Latest Contributions by LukeSturtevant

‎11-12-2015

As Darren said, if you are running your tool from ArcMap or ArcCatalog you do not need to explicitly set arcpy.env.workspace to a path. It will automatically be set to your default geodatabase as I mentioned earlier.

‎11-12-2015

Are you talking about referencing your current workspace that is set in your environment settings in ArcMap? If so then arcpy.env.workspace would correctly reference your default geodatabase.

‎11-10-2015

Kon, If you are still having this issue, maybe you would consider uploading one of your tins so that we can better understand what is happening here.

‎11-09-2015

Madan, If you have decided to go with Model Builder then you should use the Iterate Rasters function to have the tool step through the rasters in your workspace. Then string all the processing tools together as you have mentioned in your steps above. Then you could use the Collect Values function to collect all the output rasters from step 5 (subtract each cell of each raster to Value N). You can then use Cell Statistics with the 'SUM' stats type to sum up all the rasters. Let us know how you are doing!

‎11-06-2015

Are you set on using python? You could also just use model builder to string together the tools which maybe an easier way to start. The you could export the model builder tool to python script and work from there.

‎11-06-2015

The computed cellsize the tool picks is the same computation it does as if you were to pick OBSERVATION. It just does a division by 250. You can arbitrarily decrease the cellsize it picks for you because it will interpolate the cell values but you just have a better resolution raster for your viewshed. What you wouldn't want to do is increase the cellsize because then you will lose data in the raster interpolation.

‎11-06-2015

This is a very helpful code thanks Dan and Toni. Not that it is needed Toni, but if all your raster are in a single workspace and you would like to compress your code you could have something like this: import arcpy,numpy from arcpy import env env.workspace = "C:\\RasterFolder" rasters = [arcpy.RasterToNumPyArray(arcpy.Describe(ras).catalogPath,nodata_to_value=10000000000) for ras in arcpy.ListRasters()] marray = numpy.ma.masked_values(rasters,10000000000) per = 20 percentile_val = numpy.percentile(marray.compressed(),(per)) print percentile_val You could use a wild card to filter raster names and folders as well for the ListRasters function.

‎11-06-2015

Yes you should change it to CELLSIZE 10 if you want 10 meter cells. Unless your spatial reference has a linear unit of feet then that will produce a cell size of 10 feet.

‎11-06-2015

What spatial reference are you using? Did you set the sample distance to 'CELLSIZE 10' ?

‎11-05-2015

Litan, is this problem in addition to the one you posted here: i am applying interpolation tool, that time i am defining the output cell size 75, here the processing is completed but… Are you getting a raster output, but not masked to your specified extent?

‎11-05-2015

Typically this happens when your spatial reference is a geographic coordinate system the horizontal unit is an arc second or degree. If you set the output coordinate system in the environment settings to a projected coordinate system with a horizontal unit of meter it should work such as this:

‎11-03-2015

If I'm understanding your question you are looking to get a total count of layers represented in the legend? If so then you could use: mxd = arcpy.mapping.MapDocument("Current") df = mxd.activeDataFrame legendCount = len(arcpy.mapping.ListLayoutElements(mxd, "LEGEND_ELEMENT", "Legend")[0].items)

‎11-03-2015

I'm not sure if you are looking to do this with python or not, but you might look into using various List functions and arcpy.da.walk().

‎11-03-2015

Oh Okay I must have misunderstood what Siyang was trying to get at. Darren your code definitely works if he is searching for shapefiles of the same name throughout the workspace.

‎11-03-2015

I definitely agree with Darren, but as Siyang has is set up right now he will produce a dictionary with keys for every feature class name found and the only item entry will be the same feature class path. This will essential attempt to merge a single feature class for each dictionary key. If Siyang's folder structure is set up with a workspace folder with multiple folders containing many shapefiles of the same type then he could modify his code like this: import arcpy, os workspace = r'C:\Users\xxx\Desktop\boundary' output_folder = r'C:\Users\xxx\Desktop\New folder' Dict = {} for root, dirs, files in os.walk(workspace,topdown=True): for dir in dirs: env.workspace = os.path.join(root,dir) Dict[dir]=[] for fc in arcpy.ListFeatureClasses(): if not fc in Dict[dir]: Dict[dir].append(os.path.join(root,dir,fc)) else: Dict[dir].append(os.path.join(root,dir,fc)) for key in Dict: output = os.path.join(output_folder,key) + '_merge' arcpy.Merge_management(Dict[key], output) print output + " created" The code I previously posted would do the same thing, but it also allows for feature type filtering and also for a folder structure with a top level workspace with multiple folders each containing multiple subfolders with their own feature classes in them. Hope this helps!

Online Status	Offline
Date Last Visited	‎05-01-2024 04:58 PM

My Ideas

Latest Contributions by LukeSturtevant

Re: How to add a path to default database in ArcPy script?

Re: How to add a path to default database in ArcPy script?

Re: TIN to Raster

Re: multiple file, task process in python

Re: multiple file, task process in python

Re: TIN to Raster

Re: Pool of raster values to calculate percentile

Re: TIN to Raster

Re: TIN to Raster

Re: when i am applying interpolation tool, then i defined the extent(interpolation boundary) of interpolation but it did not working. so what are the parameter should be defined?

Re: i am applying interpolation tool, that time i am defining the output cell size 75, here the processing is completed but the output raster layer is not visible, so what might be the proble, plz help.

Re: How to to get current legend column count using arcpy?

Re: Results from object classes and objects!

Re: Batch merge files with same name in folders, subfolders etc.

Re: Batch merge files with same name in folders, subfolders etc.

Re: how to combine two attribute table in acgis th...

Re: Draw Line.

Re: Adding extra files to a zip?

Re: Array to polygon-Python Script

Re: ArcGIS Split Tool: How to change output name p...