Несмотря на существенный возраст приборов MODIS, установленных на платформах Terra (EOS AM-1; запущена 18 декабря 1999 г.) и Aqua (EOS PM-1; запущена 4 мая 2002 г.), их данные не теряют актуальность и в настоящее время, а может быть даже и приобретают ее, ввиду наличия находящегося в общественном достоянии архива собранных данных, покрывающего более чем десять лет работы данных приборов.
Однако, использование этих данных совместно с распространенным программным обеспечением связано с рядом сложностей (отчасти связанных с дефектами метаданных, отчасти — с несколько необычной структурой используемого в программе EOS формата данных HDF-EOS.) Цель данного учебника — познакомить с некоторыми из типов результатов обработки данных MODIS, на примере взаимодействия с ними с использованием программного пакета QGIS.
Для публикуемых NASA результатов обработки данных MODIS до уровня 3 («физические величины на регулярной сетке») используется определенная номенклатура, которую рассмотрим на следующих двух примерах.
MCD12C1.A2001001.051.2012264184830.hdf
MCD
результаты совместной обработки данных MODIS/Terra и MODIS/Aqua.
12C1
код продукта MODIS — типы подстилающей поверхности.
A2001001
дата начала обрабатываемого временного ряда — 1-й день 2001 года (1 января 2001.)
051
версия обрабатывающего кода 5.1.
2012264184830
дата и время выполнения обработки — 264-й день 2012 года (20 сентября), 18:48:30 UTC.
MYD13A1.A2013249.h22v03.005.2013266052613.hdf
MYD
результаты обработки данных MODIS/Aqua.
MOD
результаты обработки данных MODIS/Terra.
13A1
код продукта MODIS — вегетационные индексы.
A2013249
дата начала обрабатываемого временного ряда — 249-й день 2013 года (9 сентября 2013.)
Чтобы получить общее представление о данных MODIS уровня 3, рассмотрим данные в системе координат «широта, долгота» на сетке 0.05°. (Данная сетка ориентирована на использование в задачах моделирования климата и, поэтому, также известна как CMG — англ. climate modelling grid.)
Загрузим () в QGIS, к примеру, растр Nadir_Reflectance_Band1 (коэффициенты спектральной яркости, приведенные к стандартной геометрии освещения—наблюдения; канал MODIS 1 — «красный») сетки MCD_CMG_BRDF_0.05Deg файла MCD43C4.A2013169.005.2013191204312.hdf.
Установим текущую область по границам слоя, а также изменим стиль на одноканальное псевдоцветное и загрузим цветовую картуreflectance-colors-10000.txt. Результат должен получиться похожим на изображение справа.
Здесь уместно упомянуть, что средствами QGIS невозможно объединить три отдельных растра для получения изображения «в истинных цветах» (каналы MODIS 1, 4, 3.) Кроме того, в описании системы координат данного набора допущена ошибка: несложно убедиться, что с точки зрения QGIS, охват слоя — 0.1° × 0.05°, а вовсе не ожидаемый 360° × 180°.
Исправить оба недостатка можно объединив отдельные растры в .vrt-файл, подобно:
С другой стороны, применить цветовую карту reflectance-colors-10000.txt к каждому из результирующих каналов средствами QGIS окажется невозможным. Впрочем, соответствующая информации также может быть внесена в .vrt-файл.
Вновь установим текущую область по границам слоя, а также изменим стиль на одноканальное псевдоцветное и загрузим цветовую картуmcd12q1-lct1-colors.txt. Ожидаемый результат представлен на изображении справа.
Наконец, рассмотрим растр CMG 0.05 Deg 16 days NDVI (нормализованный разностный вегетационный индекс; англ. NDVI) сетки MODIS_Grid_16Day_VI_CMG файла MYD13C1.A2013201.005.2013218212020.hdf.
Вновь установим текущую область по границам слоя, а также изменим стиль на одноканальное псевдоцветное и загрузим цветовую картуndvi-colors-10000.txt. Ожидаемый результат представлен на изображении справа.