OpenCV3学习（9.3）图像轮廓的多边形逼近approxPolyDP、轮廓包围框\圆\椭圆、凸包

1. 多边形逼近approxPolyDP

    当需要对图像进行形状分析时，需要使用多边形逼近一个轮廓，使得顶点数目变少，算法原理比较简单，核心就是不断找多边形最远的点加入形成新的多边形，直到最短距离小于指定的精度。OpenCV里面用函数approxPolyDP()实现。approxPolyDP()用另一条顶点较少的曲线来逼近一条曲线或者一个多边形，这样两条曲线之间的距离小于或等于指定的精度。同时也有使闭合逼近曲线的选项（那就是说，起始点和终止点相同）。

     findContours后的轮廓信息contours可能过于复杂不平滑，可以用approxPolyDP函数对该多边形曲线做适当近似，

approxPolyDP 主要功能是把一个连续光滑曲线折线化，对图像轮廓点进行多边形拟合。原理图：对比之前黑点连线，之后蓝色连线：

                                       

函数原型：

void approxPolyDP( InputArray curve,
                                OutputArray approxCurve,
                                double epsilon, bool closed );

//调用方式
vector<vector<Point>> contours_poly(contours.size());//用于存放折线点集
approxPolyDP(Mat(contours[i]), contours_poly[i], 15, true);//contours[i]是有findcontours函数得到，
drawContours(dstImg, contours_poly, i, Scalar(0, 255, 255), 2, 8);  //绘制

参数详解;

InputArray curve:一般是由图像的轮廓点组成的点集

OutputArray approxCurve：表示输出的多边形点集

double epsilon：主要表示输出的精度，就是各个轮廓点之间最大距离数，这个参数表示的是精度，越小精度越高，因为表示的意思是是原始曲线与近似曲线之间的最大距离

bool closed：表示输出的多边形是否封闭，若设为true，则认为曲线是闭合的，第一个点和最后一个点相连。

理解Douglas-Peucker算法：

1、该算法从轮廓中挑出两个最远的点，进行相连；

2、然后再原轮廓上寻找一个离线段距离最远的点，将该点加入逼近后的新轮廓，即连接着三个点形成的三角型作为轮廓；

3、选择三角形的任意一条边出发，进行步骤二，将距离最远点加入新轮廓，直至满足输出的精度要求。

 

from： https://www.cnblogs.com/k1412/p/6884484.html

2、计算轮廓周长或曲线长度（不适用于非轮廓的点集）

double arcLength( InputArray curve, bool closed );

curve 曲线点集序列或数组 
closed 表示曲线是否闭合： closed=true- 假设曲线闭合 ，则最后一点到第一点的距离计入总弧长。函数 arcLength 通过依次计算序列点之间的线段长度，并求和来得到曲线的长度

     当我们得到对象轮廓后，可用boundingRect()得到包覆此轮廓的最小正矩形，minAreaRect()得到包覆轮廓的最小斜矩形，minEnclosingCircle()得到包覆此轮廓的最小圆形，这些函式协助我们填补空隙，或者作进一步的对象辨识，boundingRect()函式返回的是正矩形，所以如果对象有倾斜的情形，返回的可能不是我们想要的结果。

3、获得矩形包围框（适合所有点集）

Rect boundingRect( InputArray points );

points：输入信息，可以为包含点的容器(vector)或是Mat。

返回包覆输入信息的最小正矩形

 4、获取最小矩形框

 通过建立凸外形并且旋转外形以寻找给定 2D 点集的最小面积的包围矩形。这个长方形很有可能是倾斜的。

RotatedRect minAreaRect( InputArray points );

points：输入信息，可以为包含点的容器(vector)或是Mat。

返回包覆输入信息的最小斜矩形。

矩形边框，用一个最小的矩形，把找到的形状包起来。还有一个带旋转的矩形，面积会更小，如下图：

                                                   

RotatedRect该类表示平面上的旋转矩形，有三个属性：

矩形中心点(质心)

边长(长和宽)

旋转角度

//定义
class CV_EXPORTS RotatedRect
{   
public:  
    RotatedRect(); //构造函数
    RotatedRect(const Point2f& center, const Size2f& size, float angle);
    RotatedRect(const CvBox2D& box);
    void points(Point2f pts[]) const;   //返回矩形的4个顶点
    Rect boundingRect() const;  //返回包含旋转矩形的最小矩形
    operator CvBox2D() const;   //转换到旧式的cvbox2d结构
    Point2f center;  //矩形的质心
    Size2f size;     //矩形的边长
    float angle;     //旋转角度，当角度为0、90、180、270等时，矩形就成了一个直立的矩形
};

//使用范例：
    Mat test_image(200, 200, CV_8UC3, Scalar(0));
    RotatedRect rRect = RotatedRect(Point2f(100,100), Size2f(100,50), 30);
    Point2f vertices[4];
    rRect.points(vertices);
    for (int i = 0; i < 4; i++)
        line(test_image, vertices[i], vertices[(i+1)%4], Scalar(0,255,0), 2);
    Rect brect = rRect.boundingRect();
    rectangle(test_image, brect, Scalar(255,0,0), 2);
    imshow("rectangles", test_image);
    waitKey(0);

5、获取最小包围圆

void minEnclosingCircle( InputArray points,
                                      CV_OUT Point2f& center, CV_OUT float& radius );

 

points：输入信息，可以为包含点的容器(vector)或是Mat。

center：包覆圆形的圆心。

radius：包覆圆形的半径。

6、获得椭圆边框

对给定的一组二维点集作椭圆的最佳拟合(最小二乘意义上的)。返回结构 RotatedRect， size 表示椭圆轴的长轴和短轴。

RotatedRect fitEllipse( InputArray points );

points ：点集的序列或数组 

7、获得轮廓的最佳拟合直线

	void fitLine(
		cv::InputArray points, // 二维点的数组或vector
		cv::OutputArray line, // 输出直线,Vec4f (2d)或Vec6f (3d)的vector
		int distType, // 距离类型
		double param, // 距离参数
		double reps, // 径向的精度参数
		double aeps // 角度精度参数
	);

第一个参数是用于拟合直线的输入点集，可以是二维点的cv::Mat数组，也可以是二维点的STL vector。

第二个参数是输出的直线，对于二维直线而言类型为cv::Vec4f，对于三维直线类型则是cv::Vec6f，输出参数的前半部分给出的是直线的方向，而后半部分给出的是直线上的一点（即通常所说的点斜式直线）。

第三个参数是距离类型，拟合直线时，要使输入点到拟合直线的距离和最小化（即下面公式中的cost最小化），可供选的距离类型如下表所示，ri表示的是输入的点到直线的距离。
           

第四个参数是距离参数，跟所选的距离类型有关，值可以设置为0，cv::fitLine()函数本身会自动选择最优化的值

第五、六两个参数用于表示拟合直线所需要的径向和角度精度，通常情况下两个值均被设定为1e-2

from:https://blog.csdn.net/guduruyu/article/details/69505487

 8、计算整个轮廓或部分轮廓的面积

double contourArea( InputArray contour, bool oriented=false );

contour:findcontour函数得到的轮廓；

oriented: 轮廓的方向影响面积的符号。因此函数也许会返回负的结果。应用函数 fabs() 得到面积的绝对值

9、获得轮廓的凸包

指一个完全包含原轮廓且只由原轮廓上的点组成的多边形。

 void convexHull( InputArray points, OutputArray hull,
                              bool clockwise=false, bool returnPoints=true );

points:输入的二维点集，Mat类型数据即可 
hull:输出参数，用于输出函数调用后找到的凸包 
clockwise:操作方向，当标识符为真时，输出凸包为顺时针方向，否则为逆时针方向。 
returnPoints:操作标识符，默认值为true，此时返回各凸包的各个点，否则返回凸包各点的指数，当输出数组时std::vector时，此标识被忽略。

10、测试点是否在多边形中

double pointPolygonTest( InputArray contour, Point2f pt, bool measureDist );

contour：输入轮廓.；                        pt针对轮廓需要测试的点。

measure_dist：如果非0，函数将估算点到轮廓最近边的距离。

       函数PointPolygonTest 决定测试点是否在轮廓内，轮廓外，还是轮廓的边上（或者共边的交点上），当measure_dist为false时，返回值对应是1, -1,0；当 measure_dist为true ，它返回一个从点到最近的边的带符号距离。

11、测试轮廓是否为凸

bool isContourConvex( InputArray contour );