[GRASS-user] does i.sentinel.mask only work with output from i.sentinel.preproc?

Hi,

I’m trying to use i.sentinel.mask on a Sentinel2 L2A.

i.sentinel.mask -r input_file=input_cloud_mask.txt cloud_mask=cloud shadow_mask=shadow cloud_threshold=25000 shadow_threshold=5000 mtd_file=MTD_TL.xml

The input_cloud_mask.txt file:

blue=L2A_T34UGD_20161206T094402_B02_10m
green=L2A_T34UGD_20161206T094402_B03_10m
red=L2A_T34UGD_20161206T094402_B04_10m
swir11=L2A_T34UGD_20161206T094402_B11_20m
nir=L2A_T34UGD_20161206T094402_B08_10m
swir12=L2A_T34UGD_20161206T094402_B12_20m
nir8a=L2A_T34UGD_20161206T094402_B8A_20m

I get the following error:

WARNING: Any rescale factor has been applied
ERROR: Raster maps must be DCELL o FCELL

I didn’t use i.sentinel.preproc because our colleagues already preprocess L1C images and provide L2A, but AFAIU here I should start from a L1C image, and preprocess it with i.sentinel.preproc otherwise I can’t apply i.sentinel.mask, is this correct?

Thanks for any pointers
Regards,

···

Margherita Di Leo

Hi Madi

El jue., 21 feb. 2019 13:20, Margherita Di Leo <diregola@gmail.com> escribió:

Hi,

I’m trying to use i.sentinel.mask on a Sentinel2 L2A.

i.sentinel.mask -r input_file=input_cloud_mask.txt cloud_mask=cloud shadow_mask=shadow cloud_threshold=25000 shadow_threshold=5000 mtd_file=MTD_TL.xml

The input_cloud_mask.txt file:

blue=L2A_T34UGD_20161206T094402_B02_10m
green=L2A_T34UGD_20161206T094402_B03_10m
red=L2A_T34UGD_20161206T094402_B04_10m
swir11=L2A_T34UGD_20161206T094402_B11_20m
nir=L2A_T34UGD_20161206T094402_B08_10m
swir12=L2A_T34UGD_20161206T094402_B12_20m
nir8a=L2A_T34UGD_20161206T094402_B8A_20m

I get the following error:

WARNING: Any rescale factor has been applied
ERROR: Raster maps must be DCELL o FCELL

I didn’t use i.sentinel.preproc because our colleagues already preprocess L1C images and provide L2A, but AFAIU here I should start from a L1C image, and preprocess it with i.sentinel.preproc otherwise I can’t apply i.sentinel.mask, is this correct?

I haven’t tried myself, but the error only says you need to provide double or float map. How do L2 products come, maybe you just need to rescale if they come as integer and that’s it. Have you tried?

Just guessing

Vero

Ciao Roberta,

···

Margherita Di Leo