Продукты MODIS (Moderate - resolution imaging spectroradiometer)

1.

Лабораторная работа №3
Продукты MODIS
(Moderate-Resolution Imaging Spectroradiometer)

2.

Продукты MODIS
Префикс MOD для данных TERRA, префикс
MYD для данных AQUA, MCD для
TERRA+AQUA;
http://en.wikipedia.org/wiki/ModerateResolution_Imaging_Spectroradiometer
Представлены на синусоидальной (250, 500 и
1000 м) и градусной (0.05 и 1) сетках;
Ежедневные продукты, 8-ми, 16-ми, 32-ми
дневные, средние за месяц и год.

3.

NDVI

4.

Синусоидальная сетка L3/L4
представления данных MODIS

5.

Задание №1 MOD13
Откройте поле NDVI (Normalized Difference Vegetation Index) из файла
MOD13Q1*.hdf
Значения массива должны лежать в диапазоне [0..1]. Если это не так, то
используйте растер→калькулятор растров, для преобразования данных к
диапазону [0..1]. Здесь коэффициент равен 10000, но это не точно.
Задайте цветовую палитру полученному слою, использую следующие параметры:
0. - 0.1 коричневый;
0.1 — 1 — градиент от белого к зеленому.
Наложите векторные слои, например Алтайский край или субъекты РФ
Постройте полученное изображение в системе координат UTM44.
Идентифицируйте ярко зеленые полосы на изображении.
Изучите атрибуты файла MOD13Q1*.hdf. Для этого откройте его в HDFview. При
нажатии ЛКМ на объект в левой панели HDFview в нижней панели
отображаются атрибуты этого объекта.
Определите в какой период времени была сделана спутниковая съемка?

6.

Задание №2 MOD09(GA/Q1)
Откройте данные каналов b01 и b02 из файла MOD09G*.hdf, у вас должно
получиться два слоя.
Значения в массивах должны лежать в диапазоне [0..1]. Если это не так, то
используйте растер→калькулятор растров, для преобразования данных
к диапазону [0..1]. Здесь коэффициент равен 10000, но это не точно.
Рассчитайте NDVI по формуле NDVI=(B02-B01)/(B02+B01), для этого
используйте растер→калькулятор растров.
Скопируйте палитру цветов из слоя, полученного в задании №1, диапазон
значений должен лежать в диапазоне [0..1].
Изучите атрибуты файла MOD09G*.hdf. Для этого откройте его в HDFview.
Определите в какой период времени была сделана спутниковая
съемка?
Сопоставьте полученный результат с полем NDVI из MOD13. Должны ли
они быть похоже друг на друга?
Постройте полученное изображение в системе координат UTM44.

7.

Задание №3 MCD12Q1
Откройте поле Land_Cover_Type_5 из файла MCD12Q1*h22v03*.hdf для 2002 и 2012 годов.
Описание данных Land_Cover_Type_5 вы сможете найти на следующем слайде. Данные каких
типов поверхности могут быть полезны при работе с природоохранными объектами на
территории Алтайского края?
Задайте палитру цветов для этих полей. Количество цветов должно быть кратно количество
категорий данных. Соответствуют ли, в общем, полученные данные типам поверхности,
которые характерны для территории Алтайского края? Для подсказки можете наложить
векторные слои, например водные объекты. Для удобства вы можете работать в UTM44.
Определите к какому типу из Land_Cover_Type_5 соответствую ленточные боры, Описание
данных Land_Cover_Type_5 вы сможете найти на следующем слайде.
Исправьте палитру цветов так, чтобы ленточные боры были обозначены, например красным
цветом, а остальные категории - белым.
У вас должно получиться 2 слоя, один для 2002 года, а второй для 2012 года. Замените в
одном из слоев красный цвет на синий.
Включите прозрачность для всех слоев равную 50%. Оцените, как изменилась площадь боров
с 2002 по 2012 годы.
Добавьте в проект файлы MCD12Q1*h22v03*.hdf и выполните для них аналогичные действия.
Изменилась ли ваша оценка площадь боров.
Выполните аналогичный анализ для пахатных земель, кустарников и лугов.

8.

Описание поля Land_Cover_Type_5
0
1
2
3
4
5
6
7
8
9
10
11
254
255
Water
Evergreen Needleleaf trees
Evergreen Broadleaf trees
Deciduous Needleleaf trees
Deciduous Broadleaf trees
Shrub
Grass
Cereal crops
Broad-leaf crops
Urban and built-up
Snow and ice
Barren or sparse vegetation
Unclassified
Fill Value

9.

FAQ
Где взять данные MODIS?
ftp://ladsftp.nascom.nasa.gov/
http://e4ftl01.cr.usgs.gov/
Где найти описание продуктов MODIS?
https://lpdaac.usgs.gov/products/modis_products_
table