diff options
| -rw-r--r-- | src/projections/ortho.cpp | 2 | ||||
| -rw-r--r-- | test/gie/builtins.gie | 10 |
2 files changed, 11 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; diff --git a/test/gie/builtins.gie b/test/gie/builtins.gie index 9e1d3339..9cebec84 100644 --- a/test/gie/builtins.gie +++ b/test/gie/builtins.gie @@ -4531,6 +4531,16 @@ tolerance 0.1 mm accept 0 -6343600 expect 0 -59.966377950099655436 +------------------------------------------------------------------------------- +# Oblique +# Test case from https://github.com/OSGeo/PROJ/issues/2844 +------------------------------------------------------------------------------- +operation +proj=ortho +ellps=WGS84 +lon_0=23.0 +lat_0=37.0 +tolerance 0.1 mm +accept 120.0 84.0 +expect 663929.0678 5118338.2423 +roundtrip 1 + ------------------------------------------------------------------------------- # Equatorial |