色々とパッケージがあるそうですが、、raster libraryを使えば簡単だということがわかりました。ただしncdf4パッケージが必要です。 rstr.brick NAMEはgdalinfoで確認できる。 以上！

levelplotと似てるとこもあるけど違うとこも。 よくごっちゃになる。。 #タイトルの追加 spplot( raster, main = list(label="title", cex = 3)) #色の変更 spplot( raster, col.regions=colorRampPalette(brewer.pal(9, 'Reds')))#値の区切り spplot( raste…

色の変更 →library(RColorBrewer)を使って col.regions=colorRampPalette(brewer.pal(9, "RdBu"))などを追加 colorRampPaletteをつけないと変な挙動になることあり（なんで？？）デフォルトであるLongtitude, Latitudeの表示を消す → xlab="", ylab="" を追…

以前過去記事（時系列ラスターデータの処理について）ではでndvitsを用いたけどもかなり古いし、ライブラリがアーカイブされてしまったので、新たな取組を。 １．複数枚のtifファイルを複数バンドのもつ一つのtifにしてstackとして読み込む。 過去記事（大量…