////////////////////////////////////////////////////////////////////////////////////////////////////// //Simple test cases for the Stereographic Projection //Spherical tests cases used the ArcInfo 'type 1' stereographic //Ellipse test used ArcInfo 'type 2' stereographic // //Rueben Schulz //June 22, 2003 set _sphere_ = GEOGCS["WGS84", DATUM["WGS84", SPHEROID["WGS84", 6370997.0, 0],TOWGS84[0,0,0,0,0,0,0]], PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]] set _WGS84_ = GEOGCS["WGS84", DATUM["WGS84", SPHEROID["WGS84", 6378137.0, 298.257223563]], PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]] set _Stereo1_ = PROJCS["North Polar Stereographic", GEOGCS["Sphere", DATUM["Sphere", SPHEROID["Sphere", 6370997.0, 0],TOWGS84[0,0,0,0,0,0,0]],PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]], PROJECTION["Polar_Stereographic"], PARAMETER["semi_major", 6370997], PARAMETER["semi_minor", 6370997], PARAMETER["central_meridian", 50.0], PARAMETER["latitude_of_origin", 90.0], PARAMETER["scale_factor", 1.0], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["metre",1.0], AXIS["x",EAST], AXIS["y",NORTH]] set _Stereo2_ = PROJCS["South Polar Stereographic", GEOGCS["Sphere", DATUM["Sphere", SPHEROID["Sphere", 6370997.0, 0],TOWGS84[0,0,0,0,0,0,0]],PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]], PROJECTION["Polar_Stereographic"], PARAMETER["semi_major", 6370997], PARAMETER["semi_minor", 6370997], PARAMETER["central_meridian", -100.0], PARAMETER["latitude_of_origin", -90.0], PARAMETER["scale_factor", 1.0], PARAMETER["false_easting", 1000.0], PARAMETER["false_northing", 500000.0], UNIT["metre",1.0], AXIS["x",EAST], AXIS["y",NORTH]] set _Stereo3_ = PROJCS["Equitorial Stereographic", GEOGCS["Sphere", DATUM["Sphere", SPHEROID["Sphere", 6370997.0, 0],TOWGS84[0,0,0,0,0,0,0]],PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]], PROJECTION["Oblique_Stereographic"], PARAMETER["semi_major", 6370997], PARAMETER["semi_minor", 6370997], PARAMETER["central_meridian", 10.0], PARAMETER["latitude_of_origin", 0.0], PARAMETER["scale_factor", 1.0], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["metre",1.0], AXIS["x",EAST], AXIS["y",NORTH]] set _Stereo4_ = PROJCS["Equitorial Stereographic", GEOGCS["Sphere", DATUM["Sphere", SPHEROID["Sphere", 6300000.0, 0],TOWGS84[0,0,0,0,0,0,0]],PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]], PROJECTION["Oblique_Stereographic"], PARAMETER["semi_major", 6300000], PARAMETER["semi_minor", 6300000], PARAMETER["central_meridian", 180.0], PARAMETER["latitude_of_origin", 0.0], PARAMETER["scale_factor", 1.0], PARAMETER["false_easting", 1640416.66666667], PARAMETER["false_northing", 0.0], UNIT["US survey foot",0.304800609601219], AXIS["x",EAST], AXIS["y",NORTH]] set _Stereo5_ = PROJCS["Oblique Stereographic", GEOGCS["Sphere", DATUM["Sphere", SPHEROID["Sphere", 6370997.0, 0],TOWGS84[0,0,0,0,0,0,0]],PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]], PROJECTION["Oblique_Stereographic"], PARAMETER["semi_major", 6370997], PARAMETER["semi_minor", 6370997], PARAMETER["central_meridian", 140.0], PARAMETER["latitude_of_origin", 60.0], PARAMETER["scale_factor", 1.0], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["metre",1.0], AXIS["x",EAST], AXIS["y",NORTH]] set _Stereo6_ = PROJCS["Oblique Stereographic", GEOGCS["Sphere", DATUM["Sphere", SPHEROID["Sphere", 6370997.0, 0],TOWGS84[0,0,0,0,0,0,0]],PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]], PROJECTION["Oblique_Stereographic"], PARAMETER["semi_major", 6370997], PARAMETER["semi_minor", 6370997], PARAMETER["central_meridian", -120.0], PARAMETER["latitude_of_origin", -50.0], PARAMETER["scale_factor", 1.0], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["metre",1.0], AXIS["x",EAST], AXIS["y",NORTH]] set _Stereo7_ = PROJCS["North Polar Stereographic", GEOGCS["WGS84", DATUM["WGS84", SPHEROID["WGS84", 6378137.0, 298.257223563]], PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]], PROJECTION["Polar_Stereographic"], PARAMETER["semi_major", 6378137.0], PARAMETER["semi_minor", 6356752.314245179], PARAMETER["central_meridian", 150.0], PARAMETER["latitude_of_origin", 90.0], PARAMETER["scale_factor", 1.0], PARAMETER["latitude_true_scale", 85.0], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["metre",1.0], AXIS["x",EAST], AXIS["y",NORTH]] set _Stereo8_ = PROJCS["North Polar Stereographic", GEOGCS["WGS84", DATUM["WGS84", SPHEROID["WGS84", 6378137.0, 298.257223563]], PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]], PROJECTION["Polar_Stereographic"], PARAMETER["semi_major", 6378137.0], PARAMETER["semi_minor", 6356752.314245179], PARAMETER["central_meridian", -67.0], PARAMETER["latitude_of_origin", 90.0], PARAMETER["scale_factor", 1.0], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["metre",1.0], AXIS["x",EAST], AXIS["y",NORTH]] set _Stereo9_ = PROJCS["South Polar Stereographic", GEOGCS["WGS84", DATUM["WGS84", SPHEROID["WGS84", 6378137.0, 298.257223563]], PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]], PROJECTION["Polar_Stereographic"], PARAMETER["semi_major", 6378137.0], PARAMETER["semi_minor", 6356752.314245179], PARAMETER["central_meridian", -177.0], PARAMETER["latitude_of_origin", -90.0], PARAMETER["scale_factor", 1.0], PARAMETER["latitude_true_scale", 71.0], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["metre",1.0], AXIS["x",EAST], AXIS["y",NORTH]] set _Stereo10_ = PROJCS["South Polar Stereographic", GEOGCS["WGS84", DATUM["WGS84", SPHEROID["WGS84", 6378137.0, 298.257223563]], PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]], PROJECTION["Polar_Stereographic"], PARAMETER["semi_major", 6378137.0], PARAMETER["semi_minor", 6356752.314245179], PARAMETER["central_meridian", 22.0], PARAMETER["latitude_of_origin", -90.0], PARAMETER["scale_factor", 1.0], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["metre",1.0], AXIS["x",EAST], AXIS["y",NORTH]] set _Stereo11_ = PROJCS["Equitorial Stereographic", GEOGCS["WGS84", DATUM["WGS84", SPHEROID["WGS84", 6378137.0, 298.257223563]], PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]], PROJECTION["Oblique_Stereographic"], PARAMETER["semi_major", 6378137.0], PARAMETER["semi_minor", 6356752.314245179], PARAMETER["central_meridian", 56.0], PARAMETER["latitude_of_origin", 0.0], PARAMETER["scale_factor", 1.0], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["metre",1.0], AXIS["x",EAST], AXIS["y",NORTH]] set _Stereo12_ = PROJCS["Equitorial Stereographic", GEOGCS["WGS84", DATUM["WGS84", SPHEROID["WGS84", 6378137.0, 298.257223563]], PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]], PROJECTION["Oblique_Stereographic"], PARAMETER["semi_major", 6378137.0], PARAMETER["semi_minor", 6356752.314245179], PARAMETER["central_meridian", -4.0], PARAMETER["latitude_of_origin", 0.0], PARAMETER["scale_factor", 0.91], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["metre",1.0], AXIS["x",EAST], AXIS["y",NORTH]] set _Stereo13_ = PROJCS["Oblique Stereographic", GEOGCS["WGS84", DATUM["WGS84", SPHEROID["WGS84", 6378137.0, 298.257223563]], PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]], PROJECTION["Oblique_Stereographic"], PARAMETER["semi_major", 6378137.0], PARAMETER["semi_minor", 6356752.314245179], PARAMETER["central_meridian", -55.0], PARAMETER["latitude_of_origin", 51.0], PARAMETER["scale_factor", 0.99], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["metre",1.0], AXIS["x",EAST], AXIS["y",NORTH]] set _Stereo14_ = PROJCS["Oblique Stereographic", GEOGCS["WGS84", DATUM["WGS84", SPHEROID["WGS84", 6378137.0, 298.257223563]], PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]], PROJECTION["Oblique_Stereographic"], PARAMETER["semi_major", 6378137.0], PARAMETER["semi_minor", 6356752.314245179], PARAMETER["central_meridian", 64.0], PARAMETER["latitude_of_origin", -32.0], PARAMETER["scale_factor", 1.0], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["metre",1.0], AXIS["x",EAST], AXIS["y",NORTH]] //epsg oblique test set _Bessel_ = GEOGCS["Bessel 1841", DATUM["Bessel", SPHEROID["Bessel", 6377397.155, 299.15281]], PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]] set _Stereo15_ = PROJCS["Oblique Stereographic EPSG", GEOGCS["Bessel 1841", DATUM["Bessel", SPHEROID["Bessel", 6377397.155, 299.15281]], PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]], PROJECTION["Oblique_Stereographic_EPSG"], PARAMETER["semi_major", 6377397.155], PARAMETER["semi_minor", 6356078.96261866], PARAMETER["central_meridian", 5.387638889], PARAMETER["latitude_of_origin", 52.15616056], PARAMETER["scale_factor", 0.9999079], PARAMETER["false_easting", 155000.0], PARAMETER["false_northing", 463000.0], UNIT["metre",1.0], AXIS["x",EAST], AXIS["y",NORTH]] //////////////////////////////////////////// //Spherical North cs_source=_sphere_ cs_target=_Stereo1_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) pt_target=(-569039.6270, 4702290.0974) pt_source=(2.33656753, 43.96499314) pt_target=(-4001482.1176, -3645739.640) pt_source=(151.283333333, -33.916666666) pt_target=(23459842.7218, 4680642.8947) pt_source=(-47.95, -15.75) pt_target=(-16670888.6766, 2328106.9805) pt_source=(139.733333333, 35.6833333333) pt_target=(6536715.3141, -30423.4742) pt_source=(18.45, -33.9166666666) pt_target=(-12517120.9473, -20386130.6194) //////////////////////////////////////////// //Spherical South cs_source=_sphere_ cs_target=_Stereo2_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) pt_target=(-13447315.1416, 32029133.6112) pt_source=(2.33656753, 43.96499314) pt_target=(29301227.6488, -5908074.5888) pt_source=(151.283333333, -33.916666666) pt_target=(-6427015.6000, -1677847.5077) pt_source=(-47.95, -15.75) pt_target=(7606887.4054, 6431688.8876) pt_source=(139.733333333, 35.6833333333) pt_target=(-21451003.9982, -12018802.8096) pt_source=(18.45, -33.9166666666) pt_target=(5968293.0902, -2733236.3811) //////////////////////////////////////////// //Spherical Equitorial cs_source=_sphere_ cs_target=_Stereo3_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) pt_target=(-10975901.6974, 17425146.9141) pt_source=(2.33656753, 43.96499314) pt_target=( -713822.4244, 5162871.3803) pt_source=(151.283333333, -33.916666666) pt_target=(18761641.7965, -20169201.7279) pt_source=(-47.95, -15.75) pt_target=(-6880387.3644, -2289409.5189) pt_source=(139.733333333, 35.6833333333) pt_target=(16554295.1489, 15458696.4976) pt_source=(18.45, -33.9166666666) pt_target=(853342.6017, -3904711.5139) cs_source=_sphere_ cs_target=_Stereo4_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) pt_target=(18313312.9974, 23071108.5406) pt_source=(2.33656753, 43.96499314) pt_target=(-2679048.8043, 102188221.1560) pt_source=(151.283333333, -33.916666666) pt_target=(-7899377.2908, -13350238.6363) pt_source=(-47.95, -15.75) pt_target=(84776405.1579, -31575568.6281) pt_source=(139.733333333, 35.6833333333) pt_target=(-11758008.0552, 14886510.7664) pt_source=(18.45, -33.9166666666) pt_target=(-49376756.6233, -108392290.3667) //////////////////////////////////////////// //Spherical Oblique cs_source=_sphere_ cs_target=_Stereo5_ test_tolerance=(1e-2, 1e-2) //pt_source=(-123.1, 49.2166666666) //pt_target=(5111518.3838, 3519921.2597) pt_source=(2.33656753, 43.96499314) pt_target=(-4626083.2478, 7709778.2077) pt_source=(151.283333333, -33.916666666) pt_target=(2239865.7214, -13571053.7552) pt_source=(-47.95, -15.75) pt_target=(5882861.1803, 30483752.1708) pt_source=(139.733333333, 35.6833333333) pt_target=(-25202.8272, -2745174.2686) pt_source=(18.45, -33.9166666666) pt_target=(-30070247.8619, 4126498.8667) cs_source=_sphere_ cs_target=_Stereo6_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) pt_target=(-536335.5224, 14976051.8319) pt_source=(2.33656753, 43.96499314) pt_target=(35106533.5793, 8734803.7398) pt_source=(151.283333333, -33.916666666) pt_target=(-7344290.0367, -3049012.9721) pt_source=(-47.95, -15.75) pt_target=(8341703.4925, 480514.1948) pt_source=(139.733333333, 35.6833333333) pt_target=(-22133987.2972, 7312269.0941) pt_source=(18.45, -33.9166666666) pt_target=(6820762.1489, -10340013.7044) //////////////////////////////////////////// //Ellipse North cs_source=_WGS84_ cs_target=_Stereo7_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) pt_target=(4734112.9430, -256390.3922) pt_source=(2.33656753, 43.96499314) pt_target=(-2897000.2339, 4576132.3583) pt_source=(151.283333333, -33.916666666) pt_target=(531572.7533, -23728661.1588) pt_source=(-47.95, -15.75) pt_target=(5156818.7159, 15918309.4169) pt_source=(139.733333333, 35.6833333333) pt_target=(-1164785.0278, -6430662.1897) pt_source=(18.45, -33.9166666666) pt_target=(-17762444.1908, 15742538.1165) cs_source=_WGS84_ cs_target=_Stereo8_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) pt_target=(-3942631.5543, -2649338.4561) pt_source=(2.33656753, 43.96499314) pt_target=(5077292.5287, -1914846.7289) pt_source=(151.283333333, -33.916666666) pt_target=(-14732828.3540, 18666157.2113) pt_source=(-47.95, -15.75) pt_target=( 5471869.9971, -15846527.7746) pt_source=(139.733333333, 35.6833333333) pt_target=(-2945434.4200, 5847866.4702) pt_source=(18.45, -33.9166666666) pt_target=(23704915.6430, -1886433.3285) //////////////////////////////////////// //Elliptical South cs_source=_WGS84_ cs_target=_Stereo9_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) pt_target=(26745521.5872, 19503169.7905) pt_source=(2.33656753, 43.96499314) pt_target=(335506.9066, -28973967.1999) pt_source=(151.283333333, -33.916666666) pt_target=(-3476092.5649, 5624613.2248) pt_source=(-47.95, -15.75) pt_target=(7283623.7845, -5908693.9040) pt_source=(139.733333333, 35.6833333333) pt_target=(-16458682.7200, 17485900.3222) pt_source=(18.45, -33.9166666666) pt_target=(-1761438.3929, -6373133.3138) cs_source=_WGS84_ cs_target=_Stereo10_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) pt_target=(-19468934.7924, -27908063.5417) pt_source=(2.33656753, 43.96499314) pt_target=(-10023169.8500, 28050023.6648) pt_source=(151.283333333, -33.916666666) pt_target=(5261171.9431, -4303664.2698) pt_source=(-47.95, -15.75) pt_target=(-9057118.4307, 3305475.2442) pt_source=(139.733333333, 35.6833333333) pt_target=(21849835.1170, -11487640.3953) pt_source=(18.45, -33.9166666666) pt_target=(-420877.4389, 6784122.5330) //////////////////////////////////////// //Elliptical Equitorial cs_source=_WGS84_ cs_target=_Stereo11_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) pt_target=(-381517.7735, 27967467.8233) pt_source=(2.33656753, 43.96499314) pt_target=(-5196633.8476, 6180449.7189) pt_source=(151.283333333, -33.916666666) pt_target=(11438665.7948, -7672420.6307) pt_source=(-47.95, -15.75) pt_target=(-15525059.5040, -4481382.3081) pt_source=(139.733333333, 35.6833333333) pt_target=(9480395.8620, 6803267.4579) pt_source=(18.45, -33.9166666666) pt_target=(-3896229.3187, -4269843.2094) cs_source=_WGS84_ cs_target=_Stereo12_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) pt_target=(-9764791.6464, 12867824.7061) pt_source=(2.33656753, 43.96499314) pt_target=(538591.9946, 4675248.9850) pt_source=(151.283333333, -33.916666666) pt_target=(16501209.0722, -26357916.0339) pt_source=(-47.95, -15.75) pt_target=(-4581578.1738, -1849327.3326) pt_source=(139.733333333, 35.6833333333) pt_target=(16269126.0087, 19618474.4606) pt_source=(18.45, -33.9166666666) pt_target=(2084362.8747, -3645505.7204) //////////////////////////////////////// //Elliptical Oblique cs_source=_WGS84_ cs_target=_Stereo13_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) pt_target=(-4407434.5194, 2082811.7883) pt_source=(2.33656753, 43.96499314) pt_target=(4298278.4273, 955057.5666) pt_source=(151.283333333, -33.916666666) pt_target=(-47167588.1622, 29071179.3031) pt_source=(-47.95, -15.75) pt_target=(1067938.1787, -8233254.2708) pt_source=(139.733333333, 35.6833333333) pt_target=(-2739234.3871, 12928389.0428) pt_source=(18.45, -33.9166666666) pt_target=(13968985.4665, -9367398.2205) cs_source=_WGS84_ cs_target=_Stereo14_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) pt_target=(21084353.9245, 77613853.8593) pt_source=(2.33656753, 43.96499314) pt_target=(-8743580.0850, 10572818.2134) pt_source=(151.283333333, -33.916666666) pt_target=(7981081.4514, -4334411.3396) pt_source=(-47.95, -15.75) pt_target=(-13608055.3148, -6384283.3272) pt_source=(139.733333333, 35.6833333333) pt_target=(11635190.0266, 8836389.2762) pt_source=(18.45, -33.9166666666) pt_target=(-4231988.1148, -1174790.7406) cs_source=_Bessel_ cs_target=_Stereo15_ test_tolerance=(1e-2, 1e-2) pt_source=(6.0,53.0) pt_target=(196105.283, 557057.739)