r语言 - 添加平面多边形到栅格可视

Question

我想创建一个结合3D水深网格数据(海洋)和2D多边形(陆地)的图形。

例如，在这样一个图(http://chenangliu.info/en/wp-content/uploads/2014/02/wireframe.jpg)中，我想在紫色区域添加一个平面，以及点，文本，…

我注意到rasterVis包，它完美地绘制了我的测深光栅。然而，我不知道如何添加陆地，就像普通的2D多边形一样。我尝试通过将陆地面积的网格值设置为零来解决这个问题，但结果并不完美，因为边界不清晰，有湖泊，……不应该是3D的

所以这个很好:

library(rasterVis)
pal<-colorRampPalette(c("darkblue","lightblue","green"))
    #----this doesn't help
     bathy[bathy>0,]=0
    #----
plot3D(bathy,col=pal)

这并不:

polygon3d(europa)  # error: Error: n > 2 is not TRUE

最后一个函数不起作用[idem for polygon()]，就像我的其他尝试一样。我想这并不难，但我找不到任何这样的例子。所以欢迎任何提示

Answer 1

您必须根据其规格使用polygon3d。由于您没有提供数据来复制您的代码，因此我使用plot3D帮助页面中的示例:

library(raster)
library(rasterVis)
library(rgl)
data(volcano)
r <- raster(volcano)
extent(r) <- c(0, 610, 0, 870)
plot3D(r)

现在你可以用rgl包中定义的函数添加一个3D多边形:

x <- c(30, 400, 400, 30)
y <- c(30, 30, 500, 500)
polygon3d(x, y, z=rep(100, 4), col = 'red')
polygon3d(x, y, z=rep(120, 4), col = 'blue')
polygon3d(x, y, z=rep(140, 4), col = 'black')