【opencv】动态背景下运动目标检测 SURF配准差分
主要思路是,读入视频,隔帧采用SURF计算匹配的特征点,进而计算两图的投影映射矩阵,做差分二值化,连通域检测,绘制目标。
如果背景是静态的采用camshift即可。
本文方法速度debug下大概2-3帧,release下8-9帧(SURF部分,不包含连通域以及绘制),后续可增加选定目标,动态模版小邻域中跟踪目标。实现对动态背景下的运动目标检测,模版跟踪速度可达150帧。
环境:opencv2.4.9 + vs2012
#include <iostream>
#include <opencv2/opencv.hpp>
#include <opencv2/nonfree/nonfree.hpp> using namespace cv;
using namespace std;void main()
{//VideoCapture capture(0);VideoCapture capture("3.mov");Mat image01,image02,imgdiff;while (true){//隔两帧配准capture >> image01;if (image01.empty()){break;}capture >> image02;capture >> image02;if (image02.empty()){break;}//GaussianBlur(image02, image02, Size(3,3), 0);double time0 = static_cast<double>(getTickCount());//开始计时//灰度图转换 Mat image1,image2; cvtColor(image01,image1,CV_RGB2GRAY); cvtColor(image02,image2,CV_RGB2GRAY); //提取特征点 SurfFeatureDetector surfDetector(2500); // 海塞矩阵阈值,高一点速度会快些vector<KeyPoint> keyPoint1,keyPoint2; surfDetector.detect(image1,keyPoint1); surfDetector.detect(image2,keyPoint2); //特征点描述,为下边的特征点匹配做准备 SurfDescriptorExtractor SurfDescriptor; Mat imageDesc1,imageDesc2; SurfDescriptor.compute(image1,keyPoint1,imageDesc1); SurfDescriptor.compute(image2,keyPoint2,imageDesc2); //获得匹配特征点,并提取最优配对 FlannBasedMatcher matcher; vector<DMatch> matchePoints; matcher.match(imageDesc1,imageDesc2,matchePoints,Mat()); sort(matchePoints.begin(),matchePoints.end()); //特征点排序 //获取排在前N个的最优匹配特征点 vector<Point2f> imagePoints1,imagePoints2; for(int i=0; i<(int)(matchePoints.size()*0.25); i++) { imagePoints1.push_back(keyPoint1[matchePoints[i].queryIdx].pt); imagePoints2.push_back(keyPoint2[matchePoints[i].trainIdx].pt); } //获取图像1到图像2的投影映射矩阵 尺寸为3*3 Mat homo=findHomography(imagePoints1,imagePoints2,CV_RANSAC); //cout<<"变换矩阵为:\n"<<homo<<endl<<endl; //输出映射矩阵 //图像配准 Mat imageTransform1,imgpeizhun,imgerzhi; warpPerspective(image01,imageTransform1,homo,Size(image02.cols,image02.rows)); //imshow("经过透视矩阵变换后",imageTransform1); absdiff(image02, imageTransform1, imgpeizhun);//imshow("配准diff", imgpeizhun); threshold(imgpeizhun, imgerzhi, 50, 255.0 , CV_THRESH_BINARY);//imshow("配准二值化", imgerzhi);//输出所需时间time0 = ((double)getTickCount()-time0)/getTickFrequency();cout<<1/time0<<endl;Mat temp,image02temp;float m_BiLi = 0.9;image02temp = image02.clone();cvtColor(imgerzhi,temp,CV_RGB2GRAY); //检索连通域Mat se=getStructuringElement(MORPH_RECT, Size(5,5));morphologyEx(temp, temp, MORPH_DILATE, se);vector<vector<Point>> contours;findContours(temp, contours, RETR_EXTERNAL, CHAIN_APPROX_NONE);if (contours.size()<1){continue;}for (int k = 0; k < contours.size(); k++){Rect bomen = boundingRect(contours[k]);//省略由于配准带来的边缘无效信息if (bomen.x > image02temp.cols * (1 - m_BiLi) && bomen.y > image02temp.rows * (1 - m_BiLi) && bomen.x + bomen.width < image02temp.cols * m_BiLi && bomen.y + bomen.height < image02temp.rows * m_BiLi){rectangle(image02temp, bomen, Scalar(255,0,255), 2, 8, 0);}}/*for (int i = 50; i < image02.rows - 100; i++){for (int j = 50; j < image02.cols - 100; j++){uchar pixel = temp.at<uchar>(i,j);if (pixel == 255){Rect bomen(j-7, i-7, 14, 14);rectangle(image02, bomen, Scalar(255,255,255),1,8,0);}}}*/imshow("检测与跟踪",image02temp);waitKey(20); }
}
检测远处运动的车辆
surf消除误匹配点
int surf2(Mat image01, Mat image02)
{Mat image1,image2; image1=image01.clone(); image2=image02.clone(); //提取特征点 SurfFeatureDetector surfDetector(2000); //hessianThreshold,海塞矩阵阈值,并不是限定特征点的个数 vector<KeyPoint> keyPoint1,keyPoint2; surfDetector.detect(image1,keyPoint1); surfDetector.detect(image2,keyPoint2); //绘制特征点 drawKeypoints(image1,keyPoint1,image1,Scalar::all(-1),DrawMatchesFlags::DEFAULT); drawKeypoints(image2,keyPoint2,image2,Scalar::all(-1),DrawMatchesFlags::DRAW_RICH_KEYPOINTS); /* imshow("KeyPoints of image1",image1); imshow("KeyPoints of image2",image2); */ //特征点描述,为下边的特征点匹配做准备 SurfDescriptorExtractor SurfDescriptor; Mat imageDesc1,imageDesc2; SurfDescriptor.compute(image1,keyPoint1,imageDesc1); SurfDescriptor.compute(image2,keyPoint2,imageDesc2); //特征点匹配并显示匹配结果 //BruteForceMatcher<L2<float>> matcher; FlannBasedMatcher matcher; vector<DMatch> matchePoints; matcher.match(imageDesc1,imageDesc2,matchePoints,Mat()); //提取强特征点 double minMatch=1; double maxMatch=0; for(int i=0;i<matchePoints.size();i++) { //匹配值最大最小值获取 minMatch=minMatch>matchePoints[i].distance?matchePoints[i].distance:minMatch; maxMatch=maxMatch<matchePoints[i].distance?matchePoints[i].distance:maxMatch; } //最大最小值输出 cout<<"最佳匹配值是: "<<minMatch<<endl; cout<<"最差匹配值是: "<<maxMatch<<endl; //获取排在前边的几个最优匹配结果 vector<DMatch> goodMatchePoints; for(int i=0;i<matchePoints.size();i++) { if(matchePoints[i].distance<minMatch+(maxMatch-minMatch)/2) { goodMatchePoints.push_back(matchePoints[i]); } } //绘制最优匹配点 Mat imageOutput; drawMatches(image01,keyPoint1,image02,keyPoint2,goodMatchePoints,imageOutput,Scalar::all(-1), Scalar::all(-1),vector<char>(),DrawMatchesFlags::NOT_DRAW_SINGLE_POINTS); imwrite("匹配图.jpg",imageOutput);return 0;
}
【opencv】动态背景下运动目标检测 SURF配准差分相关推荐
- 基于opencV的动态背景下运动目标检测及跟踪(修改版)
基于openCV的动态背景下的运动目标检测 from: http://www.mianfeiwendang.com/doc/89c6692a222a84b2ced0d502/1 摘要:介绍在动态背景下 ...
- 静态背景下运动目标检测
前言:运动对象常用在视频监控领域,目的是从序列图像中将变化区域从背景图像中提取出来,运动区域的有效检测对目标分类.跟踪.行为理解等后期处理非常重要.根据摄像机与运动目标之间的关系可分为静态背景下的运动 ...
- 静态背景下运动目标检测 matlab_动态拉伸、静态拉伸你做对了么?
↑ 点击图片获得JOINFIT健身装备 ↑ "运动乃精神之砥柱,思想之动力!" -----马库斯·图留斯·西塞罗 拉伸对于很多肌友都不算陌生,人多人都喜欢在开始训练前做一些拉伸.但 ...
- 针对无人机航拍视频中动态背景下的目标检测
目录 目录 传统目标检测技术 传统目标检测技术 1. 帧间差分 通过连续两帧相同位置像素点间的灰度差来确定目标移动.但只适用于静态背景和目标单一条件的目标检测.仅适用于无人机悬停状态下的目标检测. 2 ...
- 动态背景下目标跟踪算法目标跟踪
1.问题描述: 动态背景下目标跟踪算法目标跟踪 2.部分程序: % Programming By ALireza Fasih % Email: ar_fasih@yahoo.com % Please ...
- 【OpenCV】 车辆识别 运动目标检测
目录 一:车辆识别 运动目标检测 二:车辆识别实现 超详细步骤解析 步骤一:灰度化处理 步骤二:帧差处理 步骤三:二值化处理 步骤四:图像降噪 4-1 腐蚀处理 目的 去除白色噪点 4-2 膨胀处理 ...
- 计算机视觉与深度学习 | 动态背景下的前景目标提取
################################################ 博主github:https://github.com/MichaelBeechan 博主CSDN:h ...
- 动态背景下的运动目标检测
各种目标检测方法介绍(懒人可以直接略过) 目标检测是一个老话题了,在很多算法当中都有它的身影.目标检测要做的就两件事:检测当前图片中有没有目标?如果有的话,在哪?按照先验知识和背景运动来划分的话,目标 ...
- opencv自然背景下交通标志形状分类c++代码_前端革命时刻:前端代码是怎样智能生成的-图像分离篇
作为阿里经济体前端委员会四大技术方向之一,前端智能化项目经历了 2019 双十一的阶段性考验,交出了不错的答卷,天猫淘宝双十一会场新增模块 79.34% 的线上代码由前端智能化项目自动生成.在此期间研 ...
最新文章
- mongodb报错:E QUERY [thread1] SyntaxError: missing ; before statement @(shell):1:4
- 聊一聊Elasticsearch和MySQL的常用操作
- 用Mesos分布式架构进行工作
- java工作流引擎Jflow流程事件和流程节点事件设置
- linux下日志rorate,[转载]linux下日志分割logrotate 设置和理解
- diff算法_[VUE学习]徒手撸一个虚拟dom+diff算法
- 广联达深思6代锁_广联达|新版加密锁驱动使用文档
- IO、NIO和AIO的区别
- ORACLE11g数据库安装-刘建-专题视频课程
- 如何在Linux系统上装jdk1.8以及环境配置的方法步骤
- AnaConda下载安装完了只有anaconda prompt问题
- python游戏开发实战:黑客帝国特效
- 【FPGA 学习笔记】sof文件和jic文件的区别,程序固化(将sof文件装换位jic文件)
- Shopify卖家:EDM营销就要搭配SaleSmartly,轻松搞定转化率
- 在阿里云OSS,如何更好搭建自己的云储存
- 稍微好看点的登录页面html
- 转载国内私募界知名的操盘手刘文文的访问,收益匪浅
- Varnish反向代理服务器和 HTTP 加速器
- 火山引擎联合中国信通院发布《边缘云总体技术与测试要求》
- android不设密码装证书,安卓4.0不设置锁屏安装根证书教程