Application of keyhole imaging to interventional MRI: A simulation study to predict sequence requirements

A wide variety of techniques have been proposed recently to improve the temporal resolution of MRI. These include echo-planar imaging methods, wavelet encoding, singular value decomposition encoding, and k-space sharing methods known as ?keyhole”? imaging. In this work, we use a simulation study to investigate the phase-encoding ordering and data-sharing methods required for the application of keyhole imaging to interventional MRI (I-MRI). The advantages of keyhole imaging over other methods are its simplicity and the use of conventional phase encoding and Fourier transform reconstruction found on virtually all modern MR imagers. Our analysis has predicted that conventional keyhole methods that repeatedly acquire only the center portion of k space, and those that sequentially progress from the center of k space outward, will not meet the combination of temporal and spatial resolution required for tip localization during I-MRI needle insertion. Instead, acquisitions that acquire both high and low k-space data, in ranked order, should provide acceptable tip position and needle width accuracy in both temporal and spatial domains for use in I-MRI.     

Dynamic imaging of contrast-enhanced coronary vessels with a magnetization prepared rotated stripe keyhole acquisition

PURPOSE: To evaluate dynamic coronary imaging based on a magnetization prepared contrast-enhanced (CE) rotated stripe keyhole acquisition scheme. MATERIALS AND METHODS: Background suppression of long T(1) tissue was used so that the k-space would be selectively dominated by the contribution of the CE vessel. The phase-encoding axis was then adjusted parallel to the long axis of the vessel to sample the significant power spectrum of the vessel. The performance of this approach was evaluated by means of computer simulations and experimental studies on phantoms and a pig model instrumented with an intracoronary catheter for infusion of contrast media. RESULTS: Computer simulations and phantom studies demonstrated that by rotating the gradient axes to match the k-space pattern of the frequency spectrum, one can reduce the keyhole band to 20% of the full k-space while preserving the structure's lumen width and sharpness. In vivo studies validated those findings, and dynamic angiograms of the CE coronary arteries were obtained as rapidly as 140 msec per image, with an in-plane spatial resolution of 1.5 mm. CONCLUSION: With efficient background suppression, a rotated stripes keyhole acquisition can efficiently acquire the significant k-space of a CE vessel, and provide improved vessel definition with a reduced acquisition matrices scheme.     

Fast “real time” imaging with different k-space update strategies for interventional procedures

Interventional procedures under MR guidance require the images to be acquired with a fast acquisition strategy, a rapid reconstruction algorithm for “real-time” imaging (ie, high temporal resolution), acquisition of at least three adjacent slices to track a tool reliably, and high tissue contrast to ensure safe positioning of interventional devices. Often times, the field strength for interventional MR-imaging units is limited by the open magnet design. This complicates the trade-off between scan time and image quality, particularly when applied during low field interventional MRI procedures. To minimize the impact of some of these tradeoffs, a combination of keyhole techniques or modified k-space trajectories, in conjunction with a fluoroscopic (ie, continuous acquisition) mode and a real time reconstruction, permits rapid imaging in a low field system using standard (speed optimized) reconstruction hardware and standard gradient electronics. The purpose of this study was to design and describe different keyhole strategies that can be used in a real time mode to increase the image frame rate by a factor of up to 16. By updating the entire raw data space with our strategies, even small changes of the object could be recognized. Our results using these new strategies on two commercially available open magnet MR-imaging units (Siemens Magnetom Open 0.2T resistive magnet, Toshiba Access 0.064T permanent magnet) and a 1.5T superconductive solenoidal magnet design imager (Siemens SP) are presented to show the potential of these acquisition strategies in interventional MRI. Furthermore, these strategies may also be helpful for several other medical applications requiring high temporal resolution like contrast-enhanced breast imaging or functional brain imaging.     

Functional magnetic resonance imaging using PROPELLER-EPI

Periodically rotated overlapping parallel lines with enhanced reconstruction-echo-planar imaging (PROPELLER-EPI) is a multishot technique that samples k-space by acquisition of narrow blades, which are subsequently rotated until the entire k-space is filled. It has the unique advantage that the center of k-space, and thus the area containing the majority of functional MRI signal changes, is sampled with each shot. This continuous refreshing of the k-space center by each acquired blade enables not only sliding-window but also keyhole reconstruction. Combining PROPELLER-EPI with a fast gradient-echo readout scheme allows for high spatial resolutions to be achieved while maintaining a temporal resolution, which is suitable for functional MRI experiments. Functional data acquired with a novel interlaced sequence that samples both single-shot EPI and blades in an alternating fashion suggest that PROPELLER-EPI can achieve comparable functional MRI results. PROPELLER-EPI, however, features different spatiotemporal characteristics than single-shot EPI, which not only enables keyhole reconstruction but also makes it an interesting alternative for many functional MRI applications.     

Keyhole Dixon method for faster,perceptually equivalent fat suppression

PURPOSE: To reduce the acquisition time associated with the two-point Dixon fat suppression technique by combining a keyhole in-phase (Water + Fat) k-space data set with a full out-of-phase (Water - Fat) k-space data set and optimizing the keyhole size with a perceptual difference model. MATERIALS AND METHODS: A set of keyhole Dixon images was created by varying the number of lines in the keyhole data set. Off-resonance correction was incorporated into the image reconstruction process to improve the homogeneity of the fat suppression. A perceptual difference model (PDM) was validated with human observer experiments and used to compare the keyhole images to images from a full two-point Dixon acquisition. The PDM was used to determine the smallest keyhole width required to obtain perceptual equivalence to images obtained from the full two-point Dixon method. RESULTS: In experimental phantom studies, the keyhole Dixon image reconstructed from 96 of 192 Water + Fat k-space lines and 192 Water - Fat k-space lines was perceptually equivalent to the full (192 + 192) two-point Dixon images, resulting in a 25% reduction in scan time. Clinical images of a volunteer's knee, orbits, and abdomen created from the smallest, perceptually equivalent keyhole width resulted in a 27%-38% reduction in total scan time. CONCLUSION: This method improves the temporal efficiency of the conventional two-point Dixon technique and may prove especially useful for high-field systems where specific absorption rate (SAR) limits will constrain radiofrequency (RF)-based fat suppression techniques.     

不同场强影响介入性磁共振穿刺针成像因素的初步探讨

目的　探讨 0 .3T、1.5T场强下影响介入性磁共振穿刺针成像的因素。以便提供正确的穿刺方法及成像技术参数 ,确保介入性磁共振操作的安全性及准确性。方法　将MR相容性介入穿刺针置于琼脂模型内 ,在 0 .3T、1.5T场强下 ,采用自旋回波序列、快速自旋回波序列、梯度回波序列 ,将穿刺针平行、垂直于主磁场 ,频率编码方向平行或垂直于穿刺针的方向进行扫描 ,在工作站测量图像上穿刺针的宽度及针尖的位置 ,比较图像上和实际针尖位置的差异 ,并比较不同场强下影响穿刺针伪影的因素。结果　各序列中梯度回波序列产生的伪影较大 ,快速自旋回波序列产生的伪影较小 ,自选回波序列产生的伪影最小 ,但快速自旋回波序列与自旋回波序列产生的伪影无明显差异。快速自旋回波序列与自旋回波序列中 ,当频率编码轴垂直于穿刺针长轴时伪影较大 ;当穿刺针方向逐渐平行于主磁场方向时 ,伪影逐渐减小。所有图像上针尖位置与实际针尖位置的差异在 1cm内。结论　 1.5T场强下伪影宽度大于 0 .3T场强。在 0 .3T、1.5T场强下改变频率编码方向、脉冲序列、成像参数时 ,针尖位置的变化在 1cm内 ,穿刺针平行于主磁场方向时 ,在不同场强下明显降低伪影的宽度     

Parallel magnetic resonance imaging with adaptive radius in k-space (PARS): constrained image reconstruction using k-space locality in radiofrequency coil encoded data.

A parallel image reconstruction algorithm is presented that exploits the k-space locality in radiofrequency (RF) coil encoded data. In RF coil encoding, information relevant to reconstructing an omitted datum rapidly diminishes as a function of k-space separation between the omitted datum and the acquired signal data. The proposed method, parallel magnetic resonance imaging with adaptive radius in k-space (PARS), harnesses this physical property of RF coil encoding via a sliding-kernel approach. Unlike generalized parallel imaging approaches that might typically involve inverting a prohibitively large matrix for arbitrary sampling trajectories, the PARS sliding-kernel approach creates manageable and distributable independent matrices to be inverted, achieving both computational efficiency and numerical stability. An empirical method designed to measure total error power is described, and the total error power of PARS reconstructions is studied over a range of k-space radii and accelerations, revealing "minimal-error" conditions at comparatively modest k-space radii. PARS reconstructions of undersampled in vivo Cartesian and non-Cartesian data sets are shown and are compared selectively with traditional SENSE reconstructions. Various characteristics of the PARS k-space locality constraint (such as the tradeoff between signal-to-noise ratio and artifact power and the relationship with iterative parallel conjugate gradient approaches or nonparallel gridding approaches) are discussed.     

Averaging keyhole pulse sequence with presaturation pulses and EXORCYCLE phase cycling for dynamic contrast-enhanced MRI.

The purpose of this study was to design a keyhole pulse sequence for quantitative 2D dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) based on a spoiled gradient echo T1-weighted acquisition. Saturation recovery was applied to achieve a linear correlation between signal intensity and contrast agent concentration in an arterial input function (AIF) while simultaneously removing time-of-flight effect. To remove ghosting artifacts arising from incomplete presaturation, EXORCYCLE phase cycling with averaging was applied to the pulse sequence. RF spoiling by radiofrequency switching with the synthesizer can be combined with EXORCYCLE phase cycling. Images affected by the large difference in signal intensity before and after contrast agent administration with the keyhole technique were improved by interleaving of peripheral lines of k-space with groups of central lines. Both peripheral and central lines were renewed during the dynamic scan. AIFs were obtained from the rat abdominal aorta with this keyhole sequence.     

Transmit SENSE.

The idea of using parallel imaging to shorten the acquisition time by the simultaneous use of multiple receive coils can be adapted for the parallel transmission of a spatially-selective multidimensional RF pulse. As in data acquisition, a multidimensional RF pulse follows a certain k-space trajectory. Shortening this trajectory shortens the pulse duration. The use of multiple transmit coils, each with its own time-dependent waveform and spatial sensitivity, can compensate for the missing parts of the excitation k-space. This results in a maintained spatial definition of the pulse profile, while its duration is reduced. This work introduces the concept of parallel transmission with arbitrarily shaped transmit coils (termed "Transmit SENSE"). Results of numerical studies demonstrate the theoretical feasibility of the approach. The experimental proof of principle is provided on a commercial MR scanner. The lack of multiple independent transmit channels was addressed by combining the excitation patterns from two separate subexperiments with different transmit setups. Shortening multidimensional RF pulses could be an interesting means of making 3D RF pulses feasible even for fast T(2)(*) relaxing species or strong main field inhomogeneities. Other applications might benefit from the ability of Transmit SENSE to improve the spatial resolution of the pulse profile while maintaining the transmit duration.     

Non-Fourier-encoded parallel MRI using multiple receiver coils.

This paper describes a general theoretical framework that combines non-Fourier (NF) spatially-encoded MRI with multichannel acquisition parallel MRI. The two spatial-encoding mechanisms are physically and analytically separable, which allows NF encoding to be expressed as complementary to the inherent encoding imposed by RF receiver coil sensitivities. Consequently, the number of NF spatial-encoding steps necessary to fully encode an FOV is reduced. Furthermore, by casting the FOV reduction of parallel imaging techniques as a dimensionality reduction of the k-space that is NF-encoded, one can obtain a speed-up of each digital NF spatial excitation in addition to accelerated imaging. Images acquired at speed-up factors of 2x to 8x with a four-element RF receiver coil array demonstrate the utility of this framework and the efficiency afforded by it.     

MR coronary vessel wall imaging: comparison between radial and spiral k-space sampling

PURPOSE: To compare radial and spiral k-space sampling in navigator-gated ECG-triggered three-dimensional (3D) coronary vessel wall imaging. MATERIALS AND METHODS: The right coronary artery (RCA) vessel walls of eight healthy subjects were imaged using a modified double-inversion prepulse in concert with radial and spiral data acquisition. For data analysis, two investigators blinded to the sequence parameters subjectively assessed image quality in terms of artifacts and vessel wall visualization. Objective measures of the signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), and vessel wall definition were also determined. RESULTS: Radial k-space sampling demonstrated fewer artifacts and led to improved visualization of the coronary vessel wall compared to spiral imaging (P < 0.05). This finding was also reflected in a better vessel wall definition using radial data acquisition (P < 0.05). SNR and CNR were found to be higher when spiral k-space sampling was used (n.s.). CONCLUSION: Radial k-space sampling in concert with free-breathing navigator-gated cardiac-triggered MRI of the coronary vessel wall resulted in fewer motion artifacts and improved vessel wall definition compared to spiral k-space sampling. The proposed approach therefore appears to be preferable.     

Inherently self-calibrating non-Cartesian parallel imaging.

The use of self-calibrating techniques in parallel magnetic resonance imaging eliminates the need for coil sensitivity calibration scans and avoids potential mismatches between calibration scans and subsequent accelerated acquisitions (e.g., as a result of patient motion). Most examples of self-calibrating Cartesian parallel imaging techniques have required the use of modified k-space trajectories that are densely sampled at the center and more sparsely sampled in the periphery. However, spiral and radial trajectories offer inherent self-calibrating characteristics because of their densely sampled center. At no additional cost in acquisition time and with no modification in scanning protocols, in vivo coil sensitivity maps may be extracted from the densely sampled central region of k-space. This work demonstrates the feasibility of self-calibrated spiral and radial parallel imaging using a previously described iterative non-Cartesian sensitivity encoding algorithm.     

4D time-resolved MR angiography with keyhole (4D-TRAK): more than 60 times accelerated MRA using a combination of CENTRA, keyhole, and SENSE at 3.0T

PURPOSE: To present a new 4D method that is designed to provide high spatial resolution MR angiograms at subsecond temporal resolution by combining different techniques of view sharing with parallel imaging at 3.0T. MATERIALS AND METHODS: In the keyhole-based method, a central elliptical cylinder in k-space is repeated n times (keyhole) with a random acquisition (CENTRA), and followed by the readout of the periphery of k-space. 4D-MR angiography with CENTRA keyhole (4D-TRAK) was combined with parallel imaging (SENSE) and partial Fourier imaging. In total, a speed-up factor of 66.5 (6.25 [CENTRA keyhole] x 8 [SENSE] x 1.33 [partial Fourier imaging]) was achieved yielding a temporal resolution of 608 ms and a spatial resolution of (1.1 x 1.4 x 1.1) mm(3) with whole-brain coverage 4D-TRAK was applied to five patients and compared with digital subtraction angiography (DSA). RESULTS: 4D-TRAK was successfully completed with an acceleration factor of 66.5 in all five patients. Sharp images were acquired without any artifacts possibly created by the transition of the central cylinder and the reference dataset. MRA findings were concordant with DSA. CONCLUSION: 4D time-resolved MRA with keyhole (4D-TRAK) is feasible using a combination of CENTRA, keyhole, and SENSE at 3.0T and allows for more than 60 times accelerated MRA with high spatial resolution.     

Parallel excitation with an array of transmit coils.

Theoretical and experimental results are presented that establish the value of parallel excitation with a transmit coil array in accelerating excitation and managing RF power deposition. While a 2D or 3D excitation pulse can be used to induce a multidimensional transverse magnetization pattern for a variety of applications (e.g., a 2D localized pattern for accelerating spatial encoding during signal acquisition), it often involves the use of prolonged RF and gradient pulses. Given a parallel system that is composed of multiple transmit coils with corresponding RF pulse synthesizers and amplifiers, the results suggest that by exploiting the localization characteristics of the coils, an orchestrated play of shorter RF pulses can achieve desired excitation profiles faster without adding strains to gradients. A closed-form design for accelerated multidimensional excitations is described for the small-tip-angle regime, and its suppression of interfering aliasing lobes from coarse excitation k-space sampling is interpreted based on an analogy to sensitivity encoding (SENSE). With or without acceleration, the results also suggest that by taking advantage of the extra degrees of freedom inherent in a parallel system, parallel excitation provides better management of RF power deposition while facilitating the faithful production of desired excitation profiles. Sample accelerated and specific absorption rate (SAR)-reduced excitation pulses were designed in this study, and evaluated in experiments.     

Contrast-enhanced timing robust acquisition order with a preparation of the longitudinal signal component (CENTRA plus) for 3D contrast-enhanced abdominal imaging

PURPOSE: To investigate a new image acquisition method that enables an accurate hepatic arterial phase definition and the visualization of contrast agent uptake processes in abdominal organs like liver, spleen, and pancreas. MATERIALS AND METHODS: A 3D turbo gradient echo method where a fat suppression prepulse is followed by the acquisition of several profiles was combined with an elliptical centric k-space ordering technique and 3D dynamic elliptical centric keyhole. The new k-space ordering method (CENTRA+) was validated experimentally. In an initial clinical evaluation phase the method was employed in five patients to assess the accuracy of the hepatic arterial phase definition and the visualization of the contrast uptake processes in dynamic scanning in abdominal organs like liver, spleen, and pancreas. RESULTS: In total, five patients were evaluated using the new k-space order. Our initial results indicate that the new k-space order allows consistent capture of the hepatic arterial phase. In dynamic scanning the extreme short temporal resolution obtained with 3D elliptical centric keyhole enables contrast enhancement to be followed in organs with fast contrast uptake characteristics. CONCLUSION: The elliptical centric nature of the new image acquisition method effectively allows capture of the contrast enhancement processes with good fat suppression.     

Double average parallel steady-state free precession imaging: optimized eddy current and transient oscillation compensation.

Balanced steady-state free precession (SSFP) imaging is sensitive to off-resonance effects, which can lead to considerable artifacts during a transient phase following magnetization preparation or steady-state interruption. In addition, nonlinear k-space encoding is required if contrast-relevant k-space regions need to be acquired at specific delays following magnetization preparation or for transient artifact reduction in cardiac-gated k-space segmented CINE imaging. Such trajectories are problematic for balanced SSFP imaging due to nonconstant eddy current effects and resulting disruption of the steady state.In this work, a novel acquisition strategy for balanced SSFP imaging is presented that utilizes scan time reduction by parallel imaging for optimized "double average" eddy current compensation and artifact reduction during the transient phase following steady-state storage and magnetization preparation. Double average parallel SSFP imaging was applied to k-space segmented CINE SSFP tagging as well as nongated centrically encoded SSFP imaging. Phantom and human studies exhibit substantial reduction in steady-state storage and eddy current artifacts while maintaining spatial resolution, signal-to-noise ratio, and similar total scan time of a standard SSFP acquisition. The proposed technique can easily be extended to other acquisition schemes that would benefit from nonlinear reordering schemes and/or rely on interruption of the balanced SSFP steady state.     

T1rho imaging using magnetization-prepared projection encoding (MaPPE).

T1rho contrast weighting using a magnetization-prepared projection encoding (MaPPE) pulse sequence was investigated. Fast radial imaging was implemented by applying magnetization preparation pulses, each followed by multiple RF alpha pulses encoding radial trajectories of k-space. Acquiring multiple views per preparatory pulse imposes view-to-view variation; the resultant distortion of the point-spread function is examined. The issue of maximizing signal while preserving the intended contrast weighting is addressed. Under modification of repetition time and flip angle (alpha), three distinct behavior regimes of the sequence are identified. The utility of the pulse sequence as a quantitative relaxation measurement tool is also examined by comparing imaging and spectroscopy experiments. A mouse was imaged in vitro to demonstrate the viability of application to MR histology. These images exhibit the utility of spinlocking and projection encoding as an aftemative contrast source to both T2-weighted MaPPE images and conventional T2-weighted spin-echo images.     

Strategies to improve contrast in turboFLASH imaging: reordered phase encoding and k-space segmentation.

TurboFLASH (fast low-angle shot) sequences enable the acquisition of an image in a fraction of a second. However, unique to T1-weighted ultrafast imaging, the magnetization variation during image acquisition can produce artifacts along the phase-encoding direction. In this study, the signal behavior and nature of these artifacts were analyzed with various acquisition schemes to improve image contrast. The magnetization variation during image acquisition and its filtering effect on the image were simulated for three different approaches to T1-weighted turboFLASH imaging: standard turboFLASH with (a) monotonically ascending phase-encoding steps, (b) reordered phase encoding, and (c) k-space segmentation. Each of the modified data acquisition schemes has advantages. However, for subsecond imaging, reordered phase encoding produced improved image contrast over that of standard turboFLASH, and segmented k-space imaging gave superior tissue contrast compared with that of both standard and reordered turboFLASH, with imaging time that permits breath-hold studies.     

Description of parallel imaging in MRI using multiple coils.

A general formulation for parallel imaging using multiple coils is derived from the Fourier transform of coil sensitivity functions. This formulation provides a unified account for developed parallel imaging techniques such as subencoding, simultaneous acquisition of spatial harmonics (SMASH), and sensitivity encoding (SENSE), and indicates a guideline for coil configuration and k-space sampling in parallel imaging. The views that can be acquired simultaneously in parallel imaging have to be contained in the spatial frequency band of coil sensitivity functions.     

K-space in the clinic

Magnetic resonance imaging (MRI) sequences are characterized by both radio frequency (RF) pulses and time-varying gradient magnetic fields. The RF pulses manipulate the alignment of the resonant nuclei and thereby generate a measurable signal. The gradient fields spatially encode the signals so that those arising from one location in an excited slice of tissue may be distinguished from those arising in another location. These signals are collected and mapped into an array called k-space that represents the spatial frequency content of the imaged object. Spatial frequencies indicate how rapidly an image feature changes over a given distance. It is the action of the gradient fields that determines where in the k-space array each data point is located, with the order in which k-space points are acquired being described by the k-space trajectory. How signals are mapped into k-space determines much of the spatial, temporal, and contrast resolution of the resulting images and scan duration. The objective of this article is to provide an understanding of k-space as is needed to better understand basic research in MRI and to make well-informed decisions about clinical protocols. Four major classes of trajectories-echo planar imaging (EPI), standard (non-EPI) rectilinear, radial, and spiral-are explained. Parallel imaging techniques SMASH (simultaneous acquisition of spatial harmonics) and SENSE (sensitivity encoding) are also described.