Skip to content

HASARD-MM service specifications


icon


Service Description

The HASARD Multi Mission (HASARD-MM)_ service performs a mapping of standing waters from a co-pol backscatter image. It supports backscatter images from multiple missions. This service relies on the “HASARD”1 flood mapping algorithm that enables systematic, automatic and reliable SAR mapping of terrestrial water bodies. HASARD implements the Hierarchical Split-Based Approach (HSBA), described by Chini et al. 20172, and is specifically designed to support flooding-related disaster risk reduction at global scale. The underlying end-to-end image processing chain is based on a patented technology from the Luxembourg Institute of Science and Technology (LIST)3 that combines and fully integrates innovative hierarchical image splitting and open water backscatter modelling.

HASARD-MM exploits a single satellite image (Sigma Nought) to rapidly produce accurate standing water maps. The service employs Height Above the Nearest Drainage (HAND) auxiliary data to exclude false alarms in slope areas due to layover. HASARD-MM generates a water extent map in both raster and vector format.

Workflow

The HASARD-MM processing service applies the workflow below:

graph TB a[(Calibrated Datasets)] style a fill:#ffde86,stroke:#333,color:#282828,stroke-width:2px b[(Auxiliary Datasets)] style b fill:#ffde86,stroke:#333,color:#282828,stroke-width:2px a --> a1[input 1] a --> a2[input 2] b --> a3[input 3] b --> a4[input 4] subgraph Input style Input fill:#e8e8e8,stroke:#818181,color:#282828 a1[/Radar Dataset/] style a1 fill:#acc8ff,stroke:#333,color:#282828,stroke-width:2px a2[/Sigma0 asset name/] style a2 fill:#acc8ff,stroke:#333,color:#282828,stroke-width:2px a3[/HAND Dataset/] style a3 fill:#acc8ff,stroke:#333,color:#282828,stroke-width:2px a4[/HAND asset name/] style a4 fill:#acc8ff,stroke:#333,color:#282828,stroke-width:2px a5[/AOI as WKT/] style a5 fill:#acc8ff,stroke:#333,color:#282828,stroke-width:2px a6[/HSBA coefficients/] style a6 fill:#acc8ff,stroke:#333,color:#282828,stroke-width:2px a7[/HAND threshold/] style a7 fill:#acc8ff,stroke:#333,color:#282828,stroke-width:2px a8[/Backscatter threshold/] style a8 fill:#acc8ff,stroke:#333,color:#282828,stroke-width:2px end subgraph HASARD-MM a1--> aa1[Event backscatter <br> co-pol asset] a2--> aa1[Event backscatter <br> co-pol asset] style aa1 fill:#cfdfff,stroke:#333,color:#282828 a3--> aa2[HAND single <br> band asset] a4--> aa2[HAND single <br> band asset] style aa2 fill:#cfdfff,stroke:#333,color:#282828 aa1 --> stack[Co-located <br> stacking <br> and subset] aa2 --> stack[Co-located <br> stacking <br> and subset] a5 -.-> stack[Co-located <br> stacking <br> and subset] stack --> hasard[HSBA algorithm] a6 --> hasard[HSBA algorithm] a8 --> hasard[HSBA algorithm] hasard --> hand[HAND masking] a7 --> hand[HAND masking] hand --> cog[Convert to COG] cog --> p1(water mask) style p1 fill:#cfdfff,stroke:#87afff,color:#282828 p1 --> overview[Overview creation] overview --> p2(water mask overview) style p2 fill:#cfdfff,stroke:#87afff,color:#282828 p1 --> vec[Vectorize] vec --> p3(water mask vector) style p3 fill:#cfdfff,stroke:#87afff,color:#282828 p1 --> stac[Create STAC item] p2 --> stac[Create STAC item] p3 --> stac[Create STAC item] end subgraph Output style Output fill:#e8e8e8,stroke:#818181,color:#282828 stac --> o1[/Water mask<br>single band/] style o1 fill:#acc8ff,stroke:#87afff,color:#282828,stroke-width:2px stac --> o2[/Water mask<br>overview/] style o2 fill:#acc8ff,stroke:#87afff,color:#282828,stroke-width:2px stac --> o3[/Water mask<br>vector/] style o3 fill:#acc8ff,stroke:#87afff,color:#282828,stroke-width:2px end

Inputs

Input of the HASARD-MM service is a single band Sigma0 asset in dB from a Calibrated Dataset [CD] and the single band asset from an HAND Auxiliary Dataset. Thus to employ the service you need the following two Datasets:

  1. An event SAR Calibrated Dataset including a co-pol Sigma Nought asset in dB (e.g. s0_db_x_hh for TerraSAR-X) acquired during crisis,

  2. HAND Auxiliary Dataset including the hand single band asset providing relative elevations in m.

Warning

HASARD-MM requires as input the HAND auxiliary dataset. The hand single band asset must cover at least the area of interest defined within the footprint of the input sigma nought asset.

Warning

This service does not support SAR complex data (e.g. Sentinel-1 SLC).

Parameters

The HASARD-MM service requires a specified number of mandatory and optional parameters. All service parameters are listed in the below Table 1.

Parameter Description Required Default value
Radar calibrated dataset Input SAR Calibrated Dataset including an event Sigma0 single band asset providing backscatter values in db YES
Sigma nought asset name Name of the co-pol Sigma0 single band asset providing backscatter values in db (e.g. s0_db_x_hh for a Calibrated Dataset in X-Band) YES s0_db_c_vv
HAND auxiliary dataset Input Auxiliary Dataset including the hand single band asset providing relative elevations in meters YES
HAND asset name Name of the HAND single band asset providing relative elevations in meters YES hand
Area of Interest Area of interest expressed in WKT NO
Starting tiling level Hierarchical level of splitting. When the level is equal to 1 the image is split into 4 tiles, equal to 2 in 16 tiles, equal to 3 in 64 tiles, etc. YES 3
Check bimodal distribution (Ashman D) Ashman D coefficient to be used to check bimodal distribution YES 2.4
Minimum number of pixels per tile Minimum number of pixels per tile for checking the bimodality hypothesis YES 1000
Number of pixels for small object removal Number of pixels for small objects (blobs) removal YES 20
Threshold for HAND asset Relative elevation threshold to be considered when mapping water extent. Pixel classified as water that have height above the nearest drainage above this threshold in meters (e.g. 15m) are masked out. YES 15
Backscatter threshold Backscatter threshold to be considered when mapping water extent. The threshold must be given as a signed decimal in decibel (e.g. -20 db). YES -20

Table 1 - Service parameters for the HASARD-MM processor.

Input SAR calibrated dataset and asset name

First two mandatory parameters are dedicated to specify the input calibrated dataset and define which single band backscatter asset from this dataset shall be used for the water mapping computation. The system does not automatically select the co-pol asset from the given Calibrated Dataset to be used in the computation.

Warning

The drag and drop of the single-band asset (e.g. "s0_db_c_vv") is not possible. Users must drag and drop first a Calibrated Dataset (e.g. "[CD] SENTINEL-1A C-SAR IW VV/VH 134 2021/11/06 21:43:07") into the Radar calibrated dataset parameter and then specify the STAC key (name of the asset) in the Sigma nought asset name one.

Warning

Input single band asset specified into the Sigma nought asset name parameter must be a co-pol one. As an example from a full pol calibrated dataset in C-band you must insert either s0_db_c_hh or s0_db_c_hh as asset name.

Hint

To quickly find the names of all single band assets in a Calibrated Dataset, click on the Show assets button available near the feature title. A list with all single-band assets (CBNs) included within the Calibrated Dataset will appear under the feature title.

Image-note-001

Note

All CBNs available in the ESA Charter Mapper can be found here.

Input HAND auxiliary dataset and asset name

With the third HAND auxiliary dataset and fourth mandatory parameters, ‘HAND auxiliary datasetandHAND asset namerespectively, the user must specify the input auxiliary dataset and specify the name of the single band asset to be used for the water mapping computation. The name of the asset in the forth parameterHAND asset nameis pre-filled by default and the user must just drag and drop the auxiliary dataset in the third parameterHAND auxiliary dataset`.

Warning

The drag and drop of the single-band asset (e.g. hand) is not possible. Users must drag and drop the Auxiliary Dataset (e.g. "[CD] MERIT Hydro HAND - Hand") into the HAND auxiliary dataset parameter and then check that the STAC key (name of the asset) in the HAND asset name is set as hand.

Hint

Specifications of the hand single band asset can be found here.

AOI (optional)

The fifth parameter (optional) can be used to define the area of interest expressed as a Well-Known Text value.

Warning

The input area of interest must have an area higher than 1600 skm or 160000 ha. If a smaller area is inserted the job will not be executed and a warning will appear.

Tip

An AOI higher than 160000 ha can be easily drawn in the ESA Charter Mapper by using the Spatial Filter button available in the Navigation and Search Toolbar in the upper left side of the map. During the drawing of the rectangular area, a popup showing the size of the current area in ha will appear in the map next to the cursor.

Image-note-001

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 taken from the current search area in WKT format.

HSBA parameters

The following four mandatory parameters: Starting tiling level, Check bimodal distribution (Ashman D), Minimum number of pixels per tile, Number of pixels for small object removal represents the parameters of the retrieval algorithm. Default values given in Table 1 are the ones recommended for a general use.

HAND threshold

In the tenth mandatory parameter the user must specify a relative elevation threshold to be considered when mapping water extent. Pixels classified as water that have values of Height Above the Nearest Drainage (HAND) above this threshold in meters (e.g. 15m) are masked out. In particular, this threshold is used to derive a binary mask of potential floodplains from the HAND dataset which is then applied for masking the water extent.

Hint

To better define your HAND threshold, select under the Results tab in the left panel the HAND auxiliary dataset, open the Details of this auxiliary dataset, select the hand single band asset, and inspect its histogram. In the map you will see by default the HAND dataset with a grayscale color bar with the min max of the entire dataset. Under Expression insert the following

where((hand<=15),1,0);

and click on the underneath Apply button to visualize in the map the binary mask of potential floodplains applied in HASARD-CD. If you are not happy with the result, change the value of 15 in the expression with another one and click again on the Apply button. The map will be immediately refreshed with the updated binary map. Repeat this binarization on the fly of the HAND asset as much as you need. Once you are satisfied with the binary mask displayed in the map, you can replace the default value of 15 in the HAND threshold parameter with the one you identified.

Backscatter threshold

In the last mandatory parameter the user must insert a threshold value in decibel (e.g. -20 db) to be used for the computation of water extent from the backscatter signal.

Hint

To better define your backscatter threshold, select under the Results tab in the left panel the Radar Calibrated dataset that you defined as input of HASARD-CD (e.g. [CD] TSX-1 EEC SM HH 86 2012-09-26 17:41:16), open the Details of this calibrated dataset, select the s0_db_B_PP (where PP is the polarization [HH, HV, VH, VV] and B the SAR-band [X,C,L]). Inspect backscatter value using the Get Point Value function after a right click on the map. Employ the Get Point Value function over potential water pixel having low backscatter values. Define the suitable threshold by choosing a value in dB around the average of the left lobe values from the Sigma Nought histogram distribution. More information can be found in Chini et al., 20172.

Outputs

The HASARD-MM processor provides as output the following products:

  1. Water extent overview

  2. Water extent single band

  3. Water extent polygon vector

HASARD-MM Product Specifications can be found in the below tables.

Attribute Value / description
Long Name Water Extente RGB composite
Short Name overview-water-mask
Description Water bitmask in the blue channel including transparency
Data Type UnSigned 8-bit Integer
Band 4
Format COG
Projection UTM or EPSG:4326 - WGS84
Valid Range [1 - 255]
Fill Value 0
Attribute Value / description
Long Name Water mask from HASARD-MM
Short Name water-mask
Description Water extent bitmask: 1=Water, 0 =Not-water
Data Type 1-bit
Band 1
Format COG
Projection UTM or EPSG:4326 - WGS84
Valid Range [0 - 1]
Fill Value N/A
Attribute Value / description
Description Polygon vector converted from HASARD-MM water mask
Short Name water-mask-vector
Geospatial Data Type Vector
Vector Data Type Polygon
Feature properties DN values of source raster
Format Shapefile
Projection Native

  1. LIST, An automatic satellite image processing tool for flood hazard prevention on a global scale, available at www.list.lu

  2. M. Chini, R. Hostache, L. Giustarini and P. Matgen (2017), "A Hierarchical Split-Based Approach for Parametric Thresholding of SAR Images: Flood Inundation as a Test Case," in IEEE Transactions on Geoscience and Remote Sensing, vol. 55, no. 12, pp. 6975-6988, Dec. 2017, DOI: 10.1109/TGRS.2017.2737664

  3. LIST, HASARD developed by the Remote Sensing and Natural Resources Modelling group of the Environmental Research and Innovation Department, available at www.list.lu