SARTOPOCAL

Apply radiometric correction to compensate for terrain-slope effects


EnvironmentsPYTHON :: EASI
Quick linksDescription :: Parameters :: Parameter descriptions :: Return Value :: Details :: Example :: Related

Back to top

Description


SARTOPOCAL rescales synthetic aperture radar (SAR) intensities to compensate for variations in backscattering due to slope effects.
Back to top

Parameters


sartopocal(fili, filedem, filo)

Name Type Caption Length Value range
FILI * str Name of input SAR reference file 1 -   1 -
FILEDEM * str File name of DEM 1 -    
FILO * str Name of output corrected file 1 -   1 -

* Required parameter
Back to top

Parameter descriptions

FILI

The name of the input file containing the SAR data to use as reference.

The input can be in slant range or ground range.

FILEDEM

The name of the input file containing elevation data.

The elevations are resampled to match the input sample size, Earth model (usually WGS-84), and processing elevation of the SAR data set. If the DATUM metadata tag is not specified, mean sea level (MSL), in meters, is assumed.

FILO

The name of the output file containing the radiometrically corrected data.

The file name you specify must not already exist. The output file is written in the same projection (slant or ground range) as the input file and contains the same number of image rows, columns, and channels. Pixels that do not overlap with the DEM are marked as NODATA.

Back to top

Return Value

Returns: Program execution status

Type:  PCI_INT

The return value is 0.

This function returns only if it executes successfully; otherwise, it throws an exception.

Back to top

Details

SARTOPOCAL rescales SAR data to compensate for variations in backscattering due to slope effects. The radiometric correction is determined by combining the elevations of the resampled digital elevation model (DEM) with the known viewing geometry, nominal processing elevation, and Earth model (usually WGS-84) of the input SAR data set. Triangular terrain facets are projected to the plane perpendicular to the local slant range direction to compute the intensity due to slope effects. The original input values are normalized accordingly. If the input contains complex-valued data, the magnitude (amplitude) is rescaled, but the phase value remains unaltered. Areas not covered by the DEM are marked as NODATA.

Back to top

Example

In the following example, the input data is rescaled to compensate for terrain effects caused by the SAR viewing geometry. The input DEM file is converted automatically from MSL to WGS-84.

Note: The radiometric calibration type of the output need not match the calibration type of the input.
from pci.sartopocal import sartopocal
fili="beta_flat_earth_radiometry.pix"
filedem="n38w123.hgt"
filo="sigma_terrain_corrected_radiometry.pix"
sartopocal(fili, filedem, filo)
      

© PCI Geomatics Enterprises, Inc.®, 2026. All rights reserved.