Опубликован: 10.10.2007 | Уровень: профессионал | Доступ: платный | ВУЗ: Московский государственный университет имени М.В.Ломоносова
Лекция 7:

Алгоритмы сжатия изображений

< Лекция 6 || Лекция 7: 12345 || Лекция 8 >

Алгоритм JPEG 2000

Алгоритм JPEG-2000 разработан той же группой экспертов в области фотографии, что и JPEG. Формирование JPEG как международного стандарта было закончено в 1992 году. В 1997 стало ясно, что необходим новый, более гибкий и мощный стандарт, который и был доработан к зиме 2000 года. Основные отличия алгоритма в JPEG 2000 от алгоритма в JPEG заключаются в следующем:

  1. Лучшее качество изображения при сильной степени сжатия. Или, что то же самое, большая степень сжатия при том же качестве для высоких степеней сжатия. Фактически это означает заметное уменьшение размеров графики "Web-качества", используемой большинством сайтов.
  2. Поддержка кодирования отдельных областей с лучшим качеством. Известно, что отдельные области изображения критичны для восприятия человеком (например, глаза на фотографии), в то время как качеством других можно пожертвовать (например, задний план). При "ручной" оптимизации увеличение степени сжатия проводится до тех пор, пока не будет потеряно качество в какой-то важной части изображения. Сейчас появляется возможность задать качество в критичных областях, сжав остальные области сильнее, т.е. мы получаем еще большую окончательную степень сжатия при субъективно равном качестве изображения.
  3. Основной алгоритм сжатия заменен на wavelet. Помимо указанного повышения степени сжатия это позволило избавиться от 8-пиксельной блочности, возникающей при повышении степени сжатия. Кроме того, плавное проявление изображения теперь изначально заложено в стандарт (Progressive JPEG, активно применяемый в Интернет, появился много позднее JPEG).
  4. Для повышения степени сжатия в алгоритме используется арифметическое сжатие. Изначально в стандарте JPEG также было заложено арифметическое сжатие, однако позднее оно было заменено менее эффективным сжатием по Хаффману, поскольку арифметическое сжатие было защищено патентами. Сейчас срок действия основного патента истек, и появилась возможность улучшить алгоритм.
  5. Поддержка сжатия без потерь. Помимо привычного сжатия с потерями новый JPEG теперь будет поддерживать и сжатие без потерь. Таким образом, становится возможным использование JPEG для сжатия медицинских изображений, в полиграфии, при сохранении текста под распознавание OCR системами и т.д.
  6. Поддержка сжатия однобитных (2-цветных) изображений. Для сохранения однобитных изображений (рисунки тушью, отсканированный текст и т.п.) ранее повсеместно рекомендовался формат GIF, поскольку сжатие с использованием ДКП весьма неэффективно к изображениям с резкими переходами цветов. В JPEG при сжатии 1-битная картинка приводилась к 8-битной, т.е. увеличивалась в 8 раз, после чего делалась попытка сжимать, нередко менее чем в 8 раз. Сейчас можно рекомендовать JPEG 2000 как универсальный алгоритм.
  7. На уровне формата поддерживается прозрачность. Плавно накладывать фон при создании WWW страниц теперь можно будет не только в GIF, но и в JPEG 2000. Кроме того, поддерживается не только 1 бит прозрачности (пиксел прозрачен/непрозрачен), а отдельный канал, что позволит задавать плавный переход от непрозрачного изображения к прозрачному фону.

Кроме того, на уровне формата поддерживаются включение в изображение информации о копирайте, поддержка устойчивости к битовым ошибкам при передаче и широковещании, можно запрашивать для декомпрессии или обработки внешние средства (plug-ins), можно включать в изображение его описание, информацию для поиска и т.д.

Идея алгоритма

Базовая схема JPEG-2000 (рис. 7.3) очень похожа на базовую схему JPEG. Отличия заключаются в следующем:

  1. Вместо дискретного косинусного преобразования ( DCT ) используется дискретное вэйвлет-преобразование ( DWT ).
  2. Вместо кодирования по Хаффману используется арифметическое сжатие.
  3. В алгоритм изначально заложено управление качеством областей изображения.
  4. Не используется явно дискретизация компонент U и V после преобразования цветовых пространств, поскольку при DWT можно достичь того же результата, но более аккуратно.
Конвейер операций, используемый в алгоритме JPEG-2000

Рис. 7.3. Конвейер операций, используемый в алгоритме JPEG-2000

Рассмотрим алгоритм по шагам.

Шаг 1.

В JPEG-2000 предусмотрен сдвиг яркости ( DC level shift ) каждой компоненты ( RGB ) изображения перед преобразованием в YUV. Это делается для выравнивания динамического диапазона (приближения к 0 гистограммы частот), что приводит к увеличению степени сжатия. Формулу преобразования можно записать как:

I'(x,y) = I(x,y) - 2^{ST - 1}

Значение степени ST для как каждой компоненты R, G и В свое (определяется при сжатии компрессором). При восстановлении изображения выполняется обратное преобразование:

I'(x,y) = I(x,y) + 2^{ST - 1}

Шаг 2.

Переводим изображение из цветового пространства RGB, с компонентами, отвечающими за красную ( Red ), зеленую ( Green ) и синюю ( Blue ) составляющие цвета точки, в цветовое пространство YUV. Этот шаг аналогичен JPEG (см. матрицы преобразования в описании JPEG), за тем исключением, что кроме преобразования с потерями предусмотрено также и преобразование без потерь. Его матрица выглядит так:

\left|
  \begin{array}{c}
    Y \\
    U \\
    V \\
  \end{array}
\right|
=
\left(
  \begin{array}{c}
    {\left\lfloor {{{{\rm{R}} + {\rm{2G}} + {\rm{B}}} \over {\rm{4}}}} \right\rfloor } \\ R-G \\
    B-G \\
  \end{array}
\right)

Обратное преобразование осуществляется с помощью обратной матрицы:

\left|
  \begin{array}{c}
    R \\
    G \\
    B \\
  \end{array}
\right|
=
\left(
  \begin{array}{c}
    U+G \\ 
    Y - {\left\lfloor {{{\rm{U}} + {\rm{V}}} \over {\rm{4}}}} \right\rfloor } \\
    V+G \\
  \end{array}
\right)

Шаг 3.

Дискретное wavelet преобразование ( DWT ) также может быть двух видов - для случая сжатия с потерями и для сжатия без потерь. Его коэффициенты задаются таблицами, приведенными ниже. Для сжатия с потерями коэффициенты представлены в табл. 7.1 и табл. 7.2

Таблица 7.1.
Коэффициенты при упаковке
i Низкочастотные коэффициенты {\rm{h}}_{\rm{L}} {\rm{(i)}} Высокочастотные коэффициенты {\rm{h}}_{\rm{H}} {\rm{(i)}}
0 1.115087052456994 0.6029490182363579
±1 0.5912717631142470 -0.2668641184428723
±2 -0.05754352622849957 -0.07822326652898785
±3 -0.09127176311424948 0.01686411844287495
±4 0 0.02674875741080976
Другие i 0 0

Таблица 7.2.
Коэффициенты при распаковке
i Низкочастотные коэффициенты {\rm{g}}_{\rm{L}} {\rm{(i)}} Высокочастотные коэффициенты {\rm{g}}_{\rm{H}} {\rm{(i)}}
0 0.6029490182363579 1.115087052456994
±1 -0.2668641184428723 0.5912717631142470
±2 -0.07822326652898785 -0.05754352622849957
±3 0.01686411844287495 -0.09127176311424948
±4 0.02674875741080976 0
Другие i 0 0
< Лекция 6 || Лекция 7: 12345 || Лекция 8 >
Сергей Иванов
Сергей Иванов
Россия, г. Смоленск ул. Николаева д. 19а кв. 56
Ирина Усанова
Ирина Усанова
Россия, г. Чистополь