[发明专利]基于图像融合的光学遥感图像变化检测

权利要求书

1.一种基于图像融合的光学遥感图像变化检测方法，包括如下步骤：

(1)对同一地域不同时间获取的两幅光学遥感图像进行滤波去噪，辐射校正和几何配准的预处理，得到预处理后的两幅图像X₁,X₂；

(2)利用预处理后的两幅图像X₁,X₂，分别构造差值法差异图X_D和比值法差异图X_R：

X_D＝|X₁-X₂|，

X_R＝X₁/X₂；

(3)对构造出的差值法差异图X_D和比值法差异图X_R分别进行N层小波分解，得到每幅差异图在N个分解层上高频带和低频带上的多尺度分解的小波系数，N＝3；

(4)对高频带和低频带的小波系数用不同的融合算子进行融合处理，即对低频带小波系数采用取平均值的方法进行融合，得到低频带融合小波系数，对高频带小波系数采用选能量最小的方法进行融合，得到高频带融合小波系数；

(5)对低频带融合小波系数和高频带融合小波系数进行小波逆变换，所得的重构图像即为融合后差异图X_F；

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

2.根据权利要求1所述的基于图像融合的光学遥感图像变化检测方法，其中步骤（3）中所述的对构造出的差值法差异图X_D和比值法差异图X_R分别进行N层小波分解，按如下步骤进行：

(3a)在第一层，即N＝1，将差值法差异图X_D分解为一个低频带子图像和三个高频带子图像得到一个低频带小波系数和三个高频带小波系数将比值法差异图X_R分解为一个低频带子图像和三个高频带子图像得到一个低频带小波系数和三个高频带小波系数其中LL表示低频带，ε＝LH,HL,HH，LH表示垂直方向的高频带，HL表示水平方向的高频带，HH表示对角线方向的高频带；

(3b)在第二层，即N＝2，将差值法差异图X_D经步骤(3a)分解所得的低频带子图像继续分解为一个低频带子图像和三个高频带子图像得到一个低频带小波系数和三个高频带小波系数将比值法差异图X_R经步骤(3a)分解所得的低频带子图像继续分解为一个低频带子图像和三个高频带子图像得到一个低频带小波系数和三个高频带小波系数

(3c)在第三层，即N＝3，将差值法差异图X_D经步骤(3b)分解所得的低频带子图像继续分解为一个低频带子图像和三个高频带子图像得到一个低频带小波系数和三个高频带小波系数将比值法差异图X_R经步骤(3b)分解所得的低频带子图像继续分解为一个低频带子图像和三个高频带子图像得到一个低频带小波系数和三个高频带小波系数