sen2meta.py: Cope with non-compatible API change in GDAL 3, in which …

…long/lat can be swapped.
ubarsc · Nov 29, 2019 · 3071ca2 · 3071ca2
1 parent 75fdec7
commit 3071ca2
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/fmask/sen2meta.py b/fmask/sen2meta.py
@@ -175,6 +175,10 @@ def getCtrLongLat(self):
         srUTM.ImportFromEPSG(int(self.epsg))
         srLL = osr.SpatialReference()
         srLL.ImportFromEPSG(4326)
+        if hasattr(srLL, 'SetAxisMappingStrategy'):
+            # We are in GDAL >= 3, so guard against axis swapping
+            srLL.SetAxisMappingStrategy(osr.OAMS_TRADITIONAL_GIS_ORDER)
+            srUTM.SetAxisMappingStrategy(osr.OAMS_TRADITIONAL_GIS_ORDER)
         tr = osr.CoordinateTransformation(srUTM, srLL)
         (longitude, latitude, z) = tr.TransformPoint(ctrX, ctrY)
         return (longitude, latitude)