[MATLAB图像处理] 乳房X光检查 - 胸肌提取方法

[ooivanoo]

大家好，以下是我所面对的问题：
http://imageshack.us/photo/my-images/171/tasko.png/

我的任务：
1。图像1至5我已经做了的步骤，现在我想找到的方式来获得图像6和7。
2。我需要图像6与原始图像，以获得最终的图像（图片7）。



以下是我在网上找到最简单的“胸肌提取方法”方法,有谁可以帮我转换成matlab code吗？第一至第三我已经完成，但我不是很明白第四个步骤：

Algorithm:
Step 1
Image4 apply local threshold i.e. convert it to binary with Threshold =173. (在第1步你可以找到image5)

Step 2
Travel up to "First non-zero pixel"; Call it as POLE.

Step 3
Now travel up to "First zero pixel" and replace all travelled pixels with zero.

Step 4
Perform Step 1 to 3 for all rows but with following rules:
a) For every POLE ,it should not exceed no. of pixel travelled by previous row
b) Now if at all Rule a is violating for consecutive 5 times then decrease pole position for next row by 1 and replace all pixels with zero up to calculated POLE.



如需测试可以从这里下载原始图像：
http://peipa.essex.ac.uk/ipa/pix/mias/mdb041.pgm


谢谢。