专利名称：基于图像融合的光学遥感图像变化检测的制作方法
技术领域：
本发明属于图像处理领域，涉及光学遥感图像变化检测，可用于环境变化中的湖泊水位动态检测，农作物生长状态的动态检测，军事侦察，自然灾害评估等诸多领域中，提高光学遥感图像变化检测结果的精确度。
背景技术：
遥感图像的变化检测技术是遥感图像研究的重要组成部分，它是通过对比分析同一地域不同时刻获得的多时相遥感影像，根据图像之间的差异来得到人们所需要的地物或目标随时间发生的变化信息。变化检测技术可以检测出不同时期图像灰度值或局部纹理之间的变化，在此基础上获得感兴趣目标在形状、位置、数量及其它属性的变化情况。这些变 化可能是由图像场景的真实变化引起的，也可能是由入射角、大气条件、传感器精度、地面湿度等变化引起的。光学遥感图像的变化检测是检测同一场景的两幅或多幅光学遥感图像之间随时间变化发生的变化信息。变化检测过程主要包括三部分(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中每个位于i行j列的像素点的灰度值从左到右，从上到下都进行比值运算，构造如图4(b)所示的比值法差异图XK。步骤3，分别对差值法差异图Xd和比值法差异图Xk进行N层小波分解，得到每幅差异图在N个分解层上高频带和低频带上的多尺度分解的小波系数，N=3。(3.1)在第一层，即N= 1，将如图4(a)所示的差值法差异图Xd分解为一个低频带子图像Z1flU和三个高频带子图像Zf1,得到一个低频带小波系数Fi1和三个高频带小波系数;将如图4(b)所示的比值法差异图Xk分解为一个低频带子图像Xii和三个高频带子图像尤^，得到一个低频带小波系数好 =JP三个高频带小波系数Mf1，其中LL表示低频带，e=LH, HL, HH, LH表示垂直方向的高频带，HL表示水平方向的高频带，HH表示对角线方向的高频带；(3. 2)在第二层，即N=2，将如图4(a)所示的差值法差异图Xd经步骤(3a)分解所得的低频带子图像11」继续分解为一个低频带子图像Ifw和三个高频带子图像Ii32，得到一个低频带小波系数％:12和三个高频带小波系数;将如图4(b)所示的比值法差异图Xk经步骤(3a)分解所得的低频带子图像续分解为一个低频带子图像Xfi2和三个高频带子图像工i2，得到一个低频带小波系数Ml2和三个高频带小波系数灰(3. 3)在第三层，即N=3，将如图4(a)所示的差值法差异图Xd经步骤(3b)分解所得的低频带子图像4'2继续分解为一个低频带子图像Il3和三个高频带子图像If3，得到一个低频带小波系数》13和三个高频带小波系数,^ ;将如图4(b)所示的比值法差异图Xk经步骤(3b)分解所得的低频带子图像ZiI3继续分解为一个低频带子图像P和三个高频带子图像，得到一个低频带小波系数ML和三个高频带小波系数。步骤4，对高频带和低频带的小波系数用不同的融合算子进行融合处理，得到高频带、低频带融合小波系数。 (4. I)对低频带小波系数采用取平均值的方法进行融合，得到低频带融合小波系数通过如下公式进行wlu (/，,) = lV 丨”屮'丨、'丨、其中，D指代差值法差异图的子图像，R指代比值法差异图的子图像，F指代融合后差异图的子图像，(i，j)表示图像中像素点的位置为i行j列，LL表示低频带，力指融合后差异图像的低频带小波系数rfi在像素点(i，j)的值，指差值法差异图在第三分解层的低频带小波系数》1>3在像素点(i，j)的值，》t>3(U')指比值法差异图在第三分解层的低频带小波系数《1;在像素点(i，j)的值。由于分解后的低频带主要代表了源图像的轮廓特征，可以明确的反应两个不同时刻的源图像之间的变化信息，所以对低频带小波系数采用取平均值的方法进行融合，可以最大程度的增强变化区域的边缘特征。(4. 2)对高频带小波系数采用选能量最小的方法进行融合处理，得到高频带融合小波系数Z “ ；) — \nfN^-
[^￡.ArE￡JV(U)其中，N= 1，2，3，D指差值法差异图的子图像，R指比值法差异图的子图像，F指融合后差异图的子图像，(i，j)表示图像中像素点的位置为i行j列，4m为融合后差异图像中第N分解层的高频带小波系数<、.在像素点(i，j)的值，匕为融合后差异图像中第N分解层的高频带小波系数n 、在像素点(i，j)的值，rfvW 为融合后差异图像中第N分解层的高频带小波系数在像素点(i，j)的值，e = LH，HL，HH，LH表示垂直方向的高频带，HL表示水平方向的高频带，HH表示对角线方向的高频带；iC表示差值法差异图Xd第N分解层的局部能量系数，其中
'‘ ij
C⑷表示以(i，j)为中心的局部窗口#^中第k个小波系数值；
权利要求
1.一种基于图像融合的光学遥感图像变化检测方法，包括如下步骤 (1)对同一地域不同时间获取的两幅光学遥感图像进行滤波去噪，辐射校正和几何配准的预处理，得到预处理后的两幅图像X1, X2 ； (2)利用预处理后的两幅图像X1,X2，分别构造差值法差异图Xd和比值法差异图Xk Xd = I X1-X21 , Xr = ； (3)对构造出的差值法差异图Xd和比值法差异图Xr分别进行N层小波分解，得到每幅差异图在N个分解层上高频带和低频带上的多尺度分解的小波系数，N = 3 ; (4)对高频带和低频带的小波系数用不同的融合算子进行融合处理，即对低频带小波系数采用取平均值的方法进行融合，得到低频带融合小波系数，对高频带小波系数采用选能量最小的方法进行融合，得到高频带融合小波系数； (5)对低频带融合小波系数和高频带融合小波系数进行小波逆变换，所得的重构图像即为融合后差异图Xf; (6)运用模糊局部C均值聚类方法对融合后差异图Xf进行图像分割，生成变化检测结果图，完成对两幅光学遥感图像变化信息的最终检测。
2.根据权利要求I所述的基于图像融合的光学遥感图像变化检测方法，其中步骤(3)中所述的对构造出的差值法差异图Xd和比值法差异图Xk分别进行N层小波分解，按如下步骤进行 (3a)在第一层，即N= 1，将差值法差异图Xd分解为一个低频带子图像Ifii和三个高频带子图像,得到一个低频带小波系数和三个高频带小波系数汗^ ;将比值法差异图Xk分解为一个低频带子图像和三个高频带子图像，得到一个低频带小波系数Ml1和三个高频带小波系数，其中LL表示低频带，e = LH, HL, HH, LH表示垂直方向的高频带，HL表示水平方向的高频带，HH表示对角线方向的高频带； (3b)在第二层，即N = 2，将差值法差异图Xd经步骤(3a)分解所得的低频带子图像IfL继续分解为一个低频带子图像和三个高频带子图像，得到一个低频带小波系数『|2和三个高频带小波系数^^将比值法差异图Xk经步骤(3a)分解所得的低频带子图像11,继续分解为一个低频带子图像A和三个高频带子图像夂匕，得到一个低频带小波系数ML和三个高频带小波系数RR (3c)在第三层，即N = 3，将差值法差异图Xd经步骤(3b)分解所得的低频带子图像 继续分解为一个低频带子图像和三个高频带子图像，得到一个低频带小波系数灰左3和三个高频带小波系数将比值法差异图Xk经步骤(3b)分解所得的低频带子图像继续分解为一个低频带子图像Ii3和三个高频带子图像夂〗3，得到一个低频带小波系数灰左3和三个高频带小波系数砂=o
3.根据权利要求I所述的基于图像融合的光学遥感图像变化检测方法，其中步骤(4)中所述的对低频带小波系数采用取平均值的方法进行融合，通过如下公式进行
4.根据权利要求I所述的基于图像融合的光学遥感图像变化检测方法，其中步骤(4)中所述的对高频带小波系数采用选能量最小的方法进行融合，是选择能量值较小者的小波系数作为融合后的小波系数，通过如下公式进行
全文摘要
本发明公开一种基于图像融合的光学遥感图像变化检测方法，主要解决现有变化检测技术检测结果精度不高的问题。其实现过程为对同一地域不同时间获取的两幅光学遥感图像预处理后，分别构造差值法差异图和比值法差异图；对这两幅差异图分别进行N层小波分解后获取各分解层高频带和低频带的小波系数；对高频带和低频带的小波系数分别用不同的融合算子融合处理，得到高频带和低频带的融合小波系数；对高频带和低频带的融合小波系数进行逆变换得到融合后差异图；运用模糊局部C均值聚类法分割融合后的差异图，得到变化检测结果。本发明通过提高光学遥感图像的差异图性能，提高了检测结果的精度，可用于自然灾害评估和环境检测。
文档编号G06T7/00GK102750705SQ20121023407
公开日2012年10月24日 申请日期2012年7月8日 优先权日2012年7月8日
发明者公茂果, 李瑜, 焦李成, 王桂婷, 王爽, 翟路, 贾萌, 陈默, 马文萍, 马晶晶 申请人:西安电子科技大学