R：轮廓图正在下降点

Question

我正在绘制x，y和z维度中许多点的数据集。对于每个x和y，可以有多个z。我试图做的是将xy平面划分为50x50的网格，然后绘制等高线图中每个框的z的最大值。我认为我在网上发现的代码工作正常，但经过仔细检查后，看起来它正在丢弃点或倾斜数据。我的代码有问题吗？

df<-eff[c("Xs",Ys","Zs")]
df$x<-cut(df$Xs, breaks=50, labels=FALSE)
df$y<-cut(df$Ys, breaks=50, labels=FALSE)
df.max<-expand.grid(x=1:50, y=1:50)
df.max<-merge(df.max, aggregate(Zs~x+y, df, max), all.x=TRUE)
z<-matrix(df.max$Zs, nr=50, nc=50)
x.values <- min(df$Xs)+(0:49)*diff(range(df$Xs))/50
y.values <- min(df$Ys)+(0:49)*diff(range(df$Ys))/50
image2D(x=x.values, y=y.values, z, rasterImage = TRUE, contour = list(lwd = 2, col = jet.col(11)),
    main="Contour Plot", xlab="Xs",ylab="Ys", clab="Zs", ylim=c(0,max(df$Ys)*1.05), 
    xlim=c(0,max(df$Xs)*1.05))

根据xy平面的图（我认为是正确的），等高线图的轮廓应如下所示：

请注意，它与y轴相交约90k，并且最小值约为（25k，73k）。

来自相同数据的轮廓图如下所示：

我目前没有提供样本数据，因为我的模型正处于运行过程中，我不想全部启动它。如有必要，我可以稍后添加一个样本。

这是指向data

的链接