frangi_hessian 多尺度海森矩阵血管增强算法

多尺度海森矩阵血管增强算法,在MATLAB下实现了2维医学图像的血管增强

%求hessian矩阵的特征值（按绝对值的大小排列），并给出边缘方向（最小特征值的特征向量）
function [Lambda1Lambda2IxIy]=eig2image（DxxDxyDyy）

%计算特征向量
Dxx=double（Dxx）;
Dxy=double（Dxy）;
Dyy=double（Dyy）;
temp=sqrt（double（（Dxx-Dyy）.^2+4*Dxy.^2））;
v2_x=2*Dxy;
v2_y=Dyy-Dxx+temp;

%标准化
mag=sqrt（v2_x.^2+v2_y.^2）;
i=（mag~=0）;
v2_x（i）=v2_x（i）./mag（i）;
v2_y（i）=v2_y（i）./mag（i）;

%正交化
v1_x=-v2_y;
v1_y=v2_x;

%计算特征值
mu1=0.5*（Dxx+Dyy+temp）;
mu2=0.5*（Dxx+Dyy-temp）;

%特征值按绝对值大小排序，取特征值大的为lambda1
check=abs（mu1）>abs（mu2）;
Lambda1=mu1;
Lambda1（check）=mu2（check）;
Lambda2=mu2;
Lambda2（check）=mu1（check）;

Ix=v1_x;
Ix（check）=v2_x（check）;
Iy=v1_y;
Iy（check）=v2_y（check）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        723  2016-04-18 08:31  frangi_hessian\eig2image.m

     文件       3528  2016-04-18 14:20  frangi_hessian\frangifilter2D.m

     文件        652  2016-04-18 10:52  frangi_hessian\Hessian2D.m

     目录          0  2016-04-19 15:34  frangi_hessian

----------- ---------  ---------- -----  ----

                 4903                    4