r - 在 R 中绘制 3D 栅格(DEM)？

Question

我正在尝试使用 raster 包在 R 中绘制作为光栅导入的 DEM 的 3D 表面。

到目前为止，我的代码是:

DEM <- raster("DSM_TLS_2010_25cm_v4.tif")
DEM <- setMinMax(DEM)
col <- rainbow(20)
plot(DEM, col=col, zlim=c(0,790.22), main="Digital Elevation Model (DEM)")

这非常适合 DEM 的 2D 图，但是，当我尝试将其制作成 3D 时，使用:

plot3d(DEM)

或

surface3d(DEM)

它说 不能将类型“S4”强制转换为类型“double”的向量。

我确信答案非常简单，但我还没有设法解决我发现的类似问题。

Answer 1

raster() 函数返回 RasterLayer 对象，我怀疑 rgl 包中的 plot3d() 和 surface3d() 函数不(或不完全)支持 RasterLayer 对象。

但是，RasterVis 包中的 plot3D() 函数可以。试一试:

install.packages("rasterVis")
library(rasterVis)
plot3D(DEM)   # note: 3D not 3d