Почему фильтры Гаусса используются в качестве фильтров нижних частот при обработке изображений?

30

При обработке 1d сигнала используются многие типы фильтров нижних частот. Гауссовы фильтры почти никогда не используются.

Почему они так популярны в приложениях обработки изображений? Являются ли эти фильтры результатом оптимизации какого-либо критерия или являются специальным решением, так как «пропускная способность» изображения обычно не определена правильно.

nimrodm
источник

Ответы:

28

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

Другим важным количественным свойством гауссовых фильтров является то, что они везде неотрицательны . Это важно, потому что большинство 1D-сигналов изменяются примерно на 0 ( ) и могут иметь положительные или отрицательные значения. Изображения отличаются в том смысле, что все значения изображения неотрицательны ( ). Свертка с гауссовым ядром (фильтр) гарантирует неотрицательный результат, поэтому такая функция отображает неотрицательные значения в другие неотрицательные значения ( ). Поэтому результатом всегда является другое действительное изображение.ИксрИкср+е:р+р+

В целом, подавление частоты при обработке изображений не так критично, как в 1D сигналах. Например, в схемах модуляции ваши фильтры должны быть очень точными, чтобы отклонять другие каналы, передаваемые на разных несущих частотах, и так далее. Я не могу думать ни о чем, кроме как о проблемах обработки изображений.

Phonon
источник
21

Гауссовы фильтры используются при обработке изображений, потому что они обладают свойством, что их поддержка во временной области равна их поддержке в частотной области. Это происходит из-за Гаусса, являющегося его собственным преобразованием Фурье.

Каковы последствия этого? Что ж, если поддержка фильтра одинакова в обоих доменах, это означает, что соотношение обеих опор равно 1. Как оказалось, это означает, что фильтры Гаусса имеют «произведение минимальной ширины полосы пропускания».

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

Но если изображения имеют края в качестве высокочастотных компонентов, традиционные LPF-функции также удаляют их, и визуально это проявляется в том, что края становятся более «смазанными».

Как тогда, чтобы удалить шум, но и сохранить высокочастотные края? Введите ядро ​​Гаусса. Поскольку преобразование Фурье гауссиана также является гауссовым, у фильтра Гаусса нет резкого обрезания на некоторой частоте полосы пропускания, за пределами которой удаляются все более высокие частоты. Вместо этого он имеет изящный и естественный хвост, который становится все ниже и ниже с увеличением частоты. Это означает, что он будет действовать как фильтр нижних частот, но также позволит использовать более высокочастотные компоненты, соразмерные с тем, насколько быстро его хвост распадается. (С другой стороны, у ФНЧ будет продукт с более высокой пропускной способностью, потому что его поддержка в F-области не так велика, как у гауссовцев).

Это тогда позволяет достигнуть лучшего из обоих миров - удаление шума, плюс сохранение края.

ошалевший
источник
2
Я не уверен, что вы можете сравнить две опоры напрямую, так как одна измеряется по времени / длине, а другая по Гц / радианам. Их морфология идентична, но свойство универсального масштабирования все еще сохраняется.
Фонон
Спасибо за напоминание о продукте с минимальной пропускной способностью. Однако, как отметил Фонон, уменьшение поддержки пространственной (~ временной) области обязательно увеличивает пропускную способность. Там нет никакого способа, как вы можете иметь как подавление шума, так и сохранить края с помощью простого фильтра Гаусса. Вот почему Perona & Malik разработали анизотропную фильтрацию.
nimrodm
@Phonon Как я видел, опоры - это просто количество ненулевых записей, описывающих функцию в любом домене - я верю, что они одинаковы. (Следовательно, соотношение 1). При этом произведение ширины полосы времени измеряется как произведение дисперсии функции по времени и частоте. Чем его нормализация отличается от авторов, я видел, что он приравнивается либо к 1/2, либо к 1/4.
Spacey
1
@nimrodm «уменьшение поддержки пространственной (~ временной) области обязательно увеличивает пропускную способность.» Да, это тенденция , проистекающая из обратной зависимости частоты от времени. (Вот откуда берется частотно-временная неопределенность). Однако функция Гаусса относится к классу, который полностью минимизирует это произведение. Учитывая обратную зависимость между временем и частотой, сделать это невозможно, если он не имеет равной поддержки в обеих областях.
Spacey
@nimrodm При анизотропной диффузии ядра, которые я видел, все еще гауссовы, хотя и с ковариационными матрицами, которые зависят от градиента изображения. (Это также нелинейный метод, VS гауссово сглаживание, которое является линейным). Но гауссиан остается используемым из-за его свойств.
Spacey
11

У вас уже есть хорошие ответы, но я просто добавлю еще одно полезное свойство 2D гауссовых фильтров, заключающееся в том, что они являются разделимыми , то есть 2D фильтр можно разложить на два 1D фильтра. Это может быть важным фактором производительности для больших размеров ядра, так как разделяемый фильтр MxN может быть реализован с M+Nмногократным сложением, тогда как неразделимый фильтр MxN требует M*Nмногократного добавления.

Пол Р
источник
2
Это хороший аргумент. 2D-фильтр Гаусса является радиально-симметричным и все же разделяемым, поэтому сложность реализации значительно снижается.
Nimrodm
1
В качестве ссылки, Руководство для ученых и инженеров по DSP дает превосходное описание этого свойства в главе 24 .
Калеб Рейстер
6

Руководство по imagemagick имеет отличное объяснение того, почему фильтрация с помощью функций sinc приводит к эффектам «звонка», а гауссианы - нет. ( http://www.imagemagick.org/Usage/fourier/#blurring и http://www.imagemagick.org/Usage/fourier/#circle_spectrum ). Если у вас есть края (несплошности) в вашем изображении (что делает большинство изображений), тогда полное отключение всех высоких частот оставляет вас в ряду в пространственной области. Вы также получаете сигнал, когда фильтруете прямоугольные волны с помощью функции sinc в одном измерении.

Блуждающая логика
источник
2

Там уже были прекрасные ответы, но я добавлю зерно соли, а не другую точку зрения:

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

Для изображения классическим априором является плавность значений (например, интенсивности) относительно положения (это можно рассматривать как функцию рассеяния точки, упомянутую @Phonon). Он часто моделируется как гауссовский, поскольку это форма, которую вы получите, смешивая различные объекты с известным радиусом гладкости (это называется центральной предельной теоремой ). Это в основном полезно, когда вы хотите создать производные изображения: вместо того, чтобы дифференцировать необработанный сигнал (который будет создавать шумный выход), вы должны делать это на сглаженном изображении. Это эквивалентно применению вейвлет-подобного оператора, такого как фильтры Габора .

meduz
источник