Latest Contributions by XanderBakker

Below the code based on the thread mentioned by Dan Patterson‌ adapted to create squares or rectangles. Change lines 6 and 11 to point to your input and output (will be created) and change the sizes on lines 7 and 8. For some more snippets with Python have a look at: https://community.esri.com/docs/DOC-1927 import arcpy def main():     arcpy.env.overwriteOutput = True     # inputs     fc = r'D:\Xander\Rectangles\data.gdb\points'     width = 0.005 # with of the output rectangle     height = 0.005 # height of the output rectangle     # output     fcout = r'D:\Xander\Rectangles\data.gdb\rectangles'     sr = arcpy.Describe(fc).spatialReference     features = []     with arcpy.da.SearchCursor(fc, ("SHAPE@")) as cursor:         for row in cursor:             rectangle = createRectangleFromPoint(row[0].firstPoint, width, height, sr)             features.append(rectangle)     # write to output     arcpy.CopyFeatures_management(features, fcout) def createRectangleFromPoint(pnt, width, height, sr):     arrPnts = arcpy.Array()     pnt2 = arcpy.Point(pnt.X - width, pnt.Y - height)     arrPnts.add(pnt2)     pnt2 = arcpy.Point(pnt.X - width, pnt.Y + height)     arrPnts.add(pnt2)     pnt2 = arcpy.Point(pnt.X + width, pnt.Y + height)     arrPnts.add(pnt2)     pnt2 = arcpy.Point(pnt.X + width, pnt.Y - height)     arrPnts.add(pnt2)     pnt2 = arcpy.Point(pnt.X - width, pnt.Y - height)     arrPnts.add(pnt2)     return arcpy.Polygon(arrPnts, sr) if __name__ == '__main__':     main() ... View more

ArcGIS 10.2 uses 32 bits... ... View more

Great that you have advanced on the matter, but it sounds a little strange a company is not willing to invest in some tools that you will need to process 49000 LAS file... The additional time you will have to invest to do things in a less optimal way (since you lack the proper tools) is costly too. ... View more

If you plan is to run the tool 800 times, you will want to script this. Don't do something 800 manual when with little effort it can be scripted A very easy start is to perform the IDW manually and in the Results windows, right click on the IDW result and select Copy Python Snippet from the context sensitive menu. This will give you the python syntax of the IDW interpolation. Based on that you can start scripting (changing dataset, IDW settings, etc). If you want some help with it, just explain what you would like to do with the 800 interpolations (how they differ from each other). ... View more

If you don't have access to 3D analyst then your options inside ArcGIS will be limited, LAStools (as Mike Cusi suggests is a very good alternative). ... View more

There should be a difference (Inside should give multiple results, while Enter should only provide one): and... Source: ArcGIS Help (10.2, 10.2.1, and 10.2.2) Are you sure that it enters the GeoFence only once? If it is a narrow GeoFence following the route of a car, GPS imprecision may invoke multiple entries. ... View more

Try it without the double equal sign: SetNull("%wcs_int.tiff%","%wcs_int.tiff%","Value = 0") ... View more

Did ArcGIS work correctly before? Any changes occurred (other software installed)? You might want to re-install ArcGIS. ... View more

For that purpose you can use Cell Statistics to determine the range of a series of rasters (and min, max, etc): See: ArcGIS Help (10.2, 10.2.1, and 10.2.2) ... View more

If you use a ArcGIS tool from a toolbox, the messages will contain information about the elapsed time: Messages Executing: CostDistance Water SlopeSetNull C:\Forum\Pasture\gdb\Pasture.gdb\CostWater2 # # Start Time: Thu Jan 08 22:33:20 2015 Succeeded at Thu Jan 08 22:33:45 2015 (Elapsed Time: 24,92 seconds) Would that be enough? If you would use python to program a process you can simply determine the moment the process started and the moment the process ended and determine the elapsed time ( datetime.timedelta) For exact measurement in Python, you can execute code multiple times and obtains stats, using the timeit module. 26.6. timeit — Measure execution time of small code snippets — Python 2.7.9 documentation ... View more

Yesterday I performed some steps on the dataset provided by Andrew. See them listed below: # Create a hillshade for visualization arcpy.gp.HillShade_sa("ned10m34101g3","C:/Forum/Pasture/gdb/Pasture.gdb/HillShade","315","45","NO_SHADOWS","1") # Calculate the slope in percentage based on your DEM ("ned10m34101g3") arcpy.gp.Slope_sa("ned10m34101g3","C:/Forum/Pasture/gdb/Pasture.gdb/Slope","PERCENT_RISE","1") # Reclassify the slope in two classes 0-15% and >15% arcpy.gp.Reclassify_sa("Slope","Value","0 15 1;15 159,99526977539062 2","C:/Forum/Pasture/gdb/Pasture.gdb/SlopeClfy","DATA") # Set those areas to NoData that have a slope > 15% (value = 2) arcpy.gp.SetNull_sa("SlopeClfy","SlopeClfy","C:/Forum/Pasture/gdb/Pasture.gdb/SlopeSetNull","Value = 2") # Calculate the cost distance from the water points, using SlopeSetNull as cost raster (slopes > 15% are inaccessible) arcpy.gp.CostDistance_sa("Water","SlopeSetNull","C:/Forum/Pasture/gdb/Pasture.gdb/CostWater") The result looks something like this (accessibility calculated from water points are represented with blue area): ... View more

There are many option when converting Las to raster. What options did you use (Binning, Triangulation)? Is it posible to include a screenshot to show what your result is (refering to the -0.2)? ... View more