Automatic, whole-brain segmentation of cortical layers for the high-resolution anatomical data at 7
Poster No:
1493
Submission Type:
Abstract Submission
Authors:
Roman Belenya1, Gabriel Castrillón2, Valentin Riedl2
Institutions:
1TUM Klinikum rechts der Isar, Munich, Bavaria, 2Friedrich-Alexander University, Erlangen, Bavaria
First Author:
Co-Author(s):
Introduction:
Developments of high-resolution MR imaging at 7 Tesla have allowed researchers to study the brain in unprecedented detail. The increased magnetic field improves the signal-to-noise ratio and the resolution of the images. However, 7T imaging comes with a set of challenges. Among them are magnetic field inhomogeneities and imaging artefacts. They present a common challenge for data preprocessing. Over the years, robust and reproducible preprocessing methods have been well-optimized for standard 3T imaging. However, there are currently no standard optimized pipelines for high-resolution brain images. In this work, we develop a preprocessing pipeline for the 7T hires anatomical data. We focus on performing as few processing steps as possible while maximizing the performance of automatic whole-brain segmentation. Specifically, our end-goal is to maximize the accuracy and reliability of cortical layer segmentation for the analyses of layer fMRI data.
Methods:
We processed the anatomical data of 14 participants from the 7T Siemens MAGNETOM Terra.X machine. We collected the 3D T1-weighted MP2RAGE as well as T2-weighted TSE images. The preprocessing pipeline is similar for both modalities:
Pre-freesurfer
1. Image denoising using the DenoiseImage program in ANTs.
1a. Gibbs ringing artefacts removal with the mrdegibbs program implemented in mrtrix. We run this step only for the T2w images because they are affected by the artefact the most.
2. Bias field correction using SPM tissue-type segmentation.
3. Brain extraction with HD-BET.
Freesurfer
4. Brain surface reconstruction with freesurfer recon-all script. We run the standard processing using the high-resolution (-cm) option and an additional T2w contrast to improve the tissue type segmentation. We entered the combined brain mask from Step 3 using the expert file.
Post-Freesurfer
5. To convert brain surfaces generated by the freesurfer to a voxel space rim file required by LayNii, we use the surf-laynii script by Sriranga Kashyap. The script first expands the grey matter and shrinks the white matter surfaces using the mri_expand tool. The mri_fill tool fills the original, expanded and shrunk surfaces, resulting in a set of masks. Subtracting the smaller from a larger mask estimates the boundary between the white matter and the grey matter and the boundary between the CSF and the white matter in voxel space.
6. We use these boundaries as the input to the LayNii LN2_LAYERS program, which estimates cortical layers using an equivolumetric approach.
Pre-freesurfer
1. Image denoising using the DenoiseImage program in ANTs.
1a. Gibbs ringing artefacts removal with the mrdegibbs program implemented in mrtrix. We run this step only for the T2w images because they are affected by the artefact the most.
2. Bias field correction using SPM tissue-type segmentation.
3. Brain extraction with HD-BET.
Freesurfer
4. Brain surface reconstruction with freesurfer recon-all script. We run the standard processing using the high-resolution (-cm) option and an additional T2w contrast to improve the tissue type segmentation. We entered the combined brain mask from Step 3 using the expert file.
Post-Freesurfer
5. To convert brain surfaces generated by the freesurfer to a voxel space rim file required by LayNii, we use the surf-laynii script by Sriranga Kashyap. The script first expands the grey matter and shrinks the white matter surfaces using the mri_expand tool. The mri_fill tool fills the original, expanded and shrunk surfaces, resulting in a set of masks. Subtracting the smaller from a larger mask estimates the boundary between the white matter and the grey matter and the boundary between the CSF and the white matter in voxel space.
6. We use these boundaries as the input to the LayNii LN2_LAYERS program, which estimates cortical layers using an equivolumetric approach.
Results:
We run the pipeline on data from 14 subjects. On average, the pre-freesurfer pipeline doubles the spatial SNR and CNR within the grey and white matter masks for T1w and T2w images. Combining the preprocessed contrasts allows for more accurate placement of brain cortical surfaces and estimation of cortical layers. We demonstrate an example of processing steps and their effect on the images in Figure 1.
Conclusions:
The accuracy of brain segmentation is crucial in analyzing neuroimaging data. While manual segmentations can be most accurate, they are time-consuming, especially on the whole-brain level. The existing automatic segmentation tools are usually optimized for images from 3T scanners and cannot be used out of the box for high-resolution 7T data. Thus, it is essential to configure the preprocessing pipelines to maximize the accuracy of automatic whole-brain segmentation. Here, we describe the processing pipeline that optimized the performance of cortical layer segmentation for our anatomical data. Accurate automatic layer segmentation can dramatically increase the accessibility of high-resolution data analyses. This work will assist researchers in developing analysis methods for the high-resolution neuroimaging data.
Modeling and Analysis Methods:
Image Registration and Computational Anatomy 1
Methods Development 2
Segmentation and Parcellation
Keywords:
Acquisition
Cortical Layers
HIGH FIELD MR
MRI
Segmentation
STRUCTURAL MRI
1|2Indicates the priority used for review
Abstract Information
By submitting your proposal, you grant permission for the Organization for Human Brain Mapping (OHBM) to distribute your work in any format, including video, audio print and electronic text through OHBM OnDemand, social media channels, the OHBM website, or other electronic publications and media.
The Open Science Special Interest Group (OSSIG) is introducing a reproducibility challenge for OHBM 2025. This new initiative aims to enhance the reproducibility of scientific results and foster collaborations between labs. Teams will consist of a “source” party and a “reproducing” party, and will be evaluated on the success of their replication, the openness of the source work, and additional deliverables. Click here for more information. Propose your OHBM abstract(s) as source work for future OHBM meetings by selecting one of the following options:
Please indicate below if your study was a "resting state" or "task-activation” study.
Healthy subjects only or patients (note that patient studies may also involve healthy subjects):
Was this research conducted in the United States?
Were any human subjects research approved by the relevant Institutional Review Board or ethics panel? NOTE: Any human subjects studies without IRB approval will be automatically rejected.
Were any animal research approved by the relevant IACUC or other animal research panel? NOTE: Any animal studies without IACUC approval will be automatically rejected.
Please indicate which methods were used in your research:
For human MRI, what field strength scanner do you use?
Which processing packages did you use for your study?
Provide references using APA citation style.
Isensee F. (2019). Automated brain extraction of multi-sequence MRI using artificial neural networks. Human Brain Mapping, 1–13
Sriranga Kashyap. Surf-laynii. https://github.com/srikash/surf_laynii
Tournier, J.-D. (2019). MRtrix3: A fast, flexible and open software framework for medical image processing and visualisation. NeuroImage, 202, 116137
Tustison, N. J. (2014). Large-scale evaluation of ANTs and FreeSurfer cortical thickness measurements. NeuroImage, 99, 166–179