matlab 心脏血管中心线提取

基于Matlab程序，对心脏血管中心线进行提取，并实现心血管分割

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  心血管分割\22.bmp

     文件        524  2011-05-22 20:32  心血管分割\eig2image.m

     文件       2083  2011-05-22 20:28  心血管分割\FrangiFilter2D.m

     文件        352  2011-05-22 17:09  心血管分割\grdAng.m

     文件        433  2011-05-22 20:25  心血管分割\Hessian2D.m

     文件        455  2011-05-22 17:09  心血管分割\HessianAng.m

     文件      22500  2011-05-21 12:30  心血管分割\imgaussian.c

     文件       1118  2011-05-22 20:22  心血管分割\imgaussian.m

     文件        402  2011-05-22 16:02  心血管分割\interpll.m

     文件         69  2011-08-01 17:54  心血管分割\shuoming.txt

     文件        538  2011-08-05 11:11  心血管分割\test56.asv

     文件        559  2011-08-05 11:11  心血管分割\test56.m

     文件        724  2011-07-17 18:18  心血管分割\testcut.asv

     文件        724  2011-05-28 08:32  心血管分割\testcut.m

     文件       1088  2011-07-14 22:00  心血管分割\testpp.asv

     文件       1091  2011-08-05 11:08  心血管分割\testpp.m

     文件       1781  2011-05-22 11:22  心血管分割\ThreelinkFilter.m

     目录          0  2011-08-25 09:55  心血管分割

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

               233151                    18