フーリエ変換をして、画像の特定の周波数成分を除去するプログラムの紹介（Pythonでの実装）

本日は画像のフーリエ変換を用いて画像を周波数成分に変換し、特定の周波数成分をカットしてから元に戻すことで、画像を変化させてみます。

画像のフーリエ変換

フーリエ変換は、ある時間領域の信号を周波数領域の信号に変換する変換手法です。

画像に対しては2次元離散フーリエ変換が使われますが、これを高速に計算する手段としてFFT（Fast Fourier Transform：高速フーリエ変換）が広く用いられています。

一般的な画像は、低周波成分が多く、高周波成分が少ない傾向にあります。

よって低周波成分だけを残しても、画像の多くの特徴は残ります。よって、低周波成分だけを残して保存することで、画像のデータサイズを圧縮することができます。

また、一般にホワイトノイズと呼ばれるようなノイズは、低周波・高周波によらず全周波数帯に一様にノイズ成分が載る特徴があります。

よって、高周波成分を積極的に削減することによって、ノイズの成分のみを多く削減することができます。

このように、信号の見方を変えて周波数成分で扱うことは、画像処理において非常に強力なツールとなり得ます。

以下に画像を入力にフーリエ変換を施し、特定の周波数を除去するプログラムを紹介します。今回はグレースケール画像に対してフーリエ変換を実施します。

上記の箇所で画像に対してFFTを実施し、周波数成分へと変換します。

デフォルトだと低周波成分が四隅に現れるので、fftshiftを用いて低周波成分が画像の中央に来るようにシフトします。

上記の箇所で特定の周波数成分をマスクします。

「CUTOFF_FREQ = 50」が残す周波数成分の閾値を表しており、この値が大きいほど高周波成分まで残ります。

周波数成分をシフトさせていたのを元に戻します。

上記の箇所で周波数成分に対して逆変換を実施し、元の画像へと戻します。

Frequency Filterの白の部分が通過させる周波数、黒の部分が遮断する周波数になります。

わずかな低周波成分のみを通過させる場合は、かなり元の画像に劣化が加わっていることがわかります。

一方、ある程度低周波成分を通過させた場合には、かなり元の画像の特徴を保っていることがわかります。

このことは、画像を周波数領域に変換し、低周波数成分のみを保存することで、画像の大部分の特徴を保持しながら、データサイズを削減できることを示します。

本日は画像のフーリエ変換を行い、特定の周波数成分を遮断してみました。

低周波成分のみを残すことでノイズ除去をしたり、高周波成分のみを残すことでエッジを抽出したり、フーリエ変換はさまざまな用途に応用することができます。