基于Hessian滤波器的血管增强算法

基于Hessian滤波器的血管增强算法

function [Lambda1Lambda2IxIy]=eig2image（DxxDxyDyy）

tmp = sqrt（（Dxx - Dyy）.^2 + 4*Dxy.^2）;
v2x = 2*Dxy;
v2y = Dyy - Dxx + tmp;

mag = sqrt（v2x.^2 + v2y.^2）; 
i = （mag ~= 0）;
v2x（i） = v2x（i）./mag（i）;
v2y（i） = v2y（i）./mag（i）;

v1x = -v2y; 
v1y = v2x;

mu1 = 0.5*（Dxx + Dyy + tmp）;
mu2 = 0.5*（Dxx + Dyy - tmp）;

check=abs（mu1）>abs（mu2）;

Lambda1=mu1; 
Lambda1（check）=mu2（check）;
Lambda2=mu2; 
Lambda2（check）=mu1（check）;

Ix=v1x; 
Ix（check）=v2x（check）;
Iy=v1y; 
Iy（check）=v2y（check）;
end

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

     文件     198710  2011-05-21 13:21  skeleton blood\22.bmp

     文件        524  2011-05-22 20:32  skeleton blood\eig2image.m

     文件       2083  2011-05-22 20:28  skeleton blood\FrangiFilter2D.m

     文件        352  2011-05-22 17:09  skeleton blood\grdAng.m

     文件        433  2011-05-22 20:25  skeleton blood\Hessian2D.m

     文件        455  2011-05-22 17:09  skeleton blood\HessianAng.m

     文件      22500  2011-05-21 12:30  skeleton blood\imgaussian.c

     文件       1118  2011-05-22 20:22  skeleton blood\imgaussian.m

     文件        402  2011-05-22 16:02  skeleton blood\interpll.m

     文件        674  2016-12-13 21:04  skeleton blood\main.m

     文件         69  2011-08-01 17:54  skeleton blood\shuoming.txt

     文件        589  2016-12-13 20:07  skeleton blood\test56.m

     文件        750  2016-12-13 19:53  skeleton blood\testcut.m

     文件       1117  2016-12-13 19:53  skeleton blood\testpp.m

     文件       1781  2016-12-13 19:52  skeleton blood\ThreelinkFilter.m

     目录          0  2016-12-13 20:40  skeleton blood

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

               231557                    16