Add S2 projection (#2749)

1 files changed, 167 insertions, 0 deletions

diff --git a/test/gie/builtins.gie b/test/gie/builtins.gie
index 1f3824c8..b71254f4 100644
--- a/test/gie/builtins.gie
+++ b/test/gie/builtins.gie
@@ -5273,6 +5273,173 @@ expect  -223368.098302014 111769.110486991
 accept  -2 -1
 expect  -223368.098302014 -111769.110486991
 
+===============================================================================
+# S2
+# Input lats are converted from nice spherical values to
+# messy ellipsoidal lats e.g. 45 vs 45.1924232....
+===============================================================================
+
+-------------------------------------------------------------------------------
+operation +proj=s2   +ellps=WGS84 +lat_0=0 +lon_0=0 +UVtoST=linear
+-------------------------------------------------------------------------------
+tolerance 0.1 mm
+accept 0 0
+expect 0.5 0.5
+accept 0 45.19242321598196
+expect 0.5 1
+accept 0 -45.19242321598196
+expect 0.5 0
+accept -45 0
+expect 0 0.5
+accept 45 0
+expect 1 0.5
+accept  -45 -35.446011426401625
+expect  0 0
+accept  45 -35.446011426401625
+expect  1 0
+accept  45 35.446011426401625
+expect  1 1
+accept  -45 35.446011426401625
+expect  0 1
+accept 20 20.124006563576454
+expect 0.6819851171331012 0.6936645165744716
+accept 20 -20.124006563576454
+expect 0.6819851171331012 0.3063354834255284
+accept -20 -20.124006563576454
+expect 0.31801488286689883 0.3063354834255284
+accept -20 20.124006563576454
+expect 0.31801488286689883 0.6936645165744716
+
+direction inverse
+accept  0.5 0.5
+expect  0 0
+accept 0.5 1
+expect 0 45.19242321598196
+accept 0.5 0
+expect 0 -45.19242321598196
+accept 0 0.5
+expect -45 0
+accept 1 0.5
+expect 45 0
+accept 0 0
+expect -45 -35.446011426401625
+accept 1 0 
+expect 45 -35.446011426401625
+accept 1 1
+expect 45 35.446011426401625
+accept 0 1
+expect -45 35.446011426401625
+accept 0.6819851171331012 0.6936645165744716
+expect 20 20.124006563576454
+accept 0.6819851171331012 0.3063354834255284
+expect 20 -20.124006563576454
+accept 0.31801488286689883 0.3063354834255284
+expect -20 -20.124006563576454
+accept 0.31801488286689883 0.6936645165744716
+expect -20 20.124006563576454
+
+-------------------------------------------------------------------------------
+operation +proj=s2   +ellps=WGS84 +lat_0=0 +lon_0=90 +UVtoST=quadratic
+-------------------------------------------------------------------------------
+tolerance 0.1 mm
+accept 90 0
+expect 0.5 0.5
+accept 70 20.124006563576454
+expect 0.27682804555233764 0.7351848576118168
+accept 110 20.124006563576454
+expect 0.7231719544476624 0.7351848576118168
+
+direction inverse
+accept 0.5 0.5
+expect 90 0
+accept 0.27682804555233764 0.7351848576118168
+expect 70 20.124006563576454
+accept 0.7231719544476624 0.7351848576118168
+expect 110 20.124006563576454
+
+-------------------------------------------------------------------------------
+operation +proj=s2   +ellps=WGS84 +lat_0=90 +UVtoST=tangent
+-------------------------------------------------------------------------------
+tolerance 0.1 mm
+accept 0 90
+expect 0.5 0.5
+accept 20 70.12337013762532
+expect 0.29020309743436806 0.4211558922141421
+accept -20 70.12337013762532
+expect 0.29020309743436806 0.5788441077858579
+
+direction inverse
+accept 0.5 0.5
+expect 0 90
+accept 0.29020309743436806 0.4211558922141421
+expect 20 70.12337013762532
+accept 0.29020309743436806 0.5788441077858579
+expect -20 70.12337013762532
+
+-------------------------------------------------------------------------------
+operation +proj=s2   +ellps=WGS84 +lat_0=0 +lon_0=180 +UVtoST=none
+-------------------------------------------------------------------------------
+tolerance 0.1 mm
+accept 180 0
+expect 0 0
+accept 160 20.124006563576454
+expect -0.3873290331489431 -0.3639702342662023
+accept -160 20.124006563576454
+expect -0.3873290331489431 0.3639702342662023
+
+direction inverse
+accept 0 0
+expect 180 0
+accept -0.3873290331489431 -0.3639702342662023
+expect 160 20.124006563576454
+accept -0.3873290331489431 0.3639702342662023
+expect -160 20.124006563576454
+
+-------------------------------------------------------------------------------
+operation +proj=s2   +ellps=WGS84 +lat_0=0 +lon_0=-90
+-------------------------------------------------------------------------------
+tolerance 0.1 mm
+accept -90 0
+expect 0.5 0.5
+accept -70 20.124006563576454
+expect 0.26481514238818316 0.7231719544476624
+accept -110 20.124006563576454
+expect 0.26481514238818316 0.27682804555233764
+
+direction inverse
+accept 0.5 0.5
+expect -90 0
+accept 0.26481514238818316 0.7231719544476624
+expect -70 20.124006563576454
+accept 0.26481514238818316 0.27682804555233764
+expect -110 20.124006563576454
+
+-------------------------------------------------------------------------------
+operation +proj=s2   +ellps=WGS84 +lat_0=-90 +UVtoST=linear
+-------------------------------------------------------------------------------
+tolerance 0.1 mm
+accept 0 -90
+expect 0.5 0.5
+accept 20 -70.12337013762533
+expect 0.5622425758450019 0.6710100716628344
+accept -20 -70.12337013762533
+expect 0.4377574241549981 0.6710100716628344
+
+direction inverse
+accept 0.5 0.5
+expect 0 -90
+accept 0.5622425758450019 0.6710100716628344
+expect 20 -70.12337013762533
+accept 0.4377574241549981 0.6710100716628344
+expect -20 -70.12337013762533
+
+-------------------------------------------------------------------------------
+operation +proj=s2   +ellps=WGS84 +lat_0=0 +lon_0=0 +UVtoST=invalid
+-------------------------------------------------------------------------------
+tolerance 0.1 mm
+
+accept 0 0
+expect failure errno invalid_op_illegal_arg_value
 
 ===============================================================================
 # Sinusoidal (Sanson-Flamsteed)