interp函数在插值前后网格对应。
例如:
x = 0.5:.5:5;
y = 0.5:.5:3;
temps = [82 81 80 82 84 83 81 79 83 87;76 69 76 79 63 61 65 77 80 81;84 87 84 82 85 86 79 82 83 84;...
76 75 78 87 79 80 81 85 79 79; 79 84 83 82 80 80 84 87 80 79;79 84 83 82 80 80 84 87 80 79];
% mesh(x,y,temps);
x1 =0.8:4.8;
y1 =1:3;
[Y,X]=meshgrid(x,y);
[Y1,X1]=meshgrid(x1,y1);
zzi = griddata(Y,X,temps,Y1,X1);
zzi = interp2(Y,X,temps,Y1,X1);
目前,在线性插值(linear)条件下,griddata和interp2效果一样。二者具体差别这里不做讨论。
temps=[82 81 80 82 84 83 81 79 83 87
76 69 76 79 63 61 65 77 80 81
84 87 84 82 85 86 79 82 83 84
76 75 78 87 79 80 81 85 79 79
79 84 83 82 80 80 84 87 80 79
79 84 83 82 80 80 84 87 80 79];
插值前网格为6X10,插值后为3X5,但并不是2X2的网格求平均,而是对应插值位置进行线性运算,例如:
插值后矩阵的第一个网格(1,1)实际上是 x1=0.8, y1=1
对应来看位于插值前矩阵网格(1,2)和(2,2)之间,即在69和76之间进行线性插值,76 - 69在0.5~1之间分为5个间隔,每个间隔相差-1.4,对应到0.8的位置即为 76-1.4*3=71.8
如此可以推知插值后矩阵的结果。
插值后
zzi=[71.8 77.8 61.8 72.2 80.6
75.4 83.4 79.6 83.4 79
82 82.4 80 85.8 79.4]