7

我尝试为做中值滤波器制作python程序。我得到了这篇文章http://www.programming-techniques.com/2013/02/median-filter-using-c-and-opencv-image.html,所以我尝试将该代码转换为 python 代码。

这是python中的代码

from cv2 import * #Import functions from OpenCV
import cv2

if __name__ == '__main__':
    source = cv2.imread("Medianfilterp.png", CV_LOAD_IMAGE_GRAYSCALE)
    final = source[:]
    for y in range(len(source)):
        for x in range(y):
            final[y,x]=source[y,x]

    members=[source[0,0]]*9
    for y in range(1,len(source)-1):
        for x in range(1,y-1):
            members[0] = source[y-1,x-1]
            members[1] = source[y,x-1]
            members[2] = source[y+1,x-1]
            members[3] = source[y-1,x]
            members[4] = source[y,x]
            members[5] = source[y+1,x]
            members[6] = source[y-1,x+1]
            members[7] = source[y,x+1]
            members[8] = source[y+1,x+1]

            members.sort()
            final[y,x]=members[4]

    cv.NamedWindow('Source_Picture', cv.CV_WINDOW_AUTOSIZE)
    cv.NamedWindow('Final_Picture', cv.CV_WINDOW_AUTOSIZE)
    cv2.imshow('Source_Picture', source) #Show the image
    cv2.imshow('Final_Picture', final) #Show the image
    cv2.waitKey()

这是中值滤波前的图片: 源图片

但我得到了奇怪的结果,程序的结果: 最终图片

4

1 回答 1

18

首先,我建议你不要重新发明轮子。OpenCV 已经包含了一种执行中值滤波的方法:

final = cv2.medianBlur(source, 3)

也就是说,您的实现问题在于您的迭代界限。你的y范围是正确的。但是,for x in range(1,y-1):仅迭代到当前y值,而不是图像的整个x范围。这就解释了为什么过滤器只应用于图像左下角的三角形区域。您可以使用shape图像的字段(实际上只是一个 numpy 数组)来获取图像尺寸,然后可以对其进行迭代:

for y in range(1,source.shape[0]-1):
    for x in range(1,source.shape[1]-1):

这会将过滤器应用于整个图像:

中值过滤结果

于 2013-08-26T17:13:24.127 回答