Estoy modelizando la respuesta de un ave (alondra) a diversas variables climáticas y topográficas en el programa R. Utilizando la función predict he predicho la abundancia por píxel (ver código más abajo). Ahora me gustaría exportar mis tres rásters a ArcGIS 10.2? ¿Cómo puedo hacerlo?
¿Puedo añadir las coordenadas de referencia utilizando ArcGis o es fácil hacerlo en R? ¿Hay una forma más rápida de predecir la densidad? (Estoy feliz de mapa de densidad sin los errores inicialmente?
library(raster)
# Load the raster layers of the predictors for the entire landscape as a stack.
nsw.stack <- stack(list.files(pattern="tif$", full.names=FALSE))
names(nsw.stack) <- c("riv","lake","ndvi","temp","prec","dem")
summary(nsw.stack)
class : RasterStack
dimensions : 2129, 1905, 4055745, 6 (nrow, ncol, ncell, nlayers)
resolution : 0.008926991, 0.008926991 (x, y)
extent : 55.99709, 73.003, 43.99541, 63.00097 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0
names : riv, lake, ndvi, temp, prec, dem
A continuación, utilicé la función predecir para asignar los datos utilizando mi modelo mejor clasificado (t5)
skylark.psi <- predict(t5, type="state", newdata=nsw.stack)
plot(skylark.psi, ylab="psi - GLM")
> skylark.psi
class : RasterStack
dimensions : 2129, 1905, 4055745, 4 (nrow, ncol, ncell, nlayers)
resolution : 0.008926991, 0.008926991 (x, y)
extent : 55.99709, 73.003, 43.99541, 63.00097 (xmin, xmax, ymin, ymax)
coord. ref. : NA
names : Predicted, SE, lower, upper
min values : 2.987641e-05, 1.275741e-05, 1.293780e-05, 6.899167e-05
max values : 779.48502, 85.93247, 675.88093, 898.97033