Comparing Magnitude versus Complex Data Fitting in Liver R2* Relaxometry
Arthur Peter Wunderlich1,2, Stefan Andreas Schmidt1, Meinrad Beer1, Armin Michael Nagel1,2, and Holger Cario3

1Clinic for Diagnsotic and Interventional Radiology, Ulm University, Medical Center, Ulm, Germany, 2Section for Experimental Radiology, Ulm University, Medical Center, Ulm, Germany, 3Department of Pediatrics and Adolescent Medicine, Ulm University, Medical Center, Ulm, Germany

Synopsis

Relaxometry of patient data was performed comparing the use of magnitude versus complex data. 94 patients suspected for liver iron overload were scanned with mulit-contrast GRE-MRI at 1.5 T, involving multiple TE, TR and FA. Analysis was performed as conjoined fit incorporating effects of fat/water dephasing. One fit was based on magnitude images modeling noise as free fit parameter, the other on complex data. Magnitude fit yielded similar results, but showed superior convergence and lower result uncertainty compared to the approach involving complex data.

Purpose

There are several publications concerning liver relaxometry based on fitting of magnitude MRI data. However, the full information is contained in complex data (1). This work was performed to compare both approaches directly in identical ROIs on real patient data.

Methods

94 regular transfused patients suspected for liver iron overload were investigated with multi-contrast GRE-MRI at 1.5 T (Siemens Avanto) using multiple TE/TR/FA, cf. Table 1, obtaining five transversal slices positioned over the liver. Two slices best suited for comprising vessel-free parts of liver tissue were analyzed by placing manually ROIs of fixed size containing 45 voxels, three in each slice. R2* relaxometry of averaged ROI signal values was performed, as conjoined fit of all FA and both echo-spacings, in two manners: a) using only magnitude images, modelling noise as free fit parameter with the expected signal as sum of squares of noise and ideal signal, and b) involving complex information (1). For both fitting processes, performed at identical ROIs, modulation of GRE signal caused by fat/water-dephasing was considered according to (1). Not only R2* values, but also their uncertainties were determined by the fit algorithm. Examinations were excluded from further analysis if fit failed to converge in more than three of six ROIs. Median R2* values for ROIs were determined for each patient seperately for both methods. Correlation between results was studied, and mean relative uncertainties for all patients and both methods were calculated.

Results

While magnitude fit worked in all cases, three examinations (3.1%) had to be excluded since the complex fit did not converge sufficiently. Correlation and correspondence of R2* values was good, cf. Fig. 1, yielding an R2 value of 0.985. Slope was calculated as 0.937, indicating that magnitude fit returned slightly larger values in some cases. The mean relative uncertainty of R2* values was 2.4 ± 1.8 % for magnitude and 12.8 ± 9.4 % for complex fit. Uncertainties as a function of resulting R2* values for both approaches are shown in Figs. 2 and 3.

Discussion

Multipeak fat-corrected R2* relaxometry worked on data acquired in a number of patients. ROI based analysis was employed to reduce uncertainty of measured data, assuming this would also diminishing result uncertainy. For three examinations, complex fit failed to converge, probably due to phase inconsistencies within ROIs, whereas magnitude fit worked well on all patients.

R2* values are comparable, however, slightly larger in some cases with the magnitude approach, yielding a tendency similar to that demonstrated in (1).

Uncertainty was lower for the magnitude fit by about a factor of five. This is surprising at the first sight since not magnitude, but complex fit considers complete information contained in MRI data. However, larger uncertainties for complex fit may again be caused by the ROI based approach. Inconsistent phase evolution, i.e. dephasing between different voxels in one ROI, leads to signal reduction for the averaged signal and may cause observed uncertainties.

Multiple breathholds with different FA were performed to get sufficient data for reliable fits in patients showing large R2*, which were present in our cohort of regular transfused patientes. These multiple acquisitions improve magnitude fit, but may impair complex fit because of scanner phase instabilities. To evaluate influence of phase incongruities within ROIs, voxel-wise fit should be studied. Though, in our preliminary experience, single-voxel fits frequently failed to converge especially for large R2* and showed large uncertainty of results. Both findings are probably due to larger noise in one voxel than in averaged ROI data. Furthermore, single-voxel fit is time-consuming and therefore less suited for routine application.

Conclusion

In our implementation of ROI based analysis, fit of magnitude data yielded R2* values comparable to fit involving complex data, at improved fit stability and reduced R2* uncertainty.

Acknowledgements

We acknowledge G. Glatting for help with the fit procedure.

References

1. Hernando D, Kramer JH, Reeder SB. Multipeak fat-corrected complex R2* relaxometry: theory, optimization, and clinical validation. Magnetic resonance in medicine 2013;70(5):1319-1331.

Figures

Table 1. Acquisition parameters for multi-contrast GRE-Sequence

Figure 1. Correlation of R2* values obtained with magnitude and complex fit. The displayed regression line is close to identity, with R2 of 0.99. However, slope of 0.94 indicates lack of congruence, with magnitude fit yielding slightly larger results in some cases.

Fig. 2. Uncertainties vs. R2* results for magnitude fit.

Fig. 3. Uncertainties vs. R2* results for complex fit.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
0360