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

Question

我有一组3d数据点，每个数据点都与某个特定数量的值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循环。