Rで解析：ラスタデータをデータキューブで確認「cubeview」パッケージ

ラスタデータをインタラクティブに3Dデータキューブで確認できるパッケージの紹介です。plot(ラスターデータ)で出力されるグラフを3Dでインタラクティブに確認することができます。

パッケージバージョンは0.2.0。実行コマンドはRStudioのRStudio Desktop 2021.09.2+382、windows 11のR version 4.1.3で確認しています。

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

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

#&#12497;&#12483;&#12465;&#12540;&#12472;&#12398;&#12452;&#12531;&#12473;&#12488;&#12540;&#12523;
install.packages("cubeview")

実行コマンド

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

#&#12497;&#12483;&#12465;&#12540;&#12472;&#12398;&#12452;&#12531;&#12473;&#12488;&#12540;&#12523;
library("cubeview")

###&#12487;&#12540;&#12479;&#20363;&#12398;&#20316;&#25104;_1#####
#R&#12525;&#12468;&#12398;jpeg&#30011;&#20687;&#12434;&#35501;&#12415;&#36796;&#12415;
Rlogo_data <- system.file("img", "Rlogo.jpg", package = "jpeg")
#raster::stack&#12395;&#22793;&#25563;
Rlogo_Stack <- raster::stack(Rlogo_data)
########

##&#12487;&#12540;&#12479;&#20363;&#12398;&#20316;&#25104;_1
#raster&#12487;&#12540;&#12479;&#12461;&#12517;&#12540;&#12502;&#12434;&#20316;&#25104;:cubeview&#12467;&#12510;&#12531;&#12489;
cubeview(Rlogo_Stack, col.regions = viridisLite::viridis)


###&#12487;&#12540;&#12479;&#20363;&#12398;&#20316;&#25104;_2#####
#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
########

##&#12487;&#12540;&#12479;&#20363;&#12398;&#20316;&#25104;_2
#raster&#12487;&#12540;&#12479;&#12461;&#12517;&#12540;&#12502;&#12434;&#20316;&#25104;:cubeview&#12467;&#12510;&#12531;&#12489;
cubeView(TESTRaster, col.regions = viridisLite::viridis)