Apply Radiometric Terrain Flattening (Image Analyst)

Available with Image Analyst license.

Summary

Corrects the input synthetic aperture radar (SAR) data for radiometric distortions due to topography.

Due to the side-looking nature of SAR sensors, features facing the sensor appear artificially brighter and features facing away from the sensor appear artificially darker. Radiometric terrain flattening normalizes the backscatter values so that value variations will be due to surface scattering properties.

Radiometric terrain flattening is necessary to obtain meaningful backscatter that can be related directly to the surface scattering properties of features in a SAR image over any terrain.

Usage

  • The input SAR data must be calibrated to beta nought.

  • Use the Apply Radiometric Calibration tool to calibrate the SAR data to beta nought.

  • This tool does not support a geodatabase as an output location.

  • If the input DEM does not span the entire SAR dataset, the tool will output NoData values for the pixels outside of the DEM extent for gamma nought and sigma nought.

  • The input DEM must be in the WGS84 (EPSG:4326) geographic coordinate system.

  • The geometric distortion mask is an optional output that can provide insight into the various geometric distortions caused by terrain in the input radar data. The output distortions provided are foreshortening, lengthening, layover, and shadow.

    Geometric distortions

    In the illustration above, the blue slope facing the sensor and the magenta slope facing away from the sensor are the same length on the ground, but in the SAR image, the blue foreshortening region appears to be shorter than the magenta lengthening region. This is because the slope facing away from the sensor has more pixels than the slope that is facing the sensor.

    Layover occurs when the radar signal reaches the top of a tall feature before it reaches the base. The green section of the steep mountain is an example of layover; it appears in the same pixel as the ground surface. The brown section is a combination of layover and shadow, which appears to the right of the upper section in the SAR image, despite being located to the left of it on the ground.

    Shadow occurs when an object blocks the radar signal. The yellow slope facing away from the sensor is not illuminated. Since radar illumination is not scattered in the atmosphere, the shadows in a SAR image appear black.

Parameters

LabelExplanationData Type
Input Radar Data

The input radar data.

The data that will have radiometric terrain flattening applied. The data must be radiometrically calibrated to beta nought.

Raster Dataset; Raster Layer
Output Radar Data

The radiometrically terrain-flattened radar data.

Raster Dataset
DEM Raster

The input DEM.

The DEM that will be used to estimate the local illuminated area and the local incidence angle.

Mosaic Layer; Raster Layer
Apply geoid correction
(Optional)

Specifies whether the vertical reference system of the input DEM will be transformed to ellipsoidal height. Most elevation datasets are referenced to sea level orthometric height, so a correction is required in these cases to convert to ellipsoidal height.

  • Checked—A geoid correction will be made to convert orthometric height to ellipsoidal height (based on EGM96 geoid). This is the default.
  • Unchecked—No geoid correction will be made. Use this option only if the DEM is provided in ellipsoidal height.

Boolean
Polarization Bands
(Optional)

The polarization bands that will be radiometrically terrain flattened.

The first band is selected by default.

String
Calibration Type
(Optional)

Specifies whether the output will be terrain flattened using sigma nought or gamma nought.

  • Gamma nought The beta nought backscatter will be corrected using an accurate computation of an area using a DEM. This is the default.
  • Sigma nought The beta nought backscatter will be corrected using the unit area of a plane that is locally tangent to the DEM.
String
Output Scattering Area
(Optional)

The scattering area radar dataset.

Raster Dataset
Output Geometric Distortion
(Optional)

The 4-band geometric distortion radar dataset. The first band is the terrain slope, the second band is look angle, the third band is the foreshortening ratio, and the fourth band is the local incidence angle.

Raster Dataset
Output Geometric Distortion Mask
(Optional)

The 1-band geometric distortion mask radar dataset. The pixels are classified using six unique values, one for each distortion type:

  • Undetermined Red, Green, Blue = 160, 177, 186—Value of 0
  • Foreshortening Red, Green, Blue = 0, 154, 222—Value of 1
  • Lengthening Red, Green, Blue = 175, 88, 186—Value of 2
  • Shadow Red, Green, Blue = 255, 198, 30—Value of 3
  • Layover Red, Green, Blue = 0, 205, 108—Value of 4
  • Layover and shadow Red, Green, Blue = 166, 118, 29—Value of 5

Raster Dataset

ApplyRadiometricTerrainFlattening(in_radar_data, out_radar_data, in_dem_raster, {geoid}, {polarization_bands}, {calibration_type}, {out_scattering_area}, {out_geometric_distortion}, {out_geometric_distortion_mask})
NameExplanationData Type
in_radar_data

The input radar data.

The data that will have radiometric terrain flattening applied. The data must be radiometrically calibrated to beta nought.

Raster Dataset; Raster Layer
out_radar_data

The radiometrically terrain-flattened radar data.

Raster Dataset
in_dem_raster

The input DEM.

The DEM that will be used to estimate the local illuminated area and the local incidence angle.

Mosaic Layer; Raster Layer
geoid
(Optional)

Specifies whether the vertical reference system of the input DEM will be transformed to ellipsoidal height. Most elevation datasets are referenced to sea level orthometric height, so a correction is required in these cases to convert to ellipsoidal height.

  • GEOIDA geoid correction will be made to convert orthometric height to ellipsoidal height (based on EGM96 geoid). This is the default.
  • NONENo geoid correction will be made. Use this option only if the DEM is provided in ellipsoidal height.
Boolean
polarization_bands
[polarization_bands,...]
(Optional)

The polarization bands that will be radiometrically terrain flattened.

The first band is selected by default.

String
calibration_type
(Optional)

Specifies whether the output will be terrain flattened using sigma nought or gamma nought.

  • GAMMA_NOUGHT The beta nought backscatter will be corrected using an accurate computation of an area using a DEM. This is the default.
  • SIGMA_NOUGHT The beta nought backscatter will be corrected using the unit area of a plane that is locally tangent to the DEM.
String
out_scattering_area
(Optional)

The scattering area radar dataset.

Raster Dataset
out_geometric_distortion
(Optional)

The 4-band geometric distortion radar dataset. The first band is the terrain slope, the second band is look angle, the third band is the foreshortening ratio, and the fourth band is the local incidence angle.

Raster Dataset
out_geometric_distortion_mask
(Optional)

The 1-band geometric distortion mask radar dataset. The pixels are classified using six unique values, one for each distortion type:

  • Undetermined Red, Green, Blue = 160, 177, 186—Value of 0
  • Foreshortening Red, Green, Blue = 0, 154, 222—Value of 1
  • Lengthening Red, Green, Blue = 175, 88, 186—Value of 2
  • Shadow Red, Green, Blue = 255, 198, 30—Value of 3
  • Layover Red, Green, Blue = 0, 205, 108—Value of 4
  • Layover and shadow Red, Green, Blue = 166, 118, 29—Value of 5

Raster Dataset

Code sample

ApplyRadiometricTerrainFlattening example 1 (Python window)

This example corrects a cross polarized radar dataset using a DEM and gamma nought parameters.

import arcpy
arcpy.env.workspace = r"C:\Data\SAR"
outRadar = arcpy.ia.ApplyRadiometricTerrainFlattening(
     "IW_manifest_TNR_CalB0_Dspk.crf", r"C:\Data\DEM\dem.tif", "GEOID", 
     "VH;VV", "GAMMA_NOUGHT")
outRadar.save("IW_manifest_TNR_CalB0_Dspk_RTFG0.crf")
ApplyRadiometricTerrainFlattening example 2 (stand-alone script)

This example corrects a cross polarized radar dataset using a DEM and gamma nought parameters.

# Import system modules and check out ArcGIS Image Analyst extension license
import arcpy
arcpy.CheckOutExtension("ImageAnalyst")
from arcpy.ia import *

# Set local variables
in_radar = r"C:\Data\SAR\IW_manifest_TNR_CalB0_Dspk.crf"
out_radar = r"C:\Data\SAR\IW_manifest_TNR_CalB0_Dspk_RTFG0.crf"
in_dem_raster = r"C:\Data\DEM\dem.tif"
ApplyGeoid = "GEOID" 
polarization = "VH;VV"
calibration_type = "GAMMA_NOUGHT"

# Execute 
outRadar = arcpy.ia.ApplyRadiometricTerrainFlattening(
     in_radar, in_dem_raster, ApplyGeoid, polarization, calibration_type)
outRadar.save(out_radar)

Related topics