我正在尝试使用 rasterVis::plot3D 绘制一些 3D 数字高程模型。
DEM最初是geo tifs,它们都是同一位置的所有DEM,只是该位置的不同区域,它们都具有相同的坐标参考系(“+proj=utm +zone=21 +south +datum=WGS84 +units= m +no_defs"),唯一不同的是空间范围和分辨率。一个是面积约。 2 米分辨率下为 200 x 700 米,其余为 1 厘米分辨率下的 1 x 1 米。 x 和 y 值是坐标,z 是高程。当以 2D 方式绘制它们时,它们都绘制得很好并且看起来很正常。
但是,当我使用 rasterVis::plot3D 绘制更大面积的 2m 分辨率 DEM 时,我得到了一个漂亮的 3D 图 正如预期的那样,但是当我绘制任何更小、更高分辨率的 DEM 时,我得到了一个非常奇怪的输出 .
我不知道是什么导致了这个奇怪的输出。
Onedrive 链接下载 geo tif 文件 这就是 Q1cm.tif 在 2D 中绘制时的样子(正方形是塑料样方 - 我们在样方内进行了植被调查)。
library(raster)
library(rasterVis)
library(rgl)
dem2m<-raster("lidar.tif")
dem1cm<-raster("Q1cm.tif")
dem2m
class : RasterLayer
dimensions : 197, 292, 57524 (nrow, ncol, ncell)
resolution : 2, 2 (x, y)
extent : 360047, 360631, 3080519, 3080913 (xmin, xmax, ymin, ymax)
crs : +proj=utm +zone=21 +south +datum=WGS84 +units=m +no_defs
source : lidar.tif
names : lidar
values : 18.7875, 59.18345 (min, max)
dem1cm
class : RasterLayer
dimensions : 314, 336, 105504 (nrow, ncol, ncell)
resolution : 0.01, 0.01 (x, y)
extent : 360091.7, 360095, 3080597, 3080601 (xmin, xmax, ymin, ymax)
crs : +proj=utm +zone=21 +south +datum=WGS84 +units=m +no_defs
source : memory
names : Q1cm
values : 54.84027, 55.15403 (min, max)
rgl::open3d()
plot3D(dem1cm)
关于正在发生的事情或如何解决这个问题有什么想法吗?
提前致谢!