QGIS kullanarak birden çok satırdan “ortalama” satır oluşturma

Çeşitli GPS hassasiyeti ile benzer yolu (bu durumda dağ izleri) temsil eden birden fazla çizgi özelliğine (aynı zamanda birden fazla noktaya sahip ...) sahibim, bu da birbirine yakın birçok çizgiye neden oluyor, ancak mükemmel bir şekilde çakışmıyor.

Bu projenin amacı için, hepsinden bir "ortalama" çizgi hesaplamanın ve patikanın en olası pozisyonunun bir sonucunu üretmenin bir yolunu arıyorum.

QGIS veya başka bir araç kullanarak nasıl ilerlersiniz (OGR'yi düşündüm ...)?

Test verisi:

• QGIS 2.18.16, ÇİM CBS 7
• 4 GPS Parkuru
• 1x1km'lik bir ızgara içinde

BEN.)

QGIS Eklentisi Locate Points Along Lines( https://plugins.qgis.org/plugins/LocatePoints/ ) ile GPS Parkurlarınızda puan oluşturun . Örneklemimde 5 metrelik bir aralık kullandım.

II.)

Bir oluşturma Concave Hullile Processing > Toolbox > QGIS geoalgorithms >Vector geometry tools > Concave hull. Örneğimde bir eşik kullandım 0.1. Eşik düşükse, çıkış poligonunda delikler olabilir.

III.)

Şimdi, iskelet algoritması yardımıyla "ortalama" çizgiyi hesaplayabilirsiniz. İskeleti arayın Processing Toolbox. v.voronoi.skeletonAracını kullanın GRASS GIS 7 commands.

Isı haritası yaklaşımı:

İşleme aşamalarının hesaplanma süresi nedeniyle yaklaşım daha fazla zaman alır. Belki daha genel bir çözüme yaklaşmak bir fikir olarak görülebilir.

Test verisi:

• QGIS 2.18.16, ÇİM CBS 7
• 4 GPS Parkuru
• 1x1km'lik bir ızgara içinde

BEN.)

QGIS Eklentisi ile GPS Parkurlarınız boyunca noktalar Boyunca Nokta Bulun ( https://plugins.qgis.org/plugins/LocatePoints/ ). Isı haritası yaklaşımı için 2m aralık kullandım .

II.)

QGIS ısı haritası eklentisi ile bir ısı haritası oluşturun. 40 m yarıçapı kullandım. Çıkış rasterinde delik kalmayana kadar yarıçapı artırıyorum. Bunu farklı yarıçap değerleriyle denemelisiniz.

III.)

DÜZENLENMİŞ Isı haritası tarama değerinin tam olarak tutulması gerekmez.

Şimdi rasterini en fazla puanın yoğunlaştığı "beyaz" alanlara inceltmek istiyorum. Bunun için çıkış rasterini yeniden hesaplıyorum. Çıkış rasterinin min / maks değerleri 0ve 89.7935. Sadece yukarıdaki değerleri kullanıyorum 44. Bu yüzden bir "kural" kullandım. Maksimum değeri yuvarlayın ve ikiye bölün. Bu değeri başka bir zaman aşağı yuvarlayın. 89/2 = 44,5-> 44. Ben OSGeo4W Shell kullandım: gdal_calc -A heatmap.tif --calc="A>=44" --NoDataValue=0 --outfile=heatmap_44_NoData.tif.

IV.)

REDAKTE

a) Yeniden hesaplanan ısı haritasını ile poligonize edin Raster > Conversion > Polygonize ...

b) Çokgeni basitleştirin Vector > Geometry Tools > Simplify geometries. Bir hoşgörü kullandım 2. Daha basit bir çokgen iskeletlerin işlem süresini azaltır.

c) İskeletleri hesaplayın: İşleme Araç Kutusunda iskeleti arayın. v.voronoi.skeletonGRASS GIS 7 komutlarındaki aracı kullanın .

Ortaya çıkan çizginin, ilk cevabımdakinden daha fazla yolun en olası konumunu temsil ettiğini görebilirsiniz. Özellikle kuzeydeki kıvrım için ortalama çizgi birbirine daha yakın olan üç pisti takip eder. Doğudaki viraj için de aynı şey geçerli.

Yaklaşımın avantajları:

• sadece QGIS kullanarak makul iyi sonuçlar

Dezavantajları:

• büyük veri kümeleri için işlem süresi
• parametreleri bir priori denemelisiniz (ısı haritası yarıçapı, min / maks değerler)
• işlem adımlarını otomatikleştirmek zor
• Dar virajlar / eğriler ve gerçekten çizginin dışına çıkan parçalar için test edilmedi

Birisi işleme adımlarını optimize edebilirse, hoş geldiniz!