Degree of Polarization (dopdp)#

polsartools.dopdp(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]

Computes Barakat degree of polarization (DOP) from the input dual-polarization (dual-pol) C2 matrix data, and writes the output in various formats (GeoTIFF or binary). The computation is performed in parallel for efficiency.

Example:#

>>> dopdp("path_to_C2_folder", window_size=5, outType="tif", cog_flag=True)
This will compute Barakat degree of polarization from the C2 matrix in the specified folder,
generating output in Geotiff format with Cloud Optimized GeoTIFF settings enabled.

Parameters:#

infolderstr

Path to the input folder containing C2 matrix data.

window_sizeint, optional

Size of the processing window (default is 1).

outTypestr, optional

Output format of the files; can be “tif” (GeoTIFF) or “bin” (binary) (default is “tif”).

cog_flagbool, optional

If True, outputs Cloud Optimized GeoTIFF (COG) (default is False).

cog_overviewslist of int, optional

List of overview levels to be used for COGs (default is [2, 4, 8, 16]).

write_flagbool, optional

Whether to write the computed output files (default is True).

max_workersint, optional

Number of parallel workers for processing (default is None, which uses one less than the number of available CPU cores).

block_sizetuple of int, optional

Size of each processing block (default is (512, 512)), defining the spatial chunk dimensions used in parallel computation.

Returns:#

None

The function writes the computed Barakat degree of polarization to the specified output format.

Output Files:#

  • “dopdp.tif” or “dopdp.bin”

\[\text{DOP}_{dp}=\sqrt{1-\frac{4\times\text{det([C2])}}{\text{(Trace[C2])}^2}}\]

where, \(\text{[C2]}\) is co-variance matrix. Further details on the Barakat Degree of polarization can be found in [[10]](#10)

This Page