Journal Papers

2015

Gholipour, Ali, Onur Afacan, Iman Aganj, Benoit Scherrer, Sanjay Prabhu, Mustafa Sahin, and Simon Warfield. (2015) 2015. “Super-Resolution Reconstruction in Frequency, Image, and Wavelet Domains to Reduce Through-Plane Partial Voluming in MRI”. Med Phys 42 (12): 6919-32. https://doi.org/10.1118/1.4935149.
PURPOSE: To compare and evaluate the use of super-resolution reconstruction (SRR), in frequency, image, and wavelet domains, to reduce through-plane partial voluming effects in magnetic resonance imaging. METHODS: The reconstruction of an isotropic high-resolution image from multiple thick-slice scans has been investigated through techniques in frequency, image, and wavelet domains. Experiments were carried out with thick-slice T2-weighted fast spin echo sequence on the Academic College of Radiology MRI phantom, where the reconstructed images were compared to a reference high-resolution scan using peak signal-to-noise ratio (PSNR), structural similarity image metric (SSIM), mutual information (MI), and the mean absolute error (MAE) of image intensity profiles. The application of super-resolution reconstruction was then examined in retrospective processing of clinical neuroimages of ten pediatric patients with tuberous sclerosis complex (TSC) to reduce through-plane partial voluming for improved 3D delineation and visualization of thin radial bands of white matter abnormalities. RESULTS: Quantitative evaluation results show improvements in all evaluation metrics through super-resolution reconstruction in the frequency, image, and wavelet domains, with the highest values obtained from SRR in the image domain. The metric values for image-domain SRR versus the original axial, coronal, and sagittal images were PSNR = 32.26 vs 32.22, 32.16, 30.65; SSIM = 0.931 vs 0.922, 0.924, 0.918; MI = 0.871 vs 0.842, 0.844, 0.831; and MAE = 5.38 vs 7.34, 7.06, 6.19. All similarity metrics showed high correlations with expert ranking of image resolution with MI showing the highest correlation at 0.943. Qualitative assessment of the neuroimages of ten TSC patients through in-plane and out-of-plane visualization of structures showed the extent of partial voluming effect in a real clinical scenario and its reduction using SRR. Blinded expert evaluation of image resolution in resampled out-of-plane views consistently showed the superiority of SRR compared to original axial and coronal image acquisitions. CONCLUSIONS: Thick-slice 2D T2-weighted MRI scans are part of many routine clinical protocols due to their high signal-to-noise ratio, but are often severely affected by through-plane partial voluming effects. This study shows that while radiologic assessment is performed in 2D on thick-slice scans, super-resolution MRI reconstruction techniques can be used to fuse those scans to generate a high-resolution image with reduced partial voluming for improved postacquisition processing. Qualitative and quantitative evaluation showed the efficacy of all SRR techniques with the best results obtained from SRR in the image domain. The limitations of SRR techniques are uncertainties in modeling the slice profile, density compensation, quantization in resampling, and uncompensated motion between scans.
Kurugol, Sila, Moti Freiman, Onur Afacan, Liran Domachevsky, Jeannette Perez-Rossello, Michael Callahan, and Simon Warfield. (2015) 2015. “Motion Compensated Abdominal Diffusion Weighted MRI by Simultaneous Image Registration and Model Estimation (SIR-ME)”. Med Image Comput Comput Assist Interv 9351: 501-9. https://doi.org/10.1007/978-3-319-24574-4_60.
Non-invasive characterization of water molecule's mobility variations by quantitative analysis of diffusion-weighted MRI (DW-MRI) signal decay in the abdomen has the potential to serve as a biomarker in gastrointestinal and oncological applications. Accurate and reproducible estimation of the signal decay model parameters is challenging due to the presence of respiratory, cardiac, and peristalsis motion. Independent registration of each b-value image to the b-value=0 s/mm(2) image prior to parameter estimation might be sub-optimal because of the low SNR and contrast difference between images of varying b-value. In this work, we introduce a motion-compensated parameter estimation framework that simultaneously solves image registration and model estimation (SIR-ME) problems by utilizing the interdependence of acquired volumes along the diffusion weighting dimension. We evaluated the improvement in model parameters estimation accuracy using 16 in-vivo DW-MRI data sets of Crohn's disease patients by comparing parameter estimates obtained using the SIR-ME model to the parameter estimates obtained by fitting the signal decay model to the acquired DW-MRI images. The proposed SIR-ME model reduced the average root-mean-square error between the observed signal and the fitted model by more than 50%. Moreover, the SIR-ME model estimates discriminate between normal and abnormal bowel loops better than the standard parameter estimates.
Stamm, Aymeric, Jolene Singh, Onur Afacan, and Simon Warfield. (2015) 2015. “Analytic Quantification of Bias and Variance of Coil Sensitivity Profile Estimators for Improved Image Reconstruction in MRI”. Med Image Comput Comput Assist Interv 9350: 684-91. https://doi.org/10.1007/978-3-319-24571-3_82.
Magnetic resonance (MR) imaging provides a unique in-vivo capability of visualizing tissue in the human brain non-invasively, which has tremendously improved patient care over the past decades. However, there are still prominent artifacts, such as intensity inhomogeneities due to the use of an array of receiving coils (RC) to measure the MR signal or noise amplification due to accelerated imaging strategies. It is critical to mitigate these artifacts for both visual inspection and quantitative analysis. The cornerstone to address this issue pertains to the knowledge of coil sensitivity profiles (CSP) of the RCs, which describe how the measured complex signal decays with the distance to the RC. Existing methods for CSP estimation share a number of limitations: (i) they primarily focus on CSP magnitude, while it is known that the solution to the MR image reconstruction problem involves complex CSPs and (ii) they only provide point estimates of the CSPs, which makes the task of optimizing the parameters and acquisition protocol for their estimation difficult. In this paper, we propose a novel statistical framework for estimating complex-valued CSPs. We define a CSP estimator that uses spatial smoothing and additional body coil data for phase normalization. The main contribution is to provide detailed information on the statistical distribution of the CSP estimator, which yields automatic determination of the optimal degree of smoothing for ensuring minimal bias and provides guidelines to the optimal acquisition strategy.
Scherrer, Benoit, Onur Afacan, Maxime Taquet, Sanjay Prabhu, Ali Gholipour, and Simon Warfield. (2015) 2015. “Accelerated High Spatial Resolution Diffusion-Weighted Imaging”. Inf Process Med Imaging 24: 69-81. https://doi.org/10.1007/978-3-319-19992-4_6.
Acquisition of a series of anisotropically oversampled acquisitions (so-called anisotropic "snapshots") and reconstruction in the image space has recently been proposed to increase the spatial resolution in diffusion weighted imaging (DWI), providing a theoretical 8x acceleration at equal signal-to-noise ratio (SNR) compared to conventional dense k-space sampling. However, in most works, each DW image is reconstructed separately and the fact that the DW images constitute different views of the same anatomy is ignored. In addition, current approaches are limited by their inability to reconstruct a high resolution (HR) acquisition from snapshots with different subsets of diffusion gradients: an isotropic HR gradient image cannot be reconstructed if one .of its anisotropic snapshots is missing, for example due to intra-scan motion, even if other snapshots for this gradient were successfully acquired. In this work, we propose a novel multi-snapshot DWI reconstruction technique that simultaneously achieves HR reconstruction and local tissue model estimation while enabling reconstruction from snapshots containing different subsets of diffusion gradients, providing increased robustness to patient motion and potential for acceleration. Our approach is formalized as a joint probabilistic model with missing observations, from which interactions between missing snapshots, HR reconstruction and a generic tissue model naturally emerge. We evaluate our approach with synthetic simulations, simulated multi-snapshot scenario and in vivo multi-snapshot imaging. We show that (1) our combined approach ultimately provides both better HR reconstruction and better tissue model estimation and (2) the error in the case of missing snapshots can be quantified. Our novel multi-snapshot technique will enable improved high spatial characterization of the brain connectivity and microstructure in vivo.
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.
Velasco-Annis, Clemente, Ali Gholipour, Onur Afacan, Sanjay Prabhu, Judy Estroff, and Simon Warfield. (2015) 2015. “Normative Biometrics for Fetal Ocular Growth Using Volumetric MRI Reconstruction”. Prenat Diagn 35 (4): 400-8. https://doi.org/10.1002/pd.4558.
OBJECTIVE: To determine normative ranges for fetal ocular biometrics between 19 and 38 weeks gestational age (GA) using volumetric MRI reconstruction. METHOD: The 3D images of 114 healthy fetuses between 19 and 38 weeks GA were created using super-resolution volume reconstructions from MRI slice acquisitions. These 3D images were semi-automatically segmented to measure fetal orbit volume, binocular distance (BOD), interocular distance (IOD), and ocular diameter (OD). RESULTS: All biometry correlated with GA (Volume, Pearson's correlation coefficient (CC) = 0.9680; BOD, CC = 0.9552; OD, CC = 0.9445; and IOD, CC = 0.8429), and growth curves were plotted against linear and quadratic growth models. Regression analysis showed quadratic models to best fit BOD, IOD, and OD and a linear model to best fit volume. CONCLUSION: Orbital volume had the greatest correlation with GA, although BOD and OD also showed strong correlation. The normative data found in this study may be helpful for the detection of congenital fetal anomalies with more consistent measurements than are currently available. © 2015 John Wiley & Sons, Ltd.
Kurugol, Sila, Moti Freiman, Onur Afacan, Liran Domachevsky, Jeannette M Perez-Rossello, Michael J Callahan, and Simon K Warfield. (2015) 2015. “Motion Compensated Abdominal Diffusion Weighted MRI by Simultaneous Image Registration and Model Estimation (SIR-ME).”. Medical Image Computing and Computer-Assisted Intervention : MICCAI . International Conference on Medical Image Computing and Computer-Assisted Intervention 9351: 501-9. https://doi.org/10.1007/978-3-319-24574-4_60.

Non-invasive characterization of water molecule's mobility variations by quantitative analysis of diffusion-weighted MRI (DW-MRI) signal decay in the abdomen has the potential to serve as a biomarker in gastrointestinal and oncological applications. Accurate and reproducible estimation of the signal decay model parameters is challenging due to the presence of respiratory, cardiac, and peristalsis motion. Independent registration of each b-value image to the b-value=0 s/mm(2) image prior to parameter estimation might be sub-optimal because of the low SNR and contrast difference between images of varying b-value. In this work, we introduce a motion-compensated parameter estimation framework that simultaneously solves image registration and model estimation (SIR-ME) problems by utilizing the interdependence of acquired volumes along the diffusion weighting dimension. We evaluated the improvement in model parameters estimation accuracy using 16 in-vivo DW-MRI data sets of Crohn's disease patients by comparing parameter estimates obtained using the SIR-ME model to the parameter estimates obtained by fitting the signal decay model to the acquired DW-MRI images. The proposed SIR-ME model reduced the average root-mean-square error between the observed signal and the fitted model by more than 50%. Moreover, the SIR-ME model estimates discriminate between normal and abnormal bowel loops better than the standard parameter estimates.

2014

Gholipour, Ali, Judith Estroff, Carol Barnewolt, Richard Robertson, Ellen Grant, Borjan Gagoski, Simon Warfield, et al. (2014) 2014. “Fetal MRI: A Technical Update With Educational Aspirations”. Concepts Magn Reson Part A Bridg Educ Res 43 (6): 237-66. https://doi.org/10.1002/cmr.a.21321.
Fetal magnetic resonance imaging (MRI) examinations have become well-established procedures at many institutions and can serve as useful adjuncts to ultrasound (US) exams when diagnostic doubts remain after US. Due to fetal motion, however, fetal MRI exams are challenging and require the MR scanner to be used in a somewhat different mode than that employed for more routine clinical studies. Herein we review the techniques most commonly used, and those that are available, for fetal MRI with an emphasis on the physics of the techniques and how to deploy them to improve success rates for fetal MRI exams. By far the most common technique employed is single-shot T2-weighted imaging due to its excellent tissue contrast and relative immunity to fetal motion. Despite the significant challenges involved, however, many of the other techniques commonly employed in conventional neuro- and body MRI such as T1 and T2*-weighted imaging, diffusion and perfusion weighted imaging, as well as spectroscopic methods remain of interest for fetal MR applications. An effort to understand the strengths and limitations of these basic methods within the context of fetal MRI is made in order to optimize their use and facilitate implementation of technical improvements for the further development of fetal MR imaging, both in acquisition and post-processing strategies.
Akhondi-Asl, Alireza, Onur Afacan, Robert Mulkern, and Simon Warfield. (2014) 2014. “T(2)-Relaxometry for Myelin Water Fraction Extraction Using Wald Distribution and Extended Phase Graph”. Med Image Comput Comput Assist Interv 17 (Pt 3): 145-52. https://doi.org/10.1007/978-3-319-10443-0_19.
Quantitative assessment of myelin density in the white matter is an emerging tool for neurodegenerative disease related studies such as multiple sclerosis and Schizophrenia. For the last two decades, T2 relaxometry based on multi-exponential fitting to a single slice multi-echo sequence has been the most common MRI technique for myelin water fraction (MWF) mapping, where the short T2 is associated with myelin water. However, modeling the spectrum of the relaxations as the sum of large number of impulse functions with unknown amplitudes makes the accuracy and robustness of the estimated MWF's questionable. In this paper, we introduce a novel model with small number of parameters to simultaneously characterize transverse relaxation rate spectrum and B1 inhomogeneity at each voxel. We use mixture of three Wald distributions with unknown mixture weights, mean and shape parameters to represent the distribution of the relative amount of water in between myelin sheets, tissue water, and cerebrospinal fluid. The parameters of the model are estimated using the variable projection method and are used to extract the MWF at each voxel. In addition, we use Extended Phase Graph (EPG) method to compensate for the stimulated echoes caused by B1 inhomogeneity. To validate our model, synthetic and real brain experiments were conducted where we have compared our novel algorithm with the non-negative least squares (NNLS) as the state-of-the-art technique in the literature. Our results indicate that we can estimate MWF map with substantially higher accuracy as compared to the NNLS method.

2013

Freiman, Moti, Onur Afacan, Robert Mulkern, and Simon Warfield. (2013) 2013. “Improved Multi B-Value Diffusion-Weighted MRI of the Body by Simultaneous Model Estimation and Image Reconstruction (SMEIR)”. Med Image Comput Comput Assist Interv 16 (Pt 3): 1-8. https://doi.org/10.1007/978-3-642-40760-4_1.
Diffusion-weighted MRI images acquired with multiple b-values have the potential to improve diagnostic accuracy by increasing the conspicuity of lesions and inflammatory activity with background suppression. Unfortunately, the inherently low signal-to-noise ratio (SNR) of DW-MRI reduces enthusiasm for using these images for diagnostic purposes. Moreover, lengthy acquisition times limit our ability to improve the quality of multi b-value DW-MRI images by multiple excitations acquisition and signal averaging at each b-value. To offset these limitations, we propose the Simultaneous Model Estimation and Image Reconstruction (SMEIR) for DW-MRI, which substantially improves the quality of multi b-value DW-MRI images without increasing acquisition times. Our model introduces the physiological signal decay model of DW-MRI as a constraint in the reconstruction of the DW-MRI images. An in-vivo experiment using 6 low-quality DW-MRI datasets of a healthy subject showed that SMEIR reconstruction of low-quality data improved SNR by 55% in the liver and by 41% in the kidney without increasing acquisition times. We also demonstrated the clinical impact of our SMEIR reconstruction by increasing the conspicuity of inflamed bowel regions in DW-MRI of 12 patients with Crohn's disease. The contrast-to-noise ratio (CNR) of the inflamed regions in the SMEIR images was higher by 12.6% relative to CNR in the original DW-MRI images.