[Insight-users] Mean Curvature Calculation in Level Set

[Joshua Cates]

Hello Sah,

The default curvature calculation in the Level Set function is mean
curvature.  I cannot find the line of code that you reference.  Look in
the file itkLevelSetFunction.txx, method ComputeMeanCurvature().  The
terms for the calculation come from the precomputed Hessian matrix.

Josh.

______________________________
 Josh Cates			
 Scientific Computing and Imaging Institute
 University of Utah
 Email: cates at sci . utah . edu
 Phone: (801) 587-7697
 URL:   http://www . sci . utah . edu/~cates


On Fri, 14 Nov 2003, Sah Rayman wrote:

> I am looking at function:
> 
> template <class TImageType>
> typename LevelSetFunction<TImageType>::ScalarValueType
> LevelSetFunction<TImageType>::ComputeMeanCurvature()
> 
> in itkLevelSetFunction.txx, line 158 (ITK 1.4)
> 
> When I compared that to Sethian's "Level Set Methods
> and Fast Marching Methods", 2ndEd, pp. 70, Eq. 6.36, I
> found that the C++ implementation contain some extra
> terms in nominator:
> 
> -dx[0]*dx[0]*dxy[0][0]-dx[1]*dx[1]*dxy[1][1]-...
> 
> Can someone tell me what are these additional terms
> for? Is the result still mean curvature? Any reference
> for that?
> 
> 
> __________________________________
> Do you Yahoo!?
> Protect your identity with Yahoo! Mail AddressGuard
> http://antispam . yahoo . com/whatsnewfree
> _______________________________________________
> Insight-users mailing list
> Insight-users at itk . org
> http://www . itk . org/mailman/listinfo/insight-users
>