SDO_RELATE ANYINTERACT 不适用于二维地理坐标中的点

Question

最初的问题是，我试图找出集合中的一个点是否在给定的多边形内。我没有得到任何结果，所以我把它简化为一个最简单的案例，但仍然没有结果。

多边形几何形状如下（坐标中心周围的小方块）：

geom1 := SDO_GEOMETRY(2003, 8307, null, SDO_ELEM_INFO_ARRAY(1,1003,1),     
SDO_ORDINATE_ARRAY(10,10, 10,-10, -10,-10, -10,10))

现在我试图找出点 [0,0] 是否在使用 SDO_RELATE 的多边形内部。

sdo_relate(geom1, 
        SDO_GEOMETRY('MULTIPOINT((0 0))',8307),
       'mask=anyinteract')

我用过anyinteract，因为根据Oracle documentation：

ANYINTERACT: The objects are non-disjoint.

多边形内的点显然与其不相交。所以我花了大约一个小时尝试以不同的方式初始化点并检查坐标和所有内容，然后我尝试将 CONTAINS 而不是 ANYINTERACT 放在那里并最终获得所需的输出。

所以我的问题是：

这是一个错误吗？ contains 参数显然比 anyinteract
检查集合中的任何点是否位于多边形内的最佳方法是什么？包含在这里没有帮助，因为如果任何点位于外部，则结果 od SDO_RELATE 为假。没有比 anyinteract 更合适的交集类型了，但这是行不通的。

Answer 1

首先你的多边形是畸形的：它没有关闭。对于多边形，最后一个点必须与第一个点匹配。您可以通过执行以下操作来检测错误：

SQL> select sdo_geom.validate_geometry_with_context (SDO_GEOMETRY(2003, 8307, null, SDO_ELEM_INFO_ARRAY(1,1003,1),SDO_ORDINATE_ARRAY(10,10, 10,-10, -10,-10, -10,10)), 0.05) from dual;

13348 [Element <1>] [Ring <1>]    
1 row selected.

ORA-13348: polygon boundary is not closed

一旦你更正它，你会得到一个不同的错误：

SQL> select sdo_geom.validate_geometry_with_context (SDO_GEOMETRY(2003, 8307, null, SDO_ELEM_INFO_ARRAY(1,1003,1),SDO_ORDINATE_ARRAY(10,10, 10,-10, -10,-10, -10,10, 10,10)), 0.05) from dual;

13367 [Element <1>] [Ring <1>]
1 row selected.

ORA-13367: wrong orientation for interior/exterior rings

对于多边形，点必须处于逆时针方向。更正后，多边形就正确了：

SQL> select sdo_geom.validate_geometry_with_context (SDO_GEOMETRY(2003, 8307, null, SDO_ELEM_INFO_ARRAY(1,1003,1),SDO_ORDINATE_ARRAY(10,10, -10,10, -10,-10, 10,-10, 10,10)), 0.05) from dual;

TRUE
1 row selected.

并且查询工作正常。这是一个例子：

create table t1 (id number, note varchar2(20), geom sdo_geometry);

insert into t1 (id, note, geom) 
values (
 1,
 'Not closed',
 SDO_GEOMETRY(2003, 8307, null, SDO_ELEM_INFO_ARRAY(1,1003,1), SDO_ORDINATE_ARRAY(10,10, 10,-10, -10,-10, -10,10))
);

insert into t1 (id, note, geom) 
values (
 2,
 'Wrong orientation',
 SDO_GEOMETRY(2003, 8307, null, SDO_ELEM_INFO_ARRAY(1,1003,1), SDO_ORDINATE_ARRAY(10,10, 10,-10, -10,-10, -10,10, 10,10))
);

insert into t1 (id, note, geom) 
values (
 3,
 'Valid',
 SDO_GEOMETRY(2003, 8307, null, SDO_ELEM_INFO_ARRAY(1,1003,1), SDO_ORDINATE_ARRAY(10,10, -10,10, -10,-10, 10,-10, 10,10))
);

commit;

insert into user_sdo_geom_metadata (table_name, column_name, diminfo, srid)
values (
  'T1',
  'GEOM', 
  sdo_dim_array (
    sdo_dim_element ('Long',-180,180,0.5),
    sdo_dim_element ('Lat',-90,90,0.5)
  ),
  8307
);
commit;

create index t1_sx on t1 (geom) indextype is mdsys.spatial_index;

测试 table 包含多边形的三个变体。让我们尝试查询：

select *
from t1
where sdo_relate(
  geom, 
  SDO_GEOMETRY('MULTIPOINT((0 0))',8307),
  'mask=anyinteract'
) = 'TRUE';

哪个 returns 正确的结果：只有有效的多边形被正确识别为包含您的点：

        ID NOTE                 GEOM(SDO_GTYPE, SDO_SRID, SDO_POINT(X, Y, Z), SDO_ELEM_INFO, SDO_ORDINATES)
---------- -------------------- -------------------------------------------------------------------------------
         3 Valid                SDO_GEOMETRY(2003, 8307, NULL, SDO_ELEM_INFO_ARRAY(1, 1003, 1), SDO_ORDINATE_ARRAY(10, 10, -10, 10, -10, -10, 10, -10, 10, 10))

1 row selected.

至于多边形中点的 ANYINTERACT 和 CONTAINS 之间的区别，它会影响那些 在多边形边界 上的点（或者更准确地说在边界）。 ANYINTERACT 搜索 returns 那些。 CONTAINS/INSIDE 搜索没有。

这可能很重要，例如，当您想要将点分布到多边形（如销售区域中的客户点位置）并且一些客户位置点恰好落在相邻区域之间的边界上时：ANYINTERACT 搜索将报告那些地点位于两个地区。 INSIDE 搜索将报告它们不在中。

最后，使用更简单的查询语法：

select *
from t1
where sdo_anyinteract(
  geom, 
  SDO_GEOMETRY('MULTIPOINT((0 0))',8307)
) = 'TRUE';

有关如何获取与错误代码（如 13348）对应的消息的一些信息。我通常使用这种技术：

set serveroutput on
exec dbms_output.put_line(sqlerrm(-13348))

这将打印出完整的错误消息。请注意，您必须传递带有减号的错误。

SDO_RELATE ANYINTERACT 不适用于二维地理坐标中的点

SDO_RELATE ANYINTERACT doesn't work with points in 2d geocoordinates

oracle

oracle-spatial