diff options
| author | Even Rouault <even.rouault@spatialys.com> | 2021-09-15 18:08:05 +0200 |
|---|---|---|
| committer | Even Rouault <even.rouault@spatialys.com> | 2021-09-15 18:08:05 +0200 |
| commit | c9b506005129cbd4106a332d33cc9b926291d2b9 (patch) | |
| tree | 21307d889df01c298ba7ee15aa35ff0cf9e67906 /src/projections | |
| parent | e6e4d5aa20cfc64b1b41eb6dd6d4516e0ea4807e (diff) | |
| download | PROJ-c9b506005129cbd4106a332d33cc9b926291d2b9.tar.gz PROJ-c9b506005129cbd4106a332d33cc9b926291d2b9.zip | |
Inverse ortho ellipsoidal oblique: address a few remarks from https://github.com/OSGeo/PROJ/issues/2844#issuecomment-920138371
Diffstat (limited to 'src/projections')
| -rw-r--r-- | src/projections/ortho.cpp | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/src/projections/ortho.cpp b/src/projections/ortho.cpp index c334f398..9dae0980 100644 --- a/src/projections/ortho.cpp +++ b/src/projections/ortho.cpp @@ -263,13 +263,13 @@ static PJ_LP ortho_e_inverse (PJ_XY xy, PJ *P) { /* Ellipsoidal, inver lp.phi += dphi; if( lp.phi > M_PI_2) { - lp.phi = M_PI_2 - (M_PI_2 - lp.phi); - lp.lam = adjlon(lp.lam + 180); + lp.phi = M_PI_2 - (lp.phi - M_PI_2); + lp.lam = adjlon(lp.lam + M_PI); } else if( lp.phi < -M_PI_2) { lp.phi = -M_PI_2 + (-M_PI_2 - lp.phi); - lp.lam = adjlon(lp.lam + 180); + lp.lam = adjlon(lp.lam + M_PI); } lp.lam += dlam; if( fabs(dphi) < 1e-12 && fabs(dlam) < 1e-12 ) |