diff options
| author | Even Rouault <even.rouault@spatialys.com> | 2021-09-16 10:20:14 +0200 |
|---|---|---|
| committer | GitHub <noreply@github.com> | 2021-09-16 10:20:14 +0200 |
| commit | 1cedab9d8d48302c61a2ae1f3c139332560f9b73 (patch) | |
| tree | 21307d889df01c298ba7ee15aa35ff0cf9e67906 /src | |
| parent | e6e4d5aa20cfc64b1b41eb6dd6d4516e0ea4807e (diff) | |
| parent | c9b506005129cbd4106a332d33cc9b926291d2b9 (diff) | |
| download | PROJ-1cedab9d8d48302c61a2ae1f3c139332560f9b73.tar.gz PROJ-1cedab9d8d48302c61a2ae1f3c139332560f9b73.zip | |
Merge pull request #2855 from rouault/ortho_e_followup
Inverse ortho ellipsoidal oblique: address a few remarks...
Diffstat (limited to 'src')
| -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 ) |