Форумы GIS-Lab.info

Здравствуйте, уважаемое сообщество! Я на сайте новичок, поэтому простите если мой вопрос будет уж слишком глупый.

Задача такая у меня есть hdf файл( Modis Level 1B) мне нужно сделать с него JPEG файл. Как я понимаю это можно сделать с помощью gdal_translate. Но в файле же содержится несколько датасетов. Вот они:
Subdatasets:
SUBDATASET_1_NAME=HDF4_EOS:EOS_SWATH:"/home/alexandr/Загрузки/test.hdf":MODIS_SWATH_Type_L1B:EV_500_RefSB
SUBDATASET_1_DESC=[5x4060x2708] EV_500_RefSB MODIS_SWATH_Type_L1B (16-bit unsigned integer)
SUBDATASET_2_NAME=HDF4_EOS:EOS_SWATH:"/home/alexandr/Загрузки/test.hdf":MODIS_SWATH_Type_L1B:EV_500_RefSB_Uncert_Indexes
SUBDATASET_2_DESC=[5x4060x2708] EV_500_RefSB_Uncert_Indexes MODIS_SWATH_Type_L1B (8-bit unsigned integer)
SUBDATASET_3_NAME=HDF4_EOS:EOS_SWATH:"/home/alexandr/Загрузки/test.hdf":MODIS_SWATH_Type_L1B:EV_250_Aggr500_RefSB
SUBDATASET_3_DESC=[2x4060x2708] EV_250_Aggr500_RefSB MODIS_SWATH_Type_L1B (16-bit unsigned integer)
SUBDATASET_4_NAME=HDF4_EOS:EOS_SWATH:"/home/alexandr/Загрузки/test.hdf":MODIS_SWATH_Type_L1B:EV_250_Aggr500_RefSB_Uncert_Indexes
SUBDATASET_4_DESC=[2x4060x2708] EV_250_Aggr500_RefSB_Uncert_Indexes MODIS_SWATH_Type_L1B (8-bit unsigned integer)
SUBDATASET_5_NAME=HDF4_EOS:EOS_SWATH:"/home/alexandr/Загрузки/test.hdf":MODIS_SWATH_Type_L1B:EV_250_Aggr500_RefSB_Samples_Used
SUBDATASET_5_DESC=[2x4060x2708] EV_250_Aggr500_RefSB_Samples_Used MODIS_SWATH_Type_L1B (8-bit integer)


Не очень понимаю что эти данные означают, и какой из них брать чтобы получить нормальное RGB изображения.
Заранее спасибо за помощь!

Какие именно данные вы импортируете? MOD02HKM? Как выглядит имя исходного непереименованного файла? Откуда вы его брали?

Если вы новичок, вам возможно будет проще использовать MRT Swath, к нему есть графический интерфейс пользователя.

Описания для командной строки и спец. инструмента можно найти по ссылкам:
http://gis-lab.info/qa/modis-l1-import.html
http://gis-lab.info/qa/modis-l1-import-tool.html

Пример для gdal будет выглядеть примерно так (это другие данные но идея должна быть понятна):

данные брал от сюда http://ladsweb.nascom.nasa.gov/data/search.html.
Данные MOD02HKM.

Как я понял файл нельзя переименовывать верно?

если вы не хотите, чтобы мы играли в телепатов пытаясь угадать тип данных, то не стоит
для всего остального значения это не имеет

Спасибо! В следующий раз учту!

Все равно немножко не понял.. Вот я вытаскиваю допустим из hdf все датасеты в отдельные tif файлы.. как мне потом сделать из этого нормальную RGB картинку?

Раз, два

Спасибо за помощь. Разобрался.

Теперь не могли бы вы сказать как осуществить это все что бы у меня еще была геопривязка!

А разве при использовании GDAL привязка не сохраняется? Посмотрите свойства ваших растров при помощи gdalinfo. Например, если растр называется b1.tif, команда будет следующая

импортировать gdal'ом MOD02QKM я бы не рекомендовал, последний раз когда я смотрел там было слишком много глюков, лучше MRT Swath

но если уж решите gdal'ом, то см. ключ -geolocate

Спасибо!!!

А можете все таки подсказать как с помощью MRT Swath выгрузить правильно с гео-привязкой.

Еще вопрос. Вытаскиваю я данные МРТ нужные мне три канала для РГБ. Потом с помощью gdal_merge.py их склеиваю, но изображение получается все равно серое. Знаю есть много еще средств чтобы это сделать. Но мне нужно только gdal_merge.py чтобы потом это все заскриптовать. Так как все это сделать, чтобы сохранить геопривязку и создать коректный РГБ файл.

Покажите команду, которой выполняете склейку.