TimeSeries/s2theia.py

@@ -176,7 +176,12 @@ class S2TheiaTilePipeline:
             er.SetParameterString('in', bnd)
             er.SetParameterString('mode', 'fit')
             er.SetParameterString('mode.fit.vect', roi)
-            er.Execute()
+            try:
+                er.Execute()
+            except Exception as e:
+                if not 'Extraction Region not consistent with output image' in str(e):
+                    print(e)
+                return False
             arr = er.GetImageAsNumpyArray('out')
             if (np.sum(arr != cls.NDT) / (arr.shape[0]*arr.shape[1])) <= min_surf:
                 return False
@@ -259,6 +264,9 @@ class S2TheiaTilePipeline:
     def is_empty(self):
         return len(self.image_list) == 0
 
+    def is_empty(self):
+        return len(self.image_list) == 0
+
     def reset(self):
         self.pipe = []
         self.files = []
@@ -797,6 +805,7 @@ class S2TheiaPipeline:
         self.temp_fld = temp_fld
         self.input_date_interval = input_date_interval
         self.max_clouds_percentage = max_clouds_percentage
+        self.roi = roi
         self.tile_list = set()
         self.roi = roi
         img_list = [os.path.abspath(x) for x in glob.glob(os.path.join(self.folder, self.S2TilePipeline.PTRN_dir))
@@ -846,6 +855,7 @@ class S2TheiaPipeline:
             for t in self.tiles:
                 t.preprocess()
                 if self.roi is not None:
+                    print('Clipping')
                     t.clip(self.roi)
                 if align:
                     if type(align_to) == str and os.path.exists(align_to):