[Bug]使用ST_Geometry查询结果不准确的解决方法
今天有个用户咨询，说是使用ST_Geometry查询，主要是面状图层查询点状图层，当面状图层的面积达到数平方公里时，就会出现不该查询到的信息却查了出来。用户是使用了sde.st_intersects，而且问题出现在为了做前端分页，使用了RowNum，那么我就找了一个数据量特别大的点图层和同投影的面图层进行测试，发现了一个问题，也就是我通过桌面叠加和ST_Geometry查询不一致的情况出现。
相关环境：
ArcSDE10 SP5 64Bit
ArcGIS Desktop 10.1
Oracle 11.2.0.1 64Bit
点状图层：ccc 有93674个点
面状图层：myarea 有3个面
从上图我们可以明显的看出ObjectiD=18的面状要素包含9个点，但是我使用ST_Geometry查询：
SQL> select ccc.objectid from ccc,myarea where myarea.objectid=18 and sde.st_intersects(ccc.shape,myarea.shape)=1 ;
OBJECTID
----------
93767
93766
93769
93765
93768
93772
已选择6行。
明显有问题，因为投影一样，所以这个因素排除，使用空间SQL查询最大的问题就是索引，将点状要素类的空间索引删除掉
SQL> select index_name from user_indexes where table_name='CCC'; INDEX_NAME ------------------------------ SYS_IL0000075699C00028$$ R93_SDE_ROWID_UK
然后我们再进行空间查询
SQL> select ccc.objectid from ccc,myarea where myarea.objectid=18 and sde.st_intersects(ccc.shape,myarea.shape)=1 ;
OBJECTID
----------
93761
93763
93765
93766
93767
93768
93769
93772
93773
已选择9行。
可见，空间索引的问题，我测试了将该索引删除重建还是不行，我就将我的数据导出来，删除SDE里面的数据，再导入进去，还不行。
得出的结论：数据有问题，至于什么问题不详！
另外，我又将图上所有点导出一个单独的要素类数据
SQL> select objectid from ccd;
OBJECTID
----------
93764
93769
93771
93766
93768
93770
93767
93772
93773
93763
93761
OBJECTID
----------
93765
93762
已选择13行。同样的SQL语句
删除空间索引查询出来的结果也是9条
SQL> select index_name from user_indexes where table_name='CCD';
INDEX_NAME
------------------------------
R276_SDE_ROWID_UK
SYS_IL0000080306C00029$$
SQL> select ccd.objectid from ccd,myarea where myarea.objectid=18 and sde.st_intersects(ccd.shape,myarea.shape)=1
;
OBJECTID
----------
93769
93766
93768
93767
93772
93773
93763
93761
93765
已选择9行。如果添加上空间索引
SQL> select index_name from user_indexes where table_name='CCD';
INDEX_NAME
------------------------------
A241_IX1
R276_SDE_ROWID_UK
SYS_IL0000080306C00029$$
SQL> select ccd.objectid from ccd,myarea where myarea.objectid=18 and sde.st_intersects(ccd.shape,myarea.shape)=1
;
OBJECTID
----------
93766
93767
93765
93769
93768
93772
93773
93761
已选择8行。我已经晕了，就当是数据的原因吧，我是无语了............
可见ArcSDE的ST_Geometry的问题还是很多的.
-----------------------------------------------------------------
经过跟美国Esri的工程师的沟通，他们怀疑是一个Bug，不过可以将ArcCatalog自己生成的grid值进行自定义的修改，这样就可以得到正确的结果了。如果用户遇到类似的问题，可以有两种方式
1：重建一下空间索引
2：重新修改一下空间索引的Grid值
感兴趣的用户可以下载测试数据
最后裁定是一个bug
===================================
[#NIM086005 ST_Intersects function yields different results when using auto-calculated spatial index compared to no spatial index. ]
===================================
-------------------------------------------------------------------------------------------------------
版权所有，文章允许转载，但必须以链接方式注明源地址，否则追究法律责任!
------------------------------------------------------------------------------------------------------
Bug NIM-086005
| Bug Number | NIM086005 |
|---|---|
| Submitted | Oct 25, 2012 9:34 PM |
| Severity | Low |
| Applies To | ArcSDE/Enterprise Geodatabase |
| Version Found | 10.0 SP5 |
| Prog Language | N/A |
| Server Platform | |
| Client Platform | |
| Database | Oracle |
| Locale | N/A |
| Status | |
| Version Fixed | 10.1 |
| SP Fixed | N/A |
Synopsis
ST_Intersects function yields different results when using auto-calculated spatial index compared to no spatial index.
Additional Status Information
Behavior cannot be reproduced using 10.1.
Alternate Solution
N/A