由点集序列或数组创建凸边形

2019-04-13 15:52发布

站内文章 / 模拟电子
9446 0 1195
下面是一个有点集序列或数组创建凸边形的简单示例 #include #include #include #define ARRAY 0 int main() { IplImage* img = cvCreateImage (cvSize(500, 500), 8, 3); cvNamedWindow ("hull", 1); #if ! ARRAY CvMemStorage* storage = cvCreateMemStorage (0); #endif int i; int contour = rand() % 100 + 1; int hullcontour; CvPoint pt0; #if ! ARRAY CvSeq* ptseq = cvCreateSeq (CV_SEQ_KIND_GENERIC | CV_32SC2, sizeof(CvContour), sizeof(CvPoint), storage); CvSeq* hull; for (i = 0; i < contour; i++) { pt0.x = rand() % (img->width / 2) + img->width / 4; pt0.y = rand() % (img->height / 2) + img->height / 4; cvSeqPush (ptseq, &pt0); } hull = cvConvexHull2 (ptseq, 0, CV_CLOCKWISE, 0); /*hull = cvConvexHull2 (ptseq, 0, CV_CLOCKWISE, 1);*/ hullcontour = hull->total; #else CvPoint* points = (CvPoint*)malloc(contour * sizeof(points[0])); int* hull = (int*)malloc(contour * sizeof(hull[0])); CvMat point_mat = cvMat (1, contour, CV_32SC2, points); CvMat hull_mat = cvMat (1, contour, CV_32SC1, hull); for (i = 0; i < contour; i++) { pt0.x = rand() % (img->width / 2) + img->width / 4; pt0.y = rand() % (img->height / 2) + img->height / 4; points[i] = pt0; } cvConvexHull2 (&point_mat, &hull_mat, CV_CLOCKWISE, 0); hullcontour = hull_mat.cols; #endif cvZero (img); for (i = 0; i < contour; i++) { #if ! ARRAY pt0 = * CV_GET_SEQ_ELEM(CvPoint, ptseq, i); /*pt0 = ** CV_GET_SEQ_ELEM(CvPoint*, ptseq, i);*/ #else pt0 = points[i]; #endif cvCircle (img, pt0, 2, CV_RGB(255, 0, 0), CV_FILLED); } #if ! ARRAY pt0 = ** CV_GET_SEQ_ELEM(CvPoint*, hull, hullcontour - 1); /*pt0 = * CV_GET_SEQ_ELEM(CvPoint, hull, hullcontour - 1);*/ #else pt0 = points[hull[hullcontour - 1]]; #endif for (i = 0; i < hullcontour; i++) { #if ! ARRAY CvPoint pt = **CV_GET_SEQ_ELEM(CvPoint*, hull, i); /*CvPoint pt = *CV_GET_SEQ_ELEM(CvPoint, hull, i);*/ #else CvPoint pt = points[hull[i]]; #endif cvLine (img, pt0, pt, CV_RGB(255, 0, 0)); pt0 = pt; } cvShowImage ("hull", img); cvWaitKey (0); #if ! ARRAY cvClearMemStorage (storage); #else free (points); free (hull); #endif return 0; }

代码中有3条注释语句，这是我在学习这个程序时遇到的问题，自己写的测试代码。这个代码相对来说还是比较 好理解的，主要就是开始对程序中CV_GET_SEQ_ELEM宏的使用有些不太理解。这个宏是用来在序列中提取元 素的，本程序中共有3处用到了。 第一处： pt0 = * CV_GET_SEQ_ELEM(CvPoint, ptseq, i);
第二处： pt0 = ** CV_GET_SEQ_ELEM(CvPoint*, hull, hullcontour - 1);
第三处： CvPoint pt = **CV_GET_SEQ_ELEM(CvPoint*, hull, i);
显然，第一处的调用方式和后面两处是不同的。但这里pt0和pt的类型都是CvPoint，ptseq和hull都是CvSeq*类 型的，唯一不同的是宏里的第一个参数，这个参数表上宏返回的类型+*，比如：第一处这个参数是CvPoint,则 返回CvPoint*，以此类推。这样的话这两个宏最后结果都是返回一个CvPoint类型值。这里对宏的使用是不是有 点类型于函数重载呢！既然如此，我想干脆都用相同的方式调用这个宏，于是我把第一处的代码改为 pt0 = ** CV_GET_SEQ_ELEM(CvPoint*, ptseq, i);
运行后程序中止了，于是我把这里又改回来，然后改动后两处，改为 pt0 = * CV_GET_SEQ_ELEM(CvPoint, hull, hullcontour - 1); CvPoint pt = *CV_GET_SEQ_ELEM(CvPoint, hull, i);
运行后程序依然中止啦。这样起码说明了在这里关于这个宏的调用方式是不可以互换的。但具体为什么，确实 纠结了很久，当就要绝望时，在Opencv中文论坛上看到了关于使用cvConvexHull2函数得到的凸包的问题，里 面有句话对我很有帮助：当return_points=0时，用cvConvexHull2函数得到的是凸外形，包含的是轮廓的定点 的指针或下标，而当return_points非0时，得到的是外形点本身。cvConvexHull2函数调用代码如下 hull = cvConvexHull2 (ptseq, 0, CV_CLOCKWISE, 0);
这里的return_points被设置为0了，说明hull并不是直接指向轮廓点的序列，而是指向指向轮廓点序列的指针， 这就有点像hull是个二级指针啦!所以在后两处对宏的调用第一个参数设置为CvPoint*，这样的话就返回 CvPoint**类型。到目前为止，这还只是我个人推测，还需要进一步的验证。于是我把上面代码中最后一个参数 return_points改为1，运行。一切正常啦！注：这里后面两处对宏的调用已经在前面改了，所以就不必在修改 了。以上所有的分析和理解都只是我个人的看法，不一定就是对的。希望大家能够一起交流下这个问题，也希 望各位大神能够不吝赐教！！！