一种面向对象的基于CRF的多时相多光谱遥感影像变化检测方法技术

本发明专利技术公开了一种面向对象的基于CRF的多时相多光谱遥感影像变化检测方法，包括以下步骤：输入同一区域、不同时相的两幅高分辨率光学遥感影像；利用遥感软件ENVI对X1和X2进行影像配准，分为粗校正和精校；利用多元变化检测MAD方法对X1和X2进行辐射归一化校正；基于SLIC算法进行分割；计算多时相影像的差异影像，利用FCM算法进行聚类；根据SLIC的分割图像构建CRF的多邻域系统，构建CRF的二元项能量；根据循环信度传播LBP优化算法对CRF能量项U进行优化，获取最终的变化检测结果。本发明专利技术使得变化检测的结果更加可靠，也更加具有稳健性。

An Object-Oriented Multi-temporal and Multi-spectral Remote Sensing Image Change Detection Method Based on CRF

The invention discloses an object-oriented multi-temporal and multi-spectral remote sensing image change detection method based on CRF, which includes the following steps: input two high-resolution optical remote sensing images of the same region and different time phases; use remote sensing software ENVI to register X1 and X2 images, and divide them into rough correction and correction; use multi-variate change detection MAD method to carry out radiation normalization correction for X1 and X2. Segmentation based on SLIC algorithm; Computing difference images of multi-temporal images, clustering using FCM algorithm; Constructing multi-neighborhood system of CRF based on SLIC segmentation images, building binary energy of CRF; Optimizing CFR energy item U based on cyclic reliability propagation LBP optimization algorithm to obtain final change detection results. The invention makes the result of change detection more reliable and more robust.

【技术实现步骤摘要】
一种面向对象的基于CRF的多时相多光谱遥感影像变化检测方法
本专利技术涉及一种面向对象的基于条件随机场CRF的多时相多光谱高分辨率遥感影像非监督变化检测方法，属于光学遥感影像变化检测

技术介绍
随着多时相遥感数据的不断积累以及空间数据库的相继建立，如何从这些遥感数据中提取和检测变化信息已成为遥感科学和地理信息科学的重要研究课题。根据同一区域不同时相的遥感影像，可以提取城市、环境等动态变化的信息，为资源管理与规划、环境保护等部门提供科学决策的依据。遥感影像的变化检测就是从不同时期的遥感数据中，定量地分析和确定地表变化的特征与过程。各国学者从不同的角度和应用研究提出了许多有效的检测算法，总体来说，根据检测过程中是否需要训练样本，可以将变化检测分为三个大类：非监督变化检测算法、半监督变化检测算法和监督变化检测算法。非监督变化检测一般先通过构造差分影像，进而再通过预先设定的阈值将差分影像分为变化和非变化区域。传统的非监督变化检测主要有算法运算法、CVA法及其改进的算法、主成分变换(PrincipleComponentAnalysis，PCA)法、影像回归法等。上述传统的变化检测方法大多是基于像素的方法，适用于中-低分辨率的遥感影像，且这些方法都是暗含像素在空间上具有独立性。针对高分辨率影像，由于其结构、纹理等信息更加突出，面向对象的变化检测更有利于知识的结合利用，也更能有效地利用高分辨率影像所具有的多特征优势。另外一方面，变化检测建模过程中，像素和其邻域具有高度相关性，MRF和CRF两种概率图模型均可以刻画这种邻域关系。由于CRF能够同时考虑观察场与标记场的空间上下文信息，比MRF建模更加灵活，其在多时相变化检测领域应用具有优势。如李民等提出的基于非平稳分析与条件随机场的SAR图像变化检测方法(2015，中国专利技术专利，西安电子科技大学)，但是该专利技术专利技术是监督类变化检测，实际应用中需消耗大量的人力成本等去构建训练样本。Cao等提出的基于CRF的非监督变换检测方法应用于多光谱多时相遥感图像检测过程中，提高了检测的精度。但是该方法的是基于像素的，检测的效率和精度有一定的局限性。B.Desclée等提出的面向对象的基于卡方变换的检测方法(OBCT)，可以克服基于像素的变化检测的局限性，但是该面向对象的方法没有考虑对象之间的相关性。
技术实现思路
针对现有技术存在的不足，本专利技术目的是提供一种面向对象的基于条件随机场CRF的多时相多光谱高分辨率遥感影像非监督变化检测方法，本专利技术首先利用简单线性迭代聚类SLIC算法对多时相遥感影像进行分割，然后对差异影像的均值应用FCM算法算法获取CRF的一元能量项；其次根据图像的每一个节点(影像的对象)的邻域和差异影像的幅值，获取CRF的二元能量项；最后通过循环信度传播LBP算法对CRF的能量最小化获取最终的变化检测结果，使得变化检测的结果更加可靠，也更加具有稳健性。为了实现上述目的，本专利技术是通过如下的技术方案来实现：本专利技术的一种面向对象的基于CRF的多时相多光谱遥感影像变化检测方法，包括以下步骤：步骤1：输入同一区域、不同时相的两幅高分辨率光学遥感影像，分别记为：X1和X2；步骤2：利用遥感软件ENVI对X1和X2进行影像配准，分为粗校正和精校；步骤3：利用多元变化检测MAD方法对X1和X2进行辐射归一化校正；步骤4：对经过预处理的多时相高分辨率影像进行叠合，然后基于SLIC算法进行分割；步骤5：计算多时相影像的差异影像，对所述差异影像的每一个波段，计算每一个对象的均值；记各波段各对象均值所组成的影像为XM，继而利用FCM算法进行聚类；步骤6：根据SLIC的分割图像构建CRF的多邻域系统，并结合步骤5的差异影像，构建CRF的二元项能量；步骤7：根据循环信度传播LBP优化算法对CRF能量项U进行优化，获取最终的变化检测结果。步骤2中，对于所述粗校正，具体步骤如下：2.1)显示基准影像和待校正影像；2.2)采集地面控制点GCPs，所述GCPs均匀分布在整幅图像内，所述GCPs的数目至少大于等于25；2.3)计算均方根误差；2.4)选择多项式校正模型；2.5)采用双线性插值进行重采样输出。上述双线性插值具体方法如下：求未知函数f在坐标点P＝(x,y)的值，x是横坐标，y是纵坐标；基于预设的已知函数f在Q11＝(x1,y1),Q12＝(x1,y2),Q21＝(x2,y1)及Q22＝(x2,y2)四个坐标点的值(x1，x2是横坐标，y1，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(1)步骤2中，所述精校正即将经过几何粗校正的多光谱遥感影像数据，利用自动匹配与三角剖分法进行几何精校正。上述三角剖分法即采用逐点插入法构建Delaunay三角网，对每一个三角形，利用其三个顶点的行列号与其对应的基准影像同名点的地理坐标来确定该三角形内部的仿射变换模型参数，对待校正影像进行纠正，得到校正后的遥感影。步骤3中，所述多元变化检测MAD相对辐射归一化校正具体方法如下：首先，通过典型相关分析找到两期影像各波段亮度值的一个线性组合，得到变化信息增强的差异影像，通过对归一化的MAD变量平方和取阈值确定变化和未变化区域；然后，通过所述未变化区域对应的两时相像元对的映射方程，完成相对辐射校正。步骤4，具体实现步骤如下：4.1)初始化种子点即聚类中心：按照设定的超像素个数，在图像内均匀的分配种子点；基于预设的图片总共有N个像素点，预分割为K个相同尺寸的超像素，那么每个超像素的大小为N/K，则相邻种子点的距离即步长为4.2)在种子点的z×z邻域内重新选择种子点；4.3)在每个种子点周围的邻域内为每个像素点分配类标签，搜索范围限制为2S×2S；4.4)对于每个搜索到的像素点，分别计算它和种子点的颜色距离和空间距离；取最小值对应的种子点作为像素点的聚类中心；4.5)对步骤4.1)-步骤4.4)进行迭代直到误差收敛。步骤4.2)中，所述种子点具体选择方法为：计算邻域内所有像素点的梯度值，将种子点移到邻域内梯度最小的地方。步骤5中，基于预设的加权指数m＝2，聚类数目C＝2，具体方法如下：5.1)基于预设的构建FCM：式中，C是聚类数目，K是对象样本的总数，k表示样本的序号，j表示类别中心的序号，XM表示多时相差异影像各波段所有对象均值所组成的影像，vj表示第j聚类中心，dkj表示第k个样本和第j聚类中心的欧式距离，表示第k样本对于第j类聚类中心的模糊隶属度，m为隶属度的加权指数，ujk∈[0,1]且5.2)公式(2)的目标函数最小化用下述的公式交替进行：式中，dki表示第k个样本和第i聚类中心的欧式距离、dkj表示第k个样本和第j聚类中心的欧式距离；5.3)FCM算法收敛后，获得CRF的一元能量项Uj＝ujk。步骤6)中，构建CRF的二元项能量如下：式中，i,j表示图像相邻的两个节点即相邻的对象(本专利技术中，将每一个对象看作图模型中的一个节点)；ci,cj表示两个节点i,j处的类别标签；两种类别为变化和非变化本文档来自技高网...

【技术保护点】
1.一种面向对象的基于CRF的多时相多光谱遥感影像变化检测方法，其特征在于，包括以下步骤：步骤1：输入同一区域、不同时相的两幅高分辨率光学遥感影像，分别记为：X1和X2；步骤2：利用遥感软件ENVI对X1和X2进行影像配准，分为粗校正和精校；步骤3：利用多元变化检测MAD方法对X1和X2进行辐射归一化校正；步骤4：对经过预处理的多时相高分辨率影像进行叠合，然后基于SLIC算法进行分割；步骤5：计算多时相影像的差异影像，对所述差异影像的每一个波段，计算每一个对象的均值；记各波段各对象均值所组成的影像为XM，继而利用FCM算法进行聚类；步骤6：根据SLIC的分割图像构建CRF的多邻域系统，并结合步骤5的差异影像，构建CRF的二元项能量；步骤7：根据循环信度传播LBP优化算法对CRF能量项U进行优化，获取最终的变化检测结果。

【技术特征摘要】
1.一种面向对象的基于CRF的多时相多光谱遥感影像变化检测方法，其特征在于，包括以下步骤：步骤1：输入同一区域、不同时相的两幅高分辨率光学遥感影像，分别记为：X1和X2；步骤2：利用遥感软件ENVI对X1和X2进行影像配准，分为粗校正和精校；步骤3：利用多元变化检测MAD方法对X1和X2进行辐射归一化校正；步骤4：对经过预处理的多时相高分辨率影像进行叠合，然后基于SLIC算法进行分割；步骤5：计算多时相影像的差异影像，对所述差异影像的每一个波段，计算每一个对象的均值；记各波段各对象均值所组成的影像为XM，继而利用FCM算法进行聚类；步骤6：根据SLIC的分割图像构建CRF的多邻域系统，并结合步骤5的差异影像，构建CRF的二元项能量；步骤7：根据循环信度传播LBP优化算法对CRF能量项U进行优化，获取最终的变化检测结果。2.根据权利要求1面向对象的基于CRF的多时相多光谱遥感影像变化检测方法，其特征在于，步骤2中，对于所述粗校正，具体步骤如下：2.1)显示基准影像和待校正影像；2.2)采集地面控制点GCPs，所述GCPs均匀分布在整幅图像内，所述GCPs的数目至少大于等于25；2.3)计算均方根误差；2.4)选择多项式校正模型；2.5)采用双线性插值进行重采样输出。3.根据权利要求2面向对象的基于CRF的多时相多光谱遥感影像变化检测方法，其特征在于，所述双线性插值具体方法如下：求未知函数f在坐标点P＝(x,y)的值，x和y分别表示横坐标和纵坐标；基于预设的已知函数f在Q11＝(x1,y1),Q12＝(x1,y2),Q21＝(x2,y1)及Q22＝(x2,y2)四个坐标点的值，x1，x2是横坐标，y1，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(1)。4.根据权利要求1面向对象的基于CRF的多时相多光谱遥感影像变化检测方法，其特征在于，步骤2中，所述精校正即将经过几何粗校正的多光谱遥感影像数据，利用自动匹配与三角剖分法进行几何精校正。5.根据权利要求4面向对象的基于CRF的多时相多光谱遥感影像变化检测方法，其特征在于，所述三角剖分法即采用逐点插入法构建Delaunay三角网，对每一个三角形，利用其三个顶点的行列号与其对应的基准影像同名点的地理坐标来确定该三角形内部的仿射变换模型参数，对待校正影像进行纠正，得到校正后的遥感影。6.根据权利要求1面向对象的基于CRF的多时相多光谱遥感影像变化检测方法，其特征在于，步骤3中，所述多元变化检测MAD相对辐射归一化校正具体方法如下：...

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