%求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