select by location and calculate field from update cursor row: help

Question

Hi,
I'm trying to use update.cursor to do a spatial selection of another layer, then populate the selection layer with an attribute from the cursor layer. More specifically, I'm trying to select "parcLyr" by each feature in "JurLyr", then populate "parcLyr" field "Dist" with the value in "JurLyr" field "Jur_Name".arcpy.MakeFeatureLayer_management(Jurs_dis, "JurLyr")
rows = arcpy.SearchCursor("JurLyr", "", "", "Jur_Name", "")
arcpy.MakeFeatureLayer_management(parcels, "parcLyr")
for row in rows:
    arcpy.SelectLayerByLocation_management("parcLyr", "HAVE_THEIR_CENTER_IN", "JurLyr", "", "NEW_SELECTION")
    arcpy.CalculateField_management("parcLyr", "Dist", "%s") %(row.Jur_Name)

I'm getting an "invalid character" error. I feel like I'm close.

Any help would be greatly appreciated.

Thanks!

Rich

JenniferHorsman · Answer

Thank you Kim for your first suggestion! I have a similar task and was trying to write a script, but when I read your suggestion, I gave my task some more thought and realized I could do what I needed using various selections, Spatial Join, and Join by Attributes with no need for a script except to put it altogether in one place. Whew!

KimOllivier · Answer

Forget any idea of iterating over each row. Better analysis tools do that for you.

Why not just use an overlay tool between the two featureclasses?
If you want to find if a parcel centroid is inside some other polygon featureclass then create a point featureclass of parcel centroids and use Identity_analysis.

Performance in my terms means finishing in preferably 30 seconds for something so simple, so you can run it several times until you get it right before morning tea. 😉 import arcpy
import datetime
start = datetime.datetime.now()
arcpy.env.overwriteOutput = True
arcpy.env.workspace = "e:/crs/customer/waiheke/mobile.gdb"
arcpy.FeatureToPoint_management("parcel", "in_memory/parcel_center","INSIDE")
arcpy.management.MakeFeatureLayer("E:/data/census2006/census2006.gdb/nztm/mb2006","jurs_lay","TA06V2 = '007'")
arcpy.Identity_analysis("in_memory/parcel_center", "jurs_lay", "in_memory/parcel_jurs")
# join back district to parcels
arcpy.JoinField_management("parcel", "objectid", "in_memory/parcel_jurs", "ORIG_FID", ["MB06"])
print "Done",datetime.datetime.now() - start, " seconds"

# >>> Done 0:00:26.953 seconds

RichardThurau · Answer

kimo
Thanks for the reply.
I'm not having any performance issues, which, once I get this type of script to run, perhaps i will. At that point, I would submit a thread to fix.

Your other suggestions sound interesting, but a little beyond me. I'm trying to learn all I can, but at the same time not get fired for spending too much time learning and not enough doing.

If I am "identifying" parcels that are within a spatial domain of another selected layer, it seems redundant to create a 'dictionary' of this relationship anywhere except in the parcels attribute table. This is a very simple process to do in ArcMap, but as i said, trying to learn.

Thanks to another suggestion, I am able to get the code to run using the following:arcpy.MakeFeatureLayer_management(Jurs_dis, "JurLyr")
rows = arcpy.SearchCursor("JurLyr", "", "", "Jur_Name", "")
arcpy.MakeFeatureLayer_management(parcels, "parcLyr")
for row in rows:
    arcpy.SelectLayerByLocation_management("parcLyr", "HAVE_THEIR_CENTER_IN", "JurLyr", "", "NEW_SELECTION")
    arcpy.CalculateField_management("parcLyr", "Dist_amec", '"%s" %(row.Jur_Name)', "PYTHON")

However, this code populates all the parcels with the first jurName in the JurLyr. Instead of just populating the correct parcels within each jurName. Since there is one feature for each jurName, i assumed the updatecursor would loop through the attribute table. But, maybe I need something to make this iterate?

Thanks for any input you can provide.

Rich