Step by Step: Recommended Practice Flood Mapping
Step by Step: Recommended Practice Flood Mapping jan.stefanskiSAR data from Sentinel-1, a SAR mission from ESA, is available free of charge. To download Sentinel-1 data, register at the Copernicus Data Space Ecosystem. Then a search option becomes available which can be used to specify the data need such as region of interest, product type, sensor mode, sensing period, among others. In the example below, we used Level-1 Ground Range Detected (GRD) Sentinel-1 data, which incorporates already some basic preprocessing.
Other SAR imagery is not freely available; these include: Radarsat-2, TerraSAR-X, and Cosmo-SkyMed. But archived SAR images from past missions can be freely obtained. For example, SAR images from Envisat/ASAR can be obtained from the ESA Category-1 programme (http://eopi.esa.int/). Also, for commercial images (Radarsat-2, TerraSAR-X, and Cosmo-SkyMed) there are opportunities for scientific use of data for free or reduced rates.
Processing and screenshots are based on ESA's Sentinel Application Platform (SNAP) version 8.0 64-bit.
Step 1: Pre-processing - Calibration
Step 2: Pre-processing - Speckle filtering
Step 4: Post-processing - Geometric correction
SNAP supports all major data formats that are used to store SAR images.
0.1 To read the data set point to File --> Open Product.
0.2 Select the Sentinel-1 image, either directly the *.zip file or the manifest.safe from the unzipped folder.
0.3 The Product Explorer on the left shows relevant information on the product. This includes: Metadata (different SAR parameters on the orbit and image); Tie-points grid (interpolated latitude, longitude, incident angle and slant range time values); Bands (actual image bands). By right clicking on the Product, Properties can be opened which include information on mission, acquisition date, pass, etc.
0.4 For each polarization recorded there are two bands: Amplitude and Intensity. (The Intensity band is a virtual one. It is the square of the amplitude). Double click on either Amplitude or Intensity to view the image. In the bottom left corner, the WorldWind View shows the footprint of the selected scene.
0.5 If you need a subset, select on the Menu panel Raster --> Subset. Specify parameters of the region of interest and click OK button. SNAP will open the product with all metadata relevant to the subset.
0.6 The created subset is added as a new product, which is handled simultaneously as described in step 0.4. The footprint of the subset is also shown in the WorldWind View.
1. Pre-processing - Calibration
1.1 Calibration: select on the Menu panel Radar --> Radiometric --> Calibrate.
1.2 A Calibration window will open. Select the tab: Processing Parameters. Select the Polarizations you wish to process. This will create a new product with calibrated values of the backscatter coefficient.
2. Pre-processing - Speckle filtering
2.1 Next step is filtration. Select Radar --> Speckle Filtering --> Single Product Speckle Filter.
2.2 A Speckle filtering window will appear. Select the tab: Processing Parameters. Select Sigma0_VV and select Lee filter with filter size 7 by 7. Click the Run button.
2.3 A new product will be created in the Product Explorer. Open the band in the newly created product.
2.4 Compare the non-filtered and filtered versions.
3.1 To separate water from non-water a threshold can be selected. For this, we will analyse the histogram of the filtered backscatter coefficient. On the left side panel select the Colour Manipulation tab. The histogram of the backscatter coefficient will show up and one might need to use the logarithmic display. The histogram will show one or more peaks of different magnitude depending on the data. Low values of the backscatter will correspond to the water, and high values will correspond to the non-water class. We need to select the value that will separate water from non-water. In our case, the threshold value will be 2.22E-2.
3.2 To segment or binarize the image we will apply band arithmetic. For this, go to Raster --> Band Maths…
3.3 A window will open. Type a name for the new band, for example, water. Remove checkbox from the virtual band (the virtual band is in the memory but not physically on the disc). Click Edit Expression… button. The band math utility allows one to make any sorts of mathematical and logical expressions on the bands available in the product. Our goal is to create a new image which will be, for example, 255, for water objects (for backscatter lower our threshold of 2.22E-2), and 0 for higher values. To do this, enter into the Expression field the following expression: 255*(Sigma0_VV<2.22E-2). The expression (Sigma0_VV<2.22E-2) will return the logical value: true (or 1) for values less than 2.22E-2, and false (or 0) for higher values. Then we just multiply by 255. Click OK button.
Alternatively, another Band Math equation could be
if Sigma0_VV < 2.22E-2 then 1 else 0
3.4 A new band named water will be added to the product. Double click this band to view the image.
4. Post-processing - Geometric correction
4.1 The obtained image is in the geometry of the sensor. We need to reproject it to the geographic projection. To do this, select Radar --> Geometric --> Terrain Correction --> Range-Doppler Terrain Correction.
4.2 A window with parameters will open. Select the Processing Parameters tab. In the Source Bands select only water; Digital Elevation Model - SRTM 3Sec (Auto Download), the DEM over the region that SAR image covers will be automatically downloaded, use SRTM 1Sec (Auto Download) if the tool does not run properly; DEM Resampling Method - BILINEAR_INTERPOLATION; Image Resampling Method - NEAREST_NEIGHBOUR; Pixel Spacing - 10m (which depends on the sensor and its acquisition mode); Map projection - WGS84(DD), if UTM coordinates are needed, use UTM/WGS84 (Automatic), SNAP will automatically select the zones. Click OK button.
4.3 A new product will be created and appear in the Product Explorer. The product will be created in the DIM format and the image will be in img format. Double click the band in the newly created product. Notice that the geometry of the image changed to the geographic projection. This file can now be opened in a GIS (for example QGIS) to visualize and create a map.
5. Visualization in Google Earth
5.1 In order to display solely the water pixels, go to the left side panel, select the Colour Manipulation tab and the Table editor. Click on the boxes under Colour and set the background colours to None to make it transparent; optionally, you can change the water pixels (white colour) to another colour.
5.2 To visualize the water band in Google Earth, export it as a KMZ file; right click on the opened image --> Export View as Google Earth KMZ. Another way to export the flood mask is through File --> Export --> select file format to export.
5.3 Import the KMZ file to Google Earth.
Example of the resulting product
To assess the quality of the resulting flood mask, we compare it with the flood mask created by the Copernicus Emergency Management Service (EMS). The Copernicus EMS website for the flood event in Malawi (January 2015) can be found here and includes the results in vector format. These files can be imported and analyzed in a GIS (e.g., Quantum GIS).
The following figure shows both flood masks in comparison. Due to the thresholding without manual post-editing or other post-processing techniques in the Recommended Practices Flood Mapping, several pixels around the main water bodies are found as water, whereby most of them might be false positives (a wrong classification of water). Apart from that and some river structures, there is an overall concordance between both flood masks.
In the next figure Sentinel-1 subsets from before and during the flood event as well as the Recommended Practies and EMS flood masks are shown in detail. As an example, the red markings indicate false negatives meaning that water pixels were not classified as water. The orange marking of the EMS flood mask indicates an example for a false positive - a wrong classification of water.
