java gis 最短路径_用Postgis算最短路径(在任意位置选择起点终点)(下)
今天将上一次路径绘制的坑(上次绘制的是子路径端点的连线)给填了,本文中我将实现一个输入起点和终点经纬度,输出完整路径(包含子路径端点间的真实弯曲情况,完美贴合真实路径)的函数,能够计算准确的路径长度。
经过一番sql的洗礼我终于能够获取完整的规划路径了:
完整路径,完美贴合
路径端点连线,不完美贴合
并且我们将路径规划做成一个Sql函数,要规划路线的话只需要调用函数即可。
select whz_getPathBeta(106.7267,26.5690,106.8046,26.56904,1000);
废话不多说,上sql。
首先是生成路径的sql,参数填写起点终点经纬度和poi搜索范围。函数中首先将用户的起点终点更新入表中,然后获取路径段的数目,最后获取完整的路径图形。
-- 路网表daolu_simp
-- POI表pointsOfInterest2
CREATE OR REPLACE FUNCTION whz_getPathBeta(x1 float8, y1 float8, x2 float8, y2 float8, roadWithin numeric) RETURNS setof geometry AS
$$
declare poiId1 integer;
declare poiId2 integer;
declare segCount integer;
begin
poiId1=whz_updatePoi(x1,y1,roadWithin);
poiId2=whz_updatePoi(x2,y2,roadWithin);
select count(*) into segCount from pgr_withPoints(
'SELECT gid as id, source, target, cost_len as cost, rcost_len as reverse_cost FROM daolu_simp ORDER BY gid',
'SELECT pid as pid, edge_id, fraction, side from pointsOfInterest2',
-poiId1, -poiId2);
return query SELECT whz_getSegGeom(wp.seq,wp.node::integer,wp.edge::integer, segCount,-poiId1,-poiId2) as the_geom FROM pgr_withPoints(
'SELECT gid as id, source, target, cost_len as cost, rcost_len as reverse_cost FROM daolu_simp ORDER BY gid',
'SELECT pid as pid, edge_id, fraction, side from pointsOfInterest2',
-poiId1, -poiId2) as wp;
end;
$$
LANGUAGE plpgsql;
;
然后是往里面插入POI(起点、终点、途径点)的sql:
-- POI表pointsOfInterest2
-- 指定SRID=4326
-- 指定改点在路段的两侧
CREATE OR REPLACE FUNCTION whz_updatePoi(xo float8, yo float8, roadWithin numeric) RETURNS integer AS
$$
declare idRet numeric;
begin
INSERT INTO pointsOfInterest2(x, y, edge_id, side, fraction, the_geom, newPoint)
(SELECT xo, yo,gid,'b',ST_LineLocatePoint(geom, concat('SRID=4326;POINT(',xo,' ',yo,')')::geometry) as fraction,
st_makePoint(xo,yo),
ST_LineInterpolatePoint(geom, ST_LineLocatePoint(geom, concat('SRID=4326;POINT(',xo,' ',yo,')')::geometry))
FROM daolu_simp
WHERE ST_DWithin(geom, concat('SRID=4326;POINT(',xo,' ',yo,')')::geometry, roadWithin)
order by geom concat('SRID=4326;POINT(',xo,' ',yo,')')::geometry
limit 1)
RETURNING pid into idRet;
return idRet;
end;
$$
LANGUAGE plpgsql;
;
最后是获取每一段子路径的完整路径的sql:
-- POI表pointsofinterest2
-- 路网节点表daolu_simp_vertices_pgr
-- 路网表daolu_simp
-- 本方法试用于起点终点都在poi,且一起点一终点
CREATE OR REPLACE FUNCTION whz_getSegGeom(myseq integer, mynode integer, myedge integer, node_num integer,poiId1 integer,poiId2 integer) RETURNS geometry AS
$$
declare geomRet geometry;
declare startFraction float8;
declare endFraction float8;
declare edgeidForPoi integer;
declare pointTmp geometry;
declare nodeTmp integer;
declare edgeTmp geometry;
begin
IF(mynode<0 AND myseq=1)THEN
SELECT node into nodeTmp FROM pgr_withPoints(
'SELECT gid as id, source, target, cost_len as cost, rcost_len as reverse_cost FROM daolu_simp ORDER BY gid',
'SELECT pid as pid, edge_id, fraction, side from pointsOfInterest2',
poiId1, poiId2) as wp1
where wp1.seq= 2;
select pointsofinterest2.fraction into endFraction
from pointsofinterest2
where pointsofinterest2.pid=-mynode
limit 1;
select pointsofinterest2.edge_id into edgeidForPoi
from pointsofinterest2
where pointsofinterest2.pid=-mynode
limit 1;
select the_geom into pointTmp
from daolu_simp_vertices_pgr
where daolu_simp_vertices_pgr.id= nodeTmp
limit 1;
select geom into edgeTmp
from daolu_simp
where daolu_simp.gid=edgeidForPoi
limit 1;
startFraction=st_linelocatepoint(edgeTmp, pointTmp);
if(startFraction>=endFraction) THEN
geomRet=ST_LineSubstring(edgeTmp, endFraction,startFraction);
ELSE geomRet=ST_LineSubstring(edgeTmp, startFraction,endFraction);
end if;
ELSIF(mynode<0 AND myseq=node_num)THEN
SELECT node into nodeTmp FROM pgr_withPoints(
'SELECT gid as id, source, target, cost_len as cost, rcost_len as reverse_cost FROM daolu_simp ORDER BY gid',
'SELECT pid as pid, edge_id, fraction, side from pointsOfInterest2',
poiId1, poiId2) as wp2
where wp2.seq= myseq-1;
select edge_id into edgeidForPoi
from pointsofinterest2
where pointsofinterest2.pid=-mynode
limit 1;
select fraction into endFraction
from pointsofinterest2
where pointsofinterest2.pid=-mynode
limit 1;
select the_geom into pointTmp
from daolu_simp_vertices_pgr
where daolu_simp_vertices_pgr.id= nodeTmp
limit 1;
select geom into edgeTmp
from daolu_simp
where daolu_simp.gid=edgeidForPoi
limit 1;
startFraction=st_linelocatepoint(edgeTmp, pointTmp);
if(startFraction>=endFraction) THEN
geomRet=ST_LineSubstring(edgeTmp, endFraction,startFraction);
ELSE geomRet=ST_LineSubstring(edgeTmp, startFraction,endFraction);
end if;
ELSIF(mynode>0 AND myseq
select geom into geomRet
from daolu_simp
where daolu_simp.gid=myedge;
END IF;
return geomRet;
end;
$$
LANGUAGE plpgsql;
注意,这个sql里面的表名和相关参数是写死的,实际使用的时候需要根据情况调整。
实际应用
例如可以借此实现Cesium的路径规划。在Web端的业务中,使用AJAX请求后端(例如GeoServer),后端再使用Mybatis或者Sql视图(Sql视图在GeoServer中配置),调用数据库中写好的该路径规划函数,函数将返回路径的几何图形,后端再将这些几何图形返回给前端进行渲染。
在倾斜摄影上渲染路径规划结果
java gis 最短路径_用Postgis算最短路径(在任意位置选择起点终点)(下)相关推荐
- 用Postgis算最短路径(在任意位置选择起点终点)
1.前言 阅读本文需要知道什么是shapfile,什么是路径分析,什么是GIS.相比Arcgis的路径分析功能,本文介绍的方法稍微复杂,需要注意的细节更多,但却是完全免费的.PostGis+QGIS+ ...
- java gis项目_记录:GIS常用开发工具
一.GIS软件 ArcGIS Erdas Envi QGIS:开源 国产:超图.中地MapGIS 二.二次开发编程语言 页面前端:html.css.JavaScript .JQuery.Bootstr ...
- python贪心算法最短路径_贪心算法之最短路径
1.Dijkstra算法 Dijkstra算法(单源最短路径) 单源最短路径问题,即在图中求出给定顶点到其它任一顶点的最短路径.在弄清楚如何求算单源最短路径问题之前,必须弄清楚最短路径的最优子结构性质 ...
- Java用蚁群算法求最短路径_蚁群算法最短路径java
改进蚁群算法求解最短路径问题 袁亚博,刘羿,吴斌 [摘要]摘要:针对蚁群算法在求解最短路径问题时存在容易陷入局部最优 解的问题,对经典蚁群算法提出三方面改进.... 网络出版时间:2016-11-21 ...
- 迪杰斯特拉算法求经纬度坐标的最短路径_【图的最短路径】迪杰斯特拉算法求图的最短路径...
#include using namespace std; const int INFINITY=23678; const int M=3; /*typedef struct G { int ver[ ...
- java公网对讲_【对讲机的那点事】选择公网对讲机你必须要知道使用的网络信号!...
公网对讲机对于广大的用户来说已经不再陌生,公网对讲机的构成主要由硬件部分.IP网络运营平台.运营商提供的流量卡三大要素组成.利用移动通信的数据通道,将话音数字化,压缩,然后经现有的公众移动数据网络发送 ...
- java高斯分布随机数_生成符合高斯分布或者其他任意分布的随机数
在一些情况下经常需要用到随机数,而高斯随机数又是最常用到的.这一篇讲一下如何编程生成符合正态分布的高斯随机数,甚至任何其他分布的随机数. 我们知道C语言的标准库函数可以生成符合均匀分布的伪随机数.那么 ...
- Hdu1874 最短路径_畅通工程续
Problem Description 某省自从实行了很多年的畅通工程计划后,终于修建了很多路.不过路多了也不好,每次要从一个城镇到另一个城镇时,都有许多种道路方案可以选择,而某些方案要比另一些方案行 ...
- postgis+geoserver最短路径
postgis+geoserver最短路径 1 安装软件 2 数据预处理 3 操作postgres 4 操作Geoserver 5 计算最短路径 6 问题记录 1 安装软件 安装PostgreSQL与 ...
最新文章
- 第11章 UART 串口通信(手把手教你学51单片机pdf部分)
- excel设置单元格为文本
- 51CTO,什么情况?
- Python编程语言学习:python中与数字相关的函数(取整等)、案例应用之详细攻略
- 简单的ASP.NET无刷新分页
- git命令下载项目,上传android项目到github步骤,以及取消项目关联Git,设置git忽略文件
- datatable中某一列最小值_Asp.net中获取DataTable选择第一行某一列值
- 创业经验点滴 五(转) 关于淘宝
- [原] CentOS 7 安装 nginx, php mysql 套件
- 实战项目——小王优品铺
- Soap+xml实现webservice 调用
- 滑动窗口协议如何实现流量控制
- linux syslog原理,Linux日志管理之詳解 syslog/vsyslog
- mysql启动失败LSB_mysql启动报错:Failed to start LSB: start and stop MySQL
- 想成为牛逼网页设计师吗?
- QT报make: Circular all <- first dependency dropped.错误解决。
- 手撕自动驾驶算法——无迹卡尔曼滤波(UKF)
- 图文并茂的电脑组装过程讲解
- 计算机课程考试质量分析,计算机考试质量分析.pdf
- 软件构造期末复习总结【一】:关于SCM和Git的一些概念整理