R绘图曲面，包含NA值

我想在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)

#需要地形表面图！
图书馆（rgl）
x除非您试图生成三维热图，否则看起来您的值太多了。如果要生成曲面图，Z应该是数组，而不是矩阵。加上现在，Z输出NA值。这是因为线
z[8:10,8:9]您可以尝试
persp
以三维透视图绘制曲面——它还允许
NA
s和不均匀分布的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')

#θ和φ控制视角
persp（x=y，y=x，z=z，θ=50，φ=40，zlim=c（0,15））
#对于花式颜色：
天然橡胶