Преобразование линий в точки

[KartaBY]

Добрый день.
Подскажите как можно преобразовать линию в точку?
Существующий инструмент "извлечение узлов" создает точки на концах линии, а мне надо создать одну точку на центре линии.

Спасибо.

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

В "чистом" QGIS, боюсь, никак. Можно с помощью GRASS Tools (модуль v.to.points).

[KartaBY]

Н-да. Не хочет ГРАСС импортировать шейпы хоть убей.
Но это уже совсем другая история (ветка форума).

[rhot]

спрашивайте, поможем =)

[KartaBY]

Да собственно и спрашивать пока нечего.
Выдает ошибки, но прочитать их нет возможности, одни кракозябры.
Надо с ними бороться, так что курю пока форум.

[KartaBY]

Вот описана моя проблема.
Решения пока нет. (
viewtopic.php?f=29&t=8210

[Voltron]

У вас все линии содержат только по 2 узла и прямые? В таком случае задача решается с использованием школьного курса геометрии.

Если же у вас ломанные, то не совсем понятно как вы собрались находить «центр линии».

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

Если же у вас ломанные, то не совсем понятно как вы собрались находить «центр линии».

Поделить длину линии пополам и отступить это расстояние вдоль линии от одной из конечных точек?

[Voltron]

Поделить длину линии пополам и отступить это расстояние вдоль линии от одной из конечных точек?

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

[KartaBY]

У вас все линии содержат только по 2 узла и прямые?

да

ps. Школьный курс геометрии не говорит как этот процесс сделать в QGIS, а в ручную 150000 линий я буду до второго пришествия преобразовывать.

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

KartaBY, если возможно выложить куда-то данные, могу помочь с точками (в GRASS). В QGIS без спец. программирования вы это не сделаете, как я понял.

[Voltron]

ps. Школьный курс геометрии не говорит как этот процесс сделать в QGIS

Школьный курс говорит как найти середину отрезка. Внезапно.

а в ручную 150000 линий я буду до второго пришествия преобразовывать.

Две минуты, включая набор и проверку на тестовом слое

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

from qgis.core import *
layer = qgis.utils.iface.mapCanvas().currentLayer()
vl = QgsVectorLayer("Point", "temporary_points", "memory")
pr = vl.dataProvider()
layer.select()
ft = QgsFeature()
while layer.nextFeature(ft):
  tmp = ft.geometry().asPolyline()
  x = (tmp[0].x() + tmp[1].x())/2.0
  y = (tmp[0].y() + tmp[1].y())/2.0
  f = QgsFeature()
  f.setGeometry(QgsGeometry.fromPoint(QgsPoint(x, y)))
  pr.addFeatures([f])

QgsMapLayerRegistry.instance().addMapLayer( vl )

Надеюсь, вы знаете что такое Консоль Python в QGIS

[KartaBY]

KartaBY, если возможно выложить куда-то данные, могу помочь с точками (в GRASS). В QGIS без спец. программирования вы это не сделаете, как я понял.

Вот для примера один из кусков.
Мне требуется преобразовать каждый линейный объект в отдельную точку ( в середине линии) с сохранением всех атрибутов.

[KartaBY]

Две минуты, включая набор и проверку на тестовом слое

Спасибо что уделили время для моей проблемы,
но ... я не смог получить правильный точечный слой с помощью вашего кода.
Выше приложен пример моих шейпов, может дело в нем?
ps. Увы, в программировании я не силен.

[Voltron]

Выше приложен пример моих шейпов, может дело в нем?

Да, проблема в исходных данных. У вас в слое часть объектов имеет тип LineString, а часть — MultiLinestring. Прогоните слой через Multipart to singleparts и будем вам счастье