【R】persp

perspは、xy平面上の透視投影図を描くbase-r標準の関数なのですが、今さら知ったので、使ってみた。

おなじみのvolcanoデータです。

image(volcano)

等高線で表示してみます。

contour(volcano)

で、persp関数を使って、こんな感じで表現できます。

persp(volcano,  expand=0.4, shade=0.2, theta=12.2, phi=23.4)

色も付けられます。

persp(volcano, col = rainbow(30), expand=0.4, shade=0.2, theta=12.2, phi=23.4)

par(bg = "white")
x <- seq(-1.95, 1.95, length = 30)
y <- seq(-1.95, 1.95, length = 35)
z <- outer(x, y, function(a, b) a*b^3)
nrz <- nrow(z)
ncz <- ncol(z)
jet.colors <- colorRampPalette( c("green", "white", "yellow") )
nbcol <- 20
color <- jet.colors(nbcol)
zfacet <- z[-1, -1] + z[-1, -ncz] + z[-nrz, -1] + z[-nrz, -ncz]
facetcol <- cut(zfacet, nbcol)
persp(x, y, z, col = color[facetcol], phi = 25, theta = -20)