[发明专利]一种雨杂波信号的处理方法、装置及设备

权利要求书

1.一种雨杂波信号的处理方法，其特征在于，所述方法包括：

获取雷达的回波信号；所述雷达的回波信号为网格化的回波数据，所述回波数据表示为s(m，n)，其中，m、n分别为所述回波数据矩阵的行数和列数；

对所述雷达的回波信号进行预处理，得到所述回波信号的时域回波信号；

对所述时域回波信号进行梯度计算处理，得到回波梯度；

根据所述回波梯度进行时域回波信号的连通域处理，得到连通域中的梯度方差；

根据所述连通域中的梯度方差和预设方差阈值，确定所述回波信号中的雨杂波信号；

其中，对所述雷达的回波信号进行预处理，得到所述回波信号的时域回波信号，包括：

对所述雷达的回波信号进行傅里叶变化得到频域回波信号S；

通过二阶高通滤波器对所述频域回波信号S进行高通滤波处理，得到滤波处理后的频域回波信号；

对所述滤波处理后的频域回波信号进行反傅里叶变换，得到所述滤波处理后的频域回波信号的时域回波信号 S^’；

通过二阶高通滤波器对所述频域回波信号S进行高通滤波处理，得到滤波处理后的频域回波信号，包括：

将所述频域回波信号S与二阶高通滤波器系数h相乘，得到滤波处理后的频域回波信号；

其中，所述二阶高通滤波器的系数可以通过以下公式获得：

，其中，h为二阶高通滤波器的系数，为d₀为滤波器截止频率，t为滤波器阶数，用于控制滤波边沿陡峭程度，d表示频谱图中频点位置至中心位置距离，由于d受数据矩阵维度影响，考虑自适应数据维度，设计数据维度缩放因子为zoomScaling=m/n，则d可以依据公式：获得，

其中m，n为所述雷达的回波信号对应的回波数据矩阵的行数和列数；i=1，2，3，...，m；j=1，2，3，...，n。

2.根据权利要求1所述的雨杂波信号的处理方法，其特征在于，对所述时域回波信号进行梯度计算处理，得到回波梯度，包括：

根据公式计算所述时域回波信号 的横向角度梯度分量；

根据公式计算所述时域回波信号 的纵向角度梯度分量；

根据所述横向角度梯度分量以及所述纵向角度梯度分量，获得所述时域回波信号 的回波梯度；

其中，m，n为所述时域回波信号 对应的回波矩阵数据的下标。

3.根据权利要求1所述的雨杂波信号的处理方法，其特征在于，根据所述回波梯度进行时域回波信号的连通域处理，得到连通域中的梯度方差，包括：

利用连通域统计方法统计所述回波梯度的多个梯度连通域bwlabel(L)，其中L为连通域数量；

遍历所述多个梯度连通域bwlabel(L)，获取所述多个梯度连通域bwlabel(L)中每个连通域的梯度方差Var(T)，其中T=1，2，3，...，L。

4.根据权利要求1所述的雨杂波信号的处理方法，其特征在于，根据所述连通域中的梯度方差和预设方差阈值，确定所述回波信号中的雨杂波信号，包括：

将所述连通域中的梯度方差小于或者等于预设方差阈值的回波信号确定为雨杂波信号。

5.根据权利要求4所述的雨杂波信号的处理方法，其特征在于，还包括：

当所述雨杂波信号大于检测门限时，对所述雨杂波信号进行置零。