Equilibrium signal intensity mapping, an MRI method for fast mapping of longitudinal relaxation rates and for image enhancement

Abstract

Introduction: Inhomogeneity of magnetic fields, both B(0) and B(1), has been a major challenge in magnetic resonance imaging (MRI). Field inhomogeneity leads to image artifacts and unreliability of signal intensity (SI) measurements. This work proposes and shows the feasibility of generating equilibrium signal intensity (SI(Eq)) maps that can be utilized either to speed up relaxation-rate measurement or to enhance image quality and relaxation-rate-based weighting in various applications.

Methods: A 1.5-T MRI scanner was used. In canines (n=4), myocardial infarction was induced, and 48 h after the administration of 0.05 mmol kg(-1) Gd(ABE-DTTA), a contrast agent with slow tissue kinetics, in vivo R(1) mapping was carried out using an inversion recovery (IR)-prepared, fast gradient-echo sequence with varying inversion times (TIs). To test the SI(Eq) mapping method without the confounding effects of motion and blood flow, we carried out ex vivo R(1) mapping after the administration of 0.2 mmol kg(-1) Gd(DTPA) using an IR-prepared, fast spin-echo sequence in another group of dogs (n=2). R(1,full) maps and SI(Eq) maps were generated from the data from both sequences by three-parameter nonlinear curve fitting of the SI versus TI dependence. R(1,full) maps served as the reference standard. Raw IR images were then divided by the SI(Eq) maps, yielding corrected SI maps (COSIMs). Additionally, R(1) values were calculated from each single-TI image separately, using the SI(Eq) value and a one-parameter curve-fitting procedure (R(1,single)). Voxelwise correlation analysis was carried out for the COSIMs and the R(1,single) maps, both versus the standard R(1,full) maps. Deviations of R(1,single) from R(1,full) were statistically evaluated.

Results: In vivo, COSIM versus R(1,full) showed significantly (P<.05) better correlation [correlation coefficient (CC)=0.95] than SI versus R(1,full) with a TI=700-800 ms, which is 200-300 ms longer than the tau(null) (500 ms) of viable myocardium. With such TIs, SI versus R(1,full) yielded CCs of 0.86-0.88. R(1,single) versus R(1,full) yielded a peak CC of 0.96 at TI=700-900 ms. Mean deviations of R(1,single) from R(1,full) were below 5% for TIs between 500 and 1000 ms. Ex vivo, where tau(null) was 300 ms, the advantage of correction with SI(Eq) was not in the improvement of linear correlation but more in the reduction of scatter. Peak CCs for SI versus R(1,full) and COSIM versus R(1,full) at TI=500 ms were 0.96 for both. The ex vivo CC for R(1,single) versus R(1,full) at TI=500 ms was 0.98. Mean deviations of R(1,single) from R(1,full) were below 5% for TIs between 400 and 700 ms.

Conclusions: Once the corresponding SI(Eq) map is obtained from a control stack, R(1) can be obtained accurately, using only a single IR image and without the need for a stack of TI-varied images. This approach could be applied in various dynamic MRI studies where short measurement time, once the dynamics has started, is of essence. When using this method with IR-prepared T(1)-weighted images, it is essential that the single TI be chosen such that the longitudinal relaxation in all voxels of interest would have passed tau(null). SI(Eq) maps are also useful in eliminating confounders from MR images to allow obtaining SI values that reflect more faithfully the relaxation parameter (R(1)) sought.