Получение высоты точки по координатам

[Soitanen]

Помогите решить задачу:
Имеются координаты точки, имеется 3-х секундный SRTM *.hgt
Каким инструментом можно получить высоту точки с заданными координатами по имеющемуся SRTM файлу? Использую Ubuntu, желательно, чтобы это можно было сделать из командной строки, так как в итоге будет скрипт, и огромное количество таких вот точек.

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

При всём перечисленном выше вам сам Бог велел использовать GRASS. Автоматизировать при умении можно будет буквально всё.

Алгоритм действий:
- импорт векторных точек (v.in.ogr); м.б., понадобится ещё создать таблицу атрибутов (v.db.addtable) или поле в таблице атрибутов для высот (v.db.addcol);
- импорт данных SRTM в GRASS (r.in.srtm), а ещё лучше брать SRTM v4 в геотифф (тогда импорт через r.in.gdal);
- задаём текущий регион по нужному растру:

Код: Выделить всё

g.region rast=elevation

- получаем высоты точек, попадающих в охват заданного растра (v.what.rast):

Код: Выделить всё

v.what.rast vect=pnts raster=elevation column=heights

[Soitanen]

Ух ты! Прочитал 3 раза, нифига не понял, так как GRASS раньше никогда не использовал, пошел искать статьи Вечером буду мучить скрипт тогда. Спасибо!

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

влезу с не-GRASS решением)

http://gis-lab.info/qa/extract-values-rasters.html

растры перед этим объединить в vrt если нужно.

[Soitanen]

влезу с не-GRASS решением)

http://gis-lab.info/qa/extract-values-rasters.html

растры перед этим объединить в vrt если нужно.

А скачанный SRTM в виде *.hgt файлов ведь не является растром? У меня скрипт качает именно такие данные (они потом конвертируются в FlightGear).

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

А скачанный SRTM в виде *.hgt файлов ведь не является растром?

SRTM HGT — растровый формат. Для исп-я упомянутого Максимом скрипта можно конвертировать HGT в GeoTIFF (с помощью GDAL).

[Soitanen]

У меня точки не в shp формате, они текстовым файлом идут, так что вариант прикоснуться к GRASS будет более правильным.

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

У меня точки не в shp формате, они текстовым файлом идут

На всякий: Конвертация данных из CSV в SHP и обратно с OGR

[Дмитрий Барышников]

Еще можно cs2cs использовать, но этот инструмент выдает высоту по гравитационной модели - ему рельеф не нужен.