从粗分辨率网格插入到更精细分辨率网格

Question

我试图将温度 (Tas) 从粗网格 (2.5° x 2.5°) 插入到更精细的分辨率网格 (0.5° x 0.5°)。 我的数据采用以下形式：

COARSE = matrix(c(rep(seq(-180, -170, by = 2.5), times = 5), 
            rep(seq(40, 50, by = 2.5), each = 5), 
            rep(seq(6, 10, by = 1), times = 5)),
            nrow = 25, ncol = 3)
colnames(COARSE) = c("Lon", "Lat", "Tas")

FINE = matrix(c(rep(seq(-180, -170.5, by = 0.5), times = 20), 
              rep(seq(40, 49.5, by = 0.5), each = 20), 
            rep(seq(5, 9, by = 1), times = 80)),
            nrow = 400, ncol = 3)
colnames(FINE) = c("Lon", "Lat", "Tas")

输出应该是一个矩阵，前两列中包含较细网格的 Lon 和 Lat，然后第三列应表示插值到较细网格上的粗网格 Tas 值。

我曾尝试在 Akima 包中使用双线性函数，但我无法理解它如何工作得足够好以适应我的情况。 当然，上面的例子真的很容易插值吗？ 在过去的几天里，我搜索了很多例子，并被难住了。 任何建议都非常感谢。

Answer 1

也许您可以尝试这样的插值，其中使用包pracma的二维插值interp2 ：

library(pracma)
FINE[,3] <- interp2(unique(COARSE[,1]),
                    unique(COARSE[,2]),
                    matrix(COARSE[,3],nrow = length(unique(COARSE[,1]))),
                    unique(FINE[,1]),
                    unique(FINE[,2]))

以至于

> head(FINE,20)
         Lon Lat Tas
 [1,] -180.0  40 6.0
 [2,] -179.5  40 6.2
 [3,] -179.0  40 6.4
 [4,] -178.5  40 6.6
 [5,] -178.0  40 6.8
 [6,] -177.5  40 7.0
 [7,] -177.0  40 7.2
 [8,] -176.5  40 7.4
 [9,] -176.0  40 7.6
[10,] -175.5  40 7.8
[11,] -175.0  40 8.0
[12,] -174.5  40 8.2
[13,] -174.0  40 8.4
[14,] -173.5  40 8.6
[15,] -173.0  40 8.8
[16,] -172.5  40 9.0
[17,] -172.0  40 9.2
[18,] -171.5  40 9.4
[19,] -171.0  40 9.6
[20,] -170.5  40 9.8