一种两时相遥感影像变化检测方法技术

本发明专利技术公开了一种两时相遥感影像变化检测方法，包括：获取两时相的原始光学遥感影像；对所述原始光学遥感影像进行影像配准；利用直方图调整方法对影像配准后的遥感影像进行辐射归一化校正；对辐射归一化校正后的遥感影像，计算每个相应波段的波段差；对每一个波段差，分别计算扩展属性剖面；将每一个波段差作为FCM的输入，并将扩展属性剖面作为空间邻域信息嵌入到FCM模型中，进而确定影像的变化区域和非变化区域。本发明专利技术可以提高变化检测的精度。

【技术实现步骤摘要】

本专利技术涉及遥感影像处理
，尤其涉及一种两时相遥感影像变化检测方法。
技术介绍
遥感影像的变化检测就是从不同时期的遥感数据中，定量地分析和确定地表变化的特征与过程。各国学者从不同的角度和应用研究提出了许多有效的检测算法，如变化矢量分析法(ChangeVectorAnalysis，CVA)、基于FuzzyC-means(FCM)的聚类方法等。其中，传统的基于FCM的多时相光学遥感变化检测，多先进行CVA变换，然后对变化矢量的幅值进行FCM聚类，进而得到变化检测结果。该类技术中，使用FCM的不足是仅仅使用变化矢量的幅值，使得原始多波段信息没有得到充分的挖掘。基于FCM的变化检测算法中，由于存在空间相关性，图像中的像素点所属类别(变化类与非变化类)不仅与当前像素点的灰度值有关，还受到其邻域内像素点灰度值的影响。目前在FCM算法的基础上，目前多以局部空间邻域均值或中值作为空间约束条件嵌入到FCM的目标函数中。然而这些方法使用固定的窗口得到空间信息，导致了尺寸比例的选择问题。另外，由于多时相光学遥感影像检测环境的复杂化以及目标先验信息匮乏等，导致这些算法都存在着一定的局限性。为此，有必要研究新的光学遥感图像变化检测技术来有效克服上述难点。
技术实现思路
本专利技术所要解决的技术问题在于，提供一种两时相遥感影像变化检测方法，基于扩展属性滤波(ExtendedAttributeProfile，EAP)和FCM(FuzzyC-Means)，提高了变换检测精度。为了解决上述技术问题，本专利技术提供了一种两时相遥感影像变化检测方法，包括：获取两时相的原始光学遥感影像；对所述原始光学遥感影像进行影像配准；利用直方图调整方法对影像配准后的遥感影像进行辐射归一化校正；对辐射归一化校正后的遥感影像，计算每个相应波段的波段差；对每一个波段差，分别计算扩展属性剖面；将每一个波段差作为FCM的输入，并将扩展属性剖面作为空间邻域信息嵌入到FCM模型中，确定影像的变化区域和非变化区域。实施本专利技术，具有如下有益效果：1)采用EAP提取影像的纹理特征，其滤波窗口尺寸的选择不是固定的，而是自适应变化的。2)在FCM目标函数中，利用EAP滤波影像作为空间约束，提高了影像的变化检测精度。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1是本专利技术提供的两时相遥感影像变化检测方法的一个实施例的流程示意图；图2是原始光学遥感影像；图3是本专利技术方法与其他方法的处理结果对比图。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。图1是本专利技术提供的两时相遥感影像变化检测方法的一个实施例的流程示意图，如图1所示，包括：S101、获取两时相的原始光学遥感影像X1和X2。其中，X1、X2是同一区域不同时相的两幅光学遥感影像。S102、对所述原始光学遥感影像X1和X2进行影像配准。具体的，步骤S102具体包括步骤：S1021、采用ENVI14.8遥感软件对光学遥感影像X1和X2进行几何粗校正。几何粗校正具体操作步骤为：(1)显示基准影像和待校正影像；(2)采集地面控制点GCPs；GCPs应均匀分布在整幅图像内，GCPs的数目至少大于等于9；(3)计算误差；(4)选择多项式模型；(5)采用双线性插值进行重采样输出。其中的双线性差值法为：若求未知函数f1在点P＝(x,y)的值，假设我们已知函数f1在Q11＝(x1,y1),Q12＝(x1,y2),Q21＝(x2,y1),及Q22＝(x2,y2)四个点的值。如果选择一个坐标系统使得这四个点的坐标分别为(0,0)、(0,1)、(1,0)和(1,1)，那么双线性插值公式就可以表示为：f1(x,y)≈f1(0,0)(1-x)(1-y)+f1(1,0)x(1-y)+f(0,1)(1-x)y+f1(1,1)xy。S1022、利用自动匹配与三角剖分法对几何粗校正后的遥感影像X1和X2进行几何精校正。其中，三角剖分法为，采用逐点插入法构建Delaunay三角网，对每一个三角形，利用其三个顶点的行列号与其对应的基准影像同名点的地理坐标来确定该三角形内部的仿射变换模型参数，对待校正影像进行纠正，得到校正后的遥感影像。S103利用直方图调整方法对影像配准后的遥感影像进行辐射归一化校正。具体的，步骤S103包括步骤：将影像配准后的遥感影像的像素灰度值分别调整为从而实现辐射归一化校正，其中，影像配准后的遥感影像各个波段服从高斯分布，μf、σf分别为参考影像的像素均值和标准差，gr、μr、σr分别为影像配准后的遥感影像的像素灰度值、均值和标准差。本步骤中按照公式分别对影像配准后的遥感影像X1和X2进行辐射归一化校正。S104、对辐射归一化校正后的遥感影像，计算每个相应波段的波段差。其中，计算到的波段差为Di，i＝1,2,…,B，B表示影像的波段数目。S105、对每一个波段差，分别计算扩展属性剖面。具体的，步骤S105包括步骤：S1051、对每一个波段差Di，设置其灰度级为L+1，其层级定义为Xl(Di)＝{p∈Di|Di(p)≥l本文档来自技高网...

【技术保护点】
一种两时相遥感影像变化检测方法，其特征在于，包括：获取两时相的原始光学遥感影像；对所述原始光学遥感影像进行影像配准；利用直方图调整方法对影像配准后的遥感影像进行辐射归一化校正；对辐射归一化校正后的遥感影像，计算每个相应波段的波段差；对每一个波段差，分别计算扩展属性剖面；将每一个波段差作为FCM的输入，并将扩展属性剖面作为空间邻域信息嵌入到FCM模型中，确定影像的变化区域和非变化区域。

【技术特征摘要】
1.一种两时相遥感影像变化检测方法，其特征在于，包括：
获取两时相的原始光学遥感影像；
对所述原始光学遥感影像进行影像配准；
利用直方图调整方法对影像配准后的遥感影像进行辐射归一化校正；
对辐射归一化校正后的遥感影像，计算每个相应波段的波段差；
对每一个波段差，分别计算扩展属性剖面；
将每一个波段差作为FCM的输入，并将扩展属性剖面作为空间邻域信息嵌
入到FCM模型中，确定影像的变化区域和非变化区域。
2.如权利要求1所述的两时相遥感影像变化检测方法，其特征在于，所述
对所述原始光学遥感影像进行影像配准，具体包括：
采用ENVI14.8遥感软件对原始光学遥感影像进行几何粗校正；
利用自动匹配与三角剖分法对几何粗...

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