У меня есть SpatialPolygonsDataFrame со 120 полигонами и некоторыми связанными данными. Теперь я хотел бы извлечь среднее значение значений в растре каждого полигона отдельно. Мне удалось построить отдельные полигоны с помощью:
plot(SpatialPolygons(SPdataframe@polygons)[i])
Но не удалось извлечь значения таким же образом:
extract(raster, SpatialPolygons(SPdataframe@polygons)[i],fun="mean",na.rm=TRUE,method="simple")
Может кто-нибудь объяснить разницу между использованием одной и той же индексации в этих двух случаях? Каков официальный способ выбора отдельных полигонов элемента SpatialPolygonsDataFrame с индексами?
Заранее большое спасибо за помощь!