Fixed invalid polygons. Fixed wrong call within run_gms.sh

bin/run_gms.sh

@@ -5,4 +5,4 @@
 #       variable points to the correct Python interpreter.
 
 #ipython --colors='NoColor' run_gms.py jobid "$@"     # NoColor must be active for ipython because gwt-log cannot interpret ANSI console colors
-python run_gms.py jobid "$@"
+run_gms.py jobid "$@"

gms_preprocessing/misc/helper_functions.py

@@ -473,13 +473,20 @@ def scene_ID_to_shapelyPolygon(scene_ID):
     """
     Returns a LonLat shapely.Polygon() object corresponding to the given scene_ID.
     """
-
-    return Polygon(reorder_CornerLonLat(sceneID_to_trueDataCornerLonLat(scene_ID)))
+    poly = Polygon(reorder_CornerLonLat(sceneID_to_trueDataCornerLonLat(scene_ID)))
+    if not poly.is_valid:
+        poly = poly.buffer(0)
+        assert poly.is_valid
+    return poly
 
 
 def CornerLonLat_to_shapelyPoly(CornerLonLat):
     """Returns a shapely.Polygon() object based on the given coordinate list. """
-    return Polygon(reorder_CornerLonLat(CornerLonLat))
+    poly = Polygon(reorder_CornerLonLat(CornerLonLat))
+    if not poly.is_valid:
+        poly = poly.buffer(0)
+        assert poly.is_valid
+    return poly
 
 
 def find_in_xml_root(namespace, xml_root, branch, *branches, findall=None):

gms_preprocessing/misc/spatial_index_mediator.py

@@ -324,4 +324,10 @@ class Scene:
         self.bounds = bounds
         tempList = list(bounds) + [None] * 2
         self.coordsLonLat = [tempList[n:n + 2] for n in range(0, len(bounds), 2)]
-        self.polyLonLat = Polygon(self.coordsLonLat)
+
+        # set validated (!) polygon
+        poly = Polygon(self.coordsLonLat)
+        if not poly.is_valid:
+            poly = poly.buffer(0)
+            assert poly.is_valid
+        self.polyLonLat = poly