Clutter detection - Satellite Cloud Images

Contents

Clutter detection - Satellite Cloud Images#

In this notebook we show howto identify clutter based on collocated cloudtype.

import matplotlib.pyplot as plt
import numpy as np
import wradlib as wrl
import wradlib_data
import xarray as xr
import xradar as xd
from IPython.display import display
from osgeo import osr

/home/docs/checkouts/readthedocs.org/user_builds/wradlib/conda/latest/lib/python3.13/site-packages/tqdm/auto.py:21: TqdmWarning: IProgress not found. Please update jupyter and ipywidgets. See https://ipywidgets.readthedocs.io/en/stable/user_install.html
  from .autonotebook import tqdm as notebook_tqdm

Estimate localisation errors#

timelag = 9 * 60
wind = 10
error = np.absolute(timelag) * wind

Plot the results#

fig = plt.figure(figsize=(16, 8))

tilt = 0

ax = fig.add_subplot(131)
pm = vol.DBZH[tilt].wrl.vis.plot(x="xp", y="yp", ax=ax)
# plt.colorbar(pm, shrink=0.5)
plt.title("Radar reflectivity")

ax = fig.add_subplot(132)
pm = vol.CT[tilt].wrl.vis.plot(x="xp", y="yp", ax=ax)
# plt.colorbar(pm, shrink=0.5)
plt.title("Satellite cloud classification")

ax = fig.add_subplot(133)
pm = vol.CMAP[tilt].wrl.vis.plot(x="xp", y="yp", ax=ax)
# plt.colorbar(pm, shrink=0.5)
plt.title("Detected clutter")

fig.tight_layout()

../../../_images/b95a433832d5926407b4941bebad4840fadc96915b1ac5d58c831b4c24880768.png