[发明专利]基于轮廓波变换的光学遥感图像变化检测方法

权利要求书

1.一种基于轮廓波变换的光学遥感图像变化检测方法，其特征是，包括以下步骤：

（1）对同一地域不同时间获取的两幅光学遥感图像进行滤波去噪，辐射校正和几何配准的预处理，得到预处理后的两幅图像X_A、X_B；

（2）利用预处理后的两幅图像构造比值法差异图X_l；

（3）利用预处理后的两幅图像构造均值比法差异图X_m；

（4）分别对比值法差异图X_l、均值比法差异图X_m进行N层轮廓波分解，得到每幅差异图在N个分解层上高频带和低频带上的多方向多尺度分解的轮廓波系数,N=5；

（5）对高频带和低频带的轮廓波系数用不同的融合算子进行融合处理，即对低频带轮廓波系数采用取平均值的方法进行融合，得到低频带融合轮廓波系数，对高频带轮廓波系数采用选邻域能量最小的方法进行融合处理，得到高频带融合轮廓波系数；

（6）由轮廓波生成融合后图像的轮廓波系数Y_f，Y_f由低频系数Y_f{1}和高频系数构成，其中k＝2,3,4,5,6，N＝1,2,3,4,5，t＝1,2,...,2^N，让k，N，t遍历这些取值，即可得到整体轮廓波系数Y_f；

（7）将融合后得到的轮廓波系数Y_f进行逆轮廓波变换，也称为图像重构，得到融合后差异图X_F；

（8）运用模糊局部C均值聚类方法对融合后差异图X_F进行图像分割，生成变化检测结果图，完成对两幅光学遥感图像变化信息的最终检测。

2.如权利要求1所述一种基于轮廓波变换的光学遥感图像变化检测方法，其特征是，其中步骤（3）所述的利用预处理后的两幅图像构造均值比法差异图X_m，构造方式为：

将图像X_A中位于i行j列的像素点(i,j)的灰度值X_A(i,j)和相对应的图像X_B中位于i行j列的像素点的灰度值X_B(i,j)，通过均值比值运算其中μ_A(i,j)，μ_B(i,j)分别为图像A，B的邻域像素平均值，得到均值比值法差异图X_m中位于i行j列的像素点(i,j)的灰度值X_m(i,j)，如果X_m(i,j)接近于0,表示图像X_A中位于i行j列的像素点没有随时间的变化而发生变化，反之，若X_m(i,j)远大于0，则认为该像素点发生了变化；对图像X_A和图像X_B中每个位于i行j列的像素点的灰度值从左到右，从上到下都进行此运算，构造均值比法差异图X_m。