【图像增强】基于Frangi滤波器实现血管图像增强附matlab代码

1 简介

分析视网膜血管结构的变化是诊断和检测糖尿病,高血压等血管类相关疾病的最重要步骤.针对眼底视网膜图像对比度低,受病变区域边界干扰,很难正确提取血管细节的问题提出了一种基于Frangi滤波器的视网膜血管分割的方法,仿真结果表明上述方法对细小血管的提取表现出良好的效果,具备很强的实用价值.

2 部分代码

function [Dxx,Dxy,Dyy] = Hessian2D(I,Sigma)

% This function Hessian2 Filters the image with 2nd derivatives of a

% Gaussian with parameter Sigma.

%

% [Dxx,Dxy,Dyy] = Hessian2(I,Sigma);

%

% inputs,

% I : The image, class preferable double or single

% Sigma : The sigma of the gaussian kernel used

%

% outputs,

% Dxx, Dxy, Dyy: The 2nd derivatives

%

% example,

% I = im2double(imread('moon.tif'));

% [Dxx,Dxy,Dyy] = Hessian2(I,2);

% figure, imshow(Dxx,[]);

%

% Function is written by D.Kroon University of Twente (June )

if nargin < 2, Sigma = 1; end

% Make kernel coordinates

[X,Y] = ndgrid(-round(3*Sigma):round(3*Sigma));

% Build the gaussian 2nd derivatives filters

DGaussxx = 1/(2*pi*Sigma^4) * (X.^2/Sigma^2 - 1) .* exp(-(X.^2 + Y.^2)/(2*Sigma^2));

DGaussxy = 1/(2*pi*Sigma^6) * (X .* Y).* exp(-(X.^2 + Y.^2)/(2*Sigma^2));

DGaussyy = DGaussxx';

Dxx = imfilter(I,DGaussxx,'conv');

Dxy = imfilter(I,DGaussxy,'conv');

Dyy = imfilter(I,DGaussyy,'conv');

3 仿真结果

4 参考文献

[1]袁盼, 陈以. 基于多尺度Frangi滤波器的视网膜血管分割[J]. 现代信息科技, .

[2]佘黎煌, 郭一蓉, 张石. 基于方向分数和Frangi滤波器的视网膜血管分割算法[J]. 东北大学学报：自然科学版, , 41(2):6.​

博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机、雷达通信、无线传感器等多种领域的Matlab仿真，相关matlab代码问题可私信交流。

部分理论引用网络文献，若有侵权联系博主删除。