wradlib.dp.system_phidp_window

wradlib.dp.system_phidp_window#

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

Estimate the system differential phase (\(\Phi_{DP}^{sys}\)) from the range interval with the highest valid-data coverage.

For each ray, a rolling window of length rng is evaluated and the interval containing the maximum number of valid \(\Phi_{DP}\) bins is selected. The median \(\Phi_{DP}\) within this interval 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.

Unlike system_phidp_block, this method does not require all bins within the selected interval to be valid. It therefore provides estimates for rays that contain gaps in \(\Phi_{DP}\) coverage.

Parameters:

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

  • rng (float) – Physical range length (m) used to determine the rolling window size.

  • n_lowest_rays (int, optional) – Number of lowest ray-wise 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 interval.

  • stop_range : Stop range of the selected interval.

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

Notes

The selected interval maximizes the count of valid \(\Phi_{DP}\) bins within the specified window length. No minimum coverage threshold is applied.

Examples

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