データで遊ぶ

ステップバイステップ

片栗粉Ｒ～ラピエ・ドリーネ・ポドゾル！～

片栗粉Ｒとは

片栗粉Ｒは、フリーの統計解析環境「R」を使って、地形データ（ハイトマップ）を立体感のある角度や光の当てかたで描画する方法です。この絶妙な設定値にたどりつくまでにかなり試行錯誤しました。

PCゲーム「A列車で行こう9」では、立体感をともなってマップの全景を見る機能がないため「広さ」が実感しにくいという課題がある

セーブ・ロードの画面から256×256のビットマップをキャプチャし、各ピクセルのRGB値から高度（標高）を読み取る

サンプルデータ	ひしめきＭ
	ｴｱﾎﾟｰﾄＭ
	ﾒﾄﾛﾎﾟﾘﾀﾝＭ
	タ陽Ｍ
	あかりＭ
	ｴｷｽﾊﾟｧＭ．
	三日月島ＭＶ

その他のおもなシナリオマップについてはこちらから。
マップコンストラクションで自動生成した地形についてはこちら。

「片栗粉Ｒ」のRスクリプト
myakari=data.matrix(read.table("clipboard",h=0)) # 256x256 z <- 0.5 * myakari x <- 10 * (1:nrow(z)) y <- 10 * (1:ncol(z)) par(bg = "black") # 背景色を黒にする persp(x, y, z, theta = 115, phi = 22.5, col = "white", scale = FALSE, ltheta = 283.95 + 90, lphi = 21.76, shade = 0.75, border = NA, box = FALSE) # プロットされた図を右クリックしてコピーや保存をします # （高精細な画像にするには図のウィンドウを最大化してから）

「片栗粉Ｒ」のRスクリプト

myakari=data.matrix(read.table("clipboard",h=0)) # 256x256
z <- 0.5 * myakari
x <- 10 * (1:nrow(z))
y <- 10 * (1:ncol(z))
par(bg = "black") # 背景色を黒にする
persp(x, y, z, theta = 115, phi = 22.5, col = "white", scale = FALSE, ltheta = 283.95 + 90, lphi = 21.76, shade = 0.75, border = NA, box = FALSE)
# プロットされた図を右クリックしてコピーや保存をします
# （高精細な画像にするには図のウィンドウを最大化してから）

※諸般の事情に鑑み、マップによらず変数名を「akari」にしています。Windowsでの実行を念頭にしていますので入出力に「クリップボード」を使用しています。

統計解析環境「R」に読み込んでプロットした結果を以下に示す

ゲーム中には気づかなかった大域的な特徴や細い谷筋などの微地形、整地可能な高台などが発見しやすくなることがわかった
マップの全景から受ける感覚的な「広さ」（距離感）が地形表現の詳細さ（きめ細かさ）や水面の分布（ひとまとまりの水面の面積や位置、方角）に依存していることが示唆される
このゲームではマップの広さを「10km四方」と公称しているが、ここでは地形および学習内容にあわせて適当な距離感をじぶんで定義することとする

【注意】PCゲーム「A列車で行こう9」は「地球の丸み」を表現しない。「A列車で行こう9」の地形データをGISや数値地図の可視化アプリなどに読み込ませるのは妥当ではない。Rのように愚直な3次元プロットを行なうのが正しい。あなたの学校や職場で数値計算や統計解析の商用ソフトウェアへのアクセスがあるなら、使い勝手やマニュアルが優れているのでまずはそちらを使ってほしい。なければRを使ってほしい。

ここでいう片栗粉とは、~~ガイドさんとあるく~~柱状節理を説明する演示実験の愛称です。東京・渋谷のカタクリコホットライン・ゴールドラッシュとは無関係です。

スーパーカタクリコＲ

スーパーカタクリコＲとは

スーパーカタクリコＲは、片栗粉Ｒのカラー版です。地形データの標高に応じて色を着けます。不必要に凝ったことはせず、フリーの統計解析環境「R」に組み込まれた標準の色を使います。

~~熱した溶岩で溶岩焼きステーキ…あつっ！！~~白黒の画像では味気ない・~~お能をなさっておられたとか~~厚みを表現しないと屋台のお面（の裏側）のようできもちわるい

~~バウムクーヘンに生クリームと栗のジャムをつけ~~着色する

「スーパーカタクリコＲ」のRスクリプト
myakari=data.matrix(read.table("clipboard",h=0)) # 256x256 install.packages('plot3D') # これは最初に1回だけ必要な操作です library(plot3D) # 起動するたびに必要 z <- 0.5 * myakari x <- 10 * (1:nrow(z)) y <- 10 * (1:ncol(z)) par(bg = "white") # 背景色を白にする persp3D(x, y, z, zlim = c(-80, 240), theta = 115, phi = 22.5, scale = FALSE, ltheta = 283.95 + 90, lphi = 21.76, shade = 0.75, border = NA, box = FALSE, clim = c(-320, 260), curtain = TRUE) # プロットされた図を右クリックしてコピーや保存をします

「スーパーカタクリコＲ」のRスクリプト

myakari=data.matrix(read.table("clipboard",h=0)) # 256x256
install.packages('plot3D') # これは最初に1回だけ必要な操作です
library(plot3D) # 起動するたびに必要
z <- 0.5 * myakari
x <- 10 * (1:nrow(z))
y <- 10 * (1:ncol(z))
par(bg = "white") # 背景色を白にする
persp3D(x, y, z, zlim = c(-80, 240), theta = 115, phi = 22.5, scale = FALSE, ltheta = 283.95 + 90, lphi = 21.76, shade = 0.75, border = NA, box = FALSE, clim = c(-320, 260), curtain = TRUE)
# プロットされた図を右クリックしてコピーや保存をします

	統計数理研究所のミラーサイトから「R」を入手しよう統計数理研究所とは