4次元分布における複数重心点座標の取得について

Question

ryo tanaka 2019-10-17

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/485904-4

评论： Kenta 2019-10-21

采纳的回答： Kenta

在 MATLAB Online 中打开

添付図のように4次元分布を作成しました。

4次元分布上に5個の丸い分布があります。

この5個の分布それぞれの重心点を取得したいです。

分布上から1点の重心座標は下記のコードで取得できると思いますが、1つの分布上から任意の個数だけ重心点を取得する方法がありましたら

教えて頂きたいです。

v = squeeze(v); 
BW = true(size(v));
W = regionprops3(BW,v,'WeightedCentroid');
WC = W{1,{'WeightedCentroid'}};

9 个评论
显示 7更早的评论隐藏 7更早的评论

ryo tanaka 2019-10-18

在 MATLAB Online 中打开

ご返信ありがとうございます。

添付の仕方を教えて頂きありがとうございます。

コードは下記に記載しました。

しかし、mat fileの添付は下記の理由でできないようです。

初めてこのような文がでてきたので、よくわかりませんが、明日以降にファイルの添付をすることになるかと思います。申し訳ありません。

また、見分ける基準は各塊のピーク値(付近)に重心点が置ければいいかなと思います。

（あくまで、塊の中心ではなく塊のピーク点です。）

次の理由でこのファイルを添付できません:

アップロードは 1 日 10 件に制限されています。追加ファイルをアップロードする必要がある場合、今すぐ 1 つ以上のファイルを削除するか、または 24 時間が経過してからアップロードしてください。

load X.mat
load Y.mat
load Z.mat
load V.mat
xslice = [-1/2 12];
yslice = 15;
zslice = [519/2 349];
figure
slice(X, Y, Z, v, xslice, yslice, zslice)                    % display the slices
  
shading faceted
ylim([ymin-1 ymax+1])
view(-34,24)
% set(gcf, 'Color', 'none');
cb = colorbar;                                               % create and label the colorbar
cb.Label.String = 'Temperature, C';
% 重心位置の決定
v = squeeze(v); 
BW = true(size(v));
W = regionprops3(BW,v,'WeightedCentroid');
WC = W{1,{'WeightedCentroid'}};

ryo tanaka 2019-10-18

attachmentsをクリックしファイルを選択してもやはりエラーが出てきてしまいます。。

スタッフの方にご連絡したいと思います。

ryo tanaka 2019-10-19

すみません、遅くなりました。

ようやくファイルを添付できましたので、お送りいたします。

请先登录，再进行评论。

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

Answer 1

Kenta 2019-10-19

6
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/485904-4#answer_397143

在 MATLAB Online 中打开

こんにちは。データを送っていただきありがとうございます。

X,Y,Z,ｖとありますが、X-Zは、ｖのボクセル座標を定めるためのグリッドで、それぞれのグリッド内の値が

そのX（やY・Z)座標に対応している、そして、ｖが実際にそのプロットの持つ強度である、ということでしょうか。

コメントにて、コードも添付いただきありがとうございました。申し訳ございませんが、こちらで改めてコードを書きました。

質問者様は、上のようにボクセルを定義して、各重心を求める方針だったと読み取りましたが、

下では、ポイントクラウドで処理しています。下のように書けば重心を求めるまでかなり短く書くことができます。

重心は一般的な３D上でのものを計算しています。本来はこの過程でｖの強度を利用して、重みつきの重心を定義するのですかね？また、4次元、というキーワードがありますが、これを４D画像としてとらえるなら、ｖの値を利用せず重心を求めているのでもう少し変更する必要があるとは思います。適宜変更いただければと思います。

詳細は適宜コード内のコメントをご参照ください。よろしくお願いいたします。

clear;clc;close all
% data loading
load V.mat
% find out the points with non-zero intensity
idx=find(v~=0);
% extract the XYZ coordinate
[subx, suby, subz]=ind2sub(size(v),idx);
% acquire the intensity, it was not used, though
int = v(idx);
% convert the xyz into pointcloud variable
ptCloud=pointCloud([subx, suby, subz]);
% figure;pcshow(ptCloud,'MarkerSize',50)
% closest pair of points whose euclidean distance is less than "minDistance"
% was assigned into the same cluster
minDistance = 3;
[labels,numClusters] = pcsegdist(ptCloud,minDistance);
% note that each cluster could be segmented based on the distance
figure;pcshow(ptCloud.Location,labels+1,'MarkerSize',50)
colormap([hsv(numClusters);[0 0 0]])
title('Segmented Cloud Clusters');hold on
% the centroid of each cluster was caclulated and plotted on the 3D space
% note that any weighting  based on the intensity was not applied 
s = countcats(categorical(labels));
[~, c_label] = sort(s,'descend');
for i=1:5
    roi_idx=find(labels==c_label(i));
    roi = ptCloud.Location(roi_idx,:);
    xyz=mean(roi,1);
    plot3(xyz(1),xyz(2),xyz(3),'o','MarkerSize',20,'MarkerFaceColor','w');hold on 
end

6 个评论
显示 4更早的评论隐藏 4更早的评论

ryo tanaka 2019-10-20

回答ありがとうございます。

＞X,Y,Z,ｖとありますが、X-Zは、ｖのボクセル座標を定めるためのグリッドで、それぞれグ＞リッド内の値がそのX（やY・Z)座標に対応している、そして、ｖが実際にそのプロットの＞持つ強度である、ということでしょうか。

その通りですね。Vがプロットのもつ強度になります。

＞質問者様は、上のようにボクセルを定義して、各重心を求める方針だったと読み取りました＞が、下では、ポイントクラウドで処理しています。下のように書けば重心を求めるまでかな＞り短く書くことができます。

強度をもつ3次元座標（X,Y,Z）をボクセル上に投影させて、座標点にボリュームをもたせて、そのボクセル上から重心点を求めようとしていました。

ポイントクラウドというのは、ボクセルではなく、各座標点をそのまま点として格納する方法でしょうか？

＞重心は一般的な３D上でのものを計算しています。本来はこの過程でｖの強度を利用して、＞重みつきの重心を定義するのですかね？

そうですね。強度値から重心を求めようとしていました。

＞また、4次元、というキーワードがありますが、これを４D画像としてとらえるなら、ｖの＞値を利用せず重心を求めているのでもう少し変更する必要があるとは思います。適宜変更い＞ただければと思います。

自分がやってる方法だとｖの値は利用されていなかったのですね、勉強不足でした。

参照コードの添付ありがとうございます。図に示してある白い丸が重心点でしょうか？

この5つの丸の座標はどのように出すのでしょうか？

Kenta 2019-10-20

こんばんは、ご丁寧に返信ありがとうございます。

重心点を得られたとのことでよかったです。

>>できれば強度の値から重心位置を取得したいと思いますので、最後の重心位置を求めるコー>>ド部分を変更しないといけませんね

はい、そうですね。またいろいろと試したうえで不明点があれば新たに質問いただければと思います。

>>labels+1と[0 0 0]がどういった意味をしめしているのかわからない

詳しく目を通していただきありがとうございます。特に両者とも特に意味はありません。

labelsの値が１以上になってほしかったので、こうしましたが、labelsはすべて１以上の値が入っています。特に意図せず書いていました。

[0 0 0]のほうですが、[0 0 0]だけ独立してあるのではなくて、

hsv(numclusters)と縦方向に積み重ねています。; マークが縦に積むということを意味しています。つまり、カラーマップの終点を[0 0 0]（黒）にしています。こうするとFigure中の点の色が一部変更されます。特に大きな意味はないです。下で、ほかのいろのパターンも探せます。

https://jp.mathworks.com/help/matlab/ref/colormap.html#buc3wsn-1-map

ryo tanaka 2019-10-21

ご返信ありがとうございます。

＞はい、そうですね。またいろいろと試したうえで不明点があれば新たに質問いただければ

＞と思います。

ありがとうございます。とても助かります。

質問に回答していただきありがとうございます。

2つとも特に気にすることはないようで安心しました。

また、何か不明な点がありましたら、ご連絡させて頂きます。

Kenta 2019-10-21

はい、承知いたしました。解決してよかったです。

请先登录，再进行评论。

4次元分布における複数重心点座標の取得について

9 个评论
显示 7更早的评论隐藏 7更早的评论

采纳的回答

6 个评论
显示 4更早的评论隐藏 4更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

4次元分布における複数重心点座標の取得について

9 个评论 显示 7更早的评论隐藏 7更早的评论

采纳的回答

6 个评论 显示 4更早的评论隐藏 4更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

9 个评论
显示 7更早的评论隐藏 7更早的评论

6 个评论
显示 4更早的评论隐藏 4更早的评论