Assessment of intraoperative diffusion EPI distortion and its impact on estimation of supratentorial white matter tract positions in pediatric epilepsy surgery

PyHySCO: GPU-enabled susceptibility artifact distortion correction in seconds

Over the past decade, reversed gradient polarity (RGP) methods have become a popular approach for correcting susceptibility artifacts in echo-planar imaging (EPI). Although several post-processing tools for RGP are available, their implementations do not fully leverage recent hardware, algorithmic, and computational advances, leading to correction times of several minutes per image volume. To enable 3D RGP correction in seconds, we introduce PyTorch Hyperelastic Susceptibility Correction (PyHySCO), a user-friendly EPI distortion correction tool implemented in PyTorch that enables multi-threading and efficient use of graphics processing units (GPUs). PyHySCO uses a time-tested physical distortion model and mathematical formulation and is, therefore, reliable without training. An algorithmic improvement in PyHySCO is its use of the one-dimensional distortion correction method by Chang and Fitzpatrick to initialize the non-linear optimization. PyHySCO is published under the GNU public license and can be used from the command line or its Python interface. Our extensive numerical validation using 3T and 7T data from the Human Connectome Project suggests that PyHySCO can achieve accuracy comparable to that of leading RGP tools at a fraction of the cost. We also validate the new initialization scheme, compare different optimization algorithms, and test the algorithm on different hardware and arithmetic precisions.

Fibre orientation atlas guided rapid segmentation of white matter tracts

Fibre tract delineation from diffusion magnetic resonance imaging (MRI) is a valuable clinical tool for neurosurgical planning and navigation, as well as in research neuroimaging pipelines. Several popular methods are used for this task, each with different strengths and weaknesses making them more or less suited to different contexts. For neurosurgical imaging, priorities include ease of use, computational efficiency, robustness to pathology and ability to generalise to new tracts of interest. Many existing methods use streamline tractography, which may require expert neuroimaging operators for setting parameters and delineating anatomical regions of interest, or suffer from as a lack of generalisability to clinical scans involving deforming tumours and other pathologies. More recently, data-driven approaches including deep-learning segmentation models and streamline clustering methods have improved reproducibility and automation, although they can require large amounts of training data and/or computationally intensive image processing at the point of application. We describe an atlas-based direct tract mapping technique called 'tractfinder', utilising tract-specific location and orientation priors. Our aim was to develop a clinically practical method avoiding streamline tractography at the point of application while utilising prior anatomical knowledge derived from only 10-20 training samples. Requiring few training samples allows emphasis to be placed on producing high quality, neuro-anatomically accurate training data, and enables rapid adaptation to new tracts of interest. Avoiding streamline tractography at the point of application reduces computational time, false positives and vulnerabilities to pathology such as tumour deformations or oedema. Carefully filtered training streamlines and track orientation distribution mapping are used to construct tract specific orientation and spatial probability atlases in standard space. Atlases are then transformed to target subject space using affine registration and compared with the subject's voxel-wise fibre orientation distribution data using a mathematical measure of distribution overlap, resulting in a map of the tract's likely spatial distribution. This work includes extensive performance evaluation and comparison with benchmark techniques, including streamline tractography and the deep-learning method TractSeg, in two publicly available healthy diffusion MRI datasets (from TractoInferno and the Human Connectome Project) in addition to a clinical dataset comprising paediatric and adult brain tumour scans. Tract segmentation results display high agreement with established techniques while requiring less than 3 min on average when applied to a new subject. Results also display higher robustness than compared methods when faced with clinical scans featuring brain tumours and resections. As well as describing and evaluating a novel proposed tract delineation technique, this work continues the discussion on the challenges surrounding the white matter segmentation task, including issues of anatomical definitions and the use of quantitative segmentation comparison metrics.

Distortion correction of functional MRI without reverse phase encoding scans or field maps

Functional magnetic resonance images (fMRI) acquired using echo planar sequences typically suffer from spatial distortions due to susceptibility induced off-resonance fields, which may cause geometric mismatch with structural images and affect subsequent quantification and localization of brain function. State-of-the art distortion correction methods (for example, using FSL's topup or AFNI's 3dQwarp algorithms) require the collection of additional scans - either field maps or images with reverse phase encoding directions (i.e., blip-up/blip-down acquisitions) - to estimate and correct distortions. However, not all imaging protocols acquire these additional data and thus cannot take advantage of these post-acquisition corrections. In this study, we aim to enable state-of-the art processing of historical or limited datasets that do not include specific sequences for distortion correction by using only the acquired functional data and a single commonly acquired structural image. To achieve this, we synthesize an undistorted image with contrast similar to the fMRI data and use the non-distorted synthetic image as an anatomical target for distortion correction. We evaluate the efficacy of this approach, named SynBOLD-DisCo (Synthetic BOLD contrast for Distortion Correction), and show that this distortion correction process yields fMRI data that are geometrically similar to non-distorted structural images, with distortion correction virtually equivalent to acquisitions that do contain both blip-up/blip-down images. Our method is available as a Singularity container, source code, and an executable trained model to facilitate evaluation and integration into existing fMRI preprocessing pipelines.

Quantitative Validation of the Correlation Between Optimized Pyramidal Tract Delineation After Brain Shift Compensation and Direct Electrical Subcortical Stimulation During Brain Tumor Surgery

It remains unclear whether tractography of pyramidal tracts is correlated with the intraoperative direct electrical subcortical stimulation (DESS), and brain shift further complicates the issue. The objective of this research is to quantitatively verify the correlation between optimized tractography (OT) of pyramidal tracts after brain shift compensation and DESS during brain tumor surgery. OT was performed for 20 patients with lesions in proximity to the pyramidal tracts based on preoperative diffusion-weighted magnetic resonance imaging. During surgery, tumor resection was guided by DESS. A total of 168 positive stimulation points and their corresponding stimulation intensity thresholds were recorded. Using the brain shift compensation algorithm based on hierarchical B-spline grids combined with a Gaussian resolution pyramid, we warped the preoperative pyramidal tract models and used receiver operating characteristic (ROC) curves to investigate the reliability of our brain shift compensation method based on anatomic landmarks. Additionally, the minimum distance between the DESS points and warped OT (wOT) model was measured and correlated with DESS intensity threshold. Brain shift compensation was achieved in all cases, and the area under the ROC curve was 0.96 in the registration accuracy analysis. The minimum distance between the DESS points and the wOT model was found to have a significantly high correlation with the DESS stimulation intensity threshold (r = 0.87, P < 0.001), with a linear regression coefficient of 0.96. Our OT method can provide comprehensive and accurate visualization of the pyramidal tracts for neurosurgical navigation and was quantitatively verified by intraoperative DESS after brain shift compensation.

SynBOLD-DisCo: Synthetic BOLD images for distortion correction of fMRI without additional calibration scans

The blood oxygen level dependent (BOLD) signal from functional magnetic resonance imaging (fMRI) is a noninvasive technique that has been widely used in research to study brain function. However, fMRI suffers from susceptibility-induced off resonance fields which may cause geometric distortions and mismatches with anatomical images. State-of-the-art correction methods require acquiring reverse phase encoded images or additional field maps to enable distortion correction. However, not all imaging protocols include these additional scans and thus cannot take advantage of these susceptibility correction capabilities. As such, in this study we aim to enable state-of-the-art distortion correction with FSL's topup algorithm of historical and/or limited fMRI data that include only a structural image and single phase encoded fMRI. To do this, we use 3D U-net models to synthesize undistorted fMRI BOLD contrast images from the structural image and use this undistorted synthetic image as an anatomical target for distortion correction with topup. We evaluate the efficacy of this approach, named SynBOLD-DisCo (synthetic BOLD images for distortion correction), and show that BOLD images corrected using our approach are geometrically more similar to structural images than the distorted BOLD data and are practically equivalent to state-of-the-art correction methods which require reverse phase encoded data. Future directions include additional validation studies, integration with other preprocessing operations, retraining with broader pathologies, and investigating the effects of spin echo versus gradient echo images for training and distortion correction. In summary, we demonstrate SynBOLD-DisCo corrects distortion of fMRI when reverse phase encoding scans or field maps are not available.