椭圆拟合
看了網(wǎng)上光宇橢圓擬合的各種算法,都差不多,但是寫的條理很亂,借鑒別人的成果,修改了下程序,下邊是我梳理后的過程。
#include "cv.h"
#include "highgui.h"
int main()
{ ? //定義圖像
IplImage *img3 = cvLoadImage("spot.JPG",0);
IplImage *img2 = cvCloneImage(img3);
//將原圖img3閾值化為img2;
cvThreshold(img3, img2, 50, 255, CV_THRESH_BINARY);
IplImage *img4 = cvCloneImage(img3);
cvNamedWindow("Source", 1);
cvNamedWindow("Result", 1);
cvShowImage("Source", img3);
//設(shè)置橢圓擬合用到的幾個參數(shù)
CvMemStorage *stor = cvCreateMemStorage(0);
CvSeq *cont = cvCreateSeq(CV_SEQ_ELTYPE_POINT, sizeof(CvSeq), sizeof(CvPoint), stor);
CvBox2D32f *box;
CvPoint *PointArray;
CvPoint2D32f *PointArray2D32f;
? ? //檢索橢圓的輪廓,存儲在stor
cvFindContours(img2, stor, &cont, sizeof(CvContour),
CV_RETR_LIST, CV_CHAIN_APPROX_NONE, cvPoint(0, 0));
//清空圖像img2,img4;
cvZero(img2);
cvZero(img4);
//繪制所有輪廓并用橢圓擬合
for (; cont; cont = cont->h_next)
{
int i;
int count = cont->total;//輪廓個數(shù)
CvPoint center;
CvSize size;
/*個數(shù)必須大于6,這是cvFitEllipse_32f的要求*/
if (count<6)
{
continue;
}
//分配內(nèi)存給點集
PointArray = (CvPoint *)malloc(count*sizeof(CvPoint));
PointArray2D32f = (CvPoint2D32f*)malloc(count*sizeof(CvPoint2D32f));
//分配內(nèi)存給橢圓數(shù)據(jù)
box = (CvBox2D32f *)malloc(sizeof(CvBox2D32f));
//得到點集(這個方法值得借鑒)
cvCvtSeqToArray(cont, PointArray, CV_WHOLE_SEQ);
//將CvPoint點集轉(zhuǎn)化為CvBox2D32f集合
for (i = 0; i<count; i++)
{
PointArray2D32f[i].x = (float)PointArray[i].x;
PointArray2D32f[i].y = (float)PointArray[i].y;
}
//擬合當(dāng)前輪廓
cvFitEllipse(PointArray2D32f, count, box);
//繪制當(dāng)前輪廓
cvDrawContours(img4, cont, CV_RGB(255, 255, 255), CV_RGB(255, 255, 255),
0, 1, 8, cvPoint(0, 0));
//將橢圓數(shù)據(jù)從浮點轉(zhuǎn)化為整數(shù)表示
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;
//畫橢圓
cvEllipse(img4, center, size, box->angle, 0, 360, CV_RGB(0, 0, 255), 1, CV_AA, 0);
free(PointArray);
free(PointArray2D32f);
free(box);
}
cvShowImage("Result", img4);
cvWaitKey(0);
cvSaveImage("1.JPG", img4);
cvReleaseImage(&img2);
cvReleaseImage(&img3);
cvDestroyWindow("Source");
cvDestroyWindow("Result");
return 0;
}
總結(jié)
- 上一篇: tp5 php 使用array_colu
- 下一篇: 苹果开发者中心上传APP屏幕快照