Rで解析：ラスタデータの表現に便利です！「rasterVis」パッケージ

ラスタデータは海面温度やデジタル写真、データ分布などの表現に使用される座標と数値をもつデータ形式です。座標と数値をもつ表現にはヒートマップなどがあります。ラスタデータの適切な表現は解釈や物事を伝える上で非常に有用な方法かと思います。

ラスタデータを表現するのに有用な「rasterVis」パッケージを紹介します。海面温度のデータを用意し、今話題のエルニーニョなどを表現してみてはいかがでしょうか。plot3Dコマンドはインタラクティブに回転などが可能です。ぜひ試してみてください。

パッケージのバージョンは0.51.2。実行コマンドはwindows 11のR version 4.1.3で確認しています。

パッケージのインストール

下記コマンドを実行してください。

#パッケージのインストール
install.packages("rasterVis")

実行コマンドの紹介

詳細はコマンド、パッケージヘルプを確認してください。

#&#12497;&#12483;&#12465;&#12540;&#12472;&#12398;&#35501;&#12415;&#36796;&#12415;
library("rasterVis")

###&#12487;&#12540;&#12479;&#20363;&#12398;&#20316;&#25104;#####
#raster&#12497;&#12483;&#12465;&#12540;&#12472;&#12364;&#12394;&#12369;&#12428;&#12400;&#12452;&#12531;&#12473;&#12488;&#12540;&#12523;
if(!require("raster", quietly = TRUE)){
  install.packages("raster");require("raster")
}
#raster&#12524;&#12452;&#12516;&#12540;&#12434;&#20316;&#25104;
TESTRaster <- raster(nrow = 100, ncol = 100, 
                     xmn = 0, xmx = 900, ymn = 0, ymx = 900)
#raster&#12524;&#12452;&#12516;&#12540;&#12395;&#12487;&#12540;&#12479;&#12434;&#20184;&#19982;
set.seed(1234)
values(TESTRaster) <- sample(c(1:100, 1500, 5000), 100 * 100,
                             replace = TRUE)
#&#20999;&#12426;&#25244;&#12365;&#31684;&#22258;&#12434;&#25351;&#23450;
CellExtent <- cellsFromExtent(TESTRaster, extent(200, 700,
                                                 250, 630))
#&#31684;&#22258;&#12434;NA&#12391;&#32622;&#12365;&#25563;&#12360;
TESTRaster[CellExtent] <- NA
########

#&#12521;&#12473;&#12479;&#12540;&#12487;&#12540;&#12479;&#12398;&#12503;&#12525;&#12483;&#12488;:levelplot&#12467;&#12510;&#12531;&#12489;
#&#19979;&#37096;&#12398;&#12473;&#12465;&#12540;&#12523;&#34920;&#31034;:scales&#12458;&#12503;&#12471;&#12519;&#12531;
#&#12473;&#12465;&#12540;&#12523;&#20301;&#32622;:margin&#12458;&#12503;&#12471;&#12519;&#12531;
#&#12475;&#12523;&#33394;&#12398;&#35373;&#23450;:par.settings&#12458;&#12503;&#12471;&#12519;&#12531;
#rasterTheme, RdBuTheme, BuRdTheme, GrTheme, BTCTheme, PuOrTheme, streamTheme
#&#12362;&#12424;&#12403;rasterTheme&#12398;&#20351;&#29992;&#12364;&#21487;&#33021;
levelplot(TESTRaster, xlab = "TEST1", ylab = "TEST2",
          colorkey = TRUE, margin = TRUE,
          par.settings = rasterTheme(pch = 19, cex = 0.7,
                                     region = c("#00bfd4", "#6e5f72")))

#&#12521;&#12473;&#12479;&#12540;&#12487;&#12540;&#12479;&#12398;3D&#12503;&#12525;&#12483;&#12488;:plot3D&#12467;&#12510;&#12531;&#12489;
plot3D(TESTRaster, col = c("#00bfd4", "#6e5f72"))

#&#12521;&#12473;&#12479;&#12540;&#12487;&#12540;&#12479;&#12398;&#20986;&#29694;&#25968;&#12434;&#12503;&#12525;&#12483;&#12488;:histogram&#12467;&#12510;&#12531;&#12489;
histogram(TESTRaster)

出力例

・levelplotコマンド

・plot3Dコマンド

・histogramコマンド

少しでも、あなたの解析が楽になりますように！！