如何使用PostGIS计算连接所有点的最小网络？

我有一组postgis脚本，它们生成两个表-一个是一组点，第二个是围绕它们的道路。所有数据都在同一投影中，并且两个输出都存储在带有postgis 2.1的postgres 9.2表中

道路网的注浆拓扑已创建，点表中的列包含最近的路段。

然后，我想生成道路网的一个子集，该子集代表使用最小生成树之类的东西连接所有点的最小网络。道路网络是无方向性的，成本只是路线长度。

我可以使用v.net系列模块在QGIS / Grass中执行此操作，但理想情况下，我也希望在SQL中保留最后一步。

我已经看过新的apspWarshall postgis函数，但是对于如何鼓励它将精力集中在连接点而不是整个网络上，我感到困惑。

这是我编写的简短脚本，试图创建一个框架来解决此问题，但我看不到在哪里可能集中功能从边缘的子集开始。

SELECT seq, id1 AS node, id2 AS edge, cost, the_geom
FROM   pgr_apspWarshall('SELECT gid AS id, 
                                source, 
                                target, 
                                st_length(the_geom) AS cost 
                         FROM   road_network
                        ',
                        false, false
                       ) AS tree
JOIN   road_network As roads
ON     tree.id2 = roads.gid

在单路径最短路径问题中，该功能要求起点和终点，但显然不是在所有问题上都要求。同样在草丛中，v.net.spanningtree和v.net.steiner期望将一组点和线作为组合网络使用。

有没有人对如何在PostGIS中做到这一点提出建议？

该答案尚未完成或未通过测试，但请尝试以下操作：

根据问题/ 39210：

with index_query as (
SELECT
        ,ST_Distance(i.geom, i.b_geom) AS dist
        ,ST_MakeLine(i.geom, i.b_geom) as geom
FROM(
SELECT
        ,a.geom
        ,b.geom AS b_geom
        ,rank() OVER (PARTITION BY a.id ORDER BY ST_Distance(a.centroid_geom, b.the_geom)) AS pos
FROM points a, points b 
WHERE a.id <> b.id) i
WHERE pos = 1
) select ST_Union(geom) from index_query;

我认为这不是很有效。

@Adrian，我真的不了解灌浆的结果，但是文档非常详细。我的答案基于一个两步函数，该函数在SQL中将是非常无用的，但是[可能]会产生结果。此[未测试]解决方案不会优化哪个是最佳起点，但是会将整个路线网络缩小到仅连接所有站点的边缘，然后有效地路由到所有站点。

第1步（连接所有站点的道路网络子集的子选择）这使用多目标（K Dijkstr路径）路由功能返回一组路径（当成本<> -1时）确实连接了所有停止。

SELECT id1 as path, st_astext(st_linemerge(st_union(b.the_geom))) as the_geom
FROM pgr_kdijkstraPath(
’SELECT id, source, target, cost FROM edge_table’,
min(all_your_stop_ids), [array_of_all_your_stop_ids], false, false
) a,
edge_table b
WHERE a.id3=b.id
GROUP by id1
ORDER by id1

我这里的问题是从停止表组装数组的语法，因为在问题中并未真正描述。但是，我们假设SQL语法可以汇编该数组，并且最小id停止应该是到其余目标停止的所有K条路径的起点。

第2步（根据连接所有站点的上述子集道路网络路径，最终选择最小路径）这实际上是您开始的工作，但是我建议您将道路网络与id1（路径）的初始结果等价加入这样在最终的“实战”路线中仅使用道路的一部分：

SELECT seq, id1 AS node, id2 AS edge, cost, the_geom
FROM   pgr_apspWarshall('SELECT R.gid AS id, 
                                R.source, 
                                R.target, 
                                st_length(R.the_geom) AS cost 
             FROM   road_network AS R JOIN
                   (SELECT id1 as path
                     FROM pgr_kdijkstraPath(
                            ’SELECT id, source, target, cost FROM edge_table’,
                            min(all_your_stop_ids), 
                            [array_of_all_your_stop_ids], false, false
                           ) a,
                     edge_table b
                    WHERE a.id3=b.id
                    GROUP by id1
                    ORDER by id1
                        ',
                        false, false
                  ) AS  Just_K_Paths
         on R.id1 = just_K_paths.id1',       /* this join reduces R down to K paths */
         false, false
        ) AS tree
  JOIN   road_network As roads
  ON     tree.id2 = roads.gid

因此，总而言之...内部的k_dijkstra_path路由查询将整个道路网络减少到仅连接所有站点的路径，然后外部的fField_Warshal路由仅使用那些边缘ID来解决路径优化查询....也许。