使用pgrouting进行最短路径搜索

最新推荐文章于 2024-11-18 16:25:08 发布

转载最新推荐文章于 2024-11-18 16:25:08 发布 · 872 阅读

地理几何专栏收录该内容

3 篇文章

订阅专栏

PgRouting是基于PostGIS用于网络分析的扩展模块，依托PostGIS优势保障网络分析效果。本文介绍其使用，包括环境搭建（安装PostgreSQL和PostGIS）、启用功能、路网数据处理、拓扑数据计算、查询最短路径及参数化查询，还提及在geoserver和ol中显示导航路径。

转载 https://www.cnblogs.com/weiweictgu/p/8280285.html

PgRouting是基于开源空间数据库PostGIS用于网络分析的扩展模块，最初它被称作pgDijkstra，因为它只是利用Dijkstra算法实现最短路径搜索，之后慢慢添加了其他的路径分析算法，如A算法，双向A算法，Dijkstra算法，双向Dijkstra算法，tsp货郎担算法等，然后被更名为pgRouting。该扩展库依托PostGIS自身的gist索引，丰富的坐标系与图形类型，强大的几何处理能力，如空间查询，空间处理，线性参考等优势，能保障在较大数据级别下的网络分析效果更快更好。

1、环境搭建

需要先安装PostgreSQL9.6，再安装PostGIS，注意版本号PostGIS2.4后已经集成了pgRouting，无需单独安装

PostGIS 官网 http://www.postgis.org/

pgRouting官网http://pgrouting.org/documentation.html

2、启用pgRouting功能，连接到对应postgres数据库后执行下面的sql语句

CREATE EXTENSION postgis;
CREATE EXTENSION pgrouting;

3、路网数据处理

由于路网分析的特殊性，只支持LineString类型，不支持MultiLineString类型，如果路网数据为MultiLineString类型导入前需加工处理，可使用arcgis在交叉点处使用打断相交线来分割线，参见

http://resources.arcgis.com/zh-cn/help/main/10.1/index.html#//01m8，最终可用的数据必须为路网不自相交，导入时选择简单的线类型，否则网络拓扑计算数据会错误

4、路网拓扑数据计算处理，执行成功后，执行成功后会生产一个**_vertices_pgr的表，里面包含路网相交点的空间数据

复制代码

 alter table road add column source int;

alter table road add column target int;

create index road_source_idx on road("source");
create index road_target_idx on road("target");


ALTER TABLE road  ADD COLUMN length double precision;  


SELECT pgr_createTopology('road',0.00001, 'geom', 'gid');  

update road set length =st_length(geom);

复制代码

5、查询最短路径sql语句

复制代码

SELECT seq, id1 AS node, id2 AS edge, cost,geom  FROM pgr_dijkstra('  
SELECT gid AS id,                      
source::integer,                         
target::integer,                        
length::double precision AS cost  
FROM xmpark_road',  
1, 10, false, false) as di  
join xmpark_road pt  
on di.id2 = pt.gid

复制代码

6、参数化的查询sql语句，可以用下面的语句，在geoserver中使用发布sql视图图层来显示导航搜索图层，见下图

复制代码

SELECT seq, id1 AS node, id2 AS edge, cost,geom  FROM pgr_dijkstra('  
SELECT gid AS id,                      
source::integer,                         
target::integer,                        
length::double precision AS cost  
FROM xmpark_road',  
%a%, %b%, false, false) as di  
join xmpark_road pt  
on di.id2 = pt.gid

复制代码