Trasferimento di flussi (connessioni + valori) tra poligoni

In QGIS ci sono due shapefile che rappresentano i dati in movimento tra le celle e un livello aggiuntivo, vedi immagine sotto

Spostamento dei dati definito da:

• Poligono "LayerA"(quadrati trasparenti con contorno rosso). Inoltre si riferisce anche ai cerchi che rappresentano i movimenti all'interno delle cellule, visualizzati sulla posizione dei "LayerA"geocentroidi.

  

• Lo strato di polilinea "Flows"(frecce gialle / grigie) trasmette valori tramite connessioni tra geocentroidi di "LayerA"caratteristiche

  

Livello target:

• Poligono "LayerB"(lineamenti lilla chiaro con contorno grigio scuro).

  

Inoltre, ho già trasferito "FLUX"e valori di movimento all'interno delle cellule da "LayerA"in "LayerB"poligoni, vedere la mia domanda precedente: Ereditato valori compresi tra i poligoni in QGIS? . È stato fatto usando il %di$area calcolo.

Potrebbe esserci una soluzione / approccio significativo per trasferire / trasmettere / trasformare connessioni di flusso rappresentate da "Flows"e i suoi valori da relazioni di "LayerA"a relazioni di"LayerB" .

Come posso ottenere tali connessioni come polilinee?

Inoltre, i nuovi flussi erediteranno uno stile simile a "Flows" .

Su richiesta, posso fornire un campione di dati.

I flussi non esistono tra le funzionalità di "LayerA", ma tra le funzionalità di "LayerB" . L'obiettivo principale è quello di ottenere l'attributo "FLUX"(cioè da / a) per le connessioni tra "LayerB"possibili come tabella / matrice origine-destinazione.

Ci sono alcuni requisiti / criteri che dovrebbero essere rispettati:

1. Non ci sono connessioni di flusso tra le parti delle caratteristiche (selezionate in giallo) nella stessa cella

2. Non ci sono connessioni tra la stessa funzione anche se le sue parti sono in celle diverse

3. Esistono connessioni tra parti di funzioni "LayerB"(basate "Union"sull'output) se si trovano interamente in due "LayerA"funzioni di celle distinte

4. Nuovo"FLUX" valore trasmesso verrà calcolato come mostrato nell'immagine seguente.

Ad esempio, v'è una connessione tra due celle Ie II, dove "FLUX"è 100. Supponendo altri valori, "NEW_FLUX"tra A'e B''sarà in giro 1.5625. 100è solo un singolo esempio.

Riferimenti:

• Valori ereditati tra poligoni in QGIS?
• Regolazione degli stili delle mappe di flusso in QGIS 2
• Calcola l'area proporzionale del poligono all'interno del poligono di un altro livello
• Rambling GIS gratuiti e open source | Dati di movimento in GIS: problemi e idee
• Preservare la geometria con Virtual Layer in QGIS 2

Con i livelli virtuali, in teoria, è possibile (con gli shapefile, il processo sarà molto lungo, ma se i livelli si trovano in un database spaziale, penso che sia molto più veloce).

Ecco il codice:

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

L'output grafico sarà simile

Il risultato è stato testato manualmente. La differenza di "FLUX"valori è trascurabile.

L'output finale erediterà "Flow"e assomiglierà agli stili

Consiglio di testarlo con pochi dati e, se impiega troppo tempo per set di dati di grandi dimensioni, esegui le query passo dopo passo ("inter_ab" , "new_flux") e salvi il risultato ed esegui la query successiva.

È possibile eseguire un join tra i tre livelli, quindi aggregare per layerB. È possibile utilizzare probabilmente i livelli virtuali. Non sono sicuro che i dati importanti si trovino nel layerAo nel flowlayer. Ecco una possibilità (non testata):

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