[发明专利]用于天气雷达图像径向干扰回波消除的剪切波方法

权利要求书

1.一种用于天气雷达图像径向干扰回波消除的剪切波方法，其特征是，对天气雷达图像中的径向干扰回波进行标识，得到标识矩阵，进而得到标识图像，通过RGB三个通道分别进行剪切波方法迭代处理，再进行整合，从而消除天气雷达图像中的径向干扰回波。

2.如权利要求1所述的用于天气雷达图像径向干扰回波消除的剪切波方法，其特征是，具体步骤如下：

步骤1：输入一幅含有径向干扰回波的天气雷达图像f；

步骤2：标识图像f中的径向干扰回波，得到标识矩阵M；

步骤3：将标识矩阵M作用于待处理图像f，按照公式f_m＝f*M进行计算，得到标识图像f_m；

步骤4：将标识后的彩色天气雷达图像分成RGB三个通道；

步骤5：分别选取R通道、G通道和B通道的图像，采用剪切波方法迭代；

步骤6：将分别处理后的RGB三个通道图像整合，输出处理后的图像。

3.如权利要求2所述的用于天气雷达图像径向干扰回波消除的剪切波方法，其特征是，其中步骤2的标识过程包括以下步骤：

步骤2-1：输入一幅含有径向干扰回波的天气雷达图像f；

步骤2-2：定义一个大小和图像f尺寸相等的标识矩阵M，将图像中的径向干扰回波部分标识为0，将图像中的其余部分标识为1,标识矩阵M的具体形式为：

其中下标i和j代表图像像素的位置；

步骤5的迭代过程包括以下步骤：

步骤5-1：输入R通道的待处理图像f_mR；

步骤5-2：设定初值和参数，设置迭代初值：f₀为原始待处理图像f_mR，初始循环次数n＝0；设置参数：阈值thr、剪切级数[sc(1),sc(2),sc(3),sc(4),sc(5)]、迭代次数N，其中，阈值用于图像剪切波变换硬阈值处理，剪切级数用于构建剪切波变换，迭代次数用于确定图像更新的次数；

步骤5-3：对图像f₀进行剪切波变换得到SH(f₀)，计算公式为：

SH(f₀)(a,s,t)＝<f₀,ψ_a,s,t>，

其中ψ_a,s,t为剪切波，具体形式为

其中为尺度矩阵，为剪切矩阵；

步骤5-4：对SH(f₀)进行硬阈值处理，得到G₀；

步骤5-5：对G₀进行剪切波逆变换，得到

步骤5-6：更新f_mR＝f_n+1，重复步骤5-3到步骤5-5；

步骤5-7：当n达到设置的最大迭代次数N时停止迭代，输出此时的f_mR；

步骤5-8：输入G通道的待处理图像f_mG，重复步骤5-2到步骤5-7；

步骤5-9：输入B通道的待处理图像f_mB，重复步骤5-2到步骤5-7。