opencv 提取曲线并用矩形框框出来
老师提出写个提取裂缝,并用矩形框框出来的小程序,不用多好,能跑起来看看就行,由于裂纹没有现成的,就随便划了 黑线做个模仿,由于水平比较菜,就在网上抄些程序改改,拼拼,好歹算搞定了,也算个demo吧,嘿嘿~~~希望对大家有帮助
原图在这:,,,额 ,真的好大,算了 ,不弄了,就这样了
下面直接贴代码喽~~~
#include "stdafx.h"
#include "stdio.h"
#include "cv.h"
#include "highgui.h"
#include
#include
using namespace std;
using namespace cv;
void sum_rgb( IplImage* src, IplImage* dst) //该函数的作用就是转灰度后二值化
{
IplImage* r=cvCreateImage(cvGetSize(src), IPL_DEPTH_8U,1);
IplImage* g=cvCreateImage(cvGetSize(src), IPL_DEPTH_8U,1);
IplImage* b=cvCreateImage(cvGetSize(src), IPL_DEPTH_8U,1);//函数 cvSplit 分割多通道数组成分离的单通道数组d
cvSplit(src,r,g,b,NULL);
IplImage* s=cvCreateImage(cvGetSize(src), IPL_DEPTH_8U,1);cvAddWeighted(r, 1./3, g, 1./3, 0.0, s);
cvAddWeighted(s, 2./3, b, 1./3, 0.0, s);
cvThreshold(s,dst,50,255,CV_THRESH_BINARY );cvReleaseImage(&r);
cvReleaseImage(&g);
cvReleaseImage(&b);
cvReleaseImage(&s);}
int main(int argc,char** argv)
{
CvSeq * contour=NULL;
CvSeq * contourMax=NULL;
CvMemStorage * storage=cvCreateMemStorage();
cvNamedWindow("aaa",1);
IplImage* src=cvLoadImage("1.png");
IplImage* dst=cvCreateImage(cvGetSize(src),src->depth,1);
IplImage * tempdst1=cvCreateImage(cvGetSize(src),src->depth,3);
IplImage *pCannyImg=cvCreateImage(cvGetSize(src),src->depth,1);
sum_rgb(src,dst);cvErode(dst, dst, 0, 3);
cvZero(tempdst1);
int maxh = 0,maxw=0;
int minh=100000,minw=100000;
//取边缘 以便轮廓提取
cvCanny(dst,pCannyImg,0,10,3);
//轮廓提取
cvFindContours(pCannyImg,storage,&contour,sizeof(CvContour),CV_RETR_CCOMP ,CV_CHAIN_APPROX_SIMPLE);
cvZero(tempdst1);
double areaMax=cvContourArea( contour,CV_WHOLE_SEQ );
for( ; contour != NULL; contour = contour->h_next)
{ //轮廓的方向影响面积的符号。因此函数也许会返回负的结果。应用函数 fabs() 得到面积的绝对值。
double area = cvContourArea( contour,CV_WHOLE_SEQ );
//计算整个轮廓或部分轮廓的面积
if(fabs(area) >= areaMax)
{
areaMax=area;
contourMax=contour;
}
}//把黑线框出来
CvScalar color = CV_RGB(255, 0, 0);
CvRect rect = cvBoundingRect(contourMax,0);
cvRectangle(tempdst1, cvPoint(rect.x, rect.y), cvPoint(rect.x + rect.width, rect.y + rect.height),CV_RGB(255,0,0), 3, 8, 0);
IplImage* temp2 = cvCreateImage(cvGetSize(tempdst1), IPL_DEPTH_8U, 1); //创建目标图像
cvCvtColor(tempdst1,temp2,CV_BGR2GRAY); //cvCvtColor(src,des,CV_BGR2GRAY)
cvAddWeighted(src,1,tempdst1,0.9,0,src);
cvReleaseMemStorage(&storage);
storage=NULL;
cvShowImage("aaa",src);
cvWaitKey(0);
cvDestroyAllWindows();
cvReleaseImage(&src);
cvReleaseImage(&dst);
return 0;
}
演示下结果:
嗯,应该算完成了任务了吧~~~嘎嘎~~~
本文来自互联网用户投稿,文章观点仅代表作者本人,不代表本站立场,不承担相关法律责任。如若转载,请注明出处。 如若内容造成侵权/违法违规/事实不符,请点击【内容举报】进行投诉反馈!
