Compute eigen vectors from 3D Hessian matrix in Python

Hi all,

I saw the function itk.SymmetricEigenAnalysisImageFilter can compute the eigen values from a hessian image, but I didn’t see a function that could compute te eigen vectors from a hessian image. Could you help me with this question? I’d like to compute both eigen values and eigen vectors from a hessian image for my task. My data is 3D data and I’m using Python.

Thank you very much!

ITK has a specific vesselness filter, which uses eigenalysis of the Hessian but specifically for tubular structures. There’s an example at
https://examples.itk.org/src/filtering/imagefeature/segmentbloodvessels/documentation
Hope that helps.