[发明专利]基于第二代曲线波变换的SAR图像相干斑抑制方法

权利要求书

1.一种基于第二代曲线波变换的SAR图像相干斑抑制方法，包括如下步骤：

(1)对选取的测试图像I进行第二代曲线波变换，将其分成5层子带；

(2)保持第1层的系数不变，把第5层的系数全部置零；

(3)用EM方法分别估计第2～4层的混合高斯模型的参数向量θ；

(4)对第2～4层的系数进行标记，标记为1的系数保持不变，标记为0的系数置零，标记为2的系数用贝叶斯收缩因子进行收缩；

(5)对经过步骤(2)和步骤(4)处理后的系数进行第二代曲线波反变换得到重构图像R；

(6)检测重构图像R的边缘，并对其均匀区域进行均值滤波，得到滤波后图像U；

(7)将原始图像和滤波后的图像相减得到差值图像：V＝I-U，对该差值图像进行非线性各向异性扩散k次迭代到v^k；

(8)将v^k与滤波后图像U相加，得到相干斑抑制后的图像D＝U+v^k。

2.根据权利要求1所述的SAR图像相干斑抑制方法，其中步骤(3)所述的用EM方法分别估计第2～4层的混合高斯模型的参数向量θ，按如下步骤进行：

(2a)设Y为第S层的系数矩阵，S＝2、3、4，用两个零均值的混合高斯函数逼近Y的系数直方图，得到混合高斯函数为：

f(yn)=Σm=1m=2ϵmfm(yn)]]>

其中分别是f_m(y_n)的权值和方差，

将参数记为混合高斯模型的参数向量θ，即

(2b)选择参数的初始值为：ε₁＝ε₂＝1/2，σ22=2(max-min)/3;]]>

(2c)计算每个样本y_n高斯概率密度函数的后验概率：

E(hnm)=fm(yn)ϵmΣm=1m=2fm(yn)ϵm;]]>

(2d)用公式对参数向量θ进行迭代训练，N为Y中的样本个数；

(2e)当时，停止迭代，得到参数向量θ，表示第i次迭代得到的参数向量的第j个元素。