Ground clutter characterization and elimination in mountainousterrain -M. Gabella and R. Notarpietro，2002

算法基础：地物和异常回波的空间变化程度大于气象回波。

数据准备：单仰角雷达反射率，也可拓展为组合反射率。

技术实现：

Step1：“spatial-proximity” filter 空间连续性滤波

选取5*5网格点，判断中间点和周围点的阈值（tr_var）差异，如异常给予异常标记；

tr_var选取：代表降水回波场的空间连续性，分析其多尺度的统计特征（例如空间傅里叶功率谱spatial Fourier power spectra，普适的结构函数generalized structure function，moment-scale analysis，Harris et al., 2001）

选取variogram-based by Germann and Joss (2001)，认为即使是1小时的中尺度对流系统，代表1km分辨率的雷达距离库Bin，周边分布8个像素点，其变化的反射率值不超过6dBZ。

采用5*5的网格点（好于3*3的网格点），选择合理的阈值(例如3 dBZ< tr1<9 dBZ, 6<np<10)对于滤波的结果不会剧烈起伏。（Np表示像素点多少）

Step2: test of compactness

判断临近的非零值的点，因为其可能是受旁边或对角的值（地物）影响。

比值R：选取范围内的总像素  /  目标像素及其评估范围内的像素数目的比值

地物的 R常接近于1，阈值 tr2 选择稍微大于1，可以去除地物杂波。

对于小的孤立的强回波单体Cell，算法会误判，所以阈值选取1.3，Np>11。

Clutter detection using the Gabella approach

http://wradlib.org/wradlib-docs/latest/notebooks/classify/wradlib_clutter_gabella_example.html