The aim of the project is to familiarize students with calculating Normalized satellite indexes, electromagnetic spectrum bands, visualization of map data. Besides the development of software, students should also work on a literature review of the topic and prepare a demonstration of Normalized Satelite Indexes (NSIs) for their final presentation.
The result of the project should be an interface that could demonstrate differences between NSIs. The interface could be graphical (GUI) or at least command line-based (CLI) capable of generating adequate visualizations.
Input data of the project should be data downloaded from any of the services that provide the Sentinel-2 dataset. Usage of other satellite missions is allowed if students present sufficient arguments for the change in the assignment. The toy dataset created in this project might contain a combination of data from several missions – for instance, Sentinel-2 primarily, and for situations when Sentinel-2 data are not available, MODIS.
Some steps of the work
The project should have roughly the following parts:
- Review of existing normalized satellite indexes/spectral indices
- Mapping the satellite indexes to Sentinel-2 bands or bands of the products used to construct the toy dataset
- Table of used normalized satellite indexes/spectral indexes with references to literature, notes, examples, etc.
- Approximately 5-page review.
- Construction of a toy dataset that includes samples for different times of the year at the same location, different geographical locations, known fire areas, etc. Students should develop software, which handles the downloading of the data, for instance from the Copernicus Open Access Hub.
- Creating or adapting a user interface capable of displaying the maps and raster overlays that supports quick “switching” and comparison of the various satellite indexes. The user interface should use adequate colormap for each of the used satellite indexes. Ideally, all data in the toy dataset should be shown – in a similar way to the Google Maps’ satellite view.
- Preparation of a presentation of the developed software, including demonstrations of the satellite normalized indexes on the data from the toy dataset constructed in this project. There should be a suitable example for each of the indices calculated in this project.
Computational resources
A remote computer for long-running code can be provided to students after discussion.
Semester 2020/2021
This project is assigned to Gabriel Gorta and Michal Marcinčin.
Links/Resources
- Detection of water bodies based on NDWI > 0
- RUS Webinar: Mapping waterbodies from space – HYDR01 (PDF, Tereza Šmejkalová, 3 April 2018)
- Patrick Gray’s example – chapter_5_classification
- Simple NDVI example by CSC Finland – Raster map algebra — Intro to Python GIS CSC documentation
- Spectral indexes
- Table of normalized satellite indexes: https://www.usna.edu/Users/oceano/pguth/md_help/html/norm_sat.htm
NDWI (Normalized Difference Water Index)
- Water value is positive
- McFeeters S.K. (1996)
- Monitor changes related to water content in water bodies
- NDWI = (G – NIR) / (G + NIR)
- In Sentinel-2:
- NDWI = (B03 – B08) / (B03 + B08)
- S. K. McFEETERS (1996) The use of the Normalized Difference Water Index (NDWI) in the delineation of open water features, International Journal of Remote Sensing, 17:7, 1425-1432, DOI: 10.1080/01431169608948714
- B.C. Gao. (1996)
- Monitor changes in water content of leaves
- (NIR – SWIR) / (NIR + SWIR)
- Wavelengths in the paper:
[ ⍴(860nm) – ⍴(1240nm) ] / [ ⍴(860nm) – ⍴(1240nm) ]- Band ⍴(1240nm) is a problematic
- Not present in to Sentinel-2, LANDSAT
- Applicable to MODIS data
- For Sentinel-2 some sources* use B11 and B12
- B.C. GAO. (1996) “NDWI—A normalized difference water index for remote sensing of vegetation liquid water from space.” Remote Sensing of Environment, 58:3, 257-266, DOI: 10.1016/S0034-4257(96)00067-3
- Xu H. (2006) – “Modified Normalized Difference Water Index”
- MNDWI = (Green − MIR)/(Green + MIR)
- MIR – Mid-Infrared, 300–5000 nm
- Sentinel-2: B11 and B12 are in MIR
- Hanqiu Xu (2006) Modification of normalised difference water index (NDWI) to enhance open water features in remotely sensed imagery, International Journal of Remote Sensing, 27:14, 3025-3033, DOI: 10.1080/01431160600589179
- MIR: LANDSAT TM band 5 ≈ Sentinel-2 band 11
Landsat Surface Reflectance-Derived Spectral Indices
- https://www.usgs.gov/core-science-systems/nli/landsat/landsat-surface-reflectance-derived-spectral-indices
- Normalized Difference Vegetation Index (NDVI)
- (NIR – R) / (NIR + R)
- In Landsat 8: NDVI = (Band 5 – Band 4) / (Band 5 + Band 4)
- In Sentinel-2:
- NDVI = (B08 – B04) / (B08 + B04)
- NDVI = (B8a – B04) / (B8a + B04)
- https://www.usgs.gov/core-science-systems/nli/landsat/landsat-normalized-difference-vegetation-index
- https://eos.com/ndvi/
- a measure of the state of plant health
- Based on reflectance of light at certain frequencies by a plant
- it helps to compensate for changes in lighting conditions, surface slope, exposure, and other external factors → preferable for global vegetation monitoring
- difference in the intensities of reflected light in the red and infrared range divided by the sum of these intensities
- values from -1.0 to 1.0
- negative values – mainly formed from clouds, water and snow
- values close to zero – primarily formed from rocks and bare soil
- Very small values (0.1 or less) – empty areas of rocks, sand or snow
- Moderate values (from 0.2 to 0.3) – shrubs and meadows
- large values (from 0.6 to 0.8) – temperate and tropical forests
- a measure of the state of plant health
- See: NDVI products provided by Copernicus LMS https://land.copernicus.eu/global/products/ndvi
- Enhanced Vegetation Index (EVI)
- EVI = G * ((NIR – R) / (NIR + C1 * R – C2 * B + L))
- https://www.usgs.gov/core-science-systems/nli/landsat/landsat-enhanced-vegetation-index
- Soil Adjusted Vegetation Index (SAVI)
- ((NIR – R) / (NIR + R + L)) * (1 + L)
- In Landsat 8: SAVI = ((Band 5 – Band 4) / (Band 5 + Band 4 + 0.5)) * (1.5)
- https://www.usgs.gov/core-science-systems/nli/landsat/landsat-soil-adjusted-vegetation-index
- Modified Soil Adjusted Vegetation Index (MSAVI)
- MSAVI = (2 * NIR + 1 – sqrt ((2 * NIR + 1)2 – 8 * (NIR – R))) / 2
- In Landsat 8: MSAVI = (2 * Band 5 + 1 – sqrt ((2 * Band 5 + 1)2 – 8 * (Band 5 – Band 4))) / 2
- https://www.usgs.gov/core-science-systems/nli/landsat/landsat-modified-soil-adjusted-vegetation-index
- Normalized Difference Moisture Index (NDMI)
- (NIR – SWIR) / (NIR + SWIR)
- In Landsat 8: NDMI = (Band 5 – Band 6) / (Band 5 + Band 6)
- Sentinel-2:
- NDMI = (B8A – B11) / (B8A + B11)
- NDMI = (B08 – B11) / (B08 + B11)
- https://www.usgs.gov/core-science-systems/nli/landsat/normalized-difference-moisture-index
- Normalized Burn Ratio (NBR)
- (NIR – SWIR) / (NIR + SWIR)
- In Landsat 8: NBR = (Band 5 – Band 7) / (Band 5 + Band 7)
- Sentinel-2:
- NBR = (B8A – B12) / (B8A + B12)
- NBR = (B08 – B12) / (B08 + B12)
- https://www.usgs.gov/core-science-systems/nli/landsat/landsat-normalized-burn-ratio
- Normalized Burn Ratio 2 (NBR2)
- (SWIR1 – SWIR2) / (SWIR1 + SWIR2)
- In Landsat 8: NBR2 = (Band 6 – Band 7) / (Band 6 + Band 7)
- Sentinel-2: NBR2 = (B11 – B12) / (B11 + B12)
- https://www.usgs.gov/core-science-systems/nli/landsat/landsat-normalized-burn-ratio-2
- Normalized Difference Vegetation Index (NDVI)
6 Spectral Indexes To Make Vegetation Analysis Complete
https://eos.com/blog/6-spectral-indexes-on-top-of-ndvi-to-make-your-vegetation-analysis-complete
- The Soil Adjusted Vegetation Index (SAVI)
- soil adjustment factor L to the equation of NDVI in order to correct for soil noise effects
- https://doi.org/10.1016/0034-4257(88)90106-X
- Atmospherically Resistant Vegetation Index (ARVI)
- notrelatively prone to atmospheric factors (such as aerosol)
- NDVI corrected for atmospheric scattering effects in the red reflectance spectrum by using the measurements in blue wavelengths
- more robust to topographic effects
- A highly effective monitoring tool for tropical mountainous regions
- ARVI = (NIR – (2 * Red) + Blue) / (NIR + (2 * Red) + Blue)
- Enhanced Vegetation Index (EVI)
- EVI = 2.5 * ((NIR – Red) / ((NIR) + (C1 * Red) – (C2 * Blue) + L))
- Aim to correct NDVI results for atmospheric influences and soil background signals
- coefficients
- C1 and C2 to correct for aerosol scattering present in the atmosphere,
- L to adjust for soil and canopy background
- NASA’s MODIS: C1=6, C2=7.5, and L=1
- Green Chlorophyll Index (GCI)
- GCI = (NIR) / (Green) – 1
- content of leaf chlorophyll in various species of plants
- reflects the physiological state of vegetation
- decreases in stressed plants
- Structure Insensitive Pigment Index (SIPI)
- SIPI = (NIR – Blue) / (NIR – Red)
- analysis of vegetation with the variable canopy structure
- Normalized Burn Ratio (NBR)
- Highlight burned areas following the fire
- Healthy vegetation shows high reflectance in NIR spectrum
- recently burned areas of vegetation reflect highly in the SWIR spectrum
- A common practice to assess burn extent and severity with the relative differenced NBR
References
- Table of normalized satellite indexes: https://www.usna.edu/Users/oceano/pguth/md_help/html/norm_sat.htm
- NDWI: by European Drought Observatory https://edo.jrc.ec.europa.eu/documents/factsheets/factsheet_ndwi.pdf
- NDVI: by NASA earth observatory
https://earthobservatory.nasa.gov/features/MeasuringVegetation/measuring_vegetation_2.php - Comparison of band numbers: