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.

  • 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.

Parameters

LabelExplanationData Type
Input Radar Data

The input radar data.

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 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 expressed 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
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

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 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 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 expressed 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 = "D:\Data\SAR\S1\20181014"
outRadar = arcpy.ia.ApplyRadiometricTerrainFlattening(
     "IW_manifest_TNR_CalB0_Dspk.crf", "D:\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 = "D:\Data\SAR\S1\20181014\IW_manifest_TNR_CalB0_Dspk.crf"
out_radar = "D:\Data\SAR\S1\20181014\IW_manifest_TNR_CalB0_Dspk_RTFG0.crf"
in_dem_raster = "D:\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)

Licensing information

  • Basic: Requires Image Analyst
  • Standard: Requires Image Analyst
  • Advanced: Requires Image Analyst

Related topics