# Reflectivity clutter correction

**URL:** <https://openradar.discourse.group/t/reflectivity-clutter-correction/86>\
**Category:** Corrections/Data Quality\
**Created:** [November 3, 2022, 11:23am UTC](https://openradar.discourse.group/t/reflectivity-clutter-correction/86 "2022-11-03T11:23:52Z")\
**Posts on this page:** 2\
**Page:** 1

<div class="post-metadata">

**Author:** ![adithiyrnair](https://yyz2.discourse-cdn.com/free1/user_avatar/openradar.discourse.group/adithiyrnair/32/84_2.png) [@adithiyrnair](https://openradar.discourse.group/u/adithiyrnair)\
**Post date:** [November 3, 2022, 11:23am UTC](https://openradar.discourse.group/t/reflectivity-clutter-correction/86/1 "2022-11-03T11:23:52Z")

</div>

Hello,  
I am working on CFradial2 type nc radar data. I am actually trying to apply glabella approach clutter filter and vulpiani correction on my radar data to produce MAXZ plot. my radar data have DBZ, VEL, Spec. WIdth, ZDR,RHOHV and PHI.

Please help

---

<div class="post-metadata">

**Author:** ![syedhamidali](https://yyz2.discourse-cdn.com/free1/user_avatar/openradar.discourse.group/syedhamidali/32/377_2.png) [@syedhamidali](https://openradar.discourse.group/u/syedhamidali)\
**Post date:** [November 14, 2022, 7:32am UTC](https://openradar.discourse.group/t/reflectivity-clutter-correction/86/2 "2022-11-14T07:32:52Z")

</div>

Hi Adithiy,

You should be able to use `wradlib.io.open_cfradial2_dataset("RCTLS_30NOV2017_000100_L2B_STD.nc", group = "sweep_0")` to read the data, and then you can select the values of all the products, where PHIDP \< 179, or RHOHV \< 0.9.

```auto
vol2 = wrl.io.open_cfradial2_dataset("RCTLS_30NOV2017_000100_L2B_STD.nc", group = "sweep_0")
vol2.RHOHV.max().values
vol2.PHIDP.max().values, vol2.PHIDP.min().values, vol2.PHIDP.median().values, vol2.PHIDP.std().values
swp = vol2.where((vol2.PHIDP<=179) & (vol2.RHOHV<=0.99)).pipe(wrl.georef.georeference_dataset)
swp['x'] /= 1000
swp['y'] /= 1000
swp['z'] /= 1000
fig = plt.figure(figsize = [10,4])
ax = plt.subplot(121)
swp.PHIDP.plot(x="x", y = "y")
ax = plt.subplot(122)
swp.DBZ.plot(x="x", y = "y", vmin = -10, vmax = 70, cmap = cmap)

```

You can have a look at the below-attached notebook, where I tried pyart, wradlib and lrose as well.  
`https://github.com/syedhamidali/test_scripts/blob/test/for_Adithiy.ipynb`

Best,  
Hamid