一种高光谱图像特征的并行提取与分类方法
【专利摘要】本发明提出了一种高光谱图像特征的并行提取与分类方法，通过三维空间域以及频域联合分析，对三维和高光谱图像进行并行特征提取，然后对这些特征进行增强和融合，可以充分利用高维数据的丰富信息以及这些特征之间的结构关系，大大减少噪声对准确率的影响。本发明对各种相位、方向、频域以及三维空间编码方法进行深入研究，然后从提取各Gabor特征在方向与频率空间结构关系出发，在联合提取空间?光谱域和频域的信号变化特征同时，开拓性地将智能算法引入波段选择，从而可以在保证识别准确率的前提下减少冗余信息，提高识别效率，具有较广泛的应用和推广前景。
【专利说明】
-种高光谱图像特征的并行提取与分类方法
技术领域
[0001] 本发明属于图像处理技术领域，具体设及一种高光谱图像特征的并行提取与分类 方法。
【背景技术】
[0002] 高光谱成像技术是近二十年来发展起来的基于非常多窄波段的影像数据技术，其 最突出的应用是遥感探测领域，并在越来越多的民用领域有着广泛的应用前景。高光谱成 像主要通过不同波长的光源对物体进行拍照，获取物体在不同波段下的照片，从而提取更 丰富的有用信息。彩色成像通过红、蓝、绿=个通道获得多光谱图像，而高光谱技术通常将 可见光、近红外分割成上百个波段活动高达数百张的高光谱图像。高光谱遥感的光谱分辨 率一般达到纳米级，光谱波段数为几十到数百，提供了窄波段连续光谱信息，对物体的测谱 能力非常强。
[0003] 随着高光谱技术在遥感领域的成功应用，现在该技术也开始在食品安全，医学诊 断、人脸、掌纹识别W及材料探测等领域广泛运用。食品安全和材料检测主要通过高光谱图 像分析食品、物体的光谱，来判断是否有异常物质的加入，从而确定食品是否安全W及材料 的种类。医学诊断可W通过高光谱图像获得人体组织的光谱，从而判断该组织是否存在病 变。鉴于不同人的皮肤所含的胡萝h素、黑色素和血色素各不相同，在不同光谱下会呈现出 不同的反射特性，研究人员也已经开始通过高光谱图像提取特征对人脸和掌纹进行识别。
[0004] 目前高光谱图像识别分类的研究还处于起步阶段，研究工作主要集中为遥感图像 的地物分类领域。如何从运些海量遥感数据中抽取代表性特征，从而实现准确的地物分类， 是目前高光谱遥感影像研究的热点。但在高光谱人脸、掌纹等其他领域，相关研究工作还比 较少。目前开展高光谱图像识别分类的研究，可W大大促进该技术在诸如高光谱人脸、高光 谱掌纹W及医学诊断领域的推广，从而提升基于高光谱信息的物体识别分类技术进步。
[0005] 在高光谱遥感影像中，像元的光谱是最简单、可用于表征地物类别的特征。但是由 于光谱数据中噪声比较多，基于光谱数据的分类效果往往比较差。为抑制光谱数据中的噪 声，主成分分析（Principal Component Analysis,PCA)，线性判别分析化inear Discriminant Analysis , LDA似及独立成分分析（Independent Component Analysis , ICA)等子空间方法都被用来对数据进行分类前的预处理。除了光谱，研究者们也开始提取 诸如纹理、形状之类的空间域特征，从而提高不同地物之间的可区分性。实验表明，纹理和 形状等空间域特征可有效地和光谱形成互补，从而大大增加地物分类准确率。鉴于Gabor小 波在二维图像特征提取领域的成功，不少学者在二维切片/帖上提取二维Gabor特征，用于 =维图像的分析与识别。由于只是提取了二维空间域信息，该方法显然没有充分利用=维 数据所包含的所有信息。除了在空间域和光谱域分别提取特征，一种更直接有效的方法是 提取空间-光谱域S维空间的联合变化特征。非专利文献1 (S. Jia，L.化en，and Q丄i，" Gabor feature-based collaborative representation for hyperspectral imagery classification,"IEEE Transactions on Geoscience and Remote Sensing,vol.53, PP. 1118-1129，2015)通过S维Gabor小波直接在S维空间提取信号变化信息，率先将S维 Gabor特征应用于高光谱遥感图像分类。该方法将S维数据和S维Gabor小波卷积后的幅值 结果向量化，然后拼接成为一个长特征向量。实验表明，直接从S维数据提取的S维Gabor 特征可W包含更丰富信息，从而获得比二维Gabor特征高得多的准确率。Jia和化en等人将 向量连接形式构造的S维Gabor特征分别和基于LlW及L2优化准则的稀释表示框架进行了 结合，用于高光谱遥感图像分类。实验结果表明，和采用光谱特征时不同，在采用区分度高 的S维Gabor特征时，依据L2优化准则的协同表示方法能够在大大提高分类效率的同时，取 得和采用Ll稀疏表示方法相似的准确率。从上述文献可W看出，目前WS维Gabor小波为代 表的空间-光谱域联合变化特征提取方法已经在高光谱遥感图像分类中开始广泛应用，各 种实验对比表明，该特征可W取得比光谱和二维纹理好得多的准确性能。
[0006] 相对于可见光人脸识别，高光谱人脸识别研究的起步较晚、研究成果也较少。美国 Rochester技术学院的Munsell色彩科学实验室在1999年率先对高光谱人脸数据库的采集 进行了尝试。随后，卡内基梅隆大学在2002年采集了 18人的数据库用于人脸非对称性的研 究。加州大学埃尔文分校计算机视觉实验室2003年在700到1000纳米的近红外波段采集了 200个人的高光谱人脸数据，通过前额、左右脸颊和嘴唇等采样区域的平均光谱特性进行识 另IJ，实验结果表明高光谱成像在姿态和表情变化下仍能取得令人满意的识别效果。该实验 室的研究团队2007年进一步对室外光照情况进行了研究，在原型集(Gallery)和测试集 (Probe)的光照存在巨大变化的情况下，高光谱方法同样取得了较高的准确率。该方法只采 用了人脸在谱域的信号变化特性，没有充分利用空域信息。
[0007] 基于高光谱图像的掌纹识别也开始起步。Guo等人在非专利文献2 (Z. Guo, D.Zhang，L.Zhang,and w.Liu，"Feature band selection for online multispectral palmprint recognition, IEEE Transactions on Information Forensics and Security，vol.7，pp. 1094-1099,2012.)中率先开展相关数据库的采集和研究。该方法主要 通过提取选择波段上的二维特征来进行识别。调研表明，相对于数目众多的可见光人脸、掌 纹识别研究团队和成果，目前高光谱人脸/掌纹识别方面的研究还处于起步阶段，需要大量 理论和应用上的尝试。通过S维Gabor小波提取空间-光谱域联合变化特征，是未来的一个 发展方向。
[0008] 从高光谱图像中提取的S维空间-光谱域Gabor特征维度巨大。例如，对于大小为 100 X 100 X 100的=维数据，52个小波抽取的特征维度可高达52,000,000。如此高维特征会 给后续分类带来巨大的计算量。例如在稀疏表示框架下，原始S维Gabor特征由于维度太大 而无法让现有的求解算法获得稀疏解，必须要通过主成分分析来降维。相对二维Gabor特 征，采用向量化拼接方式获得的S维Gabor特征在维度上增加了数百倍。目前对运些特征向 量进行选择是主要的有效手段，还未能见到有更有效的相关融合方法研究。一个可行的研 究方向是将每个S维Gabor小波提取的特征分别进行处理，从而大大降低每个任务所要处 理的特征维度和计算复杂度，后续可W通过多任务联合稀疏表示框架或者多视角张量结构 把每个任务的结果/特征融合起来。
[0009] 相对二维可见光图像的广泛应用，高光谱图像无论在基础理论、还是应用领域的 研究都还处于起步阶段。目前高光谱图像的特征提取正在从光谱、二维纹理特征的融合向 空间-光谱域联合变化特征的直接提取过度。在运个过程中，采用S维Gabor小波提取的S 维Gabor特征在高光谱遥感图像分类、高光谱人脸/掌纹识别应用中表现出了优异的准确性 能，正在日益受到相关学者的关注。虽然目前的S维Gabor特征性能好，但是也存在着特征 维度大、冗余性强等缺点。不同Gabor特征的融合目前主要采取向量化和向量拼接方式，急 需对能够保持方向域、频域W及=维空间结构关系的融合方法进行研究。该发明对高维 Gabor特征选择、不同小波提取的Gabor特征融合方法进行深入研究，不仅可W大大丰富和 完善空间-光谱域联合变化特征提取理论，还可W有力推动该工具和日臻完善的高光谱成 像技术的结合，从而促进基于高光谱数据的生物特征识别、医学诊断和遥感监测领域的技 术进步。
【发明内容】
[0010] 为解决现有技术中存在的问题，本发明提出了一种高光谱图像特征的并行提取与 分类方法，大大减少噪声对准确率的影响，在保证识别准确率的前提下减少冗余f目息，提局 识别效率，具有较广泛的应用和推广前景。
[0011] 本发明具体通过如下技术方案实现：
[0012] -种高光谱图像特征的并行提取与分类方法，包括W下步骤：
[0013] A.获取高光图图像；
[0014] B.对所述高光图图像进行基于S维Gabor小波的空间-光谱域联合变化的并行特 征提取;具体为：
[0015] 给定S维高光谱数据Wx，y，A)，S维小波分析通过S维Gabor小波函数Wf,e,*(x, y，A)和V(x，y，A)的内积来同时获得信号的空域和频域信息，其中，参数f为小波的中屯、频 率，d)为该频率在S维频域与W轴的夹角，0为该频率在S维频域的U-V平面内投影与U轴的 夹角；在采用Gabor小波抽取特征时，需要预先设计多个不同中屯、频率fi W及方向(0,，&)的 小波
I然后将其和V(x，y，A)进行卷积操作，卷积结 果
代表了在S维空间（x，y，A)、频率为方向为(巧，稱)的信号变 化信息；给定N个S维Gabor小波Wn(n= 1，2，-'N)，通过卷积获得N个同样大小的S维数据 F 0聲。fe，j，:躬:;从S维空间、频域W及方向域S个角度对所述N个S维Gabor特征进行空 间编码W及频率、方向编码后，V(x，y，A)获得K个编码后的特征图像gi(x，y，A)，…gk(x，y, 入)..心(又山入）
[0016] C.结合波段选择对所述步骤B提取的S维Gabor特征进行筛选，通过所述波段选择 去除高光谱数据中冗余信息；
[0017] D.构造并学习基于多视角的高阶Gabor张量，实现S维Gabor特征融合；
[0018] E.构造核函数实现基于支持向量机的图像分类。
[0019] 进一步地，所述空间编码的主要思想是将中屯、像素的响应与相邻像素的响应进行 差异性比较，然后对大小关系进行0/1编码，最后再转换成整数，通过区域直方图统计对 Gabor特征进行增强。
[0020] 进一步地，所述频率、方向编码的核屯、思想是将中屯、像素的Gabor响应和相邻频 率、方向的响应进行差异性比较，然后对大小关系进行0/1编码，最后再转换成整数，通过区 域直方图统计对Gabor特征进行增强。
[0021] 进一步地，所述步骤C具体包括：
[0022] Stepl:设计L个不同中屯、频率和方向的S维Gabor小波{ Wi, 1 = 1，..，L}，将其和B 个波段的多光谱数据Wx，y，A)卷积，获得BXL个卷积图像：|Mi(x，y，Ai)，…化(x，y，Ai)…Mi (x,y,Ab)...]V[L(x,y,^b)...Mi(x,y,AB)...ML(x,y,AB)};
[0023] Step2:依次提取每个波段Ab的L个卷积图像{Mi(x，y，Ab)…化^，7，人〇}，提取维度 较低的特征，设计评价函数进行排序，找出最优波段;将该最优波段从候选集合移出，加入 到已选集合；
[0024] Step3:从候选集合选取一个波段加入到已选集合，使得已选集合的评价函数值得 到最大的增加；
[0025] Step4:重复上述过程Stepl-3,直到增加量小于设定阔值或者迭代次数达到预先 设定值。
[0026] 进一步地，所述步骤D具体为:波段选择后提取的S维Gabor特征图像表示为K个S 阶张量X={xl，X2…XK}，其中义e臾WxS利用稀疏非负多视角张量分解算法，同时学习K个 S阶张量的共同因子和各个视角上的潜在因子:一个S阶张量Xk的CP分解优化模型写为
[0027]
[002引为简便起见，脚("U(2)U(3嗦示
，其中啤为UW的第j个列向量，u(。 通常称为第i模态的因子矩阵，拟优化的目标函数如下：
[0029]
[0030]
[0031] 其中，巧々为第i个视角下第j个模态上的因子矩阵，[V阳1]表示第一模态上的因 子矩阵，护为共同因子矩阵(运里也可W根据实际情况考虑其他模态）；M ||2,1为^2.1范数， 用于特征的稀疏化学习；所述目标函数首先根据范数的分析方法，将其转换为矩阵迹的 表达形式，进而利用交替最小二乘和半二次法结合进行优化求解得到因子矩阵护和1/,W，即 为多视角张量分解提取的潜在因子矩阵。
【附图说明】
[0032] 图1是S维频域的示意图；
[0033] 图2是Gabor小波频域示意图；
[0034] 图3是尺度、方向编码示意图；
[0035] 图4是多视角张量分解示意图。
【具体实施方式】
[0036] 下面结合【附图说明】及【具体实施方式】对本发明进一步说明。
[0037] 本发明主要研究如何从高光谱影像中提取空间-光谱域联合变化特征，并通过智 能优化算法对运些特征进行筛选与融合，去除冗余特征，最终实现准确和高效的识别分类 效果。本发明重点关注基于S维Gabor小波的空间-光谱域联合变化特征提取，具体如下：
[0038] (1)基于S维Gabor小波的空间-光谱域联合变化的并行特征提取
[0039] 本发明在深入研究=维空间域和频域同时提取信号变化特性机理的同时，对运些 编码方法进行优化，重点关注能够提取Gabor小波在方向、频率W及S维空间结构关系的编 码方法。
[0040] 目前S维Gabor特征主要W提取高光谱数据和Gabor小波的卷积响应幅值为主，将 =维空间各个点的小波响应直接连接成一个向量，用于后续的分类识别。直接向量化和简 单连接将各个特征独立对待，往往丢失了各个特征在方向、频率W及=维空间的结构关系， 从而降低特征的识别精度。为此，一个更好的办法是对各个Gabor特征在方向、频率W及S 维空间的结构关系进行编码，从而获得更具分辨能力的编码特征。
[0041] =维空间编码的主要思想是将中屯、像素的响应与诸如3 X 3 X 3邻域内的相邻像素 的响应进行差异性比较，然后对大小关系进行0/1编码，频率、方向编码的核屯、思想是将中 屯、像素的Gabor响应和相邻频率(尺度）、方向的响应进行差异性比较，然后对大小关系进行 0/1编码。本发明将对不同空间领域，W及频率、方向的空间排列关系进行探索，最后找到最 佳编码增强方式。
[0042] (2)结合波段选择的S维Gabor特征筛选
[0043] 在总结高光谱人脸数据的特点的基础上，深入研究在空域和谱域同时提取信号变 化特性的方法，从理论上分析比较运些方法的优缺点并对其进行优化，重点关注S维Gabor 小波在该问题上的应用。
[0044] 探索通过波段选择去除高光谱数据中冗余信息的方法，量化比较各种选择方案的 优劣，尝试引入生物启发算法对选择方案进行优化，减少其对特征提取的负面影响；比较空 域、谱域联合分析和波段选择的不同结合方案，使得在重要波段上提取的特征能够获得最 佳的识别性能。
[0045] (3)基于多视角张量的S维Gabor特征融合
[0046] 在获得各个小波的增强编码特征之后，需要对运些特征的融合方法进行深入研 究。现有的融合方法是将各个Gabor编码特征转换成为一维向量，然后连接成为一个长特征 向量。对于大小为100 X 100 X 100的S维数据，52个小波抽取的特征维度可高达52，000， 000。如此高维特征会给后续分类带来巨大的计算量。除了采用特征选择方法减少特征维 度，通过张量表示可W在融合运些特征的同时，保持不同Gabor特征结构关系。
[0047] 本发明将首先研究目前最新的张量理论与方法，比较包括CP分解在内的不同张量 分解方法，重点关注多视角张量数据构造W及基于张量CP分解和因子分解的张量模型，结 合Gabor特征特性设计Gabor张量表示，发掘运些多视角张量之间的结构关系，构造基于核 函数的相似性度量W及相应支持向量机分类算法。
[004引探索基于张量的S维Gabor特征表示方法，将多个Gabor响应或编码结果表示成S 阶张量，根据Gabor张量的多视角特性设计相应张量分解算法，同时学习运些S阶张量的共 同因子和各个视角上的潜在因子，获得各Gabor特征之间的结构关系，最后为每一个Gabor 视角上提取的因子矩阵设计相应的结构保距核函数，通过支持张量机实现分类。
[0049]本发明需要解决的关键问题及方案
[0050] (I)基于方向、频率和S维空间的Gabor特征编码问题
[0051] 在将S维数据和不同的S维Gabor小波卷积获得不同响应之后，如何在方向、频率 和S维空间内对运些响应进行编码，提取运些响应之间的联系，从而增强S维Gabor特征对 噪声的鲁棒性，是本发明所要解决的第一个关键问题。
[0052] 给定高光谱数据V(x，y，A)，S维小波分析通过小波函数口巧日V(x，y，A)的内积 来同时获得信号的空域和频域信息。
[0053] 假定中屯、为(Xc，yc，、）的对称小波函数口在立维空间和频域的覆盖范围分别为Xc ± A x、yc± A y、Ac± A AW及fc± A f，S维小波分析则可W提取信号在[Xe- A x,Xc+ A X] X [ye-A y ,yc+A y] X [Ac-A A,、+A A] X [fc-A f ,fc+A f ]范围内的变化信息。在众多小波函 数中，S维Gabor小波可W使A X、A y、A AW及A f的乘积达到最小，即达到分析信号的最佳 分辨率。S维Gabor函数主要由S维高斯函数和实虚部分别为正余弦的复函数乘积组成，具 体定义如下：
[0化4] (1)
[00对公式中S为归一化常量；（Xe,yc，、)为Gabor小波的中屯、位置，类似于神经元感受野 中屯、；Ox、Oy W及OA决定小波在X、y和A方向的高斯包络大小；参数f为小波的中屯、频率，4为 该频率在S维频域与W轴的夹角，0为该频率在S维频域的U-V平面内投影与U轴的夹角（参 见附图1);R为将高斯包络主轴和中屯、频率振荡方向调整为一致的旋转矩阵。
[0化6] 在采用Gabor小波抽取特征时，需要预先设计多个不同中屯、频率fi W及方向(0,，A) 的小波
，然后将其和=维图像数据V(x，y，A)进行 卷积操作。卷积结果
t表了在=维空间（x，y，A)、频率为。、方向 为(0,，私)的信号变化信息。附图2显示了多个S维Gabor小波在频域的示意图，运些小波覆盖 了频域中多个不同频率和方向。
[0化7] 给定N个S维Gabor小波= 2,…N)，可W通过卷积获得N个同样大小的S维 数据F @ 7, O本发明从S维空间、频域W及方向域S个角度对运N个S维Gabor特 征进行编码，W获得更具辨识力和抗干扰性能的特征。空间编码的主要思想是将中屯、像素 的响应与相邻像素的响应进行差异性比较，然后对大小关系进行0/1编码，最后再转换成整 数，通过区域直方图统计对Gabor特征进行增强。在S维空间，中屯、像素的3X3X3邻域中有 6邻域，18邻域和26邻域=种。本发明将为运=种邻域关系分别设计相应的0/1编码方式，并 测试用不同方法增强后Gabor特征的鲁棒性，从而确定最佳邻域和编码方式。
[0058] 频率、方向编码的核屯、思想是将中屯、像素的Gabor响应和相邻频率(尺度）、方向的 响应进行差异性比较，然后对大小关系进行0/1编码，最后再转换成整数，通过区域直方图 统计对Gabor特征进行增强。本发明将相邻尺度、方向的Gabor响应排列成如附图3所示的虚 拟空间位置，然后通过该空间的邻域编码来实现特征增强。
[0059] 在对每个S维Gabor特征进行空间编码W及多个Gabor特征进行频率、方向编码 后，一幅S维图像Wx，y，A)可W获得K个编码后的特征图像gi(x，y，A)，…gk(x，y，A)，，，gK(x, y，入)。
[0060] (2)特征提取和波段选择的结合问题
[0061] 如何在不影响特征提取效果的前提下，去除冗余信息，从而减少后续处理的数据 量;如何使选择的波段更有利于特征提取，从而获得最佳识别性能。
[0062] 对于图像大小为64X64、波段数为30的多光谱图像，52个小波卷积后的数据量高 达64 X 64 X 30 X 52 = 6，389，760。对运些数据进行分类处理需要大量的计算时间，因此很有 必要去除那些信噪比低、鉴别性差的波段数据，从而减少冗余信息W及后续处理复杂度。本 发明通过波段和特征选择相接合的方法进行，主要步骤如下：
[0063] Stepl:设计L个不同中屯、频率和方向的S维Gabor小波{ Wi, 1 = 1，..，L}，将其和B 个波段的多光谱数据Wx，y，A)卷积，获得BXL个卷积图像：|Mi(x，y，Ai)，…化(x，y，Ai)…Mi (x,y,Ab)...]V[L(x,y,^b)...Mi(x,y,AB)...ML(x,y,AB)};
[0064] Step2:依次提取每个波段Ab的L个卷积图像{Mi(x，y，Ab)…化^，7，人〇}，通过下采 样等手段提取维度较低的特征，设计评价函数进行排序，找出最优波段。将该波段从候选集 合移出，加入到已选集合；
[0065] Step3:从候选集合选取一个波段加入到已选集合，使得已选集合的评价函数值得 到最大的增加；
[0066] Step4:重复上述过程Stepl-3,直到增加量小于设定阔值或者迭代次数达到预先 设定值。
[0067] 对波段评价函数的直接选择是识别准确率，但是每次迭代对每个个体进行准确率 测试会耗费大量计算时间，比较折中的选择是采用诸如Fisher准则之类的可分性指标。类 间距和类内距的比值越大，特征越好。本发明将对各种评价函数进行测试，最终找到合理的 方案，考虑计算复杂度，实现高效、准确的多光谱目标识别方法。
[0068] (3)多视角Gabor张量的构造与学习问题
[0069] 在获取经过选择的S维Gabor特征之后，运些特征包含了高光谱图像不同方向、频 率的信号变化信息，如何将运些多视角特征通过张量表示，在后续张量分解W及分类过程 中提取运些Gabor特征在方向、频率和S维空间内的结构关系，W及如何构造核函数实现基 于支持向量机的分类，是另一个关键问题。
[0070] 波段选择后提取的S维Gabor特征图像可表示为K个S阶张量X= {xi，X2…XkK其 中义j e沢/'xjzws。传统的张量学习方法往往将其整合为一个更大的四阶张量进行张量分解学 习共同的潜在因子，或者针对多个=阶张量进行独立的张量分解学习各个视角上的潜在因 子，然后进行后续的识别过程。本发明针对该Gabor张量的多视角特性W及编码数据本身的 非负性和冗余性，设计稀疏非负多视角张量分解算法，同时学习 K个=阶张量的共同因子和 各个视角上的潜在因子(如附图4所示）。
[0071] -般地，一个S阶张量Xk的CP分解优化模型可W写为：
[007^
口）
[007引为简便起见，用1](1￥2)1](3)表示^"!.。。"^。"^，其中吟为妒喻第件列向量，^， M- 通常称为第i模态的因子矩阵。基于此，本发明拟优化的目标函数如下：
[0074]
[0075]
[0076] 其中，巧W为第i个视角下第j个模态上的因子矩阵，1巧]表示第一模态上的因 子矩阵，护为共同因子矩阵(运里也可W根据实际情况考虑其他模态）；M I kl为iy范数， 用于特征的稀疏化学习。目标函数(3)可首先根据《2,,范数的分析方法，将其转换为矩阵迹 的表达形式，进而利用交替最小二乘（alternating least squares)和半二次法化alf- qua化atic way)结合进行优化求解。即可迭代优化如下目标函数
[0077]
[007引
[0079] 其中矿和势^'>分别为与护和巧"相关的辅助变量。
[0080] 通过求解上述公式可得到因子矩阵护和巧A，:即为多视角张量分解提取的潜在因 子矩阵。需要说明的是，每一个视角上提取的因子矩阵仍然是CP分解形式，可有效地近似表 示原始视角张量。在此基础上，可W根据基于CP分解构造的张量结构保距核函数W及多核 学习理论，构造多视角张量核函数，进而结合第二类张量机进行分类学习。
[0081] 相对可见光二维图像的广泛应用，高光谱图像识别分类无论在基础理论、还是应 用领域的研究都还处于起步阶段。针对基于S维Gabor小波的空间-光谱域联合变化特征提 取W及相关的=个关键问题，本发明的特色与创新之处在于：
[0082] (1)高光谱图像识别本身是模式识别领域一个很新的方向，通过=维空间域W及 频域联合分析，对=维和高光谱图像进行并行特征提取，然后对运些特征进行增强和融合， 可W充分利用高维数据的丰富信息W及运些特征之间的结构关系，大大减少噪声对准确率 的影响，本发明技术内容是该领域理论和应用上的全新尝试；
[0083] (2)目前S维Gabor特征主要基于幅值响应W及区域统计为主，业界并没有对各种 相位、方向、频域W及S维空间编码方法进行深入研究。本发明从提取各Gabor特征在方向 与频率空间结构关系出发，对不同编码方法进行研究，可大大加强高维Gabor特征的区分能 力；
[0084] (3)在联合提取空间-光谱域和频域的信号变化特征同时，本发明开拓性地将智能 算法引入波段选择，从而可W在保证识别准确率的前提下减少冗余信息，提高识别效率，具 有较广泛的应用和推广前景；
[0085] (4)传统的向量化和简单拼接方式极其容易造成S维Gabor特征的维数灾难问题， 从而给后续分类带来非常大的计算复杂度。本发明从张量角度出发，国内外首次提出基于 多视角的高阶Gabor张量构造、学习W及分类方法，可W在提取不同Gabor特征结构关系的 同时，对运些特征进行选择和融合，从而进一步提高高维图像识别系统的准确性和鲁棒性。
[0086] W上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定 本发明的具体实施只局限于运些说明。对于本发明所属技术领域的普通技术人员来说，在 不脱离本发明构思的前提下，还可W做出若干简单推演或替换，都应当视为属于本发明的 保护犯i围。
【主权项】
1. 一种高光谱图像特征的并行提取与分类方法，其特征在于，所述方法包括以下步骤： A. 获取高光图图像； B. 对所述高光图图像进行基于三维Gabor小波的空间-光谱域联合变化的并行特征提 取;具体为： 给定三维高光谱数据V(x，y，λ)，三维小波分析通过三维Gabor小波函数ΨΜ,φ(χ，γ，λ) 和V(x，y，A)的内积来同时获得信号的空域和频域信息，其中，参数f为小波的中心频率，φ 为该频率在三维频域与w轴的夹角，Θ为该频率在三维频域的u-v平面内投影与u轴的夹角； 在采用Gabor小波抽取特征时，需要预先设计多个不同中心频率fi以及方向(0,,?)的小波然后将其和VU，y，A)进行卷积操作，卷积结果 心,#，>^) = (?，0平,,>,>'3)代表了在三维空间（1，7，\)、频率为匕方向为吆说)的信号变化 信息；给定N个三维Gabor小波Ψη(η=1，2，··· N)，通过卷积获得N个同样大小的三维数据 r ? Ψλ(α乃2);从三维空间、频域以及方向域三个角度对所述N个三维Gabor特征进行空 间编码以及频率、方向编码后，V(x，y，A)获得K个编码后的特征图像 gHx.yA)，.Hgk(XjJ)…gK(x，y，A); C. 结合波段选择对所述步骤B提取的三维Gabor特征进行筛选，通过所述波段选择去除 高光谱数据中冗余信息； D. 构造并学习基于多视角的高阶Gabor张量，实现三维Gabor特征融合； E. 构造核函数实现基于支持向量机的图像分类。2. 根据权利要求1所述的方法，其特征在于:所述空间编码的主要思想是将中心像素的 响应与相邻像素的响应进行差异性比较，然后对大小关系进行0/1编码，最后再转换成整 数，通过区域直方图统计对Gabor特征进行增强。3. 根据权利要求1所述的方法，其特征在于:所述频率、方向编码的核心思想是将中心 像素的Gabor响应和相邻频率、方向的响应进行差异性比较，然后对大小关系进行0/1编码， 最后再转换成整数，通过区域直方图统计对Gabor特征进行增强。4. 根据权利要求1所述的方法，其特征在于:所述步骤C具体包括： Step 1:设计L个不同中心频率和方向的三维Gabor小波{Ψι，1 = 1，..，L}，将其和B个波 段的多光谱数据￥(1，7，人)卷积，获得8\1^个卷积图像：{11(1，7，人1)，"低(1，7，人1；^"]/[1(1，7， Xb)...ML(x，y，Xb)...Mi(x，y，λΒ)···ΜιΧχ，γ，λβ)}; Step2:依次提取每个波段AdtlL个卷积图像{Μι (χ，y，Ab)…Ml(χ，y，Ab)}，提取维度较低 的特征，设计评价函数进行排序，找出最优波段;将该最优波段从候选集合移出，加入到已 选集合； Step3:从候选集合选取一个波段加入到已选集合，使得已选集合的评价函数值得到最 大的增加； Step4:重复上述过程Stepl-3,直到增加量小于设定阈值或者迭代次数达到预先设定 值。5. 根据权利要求1所述的方法，其特征在于:所述步骤D具体为:波段选择后提取的三维 Gabor特征图像表示为K个三阶张量X= {xi,X2…Χκ}，其中Z e货/1></2><13 ;利用稀疏非负多视角 张量分解算法，同时学习 K个三阶张量的共同因子和各个视角上的潜在因子:一个三阶张量 Xk的CP分解优化模型写为为简便起见，用U(1)u(2)u(34丨第j个列向量，u(i)通常 称为第i模态的因子矩阵，拟优化的目标函数如下：其中，?/Ρ为第i个视角下第j个模态上的因子矩阵，[V Pm]表示第一模态上的因子矩 阵，#为共同因子矩阵(这里也可以根据实际情况考虑其他模态）；I I I I24SA4范数，用于 特征的稀疏化学习；所述目标函数首先根据:范数的分析方法，将其转换为矩阵迹的表达 形式，进而利用交替最小二乘和半二次法结合进行优化求解得到因子矩阵If和?/f，即为多 视角张量分解提取的潜在因子矩阵。
【文档编号】G06K9/62GK106022391SQ201610380345
【公开日】2016年10月12日
【申请日】2016年5月31日
【发明人】王岢, 张海军, 李旭涛, 叶允明, 徐晓飞
【申请人】哈尔滨工业大学深圳研究生院