绘制多边形或者进行形状分析的时候,通常需要使用多边形逼近一个轮廓。
ApproxPoly
用指定精度逼近多边形曲线
CvSeq* cvApproxPoly( const void* src_seq, int header_size, CvMemStorage* storage, int method, double parameter, int parameter2=0 );
-
src_seq
- 点集数组序列 header_size
- 逼近曲线的头尺寸 storage
- 逼近轮廓的容器。如果为 NULL, 则使用输入的序列 method
- 逼近方法。目前仅支持 CV_POLY_APPROX_DP , 对应 Douglas-Peucker 算法. parameter
- 方法相关参数。对 CV_POLY_APPROX_DP 它是指定的逼近精度 parameter2
- 如果 src_seq 是序列,它表示要么逼近单个序列,要么在 src_seq 的同一个或低级层次上逼近所有序列 (参考 cvFindContours 中对轮廓继承结构的描述). 如果 src_seq 是点集的数组 (CvMat*) , 参数指定曲线是闭合 (parameter2!=0) 还是非闭合 (parameter2=0).
函数 cvApproxPoly 逼近一个或多个曲线,并返回逼近结果。对多个曲线的逼近,生成的树将与输入的具有同样的结构。(1:1 的对应关系).
用到的另一个函数的定义是:
#define cvContourPerimeter(contour)
这个函数作用于一个轮廓并返回其长度,事实上,此函数是一个调用通用函数cvArcLength()的宏
ArcLength
计算轮廓周长或曲线长度
double cvArcLength( const void* curve, CvSlice slice=CV_WHOLE_SEQ, int is_closed=-1 );
-
curve
- 曲线点集序列或数组 slice
- 曲线的起始点,缺省是计算整个曲线的长度 is_closed
-
表示曲线是否闭合,有三种情况:
- is_closed=0 - 假设曲线不闭合
- is_closed>0 - 假设曲线闭合
- is_closed<0 - 若曲线是序列,检查 ((CvSeq*)curve)->flags 中的标识 CV_SEQ_FLAG_CLOSED 来确定曲线是否闭合。否则 (曲线由点集的数组 (CvMat*) 表示) 假设曲线不闭合。
函数 cvArcLength 通过依次计算序列点之间的线段长度,并求和来得到曲线的长度。
函数代码:
/*
* 多边形逼近
*/
#include "highgui.h"
#include "cv.h"
void doApporx(IplImage* img)
{
IplImage* s =cvCreateImage(cvGetSize(img), IPL_DEPTH_8U, 1);
IplImage* dst = cvCreateImage(cvGetSize(img), img->depth, img->nChannels);
IplImage* dst2 = cvCreateImage(cvGetSize(img), img->depth, img->nChannels);
cvZero(dst);
cvZero(dst2);
CvMemStorage* storage = cvCreateMemStorage(0);
CvMemStorage* storage1 = cvCreateMemStorage(0);
CvSeq* contour = 0;
CvSeq* con ;
CvSeq* mcon = NULL;
cvCvtColor(img, s, CV_BGR2GRAY);
cvThreshold(s, s, 100, 200, CV_THRESH_BINARY);
cvFindContours(s,storage,&contour,sizeof(CvContour),CV_RETR_LIST, CV_CHAIN_APPROX_SIMPLE);
// 不使用多边形逼近的轮廓
cvDrawContours(dst2, contour, CV_RGB(0,0,255), CV_RGB(0,0,255), 2,2,8, cvPoint(0,0));
cvNamedWindow("cont",1);
cvShowImage("cont",dst2);
while(contour != 0)
{
// 使用对变形逼近
mcon = cvApproxPoly(contour, sizeof(CvContour), storage1, CV_POLY_APPROX_DP, cvContourPerimeter(contour)*0.02, 0);
cvDrawContours(dst, mcon, CV_RGB(0,0,255), CV_RGB(0,255,0), 2,2,8,cvPoint(0,0));
contour = contour->h_next;
}
cvNamedWindow("contour",1);
cvShowImage("contour",dst);
cvWaitKey(0);
}
显示结果:
原始图像:
结果显示:
参考博客:
http://buaagc.blog.163.com/blog/static/727883942013109145356/
http://blog.youkuaiyun.com/augusdi/article/details/9021939
图片来源: