Latest Contributions by BrookeHodge

Hi all, I'm hoping someone can help me a problem I can't figure out.  Below is my code.  I keep getting a RuntimeError: cannot open 'file I'm writing to' using a arcpy.da.InsertCursor.  It keeps locking for some reason, I do not have arc open when running (running it from PyCharm), and it will write to the file thousands of times... I am writing line segments to a feature class and it will do it thousands of times, but then seemingly randomly say it can't access it.  I have run it many times and it get's locked up in different places each time.  The only info I can find is that people need to close their cursors, but all my cursors (both search and insert cursors) are closed after each use.  I'm not sure what else to do.  Any ideas?  Description of what code is trying to do below if that is helpful at all as well as a screen shot of the error. import arcpy # Set environments / workspace. Workspace should point to a gdb that has all the effort point feature classes arcpy.env.workspace = r"C:\Users\bhodge\Dropbox (New England Aquarium)\MonumentsWork\MonumentSpeciesDiversity\DataFromBob\DerivedData.gdb" arcpy.env.overwriteOutput = True # This is the input points fc = r"C:\Users\bhodge\Dropbox (New England Aquarium)\MonumentsWork\MonumentSpeciesDiversity\DataFromBob\DerivedData.gdb\EffortPoints_Prj_TEST_B" # This is the output lines sr = arcpy.Describe(fc).spatialReference outputFC = arcpy.management.CreateFeatureclass(r"C:\Users\bhodge\Dropbox (New England Aquarium)\MonumentsWork\MonumentSpeciesDiversity\DataFromBob\DerivedData.gdb", "EffortLines_TEST", "POLYLINE", "", "DISABLED", "DISABLED", sr) # Add desired fields to outputFC arcpy.AddField_management(outputFC, "MONTH", "TEXT") arcpy.AddField_management(outputFC, "YEAR", "SHORT") arcpy.AddField_management(outputFC, "FILEID", "TEXT") arcpy.AddField_management(outputFC, "EVENTNO", "LONG") # Make a feature Layer of the input points to work off. This layer is temporary and will not persist after session ends unless saved arcpy.MakeFeatureLayer_management(fc, "fc_lyr") # Sort the table by fileid and eventno and create a new fc called fc_sort sort_fields = [["FILEID", "ASCENDING"], ["EVENTNO", "ASCENDING"]] fc_Sorted = arcpy.Sort_management("fc_lyr", "fc_sort", sort_fields) # Create the search cursor fields = ['FILEID', 'beauf_num', 'LEGTYPE', 'LEGSTAGE', 'vis_num', 'olviz_num', 'alt_num', 'LAND', 'MONTH', 'YEAR', "EVENTNO", "SHAPE@"] with arcpy.da.SearchCursor(fc_Sorted, fields) as cur1: for row in cur1: # Name variables and assign values starting on first record of table fileid1 = row[0] beaufort = row[1] legtype = row[2] legstage = row[3] visiblty = row[4] oldvis = row[5] alt = row[6] land = row[7] month = row[8] year = row[9] eventno = row[10] # Get the geometry of the point for that record geom1 = row[11].getPart() # Go to the next row, get the fileid and geometry next_row = cur1.next() fileid2 = next_row[0] land2 = next_row[7] geom2 = next_row[11].getPart() # Create an Array array = arcpy.Array() # Test the parameters to see if it should be on-effort if ((fileid1 == fileid2) and (beaufort != None and beaufort <= 4) and ((visiblty != None and visiblty >= 2) or (oldvis == 1)) and \ (land != 1 and land2 != 1) and (alt != None and alt <= 366) and \ (not legtype == '0') and \ (not (legtype == '7' and legstage == ' ')) and \ (not (legtype == '9' and legstage == ' ')) and \ (not (legtype == '5' and legstage == ' '))): array.add(geom1) array.add(geom2) polyline = arcpy.Polyline(array, sr) with arcpy.da.InsertCursor(outputFC, ['SHAPE@', "MONTH", "YEAR", "FILEID", "EVENTNO"]) as insertCursor1: insertCursor1.insertRow([polyline, month, year, fileid1, eventno]) del insertCursor1 array.removeAll() else: pass del cur1 print("first with loop completed") with arcpy.da.SearchCursor(fc_Sorted, fields) as cur2: cur2.next() for row in cur2: # Name variables and assign values starting on first record of table fileid1 = row[0] beaufort = row[1] legtype = row[2] legstage = row[3] visiblty = row[4] oldvis = row[5] alt = row[6] land = row[7] month = row[8] year = row[9] eventno = row[10] # Get the geometry of the point for that record geom1 = row[11].getPart() # Go to the next row, get the fileid and geometry next_row = cur2.next() fileid2 = next_row[0] land2 = next_row[7] geom2 = next_row[11].getPart() # Create an Array array = arcpy.Array() # Test the parameters to see if it should be on-effort if ((fileid1 == fileid2) and (beaufort != None and beaufort <= 4) and ((visiblty != None and visiblty >= 2) or (oldvis == 1)) and \ (land != 1 and land2 != 1) and (alt != None and alt <= 366) and \ (not legtype == '0') and \ (not (legtype == '7' and legstage == ' ')) and \ (not (legtype == '9' and legstage == ' ')) and \ (not (legtype == '5' and legstage == ' '))): array.add(geom1) array.add(geom2) polyline = arcpy.Polyline(array, sr) with arcpy.da.InsertCursor(outputFC, ['SHAPE@', "MONTH", "YEAR", "FILEID", "EVENTNO"]) as insertCursor2: insertCursor2.insertRow([polyline, month, year, fileid1, eventno]) del insertCursor2 array.removeAll() else: pass del cur2 print("Done")‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ This is the Error I get, but this changes depending on the run.  Sometimes it's in in insertCursor1 sometimes in insertCursor2 depending on how far it get's that run. "C:\Program Files\ArcGIS\Pro\bin\Python\envs\arcgispro-py3\python.exe" "C:/Users/bhodge/Dropbox (New England Aquarium)/MonumentsWork/MonumentSpeciesDiversity/MonumentsPython/Trackline_Fields.py" Traceback (most recent call last): File "C:/Users/bhodge/Dropbox (New England Aquarium)/MonumentsWork/MonumentSpeciesDiversity/MonumentsPython/Trackline_Fields.py", line 59, in <module> with arcpy.da.InsertCursor(outputFC, ['SHAPE@', "MONTH", "YEAR", "FILEID", "EVENTNO"]) as insertCursor1: RuntimeError: cannot open 'C:\Users\bhodge\Dropbox (New England Aquarium)\MonumentsWork\MonumentSpeciesDiversity\DataFromBob\DerivedData.gdb\EffortLines_TEST' Process finished with exit code 1‍‍‍‍‍‍‍ Code description: Basically I'm running through a point table, looking at variables in the first row, and some in the next row and if those parameters meet my requirements, I draw a line to connect the points and save that to a featureclass I created along with some variables I need attached to those line segments.  I run through two searchCursors because of how I need to access the first and second line, it skips one iteration, (basically the first one compares points 1&2, 3&4, etc.) so I run it through a second searchCursor accessing the opposites (comparing 2&3, 4&5, etc.).  It's somehow locking the file I'm trying to write to in the insertCursor, but it does so randomly and I don't know why or how to fix it. Thanks in advance for any suggestions! ... View more

I'm going crazy and hoping someone can help.  I'm trying to bring in a point feature class, go row by row (using a .da.searchcursor), add two points into an array, and create a polyline from them. Then use an da.update cursor to add that line segment into a line feature class.  Then go to the next row and do the same until I have made a number of line segments from the point features.  I have tried this a ton of different ways, all with the same result, an updated line feature class that has the added rows, but NO lines draw and the length field for all segments are 0.  I've altered this to just run the points and add the points to a new point feature and that works just fine.  This code has also been updated to incorporate the .da.cursors (I have similar code running off the old cursors that work, i'm just trying to update it and the .polyline is not working right).  Does anyone have any ideas or things to try?  I'm going crazy. TIA! import arcpy # Set environments / workspace. Workspace will be a parameter in tool. This parameter should point to a gdb that has all the effort point feature classes arcpy.env.workspace = r"C:\Users\bhodge\Dropbox (New England Aquarium)\MonumentsWork\MonumentSpeciesDiversity\DataProcessing\MichelleCData.gdb" arcpy.env.overwriteOutput = True # input point fc fc = r"C:\Users\bhodge\Dropbox (New England Aquarium)\MonumentsWork\MonumentSpeciesDiversity\DataProcessing\MichelleCData.gdb\SurveyPointsUTM19_TestSet_Short_6" arcpy.MakeFeatureLayer_management(fc, "fc_lyr") # Sort the table by fileid and eventno and create a new fc called fc_sort sort_fields = [["FILEID", "ASCENDING"], ["EVENTNO", "ASCENDING"]] fc_Sorted = arcpy.Sort_management("fc_lyr", "fc_sort", sort_fields) # Create the search cursor fields = ["SHAPE@"] with arcpy.da.SearchCursor(fc_Sorted, fields) as cur: for row in cur: geom1 = row[0].getPart() cur.next() geom2 = row[0].getPart() array = arcpy.Array() array.add(geom1) array.add(geom2) polyline = arcpy.Polyline(array, arcpy.SpatialReference(26919)) with arcpy.da.InsertCursor(r"C:\Users\bhodge\Dropbox (New England Aquarium)\MonumentsWork\MonumentSpeciesDiversity\DataProcessing\MichelleCData.gdb\OnEffortLinesTEST", 'SHAPE@') as insertCursor: insertCursor.insertRow([polyline]) array.removeAll() del insertCursor del cur The result is all the points used draw and are correct (the fc_Sorted /fc_sort), and the line feature class seems to updated in that rows are added, but no lines draw and the Length field is 0. ... View more