在MATLAB里求解最大连通域用bwlabel就行了,在OpenCV中没有看到现成的函数,在网上
找了一些程序,看起来头大,于是自己编了一个,抛砖引玉了这里
本程序的主要结构是:
1.读彩色bmp文件,提取Green通道
2.中值滤波,Otsu二值化
3.形态学操作去除小细节
4.标注连通域
找了一些程序,看起来头大,于是自己编了一个,抛砖引玉了这里
本程序的主要结构是:
1.读彩色bmp文件,提取Green通道
2.中值滤波,Otsu二值化
3.形态学操作去除小细节
4.标注连通域
5.求解最大连通域
#include <cv.h>
#include <highgui.h>
#include <cxcore.h>
#include <vector>
#include <algorithm>
using namespace std;
///
void CMy20111114qian的matlab程序转opencvDlg::OnBnClickedButton1()
{
IplImage *src = cvLoadImage("before.bmp", CV_LOAD_IMAGE_COLOR);//读取图像文件
cvNamedWindow("origin");
cvShowImage("origin", src);//显示原始图像
//摄像机保存的图像是32位的,有R、G、B和Alpha通道
//图像中实际存储顺序是B、G、R
IplImage *blue = cvCreateImage(cvGetSize(src), IPL_DEPTH_8U, 1);//制作一个单通道图像
IplImage *green = cvCreateImage(cvGetSize(src), IPL_DEPTH_8U, 1);//制作一个单通道图像
IplImage *red = cvCreateImage(cvGetSize(src), IPL_DEPTH_8U, 1);//制作一个单通道图像
cvSplit(src, blue, green, red, NULL);//分割颜色通道
cvSmooth(green, green, CV_MEDIAN, 7);//7*7中值滤波
cvThreshold(green, green, 0.0, 255.0, CV_THRESH_BINARY | CV_THRESH_OTSU);//OTSU法二值化
{
IplConvKernel *element = cvCreateStructuringElementEx(5, 5, 0, 0, CV_SHAPE_ELLIPSE);//定义形态学结构指针
cvMorphologyEx(green, green, NULL, element, CV_MOP_OPEN);//开运算,去除比结构元素小的亮点
cvReleaseStructuringElement(&element);
}
cvNamedWindow("binary");
cvShowImage("binary", green);//显示二值化图像
{
int color = 254;// 从254开始,因此连通域不能多于253个
CvSize sz = cvGetSize(green);
int w;
int h;
for (w=0; w<sz.width; w++)
{
for (h=0; h<sz.height; h++)
{
if (color > 0)
{
if (CV_IMAGE_ELEM(green, unsigned char, h, w) == 255)
{
cvFloodFill(green, cvPoint(w,h), CV_RGB( color,color,color));//把各连通域标记上颜色
color--;
}
}
}
}
cvNamedWindow("labeled");
cvShowImage("labeled", green);//显示标记后的图像
int colorsum[255] = {0};
for (w=0; w<sz.width; w++)
{
for (h=0; h<sz.height; h++)
{
if (CV_IMAGE_ELEM(green, unsigned char, h, w) > 0)//不对0值计数,不可能为255
{
colorsum[CV_IMAGE_ELEM(green, unsigned char, h, w)]++;//统计每种颜色的数量
}
}
}
vector<int> v1(colorsum, colorsum+255);//用数组初始化vector
int maxcolorsum = max_element(v1.begin(), v1.end()) - v1.begin();//求出最多数量的颜色
for (w=0; w<sz.width; w++)
{
for (h=0; h<sz.height; h++)
{
if (CV_IMAGE_ELEM(green, unsigned char, h, w) == maxcolorsum)
{
CV_IMAGE_ELEM(green, unsigned char, h, w) = 255;//只把最多数量的颜色标为255
}
else
{
CV_IMAGE_ELEM(green, unsigned char, h, w) = 0;//其他标为0
}
}
}
cvNamedWindow("最大连通域");
cvShowImage("最大连通域", green);//显示最大连通域
}
cvReleaseImage(&src);
cvReleaseImage(&blue);
cvReleaseImage(&green);
cvReleaseImage(&red);
}