การจัดแนวหลายจุดให้เป็นบรรทัดใน QGIS?

ฉันต้องการจัดตำแหน่งหลายจุดให้เป็นหนึ่งบรรทัดหรือหลายบรรทัดภายในเลเยอร์โดยใช้เกณฑ์ความคลาดเคลื่อนหรือบัฟเฟอร์รอบวัตถุเส้นตรง โปรดดูตัวอย่างร่างที่แนบมา

เพื่อให้เป็นไปตามตัวอย่างนี้จุดที่อยู่ใกล้กับเส้นมากที่สุดในภาพก่อนหน้านั้นจะอยู่ในระยะ 5 หน่วยแผนที่ของเส้นในขณะที่จุดนอกสุดอยู่ห่างจากแผนที่มากกว่า 10 หน่วย ฉันต้องการถ่ายจุดที่ใกล้ที่สุดลงบนเส้นที่ใกล้ที่สุดโดยใช้ความอดทน 5 หน่วยแผนที่เพื่อให้ได้ผลลัพธ์ในภาพหลัง

มีเครื่องมือในตัวให้ทำเช่นนี้ในเวอร์ชัน QGIS 3.0 (ไม่ได้เผยแพร่) คุณสามารถรับภาพรวมจากเว็บไซต์ QGIS ทุกคืนเพื่อทดสอบล่วงหน้า

เพื่อทำสิ่งนี้:

1. เรียกใช้อัลกอริทึมการประมวลผล "Snap geometries to layer"
2. เลือกเลเยอร์คะแนนของคุณเป็น "เลเยอร์อินพุต"
3. เลือกเลเยอร์บรรทัดเป็น "เลเยอร์อ้างอิง"
4. ป้อนความคลาดเคลื่อนที่เหมาะสม (ระยะทางสูงสุดในการเคลื่อนย้ายจุดขณะถ่ายภาพ)
5. เปลี่ยนพฤติกรรมเป็น "ชอบจุดที่ใกล้ที่สุด"

นี่คือผลลัพธ์โดยแสดงจุดเดิมเป็น "x" และจุดที่จัดชิดเป็นจุดสีเขียว ฉันใช้ความอดทนที่นี่เพื่อให้มีเพียงบางจุดอินพุตเท่านั้น

นี้สามารถ afforded กับPyQGIS สำหรับสถานการณ์ต่อไป:

รหัสต่อไปนี้เมื่อพิจารณาถึงความทนทานต่อแผนที่ 5 หน่วยถูกเรียกใช้ที่ Python Console ของ 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)

มันผลิตเลเยอร์หน่วยความจำที่มีการหักคะแนนตามความคลาดเคลื่อนที่พิจารณาก่อนหน้านี้ของ 5 หน่วยแผนที่และเส้นทางตั้งฉากกับเส้น

คุณยังสามารถทำสิ่งนี้ได้ใน Field Calculator ด้วยปลั๊กอิน refFunctions คุณสามารถใช้เครื่องคำนวณภาคสนามเพื่ออัปเดตรูปทรงเรขาคณิตของเลเยอร์และฟิลด์ refFunctions ให้ฟังก์ชัน "geomdistance" เพื่อค้นหาบรรทัดที่ใกล้ที่สุดภายในระยะทางที่กำหนด (หรือ "geomnearest" หากคุณไม่ต้องการขีด จำกัด ) และจะส่งคืนแอตทริบิวต์หรือเรขาคณิตและฟังก์ชัน "closest_point" จะค้นหาสิ่งที่ใกล้เคียงที่สุด ชี้ไปที่รูปทรงเรขาคณิตที่กำหนด สตริงพวกมันเข้าด้วยกันอย่างนั้นเพื่อคำนวณรูปทรงเรขาคณิตใหม่สำหรับเลเยอร์จุดของคุณ:

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

แทนที่จะอัพเดทรูปทรงเรขาคณิตโดยตรงคุณสามารถคำนวณเขตข้อมูลด้วยรูปทรงเรขาคณิตที่จัดชิดแทน ฉันเก็บรูปทรงหลายรูปแบบเพื่อจับจุดระบายน้ำไปยังเลเยอร์สตรีมที่แตกต่างกันและฉันสามารถอัปเดตรูปทรงเรขาคณิตของจุดในฟิลด์เครื่องคิดเลขได้อย่างง่ายดายโดยขึ้นอยู่กับว่าฉันต้องการใช้สายข้อมูลใด

มีข้อ จำกัด บางประการเช่นนี้เลเยอร์ทั้งคู่ต้องเป็น CRS เดียวกันและฟังก์ชัน geomdistance จะให้ข้อผิดพลาดหากคุณมีมากกว่า 100,000 คะแนน แต่คุณสามารถเปลี่ยนข้อ จำกัด นี้ได้หากคุณแก้ไขไฟล์ปลั๊กอิน refFunctions