diff --git a/sen2chain/time_series.py b/sen2chain/time_series.py index b8c6d91b684a62130963f0e1c476b49e4d173134..389b276ce348aebb0ce92add45f8a6d49fc48ad9 100644 --- a/sen2chain/time_series.py +++ b/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