Matlab的4D分散数据插值

Question

我们有一个195x4（双精度）矩阵：A=[X Y Z temp]，当用scatter3(A(:,1), A(:,2), A(:,3),30, A(:,4), 'filled' )绘制时，会得到如下所示：

现在，我们要生成一个带有temp=A(:,4)矢量数据插值的“立方体”。 到目前为止，我们已经尝试过interp3

% Base Grid
[Xm Ym Zm] = meshgrid(A(:,1), A(:,2), A(:,3)); 
% Grid Refinement
[Xq,Yq,Zq] = meshgrid(xmin:dx:xmax, ymin:dy:ymax, zmin:dz:zmax);
Aq = interp3(Xm,Ym,Zm,A(:,4),Xq,Yq,Zq);

返回以下错误：

Error using griddedInterpolant
The number of input coordinate arrays does not equal the number of dimensions
(NDIMS) of these arrays.

Error in interp3 (line 144)
            F = griddedInterpolant(X, Y, Z, V, method,extrap);

Error in PDGEM_MT (line 112)
Aq = interp3(Xm,Ym,Zm,A(:,4),Xq,Yq,Zq);

因此，我认为这可能是错误的实现方式和/或对问题的错误解释。 如何生成用A(:,4)的体积插值着色的那个空间的“立方体”？

谢谢。

Answer 1

您有一个分散的数据集。

interp3仅在您的数据点为网格格式read this时有效。此功能的简短描述是Interpolation for 3-D gridded data in meshgrid format

您可以使用griddata来处理分散的数据read this。简短描述为Interpolate 2-D or 3-D scattered data。

示例：

X = rand(100,1);
Y = rand(100,1);
Z = rand(100,1);
C = rand(100,1);

figure
scatter3(X, Y, Z,30, C, 'filled' )

[Xm, Ym, Zm] = meshgrid(min(X):.01:max(X), min(Y):.01:max(Y), min(Z):.01:max(Z)); 

Cm = griddata(X,Y,Z,C,Xm,Ym,Zm);

figure
scatter3(Xm(:), Ym(:), Zm(:), 30, Cm(:), 'filled' )

数据点：

插值：