python进行图像滤波_计算机视觉实验_图像滤波处理

最新推荐文章于 2023-02-10 15:25:40 发布

weixin_40000999

最新推荐文章于 2023-02-10 15:25:40 发布

阅读量773

点赞数

文章标签： python进行图像滤波

本文链接：https://blog.csdn.net/weixin_40000999/article/details/111515205

版权

本文介绍了使用Python进行图像滤波的实验，包括高斯滤波和平滑处理。首先，详细阐述了如何通过调整高斯函数标准差实现图像的高斯滤波，讨论了滤波窗口大小的计算方法。接着，探讨了联合双边滤波的实现，涉及输入图像的缩放、双线性插值以及滤波器核的计算。实验展示了不同参数设置下滤波效果的变化，对图像处理初学者具有指导意义。

摘要由CSDN通过智能技术生成

实验内容

2-1：实现图像的高斯滤波处理

实验要求：

1)通过调整高斯函数的标准差(sigma)来控制平滑程度；

给定函数：void Gaussian(const MyImage &input, MyImage &output, double sigma);

2)滤波窗口大小取为[6sigma-1]/22+1，[.]表示取整；

3)利用二维高斯函数的行列可分离性进行加速；

先对每行进行一维高斯滤波，再对结果的每列进行同样的一维高斯滤波；

空间滤波=图像卷积；

高斯滤波=以高斯函数为卷积核的图像卷积。

2-2：实现图像的联合双边滤波处理

实验要求：

给定函数：function im = jbf(D,C,w, sigma_f, sigma_g)

其中：D为输入图像；C为引导图像；W为滤波窗口大小；

sigma_f 为spatial kernel标准差；

sigma_g为range kernel 标准差；

给定公式：

其中：p是q的邻域中的一个像素。 f和g是空间和距离内核，通常以高斯的形式表示，I是输入图像。

f 是空间滤波器内核，定义为：exp([-d_f2]/[2*sigma_f2])

其中：d_f 为输入图像I的像素位置差，sigma_f为空间滤波核函数的标准差

g是距离滤波器内核，定义为：exp([-d_g2]/[2*sigma_g2])

其中：d_g为引导图像的像素灰度值差，sigma_g为滤波核函数的标准差

输入图像I下采样(1/2)得到LR低分辨率图像，再由LR图像上采样2倍得到引导图像。

注：图像缩放采用双线性插值。

实验过程：

2.1 高斯滤波处理

使用cv2做高斯模糊，只要一行代码调用GaussianBlur函数，给出高斯矩阵的尺寸和标准差就可以，也可以自己构造高斯核，相关函数为：cv2.GaussianKernel()。

根据要求sigma自己规定，窗口大小用[6*sigma-1]/2 *2+1计算，刚开始在想这个为什么要先除2再乘2，结果在python中因为自动转换类型没有区别，但是在C++里就不一样了，他会使这个式子永远是基数，所以在python中跑出来报错，输入参数只能是大于0 的奇数；修改后可以正常运行，代码如下：

img = cv2.imread(r"F:\CVcode\lab21\a.jpg")

sigma = float(input())

ksize = int((int(6*sigma - 1)) /2 )*2 + 1

if sigma <= 0: #扩展的公式，如果 sigma 非法

sigma = int(0.3*((ksize-1)*0.5 - 1) + 0.8)

blur = cv2.GaussianBlur(img,(ksize,ksize),sigma)

while(1):

cv2.imshow('image',img)

cv2.imshow('blur',blur)

k=cv2.waitKey(1)

if k == ord('q'):#按q键退出

break

cv2.destroyAllWindows()1

当Alpha = 0.8时ksize = 3，变化不是很大：

当Alpha = 1.5 时ksize = 9，有明显模糊感：

趋势是高斯矩阵尺寸和标准差越大，处理后的图片越模糊；

接下来利用二维高斯函数的行列可分离性进行加速(高斯卷积算子是可分离卷积核)：

高斯卷积核可以分成一维水平方向上的高斯核和一维垂直方向上的高斯核，或者反过来，即满足交换律。使用Mat getGaussianKernel(int ksize,doable sigma,int ktype = CV_64F)返回ksize*1的垂直方向上的高斯核，而对于一位数水平方向上的高斯核，只需要对垂直方向上的转置处理即可。

def gaussBlur(image,sigma,H,W,_boundary = 'fill',_fillvalue = 0):

'''

通过内置函数得到高斯算子

'''

gaussKenrnel_x = np.array(cv2.getGaussianKernel(sigma,1,W))

gaussBlur_x = signal.convolve2d(image,gaussKenrnel_x,mode='same',boundary = _boundary,fillvalue=_fillvalue)

gaussKenrnel_y = np.array(cv2.getGaussianKernel(sigma,H,1))

gaussBlur_xy = signal.convolve2d(gaussBlur_x,gaussKenrnel_y,mode='same',boundary = _boundary,fillvalue=_fillvalue)

return gaussBlur_xy1

但根据老师给的实验要求好像是要自己写函数得到高斯平滑算子，下面是代码：

def getGaussKernel(sigma,H,W):

# 第一步：构建高斯矩阵gaussMatrix