ポイントレイヤーとラインレイヤーの間の垂線の自動作成

QGISを使用していて、ポイントレイヤーから別のラインレイヤーのラインに多数の垂直線を作成できるスクリプトまたはプラグインを探しています。

これまで、MMQGISのハブ距離関数（ラインをポイントに変換し、ポイントを最も近いハブに接続する）とQGIS地理アルゴリズムのそれぞれのツールを使用してみました。どちらもうまくいきませんでした。どちらも2時間以上かかり、すべてのレイヤーの線、またはポイントに垂直または接続されていない線を作成します。

写真では、プロジェクトの現在のステータスを確認できます。垂直線は、ポイントから最も近い線まで延びる必要があります。最後に、ポイントと国境を結ぶ線との交点を使用して、2ポリゴンの深さである4辺ポリゴンのバッファを作成します。これを行う簡単な方法がある場合に備えて、これについて触れます。

垂線を作成する方法についていくつかの記事があることは知っていますが、どちらも私の問題を解決しませんでした。

次のスクリプトは、ポイントレイヤーとラインレイヤーの間の垂線の作成を自動化しました。作成された垂直セグメント（メモリレイヤーのフィーチャ）は、ポイントからラインレイヤーの最も近いフィーチャまで続きます。

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)

私は質問で提示されたのと非常に似た状況でそれを試しました：

QGISのPythonコンソールでコードを実行した後、コードが取得されました。

以下は、Python 3.x（または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)