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

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

下面是一个有点集序列或数组创建凸边形的简单示例

[cpp]  view plain copy print ?
  1. #include<cv.h>  
  2. #include<highgui.h>  
  3. #include<stdlib.h>  
  4.   
  5. #define ARRAY 0  
  6.   
  7. int main()  
  8. {  
  9.     IplImage* img = cvCreateImage (cvSize(500, 500), 8, 3);  
  10.     cvNamedWindow ("hull", 1);  
  11.   
  12. #if ! ARRAY  
  13.     CvMemStorage* storage = cvCreateMemStorage (0);  
  14. #endif  
  15.     int i;  
  16.     int contour = rand() % 100 + 1;  
  17.     int hullcontour;  
  18.     CvPoint pt0;  
  19. #if ! ARRAY  
  20.     CvSeq* ptseq = cvCreateSeq (CV_SEQ_KIND_GENERIC | CV_32SC2, sizeof(CvContour), sizeof(CvPoint), storage);  
  21.     CvSeq* hull;  
  22.     for (i = 0; i < contour; i++)  
  23.     {  
  24.         pt0.x = rand() % (img->width / 2) + img->width / 4;  
  25.         pt0.y = rand() % (img->height / 2) + img->height / 4;  
  26.         cvSeqPush (ptseq, &pt0);  
  27.     }  
  28.     hull = cvConvexHull2 (ptseq, 0, CV_CLOCKWISE, 0);  
  29.     /*hull = cvConvexHull2 (ptseq, 0, CV_CLOCKWISE, 1);*/  
  30.     hullcontour = hull->total;  
  31. #else  
  32.     CvPoint* points = (CvPoint*)malloc(contour * sizeof(points[0]));  
  33.     int* hull = (int*)malloc(contour * sizeof(hull[0]));  
  34.     CvMat point_mat = cvMat (1, contour, CV_32SC2, points);  
  35.     CvMat hull_mat = cvMat (1, contour, CV_32SC1, hull);  
  36.     for (i = 0; i < contour; i++)  
  37.     {  
  38.         pt0.x = rand() % (img->width / 2) + img->width / 4;  
  39.         pt0.y = rand() % (img->height / 2) + img->height / 4;  
  40.         points[i] = pt0;  
  41.     }  
  42.     cvConvexHull2 (&point_mat, &hull_mat, CV_CLOCKWISE, 0);  
  43.     hullcontour = hull_mat.cols;  
  44. #endif  
  45.     cvZero (img);  
  46.     for (i = 0; i < contour; i++)  
  47.     {  
  48. #if ! ARRAY  
  49.         pt0 = * CV_GET_SEQ_ELEM(CvPoint, ptseq, i);  
  50.         /*pt0 = ** CV_GET_SEQ_ELEM(CvPoint*, ptseq, i);*/  
  51. #else  
  52.         pt0 = points[i];  
  53. #endif  
  54.         cvCircle (img, pt0, 2, CV_RGB(255, 0, 0), CV_FILLED);  
  55.     }  
  56. #if ! ARRAY  
  57.     pt0 = ** CV_GET_SEQ_ELEM(CvPoint*, hull, hullcontour - 1);  
  58.     /*pt0 = * CV_GET_SEQ_ELEM(CvPoint, hull, hullcontour - 1);*/  
  59. #else  
  60.     pt0 = points[hull[hullcontour - 1]];  
  61. #endif  
  62.   
  63.     for (i = 0; i < hullcontour; i++)  
  64.     {  
  65. #if ! ARRAY  
  66.         CvPoint pt = **CV_GET_SEQ_ELEM(CvPoint*, hull, i);  
  67.         /*CvPoint pt = *CV_GET_SEQ_ELEM(CvPoint, hull, i);*/  
  68. #else  
  69.         CvPoint pt = points[hull[i]];  
  70. #endif  
  71.         cvLine (img, pt0, pt, CV_RGB(255, 0, 0));  
  72.         pt0 = pt;  
  73.     }  
  74.   
  75.     cvShowImage ("hull", img);  
  76.   
  77.     cvWaitKey (0);  
  78.   
  79. #if ! ARRAY   
  80.     cvClearMemStorage (storage);  
  81. #else  
  82.     free (points);  
  83.     free (hull);  
  84. #endif  
  85.     return 0;  
  86.   
  87. }  


代码中有3条注释语句,这是我在学习这个程序时遇到的问题,自己写的测试代码。这个代码相对来说还是比较

好理解的,主要就是开始对程序中CV_GET_SEQ_ELEM宏的使用有些不太理解。这个宏是用来在序列中提取元

素的,本程序中共有3处用到了。

第一处:

[cpp]  view plain copy print ?
  1. pt0 = * CV_GET_SEQ_ELEM(CvPoint, ptseq, i);  

第二处:

[cpp]  view plain copy print ?
  1. pt0 = ** CV_GET_SEQ_ELEM(CvPoint*, hull, hullcontour - 1);  

第三处:

[cpp]  view plain copy print ?
  1. CvPoint pt = **CV_GET_SEQ_ELEM(CvPoint*, hull, i);  

显然,第一处的调用方式和后面两处是不同的。但这里pt0pt的类型都是CvPoint,ptseqhull都是CvSeq*

型的,唯一不同的是宏里的第一个参数,这个参数表上宏返回的类型+*,比如:第一处这个参数是CvPoint,则

返回CvPoint*,以此类推。这样的话这两个宏最后结果都是返回一个CvPoint类型值。这里对宏的使用是不是有

点类型于函数重载呢!既然如此,我想干脆都用相同的方式调用这个宏,于是我把第一处的代码改为

[cpp]  view plain copy print ?
  1. pt0 = ** CV_GET_SEQ_ELEM(CvPoint*, ptseq, i);  

运行后程序中止了,于是我把这里又改回来,然后改动后两处,改为

[cpp]  view plain copy print ?
  1. pt0 = * CV_GET_SEQ_ELEM(CvPoint, hull, hullcontour - 1);  
  2.   
  3. CvPoint pt = *CV_GET_SEQ_ELEM(CvPoint, hull, i);  

运行后程序依然中止啦。这样起码说明了在这里关于这个宏的调用方式是不可以互换的。但具体为什么,确实

纠结了很久,当就要绝望时,在Opencv中文论坛上看到了关于使用cvConvexHull2函数得到的凸包的问题,里

面有句话对我很有帮助:当return_points=0时,用cvConvexHull2函数得到的是凸外形,包含的是轮廓的定点

的指针或下标,而当return_points非0时,得到的是外形点本身。cvConvexHull2函数调用代码如下

[cpp]  view plain copy print ?
  1. hull = cvConvexHull2 (ptseq, 0, CV_CLOCKWISE, 0);  

这里的return_points被设置为0了,说明hull并不是直接指向轮廓点的序列,而是指向指向轮廓点序列的指针,

这就有点像hull是个二级指针啦!所以在后两处对宏的调用第一个参数设置为CvPoint*,这样的话就返回

CvPoint**类型。到目前为止,这还只是我个人推测,还需要进一步的验证。于是我把上面代码中最后一个参数

return_points改为1,运行。一切正常啦!注:这里后面两处对宏的调用已经在前面改了,所以就不必在修改

了。以上所有的分析和理解都只是我个人的看法,不一定就是对的。希望大家能够一起交流下这个问题,也希

望各位大神能够不吝赐教!!!

你可能感兴趣的:(由点集序列或数组创建凸边形)