沿2D矩阵的一维插值

Question

我正在尝试对M*N矩阵进行插值。 例如，矩阵T代表温度，矩阵B代表位置：

T = [1 3 5; ...
     2 4 6; ...
     1 2 3];
B = [0.1 0.2 0.3; ...
     0.1 0.2 0.3; ...
     0.1 0.2 0.3];

我正在尝试获取位置0.15和0.25处的温度信息。 因此，矩阵B将类似于：

New_B = [0.1 0.15 0.2 0.25 0.3; ...
         0.1 0.15 0.2 0.25 0.3; ...
         0.1 0.15 0.2 0.25 0.3];

矩阵T的预期结果将是：

T = [1 2 3 4 5; ...
     2 3 4 5 6; ...
     1 1.5 2 2.5 3];

我试过了：

New_T = interp2(T, B, New_B);

但这行不通。 我该如何进行插值？

Answer 1

您将必须使用for循环，因为interp1仅适用于单行。

T=[ 1 3 5; 2 4 6; 1 2 3];
B=[0.1 0.2 0.3; 0.1 0.2 0.3; 0.1 0.2 0.3];


B_new = [0.1 0.15 0.2 0.25 0.3; 0.1 0.15 0.2 0.25 0.3; 0.1 0.15 0.2 0.25 0.3];

T_new = cell(2,1);

for k=1:size(B,1)

    T_new{k} =  interp1(B(1,:),T(k,:),B_new(k,:)); 

end

T_new = cell2mat(T_new)

输出：

T_new =

    1.0000    2.0000    3.0000    4.0000    5.0000
    2.0000    3.0000    4.0000    5.0000    6.0000
    1.0000    1.5000    2.0000    2.5000    3.0000

Answer 2

您可以将interp1应用于矩阵，在这种情况下，它将沿每一列运行。 由于要沿每一行插值，因此必须转置输入和输出。 您还只需要B和New_B每一行：

New_T = interp1(B(1, :).', T.', New_B(1, :).').';

New_T =

    1.0000    2.0000    3.0000    4.0000    5.0000
    2.0000    3.0000    4.0000    5.0000    6.0000
    1.0000    1.5000    2.0000    2.5000    3.0000

如果您感到好奇，则必须指定其他行网格点才能使用interp2 ：

New_T = interp2(B(1, :), (1:size(B, 1)).', T, New_B(1, :), (1:size(New_B, 1)).');