[发明专利]一种基于光谱指纹特征的弱小目标精确跟踪方法

权利要求书

1.一种基于光谱指纹特征的弱小目标跟踪方法，包括以下三个步骤：

步骤一：目标光谱稳定吸收特征描述

通过画幅式成像光谱仪获得前景高光谱图像，并对高光谱图像进行数据预处理，用以消除设备噪声；为了准确高效地描述目标的光谱特性，使用一系列光谱统计数据来准确的反应每个像素点上所包含的数百个波段的值，其具体的步骤如下：

(1)在特征库中引入目标二值化光谱曲线；

(2)构造一组特征向量来表述目标的强光谱特征，强光谱特征包括目标物体的固有稳定光谱吸收特征，目标与背景反差明显的光谱波段区间和近红外波段区间；

通过上述特征提取，获得一组代表目标稳定特征的数据，这一组特征数据称之为光谱指纹特征；在光谱跟踪算法中，把特征数据引入用光谱学特征描述目标形态的数学模型，通过扩展核函数的特征直方图维数，将这些特征值变换到直方图分布函数中，建立概率直方图特征模型；

步骤二：建立基于高光谱影像的运动目标跟踪数学模型

结合基本的可见光目标运动模型和光谱吸收特征，建立基于高光谱影像的运动目标跟踪数学模型；设m帧的高光谱图像序列为像素点x_i下的n维特征空间Hⁿ，包括步骤一中的光谱吸收特征：二值化曲线跳变位置h_b，近红外积分特征h_inf、大于均值的波峰数h_aMax、大于均值的波谷数h_aMin、小于均值的波峰数h_bMax、小于均值的波谷数h_bMin、上升趋势点数h_uGdt、下降趋势点数h_dGdt、波峰位置h_max和波谷h_min位置；对于u∈{1...m}，第u帧图像中的目标t_u可以表示为：

tu=1CsΣi=1nkK(||xi-x0h||2)δ(bs(f(xi)),u)---(1)]]>

其中，{x_i}_i＝1...n是像素点在图像中的位置，C_s是核直方图的归一化常数，h为跟踪窗口半径，函数K和δ的定义如上述式(1)所述；函数f(x_i)是一个特征空间函数，其作用在于获取像素点x_i下的n维特征空间Hⁿ中的10个特征值{h_b，h_inf，h_aMax，h_aMin，h_bMax，h_bMin，h_uGdt，h_dGdt，h_max，h_min}∈Hⁿ，并将这些特征值变换到[0，255]的区间上后返回到直方图分布函数b_s(x)中；根据获得的10个特征值，b_s(x)将该像素点分配到一个量化级数为λ＝16¹⁰的直方图；最后这个直方图被转换成概率直方图；

假设弱小目标的跟踪窗口为矩形区域，目标的状态模型可以用如下向量表示：

Sk=(CkT,WkT)T---(2)]]>

S_k表示目标第k个运动状态，C表示目标的位置属性，W表示目标的大小属性；

以跟踪窗口的概率直方图为特征，搜索最相似的候选区域S_k；假设S_k和目标参考状态S^*都具有λ个分量的核直方图，通过Bhattacharrya系数ρ(S^*，S_k)来度量两个离散概率分布的可分离误差，即

ρ(S*,Sk)=Σk=1mSkSk-1---(3)]]>

完成相似度计算后，通过启发式策略在候选空间中搜索最佳匹配目标，避免全局搜索；

上式在S_k(C_k)处泰勒展开，得到

ρ[S*(C‾ki),Sk)≈12Σk=1mS*C0Sk+λ2Σi=1nkωiK(||C‾ki-Cih||2)---(4)]]>

对于式(4)中的右边第二项，采用均值偏移算法迭代寻找其最大值，从而得到起始候选区域移向真实目标区域的平均位移：

m(C‾ki)=Σj=1mCjωjg(||C‾ki-Cjh||2)Σj=1mωg(||C‾ki-Cjh||2)-C‾ki---(5)]]>

其中，{C_j}_{j＝1，2，...，m}是状态在矩形区域像素坐标点；ω_j是在位置C_j处像素值的相似权重，是参考概率直方图与概率直方图二进制比率的平方根；g(.)是一个核函数；h为归一化后的窗口半径；在式(5)中，g(x)＝-K′(x)；该算法沿偏移向量的方向不断向核函数中心位置移动，这个过程已经被证明是收敛的，通过5次重复迭代基本可以达到准确跟踪目标的效果；

步骤三：针对高光谱图像的层次化粒子滤波跟踪

其具体流程如下：

(1)在跟踪区间内，选择N个随机不重合的目标点，N即为初始化的粒子个数；

(2)对每个目标点独立同分布采样得到初始粒子集

(3)计算k状态的粒子的数学期望：同时得到该粒子的权重：ωki=p(Ok|Ski);]]>

(4)将粒子集随机分成相等的两个子集，即

(5)对其中一个子集中的每个粒子，以精确跟踪模型进行迭代运算，完成目标的搜索；在每次迭代过程中，首先选择出加权值最大的粒子，然后使用该子集中的所有粒子，来计算该粒子状态中心点的平均位移，如公式(5)所示；

(6)经过一次迭代得到的子集与原有不变的子集合并，得到新的粒子集同时对两个子集的权值ω_j进行合并，使得

(7)重复(5)和(6)的迭代过程t次，t＝5；获得第k帧中车辆目标的状态向量估计值为：

S~k=Σj=1NSkjωj;]]>

综上三个步骤，最终算法输出目标的光谱指纹特征集、当前的即运动状态向量以及目标的运动轨迹。

2.根据权利要求1所述的一种基于光谱指纹特征的弱小目标跟踪方法，其特征在于：针对目标物体的固有稳定光谱吸收特征、目标与背景反差明显的光谱波段区间和近红外波段区间这三种特征，分别采用对应的提取方法进行：

1)物质的固有稳定光谱吸收特征反映了目标材质的特点，针对不同材质的目标，使用光谱数据库中对应的吸收波段知识对目标进行基本的提取与识别；

2)对目标与背景反差明显的光谱波段区间以及近红外波段区间的反射率数据进行变系数积分；而后，通过一个线性变换函数将数据放大，使得目标特征更加明显；

3)对于光谱数据中噪声相对较大，不适宜作为主参考特征源的数据，构造一组相应的弱光谱特征以强化参考特征数据；这些特征包括大于均值的波峰数、大于均值的波谷数、小于均值的波峰数、小于均值的波谷数、上升趋势点数、下降趋势点数、波峰位置和波谷位置。