1
|
Geldschläger O, Bosch D, Henning A. OTUP workflow: target specific optimization of the transmit k-space trajectory for flexible universal parallel transmit RF pulse design. NMR IN BIOMEDICINE 2022; 35:e4728. [PMID: 35297104 DOI: 10.1002/nbm.4728] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/26/2021] [Revised: 02/09/2022] [Accepted: 03/14/2022] [Indexed: 06/14/2023]
Abstract
PURPOSE To optimize transmit k-space trajectories for a wide range of excitation targets and to design "universal pTx RF pulses" based on these trajectories. METHODS Transmit k-space trajectories (stack of spirals and SPINS) were optimized to best match different excitation targets using the parameters of the analytical equations of spirals and SPINS. The performances of RF pulses designed based on optimized and non-optimized trajectories were compared. The optimized trajectories were utilized for universal pulse design. The universal pulse performances were compared with subject specific tailored pulse performances. The OTUP workflow (optimization of transmit k-space trajectories and universal pulse calculation) was tested on three test target excitation patterns. For one target (local excitation of a central area in the human brain) the pulses were tested in vivo at 9.4 T. RESULTS The workflow produced appropriate transmit k-space trajectories for each test target. Utilization of an optimized trajectory was crucial for the pulse performance. Using unsuited trajectories diminished the performance. It was possible to create target specific universal pulses. However, not every test target is equally well suited for universal pulse design. There was no significant difference in the in vivo performance between subject specific tailored pulses and a universal pulse at 9.4 T. CONCLUSIONS The proposed workflow further exploited and improved the universal pulse concept by combining it with gradient trajectory optimization for stack of spirals and SPINS. It emphasized the importance of a well suited trajectory for pTx RF pulse design. Universal and tailored pulses performed with a sufficient degree of similarity in simulations and a high degree of similarity in vivo. The implemented OTUP workflow and the B0 /B1+ map data from 18 subjects measured at 9.4 T are available as open source (https://github.com/ole1965/workflow_OTUP.git).
Collapse
Affiliation(s)
- Ole Geldschläger
- High-Field Magnetic Resonance Center, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
| | - Dario Bosch
- High-Field Magnetic Resonance Center, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
- Biomedical Magnetic Resonance, University Hospital Tübingen, Tübingen, Germany
| | - Anke Henning
- High-Field Magnetic Resonance Center, Max Planck Institute for Biological Cybernetics, Tübingen, Germany
- Advanced Imaging Research Center, University of Texas Southwestern Medical Center, Dallas, TX, USA
| |
Collapse
|
2
|
Engel M, Kasper L, Wilm B, Dietrich B, Vionnet L, Hennel F, Reber J, Pruessmann KP. T-Hex: Tilted hexagonal grids for rapid 3D imaging. Magn Reson Med 2020; 85:2507-2523. [PMID: 33270941 DOI: 10.1002/mrm.28600] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2020] [Revised: 10/20/2020] [Accepted: 10/22/2020] [Indexed: 01/09/2023]
Abstract
PURPOSE The purpose of this work is to devise and demonstrate an encoding strategy for 3D MRI that reconciles high speed with flexible segmentation, uniform k-space density, and benign T 2 ∗ effects. METHODS Fast sampling of a 3D k-space is typically accomplished by 2D readouts per shot using EPI trains or spiral readouts. Tilted hexagonal (T-Hex) sampling is a way of acquiring more k-space volume per excitation while maintaining uniform sampling density and a smooth T 2 ∗ filter. The k-space volume covered per shot is controlled by the tilting angle. Image reconstruction is performed with a 3D extension of the iterative SENSE approach, incorporating actual field dynamics and static off-resonance. T-Hex imaging is compared with established 3D schemes in terms of speed and noise performance. RESULTS Tilted hexagonal acquisition is found to achieve greater imaging speed than known alternatives, particularly in combination with spiral trajectories. The interplay of the proposed 3D trajectories, array detection, and off-resonance is successfully addressed by iterative inversion of the full signal model. Enhanced coverage per shot is of greatest utility for high speed in an intermediate resolution regime of 1 to 4 mm. T-Hex EPI combines the benefits of extended coverage per shot with increased robustness against off-resonance effects. CONCLUSION Sampling of tilted hexagonal grids is a feasible means of gaining 3D imaging speed with near-optimal SNR efficiency and benign depiction properties. It is a particularly promising technique for time-resolved applications such as fMRI.
Collapse
Affiliation(s)
- Maria Engel
- Institute for Biomedical Engineering, ETH Zurich and University of Zurich, Zurich, Switzerland
| | - Lars Kasper
- Institute for Biomedical Engineering, ETH Zurich and University of Zurich, Zurich, Switzerland.,Translational Neuromodeling Unit, IBT, University of Zurich and ETH Zurich, Zurich, Switzerland
| | - Bertram Wilm
- Institute for Biomedical Engineering, ETH Zurich and University of Zurich, Zurich, Switzerland
| | - Benjamin Dietrich
- Institute for Biomedical Engineering, ETH Zurich and University of Zurich, Zurich, Switzerland
| | - Laetitia Vionnet
- Institute for Biomedical Engineering, ETH Zurich and University of Zurich, Zurich, Switzerland
| | - Franciszek Hennel
- Institute for Biomedical Engineering, ETH Zurich and University of Zurich, Zurich, Switzerland
| | - Jonas Reber
- Institute for Biomedical Engineering, ETH Zurich and University of Zurich, Zurich, Switzerland
| | - Klaas P Pruessmann
- Institute for Biomedical Engineering, ETH Zurich and University of Zurich, Zurich, Switzerland
| |
Collapse
|
3
|
Sharma S, Coutino M, Chepuri SP, Leus G, Hari KVS. Towards a general framework for fast and feasible k-space trajectories for MRI based on projection methods. Magn Reson Imaging 2020; 72:122-134. [PMID: 32668272 DOI: 10.1016/j.mri.2020.06.016] [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: 02/20/2020] [Revised: 06/11/2020] [Accepted: 06/24/2020] [Indexed: 10/23/2022]
Abstract
The design of feasible trajectories to traverse the k-space for sampling in magnetic resonance imaging (MRI) is important while considering ways to reduce the scan time. Over the recent years, non-Cartesian trajectories have been observed to result in benign artifacts and being less sensitive to motion. In this paper, we propose a generalized framework that encompasses projection-based methods to generate feasible non-Cartesian k-space trajectories. This framework allows to construct feasible trajectories from both random or structured initial trajectories, e.g., based on the traveling salesman problem (TSP). We evaluate the performance of the proposed methods by simulating the reconstruction of 128 × 128 and 256 × 256 phantom and brain MRI images in terms of structural similarity (SSIM) index and peak signal-to-noise ratio (PSNR) using compressed sensing techniques. It is observed that the TSP-based trajectories from the proposed projection method with constant acceleration parameterization (CAP) result in better reconstruction compared to the projection method with constant velocity parameterization (CVP) and this for a similar read-out time. Further, random-like trajectories are observed to be better than TSP-based trajectories as they reduce the read-out time while providing better reconstruction quality. A reduction in read-out time by upto 67% is achieved using the proposed projection with permutation (PP) method.
Collapse
Affiliation(s)
- Shubham Sharma
- Department of Electrical Communication Engineering, Indian Institute of Science, Bangalore, India.
| | - Mario Coutino
- Department of Microelectronics, Delft University of Technology, Netherlands
| | | | - Geert Leus
- Department of Microelectronics, Delft University of Technology, Netherlands
| | - K V S Hari
- Department of Electrical Communication Engineering, Indian Institute of Science, Bangalore, India
| |
Collapse
|
4
|
Wang F, Hennig J, LeVan P. Time-domain principal component reconstruction (tPCR): A more efficient and stable iterative reconstruction framework for non-Cartesian functional MRI. Magn Reson Med 2020; 84:1321-1335. [PMID: 32068309 DOI: 10.1002/mrm.28208] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2019] [Revised: 12/27/2019] [Accepted: 01/19/2020] [Indexed: 12/20/2022]
Abstract
PURPOSE To improve the reconstruction efficiency (i.e., computational load) and stability of iterative reconstruction for non-Cartesian fMRI when using high undersampling rates and/or in the presence of strong off-resonance effects. THEORY AND METHODS The magnetic resonance encephalography (MREG) sequence with 3D non-Cartesian trajectory and 0.1s repetition time (TR) was applied to acquire fMRI datasets. Different from a conventional time-point-by-time-point sequential reconstruction (SR), the proposed time-domain principal component reconstruction (tPCR) performs three steps: (1) decomposing the k-t-space fMRI datasets into time-domain principal component space using singular value decomposition, (2) reconstructing each principal component with redistributed computation power according to their weights, and (3) combining the reconstructed principal components back to image-t-space. The comparison of reconstruction accuracy was performed by simulation experiments and then verified in real fMRI data. RESULTS The simulation experiments showed that the proposed tPCR was able to significantly reduce reconstruction errors, and subsequent functional activation errors, relative to SR at identical computational cost. Alternatively, at fixed reconstruction accuracy, computation time was greatly reduced. The improved performance was particularly obvious for L1-norm nonlinear reconstructions relative to L2-norm linear reconstructions and robust to different regularization strength, undersampling rates, and off-resonance effects intensity. By examining activation maps, tPCR was also found to give similar improvements in real fMRI experiments. CONCLUSION The proposed proof-of-concept tPCR framework could improve (1) the reconstruction efficiency of iterative reconstruction, and (2) the reconstruction stability especially for nonlinear reconstructions. As a practical consideration, the improved reconstruction speed promotes the application of highly undersampled non-Cartesian fast fMRI.
Collapse
Affiliation(s)
- Fei Wang
- Department of Radiology, Medical Physics, Faculty of Medicine, Medical Center - University of Freiburg, Freiburg, Germany.,Center for Basics in NeuroModulation (NeuroModul Basics), Faculty of Medicine, University of Freiburg, Freiburg, Germany
| | - Jürgen Hennig
- Department of Radiology, Medical Physics, Faculty of Medicine, Medical Center - University of Freiburg, Freiburg, Germany.,Center for Basics in NeuroModulation (NeuroModul Basics), Faculty of Medicine, University of Freiburg, Freiburg, Germany
| | - Pierre LeVan
- Department of Radiology, Medical Physics, Faculty of Medicine, Medical Center - University of Freiburg, Freiburg, Germany.,Departments of Radiology and Paediatrics, Cumming School of Medicine, University of Calgary, Calgary, Canada.,Alberta Children's Hospital Research Institute and Hotchkiss Brain Institute, University of Calgary, Calgary, Canada
| |
Collapse
|
5
|
Li Q, Liao C, Ye H, Chen Y, Cao X, Yuan L, He H, Zhong J. Squeezed Trajectory Design for Peak RF and Integrated RF Power Reduction in Parallel Transmission MRI. IEEE TRANSACTIONS ON MEDICAL IMAGING 2018; 37:1809-1821. [PMID: 29993630 DOI: 10.1109/tmi.2018.2828112] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
High peak RF amplitude and excessive specific absorption rate (SAR) are two critical concerns for hardware implementation and patient safety in scientific and clinical research for high field MRI using parallel transmissions (pTX). In this paper, we introduce a squeezing strategy to reduce peak RF amplitude and integrated RF power via direct reshaping of the k-space trajectory. In the existing peak RF / integrated RF power optimization methods gradient amplitude or slew rate is reduced, but the k-space trajectory remains unchanged. Unlike these traditional methods, we worked directly in the excitation k-space to reshape k-space traversal by a squeezing vector in order to achieve peak RF and total RF power optimization, using a particle swarm optimization algorithm. The squeezing strategy was applied to the conventional variable density spiral (CVDS) and the variable rate selective excitation (VERSE) trajectories, dubbed SVDS (squeezed variable density spiral) and SVERSE (squeezing trajectory with VERSE), respectively, for different excitation profiles of small or large tip angles. Pulse acceleration and off-resonance effects were evaluated for an 8-ch pTX via Bloch simulation. CVDS, VERSE, SVDS, and SVERSE pulses were implemented on a 3T scanner with a 2-ch pTX. Phantom and in vivo experiments were performed for reduced FOV (rFOV) imaging. The results show that SVDS pulses simultaneously reduce integrated RF power and peak RF by about 30% on average compared to CVDS pulses for a square pattern ( $80\times80$ mm2) with flip angles of 30°, 90°, and 180°. Compared with the VERSE method under the same peak RF constraints, the SVDS method reduces integrated RF power by an average of 20% for small tip excitations for profiles of slice, rectangular, square, and circle, and has slightly reduced excitation accuracy slightly (about 0.6%, from 6.8% to 7.4%). The SVERSE method shortens the duration of the VERSE pulse by 12.8% at large ti p angle (180°). Feasibility for rFOV imaging was demonstrated with phantom and in vivo experiments with squeezed pulses.
Collapse
|
6
|
Davids M, Schad LR, Wald LL, Guérin B. Fast three-dimensional inner volume excitations using parallel transmission and optimized k-space trajectories. Magn Reson Med 2016; 76:1170-82. [PMID: 26527590 PMCID: PMC4854802 DOI: 10.1002/mrm.26021] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2015] [Revised: 09/28/2015] [Accepted: 09/29/2015] [Indexed: 11/05/2022]
Abstract
PURPOSE To design short parallel transmission (pTx) pulses for excitation of arbitrary three-dimensional (3D) magnetization patterns. METHODS We propose a joint optimization of the pTx radiofrequency (RF) and gradient waveforms for excitation of arbitrary 3D magnetization patterns. Our optimization of the gradient waveforms is based on the parameterization of k-space trajectories (3D shells, stack-of-spirals, and cross) using a small number of shape parameters that are well-suited for optimization. The resulting trajectories are smooth and sample k-space efficiently with few turns while using the gradient system at maximum performance. Within each iteration of the k-space trajectory optimization, we solve a small tip angle least-squares RF pulse design problem. Our RF pulse optimization framework was evaluated both in Bloch simulations and experiments on a 7T scanner with eight transmit channels. RESULTS Using an optimized 3D cross (shells) trajectory, we were able to excite a cube shape (brain shape) with 3.4% (6.2%) normalized root-mean-square error in less than 5 ms using eight pTx channels and a clinical gradient system (Gmax = 40 mT/m, Smax = 150 T/m/s). This compared with 4.7% (41.2%) error for the unoptimized 3D cross (shells) trajectory. Incorporation of B0 robustness in the pulse design significantly altered the k-space trajectory solutions. CONCLUSION Our joint gradient and RF optimization approach yields excellent excitation of 3D cube and brain shapes in less than 5 ms, which can be used for reduced field of view imaging and fat suppression in spectroscopy by excitation of the brain only. Magn Reson Med 76:1170-1182, 2016. © 2015 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Mathias Davids
- Computer Assisted Clinical Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim, BW, Germany.
- A. A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, United States.
| | - Lothar R Schad
- Computer Assisted Clinical Medicine, Medical Faculty Mannheim, Heidelberg University, Mannheim, BW, Germany
| | - Lawrence L Wald
- A. A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, United States
- Harvard-MIT, Division of Health Sciences and Technology, Cambridge, Massachusetts, United States
- Harvard Medical School, Boston, Massachusetts, United States
| | - Bastien Guérin
- A. A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, Massachusetts, United States
- Harvard Medical School, Boston, Massachusetts, United States
| |
Collapse
|
7
|
Chauffert N, Weiss P, Kahn J, Ciuciu P. A Projection Algorithm for Gradient Waveforms Design in Magnetic Resonance Imaging. IEEE TRANSACTIONS ON MEDICAL IMAGING 2016; 35:2026-2039. [PMID: 27019479 DOI: 10.1109/tmi.2016.2544251] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Collecting the maximal amount of information in a given scanning time is a major concern in magnetic resonance imaging (MRI) to speed up image acquisition. The hardware constraints (gradient magnitude, slew rate, etc.), physical distortions (e.g., off-resonance effects) and sampling theorems (Shannon, compressed sensing) must be taken into account simultaneously, which makes this problem extremely challenging. To date, the main approach to design gradient waveform has consisted of selecting an initial shape (e.g., spiral, radial lines, etc.) and then traversing it as fast as possible using optimal control. In this paper, we propose an alternative solution which first consists of defining a desired parameterization of the trajectory and then of optimizing for minimal deviation of the sampling points within gradient constraints. This method has various advantages. First, it better preserves the density of the input curve which is critical in sampling theory. Second, it allows to smooth high curvature areas making the acquisition time shorter in some cases. Third, it can be used both in the Shannon and CS sampling theories. Last, the optimized trajectory is computed as the solution of an efficient iterative algorithm based on convex programming. For piecewise linear trajectories, as compared to optimal control reparameterization, our approach generates a gain in scanning time of 10% in echo planar imaging while improving image quality in terms of signal-to-noise ratio (SNR) by more than 6 dB. We also investigate original trajectories relying on traveling salesman problem solutions. In this context, the sampling patterns obtained using the proposed projection algorithm are shown to provide significantly better reconstructions (more than 6 dB) while lasting the same scanning time.
Collapse
|
8
|
Sun H, Fessler JA, Noll DC, Nielsen JF. Joint Design of Excitation k-Space Trajectory and RF Pulse for Small-Tip 3D Tailored Excitation in MRI. IEEE TRANSACTIONS ON MEDICAL IMAGING 2016; 35:468-79. [PMID: 26390450 PMCID: PMC4792784 DOI: 10.1109/tmi.2015.2478880] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
We propose a new method for the joint design of k-space trajectory and RF pulse in 3D small-tip tailored excitation. Designing time-varying RF and gradient waveforms for a desired 3D target excitation pattern in MRI poses a non-linear, non-convex, constrained optimization problem with relatively large problem size that is difficult to solve directly. Existing joint pulse design approaches are therefore typically restricted to predefined trajectory types such as EPI or stack-of-spirals that intrinsically satisfy the gradient maximum and slew rate constraints and reduce the problem size (dimensionality) dramatically, but lead to suboptimal excitation accuracy for a given pulse duration. Here we use a 2nd-order B-spline basis that can be fitted to an arbitrary k-space trajectory, and allows the gradient constraints to be implemented efficiently. We show that this allows the joint optimization problem to be solved with quite general k-space trajectories. Starting from an arbitrary initial trajectory, we first approximate the trajectory using B-spline basis, and then optimize the corresponding coefficients. We evaluate our method in simulation using four different k-space initializations: stack-of-spirals, SPINS, KT-points, and a new method based on KT-points. In all cases, our approach leads to substantial improvement in excitation accuracy for a given pulse duration. We also validated our method for inner-volume excitation using phantom experiments. The computation is fast enough for online applications.
Collapse
Affiliation(s)
- Hao Sun
- Department of Electrical Engineering and Computer Science, University of Michigan, Ann Arbor, MI 48109 USA
| | - Jeffrey A. Fessler
- Department of Electrical Engineering and Computer Science, University of Michigan, Ann Arbor, MI 48109 USA
| | - Douglas C. Noll
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI 48109 USA
| | - Jon-Fredrik Nielsen
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI 48109 USA
| |
Collapse
|