opencv#30 线性滤波

均值滤波原理

均值滤波步骤

Step1:求和。

Step2:计算平均值。

所谓均值滤波,就是求平均值的意思。我们假设在一个3*3的范围内有一个图像,其中这个图像每一个像素可能含有噪声,也可能不含噪声,我们是不知道的,因此通过均值滤波的方式,对图像中所有像素进行求和,并除以像素个数,得到的结果就是滤波后的结果,将3*3区域的中心位置更改为滤波后的平均值,这样的过程就是均值滤波。

均值滤波和图像的卷积操作相类似,它们都是在原图像中将每一个像素分别覆盖滤波模板,也可称为滤波器,覆盖原图像后求取平均值,放入中心区域作为滤波结果,然后移动滤波器,使滤波器能够覆盖图像的每一个位置。对于上图3*3的矩阵来说,其中每一个位置都乘了1,也就是说它进行卷积的模板是一个3*3的矩阵,此矩阵中的值都为1,之后又除以了和,因此它的卷积模板(滤波器)就是一个1/9乘全为1的矩阵(3*3)。

均值滤波函数

blur()

cpp 复制代码
void cv::blur(InuputArray   src,
              OutputArray   dst,
              Size          ksize,
              Point         anchor = Point(-1,-1),
              int           borderType = BORDER_DEFAULT
             )

**·**src:待均值滤波的图像,图像数据类型必须是CV_8U,CV_16U,CV_16S,CV_32F和CV_64F这五种数据类型之一。均值滤波输出的结果的数据类型与输入图像的数据类型保持一致。

**·**dst:均值滤波后的图像,与输入图像具有相同的尺寸和数据类型。

**·**ksize:卷积核尺寸(滤波器尺寸:常用的是奇数正方形)。因为均值滤波中,滤波器中的每一位参数值是恒定的,因此这里面我们只需给出其尺寸即可。

· anchor:内核的基准点(锚点),其默认值为(-1,-1)代表内核基准点位于kernel的中心位置。**·**borderType:像素外推法选择标志。

方框滤波函数

boxFilter()

cpp 复制代码
void cv::boxFilter(InputArray   src,
                   OutputArray  dst,
                   int          ddepth,
                   Size         ksize,
                   Point        anchor = Point(-1,-1),
                   bool         normalize = true,
                   int          borderType = BORDER_DEFAULT
                  )

方框滤波函数与均值滤波其实是等价的,只不过方框滤波运行滤波器不进行归一化,但是在均值滤波中是强制进行归一化的。 由于方框滤波是一个求和的参数,可能会出现255+255+255+....很多个255相加,最终超出255范围的情况,因此第三个参数运行用户设定输出图像的数据类型,也就是输出图像的数据类型可以与输入图像不同,但是尺寸和通道数是一致的。

·src:输入图像。

·dst:输出图像,与输入图像具有相同的尺寸和通道数。

·ddepth:输出图像的数据类型(深度),根据输入图像的数据类型不同拥有不同的取值范围。

·ksize:卷积核尺寸。

·anchor:内核的基准点(锚点),其默认值为(-1,-1)代表内核基准点位于kernel的中心位置。

·normalize:是否将卷积核进行归一化的标志,默认参数为true,表示进行归一化(在矩阵前面乘所有数值之和)。

·borderType:像素外推法选择标志。

高斯滤波原理

高斯滤波也就是滤波器是一个高斯分布的形式,如上图所示,滤波器中中间的像素值较大,而四周像素值较小的情况,是一个中心对称的滤波器。若以这样的滤波器作为卷积核与原图像进行卷积操作的时候,得到的就是一个高斯滤波之后的结果 。

高斯滤波函数

GaussianBlur()

cpp 复制代码
void cv::GaussianBlur(InputArray    src,
                      OutputArray   dst,
                      Size          xsize,
                      double        sigmaX,
                      double        sigmaY = 0,
                      int           borderType = BORDER_DEFAULT
                     )

·src:待高斯滤波图像,数据类型必须为CV_8U,CV_16U,CV_16S,CV_32F或CV_64F。

·dst:输出图像,与输入图像src具有相同的尺寸,通道数和数据类型。

·ksize:高斯滤波器的尺寸,滤波器可以不为正方形,但是必须是正奇数。如果尺寸为0,则由标准偏差计算尺寸。高斯滤波器不需要人为的输入,因为高斯滤波器的高斯函数分布形式是固定的,只需给出滤波器尺寸即可。

·sigmaX:X方向的高斯滤波器标准偏差。

·sigmaY:Y方向的高斯滤波器标准偏差。如果输入量为0,则将其设置为等于sigmaX,如果两个轴的标准差均为0,则根据输入的高斯滤波器尺寸计算标准差。

ksize ,sigmaX,sigmaY不允许同时为0,若某个参数为0,可以根据对应关系计算出参数。

·borderType:像素外推法选择标志。

示例
cpp 复制代码
#include <opencv2/opencv.hpp>
#include <iostream>

using namespace cv; //opencv的命名空间
using namespace std;



//主函数
int main()
{
	Mat lenaGray = imread("E:/opencv/opencv-4.6.0-vc14_vc15/opencv/lenaGray.png");
	Mat equalLena_Salt = imread("E:/opencv/opencv-4.6.0-vc14_vc15/opencv/equalLena_S.png");
	Mat equalLena_Gauss = imread("E:/opencv/opencv-4.6.0-vc14_vc15/opencv/equalLena_G.png");

	if (lenaGray.empty() || equalLena_Salt.empty() || equalLena_Gauss.empty())
	{
		cout << "请确认图像名称是否正确" << endl;
		return -1;
	}

	Mat result_3, result_9; //存放不含噪声滤波结果,后面的数字代表滤波器尺寸
	Mat result_3gauss, result_9gauss; //存放含有高斯噪声滤波结果,后面数字代表滤波器尺寸
	Mat result_3salt, result_9salt; //存放含有椒盐噪声滤波结果,后面的数字代表滤波器尺寸

	//调用均值滤波函数blur()进行滤波
	blur(lenaGray, result_3, Size(3, 3));
	blur(lenaGray, result_9, Size(9, 9));
	blur(equalLena_Salt, result_3salt, Size(3, 3));
	blur(equalLena_Salt, result_9salt, Size(9, 9));
	blur(equalLena_Gauss, result_3gauss, Size(3, 3));
	blur(equalLena_Gauss, result_9gauss, Size(9, 9));


	//显示不含噪声图像
	imshow("lenaGray", lenaGray);
	imshow("result_3", result_3);
	imshow("result_9", result_9);
	//显示不含噪声图像
	imshow("equalLena_Salt", equalLena_Salt);
	imshow("result_3salt", result_3salt);
	imshow("result_9salt", result_9salt);
	//显示不含噪声图像
	imshow("equalLena_Gauss", equalLena_Gauss);
	imshow("result_3gauss", result_3gauss);
	imshow("result_9gauss", result_9gauss);

	cout << "下面是方框滤波" << endl;
	waitKey(0);//等待函数用于显示图像,按下键盘任意键后退出
	
	Mat resultNorm, result;
	//方框滤波boxFilter()和sqrBoxFilter()
	boxFilter(lenaGray, resultNorm, -1, Size(3, 3), Point(-1, -1), true); //进行归一化
	boxFilter(lenaGray, result, -1, Size(3, 3), Point(-1, -1), false); //进行归一化

	//显示处理结果
	imshow("resultNorm", resultNorm);
	imshow("result", result);
	cout << "下面是高斯滤波" << endl;
	waitKey(0);//等待函数用于显示图像,按下键盘任意键后退出

	Mat result_5_G, result_9_G;//存放不含噪声滤结果,后面数字代表滤波器尺寸
	Mat result_5Gauss_G, result_9Gauss_G;//存放含高斯噪声滤结果,后面数字代表滤波器尺寸
	Mat result_5Salt_G, result_9Salt_G;//存放含椒盐噪声滤结果,后面数字代表滤波器尺寸

	GaussianBlur(lenaGray, result_5_G, Size(5, 5), 10, 20);
	GaussianBlur(lenaGray, result_5_G, Size(5, 5), 10, 20);
	GaussianBlur(equalLena_Gauss, result_5Gauss_G, Size(5, 5), 10, 20);
	GaussianBlur(equalLena_Gauss, result_9Gauss_G, Size(9, 9), 10, 20);
	GaussianBlur(equalLena_Salt, result_5Salt_G, Size(5, 5), 10, 20);
	GaussianBlur(equalLena_Salt, result_9Salt_G, Size(9, 9), 10, 20);


	//显示不含噪声图像
	imshow("lenaGray", lenaGray);
	imshow("result_5_G", result_5_G);
	imshow("result_9_G", result_9_G);
	//显示不含噪声图像
	imshow("equalLena_Salt", equalLena_Salt);
	imshow("result_5Gauss_G", result_5Gauss_G);
	imshow("result_9Gauss_G", result_9Gauss_G);
	//显示不含噪声图像
	imshow("equalLena_Gauss", equalLena_Gauss);
	imshow("result_5Salt_G", result_5Salt_G);
	imshow("result_9Salt_G", result_9Salt_G);
	waitKey(0);//等待函数用于显示图像,按下键盘任意键后退出

	return 0;

}

滤波完成后:

在均值滤波中:对于椒盐噪声和高斯噪声的滤波结果,尺寸越大的滤波器,对于图像来说,滤波后的结果也越模糊,噪声产生的影响变小。如果想最大程度的去除噪声,采用均值滤波的形式的话,那么需要一个很大的尺寸,同时对图像的整体质量有较大的影响。

对椒盐噪声,由于我们的均值滤波是每一个像素值都贡献了相同作用,所以对于椒盐噪声来说,仍然可以看出图像中含有椒盐噪声,只不过椒盐噪声的颗粒没有那么亮,没有那么明显,但是椒盐噪声的颗粒变得更粗了。

在方框滤波中:方框滤波中如果进行归一化操作,所得到的结果与均值滤波的结果一致,若不进行归一化操作,那么有可能出现图像中数值远远大于255,但是在opencv中,会强制为255,所以导致整个图像中,较大的区域都为白色,而有纹理的恰恰是原图像中较暗的区域。

在高斯滤波中:顾名思义对高斯噪声的滤波更为理想。结果也会变得模糊,但是处理结果也较好。若图像中含有椒盐噪声,采用线性滤波的方式得到的结果不太理想。

相关推荐
边缘计算社区31 分钟前
首个!艾灵参编的工业边缘计算国家标准正式发布
大数据·人工智能·边缘计算
游客52042 分钟前
opencv中的各种滤波器简介
图像处理·人工智能·python·opencv·计算机视觉
一位小说男主42 分钟前
编码器与解码器:从‘乱码’到‘通话’
人工智能·深度学习
小俊俊的博客1 小时前
海康RGBD相机使用C++和Opencv采集图像记录
c++·opencv·海康·rgbd相机
KeyPan1 小时前
【IMU:视觉惯性SLAM系统】
计算机视觉
7yewh1 小时前
嵌入式Linux QT+OpenCV基于人脸识别的考勤系统 项目
linux·开发语言·arm开发·驱动开发·qt·opencv·嵌入式linux
深圳南柯电子1 小时前
深圳南柯电子|电子设备EMC测试整改:常见问题与解决方案
人工智能
Kai HVZ1 小时前
《OpenCV计算机视觉》--介绍及基础操作
人工智能·opencv·计算机视觉
biter00881 小时前
opencv(15) OpenCV背景减除器(Background Subtractors)学习
人工智能·opencv·学习
吃个糖糖1 小时前
35 Opencv 亚像素角点检测
人工智能·opencv·计算机视觉