Mammogram - pectoral muscle removal

3 次查看(过去 30 天)
ws
ws 2011-10-8
Hi all
I have a problem with removing pectoral muscle on mammogram. After identified pectoral muscle, what should I do next? I found on some paper, they are using loop function to travel up to First zero pixel and then replace all travelled pixels with zero and in some condition is violating for consecutive 5 times, then decrease pole position for next row by 1 and replace all pixels with zero.
Below is my sample code to identified pectoral muscle area:
pect = I;
[x y] = size(I);
thres=173;
for i=1:x
for j=1:y
if (pect(i,j)<thres)
pect(i,j)=0;
end
end
end
*Note: I is the mammogram which I have removed "white text" area.
Thank you.
  2 个评论
krishnasai gudavalli
How to identify the pectoral muscle?
Image Analyst
Image Analyst 2015-3-20
Post your image in a new question if my answer below about using thresholding does not work for you.

请先登录,再进行评论。

回答(2 个)

Image Analyst
Image Analyst 2011-10-8
I think you should just do
pect(pect <= thres) = 0;
instead of the loops.
  12 个评论
Image Analyst
Image Analyst 2011-10-9
http://www.mathworks.com/matlabcentral/answers/7924-where-can-i-upload-images-and-files-for-use-on-matlab-answers
ws
ws 2011-10-9
This is the image link: http://imageshack.us/photo/my-images/171/tasko.png/
My task:
1. Image 1 until image 5 are the steps I already done, now I want to find the way to get image 6 and 7.
2. I need image 6 to work with original image to get the final image (image 7).
3. After done this task I will use image7 for segmentation part.
Original image for testing purpose can get it from here: http://peipa.essex.ac.uk/ipa/pix/mias/mdb041.pgm
Thank you.

请先登录,再进行评论。


Image Analyst
Image Analyst 2011-10-9
You could threshold your image 5 to get a binary image. Then call imclearborder() to get the image of just the breast. Subtract that from the first binary image to get just the pect, which is touching the border. If you have cases where the breast also touches the border you'll have to do something like find the std dev of the regions and assume that the pect is the lower std dev.
  2 个评论
ws
ws 2011-10-9
Thank you, I will try on it.
ws
ws 2011-10-13
This method is looking for the opposite, I want pectorial muscle, but it removed it.

请先登录,再进行评论。

Community Treasure Hunt

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

Start Hunting!

Translated by