环境： Win10 ArcMap10.4（用于数据处理） postgresql9.4 postgis2.2.3 pgRouting2.3（postgresql插件）说明：做爆管分析的第一步，需要先将数据做拓扑处理，以下在构建拓扑数据成功的基础上继续（保证gid、source、target字段）步骤： 1、用户设定的爆管点（startx，starty），会与真实管网位置有差距，这里设置15米容差：假设爆管点距离15米内，最近的管网为爆点管段 execute 'select geom, source, target, ST_StartPoint(geom) as startpoint,ST_EndPoint(geom) as endpoint from ' ||tbl|| ' where ST_DWithin(geom,ST_Geometryfromtext(''point('|| startx ||' ' || starty ||')'',3857),15) order by ST_Distance(geom,ST_GeometryFromText(''point('|| startx ||' '|| starty ||')'',3857)) limit 1' into v_startLine, v_startSource ,v_startTarget, v_statpoint ,v_endpoint; 2、开始查找，从爆管开始 --循环第一次，将爆管放入查询条件，查询爆管上是否有阀门 IF(v_up_idx = 0) THEN SELECT array_append(v_up_where, v_startSource) into v_up_where; ELSE --v_up_where = null; END IF; --循环开始 FOR up_temprow IN select zy1.gid,zy1.source,zy1.target from zy zy1 where source = any(v_up_where) or target = any(v_up_where) LOOP --查询管网上的点 select t.gid,t.geom from fm t where t.gid in ( select a.gid from fm a,(select c.* from zy c where c.gid = up_temprow.gid) b where ST_intersects(a.geom,b.geom) ) into v_uptap_gid, v_uptap_geom; 3、如爆管没有阀门，则继续循环下一层级，往下找与爆管相接的管段；如果有阀门，则返回阀门，从下一次循环中剔除（不再找与他相接的管段） --如果没查找到阀门，则继续往下查 IF(v_uptap_gid is null) then --source去重，判断如果数组中已有，则不添加 IF (v_up_where @> ARRAY[up_temprow.source::integer] OR v_all_where @> ARRAY[up_temprow.source::integer]) THEN ELSE SELECT array_append(v_up_where,up_temprow.source) into v_up_where; SELECT array_append(v_all_where,up_temprow.source) into v_all_where; END IF; --target去重，判断如果数组中已有，则不添加 IF (v_up_where @> ARRAY[up_temprow.target::integer] OR v_all_where @> ARRAY[up_temprow.target::integer]) THEN ELSE SELECT array_append(v_up_where,up_temprow.target) into v_up_where; SELECT array_append(v_all_where,up_temprow.target) into v_all_where; END IF; ELSE 4、如果有阀门，则返回阀门gid和geom --执行返回结果 return query select v_uptap_gid as res_uptap_gid,v_uptap_geom as res_uptap_geom ; 5、附上全部存储过程 -- Function: test_getpoint5(character varying, double precision, double precision) -- DROP FUNCTION test_getpoint5(character varying, double precision, double precision); CREATE OR REPLACE FUNCTION test_getpoint5( IN tbl character varying, IN startx double precision, IN starty double precision) RETURNS TABLE(v_gid integer, v_res geometry) AS $BODY$ declare v_startLine geometry;--离起点最近的线 v_startTarget integer;--距离起点最近线的终点 v_startSource integer; v_statpoint geometry;--在v_startLine上距离起点最近的点 v_endpoint geometry;--在v_endLine上距离终点最近的点 v_up_source integer;--游标，记录是否有记录 v_up_idx integer;--记录遍历到多少层级 v_uptap_gid integer;--上游阀门gid v_uptap_geom geometry;--上游阀门要素 v_all_where integer[];--记录所有查询过的管段 v_up_where integer[];--where条件，将遍历到阀门的管段gid排除 up_temprow record ; test integer; begin --查询离起点最近的线 --3857坐标系 --找起点15米范围内的最近线 execute 'select geom, source, target, ST_StartPoint(geom) as startpoint,ST_EndPoint(geom) as endpoint from ' ||tbl|| ' where ST_DWithin(geom,ST_Geometryfromtext(''point('|| startx ||' ' || starty ||')'',3857),15) order by ST_Distance(geom,ST_GeometryFromText(''point('|| startx ||' '|| starty ||')'',3857)) limit 1' into v_startLine, v_startSource ,v_startTarget, v_statpoint ,v_endpoint; --查找上游阀门 v_up_idx = 0; v_up_source = 1; test = 0; SELECT array_append(v_up_where, v_startSource) into v_up_where; WHILE array_length(v_up_where,1) > 0 LOOP --游标归零 v_up_source = 0; --IF(v_up_idx = 0) THEN --SELECT array_append(v_up_where, v_startSource) into v_up_where; --SELECT array_append(v_up_where, v_startTarget) into v_up_where; --ELSE --v_up_where = null; --END IF; --记录层级 v_up_idx = v_up_idx + 1; --获取当前层级节点 FOR up_temprow IN select zy1.gid,zy1.source,zy1.target from zy zy1 where source = any(v_up_where) or target = any(v_up_where) LOOP test = test +1; --清空需要查的点 IF(v_up_source = 0) THEN v_up_where = null; END IF; --清空初始执行节点 v_startSource = 0; --标志执行有数据 v_up_source = 1; --查询管网上的点 select t.gid,t.geom from fm t where t.gid in ( select a.gid from fm a,(select c.* from zy c where c.gid = up_temprow.gid) b where ST_intersects(a.geom,b.geom) ) into v_uptap_gid, v_uptap_geom; --如果没查找到阀门，则继续往下查 IF(v_uptap_gid is null) then --source去重，判断如果数组中已有，则不添加 IF (v_up_where @> ARRAY[up_temprow.source::integer] OR v_all_where @> ARRAY[up_temprow.source::integer]) THEN ELSE SELECT array_append(v_up_where,up_temprow.source) into v_up_where; SELECT array_append(v_all_where,up_temprow.source) into v_all_where; END IF; --target去重，判断如果数组中已有，则不添加 IF (v_up_where @> ARRAY[up_temprow.target::integer] OR v_all_where @> ARRAY[up_temprow.target::integer]) THEN ELSE SELECT array_append(v_up_where,up_temprow.target) into v_up_where; SELECT array_append(v_all_where,up_temprow.target) into v_all_where; END IF; ELSE raise notice '%' , v_uptap_gid ||'---'||cast(test as text); --执行返回结果 return query select v_uptap_gid as res_uptap_gid,v_uptap_geom as res_uptap_geom ; END IF; --return next; END LOOP; END LOOP; end; $BODY$ LANGUAGE plpgsql VOLATILE STRICT COST 100 ROWS 1000; ALTER FUNCTION test_getpoint5(character varying, double precision, double precision) OWNER TO postgres; 结尾：这里本想将上游阀门和下游阀门分开，但是我们建立的拓扑中并没有方向，所以改成了查询出所有的影响阀门。后续将继续研究，把方向数据放进去，实现上游阀门、下游阀门、精确找出总阀门的功能

PostGIS 爆管分析之根据爆点找出所有影响阀门

vip168

环境：
Win10
ArcMap10.4（用于数据处理）
postgresql9.4
postgis2.2.3
pgRouting2.3（postgresql插件）

说明：
做爆管分析的第一步，需要先将数据做拓扑处理，以下在构建拓扑数据成功的基础上继续（保证gid、source、target字段）

步骤：
1、用户设定的爆管点（startx，starty），会与真实管网位置有差距，这里设置15米容差：假设爆管点距离15米内，最近的管网为爆点管段

execute 'select geom, source, target, ST_StartPoint(geom) as startpoint,ST_EndPoint(geom) as endpoint from ' ||tbl|| ' where ST_DWithin(geom,ST_Geometryfromtext(''point('|| startx ||' ' || starty ||')'',3857),15)
           order by ST_Distance(geom,ST_GeometryFromText(''point('|| startx ||' '|| starty ||')'',3857))  limit 1' 
           into v_startLine, v_startSource ,v_startTarget, v_statpoint ,v_endpoint;

2、开始查找，从爆管开始

      --循环第一次，将爆管放入查询条件，查询爆管上是否有阀门
      IF(v_up_idx = 0) THEN
        SELECT array_append(v_up_where, v_startSource) into v_up_where;
      ELSE
        --v_up_where = null;
      END IF;  
   

  --循环开始
  FOR up_







此内容 登录 或 注册 后可见