[xbroyw]

I made a simple hessian matrix using Matlab 704.

syms a b c x
f=a*exp(-(x-b)^2/2/c^2);
Hf = jacobian(jacobian(f));
Hf = simple(Hf); 

however, this is only for scalar values, ie,
subs([Hf,{a,b,c,x},{0,0,0,1})

how can I process my 3D image files using hessian matrix.

This post has been edited by xbroyw: 28 December 2008 - 11:34 PM

[xbroyw]

xbroyw, on 28 Dec, 2008 - 10:30 PM, said:

I made a simple hessian matrix using Matlab 704.

syms a b c x
f=a*exp(-(x-b)^2/2/c^2);
Hf = jacobian(jacobian(f));
Hf = simple(Hf); 

however, this is only for scalar values, ie,
subs([Hf,{a,b,c,x},{0,0,0,1})

how can I process my 3D image files using hessian matrix.

Ok, I solved my question.

I used diff(fxxx),..., diff(fyyy) as mentioned http://www.dreaminco...wtopic53310.htm
then
convolved to my image files.

However, I'm not sure, how to get maximum or minimum eigen value image from hessian matrix convolution.

Thanks to all.

This post has been edited by xbroyw: 29 December 2008 - 03:16 AM