Hot answers tagged raster
In this test raster (20x20), with values between 1 and 50, I replaced 19 value by -32768 value. In the following raster, the -32768 values were replaced by -3.40282e+38 values (see window of Value Tool Plugin): The used expression is in the window of "Raster calculator expression" (below image and code sample): ...
Both variables are zonal means. The average distance to the nearest facility is the zonal mean of the Euclidean distance grid (based on the facilities). The average number of facilities is the zonal mean of a one-kilometer radius focal sum of the facilities grid. (This is merely a grid whose cell values count the number of facilities within each cell. ...
You're trying to do two things at once: reclassify the values of r2 and then multiply those by r1. Instead, do them separately: d <- 1:3/4 out = overlay(r1, calc(r2, fun=function(i) d[i]), fun="*") It is, of course, your responsibility to ensure that the values of r2 are all valid indexes into array d.
whuber showed the way, but here's another way to get there (using the functions that match the operations he suggests) d <- c(.25,.50,.75) m <- cbind(1:3, d) r3 <- reclassify(r2, m) out <- r1 * r3
Ok Here is what I would do. Use the "Extract Values to Points" tool under "Extraction" which is in "Spatial Analyst Tools". You will use your point shapefile as the Input Point Feature and your continuous raster as the input Raster. I think that may be all you need to do. BTW instead of a 1-word comment, sometimes it is good policy to edit your question. ...
I see two fixes: For the first error: ERROR 000161: The length of the grid name must not exceed 13 characters Convert your output to tiff format, which does not have the character length limitations that the Esri grid raster format does. mamcnty_out = os.path.join(output_workSpace, "cnty_" + species_str + ".tif") For the second error: ERROR ...
Only top voted, non community-wiki answers of a minimum length are eligible