K-space Lines Research Articles

0057: Influence of the model used to combine data acquired during multiple heart beats on temporal resolution of cardiac Cine MRI: is high temporal resolution achievable with children and young adults?

Cardiac Cine-MRI combines data acquired in different cycles of variable lengths to reconstruct images. This study aims at assessing the temporal misalignment induced by different methods for combining k-space lines. The durations of the six cardiac stages of 306 cardiac cycles were assessed with Tissue-Doppler Imaging in a population of 7 children and young adults. Different models from the literature (Chung and Feinstein), and adaptive two-stage and six-stage models were used to combine these cardiac cycles. Temporal shift between the modelized and the real positions within the cycles were recorded. The averages of the 95% confidence limits of temporal misalignments caused by two-stage models were between 20 and 30ms but during early diastole they reached 40-50 ms. Chung's model behaved slightly better than Feinstein's (26 ms vs 31 ms, p<0.001). The adaptive models significantly reduced time misalignments (22 ms vs 26 ms, p<0.001 for the 2-stage model and 18ms vs 22ms, p<0.001 for the six-stage model). There is a theoretical limitation to high temporal resolution cardiac acquisitions due to time misalignments during cine reconstruction, at least within a pediatric population. Simple personalized adaptive cardiac models may reduce this limitation. Abstract 0057 - Figure: Time missallignments when combining cardiac cycles

Lipid elimination with an echo‐shifting N/2‐ghost acquisition (LEENA) MRI

The Dixon techniques provide uniform water-fat separation but require multiple image sets, which extend the overall acquisition time. Here, an alternative rapid single acquisition method, lipid elimination with an echo-shifting N/2-ghost acquisition (LEENA), was introduced. The LEENA method utilized a fast imaging with steady-state free precession sequence to obtain a single k-space dataset in which successive k-space lines are acquired to allow the fat magnetization to precess 180°. The LEENA data were then unghosted using either image-domain (LEENA-S) or k-space domain (LEENA-G) parallel imaging techniques to reconstruct water-only and fat-only images. An off-resonance correction technique was incorporated to improve the uniformity of the water-fat separation. Uniform water-fat separation was achieved for both the LEENA-S and LEENA-G methods for phantom and human body and leg imaging applications at 1.5T and 3T. The resultant water and fat images were qualitatively similar to conventional 2-point Dixon and fat-suppressed images. The LEENA-S and LEENA-G methods provide uniform water and fat images from a single MRI acquisition. These straightforward methods can be adapted to 1.5T and 3T clinical MRI scanners and provide comparable fat/water separation with conventional 2-point Dixon and fat-suppression techniques.

An improved method for retrospective motion correction in quantitative T2* mapping

A new method for motion correction of T2*-weighted data and resulting quantitative T2* maps is presented. For this method, additional data sets with a reduced number of phase encoding steps covering the k-space centre are acquired. Motion correction is based on a 3-step procedure: (1) calculation of improved input data sets with reduced artefact levels from the original data, (2) creation of a target data set free of movement artefacts on the basis of the improved input data sets, and (3) fitting of original data to the target data set, yielding an optimum combination of acquired k-space data which suppresses lines affected by movement. The method was tested on healthy subjects performing pre-trained movement. Motion correction was successful unless the same k-space line was affected by movement in all data sets acquired on a specific subject. The method was applied to patients suffering from subarachnoid haemorrhage (group 1) or tumours (group 2) with accompanying edema in the brain. Motion correction improved the interpretability of T2*-weighted patient data and resulting quantitative T2* maps considerably by allowing a clear delineation between ventricle and edema and a clear localisation of haemorrhage (group 1) or a clear delineation of tumour accompanying edema (group 2) which was not possible in data affected by movement.

Accelerated myocardial T1 mapping using SMART1Map and compressed sensing with temporal PCA

Background Myocardial T1 mapping has been proposed for characterizing a wide range of diffuse and focal fibrotic pathologies. Due to the demands of requiring multiple sample points along the signal relaxation curve in a breath-held ECG-gated scan, current T1 mapping methods have focused on single-shot acquisitions. Collecting all data in a quiescent time window has limited T1 mapping to a reduced number of k-space lines, often with parallel imaging. These compromises lead to low spatial resolution and/or low signal-to-noise (SNR) images. This work explores the initial feasibility of accelerating the recently proposed SMART1Map method[1] using Compressed Sensing with temporal Principle Component Analysis (CS-tPCA)[2] to obtain higher spatial resolution and higher SNR T1 maps.

Towards a new method for cardiac tissue velocity measurements using MRI, comparison with echocardiography

Methods Eighteen healthy volunteers underwent cardiac MRI examination on a Signa HDxt 3T scanner (General Electric, Waukesha). A phase contrast sequence modified to acquire only the central k-space line was used in a midventricular small axis scan plane with through-plane velocity encoding during 128 heartbeats. The frequency encoding direction was angulated to encompass the projection of the left ventricular outflow tract plane (Figure 1). Relevant parameters were: 256-point frequency k-space size, 35 cm FOV, 1 view per segment, 150 cm/s velocity encoding and 7 ms TR. Data were Fourier transformed into 1D+t image space, subjected to SVD decomposition for automatic spatial segmentation and finally automatic temporal domain peak detection. This process lead to the measurement of four durations for each heartbeat (RR): IsoVolumic Contraction (IVC), Systolic Ejection (SE), IsoVolumic Relaxation (IVR), and Diastolic Diastasis (DD). For every cardiac cycle, the results of the automatic detection were quality controlled by a physician who could either discard the data, correct or accept it. Myocardiac Performance Index (MPI) [4] was computed as (IVC+IVR)/SE as well as normalized systole duration (IVC+SE)/RR and normalized diastasis duration DD/RR. These values were checked

Frequency correction method for improved spatial correlation of hyperpolarized 13C metabolites and anatomy

Blip-reversed echo-planar imaging (EPI) is investigated as a method for measuring and correcting the spatial shifts that occur due to bulk frequency offsets in (13)C metabolic imaging in vivo. By reversing the k-space trajectory for every other time point, the direction of the spatial shift for a given frequency is reversed. Here, mutual information is used to find the 'best' alignment between images and thereby measure the frequency offset. Time-resolved 3D images of pyruvate/lactate/urea were acquired with 5 s temporal resolution over a 1 min duration in rats (N = 6). For each rat, a second injection was performed with the demodulation frequency purposely mis-set by +35 Hz, to test the correction for erroneous shifts in the images. Overall, the shift induced by the 35 Hz frequency offset was 5.9 ± 0.6 mm (mean ± standard deviation). This agrees well with the expected 5.7 mm shift based on the 2.02 ms delay between k-space lines (giving 30.9 Hz per pixel). The 0.6 mm standard deviation in the correction corresponds to a frequency-detection accuracy of 4 Hz. A method was presented for ensuring the spatial registration between (13)C metabolic images and conventional anatomical images when long echo-planar readouts are used. The frequency correction method was shown to have an accuracy of 4 Hz. Summing the spatially corrected frames gave a signal-to-noise ratio (SNR) improvement factor of 2 or greater, compared with the highest single frame.

ViP MRI: virtual phantom magnetic resonance imaging

The ability to generate reference signals is of great benefit for quantitation of the magnetic resonance (MR) signal. The aim of the present study was to implement a dedicated experimental set-up to generate MR images of virtual phantoms. Virtual phantoms of a given shape and signal intensity were designed and the k-space representation was generated. A waveform generator converted the k-space lines into a radiofrequency (RF) signal that was transmitted to the MR scanner bore by a dedicated RF coil. The k-space lines of the virtual phantom were played line-by-line in synchronization with the magnetic resonance imaging data acquisition. Virtual phantoms of complex patterns were reproduced well in MR images without the presence of artifacts. Time-series measurements showed a coefficient of variation below 1% for the signal intensity of the virtual phantoms. An excellent linearity (coefficient of determination r (2) = 0.997 as assessed by linear regression) was observed in the signal intensity of virtual phantoms. Virtual phantoms represent an attractive alternative to physical phantoms for providing a reference signal. MR images of virtual phantoms were here generated using a stand-alone, independent unit that can be employed with MR scanners from different vendors.

Highly accelerated projection imaging with coil sensitivity encoding for rapid MRI

Rapid magnetic resonance imaging (MRI) acquisition is typically achieved by acquiring all or most lines of k-space after one radio frequency (RF) excitation. Parallel imaging techniques can further accelerate data acquisition by acquiring fewer phase-encoded lines and utilizing the spatial sensitivity information of the RF coil arrays. The goal of this study was to develop a new MRI data acquisition and reconstruction technique that is capable of reconstructing a two-dimensional (2D) image using highly undersampled k-space data without any special hardware. Such a technique would be very efficient, as it would significantly reduce the time wasted during multiple RF excitations or phase encoding and gradient switching periods. The essence of this new technique is to densely sample a small number of projections, which should be acquired at an angle other than 0° or multiples of 45°. This results in multiple rays passing through a voxel and provides new and independent measurements for each voxel. Then the images are reconstructed using the unique information coming from these projections combined with RF coil sensitivity profiles. The feasibility of this new technique was investigated with realistic simulations and experimental studies using a phantom and compared with conventional nonuniform fast Fourier transform technique. Eigenvalue analysis and error calculations were conducted to find optimal projection angles and minimum requirements for dense sampling. Reconstruction of 64 × 64 images was done using a single projection from simulated data under different noise levels. Simulated reconstruction was also tested with two projections to assess the improvement. Experimental phantom images were reconstructed at higher resolution using 4, 8, and 16 projections. Cross-sectional profiles illustrate that the new technique resolved compartment boundaries clearly. Simulations demonstrated that only a single k-space line might be sufficient to reconstruct a 2D image using this new technique. Experimental results showed that this is a promising new technique for fast imaging. Using the information from the simulations and fast imaging parameters published in the literature, it could be predicted that a two-dimensional image could be acquired in about 10 ms. One of the major advantages of this new technique is that it does not require any additional hardware and can be implemented on a conventional scanner with an eight-channel coil.

Temporally constrained respiratory gating improves continuously moving table MRI during free breathing

To evaluate a novel breathing motion correction algorithm for continuously moving table magnetic resonance imaging (CMT-MRI) that optimizes motion consistency in a fixed time span. In 22 patients CMT-MRI was performed during free breathing. During a preparatory phase (constant) or continuously during the scan (adaptive) gating thresholds were computed from breathing states that should allow for motion consistent k-space sampling. After data from a first k-space traversal was acquired irrespective of breathing motion, subsequently k-space lines with discordant breathing states were reacquired below the gating threshold. Time constraints of CMT-MRI were respected, because a fixed time was allocated for reacquisition. Image quality and lesion depiction were evaluated on images reconstructed from the first traversal and motion-corrected images. Compared to constant thresholds, gating with adaptive thresholds led to a higher number of reacquired k-space lines (60.1%/41.7%) and a larger fraction of motion consistent final k-space lines (96.6%/78.8%). Adaptive gating induced a significant increase in image quality for all regions affected by breathing motion. Only one of 22 lesions was not depicted on the adaptively corrected images, whereas 15 were readily appreciable. Temporally constrained respiratory gating with adaptive thresholds allows for fully sampled, motion-corrected CMT-MRI acquisitions during free breathing.

Utility of respiratory‐navigator‐rejected k‐space lines for improved signal‐to‐noise ratio in three‐dimensional cardiac MR

To develop and evaluate a technique that uses the k-space lines rejected by prospective respiratory navigator (NAV) to improve the signal-to-noise ratio (SNR) without increasing the scan time. In conventional image reconstruction, the motion-corrupted k-space lines rejected by the NAV are not used. In this study, a set of translational motion parameters for the NAV-rejected lines and a phase-corrected average for the k-space line are estimated jointly using a maximum-likelihood approach and the information from the corresponding accepted k-space lines. Left coronary artery images were acquired in 10 healthy adult subjects, and the proposed approach incorporating the NAV-rejected lines was compared with the conventional dataset with NAV-accepted lines only, as well as a simple average of all k-space lines, in terms of SNR, normalized vessel sharpness and qualitative image scores on a four-point scale (1 = poor, 4 = excellent). Late gadolinium enhancement images of the left atrium were also acquired in 21 patients with atrial fibrillation pre- or post-pulmonary vein isolation. Images reconstructed with the proposed, conventional, and simple averaging methods were compared in terms of SNR, and subjective image quality on a four-point scale. For coronary MRI, there was a significant improvement in SNR with the proposed technique, but no significant difference in normalized vessel sharpness or qualitative image scores were observed with respect to the conventional method. Simple averaging resulted in an SNR gain, but significant loss in vessel sharpness and image quality. For late gadolinium enhancement, there was a significant increase in SNR, but no significant differences were observed in subjective image quality scores between the proposed and conventional methods. There was an SNR gain, but image quality loss for simple averaging, when compared with the conventional technique. In both coronary MRI and late gadolinium enhancement, the SNR gain of the proposed method was not significantly different than the maximum theoretical SNR gain. The proposed technique improves SNR using the additional information from NAV-rejected k-space lines, while providing similar image quality to standard reconstruction using motion-free k-space data only, with no increase in scan time.

EFFECT OF PHASE-ENCODING REDUCTION ON GEOMETRIC DISTORTION AND BOLD SIGNAL CHANGES IN FMRI

Introduction Echo-planar imaging (EPI) is a group of fast data acquisition methods commonly used in fMRI studies. It acquires multiple image lines in k-space after a single excitation, which leads to a very short scan time. A well-known problem with EPI is that it is more sensitive to distortions due to the used encoding scheme. Source of distortion is inhomogeneity in the static B0 field that causes more geometric distortion in phase encoding direction. This inhomogeneity is induced mainly by the magnetic susceptibility differences between various structures within the object placed inside the scanner, often at air-tissue or bone-tissue interfaces. Methods of reducing EPI distortion are mainly based on decreasing steps of the phase encoding. Reducing steps of phase encoding can be applied by reducing field of view, slice thickness, and/or the use of parallel acquisition technique. Materials and Methods We obtained three data acquisitions with different FOVs including: conventional low resolution, conventional high resolution, and zoomed high resolution EPIs. Moreover we used SENSE technique for phase encoding reduction. All experiments were carried out on three Tesla scanners (Siemens, TIM, and Germany) equipped with 12 channel head coil. Ten subjects participated in the experiments. Results The data were processed by FSL software and were evaluated by ANOVA. Distortion was assessed by obtaining low displacement voxels map, and calculated from a field map image. Conclusion We showed that image distortion can be reduced by decreasing slice thickness and phase encoding steps. Distortion reduction in zoomed technique resulted the lowest level, but at the cost of signal-to-noise loss. Moreover, the SENSE technique was shown to decrease the amount of image distortion, efficiently.

High temporal resolution in vivo blood oximetry via projection-basedT2measurement

Measuring venous oxygen saturation (HbO2) in large blood vessels can provide important information about oxygen delivery and its consumption in vital organs. Quantification of blood's T2 value via MR can be utilized to determine HbO2 noninvasively. We propose a fast method for in vivo blood T2 quantification via computing the complex difference of velocity-encoded projections. As blood flows continuously, its signal can be robustly isolated from the surrounding tissue by computing the complex difference of two central k-space lines with different velocity encodings. This resultant signal can then be measured as a function of echo time for rapidly quantifying T2 of blood. We applied the method to quantify HbO2 in three cerebral veins at rest and in one of the veins in response to hypercapnia. Average HbO2 measurements in superior sagittal sinus (SSS), straight sinus and internal jugular vein in the group were 63 ± 3%, 68 ± 4% and 65 ± 4%, respectively. Average HbO2 values in SSS during baseline, hypercapnia, and recovery were 63 ± 2%, 79 ± 5%, and 61 ± 3%, respectively. When compared with standard T2 quantification techniques, the proposed method is fast, reliable, and robust against partial volume effects.

Center‐out echo‐planar spectroscopic imaging with correction of gradient‐echo phase and time shifts

A procedure to prevent the formation of image and spectral Nyquist ghosts in echo-planar spectroscopic imaging is introduced. It is based on a novel Cartesian center-out echo-planar spectroscopic imaging trajectory, referred to as EPSICO, and combined with a correction of the gradient-echo phase and time shifts. Processing of homogenous sets of forward and reflected echoes is no longer necessary, resulting in an optimized spectral width. The proposed center-out trajectory passively prevents the formation of Nyquist ghosts by privileging the acquisition of the center k-space line with forward echoes at the beginning of an echo-planar spectroscopic imaging dwell time and by ensuring that all k-space lines and their respective complex conjugates are acquired at equal time intervals. With the proposed procedure, concentrations of N-acetyl aspartate, creatine, choline, glutamate, and myo-inositol were reliably determined in human white matter at 3 T.

Keyhole Imaging기법을 적용한 위상대조도 자기공명 혈관조영기법

Phase Contrast MR Angiography(PC MRA) is excellent MRA technique for measuring the velocity of vessels in the human body. PC MRA need to at least four images for angiogram reconstruction and it caused longer scan time. Therefore, we used keyhole imaging combined PC MRA to reduce the scan time. However, keyhole imaging can lead the erroneous effects as loss of phase information or frequency discontinuous. In this study, we applied the keyhole imaging combined 2D PC MRA for improving the temporal resolution and also measured the velocity to evaluate the accuracy of phase information. We used 0.32T MRI scanner(Magfinder II, Scimedix, Korea). Using the 2D PC MRA pulse sequence, the vascular images for a human brain targeted on the Superior Sagittal Sinus(SSS) were obtained. We applied tukey window function for keyhole images to minimize the ringing artifact and erroneous factors that are induced frequency discontinuous and phase information loss. We also applied zero-padded algorithm to peripheral missing k-space lines to compare keyhole imaging results and the artifact power(AP) value was measured on the complex difference images to validate the image quality. Consider as based on our results, heavy image distortions and artifacts were shown until using at least 50% keyhole factor. Using above the 50% keyhole factors are shown well reconstructed and matched for magnitude images and velocity information measurements. In conclusion, we confirmed the image quality and velocity information of keyhole technique combined 2D PC MRA. Especially, measured velocity information through the keyhole imaging combination was similar to the velocity information of full sampled k-space image despite of frequency discontinuous and phase information loss in the keyhole imaging reconstruction process. Consequently, the keyhole imaging combined 2D PC MRA will give some clinical usefulness and advantages as improving the temporal resolution and measuring the velocity information via selecting the appropriate keyhole factor at low tesla MRI system.

Improving diffusion MRI using simultaneous multi-slice echo planar imaging

In diffusion MRI, simultaneous multi-slice single-shot EPI acquisitions have the potential to increase the number of diffusion directions obtained per unit time, allowing more diffusion encoding in high angular resolution diffusion imaging (HARDI) acquisitions. Nonetheless, unaliasing simultaneously acquired, closely spaced slices with parallel imaging methods can be difficult, leading to high g-factor penalties (i.e., lower SNR). The CAIPIRINHA technique was developed to reduce the g-factor in simultaneous multi-slice acquisitions by introducing inter-slice image shifts and thus increase the distance between aliased voxels. Because the CAIPIRINHA technique achieved this by controlling the phase of the RF excitations for each line of k-space, it is not directly applicable to single-shot EPI employed in conventional diffusion imaging. We adopt a recent gradient encoding method, which we termed "blipped-CAIPI", to create the image shifts needed to apply CAIPIRINHA to EPI. Here, we use pseudo-multiple replica SNR and bootstrapping metrics to assess the performance of the blipped-CAIPI method in 3× simultaneous multi-slice diffusion studies. Further, we introduce a novel image reconstruction method to reduce detrimental ghosting artifacts in these acquisitions. We show that data acquisition times for Q-ball and diffusion spectrum imaging (DSI) can be reduced 3-fold with a minor loss in SNR and with similar diffusion results compared to conventional acquisitions.

Integrated variable projection approach (IVAPA) for parallel magnetic resonance imaging

Parallel magnetic resonance imaging (pMRI) is a fast method which requires algorithms for the reconstructing image from a small number of measured k-space lines. The accurate estimation of the coil sensitivity functions is still a challenging problem in parallel imaging. The joint estimation of the coil sensitivity functions and the desired image has recently been proposed to improve the situation by iteratively optimizing both the coil sensitivity functions and the image reconstruction. It regards both the coil sensitivities and the desired images as unknowns to be solved for jointly. In this paper, we propose an integrated variable projection approach (IVAPA) for pMRI, which integrates two individual processing steps (coil sensitivity estimation and image reconstruction) into a single processing step to improve the accuracy of the coil sensitivity estimation using the variable projection approach. The method is demonstrated to be able to give an optimal solution with considerably reduced artifacts for high reduction factors and a low number of auto-calibration signal (ACS) lines, and our implementation has a fast convergence rate. The performance of the proposed method is evaluated using a set of in vivo experiment data.

Using learned under-sampling pattern for increasing speed of cardiac cine MRI based on compressive sensing principles

Abstract This article presents a compressive sensing approach for reducing data acquisition time in cardiac cine magnetic resonance imaging (MRI). In cardiac cine MRI, several images are acquired throughout the cardiac cycle, each of which is reconstructed from the raw data acquired in the Fourier transform domain, traditionally called k-space. In the proposed approach, a majority, e.g., 62.5%, of the k-space lines (trajectories) are acquired at the odd time points and a minority, e.g., 37.5%, of the k-space lines are acquired at the even time points of the cardiac cycle. Optimal data acquisition at the even time points is learned from the data acquired at the odd time points. To this end, statistical features of the k-space data at the odd time points are clustered by fuzzy c-means and the results are considered as the states of Markov chains. The resulting data is used to train hidden Markov models and find their transition matrices. Then, the trajectories corresponding to transition matrices far from an identity matrix are selected for data acquisition. At the end, an iterative thresholding algorithm is used to reconstruct the images from the under-sampled k-space datasets. The proposed approaches for selecting the k-space trajectories and reconstructing the images generate more accurate images compared to alternative methods. The proposed under-sampling approach achieves an acceleration factor of 2 for cardiac cine MRI.

Single‐shot echo‐planar imaging with Nyquist ghost compensation: Interleaved dual echo with acceleration (IDEA) echo‐planar imaging (EPI)

Echo planar imaging (EPI) is most commonly used for blood oxygen level-dependent fMRI, owing to its sensitivity and acquisition speed. A major problem with EPI is Nyquist (N/2) ghosting, most notably at high field. EPI data are acquired under an oscillating readout gradient and hence vulnerable to gradient imperfections such as eddy current delays and off-resonance effects, as these cause inconsistencies between odd and even k-space lines after time reversal. We propose a straightforward and pragmatic method herein termed "interleaved dual echo with acceleration (IDEA) EPI": two k-spaces (echoes) are acquired under the positive and negative readout lobes, respectively, by performing phase encoding blips only before alternate readout gradients. From these two k-spaces, two almost entirely ghost free images per shot can be constructed, without need for phase correction. The doubled echo train length can be compensated by parallel imaging and/or partial Fourier acquisition. The two k-spaces can either be complex averaged during reconstruction, which results in near-perfect cancellation of residual phase errors, or reconstructed into separate images. We demonstrate the efficacy of IDEA EPI and show phantom and in vivo images at both 3 T and 7 T.

Cardiac arterial spin labeling using segmented ECG‐gated Look‐Locker FAIR: Variability and repeatability in preclinical studies

MRI is important for the assessment of cardiac structure and function in preclinical studies of cardiac disease. Arterial spin labeling techniques can be used to measure perfusion noninvasively. In this study, an electrocardiogram-gated Look-Locker sequence with segmented k-space acquisition has been implemented to acquire single slice arterial spin labeling data sets in 15 min in the mouse heart. A data logger was introduced to improve data quality by: (1) allowing automated rejection of respiration-corrupted images, (2) providing additional prospective gating to improve consistency of acquisition timing, and (3) allowing the recombination of uncorrupted k-space lines from consecutive data sets to reduce respiration corruption. Finally, variability and repeatability of perfusion estimation within-session, between-session, between-animal, and between image rejection criteria were assessed in mice. The criterion used to reject images from the T(1) fit was shown to affect the perfusion estimation. These data showed that the between-animal coefficient of variability (24%) was greater than the between-session variability (17%) and within-session variability (11%). Furthermore, the magnitude of change in perfusion required to detect differences was 30% (within-session) and 55% (between-session) according to Bland-Altman repeatability analysis. These technique developments and repeatability statistics will provide a platform for future preclinical studies applying cardiac arterial spin labeling.

Improved navigator-gated motion compensation in cardiac MR using additional constraint of magnitude of motion-corrupted data

Background In conventional prospective respiratory navigator (NAV) acquisitions, 40-60% of the acquired data are discarded resulting in low efficiency and long scan times [1,2]. Compressed-sensing Motion Compensation (CosMo) has a shorter fixed scan time by acquiring the full inner k-space and estimating the NAV-rejected outer k-space lines [3]. Respiratory motion will mainly manifest itself as phase variation in the acquired k-space data. We sought to determine if the addition of the magnitude of the rejected k-space lines as a constraint in image reconstruction will improve the performance of CosMo.