Efficient Algorithms for Parallel Excitation and Parallel Imaging with Large Arrays

During the past two decades, techniques and devices were developed to transmit and receive signals with a phased array instead of a single coil in the MRI (Magnetic Resonance Imaging) system. The two techniques to simultaneously transmit and receive RF signals using phased arrays are called parallel excitation (pTx) and parallel imaging (PI), respectively. These two techniques lead to shorter transmit pulses for higher imaging quality and faster data acquisition correspondingly. This dissertation focuses on improving the efficiency of the pTx pulse design and the PI reconstruction in MRI. Both PI and pTx benefit from the increased number of elements of the array. However, efficiency concerns may arise which include: (1) In PI, the computation cost of the reconstructions and the achievable acceleration factors and (2) in pTx, the pulse design speed and memory cost. The work presented in this dissertation addresses these issues. First, a correlation based channel reduction algorithm is developed to reduce the computation cost of PI reconstruction. In conventional k-domain methods, the individual channel data is reconstructed via linear interpolation of the neighbourhood data from all channels. In this proposed algorithm, we choose only a subset of the channels based on the spatial correlation. The results have shown that the computation cost can be significantly reduced with similar or higher reconstruction accuracy. Then, a new parallel imaging method named parallel imaging using localized receive arrays with Sinc interpolation(PILARS) is proposed to improve the actual acceleration factor and to reduce the computation cost. It employs the local support of individual coils and pre-determines the magnitude of the reconstruction coefficients. Thus, it requires much less auto-calibration signals (ACS) data and achieves higher acceleration factors. The results show that this method can increase the acceleration factor and the reconstruction speed while achieving the same level of reconstruction quality. Finally, a fast pTx pulse design method is proposed to accelerate the design speed. This method is based on the spatial domain pulse design method and can be used to accelerate similar methods. We substitute the two computational expensive matrix- vector multiplications in the conjugate gradient (CG) solver with gridding and fast Fourier transform (FFT). Theoretical and simulation results have shown that the design speed can be improved by 10 times. Meanwhile, the memory cost is reduced by 103 times. This breaks the memory burden of implementing pulse designs on GPU which enables further accelerations

Acquisition strategies for fat/water separated MRI

This thesis focuses on new ways to more efficiently acquire the signal for fat/water separated MRI, also known as Dixon methods. In paper I, the concept of dual bandwidths was introduced to improve SNR efficiency by removing dead times in a spin echo PROPELLER sequence. By correcting for the displacement of fat, we were able to improve the motion correction. This required additional considerations during reconstruction in order to avoid noise amplification, which was solved with a noise-whitening Tikhonov regularization. Paper II explores the combination of fat/water separation in k-space with partially acquired data, i.e. partial Fourier sampling. With reduced sampling coverage comes the ability of increased spatial resolution, which is often limited in fat/water imaging, particularly in gradient echo sequences. A modified POCS routine was also developed with real-valued estimates, exploiting Hermitian symmetry to improve the inverse problem conditioning in the fully sampled region. A single-TR dual-bandwidth RARE (fast/turbo spin echo) sequence without dead times was developed in Paper III, which uses partial Fourier sampling with late and early echoes to improve the chemical shift encoding. The proposed sequence can acquire images with 0.5 mm in-plane resolution without dead times, with image quality exceeding current state-of-the-art techniques. An automated selection of gradient waveforms based on Cramér-Rao bounds was implemented on the scanner. In Paper IV, the dual-bandwidth concept was generalized to continuous bandwidths. Instead of the conventional shift of a trapezoidal readout gradient, we describe a new method of encoding chemical shift by using asymmetric readout waveforms. Asymmetric readouts were implemented in a RARE sequence to completely remove dead times from multi-TR acquisitions, with typical scan time reductions of 25 %. The developed methods enable fat/water imaging with reduced scan times and increased spatial resolution, which has previously limited their use

Rapid 3D Phase Contrast Magnetic Resonance Angiography through High-Moment Velocity Encoding and 3D Parallel Imaging

abstract: Phase contrast magnetic resonance angiography (PCMRA) is a non-invasive imaging modality that is capable of producing quantitative vascular flow velocity information. The encoding of velocity information can significantly increase the imaging acquisition and reconstruction durations associated with this technique. The purpose of this work is to provide mechanisms for reducing the scan time of a 3D phase contrast exam, so that hemodynamic velocity data may be acquired robustly and with a high sensitivity. The methods developed in this work focus on the reduction of scan duration and reconstruction computation of a neurovascular PCMRA exam. The reductions in scan duration are made through a combination of advances in imaging and velocity encoding methods. The imaging improvements are explored using rapid 3D imaging techniques such as spiral projection imaging (SPI), Fermat looped orthogonally encoded trajectories (FLORET), stack of spirals and stack of cones trajectories. Scan durations are also shortened through the use and development of a novel parallel imaging technique called Pretty Easy Parallel Imaging (PEPI). Improvements in the computational efficiency of PEPI and in general MRI reconstruction are made in the area of sample density estimation and correction of 3D trajectories. A new method of velocity encoding is demonstrated to provide more efficient signal to noise ratio (SNR) gains than current state of the art methods. The proposed velocity encoding achieves improved SNR through the use of high gradient moments and by resolving phase aliasing through the use measurement geometry and non-linear constraints.Dissertation/ThesisPh.D. Bioengineering 201

Doctor of Philosophy

dissertationDynamic contrast enhanced magnetic resonance imaging (DCE-MRI) is a powerful tool to detect cardiac diseases and tumors, and both spatial resolution and temporal resolution are important for disease detection. Sampling less in each time frame and applying sophisticated reconstruction methods to overcome image degradations is a common strategy in the literature. In this thesis, temporal TV constrained reconstruction that was successfully applied to DCE myocardial perfusion imaging by our group was extended to three-dimensional (3D) DCE breast and 3D myocardial perfusion imaging, and the extension includes different forms of constraint terms and various sampling patterns. We also explored some other popular reconstruction algorithms from a theoretical level and showed that they can be included in a unified framework. Current 3D Cartesian DCE breast tumor imaging is limited in spatiotemporal resolution as high temporal resolution is desired to track the contrast enhancement curves, and high spatial resolution is desired to discern tumor morphology. Here temporal TV constrained reconstruction was extended and different forms of temporal TV constraints were compared on 3D Cartesian DCE breast tumor data with simulated undersampling. Kinetic parameters analysis was used to validate the methods

Aggregated motion estimation for real-time MRI reconstruction

Real-time magnetic resonance imaging (MRI) methods generally shorten the measuring time by acquiring less data than needed according to the sampling theorem. In order to obtain a proper image from such undersampled data, the reconstruction is commonly defined as the solution of an inverse problem, which is regularized by a priori assumptions about the object. While practical realizations have hitherto been surprisingly successful, strong assumptions about the continuity of image features may affect the temporal fidelity of the estimated images. Here we propose a novel approach for the reconstruction of serial real-time MRI data which integrates the deformations between nearby frames into the data consistency term. The method is not required to be affine or rigid and does not need additional measurements. Moreover, it handles multi-channel MRI data by simultaneously determining the image and its coil sensitivity profiles in a nonlinear formulation which also adapts to non-Cartesian (e.g., radial) sampling schemes. Experimental results of a motion phantom with controlled speed and in vivo measurements of rapid tongue movements demonstrate image improvements in preserving temporal fidelity and removing residual artifacts.Comment: This is a preliminary technical report. A polished version is published by Magnetic Resonance in Medicine. Magnetic Resonance in Medicine 201

Doctor of Philosophy

dissertationThis dissertation presents original research that improves the ability of magnetic resonance imaging (MRI) to measure temperature in aqueous tissue using the proton resonance frequency (PRF) shift and T1 measurements in fat tissue in order to monitor focused ultrasound (FUS) treatments. The inherent errors involved in measuring the longitudinal relaxation time T1 using the variable flip angle method with a two-dimensional (2D) acquisition are presented. The edges of the slice profile can contribute a significant amount of signal for large flip angles at steady state, which causes significant errors in the T1 estimate. Only a narrow range of flip angle combinations provided accurate T1 estimates. Respiration motion causes phase artifacts, which lead to errors when measuring temperature changes using the PRF method. A respiration correction method for 3D imaging temperature of the breast is presented. Free induction decay (FID) navigators were used to measure and correct phase offsets induced by respiration. The precision of PRF temperature measurements within the breast was improved by an average factor of 2.1 with final temperature precision of approximately 1 °C. Locating the position of the ultrasound focus in MR coordinates of an ultrasound transducer with multiple degrees of freedom can be difficult. A rapid method for predicting the position using 3 tracker coils with a special MRI pulse iv sequence is presented. The Euclidean transformation of the coil's current positions to their calibration positions was used to predict the current focus position. The focus position was predicted to within approximately 2.1 mm in less than 1 s. MRI typically has tradeoffs between imaging field of view and spatial and temporal resolution. A method for acquiring a large field of view with high spatial and temporal resolution is presented. This method used a multiecho pseudo-golden angle stack of stars imaging sequence to acquire the large field of view with high spatial resolution and k-space weighted image contrast (KWIC) to increase the temporal resolution. The pseudo-golden angle allowed for removal of artifacts introduced by the KWIC reconstruction algorithm. The multiple echoes allowed for high readout bandwidth to reduce blurring due to off resonance and chemical shift as well as provide separate water/fat images, estimates of the initial signal magnitude M(0), T2 * time constant, and combination of echo phases. The combined echo phases provided significant improvement to the PRF temperature precision, and ranged from ~0.3-1.0 °C within human breast. M(0) and T2 * values can possibly be used as a measure of temperature in fat