朱莉娅中的 3D 网格数据插值

Question

我正在努力将一些MATLAB代码转换为Julia。我有一些 3D 网格数据(温度随二维和随时间变化(，并希望从 (x，y，t( 网格更改为更松散的 (xi，yi，ti( 网格。在 MATLAB 中，这将是一个简单的interp(x,y,t,T,xi,yi,ti).

我尝试使用Interpolations，Dierckx，但两者似乎都只适用于2D网格数据。我做错了什么吗？我对朱莉娅编程很陌生...

我已经在考虑通过具有一些NumPy/SciPy功能的PyCall解决问题的可能性。

谢谢！

Answer 1

是什么让你相信Interpolations.jl只适用于二维数据？

julia> a = rand(1:100, 10, 10, 10);
julia> using Interpolations
julia> itp = interpolate(a, BSpline(Linear()));
julia> v = itp(1.4, 2.3, 3.7)
55.24