有几个栅格文件具有不同的空间范围。需要在一个shapefile中将平均值提取到多个多边形。输出应具有面属性和提取的栅格文件名作为单独的列。这是经过测试的代码,但无法在提取的输出中将每个属性作为一列获取。
library(rgdal)
library(raster)
library(rasterVis)
## Create a list with raster with different spatial extent
stack_list <- list.files(pattern=".tif$", full.names=T)
stack_list
## Loop through each raster and extract the mean polygon value
for (i in 1:length(stack_list)) {
ras <- raster(stack_list[i])
#print(ras@crs)
plot(ras)
poly_d <- extract(ras, poly, weights=FALSE, fun=mean)
}