Matlab的4D分散数据插值

Question

当使用scatter3(A(:,1), A(:,2), A(:,3),30, A(:,4), 'filled' )绘制时，我们有一个195x4（双）矩阵： A=[XYZ 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如果您的数据点在meshgrid格式并只能工作阅读 。 此功能的简短描述是Interpolation for 3-D gridded data in meshgrid format

相反，您可以使用适用于分散数据的griddata来阅读此内容 。 简短描述是内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' )

数据点：

插值。