错误：增加 persp {fields} 中预期的“x”和“y”值

Question

我在lat, lng 的空间中有 500 个点：

lng=runif(100,1,4)
lat=runif(100,40,40.1)

假设我们有（密度）作为

z=rnorm(1000)

我尝试使用 persp 绘制 3d 曲面

persp(lng,lat,z)

但我有这个错误：

Error : increasing 'x' and 'y' values expected

我使用 x 和 y 排序

data=data.frame(lng,lat)
data=data[ order(data[,1],data[,2]), ]

我又试了一次

persp(data[,1],data[,2],z)

但我还是有同样的问题。

我该怎么做？

谢谢。

Answer 1

问题在于参数z 以及x 和y 的顺序。正如阴影在 cmets 中提到的，这需要是 dim(z) = length(x) x length(y) 在这种情况下为 100x100 的矩阵。显然 x 和 y 需要both同时排序，就像上面提到的阴影一样，这在persp.default as.well 的源代码中可以看到：

if (any(diff(x) <= 0) || any(diff(y) <= 0)) #both need to be ordered
        stop("increasing 'x' and 'y' values expected")

所以下面的工作：

lng=runif(100,1,4)
lat=runif(100,40,40.1)
z=matrix(rnorm(10000),ncol=100,nrow=100) #use the correct z as matrix
data=data.frame(lng,lat)
data=data[ order(data[,1],data[,2]), ]   #just order lng as you do 
data$lat <- sort(data$lat)               #you need to sort lat too

persp(data[,1] ,data[,2], z) 

现在可以了。

编辑

我认为您真正想要做的是让您在 z 中获得 lan/lng 坐标并绘制它们。 z 参数需要携带坐标。 x 和 y 参数分别只携带 x 和 y 轴的值。您可以使用默认的 x 和 y 参数（即，将它们留空）或仅对 lan 和 lng 分别排序以提供实际值。因此，您正在寻找的解决方案是：

lng=runif(100,1,4)
lat=runif(100,40,40.1)
z= as.matrix(data.frame(lng=lng,lat=lat))
z=matrix(rnorm(10000),ncol=100,nrow=100) #z contains lan and lng

#plot z with default x and y values (or use sorted lng and lat instead)
persp( z = z)  

哪些输出：