R绘图面，包含NA值

Question

我想在R中绘制一个3d曲面。我的数据包含NA。

# Topographic Surface Plot wanted!
library(rgl)
x<-seq(from = 0,to = 4000,by = 500)
y<-seq(from = 10,to = 580,by = 60)
z<-matrix(data = rbinom(n = 9*10,size = 10,prob = 0.7),nrow = length(y))
# Everybody can plot this in 3d, but my real world data contains empty space NA values like here:
z[8:10,8:9]<-NA
# Also, the y-axis is not homogeneous, it looks more like this:
y<-y+rbinom(n = length(y),size = 15,prob = 0.4)

需要帮助:如何绘制xyz地形表面？

提前感谢！

Answer 1

除非您正在尝试生成3D热图，否则看起来您的值太多了。如果要生成曲面图，则Z应为数组，而不是矩阵。再加上现在的情况，Z是以NA值输出的。这是因为z[8:10,8:9]<-NA这一行。

试试这里的建议：R: Plotting a 3D surface from x, y, z

Answer 2

您可以尝试使用persp在3D透视图中绘制曲面--它还允许NAs和不均匀分布的x和y值：

# theta and phi control view angle    
persp(x=y,y=x,z=z, theta=50, phi = 40, zlim=c(0,15))

# for fancy colours:
nr <- nrow(z)
nc <- ncol(z)

# Calculate value at center of each cell
zfacet <- (z[-1, -1] + z[-1, -nc] + z[-nr, -1] + z[-nr, -nc])/4

# Generate the desired colors
cols = heat.colors(10)

# Cut matrix values into 10 bins by manual breaks
zbinned <- cut(zfacet, breaks=10)

persp(x=y,y=x,z=z, theta=30, phi = 30, r=2, zlim=c(-15,15), col=cols[zbinned], ticktype='detailed')