在r中将地图叠加在3D曲面地图之上

我使用 rgl.surface() 创建了一个 3d 地图，主要遵循 Shane 在中的回答this https://stackoverflow.com/questions/1896419/plotting-a-3d-surface-plot-with-contour-map-overlay-using-r邮政。使用我自己的数据，我得到了这张地图

在此表面地图的顶部，我想添加植被密度地图，以便我获得类似这样的东西（通过软件 Surfer 获得）：

是否可以使用 rgl 来做到这一点，或者 r 中的任何其他包，或者是像 Shane 的答案那样拥有两个地图的唯一解决方案？

谢谢。

Edit:

根据@gsk3的请求，以下是该地图的代码：

library(rgl)

# Read the z (i.e. elevation) dimension from file
z1 = matrix(scan("myfile.txt"),nrow=256, ncol=256, byrow=TRUE)
#create / open x y (i.e. easting and northing coordinates) dimensions 
y=8*(1:ncol(z)) # Each point is 8 m^2
x=8*(1:nrow(z))

# See https://stackoverflow.com/questions/1896419/plotting-a-3d-surface-plot-with-contour-map-overlay-using-r for details of code below
zlim <- range(z)
zlen <- zlim[2] - zlim[1] + 1
colorlut <- terrain.colors(zlen,alpha=0) # height color lookup table
col <- colorlut[ z-zlim[1]+1 ] # assign colors to heights for each point
open3d()
rgl.surface(x,y,z)

我无法发布海拔代码，因为有 65536 个（即 x*y=256*256）点，但它是一个看起来像这样的矩阵

            [,1]     [,2]     [,3]     [,4]     [,5] 
[1,]    1513.708 1513.971 1514.067 1513.971 1513.875 
[2,]    1513.622 1513.524 1513.578 1513.577 1513.481

等等。 植被密度图也是如此，它的格式完全相同，并且每个 x*y 点都有一个值。我希望这能让事情变得更清楚......？

编辑2，最终版本

这是我用 R 制作的地图。我还没有在上面添加图例，但这是我稍后会做的事情。

最终的代码是

library(rgl)
z1 = matrix(scan("myfile.txt"),nrow=256, ncol=256, byrow=TRUE)
# Multiply z by 2 to accentuate the relief otherwise it looks a little bit flat.
z= z1*2

#create / open x y dimensions
y=8*(1:ncol(z))
x=8*(1:nrow(z))

trn = matrix(scan("myfile.txt"),nrow=256, ncol=256, byrow=TRUE)
fv = trn*100
trnlim = range(fv)

fv.colors = colorRampPalette(c("white","tan4","darkseagreen1","chartreuse4")) ## define the color ramp
colorlut =fv.colors(100)c(1,seq(35,35,length.out=9),seq(35,75,length.out=30),seq(75,100,length.out=61))] 

# Assign colors to fv for each point
col = colorlut[fv-trnlim[1]+1 ] 
open3d()
rgl.surface(x,y,z,color=col) 

非常感谢@gsk3 和@nullglobthis https://stackoverflow.com/questions/6981631/sharp-break-in-color-palette-in-r发帖寻求他们的帮助。希望这篇文章能帮助更多人！

修改上面的代码来给出答案。请注意，地形应该是与高程矩阵格式相同的矩阵。我添加了一个,color函数调用的参数，因此它实际上使用您创建的颜色矩阵。

library(rgl)

# Read the z (i.e. elevation) dimension from file
z1 = matrix(scan("myfile.txt"),nrow=256, ncol=256, byrow=TRUE)
#create / open x y (i.e. easting and northing coordinates) dimensions 
y=8*(1:ncol(z)) # Each point is 8 m^2
x=8*(1:nrow(z))

# Read the terrain types from a file
trn = matrix(scan("terrain.txt"),nrow=256, ncol=256, byrow=TRUE)

# See http://stackoverflow.com/questions/1896419/plotting-a-3d-surface-plot-with-contour-map-overlay-using-r for details of code below
trnlim <- range(trn)
trnlen <- trnlim[2] - trnlim[1] + 1
colorlut <- terrain.colors(trnlen,alpha=0) # height color lookup table
col <- colorlut[ trn-trnlim[1]+1 ] # assign colors to heights for each point
open3d()
rgl.surface(x,y,z,color=col)