K-space Lines Research Articles

The PRESTO technique for fMRI

In the early days of BOLD fMRI, the acquisition of T(2)(*) weighted data was greatly facilitated by rapid scan techniques such as EPI. The latter, however, was only available on a few MRI systems that were equipped with specialized hardware that allowed rapid switching of the imaging gradients. For this reason, soon after the invention of fMRI, the scan technique PRESTO was developed to make rapid T(2)(*) weighted scanning available on standard clinical scanners. This method combined echo shifting, which allows for echo times longer than the sequence repetition time, with acquisition of multiple k-space lines per excitation. These two concepts were combined in order to achieve a method fast enough for fMRI, while maintaining a sufficiently long echo time for optimal contrast. PRESTO has been primarily used for 3D scanning, which minimized the contribution of large vessels due to inflow effects. Although PRESTO is still being used today, its appeal has lessened somewhat due to increased gradient performance of modern MRI scanners. Compared to 2D EPI, PRESTO may have somewhat reduced temporal stability, which is a disadvantage for fMRI that may not outweigh the advantage of reduced inflow effects provided by 3D scanning. In this overview, the history of the development of the PRESTO is presented, followed by a qualitative comparison with EPI.

A Technique to Reduce Motion Artifact for Externally Triggered Cine-MRI(EC-MRI) Based on Detecting the Onset of the Articulated Word with Spectral Analysis

One issue in externally triggered cine-magnetic resonance imaging (EC-MRI) for the dynamic observation of speech organs is motion artifact in the phase-encoding direction caused by unstable repetitions of speech during data acquisition. We propose a technique to reduce such artifact by rearranging the k-space data used to reconstruct MR images based on the analysis of recorded speech sounds. We recorded the subject's speech sounds during EC-MRI and used post hoc acoustical processing to reduce scanning noise and detect the onset of each utterance based on analysis of the recorded sounds. We selected each line of k-space from several data acquisition sessions and rearranged them to reconstruct a new series of dynamic MR images according to the analyzed time of utterance onset. Comparative evaluation showed significant reduction in motion artifact signal in the dynamic MR images reconstructed by the proposed method. The quality of the reconstructed images was sufficient to observe the dynamic aspects of speech production mechanisms.

Contrast enhancement in TOF cerebral angiography at 7 T using saturation and MT pulses under SAR constraints: Impact of VERSE and sparse pulses

Cerebral three-dimensional time of flight (TOF) angiography significantly benefits from ultrahigh fields, mainly due to higher signal-to-noise ratio and to longer T(1) relaxation time of static brain tissues; however, specific absorption rate (SAR) significantly increases with B(0). Thus, additional radiofrequency pulses commonly used at lower field strengths to improve TOF contrast such as saturation of venous signal and improved background suppression by magnetization transfer typically cannot be used at higher fields. In this work, we aimed at reducing SAR for each radiofrequency pulse category in a TOF sequence. We use the variable-rate selective excitation principle for the slab selective TOF excitation as well as the venous saturation radiofrequency pulses. In addition, magnetization transfer pulses are implemented by sparsely applying the pulses only during acquisition of the central k-space lines to limit their SAR contribution. Image quality, angiographic contrast, and SAR reduction were investigated as a function of variable-rate selective excitation parameters and of the total number of magnetization transfer pulses applied. Based on these results, a TOF protocol was generated that increases the angiographic contrast by more than 50% and reduces subcutaneous fat signal while keeping the resulting SAR within regulatory limits.

Use of 2D Sensitivity Encoding for Slow-Infusion Contrast-Enhanced Isotropic 3-T Whole-Heart Coronary MR Angiography

The purpose of this study was to improve the blood-pool signal-to-noise ratio (SNR) and blood-myocardium contrast-to-noise ratio (CNR) of slow-infusion 3-T whole-heart coronary MR angiography (MRA). In 2D sensitivity encoding (SENSE), the number of acquired k-space lines is reduced, allowing less radiofrequency excitation per cardiac cycle and a longer TR. The former can be exploited for signal enhancement with a higher radiofrequency excitation angle, and the latter leads to noise reduction due to lower data-sampling bandwidth. Both effects contribute to SNR gain in coronary MRA when spatial and temporal resolution and acquisition time remain identical. Numeric simulation was performed to select the optimal 2D SENSE pulse sequence parameters and predict the SNR gain. Eleven patients underwent conventional unenhanced and the proposed 2D SENSE contrast-enhanced coronary MRA acquisition. Blood-pool SNR, blood-myocardium CNR, visible vessel length, vessel sharpness, and number of side branches were evaluated. Consistent with the numeric simulation, using 2D SENSE in contrast-enhanced coronary MRA resulted in significant improvement in aortic blood-pool SNR (unenhanced vs contrast-enhanced, 37.5 ± 14.7 vs 121.3 ± 44.0; p < 0.05) and CNR (14.4 ± 6.9 vs 101.5 ± 40.8; p < 0.05) in the patient sample. A longer length of left anterior descending coronary artery was visualized, but vessel sharpness, coronary artery coverage, and image quality score were not improved with the proposed approach. In combination with contrast administration, 2D SENSE was found effective in improving SNR and CNR in 3-T whole-heart coronary MRA. Further investigation of cardiac motion compensation is necessary to exploit the SNR and CNR advantages and to achieve submillimeter spatial resolution.

Compressed-sensing motion compensation (CosMo): a joint prospective-retrospective respiratory navigator for coronary MRI.

Prospective right hemidiaphragm navigator (NAV) is commonly used in free-breathing coronary MRI. The NAV results in an increase in acquisition time to allow for resampling of the motion-corrupted k-space data. In this study, we are presenting a joint prospective-retrospective NAV motion compensation algorithm called compressed-sensing motion compensation (CosMo). The inner k-space region is acquired using a prospective NAV; for the outer k-space, a NAV is only used to reject the motion-corrupted data without reacquiring them. Subsequently, those unfilled k-space lines are retrospectively estimated using compressed sensing reconstruction. We imaged right coronary artery in nine healthy adult subjects. An undersampling probability map and sidelobe-to-peak ratio were calculated to study the pattern of undersampling, generated by NAV. Right coronary artery images were then retrospectively reconstructed using compressed-sensing motion compensation for gating windows between 3 and 10 mm and compared with the ones fully acquired within the gating windows. Qualitative imaging score and quantitative vessel sharpness were calculated for each reconstruction. The probability map and sidelobe-to-peak ratio show that the NAV generates a random undersampling k-space pattern. There were no statistically significant differences between the vessel sharpness and subjective score of the two reconstructions. Compressed-sensing motion compensation could be an alternative motion compensation technique for free-breathing coronary MRI that can be used to reduce scan time.

Increasing spoiling efficiency in RF‐spoiled gradient echo sequences by averaging of RF phase‐cycle‐adaptedk‐spaces

A method is presented that improves the effectiveness of RF-spoiling whenever an image consists of the averaging of several identical single acquisitions. Using this method, it becomes possible to reduce the necessary moment of the spoiler gradient by a factor equal to the number of averages. The idea is to perform the collection of k-space lines such that the averaging takes over a part of the spoiling process. The number of sequence cycles played out between two acquisitions of a certain k-space line is adjusted based on both the number of total averages of this line and the phase-cycling pattern of the RF-pulse train of a RF-spoiled gradient echo sequence. In this way, ghost artifacts visible in the single images add to zero in the averaged image while the depiction of the object experiences normal averaging with enhanced signal-to-noise ratio.

Quantitative T*2‐mapping based on multi‐slice multiple gradient echo flash imaging: Retrospective correction for subject motion effects

Numerous clinical and research applications for quantitative mapping of the effective transverse relaxation time T*(2) have been described. Subject motion can severely deteriorate the quality and accuracy of results. A correction method for T*(2) maps acquired with multi-slice multiple gradient echo FLASH imaging is presented, based on acquisition repetition with reduced spatial resolution (and consequently reduced acquisition time) and weighted averaging of both data sets, choosing weighting factors individually for each k-space line to reduce the influence of motion. In detail, the procedure is based on the fact that motion artifacts reduce the correlation between acquired and exponentially fitted data. A target data set is constructed in image space, choosing the data yielding best correlation from the two acquired data sets. The k-space representation of the target is subsequently approximated as linear combination of original raw data, yielding the required weighting factors. As this method only requires a single acquisition repetition with reduced spatial resolution, it can be employed on any clinical system offering a suitable sequence with export of modulus and phase images. Experimental results show that the method works well for sparse motion, but fails for strong motion affecting the same k-space lines in both acquisitions.

Motion correction using coil arrays (MOCCA) for free‐breathing cardiac cine MRI

In this study, we present a motion correction technique using coil arrays (MOCCA) and evaluate its application in free-breathing respiratory self-gated cine MRI. Motion correction technique using coil arrays takes advantages of the fact that motion-induced changes in k-space signal are modulated by individual coil sensitivity profiles. In the proposed implementation of motion correction technique using coil arrays self-gating for free-breathing cine MRI, the k-space center line is acquired at the beginning of each k-space segment for each cardiac cycle with 4 repetitions. For each k-space segment, the k-space center line acquired immediately before was used to select one of the 4 acquired repetitions to be included in the final self-gated cine image by calculating the cross correlation between the k-space center line with a reference line. The proposed method was tested on a cohort of healthy adult subjects for subjective image quality and objective blood-myocardium border sharpness. The method was also tested on a cohort of patients to compare the left and right ventricular volumes and ejection fraction measurements with that of standard breath-hold cine MRI. Our data indicate that the proposed motion correction technique using coil arrays method provides significantly improved image quality and sharpness compared with free-breathing cine without respiratory self-gating and provides similar volume measurements compared with breath-hold cine MRI.

Improved cardiac arterial spin labelling in the mouse heart by optimisation of acquisition and analysis

Methods Perfusion maps are generated by pixel-wise comparison of T1 maps from slice-selective and global inversion recoveries. T1 was quantified using ECG-gated LookLocker approach [3], with 4 lines of k-space acquired every heartbeat to speed up acquisition by approximately four times compared with previous approaches (eg. [2],[3]). A data logger (CED) was used to record ECG, RF and respiration events and was programmed to gate on each ECG peak detected as well as to interpolate peaks missed by physiological monitoring software. Using the information recorded, lines of k-space were binned as ‘acquired during respiration’ or not. Images with many lines or centre lines of k-space acquired during respiration were rejected. ASL data were acquired on a 9.4T Varian system with 35mm-long volume RF coil (Rapid) (TE=1.18ms, flip=5°, resolution=200μm, slice thickness=1.5mm) and image analysis was performed in Matlab.

Two-Axis Acceleration of Functional Connectivity Magnetic Resonance Imaging by Parallel Excitation of Phase-Tagged Slices and Half k-Space Acceleration

Whole brain functional connectivity magnetic resonance imaging requires acquisition of a time course of gradient-recalled (GR) volumetric images. A method is developed to accelerate this acquisition using GR echo-planar imaging and radio frequency (RF) slice phase tagging. For N-fold acceleration, a tailored RF pulse excites N slices using a uniform-field transmit coil. This pulse is the Fourier transform of the profile for the N slices with a predetermined RF phase tag on each slice. A multichannel RF receive coil is used for detection. For n slices, there are n/N groups of slices. Signal-averaged reference images are created for each slice within each slice group for each member of the coil array and used to separate overlapping images that are simultaneously received. The time-overhead for collection of reference images is small relative to the acquisition time of a complete volumetric time course. A least-squares singular value decomposition method allows image separation on a pixel-by-pixel basis. Twofold slice acceleration is demonstrated using an eight-channel RF receive coil, with application to resting-state functional magnetic resonance imaging in the human brain. Data from six subjects at 3 T are reported. The method has been extended to half k-space acquisition, which not only provides additional acceleration, but also facilitates slice separation because of increased signal intensity of the central lines of k-space coupled with reduced susceptibility effects.

Self-calibrating Coil Sensitivity Profiles for Parallel Imaging Based on Anisotropic Diffusion

Calibration of the spatial sensitivity functions of coil arrays is a crucial element in parallel magnetic resonance imaging (PMRI).In order to overcome the sensitivity miscalibration errors introduced by pre-scan calibration method and decrease the total examination time , a new self-calibrating method for estimating the coil sensitivity profiles was proposed i n th e paper . F rom a fully - sampled central region of a variable-density k -space acquisition, t he sensitivity calibration images were directly extracted , and then were adaptively smoothed by a new anisotropic diffusion scheme to extract encoding effects of coil sensitivities. When the estimated coil sensitivity profiles were applied to reconstruct the full Field-Of-View image from the uniformly under-sampled MR data when under-sampl e rate was greater than 2 , the experimental results showed that the quality of reconstruct ion image was evidently improv ed. Therefore, t h is self-calibrating method for calculating sensitivity profiles is suitable for much faster imaging in parallel magnetic imaging. Meanwhile, the paper suggest that the suitable number of central k-space lines for generating the internal reference images should be chosen between 8 to 20.

Time‐resolved absolute velocity quantification with projections

Quantitative information on time-resolved blood velocity along the femoral/popliteal artery can provide clinical information on peripheral arterial disease and complement MR angiography as not all stenoses are hemodynamically significant. The key disadvantages of the most widely used approach to time-resolve pulsatile blood flow by cardiac-gated velocity-encoded gradient-echo imaging are gating errors and long acquisition time. Here, we demonstrate a rapid nontriggered method that quantifies absolute velocity on the basis of phase difference between successive velocity-encoded projections after selectively removing the background static tissue signal via a reference image. The tissue signal from the reference image's center k-space line is isolated by masking out the vessels in the image domain. The performance of the technique, in terms of reproducibility and agreement with results obtained with conventional phase contrast-MRI was evaluated at 3 T field strength with a variable-flow rate phantom and in vivo of the triphasic velocity waveforms at several segments along the femoral and popliteal arteries. Additionally, time-resolved flow velocity was quantified in five healthy subjects and compared against gated phase contrast-MRI results. To illustrate clinical feasibility, the proposed method was shown to be able to identify hemodynamic abnormalities and impaired reactivity in a diseased femoral artery. For both phantom and in vivo studies, velocity measurements were within 1.5 cm/s, and the coefficient of variation was less than 5% in an in vivo reproducibility study. In five healthy subjects, the average differences in mean peak velocities and their temporal locations were within 1 cm/s and 10 ms compared to gated phase contrast-MRI. In conclusion, the proposed method provides temporally resolved arterial velocity with a temporal resolution of 20 ms with minimal post processing.

Intrinsic detection of motion in segmented sequences

While many motion correction techniques for MRI have been proposed, their use is often limited by increased patient preparation, decreased patient comfort, additional scan time, or the use of specialized sequences not available on many commercial scanners. For this reason, we propose a simple self-navigating technique designed to detect motion in segmented sequences. We demonstrate that comparing two segments containing adjacent sets of k-space lines results in an aliased error function. A global shift of the aliased error function indicates the presence of in-plane rigid-body translation, while other types of motion are evident in the dispersion or breadth of the error function. Since segmented sequences commonly acquire data in sets of adjacent k-space lines, this method provides these sequences with an inherent method of detecting object motion. Motion corrupted data can then be reacquired proactively or in some cases corrected or removed retrospectively.

CEST‐FISP: A novel technique for rapid chemical exchange saturation transfer MRI at 7 T

Chemical exchange saturation transfer (CEST) and magnetization transfer techniques provide unique and potentially quantitative contrast mechanisms in multiple MRI applications. However, the in vivo implementation of these techniques has been limited by the relatively slow MRI acquisition techniques, especially on high-field MRI scanners. A new, rapid CEST-fast imaging with steady-state free precession technique was developed to provide sensitive CEST contrast in ∼20 sec. In this study at 7 T with in vitro bovine glycogen samples and initial in vivo results in a rat liver, the CEST-fast imaging with steady-state free precession technique was shown to provide equivalent CEST sensitivity in comparison to a conventional CEST-spin echo acquisition with a 50-fold reduction in acquisition time. The sensitivity of the CEST-fast imaging with steady-state free precession technique was also shown to be dependent on k-space encoding with centric k-space encoding providing a 30-40% increase in CEST sensitivity relative to linear encoding for 256 or more k-space lines. Overall, the CEST-fast imaging with steady-state free precession acquisition technique provides a rapid and sensitive imaging platform with the potential to provide quantitative CEST and magnetization transfer imaging data.

Analysis of time reduction methods for magnetic resonance elastography of the brain

Magnetic resonance elastography (MRE) uses a phase-contrast MRI technique to image shear wave propagation in tissue followed by the mathematical inversion of the equations of motion governing tissue mechanics to noninvasively image tissue stiffness. This work investigates the impact of various MR sampling strategies designed to reduce acquisition times on the accuracy of MRE inversions. The results indicate that brain MRE data can be significantly truncated while maintaining a mean global stiffness error less than 10%. The results also indicate that brain MRE data can be collected in as few as eight lines of k-space. This degree of data truncation is possible due to the relatively low spatial frequency content and low amplitude of the shear waves observed during brain MRE exams and will facilitate the design of rapid brain MRE protocols for future clinical investigations.

Optimized parallel imaging for dynamic PC‐MRI with multidirectional velocity encoding

Phase contrast MRI with multidirectional velocity encoding requires multiple acquisitions of the same k-space lines to encode the underlying velocities, which can considerably lengthen the total scan time. To reduce scan time, parallel imaging is often applied. In dynamic phase contrast MRI using standard generalized autocalibrating partially parallel acquisitions (GRAPPA), several central k-spaces for autocalibration of the reconstruction (autocalibrating signal lines (ACS)) are typically acquired, separately for each velocity direction and each cardiac timeframe, for calculating the reconstruction weights. To further accelerate data acquisition, we developed two methods, which calculated weights with a substantially reduced number of ACSl lines. The effects on image quality and flow quantification were compared to fully sampled data, standard GRAPPA, and time-interleaved sampling scheme in combination with generalized autocalibrating partially parallel acquisitions (TGRAPPA). The results show that the two proposed methods can clearly improve scan efficiency while maintaining image quality and accuracy of measured flow or myocardial tissue velocities. Compared to TGRAPPA, the proposed methods were more accurate in evaluating flow velocity. In conclusion, the proposed reconstruction strategies are promising for dynamic multidirectionally encoded acquisitions and can easily be implemented using the standard GRAPPA reconstruction algorithm.

Model‐based reconstruction for cardiac cine MRI without ECG or breath holding

This paper describes an acquisition and reconstruction strategy for cardiac cine MRI that does not require the use of electrocardiogram or breath holding. The method has similarities with self-gated techniques as information about cardiac and respiratory motion is derived from the imaging sequence itself; here, by acquiring the center k-space line at the beginning of each segment of a balanced steady-state free precession sequence. However, the reconstruction step is fundamentally different: a generalized reconstruction by inversion of coupled systems is used instead of conventional gating. By correcting for nonrigid cardiac and respiratory motion, generalized reconstruction by inversion of coupled systems (GRICS) uses all acquired data, whereas gating rejects data acquired in certain motion states. The method relies on the processing and analysis of the k-space central line data: local information from a 32-channel cardiac coil is used in order to automatically extract eigenmodes of both cardiac and respiratory motion. In the GRICS framework, these eigenmodes are used as driving signals of a motion model. The motion model is defined piecewise, so that each cardiac phase is reconstructed independently. Results from six healthy volunteers, with various slice orientations, show improved image quality compared to combined respiratory and cardiac gating.

Optimization of multiply acquired magnetic flux density Bz using ICNE-Multiecho train in MREIT

The aim of magnetic resonance electrical impedance tomography (MREIT) is to visualize the electrical properties, conductivity or current density of an object by injection of current. Recently, the prolonged data acquisition time when using the injected current nonlinear encoding (ICNE) method has been advantageous for measurement of magnetic flux density data, Bz, for MREIT in the signal-to-noise ratio (SNR). However, the ICNE method results in undesirable side artifacts, such as blurring, chemical shift and phase artifacts, due to the long data acquisition under an inhomogeneous static field. In this paper, we apply the ICNE method to a gradient and spin echo (GRASE) multi-echo train pulse sequence in order to provide the multiple k-space lines during a single RF pulse period. We analyze the SNR of the measured multiple Bz data using the proposed ICNE-Multiecho MR pulse sequence. By determining a weighting factor for Bz data in each of the echoes, an optimized inversion formula for the magnetic flux density data is proposed for the ICNE-Multiecho MR sequence. Using the ICNE-Multiecho method, the quality of the measured magnetic flux density is considerably increased by the injection of a long current through the echo train length and by optimization of the voxel-by-voxel noise level of the Bz value. Agarose-gel phantom experiments have demonstrated fewer artifacts and a better SNR using the ICNE-Multiecho method. Experimenting with the brain of an anesthetized dog, we collected valuable echoes by taking into account the noise level of each of the echoes and determined Bz data by determining optimized weighting factors for the multiply acquired magnetic flux density data.

SE_MC sequence improves image quality of carotid arteries and atherosclerotic plaques

Purpose FSE images suffer from blurring because of non-uniform T2-weighting of k-space [Constable, Gore.MRM.1992.28(1):9-24], which is stronger on small structures with short T2 times, such as carotid arteries. To avoid this, we propose to generate T1-PD-T2-weighted images using Spin-Echo Multi-Contrast (SE_MC) sequence. FSE and SE_MC are both CPMG multi-echo sequences; however FSE acquires different lines of the same k-space at different TEs, whereas SE_MC samples multiple k-spaces (with constant TE) at different TEs.

Navigator accuracy requirements for prospective motion correction

Prospective motion correction in MRI is becoming increasingly popular to prevent the image artifacts that result from subject motion. Navigator information is used to update the position of the imaging volume before every spin excitation so that lines of acquired k-space data are consistent. Errors in the navigator information, however, result in residual errors in each k-space line. This paper presents an analysis linking noise in the tracking system to the power of the resulting image artifacts. An expression is formulated for the required navigator accuracy based on the properties of the imaged object and the desired resolution. Analytical results are compared with computer simulations and experimental data.