xyz座標のベクトル​​データ(散布データ​)​から3次元メッシ​ュプ​ロットや表面積​の出し方

29 次查看（过去 30 天）
和大 坂田
和大 坂田 2020-12-20
编辑： Atsushi Ueno 2021-10-17
Excelやテキストファイルにx, y,z のデータがそれぞれベクトルデータ（散布データ）として定義されています。
このデータを使用して3次元メッシュプロットを描き表面積のを求める方法を教えてください。
  2 个评论
源樹 上林
源樹 上林 2020-12-22
delaunay と polyarea あたりで、がんばればできる…ような気がする。
源樹 上林
源樹 上林 2021-10-15
convexHull なんて便利なのがあったのか。例題はyのデータがゼロだから polyarea で面積出せそうですね。
もしyにも凹凸がある場合は、面積計算は cross を利用できます。外積の半分で三角形の面積になるので delaunay で作った三角形の面積を全部足し合わせれば表面積が出ます。

请先登录，再进行评论。

请先登录，再回答此问题。

回答（2 个）

Kenta
Kenta 2020-12-27
clear;clc;close all
ptCloud = pcread('teapot.ply');
pcshow(ptCloud);
ptCloud=pcdownsample(ptCloud,"random",0.1);
P=double(ptCloud.Location);
% k = boundary(P);
% trisurf(k,P(:,1),P(:,2),P(:,3),'Facecolor','red','FaceAlpha',0.1)
DT = delaunayTriangulation(P(:,1),P(:,2),P(:,3));
tetramesh(DT,'FaceAlpha',0.3);
[K,v] = convexHull(DT);
trisurf(K,DT.Points(:,1),DT.Points(:,2),DT.Points(:,3))
こんにちは、源樹 上林さまのおっしゃるように、delauney三角形を作成すれば、それらしいものはできます。
試したことはないのですが、以下のようなツールもあるようです。いろいろと試してみていい感じのメッシュができたら面積計算を試みてはどうでしょうか
http://iso2mesh.sourceforge.net/cgi-bin/index.cgi?Doc/Installation
Atsushi Ueno
Atsushi Ueno 2021-10-16
编辑：Atsushi Ueno 2021-10-17
Ran in:
2 次元および 3 次元の点で構成される多角形および多面体 - MATLAB - MathWorks 日本
3 次元アルファ形状の表面積 - MATLAB surfaceArea - MathWorks 日本
回答に際し今知った関数を使ってみたのですが、理解が浅いためか使い勝手が悪いです。
  • アルファ形状の作成時に指定する「アルファ半径」が謎
  • 何故か面積の値が倍になります。表と裏があるって事でしょうか (追記)見た目から明らかに「アルファ形状」は三角メッシュです。三角形1個を構成する2辺の外積を求めているのであれば、2で割って三角形の面積にしなければならないという事なのでしょうか？surfaceAreaの仕様がよくわかりません
Excelやテキストファイルに定義されたベクトルデータ（散布データ）はこういう事ですよね
n = 10; % xyz:ベクトルデータ(散布データ)の例を作成
xyz(:, 1) = repmat((0:n)', [n+1 1]);
xyz(:, 2) = repelem((0:n)', n+1);
xyz(:, 3) = rand((n+1)^2, 1);
plot(xyz) % xyz:ベクトルデータ(散布データ)の例
shp = alphaShape(xyz, 1.9); % アルファ形状を作成
totalsurfarea = surfaceArea(shp) % アルファ形状の表面積を計算
totalsurfarea = 221.9298
plot(shp);
xyz(:, 3) = rand((n+1)^2, 1).* eps; % xyz:ベクトルデータ(散布データ)をまっすぐに伸ばす
shp = alphaShape(xyz, 1.0); % アルファ形状を作成
totalsurfarea = surfaceArea(shp) % アルファ形状の表面積を計算
totalsurfarea = 200
plot(shp);

类别

在 Help Center 和 File Exchange 中查找有关 表面プロットとメッシュ プロット 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!