Skip to content

DInSAR tutorial


image

This InSAR service derives line of sight displacement, interferometric phase and coherence from a pair of SAR complex images. It supports Sentinel-1 SLC data in TOPSAR (IW and EW) mode and SAOCOM SLC data in Stripmap (SM) mode only

📕 DInSAR service description and specifications are available in this section.


Use case 1: DInSAR with Sentinel-1 data

Select the processing service

After the opening of the activation workspace, in the right panel of the interface, open the Processing Services tab and select the processing service DInSAR Displacement Mapping - Sentinel-1 (DInSAR-S1).

Image-001

The "DInSAR Displacement Mapping - Sentinel-1 (DInSAR-S1)" panel is displayed with parameters values to be filled-in.

Image-002


Find the data using multiple filter criteria

Choose an area in which you want to focus your analysis. As an example the Menyuan County, Qinghai Province near the border with Gansu Province, China.

From the Navigation and Search toolbar (located in the upper left side of the map), click on “Spatial Filter” and draw a square around your Area of Interest. This spatial filter allows you to select only the EO data acquired over this area.

From the top of the left panel, use Filter Criterias to search for Radar and sentinel-1 data collections from the list.

To further filter results by selecting only SLC products you can use the Text filter by inserting SLC in the Search Terms text box always visible on the top left of the map.

After the query, the Results list will include only the features satisfying all the above defined filters.

Now it is possible to choose a pair of "Reference" and "Secondary" SLC Datasets to be used for DInSAR. This pair must come from the same sensor and shall have the same radar geometry (incident angle, orbit path).

Fill the parameters

After the definition of all the filters, you can employ DInSAR by importing at least a pair of Sentinel-1 SLC Datasets.

Job name

Insert as job name:

DInSAR - Sentinel-1 Earthquake in Jinchang China Jan 2022

Input Reference and Secondary datasets

Drag and Drop in the Reference SLC dataset and in the Secondary SLC dataset fields the following datasets:

  1. [CD] Sentinel-1A SLC VV/VH 33 29/12/2020 23:19:26,

  2. [CD] Sentinel-1A SLC VV/VH 33 10/01/2021 23:19:26.

Image-003

!! Tip For convention the reference is the oldest SAR acquisition and you can check how the motion was evolved by inserting the last acquired image as secondary.

Note

This service supports only SAR complex imagery.

Warning

DInSAR currently supports only SLC Datasets from the Sentinel-1 EO mission.

Warning

Drag and drop the Dataset (e.g. "[CD] Sentinel-1A SLC VV/VH 33 29/12/2020 23:19:26") and not the single-band asset (e.g. "magnitude-vh-iw1-001") into the Reference/Secondary SLC datasets fields.

Area of interest expressed as Well-known text

The “Area of interest as Well Known Text” can be defined by using the drawn polygon defined with the area filter.

Tip

In the definition of “Area of interest as Well Known Text” it is possible to apply as AOI the drawn polygon defined with the area filter. To do so, click on the button in the left side of the "Area of interest expressed as Well-known text" box and select the option AOI from the list. The platform will automatically fill the parameter value with the rectangular bounding box which is taken from the current search area in WKT format.

Note

This parameter is optional.

Once all parameters are inserted the processing service panel shall appear as the one shown in the figure below.

Image-004

Run the job

Click on the button Run Job and see the Running Job. You can monitor job progress through the progress bar.

Results: download and visualization

Once the job is completed, click on the button Show results at the bottom of the processing service panel.

Image-005

Tip

You can also save the parameters employed in this job by clicking on the Export params button in the right panel. This allows you to copy all your entries to the clipboard. This is meant to be used for a quick re-submission of a similar job after a fine tuning of the parameters (e.g. to add a colour formula later).

Below is reported the syntax which includes all the parameters employed in this example.

{
   "reference": "https://catalog.disasterscharter.org//charter/cat/[chartercalibrateddataset,{callid932}]/search?format=json&uid=call932_S1A_IW_SLC__1SDV_20211229T231926_20211229T231953_041230_04E66A_3DBE-calibrated",
   "secondary": "https://catalog.disasterscharter.org//charter/cat/[chartercalibrateddataset,{callid932}]/search?format=json&uid=call932_S1A_IW_SLC__1SDV_20220110T231926_20220110T231953_041405_04EC57_103E-calibrated",
   "aoi": "POLYGON((100.069 37.257,100.069 38.031,101.739 38.031,101.739 37.257,100.069 37.257))"
}

Visualization

See the result on the map. The preview appears within the area defined in the spatial filter.

Image-006

To get more information about the product just click on the preview in the map, a bubble showing the name of the layer “DInSAR Displacement Mapping from call932_S1A_IW_SLC__1SDV_20211229T231926_20211229T231953_041230_04E66A_3DBE-calibrated and call932_S1A_IW_SLC__1SDV_20220110T231926_20220110T231953_041405_04EC57_103E-calibrated” will appear and then click on the Show details button.

In the left panel under the result Details is possible to customize the overview on the fly by clicking on Layer Styling and Select Asset or Combine Assets. As an example under Select Assets you can select the asset displacement_c_vv_20211229_20220110 band change the color map and set the desired opacity level.

Image-007

It is also possible to make RGB composites on the fly using the single band assets generated by the service. To do so click on Combine Assets and specify the assets in the RGB channels.

Warning

When combining single-band assets under Combine assets please consider the valid range associated with the specific asset. As an example consider [0,1] for coherence and [-3.14,3.14] for the interferometric phase in radians. Valid ranges for all single-band products can be found here. Min and Max values shall be inserted manually in the Channel histogram min/max boxes.

Download

In this first use case the DInSAR-S1 service returned as output the following 6 products:

  • coh_c_vv_20210311_20210323: single-band asset coh_c_vv product for Coherence in C-band and VV polarization from the interferometric pair 20211229_20220110 as single band GeoTIFF in COG format,
  • phase_c_vv_20210311_20210323: single-band asset phase_c_vv product for the Interferometric Phase in C-band and VV polarization from the interferometric pair 20211229_20220110 as single band GeoTIFF in COG format,
  • displacement_c_vv_20210311_20210323: single-band asset displacement_c_vv product for LOS displacement in C-band and VV polarization from the interferometric pair 20211229_20220110 as single band GeoTIFF in COG format.
  • overview-coherence: multi-band visual asset derived from the Interferometric Coherence product using the "Gray" color map given as 4-band 8-bit GeoTIFF in COG format,
  • overview-phase: multi-band visual asset derived from the Interferometric Phase product in radians using the "Jet" color map given as 4-band 8-bit GeoTIFF in COG format,
  • overview-displacement: multi-band visual asset derived from the Line of Sight Displacement in m product using a "RdBu" color map given as 4-band 8-bit GeoTIFF in COG format,

To download these products click on the Download button located at the bottom of the Product Details tab in the left panel. Users can also dowloand the legend of each overview asset in PNG format.


Use case 2: DInSAR with SAOCOM-1 data

Select the processing service

After the opening of the activation workspace, in the right panel of the interface, open the Processing Services tab and select the processing service DInSAR Displacement Mapping - Stripmap (DInSAR-SM).

Image-008

The "DInSAR Displacement Mapping - Stripmap (DInSAR-SM)" panel is displayed with parameters values to be filled-in.

Image-009


Find the data using multiple filter criteria

Choose an area in which you want to focus your analysis. As an example near San Juan in Argentina (EO data credits: CONAE)1.

From the Navigation and Search toolbar (located in the upper left side of the map), click on “Spatial Filter” and draw a square around your Area of Interest. This spatial filter allows you to select only the EO data acquired over this area.

From the top of the left panel, use Filter Criterias to search for Radar and saocom-1 data collections from the list.

To further filter results by selecting only SLC products you can use the Text filter by inserting SLC in the Search Terms text box always visible on the top left of the map.

After the query, the Results list will include only the features satisfying all the above defined filters.

Now it is possible to choose a pair of "Reference" and "Secondary" SLC Datasets to be used for DInSAR. This pair must come from the same sensor and shall have the same radar geometry (incident angle, orbit path).

Fill the parameters

After the definition of all the filters, you can employ DInSAR by importing at least a pair of Sentinel-1 SLC Datasets.

Job name

Insert as job name:

DInSAR Stripmap SAO DEM Argentina CONAE pair

Input Reference and Secondary datasets

Drag and Drop in the Reference SLC dataset and in the Secondary SLC dataset fields the following datasets:

  1. [CD] SAOCOM-1A SLC SM VV/VH 2019-03-04 21:50:08,

  2. [CD] SAOCOM-1A SLC SM VV/VH 2019-02-16 21:50:07.

Image-010

!! Tip For convention the reference is the oldest SAR acquisition and you can check how the motion was evolved by inserting the last acquired image as secondary.

Note

This service supports only SAR complex imagery.

Warning

DInSAR-SM currently supports only SLC Datasets in Stripmap mode from the Saocom-1 EO mission.

Warning

Drag and drop the Dataset (e.g. [CD] SAOCOM-1A SLC SM VV/VH 2019-03-04 21:50:08) and not the single-band asset (e.g. amplitude-vv) into the Reference/Secondary SLC datasets fields.

Processing type

With this parameter you shall choose among the 2 possible processing types for the generation of the interferogram:

  1. DEM - a DEM assisted co-registration is employed for the interferogram generation,

  2. XCORR - the interferogram is generated without the DEM through a sequence of stacking, cross-correlation, and warping of reference and secondary.

Image-011

Leave the default option DEM.

Area of interest expressed as Well-known text

The “Area of interest as Well Known Text” can be defined by using the drawn polygon defined with the area filter.

Tip

In the definition of “Area of interest as Well Known Text” it is possible to apply as AOI the drawn polygon defined with the area filter. To do so, click on the button in the left side of the "Area of interest expressed as Well-known text" box and select the option AOI from the list. The platform will automatically fill the parameter value with the rectangular bounding box which is taken from the current search area in WKT format.

Note

This parameter is optional.

In this tutorial AOI is not defined.

Run the job

Once all parameters are inserted the processing service panel shall appear as the one shown in the figure below. Click on the button Run Job and see the Running Job. You can monitor job progress through the progress bar.

Image-012

Results: download and visualization

Once the job is completed, click on the button Show results at the bottom of the processing service panel.

Image-013

Tip

You can also save the parameters employed in this job by clicking on the Export params button in the right panel. This allows you to copy all your entries to the clipboard. This is meant to be used for a quick re-submission of a similar job after a fine tuning of the parameters (e.g. to add a colour formula later).

Below is reported the syntax which includes all the parameters employed in this example.

{
   "reference": "https://catalog.disasterscharter.org//charter/cat/[chartercalibrateddataset,!excluded,{callid969}]/search?format=json&uid=call969_S1A_OPER_SAR_EOSSP__CORE_L1A_OLF_20200221T122503-calibrated",
   "secondary": "https://catalog.disasterscharter.org//charter/cat/[chartercalibrateddataset,!excluded,{callid969}]/search?format=json&uid=call969_S1A_OPER_SAR_EOSSP__CORE_L1A_OLF_20200221T122606-calibrated",
   "processing_type": "DEM"
}

Visualization

See the result on the map. The preview appears within the area defined in the spatial filter.

Image-014

To get more information about the product just click on the preview in the map, a bubble showing the name of the layer “DInSAR Displacement Mapping from S1A_OPER_SAR_EOSSP__CORE_L1A_OLF_20200221T122503-calibrated and S1A_OPER_SAR_EOSSP__CORE_L1A_OLF_20200221T122606-calibrated” (Feature under the Results list with the tag DInSAR StripMap SAO DEM Argentina CONAE pair) will appear and then click on the Show details button.

In the left panel under the result Details is possible to customize the overview on the fly by clicking on Layer Styling and Select Asset or Combine Assets. As an example under Select Assets you can select another asset, as the overview-coherence one, and visualize it in the geobrowser together with the legend displayed in the bottom righter corner of the map.

Image-015

Same for a single band asset (e.g. by selecting coh_l_vv_20190216_20190304).

Image-016

It is also possible to make RGB composites on the fly using the single band assets generated by the service. To do so click on Combine Assets and specify the assets in the RGB channels.

Warning

When combining single-band assets under Combine assets please consider the valid range associated with the specific asset. As an example consider [0,1] for coherence and [-3.14,3.14] for the interferometric phase in radians. Valid ranges for all single-band products can be found here. Min and Max values shall be inserted manually in the Channel histogram min/max boxes.

Download

In this second use case the DInSAR-SM service returned as output the following 6 products:

  • coh_l_vv_20190216_20190304: single-band asset coh_l_vv product for Coherence in L-band and VV polarization from the interferometric pair 20190216_20190304 as single band GeoTIFF in COG format,
  • phase_c_vv_20190216_20190304: single-band asset phase_l_vv product for the Interferometric Phase in L-band and VV polarization from the interferometric pair 20190216_20190304 as single band GeoTIFF in COG format,
  • displacement_c_vv_20190216_20190304: single-band asset displacement_l_vv product for LOS displacement in L-band and VV polarization from the interferometric pair 20190216_20190304 as single band GeoTIFF in COG format.
  • overview-coherence: multi-band visual asset derived from the Interferometric Coherence product using the "Gray" color map given as 4-band 8-bit GeoTIFF in COG format,
  • overview-phase: multi-band visual asset derived from the Interferometric Phase product in radians using the "Jet" color map given as 4-band 8-bit GeoTIFF in COG format,
  • overview-displacement: multi-band visual asset derived from the Line of Sight Displacement in m product using a "RdBu" color map given as 4-band 8-bit GeoTIFF in COG format,

To download these products click on the Download button located at the bottom of the Product Details tab in the left panel.

Image-017

Users can also dowloand the legend of each overview asset in PNG format.


  1. CONAE, SAOCOM 1A Test Products for satellite information users, Interferometric pair - Stripmap Dual Pol products, Area of El Leoncito National Park, San Juan Province, Argentina. Available at: https://catalogos.conae.gov.ar