Làm cách nào để thêm Hướng và Khoảng cách vào bảng thuộc tính?

bất cứ ai có thể giúp đỡ. tôi chỉ muốn thêm hướng (Mang: tức là N 25 35 E) và Khoảng cách (Chiều dài: 125 mét) làm trường mới của tôi trong dữ liệu đa tuyến / dòng. Có một trình cắm để tạo các trường này không? tôi đã cố sử dụng "xuất / thêm cột hình học" trong dữ liệu dòng của mình, nhưng chỉ thêm giá trị "Độ dài".

Bạn có thể tính toán mang trong Máy tính Trường trong QGIS. Điều này hoạt động theo tọa độ UTM (số liệu) trên khoảng cách nhỏ (hàng trăm km), nhưng cần một cái gì đó tinh vi hơn cho khoảng cách lớn hoặc cho độ thập phân.

Mở bảng thuộc tính cho lớp dòng của bạn, chuyển đổi chỉnh sửa và nhấp vào nút Máy tính trường để mở hộp thoại:

Tạo một trường mới dưới dạng thập phân với độ chính xác 1 hoặc 2.

Dán mã này vào hộp "Biểu thức" và nhấp vào "OK": (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)))

Phần đầu tiên tính toán tiếp tuyến nghịch đảo của sự khác biệt x và y và chuyển đổi nó thành độ (180 / pi). Phần thứ hai thêm 180 hoặc 360 vào con số kết quả để mang ổ đỡ từ 0-360 °.

Bạn không cần một plugin. Tất cả mọi thứ đều nằm trong lớp QssPoint của PyQGIS

Nếu bạn kiểm tra nội dung của lớp điểm QGIS với hàm dir () tích hợp sẵn trong Python Console.

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']

Bạn có thể thấy có các hàm phương vị và sqrDist và sau một vài lần thử:

- 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)

Vấn đề

Vì vậy, trong bảng điều khiển 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()

bây giờ xy chứa tất cả các nút (điểm) của dòng

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

Sử dụng tất cả các điểm nút của dòng:

1) góc phương vị i đến điểm i + 1 (+/- 180 °) (các nút của một dòng)

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) khoảng cách euclide giữa điểm i và điểm 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
[...]

Sau đó, không khó để thêm các giá trị này vào bảng thuộc tính.

Tôi sử dụng kỹ thuật này để phân tích các dòng dõi (địa chất) với matplotlib và plugin Script Runner

Giải pháp được cung cấp bởi @Simbamangu khá hiệu quả nhưng không bao gồm tất cả các trường hợp. Ví dụ: áp dụng công thức với chuyển vị ngang sẽ NULL kết quả, do đó bạn phải sử dụng công thức này trong Máy tính trường của 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