Postgresql+Geoserver+Openlayers+pgRouting学习1
pgRouting 学习 最短路径算法
1、数据预处理
数据的投影需要定义，空间数据库中矢量图层的数据投影需要定义，而且需要同geoserver以及openlayer中投影定义一致，我使用的是EPSG:4326
编辑矢量图层数据
最短路径使用的 Dijkstra算法，了解这个算法对于数据的要求是矢量数据要是一张图，意味着如果相邻的两个节点之间连接的线段是独立的，就是说如果有一条街道和别的街道相交就需要把它拆分。
将矢量数据导入postgresql
postgresql处理线表
--ShortestPath为表名
--source为线表起点字段名称
--target为线表终点字段名称
--source和target可以在导入数据库前直接先自己定义
--如果不自己定义起点终点，则自动生成起点和终点点号
ALTER TABLE shortestpath ADD COLUMN "source" integer;
ALTER TABLE shortestpath ADD COLUMN "target" integer;
--创建拓扑
SELECT pgr_createTopology('shortestpath', 0.00001, 'geom', 'id');
--为起点号终点号加空间索引
CREATE INDEX source_idx ON ShortestPath("source");
CREATE INDEX target_idx ON ShortestPath("target");
--添加长度字段、并计算赋值
ALTER TABLE shortestpath ADD COLUMN length double precis;
update shortestpath set length =st_length(geom);
--将长度值赋给reverse_cost，作为路线选择标准
ALTER TABLE shortestpath ADD COLUMN reverse_cost double precision;
UPDATE shortestpath SET reverse_cost = length;
构建路径函数
--删除已存在的函数
DROP FUNCTION pgr_fromAtoB(tbl varchar,startx float, starty float,endx float,endy float);
--tbl路网表名
--startx起点经度
--starty起点纬度
--endx终点经度
--endy终点纬度
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_startSource integer;
v_endSource integer;--距离终点最近线的起点
v_endTarget integer;
v_statpoint geometry;--在v_startLine上距离起点最近的点
v_endpoint geometry;--在v_endLine上距离终点最近的点
v_res geometry;--最短路径分析结果
v_res_a geometry;
v_res_b geometry;
v_res_c geometry;
v_res_d geometry;
v_perStart float;--v_statpoint在v_res上的百分比
v_perEnd float;--v_endpoint在v_res上的百分比
v_shPath_se geometry;--开始到结束
v_shPath_es geometry;--结束到开始
v_shPath geometry;--最终结果
tempnode float;
begin
--查询离起点最近的线
--4326坐标系
--找起点15米范围内的最近线
execute 'select geom, source, target from ' ||tbl||
' where ST_DWithin(geom,ST_Geometryfromtext(''point('|| startx ||' ' || starty||')'',4326),15)
order by ST_Distance(geom,ST_GeometryFromText(''point('|| startx ||' '|| starty ||')'',4326)) limit 1'
into v_startLine, v_startSource ,v_startTarget;
--查询离终点最近的线
--找终点15米范围内的最近线
execute 'select geom, source, target from ' ||tbl||
' where ST_DWithin(geom,ST_Geometryfromtext(''point('|| endx || ' ' || endy ||')'',4326),15)
order by ST_Distance(geom,ST_GeometryFromText(''point('|| endx ||' ' || endy ||')'',4326)) limit 1'
into v_endLine, v_endSource,v_endTarget;
--如果没找到最近的线，就返回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 ||')',4326)) into v_statpoint;
select ST_ClosestPoint(v_endLine, ST_GeometryFromText('point('|| endx ||' ' || endy ||')',4326)) into v_endpoint;
-- ST_Distance
--从开始的起点到结束的起点最短路径
execute 'SELECT st_linemerge(st_union(b.geom)) ' ||
'FROM pgr_kdijkstraPath(
''SELECT gid as id, source, target, length as cost FROM ' || tbl ||''','
||v_startSource || ', ' ||'array['||v_endSource||'] , false, false
) a, '
|| tbl || ' b
WHERE a.id3=b.gid
GROUP by id1
ORDER by id1' into v_res ;
--从开始的终点到结束的起点最短路径
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_b ;
--从开始的起点到结束的终点最短路径
execute 'SELECT st_linemerge(st_union(b.geom)) ' ||
'FROM pgr_kdijkstraPath(
''SELECT gid as id, source, target, length as cost FROM ' || tbl ||''','
||v_startSource || ', ' ||'array['||v_endTarget||'] , false, false
) a, '
|| tbl || ' b
WHERE a.id3=b.gid
GROUP by id1
ORDER by id1' into v_res_c ;
--从开始的终点到结束的终点最短路径
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_endTarget||'] , false, false
) a, '
|| tbl || ' b
WHERE a.id3=b.gid
GROUP by id1
ORDER by id1' into v_res_d ;
if(ST_Length(v_res) > ST_Length(v_res_b)) then
v_res = v_res_b;
end if;
if(ST_Length(v_res) > ST_Length(v_res_c)) then
v_res = v_res_c;
end if;
if(ST_Length(v_res) > ST_Length(v_res_d)) then
v_res = v_res_d;
end if;
--如果找不到最短路径，就返回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;
2、geoserver 发布图层
geoserver可以发布矢量图层作为地图，sql视图图层作为分析图层
SELECT * FROM pgr_fromAtoB('shortestpath', %x1%, %y1%, %x2%, %y2%)
从sql语句来看，实则调用空间数据库pgr_fromAtoB函数，将起点和终点坐标作为参数。
3、前端调用
测试可使用url测试，如返回图片则成功；如果出现错误，则在geoserver控制台查看错误信息，我遇到的错误有：
-
混合投影
解决办法：定义空间数据库表格投影
-
数据格式不一致
解决办法：修改空间数据库表格类型
http://localhost:8082/geoserver/jinmen/wms?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&FORMAT=image%2Fpng&TRANSPARENT=true&LAYERS=jinmen%3AShortest%20Path&viewparams=x1%3A120.53052%3By1%3A30.637777%3Bx2%3A120.55427%3By2%3A30.62243&CRS=EPSG%3A4326&STYLES=&FORMAT_OPTIONS=dpi%3A124&WIDTH=789&HEIGHT=616&BBOX=30.61324977874755%2C120.52187555486509%2C30.65170192718505%2C120.5711267644709
