Форумы GIS-Lab.info

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

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

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

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


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

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

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

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

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

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

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

Я открывал их просто при помощи 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, хотя композиты можно просматривать и без их записи.

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


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

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

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

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

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

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

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

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

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

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

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

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

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