Polarimetric Radar Vegetation Index (prvifp )#

polsartools.prvifp(infolder, window_size=1, outType='tif', cog_flag=False, cog_overviews=[2, 4, 8, 16], write_flag=True, max_workers=None, block_size=(512, 512), progress_callback=None)[source]

Calculate Polarimetric Radar Vegetation Index (PRVI) from full-pol SAR data.

This function computes the Polarimetric Radar Vegetation Index (PRVI) using full-polarimetric SAR data. PRVI is specifically designed to assess vegetation density and biomass by utilizing the complete polarimetric information available in the coherency (T3) or covariance (C3) matrix.

Examples

>>> # Basic usage with default parameters
>>> prvifp("/path/to/fullpol_data")
>>> # Advanced usage with custom parameters
>>> prvifp(
...     infolder="/path/to/fullpol_data",
...     window_size=5,
...     outType="tif",
...     cog_flag=True,
...     block_size=(1024, 1024)
... )
Parameters:
  • infolder (str) – Path to the input folder containing full-pol T3 or C3 matrix files.

  • window_size (int, default=1) – Size of the spatial averaging window. Larger windows reduce speckle noise but decrease spatial resolution.

  • outType ({'tif', 'bin'}, default='tif') – Output file format: - ‘tif’: GeoTIFF format with georeferencing information - ‘bin’: Raw binary format

  • cog_flag (bool, default=False) – If True, creates a Cloud Optimized GeoTIFF (COG) with internal tiling and overviews for efficient web access.

  • cog_overviews (list[int], default=[2, 4, 8, 16]) – Overview levels for COG creation. Each number represents the decimation factor for that overview level.

  • write_flag (bool, default=True) – If True, writes results to disk. If False, only processes data in memory.

  • max_workers (int | None, default=None) – Maximum number of parallel processing workers. If None, uses CPU count - 1 workers.

  • block_size (tuple[int, int], default=(512, 512)) – Size of processing blocks (rows, cols) for parallel computation. Larger blocks use more memory but may be more efficient.

Returns:

Writes one output file to disk: - ‘prvi_fp.tif’ or ‘prvi_fp.bin’: PRVI values

Return type:

None

The formlation of PRVI interms of degree of polarization and cross-pol backscatter intensity can be expressed as follows:

\[\text{PRVI}_{fp}=(1-\text{DOP}_{fp})\sigma^\circ_{\text{XY}}\]

where, \(\text{DOP}_{fp}\) 3D Barakt degree of polarization and can be expressed as shown below. Further details on the PRVI can be found in [[1]](#1)

\[\text{DOP}_{fp}=\sqrt{1-\frac{27\times\text{det([T3])}}{\text{(Trace[T3])}^3}}\]