Latest Contributions by DarrenWiens2

The way it is now, you've listed some feature classes in the variable 'fc', listed some feature classes in your workspace into the variable 'vertices', destroyed the 'fc' list, and are looping through everything in 'vertices' in a new variable called 'fc'. I think you may want: import arcpy, os #array of files to be processed fcs = ['D:\\Broomfield.Lidar\\HARN_Vertices.gdb\\VertexN2W187North', 'D:\\Broomfield.Lidar\\HARN_Vertices.gdb\\VertexN2W188', 'D:\\Broomfield.Lidar\\HARN_Vertices.gdb\\VertexN3E180', 'D:\\Broomfield.Lidar\\HARN_Vertices.gdb\\VertexN3E190', 'D:\\Broomfield.Lidar\\HARN_Vertices.gdb\\VertexN3E191'] outputFolder = r'C:\junk' # or whatever you choose for fc in fcs:      fcName = os.path.basename(fc) # returns 'VertexN2W187North', etc.     arcpy.MinimumBoundingGeometry_management(fcName, os.path.join(outputFolder, fcName + "_Poly"), "RECTANGLE_BY_AREA", "ALL", "", "NO_MBG_FIELDS") ... View more

Just curious, are you using 32- or 64-bit Python? I can't seem to shake these Memory Errors with 32-bit. Any more than 10,000 x 10,000 comparisons returns an error, but I may be sloppy managing memory. edit: or when you say you worked around it, do you mean you processed chunks at a time? ... View more

Oh, you may want to check Tabulate Area. With both of these, if you have overlapping polygons, you will have to process zone features individually, since the zones are converted to raster (flattened), internally. ... View more

Zonal Statistics as Table should do this for you. ... View more

This thread is getting too long to read, so sorry if this has been mentioned and/or ruled out. You may be able to save yourself some time by Merging (that is, edit merge) or Dissolving all lines into one massive multipart polygon, then Exploding or running Multipart to Singlepart to separate into touching clusters. Finally, Spatial Join named segments to clusters. You will have to manually inspect clusters that should have more than one name. ... View more

Sorry to say, but I can't reproduce your error. I'm using version 6.0 for 10.2. I don't see anything wrong with your syntax - the backslash situation in the call to ETS_GPIDW raised an eyebrow, but I could also save to \Sal.img and \pSal.img, so I don't think that's it. This works for me: import arcpy arcpy.ImportToolbox(r'C:\Program Files (x86)\ET SpatialTechniques\ET Surface 6.0 for ArcGIS 10.2\ET Surface.tbx') arcpy.ETS_GPIDW("C:\junk\csv_points.shp", "C:\junk\idw.img", "Elev", 25, 2, 12, 2000) ... View more

If you're feeling brave, install PySAL and use the natural breaks function once you have a list of values (get using a SearchCursor​, or better yet FeatureClassToNumPyArray). I'd consider this an advanced manoeuvre which would probably spawn many more questions, but worth noting. ... View more

Probably not possible in batch processing. Possibly possible in model builder using an iterator. Definitely possible using python. ... View more

Hmmm, well my numpyNear() function gives 'MemoryError' at around 10,000 points. Whether that's due to SciPy or my setup, I don't know. edit: I see you said comparisons, so I guess that's about in line with 10,000 points compared to itself: 100,000,000 distance comparisons. ... View more

No, you're right, it isn't apples-to-apples (I meant to say so). stat_points() is much more fully functional, but that comes at a cost. My main observation was that your pure numpy method is likely done behind the scenes using scipy, which comes at its own cost. ... View more

Just ran a quick test between stat_point() and the following function for 1,000 points and the results are: stat_point() - 29s, numpyNear() - 0.5s. Almost of the stat_point() time is due to calls to arcpy.Describe (14.1s) and arcpy.AddField (14.7s), so the numpy stuff isn't the problem. ... def numpyNear(fc): ...     npArray = arcpy.da.FeatureClassToNumPyArray(fc,["OID@","SHAPE@XY"]) ...     distArray = scipy.spatial.distance.pdist(npArray['SHAPE@XY'],'euclidean') ...     sqDistArray = scipy.spatial.distance.squareform(distArray) ...     nearFID = npArray['OID@'][numpy.argsort(sqDistArray)[:,1]].transpose() ...     npAppend = numpy.lib.recfunctions.append_fields(npArray,'NearFID',nearFID) ...     nearDist = numpy.sort(sqDistArray)[:,1].transpose() ...     npAppend = numpy.lib.recfunctions.append_fields(npAppend,'NearDist',nearDist) ...     outFC = r'in_memory\outPts'+str(random.random()*1000) ...     outFeat = arcpy.da.NumPyArrayToFeatureClass(npAppend,outFC,"SHAPE@XY") ...     arcpy.CopyFeatures_management(outFC,r'in_memory\numpy_Pts') ... View more

I should probably read this more closely, but I've used scipy.spatial.distance.pdist for approximately what I think may be happening in lines 43-45 (pairwise distances?). ... View more

See this help page for how to use the Check Values button to generate output file names, based on inputs: ... View more