Biomedical subjects

Jeffrey Tsao

Publications and source records attributed to Jeffrey Tsao.

16 recordsLinked to original sources

PubMed2006

k-t BLAST reconstruction from non-Cartesian k-t space sampling.

Current implementations of k-t Broad-use Linear Acqusition Speed-up Technique (BLAST) require the sampling in k-t space to conform to a lattice. To permit the use of k-t BLAST with non-Cartesian sampling, an iterative reconstruction approach is proposed in this work. This method, which is based on the conjugate gradient (CG) method and gridding reconstruction principles, can efficiently handle data that are sampled along non-Cartesian trajectories in k-t space. The approach is demonstrated on prospectively gated radial and retrospectively gated Cartesian imaging. Compared to a sliding window (SW) reconstruction, the resulting image series exhibit lower artifact levels and improved temporal fidelity. The proposed approach thus allows investigators to combine the specific advantages of non-Cartesian imaging or retrospective gating with the acceleration provided by k-t BLAST.

Artifacts↗

PubMed2006

Lattice permutation for reducing motion artifacts in radial and spiral dynamic imaging.

Radial and spiral trajectories exhibit favorable characteristics for dynamic imaging. Nevertheless, changes in image contents during acquisition lead to inconsistencies in the k-space data, which are manifested as streaks or spiral artifacts, respectively. This work proposes the concept of lattice permutation to reorder the data segments for artifact suppression. This acts to reshuffle the alias pattern along the temporal frequency axis. The proposed approach is well suited to sliding window reconstruction, although more sophisticated methods are also possible. For typical image series where the signal energies are concentrated in the low temporal frequencies, the permutation displaces most of the aliased signals from the low temporal frequencies to the high temporal frequencies, where they are attenuated by sliding window reconstruction, while the signals in the low temporal frequencies are mostly contaminated by aliasing from the much weaker signals in the higher temporal frequencies. This results in considerably reduced artifacts without any increase in scan time. In practice, lattice permutation achieves similar artifact suppression as the bit-reversed order, but with a less stringent restriction on the number of segments. At the same time, it provides a more powerful approach to controlling the alias pattern exactly. Results from real-time cardiac imaging are demonstrated.

Artifacts↗

PubMed2006

Minimum-norm reconstruction for sensitivity-encoded magnetic resonance spectroscopic imaging.

In this work we propose minimum-norm reconstruction as a means to enhance the spatial response behavior in parallel spectroscopic MRI. By directly optimizing the shape of the spatial response function (SRF), the new method accounts for coil sensitivity variation across individual voxels and their side lobes. In this fashion, it mitigates the signal contamination and side-lobe aliasing, to which previous techniques are susceptible at low resolution. Although the computational burden is higher, minimum-norm reconstruction is shown to be feasible using an iterative algorithm. Benefits in terms of SRF shape and artifact suppression are demonstrated.

Algorithms↗

PubMed2005

Linear response equilibrium.

A new periodic pulse sequence employing weak excitation is presented. This type of sequence drives the system into a steady-state with periodic time evolution from which the data can be reconstructed to a spectrum. It is demonstrated that the frequency response of such a sequence can be analyzed using perturbation methods and linear system analysis. A mathematical framework is proposed allowing the frequency response to be tailored by weighting a periodic flip function. The weak excitation level used implies very low specific absorption rates while generating a highly frequency selective signal in the order of 1/T2 with signal strengths comparable to those obtainable with conventional large flip angle balanced steady-state free precession techniques. The concept is illustrated with phantom experiments and in vivo feasibility of water fat separation is shown on human knee images.

Adipose Tissue↗

PubMed2005

Optimizing spatiotemporal sampling for k-t BLAST and k-t SENSE: application to high-resolution real-time cardiac steady-state free precession.

In k-t BLAST and k-t SENSE, data acquisition is accelerated by sparsely sampling k-space over time. This undersampling in k-t space causes the object signals to be convolved with a point spread function in x-f space (x = spatial position, f = temporal frequency). The resulting aliasing is resolved by exploiting spatiotemporal correlations within the data. In general, reconstruction accuracy can be improved by controlling the k-t sampling pattern to minimize signal overlap in x-f space. In this work, we describe an approach to obtain generally favorable patterns for typical image series without specific knowledge of the image series itself. These optimized sampling patterns were applied to free-breathing, untriggered (i.e., real-time) cardiac imaging with steady-state free precession (SSFP). Eddy-current artifacts, which are otherwise increased drastically in SSFP by the undersampling, were minimized using alternating k-space sweeps. With the synergistic combination of the k-t approach with optimized sampling and SSFP with alternating k-space sweeps, it was possible to achieve a high signal-to-noise ratio, high contrast, and high spatiotemporal resolutions, while achieving substantial immunity against eddy currents. Cardiac images are shown, demonstrating excellent image quality and an in-plane resolution of approximately 2.0 mm at >25 frames/s, using one or more receiver coils.

Artifacts↗

PubMed2005

Accelerating cine phase-contrast flow measurements using k-t BLAST and k-t SENSE.

Conventional phase-contrast velocity mapping in the ascending aorta was combined with k-t BLAST and k-t SENSE. Up to 5.3-fold net acceleration was achieved, enabling single breath-hold acquisitions. A standard phase-contrast (PC) sequence with interleaved acquisition of the velocity-encoded segments was modified to collect data in 2 stages, a high-resolution under sampled and a low-resolution fully sampled training stage. In addition, a modification of the k-t reconstruction strategy was tested. This strategy, denoted as "plug-in," incorporates data acquired in the training stage into the final reconstruction for improved data consistency, similar to conventional keyhole. "k-t SENSE plug-in" was found to provide best image quality and most accurate flow quantification. For this strategy, at least 10 training profiles are required to yield accurate stroke volumes (relative deviation <5%) and good image quality. In vivo 2D cine velocity mapping was performed in 6 healthy volunteers with 30-32 cardiac phases (spatial resolution 1.3 x 1.3 x 8-10 mm(3), temporal resolution of 18-38 ms), yielding relative stroke volumes of 106 +/- 18% (mean +/- 2*SD) and 112 +/- 15% for 3.8 x and 5.3 x net accelerations, respectively. In summary, k-t BLAST and k-t SENSE are promising approaches that permit significant scan-time reduction in PC velocity mapping, thus making high-resolution breath-held flow quantification possible.

Algorithms↗

PubMed2004

Accelerated dynamic Fourier velocity encoding by exploiting velocity-spatio-temporal correlations.

OBJECTIVE: To describe how the information content in a Fourier velocity encoding (FVE) scan can be transformed into a very sparse representation and to develop a method that exploits the compactness of the data to significantly accelerate the acquisition. MATERIALS AND METHODS: For validation, fully sampled FVE datasets were acquired in phantom and in vivo experiments. Fivefold and eightfold acceleration was simulated by using only one fifth or one eighth of the data for reconstruction in the proposed method based on the k-t BLAST framework. Reconstructed images were compared quantitatively to those from the fully sampled data. RESULTS: Velocity spectra in the accelerated datasets were comparable to the spectra from fully sampled datasets. The detected peak velocities remained accurate even at eightfold acceleration, and the overall shape of the spectra was well preserved. Slight temporal smoothing was seen in the accelerated datasets. CONCLUSION: A novel technique for accelerating time-resolved FVE scan is presented. It is possible to accelerate FVE to acquisition speeds comparable to a standard time-resolved phase-contrast scan.

Algorithms↗

PubMed2004

Peak-combination HARP: a method to correct for phase errors in HARP.

PURPOSE: To introduce a method to correct phase errors (e.g., from B0 inhomogeneity) in tagging images, which may affect harmonic phase (HARP) evaluation. MATERIALS AND METHODS: The phase images corresponding to the negative and positive harmonic peaks in k-space are combined before HARP evaluation to eliminate any spurious phase. To validate in vivo, two complementary spatial modulation of magnetization (CSPAMM) data sets were collected for each volunteer and evaluated with conventional HARP, using either the positive or the negative harmonic peak, and with peak-combination HARP. RESULTS: Elimination of phase distortion by peak combination was observed in vitro and in vivo. Improved reproducibility of motion parameters was found with peak-combination HARP. CONCLUSION: With peak-combination HARP, reproducibility of contractile parameters is improved, and consequently, the number of subjects needed to detect statistically significant changes in contractile function can be reduced to one third compared to conventional HARP evaluation.

Adult↗

PubMed2004

Accelerating cardiac cine 3D imaging using k-t BLAST.

By exploiting spatiotemporal correlations in cardiac acquisitions using k-t BLAST, gated cine 3D acquisitions of the heart were accelerated by a net factor of 4.3, making single breathhold acquisitions possible. Sparse sampling of k-t space along a sheared grid pattern was implemented into a cine 3D SSFP sequence. The acquisition of low-resolution training data, which was required to resolve aliasing in the k-t BLAST method, was either interleaved into the sampling process or obtained in a separate prescan to allow for shorter breathhold durations in patients with heart disease. Volumetric datasets covering the heart with 20 slices at a spatial resolution of 2 x 2 x 5 mm3 were recorded with 20 cardiac phases in a total breathhold duration of 25-27 sec, or 18 sec if partial Fourier sampling was additionally employed. The feasibility of the method was demonstrated on healthy volunteers and on patients. The comparison of endocardial area derived from single slices of the 3D dataset with values extracted from separate single-slice acquisitions showed no significant differences. By shortening the acquisition substantially, k-t BLAST may greatly facilitate volumetric imaging of the heart for evaluation of regional wall motion and the assessment of ventricular volume and ejection fraction.

Computer Simulation↗

PubMed2004

On the influence of training data quality in k-t BLAST reconstruction.

This work investigated how the quality of prior information (i.e., data acquired during the training stage) influences k-t BLAST reconstruction. The impact of several factors, such as the amount of training data, the presence of spatial misregistration in the training data, and the effects of filtering, was investigated with simulations and in vivo data. It is shown that k-t BLAST outperforms sliding window reconstruction, even with very limited training data. By increasing the amount of training data, reconstruction error continues to decrease, albeit by a diminishing amount. However, an increased amount of training data also increases susceptibility to misregistration of the training data. Filtering of the training data with the goal of reducing truncation artifacts had only minor impact on reconstruction errors. Considering the balance among obtaining the most benefit from the training data, minimizing susceptibility to misregistration, and keeping data acquisition to a minimum, it is concluded that in cardiac imaging the training datasets should be limited to 10-20 profiles in k-space for a typical field of view. The training data may be acquired in a separate breathhold without much penalty, if care is taken to minimize misregistration, such as with a navigator.

Adult↗

PubMed2004

Equi-ripple design of quadratic-phase RF pulses.

An improved strategy for the design of quadratic-phase RF pulses with high selectivity and broad bandwidths using the Shinnar-Le Roux (SLR) transformation is proposed. Unlike previous implementations, the required quadratic-phase finite impulse response (FIR) filters are generated using the complex Remez exchange algorithm, which ensures an equi-ripple deviation from the ideal response function. It is argued analytically that quadratic-phase pulses are near-optimal in terms of minimising the B1-amplitude for a given bandwidth and flip angle. Furthermore, several parameter relations are derived, providing practical design guidelines. The effectiveness of the proposed design method is demonstrated by examples of excitation and saturation pulses applied in vitro and in vivo.

Algorithms↗

PubMed2004

Reduced data acquisition methods in cardiac imaging.

In cardiac imaging, acquisition speed is of primary importance. While improved performance has mainly been achieved through improvements in gradient hardware in the past, further developments along this direction are limited due to physiological constraints such as the risk of peripheral nerve stimulation. With the introduction of parallel imaging, alternative means for increasing acquisition speed have become available. Using information from multiple receiver coils, images can be reconstructed from a sparsely sampled set of data. In practice, parallel imaging allows for 2- to 3-fold acceleration of the imaging process in typical cardiac applications. Further increases in acquisition speed are, however, difficult to achieve for current clinical field strengths and typical field of views. To address the limited gain in acquisition speed achievable with parallel imaging, a new set of methods has been proposed to take into account the similarity of image information at different time points during a dynamic series. Using these methods, 5- to 8-fold acceleration can be achieved in cardiac imaging. It is the purpose of this paper to review cardiac applications of reduced data acquisition methods with focus on parallel imaging and the recently developed k-t BLAST and k-t SENSE techniques.

Heart↗

PubMed2003

Parallel spectroscopic imaging with spin-echo trains.

A reduction in scan time in spectroscopic imaging (SI) can be achieved by both fast and reduced k-space sampling. This work presents an ultrafast SI technique that combines the two approaches. The synergy of multiple spin-echo (MSE) acquisition and sensitivity encoding (SENSE) enables high-resolution SI to be performed within a clinically acceptable scan time. MSE-SENSE-SI with echo train lengths ranging from one to four echoes is evaluated with respect to SNR and spatial response function by means of in vitro experiments. It is shown that acquiring two spin-echoes (SEs) per acquisition yields a good practical trade-off among scan time, SNR, and spatial response. The clinical feasibility of the technique is demonstrated in a patient with an astrocytoma, and SI data are obtained with an image matrix of 24 x 24 in just over 2 min.

Astrocytoma↗

PubMed2003

k-t BLAST and k-t SENSE: dynamic MRI with high frame rate exploiting spatiotemporal correlations.

Dynamic images of natural objects exhibit significant correlations in k-space and time. Thus, it is feasible to acquire only a reduced amount of data and recover the missing portion afterwards. This leads to an improved temporal resolution, or an improved spatial resolution for a given amount of acquisition. Based on this approach, two methods were developed to significantly improve the performance of dynamic imaging, named k-t BLAST (Broad-use Linear Acquisition Speed-up Technique) and k-t SENSE (SENSitivity Encoding) for use with a single or multiple receiver coils, respectively. Signal correlations were learned from a small set of training data and the missing data were recovered using all available information in a consistent and integral manner. The general theory of k-t BLAST and k-t SENSE is applicable to arbitrary k-space trajectories, time-varying coil sensitivities, and under- and overdetermined reconstruction problems. Examples from ungated cardiac imaging demonstrate a 4-fold acceleration (voxel size 2.42 x 2.52 mm(2), 38.4 fps) with either one or six receiver coils. k-t BLAST and k-t SENSE are applicable to many areas, especially those exhibiting quasiperiodic motion, such as imaging of the heart, the lungs, the abdomen, and the brain under periodic stimulation.

Heart↗

PubMed2003

Interpolation artifacts in multimodality image registration based on maximization of mutual information.

Mutual information (MI) is an increasingly popular match metric for multimodality image registration. However, its value is affected by interpolation, which may limit registration accuracy. The purpose of this study was to characterize the artifacts from eight interpolators and to investigate efficient strategies to overcome these artifacts. The interpolators were: 1) nearest neighbor; 2) linear; 3) cubic Catmull-Rom; 4) Hamming-windowed sinc; 5) partial volume; 6) NN with jittered sampling (JIT); 7) NN with histogram blurring (BLUR); and 8) NN with JIT and BLUR. The impact of interpolation on MI was evaluated in two dimensions over different translational and rotational misregistration. Interpolation caused spurious fluctuations in MI whenever the voxel grids had coinciding periodicities and were nearly aligned. The artifacts did not lessen by using intensity interpolators with wider support (e.g., cubic Catmull-Rom, Hamming-windowed sinc). PV could lead to either arch artifacts or inverted-arch artifacts, depending on the relative voxel sizes. Several strategies reduced artifacts and improved registration robustness: JIT, BLUR, avoiding an extreme number of intensity bins, and resampling the images in a rotated orientation with different relative voxel sizes (e.g., pi/3). These findings also apply to related methods, including normalized MI, joint entropy, and Hill's third moment.

Algorithms↗

PubMed2002

On the UNFOLD method.

A graphical formalism is presented, showing that "UNaliasing by Fourier-encoding the Overlaps Using the temporaL Dimension" (UNFOLD) is equivalent to sampling k-t-space in a sheared grid pattern. Discrete regular sampling in k-t-space leads to periodic replication of the support region in x-f-space. Thus, the maximum acceleration achievable by UNFOLD is equivalent to the maximum packing of support regions in x-f-space. When the support region is separable along the x and f axes, the reconstruction can be performed separately for each k. UNFOLD can be combined with SiMultaneous Acquisition of Spatial Harmonics (SMASH) to further accelerate acquisition. However, a straightforward combination of the methods has been shown to result in a size restriction, which limits the portion of the field of view (FOV) with a larger temporal bandwidth to only a quarter of the FOV. Two solutions are presented to overcome this restriction.

Algorithms↗