Страница 1 из 1

MODIS - откалиброваный radiance

Добавлено: 21 июл 2016, 15:47
rastr
Добрый день,

Не нашел подобную тему на форуме. Собственно вопрос:

Возможно ли получить значения radiance сенсора MODIS на основе данных портала https://ladsweb.nascom.nasa.gov/data/search.html ?

Интересует radiance отдельных каналов в текстовом формате ASCII. Без использования пропиетарного ПО (Image Processor).

С уважением

Re: MODIS - откалиброваный radiance

Добавлено: 21 июл 2016, 16:10
rastr
UPD.

Возможно ли это осуществить в ENVI средствами модуля band math?

Re: MODIS - откалиброваный radiance

Добавлено: 22 июл 2016, 20:23
pendduduk
rastr писал(а):озможно ли это осуществить в ENVI средствами модуля band math?
То есть вы считаете что ENVI удовлетворяет
rastr писал(а):Без использования пропиетарного ПО
?

HDF View откройте продукт, найдите коэффициенты для пересчета DN в radiance и преобразуйте в любой считалке свой текстовый файл (хоть в Ёкселе ).

Re: MODIS - откалиброваный radiance

Добавлено: 23 июл 2016, 10:45
bolotoved
rastr писал(а): Возможно ли получить значения radiance сенсора MODIS на основе данных портала https://ladsweb.nascom.nasa.gov/data/search.html ?
1. Обратите внимание, что многие продукты MODIS распространяются уже в radiance. Т.е. вам остается выбрать продукт включающий необходимые вам каналы и имеющий необходимое вам пространственное/временное разрешение.
2. MODIS распространяется в виде HDF - это контейнер в котором лежат сжатые слои и их метаданные. Формат распространенный, открытый, многое ПО работает с ним напрямую.
rastr писал(а): Интересует radiance отдельных каналов в текстовом формате ASCII.
3. Когда речь заходит о конвертации форматов (если это конечно не какой-то закрытый или узкоспециальный формат), в первую очередь нужно смотреть на GDAL.

4. Особенность работы GDAL c HDF-форматом:
Вначале, нужно посмотреть, какие поднаборы (subdatasets) находятся в вашем HDF
Потом экспортировать нужный поднабор в нужный формат.

5. Смотрим поднаборы

Код: Выделить всё

gdalinfo ваш.hdf
или, если вы работаете в Linux, то можно сократить вывод:

Код: Выделить всё

gdalinfo ваш.hdf  | grep SUBDATASET
Эта команда выдаст что-то вроде этого:

Код: Выделить всё

  SUBDATASET_1_NAME=HDF4_EOS:EOS_GRID:"MOD17A3.A2010001.h09v02.055.2011276120902.hdf":MOD_Grid_MOD17A3:Gpp_1km
  SUBDATASET_1_DESC=[1200x1200] Gpp_1km MOD_Grid_MOD17A3 (16-bit unsigned integer)
  SUBDATASET_2_NAME=HDF4_EOS:EOS_GRID:"MOD17A3.A2010001.h09v02.055.2011276120902.hdf":MOD_Grid_MOD17A3:Npp_1km
  SUBDATASET_2_DESC=[1200x1200] Npp_1km MOD_Grid_MOD17A3 (16-bit unsigned integer)
  SUBDATASET_3_NAME=HDF4_EOS:EOS_GRID:"MOD17A3.A2010001.h09v02.055.2011276120902.hdf":MOD_Grid_MOD17A3:Gpp_Npp_QC_1km
  SUBDATASET_3_DESC=[1200x1200] Gpp_Npp_QC_1km MOD_Grid_MOD17A3 (8-bit unsigned integer)
У нас есть 3 поднабора данных с разрешением 1 км: Gpp_1km, Npp_1km и Gpp_Npp_QC_1km


6. Экспортируем поднабор Npp_1km

Код: Выделить всё

gdal_translate -of AAIGrid HDF4_EOS:EOS_GRID:"MOD17A3.A2010001.h09v02.055.2011276120902.hdf":MOD_Grid_MOD17A3:Npp_1km npp.asc
Где -of AAIGrid -- output format, в вашем случае ASCI
Повторяем процедуру для всех нужных каналов.