Инверсная фильтрация

Поделим равенство (4) на H(u, v) и получаем следующую оценку исходного изображения:

. (5)

Это называется инверсной фильтрацией, но на практике практически никогда не работает. Посмотрим на последнее слагаемое в формуле (5) - если функция H(u,v) принимает значения близкие к нулю или нулевые, то вклад этого слагаемого будет доминирующим. Это практически всегда встречается в реальных примерах - для объяснения этого вспомним, как выглядит спектр после преобразования Фурье.

Берем исходное изображение, преобразуем его в полутоновое и используя Matlab, получаем спектр:

В результате получаем две компоненты: амплитудный и фазовый спектры. Амплитудный спектр показан в логарифмической шкале, т.к. его значения варьируются очень сильно, на несколько порядков. В центре находятся максимальные значения (порядка миллионов), которые быстро убывают по мере удаления от центра практически до нулевых. Именно из-за этого, инверсная фильтрация будет работать только при нулевых или практически нулевых значениях шума. Продемонстрируем это на практике с помощью следующего скрипта:

noise_var = 0.0000001 noise_var = 0.000005

Хорошо видно, что добавление даже очень небольшого шума приводит к значительным помехам, что сильно ограничивает практическое применение метода.

Лучшие статьи по информатике

Проектирование микропроцессорного устройства
Спроектировать микропроцессорное устройство содержащее МП, системный контроллер, адресные буферы, ОЗУ, ПЗУ, порт ввода/вывода, адресный дешифратор. ...

Расчет усилительного каскада на биполярном транзисторе
Транзистор - это полупроводниковый прибор с двумя или несколькими р-n-переходами, позволяющий усиливать электрические сигналы и имеющий три вывода или более ...

Разработка цикловой системы управления промышленным роботом
Электроника - это область науки и техники, которая занимается изучением физических основ функционирования, исследованием, разработкой и применением приборо ...