[发明专利]基于非下采样轮廓波的合成孔径雷达图像去噪方法

权利要求书

1.一种基于非下采样轮廓波的合成孔径雷达图像去噪方法，包括如下步骤：

(1)输入SAR图像X，并对其进行L层非下采样轮廓波变换，得到低频子带D_L(a，b)和不同尺度上的高频方向子带{C_l，i(a，b)，0≤l≤L-1，1≤i ≤k_l}，k_l表示图像在尺度2^-l上的高频方向子带的数目，C_l，i(a，b)表示SAR图像在尺度2^-l上的第i个高频方向子带的系数，L为3～5；

(2)计算不同尺度各个高频方向子带的斑点噪声方差 

2a)计算对SAR图像X进行非下采样轮廓波变换时，施加在不同尺度各个高频方向子带的等价滤波器操作算子的二阶累积量 

其中， 表示第i级塔形滤波的非下采样拉普拉斯滤波器系数， 表示第j级方向滤波的非下采样方向滤波器系数；

2b)计算不同尺度各个高频方向子带的均值μ_X和方差 

2c)计算不同尺度各个高频方向子带的斑点噪声方差 

其中，G_Z是输入SAR图像X的规范化标准差，且对于强度SAR图像有 对于幅度SAR图像有 R为已知输入SAR图像的视数，R≥1；

(3)利用高频方向子带系数C_l，i(a，b)的局部均值mean[C_l，i(a，b)]和高频方向子带的斑点噪声方差 将高频方向子带系数C_l，i(a，b)区分为信号或是噪声的变换系数，判定准则如下：

其中，k在1～5之间取值，mean[C_l，i(a，b)]是以高频方向子带系数C_l，i(a，b)为中心的5×5大小的局部区域内系数的均值；

(4)对判定后的高频方向子带系数C_l，i(a，b)进行如下去噪操作：

得到去噪后的高频方向子带系数 

(5)对步骤(1)中得到的低频子带D_L(a，b)和去噪后的高频方向子带系数 作非下采样轮廓波逆变换，得到去噪后的SAR图像X^*。

2.一种基于非下采样轮廓波的合成孔径雷达图像去噪虚拟装置，包括：

变换装置，用于输入SAR图像X，并对其进行L层非下采样轮廓波变换，得到低频子带D_L(a，b)和不同尺度上的高频方向子带{C_l，i(a，b)，0≤l≤L-1，1≤i≤k_l}，C_l，i(a，b)表示SAR图像在尺度2^-l上的第i个高频方向子带的系数，L为3～5，k_l表示图像在尺度2^-l上的高频方向子带的数目；

计算装置，用于计算不同尺度各个高频方向子带的斑点噪声方差 

判定装置，用于利用高频方向子带系数C_l，i(a，b)的局部均值mean[C_l，i(a，b)]和高频方向子带的斑点噪声方差 将高频方向子带系数C_l，i(a，b)区分为信号或是噪声的变换系数，判定准则如下：

其中，k在1～5之间取值，mean[C_l，i(a，b)]是以高频方向子带系数C_l，i(a，b)为中心的5×5大小的局部区域内系数的均值；

去噪装置，用于对判定后的高频方向子带系数C_l，i(a，b)进行如下去噪操作：

得到去噪后的高频方向子带系数 

逆变换装置，用于对低频子带D_L(a，b)和去噪后的高频方向子带系数 作非下采样轮廓波逆变换，得到去噪后的SAR图像X^*；

所述的计算装置，进一步包括：

二阶累积量计算子装置，用于计算对SAR图像X进行非下采样轮廓波变换时，施加在不同尺度各个高频方向子带的等价滤波器操作算子的二阶累积量 

其中， 表示第i级塔形滤波的非下采样拉普拉斯滤波器系数， 表示第j级方向滤波的非下采样方向滤波器系数；

高频方向子带均值和方差计算子装置，用于计算不同尺度各个高频方向子带的均值 μ_X和方差 

斑点噪声方差计算子装置，用于计算不同尺度各个高频方向子带的斑点噪声方差 

其中，G_Z是输入SAR图像X的规范化标准差，且对于强度SAR图像有 对于幅度SAR图像有 R为已知输入SAR图像的视数，R≥1。