Импорт грида в GRASS

[eucrariano]

Здравствуйте! Говорят, все проблемы - от незнания матчасти, и это правда. Решил воссоздать пример, реализованый в теме viewtopic.php?f=1&t=4579 и столкнулся с проблемой импорта данных. В сёрфере 9 DGM_30m_Mt.St.Helens_SRTM.grd прекрасно открывается, легко строится 3Д и т.п. А вот если экспортировать, а пробовал я по всякому, то в грассе не могу ничего открыть. То-есть, кое-что могу. Например, горизонтали в 2д представлении импортируются через mif, а вот атрибутов высоты там нет.
Вопрос, как импортировать в ГРАСС чистый грид, или во что его преобразовать, чтобы можно было работать далее?
Спасибо.

[KolesovDmitry]

Сразу скажу, что серфер не знаю, но советы даю.

GRASS может импортировать множество форматов, в первую очередь, смотрите модуль r.in.gdal, в документации к которому перечислены форматы, поддерживаемые gdal.

Потом посмотрите, есть ли общий обменный формат, в который может импортнуть серфер, а грасс прочитать.

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

Говорят, все проблемы - от незнания матчасти

К сожалению, далеко не все..

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 КБ) 6301 просмотр

Интересно другое. Команда 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). Впрочем, на импорте туда файла это никак не сказывается..

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

Как оказалось, вулкан Сент-Хеленс находится вообще-то в северной 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

Полученная область:

GRASS__DGM_30m_Mt.St.Helens_SRTM__UTM_10N.tar.gz

(500.36 КБ) 808 скачиваний

[eucrariano]

Спасибо. Снимаю шляпу !!!