提取具有与多边形不同的空间范围的栅格数据

时间:2019-01-21 03:03:34

标签: r loops extract polygon raster

有几个栅格文件具有不同的空间范围。需要在一个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)



}

0 个答案:

没有答案