diff options
| author | Even Rouault <even.rouault@spatialys.com> | 2021-09-15 05:42:17 +0200 |
|---|---|---|
| committer | Even Rouault <even.rouault@spatialys.com> | 2021-09-15 06:41:39 +0200 |
| commit | 3bd64e8532521d61119df0a6a365a0a32b5c9b15 (patch) | |
| tree | d959f1c027722e532cc0c088de35e95dccb812d1 /src/projections/ortho.cpp | |
| parent | 47fb85bdf0e45bf5660e0b1dbde7b76f8824a87a (diff) | |
| download | PROJ-3bd64e8532521d61119df0a6a365a0a32b5c9b15.tar.gz PROJ-3bd64e8532521d61119df0a6a365a0a32b5c9b15.zip | |
Fix error in implementation of Inverse ellipsoidal orthographic projection (oblique case) that cause convergence to sometimes fail (fixes #2844)
Diffstat (limited to 'src/projections/ortho.cpp')
| -rw-r--r-- | src/projections/ortho.cpp | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/src/projections/ortho.cpp b/src/projections/ortho.cpp index 908f283d..d2d113ef 100644 --- a/src/projections/ortho.cpp +++ b/src/projections/ortho.cpp @@ -254,7 +254,7 @@ static PJ_LP ortho_e_inverse (PJ_XY xy, PJ *P) { /* Ellipsoidal, inver const double J11 = -rho * sinphi * sinlam; const double J12 = nu * cosphi * coslam; const double J21 = rho * (cosphi * Q->cosph0 + sinphi * Q->sinph0 * coslam); - const double J22 = nu * Q->sinph0 * Q->cosph0 * sinlam; + const double J22 = nu * Q->sinph0 * cosphi * sinlam; const double D = J11 * J22 - J12 * J21; const double dx = xy.x - xy_new.x; const double dy = xy.y - xy_new.y; |