[发明专利]基于轮廓波变换的SAR图像变化检测差异图生成方法

权利要求书

1.一种基于轮廓波变换的SAR图像变化检测差异图生成方法，其特征在于：基于轮廓波变换的SAR图像变化检测差异图生成具体实现步骤包括有：

(1)输入两幅不同时间相同地域的经过滤波去噪，辐射校正与几何配准的预处理SAR图像X_A，X_B；

(2)对该两幅图像X_A，X_B构造对数比值图其中X_l为生成的对数比值差异图，其中X_l(i，j)，X_A(i，j)，X_B(i，j)分别为对数比值图X_l，图像X_A，图像X_B在图像中坐标为(i，j)的像素点的像素值；

(3)对该两幅图像X_A、X_B，构造均值比值图其中μ_A(i，j)，μ_B(i，j)分别为图像X_A，X_B中以坐标为(i，j)的像素点为中心邻域像素平均均值，X_m(i，j)为均值比值图X_m图像中坐标为(i，j)的像素点的像素值；

(4)分别对对数比值图像X_l、均值比值图像X_m进行Contourlet变换，得到对数比值图的Contourlet系数Y_l和均值比值图的Contourlet系数Y_m；

(5)对得到的对数比值图的Contourlet系数Y_l和均值比值图的Contourlet系数Y_m分为低频部分和高频部分进行融合：

5.1获取融合后图像低频系数，按取平均规则融合得到低频系数，即：其中和分别为Y_m和Y_l的低频部分；

5.2 获取融合后图像高频系数，按取邻域能量最小原则规则融合得到高频系数，即：YN,tf{k}(i,j)=YN,tl{k}(i,j),DN,tl{k}(i,j)≤DN,tm{k}(i,j)YN,tm{k}(i,j),DN,tl{k}(i,j)≥DN,tm{k}(i,j),]]>其中代表坐标位于(i，j)像素的contourlet系数M×N邻域内的能量，和分别为Y_m和Y_l的第N层第t个方向的的高频部分，k＝N+1；

(5)由Contourlet生成融合后图像的contourlet系数Y_f，Y_f由低频系数Y_f{1}和高频系数构成；

(6)将融合后得到的Contourlet系数Y_f进行逆Contourlet变换，也称为图像重构，得到目标图像X_F，即SAR图像X_A、X_B的差异信息图。