来自X，Y，Z，V的Isosurface在Matlab / Octave中排序数据

Question

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

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）。数据可以下载here。

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循环。