ArcSDE vs. Oracle Spatial 13
I. 空间关系运算
以上两个章节是必不可少的基础，涉及到空间数据在数据库中的存储和通过索引加速空间数据的获取。从这句话也可以看出，后续的章节并不是必须的，也就是说某些基于空间数据的应用可能并不需要诸如空间关系判断、几何对象处理等功能。这并不是说这些功能就用不着了，而是这些功能并不一定需要在数据库端执行。比如基于ArcSDE，这些空间算法和功能在ArcGIS的产品线中无处不在，很多时候都不会把这些功能放到数据库上去。
· Oracle Spatial
在Oracle Spatial中，主要的空间关系操作在下表中列出：
| 空间操作 | 描述 |
| SDO_FILTER | 主过滤，判断哪些几何对象可能和给定的几何对象相交 |
| SDO_JOIN | 基于一定空间关系的表连接 |
| SDO_NN | 查询离某个几何对象最近的几何对象 |
| SDO_NN_DISTANCE | 查询离某个几何对象最近的几何对象与当前对象的距离 |
| SDO_RELATE | 判断两个几何对象是否满足某种空间关系 |
| SDO_WITHIN_DISTANCE | 判断两个几何对象间距离是否小于某给定值 |
其中SDO_RELATE操作比较特殊，它通过mask参数来判断几何对象间是否满足某种空间关系，但是事实上如果都用参数来指定空间关系的话使用上就不太方便，因此，Oracle Spatial还提供了基于SDO_RELATE不同参数代表的不同空间关系的一些操作，如下表：
| 衍生自SDO_RELATE的空间操作 | 描述 |
| SDO_ANYINTERACT | 任意部分有相交 |
| SDO_CONTAINS | |
| SDO_COVEREDBY | |
| SDO_COVERS | |
| SDO_EQUAL | |
| SDO_INSIDE | |
| SDO_ON | |
| SDO_OVERLAPBDYDISJOINT | |
| SDO_OVERLAPBDYINTERSECT | |
| SDO_OVERLAPS | |
| SDO_TOUCH |
下面我们看一些空间关系操作的大概用法，比如我用一个三角形去查询所有满足“INSIDE”关系的几何对象，我们可以通过以下两种等价的方式操作：
SQL> select continent from continent where sdo_inside(geom, SDO_GEOMETRY(2003,NULL,NULL,SDO_ELEM_INFO_ARRAY(1,1003,1),SDO_ORDINATE_ARRAY(0,0,180,-90,180,90,0,0)))='TRUE';
CONTINENT
---------------------------------------
Australia
SQL> select continent from continent where sdo_relate(geom, SDO_GEOMETRY(2003,NULL,NULL,SDO_ELEM_INFO_ARRAY(1,1003,1),SDO_ORDINATE_ARRAY(0,0,180,-90,180,90,0,0)), 'mask=INSIDE')='TRUE';
CONTINENT
---------------------------------------
Australia
图 13 判断和图示三角形满足一定空间关系的几何对象
虽然有一些比较方便的SDO_RELATE的衍生操作，不过SDO_RELATE支持更灵活的空间关系判断，比如以下的SQL语句可以查询和上面的三角形满足“OVERLAPBDYDISJOINT”关系、“INSIDE”关系及同时满足这两种关系的几何对象：
SQL> select continent from continent where sdo_relate(geom, SDO_GEOMETRY(2003,NULL,NULL,SDO_ELEM_INFO_ARRAY(1,1003,1),SDO_ORDINATE_ARRAY(0,0,180,-90,180,90,0,0)), 'mask=OVERLAPBDYDISJOINT')='TRUE';
CONTINENT
---------------------------------------
North America
SQL> select continent from continent where sdo_relate(geom, SDO_GEOMETRY(2003,NULL,NULL,SDO_ELEM_INFO_ARRAY(1,1003,1),SDO_ORDINATE_ARRAY(0,0,180,-90,180,90,0,0)), 'mask=INSIDE')='TRUE';
CONTINENT
---------------------------------------
Australia
SQL> select continent from continent where sdo_relate(geom, SDO_GEOMETRY(2003,NULL,NULL,SDO_ELEM_INFO_ARRAY(1,1003,1),SDO_ORDINATE_ARRAY(0,0,180,-90,180,90,0,0)), 'mask=INSIDE+OVERLAPBDYDISJOINT')='TRUE';
CONTINENT
---------------------------------------
North America
Australia