Performance evaluation of matrix gradient coils.

OBJECTIVE: In this paper, we present a new performance measure of a matrix coil (also known as multi-coil) from the perspective of efficient, local, non-linear encoding without explicitly considering target encoding fields. MATERIALS AND METHODS: An optimization problem based on a joint optimization for the non-linear encoding fields is formulated. Based on the derived objective function, a figure of merit of a matrix coil is defined, which is a generalization of a previously known resistive figure of merit for traditional gradient coils. RESULTS: A cylindrical matrix coil design with a high number of elements is used to illustrate the proposed performance measure. The results are analyzed to reveal novel features of matrix coil designs, which allowed us to optimize coil parameters, such as number of coil elements. A comparison to a scaled, existing multi-coil is also provided to demonstrate the use of the proposed performance parameter. CONCLUSIONS: The assessment of a matrix gradient coil profits from using a single performance parameter that takes the local encoding performance of the coil into account in relation to the dissipated power.

Monoplanar gradient system for imaging with nonlinear gradients.

OBJECT: In this paper we present a monoplanar gradient system capable of imaging a volume comparable with that covered by linear gradient systems. Such a system has been designed and implemented. MATERIALS AND METHODS: Building such a system was made possible by relaxing the constraint of global linearity and replacing it with a requirement for local orthogonality. A framework was derived for optimization of local orthogonality within the physical boundaries and geometric constraints. Spatial encoding of magnetic fields was optimized for their local orthogonality over a large field of view. RESULTS: A coil design consisting of straight wire segments was optimized, implemented, and integrated into a 3T human scanner to show the feasibility of this approach. Initial MR images are shown and further applications of the derived optimization method and the nonlinear planar gradient system are discussed. CONCLUSION: Encoding fields generated by the prototype encoding system were shown to be locally orthogonal and able to encode a cylindrical volume sufficient for some abdomen imaging applications for humans.

A g-factor metric for k-t-GRAPPA- and PEAK-GRAPPA-based parallel imaging.

PURPOSE: The aim of this work is to derive a theoretical framework for quantitative noise and temporal fidelity analysis of time-resolved k-space-based parallel imaging methods. THEORY: An analytical formalism of noise distribution is derived extending the existing g-factor formulation for nontime-resolved generalized autocalibrating partially parallel acquisition (GRAPPA) to time-resolved k-space-based methods. The noise analysis considers temporal noise correlations and is further accompanied by a temporal filtering analysis. METHODS: All methods are derived and presented for k-t-GRAPPA and PEAK-GRAPPA. A sliding window reconstruction and nontime-resolved GRAPPA are taken as a reference. Statistical validation is based on series of pseudoreplica images. The analysis is demonstrated on a short-axis cardiac CINE dataset. RESULTS: The superior signal-to-noise performance of time-resolved over nontime-resolved parallel imaging methods at the expense of temporal frequency filtering is analytically confirmed. Further, different temporal frequency filter characteristics of k-t-GRAPPA, PEAK-GRAPPA, and sliding window are revealed. CONCLUSION: The proposed analysis of noise behavior and temporal fidelity establishes a theoretical basis for a quantitative evaluation of time-resolved reconstruction methods. Therefore, the presented theory allows for comparison between time-resolved parallel imaging methods and also nontime-resolved methods. Magn Reson Med 74:125-135, 2015. © 2014 Wiley Periodicals, Inc.

Parallel imaging with phase scrambling.

PURPOSE: Most existing methods for accelerated parallel imaging in MRI require additional data, which are used to derive information about the sensitivity profile of each radiofrequency (RF) channel. In this work, a method is presented to avoid the acquisition of separate coil calibration data for accelerated Cartesian trajectories. METHODS: Quadratic phase is imparted to the image to spread the signals in k-space (aka phase scrambling). By rewriting the Fourier transform as a convolution operation, a window can be introduced to the convolved chirp function, allowing a low-resolution image to be reconstructed from phase-scrambled data without prominent aliasing. This image (for each RF channel) can be used to derive coil sensitivities to drive existing parallel imaging techniques. As a proof of concept, the quadratic phase was applied by introducing an offset to the x(2) - y(2) shim and the data were reconstructed using adapted versions of the image space-based sensitivity encoding and GeneRalized Autocalibrating Partially Parallel Acquisitions algorithms. RESULTS: The method is demonstrated in a phantom (1 × 2, 1 × 3, and 2 × 2 acceleration) and in vivo (2 × 2 acceleration) using a 3D gradient echo acquisition. CONCLUSION: Phase scrambling can be used to perform parallel imaging acceleration without acquisition of separate coil calibration data, demonstrated here for a 3D-Cartesian trajectory. Further research is required to prove the applicability to other 2D and 3D sampling schemes.

Image reconstruction in k-space from MR data encoded with ambiguous gradient fields.

PURPOSE: In this work, the limits of image reconstruction in k-space are explored when non-bijective gradient fields are used for spatial encoding. THEORY: The image space analogy between parallel imaging and imaging with non-bijective encoding fields is partially broken in k-space. As a consequence, it is hypothesized and proven that ambiguities can only be resolved partially in k-space, and not completely as is the case in image space. METHODS: Image-space and k-space based reconstruction algorithms for multi-channel radiofrequency data acquisitions are programmed and tested using numerical simulations as well as in vivo measurement data. RESULTS: The hypothesis is verified based on an analysis of reconstructed images. It is found that non-bijective gradient fields have the effect that densely sampled autocalibration data, used for k-space reconstruction, provide less information than a separate scan of the receiver coil sensitivity maps, used for image space reconstruction. Consequently, in k-space only the undersampling artifact can be unfolded, whereas in image space, it is also possible to resolve aliasing that is caused by the non-bijectivity of the gradient fields. CONCLUSION: For standard imaging, reconstruction in image space and in k-space is nearly equivalent, whereas there is a fundamental difference with practical consequences for the selection of image reconstruction algorithms when non-bijective encoding fields are involved.

MR image reconstruction from generalized projections.

PURPOSE: Currently, the time required for image reconstruction is prohibitively long if data are acquired using multidimensional imaging trajectories that make use of multichannel systems equipped with nonlinear gradients. Methods are presented that reduce the computational complexity of the iterative time-domain reconstruction algorithm down from O(N(4)) to O(N(3)). THEORY: For generalized projections, a large class of multidimensional imaging trajectories, the encoding matrix can be focused to sparse bands by introducing an appropriate filter function along the frequency-encoding direction. The reconstruction can be speeded up by ignoring values below a predefined threshold level. METHODS: Two methods are presented that differ in how the filter is incorporated into the reconstruction algorithm. The first method represents, without implementation of a threshold, a weighted version of the time-domain method, while the second method is equivalent to it. RESULTS: Simulation and measurement results show that image reconstruction from high-resolution imaging data can be speeded up by up to two orders of magnitude. While the weighted reconstruction requires more iterations to reach an optimum than the second method, it is less sensitive to thresholding. CONCLUSION: For complex spatial encoding strategies that involve nonlinear gradient fields, fast and accurate image reconstruction methods are provided that are particularly efficient for high-resolution anatomical imaging.

Stages: sub-Fourier dynamic shim updating using nonlinear magnetic field phase preparation.

Heterogeneity of the static magnetic field in magnetic resonance imaging may cause image artifacts and degradation in image quality. The field heterogeneity can be reduced by dynamically adjusting shim fields or dynamic shim updating, in which magnetic field homogeneity is optimized for each tomographic slice to improve image quality. A limitation of this approach is that a new magnetic field can be applied only once for each slice, otherwise image quality would improve somewhere to its detriment elsewhere in the slice. The motivation of this work is to overcome this limitation and develop a technique using nonlinear magnetic fields to dynamically shim the static magnetic field within a single Fourier-encoded volume or slice, called sub-Fourier dynamic shim updating. However, the nonlinear magnetic fields are not used as shim fields; instead, they impart a strong spatial dependence to the acquired MR signal by nonlinear phase preparation, which may be exploited to locally improve magnetic field homogeneity during acquisition. A theoretical description of the method is detailed, simulations and a proof-of-principle experiment are performed using a magnet coil with a known field geometry. The method is shown to remove artifacts associated with magnetic field homogeneity in balanced steady-state free-precession pulse sequences. We anticipate that this method will be useful to improve the quality of magnetic resonance images by removing deleterious artifacts associated with a heterogeneous static magnetic field.

Local field of view imaging for alias-free undersampling with nonlinear spatial encoding magnetic fields.

PURPOSE: Nonlinear spatial encoding magnetic fields result in an inhomogeneous image resolution. Within this study, this characteristic property of nonlinear encoding is investigated with regard to its potential to accelerate MRI acquisitions. THEORY: A dependency between k-space coverage and local resolvability of the image causes k-space samples to have a spatially localized contribution to the reconstruction of the spin density. On the basis of this observation, a concept for alias-free data undersampling is developed, which is referred to as the local field of view concept. METHODS: On the basis of this concept, a fast sampling trajectory is developed. It is evaluated with simulations and experiments (both using a phantom and in vivo) for MRI with, as an example, pure quadrupolar encoding fields. To demonstrate that the concept is only applicable to (spatially) nonlinear encoding, a comparison with linear encoding is provided. RESULTS: Application of the local field of view concept results in a localized adaptation of the image resolution by undersampling higher frequency k-space samples without introducing aliasing. CONCLUSIONS: A new effect of nonlinear spatial encoding magnetic fields was found, which allows more efficient data sampling and at the same time counterbalancing the natural variation in image resolution.

Excitation and geometrically matched local encoding of curved slices.

In this work, the concept of excitation and geometrically matched local in-plane encoding of curved slices (ExLoc) is introduced. ExLoc is based on a set of locally near-orthogonal spatial encoding magnetic fields, thus maintaining a local rectangular shape of the individual voxels and avoiding potential problems arising due to highly irregular voxel shapes. Unlike existing methods for exciting curved slices based on multidimensional radiofrequency-pulses, excitation and geometrically matched local encoding of curved slices does not require long duration or computationally expensive radiofrequency-pulses. As each encoding field consists of a superposition of potentially arbitrary (spatially linear or nonlinear) magnetic field components, the resulting field shape can be adapted with high flexibility to the specific region of interest. For extended nonplanar structures, this results in improved relevant volume coverage for fewer excited slices and thus increased efficiency. In addition to the mathematical description for the generation of dedicated encoding fields and data reconstruction, a verification of the ExLoc concept in phantom experiments and examples for in vivo curved single and multislice imaging are presented.

Reconstruction of undersampled radial PatLoc imaging using total generalized variation.

In the case of radial imaging with nonlinear spatial encoding fields, a prominent star-shaped artifact has been observed if a spin distribution is encoded with an undersampled trajectory. This work presents a new iterative reconstruction method based on the total generalized variation, which reduces this artifact. For this approach, a sampling operator (as well as its adjoint) is needed that maps data from PatLoc k-space to the final image space. It is shown that this can be realized as a type-3 nonuniform fast Fourier transform, which is implemented by a combination of a type-1 and type-2 nonuniform fast Fourier transform. Using this operator, it is also possible to implement an iterative conjugate gradient SENSE based method for PatLoc reconstruction, which leads to a significant reduction of computation time in comparison to conventional PatLoc image reconstruction methods. Results from numerical simulations and in vivo PatLoc measurements with as few as 16 radial projections are presented, which demonstrate significant improvements in image quality with the total generalized variation-based approach.

Practical considerations for in vivo MRI with higher dimensional spatial encoding.

OBJECT: This work seeks to examine practical aspects of in vivo imaging when spatial encoding is performed with three or more encoding channels for a 2D image. MATERIALS AND METHODS: The recently developed 4-Dimensional Radial In/Out (4D-RIO) trajectory is compared in simulations to an alternative higher-order encoding scheme referred to as O-space imaging. Direct comparison of local k-space representations leads to the proposal of a modification to the O-space imaging trajectory based on a scheme of prephasing to improve the reconstructed image quality. Data were collected using a 4D-RIO acquisition in vivo in the human brain and several image reconstructions were compared, exploiting the property that the dense encoding matrix, after a 1D or 2D Fourier transform, can be approximated by a sparse matrix by discarding entries below a chosen magnitude. RESULTS: The proposed prephasing scheme for the O-space trajectory shows a marked improvement in quality in the simulated image reconstruction. In experiments, 4D-RIO data acquired in vivo in the human brain can be reconstructed to a reasonable quality using only 5 % of the encoding matrix--massively reducing computer memory requirements for a practical reconstruction. CONCLUSION: Trajectory design and reconstruction techniques such as these may prove especially useful when extending generalized higher-order encoding methods to 3D images.

Reconstruction of MRI data encoded by multiple nonbijective curvilinear magnetic fields.

Parallel imaging technique using localized gradients (PatLoc) uses the combination of surface gradient coils generating nonbijective curvilinear magnetic fields for spatial encoding. PatLoc imaging using one pair of multipolar spatial encoding magnetic fields (SEMs) has two major caveats: (1) The direct inversion of the encoding matrix requires exact determination of multiple locations which are ambiguously encoded by the SEMs. (2) Reconstructed images have a prominent loss of spatial resolution at the center of field-of-view using a symmetric coil array for signal detection. This study shows that a PatLoc system actually has a higher degree of freedom in spatial encoding to mitigate the two challenges mentioned above. Specifically, a PatLoc system can generate not only multipolar but also linear SEMs, which can be used to reduce the loss of spatial resolution at the field-of-view center. Here, we present an efficient and generalized image reconstruction method for PatLoc imaging using multiple SEMs without explicitly identifying the locations where SEM encoding is not unique. Reconstructions using simulations and empirical experimental data are compared with those using conventional linear gradients to demonstrate that the general combination of SEMs can improve image reconstructions.

Localization by nonlinear phase preparation and k-space trajectory design.

A technique is described to localize MR signals from a target volume using nonlinear pulsed magnetic fields and spatial encoding trajectories designed using local k-space theory. The concept of local k-space is outlined theoretically, and this principle is applied to simulated phantom and cardiac MRI data in the presence of surface and quadrupolar gradient coil phase modulation. Phantom and in vivo human brain images are obtained using a custom, high-performance quadrupolar gradient coil integrated with a whole-body 3-T MRI system to demonstrate target localization using three-dimensional T 2*-weighted spoiled gradient echo, two-dimensional segmented, multiple gradient encoded spin echo, and three-dimensional balanced steady-state free precession acquisitions. This method may provide a practical alternative to selective radiofrequency excitation at ultra-high-field, particularly for steady-state applications where repetition time (TR) must be minimized and when the amount of energy deposited in human tissues is prohibitive. There are several limitations to the approach including the spatial variation in resolution, high frequency aliasing artifacts, and spatial variation in echo times and contrast.

Radial imaging with multipolar magnetic encoding fields.

We present reconstruction methods for radial magnetic resonance imaging (MRI) data which were spatially encoded using a pair of orthogonal multipolar magnetic fields for in-plane encoding and parallel imaging. It is shown that a direct method exists in addition to iterative reconstruction. Standard direct projection reconstruction algorithms can be combined with a previously developed direct reconstruction for multipolar encoding fields acquired with Cartesian trajectories. The algorithm is simplified by recasting the reconstruction problem into polar coordinates. In this formulation distortion and aliasing become separate effects. Distortion occurs only along the radial direction and aliasing along the azimuthal direction. Moreover, aliased points are equidistantly distributed in this representation, and, consequently, Cartesian SENSE is directly applicable with highly effective unfolding properties for radio-frequency coils arranged with a radial symmetry. The direct and iterative methods are applied to simulated data to analyze basic properties of the algorithms and for the first time also measured in vivo data are presented. The results are compared to linear spatial encoding using a radial trajectory and quadrupolar encoding using a Cartesian trajectory. The direct reconstruction gives good results for fully sampled datasets. Undersampled datasets, however, show star-shaped artifacts, which are significantly reduced with the iterative reconstruction.

Simultaneously driven linear and nonlinear spatial encoding fields in MRI.

Spatial encoding in MRI is conventionally achieved by the application of switchable linear encoding fields. The general concept of the recently introduced PatLoc (Parallel Imaging Technique using Localized Gradients) encoding is to use nonlinear fields to achieve spatial encoding. Relaxing the requirement that the encoding fields must be linear may lead to improved gradient performance or reduced peripheral nerve stimulation. In this work, a custom-built insert coil capable of generating two independent quadratic encoding fields was driven with high-performance amplifiers within a clinical MR system. In combination with the three linear encoding fields, the combined hardware is capable of independently manipulating five spatial encoding fields. With the linear z-gradient used for slice-selection, there remain four separate channels to encode a 2D-image. To compare trajectories of such multidimensional encoding, the concept of a local k-space is developed. Through simulations, reconstructions using six gradient-encoding strategies were compared, including Cartesian encoding separately or simultaneously on both PatLoc and linear gradients as well as two versions of a radial-based in/out trajectory. Corresponding experiments confirmed that such multidimensional encoding is practically achievable and demonstrated that the new radial-based trajectory offers the PatLoc property of variable spatial resolution while maintaining finite resolution across the entire field-of-view.

Reconstruction of MRI data encoded with arbitrarily shaped, curvilinear, nonbijective magnetic fields.

A basic framework for image reconstruction from spatial encoding by curvilinear, nonbijective magnetic encoding fields in combination with multiple receivers is presented. The theory was developed in the context of the recently introduced parallel imaging technique using localized gradients (PatLoc) approach. In this new imaging modality, the linear gradient fields are generalized to arbitrarily shaped, nonbijective spatial encoding magnetic fields, which lead to ambiguous encoding. Ambiguities are resolved by adaptation of concepts developed for parallel imaging. Based on theoretical considerations, a practical algorithm for Cartesian trajectories is derived in the case that the conventional gradient coils are replaced by coils for PatLoc. The reconstruction method extends Cartesian sensitivity encoding (SENSE) reconstruction with an additional voxelwise intensity-correction step. Spatially varying resolution, signal-to-noise ratio, and truncation artifacts are described and analyzed. Theoretical considerations are validated by two-dimensional simulations based on multipolar encoding fields and they are confirmed by applying the reconstruction algorithm to initial experimental data.

Parallel imaging in non-bijective, curvilinear magnetic field gradients: a concept study.

OBJECTIVES: The paper presents a novel and more generalized concept for spatial encoding by non-unidirectional, non- bijective spatial encoding magnetic fields (SEMs). In combination with parallel local receiver coils these fields allow one to overcome the current limitations of neuronal nerve stimulation. Additionally the geometry of such fields can be adapted to anatomy. MATERIALS AND METHODS: As an example of such a parallel imaging technique using localized gradients (PatLoc)- system, we present a polar gradient system consisting of 2 x 8 rectangular current loops in octagonal arrangement, which generates a radial magnetic field gradient. By inverting the direction of current in alternating loops, a near sinusoidal field variation in the circumferential direction is produced. Ambiguities in spatial assignment are resolved by use of multiple receiver coils and parallel reconstruction. Simulations demonstrate the potential advantages and limitations of this approach. RESULTS AND CONCLUSIONS: The exact behaviour of PatLoc fields with respect to peripheral nerve stimulation needs to be tested in practice. Based on geometrical considerations SEMs of radial geometry allow for about three times faster gradient switching compared to conventional head gradient inserts and even more compared to whole body gradients. The strong nonlinear geometry of the fields needs to be considered for practical applications.