2012-08-31 14:08:52 +08:00
|
|
|
#include "opencv2/objdetect/objdetect.hpp"
|
|
|
|
#include "opencv2/highgui/highgui.hpp"
|
|
|
|
#include "opencv2/imgproc/imgproc.hpp"
|
|
|
|
#include "opencv2/ocl/ocl.hpp"
|
2013-04-08 02:45:38 +08:00
|
|
|
|
|
|
|
#include "opencv2/highgui/highgui_c.h"
|
|
|
|
|
2012-08-31 14:08:52 +08:00
|
|
|
#include <iostream>
|
|
|
|
#include <stdio.h>
|
|
|
|
|
2013-04-12 18:12:12 +08:00
|
|
|
|
2012-08-31 14:08:52 +08:00
|
|
|
using namespace std;
|
|
|
|
using namespace cv;
|
2013-06-04 22:31:51 +08:00
|
|
|
#define LOOP_NUM 10
|
2013-05-24 15:52:33 +08:00
|
|
|
|
|
|
|
const static Scalar colors[] = { CV_RGB(0,0,255),
|
2013-06-19 16:36:23 +08:00
|
|
|
CV_RGB(0,128,255),
|
|
|
|
CV_RGB(0,255,255),
|
|
|
|
CV_RGB(0,255,0),
|
|
|
|
CV_RGB(255,128,0),
|
|
|
|
CV_RGB(255,255,0),
|
|
|
|
CV_RGB(255,0,0),
|
|
|
|
CV_RGB(255,0,255)
|
|
|
|
} ;
|
|
|
|
|
2012-08-31 14:08:52 +08:00
|
|
|
|
2013-05-24 15:52:33 +08:00
|
|
|
int64 work_begin = 0;
|
|
|
|
int64 work_end = 0;
|
2013-06-19 16:36:23 +08:00
|
|
|
string outputName;
|
2013-05-24 15:52:33 +08:00
|
|
|
|
2013-06-04 22:31:51 +08:00
|
|
|
static void workBegin()
|
|
|
|
{
|
2013-05-24 15:52:33 +08:00
|
|
|
work_begin = getTickCount();
|
|
|
|
}
|
|
|
|
static void workEnd()
|
2012-08-31 14:08:52 +08:00
|
|
|
{
|
2013-05-24 15:52:33 +08:00
|
|
|
work_end += (getTickCount() - work_begin);
|
2012-08-31 14:08:52 +08:00
|
|
|
}
|
2013-04-12 18:12:12 +08:00
|
|
|
|
2013-06-19 16:36:23 +08:00
|
|
|
static double getTime()
|
|
|
|
{
|
2013-05-24 15:52:33 +08:00
|
|
|
return work_end /((double)cvGetTickFrequency() * 1000.);
|
|
|
|
}
|
|
|
|
|
2013-06-04 22:31:51 +08:00
|
|
|
void detect( Mat& img, vector<Rect>& faces,
|
2013-07-22 17:53:56 +08:00
|
|
|
ocl::OclCascadeClassifier& cascade,
|
2013-06-19 16:36:23 +08:00
|
|
|
double scale, bool calTime);
|
|
|
|
|
2012-08-31 14:08:52 +08:00
|
|
|
|
2013-06-04 22:31:51 +08:00
|
|
|
void detectCPU( Mat& img, vector<Rect>& faces,
|
2013-06-19 16:36:23 +08:00
|
|
|
CascadeClassifier& cascade,
|
|
|
|
double scale, bool calTime);
|
2013-05-24 15:52:33 +08:00
|
|
|
|
|
|
|
void Draw(Mat& img, vector<Rect>& faces, double scale);
|
|
|
|
|
2013-06-19 16:36:23 +08:00
|
|
|
|
2013-05-24 15:52:33 +08:00
|
|
|
// This function test if gpu_rst matches cpu_rst.
|
|
|
|
// If the two vectors are not equal, it will return the difference in vector size
|
|
|
|
// Else if will return (total diff of each cpu and gpu rects covered pixels)/(total cpu rects covered pixels)
|
2013-06-19 16:36:23 +08:00
|
|
|
double checkRectSimilarity(Size sz, vector<Rect>& cpu_rst, vector<Rect>& gpu_rst);
|
|
|
|
|
2012-08-31 14:08:52 +08:00
|
|
|
|
|
|
|
int main( int argc, const char** argv )
|
|
|
|
{
|
2013-05-24 15:52:33 +08:00
|
|
|
const char* keys =
|
2013-07-22 17:53:56 +08:00
|
|
|
"{ h help | false | print help message }"
|
|
|
|
"{ i input | | specify input image }"
|
|
|
|
"{ t template | haarcascade_frontalface_alt.xml |"
|
2013-06-19 16:36:23 +08:00
|
|
|
" specify template file path }"
|
2013-07-22 17:53:56 +08:00
|
|
|
"{ c scale | 1.0 | scale image }"
|
|
|
|
"{ s use_cpu | false | use cpu or gpu to process the image }"
|
|
|
|
"{ o output | facedetect_output.jpg |"
|
2013-06-19 16:36:23 +08:00
|
|
|
" specify output image save path(only works when input is images) }";
|
2013-05-24 15:52:33 +08:00
|
|
|
|
|
|
|
CommandLineParser cmd(argc, argv, keys);
|
|
|
|
if (cmd.get<bool>("help"))
|
2012-10-17 07:18:30 +08:00
|
|
|
{
|
2013-05-24 15:52:33 +08:00
|
|
|
cout << "Avaible options:" << endl;
|
|
|
|
return 0;
|
2012-10-17 07:18:30 +08:00
|
|
|
}
|
2013-05-24 15:52:33 +08:00
|
|
|
CvCapture* capture = 0;
|
|
|
|
Mat frame, frameCopy, image;
|
2012-10-17 07:18:30 +08:00
|
|
|
|
2013-05-24 15:52:33 +08:00
|
|
|
bool useCPU = cmd.get<bool>("s");
|
|
|
|
string inputName = cmd.get<string>("i");
|
2013-06-19 16:36:23 +08:00
|
|
|
outputName = cmd.get<string>("o");
|
2013-05-24 15:52:33 +08:00
|
|
|
string cascadeName = cmd.get<string>("t");
|
|
|
|
double scale = cmd.get<double>("c");
|
2013-07-22 17:53:56 +08:00
|
|
|
ocl::OclCascadeClassifier cascade;
|
2013-05-24 15:52:33 +08:00
|
|
|
CascadeClassifier cpu_cascade;
|
|
|
|
|
|
|
|
if( !cascade.load( cascadeName ) || !cpu_cascade.load(cascadeName) )
|
2012-10-17 07:18:30 +08:00
|
|
|
{
|
|
|
|
cerr << "ERROR: Could not load classifier cascade" << endl;
|
|
|
|
return -1;
|
|
|
|
}
|
|
|
|
|
2013-05-24 15:52:33 +08:00
|
|
|
if( inputName.empty() )
|
2012-10-17 07:18:30 +08:00
|
|
|
{
|
2013-05-24 15:52:33 +08:00
|
|
|
capture = cvCaptureFromCAM(0);
|
2013-06-04 22:31:51 +08:00
|
|
|
if(!capture)
|
2013-05-24 15:52:33 +08:00
|
|
|
cout << "Capture from CAM 0 didn't work" << endl;
|
2012-10-17 07:18:30 +08:00
|
|
|
}
|
|
|
|
else if( inputName.size() )
|
|
|
|
{
|
|
|
|
image = imread( inputName, 1 );
|
|
|
|
if( image.empty() )
|
|
|
|
{
|
|
|
|
capture = cvCaptureFromAVI( inputName.c_str() );
|
2013-06-04 22:31:51 +08:00
|
|
|
if(!capture)
|
2013-05-24 15:52:33 +08:00
|
|
|
cout << "Capture from AVI didn't work" << endl;
|
|
|
|
return -1;
|
2012-10-17 07:18:30 +08:00
|
|
|
}
|
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
|
|
|
image = imread( "lena.jpg", 1 );
|
2013-06-04 22:31:51 +08:00
|
|
|
if(image.empty())
|
2013-05-24 15:52:33 +08:00
|
|
|
cout << "Couldn't read lena.jpg" << endl;
|
|
|
|
return -1;
|
2012-10-17 07:18:30 +08:00
|
|
|
}
|
|
|
|
|
2013-06-19 16:36:23 +08:00
|
|
|
|
2012-10-17 07:18:30 +08:00
|
|
|
cvNamedWindow( "result", 1 );
|
2013-06-19 16:36:23 +08:00
|
|
|
vector<ocl::Info> oclinfo;
|
|
|
|
int devnums = ocl::getDevice(oclinfo);
|
2013-05-24 15:52:33 +08:00
|
|
|
if( devnums < 1 )
|
2012-10-17 07:18:30 +08:00
|
|
|
{
|
|
|
|
std::cout << "no device found\n";
|
|
|
|
return -1;
|
|
|
|
}
|
|
|
|
//if you want to use undefault device, set it here
|
|
|
|
//setDevice(oclinfo[0]);
|
2013-05-24 15:52:33 +08:00
|
|
|
ocl::setBinpath("./");
|
2012-10-17 07:18:30 +08:00
|
|
|
if( capture )
|
|
|
|
{
|
|
|
|
cout << "In capture ..." << endl;
|
|
|
|
for(;;)
|
|
|
|
{
|
|
|
|
IplImage* iplImg = cvQueryFrame( capture );
|
2013-03-30 04:47:22 +08:00
|
|
|
frame = cv::cvarrToMat(iplImg);
|
2013-05-24 15:52:33 +08:00
|
|
|
vector<Rect> faces;
|
2012-10-17 07:18:30 +08:00
|
|
|
if( frame.empty() )
|
|
|
|
break;
|
|
|
|
if( iplImg->origin == IPL_ORIGIN_TL )
|
|
|
|
frame.copyTo( frameCopy );
|
|
|
|
else
|
|
|
|
flip( frame, frameCopy, 0 );
|
2013-06-19 16:36:23 +08:00
|
|
|
if(useCPU)
|
|
|
|
{
|
2013-05-24 15:52:33 +08:00
|
|
|
detectCPU(frameCopy, faces, cpu_cascade, scale, false);
|
|
|
|
}
|
2013-06-19 16:36:23 +08:00
|
|
|
else
|
|
|
|
{
|
2013-06-04 22:31:51 +08:00
|
|
|
detect(frameCopy, faces, cascade, scale, false);
|
2013-05-24 15:52:33 +08:00
|
|
|
}
|
|
|
|
Draw(frameCopy, faces, scale);
|
2012-10-17 07:18:30 +08:00
|
|
|
if( waitKey( 10 ) >= 0 )
|
|
|
|
goto _cleanup_;
|
|
|
|
}
|
|
|
|
|
2013-06-19 16:36:23 +08:00
|
|
|
|
2012-10-17 07:18:30 +08:00
|
|
|
waitKey(0);
|
2012-08-31 14:08:52 +08:00
|
|
|
|
2013-06-19 16:36:23 +08:00
|
|
|
|
2012-08-31 14:08:52 +08:00
|
|
|
_cleanup_:
|
2012-10-17 07:18:30 +08:00
|
|
|
cvReleaseCapture( &capture );
|
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
|
|
|
cout << "In image read" << endl;
|
2013-05-24 15:52:33 +08:00
|
|
|
vector<Rect> faces;
|
|
|
|
vector<Rect> ref_rst;
|
|
|
|
double accuracy = 0.;
|
2013-06-19 16:36:23 +08:00
|
|
|
for(int i = 0; i <= LOOP_NUM; i ++)
|
2012-10-17 07:18:30 +08:00
|
|
|
{
|
2013-05-24 15:52:33 +08:00
|
|
|
cout << "loop" << i << endl;
|
2013-06-19 16:36:23 +08:00
|
|
|
if(useCPU)
|
|
|
|
{
|
2013-06-04 22:31:51 +08:00
|
|
|
detectCPU(image, faces, cpu_cascade, scale, i==0?false:true);
|
2013-05-24 15:52:33 +08:00
|
|
|
}
|
2013-06-19 16:36:23 +08:00
|
|
|
else
|
|
|
|
{
|
2013-05-24 15:52:33 +08:00
|
|
|
detect(image, faces, cascade, scale, i==0?false:true);
|
2013-06-19 16:36:23 +08:00
|
|
|
if(i == 0)
|
|
|
|
{
|
2013-05-24 15:52:33 +08:00
|
|
|
detectCPU(image, ref_rst, cpu_cascade, scale, false);
|
|
|
|
accuracy = checkRectSimilarity(image.size(), ref_rst, faces);
|
2012-10-17 07:18:30 +08:00
|
|
|
}
|
2013-05-24 15:52:33 +08:00
|
|
|
}
|
|
|
|
if (i == LOOP_NUM)
|
2012-10-17 07:18:30 +08:00
|
|
|
{
|
2013-05-24 15:52:33 +08:00
|
|
|
if (useCPU)
|
|
|
|
cout << "average CPU time (noCamera) : ";
|
|
|
|
else
|
|
|
|
cout << "average GPU time (noCamera) : ";
|
|
|
|
cout << getTime() / LOOP_NUM << " ms" << endl;
|
|
|
|
cout << "accuracy value: " << accuracy <<endl;
|
2012-10-17 07:18:30 +08:00
|
|
|
}
|
|
|
|
}
|
2013-05-24 15:52:33 +08:00
|
|
|
Draw(image, faces, scale);
|
|
|
|
waitKey(0);
|
2012-10-17 07:18:30 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
cvDestroyWindow("result");
|
|
|
|
return 0;
|
2012-08-31 14:08:52 +08:00
|
|
|
}
|
|
|
|
|
2013-06-04 22:31:51 +08:00
|
|
|
void detect( Mat& img, vector<Rect>& faces,
|
2013-07-22 17:53:56 +08:00
|
|
|
ocl::OclCascadeClassifier& cascade,
|
2013-06-19 16:36:23 +08:00
|
|
|
double scale, bool calTime)
|
2012-08-31 14:08:52 +08:00
|
|
|
{
|
2013-06-19 16:36:23 +08:00
|
|
|
ocl::oclMat image(img);
|
|
|
|
ocl::oclMat gray, smallImg( cvRound (img.rows/scale), cvRound(img.cols/scale), CV_8UC1 );
|
2013-05-24 15:52:33 +08:00
|
|
|
if(calTime) workBegin();
|
2013-06-25 18:08:23 +08:00
|
|
|
ocl::cvtColor( image, gray, COLOR_BGR2GRAY );
|
2013-06-19 16:36:23 +08:00
|
|
|
ocl::resize( gray, smallImg, smallImg.size(), 0, 0, INTER_LINEAR );
|
|
|
|
ocl::equalizeHist( smallImg, smallImg );
|
2012-10-17 07:18:30 +08:00
|
|
|
|
2013-05-24 15:52:33 +08:00
|
|
|
cascade.detectMultiScale( smallImg, faces, 1.1,
|
2013-06-19 16:36:23 +08:00
|
|
|
3, 0
|
2013-07-22 17:53:56 +08:00
|
|
|
|CASCADE_SCALE_IMAGE
|
2013-06-19 16:36:23 +08:00
|
|
|
, Size(30,30), Size(0, 0) );
|
2013-05-24 15:52:33 +08:00
|
|
|
if(calTime) workEnd();
|
|
|
|
}
|
|
|
|
|
2013-06-04 22:31:51 +08:00
|
|
|
void detectCPU( Mat& img, vector<Rect>& faces,
|
2013-06-19 16:36:23 +08:00
|
|
|
CascadeClassifier& cascade,
|
|
|
|
double scale, bool calTime)
|
2013-05-24 15:52:33 +08:00
|
|
|
{
|
|
|
|
if(calTime) workBegin();
|
|
|
|
Mat cpu_gray, cpu_smallImg( cvRound (img.rows/scale), cvRound(img.cols/scale), CV_8UC1 );
|
2013-07-22 17:53:56 +08:00
|
|
|
cvtColor(img, cpu_gray, COLOR_BGR2GRAY);
|
2013-05-24 15:52:33 +08:00
|
|
|
resize(cpu_gray, cpu_smallImg, cpu_smallImg.size(), 0, 0, INTER_LINEAR);
|
|
|
|
equalizeHist(cpu_smallImg, cpu_smallImg);
|
|
|
|
cascade.detectMultiScale(cpu_smallImg, faces, 1.1,
|
2013-07-22 17:53:56 +08:00
|
|
|
3, 0 | CASCADE_SCALE_IMAGE,
|
2013-06-19 16:36:23 +08:00
|
|
|
Size(30, 30), Size(0, 0));
|
2013-06-04 22:31:51 +08:00
|
|
|
if(calTime) workEnd();
|
2013-05-24 15:52:33 +08:00
|
|
|
}
|
|
|
|
|
2013-06-19 16:36:23 +08:00
|
|
|
|
2013-05-24 15:52:33 +08:00
|
|
|
void Draw(Mat& img, vector<Rect>& faces, double scale)
|
|
|
|
{
|
|
|
|
int i = 0;
|
2012-10-17 07:18:30 +08:00
|
|
|
for( vector<Rect>::const_iterator r = faces.begin(); r != faces.end(); r++, i++ )
|
|
|
|
{
|
|
|
|
Point center;
|
|
|
|
Scalar color = colors[i%8];
|
|
|
|
int radius;
|
|
|
|
center.x = cvRound((r->x + r->width*0.5)*scale);
|
|
|
|
center.y = cvRound((r->y + r->height*0.5)*scale);
|
|
|
|
radius = cvRound((r->width + r->height)*0.25*scale);
|
|
|
|
circle( img, center, radius, color, 3, 8, 0 );
|
|
|
|
}
|
2013-06-19 16:36:23 +08:00
|
|
|
imwrite( outputName, img );
|
2013-06-28 15:08:39 +08:00
|
|
|
if(abs(scale-1.0)>.001)
|
|
|
|
{
|
2013-06-28 16:23:01 +08:00
|
|
|
resize(img, img, Size((int)(img.cols/scale), (int)(img.rows/scale)));
|
2013-06-28 15:08:39 +08:00
|
|
|
}
|
|
|
|
imshow( "result", img );
|
|
|
|
|
2012-08-31 14:08:52 +08:00
|
|
|
}
|
2013-05-24 15:52:33 +08:00
|
|
|
|
2013-06-19 16:36:23 +08:00
|
|
|
|
|
|
|
double checkRectSimilarity(Size sz, vector<Rect>& ob1, vector<Rect>& ob2)
|
2013-05-24 15:52:33 +08:00
|
|
|
{
|
|
|
|
double final_test_result = 0.0;
|
|
|
|
size_t sz1 = ob1.size();
|
|
|
|
size_t sz2 = ob2.size();
|
|
|
|
|
|
|
|
if(sz1 != sz2)
|
2013-06-19 16:36:23 +08:00
|
|
|
{
|
2013-05-24 15:52:33 +08:00
|
|
|
return sz1 > sz2 ? (double)(sz1 - sz2) : (double)(sz2 - sz1);
|
2013-06-19 16:36:23 +08:00
|
|
|
}
|
2013-05-24 15:52:33 +08:00
|
|
|
else
|
|
|
|
{
|
2013-06-19 16:36:23 +08:00
|
|
|
if(sz1==0 && sz2==0)
|
|
|
|
return 0;
|
|
|
|
Mat cpu_result(sz, CV_8UC1);
|
2013-05-24 15:52:33 +08:00
|
|
|
cpu_result.setTo(0);
|
|
|
|
|
|
|
|
for(vector<Rect>::const_iterator r = ob1.begin(); r != ob1.end(); r++)
|
2013-06-04 22:31:51 +08:00
|
|
|
{
|
2013-06-19 16:36:23 +08:00
|
|
|
Mat cpu_result_roi(cpu_result, *r);
|
2013-05-24 15:52:33 +08:00
|
|
|
cpu_result_roi.setTo(1);
|
|
|
|
cpu_result.copyTo(cpu_result);
|
|
|
|
}
|
2013-06-19 16:36:23 +08:00
|
|
|
int cpu_area = countNonZero(cpu_result > 0);
|
2013-05-24 15:52:33 +08:00
|
|
|
|
|
|
|
|
2013-06-19 16:36:23 +08:00
|
|
|
Mat gpu_result(sz, CV_8UC1);
|
2013-05-24 15:52:33 +08:00
|
|
|
gpu_result.setTo(0);
|
|
|
|
for(vector<Rect>::const_iterator r2 = ob2.begin(); r2 != ob2.end(); r2++)
|
|
|
|
{
|
|
|
|
cv::Mat gpu_result_roi(gpu_result, *r2);
|
|
|
|
gpu_result_roi.setTo(1);
|
|
|
|
gpu_result.copyTo(gpu_result);
|
|
|
|
}
|
|
|
|
|
2013-06-19 16:36:23 +08:00
|
|
|
Mat result_;
|
2013-05-24 15:52:33 +08:00
|
|
|
multiply(cpu_result, gpu_result, result_);
|
2013-06-19 16:36:23 +08:00
|
|
|
int result = countNonZero(result_ > 0);
|
|
|
|
if(cpu_area!=0 && result!=0)
|
|
|
|
final_test_result = 1.0 - (double)result/(double)cpu_area;
|
|
|
|
else if(cpu_area==0 && result!=0)
|
|
|
|
final_test_result = -1;
|
2013-05-24 15:52:33 +08:00
|
|
|
}
|
|
|
|
return final_test_result;
|
2013-07-22 17:53:56 +08:00
|
|
|
}
|