GRASS HDF файлы

[Rumato]

Добрый день, сразу извиняюсь за наверное не сильно умный вопрос.

Мне нужно открыть hdf-файл, вида MOD09GA.A2011185.h22v03.005.2011210213142.hdf, прибора MODIS, и потом данне из его объединить с другими данными, т.е. слой от модиса, слой всякие районы и т.д.

Теперь вопрос: как открыть файл HDF и использовать его данные как слои? Я извиняюсь за глупость, но только начал разбираться в этой тематике.

Как я понимаю нужно определить максимальный георафический обхват, нужно определить спектральные каналы и т.д. Если не сложно, подскажите пожалуйста, как пошагово работать с этими данными в GRASS.


Заранее спасибо за помощь.

[Александр Мурый]

Так, давайте по очереди.
Имеем дело с продуктом MODIS MOD09GA (Surface Reflectance Daily L2G Global 1km and 500m).

Теперь вопрос: как открыть файл HDF и использовать его данные как слои?

Для начала нужно вытащить данные из HDF в более понятный для ГИС растровый формат и более употребимую проекцию (например, GeoTIFF в WGS84). Сделать это можно, например, с помощью MODIS Reprojection Tool. Далее уже можно импортировать в GRASS и работать там с растрами.

[Rumato]

Я пользовался MODIS Reprojection Tool, но меня смущает то что геотифовские изображения получаются чёрно-белые, их можно как-нибудь сделать в виде того же RGB слоя? как-бы придать какое-либо изображение?

[Александр Мурый]

Я пользовался MODIS Reprojection Tool, но меня смущает то что геотифовские изображения получаются чёрно-белые, их можно как-нибудь сделать в виде того же RGB слоя? как-бы придать какое-либо изображение?

В каких программах вы просматриваете геотифы? Они сами по себе не станут "цветными". Надо для начала вытащить из них данные в нормальном виде (сейчас они там в виде целых чисел, насколько я понимаю). И вообще, выложите образцы HDF и GeoTIFF-файлов, пож-та.

[Rumato]

Я открывал их просто при помощи GIMP-a или при помощи QGIS и там и там чёрно-белые, вот сами файлы
http://scienceview.ru/temp/MOD09GA.A201 ... 213142.hdf

http://scienceview.ru/temp/data.zip

В архиве tif-ы что я наконвертировал.

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

склейте получившиеся геотифы в композиты примерно так как рассказано здесь
потом показывайте как RGB

http://gis-lab.info/qa/qgis-landsat-merge.html

[Александр Мурый]

Если в GRASS, то для создания RGB-композитов используйте модуль r.composite, хотя композиты можно просматривать и без их записи.

[Rumato]

Я сделал всё как по руководству, что посоветовал Максим Дубинин , но вот картинка получилась не совсем хорошей, т.е. http://scienceview.ru/temp/test.jpg


Как-нибудь это можно побороть? Я вроде пытался править и настраивать показ, всё равно примерно так же с минимальными изменениями.

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

Заранее спасибо за помощь.

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

версия QGIS? выложите еще раз композит посмотреть.

[Rumato]

Версия 1.8.0 композит это что?(не ругайтесь) если начальные файлы, то всё всё так же по темже ссылкам.

[Александр Мурый]

Версия 1.8.0 композит это что?(не ругайтесь) если начальные файлы, то всё всё так же по темже ссылкам.

Композит — это в данном контексте многоканальное изображение.

Попробуйте в свойствах растра на вкладке "Стиль" изменить "Улучшение контраста" на "Stretch to MinMax". Картинка, приведённая раньше, получается, если стоит "No enhancement".

Вот пример комбинации каналов 1-2-3:

hdf_rgb.png (423.41 КБ) 12286 просмотров

[Rumato]

Александр, Максим, большое спасибо, Вы мне очень помогли! И вопрос как бы в догонку: а вот это всё, объединение файлов tif, затем выбор каналов, перепроецирование(если нужно), можно ли это в консольном режиме проделать с QGIS, чтобы была возможность запускать это и получать результаты в автоматическом режиме. Возможно так?

[Александр Мурый]

а вот это всё, объединение файлов tif, затем выбор каналов, перепроецирование(если нужно), можно ли это в консольном режиме проделать с QGIS, чтобы была возможность запускать это и получать результаты в автоматическом режиме.

Все процедуры в QGIS из меню "Растр" — интерфейс к утилитам GDAL. Когда вы делаете что-то, вы можете видеть в нижней части окна соответствующую команду для запуска в консоли или в скрипте. Объединение в композит — gdal_merge, выбор каналов с созданием нового растра — gdal_translate, перепроецирование — gdalwarp. Смотрите примеры в этой статье.

[Rumato]

Большое спасибо, буду разбираться.