Texture Feature Extraction from a mammography Image
6 次查看(过去 30 天)
显示 更早的评论
I need to extract Texture Features of ROI to be classified whether as benign(normal) or malignant(abnormal).
I used a code, but it's results is not good using different well-known classifiers (SVM kernels for example)
I uploaded a picture that illustrates the following:
1- The first picture titled as('Removed Pectoral'): Breast Image after removing (pectoral muscle)
2- The second picture titled as('ROI'):The Extracted ROI from the removed pectoral muscle image
3- The third picture titled as('segmented ROI): Is the Segmented ROI (Segmentation of ROI image)
My main aim is to obtain the features I uploaded (Additional features in the code)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/166408/image.png)
I used this code for obtaining the features
% Use the Extracted ROI named as "ROI" in the uploaded Images
% Use the following angles "offset":(0,45,90,135)
GLCM = graycomatrix(ROI,'Offset',[0 1;-1 1;-1 0;-1 -1]);
% Calculate the four built-in MATLAB features
stats=graycoprops(GLCM,{'Contrast','Homogeneity','Correlation','Energy'});
contrast=(stats.Contrast);
en=(stats.Energy);
co=(stats.Correlation);
hom=(stats.Homogeneity);
% Calculate mean and standard deviation
m=mean(mean(ROI));
s=std2((ROI));
% The first feature vector
f1=[m s hom co en contrast];
% Calculate the texture Features
% Picture used is" segmented ROI" in the uploaded images
% This picture is the segmentation of " ROI " picture used at the
% beginning of the code
% Features to be extracted " Area, Perimeter, compactness, smoothness, entropy,..."
J2=SegmentedROI;
I=J2;
J2 = uint8(255 * mat2gray(J2));
% Detecting the edges of the ROI
J3=edge(J2,'log');
bw=bwareaopen(J2,150);
bwp=edge(bw,'sobel');
geometric.area=sum(sum(bw));
geometric.peri=sum(sum(bwp));
geometric.compact=geometric.peri^2/geometric.area;
count=1;
roih=0;
for i=1:size(bw,1);
for j=1:size(bw,2);
if bw(i,j)
roih(count,1)=double(SegI(i,j));
count=count+1;
end
end
end
nh=roih/geometric.area;
texture.mean=mean(roih);
texture.globalmean=texture.mean/mean(mean(double(I)));
texture.std=std(roih);
texture.smoothness=1/(1+texture.std^2);
texture.uniformity=sum(nh.^2);
texture.entropy=sum(nh.*log10(nh));
texture.skewness=skewness(roih);
texture.correlation=sum(nh'*roih)-texture.mean/texture.std;
yy=[texture.mean,texture.globalmean,texture.std,texture.smoothness,texture.uniformity,...
texture.entropy,texture.skewness,texture.correlation];
y1=[y1;yy];
xy=[geometric.area,geometric.peri,geometric.compact];
x1=[x1;xy];
f2=[x1 y1];
Totalfeaures=[f1,f1];
2 个评论
Walter Roberson
2017-7-29
编辑:Walter Roberson
2018-11-10
This should have been posted to https://www.mathworks.com/matlabcentral/answers/350414-mammography-image-feature-extraction
Ankur Pruthi
2021-2-22
Hi mohamed ahmed. Can u please share the code for pectoral muscle removal from the image. It would be of great help to me. Mail id- ankpru@gmail.com
回答(8 个)
fatin rasyidah
2017-11-26
hi Mohamed Ahmed, may I know what method that u use to extract ROI and segment the ROI. I hope u can share it, because i am currently doing the same thing like u have did. thanks in advance
5 个评论
Hoang Duy
2019-12-15
Hi Mohamed Ahmed can u send me the MATLAB code for ROI extraction and Segment the ROI with thanks my email duyhm1998@gmail.com
sathya lakshmanan
2018-10-16
Hi Mohamed Ahmed can u send me the MATLAB code for ROI extraction and Segment the ROI with thanks my email sathyass660@gmail.com
0 个评论
hamid hosseinabadi
2018-12-17
Hi Mohamed Ahmed plz send to me your MATLAB code hosseinabadi_hamid@yahoo.com
0 个评论
Loai Kinani
2019-4-17
Hi Mohammad,
can you share the the matlab code of pectoral muscle removal.
Thank you
0 个评论
ahmed nageeb
2019-6-22
hello every one
when use SegI , give me error
Undefined function or variable 'SegI'.
can help me
1 个评论
Walter Roberson
2019-6-23
You are correct, the code that was posted does not define several of the variables.
SegI is not a function: it is a variable the user defined but did not show us the construction of.
mary mwadulo
2019-9-25
hi Mohamed Ahmed,can you kindly share with me the MATLAB code for pectoral muscle removal, extracted ROI and Segmented ROI.Am currently working on the same. my email is mwalowe@yahoo.com thank you.
1 个评论
manish kumari kumari
2020-4-15
hi Mohamed Ahmed,can you kindly share with me the MATLAB code for pectoral muscle removal, extracted ROI and Segmented ROI.Am currently working on the same. my email is mk2634347@gmil.com
Amelougou Abderrahman
2022-5-19
编辑:Amelougou Abderrahman
2022-5-19
hello everyone please I need your help urgently...I worked at first with the Matlab FirstCode code which allows me to delete the pectoral muscle after the background deletion but I would like to see my final image of the same original image size that's why I tried to work on the same code with modifications: 1. orientation of the images right to left and 2. do not remove the background 3 .determine the first and the last non-black pixel in the image i.e. I determine my region of interest ROI 4.application of pectoral muscle removal algorithm .....but unfortunately does not give results requested (muscle removal).....please help
0 个评论
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!