Mentransfer aliran (koneksi + nilai) antara poligon

Dalam QGIS ada dua shapefile yang mewakili data bergerak antara sel dan satu lapisan tambahan, lihat gambar di bawah ini

Memindahkan data yang didefinisikan oleh:

• Poligon "LayerA"(kotak transparan dengan garis merah). Selain itu juga berhubungan dengan lingkaran yang mewakili gerakan di dalam sel, divisualisasikan pada posisi "LayerA"geocentroid.

  

• Lapisan polyline "Flows"(panah kuning / abu-abu), menyampaikan nilai melalui koneksi antara "LayerA"fitur geocentroid

  

Lapisan target:

• Polygon "LayerB"(fitur ungu muda dengan garis abu-abu gelap).

  

Selain itu, saya telah mentransfer "FLUX"dan memindahkan nilai dalam sel dari "LayerA"ke dalam "LayerB"poligon, lihat pertanyaan saya sebelumnya: Nilai yang diwariskan antara poligon di QGIS? . Itu dilakukan menggunakan% dari $areaperhitungan.

Mungkin ada solusi / pendekatan yang berarti untuk mentransfer / mentransmisikan / mengubah koneksi aliran yang diwakili oleh"Flows" dan nilai-nilainya dari hubungan "LayerA"ke dalam hubungan "LayerB".

Bagaimana saya bisa mencapai koneksi tersebut sebagai polyline?

Selain itu, aliran baru akan mewarisi gaya yang mirip dengan "Flows" .

Dengan permintaan, saya dapat memberikan sampel data.

Aliran tidak akan ada antara fitur "LayerA", tetapi antara fitur "LayerB" . Tujuan utamanya adalah untuk mencapai atribut "FLUX"(yaitu dari / ke) untuk koneksi antara "LayerB"kemungkinan sebagai tabel / Asal-Tujuan Matriks.

Ada beberapa persyaratan / kriteria yang harus dipatuhi:

1. Tidak ada koneksi aliran antara bagian fitur (dipilih dengan warna kuning) di sel yang sama

2. Tidak ada koneksi antara fitur yang sama bahkan bagian-bagiannya ada di sel yang berbeda

3. Koneksi ada di antara bagian-bagian fitur "LayerB"(berdasarkan "Union"output) jika mereka sepenuhnya dalam dua "LayerA"fitur sel yang berbeda

4. Baru"FLUX" yang disampaikan, akan dihitung seperti yang ditunjukkan pada gambar di bawah ini.

Misalnya, ada hubungan antara dua sel Idan II, di mana "FLUX"ada 100. Dengan asumsi nilai-nilai lain, "NEW_FLUX"antara A'dan B''akan ada 1.5625. 100hanya satu contoh.

Referensi:

• Nilai-nilai yang diwariskan antara poligon di QGIS?
• Menyesuaikan gaya peta aliran di QGIS 2
• Hitung luas proporsional poligon dalam poligon lapisan lain
• Sumber GIS Rambling Gratis dan Terbuka | Memindahkan data dalam GIS: masalah & ide
• Mempertahankan geometri dengan Virtual Layer di QGIS 2

Dengan Virtual Layers, secara teoritis, itu mungkin (dengan shapefile, prosesnya akan ekstra panjang, tetapi jika lapisannya ada dalam Database Spasial, saya pikir ini jauh lebih cepat).

Berikut kodenya:

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

Output grafis akan terlihat seperti

Hasilnya diuji secara manual. Perbedaan "FLUX"nilai diabaikan.

Hasil akhir akan mewarisi gaya dari "Flow"dan terlihat seperti

Saya merekomendasikan untuk mengujinya dengan beberapa data, dan jika terlalu lama untuk set data yang besar, jalankan langkah demi langkah kueri ( "inter_ab", "new_flux") dan simpan hasilnya dan jalankan kueri berikutnya.

Anda bisa melakukan gabungan antara tiga lapisan, lalu digabungkan dengan layerB. Lapisan virtual kemungkinan dapat digunakan. Saya tidak yakin apakah data penting ada di layerAatau di dalam flowlapisan .. Berikut ini adalah kemungkinan (yang belum diuji):

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