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