Sentinel 2 biodiversity
Author(s) | Marie Josse Yvan Le Bras |
OverviewQuestions:Objectives:
How to get spectral and biodiversity indicators from remote sensing data ?
Which kind of ecosystem are you studying ? global or canopy data ?
Requirements:
Getting Sentinel 2 data and reformatting them in a generalized way
Computing spectral indices such as the NDVI
Calculating and visualizing biodiversity indicators
Comparing with in-situ data
Time estimation: 48 hoursSupporting Materials:Published: Apr 13, 2023Last modification: Jun 27, 2024License: Tutorial Content is licensed under Creative Commons Attribution 4.0 International License. The GTN Framework is licensed under MITpurl PURL: https://gxy.io/GTN:T00333version Revision: 5
This tutorial will guide you on getting Sentinel 2 data and processing them in order to calculate and visualize biodiversity indicators. This workflow made of 6 tools will allow you to explore Sentinel 2 data in the view of making biodiversity analyses.
Spatial diversity measurements should not replace in situ biodiversity data, but rather complement existing data and approaches. Spatial diversity estimates are currently based on long time scales, allowing more general predictions about rates of change in diversity. In practice, spatial data incorporate information on surface properties, including functional aspects, taxonomy, phylogeny and genetic diversity.
The tools explained here are useful for observing variations in spatial and temporal ecosystem properties given the intrinsic relationship between spatial variations in ecosystems and pixel values of spectral signals. A single measurement cannot provide a complete description of all the different aspects of ecosystem heterogeneity. Therefore, the combination of multiple tools in a Galaxy-Ecology workflow offers multiple approaches to unravel the complexity of ecosystem heterogeneity in space and time.
So, we will compute biodiversity and spectral indices mainly using reflectance information.
The reflectance is a proportion on reflected light on an area. It’s the ratio between the electromagnetic incident wave on the area and the reflected wave. It’s often a percentage between reflected intensity and incident intensity assumed as energy quantity.
Each part of this workflow has elementary steps :
- A first step to preprocess Sentinel 2 data:
- Preprocess Sentinel 2 data
- A second step to compute biodiversity indicators
- Global overview
- For Canopy
- A third step to compute spectral indices:
- Spectral indices
- EBV
Spectral indices are used to highlight particular features or properties of the earth’s surface, e.g. vegetation, soil, water. They are developed on the basis of the spectral properties of the object of interest.
Knowledge of the leaf cell, plant structure, state, condition and spectral properties is essential to perform vegetation analysis using remote sensing data. Spectral indices dedicated to vegetation analysis are developed on the basis that healthy vegetation reflects strongly in the near infrared (NIR) spectrum while absorbing strongly in the visible red.
In this tutorial, we’ll be working on Sentinel 2 data extracted from the Theia Land portal. First those data will be prepared. After pre-processing to fit the input format of the tools, we’ll see how to calculate biodiversity metrics.
AgendaIn this tutorial, we will cover:
Warning: It can take a bit of timeAs your processing heavy data some of the steps can take some time. Notably, the step Create Biodiversity maps can run for 2 days. We advise to let it run and proceed with the tutorial. Once, the pre-processing part is done all the steps can be done separately. Thus, no need for you to wait the end of each tool before running another one.
Upload and pre-processing of the data
This first step consist of downloading and properly prepare the data to use it in Sentinel 2 toolsuite.
Hands-on: Data upload
Create a new history for this tutorial and give it a name (example: “Sentinel 2 data for biodiversity tutorial”) for you to find it again later if needed.
To create a new history simply click the new-history icon at the top of the history panel:
Download the files from Scihub, PEPS or Theia :
You will have to to create an account for either of these platform. Select Reflectance, Sentinel 2 and “Niveau 2A” (level 2A).
This an example of the Theia land portal. You need to download a zip folder. Keep it that way.
Upload the zip folder
- Click on Upload Data on the top of the left panel
- Click on Choose local file and select the files or drop the files in the Drop files here part
- Click on Start
- Click on Close
You can rename galaxy-pencil the dataset,
sentinel_2_data.zip
for example, and to keep informations about the original nameSENTINEL2A.....zip
- Click on the galaxy-pencil pencil icon for the dataset to edit its attributes
- In the central panel, change the Name field
- Click the Save button
Preparing data
Using preprocS2 R package, this step provides a unique tool to read, crop, resample the original image directory, and write it as a raster stack
Preprocessing sentinel 2 data
Hands-on: Preprocess
- Preprocessing sentinel 2 data ( Galaxy version 0.0.1) with the following parameters:
- param-file “Input data”:
sentinel_2_data.zip
(Input dataset)- param-select “Where does your data come from ?”: ‘From Theia’
Click on Execute
CommentThe interesting output is the ENVI image format which is a binary raster file with an accompanying header file. The data are stored as a binary stream of bytes in a BIL file without extension and the metadata are stored in the .hdr file. These data are in the output Reflectance. You can directly use the output Reflectance for the rest of the analysis.
Question
- What are the files you are interested in for the following tools ?
- The 2 files in the Reflectance folder that finish by “_Refl” and “_Refl.hdr”
Producing biodiversity indicators
You can choose to compute spectral and biodiversity indicators either for global remote sensing data or for a canopy.
Hands-on: Choose Your Own TutorialThis is a "Choose Your Own Tutorial" section, where you can select between multiple paths. Click one of the buttons below to select how you want to follow the tutorial
Here you can choose which tutorial you want to folllow according to if your are more interested about studying canopy remote sensing data or more global ones
Compute a PCA
Hands-on: Principal components analysis for remote sensing data
- Compute a PCA ( Galaxy version 0.0.1) with the following parameters:
- param-select “In which format are your data ?”: ‘The data you are using are in a zip folder Reflectance’
- param-file “Input data”:
Reflectance
(output of Preprocessing sentinel 2 data tool)- param-select “Do you want to do a PCA or a SPCA ?”: ‘PCA’
Check that the “Input raster” datatype is bil and that “Input raster header” datatype is hdr
- Go on your raster data
- Click on galaxy-pencil to edit it
- Click on galaxy-chart-select-data Datatypes
- On “New Type” Select bil
Press Save
CommentDo the same for the raster header with the datatype hdr
Compute biodiversity indices
Hands-on: Biodiversity indicators for global remote sensing data
- Compute biodiversity indices ( Galaxy version 0.0.1) with the following parameters:
- param-select “In which format are your data ?”: ‘The data you are using are in a zip folder Reflectance’
- param-file “Input data”:
Reflectance
(output of Preprocessing sentinel 2 data tool)- param-text “Write a number of the value of alpha”: ‘1’
Check that the “Input raster” datatype is bil and that “Input raster header” datatype is hdr
- Go on your raster data
- Click on galaxy-pencil to edit it
- Click on galaxy-chart-select-data Datatypes
- On “New Type” Select bil
Press Save
CommentDo the same for the raster header with the datatype hdr
- param-select “In which format are your data ?”: ‘Your already have the files in ENVI BIL format’
- param-file “Input raster”:
PCA raster
(output of Compute a PCA tool)- param-file “Input header”:
PCA header
(output of Compute a PCA tool)- param-text “Write a number of the value of alpha”: ‘1’
Here again check that the “Input raster” datatype is bil and that “Input raster header” datatype is hdr
These 2 pictures are only exemple it is normal if you don’t have exactly the same output.
Question
- How many biodiversity indicators do you have ?
- You should have 7 of them (Shannon, Renyi, Prao, Pielou, Hill, CRE, Berger-Parker). If you have only 5 of them, no problem, it just means you data are too small to compute CRE and Pielou but you can still continue your analysis.
Compute a PCA
Hands-on: Principal components analysis for remote sensing data
- Compute a PCA ( Galaxy version 0.0.1) with the following parameters:
- param-select “In which format are your data ?”: ‘The data you are using are in a zip folder Reflectance’
- param-file “Input data”:
Reflectance
(output of Preprocessing sentinel 2 data tool)- param-select “Do you want to do a PCA or a SPCA ?”: ‘PCA’
Check that the “Input raster” datatype is bil and that “Input raster header” datatype is hdr
- Go on your raster data
- Click on galaxy-pencil to edit it
- Click on galaxy-chart-select-data Datatypes
- On “New Type” Select bil
Press Save
CommentDo the same for the raster header with the datatype hdr
Create Biodiversity maps
Hands-on: Biodiversity indicators for canopy remote sensing data
- Mapping diversity ( Galaxy version 0.0.1) with the following parameters:
- param-select “In which format are your data ?”: ‘The data you are using are in a zip folder Reflectance’
- param-file “Input data”:
Reflectance
(output of Preprocessing sentinel 2 data tool)- param-select “Alpha, beta, functional diversity and comparison plot and map”: ‘All of the above’
Check that the “Input raster” datatype is bil and that “Input raster header” datatype is hdr
- Go on your raster data
- Click on galaxy-pencil to edit it
- Click on galaxy-chart-select-data Datatypes
- On “New Type” Select bil
Press Save
CommentDo the same for the raster header with the datatype hdr
Processing remote sensing data
Hands-on: Comparing biodiversity indicators for canopy
- Comparing remote sensing data ( Galaxy version 0.0.1) with the following parameters:
- param-select “In which format are your data ?”: ‘The data you are using are in a zip folder Reflectance’
- param-file “Input data”:
Reflectance
(output of Preprocessing sentinel 2 data tool)- param-file “Plots folder zip”:
output
(Input dataset)Comment: ShapefilesHere you must provide your folder of shapefiles (at least 2 in order to have the beta diversity).
Comment: Check datatype if you use your own ENVI BIL filesSame as the compute spectral indices make sure you have the right datatypes bil and hdr.
Question
- What kind of data do you need to use these tools ?
- Do you need a shapefile for mapping the diversity ?
- Why do you need multiple locations for comparing biodiversity ?
- This analisys is for data on forest, it’s a canopy study.
- No, only for the comparison with in situ data.
- The Bray curtis table compare the diversity between locations.
Spectral indices
Hands-on: Compute spectral indices
- Compute spectral indices ( Galaxy version 0.0.1) with the following parameters:
- param-select “In which format are your data ?”: ‘The data you are using are in a zip folder Reflectance’
- param-file “Input data”:
Reflectance
(output of Preprocessing sentinel 2 data tool)param-select “Input the type of indice you want”: ‘NDVI’
- param-select “Do you want the raster layer of the indice as an output ?”: ‘No’
Check that the “Input raster” datatype is bil and that “Input raster header” datatype is hdr
- Go on your raster data
- Click on galaxy-pencil to edit it
- Click on galaxy-chart-select-data Datatypes
- On “New Type” Select bil
Press Save
CommentDo the same for the raster header with the datatype hdr
CommentYou can choose whichever indice you want
Remotely sensed diversity is consistent with most of the essential spatially constrained biodiversity variables proposed by Skidmore et al. (2015). This highlights the need for increased dialogue and collaboration between the biodiversity monitoring community and the remote sensing community to make satellite remote sensing a tool of choice for conservation. Increased dialogue is also essential within the biodiversity monitoring community to achieve this. From this point of view multiple Satellite Remote Sensing EBV (SRS EBV) were created. Some of the indices proposed here will allow you to compute SRS EBV. For instance it allows you to compute one of GEO BON EBV Canopy Chlorophyll Content. This EBV is computed by GEO BON on the Netherlands, here you can compute it on which ever Sentinel 2 data you want by chosing to calculate the indice CCCI.
Question
- What’s the difference between biodiversity indicators and spectral indices ?
- Biodiversity indicators give us informations on the heterogeneity of the landscape whereas spectral indices inform us on the well being of the vegetation.
Conclusion
You are now all set to use your remote sensing data in order to do a biodiversity analysis. Before you go on one last reflexion.
Question
- Should remote sensing replace in-situ data ?
- NO ! remote sensing and in-situ data should come and complete one another to have the most complete view of the state of biodiversity.