Breast Density in Mammography Dicom Images

4 次查看(过去 30 天)
Is there a way to extract the breast density of a mammography through Matlab code?
  1 个评论
Adam Danz
Adam Danz 2019-5-15
编辑:Adam Danz 2019-5-15
I've seen this come up several times in the forum. Browse the question/answers that appear in the link above.

请先登录,再进行评论。

采纳的回答

Said Pertuz
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 个评论
Ann G
Ann G 2019-11-15
I have already found a software for breast density, but I will also like to try yours in order to compare the results.
I will let you know of the outcome.
Thank you very much!

请先登录,再进行评论。

更多回答(1 个)

Doaa
Doaa 2024-7-27
  2 个评论
Doaa
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
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!

Translated by