图片噪声去除

# -*- coding: utf-8 -*-
"""
Created on Wed Mar 29 15:32:46 2017

@author: Joker

"""
import matplotlib.pyplot as plt
from skimage import io
im=io.imread('b.jpg')
#io.imshow(img)

a = im.shape[0]
b = im.shape[1]

im_copy = [[0 for x in range(a)] for b in range(b)]
for i in range(0,im.shape[0]):
    for j in range(0,im.shape[1]):
        im_copy[i][j]=im[i][j]
#用3*3的中值滤波器
step=3
def m_filter(x,y):
    sum_s=[]
    for k in range(-int(step/2),int(step/2)+1):
        for m in  range(-int(step/2),int(step/2)+1):
            sum_s.append(im[x+k][y+m])
    sum_s.sort(key=lambda x:x[1])
    sum_s.sort(key=lambda x:x[0])
    #print('sum_s:',sum_s)
    return sum_s[(int(step*step/2)+1)]

for i in range(int(step/2),im.shape[0]-int(step/2)):
    for j in range(int(step/2),im.shape[1]-int(step/2)):
        im_copy[i][j]=m_filter(i,j)

#io.imshow(im_copy)
plt.imshow(im_copy)
plt.show()