低通滤镜可模糊图像

Question

我正在尝试通过fft来模糊图像，方法是创建一个低通滤波器，但输出会变成充满灰色噪声的图像。我只是在尝试遵循这里的基础知识，但是我的实现似乎有问题：

from scipy import fftpack
import numpy as np
import imageio
from PIL import Image, ImageDraw

image1 = imageio.imread('image.jpg',as_gray=True)

#convert image to numpy array
image1_np=np.array(image)

#fft of image
fft1 = fftpack.fftshift(fftpack.fft2(image1_np))

#Create a low pass filter image
x,y = image1_np.shape[0],image1_np.shape[1]
#size of circle
e_x,e_y=50,50
#create a box 
bbox=((x/2)-(e_x/2),(y/2)-(e_y/2),(x/2)+(e_x/2),(y/2)+(e_y/2))

low_pass=Image.new("L",(image1_np.shape[0],image1_np.shape[1]),color=0)

draw1=ImageDraw.Draw(low_pass)
draw1.ellipse(bbox, fill=255)

low_pass_np=np.array(low_pass)
low_pass_fft=fftpack.fftshift(fftpack.fft2(low_pass))

#multiply both the images
filtered=np.multiply(fft1,low_pass_fft)

#inverse fft
ifft2 = abs(fftpack.ifft2(fftpack.ifftshift(filtered)))

#save the image
imageio.imsave('fft-then-ifft.png', ifft2.astype(np .uint8))

Answer 1

如Cris Luengo的评论中所述，有几件事需要更正：

• 为低通滤波器提供的椭圆形在频域中很有意义，因此您不应该计算其FFT。

• 255的滤波器大小将结果缩放相同的数量。当您存储如此大的值时，uint8类型会回绕以仅保留8个最低有效位，从而导致看起来像噪声。只需更改过滤器的值即可解决此问题：

  draw1.ellipse(bbox, fill=1)
  

• 重新调整缩放比例后，在图像的某些区域中，计算出的filtered可能仍略微超出所需的0-255范围。这会创建环绕点（由白色像素包围的区域中的黑色区域，由黑色像素包围的区域中的白色区域，甚至是图像从白色到黑色再到白色的渐变带）。为避免这种情况，通常使用以下方法将值裁剪到0-255范围：

  ifft2 = np.real(fftpack.ifft2(fftpack.ifftshift(filtered)))
  ifft2 = np.maximum(0, np.minimum(ifft2, 255))
  

进行这些更正后，您应该具有以下代码：

from scipy import fftpack
import numpy as np
import imageio
from PIL import Image, ImageDraw

image1 = imageio.imread('image.jpg',as_gray=True)

#convert image to numpy array
image1_np=np.array(image1)

#fft of image
fft1 = fftpack.fftshift(fftpack.fft2(image1_np))

#Create a low pass filter image
x,y = image1_np.shape[0],image1_np.shape[1]
#size of circle
e_x,e_y=50,50
#create a box 
bbox=((x/2)-(e_x/2),(y/2)-(e_y/2),(x/2)+(e_x/2),(y/2)+(e_y/2))

low_pass=Image.new("L",(image1_np.shape[0],image1_np.shape[1]),color=0)

draw1=ImageDraw.Draw(low_pass)
draw1.ellipse(bbox, fill=1)

low_pass_np=np.array(low_pass)

#multiply both the images
filtered=np.multiply(fft1,low_pass_np)

#inverse fft
ifft2 = np.real(fftpack.ifft2(fftpack.ifftshift(filtered)))
ifft2 = np.maximum(0, np.minimum(ifft2, 255))

#save the image
imageio.imsave('fft-then-ifft.png', ifft2.astype(np .uint8))

以及以下经过过滤的图像：