ГУАП

КАФЕДРА № 41

ОТЧЕТ ЗАЩИЩЕН С ОЦЕНКОЙ

ПРЕПОДАВАТЕЛЬ

доц., к.т.н., доц.				О.О. Жаринов
должность, уч. степень, звание		подпись, дата		инициалы, фамилия

ОТЧЕТ О ЛАБОРАТОРНОЙ РАБОТЕ №6

Основы цифровой фильтрации изображений. Метод фильтрации в спектральном пространстве

по курсу: МУЛЬТИМЕДИА ТЕХНОЛОГИИ

РАБОТУ ВЫПОЛНИЛ

СТУДЕНТ гр. №	4116
			подпись, дата		инициалы, фамилия

Санкт-Петербург 2025

Цель работы: изучить основы обработки изображений на примере методов фильтрации, реализуемых в спектральном пространстве.

Краткие теоретические сведения

Фильтрация является одним из широко применяемых методов обработки изображений. Концептуальные идеи, положенные в основу методов фильтрации изображений, опираются на методологическую базу обработки одномерных сигналов, таких как аудио. В отличие от задач фильтрации одномерных оцифрованных сигналов, которые представляют собой последовательность числовых значений, упорядоченную по времени появления каждого из них, последовательность значений пикселей изображения упорядочена лишь по возрастанию индексов строк и столбцов соответствующих матриц. Соответственно, если для временных сигналов период измеряется в секундах, а частота в Герцах ([Гц]=[1/c]), то применительно к изображениям период повторения группы пикселей и соответствующая периоду повторения пространственная частота (англ. space frequency) выражаются в условных безразмерных единицах. Для обработки изображений также успешно применяются подходы, связанные с избирательным воздействием на определенные частотные компоненты с целью целенаправленного изменения свойств изображения.

Методы фильтрации чаще всего используются для решения задачи улучшения качества изображений (увеличения резкости расфокусированных и смазанных изображений). Кроме того, фильтрация часто является одним из этапов более сложных, многоступенчатых методов обработки изображений, таких, например, как методы выделения границ объектов на изображениях.

Фильтрация всегда основана на наличии существенных различий в спектральном составе различных компонентов, составляющих исходное изображение. Для получения информации о частотном составе цифровых изображений существует специальный метод и алгоритм – двумерное быстрое преобразование Фурье (2D-БПФ). В результате его применения к массиву { } образуется массив значений его спектра { }. Существует также и алгоритм обратного действия (обратное 2D-БПФ): по массиву значений спектра { } можно вычислить массив значений пикселей изображения { }. В этом смысле концепция выглядит так же, как и в случае одномерного БПФ для обработки аудиосигналов. Расчет двумерного БПФ реализован в методе fft2() библиотеки numpy.

При реализации фильтрации изображений часто применяют метод фильтрации в спектральном пространстве (метод Фурье-фильтрации). Идея метода заключается в вычислении двумерного спектра { } исходного изображения { }, последующей модификации спектра в соответствии заданными частотными характеристиками пространственного фильтра и формировании выходного изображения посредством вычисления обратного двумерного преобразования Фурье модифицированного спектра { }.

Для модификации спектра исходного изображения необходимо задать массив отсчетов двумерной дискретной передаточной функции фильтра (его еще называют частотной маской), после чего выполнить преобразование

= для всех n и m. (6.1)

При задании массива исходят из требований к свойствам, которые необходимо придать выходному изображению.

Фундаментом теории частотной фильтрации изображений является общая идея о том, что относительно крупным объектам на изображении соответствуют более низкие пространственные частоты, а мелким объектам и деталям – высокие частоты. Резкая граница между относительно ярким объектом и более темным окружающим фоном по своим частотным свойствам аналогична фронту одномерного импульсного сигнала, спектр такого фрагмента аудио имеет широкий спектр, в котором много высокочастотных составляющих.

Свойство комплексно-сопряженной симметрии двумерного БПФ требует обязательного выполнения такого же условия симметрии для массива . В двумерном случае симметрия должна быть относительно точки, расположенной в центре, что программно реализуется сложнее, чем в одномерном случае. Для облегчения реализации принципа симметрии в библиотечном модуле numpy имеется специальный метод центрирования спектра fftshift(). При этом осуществляется двумерный циклический сдвиг спектра, таким образом, что компоненты спектра с нулевыми и/или максимальными индексами сдвигаются в центр матрицы спектра. Тогда при необходимости остаются в изображении только низкие частоты (т.е. остаются неизмененными протяженные объекты и уменьшается заметность мелких объектов), необходимо оставлять в спектре только двумерную область, близкую к центру спектра, симметричную относительно центра (очевидно, что нужно быть очень аккуратным, т.к. в программе важно предусмотреть различные ветки кода для четных и нечетных размеров изображений по каждому размеру). Если же, наоборот, требуется выделить края (edges) объектов, а сами объекты сделать незаметными, то нужно оставить только высокие частоты, для чего следует уменьшить или вовсе сделать равными нулю значения спектра, находящиеся близко к центру матрицы. По окончании модификаций центрированного спектра обязательно осуществляется его обратный сдвиг (методом ifftshift()) и вычисляется двумерное обратное БПФ, в результате чего формируется отфильтрованное изображение, свойства которого нужным образом отличаются от исходного.