使用pgrouting求任意两点的最短路径

要利用rgrouting实现像QGIS那样求任意两点间的最短路径，可以按照以下步骤使用pl/pgsql进行自定义函数：

1 函数的参数为：_myShortPath(startxfloat,starty float,endx float,endy float,costfile varchar)，前四个参数为起点和终点的坐标，最后一个参数为最短路径分析时的cost字段；

2 分别找出距离起点和终点最近的两条线，可以使用：

ST_DWithin(geography gg1,geography gg2,double precisiondistance_meters);

在距离gg2在distance_meters范围内的gg1的线，之后使用ST_distance(geometryg1,geometry g2);进行排序，并取第一个。

3 找出在上一步求出最近线上距离起点（终点）最近的点。使用：

ST_ClosestPoint(geometryg1,geometry g2)

4 根据第2步求出的离起点最近先的target和离终点最近线的source，求最短路径。使用

SELECTst_linemerge(st_union(b.geom)) into v_res

FROMpgr_kdijkstraPath(

'SELECTid,source,target,' || costfile ||' as cost FROM route',

v_startTarget,array[v_endSource],true,false

)a,

routeb

WHEREa.id3=b.id

GROUPby id1

ORDERby id1;

5 将求出的路径与第2步找出的两条线进行合并，可以使用：

st_linemerge(ST_Union(array[v_res,v_startLine,v_endLine]))

6 求出第3步找出的点在上一步中线的上的百分比，可以使用：

ST_Line_Locate_Point(v_res,v_statpoint)

7 使用上一步求出的百分比来截取第5步的结果：

ST_Line_SubString(v_res,v_perStart,v_perEnd)

8 最后将截取的结果返回；

具体代码：

create function _myShortPath(startx float,costfile varchar) 
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;--最终结果
begin	
        
	--查询离起点最近的线
	select geom,target into v_startLine,v_startTarget from  route where
	ST_DWithin(geom,ST_Geometryfromtext('point('|| startx ||' ' || starty ||')'),0.015)
	order by ST_distance(geom,ST_GeometryFromText('point('|| startx ||' '|| starty ||')'))  limit 1;
	
	--查询离终点最近的线
	select geom,source into v_endLine,v_endSource from  route where
	ST_DWithin(geom,ST_Geometryfromtext('point('|| endx || ' ' || endy ||')'),ST_GeometryFromText('point('|| endx ||' ' || endy ||')'))  limit 1;

	--如果没找到最近的线，就返回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;

	
        --最短路径
	SELECT st_linemerge(st_union(b.geom)) into v_res
	FROM pgr_kdijkstraPath(
	'SELECT id,v_startTarget,false
	) a,route b
	WHERE a.id3=b.id
	GROUP by id1
	ORDER by id1;

	--如果找不到最短路径，就返回null
	if(v_res is null) then
		return null;
	end if;
	
	--将v_res,v_endLine进行拼接
	select  st_linemerge(ST_Union(array[v_res,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;

	--截取v_res
	SELECT ST_Line_SubString(v_res,v_perEnd) into v_shPath;
	 
	return v_shPath;

	
	
end;
$body$
LANGUAGE plpgsql VOLATILE STRICT

select_myShortPath(119.274,26.0649,119.258,26.019,'length') intoshortPath；

在QGIS中打开shortPath，可以看出结果与在QGIS执行最短路径结果一致：

使用pgrouting求任意两点的最短路径

相关推荐