Векторные данные для GRASS и решение задачи Штайнера

[KolesovDmitry]

Такой еще вопрос, в примере координаты точек в метрах даны? думаю, можно ли тоже самое с широтой и долготой? и как оценить погрешность при пересчете из географических координат в декартовы?

Думаю, что можно и в градусах. Уверен почти на 100% (почти - потому, что сам еще не пробовал).
Как оценить погрешность пересчета в декартову систему (видимо, вы имеете в виду UTM?)? Хорошего способа не знаю, могу предложить только взять и посчитать несколькими разными программами и сравнить. Но вообще говоря, погрешность не должна превышать нескольких метров, уж точно не сильно больше чем десять. Наверное, по погрешности пересчета есть смысл задать вопрос в отдельной теме - вам ответят более знающие люди.

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

как оценить погрешность при пересчете из географических координат в декартовы?

субметровая

[Callisto]

Sim, это, я так понимаю, так погрешность называется. а как определить ее? какие-то определенные методы есть или как предложил KolesovDmitry?
Дмитрий, спасибо. я попробую с градусами. Но если вы тоже проверите, то напишите сюда, пожалуйста, если получится. Только если будете пробовать, то пишите прямо как в тексте программы.
Спасибо за помощь всем!

[KolesovDmitry]

Но если вы тоже проверите, то напишите сюда, пожалуйста, если получится. Только если будете пробовать, то пишите прямо как в тексте программы.

Попробовал в градусах, абсолютно никаких отличий - все как и ожидалось "Программу" тоже выкладываю, хотя она и не отличается от того, что уже было приведено в ссылке выше. Единственное что по-другому - я не стал заморачиваться с разными весами ребер для разных типов дорог. Регион создан в формате широта-долгота на wgs84

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

# Импортирую железные дороги из vmap0
v.in.ogr -o dsn=/home/Dima/laboro/GIS/DATA/vmap0/RF/Транспорт/ layer=trans-railroad-l out=RailRoads

#
# При помощи v.digit создаю несколько (6) точек на карте,
# тыкая наугад, но стараясь попасть в дорогу
#

# Создаю сеть
v.net RailRoads points=centroj out=RailRoads_net op=connect thresh=0.01
# Поиск дерева
v.net.steiner RailRoads_net tcat=1-6 out=test

Считал долго, около получаса, так что, если будете экспериментировать, берите дороги не на всю Россию, как я, а на свой кусок