Applied-GeoSolutions / gips

Geospatial Image Processing System
GNU General Public License v3.0
17 stars 5 forks source link

S2 gips_project mosaicking error #427

Open lledoux opened 6 years ago

lledoux commented 6 years ago

So there are products for these dates (as can be seen in each tile's respective directory), but for some reason gips_project doesn't want to mosaic them . I am running this in my virtual environment /scratch/lledoux./cmsidno/test-venv/ on rio:

(test-venv) lmelendy@rio:/projects/ForestDisturbance/data/Cambodia$ gips_project sentinel2 -s /scratch/lledoux/forest_disturbance/2017_cambodia_mangroves/cambodia_mangroves_landsat_tiles_dis.shp -d 2015-12-01,2016-03-31 -p evi ndvi lswi satvi ref cfmask --res 30 30 --outdir . -v4

The errors I receive are:

Creating mosaic project ./cambodia_mangroves_landsat_tiles_dis_30.0x30.0_sentinel2/0
  Dates: 2 dates (2015-12-01 - 2015-12-11)
  Products: cfmask evi lswi ndvi ref satvi
Error mosaicking ./cambodia_mangroves_landsat_tiles_dis_30.0x30.0_sentinel2/0/2015335_S2A_cfmask.tif:
Traceback (most recent call last):
  File "/home/lmelendy/code/gips/gips/utils.py", line 474, in cli_error_handler
    yield
  File "/home/lmelendy/code/gips/gips/tiles.py", line 94, in mosaic
    filenames = [self.tiles[t].filenames[(sensor, product)] for t in self.tiles]
KeyError: ('S2A', 'cfmask')
Error mosaicking ./cambodia_mangroves_landsat_tiles_dis_30.0x30.0_sentinel2/0/2015335_S2A_evi.tif:
Traceback (most recent call last):
  File "/home/lmelendy/code/gips/gips/utils.py", line 474, in cli_error_handler
    yield
  File "/home/lmelendy/code/gips/gips/tiles.py", line 94, in mosaic
    filenames = [self.tiles[t].filenames[(sensor, product)] for t in self.tiles]
KeyError: ('S2A', 'evi')

And these are the inputs of a directory where the products are properly mosaicked:

lmelendy@lmelendydt:/titan/data/sentinel2/tiles/48PTT/2017004$ ls -lrth
total 1.2G
-rw-rw-r-- 1 lmelendy staff 583M Oct 13 16:46 S2A_MSIL1C_20170104T033122_N0204_R018_T48PTT_20170104T033834.zip
-rw-rw-r-- 1 lmelendy staff  14K Oct 16 07:53 S2A_MSIL1C_20170104T033122_N0204_R018_T48PTT_20170104T033834.zip.index
-rw-rw-r-- 1 lmelendy staff 345M Oct 16 08:03 48PTT_2017004_S2A_ref.tif
-rw-rw-r-- 1 lmelendy staff  29M Oct 16 08:03 48PTT_2017004_S2A_cfmask.tif
-rw-rw-r-- 1 lmelendy staff  58M Oct 16 08:03 48PTT_2017004_S2A_satvi.tif
-rw-rw-r-- 1 lmelendy staff  58M Oct 16 08:03 48PTT_2017004_S2A_ndvi.tif
-rw-rw-r-- 1 lmelendy staff  58M Oct 16 08:03 48PTT_2017004_S2A_lswi.tif
-rw-rw-r-- 1 lmelendy staff  58M Oct 16 08:03 48PTT_2017004_S2A_evi.tif

And where the products are not able to be properly mosaicked:

lmelendy@lmelendydt:/titan/data/sentinel2/tiles/48PTT/2015335$ ls -lrth
total 6.9G
-rw-rw-r-- 14 lmelendy staff 6.3G Oct 13 15:59 48PTT_S2A_OPER_PRD_MSIL1C_PDMC_20170517T123855_R018_V20151201T033112_20151201T033112.zip
-rw-rw-r--  1 lmelendy staff 272K Oct 16 06:45 48PTT_S2A_OPER_PRD_MSIL1C_PDMC_20170517T123855_R018_V20151201T033112_20151201T033112.zip.index
-rw-rw-r--  1 lmelendy staff 345M Oct 16 06:53 48PTT_2015335_S2A_ref.tif
-rw-rw-r--  1 lmelendy staff  29M Oct 16 06:53 48PTT_2015335_S2A_cfmask.tif
-rw-rw-r--  1 lmelendy staff  58M Oct 16 06:53 48PTT_2015335_S2A_satvi.tif
-rw-rw-r--  1 lmelendy staff  58M Oct 16 06:53 48PTT_2015335_S2A_ndvi.tif
-rw-rw-r--  1 lmelendy staff  58M Oct 16 06:53 48PTT_2015335_S2A_lswi.tif
-rw-rw-r--  1 lmelendy staff  58M Oct 16 06:53 48PTT_2015335_S2A_evi.tif