我有一个形状为 f[nz,ny,nx] 的某个数量 f 的 3D 数据集，其中 nz,ny,nx 是对应的笛卡尔网格 z[nz]、y[ny]、x[nx] 的索引。在我的特殊情况下，nz=ny=nz=256。这个笛卡尔盒子的中心是（0,0,0）。在空间的同一点，我想在笛卡尔盒子内放置一个半径为 r=1 的球体。现在，我想在球体表面对 f 进行插值，最终得到 f(phi,theta,r=1)，其中 phi 是经度，theta 是纬度。
我已经看到了许多不同的方法来解决与此类似的问题（主要是在 Python 和 Matlab 中），但还没有能够做到。有人做过这个或知道如何做吗？谢谢！