Глава 5. АДАПТИВНЫЕ ОРТОГОНАЛЬНЫЕ ПРЕОБРАЗОВАНИЯ

Вейвлет-преобразование сигнала является сигнально-независимым. Октавополосное разбиение спектра, производимое им, подходит для большинства, но не для всех реальных сигналов. Желательно было бы иметь преобразование, адаптированное к сигналу, подобно ПКЛ, но имеющее быстрый алгоритм выполнения. Это эквивалентно тому, что преобразование было бы способно произвольно менять структуру разбиения частотно-временной плоскости в зависимости от сигнала. Каскадно соединенные блоки вейвлет-фильтров позволяют достичь этого. В данной главе будут рассмотрены адаптивные преобразования, которые на основе введенной функции стоимости реализуют произвольное разбиение частотно-временной плоскости сигнала.

В разделе 5.1 рассмотрены так называемые пакеты вейвлетов, или адаптация в частотной области. В разделе 5.2 рассмотрен так называемый алгоритм двойного дерева, или адаптация базиса разложения как в частотной, так и в пространственной областях. Дальнейшее развитие этих идей приведено в разделе 5.3. В разделе 5.4 обсуждаются вопросы размерности библиотеки базисов для всех преобразований и их вычислительная сложность. Раздел 5.5 содержит выводы по данной главе.

5.1. Пакеты вейвлетов (алгоритм одиночного дерева)

Итак, вейвлет-преобразование сигнала выполняется путем его пропускания через каскадно соединенные двухканальные схемы А-С (см. рис. 3.2). При этом каскадирование производится по низкочастотной области. Причина этого в неявном предположении, что эта область содержит больше информации об исходном сигнале. В результате получается «однобокое» дерево (рис. 5.1(а)). Данное предположение оправданно для многих реальных сигналов. В самом деле, оно означает, что наш сигнал является низкочастотным на большом интервале времени, а высокочастотные составляющие появляются на коротком интервале. Однако для некоторых сигналов это предположение не выполняется. Метод пакетов вейвлетов основан на определении того, по какой области на данном уровне выгоднее производить каскадирование. Для этого вначале производится каскадирование по обеим субполосам. В результате получается так называемое «полное», «сбалансированное» дерево (рис.5.1(б)), напоминающее дерево, присущее кратковременному преобразованию Фурье. Далее, на основе введенной функции стоимости определяется наилучший путь по этому дереву (рис. 5.1(в)). Если исходный блок вейвлет-фильтров

Рис. 5.1. Разбиение частотно-временной плоскости при помощи пакетов вейвлетов: (а) вейвлет-декомпозиция; (б) полная, аналогичная STFT декомпозиция; (в) пример декомпозиции при помощи пакетов вейвлетов

был ортогональным, то и схема, соответствующая любой конфигурации дерева, будет ортогональной, так как она есть не что иное, как каскадное соединение ортогональных блоков.

Таким образом, получается базис, адаптированный к сигналу. Отметим, что эта адаптация не требует обучения или знания статистических свойств сигнала. Вейвлет-преобразование (DWT), как и STFT, является частным случаем этого базиса. Адаптивность достигается за счет увеличения вычислительной стоимости. К счастью, разработан быстрый алгоритм поиска наилучшего базиса.

Пакеты вейвлетов были разработаны и исследованы Р.Койфманом и М.Викерхаузером. В качестве функции стоимости они использовали энтропию, понимаемую ими, как «концентрацию» числа коэффициентов M, требующихся для описания сигнала. Данная функция будет большой, если коэффициенты примерно одной величины, и малой, если все, кроме нескольких коэффициентов, близки к нулю. Таким образом, любое усреднение приводит к увеличению энтропии. Функция стоимости должна быть аддитивной. Это означает, что

Под энтропией в данном контексте понимается величина

где

Эта энтропия вычисляется для каждого узла полного дерева пакета вейвлетов. Далее сравнивается сумма энтропии двух потомков и энтропия их предка на дереве. Если энтропия предка оказалась меньше, отказываемся от его декомпозиции, то есть «обрезаем» дерево. Алгоритм рекурсивно продолжается до достижения вершины дерева. Доказано, что данный алгоритм приводит к наилучшему базису относительно

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

Рис. 5.2. Принятие решения в алгоритме одиночного дерева

Данный алгоритм получил название алгоритма одиночного (частотного) дерева. Он был применен для кодирования изображений. При этом на каждом этапе изображение делилось на четыре субполосы (структура, называемая квадродеревом вейвлет-коэффициентов).

При применении этого алгоритма для целей сжатия нельзя забывать о необходимости передавать декодеру информацию о структуре дерева. Одним из методов может быть посылка декодеру одного бита, указывающего, производилась или нет декомпозиция исходного изображения. Если - да, то посылаем еще четыре бита, указывающих решение по разбиению каждой из субполос. Легко показать, что для дерева максимальной глубины число дополнительных бит не превышает Скажем, при 4- уровневом разбиении изображения размером 512 х 512 потребуется 85 бит или примерно 0.000324 бит/пиксел, что совершенно незначительно.