Package 'freesurfer'

Description

Converts Freesurfer aparcs table to brainGraph naming convention, relying on aparcstats2table

Usage

aparcs_to_bg(subjects, measure, clean_col_names = TRUE, ...)

Arguments

Value

Long data.frame

Examples

## Not run: 
fs_subj_dir()
df = aparcs_to_bg(subjects = "bert", measure = "thickness")
print(head(df))

## End(Not run)

Description

Ensures the output to be a character filename (or vector) from an input image or nifti to have .mnc extension and be converted to MNC when necessary

Usage

checkmnc(file, ...)

## S4 method for signature 'nifti'
checkmnc(file, ...)

## S4 method for signature 'character'
checkmnc(file, ...)

## S4 method for signature 'list'
checkmnc(file, ...)

checkmnc(file, ...)

Arguments

Value

Character filename of mnc image

Author(s)

John Muschelli [email protected]

Description

This function copies files specified by the types of data, determined by the folder Freesurfer put them in, into a temporary directory for easier separation of data and different structuring of data.

Usage

construct_subj_dir(
  label = NULL,
  mri = NULL,
  stats = NULL,
  surf = NULL,
  touch = NULL,
  subj = NULL,
  subj_root_dir = tempdir(check = TRUE)
)

Arguments

Value

List with the subject name, the SUBJECTS_DIR to use (the directory that contains the subject name), and the types of objects copied

Examples

## Not run: 
library(freesurfer)
label = "/Applications/freesurfer/subjects/bert/label/aparc.annot.a2009s.ctab"
mri = c(
  "/Applications/freesurfer/subjects/bert/mri/aparc.a2009s+aseg.mgz",
  "/Applications/freesurfer/subjects/bert/mri/aseg.auto.mgz")
stats = c("/Applications/freesurfer/subjects/bert/stats/lh.aparc.stats",
          "/Applications/freesurfer/subjects/bert/stats/aseg.stats")
surf = "/Applications/freesurfer/subjects/bert/surf/lh.thickness"
construct_subj_dir(
  label = label,
  mri = mri,
  stats = stats,
  surf = surf)

## End(Not run)

Description

Reads in a surface file from Freesurfer and separates into vertices and faces

Usage

convert_surface(infile, ...)

Arguments

Value

List of 3 elements: a header indicating the number of vertices and faces, the vertices, and the faces

Note

This was adapted from the gist: https://gist.github.com/mm–/4a4fc7badacfad874102

Examples

## Not run: 
infile = file.path(fs_subj_dir(),
                   "bert", "surf", "rh.pial")
res = convert_surface(infile = infile)

## End(Not run)

Description

Reads a Freesurfer curvature file according to the FREESURFER_HOME/matlab/read_curv.m file.

Usage

freesurfer_read_curv(file)

Arguments

Value

Numeric vector

Examples

bert_dir = file.path(fs_subj_dir(), "bert", "surf")
file = file.path(bert_dir, "lh.thickness")
fid = file(file, open = "rb")
out = freesurfer_read_curv(file)

Description

Reads a Freesurfer Surface file from the surf/ directory from recon-all

Usage

freesurfer_read_surf(file, verbose = get_fs_verbosity())

Arguments

Value

List of length 2: vertices and faces are the elements

Examples

fname = file.path(fs_subj_dir(), "bert", "surf", "lh.inflated")
out = freesurfer_read_surf(fname)

Description

fs_cmd() is the core wrapper function that executes FreeSurfer command-line tools from within R. It handles file path validation, command construction, execution, and optional return of neuroimaging data as nifti objects.

Usage

fs_cmd(
  func,
  file,
  outfile = NULL,
  retimg = TRUE,
  reorient = FALSE,
  verbose = get_fs_verbosity(),
  intern = verbose,
  opts = "",
  opts_after_outfile = FALSE,
  frontopts = "",
  bin_app = "bin",
  timeout_seconds = NULL,
  validate_inputs = TRUE,
  ...
)

Arguments

Details

This function provides a unified interface to FreeSurfer's command-line tools, with several key features:

• Automatic validation of input/output files

• Support for nifti objects as inputs

• Optional timeout for long-running commands

• Flexible command option placement

• Automatic result verification

• Integration with FreeSurfer environment setup

Overwriting Files

To overwrite the input file, set outfile = file. A warning will be displayed for safety. The function checks file existence before and after command execution to verify successful completion.

Command Construction

The command is constructed in the following order:

1. FreeSurfer environment setup (from get_fs())

2. Command name (func)

3. Front options (frontopts)

4. Input file path (quoted)

5. Output file path (quoted, if provided)

6. Additional options (opts)

Use opts_after_outfile = TRUE to place opts after the output file.

Value

If retimg = TRUE: A nifti object containing the output image, or NULL if output creation failed. If retimg = FALSE: The result from base::system(), typically the exit status or captured output.

See Also

get_fs() for FreeSurfer environment setup, mri_convert() for format conversion, neurobase::readnii() for reading nifti files

Examples

## Not run: 
# Basic usage: convert MGZ to NIfTI
fs_cmd(
  func = "mri_convert",
  file = "input.mgz",
  outfile = "output.nii.gz",
  opts = "--conform"
)

# Return as nifti object
img <- fs_cmd(
  func = "mri_convert",
  file = "input.mgz",
  retimg = TRUE
)

# Use nifti object as input
library(oro.nifti)
img <- nifti(array(rnorm(10*10*10), dim = c(10, 10, 10)))
result <- fs_cmd(
  func = "mri_convert",
  file = img,
  outfile = "output.nii.gz"
)

## End(Not run)

Description

Functions to retrieve FreeSurfer's installation and subjects directories. Multiple aliases are provided for convenience and backward compatibility.

Usage

freesurferdir()

freesurfer_dir()

fs_dir()

fs_subj_dir()

Value

Character; path to the FreeSurfer home or subjects directory.

Functions

• freesurferdir(): Get FreeSurfer installation directory

• freesurfer_dir(): Alias for freesurferdir

• fs_dir(): Short alias for freesurferdir

• fs_subj_dir(): Get FreeSurfer subjects directory

See Also

get_fs_home() for more detailed information, have_fs() to check if FreeSurfer is accessible

Examples

# All return the same value
freesurferdir()
freesurfer_dir()
fs_dir()

# Get subjects directory
fs_subj_dir()

Description

This function takes in the function and returns the help from FreeSurfer for that function with simple validation.

Usage

fs_help(
  func_name,
  help_arg = "--help",
  extra_args = "",
  timeout_seconds = 30,
  display = TRUE,
  warn = TRUE,
  bin_app = "bin",
  ...
)

Arguments

Value

Prints help output and returns output as character vector

Examples

fs_help("mri_watershed")

# For MNI tools
fs_help("nu_correct", help_arg = "-help", bin_app = "mni/bin")

Description

Runs get_fs_output() to extract the Freesurfer output type and then gets the corresponding file extension (such as .nii.gz).

Usage

fs_imgext()

Value

Character string representing the file extension for the output type.

Examples

fs_imgext()

Description

A data.frame with the index, label, and RGBA (red, blue, green, alpha) specification for the segmentations

Usage

fs_lut

Format

An object of class data.frame with 1266 rows and 6 columns.

Description

Get a situation report on your current FreeSurfer installation and settings within R. This function helps diagnose problems by showing how FreeSurfer paths and options are determined.

Usage

fs_sitrep(clear_cache = FALSE, test_commands = TRUE)

Arguments

Examples

# Report all FreeSurfer settings
fs_sitrep()

Description

Finds the Freesurfer version from FREESURFER_HOME/build-stamp.txt

Usage

fs_version()

Value

If the version file does not exist, it will throw a warning, but it will return an empty string. Otherwise it will be a string of the version.

Note

This will use fs_dir() to get the directory of FREESURFER

Examples

fs_version()

Description

This function generates a bash command string to set up the environment for using FreeSurfer. It ensures the required FreeSurfer installation, license, and environment setup files are validated and included in the command string. The function handles different FreeSurfer binaries like bin, mni/bin, and others, while ensuring proper initialization of the MNI environment if required.

Usage

get_fs(bin_app = c("bin", "mni/bin"), fs_home = get_fs_home(simplify = FALSE))

Arguments

Value

character A bash command string that includes environment setup for FreeSurfer. If the FreeSurfer environment or required configurations cannot be initialized, the function throws an error or issues a warning. On success, the returned string can be used directly in shell operations to load the FreeSurfer environment.

See Also

get_fs_home(), get_fs_license(), get_fs_output()

Examples

# Generate a shell command to set up FreeSurfer with the default `bin`
get_fs(bin_app = "bin")

# Generate a shell command to include MNI environment setup
get_fs(bin_app = "mni/bin")

Description

These functions retrieve FreeSurfer configuration settings using a hierarchical lookup system: R options → environment variables → defaults.

Usage

get_fs_setting(env_var, opt_var, defaults = NULL, is_path = TRUE)

get_fs_home(simplify = TRUE)

get_fs_license(fs_home = get_fs_home(), simplify = TRUE)

get_fs_subdir(fs_home = get_fs_home(), simplify = TRUE)

get_fs_source(fs_home = get_fs_home(), simplify = TRUE)

get_fs_verbosity(simplify = TRUE)

get_fs_output(simplify = TRUE)

get_mni_bin(fs_home = get_fs_home(), simplify = TRUE)

Arguments

Details

Lookup Hierarchy

Settings are resolved in the following order (first match wins):

1. R options (set with options())

2. Environment variables (set with Sys.setenv())

3. Default paths (platform-specific)

Setting Options

You can set R options in your .Rprofile:

options(
  freesurfer.home = "/usr/local/freesurfer",
  freesurfer.subj_dir = "~/freesurfer_subjects",
  freesurfer.verbose = TRUE
)

Value

When simplify = TRUE: Character string with the setting value, or NA. When simplify = FALSE: List with value, source, and exists components.

Functions

• get_fs_setting(): Core function for retrieving settings

• get_fs_home(): Retrieve FreeSurfer installation directory

• get_fs_license(): Retrieve FreeSurfer license file path

• get_fs_subdir(): Retrieve FreeSurfer subjects directory

• get_fs_source(): Retrieve FreeSurfer source script path

• get_fs_verbosity(): Retrieve FreeSurfer verbosity setting

• get_fs_output(): Retrieve FreeSurfer output format

• get_mni_bin(): Retrieve MNI tools directory

See Also

have_fs() to check if FreeSurfer is accessible, fs_sitrep() for comprehensive diagnostics, get_fs() to generate environment setup commands

Examples

# Get FreeSurfer home directory
fs_home <- get_fs_home()

# Get detailed information
fs_home_info <- get_fs_home(simplify = FALSE)

# Check license file
license <- get_fs_license()

# Get subjects directory
subj_dir <- get_fs_subdir()

# Check output format
format <- get_fs_output()  # Returns "nii.gz", "nii", "mgz", etc.

# Check verbosity
verbose <- get_fs_verbosity()  # Returns TRUE/FALSE

## Not run: 
# Set custom options
options(freesurfer.home = "/custom/path/freesurfer")
get_fs_home()  # Returns: "/custom/path/freesurfer"

## End(Not run)

Description

Checks if FreesurferDIR is accessible and optionally if a license file exists.

Usage

have_fs(check_license = TRUE)

Arguments

Value

Logical TRUE if Freesurfer is accessible and license (if checked) is found, FALSE otherwise.

Examples

have_fs()

Description

This function calls mnc2nii to convert MNC files to NIfTI

Usage

mnc2nii(file, outfile = NULL, ...)

mnc2nii.help()

Arguments

Value

Character filename of output

Functions

• mnc2nii.help(): Display information about mnc2nii command

Examples

img = oro.nifti::nifti(array(rnorm(5*5*5), dim = c(5,5,5)))
mnc = nii2mnc(img)
img_file = mnc2nii(mnc, outfile = temp_file(fileext = ".nii"))
neurobase::readnii(img_file, verbose = get_fs_verbosity())

Description

Calls FreeSurfer's mri_convert to convert between different medical imaging formats (NIfTI, MGZ, MINC, etc.) with enhanced format validation and error handling.

Usage

mri_convert(
  file,
  outfile,
  opts = "",
  format_check = TRUE,
  timeout_seconds = 300,
  ...
)

mri_convert.help(...)

Arguments

Details

Runtime FreeSurfer CLI help is available via the helper function mri_convert.help(). When called, that helper will attempt to fetch and display the underlying FreeSurfer command-line help if FreeSurfer is installed on the system. The package does not embed the CLI help into the installed Rd at package build time.

This function provides format validation and informative messages about conversions between different image types.

Value

Result of base::system() command, typically exit status.

Functions

• mri_convert.help(): Display FreeSurfer help for mri_convert

FreeSurfer Command Help

When FreeSurfer is installed and available, detailed command-line help for the underlying mri_convert command can be accessed via mri_convert.help(). This provides comprehensive information about FreeSurfer's native command options, file format support, and usage examples.

To see FreeSurfer CLI help: mri_convert.help()

See Also

fs_cmd() for the underlying command wrapper

Examples

# Convert nifti object to MGZ format
img <- oro.nifti::nifti(array(rnorm(5*5*5), dim = c(5, 5, 5)))
res <- mri_convert(img, outfile = temp_file(fileext = ".mgz"))

## Not run: 
# Convert MGZ to NIfTI
mri_convert("brain.mgz", "brain.nii.gz")

# With additional options
mri_convert("brain.mgz", "brain.nii.gz", opts = "--conform")

## End(Not run)

Description

This calls Freesurfer's mri_deface

Usage

mri_deface(file, brain_template = NULL, face_template = NULL, ...)

Arguments

Value

Result of fs_cmd, which type depends on arguments to ...

Note

If brain_template orface_template is NULL, they will be downloaded.

Examples

## Not run: 
base_url = "https://surfer.nmr.mgh.harvard.edu/pub/dist/mri_deface"
url = file.path(base_url, "sample_T1_input.mgz")
x = temp_file(fileext = ".mgz")
out = try({
  utils::download.file(url, destfile = x)
})
if (!inherits(out, "try-error")) {
  noface = mri_deface(x)
} else {
  url = paste0(
    "https://raw.githubusercontent.com/muschellij2/kirby21.t1/master/",
    "inst/visit_1/113/113-01-T1.nii.gz")
   x = temp_file(fileext = ".nii.gz")
   out = try({
     utils::download.file(url, destfile = x)
   })
   noface = mri_deface(x)
 }

## End(Not run)

Description

This calls Freesurfer's mri_info

Usage

mri_info(file, ...)

Arguments

Value

Result of fs_cmd, which type depends on arguments to ...

Examples

## Not run: 
img = oro.nifti::nifti(array(rnorm(5*5*5), dim = c(5,5,5)))
mri_info(img)

## End(Not run)

Description

This calls Freesurfer's mri_info help

Usage

mri_info.help(...)

Arguments

Description

This function calls mri_mask to mask an image

Usage

mri_mask(file, mask, outfile = NULL, retimg = TRUE, opts = "", ...)

Arguments

Value

Character or nifti depending on retimg

Examples

## Not run: 
img = oro.nifti::nifti(array(rnorm(5*5*5), dim = c(5,5,5)))
mask = img > 1
res = mri_mask(img, mask)

## End(Not run)

Description

This calls Freesurfer's mri_mask help

Usage

mri_mask.help(...)

Arguments

Description

This function calls mri_normalize to normalize the values of the image, with white matter voxels around 110.

Usage

mri_normalize(file, outfile = NULL, retimg = TRUE, opts = "", ...)

Arguments

Value

Character or nifti depending on retimg

Examples

## Not run: 
mri_normalize("/path/to/T1.nii.gz")

## End(Not run)

Description

This calls Freesurfer's mri_normalize help

Usage

mri_normalize.help(...)

Arguments

Value

Result of fs_help

Description

This function calls mri_segment

This calls Freesurfer's mri_segment help

Usage

mri_segment(file, outfile = NULL, retimg = TRUE, opts = "", ...)

mri_segment.help(...)

Arguments

Value

Character or nifti depending on retimg

Result of fs_help

Functions

• mri_segment.help(): Display FreeSurfer help for mri_segment

Note

NOT COMPLETE

Description

Calls FreeSurfer's mri_surf2surf to resample one cortical surface onto another, enabling comparison of surface data across subjects or atlases.

Usage

mri_surf2surf(
  subject = NULL,
  target_subject = NULL,
  trg_type = c("curv", "w", "mgh", "nii"),
  src_type = c("curv", "w"),
  outfile = NULL,
  hemi = c("lh", "rh"),
  sval = c("thickness"),
  subj_dir = NULL,
  opts = "",
  verbose = get_fs_verbosity(),
  ...
)

mri_surf2surf.help(...)

Arguments

Details

This function is commonly used to:

• Project individual subject data onto a template (e.g., fsaverage)

• Resample between different surface resolutions

• Convert between surface file formats

The output filename is automatically prefixed with the hemisphere (e.g., "lh.output.mgz").

Value

Character; path to the output file (prefixed with hemisphere).

Functions

• mri_surf2surf.help(): Display FreeSurfer help for mri_surf2surf

See Also

fs_help() for FreeSurfer command documentation

Examples

## Not run: 
out <- mri_surf2surf(
  subject = "bert",
  target_subject = "fsaverage",
  trg_type = "curv",
  src_type = "curv",
  hemi = "rh",
  sval = "thickness"
)

## End(Not run)

Description

This function calls mri_mask to mask an image

Usage

mri_synthstrip(
  file,
  outfile = NULL,
  retimg = TRUE,
  maskfile = NULL,
  opts = "",
  ...
)

synthstrip(
  file,
  outfile = NULL,
  retimg = TRUE,
  maskfile = NULL,
  opts = "",
  ...
)

Arguments

Value

Character or nifti depending on retimg

Examples

## Not run: 
mock_nifti = array(rnorm(5*5*5), dim = c(5,5,5))
img = oro.nifti::nifti(mock_nifti)
res = mri_synthstrip(img)

## End(Not run)

Description

This calls Freesurfer's mri_mask help

Usage

mri_synthstrip.help(...)

Arguments

Value

Result of fs_help

Description

This function calls mri_watershed to extract a brain from an image, usually for skull stripping.

This calls Freesurfer's mri_watershed help

Usage

mri_watershed(file, outfile = NULL, retimg = TRUE, opts = "", ...)

mri_watershed.help(...)

Arguments

Value

Character or nifti depending on retimg

Result of fs_help

Functions

• mri_watershed.help(): Display FreeSurfer help for mri_watershed

Examples

## Not run: 
mri_watershed("/path/to/T1.nii.gz")

## End(Not run)

Description

These functions call FreeSurfer's mris_convert to convert between cortical surface file formats. The base function mris_convert() provides general conversion, while specialized variants handle specific data types.

Usage

mris_convert(
  infile,
  outfile = NULL,
  ext = ".asc",
  opts = "",
  verbose = get_fs_verbosity(),
  ...
)

mris_convert_annot(annot, opts = "", ...)

mris_convert_curv(curv, opts = "", ...)

mris_convert_normals(opts = "", ...)

mris_convert_vertex(opts = "", ...)

mris_convert.help(...)

Arguments

Details

FreeSurfer's mris_convert is a general conversion program for cortical surface file formats. It can convert between binary and ASCII formats, extract specific data types, and transform between coordinate systems.

Value

Character; path to the output file. For mris_convert(), the output has a "separator" attribute indicating the file's field separator.

Functions

• mris_convert_annot(): Convert with annotation or gifti label data

• mris_convert_curv(): Convert with scalar curvature overlay

• mris_convert_normals(): Extract surface normals

• mris_convert_vertex(): Extract vertex coordinates

• mris_convert.help(): Display FreeSurfer help for mris_convert

Note

For curvature conversions, the output filename may be modified by FreeSurfer. You may need to prepend the hemisphere prefix to get the correct path.

Examples

## Not run: 
bert_dir <- file.path(fs_subj_dir(), "bert")
bert_surf <- file.path(bert_dir, "surf", "lh.white")

# Basic conversion to ASCII
asc_file <- mris_convert(infile = bert_surf)

# Convert with annotation data
gii_file <- mris_convert_annot(
  infile = bert_surf,
  annot = file.path(bert_dir, "label", "lh.aparc.annot"),
  ext = ".gii"
)

# Convert with curvature overlay
curv_file <- mris_convert_curv(
  infile = bert_surf,
  curv = file.path(bert_dir, "surf", "lh.thickness")
)

# Extract surface normals
normals_file <- mris_convert_normals(infile = bert_surf)

# Extract vertex coordinates
vertex_file <- mris_convert_vertex(infile = bert_surf)

## End(Not run)

Description

This function calls mris_euler_number to calculate the Euler Number with improved error handling and output management.

Usage

mris_euler_number(
  file,
  outfile = NULL,
  opts = "",
  cleanup_temp = TRUE,
  timeout_seconds = 120,
  validate_format = TRUE,
  ...
)

Arguments

Value

Character vector containing the Euler number results

Examples

## Not run: 
img = oro.nifti::nifti(array(rnorm(5*5*5), dim = c(5,5,5)))
res = mris_euler_number(img, outfile = temp_file(fileext = ".mgz"))

## End(Not run)

Description

This calls FreeSurfer's mris_euler_number help

Usage

mris_euler_number.help(...)

Arguments

Description

This function calls nii2mnc to convert NIfTI to MNC files

Usage

nii2mnc(file, outfile = tempfile(fileext = ".mnc"), ...)

nii2mnc.help()

Arguments

Value

Character filename of output

Functions

• nii2mnc.help(): Display information about nii2mnc command

Examples

img = oro.nifti::nifti(array(rnorm(5*5*5), dim = c(5,5,5)))
mnc = nii2mnc(img)
img_file = mnc2nii(mnc)

Description

This function calls nu_correct to correct for intensity non-uniformity in MRI images using N3 (nonparametric nonuniform normalization).

Usage

nu_correct(file, mask = NULL, opts = "", verbose = get_fs_verbosity(), ...)

nu_correct.help(...)

Arguments

Details

The nu_correct command performs N3 (nonparametric nonuniform intensity normalization) to correct for intensity inhomogeneity in MRI images. This is particularly useful for T1-weighted images where intensity varies across the image due to magnetic field inhomogeneities.

Common options:

• ⁠-mask <file>⁠: Use a binary mask

• ⁠-distance <value>⁠: N3 spline distance (default 200mm)

• ⁠-iterations <value>⁠: Maximum iterations (default 4)

• ⁠-stop <value>⁠: Stopping criterion (default 0.01)

Note: This is an MNI tool and uses different help syntax (-help instead of --help).

Value

Object of class nifti.

Functions

• nu_correct.help(): Display information about nu_correct

See Also

mri_normalize() for FreeSurfer's normalization, nu_correct.help() for command information

Examples

## Not run: 
# Basic usage
nu_correct("/path/to/T1.nii.gz")

# With mask
nu_correct("/path/to/T1.nii.gz", mask = "/path/to/mask.nii.gz")

## End(Not run)

Description

Reads Freesurfer binary annotation files that contain information on vertex labels and colours for use in analyses and brain area lookups.

Usage

read_annotation(path, verbose = get_fs_verbosity())

Arguments

Details

This function is heavily based on Freesurfer's read_annotation.m Original Author: Bruce Fischl CVS Revision Info: $Author: greve $ $Date: 2014/02/25 19:54:10 $ $Revision: 1.10 $

Value

list of 3 with vertices, labels, and colortable

Examples

bert_dir = file.path(fs_subj_dir(), "bert")
annot_file = file.path(bert_dir, "label", "lh.aparc.annot")
res = read_annotation(annot_file)

Description

Reads and parses an aseg.stats file produced by FreeSurfer's anatomical segmentation pipeline. Returns both global brain measures and structure- specific statistics.

Usage

read_aseg_stats(file, lowercase = TRUE)

Arguments

Details

The aseg.stats file contains volumetric and morphometric statistics for subcortical structures and global brain measures computed by FreeSurfer's recon-all pipeline.

Output Structure

The function returns a list with two data frames:

measures: Global brain measures including:

• Brain segmentation volume

• Estimated total intracranial volume (eTIV)

• Normalization factors

structures: Structure-specific measures including:

• Structure index and name

• Number of voxels

• Volume (mm³)

• Mean intensity and standard deviation

• Minimum, maximum, and range of intensities

File Location

For a subject processed with recon-all, the file is typically located at: ⁠$SUBJECTS_DIR/<subject_id>/stats/aseg.stats⁠

Value

A list with two components:

• measures: A data frame with global brain measures, with columns:

  • measure: Short name of the measure

  • measure_long: Long descriptive name

  • meaning: Description of what the measure represents

  • value: Numeric value of the measure

  • units: Units of measurement

• structures: A data frame with structure-specific statistics, with columns:

  • Index: Structure index

  • SegId: Segmentation ID

  • NVoxels: Number of voxels

  • Volume_mm3: Volume in cubic millimeters

  • StructName: Name of the anatomical structure

  • normMean: Normalized mean intensity

  • normStdDev: Normalized standard deviation

  • normMin: Normalized minimum

  • normMax: Normalized maximum

  • normRange: Normalized range

See Also

• read_fs_table() for reading other FreeSurfer tables

• recon_all() for running the FreeSurfer pipeline

• stats2table() for combining stats from multiple subjects

Examples

# Read stats for the "bert" example subject
file <- file.path(fs_subj_dir(), "bert", "stats", "aseg.stats")

if (file.exists(file)) {
  out <- read_aseg_stats(file)

  # Examine global measures
  print(head(out$measures))

  # Get estimated total intracranial volume
  etiv <- out$measures$value[out$measures$measure == "estimatedtotalintracranialvol"]
  cat("eTIV:", etiv, "mm³\n")

  # Examine structure-specific stats
  print(head(out$structures))

  # Get hippocampal volumes
  hippo <- out$structures[grepl("Hippocampus", out$structures$StructName), ]
  print(hippo[, c("StructName", "Volume_mm3")])
}

Description

Reads an label file from an individual subject

Usage

read_fs_label(file)

Arguments

Value

data.frame with 5 columns:

vertex_num:

Vertex Number

r_coord:

Coordinate in RL direction

a_coord:

Coordinate in AP direction

s_coord:

Coordinate in SI direction

value:

Value of label (depends on file)

Examples

file = file.path(fs_subj_dir(), "bert", "label", "lh.BA1.label")
 if (!file.exists(file)) {
 file = file.path(fs_subj_dir(), "bert", "label", "lh.BA1_exvivo.label")
 }
 out = read_fs_label(file)

Description

This function reads output from a FreeSurfer table command, e.g. aparcstats2table, asegstats2table

Usage

read_fs_table(
  file,
  sep = NULL,
  header = TRUE,
  validate_format = TRUE,
  stringsAsFactors = FALSE,
  ...
)

read_stats_table(
  file,
  sep = NULL,
  header = TRUE,
  validate_format = TRUE,
  stringsAsFactors = FALSE,
  ...
)

Arguments

Value

data.frame from the file

Examples

outfile = aparcstats2table(
   subjects = "bert",
   hemi = "lh",
   meas = "thickness"
)
df = read_fs_table(outfile)
seg_outfile = asegstats2table(subjects = "bert", meas = "mean")
df_seg = read_fs_table(seg_outfile)

Description

This function calls mri_convert to convert MGH/MGZ files to NIfTI, then reads it in using readnii with enhanced validation.

Usage

read_mgz(file, validate_format = TRUE, cleanup_temp = TRUE, ...)

read_mgh(file, validate_format = TRUE, cleanup_temp = TRUE, ...)

Arguments

Value

Object of class nifti

Description

This function calls mnc2nii to convert MNC files to NIfTI, then reads it in using readnii

Usage

read_mnc(file)

Arguments

Value

Object of class nifti

Description

Functions to run FreeSurfer's cortical reconstruction pipeline with varying levels of control and customization.

Usage

recon_all(
  subjid = NULL,
  infile = NULL,
  outdir = NULL,
  verbose = get_fs_verbosity(),
  opts = "-all",
  ...
)

recon(
  subjid = NULL,
  infile = NULL,
  outdir = NULL,
  opts = "",
  options = recon_steps(),
  verbose = get_fs_verbosity()
)

recon_steps()

Arguments

Details

Pipeline Overview

The reconstruction pipeline performs:

1. Motion correction and intensity normalization

2. Skull stripping and subcortical segmentation

3. White matter segmentation

4. Surface generation (white and pial)

5. Cortical parcellation and thickness calculation

Processing Time

A full reconstruction typically takes 6-24 hours. Use opts = "-parallel" for parallel processing if available.

Output Location

Results stored in ⁠$SUBJECTS_DIR/<subjid>/⁠ with subdirectories: ⁠mri/⁠, ⁠surf/⁠, ⁠label/⁠, ⁠stats/⁠

Value

Result from base::system() call, typically exit status (0 = success).

Functions

• recon_all(): Run complete reconstruction pipeline with default settings

• recon(): Run reconstruction with step-by-step control

• recon_steps(): Get named vector of available reconstruction steps

References

Dale et al. (1999) Neuroimage 9:179-194. Fischl et al. (1999) Neuroimage 9:195-207.

See Also

read_aseg_stats() to read segmentation statistics reconner() for low-level control of recon-all

Examples

## Not run: 
# Full reconstruction with recon_all
recon_all(
  infile = "T1_scan.nii",
  subjid = "subject01",
  outdir = "~/freesurfer_output"
)

# Step-by-step control with recon
steps <- recon_steps()
steps["normalization"] <- FALSE
recon(
  infile = "T1_scan.nii",
  subjid = "subject02",
  options = steps
)

# Low-level control with reconner
reconner(
  infile = "T1_scan.nii",
  subjid = "subject03",
  opts = "-autorecon1 -parallel"
)

## End(Not run)

Description

A set of helper functions for running different stages of Freesurfer's recon-all pipeline: autorecon1, autorecon2, and autorecon3. These functions are wrappers around the reconner function, simplifying the reconstruction process by specifying pre-defined options for each stage.

Usage

recon_con1(infile, subjid, outdir = NULL, verbose = get_fs_verbosity())

autorecon1(infile, subjid, outdir = NULL, verbose = get_fs_verbosity())

recon_con2(infile, subjid, outdir = NULL, verbose = get_fs_verbosity())

autorecon2(infile, subjid, outdir = NULL, verbose = get_fs_verbosity())

recon_con3(infile, subjid, outdir = NULL, verbose = get_fs_verbosity())

autorecon3(infile, subjid, outdir = NULL, verbose = get_fs_verbosity())

Arguments

Details

Freesurfer is a widely used software suite for processing and analyzing brain MRI images. The recon-all tool provides a stepwise reconstruction process:

• -autorecon1: Motion correction to skull stripping (steps 1-5)

• -autorecon2: Further processing, including intensity normalization

• -autorecon3: Final surface generation and QC measures

These helpers call specific stages of the recon-all pipeline to provide finer control over processing.

Value

The result of the system call to recon-all, either the command's output or a diagnostic message if an issue occurs.

Functions

• recon_con1(): Performs -autorecon1, which includes steps 1-5 of Freesurfer's reconstruction process.

• autorecon1(): Helper function for recon_con1 with the same functionality.

• recon_con2(): Performs -autorecon2, which includes steps 6+ of Freesurfer's reconstruction process.

• autorecon2(): Helper function for recon_con2 with the same functionality.

• recon_con3(): Performs the final stage of Freesurfer's reconstruction process, which includes generating cortical surfaces,

• autorecon3(): Helper function for recon_con3 with the same functionality.

Note

See the official Freesurfer documentation for additional details on each autorecon stage: https://surfer.nmr.mgh.harvard.edu/fswiki/recon-all.

If infile is set to NULL, the -i flag (input file flag) in recon-all will be omitted, which may or may not be appropriate depending on your use case.

See Also

reconner, recon, recon_all

Examples

## Not run: 
# Example: Perform Step 1-5 reconstruction (Motion Correction to Skull Strip)
recon_con1(
  infile = "subject_001.nii",
  outdir = "/output_dir",
  subjid = "subj01",
  verbose = TRUE
)

# Example: Run autorecon1 with equivalent helper function
autorecon1(
  infile = "subject_001.nii",
  outdir = "/output_dir",
  subjid = "subj01",
  verbose = TRUE
)

# Example: Perform further processing for Step 6+ using autorecon2
autorecon2(
  infile = "subject_002.nii",
  outdir = "/output_dir",
  subjid = "subj02",
  verbose = TRUE
)

# Example: Complete reconstruction with autorecon3
autorecon3(
  infile = "subject_003.nii",
  outdir = "/output_dir",
  subjid = "subj03",
  verbose = FALSE
)

## End(Not run)

Description

Wrapper around FreeSurfer's recon-all command for brain surface reconstruction from MRI data. Handles input processing, subject directory creation, and command flag management.

Usage

reconner(
  subjid = NULL,
  infile = NULL,
  outdir = NULL,
  opts = "-all",
  force = FALSE,
  verbose = get_fs_verbosity()
)

Arguments

Details

FreeSurfer's recon-all performs cortical reconstruction and volumetric segmentation. This function simplifies usage by:

• Automatically deriving subject ID from input filename if not provided

• Managing subject directory paths

• Providing force option for re-running on existing subjects

Value

Result of base::system() call, typically exit status (0 = success).

See Also

tracker() for diffusion tractography pipeline

Examples

## Not run: 
reconner(infile = "input.nii", outdir = "/output_dir", subjid = "subj01")
reconner(infile = "input.nii", outdir = "/output_dir")
reconner(infile = "input.nii", opts = "-autorecon2", force = TRUE)

## End(Not run)

Description

This function serves as a flexible abstraction to execute the FreeSurfer commands asegstats2table or aparcstats2table, which are primarily used to convert parcellation and segmentation statistics to tabular formats. The function dynamically handles shared logic, constructs command-line arguments, and runs the appropriate FreeSurfer command. Users can specify the type of input (subjects or input file paths) and various configuration options. With the appropriate parameters, it constructs and executes the matching system command.

Usage

stats2table(
  type = c("aseg", "aparc"),
  input,
  measure,
  input_type = c("subjects", "inputs"),
  outfile = NULL,
  delim = c("tab", "space", "comma", "semicolon"),
  skip = TRUE,
  subj_dir = NULL,
  opts = "",
  verbose = get_fs_verbosity(),
  ...
)

asegstats2table(
  subjects = NULL,
  inputs = NULL,
  measure = c("volume", "mean", "std"),
  ...
)

aparcstats2table(
  subjects,
  hemi = c("lh", "rh"),
  measure = c("area", "volume", "thickness", "thicknessstd", "meancurv", "gauscurv",
    "foldind", "curvind"),
  parc = c("aparc", "aparc.a2009s"),
  opts = "",
  ...
)

aparcstats2table.help(...)

asegstats2table.help(...)

Arguments

Value

A character string: the path to the output file with its delimiter stored as an attribute.

Functions

• asegstats2table(): Converts subcortical segmentation statistics into tabular format by calling the FreeSurfer asegstats2table command. ... is passed to stats2table for additional options.

• aparcstats2table(): Converts cortical parcellation statistics into tabular format by calling the FreeSurfer aparcstats2table command. ... is passed to stats2table for additional options.'

• aparcstats2table.help(): Display FreeSurfer help for aparcstats2table

• asegstats2table.help(): Display FreeSurfer help for asegstats2table

Examples

## Not run: 
# Example for asegstats2table
outfile_aseg <- asegstats2table(
    subjects = "bert",
    measure = "mean",
    delim = "tab"
)
print(outfile_aseg)

# Example for aparcstats2table
outfile_aparc <- aparcstats2table(
  subjects = "bert",
  hemi = "lh",
  measure = "thickness",
  delim = "tab",
  opts = "--etiv --scale=1.0"
)
print(outfile_aparc)

## End(Not run)

Description

This function calls mri_convert to convert a measure from surfaces to an ASCII file and reads it in.

Usage

surf_convert(file, outfile = NULL, ...)

Arguments

Value

data.frame

Examples

## Not run: 
   fname = file.path(fs_subj_dir(), "bert", "surf", "lh.thickness")
   out = surf_convert(fname)

## End(Not run)

Description

Reads in a surface file from Freesurfer and converts it to a Wavefront OBJ file

Usage

surface_to_obj(infile, outfile = NULL, ...)

Arguments

Value

Character filename of output file

Examples

## Not run: 
infile = file.path(fs_subj_dir(),
                   "bert", "surf", "rh.pial")
res = surface_to_obj(infile = infile)

## End(Not run)

Description

Reads in a surface file from Freesurfer and converts it into triangles

Usage

surface_to_triangles(infile, ...)

Arguments

Value

Matrix of triangles with the number of rows equal to the number of faces (not the triplets - total faces)

Examples

## Not run: 
infile = file.path(fs_subj_dir(),
                   "bert", "surf", "rh.pial")
right_triangles = surface_to_triangles(infile = infile)
infile = file.path(fs_subj_dir(),
                   "bert", "surf", "lh.pial")
left_triangles = surface_to_triangles(infile = infile)
rgl::open3d()
rgl::triangles3d(right_triangles,
color = rainbow(nrow(right_triangles)))
rgl::triangles3d(left_triangles,
color = rainbow(nrow(left_triangles)))

infile = file.path(fs_subj_dir(),
                   "bert", "surf", "rh.inflated")
right_triangles = surface_to_triangles(infile = infile)
infile = file.path(fs_subj_dir(),
                   "bert", "surf", "lh.inflated")
left_triangles = surface_to_triangles(infile = infile)
rgl::open3d()
rgl::triangles3d(left_triangles,
color = rainbow(nrow(left_triangles)))
rgl::triangles3d(right_triangles,
color = rainbow(nrow(right_triangles)))

## End(Not run)

Description

Create a Temporary File with a Newly Created Directory

Usage

temp_file(tmpdir = tempdir(check = TRUE), ...)

Arguments

Value

Full file path to temporary file with created directory

Description

These functions provide wrappers for FreeSurfer's trac-all command, which performs automated reconstruction of diffusion pathways.

Usage

tracker(
  infile = NULL,
  outdir = NULL,
  subjid = NULL,
  verbose = get_fs_verbosity(),
  opts = ""
)

trac_all(
  infile = NULL,
  outdir = NULL,
  subjid = NULL,
  verbose = get_fs_verbosity(),
  opts = ""
)

trac_prep(infile, outdir = NULL, subjid, verbose = get_fs_verbosity())

trac_bedpost(infile, outdir = NULL, subjid, verbose = get_fs_verbosity())

trac_path(infile, outdir = NULL, subjid, verbose = get_fs_verbosity())

tracker.help(...)

Arguments

Details

All functions call the same underlying FreeSurfer trac-all command with different option flags.

Value

Result of base::system() call, typically exit status (0 = success).

Functions

• tracker(): Low-level wrapper for trac-all command with custom options

• trac_all(): High-level wrapper running complete tractography pipeline

• trac_prep(): Run pre-processing step (step 1: image corrections, registration)

• trac_bedpost(): Run bedpost step (step 2: ball-and-stick model fitting)

• trac_path(): Run pathway reconstruction step (step 3: probabilistic tractography)

• tracker.help(): Display FreeSurfer help for trac-all

See Also

recon_all() for structural reconstruction pipeline

Examples

## Not run: 
# Run full tractography pipeline
trac_all(
  infile = "dwi.nii",
  subjid = "subject01",
  outdir = "~/tractography_output"
)

# Run step-by-step
trac_prep(infile = "dwi.nii", subjid = "subject02")
trac_bedpost(infile = "dwi.nii", subjid = "subject02")
trac_path(infile = "dwi.nii", subjid = "subject02")

# Run with custom options
tracker(
  infile = "dwi.nii",
  subjid = "subject03",
  opts = "-prep -bedp"
)

## End(Not run)