5

点群のすべての点にまたがる3Dサーフェスのプロットを作成したいと思います[X,Y,Z]。たとえば、これは私の点群の散布図です。

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一致しないためです。Xbtri

正しい補間された表面色を確実に取得するにはどうすればよいですか？

4

1 に答える 1

4

freeBoundaryを呼び出して、プロットされた三角形分割のポイント数を変更しました。サーフェスポイントのみが残り、内側のポイントはサーフェスに属しません。したがって、これらのポイントに対応するC値を抽出する必要があります。'intersect（...、' rows'）'を使用して、サーフェスポイント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');
于 2012-09-13T12:47:01.300 に答える