形态学概念介绍
形态学现在学完基本的几个了,但我还是不知道什么是形态学!原理其实就是和“卷积”在图像处理中的应用一样,就是一个“内核”遍历图像之后进行处理,内核的不同使得处理得到的图像效果也是不同的。下面介绍几种形态学滤波原理你就懂了:
注->RGB:0-255,0代表黑色,255代表白色
腐蚀:腐蚀的顾名思义就是一个东西变黑变坏变烂了,那么简单的理解就是把大于0的像素都都变得接近0就行了啊!
那对应的图像处理:
膨胀:一个人膨胀了的样子怎样的?变得越来越耀眼、越来越明亮!那么对应的像素就是像素<255的就越来越接近255啊。
那对应的图像处理:
开运算:从名字记忆是图像打开,既然是打开那就是最后的结果是膨胀-------------------->先腐蚀后膨胀
闭运算:和开运算相对,从名字记忆是图像关闭,既然是打开那就是最后的结果是腐蚀----->先膨胀后腐蚀
形态学梯度:梯度就是一个阶梯的长度,对应于图像那就是像素的差值,膨胀—原图、原图—腐蚀、膨胀—腐蚀、X/Y等方向的
顶帽:不解释了---->原图—开运算
黑帽:------------>原图—闭运算
---实例分析---
注意点:形态学滤波一般运用在二值化的图像上,对于那些彩色的图像运用不明显(用过之后很难看),看了很多书本的介绍都是随便找个例子,这是在课程中看到的,感觉按照下面的步骤学习形态学真的很简单而且实用!
例一:腐蚀的作用
原图如下,去除图片上的小白点。
用内核大小3X3进行的图片:小的白点已经没有了,但是稍微大点的杂点还是没去除!
这是实用15X15的内核进行的图片:图片的白点完全去除了。
这是不是完成了我们的要求了呢?仔细的看会发现,我们想要的大白色区域变小了,这是什么原因呢?
从我们上面的原理分析可以得知:腐蚀会把目标区域给变小的,请看下面的图片->>>红色区域是内核,
一号区域->黑色
二号区域->白色
三号区域->黑色
二号区域->黑色
五号区域->黑色
所以图片缩小的区域就是内核的大小,每个边都会缩小!
在想一下,如果我们用膨胀处理经过腐蚀的图片会怎么样?由上面的分析可以很快得到结论,就是恢复我们目标区域的原始尺寸
看下面的效果图:
代码比较简单,就是几行API,但是如何运用,为什么这么运用,这才是关键:
1 int main(int argc,char**argv)2 {3 Mat input_image = imread("1.jpg");4 if (input_image.data==NULL) {5 return -1; cout << "can't open image.../";6 }7 imshow("Sourse image", input_image);8 Mat output_image;9 Mat kernel = getStructuringElement(MORPH_RECT,Size(15,15));
10 erode(input_image,output_image,kernel);
11 dilate(output_image, output_image,kernel);
12 imshow("Destinate image",output_image);
13 waitKey(0);
14 return 0;
15 }
例二:提取行和列
要求提取其中的行线段-->>
经过处理的线段图片:
这里没给其他的特殊照片,看代码直接改一下就可以了。
注意点: getStructuringElement()获得内核的一些参数->核大小、核形状、核锚点等。其中控制核的大小可以滤波不同的噪点:
我要滤去下面的三个大噪点,保留上面的大白色区域,其实滤波的核定义成红色的大小就可以了,不一定是正方形,矩形就可以了。
int main(int argc,char*=*argv)2 {3 Mat input_image = imread("1.jpg");4 if (input_image.data==NULL) {5 return -1; cout << "can't open image.../";6 }7 imshow("Sourse image", input_image);8 Mat output_image;9 Mat kernel = getStructuringElement(MORPH_RECT,Size(input_image.cols/30,1));//这个input_image.cols/30,是定义核的长度是图片长度的三十分之一,如果直 接给定一个数200也可以,但是你不知道200在图像上是多大啊。
10 erode(input_image,output_image,kernel);
11 dilate(output_image, output_image,kernel);
12 imshow("Destinate image",output_image);
13 waitKey(0);
14 return 0;
15 }
例三:简单的提取字母
目的是提取图片中的字母
阈值化:
形态学滤波:
取反之后:
1 int main(int argc,char**argv)2 {3 Mat input_image = imread("1.jpg");4 if (input_image.data==NULL) {5 return -1; cout << "can't open image.../";6 }7 imshow("Sourse image", input_image);8 cvtColor(input_image,input_image,CV_RGB2GRAY);9 imshow("Sourse1 image", input_image);
10 //adaptiveThreshold(input_image,input_image,255,ADAPTIVE_THRESH_GAUSSIAN_C,THRESH_BINARY_INV,171,0);
11 threshold(input_image,input_image,15,255,THRESH_BINARY);
12 imshow("Sourse2 image", input_image);
13 Mat output_image;
14 Mat kernel = getStructuringElement(MORPH_RECT,Size(3,3));
15 erode(input_image,output_image,kernel);
16 dilate(output_image, output_image,kernel);
17 imshow("Destinate image",output_image);
18 bitwise_not(output_image,output_image);
19 imshow("Destinate2 image", output_image);
20 Mat my_kernel = (Mat_<uchar>(3, 3) << 0, -1, 0, 5, -1, 0, -1, 0);
21 filter2D(output_image,output_image,output_image.depth(),my_kernel);//加强显示
22 imshow("Destinate3 image", output_image);
23 waitKey(0);
24 return 0;
25 }
例四:稍微困难的提取字母
这个图像对我来说有点麻烦的,形态学滤波不行的,而且形态学操作之后留下很多噪点。。。。。
形态学操作之后:
轮廓检测去除噪点:
霍夫变换去除粗实线:
这个图的小噪点用上面的步骤可以去除,这个就没再继续了
上代码:
1 Mat input_image = imread("2.jpg");2 if (input_image.data==NULL) {3 return -1; cout << "can't open image.../";4 }5 imshow("Sourse image", input_image);6 cvtColor(input_image,input_image,CV_RGB2GRAY);7 //adaptiveThreshold(input_image,input_image,255,ADAPTIVE_THRESH_GAUSSIAN_C,THRESH_BINARY_INV,171,0);8 threshold(input_image,input_image,0,255,THRESH_BINARY|THRESH_OTSU);9 //-----------------------去除细实线------------------------//
10 Mat output_image;
11 Mat kernel1 = getStructuringElement(MORPH_RECT,Size(3,3));
12 morphologyEx(input_image, input_image, MORPH_CLOSE, kernel1);
13 bitwise_not(input_image, input_image);
14 output_image = input_image.clone();
15 imshow("DeleteThick image", input_image);
16 //----------------------去除形态学不能去除的噪点----------------------//
17 vector<vector<Point> > contours;
18 vector<Vec4i> hierarchy;
19 findContours(output_image, contours, hierarchy, RETR_TREE, CHAIN_APPROX_SIMPLE);
20 Mat contours_image = Mat::zeros(input_image.size(), input_image.type());
21 for (size_t i = 0; i < contours.size(); i++)
22 {
23 double Area = contourArea(contours[i]);
24 if (Area > 50) continue;
25 drawContours(contours_image, contours, static_cast<int>(i),Scalar(255,255,255),1);
26 }
27 input_image = input_image - contours_image;
28 morphologyEx(input_image, input_image, MORPH_OPEN, kernel1);
29 imshow("contours image", input_image);
30 //------------------去除粗实线--------------------//
31 vector<Vec4i> lines;
32 HoughLinesP(input_image,lines,1,CV_PI/180,100,0,200);
33 Mat line_image = Mat::zeros(input_image.size(), input_image.type());
34 for (size_t i = 0; i < lines.size(); i++)
35 {
36 line(line_image, Point(lines[i][0], lines[i][1]),
37 Point(lines[i][2], lines[i][3]), Scalar(255, 255, 255), 1, 8);
38 }
39 bitwise_not(input_image, input_image);
40 input_image = line_image + input_image;
41 morphologyEx(input_image, input_image, MORPH_CLOSE, kernel1);
42 imshow("Last image", input_image);