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

128

ISSN 2223-1560. Известия Юго-Западного государственного университета. 2015. № 4 (61).

Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

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

1. Первичная обработка данных ДЗЗ

В настоящее время наблюдается постоянный рост использования данных дистанционного зондирования Земли (ДЗЗ) при решении прикладных задач в различных сферах: в сельском хозяйстве [1,2], образовании, муниципальном и административном управлении [3], что, в свою очередь, требует разработки методов и средств фотограмметрической обработки данных ДЗЗ для удовлетворения потребностей в оперативности и объёмах получения этих данных [4, 5].

Первичная обработка данных ДЗЗ учитывает конструктивные особенности специальной бортовой аппаратуры и может носить уникальный характер для материалов, получаемых от различных космических аппаратов [4]. Данный вид обработки осуществляется, как правило, в наземных центрах приёма информации. Методы и способы сокращения времени её выполнения подробно изложены в [4]. Одной из ключевых операций при фотограмметрической обработке данных ДЗЗ является трансформирование космического изображения местности [4, 6]. Данная операция присутствует в тракте стандартной обработки всех космических изображений [7] второго уровня и выше, осуществляется по типовым моделям [6] и представляет интерес в части сокращения времени её выполнения [8].

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

- свойства данных ДЗЗ (формат и внутренняя структура исходных растровых данных);

- свойства трансформированных изображений местности (формат и внутренняя структура результирующих растровых данных);

- используемые методы и параметры трансформирования изображения местности (типы полиномов трансформации и алгоритмов сглаживания);

- параметры аппаратной платформы (тип дискового массива и размер кластера файловой системы).

В настоящей статье рассмотрены результаты оценки влияния указанных факторов на время трансформирования космических изображений местности. Сформированы рекомендации по заданию параметров для выполнения данной операции.

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

1) исследования проводятся на цифровой фотограмметрической станции (ЦФС) на следующей платформе: процессор IntelXeon E5-2650, 2.00 ГГц; оперативная память DDR3-1600 DDR3 SDRAM (65484 Мб); дисковый массив SSD 222ГБ (RAID 1 и RAID 0); операционная система Microsoft Windows 7;

2) трансформирование изображений местности осуществляется аппроксимационными (полиномиальными) методами, хотя полученные результаты могут быть распространены и на другие методы [8];

3) размер пиксела исходного изображения местности равен размеру пиксела результирующего изображения;

4) трансформирование изображений осуществляется в одну и ту же систему координат (WGS-84);

5) при трансформировании применяется параллельный алгоритм трансформирования изображений (рис. 1).

Рассмотрим влияние указанных факторов на время трансформирования космических изображений, которое будем оценивать через производительность ЦФС в соответствующих условиях. Исходя из размера изображения местности (в Гбайт) и производительность ЦФС (в Гбайт/мин), рассчитываем время его трансформирования (мин).

Рис. 1. Блок-схема алгоритма трансформирования изображений

изображение местность дистанционный зондирование

1. Влияние свойств данных ДЗЗ

и трансформированных

изображений местности

на производительность ЦФС

Космические изображения местности (как исходные, так и полученные в результате трансформирования) обладают различными свойствами [9], основными из которых являются: размеры изображения, размер пиксела, формат файла изображения, организационная структура файла (размер блоков изображения), степень и применяемый алгоритм сжатия изображения.

Характер зависимости времени трансформирования изображения от его размеров очевиден и поэтому в рамках настоящей статьи не рассматривается.

Рассмотрим влияние сжатия изображения на производительность ЦФС при прочих равных условиях. Для этого используем в качестве исходных изображение без компрессии в формате GeoTIFF и изображения, обработанные различными алгоритмами сжатия, в форматах JPEG2000, ECW, JPEG. Зависимости производи-тельности ЦФС от размера результирующего изображения приведены на рисунке 2.

Рис. 2. Зависимость производительности ЦФС от размера результирующего изображения при различных методах сжатия

Сжатые космические изображения занимают меньший объём дискового массива, чем не сжатые, однако при их трансформировании требуется выполнение операции декомпрессии, что на порядок увеличивает время на их обработку (рис. 2).

Рассмотрим влияние организационной структуры файла изображения различных форматов на производительность ЦФС, обусловленное хранением растровых данных в файле в виде блоков, которое определяется размерами этих блоков. Одним из распространённых форматов, в котором данные ДЗЗ поступают на трансформирование, является GeoTIFF со «строчным» блоком (размеры блока равны одной строке изображения). Другим форматом является IMG, блок которого является квадратом с размерами стороны от 64 до 2048 пикселов. Изображения в формате GeoTIFF могут иметь такую же блочную структуру, если заданы соответствующие параметры при их формировании. Далее на рисунках 3 - 6 приводятся зависимости производительности ЦФС от размера входного изображения при различных вариациях размеров блоков двух указанных форматов.

Рис. 3. Зависимость производительности ЦФС от размера входного изображения при трансформировании (формат исходного и выходного изображений - GeoTIFF, тип блока совпадает)

Рис. 4. Зависимость производительности ЦФС от размера входного изображения при трансформировании для различных размеров блока (формат исходного изображения GeoTIFF с квадратным блоком, выходного - IMG, размеры блока совпадают)

Рис. 5. Зависимость производительности ЦФС от размера обрабатываемого изображения при трансформировании для различных размеров блока результирующего изображения (формат исходного изображения GeoTIFF со «строчным» блоком, выходного - IMG)

Рис. 6. Зависимость производительности ЦФС от размера обрабатываемого изображения при трансформировании для различных размеров блока (формат входного и выходного изображения - IMG, размеры блока совпадают)

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

Из анализа рисунков 3 - 6 видно, что максимальная производительность ЦФС при указанной выше конфигурации аппаратной платформы достигается при размере стороны блока выходного изображения от 1024 до 2048 пикселов и объёме изображения от 0,5 Гб до 8 ГБ.

2. Влияние используемых методов и параметров трансформирования изображений местности на производительность ЦФС

Для трансформирования изображений местности применялись аппроксимационные методы [6], в соответствии с которыми связь точек изображения с точками местности задаётся полиномами различной степени. В данной работе исследовались связи на основе полиномов 1, 2 и 3 степени, а также сплайнов. Чем выше порядок полинома при трансформировании изображения местности, тем более сложные искажения могут быть скорректированы и тем выше общая и локальная точность формируемого изображения. Трансформирование изображений на основе полиномов высокого порядка требует большего количества опорных точек и, следовательно, больше времени на обработку. Зависимость производительности ЦФС от размера обрабатываемого изображения при трансформировании для различных типов полиномов приведена на рисунке 7.

Рис. 7. Зависимость производительности ЦФС от размера входного изображения при трансформировании на основе полиномов различных типов

Как видно из рисунка 7, ЦФС демонстрирует наименьшую производительность при трансформировании изображения на основе сплайнов, однако невязка в опорных точках изображения в этом случае стремится к нулю, что дает предпосылки формирования результирующего изображения с лучшей точностью, чем при использовании полиномов других типов [6].

Поскольку в процессе трансформирования изображения изменяется его фотометрическая структура, то для получения значения яркости пиксела результирующего изображения по значениям яркости ограниченного набора пикселов исходного изображения применяется их интерполяция. Суть интерполяции заключается в использовании имеющихся данных для получения ожидаемых значений в неизвестных точках изображения местности [10]. Зависимость производительности ЦФС от размера обрабатываемого изображения при трансформировании для различных типов интерполяции пикселов приведена на рисунке 8.

Рис. 8. Зависимость производительности ЦФС от размера обрабатываемого изображения при трансформировании для различных типов интерполяции пикселов, тип дискового массива -Raid0

Следует отметить, что вид интерполяции пикселов влияет на размытие резких переходов яркости изображения и наличие алиасинга (ступенчатости контуров отображаемых объектов), то есть - на визуальное восприятие изображения местности. В задачах космического мониторинга стационарных объектов и участков поверхности Земли возникает необходимость четкого прослеживания контрастных границ этих объектов. Для решения подобных задач применяются различные виды кубической интерполяции и Ланцоша. Выполнение данных видов интерполяции является ресурсоемкой процедурой и их применение приводит к снижению производительности ЦФС, что видно из рисунка 8.

3. Влияние параметров аппаратной части ЦФС на её производительность

В данной работе не рассматривалось влияние аппаратной платформы ЦФС (частота работы процессора, количество ядер процессора, объем оперативного запоминающего устройства (ОЗУ) и пр.) на её производительность, поскольку оператор не может штатно повлиять на эти параметры (кроме замены комплектующих). Однако существует ряд параметров ЦФС, установка которых позволит повлиять на её производительность. В качестве таких параметров исследовались: тип дискового массива и размер кластера - минимального объема дискового пространства, который может быть выделен для размещения файла.

Из анализа зависимостей, изображенных на рисунках 8 и 9, следует, что использование Raid0 (объединение дисков) позволяет повысить производительность ЦФС в среднем на 10 - 15 % по сравнению с Raid1 (дублирование дисков) при прочих равных условиях.

Результаты проведенных исследований показали, что увеличение размера кластера дискового массива существенно не влияет на производительность ЦФС (рис. 10). Влияние изменения размера кластера при изменении размера блока изображения на производительность ЦФС не выявлено.

4. Формирование рекомендаций для определения параметров выполнения трансформирования космических изображений местности

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

- максимальная производительность ЦФС при трансформировании изображений достигается при размере блока от 1024 пиксела и более;

- изображения местности в формате GeoTIFF со «строчным» блоком рекомендуется использовать в случае, если их размер составляет не более 50 % от объема оперативного запоминающего устройства ЦФС в виду возможного падения её производительности при трансформировании из-за многократного обращения к блокам изображения, находящимся на дисковом массиве;

Рис. 9. Зависимость производительности ЦФС от размера обрабатываемого изображения при трансформировании для различных типов интерполяции пикселов, тип дискового массива - Raid1

Рис. 10. Зависимость производительности ЦФС от размера кластера дискового массива при трансформировании изображений с различным размером блока

- в качестве выходного формата при трансформировании космических изображений местности предпочтительно использовать GeoTIFF с размером блока 2048 пикселов;

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

- при трансформировании изображений местности рекомендуется использовать полиномы 2-й и 3-й степени, поскольку они позволяют сформировать изображение с меньшими геометрическими искажениями объектов по сравнению с полиномом 1-й степени при сопоставимом времени трансформирования;

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

- размер кластера дискового массива на время обработки изображений местности влияет незначительно.

Выводы

В работе проведен анализ влияния свойств данных ДЗЗ, свойств трансформированных изображений местности, методов и параметров процесса трансформирования и параметров аппаратной платформы на производительность цифровой фотограмметрической станции. Сформулированы рекомендации по выбору указанных параметров для повышения производительности ЦФС при сохранении точности и визуальных свойств результирующих изображений

Список литературы

1. Мирошниченко С.Ю. Выделение теней для определения высот объектов местности на космических изображениях // Известия Юго-Западного государственного университета. Серия: Управление, вычислительная техника, информатика. Медицинское приборостроение. - 2012. - №2. Ч. 2. - С. 170 - 173.

2. Ященко А.А., Мирошниченко С.Ю. Метод распознавания объектов дорожной сети на космических снимках // Известия Юго-Западного государственного университета. Серия: Управление, вычислительная техника, информатика. Медицинское приборостроение. - 2012. - №2. Ч. 1. - С. 158 - 162.

3. Панарин В.А., Кубовская О.М. Методология создания муниципальной информационной системы на основе требований к ИСОГД // Геоматика. - 2015. - № 4. - С. 58 - 66.

4. Злобин В.К., Еремеев В.В. Обработка аэрокосмических изображений. - М.: Физматлит, 2006. - 288 с.

5. Бекренев О.В. [и др.]. Возможности приема и обработки данных ДЗЗ в Консорциуме УНИГЕО // Ракетно-космическое приборостроение и информационные системы. - 2014. - Т. 1, вып. 3. - C. 66-72.

6. Шовенгердт Р.А. Дистанционное зондирование. Модели и методы обработки изображений. - М.: Техносфера, 2010. - 560 с.

7. Беленов А.В. Стандартные уровни обработки и форматы представления данных ДЗЗ из космоса. Мировой опыт // Геоматика. - 2009. - № 4. - С. 18 - 20.

8. Сонюшкин А.В. Совершенствование технологии создания ортофотопланов по космическим изображениям высокого разрешения: дис. ... канд. техн. наук / Московский государственный университет геодезии и картографии. - М., 2015.

9. Шуклин И.И., Ющенко С.П. Адаптивная распределенная система фотограмметрической обработки данных дистанционного зондирования Земли из космоса в ЕТРИС [Электронный ресурс] // Тезисы докладов 13-й Всероссийской открытой конференции «Современные проблемы дистанционного зондирования Земли из космоса» (16 - 20 ноября 2015 г.) / Институт космических исследований РАН. - М., 2015. - 1 электрон. опт. диск (CD-ROM) : зв., цв. ; 12 см.

10. Евстратова Л.Г. Трансформирование космических снимков с использованием программного комплекса ENVI: учеб. пособие. - Новосибирск: СГГА, 2008. - 53 с.

Размещено на Allbest.ru