PostGIS拓扑:TOPOLOGY

数据准备:建立架构waterdataset,导入样例供水管网数据shp文件,pipesectionmain管线表,fitting节点表。

 

1)创建与删除拓扑

创建拓扑,CreateTopology(topo名称,坐标系srid,容差值)

示例:通过现有空间数据表(读取现有空间数据的坐标系)创建拓扑

SELECT topology.CreateTopology('topo', find_srid('waterdataset','fitting','geom'),0.001);

删除拓扑,DropTopology(topo名称)

示例:删除建立的topo

SELECT topology.DropTopology('ma_topo');

示例:删除建立的TopoGeometry字段

语法,DropTopoGeometryColumn(架构名,表名,字段名);

SELECT topology.DropTopoGeometryColumn('waterdataset', 'pipesectionmain', 'topo_geom');

2)往空拓扑中加载数据

ST_CreateTopoGeo(已创建的topo名称,几何集合)

例:通过现有的空间表数据填充拓扑,该方法只在空topo的时候可以使用,导入过一次后该方法不可再用。

注:需导入线数据,导入数据后会填充edge_data表,并用线的端点填充node表。

SELECT topology.ST_CreateTopoGeo('topo',(SELECT ST_Collect(geom) FROM waterdataset.pipesectionmain) );

3)建立TopoGeometry列

可在原空间数据表上加入,也可以单独新建表加入。AddTopoGeometryColumn(topo名称,表空间名,表名,字段名,几何类型) ;

单独建新表示例:

//创建表空间

CREATE SCHEMA watertopo_feature;

//创建要素表

CREATE TABLE watertopo_feature.fitting(gid serial, fid integer PRIMARY KEY);

在原空间数据管线表添加示例:

SELECT topology.AddTopoGeometryColumn('hkgwater_topo', 'waterdataset', ' pipesectionmain ', 'topo_geom', 'LINESTRING');

4)给空间数据表中的TopoGeometry列填充数据

填充数据会自动建立起几何对象和topo对象的关系(会关联更新relation表)。

示例:更新管线表中的TopoGeometry列

toTopoGeom(空间数据表几何字段geom,自己建立的topo表空间,更新的layer_id,容差)

UPDATE waterdataset.pipesectionmain SET topo_geom = topology.toTopoGeom(geom, 'hkgwater_topo', 4, 0.001);

TopoGeometry存储方式,(topology_id,layer_id ,id, type)

topology_id:建立的拓扑的id,存储在topology表中

layer_id:建立的拓扑图层的id,存储在layer表中

id:对应空间数据表的gid

type:空间数据类型(1:[multi]point, 2:[multi]line, 3:[multi]poly, 4:collection),我们管网中常用的简单点为1,简单线为2。

 

5)重复3~4步骤,给fitting节点表建立TopoGeometry列并填充数据,fitting表会建立起与node表的对应关系。

 

参考资料

官方文档:http://postgis.net/docs/manual-dev/Topology.html#Topology_Processing

WIKI:http://trac.osgeo.org/postgis/wiki/UsersWikiPostgisTopology

示例:http://blog.mathieu-leplatre.info/use-postgis-topologies-to-clean-up-road-networks.html

posted @ 2019-09-04 16:19  GoodGF  阅读(2279)  评论(0编辑  收藏  举报