树莓派摄像头 C++ OpenCV YoloV3 实现实时目标检测
本文将实现树莓派摄像头 C++ OpenCV YoloV3 实现实时目标检测,我们会先实现树莓派对视频文件的逐帧检测来验证算法流程,成功后,再接入摄像头进行实时目标检测。
先声明一下笔者的主要软硬件配置:
树莓派4B 4GB内存
CSI 摄像头
Ubuntu 20.04
OpenCV 的安装
不多讲,参考 Ubuntu 18.04 安装OpenCV C++ 。
准备YoloV3模型权重文件和视频文件
模型配置文件和权重、COCO数据集名称文件
我们先下载作者官方发布的 YoloV3 模型配置文件、权重文件:
wget https://pjreddie.com/media/files/yolov3.weights
wget https://github.com/pjreddie/darknet/blob/master/cfg/yolov3.cfg?raw=true -O ./yolov3.cfg
上面是比较大的网络,由于我们的树莓派算力比较一遍,所以建议使用轻量型的网络 yolov3-tiny:
wget https://pjreddie.com/media/files/yolov3-tiny.weights
wget https://github.com/pjreddie/darknet/blob/master/cfg/yolov3-tiny.cfg?raw=true -O ./yolov3-tiny.cfg
另外,由于模型权重是在 COCO 数据集上进行预训练的,所以我们还要准备 COCO 的类别名称文件,方便在模型输出检测结果后进行后处理,将类别显示在检测结果框上。
wget https://github.com/pjreddie/darknet/blob/master/data/coco.names?raw=true -O ./coco.names
注:如果上面的 github 中的配置文件在命令行下载的比较慢的话,可以直接去 github 网页复制粘贴下来。
准备视频文件
我们会先对一个视频文件进行逐帧检测来验证算法的流程,在之后再使用摄像头进行实时检测。
我们直接通过 you-get 工具去B站下载视频文件并改个名:
pip install you-get
you-get https://www.bilibili.com/video/av32184680
rm fileName.cmt.xml
mv fileName.mp4 demo.mp4
如果是 flv 文件,可以用 ffmpeg 转为 mp4 文件:
ffmpeg -i input.flv output.mp4
视频文件的检测
一切准备就绪我们开始先测试一下视频文件的检测,我们先讲解一遍代码,在最后会给出整个源码。
1 初始化参数
YOLOv3算法的预测结果就是边界框。每一个边界框都旁随着一个置信值。第一阶段中,全部低于置信度阀值的都会排除掉。 对剩余的边界框执行非最大抑制算法,以去除重叠的边界框。非最大抑制由一个参数 nmsThrehold 控制。读者可以尝试改变这个数值,观察输出的边界框的改变。 接下来,设置输入图片的宽度inpWidth 和高度 inpHeight 。这里设置为416。如果想要更快的速度,可以把宽度和高度设置为320。如果想要更准确的结果,可改变为608。
#include <iostream>
#include <fstream>
#include <vector>
#include <string>
#include <opencv.hpp>
using namespace std;
float confThreshold = 0.5;
float nmsThreshold = 0.4;
int inpWidth = 416;
int inpHeight = 416;
2 读取模型和COCO类别名
接下来我们读入COCO 类别名并存入 classes 容器。并加载模型与权重文件 yolov3-tiny.cfg 和 yolov3-tiny.weights 。这里用到的几个文件就是我们刚才下载好的,读者需要改为自己的路径。最后把DNN的后端设置为OpenCV,目标设置CPU。这里我们树莓派没有GPU等加速推理硬件,就用CPU,如果有GPU,可改为OpenCL、CUDA等
vector<string> classes;
string classesFile = "./coco.names";
ifstream ifs(classesFile.c_str());
string line;
while(getline(ifs,line))classes.push_back(line);
cv::String modelConfiguration = "./yolov3-tiny.cfg";
cv::String modelWeights = "./yolov3-tiny.weights";
cv::dnn::Net net = cv::dnn::readNetFromDarknet(modelConfiguration,modelWeights);
net.setPreferableBackend(cv::dnn::DNN_BACKEND_OPENCV);
net.setPreferableTarget(cv::dnn::DNN_TARGET_OPENCL);
3 读取输入
这里我们先读入下载好的视频文件,一会儿再使用本地摄像头测试。这里如果是树莓派外接显示器的话可以用创建GUI窗口来查看,但是我们通常是命令行SSH连接树莓派,这时我们就直接将每一帧的检测结果保存为图像文件查看:
string str, outputFile;
cv::VideoCapture cap("./demo.mp4");
cv::VideoWriter video;
cv::Mat frame,blob;
4 循环处理每一帧
OpenCV中,输入到神经网络的图像需要以一种叫 bolb 的格式保存。 读取了输入图片或者视频流的一帧图像后,这帧图像需要经过bolbFromImage() 函数处理为神经网络的输入类型 bolb。在这个过程中,图像像素以一个 1/255 的比例因子,被缩放到0到1之间。同时,图像在不裁剪的情况下,大小调整到 416x416。注意我们没有降低图像平均值,因此传递 [0,0,0] 到函数的平均值输入,保持swapRB 参数到默认值1。 输出的bolb传递到网络,经过网络正向处理,网络输出了所预测到的一个边界框清单。这些边界框通过后处理,滤除了低置信值的。我们随后再详细的说明后处理的步骤。我们在每一帧的左上方打印出了推断时间。伴随着最后的边界框的完成,图像保存到硬盘中,之后可以作为图像输入或者通过 VideoWriter 作为视频流输入。
while(cv::waitKey(1)<0){
cap>>frame;
if(frame.empty()){
break;
}
cv::dnn::blobFromImage(frame,blob,1/255.0,cv::Size(inpWidth,inpHeight));
net.setInput(blob);
vector<cv::Mat> outs;
net.forward(outs,getOutputNames(net));
postprocess(frame,outs);
vector<double> layersTimes;
double freq = cv::getTickFrequency()/1000;
double t=net.getPerfProfile(layersTimes)/freq;
string label = cv::format("Infercence time for a frame:%.2f ms",t);
cv::putText(frame,label,cv::Point(0,15),cv::FONT_HERSHEY_SIMPLEX,0.5,cv::Scalar(0,255,255));
cout << "Frame: " << frame_cnt++ << ", time: " << t << "ms" << endl;
cv::imwrite("output.jpg",frame);
}
5-1 得到输出层的名字
第五步我们给出几个用到的函数的实现。
OpenCV 的网络类中的前向功能需要结束层,直到它在网络中运行。因为我们需要运行整个网络,所以我们需要识别网络中的最后一层。我们通过使用 getUnconnectedOutLayers() 获得未连接的输出层的名字,该层基本就是网络的最后层。然后我们运行前向网络,得到输出,如前面的代码片段 net.forward(getOutputsNames(net)) 。
vector<cv::String> getOutputNames(const cv::dnn::Net& net){
static vector<cv::String> names;
if(names.empty()){
vector<int> outLayers = net.getUnconnectedOutLayers();
vector<cv::String> layersNames = net.getLayerNames();
names.resize(outLayers.size());
for(size_t i =0;i<outLayers.size();i++){
names[i] = layersNames[outLayers[i]-1];
}
}
return names;
}
5-2 后处理
网络输出的每个边界框都分别由一个包含着类别名字和5个元素的向量表示。 头四个元素代表center_x, center_y, width, height 。第五个元素表示包含着目标的边界框的置信度。 其余的元素是和每个类别(如目标种类)有关的置信度。边界框分配给最高分数对应的那一种类。 一个边界框的最高分数也叫做它的置信度 confidence 。如果边界框的置信度低于规定的阀值,算法上不再处理这个边界框。 置信度大于或等于置信度阀值的边界框,将进行非最大抑制。这会减少重叠的边界框数目。
void postprocess(cv::Mat& frame,const vector<cv::Mat>& outs){
vector<int> classIds;
vector<float> confidences;
vector<cv::Rect> boxes;
for(size_t i=0;i<outs.size();i++){
float* data = (float*)outs[i].data;
for(int j=0;j<outs[i].rows;j++,data+=outs[i].cols){
cv::Mat scores = outs[i].row(j).colRange(5,outs[i].cols);
cv::Point classIdPoint;
double confidence;
cv::minMaxLoc(scores,0,&confidence,0,&classIdPoint);
if(confidence>confThreshold){
int centerX = (int)(data[0]*frame.cols);
int centerY = (int)(data[1]*frame.rows);
int width = (int)(data[2]*frame.cols);
int height = (int)(data[3]*frame.rows);
int left = centerX-width/2;
int top = centerY-height/2;
classIds.push_back(classIdPoint.x);
confidences.push_back((float)confidence);
boxes.push_back(cv::Rect(left, top, width, height));
}
}
}
vector<int> indices;
cv::dnn::NMSBoxes(boxes,confidences,confThreshold,nmsThreshold,indices);
for(size_t i=0;i<indices.size();i++){
int idx = indices[i];
cv::Rect box = boxes[idx];
drawPred(classIds[idx],confidences[idx],box.x,box.y,
box.x+box.width,box.y+box.height,frame);
}
}
5-3 画出边界框
最后,经过非最大抑制后,得到了边界框。我们把边界框在输入帧上画出,并标出种类名和置信值。
void drawPred(int classId,float conf,int left,int top,int right,int bottom,cv::Mat& frame){
cv::rectangle(frame,cv::Point(left,top),cv::Point(right,bottom),cv::Scalar(255,178,50),3);
string label = cv::format("%.2f",conf);
if(!classes.empty()){
CV_Assert(classId < (int)classes.size());
label = classes[classId]+":"+label;
}
int baseLine;
cv::Size labelSize = cv::getTextSize(label,cv::FONT_HERSHEY_SIMPLEX,0.5,1,&baseLine);
top = max(top,labelSize.height);
cv::rectangle(frame,cv::Point(left,top-round(1.5*labelSize.height)),cv::Point(left+round(1.5*labelSize.width),top+baseLine),cv::Scalar(255,255,255),cv::FILLED);
cv::putText(frame, label,cv::Point(left, top), cv::FONT_HERSHEY_SIMPLEX, 0.75,cv::Scalar(0, 0, 0), 1);
}
文件全部源码在文末。
写好之后我们编译执行即可,关于 OpenCV 众多头文件包含、链接库链接时、运行时的链接,对初学者来说可能会遇到一些问题,可以参考:
Linux下编译、链接、加载运行C++ OpenCV的两种方式及常见问题的解决
Linux下C/C++程序编译链接加载过程中的常见问题及解决方法
可以从左上角和标准输出看到,每帧的检测时间大概在 280ms,速度还可以,精度大体也不错。但是由于模型较小,性能受限,对于一些边缘小物体会有误差,如上图中右侧的小车。
树莓派摄像头实时检测
树莓派摄像头调试参考:树莓派摄像头基础配置及测试 。
在视频文件的检测顺利完成后,实时树莓派的检测就很简单了,只需要将读取输入部分从视频文件改为本地摄像头即可。
主要就是这一行修改:
cv::VideoCapture cap(0);
另外,我们需要设置一些 OpenCV 读取摄像头输入的尺寸大小,否则笔者亲测是有一些小bug:
cap.set(cv::CAP_PROP_FRAME_WIDTH, 640);
cap.set(cv::CAP_PROP_FRAME_HEIGHT, 480);
在笔者实验室中实测,速度和精度也都还可以。
全部代码
这里给出摄像头实时测试的全部代码
#include <iostream>
#include <fstream>
#include <vector>
#include <string>
#include <opencv.hpp>
using namespace std;
vector<string> classes;
float confThreshold = 0.5;
float nmsThreshold = 0.4;
int inpWidth = 320;
int inpHeight = 320;
void postprocess(cv::Mat& frame,const vector<cv::Mat>& out);
void drawPred(int classId,float conf,int left,int top,int right,int bottom,cv::Mat& frame);
vector<cv::String> getOutputNames(const cv::dnn::Net& net);
int main(int argc, char const *argv[])
{
string classesFile = "./coco.names";
ifstream ifs(classesFile.c_str());
string line;
while(getline(ifs,line))classes.push_back(line);
cv::String modelConfiguration = "./model_file/yolov3-tiny.cfg";
cv::String modelWeights = "./model_file/yolov3-tiny.weights";
cv::dnn::Net net = cv::dnn::readNetFromDarknet(modelConfiguration,modelWeights);
net.setPreferableBackend(cv::dnn::DNN_BACKEND_OPENCV);
net.setPreferableTarget(cv::dnn::DNN_TARGET_CPU);
string str, outputFile;
cv::Mat frame, blob;
cv::VideoCapture cap(0);
cap.set(cv::CAP_PROP_FRAME_WIDTH, 640);
cap.set(cv::CAP_PROP_FRAME_HEIGHT, 480);
int frame_cnt = 0;
while(cv::waitKey(1)<0){
cap >> frame;
if(frame.empty()){
break;
}
cv::dnn::blobFromImage(frame,blob,1/255.0,cv::Size(inpWidth,inpHeight));
net.setInput(blob);
vector<cv::Mat> outs;
net.forward(outs,getOutputNames(net));
postprocess(frame,outs);
vector<double> layersTimes;
double freq = cv::getTickFrequency()/1000;
double t=net.getPerfProfile(layersTimes)/freq;
string label = cv::format("Infercence time for a frame:%.2f ms",t);
cv::putText(frame,label,cv::Point(0,15),cv::FONT_HERSHEY_SIMPLEX,0.5,cv::Scalar(0,255,255));
cv::Mat detecteFrame;
frame.convertTo(detecteFrame,CV_8U);
cout << "Frame: " << frame_cnt++ << ", time: " << t << "ms" << endl;
cv::imwrite("output.jpg",frame);
}
cap.release();
return 0;
}
void postprocess(cv::Mat& frame,const vector<cv::Mat>& outs){
vector<int> classIds;
vector<float> confidences;
vector<cv::Rect> boxes;
for(size_t i=0;i<outs.size();i++){
float* data = (float*)outs[i].data;
for(int j=0;j<outs[i].rows;j++,data+=outs[i].cols){
cv::Mat scores = outs[i].row(j).colRange(5,outs[i].cols);
cv::Point classIdPoint;
double confidence;
cv::minMaxLoc(scores,0,&confidence,0,&classIdPoint);
if(confidence>confThreshold){
int centerX = (int)(data[0]*frame.cols);
int centerY = (int)(data[1]*frame.rows);
int width = (int)(data[2]*frame.cols);
int height = (int)(data[3]*frame.rows);
int left = centerX-width/2;
int top = centerY-height/2;
classIds.push_back(classIdPoint.x);
confidences.push_back((float)confidence);
boxes.push_back(cv::Rect(left, top, width, height));
}
}
}
vector<int> indices;
cv::dnn::NMSBoxes(boxes,confidences,confThreshold,nmsThreshold,indices);
for(size_t i=0;i<indices.size();i++){
int idx = indices[i];
cv::Rect box = boxes[idx];
drawPred(classIds[idx],confidences[idx],box.x,box.y,
box.x+box.width,box.y+box.height,frame);
}
}
void drawPred(int classId,float conf,int left,int top,int right,int bottom,cv::Mat& frame){
cv::rectangle(frame,cv::Point(left,top),cv::Point(right,bottom),cv::Scalar(255,178,50),3);
string label = cv::format("%.2f",conf);
if(!classes.empty()){
CV_Assert(classId < (int)classes.size());
label = classes[classId]+":"+label;
}
int baseLine;
cv::Size labelSize = cv::getTextSize(label,cv::FONT_HERSHEY_SIMPLEX,0.5,1,&baseLine);
top = max(top,labelSize.height);
cv::rectangle(frame,cv::Point(left,top-round(1.5*labelSize.height)),cv::Point(left+round(1.5*labelSize.width),top+baseLine),cv::Scalar(255,255,255),cv::FILLED);
cv::putText(frame, label,cv::Point(left, top), cv::FONT_HERSHEY_SIMPLEX, 0.75,cv::Scalar(0, 0, 0), 1);
}
vector<cv::String> getOutputNames(const cv::dnn::Net& net){
static vector<cv::String> names;
if(names.empty()){
vector<int> outLayers = net.getUnconnectedOutLayers();
vector<cv::String> layersNames = net.getLayerNames();
names.resize(outLayers.size());
for(size_t i =0;i<outLayers.size();i++){
names[i] = layersNames[outLayers[i]-1];
}
}
return names;
}
如有错误或遗漏,欢迎留言指正。
Ref:
https://blog.csdn.net/cuma2369/article/details/107666559
https://ryanadex.github.io/2019/01/15/opencv-yolov3/
|