opencv 提取曲线并用矩形框框出来

来源：互联网发布：手机电子书软件编辑：程序博客网时间：2024/04/28 11:36

老师提出写个提取裂缝，并用矩形框框出来的小程序，不用多好，能跑起来看看就行，由于裂纹没有现成的，就随便划了黑线做个模仿，由于水平比较菜，就在网上抄些程序改改，拼拼，好歹算搞定了，也算个demo吧，嘿嘿~~~希望对大家有帮助

原图在这：，，，额，真的好大，算了，不弄了，就这样了

下面直接贴代码喽~~~

<pre name="code" class="html">#include "stdafx.h"#include "stdio.h"#include "cv.h"#include "highgui.h"#include <iostream>#include <string>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 分割多通道数组成分离的单通道数组dcvSplit(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;}

演示下结果：

嗯，应该算完成了任务了吧~~~嘎嘎~~~

0 0