Latest Contributions by RachaelJohnson

I've been developing my program with just one dataset and I got all the bugs worked out from that one dataset.  Now, I am trying it on some other datasets to make sure everything is working as expected. I am having a problem in the setting-up-the-environments stage of my program.  The user inputs a polygon shapefile that contains a polygon representing the parcel outline.  This polygon is used to create a mask for future raster calculations.  The basic outline of this part of the program is as follows: Buffer the outline polygon by 50 ft Convert the buffered polygon into a raster Set the raster as the mask and the snapraster for the remainder of the program An empty raster is created when converting the buffered polygon to a raster.  I tried doing it manually in GIS and had the same problem.  I added an integer field to my buffered shapefile and converted to a raster based on that integer field (rather than the ID or FID which were both 0.)  That didn't fix it.  After setting the processing extent to the same extent as my DEM, it worked.  I made those changes to my code and it still is not producing the expected results.  It's not giving me an error, but it tells me that the raster is empty and then keeps running through my code.  I've exhausted my troubleshooting abilities.  I also want to repeat: it works when I do it by hand in GIS but it isn't working via code for some reason. My code is printing out the following: Preparing necessary files... Checking for and creating output folders for GIS files... Checking for and creating output space for tables... Setting cell size to DEM cell size: 5.0 ft... Creating an environment mask from the site outline shapefile... All cells in grid k:\gradwork\gis\colleg~1\cp\scratch\rastermask have NODATA value. VAT will not be built. (VatBldNew) My code is below. Line 84 is where things are going wrong. #### _______________IMPORT MODULES_________________### print("Preparing necessary files...") import os import arcpy import copy ### ______________INPUT FILES___________________### outline = r"K:\GradWork\GIS\CollegePark\CP_outline.shp" DA = r"K:\GradWork\GIS\CollegePark\CPDrainage.shp" DAID = "Id"  # field in DA shapefile where unique DA values exist soil = r"C:\Users\Rachael J\Documents\GradWork\Project Files\BMPGIS\soils\VB_Soils.shp" WTin = r"K:\GradWork\GIS\CollegePark\CPGW_adj1" DEMin = r"K:\GradWork\GIS\CollegePark\cpelev" MapLoc = r"K:\GradWork\GIS\CollegePark\CP.mxd" WT = arcpy.Raster(WTin) DEM = arcpy.Raster(DEMin) ### ________________SET ENVIRONMENTS__________________### # Check out extension and overwrite outputs arcpy.CheckOutExtension("spatial") arcpy.env.overwriteOutput = True # Set Map Document mxd = arcpy.mapping.MapDocument(MapLoc) # Create project folder and set workspace print("Checking for and creating output folders for GIS files...") WorkPath = MapLoc[:-4] if not os.path.exists(WorkPath):     os.makedirs(WorkPath) arcpy.env.workspace = WorkPath # Create scratch workspace ScratchPath = str(WorkPath) + r"\scratch" if not os.path.exists(ScratchPath):     os.makedirs(ScratchPath) arcpy.env.scratchWorkspace = ScratchPath # Create GDB path, filename = os.path.split(MapLoc) GDB = filename[:-4] + ".gdb" GDBpath = MapLoc[:-4] + ".gdb" if not os.path.exists(GDBpath):     arcpy.CreateFileGDB_management(path, GDB) # Create output table folder if it does not exist and create project folder print("Checking for and creating output space for tables...") TabPath = r"C:\outputtable" "\\" ProjFolder = TabPath + filename[:-4] if not os.path.exists(TabPath):     os.makedirs(TabPath) if not os.path.exists(ProjFolder):     os.makedirs(ProjFolder) ### _____________SET PROCESSING EXTENTS____________ ### # Set cell size description = arcpy.Describe(DEM) cellsize = description.children[0].meanCellHeight print("Setting cell size to DEM cell size: " + str(cellsize) + " ft...") arcpy.env.cellSize = cellsize arcpy.env.extent = arcpy.Describe(DEM).extent # Create buffer around outline to use as mask # Buffer distance is in feet print("Creating an environment mask from the site outline shapefile...") maskshp = arcpy.Buffer_analysis(outline, ScratchPath + r"\outline_buff", "50 Feet", "", "", "ALL",) arcpy.AddField_management(maskshp, "RID", "SHORT") with arcpy.da.UpdateCursor(maskshp, ["RID"]) as cursor:     for row in cursor:         row[0] = 8         cursor.updateRow(row) # Convert buffer to raster mask = arcpy.Raster(arcpy.PolygonToRaster_conversion(maskshp, "RID", ScratchPath + r"\rastermask")) mask.save(ScratchPath + r"\rastermask") # Set raster mask and snap raster print("Setting raster mask and snap raster for project...") arcpy.env.mask = mask arcpy.env.snapRaster = mask EDIT:  Changing line 69 to "arcpy.env.extent = DEM.extent" did not solve the problem. EDIT 2:  My files aren't in the same coordinate system! ... View more

I am comparing the values in one table to the values in another table.  I converted the values of one table into a dictionary to avoid using nested cursors.  Ultimately, I want to create a table of items that pass each comparison test.  I did another version of this where I have separate tables for each area I am testing, start with a full table, and remove rows from the table as each row fails a test.  This results in many tables, which is fine for a handful of areas but messy for more than a few.  I am trying to build a version that generates all area outputs in the same table. I think it is easier to keep track of those that fail tests than those that pass.  I copy the main dictionary, D, to a new dictionary.  I create a blank dictionary for the failed tests.  Each time something fails a test, I record its key into the Failed Test dictionary and store the reason for failure as a string value.  For subsequent tests, if the key already exists in the Failed test dictionary, I add on to the existing value as a pseudo-list to prevent test failures from being overwritten.  When all tests are complete, I want to compare the keys in the Failed Test dictionary to the copy of D.  If the keys are a match, I want that key and it's associated value removed from the copy of D.  What should remain in the copy of D are the BMPs that pass all 4 tests.  Then I want to use the values stored in the copy of D to write to a table.  Maybe this logic is all sorts of messed up and there's a better way to do it.  I am open to suggestions.  I have been playing with this code for a while to figure out what the problem is and it seems I have found the first source of error in my code.  There appears to be a problem with copying the dictionary, D.  Here's the code I have so far.  It is missing the write-to-table part because I am not getting the outputs I would expect. print("Comparing site values to constraints...") # Compare the lumped parameters to the constraint dictionary for row in arcpy.da.SearchCursor(combo_table, ["DARAS", "HSG", "MEDSLOPE", "MEDWT"]):     BMP = {}     # Temporarily store the area of each DA for later comparison     for r in arcpy.da.SearchCursor(DA_area, [DAID, "SUM_AREA"]):     # NOTE: THE DAID/DARAS MISMATCH COULD BE A PROBLEM LATER. NEED TO TEST.         if r[0] == row[0]:             totalArea = r[1]     # Duplicate criteria dictionary that can be amended throughout the loop     print D     BMP = D     print D     # Initialize empty dictionary to store BMPs that fail each test     NoBMP = {}     # Compare lumped values in each DA/HSG pair to those in the constraint table     print(row[0], row[1])     for k, v in D.items():         # Test if soil type is incorrect for each BMP         if row[1] not in v:             # Add failed BMP to the failed BMP dictionary, NoBMP             NoBMP = "Soil Type Mismatch"         # Compare keys in BMP and NoBMP dictionaries. Remove matching pairs from the BMP dictionary.     for key in BMP.keys():         if key in NoBMP.keys():             del BMP[key] The main takeaway from the following output information is that D changes sizes from one iteration of the Cursor to the next iteration of the Cursor; D does not change size in future iterations of the loop.  I think it has something to do with me changing the size of BMP at the end of the code by deleting the matching keys.  The keys that are removed are the keys that do not show up in the next iteration.  However, I reinitialize BMP as a blank dictionary at the beginning of the Cursor, so I am not sure what else to do.   I have also tried BMP = dict() to reinitialize and gotten the same result.  Comparing site values to constraints... {u'WPGW1': u'C,D', u'WPGW2': u'C,D', u'VFA': u'A', u'SDCD': u'C,D', u'BRE2': u'A,B,C,D', u'BRE1': u'A,B,C,D', u'SDAB': u'A,B', u'EDP2': u'A,B,C,D', u'WP1': u'C,D', u'CW1': u'C,D', u'CW2': u'C,D', u'DS2': u'A,B,C,D', u'DS1': u'A,B,C,D', u'GCCD': u'C,D', u'F2': u'A,B,C,D', u'F1': u'A,B,C,D', u'CAAB': u'A,B', u'WP2': u'C,D', u'GCAB': u'A,B', u'CACD': u'C,D', u'VFSA': u'B,C,D', u'SI2': u'A,B', u'SI1': u'A,B', u'EDP1': u'A,B,C,D', u'SDSA': u'C,D', u'WS1': u'C,D', u'CI2': u'A,B', u'CI1': u'A,B', u'WS2': u'C,D', u'GCSA': u'C,D', u'MI1': u'A,B', u'MI2': u'A,B'} {u'WPGW1': u'C,D', u'WPGW2': u'C,D', u'VFA': u'A', u'SDCD': u'C,D', u'BRE2': u'A,B,C,D', u'BRE1': u'A,B,C,D', u'SDAB': u'A,B', u'EDP2': u'A,B,C,D', u'WP1': u'C,D', u'CW1': u'C,D', u'CW2': u'C,D', u'DS2': u'A,B,C,D', u'DS1': u'A,B,C,D', u'GCCD': u'C,D', u'F2': u'A,B,C,D', u'F1': u'A,B,C,D', u'CAAB': u'A,B', u'WP2': u'C,D', u'GCAB': u'A,B', u'CACD': u'C,D', u'VFSA': u'B,C,D', u'SI2': u'A,B', u'SI1': u'A,B', u'EDP1': u'A,B,C,D', u'SDSA': u'C,D', u'WS1': u'C,D', u'CI2': u'A,B', u'CI1': u'A,B', u'WS2': u'C,D', u'GCSA': u'C,D', u'MI1': u'A,B', u'MI2': u'A,B'} (1, u'D') {u'WPGW1': u'C,D', u'WPGW2': u'C,D', u'SDCD': u'C,D', u'BRE2': u'A,B,C,D', u'BRE1': u'A,B,C,D', u'EDP2': u'A,B,C,D', u'WP1': u'C,D', u'CW1': u'C,D', u'CW2': u'C,D', u'DS2': u'A,B,C,D', u'DS1': u'A,B,C,D', u'GCCD': u'C,D', u'F2': u'A,B,C,D', u'F1': u'A,B,C,D', u'WP2': u'C,D', u'CACD': u'C,D', u'VFSA': u'B,C,D', u'EDP1': u'A,B,C,D', u'SDSA': u'C,D', u'WS1': u'C,D', u'WS2': u'C,D', u'GCSA': u'C,D'} {u'WPGW1': u'C,D', u'WPGW2': u'C,D', u'SDCD': u'C,D', u'BRE2': u'A,B,C,D', u'BRE1': u'A,B,C,D', u'EDP2': u'A,B,C,D', u'WP1': u'C,D', u'CW1': u'C,D', u'CW2': u'C,D', u'DS2': u'A,B,C,D', u'DS1': u'A,B,C,D', u'GCCD': u'C,D', u'F2': u'A,B,C,D', u'F1': u'A,B,C,D', u'WP2': u'C,D', u'CACD': u'C,D', u'VFSA': u'B,C,D', u'EDP1': u'A,B,C,D', u'SDSA': u'C,D', u'WS1': u'C,D', u'WS2': u'C,D', u'GCSA': u'C,D'} (2, u'D') ... View more