Creating binary 3D mask from MR volume

조회 수: 34 (최근 30일)
banikr
banikr 2020년 4월 26일
답변: banikr 2020년 4월 29일
Hello All,
I am trying to create a 3D binary mask for MR head volume. The image has 0 intensity for background and the tissue inside the image has integers. Morphological operations provide me very noisy output.
This is after bwconncomp, thresholding and imerode options.
Are there better ways to get volumetric masks?
  댓글 수: 8
Image Analyst
Image Analyst 2020년 4월 27일
I don't think the noise on the mask is the issue. Why did your classification routine quit classifying stuff below that line other than the blue stuff? For some reason, you have a mask that just quits below that line, and that has nothing to do with little noisy blips on it.
banikr
banikr 2020년 4월 27일
Yeah you are right, the noise has nothing to do with the classification.
The classification comes from CT image which only exists for the segmented section. Now for MR images there is no specific intensity scaling mechanism like CT. For my experiment, I also need the inferior section of the head labelled in MR space.
CT image was registered on MR space.

댓글을 달려면 로그인하십시오.

답변 (2개)

Mrutyunjaya Hiremath
Mrutyunjaya Hiremath 2020년 4월 27일
Hello banikr,
I think this is what you want ... if I am not wrong.
load mri;
% load mriVolume; %3D MRI data with data variable name 'mriVolume'
mriVolume = D;
% Assuming background pixels are 0
bwMask = false(size(mriVolume));
pixIndex = find(mriVolume > 0);
bwMask(pixIndex) = 1;
  댓글 수: 3
Mrutyunjaya Hiremath
Mrutyunjaya Hiremath 2020년 4월 27일
@ banikr,
ok got it, if MRI data is possible to share, please share. love to work on it.
banikr
banikr 2020년 4월 27일

댓글을 달려면 로그인하십시오.


banikr
banikr 2020년 4월 29일
I just wanted to update the finding I have(following mask image).
I used the intensity 10 for thresholding and removed objects with area under 1000.
bi = (mr>10);
% volumeViewer(bi);
stat = regionprops(bi, 'Area', 'PixelIdxList');
for nn=1:length(stat)
s = stat(nn);
if s.Area<1000
continue;
end % remove small objects
bi = zeros(sz);
bi(s.PixelIdxList) = 1;
end
volumeViewer(bi);
I tried to explain the value 10.
J = mr(mr~=0);
>> SkinMusclePrc = prctile(J,[20, 100], 'all')
SkinMusclePrc =
2×1 int16 column vector
10
626
The value 10 is on 20 percentile of the distribution.
let me know if that helps your search also.

Community Treasure Hunt

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

Start Hunting!

Translated by