use*_*151 2 python image image-processing noise noise-reduction
我尝试实现以下算法,但生成的图像看起来相同。
第 1 步:读取噪声图像。
步骤 2:选择大小为 3x3 的 2D 窗口,中心元素作为处理像素。假设正在处理的像素是 P ij 。
步骤3:如果P ij 是未损坏的像素(即0< P ij <255),则其值保持不变。
步骤4:如果P ij = 0 或P ij = 255,则P ij 是损坏的像素。
第 5 步:如果所选窗口中有 3/4 或更多像素有噪声,则将窗口大小增加到 5x5。步骤 6:如果所选窗口中的所有元素都是 0?s 和 255?s,则将 Pij 替换为窗口中元素的均值,否则转到步骤 7。
步骤 7:从所选窗口中消除 0?s 和 255?s,并找到剩余元素的中值。用中值替换 Pij。
步骤8:重复步骤2到6,直到处理完整幅图像中的所有像素。
这是我的代码。请提出改进建议。
import Image
im=Image.open("no.jpg")
im = im.convert('L')
for i in range(2,im.size[0]-2):
for j in range(2,im.size[1]-2):
b=[]
if im.getpixel((i,j))>0 and im.getpixel((i,j))<255:
pass
elif im.getpixel((i,j))==0 or im.getpixel((i,j))==255:
c=0
for p in range(i-1,i+2):
for q in range(j-1,j+2):
if im.getpixel((p,q))==0 or im.getpixel((p,q))==255:
c=c+1
if c>6:
c=0
for p in range(i-2,i+3):
for q in range(j-2,j+3):
b.append(im.getpixel((p,q)))
if im.getpixel((p,q))==0 or im.getpixel((p,q))==255:
c=c+1
if c==25:
a=sum(b)/25
print a
im.putpixel((i,j),a)
else:
p=[]
for t in b:
if t not in (0,255):
p.append(t)
p.sort()
im.putpixel((i,j),p[len(p)/2])
else:
b1=[]
for p in range(i-1,i+2):
for q in range(j-1,j+2):
b1.append(im.getpixel((p,q)))
im.putpixel((i,j),sum(b1)/9)
im.save("nonoise.jpg")
Run Code Online (Sandbox Code Playgroud)