如何在属性表中添加方向和距离？

任何可以提供帮助的人。我只想在多义线/线数据中添加方向（轴承：即N 25 35 E）和距离（长度：125米）作为我的新字段。是否有插件来生成这些字段？我试图在行数据中使用“导出/添加几何列”，但仅添加了“长度”值。

您可以在QGIS的字段计算器中计算方位角。这适用于小距离（几百公里）内的UTM（公制）坐标，但是对于大距离或十进制度数，则需要更复杂的方法。

打开您的线图层的属性表，切换编辑，然后单击“字段计算器”按钮以打开对话框：

以1或2的精度创建一个新的字段作为十进制。

将此代码粘贴到“表达式”框中，然后单击“确定”： (atan((xat(-1)-xat(0))/(yat(-1)-yat(0)))) * 180/3.14159 + (180 *(((yat(-1)-yat(0)) < 0) + (((xat(-1)-xat(0)) < 0 AND (yat(-1) - yat(0)) >0)*2)))

第一部分计算x和y差的反正切，并将其转换为度数（180 / pi）。第二部分将180或360添加到所得图形中，以提供0-360°的方位角。

您不需要插件。一切都在PyQGIS的QgsPoint类中

如果使用Python控制台中的Python内置函数dir（）检查QGIS点类的内容。

dir(point])
['__class__', '__delattr__', '__dict__', '__doc__', '__eq__', '__format__', '__ge__', '__getattribute__'
, '__getitem__', '__gt__', '__hash__', '__init__', '__le__', '__len__', '__lt__', '__module__', 
'__ne__', '__new__', '__reduce__', '__reduce_ex__', '__repr__', '__setattr__', '__sizeof__', 
'__str__', '__subclasshook__', '__weakref__', 'azimuth', 
'multiply', 'set', 'setX', 'setY', 'sqrDist', 'sqrDistToSegment', 'toDegreesMinutesSeconds', 'toString', 'wellKnownText', 'x', 'y']

您可以看到有方位角和sqrDist函数，并且经过几次尝试：

- xy[0].azimuth(xy[1]) or xy[1].azimuth(xy[0]) gives the azimuth direction between two points(in degrees, +/- 180°)
- xy[0].sqrDist(xy[1]) give the square distance between two points (in the unit of the project)

问题

因此，在Python控制台中

def select_all(layer):
     layer.select([])
     layer.setSelectedFeatures([obj.id() for obj in layer])

myline = qgis.utils.iface.activeLayer()
select_all(myline)
for elem in myline.selectedFeatures():
      xy = elem.geometry().asPolyline()

现在xy包含该行的所有节点（点）

# first point
print "x=%2d y=%2d" % (xy[0].x(),xy[0].y())
x=112935 y=117784
# and others...

使用线的所有节点：

1）方位角i至点i + 1（+/- 180°）（一条线的节点）

for i in range(len(xy)-1):
     print "x=%2d y=%2d azim=%6.1f azim2=%6.1f" % (xy[i].x(), xy[i].y(), xy[i].azimuth(xy[i+1]), xy[i+1].azimuth(xy[i]))

x=112935 y=117784 azim= 168.4 azim2= -11.6
x=113032 y=117312 azim=-167.5 azim2=  12.5
x=112926 y=116835 azim= 177.3 azim2=  -2.7
x=112943 y=116472 azim= 145.1 azim2= -34.9
[...]

2）点i与点i + 1之间的欧式距离

for i in range(len(xy)-1):
     print "x=%2d y=%2d dist=%6.1f" % (xy[i].x(), xy[i].y(), xy[i].sqrDist(xy[i+1]))

x=112935 y=117784 dist=232533.9
x=113032 y=117311 dist=238243.6
x=112926 y=116835 dist=131839.8
x=112943 y=116472 dist=209268.1
[...]

之后，将这些值添加到属性表并不是很困难。

我使用这项技术通过matplotlib和Script Runner插件来分析线条（地质）

@Simbamangu提供的解决方案非常有效，但不能涵盖所有情况。例如，将公式与水平位移一起应用将使结果为NULL，因此必须在QGIS的字段计算器中使用此公式

case
when yat(-1)-yat(0) < 0 or yat(-1)-yat(0) > 0 then 
(atan((xat(-1)-xat(0))/(yat(-1)-yat(0)))) * 180/3.14159 + 
(180 *
(((yat(-1)-yat(0)) < 0) + 
(((xat(-1)-xat(0)) < 0 AND (yat(-1) - yat(0)) >0)*2)
))
when ((yat(-1)-yat(0)) = 0 and (xat(-1) - xat(0)) >0) then 90
when ((yat(-1)-yat(0)) = 0 and (xat(-1) - xat(0)) <0) then 270
end