КОНТРОЛЬ СТЕПЕНИ СЖАТИЯ ПРИ ИЕРАРХИЧЕСКОЙ КОМПРЕССИИ ИЗОБРАЖЕНИЙ
Гашников М.В., Глумов Н.И., Сергеев В.В. Самарский государственный аэрокосмический университет
Аннотация
Предлагается новый метод обеспечения заданной степени сжатия при иерархической компрессии изображений. Основная идея метода заключается в использовании различных значений управляющего параметра (максимальной погрешности) для различных масштабных уровней представления изображения и коррекции значения погрешности на последнем уровне в зависимости от получаемого объема сжатых данных.
Введение
Основная идея иерархической компрессии [1-4,6] заключается в иерархическом прореживании сжимаемого изображения, интерполяции пропущенных отсчетов и последующем энтропийном кодировании [5] ошибок интерполяции (постинтерполяционных остатков). Контролируемым показателем качества при этом обычно является максимальная погрешность
етах = тах |х (т,п)- x (m,n)| ,
((m, n)
где x (m, n) и x (m, n) отсчеты исходного и декомпрессированного изображения соответственно. Очевидно, что невозможно одновременно обеспечить требуемые значения качества и степени сжатия, поскольку для различных изображений при одинаковом качестве достигаются различные степени сжатия.
Существует ряд важных прикладных задач, в которых требуется контролировать именно степень сжатия. Например, при компрессии изображений в системах дистанционного зондирования Земли требование равномерной загрузки канала связи с фиксированной пропускной способностью приводит к необходимости стабилизации скорости формирования выходного потока сжатых данных. В свою очередь, стабилизация выходного потока может быть достигнута только путем контроля степени сжатия при выборе управляющих параметров реализуемого алгоритма сжатия. Эти параметры должны задаваться для каждой очередной "порции" регистрируемого изображения, содержащей одну или несколько строк фиксированной длины (блок изображения).
Обычно для стабилизации скорости формирования выходного потока при компрессии используется буферная память, которая сглаживает колебания выходного объема данных. При этом значения управляющих параметров алгоритма сжатия выбираются исключительно исходя из заполненности буфера.
Однако в работах [3, 4] авторами было показано, что такой подход для иерархической компрессии малоэффективен, поскольку на нестационарных изображениях ведет к значительным колебаниям степени сжатия и заполненности буфера. Был предложен другой подход, основанный на использовании статистических характеристик сжимаемого изображения. Предложенный подход обеспечивал меньшую, в среднем, погрешность и более плавное изменение качества восстановленной видеоинфор-
мации, однако, также как и традиционный подход, требовал значительного объема буферной памяти, поскольку допускал довольно сильное изменение степени сжатия для различных блоков изображения.
В данной работе предлагается новый, более эффективный метод стабилизации скорости, практически не требующий буферной памяти за счет гораздо более точного обеспечения задаваемой степени сжатия для каждого блока изображения.
Описание метода стабилизации
При иерархической компрессии изображение X = {х (т,п)} представляется [2-4] в виде набора
непересекающихся иерархических уровней Х1: ь-1
X = и Х/,
/=0
Хь-1 = {хь (т п)}
Х/ = { (т,п)}\ {+1 (т,п)},0 < / < Ь - 2, где Х1 (т, п) - прореженная матрица отсчетов изображения, взятых с шагом 2/ по каждой координате:
XI (т, п ) = х (2} т,21 п) .
Данные иерархических уровней (квантованные постинтерполяционные остатки отсчетов уровней) сжимаются последовательно, начиная с наиболее прореженного.
В [3, 4] для каждого уровня предлагалось использовать одинаковую максимальную погрешность, определяемую в результате оценок статистических характеристик изображения. В данной работе предлагается использовать переменную по уровням погрешность. Сначала следует вычислить оценку необходимой погрешности для обеспечения требуемой степени сжатия. Однако в отличие от ранее использованного алгоритма в новом алгоритме используются другие, более информативные признаки. Далее на каждом уровне предлагается использовать погрешности, не превышающую вычисленную оценку. Однако такой способ не гарантирует достаточную точность достигаемой степени сжатия всего изображения. Следовательно, на последнем (наименее прореженном) необходимо уточнение значения максимальной погрешности. К моменту сжатия последнего уровня уже известны объемы сжатых данных, полученных для всех
остальных уровней. Соответственно, известен объем сжатых данных, который должен быть получен на последнем уровне для обеспечения заданной степени сжатия. Следовательно, максимальную погрешность на последнем уровне можно подбирать исходя из того, чтобы объем сжатых данных был равен требуемому.
Опишем подробнее алгоритм выбора погрешности. Сначала вычисляется признаки:
.-0—^, -,
-2 — ^, -3 —^,
где В - требуемая степень сжатия в бит/отсчет,
) - дисперсия постинтерполяционных остатков в отсчетах иерархического уровня X/ с координатами
2 N (д, е) (N (д, е)
(21 (2т +1),27 (2п +1)) ,
вида
) - дисперсия постинтерполяционных остатков в остальных отсчетах иерархического уровня Х[.
Оценка максимальной погрешности вычисляется как квадратичная форма:
г ] к г аг,у,к,г -0 -1 -2 -3
етах — X
0<г+у+к+г <2
где аг,у,к,г - коэффициенты, рассчитываемые путем параметрической настройки алгоритма сжатия для предполагаемого класса изображений.
На всех уровнях кроме младшего используются следующие значения максимальной погрешности:
А
етах — Ь1 етах, Ь -1 > I > 1,
где - коэффициенты, обеспечивающие наибольшую эффективность сжатия в координатах «качество-степень сжатия» и подбираемые экспериментально при настройке алгоритма для заданного класса изображений.
Для выбора погрешности на младшем уровне сначала вычисляется максимально допустимый объем данных (в битах), который можно сгенерировать, не нарушив ограничение на степень сжатия всего изображения:
V(0) — MNB - V(1), где MxN - размер изображения,
а V*-1-* - объем данных в битах, полученный при компрессии уровней с номерами от (Ь -1) до 0.
Затем для всех возможных значений максимальной погрешности вычисляется оценка степени сжатия, которая будет получена при использовании этого значения погрешности на последнем уровне. В качестве этой оценки используется энтропия [5] квантованных постинтерполяционных остатков на последнем уровне:
д — -2 ^
N
й
где NQ - общее количество квантованных постинтерполяционных остатков на последнем уровне; N (д, е) - количество квантованных постинтерполяционных остатков, равных д при квантовании с максимальной погрешностью, равной е; а 2 - максимальное значение отсчета изображения
г — 2Ь -1,
где Ь - количество бит для хранения каждого отсчета исходного (несжатого) изображения.
Собственно значение максимальной погрешности на последнем уровне выбирается исходя из того, чтобы эта оценка не превышала степень сжатия, допустимую на последнем уровне:
V (°) 1
е(0) — .
! Н (е): Н (е) < •
0<е<г | W V ' MN
Описанный способ выбора погрешности позволяет обеспечивать заданную степень сжатия с высокой точностью.
Вычислительный эксперимент
Для исследования эффективности разработанного метода контроля степени сжатия был проведен ряд вычислительных экспериментов. На рис. 1 показано 10-битное тестовое изображение, а на рис. 2 -то же самое изображение после компрессии и декомпрессии со степенью сжатия В — 1 бит/отсчет.
Рис. 1. Исходное изображение (10 бит/отсчет)
Рис. 2. Декомпрессированное изображение (1 бит/отсчет)
График на рис. 3 построен в координатах "степень сжатия - погрешность" и характеризует эффективность иерархического метода сжатия в целом. На
этом графике показана зависимость относительной квадратичной погрешности восстановления изображения Derror / Dimage от реально достигнутой степени сжатия Breai.
^error ^image 0,06
0,04
0,02
real
Рис. 3. Зависимость относительной квадратичной погрешности от фактической степени сжатия
Однако для исследования эффективности собственно разработанного метода контроля степени сжатия более информативным является другой параметр: отклонение фактически достигнутой степени сжатия Бгес/1 от требуемой степени сжатия В , заданной при компрессии: ДВ = Вгеа/ - В .
0 1 2 3 4 5 6 7 В
Рис. 4. Точность работы метода контроля степени сжатия
Зависимость этого параметра от заданной степени сжатия B показана на рис. 4. Как видно из этого рисунка, отклонения фактической степени сжатия от требуемой все-таки возможны, но чрезвычайно малы, что позволяет сделать вывод о высокой эффективности и перспективности использования разработанного метода контроля степени сжатия при иерархической компрессии изображений.
Благодарности Работа выполнена при поддержке Министерства образования РФ, Администрации Самарской области и Американского фонда гражданских исследований и развития (CRDF Project SA-014-02) в рамках российско-американской программы «Фундаментальные исследования и высшее образование» (BRHE); Российского фонда фундаментальных исследований (РФФИ), грант 04-01-96507.
Литература
1. Александров В. В., Горский Н. Д. Представление и обработка изображений // Рекурсивный подход. - Л.: Наука, 1985, - 192 с.
2. Методы компьютерной обработки изображений // Под редакцией В.А.Сойфера. М.: Физматлит, 2001. - 784 с.
3. Гашников М.В., Глумов Н.И., Сергеев В.В. Метод компрессии изображений в системах оперативного дистанционного зондирования // Доклады IX Всероссийской конференции «Математические методы распознавания образов» (ММРО-9-99) С.160-163.
4. Гашников М. В., Глумов Н. И., Сергеев В. В. Информационная технология компрессии изображений для систем реального времени // Труды VI международной конференции «Распознавание образов и анализ изображений» (РОАИ-6-2002), Великий Новгород, 2126 октября 2002. Т. 1. С.138-141.
5. Колесник В. Д., Полтырев Г. Ш. Курс теории информации ,М.: Наука, 1982. 416 с.
6. Kortman C. M. Redundancy Reduction - A Practical Method of Data Compression // IEEE, 1967. Vol. 55, No 3. Р. 253-263.