Zonal Estadísticas para los Polígonos con R

Question

Me gustaría hacer zonal statistics / análisis de terreno con el polígono superposiciones - he logrado conseguir min / max así:

library(raster)
library(rgdal)

alt <- getData('alt', country = "AT")
gadm <- getData('GADM', country = "AT", level = 2)
gadm_sub <- gadm[1:3, ]

plot(alt)
plot(gadm_sub, add=T)

extract(alt, gadm_sub, fun = max, na.rm = T, small = T, df = T)

Pero para la pendiente / aspecto me sale este error que yo no entender cómo prevenir:

extract(alt, gadm_sub, fun = function(x) terrain(x, opt = "aspect"), na.rm = T, small = T, df = T)
# > Fehler in FUN(X[[1L]], ...) : unbenutzte(s) Argument(e) (na.rm = TRUE)

Answer 1

Sry que era demasiado fácil.. Para el registro lo voy a publicar aquí:

    library(raster)
    library(rgdal)

    alt <- getData('alt', country = "AT")
    gadm <- getData('GADM', country = "AT", level = 2)
    gadm_sub <- gadm[1:3, ]

    plot(alt)
    plot(gadm_sub, add=T)

    asp <- terrain(alt, opt = "aspect", unit = "degrees", df = F)
    slo <- terrain(alt, opt = "slope", unit = "degrees", df = F)

    extract(slo, gadm_sub, fun = mean, na.rm = T, small = T, df = T)
    extract(asp, gadm_sub, fun = mean, na.rm = T, small = T, df = T)