伙伴云客服论坛»论坛 S区 S软件开发 查看内容

0 评论

0 收藏

分享

OpenCV使用稀疏光流实现视频对象跟踪的方法详解

目录

    1、概述2、代码示例3、图像演示


1、概述

案例:使用稀疏光流实现对象跟踪
稀疏光流API介绍:
  1. calcOpticalFlowPyrLK( InputArray prevImg, InputArray nextImg,
  2.                                         InputArray prevPts, InputOutputArray nextPts,
  3.                                         OutputArray status, OutputArray err,
  4.                                         Size winSize = Size(21,21), int maxLevel = 3,
  5.                                         TermCriteria criteria = TermCriteria(TermCriteria::COUNT+TermCriteria::EPS, 30, 0.01),
  6.                                         int flags = 0, double minEigThreshold = 1e-4 );
复制代码
    prevImg:视频前一帧图像/金字塔,单通道CV_8UC1nextImg:视频后一帧图像/金字塔,单通道CV_8UC1preVPts:前一帧图像的特征向量(输入)需要找到流的2D点的矢量(vector of 2D points for which the flow needs to be found;);点坐标必需是单精度浮点数nextPts:后一帧图像的特征向量(输出),输出二维点的矢量(具有单精度浮点坐标),包含第二图像中输入特征的计算新位置;当传送OPTFLOW_USE_INITIAL_FLOW标志时,向量必需与输入中的大小相同。status:输出状态向量(无符号字符);假设找到相应特征的流,则向量的每个元素设置为1,否则设置为0err:输出错误的矢量; 向量的每个元素都设置为相应特征的错误,错误度量的类型可以在flags参数中设置; 假设未找到流,则未定义错误(使用status参数查找此类情况)winSize:每个金字塔等级的搜索窗口的winSize大小maxLevel:基于0的最大金字塔等级数;假设设置为0,则不使用金字塔(单级),假设设置为1,则使用两个级别,依此类推;假设将金字塔传送给输入,那么算法将使用与金字塔一样多的级别,但不超越maxLevelcriteria:停止条件,指定迭代搜索算法的终止条件(在指定的最大迭代次数criteria.maxCount之后或当搜索窗口挪动小于criteria.epsilon时)。flags:操作标志,OPTFLOW_USE_INITIAL_FLOW使用初始估计,存储在nextPts中;假设未设置标志,则将prevPts复制到nextPts并将其视为初始估计。OPTFLOW_LK_GET_MIN_EIGENVALS使用最小特征值作为误差丈量(参见minEigThreshold描绘);假设没有设置标志,则将原稿四周的色块和挪动点之间的L1间隔除以窗口中的像素数,用作误差丈量minEigThreshold:算法计算光流方程的2x2正常矩阵的最小特征值,除以窗口中的像素数;假设此值小于minEigThreshold,则过滤掉相应的功能并且不处置其流程,因而它允许删除坏点并获得性能提升
算法实现步骤:
1.实例化VideoCapture
2.循环读取视频数据
3.视频帧灰度转换
4.执行角点检测
5.保管角点检测的特征数据
6.初始化时假设检测到前一帧为空,把当前帧的灰度图像给前一帧
7.执行光流跟踪,并输出跟踪后的特征向量
8.遍历光流跟踪的输出特征向量,并得到间隔和状态都符合预期的特征向量。让后将其重新填充到fpts[1]中备用
9.重置集合大小
10.绘制光流线
11.交换特征向量的输入和输出
12.将用于跟踪的角点绘制出来
13.展示最终的跟踪效果
14.循环3~13步骤
15.完毕

2、代码示例
  1. KLT_Object_Tracking::KLT_Object_Tracking(QWidget *parent)
  2.     : MyGraphicsView{parent}
  3. {
  4.     isShowLine = false;
  5.     this->setWindowTitle("KLT稀疏光流实现对象跟踪");
  6.     QPushButton *btn = new QPushButton(this);
  7.     btn->setText("选择视频");
  8.     connect(btn,&QPushButton::clicked,[=](){
  9.         //选择视频
  10.         path = QFileDialog::getOpenFileName(this,"请选择视频","/Users/yangwei/Downloads/",tr("Image Files(*.mp4 *.avi)"));
  11.         qDebug()<<"视频途径:"<<path;
  12.         startKltTracking(path.toStdString().c_str());
  13.     });
  14.     //
  15.     QButtonGroup * group = new QButtonGroup(this);
  16.     QRadioButton * radioNo = new QRadioButton(this);
  17.     radioNo->setText("否");
  18.     radioNo->setChecked(true);
  19.     QRadioButton *radioYes = new QRadioButton(this);
  20.     radioYes->setText("是");
  21.     group->addButton(radioNo,0);
  22.     group->addButton(radioYes,1);
  23.     radioNo->move(0,btn->y()+btn->height()+20);
  24.     radioYes->move(radioNo->x()+radioNo->width()+20,btn->y()+btn->height()+20);
  25.     connect(radioNo,&QRadioButton::clicked,[=](){
  26.         isShowLine = false;//显示光流线
  27.     });
  28.     connect(radioYes,&QRadioButton::clicked,[=](){
  29.         isShowLine = true;//不显示光流线
  30.     });
  31. }
  32. void KLT_Object_Tracking::startKltTracking(const char* filePath){
  33.     //【1】实例化VideoCapture并翻开视频
  34.     VideoCapture capture;//实例化视频捕获器
  35.     capture.open(filePath);//翻开视频文件(或摄像头)
  36.     if(!capture.isOpened()){//检测文件是否翻开,假设没翻开直接退出
  37.         qDebug()<<"无法翻开视频";
  38.         return;
  39.     }
  40.     Mat frame,gray;
  41.     vector<Point2f> features;//检测出来的角点集合
  42.     vector<Point2f> inPoints;//这个主要是为了画线用的
  43.     vector<Point2f> fpts[2];//[0],存入的是是二维特征向量,[1]输出的二维特征向量
  44.     Mat pre_frame,pre_gray;
  45.     vector<uchar> status;//光流输出状态
  46.     vector<float> err;//光流输出错误
  47.     //【2】循环读取视频
  48.     while(capture.read(frame)){//循环读取视频中每一帧的图像
  49.         //【3】将视频帧图像转为灰度图
  50.         cvtColor(frame,gray,COLOR_BGR2GRAY);//ps:角点检测输入要求单通道
  51.         //【4】假设特征向量(角点)小于40个我们就重新执行角点检测
  52.         if(fpts[0].size()<40){//假设小于40个角点就重新开端执行角点检测
  53.             //执行角点检测
  54.             goodFeaturesToTrack(gray,features,5000,0.01,10,Mat(),3,false,0.04);
  55.             //【5】将检测到的角点放入fpts[0]中作为,光流跟踪的输入特征向量
  56.             //将检测到的角点插入vector
  57.             fpts[0].insert(fpts[0].begin(),features.begin(),features.end());
  58.             inPoints.insert(inPoints.end(),features.begin(),features.end());
  59.             qDebug()<<"角点检测执行完成,角点个数为:"<<features.size();
  60.         }else{
  61.             qDebug()<<"正在跟踪...";
  62.         }
  63.         //【6】初始化的时候假设检测到前一帧为空,这个把当前帧的灰度图像给前一帧
  64.         if(pre_gray.empty()){//假设前一帧为空就给前一帧赋值一次
  65.             gray.copyTo(pre_gray);
  66.         }
  67.         //执行光流跟踪
  68.         qDebug()<<"开端执行光流跟踪";
  69.         //【7】执行光流跟踪,并将输出的特征向量放入fpts[1]中
  70.         calcOpticalFlowPyrLK(pre_gray,gray,fpts[0],fpts[1],status,err);
  71.         qDebug()<<"光流跟踪执行完毕";
  72.         //【8】遍历光流跟踪的输出特征向量,并得到间隔和状态都符合预期的特征向量。让后将其重新填充到fpts[1]中备用
  73.         int k =0;
  74.         for(size_t i=0;i<fpts[1].size();i++){//循环遍历二维输出向量
  75.             double dist = abs(fpts[0][i].x - fpts[1][i].x) + abs(fpts[0][i].y - fpts[1][i].y);//特征向量挪动间隔
  76.             if(dist>2&&status[i]){//假设间隔大于2,status=true(正常)
  77.                 inPoints[k] = inPoints[i];
  78.                 fpts[1][k++] = fpts[1][i];
  79.             }
  80.         }
  81.         //【9】重置集合大小(由于有错误/不符合条件的输出特征向量),只拿状态正确的
  82.         //重新设置集合大小
  83.         inPoints.resize(k);
  84.         fpts[1].resize(k);
  85.         //【10】绘制光流线,这一步要不要都行
  86.         //绘制光流线
  87.         if(isShowLine){
  88.             for(size_t i = 0;i<fpts[1].size();i++){
  89.                 line(frame,inPoints[i],fpts[1][i],Scalar(0,255,0),1,8,0);
  90.                 circle(frame, fpts[1][i], 2, Scalar(0, 0, 255), 2, 8, 0);
  91.             }
  92.         }
  93.         qDebug()<<"特征向量的输入输出交换数据";
  94.         //【11】交换特征向量的输入和输出,(循环往复/进入下一个循环),此时特征向量的值会递减
  95.         std::swap(fpts[1],fpts[0]);//交换特征向量的输入和输出,此处焦点的总数量会递减
  96.         //【12】将用于跟踪的角点绘制出来
  97.         //将角点绘制出来
  98.         for(size_t i = 0;i<fpts[0].size();i++){
  99.             circle(frame,fpts[0][i],2,Scalar(0,0,255),2,8,0);
  100.         }
  101.         //【13】重置前一帧图像(每一个循环都要刷新)
  102.         gray.copyTo(pre_gray);
  103.         frame.copyTo(pre_frame);
  104.         //【14】展示最终的效果
  105.         imshow("frame",frame);
  106.         int keyValue = waitKey(100);
  107.         if(keyValue==27){//假设用户按ese键退出播放
  108.             break;
  109.         }
  110.     }
  111. }
复制代码
3、图像演示

OpenCV使用稀疏光流实现视频对象跟踪的方法详解-1.jpg

OpenCV使用稀疏光流实现视频对象跟踪的方法详解-2.jpg

到此这篇关于OpenCV使用稀疏光流实现视频对象跟踪的方法详解的文章就介绍到这了,更多相关OpenCV视频对象跟踪内容请搜索网站以前的文章或继续阅读下面的相关文章希望大家以后多多支持网站!

回复

举报 使用道具

全部回复
暂无回帖,快来参与回复吧
本版积分规则 高级模式
B Color Image Link Quote Code Smilies

萌晓
注册会员
主题 18
回复 22
粉丝 0
|网站地图
快速回复 返回顶部 返回列表