/dcs-rss-snap-s1-insar

Primary LanguageShellGNU General Public License v3.0GPL-3.0

Interferometric SAR processing of Sentinel 1 TOPSAR IW images with SNAP

Overview

This repository contains the application files and scripts to process a pair (Master and Slave) of Sentinel 1 TOPSAR IW data with SNAP (SeNtinel’s Application Platform), which is the common architecture for all Sentinel Toolboxes jointly developed by Brockmann Consult, Array Systems Computing and C-S. The interferometric processing chain for Sentinel 1 TOPSAR IW images is implemented through the tools contained in the Sentinel-1 Toolbox.

The interferometric SAR technique

Interferometric synthetic aperture radar (InSAR) exploits the phase difference between two complex radar SAR observations taken from slightly different sensor positions and extracts information about the earth’s surface. A SAR signal contains amplitude and phase information. The amplitude is the strength of the radar response and the phase is the fraction of one complete sine wave cycle (a single SAR wavelength). The phase of the SAR image is determined primarily by the distance between the satellite antenna and the ground targets. By combining the phase of these two images after coregistration, an interferogram can be generated whose phase is highly correlated to the terrain topography. The InSAR technique can potentially measure millimetre-scale changes in deformation over spans of days to years. It has applications for geophysical monitoring of natural hazards, for example earthquakes, volcanoes and landslides, and in structural engineering, in particular monitoring of subsidence and structural stability.

Sentinel-1 Interferometric Wide Swath Products

The Interferometric Wide (IW) swath mode is the main acquisition mode over land for Sentinel-1. It acquires data with a 250 km swath at 5 m by 20 m spatial resolution (single look). IW mode captures three sub-swaths using Terrain Observation with Progressive Scans SAR (TOPSAR). With the TOPSAR technique, in addition to steering the beam in range as in ScanSAR, the beam is also electronically steered from backward to forward in the azimuth direction for each burst, avoiding scalloping and resulting in homogeneous image quality throughout the swath. TOPSAR mode replaces the conventional ScanSAR mode, achieving the same coverage and resolution as ScanSAR, but with a nearly uniform SNR (Signal-to-Noise Ratio) and DTAR (Distributed Target Ambiguity Ratio). IW SLC products contain one image per sub-swath and one per polarisation channel, for a total of three (single polarisation) or six (dual polarisation) images in an IW product. Each sub-swath image consists of a series of bursts, where each burst has been processed as a separate SLC image. The individually focused complex burst images are included, in azimuth time order, into a single sub-swath image with black-fill demarcation in between, similar to ENVISAT ASAR Wide ScanSAR SLC products.

Quick link

Getting Started

To run this application you will need a Developer Cloud Sandbox, that can be either requested from:

A Developer Cloud Sandbox provides Earth Sciences data access services, and helper tools for a user to implement, test and validate a scalable data processing application. It offers a dedicated virtual machine and a Cloud Computing environment. The virtual machine runs in two different lifecycle modes: Sandbox mode and Cluster mode. Used in Sandbox mode (single virtual machine), it supports cluster simulation and user assistance functions in building the distributed application. Used in Cluster mode (a set of master and slave nodes), it supports the deployment and execution of the application with the power of distributed computing for data processing over large datasets (leveraging the Hadoop Streaming MapReduce technology).

Installation

Dependencies

SNAP

SNAP is the common architecture for all Sentinel Toolboxes. It should be automatically installed together with this application installation. SNAP software package is available through the Terradue Cloud Platform software repository. It can be manually installed through this command:

sudo yum install snap -y

SNAP is licensed under GNU GPL v3. For more information about SNAP go to SNAP.

SNAPHU

SNAPHU is the phase unwrapper processor developed by the Stanford Radar Interferometry Research Group. It should be automatically installed together with this application installation. SNAPHU software package is available through the Terradue Cloud Platform software repository. It can be manually installed through this command:

sudo yum install snaphu -y

Copyright information available at SNAPHU distribution README file.

Using the releases

Log on the developer cloud sandbox. Download the RPM package from https://github.com/geohazards-tep/dcs-rss-snap-s1-insar/releases. Install the downloaded package by running these commands in a shell:

sudo yum -y install rss-snap-s1-insar-<version>.noarch.rpm

Using the development version

Log on the developer sandbox and run these commands in a shell:

cd
git clone https://github.com/geohazards-tep/dcs-rss-snap-s1-insar.git
cd rss-snap-s1-insar
mvn install

Submitting the workflow

Invoke the Web Processing Service via the Sandbox dashboard or the Geohazards Thematic Exploitation platform providing the following items:

Master and Slave products' reference:

These are the URLs of the master and slave products to be processed, for example to have an interferogram related to the center Italy earthquake of August 2016:

Product polarisation

Define the polarization to be processed. Note that such polarization must be contained in the master-slave couple. For example VV is not present in "DH" or "SH" products. The following values can be chosen:

  • VV (default)
  • VH
  • HH
  • HV

Orbit type

Define the orbit source for the Orbit Correction:

  • Sentinel Precise
  • Sentinel Restituted (default)

NOTE: Precise orbit is available after 20 days of the data acquisition time.

DEM type

Define the DEM source for the Back-Geocoding Coregistration, Topographic Phase Removal and Terrain Correction processing:

  • SRTM 3Sec (default)

NOTE: SRTM valid in the [-56 deg,+60 deg] range of latitudes.

Azimuth coherence window size

Define the coherence estimation azimuth window size for the Interferogram processing [integer number of pixels]: 6 is used as default.

Range coherence window size

Define the coherence estimation range window size for the Interferogram processing [integer number of pixels]: 20 is used as default.

Azimuth Multilook factor

Define the multilook factor applied for Azimuth direction in the Multilooking processing [integer]: 1 is used as default.

Range Multilook factor

Define the multilook factor applied for Range direction in the Multilooking processing [integer]: 4 is used as default.

Perform phase unwrapping

Select if the phase unwrapping (through SNAPHU) has to be performed:

  • true
  • false (default)

Subset Bounding Box for Unwrapping

Define the area of interest for the Phase Unwrapping processing [lon_min,lat_min,lon_max,lat_max].

NOTE (1): This parameter is mandatory even if the Phase Unwrapping is not requested (i.e. when "Perform phase unwrapping"=false). Please put "-180.0,-56.0,180.0,60.0" as bounding box value in case of wrapped interferogram generation only.

NOTE (2): Area is limited to 0.25x0.25 degrees around the center of the selected area of interest.

Pixel spacing in meters

Define the pixel spacing for the Terrain-Correction processing [meters]: 15.0 is used as default.

To learn more and find information go to

Authors (alphabetically)

  • RSS Team

Questions, bugs, and suggestions

Please file any bugs or questions as issues or send in a pull request.

License

This application is licensed under the GNU GPL v3.