CN107341820A - 一种融合Cuckoo搜索和KCF的突变运动目标跟踪方法 - Google Patents

一种融合Cuckoo搜索和KCF的突变运动目标跟踪方法 Download PDF

Info

Publication number
CN107341820A
CN107341820A CN201710531498.8A CN201710531498A CN107341820A CN 107341820 A CN107341820 A CN 107341820A CN 201710531498 A CN201710531498 A CN 201710531498A CN 107341820 A CN107341820 A CN 107341820A
Authority
CN
China
Prior art keywords
mrow
target
tracking
kcf
msup
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710531498.8A
Other languages
English (en)
Other versions
CN107341820B (zh
Inventor
张焕龙
张建伟
陈宜滨
张秀娇
吴青娥
蒋斌
钱晓亮
贺振东
王延峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhengzhou University of Light Industry
Original Assignee
Zhengzhou University of Light Industry
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhengzhou University of Light Industry filed Critical Zhengzhou University of Light Industry
Priority to CN201710531498.8A priority Critical patent/CN107341820B/zh
Publication of CN107341820A publication Critical patent/CN107341820A/zh
Application granted granted Critical
Publication of CN107341820B publication Critical patent/CN107341820B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/20Analysis of motion
    • G06T7/277Analysis of motion involving stochastic approaches, e.g. using Kalman filters
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10016Video; Image sequence

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开了一种融合Cuckoo搜索和KCF的突变运动目标跟踪方法,步骤如下:初始化目标状态参数和优化方法初始参数;采用KCF跟踪方法获得目标前几帧内的最大响应值,计算置信度阈值初始值;依据当前帧最大响应值与置信度阈值的关系,确定不同的基样本图像产生方式:当大于置信度阈值,随机选取基图像样本,执行KCF方法跟踪目标;当小于置信度阈值,采用Cuckoo搜索机制获取全局最优的目标预测状态,产生新的基图像样本,执行KCF方法跟踪目标;动态更新置信度阈值,重复上述步骤实现目标跟踪。本发明有效地实现动态摄像机下目标存在帧间突变运动的持续性跟踪,实现帧间存在运动突变目标的精确跟踪,提高了在复杂场景下跟踪方法的适应能力。

Description

一种融合Cuckoo搜索和KCF的突变运动目标跟踪方法
技术领域
本发明属于目标跟踪的技术领域,具体涉及一种融合Cuckoo搜索和KCF的突变运动目标跟踪方法,实现目标的持续性跟踪,尤其涉及目标在视频相邻图像帧间产生运动突变情况下的状态空间搜索机制和目标持续跟踪能力。
背景技术
受到跟踪环境、目标自身和摄像机成像等因素的影响,兴趣目标在被跟踪过程中往往会出现运动突变的异常现象,这使传统基于运动平滑性假设条件提出的众多跟踪算法难以适应而导致失败。目标的突变运动已经成为一个算法能否实现持续性跟踪的重要因素之一,必须在传统目标跟踪框架下设计能够自适应运动突变问题的新方法,保证目标跟踪算法的鲁棒性。
视频目标跟踪算法研究主要分为两类:基于概率方法和基于决策方法。基于概率方法将跟踪过程视为贝叶斯框架下的状态评估问题,其中代表性算法是基于粒子滤波的视频目标跟踪方法,其能够处理非线性和多模态等问题,获得好的跟踪效果。然而,当面临目标运动突变时,基于概率方法往往通过扩大粒子数目以完全覆盖状态空间的方式提高跟踪的性能,这将导致较大的计算代价,难以满足目标跟踪的实时性要求。基于决策方法通过在每帧图像中采用迭代搜索方式寻找与真实目标最近似的区域实现目标定位,往往将跟踪过程转换为求解最优化问题来解决。为了使算法适应运动突变的问题,全局优化方法被引入视频目标跟踪框架,如退火粒子滤波方法、基于粒子群最优化的跟踪方法,基于马尔可夫蒙特卡洛采样的跟踪方法等。上述方法存在两个明显的问题:(1)在每帧内实现最优化过程,均需要全局搜索以获得目标定位,显然算法运行效率难以提高;(2)启发式优化策略包含较多参数调节,这使得跟踪算法针对较强,难以泛化。
因此,必须寻找一种既能利用全局优化方法适应目标运动突变问题,又能保证算法的运行效率,同时,优化策略应具有较少的模型调节参数以提高算法的通用性。
发明内容
针对现有运动目标跟踪方法运行效率不高、通用性差的技术问题,本发明提出一种融合Cuckoo搜索和KCF的突变运动目标跟踪方法,将布谷鸟(Cuckoo)搜索策略引入到核相关滤波跟踪器(KCF)设计中,依据置信度阈值确定产生基样本图像的方式,以增强评估突变运动状态的能力;将空间域中的回归问题等效为相关滤波操作,将耗时的卷积运算用点乘运算代替,提高了算法的运行效率,最终实现能够适应突变运动的视频目标跟踪方法。
为了解决上述的技术问题,本发明的技术方案是:一种融合Cuckoo搜索和KCF的突变运动目标跟踪方法,其步骤如下:
步骤一:初始化目标状态参数、Cuckoo搜索机制和KCF跟踪方法的初始参数;
步骤二:利用KCF跟踪方法获得前k帧的k个最大响应值,计算置信度阈值初始阈值;
步骤三:依据当前帧最大响应值与置信度阈值的关系,确定不同的基样本图像产生方式:当大于置信度阈值,随机选取基图像样本,执行KCF方法跟踪目标;当小于置信度阈值,采用Cuckoo搜索机制获取全局最优的目标预测状态,产生新的基图像样本,执行KCF方法跟踪目标;
步骤四:采用新的最大响应值动态更新置信度阈值,重复上述操作实现突变运动目标的有效跟踪。
所述目标状态参数为:读取第一帧图像的数据信息,确定目标在第一帧图像内的状态参数[px,py,width,high],其中,(px,py)为目标左上角像素点的坐标值,width为目标宽度,high为目标高度;Cuckoo搜索机制的参数为:初始图像块位置数num,位置发现概率pa,迭代优化次数iter;KCF跟踪方法的参数为:候选区域与目标区域的比例因子γ,设置计算置信度阈值的集合元素长度Len。
所述KCF跟踪方法实现的步骤为:
A、依据基样本图像构造循环矩阵X;
候选区域的宽和高分别为γ*width和γ*high,γ>1,将目标拉成列向量x=[x1,x2,…xn]T作为基样本向量,n=width*high,依据基样本向量的循环移位操作产生其他的候选样本,基样本向量和其产生的n-1个候选样本形成循环矩阵:
依据循环矩阵能够通过离散傅里叶变换实现对角化的性质,循环矩阵等价为:其中,F是离散傅里叶矩阵,是对基样本向量x进行离散傅里叶变换diag表示向量对角化操作;
B、跟踪问题的时域-频域转化,实现分类问题的求解;
xi为第i个训练样本,yi是其对应的回归值,训练的目标是求最小化平方误差下的权值w,从而获得决策函数f(z)=wTz,z是观测样本或称为测试样本;将跟踪的目标样本从候选样本中分离,权值w通过下式获得:
其中,λ控制过度拟合的正则化因子;在复数域下权值w解的向量描述形式为:
w=(XHX+λI)-1XHy,
其中,循环矩阵X的每行由训练样本xi形成,y是回归值yi的其中任意一个,I为识别矩阵,XH是共轭转换XH=(X*)T,X*是循环矩阵X的复共轭矩阵;
将循环矩阵的频域表示形式带入线性求解公式:
则上式改写成: 为逐元素的点乘操作,代入到权值w的向量求解表达式中获得频率域内的表达式如下:
其中,表示y的傅里叶变换;
利用核技巧将一个线性问题的输入映射到非线性的特征空间上,w在特征空间上被表示为问题转换成求解系数αi和非线性映射的关系,获得αi的系数组成的向量α的频率域表示:
其中,kxx是核循环矩阵k=C(kxx)的第一行元素,^代表向量的傅里叶变换;
C、获得最大响应值,确定跟踪目标;
通过上式获得向量α的值,得到线性回归问题中的权值w,然后通过分类器的分类方式实现目标检测,选取候选样本z,则:
其中,kxz是基样本向量x和观测样本z的核相关性,f(z)是一个向量、维度和回归值y一样,将所有候选样本和目标样本依据相关性大小排序,选择最大的值作为目标在该帧图像中获得的最大响应值response,其相对应的候选目标被确定为跟踪结果;
D、依据前k帧的最大响应值,计算置信度阈值初值:
如果k<Len,则返回步骤二重新执行;
如果k=Len,依据获得的k个最大响应值,计算初始置信度阈值如下:
其中,median(·)表示对集合的元素取中间值;
如果k>Len,依据第k,k-1,k-2,k-3,k-4帧获得的最大响应值动态更新集合Rk={responsek,responsek-1…responsek-4},则第k帧的置信度阈值为:
所述确定不同的基样本图像产生方式的方法是:当k+1帧图像到来时,如果采用布谷鸟全局搜索获得基样本图像,确定基样本图像块区域,构造循环矩阵X,采用KCF跟踪方法获得最大响应值,实现突变目标的跟踪;如果则依据第k帧目标在k+1帧内的映射位置附近随机获得基样本图像块,执行KCF跟踪方法对运动平滑的目标跟踪。
所述布谷鸟搜索的策略为:
A、产生初始位置,构造目标函数:随机选取num个初始化鸟窝,在各个位置上选取图像块,大小和跟踪目标相同,构造优化过程采用的目标函数;目标图像块和候选图像块分别用TX和CY表示,将它们作为随机变量,计算它们的相似性度如下:
其中,D(·)表示方差,Cov(·)表示协方差,M和N分别是目标图像块TX和候选图像块CY的HOG特征,目标函数定义为:
E=2+2*ρ(TX,CY);
B、预测位置的产生机制:布谷鸟搜索方法采用莱维飞行准则对每个初始化的位置产生一一对应的新位置,Px=[px,py]位置的更新公式为:
其中,分别表示第i个图像块在第t+1和t次迭代时的位置,β表示步长控制量,代表点乘操作,Levy(λ)表示莱维飞行随机搜索的跳跃式路径,β按下式计算:
β=u/|v|2/3
其中,u和v分别是服从正态分布的随机数;
依据莱维飞行规则对初始化的num个起始位置产生新的位置Ssum=(S1,S2…Snum),计算新位置上图像块与目标图像的相似度量函数值,保留函数值最大的新位置Smax,然后比较新解、旧解位置上图像块的能量函数值,将目标函数E值大的解作为新位置保留,确定新位置的集合元素;
C、依据发现概率pa,更新新的位置:依据布谷鸟搜索理论,如果新解被发现rand≥pa,则预测的位置被废弃,需要按照高斯噪声形式随机产生新位置更新原来的位置;如果rand<pa,则预测的位置继续保留,这里rand表示随机数产生函数;
D、依据所有新位置和目标计算获得的能量值,确定最优的位置作为本次迭代的最优输出值
E、一次迭代完成,判断是否达到最大的迭代次数iter,达到终止条件结束,输出最优的位置作为基样本图像块;否则,进入下次迭代转到B步。
本发明的有益效果:采用全局最优的运动状态预测机制,获得有效的基样本图像,保证循环矩阵能够覆盖目标状态空间,以弥补传统跟踪方法无法适应运动突变的问题;利用空间域和频率域处理数据的特征,在频率域实现候选目标样本响应值的确定,提高了运行效率;依据置信度阈值自适应选择不同跟踪模式,能够适应运动突变问题的同时也兼顾运行效率,避免传统方法因采用全局优化方法解决运动问题而带来计算代价过大的现象,能够很好适应目标的运动突变问题,后续对跟踪目标的识别、理解和分析具有重要意义。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图理解算法的本质。
图1为本发明的流程图。
图2为本发明与其他算法的距离精度比较示意图。
图3为本发明与其他算法的重叠率比较示意图。
图4为本发明与其他算法的跟踪效果示意图。
具体实施方式
下面将结合本发明实施中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
一种融合Cuckoo搜索和KCF的突变运动目标跟踪方法,将布谷鸟(Cuckoo)搜索机制引入KCF跟踪问题中,提出能够适应平滑和突变两种运动模式的新型跟踪框架;提出全局性的运动评估机制,能够对目标的不确定性运动状态进行预测,提高基样本图像的获取能力;利用置信度阈值切换不同的模型,通过信息的交互提高循环矩阵的可靠性,能够适应目标在图像帧间存在的运动突变问题。
本发明的思想在于:(1)将布谷鸟搜索技术引入到KCF跟踪问题中,给突变运动目标跟踪方法研究指出一个新的研究思路;(2)基于布谷鸟搜索机制能够利用多元启发提高搜索效率和利用莱维飞行避免局部陷入问题的优点,提出能够有效评估目标突变运动的运动模型;(3)建立一种新的跟踪算法框架,能够兼顾平滑和突变两种运动模式下的目标持续性跟踪。
本发明用于实施的硬件环境为:Intel(R)Core(TM)i3CPU 3.2G计算机、4GB内存、1G显卡,运行的软件环境是:Matlab R20012a和Windows 7,使用的数据为实验室拍摄的视频和吴毅老师公布的数据库。
如图1所示,本发明实例提供了一种融合Cuckoo搜索和KCF的突变运动目标跟踪方法,其步骤如下:
步骤一:初始化目标状态参数和Cuckoo搜索方法的参数。
读取第一帧图像的数据信息,确定目标在第一帧图像内的状态参数[px,py,width,high],其中,(px,py)为目标左上角像素点的坐标值,width为目标宽度,high为目标高度。设置Cuckoo搜索的初始图像块位置数num,位置发现概率pa,迭代优化次数iter;设置KCF跟踪方法中候选区域与目标区域的比例因子γ=1.3,设置计算置信度阈值的集合元素长度Len=5。Cuckoo搜索初始鸟窝数目即图像块位置数num=250,位置发现概率pa=0.5,迭代优化次数iter=100。
步骤二:利用KCF跟踪方法获得最大响应值,确定跟踪目标。
当新的图像到来时,采用KCF跟踪方法获得最大响应值,依据最大响应值定位目标。KCF方法实现主要分为三步:
A、依据基样本图像构造循环矩阵X。
基样本图像依据上一帧目标图像在当前帧映射的位置进行随机性采样。目标的宽和高分别为width和high,候选样本图像块和目标的尺寸一样。候选区域的宽和高分别为γ*width和γ*high,这里γ>1。将目标拉成列向量x=[x1,x2,…xn]T,n=width*high,将其作为基样本向量,依据基样本向量的循环移位操作产生其他的候选样本,候选样本集合包含了基样本向量和其产生的n-1个候选样本,形成的循环矩阵如下:
循环矩阵X的第一行是基样本向量的转置xT。依据循环矩阵能够通过离散傅里叶变换实现对角化的性质,循环矩阵可以等价为:
其中,F是离散傅里叶矩阵,是个常量矩阵,是对基样本向量x进行离散傅里叶变换diag表示向量对角化操作。对于任意一个向量有可以获得离散傅里叶矩阵F。
B、跟踪问题的时域-频域转化,实现分类问题的求解。
目标跟踪过程可以视为目标样本和候选样本持续分类的操作。已知xi为第i个训练样本,yi是其对应的回归值,训练的目标就是求最小化平方误差下的权值w,从而获得决策函数f(z)=wTz,z是观测样本或称为测试样本。将跟踪的目标样本从候选样本中分离,权值w通过下式获得:
其中,λ控制过度拟合的正则化因子。这个线性回归问题具有封闭解,在复数域下解的向量描述形式为:
w=(XHX+λI)-1XHy
其中,循环矩阵X的每行由基向量和其经过移位操作形成的所有候选样本的向量组成,xi表示任意一个候选样本,i表示候选样本的索引。y是回归值yi的其中任意一个值。I为识别矩阵,XH是共轭转换XH=(X*)T,X*是循环矩阵X的复共轭矩阵。对于实数的情况XH=(X)T
为了获得权值w,上式需要矩阵求逆,在大规模视频运算中,复杂度很大,运算耗时,将时域复杂运算转换到频域处理以提高运行效率。
将循环矩阵的频域表示形式带入线性求解公式:
定义逐元素的点乘操作符号为⊙,则上式改写成:
代入到权值w的向量求解表达式中,获得其频率域内的表达式如下:
其中,表示y的傅里叶变换。
时域中的求逆转换成频率域中逐元素除法效率大大提高,通过反傅里叶变换很容易获得权值w。实际跟踪中很多问题属于非线性求解问题,利用核技巧用x的非线性组合来表示权值问题转换成求解系数αi和非线性映射的关系。类似线性情况下的数学推导,获得α的频率域表示:
其中,α是αi的系数组成的向量。kxx是核循环矩阵k=C(kxx)的第一行元素,^代表向量的傅里叶变换。于是目标跟踪的分类问题转换到频率域中对向量α的求解,然后经反傅里叶变换获得时域中决策函数形式,实现高效的目标跟踪过程。
C、获得最大响应值,确定跟踪目标。
通过上式获得向量α的值,相当于在线性回归问题中获得了权值w,然后通过分类器的分类方式实现目标检测。选取候选样本z,其和目标样本x的维度相同,依据公式:
其中,kxz是基样本向量x和观测样本z的核相关性,f(z)是一个向量、维度和回归值y一样,获得最大响应值response的候选样本,确定为跟踪结果,即目标被再次定位。
步骤三:依据目标的最大响应值,确定置信度阈值。
如果k<5,则返回步骤二重新执行;
如果k=5,则按如下规则确定初始置信度阈值
依据获得的k个最大响应值,构造集合Rk={response1,response2…responsek},计算初始置信度阈值如下:
其中,median(·)表示对集合的元素取中间值;
如果k>5,依据第k,k-1,k-2,k-3,k-4帧获得的最大响应值动态更新集合Rk={responsek,responsek-1…responsek-4},则第k帧的置信度阈值为:
步骤四:依据当前帧的最大响应值与置信度阈值的关系,确定不同的基样本图像产生方式。
当k+1帧图像到来时,依据最大响应值responsek+1与置信度阈值比较,以确定基样本图像的产生形式。
如果采用布谷鸟全局搜索获得基样本图像,确定基样本图像块区域,结合KCF方法实现目标跟踪。布谷鸟搜索策略如下:
A、产生初始位置,构造目标函数。随机选取250个初始化鸟窝,即并行迭代优化方法的初始位置。在各个位置上选取图像块,大小和跟踪目标相同,构造优化过程采用的目标函数。在本发明中,目标图像块和候选图像块分别用TX和CY表示,将它们作为随机变量,计算它们的相似性度如下:
其中,D(·)表示方差,Cov(·)表示协方差,M和N分别是目标图像块TX和候选图像块CY的HOG特征。目标函数被定义如下:
E=2+2*ρ(TX,CY)。
目标函数的函数值即能量值表征两个图像块的相似性,其值越大说明越相似。
B、预测位置的产生机制。布谷鸟搜索方法采用莱维(Levy)飞行准则对每个初始化的位置产生一一对应的新位置,也就是沿着能量值降低的方向进行迭代,Px=[px,py]位置的更新公式为:
其中,分别表示第i个图像块在第t+1和t次迭代时的位置,β表示步长控制量,代表点乘操作,Levy(λ)表示Levy飞行随机搜索的跳跃式路径,β按下式计算:
β=u/|v|2/3
其中,u和v分别是服从正态分布的随机数。依据莱维飞行规则对初始化的num个起始位置产生新的位置Ssum=(S1,S2…Snum),计算新位置上图像块与目标图像的相似度量函数值,保留函数值最大的新位置Smax,然后比较新解、旧解位置上图像块的能量函数值,将目标函数值E大的解作为新位置保留,确定新解的集合元素。依据莱维飞行规则获得250个新解的位置,分别与迭代前位置上的图像块进行能量函数值E比较,保留能量值的位置作为最终的新位置,即重新确定250个较优的新位置。
C、依据位置发现概率,更新新的位置。依据布谷鸟搜索理论,如果新解被发现rand≥pa,则新解被废弃,需要按照高斯噪声形式随机产生新位置更新原来的位置;如果rand<pa,则预测的位置继续保留,这里rand表示随机数产生函数。
D、依据所有新位置和目标计算获得的能量值,确定一个最优的位置作为本次迭代的最优输出值
E、一次迭代完成,判断是否达到最大的迭代次数,达到终止条件结束,输出最优的位置作为基样本图像块。否则,进入下次迭代转到B步。
如果则依据第k帧目标在k+1帧内的映射位置附近随机获得基样本图像块,执行KCF方法对运动平滑的目标进行有效跟踪。
使用更新的基样本图像块构造循环矩阵X,采用KCF跟踪方法获得最大响应值,实现对突变目标的有效跟踪。
步骤五:依据最大响应值,确定新的基图像样本,返回步骤二,确定跟踪目标。
重复二-五步实现突变运动目标的跟踪。
本发明的实施步骤如下:初始化目标状态参数和优化模型参数;采用KCF跟踪方法获得目标前几帧内的最大响应值,计算置信度阈值初始值;依据当前帧最大响应值与置信度阈值的关系,确定不同的基样本图像产生方式:当大于置信度阈值,随机选取基图像样本,执行KCF方法跟踪目标,当小于置信度阈值,采用Cuckoo搜索机制获取全局最优的目标预测状态,产生新的基图像样本,执行KCF方法跟踪目标;采用新的最大响应值动态更新置信度阈值,重复上述操作实现突变运动目标的有效跟踪。本发明对动态摄像机下目标存在帧间突变运动时能够有效地实现持续性跟踪,实现了帧间存在运动突变的目标的精确跟踪,提高了在复杂场景下跟踪方法的适应能力。
采用定性和定量两种评价方式来评价本发明的有效性。定性评价方式采用本发明方法和各种代表性方法在部分帧内的跟踪效果图如图4所示。图4分别是Deer、Face2、FHC和ZT四个视频在部分帧内的跟踪效果图,其中,第一个和第二个视频中,目标在水平方向的帧间运动位移分别达到38和88个像素的大位移运动,第三个和第四个视频中目标在垂直方向上的帧间运动位移达到188和256个像素,这些突变运动而产生的帧间大位移使得其它跟踪算法表现出不适应,本发明提出的方法获得了好的跟踪效果。定量评价方式采用中心位置误差率和目标重叠率来评价,距离精度DP是指依据设定的阈值图像序列中能够成功跟踪目标的帧数与总帧数的比值。其中,阈值是通过跟踪结果与真实结果的中心误差值确定的,本发明设置阈值为0.5。中心误差是通过计算跟踪目标的中心位置和真实位置之间的欧氏距离,一般其值越小说明跟踪结果比较优。目标重叠率OP是指跟踪结果与真实目标区域面积与操作和并操作的比值,其值越大,说明跟踪结果越好。图2展示了本发明与代表性跟踪算法CSK、DSST、FCT、KCF和STC的距离精度DP值比较结果,图3展示相应跟踪算法的目标重叠率OP值的比较结果。综合图4结果表示,本发明提供的目标跟踪方法能够很好地解决目标突变的运动问题,获得了较好地跟踪表现。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (5)

1.一种融合Cuckoo搜索和KCF的突变运动目标跟踪方法,其特征在于,其步骤如下:
步骤一:初始化目标状态参数、Cuckoo搜索机制和KCF跟踪方法的初始参数;
步骤二:利用KCF跟踪方法获得前k帧的k个最大响应值,计算置信度阈值初始阈值;
步骤三:依据当前帧最大响应值与置信度阈值的关系,确定不同的基样本图像产生方式:当大于置信度阈值,随机选取基图像样本,执行KCF方法跟踪目标;当小于置信度阈值,采用Cuckoo搜索机制获取全局最优的目标预测状态,产生新的基图像样本,执行KCF方法跟踪目标;
步骤四:采用新的最大响应值动态更新置信度阈值,重复上述操作实现突变运动目标的有效跟踪。
2.根据权利要求1所述的融合Cuckoo搜索和KCF的突变运动目标跟踪方法,其特征在于,所述目标状态参数为:读取第一帧图像的数据信息,确定目标在第一帧图像内的状态参数[px,py,width,high],其中,(px,py)为目标左上角像素点的坐标值,width为目标宽度,high为目标高度;Cuckoo搜索机制的参数为:初始图像块位置数num,位置发现概率pa,迭代优化次数iter;KCF跟踪方法的参数为:候选区域与目标区域的比例因子γ,设置计算置信度阈值的集合元素长度Len。
3.根据权利要求1所述的融合Cuckoo搜索和KCF的突变运动目标跟踪方法,其特征在于,所述KCF跟踪方法实现的步骤为:
A、依据基样本图像构造循环矩阵X;
候选区域的宽和高分别为γ*width和γ*high,γ>1,将目标拉成列向量x=[x1,x2,…xn]T作为基样本向量,n=width*high,依据基样本向量的循环移位操作产生其他的候选样本,基样本向量和其产生的n-1个候选样本形成循环矩阵:
依据循环矩阵能够通过离散傅里叶变换实现对角化的性质,循环矩阵等价为:其中,F是离散傅里叶矩阵,是对基样本向量x进行离散傅里叶变换diag表示向量对角化操作;
B、跟踪问题的时域-频域转化,实现分类问题的求解;
xi为第i个训练样本,yi是其对应的回归值,训练的目标是求最小化平方误差下的权值w,从而获得决策函数f(z)=wTz,z是观测样本或称为测试样本;将跟踪的目标样本从候选样本中分离,权值w通过下式获得:
<mrow> <munder> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> <mi>w</mi> </munder> <munder> <mo>&amp;Sigma;</mo> <mi>i</mi> </munder> <msup> <mrow> <mo>(</mo> <mi>f</mi> <mo>(</mo> <msup> <mi>x</mi> <mi>i</mi> </msup> <mo>)</mo> <mo>-</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mi>&amp;lambda;</mi> <mo>|</mo> <mo>|</mo> <mi>w</mi> <mo>|</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> <mo>,</mo> </mrow>
其中,λ控制过度拟合的正则化因子;在复数域下权值w解的向量描述形式为:
w=(XHX+λI)-1XHy,
其中,循环矩阵X的每行由训练样本xi形成,y是回归值yi的其中任意一个,I为识别矩阵,XH是共轭转换XH=(X*)T,X*是循环矩阵X的复共轭矩阵;
将循环矩阵的频域表示形式带入线性求解公式:
<mrow> <msup> <mi>X</mi> <mi>H</mi> </msup> <mi>X</mi> <mo>=</mo> <mi>F</mi> <mo>&amp;CenterDot;</mo> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mrow> <mo>(</mo> <msup> <mover> <mi>x</mi> <mo>^</mo> </mover> <mo>*</mo> </msup> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msup> <mi>F</mi> <mi>H</mi> </msup> <mi>F</mi> <mo>&amp;CenterDot;</mo> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>g</mi> <mrow> <mo>(</mo> <mover> <mi>x</mi> <mo>^</mo> </mover> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msup> <mi>F</mi> <mi>H</mi> </msup> <mo>;</mo> </mrow>
则上式改写成:⊙为逐元素的点乘操作,代入到权值w的向量求解表达式中获得频率域内的表达式如下:
其中,表示y的傅里叶变换;
利用核技巧将一个线性问题的输入映射到非线性的特征空间上,w在特征空间上被表示为问题转换成求解系数αi和非线性映射的关系,获得αi的系数组成的向量α的频率域表示:
<mrow> <mover> <mi>&amp;alpha;</mi> <mo>^</mo> </mover> <mo>=</mo> <mfrac> <mover> <mi>y</mi> <mo>^</mo> </mover> <mrow> <msup> <mover> <mi>k</mi> <mo>^</mo> </mover> <mrow> <mi>x</mi> <mi>x</mi> </mrow> </msup> <mo>+</mo> <mi>&amp;lambda;</mi> </mrow> </mfrac> <mo>,</mo> </mrow>
其中,kxx是核循环矩阵k=C(kxx)的第一行元素,^代表向量的傅里叶变换;
C、获得最大响应值,确定跟踪目标;
通过上式获得向量α的值,得到线性回归问题中的权值w,然后通过分类器的分类方式实现目标检测,选取候选样本z,则:
其中,kxz是基样本向量x和观测样本z的核相关性,f(z)是一个向量、维度和回归值y一样,将所有候选样本和目标样本依据相关性大小排序,选择最大的值作为目标在该帧图像中获得的最大响应值response,其相对应的候选目标被确定为跟踪结果;
D、依据前k帧的最大响应值,计算置信度阈值初值:
如果k<Len,则返回步骤二重新执行;
如果k=Len,依据获得的k个最大响应值,计算初始置信度阈值如下:
<mrow> <msubsup> <mi>C</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>p</mi> </mrow> <mn>0</mn> </msubsup> <mo>=</mo> <mi>m</mi> <mi>e</mi> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mi>response</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>response</mi> <mn>2</mn> </msub> <mo>...</mo> <msub> <mi>response</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>,</mo> </mrow>
其中,median(·)表示对集合的元素取中间值;
如果k>Len,依据第k,k-1,k-2,k-3,k-4帧获得的最大响应值动态更新集合Rk={responsek,responsek-1…responsek-4},则第k帧的置信度阈值为:
<mrow> <msubsup> <mi>C</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>p</mi> </mrow> <mi>k</mi> </msubsup> <mo>=</mo> <mi>m</mi> <mi>e</mi> <mi>d</mi> <mi>i</mi> <mi>a</mi> <mi>n</mi> <mo>{</mo> <msub> <mi>response</mi> <mi>k</mi> </msub> <mo>,</mo> <msub> <mi>response</mi> <mrow> <mi>k</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>...</mo> <msub> <mi>response</mi> <mrow> <mi>k</mi> <mo>-</mo> <mn>4</mn> </mrow> </msub> <mo>}</mo> <mo>.</mo> </mrow>
4.根据权利要求1所述的融合Cuckoo搜索和KCF的突变运动目标跟踪方法,其特征在于,所述确定不同的基样本图像产生方式的方法是:当k+1帧图像到来时,如果采用布谷鸟全局搜索获得基样本图像,确定基样本图像块区域,构造循环矩阵X,采用KCF跟踪方法获得最大响应值,实现突变目标的跟踪;如果则依据第k帧目标在k+1帧内的映射位置附近随机获得基样本图像块,执行KCF跟踪方法对运动平滑的目标跟踪。
5.根据权利要求4所述的融合Cuckoo搜索和KCF的突变运动目标跟踪方法,其特征在于,所述布谷鸟搜索的策略为:
A、产生初始位置,构造目标函数:随机选取num个初始化鸟窝,在各个位置上选取图像块,大小和跟踪目标相同,构造优化过程采用的目标函数;目标图像块和候选图像块分别用TX和CY表示,将它们作为随机变量,计算它们的相似性度如下:
<mrow> <mi>&amp;rho;</mi> <mrow> <mo>(</mo> <mi>T</mi> <mi>X</mi> <mo>,</mo> <mi>C</mi> <mi>Y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>C</mi> <mi>o</mi> <mi>v</mi> <mrow> <mo>(</mo> <mi>M</mi> <mo>,</mo> <mi>N</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msqrt> <mrow> <mi>D</mi> <mrow> <mo>(</mo> <mi>M</mi> <mo>)</mo> </mrow> </mrow> </msqrt> <msqrt> <mrow> <mi>D</mi> <mrow> <mo>(</mo> <mi>N</mi> <mo>)</mo> </mrow> </mrow> </msqrt> </mrow> </mfrac> <mo>,</mo> </mrow>
其中,D(·)表示方差,Cov(·)表示协方差,M和N分别是目标图像块TX和候选图像块CY的HOG特征,目标函数定义为:
E=2+2*ρ(TX,CY);
B、预测位置的产生机制:布谷鸟搜索方法采用莱维飞行准则对每个初始化的位置产生一一对应的新位置,Px=[px,py]位置的更新公式为:
<mrow> <msubsup> <mi>Px</mi> <mi>i</mi> <mrow> <mi>t</mi> <mo>+</mo> <mn>1</mn> </mrow> </msubsup> <mo>=</mo> <msubsup> <mi>Px</mi> <mi>i</mi> <mi>t</mi> </msubsup> <mo>+</mo> <mi>&amp;beta;</mi> <mo>&amp;CirclePlus;</mo> <mi>L</mi> <mi>e</mi> <mi>v</mi> <mi>y</mi> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> <mo>,</mo> </mrow>
其中,分别表示第i个图像块在第t+1和t次迭代时的位置,β表示步长控制量,代表点乘操作,Levy(λ)表示莱维飞行随机搜索的跳跃式路径,β按下式计算:
β=u/|v|2/3
其中,u和v分别是服从正态分布的随机数;
依据莱维飞行规则对初始化的num个起始位置产生新的位置Ssum=(S1,S2…Snum),计算新位置上图像块与目标图像的相似度量函数值,保留函数值最大的新位置Smax,然后比较新解、旧解位置上图像块的能量函数值,将目标函数E值大的解作为新位置保留,确定新位置的集合元素;
C、依据发现概率pa,更新新的位置:依据布谷鸟搜索理论,如果新解被发现rand≥pa,则预测的位置被废弃,需要按照高斯噪声形式随机产生新位置更新原来的位置;如果rand<pa,则预测的位置继续保留,这里rand表示随机数产生函数;
D、依据所有新位置和目标计算获得的能量值,确定最优的位置作为本次迭代的最优输出值
E、一次迭代完成,判断是否达到最大的迭代次数iter,达到终止条件结束,输出最优的位置作为基样本图像块;否则,进入下次迭代转到B步。
CN201710531498.8A 2017-07-03 2017-07-03 一种融合Cuckoo搜索和KCF的突变运动目标跟踪方法 Active CN107341820B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710531498.8A CN107341820B (zh) 2017-07-03 2017-07-03 一种融合Cuckoo搜索和KCF的突变运动目标跟踪方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710531498.8A CN107341820B (zh) 2017-07-03 2017-07-03 一种融合Cuckoo搜索和KCF的突变运动目标跟踪方法

Publications (2)

Publication Number Publication Date
CN107341820A true CN107341820A (zh) 2017-11-10
CN107341820B CN107341820B (zh) 2018-04-20

Family

ID=60218394

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710531498.8A Active CN107341820B (zh) 2017-07-03 2017-07-03 一种融合Cuckoo搜索和KCF的突变运动目标跟踪方法

Country Status (1)

Country Link
CN (1) CN107341820B (zh)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108245194A (zh) * 2017-12-21 2018-07-06 四川省人民医院 基于布谷鸟优化策略的超声图像心脏流场运动估计方法
CN108446383A (zh) * 2018-03-21 2018-08-24 吉林大学 一种基于地理分布式数据查询的数据任务重分配方法
CN108734177A (zh) * 2018-05-17 2018-11-02 中国人民解放军陆军工程大学 双步相关滤波目标跟踪方法
CN108805909A (zh) * 2018-04-28 2018-11-13 哈尔滨工业大学深圳研究生院 在相关滤波框架下基于粒子滤波重检测的目标跟踪方法
CN109409235A (zh) * 2018-09-27 2019-03-01 Oppo广东移动通信有限公司 图像识别方法和装置、电子设备、计算机可读存储介质
CN110197472A (zh) * 2018-02-26 2019-09-03 四川省人民医院 一种用于超声造影图像稳定定量分析的方法和系统
CN110378260A (zh) * 2019-07-05 2019-10-25 国网安徽省电力有限公司铜陵供电公司 基于kcf的实时刀闸状态跟踪方法及系统
CN110490899A (zh) * 2019-07-11 2019-11-22 东南大学 一种结合目标跟踪的可变形施工机械的实时检测方法
CN110717930A (zh) * 2019-08-28 2020-01-21 河南中烟工业有限责任公司 一种基于扩展Wang-landau蒙特卡罗和KCF的突变运动目标跟踪方法
CN110751096A (zh) * 2019-10-21 2020-02-04 陕西师范大学 一种基于kcf轨迹置信度的多目标跟踪方法
CN110751670A (zh) * 2018-07-23 2020-02-04 中国科学院长春光学精密机械与物理研究所 一种基于融合的目标跟踪方法
CN111739055A (zh) * 2020-06-10 2020-10-02 新疆大学 一种红外点状目标跟踪方法
CN111784737A (zh) * 2020-06-10 2020-10-16 中国人民解放军军事科学院国防科技创新研究院 一种基于无人机平台的目标自动跟踪方法及系统
CN113239860A (zh) * 2021-05-28 2021-08-10 中建材信息技术股份有限公司 一种基于视频的烟火检测方法
CN113658216A (zh) * 2021-06-24 2021-11-16 北京理工大学 基于多级自适应kcf的遥感目标跟踪方法及电子设备
CN113239860B (zh) * 2021-05-28 2024-05-24 中建材信息技术股份有限公司 一种基于视频的烟火检测方法

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109461166A (zh) * 2018-10-26 2019-03-12 郑州轻工业学院 一种基于kcf混合mfo的快速运动目标跟踪方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106557774A (zh) * 2015-09-29 2017-04-05 南京信息工程大学 多通道核相关滤波的实时跟踪方法
CN106570486A (zh) * 2016-11-09 2017-04-19 华南理工大学 基于特征融合和贝叶斯分类的核相关滤波目标跟踪方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106557774A (zh) * 2015-09-29 2017-04-05 南京信息工程大学 多通道核相关滤波的实时跟踪方法
CN106570486A (zh) * 2016-11-09 2017-04-19 华南理工大学 基于特征融合和贝叶斯分类的核相关滤波目标跟踪方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
GURJIT SINGH WALIA 等: "Intelligent video target tracking using an evolutionary particle filter based upon improved cuckoo search", 《EXPERTS SYSTEMS WITH APPLICATIONS》 *
MING-LIANG GAO 等: "Visual tracking method based on cuckoo search algorithm", 《OPTICAL ENGINEERING》 *
TARIK LJOUAD 等: "A hybrid mobile object tracker based on the modified Cuckoo Search algorithm and the Kalman filter", 《PATTERN RECOGNITION》 *
刘剑 等: "一种基于多样性优化的视频目标跟踪方法", 《计算应用与软件》 *

Cited By (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108245194A (zh) * 2017-12-21 2018-07-06 四川省人民医院 基于布谷鸟优化策略的超声图像心脏流场运动估计方法
CN108245194B (zh) * 2017-12-21 2020-06-19 四川省人民医院 基于布谷鸟优化策略的超声图像心脏流场运动估计方法
CN110197472B (zh) * 2018-02-26 2023-08-08 四川省人民医院 一种用于超声造影图像稳定定量分析的方法和系统
CN110197472A (zh) * 2018-02-26 2019-09-03 四川省人民医院 一种用于超声造影图像稳定定量分析的方法和系统
CN108446383A (zh) * 2018-03-21 2018-08-24 吉林大学 一种基于地理分布式数据查询的数据任务重分配方法
CN108446383B (zh) * 2018-03-21 2021-12-10 吉林大学 一种基于地理分布式数据查询的数据任务重分配方法
CN108805909A (zh) * 2018-04-28 2018-11-13 哈尔滨工业大学深圳研究生院 在相关滤波框架下基于粒子滤波重检测的目标跟踪方法
CN108805909B (zh) * 2018-04-28 2022-02-11 哈尔滨工业大学深圳研究生院 在相关滤波框架下基于粒子滤波重检测的目标跟踪方法
CN108734177B (zh) * 2018-05-17 2021-06-29 中国人民解放军陆军工程大学 双步相关滤波目标跟踪方法
CN108734177A (zh) * 2018-05-17 2018-11-02 中国人民解放军陆军工程大学 双步相关滤波目标跟踪方法
CN110751670B (zh) * 2018-07-23 2022-10-25 中国科学院长春光学精密机械与物理研究所 一种基于融合的目标跟踪方法
CN110751670A (zh) * 2018-07-23 2020-02-04 中国科学院长春光学精密机械与物理研究所 一种基于融合的目标跟踪方法
CN109409235A (zh) * 2018-09-27 2019-03-01 Oppo广东移动通信有限公司 图像识别方法和装置、电子设备、计算机可读存储介质
CN110378260A (zh) * 2019-07-05 2019-10-25 国网安徽省电力有限公司铜陵供电公司 基于kcf的实时刀闸状态跟踪方法及系统
CN110378260B (zh) * 2019-07-05 2024-02-06 国网安徽省电力有限公司铜陵供电公司 基于kcf的实时刀闸状态跟踪方法及系统
CN110490899A (zh) * 2019-07-11 2019-11-22 东南大学 一种结合目标跟踪的可变形施工机械的实时检测方法
CN110717930A (zh) * 2019-08-28 2020-01-21 河南中烟工业有限责任公司 一种基于扩展Wang-landau蒙特卡罗和KCF的突变运动目标跟踪方法
CN110717930B (zh) * 2019-08-28 2024-02-20 河南中烟工业有限责任公司 一种基于扩展Wang-landau蒙特卡罗和KCF的突变运动目标跟踪方法
CN110751096A (zh) * 2019-10-21 2020-02-04 陕西师范大学 一种基于kcf轨迹置信度的多目标跟踪方法
CN110751096B (zh) * 2019-10-21 2022-02-22 陕西师范大学 一种基于kcf轨迹置信度的多目标跟踪方法
CN111784737A (zh) * 2020-06-10 2020-10-16 中国人民解放军军事科学院国防科技创新研究院 一种基于无人机平台的目标自动跟踪方法及系统
CN111739055A (zh) * 2020-06-10 2020-10-02 新疆大学 一种红外点状目标跟踪方法
CN111739055B (zh) * 2020-06-10 2022-07-05 新疆大学 一种红外点状目标跟踪方法
CN113239860A (zh) * 2021-05-28 2021-08-10 中建材信息技术股份有限公司 一种基于视频的烟火检测方法
CN113239860B (zh) * 2021-05-28 2024-05-24 中建材信息技术股份有限公司 一种基于视频的烟火检测方法
CN113658216A (zh) * 2021-06-24 2021-11-16 北京理工大学 基于多级自适应kcf的遥感目标跟踪方法及电子设备

Also Published As

Publication number Publication date
CN107341820B (zh) 2018-04-20

Similar Documents

Publication Publication Date Title
CN107341820B (zh) 一种融合Cuckoo搜索和KCF的突变运动目标跟踪方法
Ebert et al. Robustness via retrying: Closed-loop robotic manipulation with self-supervised learning
CN106169188B (zh) 一种基于蒙特卡洛树搜索的对象跟踪方法
CN113110509B (zh) 一种基于深度强化学习的仓储系统多机器人路径规划方法
CN105405151B (zh) 基于粒子滤波和加权Surf的抗遮挡目标跟踪方法
CN112348849A (zh) 一种孪生网络视频目标跟踪方法及装置
CN107424177A (zh) 基于连续相关滤波器的定位修正长程跟踪算法
Zhu et al. Off-road autonomous vehicles traversability analysis and trajectory planning based on deep inverse reinforcement learning
CN103259962A (zh) 一种目标追踪方法和相关装置
CN102708382B (zh) 一种基于变处理窗口和变坐标系的多目标跟踪方法
CN109940614B (zh) 一种融合记忆机制的机械臂多场景快速运动规划方法
CN103985120A (zh) 一种遥感图像多目标关联的方法
CN106447696A (zh) 一种基于双向sift流运动评估的大位移目标稀疏跟踪方法
CN109544600A (zh) 一种基于上下文相关和判别相关滤波器的目标跟踪方法
CN110189362B (zh) 基于多分支自编码对抗网络的高效目标跟踪方法
CN112651998A (zh) 基于注意力机制和双流多域卷积神经网络的人体跟踪算法
CN104299248A (zh) 利用前视声呐图像对水下多个动态目标运动预测的方法
CN110070565A (zh) 一种基于图像叠加的船舶轨迹预测方法
CN102592135B (zh) 融合目标空间分布和时序分布特征子空间的视觉跟踪方法
CN113344973A (zh) 一种基于时空正则化和特征可靠性评估的目标跟踪方法
CN111462180B (zh) 基于与或图aog的目标跟踪方法
CN109461166A (zh) 一种基于kcf混合mfo的快速运动目标跟踪方法
CN113920170B (zh) 结合场景上下文和行人社会关系的行人轨迹预测方法、系统及存储介质
Zhao et al. Gaussian processes for flow modeling and prediction of positioned trajectories evaluated with sports data
CN106127766A (zh) 基于空间耦合关系和历史模型的目标跟踪方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Zhang Huanlong

Inventor after: Wang Yanfeng

Inventor after: Zheng Yubin

Inventor after: Zhang Jianwei

Inventor after: Li Zuhe

Inventor after: Chen Yibin

Inventor after: Zhang Xiujiao

Inventor after: Wu Qinge

Inventor after: Jiang Bin

Inventor after: Qian Xiaoliang

Inventor after: He Zhendong

Inventor before: Zhang Huanlong

Inventor before: Zhang Jianwei

Inventor before: Chen Yibin

Inventor before: Zhang Xiujiao

Inventor before: Wu Qinge

Inventor before: Jiang Bin

Inventor before: Qian Xiaoliang

Inventor before: He Zhendong

Inventor before: Wang Yanfeng

GR01 Patent grant
GR01 Patent grant