如何使用Matlab提取部分图像

Question

我没有太多使用matlab，我需要从给定的心脏图像中提取心脏的左冠状动脉和右冠状动脉的部分。

这是我的头像，

基于形态学运算，这就是我想出的，

f=imread('heart.jpg');
diam=strel('diamond',19);
line=strel('line',10,90);
linef=imclose(f,line);
line120=strel('line',10,120);
line120f= imclose(f,line120);
bothline=linef+line120f;
diamf=imclose(f,diam);
arterybm=diamf-bothline;
binaryartery= im2bw(arterybm,0);
mask=cast(binaryartery,class(f));
redPlane=f(:,:,1);
greenPlane=f(:,:,2);
bluePlane=f(:,:,3);
maskedRed=redPlane.*mask;
maskedGreen=greenPlane.*mask;
maskedBlue=bluePlane.*mask;
maskedRGBImage=cat(3,maskedRed,maskedGreen,maskedBlue);
subplot(2,3,1);imshow(f);title('Input Image');subplot(2,3,2);imshow(diamf);title('imclose with Diamond Mask');subplot(2,3,3);imshow(bothline);title('imclose with Line 120 and 90 mask');subplot(2,3,4);imshow(arterybm);title('Difference of line and diamond');subplot(2,3,5);imshow(binaryartery);title('Convert to binary image');subplot(2,3,6);imshow(maskedRGBImage);title('Apply mask to input image');

还有更好的方法吗？

Answer 1

这项任务是一项相当困难的任务，如果您发现解决方案在大多数情况下都能完美运行，则值得撰写学术文章。我的建议：搜索有关该主题的文章，并尝试“Matlab 文件交换”(http://www.mathworks.com/matlabcentral/fileexchange/)。如果你很幸运，有人可能已经解决了这个问题并发布了解决方案。

Answer 2

看看 Frangi 过滤器（又名脊线检测过滤器），它的目的是检测血管。 文件交换上有一个可用的实现：

http://www.mathworks.com/matlabcentral/fileexchange/24409-hessian-based-frangi-vesselness-filter