Sắp xếp nhiều điểm để xếp hàng trong QGIS?

Tôi muốn căn chỉnh nhiều điểm cho một dòng hoặc các dòng trong một lớp, sử dụng dung sai hoặc bộ đệm được chỉ định xung quanh các đối tượng dòng. Vui lòng tham khảo bản phác thảo mẫu đính kèm.

Vì lợi ích của ví dụ này, các điểm gần đường thẳng nhất trong hình TRƯỚC là trong 5 đơn vị bản đồ của đường, trong khi các điểm ngoài cùng cách xa hơn 10 đơn vị bản đồ. Tôi muốn chụp các điểm gần nhất vào đường gần nhất, sử dụng dung sai 5 đơn vị bản đồ để đạt được kết quả trong ảnh SAU.

Có một công cụ tích hợp để thực hiện việc này trong phiên bản QGIS 3.0 (chưa được phát hành). Bạn có thể nhận được một ảnh chụp nhanh hàng đêm từ trang web của QGIS để kiểm tra trước điều này.

Để làm điều này:

1. Chạy thuật toán xử lý "Snap geometries to layer"
2. Chọn lớp điểm của bạn làm "lớp đầu vào"
3. Chọn lớp dòng là "lớp tham chiếu"
4. Nhập dung sai phù hợp (khoảng cách tối đa để di chuyển điểm trong khi chụp)
5. Thay đổi hành vi thành "Thích điểm gần nhất"

Đây là kết quả, hiển thị các điểm ban đầu là "x" và các điểm được ngắt là các chấm màu lục. Tôi đã sử dụng một dung sai ở đây để chỉ một số điểm đầu vào bị ngắt.

Điều này có thể được cung cấp với PyQGIS . Đối với tình huống tiếp theo:

mã sau, xem xét dung sai của 5 đơn vị bản đồ, đã được chạy tại Bảng điều khiển Python của QGIS:

from math import sqrt

registry = QgsMapLayerRegistry.instance()

points = registry.mapLayersByName('points')
line = registry.mapLayersByName('line')

feat_points = [ feat for feat in points[0].getFeatures() ]
feat_line = line[0].getFeatures().next()

new_points = []

for feat in feat_points:
    pt = feat.geometry().asPoint()
    sqrdist, point, vertex = feat_line.geometry().closestSegmentWithContext(pt)
    if sqrt(sqrdist) <= 5:
        new_points.append(point)
    else:
        new_points.append(pt)

epsg = points[0].crs().postgisSrid()

uri = "Point?crs=epsg:" + str(epsg) + "&field=id:integer""&index=yes"

mem_layer = QgsVectorLayer(uri,
                           'new_points',
                           'memory')

prov = mem_layer.dataProvider()

feats = [ QgsFeature() for i in range(len(new_points)) ]

for i, feat in enumerate(feats):
    feat.setAttributes([i])
    feat.setGeometry(QgsGeometry.fromPoint(new_points[i]))

prov.addFeatures(feats)

QgsMapLayerRegistry.instance().addMapLayer(mem_layer)

Nó được tạo ra một lớp bộ nhớ trong đó các điểm được chụp theo dung sai được coi là trước đây của 5 đơn vị bản đồ và đường đi vuông góc với đường thẳng.

Bạn cũng có thể thực hiện việc này trong Máy tính trường với plugin refFifts. Bạn có thể sử dụng Máy tính trường để cập nhật hình dạng lớp cũng như các trường. refFifts cung cấp cho bạn hàm "geomdistance" để tìm dòng gần nhất trong một khoảng cách nhất định (hoặc "geomgầnest" nếu bạn không muốn có ngưỡng) và sẽ trả về một thuộc tính hoặc hình học, và hàm "near_point" sẽ tìm thấy gần nhất điểm trên một hình học nhất định. Chuỗi chúng lại với nhau như vậy để tính toán hình học mới cho lớp điểm của bạn:

closest_point(geom_from_wkt(geomdistance('snap_lines','$geometry',10)) , $geometry)

Thay vì trực tiếp cập nhật hình học, bạn có thể tính toán một trường với hình dạng được chụp thay thế. Tôi lưu trữ nhiều hình học để chụp các điểm cống vào các lớp luồng khác nhau và tôi có thể dễ dàng cập nhật hình dạng điểm trong Máy tính Trường tùy thuộc vào dòng dòng tôi cần sử dụng.

Có một số hạn chế với điều này, cả hai lớp cần phải là cùng một CRS và chức năng phân tích địa lý sẽ gây ra lỗi nếu bạn có hơn 100.000 điểm nhưng bạn có thể thay đổi giới hạn này nếu bạn chỉnh sửa tệp plugin refFifts.