使用QGIS在沿线的指定距离内创建点？

我想使用QGIS沿直线创建点，但要指定距离。我在SEXTANTE工具箱中找到了两个GRASS工具：

• v.to.points
• 段

我尝试了v.to.points，将地图单位中的点之间的最大距离设置为100，得到的是：

第二张照片显示了上部的近距离观察。结果是我想要的，但是我想获得更少的分数，但是即使我将最大距离设置为1000或1000000，结果也是一样。

我尝试使用v.segment，但是我没有包含细分规则的文件。

Underdark给了我一些建议，并且我尝试给定间隔一段时间使Densify几何形状成为可能。首先，它给了我与输入相同的折线shapefile，然后使用提取节点尝试获取这些点。但是无论设置距离间隔如何，结果都与使用v.to.points获得的结果相同。

如何沿直线并在指定距离内创建点？

在Python控制台中：

• 您可以将 Shapely模块（如“ 如何在QGIS中创建等距点”中所述）一起使用

  点= line.interpolate（currentdistance）

• 主版本（1.9）的新Python API具有等效的命令（如Nathan Woodrow的QGIS中的生成里程（距离）节点）

  点= 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后，可以使用“ 提取节点”工具提取点。

您可以从Github获得它，安装说明在我的博客上。

如果您不想使用python脚本，则可以简单地安装并使用“来自行的配置文件”插件，而忽略/删除栅格值列。桩号将是您设置的采样间隔。

不太确定这是否在稳定版本中，但是在Geoalgorithims-> Vector下1.9 alpha的sextante工具箱中是“将线转换为点”选项。很有效，但如果它也沿着直线增加了一个距离字段，那就太好了。

有趣的是，确实从我的矢量层添加了NAME字段。

我试过让Nathan Woodrow来工作的Python脚本，但是我对python很烂。而且一般来说，代码似乎。

我花了很长时间在这个问题上苦苦挣扎，但最终还是得到了我使用QChainage寻找的东西。我想分享的是可以帮助他人的东西。在生产线上使用SAGA地理处理分解工具时，QChainage会产生错误的结果。使用其他溶解工具可以使其正常工作。