[ITK-users] ComputeMeanCurvature() in itkLevelSetFunction.hxx

[keepdash]

>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.