From 35bcec8211733e45d31d06ef24c4e3bac0707db0 Mon Sep 17 00:00:00 2001 From: Raffaele Gaetano <raffaele.gaetano@cirad.fr> Date: Thu, 23 Nov 2023 16:10:03 +0100 Subject: [PATCH] FIX: bug on time series mosaicking. --- TimeSeries/s2theia.py | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/TimeSeries/s2theia.py b/TimeSeries/s2theia.py index cc5ade7..8c3f1fd 100644 --- a/TimeSeries/s2theia.py +++ b/TimeSeries/s2theia.py @@ -852,9 +852,11 @@ class S2TheiaPipeline: out.extend(t.write_outputs(out_fld)) t.reset() if len(self.tiles) > 1 and mosaicking == 'vrt': + n_dates = len(self.tiles[0].output_dates) + out = [out[i:i+n_dates] for i in range(0,len(out),n_dates)] out_mos = [] vrtopt = gdal.BuildVRTOptions() - for i in range(len(self.tiles[0].output_dates)): + for i in range(n_dates): fn = out_fld + os.sep + 'SENTINEL2_MOSAIC_GAPFILL_' + self.tiles[0].output_dates[i] + '_' + stack_name + '.vrt' to_mosaic = [x[i] for x in out] gdal.BuildVRT(fn, to_mosaic, options=vrtopt) -- GitLab