Fix comparison of GeodeticRefrenceFrame vs DynamicGeodeticReferenceFrame

If comparing a DynamicGeodeticReferenceFrame object and its export to
WKT1, which is a simple DATUM object, currently in non-strict comparison
mode, we'd consider the datum to be equivalent to the dynamic datum, but
not the reverse, which breaks the symmetric property of the
isEquivalentTo() operation. So fix this, to consider both equivalent
whatever the operand order.
(in strict mode, the objects will be considered different of course)

Spotted in the GDAL GeoTIFF CRS reader code:
https://github.com/OSGeo/gdal/blob/f9d48bdcc8c90df20e53b5af5785f1e5d78910db/frmts/gtiff/gt_wkt_srs.cpp#L832

Do same change for vertical datum vs dynamic vertical datum.

1 files changed, 53 insertions, 8 deletions

diff --git a/src/iso19111/datum.cpp b/src/iso19111/datum.cpp
index 7c76061e..758eaec3 100644
--- a/src/iso19111/datum.cpp
+++ b/src/iso19111/datum.cpp
@@ -1414,7 +1414,8 @@ void GeodeticReferenceFrame::_exportToJSON(
 // ---------------------------------------------------------------------------
 
 //! @cond Doxygen_Suppress
-bool GeodeticReferenceFrame::_isEquivalentTo(
+
+bool GeodeticReferenceFrame::isEquivalentToNoExactTypeCheck(
     const util::IComparable *other, util::IComparable::Criterion criterion,
     const io::DatabaseContextPtr &dbContext) const {
     auto otherGRF = dynamic_cast<const GeodeticReferenceFrame *>(other);
@@ -1427,6 +1428,19 @@ bool GeodeticReferenceFrame::_isEquivalentTo(
            ellipsoid()->_isEquivalentTo(otherGRF->ellipsoid().get(), criterion,
                                         dbContext);
 }
+
+// ---------------------------------------------------------------------------
+
+bool GeodeticReferenceFrame::_isEquivalentTo(
+    const util::IComparable *other, util::IComparable::Criterion criterion,
+    const io::DatabaseContextPtr &dbContext) const {
+    if (criterion == Criterion::STRICT &&
+        !util::isOfExactType<GeodeticReferenceFrame>(*other)) {
+        return false;
+    }
+    return isEquivalentToNoExactTypeCheck(other, criterion, dbContext);
+}
+
 //! @endcond
 
 // ---------------------------------------------------------------------------
@@ -1558,11 +1572,20 @@ DynamicGeodeticReferenceFrame::deformationModelName() const {
 bool DynamicGeodeticReferenceFrame::_isEquivalentTo(
     const util::IComparable *other, util::IComparable::Criterion criterion,
     const io::DatabaseContextPtr &dbContext) const {
-    auto otherDGRF = dynamic_cast<const DynamicGeodeticReferenceFrame *>(other);
-    if (otherDGRF == nullptr ||
-        !GeodeticReferenceFrame::_isEquivalentTo(other, criterion, dbContext)) {
+    if (criterion == Criterion::STRICT &&
+        !util::isOfExactType<DynamicGeodeticReferenceFrame>(*other)) {
         return false;
     }
+    if (!GeodeticReferenceFrame::isEquivalentToNoExactTypeCheck(
+            other, criterion, dbContext)) {
+        return false;
+    }
+    auto otherDGRF = dynamic_cast<const DynamicGeodeticReferenceFrame *>(other);
+    if (otherDGRF == nullptr) {
+        // we can go here only if criterion != Criterion::STRICT, and thus
+        // given the above check we can consider the objects equivalent.
+        return true;
+    }
     return frameReferenceEpoch()._isEquivalentTo(
                otherDGRF->frameReferenceEpoch(), criterion) &&
            metadata::Identifier::isEquivalentName(
@@ -2101,7 +2124,7 @@ void VerticalReferenceFrame::_exportToJSON(
 // ---------------------------------------------------------------------------
 
 //! @cond Doxygen_Suppress
-bool VerticalReferenceFrame::_isEquivalentTo(
+bool VerticalReferenceFrame::isEquivalentToNoExactTypeCheck(
     const util::IComparable *other, util::IComparable::Criterion criterion,
     const io::DatabaseContextPtr &dbContext) const {
     auto otherVRF = dynamic_cast<const VerticalReferenceFrame *>(other);
@@ -2121,6 +2144,19 @@ bool VerticalReferenceFrame::_isEquivalentTo(
     }
     return true;
 }
+
+// ---------------------------------------------------------------------------
+
+bool VerticalReferenceFrame::_isEquivalentTo(
+    const util::IComparable *other, util::IComparable::Criterion criterion,
+    const io::DatabaseContextPtr &dbContext) const {
+    if (criterion == Criterion::STRICT &&
+        !util::isOfExactType<VerticalReferenceFrame>(*other)) {
+        return false;
+    }
+    return isEquivalentToNoExactTypeCheck(other, criterion, dbContext);
+}
+
 //! @endcond
 
 // ---------------------------------------------------------------------------
@@ -2195,11 +2231,20 @@ DynamicVerticalReferenceFrame::deformationModelName() const {
 bool DynamicVerticalReferenceFrame::_isEquivalentTo(
     const util::IComparable *other, util::IComparable::Criterion criterion,
     const io::DatabaseContextPtr &dbContext) const {
-    auto otherDGRF = dynamic_cast<const DynamicVerticalReferenceFrame *>(other);
-    if (otherDGRF == nullptr ||
-        !VerticalReferenceFrame::_isEquivalentTo(other, criterion, dbContext)) {
+    if (criterion == Criterion::STRICT &&
+        !util::isOfExactType<DynamicVerticalReferenceFrame>(*other)) {
         return false;
     }
+    if (!VerticalReferenceFrame::isEquivalentToNoExactTypeCheck(
+            other, criterion, dbContext)) {
+        return false;
+    }
+    auto otherDGRF = dynamic_cast<const DynamicVerticalReferenceFrame *>(other);
+    if (otherDGRF == nullptr) {
+        // we can go here only if criterion != Criterion::STRICT, and thus
+        // given the above check we can consider the objects equivalent.
+        return true;
+    }
     return frameReferenceEpoch()._isEquivalentTo(
                otherDGRF->frameReferenceEpoch(), criterion) &&
            metadata::Identifier::isEquivalentName(