diff options
| author | Even Rouault <even.rouault@spatialys.com> | 2021-09-08 12:50:37 +0200 |
|---|---|---|
| committer | Even Rouault <even.rouault@spatialys.com> | 2021-09-08 17:05:45 +0200 |
| commit | 85733181ee7c2777139f5d1db94f2beabb737e96 (patch) | |
| tree | 2936c5f142c390bae34388925fbed74d0ad520a6 /test/cli | |
| parent | 5527b10ed140e20fac8e183317514fd59e4c8b99 (diff) | |
| download | PROJ-85733181ee7c2777139f5d1db94f2beabb737e96.tar.gz PROJ-85733181ee7c2777139f5d1db94f2beabb737e96.zip | |
createOperations(): deal with spherical planetocentric geodetic CRS
This also fixes conversion between geocentric latlong and geodetic latlong
with cs2cs. This was dealt with in PR 1093, but in the wrong direction
(the geocentric latitude must be <= in absolute value to the geodetic one)
The issue here was linked to the semantics of the +geoc specifier, which
affects the semantics of the input coordinates in the forward direction
(+geoc means that the input coordinate is is a geocentric latitude),
which defeats the logic of doing A to B by using the inverse path of A
and the forward path of B.
Diffstat (limited to 'test/cli')
| -rwxr-xr-x | test/cli/testvarious | 12 | ||||
| -rw-r--r-- | test/cli/tv_out.dist | 4 |
2 files changed, 8 insertions, 8 deletions
diff --git a/test/cli/testvarious b/test/cli/testvarious index 6e9cd43c..80340369 100755 --- a/test/cli/testvarious +++ b/test/cli/testvarious @@ -144,10 +144,10 @@ $EXE +proj=geocent +datum=WGS84 \ EOF # echo "#############################################################" >> ${OUT} -echo Test conversion from geocentric latlong to geodetic latlong >> ${OUT} +echo Test conversion from geodetic latlong to geocentric latlong >> ${OUT} # -$EXE +proj=latlong +datum=WGS84 +geoc \ - +to +proj=latlong +datum=WGS84 \ +$EXE +proj=latlong +datum=WGS84 \ + +to +proj=latlong +datum=WGS84 +geoc \ -E >>${OUT} <<EOF 0d00'00.000"W 0d00'00.000"N 0.0 79d00'00.000"W 45d00'00.000"N 0.0 @@ -156,10 +156,10 @@ $EXE +proj=latlong +datum=WGS84 +geoc \ EOF # echo "#############################################################" >> ${OUT} -echo Test conversion from geodetic latlong to geocentric latlong >> ${OUT} +echo Test conversion from geocentric latlong to geodetic latlong >> ${OUT} # -$EXE +proj=latlong +datum=WGS84 \ - +to +proj=latlong +datum=WGS84 +geoc \ +$EXE +proj=latlong +datum=WGS84 +geoc \ + +to +proj=latlong +datum=WGS84 \ -E >>${OUT} <<EOF 0d00'00.000"W 0d00'00.000"N 0.0 79d00'00.000"W 44d48'27.276"N 0.000 diff --git a/test/cli/tv_out.dist b/test/cli/tv_out.dist index 9bb85933..06aa1a2d 100644 --- a/test/cli/tv_out.dist +++ b/test/cli/tv_out.dist @@ -43,13 +43,13 @@ Test geocentric x/y/z consumption. 861996.98 -4434590.01 4487348.41 79dW 45dN 0.001 0.00 -0.00 6356752.31 0dE 90dN -0.004 ############################################################# -Test conversion from geocentric latlong to geodetic latlong +Test conversion from geodetic latlong to geocentric latlong 0d00'00.000"W 0d00'00.000"N 0.0 0dE 0dN 0.000 79d00'00.000"W 45d00'00.000"N 0.0 79dW 44d48'27.276"N 0.000 12d00'00.000"W 45d00'00.000"N 0.0 12dW 44d48'27.276"N 0.000 0d00'00.000"W 90d00'00.000"N 0.0 0dE 90dN 0.000 ############################################################# -Test conversion from geodetic latlong to geocentric latlong +Test conversion from geocentric latlong to geodetic latlong 0d00'00.000"W 0d00'00.000"N 0.0 0dE 0dN 0.000 79d00'00.000"W 44d48'27.276"N 0.000 79dW 45dN 0.000 12d00'00.000"W 44d48'27.276"N 0.0 12dW 45dN 0.000 |