PostGIS 爆管分析实战记录之更新日志&&BUG修改

随着爆管分析运用场景越来越多，算法的问题也逐渐突显。在此新开一贴，专门记录算法的优化以及bug的更正（根据进程不断更新此贴）

问题一：寻找上游阀门性能优化

参考：https://www.cnblogs.com/giser-s/p/12091219.html

问题二：查找周围所有阀门：遍历死角

遇到这样的情况：爆点所在管段的source和target端都有阀门，但是算法给出的结果是另一个管段上的阀门。
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-2ol7f00K-1589766952086)(https://img2018.cnblogs.com/blog/365445/202002/365445-20200226145536128-1839740167.png)]
经检查发现：原来算法中，以source为起始点开始遍历，寻找与其相接管段（其他source或target等于爆管source的），由此找到了1、2管段

接着判断1、2管段上是否有阀门，有阀门就返回，没阀门就继续往下…

这里出现了问题，情况是：1、2管段被找出，理应查到a、b、c三个阀门，但是在查找管段上的阀门时，并没有循环所有被查出的阀门，故a被查丢了（其实丢a丢b是随机的）。
[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-ufmbfStJ-1589766952093)(https://img2018.cnblogs.com/blog/365445/202002/365445-20200226151643865-1645961507.png)]
出现这种情况的原因是，我在查找管段上阀门，考虑到广度遍历时每次循环只判断下游节点，只会有一个结果，因此没有对查出的阀门结果再循环。但是恰恰存在爆点所在管段上下都有阀门的情况（即第一次）。

所以打了个补丁，第一次的时候判断，如果爆管管段上下都有阀门，则直接返回，不再执行查找。

	  tap_count = 0;
	  FOR tap_temprow IN
	    select t.gid as v_uptap_gid,t.geom as v_uptap_geom from fm t where t.gid in (
          select a.gid from fm a,(select c.* from zy c where c.gid = v_startGid) b where ST_intersects(a.geom,b.geom) 
        ) 
	  LOOP
	    tap_count=tap_count+1;
	    --执行返回结果
        return query
        select tap_temprow.v_uptap_gid as res_uptap_gid,tap_temprow.v_uptap_geom as res_uptap_geom ;
	  END LOOP;
	  raise notice '%' , tap_count ;
	  IF(tap_count >=2) THEN
	    return;
	  END IF;

问题三：查找周围所有阀门：遍历时无法判断起始方向

遇到这样的一个问题，如果爆管点如下图黄点，返回的结果是3个点。从图中可以分析出，最上面那个点是错误阀门，不应该被返回

经调试发现，原来我查到离爆点最近的管段时，直接就去找与其source相接的管段了，爆点所在管段与他相接的上下游管段一起成为判断条件了。
故在遍历之前，先判断爆点所在管段两端是否有阀门，有阀门的话从另一侧开始遍历

 	  FOR tap_temprow_source IN 
		 select a.gid,a.geom,count(a.gid) as cnt from fm a,(select c.* from zy c where (c.source = v_startSource or c.target = v_startSource) AND c.gid <> v_startGid) b where ST_intersects(a.geom,b.geom) group by a.gid
 	  LOOP
		IF(tap_temprow_source.cnt > 0) THEN
		  IF (select ST_intersects(tap_temprow_source.geom,c.geom) from zy c where c.gid = v_startGid) THEN
            SELECT array_append(v_up_where, v_startTarget) into v_up_where;  
		    --v_up_where = v_startTarget;
		    tap_count_source = tap_count_source + 1;
		  END IF;
		END IF;
 	  END LOOP;
	  
 	  FOR tap_temprow_target IN 	 
		 select a.gid,a.geom,count(a.gid) as cnt from fm a,(select c.* from zy c where (c.source = v_startTarget or c.target = v_startTarget) AND c.gid <> v_startGid) b where ST_int