Matlab Delaunay点云的三角剖分 - 颜色矩阵

Question

我想创建一个横跨[X,Y,Z]点云所有点的3D表面图。例如，这是我的点云的散点图：

scatter3(X,Y,Z,5,C)

正如您所看到的，每个数据点的强度值都为C。

我现在进行三角测量

dt      = DelaunayTri(X,Y,Z); 
[tri Xb]= freeBoundary(dt); 

我得到了三角形表面

figure 
trisurf(tri,Xb(:,1),Xb(:,2),Xb(:,3), 'FaceColor', 'cyan', 'faceAlpha', 0.8);

然而，当我尝试使用

设置表面的颜色时

trisurf(tri,Xb(:,1),Xb(:,2),Xb(:,3),C,'EdgeAlpha',0,'FaceColor','interp')

我收到错误消息：“警告：颜色数据未设置为插值着色”，这是因为C的大小与Xb或{{1}不匹配}。

如何确保获得正确的插值表面颜色？

Answer 1

您通过调用freeBoundary更改了绘制三角剖分中的点数：仅剩下曲面点，内部点不属于曲面。因此，您必须提取与这些点对应的C值。您可以使用＆＃39; cross（...，＆＃39; rows＆＃39;）＆＃39;将表面点Xb映射到原始点集XYZ上。基于此映射，您可以从C中提取所需的值。下面的代码执行此操作。

clear all;

XYZ = rand(100,3);
X=XYZ(:,1);
Y=XYZ(:,2);
Z=XYZ(:,3);
C=rand(size(X));

scatter3(X, Y, Z, 5,C);

dt = DelaunayTri(X, Y, Z);
[tri Xb]=freeBoundary(dt);

% map Xb onto XYZ
[~,IA,IB]=intersect(XYZ, Xb, 'rows');

% extract the needed colors using the IA map
Cn      = C(IA);

% permute the surface triangulation points using IB map
Xbn     = Xb(IB,:);

% map the point numbers used in triangle definitions
% NOTE: for that you need inverse map
iIB(IB) = 1:length(IB);
trin    = iIB(tri);

trisurf(trin,Xbn(:,1),Xbn(:,2),Xbn(:,3),Cn,'EdgeAlpha',0,'FaceColor','interp');