Банк фильтров в матричном определении преобразования Хаара
Матричное определение преобразования Хаара будет использовано для введения понятия банка фильтров. Преобразование Хаара можно интерпретировать как банк, состоящий из двух фильтров: один пропускает низкие частоты, а другой - высокие. Объяснение термину «фильтр», как простая идея банка фильтров ложится в основу концепции под диапазонных преобразований.
Конечно, преобразование Хаара является простейшим вейвлетным преобразованием, которое здесь употребляется для иллюстрации новых идей. Однако, его использование в качестве банка фильтров не очень эффективно. В конкретных приложениях используются более сложные множества фильтров, однако общая идея остается без изменений.
Фильтром называется линейный оператор, определяемый с помощью коэффициентов фильтра. Этот оператор применяется к входному вектору, в результате чего получается выходной вектор. Заметим, что пределы суммирования зависят от выбора последовательностей. Независимой переменной является время, поэтому удобно считать, что входная и выходная последовательности заданы при всех временах.
На практике входная последовательность всегда является конечной, поэтому считается, что бесконечный вектор имеет лишь конечное число ненулевых компонентов. Глубже вникнуть в работу фильтра можно с помощью простейшего входного сигнала. Эта последовательность равна нулю всюду кроме момента. Она называется единичным импульсом.
Несмотря на то, что в сумме, задающей свертку, не заданы пределы суммирования, легко видеть, что при любом имеется всего одно ненулевое слагаемое. Будем говорить, что выходной сигнал является откликом в момент времени на единичный импульс. Поскольку число коэффициентов фильтра конечно, то фильтр называется конечным импульсным откликом (FIR, finite impulse response).
Входной сигнал может быть одномерным (вектором из вещественных компонентов, как это предполагается в этом параграфе) или двумерным, то есть, изображением. Элементы подаются на вход фильтров один за другим, и каждый фильтр вычисляет и выдает на выход один сигнальный отклик. Число откликов в два раза больше числа входных сигналов (так как рассматривается два фильтра).
Разочаровывающий результат, поскольку мы хотим получить сжатие. Для исправления этого обстоятельства после прохождения через фильтр делается прореживание, при котором выбрасываются отклики с нечетными номерами. Эту операцию также называют децимацией.
После применения децимации число элементов на выходе равно числу элементов на входе. Пример: Легко построить банк фильтров, в котором низкочастотный фильтр вычисляет средние значения, а высокочастотный фильтр вычисляет полуразности, то есть, делается преобразование. Низкочастотный фильтр использует свертку для удаления из сигнала высокочастотной составляющей. Он пропускает низкие частоты. А высокочастотный фильтр делает все наоборот. Он удаляет из сигнала высокие частоты. Вместе они разделяют входной сигнал на поддиапазоны частот.
Вычисления четырех восстановленных значений. Затем эти восемь значений перемежаются для получения окончательной реконструкции исходной последовательности. Банки фильтров позволяют взглянуть на преобразование Хаара с общих позиций. Они же дают возможность построить другие, более изощренные вейвлетные преобразования. Использование банка фильтров имеет существенное преимущество по сравнению с одним фильтром, так как позволяет использовать свойства, недоступные одному фильтру. Прежде всего, это возможность восстанавливать исходный сигнал из выходных сигналов после их децимации (прореживания). Децимация не является трансляционно инвариантной операцией. После ее применения выходная последовательность состоит из четных членов, но если сделать задержку входа на один отсчет времени, то выходом будет служить последовательность, которая полностью отличается от подлинного выхода. Эти две последовательности являются двумя фазами выходного сигнала.
Выходы банка анализа называются коэффициентами поддиапазона. Их можно квантовать (если допустима частичная потеря информации), а затем сжимать с помощью RLE, методом Хаффмана, арифметическим кодированием или любым другим методом. В конечном счете, они подаются на вход банка синтеза, где к ним сначала добавляются нули (на место отброшенных членов), после чего они пропускаются через обратные фильтры, где из последовательностей формируется выходной вектор.
Децимация означает потерю данных. Обратная процедура не может компенсировать эту потерю, поскольку она лишь добавляет нули. Для того, чтобы достигнуть полного восстановления исходного сигнала, необходимо конструировать фильтры так, чтобы они могли компенсировать эти потери. Первое свойство, которое часто применяется при проектировании фильтров - это свойство ортогональности.
Банк анализа преобразования Хаара состоит из коэффициентов низкочастотного фильтра и из коэффициентов высокочастотного фильтра. Скалярное произведение этих двух векторов равно нулю. То есть, они ортогональны. Аналогично, банк синтеза состоит из двух ортогональных векторов коэффициентов фильтров.
Множество ортогональных фильтров размера. Эти фильтры ортогональны, поскольку скалярное произведение векторов их коэффициентов равно нулю: = 0. Примером служит фильтр Добеши (Daubechies) D4. Его коэффициенты приведены в уравнении. Если с помощью этого фильтра вручную посчитать несколько примеров, то будет видно, что реконструированный сигнал идентичен исходному входному сигналу, но отстает от него на три отсчета по времени.
Банк фильтров может быть биортогональным, менее ограничительным типом фильтров. Обратите внимание на схожесть фильтров. Нам уже известно, что выходы фильтра низких частот обычно пропускают через фильтры анализа несколько раз, при этом образуются все более короткие выходные сигналы.
Эту рекурсивную процедуру можно изобразить в виде дерева. Поскольку каждый узел этого дерева вдвое сокращает число выходных символов своего предшественника, то это дерево называется логарифмическим. Масштабирующая функция и вейвлетная функция получаются в пределе из логарифмического дерева.
Здесь обнаруживается связь дискретного вейвлетного преобразования (использующего банк фильтров) и континуального (непрерывного) вейвлетного преобразования. Если «взбираться» вверх по логарифмическому дереву с уровня на уровень, то одновременно вычисляется новое среднее с помощью новой масштабирующей функции, имеющей большую частоту, а также новые детали с помощью новой вейвлетной функции. Каждый уровень дерева соответствует удвоению частоты (или разрешения) по отношению к предыдущему уровню, поэтому логарифмическое дерево еще называют деревом с мультиразрешением.
Тот, кто профессионально работает со звуком и музыкой, знает, что два тона, имеющие частоты; звучат как один тон, но разной высоты. Частотный интервал между и делится на 12 подынтервалов (называемых хроматической гаммой), однако в западной традиции отдается предпочтение 8 из 12 тонов этого разделения (диатоническая гамма, которая состоит из 7 нот, а 8 нота называется «октавой»).
По этой причине основной частотный интервал, используемый в музыке, называется октавой. Следуя этой аналогии, можно сказать, что соседние уровни дерева мультиразрешения различаются на октаву частот. Основное преимущество под диапазонных преобразований состоит в том, что они выделяют различные частоты из входного сигнала, после чего можно точно контролировать часть сохраняемой (и отбрасываемой) информации из каждого поддиапазона частот.
На практике такое преобразование разлагает образ на несколько поддиапазонов, отвечающих разным частотным областям данного образа, после чего каждый поддиапазон можно квантовать независимо от остальных поддиапазонов
Главным недостатком таких преобразований является появление искусственных артефактов (то есть элементов, которых не было в исходном образе, таких как наложения, затухания, «звоны») в реконструированном образе.
По этой причине преобразование Хаара не является удовлетворительным, и основные исследования в этой области прежде всего направлены на поиски улучшенных фильтров.
Общая схема банка фильтров, включающая N частотных фильтров и 3 стадии процесса преобразования. Обратите внимание, что выход частотного фильтра Но каждой стадии направляется на вход следующей стадии для дальнейшего разделения, а комбинированный выход банка синтеза 3-ей и 2-ой стадий посылается наверх обратного фильтра банка синтеза предшествующей стадии.
- RSS
Наши услуги: