R/antsImageStats.R

In stnava/ANTsR: ANTs in R: Quantification Tools for Biomedical Images

antsImageStats <- function(image, mask = NULL, na.rm = FALSE) {
  if (!is.null(mask)) {
    mask <- check_ants(mask)
  }
  image <- antsImageClone(image, out_pixeltype = "double")
  res <- ANTsRCore::antsImageStats(image, mask, na.rm)
  check <- !any(sapply(res, is.raw))
  if (!check) {
    stop("raw output form antsImageStats, please report this bug")
  }
  res <- lapply(res, drop)
  return(res)
}