Создание точек на указанном расстоянии вдоль линии с использованием QGIS？

Я хочу использовать QGIS для создания точек вдоль линии, но на указанном расстоянии. Я нашел два инструмента GRASS в SEXTANTE Toolbox ：

• v.to.points
• v.segment

Я попытался v.to.points , установив максимальное расстояние между точками в единицах карты до 100 и получил это:

Второе фото показывало близкий взгляд с верхнего. Результатом было то, что я хочу, но я хочу получить меньше очков , но даже я установил максимальное расстояние 1000 или 1000000, результат был таким же.

Я пытался использовать v.segment , но у меня не было файла, содержащего правила сегмента.

Подземье дало мне несколько советов, и я попробовал геометрию Densify с заданным интервалом . Сначала он дал мне шейп-файл ломаной линии, такой же, как я, и затем я использовал узлы извлечения, чтобы попытаться получить эти точки. Но результат был таким же, как я использовал v.to.points , независимо от того, какой интервал расстояния я установил.

Как я могу создать точки вдоль линии и на указанном расстоянии？

В консоли Python:

• Вы можете использовать модуль Shapely (как в Как создать равноудаленные точки в QGIS? ) с функцией

  point = line.interpolate (currentdistance)

• Новый Python API основной версии (1.9) имеет эквивалентную команду (как в Генерировании узлов (расстояний) в QGIS Натана Вудроу)

  point = geom.interpolate (currentdistance)

• или вы можете использовать векторную алгебру и направляющие косинусы (как в PyQGIS: геометрии, векторы, векторная алгебра или направляющие косинусы, примеры применения на французском языке)

1) создать класс слоя памяти (полный класс (без атрибутов) в crea_mem_layer.py )

class Create_vlayer(object):
    '''creation of a virtual layer''' 
     def __init__(self,nom,type):
         self.type=type
         self.name = nom
         self.layer =  QgsVectorLayer(self.type, self.name , "memory")
         self.pr =self.layer.dataProvider() 
     def create_point(self,geometry):
         # add point to the layer
         self.seg = QgsFeature()
         self.seg.setGeometry(QgsGeometry.fromPoint(geometry))
         self.pr.addFeatures([self.seg])
         self.layer.updateExtents()
     @property
     def display_layer(self):
         #end of layer and display layer 
         QgsMapLayerRegistry.instance().addMapLayers([self.couche])

2) функции векторной алгебры (из algèbre_vect_PyQGIS.py )

import math 
def mag(point):
    # magnitude of a vector
    return math.sqrt(point.x()**2 + point.y()**2)
def diff(point2, point1):
    # substraction betwen two vector
    return QgsPoint(point2.x()-point1.x(), point2.y() - point1.y())
def length(point1,point2):
    # with PyQGIS: sqrDist
    return math.sqrt(point1.sqrDist(point2))

3) направляющие косинусы

def dircos(point):
    cosa = point.x() / mag(point)
    cosb = point.y()/ mag(point)
    return cosa,cosb

4) технологическая линия или отрезки

def pairs(list):
    # list pairs iteration 
    for i in range(1, len(list)):
    yield list[i-1], list[i]

layer = qgis.utils.iface.activeLayer()
# interval between points
interval = 5 m
# create virtual layer
gen_pt  = Create_vlayer("mid5", "Point")

for elem in layer():
    line = elem.geometry()
    for seg_start, seg_end in pairs(line.asPolyline()):
       line_start = QgsPoint(seg_start)
       line_end = QgsPoint(seg_end)
       # mid point = vector coordinates [x2-x1,y2-y1]
       pointm =diff(line_end, line_start)
       # direction cosines of the segment
       cosa,cosb = dircos(pointm)
       # length of the segment
       lg = length(line_end, line_start)
       # generate and add points to the virtual layer 
       for i in range(interval,lg,interval):
           gen_pt.create_point(QgsPoint(line_start.x()  + (i * cosa), line_start.y() + (i*cosb)))

# display layer
gen_pt.display_layer

Результаты

полилиния, равноудаленная точка с Shapely или PyQGIS2, с направленными косинусами

Затем просто настройте интервал

Теперь есть плагин под названием QChainage, который делает именно то, что вы хотите. Вы указываете интервальное расстояние, и плагин создает точки вдоль вашей линии с указанным интервалом.

Я написал скрипт, который меняет инструмент геометрии Sextante Densify, чтобы принять определенное расстояние. Это называется Densify геометрией с заданным интервалом .

После запуска Densify, вы можете извлечь точки с помощью инструмента Извлечь узлы .

Вы можете получить его от Github и инструкции по установке находятся в моем блоге .

Если вы не хотите использовать скрипты Python, вы можете просто установить и использовать плагин «Профиль из линии» и игнорировать / удалить столбец значений растра. Цепочка будет тем, что вы установили для интервала выборки.

Ну, не уверен, что это в стабильном выпуске, но в наборе инструментов sextante на 1.9 alpha под Geoalgorithims-> Vector есть опция «Преобразовать линии в точки». Работает угощение, было бы здорово, если бы оно добавило поле для расстояния вдоль линии, хотя.

Интересно, что оно добавило поле NAME из моего векторного слоя.

Я пытался заставить работать этот Python Script Нейтана Вудроу, но я сосу с Python. А с кодом вообще похоже.

Я долго боролся с этой проблемой, но в итоге получил более или менее то, что искал, используя QChainage. Я хотел бы поделиться тем, что может помочь другим. При использовании инструмента растворения геообработки SAGA на линии QChainage дает неверные результаты. Использование другого инструмента растворения заставляет его работать нормально.