////////////////////////////////////////////////////////////////////////////////////////////////////// //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["Sphere", DATUM["Sphere", 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 _sphere2_ = GEOGCS["Sphere", DATUM["Sphere", SPHEROID["WGS84", 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]] 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_Series"], 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["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["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["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["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["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["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["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["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"], 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=_sphere2_ 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) //a test point from the EPSG guide book 7 cs_source=_Bessel_ cs_target=_Stereo15_ test_tolerance=(1e-2, 1e-2) pt_source=(6.0,53.0) pt_target=(196105.283, 557057.739) ////////////////////////////////////////////////////////////////////////////////////////////////////// //More tests for the Stereographic Projection //These test the EPSG form of the equations and were //calculated using the 'Double Stereographic' projection in ArcGIS 8.3 //Note: there are some differences beteen what arcGIS calculates and //what we calculate. // //Rueben Schulz //Feb 17, 2004 set _DStereo1_ = 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", 0.0], PARAMETER["latitude_of_origin", 70.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 _DStereo2_ = 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", -110.0], PARAMETER["latitude_of_origin", -50.0], PARAMETER["scale_factor", 0.96], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["metre",1.0], AXIS["x",EAST], AXIS["y",NORTH]] set _DStereo3_ = 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", 20.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 _DStereo4_ = PROJCS["Polar Stereographic (North)", 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", 140.0], PARAMETER["latitude_of_origin", 90.0], PARAMETER["scale_factor", 0.8], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["metre",1.0], AXIS["x",EAST], AXIS["y",NORTH]] set _DStereo5_ = PROJCS["Polar Stereographic (South)", 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", 0.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]] //////////////////////////////////////// //Elliptical Oblique cs_source=_WGS84_ cs_target=_DStereo1_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) pt_target=(-4398811.24918829, 4780023.58535872) pt_source=(2.33656753, 43.96499314) pt_target=(197487.535049174, -2946384.25777934) //pt_source=(151.283333333, -33.9166666666) //Error: The operation was attempted on an empty geometry.::"esri.Point" pt_source=(-47.95, -15.75) pt_target=(-9417083.95801759, -9170019.17553014) pt_source=(139.733333333, 35.6833333333) pt_target=(5013546.19775259, 7475371.30940516) pt_source=(18.45, -33.9166666666) pt_target=(4486418.53520613, -15840660.7109299) pt_source=(0, 90) pt_target=(0, 2255899.32774423) //pt_source=(0, -90) //Error: The operation was attempted on an empty geometry.::"esri.Point" cs_source=_WGS84_ cs_target=_DStereo2_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) pt_target=(-2187792.49332045, 14311891.5579733) //pt_source=(2.33656753, 43.96499314) //Error: The operation was attempted on an empty geometry.::"esri.Point" //error pt_source=(151.283333333, -33.9166666666) pt_target=(-7460003.39378955, -4124311.05369635) pt_source=(-47.95, -15.75) pt_target=(6949231.62438508, 1382215.87566019) //error pt_source=(139.733333333, 35.6833333333) pt_target=(-24924065.4478382, 5124122.41314095) pt_source=(18.45, -33.9166666666) pt_target=(7265980.56818263, -8446857.93987766) //pt_source=(0, 90) //Error: The operation was attempted on an empty geometry.::"esri.Point" pt_source=(0, -90) pt_target=(3.9914448559467E-10, -4464514.77931394) //////////////////////////////////////// //Elliptical Equitorial cs_source=_WGS84_ cs_target=_DStereo3_ test_tolerance=(1e-2, 1e-2) //error pt_source=(-123.1, 49.2166666666) pt_target=(-10422471.5923968, 20303079.0429489) //error pt_source=(2.33656753, 43.96499314) pt_target=(-1767825.34583943, 5241329.09731566) //error pt_source=(151.283333333, -33.9166666666) pt_target=(17436880.6584282, -15376813.6457952) //error pt_source=(-47.95, -15.75) pt_target=(-8556259.14610263, -2569702.00752923) //error pt_source=(139.733333333, 35.6833333333) pt_target=(14910741.4194107, 12198910.9989865) //error pt_source=(18.45, -33.9166666666) pt_target=(-278548.416542974, -3869437.88697763) pt_source=(0, 90) pt_target=(-2.82499596765061E-10, 12713504.6284904) pt_source=(0, -90) pt_target=(-2.82499596765061E-10, -12713504.6284904) //////////////////////////////////////// //Elliptical Polar cs_source=_WGS84_ cs_target=_DStereo4_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) pt_target=(3772547.94163496, 456528.463621665) pt_source=(2.33656753, 43.96499314) pt_target=(-2923662.83596991, 3208946.80055835) pt_source=(151.283333333, -33.9166666666) pt_target=(3722230.51407769, -18656185.5710066) pt_source=(-47.95, -15.75) pt_target=(1854960.95732183, 13282829.3019208) pt_source=(139.733333333, 35.6833333333) pt_target=(-24379.6291743415, -5238149.13743781) pt_source=(18.45, -33.9166666666) pt_target=(-16211849.3469659, 9954104.7217239) pt_source=(0, 90) pt_target=(0, 0) //pt_source=(0, -90) //Error: The operation was attempted on an empty geometry.::"esri.Point" //######################################################################## //the following tests all failed. The commented out points are from //arcGIS double stereographic. The new values are from the arcInfo stereographic. //Note that they should be the same for the polar case (I think???) //Points calculated using double stereographic and stereographic are equal //for the north polar case above (as expected). cs_source=_WGS84_ cs_target=_DStereo5_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) //Error: The operation was attempted on an empty geometry.::"esri.Point" pt_target=(-28505826.6173, -18582714.5881) pt_source=(2.33656753, 43.96499314) //pt_target=(1230808.15179152, 30164338.8749791) pt_target=(1214402.5538, 29762274.5749) pt_source=(151.283333333, -33.9166666666) //pt_target=(3310012.24359201, -6041688.15332015) pt_target=(3265892.6704, -5961157.7254) pt_source=(-47.95, -15.75) //pt_target=(-7256078.34586213, 6544877.05627666) pt_target=(-7159361.1570, 6457639.5596) pt_source=(139.733333333, 35.6833333333) //pt_target=(16171014.0080265, -19090709.4904612) pt_target=(15955468.5106, -18836247.0015) pt_source=(18.45, -33.9166666666) //pt_target=(2180206.51455429, 6534896.81239062) pt_target=(2151146.2653, 6447792.3436) //pt_source=(0, 90) //Error: The operation was attempted on an empty geometry.::"esri.Point" pt_source=(0, -90) pt_target=(0, 0) set _sphere3_ = GEOGCS["Sphere", DATUM["Sphere", SPHEROID["Sphere", 6371000.0, 0]], PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]] set _DStereo6_ = PROJCS["Polar Stereographic (South)", GEOGCS["Sphere", DATUM["Sphere", SPHEROID["Sphere", 6371000.0, 0]], PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]], PROJECTION["Polar_Stereographic"], PARAMETER["semi_major", 6371000.0], PARAMETER["semi_minor", 6371000.0], PARAMETER["central_meridian", 0.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 _DStereo7_ = PROJCS["Equitorial Stereographic", GEOGCS["Sphere", DATUM["Sphere", SPHEROID["Sphere", 6371000.0, 0]], PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]], PROJECTION["Polar_Stereographic"], PARAMETER["central_meridian", -45.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 _DStereo8_ = PROJCS["Oblique Stereographic", GEOGCS["Sphere", DATUM["Sphere", SPHEROID["Sphere", 6371000.0, 0]], PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]], PROJECTION["Oblique_Stereographic"], PARAMETER["central_meridian", -55.0], PARAMETER["latitude_of_origin", 70.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 _DStereo9_ = PROJCS["Oblique Stereographic", GEOGCS["Sphere", DATUM["Sphere", SPHEROID["Sphere", 6371000.0, 0]], PRIMEM["Greenwich", 0.0], UNIT["degree",0.017453292519943295], AXIS["Longitude",EAST], AXIS["Latitude",NORTH]], PROJECTION["Oblique_Stereographic"], PARAMETER["central_meridian", 170.0], PARAMETER["latitude_of_origin", -22.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]] //////////////////////////////////////// //Spherical Polar cs_source=_sphere3_ cs_target=_DStereo6_ test_tolerance=(1e-2, 1e-2) //pt_source=(-123.1, 49.2166666666) //Error: The operation was attempted on an empty geometry.::"esri.Point" pt_source=(2.33656753, 43.96499314) pt_target=(1222790.97629403, 29967855.9396212) pt_source=(151.283333333, -33.9166666666) pt_target=(3260975.84097009, -5952183.1512236) pt_source=(-47.95, -15.75) pt_target=(-7162325.06507552, 6460312.9615793) pt_source=(139.733333333, 35.6833333333) pt_target=(16053723.026356, -18952241.4849415) pt_source=(18.45, -33.9166666666) pt_target=(2147907.69612739, 6438085.13690347) //pt_source=(0, 90) //Error: The operation was attempted on an empty geometry.::"esri.Point" pt_source=(0, -90) pt_target=(0, 0) //////////////////////////////////////// //Spherical Equitorial cs_source=_sphere3_ cs_target=_DStereo7_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) pt_target=(-7177457.04562655, 8502789.26341092) pt_source=(2.33656753, 43.96499314) pt_target=(4532964.6771832, 5945604.3806171) //pt_source=(151.283333333, -33.9166666666) //Error: The operation was attempted on an empty geometry.::"esri.Point" pt_source=(-47.95, -15.75) pt_target=(-321816.428830414, -1763578.32357749) //pt_source=(139.733333333, 35.6833333333) //Error: The operation was attempted on an empty geometry.::"esri.Point" pt_source=(18.45, -33.9166666666) pt_target=(6899624.53423796, -5186179.56202058) //the long. of the the next 2 points were changed from arcGIS output (was 0) //we have differing views about what the long. of a pole should be. pt_source=(-45, 90) pt_target=(5.5168238290053E-10, 12742000) pt_source=(-45, -90) pt_target=(5.5168238290053E-10, -12742000) cs_source=_DStereo7_ cs_target=_sphere3_ test_tolerance=(1e-5, 1e-5) pt_source=(5.5168238290053E-10, 12742000) pt_target=(-45, 90) pt_source=(5.5168238290053E-10, -12742000) pt_target=(-45, -90) //////////////////////////////////////// //Spherical Oblique cs_source=_sphere3_ cs_target=_DStereo8_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) pt_target=(-4302563.66689065, 213187.830274701) pt_source=(2.33656753, 43.96499314) pt_target=(4324891.73045712, -910713.807515885) //pt_source=(151.283333333, -33.9166666666) //Error: The operation was attempted on an empty geometry.::"esri.Point" pt_source=(-47.95, -15.75) pt_target=(1404583.27391585, -11776410.5577183) pt_source=(139.733333333, 35.6833333333) pt_target=(-2057243.67777944, 9338231.09361794) pt_source=(18.45, -33.9166666666) pt_target=(18213225.8877641, -9455459.3973766) pt_source=(0, 90) pt_target=(3.29485054216154E-10, 2246758.38818726) //pt_source=(0, -90) //Error: The operation was attempted on an empty geometry.::"esri.Point" cs_source=_sphere3_ cs_target=_DStereo9_ test_tolerance=(1e-2, 1e-2) pt_source=(-123.1, 49.2166666666) pt_target=(8025167.41865162, 10659462.713673) //pt_source=(2.33656753, 43.96499314) //Error: The operation was attempted on an empty geometry.::"esri.Point" pt_source=(151.283333333, -33.9166666666) pt_target=(-1751021.5482452, -1465895.97287927) pt_source=(-47.95, -15.75) pt_target=(18949003.1287529, -17159231.214752) pt_source=(139.733333333, 35.6833333333) pt_target=(-3642974.56059216, 7151031.37033279) pt_source=(18.45, -33.9166666666) pt_target=(-9459391.5721868, -18919110.2925948) pt_source=(0, 90) pt_target=(-2.16631219873976E-10, 18890791.8607893) pt_source=(0, -90) pt_target=(-9.85589165790218E-11, -8594587.5216062)