在opencv中,通過函數convexHulll能很容易的得到一系列點的凸包,比如由點組成的輪廓,通過convexHull函數,我們就能得到輪廓的凸包。下面的圖就是一些點集的凸包。
求凸包的代碼如下:
int main( int /*argc*/, char** /*argv*/ )
{
Mat img(500, 500, CV_8UC3);
RNG& rng = theRNG();
cout << "\n這個程序演示了凸包函數的使用,任意給定一些點,求出包圍這些點的凸包\n" <<endl;
for(;;)
{
char key;
int i, count = (unsigned)rng%100 + 1;
vector<Point> points;
//隨機在1-100個點,這些點位於圖像中心3/4處。
for( i = 0; i < count; i++ )
{
Point pt;
pt.x = rng.uniform(img.cols/4, img.cols*3/4);
pt.y = rng.uniform(img.rows/4, img.rows*3/4);
points.push_back(pt);
}
//計算凸包
vector<int> hull;
convexHull(Mat(points), hull, true);
//畫隨即點
img = Scalar::all(0);
for( i = 0; i < count; i++ )
circle(img, points[i], 3, Scalar(0, 0, 255), CV_FILLED, CV_AA);
int hullcount = (int)hull.size();
Point pt0 = points[hull[hullcount-1]];
//畫凸包
for( i = 0; i < hullcount; i++ )
{
Point pt = points[hull[i]];
line(img, pt0, pt, Scalar(0, 255, 0), 1, CV_AA);
pt0 = pt;
}
imshow("hull", img);
key = (char)waitKey();
if( key == 27 || key == 'q' || key == 'Q' ) // 'ESC'
break;
}
return 0;
}
convexHull(Mat(points), hull, true);
convex
convexHull第一個參數是要求凸包的點集,第二個參數是輸出的凸包點,第三個參數是一個bool變量,表示求得的凸包是順時針方向還是逆時針方向,true是順時針方向。注意:第二個參數可以為vector<int>,此時返回的是凸包點在原輪廓點集中的索引,也可以為vector<Point>,此時存放的是凸包點的位置。
OpenCV中求凸包的算法參見paper:
Sklansky, J., Finding the Convex Hull of a Simple Polygon. PRL 1 $number, pp 79-83 (1982)
程序代碼:參考工程FirstOpenCV24