Project

General

Profile

Crocus Development in progress ยป

Download Sentinel2 data for SURFEX-Crocus experiments evaluation

1. Create a Theia account : https://sso.theia-land.fr/theia/register/register.xhtml

2. Get download script : https://github.com/olivierhagolle/theia_download

git clone https://github.com/olivierhagolle/theia_download.git
cd theia_downlad

3. Set configuration file : config_theia.cfg
serveur = https://theia.cnes.fr/atdistrib
resto = resto2
token_type = text
login_theia =
password_theia = password_theia

DO NOT FORGET TO SECURE config_theia.cfg :

chmod 600 config_theia.cfg

4. Set MF proxy :

export http_proxy=http://vernaym::11011

5. Extract tiles ( https://umap.openstreetmap.fr/fr/map/theias-sentinel-2-snow-tiles_156646#10/45.1699/6.1345)

p theia_download.py -t T31TGL -c Snow -a config_theia.cfg -d 2021-09-01 -f 2022-08-31 --snow_level L3B-SNOW
p theia_download.py -t T31TGK -c Snow -a config_theia.cfg -d 2021-09-01 -f 2022-08-31 --snow_level L3B-SNOW

Products documentation : https://www.theia-land.fr/product/neige/
Reference paper : https://essd.copernicus.org/articles/11/493/2019/
NB : Snow cover duration computed between 1 september and 31 august

6. Unzip extracted files :

unzip MULTISAT_20210901-000000-000_L3B-SNOW_T31TGL_D.zip
unzip MULTISAT_20210901-000000-000_L3B-SNOW_T31TGK_D.zip

7. Concatenation + projection + domain extraction : https://gdal.org/programs/gdalwarp.html

gdalwarp -t_srs EPSG:2154 -te 937948 6439100 973434 6464196 -tap -tr 250 250 -r average MULTISAT_20210901-000000-000_L3B-SNOW_T31TGK_C_V1-0/MULTISAT_20210901-000000-000_L3B-SNOW_T31TGK_C_V1-0_SCD_R2.tif MULTISAT_20210901-000000-000_L3B-SNOW_T31TGL_C_V1-0/MULTISAT_20210901-000000-000_L3B-SNOW_T31TGL_C_V1-0_SCD_R2.tif 20210901_L3B-SNOW_SCD_R2.tif

gdalwarp -t_srs EPSG:2154 -te 937948 6439100 973434 6464196 -tap -tr 250 250 -r average MULTISAT_20210901-000000-000_L3B-SNOW_T31TGK_C_V1-0/MULTISAT_20210901-000000-000_L3B-SNOW_T31TGK_C_V1-0_SMD_R2.tif MULTISAT_20210901-000000-000_L3B-SNOW_T31TGL_C_V1-0/MULTISAT_20210901-000000-000_L3B-SNOW_T31TGL_C_V1-0_SMD_R2.tif 20210901_L3B-SNOW_SMD_R2.tif

8. Convert to netcdf

gdal_translate -of NETCDF 20210901_L3B-SNOW_SCD_R2.tif 20210901_L3B-SNOW_SCD_R2.nc
gdal_translate -of NETCDF 20210901_L3B-SNOW_SMD_R2.tif 20210901_L3B-SNOW_SMD_R2.nc