From 4578e9c7ae264ad706876cc2bb6a412ff1681fe0 Mon Sep 17 00:00:00 2001 From: Kristian Evers Date: Sun, 1 Apr 2018 20:48:22 +0200 Subject: Add tests extra test cases for the gnomonic projection --- test/gie/builtins.gie | 110 +++++++++++++++++++++++++++++++++++++++++--------- 1 file changed, 92 insertions(+), 18 deletions(-) (limited to 'test') diff --git a/test/gie/builtins.gie b/test/gie/builtins.gie index 9cfb7cce..27d5c143 100644 --- a/test/gie/builtins.gie +++ b/test/gie/builtins.gie @@ -1366,28 +1366,102 @@ Gnomonic =============================================================================== ------------------------------------------------------------------------------- -operation +proj=gnom +a=6400000 +lat_1=0.5 +lat_2=2 +Test material from Snyder p. 168, table 26. +Tests the equatorial aspect of the projection. +------------------------------------------------------------------------------- +operation +proj=gnom +R=1 ------------------------------------------------------------------------------- tolerance 0.1 mm -accept 2 1 -expect 223492.924747185 111780.509206593 -accept 2 -1 -expect 223492.924747185 -111780.509206593 -accept -2 1 -expect -223492.924747185 111780.509206593 -accept -2 -1 -expect -223492.924747185 -111780.509206593 +accept 0 0 +expect 0 0 +roundtrip 100 +accept 10 80 +expect 0.1763 5.7588 +roundtrip 100 +accept 20 70 +expect 0.3640 2.9238 +roundtrip 100 +accept 30 60 +expect 0.5774 2.0000 +roundtrip 100 +accept 40 50 +expect 0.8391 1.5557 +roundtrip 100 +accept 50 40 +expect 1.1918 1.3054 +roundtrip 100 +accept 60 30 +expect 1.7321 1.1547 +roundtrip 100 +accept 70 20 +expect 2.7475 1.0642 +roundtrip 100 +accept 80 10 +expect 5.6713 1.0154 +roundtrip 100 +accept 80 80 +expect 5.6713 32.6596 +roundtrip 100 +accept 0 90 +expect failure errno tolerance_condition -direction inverse -accept 200 100 -expect 0.001790493 0.000895247 -accept 200 -100 -expect 0.001790493 -0.000895247 -accept -200 100 -expect -0.001790493 0.000895247 -accept -200 -100 -expect -0.001790493 -0.000895247 +# test that extreme northings are mapped to the sphere +direction inverse +accept 0 1e8 +expect 0 90 + + +------------------------------------------------------------------------------- +Test the northern polar aspect of the gnonomic projection +------------------------------------------------------------------------------- +operation +proj=gnom +R=1 +lat_0=90 +------------------------------------------------------------------------------- +tolerance 0.1 mm +accept 0 90 +expect 0 0 +roundtrip 100 +accept 45 45 +expect 0.7071 -0.7071 +roundtrip 100 +accept 0 0 +expect failure errno tolerance_condition +accept 90 0 +expect failure errno tolerance_condition +------------------------------------------------------------------------------- +Test the sourthern polar aspect of the gnonomic projection +------------------------------------------------------------------------------- +operation +proj=gnom +R=1 +lat_0=-90 +------------------------------------------------------------------------------- +tolerance 0.1 mm +accept 0 -90 +expect 0 0 +roundtrip 100 +accept 45 -45 +expect 0.7071 0.7071 +roundtrip 100 +accept 0 0 +expect failure errno tolerance_condition +accept 90 0 +expect failure errno tolerance_condition + +------------------------------------------------------------------------------- +Test the oblique aspect of the gnonomic projection +------------------------------------------------------------------------------- +operation +proj=gnom +R=1 +lat_0=45 +------------------------------------------------------------------------------- +tolerance 0.1 mm +accept 0 45 +expect 0 0 +roundtrip 100 +accept 0 0 +expect 0 -1 +roundtrip 100 +accept 0 90 +expect 0 1 +roundtrip 100 +accept 0 -45 +expect failure errno tolerance_condition =============================================================================== Goode Homolosine -- cgit v1.2.3 From 1777869afa75e3e9ee3cdc4fd0b777d05b4f51d6 Mon Sep 17 00:00:00 2001 From: Kristian Evers Date: Mon, 2 Apr 2018 00:28:05 +0200 Subject: Add test cases for aeqd projection --- test/gie/builtins.gie | 252 +++++++++++++++++++++++++++++++++++++++++++------- 1 file changed, 217 insertions(+), 35 deletions(-) (limited to 'test') diff --git a/test/gie/builtins.gie b/test/gie/builtins.gie index 27d5c143..4d744135 100644 --- a/test/gie/builtins.gie +++ b/test/gie/builtins.gie @@ -71,51 +71,233 @@ Azimuthal Equidistant =============================================================================== ------------------------------------------------------------------------------- -operation +proj=aeqd +ellps=GRS80 +lat_1=0.5 +lat_2=2 +Test equatorial aspect of the spherical azimuthal equidistant. Test data from +Snyder pp. 196-197, table 30. +------------------------------------------------------------------------------- +operation +proj=aeqd +R=1 +lat_0=0 ------------------------------------------------------------------------------- tolerance 0.1 mm -accept 2 1 -expect 222616.522190052 110596.996549550 -accept 2 -1 -expect 222616.522190052 -110596.996549550 -accept -2 1 -expect -222616.522190052 110596.996549550 -accept -2 -1 -expect -222616.522190052 -110596.996549550 +accept 0 0 +expect 0 0 +roundtrip 100 +accept 0 90 +expect 0 1.57080 +roundtrip 100 +accept 10 80 +expect 0.04281 1.39829 +roundtrip 100 +accept 40 30 +expect 0.62896 0.56493 +roundtrip 100 +accept 90 0 +expect 1.57080 0 +roundtrip 100 +accept 90 90 +expect 0 1.57080 +roundtrip 100 -direction inverse -accept 200 100 -expect 0.001796631 0.000904369 -accept 200 -100 -expect 0.001796631 -0.000904369 -accept -200 100 -expect -0.001796631 0.000904369 -accept -200 -100 -expect -0.001796631 -0.000904369 +# point opposite projection center is undefined +accept 180 0 +expect failure errno tolerance_condition ------------------------------------------------------------------------------- -operation +proj=aeqd +R=6400000 +lat_1=0.5 +lat_2=2 +Test equatorial aspect of the ellipsoidal azimuthal equidistant. Test data from +Snyder pp. 196-197, table 30. +------------------------------------------------------------------------------- +operation +proj=aeqd +ellps=GRS80 +lat_0=0 ------------------------------------------------------------------------------- tolerance 0.1 mm -accept 2 1 -expect 223379.456047271 111723.757570854 -accept 2 -1 -expect 223379.456047271 -111723.757570854 -accept -2 1 -expect -223379.456047271 111723.757570854 -accept -2 -1 -expect -223379.456047271 -111723.757570854 +accept 0 90 +expect 0 10001965.7292 +roundtrip 100 +accept 0 0 +expect 0 0 +roundtrip 100 +accept 90 0 +expect 10_018_754.1714 0 +roundtrip 100 +accept 90 0 +expect 10_018_754.1714 0 +roundtrip 100 +accept 45 45 +expect 3_860_398.3783 5_430_089.0490 +roundtrip 100 + +------------------------------------------------------------------------------- +Test equatorial aspect of the ellipsoidal azimuthal equidistant. Test data from +Snyder pp. 196-197, table 30. +------------------------------------------------------------------------------- +operation +proj=aeqd +ellps=GRS80 +guam +------------------------------------------------------------------------------- +tolerance 1 m +accept 0 0 +expect 0.0000 0.0000 +roundtrip 100 +accept 90 90 +except 0.0000 10_001_965.7292 +roundtrip 100 +accept 0 90 +expect 0.0000 10_001_965.7292 +roundtrip 100 +accept 90 90 +expect 0.0000 10_001_965.7292 +roundtrip 100 +accept 45 45 +expect 3548107.5793 5970183.542 +#roundtrip 100 +accept -45 -45 +expect -3548107.5793 -5970183.542 +#roundtrip 100 + +------------------------------------------------------------------------------- +Test northern polar aspect of the ellipsoidal azimuthal equidistant. Test data +from Snyder p. 198, table 31. +------------------------------------------------------------------------------- +operation +proj=aeqd +ellps=intl +lat_0=90 +------------------------------------------------------------------------------- +tolerance 0.1 m +accept 0 90 +expect 0 0 +roundtrip 100 +accept 0 85 +expect 0 -558_485.4 +roundtrip 100 +accept 0 80 +expect 0 -1_116_885.2 +roundtrip 100 +accept 0 70 +expect 0 -2_233_100.9 +roundtrip 100 + +------------------------------------------------------------------------------- +Test sourthern polar aspect of the ellipsoidal azimuthal equidistant. Test data +from Snyder p. 198, table 31. +------------------------------------------------------------------------------- +operation +proj=aeqd +ellps=intl +lat_0=-90 +------------------------------------------------------------------------------- +tolerance 0.1 m +accept 0 -90 +expect 0 0 +roundtrip 100 +accept 0 -85 +expect 0 558_485.4 +roundtrip 100 +accept 0 -80 +expect 0 1_116_885.2 +roundtrip 100 +accept 0 -70 +expect 0 2_233_100.9 +roundtrip 100 + +------------------------------------------------------------------------------- +Test northern polar aspect of the spherical azimuthal equidistant. +------------------------------------------------------------------------------- +operation +proj=aeqd +R=1 +lat_0=90 +------------------------------------------------------------------------------- +tolerance 0.1 m +accept 0 0 +expect 0 -1.5708 +roundtrip 100 +accept 0 90 +expect 0 0 +roundtrip 100 +accept 90 90 +expect 0 0 +roundtrip 100 +accept 90 0 +expect 1.5708 0 +roundtrip 100 +accept 45 45 +expext 0.5554 -0.5554 +roundtrip 100 + +#point opposite of projection center is undefined +accept 0 -90 +expect failure errno tolerance_condition direction inverse -accept 200 100 -expect 0.001790493 0.000895247 -accept 200 -100 -expect 0.001790493 -0.000895247 -accept -200 100 -expect -0.001790493 0.000895247 -accept -200 -100 -expect -0.001790493 -0.000895247 +accept 0 5 +expect failure errno tolerance_condition + +accept 0 3.14159265359 +expect 180 -90 +------------------------------------------------------------------------------- +Test sourthnern polar aspect of the spherical azimuthal equidistant. +------------------------------------------------------------------------------- +operation +proj=aeqd +R=1 +lat_0=-90 +------------------------------------------------------------------------------- +tolerance 0.1 m +accept 0 0 +expect 0 1.5708 +roundtrip 100 +accept 0 -90 +expect 0 0 +roundtrip 100 +accept 90 -90 +expect 0 0 +roundtrip 100 +accept 90 0 +expect 1.5708 0 +roundtrip 100 +accept 45 -45 +expext 0.5554 -0.5554 +roundtrip 100 + +#point opposite of projection center is undefined +accept 0 90 +expect failure errno tolerance_condition + + +------------------------------------------------------------------------------- +Test oblique aspect of the spherical azimuthal equidistant. +------------------------------------------------------------------------------- +operation +proj=aeqd +R=1 +lat_0=45 +------------------------------------------------------------------------------- +tolerance 0.1 m +accept 0 0 +expect 0.0000 -0.7854 +roundtrip 100 +accept 0 45 +expect 0.0000 0.0000 +roundtrip 100 +accept 0 90 +expect 0.0000 0.7854 +roundtrip 100 +accept 90 0 +expect 1.5708 -0.0000 +roundtrip 100 +accept 90 45 +expect 0.8550 0.6046 +#roundtrip 100 # roundtrip performs badly for this test on some platforms +accept 90 90 +expect 0.0000 0.7854 +roundtrip 100 + +------------------------------------------------------------------------------- +Test oblique aspect of the ellipsoidal azimuthal equidistant. +------------------------------------------------------------------------------- +operation +proj=aeqd +ellps=GRS80 +lat_0=45 +------------------------------------------------------------------------------- +tolerance 0.1 m +accept 0 0 +expect 0.0000 -4984944.3779 +roundtrip 100 +accept 0 45 +expect 0.0000 0.0000 +roundtrip 100 +accept 0 90 +expect 0.0000 5017021.3514 +roundtrip 100 +accept 90 0 +expect 10010351.5666 26393.3781 +roundtrip 100 +accept 90 45 +expect 5461910.9128 3863514.7047 +roundtrip 100 +accept 90 90 +expect 0.0000 5017021.3514 +roundtrip 100 =============================================================================== Airy -- cgit v1.2.3 From 3a2ddb6c6efbccf388ea89e177ca51fd25946ecf Mon Sep 17 00:00:00 2001 From: Even Rouault Date: Sat, 7 Apr 2018 13:41:12 +0200 Subject: Add webmerc projection --- test/gie/4D-API_cs2cs-style.gie | 13 +++++++++++++ 1 file changed, 13 insertions(+) (limited to 'test') diff --git a/test/gie/4D-API_cs2cs-style.gie b/test/gie/4D-API_cs2cs-style.gie index bbd9ee39..6eed4faa 100644 --- a/test/gie/4D-API_cs2cs-style.gie +++ b/test/gie/4D-API_cs2cs-style.gie @@ -215,6 +215,19 @@ accept 487147.594520173 4934316.46263998 0 expect -10370728.80 5552839.74 0 ------------------------------------------------------------------------------- +------------------------------------------------------------------------------- +Test that Google's Web Mercator with +proj=webmerc +------------------------------------------------------------------------------- +operation proj=pipeline step init=epsg:26915 inv step proj=webmerc datum=WGS84 +------------------------------------------------------------------------------- +tolerance 20 cm +accept 487147.594520173 4934316.46263998 +expect -10370728.80 5552839.74 + +accept 487147.594520173 4934316.46263998 0 +expect -10370728.80 5552839.74 0 +------------------------------------------------------------------------------- + ------------------------------------------------------------------------------- Test that +datum parameters are handled correctly in pipelines. See #872 for details. -- cgit v1.2.3 From b7c59722cb0058391e47d2899d510961a0985c59 Mon Sep 17 00:00:00 2001 From: Kristian Evers Date: Mon, 9 Apr 2018 22:28:20 +0200 Subject: Add more test cases for airy projection --- test/gie/builtins.gie | 68 +++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 68 insertions(+) (limited to 'test') diff --git a/test/gie/builtins.gie b/test/gie/builtins.gie index 4d744135..84d5c832 100644 --- a/test/gie/builtins.gie +++ b/test/gie/builtins.gie @@ -318,6 +318,74 @@ expect -189109.886908621 94583.752387504 accept -2 -1 expect -189109.886908621 -94583.752387504 +------------------------------------------------------------------------------- +Test north polar aspect +------------------------------------------------------------------------------- +operation +proj=airy +R=1 +lat_0=90 +------------------------------------------------------------------------------- +tolerance 0.1 mm +accept 0 0 +expect 0 -1.3863 +accept 0 90 +expect 0 0 +accept 0 -90 +expect failure errno tolerance_condition + + +------------------------------------------------------------------------------- +Test south polar aspect +------------------------------------------------------------------------------- +operation +proj=airy +R=1 +lat_0=-90 +------------------------------------------------------------------------------- +tolerance 0.1 mm +accept 0 0 +expect 0 1.3863 +accept 0 -90 +expect 0 0 +accept 0 90 +expect failure errno tolerance_condition + +------------------------------------------------------------------------------- +Test oblique aspect +------------------------------------------------------------------------------- +operation +proj=airy +R=1 +lon_0=45 +lat_0=45 +------------------------------------------------------------------------------- +tolerance 0.1 mm +accept 45 45 +expect 0 0 +accept 0 0 +expext -0.7336 -0.5187 +accept -45 -45 +expect failure errno tolerance_condition + +------------------------------------------------------------------------------- +Test that coordinates on the opposing hemisphere are projected when using ++no_cut. +------------------------------------------------------------------------------- +operation +proj=airy +R=1 +lat_0=-90 +no_cut +------------------------------------------------------------------------------- +tolerance 0.1 mm +accept 0 10 +expect 0 1.5677 + + +------------------------------------------------------------------------------- +Test the +lat_b parameter +------------------------------------------------------------------------------- +operation +proj=airy +R=1 +lat_b=89.99999999 # check tolerance +------------------------------------------------------------------------------- +tolerance 0.1 mm +accept 0 0 +expect 0 0 +------------------------------------------------------------------------------- +operation +proj=airy +R=1 +lat_b=30 +------------------------------------------------------------------------------- +tolerance 0.1 mm +accept 0 0 +expect 0 0 +accept 25 25 +expect 0.3821 0.4216 + =============================================================================== Aitoff -- cgit v1.2.3 From 8bd2fc699e61eb3665b7d16d450886a914258933 Mon Sep 17 00:00:00 2001 From: Kristian Evers Date: Wed, 11 Apr 2018 00:03:03 +0200 Subject: Add more tests for laea --- test/gie/builtins.gie | 151 +++++++++++++++++++++++++++++++++++++++++++++++++- 1 file changed, 149 insertions(+), 2 deletions(-) (limited to 'test') diff --git a/test/gie/builtins.gie b/test/gie/builtins.gie index 4d744135..47d63bf1 100644 --- a/test/gie/builtins.gie +++ b/test/gie/builtins.gie @@ -2143,7 +2143,7 @@ Lambert Azimuthal Equal Area =============================================================================== ------------------------------------------------------------------------------- -operation +proj=laea +ellps=GRS80 +lat_1=0.5 +lat_2=2 +operation +proj=laea +ellps=GRS80 ------------------------------------------------------------------------------- tolerance 0.1 mm accept 2 1 @@ -2155,6 +2155,9 @@ expect -222602.471450095 110589.827224410 accept -2 -1 expect -222602.471450095 -110589.827224409 +accept 180 0 +expect failure errno tolerance_condition + direction inverse accept 200 100 expect 0.001796631 0.000904369 @@ -2166,7 +2169,7 @@ accept -200 -100 expect -0.001796631 -0.000904369 ------------------------------------------------------------------------------- -operation +proj=laea +R=6400000 +lat_1=0.5 +lat_2=2 +operation +proj=laea +R=6400000 ------------------------------------------------------------------------------- tolerance 0.1 mm accept 2 1 @@ -2178,6 +2181,9 @@ expect -223365.281370125 111716.668072916 accept -2 -1 expect -223365.281370125 -111716.668072916 +accept 180 0 +expect failure errno tolerance_condition + direction inverse accept 200 100 expect 0.001790493 0.000895247 @@ -2188,6 +2194,147 @@ expect -0.001790493 0.000895247 accept -200 -100 expect -0.001790493 -0.000895247 +------------------------------------------------------------------------------- +Test oblique aspect of the spherical form +------------------------------------------------------------------------------- +operation +proj=laea +R=1 +lat_0=45 +------------------------------------------------------------------------------- +tolerance 0.1 mm +accept 0 45 +expect 0 0 +accept 0 0 +expect 0 -0.7654 +accept 0 90 +expect 0 0.7654 +accept 0 -45 +expect 0 -1.4142 +accept 45 45 +expect 0.5194 0.1521 + +tolerance 0.1 mm +accept 45 45 +roundtrip 100 + +# error when waaay outside the sphere +direction inverse +accept 0 10 +expect failure errno tolerance_condition + +------------------------------------------------------------------------------- +Test oblique aspect of the ellipsoidal form +------------------------------------------------------------------------------- +operation +proj=laea +ellps=GRS80 +lat_0=45 +------------------------------------------------------------------------------- +tolerance 0.1 mm +accept 0 45 +expect 0 0 +accept 0 0 +expect 0 -4860248.8602 +accept 0 90 +expect 0 4886594.2207 +accept 0 -45 +expect 0 -8984728.0442 +accept 45 45 +expect 3318800.8682 968788.2336 + +tolerance 10 cm +accept 45 45 +roundtrip 100 + +# test rho Date: Wed, 11 Apr 2018 23:11:41 +0200 Subject: Add more tests for nsper --- test/gie/builtins.gie | 54 +++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 54 insertions(+) (limited to 'test') diff --git a/test/gie/builtins.gie b/test/gie/builtins.gie index 4f6d6a44..fb8a53c3 100644 --- a/test/gie/builtins.gie +++ b/test/gie/builtins.gie @@ -3347,6 +3347,60 @@ expect -0.001790493 0.000895247 accept -200 -100 expect -0.001790493 -0.000895247 +------------------------------------------------------------------------------- +Test north polar aspect +------------------------------------------------------------------------------- +operation +proj=nsper +R=1 +h=3 +lat_0=90 +------------------------------------------------------------------------------- +tolerance 0.1 mm +accept 0 90 +expect 0 0 +accept 45 45 +expect 0.4555 -0.4555 +roundtrip 100 +accept 0 0 +expect failure errno tolerance_condition + +direction inverse +accept 0 0 +expect 0 90 +accept 0 2 # projected coordinate is outside the sphere +expect failure errno tolerance_condition + +------------------------------------------------------------------------------- +Test south polar aspect +------------------------------------------------------------------------------- +operation +proj=nsper +R=1 +h=3 +lat_0=-90 +------------------------------------------------------------------------------- +tolerance 0.1 mm +accept 0 -90 +expect 0 0 +accept -45 -45 +expect -0.4555 0.4555 +roundtrip 100 +accept 0 0 +expect failure errno tolerance_condition + +------------------------------------------------------------------------------- +operation +proj=nsper +R=1 +h=3 +lat_0=45 +------------------------------------------------------------------------------- +tolerance 0.1 mm +accept 0 45 +expect 0 0 +accept 0 90 +expect 0 0.6442 +accept 0 0 +expect 0 -0.6442 +accept 45 45 +expect 0.4767 0.1396 +roundtrip 100 + + +------------------------------------------------------------------------------- +operation +proj=nsper +R=1 +h=0 +------------------------------------------------------------------------------- +expect failure errno h_less_than_zero + =============================================================================== New Zealand Map Grid -- cgit v1.2.3 From 299135153e33869157ecd432dd0194743e959170 Mon Sep 17 00:00:00 2001 From: Kristian Evers Date: Fri, 13 Apr 2018 22:42:54 +0200 Subject: Fix segfault in deformation When transforming coordinates outside the grid model the deformation operation failed spectatularly. This is now fixed by checking that the coordinate is inside the grid. If it isn't an error is returned and a debugging log message is issued. Closes #934 --- test/gie/deformation.gie | 8 ++++++++ 1 file changed, 8 insertions(+) (limited to 'test') diff --git a/test/gie/deformation.gie b/test/gie/deformation.gie index 55283b84..b6ca3e0f 100644 --- a/test/gie/deformation.gie +++ b/test/gie/deformation.gie @@ -23,6 +23,14 @@ accept -3004295.5882503074 -1093474.1690603832 5500477.1338251457 expect -3004295.7025 -1093474.2106 5500477.3444 roundtrip 5 +# Test that errors are reported for coordinates outside the grid. +# Here we test 120W 40N which is well outside the alaska grid. +accept -2446353.8001 -4237209.0750 4077985.572 +expect failure errno grid_area +accept -2446353.8001 -4237209.0750 4077985.572 2000 +expect failure errno grid_area + + ------------------------------------------------------------------------------- Test using both horizontal and vertical grids ------------------------------------------------------------------------------- -- cgit v1.2.3 From 111ce96fcadd9bd700c7000d6f22fea60211d8c7 Mon Sep 17 00:00:00 2001 From: Kristian Evers Date: Sun, 15 Apr 2018 18:35:38 +0200 Subject: Update +proj=aeqd +guam test case with data from EPSG Note 7-2 --- test/gie/builtins.gie | 34 ++++++++++++---------------------- 1 file changed, 12 insertions(+), 22 deletions(-) (limited to 'test') diff --git a/test/gie/builtins.gie b/test/gie/builtins.gie index fb8a53c3..63e56967 100644 --- a/test/gie/builtins.gie +++ b/test/gie/builtins.gie @@ -124,30 +124,20 @@ expect 3_860_398.3783 5_430_089.0490 roundtrip 100 ------------------------------------------------------------------------------- -Test equatorial aspect of the ellipsoidal azimuthal equidistant. Test data from -Snyder pp. 196-197, table 30. +# Test the azimuthal equidistant modified for Guam. Test data from the EPSG +Guidance Note 7 part 2, September 2016, p. 85 ------------------------------------------------------------------------------- -operation +proj=aeqd +ellps=GRS80 +guam +operation +proj=aeqd +guam +ellps=clrk66 +x_0=50000.00 +y_0=50000.00 + +lon_0=144.74875069444445 +lat_0=13.47246633333333 ------------------------------------------------------------------------------- -tolerance 1 m -accept 0 0 -expect 0.0000 0.0000 -roundtrip 100 -accept 90 90 -except 0.0000 10_001_965.7292 -roundtrip 100 -accept 0 90 -expect 0.0000 10_001_965.7292 -roundtrip 100 -accept 90 90 -expect 0.0000 10_001_965.7292 -roundtrip 100 -accept 45 45 -expect 3548107.5793 5970183.542 -#roundtrip 100 -accept -45 -45 -expect -3548107.5793 -5970183.542 -#roundtrip 100 +tolerance 1 cm +accept 144.635331291666660 13.33903846111111 +expect 37712.48 35242.00 +roundtrip 100 + +direction inverse +accept 37712.48 35242.00 +expect 144.635331291666660 13.33903846111111 ------------------------------------------------------------------------------- Test northern polar aspect of the ellipsoidal azimuthal equidistant. Test data -- cgit v1.2.3