Find Center of Gravity of Image, after thresholding

cog(img, thresh = 0, ceil = FALSE, warn = TRUE)

Arguments

img

Object of class nifti

thresh

threshold for image, will find img > 0

ceil

Run ceiling to force integers (usu for plotting)

warn

Produce a warning if the image is empty after thresholding

Value

Vector of length 3

Examples

dims = rep(20, 3)
x = array(rnorm(prod(dims)), dim = dims)
img = nifti(x, dim= dims, 
datatype = convert.datatype()$FLOAT32, cal.min = min(x), 
cal.max = max(x), pixdim = rep(1, 4))
cog(img)
#>     dim1     dim2     dim3 
#> 10.54665 10.50265 10.36637