diff options
| -rw-r--r-- | src/projections/aeqd.cpp | 2 | ||||
| -rw-r--r-- | test/gie/builtins.gie | 23 |
2 files changed, 24 insertions, 1 deletions
diff --git a/src/projections/aeqd.cpp b/src/projections/aeqd.cpp index a187edb7..e0fa5c37 100644 --- a/src/projections/aeqd.cpp +++ b/src/projections/aeqd.cpp @@ -151,7 +151,7 @@ oblcon: return xy; } else - xy.x = xy.y = 0.; + return aeqd_e_forward(lp, P); else { xy.y = acos(xy.y); xy.y /= sin(xy.y); diff --git a/test/gie/builtins.gie b/test/gie/builtins.gie index f86b7e79..676086e0 100644 --- a/test/gie/builtins.gie +++ b/test/gie/builtins.gie @@ -136,6 +136,29 @@ accept 45 45 expect 3_860_398.3783 5_430_089.0490 roundtrip 100 + +# Test oblique aeqd with point very close lon_0, lat_0, on a perfect sphere +operation +proj=aeqd +a=6371008.771415 +b=6371008.771415 +lat_0=30.2345 +lon_0=-120.2345 +tolerance 1 mm +accept -120.234501 30.234501 +expect -0.096 0.111 +roundtrip 1 + +accept -120.2345 30.2345 +expect 0.000 0.000 +roundtrip 1 + +# Same on an ellipsoid very close to the sphere +operation +proj=aeqd +a=6371008.771415 +b=6371008.771414 +lat_0=30.2345 +lon_0=-120.2345 +tolerance 1 mm +accept -120.234501 30.234501 +expect -0.096 0.111 +roundtrip 1 + +accept -120.2345 30.2345 +expect 0.000 0.000 +roundtrip 1 + ------------------------------------------------------------------------------- # Test the Modified Azimuthal Equidistant / EPSG 9832. Test data from the EPSG Guidance Note 7 part 2, April 2018, p. 85 |