由點集序列或數組創建凸邊形

[cpp] 
#include<cv.h> 
#include<highgui.h> 
#include<stdlib.h> 
 
#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處用到了。
第一處：
[cpp] 
pt0 = * CV_GET_SEQ_ELEM(CvPoint, ptseq, i); 

第二處：
[cpp] ll, hullcontour - 1); 

第三處：
[cpp] 
CvPoint pt = **CV_GET_SEQ_ELEM(CvPoint*, hull, i); 

顯然，第一處的調用方式和後面兩處是不同的。但這裡pt0和pt的類型都是CvPoint，ptseq和hull都是CvSeq*類
型的，唯一不同的是宏裡的第一個參數，這個參數表上宏返回的類型+*，比如：第一處這個參數是CvPoint,則
返回CvPoint*，以此類推。這樣的話這兩個宏最後結果都是返回一個CvPoint類型值。這裡對宏的使用是不是有
點類型於函數重載呢！既然如此，我想干脆都用相同的方式調用這個宏，於是我把第一處的代碼改為
[cpp] 
pt0 = ** CV_GET_SEQ_ELEM(CvPoint*, ptseq, i); 

運行後程序中止了，於是我把這裡又改回來，然後改動後兩處，改為
[cpp]
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函數調用代碼如下
[cpp] 
hull = cvConvexHull2 (ptseq, 0, CV_CLOCKWISE, 0); 

這裡的return_points被設置為0了，說明hull並不是直接指向輪廓點的序列，而是指向指向輪廓點序列的指針，
這就有點像hull是個二級指針啦!所以在後兩處對宏的調用第一個參數設置為CvPoint*，這樣的話就返回
CvPoint**類型。到目前為止，這還只是我個人推測，還需要進一步的驗證。於是我把上面代碼中最後一個參數
return_points改為1，運行。一切正常啦！注：這裡後面兩處對宏的調用已經在前面改了，所以就不必在修改
了。以上所有的分析和理解都只是我個人的看法，不一定就是對的。