image curvature for 2D or 3D

uses the weingarten map to estimate image mean or gaussian curvature

weingartenImageCurvature(image, sigma = 1, opt = "mean")

Arguments

image	antsImage
sigma	smoothing parameter
opt	mean by default, otherwise use string gaussian or characterize

Value

image

References

Avants, B, J. Gee, and B. Avants. Shape operator for differential image analysis Information Processing in Medical Imaging, 2003.

Examples

img = makeImage( c(10,10,10) , rnorm( 1000 ) )
fik <- weingartenImageCurvature( img )
if (FALSE) {
fi <- antsImageRead( getANTsRData( "mni" ) )
fik <- weingartenImageCurvature( fi )
}