質問: interp2 による内挿と、griddataにおける内挿との結果が異なります 正しい動作でしょうか?
回答: 補間アルゴリズムが異なるため、計算結果は異なります。 特に、griddata は、三角形分割に基づく内挿補間をしています。 サンプル例を示しています。 理解しやすい線形補間の例です。 rng(100) [x,y] = meshgrid(-3:3); v = reshape(randperm(49),[7 7]); [xq,yq] = meshgrid(-3:0.5:3); vqi = interp2(x,y,v,xq,yq,'linear'); vqg = griddata(x,y,v,xq,yq,'linear'); figure diffdata = vqg - vqi; mesh(xq,yq,diffdata) % 確かに下の例で、中央の数値が異なっています。 vqi(1:3,1:3) vqg(1:3,1:3) % interp2においては、メッシュの4頂点の平均 isequal((vqi(1,1) + vqi(3,1) + vqi(1,3) + vqi(3,3))/4,vqi(2,2)) % griddataは、斜めに対角線を引くので、(3,1) と(1,3)と平均 isequal((vqi(3,1) + vqi(1,3))/2,vqg(2,2))