MATLAB（或Python或C）中的3D数据插值和导数
我有三个三维矩阵 X、Y 和 Z，它们定义了一个与它们大小相同的矩阵 V，这个矩阵在某些三维空间中。这个矩阵的点是均匀分布的。现在，我想进行插值计算，并且还想计算 V 的空间偏导数，也就是说，对于每一个像素，我想知道 V 在 x、y 和 z 方向上是如何变化的。我听说用样条插值（splines）来进行插值和计算导数效果很好。比如，我之前用过 splinefit 和 ppdiff（http://www.mathworks.com/matlabcentral/fileexchange/13812-splinefit）。
我该如何使用样条插值来处理我手头的数据集呢？有没有现成的代码，最好是 MATLAB 的（Python 和 C 也可以）来进行这些计算？
假设我只想在 X、Y 和 Z 定义的采样位置上得到导数，我是否可以对每个维度做一维样条近似，然后以这种方式计算偏导数？也许这个问题更适合数学交流。虽然可能需要一些时间，但这样应该是可行的，对吧？
谢谢你的帮助！