[发明专利]一种面向精确打击的地面运动目标检测与识别方法

权利要求书

1.一种面向精确打击的地面运动目标检测与识别方法，其特征在于，包括以下步骤：

步骤一：联合时空信息提取目标候选区域：先利用时域信息获取运动图像，根据运动图像计算前向运动历史图和后向运动历史图，将前向运动历史图和后向运动历史图融合，按像素取最小值得到前后向运动历史图，作为输入图像；通过自适应阈值处理和连通域提取获得候选运动区域；再利用空域信息，计算所述候选运动区域的物质性分值，通过阈值处理获得目标候选区域；

其中，利用时域信息获取运动图像的具体方法包括：

(1.1)背景运动估计：采用KLT特征点跟踪器结合RANSAC算法估计背景运动，给定相邻两帧图像第τ帧图像I(τ)和第τ+1帧图像I(τ+1)，在第τ帧图像I(τ)上均匀取点，采用KLT特征点跟踪器提取第τ+1帧图像I(τ+1)上的匹配特征点，再利用RANSAC算法去除离群值，用得到的特征点对拟合8参数的平面投影变换，获得的单应性矩阵即为第τ帧图像到第τ+1帧图像的背景运动估计；

(1.2)采用帧差法，每N帧图像计算一幅运动图像；设第τ帧图像的背景运动补偿图像为I′(τm(N-1))，则根据所述背景运动估计对图像进行背景运动补偿：

(1.3)设运动图像为D(τ)，通过第τ帧图像与背景运动补偿图像I′(τm(N-1))的绝对差分得到D(τ)：

D(τ)＝|I(τ)-I′(τ m(N-1))| (4)

其中，“-”代表前向差分，得到前向运动图像D_F(τ)；“+”代表后向差分，得到后向运动图像D_B(τ)；

计算前后向运动历史图的具体方法包括：

(2.1)设第τ帧图像前向运动历史图为H_F(τ)，设第τ帧图像后向运动历史图为H_B(τ)；

(2.2)融合多层前向运动图像获得H_F(τ)，即将第τ帧图像前向运动历史图H_F(τ)表示成第τ-1帧图像前向运动历史图H_F(τ-1)与第τ帧图像前向运动图像D_F(τ)的函数：

其中，T为设定的阈值，d＝255/L为衰减项，L为前向运动历史图中包含的前向运动图像的有效层数；

(2.3)融合多层后向运动图像获得H_B(τ)，后向运动历史图包含目标的未来运动信息，表示成类似式(6)的递推形式：

不同的是，H_F(τ)由H_F(τ-1)递推一次即可得到，而H_B(τ)要由H_B(τ+L)递推L次才能得到；

(2.4)融合前向运动历史图H_F(τ)和后向运动历史图H_B(τ)，获得前后向运动历史图H_FB(τ)：

H_FB(τ)＝min(blur(H_F(τ))，blur(H_B(τ))) (7)

其中，blur是指平滑滤波器；

提取候选运动区域的具体方法包括：

(3.1)自适应阈值计算：设定一个灰度级阈值，当前后向运动历史图中灰度级的数量大于所述灰度级阈值时，采用大津法计算双阈值，选取较小的阈值作为自适应阈值；当图像中灰度级数量少于所述灰度级阈值时，采用大津法求单阈值，作为自适应阈值；

(3.2)候选运动区域提取：采用步骤(3.1)中得到的自适应阈值对所述输入图像进行阈值处理、二值化后，再进行腐蚀、膨胀的形态学处理，最后通过区域生长法提取连通域，得到候选运动区域；

获得目标候选区域的具体方法包括：

(4.1)设候选运动区域为R，候选运动区域的边框为候选框，以R的中心为中心，将候选框向外扩大θ倍，得到矩形环的外边框，以候选框作为矩形环的内边框，由此确定矩形环区域R_S(R，θ)，满足：

其中，外接矩形环参数θ为经验值，可以采用贝叶斯框架学习得到；

(4.2)计算中心环绕直方图，即分别计算候选运动区域的RGB/LAB颜色直方图H与矩形环区域的RGB/LAB颜色直方图H_S；然后计算中心环绕直方图的卡方距离，作为候选运动区域R在参数θ下的颜色对比度，设颜色对比度为CC(R，θ)：

CC(R，θ)＝x²(H(R)，H_S(R_S(R，θ))) (9)

其中，n_i为直方图中箱子的数量，H与H_S的箱子数量是相同的；i表示直方图中第i个箱子；Hⁱ、分别为直方图H与H_S第i个箱子的数值；

(4.3)RGB/LAB颜色直方图包括三个通道，对不同通道的颜色对比度采用最大值操作符进行区域物质性分值计算：

其中，c₁，c₂，c₃为所述三个通道；objectness为物质性分值；

(4.4)设定物质性分值阈值，滤除物质性分值小于所述物质性分值阈值的候选运动区域，得到目标候选区域；

步骤二：目标候选区域识别：通过局部低约束编码完成特征学习，并利用线性支持向量机实现对目标候选区域的识别，保留被识别为特定目标的候选区域，得到最终的检测结果；

所述目标候选区域识别分为离线训练和在线检测两个阶段，在离线训练阶段，构建典型底面运动目标和背景样本集，将样本映射到高维特征空间，训练得到一个多类分类器；在检测阶段，用与训练阶段相同的方法对候选区域进行特征提取，再对得到的高维特征进行分类，得到类标签和置信度；

其中，目标候选区域识别的具体方法包括：

(5.1)字典构造：

采用交替优化的方式训练字典，目标函数如式(12)：

其中，为字典训练过程中待优化的变量，即字典的估计和局部特征的重构系数的估计，局部特征的重构系数即编码向量；x_i为待编码向量；为用于训练的局部特征的数量；为基矩阵，基矩阵中每一列表示一个基向量，所有列构成的集合为字典，为字典中每个原子的维度，为字典中原子的数量；是B的凸可行集；是局部特征的重构系数；为对角正定权值矩阵，用来控制编码的局部性，根据每一个基向量与待编码向量x_i的相似度为其分配不同的权重，待编码向量离基向量越近，重构系数越大，权值d_ij为距离函数，例如：或σ用于调节距离增大时权重d_ij的衰减速度；λ为拉格朗日乘子，用于折中目标函数的重构误差项与局部性项的重要程度；

(5.2)特征提取：

首先对目标候选区域进行滑动网格稠密采样，得到均匀分布的图像块，设目标候选区域尺寸为图像块尺寸为patchsize×patchsize，步长为s，则得到的图像块的个数为：

其中，floor(·)为向下取整运算符；

其次，对每一个图像块提取SIFT特征描述子：先将图像块分为4×4的子区域，再对每一个区域进行8方向的梯度直方图统计，最后将每个子区域的归一化梯度直方图拼接起来得到该图像块的SIFT特征描述子，其维度为：

目标候选区域局部特征提取结果可以表示为式(15)：

采用k近邻搜索，选取所述字典中距离待编码描述子最近的k个基，形成子字典；然后，使用子字典重构特征描述子；最后，将得到的编码在非临近基的位置补0，得到编码向量

对已提取的目标候选区域局部特征进行局部约束线性编码的结果表示为式(16)：

采用空间金字塔池化操作来表示目标候选区域特征：

(6.1)对目标候选区域进行金字塔划分，将其分为L层，每层细分为不同分辨率的网格：第l层划分为2^l×2^l个网格，l＝0，1，...，L-1；

(6.2)进行池化操作，即对每一个网格中的局部约束线性编码特征进行最大池化，得到每个网格的池化特征；

(6.3)将不同分辨率下每个网格的池化特征级联起来得到输入图像最终的特征表达，其维度为：

其中是每个网格的池化特征的维度：