计算和绘制两个点云之间的差异

Question

我有两个不同维度的点云（XYZ 坐标），并希望能够计算它们之间的差异（结果为 XYZ 数组，Z 是它们之间的距离），然后将两者绘制为表面连同作为不同颜色的差异。

这个问题似乎很相关，但不是我想要的：Subtract two trisurf plots from one another

这里是一些示例数据和代码。我不知道如何将一个数据集插入另一个数据集。

%Point Cloud 1
X1 = randn(100,1);
Y1 = randn(100,1);
Z1 =(exp(-X1.^2-Y1.^2));
% Point Cloud 2
X2 = randn(107,1);
Y2 = randn(107,1);
Z2 = (exp(-X2.^2-Y2.^2));

tri1 = delaunay(X1, Y1);
tri2 = delaunay(X2, Y2);

trisurf(tri1, X1, Y1, Z1, 1)
hold on
trisurf(tri2, X2, Y2, Z2, 100)
hold off

我上面提到的问题指向这里：How Do I Generate a 3-D Surface From Isolines? 但我是在 Matlab 中插值 3d 数据的新手，似乎无法弄清楚。任何帮助，将不胜感激。谢谢。

Answer 1

我认为这会有所帮助，实际上它只是您没有的插值方法。

%Point Cloud 1
X1 = randn(100,1);
Y1 = randn(100,1);
Z1 =(exp(-X1.^2-Y1.^2));
% Point Cloud 2
X2 = randn(107,1);
Y2 = randn(107,1);
Z2 = (exp(-X2.^2-Y2.^2));

% Mesh for interpolation
x=linspace(min([X1;X2]),max([X1;X2]),40);
y=linspace(min([Y1;Y2]),max([Y1;Y2]),40);
[X,Y]=meshgrid(x,y);

% Calculate interpolants
V1=TriScatteredInterp(X1,Y1,Z1);
V2=TriScatteredInterp(X2,Y2,Z2);
F1=V1(X,Y);
F2=V2(X,Y);

% Plot to check results!
figure(1)
scatter3(X1,Y1,Z1)
hold on
mesh(X,Y,F1)
hold off
figure(2)
scatter3(X2,Y2,Z2)
hold on
mesh(X,Y,F2)
hold off

你也可以使用griddata:

% Calculate interpolants
F1=griddata(X1,Y1,Z1,X,Y);
F2=griddata(X2,Y2,Z2,X,Y);