Homeworks for the course Earth Observation Data Analysis, 2020, Sapienza University of Rome
This repository contains three pdf files, one for each project undertaken during the course Earth Observation Data Analysis (Control Engineering, DIAG, Sapienza University of Rome).
All the homeworks have been done through the use of the Sentinel Application Platform (SNAP) toolbox, an open source software which is used for Earh Observation processing and analysis. You can find more information here and you can download it here.
The satellite image data are available from the following sources:
- MODIS data from Terra and Aqua satellites from NASA LAADS DAAC
- MSI Sentinel-2 data from Copernicus
- SAR Sentinel-1 data from Copernicus
Modis images from NASA LAADS DAAC. The objective of this homework was the quantitative detection of the vegetation coverage between two MODIS images representing the same region of interest (ROI) in the winter and summer seasons.
- Download of a MODIS image at 1-km resolution over a ROI during summer season.
- Quality check
- Data analysis (spectrum, histogram, profile tools)
- Channel data correlation of whole image and selected ROI
- Principal Component Analysis computation
- Unsupervised and supervised classification with sea, land and cloud classes
- Application of the formulas of 2-band and 3-band vegetation index (VI)
- Unsupervised classification using a VI index in place of a MODIS channel
- Application of VI indices to the second MODIS image in the winter season
- Quantitative change detection of the vegetation coverage class by reprojecting the two MODIS images in July and January.
Here is the vegetation coverage of Italy in July and January seasons:
Quantitative change detection:
## Homework 02: Surface mapping from MSI Sentinel-2 Data 🌿 💧Sentinel-2 MSI images downloaded from Copernicus.
The objective was the estimation of vegetation cover, inland water and chlorophyll-a sea concentration of Sentinel-2 MSI data and supervised classification within a ROI.
- Data exploration: download of two summer and winter MSI images choosing a coastal target area (Tiber river estuary in Italy).
Winter images of Ronciglione lake:
Summer images of Ronciglione lake:
-
Combination of the B2, B3, B4 channels to create an RGB image
-
L2 MSI products for scene classification masks: cloud, vegetation, soil, water.
-
SNAP Graph Builder tool applied to a subset around a lake or a river, resampling of the channels at 10 meters and their visualization.
-
Water and vegetation normalized index: NDWI and NDWI2 on winter product.
- Apply NDWI to highlight water content in the leaves and NDWI2 to build a water mask estimate chlorophyll-a (Chl-a) and total suspended sediments (TSS).
- Computation of NDVI on winter and summer MSI resampled products.
February and August water masks and difference water mask:
- Sea Chlorophyll-a and sediment estimation of an estuary/delta
- EmpReg regressive algorithm to Chl-a and TSS estimation
- Comparison with MCI (Maximum Chlorophyll Index) plugged-in algorithm
- Comparison with Chl-a and TSS products (from Copernicus Marine Service)
February and August Tiber river images:
Chlorophyll-a and Total Suspended Sediments differences between the two seasons:
- Scene classification
- Scene classification masks over an L2 MSI summer image defining vegetation, water, urban and bare soil as training areas.
- Reprojection and SNAP supervised classification algorithms (Maximum Likelihood and Random Forest).
- Comparison between the classifiers
Images downloaded from Copernicus and provided by ESA.
- Earthquake detection by SAR differential interferometry (Amatrice Earthquake in 2016, Italy) :world_map:
- DInSAR processing steps
Displacement map:
Displacement map terrain corrected and displacement map in Google Earth:
- Ship detection by SAR backscattering (Messina strait in 2020, Italy) 🚢
- SAR processing steps
Detected ship in the Messina strait:
- Flood detection by SAR backscattering (Flood in Mozambico, 2019) ⛈️
- DInSAR processing steps
Pre-event and post-event acquisitions:
Change detection after color manipulation: