Here's the approach I was experimenting with that uses a spatial join (and field mapping).
points = "GeoWellTemp"
parcels = "WellTaxlots"
pointsFlds = ['OBJECTID']
parcelFlds = ['WellID', 'W_DEPTH', 'Geothermal', 'GeoWaterUs', 'TEMP', 'TEMP_DATE']
fieldmappings = arcpy.FieldMappings()
for f in pointsFlds:
fm = arcpy.FieldMap()
fm.addInputField(points, f)
fieldmappings.addFieldMap(fm)
for f in parcelFlds:
try:
fm = arcpy.FieldMap()
fm.addInputField(points, f)
outname = fm.outputField
outname.name = outname.name + "_99"
outname.aliasName = outname.aliasName + "_99"
fm.outputField = outname
fieldmappings.addFieldMap(fm)
except:
pass
for f in parcelFlds:
fm = arcpy.FieldMap()
fm.addInputField(parcels, f)
fieldmappings.addFieldMap(fm)
join = 'in_memory\sj_test'
if arcpy.Exists(join):
arcpy.Delete_management(join)
arcpy.SpatialJoin_analysis(target_features = points, join_features = parcels, out_feature_class = join,
join_operation="JOIN_ONE_TO_MANY", join_type="KEEP_ALL", field_mapping = fieldmappings,
match_option="INTERSECT", search_radius="", distance_field_name="")
joinFlds = pointsFlds + parcelFlds
valueDict = {r[0]:(r[1:]) for r in arcpy.da.SearchCursor(join, joinFlds)}
with arcpy.da.UpdateCursor(points, joinFlds) as updateRows:
for updateRow in updateRows:
keyValue = updateRow[0]
if keyValue in valueDict:
for n in range (1,len(joinFlds)):
updateRow[n] = valueDict[keyValue][n-1]
updateRows.updateRow(updateRow)
if arcpy.Exists(join):
arcpy.Delete_management(join)
del valueDict
The field mappings took some experimenting. I first tried a spatial join and examined the code the tool creates (copy as python snippet). This provided a way to examine the default field map the tool generates. The points and parcels layers had a number of fields that used the same name. Normally, the target feature (points, in this case) would use the original field name, and the join feature (parcels) would append an underscore with a number to rename any fields sharing a name in the target feature. As I wanted the original field names for the join feature, I renamed the target feature's fields when setting up the field mapping. I think the identical names caused some issues for the spatial join tool, as it would use the values of the target feature when I selected only the fields from the join feature.
I ended up performing three steps to assemble the field map. First, I selected the object id field from the target feature to be used for linking to the points feature. The second step was to select and rename any fields in the target feature that had identical names in the join feature that were to be used. The final step was to select the fields in the join feature that were used to update the points feature.
I used an in_memory location for the spatial join; it needed to be a temporary join anyway. I used a 'one to many' join to have access to the JOIN_FID field that is created which would assist in linking the two layers together for any debugging necessary.
The last steps were to load the join layer's data into a dictionary and to use it in an update cursor.
I think this approach is a bit faster than the one in my previous post. But it took much more time to get debugged and working. Hope this helps.
