ディジタル画像処理~pythonによる空間フィルタリングpart2~

今回はフーリエ変換を用いた周波数領域でのフィルタリングをpythonで行ってみました。

ローパスフィルタをかけたimg01.png

上絵のようにローパスフィルタをかけた画像を得ることができました。先述通り、平滑化されてぼんやりとした画像が得られます。これらのコードを図式化すると下絵のようになります。

2. ハイパスフィルタ

次にハイパスフィルタをかけていきます。ローパスフィルタと同様にフーリエ変換した画像にデザインしたフィルタをかけます。フーリエ変換までの過程は省略して、フィルタの設計から書きます。

size = gray.shape
mask = np.zeros(size)
length = size[0]
centery = size[0]/2
R = 10
for x in range(0,length):
    for y in range(0,length):
        if (x- centery)**2 +(y- centery)**2 >R**2:
            mask[x,y]=1

Design filter.py

これをフーリエ変換した先ほどの画像にかけていきます。

img_fft2_mask = img_fft2*mask
img_fft2_mask= np.fft.fftshift(img_fft2_mask) 
ifimage = np.fft.ifft2(img_fft2_mask) 
ifimage = ifimage.real
Image.fromarray(np.uint8(ifimage)).show()

Apply filter.py

これによって、ハイパスフィルタをかけた画像を得ることができます。