Извлечение растровых данных с разным пространственным экстентом в полигон - PullRequest
Новая
       36

Извлечение растровых данных с разным пространственным экстентом в полигон

0 голосов
/ 21 января 2019

Существует несколько растровых файлов с разным пространственным экстентом.Необходимо извлечь средние значения для нескольких полигонов в одном шейп-файле.Выходные данные должны иметь атрибуты многоугольника и имя извлеченного растрового файла в виде отдельных столбцов.Это проверенный код, но он не может получить каждый атрибут в виде столбца в извлеченном выводе. 

   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)



}
...