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

mit den Daten arbeiten

Das ist ein Datensatz im SAFE-Format. Wo sind die Daten?

GDAL

Was ist GDAL?

GDAL is a translator library for raster and vector geospatial data formats

GDAL bietet viele verschiedene allgemeine Werkzeuge und Treiber für Sentinel-2.


    ############################################################################
    # Prepare data
    ############################################################################

    # prepare working directory
    WORK_DIR=/home/user/Desktop/open_data_mit_open_source/
    mkdir -p ${WORK_DIR}
    cd ${WORK_DIR}

    # extract data
    unzip S2A_MSIL1C_20170619T103021_N0205_R108_T32ULB_20170619T103021.zip
    unzip S2A_MSIL2A_20170619T103021_N0205_R108_T32ULB_20170619T103021.zip

    # set Level 1C
    L1C=S2A_MSIL1C_20170619T103021_N0205_R108_T32ULB_20170619T103021.SAFE
    echo $L1C

    # set Level 2A
    L2A=S2A_MSIL2A_20170619T103021_N0205_R108_T32ULB_20170619T103021.SAFE
    echo $L2A

    # set single band
    L2A_B04=${L2A}/GRANULE/L2A_T32ULB_A010401_20170619T103021/IMG_DATA/R10m/L2A_T32ULB_20170619T103021_B04_10m.jp2


    ############################################################################
    # Retrieve information with GDAL
    ############################################################################

    # set gdal cache higher to calculate faster
    export GDAL_CACHEMAX=2000

    gdalinfo --version

    # retrieve information about full Level 1C Product
    gdalinfo ${L1C}/MTD_MSIL1C.xml

    # retrieve information about 10m Level 1C Products
    cd ${L1C}
    gdalinfo "SENTINEL2_L1C:MTD_MSIL1C.xml:10m:EPSG_32632"
    cd ..

    # retrieve information about full Level 2A Product
    gdalinfo ${L2A}/MTD_MSIL2A.xml

    # retrieve information about single Level 2A Band
    gdalinfo ${L2A_B04}


    ############################################################################
    # Manipulate data with GDAL
    ############################################################################

    # Create Virtual Raster. -projwin ulx uly lrx lry
    gdal_translate -of VRT -projwin 350000 5640000 378000 5615000 ${L2A_B04} L2A_B04_subset.vrt

Impressum - Datenschutzerklärung