r - How to extract raster values from a RasterStackTS object using a shapefile? -

i have rasterstack , need extract values every raster using shapefile.

i know it's easy "normal" rasterstack, have rasterstackts object. unfortunately must manage rasterstackts object, because i've obtained if previous computation.

here reproducible example, taken babak naimi's website.

library(rts)  # location of files path <- system.file("external", package="rts")  # list of raster files: lst <- list.files(path=path,pattern='.asc$',full.names=true)  lst  #create random shapefile on raster extent r_4_shp<- raster(lst[1]) shp <- rastertopolygons(r_4_shp, fun=function(x){x>0.77})  #visualize shp on raster plot(r_4_shp) plot(shp, add=t)  # creating rasterstack object r <- stack(lst)  # corresponding dates 4 rasters d <- c("2000-02-01","2000-03-01","2000-04-01","2000-05-01")   d <- as.date(d)    # creating rasterstackts object: rt <- rts(r,d)  rt  #try extract raster values shapefile sp_basin_mean <- extract(rt, shp) 

i found answer here. of raster functions cannot work rasterstackts object. raster object kept in slot @raster. indeed line runs, using data example i've posted before:

sp_basin_mean <- extract(rt@raster, shp)