Accelerated liver water T1 mapping using single-shot continuous inversion-recovery spiral imaging

PURPOSE

Liver T1 mapping techniques typically require long breath holds or long scan time in free-breathing, need correction for B 1 + inhomogeneities and process composite (water and fat) signals. The purpose of this work is to accelerate the multi-slice acquisition of liver water selective T1 (wT1) mapping in a single breath hold, improving the k-space sampling efficiency.

METHODS

The proposed continuous inversion-recovery (IR) Look-Locker methodology combines a single-shot gradient echo spiral readout, Dixon processing and a dictionary-based analysis for liver wT1 mapping at 3 T. The sequence parameters were adapted to obtain short scan times. The influence of fat, B 1 + inhomogeneities and TE on the estimation of T1 was first assessed using simulations. The proposed method was then validated in a phantom and in 10 volunteers, comparing it with MRS and the modified Look-Locker inversion-recovery (MOLLI) method. Finally, the clinical feasibility was investigated by comparing wT1 maps with clinical scans in nine patients.

RESULTS

The phantom results are in good agreement with MRS. The proposed method encodes the IR-curve for the liver wT1 estimation, is minimally sensitive to B 1 + inhomogeneities and acquires one slice in 1.2 s. The volunteer results confirmed the multi-slice capability of the proposed method, acquiring nine slices in a breath hold of 11 s. The present work shows robustness to B 1 + inhomogeneities (wT 1 , No B 1 + = 1.07 wT 1 , B 1 + - 45.63 , R 2 = 0.99 ) , good repeatability (wT 1 , 2 ° = 1 . 0 wT 1 , 1 ° - 2.14 , R 2 = 0.96 ) and is in better agreement with MRS (wT 1 = 0.92 wT 1 MRS + 103.28 , R 2 = 0.38 ) than is MOLLI (wT 1 MOLLI = 0.76 wT 1 MRS + 254.43 , R 2 = 0.44 ) . The wT1 maps in patients captured diverse lesions, thus showing their clinical feasibility.

CONCLUSION

A single-shot spiral acquisition can be combined with a continuous IR Look-Locker method to perform rapid repeatable multi-slice liver water T1 mapping at a rate of 1.2 s per slice without a B 1 + map. The proposed method is suitable for nine-slice liver clinical applications acquired in a single breath hold of 11 s.

AI-based motion artifact severity estimation in undersampled MRI allowing for selection of appropriate reconstruction models

BACKGROUND

Magnetic Resonance acquisition is a time consuming process, making it susceptible to patient motion during scanning. Even motion in the order of a millimeter can introduce severe blurring and ghosting artifacts, potentially necessitating re-acquisition. Magnetic Resonance Imaging (MRI) can be accelerated by acquiring only a fraction of k-space, combined with advanced reconstruction techniques leveraging coil sensitivity profiles and prior knowledge. Artificial intelligence (AI)-based reconstruction techniques have recently been popularized, but generally assume an ideal setting without intra-scan motion.

PURPOSE

To retrospectively detect and quantify the severity of motion artifacts in undersampled MRI data. This may prove valuable as a safety mechanism for AI-based approaches, provide useful information to the reconstruction method, or prompt for re-acquisition while the patient is still in the scanner.

METHODS

We developed a deep learning approach that detects and quantifies motion artifacts in undersampled brain MRI. We demonstrate that synthetically motion-corrupted data can be leveraged to train the convolutional neural network (CNN)-based motion artifact estimator, generalizing well to real-world data. Additionally, we leverage the motion artifact estimator by using it as a selector for a motion-robust reconstruction model in case a considerable amount of motion was detected, and a high data consistency model otherwise.

RESULTS

Training and validation were performed on 4387 and 1304 synthetically motion-corrupted images and their uncorrupted counterparts, respectively. Testing was performed on undersampled in vivo motion-corrupted data from 28 volunteers, where our model distinguished head motion from motion-free scans with 91% and 96% accuracy when trained on synthetic and on real data, respectively. It predicted a manually defined quality label ('Good', 'Medium' or 'Bad' quality) correctly in 76% and 85% of the time when trained on synthetic and real data, respectively. When used as a selector it selected the appropriate reconstruction network 93% of the time, achieving near optimal SSIM values.

CONCLUSIONS

The proposed method quantified motion artifact severity in undersampled MRI data with high accuracy, enabling real-time motion artifact detection that can help improve the safety and quality of AI-based reconstructions.

Digital reference objects for evaluating algorithm performance in MR image formation

PURPOSE

Digital Reference Objects (DROs) are mathematical phantoms that can serve as a basis for evaluating MR image quality (IQ) in an objective way. Their main purpose is to facilitate the establishment of fully automated and perfectly reproducible IQ metrics to objectively compare different algorithms in MR image formation in a standardized manner. They also allow to re-build parts of standard phantoms.

METHODS

We sample DROs directly in k-space, using analytical formulas for the continuous Fourier transform of primitive shapes. We demonstrate this DRO approach by applying a state-of-the-art CNN-based denoising algorithm that is robust to varying noise levels to noisy images of the resolution section of the well-known ACR phantom for IQ assessment, reconstructed from both measured and simulated k-space data.

RESULTS

Applying the CNN-based denoising algorithm to the measured and simulated version of the ACR phantom resolution section produced virtually identical results, as confirmed by visual and quantitative comparison.

CONCLUSIONS

DROs can help guide technology selection during the development of new algorithms in MR image formation, e.g., via deep learning. This could be an important step towards reproducible MR image formation.

High-fidelity direct contrast synthesis from magnetic resonance fingerprinting

PURPOSE

This work was aimed at proposing a supervised learning-based method that directly synthesizes contrast-weighted images from the Magnetic Resonance Fingerprinting (MRF) data without performing quantitative mapping and spin-dynamics simulations.

METHODS

To implement our direct contrast synthesis (DCS) method, we deploy a conditional generative adversarial network (GAN) framework with a multi-branch U-Net as the generator and a multilayer CNN (PatchGAN) as the discriminator. We refer to our proposed approach as N-DCSNet. The input MRF data are used to directly synthesize T1-weighted, T2-weighted, and fluid-attenuated inversion recovery (FLAIR) images through supervised training on paired MRF and target spin echo-based contrast-weighted scans. The performance of our proposed method is demonstrated on in vivo MRF scans from healthy volunteers. Quantitative metrics, including normalized root mean square error (nRMSE), peak signal-to-noise ratio (PSNR), structural similarity (SSIM), learned perceptual image patch similarity (LPIPS), and Fréchet inception distance (FID), were used to evaluate the performance of the proposed method and compare it with others.

RESULTS

In-vivo experiments demonstrated excellent image quality with respect to that of simulation-based contrast synthesis and previous DCS methods, both visually and according to quantitative metrics. We also demonstrate cases in which our trained model is able to mitigate the in-flow and spiral off-resonance artifacts typically seen in MRF reconstructions, and thus more faithfully represent conventional spin echo-based contrast-weighted images.

CONCLUSION

We present N-DCSNet to directly synthesize high-fidelity multicontrast MR images from a single MRF acquisition. This method can significantly decrease examination time. By directly training a network to generate contrast-weighted images, our method does not require any model-based simulation and therefore can avoid reconstruction errors due to dictionary matching and contrast simulation (code available at:https://github.com/mikgroup/DCSNet).

Quantitative longitudinal mapping of radiation-treated prostate cancer using MR fingerprinting with radial acquisition and subspace reconstruction

MR fingerprinting (MRF) enables fast multiparametric quantitative imaging with a single acquisition and has been shown to improve diagnosis of prostate cancer. However, most prostate MRF studies were performed with spiral acquisitions that are sensitive to B0 inhomogeneities and consequent blurring. In this work, a radial MRF acquisition with a novel subspace reconstruction technique was developed to enable fast T1/T2 mapping in the prostate in under 4 min. The subspace reconstruction exploits the extensive temporal correlations in the MRF dictionary to pre-compute a low dimensional space for the solution and thus reduce the number of radial spokes to accelerate the acquisition. Iterative reconstruction with the subspace model and additional regularization of the signal representation in the subspace is performed to minimize the number of spokes and maintain matching quality and SNR. Reconstruction accuracy was assessed using the ISMRM NIST phantom. In-vivo validation was performed on two healthy subjects and two prostate cancer patients undergoing radiation therapy. The longitudinal repeatability was quantified using the concordance correlation coefficient (CCC) in one of the healthy subjects by repeated scans over 1 year. One prostate cancer patient was scanned at three time points, before initiating therapy and following brachytherapy and external beam radiation. Changes in the T1/T2 maps obtained with the proposed method were quantified. The prostate, peripheral and transitional zones, and visible dominant lesion were delineated for each study, and the statistics and distribution of the quantitative mapping values were analyzed. Significant image quality improvements compared with standard reconstruction methods were obtained with the proposed subspace reconstruction method. A notable decrease in the spread of the T1/T2 values without biasing the estimated mean values was observed with the subspace reconstruction and agreed with reported literature values. The subspace reconstruction enabled visualization of small differences in T1/T2 values in the tumor region within the peripheral zone. Longitudinal imaging of a volunteer subject yielded CCC of 0.89 for MRF T1, and 0.81 for MRF T2 in the prostate gland. Longitudinal imaging of the prostate patient confirmed the feasibility of capturing radiation treatment related changes. This work is a proof-of-concept for a high resolution and fast quantitative mapping using golden-angle radial MRF combined with a subspace reconstruction technique for longitudinal treatment response assessment in subjects undergoing radiation treatment.

Comparing myelin-sensitive magnetic resonance imaging measures and resulting g-ratios in healthy and multiple sclerosis brains

The myelin concentration and the degree of myelination of nerve fibers can provide valuable information on the integrity of human brain tissue. Magnetic resonance imaging (MRI) of myelin-sensitive parameters can help to non-invasively evaluate demyelinating diseases such as multiple sclerosis (MS). Several different myelin-sensitive MRI methods have been proposed to determine measures of the degree of myelination, in particular the g-ratio. However, variability in underlying physical principles and different biological models influence measured myelin concentrations, and consequently g-ratio values. We therefore investigated similarities and differences between five different myelin-sensitive MRI measures and their effects on g-ratio mapping in the brains of both MS patients and healthy volunteers. We compared two different estimates of the myelin water fraction (MWF) as well as the inhomogeneous magnetization transfer ratio (ihMTR), magnetization transfer saturation (MTsat), and macromolecular tissue volume (MTV) in 13 patients with MS and 14 healthy controls. In combination with diffusion-weighted imaging, we derived g-ratio parameter maps for each of the five different myelin measures. The g-ratio values calculated from different myelin measures varied strongly, especially in MS lesions. While, compared to normal-appearing white matter, MTsat and one estimate of the MWF resulted in higher g-ratio values within lesions, ihMTR, MTV, and the second MWF estimate resulted in lower lesion g-ratio values. As myelin-sensitive measures provide rough estimates of myelin content rather than absolute myelin concentrations, resulting g-ratio values strongly depend on the utilized myelin measure and model used for g-ratio mapping. When comparing g-ratio values, it is, thus, important to utilize the same MRI methods and models or to consider methodological differences. Particular caution is necessary in pathological tissue such as MS lesions.

High-Resolution Free-Breathing Quantitative First-Pass Perfusion Cardiac MR Using Dual-Echo Dixon With Spatio-Temporal Acceleration

Introduction

To develop and test the feasibility of free-breathing (FB), high-resolution quantitative first-pass perfusion cardiac MR (FPP-CMR) using dual-echo Dixon (FOSTERS; Fat-water separation for mOtion-corrected Spatio-TEmporally accelerated myocardial peRfuSion).

Materials and Methods

FOSTERS was performed in FB using a dual-saturation single-bolus acquisition with dual-echo Dixon and a dynamically variable Cartesian k-t undersampling (8-fold) approach, with low-rank and sparsity constrained reconstruction, to achieve high-resolution FPP-CMR images. FOSTERS also included automatic in-plane motion estimation and T2* correction to obtain quantitative myocardial blood flow (MBF) maps. High-resolution (1.6 x 1.6 mm²) FB FOSTERS was evaluated in eleven patients, during rest, against standard-resolution (2.6 x 2.6 mm²) 2-fold SENSE-accelerated breath-hold (BH) FPP-CMR. In addition, MBF was computed for FOSTERS and spatial wavelet-based compressed sensing (CS) reconstruction. Two cardiologists scored the image quality (IQ) of FOSTERS, CS, and standard BH FPP-CMR images using a 4-point scale (1–4, non-diagnostic – fully diagnostic).

Results

FOSTERS produced high-quality images without dark-rim and with reduced motion-related artifacts, using an 8x accelerated FB acquisition. FOSTERS and standard BH FPP-CMR exhibited excellent IQ with an average score of 3.5 ± 0.6 and 3.4 ± 0.6 (no statistical difference, p > 0.05), respectively. CS images exhibited severe artifacts and high levels of noise, resulting in an average IQ score of 2.9 ± 0.5. MBF values obtained with FOSTERS presented a lower variance than those obtained with CS.

Discussion

FOSTERS enabled high-resolution FB FPP-CMR with MBF quantification. Combining motion correction with a low-rank and sparsity-constrained reconstruction results in excellent image quality.

A unifying view on extended phase graphs and Bloch simulations for quantitative MRI

Quantitative MRI methods and learning-based algorithms require exact forward simulations. One critical factor to correctly describe magnetization dynamics is the effect of slice-selective RF pulses. While contemporary simulation techniques correctly capture their influence, they only provide final magnetization distributions, require to be run for each parameter set separately, and make it hard to derive general theoretical conclusions and to generate a fundamental understanding of echo formation in the presence of slice-profile effects. This work aims to provide a mathematically exact framework, which is equally intuitive as extended phase graphs (EPGs), but also considers slice-profiles through their natural spatial representation. We show, through an analytical, hybrid Bloch-EPG formalism, that the spatially-resolved EPG approach allows to exactly predict the signal dependency on off-resonance, spoiling moment, microscopic dephasing, and echo time. We also demonstrate that our formalism allows to use the same phase graph to simulate both gradient-spoiled and balanced SSFP-based MR sequences. We present a derivation of the formalism and identify the connection to existing methods, i.e. slice-selective Bloch, slice-selective EPG, and the partitioned EPG. As a use case, the proposed hybrid Bloch-EPG framework is applied to MR Fingerprinting.

A pilot study of magnetic resonance fingerprinting in Parkinson's disease

Parkinson's disease (PD) affects more than six million people, but reliable MRI biomarkers with which to diagnose patients have not been established. Magnetic resonance fingerprinting (MRF) is a recent quantitative technique that can provide relaxometric maps from a single sequence. The purpose of this study is to assess the potential of MRF to identify PD in patients and their disease severity, as well as to evaluate comfort during MRF. Twenty-five PD patients and 25 matching controls underwent 3 T MRI, including an axial 2D spoiled gradient echo MRF sequence. T1 and T2 maps were generated by voxel-wise matching the measured MRF signal to a precomputed dictionary. All participants also received standard inversion recovery T1 and multi-echo T2 mapping. An ROI-based analysis of relaxation times was performed. Differences between patients and controls as well as techniques were determined by logistic regression, Spearman correlation and t-test. Patients were asked to estimate the subjective comfort of the MRF sequence. Both MRF-based T1 and T2 mapping discriminated patients from controls: T1 relaxation times differed most in cortical grey matter (PD 1337 ± 38 vs. control 1386 ± 37 ms; mean ± SD; P = .0001) and, in combination with normal-appearing white matter, enabled correct discrimination in 85.7% of cases (sensitivity 83.3%; specificity 88.0%; receiver-operating characteristic [ROC]) area under the curve [AUC] 0.87), while for T2 mapping the left putamen was the strongest classifier (40.54 ± 6.28 vs. 34.17 ± 4.96 ms; P = .0001), enabling differentiation of groups in 84.0% of all cases (sensitivity 80.0%; specificity 88.0%; ROC AUC 0.87). Relaxation time differences were not associated with disease severity. Standard mapping techniques generated significantly different relaxation time values and identified other structures as different between groups other than MRF. Twenty-three out of 25 PD patients preferred the MRF examination instead of a standard MRI. MRF-based mapping can identify PD patients with good comfort but needs further assessment regarding disease severity identification and its potential for comparability with standard mapping technique results.

Confounding factors in breast magnetic resonance fingerprinting: B 1 + , slice profile, and diffusion effects

PURPOSE

Magnetic resonance fingerprinting (MRF) offers rapid quantitative imaging but may be subject to confounding effects (CE) if these are not included in the model-based reconstruction. This study characterizes the influence of in-plane B 1 + , slice profile and diffusion effects on T₁ and T₂ estimation in the female breast at 1.5T.

METHODS

Simulations were used to predict the influence of each CE on the accuracy of MRF and to investigate the influence of electronic noise and spiral aliasing artefacts. The experimentally observed bias in regions of fibroglandular tissue (FGT) and fatty tissue (FT) was analyzed for undersampled spiral breast MRF data of 6 healthy volunteers by performing MRF reconstruction with and without a CE.

RESULTS

Theoretic analysis predicts T₁ under-/T₂ overestimation if the nominal flip angles are underestimated and inversely, T₁ under-/T₂ overestimation if omitting slice profile correction, and T₁ under-/T₂ underestimation if omitting diffusion in the signal model. Averaged over repeated signal simulations, including spiral aliasing artefacts affected precision more than accuracy. Strong in-plane B 1 + effects occurred in vivo, causing T₂ left-right inhomogeneity between both breasts. Their correction decreased the T₂ difference from 29 to 5 ms in FGT and from 29 to 9 ms in FT. Slice profile correction affected FGT T₂ most strongly, resulting in -22% smaller values. For the employed spoiler gradient strengths, diffusion did not affect the parameter maps, corresponding well with theoretic predictions.

CONCLUSION

Understanding CEs and their relative significance for an MRF sequence is important when defining an MRF signal model for accurate parameter mapping.

T1, T2, and Fat Fraction Cardiac MR Fingerprinting: Preliminary Clinical Evaluation

BACKGROUND

Dixon cardiac magnetic resonance fingerprinting (MRF) has been recently introduced to simultaneously provide water T₁ , water T₂ , and fat fraction (FF) maps.

PURPOSE

To assess Dixon cardiac MRF repeatability in healthy subjects and its clinical feasibility in a cohort of patients with cardiovascular disease.

POPULATION

T1MES phantom, water-fat phantom, 11 healthy subjects and 19 patients with suspected cardiovascular disease.

STUDY TYPE

Prospective.

FIELD STRENGTH/SEQUENCE

1.5T, inversion recovery spin echo (IRSE), multiecho spin echo (MESE), modified Look-Locker inversion recovery (MOLLI), T₂ gradient spin echo (T₂ -GRASE), 6-echo gradient rewound echo (GRE), and Dixon cardiac MRF.

ASSESSMENT

Dixon cardiac MRF precision was assessed through repeated scans against conventional MOLLI, T₂ -GRASE, and PDFF in phantom and 11 healthy subjects. Dixon cardiac MRF native T₁ , T₂ , FF, postcontrast T₁ and synthetic extracellular volume (ECV) maps were assessed in 19 patients in comparison to conventional sequences. Measurements in patients were performed in the septum and in late gadolinium enhanced (LGE) areas and assessed using mean value distributions, correlation, and Bland-Altman plots. Image quality and diagnostic confidence were assessed by three experts using 5-point scoring scales.

STATISTICAL TESTS

Paired Wilcoxon rank signed test and paired t-tests were applied. Statistical significance was indicated by *(P < 0.05).

RESULTS

Dixon cardiac MRF showed good overall precision in phantom and in vivo. Septal average repeatability was ~23 msec for T₁ , ~2.2 msec for T₂ , and ~1% for FF. Biases in healthy subjects/patients were measured at +37 msec*/+60 msec* and -8.8 msec*/-8 msec* when compared to MOLLI and T₂ -GRASE, respectively. No statistically significant differences in postcontrast T₁ (P = 0.17) and synthetic ECV (P = 0.19) measurements were observed in patients.

DATA CONCLUSION

Dixon cardiac MRF attained good overall precision in phantom and healthy subjects, while providing coregistered T₁ , T₂ , and fat fraction maps in a single breath-hold scan with similar or better image quality than conventional methods in patients.

LEVEL OF EVIDENCE

2.

TECHNICAL EFFICACY STAGE

2.

3D free-breathing cardiac magnetic resonance fingerprinting

PURPOSE

To develop a novel respiratory motion compensated three-dimensional (3D) cardiac magnetic resonance fingerprinting (cMRF) approach for whole-heart myocardial T₁ and T₂ mapping from a free-breathing scan.

METHODS

Two-dimensional (2D) cMRF has been recently proposed for simultaneous, co-registered T₁ and T₂ mapping from a breath-hold scan; however, coverage is limited. Here we propose a novel respiratory motion compensated 3D cMRF approach for whole-heart myocardial T₁ and T₂ tissue characterization from a free-breathing scan. Variable inversion recovery and T₂ preparation modules are used for parametric encoding, respiratory bellows driven localized autofocus is proposed for beat-to-beat translation motion correction and a subspace regularized reconstruction is employed to accelerate the scan. The proposed 3D cMRF approach was evaluated in a standardized T₁ /T₂ phantom in comparison with reference spin echo values and in 10 healthy subjects in comparison with standard 2D MOLLI, SASHA and T₂ -GraSE mapping techniques at 1.5 T.

RESULTS

3D cMRF T₁ and T₂ measurements were generally in good agreement with reference spin echo values in the phantom experiments, with relative errors of 2.9% and 3.8% for T₁ and T₂ (T₂ < 100 ms), respectively. in vivo left ventricle (LV) myocardial T₁ values were 1054 ± 19 ms for MOLLI, 1146 ± 20 ms for SASHA and 1093 ± 24 ms for the proposed 3D cMRF; corresponding T₂ values were 51.8 ± 1.6 ms for T2-GraSE and 44.6 ± 2.0 ms for 3D cMRF. LV coefficients of variation were 7.6 ± 1.6% for MOLLI, 12.1 ± 2.7% for SASHA and 5.8 ± 0.8% for 3D cMRF T₁ , and 10.5 ± 1.4% for T2-GraSE and 11.7 ± 1.6% for 3D cMRF T₂ .

CONCLUSION

The proposed 3D cMRF can provide whole-heart, simultaneous and co-registered T₁ and T₂ maps with accuracy and precision comparable to those of clinical standards in a single free-breathing scan of about 7 min.

SURE-based automatic parameter selection for ESPIRiT calibration

PURPOSE

ESPIRiT is a parallel imaging method that estimates coil sensitivity maps from the auto-calibration region (ACS). This requires choosing several parameters for the optimal map estimation. While fairly robust to these parameter choices, occasionally, poor selection can result in reduced performance. The purpose of this work is to automatically select parameters in ESPIRiT for more robust and consistent performance across a variety of exams.

METHODS

By viewing ESPIRiT as a denoiser, Stein's unbiased risk estimate (SURE) is leveraged to automatically optimize parameter selection in a data-driven manner. The optimum parameters corresponding to the minimum true squared error, minimum SURE as derived from densely sampled, high-resolution, and non-accelerated data and minimum SURE as derived from ACS are compared using simulation experiments. To avoid optimizing the rank of ESPIRiT's auto-calibrating matrix (one of the parameters), a heuristic derived from SURE-based singular value thresholding is also proposed.

RESULTS

Simulations show SURE derived from the densely sampled, high-resolution, and non-accelerated data to be an accurate estimator of the true mean squared error, enabling automatic parameter selection. The parameters that minimize SURE as derived from ACS correspond well to the optimal parameters. The soft-threshold heuristic improves computational efficiency while providing similar results to an exhaustive search. In-vivo experiments verify the reliability of this method.

CONCLUSIONS

Using SURE to determine ESPIRiT parameters allows for automatic parameter selections. In-vivo results are consistent with simulation and theoretical results.

Free-running cardiac magnetic resonance fingerprinting: Joint T1/T2 map and Cine imaging

PURPOSE

To develop and evaluate a novel non-ECG triggered 2D magnetic resonance fingerprinting (MRF) sequence allowing for simultaneous myocardial T₁ and T₂ mapping and cardiac Cine imaging.

METHODS

Cardiac MRF (cMRF) has been recently proposed to provide joint T₁/T₂ myocardial mapping by triggering the acquisition to mid-diastole and relying on a subject-dependent dictionary of MR signal evolutions to generate the maps. In this work, we propose a novel "free-running" (non-ECG triggered) cMRF framework for simultaneous myocardial T₁ and T₂ mapping and cardiac Cine imaging in a single scan. Free-running cMRF is based on a transient state bSSFP acquisition with tiny golden angle radial readouts, varying flip angle and multiple adiabatic inversion pulses. The acquired data is retrospectively gated into several cardiac phases, which are reconstructed with an approach that combines parallel imaging, low rank modelling and patch-based high-order tensor regularization. Free-running cMRF was evaluated in a standardized phantom and ten healthy subjects. Comparison with reference spin-echo, MOLLI, SASHA, T₂-GRASE and Cine was performed.

RESULTS

T₁ and T₂ values obtained with the proposed approach were in good agreement with reference phantom values (ICC(A,1) > 0.99). Reported values for myocardium septum T₁ were 1043 ± 48 ms, 1150 ± 100 ms and 1160 ± 79 ms for MOLLI, SASHA and free-running cMRF respectively and for T₂ of 51.7 ± 4.1 ms and 44.6 ± 4.1 ms for T₂-GRASE and free-running cMRF respectively. Good agreement was observed between free-running cMRF and conventional Cine 2D ejection fraction (bias = -0.83%).

CONCLUSION

The proposed free-running cardiac MRF approach allows for simultaneous assessment of myocardial T₁ and T₂ and Cine imaging in a single scan.

Fast multi-component analysis using a joint sparsity constraint for MR fingerprinting

PURPOSE

To develop an efficient algorithm for multi-component analysis of magnetic resonance fingerprinting (MRF) data without making a priori assumptions about the exact number of tissues or their relaxation properties.

METHODS

Different tissues or components within a voxel are potentially separable in MRF because of their distinct signal evolutions. The observed signal evolution in each voxel can be described as a linear combination of the signals for each component with a non-negative weight. An assumption that only a small number of components are present in the measured field of view is usually imposed in the interpretation of multi-component data. In this work, a joint sparsity constraint is introduced to utilize this additional prior knowledge in the multi-component analysis of MRF data. A new algorithm combining joint sparsity and non-negativity constraints is proposed and compared to state-of-the-art multi-component MRF approaches in simulations and brain MRF scans of 11 healthy volunteers.

RESULTS

Simulations and in vivo measurements show reduced noise in the estimated tissue fraction maps compared to previously proposed methods. Applying the proposed algorithm to the brain data resulted in 4 or 5 components, which could be attributed to different brain structures, consistent with previous multi-component MRF publications.

CONCLUSIONS

The proposed algorithm is faster than previously proposed methods for multi-component MRF and the simulations suggest improved accuracy and precision of the estimated weights. The results are easier to interpret compared to voxel-wise methods, which combined with the improved speed is an important step toward clinical evaluation of multi-component MRF.

High Acceleration Three-Dimensional T1-Weighted Dual Echo Dixon Hepatobiliary Phase Imaging Using Compressed Sensing-Sensitivity Encoding: Comparison of Image Quality and Solid Lesion Detectability with the Standard T1-Weighted Sequence

Objective

To compare a high acceleration three-dimensional (3D) T1-weighted gradient-recalled-echo (GRE) sequence using the combined compressed sensing (CS)-sensitivity encoding (SENSE) method with a conventional 3D GRE sequence using SENSE, with respect to image quality and detectability of solid focal liver lesions (FLLs) in the hepatobiliary phase (HBP) of gadoxetic acid-enhanced liver MRI.

Materials and Methods

A total of 217 patients with gadoxetic acid-enhanced liver MRI at 3T (54 in the preliminary study and 163 in the main study) were retrospectively included. In the main study, HBP imaging was done twice using the standard mDixon-3D-GRE technique with SENSE (acceleration factor [AF]: 2.8, standard mDixon-GRE) and the high acceleration mDixon-3D GRE technique using the combined CS-SENSE technique (CS-SENSE mDixon-GRE). Two abdominal radiologists assessed the two MRI data sets for image quality in consensus. Three other abdominal radiologists independently assessed the diagnostic performance of each data set and its ability to detect solid FLLs in 117 patients with 193 solid nodules and compared them using jackknife alternative free-response receiver operating characteristics (JAFROC).

Results

There was no significant difference in the overall image quality. CS-SENSE mDixon-GRE showed higher image noise, but lesser motion artifact levels compared with the standard mDixon-GRE (all p < 0.05). In terms of lesion detection, reader-averaged figures-of-merit estimated with JAFROC was 0.918 for standard mDixon-GRE, and 0.953 for CS-SENSE mDixon-GRE (p = 0.142). The non-inferiority of CS-SENSE mDixon-GRE over standard mDixon-GRE was confirmed (difference: 0.064 [−0.012, 0.081]).

Conclusion

The CS-SENSE mDixon-GRE HBP sequence provided comparable overall image quality and non-inferior solid FFL detectability compared with the standard mDixon-GRE sequence, with reduced acquisition time.

Water-fat Dixon cardiac magnetic resonance fingerprinting

Purpose

Cardiac magnetic resonance fingerprinting (cMRF) has been recently introduced to simultaneously provide T₁, T₂, and M₀ maps. Here, we develop a 3‐point Dixon‐cMRF approach to enable simultaneous water specific T₁, T₂, and M₀ mapping of the heart and fat fraction (FF) estimation in a single breath‐hold scan.

Methods

Dixon‐cMRF is achieved by combining cMRF with several innovations that were previously introduced for other applications, including a 3‐echo GRE acquisition with golden angle radial readout and a high‐dimensional low‐rank tensor constrained reconstruction to recover the highly undersampled time series images for each echo. Water–fat separation of the Dixon‐cMRF time series is performed to allow for water‐ and fat‐specific T₁, T₂, and M₀ estimation, whereas FF estimation is extracted from the M₀ maps. Dixon‐cMRF was evaluated in a standardized T₁–T₂ phantom, in a water–fat phantom, and in healthy subjects in comparison to current clinical standards: MOLLI, SASHA, T₂‐GRASE, and 6‐point Dixon proton density FF (PDFF) mapping.

Results

Dixon‐cMRF water T₁ and T₂ maps showed good agreement with reference T₁ and T₂ mapping techniques (R² > 0.99 and maximum normalized RMSE ~5%) in a standardized phantom. Good agreement was also observed between Dixon‐cMRF FF and reference PDFF (R² > 0.99) and between Dixon‐cMRF water T₁ and T₂ and water selective T₁ and T₂ maps (R² > 0.99) in a water–fat phantom. In vivo Dixon‐cMRF water T₁ values were in good agreement with MOLLI and water T₂ values were slightly underestimated when compared to T₂‐GRASE. Average myocardium septal T₁ values were 1129 ± 38 ms, 1026 ± 28 ms, and 1045 ± 32 ms for SASHA, MOLLI, and the proposed water Dixon‐cMRF. Average T₂ values were 51.7 ± 2.2 ms and 42.8 ± 2.6 ms for T₂‐GRASE and water Dixon‐cMRF, respectively. Dixon‐cMRF FF maps showed good agreement with in vivo PDFF measurements (R² > 0.98) and average FF in the septum was measured at 1.3%.

Conclusion

The proposed Dixon‐cMRF allows to simultaneously quantify myocardial water T₁, water T₂, and FF in a single breath‐hold scan, enabling multi‐parametric T₁, T₂, and fat characterization. Moreover, reduced T₁ and T₂ quantification bias caused by water–fat partial volume was demonstrated in phantom experiments.

MR fingerprinting as a diagnostic tool in patients with frontotemporal lobe degeneration: A pilot study

Several very rare forms of dementia are associated with characteristic focal atrophy predominantly of the frontal and/or temporal lobes and currently lack imaging solutions to monitor disease. Magnetic resonance fingerprinting (MRF) is a recently developed technique providing quantitative relaxivity maps and images with various tissue contrasts out of a single sequence acquisition. This pilot study explores the utility of MRF-based T1 and T2 mapping to discover focal differences in relaxation times between patients with frontotemporal lobe degenerative dementia and healthy controls. 8 patients and 30 healthy controls underwent a 3 T MRI including an axial 2D spoiled gradient echo MRF sequence. T1 and T2 relaxation maps were generated based on an extended phase graphs algorithm-founded dictionary involving inner product pattern matching. A region of interest (ROI)-based analysis of T1 and T2 relaxation times was performed with FSL and ITK-SNAP. Depending on the brain region analyzed, T1 relaxation times were up to 10.28% longer in patients than in controls reaching significant differences in cortical gray matter (P = .047) and global white matter (P = .023) as well as in both hippocampi (P = .001 left; P = .027 right). T2 relaxation times were similarly longer in the hippocampus by up to 19.18% in patients compared with controls. The clinically most affected patient had the most control-deviant relaxation times. There was a strong correlation of T1 relaxation time in the amygdala with duration of the clinically manifest disease (Spearman Rho = .94; P = .001) and of T1 relaxation times in the left hippocampus with disease severity (Rho = .90, P = .002). In conclusion, MRF-based relaxometry is a promising and time-saving new MRI tool to study focal cerebral alterations and identify patients with frontotemporal lobe degeneration. To validate the results of this pilot study, MRF is worth further exploration as a diagnostic tool in neurodegenerative diseases.

Dictionary-based electric properties tomography

PURPOSE

To develop and validate a new algorithm called "dictionary-based electric properties tomography" (dbEPT) for deriving tissue electric properties from measured B₁ maps.

METHODS

Inspired by Magnetic Resonance fingerprinting, dbEPT uses a dictionary of local patterns ("atoms") of B₁ maps and corresponding electric properties distributions, derived from electromagnetic field simulations. For reconstruction, a pattern from a measured B₁ map is compared with the B₁ atoms of the dictionary. The B₁ atom showing the best match with the measured B₁ pattern yields the optimum electric properties pattern that is chosen for reconstruction. Matching was performed through machine learning algorithms. Two dictionaries, using transmit and transceive phases, were evaluated. The spatial distribution of local matching distance between optimal atom and measured pattern yielded a reconstruction reliability map. The method was applied to reconstruct conductivity of 4 volunteers' brains. A conventional, Helmholtz-based Electric properties tomography (EPT) reconstruction was performed for reference. Noise performance was studied through phantom simulations.

RESULTS

Quantitative values of conductivity agree with literature values. Results of the 2 dictionaries exhibit only minor differences. Somewhat larger differences are visible between dbEPT and Helmholtz-based EPT. Quantified by the correlation between conductivity and anatomic images, dbEPT depicts brain details more clearly than Helmholtz-based EPT. Matching distance is minimal in homogeneous brain ventricles and increases with tissue heterogeneity. Central processing unit time was approximately 2 minutes per dictionary training and 3 minutes per brain conductivity reconstruction using standard hardware equipment.

CONCLUSION

A new, dictionary-based approach for reconstructing electric properties is presented. Its conductivity reconstruction is able to overcome the EPT transceive-phase problem.

Towards predicting the encoding capability of MR fingerprinting sequences

Sequence optimization and appropriate sequence selection is still an unmet need in magnetic resonance fingerprinting (MRF). The main challenge in MRF sequence design is the lack of an appropriate measure of the sequence's encoding capability. To find such a measure, three different candidates for judging the encoding capability have been investigated: local and global dot-product-based measures judging dictionary entry similarity as well as a Monte Carlo method that evaluates the noise propagation properties of an MRF sequence. Consistency of these measures for different sequence lengths as well as the capability to predict actual sequence performance in both phantom and in vivo measurements was analyzed. While the dot-product-based measures yielded inconsistent results for different sequence lengths, the Monte Carlo method was in a good agreement with phantom experiments. In particular, the Monte Carlo method could accurately predict the performance of different flip angle patterns in actual measurements. The proposed Monte Carlo method provides an appropriate measure of MRF sequence encoding capability and may be used for sequence optimization.

Magnetic Resonance Fingerprinting with short relaxation intervals

PURPOSE

The aim of this study was to investigate a technique for improving the performance of Magnetic Resonance Fingerprinting (MRF) in repetitive sampling schemes, in particular for 3D MRF acquisition, by shortening relaxation intervals between MRF pulse train repetitions.

MATERIAL AND METHODS

A calculation method for MRF dictionaries adapted to short relaxation intervals and non-relaxed initial spin states is presented, based on the concept of stationary fingerprints. The method is applicable to many different k-space sampling schemes in 2D and 3D. For accuracy analysis, T₁ and T₂ values of a phantom are determined by single-slice Cartesian MRF for different relaxation intervals and are compared with quantitative reference measurements. The relevance of slice profile effects is also investigated in this case. To further illustrate the capabilities of the method, an application to in-vivo spiral 3D MRF measurements is demonstrated.

RESULTS

The proposed computation method enables accurate parameter estimation even for the shortest relaxation intervals, as investigated for different sampling patterns in 2D and 3D. In 2D Cartesian measurements, we achieved a scan acceleration of more than a factor of two, while maintaining acceptable accuracy: The largest T₁ values of a sample set deviated from their reference values by 0.3% (longest relaxation interval) and 2.4% (shortest relaxation interval). The largest T₂ values showed systematic deviations of up to 10% for all relaxation intervals, which is discussed. The influence of slice profile effects for multislice acquisition is shown to become increasingly relevant for short relaxation intervals. In 3D spiral measurements, a scan time reduction of 36% was achieved, maintaining the quality of in-vivo T1 and T2 maps.

CONCLUSIONS

Reducing the relaxation interval between MRF sequence repetitions using stationary fingerprint dictionaries is a feasible method to improve the scan efficiency of MRF sequences. The method enables fast implementations of 3D spatially resolved MRF.

Highly efficient respiratory motion compensated free-breathing coronary MRA using golden-step Cartesian acquisition

PURPOSE

To develop an efficient 3D affine respiratory motion compensation framework for Cartesian whole-heart coronary magnetic resonance angiography (MRA).

MATERIALS AND METHODS

The proposed method achieves 100% scan efficiency by estimating the affine respiratory motion from the data itself and correcting the acquired data in the reconstruction process. For this, a golden-step Cartesian sampling with spiral profile ordering was performed to enable reconstruction of respiratory resolved images at any breathing position and with different respiratory window size. Affine motion parameters were estimated from image-based registration of 3D undersampled respiratory resolved images reconstructed with iterative SENSE and motion correction was performed directly in the reconstruction using a multiple-coils generalized matrix formulation method. This approach was tested on healthy volunteers and compared against a conventional diaphragmatic navigator-gated acquisition using quantitative and qualitative image quality assessment.

RESULTS

The proposed approach achieved 47 ± 12% and 59 ± 6% vessel sharpness for the right (RCA) and left (LAD) coronary arteries, respectively. Also, good quality visual scores of 2.4 ± 0.74 and 2.44 ± 0.86 were observed for the RCA and LAD (scores from 0, no to 4, excellent coronary vessel delineation). A not statically significant difference (P = 0.05) was found between the proposed method and an 8-mm navigator-gated and tracked scan, although scan efficiency increased from 61 ± 10% to 100%.

CONCLUSION

We demonstrate the feasibility of a new 3D affine respiratory motion correction technique for Cartesian whole-heart CMRA that achieves 100% scan efficiency and therefore a predictable acquisition time. This approach yields image quality comparable to that of an 8-mm navigator-gated acquisition with lower scan efficiency. Further evaluation of this technique in patients is now warranted to determine its clinical use.

Regional differences in diabetic patients' pharmacotherapy in Bulgaria

BACKGROUND

The regional analyses play an important role in understanding a state of diabetes mellitus management and to support informed policy options. They need to be explored in more details in order to ensure an equal patients' access to health care services of the same value and quality.

AIM

The aim of this study is to analyze regional differences in a cost of diabetes therapy for insulin users in Bulgaria.

MATERIALS AND METHODS

It is a combined prospective and retrospective observational study with duration of 6 months. Diabetic patients on insulin therapy were recruited by 35 endocrinologists. Information about the health care resources used was collected within 3-prospective and 3 retrospective months in 2010 and 2011. The regional cost of illness analysis was performed by calculating the average cost attributable to each individual patient despite the fact that some might not use a particular health care resource. Subgroup analysis was performed for hospitalized patients.

RESULTS

A detailed analysis revealed cost differences in the regions, especially with more vulnerable population like Burgas and Pleven regions. Another reason for the cost differences is the type of insulin or type of therapy. Our study confirms the fact that the hospitalizations are the major cost driver. Rising diabetes prevalence and associated costs, including hospitalizations and management of diabetes complications, are a growing concern. The last possible reason for regional differences is the patients' characteristics and therapy differences. We add evidence demonstrating that the patients on insulin and OAD consume more resources including hospitalizations and suffer from more complications of diabetes.

CONCLUSIONS

Reasons for regional differences might have different origin as there are various population characteristics, type of therapy, socio economic status and others.

Compressed sensing reconstruction for magnetic resonance parameter mapping

Compressed sensing (CS) holds considerable promise to accelerate the data acquisition in magnetic resonance imaging by exploiting signal sparsity. Prior knowledge about the signal can be exploited in some applications to choose an appropriate sparsifying transform. This work presents a CS reconstruction for magnetic resonance (MR) parameter mapping, which applies an overcomplete dictionary, learned from the data model to sparsify the signal. The approach is presented and evaluated in simulations and in in vivo T(1) and T(2) mapping experiments in the brain. Accurate T(1) and T(2) maps are obtained from highly reduced data. This model-based reconstruction could also be applied to other MR parameter mapping applications like diffusion and perfusion imaging.

Compressed sensing for chemical shift-based water-fat separation

Multi echo chemical shift-based water-fat separation methods allow for uniform fat suppression in the presence of main field inhomogeneities. However, these methods require additional scan time for chemical shift encoding. This work presents a method for water-fat separation from undersampled data (CS-WF), which combines compressed sensing and chemical shift-based water-fat separation. Undersampling was applied in the k-space and in the chemical shift encoding dimension to reduce the total scanning time. The method can reconstruct high quality water and fat images in 2D and 3D applications from undersampled data. As an extension, multipeak fat spectral models were incorporated into the CS-WF reconstruction to improve the water-fat separation quality. In 3D MRI, reduction factors of above three can be achieved, thus fully compensating the additional time needed in three-echo water-fat imaging. The method is demonstrated on knee and abdominal in vivo data.

Automatic coil selection for channel reduction in SENSE-based parallel imaging

OBJECTIVE

Coil arrays with large number of receive elements allow improved imaging performance and higher signal-to-noise-ratio. The MR systems supporting these arrays have to handle an increased amount of data and higher reconstruction burden. To overcome these problems, data reduction techniques need to be applied, realized either by linear combination of the original coil data prior to reconstruction or by discarding particular data from unimportant coil elements.

MATERIALS AND METHODS

This work focuses on the latter approach and presents an efficient algorithm for automatic coil selection applicable to SENSE imaging. A singular value decomposition (SVD)-based coil selection is proposed that performs a coil element ranking quantifying the contribution of each coil element to the image reconstruction allowing appropriate coil selection. This approach makes use of the coil sensitivity information and takes reduction factor and phase encoding direction into account.

RESULTS

Simulations, phantom and in vivo experiments were performed to validate the SVD-based coil selection algorithm. The proposed approach proved to be computationally efficient without remarkable image quality degradation.

CONCLUSION

The SVD-based approach offers the opportunity for fast automatic coil selection. This could simplify clinical workflow and may, furthermore, pave the way for various 2D real-time and interventional applications.