Unit 20 - Sentinel downloader

There are plenty of libraries or tools which allow downloading Sentinel products from Copernicus Open Access Hub.

For GRASS GIS there is the i.sentinel toolbox. It consists of six GRASS addon modules:

Addons modules are not internal part of GRASS installation. Addons can be easily installed by g.extension (Settings ‣ Addons extensions ‣ Install extension from addons).

g.extension extension=i.sentinel
../_images/g-extension1.png

Fig. 111 Install i.sentinel GRASS Addon.

Note that i.sentinel.download requires also sentinelsat library to be installed.

Note

Sentinelsat library can be easily installed from the Console tab by entering the command below.

python3 -m pip install sentinelsat pandas

Switch to jena-region location.

Download data

Important

Pre-downloaded Sentinel scenes are available in the sample dataset jena-sample-data-sentinel.7z (1.7GB) geodata/sentinel. Readers can continue with importing sample data.

Let’s download suitable Sentinel products for our area of interest (AOI) and perform the NDVI calculation as described in Unit 05 - Raster processing (implemented as a model in Unit 08 - Modeler or as a Python script in Unit 11 - PyGRASS scripting). AOI region is defined by Jena city region created in Unit 03 - Data Management.

Sentinel-2 L2A products will be used to avoid computing atmospheric corrections. Let’s search for the latest available product by means of i.sentinel.download. Setting the -l flag, the result will only be printed. The download procedure will be performed later. In order to search and download Sentinel products from the Copernicus Open Access Hub, you have to create an account first. See the manual page of i.sentinel.download module for details. Create a new text file sentinel.txt containing two lines (username and password).

Note

To get username and password you need to register at the Copernicus Open Access Hub, see Register new account page for signing up.

i.sentinel.download -l map=jena_boundary producttype=S2MSI2A settings=sentinel.txt
11 Sentinel product(s) found
37438a05-f795-4e77-b68e-3a3c0ec5f04b ... 2022-01-06T10:23:19Z  4% S2MSI2A 1.08 GB
b4bc9d3b-1e1a-4cfd-bd01-14c21a596e8a ... 2022-02-05T10:21:19Z 64% S2MSI2A 1.03 GB
f2a604a6-8e50-4c62-b98b-22077365df7a ... 2022-01-21T10:23:31Z 65% S2MSI2A 933.01 MB
...

By default the module returns products for the last 60 days. Let’s change the search period setting start and end options. To be sure that our AOI is fully covered by a Sentinel product we also set area_relation option. We will also limit products by clouds coverage percentage threshold

i.sentinel.download -l map=jena_boundary producttype=S2MSI2A settings=sentinel.txt \
start=2021-04-01 end=2021-10-01 area_relation=Contains clouds=10
6 Sentinel product(s) found
a844500a-049f-46a3-92de-bcda2c38fc3c ... 2021-05-31T10:15:59Z  2% S2MSI2A 1.09 GB
d5b73db9-0acf-401d-9bf4-a6f199df1119 ... 2021-09-08T10:15:59Z  3% S2MSI2A 1.09 GB
b00d5dfd-9cce-48c6-a011-fd46b85de814 ... 2021-09-03T10:20:21Z  3% S2MSI2A 1.09 GB
...

Let’s download the desired product(s). Just remove the -l flag and add the output option in order to define the path to the output directory where data should be saved.

i.sentinel.download map=jena_boundary producttype=S2MSI2A settings=sentinel.txt \
start=2021-04-01 end=2021-10-01 area_relation=Contains clouds=10 \
limit=1 output=/home/user/geodata/sentinel/2021

Note

Note all products are available online. In this case the module fails with Product a844500a-049f-46a3-92de-bcda2c38fc3c is not online. Triggering retrieval from long term archive. error message. In this case try to download affected product in the next days.

Import data

Before importing or linking Sentinel data try to print a list of filtered raster files including projection match (second column, 1 for match otherwise 0). If the CRS of input data differs from the current location consider reprojection (-r) or creating a new location for import.

Data will be imported into the GRASS location by means of the i.sentinel.import tool. The command will import all Sentinel bands from input directory recursively. Before importing data let’s check content of the input directory by -p flag. The import procedure will be limited to the 4th and 8th bands in 10m spatial resolution by pattern option.

i.sentinel.import -p input=/home/user/geodata/sentinel/2019 pattern="20190626T102031_B0[4|8]_10m"
...L2A_T32UPB_A020940_20190626T102028/IMG_DATA/R10m/T32UPB_20190626T102031_B08_10m.jp2 1 (EPSG: 32632)
...L2A_T32UPB_A020940_20190626T102028/IMG_DATA/R10m/T32UPB_20190626T102031_B04_10m.jp2 1 (EPSG: 32632)

By default, input data are imported into GRASS data format. Alternatively, data can be linked if -l is given. It is also useful to import cloud mask vector features by -c flag. We also use register_output option to produce a timestamp plain text file which will be used in Unit 21 - Sentinel spatio-temporal.

i.sentinel.import -l -c input=/home/user/geodata/sentinel/2019 pattern="20190626T102031_B0[4|8]_10m"

Note

Cloud mask is computed by i.sentinel.import, the module doesn’t use cloud mask product stored in a SAFE directory.

Todo

Fix i.sentinel.import -c flag to respect pattern parameter.

Now launch the NDVI sample script created in Unit 12 - Script User Interface (ndvi-v4.py) in order to compute NDVI classes.

../_images/run-script1.png

Fig. 112 Run script to compute NDVI classes.

../_images/ndvi-classes.png

Fig. 113 Computed NDVI classes for given AOI.