Small BAseline Subset – SBAS inSAR (Beta) for Sentinel-1

USER MANUAL

Small BAseline Subset – SBAS inSAR (Beta) for Sentinel-1

Owner: IREA CNR

This is a short user manual for the SBAS-InSAR chain developed on the EarthConsole® P-PRO ON DEMAND infrastructure to help the users to properly exploit the service.

This manual is addressed to users already familiar with InSAR processing, analysis and products.

Introduction

Differential Synthetic Aperture Radar Interferometry (InSAR) is a technique for generating displacement maps of single events by calculating the phase difference (interferogram) of two SAR images acquired over the same area and from temporally separated acquisition orbits. By exploiting a large number of differential interferograms, several Advanced InSAR techniques have been developed, allowing to follow the temporal evolution of the ground displacements. From now on we particularly address the Advanced InSAR technique referred to as Small BAseline Subset (SBAS) algorithm [1].

Following up the collaboration between CNR-IREA [2] and ESA-RSS [3] in the framework of the ESA’s Grid Processing, the Parallel-SBAS (P-SBAS) [4] algorithm, which takes benefit from High Performance Computing resources, has been integrated within the EarthConsole® P-PRO On Demand [5] environment to provide users with a service for the generation of Earth’s surface deformation time series and, more generally, interferometric products, in unsupervised way. The P-PRO ON DEMAND flexibility and an user-friendly interface allows users to perform a SBAS-InSAR processing in a reasonable time. The aim of this activity is to provide the Earth Observation community with an automatic tool for improving the surface deformation analyses and thus making advances in understanding several geophysical phenomena.

The main user actions are the following:

select the SAR data to be processed through a defined computing element;
set input parameters for SBAS-InSAR processing;
obtain SBAS-InSAR geocoded (Lat/Lon and/or UTM WGS84) results, such as mean deformation velocity map and deformation time series.

The provided service performs the full SBAS InSAR chain from L1 data focusing to displacement time series generation.

The service provides in output SBAS-InSAR geocoded (Lat/Lon and/or UTM WGS84) results, such as mean deformation velocity map and deformation time series.

Search Data Parameters

Name	Value	Description
Start/Stop date		The user can select the temporal range of interest
BBOX		The user can either select the Area of Interest drawing it on the map or specifying it manually into this box
Number of results returned per page		Enter an appropriate number of results returned per page (e.g. 100)
Relative orbit		It allows the user to filter the data by relative orbit
Acquisition	Ascending/Descending	It allows the user to filter the data by the acquisition type

Note that this step must be carried out with particular care, since a wrong data selection can result to an unfeasible processing.

Recommendations:

The bounding box must have an extension in latitude between 0.7 and 2.5 degrees.
It is important that the user select images related to the same relative orbit.
The user should avoid temporal gaps larger than one year between two consecutive images in the data set selection.
It is preferable to select the reference point in this tab by clicking on “Edit layers” on the left (1) and then move the pin point with the mouse on the desired position (2). The reference point should be on land and included within the selected area of interest. By default the coordinates are set in the middle of the selected area of interest. As a suggestion, urbanized areas are usually well suitable to locate the reference point.

Processing Parameters

Name	Value	Description
Name		Name of the processing task
Results Publish Server	Default: local	The user can specify where the results will be published, inserting the publish server address. Username and password shall be specified respectively in the fields Publish server username and Publish server password. By default the results are published in the EarthConsole publish server.
Lat point reference Lon point reference	The reference point should be on land and included within the selected area of interest. It should be located in a stable area or its deformation behaviour shall be known As a suggestion, urbanized areas are usually well suitable to locate the reference point. To help the PhU, the algorithm automatically refines the reference point and selects the one with best coherence conditions close to the one selected by the user. By default the coordinates are set in the middle of the selected area of interest.	Point used as reference for the displacement measurements generated by the SBAS-InSAR algorithm. It can be set either manually inserting the coordinates into the boxes or automatically moving the Pin point over the map. To move the pin point on the map select “Edit layers” on the left (1) and then move the pin point with the mouse on the desired position (2). The reference point coordinates should update automatically. To finish moving the pin, click on “Finish” on the map.
Cut data over selected AOI	True/false	If checked, the data are cut over selected AOI. In this case acquisitions that do not cover the selected area are automatically discarded.

Output Format

Conventions and assumptions

Results are provided in the satellite Line Of Sight (LOS). Positive values indicate that the target moves toward the satellite. Processing results are provided according to the EPOS-IP project specifications along with the corresponding metadata.

Published Results

The main outputs of the MTA mode are the:

Displacement Time series in LOS;
Mean LOS velocity;
Temporal coherence;
Pixel location (coordinates);
LOS unit vectors.
Information is organized in a CSV ASCII table according to the following figure:

Provided information consists, per each pixel considered reliable, in:

Unique pixel identifier (ID);
WGS84 Latitude in degree (Lat);
WGS84 Longitude in degree (Lon);
Topography above the ellipsoid (Topo);
Mean Velocity in cm/year, as linear regression of the displacement Time series (Vel);
Temporal coherence (Coer);
Components of LOS unit vector along the North, East and Vertical directions (cosN, cosE, cosU);
LOS displacement time-series in cm (TS): the length of this field depends on the number of acquisitions used in the time series generation.

Citation and Feedback

Commercial use of any part of this service is not allowed without express permission from the CNR-IREA Institute and EarthConsole®.

Users are encouraged to use the InSAR service here described for any scientific aim. Please recognize the effort made by the authors by citing:

Casu, S. Elefante, P. Imperatore, I. Zinno, M. Manunta, C. De Luca and R. Lanari, “SBAS-DInSAR Parallel Processing for Deformation Time-Series Computation,” IEEE JSTARS, vol. 7, no. 8, pp. 3285-3296, 2014, doi: 10.1109/JSTARS.2014.2322671

in relevant talks and publications prepared by using InSAR results generated by this service.

Users are also kindly invited to report any issue and problem encountered during the use of this service to support@earthconsole.eu and sbas-help@irea.cnr.it. Moreover, suggestions and comments are warmly appreciated to keep the service as much as possible appealing, effective and efficient.

Acronyms & References

ACRONYM	DESCRIPTION
APS	Atmosferic Phase Screen
CNR	National Research Council of Italy
DEM	Digital Elevation Model
EMCF	Extended Minimum Cost Flow
ENVISAT	Environment Satellite
ERS	European Remote Sensing Satellite
ESA	European Space Agency
GIS	Geographic Information System
InSAR	Synthetic Aperture Radar Interferometry
IREA	Institute for the Electromagnetic Sensing of Environment
MCF	Minimum Cost Flow
P-PRO ON DEMAND	Parallel Processing On Demand
P-SBAS	Parallel Small Baseline Subset
SBAS	Small BAseline Subset
SRTM	Shuttle Radar Topography Mission
SVD	Singular Value Decomposition
TEP	Thematic Exploitation Platform
UTM	Universal Transverse Mercator
VA4	Virtual Archive 4

References

[1] P. Berardino, G. Fornaro, R. Lanari, and E. Sansosti, “A new algorithm for surface deformation monitoring based on small baseline differential SAR interferograms,” IEEE Trans. Geosci. Remote Sens., vol. 40, no. 11, pp. 2375–2383, Nov. 2002
[2] http://www.irea.cnr.it/
[3] http://wiki.services.eoportal.org/rss-portal.php
[4] F. Casu, S. Elefante, P. Imperatore, I. Zinno, M. Manunta, C. De Luca and R. Lanari, “SBAS-DInSAR Parallel Processing for Deformation Time-Series Computation,” IEEE JSTARS, 2014, doi: 10.1109/JSTARS.2014.2322671
[5] https://ui-ppro.earthconsole.eu/
[6] R. M. Goldstein and C. L. Werner, “Radar interferogram filtering for geophysical applications,” Geophys. Res. Lett., vol. 25, no. 21,pp. 4035–4038, 1998.
[7] E. Sansosti, P. Berardino, M. Manunta, F. Serafino and G. Fornaro, “Geometrical SAR image registration,” IEEE Trans. Geosci. Remote Sens., vol. 44, no. 10, pp. 2861–2870, Oct. 2006
[8] A. Pepe and R. Lanari, “On the extension of the minimum cost flow algorithm for phase unwrapping of multitemporal differential SAR interferograms,” IEEE Trans. Geosci. Remote Sens., vol. 44, no. 9, pp. 2374–2383, Sep. 2006

Small BAseline Subset – SBAS inSAR (Beta) for Sentinel-1

Resources

Headquarters