绘制网格的 3D 地形图

Question

我想使用 R 复制类似下图的内容，可能使用 ggplot（尽管我怀疑这是可能的，因为 AFAIK，它没有 3D 功能）。我拥有的数据通常是栅格包中的栅格文件，但我可以将其转换为最合适的格式。

剧情取自： "Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change"，图1.14。我不知道那个情节是哪个软件制作的。

我想唯一的办法就是使用 lattice::cloud（线框）或类似的东西？我似乎找不到任何方法来强制线框使用条形图而不是曲面图；此外，基于海上高度的着色，网格保持水平，可能是不可能的..

Answer 1

我抽空看了看这个。最后我采用了以下方法（链接的 R 代码）：

http://pastebin.com/dA2nNNS0

感谢评论中为我指出正确方向的人。

我在保存文件时仍然遇到问题，但这是 material 另一个问题。

来自 pastebin 的代码：

library(raster) 
library(rgl)

r <- raster("topo12.nc", varname="topo") #Read file
cols <- terrain.colors(3100) #Define colors

binplot.3d <- function(x,y,z,alpha=1,topcol="#ff0000",sidecol="#aaaaaa"){ #Binplotting function
  save <- par3d(skipRedraw=TRUE)
  on.exit(par3d(save))

  x1<-c(rep(c(x[1],x[2],x[2],x[1]),3),rep(x[1],4),rep(x[2],4))
  z1<-c(rep(0,4),rep(c(0,0,z,z),4))
  y1<-c(y[1],y[1],y[2],y[2],rep(y[1],4),rep(y[2],4),rep(c(y[1],y[2],y[2],y[1]),2))
  x2<-c(rep(c(x[1],x[1],x[2],x[2]),2),rep(c(x[1],x[2],rep(x[1],3),rep(x[2],3)),2))
  z2<-c(rep(c(0,z),4),rep(0,8),rep(z,8) )
  y2<-c(rep(y[1],4),rep(y[2],4),rep(c(rep(y[1],3),rep(y[2],3),y[1],y[2]),2) )
  rgl.quads(x1,z1,y1,col=rep(sidecol,each=4),alpha=alpha)
  rgl.quads(c(x[1],x[2],x[2],x[1]),rep(z,4),c(y[1],y[1],y[2],y[2]),
              col=rep(topcol,each=4),alpha=1) 
  rgl.lines(x2,z2,y2,col="#000000")
}

cat("Row ( of", dim(r)[1],"):")
for (row in 1:dim(r)[1]) { #Plotting loop
    for (col in 1:dim(r)[2]) {
    if (round(r[row, col]) < 1) {
    binplot.3d(c(col-1,col), c(row-1,row), r[row, col]/500, alpha=1, topcol="cadetblue3")
    } else {
    binplot.3d(c(col-1,col), c(row-1,row), r[row, col]/500, alpha=1, topcol=cols[round(r[row, col])])
#   cat(round(r[row, col]), "\t")
    }
    }
    cat(row, "")
}

绘制网格的 3D 地形图

Plot 3D orography of a grid

plot

r

ggplot2