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计算是通过后台程序实现的。

转载自:https://blog.csdn.net/GISShiXiSheng/article/details/49642745

You may also like...