How to combine all the segmented images in one figure

5 ビュー (過去 30 日間)
mohd akmal masud
mohd akmal masud 2021 年 10 月 21 日
回答済み: yanqi liu 2021 年 10 月 26 日
Hi all, I have code below.
Anyone know how to combine all the segmented images in one figure?
clc
clear all
[spect map]=dicomread('i-131-phantom validation 04112019 256x256');
info = dicominfo('i-131-phantom validation 04112019 256x256');
%gp=info.SliceThickness;
spect=(squeeze(spect));%smooth3
aa=size(spect);aa=aa(3);
% NI KALAU VOLUME1 en amin prepapre 10:1
% BCKG 468 MAX 381
seedR1 = 111; seedC1 = 150; seedP1 = 88;
W1 = graydiffweight(spect, seedC1, seedR1, seedP1 , 'GrayDifferenceCutoff', 1000000);
thresh1 = 0.004349;
[BW1, D1] = imsegfmm(W1, seedC1, seedR1, seedP1, thresh1);
% figure, imshow3D(BW1)
T1 = regionprops('table', BW1,'Area','Centroid')
% NI KALAU VOLUME2 en amin prepapre 10:1
% BCKG 468 MAX 221
seedR2 = 135; seedC2 = 149; seedP2 = 88;
W2 = graydiffweight(spect, seedC2, seedR2, seedP2 , 'GrayDifferenceCutoff', 1000000);
thresh2 = 0.001399;
[BW2, D2] = imsegfmm(W2, seedC2, seedR2, seedP2, thresh2);
% figure, imshow3D(BW2)
T2 = regionprops('table', BW2,'Area','Centroid')
% NI KALAU VOLUME3 en amin prepapre 10:1
% BCKG 468 MAX 143
seedR3 = 147; seedC3 = 129; seedP3 = 88;
W3 = graydiffweight(spect, seedC3, seedR3, seedP3 , 'GrayDifferenceCutoff', 1000000);
thresh3 = 0.000319;
[BW3, D3] = imsegfmm(W3, seedC3, seedR3, seedP3, thresh3);
% figure, imshow3D(BW3)
T3 = regionprops('table', BW3,'Area','Centroid')
% NI KALAU VOLUME4 en amin prepapre 10:1
% BCKG 468 MAX 109
seedR4 = 135; seedC4 = 108; seedP4 = 88;
W4 = graydiffweight(spect, seedC4, seedR4, seedP4 , 'GrayDifferenceCutoff', 1000000);
thresh4 = 0.000091;
[BW4, D4] = imsegfmm(W4, seedC4, seedR4, seedP4, thresh4);
% figure, imshow3D(BW4)
T4 = regionprops('table', BW4,'Area','Centroid')
imshow3D(BW1, BW2, BW3,BW4)
If I did seperately, can. But I do not seperately.
This one is BW1
This one is BW2
This one is BW3
This one is BW4
Should be all the sphere in one figure as example below

採用された回答

yanqi liu
yanqi liu 2021 年 10 月 26 日
use
BW = logical(BW1+BW2+BW3+BW4);

その他の回答 (0 件)

Community Treasure Hunt

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

Start Hunting!

Translated by