文档章节

图像的膨胀与腐蚀——OpenCV与C++的具体实现

o
 osc_w9s1w4o0
发布于 2019/03/31 23:01
字数 1142
阅读 12
收藏 0

精选30+云产品,助力企业轻松上云!>>>

[toc]

1. 膨胀与腐蚀的原理

膨胀与腐蚀是数学形态学在图像处理中最基础的操作。在笔者之前的文章《图像的卷积(滤波)运算(一)——图像梯度》《图像的卷积(滤波)运算(二)——高斯滤波》具体介绍了图像卷积\滤波的具体的概念与操作,图像的膨胀与腐蚀其实也是一种类似的卷积操作。其卷积操作非常简单,对于图像的每个像素,取其一定的邻域,计算最大值/最小值作为新图像对应像素位置的像素值。其中,取最大值就是膨胀,取最小值就是腐蚀。

2. 膨胀的具体实现

1) OpenCV实现

在OpenCV中实现了图像膨胀的函数dilate(),可以直接调用:

Mat img = imread(imagename, IMREAD_GRAYSCALE);
if (img.empty())
{
	fprintf(stderr, "Can not load image %s\n", imagename);
	return -1;
}

//OpenCV方法
Mat dilated_cv;
dilate(img, dilated_cv, Mat());

dilate()函数第一个参数表示输入影像,第二个参数表示输出影像,第三个表示一个默认的核,在3X3的范围内寻找最大值。

2) C/C++实现

在一般的图像处理时,图像读写是由专门的组件进行读取的。这这里仍然使用OpenCV进行读取,以免增加复杂性。而在CV::Mat类中,提供了at()函数访问某一行某一列的像素值,可以通过at()函数去访问每一个像素的领域。

与之前OpenCV实现的一样,对于每一个像素,遍历以其像素位置为中心的3X3邻域,取最大值作为新图像对应位置的像素值。 其具体实现如下:

//从文件中读取成灰度图像
const char* imagename = "D:\\Data\\imgDemo\\lena.jpg";
Mat img = imread(imagename, IMREAD_GRAYSCALE);
if (img.empty())
{
	fprintf(stderr, "Can not load image %s\n", imagename);
	return -1;
}

//自定义方法
Mat dilated_my;
dilated_my.create(img.rows, img.cols, CV_8UC1);
for (int i = 0; i < img.rows; ++i)
{
	for (int j = 0; j < img.cols; ++j)
	{	
		//uchar minV = 255;
		uchar maxV = 0;
	
		//遍历周围最大像素值
		for (int yi = i-1; yi <= i+1; yi++)
		{
			for (int xi = j-1; xi <= j+1; xi++)
			{					
				if (xi<0||xi>= img.cols|| yi<0 || yi >= img.rows)
				{
					continue;
				}					
				//minV = (std::min<uchar>)(minV, img.at<uchar>(yi, xi));
				maxV = (std::max<uchar>)(maxV, img.at<uchar>(yi, xi));			
			}
		}
		dilated_my.at<uchar>(i, j) = maxV;
	}
}	

3) 验证与结果

为了验证自己的算法是否正确,可以通过把两者膨胀的结果通过compare()函数进行比较。compare()函数会逐个比较两者的像素值,如果相同就会返回255(白色),如果不相同就会返回0(黑色)。整个过程的具体实现如下:

#include <iostream>
#include <algorithm>
#include <opencv2\opencv.hpp>

using namespace cv;
using namespace std;

int main()
{
	//从文件中读取成灰度图像
	const char* imagename = "D:\\Data\\imgDemo\\lena.jpg";
	Mat img = imread(imagename, IMREAD_GRAYSCALE);
	if (img.empty())
	{
		fprintf(stderr, "Can not load image %s\n", imagename);
		return -1;
	}
	
	//OpenCV方法
	Mat dilated_cv;
	dilate(img, dilated_cv, Mat());

	//自定义方法
	Mat dilated_my;
	dilated_my.create(img.rows, img.cols, CV_8UC1);
	for (int i = 0; i < img.rows; ++i)
	{
		for (int j = 0; j < img.cols; ++j)
		{	
			//uchar minV = 255;
			uchar maxV = 0;
		
			//遍历周围最大像素值
			for (int yi = i-1; yi <= i+1; yi++)
			{
				for (int xi = j-1; xi <= j+1; xi++)
				{					
					if (xi<0||xi>= img.cols|| yi<0 || yi >= img.rows)
					{
						continue;
					}					
					//minV = (std::min<uchar>)(minV, img.at<uchar>(yi, xi));
					maxV = (std::max<uchar>)(maxV, img.at<uchar>(yi, xi));			
				}
			}
			dilated_my.at<uchar>(i, j) = maxV;
		}
	}	

	//比较两者的结果
	Mat c;
	compare(dilated_cv, dilated_my, c, CMP_EQ);

	//显示
	imshow("原始", img);
	imshow("膨胀_cv", dilated_cv);
	imshow("膨胀_my", dilated_my);
	imshow("比较结果", c);
		
	waitKey();
	
	return 0;
}

其运行结果如下所示。可以发现最后的比较结果是一张白色的图像,说明自己实现的算法是正确的。

3. 腐蚀的具体实现

同样的办法可以实现图像腐蚀的过程,只要将求局部最大值改成局部最小值就可以了。具体实现过程如下:

#include <iostream>
#include <algorithm>
#include <opencv2\opencv.hpp>

using namespace cv;
using namespace std;

int main()
{
	//从文件中读取成灰度图像
	const char* imagename = "D:\\Data\\imgDemo\\lena.jpg";
	Mat img = imread(imagename, IMREAD_GRAYSCALE);
	if (img.empty())
	{
		fprintf(stderr, "Can not load image %s\n", imagename);
		return -1;
	}
	
	//OpenCV方法
	Mat eroded_cv;
	erode(img, eroded_cv, Mat());

	//自定义方法
	Mat eroded_my;
	eroded_my.create(img.cols, img.rows, CV_8UC1);
	for (int i = 0; i < img.rows; ++i)
	{
		for (int j = 0; j < img.cols; ++j)
		{	
			uchar minV = 255;
			//uchar maxV = 0;
		
			//遍历周围最大像素值
			for (int yi = i-1; yi <= i+1; yi++)
			{
				for (int xi = j-1; xi <= j+1; xi++)
				{					
					if (xi<0||xi>= img.cols|| yi<0 || yi >= img.rows)
					{
						continue;
					}					
					minV = (std::min<uchar>)(minV, img.at<uchar>(yi, xi));
					//maxV = (std::max<uchar>)(maxV, img.at<uchar>(yi, xi));			
				}
			}
			eroded_my.at<uchar>(i, j) = minV;
		}
	}	

	//比较两者的结果
	Mat c;
	compare(eroded_cv, eroded_my, c, CMP_EQ);

	//显示
	imshow("原始", img);
	imshow("膨胀_cv", eroded_cv);
	imshow("膨胀_my", eroded_my);
	imshow("比较结果", c);
		
	waitKey();
	
	return 0;
}

其运行结果如下:

o
粉丝 0
博文 500
码字总数 0
作品 0
私信 提问
加载中
请先登录后再评论。
图像的膨胀与腐蚀——OpenCV与C++的具体实现

目录 1. 膨胀与腐蚀的原理 2. 膨胀的具体实现 3. 腐蚀的具体实现 1. 膨胀与腐蚀的原理 膨胀与腐蚀是数学形态学在图像处理中最基础的操作。在笔者之前的文章《图像的卷积(滤波)运算(一)——图...

charlee44
2019/03/31
0
0
[Python图像处理] 十二.图像几何变换之图像仿射变换、图像透视变换和图像校正

版权声明:本文为博主原创文章,转载请注明CSDN博客源地址!共同学习,一起进步~ https://blog.csdn.net/Eastmount/article/details/88679772 该系列文章是讲解Python OpenCV图像处理知识,前...

Eastmount
2019/03/20
0
0
[Python图像处理] 九.形态学之图像开运算、闭运算、梯度运算

版权声明:本文为博主原创文章,转载请注明CSDN博客源地址!共同学习,一起进步~ https://blog.csdn.net/Eastmount/article/details/83651172 该系列文章是讲解Python OpenCV图像处理知识,前...

Eastmount
2018/11/02
0
0
【全解】基于OpenCv的SVM实现车牌检测与识别

都说深度学习的出现极力地打压着传统机器学习算法的地位,作为一个二刷机器学习经典算法的小伙伴告诉你:还真多半是这样,咳,毕竟还是差距较大,深度学习处理真实世界多维度的问题更权威!不...

周小夏(cv调包侠)
05/14
0
0
[Python图像处理] 八.图像腐蚀与图像膨胀

版权声明:本文为博主原创文章,转载请注明CSDN博客源地址!共同学习,一起进步~ https://blog.csdn.net/Eastmount/article/details/83581277 该系列文章是讲解Python OpenCV图像处理知识,前...

Eastmount
2018/10/31
0
0

没有更多内容

加载失败,请刷新页面

加载更多

自从尝了 Rust,Java 突然不香了

Rust 是软件行业中相对而言比较新的一门编程语言,如果从语法上来比较,该语言与 C++ 其实非常类似,但从另一方面而言,Rust 能更高效地提供许多功能来保证性能和安全。而且,Rust 还能在无需...

osc_k3vwonkw
16分钟前
10
0
Java 高级 面试题 及 参考答案

一、面试题基础总结 1、 JVM结构原理、GC工作机制详解 答:具体参照:JVM结构、GC工作机制详解 ,说到GC,记住两点:1、GC是负责回收所有无任何引用对象的内存空间。 注意:垃圾回收回收的是无...

FH-Admin
17分钟前
14
0
机器学习中的AUC-ROC曲线

作者|ANIRUDDHA BHANDARI 编译|VK 来源|Analytics Vidhya AUC-ROC曲线 你已经建立了你的机器学习模型-那么接下来呢?你需要对它进行评估,并验证它有多好(或有多坏),这样你就可以决定是否...

osc_bg8v9gvf
18分钟前
8
0
音视频(消息)应用场景 :连麦交友例子

实现一个小例子: 效果类似唱吧APP里的 连麦交友功能,音视频,IM 及音视频 SDK参考融云服务商。 没有印象的可以搜索 ’连麦’ 关键字在 应用商店下载一款 连麦的软件 体验下 业务方面的需求...

T型人才追梦者
19分钟前
11
0
逛淘宝天猫想到SSO单点登录

我的原文地址:https://mp.weixin.qq.com/s/77xukPDlgkKnYpwu4LrqaA

osc_yy65eb2q
19分钟前
11
0

没有更多内容

加载失败,请刷新页面

加载更多

返回顶部
顶部