Opencv椭圆拟合

一 轮廓检测

在进行椭圆拟合时需要先检测出轮廓中的点,这就需要用到cvFindContous命令。

cvFindContours(CvArr* image, CvSeq** first_contour,int header_size CV_DEFAULT(sizeof(CvContour)),int mode CV_DEFAULT(CV_RETR_LIST),int method CV_DEFAULT(CV_CHAIN_APPROX_SIMPLE),CvPoint offset CV_DEFAULT(cvPoint(0,0)))

image: 是一个二值图像,一般由cvThreShold得到
storage: 返回轮廓的容器
first_contour: 输出函数, 指向第一个轮廓的地址
header_size: 序列的尺寸,如果method=CV_CHAIN_CODE, 则
header_size=sizeof(CvChain),否则为sizeof(CvContour)
mode: CV_RETR_EXTERNAL, 只检索最外轮廓
CV_RETR_LIST, 采用横向列表方式检索所用的轮廓
method: CV_CHAIN_CODE - Freeman 链码的输出轮廓. 其它方法输出多边形(定点序列).
CV_CHAIN_APPROX_NONE - 将所有点由链码形式翻译(转化）为点序列形式
CV_CHAIN_APPROX_SIMPLE - 压缩水平、垂直和对角分割，即函数只保留末端的象素点
CV_CHAIN_APPROX_TC89_L1
CV_CHAIN_APPROX_TC89_KCOS - 应用 Teh-Chin 链逼近算法. CV_LINK_RUNS - 通过连接为 1 的水平碎片使用完全不同的轮廓提取算法。仅有 CV_RETR_LIST 提取模式可以在本方法中应用
offset: 偏移量，用于移动所有轮廓点

二 椭圆拟合

cvFitEllipse( const CvPoint2D32f* points, int count, CvBox2D* box )

三 程序操作


CvMemStorage*stor;
CvSeq*cont;
CvFont font;
CvBox2D32f*box;
CvPoint*PointArray;
CvPoint2D32f*PointArray2D32f;
double minarea=10,maxarea=500;
stor=cvCreateMemStorage(0);
cont=cvCreateSeq(CV_SEQ_ELTYPE_POINT,sizeof(CvSeq),sizeof(CvPoint),stor);
cvThreshold(image03,image02,slider_pos,255,CV_THRESH_BINARY);
cvFindContours(image02,stor,&cont,sizeof(CvContour),
CV_RETR_LIST,CV_CHAIN_APPROX_NONE,cvPoint(0,0));
for(;cont;cont=cont->h_next)
{
int i;
int count=cont->total;//轮廓个数
CvPoint center;
CvSize size;
/*个数必须大于6，这是cvFitEllipse_32f的要求*/
if(count<6)
{
continue;
}
double a=abs(cvContourArea(cont));
cvInitFont(&font,CV_FONT_HERSHEY_COMPLEX,0.5f,0.5f,0,1,8); //初始化字体
if (a>=minarea&&a<=maxarea)
{
//分配内存给点集
PointArray=(CvPoint*)malloc(count*sizeof(CvPoint));
PointArray2D32f= (CvPoint2D32f*)malloc(count*sizeof(CvPoint2D32f));
//分配内存给椭圆数据
box=(CvBox2D32f*)malloc(sizeof(CvBox2D32f));
//得到点集
cvCvtSeqToArray(cont,PointArray,CV_WHOLE_SEQ);
//将CvPoint点集转化为CvBox2D32f集合
for(i=0;i<count;i++)
{
PointArray2D32f[i].x=(float)PointArray[i].x;
PointArray2D32f[i].y=(float)PointArray[i].y;
}
//拟合当前轮廓
cvFitEllipse(PointArray2D32f,count,box);
//绘制当前轮廓 cvDrawContours(image04,cont,CV_RGB(255,255,255),CV_RGB(255,255,255),0,1,8,cvPoint(0,0));
//将椭圆数据从浮点转化为整数表示
center.x=cvRound(box->center.x);
center.y=cvRound(box->center.y);
size.width=cvRound(box->size.width*0.5);
size.height=cvRound(box->size.height*0.5);
box->angle=-box->angle;
centerfile<<box->center.x<<" "<<box->center.y<<endl;
//画椭圆
cvEllipse(image04,center,size,box- >angle,0,360,CV_RGB(0,0,255),1,CV_AA,0);
free(PointArray);
free(PointArray2D32f);
free(box);
}
}