Merge pull request #1655 from rouault/aeqd_obliq_sphere

aeqd: for spherical forward path, go to higher precision ellipsoidal case when the point coordinates are super close to the origin (fixes #1654)

1 files changed, 1 insertions, 1 deletions

diff --git a/src/projections/aeqd.cpp b/src/projections/aeqd.cpp
index a187edb7..e0fa5c37 100644
--- a/src/projections/aeqd.cpp
+++ b/src/projections/aeqd.cpp
@@ -151,7 +151,7 @@ oblcon:
                 return xy;
             }
             else
-                xy.x = xy.y = 0.;
+                return aeqd_e_forward(lp, P);
         else {
             xy.y = acos(xy.y);
             xy.y /= sin(xy.y);