[发明专利]基于NSCT域局部高斯模型的SAR图像降斑方法

摘要

本发明公开了一种基于NSCT域局部高斯模型的SAR图像降斑方法，它涉及SAR图像处理技术领域，包括如下步骤：1.输入待降斑SAR图像；2.非下采样Contourlet变换；3.高频系数收缩；4.对低频子带系数和估计的无噪图像Contourlet高频各方向子带系数进行非下采样Contourlet逆变换，得到降斑后SAR图像；5.输出降斑后SAR图像。本发明对SAR图像同质区域抑斑效果优良，同时能够有效保持图像结构和纹理信息以及SAR图像的辐射特性。

主权项

一种基于NSCT域局部高斯模型的SAR图像降斑方法，包括如下步骤：(1)输入待降斑SAR图像；(2)非下采样Contourlet变换2a)按照下式，将SAR图像的乘性噪声转化为加性噪声：I＝RX＝X+(R‑1)X＝X+N其中，I为所观测到的被噪声污染的图像强度，R为相干斑，其均值为1，X为地物的真实后向散射强度，N为将要滤除的加性噪声；2b)进行非下采样Contourlet变换，获得Contourlet变换低频子带系数和高频方向子带系数；(3)高频系数收缩3a)保持低频子带系数不变；3b)利用降斑后图像与无噪图像之间的最小均方误差，按下式求解权重系数；ak＝[θk(y)Tθk(y)]‑1{θk(y)Ty‑σ2div[θk(y)]}其中，ak为权重系数，[θk(y)Tθk(y)]‑1为逆运算符号，θk(y)T为θk(y)的转置，θk(y)为阈值函数，k＝1，2……K，K的取值范围为1～100，y为输入的待降斑图像Contourlet变换高频方向子带系数，σ为噪声标准差，div为求导运算；3c)按照下式对高频方向子带系数进行收缩，得到估计的无噪图像Contourlet高频各方向子带系数： F ( y ) = Σ k = 1 K a k θ k ( y ) 其中，F(y)为估计的无噪图像Contourlet高频各方向子带系数，∑为求和运算，ak为权重系数，θk(y)为阈值函数，k＝1，2……K，K的取值范围为1～100；(4)对低频子带系数和估计的无噪图像Contourlet高频各方向子带系数进行非下采样Contourlet逆变换，得到降斑后SAR图像；(5)输出降斑后SAR图像。