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：对象是不相交的。

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

所以我的问题是：

这是一个错误吗？ contains参数显然比anyinteract
检查集合中的任何点是否位于多边形中的最佳方法是什么？包含在这里没有帮助，因为如果任何一点位于外面，则结果od SDO_RELATE为false。没有比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;

测试表包含多边形的三种变体。让我们尝试一下查询：

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

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

        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搜索返回那些。 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不适用于2d地理坐标中的点数

1 个答案: