ModelBuilder how to compare several fields between 2 FCs
by
1 Solution
Accepted Solutions
by
by
by
XanderBakker
To help you a bit, have a look at the tool and code below:
Parameters:
- 0: FeatureLayer (1)
- 1: Field (multivalue, derived from featurelayer1)
- 2: Featurelayer (2)
- 3: String (direction output) additions
- 4: String (direction output) deletes
What you choose to do with the resulting lists defines how they should be defined and formatted as output.
Result:
Code:
def main():
import arcpy
fc1 = arcpy.GetParameterAsText(0)
selected_flds = arcpy.GetParameterAsText(1)
fc2 = arcpy.GetParameterAsText(2)
flds = selected_flds.split(';')
# check if selected fields in fc2
arcpy.AddMessage('Check if selected fields in fc2')
for fld in flds:
if len(arcpy.ListFields(fc2, fld)) == 0:
arcpy.AddError('ERROR: Field {0} not in {1}'.format(fld, fc2))
fld_oid1 = arcpy.Describe(fc1).OIDFieldName
fld_oid2 = arcpy.Describe(fc2).OIDFieldName
flds1 = [fld_oid1] + flds
flds2 = [fld_oid2] + flds
# validation: are the three fields combined unique?
arcpy.AddMessage('Validation: are the three fields combined unique?')
if not isUnique(fc1, flds1):
arcpy.AddError("The combination of the 3 field for FC1 is not unique")
if not isUnique(fc2, flds2):
arcpy.AddError("The combination of the 3 field for FC2 is not unique")
# create dictionary
arcpy.AddMessage('Create dictionaries')
dct1 = generateDict(fc1, flds1)
dct2 = generateDict(fc2, flds2)
# create sets
arcpy.AddMessage('Create sets')
set1 = set(dct1.keys())
set2 = set(dct2.keys())
# check additions (in FC2, but not in FC1)
arcpy.AddMessage('Check additions (in FC2, but not in FC1)')
adds = [dct2 for a in list(set2 - set1)]
if len(adds) > 0:
arcpy.AddMessage('OIDs added: {0}'.format(adds))
arcpy.SetParameterAsText(3, adds)
# check deletes (in FC1, but not in FC2)
arcpy.AddMessage('Check deletes (in FC1, but not in FC2)')
dels = [dct1 for a in list(set1 - set2)]
if len(dels) > 0:
arcpy.AddMessage('OIDs deleted: {0}'.format(dels))
arcpy.SetParameterAsText(4, list(dels))
def generateDict(fc, flds):
return {u'{0}#{1}#{2}'.format(r[1], r[2], r[3]): r[0] for r in arcpy.da.SearchCursor(fc, (flds))}
def isUnique(fc, flds):
lst1 = [u'{0}#{1}#{2}'.format(r[1], r[2], r[3]) for r in arcpy.da.SearchCursor(fc, (flds))]
return len(lst1) == len(list(set(lst1)))
if __name__ == '__main__':
main()
