Страница 2 из 2
Re: Векторные данные для GRASS и решение задачи Штайнера
Добавлено: 21 апр 2009, 08:50
KolesovDmitry
Callisto писал(а):
Такой еще вопрос, в примере координаты точек в метрах даны? думаю, можно ли тоже самое с широтой и долготой? и как оценить погрешность при пересчете из географических координат в декартовы?
Думаю, что можно и в градусах. Уверен почти на 100% (почти - потому, что сам еще не пробовал).
Как оценить погрешность пересчета в декартову систему (видимо, вы имеете в виду UTM?)? Хорошего способа не знаю, могу предложить только взять и посчитать несколькими разными программами и сравнить. Но вообще говоря, погрешность не должна превышать нескольких метров, уж точно не сильно больше чем десять. Наверное, по погрешности пересчета есть смысл задать вопрос в отдельной теме - вам ответят более знающие люди.
Re: Векторные данные для GRASS и решение задачи Штайнера
Добавлено: 21 апр 2009, 18:00
Максим Дубинин
как оценить погрешность при пересчете из географических координат в декартовы?
субметровая
Re: Векторные данные для GRASS и решение задачи Штайнера
Добавлено: 21 апр 2009, 21:14
Callisto
Sim, это, я так понимаю, так погрешность называется. а как определить ее? какие-то определенные методы есть или как предложил KolesovDmitry?
Дмитрий, спасибо. я попробую с градусами. Но если вы тоже проверите, то напишите сюда, пожалуйста, если получится. Только если будете пробовать, то пишите прямо как в тексте программы.
Спасибо за помощь всем!
Re: Векторные данные для GRASS и решение задачи Штайнера
Добавлено: 22 апр 2009, 11:10
KolesovDmitry
Callisto писал(а):Но если вы тоже проверите, то напишите сюда, пожалуйста, если получится. Только если будете пробовать, то пишите прямо как в тексте программы.
Попробовал в градусах, абсолютно никаких отличий - все как и ожидалось
"Программу" тоже выкладываю, хотя она и не отличается от того, что уже было приведено в ссылке выше. Единственное что по-другому - я не стал заморачиваться с разными весами ребер для разных типов дорог. Регион создан в формате широта-долгота на 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
Считал долго, около получаса, так что, если будете экспериментировать, берите дороги не на всю Россию, как я, а на свой кусок