bwlabel doesn't work as expected and what would be ideal method to segment disconnected regions like this example
1 visualización (últimos 30 días)
Mostrar comentarios más antiguos
Gopichandh Danala
el 10 de Feb. de 2017
Editada: Omkar Patil
el 16 de Feb. de 2017
Hi,
This is a segmented image of a brain from thresholding segmentation, Now I want to get A, B, and C regions, which are all brain regions.
when I apply some kind of operation. For now, i am using bwlabel to find the blob areas and correspondingly extract bigger or circular blobs of some kind and get blobs A, B,and C. but the problem is when I do bwlabel I don't get A, B blobs.
Here is the code and attached images before and after bwlabel:
temp = before_dicomFiles{1,15};
%figure, imshow(temp, [LOW1 HIGH1]);
temp = double(temp);
newImg = temp; % this is the actual brain image used for thresholding
newImg(newImg <0 ) = 0; % set all negative pixels to zeros
newImg(newImg >1080) = 0; % remove anything over this range
figure, imshow(newImg, [LOW1 HIGH1]); % thresholded image
BW = bwlabel(newImg,8) ;
figure, imshow(BW, []) % Labelled image missing A, Blob areas
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/160553/image.jpeg)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/160555/image.jpeg)
Suggest me some methods to trace back regions A, B and also better segmentation methods to not miss all the brain regions
Thanks, Gopi
0 comentarios
Respuesta aceptada
Image Analyst
el 16 de Feb. de 2017
Try this:
binaryImage = temp > 1080;
% Fill holes (depends on skull having a hole through it.
binaryImage = imfill(binaryImage, 'holes');
% Extract 3 largest blobs.
binaryImage = bwareafilt(binaryImage, 3);
labeledImage = bwlabel(binaryImage);
props = regionprops(labeledImage, temp);
4 comentarios
Image Analyst
el 16 de Feb. de 2017
You might try using bwconvhull to get the convex hull. Then use imerode() to shrink it down. Then use that as a mask to "erase" the outer parts (skull, etc.);
mask = bwconvhull(binaryImage);
mask = imerode(mask, true(15));
binaryImage(~mask) = false;
% Then the previous code I gave.
Omkar Patil
el 16 de Feb. de 2017
Editada: Omkar Patil
el 16 de Feb. de 2017
Thanks, I will try this and give you an update next time I go to lab, Oops i am commenting from my friends computer
Más respuestas (1)
Saurabh Gupta
el 15 de Feb. de 2017
You can use the Image Segmenter App to perform the segmentation, then export it as a MATLAB script to observe the methodology used to implement further tasks in an automated manner.
1 comentario
Ver también
Categorías
Más información sobre Convert Image Type en Help Center y File Exchange.
Productos
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!