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

[rastr]

Добрый день,

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

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

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

С уважением

[rastr]

UPD.

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

[pendduduk]

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

То есть вы считаете что ENVI удовлетворяет

Без использования пропиетарного ПО

?

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

[bolotoved]

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

1. Обратите внимание, что многие продукты MODIS распространяются уже в radiance. Т.е. вам остается выбрать продукт включающий необходимые вам каналы и имеющий необходимое вам пространственное/временное разрешение.
2. MODIS распространяется в виде HDF - это контейнер в котором лежат сжатые слои и их метаданные. Формат распространенный, открытый, многое ПО работает с ним напрямую.

Интересует 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
Повторяем процедуру для всех нужных каналов.