Merge branch '2051_CRS_equality_tests' into 'develop'

Fix spatial reference equality tests

Closes #2051 and #2055

See merge request orfeotoolbox/otb!720

Modules/Learning/Sampling/include/otbImageSampleExtractorFilter.hxx

@@ -117,6 +117,9 @@ void PersistentImageSampleExtractorFilter<TInputImage>::GenerateOutputInformatio
 #if GDAL_VERSION_NUM >= 3000000 // importFromWkt is const-correct in GDAL 3
     const char* projWktCstr = projectionRefWkt.c_str();
     OGRErr      err         = imgSRS.importFromWkt(&projWktCstr);
+    // Use the same mapping strategy as the one in the datasource.
+    imgSRS.SetAxisMappingStrategy(
+      this->GetOGRData()->GetLayer(this->GetLayerIndex()).GetSpatialRef()->GetAxisMappingStrategy ());
 #else
     const char* projWktCstr    = projectionRefWkt.c_str();
     char**      projWktPointer = const_cast<char**>(&projWktCstr);

Modules/Segmentation/Conversion/include/otbPersistentImageToOGRLayerFilter.hxx

@@ -89,6 +89,14 @@ void PersistentImageToOGRLayerFilter<TImage>::Initialize()
 
   oSRSESRI.morphToESRI();
   oSRSESRI.morphFromESRI();
+
+#if GDAL_VERSION_NUM >= 3000000 // importFromWkt is const-correct in GDAL 3
+    // Use the same mapping strategy as the one in the datasource.
+    auto mappingStrategy = m_OGRLayer.GetSpatialRef()->GetAxisMappingStrategy();
+    oSRS.SetAxisMappingStrategy(mappingStrategy);
+    oSRSESRI.SetAxisMappingStrategy(mappingStrategy);
+#endif
+
   if (m_OGRLayer.GetSpatialRef() && (!oSRS.IsSame(m_OGRLayer.GetSpatialRef()) && !oSRSESRI.IsSame(m_OGRLayer.GetSpatialRef())))
   {
     if ((oSRS.Validate() != OGRERR_NONE) && (oSRSESRI.Validate() != OGRERR_NONE))