Извлечение среднего значения в пределах многоугольника для стека растров - PullRequest
0 голосов
/ 18 сентября 2018

У меня есть набор растров, которые я хотел бы перебрать, чтобы извлечь среднее значение в пределах многоугольника для каждого растра.Тем не менее, это займет навсегда.Могу ли я создать растровый кирпич из стека, чтобы сократить время обработки?Любые другие предложения?У меня есть 69 растров ... но я могу сократить их до 14, если мне действительно нужно.

iterate_raster <- function(grids_location, shp_location, func){
        #list files (in this case raster TIFFs)
        grids <- list.files(grids_location, pattern = "*.tif$")
        #read-in the polygon shapefile
        shp <- readShapePoly(shp_location)
        #create a raster stack
        st <- paste0(grids_location, "/")
        s <- stack(paste0(st, grids))
        ex <- extract(s, shp, fun=func, na.rm=TRUE, df=TRUE, weights = TRUE)
        return(data.frame(ex))
}

df <- iterate_raster(dep, spsh, mean)
...