シェープファイルを使用してrasterstackすべてのラスターから値を抽出する必要があります。
「通常の」ラスタースタックを使用すると非常に簡単なことはわかっていますが、RasterStackTSオブジェクトがあります。残念ながらRasterStackTS、以前の計算から if を取得したため、オブジェクトを管理する必要があります。
これは、Babak Naimi のウェブサイトから取った再現可能な例です。
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 a random shapefile over the raster extent
r_4_shp<- raster(lst[1])
shp <- rasterToPolygons(r_4_shp, fun=function(x){x>0.77})
#visualize the shp over the raster
plot(r_4_shp)
plot(shp, add=T)
# creating a RasterStack object
r <- stack(lst)
# corresponding dates to 4 rasters
d <- c("2000-02-01","2000-03-01","2000-04-01","2000-05-01")
d <- as.Date(d)
# creating a RasterStackTS object:
rt <- rts(r,d)
rt
#try to extract raster values with the shapefile
sp_basin_mean <- extract(rt, shp)