Calculating Eigenvectors using Python and ITK


(Veronika) #1

Hey,
I’m trying to calculate eigenvalues and -vectors of a hessian matrix I received with the HessianRecursiveGaussianImageFilter using Python. So far I found neither in SimpleITK nor in ITK an algorithm that would work.
For the SymmetricEigenAnalysisImageFilter and similar filters there seems to be no possibility to access the functionalities via Python.
So are there any algorithms in ITK or SimpleITK I can use to calculate my 3D images’ eigenvectors voxelwise? I need the eigenvectors to compare neighboring voxels to get clues about the structure they are part of.
Thanks, Vero :slight_smile:


(Matt McCormick) #2

Hi Vero :-),

Welcome to ITK!

I will look into wrapping the itk.SymmetricEigenAnalysisImageFilter for you.

Cheers,
Matt