From 3bd64e8532521d61119df0a6a365a0a32b5c9b15 Mon Sep 17 00:00:00 2001 From: Even Rouault Date: Wed, 15 Sep 2021 05:42:17 +0200 Subject: Fix error in implementation of Inverse ellipsoidal orthographic projection (oblique case) that cause convergence to sometimes fail (fixes #2844) --- src/projections/ortho.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) (limited to 'src') 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; -- cgit v1.2.3