iCOR sentinel2 AOT issue


#1

Hi,
I am using the icor.py script from command line on a Centos7 machine with python 2.7.5 and I installed the “ICOR UPDATE FOR SNAP6.0” patch.
When I run the script with landsat8 products everything works fine, but with sentinel2 products I am encountering some issues during the Aerosol Optical Thickness phase. The process finishes successfully only if I set the aot parameter to false.
I tried using both sentinel2’s products released before and after the 2017 naming convention change. Tried both S2A and S2B.
I tried to shorten the s2 product’s path as suggested in iCOR Landsat 8 error

This is the command I am running

/usr/bin/python /opt/vito/icor/src/icor.py \
--keep_intermediate false \
 --cloud_average_threshold 0.19  \
 --cloud_low_band B01  \
 --cloud_low_threshold 0.25  \
 --cirrus true  \
 --aot true  \
 --aot_window_size 100  \
 --simec false  \
 --watervapor false  \
 --bg_window 1 \
 --cirrus_threshold 0.01  \
 --aot_override 0.1  \
 --ozone_override 0.33  \
 --wv_override 2.0  \
 --water_band B08  \
 --water_threshold 0.05  \
 --data_type S2  \
 --output_file /tmp/icorTest/output/S2A_*_processed.tif  \ 
/tmp/data/S2A_MSIL1C_20170101T121632_N0204_R123_T24FWG_20170101T121634.SAFE/GRANULE/L1C_T24FWG_A007985_20170101T121634/IMG_DATA/T24FWG_20170101T121632_B08.jp2

and this is the error I am getting.

======================================================================
 62% Estimate the Aerosol Optical Thickness (AOT)
======================================================================
Running ac_runner /tmp/atcconfidXCFh.conf
Running ac_runner /tmp/atcconftNMb1e.conf


0%   10   20   30   40   50   60   70   80   90   100%
|----|----|----|----|----|----|----|----|----|----|
***************************************************

Failed!
Process output:
 +---------------------------------------------------------------------
 | :: load config file /tmp/atcconftNMb1e.conf
 | :: open image : /tmp/icor_sV3oZU_proc/S2A_OPER_MSI_L1C_TL_MTI__20170103T052120_A007985_T24FWG_N02_04/S2A_OPER_MSI_L1C_TL_MTI__20170103T052120_A007985_T24FWG_N02_04_AOT_60M.tif
 | :: image does not contain valid values
 +---------------------------------------------------------------------

Thank you in advance

Blasco


#2

Hello,

what is the surface content of the images ? Only water, clouds … ? It appears that the automatic AOT retrieval fails due to the absence of land (vegetation and soil) inside the image. Also, when the images are too cloudy, the process might fail.

greetings,
Stefan


#3

Thank you for your prompt reply.
The issue was indeed the image content which only had clouds.

Thank you very much!
greetings

Blasco