Latest Contributions by BlakeTerhune

After getting a lot of requests to get latitude and longitude coordinates for features, I created a generic script that creates two new fields in the existing feature feature class and and calcs the latitude and longitude of the feature's "true centroid." The only problem is that you would have to run this script any time a feature's geometry was changed. I did some searching for something that already existed and didn't find anything that was simple and generic enough, so please let me know if you can do it better. import arcpy fc = r"C:\temp\working.gdb\temp_fc" # Define and create longitude/latitude fields ## Longitude is X, Latitude is Y fieldsNew = ["GCS_WGS_1984_X", "GCS_WGS_1984_Y"] fieldsExisting = {f.name: f.type for f in arcpy.ListFields(fc)} for fNew in fieldsNew:     # Add field if it doens't already exist     if fNew not in fieldsExisting.keys():         arcpy.AddField_management(             fc,  ## in_table             fNew,  ## field_name             "DOUBLE",  ## field_type         )     # If data type of existing field is not a Double, throw Exception     elif fieldsExisting[fNew] != "Double":         raise Exception(             "Existing {} field is type {}. Must be 'Double'".format(                 fNew, fieldsExisting[fNew]             )         ) # Update fields with GCS_WGS_1984 lat/long print("Updating Field values") fieldsNew.append("SHAPE@") GCS_WGS_1984 = arcpy.SpatialReference(4326) with arcpy.da.UpdateCursor(fc, fieldsNew, "", GCS_WGS_1984) as u_cursor:     for row in u_cursor:         point = row[2].trueCentroid         row[0] = point.X         row[1] = point.Y         u_cursor.updateRow(row) python snippets​ ... View more

I created this script to help with some geodatabase house cleaning tasks. Maybe someone else will find it useful or have an idea to improve it. import arcpy import os # Set workspace myGDB = r"C:\temp\working.gdb" # Get domains that are assigned to a field domains_used = [] for dirpath, dirnames, filenames in arcpy.da.Walk(myGDB, datatype=["FeatureClass", "Table"]):     for filename in filenames:         print "Checking {}".format(os.path.join(dirpath, filename))         try:             ## Check for normal field domains             for field in arcpy.ListFields(os.path.join(dirpath, filename)):                 if field.domain:                     domains_used.append(field.domain)             ## Check for domains used in a subtype field             subtypes = arcpy.da.ListSubtypes(os.path.join(dirpath, filename))             for stcode, stdict in subtypes.iteritems():                 if stdict["SubtypeField"] != u'':                     for field, fieldvals in stdict["FieldValues"].iteritems():                         if not fieldvals[1] is None:                             domains_used.append(fieldvals[1].name)         except Exception, err:             print "Error:", err # Get domains that exist in the geodatabase domains_existing = [dom.name for dom in arcpy.da.ListDomains(myGDB)] # Find existing domains that are not assigned to a field domains_unused = set(domains_existing) ^ set(domains_used) print "{} unused domains in {}".format(len(domains_unused), myGDB) for domain in domains_unused:     arcpy.DeleteDomain_management(myGDB, domain)     print "{} deleted".format(domain) python snippets​ ... View more