上一篇教程中，我们学习了如何计算轮廓的凸包，其实对一个轮廓而言，可能它的凸包和它本身是重合的，也有可能不是重合的。比如下面左边图像的轮廓本身就是凸包，而右边图像的轮廓则不是。我们可以通过函数bool isContourConvex(InputArray contour)，来判定一个轮廓是否是凸包，是的话返回true，否则false[注意测试的轮廓必须是简单轮廓，没有自交叉之类的]。
对一个非凸包的轮廓而言，它包括一系列的凹陷区域，这些区域称作defect，比如下面手轮廓中，包括6个defect区域。在OpenCV中，我们用下面的结构来定义defect。
struct CvConvexityDefect { CvPoint* start; // 轮廓中defect的起点 CvPoint* end; // 轮廓中defect的终点 CvPoint* depth_point; // defect中到凸包最远的点 float depth; // 最远点和凸包之间的距离};
在OpenCV中，我们通过函数
void convexityDefects(InputArray contour, InputArray convexhull, OutputArray convexityDefects)
得到轮廓的凸包，其中第一个参数和第二个参数是轮廓以及轮廓对应的凸包，注意凸包应该使用vector<int>这样的索引方式表示。第三个参数为返回的defect点集。
下面我们看下检测轮廓defects的代码：
#include "opencv2/highgui/highgui.hpp"
#include "opencv2/imgproc/imgproc.hpp"
#include <iostream>
#include <stdio.h>
#include <stdlib.h>
using namespace cv;
using namespace std;
Mat src; Mat src_gray;
RNG rng(12345);
int main( int argc, char** argv )
{
//装入图像
src = imread("../hand1.jpg", 1 );
//转化为灰度图像
cvtColor( src, src_gray, CV_BGR2GRAY );
//blur( src_gray, src_gray, Size(3,3) );
namedWindow( "image");
imshow( "image", src_gray );
Mat threshold_output;
vector<vector<Point> > contours;
vector<Vec4i> hierarchy;
//得到二值图
threshold( src_gray, threshold_output, 200, 255, THRESH_BINARY );
//查找轮廓
findContours( threshold_output, contours, hierarchy, CV_RETR_TREE, CV_CHAIN_APPROX_SIMPLE, Point(0, 0) );
printf("轮廓数目:%d\n", contours.size());
/// Find the convex hull object for each contour
vector<vector<Point> >hull( contours.size() );
for( int i = 0; i < contours.size(); i++ )
{ convexHull( Mat(contours[i]), hull[i], false ); }
/// Draw contours + hull results
Mat drawing = Mat::zeros( threshold_output.size(), CV_8UC3 );
int area = 0; //轮廓索引
int k = 0;
int i;
for(i = 0; i< contours.size(); i++ )
{
Scalar color1 = Scalar( rng.uniform(0, 255), rng.uniform(0,255), rng.uniform(0,255) );
drawContours( drawing, contours, i, color1, 1, 8, vector<Vec4i>(), 0, Point() );
Scalar color2 = Scalar( rng.uniform(0, 255), rng.uniform(0,255), rng.uniform(0,255) );
drawContours( drawing, hull, i, color2, 1, 8, vector<Vec4i>(), 0, Point() );
int tt = contourArea(contours[i]);
printf("轮廓面积%d = %d\n", i, tt);
if( tt > area)
{
area = contourArea(contours[i]);
k = i;
}
}
vector<Point> hull1;
hull1 = hull[1];
for(i = 0; i< hull1.size(); i++ )
{
printf("point %d, %d, %d\n", i, hull1[i].x, hull1[i].y);
circle(drawing, hull1[i], 6, Scalar(255,0,0), 3, CV_AA);
}
int j;
for(j=0; j< contours.size(); j++)
{
//如果没有defects或者轮廓小于三个点，则continue
if( isContourConvex(contours[j])|| contours[j].size()<3) continue;
vector<int> convexHull_IntIdx;
vector<Vec4i> defects;
if (contours[j].size() >3 )
{
convexHull(contours[j], convexHull_IntIdx, true);
convexityDefects(contours[j], convexHull_IntIdx, defects);
}
for(i=0;i < defects.size();++i)
{
Matx<int,4,1> defectVector = defects[i];
vector<Point> contours1 =contours[j];
Point point1 = contours1[defectVector.val[0]];//开始点
Point point2 = contours1[defectVector.val[1]];//结束点
Point point3 = contours1[defectVector.val[2]];//深度点
float dist = defectVector.val[3];
printf("dist: %f \n", dist);
//if ( defectVector.val[3] <= 1000 ) { continue; } // skip defects that are shorter than 100 pixel
circle(drawing, point1, 3, Scalar(255,255,0), 2, CV_AA);
circle(drawing, point2, 8, Scalar(0,255,0), 2, CV_AA);
circle(drawing, point3, 3, Scalar(0,255,255), 2, CV_AA);
}
}
/// Show in a window
namedWindow( "Hull demo");
imshow( "Hull demo", drawing );
waitKey(0);
return(0);
}
程序执行之后界面如下，注意左下有图中
程序代码：工程FirstOpenCV25
