sde用sql实现erase

概述：

本文讲述基于Arc SDE forOracle实现erase空间分析计算。

实现流程：

1、叠加计算

判断叠加，非叠加部分即为一部分所要结果，叠加部分进入第二步；

2、合并计算

根据objectid进行union计算；

3、差异计算

用原始数据data1和合并计算后的结果进行差异计算，所得的结果即为另一部分所要结果。

相关sde函数介绍：

1、sde.st_intersection(st_geometryshape1,  st_geometry shape2)

         语法：

参数：(st_geometry shape1,  st_geometryshape2)

                   返回值：st_geometry shape

         解释：

Returns a geometry that represents the shared portion of shape1 andshape2.

示意：

2、sde.st_union(st_geometryshape1,  st_geometry shape2)

         语法：

参数：(st_geometry shape1,  st_geometryshape2)

                   返回值：st_geometry shape

         解释：

Returns a geometry that represents the point set union of theGeometries.

         示意：

3、sde.st_difference(st_geometryshape1,  st_geometry shape2)

         语法：

参数：(st_geometry shape1,  st_geometryshape2)

                   返回值：st_geometry shape

         解释：

                   Returns a geometry thatrepresents that part of geometry A that does not intersect with geometry B.

         示意：

说明：

看到图3，很多人就开始问了：erase的效果不就是图3的效果吗，为什么还要那么多步？是的，对于两两的geometry来说，difference结果即为我们想要的erase的结果，但是，对于两个图层来说，difference后的结果是整个data2的结果，并且结果中会有一些叠加与重复。为什么会出现这样的结果呢，我想sde的算法实现中，也是两两做的计算，因此，计算的结果就是整个data2的结果。举例：A为data1的一个要素，B、C分别为data2的两个要素，在sde中，difference计算的结果为图4,1和2分别为计算后的结果，并重叠在一起显示，但是实际中，我们想要的及全国是图5。

图4

图5

实现sql实例：

1、计算叠加

insert intohx2q_risk_result_temp(objid,shape)

  select *

    from (select a.objectid,sde.st_intersection(a.shape, b.shape) as shape

            from hx2q_project a, hx2q_landnormb

           where a.ptcode = 'm01'

             and b.ptcode = 'm01'

             and a.status = '开工'

             and b.confirm = '已取得'

             and a.flag = 0

             and b.flag = 0)

   where sde.st_area(shape) > 0

2、计算union

3、计算difference

insert intoHX2Q_RISK_RESULT

  (objid, ptcode, Shape, objectid)

select a.objectid,

       a.ptcode,

       sde.st_difference(a.shape, b.shape) asshape,

       SDE.version_user_ddl.NEXT_ROW_ID('HX2Q','36') as objectid

  from hx2q_project a, hx2q_risk_result_temp b

 where a.objectid = b.objid

4、获取非叠加

insert intoHX2Q_RISK_RESULT

  (objid, ptcode, Shape, objectid)

select objectid,

       ptcode,

       shape,

       SDE.version_user_ddl.NEXT_ROW_ID('HX2Q','36')

  from hx2q_project

 where objectid not in (select distinct objidfrom hx2q_risk_result_temp)

 and ptcode = 'm01' and flag = 0 and status = '开工'

说明：

1、difference计算中，where后的条件是必须的；

2、union计算是通过后台程序实现的。