wradlib.dp.system_phidp_block

wradlib.dp.system_phidp_block#

wradlib.dp.system_phidp_block(phidp, rng, n_lowest_rays=30)[source]#

Estimate the system differential phase (\(\Phi_{DP}^{sys}\)) from contiguous valid \(\Phi_{DP}\) segments along radar rays.

The algorithm searches each ray for a sequence of N consecutive valid \(\Phi_{DP}\) bins, where N is derived from the requested range length (rng). For each ray, the median \(\Phi_{DP}\) within the identified segment is used as a ray-wise \(\Phi_{DP}^{sys}\) estimate. The final \(\Phi_{DP}^{sys}\) estimate is computed as the median of the n_lowest_rays smallest ray-wise estimates.

Parameters:

  • phidp (xarray.DataArray) – Differential phase field containing a range dimension.

  • rng (float) – Range length (m) used to determine the required number of consecutive valid bins.

  • n_lowest_rays (int, optional) – Number of lowest ray-wise \(\Phi_{DP}\) estimates used to derive the final \(\Phi_{DP}^{sys}\) estimate. Default is 30.

Returns:

xarray.Dataset – Dataset containing:

  • sysphi_ray : ray-wise \(\Phi_{DP}^{sys}\) estimate.

  • sysphi : global \(\Phi_{DP}^{sys}\) estimate.

  • start_range : Start range of the selected valid segment.

  • stop_range : Stop range of the selected valid segment.

  • valid_bins : Number of valid \(\Phi_{DP}\) bins within the selected interval.

Notes

Rays that do not contain N consecutive valid \(\Phi_{DP}\) bins receive NaN values for the corresponding outputs.

Examples

See Core Features - Dual-Pol - System Differential Phase.