POST
|
Hi again, I have just figured out that if I use the parameter spatial_grid_1=1 instead of spatial_grid_1=0 in my python script, my SQL Query works. So it seems that I need to set the grid values explicitly. Thanks a lot for helping me so far! Now I hope this will help me solving my original problem, which is a bit more complicated. - Morten
... View more
10-31-2013
01:05 AM
|
0
|
0
|
1231
|
POST
|
Hi Vince, Thanks for your reply. I have tried both with and without data => same problem As I understand a spatial index is automatically added. I do not add a spatial index explicitly. (also see console output below) My parameters are spatial_grid_1=0, spatial_grid_2=0, and spatial_grid_3=0, so I leave it to ArcGIS to determine the values. If I run the python code below:
deconn = "c:\\arcgisconfig\\oracle.sde"
prj_file_path = "c:/Program Files (x86)/ArcGIS/Desktop10.1/Reference Systems/World Time Zones.prj"
fc = arcpy.CreateFeatureclass_management(sdeconn, "MYFEATURECLASS", "POLYGON", spatial_reference=prj_file_path, spatial_grid_1=0, spatial_grid_2=0, spatial_grid_3=0)
array = arcpy.Array([arcpy.Point(45, 50), arcpy.Point(47, 50), arcpy.Point(47, 49), arcpy.Point(45, 50)])
polygon = arcpy.Polygon(array)
cursor = arcpy.da.InsertCursor(fc, ("SHAPE@"))
cursor.insertRow((polygon,))
del cursor
desc = arcpy.Describe(fc)
print "Feature Type: " + desc.featureType
print "Shape Type : " + desc.shapeType
print "Spatial Index: " + str(desc.hasSpatialIndex)
The output from my console is:
Feature Type: Simple
Shape Type : Polygon
Spatial Index: True
And I still get the error from my SQL Query. - Morten
... View more
10-31-2013
12:52 AM
|
0
|
0
|
1231
|
POST
|
Thanks Vince and FJ Yes, I could open a tech support incident. It seems I can reproduce the error by creating a simple polygon feature class using the python script below, and then afterwards run the SQL query. PYTHON:
import arcpy
import sys
sdeconn = "c:\\arcgisconfig\\oracle.sde"
prj_file_path = "c:/Program Files (x86)/ArcGIS/Desktop10.1/Reference Systems/World Time Zones.prj"
arcpy.CreateFeatureclass_management(sdeconn, "MYFEATURECLASS", "POLYGON", spatial_reference=prj_file_path, spatial_grid_1=0, spatial_grid_2=0, spatial_grid_3=0)
SQL:
SELECT COUNT(objectid) FROM MYFEATURECLASS WHERE sde.st_intersects(shape, sde.st_geomfromtext('POINT(-50.7 64.6)', 4326)) = 1
It makes no difference whether I include or exclude the "spatial_grid_1=0, spatial_grid_2=0, spatial_grid_3=0" parameters. - Morten
... View more
10-30-2013
11:57 AM
|
1
|
0
|
1231
|
POST
|
Sorry, It was a typing mistake. I meant Oracle version 11.2.0.1. I have also tested it on 11.2.0.3, since I suspected some bug fixes in the Oracle patch to be related to my problem. However it didn't make any difference. Thanks, Morten (boedtkjer=mbhansen)
... View more
10-29-2013
12:03 PM
|
0
|
0
|
1232
|
Title | Kudos | Posted |
---|---|---|
1 | 10-30-2013 11:57 AM |
Online Status |
Offline
|
Date Last Visited |
11-11-2020
02:23 AM
|