image curvature for 2D or 3D

uses the weingarten map to estimate image mean or gaussian curvature

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

Arguments

image

antsImage

sigma

smoothing parameter

opt

mean by default, otherwise use string gaussian or characterize

labeled

boolean if TRUE then assume image contains labels which denote the surface on which curvature is to be calculated. in this case, the surface should hold integer value 1; the background should be zero; the interior of the surface should be greater than 1 (e.g. 2).

Value

image

References

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

Author

Brian B. Avants

Examples

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