基于光谱抽样直方图的超光谱降维匹配方法及系统的制作方法
【专利摘要】本发明提供一种基于光谱抽样直方图的超光谱降维匹配方法及系统，包括对待匹配光谱和光谱库中的所有光谱分别进行归一化处理，分别获取归一化后的待匹配光谱和光谱库中所有光谱的抽样直方图，计算待匹配光谱的抽样直方图与光谱库中所有光谱的抽样直方图的欧氏距离，在光谱库中选取与待匹配光谱抽样直方图欧氏距离最小的一条光谱作为匹配对象。本发明通过对归一化后的光谱使用等间距的窄带进行抽样，从而获得维数远小于原始光谱的抽样直方图，完成了光谱的降维，然后使用降维后的抽样直方图代替原始光谱进行匹配，显著降低了后续匹配时的运算量，同时在抽样时利用分段提取的方法保留了光谱图中的相对位置信息，提高了匹配的精度。
【专利说明】基于光谱抽样直方图的超光谱降维匹配方法及系统
【技术领域】
[0001]本发明涉及超光谱降维匹配【技术领域】，具体地说，本发明涉及一种基于光谱抽样直方图的超光谱降维匹配方法及系统。
【背景技术】
[0002]光谱遥感仪器正向着更高的光谱分辨率、空间分辨率以及时间分辨率的方向发展，测量得到的光谱数据维数可达数千维甚至数万维。在得到超光谱数据后，最常见也是最基础的需求是鉴别物质的种类，而鉴别种类就需要进行光谱的匹配。简单的说，光谱匹配是通过计算逐一比较光谱数据库中的所有参考光谱与待测光谱，选择最近似的一条光谱作为待测光谱的匹配结果。匹配结果可用于确定待测光谱的种类，或者作为对照信息供进一步分析时参考。
[0003]然而，传统的光谱匹配算法很难适应超光谱数据维数巨大这一特性，计算复杂度大而且耗时长，无法达到实时性的要求。解决这一问题的关键在于使用适当的光谱降维算法对超光谱数据进行降维处理，即通过大幅降低光谱维数来减小运算量，提高实时性，同时保证光谱匹配的准确性尽量不受影响。但从研究现状来看，还没有专门用于超光谱数据的降维匹配算法，目前的光谱降维算法都是为几百维的高光谱数据而设计的，若直接用于超光谱数据降维然后再匹配，其性能不佳。
【发明内容】
[0004]为克服相应技术缺陷，本发明提出了一种基于光谱抽样直方图的超光谱降维匹配技术方案。
[0005]本发明技术方案提供一种基于光谱抽样直方图的超光谱降维匹配方法，包括以下步骤:
[0006]步骤I，光谱归一化，包括对待匹配光谱和光谱库中的所有光谱分别进行归一化处理，归一化公式如下，
[0007]Pn= (P-P min) / ( P max- P min)
[0008]其中，Pn为归一化后的光谱辐射值，P为原始光谱的辐射值，Pmax为该条原始光谱中辐射值的最大值，Pmin为该条原始光谱中辐射值的最小值；
[0009]步骤2，获取抽样直方图，包括分别获取归一化后的待匹配光谱和光谱库中所有光谱的抽样直方图，实现方式为对归一化后的待匹配光谱和归一化后的光谱库中任一光谱，按照相应归一化后的光谱图，分别执行以下操作，
[0010]将归一化后的光谱图在整个超光谱波段范围平均划分为s段，s为预设参数；
[0011]在归一化后的光谱图上作M条互相平行且等间隔的宽度为2d的水平窄带，M和d为预设参数，则M条窄带在s个不同波段内与归一化后的光谱图中光谱曲线相交的次数组成抽样直方图，用(sXM)维向量{ α η, α 12，…α 1Μ,…，a sl,...，a sM}表不,其中a 为第i个波段中第j条窄带与光谱曲线相交的次数，i = I，2,..., s, j = 1,2,..., M ；
[0012]步骤3，计算欧氏距离，包括计算待匹配光谱的抽样直方图与光谱库中所有光谱的抽样直方图的欧氏距离，设光谱库中有N条光谱，则分别计算得到待测光谱与光谱库中每条光谱抽样直方图的欧氏距离EDk，k = 1，…，N，实现方式为对待匹配光谱的抽样直方图与光谱库中任一光谱的抽样直方图计算欧氏距离如下，
[0013]设A' = {α ' η, α ' 12，…α , 1Μ，…α , sl，…α , sM}为待测光谱的抽样直方图，a' u为待测光谱相应第i个波段中第j条窄带与光谱曲线相交的次数，A"={α " η, α " 12，…α " 1Μ,…α " sl,…α " sM}为光谱库中任意一条光谱的抽样直方图，α " 为光谱库中任意一条光谱相应第i个波段中第j条窄带与光谱曲线相交的次数，
[0014]计算欧氏距离的公式为，
[0015]ED{A',A") = Jih^-aD2
Y /-1/=I
[0016]步骤4，光谱匹配，包括在光谱库中选取与待匹配光谱抽样直方图欧氏距离最小的一条光谱作为匹配对象。
[0017]本发明还相应提供一种基于光谱抽样直方图的超光谱降维匹配系统，包括以下模块:
[0018]光谱归一化模块，用于对待匹配光谱和光谱库中的所有光谱分别进行归一化处理，归一化公式如下，
[0019]Pn= (P-P min) / ( P max- P min)
[0020]其中，Pn为归一化后的光谱辐射值，P为原始光谱的辐射值，Pmax为该条原始光谱中辐射值的最大值，Pmin为该条原始光谱中辐射值的最小值；
[0021]抽样直方图获取模块，用于分别获取归一化后的待匹配光谱和光谱库中所有光谱的抽样直方图，实现方式为对归一化后的待匹配光谱和归一化后的光谱库中任一光谱，按照相应归一化后的光谱图，分别执行以下操作，
[0022]将归一化后的光谱图在整个超光谱波段范围平均划分为s段，s为预设参数；
[0023]在归一化后的光谱图上作M条互相平行且等间隔的宽度为2d的水平窄带，M和d为预设参数，则M条窄带在s个不同波段内与归一化后的光谱图中光谱曲线相交的次数组成抽样直方图，用(sXM)维向量{ α η, α 12，…α 1Μ,…，a sl,...，a sM}表不,其中a 为第i个波段中第j条窄带与光谱曲线相交的次数，i = I，2,..., s, j = 1,2,..., M ；
[0024]欧氏距离计算模块，用于计算待匹配光谱的抽样直方图与光谱库中所有光谱的抽样直方图的欧氏距离，设光谱库中有N条光谱，则分别计算得到待测光谱与光谱库中每条光谱抽样直方图的欧氏距离EDk，k = 1，…，N，实现方式为对待匹配光谱的抽样直方图与光谱库中任一光谱的抽样直方图计算欧氏距离如下，
[0025]设A' = {α ' η, α ' 12，…α , 1Μ,…α , sl,...α , sM}为待测光谱的抽样直方图，a' u为待测光谱相应第i个波段中第j条窄带与光谱曲线相交的次数，A"={α " η, α " 12，…α " 1Μ,…α " sl,…α " sM}为光谱库中任意一条光谱的抽样直方图，α " 为光谱库中任意一条光谱相应第i个波段中第j条窄带与光谱曲线相交的次数，
[0026]计算欧氏距离的公式为，
[0027]ED(ArfAn) = JziK-aP2
V /=1 j=i
[0028]光谱匹配模块，用于在光谱库中选取与待匹配光谱抽样直方图欧氏距离最小的一条光谱作为匹配对象。
[0029]本发明分析超光谱的抽样直方图信息，通过对归一化后的光谱使用等间距的窄带进行抽样，从而获得维数远小于原始光谱的抽样直方图，完成了光谱的降维，显著降低了后续匹配时的运算量，同时在抽样时利用分段提取的方法保留了光谱图中的相对位置信息，提高了匹配的精度。本发明的技术方案可以抑制不同级别的噪声对光谱图抽样的影响，具有良好的抗噪性和鲁棒性，相对传统方法可以显著降低匹配运算量，实时性好，具有良好的抗噪性和鲁棒性，匹配精度高。
【专利附图】
【附图说明】
[0030]图1是本发明实施例的流程图。
[0031]图2是本发明实施例的分波段抽样示意图。
【具体实施方式】
[0032]下面结合附图和实施例对本发明进行进一步的说明。
[0033]参照附图1，本发明主要由4个步骤组成:光谱归一化，获取光谱的抽样直方图，计算欧氏距离，取最小欧氏距离完成匹配。实施例选取含有1432种物质光谱的光谱库，每种物质在光谱库中只有一条光谱数据。光谱分辨率为Λ σ = 0.1cnT1，波长范围为2-14 μ m，对应波数范围为5000-714(^'在此光谱范围内共有N = 42861个采样点，即原始光谱为42861 维。
[0034]具体实施时，本发明技术方案可采用计算机软件技术实现自动运行流程。实施例执行步骤如下:
[0035](I)光谱归一化步骤包括，分别对待匹配光谱和光谱库中的所有光谱进行归一化处理，归一化公式如下:
[0036]Pn= (P-P min) / ( P max- P min)
[0037]其中，Pn为归一化后的光谱辐射值，P为原始光谱的辐射值，Pmax为该条原始光谱中辐射值的最大值，P min为该条原始光谱中辐射值的最小值。
[0038]将待匹配光谱和光谱库中的所有光谱分别作为原始光谱，按以上归一化公式进行处理，得到相应归一化后的待匹配光谱和光谱库中所有光谱。
[0039](2)获取抽样直方图步骤包括，分别获取归一化后的待匹配光谱和光谱库中所有光谱的抽样直方图。
[0040]实施例中，对归一化后的待匹配光谱和归一化后的光谱库中任一光谱，按照相应归一化后的光谱图，分别执行包括以下操作:
[0041](2.1)划分波段，包括将归一化后的光谱图在整个超光谱波段范围平均划分为s段，S为预设参数，具体实施时本领域技术人员可自行预设的s取值，本实施例中S取3 ;
[0042](2.2)在每个波段分别抽样，联合起来组成光谱的抽样直方图:在归一化后的光谱图上作M条互相平行且等间隔的宽度为2d的水平窄带，M和d为预设参数，具体实施时本领域技术人员可自行预设的M、d取值，则M条窄带在s个不同波段内与归一化后的光谱图中光谱曲线相交的次数组成光谱的抽样直方图，可以用(sXM)维向量{αη，α12^..α1Μ)-, asl，…，asM}来表示，其中a u为第i个波段中第j条窄带与光谱曲线相交的次数，i = l,2,..., s, j = I, 2,...,Μ。
[0043]实施例中d为0.001，M取3，则3条窄带在3个不同波段内与光谱曲线相交的次数组成光谱抽样直方图，可以用(sXM)=维向量{αη, α 12，α 13，…，α 31，α 32，α 33}来表不，其中a u为第i个波段中第j条窄带与光谱曲线相交的次数，i = 1,2,3, j = 1，2，3;附图2所示光谱的抽样直方图为{1，I, 4，I, I, O, O, O, 0}。
[0044](3)计算欧氏距离步骤包括，计算待匹配光谱的抽样直方图与光谱库中的所有光谱的抽样直方图的欧氏距离，计算公式如下:
[0045]ED{A\An)=\ ΣΣΚ-Ο2
Yi=l J=I
[0046]其中，设A' = {a ' n, a / 12,…a ' 1M，…a ' sl,…a ' sM}为待测光谱的抽样直方图，a ' u为待测光谱相应第i个波段中第j条窄带与光谱曲线相交的次数，i = 1，2，...，s，j = l，2，...，M;A" = {a" n, a " 12，…a " 1M，…a " sl，…a " sM}为光谱库中任意一条光谱的抽样直方图，a " u为光谱库中任意一条光谱相应第i个波段中第j条窄带与光谱曲线相交的次数，i = 1，2，…，s，j = l，2，…，Μ;设光谱库中有N条光谱，则分别计算得到待测光谱与光谱库中每条光谱抽样直方图的欧氏距离EDk，k = 1，…，N。
[0047]实施例的光谱库中有1432条光谱，则分别计算得到待测光谱与光谱库中每条光谱的抽样直方图的欧氏距离EDk，k = 1，…，1432。
[0048](4)光谱匹配步骤包括，在光谱库中选取与待匹配光谱抽样直方图欧氏距离最小的一条光谱作为匹配对象，可表达如下:
[0049]AmiJEDmin = min{EDk}
[0050]上式的意思为Amin为光谱数据库中与待匹配光谱的抽样直方图欧氏距离最小的抽样直方图，则Aniin所对应的光谱库中的原始光谱即为待匹配光谱的匹配对象。
[0051]本实施例中M和s均取3，仅为了便于理解示范，一般这两个参数均取20左右为佳，考虑运算时间也不宜太大。
[0052]本发明还相应提供一种基于光谱抽样直方图的超光谱降维匹配系统，包括以下模块:
[0053]光谱归一化模块，用于对待匹配光谱和光谱库中的所有光谱分别进行归一化处理，归一化公式如下，
[0054]Pn= (P-P min) / ( P P min)
[0055]其中，Pn为归一化后的光谱辐射值，P为原始光谱的辐射值，Pmax为该条原始光谱中辐射值的最大值，Pmin为该条原始光谱中辐射值的最小值；
[0056]抽样直方图获取模块，用于分别获取归一化后的待匹配光谱和光谱库中所有光谱的抽样直方图，实现方式为对归一化后的待匹配光谱和归一化后的光谱库中任一光谱，按照相应归一化后的光谱图，分别执行以下操作，
[0057]将归一化后的光谱图在整个超光谱波段范围平均划分为s段，s为预设参数；
[0058]在归一化后的光谱图上作M条互相平行且等间隔的宽度为2d的水平窄带，M和d为预设参数，则M条窄带在s个不同波段内与归一化后的光谱图中光谱曲线相交的次数组成抽样直方图，用(sXM)维向量{ α η, α 12，…α 1Μ,…，a sl,...，a sM}表不,其中a 为第i个波段中第j条窄带与光谱曲线相交的次数，i = I，2,..., s, j = 1,2,..., M ；
[0059]欧氏距离计算模块，用于计算待匹配光谱的抽样直方图与光谱库中所有光谱的抽样直方图的欧氏距离，设光谱库中有N条光谱，则分别计算得到待测光谱与光谱库中每条光谱抽样直方图的欧氏距离EDk，k = 1，…，N，实现方式为对待匹配光谱的抽样直方图与光谱库中任一光谱的抽样直方图计算欧氏距离如下，
[0060]设A' = {α ' η, α ' 12，…α , 1Μ，…α , sl，…α , sM}为待测光谱的抽样直方图，a' u为待测光谱相应第i个波段中第j条窄带与光谱曲线相交的次数，A"={α " η, α " 12，…α " 1Μ,…α " sl,…α " sM}为光谱库中任意一条光谱的抽样直方图，α " 为光谱库中任意一条光谱相应第i个波段中第j条窄带与光谱曲线相交的次数，
[0061]计算欧氏距离的公式为，
I s ~M ~Z ~
[0062]ED(A'，A") =
Vf=l J=I
[0063]光谱匹配模块，用于在光谱库中选取与待匹配光谱抽样直方图欧氏距离最小的一条光谱作为匹配对象。
[0064]应当理解的是，本说明书未详细阐述的部分均属于现有技术。
[0065]应当理解的是，上述针对实施例的描述较为详细，并不能因此而认为是对本发明专利保护范围的限制，本领域的普通技术人员在本发明的启示下，在不脱离本发明权利要求所保护的范围情况下，还可以做出替换或变形，均落入本发明的保护范围之内，本发明的请求保护范围应以所附权利要求为准。
[0066]总体来说，本发明所提出的基于抽样直方图的超光谱降维匹配方法，在保证了很高匹配成功率的前提下，极大降低了匹配时间。此外，一些具体参数的取值可根据不同情况灵活选择，在实时性和准确性间取得平衡，满足实际应用的需求。
【权利要求】
1.一种基于光谱抽样直方图的超光谱降维匹配方法，其特征在于，包括以下步骤: 步骤1，光谱归一化，包括对待匹配光谱和光谱库中的所有光谱分别进行归一化处理，归一化公式如下，
Pn= (P-P min) / ( P P min) 其中，Pn为归一化后的光谱辐射值，P为原始光谱的辐射值，Pmax为该条原始光谱中辐射值的最大值，Pmin为该条原始光谱中辐射值的最小值； 步骤2，获取抽样直方图，包括分别获取归一化后的待匹配光谱和光谱库中所有光谱的抽样直方图，实现方式为对归一化后的待匹配光谱和归一化后的光谱库中任一光谱，按照相应归一化后的光谱图，分别执行以下操作， 将归一化后的光谱图在整个超光谱波段范围平均划分为S段，S为预设参数； 在归一化后的光谱图上作M条互相平行且等间隔的宽度为2d的水平窄带，M和d为预设参数，则M条窄带在S个不同波段内与归一化后的光谱图中光谱曲线相交的次数组成抽样直方图，用(sXM)维向量{ α η, α 12，…a m,…，a sl,...，a sM}表不,其中a Ij为第i个波段中第j条窄带与光谱曲线相交的次数，i = I，2,..., s, j = 1,2,..., M ； 步骤3，计算欧氏距离，包括计算待匹配光谱的抽样直方图与光谱库中所有光谱的抽样直方图的欧氏距离，设光谱库中有N条光谱，则分别计算得到待测光谱与光谱库中每条光谱抽样直方图的欧氏距离EDk，k= 1，…，N，实现方式为对待匹配光谱的抽样直方图与光谱库中任一光谱的抽样直方图计算欧氏距离如下， 设A' = {α ' η, α ' 12，…α ' 1Μ，…α ' sl，…α ' sM}为待测光谱的抽样直方图，a' u为待测光谱相应第i个波段中第j条窄带与光谱曲线相交的次数，A"={α " η, α " 12，…α " 1Μ,…α " sl,…α " sM}为光谱库中任意一条光谱的抽样直方图，α " 为光谱库中任意一条光谱相应第i个波段中第j条窄带与光谱曲线相交的次数，计算欧氏距离的公式为，
ED(A',A") = 步骤4，光谱匹配，包括在光谱库中选取与待匹配光谱抽样直方图欧氏距离最小的一条光谱作为匹配对象。
2.一种基于光谱抽样直方图的超光谱降维匹配系统，其特征在于，包括以下模块: 光谱归一化模块，用于对待匹配光谱和光谱库中的所有光谱分别进行归一化处理，归一化公式如下，
Pn= (P-P min) / ( P P min) 其中，Pn为归一化后的光谱辐射值，P为原始光谱的辐射值，Pmax为该条原始光谱中辐射值的最大值，Pmin为该条原始光谱中辐射值的最小值； 抽样直方图获取模块，用于分别获取归一化后的待匹配光谱和光谱库中所有光谱的抽样直方图，实现方式为对归一化后的待匹配光谱和归一化后的光谱库中任一光谱，按照相应归一化后的光谱图，分别执行以下操作， 将归一化后的光谱图在整个超光谱波段范围平均划分为S段，S为预设参数； 在归一化后的光谱图上作M条互相平行且等间隔的宽度为2d的水平窄带，M和d为预设参数，则M条窄带在s个不同波段内与归一化后的光谱图中光谱曲线相交的次数组成抽样直方图，用(sXM)维向量{ α η, α 12，…a m,…，a sl,...，a sM}表不,其中a Ij为第i个波段中第j条窄带与光谱曲线相交的次数，i = I，2,..., s, j = 1,2,..., M ； 欧氏距离计算模块，用于计算待匹配光谱的抽样直方图与光谱库中所有光谱的抽样直方图的欧氏距离，设光谱库中有N条光谱，则分别计算得到待测光谱与光谱库中每条光谱抽样直方图的欧氏距离EDk，k = 1，…，N，实现方式为对待匹配光谱的抽样直方图与光谱库中任一光谱的抽样直方图计算欧氏距离如下， 设A' = {α ' η, α ' 12，…α ' 1Μ，…α ' sl，…α ' sM}为待测光谱的抽样直方图，a' u为待测光谱相应第i个波段中第j条窄带与光谱曲线相交的次数，A"={α " η, α " 12，…α " 1Μ,…α " sl,…α " sM}为光谱库中任意一条光谱的抽样直方图，α " 为光谱库中任意一条光谱相应第i个波段中第j条窄带与光谱曲线相交的次数，计算欧氏距离的公式为， ^^?=^ΣΣ(^/-^)2 光谱匹配模块，用于在光谱库中选取与待匹配光谱抽样直方图欧氏距离最小的一条光谱作为匹配对象。
【文档编号】G06T7/00GK104268896SQ201410584276
【公开日】2015年1月7日 申请日期:2014年10月27日 优先权日:2014年10月27日
【发明者】黄珺, 马佳义 申请人:武汉大学