Discussion:
[ITK-users] ComputeMeanCurvature() in itkLevelSetFunction.hxx
keepdash
2017-06-06 07:57:41 UTC
Permalink
From level set papers, we can find the curvature (2D) is:

K = (fxx*fy*fy + fyy*fx*fx - 2*fx*fy*fxy) / (fx*fx+fy*fy)^(3/2)

I compared the equation with the code in "itkLevelSetFunction.hxx", all same
except the normalization, seems in ITK, the curvature (2D) is computed as:

K = (fxx*fy*fy + fyy*fx*fx - 2*fx*fy*fxy) / (fx*fx+fy*fy)

which can be found in function ComputeMeanCurvature(), line.179:

return ( curvature_term / gd->m_GradMagSqr );

where the "m_GradMagSqr" is the fx*fx+fy*fy from line.332.

Then, why ITK use this way, is it better?

Thank you.



--
View this message in context: http://itk-insight-users.2283740.n2.nabble.com/ITK-users-ComputeMeanCurvature-in-itkLevelSetFunction-hxx-tp7590001.html
Sent from the ITK Insight Users mailing list archive at Nabble.com.
_____________________________________
Powered by www.kitware.com

Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html

Kitware offers ITK Training Courses, for more information visit:
http://www.kitware.com/products/protraining.php

Please keep messages on-topic and check the ITK FAQ at:
http://www.itk.org/Wiki/ITK_FAQ

Follow this link to subscribe/unsubscribe:
http://public.kitware.com/mailman/listinfo/insight-users

Loading...