aboutsummaryrefslogtreecommitdiff
path: root/test/gie
diff options
context:
space:
mode:
authorEven Rouault <even.rouault@spatialys.com>2021-08-16 18:03:39 +0200
committerEven Rouault <even.rouault@spatialys.com>2021-08-16 18:03:39 +0200
commited16a27e1f783d1b929aaa9fb9f41bdd2b63fbc1 (patch)
treebab66c6165c8e9a4db913e50a4b4f5ef3c4e27ee /test/gie
parent9190f2b7b0dd7c9b223c171c943963f3bb37872f (diff)
downloadPROJ-ed16a27e1f783d1b929aaa9fb9f41bdd2b63fbc1.tar.gz
PROJ-ed16a27e1f783d1b929aaa9fb9f41bdd2b63fbc1.zip
test: more testing of Polar Stereographic variants
Diffstat (limited to 'test/gie')
-rw-r--r--test/gie/builtins.gie45
1 files changed, 45 insertions, 0 deletions
diff --git a/test/gie/builtins.gie b/test/gie/builtins.gie
index 70d99a01..9e1d3339 100644
--- a/test/gie/builtins.gie
+++ b/test/gie/builtins.gie
@@ -5615,6 +5615,7 @@ accept 20 70
expect 748315.3282 -2055979.4669
roundtrip 1
+# Polar Stereographic Variant B, ellipsoidal
-------------------------------------------------------------------------------
operation +proj=stere +ellps=GRS80 +lat_0=-90 +lat_ts=-70
-------------------------------------------------------------------------------
@@ -5628,6 +5629,50 @@ accept 20 -70
expect 748315.3282 2055979.4669
roundtrip 1
+# Polar Stereographic Variant A, ellipsoidal
+-------------------------------------------------------------------------------
+operation +proj=stere +ellps=GRS80 +lat_0=-90 +k_0=0.97
+-------------------------------------------------------------------------------
+tolerance 1e-15m
+accept 0 -90
+expect 0 0
+roundtrip 1
+
+tolerance 0.1 mm
+accept 20 -70
+expect 748424.7446 2056280.0858
+roundtrip 1
+
+# Polar Stereographic Variant B, spherical
+-------------------------------------------------------------------------------
+operation +proj=stere +R=6378137 +lat_0=-90 +lat_ts=-70
+-------------------------------------------------------------------------------
+tolerance 1e-15m
+accept 0 -90
+expect 0 0
+roundtrip 1
+
+tolerance 0.1 mm
+accept 20 -70
+expect 746100.2968 2049893.7182
+roundtrip 1
+
+
+# Polar Stereographic Variant A, spherical
+# k_0 = (1 + math.sin(abs(lat_ts=-70) / 180. * math.pi)) / 2.0
+-------------------------------------------------------------------------------
+operation +proj=stere +R=6378137 +lat_0=-90 +k_0=0.9698463103929542
+-------------------------------------------------------------------------------
+tolerance 1e-15m
+accept 0 -90
+expect 0 0
+roundtrip 1
+
+tolerance 0.1 mm
+accept 20 -70
+expect 746100.2968 2049893.7182
+roundtrip 1
+
===============================================================================
# Oblique Stereographic Alternative
# Azimuthal, Sph&Ell