In the clinical MRI practice, it is common to assess liver iron overload by T2* multi-echo gradient-echo images. However, there is no full consensus about the best image analysis approach for the T2* measurements. The currently used methods involve manual drawing of a region of interest (ROI) within MR images of the liver. Evaluation of a representative liver T2* value is done by fitting an appropriate model to the signal decay within the ROIs vs. the echo time. The resulting T2* value may depend on both ROI placement and choice of the signal decay model.
The aim of this study was to understand how the choice of the analysis methodology may affect the accuracy of T2* measurements. A software model of the iron overloaded liver was inferred from MR images acquired from 40 thalassemia major patients. Different image analysis methods were compared exploiting the developed software model. Moreover, a method for global semiautomatic T2* measurement involving the whole liver was developed. The global method included automatic segmentation of parenchyma by an adaptive fuzzyclustering algorithm able to compensate for signal inhomogeneities. Global liver T2* value was evaluated using a pixel-wise technique and an optimized signal decay model. The global approach was compared with the ROI-based approach used in the clinical practice.
For the ROI-based approach, the intra-observer and inter-observer coefficients of variation (CoVs) were 3.7% and 5.6%, respectively. For the global analysis, the CoVs for intra-observers and inter-observers reproducibility were 0.85% and 2.87%, respectively. The variability shown by the ROI-based approach was acceptable for use in the clinical practice; however, the developed global method increased the accuracy in T2* assessment and significantly reduced the operator dependence and sampling errors. This global approach could be useful in the clinical arena for patients with borderline liver iron overload and/or requiring follow-up studies.