1
|
Liu Q, Lin Z, Shen Y, Zhu J, Song J, Zhang C, Lu Y, Xu J. Use of Compressed Sensing Accelerated, Low-Velocity Encoded, Isotropic Resolution, Phase Contrast Magnetic Resonance Angiography for SEEG Electrode Implantation. World Neurosurg 2024; 181:e18-e28. [PMID: 36791880 DOI: 10.1016/j.wneu.2023.02.030] [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: 12/06/2022] [Revised: 02/05/2023] [Accepted: 02/06/2023] [Indexed: 02/15/2023]
Abstract
OBJECTIVE We assessed the feasibility of using compressed sensing accelerated, low-velocity encoded, isotropic resolution phase contrast (CLIP) magnetic resonance angiography (MRA) for avascular trajectory planning of stereoelectroencephalography. METHODS Ten healthy subjects (1 woman and 9 men; age, 33.6 ± 9.0 years) and 20 consecutive patients (12 female patients; age, 22 ± 13.6 years) were enrolled in the present study. The healthy subjects underwent CLIP-MRA, and 3 other phase contrast MRA protocols with conventional parallel imaging (PI) acceleration, including low resolution with twofold PI (PI2), high resolution (HR) with fivefold PI (PI5), and HR-PI2. The patients underwent CLIP-MRA and computed tomography angiography (CTA). The image qualities were evaluated. The numbers and locations of trajectory-vessel conflict detected using CLIP-MRA were noted. RESULTS With similar scan durations, CLIP-MRA achieved higher spatial resolution compared with low resolution with PI2 and detected significantly more branches compared with HR-PI5. With the same spatial resolution, the signal/noise and contrast/noise ratios of CLIP-MRA were higher than those with HR-PI2 with a shorter scan duration. For the 12 adult patients (10 female patients; 28.8 ± 12.7 years), CLIP-MRA had better signal/noise and contrast/noise ratios than CTA. The trajectory had required modification for 14 of the 20 patients (70%), with a proportion of trajectory modification of 10.7% (23 of 215 electrodes). The middle meningeal artery, cortical vessel, and skull vessel were the main vessels with conflict (n = 11, n = 7, and n = 5, respectively). CONCLUSIONS In the present study, CLIP-MRA provided a clear cortical vascular display noninvasively without intravascular contrast and radiation. The middle meningeal artery and diploic and emissary veins were the main conflict vessels and could be clearly displayed using CLIP-MRA but not CTA.
Collapse
Affiliation(s)
- Qiangqiang Liu
- Clinical Neuroscience Center Comprehensive Epilepsy Unit, Department of Neurosurgery, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, People's Republic of China; Clinical Neuroscience Center, Ruijin Hospital Luwan Branch, Shanghai Jiao Tong University School of Medicine, Shanghai, People's Republic of China
| | - Zengping Lin
- United Imaging Healthcare Group, Shanghai, People's Republic of China
| | - Yiwen Shen
- Department of Radiology, Ruijin Hospital Luwan Branch, Shanghai Jiao Tong University School of Medicine, Shanghai, People's Republic of China
| | - Jiachen Zhu
- United Imaging Healthcare Group, Shanghai, People's Republic of China
| | - Jian Song
- Wuhan United Imaging Healthcare Surgical Technology Co., Ltd., Wuhan, People's Republic of China
| | - Chencheng Zhang
- Department of Neurosurgery, Center for Functional Neurosurgery, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, People's Republic of China
| | - Yong Lu
- Clinical Neuroscience Center, Ruijin Hospital Luwan Branch, Shanghai Jiao Tong University School of Medicine, Shanghai, People's Republic of China
| | - Jiwen Xu
- Clinical Neuroscience Center Comprehensive Epilepsy Unit, Department of Neurosurgery, Ruijin Hospital, Shanghai Jiao Tong University School of Medicine, Shanghai, People's Republic of China; Clinical Neuroscience Center, Ruijin Hospital Luwan Branch, Shanghai Jiao Tong University School of Medicine, Shanghai, People's Republic of China.
| |
Collapse
|
2
|
Haldar JP. On Ambiguity in Linear Inverse Problems: Entrywise Bounds on Nearly Data-Consistent Solutions and Entrywise Condition Numbers. IEEE TRANSACTIONS ON SIGNAL PROCESSING : A PUBLICATION OF THE IEEE SIGNAL PROCESSING SOCIETY 2023; 71:1083-1092. [PMID: 37383695 PMCID: PMC10299746 DOI: 10.1109/tsp.2023.3257989] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/30/2023]
Abstract
Ill-posed linear inverse problems appear frequently in various signal processing applications. It can be very useful to have theoretical characterizations that quantify the level of ill-posedness for a given inverse problem and the degree of ambiguity that may exist about its solution. Traditional measures of ill-posedness, such as the condition number of a matrix, provide characterizations that are global in nature. While such characterizations can be powerful, they can also fail to provide full insight into situations where certain entries of the solution vector are more or less ambiguous than others. In this work, we derive novel theoretical lower- and upper-bounds that apply to individual entries of the solution vector, and are valid for all potential solution vectors that are nearly data-consistent. These bounds are agnostic to the noise statistics and the specific method used to solve the inverse problem, and are also shown to be tight. In addition, our results also lead us to introduce an entrywise version of the traditional condition number, which provides a substantially more nuanced characterization of scenarios where certain elements of the solution vector are less sensitive to perturbations than others. Our results are illustrated in an application to magnetic resonance imaging reconstruction, and we include discussions of practical computation methods for large-scale inverse problems, connections between our new theory and the traditional Cramér-Rao bound under statistical modeling assumptions, and potential extensions to cases involving constraints beyond just data-consistency.
Collapse
Affiliation(s)
- Justin P Haldar
- Signal and Image Processing Institute, Ming Hsieh Department of Electrical and Computer Engineering, University of Southern California, Los Angeles, CA, 90089 USA
| |
Collapse
|
3
|
Kim TH, Haldar JP. Efficient iterative solutions to complex-valued nonlinear least-squares problems with mixed linear and antilinear operators. OPTIMIZATION AND ENGINEERING 2022; 23:749-768. [PMID: 35656362 PMCID: PMC9159680 DOI: 10.1007/s11081-021-09604-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/14/2020] [Revised: 02/11/2021] [Accepted: 02/12/2021] [Indexed: 06/15/2023]
Abstract
We consider a setting in which it is desired to find an optimal complex vector x ∈ C N that satisfies A (x) ≈ b in a least-squares sense, where b ∈ C M is a data vector (possibly noise-corrupted), and A (·) : C N → C M is a measurement operator. If A (·) were linear, this reduces to the classical linear least-squares problem, which has a well-known analytic solution as well as powerful iterative solution algorithms. However, instead of linear least-squares, this work considers the more complicated scenario where A (·) is nonlinear, but can be represented as the summation and/or composition of some operators that are linear and some operators that are antilinear. Some common nonlinear operations that have this structure include complex conjugation or taking the real-part or imaginary-part of a complex vector. Previous literature has shown that this kind of mixed linear/antilinear least-squares problem can be mapped into a linear least-squares problem by considering x as a vector in R 2N instead of C N . While this approach is valid, the replacement of the original complex-valued optimization problem with a real-valued optimization problem can be complicated to implement, and can also be associated with increased computational complexity. In this work, we describe theory and computational methods that enable mixed linear/antilinear least-squares problems to be solved iteratively using standard linear least-squares tools, while retaining all of the complex-valued structure of the original inverse problem. An illustration is providedtodemonstratethatthisapproachcansimplifytheimplementationandreduce the computational complexity of iterative solution algorithms.
Collapse
Affiliation(s)
- Tae Hyung Kim
- Department of Electrical and Computer Engineering, University of Southern California, Los Angeles, CA 90089, USA
| | - Justin P. Haldar
- Department of Electrical and Computer Engineering, University of Southern California, Los Angeles, CA 90089, USA
| |
Collapse
|
4
|
Chan CC, Haldar JP. Local perturbation responses and checkerboard tests: Characterization tools for nonlinear MRI methods. Magn Reson Med 2021; 86:1873-1887. [PMID: 34080720 DOI: 10.1002/mrm.28828] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2020] [Revised: 04/09/2021] [Accepted: 04/13/2021] [Indexed: 01/11/2023]
Abstract
PURPOSE Modern methods for MR image reconstruction, denoising, and parameter mapping are becoming increasingly nonlinear, black-box, and at risk of "hallucination." These trends mean that traditional tools for judging confidence in an image (visual quality assessment, point-spread functions (PSFs), g-factor maps, etc.) are less helpful than before. This paper describes and evaluates an approach that can help with assessing confidence in images produced by arbitrary nonlinear methods. THEORY AND METHODS We propose to characterize nonlinear methods by examining the images they produce before and after applying controlled perturbations to the measured data. This results in functions known as local perturbation responses (LPRs) that can provide useful insight into sensitivity, spatial resolution, and aliasing characteristics. LPRs can be viewed as generalizations of classical PSFs, and are are very flexible-they can be applied to arbitary nonlinear methods and arbitrary datasets across a range of different reconstruction, denoising, and parameter mapping applications. Importantly, LPRs do not require a ground truth image. RESULTS Impulse-based and checkerboard-pattern LPRs are demonstrated in image reconstruction and denoising scenarios. We observe that these LPRs provide insights into spatial resolution, signal leakage, and aliasing that are not available with other methods. We also observe that popular reference-based image quality metrics (eg, mean-squared error and structural similarity) do not always correlate with good LPR characteristics. CONCLUSIONS LPRs are a useful tool that can be used to characterize and assess confidence in nonlinear MR methods, and provide insights that are distinct from and complementary to existing quality assessments.
Collapse
Affiliation(s)
- Chin-Cheng Chan
- Electrical and Computer Engineering, University of Southern California, Los Angeles, California, USA.,Signal and Image Processing Institute, University of Southern California, Los Angeles, California, USA
| | - Justin P Haldar
- Electrical and Computer Engineering, University of Southern California, Los Angeles, California, USA.,Signal and Image Processing Institute, University of Southern California, Los Angeles, California, USA
| |
Collapse
|
5
|
Lee Y, Wilm BJ, Brunner DO, Gross S, Schmid T, Nagy Z, Pruessmann KP. On the signal-to-noise ratio benefit of spiral acquisition in diffusion MRI. Magn Reson Med 2020; 85:1924-1937. [PMID: 33280160 DOI: 10.1002/mrm.28554] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2020] [Revised: 09/15/2020] [Accepted: 09/22/2020] [Indexed: 12/20/2022]
Abstract
PURPOSE Spiral readouts combine several favorable properties that promise superior net sensitivity for diffusion imaging. The purpose of this study is to verify the signal-to-noise ratio (SNR) benefit of spiral acquisition in comparison with current echo-planar imaging (EPI) schemes. METHODS Diffusion-weighted in vivo brain data from three subjects were acquired with a single-shot spiral sequence and several variants of single-shot EPI, including full-Fourier and partial-Fourier readouts as well as different diffusion-encoding schemes. Image reconstruction was based on an expanded signal model including field dynamics obtained by concurrent field monitoring. The effective resolution of each sequence was matched to that of full-Fourier EPI with 1 mm nominal resolution. SNR maps were generated by determining the noise statistics of the raw data and analyzing the propagation of equivalent synthetic noise through image reconstruction. Using the same approach, maps of noise amplification due to parallel imaging (g-factor) were calculated for different acceleration factors. RESULTS Relative to full-Fourier EPI at b = 0 s/mm2 , spiral acquisition yielded SNR gains of 42-88% and 40-89% in white and gray matter, respectively, depending on the diffusion-encoding scheme. Relative to partial-Fourier EPI, the gains were 36-44% and 34-42%. Spiral g-factor maps exhibited less spatial variation and lower maxima than their EPI counterparts. CONCLUSION Spiral readouts achieve significant SNR gains in the order of 40-80% over EPI in diffusion imaging at 3T. Combining systematic effects of shorter echo time, readout efficiency, and favorable g-factor behavior, similar benefits are expected across clinical and neurosciences uses of diffusion imaging.
Collapse
Affiliation(s)
- Yoojin Lee
- Institute for Biomedical Engineering, ETH Zurich and University of Zurich, Zurich, Switzerland.,Laboratory for Social and Neural Systems Research, University of Zurich, Zurich, Switzerland
| | - Bertram J Wilm
- Institute for Biomedical Engineering, ETH Zurich and University of Zurich, Zurich, Switzerland
| | - David O Brunner
- Institute for Biomedical Engineering, ETH Zurich and University of Zurich, Zurich, Switzerland
| | - Simon Gross
- Institute for Biomedical Engineering, ETH Zurich and University of Zurich, Zurich, Switzerland
| | - Thomas Schmid
- Institute for Biomedical Engineering, ETH Zurich and University of Zurich, Zurich, Switzerland
| | - Zoltan Nagy
- Laboratory for Social and Neural Systems Research, University of Zurich, Zurich, Switzerland
| | - Klaas P Pruessmann
- Institute for Biomedical Engineering, ETH Zurich and University of Zurich, Zurich, Switzerland
| |
Collapse
|
6
|
Kettinger AO, Setsompop K, Kannengiesser SAR, Breuer FA, Vidnyanszky Z, Blaimer M. Full utilization of conjugate symmetry: combining virtual conjugate coil reconstruction with partial Fourier imaging for g-factor reduction in accelerated MRI. Magn Reson Med 2019; 82:1073-1090. [PMID: 31081561 DOI: 10.1002/mrm.27799] [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/2018] [Revised: 02/27/2019] [Accepted: 04/13/2019] [Indexed: 11/10/2022]
Abstract
PURPOSE In this study we propose a method to combine the parallel virtual conjugate coil (VCC) reconstruction with partial Fourier (PF) acquisition to improve reconstruction conditioning and reduce noise amplification in accelerated MRI where PF is used. METHODS Accelerated measurements are reconstructed in k-space by GRAPPA, with a VCC reconstruction kernel trained and applied in the central, symmetrically sampled part of k-space, while standard reconstruction is performed on the asymmetrically sampled periphery. The two reconstructed regions are merged to form a full reconstructed dataset, followed by PF reconstruction. The method is tested in vivo using T1-weighted spin-echo and T2*-weighted gradient-echo echo planar imaging (EPI) sequences, using both in-plane and simultaneous multislice (SMS) acceleration, at 1.5T and 3T field strengths. Noise amplification is estimated with theoretical calculations and pseudo-multiple-replica computations, for different PF factors, using zero-filling, homodyne, and projection onto convex sets (POCS) PF reconstruction. RESULTS Depending on the PF algorithm and the inherent benefit of VCC reconstruction without PF, approximately 35% to 80%, 15% to 60%, and 5% to 30% of that intrinsic SNR gain can be retained for PF factors 7/8, 6/8, and 5/8, respectively, by including the VCC signals in the reconstruction. Compared with VCC-reconstructed acquisitions of higher acceleration, without PF, but having the same net acceleration, the combined method can provide a higher SNR if the inherent benefit of VCC is low or moderate. CONCLUSION The proposed technique enables the partial application of VCC reconstruction to measurements with PF using either in-plane or SMS acceleration, and therefore can reduce the noise amplification of such acquisitions.
Collapse
Affiliation(s)
- Adam O Kettinger
- Brain Imaging Centre, Research Centre for Natural Sciences, Hungarian Academy of Sciences, Budapest, Hungary.,Department of Nuclear Techniques, Budapest University of Technology and Economics, Budapest, Hungary.,Siemens Healthcare GmbH, Erlangen, Germany
| | - Kawin Setsompop
- Athinoula A. Martinos Center for Biomedical Imaging, Charlestown, Massachusetts.,Department of Radiology, Harvard Medical School, Boston, Massachusetts.,Harvard-MIT Health Sciences and Technology, MIT, Cambridge, Massachusetts
| | | | - Felix A Breuer
- Magnetic Resonance and X-ray Imaging Department, Fraunhofer Development Center X-ray Technology (EZRT), Würzburg, Germany
| | - Zoltan Vidnyanszky
- Brain Imaging Centre, Research Centre for Natural Sciences, Hungarian Academy of Sciences, Budapest, Hungary
| | - Martin Blaimer
- Magnetic Resonance and X-ray Imaging Department, Fraunhofer Development Center X-ray Technology (EZRT), Würzburg, Germany
| |
Collapse
|
7
|
Walheim J, Gotschy A, Kozerke S. On the limitations of partial Fourier acquisition in phase-contrast MRI of turbulent kinetic energy. Magn Reson Med 2018; 81:514-523. [PMID: 30265753 DOI: 10.1002/mrm.27397] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2017] [Revised: 05/04/2018] [Accepted: 05/20/2018] [Indexed: 11/11/2022]
Abstract
PURPOSE To investigate limitations of partial Fourier acquisition in phase-contrast MRI of turbulent kinetic energy (TKE). METHODS To assess the validity of partial Fourier reconstruction of TKE and phase images, computational fluid dynamics data of mean and turbulent velocities in a stenotic U-bend phantom was used. Partial Fourier acquisition with 75% k-space coverage was simulated and TKE data were reconstructed using zero-filling, homodyne reconstruction, and the method of projections onto convex sets (POCS). Results were compared to data from fully sampled k-space and 75% symmetric sampling. In addition, compressed sensing (CS) reconstruction was compared for a standard variable density sampling pattern and a variable density sampling pattern combined with 75% partial Fourier. For illustration purposes, in vivo examples of velocity magnitude and TKE maps of aortic flow reconstructed with the different methods are provided. RESULTS In accordance with theory, partial Fourier reconstruction of TKE maps from phase-contrast data results in artifacts relative to fully sampled data. It is demonstrated that neither homodyne reconstruction nor POCS can improve reconstruction of TKE data with respect to zero-filling reconstruction when compared to ground-truth (RMS error: 4.70%, 4.34%, and 2.45% for homodyne, POCS, and zero-filling reconstruction of in vivo data, respectively). CS reconstruction from data acquired with partial Fourier did not recover the resolution loss incurred by partial Fourier sampling. CONCLUSION Partial Fourier reconstruction of TKE maps from phase-contrast data does not yield a benefit over zero-filling reconstruction. In consequence, symmetric sampling is preferred over partial Fourier acquisition for a given number of phase-encodes in phase-contrast MRI.
Collapse
Affiliation(s)
- Jonas Walheim
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | - Alexander Gotschy
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland.,Department of Cardiology, University Hospital Zurich, Zurich, Switzerland
| | - Sebastian Kozerke
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| |
Collapse
|
8
|
Aiello M, Cavaliere C, Marchitelli R, d'Albore A, De Vita E, Salvatore M. Hybrid PET/MRI Methodology. INTERNATIONAL REVIEW OF NEUROBIOLOGY 2018; 141:97-128. [PMID: 30314608 DOI: 10.1016/bs.irn.2018.07.026] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/04/2022]
Abstract
The hybrid PET/MR scanner represents the first implementation of the effective integration of two modalities allowing truly synchronous/simultaneous acquisition of their imaging signals. This integration, resulting from the innovation and development of specific hardware components has paved the way for new approaches in the study of neurodegenerative diseases. This chapter will describe the hardware development that has led to the availability of different clinical solutions for PET/MR imaging as well as the still-open technological challenges and opportunities related to the processing and exploitation of the simultaneous acquisition in neurological studies.
Collapse
Affiliation(s)
| | | | | | | | - Enrico De Vita
- Department of Biomedical Engineering, School of Biomedical Engineering & Imaging Sciences, King's College London, King's Health Partners, St Thomas' Hospital, London, United Kingdom
| | | |
Collapse
|
9
|
Ong F, Cheng J, Lustig M. General phase regularized reconstruction using phase cycling. Magn Reson Med 2018; 80:112-125. [PMID: 29159989 PMCID: PMC5876131 DOI: 10.1002/mrm.27011] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2017] [Revised: 10/23/2017] [Accepted: 10/24/2017] [Indexed: 12/21/2022]
Abstract
PURPOSE To develop a general phase regularized image reconstruction method, with applications to partial Fourier imaging, water-fat imaging and flow imaging. THEORY AND METHODS The problem of enforcing phase constraints in reconstruction was studied under a regularized inverse problem framework. A general phase regularized reconstruction algorithm was proposed to enable various joint reconstruction of partial Fourier imaging, water-fat imaging and flow imaging, along with parallel imaging (PI) and compressed sensing (CS). Since phase regularized reconstruction is inherently non-convex and sensitive to phase wraps in the initial solution, a reconstruction technique, named phase cycling, was proposed to render the overall algorithm invariant to phase wraps. The proposed method was applied to retrospectively under-sampled in vivo datasets and compared with state of the art reconstruction methods. RESULTS Phase cycling reconstructions showed reduction of artifacts compared to reconstructions without phase cycling and achieved similar performances as state of the art results in partial Fourier, water-fat and divergence-free regularized flow reconstruction. Joint reconstruction of partial Fourier + water-fat imaging + PI + CS, and partial Fourier + divergence-free regularized flow imaging + PI + CS were demonstrated. CONCLUSION The proposed phase cycling reconstruction provides an alternative way to perform phase regularized reconstruction, without the need to perform phase unwrapping. It is robust to the choice of initial solutions and encourages the joint reconstruction of phase imaging applications. Magn Reson Med 80:112-125, 2018. © 2017 International Society for Magnetic Resonance in Medicine.
Collapse
Affiliation(s)
- Frank Ong
- Department of Electrical Engineering and Computer Sciences, University of California, Berkeley, California
| | - Joseph Cheng
- Department of Radiology, Stanford University, Stanford, California
| | - Michael Lustig
- Department of Electrical Engineering and Computer Sciences, University of California, Berkeley, California
| |
Collapse
|
10
|
Derakhshan JJ, McDonald ES, Siegelman ES, Schnall MD, Wehrli FW. Characterizing and eliminating errors in enhancement and subtraction artifacts in dynamic contrast-enhanced breast MRI: Chemical shift artifact of the third kind. Magn Reson Med 2018; 79:2277-2289. [PMID: 28840613 PMCID: PMC5811365 DOI: 10.1002/mrm.26879] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2017] [Revised: 07/27/2017] [Accepted: 07/30/2017] [Indexed: 12/29/2022]
Abstract
PURPOSE To characterize errors in enhancement in breast dynamic contrast-enhanced (DCE) MRI studies as a function of echo time and determine the source of dark band artifacts in clinical subtraction images. METHODS Computer simulations, oil and water substitute (methylene chloride), as well as an American College of Radiology quality control phantom were tested. Routine clinical DCE breast MRI study was bracketed with (accelerated) in-phase DCE acquisitions in five patients. RESULTS Simulation results demonstrated up to -160% suppression of the expected enhancement caused by differential enhancement of fat and water. Two-dimensional gradient-recalled echo and fat-suppressed 3D GRE phantom imaging confirmed the simulation results and showed that fat suppression does not eliminate the artifact. In vivo in-phase DCE images showed increased enhancement consistent with predictions and also confirmed increased spatial blurring on in-phase 3D gradient-recalled echo images. Combined multi-dimensional partial Fourier and parallel imaging provided a time-equivalent in-phase DCE MRI acquisition. CONCLUSION Errors in expected enhancement occur in DCE breast MRI subtraction images because of differential enhancement of fat and water and incomplete fat signal suppression. These errors can lead to artificial suppression of enhancement as well as dark band artifacts on subtraction images. These artifacts can be eliminated with a time-equivalent in-phase fat-suppressed 3D gradient-recalled echo sequence. Understanding chemical shift artifact of the third kind, a unique artifact of artificial enhancement suppression in the presence of intravoxel fat and water signal, will aid DCE breast MRI image interpretation. In-phase acquisitions (combined with simultaneous minimum echo time or opposed-phase echoes) may facilitate qualitative, quantitative and longitudinal analysis of contrast enhancement. Magn Reson Med 79:2277-2289, 2018. © 2017 International Society for Magnetic Resonance in Medicine.
Collapse
Affiliation(s)
- Jamal J Derakhshan
- Department of Radiology, University of Pennsylvania, Philadelphia, Pennsylvania, USA
| | - Elizabeth S McDonald
- Department of Radiology, University of Pennsylvania, Philadelphia, Pennsylvania, USA
| | - Evan S Siegelman
- Department of Radiology, University of Pennsylvania, Philadelphia, Pennsylvania, USA
| | - Mitchell D Schnall
- Department of Radiology, University of Pennsylvania, Philadelphia, Pennsylvania, USA
| | - Felix W Wehrli
- Department of Radiology, University of Pennsylvania, Philadelphia, Pennsylvania, USA
| |
Collapse
|
11
|
Olafsson VT, Noll DC, Fessler JA. Fast Spatial Resolution Analysis of Quadratic Penalized Least-Squares Image Reconstruction With Separate Real and Imaginary Roughness Penalty: Application to fMRI. IEEE TRANSACTIONS ON MEDICAL IMAGING 2018; 37:604-614. [PMID: 29408788 PMCID: PMC5804832 DOI: 10.1109/tmi.2017.2768825] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Penalized least-squares iterative image reconstruction algorithms used for spatial resolution-limited imaging, such as functional magnetic resonance imaging (fMRI), commonly use a quadratic roughness penalty to regularize the reconstructed images. When used for complex-valued images, the conventional roughness penalty regularizes the real and imaginary parts equally. However, these imaging methods sometimes benefit from separate penalties for each part. The spatial smoothness from the roughness penalty on the reconstructed image is dictated by the regularization parameter(s). One method to set the parameter to a desired smoothness level is to evaluate the full width at half maximum of the reconstruction method's local impulse response. Previous work has shown that when using the conventional quadratic roughness penalty, one can approximate the local impulse response using an FFT-based calculation. However, that acceleration method cannot be applied directly for separate real and imaginary regularization. This paper proposes a fast and stable calculation for this case that also uses FFT-based calculations to approximate the local impulse responses of the real and imaginary parts. This approach is demonstrated with a quadratic image reconstruction of fMRI data that uses separate roughness penalties for the real and imaginary parts.
Collapse
|
12
|
Ahmad R, Hu HH, Krishnamurthy R, Krishnamurthy R. Reducing sedation for pediatric body MRI using accelerated and abbreviated imaging protocols. Pediatr Radiol 2018; 48:37-49. [PMID: 29292482 DOI: 10.1007/s00247-017-3987-6] [Citation(s) in RCA: 51] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/03/2017] [Revised: 06/13/2017] [Accepted: 09/12/2017] [Indexed: 12/01/2022]
Abstract
Magnetic resonance imaging (MRI) is an established diagnostic imaging tool for investigating pediatric disease. MRI allows assessment of structure, function, and morphology in cardiovascular imaging, as well as tissue characterization in body imaging, without the use of ionizing radiation. For MRI in children, sedation and general anesthesia (GA) are often utilized to suppress patient motion, which can otherwise compromise image quality and diagnostic efficacy. However, evidence is emerging that use of sedation and GA in children might have long-term neurocognitive side effects, in addition to the short-term procedure-related risks. These concerns make risk-benefit assessment of sedation and GA more challenging. Therefore, reducing or eliminating the need for sedation and GA is an important goal of imaging innovation and research in pediatric MRI. In this review, the authors focus on technical and clinical approaches to reducing and eliminating the use of sedation in the pediatric population based on image acquisition acceleration and imaging protocols abbreviation. This paper covers important physiological and technical considerations for pediatric body MR imaging and discusses MRI techniques that offer the potential of recovering diagnostic-quality images from accelerated scans. In this review, the authors also introduce the concept of reporting elements for important indications for pediatric body MRI and use this as a basis for abbreviating the MR protocols. By employing appropriate accelerated and abbreviated approaches based on an understanding of the imaging needs and reporting elements for a given clinical indication, it is possible to reduce sedation and GA for pediatric chest, cardiovascular and abdominal MRI.
Collapse
Affiliation(s)
- Rizwan Ahmad
- Department of Biomedical Engineering, The Ohio State University, Columbus, OH, USA
| | - Houchun Harry Hu
- Department of Radiology, Nationwide Children's Hospital, 700 Children's Drive, Columbus, OH, 43205, USA
| | - Ramkumar Krishnamurthy
- Department of Radiology, Nationwide Children's Hospital, 700 Children's Drive, Columbus, OH, 43205, USA
| | - Rajesh Krishnamurthy
- Department of Radiology, Nationwide Children's Hospital, 700 Children's Drive, Columbus, OH, 43205, USA.
| |
Collapse
|
13
|
Kettinger AO, Kannengiesser SAR, Breuer FA, Vidnyanszky Z, Blaimer M. Controlling the object phase for g‐factor reduction in phase‐Constrained parallel MRI using spatially selective RF pulses. Magn Reson Med 2017; 79:2113-2125. [DOI: 10.1002/mrm.26890] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2017] [Revised: 08/05/2017] [Accepted: 08/07/2017] [Indexed: 12/13/2022]
Affiliation(s)
- Adam O. Kettinger
- Department of Nuclear TechniquesBudapest University of Technology and EconomicsBudapest Hungary
- Brain Imaging Centre, Research Centre for Natural Sciences, Hungarian Academy of SciencesBudapest Hungary
| | | | - Felix A. Breuer
- Magnetic Resonance and X‐ray Imaging DepartmentFraunhofer Development Center X‐ray Technology (EZRT)Würzburg Germany
| | - Zoltan Vidnyanszky
- Brain Imaging Centre, Research Centre for Natural Sciences, Hungarian Academy of SciencesBudapest Hungary
| | - Martin Blaimer
- Magnetic Resonance and X‐ray Imaging DepartmentFraunhofer Development Center X‐ray Technology (EZRT)Würzburg Germany
| |
Collapse
|
14
|
Hamilton J, Franson D, Seiberlich N. Recent advances in parallel imaging for MRI. PROGRESS IN NUCLEAR MAGNETIC RESONANCE SPECTROSCOPY 2017; 101:71-95. [PMID: 28844222 PMCID: PMC5927614 DOI: 10.1016/j.pnmrs.2017.04.002] [Citation(s) in RCA: 106] [Impact Index Per Article: 15.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/30/2017] [Revised: 04/09/2017] [Accepted: 04/17/2017] [Indexed: 05/22/2023]
Abstract
Magnetic Resonance Imaging (MRI) is an essential technology in modern medicine. However, one of its main drawbacks is the long scan time needed to localize the MR signal in space to generate an image. This review article summarizes some basic principles and recent developments in parallel imaging, a class of image reconstruction techniques for shortening scan time. First, the fundamentals of MRI data acquisition are covered, including the concepts of k-space, undersampling, and aliasing. It is demonstrated that scan time can be reduced by sampling a smaller number of phase encoding lines in k-space; however, without further processing, the resulting images will be degraded by aliasing artifacts. Nearly all modern clinical scanners acquire data from multiple independent receiver coil arrays. Parallel imaging methods exploit properties of these coil arrays to separate aliased pixels in the image domain or to estimate missing k-space data using knowledge of nearby acquired k-space points. Three parallel imaging methods-SENSE, GRAPPA, and SPIRiT-are described in detail, since they are employed clinically and form the foundation for more advanced methods. These techniques can be extended to non-Cartesian sampling patterns, where the collected k-space points do not fall on a rectangular grid. Non-Cartesian acquisitions have several beneficial properties, the most important being the appearance of incoherent aliasing artifacts. Recent advances in simultaneous multi-slice imaging are presented next, which use parallel imaging to disentangle images of several slices that have been acquired at once. Parallel imaging can also be employed to accelerate 3D MRI, in which a contiguous volume is scanned rather than sequential slices. Another class of phase-constrained parallel imaging methods takes advantage of both image magnitude and phase to achieve better reconstruction performance. Finally, some applications are presented of parallel imaging being used to accelerate MR Spectroscopic Imaging.
Collapse
Affiliation(s)
- Jesse Hamilton
- Biomedical Engineering, Case Western Reserve University, Cleveland, OH, USA.
| | - Dominique Franson
- Biomedical Engineering, Case Western Reserve University, Cleveland, OH, USA.
| | - Nicole Seiberlich
- Biomedical Engineering, Case Western Reserve University, Cleveland, OH, USA; Radiology, University Hospitals Cleveland Medical Center, Cleveland, OH, USA.
| |
Collapse
|
15
|
Bydder M, Rapacchi S, Girard O, Guye M, Ranjeva JP. Trimmed autocalibrating k-space estimation based on structured matrix completion. Magn Reson Imaging 2017; 43:88-94. [PMID: 28716683 DOI: 10.1016/j.mri.2017.07.015] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2017] [Revised: 07/07/2017] [Accepted: 07/13/2017] [Indexed: 11/28/2022]
Abstract
PURPOSE Parallel imaging allows the reconstruction of undersampled data from multiple coils. This provides a means to reject and regenerate corrupt data (e.g. from motion artefact). The purpose of this work is to approach this problem using the SAKE parallel imaging method. THEORY AND METHODS Parallel imaging methods typically require calibration by fully sampling the center of k-space. This is a challenge in the presence of corrupted data, since the calibration data may be corrupted which leads to an errors-in-variables problem that cannot be solved by least squares or even iteratively reweighted least squares. The SAKE method, based on matrix completion and structured low rank approximation, was modified to detect and trim these errors from the data. RESULTS Simulated and actual corrupted datasets were reconstructed with SAKE, the proposed approach and a more standard reconstruction method (based on solving a linear equation) with a data rejection criterion. The proposed approach was found to reduce artefacts considerably in comparison to the other two methods. CONCLUSION SAKE with data trimming improves on previous methods for reconstructing images from grossly corrupted data.
Collapse
Affiliation(s)
- Mark Bydder
- Aix-Marseille Université, Centre de Résonance Magnétique Biologique et Médicale, UMR CNRS 7339, Marseille, France.
| | - Stanislas Rapacchi
- Aix-Marseille Université, Centre de Résonance Magnétique Biologique et Médicale, UMR CNRS 7339, Marseille, France
| | - Olivier Girard
- Aix-Marseille Université, Centre de Résonance Magnétique Biologique et Médicale, UMR CNRS 7339, Marseille, France
| | - Maxime Guye
- Aix-Marseille Université, Centre de Résonance Magnétique Biologique et Médicale, UMR CNRS 7339, Marseille, France; Assistance Publique - Hôpitaux de Marseille, CEMREM, Pôle d'Imagerie Médicale, CHU Timone, Marseille, France
| | - Jean-Philippe Ranjeva
- Aix-Marseille Université, Centre de Résonance Magnétique Biologique et Médicale, UMR CNRS 7339, Marseille, France
| |
Collapse
|
16
|
Ramos-Llorden G, den Dekker AJ, Sijbers J. Partial Discreteness: A Novel Prior for Magnetic Resonance Image Reconstruction. IEEE TRANSACTIONS ON MEDICAL IMAGING 2017; 36:1041-1053. [PMID: 28026759 DOI: 10.1109/tmi.2016.2645122] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2023]
Abstract
An important factor influencing the quality of magnetic resonance (MR) images is the reconstruction method that is employed, and specifically, the type of prior knowledge that is exploited during reconstruction. In this work, we introduce a new type of prior knowledge, partial discreteness (PD), where a small number of regions in the image are assumed to be homogeneous and can be well represented by a constant magnitude. In particular, we mathematically formalize the partial discreteness property based on a Gaussian Mixture Model (GMM) and derive a partial discreteness image representation that characterizes the salient features of partially discrete images: a constant intensity in homogeneous areas and texture in heterogeneous areas. The partial discreteness representation is then used to construct a novel prior dedicated to the reconstruction of partially discrete MR images. The strength of the proposed prior is demonstrated on various simulated and real k-space data-based experiments with partially discrete images. Results demonstrate that the PD algorithm performs competitively with state-of-the-art reconstruction methods, being flexible and easy to implement.
Collapse
|
17
|
Abstract
We proposed and assessed a modified asymmetric Fourier imaging (AFI) technique named real-based polarity-preserving AFI (RepAFI), in which the low-pass filter kernel for background phase estimation in AFI is optimized to preserve the magnetization polarity information for blood vessels and cerebrospinal fluid (CSF) even for data obtained using phase-sensitive inversion-recovery spin-echo-based (PSIR-SE) sequences with asymmetrical sampling in the k-space. Our proposed RepAFI technique achieves a practical balance of image quality and simplicity to provide better performance than conventional AFI methods.
Collapse
Affiliation(s)
- Tokunori Kimura
- Clinical Application Research and Development Department, Center for Medical Research and Development, Toshiba Medical Systems Corporation
| | | |
Collapse
|
18
|
Sparse Parallel MRI Based on Accelerated Operator Splitting Schemes. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2016; 2016:1724630. [PMID: 27746824 PMCID: PMC5056009 DOI: 10.1155/2016/1724630] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/26/2016] [Accepted: 08/29/2016] [Indexed: 11/17/2022]
Abstract
Recently, the sparsity which is implicit in MR images has been successfully exploited for fast MR imaging with incomplete acquisitions. In this paper, two novel algorithms are proposed to solve the sparse parallel MR imaging problem, which consists of l1 regularization and fidelity terms. The two algorithms combine forward-backward operator splitting and Barzilai-Borwein schemes. Theoretically, the presented algorithms overcome the nondifferentiable property in l1 regularization term. Meanwhile, they are able to treat a general matrix operator that may not be diagonalized by fast Fourier transform and to ensure that a well-conditioned optimization system of equations is simply solved. In addition, we build connections between the proposed algorithms and the state-of-the-art existing methods and prove their convergence with a constant stepsize in Appendix. Numerical results and comparisons with the advanced methods demonstrate the efficiency of proposed algorithms.
Collapse
|
19
|
Takatsu Y, Ueyama T, Miyati T, Yamamura K. Simulation of the modulation transfer function dependent on the partial Fourier fraction in dynamic contrast enhancement magnetic resonance imaging. AUSTRALASIAN PHYSICAL & ENGINEERING SCIENCES IN MEDICINE 2016; 39:825-831. [DOI: 10.1007/s13246-016-0474-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/28/2016] [Accepted: 08/15/2016] [Indexed: 02/03/2023]
|
20
|
Kim TH, Setsompop K, Haldar JP. LORAKS makes better SENSE: Phase-constrained partial fourier SENSE reconstruction without phase calibration. Magn Reson Med 2016; 77:1021-1035. [PMID: 27037836 DOI: 10.1002/mrm.26182] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2015] [Revised: 02/03/2016] [Accepted: 02/04/2016] [Indexed: 12/13/2022]
Abstract
PURPOSE Parallel imaging and partial Fourier acquisition are two classical approaches for accelerated MRI. Methods that combine these approaches often rely on prior knowledge of the image phase, but the need to obtain this prior information can place practical restrictions on the data acquisition strategy. In this work, we propose and evaluate SENSE-LORAKS, which enables combined parallel imaging and partial Fourier reconstruction without requiring prior phase information. THEORY AND METHODS The proposed formulation is based on combining the classical SENSE model for parallel imaging data with the more recent LORAKS framework for MR image reconstruction using low-rank matrix modeling. Previous LORAKS-based methods have successfully enabled calibrationless partial Fourier parallel MRI reconstruction, but have been most successful with nonuniform sampling strategies that may be hard to implement for certain applications. By combining LORAKS with SENSE, we enable highly accelerated partial Fourier MRI reconstruction for a broader range of sampling trajectories, including widely used calibrationless uniformly undersampled trajectories. RESULTS Our empirical results with retrospectively undersampled datasets indicate that when SENSE-LORAKS reconstruction is combined with an appropriate k-space sampling trajectory, it can provide substantially better image quality at high-acceleration rates relative to existing state-of-the-art reconstruction approaches. CONCLUSION The SENSE-LORAKS framework provides promising new opportunities for highly accelerated MRI. Magn Reson Med 77:1021-1035, 2017. © 2016 International Society for Magnetic Resonance in Medicine.
Collapse
Affiliation(s)
- Tae Hyung Kim
- Department of Electrical Engineering, University of Southern California, Los Angeles, California, USA
| | - Kawin Setsompop
- Department of Radiology, Harvard Medical School, Boston, Massachusetts, USA
| | - Justin P Haldar
- Department of Electrical Engineering, University of Southern California, Los Angeles, California, USA
| |
Collapse
|
21
|
Brodsky EK, Bultman EM, Johnson KM, Horng DE, Schelman WR, Block WF, Reeder SB. High-spatial and high-temporal resolution dynamic contrast-enhanced perfusion imaging of the liver with time-resolved three-dimensional radial MRI. Magn Reson Med 2015; 71:934-41. [PMID: 23519837 DOI: 10.1002/mrm.24727] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
PURPOSE Detection, characterization, and monitoring the treatment of hepatocellular carcinomas (HCC) in patients with cirrhosis is challenging because of their variable and rapid arterial enhancement. Multiphase dynamic contrast-enhanced MRI is used clinically for HCC assessment; however, the method suffers from limited temporal resolution and difficulty in coordinating imaging and breath-hold timing within a narrow temporal window of interest. In this article, a volumetric, high-spatial resolution, and high-temporal resolution dynamic contrast-enhanced liver imaging method for improved detection and characterization of HCC is demonstrated. METHODS A time-resolved three-dimensional radial acquisition with iterative sensitivity-encoding reconstruction images the entire abdomen and thorax with high spatial and temporal resolution, using real-time three-dimensional fluoroscopy to match the breath hold to contrast arrival. The sequence was tested on 17 subjects, including eight patients with HCC or other hypervascular focal lesions. RESULTS This technique was successful in acquiring volumetric imaging of the entire liver with 2.1-mm isotropic spatial and true 4-s temporal resolution. CONCLUSION This technique may be suitable for detecting, characterizing, and monitoring the treatment of HCC. It also holds significant potential for perfusion modeling, which may provide a noninvasive means to rapidly determine the efficacy of chemotherapeutic agents in these tumors over the entire liver volume.
Collapse
Affiliation(s)
- Ethan K Brodsky
- Department of Radiology, University of Wisconsin, Madison, Wisconsin, USA; Department of Medical Physics, University of Wisconsin, Madison, Wisconsin, USA; Department of Biomedical Engineering, University of Wisconsin, Madison, Wisconsin, USA
| | | | | | | | | | | | | |
Collapse
|
22
|
Kimura T, Shigeta T. Magnitude-based Asymmetric Fourier Imaging (MagAFI). Magn Reson Med Sci 2015; 15:94-104. [PMID: 26346401 DOI: 10.2463/mrms.2014-0147] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022] Open
Abstract
PURPOSE We propose and assess 2 novel asymmetric Fourier imaging (AFI) techniques, magnitude-based AFI (MagAFI) and MagAFI combined with projection on to convex sets (POCS) (MagAFI+POCS). MagAFI does not require phase information because it uses only the magnitude image with zero-filling. MagAFI+POCS requires phase information but further reduces image errors. MATERIALS AND METHODS We initially compared phase maps obtained using asymmetrically sampled data for the whole of the k-space and symmetrically sampled data for the low frequency part of the k-space. We used one-dimensional simulation data and 3-dimensional gradient echo data for 2 different echo times (TEs) of the brains of volunteers and assessed the differences between the image reconstructed from the full k-space data and AFI images reconstructed from truncated k-space data. We generated AFI images in this study using the zero-filling, Margosian (homodyne), Margosian+POCS (standard POCS), MagAFI, and MagAFI+POCS techniques. RESULTS We confirmed the assumption of smaller phase errors for the full k-space data than for the symmetric low frequency k-space data. Our proposed MagAFI technique provides images with smaller phase-induced errors than those obtained using conventional methods, including standard POCS methods, which have been regarded as the best methods. MagAFI+POCS improves image quality as well as robustness. CONCLUSION Our proposed MagAFI technique achieves a practical balance of image quality and simplicity to perform better than conventional methods using only the 0-filled magnitude image. Combined with POCS, this technique can produce images of even better quality.
Collapse
Affiliation(s)
- Tokunori Kimura
- Clinical Application Research and Development Department, Center for Medical Research and Development, Toshiba Medical Systems Corporation
| | | |
Collapse
|
23
|
Narsude M, Gallichan D, van der Zwaag W, Gruetter R, Marques JP. Three-dimensional echo planar imaging with controlled aliasing: A sequence for high temporal resolution functional MRI. Magn Reson Med 2015; 75:2350-61. [PMID: 26173572 DOI: 10.1002/mrm.25835] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2014] [Revised: 06/16/2015] [Accepted: 06/16/2015] [Indexed: 11/06/2022]
Abstract
PURPOSE In this work, we combine three-dimensional echo planar imaging (3D-EPI) with controlled aliasing to substantially increase temporal resolution in whole-brain functional MRI (fMRI) while minimizing geometry-factor (g-factor) losses. THEORY AND METHODS The study was performed on a 7 Tesla scanner equipped with a 32-channel receive coil. The proposed 3D-EPI-CAIPI sequence was evaluated for: (i) image quality, compared with a conventionally undersampled parallel imaging acquisition; (ii) temporal resolution, the ability to sample and remove physiological signal fluctuations from the fMRI signal of interest and (iii) the ability to distinguish small changes in hemodynamic responses in an event-related fMRI experiment. RESULTS Whole-brain fMRI data with a voxel size of 2 × 2 × 2 mm(3) and temporal resolution of 371 ms could be acquired with acceptable image quality. Ten-fold parallel imaging accelerated 3D-EPI-CAIPI data were shown to lower the maximum g-factor losses up to 62% with respect to a 10-fold accelerated 3D-EPI dataset. FMRI with 400 ms temporal resolution allowed the detection of time-to-peak variations in functional responses due to multisensory facilitation in temporal, occipital and frontal cortices. CONCLUSION 3D-EPI-CAIPI allows increased temporal resolution that enables better characterization of physiological noise, thus improving sensitivity to signal changes of interest. Furthermore, subtle changes in hemodynamic response dynamics can be studied in shorter scan times by avoiding the need for jittering. Magn Reson Med 75:2350-2361, 2016. © 2015 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Mayur Narsude
- Laboratory for Functional and Metabolic Imaging, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland.,Department of Radiology, University of Lausanne, Lausanne, Switzerland
| | - Daniel Gallichan
- Laboratory for Functional and Metabolic Imaging, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Wietske van der Zwaag
- Centre d'Imagerie BioMédicale, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Rolf Gruetter
- Laboratory for Functional and Metabolic Imaging, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland.,Department of Radiology, University of Lausanne, Lausanne, Switzerland.,Centre d'Imagerie BioMédicale, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland.,Department of Radiology, University of Geneva, Geneva, Switzerland
| | - José P Marques
- Centre d'Imagerie BioMédicale, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland.,Donders Institute for Brain, Cognition and Behaviour, Radboud University, Nijmegen, The Netherlands
| |
Collapse
|
24
|
Untenberger M, Tan Z, Voit D, Joseph AA, Roeloffs V, Merboldt KD, Schätz S, Frahm J. Advances in real-time phase-contrast flow MRI using asymmetric radial gradient echoes. Magn Reson Med 2015; 75:1901-8. [PMID: 26096085 DOI: 10.1002/mrm.25696] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2015] [Revised: 02/17/2015] [Accepted: 02/24/2015] [Indexed: 11/11/2022]
Abstract
PURPOSE To provide multidimensional velocity compensation for real-time phase-contrast flow MRI. METHODS The proposed method introduces asymmetric gradient echoes for highly undersampled radial FLASH MRI with phase-sensitive image reconstruction by regularized nonlinear inversion (NLINV). Using an adapted gradient delay correction the resulting image quality was analyzed by simulations and experimentally validated at 3 Tesla. For real-time flow MRI the reduced gradient-echo timing allowed for the incorporation of velocity-compensating waveforms for all imaging gradients at even shorter repetition times. RESULTS The results reveal a usable degree of 20% asymmetry. Real-time flow MRI with full velocity compensation eliminated signal void in a flow phantom, confirmed flow parameters in healthy subjects and demonstrated signal recovery and phase conservation in a patient with aortic valve insufficiency and stenosis. Exemplary protocols at 1.4-1.5 mm resolution and 6 mm slice thickness achieved total acquisition times of 33.3-35.7 ms for two images (7 spokes each) with and without flow-encoding gradient. CONCLUSION Asymmetric gradient echoes were successfully implemented for highly undersampled radial trajectories. The resulting temporal gain offers full velocity compensation for real-time phase-contrast flow MRI which minimizes false-positive contributions from complex flow and further enhances the temporal resolution compared with acquisitions with symmetric echoes.
Collapse
Affiliation(s)
- Markus Untenberger
- Biomedizinische NMR Forschungs GmbH am Max-Planck-Institut für biophysikalische Chemie, Göttingen, Germany
| | - Zhengguo Tan
- Biomedizinische NMR Forschungs GmbH am Max-Planck-Institut für biophysikalische Chemie, Göttingen, Germany
| | - Dirk Voit
- Biomedizinische NMR Forschungs GmbH am Max-Planck-Institut für biophysikalische Chemie, Göttingen, Germany
| | - Arun A Joseph
- Biomedizinische NMR Forschungs GmbH am Max-Planck-Institut für biophysikalische Chemie, Göttingen, Germany.,DZHK, German Center for Cardiovascular Research, partner site Göttingen, Germany
| | - Volkert Roeloffs
- Biomedizinische NMR Forschungs GmbH am Max-Planck-Institut für biophysikalische Chemie, Göttingen, Germany
| | - K Dietmar Merboldt
- Biomedizinische NMR Forschungs GmbH am Max-Planck-Institut für biophysikalische Chemie, Göttingen, Germany
| | - Sebastian Schätz
- Biomedizinische NMR Forschungs GmbH am Max-Planck-Institut für biophysikalische Chemie, Göttingen, Germany
| | - Jens Frahm
- Biomedizinische NMR Forschungs GmbH am Max-Planck-Institut für biophysikalische Chemie, Göttingen, Germany.,DZHK, German Center for Cardiovascular Research, partner site Göttingen, Germany
| |
Collapse
|
25
|
Haldar JP, Zhuo J. P-LORAKS: Low-rank modeling of local k-space neighborhoods with parallel imaging data. Magn Reson Med 2015; 75:1499-514. [PMID: 25952136 DOI: 10.1002/mrm.25717] [Citation(s) in RCA: 97] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/21/2015] [Revised: 02/25/2015] [Accepted: 03/13/2015] [Indexed: 11/06/2022]
Abstract
PURPOSE To propose and evaluate P-LORAKS a new calibrationless parallel imaging reconstruction framework. THEORY AND METHODS LORAKS is a flexible and powerful framework that was recently proposed for constrained MRI reconstruction. LORAKS was based on the observation that certain matrices constructed from fully sampled k-space data should have low rank whenever the image has limited support or smooth phase, and made it possible to accurately reconstruct images from undersampled or noisy data using low-rank regularization. This paper introduces P-LORAKS, which extends LORAKS to the context of parallel imaging. This is achieved by combining the LORAKS matrices from different channels to yield a larger but more parsimonious low-rank matrix model of parallel imaging data. This new model can be used to regularize the reconstruction of undersampled parallel imaging data, and implicitly imposes phase, support, and parallel imaging constraints without needing to calibrate phase, support, or sensitivity profiles. RESULTS The capabilities of P-LORAKS are evaluated with retrospectively undersampled data and compared against existing parallel MRI reconstruction methods. Results show that P-LORAKS can improve parallel imaging reconstruction quality, and can enable the use of new k-space trajectories that are not compatible with existing reconstruction methods. CONCLUSION The P-LORAKS framewok provides a new and effective way to regularize parallel imaging reconstruction.
Collapse
Affiliation(s)
- Justin P Haldar
- Department of Electrical Engineering, University of Southern California, Los Angeles, California, USA
| | - Jingwei Zhuo
- Department of Electronic Engineering, Tsinghua University, Beijing, China
| |
Collapse
|
26
|
Blaimer M, Heim M, Neumann D, Jakob PM, Kannengiesser S, Breuer FA. Comparison of phase-constrained parallel MRI approaches: Analogies and differences. Magn Reson Med 2015; 75:1086-99. [PMID: 25845973 DOI: 10.1002/mrm.25685] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2014] [Revised: 02/12/2015] [Accepted: 02/13/2015] [Indexed: 12/12/2022]
Abstract
PURPOSE Phase-constrained parallel MRI approaches have the potential for significantly improving the image quality of accelerated MRI scans. The purpose of this study was to investigate the properties of two different phase-constrained parallel MRI formulations, namely the standard phase-constrained approach and the virtual conjugate coil (VCC) concept utilizing conjugate k-space symmetry. METHODS Both formulations were combined with image-domain algorithms (SENSE) and a mathematical analysis was performed. Furthermore, the VCC concept was combined with k-space algorithms (GRAPPA and ESPIRiT) for image reconstruction. In vivo experiments were conducted to illustrate analogies and differences between the individual methods. Furthermore, a simple method of improving the signal-to-noise ratio by modifying the sampling scheme was implemented. RESULTS For SENSE, the VCC concept was mathematically equivalent to the standard phase-constrained formulation and therefore yielded identical results. In conjunction with k-space algorithms, the VCC concept provided more robust results when only a limited amount of calibration data were available. Additionally, VCC-GRAPPA reconstructed images provided spatial phase information with full resolution. CONCLUSIONS Although both phase-constrained parallel MRI formulations are very similar conceptually, there exist important differences between image-domain and k-space domain reconstructions regarding the calibration robustness and the availability of high-resolution phase information.
Collapse
Affiliation(s)
- Martin Blaimer
- Research Center Magnetic-Resonance-Bavaria (MRB), Würzburg, Germany
| | - Marius Heim
- Department of Experimental Physics 5, University of Würzburg, Würzburg, Germany
| | - Daniel Neumann
- Research Center Magnetic-Resonance-Bavaria (MRB), Würzburg, Germany
| | - Peter M Jakob
- Research Center Magnetic-Resonance-Bavaria (MRB), Würzburg, Germany.,Department of Experimental Physics 5, University of Würzburg, Würzburg, Germany
| | | | - Felix A Breuer
- Research Center Magnetic-Resonance-Bavaria (MRB), Würzburg, Germany
| |
Collapse
|
27
|
Shin T, Kligerman SJ, Crawford RS, Rajagopalan S, Gullapalli RP. Noncontrast-enhanced peripheral venography using velocity-selective magnetization preparation and transient balanced SSFP. Magn Reson Med 2015; 75:653-64. [PMID: 25824323 DOI: 10.1002/mrm.25623] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2014] [Revised: 12/30/2014] [Accepted: 12/31/2014] [Indexed: 12/31/2022]
Abstract
PURPOSE To develop a three-dimensional (3D) noncontrast-enhanced (NCE) peripheral magnetic resonance venography (MRV) method and demonstrate its feasibility in vivo. METHODS The proposed MRV pulse sequence consisted of a velocity-selective (VS) inversion preparation module, inversion delay time (TI), fat inversion pulse, and 3D balanced steady-state free precession (bSSFP) dummy excitations and readout. The VS preparation module inverted arterial blood, which recovered close to zero magnetization during TI. The TI and the number of dummy excitations (Nnum ) were numerically optimized for maximizing vein-to-background contrast and tested in a healthy subject. The proposed MRV of the entire peripheral system, using four-station acquisition, was performed in six healthy subjects and three peripheral artery patients. RESULTS The numerical optimization yielded TI = 350 ms and Ndum = 40, which was supported by the largest vein contrast among the parameters chosen around the optima on in vivo venograms. Four-station peripheral MRV using the optimized parameters well visualized all major deep veins with high vein-to-background contrast. The relative vein contrast ratios were 0.80 ± 0.08, 0.75 ± 0.07, and 0.84 ± 0.06 against the arteries, muscle, and fat, respectively. CONCLUSION The proposed NCE MRV using VS preparation and transient bSSFP can generate high-contrast peripheral venograms directly with a single acquisition.
Collapse
Affiliation(s)
- Taehoon Shin
- Department of Diagnostic Radiology and Nuclear Medicine, University of Maryland, Baltimore, Maryland, USA
| | - Seth J Kligerman
- Department of Diagnostic Radiology and Nuclear Medicine, University of Maryland, Baltimore, Maryland, USA
| | - Robert S Crawford
- Division of Vascular and Endovascular Surgery, University of Maryland, Baltimore, Maryland, USA
| | - Sanjay Rajagopalan
- Division of Cardiovascular Medicine, University of Maryland, Baltimore, Maryland, USA
| | - Rao P Gullapalli
- Department of Diagnostic Radiology and Nuclear Medicine, University of Maryland, Baltimore, Maryland, USA
| |
Collapse
|
28
|
Li G, Hennig J, Raithel E, Büchert M, Paul D, Korvink JG, Zaitsev M. An L1-norm phase constraint for half-Fourier compressed sensing in 3D MR imaging. MAGNETIC RESONANCE MATERIALS IN PHYSICS BIOLOGY AND MEDICINE 2015; 28:459-72. [PMID: 25712732 DOI: 10.1007/s10334-015-0482-7] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/11/2014] [Revised: 12/15/2014] [Accepted: 01/29/2015] [Indexed: 12/22/2022]
Abstract
OBJECTIVE In most half-Fourier imaging methods, explicit phase replacement is used. In combination with parallel imaging, or compressed sensing, half-Fourier reconstruction is usually performed in a separate step. The purpose of this paper is to report that integration of half-Fourier reconstruction into iterative reconstruction minimizes reconstruction errors. MATERIALS AND METHODS The L1-norm phase constraint for half-Fourier imaging proposed in this work is compared with the L2-norm variant of the same algorithm, with several typical half-Fourier reconstruction methods. Half-Fourier imaging with the proposed phase constraint can be seamlessly combined with parallel imaging and compressed sensing to achieve high acceleration factors. RESULTS In simulations and in in-vivo experiments half-Fourier imaging with the proposed L1-norm phase constraint enables superior performance both reconstruction of image details and with regard to robustness against phase estimation errors. CONCLUSION The performance and feasibility of half-Fourier imaging with the proposed L1-norm phase constraint is reported. Its seamless combination with parallel imaging and compressed sensing enables use of greater acceleration in 3D MR imaging.
Collapse
Affiliation(s)
- Guobin Li
- Department of Radiology, Medical Physics, University Medical Center Freiburg, Breisacher Straße 60a, 79106, Freiburg, Germany.
| | - Jürgen Hennig
- Department of Radiology, Medical Physics, University Medical Center Freiburg, Breisacher Straße 60a, 79106, Freiburg, Germany
| | | | - Martin Büchert
- Department of Radiology, Medical Physics, University Medical Center Freiburg, Breisacher Straße 60a, 79106, Freiburg, Germany
| | | | - Jan G Korvink
- Department of Microsystems Engineering-IMTEK, University of Freiburg, Freiburg, Germany
| | - Maxim Zaitsev
- Department of Radiology, Medical Physics, University Medical Center Freiburg, Breisacher Straße 60a, 79106, Freiburg, Germany
| |
Collapse
|
29
|
Blaimer M, Jakob PM, Breuer FA. Regularization method for phase-constrained parallel MRI. Magn Reson Med 2013; 72:166-71. [PMID: 23904349 DOI: 10.1002/mrm.24896] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2013] [Revised: 06/12/2013] [Accepted: 06/29/2013] [Indexed: 11/06/2022]
Abstract
PURPOSE To implement a regularization method for the phase-constrained generalized partially parallel acquisitions (GRAPPA) algorithm to reduce image artifacts caused by data inconsistencies. METHODS Phase-constrained GRAPPA reconstructions are implemented through the use of virtual coils. To that end, synthetic virtual coils are generated by using complex conjugate symmetric signals from the actual coils. Regularization is achieved by applying coefficient-based penalty factors during the GRAPPA calibration procedure. Different penalizing factors are applied for the actual and virtual coils. The method is tested in vivo using T2-weighted turbo spin echo (TSE) images. RESULTS T2 signal decay perturbs conjugate k-space symmetry and produces artifacts in phase-constrained parallel MRI reconstructions of T2-weighted TSE images. Using the proposed regularization method, artifacts are suppressed at the cost of noise amplification. However, there is still a significant SNR gain compared with conventional GRAPPA. CONCLUSION The proposed regularization method is an efficient approach for artifact suppression and maintains the SNR benefit of phase-constrained parallel MRI over conventional parallel MRI.
Collapse
Affiliation(s)
- Martin Blaimer
- Research Center Magnetic-Resonance-Bavaria (MRB), Würzburg, Germany
| | | | | |
Collapse
|
30
|
Schmidt JFM, Wissmann L, Manka R, Kozerke S. Iterative k-t principal component analysis with nonrigid motion correction for dynamic three-dimensional cardiac perfusion imaging. Magn Reson Med 2013; 72:68-79. [PMID: 23913550 DOI: 10.1002/mrm.24894] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2012] [Revised: 06/25/2013] [Accepted: 06/27/2013] [Indexed: 11/08/2022]
Abstract
PURPOSE In this study, an iterative k-t principal component analysis (PCA) algorithm with nonrigid frame-to-frame motion correction is proposed for dynamic contrast-enhanced three-dimensional perfusion imaging. METHODS An iterative k-t PCA algorithm was implemented with regularization using training data corrected for frame-to-frame motion in the x-pc domain. Motion information was extracted using shape-constrained nonrigid image registration of the composite of training and k-t undersampled data. The approach was tested for 10-fold k-t undersampling using computer simulations and in vivo data sets corrupted by respiratory motion artifacts owing to free-breathing or interrupted breath-holds. Results were compared to breath-held reference data. RESULTS Motion-corrected k-t PCA image reconstruction resolved residual aliasing. Signal intensity curves extracted from the myocardium were close to those obtained from the breath-held reference. Upslopes were found to be more homogeneous in space when using the k-t PCA approach with motion correction. CONCLUSIONS Iterative k-t PCA with nonrigid motion correction permits correction of respiratory motion artifacts in three-dimensional first-pass myocardial perfusion imaging.
Collapse
Affiliation(s)
- Johannes F M Schmidt
- Institute for Biomedical Engineering, University and ETH Zurich, Zurich, Switzerland
| | | | | | | |
Collapse
|
31
|
Blaimer M, Choli M, Jakob PM, Griswold MA, Breuer FA. Multiband phase-constrained parallel MRI. Magn Reson Med 2013; 69:974-80. [PMID: 23440994 DOI: 10.1002/mrm.24685] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2012] [Accepted: 01/18/2013] [Indexed: 11/08/2022]
Abstract
PURPOSE Parallel MRI methods are typically associated with a degradation of the signal-to-noise ratio (SNR). High scan time reduction factors are therefore restricted to applications with high intrinsic SNR. One possibility to increase the intrinsic SNR is to simultaneously excite several slices by means of multiband radio-frequency (RF) pulses and subsequently separate the slices by parallel MRI reconstruction algorithms. However, the separation of closely spaced slices may suffer from severe noise amplification when there is insufficient coil sensitivity variation along the slice direction. The purpose of this work is to apply a phase-constrained reconstruction for multiband experiments to minimize the noise amplification. METHODS Pre-defined phase differences between neighboring slices are induced and slice separation is performed by a phase-constrained parallel MRI reconstruction. Phase differences between neighboring slices are tailored to achieve optimal slice separation with minimized noise amplification. The potential of the method is demonstrated through multiband in-vivo experiments. RESULTS Noise amplification in multiband phase-constrained reconstructions is significantly reduced in comparison to standard multiband reconstruction when the phase difference between neighboring slices (distance = 12 mm) is 90°. CONCLUSIONS Multiband phase-constrained parallel MRI has the potential for accelerated multi-slice imaging with an improved SNR performance.
Collapse
Affiliation(s)
- Martin Blaimer
- Research Center Magnetic-Resonance-Bavaria (MRB), Würzburg, Germany.
| | | | | | | | | |
Collapse
|
32
|
Haldar JP, Wedeen VJ, Nezamzadeh M, Dai G, Weiner MW, Schuff N, Liang ZP. Improved diffusion imaging through SNR-enhancing joint reconstruction. Magn Reson Med 2013; 69:277-89. [PMID: 22392528 PMCID: PMC3407310 DOI: 10.1002/mrm.24229] [Citation(s) in RCA: 72] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2011] [Revised: 12/20/2011] [Accepted: 02/06/2012] [Indexed: 11/09/2022]
Abstract
Quantitative diffusion imaging is a powerful technique for the characterization of complex tissue microarchitecture. However, long acquisition times and limited signal-to-noise ratio represent significant hurdles for many in vivo applications. This article presents a new approach to reduce noise while largely maintaining resolution in diffusion weighted images, using a statistical reconstruction method that takes advantage of the high level of structural correlation observed in typical datasets. Compared to existing denoising methods, the proposed method performs reconstruction directly from the measured complex k-space data, allowing for gaussian noise modeling and theoretical characterizations of the resolution and signal-to-noise ratio of the reconstructed images. In addition, the proposed method is compatible with many different models of the diffusion signal (e.g., diffusion tensor modeling and q-space modeling). The joint reconstruction method can provide significant improvements in signal-to-noise ratio relative to conventional reconstruction techniques, with a relatively minor corresponding loss in image resolution. Results are shown in the context of diffusion spectrum imaging tractography and diffusion tensor imaging, illustrating the potential of this signal-to-noise ratio-enhancing joint reconstruction approach for a range of different diffusion imaging experiments.
Collapse
Affiliation(s)
- Justin P Haldar
- Department of Electrical and Computer Engineering, Beckman Institute for Advanced Science and Technology, University of Illinois, Urbana, Illinois, USA.
| | | | | | | | | | | | | |
Collapse
|
33
|
Zhao F, Noll DC, Nielsen JF, Fessler JA. Separate magnitude and phase regularization via compressed sensing. IEEE TRANSACTIONS ON MEDICAL IMAGING 2012; 31:1713-23. [PMID: 22552571 PMCID: PMC3545284 DOI: 10.1109/tmi.2012.2196707] [Citation(s) in RCA: 45] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
Compressed sensing (CS) has been used for accelerating magnetic resonance imaging acquisitions, but its use in applications with rapid spatial phase variations is challenging, e.g., proton resonance frequency shift (PRF-shift) thermometry and velocity mapping. Previously, an iterative MRI reconstruction with separate magnitude and phase regularization was proposed for applications where magnitude and phase maps are both of interest, but it requires fully sampled data and unwrapped phase maps. In this paper, CS is combined into this framework to reconstruct magnitude and phase images accurately from undersampled data. Moreover, new phase regularization terms are proposed to accommodate phase wrapping and to reconstruct images with encoded phase variations, e.g., PRF-shift thermometry and velocity mapping. The proposed method is demonstrated with simulated thermometry data and in vivo velocity mapping data and compared to conventional phase corrected CS.
Collapse
Affiliation(s)
- Feng Zhao
- Biomedical Engineering Department, The University of Michigan, Ann Arbor, MI 48109 USA
| | - Douglas C. Noll
- Biomedical Engineering Department, The University of Michigan, Ann Arbor, MI 48109 USA
| | - Jon-Fredrik Nielsen
- Biomedical Engineering Department, The University of Michigan, Ann Arbor, MI 48109 USA
| | - Jeffrey A. Fessler
- Department of Electrical Engineering and Computer Science, The University of Michigan, Ann Arbor, MI 48109 USA
| |
Collapse
|
34
|
Yankam Njiwa J, Baltes C, Rudin M. Analytic image concept combined to SENSE reconstruction. MAGMA (NEW YORK, N.Y.) 2011; 24:305-313. [PMID: 21833790 DOI: 10.1007/s10334-011-0274-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/04/2010] [Revised: 06/01/2011] [Accepted: 07/13/2011] [Indexed: 05/31/2023]
Abstract
OBJECT Two approaches of reconstructing undersampled partial k-space data, acquired with multiple coils are compared: homodyne detection combined with SENSE (HM_SENSE) and analytic image reconstruction combined with SENSE (AI_SENSE). The latter overcomes limitations of HM_ SENSE by considering aliased images as analytic thus avoiding the need for phase correction required for HM_SENSE. MATERIALS AND METHODS In vivo imaging experiments were carried out in male Lewis rats using both gradient echo and spin echo sequences. Accelerated images obtained by using the various reconstruction algorithms were compared to fully sampled reference images both qualitatively and quantitatively. RESULTS For the various sampling patterns evaluated, both HM_SENSE and AI_SENSE were found to yield robust image reconstruction with small deviations from the reference image. Even for high acceleration factors AI_SENSE still provided useful results and was found superior compared to HM_SENSE. CONCLUSION Combination of partial k-space sampling and parallel image acquisition allows for further acceleration of data acquisition as compared to each method alone. Image reconstruction from undersampled data sets using the AI_SENSE algorithm was found to considerably reduce reconstruction errors and artifacts observed for HM_SENSE reconstruction caused by errors in phase estimation.
Collapse
Affiliation(s)
- Josiane Yankam Njiwa
- Institute for Biomedical Engineering, University and ETH Zurich, AIC-ETH, HIT E22, Wolfgang Pauli-Strasse 27, 8093, Zurich, Switzerland.
| | | | | |
Collapse
|
35
|
Jung Y, Samsonov AA, Bydder M, Block WF. Self-calibrated multiple-echo acquisition with radial trajectories using the conjugate gradient method (SMART-CG). J Magn Reson Imaging 2011; 33:980-7. [PMID: 21448967 DOI: 10.1002/jmri.22482] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022] Open
Abstract
PURPOSE To remove phase inconsistencies between multiple echoes, an algorithm using a radial acquisition to provide inherent phase and magnitude information for self correction was developed. The information also allows simultaneous support for parallel imaging for multiple coil acquisitions. MATERIALS AND METHODS Without a separate field map acquisition, a phase estimate from each echo in multiple echo train was generated. When using a multiple channel coil, magnitude and phase estimates from each echo provide in vivo coil sensitivities. An algorithm based on the conjugate gradient method uses these estimates to simultaneously remove phase inconsistencies between echoes, and in the case of multiple coil acquisition, simultaneously provides parallel imaging benefits. The algorithm is demonstrated on single channel, multiple channel, and undersampled data. RESULTS Substantial image quality improvements were demonstrated. Signal dropouts were completely removed and undersampling artifacts were well suppressed. CONCLUSION The suggested algorithm is able to remove phase cancellation and undersampling artifacts simultaneously and to improve image quality of multiecho radial imaging, the important technique for fast three-dimensional MRI data acquisition.
Collapse
Affiliation(s)
- Youngkyoo Jung
- Department of Radiology, Wake Forest University School of Medicine, Winston-Salem, North Carolina, USA.
| | | | | | | |
Collapse
|
36
|
Wilm BJ, Barmet C, Pavan M, Pruessmann KP. Higher order reconstruction for MRI in the presence of spatiotemporal field perturbations. Magn Reson Med 2011; 65:1690-701. [PMID: 21520269 DOI: 10.1002/mrm.22767] [Citation(s) in RCA: 115] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2010] [Revised: 11/01/2010] [Accepted: 11/24/2010] [Indexed: 11/11/2022]
Abstract
Despite continuous hardware advances, MRI is frequently subject to field perturbations that are of higher than first order in space and thus violate the traditional k-space picture of spatial encoding. Sources of higher order perturbations include eddy currents, concomitant fields, thermal drifts, and imperfections of higher order shim systems. In conventional MRI with Fourier reconstruction, they give rise to geometric distortions, blurring, artifacts, and error in quantitative data. This work describes an alternative approach in which the entire field evolution, including higher order effects, is accounted for by viewing image reconstruction as a generic inverse problem. The relevant field evolutions are measured with a third-order NMR field camera. Algebraic reconstruction is then formulated such as to jointly minimize artifacts and noise in the resulting image. It is solved by an iterative conjugate-gradient algorithm that uses explicit matrix-vector multiplication to accommodate arbitrary net encoding. The feasibility and benefits of this approach are demonstrated by examples of diffusion imaging. In a phantom study, it is shown that higher order reconstruction largely overcomes variable image distortions that diffusion gradients induce in EPI data. In vivo experiments then demonstrate that the resulting geometric consistency permits straightforward tensor analysis without coregistration.
Collapse
Affiliation(s)
- Bertram J Wilm
- Institute for Biomedical Engineering, ETH Zurich and University of Zurich, Zurich, Switzerland
| | | | | | | |
Collapse
|
37
|
Bydder M, Yokoo T, Yu H, Carl M, Reeder SB, Sirlin CB. Constraining the initial phase in water-fat separation. Magn Reson Imaging 2010; 29:216-21. [PMID: 21159457 DOI: 10.1016/j.mri.2010.08.011] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2010] [Revised: 08/03/2010] [Accepted: 08/27/2010] [Indexed: 10/18/2022]
Abstract
An algorithm is described for use in chemical shift-based water-fat separation to constrain the phase of both species to be equal at an echo time of zero. This constraint is physically reasonable since the initial phase should be a property of the excitation pulse and receiver coil only. The advantages of phase constrained water-fat separation, namely, improved noise performance and/or reduced data requirements (fewer echos), are demonstrated in simulations and experiments.
Collapse
Affiliation(s)
- Mark Bydder
- Department of Radiology, University of California San Diego, San Diego, CA, USA.
| | | | | | | | | | | |
Collapse
|
38
|
Samsonov AA, Velikina J, Jung Y, Kholmovski EG, Johnson CR, Block WF. POCS-enhanced correction of motion artifacts in parallel MRI. Magn Reson Med 2010; 63:1104-10. [PMID: 20373413 DOI: 10.1002/mrm.22254] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
A new method for correction of MRI motion artifacts induced by corrupted k-space data, acquired by multiple receiver coils such as phased arrays, is presented. In our approach, a projections onto convex sets (POCS)-based method for reconstruction of sensitivity encoded MRI data (POCSENSE) is employed to identify corrupted k-space samples. After the erroneous data are discarded from the dataset, the artifact-free images are restored from the remaining data using coil sensitivity profiles. The error detection and data restoration are based on informational redundancy of phased-array data and may be applied to full and reduced datasets. An important advantage of the new POCS-based method is that, in addition to multicoil data redundancy, it can use a priori known properties about the imaged object for improved MR image artifact correction. The use of such information was shown to improve significantly k-space error detection and image artifact correction. The method was validated on data corrupted by simulated and real motion such as head motion and pulsatile flow.
Collapse
Affiliation(s)
- Alexey A Samsonov
- Department of Radiology, University of Wisconsin, Madison, Wisconsin, USA.
| | | | | | | | | | | |
Collapse
|
39
|
Du J. Contrast-enhanced MR angiography using time resolved interleaved projection sampling with three-dimensional cartesian phase and slice encoding (TRIPPS). Magn Reson Med 2009; 61:918-24. [DOI: 10.1002/mrm.21805] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
40
|
Improved matrix inversion in image plane parallel MRI. Magn Reson Imaging 2009; 27:942-53. [PMID: 19269768 DOI: 10.1016/j.mri.2009.01.017] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2008] [Revised: 11/11/2008] [Accepted: 01/15/2009] [Indexed: 11/24/2022]
Abstract
A new 3D parallel magnetic resonance imaging (MRI) method named Generalized Unaliasing Incorporating Support constraint and sensitivity Encoding (GUISE) is presented. GUISE allows direct image recovery from arbitrary Cartesian k-space trajectories. However, periodic k-space sampling patterns are considered for reconstruction efficiency. Image recovery methods such as 2D SENSE (SENSitivity Encoding) and 2D CAIPIRINHA (Controlled Aliasing In Parallel Imaging Results IN Higher Acceleration) are special instances of GUISE where specific restrictions are placed on the k-space sampling patterns used. It is shown that the sampling pattern has large impacts on the image reconstruction error due to noise. An efficient sampling pattern design method that incorporates prior knowledge of object support and coil sensitivity profile is proposed. It requires no experimental trials and could be used in clinical imaging. Comparison of the proposed sampling pattern design method with 2D SENSE and 2D CAIPIRINHA are made based on both simulation and experiment results. It is seen that this new adaptive sampling pattern design method results in a lower noise level in reconstructions due to better exploitation of the coil sensitivity variation and object support constraint. In addition, elimination of the non-object region from reconstruction potentially allows an acceleration factor higher than the number of receiver coils used.
Collapse
|
41
|
Han M, Daniel BL, Hargreaves BA. Accelerated bilateral dynamic contrast-enhanced 3D spiral breast MRI using TSENSE. J Magn Reson Imaging 2009; 28:1425-34. [PMID: 19025951 DOI: 10.1002/jmri.21427] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022] Open
Abstract
PURPOSE To assess the ability of adaptive sensitivity encoding incorporating temporal filtering (TSENSE) to accelerate bilateral dynamic contrast-enhanced (DCE) 3D breast MRI. MATERIALS AND METHODS Bilateral DCE breast magnetic resonance imaging (MRI) exams were performed using a dual-band water-only excitation and a "stack-of-spirals" imaging trajectory. TSENSE was applied in the slab direction with an acceleration factor of 2. Four different techniques for sensitivity map calculation were compared by analyzing resultant contrast uptake curves qualitatively and quantitatively for 10 patient datasets. In addition, image quality and temporal resolution were compared between unaccelerated and TSENSE images. RESULTS TSENSE can increase temporal resolution by a factor of 2 in DCE imaging, providing better depiction of contrast uptake curves and good image quality. Of the different methods tested, calculation of static sensitivity maps by averaging late postcontrast frames yields the lowest aliasing artifact level based on ROI analysis. CONCLUSION TSENSE acceleration combined with 3D spiral imaging is very time-efficient, providing 11-second temporal resolution and 1.1 x 1.1 x 3 mm(3) spatial resolution over a 20 x 20 x 10 cm(3) field of view for each breast.
Collapse
Affiliation(s)
- Misung Han
- Lucas MRS/I Center, Department of Radiology, Stanford University, 1201 Welch Rd., Stanford, CA 94305-5488, USA.
| | | | | |
Collapse
|
42
|
Zhang L, Kholmovski EG, Guo J, Choi SEK, Morrell GR, Parker DL. HASTE sequence with parallel acquisition and T2 decay compensation: application to carotid artery imaging. Magn Reson Imaging 2009; 27:13-22. [PMID: 18650045 PMCID: PMC2630054 DOI: 10.1016/j.mri.2008.05.009] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2007] [Revised: 04/06/2008] [Accepted: 05/20/2008] [Indexed: 10/21/2022]
Abstract
T2-weighted carotid artery images acquired using the turbo spin-echo (TSE) sequence frequently suffer from motion artifacts due to respiration and blood pulsation. The possibility of using HASTE sequence to achieve motion-free carotid images was investigated. The HASTE sequence suffers from severe blurring artifacts due to signal loss in later echoes due to T2 decay. Combining HASTE with parallel acquisition (PHASTE) decreases the number of echoes acquired and thus effectively reduces the blurring artifact caused by T2 relaxation. Further improvement in image sharpness can be achieved by performing T2 decay compensation before reconstructing the PHASTE data. Preliminary results have shown successful suppression of motion artifacts with PHASTE imaging. The image quality was enhanced relative to the original HASTE image, but was still less sharp than a non-motion-corrupted TSE image.
Collapse
Affiliation(s)
- Ling Zhang
- Department of Electrical and Computer Engineering, University of Utah, Salt Lake City, UT 84108, USA.
| | | | | | | | | | | |
Collapse
|
43
|
Blaimer M, Gutberlet M, Kellman P, Breuer FA, Köstler H, Griswold MA. Virtual coil concept for improved parallel MRI employing conjugate symmetric signals. Magn Reson Med 2008; 61:93-102. [PMID: 19097211 DOI: 10.1002/mrm.21652] [Citation(s) in RCA: 70] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Martin Blaimer
- Department of Radiology, University Hospitals of Cleveland, Case Western Reserve University, Cleveland, Ohio, USA.
| | | | | | | | | | | |
Collapse
|
44
|
Chen NK, Oshio K, Panych LP. Improved image reconstruction for partial Fourier gradient-echo echo-planar imaging (EPI). Magn Reson Med 2008; 59:916-24. [PMID: 18383294 DOI: 10.1002/mrm.21529] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The partial Fourier gradient-echo echo planar imaging (EPI) technique makes it possible to acquire high-resolution functional MRI (fMRI) data at an optimal echo time. This technique is especially important for fMRI studies at high magnetic fields, where the optimal echo time is short and may not be achieved with a full Fourier acquisition scheme. In addition, it has been shown that partial Fourier EPI provides better anatomic resolvability than full Fourier EPI. However, the partial Fourier gradient-echo EPI may be degraded by artifacts that are not usually seen in other types of imaging. Those unique artifacts in partial Fourier gradient-echo EPI, to our knowledge, have not yet been systematically evaluated. Here we use the k-space energy spectrum analysis method to understand and characterize two types of partial Fourier EPI artifacts. Our studies show that Type 1 artifact, originating from k-space energy loss, cannot be corrected with pure postprocessing, and Type 2 artifact can be eliminated with an improved reconstruction method. We propose a novel algorithm, that combines images obtained from two or more reconstruction schemes guided by k-space energy spectrum analysis, to generate partial Fourier EPI with greatly reduced Type 2 artifact. Quality control procedures for avoiding Type 1 artifact in partial Fourier EPI are also discussed.
Collapse
Affiliation(s)
- Nan-kuei Chen
- Brain Imaging and Analysis Center and Department of Radiology, Duke University, Durham, North Carolina 27710, USA.
| | | | | |
Collapse
|
45
|
Lew C, Pineda AR, Clayton D, Spielman D, Chan F, Bammer R. SENSE phase-constrained magnitude reconstruction with iterative phase refinement. Magn Reson Med 2008; 58:910-21. [PMID: 17969127 DOI: 10.1002/mrm.21284] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
Conventional sensitivity encoding (SENSE) reconstruction is based on equations in the complex domain. However, for many MRI applications only the magnitude is relevant. If there exists an estimate of the underlying phase information, a magnitude-only phase-constrained reconstruction can help to improve the conditioning of the SENSE reconstruction problem. Consequently, this reduces g-factor-related noise enhancement. In previous attempts at phase-constrained SENSE reconstruction, image quality was hampered by strong aliasing artifacts resulting from inadequate phase estimates and high sensitivity to phase errors. If a full-resolution phase image is used, a significant reduction in aliasing errors and better noise properties compared to SENSE can be obtained. An iterative scheme that improves the phase estimate to better approximate the phase is presented. The mathematical framework of the new approach is provided together with comparisons of conventional SENSE, phase-constrained SENSE, and the new phase-refinement method. Both theory and experimental verification demonstrate significantly better noise performance at high reduction factors, i.e., close to the theoretical limit. For applications that need only magnitude data, an iterative phase-constrained SENSE reconstruction can provide substantial SNR improvement over SENSE reconstruction and less artifacts than phase-constrained SENSE.
Collapse
Affiliation(s)
- Calvin Lew
- Lucas MRS/I Center, Stanford University, Stanford, California 94035, USA.
| | | | | | | | | | | |
Collapse
|
46
|
Abstract
Parallel imaging has been the single biggest innovation in magnetic resonance imaging in the last decade. The use of multiple receiver coils to augment the time consuming Fourier encoding has reduced acquisition times significantly. This increase in speed comes at a time when other approaches to acquisition time reduction were reaching engineering and human limits. A brief summary of spatial encoding in MRI is followed by an introduction to the problem parallel imaging is designed to solve. There are a large number of parallel reconstruction algorithms; this article reviews a cross-section, SENSE, SMASH, g-SMASH and GRAPPA, selected to demonstrate the different approaches. Theoretical (the g-factor) and practical (coil design) limits to acquisition speed are reviewed. The practical implementation of parallel imaging is also discussed, in particular coil calibration. How to recognize potential failure modes and their associated artefacts are shown. Well-established applications including angiography, cardiac imaging and applications using echo planar imaging are reviewed and we discuss what makes a good application for parallel imaging. Finally, active research areas where parallel imaging is being used to improve data quality by repairing artefacted images are also reviewed.
Collapse
Affiliation(s)
- David J Larkman
- The Imaging Sciences Department, Clinical Sciences Centre, Faculty of Medicine, Imperial College London, Hammersmith Hospital Campus, Du Cane Road, London W12 0NN, UK.
| | | |
Collapse
|
47
|
Bydder M, Perthen JE, Du J. Optimization of sensitivity encoding with arbitrary k-space trajectories. Magn Reson Imaging 2007; 25:1123-9. [PMID: 17905244 DOI: 10.1016/j.mri.2007.01.003] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2006] [Revised: 12/24/2006] [Accepted: 01/05/2007] [Indexed: 11/25/2022]
Abstract
Sensitivity encoding (SENSE) is a magnetic resonance technique that unifies gradient and receive coil encoding. SENSE reconstructs the image by solving a large, ill-conditioned inverse problem, which generally requires regularization and preconditioning. The present study suggests a simple heuristic for determining the regularization parameter. Also discussed are the use of density weighting and intensity correction as preconditioners and the role that coil sensitivity estimation has in regularization. A modification to the intensity correction is proposed for use with a phase constraint.
Collapse
Affiliation(s)
- Mark Bydder
- Department of Radiology, MRI3 Research, University of California San Diego, San Diego, CA 92103-8226, USA.
| | | | | |
Collapse
|
48
|
Lin FH, Wald LL, Ahlfors SP, Hämäläinen MS, Kwong KK, Belliveau JW. Dynamic magnetic resonance inverse imaging of human brain function. Magn Reson Med 2007; 56:787-802. [PMID: 16964616 DOI: 10.1002/mrm.20997] [Citation(s) in RCA: 79] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
MRI is widely used for noninvasive hemodynamic-based functional brain imaging. In traditional spatial encoding, however, gradient switching limits the temporal resolution, which makes it difficult to unambiguously identify possible fast nonhemodynamic changes. In this paper we propose a novel reconstruction approach, called dynamic inverse imaging (InI), that is capable of providing millisecond temporal resolution when highly parallel detection is used. To achieve an order-of-magnitude speedup in generating time-resolved contrast estimates and dynamic statistical parametric maps (dSPMs), the spatial information is derived from an array of detectors rather than by time-consuming gradient-encoding methods. The InI approach was inspired by electroencephalography (EEG) and magnetoencephalography (MEG) source localization techniques. Dynamic MR InI was evaluated by means of numerical simulations. InI was also applied to measure BOLD hemodynamic time curves at 20-ms temporal resolution in a visual stimulation experiment using a 90-channel head array. InI is expected to improve the time resolution of MRI and provide increased flexibility in the trade-off between spatial and temporal resolution for studies of dynamic activation patterns in the human brain.
Collapse
Affiliation(s)
- Fa-Hsuan Lin
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, Massachusetts 02129, USA.
| | | | | | | | | | | |
Collapse
|
49
|
Agarwal H, Abd-elmoniem K, Prince J. FAST MOTION IMAGING USING REDUCED FIELD OF VIEW PARTIAL FOURIER MRI. 2007 4TH IEEE INTERNATIONAL SYMPOSIUM ON BIOMEDICAL IMAGING: FROM NANO TO MACRO 2007. [DOI: 10.1109/isbi.2007.356928] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/01/2023]
|
50
|
Luo J, Zhu Y. Magnetic resonance image reconstruction using equation system solving-based singularity function analysis. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2007; 2007:751-754. [PMID: 18002065 DOI: 10.1109/iembs.2007.4352399] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/25/2023]
Abstract
This paper presents a method of reconstructing magnetic resonance (MR) images from partial k-space using a so-called singularity function model. To estimate the parameters of the model, we propose a new technique based on equation system solving. The proposed reconstruction method is evaluated on both simulated and real brain MR data. It is also compared with existing techniques. The results show that this new strategy exhibits a great simplicity while offering a good reconstruction quality.
Collapse
Affiliation(s)
- Jianhua Luo
- College of Life Science and Technology, Shanghai, Jiaotong University, Shanghai, P.R.China
| | | |
Collapse
|