時系列ラスターデータの読み込み（R）

以前過去記事（時系列ラスターデータの処理について）ではでndvitsを用いたけどもかなり古いし、ライブラリがアーカイブされてしまったので、新たな取組を。

１．複数枚のtifファイルを複数バンドのもつ一つのtifにしてstackとして読み込む。
過去記事（大量のラスターデータを一枚のtifにまとめる）参照
ただし、GRASSGISを使ってます。

２．複数枚のtifファイルをそのままRでstack化する。
setwd("~/GIMMS/TIF/")
ls=grep(".tif",dir("~/GIMMS/TIF/"),value=T)
rasStack=stack(ls)
ただし、ls変数内のファイルの順序に注意。

簡単！！