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

[Callisto]

Скажите, пожалуйста, векторные данные для России, которые здесь http://gis-lab.info/qa/vmap0.html
подходят для использования в GRASS?

Понимаю, что вопрос, очевидно для вас глупый, просто я еще совсем новичок в работе с любыми ГИС, не очень понимаю, как и где мне искать данным для работы. вроде начиталась разной информации уже здесь, а все равно осталось много пробелов.
Главным образом ищу карты дорог России, которые мне нужны для решение задачи Штайнера.

нашла примерное решение здесь http://grass.itc.it/gdp/html_grass63/v.net.steiner.html
но пока не разбиралась, просто интересуюсь, занимался ли кто-то в GRASS решением этой задачи?

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

насчет качеству данных: подходят, но их надо импортировать командой v.in.ogr в свой векторный формат GRASS
насчет качества этих данных для ваших задач - не уверен, если разберетесь, не забудьте рассказать.

[Callisto]

спасибо.
обязательно буду разбираться. может еще кто-то подскажет

[Callisto]

программа мне пишет, что v.in.ogr не является внутренней или внешней командой, исполняемой программой или пакетным файлом
вы не знаете, что бы это могло значить?
это у меня произошло как раз при попытке импортировать Vmap0 для России

[Callisto]

Читаю в пособии по GRASS про импортирование данных с помощью v.in.ogr:
новую область можно создать из SHAPE-файла, если имеется "правильный" файл проекции (.prj)(к сожалению это большая редкость)
Искала в своей папке Vmap0 этот файл. Не нашла. Значит ли это, что данные эти я никак не смогу импортировать в GRASS и что v.in.ogr мне тут никак не поможет?

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

попробуйте поискать в нашей папке, ссылка в начале
http://gis-lab.info/qa/vmap0.html

[KolesovDmitry]

Искала в своей папке Vmap0 этот файл. Не нашла. Значит ли это, что данные эти я никак не смогу импортировать в GRASS и что v.in.ogr мне тут никак не поможет?

Конечно не значит. Это значит, что вам придется создать область или по другому файлу, или руками (это тоже очень просто). Удобно создать область в UTM (все координаты в метрах - удобно для ваших расчетов плюс легко создается в GRASS даже новичком).
После того, как вы создадите область в UTM вам перед импортом нужно будет перепроецировать векторные данные Vmap0 из системы широта/долгота в новую систему координат (UTM).

[Callisto]

спасибо. как создать область по-другому уже разобралась.
а вот про область в UTM мне интересно. попробую сделать.

[Callisto]

пытаюсь реализовать пример http://grass.itc.it/gdp/html_grass63/v.net.steiner.html

который в конце страницы там. но со строчки v.db.update myroads col=cost val=1/75 where="label='interstate'"
программа начинает постоянно спотыкаться, все говорит, что невозможно осуществить команду...

может там какая-то ошибка в строках, которую я по незнанию не замечаю?

пробовала еще сохранить весь текст в документе типа s.bat и просто выполнить, как команду черед ГРАСС, но после выполнения второй строчки, программа просто замолкает и показывает просто пустую (готовую к работе), командную строку.

помогите разобраться, пожалуйста.

GRASS в моих неумелых руках явно работать не хочет.

[KolesovDmitry]

пытаюсь реализовать пример http://grass.itc.it/gdp/html_grass63/v.net.steiner.html
который в конце страницы там. но со строчки v.db.update myroads col=cost val=1/75 where="label='interstate'"
программа начинает постоянно спотыкаться, все говорит, что невозможно осуществить команду...

Я сейчас попробовал прогнать этот пример - все получается, никаких ошибок не выдает. И дерево получилось как на картинке у них...
А что именно у вас пишет GRASS, какая ошибка?

[Callisto]

ну, когда я просто файлом пробую загружать, то выходит
s.bat
g.copy vect=roads,myroads
Copy vector <roads@PERMANRNT> to current mapset as <myroads>
v.db.addcol myroads col="cost double precision"
_

все. дальше просто ничего не выполняет. просто пустая строка.

когда пытаюсь по строкам выполнять, получается так (пишу полностью по строкам, как в окне вывода):
g.copy vect=roads,myroads
copy vector <roads@PERMANENT> to current mapset as <myroads>
v.db.addcol map=myroads {columns=cost double precision}
v.db.update map=myroads column=cost value=1/50
v.db.update map=myroads column=cost value=1/75 where=label='interstate'
BMI-DBF driver error:
Error while executing: 'UPDATE myroads SET cost=1/75 WHERE label=interstate
'
в текстовом режиме тоже пробовала. заполняла в точности, как в примере и та же ошибка вылезает.

[Callisto]

Подскажите, пожалуйста, как после решения задачи, найти длину всего построенного пути в метрах или километрах?
http://grass.itc.it/gdp/html_grass63/v.net.steiner.html

[Callisto]

для v.net.path интересует тот же вопрос, но думаю, что для обоих случаев ищется одинаково. только вот как?

[KolesovDmitry]

Подскажите, пожалуйста, как после решения задачи, найти длину всего построенного пути в метрах или километрах?
http://grass.itc.it/gdp/html_grass63/v.net.steiner.html

Самый простой способ - воспользоваться командой v.report
Например, для приведенного вами примера можно подсчитать длину сегментов (по соотв. категориям дорог в отдельности):

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

v.report mysteiner opt=length units=kilometers

PS Скажите, а с ошибкой "BMI-DBF driver error" вы разобрались?

[Callisto]

Да. разобралась. пришлось по-другому вводить текст. тот же файл s.bat созранила просто как s, без bat. и запустила с помощью команты sh s
тоже немного прога ругалась в середине выполнения, так и не поняла над чем конкретно, но вроде выполняет все правильно.

спасибо за подсказку, сегодня попробую.

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