通过参看文献,选择门要素、房间要素中点、不通行要素中点等作为节点在ArcGIS中创建点要素
连接不同点要素在ArcGIS中创建线要素,并删除穿越过无法通行要素的路径。
在ArcGIS中通过使用拐点分割线的工具将线要素进行切割。并生成最终要使用的路径图层
路网数据的构建分为两种,一种在ArcGIS中构建、一种在postGIS中构建。
将点要素与线要素放在同一文件夹,在线要素上右键选择创建网络数据集。
创建postGIS拓展
使用SHP2SQL工具,将shp导入数据库
其中
点击“View connection details…”在弹出的窗口中填入PostgreSQL的账号和密码,以及Database。
连接成功后,需要设置一下“Optionns…”
需要使用GBK编码,并勾选最下面一个选项。
添加路网数据,并设置SRID为:4326
(注意路径不能出现中文)
--添加起点id
ALTER TABLE public.point2line ADD COLUMN source integer;
--添加起点id
ALTER TABLE public.point2line ADD COLUMN target integer;
--添加道路权重值
ALTER TABLE public.point2line ADD COLUMN length double precision;
--添加路的起终点直角坐标坐标
alter table public.point2line add column source_x double precision;
alter table public.point2line add column source_y double precision;
alter table public.point2line add column target_x double precision;
alter table public.point2line add column target_y double precision;
--为sampledata表创建拓扑布局,即为source和target字段赋值,其中0.000001为容忍度,大尺度小范围的图需要更低的容忍度
SELECT pgr_createTopology('public.point2line',0.0000001, 'geom', 'gid');
--为source和target字段创建索引
CREATE INDEX source_idx ON point2line("source");
CREATE INDEX target_idx ON point2line("target");
--为length赋值
update point2line set length =st_length(geom);
--为road_xblk表添加reverse_cost字段并用length的值赋值
ALTER TABLE point2line ADD COLUMN reverse_cost double precision;
UPDATE point2line SET reverse_cost =length;
update public.point2line set source_x =st_x(ST_startpoint(geom));
update public.point2line set source_y =ST_Y(ST_STARTpoint(geom));
update public.point2line set target_x =st_x(ST_ENDpoint(geom));
update public.point2line set target_y =ST_Y(ST_endpoint(geom));
-- 使用dijkstra查询最短路径,
select b.edge as gid,
(select the_geom from point2line_vertices_pgr where id=b.node) as node_geom
, b.node ,(select geom from point2line where gid=b.edge) as geom
from pgr_dijkstra(
'select gid as id, source, target, length as cost FROM point2line',2, 11,false
) as b;
从点2-11,并将结果中的节点与边,作为条件,对节点表和边表进行查询,获得geom并可视化。
select b.edge as gid,agg_cost,
(select the_geom from point2line_vertices_pgr where id=b.node)
as node_geom, b.node ,(select geom from point2line where gid=b.edge)
as geom from pgr_astar( 'select gid as id, source, target,
length as cost ,source_x as x1,source_y as y1,
target_x as x2,target_y as y2 FROM point2line',
80, 74,false,2 ) as b;
大概是 使用了起始点的方位作为启发函数
QGIS中layer菜单,选择对应数据库,输入密码,点击连接测试成功后
点击连接获得数据库中图层
选择要导入的图层。
获得图层后,选择数据库菜单database
选择postgis,按F2调出sql操作界面
在执行框中执行查询语句。并将数据中包含的geon数据导出至QGIS界面。
导入矢量文件,六边形格网
调整样式