Kuan J. Lee

Compressed sensing in hyperpolarized 3He Lung MRI

In this work, the application of compressed sensing techniques to the acquisition and reconstruction of hyperpolarized 3He lung MR images was investigated. The sparsity of 3He lung images in the wavelet domain was investigated through simulations based on fully sampled Cartesian two‐dimensional and three‐dimensional 3He lung ventilation images, and the k‐spaces of 2D and 3D images were undersampled randomly and reconstructed by minimizing the L1 norm. The simulation results show that temporal resolution can be readily improved by a factor of 2 for two‐dimensional and 4 to 5 for three‐dimensional ventilation imaging with 3He with the levels of signal to noise ratio (SNR) (∼19) typically obtained. The feasibility of producing accurate functional apparent diffusion coefficient (ADC) maps from undersampled data acquired with fewer radiofrequency pulses was also demonstrated, with the preservation of quantitative information (mean ADCcs ∼ mean ADCfull ∼ 0.16 cm2 sec−1). Prospective acquisition of 2‐fold undersampled two‐dimensional 3He images with a compressed sensing k‐space pattern was then demonstrated in a healthy volunteer, and the results were compared to the equivalent fully sampled images (SNRcs = 34, SNRfull = 19). Magn Reson Med 63:1059–1069, 2010.

Use of forward projection to correct patient motion during SPECT imaging

We have developed an iterative method to correct axial and tangential patient motion occurring during tomographic acquisition. The method uses axial images reconstructed from the uncorrected projection images, which are then forward projected to form a basis for registering the original planar images and, in the process, directly seeks to establish a consistent data set. Our method can be applied to all SPECT scans including myocardial and brain SPECT. We demonstrate that the method is capable of detecting and quantitatively correcting for complex motion in both axial and tangential directions. Results from phantom experiments show excellent resolution and contrast recovery after simulated movement in both the axial and tangential directions and initial results with clinical data sets are encouraging.

Image-based EPI ghost correction using an algorithm based on projection onto convex sets (POCS)

This work describes the use of a method, based on the projection onto convex sets (POCS) algorithm, for reduction of the N/2 ghost in echo‐planar imaging (EPI). In this method, ghosts outside the parent image are set to zero and a model k‐space is obtained from the Fourier transform (FT) of the resulting image. The zeroth‐ and first‐order phase corrections for each line of the original k‐space are estimated by comparison with the corresponding line in the model k‐space. To overcome problems of phase wrapping, the first‐order phase corrections for the lines of the original k‐space are estimated by registration with the corresponding lines in the model k‐space. It is shown that applying these corrections will result in a reduction of the ghost, and that iterating the process will result in a convergence towards an image in which the ghost is minimized. The method is tested on spin‐echo EPI data. The results show that the method is robust and remarkably effective, reducing the N/2 ghost to a level nearly comparable to that achieved with reference scans. Magn Reson Med 47:812–817, 2002.

Brain arteriovenous malformations: measurement of nidal volume using a combination of static and dynamic magnetic resonance angiography techniques

Arteriovenous malformations of the brain are complex vascular lesions that are an important cause of death and long-term disability. Currently, catheter angiography (CA) is the reference standard procedure for the diagnosis and follow-up of treated arteriovenous malformations (AVMs). This is an invasive procedure with potential risks. Magnetic resonance angiography (MRA) is commonly used in neurovascular imaging as a non-invasive alternative. Various MRA techniques have been used in the diagnosis and follow-up of AVMs but these have suffered from lack of temporal or spatial resolution. In this 60-patient study we describe the combination of two techniques: dynamic magnetic resonance digital subtraction angiography with a high temporal resolution, and a non-dynamic contrast-enhanced time-of-flight sequence with a high spatial resolution technique, in the assessment of AVM. The results showed an excellent correlation between MRA and CA measurement of both maximum linear dimension and AVM nidus volume.

Slice profile effects in 2D slice-selective MRI of hyperpolarized nuclei.

This work explores slice profile effects in 2D slice-selective gradient-echo MRI of hyperpolarized nuclei. Two different sequences were investigated: a Spoiled Gradient Echo sequence with variable flip angle (SPGR-VFA) and a balanced Steady-State Free Precession (SSFP) sequence. It is shown that in SPGR-VFA the distribution of flip angles across the slice present in any realistically shaped radiofrequency (RF) pulse leads to large excess signal from the slice edges in later RF views, which results in an undesired non-constant total transverse magnetization, potentially exceeding the initial value by almost 300% for the last RF pulse. A method to reduce this unwanted effect is demonstrated, based on dynamic scaling of the slice selection gradient. SSFP sequences with small to moderate flip angles (<40 degrees ) are also shown to preserve the slice profile better than the most commonly used SPGR sequence with constant flip angle (SPGR-CFA). For higher flip angles, the slice profile in SSFP evolves in a manner similar to SPGR-CFA, with depletion of polarization in the center of the slice.

Automated gamma knife radiosurgery treatment planning with image registration, data-mining, and Nelder-Mead simplex optimization

Gamma knife treatments are usually planned manually, requiring much expertise and time. We describe a new, fully automatic method of treatment planning. The treatment volume to be planned is first compared with a database of past treatments to find volumes closely matching in size and shape. The treatment parameters of the closest matches are used as starting points for the new treatment plan. Further optimization is performed with the Nelder-Mead simplex method: the coordinates and weight of the isocenters are allowed to vary until a maximally conformal plan specific to the new treatment volume is found. The method was tested on a randomly selected set of 10 acoustic neuromas and 10 meningiomas. Typically, matching a new volume took under 30 seconds. The time for simplex optimization, on a 3 GHz Xeon processor, ranged from under a minute for small volumes (<1000 cubic mm, 2-3 isocenters), to several tens of hours for large volumes (>30,000 cubic mm, >20 isocenters). In 8/10 acoustic neuromas and 8/10 meningiomas, the automatic method found plans with conformation number equal or better than that of the manual plan. In 4/10 acoustic neuromas and 5/10 meningiomas, both overtreatment and undertreatment ratios were equal or better in automated plans. In conclusion, data-mining of past treatments can be used to derive starting parameters for treatment planning. These parameters can then be computer optimized to give good plans automatically.

Parallel imaging of hyperpolarized helium‐3 with simultaneous slice excitation

Hyperpolarized (HP) gas imaging of the lungs is an ideal potential application for parallel imaging. This is due to the fact that there is limited scan time (breath hold of 20 s) and limited non‐renewable polarization. Reduced phase encode parallel imaging is demanding on hardware in that it requires multiple receivers. In this work, simultaneous parallel acquisition of hyperpolarized (HP) 3He images from multiple slices was demonstrated in phantoms and in vivo using a simultaneous slice excitation method, at a field strength of 1.5 T. The pulse sequence allows simultaneous acquisition of n slices per RF excitation, thus reducing the number of RF pulses needed to fully cover a given volume with multi‐slicing. Unlike conventional parallel imaging, this method does not require prior reference scan information, which would consume some of the finite longitudinal polarization in lung ventilation studies with HP gas. Magn Reson Med, 2006.

B1AC‐MAMBA: B1 array combined with multiple‐acquisition micro B0 array parallel magnetic resonance imaging

The combination of an in‐plane B1 sensitivity encoding (SENSE) technique with a simultaneous multiple‐slice B0 field step technique (multiple‐acquisition micro B0 array (MAMBA)) has produced high scan time reduction factors (R ≤ 8). In this study, two slices were acquired simultaneously in combination with ×2 and ×4 SENSE in‐plane encoding using a MAMBA stepped B0 field coil inside a four‐channel phased‐array coil system. Experiments were performed on a 1.5 T Infinion system (Philips Medical Systems, Cleveland, OH). The signal‐to‐noise ratio (SNR) was reduced with higher R factors, as was expected from the reduced number of acquisitions used to create the unaliased images. The combination of SENSE and MAMBA offers great promise for reducing scan times through parallel acquisition while at the same time reducing the number of RF channels required by a factor equal to the number of field steps employed. The B1 array combined with MAMBA (B1AC‐MAMBA) technique is applicable when the length of an object is much greater than its diameter, as in scanning limbs or in whole‐body screening for disease. Magn Reson Med 49:1196–1200, 2003.

Interleaved pulsed MAMBA: A new parallel slice imaging method

A method of acquiring slices in parallel is described which uses interleaved sets of pulsed B0 field coils to generate discrete regions of uniform field within the main magnetic field known as interleaved MAMBA (multiple acquisition micro B0 array). Simulations of a number of coil designs were performed using the Biot‐Savart law. A six‐step coil was built and interfaced to a 0.17 T Niche MRI system and the field steps measured using an imaging technique. Measured field steps were in good agreement with the values predicted by simulation. The coil design was then scaled up by a factor of three, interfaced to a 1.5 T whole‐body MRI system, and scans of the hands and arms of volunteers were acquired from up to four field steps using standard spin and gradient echo sequences. Images were also acquired simultaneously from two field steps with no frequency encode aliasing and one excitation. The one‐dimensional interleaved pulsed MAMBA step field technique shows great promise for enabling many slices to be acquired simultaneously along the axis of the coil for rapid volumetric studies without the need for multiple shot Hadamard encoding. Extension of interleaved coil design to two or three dimensions is feasible, which could provide full spatial coverage combined with ultra‐rapid data acquisition. Magn Reson Med 48:1043–1050, 2002.

Fast two-dimensional MR imaging by multiple acquisition with micro B0 array (MAMBA)

A new method for acquiring MR data in two dimensions is described. This is achieved by combining coils into an array so as to produce a unique local magnetic field within each coil when placed within the B(0) field of a standard MRI scanner. In this way each known location of a coil is associated with a unique resonant frequency. Each coil now represents a location of a pixel in the plane, and after Fourier transformation of the signal the resulting frequency spectrum gives immediately the spin distribution in the plane of the array. In effect, the two-dimensional spatial distribution is frequency encoded without the use of switched gradients or phase encoding. As only static fields are used, this technique offers the potential of fast imaging. Furthermore, signals from different locations would also be inherently time-registered. Initial experiments to demonstrate the principle are described, using a square array of 5 by 5 coils. The currents in the coils were determined by using a genetic algorithm. Echoes from pellet phantoms placed in the array were acquired using standard spin-echo sequences with gradients switched off. The results are promising, with the spectra showing generally good resolution between peaks, enabling localisation in up to half the pixels. Technical difficulties are discussed and possible applications are outlined.