Oracle存储过程结合ST_Geometry实现根据空间关系建立管点管线属性关联

一、问题背景

  某客户单位提供了一批管点管线GIS数据，每个管点具有一个点号字段（EXP_NO，字符类型），管线则具有起始管点、终止管点字段（S_EXP、E_EXP，字符类型。只作区分，不具方向性）。
   预处理已经为每个管点生成了点号。管线数据在空间上也大多与管点之间是有关联的，即管点刚好位于管线段的起始或终止点位置。
  现在需要对每一条管线段数据进行处理，将管线段端点上的管点的点号字段的值取出来，存放到管线段的起始管点、终止管点字段中。
  本例中，使用的范围面数据来自pjrnt图层，管线段数据来自图层piplk，管点数据来自图层mhhpt。图中黄色边框红色阴影部分为所使用的范围面（其objectid=641），蓝色线段即为管线段，黄色黑心加点为管点数据。如图可知，管点具有点号值，但管线的起始点号、终止点号均为null。

二、解决思路

  1. 选择工作区范围，进行包含查询，得到位于其中的所有管线段。
  2. 循环查询结果，取得每一条管线段的起点、终点、OBJECTID。
  3. 对起点、终点分别创建一个极小缓冲区。
  4. 在缓冲区内进行包含查询，取得位于缓冲区内的管点。
  5. 将管点的起点、终点的点号分别保存到管线段的起始点号、终止点号字段中。

三、所需工具

  PL/SQL（本例使用SYDXGH用户以管理员身份登录）
  Oracle 11g r2（包括服务器端、客户端）
 已经配置了ST_Geometry的ArcSDE环境
  已经入库的数据

四、实现代码

 DECLARE
TYPE cur IS REF CURSOR;
temp_cursor cur;
aShp sydxgh.piplk.shape%TYPE;
aOid integer;
BEGIN
OPEN temp_cursor
FOR
SELECT t1.shape AS aShp,t1.objectid AS aOid FROM sydxgh.piplk t1
WHERE sde.st_within(t1.shape,(SELECT t2.shape FROM sydxgh.pjrnt t2 WHERE t2.objectid=641))=1;
LOOP
FETCH temp_cursor INTO aShp,aOid;
EXIT WHEN temp_cursor%NOTFOUND;
BEGIN
UPDATE sydxgh.piplk t2 SET
t2.s_exp = (SELECT t.exp_no FROM sydxgh.mhhpt t WHERE sde.st_within(t.shape,(sde.st_buffer(sde.st_startpoint(aShp),0.01)))=1 AND ROWNUM=1),
t2.e_exp = (SELECT t.exp_no FROM sydxgh.mhhpt t WHERE sde.st_within(t.shape,(sde.st_buffer(sde.st_endpoint(aShp),0.01)))=1 AND ROWNUM=1)
WHERE t2.objectid=aOid;
END;
END LOOP;
CLOSE temp_cursor;
END;

代码说明：
  DECLAR中声明了一系列后面将要使用的变量，其中temp_cursor为查询位于工程范围面内的所有管线段所使用的游标，使用该游标可以获得每一条管段数据；
  aShp是游标中管线段的图形；
  aOid是游标中管线段的OBJECTID（用于确定将更新哪条管段的起始终止点号）。
  FOR循环部分，逐条取得管线段记录，使用了最简单的st_within查询，含义是查询位于OBJECTID=641的工程范围面内的管线段图层piplk内的数据。
  LOOP部分是对FOR循环的每条结果进行处理。
  FETCH行，从游标中取出一条记录的aShp和aOid。
  BEGIN和END中对查询结果进行更新，其中st_startpoint、st_endpoint可取得管线段的起点和终点，st_buffer方法建立缓冲区（半径为0.01米），外层的st_within方法查询位于缓冲区内的管线点，并取得其点号即exp_no。
执行完毕后结束循环，关闭游标。




结果验证：
 提交事务后，我们在ArcMap中检查发现，所有位于工程范围面内的管线段均已正确设置了起始点号、终止点号。

五、思考总结

  St_Geometry本身已经集成了强大的空间查询、处理功能，结合Oracle的存储过程，可以通过简单代码实现复杂的空间运算。有时候一些看似需要进行功能开发的需求，利用现成的工具可以方便地达到我们的目的，并且灵活高效。