Geometry methods

Question

When I first saw this in 10 I though that would be very cool to do some spatial processing...
Never had a chance to try it until a few days ago. As usual I'm stuggling.

I have done a whole load of searches here and found plenty of discussion re the accuracy of the integer version of the geometry so that properties are incorrectly reported (like Poly1Geom.area)
But no real info on how to do the geometry comparisons. The help files are a little thin here as well.

My idea was to load the geometry objects into a python list / dictionary or whatever, to facilitate looping and comparison with other geometry objects.
The lists load fine and the geometry objects in the list have the usual properties. But when I get to the methods they are always "False". ie Poly1.inside(Poly2) always = False, even when its not.

I am using Arc10 SP4 and all this is running python code in IDLE.

Anyone got any tips or pointers to some script examples to do this?

Thanks in advance,
Neil

NeilAyres · Answer

No worries,
and yes I do read the help files before posting here....

My approach was exactly that, use cursors to cycle through the relavent geometries.
However, I was reading the geometry objects into a python list first before testing with the various methods.

N

MarcinGasior · Answer

You access geometry using cursors and then reading SHAPE field.
Here's an example where if a polygon contains a point, some attribute from polygon layer is copied to point layer:...
updCursor = arcpy.UpdateCursor(pointLayer)
for updRow in updCursor:
    pnt = updRow.Shape

#copy field values from polygon to point layer
    searchCursor = arcpy.SearchCursor(polygonLayer)
    for row in searchCursor:
        poly = row.Shape
        if poly.contains(pnt):
            updRow.setValue(copyField, row.getValue(copyField))
            updCursor.updateRow(updRow)
...

FabianBlau · Answer

Read the examples, this should help:

http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//002z0000001t000000

http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Geometry/000v00000095000000/