Matlab仅在图像中获取血管

Question

我正在尝试使用图像处理工具箱中提供的功能获取视网膜图像（下图 1），以仅显示黑白血管（下图 2），阈值化、形态学运算等。我尝试了以下方法，但没有得到我需要的结果：

% READ IMAGES 
eye1 = imread('Retina1.png');
%display coloured version(original)
imshow(eye1);

% convert to grayscale
grayeye1 =rgb2gray(eye1);

% display grayscale
imshow(grayeye1);
    
% make darker
dgrayeye1 = imadjust(grayeye1,[0.1 0.9],[]);
imshow(a);
    
se = strel('disk',1);
    
cannyeye = edge(dgrayeye1,'canny',0.15);
figure,imshow(cannyeye);
    
dilate = imdilate(cannyeye,se);
figure, imshow(dilate)

图 1：

图片 2：

Answer 1

这不是一件容易的事。您可能需要在以下几篇参考论文中进行调查：

1. Retinal vessel extraction by matched filter with first-order derivative of Gaussian。而Matlab代码包含在这个post中。

2. An Automatic Hybrid Method for Retinal Blood Vessel Extraction。它使用带有开/关和顶帽变换的形态学方法。之后应用模糊聚类。

3. 最常用的血管检测方法来自 Frangi 等人的论文Hessian-based Multiscale Vessel Enhancement Filtering。有一个 Matlab 实现，FrangiFilter2D，适用于 2D 血管图像。

我在您的情况下尝试了此代码，希望您可以从这里开始...

options=struct('FrangiScaleRange', [1 5], 'FrangiScaleRatio', 1, 'FrangiBetaOne', 1,...
        'FrangiBetaTwo', 7.5, 'verbose',true,'BlackWhite',true);
[outIm,whatScale,Direction] = FrangiFilter2D(double(dgrayeye1), options);
imshow(outIm)

图片：