Выделение растровых контуров объектов в GRASS

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

Обсуждение статьи "Выделение растровых контуров объектов в GRASS"

http://gis-lab.info/qa/grass-gradient.html

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

Спасибо за познавательную статью:))

Сейчас как раз разбираюсь с r.mapcalc -- хочу использовать его аля-EasyTrace для автоматической векторизации изогипс рельефа с генштабовских топокарт. Чем не растровые контуры..

[KolesovDmitry]

Спасибо за познавательную статью:))

Сейчас как раз разбираюсь с r.mapcalc -- хочу использовать его аля-EasyTrace для автоматической векторизации изогипс рельефа с генштабовских топокарт. Чем не растровые контуры..

Ну как, что-нибудь получилось? Особенно интересует, как разобрались с цветоделением.

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

Честно говоря, пока никак. Почитал про технологии, используемые в EasyTrace. Напрямую в GRASS все это реализовать вряд ли удасться, да и нужно ли? а если серьезно -- с математикой я не дружу и тем более не программист. Так что с цветоделением пока глухо.

Я импортировал кусочек в *jpg
Сделал RGB на каждый канал по 3 уровня. Получился растр с 26 категориями:

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

r.composite -d -c red="map.red" green="map.green" blue="map.blue" levels=3 output="map_RGB"

Изолиниям в целом отвечают категории 4 и 5. Извлекаем их и объединяем:

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

r.mapcalc "map_RGB_4=if(map_RGB == 4, map_RGB, null())"
r.mapcalc "map_RGB_5=if(map_RGB == 5, map_RGB, null())"
r.patch input="map_RGB_4,map_RGB_5" output="map_RGB_4_5"

Утоняем растровые "изолинии" -- r.thin
Перегоняем в вектор -- r.to.vect, делаем полилинии (?) -- v.build.polylines ( как видно, толком не сработало)

ДО

map_RGB.png (111.95 КБ) 13603 просмотра

ПОСЛЕ

map_RGB_4_5_thin_polylines.png (37 КБ) 13603 просмотра

P.S. Сделал затем как в сабжевой статье -- получилась какая-то каша

[KolesovDmitry]

Спасибо за подробности, наверное, есть смыл обсуждать эту тему далее в другой ветке?