0
votes
0answers
54 views

What is GDAL equivalent to ArcGIS Lookup (Spatial Analyst) tool?

I use the Lookup (Spatial Analyst) tool to create several raster layers from a combined raster layer. I know ArcPy can do the same thing. However, I would like to know whether GDAL (or another ...
1
vote
0answers
21 views

ExtractByMask with one mask and multiple layers ArcMap

I have a multiple cost distance rasters and multiple feature layers. I need to find and join the feature layers (or associated shapefiles) that lie within a threshold of costs on the cost rasters. For ...
1
vote
2answers
70 views

Open source, Python option for Arcpy's Con tool?

I'm wondering if there is a method to perform a function similar to Arc's Con tool, using open source options in python. I have two rasters, one that shows locations of clouds and the other is just a ...
0
votes
1answer
92 views

AttributeError: 'int' object has no attribute 'save'

I am trying to use a conditional statement to generate a raster with binary values from a raster with probability values (floating point raster). This is easily done within ArcMap by batching the ...
0
votes
0answers
23 views

Using GDAL and Python to replace the Spatial Analyst TabulateArea method

Has anyone tried to phase out the use of ESRI's Spatial Analyst extension by replacing it with GDAL in Python? I have a script that only uses the TabulateArea and Statistics analysis but not much else ...
1
vote
1answer
40 views

Script to “extract by mask” a list of raster files

I'm trying to mask a group of raster files and save the new rasters in a different folder using the same name of the original files. It doesn't sound as a difficult task, but I'm having trouble ...
0
votes
0answers
24 views

slope in model different value than slope in stand alone script [duplicate]

When running slope(analysis) from model I receive a high value of somewhere around 71. After then exporting model as python snippet and changing only the output raster name (using same input raster) I ...
0
votes
1answer
66 views

Weighted Overlay - passing parameter as text error?

For some reason when hardcoding the save path for the wieghted raster just before runtime it works perfectly. Like such... arcpy.env.cellSize = "90" arcpy.env.workspace = r'G:\Xcel\Route ...
2
votes
1answer
221 views

Multiprocessing issues with ArcPy

I am using ArcGIS 'watershed' routine in a script which makes use of multiprocessing. This works fine, but during execution I get the message: Unable to remove directory. Possible causes: 1- Not ...
3
votes
1answer
64 views

Any way to perform zonal statistics on column other than value or do reclass faster?

I have to perform ArcGIS zonal statistics on a raster, but the column I am interested in is not the 'value' column. Now, I can reclass the raster so that the column of interest becomes the 'value' ...
1
vote
2answers
49 views

Fill Sinks of multiple rasters in a Python expression

I have 31 rasters that represent different watersheds and I would like to use the Fill (Spatial Analysis) tool to fill the sinks of all rasters at once in Python. All the rasters have the same name ...
1
vote
0answers
180 views

Watershed delineation; how to program

I've been trying to roll my own watershed delineation program, taking a flow direction grid and pour points as input. I coded up a quick search algorithm in Python which inspects the flow direction ...
1
vote
1answer
41 views

How to build and pass remap for floating point raster without getting Reclassify ERROR 010255?

I am receiving the following message when executing my python script: ERROR 010255: String remap must be used with a string field. My input raster is floating point and I would like to ...
0
votes
0answers
77 views

Loop - convert layer to raster coming up blank

I am looping through files, selecting out building footprints one by one and converting the selected building footprint to a raster. The script worked fine for about 2 hours this morning then stopped ...
3
votes
1answer
228 views

How to exclude stations from IDW interpolation?

I have the following problem: I need to calculate an average precipitation amount on a catchment using an IDW interpolation. The interpolation is performed in a loop since each precipitation station ...
1
vote
0answers
33 views

Problem in analyzing data in a stream segment

I have to find out the head difference along the stream segment. For that I used Spatial Analyst Tool>Extraction>Sample. With that I got the elevation values in points. But after doing so, i found ...
2
votes
2answers
226 views

Tabulate Area crashes Python and ArcGIS

I have a script that runs Tabulate Area on a single feature class multiple times with different rasters. In total Tabulate Area runs 5 times. 2 of the input rasters are provided as arguments to the ...
0
votes
0answers
48 views

How specify a corridor minwidth?

I’m trying to make a corridor between two substations with ArcGIS, a macro corridor and a micro corridor. I have two questions in this respect and would be grateful for your suggestions. a) How can ...
1
vote
0answers
159 views

Problems running zonal statistics - not all zones represented

I have been working on a Python script to automate the steps illustrated by whuber here:Create points around the edge of viewshed results I am trying to create a horizon graph illustrating the ...
3
votes
1answer
377 views

Why is my simple arcpy script failing? (sets raster cell value to null if below threshold value, ArcGIS 10.1)

As mentioned in the title, simple script to set null cell values below three standard deviations beyond the mean: import arcpy, sys, os #Accept input and output raster inRas = sys.argv[1] outRas = ...
5
votes
2answers
199 views

Is there a way to increase the search radius of nearest neighbors for the Nibble tool?

I am working with a raster map vegetation types. I have about 20 types with unique raster code values that can be sorted into three classes: 1 - vegetation types I want to expand using the nibble ...
3
votes
2answers
548 views

Majority filter (minimum mapping unit) in ArcGIS with larger window size

I want to run the majority filter in ArcGIS with a window size of 4 or 5 i.e. greater than the default 3. Is it possible? Or some other command (or tool) which might help me achieve that? thanks!
5
votes
1answer
496 views

Work-around for slow processing times in ArcGIS 9.3 with Spatial Analyst functions

I am batch processing a large amount of distance functions (~50,000) using ArcGIS 9.3. In this version, using spatial analyst functions, it seems to be the case that users cannot parallelize processes ...
2
votes
4answers
326 views

Error When Trying to Run Slope Analysis Script

I have written a script which is supposed to loop through all the rasters in a folder and perform a slope analysis on each raster and save the results in another folder. I keep getting the following ...
1
vote
1answer
85 views

sa.fill creates NoData raster in standalone script ArcMap 10.0

I am trying to put together a geoprocessing workflow in Python, and I've hit a stumbling block with sa.fill. When run as a standalone script, the tool produces a raster composed entirely of NoData ...
1
vote
2answers
90 views

Distribute gridded data to corresponding subbasin boundaries

I have one ArcGIS shapefile of my river basin consisting of many subbasin. Another shapefile showing the location of gridded points over the whole basin (Its latitude and longitude are used as file ...
5
votes
3answers
2k views

Extract raster by mask output does not completely overlap the mask

I have a flow accumulation grid and a shapefile containing catchment areas. Now I want to do some calculations on them so I need to extract data from cells in the raster adjacent to certain catchment ...
2
votes
0answers
59 views

How do I make a path cross the most possible classes in a raster?

I want to create a function and possibly use Arcmap tools to create something similar to Cost Path. Instead of just choosing the LCP, it must cross through the most different integers. Let's say I ...
4
votes
1answer
332 views

Trouble Tabulate Area

I am using Kansas Cropland raster data from NASS with tabulate area to find out acreage of different crop types. The python I used works fine for 2006-2009 (56 meter cell size), except 2010 & ...
1
vote
2answers
619 views

arcpy.sa.Viewshed with OFFSETA

This is using ArcGIS10 arcpy with spatial analyst. I need to call Viewshed, but with an observer height set. I have tried using a Z-aware point feature set, but that only throws an error. I am ...
2
votes
1answer
396 views

How can you make RasterToASCII_conversion observe environment setting output cell size?

I have two rasters (an ESRI grid and an IMAGINE image) that differ in projection, extent and cell size. I want to convert the IMAGINE to an ASCII, and want it to have the extent, cell size and ...
2
votes
0answers
544 views

Spatial Analyst Cell Statistics seems to give the wrong answer in ArcGIS 10.0

I have a script that calculates standard deviation for a list of grids using the Cell Statistics function in Spatial Analyst for ArcGIS 10.0 The syntax is this: stdevGrid = ...
4
votes
1answer
5k views

Thiessen Polygon delineation within a Feature

I am trying to create Thiessen (Voronoi) polygons, based on points within a defined shapefile extent. The output should be inclusive to the feature. The output I receive now is the rectangular area ...
4
votes
1answer
1k views

Processing time using ArcGIS Hydrology toolbox tools in standalone Python script vs. ArcCatalog

My question is: what might cause the ArcGIS 9.3 Hydrology Toolbox Snap Pour Point and Watershed tools to run significantly slower in a Python script versus manually running the tools in ArcCatalog ...
7
votes
2answers
815 views

Process a raster with multiple geoprocessing tools in arcgis 10 using python - flow between processes

I'm writing a short python script that includes processing rasters with different geoprocessing tools in ArcGis 10, for example, the three processes below, is this best way to use multiple tools ...
5
votes
2answers
858 views

Using an integer variable in where-clause in arcgis python script

I've seen some questions on this, but there seem to be many different answers... and none have seemed to work. I'm trying to include an integer variable in the SQL where-clause of a tool like this ...