УДК 621.397
МЕТОД СЖАТИЯ ИЗОБРАЖЕНИЙ, ОСНОВАННЫЙ НА РАЗЛОЖЕНИИ КВАЗИЦИКЛИЧЕСКИХ КОМПОНЕНТ ИЗОБРАЖЕНИЯ ПО СОБСТВЕННЫМ ВЕКТОРАМ СООТВЕТСТВУЮЩИХ СУБПОЛОСНЫХ МАТРИЦ
Е.Г. ЖИЛЯКОВ А.А. ЧЕРНОМОРЕЦ ВА. ГОЛОЩАПОВА
В работе изложен метод сжатия изображений, основанный на разложении квазициклических компонент изображения по собственным векторам соответствующих субпо-лосных матриц, который позволяет достичь высоких коэффициентов сжатия для отдельного класса изображений.
Белгородский государственный национальный исследовательский университет
Ключевые слова: изображение, частотный субинтервал, квазициклические компоненты, информационные классы, субинтервальные преобразования.
e-mail:
Zhilyakov@bsu.edu.ru Chernomorets@bsu.edu.ru VGoloschapova@ bsu.edu.ru
Современное развитие информационно-телекоммуникационных систем, основанное на использовании все более мощных компьютеров и наличии прогресса в технологии производства цифровых камер, сканеров и принтеров, приводит к широкому использованию изображений в цифровом виде. В связи с этим растет интерес к улучшению алгоритмов сжатия изображений. Сжатие данных актуально для достижения больших скоростей их передачи и эффективности хранения. Для решение проблемы сжатия изображения используются достижения во многих областях техники и математики. Одной из основных проблем является выделение и сокращение избыточности информации на изображении. В данной работе представлен новый метод сжатия изображений, основанный на разложении квазициклических компонент по собственным векторам соответствующих субполосных матриц, который позволяет достичь высоких коэффициентов сжатия для отдельного класса изображений.
Предлагаемый метод сжатия/восстановления изображений в цифровом виде, задаваемых в виде матрицы Ф = (_//к), / = 1,2,...,N1, к = 1,2,...,М2, значения элементов которой соответствуют яркости пикселей изображения, заключается в следующем.
На первом этапе необходимо определить доли энергии Р , г 1 = 1,2,..., Я1, г2=1,2,...,Я2, компонент изображения Ф = (//к), / = 1,2,...,М1, к = 1,2,...,М2, в раз-
личных частотных субинтервалах О [1],
гг (ЛпГ □ Ф □ □ Ф т
(1)
гг (ФФ )
где Ф - матрица исходного изображения, Аг =(а'|. ), /1,/2 = 1,..,М1, и АГ2 = (а^к2),
к1,к2 =1,..,М2 - субполосные матрицы, значения элементов которых вычисляются на основании следующих выражений,
н
М
□
□
г □■
«4 = □
2Свя
01 (2г! -1)(/! - /г)а1(/1 -/2)
2
п(*1 - /2)
□ , *1 *2,
□ п
□
□
г □■ «ккг = □
п(к1 - к 2 )
□ —, к1 = к2,
□ п 1 2
*1 Ф *2
2СоА ° 2 (2Г2 1)(к1 к 2) Яп °2(к. - к 2)
2 2
п
п
Я2
соответствующие частотному субинтервалу О ,
0 ^01(и □Ь ^,а гЛ V □ [рг2, Рг2+1 ] и и (и □[а Г1,а ^
и (и □ [
и □ -аг +1 ,-а
г+1 ’
□ V □
V □[-в Г2 + 1-Р г2] и
в г2+1,-в г 2] и и (и □ [-а Г1+1,-а Л V □ [вг2,в г2+1])} >
(2)
а г=(г1-1) ^ ’ в г2=(г2-1) ^ ,
г1 =1,2,...,Я, г =1,2,...,Я2 .
Далее на основе анализа значений указанных долей энергии определяется частотная концентрация Ст [2] для различных пороговых значений доли т суммарной энергии
1т
Ст =
Я1Я2
(3)
где 1т - минимальное количество частотных субинтервалов, в которых сосредоточена заданная доля т энергии изображения. Вычисленная частотная концентрация позволяет оценить долю суммарной площади частотных субинтервалов, в которых сосредоточена подавляющая доля энергии анализируемого изображения. Ранее [2] было выявлено, что для большинства изображений значение величины частотной концентрации указывает на сосредоточенность энергии изображений в узком частотном диапазоне.
Затем, на основе анализа значений долей энергии изображения, компоненты изображения, соответствующие различным субинтервалам, разбиваются на информационные классы V*, / = 1,2,... [3], информация о которых хранится в матрице М , содержащей номера информационных классов компонент, соответствующих частотному субинтервалу О , г1 =1,2,..., Я1, г2 =1,2,..., Я2.
2
и
Кф
Далее выполняются субинтервальные преобразования на основе их разложения [4] по собственным векторам субполосных матриц только в тех частотных субинтервалах О , которые соответствуют вычисленной ранее частотной концентрации
Ст для заданного значения доли т суммарной энергии
Ж г= а^Фа*, (4)
г = 1,2,..., Я , г2=1,2,..., Я2,
где 2 , аг2 - матрицы, столбцами которых являются собственные векторы матрицы ЛГ1 и Аг2
2п=(^, яЧ,., ^), аГ2=(Гг2, Я22,..., я X).
Как известно, основная потеря информации при реализации большинства методов сжатия изображений происходит на этапе квантования по уровням результатов различных преобразований над значениями яркости пикселей изображения. В данном методе используется метод квантования по уровням значений субинтер-вальных преобразований, который основан на выборе различного количества уровней квантования компонент изображения в зависимости от величины доли энергии данных компонент изображения.
Значения полученных субинтервальных преобразований Ж1, с учетом принадлежности М , г =1,2,..., Я1, г2 =1,2,...,Я2 некоторому информационному классу V*, / = 1,2,..., квантуются по соответствующему количеству уровней квантования , которые определяются как
N22 = 2м-2 , (5)
если М =0, то субинтервальные преобразования при данных г1 и г2 отбрасываются и не сохраняются в дальнейшем. В результате квантования получаем матрицы соответствующих квантованных значений.
Далее квантованные значения оставшихся субинтервальных преобразований Жц записываются в некоторый файл, запись осуществляется в порядке обхода элементов матриц Жц" по «змейке» («зигзаг-сканирование»). Данные в полученных
файлах сжимаются на основе кодирования по методу Хаффмана с фиксированной таблицей [5] или одним из других известных методов сжатия без потерь.
В результате получаем файл, который представляет собой информацию об изображении в сжатом виде, которой достаточно для его восстановления.
Для восстановления изображения на основе субинтервального синтеза последовательно выполняем следующие операции в обратном порядке: декодирование по методу Хаффмана, деквантование с учетом различного количества уровней квантования для различных компонент изображения. Затем выполняются обратные субин-тервальные преобразования
О:
~Л= 21/ цЖ 1гг ь 1г2 т, (6)
где Ь. 1, Ь - диагональные матрицы, составленные из собственных чисел матриц Л. 1
и Лг
г2
Ь п = diag(UГ1,Ц,...Ml),
Ь = diag(\ll2,\г2,...,\г),
собственные числа упорядочены по убыванию
^>^>...>4,
^>^>...>^5 .
На последнем этапе вычисленные для каждого частотного субинтервала результаты суммируются для получения восстановленного изображения
~ Я* Я\ ~
Ф = Ц ^ • (7)
г1=1 г2 =1
Данный метод в ходе вычислительных экспериментов продемонстрировал высокую эффективность сжатия изображений, энергия которых сосредоточена в незначительном количестве частотных интервалов.
Изложенный в статье метод сжатия изображений, основанный на разложении квазициклических компонент изображения по собственным векторам соответствующих субполосных матриц, является методом сжатия с потерями изображений на основе вариационного субинтервального частотного анализа/синтеза и может быть использован для минимизации затрат ресурсов информационнотелекоммуникационных систем на их хранение и передачу.
Работа выполнена при поддержке ФЦП «Научные и научно-педагогические кадры для инновационной России» на 2009-2013 годы, гос. контракт № 14.740.11.0390.
Список литературы
1. Жиляков Е.Г., Об эффективности метода оценивания значений долей энергии изображений на основе частотных представлений / Е.Г. Жиляков, А.А. Черноморец, А.Н. Зали-вин // Известия ОрелГТУ. Информационные системы и технологии. - № 2/52 (563). - 2009.
- С. 12-22.
2. Черноморец А.А., О частотной концентрации энергии изображений / А.А. Черноморец, В.А. Голощапова, И.В. Лысенко, Е.В. Болгова // Научные ведомости БелГУ. Сер. История. Политология. Экономика. Информатика. - 2011. - № 1 (96). - Вып. 17/1. - С. 146-151.
3. Черноморец А.А., Метод разбиения частотных субинтервалов на классы в задачах частотного анализа изображений / А.А. Черноморец // Информационные системы и технологии. - № 4(66). - 2011. - С. 31-38.
4. Жиляков Е.Г., Вариационные алгоритмы анализа и обработки изображений на основе частотных представлений / Е.Г. Жиляков, А.А. Черноморец // Белгород: Изд-во 000 «ГиК», 2009. - 146 с.
5. Миано Дж., Форматы и алгоритмы сжатия изображений в действии [Текст] / Дж. Миано. - М.: Издательство Триумф, 2003. - 336 с.
IMAGE COMPRESSION METHOD BASED ON THE EXPANSION OF QUASI-CYCLIC IMAGE COMPONENTS IN EIGENVECTORS SUBBAND MATRIX
E.G.ZHILYAKOV
A.A. CHERNOMORETSV.A. GOLOSCHAPOVA
BelgorodNational Research University
e-mail: Zhilyakov@bsu.edu.ru Chernomorets@bsu.edu.ru VGoloschapova@ bsu.edu.ru
In this article we present a image compression method based on the expansion of quasi-cyclic component of the image on the eigenvectors of the subband matrix, which achieves high compression ratios for a particular class of images
Key words: image, frequency subinterval, quasicyclic components, informational classes, subband conversion.
El м