Latest Contributions by RachaelJohnson

As the title says, I am using arcpy.da.UpdateCursor() and I can't get it to update the rows.  My code is below.  (Yes, I know I am still working with a nested cursor... I haven't decided if a dictionary is the right way to go for my purposes yet.  That's for another topic.) fields = ["SOIL", "MAX_SLOPE", "MIN_CDA", "MAX_CDA", "WT_SEP", "WT_RELAX",              "COAST_SEP", "MIN_DEPTH", "DEPTH_RELA", "COAST_MIN_", "BMP", "MOD"] for row in arcpy.da.SearchCursor(combo_table, ["DARAS", "HSG", "MEDSLOPE", "MEDWT"]):      print("Saving constraint table for DA " + str(row[0]) + " and HSG " + str(row[1]) + "...")      arcpy.Copy_management(compare, BMPFold + "\\" + "DA" + str(row[0]) + "HSG" + str(row[1]))      arcpy.AddField_management(BMPFold + "\\" + "DA" + str(row[0]) + "HSG" + str(row[1]), "MOD", "TEXT", "", "", 25)      with arcpy.da.UpdateCursor(BMPFold + "\\" + "DA" + str(row[0]) + "HSG" + str(row[1]), fields) as cursor:              for r in cursor:                  # Compare water table depths                  if r[4] + r[7] > row[3]:                       if r[5] == 1:                            coast_WT = r[6]                       else:                            coast_WT = r[4]                        if r[8] == 1:                            coast_depth = r[9]                        else:                           coast_depth = r[7]                        if coast_WT + coast_depth < row[3]:                             if r[5] == 1 and r[8] == 1:                                    r[11] = "Separation and Depth"                                    print("Modifying " + str(r[11]) + "allows for " + str(r[10]) + "...")                             elif r[5] == 1:                                     r[11] = "WT Separation"                                     print("Modifying " + str(r[11]) + "allows for " + str(r[10]) + "...")                             elif r[8] == 1:                                      r[11] = "Practice Depth"                                       print("Modifying " + str(r[11]) + "allows for " + str(r[10]) + "...")                            cursor.updateRow(r)                        if coast_WT + coast_depth > row[3]:                             print("Removing " + str(r[10]) + " from DA " + str(row[0]) +                                       " HSG " + str(row[1]) +" based on water table proximity...")                             cursor.deleteRow() As the title states, I am not getting any errors when I run this code.  When I open my file, the "MOD" field is blank for all rows still.  My print statements are not printing before or after the r[11] = "TEXT" line even when they are modified to exclude calling r[11].  They worked at one point while I was getting errors with cursor.updateRow() but now there's nothing.  My last print statement works, "Removing ..." though.  ... View more

I am finding it hard to articulate what I am trying to do, so here's what I have and here's where I'm going. I have: 4 raster datasets.  One is a set of drainage zones, one is the soil type, one is slope, and one is depth below surface. I also have a large table of stormwater practices and their associated physical constraints. Desired outcome: For each drainage zone, I want a list of stormwater practices that can go on the site. What I've done:  For each drainage zone, I have found the median slope and median depth below surface using Zonal Statistics.  I have combined the median slope, median depth, drainage zone, and soil type into a single raster dataset using the Combine function.  I have used the Build Raster Attribute Table to have all of those values available to me in a single table without converting it back to a shapefile.  I have a Cursor set up to duplicate the constraint table for each unique combination of drainage zone and soil type. What I want to do:  I want to compare the Combination table to the table of constraints and delete rows from the constraint table when certain conditions are met.  For example, some stormwater practices only work on soil types A and B, so if the soil type in the constraint table is not found in the row for the combination table, I want the row in the constraint table to be deleted.  Here is the code I have come up with to help you understand my goal. # Combine rasters to give unique combinations combo = arcpy.sa.Combine([DAras, EffHSG, med_WT, med_slope]) combo.save(WorkPath + r"\combo") combo_table = arcpy.BuildRasterAttributeTable_management(combo) # Convert integer HSG code to string type arcpy.AddField_management(combo_table, "HSG", "TEXT", "", "", 6) with arcpy.da.UpdateCursor(combo_table, ["EFFECTIVEHSG", "HSG"]) as cursor:      for row in cursor:           if row[0] == 1:             row[1] = "A"         if row[0] == 2:             row[1] = "B"         if row[0] == 3:             row[1] = "C"         if row[0] == 4:             row[1] = "D"         cursor.updateRow(row) ### _____________COMPARE CRITERIA_____________________ ### print("Creating BMP output folder...") # Define location of constraint database and establish new output location CRIT = TabPath + "constraints.xlsx" BMPFold = ProjFolder + r"\BMP__Tables" if not os.path.exists(BMPFold):   os.makedirs(BMPFold) print("Loading constraint database...") # Convert Excel constraint file to GIS table compare = arcpy.ExcelToTable_conversion(CRIT, ProjFolder + r"\constraints") fields = ["SOIL", "MAX_SLOPE", "MIN_CDA", "MAX_CDA", "WT_SEP", "WT_RELAX", "COAST_SEP", "MIN_DEPTH", "DEPTH_RELA", "COAST_MIN_"] # Create output table for each drainage area for row in arcpy.da.SearchCursor(combo_table, ["DARAS", "HSG", "MEDSLOPE", "MEDWT"]):        arcpy.Copy_management(compare, BMPFold + "\\" + "DA" + str(row[0]) + "HSG" + str(row[1]))        DAmask = arcpy.sa.SetNull(DAras, DAras, "VALUE <> " + str(row[0]))        arcpy.env.mask = DAmask        arcpy.env.snapRaster = DAmask        arcpy.AddField_management(BMPFold + "\\" + "DA" + str(row[0]), "MOD", "TEXT", "", "", 14)        for r in arcpy.da.UpdateCursor(BMPFold + "\\" + "DA" + str(row[0]) + "HSG" + str(row[1]), fields):              if row[1] not in r[0]:                 del(row) The problem:  I get an error saying that "NameError: name 'row' is not defined" for line 41, "if row[1] not in r[0]:".  Once the code enters the nested cursor, "row" gets wiped out and is no longer a valid parameter.  What is a workaround for this? EDIT:  The python syntax converter works very poorly on my computer for some reason.  Sorry.  I tried to update it so the indentations are in the right place.  If they do not appear in the correct spot, I promise they are correct in my actual code and that is not the problem. ... View more

In case you're interested: I got the first part of my program to run and give good outputs! Yay! Thanks for all your help, everyone.  It took me a while to write these basic lines of code but I learned a lot in the process and I am eager to apply what I learned to the hard part of my program next.  But first... My problem: I want to script a code that will delete all the working directories created by my program. I've got this code: # This script is meant to be used to clear the working directories for a particular # map file. The user selects the map file they are trying to use the program with # and this program deletes the working folders bearing the name of the .MXD file. # This is a tool meant to clean the slate in case write errors are encountered during # tool operation. #### _______________IMPORT MODULES_________________### import os import shutil ### ________________SET ENVIRONMENTS__________________### # Set Map Document MapLoc = r"C:\Users\Rachael Johnson\Documents\GradWork\GIS\LakeRidge\LRpythontest.mxd" # Establish paths WorkPath = MapLoc[:-4] path, filename = os.path.split(MapLoc) GDBpath = MapLoc[:-4] + ".gdb" TabPath = r"C:\outputtable" "\\" ProjFolder = TabPath + filename[:-4] ### __________________DELETE PATHS ASSOCIATED WITH MXD___________ ### print("Time to delete") if os.path.exists(WorkPath):    print("deleting...")    try:   shutil.rmtree(WorkPath)    except:    raise else:    print("No such path:" + WorkPath) print("Deleted working and scratch folders in the MXD location:" + path) shutil.rmtree(ProjFolder) print("Deleted the project table output folder in" + TabPath) shutil.rmtree(GDBpath) print("Deleted the project geodatabase in the MXD location:" + path) I've also tried this code without the if statement--just shutil.rmtree(path)-- and I have the same problem.  I even tried "ignore_errors=True".  Same thing.  I have two more directories I'd like to remove but I guess I should get just one working first. I run this code and don't get any errors but the folder deletion doesn't happen either.  I added the if statement and print statements to see where the hang-up was and it doesn't even print "Time to delete", which tells me it isn't even reading the code for some reason??  I'm at a loss for what to do. EDIT: I am running Windows 8.1.  My code is in Python 2.7 or whatever version is currently used for ArcGIS 10.3.  I am running my code from PyCharm Professional v 4.0.6 ... View more

This is less of a python question and more of a general logic question, I suppose. I have a shapefile that represents the intersection of soil types and drainage boundaries.  I want to sum the total area of each soil type within each drainage boundary.  There are 4 soil types.  There are a potentially unlimited number of drainage areas. Soil type (HSG) field: "GRIDCODE" Drainage area field: "DA" Area: "Area" I have a code that extracts unique values from the drainage area field.  My problem is that I am not sure how to implement a loop sequence that iterates through each unique drainage value to find the instances of the soil types and sum the Area without manually coding each one. Earlier, I summed the areas of all the soil types with the following code: with arcpy.da.SearchCursor(SiteHSG, ["gridcode", "Area"]) as cursor:    # Initialize variables    SiteArea_A = 0.000    SiteArea_B = 0.000    SiteArea_C = 0.000    SiteArea_D = 0.000    # Calculate area for each HSG    for row in cursor:    if row[0] == 1:   SiteArea_A = SiteArea_A + row[1]    if row[0] == 2:   SiteArea_B = SiteArea_B + row[1]    if row[0] == 3:   SiteArea_C = SiteArea_C + row[1]    if row[0] == 4:   SiteArea_D = SiteArea_D + row[1] I have no idea if there was a better way to calculate those areas or not and I have no problem using a similar code over and over, basically, but I do not know how to tell the computer that for each unique DA value, find soil type, and sum those soil types, then move on and search for the next unique DA value and iterate through the grid codes and sum the area. I am kind of thinking a "while" loop may be able to accomplish what I need, but then I am not really sure how to implement a while loop in this context. I am kind of thinking of something like this so far: DAvalues = [row[0] for row in arcpy.da.SearchCursor(DAHSG, DAID)] DAunique = set(DAvalues) print(DAunique) with arcpy.da.SearchCursor(DAHSG, ["DA", "gridcode", "Area"]) as cursor: for row[0] == DAunique[0] if row[1] == 1:   SiteArea_A = SiteArea_A + row[3]    if row[1] == 2:   SiteArea_B = SiteArea_B + row[3]    if row[1] == 3:   SiteArea_C = SiteArea_C + row[3]    if row[1] == 4:   SiteArea_D = SiteArea_D + row[3] for row[0] == DAunique[1] if row[1] == 1:   SiteArea_A = SiteArea_A + row[3]    if row[1] == 2:   SiteArea_B = SiteArea_B + row[3]    if row[1] == 3:   SiteArea_C = SiteArea_C + row[3]    if row[1] == 4:   SiteArea_D = SiteArea_D + row[3] Etc.  I know there's a better way to do this than making a bunch of different for loops, especially because the number of unique values are unlimited.  It's impractical for me to code it this way but I'm not sure where to go. ... View more