[Insight-users] itkSymmetricEigenAnalysisImageFilter

[Ruben Schilling]

Hi Jing,

actually it is quite easy to compute the eigenvectors, that you need.  
Since the Hessian matrix is a symmetric tensor of rank 2 (assuming,  
that your derivatives are continuous) you should find what you need  
by including these filters (and maybe some additional smoothing filters)

#include "itkSymmetricSecondRankTensor.h"
#include "itkSymmetricEigenAnalysisImageFilter.h"
#include "itkHessianRecursiveGaussianImageFilter.h"
#include "itkGradientRecursiveGaussianImageFilter.h"

You can then easily do what you proposed.


Best regards
Ruben



Am 20.06.2007 um 15:03 schrieb jing xu:

> II use ITK to extract the centreline of 3D vessel ,The basic theory  
> is to calculate the Hessian matrix of  image at scale
> Then
> lamda1<lamda2<lamda3  and  lamda1<lamda2<0
> The scale product of eigenvector1 and gradient of image is equal to  
> zero
> The scale product of eigenvector2 and gradient of image is equal to  
> zero
> But in ITK software, I only find function  
> itkSymmetricEigenAnalysisImageFilter to deal with 3D  Hessian  
> matrix and there is no method to get the eigenvector
> I want to know there is any other fuction to get the eigenvector of  
> 3D Hessian matrix
> Thank you
> _______________________________________________
> Insight-users mailing list
> Insight-users at itk.org
> http://www.itk.org/mailman/listinfo/insight-users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://public.kitware.com/pipermail/insight-users/attachments/20070620/624be440/attachment.htm