diff options
| author | Even Rouault <even.rouault@spatialys.com> | 2020-09-26 14:40:29 +0200 |
|---|---|---|
| committer | Even Rouault <even.rouault@spatialys.com> | 2020-09-26 18:55:25 +0200 |
| commit | 3d88b6fc89f95803bfd4d59c47eef1c05bda710c (patch) | |
| tree | a67c595277c9cdce18dac60b771fdaad7ff089eb /test/gie | |
| parent | 2e104e092578347de11208e9a3a80a3bf711265d (diff) | |
| download | PROJ-3d88b6fc89f95803bfd4d59c47eef1c05bda710c.tar.gz PROJ-3d88b6fc89f95803bfd4d59c47eef1c05bda710c.zip | |
Ortho: add visibility condition for ellipsoidal case. Credits to @cffk
Diffstat (limited to 'test/gie')
| -rw-r--r-- | test/gie/builtins.gie | 41 |
1 files changed, 41 insertions, 0 deletions
diff --git a/test/gie/builtins.gie b/test/gie/builtins.gie index 48afc00a..9d5b6644 100644 --- a/test/gie/builtins.gie +++ b/test/gie/builtins.gie @@ -4434,13 +4434,23 @@ 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 89.99 0 +expect 6378136.9029 0 +roundtrip 1 + +accept -89.99 0 +expect -6378136.9029 0 +roundtrip 1 + accept 0 89.99 expect 0 6356752.2167 roundtrip 1 @@ -4449,12 +4459,33 @@ accept 0 -89.99 expect 0 -6356752.2167 roundtrip 1 +# Point not visible from the projection plane +accept 90.00001 0 +expect failure errno tolerance_condition + +# Point not visible from the projection plane +accept -90.00001 0 +expect failure errno tolerance_condition + +# Consistant with WGS84 semi-major axis +# The inverse transformation doesn't converge due to properties of the projection +accept 90 0 +expect 6378137 0 + +accept -90 0 +expect -6378137 0 + # 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 +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 @@ -4465,12 +4496,18 @@ accept 1 89.999 expect 1.9493 -111.6770 roundtrip 1 +# Point not visible from the projection plane +accept 0 -0.0000001 +expect failure errno tolerance_condition + # 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 @@ -4481,6 +4518,10 @@ accept 1 -89.999 expect 1.9493 111.6770 roundtrip 1 +# Point not visible from the projection plane +accept 0 0.0000001 +expect failure errno tolerance_condition + # Consistant with WGS84 semi-major axis # The inverse transformation doesn't converge due to properties of the projection accept 0 0 |