当x / y包含重复元素时，插值（akima）省略部分数据

Question

我正在创建一个接收三个向量的函数，使用akima对它们进行插值，并使用plot_ly（）绘制它们。虽然通用代码有效，但我遇到了interp（）输出的z矩阵缩放问题。

让我给你举个例子：

• x是包含一些重复值的非NA数字。
• y是包含一些重复值的非NA数字。
• z是非NA连续向量

一些摘要统计：

> unique(x)
[1] 60 48 36 32 18 24 30 15 12 28 21 19 54 20 16 27 10 39 14 17  9  6 50  8 13
> range(x)
[1]  6 60

> unique(y)
[1] 10.00 10.50 13.50 12.50 14.00 12.00 11.00  9.00 11.50  9.25 13.00 10.25 6.50 6.75 8.25 9.50
[17]  8.00 8.85  9.75  7.90  7.00  8.60  8.75  7.50  8.90  8.50  7.49  7.40  5.50  7.60  7.25  8.35
[33]  6.00  5.00  7.75  7.35  6.30  4.50  5.75  8.40  5.60  5.90  7.74  9.90  6.20  5.80
> range(y)
[1]  4.5 14.0


> head(z)
[1] 2.877272 3.267328 3.175478 3.843326 4.809792 2.827825
> range(z)
[1]  2.316529 28.147808

我在下面实现基线功能：

labs = list(x = 'x', y = 'y', z = 'z')

mat = interp(x, y, z, duplicate = 'mean', extrap = T, xo = sort(unique(x)))

plot_ly(x = mat$x, y = mat$y, z = mat$z, type = 'surface') %>%
    layout(title = title,
           scene = list(xaxis = list(title = labs$x),
                        yaxis = list(title = labs$y),
                        zaxis = list(title = labs$z)))

当我运行它时，输出如下：

问题是此图片中未涵盖部分数据。例如，存在大约x> 50，y <11的相当大的数据部分，其被插值省略（因此未被绘制）。

length(x[x > 50])
[1] 304
> length(y[x > 50 & y < 11])
[1] 290
> length(z[x > 50 & y < 11])
[1] 290

我怀疑这与重复的x值有关。因此，我在interp（）中配置了xo参数，以便：

mat = interp(x, y, z, duplicate = 'mean', xo = sort(unique(x)), decreasing = T)

在这种情况下，部分地绘制了先前省略的区域。它看起来如下：

尽管如此，x和y轴仍然不对应于它们各自的数据范围（尽管数据可用）。底线：如何调整函数，使曲面始终延伸到x和y的整个范围？

最好

Answer 1

事实证明，错误来自plot_ly（）。显然，z-matrix不能直接从interp（）传递到plot_ly（），因为轴会错误地传递给图形。因此，需要对内插的z矩阵进行变换。

如果组合使用这两个函数，请确保执行z的转换，如下所示：

mat = interp(x,y,z, duplicate = 'mean')

x = mat$x
y = mat$y
z = matrix(mat$z, nrow = length(mat$y), byrow = TRUE)

plot_ly(x, y ,z, type = 'surface')