Needle artifact localization in 3T MR images.

This work explores an image-based approach for localizing needles during MRI-guided interventions, for the purpose of tracking and navigation. Susceptibility artifacts for several needles of varying thickness were imaged, in phantoms, using a 3 tesla MRI system, under a variety of conditions. The relationship between the true needle positions and the locations of artifacts within the images, determined both by manual and automatic segmentation methods, have been quantified and are presented here.

Signal-to-noise ratio comparison of encoding methods for hyperpolarized noble gas MRI.

Some non-Fourier encoding methods such as wavelet and direct encoding use spatially localized bases. The spatial localization feature of these methods enables optimized encoding for improved spatial and temporal resolution during dynamically adaptive MR imaging. These spatially localized bases, however, have inherently reduced image signal-to-noise ratio compared with Fourier or Hadamad encoding for proton imaging. Hyperpolarized noble gases, on the other hand, have quite different MR properties compared to proton, primarily the nonrenewability of the signal. It could be expected, therefore, that the characteristics of image SNR with respect to encoding method will also be very different from hyperpolarized noble gas MRI compared to proton MRI. In this article, hyperpolarized noble gas image SNRs of different encoding methods are compared theoretically using a matrix description of the encoding process. It is shown that image SNR for hyperpolarized noble gas imaging is maximized for any orthonormal encoding method. Methods are then proposed for designing RF pulses to achieve normalized encoding profiles using Fourier, Hadamard, wavelet, and direct encoding methods for hyperpolarized noble gases. Theoretical results are confirmed with hyperpolarized noble gas MRI experiments.

Theoretical comparison of Fourier and wavelet encoding in magnetic resonance imaging.

A theoretical comparison of the wavelet and Fourier encoding methods is made with respect to resolution, sensitivity to artifact, and signal-to-noise ratios (SNR's). A general mathematical description is developed in which magnetic resonance (MR) image encoding is represented by a projection of the function representing the "MR signal density" onto an approximation subspace of the finite energy functions. Characteristics of the subspace are used to define a "generalized" point-spread function for space-variant systems. Using the formal model of MR image encoding it is shown that wavelet encoding approaches the resolution limit defined by Fourier encoding. Artifact is treated according to whether or not the source of the variation in the measured data is stationary. Nonstationary imperfections perturb the projection operation and result in encoding method-dependent effects which can be modeled by a distortion matrix suitable for treating shift variant systems. A framework is developed in which to derive expressions for SNR's applicable to a general class of MR encoding methods.

Spin-Echo planar spectroscopic imaging for fast lipid characterization in bone marrow.

Lipid characterization of bone marrow in vivo with proton magnetic resonance spectroscopy was performed using Spin-Echo Planar Spectroscopic Imaging sequences. The methods are shown capable of rapidly generating two-dimensional chemical shift imaging data sets suitable for measuring lipid indices that reflect unsaturation levels among triglycerides, as demonstrated in oil phantoms and bone marrow from a healthy volunteer. The volume coverage, spatial resolution, acquisition speed, and spectral characteristics of Spin-Echo Planar Spectroscopic Imaging should make it attractive for clinical studies of diseases affecting normal lipid chemical composition.

An efficient region of interest acquisition method for dynamic magnetic resonance imaging.

Motivated by work in the area of dynamic magnetic resonance imaging (MRI), we develop a new approach to the problem of reduced-order MRI acquisition. Efforts in this field have concentrated on the use of Fourier and singular value decomposition (SVD) methods to obtain low-order representations of an entire image plane. We augment this work to the case of imaging an arbitrarily-shaped region of interest (ROI) embedded within the full image. After developing a natural error metric for this problem, we show that determining the minimal order required to meet a prescribed error level is in general intractable, but can be solved under certain assumptions. We then develop an optimization approach to the related problem of minimizing the error for a given order. Finally, we demonstrate the utility of this approach and its advantages over existing Fourier and SVD methods on a number of MRI images.

Investigation of long-term reproducibility of intrinsic connectivity network mapping: a resting-state fMRI study.

BACKGROUND AND PURPOSE: Connectivity mapping based on resting-state fMRI is rapidly developing, and this methodology has great potential for clinical applications. However, before resting-state fMRI can be applied for diagnosis, prognosis, and monitoring treatment for an individual patient with neurologic or psychiatric diseases, it is essential to assess its long-term reproducibility and between-subject variations among healthy individuals. The purpose of the study was to quantify the long-term test-retest reproducibility of ICN measures derived from resting-state fMRI and to assess the between-subject variation of ICN measures across the whole brain. MATERIALS AND METHODS: Longitudinal resting-state fMRI data of 6 healthy volunteers were acquired from 9 scan sessions during >1 year. The within-subject reproducibility and between-subject variation of ICN measures, across the whole brain and major nodes of the DMN, were quantified with the ICC and COV. RESULTS: Our data show that the long-term test-retest reproducibility of ICN measures is outstanding, with >70% of the connectivity networks showing an ICC > 0.60. The COV across 6 healthy volunteers in this sample was >0.2, suggesting significant between-subject variation. CONCLUSIONS: Our data indicate that resting-state ICN measures (eg, the correlation coefficients between fMRI signal-intensity profiles from 2 different brain regions) are potentially suitable as biomarkers for monitoring disease progression and treatment effects in clinical trials and individual patients. Because between-subject variation is significant, it may be difficult to use quantitative ICN measures in their current state as a diagnostic tool.

Abnormal thalamic volume in treatment-naïve boys with Tourette syndrome.

OBJECTIVE: Thalamic abnormality has been implicated in the pathophysiology of Tourette's syndrome (TS). We examined the presence of aberrant thalamic volume from the treatment-naïve boys with TS using magnetic resonance imaging (MRI). METHOD: Volumetric MRI was performed on 18 treatment-naïve boys with TS, aged 7-14 years, and 16 healthy comparison subjects. The anatomical boundaries were then manually parcellated to measure the thalamic volume. RESULTS: Tourette's syndrome subjects had a significantly larger left thalamus in comparison with those of healthy subjects. On the contrary, no group difference was observed from the right thalamic volume. TS subjects also showed a significant reduction in rightward asymmetry in thalamic volume compared with the healthy subjects. CONCLUSION: Our findings provide new evidence of abnormal thalamic volume in pediatric TS.

A dynamically adaptive imaging algorithm for wavelet-encoded MRI.

A new adaptive algorithm based on wavelet-encoded MRI is presented for application in dynamic imaging. This algorithm is adaptive because the strategy for updating image data in the dynamic series of images is determined by the processing of the most recently acquired data. The spatially selective multi-resolution properties of the wavelet transform are exploited to selectively update only those regions of the field of view where change is actually occurring. A theoretical imaging model is presented to motivate use of the adaptive algorithm, and simulation results using both artificial and experimental wavelet-encoded data are presented.

Selection of high-definition 2D virtual profiles with multiple RF pulse excitations along interleaved echo-planar k-space trajectories.

A method for spatially selective excitation of 2D RF profiles is reported. The method makes use of multiple shots to traverse interleaved echo-planar trajectories in 2D k space during each RF pulse excitation. Results from each of the interleaved excitations are summed, with the net effect being the excitation of a virtual profile. The method allows for the excitation of high-definition 2D profiles with standard gradient hardware. Signal to noise is enhanced by a factor equal to the square root of the number of interleaved excitations, compared with a single-shot excitation. Potential applications for volume-localized spectroscopy, functional MRI, and high-resolution reduced-field-of-view imaging are discussed.

Dynamically adaptive MRI with encoding by singular value decomposition.

A new MRI spatial encoding method based upon the singular value decomposition (SVD) and using spatially selective RF excitation is described. This encoding technique is particularly applicable to dynamic adaptive MRI, because it provides a near minimal set of spatial encoding profiles computed using an image estimate that is determined from a previously obtained image. Experimental results are presented for two cases, which exemplify its potential use in different dynamic imaging tasks. SVD-encoded MRI has demonstrated to be a highly efficient encoding scheme.

Implementation of wavelet-encoded MR imaging.

Reconstructions of images from wavelet-encoded data are shown. The method of MR wavelet encoding in one dimension was proposed previously by Weaver and Healy. The technique relies on selective excitation with wavelet-shaped profiles generated by special radio-frequency waveforms. The result of the imaging sequence is a set of inner products of the image with orthogonal functions of the wavelet basis. Inversion of the wavelet data is accomplished with an efficient algorithm with processing times comparable with those of a fast Fourier transform. The experiments show that wavelet encoding by selective excitation of wavelet-shaped profiles is feasible. Wavelet-encoded images are compared with phase-encoded images that have a similar signal-to-noise ratio, and there is no discernible degradation in image quality due to the wavelet encoding. Potential benefits of wavelet encoding are briefly discussed.

Multiresolution data acquisition and detection in functional MRI.

In an investigation of a multiresolution and multistaged approach in functional MRI, the relationship between spatial resolution and detection of functional activation is examined. The difference between functional detection and mapping is defined, and a multiresolution approach to functional detection is analyzed by constructing simple theoretical and experimental models simulating variations of in-plane resolution. Experimentally measured blood oxygenation level-dependent (BOLD) signal changes as well as BOLD contrast-to-noise ratio (CNR) with respect to different spatial resolutions are compared with results from theoretical predictions and simulation. From both an experimental and a theoretical perspective, it is shown that BOLD CNR and, thus, the concomitant detection of the functional activation are maximized when the resolution matches the size of activation.

Applicability and efficiency of near-optimal spatial encoding for dynamically adaptive MRI.

Adaptive near-optimal MRI spatial encoding entails, for the acquisition of each image update in a dynamic series, the computation of encodes in the form of a linear algebra-derived orthogonal basis set determined from an image estimate. The origins of adaptive encoding relevant to MRI are reviewed. Sources of error of this approach are identified from the linear algebraic perspective where MRI data acquisition is viewed as the projection of information from the field-of-view onto the encoding basis set. The definitions of ideal and non-ideal encoding follow, with nonideal encoding characterized by the principal angles between two vector spaces. An analysis of the distribution of principal angles is introduced and applied in several example cases to quantitatively describe the suitability of a basis set derived from a specific image estimate for the spatial encoding of a given field-of-view. The robustness of adaptive near-optimal spatial encoding for dynamic MRI is favorably shown by results computed using singular value decomposition encoding that simulates specific instances of worst case data acquisition when all objects have changed or new objects have appeared in the field-of-view. The mathematical analysis and simulations presented clarify the applicability and efficiency of adaptively determined near-optimal spatial encoding throughout a range of circumstances as may typically occur during use of dynamic MRI.

Functional magnetic resonance imaging using non-Fourier, spatially selective radiofrequency encoding.

A new method for functional magnetic resonance imaging (fMRI) employing non-Fourier encoding using spatially selective radiofrequency (RF) excitation is presented. The method uses manipulation of spatially selective RF pulses to encode spins in the slice-select direction. The method has several advantages over standard multislice approaches. It provides a simple means for monitoring irregularly distributed sections throughout a volume without the need to encode the whole volume. It offers the potential for increased signal-to-noise ratio if an appropriate basis is used for encoding. With a unique design of excitation pulses, it also appears possible to significantly reduce in-flow effects. An interleaved echo-planar imaging (EPI) sequence was adapted for non-Fourier encoding in the slice-select direction and was implemented on a conventional 1.5-Telsa system. The method was then used for functional mapping of the visual and motor areas where significant reduction of in-flow effect was demonstrated. This approach can be adapted to other imaging sequences that are used for fMRI, such as single-shot EPI.

A simple noninvasive stereotactic device for routine MR head examinations.

A simple device to define reproducible reference points in MR examinations was designed. The device can be worn like a pair of eyeglasses and has reference structure visible in MR images. Combined with a specially designed algorithm, the imaging plane can be aligned to a predefined coordinate system in one step, using a single slice through the device. The device can be produced at low cost and is suitable for routine examinations as well as surgical planning.

Automation of the seizure investigation unit at the University of British Columbia Health Sciences Centre Hospital.

At the University of British Columbia Health Sciences Centre Hospital we have constructed an automated Seizure Investigation Unit for long term monitoring of epileptic patients. A central component of this system is a new device, developed at the University of British Columbia, which prevents video and EEG records of seizures recorded on video tape from being over-recorded. Computer technology is relied upon to a considerable degree in our unit. Computers are seen, in this phase of development of the SIU, as a means of making patient monitoring less dependent on supervision. They can help to redirect human energy towards complex analysis rather than time consuming and simple monitoring tasks.

Practical digital filters for reducing EMG artefact in EEG seizure recordings.

In long-term scalp EEG monitoring of epileptic patients it is virtually impossible, in the present state of the technology, to avoid movement-related artefacts. These often obscure EEG information about the location of the seizure focus. One important example of such artefact is EMG activity. Its removal or suppression is sometimes enough to make otherwise useless EEG traces readable. Different methods of filtering have been applied towards that end. We routinely use a 15 Hz setting on our polygraph in obtaining EEG seizure printouts. We have recently examined digital filters which attenuate EMG beyond what is possible with the 15 Hz filter. A concern has been that the filters are practical in that they run in real time on a simple microprocessor and cause a minimum of confusion between smoothed artefact and actual brain activity.

Comparing the FAISE method with conventional dual-echo sequences.

The FAISE (fast-acquisition interleaved spin-echo) technique consists of a hybrid rapid-acquisition relaxation-enhanced (RARE) sequence combined with a specific phase-encode reordering method. Implemented on a 1.5-T unit, this multisection, high-resolution technique permits convenient contrast manipulation similar to that of spin-echo imaging, with selection of a pseudo-echo-time parameter and a TR interval. With a TR of 2 seconds, eight 256 x 256 images are obtained in 34 seconds with either T2 or proton-density weighting. A direct comparison between FAISE and spin echo for obtaining T2-weighted head images in healthy subjects indicates that FAISE and spin-echo images are qualitatively and quantitatively similar. Image artifacts are more pronounced on "proton-density" FAISE images than on the T2-weighted FAISE images. T1 contrast can be obtained with inversion recovery and short TR FAISE images. Preliminary temperature measurements in saline phantoms do not indicate excessive temperature increases with extended FAISE acquisitions. However, extensive studies of radio-frequency power deposition effects should be performed if the FAISE technique is to be fully exploited.

Dynamic imaging with multiple resolutions along phase-encode and slice-select dimensions.

An implementation is reported of an imaging method to obtain MUltiple Resolutions along Phase-encode and Slice-select dimensions (MURPS), which enables dynamic imaging of focal changes using a graded, multiresolution approach. MURPS allows one to trade spatial resolution in part of the volume for improved temporal resolution in dynamic imaging applications. A unique method of Hadamard slice encoding is used, enabling the varying of the phase encode and slice resolution while maintaining a constant effective TR throughout the entire 3-D volume. MURPS was implemented using a gradient-recalled echo sequence, and its utility was demonstrated for MR temperature monitoring. In this preliminary work, it has been shown that changes throughout a large volume can be effectively monitored in times that would normally only permit dynamic imaging in one or a very few slices.

Non-Fourier encoding with multiple spin echoes.

The advantages and limitations of multiple spin-echo sequences for non-Fourier encoding are investigated. Complications caused by improper encoding of alternate magnetization pathways due to imperfect refocusing pulses are analyzed. It is shown that mirror image ghosts result if the encoding RF pulse matrix is real-valued. These ghosts can be avoided as long as the rows of the RF pulse matrix are conjugate symmetric, which implies that spatial profiles are real valued. Non-Fourier encoding using bases derived from wavelet, Hadamard, and other real-valued orthogonal functions does not result in a mirror ghost artifact. A RARE sequence for non-Fourier encoding has been implemented on a clinical imaging system and successfully applied for brain imaging.