matlab - Matlab / Octave中X，Y，Z，V有序数据的等值面

Question

我有一组 3d 数据点，每个数据点与某个数量的值 V 相关联。数据组织在有序列 XYZV 中，其中空间坐标分布在网格上。我在每个方向都有 21 个点，所以几列的长度是 21x21x21=9261。如何将数据转换为可以isosurface在 Octave/Matlab 中使用的网格？我的代码如下

a=load("data.txt");
X=reshape(a(:,1), 21,21,21);
Y=reshape(a(:,2), 21,21,21);
Z=reshape(a(:,3), 21,21,21);
V=reshape(a(:,2), 21,21,21);
fv=isosurface (X,Y,Z,V,0.9);
patch(fv)

但结果没有意义（我得到两个平面位于 x=0.9 和 1。）。数据可以在这里下载。

Answer 1

这是一种为您的数据创建适当网格网格的方法：

a = load('data.txt');
[X, Y, Z] = meshgrid(unique(a(:, 1)), unique(a(:, 2)), unique(a(:, 3)));
V = zeros(21, 21, 21);
for i = 1:numel(V)
  idx = find(a(:, 1) == X(i) & a(:, 2) == Y(i) & a(:, 3) == Z(i));
  V(i) = a(idx, 4);
end
fv = isosurface (X,Y,Z,V,0.9);
p = patch(fv);
set(p,'FaceColor','red','EdgeColor','none');
camlight;
lighting gouraud;
xlabel('x');
ylabel('y');
zlabel('z');

如果你有更大的数据而且速度太慢，我可能会想出一种重塑原始数据的方法，以避免上面使用的 for 循环。