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 1 and 2.5 degrees.
It is important that the user select images related to the same relative orbit.
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. As a suggestion, urbanized areas are usually well suitable to locate the reference point. 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. Acquisitions that do not cover the selected area are automatically discarded.

Output Format

Interferogram Generation Processing Mode:
a .tgz folder related to each interferogram pair (in the format _.tgz) is published.
Each folder contains:
- fr.dat: Interferogram as binary file in SAR coordinates;
- coh.dat: Coherence map as binary file in SAR coordinates;
- fr.dat_GOLD: Interferogram filtered with Goldstein kernel as binary file in SAR coordinates;
- fr.dat_GOLD_GEO_x: Geocoded Goldstein filtered interferogram as binary file;
- coh.dat_GEO_x: Geocoded Coherence map as binary file;
- fr.dat___GOLD_GEO.tif.kmz: Quick-look of the geocoded Interferogram importable in Google Earth.

Multitemporal Analysis Processing Mode:
a .tgz folder named “result_final” is published, which contains:
- One ASCII file containing all the outputs of the SBAS-InSAR chain processing.The first record contains, respectively: the number of field of each record; the number of acquisitions used in the time series and the measure unit for the mean displacement velocity.The second record contains the list of the epochs of acquisitions used in the time series, given as fraction of year.All the others records contain, per each pixel:
  - East and North UTM WGS84 coordinates (in meters);
  - Temporal coherence value;
  - Mean displacement velocity, calculated as linear fit of the relevant time series (in centimetres per year);
  - SAR coordinates (in pixels);
  - geographic WGS84 coordinates (latitude and longitude, in degrees);
  - residual topography w.r.t. the used DEM (in meters);
  - deformation time series (in centimetres).
    The basic scheme of the output ASCII file is reported in Figure 1.
    
    The file name is: SBAS_output4GIS_
    
    _.txt, where the two date indicate the starting and ending acquisitions used for generating the displacement time series.
- One .kmz file containing the Quick-look of the retrieved mean deformation velocity importable in Google Earth.

Conventions and Assumptions

Note that all the published binary files (.dat) are in 32-bit floating point raster format (if not differently specified) without any header information.
The origin of each binary matrix is located in the Lower Left corner.
Results provided in SAR coordinates have the Azimuth direction along the columns (X) and Range direction along the rows (Y).
Results provided in GEO coordinates (Lat-Lon WGS84) have the Longitude direction along the columns (X) and the Latitude direction along the rows (Y).

In the naming convention, the represents the day, month and year of a specific date while is a three character code referring to the used sensor (e.g, ERS, ENV, …).

SBAS Time Series

Figure 1 Format of the ASCII file, which encloses the outputs of the SBAS-InSAR processing chain

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