基于核互模态因素分析核融合的多光谱图像变化检测方法技术

本发明专利技术公开了基于核互模态因素分析核融合的多光谱图像变化检测方法，主要解决现有差值核函数不能充分挖掘及利用图像多种特征之间的相关性而导致检测精度低下的问题。其实现过程是：1)将两时刻三通道光谱图像变换至HIS空间，得到两时刻图像的色调、饱和度、强度信息；2)根据1)的结果提取两时刻图像颜色特征和强度纹理特征；3)利用所提特征求取颜色差值核函数矩阵和强度纹理差值核函数矩阵；4)对3)获得的两个矩阵加权融合，构造合成核函数矩阵，并自适应选取最优的加权系数；5)将合成核函数矩阵输入支持向量机SVM进行检测，得到变化检测结果。本发明专利技术检测精度高，结果稳定，运算量低，可用于多光谱图像变化检测。

【技术实现步骤摘要】
基于核互模态因素分析核融合的多光谱图像变化检测方法
本专利技术属于图像处理
，特别涉及一种图像变化检测方法，可用于对多光谱图像进行地物状态变化监测与评估。
技术介绍
多光谱图像是遥感卫星获取的在多个波段上的对地观测数字图像。多光谱图像变化检测通过分析同一地区不同时刻观测到的多幅多光谱图像来确定地物变化的特征及过程。多光谱图像具有丰富的光谱信息，能够很好地反映地物特性。多光谱图像变化检测在军事和民事领域中具有突出的作用。最常用的多光谱图像变化检测方法是变化矢量分析法CVA，该方法将多通道的原始两时刻图像定义为向量，将两时刻的向量作差得到变化向量，再对变化向量进行分析从而检测出两时刻图像是否发生变化。CVA方法在多光谱图像变化检测中得到了成功的应用，近年各国学者对其进行了多方面的改进。JChen定义了后验概率空间中的CVA算法，该算法降低了传统CVA对几何配准和辐射校正的严格要求，不同变化所得到的结果都处于同一尺度空间，对单门限检测法的鲁棒性更好。FBovolo提出了极坐标空间中的CVA算法，该算法从理论上分析了极坐标空间中变化与非变化像素点的分布及其统计特性，给出了极坐标空间中变化矢量的严格定义及分析。此后，FBovolo又提出了2维CVA和多级CVA算法。但是上述的CVA方法及其改进方法，不能够实现对多光谱图像多种特征的自适应融合，不能充分利用多种特征中的信息，导致运算量比较大，计算效率比较低。除了CVA方法，其它方法也在多光谱图像变化检测中得到了成功的应用。例如，Hichri等在2013年提出了一种实现交式互分割的多光谱图像变化检测方法，该方法使用支持向量机SVM对已知标记的像素进行训练，使用训练结果指导图像分割得到初始变化检测图像，最后使用马尔科夫随机场MRFs调整初始变化检测结果图，得到最终的变化检测结果图。但是该方法实现时需要两个步骤，即先对已知标签样本进行训练得到初分结果，再进行优化，使得该算法流程比较复杂。近年来新发展起来的是基于核方法的多光谱图像变化检测算法。GustavoCamps-Valls在2008年首先提出将核方法应用于遥感图像变化检测，通过提取图像的强度信息和纹理信息并构造强度纹理差值核函数实现了遥感图像变化检测。但是该方法处理多光谱图像时，没有考虑多种特性之间的相关性，不能充分挖掘和融合多种图像特征，因而其检测精度有待进一步提高。
技术实现思路
本专利技术的目的是在针对上述已有算法的不足，提出一种基于核互模态因素分析核融合的多光谱图像变化检测方法，通过分析图像多种特征之间的相关性，充分挖掘和融合图像特征，从而降低运算量，提高运算效率，提高检测精度。为实现上述目的，本专利技术的实现步骤包括如下：(1)给定原始的预处理后的三通道的相同尺寸两时刻光谱图像，记为{Rt,Gt,Bt}，t＝t0,t1，其中，Rt,Gt,Bt分别表示原RGB空间中的三通道光谱图像，使用RGB-HIS变换将该三通道图像变换至HIS空间，得到HIS通道的两时刻光谱图像{Ht,It,St}，其中Ht为t时刻光谱图像的色调分量，St为t时刻光谱图像的饱和度分量，It为t时刻光谱图像的强度分量；(2)利用t0,t1两时刻光谱图像的色调饱和度和强度分量，提取光谱图像特征；2a)分别提取分量光谱图像的灰度值和分量光谱图像的灰度值，将其堆砌成两时刻光谱图像的颜色特征和2b)提取强度分量光谱图像的灰度特征和纹理特征将其堆砌成强度纹理特征(3)利用两时刻光谱图像的颜色特征求取颜色差值核函数矩阵KCol；(4)利用两时刻光谱图像的强度纹理特征求取强度纹理差值核函数矩阵KST；(5)利用基于核互模态因素分析算法KCFA分别求取颜色差值核函数矩阵KCol的权重系数wc和强度纹理差值核函数矩阵KST的权重系数ws；(6)利用步骤5)、步骤4)和步骤3)的结果，构造合成核函数矩阵KF：KF＝wc·KCol+ws·KST；(7)将步骤6)得到的合成核函数矩阵KF输入到支持向量机SVM中进行检测，得到二值变化检测结果，即完成对多光谱图像的变化检测。本专利技术具有如下优点：1)本专利技术由于利用了核方法，可以将实际多光谱图像处理中存在的非线性不可分问题转化为高维空间中的线性可分问题。2)本专利技术通过将原RGB空间的伪彩色图像变换至HIS空间，充分提取图像的颜色特征和强度纹理特征，有效突出了变化区域及其边界；并将颜色信息和强度纹理信息进行融合，最大限度挖掘两种特征之间的相关性，有效实现了核空间中的特征融合。3)本专利技术能全面地提取图像特征，可以显著提高变化检测结果的精度。仿真结果表明，本专利技术与现有的差值核方法相比，运算量低，检测精度高，检测性能更为优越。附图说明图1是本专利技术的实现流程图；图2是用本专利技术对日本地震灾害实测多光谱图像的变化检测结果图；图3是用本专利技术对格洛斯特洪水灾害实测多光谱图像的变化检测结果图；图4是用本专利技术对哈萨克斯坦海岸受流沙冲击前后实测多光谱图像的变化检测结果图；图5是用本专利技术对塔霍湖火灾实测多光谱图像的变化检测结果图。具体实施方式参照图1，本专利技术的具体实施过程如下：步骤1.将预处理后的光谱图像变换至HIS空间。给定原始的预处理后的三通道的相同尺寸两时刻光谱图像，记为{Rt,Gt,Bt}，t＝t0,t1，其中，Rt,Gt,Bt分别表示原RGB空间中的三通道光谱图像；使用RGB-HIS变换将该三通道光谱图像变换至HIS空间，得到HIS通道的两时刻光谱图像{Ht,It,St}，其中Ht为t时刻光谱图像的色调分量，St为t时刻光谱图像的饱和度分量，It为t时刻光谱图像的强度分量。其中，RGB空间中的R代表光谱图像的红色分量，G代表光谱图像的绿色分量，B代表光谱图像的蓝色分量；HIS空间中的H代表光谱图像的色调分量，S代表光谱图像的饱和度分量，I光谱图像的代表强度分量。步骤2.提取光谱图像的颜色特征和强度纹理特征。利用t0,t1两时刻光谱图像的色调饱和度和强度分量，提取光谱图像特征，其步骤如下：2a)提取t0时刻光谱图像分量的灰度值和分量的灰度值用这两个灰度值堆砌构成t0时刻光谱图像的颜色特征提取t1两时刻光谱图像分量的灰度值和分量的灰度值用这两个灰度值堆砌构成t1时刻光谱图像的颜色特征表示如下：2b)提取t0时刻光谱图像的强度分量的灰度值和纹理特征用该灰度值和纹理特征堆砌构成t0时刻光谱图像强度纹理特征提取t1时刻光谱图像分量灰度值和纹理特征用该灰度值和纹理特征堆砌构成t1时刻光谱图像的强度纹理特征表示如下：步骤3.构造颜色差值核函数矩阵KCol。利用t0时刻光谱图像的颜色特征和t1时刻光谱图像的颜色特征构造颜色差值核函数矩阵KCol，其步骤如下：3a)从t0时刻光谱图像中任意选取两个像素点和从t1时刻光谱图像中选取与和位置相同的两个像素点，记为和3b)从中找出像素点对应的颜色特征和像素点对应的颜色特征从中找出像素点对应的颜色特征和像素点对应的颜色特征3c)利用3b)获得的颜色特征求取颜色插值核函数的值：其中，K(·)表示高斯径向基RBF核函数，ζ代表核宽参数，exp为指数函数；3d)在t0时刻光谱图像中找出没有被选取过的任意两个像素以及中相同位置的像素点按照与3b)相同的方法找到与对应的颜色特征再按照与3c)相同的方法计算颜色差值核函本文档来自技高网...

【技术保护点】
一种基于核互模态因素分析核融合的多光谱图像变化检测方法，包括如下步骤：1)给定原始的预处理后的三通道的相同尺寸两时刻光谱图像，记为{R

【技术特征摘要】
1.一种基于核互模态因素分析核融合的多光谱图像变化检测方法，包括如下步骤：1)给定原始的预处理后的三通道的相同尺寸两时刻光谱图像，记为{Rt,Gt,Bt}，t＝t0,t1，其中，Rt,Gt,Bt分别表示原RGB空间中的三通道光谱图像，使用RGB-HIS变换将该三通道图像变换至HIS空间，得到HIS通道的两时刻光谱图像{Ht,It,St}，其中Ht为t时刻光谱图像的色调分量，St为t时刻光谱图像的饱和度分量，It为t时刻光谱图像的强度分量；2)利用t0,t1两时刻光谱图像的色调饱和度和强度分量，提取光谱图像特征；2a)分别提取分量光谱图像的灰度值和分量光谱图像的灰度值，将其堆砌成两时刻光谱图像的颜色特征和2b)提取强度分量光谱图像的灰度特征和纹理特征将其堆砌成强度纹理特征3)利用两时刻光谱图像的颜色特征求取颜色差值核函数矩阵KCol；4)利用两时刻光谱图像的强度纹理特征求取强度纹理差值核函数矩阵KST；5)利用基于核互模态因素分析算法KCFA分别求取颜色差值核函数矩阵KCol的权重系数wc和强度纹理差值核函数矩阵KST的权重系数ws；6)利用步骤5)、步骤4)和步骤3)的结果，构造合成核函数矩阵KF：KF＝wc·KCol+ws·KST；7)将步骤6)得到的合成核函数矩阵KF输入到支持向量机SVM中进行检测，得到二值变化检测结果，即完成对多光谱图像的变化检测。2.根据权利要求1所述的方法，其中步骤2a)中堆砌t0时刻颜色特征和t1时刻颜色特征表示如下：其中，和分别为提取出的t0时刻图像色调分量的灰度值和光谱图像饱和度分量的灰度值；和分别为提取出的t1时刻光谱图像色调分量的灰度值和光谱图像饱和度分量的灰度值。3.根据权利要求1所述的方法，其中步骤2b)中堆砌t0时刻强度纹理特征和t1时刻强度纹理特征表示如下：其中和分别表示t0时刻光谱图像强度分量的灰度值和t1时刻光谱图像强度分量的灰度值；和分别表示t0时刻光谱图像的强度分量的纹理特征和t1时刻光谱图像的强度分量的纹理特征。4.根据权利要求1所述方法，其中步骤3)中利用两时刻光谱图像的颜色特征构造颜色差值核函数矩阵KCol，按如下步骤进行：3a)从t0时刻光谱图像中任意选取两个像素点和从t1时刻光谱图像中选取与和位置相同的两个像素点，记为和3b)从中找出像素点对应的颜色特征和像素点对应的颜色特征从中找出像素点对应的颜色特征和像素点对应的颜色特征3c)利用从3b)中获得的颜色特征求取颜色插值核函...

【专利技术属性】
技术研发人员：李明，谭啸峰，张鹏，贾璐，吴艳，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西,61