向rasterVis中添加平面多边形

Question

我想要创建一个图形，将三维测深网格数据(sea)与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

最后一个函数对多边形()不起作用，就像我的其他尝试一样。我想这不是太难，但我不能详细说明这方面的任何例子。所以任何暗示都是值得欢迎的。

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')