1
|
Tang X, Gao J, Aburas A, Wu D, Chen Z, Chen H, Hu C. Accelerated multi-b-value multi-shot diffusion-weighted imaging based on EPI with keyhole and a low-rank tensor constraint. Magn Reson Imaging 2024; 110:138-148. [PMID: 38641211 DOI: 10.1016/j.mri.2024.04.015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2024] [Revised: 04/08/2024] [Accepted: 04/10/2024] [Indexed: 04/21/2024]
Abstract
PURPOSE Multi-Shot (MS) Echo-Planar Imaging (EPI) may improve the in-plane resolution of multi-b-value DWI, yet it also considerably increases the scan time. Here we explored the combination of EPI with Keyhole (EPIK) and a calibrationless reconstruction algorithm for acceleration of multi-b-value MS-DWI. METHODS We firstly analyzed the impact of nonuniform phase accrual in EPIK on the reconstructed image. Based on insights gained from the analysis, we developed a calibrationless reconstruction algorithm based on a Space-Contrast-Coil Locally Low-Rank Tensor (SCC-LLRT) constraint for reconstruction of EPIK-acquired data. We compared the algorithm with a modified SPatial-Angular Locally Low-Rank (SPA-LLR) algorithm through simulations, phantoms, and in vivo study. We then compared EPIK with uniformly undersampled EPI for accelerating multi-b-value DWI in 6 healthy subjects. RESULTS Through theoretical derivations, we found that the reconstruction of EPIK with a SENSE-encoding-based algorithm, such as SPA-LLR, may cause additional aliasing artifacts due to the frequency-dependent distortion of the coil sensitivity. Results from simulations, phantoms, and in vivo study verified the theoretical finding by showing that the calibrationless SCC-LLRT algorithm reduced aliasing artifacts compared with SPA-LLR. Finally, EPIK with SCC-LLRT substantially reduced the ghosting artifacts compared with uniform undersampled multi-b-value DWI, decreasing the fitting errors in ADC (0.05 ± 0.01 vs 0.10 ± 0.01, P < 0.001) and IVIM mapping (0.026 ± 0.004 vs 0.06 ± 0.006, P < 0.001). CONCLUSION The SCC-LLRT algorithm reduced the aliasing artifacts of EPIK by using a calibrationless modeling of the multi-coil data. The dense sampling of k-space center offers EPIK a potential to improve image quality for acceleration of multi-b-value MS-DWI.
Collapse
Affiliation(s)
- Xin Tang
- National Engineering Research Center of Advanced Magnetic Resonance Technologies for Diagnosis and Therapy, School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai, China; United Imaging Healthcare Co. Ltd, Shanghai, China
| | - Juan Gao
- National Engineering Research Center of Advanced Magnetic Resonance Technologies for Diagnosis and Therapy, School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai, China
| | - Ahmed Aburas
- National Engineering Research Center of Advanced Magnetic Resonance Technologies for Diagnosis and Therapy, School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai, China
| | - Dan Wu
- Department of Biomedical Engineering, College of Biomedical Engineering & Instrument Science, Zhejiang University, Hangzhou, Zhejiang, China
| | - Zhuo Chen
- National Engineering Research Center of Advanced Magnetic Resonance Technologies for Diagnosis and Therapy, School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai, China
| | - Hao Chen
- National Engineering Research Center of Advanced Magnetic Resonance Technologies for Diagnosis and Therapy, School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai, China
| | - Chenxi Hu
- National Engineering Research Center of Advanced Magnetic Resonance Technologies for Diagnosis and Therapy, School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai, China.
| |
Collapse
|
2
|
Radhakrishnan H, Zhao C, Sydnor VJ, Baller EB, Cook PA, Fair DA, Giesbrecht B, Larsen B, Murtha K, Roalf DR, Rush‐Goebel S, Shinohara RT, Shou H, Tisdall MD, Vettel JM, Grafton ST, Cieslak M, Satterthwaite TD. A practical evaluation of measures derived from compressed sensing diffusion spectrum imaging. Hum Brain Mapp 2024; 45:e26580. [PMID: 38520359 PMCID: PMC10960521 DOI: 10.1002/hbm.26580] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2023] [Revised: 12/13/2023] [Accepted: 12/20/2023] [Indexed: 03/25/2024] Open
Abstract
Diffusion Spectrum Imaging (DSI) using dense Cartesian sampling of q-space has been shown to provide important advantages for modeling complex white matter architecture. However, its adoption has been limited by the lengthy acquisition time required. Sparser sampling of q-space combined with compressed sensing (CS) reconstruction techniques has been proposed as a way to reduce the scan time of DSI acquisitions. However prior studies have mainly evaluated CS-DSI in post-mortem or non-human data. At present, the capacity for CS-DSI to provide accurate and reliable measures of white matter anatomy and microstructure in the living human brain remains unclear. We evaluated the accuracy and inter-scan reliability of 6 different CS-DSI schemes that provided up to 80% reductions in scan time compared to a full DSI scheme. We capitalized on a dataset of 26 participants who were scanned over eight independent sessions using a full DSI scheme. From this full DSI scheme, we subsampled images to create a range of CS-DSI images. This allowed us to compare the accuracy and inter-scan reliability of derived measures of white matter structure (bundle segmentation, voxel-wise scalar maps) produced by the CS-DSI and the full DSI schemes. We found that CS-DSI estimates of both bundle segmentations and voxel-wise scalars were nearly as accurate and reliable as those generated by the full DSI scheme. Moreover, we found that the accuracy and reliability of CS-DSI was higher in white matter bundles that were more reliably segmented by the full DSI scheme. As a final step, we replicated the accuracy of CS-DSI in a prospectively acquired dataset (n = 20, scanned once). Together, these results illustrate the utility of CS-DSI for reliably delineating in vivo white matter architecture in a fraction of the scan time, underscoring its promise for both clinical and research applications.
Collapse
Affiliation(s)
- Hamsanandini Radhakrishnan
- Lifespan Informatics and Neuroimaging CenterUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
- Department of Psychiatry, Perelman School of MedicineUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
| | - Chenying Zhao
- Lifespan Informatics and Neuroimaging CenterUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
- Department of Psychiatry, Perelman School of MedicineUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
- Lifespan Brain Institute, Perelman School of MedicineUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
- Department of Bioengineering, School of Engineering and Applied ScienceUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
| | - Valerie J. Sydnor
- Lifespan Informatics and Neuroimaging CenterUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
- Department of Psychiatry, Perelman School of MedicineUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
| | - Erica B. Baller
- Lifespan Informatics and Neuroimaging CenterUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
- Department of Psychiatry, Perelman School of MedicineUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
| | - Philip A. Cook
- Department of Radiology, Perelman School of MedicineUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
| | - Damien A. Fair
- Masonic Institute for the Developing BrainUniversity of MinnesotaMinneapolisMinnesotaUSA
| | - Barry Giesbrecht
- Department of Psychological and Brain SciencesUniversity of CaliforniaSanta BarbaraCaliforniaUSA
| | - Bart Larsen
- Lifespan Informatics and Neuroimaging CenterUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
- Department of Psychiatry, Perelman School of MedicineUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
| | - Kristin Murtha
- Lifespan Informatics and Neuroimaging CenterUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
- Department of Psychiatry, Perelman School of MedicineUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
| | - David R. Roalf
- Department of Psychiatry, Perelman School of MedicineUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
- Lifespan Brain Institute, Perelman School of MedicineUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
| | - Sage Rush‐Goebel
- Lifespan Informatics and Neuroimaging CenterUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
- Department of Psychiatry, Perelman School of MedicineUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
| | - Russell T. Shinohara
- Department of Biostatistics, Epidemiology and InformaticsUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
- Center for Biomedical Image Computing & AnalyticsUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
| | - Haochang Shou
- Department of Biostatistics, Epidemiology and InformaticsUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
- Center for Biomedical Image Computing & AnalyticsUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
| | - M. Dylan Tisdall
- Department of Radiology, Perelman School of MedicineUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
| | - Jean M. Vettel
- Department of Psychological and Brain SciencesUniversity of CaliforniaSanta BarbaraCaliforniaUSA
- U.S. Army Research LaboratoryAberdeen Proving GroundAberdeenMarylandUSA
| | - Scott T. Grafton
- Department of Psychological and Brain SciencesUniversity of CaliforniaSanta BarbaraCaliforniaUSA
| | - Matthew Cieslak
- Lifespan Informatics and Neuroimaging CenterUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
- Department of Psychiatry, Perelman School of MedicineUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
| | - Theodore D. Satterthwaite
- Lifespan Informatics and Neuroimaging CenterUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
- Department of Psychiatry, Perelman School of MedicineUniversity of PennsylvaniaPhiladelphiaPennsylvaniaUSA
| |
Collapse
|
3
|
Radhakrishnan H, Zhao C, Sydnor VJ, Baller EB, Cook PA, Fair D, Giesbrecht B, Larsen B, Murtha K, Roalf DR, Rush-Goebel S, Shinohara R, Shou H, Tisdall MD, Vettel J, Grafton S, Cieslak M, Satterthwaite T. Establishing the Validity of Compressed Sensing Diffusion Spectrum Imaging. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.02.22.529546. [PMID: 36865219 PMCID: PMC9980087 DOI: 10.1101/2023.02.22.529546] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/25/2023]
Abstract
Diffusion Spectrum Imaging (DSI) using dense Cartesian sampling of q-space has been shown to provide important advantages for modeling complex white matter architecture. However, its adoption has been limited by the lengthy acquisition time required. Sparser sampling of q-space combined with compressed sensing (CS) reconstruction techniques has been proposed as a way to reduce the scan time of DSI acquisitions. However prior studies have mainly evaluated CS-DSI in post-mortem or non-human data. At present, the capacity for CS-DSI to provide accurate and reliable measures of white matter anatomy and microstructure in the living human brain remains unclear. We evaluated the accuracy and inter-scan reliability of 6 different CS-DSI schemes that provided up to 80% reductions in scan time compared to a full DSI scheme. We capitalized on a dataset of twenty-six participants who were scanned over eight independent sessions using a full DSI scheme. From this full DSI scheme, we subsampled images to create a range of CS-DSI images. This allowed us to compare the accuracy and inter-scan reliability of derived measures of white matter structure (bundle segmentation, voxel-wise scalar maps) produced by the CS-DSI and the full DSI schemes. We found that CS-DSI estimates of both bundle segmentations and voxel-wise scalars were nearly as accurate and reliable as those generated by the full DSI scheme. Moreover, we found that the accuracy and reliability of CS-DSI was higher in white matter bundles that were more reliably segmented by the full DSI scheme. As a final step, we replicated the accuracy of CS-DSI in a prospectively acquired dataset (n=20, scanned once). Together, these results illustrate the utility of CS-DSI for reliably delineating in vivo white matter architecture in a fraction of the scan time, underscoring its promise for both clinical and research applications.
Collapse
Affiliation(s)
- Hamsanandini Radhakrishnan
- Lifespan Informatics and Neuroimaging Center, University of Pennsylvania, Philadelphia, PA, USA
- Department of Psychiatry, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Chenying Zhao
- Lifespan Informatics and Neuroimaging Center, University of Pennsylvania, Philadelphia, PA, USA
- Department of Psychiatry, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
- Lifespan Brain Institute, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
- Department of Bioengineering, School of Engineering and Applied Science, University of Pennsylvania, Philadelphia, PA, USA
| | - Valerie J. Sydnor
- Lifespan Informatics and Neuroimaging Center, University of Pennsylvania, Philadelphia, PA, USA
- Department of Psychiatry, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Erica B. Baller
- Lifespan Informatics and Neuroimaging Center, University of Pennsylvania, Philadelphia, PA, USA
- Department of Psychiatry, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Philip A. Cook
- Department of Radiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Damien Fair
- Masonic Institute for the Developing Brain, University of Minnesota, Minneapolis, MN, USA
| | - Barry Giesbrecht
- Department of Psychological and Brain Sciences, University of California, Santa Barbara, CA, USA
| | - Bart Larsen
- Lifespan Informatics and Neuroimaging Center, University of Pennsylvania, Philadelphia, PA, USA
- Department of Psychiatry, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Kristin Murtha
- Lifespan Informatics and Neuroimaging Center, University of Pennsylvania, Philadelphia, PA, USA
- Department of Psychiatry, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - David R. Roalf
- Department of Psychiatry, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
- Lifespan Brain Institute, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Sage Rush-Goebel
- Lifespan Informatics and Neuroimaging Center, University of Pennsylvania, Philadelphia, PA, USA
- Department of Psychiatry, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Russell Shinohara
- Department of Biostatistics, Epidemiology and Informatics, University of Pennsylvania, Philadelphia, PA, USA
- Center for Biomedical Image Computing & Analytics, University of Pennsylvania, Philadelphia, PA, USA
| | - Haochang Shou
- Department of Biostatistics, Epidemiology and Informatics, University of Pennsylvania, Philadelphia, PA, USA
- Center for Biomedical Image Computing & Analytics, University of Pennsylvania, Philadelphia, PA, USA
| | - M. Dylan Tisdall
- Department of Radiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Jean Vettel
- Department of Psychological and Brain Sciences, University of California, Santa Barbara, CA, USA
- U.S. Army Research Laboratory, Aberdeen Proving Ground, MD, USA
| | - Scott Grafton
- Department of Psychological and Brain Sciences, University of California, Santa Barbara, CA, USA
| | - Matthew Cieslak
- Lifespan Informatics and Neuroimaging Center, University of Pennsylvania, Philadelphia, PA, USA
- Department of Psychiatry, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Theodore Satterthwaite
- Lifespan Informatics and Neuroimaging Center, University of Pennsylvania, Philadelphia, PA, USA
- Department of Psychiatry, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| |
Collapse
|
4
|
Wang Z, She H, Zhang Y, Du YP. Parallel non-Cartesian spatial-temporal dictionary learning neural networks (stDLNN) for accelerating 4D-MRI. Med Image Anal 2023; 84:102701. [PMID: 36470148 DOI: 10.1016/j.media.2022.102701] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2021] [Revised: 10/02/2022] [Accepted: 11/21/2022] [Indexed: 11/25/2022]
Abstract
Dynamic magnetic resonance imaging (MRI) acquisitions are relatively slow due to physical and physiological limitations. The spatial-temporal dictionary learning (DL) approach accelerates dynamic MRI by learning spatial-temporal correlations, but the regularization parameters need to be manually adjusted, the performance at high acceleration rate is limited, and the reconstruction can be time-consuming. Deep learning techniques have shown good performance in accelerating MRI due to the powerful representational capabilities of neural networks. In this work, we propose a parallel non-Cartesian spatial-temporal dictionary learning neural networks (stDLNN) framework that combines dictionary learning with deep learning algorithms and utilizes the spatial-temporal prior information of dynamic MRI data to achieve better reconstruction quality and efficiency. The coefficient estimation modules (CEM) are designed in the framework to adaptively adjust the regularization coefficients. Experimental results show that combining dictionary learning with deep neural networks and using spatial-temporal dictionaries can obviously improve the image quality and computational efficiency compared with the state-of-the-art non-Cartesian imaging methods for accelerating the 4D-MRI especially at high acceleration rate.
Collapse
Affiliation(s)
- Zhijun Wang
- School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai 200030, China
| | - Huajun She
- School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai 200030, China.
| | - Yufei Zhang
- School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai 200030, China
| | - Yiping P Du
- School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai 200030, China
| |
Collapse
|
5
|
Jones R, Maffei C, Augustinack J, Fischl B, Wang H, Bilgic B, Yendiki A. High-fidelity approximation of grid- and shell-based sampling schemes from undersampled DSI using compressed sensing: Post mortem validation. Neuroimage 2021; 244:118621. [PMID: 34587516 PMCID: PMC8631240 DOI: 10.1016/j.neuroimage.2021.118621] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2021] [Revised: 09/02/2021] [Accepted: 09/24/2021] [Indexed: 12/31/2022] Open
Abstract
While many useful microstructural indices, as well as orientation distribution functions, can be obtained from multi-shell dMRI data, there is growing interest in exploring the richer set of microstructural features that can be extracted from the full ensemble average propagator (EAP). The EAP can be readily computed from diffusion spectrum imaging (DSI) data, at the cost of a very lengthy acquisition. Compressed sensing (CS) has been used to make DSI more practical by reducing its acquisition time. CS applied to DSI (CS-DSI) attempts to reconstruct the EAP from significantly undersampled q-space data. We present a post mortem validation study where we evaluate the ability of CS-DSI to approximate not only fully sampled DSI but also multi-shell acquisitions with high fidelity. Human brain samples are imaged with high-resolution DSI at 9.4T and with polarization-sensitive optical coherence tomography (PSOCT). The latter provides direct measurements of axonal orientations at microscopic resolutions, allowing us to evaluate the mesoscopic orientation estimates obtained from diffusion MRI, in terms of their angular error and the presence of spurious peaks. We test two fast, dictionary-based, L2-regularized algorithms for CS-DSI reconstruction. We find that, for a CS acceleration factor of R=3, i.e., an acquisition with 171 gradient directions, one of these methods is able to achieve both low angular error and low number of spurious peaks. With a scan length similar to that of high angular resolution multi-shell acquisition schemes, this CS-DSI approach is able to approximate both fully sampled DSI and multi-shell data with high accuracy. Thus it is suitable for orientation reconstruction and microstructural modeling techniques that require either grid- or shell-based acquisitions. We find that the signal-to-noise ratio (SNR) of the training data used to construct the dictionary can have an impact on the accuracy of CS-DSI, but that there is substantial robustness to loss of SNR in the test data. Finally, we show that, as the CS acceleration factor increases beyond R=3, the accuracy of these reconstruction methods degrade, either in terms of the angular error, or in terms of the number of spurious peaks. Our results provide useful benchmarks for the future development of even more efficient q-space acceleration techniques.
Collapse
Affiliation(s)
- Robert Jones
- Department of Radiology, Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital & Harvard Medical School, Charlestown, MA 02129, USA.
| | - Chiara Maffei
- Department of Radiology, Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital & Harvard Medical School, Charlestown, MA 02129, USA
| | - Jean Augustinack
- Department of Radiology, Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital & Harvard Medical School, Charlestown, MA 02129, USA
| | - Bruce Fischl
- Department of Radiology, Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital & Harvard Medical School, Charlestown, MA 02129, USA; Computer Science and Artificial Intelligence Laboratory, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Hui Wang
- Department of Radiology, Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital & Harvard Medical School, Charlestown, MA 02129, USA
| | - Berkin Bilgic
- Department of Radiology, Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital & Harvard Medical School, Charlestown, MA 02129, USA; Harvard-MIT Division of Health Sciences and Technology, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Anastasia Yendiki
- Department of Radiology, Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital & Harvard Medical School, Charlestown, MA 02129, USA
| |
Collapse
|
6
|
Grisot G, Haber SN, Yendiki A. Diffusion MRI and anatomic tracing in the same brain reveal common failure modes of tractography. Neuroimage 2021; 239:118300. [PMID: 34171498 PMCID: PMC8475636 DOI: 10.1016/j.neuroimage.2021.118300] [Citation(s) in RCA: 35] [Impact Index Per Article: 11.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2021] [Revised: 04/29/2021] [Accepted: 06/21/2021] [Indexed: 12/15/2022] Open
Abstract
Anatomic tracing is recognized as a critical source of knowledge on brain circuitry that can be used to assess the accuracy of diffusion MRI (dMRI) tractography. However, most prior studies that have performed such assessments have used dMRI and tracer data from different brains and/or have been limited in the scope of dMRI analysis methods allowed by the data. In this work, we perform a quantitative, voxel-wise comparison of dMRI tractography and anatomic tracing data in the same macaque brain. An ex vivo dMRI acquisition with high angular resolution and high maximum b-value allows us to compare a range of q-space sampling, orientation reconstruction, and tractography strategies. The availability of tracing in the same brain allows us to localize the sources of tractography errors and to identify axonal configurations that lead to such errors consistently, across dMRI acquisition and analysis strategies. We find that these common failure modes involve geometries such as branching or turning, which cannot be modeled well by crossing fibers. We also find that the default thresholds that are commonly used in tractography correspond to rather conservative, low-sensitivity operating points. While deterministic tractography tends to have higher sensitivity than probabilistic tractography in that very conservative threshold regime, the latter outperforms the former as the threshold is relaxed to avoid missing true anatomical connections. On the other hand, the q-space sampling scheme and maximum b-value have less of an impact on accuracy. Finally, using scans from a set of additional macaque brains, we show that there is enough inter-individual variability to warrant caution when dMRI and tracer data come from different animals, as is often the case in the tractography validation literature. Taken together, our results provide insights on the limitations of current tractography methods and on the critical role that anatomic tracing can play in identifying potential avenues for improvement.
Collapse
Affiliation(s)
| | - Suzanne N Haber
- Department of Pharmacology and Physiology, University of Rochester, Rochester, NY, United States; McLean Hospital, Belmont, MA, United States
| | - Anastasia Yendiki
- Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital and Harvard Medical School, Charlestown, MA, United States.
| |
Collapse
|
7
|
Mani M, Magnotta VA, Jacob M. qModeL: A plug-and-play model-based reconstruction for highly accelerated multi-shot diffusion MRI using learned priors. Magn Reson Med 2021; 86:835-851. [PMID: 33759240 PMCID: PMC8076086 DOI: 10.1002/mrm.28756] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2020] [Revised: 02/08/2021] [Accepted: 02/09/2021] [Indexed: 01/19/2023]
Abstract
PURPOSE To introduce a joint reconstruction method for highly undersampled multi-shot diffusion weighted (msDW) scans. METHODS Multi-shot EPI methods enable higher spatial resolution for diffusion MRI, but at the expense of long scan-time. Highly accelerated msDW scans are needed to enable their utilization in advanced microstructure studies, which require high q-space coverage. Previously, joint k-q undersampling methods coupled with compressed sensing were shown to enable very high acceleration factors. However, the reconstruction of this data using sparsity priors is challenging and is not suited for multi-shell data. We propose a new reconstruction that recovers images from the combined k-q data jointly. The proposed qModeL reconstruction brings together the advantages of model-based iterative reconstruction and machine learning, extending the idea of plug-and-play algorithms. Specifically, qModeL works by prelearning the signal manifold corresponding to the diffusion measurement space using deep learning. The prelearned manifold prior is incorporated into a model-based reconstruction to provide a voxel-wise regularization along the q-dimension during the joint recovery. Notably, the learning does not require in vivo training data and is derived exclusively from biophysical modeling. Additionally, a plug-and-play total variation denoising provides regularization along the spatial dimension. The proposed framework is tested on k-q undersampled single-shell and multi-shell msDW acquisition at various acceleration factors. RESULTS The qModeL joint reconstruction is shown to recover DWIs from 8-fold accelerated msDW acquisitions with error less than 5% for both single-shell and multi-shell data. Advanced microstructural analysis performed using the undersampled reconstruction also report reasonable accuracy. CONCLUSION qModeL enables the joint recovery of highly accelerated multi-shot dMRI utilizing learning-based priors. The bio-physically driven approach enables the use of accelerated multi-shot imaging for multi-shell sampling and advanced microstructure studies.
Collapse
Affiliation(s)
- Merry Mani
- Department of Radiology, University of Iowa, Iowa City, Iowa
| | | | - Mathews Jacob
- Department of Electrical and Computer Engineering, University of Iowa, Iowa City, Iowa
| |
Collapse
|
8
|
Zhang C, Arefin TM, Nakarmi U, Lee CH, Li H, Liang D, Zhang J, Ying L. Acceleration of three-dimensional diffusion magnetic resonance imaging using a kernel low-rank compressed sensing method. Neuroimage 2020; 210:116584. [PMID: 32004717 DOI: 10.1016/j.neuroimage.2020.116584] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2019] [Revised: 12/23/2019] [Accepted: 01/23/2020] [Indexed: 12/13/2022] Open
Abstract
Diffusion Magnetic Resonance Imaging (dMRI) has shown great potential in probing tissue microstructure and structural connectivity in the brain but is often limited by the lengthy scan time needed to sample the diffusion profile by acquiring multiple diffusion weighted images (DWIs). Although parallel imaging technique has improved the speed of dMRI acquisition, attaining high resolution three dimensional (3D) dMRI on preclinical MRI systems remained still time consuming. In this paper, kernel principal component analysis, a machine learning approach, was employed to estimate the correlation among DWIs. We demonstrated the feasibility of such correlation estimation from low-resolution training DWIs and used the correlation as a constraint to reconstruct high-resolution DWIs from highly under-sampled k-space data, which significantly reduced the scan time. Using full k-space 3D dMRI data of post-mortem mouse brains, we retrospectively compared the performance of the so-called kernel low rank (KLR) method with a conventional compressed sensing (CS) method in terms of image quality and ability to resolve complex fiber orientations and connectivity. The results demonstrated that the KLR-CS method outperformed the conventional CS method for acceleration factors up to 8 and was likely to enhance our ability to investigate brain microstructure and connectivity using high-resolution 3D dMRI.
Collapse
Affiliation(s)
- Chaoyi Zhang
- Electrical Engineering, University at Buffalo, State University of New York, Buffalo, NY, United States
| | - Tanzil Mahmud Arefin
- Radiology, New York University School of Medicine, New York City, NY, United States
| | - Ukash Nakarmi
- Radiology, Stanford University, Stanford, CA, United States
| | - Choong Heon Lee
- Radiology, New York University School of Medicine, New York City, NY, United States
| | - Hongyu Li
- Electrical Engineering, University at Buffalo, State University of New York, Buffalo, NY, United States
| | - Dong Liang
- Paul C. Lauterbur Research Center for Biomedical Imaging, Shenzhen Institutes of Advanced Technology, Shenzhen, People's Republic of China
| | - Jiangyang Zhang
- Radiology, New York University School of Medicine, New York City, NY, United States
| | - Leslie Ying
- Electrical Engineering, University at Buffalo, State University of New York, Buffalo, NY, United States; Biomedical Engineering, University at Buffalo, State University at New York, Buffalo, NY, United States.
| |
Collapse
|
9
|
Lingala SG, Guo Y, Bliesener Y, Zhu Y, Lebel RM, Law M, Nayak KS. Tracer kinetic models as temporal constraints during brain tumor DCE-MRI reconstruction. Med Phys 2019; 47:37-51. [PMID: 31663134 PMCID: PMC6980286 DOI: 10.1002/mp.13885] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2019] [Revised: 10/17/2019] [Accepted: 10/17/2019] [Indexed: 12/11/2022] Open
Abstract
Purpose To apply tracer kinetic models as temporal constraints during reconstruction of under‐sampled brain tumor dynamic contrast enhanced (DCE) magnetic resonance imaging (MRI). Methods A library of concentration vs time profiles is simulated for a range of physiological kinetic parameters. The library is reduced to a dictionary of temporal bases, where each profile is approximated by a sparse linear combination of the bases. Image reconstruction is formulated as estimation of concentration profiles and sparse model coefficients with a fixed sparsity level. Simulations are performed to evaluate modeling error, and error statistics in kinetic parameter estimation in presence of noise. Retrospective under‐sampling experiments are performed on a brain tumor DCE digital reference object (DRO), and 12 brain tumor in‐vivo 3T datasets. The performances of the proposed under‐sampled reconstruction scheme and an existing compressed sensing‐based temporal finite‐difference (tFD) under‐sampled reconstruction were compared against the fully sampled inverse Fourier Transform‐based reconstruction. Results Simulations demonstrate that sparsity levels of 2 and 3 model the library profiles from the Patlak and extended Tofts‐Kety (ETK) models, respectively. Noise sensitivity analysis showed equivalent kinetic parameter estimation error statistics from noisy concentration profiles, and model approximated profiles. DRO‐based experiments showed good fidelity in recovery of kinetic maps from 20‐fold under‐sampled data. In‐vivo experiments demonstrated reduced bias and uncertainty in kinetic mapping with the proposed approach compared to tFD at under‐sampled reduction factors >= 20. Conclusions Tracer kinetic models can be applied as temporal constraints during brain tumor DCE‐MRI reconstruction. The proposed under‐sampled scheme resulted in model parameter estimates less biased with respect to conventional fully sampled DCE MRI reconstructions and parameter estimation. The approach is flexible, can use nonlinear kinetic models, and does not require tuning of regularization parameters.
Collapse
Affiliation(s)
- Sajan Goud Lingala
- Roy J Carver Department of Biomedical Engineering, University of Iowa, Iowa City, IA, USA
| | - Yi Guo
- Snap Inc., San Francisco, CA, USA
| | - Yannick Bliesener
- Ming Hsieh Department of Electrical and Computer Engineering, University of Southern California, Los Angeles, CA, USA
| | | | - R Marc Lebel
- GE Healthcare Applied Sciences Laboratory, Calgary, Canada
| | - Meng Law
- Department of Neuroscience, Monash University, Melbourne, Australia
| | - Krishna S Nayak
- Ming Hsieh Department of Electrical and Computer Engineering, University of Southern California, Los Angeles, CA, USA
| |
Collapse
|
10
|
Vemuri BC, Sun J, Banerjee M, Pan Z, Turner SM, Fuller DD, Forder JR, Entezari A. A geometric framework for ensemble average propagator reconstruction from diffusion MRI. Med Image Anal 2019; 57:89-105. [DOI: 10.1016/j.media.2019.06.012] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2018] [Revised: 03/13/2019] [Accepted: 06/24/2019] [Indexed: 01/10/2023]
|
11
|
Exploiting structural redundancy in q-space for improved EAP reconstruction from highly undersampled (k, q)-space in DMRI. Med Image Anal 2019; 54:122-137. [DOI: 10.1016/j.media.2019.02.014] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2018] [Revised: 11/22/2018] [Accepted: 02/26/2019] [Indexed: 11/22/2022]
|
12
|
Ye C, Prince JL. Dictionary-based fiber orientation estimation with improved spatial consistency. Med Image Anal 2018; 44:41-53. [PMID: 29190575 PMCID: PMC5771867 DOI: 10.1016/j.media.2017.11.010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2017] [Revised: 11/19/2017] [Accepted: 11/22/2017] [Indexed: 10/18/2022]
Abstract
Diffusion magnetic resonance imaging (dMRI) has enabled in vivo investigation of white matter tracts. Fiber orientation (FO) estimation is a key step in tract reconstruction and has been a popular research topic in dMRI analysis. In particular, the sparsity assumption has been used in conjunction with a dictionary-based framework to achieve reliable FO estimation with a reduced number of gradient directions. Because image noise can have a deleterious effect on the accuracy of FO estimation, previous works have incorporated spatial consistency of FOs in the dictionary-based framework to improve the estimation. However, because FOs are only indirectly determined from the mixture fractions of dictionary atoms and not modeled as variables in the objective function, these methods do not incorporate FO smoothness directly, and their ability to produce smooth FOs could be limited. In this work, we propose an improvement to Fiber Orientation Reconstruction using Neighborhood Information (FORNI), which we call FORNI+; this method estimates FOs in a dictionary-based framework where FO smoothness is better enforced than in FORNI alone. We describe an objective function that explicitly models the actual FOs and the mixture fractions of dictionary atoms. Specifically, it consists of data fidelity between the observed signals and the signals represented by the dictionary, pairwise FO dissimilarity that encourages FO smoothness, and weighted ℓ1-norm terms that ensure the consistency between the actual FOs and the FO configuration suggested by the dictionary representation. The FOs and mixture fractions are then jointly estimated by minimizing the objective function using an iterative alternating optimization strategy. FORNI+ was evaluated on a simulation phantom, a physical phantom, and real brain dMRI data. In particular, in the real brain dMRI experiment, we have qualitatively and quantitatively evaluated the reproducibility of the proposed method. Results demonstrate that FORNI+ produces FOs with better quality compared with competing methods.
Collapse
Affiliation(s)
- Chuyang Ye
- Brainnetome Center, Institute of Automation, Chinese Academy of Sciences, Beijing, China; National Laboratory of Pattern Recognition, Institute of Automation, Chinese Academy of Sciences, Beijing, China.
| | - Jerry L Prince
- Department of Electrical and Computer Engineering, Johns Hopkins University, Baltimore, MD, USA
| |
Collapse
|
13
|
Fang L, Li S, Kang X, Izatt JA, Farsiu S. 3-D Adaptive Sparsity Based Image Compression With Applications to Optical Coherence Tomography. IEEE TRANSACTIONS ON MEDICAL IMAGING 2015; 34:1306-20. [PMID: 25561591 PMCID: PMC4490145 DOI: 10.1109/tmi.2014.2387336] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/18/2023]
Abstract
We present a novel general-purpose compression method for tomographic images, termed 3D adaptive sparse representation based compression (3D-ASRC). In this paper, we focus on applications of 3D-ASRC for the compression of ophthalmic 3D optical coherence tomography (OCT) images. The 3D-ASRC algorithm exploits correlations among adjacent OCT images to improve compression performance, yet is sensitive to preserving their differences. Due to the inherent denoising mechanism of the sparsity based 3D-ASRC, the quality of the compressed images are often better than the raw images they are based on. Experiments on clinical-grade retinal OCT images demonstrate the superiority of the proposed 3D-ASRC over other well-known compression methods.
Collapse
Affiliation(s)
- Leyuan Fang
- College of Electrical and Information Engineering, Hunan University, Changsha, 410082, China
| | - Shutao Li
- College of Electrical and Information Engineering, Hunan University, Changsha, 410082, China
| | - Xudong Kang
- College of Electrical and Information Engineering, Hunan University, Changsha, 410082, China
| | - Joseph A. Izatt
- Biomedical Engineering, Duke University, Durham, NC 27708 USA and also with the Department of Ophthalmology, Duke University Medical Center, Durham, NC 27710 USA
| | - Sina Farsiu
- Departments of Biomedical Engineering, and Electrical and Computer Engineering, and Computer Science Duke University, Durham, NC 27708 USA , and also with the Department of Ophthalmology, Duke University Medical Center, Durham, NC 27710 USA
| |
Collapse
|
14
|
Knoll F, Raya JG, Halloran RO, Baete S, Sigmund E, Bammer R, Block T, Otazo R, Sodickson DK. A model-based reconstruction for undersampled radial spin-echo DTI with variational penalties on the diffusion tensor. NMR IN BIOMEDICINE 2015; 28:353-66. [PMID: 25594167 PMCID: PMC4339452 DOI: 10.1002/nbm.3258] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/09/2014] [Revised: 12/08/2014] [Accepted: 12/17/2014] [Indexed: 05/04/2023]
Abstract
Radial spin-echo diffusion imaging allows motion-robust imaging of tissues with very low T2 values like articular cartilage with high spatial resolution and signal-to-noise ratio (SNR). However, in vivo measurements are challenging, due to the significantly slower data acquisition speed of spin-echo sequences and the less efficient k-space coverage of radial sampling, which raises the demand for accelerated protocols by means of undersampling. This work introduces a new reconstruction approach for undersampled diffusion-tensor imaging (DTI). A model-based reconstruction implicitly exploits redundancies in the diffusion-weighted images by reducing the number of unknowns in the optimization problem and compressed sensing is performed directly in the target quantitative domain by imposing a total variation (TV) constraint on the elements of the diffusion tensor. Experiments were performed for an anisotropic phantom and the knee and brain of healthy volunteers (three and two volunteers, respectively). Evaluation of the new approach was conducted by comparing the results with reconstructions performed with gridding, combined parallel imaging and compressed sensing and a recently proposed model-based approach. The experiments demonstrated improvements in terms of reduction of noise and streaking artifacts in the quantitative parameter maps, as well as a reduction of angular dispersion of the primary eigenvector when using the proposed method, without introducing systematic errors into the maps. This may enable an essential reduction of the acquisition time in radial spin-echo diffusion-tensor imaging without degrading parameter quantification and/or SNR.
Collapse
Affiliation(s)
- Florian Knoll
- Bernard and Irene Schwartz Center for Biomedical Imaging, New York University School of Medicine, New York, New York, USA
- Correspondence to: Florian Knoll, PhD, New York University School of Medicine, Center for Biomedical Imaging, 660 First Avenue, 4th Floor, New York, NY 10016, Phone: 212-263-0335,
| | - José G Raya
- Bernard and Irene Schwartz Center for Biomedical Imaging, New York University School of Medicine, New York, New York, USA
| | - Rafael O Halloran
- Department of Radiology, Stanford University, Stanford, California, USA
| | - Steven Baete
- Bernard and Irene Schwartz Center for Biomedical Imaging, New York University School of Medicine, New York, New York, USA
| | - Eric Sigmund
- Bernard and Irene Schwartz Center for Biomedical Imaging, New York University School of Medicine, New York, New York, USA
| | - Roland Bammer
- Department of Radiology, Stanford University, Stanford, California, USA
| | - Tobias Block
- Bernard and Irene Schwartz Center for Biomedical Imaging, New York University School of Medicine, New York, New York, USA
| | - Ricardo Otazo
- Bernard and Irene Schwartz Center for Biomedical Imaging, New York University School of Medicine, New York, New York, USA
| | - Daniel K Sodickson
- Bernard and Irene Schwartz Center for Biomedical Imaging, New York University School of Medicine, New York, New York, USA
| |
Collapse
|
15
|
Paquette M, Merlet S, Gilbert G, Deriche R, Descoteaux M. Comparison of sampling strategies and sparsifying transforms to improve compressed sensing diffusion spectrum imaging. Magn Reson Med 2014; 73:401-16. [PMID: 24478106 DOI: 10.1002/mrm.25093] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2013] [Revised: 11/21/2013] [Accepted: 12/02/2013] [Indexed: 11/08/2022]
Abstract
PURPOSE Diffusion Spectrum Imaging enables to reconstruct the ensemble average propagator (EAP) at the expense of having to acquire a large number of measurements. Compressive sensing offers an efficient way to decrease the required number of measurements. The purpose of this work is to perform a thorough experimental comparison of three sampling strategies and six sparsifying transforms to show their impact when applied to accelerate compressive sensing-diffusion spectrum imaging. METHODS We propose a novel sampling scheme that assures uniform angular and random radial q-space samples. We also compare and implement six discrete sparse representations of the EAP and thoroughly evaluate them on synthetic and real data using metrics from the full EAP, kurtosis, and orientation distribution function. RESULTS The discrete wavelet transform with Cohen-Daubechies-Feauveau 9/7 wavelets and uniform angular sampling in combination with random radial sampling showed to be better than other tested techniques to accurately reconstruct the EAP and its features. CONCLUSION It is important to jointly optimize the sampling scheme and the sparsifying transform to obtain accelerated compressive sensing-diffusion spectrum imaging. Experiments on synthetic and real human brain data show that one can robustly recover both radial and angular EAP features while undersampling the acquisition to 64 measurements (undersampling factor of 4).
Collapse
Affiliation(s)
- Michael Paquette
- Department of Computer Science, Sherbrooke Connectivity Imaging Laboratory, Université de Sherbrooke, Sherbrooke, Quebec, Canada
| | - Sylvain Merlet
- Athena Project-Team, INRIA Sophia Antipolis-Méditerranée, Sophia-Antipolis Cedex, France
| | | | - Rachid Deriche
- Athena Project-Team, INRIA Sophia Antipolis-Méditerranée, Sophia-Antipolis Cedex, France
| | - Maxime Descoteaux
- Department of Computer Science, Sherbrooke Connectivity Imaging Laboratory, Université de Sherbrooke, Sherbrooke, Quebec, Canada
| |
Collapse
|