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

Ответить
Аватара пользователя
eucrariano
Активный участник
Сообщения: 104
Зарегистрирован: 26 авг 2009, 12:09
Репутация: 0
Откуда: Ukraine, Chernivtsy

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

Сообщение eucrariano » 07 мар 2011, 16:34

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

KolesovDmitry
Гуру
Сообщения: 810
Зарегистрирован: 22 авг 2007, 14:58
Статьи: 18
Проекты: 4
Репутация: 122
Откуда: Казань

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

Сообщение KolesovDmitry » 07 мар 2011, 17:22

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

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

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

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Статьи: 3
Проекты: 5/1
Репутация: 786
Ваше звание: званий не имею
Откуда: Москва

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

создала у меня область/набор с "гридом".
2011-03-07-171633_1024x768_scrot.png
"грид" с цветовой картой haxby и с теневой отмывкой
2011-03-07-171633_1024x768_scrot.png (181.86 КБ) 3456 просмотров
Интересно другое. Команда 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). Впрочем, на импорте туда файла это никак не сказывается..
Редактор материалов, модератор форума

Александр Мурый
Гуру
Сообщения: 5173
Зарегистрирован: 26 сен 2009, 16:26
Статьи: 3
Проекты: 5/1
Репутация: 786
Ваше звание: званий не имею
Откуда: Москва

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


Полученная область:
GRASS__DGM_30m_Mt.St.Helens_SRTM__UTM_10N.tar.gz
(500.36 КБ) 454 скачивания
Редактор материалов, модератор форума

Аватара пользователя
eucrariano
Активный участник
Сообщения: 104
Зарегистрирован: 26 авг 2009, 12:09
Репутация: 0
Откуда: Ukraine, Chernivtsy

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

Сообщение eucrariano » 08 мар 2011, 13:54

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

Ответить

Вернуться в «GRASS»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 1 гость