Çokgenler arasındaki akışları (bağlantılar + değerler) aktarma

QGIS'de hücreler ve bir ek katman arasındaki hareketli verileri temsil eden iki şekil dosyası vardır, aşağıdaki resme bakın

Verileri tanımlayan taşıma:

• Çokgen "LayerA"(kırmızı anahat şeffaf kareler). Ayrıca, "LayerA"jeosentroidlerin konumu üzerinde görselleştirilen, hücreler içindeki hareketleri temsil eden daireler ile de ilgilidir .

  

• Polyline katmanı "Flows"(sarı / gri oklar), "LayerA"özelliklerin jeocentroidleri arasındaki bağlantılar aracılığıyla değerleri iletir

  

Hedef katman:

• Çokgen "LayerB"(koyu gri çerçeveli açık leylak özellikleri).

  

Ayrıca, ben zaten aktardıktan "FLUX"gelen hücrelerin içinde hareket değerleri ve "LayerA"içine "LayerB"çokgen, benim önceki soru bkz: QGIS çokgen arasındaki Kalıtsal değerler? . Bu kullanılarak yapıldı %ait $areahesaplanması.

Aktarılması anlamlı çözüm / yaklaşım da olabilir / verici / dönüştürme akış bağlantıları ile temsil "Flows"ilişkilerinden ve değerlerini "LayerA"ilişkilerine doğru "LayerB".

Bu bağlantıları çoklu hat olarak nasıl başarabilirim?

Ek olarak, yeni akışlar benzer bir stili miras alacaktır "Flows".

İstek üzerine, verilerin bir örneğini sağlayabilirim.

Akışlar özellikleri arasında değil var olacak "LayerA", fakat özellikleri arasında "LayerB" . Ana amaç, tablo / Köken-Hedef Matrisi olarak olası "FLUX"bağlantılar arasındaki özniteliği (yani / den) elde etmektir "LayerB".

Uyulması gereken bazı gereksinimler / kriterler vardır:

1. Aynı hücredeki özelliklerin parçaları (sarı ile seçilen) arasında akış bağlantısı yok

2. Parçaları farklı hücrelerde olsa bile aynı özellik arasında bağlantı yoktur.

3. Tamamen iki farklı hücre özelliği içindeyse , özelliklerin bölümleri "LayerB"( "Union"çıktıya dayalı ) arasında bağlantılar vardır"LayerA"

4."FLUX" İletilen yeni değer, aşağıdaki resimde gösterildiği gibi hesaplanacaktır.

Örneğin, iki hücre Ive IInerede "FLUX"olduğu arasında bir bağlantı vardır 100. Diğer değerleri varsayarsak, "NEW_FLUX"aralarında A've B''etrafında olacaktır 1.5625. 100sadece tek bir örnektir.

Referanslar:

• QGIS'de çokgenler arasında miras alınan değerler?
• QGIS 2'de akış haritası stillerini ayarlama
• Başka bir katmanın çokgeni içindeki çokgenin oransal alanını hesaplama
• Özgür ve Açık Kaynak GIS Ramblings | CBS'de hareket verileri: sorunlar ve fikirler
• QGIS 2'de Sanal Katman ile geometriyi koruma

Sanal Katmanlar ile teorik olarak mümkündür (şekil dosyaları ile işlem çok uzun olacaktır, ancak katmanlar bir Uzamsal Veritabanındaysa, bence çok daha hızlıdır).

İşte kod:

WITH inter_ab AS ( 
--create intersection between LayerA and LayerB 
SELECT LayerA.id || '_' || LayerB.FLAECHEID AS id, 
LayerA.id AS id_a, 
ST_AREA(LayerA.geometry) AS area_a, 
LayerB.FLAECHEID AS id_b, 
ST_INTERSECTION(LayerB.geometry, LayerA.geometry) AS geom 
FROM LayerA, LayerB 
WHERE ST_INTERSECTION(layerB.geometry, layerA.geometry) IS NOT NULL 
),

--calculation of the new flux value 
new_flux AS (SELECT t1.id_b AS origine, 
t2.id_b AS dest, 
SUM(Flows.flux * ST_AREA(t1.geom) / t1.area_a * ST_AREA(t2.geom) / t2.area_a) AS value  
FROM inter_ab t1, inter_ab t2, flows 
-- no connection between the same feature 
WHERE t1.id <> t2.id 
-- rule 1 
AND t1.id_a <> t2.id_a 
-- rule 2 
AND t1.id_b <> t2.id_b 
-- get flow data 
AND flows.origine = t1.id_a 
AND flows.dest = t2.id_a 
GROUP BY t1.id_b, t2.id_b
)

--create flows between original layerB features
SELECT new_flux.origine,
new_flux.dest,
new_flux.value AS flux,
make_line(ST_CENTROID(t3.geometry), ST_CENTROID(t4.geometry)) AS geom --ST_MakeLine under postGIS
FROM LayerB t3,
LayerB t4,
new_flux
WHERE t3.FLAECHEID = new_flux.origine
AND t4.FLAECHEID = new_flux.dest

Grafiksel çıktı

Sonuç manuel olarak test edildi. "FLUX"Değerler arasındaki fark ihmal edilebilir düzeydedir.

Son çıktı stilleri miras alır "Flow"ve

Birkaç veri ile test etmenizi öneririm ve büyük veri kümeleri için çok uzun sürerse, sorguları adım adım ( "inter_ab", "new_flux") yürütün ve sonucu kaydedin ve sonraki sorguyu yürütün.

Üç katman arasında birleştirme yapabilir ve ardından toplayabilirsiniz layerB. Sanal katmanlar kullanılabilir. Önemli verilerin katman layerAveya flowkatman içinde olup olmadığından emin değilim .. İşte (denenmemiş) bir olasılık:

SELECT b.id, b.geometry, sum(a.myVar)
FROM layerB b
LEFT JOIN flow f
   ON ST_Intersects(ST_EndPoint(f.geometry),b.geometry)
 JOIN layerA a
   ON ST_Intersects(ST_StartPoint(f.geometry),a.geometry)
GROUP BY b.id