CC BY
27
ПОВЫШЕНИЕ СТЕПЕНИ СЖАТИЯ И ВИЗУАЛЬНОГО КАЧЕСТВА ПРИ ИЕРАРХИЧЕСКОЙ КОМПРЕССИИ ИЗОБРАЖЕНИЙ ЗА СЧЕТ ПРЕДВАРИТЕЛЬНОЙ ФИЛЬТРАЦИИ
М.В. Гашников, Н.И. Глумов Институт систем обработки изображений Самарский государственный аэрокосмический университет имени академика С.П. Королева
Аннотация
В данной работе предложен алгоритм повышения степени сжатия и визуального качества при иерархической компрессии. Алгоритм основан на использовании предварительной фильтрации, приводящей изображение к виду, более «приспособленному» для иерархического сжатия. Проведено экспериментальное исследование разработанного алгоритма, показано преимущество схемы с предфильтрацией над базовым методом по коэффициенту сжатия и по визуальному качеству.
Введение
Иерархическая компрессия изображений [1-4] основана на использовании прореженных копий изображения для интерполяции промежуточных отсчетов и последующем энтропийном кодировании [5] ошибок интерполяции (постинтерполяционных остатков). В данной работе в качестве иерархического метода сжатия рассматривается метод иерархической сеточной интерполяции (ИСИ) [3-4]. Метод обладает чрезвычайно низкой вычислительной сложностью, высоким коэффициентом сжатия, позволяет контролировать погрешность и обладает многими другими достоинствами, характерными для иерархических методов компрессии.
Одним из недостатков метода ИСИ является недостаточная эффективность и снижение визуального качества декомпрессированных изображений при больших коэффициентах сжатия, обусловленное использованием равномерной шкалы квантования [5] для постинтерполяционных остатков. В данной работе предлагается способ уменьшения этих эффектов за счет предварительной фильтрации изображения.
Иерархическая компрессия изображений
Иерархическая компрессия основана на специальном представлении изображения. Пусть X = { х (т, п)} - исходное изображение. Представим
изображения X в виде объединения иерархических (масштабных) уровней X/:
Ь-1
х = их/,
/=0
ХЬ_! ={ хь_! (т,п)} ,
хі ={ хі (тп)} \ { хі+і (тп)} 1 є[°--£ _ 2] ,
где Ь - количество уровней,
а{ хі (т,п)} - массив отсчетов изображения, взятых -і
с шагом 2 по каждой координате.
При компрессии уровни обрабатываются последовательно, начиная со старшего уровня Xь _і. Старший уровень сжимается только статистическим кодером. Процедура компрессии любого из следую-
щих младших уровней X/, / е [0, Ь - 2] включает следующие этапы (см. также рис. 1):
С Исходное изображение^/
±
Интерполяция
\
Вычисление постинтерполяционных остатков
I
Квантование постинтерполяционных остатков
Статистическое
кодирование
Вычисление
восстановленных
значений
отсчетов
С Сжатое В
изображение •у
Восстановленное
изображение
О
Рис. 1. Иерархическая компрессия
1) Интерполяция
Производится интерполяция отсчетов { х/ (тп)}
уровня X/ на основе отсчетов { х% (т, п) , к > /}
уровней {, к > /}, уже прошедших компрессию и
восстановление. В результате вычисляются интерполирующие значения отсчетов уровня %1 (т, п). В
данной работе интерполяция осуществляется простым усреднением по ближайшим уже обработанным отсчетам более старшего уровня [4].
2) Вычисление постинтерполяционных остатков
Вычисляется разностный сигнал (массив «по-стинтерполяционных остатков»), т. е. массив разностей истинных и интерполирующих значений:
// (m, п) = X/ (т п)- X/ (т п) .
3) Квантование постинтерполяционных остатков
Выполняется квантование постинтерполяцион-ных остатков. В данной работе используется квантователь [5] с контролем максимальной погрешности етах
|// (т п)| +ета.
ді (т, п) = 8Ї%п (/і (т, п))
2р +1
іьшах Т1
где [..] обозначает выделение целой части.
|х(т,п) - х(т,п)| < етах, Предварительная
фильтрация
Этот квантователь гарантирует, что декомпрессированное изображение будет отличаться от исходного не более, чем на етах
(т, п)
4) Вычисление восстановленных значений отсчетов
По квантованным значениям постинтерполяци-онных остатков вычисляются восстановленные значения х/ (т, п) отсчетов уровня. Очевидно, что для
этого необходимо сначала вычислить восстановленные значения постинтерполяционных остатков:
// (т п) = д/ (т п) (тах + ^ а затем осуществить собственно восстановление:
х/ (т,п) = // (т,п) + X/ (т,п).
Эти восстановленные отсчеты будут нужны для интерполяции более младших иерархических уровней ^к , к < /} .
5) Статистическое кодирование
Осуществляется статистическое кодирование квантованных постинтерполяционных остатков. Поскольку распределение их вероятностей, как правило, является существенно неравномерным, в результате кодирования можно достичь значительного сокращения объема данных. На этом описание процедуры компрессии иерархического уровня закончено.
Таким образом, для каждого уровня формируется своя последовательность квантованных значений, которая отдельно сжимается статистическим кодером.
Предварительная фильтрация
Под предварительной фильтрацией в задаче компрессии будем понимать фильтрацию исходного изображения, которая производится до собственно компрессии (см. схему на рис. 2). Такая фильтрация, естественно, вносит погрешность, но, при правильном выборе фильтра, фильтрованное изображение может стать более «удобным» для данного метода сжатия, за счет чего схема обработки с предфильт-ром в результате может оказаться более эффективной (в координатах «коэффициент сжатия - погрешность восстановления»).
Как уже было отмечено выше, предфильтрация в данной работе применяется для борьбы со снижением качества декомпрессированных изображений при больших коэффициентах сжатия. Снижение качества в свою очередь обусловлено использованием равномерной шкалы квантования для постинтерполяцион-ных остатков.
Использование равномерной шкалы приводит к тому, что любая точка восстановленного изображения может изменить свою яркость на величину интервала квантования шкалы, а на «гладких» участках (участках с плавно изменяющейся яркостью) это приводит к появлению «выбитых» пикселов и небольших пятен (размером в несколько пикселов), яркость которых заметно отличается от яркости окружающих отсчетов.
С Исходное N
изображение /
Компрессия методом ИСИ
С Сжатое изображение >/
Рис. 2. Схема ИСИ-компрессии с предварительной фильтрацией
Такой эффект имеет место в тех отсчетах, в которых интерполяция производится неточно. В свою очередь интерполяция на гладких участках будет неточна вокруг пикселов, яркость которых сильно отличается от окружающих, то есть, прежде всего, на шумовых пикселах.
Такой вид погрешности метода ИСИ обуславливает требования к предфильтру: предфильтр должен подавлять малые шумы на участках плавно меняющейся яркости, при этом, по возможности оставляя нетронутой полезную информацию изображения, в первую очередь контурные линии. Такими свойствами обладает, в частности, сигма-фильтр [5]:
М N
X X х(т +1, п + ])и(^)1 (т +1, п + ])
,1 \ 1=-М I=^
х (т, п) =-------------------------------------
М
N
где
и(а) (к I)={1, если х(к,/)- х(m,п)<а
[0, иначе
а - параметр фильтра,
(2М +1) х (2 N +1) - размеры окна обработки,
а х’(т,п) - отфильтрованное изобра-
жение, которое подвергается собственно сжатию методом ИСИ.
Нетрудно видеть, что этот фильтр представляет собой усреднение по тем пикселам заданного окна обработки, значение яркости которых отличается от яркости центрального пиксела окна не более чем на заданную величину параметра алгоритма а .
Вычислительный эксперимент Для исследования эффективности схемы с пред-фильтрацией производились вычислительные эксперименты по ее сравнению с базовым методом ИСИ. В качестве меры эффективности при этом использовался известный показатель Р8Ж (реак-Ю-реак
nal-to-noise ratio) [5], который для байтовых изображений с диапазоном значений 0..255 имеет вид:
PSNR = 10log10
SmSn 2552
Sm-1Sn-1
z z (x(m, n )-x(m, n ))
m=0 n=0
где 8т х Бп - размеры исходного изображения, а х(т, п) и х(т, п) - отсчеты исходного и декомпрессированного изображений, соответственно.
Для схемы с предфильтрацией и базового метода ИСИ строилась зависимость Р8МЯ от коэффициента сжатия. Типичный вид этой зависимости, полученный для изображения «Лена» (см. рис. 3), показан на рис. 4. Так как в координатах «Р8МЯ-коэффициент сжатия» более эффективному методу соответствует более высоко расположенная кривая, то нетрудно видеть, что схема с предфильтрацией имеет выигрыш.
Рис. З. Тестовое изображение «Лена»
PSNR 38
37
36
35
34
■ Предфилътр+ИСИ Базовый ИСИ
10
12
14
16
18
Рис. 4. Зависимость от коэффициента сжатия К
для метода ИСИ с предфильтрацией и базового метода ИСИ
Производилось также сравнение визуального качества декомпрессированных изображений. На рис. 5 и рис. 6 показаны фрагменты восстановленных изображений для базового метода ИСИ и схемы с пред-
фильтром, соответственно. Изображения имеют практически одинаковую погрешность по сравнению с исходным изображением «Лена», причем изображению для схемы с предфильтром соответствует немного больший коэффициент компрессии.
При этом нетрудно видеть, что изображение для схемы с предфильтром лучше по визуальному качеству, так как на нем гораздо меньше выражены артефакты в виде «выбитых» пикселов и небольших пятен.
Поученные результаты позволяют сделать вывод о том, что использование сигма-предфильтрации является перспективным при иерархической компрессии изображений.
Рис. 5. Фрагмент декомпрессированного изображения для базового метода ИСИ
Рис. 6. Фрагмент декомпрессированного изображения для метода ИСИ с предфильтрацией
Заключение
В данной работе предложен алгоритм предварительной фильтрации для иерархической компрессии изображений на основе сигма-фильтра. Алгоритм предназначен для повышения коэффициента сжатия и улучшения визуального качества декомпрессиро-
ванных изображений. Проведено экспериментальное исследование разработанного алгоритма, показано преимущество схемы с предфильтрацией над базовым методом в координатах «PSNR-коэффициент сжатия» и по визуальному качеству.
Благодарности Работа выполнена при поддержке Министерства науки и образования РФ, Администрации Самарской области и Американского фонда гражданских исследований и развития (CRDF Project SA-014-02) в рамках российско-американской программы «Фундаментальные исследования и высшее образование» (BRHE), а также при поддержке гранта Президента РФ № НШ-1007.2003.01, и гранта РФФИ № 04-01-96507.
Литература
1. Александров В.В., Горский Н.Д. Представление и обработка изображений // Рекурсивный подход. Л.: Наука, 1985. - 192 с.
2. Kortman C.M. Redundancy Reduction - A Practical Method of Data Compression, Proc. IEEE, 55. (Marth, 67). No. 3. Р. 253-263.
3. Гашников М.В., Глумов Н.И., Сергеев В.В. Информационная технология компрессии изображений в системах оперативного дистанционного зондирования // Известия Самарского научного центра РАН, 1999. № 1. С. 99-107.
4. Gashnikov M.V., Glumov N.I., Sergeyev V.V. Control of Compression Ratios during Hierarchical Image Compression // Pattern Recognition and Image Analysis, 2005. Vol. 15. No. 1. Р. 170-171.
5. Методы компьютерной обработки изображений // Под редакцией В.А.Сойфера. М.: Физматлит, 2001. - 784 с.
