Performance analysis for magnetic resonance imaging with nonlinear encoding fields

Use of nonlinear pulsed magnetic fields for spatial encoding in magnetic resonance imaging

This study examines the use of nonlinear magnetic field coils for spatial encoding in magnetic resonance imaging. Existing theories on imaging with such coils share a complex reconstruction process that originates from a suboptimal signal interpretation in the spatial-frequency domain (k-space). In this study, a new solution to this problem is proposed, namely a two-step reconstruction process, in which in the first step, the image signal is converted into a frequency spectrum, and in the second step, the spectrum, which represents the distorted image, is geometrically and intensity corrected to obtain an undistorted image. This theory has been verified by numerical simulations and experimentally using a straight wire as a coil model for an extremely nonlinear magnetic field. The results of this study facilitate the use of simple encoding coil designs that can feature low inductance, allowing for much faster switching times and higher magnetic field gradients.

Encoding scheme design for gradient-free, nonlinear projection imaging using Bloch-Siegert RF spatial encoding in a low-field, open MRI system

Eliminating conventional pulsed B0-gradient coils for magnetic resonance imaging (MRI) can significantly reduce the cost of and increase access to these devices. Phase shifts induced by the Bloch-Siegert shift effect have been proposed as a means for gradient-free, RF spatial encoding for low-field MR imaging. However, nonlinear phasor patterns like those generated from loop coils have not been systematically studied in the context of 2D spatial encoding. This work presents an optimization algorithm to select an efficient encoding trajectory among the nonlinear patterns achievable with a given hardware setup. Performance of encoding trajectories or projections was evaluated through simulated and experimental image reconstructions. Results show that the encodings schemes designed by this algorithm provide more efficient spatial encoding than comparison encoding sets, and the method produces images with the predicted spatial resolution and minimal artifacts. Overall, the work demonstrates the feasibility of performing 2D gradient-free, low-field imaging using the Bloch-Siegert shift which is an important step towards creating low-cost, point-of-care MR systems.

Accelerated 2D Cartesian MRI with an 8-channel local B0 coil array combined with parallel imaging

PURPOSE

In MRI, the magnetization of nuclear spins is spatially encoded with linear gradients and radiofrequency receivers sensitivity profiles to produce images, which inherently leads to a long scan time. Cartesian MRI, as widely adopted for clinical scans, can be accelerated with parallel imaging and rapid magnetic field modulation during signal readout. Here, by using an 8-channel localB 0 $$ {\mathrm{B}}_0 $$ coil array, the modulation scheme optimized for sampling efficiency is investigated to speed up 2D Cartesian scans.

THEORY AND METHODS

An 8-channel localB 0 $$ {\mathrm{B}}_0 $$ coil array is made to carry sinusoidal currents during signal readout to accelerate 2D Cartesian scans. An MRI sampling theory based on reproducing kernel Hilbert space is exploited to visualize the efficiency of nonlinear encoding in arbitrary sampling duration. A field calibration method using current monitors for localB 0 $$ {\mathrm{B}}_0 $$ coils and the ESPIRiT algorithm is proposed to facilitate image reconstruction. Image acceleration with various modulation field shapes, aliasing control, and distinct modulation frequencies are scrutinized to find an optimized modulation scheme. A safety evaluation is conducted. In vivo 2D Cartesian scans are accelerated by the localB 0 $$ {\mathrm{B}}_0 $$ coils.

RESULTS

For 2D Cartesian MRI, the optimal modulation field by this localB 0 $$ {\mathrm{B}}_0 $$ array converges to a nearly linear gradient field. With the field calibration technique, it accelerates the in vivo scans (i.e., proved safe) by threefold and eightfold free of visible artifacts, without and with SENSE, respectively.

CONCLUSION

The nonlinear encoding analysis tool, the field calibration method, the safety evaluation procedures, and the in vivo reconstructed scans make significant steps to push MRI speed further with the localB 0 $$ {\mathrm{B}}_0 $$ coil array.

Design of sparse Halbach magnet arrays for portable MRI using a genetic algorithm

Permanent magnet arrays offer several attributes attractive for the development of a low-cost portable MRI scanner for brain imaging. They offer the potential for a relatively lightweight, low to mid-field system with no cryogenics, a small fringe field, and no electrical power requirements or heat dissipation needs. The cylindrical Halbach array, however, requires external shimming or mechanical adjustments to produce B0 fields with standard MRI homogeneity levels (e.g., 0.1 ppm over FOV), particularly when constrained or truncated geometries are needed, such as a head-only magnet where the magnet length is constrained by the shoulders. For portable scanners using rotation of the magnet for spatial encoding with generalized projections, the spatial pattern of the field is important since it acts as the encoding field. In either a static or rotating magnet, it will be important to be able to optimize the field pattern of cylindrical Halbach arrays in a way that retains construction simplicity. To achieve this, we present a method for designing an optimized cylindrical Halbach magnet using the genetic algorithm to achieve either homogeneity (for standard MRI applications) or a favorable spatial encoding field pattern (for rotational spatial encoding applications). We compare the chosen designs against a standard, fully populated sparse Halbach design, and evaluate optimized spatial encoding fields using point-spread-function and image simulations. We validate the calculations by comparing to the measured field of a constructed magnet. The experimentally implemented design produced fields in good agreement with the predicted fields, and the genetic algorithm was successful in improving the chosen metrics. For the uniform target field, an order of magnitude homogeneity improvement was achieved compared to the un-optimized, fully populated design. For the rotational encoding design the resolution uniformity is improved by 95% compared to a uniformly populated design.

O-space with high resolution readouts outperforms radial imaging

PURPOSE

While O-Space imaging is well known to accelerate image acquisition beyond traditional Cartesian sampling, its advantages compared to undersampled radial imaging, the linear trajectory most akin to O-Space imaging, have not been detailed. In addition, previous studies have focused on ultrafast imaging with very high acceleration factors and relatively low resolution. The purpose of this work is to directly compare O-Space and radial imaging in their potential to deliver highly undersampled images of high resolution and minimal artifacts, as needed for diagnostic applications. We report that the greatest advantages to O-Space imaging are observed with extended data acquisition readouts.

THEORY AND METHODS

A sampling strategy that uses high resolution readouts is presented and applied to compare the potential of radial and O-Space sequences to generate high resolution images at high undersampling factors. Simulations and phantom studies were performed to investigate whether use of extended readout windows in O-Space imaging would increase k-space sampling and improve image quality, compared to radial imaging.

RESULTS

Experimental O-Space images acquired with high resolution readouts show fewer artifacts and greater sharpness than radial imaging with equivalent scan parameters. Radial images taken with longer readouts show stronger undersampling artifacts, which can cause small or subtle image features to disappear. These features are preserved in a comparable O-Space image.

CONCLUSIONS

High resolution O-Space imaging yields highly undersampled images of high resolution and minimal artifacts. The additional nonlinear gradient field improves image quality beyond conventional radial imaging.

Experimental O-space turbo spin echo imaging

PURPOSE

Turbo spin echo (TSE) imaging reduces imaging time by acquiring multiple echoes per repetition (TR), requiring fewer TRs. O-space can also require fewer TRs by using a combination of nonlinear magnetic gradient fields and surface coil arrays. Although to date, O-space has only been demonstrated for gradient echo imaging, it is valuable to combine these two techniques. However, collecting multiple O-space echoes per TR is difficult because of the different local k-space trajectories and variable T2-weighting.

THEORY AND METHODS

A practical scheme is demonstrated to combine the benefits of TSE and O-space for highly accelerated T2-weighted images. The scheme uses a modified acquisition order and filtered projection reconstruction to reduce artifacts caused by T2 decay, while retaining T2 contrast that corresponds to a specific echo time.

RESULTS

The experiments revealed that the proposed method can produce highly accelerated T2-weighted images. Moreover, the method can generate multiple images with different T2 contrasts from a single dataset.

CONCLUSIONS

The proposed O-space TSE imaging method requires fewer echoes than conventional TSE and fewer repetitions than conventional O-space imaging. It retains resilience to undersampling, clearly outperforming Cartesian SENSE at high levels of undersampling, and can generate undistorted images with a range of T2 contrast from a single acquired dataset.

Fast rotary nonlinear spatial acquisition (FRONSAC) imaging

PURPOSE

Nonlinear spatial encoding magnetic fields (SEMs) have been studied to reconstruct images from a minimum number of echoes. Previous work has also explored single shot trajectories in nonlinear SEMs. However, the search continues for optimal schemes that apply nonlinear SEMs to improve spatial encoding efficiency and image quality.

THEORY AND METHODS

We enhance the encoding efficiency of standard linear gradient trajectories by adding a rapidly rotating nonlinear SEM of moderate amplitude, the so called FRONSAC (Fast ROtary Nonlinear Spatial ACquisition) imaging. This additional gradient greatly improves the image quality of highly undersampled single-shot trajectories, including EPI, Spiral, and Rosette trajectories.

RESULTS

Our simulations, including noise and dephasing effects, test the effect of adding FRONSAC gradients, demonstrating the applicability of this approach. Performance is explained by demonstrating the additional k-space sampling the nonlinear gradient provides. Studies of the optimal amplitude and frequency of the additional FRONSAC field are presented, and the role of enhanced sampling during the readout demonstrated. Dynamic field mapping in a second-order gradient system shows the proposed gradient waveforms are feasible.

CONCLUSION

Images resulting from highly undersampled existing k-space trajectories, such as EPI, Spiral, and Rosette, are greatly enhanced simply by adding a rotating nonlinear SEM field.

Trajectory optimization based on the signal-to-noise ratio for spatial encoding with nonlinear encoding fields

PURPOSE

Multiple nonlinear gradient fields offer many potential benefits for spatial encoding including reduced acquisition time, fewer artefacts and region-specific imaging, although designing a suitable trajectory for such a setup is difficult. This work aims to optimize encoding trajectories for multiple nonlinear gradient fields based on the image signal-to-noise ratio.

THEORY AND METHODS

Image signal-to-noise ratio is directly linked to the covariance of the reconstructed pixels, which can be calculated recursively for each projection of the trajectory under a Bayesian formulation. An evolutionary algorithm is used to find the higher-dimensional projections that minimize the pixel covariance, incorporating receive coil profiles, intravoxel dephasing, and reconstruction regularization. The resulting trajectories are tested through simulations and experiments.

RESULTS

The optimized trajectories produce images with higher resolution and fewer artefacts compared with traditional approaches, particularly for high undersampling. However, higher-dimensional projection experiments strongly depend on accurate hardware and calibration.

CONCLUSION

Computer-based optimization provides an efficient means to explore the large trajectory space created by the use of multiple nonlinear encoding fields. The optimization framework, as presented here, is necessary to fully exploit the advantages of nonlinear fields. Magn Reson Med 76:104-117, 2016. © 2015 Wiley Periodicals, Inc.

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.

Pseudo-random center placement O-space imaging for improved incoherence compressed sensing parallel MRI

PURPOSE

Nonlinear spatial encoding magnetic (SEM) field strategies such as O-space imaging have previously reported dispersed artifacts during accelerated scans. Compressed sensing (CS) has shown a sparsity-promoting convex program allows image reconstruction from a reduced data set when using the appropriate sampling. The development of a pseudo-random center placement (CP) O-space CS approach optimizes incoherence through SEM field modulation to reconstruct an image with reduced error.

THEORY AND METHODS

The incoherence parameter determines the sparsity levels for which CS is valid and the related transform point spread function measures the maximum interference for a single point. The O-space acquisition is optimized for CS by perturbing the Z(2) strength within 30% of the nominal value and demonstrated on a human 3T scanner.

RESULTS

Pseudo-random CP O-space imaging is shown to improve incoherence between the sensing and sparse domains. Images indicate pseudo-random CP O-space has reduced mean squared error compared with a typical linear SEM field acquisition method.

CONCLUSION

Pseudo-random CP O-space imaging, with a nonlinear SEM field designed for CS, is shown to reduce mean squared error of images at high acceleration over linear encoding methods for a 2D slice when using an eight channel circumferential receiver array for parallel imaging.

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.

Spatial encoding using the nonlinear field perturbations from magnetic materials

PURPOSE

A proof-of-concept study was performed to assess the technical feasibility of using magnetic materials to generate spatial encoding fields.

THEORY AND METHODS

Spatially varying magnetic fields were generated by the placement of markers with different volume susceptibilities within the imaging volume. No linear gradients were used for spatial encoding during the signal acquisition. A signal-encoding model is described for reconstructing the images encoded with these field perturbations. Simulation and proof-of-concept experimental results are presented. Experiments were performed using field perturbations from a cylindrical marker as an example of the new encoding fields. Based on this experimental setup, annular rings were reconstructed from signals encoded with the new fields.

RESULTS

Simulation results were presented for different acquisition parameters. Proof-of-concept was supported by the correspondence of regions in an image reconstructed from experimental data compared to those in a conventional gradient-echo image. Experimental results showed that inclusions of dimensions 1.5 mm in size could be resolved with the experimental setup.

CONCLUSION

This study shows the technical feasibility of using magnetic markers to produce encoding fields. Magnetic materials will allow generating spatial encoding fields, which can be tailored to an imaging application with less complexity and at lower cost compared to the use of gradient inserts.

In vivo O-Space imaging with a dedicated 12 cm Z2 insert coil on a human 3T scanner using phase map calibration

Recently, spatial encoding with nonlinear magnetic fields has drawn attention for its potential to achieve faster gradient switching within safety limits, tailored resolution in regions of interest, and improved parallel imaging using encoding fields that complement the sensitivity profiles of radio frequency receive arrays. Proposed methods can broadly be divided into those that use phase encoding (Cartesian-trajectory PatLoc and COGNAC) and those that acquire nonlinear projections (O-Space, Null space imaging, radial PatLoc, and 4D-RIO). Nonlinear projection data are most often reconstructed with iterative algorithms that backproject data using the full encoding matrix. Just like conventional radial sequences that use linear spatial encoding magnetic fields, nonlinear projection methods are more sensitive than phase encoding methods to imperfect calibration of the encoding fields. In this work, voxel-wise phase evolution is mapped at each acquired point in an O-Space trajectory using a variant of chemical shift imaging, capturing all spin dynamics caused by encoding fields, eddy currents, and pulse timing. Phase map calibration is then applied to data acquired from a high-power, 12 cm, Z2 insert coil with an eight-channel radio frequency transmit-receive array on a 3T human scanner. We show the first experimental proof-of-concept O-Space images on in vivo and phantom samples, paving the way for more in-depth exploration of O-Space and similar imaging methods.

Single shot trajectory design for region-specific imaging using linear and nonlinear magnetic encoding fields

It has recently been demonstrated that nonlinear encoding fields result in a spatially varying resolution. This work develops an automated procedure to design single-shot trajectories that create a local resolution improvement in a region of interest. The technique is based on the design of optimized local k-space trajectories and can be applied to arbitrary hardware configurations that employ any number of linear and nonlinear encoding fields. The trajectories designed in this work are tested with the currently available hardware setup consisting of three standard linear gradients and two quadrupolar encoding fields generated from a custom-built gradient insert. A field camera is used to measure the actual encoding trajectories up to third-order terms, enabling accurate reconstructions of these demanding single-shot trajectories, although the eddy current and concomitant field terms of the gradient insert have not been completely characterized. The local resolution improvement is demonstrated in phantom and in vivo experiments.

The Role of Nonlinear Gradients in Parallel Imaging: A k-Space Based Analysis

Sequences that encode the spatial information of an object using nonlinear gradient fields are a new frontier in MRI, with potential to provide lower peripheral nerve stimulation, windowed fields of view, tailored spatially-varying resolution, curved slices that mirror physiological geometry, and, most importantly, very fast parallel imaging with multichannel coils. The acceleration for multichannel images is generally explained by the fact that curvilinear gradient isocontours better complement the azimuthal spatial encoding provided by typical receiver arrays. However, the details of this complementarity have been more difficult to specify. We present a simple and intuitive framework for describing the mechanics of image formation with nonlinear gradients, and we use this framework to review some the main classes of nonlinear encoding schemes.