[发明专利]一种面向对象的多光谱高分辨率遥感影像变化检测方法

权利要求书

1.一种面向对象的基于CRF的多时相多光谱遥感影像变化检测方法，其特征在于，包括以下步骤：

步骤1：输入同一区域、不同时相的两幅高分辨率光学遥感影像，分别记为：X₁和X₂；

步骤2：利用遥感软件ENVI对X₁和X₂进行影像配准，分为粗校正和精校；

步骤3：利用多元变化检测MAD方法对X₁和X₂进行辐射归一化校正；

步骤4：对经过预处理的多时相高分辨率影像进行叠合，然后基于SLIC算法进行分割；

步骤5：计算多时相影像的差异影像，对所述差异影像的每一个波段，计算每一个对象的均值；记各波段各对象均值所组成的影像为X_M，继而利用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在Q₁₁＝(x₁,y₁),Q₁₂＝(x₁,y₂),Q₂₁＝(x₂,y₁)及Q₂₂＝(x₂,y₂)四个坐标点的值，x₁，x₂是横坐标，y₁，y₂是纵坐标，如果选择一个坐标系统使得这四个点的坐标分别为(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相对辐射归一化校正具体方法如下：

首先，通过典型相关分析找到两期影像各波段亮度值的一个线性组合，得到变化信息增强的差异影像，通过对归一化的MAD变量平方和取阈值确定变化和未变化区域；

然后，通过所述未变化区域对应的两时相像元对的映射方程，完成相对辐射校正。