python 图像滤波算法_[Python图像处理]四.图像平滑中四种常用的滤波

主要讲解Python调用OpenCV实现图像平滑,包括四个算法:均值滤波、方框滤波、高斯滤波和中值滤波.

给图像增加噪声:

importcv2importnumpy as npdeftest10():

img= cv2.imread("result.jpg", cv2.IMREAD_UNCHANGED)

rows, cols, chn=img.shape#加噪声

for i in range(5000):

x=np.random.randint(0, rows)

y=np.random.randint(0, cols)

img[x, y, :]= 255cv2.imshow("demo", img)if cv2.waitKey(0) == 27:

cv2.destroyWindow("demo")

test10()

效果如下:

1108198-20200510192524244-1316286207.png

均值滤波: 均值滤波是指任意一点的像素值,都是周围N*M个像素值的均值, result = cv2.blur(图像, 核大小),其中核大小是以(宽度, 高度)表示的元组形式,常见的形式包括: (3, 3)和(5, 5)

代码如下:

importmatplotlib.pyplot as pltimportcv2deftest11():

img= cv2.imread("result.jpg")

source=cv2.cvtColor(img, cv2.COLOR_BGR2BGRA)#均值滤波

result = cv2.blur(source, (5, 5))#显示图像

title = ["demo1", "demo2"]

imgs=[source, result]for i in range(2):

plt.subplot(1, 2, i+1)

plt.imshow(imgs[i], "gray")

plt.title(title[i])

plt.xticks([])

plt.yticks([])

plt.show()

test11()

效果如下:

1108198-20200510194244397-529668009.png

将核值调大会让图像变模糊,例如设置为(20, 20),效果如下:如果设置为(1, 1)处理的结果就是原图,核中每个权重值相同,称为均值。

1108198-20200510194421165-147296074.png

方框滤波

方框滤波和均值滤波核基本一致,区别是需不需要均一化处理。OpenCV调用boxFilter()函数实现方框滤波。函数如下:

result = cv2.boxFilter(原始图像, 目标图像深度, 核大小, normalize属性) 其中,目标图像深度是int类型,通常用“-1”表示与原始图像一致;核大小主要包括(3,3)和(5,5);normalize属性表示是否对目标图像进行归一化处理。当normalize为true时需要执行均值化处理,当normalize为false时,不进行均值化处理,实际上为求周围各像素的和,很容易发生溢出,溢出时均为白色,对应像素值为255。

代码如下:

importcv2importmatplotlib.pyplot as pltdeftest12():

img= cv2.imread("result.jpg")

source=cv2.cvtColor(img, cv2.COLOR_BGR2BGRA)#方框滤波

result = cv2.boxFilter(source, -2, (5, 5), normalize=1)#显示图像

title = ["demo1", "demo2"]

imgs=[source, result]for i in range(2):

plt.subplot(1, 2, i+1), plt.imshow(imgs[i], "gray")

plt.title(title[i])

plt.xticks([])

plt.yticks([])

plt.show()

test12()

效果如下:

1108198-20200510200226005-701112295.png

如果省略参数normalize,则默认是进行归一化处理。如果normalize=0则不进行归一化处理,像素值为周围像素之和,图像更多为白色。

效果如下:

1108198-20200510200413485-921805040.png

上图很多像素为白色,因为图像求和结果几乎都是255。如果设置的是2*2矩阵,只取四个像素结果要好些。 result = cv2.boxFilter(source, -1, (2,2), normalize=0)

1108198-20200510200636960-1186387801.png

高斯滤波: 高斯滤波让临近的像素具有更高的重要度,对周围像素计算加权平均值,较近的像素具有较大的权重值。dst = cv2.GaussianBlur(src, ksize, sigmaX)其中,src表示原始图像,ksize表示核大小,sigmaX表示X方向方差。注意,核大小(N, N)必须是奇数,X方向方差主要控制权重

代码如下:

importcv2importmatplotlib.pyplot as pltdeftest12():

img= cv2.imread("result.jpg")

source=cv2.cvtColor(img, cv2.COLOR_BGR2BGRA)#高斯滤波

result = cv2.GaussianBlur(source, (3, 3), 0)#显示图像

title = ["demo1", "demo2"]

imgs=[source, result]for i in range(2):

plt.subplot(1, 2, i+1), plt.imshow(imgs[i], "gray")

plt.title(title[i])

plt.xticks([])

plt.yticks([])

plt.show()

test12()

效果如下:

1108198-20200510201406255-691682109.png

核数增加,图像会变模糊,例如设置为(17, 17)

1108198-20200510201505753-1966967340.png

中值滤波: 在使用邻域平均法去噪的同时也使得边界变得模糊。而中值滤波是非线性的图像处理方法,在去噪的同时可以兼顾到边界信息的保留。选一个含有奇数点的窗口W,将这个窗口在图像上扫描,把窗口中所含的像素点按灰度级的升或降序排列,取位于中间的灰度值来代替该点的灰度值.

OpenCV主要调用medianBlur()函数实现中值滤波,图像平滑里中值滤波的效果最好。

dst = cv2.medianBlur(src, ksize)

其中,src表示源文件,ksize表示核大小。核必须是大于1的奇数,如3、5、7等

代码如下:

importcv2deftest13():

img= cv2.imread("result.jpg")#高斯滤波

result = cv2.medianBlur(img, 5)

cv2.imshow("demo1", img)

cv2.imshow("demo2", result)if cv2.waitKey(0) == 27:

cv2.destroyWindow("demo1")

cv2.destroyWindow("demo2")

test13()

效果如下:

1108198-20200510202312112-50827314.png

转自: https://blog.csdn.net/Eastmount/article/details/82216380