基于图像融合的光学遥感图像变化检测制造技术

本发明专利技术公开一种基于图像融合的光学遥感图像变化检测方法，主要解决现有变化检测技术检测结果精度不高的问题。其实现过程为：对同一地域不同时间获取的两幅光学遥感图像预处理后，分别构造差值法差异图和比值法差异图；对这两幅差异图分别进行N层小波分解后获取各分解层高频带和低频带的小波系数；对高频带和低频带的小波系数分别用不同的融合算子融合处理，得到高频带和低频带的融合小波系数；对高频带和低频带的融合小波系数进行逆变换得到融合后差异图；运用模糊局部C均值聚类法分割融合后的差异图，得到变化检测结果。本发明专利技术通过提高光学遥感图像的差异图性能，提高了检测结果的精度，可用于自然灾害评估和环境检测。

【技术实现步骤摘要】

本专利技术属于图像处理领域，涉及光学遥感图像变化检测，可用于环境变化中的湖泊水位动态检测，农作物生长状态的动态检测，军事侦察，自然灾害评估等诸多领域中，提高光学遥感图像变化检测结果的精确度。
技术介绍
遥感图像的变化检测技术是遥感图像研究的重要组成部分，它是通过对比分析同一地域不同时刻获得的多时相遥感影像，根据图像之间的差异来得到人们所需要的地物或目标随时间发生的变化信息。变化检测技术可以检测出不同时期图像灰度值或局部纹理之间的变化，在此基础上获得感兴趣目标在形状、位置、数量及其它属性的变化情况。这些变 化可能是由图像场景的真实变化引起的，也可能是由入射角、大气条件、传感器精度、地面湿度等变化引起的。光学遥感图像的变化检测是检测同一场景的两幅或多幅光学遥感图像之间随时间变化发生的变化信息。变化检测过程主要包括三部分(I)图像预处理(2)构造差异图(3)对差异图进行分析。其中图像的预处理主要包括几何配准和辐射校正，消除图像的几何误差，以此达到同一区域不同图像的地理坐标的匹配，消除传感器自身引起的，大气辐射引起的辐射噪声。构造差异图通常采用图像差值法和比值法构造，这两种方法构造的差异图进行变化检测操作简单，对差异图进行分析是利用阈值或聚类的分割方法生成最终的变化检测结果图。差值法是对多时相图像中对应的像素点灰度值进行相减，结果代表两个图像随时间发生的变化情况。如果差值为0，则判定该像素点的信息未发生变化，否则，判定像素点的信息发生了变化。比值法是计算已配准的多时相图像对应像素点灰度值的比值，结果代表两个图像随时间发生的变化情况。如果一个像素点没有发生变化，则比值应该接近于1，反之，则远大于I或远小于I。尽管用差值法构造差异图和比值法构造差异图这两种变化检测方法都可以有效的检测出光学遥感图像的变化信息，但这两种方法仍然存在各自的不足差值法容易受图像成像质量、波谱特征等客观条件影响，易产生“伪变化”信息；不能完全反映地面地物的辐射能量变化，如灰度值从50到30与250到220的变化，地面相应处的能量变化并不等量，不能等同看待；另外利用差值法构造的差异图对光学遥感图像进行变化检测时，检测结果漏检率较高。图像比值法在增强变化信息、抑制背景信息的同时可以帮助减少大气条件造成的影响，但有时会过于夸大部分的变化，如对灰度值从200到20与从20到10的变化，比值法不能区分出来，而用差值法得到的差异却非常大；此外，用比值法构造的差异影像进行光学遥感图像变化检测时，检测结果具有较高的误检率。
技术实现思路
本专利技术的目的在于针对上述差值法构造差异图和比值法构造差异图这两种变化检测方法各自的不足，提出了一种基于图像融合的光学遥感图像变化检测方法，以全面反映地面地物的辐射能量变化，减少“伪变化”信息，增强变化信息、抑制背景信息，提高检测正确率和检测精度。实现本专利技术目的的技术思路是采用基于离散小波变换的图像融合方法，对用差值法构造的差异图和用比值法构造的差异图进行融合，实现对光学遥感图像的变化检测。具体步骤包括如下(I)对同一地域不同时间获取的两幅光学遥感图像进行滤波去噪，辐射校正和几何配准的预处理，得到预处理后的两幅图像X1, X2 ；(2)利用预处理后的两幅图像X1, X2,分别构造差值法差异图Xd和比值法差异图XeXd = I X「X21， Xe = X1ZX2 ； (3)对构造出的差值法差异图Xd和比值法差异图Xk分别进行N层小波分解，得到每幅差异图在N个分解层上高频带和低频带上的多尺度分解的小波系数，N = 3 ;(4)对高频带和低频带的小波系数用不同的融合算子进行融合处理，即对低频带小波系数采用取平均值的方法进行融合，得到低频带融合小波系数，对高频带小波系数采用选能量最小的方法进行融合处理，得到高频带融合小波系数；(5)对低频带融合小波系数和高频带融合小波系数进行小波逆变换，所得的重构图像即为融合后差异图Xf;(6)运用模糊局部C均值聚类方法对融合后差异图Xf进行图像分割，生成变化检测结果图，完成对两幅光学遥感图像变化信息的最终检测。本专利技术与现有的技术相比具有如下优点I、本专利技术由于采用了基于小波分解的图像融合技术，在时间和空间上分解频率，可以更容易的提取细节信息，可以很好的保存图像的细节信息，从而得到正确率更高的变化检测结果。2、本专利技术由于利用小波分解的方向性，针对人眼对不同方向的高频分量具有不同的分辨率这一视觉特性，可以获得视觉效果更佳的融合图像。3、本专利技术由于选择不同的合适的融合规则融合图像，可以更大程度上的抑制变化区域的背景信息，增强变化区域的变化信息，使差值法的漏检率和比值法的误检率互补，获得了较高的检测精度。附图说明图I是本专利技术的实现总流程图；图2是本专利技术中基于小波变换的图像融合子流程图；图3是本专利技术仿真使用的Mexico地区两幅光学遥感图像和一幅标准参考图；图4是本专利技术分别使用差值运算、比值运算和图像融合法构造的Mexico地区三幅光学遥感差异图；图5是用本专利技术和现有的差值法、比值法对Mexico地区光学遥感图像进行变化检测的仿真结果图； 图6是本专利技术仿真使用的意大利Sardinia地区两幅光学遥感图像和一幅标准参考图；图7是本专利技术分别使用差值运算、比值运算和图像融合法构造的意大利Sardinia地区三幅光学遥感图像差异图；图8是用本专利技术和现有的差值法、比值法对意大利Sardinia地区光学遥感图像进行变化检测的仿真结果图。具体实施例方式以下参照附图，对本专利技术的实现方案及优点进行详细描述。实施例1，参照图I，本实例的实现步骤如下 步骤I,在Mexico同一地域的不同时间获取两幅光学遥感图像,将这两幅光学遥感图像进行滤波去噪，辐射校正与几何配准的预处理，得到处理后的两幅图像X1, X2，其中预处理后得到的图像X1如图3(a)所示，预处理后得到的图像X2如图3(b)所示。通过预处理可以消除图像的几何误差，以达到对同一区域不同图像的地理坐标的匹配，消除传感器自身引起的噪声和大气辐射引起的辐射噪声。步骤2，用预处理后的两幅图像X1, X2，分别构造差值法差异图Xd和比值法差异图XK。(2. I)将如图3 (a)所示的图像X1中位于i行j列的像素点(i，j)的灰度值X1 (i, j)和相对应的图3(b)所示的图像X2中位于i行j列的像素点的灰度值X2 (i，j)，通过差值运算\(i，j) = IX1(IjJ)-X2(IjJ) I得到差值法差异图Xd中位于i行j列的像素点(i，j)的灰度值XD(i，j)，若XD(i，j)为0，则表示图像X1中该像素点没有随时间的变化而发生变化，否贝U，认为图像X1中的该像素点发生了变化；对图像X1和图像X2中每个位于i行j列的像素点的灰度值从左到右，从上到下都进行差值运算，得到如图4(a)所示的差值法差异图Xd ；(2. 2)将如图3 (a)所示的图像X1中位于i行j列的像素点(i，j)的灰度值X1 (i, j)和相对应的图3(b)所示的图像X2中位于i行j列的像素点的灰度值X2(i，j)，通过比值运算XK(i，j) =X1Q, j)/X2(i，j)，得到比值法差异图Xk中位于i行j列的像素点(i，j)的灰度值XK(i，j)，如果比值接近于1，表示图像X1中位于i行j列的像素点本文档来自技高网...

【技术保护点】
一种基于图像融合的光学遥感图像变化检测方法，包括如下步骤：(1)对同一地域不同时间获取的两幅光学遥感图像进行滤波去噪，辐射校正和几何配准的预处理，得到预处理后的两幅图像X1,X2；(2)利用预处理后的两幅图像X1,X2，分别构造差值法差异图XD和比值法差异图XR：XD＝|X1？X2|，XR＝X1/X2；(3)对构造出的差值法差异图XD和比值法差异图XR分别进行N层小波分解，得到每幅差异图在N个分解层上高频带和低频带上的多尺度分解的小波系数，N＝3；(4)对高频带和低频带的小波系数用不同的融合算子进行融合处理，即对低频带小波系数采用取平均值的方法进行融合，得到低频带融合小波系数，对高频带小波系数采用选能量最小的方法进行融合，得到高频带融合小波系数；(5)对低频带融合小波系数和高频带融合小波系数进行小波逆变换，所得的重构图像即为融合后差异图XF；(6)运用模糊局部C均值聚类方法对融合后差异图XF进行图像分割，生成变化检测结果图，完成对两幅光学遥感图像变化信息的最终检测。

【技术特征摘要】

【专利技术属性】
技术研发人员：公茂果，焦李成，翟路，马晶晶，贾萌，李瑜，陈默，王爽，王桂婷，马文萍，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：