[发明专利]一种基于特征融合粒子滤波的仿生机器鱼视觉跟踪方法

权利要求书

1.一种基于特征融合粒子滤波的仿生机器鱼视觉跟踪方法，其特征在于，包括如下步骤：

步骤1，初始化：利用矩形框手动选定跟踪目标，提取目标的HSV颜色直方图、运动边缘直方图，组成目标模板；选取目标模板中的粒子数量为N，粒子权重为

其中，运动边缘直方图的提取方法如下：

设h_k,h_k+1为相邻的两帧图像，计算这两帧的绝对差分图像D_k，对差分图像D_k求梯度，得到运动边缘图像E_k，然后利用运动边缘的幅值或方向角获得运动边缘直方图；

所述目标的HSV颜色直方图的提取方法如下：

从摄像头采集的视频图像获取矩形框区域RGB颜色模型，接着完成RGB向HSV颜色空间的转换；将HSV颜色空间的每个子空间H、S、V分别分成Nh、Ns、Nv个等级，将整个HSV空间的颜色直方图划分为NN＝Nh×Ns+Nv个颜色子空间；根据目标区域图像像素的(h,s,v)分量所属的等级，得到该像素属于哪个颜色子空间中，从而得到目标的HSV颜色直方图；

步骤2，粒子状态估计：

k时刻第i个粒子状态估计为：

其中，分别表示第i个粒子在k-2、k-1、k时刻的状态，表示k-1时刻所有粒子的状态估计值的均值；w_k-1为k-1时刻系统随机过程噪声；A₁、A₂、B为常量；

步骤3，构建以k时刻第i个粒子状态估计的位置为中心的、与目标模板相同大小的候选区域；i＝1,2，…，N；提取当前帧所有N个候选区域的HSV颜色直方图和运动边缘直方图；

步骤4，修改粒子权重值：

步骤4.1，采用式(1)计算基于颜色直方图的第i个粒子的权值更新值：

其中，分别为第k-1、k帧基于颜色直方图的第i个粒子的权值更新值，ξ_clour₁ⁱ＝1；为颜色信息的观测似然函数；σ为高斯方差；为第i个粒子对应的候选区域与目标模板颜色特征的Bhattacharyya距离，为颜色特征的Bhattacharyya系数，q＝(q_(u))_{u＝1,2,...,NN}为目标模板的颜色直方图，为第k帧第i个粒子对应的候选区域的颜色直方图；NN为颜色直方图的颜色子空间数；

步骤4.2，采用式(2)计算基于运动边缘直方图的第i个粒子的权值更新值：

其中分别为第k-1、k帧基于运动边缘直方图的第i个粒子的权值更新值，为运动边缘信息的观测似然函数；σ₁为高斯方差，为第i个粒子对应的候选区域与目标模板运动边缘特征的Bhattacharyya距离，其中为运动边缘特征的Bhattacharyya系数，γ₁＝(γ_1(u))_{u＝1,2,...,MM}为目标模板的运动边缘直方图，为第k帧第i个粒子对应的候选区域的运动边缘直方图；MM为运动边缘直方图的幅度或方向子空间数；

步骤4.3，采用式(3)进行特征融合的粒子权重更新：

其中，分别为第k-1、k帧的特征融合的第i个粒子的权重更新值；为联合似然函数：

其中，η₁,η₂为归一化权值，其中d_{1_min},d_{2_min}分别为颜色特征和运动边缘特征中所有候选区域和目标模板之间的最小Bhattacharyya距离；

归一化粒子权重

步骤5，利用步骤4获得的粒子权重得到k时刻目标的状态估计为：

步骤6，目标模板更新：

设定一个Bhattacharyya系数的阈值，若当前帧所有候选区域的颜色特征和运动边缘特征的Bhattacharyya系数中的最大值小于设定的阈值，进行目标模板更新，按照步骤1的方法建立新的目标模板；然后按照下式进行目标模板更新：

M_opt＝β×M_old+(1-β)×M_now

其中，M_opt为更新后的目标模板，M_old为更新前的目标模板，M_now为新建的目标模板，β表示更新前的目标模板在更新后的目标模板中所占的比例；

步骤7，重采样粒子：

计算粒子退化程度N_eff，并与设定阈值N_th比较，如果N_effN_th，则重采样粒子，并令新的粒子具有相等的权重值否则，不重采样粒子；

步骤8：读取下一张图像，重复执行步骤2～8。