Иерархическая Компрессия в задаче хранения гиперспектральных изображений

АНАЛИЗ ГИПЕРСПЕКТРАЛЬНЫХ ДАННЫХ

ИЕРАРХИЧЕСКАЯ КОМПРЕССИЯ В ЗАДАЧЕ ХРАНЕНИЯ ГИПЕРСПЕКТРАЛЬНЫХ ИЗОБРАЖЕНИЙ

Гашников М.В., Глумов Н.И. Самарский государственный аэрокосмический университет имени академика С.П. Королёва (национальный исследовательский университет) (СГАУ), Институт систем обработки изображений РАН Аннотация

Исследуются возможности применения иерархической компрессии в задаче хранения гиперспектральных изображений. Приводятся результаты анализа изображений спектрометров Брее-Т1Я и АУМБ. Предлагаются алгоритмы аппроксимации спектральных каналов, позволяющие повысить эффективность компрессии при сохранении возможности доступа к отдельным компонентам. Приводятся результаты вычислительных экспериментов по исследованию эффективности разработанных алгоритмов на 16-битных гиперспектральных изображениях.

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

Введение

Основной спецификой гиперспектральных [1] изображений является большое количество каналов, которое часто достигает нескольких сотен. Например, результаты съёмки гиперспектрометрами SpecTIR [2] и AVIRIS [3] содержат 360 и 224 16-битных канала (компонент) соответственно. Следствием этого является нетривиальность задачи хранения таких изображений: чрезвычайно большой объём накопленных данных влечёт за собой непомерные требования к ёмкости запоминающих устройств и неприемлемо низкую скорость доступа к этим данным. Поэтому при решении задачи хранения гиперспектральных изображений невозможно обойтись без использования компрессии.

В работе [4] для компрессии гиперспектральных изображений было предложено использовать метод на основе иерархической сеточной [5-7] интерполяции (hierarchical grid interpolation, HGI), приведены результаты вычислительных экспериментов, демонстрирующих преимущества HGI над методом JPEG [8] в этой ситуации.

Однако предложенный в указанной работе вариант реализации метода HGI недостаточно эффективен для решения задачи хранения гиперспектральных изображений, т.к. не обеспечивает возможности доступа к отдельным компонентам без декомпрессии значительной части сжатого изображения. Кроме того, в указанной работе при проведении вычислительных экспериментов разрядность гиперспектральных данных с целью корректного сравнения с известными методами компрессии была редуцирована до 8 бит, в результате нельзя считать достаточными проведённые исследования по применимости и эффективности метода HGI для компрессии гиперспектральных данных.

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

тельных экспериментов по исследованию эффективности предложенных алгоритмов на реальных 16-битных гиперспектральных изображениях.

1. Показатели качества

В качестве меры погрешности при оценке эффективности методов компрессии в данной работе использовалась максимальная

и квадратичная 1

= max x

e2 =-

SVH_ погрешности [9]

^

где xs(m, n), xs (m, n) - отсчёты

компоненты номер 5 исходного и декомпрессированного изображений соответственно, 5 - количество спектральных компонент, Vх Н - пиксельный размер компонент изображений.

2. Анализ гиперспектральных изображений

Анализ особенностей гиперспектральных изображений производился на примере следующих находящихся в открытом доступе данных:

• пяти фрагментов изображений гиперспектрометра БрееТЖ [2] (от 356 до 360 16-битных каналов размером 320x600 пикселов);

• пяти фрагментов изображений гиперспектрометра АУ1ИБ [3] (224 16-битных канала размером 614x1086 пикселов).

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

Для анализа особенностей гиперспектральных изображений были вычислены оценки внутрикомпонентных и межкомпонентных характеристик этих изображений. Некоторые результаты показаны на рис. 3-5. Анализ этих результатов позволил сделать следующие выводы:

разница между максимум и минимумом градации яркости достигает тысяч и десятков тысяч. Такие изображения нельзя преобразовать в «байтовые»; компоненты очень сильно зависимы, т. к. межкомпонентная корреляция чрезвычайно высока (выше 0,95 для 85,2 % пар соседних компонент); большинство компонент имеют высокую внутри-компонентную корреляцию (выше 0,85 для 87,4 % компонент). Можно ожидать, что использование компрессии даст заметный эффект.

г V v т

'¿Г*

¿тщ

'V •

'Ч 1' Л ^ "

Рис. 1. Компоненты №27, №269 изображения "Urban and Mixed Environment" спектрометра SpecTIR

Рис. 2. Компоненты №26, №102 изображения "Сиргие-1 "спектрометра ЛУ1ШБ

3. Требования к методу компрессии гиперспектральных изображений

Сделанные в предыдущем разделе выводы, в свою очередь, порождают следующие требования к методу сжатия гиперспектральных изображений:

1) возможность компрессии многокомпонентных изображений;

2) быстрый доступ к фрагментам сжатых изображений в различных масштабах;

3) низкая вычислительная сложность декомпрессии;

4) строгий контроль погрешности;

5)

6) 7)

высокая эффективность в режиме «без погрешности» и при малых погрешностях; использование межкомпонентных зависимостей; быстрый доступ к заданным компонентам сжатых изображений в различных масштабах; компресия 16-битных изображений.

АХ 30000

20000 10000

0 АХ

9000 6000 3000

рг\

т

100 200 300 s

Г 1

1 V

О 50 100 150 200 s

Рис. 3. Внутрикомпонентная разница DX между максимумом и минимумом в зависимости от номера компоненты s для изображений "Urban and Mixed Environment"

E, xlO2

200 s

0 50 100 150 200 s

Рис. 4. Оценки внутрикомпонентных среднего E, дисперсии D и коэффициента корреляции р в зависимости от номера компоненты s для изображения «LowAltitude» спектрометра AVIRIS

4. Компрессия на основе HGI

В работе [4] предлагается для компрессии гиперспектральных изображений использовать метод HGI [4-7]. Метод основан на использовании «сильно» прореженных наборов отсчётов изображения для интерполяции «сла-

бо» прореженных наборов отсчётов с последующим кодированием гостинтерполяционных остатков.

Метод удовлетворяет требованиям 1 - 5 предыдущего раздела, в частности, обеспечивает [5] контроль максимальной погрешности и независимость скорости доступа к фрагменту изображения от масштаба за счёт блочно-иерархического [6 - 7] формата сжатых данных. Этот формат подразумевает разбиение изображения на независимо обрабатываемые блоки при компрессии. В архиве сжатые данные блоков и сжатые данные масштабных уровней внутри блока лежат отдельно. В результате этого для декомпрессии фрагмента изображения достаточно декомпрессировать только необходимые блоки в необходимом масштабе, что позволяет добиться независимости времени доступа к сжатым данным от масштаба.

h"— i

1

0 50 100 150 200 s

Рис. 5. Оценка коэффициента корреляции р1 между соседними компонентами (компонентами номер s и номер s+1) в зависимости от номера компоненты s для изображения «LowAltitude» спектрометра AVIRIS

5. Компрессия на основе «скользящей аппроксимации компонент»

Для использования межкомпонентных зависимостей (требование 6 предыдущего раздела) в работе [3] предложено использовать зависимость между спектральными компонентами за счёт предсказания компонент. Собственно предсказание предлагается осуществлять как аппроксимацию компоненты на основе уже прошедших компрессию и восстановление других компонент. Высокая межкомпонентная корреляция должна обеспечить хорошую точность предсказания, а сжатие разности между исходной и предсказанной спектральной компонентами вместо сжатия исходной спектральной компоненты должно существенно повысить коэффициент сжатия. Опишем алгоритм формально.

Пусть { Xs ,0 < s < S} - S-компонентное гиперспектральное изображение, состоящее из двумерных спектральных компонент - однокомпонентных изображений Xs. Указанные компоненты будем сжимать последовательно, от меньших номеров к большим. При сжатии каждой компоненты Xs сначала будем вычислять её предсказанное (аппроксимирующее) значение:

N-1 _

Xs = £ k1Xs' ,0 < s < S ,

i=0

где X , 1 > 0 - уже прошедшие компрессию и восстановление предыдущие компоненты, на основе которых производится аппроксимация, X , 1 < 0 - заполненные нулями матрицы, введённые для упрощения запи-

си формул, N - количество предыдущих восстановленных спектральных компонент, используемых для аппроксимации (параметр алгоритма), {к1,0 < 1 < Щ -коэффициенты аппроксимации, которые, исходя из метода наименьших квадратов, являются решением системы линейных уравнений: Як = В ,

где к = {к1,0 < 1 < N} - искомый вектор коэффициентов аппроксимации, Я = ,0 < 1, ] < Щ - матрица коэффициентов корреляции между декомпрес-

— 3-1-1 "Т^-1

сированными компонентами X их ,

В ={Б1,0 < 1 < Щ - вектор коэффициентов корреляции между текущей предсказываемой компонентой X11 и декомпрессированной компонентой X .

Собственно компрессии методом ЫШ (с заданной максимальной погрешностью) будем подвергать не компоненту X11, а разность между X33 и предсказанной

(аппроксимирующей) компонентой X1 . После описанной процедуры компрессии компоненты X11 необходимо сразу же осуществить её декомпрессию. Это необходимо для обратной связи: восстановленная

компонента X будет использована при аппроксимации последующих спектральных компонент.

Таким образом, набор опорных компонент для каждой текущей предсказываемой компоненты представляет собой расположенное в спектральной плоскости «скользящее окно», поэтому описанный алгоритм компрессии далее будем называть «алгоритм на базе скользящей аппроксимации компонент».

Описанный подход предназначен для увеличения коэффициента сжатия за счёт использования межкомпонентных зависимостей при сохранении контроля максимальной погрешности.

6. Компрессия на основе «независимых порций компонент»

Предложенный в [4] алгоритм компрессии на основе «скользящей аппроксимации компонент» (описание см. выше) непригоден для использования в задаче хранения гиперспектральных изображений. Причина заключается в последовательном сжатии спектральных компонент, при котором «предыдущие» спектральные компоненты используются для аппроксимации «последующих». В результате для декомпрессии произвольной спектральной компоненты придётся декомпрессировать все предыдущие компоненты (их могут быть сотни), что является серьёзным недостатком при организации быстрого доступа к хранимым гиперспектральным изображениям (не выполняется требование 7 раздела 3).

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

пользовать подход, основанный на аппроксимации компонент внутри «независимых порций компонент».

При сжатии набор спектральных компонент разделяется на независимые «спектральные порции» (по N компонент в каждой, рис. 6), а внутри каждой порции используется упомянутый алгоритм «скользящей» аппроксимации спектральных компонент. Тогда для декомпрессии произвольной компоненты необходимо декомпрессировать не все предыдущие компоненты изображения, а только предыдущие компоненты соответствующей спектральной порции.

7. Компрессия на основе «общих опорных компонент»

Коэффициент сжатия алгоритма на основе «независимых порций компонент» неизбежно будет меньше, чем коэффициент сжатия алгоритма на основе «скользящей аппроксимации компонент». Причиной этого является то, что для аппроксимации каждой спектральной компоненты в алгоритме «независимых порций» будет использоваться гораздо меньшее (в среднем) число базовых компонент. Особенно велик проигрыш будет при компрессии первой компоненты каждой порции, т.к. для неё вообще нет опорных компонент, и она не будет аппроксимирована. Для уменьшения описанных эффектов в данной работе предлагается алгоритм аппроксимации на основе «общих опорных компонент» (рис. 7).

Из изображения с некоторым шагом N выбираются «общие опорные» спектральные компоненты. Также задаётся целочисленный параметр С. Общие опорные компоненты сжимаются алгоритмом на основе «скользящей аппроксимации компонент». При этом для аппроксимации каждой из «общих опорных компонент» используются ^ + С - 2) предыдущие опорные компоненты, уже прошедшие компрессию.

Затем, после сжатия «общих опорных компонент», независимо друг от друга сжимаются порции компонент, расположенных между «общими опорными». Для сжатия этих порций используется аппроксимационный алгоритм компрессии, аналогичный алгоритму на основе «скользящей аппроксимации компонент». Компоненты каждой порции сжимаются последовательно, с использованием аппроксимации на основе компонент, уже прошедших компрессию. Этими аппроксимирующими компонентами являются предыдущие компоненты той же самой порции, дополненные ближайшими «общими опорными компонентами» таким образом, чтобы общее количество аппроксимирующих компонент также равнялось (Ы+ С- 2). Таким образом, параметр С равен минимальному количеству «общих опорных компонент», используемых для аппроксимации «остальных» компонент.

Можно ожидать, что алгоритм на основе «общих опорных компонент» будет иметь коэффициент сжатия лучше, чем алгоритм на основе «независимых порций компонент». «Расплатой» за это будет являться некоторое уменьшение скорости, т. к. для декомпрессии произвольной компоненты нужно декомпрессировать не только предыдущие компоненты той же самой порции, но и необходимые «общие опорные» компоненты. Не-

гативный эффект от такого замедления будет тем меньше, чем больше компонент изображения нужно декомпрессировать, т. к. «общие опорные» компоненты распаковываются однократно при декомпрессии первой необходимой компоненты, а затем используются при декомпрессии всех последующих компонент.

Спектральная порция №1 Спектральная порция №2

__. ._, ± 0 1 2 3 V г( ♦ 4 + 5 + 6 + 7 V г(

Гиперспектральное изображение

Рис. 6. Алгоритм аппроксимации компонент на основе «независимых порций компонент» N = 4)

Набор «общих опорных» спектральных компонент

1

О

Спектральная порция №1

Спектральная порция №2

Гиперспектральное изображение

Рис. 7. Алгоритм аппроксимации компонент на основе «общих опорных» компонент N = 3, С = 1)

8. Реализация алгоритмов компрессии

Проведена реализация следующих алгоритмов ИИ-компрессии гиперспектральных изображений: алгоритма с независимой обработкой спектральных компонент, алгоритма со «скользящей аппроксимацией» компонент, алгоритма с аппроксимацией на основе «независимых порций компонент», алгоритма с аппроксимацией на основе «общих опорных компонент».

Общие параметры алгоритмов: допустимая максимальная погрешность £шж; размер Уь хНь независимо обрабатываемых блоков; количество Ь масштабных уровней [4]; количество N спектральных компонент, используемых для аппроксимации текущей компоненты. Кроме того, алгоритм «общих опорных» компонент имеет ещё параметр С - количество опорных компонент, которые используются для аппроксимации «неопорных» компонент (смотри описание алгоритма выше).

Все реализованные алгоритмы работают с 16-битными изображениями (требование 8 раздела 3).

9. Исследование влияния параметров

на эффективность компрессии

Для того, чтобы сформулировать рекомендации по выбору наилучших значений для параметров алгоритмов, были проведены вычислительные эксперименты на описанных 16-битных изображениях. Некоторые результаты для изображений спектрометра ЛУ1ИБ показаны на рис. 8 - 11.

К

3,3

2,9

2,5

——Е-

Общие опорные компоненты ....(&).... Скользящая

аппроксимация

__гт__ Независимые порции

__ж..— Независимые

компоненты

2,1

12 34 5 6789 10 N

Рис. 8. Зависимость коэффициента сжатия К от количества N опорных компонент при аппроксимации для изображения «Сирг^е-2» Ь = 5, С = 4, УьхНь = 112x614, ешах = 0)

О 2 4 б 8 10 С

Рис. 9. Зависимость коэффициента сжатия К от минимального количества С «общих опорных компонент», используемых для аппроксимации «остальных» компонент для изображения «СиргИе-1» Ь = 5, N= 7, УьхИь = 112x614, ешах= 0)

Vb=224

..... Vb=112

-д— Vb=64 -Sfc— Vb=224 #— Vb=32

120 220 320 420 520 620 Hb

Рис. 10. Зависимость коэффициента сжатия K от горизонтального размера блока Hb

К 3,51

3,48

3,45

3,42

1 2 3 4 5 6 7 8 9 L

Рис. 11. Зависимость коэффициента сжатия K от количества иерархических уровней L

для изображения «Cuprite-1» N= 7, С = 4, VbxHb = 128x128, emax = 0)

-Д- Общие опорные компоненты -<§>- Скользящая аппроксимация

р.....<§>.....в.....<§>.....<§>.....<§>.....•

Полученные результаты позволяют дать следующие рекомендации по выбору параметров.

• Результаты улучшаются при увеличении размера блоков. Скорее всего, причина заключается в уменьшении доли служебной информации блоков. Для блоков размером больше, чем 100x100, отставание от наилучших значений является незначительным. Можно выбирать любые размеры блока больше, чем 100x100, т.к. для них отставание от наилучших значений является незначительным.

• Наилучшие значения количества уровней: 4, 5, 6. Однако почти без потери эффективности для количества иерархических уровней можно брать любое значение больше трёх. При этом, исходя из соображений скорости доступа к сжатым данным в различных масштабах, количество уровней рекомендуется рассчитывать исходя из размера сжимаемых изображений: Ь = 1 + [^2шт(У,Н)], где [..] - знак выделения целой части.

• Количество опорных компонент N слабо влияет на коэффициент сжатия алгоритмов скользящей аппроксимации и общих опорных компонент, поэтому следует исходить из увеличения скорости доступа и использовать значения больше 9.

• Уменьшение количества опорных компонент N для алгоритма «независимых порций» заметно ухудшает его коэффициент сжатия. При использовании этого алгоритма необходимо выбирать количество компонент, обеспечивающее компромисс между коэффициентом сжатия и скоростью доступа. Предположительно, наилучшими значениями будут значения N между 7 и 10, т.к. после их превышения рост эффективности замедляется.

• Параметр С слабо влияет на коэффициент компрессии разработанных аппроксимационных алгоритмов. Его значение можно выбирать от 1 до 4, т. к. после превышения этих значений происходит замедление роста эффективности.

10. Оценка эффективности НС1-компрессии при хранении изображений

Для оценки эффективности компрессии были построены зависимости коэффициента компрессии от максимальной и квадратичной погрешностей для всех разработанных алгоритмов (рис. 12) на реальных 16-битных изображениях гиперспектрометров БрееИЯ и ЛУ1ШБ.

Поученные результаты позволяют сделать следующие выводы:

а) все алгоритмы показывают достаточно высокий коэффициент сжатия, чтобы рекомендовать их использование в системах хранения гиперспектральных изображений;

б) использование любой аппроксимации компонент позволяет существенно улучшить коэффициент сжатия; аппроксимационные алгоритмы можно рекомендовать к использованию;

в) алгоритмы в порядке убывания коэффициента сжатия: алгоритм общих опорных компонент, алгоритм скользящей аппроксимации, алгоритм не-

г)

д)

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

выигрыш от использования аппроксимации компонент растёт с увеличением коэффициента сжатия; при нулевой и малой погрешностях алгоритмы «скользящей аппроксимации» и «общих опорных компонент» показывают лучшие результаты. Так как эти результаты примерно одинаковы и при этом алгоритм «общих опорных компонент» позволяет уменьшить время доступа, то он является более предпочтительным из двух указанных алгоритмов (при малых погрешностях).

К (-Общие опорные компоненты^

♦...... Скользящая аппроксимация

8 и---Н5- Независимые порции

-ф----Независимые

—СЬ Общие опорные

компоненты ♦— Скользящая аппроксимация

0 2 4 6 8

Рис. 12. Усреднённый по пяти изображениям БресТШ коэффициент сжатия К в зависимости от максимальной погрешности етах и от квадратичной погрешности ^

Заключение

Разработаны новые алгоритмы аппроксимации спектральных компонент, хорошо приспособленные для использования ИИ-компрессии при решении задачи хранения гиперспектральных изображений. Проведена оценка эффективности, а также сравнение между собой указанных и ранее разработанных алгоритмов ИЯ-компрессии на реальных 16-битных изображениях гиперспектрометров. Показана перспективность использования ИИ-компрессии для решения задачи хранения гиперспектральных изображений.

Благодарности

Работа выполнена при государственной поддержке:

• Министерства образования и науки РФ в рамках реализации мероприятий Программы повышения конкурентоспособности СГАУ среди ведущих мировых научно-образовательных центров на 2013-2020 годы;

• грантов РФФИ (проекты 13-01-12080-офи_м, 13-07-97006-р_поволжье_а, 12-07-00751 -а);

• Министерства образования и науки Российской Федерации (в рамках постановления Правительства Российской Федерации от 09.04.2010 г. № 218: договор № 02.Г36.31.0001 от 12.02.2013).

Литература

1. Borengasser, M. Hyperspectral Remote Sensing - Principles and Applications / M. Borengasser [et al]. - CRC Press, 2004. - 128 p.

2. SpecTIR Data - Advanced Hyperspectral and Geospatial Solutions / Corporate Headquarters SpecTIR Remote Sensing Division // http://www.spectir.com/free-data-samples.

3. AVIRIS Data - Ordering Free AVIRIS Standard Data Products / Jet Propulsion Laboratory // http://aviris.jpl.na-sa. gov/data/ free_data. html.

4. Гашников, М.В. Иерархическая сеточная интерполяция при сжатии гиперспектральных изображений / М.В. Гашников, Н.И. Глумов // Компьютерная оптика. -2014. - Т. 38, № 1. - С. 87-93.

5. Сергеев, В. В. Метод иерархической компрессии космических изображений / М.В. Гашников, Н.И. Глумов, В.В. Сергеев // Автоматика и телемеханика. - 2010. -№ 3. - C. 147-161.

6. Gashnikov, M.V. Preparing a Common Raster Coverage for a Territory Based on Hierarchical Compressed Presentation of Orthoimages / M.V. Gashnikov, N.I. Glumov, A.V. Chernov // Pattern Recognition and Image Analysis. Advances in Mathematical Theory and Applications. - 2009. - V. 19, Issue 1. - P. 39-42.

7. Gashnikov, M.V. Hierarchical compression of the multidi-mansional data in the regional bank of samara region satellite images / M.V. Gashnikov, N.I. Glumov // Proceedings of 9-th International Conference on Pattern Recognition and Image Analysis: New Information Technologies (PRIA-9-2008). -Russian Federation, Nizhni Novgorod, September 15-19, 2008. - V. 1. - P. 159-161.

8. Wallace, G. The JPEG Still Picture Compression Standard / G. Wallace // Communications of the ACM. - 1991. -V. 34, Issue 4. - P. 30-44.

9. Методы компьютерной обработки изображений / М.В. Гашников, Н.И. Глумов, Н.Ю. Ильясова, В.В. Мясников, С.Б. Попов, В.В. Сергеев, В.А. Сойфер, А.Г. Храмов, А.В. Чернов, В.М. Чернов, М.А. Чичёва, В.А. Фурсов; под ред. В.А. Сойфера. - 2-е изд., испр. - М.: Физ-матлит, 2003. - 784 с.

References

1. Borengasser, M. Hyperspectral Remote Sensing - Principles and Applications / M. Borengasser [et al]. - CRC Press, 2004. - 128 p.

2. SpecTIR Data - Advanced Hyperspectral and Geospatial Solutions / Corporate Headquarters SpecTIR Remote Sensing Division // http://www.spectir.com/free-data-samples.

3. AVIRIS Data - Ordering Free AVIRIS Standard Data Products / Jet Propulsion Laboratory // http://aviris.jpl.na-sa. gov/data/ free_data. html.

4. Gashnikov, M.V. Hierarchical grid interpolation for hyperspectral image compression / M.V. Gashnikov, N.I. Glumov // Computer Optics. - 2014. - V. 38(1). - P. 87-93.

5. Sergeev, V.V. A hierarchical compression method for space images / M.V. Gashnikov, N.I. Glumov, V.V. Sergeev // Automation and Remote Control. - 2010. - V. 71(3). - P. 501-513.

6. Gashnikov, M.V. Preparing a Common Raster Coverage for a Territory Based on Hierarchical Compressed Presentation of Orthoimages / M.V. Gashnikov, N.I. Glumov, A.V. Chernov // Pattern Recognition and Image Analysis.

Advances in Mathematical Theory and Applications. -2009. - V. 19, Issue 1. - P. 39-42.

7. Gashnikov, M.V. Hierarchical compression of the multidi-mansional data in the regional bank of samara region satellite images / M.V. Gashnikov, N.I. Glumov // Proceedings of 9-th International Conference on Pattern Recognition and Image Analysis: New Information Technologies (PRIA-9-2008). - Russian Federation, Nizhni Novgorod, September 15-19, 2008. - V. 1. - P. 159-161.

8. Wallace, G. The JPEG Still Picture Compression Standard / G. Wallace // Communications of the ACM. - 1991. -V. 34, Issue 4. - P. 30-44.

9. Computer Image Processing, Part II: Methods and algorithms / A.V. Chernov, V.M. Chernov, M.A. Chicheva, V.A. Fursov, M.V. Gashnikov, N.I. Glumov, N.Yu. Ilyasova, A.G. Khra-mov, A.O. Korepanov, A.V. Kupriyanov, E.V. Myasnikov, V.V. Myasnikov, S.B. Popov, V.V. Sergeyev, V.A. Soifer. -VDM Verlag, 2009. - 584 p.

HIERARCHICAL COMPRESSION FOR HYPERSPECTRAL IMAGE STORAGE

M. V. Gashnikov, N.I. Glumov Samara State Aerospace University, Image Processing Systems Institute, Russian Academy of Sciences

Abstract

We investigate the possibility of using of hierarchical storage compression for problem of hy-perspectral images storage. Results of image analysis of SpecTIR and AVIRIS spectrometers shown. Approximation of the spectral channels is proposed to improve the efficiency of the method while still allowing access to individual components. Computational experiments to study the efficiency of the developed algorithms on 16-bit hyperspectral images implemented.

Key words: image compression, hyperspectral images storage, control of maximum deviation, hierarchical grid interpolation.

Сведения об авторах Гашников Михаил Валерьевич, 1975 года рождения, в 1998 году окончил Самарский государственный аэрокосмический университет имени академика С.П. Королёва (СГАУ), кандидат наук, доцент, работает старшим научным сотрудником лаборатории математических методов обработки изображений Института систем обработки изображений РАН (ИСОИ РАН). Область научных интересов: сжатие изображений, обработка космических снимков, геоинформационные технологии. E-mail: mgash @smr.ru .

Mikhael Valeryevich Gashnikov (b. 1975) graduated from S.P. Korolyov Samara State Aerospace University (SSAU), PhD and associate professor. He is the senior researcher of Mathematical Methods laboratory of the Image Processing Systems Institute of the RAS. His research interests are currently focused on image compression, space images processing, geoinformation technologies.

Глумов Николай Иванович, родился в 1962 году. В 1985 году окончил Куйбышевский авиационный институт (ныне Самарский государственный аэрокосмический университет (СГАУ)). В 1994 году защитил диссертацию на соискание степени кандидата технических наук. В настоящее время работает старшим научным сотрудником в Институте систем обработки изображений РАН. Круг научных интересов включает обработку изображений и распознавание образов, компрессию изображений, моделирование систем формирования цифровых изображений. Имеет свыше 100 публикаций, в том числе более 40 статей, две монографии (в соавторстве). E-mail: nglu@smr.ru .

Nikolay Ivanovich Glumov (b. 1962) graduated with honours (1985) from S. P. Korolyov Kuibyshev Aviation Institute (presently, S. P. Korolyov Samara State Aerospace University - SSAU). He received his Candidate in Technics (1994) degree from Samara State Aerospace University. He is the senior research scientist at the Image Processing Systems Institute of the Russian Academy of Sciences (IPSI RAS). His current research interests include image processing and pattern recognition, images compression, digital images forming systems modelling. He has more than 100 publications, including more than 40 scientific papers, 2 monographs (in co-authorship).

Поступила в редакцию 23 мая 2014 г.