Breast Density in Mammography Dicom Images
4 次查看(过去 30 天)
显示 更早的评论
Is there a way to extract the breast density of a mammography through Matlab code?
1 个评论
采纳的回答
Said Pertuz
2019-11-14
I hope is not too late an answer. Please take a look at the following tool:https://www.mathworks.com/matlabcentral/fileexchange/73360-breast-density-segmentation.
Beware that this implementation has been tested on digital mammograms (such as those from the INbreast dataset) and has not been tested on digitized mammograms (e.g. MIAS).
更多回答(1 个)
Doaa
2024-7-27
2 个评论
Doaa
2024-7-27
编辑:Walter Roberson
2024-7-27
% Load the mammogram image
img = imread('mammogram.jpg');
% Display the original image
figure;
imshow(img);
title('Original Mammogram Image');
% Convert the image to grayscale if it is not already
if size(img, 3) == 3
img = rgb2gray(img);
end
% Display the grayscale image
figure;
imshow(img);
title('Grayscale Mammogram Image');
% Apply median filter to reduce noise
filtered_img = medfilt2(img);
% Display the filtered image
figure;
imshow(filtered_img);
title('Filtered Mammogram Image');
% Binarize the image using a threshold
level = graythresh(filtered_img);
bw = imbinarize(filtered_img, level);
% Display the binary image
figure;
imshow(bw);
title('Binary Mammogram Image');
% Calculate the breast density
density = sum(bw(:)) / numel(bw) * 100;
% Display the breast density
fprintf('Breast density: %.2f%%\n', density);
% Classify the breast density
if density < 25
density_type = 'Fatty';
elseif density < 50
density_type = 'Scattered';
elseif density < 75
density_type = 'Heterogeneously dense';
else
density_type = 'Extremely dense';
end
fprintf('Breast density type: %s\n', density_type);
Doaa
2024-7-27
编辑:Walter Roberson
2024-7-27
% Load the mammogram image
img = imread('mammogram.jpg');
% Display the original image
figure;
imshow(img);
title('Original Mammogram Image');
% Convert the image to grayscale if it is not already
if size(img, 3) == 3
img = rgb2gray(img);
end
% Display the grayscale image
figure;
imshow(img);
title('Grayscale Mammogram Image');
% Apply median filter to reduce noise
filtered_img = medfilt2(img);
% Display the filtered image
figure;
imshow(filtered_img);
title('Filtered Mammogram Image');
% Binarize the image using a threshold
level = graythresh(filtered_img);
bw = imbinarize(filtered_img, level);
% Display the binary image
figure;
imshow(bw);
title('Binary Mammogram Image');
% Calculate the breast density
density = sum(bw(:)) / numel(bw) * 100;
% Display the breast density
fprintf('Breast density: %.2f%%\n', density);
% Classify the breast density
if density < 25
density_type = 'Fatty';
elseif density < 50
density_type = 'Scattered';
elseif density < 75
density_type = 'Heterogeneously dense';
else
density_type = 'Extremely dense';
end
fprintf('Breast density type: %s\n', density_type);
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!