Page Last Updated: November 13, 2025

HBCD MR Quality Control Procedures🔗

Raw MR Data QC🔗

Raw MR QC includes automated and manual checks to evaluate unprocessed MRI data. Raw data QC is performed to detect acquisition errors, image artifacts, or corrupted files early so that problematic scans are excluded from the released raw BIDS dataset and downstream processing.

Location of Raw MR QC metrics in Release Data: Scans TSV Files (*_scans.tsv) â–¸

The session folder sub-{ID}_ses-{V0X}_scans.tsv files include the following QC metrics. Note: x values for bx below include 0, 500, 1000, 2000, and 3000

Field Description
filename Relative paths to files
acq_time Acquisition time
loris_qc_status Pass/Fail mapping from UCSD QC JSON file
loris_selected Whether file is selected for further processing
site Site where the session data was collected
age Age (in years) of participant at time of the scan
age_adjusted Adjusted age (in days) based on the EDD
head_size Head size at the time of the scan
nrevManual QC Metric Number of manual QC reviewers
revdispManual QC Metric If there was disparity/disagreement between reviewers
QCManual QC Metric Overall manual QC score: 1 (pass), 0 (fail)
Note: scans with only automated, and not manual, QC performed have a QC field value of 1
notesManual QC Metric Optional notes from manual QC review
QU_motion Manual QC Metric Qualitative manual QC score for motion (sMRI & qMRI)
QU_sus Manual QC Metric Qualitative manual QC score for susceptibility artifact (dMRI, fMRI, & field maps)
QU_cutoff Manual QC Metric Qualitative manual QC score for FOV cutoff artifact (dMRI, fMRI, & field maps)
QU_line Manual QC Metric Qualitative manual QC score for line artifact (dMRI, fMRI, & field maps)
line_nframes Number of frames with line artifacts
line_<max|mean>_score Maximum/Mean line artifact score across frames
line_<max|mean>_count Maximum/Mean line artifact count across frames
cutoff Sum of dorsal and ventral cutoff scores
<dorsal|ventral>_cutoff Dorsal/Ventral cutoff score
brain_<mean|std> Mean/Standard deviation of image intensity within brain mask
brain_SNR Signal-to-noise ratio (mean/stdev) of image intensity within brain mask
brain_<min|max|median> Minimum/Maximum/Median image intensity within brain mask
brain_tSNR_<mean|median|std> Mean/Median/Standard deviation of temporal SNR in brain mask
mean_<motion|trans|rot> Average framewise displacement/translation/rotation (mm)
max_d<x|y|z> Maximum absolute x/y/z translation (mm)
max_r<x|y|z> Maximum absolute x/y/z rotation (mm)
subthresh_<02|03|04> Number of seconds with framewise displacement less than 0.2/0.3/0.4 mm
aqc_motion Automated QC motion score for sMRI
nreps Number of repetitions / frames
brainvol Volume of brain mask (mm^3)
fwhm_x Full width half max spatial smoothness in x-axis (left-right)
b0_<median|std> Median/Standard Deviation b=0 intensity in brain mask
DTerr_<mean|median|std> Mean/Median/STD across frames of RMS residual error relative to RMS signal in brain voxels
Completed Whether series has the expected number of files
NumberOfFiles Number of DICOM files
HBCD_compliant Whether series passes minimal protocol compliance check
AdditionalInfo Notes on classification and protocol compliance
MD_std Standard deviation of mean diffusivity in brain mask
b0_mean Average b=0 intensity in brain mask
FA_std Standard deviation of fractional anisotropy in brain mask
MD_<mean|median> Mean/Median MD (mean diffusivity) in brain mask
max_nbad_frames_per_slice Maximum number of outlier frames per slice
max_nbad_frames_per_frame Maximum number of outlier frames per slice/frame
FA_<mean|median> Mean/Median fractional anisotropy in brain mask
fwhm_z Full width half max spatial smoothness in z-axis (inferior-superior)
nbad_frame_slices Number of outlier frame-slices (dMRI)
nbad_frames Number of frames with outlier slices (dMRI)
nbad_slices Number of slices with outlier frames (dMRI)
fwhm_y Full width half max spatial smoothness in y-axis (anterior-posterior)
qc_status Whether review is pending, complete, or has other status
ngood_frames Number of frames without outlier slices (dMRI)
censor_thresh Threshold used for censoring outlier slices (dMRI)
nframes_bx Number of b=x frames
nbad_frame_slices_bx Number of outlier frame-slices for dMRI b=x frames
nbad_frames_bx Number of frames with outlier slices for dMRI b=x frames
nbad_slices_bx Number of slices with outlier frames for dMRI b=x frames
ngood_frames_bx Number of frames without outlier slices for dMRI b=x frames
FWHMx_bx FWHM spatial smoothness in x-axis (left-right) for b=x frames
FWHMy_bx FWHM spatial smoothness in y-axis (anterior-posterior) for b=x frames
FWHMz_bx FWHM spatial smoothness in z-axis (inferior-superior) for b=x frames
tSNR_bx Median temporal SNR in brain mask for for b=x frames
DTerr_rel_bx Median DTI RMS residual error for b=x frames across voxels relative to within-voxel RMS b=0 signal
RSIerr_rel_bx Median RSI RMS residual error for b=x frames across voxels relative to within-voxel RMS b=0 signal
DTerr_rel Median of DTI RMS residual error for all frames across voxels relative to within-voxel RMS signal
RSIerr_rel Median RSI RMS residual error for all frames across voxels relative to within-voxel RMS signal
NumberOfFilesMissing Number of DICOM files apparently missing (based on gaps in InstanceNumbers)
Num<Head|Neck|Spine>CoilElem Number of head/neck/spine coil elements
brain_nvox_max Number of voxels within brain mask at maximum image intensity
brain_fvox_max Fraction of voxels within brain mask at maximum image intensity
nonbrain_<mean|std|snr> Mean/Standard deviation/Signal-to-noise ratio (mean/stdev) image intensity outside brain mask
NumberOfFilesOrig Number of DICOM files received (before excluding non-image, corrupt, or extra files)
NumberOfFilesExtra Number of extra DICOM files received (non-image, corrupt, or extra files)

Automated QC🔗

Automated QC is performed at the HBCD Data Coordinating Center (HDCC) immediately after data upload. Following protocol compliance and completeness checks, the following automated QC metrics available in the release are calculated:

Automated QC Metrics
Modality QC Procedures
sMRI & qMRI • Estimate motion artifacts using a deep learning model
• Compute signal-to-noise ratio (SNR)
fMRI • Estimate head motion with average FDframewise displacement and data (sec) at FD thresholds of 0.2/0.3/0.4 mm (Power et al., 2012)
• Detect line artifacts and FOV cutoff
• Compute spatial smoothness (FWHM) and temporal SNR (tSNR) after motion correction (Triantafyllou et al., 2005)
dMRI • Estimate head motion (framewise displacement, FD)
• Refine motion estimates via registration to tensor-synthesized images (Hagler et al., 2009)
• Identify dark slices (caused by abrupt head movements) using RMS difference between raw and tensor-fitted data
• Calculate total slices and frames with motion artifacts
• Detect line artifacts and field-of-view (FOV) cutoff
Field Maps Detect line artifacts and field-of-view (FOV) cutoff
All Compute SNR where applicable
Protocol Compliance & Completeness Checks â–¸

Protocol Compliance

  • Extract imaging parameters from DICOM headers
  • Confirm key parameters (e.g., voxel size, TR, orientation) match the expected protocol for each scanner
  • Flag out-of-compliance series for review; sites are contacted if corrective action is needed

Completeness Checks
Completeness checks verify that all expected series are present in each imaging session. For example, dMRI and fMRI require paired EPI field maps for distortion correction. Missing data usually indicate an aborted scan or incomplete data transfer (often resolved by re-sending files). Series included in a valid session include the following:

Modality Required Series
Structural MRI (sMRI) T1w & T2w
Functional MRI (fMRI) Resting state (rs) runs 1 & 2, each accompanied by fieldmaps acquired in AP and PA phase encoding directions
Diffusion MRI (dMRI) Diffusion scans acquired in AP and PA phase encoding directions
Quantitative MRI (qMRI) 3DMagic/QALAS and B1 Map
MR Spectroscopy (MRS) MRS scan and SVS localizer

Manual Review🔗

Automated metrics flag series for manual review using multivariate prediction and Bayesian classifiers. Trained technicians then score artifact severity on a 0–3 scale (0 = none, 1 = mild, 2 = moderate, 3 = severe). Series with severe artifacts (score = 3) are rejected (QC = 0) and excluded from processing. Data are selected from the remaining series based on manual ratings, notes, and automated scores.

Manual QC Metrics
Modality QC Procedures & Scoring
sMRI • Motion artifacts (ripples, blurring), scored 0–3
• Document additional issues (e.g., intensity inhomogeneity, ghostingfaint displaced copy of anatomy due to slices outside FOV)
qMRI • Same artifact scoring (0–3)
• Inspect derived data (parametric maps, ROI analysis, quantitative checks for 3D-QALAS)
dMRI & fMRI, & field maps • Score susceptibility artifacts, FOV cutoff, and horizontal line artifacts (present in the sagittal view)
• Note susceptibility artifacts, including signal dropoutcommon in posterior occipital cortex of infant fMRI data acquired in PA phase encoding direction, signal bunching, and warping
MRS Visual inspection and overall QC only of SVS localizer (used to define spectroscopy ROI)

BrainSwipes🔗

Manual visual inspection remains the gold standard for detecting artifacts in structural, functional (e.g., XCP-D), and diffusion (e.g., QSIPrep) derivatives. To support this, derivative visual reports are integrated into BrainSwipes, a gamified, crowdsourced QC platform built on the open-source Swipes For Science framework. BrainSwipes provides an intuitive interface for large-scale studies, guiding users through a short tutorial before they evaluate images and classify them as pass or fail.

BrainSwipes QC Procedures â–¸

Example quality assessment of surface delineation in BrainSwipes

Example quality assessment of surface delineation on BrainSwipes platform (displaying brain in axial plane at level of basal ganglia/putamen).

Surface Delineation:

For structural QA, swipers are presented with image slices in coronal, axial, and sagittal planes to assess the accuracy of T1w and T2w surface delineations in differentiating gray and white matter. Images are derived from XCP-D visual reports.

Atlas Registration:

In addition to surface delineation, structural QA also includes atlas registration quality, evaluated by overlaying delineations of the subject’s image onto the atlas, and vice versa. Swipes display nine T1w slices for visual inspection, with three slices per anatomical plane. Quality is assessed based on the alignment of the outer boundaries of the overlaid contours with those of the underlying image, ensuring minimal gaps or misalignments. Images are derived from XCP-D visual reports.

Functional Registration:

Functional registration is evaluated by overlaying outlines of functional images onto structural images and vice versa. Swipes display nine slices of the same functional image for visual inspection, with three slices per anatomical plane. Quality is assessed similarly to structural atlas registration, focusing on the alignment of the overlaid contours. Additional evaluation includes checking for artifacts such as signal dropout. Images are derived from XCP-D visual reports.

Diffusion Direction Encoding (to be included in future release):

Swipes display GIFs of full-resolution T2w images as a grayscale background, with the "Direction Encoded Color" (DEC) map overlaid. These GIFs sweep through a portion of the brain across the three anatomical planes. High-quality processed DWI images exhibit bands of color that closely follow the folds and contours of the grayscale background. These visuals are derived from the QSIPrep report.

Each visual report for a given modality is independently reviewed and rated as a pass or fail, which in the outputs are scored as values of 1 and 0 respectively. BrainSwipes generates a summary of these results that includes the average score as well as number of reviewers for each visual report of each modality.

BrainSwipes QC results are provided as tabulated data: img_brainswipes_xcpd-T2w and img_brainswipes_xcpd-bold. QC scores range from 0 (Fail) to 1 (Pass), averaged across reviewers. For example, a score of 0.6 indicates that 60% of reviewers rated the image as a pass. The data includes overall average QC scores and average number of reviewers for each session-level T2w and BOLD run (summarizing across all visual reports for a given run) as well as the average QC and number of reviewers for each visual report. A Python helper function is provided below to load a BrainSwipes TSV file into a Pandas DataFrame and filter runs with average QC scores above a user-specified threshold:

Python Helper Function â–¸

import pandas as pd

def read_and_filter_tsv(file_path, threshold):
    """
    Reads an HBCD BrainSwipes TSV file into a DataFrame and keeps subject rows where 
    the overall average QC value is greater than or equal to the specified threshold.

    Parameters:
    - file_path (str): Path to the .tsv file
    - threshold (float): Threshold value for filtering

    Returns:
    - pd.DataFrame: Filtered DataFrame
    """
    df = pd.read_csv(file_path, sep='\t')

    fourth_col = df.columns[3]
    return df[df[fourth_col] >= threshold]

# Example usage:
# Filter to keep only subjects with an average structural QC of at least 0.6 (60% pass rate)
filtered_df = read_and_filter_tsv("img_brainswipes_xcpd_T2w.tsv", 0.6)
print(filtered_df.head())

References🔗

Dean III, D. C., Tisdall, M. D., Wisnowski, J. L., Feczko, E., Gagoski, B., Alexander, A. L., ... & HBCD MRI Working Group. (2024). Quantifying brain development in the HEALthy Brain and Child Development (HBCD) Study: The magnetic resonance imaging and spectroscopy protocol. Developmental Cognitive Neuroscience, 70, 101452. https://doi.org/10.1016/j.dcn.2024.101452

Gard, A. M., Hyde, L. W., Heeringa, S. G., West, B. T., & Mitchell, C. (2023). Why weight? Analytic approaches for large-scale population neuroscience data. Developmental Cognitive Neuroscience, 59, 101196. https://doi.org/10.1016/j.dcn.2023.101196

Hagler, D. J., Jr, Ahmadi, M. E., Kuperman, J., Holland, D., McDonald, C. R., Halgren, E., & Dale, A. M. (2009). Automated white-matter tractography using a probabilistic diffusion tensor atlas: Application to temporal lobe epilepsy. Human Brain Mapping, 30(5), 1535–1547. https://doi.org/10.1002/hbm.20619

Power, J. D., Barnes, K. A., Snyder, A. Z., Schlaggar, B. L., & Petersen, S. E. (2012). Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. NeuroImage, 59(3), 2142–2154. https://doi.org/10.1016/j.neuroimage.2011.10.018

Triantafyllou, C., Hoge, R. D., Krueger, G., Wiggins, C. J., Potthast, A., Wiggins, G. C., & Wald, L. L. (2005). Comparison of physiological noise at 1.5 T, 3 T and 7 T and optimization of fMRI acquisition parameters. NeuroImage, 26(1), 243–250. https://doi.org/10.1016/j.neuroimage.2005.01.007