Tự động tạo các đường vuông góc giữa một lớp điểm và một lớp đường

Tôi đang sử dụng QGIS và tôi đang tìm kiếm một tập lệnh hoặc plugin có thể tạo ra một số lượng lớn các đường vuông góc từ một lớp điểm đến các dòng trong một lớp dòng riêng biệt.

Cho đến nay, tôi đã thử sử dụng chức năng Khoảng cách Hub trong MMQGIS (biến các đường thành điểm và sau đó kết nối các điểm với trung tâm gần nhất) và công cụ tương ứng của thuật toán địa lý QGIS. Không làm việc. Cả hai mất hơn 2 giờ và tạo một trong hai dòng tất cả các lớp hoặc các đường không vuông góc hoặc kết nối với các điểm.

Trong hình, bạn có thể thấy trạng thái hiện tại của dự án. Các đường vuông góc sẽ chạy từ các điểm đến đường gần nhất. Cuối cùng, tôi muốn sử dụng các điểm giao nhau với một đường thẳng giữa các điểm và biên giới quốc gia để tạo ra một bộ đệm của đa giác 4 ​​mặt sâu hai đa giác. Tôi đề cập đến điều này trong trường hợp có một cách dễ dàng hơn để làm điều này.

Tôi biết rằng có một số bài viết về cách tạo các đường vuông góc, nhưng cả hai đều không giải quyết được vấn đề của tôi.

Kịch bản tiếp theo tự động tạo các đường vuông góc giữa một lớp điểm và một lớp đường. Các phân đoạn vuông góc (tính năng của lớp bộ nhớ) được tạo chạy từ các điểm đến tính năng gần nhất của lớp đường.

mapcanvas = iface.mapCanvas()

layers = mapcanvas.layers()

p_lyr = layers[0]
l_lyr = layers[1]

epsg = p_lyr.crs().postgisSrid()

uri = "LineString?crs=epsg:" + str(epsg) + "&field=id:integer""&field=distance:double(20,2)&index=yes"

dist = QgsVectorLayer(uri, 
                      'dist', 
                      'memory')

QgsMapLayerRegistry.instance().addMapLayer(dist)

prov = dist.dataProvider()

lines_features = [ line_feature for line_feature in l_lyr.getFeatures() ] 
points_features = [ point_feature for point_feature in p_lyr.getFeatures() ]

feats = []

for p in points_features:

    minDistPoint = min([l.geometry().closestSegmentWithContext( p.geometry().asPoint() ) for l in lines_features])[1]
    feat = QgsFeature()
    feat.setGeometry(QgsGeometry.fromPolyline([p.geometry().asPoint(), minDistPoint]))
    feat.setAttributes([points_features.index(p),feat.geometry().length()])
    feats.append(feat)

prov.addFeatures(feats)

Tôi đã thử nó với một tình huống rất giống với câu hỏi:

Sau khi chạy mã tại Bảng điều khiển Python của QGIS, nó đã được lấy:

Đây là mã giống như trong câu trả lời được chấp nhận chỉ được điều chỉnh để hoạt động với Python 3.x (hoặc QGIS v3.x):

from qgis.core import QgsProject

mapcanvas = iface.mapCanvas()

layers = mapcanvas.layers()

p_lyr = layers[0]
l_lyr = layers[1]

epsg = p_lyr.crs().postgisSrid()

uri = "LineString?crs=epsg:" + str(epsg) + "&field=id:integer""&field=distance:double(20,2)&index=yes"

dist = QgsVectorLayer(uri,
                      'dist',
                      'memory')

QgsProject.instance().addMapLayer(dist)

prov = dist.dataProvider()

lines_features = [ line_feature for line_feature in l_lyr.getFeatures() ]
points_features = [ point_feature for point_feature in p_lyr.getFeatures() ]

feats = []

for p in points_features:

    minDistPoint = min([l.geometry().closestSegmentWithContext( p.geometry().asPoint() ) for l in lines_features])[1]
    feat = QgsFeature()
    feat.setGeometry(QgsGeometry.fromPolylineXY([p.geometry().asPoint(), minDistPoint]))
    feat.setAttributes([points_features.index(p),feat.geometry().length()])
    feats.append(feat)

prov.addFeatures(feats)