基于pgrouting的任意两点间的最短路径查询函数二
来源:互联网 发布:js 数组 元素个数 编辑:程序博客网 时间:2024/05/19 20:43
在前面的博文中写过一篇查询任意两点间最短路径的函数,当时对pgrouting不熟悉,功能很low。现在对该函数进行扩展,支持用户自己输入查询的数据库表,这一点看似简单,其实意义很大,在做室内导航的时候当用户所在的楼层变化的时候最短路径函数查询的数据表名称也会发生变化,不可能一栋大楼里的道路都是一样的吧,另外进行跨楼层的最短路径规划时,需要查询从A到楼梯口的最短路径和楼梯口到B的最短路径,这些都需要进行最短路径规划的时候能够自己选择数据表。
先解释一下最短路径规划的处理步骤,首先要确定用户的出发点和目的地所在的道路,再在相应的道路上确定道路的节点,查找这两个节点之间的最短路径,最后再处理出发点和目的地到道路节点之间的路段。具体过程为:
(1)查找距离用户出发点最近的道路和该道路的终点T。
(2)查找距离用户目的地最近的道路和该道路的起点S。
(3)计算前两步找出的两点之间的最短路径。
(4)处理出发点和道路终点T以及目的去和道路起点S之间的路段。
DROP FUNCTION pgr_fromAtoB(tbl varchar,startx float, starty float,endx float,endy float);CREATE OR REPLACE function pgr_fromAtoB(tbl varchar,startx float, starty float,endx float,endy float) returns geometry as $body$ declare v_startLine geometry;--离起点最近的线 v_endLine geometry;--离终点最近的线 v_startTarget integer;--距离起点最近线的终点 v_endSource integer;--距离终点最近线的起点 v_statpoint geometry;--在v_startLine上距离起点最近的点 v_endpoint geometry;--在v_endLine上距离终点最近的点 v_res geometry;--最短路径分析结果 v_perStart float;--v_statpoint在v_res上的百分比 v_perEnd float;--v_endpoint在v_res上的百分比 v_shPath geometry;--最终结果 tempnode float;begin --查询离起点最近的线 execute 'select geom ,target from ' ||tbl||' where ST_DWithin(geom,ST_Geometryfromtext(''point('||startx ||' ' || starty||')''),15) order by ST_Distance(geom,ST_GeometryFromText(''point('|| startx ||' '|| starty ||')'')) limit 1' into v_startLine ,v_startTarget; --查询离终点最近的线 execute 'select geom,source from ' ||tbl||' where ST_DWithin(geom,ST_Geometryfromtext(''point('|| endx || ' ' || endy ||')''),15) order by ST_Distance(geom,ST_GeometryFromText(''point('|| endx ||' ' || endy ||')'')) limit 1' into v_endLine,v_endSource; --如果没找到最近的线,就返回null if (v_startLine is null) or (v_endLine is null) then return null; end if ; select ST_ClosestPoint(v_startLine, ST_Geometryfromtext('point('|| startx ||' ' || starty ||')')) into v_statpoint; select ST_ClosestPoint(v_endLine, ST_GeometryFromText('point('|| endx ||' ' || endy ||')')) into v_endpoint; --最短路径 execute 'SELECT st_linemerge(st_union(b.geom)) ' || 'FROM pgr_kdijkstraPath( ''SELECT gid as id, source, target, length as cost FROM ' || tbl ||''',' ||v_startTarget || ', ' ||'array['||v_endSource||'] , false, false ) a, ' || tbl || ' b WHERE a.id3=b.gid GROUP by id1 ORDER by id1' into v_res; --如果找不到最短路径,就返回null if(v_res is null) then return null; end if; --将v_res,v_startLine,v_endLine进行拼接 select st_linemerge(ST_Union(array[v_res,v_startLine,v_endLine])) into v_res; select ST_Line_Locate_Point(v_res, v_statpoint) into v_perStart; select ST_Line_Locate_Point(v_res, v_endpoint) into v_perEnd; if(v_perStart > v_perEnd) then tempnode = v_perStart;v_perStart = v_perEnd;v_perEnd = tempnode; end if; --截取v_res SELECT ST_Line_SubString(v_res,v_perStart, v_perEnd) into v_shPath; return v_shPath; end; $body$ LANGUAGE plpgsql VOLATILE STRICT
使用这个最短路径规划函数进行路径规划一般不会出现问题,但是有的时候会出现这样的问题,因为函数的查找思路是查找距离起点最近的线和改线的终点,查找离目的地最近的线和该线的起点,这样有可能会出现这样的情况,就是起点和终点对应的线段是相连接的,这时查找的距离起点最近线的终点和距离目的地最近线上的起点就是同一个点,这样的话这两点直接就没有最短距离。
所以需要对上述代码第59到61行进行修改,修改后的代码如下:
DROP FUNCTION pgr_fromAtoB(tbl varchar,startx float, starty float,endx float,endy float);CREATE OR REPLACE function pgr_fromAtoB(tbl varchar,startx float, starty float,endx float,endy float) returns geometry as $body$ declare v_startLine geometry;--离起点最近的线 v_endLine geometry;--离终点最近的线 v_startTarget integer;--距离起点最近线的终点 v_endSource integer;--距离终点最近线的起点 v_statpoint geometry;--在v_startLine上距离起点最近的点 v_endpoint geometry;--在v_endLine上距离终点最近的点 v_res geometry;--最短路径分析结果 v_perStart float;--v_statpoint在v_res上的百分比 v_perEnd float;--v_endpoint在v_res上的百分比 v_shPath geometry;--最终结果 tempnode float;begin --查询离起点最近的线 execute 'select geom ,target from ' ||tbl||' where ST_DWithin(geom,ST_Geometryfromtext(''point('||startx ||' ' || starty||')''),15) order by ST_Distance(geom,ST_GeometryFromText(''point('|| startx ||' '|| starty ||')'')) limit 1' into v_startLine ,v_startTarget; --查询离终点最近的线 execute 'select geom,source from ' ||tbl||' where ST_DWithin(geom,ST_Geometryfromtext(''point('|| endx || ' ' || endy ||')''),15) order by ST_Distance(geom,ST_GeometryFromText(''point('|| endx ||' ' || endy ||')'')) limit 1' into v_endLine,v_endSource; --如果没找到最近的线,就返回null if (v_startLine is null) or (v_endLine is null) then return null; end if ; select ST_ClosestPoint(v_startLine, ST_Geometryfromtext('point('|| startx ||' ' || starty ||')')) into v_statpoint; select ST_ClosestPoint(v_endLine, ST_GeometryFromText('point('|| endx ||' ' || endy ||')')) into v_endpoint; --最短路径 execute 'SELECT st_linemerge(st_union(b.geom)) ' || 'FROM pgr_kdijkstraPath( ''SELECT gid as id, source, target, length as cost FROM ' || tbl ||''',' ||v_startTarget || ', ' ||'array['||v_endSource||'] , false, false ) a, ' || tbl || ' b WHERE a.id3=b.gid GROUP by id1 ORDER by id1' into v_res ; --如果找不到最短路径,就返回null --if(v_res is null) then -- return null; --end if; --将v_res,v_startLine,v_endLine进行拼接 select st_linemerge(ST_Union(array[v_res,v_startLine,v_endLine])) into v_res; select ST_Line_Locate_Point(v_res, v_statpoint) into v_perStart; select ST_Line_Locate_Point(v_res, v_endpoint) into v_perEnd; if(v_perStart > v_perEnd) then tempnode = v_perStart;v_perStart = v_perEnd;v_perEnd = tempnode; end if; --截取v_res SELECT ST_Line_SubString(v_res,v_perStart, v_perEnd) into v_shPath; return v_shPath; end; $body$ LANGUAGE plpgsql VOLATILE STRICT
3 0
- 基于pgrouting的任意两点间的最短路径查询函数二
- 基于pgrouting的任意两点间的最短路径查询函数
- 基于pgrouting求任意两点的最短路径的函数pgr_fromAtoB
- 使用pgrouting求任意两点的最短路径
- pgrouting路径分析(任意两点之间的最短路径)
- 任意两点间的最短路径
- 求任意两点间的最短路径
- 任意两点间的最短路径---floyd_warshall算法
- floyd算法(求任意两点间的最短路径)
- Floyd算法(任意两点间的最短路径)
- POJ 2139-- 最短路径 (floyd算法,任意两点间的最短路径)
- pgrouting最短路径
- 运用Floyd算法求得带权有向图任意两点间的最短路径C/C++
- Dijkstra单源最短路径实现 及 Floyd任意两点之间的最短路径
- 任意两点之间最短路径
- 【OI杂记】求二叉树上任意两点的最短路径上的边权最大值
- 两点之间的所有最短路径
- 软件测试作业 图中任意两点的最短路径
- 自我介绍
- 斯坦福《机器学习》Lesson4感想-------2、广义线性模型
- CListCtrl 虚拟列表显示复选框
- POJ 3311 Hie with the Pie (状压DP)
- PCB电路板元器件布局的一般原则
- 基于pgrouting的任意两点间的最短路径查询函数二
- msp430学习笔记之定时器A
- grep与正则表达式
- ActiveX 属性持久性 是否是测试环境判断
- iOS-compile_debug-crash_record
- 主流图片加载框架对比,及遇到的问题
- AVL Tree(平衡二叉树)
- ANE原生代码的调试(安卓)++flex通过ANE调试原生安卓代码
- 1153 选择子序列