aboutsummaryrefslogtreecommitdiff
path: root/test/gie
diff options
context:
space:
mode:
authorEven Rouault <even.rouault@spatialys.com>2020-09-25 22:53:38 +0200
committerEven Rouault <even.rouault@spatialys.com>2020-09-26 11:09:23 +0200
commit2e104e092578347de11208e9a3a80a3bf711265d (patch)
tree821679ce3db494a531c7727ad4b5db926b005370 /test/gie
parent485509f3b3f5133f0bb0db5ef63e932615fa2f2e (diff)
downloadPROJ-2e104e092578347de11208e9a3a80a3bf711265d.tar.gz
PROJ-2e104e092578347de11208e9a3a80a3bf711265d.zip
Implement ellipsoidal formulation of +proj=ortho (fixes #397)
- Map ESRI 'Local' to +proj=ortho when Scale_Factor = 1 and Azimuth = 0 - Map ESRI 'Orthographic' to a PROJ WKT2 'Orthographic (Spherical)' which maps to +proj=ortho +f=0 to froce spherical evaluation
Diffstat (limited to 'test/gie')
-rw-r--r--test/gie/builtins.gie68
1 files changed, 67 insertions, 1 deletions
diff --git a/test/gie/builtins.gie b/test/gie/builtins.gie
index f755526e..48afc00a 100644
--- a/test/gie/builtins.gie
+++ b/test/gie/builtins.gie
@@ -4266,7 +4266,7 @@ expect -223374.577355253 -111701.072127637
===============================================================================
# Orthographic
-# Azi, Sph.
+# Azi, Sph&Ell.
===============================================================================
-------------------------------------------------------------------------------
@@ -4421,6 +4421,72 @@ expect failure errno tolerance_condition
accept 0.70710678118 0.7071067812
expect 45 0
+
+-------------------------------------------------------------------------------
+# Test the ellipsoidal formulation
+
+# Test data from Guidance Note 7 part 2, March 2020, p. 90
+-------------------------------------------------------------------------------
+operation +proj=ortho +ellps=WGS84 +lat_0=55 +lon_0=5
+-------------------------------------------------------------------------------
+tolerance 1 mm
+accept 2.12955 53.80939444444444
+expect -189011.711 -128640.567
+roundtrip 1
+
+# Equatorial
+operation +proj=ortho +ellps=WGS84
+tolerance 0.1 mm
+accept 0 0
+expect 0 0
+roundtrip 1
+
+accept 0 89.99
+expect 0 6356752.2167
+roundtrip 1
+
+accept 0 -89.99
+expect 0 -6356752.2167
+roundtrip 1
+
+# Consistant with WGS84 semi-minor axis
+# The inverse transformation doesn't converge due to properties of the projection
+accept 0 90
+expect 0 6356752.3142
+
+# North pole tests
+operation +proj=ortho +ellps=WGS84 +lat_0=90
+tolerance 0.1 mm
+accept 0 90
+expect 0 0
+roundtrip 1
+
+accept 1 89.999
+expect 1.9493 -111.6770
+roundtrip 1
+
+# Consistant with WGS84 semi-major axis
+# The inverse transformation doesn't converge due to properties of the projection
+accept 0 0
+expect 0 -6378137
+
+# South pole tests
+operation +proj=ortho +ellps=WGS84 +lat_0=-90
+tolerance 0.1 mm
+accept 0 -90
+expect 0 0
+roundtrip 1
+
+accept 1 -89.999
+expect 1.9493 111.6770
+roundtrip 1
+
+# Consistant with WGS84 semi-major axis
+# The inverse transformation doesn't converge due to properties of the projection
+accept 0 0
+expect 0 6378137
+
+
===============================================================================
# Perspective Conic
# Conic, Sph