การสร้างจุดในระยะทางที่กำหนดตามเส้นโดยใช้ QGIS？

ฉันต้องการใช้จุดสร้าง QGIS ตามแนวเส้น แต่ตามระยะทางที่กำหนด ฉันพบเครื่องมือ GRASS สองตัวในกล่องเครื่องมือ SEXTANTE：

• v.to.points
• v.segment

ฉันลองv.to.pointsตั้งค่าระยะห่างสูงสุดระหว่างจุดในหน่วยแผนที่เป็น 100 และได้สิ่งนี้

ภาพถ่ายที่สองแสดงให้เห็นอย่างใกล้ชิดจากภาพบน ผลที่ได้คือสิ่งที่ฉันต้องการ แต่ฉันต้องการได้คะแนนน้อยลงแต่แม้ฉันจะตั้งระยะทางสูงสุดเป็น 1,000 หรือ 1000000 ผลลัพธ์ก็เหมือนกัน

ฉันพยายามใช้v.segmentแต่ฉันไม่มีไฟล์ที่มีกฎของเซ็กเมนต์

Underdark ให้คำแนะนำบางอย่างและฉันพยายามรูปทรงเรขาคณิตที่บดอัดให้ช่วงเวลา ครั้งแรกมันทำให้ฉันมีรูปร่างไฟล์ polyline เหมือนกับอินพุตของฉันและจากนั้นฉันใช้แยกโหนดเพื่อพยายามหาจุดเหล่านั้น แต่ผลลัพธ์เหมือนกันกับที่ฉันใช้v.to.pointsไม่ว่าฉันจะตั้งช่วงเวลาไว้ที่ใด

ฉันจะสร้างคะแนนตามเส้นได้อย่างไรและในระยะทางที่กำหนด？

ในคอนโซล Python:

• คุณสามารถใช้ โมดูลShapely (เช่นเดียวกับในวิธีการสร้างจุดที่เท่ากันใน QGIS? ) กับฟังก์ชั่น

  จุด = line.interpolate (ปัจจุบันระยะทาง)

• Python API ใหม่ของเวอร์ชันต้นแบบ (1.9) มีคำสั่งที่เทียบเท่า (เช่นในการสร้างโหนด chainage (ระยะทาง) ใน QGISของ Nathan Woodrow)

  จุด = geom.interpolate (ระยะทางปัจจุบัน)

• หรือคุณสามารถใช้พีชคณิตเวกเตอร์และทิศทางโคไซน์ (เช่นในPyQGIS: รูปทรงเรขาคณิต, เวกเตอร์, พีชคณิตเวกเตอร์หรือทิศทางโคไซน์, ตัวอย่างการใช้งาน , ในภาษาฝรั่งเศส)

1) สร้างคลาสเลเยอร์หน่วยความจำ (คลาสที่สมบูรณ์ (ไม่มีแอตทริบิวต์) ที่crea_mem_layer.py )

class Create_vlayer(object):
    '''creation of a virtual layer''' 
     def __init__(self,nom,type):
         self.type=type
         self.name = nom
         self.layer =  QgsVectorLayer(self.type, self.name , "memory")
         self.pr =self.layer.dataProvider() 
     def create_point(self,geometry):
         # add point to the layer
         self.seg = QgsFeature()
         self.seg.setGeometry(QgsGeometry.fromPoint(geometry))
         self.pr.addFeatures([self.seg])
         self.layer.updateExtents()
     @property
     def display_layer(self):
         #end of layer and display layer 
         QgsMapLayerRegistry.instance().addMapLayers([self.couche])

2) ฟังก์ชันพีชคณิตเวกเตอร์ (จากalgèbre_vect_PyQGIS.py )

import math 
def mag(point):
    # magnitude of a vector
    return math.sqrt(point.x()**2 + point.y()**2)
def diff(point2, point1):
    # substraction betwen two vector
    return QgsPoint(point2.x()-point1.x(), point2.y() - point1.y())
def length(point1,point2):
    # with PyQGIS: sqrDist
    return math.sqrt(point1.sqrDist(point2))

3) ทิศทางโคไซน์

def dircos(point):
    cosa = point.x() / mag(point)
    cosb = point.y()/ mag(point)
    return cosa,cosb

4) สายการผลิตหรือส่วนของเส้น

def pairs(list):
    # list pairs iteration 
    for i in range(1, len(list)):
    yield list[i-1], list[i]

layer = qgis.utils.iface.activeLayer()
# interval between points
interval = 5 m
# create virtual layer
gen_pt  = Create_vlayer("mid5", "Point")

for elem in layer():
    line = elem.geometry()
    for seg_start, seg_end in pairs(line.asPolyline()):
       line_start = QgsPoint(seg_start)
       line_end = QgsPoint(seg_end)
       # mid point = vector coordinates [x2-x1,y2-y1]
       pointm =diff(line_end, line_start)
       # direction cosines of the segment
       cosa,cosb = dircos(pointm)
       # length of the segment
       lg = length(line_end, line_start)
       # generate and add points to the virtual layer 
       for i in range(interval,lg,interval):
           gen_pt.create_point(QgsPoint(line_start.x()  + (i * cosa), line_start.y() + (i*cosb)))

# display layer
gen_pt.display_layer

ผล

polyline จุดกึ่งกลางที่มี Shapely หรือ PyQGIS2 ที่มีทิศทาง cosines

จากนั้นเพียงแค่ปรับช่วงเวลา

ขณะนี้มีปลั๊กอินที่เรียกว่า QChainage ที่ทำสิ่งที่คุณต้องการ คุณระบุระยะห่างของช่วงเวลาและปลั๊กอินสร้างจุดตามเส้นของคุณในช่วงเวลาที่ระบุ

ฉันเขียนสคริปต์ที่เปลี่ยนเครื่องมือรูปทรงเรขาคณิต Sextante Densify เพื่อยอมรับระยะทางที่แน่นอน มันเรียกว่ารูปทรงเรขาคณิตที่บดอัดให้ช่วงเวลา

หลังจากรัน Densify คุณสามารถแยกคะแนนโดยใช้เครื่องมือExtract nodes

คุณจะได้รับจากGithubและติดตั้งคำแนะนำบนของบล็อก

หากคุณไม่ต้องการใช้สคริปต์ python คุณสามารถติดตั้งและใช้ปลั๊กอิน "Profile from line" และละเว้น / ลบคอลัมน์ค่าแรสเตอร์ เชนจะเป็นสิ่งที่คุณกำหนดช่วงเวลาการสุ่มตัวอย่างให้เป็น

ไม่แน่ใจว่านี่เป็นรุ่นเสถียรหรือไม่ในกล่องเครื่องมือ sextante ที่ 1.9 alpha ภายใต้ Geoalgorithims-> Vector เป็นตัวเลือก "แปลงสายเป็นคะแนน" ใช้การรักษาได้ดีถ้ามันเพิ่มระยะทางตามแนวนอนเช่นกัน

ที่น่าสนใจก็เพิ่มฟิลด์ NAME จากเลเยอร์เวกเตอร์ของฉัน

ฉันพยายามทำให้ Python Script ทำงานโดย Nathan Woodrow แต่ฉันดูดด้วย python และด้วยรหัสโดยทั่วไปดูเหมือนว่า

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