Latest Contributions by LukeWebb

That Error Message only tells us the line it occurs, you seem to have left of the actual Error from the message.  ... View more

Joshuas answer seems the best way.   Another take, is to use the tool you were doing, but only run it on 1 feature at a time, something like: src_shp = r'a_shapfile' result_fc = 'a_different_shapefile' for row in arcpy.SearchCursor(src_shp): row_oid = row.getValue('OBJECTID') row_where = """ObjectID = %s""" % row_oid active_lyr = arcpy.MakeFeatureLayer(src_shp, where_clause=row_where) if arcpy.Exists(row_result): arcpy.Delete(row_result) row_result = Feature to Polygon (active_lyr, temp_layer) append_management(row_result, result_fc) ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ I imagine Joshuas approach above will perform much faster though! ... View more

Ive always done.. Add a Unique ID Field to data "Generate Points along lines"   (Choose a resolution for the vertices, that is relevant to your data without going crazy!) "Feature Vertices To Points"    "Create Thiessan Polygons" "Dissolve"  using the UniqueID Field Works great..ish for me.   (The boundarys can get a bit zigzaggy!) ... View more

Im trying to understand what your achieving with this code: with arcpy.da.SearchCursor(feature_class, 'SHAPE@AREA') as acursor:     for row in acursor:         if row[0] > 100000000:             print("{} greater than 10K km2, splitting...".format(row[0])             tgrid = arcpy.GenerateTessellation_management(t_out, oe, "SQUARE", 100000000, sr)             print("Splitting complete")         else:             print("point extent smaller than 10K km2, not splitting.") To me, you need to do something with TGrid, otherwise you are just overwriting each iterations previous output. I think this is also what causes your issue as you are still "Referencing" the output in the python script, so it will lock the shapefile maybe. You could test this theory, by changing it to: with arcpy.da.SearchCursor(feature_class, 'SHAPE@AREA') as acursor:     for row in acursor:         if row[0] > 100000000:             print("{} greater than 10K km2, splitting...".format(row[0])             tgrid = arcpy.GenerateTessellation_management(t_out, oe, "SQUARE", 100000000, sr)             print("Splitting complete")             del(tgrid)         else:             print("point extent smaller than 10K km2, not splitting.") ... View more

Theres a slight chance if you had the layer in the map the whole time, that ArcGIS had cached the Shape column so its not loading all the data from disk permanantly. You could try dragging the layer into a new ArcMap window to be sure!   Or add a call at the end of your script to http://desktop.arcgis.com/en/arcmap/10.3/analyze/arcpy-functions/refreshactiveview.htm For me, normally if i zoom in and out, then it seems to reload the data as well. (Panning not so much) ... View more

when it says "City_Lyr is not defined", what this means is, you have not yet "defined" your city layer This is because in this line of code: for lyr in arcpy.mapping.ListLayers(mxd, "City*", df):     if lyr.name == "City Points (2010)":         city_lyr = lyr it never actually gets to run the line of code:  city_lyr = lyr Which is where you try to define it. The reason for that, I am not so sure, but the "If" statement is somehow not letting any of your layers get defined as the city lyr. Same issue regarding your other layer. ... View more

Dans Answer is correct. in_fc = 'C:/Your_spaceless_path/CoordGeom.gdb/Polygons' arr = arcpy.da.TableToNumPyArray(in_fc, 'Shape_Area')  # pick a field # nan functions account for nulls, cast to the dtype of the field using astype np.nanmax(arr.astype('float'))  155.0 If this is returning different numbers based on selection, then this is because you are not entering "in_fc" as shwon above. E.g. in_fc = 'C:/Your_spaceless_path/CoordGeom.gdb/Polygons' Is different than: in_fc = 'Polygons' One points to a layer in a map, and will return results based on selections and things, the other points to data on a disk, and will return counts of whatever is inside that dataset. The only possible thing that could go wrong using this, is if a user is inside an Edit session, then there features are not saved to the "dataset on disk" until they save Edits. ... View more

Thanks great to see the shapefile thing has changed! My ArcPro is a bit out of date, as we use a concurrent server now instead of named users, and im not able to download the update since the switch. Need to contact my Org admin to get help with that! With regards to new folders, this is specifically when navigating folders, via any geoprocessing "Select output path" interface.  In a perfect world I would create the folder before opening a tool, but I tend to forget as I always opened the tool, then created the folder for the past 15 years ... View more

Thanks Rory, I think I must have been thinking you cant Right click inside a folder to create shapefiles. (I know you can with FGDBs and FCs, but we still use shapefiles for most workflows) This is why I believe its an ESRI design choice to force us to adopt more modern and improved practises, but our ecosystem exists beyond esri and other softwares do not tend to interface natively with these FGDBs. Another example -  When choosing an output location of any geoprocess tool. In ArcMap whilst browsing locations I could simply click -- New folder, or new GDB, in ArcPro this isnt there. So I either have to close the window then go and make the GDB in catalog then reopen it, or use windows to go and make a folder.   Simple usability things that ArcGIS has these covered and works for everything I need stop me from considering upgrading. The new good things are very good and the one time I need to use them I load ArcPro, but it fails for me on the activities I typically do 100s of times a day.  (Creating shapefiles, and choosing output locations) ... View more

Here is a function I use to build where cluases from lists, I use it like this. values = ['test', 'test2'] field_name = 'query_field' sql_where = build_where_clause_from_list (field_name, values) values = [11, 27] field_name = 'query_field' sql_where = build_where_clause_from_list (field_name, values, field_type='Number') def build_where_clause_from_list(field_name, value_list, field_type="String"): """Build a SQL In Where clause using a list e.g. field_name = 'Cheese', value_list=['Edam', 'Cheddar', 'Stinky Bishop'] output: "Cheese" in ('Edam', 'Cheddar', 'Stinky Bishop') """ where = """ "%s" in (""" % field_name if field_type == 'String': where+= "'" if field_type == 'String': joinTxt = "','" else: joinTxt = "," where += joinTxt.join(map(str, value_list)) if field_type == 'String': where += "'" where += ")" return where ... View more

The short and dirty answer is, inside Pycharm / Spyder / any python IDE, there is a setting "Project Interpreter", which is pointing to a python.exe in a folder such as C:\python27\python.exe You need to update this setting to point to a place such as: C:\Program Files\ArcGIS\Pro\bin\Python\python.exe ... View more

Ive found that if I store the code in a temporary Microsoft Format. (Word, Excel etc) The single quotes, are not compatible with some of the ESRI stuff.   But if you copy the same text into notepad, then copy that into Arc, many tools will work. ... View more