Open data mit Open source - Prozessierung von Copernicus Daten mit freier (Kommandozeilen-) Software

Atmosphärenkorrektur

Aus Zeitgründen werden wir hier keine Atmosphärenkorrektur durchführen, sondern von der ESA bereitgestellte L2A Daten verwenden.

Auch hier gibt es verschiedene Tools:

• sen2cor ist ein Plugin für SNAP. Es gibt es aber auch als isolierters Kommandozeilentool
• i.atcorr - i.atcorr ist ein Modul in GRASS GIS
• arcsi - arcsi ist als Dockerimage verfügbar:
  • https://hub.docker.com/r/mundialis/arcsi
  • docker pull mundialis/arcsi
  • docker run -t mundialis/arcsi arcsi.py -h

  • 
                # define name of Sentinel-2 scene - note: omit: .SAFE
                S2IMG=S2A_MSIL1C_20170327T105021_N0204_R051_T31UFT_20170327T105021
                DEM=srtm_30m_myregion.tif
                OUTDIR=arcsi_output_AOT_inv
                TMPDIR=~/tmp/arcsi
    
                # Note:
                #   remove RAD entry below to not keep this tmp dataset
    
                cd ${S2IMG}.SAFE/
                mkdir ${TMPDIR}
    
                # simple DOS1 correction example
                arcsi.py --sensor sen2 -i MTD_MSIL1C.xml -o ${OUTDIR} \
                --tmpath ${TMPDIR} -f KEA --stats -p RAD DOSAOTSGL SREF \
                --aeroimg ${CONDA_PREFIX}/share/arcsi/WorldAerosolParams.kea \
                --atmosimg ${CONDA_PREFIX}/share/arcsi/WorldAtmosphereParams.kea \
                --dem ${DEM} --minaot 0.05 --maxaot 0.6 --simpledos