search

Zhiwei-Zhai / Lung-Vessel-Segmentation-Using-Graph-cuts

GNU General Public License v3.0
29 stars 19 forks source link

Calculation Strain Energy - Relative Hessian Strenght #2

Open Ingrid1989 opened 4 years ago

Ingrid1989 commented 4 years ago

Dear,

I've tried to make the described toolboxes in the other issue work, but unfortunately, I was not able to figure out how to get it working. Luckily, you referent to the articles in which the complete algorithm is described and therefor I'm implementing the code myself in python. According to the article "Lung vessel segmentation in CT images using Graph Cuts" a strain energy filter is used to determine the vesselness. I was wondering If you could explain to me about the measurement of the relative Hessian Strength in this equation (equation (1) from the article). Because the norm of the gradient is just a scalar while the maximum magnitude of the eigenvalues is a matrix. Also when I use this measure, the complete vesselness measure becomes random, while when I leave this measure out of the equation, the outcome actually makes sense.

Can you tell me more about de dimensions of all variables in this hessian strength term and how you calculated them?

Zhiwei-Zhai commented 4 years ago

Ah, I thought you may want to use the development from ITKtools. Sorry, For details about strain energy filter, you can check the literature 'A strain energy filter for 3D vessel enhancement with application to pulmonary CT images'. The C++ source code of strain energy filter could be found in https://github.com/ITKTools/ITKTools/blob/master/src/enhancement/itkStrainEnergyVesselnessImageFilter.h, which might be useful for you to implement it in python. I am not the developer of strain energy filter and I don't know its details.

Ingrid1989 commented 4 years ago

No Problem. I have contact with them and they Helped me to make the toolbox running. Tank you for you're reaction.