MapUnits Symbols Check script

[dvohra-igs]

This script was developed to make the user aware if any map units are symbolized by more than one unit, or shares a symbol in common with another unit. In such case, the user is told what the offending map units are along with their symbol codes. In my experience, it's been helpful for QA, particularly when dealing with compilation databases.

#This script checks whether or not each MapUnit has only one symbol code associated with it, and whether that symbol code is used by any other MapUnit.
#Idaho Geological Survey
#David Vohra
#[email protected]
#2024-10-30
import arcpy, os, sys

gdb = r'{}'.format(sys.argv[1])#geodatabase filepath goes here
feature_dataset = "GeologicMap"
lookup_table = "MapUnitPolys"
lookup_table_path = os.path.join(gdb, feature_dataset, lookup_table)

#You should not need to edit anything below this line.
def unique_values(table, field):
    with arcpy.da.SearchCursor(table, [field]) as cursor:
        return sorted({row[0] for row in cursor})

arcpy.env.OverwriteOutputs = True
arcpy.env.workspace = gdb
multisymbolmapunits = []
skippedmapunits = []
mapunithistorylog = []
sharedsymbolmapunits = []
myMapUnits = unique_values(lookup_table_path , 'MapUnit')#list of MapUnits present in the MapUnitPolys feature class.
for item in myMapUnits:
    #arcpy.SelectLayerByAttribute_management(lookup_table_path, "NEW_SELECTION", "MapUnit = '{}'".format(item))
    mup_lyr = "MUP_{}_lyr".format(item)
    arcpy.management.MakeFeatureLayer(lookup_table_path, mup_lyr, "MapUnit = '{}'".format(item))
    mySymbols = unique_values(mup_lyr, 'Symbol')#returns a list of all symbols for that one MapUnit. Hopefully the list has just one entry in it.
    sharedsymbollog = []
    for mupsymbol in mySymbols:
        if mapunithistorylog.count(mupsymbol) > 0: #if this symbol has already been encountered before,
            sharedsymbolmapunits.append(item) #log the mapunit in a permanent list for reporting later.
            sharedsymbollog.append(mupsymbol) #Also log the symbol in a temporary list specific to this mapunit.
        mapunithistorylog.append(mupsymbol) #keep a permanent log of which symbols the script has seen before. This is needed to determine if a mapunit shares the same symbol as another mapunit.
    if len(sharedsymbollog) > 0: #if the mapunit shares symbols with any other map unit,
        sharedsymbolmapunits.append(sharedsymbollog) #log the offending symbols into a permanent list for reporting later.
    if len(mySymbols) > 1: #If there are multiple symbols assigned to the same map unit,
        multisymbolmapunits.append(item) #log the mapunit, then
        multisymbolmapunits.append(mySymbols) #log the list of symbols for that mapunit.
    arcpy.Delete_management(mup_lyr)

if len(multisymbolmapunits) > 0:
    print("Map units with multiple symbols:" + str(multisymbolmapunits))
    arcpy.AddMessage("Map units with multiple symbols:" + str(multisymbolmapunits))
else:
    print("No map units assigned multiple symbols!")
    arcpy.AddMessage("No map units assigned multiple symbols!")

if len(sharedsymbolmapunits) > 0:
    print("Map units with non-unique symbols: " + str(sharedsymbolmapunits))
    arcpy.AddMessage("Map units with non-unique symbols: " + str(sharedsymbolmapunits))
else:
    print("All map units have unique symbol codes!")
    arcpy.AddMessage("All map units have unique symbol codes!")

if len(skippedmapunits) > 0:
    print("Map units skipped: " + str(skippedmapunits))
    arcpy.AddMessage("Map units skipped: " + str(skippedmapunits))