MATLAB Answers

Segmenting bitmap micro CT images using MATLAB's image processing capabilities.

9 views (last 30 days)
Ross Remesz
Ross Remesz on 10 Jul 2020
Commented: Ross Remesz on 11 Jul 2020
Hello, I am interested in getting some input about using MATLAB's image processing app. I want to obtain a binary image of my area of interest so I can perform some basic analysis (calculate area, path length, etc.).
I want to know if what I describe below is possible in MATLAB or not. I will show the process that I have gone through so far and any input would be greatly appreciated.
Progress Thus Far:
% This will be used for image segmentation of cranial sutures.
% Reading raw images
Raw = imread('3351_Coronal_CRotated0964.bmp');
figure
imshow(Raw)
title('Raw Image')
% Adjusting raw images so they can be turned into binary
Raw_Adjusted = imadjust(Raw);
figure
imshow(Raw_Adjusted)
title('Adjusted Image')
% Converting to black and white images
BW_Original = imbinarize(Raw_Adjusted,'adaptive','Sensitivity',0.6);
figure
imshow(BW_Original)
title('Binary Image')
% Inverting the black and white image
BW_Inverse = imcomplement(BW_Original);
figure
imshow(BW_Inverse)
title('Inverted Binary Image')
% Cropping the image
BW_Cropped = imcrop(BW_Inverse,[900,1000,825,500]);
figure
imshow(BW_Cropped)
title('Cropped Image')
% Segmenting
Segmented = bwpropfilt(BW_Cropped,'Area',1);
figure
imshow(Segmented)
title('Semi-Isolated Cropped Image')
Next Steps/Moving Forward
I want the final product to include only the squiggly line running horizontal, outlines below in red, and get rid of the sections that are outline in blue. (i.e. turn them black)
Any guidance would be greatly appreciated, thank you for reading my question.

Accepted Answer

Image Analyst
Image Analyst on 10 Jul 2020
If Segmented is your binary image, obtain the largest blob, which will be the one you outlined in red, using bwareafilt():
largestBlob = bwareafilt(Segmented, 1);
  3 Comments
Ross Remesz
Ross Remesz on 11 Jul 2020
Perfect, this is exactly what I was looking for.
I just implemented it into my script and it works like a dream, thank you very much for your time.

Sign in to comment.

More Answers (0)

Community Treasure Hunt

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

Start Hunting!