Maxim Zaitsev

To further improve the quality and robustness of the point-spread function (PSF) mapping method for fully automatic and accurate correction of geometric distortions in EPI at ultra high field such as 7 Tesla with high fidelity. Conventional gradient-echo EPI and corresponding PSF reference data in phantoms and in in vivo measurements of the human brain were acquired at 7 Tesla. To accurately determine and correct geometric distortions, i.e., in peripheral areas, the method calculates the local shift in the distorted phase-encoding dimension instead of the non-distorted spin-wrap encoding dimension. The results of the proposed method are compared with those originally demonstrated (Zaitsev et al. Magn Reson Med (2004) 52:1156-1166). The results show that the proposed method allows measurement and correction of the geometric distortions in EPI with high accuracy, and reduction of residual blurring. In addition, this method prevents errors in the shift map induced by PSF-ghost artifacts. More precise mapping and correction of pixel shifts and blurring is accomplished with the proposed improvements. Errors in the shift map that are caused by PSF-ghost artifacts in the multi-shot PSF acquisition, e.g., from small motion during the reference scan, do not affect the improved shift estimation method.

Multiple nonlinear gradient fields offer many potential benefits for spatial encoding including reduced acquisition time, fewer artefacts and region-specific imaging, although designing a suitable trajectory for such a setup is difficult. This work aims to optimize encoding trajectories for multiple nonlinear gradient fields based on the image signal-to-noise ratio. Image signal-to-noise ratio is directly linked to the covariance of the reconstructed pixels, which can be calculated recursively for each projection of the trajectory under a Bayesian formulation. An evolutionary algorithm is used to find the higher-dimensional projections that minimize the pixel covariance, incorporating receive coil profiles, intravoxel dephasing, and reconstruction regularization. The resulting trajectories are tested through simulations and experiments. The optimized trajectories produce images with higher resolution and fewer artefacts compared with traditional approaches, particularly for high undersampling. However, higher-dimensional projection experiments strongly depend on accurate hardware and calibration. Computer-based optimization provides an efficient means to explore the large trajectory space created by the use of multiple nonlinear encoding fields. The optimization framework, as presented here, is necessary to fully exploit the advantages of nonlinear fields. Magn Reson Med, 2015. © 2015 Wiley Periodicals, Inc.

High field MRI systems, such as 7 Tesla (T) scanners, can deliver higher signal to noise ratio (SNR) than lower field scanners and thus allow for the acquisition of data with higher spatial resolution, which is often demanded by users in the fields of clinical and neuroscientific imaging. However, high resolution scans may require long acquisition times, which in turn increase the discomfort for the subject and the risk of subject motion. Even with a cooperative and trained subject, involuntary motion due to heartbeat, swallowing, respiration and changes in muscle tone can cause image artifacts that reduce the effective resolution. In addition, scanning with higher resolution leads to increased sensitivity to even very small movements. Prospective motion correction (PMC) at 3T and 7T has proven to increase image quality in case of subject motion. Although the application of prospective motion correction is becoming more popular, previous articles focused on proof of concept studies ...

In this paper we present a monoplanar gradient system capable of imaging a volume comparable with that covered by linear gradient systems. Such a system has been designed and implemented. Building such a system was made possible by relaxing the constraint of global linearity and replacing it with a requirement for local orthogonality. A framework was derived for optimization of local orthogonality within the physical boundaries and geometric constraints. Spatial encoding of magnetic fields was optimized for their local orthogonality over a large field of view. A coil design consisting of straight wire segments was optimized, implemented, and integrated into a 3T human scanner to show the feasibility of this approach. Initial MR images are shown and further applications of the derived optimization method and the nonlinear planar gradient system are discussed. Encoding fields generated by the prototype encoding system were shown to be locally orthogonal and able to encode a cylindrica...

To accelerate dynamic imaging of the vocal tract during articulation. Five subjects were imaged by magnetic resonance imaging (MRI) while repeating the word "Partikel" at 90 beats per minute, using both a radio-frequency-spoiled radial gradient echo sequence with golden angle projection rotation and a previously applied real-time Cartesian sequence. The acquired data were reconstructed using a CG-SENSE method and Cartesian reconstruction. The images from both methods were compared by measuring distances between anatomical landmarks that are important for resonance behavior of the vocal tract. Only commonly available hardware and software were used. With the presented radial method a spatial resolution of 1.8 mm at 25 frames per second could be achieved. Overall, the morphometric measurements of the vocal tract showed less deviation for the radial sequence both across repetitions and for all subjects. Fast modifications, such as complete lip closing, could be analyzed with ...

In most half-Fourier imaging methods, explicit phase replacement is used. In combination with parallel imaging, or compressed sensing, half-Fourier reconstruction is usually performed in a separate step. The purpose of this paper is to report that integration of half-Fourier reconstruction into iterative reconstruction minimizes reconstruction errors. The L1-norm phase constraint for half-Fourier imaging proposed in this work is compared with the L2-norm variant of the same algorithm, with several typical half-Fourier reconstruction methods. Half-Fourier imaging with the proposed phase constraint can be seamlessly combined with parallel imaging and compressed sensing to achieve high acceleration factors. In simulations and in in-vivo experiments half-Fourier imaging with the proposed L1-norm phase constraint enables superior performance both reconstruction of image details and with regard to robustness against phase estimation errors. The performance and feasibility of half-Fourier ...

We aimed to demonstrate that follow-up scans in longitudinal examinations can be significantly accelerated by using images from previous scans as priors for constrained reconstruction. In this work, we propose a method for incorporating a prior image to improve the reconstruction of a new acquisition with considerable k-space undersampling, which contains a two-level registration scheme with non-parametric transformation, an adaptive synthesis procedure, and a constrained reconstruction with weighted total variation constraint. The performance of the method is evaluated using simulations, as well as results from volunteer and patient examinations. In vivo experiments with both volunteers and patients show that incorporating a prior image into the constrained reconstruction produces many fewer reconstruction errors compared to the conventional reconstruction using only the highly undersampled k-space data. The redundant information in the prior image can be efficiently adopted to imp...

A sequence for T(1) relaxation-time mapping which enables high-resolution, multislice imaging in short acquisition times is presented. The sequence is based on the Look-Locker method and employs a magnetization-preparation module prior to data acquisition with a banded k-space data collection scheme. The method was implemented on a standard clinical scanner and the accuracy of the T(1) results was evaluated against spectroscopic measurements. The accuracy of the T(1) maps validated by phantom imaging measurements is <3% for slow-relaxing compartments (T(1) congruent with 2000 ms) and is around 1% for faster-relaxing species (T(1) < or = 1200 ms). Additionally, the inherent multislice, multipoint capability of the method is demonstrated. Multislice, multipoint in vivo results of the human brain obtained using this method are presented. An acquisition time of approximately 8 min was achieved for a T(1) map, which, in principle, can provide whole-brain coverage with 25 slices, a matrix size of 256 x 256, and 12 time points. The speed of the sequence is derived through optimized interleaving of slices and time points, together with the acquisition of multiple echoes, which are used to fill a 3-segment k-space.

Most functional magnetic resonance imaging (fMRI) experiments use gradient-echo echo planar imaging (GE EPI) to detect the blood oxygenation level-dependent (BOLD) effect. This technique may fail in the presence of anatomy-related susceptibility-induced field gradients in the human head. In this work, we present a novel 3D compensation method in combination with a template-based correction that can be optimized over particular regions of interest to recover susceptibility-induced signal loss without acquisition time penalty. Based on an evaluation of B(0) field maps of eight subjects, slice-dependent gradient compensation moments are derived for maximal BOLD sensitivity in two compromised regions: the orbitofrontal cortex and the amygdala areas. A modified EPI sequence uses these additional gradient moments in all three imaging directions. The method is compared to non-compensated, template-based and subject-specific correction gradients and also in a breath-holding experiment. The slice-dependent gradient compensation method significantly improves signal intensity/BOLD sensitivity by about 35/43% in the orbitofrontal cortex and by 17/30% in the amygdala areas compared to a conventional acquisition. Template-based correction and subject-specific correction perform equally well. The BOLD sensitivity in the breath hold experiment is effectively increased in compensated regions. The new method addresses the problem of susceptibility-induced signal loss, without compromising temporal resolution. It can be used for event-related functional experiments without requiring additional subject-specific calibration or calculation time.

Subject motion is still the major source of data quality degradation in functional magnetic resonance imaging (fMRI) studies. Established methods correct motion between successive repetitions based on the acquired imaging volumes either retrospectively or prospectively. A fast, highly accurate, and prospective real-time correction method for fMRI using external optical motion tracking has been implemented. The head position is determined by means of an optical stereoscopic tracking system. The method corrects motion during the acquisition of an fMRI time series on a slice-by-slice basis by continuously updating the imaging volume position to follow the motion of the head. This method allows the measurement of fMRI data in the presence of significant motion during the acquisition of a single volume. Even without intentional motion, fMRI signal stability is maintained and higher sensitivity to detect activation is achieved without reducing specificity. With significant motion, only the proposed approach allowed detection of brain activation. The results show that the new method is superior to image-based correction methods, which fail in the case of fast or excessive motion.

Single-shot echo planar imaging (EPI) acquisitions at 7T are challenging due to increased distortions, signal dropouts, RF-power requirements, and reduced T2*. This study developed and tested pulse sequence and protocol modifications required to allow high resolution EPI for whole brain functional neuroimaging. Using geometric distortion correction methods, modified fat saturation, and parallel imaging, we acquired high resolution single-shot gradient-echo EPI data at 7T with different spatial resolution. The BOLD sensitivity was evaluated and quantified in a breath hold experiment. Single-shot EPI data with isotropic resolution from 3 to 1.1 mm were acquired in human subjects. The RF-power deposition has been reduced to allow up to 22 slices per second. In addition, acoustic noise and helium boil-off have been reduced. A reduction of the fat saturation flip angle resulted in up to 20% signal gain without compromising the fat suppression quality. For the coil used, the BOLD sensitivity is highest for 2 or 1.4 mm isotropic resolution. High resolution single-shot EPI in the whole brain can be performed at 7T with high efficiency, low signal dropout, and without major geometric distortions.

A novel fat-suppressed balanced steady-state free precession (b-SSFP) imaging method based on the transition into driven equilibrium (TIDE) sequence with variable flip angles is presented. The new method, called fat-saturated (FS)-TIDE, exploits the special behavior of TIDE signals from off-resonance spins during the flip angle ramp. As shown by simulations and experimental data, the TIDE signal evolution for off-resonant isochromats during the transition from turbo spin-echo (TSE)-like behavior to the true fast imaging with steady precession (TrueFISP) mode undergoes a zero crossing. The resulting signal notch for off-resonant spins is then used for fat suppression. The efficiency of FS-TIDE is demonstrated in phantoms and healthy volunteers on a 1.5T system. The resulting images are compared with standard TrueFISP data with and without fat suppression. It is demonstrated that FS-TIDE provides a fast and stable means for homogenous fat suppression in abdominal imaging while maintaining balanced SSFP-like image contrast and signal-to-noise ratio (SNR). The scan time of FS-TIDE is not increased compared to normal TrueFISP imaging without fat suppression and identical k-space trajectories. Because of the intrinsic fat suppression, no additional preparation is needed. Possible repetition times (TRs) are not firmly limited to special values and are nearly arbitrary.

Most existing methods for accelerated parallel imaging in MRI require additional data, which are used to derive information about the sensitivity profile of each radiofrequency (RF) channel. In this work, a method is presented to avoid the acquisition of separate coil calibration data for accelerated Cartesian trajectories. Quadratic phase is imparted to the image to spread the signals in k-space (aka phase scrambling). By rewriting the Fourier transform as a convolution operation, a window can be introduced to the convolved chirp function, allowing a low-resolution image to be reconstructed from phase-scrambled data without prominent aliasing. This image (for each RF channel) can be used to derive coil sensitivities to drive existing parallel imaging techniques. As a proof of concept, the quadratic phase was applied by introducing an offset to the x(2) - y(2) shim and the data were reconstructed using adapted versions of the image space-based sensitivity encoding and GeneRalized Autocalibrating Partially Parallel Acquisitions algorithms. The method is demonstrated in a phantom (1 × 2, 1 × 3, and 2 × 2 acceleration) and in vivo (2 × 2 acceleration) using a 3D gradient echo acquisition. Phase scrambling can be used to perform parallel imaging acceleration without acquisition of separate coil calibration data, demonstrated here for a 3D-Cartesian trajectory. Further research is required to prove the applicability to other 2D and 3D sampling schemes. Magn Reson Med 73:1407-1419, 2015. © 2014 Wiley Periodicals, Inc.

Time-dependent phenomena are of great interest, and researchers have sought to shed light on these processes with MRI, particularly in vivo. In this work, a new hybrid technique based on EPI and using the concept of keyhole imaging is presented. By sharing peripheral k-space data between images and acquiring the keyhole more frequently, it is shown that the spatial resolution of the reconstructed images can be maintained. The method affords a higher temporal resolution and is more robust against susceptibility and chemical-shift artifacts than single-shot EPI. The method, termed shared k-space echo planar imaging with keyhole (shared EPIK), has been implemented on a standard clinical scanner. Technical details, simulation results, phantom images, in vivo images, and fMRI results are presented. These results indicate that the new method is robust and may be used for dynamic MRI applications. Magn Reson Med 45:109-117, 2001.

Nonlinear spatial encoding magnetic fields allow excitation and geometrically matched local encoding of curved slices. However, the nonlinearity of the fields results in a varying slice thickness. Within this study, the technique is combined with multidimensional RF excitation for local adaptation of the slice shape. A framework originally developed for nonlinear receive encoding is applied to multidimensional excitation with nonlinear spatial encoding magnetic fields for determination of dedicated target patterns and combined with a model for assessment of minimum transmit-resolution requirements for the design of efficient transmit k-space trajectories. Cross-sections of curved slices acquired in a phantom with both locally adapted slice thickness and curvature are evaluated. In addition, resulting voxel shapes are analyzed to investigate the range of applicability of the technique. Finally, slice-thickness adaptation is applied to in vivo curved slice imaging. Local adaptation of the slice thickness is feasible both in phantom and in vivo. The technique further allows local adaptation of the slice curvature. However, its range of applicability is limited by prolonged pulse duration and voxel shape distortion. Multidimensional excitation allows imaging of curved slices with constant thickness. It also has the potential for further modification of the slice shape for increased ability to adapt to the anatomy.

Echo-planar imaging (EPI) is an ultrafast magnetic resonance (MR) imaging technique prone to geometric distortions. Various correction techniques have been developed to remedy these distortions. Here improvements of the point spread function (PSF) mapping approach are presented, which enable reliable and fully automated distortion correction of echo-planar images at high field strengths. The novel method is fully compatible with EPI acquisitions using parallel imaging. The applicability of parallel imaging to further accelerate PSF acquisition is shown. The possibility of collecting PSF data sets with total acceleration factors higher than the number of coil elements is demonstrated. Additionally, a new approach to visualize and interpret distortions in the context of various imaging and reconstruction methods based on the PSF is proposed. The reliable performance of the PSF mapping technique is demonstrated on phantom and volunteer scans at field strengths of up to 4 T.

Flow-sensitive 3-dimensional magnetic resonance imaging at 3 T and advanced 3-dimensional visualization were used to visualize local and global vascular hemodynamics in the thoracic aorta. In patients with pathological geometric alterations of the thoracic aorta, this technique revealed considerable changes in local blood flow characteristics, compared with normal volunteers. Specifically, relatively small geometric changes, such as a partially thrombosed aortic arch or a mild aneurysm of the ascending aorta, resulted in major disturbances of local blood flow patterns within and even further downstream to the pathology.