首页 / 专利库 / 图形技术 / 纹理合成 / 基于核互模态因素分析核融合的多光谱图像变化检测方法

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

阅读:553发布:2020-05-12

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

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)中获得的颜色特征 求取颜色插值核函数
的值:
其中,K(·)表示高斯径向基RBF核函数,ζ代表核宽参数,exp为指数函数;
3d)在t0时刻光谱图像 中找出没有被选取过的任意两个像素 以
及 中相同位置的像素点 按照与3b)相同的方法找到与
对应的颜色特征 再按照与3c)相同的方法计算颜色
差值核函数 的值,其中p是像素点的标号,1≤p≤N-1,N是光谱图
像像素点个数;
3e)重复步骤3d),直到所有像素点被选取过,得到颜色差值核函数的全部值,构成颜色差值核函数矩阵KCol。
5.根据权利要求1所述方法,其中步骤4)中利用两时刻光谱图像强度纹理特征构造强度纹理差值核函数矩阵KST,按如下步骤进行:
4a)从t0时刻光谱图像 中任意选取两个像素点 和 从t1时刻光谱图
像 中选取与 和 位置相同的两个像素点,记为 和
4b)从 中找出像素点 对应的强度纹理特征 和像素点 对应的强度纹理特征从 中找出像素点 对应的强度纹理特征 和像素点 对应的强度纹理特征
4c)利用从4b) 中获得的 求取强度纹理差值核函数
的值:
其中,K(·)表示高斯径向基RBF核函数,ζ代表核宽参数,exp为指数函数;
4d)在t0时刻光谱图像 中找出没有被选取过的任意两个像素 以
及 中相同位置的像素点 按照与4b)相同的方法找到与
对应的强度纹理特征 再按照与4c)相同的方法计算强
度纹理差值核函数 的值,其中q是像素点的标号,1≤q≤N-1,N是光
谱图像像素点个数;
4e)重复步骤4d),直到所有像素点被选取过,得到强度纹理差值核函数的全部值,构成强度纹理差值核函数矩阵KST。
6.根据权利要求1所述方法,其中步骤5)中求取颜色差值核函数矩阵KCol的权重系数wc和强度纹理差值核函数矩阵KST的权重系数ws,按照如下步骤进行:
5a)通过特征值分解法得到由KST*KCol特征向量构成的第一投影矩阵V和特征值构成的第一对矩阵λ,KST*KCol、V、λ这三者的关系如下:
(KST*KCol)*V=λ*V;
5b)提取第一投影矩阵V的第一列,记为α,定义颜色差值核函数矩阵KCol融合系数wc:
其中,T表示矩阵转置,KCol是颜色差值核函数矩阵;
5c)通过特征值分解法得到矩阵KCol*KST特征向量构成的第二投影矩阵U和特征值构成的第二对角矩阵λ′,KCol*KST、U、λ′这三者的关系如下:
(KCol*KST)*U=λ′*U;
5d)提取第二投影矩阵U的第一列,记为β;定义强度纹理插值核函数矩阵融合系数wc:
其中,KST是强度纹理差值核函数矩阵。
7.根据权利要求1所述方法,其中步骤6)中利用步骤5)、步骤4)以及步骤3)的结果构造合成核函数矩阵KF,按照如下公式进行:
KF=wc·KCol+ws·KST。

说明书全文

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

技术领域

[0001] 本发明属于图像处理技术领域,特别涉及一种图像变化检测方法,可用于对多光谱图像进行地物状态变化监测与评估。

背景技术

[0002] 多光谱图像是遥感卫星获取的在多个波段上的对地观测数字图像。多光谱图像变化检测通过分析同一地区不同时刻观测到的多幅多光谱图像来确定地物变化的特征及过程。多光谱图像具有丰富的光谱信息,能够很好地反映地物特性。多光谱图像变化检测在军事和民事领域中具有突出的作用。
[0003] 最常用的多光谱图像变化检测方法是变化矢量分析法CVA,该方法将多通道的原始两时刻图像定义为向量,将两时刻的向量作差得到变化向量,再对变化向量进行分析从而检测出两时刻图像是否发生变化。CVA方法在多光谱图像变化检测中得到了成功的应用,近年各国学者对其进行了多方面的改进。
[0004] J Chen定义了后验概率空间中的CVA算法,该算法降低了传统CVA对几何配准和辐射校正的严格要求,不同变化所得到的结果都处于同一尺度空间,对单限检测法的鲁棒性更好。
[0005] F Bovolo提出了极坐标空间中的CVA算法,该算法从理论上分析了极坐标空间中变化与非变化像素点的分布及其统计特性,给出了极坐标空间中变化矢量的严格定义及分析。此后,F Bovolo又提出了2维CVA和多级CVA算法。但是上述的CVA方法及其改进方法,不能够实现对多光谱图像多种特征的自适应融合,不能充分利用多种特征中的信息,导致运算量比较大,计算效率比较低。
[0006] 除了CVA方法,其它方法也在多光谱图像变化检测中得到了成功的应用。例如,Hichri等在2013年提出了一种实现交式互分割的多光谱图像变化检测方法,该方法使用支持向量机SVM对已知标记的像素进行训练,使用训练结果指导图像分割得到初始变化检测图像,最后使用尔科夫随机场MRFs调整初始变化检测结果图,得到最终的变化检测结果图。但是该方法实现时需要两个步骤,即先对已知标签样本进行训练得到初分结果,再进行优化,使得该算法流程比较复杂。
[0007] 近年来新发展起来的是基于核方法的多光谱图像变化检测算法。GustavoCamps-Valls在2008年首先提出将核方法应用于遥感图像变化检测,通过提取图像的强度信息和纹理信息并构造强度纹理差值核函数实现了遥感图像变化检测。但是该方法处理多光谱图像时,没有考虑多种特性之间的相关性,不能充分挖掘和融合多种图像特征,因而其检测精度有待进一步提高。

发明内容

[0008] 本发明的目的是在针对上述已有算法的不足,提出一种基于核互模态因素分析核融合的多光谱图像变化检测方法,通过分析图像多种特征之间的相关性,充分挖掘和融合图像特征,从而降低运算量,提高运算效率,提高检测精度。
[0009] 为实现上述目的,本发明的实现步骤包括如下:
[0010] (1)给定原始的预处理后的三通道的相同尺寸两时刻光谱图像,记为{Rt,Gt,Bt},t=t0,t1,其中,Rt,Gt,Bt分别表示原RGB空间中的三通道光谱图像,使用RGB-HIS变换将该三通道图像变换至HIS空间,得到HIS通道的两时刻光谱图像{Ht,It,St},其中Ht为t时刻光谱图像的色调分量,St为t时刻光谱图像的饱和度分量,It为t时刻光谱图像的强度分量;
[0011] (2)利用t0,t1两时刻光谱图像的色调 饱和度 和强度 分量,提取光谱图像特征;
[0012] 2a)分别提取 分量光谱图像的灰度值和 分量光谱图像的灰度值,将其堆砌成两时刻光谱图像的颜色特征 和
[0013] 2b)提取强度 分量光谱图像的灰度特征 和纹理特征 将其堆砌成强度纹理特征
[0014] (3)利用两时刻光谱图像的颜色特征 求取颜色差值核函数矩阵KCol;
[0015] (4)利用两时刻光谱图像的强度纹理特征 求取强度纹理差值核函数矩阵KST;
[0016] (5)利用基于核互模态因素分析算法KCFA分别求取颜色差值核函数矩阵KCol的权重系数wc和强度纹理差值核函数矩阵KST的权重系数ws;
[0017] (6)利用步骤5)、步骤4)和步骤3)的结果,构造合成核函数矩阵KF:KF=wc·KCol+ws·KST;
[0018] (7)将步骤6)得到的合成核函数矩阵KF输入到支持向量机SVM中进行检测,得到二值变化检测结果,即完成对多光谱图像的变化检测。
[0019] 本发明具有如下优点:
[0020] 1)本发明由于利用了核方法,可以将实际多光谱图像处理中存在的非线性不可分问题转化为高维空间中的线性可分问题。
[0021] 2)本发明通过将原RGB空间的伪彩色图像变换至HIS空间,充分提取图像的颜色特征和强度纹理特征,有效突出了变化区域及其边界;并将颜色信息和强度纹理信息进行融合,最大限度挖掘两种特征之间的相关性,有效实现了核空间中的特征融合。
[0022] 3)本发明能全面地提取图像特征,可以显著提高变化检测结果的精度。
[0023] 仿真结果表明,本发明与现有的差值核方法相比,运算量低,检测精度高,检测性能更为优越。附图说明
[0024] 图1是本发明的实现流程图
[0025] 图2是用本发明对日本地震灾害实测多光谱图像的变化检测结果图;
[0026] 图3是用本发明对格洛斯特洪灾害实测多光谱图像的变化检测结果图;
[0027] 图4是用本发明对哈萨克斯坦海岸受流沙冲击前后实测多光谱图像的变化检测结果图;
[0028] 图5是用本发明对塔霍湖火灾实测多光谱图像的变化检测结果图。

具体实施方式

[0029] 参照图1,本发明的具体实施过程如下:
[0030] 步骤1.将预处理后的光谱图像变换至HIS空间。
[0031] 给定原始的预处理后的三通道的相同尺寸两时刻光谱图像,记为{Rt,Gt,Bt},t=t0,t1,其中,Rt,Gt,Bt分别表示原RGB空间中的三通道光谱图像;
[0032] 使用RGB-HIS变换将该三通道光谱图像变换至HIS空间,得到HIS通道的两时刻光谱图像{Ht,It,St},其中Ht为t时刻光谱图像的色调分量,St为t时刻光谱图像的饱和度分量,It为t时刻光谱图像的强度分量。
[0033] 其中,RGB空间中的R代表光谱图像的红色分量,G代表光谱图像的绿色分量,B代表光谱图像的蓝色分量;HIS空间中的H代表光谱图像的色调分量,S代表光谱图像的饱和度分量,I光谱图像的代表强度分量。
[0034] 步骤2.提取光谱图像的颜色特征和强度纹理特征。
[0035] 利用t0,t1两时刻光谱图像的色调 饱和度 和强度 分量,提取光谱图像特征,其步骤如下:
[0036] 2a)提取t0时刻光谱图像 分量的灰度值 和 分量的灰度值 用这两个灰度值堆砌构成t0时刻光谱图像的颜色特征 提取t1两时刻光谱图像 分量的灰度值和 分量的灰度值 用这两个灰度值堆砌构成t1时刻光谱图像的颜色特征 表示如下:
[0037]
[0038]
[0039] 2b)提取t0时刻光谱图像的强度 分量的灰度值 和纹理特征 用该灰度值和纹理特征堆砌构成t0时刻光谱图像强度纹理特征 提取t1时刻光谱图像 分量灰度值 和纹理特征 用该灰度值和纹理特征堆砌构成t1时刻光谱图像的强度纹理特征表示如下:
[0040]
[0041]
[0042] 步骤3.构造颜色差值核函数矩阵KCol。
[0043] 利用t0时刻光谱图像的颜色特征 和t1时刻光谱图像的颜色特征 构造颜色差值核函数矩阵KCol,其步骤如下:
[0044] 3a)从t0时刻光谱图像 中任意选取两个像素点 和 从t1时刻光谱图像 中选取与 和 位置相同的两个像素点,记为 和
[0045] 3b)从 中找出像素点 对应的颜色特征 和像素点 对应的颜色特征从 中找出像素点 对应的颜色特征 和像素点 对应的颜色特征
[0046] 3c)利用3b)获得的颜色特征 求取颜色插值核函数的值:
[0047]
[0048] 其中,K(·)表示高斯径向基RBF核函数,ζ代表核宽参数,exp为指数函数;
[0049] 3d)在t0时刻光谱图像 中找出没有被选取过的任意两个像素以及 中相同位置的像素点 按照与3b)相同的方法找到与
对应的颜色特征 再按照与3c)相同的方法计算颜色
差值核函数 的值,
[0050]
[0051] 其中p是像素点的标号,1≤p≤N-1,N是光谱图像像素点个数;
[0052] 3e)重复步骤3d),直到所有像素点被选取过,得到颜色差值核函数的全部值,构成颜色差值核函数矩阵KCol。
[0053] 步骤4.构造强度纹理差值核函数矩阵KST。
[0054] 利用t0时刻光谱图像强度纹理特征 和t1时刻光谱图像强度纹理特征 构造强度纹理差值核函数矩阵KST,其步骤如下:
[0055] 4a)从t0时刻光谱图像 中任意选取两个像素点 和 从t1时刻光谱图像 中选取与 和 位置相同的两个像素点,记为 和
[0056] 4b)从 中找出像素点 对应的强度纹理特征 和像素点 对应的强度纹理特征 从 中找出像素点 对应的强度纹理特征 和像素点 对应的强度纹理特征
[0057] 4c)利用从4b)获得的 求取强度纹理差值核函数的值:
[0058]
[0059] 其中,K(·)表示高斯径向基RBF核函数,ζ代表核宽参数,exp为指数函数;
[0060] 4d)在t0时刻光谱图像 中找出没有被选取过的任意两个像素以及 中相同位置的像素点 按照与4b)相同的方法找到与
对应的强度纹理特征 再按照与4c)相同的方法计算强
度纹理差值核函数 的值,
[0061]
[0062] 其中q是像素点的标号,1≤q≤N-1,N是光谱图像像素点个数;
[0063] 4e)重复步骤4d),直到所有像素点被选取过,得到强度纹理差值核函数的全部值,构成强度纹理差值核函数矩阵KST。
[0064] 步骤5.求取权重系数wc和ws。
[0065] 利用步骤3得到的颜色差值核函数矩阵KCol和步骤4得到的强度纹理差值核函数矩阵KST,利用基于核互模态因素分析算法KCFA分别求取颜色差值核函数矩阵KCol的权重系数wc和强度纹理差值核函数矩阵KST的权重系数ws,其具体步骤如下:
[0066] 5a)通过特征值分解法得到由KST*KCol特征向量构成的第一投影矩阵V和特征值构成的第一对矩阵λ,这三者的关系如下:
[0067] (KST*KCol)*V=λ*V;
[0068] 5b)提取第一投影矩阵V的第一列,记为α,定义颜色差值核函数矩阵KCol融合系数wc:
[0069]
[0070] 其中,T表示矩阵转置,KCol是颜色差值核函数矩阵;
[0071] 5c)通过特征值分解法得到由KCol*KST特征向量构成的第二投影矩阵U和特征值构成的第二对角矩阵λ′,这三者的关系如下:
[0072] (KCol*KST)*U=λ′*U;
[0073] 5d)提取第二投影矩阵U的第一列,记为β,定义强度纹理插值核函数矩阵KST的融合系数ws:
[0074]
[0075] 其中,KST是强度纹理差值核函数矩阵。
[0076] 步骤6.构造合成核函数矩阵KF。
[0077] 利用步骤5得到的权重系数wc和ws,以及步骤3得到的颜色差值核函数矩阵KCol和步骤4得到的强度纹理差值核函数矩阵KST构造合成核函数矩阵:
[0078] KF=wc·KCol+ws·KST。
[0079] 步骤7.完成光谱图像的变化检测。
[0080] 将步骤6得到的合成核函数矩阵KF输入到支持向量机SVM中进行检测,得到二值变化检测结果,即完成对光谱图像的变化检测。
[0081] 本发明效果可以通过以下实验进一步证实:
[0082] 1.实验条件
[0083] 实验仿真环境为:MATLAB R2009b,Intel(R)Core(TM)I5-4590 CPU@3.3GHz,Window 7旗舰版。
[0084] 2.实验内容
[0085] 实验1,应用差值核方法和本发明对日本地震灾害真实多光谱图像进行变化检测,结果如图2,其中,图2(a)是实测日本地震灾害的第一时刻Landsat5TM图像,图2(b)是实测日本地震灾害的第二时刻Landsat5TM图像,图2(c)是变化检测结果参考图,图2(d)是差值核方法的变化检测结果,图2(e)是本发明的变化检测结果;
[0086] 实验2,应用差值核方法和本发明对格洛斯特洪水灾害实测多光谱图像进行变化检测,结果如图3,其中,图3(a)是实测格洛斯特洪水灾害的第一时刻SPOT-XS图像,图3(b)是实测格洛斯特洪水灾害的第二时刻SPOT-XS图像,图3(c)是变化检测结果参考图,图3(d)是差值核方法的变化检测结果,图3(e)是本发明的变化检测结果;
[0087] 实验3,应用差值核方法和本发明对哈萨克斯坦海岸受流沙冲击前后实测多光谱像进行变化检测,结果如图4,其中,图4(a)是实测哈萨克斯坦海岸受流沙冲击前的第一时刻Landsat7ETM+图像,图4(b)是实测城市受洪水哈萨克斯坦海岸受流沙冲击后的第二时刻Landsat7ETM+图像,图4(c)是变化检测结果参考图,图4(d)是差值核方法的变化检测结果,图4(e)是本发明的变化检测结果;
[0088] 实验4,应用差值核方法和本发明对塔霍湖火灾实测多光谱图像进行变化检测,结果如图5,其中,图5(a)是实测塔霍湖火灾的第一时刻Landsat5TM图像,图5(b)是实测塔霍湖火灾的第二时刻Landsat5TM图像,图5(c)是变化检测结果参考图,图5(d)是差值核方法的变化检测结果,图5(e)是本发明的变化检测结果。
[0089] 将上述四个实验的结果进行比较,如表1所示,选取检测精度、Kappa系数作为指标参数,评价本发明检测精度,验证本发明检性能。
[0090] 表1 本发明方法与差值核方法检测结果精度比较
[0091]
[0092] 将四个实验中本发明方法与差值核方法的运行时间进行比较,如表2所示,用以定量的说明本发明与差值核方法变化检测相比的速度优势。
[0093] 表2 本发明方法与差值核方法运行时间对比
[0094]
[0095] 3.实验结果分析
[0096] 从图2、图3、图4,图5和表1、表2可以看到,与基于差值核的变化检测方法相比,本发明的检测性能更为优越,在变化区域具有更好的一致性,且能突出变化边界,检测精度高,运算量低,检测速度快。这是因为本发明充分提取了图像的颜色特征和强度纹理特征,有效突出了变化区域及其边界,并最大限度的挖掘颜色特征和强度纹理两种特征之间的相关性,降低了运算量,有效实现了核空间中的特征融合。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈