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

Question

我正在制作一个接收三个向量的函数，使用 akima 对它们进行插值并使用 plot_ly() 绘制它们。尽管通用代码有效，但我遇到了 interp() 输出的 z-matrix 缩放问题。

举个例子：

x 是一个包含一些重复值的 non-NA 数字。
y 是一个包含一些重复值的 non-NA 数字。
z 是一个 non-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 矩阵不能直接从 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')

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

Interpolation (akima) omits part of data when x/y contains duplicate elements

3d

interpolation

r

plotly