Reclassify RTC Tool

Generate a water mask using a dB threshold approach

Calm surface water is generally characterized by very low radiometric returns in Synthetic Aperture Radar (SAR) imagery. This tendency can be used to identify surface water in a Radiometric Terrain Corrected (RTC) SAR product.

Case Study Application

The RTC products used in this demonstration were generated using ASF's  On Demand Processing  functionality at  Data Search - Vertex . Click the button below to search for the same two source granules used in the demonstration.

In Michigan, heavy rainfall in May 2020 resulted in flooding in Shiawassee National Wildlife Refuge. Because of the greater differentiation among dark pixels that is visible in dB scale, we will use RTC images in dB scale to determine surface water extent after flooding.

Powered by Esri
Powered by Esri

Sentinel-1 RTC imagery in dB scale before (left) and after (right) the flood. Move the bar to compare the before and after images. RTC products processed by ASF DAAC HyP3 2020 using GAMMA software. Contains modified Copernicus Sentinel data 2020, processed by ESA.

If you did not order your On Demand RTC products in dB scale, refer to the  tutorial for the Scale Conversion Tool  to learn how to convert RTC products from amplitude or power to dB scale.

Determine Appropriate Threshold Value

There is not one constant threshold value that can be applied to any dB image to identify surface water. The value will depend on the local site conditions at the time the image was required. To determine the appropriate value to apply, start by simply applying a classified symbology to the dB image.

In this example, we will use the co-polarized (VV in this case) image to generate the water mask. On a calm day, surface water will return very little VV backscatter to the sensor, because the surface of the water is smooth and water has a high dielectric constant. As a result, most of the incoming radar signal bounces off the surface and off in the opposite direction of the sensor.

On a windy day, however, the water surface can become rough enough to return some of the signal to the sensor, and it becomes much more difficult to distinguish between the water and surrounding terrain. If you must classify the water using an image collected on a windy day, you may be better off using the cross-pol (VH or HV) dataset, or a combination of the VV and VH values.

Step through the slideshow below to learn how to determine the threshold value to use for classification of surface water.

Add dB Raster to Map

Drag the dB raster onto the map to add it. Build pyramids if prompted.

Change Symbology to Classified

Right-click the dB raster and select Symbology.

Change the primary symbology from Stretch to Classify.

Set the Number of Classes to 2

We will want to have only two classes: one for water, one for everything that's not water.

Set the Non-water Color to Transparent

Click the color patch next to the second class and select No color.

Adjust the Threshold

Click on the Histogram tab, and move the arrow up and down to adjust the threshold value.

Images with a lot of surface water often have 2 peaks in the histogram with the first peak consisting mostly of water pixels. The trough between the two peaks is a good place to start. Adjust the value as necessary to capture the water extent.

Set Threshold Value

Once you determine which value best represents the water on the landscape, you can either keep the value from the Histogram, or set a specific value in the Classes tab.

Click in the cell for the Upper value number, and enter the desired break point. Click away from the cell or press tab to accept the new value.

You can check your extent map against a variety of other datasets. It can be helpful to compare the classified layer against the original dB layer in grayscale stretch. In this example, the water mask is applied to all pixels with values less than -18 dB.

Sentinel-1 RTC image in dB scale (left) and the same raster symbolized to show values less than -18 in yellow and all other pixels as transparent. Drag the swipe bar to compare the mask with the dark pixels. RTC product processed by ASF DAAC HyP3 2020 using GAMMA software. Contains modified Copernicus Sentinel data 2020, processed by ESA.

Another potential resource to use for comparison is an  RGB Decomposition  image, which combines the co- and cross-pol datasets into a false color image. Pixels with low backscatter in both co- and cross-pol datasets appear as dark blue in the decomposition. There is an RGB Decomposition Tool included in the ASF_Tools ArcGIS Toolbox. In this example, the water mask is applied to all pixels with values less than -19 dB.

Sentinel-1 RGB Decomposition of RTC data under a pink water mask generated by including all co-pol (VV) RTC pixel values less than -19 dB. Drag the swipe bar to compare the mask with the very dark blue pixels. RTC products processed by ASF DAAC HyP3 2020 using GAMMA software. Contains modified Copernicus Sentinel data 2020, processed by ESA.

Once you determine the threshold value that will work best for your specific application, you could just save the newly symbolized raster as a layer and use it for a visual overlay. In some cases, however, you may want a reclassified raster indicating surface water extent to use for analysis with other data layers.

The RTC Reclassify Tool allows you to generate a new mask raster that only includes the pixels that fall below the threshold value.

Using the Reclassify RTC Tool

The RTC Reclassify Tool in the ASF_Tools ArcGIS Toolbox generates a raster mask based on a threshold value. Any pixels in the input raster with a value below the threshold are classified with a value of 1, and all other pixels are set to No Data. It was designed to be used to generate water masks, based on a dB scale threshold value, but will work for any reclassification where the values of the pixels of interest are below a threshold value.

You can  download the toolbox  from ASF's website. Extract the zip file to the directory of your choice on your computer.

Step through the slideshow below to learn how to use the tool in ArcGIS Pro. This tool can also be used in ArcMap; the process is generally the same, but the interface looks different.

Open ArcGIS Pro and Navigate to the Toolbox

In the Catalog pane, navigate to the directory containing the toolbox, and expand the toolbox until all of the tools are displayed. Double-click on the Reclassify RTC tool to launch the Geoprocessing dialog.

Select Raster to Convert

Click the browse icon next to the first parameter field and navigate to the raster that you want to reclassify.

Set the Threshold Value

Enter the threshold value. Note that dB threshold values will often be negative; be sure to include the - sign before the number.

Run the Process

Change the default destination folder or output filename if desired.

Click the Run button to process the data. By default the new raster is added as a layer in the map.

When using the tool in ArcGIS Pro, hover over any of the blue information icons in the Geoprocessing dialog for more information about the parameters or the tool in general. The full tool information can be viewed by right-clicking the tool in the Catalog window and selecting View Metadata.

When using the tool in ArcMap, click the Show Help button in the tool dialog box to see the information about specific parameters or the tool in general in the Help pane. The full tool information can be viewed by right-clicking the tool in the Catalog window and selecting Item Description.

Important Considerations

It is important to be aware that not all surfaces with very low backscatter are water. While this approach can work well in many situations, make sure to consider the context of the image before relying on the masks.

In particular, both calm surface water and dry desert sand tend to have very low backscatter, but for opposite reasons. The high dielectric constant and smooth surface of the water result in the signal reflecting off the surface, essentially bouncing off in the opposite direction from the incoming signal.

Dry sand, on the other hand, has very low dielectric properties, allowing the signal to penetrate deeply into desert landscapes. By the time the signal is scattered, it lacks the energy necessary to return back through the medium. In both cases, there are low radiometric returns, but one is because the signal was bounced away, while the other is because the signal was absorbed by the medium.