1
|
Wang G, Fessler JA. Efficient Approximation of Jacobian Matrices Involving a Non-Uniform Fast Fourier Transform (NUFFT). IEEE TRANSACTIONS ON COMPUTATIONAL IMAGING 2023; 9:43-54. [PMID: 37090025 PMCID: PMC10118239 DOI: 10.1109/tci.2023.3240081] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
There is growing interest in learning Fourier domain sampling strategies (particularly for magnetic resonance imaging, MRI) using optimization approaches. For non-Cartesian sampling, the system models typically involve non-uniform fast Fourier transform (NUFFT) operations. Commonly used NUFFT algorithms contain frequency domain interpolation, which is not differentiable with respect to the sampling pattern, complicating the use of gradient methods. This paper describes an efficient and accurate approach for computing approximate gradients involving NUFFTs. Multiple numerical experiments validate the improved accuracy and efficiency of the proposed approximation. As an application to computational imaging, the NUFFT Jacobians were used to optimize non-Cartesian MRI sampling trajectories via data-driven stochastic optimization. Specifically, the sampling patterns were learned with respect to various model-based image reconstruction (MBIR) algorithms. The proposed approach enables sampling optimization for image sizes that are infeasible with standard auto-differentiation methods due to memory limits. The synergistic acquisition and reconstruction design leads to remarkably improved image quality. In fact, we show that model-based image reconstruction methods with suitably optimized imaging parameters can perform nearly as well as CNN-based methods.
Collapse
Affiliation(s)
- Guanhua Wang
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI 48109 USA
| | | |
Collapse
|
2
|
Yang X, Liu M, Duan J, Sun H, An J, Benkert T, Dai H, Wang C. Three-dimensional ultrashort echo time magnetic resonance imaging in assessment of idiopathic pulmonary fibrosis, in comparison with high-resolution computed tomography. Quant Imaging Med Surg 2022; 12:4176-4189. [PMID: 35919053 PMCID: PMC9338383 DOI: 10.21037/qims-21-1133] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2021] [Accepted: 05/27/2022] [Indexed: 11/20/2022]
Abstract
Background We aimed to evaluate the image quality, feasibility, and diagnostic performance of three-dimensional ultrashort echo time magnetic resonance imaging (3D UTE-MRI) to assess idiopathic pulmonary fibrosis (IPF) compared with high-resolution computed tomography (HRCT) and half-Fourier single-shot turbo spin-echo (HASTE) MRI. Methods A total of 36 patients with IPF (34 men; mean age: 62±8 years, age range: 43 to 78 years) were prospectively included and underwent HRCT and chest MRI on the same day. Chest MRI was performed with a free-breathing 3D spiral UTE pulse sequence and HASTE sequence on a 1.5 T MRI. Two radiologists independently evaluated the image quality of the HRCT, HASTE, and 3D UTE-MRI. They assessed the representative imaging features of IPF, including honeycombing, reticulation, traction bronchiectasis, and ground-glass opacities. Image quality of the 3D UTE-MRI, HASTE, and HRCT were assessed using a 5-point visual scoring method. Kappa and weighted kappa analysis were used to measure intra- and inter-observer and inter-method agreements. Sensitivity (SE), specificity (SP), and accuracy (AC) were used to assess the performance of 3D UTE-MRI for detecting image features of IPF and monitoring the extent of pulmonary fibrosis. Linear regressions and Bland-Altman plots were generated to assess the correlation and agreement between the assessment of the extent of pulmonary fibrosis made by the 2 observers. Results The image quality of HRCT was higher than that of HASTE and UTE-MRI (HRCT vs. UTE-MRI vs. HASTE: 4.9±0.3 vs. 4.1±0.7 vs. 3.0±0.3; P<0.001). Interobserver agreement of HRCT, HASTE, and 3D UTE-MRI when assessing pulmonary fibrosis was substantial and excellent (HRCT: 0.727≤ κ ≤1, P<0.001; HASTE: 0.654≤ κ ≤1, P<0.001; 3D UTE-MRI: 0.719≤ κ ≤0.824, P<0.001). In addition, reticulation (SE: 97.1%; SP: 100%; AC: 97.2%; κ =0.654), honeycombing (SE: 83.3%; SP: 100%; AC: 86.1%; κ =0.625) patterns, and traction bronchiectasis (SE: 94.1%; SP: 100%; AC: 94.4%, κ =0.640) were also well-visualized on 3D UTE-MRI, which was significantly superior to HASTE. Compared with HRCT, the sensitivity of 3D UTE-MRI to detect signs of pulmonary fibrosis (n=35) was 97.2%. The interobserver agreement in elevation of the extent of pulmonary fibrosis with HRCT and 3D UTE-MRI was R2=0.84 (P<0.001) and R2=0.84 (P<0.001), respectively. The extent of pulmonary fibrosis assessed with 3D UTE-MRI [median =9, interquartile range (IQR): 6.25 to 10.00] was lower than that from HRCT (median =12, IQR: 9.25 to 13.00; U=320.00, P<0.001); however, they had a positive correlation (R=0.72, P<0.001). Conclusions As a radiation-free non-contrast enhanced imaging method, although the image quality of 3D UTE-MRI is inferior to that of HRCT, it has high reproducibility to identify the imaging features of IPF and evaluate the extent of pulmonary fibrosis.
Collapse
Affiliation(s)
- Xiaoyan Yang
- Department of Pulmonary and Critical Care Medicine, China-Japan Friendship Hospital, Capital Medical University, Beijing, China.,National Center for Respiratory Medicine, Institute of Respiratory Medicine, Chinese Academy of Medical Sciences, National Clinical Research Center for Respiratory Diseases, Beijing, China
| | - Min Liu
- Department of Radiology, China-Japan Friendship Hospital, Beijing, China
| | - Jianghui Duan
- Department of Radiology, China-Japan Friendship Hospital, Beijing, China
| | - Haishuang Sun
- National Center for Respiratory Medicine, Institute of Respiratory Medicine, Chinese Academy of Medical Sciences, National Clinical Research Center for Respiratory Diseases, Beijing, China
| | - Jing An
- Siemens Shenzhen Magnetic Resonance Ltd., Shenzhen, China
| | - Thomas Benkert
- MR Application Predevelopment, Siemens Healthcare GmbH, Erlangen, Germany
| | - Huaping Dai
- Department of Pulmonary and Critical Care Medicine, China-Japan Friendship Hospital, Capital Medical University, Beijing, China.,National Center for Respiratory Medicine, Institute of Respiratory Medicine, Chinese Academy of Medical Sciences, National Clinical Research Center for Respiratory Diseases, Beijing, China
| | - Chen Wang
- Department of Pulmonary and Critical Care Medicine, China-Japan Friendship Hospital, Capital Medical University, Beijing, China.,National Center for Respiratory Medicine, Institute of Respiratory Medicine, Chinese Academy of Medical Sciences, National Clinical Research Center for Respiratory Diseases, Beijing, China.,Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
| |
Collapse
|
3
|
Weigel M, Haas T, Wendebourg MJ, Schlaeger R, Bieri O. Imaging of the thoracic spinal cord using radially sampled averaged magnetization inversion recovery acquisitions. J Neurosci Methods 2020; 343:108825. [PMID: 32580062 DOI: 10.1016/j.jneumeth.2020.108825] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2019] [Revised: 06/18/2020] [Accepted: 06/19/2020] [Indexed: 10/24/2022]
Abstract
BACKGROUND Spinal cord (SC) gray and white matter atrophy quantification by advanced morphometric MRI can help to better characterize the course of neurodegenerative diseases in vivo, such as e.g. lower motor neuron disorders. Imaging the lower thoracic cord - containing those motor neurons that control leg function - could be particularly informative, however, is challenging due to tissue composition, physiological motion and large field of views. NEW METHOD An "averaged magnetization inversion recovery acquisitions" (AMIRA) approach with a radial k-space acquisition scheme was developed. The method is designed for morphometric SC imaging with a focus on the thoracic SC. RESULTS In a typical setting, radial AMIRA acquires transverse slices with a high 0.50 × 0.50mm2 in-plane resolution and a pronounced positive contrast between thoracic gray and white matter, within typically 2:39 min. Additional proof-of-concept measurements in patients demonstrate that such contrast and resolving capability is indeed necessary to assess potential atrophy of the anterior horns. COMPARISON WITH EXISTING METHOD(S) Radial AMIRA utilizes two benefits of radial MRI techniques: being generally less prone to motion effects and that fold over artifacts can manifest less intrusively. These benefits are united with the original AMIRA approach which allows the contrast to be 'tuned' and improved based on the combination of five simultaneously acquired images of different tissue contrast. CONCLUSIONS Radial AMIRA is a promising approach for in vivo SC gray and white matter atrophy visualization and quantification in lower motor neuron diseases and other autoimmune or genetic diseases involving the entire (not only cervical) spinal cord.
Collapse
Affiliation(s)
- Matthias Weigel
- Division of Radiological Physics, Department of Radiology, University Hospital Basel, Petersgraben 4, 4031 Basel, Switzerland; Translational Imaging in Neurology (ThINk) Basel, Department of Medicine and Biomedical Engineering, University Hospital Basel and University of Basel, Petersgraben 4, 4031 Basel, Switzerland; Neurologic Clinic and Policlinic, Departments of Medicine, Clinical Research and Biomedical Engineering, University Hospital Basel and University of Basel, Petersgraben 4, 4031 Basel, Switzerland.
| | - Tanja Haas
- Division of Radiological Physics, Department of Radiology, University Hospital Basel, Petersgraben 4, 4031 Basel, Switzerland
| | - Maria Janina Wendebourg
- Translational Imaging in Neurology (ThINk) Basel, Department of Medicine and Biomedical Engineering, University Hospital Basel and University of Basel, Petersgraben 4, 4031 Basel, Switzerland; Neurologic Clinic and Policlinic, Departments of Medicine, Clinical Research and Biomedical Engineering, University Hospital Basel and University of Basel, Petersgraben 4, 4031 Basel, Switzerland
| | - Regina Schlaeger
- Translational Imaging in Neurology (ThINk) Basel, Department of Medicine and Biomedical Engineering, University Hospital Basel and University of Basel, Petersgraben 4, 4031 Basel, Switzerland; Neurologic Clinic and Policlinic, Departments of Medicine, Clinical Research and Biomedical Engineering, University Hospital Basel and University of Basel, Petersgraben 4, 4031 Basel, Switzerland
| | - Oliver Bieri
- Division of Radiological Physics, Department of Radiology, University Hospital Basel, Petersgraben 4, 4031 Basel, Switzerland; Department of Biomedical Engineering, University of Basel, Gewerbestrasse 14, 4123 Allschwil, Switzerland
| |
Collapse
|
4
|
Post-contrast T1-weighted sequences in pediatric abdominal imaging: comparative analysis of three different sequences and imaging approach. Pediatr Radiol 2014; 44:1258-65. [PMID: 24723237 DOI: 10.1007/s00247-014-2969-1] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/23/2013] [Revised: 02/07/2014] [Accepted: 03/10/2014] [Indexed: 10/25/2022]
Abstract
BACKGROUND Post-contrast T1-weighted imaging is an essential component of a comprehensive pediatric abdominopelvic MR examination. However, consistent good image quality is challenging, as respiratory motion in sedated children can substantially degrade the image quality. OBJECTIVE To compare the image quality of three different post-contrast T1-weighted imaging techniques-standard three-dimensional gradient-echo (3-D-GRE), magnetization-prepared gradient-recall echo (MP-GRE) and 3-D-GRE with radial data sampling (radial 3-D-GRE)-acquired in pediatric patients younger than 5 years of age. MATERIALS AND METHODS Sixty consecutive exams performed in 51 patients (23 females, 28 males; mean age 2.5 ± 1.4 years) constituted the final study population. Thirty-nine scans were performed at 3 T and 21 scans were performed at 1.5 T. Two different reviewers independently and blindly qualitatively evaluated all sequences to determine image quality and extent of artifacts. RESULTS MP-GRE and radial 3-D-GRE sequences had the least respiratory motion (P < 0.0001). Standard 3-D-GRE sequences displayed the lowest average score ratings in hepatic and pancreatic edge definition, hepatic vessel clarity and overall image quality. Radial 3-D-GRE sequences showed the highest scores ratings in overall image quality. CONCLUSIONS Our preliminary results support the preference of fat-suppressed radial 3-D-GRE as the best post-contrast T1-weighted imaging approach for patients under the age of 5 years, when dynamic imaging is not essential.
Collapse
|
5
|
Bamrungchart S, Tantaway EM, Midia EC, Hernandes MA, Srirattanapong S, Dale BM, Semelka RC. Free breathing three-dimensional gradient echo-sequence with radial data sampling (radial 3D-GRE) examination of the pancreas: Comparison with standard 3D-GRE volumetric interpolated breathhold examination (VIBE). J Magn Reson Imaging 2013; 38:1572-7. [PMID: 23417838 DOI: 10.1002/jmri.24064] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2012] [Accepted: 01/10/2013] [Indexed: 02/06/2023] Open
Abstract
PURPOSE To evaluate the diagnostic performance of free breathing three-dimensional gradient echo-sequence with radial data sampling (radial 3D-GRE) in MR imaging of the normal and diseased pancreas, using standard 3D-GRE for comparison in cooperative patients, and to perform a preliminary assessment in noncooperative patients. MATERIALS AND METHODS One hundred and eight consecutive patients underwent 1.5 Tesla MR imaging of the abdomen that included pre- and postcontrast free breathing radial 3D-GRE. The sequences were evaluated by two radiologists retrospectively, independently, and blindly. The results were compared using Wilcoxon-Mann-Whitney test. Kappa statistics were used to measure the extent of agreement between the reviewers. RESULTS The average scores indicated that the overall images quality of radial 3D-GRE was lower than 3D-GRE-VIBE in both pre- and postcontrast study (P = 0.0172 and 0.0001), however it achieved a rating that approximated good. In all patients, radial 3D-GRE had a mild extent of streak artifact, pulsation, susceptibility, and respiratory artifact. Radial 3D-GRE approximated good results for pancreatic edge sharpness and pancreatic ductal clarity, and did not differ significantly between cooperative and noncooperative patients. Respiratory artifact was worse in cooperative than in noncooperative patients (P = 0.02). Demonstration of pancreatic disease was slightly inferior with radial 3D-VIBE, but quality approximated good. CONCLUSION Free breathing radial 3D-GRE may be applicable for pancreatic MR imaging in patients who are unable to suspend respiration.
Collapse
Affiliation(s)
- Saraporn Bamrungchart
- Department of Radiology, University of North Carolina at Chapel Hill, North Carolina, USA; Department of Radiology, Lerdsin hospital, Bangkok, Thailand
| | | | | | | | | | | | | |
Collapse
|
6
|
Feng Y, Song Y, Wang C, Xin X, Feng Q, Chen W. Fast direct fourier reconstruction of radial and PROPELLER MRI data using the chirp transform algorithm on graphics hardware. Magn Reson Med 2012; 70:1087-94. [PMID: 23165973 DOI: 10.1002/mrm.24556] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2012] [Revised: 10/16/2012] [Accepted: 10/17/2012] [Indexed: 11/06/2022]
Abstract
PURPOSE To develop and test a new algorithm for fast direct Fourier transform (DrFT) reconstruction of MR data on non-Cartesian trajectories composed of lines with equally spaced points. THEORY AND METHODS The DrFT, which is normally used as a reference in evaluating the accuracy of other reconstruction methods, can reconstruct images directly from non-Cartesian MR data without interpolation. However, DrFT reconstruction involves substantially intensive computation, which makes the DrFT impractical for clinical routine applications. In this article, the Chirp transform algorithm was introduced to accelerate the DrFT reconstruction of radial and Periodically Rotated Overlapping ParallEL Lines with Enhanced Reconstruction (PROPELLER) MRI data located on the trajectories that are composed of lines with equally spaced points. The performance of the proposed Chirp transform algorithm-DrFT algorithm was evaluated by using simulation and in vivo MRI data. RESULTS After implementing the algorithm on a graphics processing unit, the proposed Chirp transform algorithm-DrFT algorithm achieved an acceleration of approximately one order of magnitude, and the speed-up factor was further increased to approximately three orders of magnitude compared with the traditional single-thread DrFT reconstruction. CONCLUSION Implementation the Chirp transform algorithm-DrFT algorithm on the graphics processing unit can efficiently calculate the DrFT reconstruction of the radial and PROPELLER MRI data.
Collapse
Affiliation(s)
- Yanqiu Feng
- School of Biomedical Engineering, Southern Medical University, Guangzhou, China
| | | | | | | | | | | |
Collapse
|
7
|
|
8
|
Free-breathing 3D T1-weighted gradient-echo sequence with radial data sampling in abdominal MRI: preliminary observations. AJR Am J Roentgenol 2011; 197:650-7. [PMID: 21862807 DOI: 10.2214/ajr.10.5881] [Citation(s) in RCA: 69] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
OBJECTIVE The purposes of this study were to evaluate the feasibility of a free-breathing 3D gradient-recalled echo sequence with radial data sampling (radial 3D GRE) in abdominal MRI compared with a standard 3D GRE volumetric interpolated breath-hold examination (VIBE) sequence for imaging of cooperative patients and to perform a preliminary assessment in imaging of noncooperative patients. MATERIALS AND METHODS Fifty-five consecutively registered patients who underwent unenhanced and contrast-enhanced abdominal MRI with the free-breathing radial 3D GRE technique constituted the study population. Two readers independently and blindly evaluated the images. RESULTS Overall image quality with the contrast-enhanced radial 3D GRE sequence was lower than but rated at least nearly as good as that with the 3D GRE VIBE sequence (p < 0.0001). Higher scores were recorded for 3D GRE VIBE images with respect to pixel graininess, streaking artifact, and sharpness (p = 0.0009 to p < 0.0001). Except for sharpness of vessels on unenhanced images, results for the radial 3D GRE sequence did not differ significantly in the comparison of cooperative and noncooperative patients (p = 0.004). For imaging of noncooperative patients, radial 3D GRE images of children had higher ratings for shading (unenhanced, p = 0.0004; contrast-enhanced, p < 0.0001) and streaking artifacts on contrast-enhanced images (p = 0.0017) than did those of adults. Overall image quality was higher for pediatric patients. In lesion analysis, use of the 3D GRE VIBE sequence was associated with significantly greater detectability, confidence, and conspicuity than was use of the radial 3D GRE sequence (p = 0.00026-0.011). CONCLUSION A free-breathing radial 3D GRE sequence is feasible for abdominal MRI and may find application in imaging of patients who are unable to suspend respiration, especially children.
Collapse
|
9
|
Jung H, Park J, Yoo J, Ye JC. Radial k-t FOCUSS for high-resolution cardiac cine MRI. Magn Reson Med 2010; 63:68-78. [PMID: 19859952 DOI: 10.1002/mrm.22172] [Citation(s) in RCA: 75] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
A compressed sensing dynamic MR technique called k-t FOCUSS (k-t FOCal Underdetermined System Solver) has been recently proposed. It outperforms the conventional k-t BLAST/SENSE (Broad-use Linear Acquisition Speed-up Technique/SENSitivity Encoding) technique by exploiting the sparsity of x-f signals. This paper applies this idea to radial trajectories for high-resolution cardiac cine imaging. Radial trajectories are more suitable for high-resolution dynamic MRI than Cartesian trajectories since there is smaller tradeoff between spatial resolution and number of views if streaking artifacts due to limited views can be resolved. As shown for Cartesian trajectories, k-t FOCUSS algorithm efficiently removes artifacts while preserving high temporal resolution. k-t FOCUSS algorithm applied to radial trajectories is expected to enhance dynamic MRI quality. Rather than using an explicit gridding method, which transforms radial k-space sampling data to Cartesian grid prior to applying k-t FOCUSS algorithms, we use implicit gridding during FOCUSS iterations to prevent k-space sampling errors from being propagated. In addition, motion estimation and motion compensation after the first FOCUSS iteration were used to further sparsify the residual image. By applying an additional k-t FOCUSS step to the residual image, improved resolution was achieved. In vivo experimental results show that this new method can provide high spatiotemporal resolution even from a very limited radial data set.
Collapse
Affiliation(s)
- Hong Jung
- Bio-Imaging and Signal Processing Laboratory, Department of Bio and Brain Engineering, Korea Advanced Institute of Science and Technology (KAIST), Republic of Korea
| | | | | | | |
Collapse
|
10
|
Viard R. Présentation des différentes avancées techniques en imagerie interventionnelle par résonance magnétique. Ing Rech Biomed 2009. [DOI: 10.1016/j.irbm.2008.12.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
11
|
Sorensen TS, Schaeffter T, Noe KO, Hansen MS. Accelerating the nonequispaced fast Fourier transform on commodity graphics hardware. IEEE TRANSACTIONS ON MEDICAL IMAGING 2008; 27:538-547. [PMID: 18390350 DOI: 10.1109/tmi.2007.909834] [Citation(s) in RCA: 51] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
We present a fast parallel algorithm to compute the nonequispaced fast Fourier transform on commodity graphics hardware (the GPU). We focus particularly on a novel implementation of the convolution step in the transform as it was previously its most time consuming part. We describe the performance for two common sample distributions in medical imaging (radial and spiral trajectories), and for different convolution kernels as these parameters all influence the speed of the algorithm. The GPU-accelerated convolution is up to 85 times faster as our reference, the open source NFFT library on a state-of-the-art 64 bit CPU. The accuracy of the proposed GPU implementation was quantitatively evaluated at the various settings. To illustrate the applicability of the transform in medical imaging, in which it is also known as gridding, we look specifically at non-Cartesian magnetic resonance imaging and reconstruct both a numerical phantom and an in vivo cardiac image.
Collapse
Affiliation(s)
- T S Sorensen
- Division of Imaging Sciences at King's College London, London WC2R 2LS, UK.
| | | | | | | |
Collapse
|
12
|
Viard R, Rousseau J. [Interventional MR imaging: state of the art and technological advances]. JOURNAL DE RADIOLOGIE 2008; 89:13-20. [PMID: 18288022 DOI: 10.1016/s0221-0363(08)70365-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/25/2023]
Abstract
Due to its excellent soft tissue contrast and lack of ionizing radiation, MR imaging is well suited for interventional procedures. MRI is being increasingly used for guidance during percutaneous procedures or surgery. Technical advances in interventional MR imaging are reviewed in this paper. Ergonomical factors with improved access to patients as well as advances in informatics, electronics and robotics largely explain this increasing role. Different elements are discussed from improved access to patients in the scanners to improved acquisition pulse sequences. Selected clinical applications and recent publications will be presented to illustrate the current status of this technique.
Collapse
Affiliation(s)
- R Viard
- INSERM, U703, ITM, Pavillon Vancostenobel, CHRU de Lille, 2 avenue Oscar Lambret, 59000 Lille Cedex.
| | | |
Collapse
|
13
|
Pulse sequences and system interfaces for interventional and real-time MRI. J Magn Reson Imaging 2008; 27:267-75. [DOI: 10.1002/jmri.21268] [Citation(s) in RCA: 40] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
14
|
Liang D, Lam EY, Fung GSK. A least squares quantization table method for direct reconstruction of MR images with non-Cartesian trajectory. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2007; 188:141-50. [PMID: 17646119 DOI: 10.1016/j.jmr.2007.06.013] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/26/2006] [Revised: 06/27/2007] [Accepted: 06/27/2007] [Indexed: 05/16/2023]
Abstract
The direct Fourier transform method is a straightforward solution with high accuracy for reconstructing magnetic resonance (MR) images from nonuniformly sampled k-space data, given that the optimal density compensation function is selected and the underlying magnetic field is sufficiently uniform. The computation however is very time-consuming, making it impractical especially for large-size images. In this paper, the least squares quantization table (LSQT) method is proposed to accelerate the direct Fourier transform computation, similar to the recently proposed methods such as using look-up table (LUT) or equal-phase-line (EPL). With LSQT, all the image pixels are first classified into several groups where the Lloyd-Max quantization scheme is used to ensure the minimal classification error. The representative value of each group is stored in a small-size LSQT in advance to reduce the computational load. The pixels in the same group receive the same contribution, which is calculated only once for each group instead of for each pixel, resulting in the reduction of computation because the number of groups is far smaller than the number of pixels. Finally, each image pixel is mapped into the nearest group and its representative value is used to reconstruct the image. The experimental results show that the LSQT method requires far smaller memory size than the LUT method and fewer multiplication operations than the LUT and EPL methods. Moreover, the LSQT method can perform large-size reconstructions that achieve comparable or higher accuracy as compared to the EPL and gridding methods when the appropriate parameters are given. The inherent parallel structure also makes the LSQT method easily adaptable to a multiprocessor system.
Collapse
Affiliation(s)
- Dong Liang
- Department of Electrical and Electronic Engineering, The University of Hong Kong, Pokfulam Road, Hong Kong
| | | | | |
Collapse
|
15
|
Santos JM, Wright GA, Pauly JM. Flexible real-time magnetic resonance imaging framework. CONFERENCE PROCEEDINGS : ... ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL CONFERENCE 2007; 2004:1048-51. [PMID: 17271862 DOI: 10.1109/iembs.2004.1403343] [Citation(s) in RCA: 65] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
The extension of MR imaging to new applications has demonstrated the limitations of the architecture of current real-time systems. Traditional real-time implementations provide continuous acquisition of data and modification of basic sequence parameters on the fly. We have extended the concept of real-time MRI by designing a system that drives the examinations from a real-time localizer and then gets reconfigured for different imaging modes. Upon operator request or automatic feedback the system can immediately generate a new pulse sequence or change fundamental aspects of the acquisition such as gradient waveforms excitation pulses and scan planes. This framework has been implemented by connecting a data processing and control workstation to a conventional clinical scanner. Key components on the design of this framework are the data communication and control mechanisms, reconstruction algorithms optimized for real-time and adaptability, flexible user interface and extensible user interaction. In this paper we describe the various components that comprise this system. Some of the applications implemented in this framework include real-time catheter tracking embedded in high frame rate real-time imaging and immediate switching between real-time localizer and high-resolution volume imaging for coronary angiography applications.
Collapse
|
16
|
Lin HY, Flask CA, Dale BM, Duerk JL. Rapid dark-blood carotid vessel-wall imaging with random bipolar gradients in a radial SSFP acquisition. J Magn Reson Imaging 2007; 25:1299-304. [PMID: 17520723 DOI: 10.1002/jmri.20821] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022] Open
Abstract
PURPOSE To investigate and evaluate a new rapid dark-blood vessel-wall imaging method using random bipolar gradients with a radial steady-state free precession (SSFP) acquisition in carotid applications. MATERIALS AND METHODS The carotid artery bifurcations of four asymptomatic volunteers (28-37 years old, mean age = 31 years) were included in this study. Dark-blood contrast was achieved through the use of random bipolar gradients applied prior to the signal acquisition of each radial projection in a balanced SSFP acquisition. The resulting phase variation for moving spins established significant destructive interference in the low-frequency region of k-space. This phase variation resulted in a net nulling of the signal from flowing spins, while the bipolar gradients had a minimal effect on the static spins. The net effect was that the regular SSFP signal amplitude (SA) in stationary tissues was preserved while dark-blood contrast was achieved for moving spins. In this implementation, application of the random bipolar gradient pulses along all three spatial directions nulled the signal from both in-plane and through-plane flow in phantom and in vivo studies. RESULTS In vivo imaging trials confirmed that dark-blood contrast can be achieved with the radial random bipolar SSFP method, thereby substantially reversing the vessel-to-lumen contrast-to-noise ratio (CNR) of a conventional rectilinear SSFP "bright-blood" acquisition from bright blood to dark blood with only a modest increase in TR (approximately 4 msec) to accommodate the additional bipolar gradients. CONCLUSION Overall, this sequence offers a simple and effective dark-blood contrast mechanism for high-SNR SSFP acquisitions in vessel wall imaging within a short acquisition time.
Collapse
Affiliation(s)
- Hung-Yu Lin
- Department of Radiology, University Hospitals of Cleveland and Case Western Reserve University, 11100 Euclid Avenue, Cleveland, OH 44195, USA
| | | | | | | |
Collapse
|
17
|
Elgort DR, Hillenbrand CM, Zhang S, Wong EY, Rafie S, Lewin JS, Duerk JL. Image-guided and -monitored renal artery stenting using only MRI. J Magn Reson Imaging 2006; 23:619-27. [PMID: 16555228 DOI: 10.1002/jmri.20554] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022] Open
Abstract
PURPOSE To demonstrate the ability of a unique interventional MR system to be used safely and effectively as the only imaging modality for all phases of MR-guided stent-supported angioplasty. MATERIALS AND METHODS An experimental disease model of renal stenosis was created in six pigs. An interventional MR system, which employed previously reported tools for real-time catheter tracking with automated scan-plane positioning, adaptive image parameters, and radial true-FISP imaging with steady-state precession (True-FISP) imaging coupled with a high-speed reconstruction technique, was then used to guide all phases of the intervention, including: guidewire and catheter insertion, stent deployment, and confirmation of therapeutic success. Pre- and postprocedural X-ray imaging was used as a gold standard to validate the experimental results. RESULTS All of the stent-supported angioplasty interventions were a technical success and were performed without complications. The average postoperative residual stenosis was 14.9%. The image guidance enabled the stents to be deployed with an accuracy of 0.98 +/- 0.69 mm. Additionally, using this interventional MRI system to guide renal artery stenting significantly reduces the procedure time, as compared to using X-ray fluoroscopy. CONCLUSION This study has clearly demonstrated the first successful treatment of renal artery stenting in an experimental animal model solely under MRI guidance and monitoring.
Collapse
Affiliation(s)
- Daniel R Elgort
- Department of Radiology, University Hospitals of Cleveland, Cleveland, Ohio 044106, USA
| | | | | | | | | | | | | |
Collapse
|
18
|
Samsonov AA, Block WF, Arunachalam A, Field AS. Advances in locally constrained k-space-based parallel MRI. Magn Reson Med 2006; 55:431-8. [PMID: 16369917 DOI: 10.1002/mrm.20757] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
In this article, several theoretical and methodological developments regarding k-space-based, locally constrained parallel MRI (pMRI) reconstruction are presented. A connection between Parallel MRI with Adaptive Radius in k-Space (PARS) and GRAPPA methods is demonstrated. The analysis provides a basis for unified treatment of both methods. Additionally, a weighted PARS reconstruction is proposed, which may absorb different weighting strategies for improved image reconstruction. Next, a fast and efficient method for pMRI reconstruction of data sampled on non-Cartesian trajectories is described. In the new technique, the computational burden associated with the numerous matrix inversions in the original PARS method is drastically reduced by limiting direct calculation of reconstruction coefficients to only a few reference points. The rest of the coefficients are found by interpolating between the reference sets, which is possible due to the similar configuration of points participating in reconstruction for highly symmetric trajectories, such as radial and spirals. As a result, the time requirements are drastically reduced, which makes it practical to use pMRI with non-Cartesian trajectories in many applications. The new technique was demonstrated with simulated and actual data sampled on radial trajectories.
Collapse
Affiliation(s)
- Alexey A Samsonov
- Department of Radiology, University of Wisconsin-Madison, Madison, WI 53792-1790, USA.
| | | | | | | |
Collapse
|
19
|
Elgort DR, Duerk JL. A review of technical advances in interventional magnetic resonance imaging. Acad Radiol 2005; 12:1089-99. [PMID: 16099690 DOI: 10.1016/j.acra.2005.06.003] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2005] [Revised: 06/01/2005] [Accepted: 06/01/2005] [Indexed: 10/25/2022]
Abstract
Initial research in the development of interventional magnetic resonance (MR) imaging in the late 1980s and early to mid-1990s focused on pulse sequences, devices, and clinical applications. This focus was largely a result of the limited number of areas in which the academic research community leading the development could provide innovation on the MR systems of the time. However, during the past decade, computational power, higher bandwidth graphical displays, faster computer networks, improved pulse sequence architectures, and improved technical specifications have accelerated the pace of development on modern MR systems. Today, it is the combination of multiple system factors that are enabling the future of interventional MR. These developments, their impact on the field, and newly emerging applications are described.
Collapse
Affiliation(s)
- Daniel R Elgort
- Department of Radiology-MRI, Case Western Reserve University and University Hospitals of Cleveland, 11100 Euclid Avenue, Cleveland, OH 44106, USA
| | | |
Collapse
|
20
|
Larson AC, Kellman P, Arai A, Hirsch GA, McVeigh E, Li D, Simonetti OP. Preliminary investigation of respiratory self-gating for free-breathing segmented cine MRI. Magn Reson Med 2005; 53:159-68. [PMID: 15690515 PMCID: PMC1939886 DOI: 10.1002/mrm.20331] [Citation(s) in RCA: 146] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Segmented cine MRI generally requires breath-holding, which can be problematic for many patients. Navigator echo techniques, particularly successful for free-breathing coronary MRA, are incompatible with the acquisition strategies and SSFP pulse sequences commonly used for cine MRI. The purpose of this work is to introduce a new self-gating technique deriving respiratory gating information directly from the raw imaging data acquired for segmented cine MRI. The respiratory self-gating technique uses interleaved radial k-space sampling to provide low-resolution images in real time during the free-breathing acquisition that are compared to target expiration images. Only the raw data-producing images with high correlation to the target images are included in the final high-resolution reconstruction. The self-gating technique produced cine series with no significant differences in quantitative image sharpness to series produced using comparable breath-held techniques. Because of the difficulties associated with breath-holding, the respiratory self-gating technique represents an important practical advance for cardiac MRI. , Inc.
Collapse
Affiliation(s)
- Andrew C Larson
- Department of Biomedical Engineering, NorthWestern University, Chicago, IL 60611, USA.
| | | | | | | | | | | | | |
Collapse
|
21
|
Dale BM, Lewin JS, Duerk JL. Optimal design of k-space trajectories using a multi-objective genetic algorithm. Magn Reson Med 2005; 52:831-41. [PMID: 15389938 DOI: 10.1002/mrm.20233] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Spiral, radial, and other nonrectilinear k-space trajectories are an area of active research in MRI due largely to their typically rapid acquisition times and benign artifact patterns. Trajectory design has commonly proceeded from a description of a simple shape to an investigation of its properties, because there is no general theory for the derivation of new trajectories with specific properties. Here such a generalized methodology is described. Specifically, a multi-objective genetic algorithm (GA) is used to design trajectories with beneficial flow and off-resonance properties. The algorithm converges to a well-defined optimal set with standard spiral trajectories on the rapid but low-quality end, and a new class of trajectories on the slower but high-quality end. The new trajectories all begin with nonzero gradient amplitude at the k-space origin, and curve gently outward relative to standard spirals. Improvements predicted in simulated imaging experiments were found to correlate well with improvements in actual experimental measures of image quality. The impact of deviations from the desired k-space trajectory is described, as is the impact of using different phantoms.
Collapse
Affiliation(s)
- Brian M Dale
- Department of Biomedical Engineering, Case Western Reserve University, Cleveland, Ohio 44106, USA
| | | | | |
Collapse
|
22
|
Peters DC, Guttman MA, Dick AJ, Raman VK, Lederman RJ, McVeigh ER. Reduced field of view and undersampled PR combined for interventional imaging of a fully dynamic field of view. Magn Reson Med 2004; 51:761-7. [PMID: 15065249 PMCID: PMC2034278 DOI: 10.1002/mrm.20037] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
Active catheter imaging was investigated using real-time undersampled projection reconstruction (PR) combined with the temporal filtering technique of reduced field of view (rFOV). Real-time rFOV processing was interactively enabled during highly undersampled catheter imaging, resulting in improved artifact suppression with better temporal resolution than that obtained by view-sharing. Imaging with 64 to 32 projections provided a resolution of 2 x 2 x 8 mm, and four to eight true frames per second. Image artifacts were reduced when rFOV processing was applied to the undersampled images. A comparison with Cartesian rFOV showed that PR image quality is less susceptible to aliasing that results from rFOV imaging with a wholly dynamic outer FOV. Simulations and MRI experiments demonstrated that PR rFOV provides significant artifact suppression, even for a fully dynamic FOV. The near doubling of temporal resolution that is possible with PR rFOV permits accurate monitoring of highly dynamic events, such as catheter movements, and arrhythmias, such as ventricular ectopy.
Collapse
Affiliation(s)
- Dana C Peters
- Laboratory of Cardiac Energetics, NHLBI, National Institutes of Health, DHHS, Bethesda, Maryland, USA.
| | | | | | | | | | | |
Collapse
|
23
|
Abstract
The need for ECG gating presents many difficulties in cardiac magnetic resonance imaging (CMRI). Real-time imaging techniques eliminate the need for ECG gating in cine CMRI, but they cannot offer the spatial and temporal resolution provided by segmented acquisition techniques. Previous MR signal-based techniques have demonstrated an ability to provide cardiac gating information; however, these techniques result in decreased imaging efficiency. The purpose of this work was to develop a new "self-gated" (SG) acquisition technique that eliminates these efficiency deficits by extracting the motion synchronization signal directly from the same MR signals used for image reconstruction. Three separate strategies are proposed for deriving the SG signal from data acquired using radial k-space sampling: echo peak magnitude, kymogram, and 2D correlation. The SG techniques were performed on seven normal volunteers. A comparison of the results showed that they provided cine image series with no significant differences in image quality compared to that obtained with conventional ECG gating techniques. SG techniques represent an important practical advance in clinical MRI because they enable the acquisition of high temporal and spatial resolution cardiac cine images without the need for ECG gating and with no loss in imaging efficiency.
Collapse
Affiliation(s)
- Andrew C Larson
- Laboratory of Cardiac Energetics, NHLBI, National Institutes of Health, Department of Health and Human Services, Bethesda, Maryland 20892-1061, USA.
| | | | | | | | | | | |
Collapse
|
24
|
Shankaranarayanan A, Herfkens R, Hargreaves BM, Polzin JA, Santos JM, Brittain JH. Helical MR: Continuously moving table axial imaging with radial acquisitions. Magn Reson Med 2003; 50:1053-60. [PMID: 14587016 DOI: 10.1002/mrm.10621] [Citation(s) in RCA: 31] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
A technique for extended field of view MRI is presented. Similar to helical computed tomography, the method utilizes a continuously moving patient table, a 2D axial slice that remains fixed relative to the MRI magnet, and a radial k-space trajectory. A fully refocused SSFP acquisition enables spatial resolution comparable to current clinical protocols in scan times that are sufficiently short to allow a reasonable breathhold duration. RF transmission and signal reception are performed using the RF body coil and the images are reconstructed in real time. Experimental results are presented that illustrate the technique's ability to resolve small structures in the table-motion direction. Simulation experiments to study the steady-state response of the fully refocused SSFP acquisition during continuous table motion are also presented. Finally, whole body images of healthy volunteers demonstrate the high image quality achieved using the helical MRI approach.
Collapse
|
25
|
Flask CA, Dale B, Lewin JS, Duerk JL. Radial alternating TE sequence for faster fat suppression. Magn Reson Med 2003; 50:1095-9. [PMID: 14587021 DOI: 10.1002/mrm.10615] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
This study describes a steady-state sequence that uses a radial k-space trajectory and alternating echo times (TEs) between even and odd k-space views. The sequence generated a single data set that was used to reconstruct images with inherent fat suppression. This fat suppression results from the fat phase variation in alternate echoes giving rise to cancellation in the central portion of k-space. This new fat-suppression method provides inherent fat suppression in half the acquisition time relative to the radial two-point Dixon method. The improvement in k-space sampling efficiency is demonstrated in phantom and clinical images, and through measured point-spread functions (PSFs). As a result, the radial alternating TE sequence offers improved temporal resolution over a radial version of the two-point Dixon sequence by requiring fewer total projections to obtain the same effective resolution in water-based tissues.
Collapse
Affiliation(s)
- Chris A Flask
- Department of Biomedical Engineering, Case Western Reserve University and University Hospitals of Cleveland, Cleveland, Ohio 44106, USA
| | | | | | | |
Collapse
|
26
|
Qian Y, Lin J, Jin D. Reconstruction of MR images from data acquired on an arbitrary k-space trajectory using the same-image weight. Magn Reson Med 2002; 48:306-11. [PMID: 12210939 DOI: 10.1002/mrm.10218] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
A sampling density compensation function denoted "same-image (SI) weight" is proposed to reconstruct MR images from the data acquired on an arbitrary k-space trajectory. An equation for the SI weight is established on the SI criterion and an iterative scheme is developed to find the weight. The SI weight is then used to reconstruct images from the data calculated on a random trajectory in a numerical phantom case and from the data acquired on interleaved spirals in an in vivo experiment, respectively. In addition, Pipe and Menon's weight (MRM 1999;41:179-186) is also used in the reconstructions to make a comparison. The images obtained with the SI weight were found to be slightly more accurate than those obtained with Pipe's weight.
Collapse
Affiliation(s)
- Yongxian Qian
- Institute of Biomedical Engineering, Huazhong University of Science & Technology, Wuhan, Hubei, China.
| | | | | |
Collapse
|
27
|
Qian Y, Lin J, Jin D. Direct reconstruction of MR images from data acquired on a non-Cartesian grid using an equal-phase-line algorithm. Magn Reson Med 2002; 47:1228-33. [PMID: 12111970 DOI: 10.1002/mrm.10165] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The equal-phase line (EPL) algorithm is proposed as a means of allowing rapid Fourier transform (FT) reconstruction of MR image data acquired on a non-Cartesian grid. The pixels on the image are grouped according to their positions. The pixels in a group have the same phase in the complex exponential function -exp[j2pi(xu + yv)] and receive the same contribution from a data point. Each group is related to an EPL in the image space. The contribution of a data point can then be distributed to the pixels along the EPLs. The described EPL algorithm enables a decrease of the reconstruction time to about 40% of the direct FT (DrFT) for the non-Cartesian data. A numerical phantom and two sets of in vivo spiral data were used to investigate an optimal number of the EPLs and to measure the reconstruction time. The EPL algorithm runs nearly as fast as the look-up table (LUT) method (Dale et al. IEEE Trans Med Imaging 2001;20:207-217), but it does not require a large memory to store the coefficients in advance, as is required in the LUT method. Thus, the EPL algorithm can be used to reconstruct images up to 512 x 512 pixels in size in a PC of limited memory, and may be more conveniently applied to a multiprocessor system.
Collapse
Affiliation(s)
- Yongxian Qian
- MRI Research Group, Institute of Biomedical Engineering, Huazhong University of Science and Technology, Wuhan, Hubei, China.
| | | | | |
Collapse
|