updated function if l2a.path does not exist

sen2chain/time_series.py

@@ -206,12 +206,10 @@ class TimeSeries:
         if cloud_proba_path:
             stats_cld_prb = zonal_stats(geom, str(cloud_proba_path), stats=["mean"])[0]
             stats['cldprb'] = stats_cld_prb["mean"]
-            
         logger.info(stats)
         
         
         
-        
         # mettre ici le cloud mask !
         
         
@@ -244,7 +242,11 @@ class TimeSeries:
                             indice_product = IndiceProduct(prod.identifier)
                             prod_path = str(indice_product.path)
                             prod_path_unmasked = prod_path.replace("_CM001", "")
-                            prod_path_cloud_proba = L2aProduct(indice_product.l2a).msk_cldprb_20m
+                            l2a = L2aProduct(indice_product.l2a)
+                            if l2a.path.exists():
+                                prod_path_cloud_proba = l2a.msk_cldprb_20m
+                            else:
+                                prod_path_cloud_proba = None
                             #~ prod_path = tile_indice_path / prod.identifier[:(-12 - len(indice))] / prod.identifier
                             #~ cloud_path = tile_obj.paths["l2a"] / (prod.identifier[:(-12 - len(indice))] + "_CLOUD_MASK.jp2")
                             #~ prod_path_unmasked = tile_indice_path / prod.identifier[:(-12 - len(indice))] / (prod.identifier[:-11] + '.jp2')
@@ -282,8 +284,11 @@ class TimeSeries:
         indice_product = IndiceProduct(proc_item[0].identifier)
         prod_path = str(indice_product.path)
         prod_path_unmasked = prod_path.replace("_CM001", "")
-        prod_path_cloud_proba = L2aProduct(indice_product.l2a).msk_cldprb_20m
-        
+        l2a = L2aProduct(indice_product.l2a)
+        if l2a.path.exists():
+            prod_path_cloud_proba = l2a.msk_cldprb_20m
+        else:
+            prod_path_cloud_proba = None
         #~ prod_path = proc_item[2] / proc_item[0].identifier[:(-12 - len(proc_item[3]))] / proc_item[0].identifier
         #~ prod_path_unmasked = proc_item[2] / proc_item[0].identifier[:(-12 - len(proc_item[3]))] / (proc_item[0].identifier[:-11] + '.jp2')
         #~ prod_path_cloud_proba = L2aProduct(proc_item[0].identifier[:(-12 - len(proc_item[3]))]).msk_cldprb_20m