Image Center of Gravity (FSL)

Source: R/fslhd.R

Find Center of Gravity of Image from FSL

fslcog(img, mm = TRUE, verbose = TRUE, ts = FALSE)

Arguments

img

Object of class nifti, or path of file
mm

Logical if the center of gravity (COG) would be in mm (default TRUE) or voxels (FALSE)
verbose

(logical) print out command before running
ts

(logical) is the series a timeseries (4D), invoking -t option

Value

Vector of length 3 unless ts option invoked

Note

FSL uses a 0-based indexing system, which will give you a different answer compared to cog, but fslcog(img, mm = FALSE) +1 should be relatively close to cog(img)

Examples

if (have.fsl()){
x = array(rnorm(1e6), dim = c(100, 100, 100))
img = nifti(x, dim= c(100, 100, 100), 
datatype = convert.datatype()$FLOAT32, cal.min = min(x), 
cal.max = max(x), pixdim = rep(1, 4))
fslcog(img)
}