Latest Contributions by DuncanHornby

Dan/Anyone else interested, So Dan showed a numpy approach which was blisteringly fast. Ouch!  But resisting the pull of numpy (because my brain does not think like that) it got me thinking about a different way of approaching this simple problem of overcoming the ever increasing slowness of the SelectLayerByLocation tool when called within a loop. My updated script which I show below burns through the same task now in 7 seconds, nice! It requires ArcPro 3.2 and above and makes use of the improved SearchCursor that can take a geometry as a spatial filter.   import arcpy, os import time arcpy.SetLogHistory(False) arcpy.SetLogMetadata(False) fp = os.path.abspath(__file__) dir = (os.path.dirname(fp)) lyrPoints = dir + r"\Test.gdb\testpoints" lyrPolygons = os.path.join(dir + r"\Test.gdb\testpoly") print(lyrPoints) # Prove to myself I actually have the full path to featureclass # Grab the polygon from layer, there is only one polygon in the layer with arcpy.da.SearchCursor(lyrPolygons, "SHAPE@") as cursor: for row in cursor: geom = row[0] print("Polygon obtained! Area is: " + str(geom.area)) d=dict() start_time = time.time() for i in range(1,3000): print("Step " + str(i)) # Use polygon as a spatial filter with arcpy.da.SearchCursor(lyrPoints, ["id"], None, None, False, None, None, spatial_filter=geom, spatial_relationship="INTERSECTS") as cursor: x = 0 for row in cursor: x += 1 # Count rows in cursor d[i] = x end_time = time.time() elapsed_time = end_time - start_time print(f" Elapsed time: {elapsed_time:.2f} seconds")       ... View more