From 13625a0c3f5a8a1b3a6ac63b4f49234a3722f62c Mon Sep 17 00:00:00 2001
From: Even Rouault <even.rouault@spatialys.com>
Date: Thu, 13 Dec 2018 16:23:01 +0100
Subject: Tune behaviour for LAEA vs LAEA Spherical or EQC vs EQC Spherical

---
 src/coordinateoperation.cpp | 27 ++++++++++++++++++++++++++-
 1 file changed, 26 insertions(+), 1 deletion(-)

(limited to 'src/coordinateoperation.cpp')

diff --git a/src/coordinateoperation.cpp b/src/coordinateoperation.cpp
index a1f305bf..ac170bbf 100644
--- a/src/coordinateoperation.cpp
+++ b/src/coordinateoperation.cpp
@@ -1503,10 +1503,35 @@ bool SingleOperation::_isEquivalentTo(const util::IComparable *other,
     const int methodEPSGCode = d->method_->getEPSGCode();
     const int otherMethodEPSGCode = otherSO->d->method_->getEPSGCode();
 
-    const bool equivalentMethods =
+    bool equivalentMethods =
         (criterion == util::IComparable::Criterion::EQUIVALENT &&
          methodEPSGCode != 0 && methodEPSGCode == otherMethodEPSGCode) ||
         d->method_->_isEquivalentTo(otherSO->d->method_.get(), criterion);
+    if (!equivalentMethods &&
+        criterion == util::IComparable::Criterion::EQUIVALENT) {
+        if ((methodEPSGCode == EPSG_CODE_METHOD_LAMBERT_AZIMUTHAL_EQUAL_AREA &&
+             otherMethodEPSGCode ==
+                 EPSG_CODE_METHOD_LAMBERT_AZIMUTHAL_EQUAL_AREA_SPHERICAL) ||
+            (otherMethodEPSGCode ==
+                 EPSG_CODE_METHOD_LAMBERT_AZIMUTHAL_EQUAL_AREA &&
+             methodEPSGCode ==
+                 EPSG_CODE_METHOD_LAMBERT_AZIMUTHAL_EQUAL_AREA_SPHERICAL) ||
+            (methodEPSGCode == EPSG_CODE_METHOD_EQUIDISTANT_CYLINDRICAL &&
+             otherMethodEPSGCode ==
+                 EPSG_CODE_METHOD_EQUIDISTANT_CYLINDRICAL_SPHERICAL) ||
+            (otherMethodEPSGCode == EPSG_CODE_METHOD_EQUIDISTANT_CYLINDRICAL &&
+             methodEPSGCode ==
+                 EPSG_CODE_METHOD_EQUIDISTANT_CYLINDRICAL_SPHERICAL)) {
+            auto geodCRS =
+                dynamic_cast<const crs::GeodeticCRS *>(sourceCRS().get());
+            auto otherGeodCRS = dynamic_cast<const crs::GeodeticCRS *>(
+                otherSO->sourceCRS().get());
+            if (geodCRS && otherGeodCRS && geodCRS->ellipsoid()->isSphere() &&
+                otherGeodCRS->ellipsoid()->isSphere()) {
+                equivalentMethods = true;
+            }
+        }
+    }
 
     if (!equivalentMethods) {
         if (criterion == util::IComparable::Criterion::EQUIVALENT) {
-- 
cgit v1.2.3