PSPHDW

Description

Estimate the power contributions for plate, helix, diplane and wire.

Parameters

psphdw(fili, filo)

Name	Type	Caption	Length	Value range
FILI *	str	Input fully polarimetric SAR image	1 -
FILO *	str	Output total power image	1 -

* Required parameter

Parameter descriptions

FILI

The name of the input polarimetric SAR data set. The input data set must contain either non-symmetrized or symmetrized fully polarimetric (quad-polarization complex) data in scattering, covariance, coherence, or Kennaugh format. The input can be either single-look or multi-look (multi-look is highly recommended). Multi-look can be achieved by applying a polarimetric filter (PSBOXCAR or PSPOLFIL) prior to computing the decomposition.

The input data must be a data set that has already been imported in PCIDSK (.pix) format by SARINGEST. Alternatively, it can be the key-file name of any GDB-supported POLSAR data set in its distribution format. For more information, and a complete list of supported POLSAR sensors and data products, see SARINGEST.

FILO

The name of the output file containing the results of the decomposition. The output file has the same dimensions as the input file, and four floating-point channels representing the contribution to the power from each of the four scattering mechanisms (plate, helix, diplane, and wire).

The specified file must not already exist.

Details

The total power for each pixel is decomposed into contributions from plate, helix, diplane, and wire.

Example

Create a decomposition-power image for a single-look RADARSAT-2 image that has already been imported in PCIDSK (.pix) format by SARINGEST.

from pci.psphdw import psphdw

fili	=	"rsat2_slc.pix"
filo	=	"rsat2_pow.pix"

psphdw( fili, filo )

Algorithm

Convert the input to symmetrized coherence matrix. Compute the helix (Ph) power as twice the magnitude of the imaginary part of the T23 element. The wire contribution (Pw) is computed as sqrt((4*real part of T12)**2 + T13**2). The surface contribution is computed as T11-Ph/2 and the double-bounce power (Pd) as T22+T33-Ph-Pw/2.