Journal Papers

S

Kurugol, Sila, Moti Freiman, Onur Afacan, Jeannette Perez-Rossello, Michael Callahan, and Simon Warfield. 2016. “Spatially-Constrained Probability Distribution Model of Incoherent Motion (SPIM) for Abdominal Diffusion-Weighted MRI”. Med Image Anal 32: 173-83. https://doi.org/10.1016/j.media.2016.03.009.
Quantitative diffusion-weighted MR imaging (DW-MRI) of the body enables characterization of the tissue microenvironment by measuring variations in the mobility of water molecules. The diffusion signal decay model parameters are increasingly used to evaluate various diseases of abdominal organs such as the liver and spleen. However, previous signal decay models (i.e., mono-exponential, bi-exponential intra-voxel incoherent motion (IVIM) and stretched exponential models) only provide insight into the average of the distribution of the signal decay rather than explicitly describe the entire range of diffusion scales. In this work, we propose a probability distribution model of incoherent motion that uses a mixture of Gamma distributions to fully characterize the multi-scale nature of diffusion within a voxel. Further, we improve the robustness of the distribution parameter estimates by integrating spatial homogeneity prior into the probability distribution model of incoherent motion (SPIM) and by using the fusion bootstrap solver (FBM) to estimate the model parameters. We evaluated the improvement in quantitative DW-MRI analysis achieved with the SPIM model in terms of accuracy, precision and reproducibility of parameter estimation in both simulated data and in 68 abdominal in-vivo DW-MRIs. Our results show that the SPIM model not only substantially reduced parameter estimation errors by up to 26%; it also significantly improved the robustness of the parameter estimates (paired Student's t-test, p < 0.0001) by reducing the coefficient of variation (CV) of estimated parameters compared to those produced by previous models. In addition, the SPIM model improves the parameter estimates reproducibility for both intra- (up to 47%) and inter-session (up to 30%) estimates compared to those generated by previous models. Thus, the SPIM model has the potential to improve accuracy, precision and robustness of quantitative abdominal DW-MRI analysis for clinical applications.
Taimouri, Vahid, Onur Afacan, Jeannette Perez-Rossello, Michael Callahan, Robert Mulkern, Simon Warfield, and Moti Freiman. (2015) 2015. “Spatially Constrained Incoherent Motion Method Improves Diffusion-Weighted MRI Signal Decay Analysis in the Liver and Spleen”. Med Phys 42 (4): 1895-903. https://doi.org/10.1118/1.4915495.
PURPOSE: To evaluate the effect of the spatially constrained incoherent motion (SCIM) method on improving the precision and robustness of fast and slow diffusion parameter estimates from diffusion-weighted MRI in liver and spleen in comparison to the independent voxel-wise intravoxel incoherent motion (IVIM) model. METHODS: We collected diffusion-weighted MRI (DW-MRI) data of 29 subjects (5 healthy subjects and 24 patients with Crohn's disease in the ileum). We evaluated parameters estimates' robustness against different combinations of b-values (i.e., 4 b-values and 7 b-values) by comparing the variance of the estimates obtained with the SCIM and the independent voxel-wise IVIM model. We also evaluated the improvement in the precision of parameter estimates by comparing the coefficient of variation (CV) of the SCIM parameter estimates to that of the IVIM. RESULTS: The SCIM method was more robust compared to IVIM (up to 70% in liver and spleen) for different combinations of b-values. Also, the CV values of the parameter estimations using the SCIM method were significantly lower compared to repeated acquisition and signal averaging estimated using IVIM, especially for the fast diffusion parameter in liver (CVIV IM = 46.61 ± 11.22, CVSCIM = 16.85 ± 2.160, p < 0.001) and spleen (CVIV IM = 95.15 ± 19.82, CVSCIM = 52.55 ± 1.91, p < 0.001). CONCLUSIONS: The SCIM method characterizes fast and slow diffusion more precisely compared to the independent voxel-wise IVIM model fitting in the liver and spleen.
Sui, Yao, Onur Afacan, Ali Gholipour, and Simon Warfield. 2020. “SLIMM: Slice Localization Integrated MRI Monitoring”. Neuroimage 223: 117280. https://doi.org/10.1016/j.neuroimage.2020.117280.
Functional MRI (fMRI) is extremely challenging to perform in subjects who move because subject motion disrupts blood oxygenation level dependent (BOLD) signal measurement. It has become common to use retrospective framewise motion detection and censoring in fMRI studies to eliminate artifacts arising from motion. Data censoring results in significant loss of data and statistical power unless the data acquisition is extended to acquire more data not corrupted by motion. Acquiring more data than is necessary leads to longer than necessary scan duration, which is more expensive and may lead to additional subject non-compliance. Therefore, it is well established that real-time prospective motion monitoring is crucial to ensure data quality and reduce imaging costs. In addition, real-time monitoring of motion allows for feedback to the operator and the subject during the acquisition, to enable intervention to reduce the subject motion. The most widely used form of motion monitoring for fMRI is based on volume-to-volume registration (VVR), which quantifies motion as the misalignment between subsequent volumes. However, motion is not constrained to occur only at the boundaries of volume acquisition, but instead may occur at any time. Consequently, each slice of an fMRI acquisition may be displaced by motion, and assessment of whole volume to volume motion may be insensitive to both intra-volume and inter-volume motion that is revealed by displacement of the slices. We developed the first slice-by-slice self-navigated motion monitoring system for fMRI by developing a real-time slice-to-volume registration (SVR) algorithm. Our real-time SVR algorithm, which is the core of the system, uses a local image patch-based matching criterion along with a Levenberg-Marquardt optimizer, all accelerated via symmetric multi-processing, with interleaved and simultaneous multi-slice acquisition schemes. Extensive experimental results on real motion data demonstrated that our fast motion monitoring system, named Slice Localization Integrated MRI Monitoring (SLIMM), provides more accurate motion measurements than a VVR based approach. Therefore, SLIMM offers improved online motion monitoring which is particularly important in fMRI for challenging patient populations. Real-time motion monitoring is crucial for online data quality control and assurance, for enabling feedback to the subject and the operator to act to mitigate motion, and in adaptive acquisition strategies that aim to ensure enough data of sufficient quality is acquired without acquiring excess data.
Afacan, Onur, Scott Hoge, Tess Wallace, Ali Gholipour, Sila Kurugol, and Simon Warfield. 2020. “Simultaneous Motion and Distortion Correction Using Dual-Echo Diffusion-Weighted MRI”. J Neuroimaging 30 (3): 276-85. https://doi.org/10.1111/jon.12708.
BACKGROUND AND PURPOSE: Geometric distortions resulting from large pose changes reduce the accuracy of motion measurements and interfere with the ability to generate artifact-free information. Our goal is to develop an algorithm and pulse sequence to enable motion-compensated, geometric distortion compensated diffusion-weighted MRI, and to evaluate its efficacy in correcting for the field inhomogeneity and position changes, induced by large and frequent head motions. METHODS: Dual echo planar imaging (EPI) with a blip-reversed phase encoding distortion correction technique was evaluated in five volunteers in two separate experiments and compared with static field map distortion correction. In the first experiment, dual-echo EPI images were acquired in two head positions designed to induce a large field inhomogeneity change. A field map and a distortion-free structural image were acquired at each position to assess the ability of dual-echo EPI to generate reliable field maps and enable geometric distortion correction in both positions. In the second experiment, volunteers were asked to move to multiple random positions during a diffusion scan. Images were reconstructed using the dual-echo correction and a slice-to-volume registration (SVR) registration algorithm. The accuracy of SVR motion estimates was compared to externally measured ground truth motion parameters. RESULTS: Our results show that dual-echo EPI can produce slice-level field maps with comparable quality to field maps generated by the reference gold standard method. We also show that slice-level distortion correction improves the accuracy of SVR algorithms as slices acquired at different orientations have different levels of distortion, which can create errors in the registration process. CONCLUSIONS: Dual-echo acquisitions with blip-reversed phase encoding can be used to generate slice-level distortion-free images, which is critical for motion-robust slice to volume registration. The distortion corrected images not only result in better motion estimates, but they also enable a more accurate final diffusion image reconstruction.
Vasylechko, Serge Didenko, Simon Warfield, Onur Afacan, and Sila Kurugol. 2022. “Self-Supervised IVIM DWI Parameter Estimation With a Physics Based Forward Model”. Magn Reson Med 87 (2): 904-14. https://doi.org/10.1002/mrm.28989.
PURPOSE: To assess the robustness and repeatability of intravoxel incoherent motion model (IVIM) parameter estimation for the diffusion-weighted MRI in the abdominal organs under the constraints of noisy diffusion signal using a novel neural network method. METHODS: Clinically acquired abdominal scans of Crohn's disease patients were retrospectively analyzed with regions segmented in the kidney cortex, spleen, liver, and bowel. A novel IVIM parameter fitting method based on the principle of a physics guided self-supervised convolutional neural network that does not require reference parameter estimates for training was compared to a conventional non-linear least squares (NNLS) algorithm, and a voxelwise trained artificial neural network (ANN). RESULTS: Results showed substantial increase in parameter robustness to the noise corrupted signal. In an intra-session repeatability experiment, the proposed method showed reduced coefficient of variation (CoV) over multiple acquisitions in comparison to conventional NLLS method and comparable performance to ANN. The use of D and f estimates from the proposed method led to the smallest misclassification error in linear discriminant analysis for characterization between normal and abnormal Crohn's disease bowel tissue. The fitting of D ∗ parameter remains to be challenging. CONCLUSION: The proposed method yields robust estimates of D and f IVIM parameters under the constraints of noisy diffusion signal. This indicates a potential for the use of the proposed method in conjunction with accelerated DW-MRI acquisition strategies, which would typically result in lower signal to noise ratio.
Sui, Yao, Onur Afacan, Camilo Jaimes, Ali Gholipour, and Simon Warfield. 2022. “Scan-Specific Generative Neural Network for MRI Super-Resolution Reconstruction”. IEEE Trans Med Imaging 41 (6): 1383-99. https://doi.org/10.1109/TMI.2022.3142610.
The interpretation and analysis of Magnetic resonance imaging (MRI) benefit from high spatial resolution. Unfortunately, direct acquisition of high spatial resolution MRI is time-consuming and costly, which increases the potential for motion artifact, and suffers from reduced signal-to-noise ratio (SNR). Super-resolution reconstruction (SRR) is one of the most widely used methods in MRI since it allows for the trade-off between high spatial resolution, high SNR, and reduced scan times. Deep learning has emerged for improved SRR as compared to conventional methods. However, current deep learning-based SRR methods require large-scale training datasets of high-resolution images, which are practically difficult to obtain at a suitable SNR. We sought to develop a methodology that allows for dataset-free deep learning-based SRR, through which to construct images with higher spatial resolution and of higher SNR than can be practically obtained by direct Fourier encoding. We developed a dataset-free learning method that leverages a generative neural network trained for each specific scan or set of scans, which in turn, allows for SRR tailored to the individual patient. With the SRR from three short duration scans, we achieved high quality brain MRI at an isotropic spatial resolution of 0.125 cubic mm with six minutes of imaging time for T2 contrast and an average increase of 7.2 dB (34.2%) in SNR to these short duration scans. Motion compensation was achieved by aligning the three short duration scans together. We assessed our technique on simulated MRI data and clinical data acquired from 15 subjects. Extensive experimental results demonstrate that our approach achieved superior results to state-of-the-art methods, while in parallel, performed at reduced cost as scans delivered with direct high-resolution acquisition.
Khawaled, Samah, Onur Afacan, Simon K Warfield, and Moti Freiman. (2025) 2025. “A Self-Attention Model for Robust Rigid Slice-to-Volume Registration of Functional MRI.”. Computerized Medical Imaging and Graphics : The Official Journal of the Computerized Medical Imaging Society 125: 102643. https://doi.org/10.1016/j.compmedimag.2025.102643.

Functional Magnetic Resonance Imaging (fMRI) is vital in neuroscience, enabling investigations into brain disorders, treatment monitoring, and brain function mapping. However, head motion during fMRI scans, occurring between shots of slice acquisition, can result in distortion, biased analyses, and increased costs due to the need for scan repetitions. Therefore, retrospective slice-level motion correction through slice-to-volume registration (SVR) is crucial. Previous studies have utilized deep learning (DL) based models to address the SVR task; however, they overlooked the uncertainty stemming from the input stack of slices and did not assign weighting or scoring to each slice. Treating all slices equally ignores the variability in their relevance, leading to suboptimal predictions. In this work, we introduce an end-to-end SVR model for aligning 2D fMRI slices with a 3D reference volume, incorporating a self-attention mechanism to enhance robustness against input data variations and uncertainties. Our SVR model utilizes independent slice and volume encoders and a self-attention module to assign pixel-wise scores for each slice. We used the publicly available Healthy Brain Network (HBN) dataset. We split the volumes into training (64%), validation (16%), and test (20%) sets. To conduct the simulated motion study, we synthesized rigid transformations across a wide range of parameters and applied them to the reference volumes. Slices were then sampled according to the acquisition protocol to generate 2,000, 500, and 200 3D volume-2D slice pairs for the training, validation, and test sets, respectively. Our experimental results demonstrate that our model achieves competitive performance in terms of alignment accuracy compared to state-of-the-art deep learning-based methods (Euclidean distance of 0.93 [mm] vs. 1.86 [mm], a paired t-test with a p-value of p<0.03). Furthermore, our approach exhibits faster registration speed compared to conventional iterative methods (0.096 s vs. 1.17 s). Our end-to-end SVR model facilitates real-time head motion tracking during fMRI acquisition, ensuring reliability and robustness against uncertainties in the inputs.

R

Coll-Font, Jaume, Onur Afacan, Scott Hoge, Harsha Garg, Kumar Shashi, Bahram Marami, Ali Gholipour, Jeanne Chow, Simon Warfield, and Sila Kurugol. 2021. “Retrospective Distortion and Motion Correction for Free-Breathing DW-MRI of the Kidneys Using Dual-Echo EPI and Slice-to-Volume Registration”. J Magn Reson Imaging 53 (5): 1432-43. https://doi.org/10.1002/jmri.27473.
BACKGROUND: Diffusion-weighted MRI (DW-MRI) of the kidneys is a technique that provides information about the microstructure of renal tissue without requiring exogenous contrasts such as gadolinium, and it can be used for diagnosis in cases of renal disease and assessing response-to-therapy. However, physiological motion and large geometric distortions due to main B0 field inhomogeneities degrade the image quality, reduce the accuracy of quantitative imaging markers, and impede their subsequent clinical applicability. PURPOSE: To retrospectively correct for geometric distortion for free-breathing DW-MRI of the kidneys at 3T, in the presence of a nonstatic distortion field due to breathing and bulk motion. STUDY TYPE: Prospective. SUBJECTS: Ten healthy volunteers (ages 29-38, four females). FIELD STRENGTH/SEQUENCE: 3T; DW-MR dual-echo echo-planar imaging (EPI) sequence (10 b-values and 17 directions) and a T2 volume. ASSESSMENT: The distortion correction was evaluated subjectively (Likert scale 0-5) and numerically with cross-correlation between the DW images at b = 0 s/mm2 and a T2 volume. The intravoxel incoherent motion (IVIM) and diffusion tensor (DTI) model-fitting performance was evaluated using the root-mean-squared error (nRMSE) and the coefficient of variation (CV%) of their parameters. STATISTICAL TESTS: Statistical comparisons were done using Wilcoxon tests. RESULTS: The proposed method improved the Likert scores by 1.1 ± 0.8 (P < 0.05), the cross-correlation with the T2 reference image by 0.13 ± 0.05 (P < 0.05), and reduced the nRMSE by 0.13 ± 0.03 (P < 0.05) and 0.23 ± 0.06 (P < 0.05) for IVIM and DTI, respectively. The CV% of the IVIM parameters (slow and fast diffusion, and diffusion fraction for IVIM and mean diffusivity, and fractional anisotropy for DTI) was reduced by 2.26 ± 3.98% (P = 6.971 × 10-2 ), 11.24 ± 26.26% (P = 6.971 × 10-2 ), 4.12 ± 12.91% (P = 0.101), 3.22 ± 0.55% (P < 0.05), and 2.42 ± 1.15% (P < 0.05). DATA CONCLUSION: The results indicate that the proposed Di + MoCo method can effectively correct for time-varying geometric distortions and for misalignments due to breathing motion. Consequently, the image quality and precision of the DW-MRI model parameters improved. LEVEL OF EVIDENCE: 2 TECHNICAL EFFICACY STAGE: 1.
Afacan, Onur, Tess Wallace, and Simon Warfield. 2020. “Retrospective Correction of Head Motion Using Measurements from an Electromagnetic Tracker”. Magn Reson Med 83 (2): 427-37. https://doi.org/10.1002/mrm.27934.
PURPOSE: To investigate the feasibility of using an electromagnetic (EM) tracker to estimate rigid body head motion parameters, and using these measurements to retrospectively reduce motion artifacts. THEORY AND METHODS: A clinically used MPRAGE sequence was modified to measure motion using the EM tracking system once per repetition time. A retrospective k-space based motion correction algorithm that corrects for phase ramps (translation in image domain) and rotation of 3D k-space (rotation in image domain) was developed, using the parameters recorded using an EM tracker. The accuracy of the EM tracker for the purpose of motion measurement and correction was tested in phantoms, volunteers, and pediatric patients. RESULTS: Position localization was accurate to the order of 200 microns compared with registration localization in a phantom study. The quality of reconstructed images was assessed by computing the root mean square error, the structural similarity metric and average edge strength. Image quality improved consistently when motion correction was applied in both volunteer scans with deliberate head motion and in pediatric patient scans. In patients, the average edge strength improved significantly with retrospective motion correction, compared with images with no correction applied. CONCLUSIONS: EM tracking was effective in measuring head motion in the MRI scanner with high accuracy, and enabled retrospective reconstruction to improve image quality by reducing motion artifacts.
Rollins, Caitlin, Cynthia Ortinau, Christian Stopp, Kevin Friedman, Wayne Tworetzky, Borjan Gagoski, Clemente Velasco-Annis, et al. 2021. “Regional Brain Growth Trajectories in Fetuses With Congenital Heart Disease”. Ann Neurol 89 (1): 143-57. https://doi.org/10.1002/ana.25940.
OBJECTIVE: Congenital heart disease (CHD) is associated with abnormal brain development in utero. We applied innovative fetal magnetic resonance imaging (MRI) techniques to determine whether reduced fetal cerebral substrate delivery impacts the brain globally, or in a region-specific pattern. Our novel design included two control groups, one with and the other without a family history of CHD, to explore the contribution of shared genes and/or fetal environment to brain development. METHODS: From 2014 to 2018, we enrolled 179 pregnant women into 4 groups: "HLHS/TGA" fetuses with hypoplastic left heart syndrome (HLHS) or transposition of the great arteries (TGA), diagnoses with lowest fetal cerebral substrate delivery; "CHD-other," with other CHD diagnoses; "CHD-related," healthy with a CHD family history; and "optimal control," healthy without a family history. Two MRIs were obtained between 18 and 40 weeks gestation. Random effect regression models assessed group differences in brain volumes and relationships to hemodynamic variables. RESULTS: HLHS/TGA (n = 24), CHD-other (50), and CHD-related (34) groups each had generally smaller brain volumes than the optimal controls (71). Compared with CHD-related, the HLHS/TGA group had smaller subplate (-13.3% [standard error = 4.3%], p < 0.01) and intermediate (-13.7% [4.3%], p < 0.01) zones, with a similar trend in ventricular zone (-7.1% [1.9%], p = 0.07). These volumetric reductions were associated with lower cerebral substrate delivery. INTERPRETATION: Fetuses with CHD, especially those with lowest cerebral substrate delivery, show a region-specific pattern of small brain volumes and impaired brain growth before 32 weeks gestation. The brains of fetuses with CHD were more similar to those of CHD-related than optimal controls, suggesting genetic or environmental factors also contribute. ANN NEUROL 2021;89:143-157.