在多边形之间传输流（连接+值）

在QGIS中，有两个shapefile表示单元之间和另一层之间的移动数据，请参见下图。

移动数据的定义：

• 多边形"LayerA"（带有红色轮廓的透明正方形）。此外，它还涉及代表细胞内部运动的圆圈，这些圆圈在质心的位置上可视化"LayerA"。

  

• 折线层"Flows"（黄色/灰色箭头），经由连接传达的geocentroids之间的值"LayerA"的特征

  

目标层：

• 多边形"LayerB"（浅紫色特征，深灰色轮廓）。

  

另外，我已经将"FLUX"单元内的值转移"LayerA"到了"LayerB"多边形中，并将其移动到多边形中，请参阅我的上一个问题：QGIS中多边形之间的继承值？。它用的是完成%的$area计算。

可能存在一种有意义的解决方案/方法，将/ 表示的流 连接"Flows"及其值从/的关系"LayerA"转换为/的关系"LayerB"。

如何将这些连接实现为折线？

此外，新流程将继承与相似的样式"Flows"。

根据要求，我可以提供数据样本。

流将不存在于特征之间"LayerA"，而存在于特征之间 "LayerB"。主要目的是实现可能的表/来源-目的地矩阵"FLUX"之间的连接属性（即，从/到）"LayerB"。

需要遵守一些要求/标准：

1.在同一单元格中，要素的各个部分（以黄色选中）之间没有流路连接

2.相同要素之间没有连接，即使其零件位于不同的单元中

3.如果要素的各个部分完全位于两个不同的像元要素中，则它们之间存在连接"LayerB"（基于"Union"输出）"LayerA"

4."FLUX"如下图所示，将计算正在传送的新值。

举例来说，有两个单元之间的连接I和II，其中"FLUX"是100。假设其他值，"NEW_FLUX"介于A'和之间B''将是1.5625。100仅是一个例子。

参考文献：

• QGIS中多边形之间的继承值？
• 在QGIS 2中调整流程图样式
• 计算另一层多边形内多边形的比例面积
• 免费和开源GIS乱流| GIS中的运动数据：问题和想法
• 在QGIS 2中使用虚拟层保留几何

从理论上讲，使用虚拟层是可行的（使用shapefile，过程会很长，但是如果这些层位于空间数据库中，我认为它会快得多）。

这里的代码：

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

图形输出看起来像

手动测试结果。"FLUX"值的差异可以忽略不计。

最终输出将继承样式"Flow"，看起来像

我建议用几个数据来测试它，如果它的时间太长，对于大数据集，执行一步一步的查询（"inter_ab"，"new_flux"），并保存结果并执行一个查询。

您可以在三层之间进行联接，然后通过进行聚合layerB。虚拟层可能会被使用。我不确定重要数据是在图层中layerA还是flow图层中。这是（未试用的）可能性：

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