[发明专利]一种用于合成孔径雷达图像的自适应滤波方法

权利要求书

1.一种用于合成孔径雷达SAR图像的自适应滤波方法，其特征在于，该方法包括如下步骤：

S1、对含斑点噪声的SAR图像数据进行对数变换，得到同态域含斑点噪声的SAR图像数据；

S2、对所述同态域含斑点噪声的SAR图像数据进行方位向和距离向的二维傅里叶变换，得到二维倒谱域含斑点噪声的SAR图像数据；

S3、利用二维矩形窗函数对所述二维倒谱域含斑点噪声的SAR图像数据进行滤波，得到二维倒谱域斑点噪声经抑制的SAR图像数据；

S4、对所述二维倒谱域斑点噪声经抑制的SAR图像数据进行方位向和距离向的二维傅里叶逆变换，得到同态域斑点噪声经抑制的SAR图像数据；

S5、对所述同态域斑点噪声经抑制的SAR图像数据进行指数变换，得到斑点噪声经抑制的SAR图像数据。

2.根据权利要求1所述的用于合成孔径雷达SAR图像的自适应滤波方法，其特征在于，所述步骤S3进一步包括如下子步骤：

S3.1、在二维倒谱域分别设定二维矩形窗函数中矩形窗的长度和宽度的初始值，并分别设定二维矩形窗函数中矩形窗的长度和宽度的收敛步长；

S3.2、利用矩形窗函数对所述二维倒谱域含斑点噪声的SAR图像数据进行滤波并计算滤波后图像数据的熵；

S3.3、分别按照二维矩形窗函数中矩形窗的长度和宽度的收敛步长减小二维矩形窗函数中矩形窗的长度和宽度；

S3.4、迭代步骤S3.2和步骤S3.3，直到滤波后图像数据的熵不再减小，得到二维倒谱域斑点噪声经抑制的SAR图像数据。

3.根据权利要求2所述的用于合成孔径雷达SAR图像的自适应滤波方法，其特征在于，所述步骤S3.1中二维矩形窗函数中矩形窗的长度和宽度的初始值分别设定为步骤S2中方位向傅里叶变换的维数的二分之一和距离向傅里叶变换的维数的二分之一。

4.根据权利要求2所述的用于合成孔径雷达SAR图像的自适应滤波方法，其特征在于，所述步骤S3.1中二维矩形窗函数中矩形窗的长度和宽度的收敛步长分别设定为步骤S2中方位向傅里叶变换的维数的十分之一和距离向傅里叶变换的维数的十分之一。

5.根据权利要求1所述的用于合成孔径雷达SAR图像的自适应滤波方法，其特征在于，步骤S2中对所述同态域含斑点噪声的SAR图像数据进行方位向和距离向的二维傅里叶变换的公式为：

G(k1,k2)=Σx=0N1-1Σy=0N2-1{log10[I(x,y)]}exp{-j2π(k1xN1+k2yN2)}=Σx=0N1-1Σy=0N2-1{log10[S(x,y)]}exp{-j2π(k1xN1+k2yN2)}+Σx=0N1-1Σy=0N2-1{log10[w(x,y)]}exp{-j2π(k1xN1+k2yN2)}]]>

公式中，x表示方位向，y表示距离向，I(x,y)为含斑点噪声的SAR图像数据，log₁₀(I(x,y))为同态域含斑点噪声的SAR图像数据，S(x,y)为不含斑点噪声的SAR图像数据，log₁₀(S(x,y))为同态域不含斑点噪声的SAR图像数据，w(x,y)为斑点噪声数据，log₁₀(w(x,y))为同态域斑点噪声数据，j表示虚数，k₁为与方位向x对应的方位向倒谱变量，k₂为与距离向y对应的距离向倒谱变量，N₁为方位向傅里叶变换的维数、N₂为距离向傅里叶变换的维数，G(k₁,k₂)为二维倒谱域含斑点噪声的SAR图像数据。