Основы фильтрации в частотной области

1. Исходное изображение умножается на , чтобы его преобразование Фурье оказалось центрированным.

2. Вычисляется прямое ДПФ изображения, полученного после шага 1.

3. Функция умножается на функцию фильтра .

4. Вычисляется обратное ДПФ от результата шага 3.

5. Выделяется вещественная часть результата шага 4.

6. Результат шага 5 умножается на .

Пусть обозначает входное изображение после шага 1, - его фурье-образ. Тогда фурье-образ выходного изображения определяется выражением:

.

Умножение функций двух переменных и осуществляется поэлементно, т.е. первый элемент функции умножается на первый элемент функции , второй элемент функции умножается на второй элемент функции и т.д. В общем случае компоненты фильтра являются комплексными величинами. На практике чаще всего используются с действительными компонентами. В этом случае и действительная и мнимая части функции умножаются на одну и ту же действительную функцию фильтра . Такие фильтры называются фильтрами нулевого фазового сдвига, поскольку не меняют фазу фурье-преобразования.

Фильтрованное изображение получается вычислением обратного преобразования Фурье от фурье-образа :

Фильтрованное изображение = F^-1 .

Искомое изображение получается выделением действительной части из последнего результата и умножения на , чтобы скомпенсировать эффект от умножения входного изображения на ту же величину. Обратное фурье-преобразование в общем случае является комплексным. Однако в случае вещественного входного изображения и вещественной функции фильтра мнимые части всех значений обратного фурье-преобразования должны равняться 0. Но на практике значения обратного фурье-преобразования, как правило, содержат паразитную мнимую составляющую, что связано с ошибками округлений при вычислениях. Этой составляющей необходимо пренебречь.

Рис.1.7. Основные этапы фильтрации в частотной области