1
|
Pala S, Paajanen A, Ristaniemi A, Nippolainen E, Afara IO, Nykänen O, Nissi MJ. Measurement of T 1ρ dispersion with compressed sensing and magnetization prepared radial balanced steady-state free precession in spontaneous human osteoarthritis. Magn Reson Med 2024; 92:2127-2139. [PMID: 38953429 DOI: 10.1002/mrm.30206] [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: 01/25/2024] [Revised: 06/14/2024] [Accepted: 06/14/2024] [Indexed: 07/04/2024]
Abstract
PURPOSE To assess the potential for accelerating continuous-wave (CW) T1ρ dispersion measurement with compressed sensing approach via studying the effect that the data reduction has on the ability to detect differences between intact and degenerated articular cartilage with different spin-lock amplitudes and to assess quantitative bias due to acceleration. METHODS Osteochondral plugs (n = 27, 4 mm diameter) from femur (n = 14) and tibia (n = 13) regions from human cadaver knee joints were obtained from commercial biobank (Science Care, USA) under Ethical permission 134/2015. MRI of specimens was performed at 9.4T with magnetization prepared radial balanced SSFP (bSSFP) readout sequence, and the CWT1ρ relaxation time maps were computed from the measured data. The relaxation time maps were evaluated in the cartilage zones for different acceleration factors. For reference, Osteoarthritis Research Society International (OARSI) grading and biomechanical measurements were performed and correlated with the MRI findings. RESULTS Four-fold acceleration of CWT1ρ dispersion measurement by compressed sensing approach was feasible without meaningful loss in the sensitivity to osteoarthritic (OA) changes within the articular cartilage. Differences were significant between intact and OA groups in the superficial and transitional zones, and CWT1ρ correlated moderately with the reference measurements (0.3 < r < 0.7). CONCLUSION CWT1ρ was able to differentiate between intact and OA cartilage even with four-fold acceleration. This indicates that acceleration of CWT1ρ dispersion measurement by compressed sensing approach is feasible with negligible loss in the sensitivity to osteoarthritic changes in articular cartilage.
Collapse
Grants
- 780598 Horizon 2020 Framework Programme
- H2020-ICT-2017-1 Horizon 2020 Framework Programme
- 358944 Research Council of Finland (Flagship of Advanced Mathematics for Sensing, Imaging and Modelling Grant)
- 315820 Research Council of Finland (Flagship of Advanced Mathematics for Sensing, Imaging and Modelling Grant)
- 324529 Research Council of Finland (Flagship of Advanced Mathematics for Sensing, Imaging and Modelling Grant)
- 324994 Research Council of Finland (Flagship of Advanced Mathematics for Sensing, Imaging and Modelling Grant)
- 325146 Research Council of Finland (Flagship of Advanced Mathematics for Sensing, Imaging and Modelling Grant)
- 348410 Research Council of Finland (Flagship of Advanced Mathematics for Sensing, Imaging and Modelling Grant)
- 352666 Research Council of Finland (Flagship of Advanced Mathematics for Sensing, Imaging and Modelling Grant)
- 354693 Research Council of Finland (Flagship of Advanced Mathematics for Sensing, Imaging and Modelling Grant)
- 357787 Research Council of Finland (Flagship of Advanced Mathematics for Sensing, Imaging and Modelling Grant)
- 240130 Sigrid Jusélius Foundation
- Olvi Foundation
- Päivikki and Sakari Sohlberg Foundation
- Instrumentarium Science foundation
- 65231459 Finnish Cultural Foundation, North-Savonia Regional Fund
Collapse
Affiliation(s)
- S Pala
- Department of Technical Physics, University of Eastern Finland, Kuopio, Finland
| | - A Paajanen
- Department of Technical Physics, University of Eastern Finland, Kuopio, Finland
| | - A Ristaniemi
- Department of Technical Physics, University of Eastern Finland, Kuopio, Finland
| | - E Nippolainen
- Department of Technical Physics, University of Eastern Finland, Kuopio, Finland
| | - I O Afara
- Department of Technical Physics, University of Eastern Finland, Kuopio, Finland
| | - O Nykänen
- Department of Technical Physics, University of Eastern Finland, Kuopio, Finland
| | - M J Nissi
- Department of Technical Physics, University of Eastern Finland, Kuopio, Finland
| |
Collapse
|
2
|
Plummer JW, Hussain R, Bdaiwi AS, Soderlund SA, Hoyos X, Lanier JM, Garrison WJ, Parra-Robles J, Willmering MM, Niedbalski P, Cleveland ZI, Walkup L. A decay-modeled compressed sensing reconstruction approach for non-Cartesian hyperpolarized 129Xe MRI. Magn Reson Med 2024; 92:1363-1375. [PMID: 38860514 PMCID: PMC11262966 DOI: 10.1002/mrm.30188] [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: 02/14/2024] [Revised: 04/15/2024] [Accepted: 05/18/2024] [Indexed: 06/12/2024]
Abstract
PURPOSE Hyperpolarized 129Xe MRI benefits from non-Cartesian acquisitions that sample k-space efficiently and rapidly. However, their reconstructions are complex and burdened by decay processes unique to hyperpolarized gas. Currently used gridded reconstructions are prone to artifacts caused by magnetization decay and are ill-suited for undersampling. We present a compressed sensing (CS) reconstruction approach that incorporates magnetization decay in the forward model, thereby producing images with increased sharpness and contrast, even in undersampled data. METHODS Radio-frequency, T1, andT 2 * $$ {\mathrm{T}}_2^{\ast } $$ decay processes were incorporated into the forward model and solved using iterative methods including CS. The decay-modeled reconstruction was validated in simulations and then tested in 2D/3D-spiral ventilation and 3D-radial gas-exchange MRI. Quantitative metrics including apparent-SNR and sharpness were compared between gridded, CS, and twofold undersampled CS reconstructions. Observations were validated in gas-exchange data collected from 15 healthy and 25 post-hematopoietic-stem-cell-transplant participants. RESULTS CS reconstructions in simulations yielded images with threefold increases in accuracy. CS increased sharpness and contrast for ventilation in vivo imaging and showed greater accuracy for undersampled acquisitions. CS improved gas-exchange imaging, particularly in the dissolved-phase where apparent-SNR improved, and structure was made discernable. Finally, CS showed repeatability in important global gas-exchange metrics including median dissolved-gas signal ratio and median angle between real/imaginary components. CONCLUSION A non-Cartesian CS reconstruction approach that incorporates hyperpolarized 129Xe decay processes is presented. This approach enables improved image sharpness, contrast, and overall image quality in addition to up-to threefold undersampling. This contribution benefits all hyperpolarized gas MRI through improved accuracy and decreased scan durations.
Collapse
Affiliation(s)
- J. W. Plummer
- Department of Biomedical Engineering, University of Cincinnati, Cincinnati, OH, United States
- Center for Pulmonary Imaging Research, Department of Pulmonary Medicine, Cincinnati Children’s Hospital Medical Center, Cincinnati, OH, United States
| | - R. Hussain
- Center for Pulmonary Imaging Research, Department of Pulmonary Medicine, Cincinnati Children’s Hospital Medical Center, Cincinnati, OH, United States
| | - A. S. Bdaiwi
- Department of Biomedical Engineering, University of Cincinnati, Cincinnati, OH, United States
- Center for Pulmonary Imaging Research, Department of Pulmonary Medicine, Cincinnati Children’s Hospital Medical Center, Cincinnati, OH, United States
| | - S. A. Soderlund
- Department of Biomedical Engineering, University of Cincinnati, Cincinnati, OH, United States
- Center for Pulmonary Imaging Research, Department of Pulmonary Medicine, Cincinnati Children’s Hospital Medical Center, Cincinnati, OH, United States
| | - X. Hoyos
- Department of Biomedical Engineering, University of Cincinnati, Cincinnati, OH, United States
- Center for Pulmonary Imaging Research, Department of Pulmonary Medicine, Cincinnati Children’s Hospital Medical Center, Cincinnati, OH, United States
| | - J. M. Lanier
- Department of Radiology, Cincinnati Children’s Hospital Medical Center, Cincinnati, OH, United States
| | - W. J. Garrison
- Center for Pulmonary Imaging Research, Department of Pulmonary Medicine, Cincinnati Children’s Hospital Medical Center, Cincinnati, OH, United States
| | - J. Parra-Robles
- Center for Pulmonary Imaging Research, Department of Pulmonary Medicine, Cincinnati Children’s Hospital Medical Center, Cincinnati, OH, United States
| | - M. M. Willmering
- Center for Pulmonary Imaging Research, Department of Pulmonary Medicine, Cincinnati Children’s Hospital Medical Center, Cincinnati, OH, United States
| | - P.J. Niedbalski
- Pulmonary, Critical Care and Sleep Medicine, Kansas University Medical Center, KS, United States
| | - Z. I. Cleveland
- Department of Biomedical Engineering, University of Cincinnati, Cincinnati, OH, United States
- Center for Pulmonary Imaging Research, Department of Pulmonary Medicine, Cincinnati Children’s Hospital Medical Center, Cincinnati, OH, United States
- Department of Radiology, Cincinnati Children’s Hospital Medical Center, Cincinnati, OH, United States
- Department of Pediatrics, University of Cincinnati, Cincinnati, OH, United States
| | - L.L. Walkup
- Department of Biomedical Engineering, University of Cincinnati, Cincinnati, OH, United States
- Center for Pulmonary Imaging Research, Department of Pulmonary Medicine, Cincinnati Children’s Hospital Medical Center, Cincinnati, OH, United States
- Department of Radiology, Cincinnati Children’s Hospital Medical Center, Cincinnati, OH, United States
- Department of Pediatrics, University of Cincinnati, Cincinnati, OH, United States
| |
Collapse
|
3
|
Qu B, Zhang J, Kang T, Lin J, Lin M, She H, Wu Q, Wang M, Zheng G. Radial magnetic resonance image reconstruction with a deep unrolled projected fast iterative soft-thresholding network. Comput Biol Med 2024; 168:107707. [PMID: 38000244 DOI: 10.1016/j.compbiomed.2023.107707] [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: 09/24/2023] [Revised: 10/31/2023] [Accepted: 11/15/2023] [Indexed: 11/26/2023]
Abstract
Radially sampling of magnetic resonance imaging (MRI) is an effective way to accelerate the imaging. How to preserve the image details in reconstruction is always challenging. In this work, a deep unrolled neural network is designed to emulate the iterative sparse image reconstruction process of a projected fast soft-threshold algorithm (pFISTA). The proposed method, an unrolled pFISTA network for Deep Radial MRI (pFISTA-DR), include the preprocessing module to refine coil sensitivity maps and initial reconstructed image, the learnable convolution filters to extract image feature maps, and adaptive threshold to robustly remove image artifacts. Experimental results show that, among the compared methods, pFISTA-DR provides the best reconstruction and achieved the highest PSNR, the highest SSIM and the lowest reconstruction errors.
Collapse
Affiliation(s)
- Biao Qu
- Pen-Tung Sah Institute of Micro-Nano Science and Technology, Xiamen University, Xiamen, China
| | - Jialue Zhang
- Pen-Tung Sah Institute of Micro-Nano Science and Technology, Xiamen University, Xiamen, China; Fujian Provincial Key Laboratory of Plasma and Magnetic Resonance, Department of Electronic Science, Xiamen University, China
| | - Taishan Kang
- Department of Radiology, Zhongshan Hospital of Xiamen University, School of Medicine, Xiamen University, Xiamen, China
| | - Jianzhong Lin
- Department of Radiology, Zhongshan Hospital of Xiamen University, School of Medicine, Xiamen University, Xiamen, China
| | - Meijin Lin
- Department of Applied Marine Physics & Engineering, College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
| | - Huajun She
- School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai, China
| | - Qingxia Wu
- Department of Medical Imaging, Henan Provincial People's Hospital, Zhengzhou, China
| | - Meiyun Wang
- Department of Medical Imaging, Henan Provincial People's Hospital, Zhengzhou, China; Laboratory of Brain Science and Brain-Like Intelligence Technology, Institute for Integrated Medical Science and Engineering, Henan Academy of Sciences, Zhengzhou, China
| | - Gaofeng Zheng
- Pen-Tung Sah Institute of Micro-Nano Science and Technology, Xiamen University, Xiamen, China.
| |
Collapse
|
4
|
Paajanen A, Hanhela M, Hänninen N, Nykänen O, Kolehmainen V, Nissi MJ. Fast Compressed Sensing of 3D Radial T 1 Mapping with Different Sparse and Low-Rank Models. J Imaging 2023; 9:151. [PMID: 37623683 PMCID: PMC10455972 DOI: 10.3390/jimaging9080151] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2023] [Revised: 07/14/2023] [Accepted: 07/21/2023] [Indexed: 08/26/2023] Open
Abstract
Knowledge of the relative performance of the well-known sparse and low-rank compressed sensing models with 3D radial quantitative magnetic resonance imaging acquisitions is limited. We use 3D radial T1 relaxation time mapping data to compare the total variation, low-rank, and Huber penalty function approaches to regularization to provide insights into the relative performance of these image reconstruction models. Simulation and ex vivo specimen data were used to determine the best compressed sensing model as measured by normalized root mean squared error and structural similarity index. The large-scale compressed sensing models were solved by combining a GPU implementation of a preconditioned primal-dual proximal splitting algorithm to provide high-quality T1 maps within a feasible computation time. The model combining spatial total variation and locally low-rank regularization yielded the best performance, followed closely by the model combining spatial and contrast dimension total variation. Computation times ranged from 2 to 113 min, with the low-rank approaches taking the most time. The differences between the compressed sensing models are not necessarily large, but the overall performance is heavily dependent on the imaged object.
Collapse
Affiliation(s)
| | | | | | | | | | - Mikko J. Nissi
- Department of Technical Physics, University of Eastern Finland, 70211 Kuopio, Finland; (A.P.); (M.H.); (N.H.); (O.N.); (V.K.)
| |
Collapse
|
5
|
Qu B, Zhang Z, Chen Y, Qian C, Kang T, Lin J, Chen L, Wu Z, Wang J, Zheng G, Qu X. A convergence analysis for projected fast iterative soft-thresholding algorithm under radial sampling MRI. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2023; 351:107425. [PMID: 37060889 DOI: 10.1016/j.jmr.2023.107425] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/30/2022] [Revised: 02/27/2023] [Accepted: 03/17/2023] [Indexed: 05/29/2023]
Abstract
Radial sampling is a fast magnetic resonance imaging technique. Further imaging acceleration can be achieved with undersampling but how to reconstruct a clear image with fast algorithm is still challenging. Previous work has shown the advantage of removing undersampling image artifacts using the tight-frame sparse reconstruction model. This model was further solved with a projected fast iterative soft-thresholding algorithm (pFISTA). However, the convergence of this algorithm under radial sampling has not been clearly set up. In this work, the authors derived a theoretical convergence condition for this algorithm. This condition was approximated by estimating the maximal eigenvalue of reconstruction operators through the power iteration. Based on the condition, an optimal step size was further suggested to allow the fastest convergence. Verifications were made on the prospective in vivo data of static brain imaging and dynamic contrast-enhanced liver imaging, demonstrating that the recommended parameter allowed fast convergence in radial MRI.
Collapse
Affiliation(s)
- Biao Qu
- Department of Instrumental and Electrical Engineering, Xiamen University, Xiamen, China
| | - Zuwen Zhang
- Department of Electronic Science, Biomedical Intelligent Cloud R&D Center, Fujian Provincial Key Laboratory of Plasma and Magnetic Resonance, Xiamen University, Xiamen, China
| | - Yewei Chen
- Department of Electronic Science, Biomedical Intelligent Cloud R&D Center, Fujian Provincial Key Laboratory of Plasma and Magnetic Resonance, Xiamen University, Xiamen, China
| | - Chen Qian
- Department of Electronic Science, Biomedical Intelligent Cloud R&D Center, Fujian Provincial Key Laboratory of Plasma and Magnetic Resonance, Xiamen University, Xiamen, China
| | - Taishan Kang
- Department of Radiology, Zhongshan Hospital of Xiamen University, School of Medicine, Xiamen University, Xiamen, China
| | - Jianzhong Lin
- Department of Radiology, Zhongshan Hospital of Xiamen University, School of Medicine, Xiamen University, Xiamen, China
| | - Lihua Chen
- Department of Radiology, Chongqing University Cancer Hospital, School of Medicine, Chongqing University, Chongqing, China
| | | | | | - Gaofeng Zheng
- Department of Instrumental and Electrical Engineering, Xiamen University, Xiamen, China.
| | - Xiaobo Qu
- Department of Electronic Science, Biomedical Intelligent Cloud R&D Center, Fujian Provincial Key Laboratory of Plasma and Magnetic Resonance, Xiamen University, Xiamen, China.
| |
Collapse
|
6
|
Nagtegaal M, Hartsema E, Koolstra K, Vos F. Multicomponent MR fingerprinting reconstruction using joint-sparsity and low-rank constraints. Magn Reson Med 2023; 89:286-298. [PMID: 36121015 PMCID: PMC9825911 DOI: 10.1002/mrm.29442] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2022] [Revised: 08/12/2022] [Accepted: 08/12/2022] [Indexed: 01/11/2023]
Abstract
PURPOSE To develop an efficient algorithm for multicomponent MR fingerprinting (MC-MRF) reconstructions directly from highly undersampled data without making prior assumptions about tissue relaxation times and expected number of tissues. METHODS The proposed method reconstructs MC-MRF maps from highly undersampled data by iteratively applying a joint-sparsity constraint to the estimated tissue components. Intermediate component maps are obtained by a low-rank multicomponent alternating direction method of multipliers (MC-ADMM) including the non-negativity of tissue weights as an extra regularization term. Over iterations, the used dictionary compression is adjusted. The proposed method (k-SPIJN) is compared with a two-step approach in which image reconstruction and multicomponent estimations are performed sequentially and tested in numerical simulations and in vivo by applying different undersampling factors in eight healthy volunteers. In the latter case, fully sampled data serves as the reference. RESULTS The proposed method shows improved precision and accuracy in simulations compared with a state-of-art sequential approach. Obtained in vivo magnetization fraction maps for different tissue types show reduced systematic errors and reduced noise-like effects. Root mean square errors in estimated magnetization fraction maps significantly reduce from 13.0% <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mo>±</mml:mo></mml:mrow> <mml:annotation>$$ \pm $$</mml:annotation></mml:semantics> </mml:math> 5.8% with the conventional, two-step approach to 9.6% <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mo>±</mml:mo></mml:mrow> <mml:annotation>$$ \pm $$</mml:annotation></mml:semantics> </mml:math> 3.9% and 9.6% <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mo>±</mml:mo></mml:mrow> <mml:annotation>$$ \pm $$</mml:annotation></mml:semantics> </mml:math> 3.2% with the proposed MC-ADMM and k-SPIJN methods, respectively. Mean standard deviation in homogeneous white matter regions reduced significantly from 8.6% to 2.9% (two step vs. k-SPIJN). CONCLUSION The proposed MC-ADMM and k-SPIJN reconstruction methods estimate MC-MRF maps from highly undersampled data resulting in improved image quality compared with the existing method.
Collapse
Affiliation(s)
- Martijn Nagtegaal
- Department of Imaging PhysicsDelft University of TechnologyDelftThe Netherlands
| | - Emiel Hartsema
- Department of Imaging PhysicsDelft University of TechnologyDelftThe Netherlands
| | - Kirsten Koolstra
- Division of Image Processing, Department of RadiologyLeiden University Medical CenterLeidenThe Netherlands
| | - Frans Vos
- Department of Imaging PhysicsDelft University of TechnologyDelftThe Netherlands,Department of RadiologyErasmus MCRotterdamThe Netherlands
| |
Collapse
|
7
|
Ramzi Z, G R C, Starck JL, Ciuciu P. NC-PDNet: A Density-Compensated Unrolled Network for 2D and 3D Non-Cartesian MRI Reconstruction. IEEE TRANSACTIONS ON MEDICAL IMAGING 2022; 41:1625-1638. [PMID: 35041598 DOI: 10.1109/tmi.2022.3144619] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Deep Learning has become a very promising avenue for magnetic resonance image (MRI) reconstruction. In this work, we explore the potential of unrolled networks for non-Cartesian acquisition settings. We design the NC-PDNet (Non-Cartesian Primal Dual Netwok), the first density-compensated (DCp) unrolled neural network, and validate the need for its key components via an ablation study. Moreover, we conduct some generalizability experiments to test this network in out-of-distribution settings, for example training on knee data and validating on brain data. The results show that NC-PDNet outperforms baseline (U-Net, Deep image prior) models both visually and quantitatively in all settings. In particular, in the 2D multi-coil acquisition scenario, the NC-PDNet provides up to a 1.2 dB improvement in peak signal-to-noise ratio (PSNR) over baseline networks, while also allowing a gain of at least 1dB in PSNR in generalization settings. We provide the open-source implementation of NC-PDNet, and in particular the Non-uniform Fourier Transform in TensorFlow, tested on 2D multi-coil and 3D single-coil k-space data.
Collapse
|
8
|
Chen Z, Chen Y, Xie Y, Li D, Christodoulou AG. Data-Consistent non-Cartesian deep subspace learning for efficient dynamic MR image reconstruction. PROCEEDINGS. IEEE INTERNATIONAL SYMPOSIUM ON BIOMEDICAL IMAGING 2022; 2022:10.1109/isbi52829.2022.9761497. [PMID: 35572068 PMCID: PMC9104888 DOI: 10.1109/isbi52829.2022.9761497] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
Non-Cartesian sampling with subspace-constrained image reconstruction is a popular approach to dynamic MRI, but slow iterative reconstruction limits its clinical application. Data-consistent (DC) deep learning can accelerate reconstruction with good image quality, but has not been formulated for non-Cartesian subspace imaging. In this study, we propose a DC non-Cartesian deep subspace learning framework for fast, accurate dynamic MR image reconstruction. Four novel DC formulations are developed and evaluated: two gradient decent approaches, a directly solved approach, and a conjugate gradient approach. We applied a U-Net model with and without DC layers to reconstruct T1-weighted images for cardiac MR Multitasking (an advanced multidimensional imaging method), comparing our results to the iteratively reconstructed reference. Experimental results show that the proposed framework significantly improves reconstruction accuracy over the U-Net model without DC, while significantly accelerating the reconstruction over conventional iterative reconstruction.
Collapse
Affiliation(s)
- Zihao Chen
- Biomedical Imaging Research Institute, Cedars-Sinai Medical Center, Los Angeles, USA
- Department of Bioengineering, UCLA, Los Angeles, USA
| | - Yuhua Chen
- Biomedical Imaging Research Institute, Cedars-Sinai Medical Center, Los Angeles, USA
- Department of Bioengineering, UCLA, Los Angeles, USA
| | - Yibin Xie
- Biomedical Imaging Research Institute, Cedars-Sinai Medical Center, Los Angeles, USA
| | - Debiao Li
- Biomedical Imaging Research Institute, Cedars-Sinai Medical Center, Los Angeles, USA
- Department of Bioengineering, UCLA, Los Angeles, USA
| | - Anthony G Christodoulou
- Biomedical Imaging Research Institute, Cedars-Sinai Medical Center, Los Angeles, USA
- Department of Bioengineering, UCLA, Los Angeles, USA
| |
Collapse
|
9
|
Koolstra K, Remis R. Learning a preconditioner to accelerate compressed sensing reconstructions in MRI. Magn Reson Med 2021; 87:2063-2073. [PMID: 34752655 PMCID: PMC9299023 DOI: 10.1002/mrm.29073] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2021] [Revised: 10/13/2021] [Accepted: 10/20/2021] [Indexed: 01/07/2023]
Abstract
Purpose To learn a preconditioner that accelerates parallel imaging (PI) and compressed sensing (CS) reconstructions. Methods A convolutional neural network (CNN) with residual connections was used to train a preconditioning operator. Training and validation data were simulated using 50% brain images and 50% white Gaussian noise images. Each multichannel training example contains a simulated sampling mask, complex coil sensitivity maps, and two regularization parameter maps. The trained model was integrated in the preconditioned conjugate gradient (PCG) method as part of the split Bregman CS method. The acceleration performance was compared with that of a circulant PI‐CS preconditioner for varying undersampling factors, number of coil elements and anatomies. Results The learned preconditioner reduces the number of PCG iterations by a factor of 4, yielding a similar acceleration as an efficient circulant preconditioner. The method generalizes well to different sampling schemes, coil configurations and anatomies. Conclusion It is possible to learn adaptable preconditioners for PI and CS reconstructions that meet the performance of state‐of‐the‐art preconditioners. Further acceleration could be achieved by optimizing the network architecture and the training set. Such a preconditioner could also be integrated in fully learned reconstruction methods to accelerate the training process of unrolled networks.
Collapse
Affiliation(s)
- Kirsten Koolstra
- Division of Image Processing, Department of Radiology, Leiden University Medical Center, Leiden, The Netherlands
| | - Rob Remis
- Circuits & Systems Group, Electrical Engineering, Mathematics and Computer Science Faculty, Delft University of Technology, Delft, The Netherlands
| |
Collapse
|
10
|
Hong K, Schiffers F, DiCarlo AL, Rigsby CK, Haji-Valizadeh H, Lee DC, Markl M, Katsaggelos AK, Kim D. Accelerating compressed sensing reconstruction of subsampled radial k-space data using geometrically-derived density compensation. Phys Med Biol 2021; 66. [PMID: 34607316 DOI: 10.1088/1361-6560/ac2c9d] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2021] [Accepted: 10/04/2021] [Indexed: 11/12/2022]
Abstract
Objective.To accelerate compressed sensing (CS) reconstruction of subsampled radial k-space data using a geometrically-derived density compensation function (gDCF) without significant loss in image quality.Approach.We developed a theoretical framework to calculate a gDCF based on Nyquist distance along the radial and circumferential directions of a discrete polar coordinate system. Our gDCF was compared against standard DCF (e.g. ramp filter) and another commonly used DCF (modified Shepp-Logan (SL) filter). The resulting image quality produced by each DCF was quantified using normalized root-mean-square-error (NRMSE), blur metric (1 = blurriest; 0 = sharpest), and structural similarity index (SSIM; 1 = perfect match; 0 = no match) compared with the reference. For filtered backprojection (FBP) of phantom data obtained at the Nyquist sampling rate, Cartesian k-space sampling was used as the reference. For CS reconstruction of subsampled cardiac magnetic resonance imaging datasets (real-time cardiac cine data with 11 projections per frame,n = 20 patients; cardiac perfusion data with 30 projections per frame,n = 19 patients), CS reconstruction without DCF was used as the reference.Main results.The NRMSE, SSIM, and blur metrics of the phantom data were good for all DCFs, confirming that our gDCF produces uniform densities at the upper limit (Nyquist). For CS reconstruction of subsampled real-time cine and cardiac perfusion datasets, the image quality metrics (SSIM, NRMSE) were significantly (p < 0.05) higher for our gDCF than other DCFs, and the reconstruction time was significantly (p < 0.05) faster for our gDCF (reference) than no DCF (11.9%-52.9% slower), standard DCF (23.9%-57.6% slower), and modified SL filter (13.5%-34.8% slower).Significance.The proposed gDCF accelerates CS reconstruction of subsampled radial k-space data without significant loss in image quality compared with no DCF as the reference.
Collapse
Affiliation(s)
- KyungPyo Hong
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, Illinois, United States of America
| | - Florian Schiffers
- Department of Computer Science, McCormick School of Engineering and Applied Science, Northwestern University, Evanston, Illinois, United States of America
| | - Amanda L DiCarlo
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, Illinois, United States of America
| | - Cynthia K Rigsby
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, Illinois, United States of America.,Department of Medical Imaging, Ann & Robert H. Lurie Children's Hospital of Chicago, Chicago, Illinois, United States of America
| | - Hassan Haji-Valizadeh
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, Illinois, United States of America
| | - Daniel C Lee
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, Illinois, United States of America.,Department of Medicine, Division of Cardiology, Northwestern University Feinberg School of Medicine, Chicago, Illinois, United States of America
| | - Michael Markl
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, Illinois, United States of America.,Department of Biomedical Engineering, McCormick School of Engineering, Northwestern University, Evanston, Illinois, United States of America
| | - Aggelos K Katsaggelos
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, Illinois, United States of America.,Department of Computer Science, McCormick School of Engineering and Applied Science, Northwestern University, Evanston, Illinois, United States of America.,Department of Electrical and Computer Engineering, McCormick School of Engineering and Applied Science, Northwestern University, Evanston, Illinois, United States of America
| | - Daniel Kim
- Department of Radiology, Northwestern University Feinberg School of Medicine, Chicago, Illinois, United States of America.,Department of Biomedical Engineering, McCormick School of Engineering, Northwestern University, Evanston, Illinois, United States of America
| |
Collapse
|
11
|
Kofler A, Haltmeier M, Schaeffter T, Kolbitsch C. An end-to-end-trainable iterative network architecture for accelerated radial multi-coil 2D cine MR image reconstruction. Med Phys 2021; 48:2412-2425. [PMID: 33651398 DOI: 10.1002/mp.14809] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2020] [Revised: 02/11/2021] [Accepted: 02/18/2021] [Indexed: 11/06/2022] Open
Abstract
PURPOSE Iterative convolutional neural networks (CNNs) which resemble unrolled learned iterative schemes have shown to consistently deliver state-of-the-art results for image reconstruction problems across different imaging modalities. However, because these methods include the forward model in the architecture, their applicability is often restricted to either relatively small reconstruction problems or to problems with operators which are computationally cheap to compute. As a consequence, they have not been applied to dynamic non-Cartesian multi-coil reconstruction problems so far. METHODS In this work, we propose a CNN architecture for image reconstruction of accelerated 2D radial cine MRI with multiple receiver coils. The network is based on a computationally light CNN component and a subsequent conjugate gradient (CG) method which can be jointly trained end-to-end using an efficient training strategy. We investigate the proposed training strategy and compare our method with other well-known reconstruction techniques with learned and non-learned regularization methods. RESULTS Our proposed method outperforms all other methods based on non-learned regularization. Further, it performs similar or better than a CNN-based method employing a 3D U-Net and a method using adaptive dictionary learning. In addition, we empirically demonstrate that even by training the network with only iteration, it is possible to increase the length of the network at test time and further improve the results. CONCLUSIONS End-to-end training allows to highly reduce the number of trainable parameters of and stabilize the reconstruction network. Further, because it is possible to change the length of the network at the test time, the need to find a compromise between the complexity of the CNN-block and the number of iterations in each CG-block becomes irrelevant.
Collapse
Affiliation(s)
- Andreas Kofler
- Physikalisch-Technische Bundesanstalt (PTB), Braunschweig, Berlin, 10587, Germany
| | - Markus Haltmeier
- Department of Mathematics, University of Innsbruck, Innsbruck, 6020, Austria
| | - Tobias Schaeffter
- Physikalisch-Technische Bundesanstalt (PTB), Braunschweig, Berlin, 10587, Germany.,School of Imaging Sciences and Biomedical Engineering, King's College London, London, SE1 7EH, UK.,Department of Biomedical Engineering, Technical University of Berlin, Berlin, 10623, Germany
| | - Christoph Kolbitsch
- Physikalisch-Technische Bundesanstalt (PTB), Braunschweig, Berlin, 10587, Germany.,School of Imaging Sciences and Biomedical Engineering, King's College London, London, SE1 7EH, UK
| |
Collapse
|
12
|
Saucedo A, Macey PM, Thomas MA. Accelerated radial echo-planar spectroscopic imaging using golden angle view-ordering and compressed-sensing reconstruction with total variation regularization. Magn Reson Med 2021; 86:46-61. [PMID: 33604944 DOI: 10.1002/mrm.28728] [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: 05/27/2020] [Revised: 12/30/2020] [Accepted: 01/20/2021] [Indexed: 11/11/2022]
Abstract
PURPOSE To implement a novel, accelerated, 2D radial echo-planar spectroscopic imaging (REPSI) sequence using undersampled radial k-space trajectories and compressed-sensing reconstruction, and to compare results with those from an undersampled Cartesian spectroscopic sequence. METHODS The REPSI sequence was implemented using golden-angle view-ordering on a 3T MRI scanner. Radial and Cartesian echo-planar spectroscopic imaging (EPSI) data were acquired at six acceleration factors, each with time-equivalent scan durations, and reconstructed using compressed sensing with total variation regularization. Results from prospectively and retrospectively undersampled phantom and in vivo brain data were compared over estimated concentrations and Cramer-Rao lower-bound values, normalized RMS errors of reconstructed metabolite maps, and percent absolute differences between fully sampled and reconstructed spectroscopic images. RESULTS The REPSI method with compressed sensing is able to tolerate greater reductions in scan time compared with EPSI. The reconstruction and quantitation metrics (i.e., spectral normalized RMS error maps, metabolite map normalized RMS error values [e.g., for total N-acetyl asparate, REPSI = 9.4% vs EPSI = 16.3%; acceleration factor = 2.5], percent absolute difference maps, and concentration and Cramer-Rao lower-bound estimates) showed that accelerated REPSI can reduce the scan time by a factor of 2.5 while retaining image and quantitation quality. CONCLUSION Accelerated MRSI using undersampled radial echo-planar acquisitions provides greater reconstruction accuracy and more reliable quantitation for a range of acceleration factors compared with time-equivalent compressed-sensing reconstructions of undersampled Cartesian EPSI. Compared to the Cartesian approach, radial undersampling with compressed sensing could help reduce 2D spectroscopic imaging acquisition time, and offers a better trade-off between imaging speed and quality.
Collapse
Affiliation(s)
- Andres Saucedo
- Department of Radiological Sciences, University of California Los Angeles, Los Angeles, California, USA.,Physics and Biology in Medicine Interdepartmental Graduate Program, University of California Los Angeles, Los Angeles, California, USA
| | - Paul M Macey
- School of Nursing, University of California, Los Angeles, Los Angeles, California, USA
| | - M Albert Thomas
- Department of Radiological Sciences, University of California Los Angeles, Los Angeles, California, USA.,Physics and Biology in Medicine Interdepartmental Graduate Program, University of California Los Angeles, Los Angeles, California, USA
| |
Collapse
|
13
|
Maier O, Baete SH, Fyrdahl A, Hammernik K, Harrevelt S, Kasper L, Karakuzu A, Loecher M, Patzig F, Tian Y, Wang K, Gallichan D, Uecker M, Knoll F. CG-SENSE revisited: Results from the first ISMRM reproducibility challenge. Magn Reson Med 2020; 85:1821-1839. [PMID: 33179826 DOI: 10.1002/mrm.28569] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2020] [Revised: 10/02/2020] [Accepted: 10/02/2020] [Indexed: 11/06/2022]
Abstract
PURPOSE The aim of this work is to shed light on the issue of reproducibility in MR image reconstruction in the context of a challenge. Participants had to recreate the results of "Advances in sensitivity encoding with arbitrary k-space trajectories" by Pruessmann et al. METHODS: The task of the challenge was to reconstruct radially acquired multicoil k-space data (brain/heart) following the method in the original paper, reproducing its key figures. Results were compared to consolidated reference implementations created after the challenge, accounting for the two most common programming languages used in the submissions (Matlab/Python). RESULTS Visually, differences between submissions were small. Pixel-wise differences originated from image orientation, assumed field-of-view, or resolution. The reference implementations were in good agreement, both visually and in terms of image similarity metrics. DISCUSSION AND CONCLUSION While the description level of the published algorithm enabled participants to reproduce CG-SENSE in general, details of the implementation varied, for example, density compensation or Tikhonov regularization. Implicit assumptions about the data lead to further differences, emphasizing the importance of sufficient metadata accompanying open datasets. Defining reproducibility quantitatively turned out to be nontrivial for this image reconstruction challenge, in the absence of ground-truth results. Typical similarity measures like NMSE of SSIM were misled by image intensity scaling and outlier pixels. Thus, to facilitate reproducibility, researchers are encouraged to publish code and data alongside the original paper. Future methodological papers on MR image reconstruction might benefit from the consolidated reference implementations of CG-SENSE presented here, as a benchmark for methods comparison.
Collapse
Affiliation(s)
- Oliver Maier
- Institute of Medical Engineering, Graz University of Technology, Graz, Austria
| | - Steven Hubert Baete
- Center for Biomedical Imaging, New York University School of Medicine, New York, NY, USA
| | - Alexander Fyrdahl
- Department of Clinical Physiology, Karolinska University Hospital, and Karolinska Institutet, Stockholm, Sweden
| | - Kerstin Hammernik
- Department of Computing, Imperial College London, London, UK.,Institute of Computer Graphics and Vision, Graz University of Technology, Graz, Austria
| | - Seb Harrevelt
- Department of Biomedical Engineering, Eindhoven University of Technology, Eindhoven, The Netherlands
| | - Lars Kasper
- Institute for Biomedical Engineering, ETH Zurich and University of Zurich, Zurich, Switzerland.,Translational Neuromodeling Unit, Institute for Biomedical Engineering, University of Zurich and ETH Zurich, Zurich, Switzerland.,Techna Institute, University Health Network, Toronto, ON, Canada
| | - Agah Karakuzu
- NeuroPoly Lab, Institute of Biomedical Engineering, Polytechnique Montréal, Montréal, QC, Canada
| | - Michael Loecher
- Department of Radiology, Stanford University, Stanford, CA, USA
| | - Franz Patzig
- Institute for Biomedical Engineering, ETH Zurich and University of Zurich, Zurich, Switzerland
| | - Ye Tian
- Utah Center for Advanced Imaging Research (UCAIR), Department of Radiology and Imaging Sciences, University of Utah, Salt Lake City, UT, USA.,Ming Hsieh Department of Electrical and Computer Engineering, Viterbi School of Engineering, University of Southern California, Los Angeles, CA, USA
| | - Ke Wang
- Department of Electrical Engineering and Computer Sciences, University of California, Berkeley, CA, USA
| | | | - Martin Uecker
- Institute for Diagnostic and Interventional Radiology, University Medical Center Göttingen, Göttingen, Germany.,German Centre for Cardiovascular Research (DZHK), Berlin, Germany.,Cluster of Excellence "Multiscale Bioimaging: from Molecular Machines to Networks of Excitable Cells" (MBExC), University of Göttingen, Göttingen, Germany.,Campus Institute Data Science (CIDAS), University of Göttingen, Göttingen, Germany
| | - Florian Knoll
- Center for Biomedical Imaging, New York University School of Medicine, New York, NY, USA
| |
Collapse
|