Ihavenothing persp
和 persp3d
接受的是一个二维矩阵,
以下面那个数据集为例,不能直接使用,要先插值拟合,不知道还有没有其他现成的插值拟合方法
data(LIDAR,package="MBA")
> head(LIDAR)
x y z
1 711000.4 5093989 466.08
21 711032.1 5093990 467.28
41 711048.0 5093988 468.81
61 711058.7 5093973 469.08
81 711063.8 5093998 469.37
101 711071.2 5093983 469.09
mba.int <- MBA::mba.surf(LIDAR, 300, 300, extend=TRUE)$xyz.est
z <- mba.int$z
nrz <- nrow(z)
ncz <- ncol(z)
jet.colors <- colorRampPalette( c("gray100","gray0") )
nbcol <- 100
color <- jet.colors(nbcol)
zfacet <- z[-1, -1] + z[-1, -ncz] + z[-nrz, -1] + z[-nrz, -ncz]
facetcol <- cut(zfacet, nbcol)
persp(mba.int, theta =90, phi = 20,
col = color[facetcol], scale = FALSE,
shade = 0.75, expand = 20,border = NA, box = FALSE)