Страница 1 из 1
Импорт грида в GRASS
Добавлено: 07 мар 2011, 16:34
eucrariano
Здравствуйте! Говорят, все проблемы - от незнания матчасти, и это правда. Решил воссоздать пример, реализованый в теме
viewtopic.php?f=1&t=4579 и столкнулся с проблемой импорта данных. В сёрфере 9 DGM_30m_Mt.St.Helens_SRTM.grd прекрасно открывается, легко строится 3Д и т.п. А вот если экспортировать, а пробовал я по всякому, то в грассе не могу ничего открыть. То-есть, кое-что могу. Например, горизонтали в 2д представлении импортируются через mif, а вот атрибутов высоты там нет.
Вопрос, как импортировать в ГРАСС чистый грид, или во что его преобразовать, чтобы можно было работать далее?
Спасибо.
Re: Импорт грида в GRASS
Добавлено: 07 мар 2011, 17:22
KolesovDmitry
Сразу скажу, что серфер не знаю, но советы даю.
GRASS может импортировать множество форматов, в первую очередь, смотрите модуль
r.in.gdal, в документации к которому перечислены форматы, поддерживаемые gdal.
Потом посмотрите, есть ли общий обменный формат, в который может импортнуть серфер, а грасс прочитать.
Re: Импорт грида в GRASS
Добавлено: 07 мар 2011, 17:24
Александр Мурый
eucrariano писал(а):Говорят, все проблемы - от незнания матчасти
К сожалению, далеко не все..
eucrariano,
r.in.gdal не помог?
Команда
Код: Выделить всё
r.in.gdal -e in=/home/amuriy/DGM_30m_Mt.St.Helens_SRTM.grd out=DGM_30m_Mt.St.Helens_SRTM location=DGM_30m_Mt.St.Helens_SRT
создала у меня область/набор с "гридом".
- "грид" с цветовой картой haxby и с теневой отмывкой
- 2011-03-07-171633_1024x768_scrot.png (181.86 КБ) 6171 просмотр
Интересно другое. Команда
gdalinfo DGM_30m_Mt.St.Helens_SRTM.grd выдаёт такое:
Код: Выделить всё
Driver: GSBG/Golden Software Binary Grid (.grd)
Files: DGM_30m_Mt.St.Helens_SRTM.grd
Size is 312, 458
Coordinate System is `'
Origin = (557955.000000000000000,5121855.000000000000000)
Pixel Size = (30.000000000000000,-30.000000000000000)
Corner Coordinates:
Upper Left ( 557955.000, 5121855.000)
Lower Left ( 557955.000, 5108115.000)
Upper Right ( 567315.000, 5121855.000)
Lower Right ( 567315.000, 5108115.000)
Center ( 562635.000, 5114985.000)
Band 1 Block=312x1 Type=Float32, ColorInterp=Undefined
Min=703.896 Max=2522.912
NoData Value=1.70141000918782799e+38
Т.е. система координат не определяется GDAL'ом. Соответственно, область в грасс создаётся
XY location (unprojected). Впрочем, на импорте туда файла это никак не сказывается..
Re: Импорт грида в GRASS
Добавлено: 07 мар 2011, 18:12
Александр Мурый
Как оказалось,
вулкан Сент-Хеленс находится вообще-то в северной
10-й зоне UTM :^)
Тогда делаем так:
Код: Выделить всё
gdal_translate -a_srs EPSG:32610 -of GTiff DGM_30m_Mt.St.Helens_SRTM.grd DGM_30m_Mt.St.Helens_SRTM.tif
, т.е. заменяем "неопознанную" систему координат/проекцию правильной (
http://spatialreference.org/ref/epsg/32610/), попутно конвертируем из сёрфер-грида в геотифф.
Теперь импортируем геотифф в GRASS из уже существующих области/набора:
Код: Выделить всё
r.in.gdal -e in=/home/amuriy/DGM_30m_Mt.St.Helens_SRTM.tif out=DGM_30m_Mt.St.Helens_SRTM location=DGM_30m_Mt.St.Helens_SRTM__UTM_10N
Полученная область:
Re: Импорт грида в GRASS
Добавлено: 08 мар 2011, 13:54
eucrariano
Спасибо. Снимаю шляпу
!!!
