CN102024153A - 高光谱图像监督分类方法 - Google Patents
高光谱图像监督分类方法 Download PDFInfo
- Publication number
- CN102024153A CN102024153A CN 201110002042 CN201110002042A CN102024153A CN 102024153 A CN102024153 A CN 102024153A CN 201110002042 CN201110002042 CN 201110002042 CN 201110002042 A CN201110002042 A CN 201110002042A CN 102024153 A CN102024153 A CN 102024153A
- Authority
- CN
- China
- Prior art keywords
- node
- classification
- pixel
- class
- mark
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种高光谱图像监督分类方法,主要解决现有方法对高光谱图像像素之间的空间信息利用不充足的问题,其分类过程为:对于待分类的高光谱图像,依据训练样本集建立二叉树模型,得到内部节点索引表;依据该索引表,依序找出本次待分类的内部节点标号和标号对应的高光谱图像像素点的位置;逐个计算相应位置处的每个像素分别属于本次内部节点的左、右子节点的似然概率值,得到高光谱图像的初始类标记;依据交替优化算法,更新初始类标记得到最终类标记;分别记录最终类标记中的两个子节点包含像素点的位置,存储到相应的子节点中;依次处理每个内部节点对应的高光谱图像数据,直至处理完所有内部节点。本发明可用于遥感图像分类与识别。
Description
技术领域
本发明属于图像处理领域,涉及一种高光谱图像监督分类方法,可用于遥感图像分类与识别。
背景技术
高光谱图像的分类与识别是遥感图像处理领域的一个新的重要分支,近年来在植被研究、矿物识别以及水文和海洋勘测等领域有着广泛的应用。高光谱遥感,是指具有很高的光谱分辨率的遥感技术,它能够在很窄的连续波段上记录地物,如地表、水域、大气的光谱反射特性,使得获取待测地物类别的细节成为现实。成像光谱仪能够为每个像素提供数十个至数百个窄波段,高光谱图像的每一帧都是对同一场景在不同波长上的成像,从而产生一条完整而连续的光谱曲线,提供所观测地物丰富的光谱信息。
针对高光谱遥感技术的上述特点,常用的高光谱图像的分类方法有:
从光谱角度出发,有SAM算法和混合像素分解技术。SAM方法通过比较目标光谱曲线和已有的样本库中的各类平均光谱曲线之间的相似程度来判定待分类光谱的类别归属;而混合像素光谱分解技术是通过将每个像素点分解为不同的基本单元即端元,来实现类别的划分。该类方法通常需要预先建立待分类图像所包含地物的光谱库,因而在实现上比较困难。
从统计模型角度出发,有基于统计距离的方法和基于概率模型的方法,最小距离分类法是基于统计距离的方法中最简单的一种,该方法通过比较待分类像素和利用训练样本计算出的各类中心点之间的距离,来实现像素点的类别划分。而基于概率模型的方法中最典型的就是高斯极大似然分类方法,该方法将像素向量作为特征向量,在假定各类服从正态分布的基础上,根据有关概率判决函数的贝叶斯准则进行分类。该类方法应用广泛,但由于算法的复杂度和数据维数的平方成正比,且类别参数由于训练样本的相对较少无法可靠估计。
还有基于决策树的分类方法,如二元层级分类器BHC,该方法假设高光谱数据中各类服从一定分布,将贝叶斯规则用在每层的决策中,不断的将一个K类分类问题转化为K-1个的二元分类问题,从而解决分类问题。该方法相比之前的方法,计算简单且易于实现,但由于没有很好的加入先验信息,忽略了图像中像素之间的相互关系,导致分类过程中对空间信息的利用不足,限制了分类精度的提高。
发明内容
本发明的目的在于克服上述已有技术的不足,提供一种树状的马尔可夫随机场模型TS-MRF的高光谱图像监督分类方法,以有效地利用图像中像素之间的空间关系,提高分类精度。
实现本发明目的的技术方案是:使用简洁高效的二叉树结构,结合马尔可夫随机场对图像像素之间的空间信息的利用,对高光谱图像进行监督分类,具体实现步骤如下:
(1)对于待分类的高光谱图像,依据训练样本集,建立二叉树模型,得到节点标号,将内部节点标号及其包含的类别组存储于索引表中:
1a)计算训练样本集中各个类别的均值向量,将其作为各类特征;
1b)使用欧式距离作为衡量各类别之间差异的度量,计算出每两个类别i1,i2之间的空间距离diff(i1,i2),所有类别之间的空间距离构成一个大小为K×K的可分性判别矩阵D,K为类别总数;
1c)找出可分性判别矩阵中最大值max(D)对应的两个类别i,j,将其作为两个类别组C1和C2的初始值,即C1={i},C2={j};逐个比较每个剩余类别i′与类别i和类别j的空间距离,若i′与类别i的空间距离diff(i′,i)小于i′与类别j的空间距离diff(i′,j),则将类别i′划到第一类别组C1中;否则将i′划到第二类别组C2中,即下式:
其中i′∈{1,2,...,K},i′≠i,j,U表示所有i′的集合;
1d)按照1c)中的方法,不断将每个类别组细分为更小的类别组,直至所有类别组都只包含一个最终类别,树模型即建立成功;树模型建立过程中的每个类别组以及最终类别统称为树的节点,其中类别组称为内部节点,最终类别称为叶子节点;给所有节点分配节点标号,树的根节点标号为1,该根节点标号1对应整个原始图像,包含图像中所有像素点的位置;对树中任一个节点t,其左子节点l(t)的节点标号为2t,右子节点r(t)的节点标号为2t+1,将内部节点标号及其包含的类别组存储于索引表中;
(2)将高光谱图像数据和树节点对应,依据树模型中的内部节点索引表,依序找出本次待分类的内部节点标号和标号对应的图像像素点的位置;
(3)根据得到的图像像素点位置,逐个计算相应位置处的每个像素分别属于本次内部节点的左、右子节点的似然概率值,将两个似然概率值中最大值对应的子节点的类别作为该像素类别,由此得到图像的初始类标记;
(4)依据交替优化算法,将得到的初始类标记进行类标记更新,得到本次的最终类标记;
(5)得到最终类标记后,分别记录本次分类后的两个子节点包含像素点的位置,存储到相应的子节点中;
(6)重复步骤(2)到步骤(5),依次处理每个内部节点对应的图像数据,直至处理完所有内部节点对应的数据,完成对高光谱图像的分类。
本发明与现有的技术相比有以下优点:
1)本发明由于采用了马尔可夫随机场,更好地利用了像素之间的空间依赖关系,使得图像的空间信息得到充分利用,提高了分类精度;
2)本发明由于采用了二叉树状的结构,使得分类的过程可以递归地进行,简洁有效;
3)本发明由于采用了一种可移植的树模型建立方法,使得基于TS-MRF模型的方法的树模型建立过程有章可依,给予树模型建立以有效的理论基础,增强该算法的适用性。
附图说明
图1是本发明的流程图;
图2是本发明所采用的高光谱图像三个波段合成的示意图;
图3是本发明所建立的二叉树模型图;
图4是本发明对图2所示高光谱图像的仿真结果图。
具体实施方式
参照图1,本发明的具体实现步骤如下:
步骤一,对于图2所示的待分类高光谱图像Y={y1,y2,...,yn},n为图像的总像素数,依据从高光谱图像中得到的带有样本标记的训练样本集,按照逐次二元划分的方式建立如图3所示的二叉树模型,得到树模型中各个节点的节点标号,将其中的内部节点标号及其包含的类别组存储于索引表中,具体步骤如下:
1a)计算训练样本集中各类别所包含的高光谱图像数据的均值向量,将其作为各类特征;
1b)使用欧式距离作为衡量各类别之间差异的度量,计算出每两个类别i1,i2之间的空间距离diff(i1,i2),所有类别之间的空间距离构成一个大小为K×K的可分性判别矩阵D,K为类别总数;
1c)找出可分性判别矩阵D中最大值max(D)对应的两个类别i,j,将其作为两个类别组C1和C2的初始值,即C1={i},C2={j};逐个比较每个剩余类别i′与类别i和类别j的空间距离,若i′与类别i的空间距离diff(i′,i)小于i′与类别j的空间距离diff(i′,j),则将类别i′划到第一类别组C1中;否则将i′划到第二类别组C2中,即下式:
其中i′∈{1,2,...,K},i′≠i,j,U表示所有i′的集合;
1d)按照1c)中的方法,不断将每个类别组细分为更小的类别组,直至所有类别组都只包含一个最终类别,树模型即建立成功;树模型建立过程中的每个类别组以及最终类别统称为树的节点,其中类别组称为内部节点,最终类别称为叶子节点;给所有节点分配节点标号,树的根节点标号为1,该根节点标号1对应整个原始图像,包含高光谱图像中所有像素点的位置S={1,2,..,n};对树中任一个节点t,其左子节点l(t)的节点标号为2t,右子节点r(t)的节点标号为2t+1,将内部节点标号及其包含的类别组存储于索引表ID中。
步骤二,将高光谱图像数据yt和树中的节点t对应,依据树模型中的内部节点索引表ID,依序找出本次待分类的内部节点标号ID(t)和标号对应的高光谱图像像素点的位置St,除根节点之外的其它内部节点对应的高光谱图像像素点的位置均为原始图像位置集的一部分,即St∈S。
步骤三,根据得到的高光谱图像像素点位置St,逐个计算其中的每个位置s处的像素分别属于本次内部节点t的左子节点l(t)和右子节点r(t)的似然概率值,将两个似然概率值中最大值对应的子节点的类别作为该像素的类别,由此得到高光谱图像的初始类标记,具体步骤如下:
其中yt和xt分别为节点t对应的像素集合和相应的类标记集合,和分别为yt和xt中位于位置s处的像素和类标记,w(t)为节点t的祖先节点集,指树中先于节点t产生的所有节点的集合,为祖先节点集w(t)的类标记估计,为中位于位置s处的类标记估计,l(t)和r(t)分别为节点t的左、右子节点;
3b)将节点的似然概率值分解为各个像素似然概率值的乘积之后,对于节点t中每一个像素点和分别表示该像素属于节点t的左子节点l(t)和右子节点r(t)的似然概率,在服从高斯分布的条件下,似然概率的计算式如下:
3c)由<3>式计算出似然概率后,对于节点t中每一个像素点,若该像素属于节点t的左子节点l(t)的似然概率大于属于右子节点r(t)的似然概率则该像素的类标记为左子节点l(t)的类别;否则,像素的类标记为右子节点r(t)的类别。
其中和分别为节点t对应数据的类标记xt和边缘罚函数βt在交替优化算法的第k次估计中的结果,为节点t中,给定祖先节点集合的类标记估计,先验分布在交替优化算法的第k次估计中的结果;为节点t中,给定类标记xt以及祖先节点集合的类标记估计似然分布在交替优化算法的第k次估计中的结果;
令在设定的经验值范围[a,b]内,以固定步长λ增加则使得<5>式等号右边的表达式值最大的即MPL估计结果;其中t为当前节点标号,St为节点t对应像素点的位置集,s为St中的一个像素点的位置,η(S)为位置s的邻域,r为邻域η(s)内任一点,为位置s的邻域中与其具有相同类标记的像素个数,L为标记空间;
使得该式最小的类标记即该像素点的类标记,由此将高光谱图像数据划分为两类,得到初始类标记;其中和分别为节点t中的像素所属的类别对应的高光谱图像数据的均值向量和协方差矩阵,为位置s的邻域中与其具有相同类标记的像素个数;
4c3)检验是否满足迭代停止条件,即类标记发生变化的像素点总数小于预先设定的阈值,若满足停止条件,则结束ICM算法;否则重复4c1)-4c3)的步骤;
本发明的效果可通过以下仿真进一步说明:
1.仿真内容
应用本发明的方法和现有的高斯极大似然分类方法,对一幅高光谱图像进行分类实验,主要从各个类别的内部一致性、类别边缘的完整性和边缘定位准确度来评价这两种方法的性能。
2仿真实验结果
用本发明方法以及高斯极大似然分类方法对图2所示的高光谱图像进行分类,分类结果如图4所示。其中图4(a)为高斯极大似然分类的分类结果,图4(b)为本发明的分类结果,图4(c)为高斯极大似然分类的标记区域分类结果,4(d)为本发明的标记区域分类结果。由图4可见,本发明中,各个类别的一致性较为良好,类别边缘完整,且边缘定位准确,分类效果良好。相比之下,高斯极大似然分类方法各个类别的一致性也较为良好,但是类别边缘的定位很不准确,出现很多错误边缘,导致分类结果不佳。
Claims (3)
1.一种高光谱图像监督分类方法,包括如下步骤:
(1)对于待分类的高光谱图像,依据训练样本集,建立二叉树模型,得到节点标号,将内部节点标号及其包含的类别组存储于索引表中:
1a)计算训练样本集中各个类别的均值向量,将其作为各类特征;
1b)使用欧式距离作为衡量各类别之间差异的度量,计算出每两个类别i1,i2之间的空间距离diff(i1,i2),所有类别之间的空间距离构成一个大小为K×K的可分性判别矩阵D,K为类别总数;
1c)找出可分性判别矩阵中最大值max(D)对应的两个类别i,j,将其作为两个类别组C1和C2的初始值,即C1={i},C2={j};逐个比较每个剩余类别i′与类别i和类别j的空间距离,若i′与类别i的空间距离diff(i′,i)小于i′与类别j的空间距离diff(i′,j),则将类别i′划到第一类别组C1中;否则将i′划到第二类别组C2中,即下式:
其中i′∈{1,2,...,K},i′≠i,j,U表示所有i′的集合;
1d)按照1c)中的方法,不断将每个类别组细分为更小的类别组,直至所有类别组都只包含一个最终类别,树模型即建立成功;树模型建立过程中的每个类别组以及最终类别统称为树的节点,其中类别组称为内部节点,最终类别称为叶子节点;给所有节点分配节点标号,树的根节点标号为1,该根节点标号1对应整个原始图像,包含图像中所有像素点的位置;对树中任一个节点t,其左子节点l(t)的节点标号为2t,右子节点r(t)的节点标号为2t+1,将内部节点标号及其包含的类别组存储于索引表中;
(2)将高光谱图像数据和树节点对应,依据数模型中的内部节点索引表,依序找出本次待分类的内部节点标号和标号对应的高光谱图像像素点的位置;
(3)根据得到的高光谱图像像素点位置,逐个计算相应位置处的每个像素分别属于本次内部节点的左、右子节点的似然概率值,将两个似然概率值中最大值对应的子节点的类别作为该像素类别,由此得到高光谱图像的初始类标记;
(4)依据交替优化算法,将得到的初始类标记进行类标记更新,得到本次的最终类标记;分别记录本次分类后的两个子节点包含像素点的位置,存储到相应的子节点中;
(5)重复步骤(2)到步骤(4),依次处理每个内部节点对应的高光谱图像数据,直至处理完所有内部节点对应的数据,完成对高光谱图像的分类。
2.根据权利要求1所述的高光谱图像分类方法,其中步骤(3),按如下步骤进行:
其中yt和xt分别为节点t对应的像素集合和相应的类标记集合,和分别为yt和xt中位于位置s处的像素和类标记,w(t)为节点t的祖先节点集,指树中先于节点t产生的所有节点的集合,为祖先节点集w(t)的类标记估计,为中位于位置s处的类标记估计,l(t)和r(t)分别为节点t的左、右子节点;
3b)将节点的似然概率值分解为各个像素似然概率值的乘积之后,对于节点t中每一个像素点和分别表示该像素属于节点t的左子节点l(t)和右子节点r(t)的似然概率,在服从高斯分布的条件下,似然概率的计算式如下:
3.根据权利要求1所述的高光谱图像分类方法,其中步骤(4)所述的依据交替优化算法,将得到的初始类标记进行类标记更新,按如下步骤进行:
其中和分别为节点t对应数据的类标记xt和边缘罚函数βt在交替优化算法的第k次估计中的结果,为节点t中,给定祖先节点集合的类标记估计先验分布在交替优化算法的第k次估计中的结果;为节点t中,给定类标记xt以及祖先节点集合的类标记估计似然分布在交替优化算法的第k次估计中的结果;
令在设定的经验值范围[a,b]内,以固定步长λ增加则使得<5>式等号右边的表达式值最大的即MPL估计结果;其中t为当前节点标号,St为节点t对应像素点的位置集,s为St中的一个像素点的位置,η(s)为位置s的邻域,r为邻域η(s)内任一点,为位置s的邻域中与其具有相同类标记的像素个数,L为标记空间;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100020425A CN102024153B (zh) | 2011-01-06 | 2011-01-06 | 高光谱图像监督分类方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100020425A CN102024153B (zh) | 2011-01-06 | 2011-01-06 | 高光谱图像监督分类方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102024153A true CN102024153A (zh) | 2011-04-20 |
CN102024153B CN102024153B (zh) | 2012-07-04 |
Family
ID=43865433
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011100020425A Active CN102024153B (zh) | 2011-01-06 | 2011-01-06 | 高光谱图像监督分类方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102024153B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102902491A (zh) * | 2012-08-29 | 2013-01-30 | 北京地拓科技发展有限公司 | 遥感影像数据写入及读取的方法及装置 |
CN103984960A (zh) * | 2014-05-30 | 2014-08-13 | 中国科学院遥感与数字地球研究所 | 一种遥感图像样本智能采集方法 |
CN104102929A (zh) * | 2014-07-25 | 2014-10-15 | 哈尔滨工业大学 | 基于深度学习的高光谱遥感数据分类方法 |
CN104298999A (zh) * | 2014-09-30 | 2015-01-21 | 西安电子科技大学 | 基于递归自动编码的高光谱特征学习方法 |
CN104778482A (zh) * | 2015-05-05 | 2015-07-15 | 西安电子科技大学 | 基于张量半监督标度切维数约减的高光谱图像分类方法 |
CN107507155A (zh) * | 2017-09-25 | 2017-12-22 | 北京奇虎科技有限公司 | 视频分割结果边缘优化实时处理方法、装置及计算设备 |
-
2011
- 2011-01-06 CN CN2011100020425A patent/CN102024153B/zh active Active
Non-Patent Citations (2)
Title |
---|
《电子学报》 20040725 汪西莉等 基于不完全分层MRF的非监督图象分割 1086-1089 1-3 第32卷, 第07期 * |
《软件学报》 20030923 汪西莉等 一种分层马尔可夫图像模型及其推导算法 1558-1563 1-3 第14卷, 第09期 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102902491A (zh) * | 2012-08-29 | 2013-01-30 | 北京地拓科技发展有限公司 | 遥感影像数据写入及读取的方法及装置 |
CN102902491B (zh) * | 2012-08-29 | 2015-09-09 | 北京地拓科技发展有限公司 | 遥感影像数据写入及读取的方法及装置 |
CN103984960A (zh) * | 2014-05-30 | 2014-08-13 | 中国科学院遥感与数字地球研究所 | 一种遥感图像样本智能采集方法 |
CN103984960B (zh) * | 2014-05-30 | 2017-12-08 | 中国科学院遥感与数字地球研究所 | 一种遥感图像样本智能采集方法 |
CN104102929A (zh) * | 2014-07-25 | 2014-10-15 | 哈尔滨工业大学 | 基于深度学习的高光谱遥感数据分类方法 |
CN104102929B (zh) * | 2014-07-25 | 2017-05-03 | 哈尔滨工业大学 | 基于深度学习的高光谱遥感数据分类方法 |
CN104298999A (zh) * | 2014-09-30 | 2015-01-21 | 西安电子科技大学 | 基于递归自动编码的高光谱特征学习方法 |
CN104298999B (zh) * | 2014-09-30 | 2017-08-25 | 西安电子科技大学 | 基于递归自动编码的高光谱特征学习方法 |
CN104778482A (zh) * | 2015-05-05 | 2015-07-15 | 西安电子科技大学 | 基于张量半监督标度切维数约减的高光谱图像分类方法 |
CN104778482B (zh) * | 2015-05-05 | 2018-03-13 | 西安电子科技大学 | 基于张量半监督标度切维数约减的高光谱图像分类方法 |
CN107507155A (zh) * | 2017-09-25 | 2017-12-22 | 北京奇虎科技有限公司 | 视频分割结果边缘优化实时处理方法、装置及计算设备 |
CN107507155B (zh) * | 2017-09-25 | 2020-02-18 | 北京奇虎科技有限公司 | 视频分割结果边缘优化实时处理方法、装置及计算设备 |
Also Published As
Publication number | Publication date |
---|---|
CN102024153B (zh) | 2012-07-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103258214B (zh) | 基于图像块主动学习的遥感图像分类方法 | |
Zhong et al. | Modeling and classifying hyperspectral imagery by CRFs with sparse higher order potentials | |
CN103996047B (zh) | 基于压缩谱聚类集成的高光谱图像分类方法 | |
CN102024153B (zh) | 高光谱图像监督分类方法 | |
CN103366367B (zh) | 基于像素数聚类的模糊c-均值灰度图像分割方法 | |
CN102073879B (zh) | 基于半监督学习的海岸海洋遥感影像特征地类的识别方法 | |
CN105808752B (zh) | 一种基于cca和2pknn的自动图像标注方法 | |
CN105956612B (zh) | 基于主动学习和半监督学习的高光谱图像分类方法 | |
CN109359525B (zh) | 基于稀疏低秩的判别谱聚类的极化sar图像分类方法 | |
CN103955709B (zh) | 基于加权合成核与tmf的极化sar图像分类方法 | |
CN103745233B (zh) | 基于空间信息迁移的高光谱图像分类方法 | |
CN103310230A (zh) | 联合解混及自适应端元提取的高光谱图像分类方法 | |
CN104376335A (zh) | 一种基于信息熵的半监督高光谱遥感影像分类方法 | |
CN102968640A (zh) | 基于Freeman分解和数据分布特征的极化SAR图像分类方法 | |
CN104036497A (zh) | 基于方差局部系数的图切交互式图像分割方法 | |
CN105160351A (zh) | 基于锚点稀疏图的半监督高光谱分类方法 | |
CN102999761A (zh) | 基于Cloude分解和K-wishart分布的极化SAR图像分类方法 | |
Polewski et al. | Active learning approach to detecting standing dead trees from ALS point clouds combined with aerial infrared imagery | |
CN113723492A (zh) | 一种改进主动深度学习的高光谱图像半监督分类方法及装置 | |
CN103365985A (zh) | 一种类别自适应的极化sar分类方法 | |
CN105069471A (zh) | 基于模糊标签的高光谱数据子空间投影和分类方法 | |
CN111222575A (zh) | 一种基于hrrp目标识别的klxs多模型融合方法及*** | |
Li et al. | A label propagation method using spatial-spectral consistency for hyperspectral image classification | |
CN105279523B (zh) | 一种结合混合像元分解与主动学习的半监督分类方法 | |
CN104732246A (zh) | 一种半监督协同训练高光谱图像分类方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |