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


30

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

Чому вони такі популярні в програмах для обробки зображень? Ці фільтри є результатом оптимізації будь-якого критерію чи є лише спеціальним рішенням, оскільки "пропускна здатність" зображення недостатньо визначена.

Відповіді:


28

Програми для обробки зображень відрізняються від прикладних програм для обробки аудіо, оскільки багато з них налаштовані на око. Маски Гаусса майже ідеально імітують оптичну розмитість (див. Також функції розподілу точок ). У будь-якому застосуванні для обробки зображень, орієнтованому на художнє виробництво, фільтри Гаусса використовуються для розмивання за замовчуванням.

Ще однією важливою кількісною властивістю фільтрів Гаусса є те, що вони скрізь є негативними . Це важливо, оскільки більшість 1D-сигналів змінюються приблизно 0 ( ) і можуть мати як позитивні, так і негативні значення. Зображення відрізняються тим сенсом, що всі значення зображення є негативними ( x R + ). Згортання з ядром Гаусса (фільтром) гарантує негативний результат, тому така функція відображає негативні значення на інші негативні значення ( f : R +R + ). Отже, результат завжди є іншим дійсним зображенням.xRxR+f:R+R+

Загалом, відхилення частоти в обробці зображень не настільки важливо, як в 1D-сигналах. Наприклад, у схемах модуляції ваші фільтри повинні бути дуже точними, щоб відхиляти інші канали, що передаються на різних несучих частотах тощо. Я нічого не можу вважати таким, що обмежує проблеми з обробкою зображень.


21

Гауссові фільтри використовуються при обробці зображень, оскільки вони мають властивість, що їх підтримка у часовій області дорівнює їх підтримці в частотній області. Це відбувається завдяки тому, що Гаусс був власною трансформацією Фур'є.

Які наслідки цього мають? Добре, якщо підтримка фільтра однакова в будь-якій області, це означає, що співвідношення обох опор дорівнює 1. Як виявляється, це означає, що фільтри Гаусса мають "продукт мінімальної пропускної здатності".

То що ви можете сказати? Ну а в обробці зображень одним дуже важливим завданням є видалення білого шуму, зберігаючи при цьому чіткі краї. Це може бути суперечливим завданням - білий шум існує на всіх частотах однаково, тоді як ребра існують у високочастотному діапазоні. (Раптові зміни просторових сигналів). У традиційному видаленні шуму за допомогою фільтрації сигнал фільтрується низькочастотним, що означає, що високочастотні компоненти вашого сигналу повністю видалені.

Але якщо зображення мають краї як високочастотні компоненти, традиційні LPF'ing також видалять їх, і візуально це проявляється тим, що краї стають все більш "розмитими".

Як тоді видалити шум, але також зберегти високі частоти ребер? Введіть ядро ​​Гаусса. Оскільки перетворення Фур'є у Гаусса також є гауссом, фільтр Гаусса не має різкого відсічення на деякій частоті пропускання смуги, за якою видаляються всі більш високі частоти. Натомість він має витончений і природний хвіст, який стає все нижчим у міру збільшення частоти. Це означає, що він буде функціонувати як фільтр низьких частот, але також допускати в компонентах більш високої частоти, пропорційних швидкості розпаду хвоста. (З іншого боку, LPF матиме більш високий часовий пропускний продукт, оскільки його підтримка в F-домені не настільки велика, як у гауссів).

Тоді це дозволяє досягти найкращого з обох світів - видалення шуму та збереження краю.


2
Я не впевнений, що ви можете порівняти дві опори безпосередньо, оскільки одна вимірюється за часом / довжиною, а інша - в Гц / радіанах. Їх морфологія ідентична, але властивість універсального масштабування все ще зберігається.
Phonon

Дякуємо, що нагадали про продукт із мінімальною пропускною здатністю в часі. Однак, як згадував Фонон, зменшення просторової (~ часової) підтримки домену обов'язково збільшує пропускну здатність. Немає способів придушити шум і зберегти краї за допомогою простого фільтра Гаусса. Ось чому Perona & Malik розробили анізотропну фільтрацію.
nimrodm

@Phonon Як я вже бачив, підтримка - це просто кількість ненульових записів, що описують функцію в будь-якому домені - я вважаю, що вони однакові. (Значить, співвідношення 1). Як говориться, продукт часової смуги вимірюється як добуток дисперсії функції у часі та частоті. Наскільки нормалізується його нормалізація від авторів, я бачив, що це дорівнює 1/2 або 1/4.
Спейсі

1
@nimrodm «зменшуючи просторову підтримку домену (~ часу) обов'язково збільшує пропускну здатність.», Так, то є тенденція , що випливають з зворотній залежності час-частота. (Звідси випливає невизначеність часу). Однак функція Гаусса є класом, який повністю мінімізує цей продукт. Зважаючи на зворотну залежність між часом і частотою, немає способу зробити це, якщо вона не має однакової підтримки в обох областях.
Спейсі

@nimrodm При анізотропній дифузії ядра, які я бачив, як і раніше є гауссовими, хоча і з матрицями коваріації, які залежать від градієнта зображення. (Це також нелінійний метод, згладжування В. С. гауссова, яке є лінійним). Але, гаусс залишається використаний завдяки своїм властивостям.
Спейсі

11

Ви вже маєте хороші відповіді, але я просто додам ще одну корисну властивість 2D-фільтрів Гаусса, тобто те, що вони є роздільними , тобто 2D-фільтр можна розкласти на два 1D-фільтри. Це може бути важливим врахуванням продуктивності для більших розмірів ядра, оскільки фільтр, що розділяється на MxN, може бути реалізований з M+Nмножинними додаваннями, тоді як нероздільний MxN-фільтр вимагає M*Nмноження-додавання.


2
Це хороший аргумент. 2D-фільтр Гаусса є радіально симетричним і все ще відокремленим, тому складність реалізації значно зменшується.
nimrodm

1
В якості довідника, Посібник для вченого та інженера щодо DSP надає чудовий опис цього властивості в главі 24 .
Caleb Reister

6

Посібник по Imagemagick чудово пояснює, чому фільтрація за допомогою функцій sinc призводить до ефекту "дзвінка", тоді як гаусианці цього не роблять. ( http://www.imagemagick.org/Usage/fourier/#blurring та http://www.imagemagick.org/Usage/fourier/#circle_spectrum ). Якщо у вас є зображення (розриви) у вашому зображенні (що робить більшість зображень), то повне відсікання всіх високих частот залишає вас з пульсаціями в просторовій області. Ви також отримуєте дзвінок, коли фільтруєте квадратні хвилі з функцією sinc в одному вимірі.


2

Відповіді вже було чудово, але я додам моє зерно солі, зовсім інше:

Фільтрування на найбільш абстрактному рівні може розглядатися як застосування деяких попередніх знань до деяких необроблених даних. Це означає, що застосовувати якийсь алгоритм фільтрації - це застосувати це до того, щоб знайти оптимальне співвідношення сигнал / шум, наприклад.

Для зображення класичним пріоритетом є плавність значень (наприклад, інтенсивності) щодо положення (це може розглядатися як функція поширення точки, згадана @Phonon). Його часто моделюють як гаусса, оскільки це форма, яку ви отримали при змішуванні різних об'єктів із відомим радіусом гладкості (це називається центральною граничною теоремою ). Це в основному корисно, коли ви хочете зробити похідні зображення: а не розмежовувати сирий сигнал (який би давав шумний вихід), ви повинні робити це на згладженому зображенні. Це еквівалентно застосувати оператор, що нагадує вейвлет, такий як фільтри Gabor .

Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.