У меня есть растровый стек, содержащий 6 слоев, каждый из которых связан с годом данных NDVI.Кроме того, у меня есть пространственная точка данных, которая состоит из точечных данных за те же годы.Я хотел бы извлечь значение NDVI для каждой точки, но я бы хотел сопоставить данные по годовой точке с соответствующим годом данных NDVI.Например, точки 2008 года должны извлекать данные только из уровня NDVI 2008 года.
Я пробовал разные условные операторы вместе с функцией извлечения {растр} ... но мне не повезло.
#####NDVI raster stack
amp08<-raster("H:/ArcGIS Files/GIS Mule Deer layers/amp2008.tif")
amp09<-raster("H:/ArcGIS Files/GIS Mule Deer layers/amp2009.tif")
amp10<-raster("H:/ArcGIS Files/GIS Mule Deer layers/amp2010.tif")
amp11<-raster("H:/ArcGIS Files/GIS Mule Deer layers/amp2011.tif")
amp12<-raster("H:/ArcGIS Files/GIS Mule Deer layers/amp2012.tif")
amp13<-raster("H:/ArcGIS Files/GIS Mule Deer layers/amp2013.tif")
amp<-stack(amp08,amp09, amp10, amp11, amp12, amp13)
####spdf
loc = cbind(xx.clean$x2_, xx.clean$y2_)
proj4string <- CRS("+proj=utm +zone=11 +ellps=GRS80 +units=m +no_defs")
deerspdf <- SpatialPointsDataFrame(loc, xx.clean, proj4string =
proj4string)
deerspdf$year<-as.factor(deerspdf$year) ###"2008", "2009", "2010"....
deerspdf
@ data : 'data.frame': 2500 obs. of 4 variables:
$id: Factor with 45 levels
$year: Factor w 6 levels
$x2_: num [1:2500]
$y2_: num [1:2500]
В качестве конечного продукта мне нужен один столбецс извлеченными значениями NDVI для каждого года.