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 T
E/T
R/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. R
2* 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 R
2* 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 R
2* 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
R
2* values was good, cf. Fig. 1, yielding an R
2 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 R
2* values was 2.4 ± 1.8 % for magnitude
and 12.8 ± 9.4 % for complex fit. Uncertainties
as a function of resulting R
2* 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 R
2* values comparable to fit involving complex data, at improved fit stability and reduced R
2* 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.