fname <- system.file("extdata", "global_mask_v4.nii", package="neuroim2")
vol <- read_vol(fname)
pmat <- perm_mat(space(vol))
vox <- c(12,12,8)
pvox <- vox %*% perm_mat(space(vol))
stopifnot(all(pvox == c(-12,12,8)))
Run the code above in your browser using DataLab