Takuya Miyashita

内挿 メモ †

1次元 †

1次元の場合は Interpolation パッケージでOK．

using Interpolations

itp = interpolate((xorg, ), Vorg, Gridded(Linear()))
Vint = itp(xint)

Interpolations.interpolate の第一引数はタプルになるらしい．

2次元 †

不均一な点群から規則的な meshgrid へ内挿する方法を諸々試そうと思ってやってみたけどどれも微妙で，結局は scipy を呼ぶのが個人的に一番楽という結論に至った．

using PyCall
scipyinterpolate = PyCall.pyimport("scipy.interpolate")

Z_grid = scipyinterpolate.griddata(hcat(vec(LON),vec(LAT)), vec(Z), (lon_grid,lat_grid), method="cubic")

ここで LON, LAT, Z は元データ，lon_grid, lat_grid はデータを内挿したい座標データ．