多时相差异影像模值计算及变化检测方法技术

本发明专利技术公开了多时相差异影像模值计算及变化检测方法，首先计算多时相多光谱影像的差异影像，在此基础上，计算影像上每一个点的CST值，根据致信水平获取阈值，得到初步的变化检测结果，然后再对该初步结果进行众数滤波(嵌入空间信息)，并根据滤波结果重新计算非变化区域的均值和方差矩阵。重复上述过程直到检测结果没有变化为止。其中检测过程中置信水平的选择是通过伪训练样本集来选择的，在最优的置信水平基础上，获取最终的变化检测结果。本发明专利技术解决了多时相多光谱遥感影像背景信息复杂、噪声干扰严重的问题。

【技术实现步骤摘要】
多时相差异影像模值计算及变化检测方法
本专利技术属于遥感影像处理
，特别涉及了多时相差异影像模值计算及变化检测方法。
技术介绍
随着多时相遥感数据的不断积累以及空间数据库的相继建立，如何从这些遥感数据中提取和检测变化信息已成为遥感科学和地理信息科学的重要研究课题。根据同一区域不同时相的遥感影像，可以提取城市、环境等动态变化的信息，为资源管理与规划、环境保护等部门提供科学决策的依据。遥感影像的变化检测就是从不同时期的遥感数据中，定量地分析和确定地表变化的特征与过程。各国学者从不同的角度和应用研究提出了许多有效的检测算法，如变化矢量分析法(ChangeVectorAnalysis，CVA)、基于FuzzyC-means(FCM)的聚类方法等。其中，传统的基于卡方变换(Chi-SquaredTransform,CST)的多时相光学遥感变化检测，先计算差异影像的均值和方差矩阵，然后再基于置信水平，确定变化检测的阈值，进而得到变化检测结果。该类技术中，使用CST的不足是仅使用多时相高分辨率差异影像的光谱信息，没有利用空间信息。另外，在利用CST进行变化检测时，置信水平的选择非常关键，常规的方法是基于试错法来选择该参数。这种方法缺少推广性。针对上述问题，有必要在检测过程中加入空间约束来提高检测精度，另一方面，对关键参数的自适应选择也是检测技术中的重点。
技术实现思路
为了解决上述
技术介绍
提出的技术问题，本专利技术旨在提供基于空间约束的卡方变换的多时相遥感影像变化检测方法，采用空间约束的卡方变换结合自适应置信水平的变化检测方法，克服了多时相多光谱遥感影像背景信息复杂、噪声干扰严重的问题。为了实现上述技术目的，本专利技术的技术方案为：基于空间约束的卡方变换的多时相遥感影像变化检测方法，包括以下步骤：(1)输入两时相的高分辨率光学遥感影像，分别记为X1和X2；(2)对X1和X2进行影像配准；(3)利用多元变化检测方法分别对X1和X2进行辐射归一化校正；(4)计算多时相差异影像DX＝X1-X2；(5)计算DX的模值XM，利用贝叶斯原理，并基于最大期望算法获取最优分割阈值T，将|XM-T|≤δ的区域作为伪训练样本集，其中δ取XM动态范围的百分比；(6)设定置信水平的搜索范围和搜索步长；(7)确定DX中的非变化区域，计算非变化区域的均值矢量和方差矩阵，并计算DX上每个点的卡方值；初始化时，将整个DX作为非变化区域；(8)在给定的置信水平A的基础上，计算检测阈值，并根据该阈值进行卡方检测，获得初步的检测结果M0；(9)给定滤波窗口大小，对M0进行众数滤波，即滤波窗口内，如果变化的像素数目大于非变化的，则窗口中心点像素为变化的，反之，窗口中心点像素为非变化的，记众数滤波的结果为M1；(10)判断M1中的非变化区域相较于M0中的非变化区域是否有改变，如果两者没有改变，将给定的置信水平A加上步骤(6)设定的搜索步长作为新的置信水平，并返回步骤(7)；如果两者改变了，则将M1中的非变化区域作为新的非变化区域，返回步骤(7)；(11)当置信水平达到步骤(6)设定的搜索范围的上界，终止循环迭代；(12)针对每一个置信水平，计算伪训练样本集的精度，选择精度最高的伪训练样本集对应的置信水平，并在此基础上，输出最终的变化检测结果。进一步地，步骤(5)的具体过程如下：(a)假设XM影像上的非变化类ωn和变化类ωc服从如下的高斯分布：ωi∈{ωn,ωc}，mi∈{mn,mc}，σi∈{σn,σc}；上式中，mn和σn分别为非变化类的均值和方差，mc和σc分别为变化类的均值和方差；(b)采用最大期望算法算法估计mn、σn、mc和σc这四个参数；(c)依据贝叶斯最小误差准则，求解XM的分割阈值T，在高斯分布的情况下，等同于求解下式：上式中，为分割阈值T的估计值，p(ωn)、p(ωc)分别为非变化类先验概率和变化类先验概率；(d)根据估计的分割阈值，构建伪训练样本集：非变化区域伪训练样本集：变化区域伪训练样本集：其中，δ取XM动态范围的15％。进一步地，在步骤(6)中，设定置信水平的搜索范围为0.95-0.999，搜索步长为0.001。进一步地，在步骤(8)中，采用下式计算检测阈值：P(Cij<χ21-α(b))＝1-α上式中，1-α是置信水平，Cij表示DX在坐标(i，j)处的卡方值，χ21-α(b)为检测阈值。进一步地，在步骤(2)中，对X1和X2进行影像配准包括几何粗校正和几何精校正，所述几何粗校正的过程：(A)选择X1和X2分别作为基准影像和待校正影像；(B)在基准影像和待校正影像上分别采集地面控制点，地面控制点的数量大于等于9，且地面控制点均匀分布在影像上；(C)计算基准影像和待校正影像各地面控制点处的均方误差；(D)采用多项式纠正法对待校正影像进行纠正；(E)采用双线性插值法对待校正影像进行重采样；所述几何精校正是将经过几何粗校正的遥感影像，利用自动匹配与三角剖分法进行校正。采用上述技术方案带来的有益效果：(1)本专利技术将众数滤波引入到CST检测结果之后，使得使用CST变换具有空间约束能力；(2)本专利技术采用迭代的方法，估计非变化区域的均值和标准差，克服了仅仅利用整个多波段影像DX计算均值和方差矩阵的不足，可以使变化检测的结果更加可靠，也更加具有稳健性；(3)在本专利技术中，最优置信水平是基于伪训练样本集获取，使得该技术通用性较好。附图说明图1是本专利技术的基本流程图；图2(a)、2(b)分别是2007年1月的沙特阿拉伯Mina地区高分辨率IKONOS图像第3波段示意图、2007年12月的沙特阿拉伯的Mina地区高分辨率IKONOS图像第3波段示意图；图3是变化检测的参考图像；图4(a)、4(b)、4(c)、4(d)分别是EM-CVA算法、ICST算法、RCST算法、本专利技术算法的检测结果示意图。具体实施方式以下将结合附图，对本专利技术的技术方案进行详细说明。如图1所示，基于空间约束的卡方变换的多时相遥感影像变化检测方法，具体步骤如下：步骤1：输入同一区域、不同时相的两幅高分辨率光学遥感影像，分别记为X1和X2。步骤2：对X1和X2进行影像配准，分为粗校正和精校正两个步骤：对于几何粗校正，利用ENVI4.8软件中的相关功能实现，具体操作步骤为：(1)显示基准影像和待校正影像；(2)采集地面控制点GCPs，GCPs应均匀分布在整幅图像内，GCPs的数目至少大于等于9；(3)计算均方误差；(4)采用多项式纠正法对待校正影像进行纠正；(5)采用双线性插值进行重采样输出，若求未知函数f在点P＝(x,y)的值，假设已知函数f在Q11＝(x1,y1),Q12＝(x1,y2),Q21＝(x2,y1),及Q22＝(x2,y2)四个点的值，如果选择一个坐标系统使得这四个点的坐标分别为(0,0)、(0,1)、(1,0)和(1,1)，那么双线性插值公式就可以表示为：f(x,y)≈f(0,0)(1-x)(1-y)+f(1,0)x(1-y)+f(0,1)(1-x)y+f(1,1)xy对于几何精校正，将经过几何粗校正的多光谱遥感影像数据，利用自动匹配与三角剖分法进行几何精校正。采用逐点插入法构建Delaunay三角网，对每一个三角形，利用其三个顶点的行列号与其对应的基准影像同名本文档来自技高网...

【技术保护点】
1.计算多时相差异影像模值的方法，其特征在于，包括(a)假设XM影像上的非变化类ωn和变化类ωc服从如下的高斯分布：

【技术特征摘要】
1.计算多时相差异影像模值的方法，其特征在于，包括(a)假设XM影像上的非变化类ωn和变化类ωc服从如下的高斯分布：ωi∈{ωn,ωc},mi∈{mn,mc},σi∈{σn,σc}；上式中，mn和σn分别为非变化类的均值和方差，mc和σc分别为变化类的均值和方差；(b)采用最大期望算法算法估计mn、σn、mc和σc这四个参数；(c)依据贝叶斯最小误差准则，求解XM的分割阈...

【专利技术属性】
技术研发人员：陈哲，石爱业，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏,32