1.3.3. Пирамида Лапласа

Один из первых методов для получения октавополосной декомпозиции был разработан и применен для кодирования изображения П.Буртом и Э.Адельсоном. Они использовали каскадно включенные гауссовские фильтры для получения избыточного представления сигнала, которое они назвали пирамидой Лапласа. Схема получения одного уровня пирамиды Лапласа (для одномерного сигнала) показана на рис. 1.5.

Сигнал пропускается через НЧ-фильтр и затем прореживается. В результате получается низкочастотная субполоса Высокочастотная субполоса формируется за счет последовательного выполнения следующих операций: интерполяции свертки с интерполирующим фильтром и вычитания результата из исходного сигнала. Реконструкция сигнала происходит путем интерполяции свертки с интерполирующим фильтром

Рис. 1.5. Схематическое изображение одного уровня пирамиды Лапласа

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

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

Для сравнения пирамиды Лапласа с другими субполосными преобразованиями представим ее как трехканальную систему А-С (см. рис.1.1), полученную путем деления на два сигнала: содержащего четные коэффициенты, и содержащего нечетные коэффициенты. Так как децимация во всех трех ветвях осуществляется в два раза, пирамидальное представление является избыточным в 3/2 раза. Фильтры системы А-С выражаются через фильтры пирамиды следующим образом:

Так как в данной системе не выполняется условие (1.18), пирамида Лапласа является неортогональным преобразованием. Пирамида Лапласа может строиться и для двумерного сигнала с использованием двумерного разделимого сглаживающего фильтра. Такая пирамида представляется как 5-полосная система А-С, в которой каждая субполоса прореживается в два раза как по горизонтали, так и по вертикали.

Пирамида Лапласа имеет несколько недостатков для кодирования изображений. Г лавным из них является то, что ошибки квантования высокочастотных коэффициентов распространяются на другие субполосы и проявляются в реконструированном изображении в виде шума. Как и в случае с преобразованием Габора, причина заключается в неортогональности преобразования. Другим недостатком является увеличение числа коэффициентов по сравнению с числом отсчетов исходного сигнала в 3/2 раза. И, наконец, двумерные базисные функции пирамиды Лапласа являются изотропными. Следовательно, они не подходят для устранения избыточности ориентированных структур, имеющихся в естественных изображениях.

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