Simple motion_correction function.

motion corrects 4D time series imaging data

.motion_correction(
  img,
  fixed = NULL,
  moreaccurate = 1,
  txtype = "Affine",
  verbose = FALSE
)

Arguments

img

4D antsImage

fixed

target fixed image

moreaccurate

0, 1 or 2 with higher values being more accurate (use 2 for real applications, 0 for testing)

txtype

Transformation type

verbose

verbosity boolean

Value

list of outputs

Note

This function may give different results on multiple runs.

Author

Avants BB

Examples

testimg <- makeImage(c(10, 10, 10, 5), rnorm(5000))
testimg <- iMath(testimg, "PadImage", 5)
mocorr <- .motion_correction(testimg)
aimg_to_array <- function(x) {
  if (is.antsImage(x)) {
    x <- as.array(x)
  }
  x
}
amocorr <- lapply(mocorr, aimg_to_array)
mocorr2 <- .motion_correction(testimg)
amocorr2 <- lapply(mocorr2, aimg_to_array)

# This function may give different results on multiple runs
# without setting the seed