主权项 |
一种基于光谱抽样直方图的超光谱降维匹配方法,其特征在于,包括以下步骤:步骤1,光谱归一化,包括对待匹配光谱和光谱库中的所有光谱分别进行归一化处理,归一化公式如下,ρ<sub>n</sub>=(ρ‑ρ<sub>min</sub>)/(ρ<sub>max</sub>‑ρ<sub>min</sub>)其中,ρ<sub>n</sub>为归一化后的光谱辐射值,ρ为原始光谱的辐射值,ρ<sub>max</sub>为该条原始光谱中辐射值的最大值,ρ<sub>min</sub>为该条原始光谱中辐射值的最小值;步骤2,获取抽样直方图,包括分别获取归一化后的待匹配光谱和光谱库中所有光谱的抽样直方图,实现方式为对归一化后的待匹配光谱和归一化后的光谱库中任一光谱,按照相应归一化后的光谱图,分别执行以下操作,将归一化后的光谱图在整个超光谱波段范围平均划分为s段,s为预设参数;在归一化后的光谱图上作M条互相平行且等间隔的宽度为2d的水平窄带,M和d为预设参数,则M条窄带在s个不同波段内与归一化后的光谱图中光谱曲线相交的次数组成抽样直方图,用s×M维向量{α<sub>11</sub>,α<sub>12</sub>,…α<sub>1M</sub>,…,α<sub>s1</sub>,…,α<sub>sM</sub>}表示,其中α<sub>ij</sub>为第i个波段中第j条窄带与光谱曲线相交的次数,i=1,2,…,s,j=1,2,…,M;步骤3,计算欧氏距离,包括计算待匹配光谱的抽样直方图与光谱库中所有光谱的抽样直方图的欧氏距离,设光谱库中有N条光谱,则分别计算得到待匹配光谱与光谱库中每条光谱抽样直方图的欧氏距离ED<sub>k</sub>,k=1,…,N,实现方式为对待匹配光谱的抽样直方图与光谱库中任一光谱的抽样直方图计算欧氏距离如下,设A′={α′<sub>11</sub>,α′<sub>12</sub>,…α′<sub>1M</sub>,…α′<sub>s1</sub>,…α′<sub>sM</sub>}为待匹配光谱的抽样直方图,α′<sub>ij</sub>为待匹配光谱相应第i个波段中第j条窄带与光谱曲线相交的次数,A″={α″<sub>11</sub>,α″<sub>12</sub>,…α″<sub>1M</sub>,…α″<sub>s1</sub>,…α″<sub>sM</sub>}为光谱库中任意一条光谱的抽样直方图,α″<sub>ij</sub>为光谱库中任意一条光谱相应第i个波段中第j条窄带与光谱曲线相交的次数,计算欧氏距离的公式为,<maths num="0001"><math><![CDATA[<mrow><mi>E</mi><mi>D</mi><mrow><mo>(</mo><msup><mi>A</mi><mo>′</mo></msup><mo>,</mo><msup><mi>A</mi><mrow><mo>′</mo><mo>′</mo></mrow></msup><mo>)</mo></mrow><mo>=</mo><msqrt><mrow><munderover><mo>Σ</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>s</mi></munderover><munderover><mo>Σ</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><msubsup><mi>α</mi><mrow><mi>i</mi><mi>j</mi></mrow><mo>′</mo></msubsup><mo>-</mo><msubsup><mi>α</mi><mrow><mi>i</mi><mi>j</mi></mrow><mrow><mo>′</mo><mo>′</mo></mrow></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt></mrow>]]></math><img file="FDA0001140248430000011.GIF" wi="790" he="167" /></maths>步骤4,光谱匹配,包括在光谱库中选取与待匹配光谱抽样直方图欧氏距离最小的一条光谱作为匹配对象。 |