Calculate Geometry - Custom Python Set Units via CSV

457
4
02-08-2021 11:26 AM
ChasCrandell
New Contributor II

ArcGIS Dektop 10.6.1

I have to repeat a process on many GDBs quarterly. Need to calculate spatial fields (e.g. ccordinateX, coordinateY, featureArea, featureAreaUom, featureLength, featureLengthUom, featurePerimeter, featurePerimeterUom) with units based on an input CSV list of feature classes. Feature classes have spatial unit attributes that may differ based on the input CSV. CSV sample:

Feature_ClassAreaUOM_AreaLineUOM_Line
Grid_AACRESAcresMILES_USUS Survey Mile
Ber_LNANAFEET_USFeet
RegA_ASQUARE_YARDSSquare YardsFEE_USFeet
AirDot_PNANANANA

 

I merged some existing scripts (I'm still a novice with python) to get the script below, which works (except I have to build out for more fields), except it sets the featureArea and featureLength units to meters based on the coordinate system of the feature classes, rather than taking units from the CSV. There is a commented out line close to the bottom "row[count] = SHAPE.LENGTH@r[3]" that is the area where I think this needs to be altered?

 

import arcpy, os, sys, shutil, csv

Repository = arcpy.GetParameterAsText(0)
InCSV = arcpy.GetParameterAsText(1)

arcpy.env.workspace = Repository

#Create list of feature classes
def getFCList(vPath):
    vFullFCList = []
    vDSList = arcpy.ListDatasets()
    for ds in vDSList:
        vFCList = arcpy.ListFeatureClasses('*', 'All', ds)
        vFullFCList = vFullFCList + vFCList
    vFCList2 = arcpy.ListFeatureClasses()
    vFullFCList = vFullFCList + vFCList2
    return vFullFCList

#Get list of feature classes
fcList = getFCList(arcpy.env.workspace)
#List of fields to be checked
geoFields = [ "coordinateX", "coordinateY", "featureLength", "featureLengthUom", "featureArea", "featureAreaUom", "featurePerimeter", "featurePerimeterUom"]

for fc in fcList:
    fields = []
    desc = arcpy.Describe(fc)
    dataSet = desc.path
    arcpy.AddMessage("Updating Attributes in {0}".format(fc))
    fieldList = arcpy.ListFields(fc)
    #If the following fields exist, adds the field name to a list to be used later.
    for field in fieldList:
            if field.name == "coordinateX":
                fields.append(field.name)
                fields.append("SHAPE@X")                    
            elif field.name == "coordinateY":
                fields.append(field.name)
                fields.append("SHAPE@Y") 
            elif field.name == "featureLength":
                fields.append(field.name)
                fields.append("SHAPE@LENGTH")
            elif field.name == "featureArea":
                fields.append(field.name)
                fields.append("SHAPE@AREA")
            elif field.name == "featurePerimeter":
                fields.append(field.name)
                fields.append("SHAPE@LENGTH")
            else:
                fields.append(field.name)

    #The tool checks if the field exists within a feature class
    #If the field exists it will populate that field with its corresponding value
    arcpy.AddMessage(fields)
    with open(InCSV, 'r') as source:
        rdr = csv.reader(source)
        UpdateName = desc.name
        for r in rdr:
            if len(fields) > 0:
                with arcpy.da.UpdateCursor(fc, fields) as cursor:
                    for row in cursor:
                        count = 0
                        while count < len(fields):
                            if fields[count] == "coordinateX":
                                idx = fields.index("SHAPE@X")
                                row[count] = row[idx] 
                            elif fields[count] == "featureLength" and UpdateName == r[0]:
                                idx = fields.index("SHAPE@LENGTH")
                                row[count] = row[idx]
                                #row[count] = SHAPE.LENGTH@r[3]
                            elif fields[count] == "featureAreaUom" and UpdateName == r[0]:
                                row[count] = r[2] 
                            elif fields[count] == "featureLengthUom" and UpdateName == r[0]:
                                row[count] = r[4]   
                            count += 1

                        cursor.updateRow(row)
sys.exit()

 

 

0 Kudos
4 Replies
DavidPike
MVP Frequent Contributor

Personally I would find it simpler to pass these into an existing GP tool such as Calculate Geometry Attributes (Data Management)—ArcGIS Pro | Documentation 

You could just loop through the for r in rdr, then use r[0], r[1]... etc as arguments to the tool.

ChasCrandell
New Contributor II

Thanks. Probably beyond my ability and time, but I see what I can do with this if something else doesn't come up.

0 Kudos
DavidPike
MVP Frequent Contributor

Your csv is well set up for this, this is an example of the process:

with open(InCSV, 'r') as source:
    rdr = csv.reader(source)
    for r in rdr:
        fc_name = r[0]
        area_unit = r[2]
        length_unit = r[4]
        geometry_fields = []
        if area_unit != 'NA':
            geometry_fields.append('AREA')
        if length_unit != 'NA':
            geometry_fields.append('LENGTH')

        arcpy.management.AddGeometryAttributes(fcName, geometry_fields, length_unit, area_unit)
0 Kudos
JeffK
by MVP Regular Contributor
MVP Regular Contributor

I don't know what IDE you are using but take some time to learn the debugger.  Step through your code and look at what the variables are and adjust as necessary.

0 Kudos