Hdf convert to JPEG

[Alex_12]

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

Задача такая у меня есть 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 будет выглядеть примерно так (это другие данные но идея должна быть понятна):

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

gdalwarp -overwrite -of GTiff HDF4_EOS:EOS_SWATH:"110601_081441.hdf":VNIR_Band1:ImageData b1.tif 

[Alex_12]

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

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

[Максим Дубинин]

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

[Alex_12]

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

[Alex_12]

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

[Voltron]

Раз, два

[Alex_12]

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

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

[Voltron]

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

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

gdalinfo b1.tif

[Максим Дубинин]

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

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

[Alex_12]

Спасибо!!!

[Alex_12]

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

[Alex_12]

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

[Voltron]

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

[Максим Дубинин]

viewtopic.php?f=30&t=9764&st=0&sk=t&sd= ... 000#p57878