Transferencia de flujos (conexiones + valores) entre polígonos

En QGIS hay dos archivos de forma que representan los datos en movimiento entre las celdas y una capa adicional, vea la imagen a continuación

Datos móviles definidos por:

• Polígono "LayerA"(cuadrados transparentes con contorno rojo). Además, también se relaciona con círculos que representan los movimientos dentro de las células, visualizados en la posición de los "LayerA"geocentroides.

  

• Capa de polilínea "Flows"(flechas amarillas / grises), transmite valores a través de conexiones entre geocentroides de "LayerA"entidades

  

Capa de destino:

• Polígono "LayerB"(características de color lila claro con contorno gris oscuro).

  

Además, ya he transferido "FLUX"y valores de movimiento dentro de las células de "LayerA"en "LayerB"polígonos, véase mi anterior pregunta: heredada valores entre polígonos en QGIS? . Se realiza mediante el %de $areacálculo.

Puede haber una solución / enfoque significativo para transferir / transmitir / transformar conexiones de flujo representadas por "Flows"y sus valores de relaciones de "LayerA"a relaciones de "LayerB".

¿Cómo puedo lograr esas conexiones como polilíneas?

Además, los nuevos flujos heredarán un estilo similar a "Flows".

Por solicitud, puedo proporcionar una muestra de los datos.

Los flujos no existirán entre las características de "LayerA", sino entre las características de "LayerB" . El objetivo principal es lograr el atributo "FLUX"(es decir, de / a) para conexiones entre "LayerB"posibles como tabla / matriz de origen-destino.

Hay algunos requisitos / criterios que deben cumplirse:

1. No hay conexiones de flujo entre las partes de las características (seleccionadas en amarillo) en la misma celda

2. No hay conexiones entre la misma característica, incluso sus partes están en celdas diferentes

3. Existen conexiones entre partes de características "LayerB"(basadas en la "Union"salida) si están completamente dentro de dos características de "LayerA"celda distintas

4. El nuevo "FLUX"valor que se transmite se calculará como se muestra en la imagen a continuación.

Por ejemplo, hay una conexión entre dos celdas Iy II, donde "FLUX"está 100. Asumiendo otros valores, el "NEW_FLUX"intermedio A'y B''estará alrededor 1.5625. 100Es solo un ejemplo.

Referencias

• ¿Valores heredados entre polígonos en QGIS?
• Ajuste de estilos de mapa de flujo en QGIS 2
• Calcular el área proporcional del polígono dentro del polígono de otra capa
• Free and Open Source GIS Ramblings | Datos de movimiento en SIG: problemas e ideas
• Preservar la geometría con Virtual Layer en QGIS 2

Con las capas virtuales, en teoría, es posible (con archivos de forma, el proceso será más largo, pero si las capas están en una base de datos espacial, creo que es mucho más rápido).

Aquí el código:

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

La salida gráfica se verá así

El resultado fue probado manualmente. La diferencia de "FLUX"valores es despreciable.

El resultado final heredará estilos de "Flow"y se verá como

Recomiendo probarlo con unos pocos datos, y si toma demasiado tiempo para grandes conjuntos de datos, ejecute las consultas paso a paso ( "inter_ab", "new_flux") y guarde el resultado y ejecute la siguiente consulta.

Puede hacer una unión entre las tres capas, luego agregar por layerB. Es probable que se puedan usar capas virtuales. No estoy seguro de si los datos importantes están en layerAo en la flowcapa. Aquí hay una posibilidad (no probada):

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