OpenCV生成点集的Delaunay剖分和Voronoi图

实现内容：

设置一副图像大小为600*600。图像像素值全为0，为黑色。

在图像中Rect(100,100,400,400)的区域随机产生20个点。并画出。

产生这些点集的Delaunay剖分和Voronoi图。并画出。

程序

#include <opencv2/imgproc/imgproc.hpp>

#include <opencv2/highgui/highgui.hpp>

#include <iostream>

using namespace cv;

using namespace std;

//随机产生一个点集

vector<Point2f> generatePointSet(int n,Rect rect)

{

    vector<Point2f> pointSet;

    for (int i = 0; i < n;i++)

    {

        Point2f fp((float)(rand() % (rect.width - 2*rect.x) + rect.x),

            (float)(rand() % (rect.height - 2*rect.y) + rect.y));

        pointSet.push_back(fp);

    }

    return pointSet;

}

//标记出点

static void drawPoint(Mat& img, Point2f fp, Scalar color)

{

    circle(img, fp, 3, color, CV_FILLED, 8, 0);

}

//画出点集

static void drawPointSet(Mat& img, vector<Point2f> pointSet, Scalar color)

{

    for (int i = 0; i < pointSet.size();i++)

    {

        drawPoint(img, pointSet[i], color);

    }

}

//画出剖分

static void drawSubdiv(Mat& img, Subdiv2D& subdiv, Scalar delaunay_color)

{

    vector<Vec6f> triangleList;

    subdiv.getTriangleList(triangleList);

    vector<Point> pt(3);

    for (size_t i = 0; i < triangleList.size(); i++)

    {

        Vec6f t = triangleList[i];

        pt[0] = Point(cvRound(t[0]), cvRound(t[1]));

        pt[1] = Point(cvRound(t[2]), cvRound(t[3]));

        pt[2] = Point(cvRound(t[4]), cvRound(t[5]));

        line(img, pt[0], pt[1], delaunay_color, 1, CV_AA, 0);

        line(img, pt[1], pt[2], delaunay_color, 1, CV_AA, 0);

        line(img, pt[2], pt[0], delaunay_color, 1, CV_AA, 0);

    }

}

//画出Voronoi图

static void paintVoronoi(Mat& img, Subdiv2D& subdiv)

{

    vector<vector<Point2f> > facets;

    vector<Point2f> centers;

    subdiv.getVoronoiFacetList(vector<int>(), facets, centers);

    vector<Point> ifacet;

    vector<vector<Point> > ifacets(1);

    for (size_t i = 0; i < facets.size(); i++)

    {

        ifacet.resize(facets[i].size());

        for (size_t j = 0; j < facets[i].size(); j++)

            ifacet[j] = facets[i][j];

        Scalar color;

        color[0] = rand() & 255;

        color[1] = rand() & 255;

        color[2] = rand() & 255;

        fillConvexPoly(img, ifacet, color, 8, 0);

        ifacets[0] = ifacet;

        polylines(img, ifacets, true, Scalar(), 1, CV_AA, 0);

        circle(img, centers[i], 3, Scalar(), CV_FILLED, CV_AA, 0);

    }

}

int main()

{

    //创建点集

    Rect rect(100,100,400,400);

    Mat img(rect.size(), CV_8UC3);

    img = Scalar::all(0);

    vector<Point2f> pointset = generatePointSet(20,rect);

    Mat img_ptst = img.clone();

    drawPointSet(img_ptst,pointset,Scalar(0,255,0));

    imshow("Point set",img_ptst);

    imwrite("pointSet.jpg",img_ptst);

    //创建Delaunay剖分

    Subdiv2D subdiv(rect);

    for (int i = 0; i < pointset.size();i++)

    {

        subdiv.insert(pointset[i]);

    }

    //画出Delaunay剖分三角形

    Mat img_delaunay = img.clone();

    drawSubdiv(img_delaunay, subdiv, Scalar(255,255,255));

    imshow("Delaunay", img_delaunay);

    imwrite("delaunay.jpg", img_delaunay);

    //画出Voronoi图

    Mat img_voronoi = img.clone();

    paintVoronoi(img_voronoi, subdiv);

    imshow("Voronoi", img_voronoi);

    imwrite("voronoi.jpg", img_voronoi);

    waitKey(0);

    return 0;

}

结果