Страница 1 из 1

Проблема с выводом данных из GMT в SHP

Добавлено: 27 фев 2010, 16:32
Василий
Здравствуйте, возникла следующая проблема.
Используется пакет GMT для отрисовки карт. Необходимо конвертировать данные в shp. В документации имеется информация о возможности вывода данных в ascii или бинарный файл похожий на csv.

В качестве примера рисую карту:
pscoast -Js0/90/3i/70 -R-180/180/70/90 -Ba30/10 -Dl -W0.1p/0/0/0 -M -W > file.csv

Файл выгляди примерно так:
# Data from the low resolution GMT shoreline, borders, and rivers database
> Shore Bin # 1, Level 1
20 80.4023804074
19.6249332418 80.5053788052
19.8123140307 80.204165713
18.9925993744 80.3334096284
19.5698481727 80.1362630655
19.259937438 80.0794995041
18.7100022889 80.2020294499
17.7376974136 80.1258869306
18.2998397803 80.0204470893
...

Далее пытаюсь привести его к виду csv добавив разделители и убрав строки начинающиеся на # и на >. Далее по аналогии со статьёй http://gis-lab.info/qa/csv2shp.html пытаюсь с помощью ogr2ogr получить shp файл, но на выходе получаю нечто очень отдалённо похожее на то чтобы я мог получить с помощью:
pscoast -Js0/90/3i/70 -R-180/180/70/90 -Ba30/10 -Dl -W0.1p/0/0/0 -W > map.ps

Пожалуйста кто-нибудь подскажите как мне провести всю эту конвертацию. GMT мануалы не особо помогли=(
Файлы прилагаю.

Re: Проблема с выводом данных из GMT в SHP

Добавлено: 03 мар 2010, 08:04
Александр Мурый
А вы пробовали напрямую через OGR?
Вроде бы, у них появился драйвер GMT:
http://www.gdal.org/ogr/drv_gmt.html

как-то так:

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

ogr2ogr -f "ESRI Shapefile" output.shp input.gmt
А вот как сохранить или где искать этот векторный файл GMT (input.gmt), я и сам толком не знаю..