3Dデータ中の球体カウント

医用画像のテクスチャ解析のために、３Ｄデータ中の球をカウントする方法を知りたいです。

対象の３Ｄデータは添付画像をイメージしています。

[michio]

添付の画像だと、求めたい球体の数は３でしょうか？

[VAMOU]

回答ありがとうございます。 添付画像の場合、求める球体の数は３になります。

[michio]

３つの球体はそれぞれが小さな球体で構成されているように見えたので質問いたしました。解析対象の元データとしては３次元の点群の座標データ（x,y,z) ということでよろしいですか？

[VAMOU]

仰るように３つの球体はそれぞれ小さな球体で構成されたかたまりとして、元データは3次元点群座標データ(x,y,z)です。 MATLAB内のヘルプにありました例「球内の乱数」に輝度値をランダムに与えて作成しております。

[michio]

詳細ありがとうございます。

Answer 1

[michio]

もし画像のように綺麗に分離している球体の形をしている点群のみであれば MATLAB の alphaShape オブジェクトでうまくいくかもしれません。いろんな形の塊が存在しその中から球体を見つける、という場合には別の方法を考える必要があります。

[x1, y1, z1] = sphere(24);
x1 = x1(:);
y1 = y1(:);
z1 = z1(:);
x2 = x1+5;
P = [x1 y1 z1; x2 y1 z1];
P = unique(P,'rows');
plot3(P(:,1),P(:,2),P(:,3),'.')
axis equal
grid on
shp = alphaShape(x,y,z);
numRegions(shp)

alphaShapeの詳細はこちら： https://jp.mathworks.com/help/matlab/ref/alphashape.alphashape.html

[VAMOU]

michioさん

無事数えることができました。ありがとうございます。

[Takuji Fukumoto]

Volume画像で数の測定だけでなく、形状の解析もされる場合には regionprops3 も便利ですので、よろしければご利用ください。 https://www.mathworks.com/help/releases/R2017b/images/ref/regionprops3.html

[VAMOU]

Fukumotoさん

参考に致します。ありがとうございます。