A robust method for assessment of iron overload in liver by magnetic resonance imaging

An automated liver segmentation in liver iron concentration map using fuzzy c-means clustering combined with anatomical landmark data

Background

To estimate median liver iron concentration (LIC) calculated from magnetic resonance imaging, excluded vessels of the liver parenchyma region were defined manually. Previous works proposed the automated method for excluding vessels from the liver region. However, only user-defined liver region remained a manual process. Therefore, this work aimed to develop an automated liver region segmentation technique to automate the whole process of median LIC calculation.

Methods

553 MR examinations from 471 thalassemia major patients were used in this study. LIC maps (in mg/g dry weight) were calculated and used as the input of segmentation procedures. Anatomical landmark data were detected and used to restrict ROI. After that, the liver region was segmented using fuzzy c-means clustering and reduced segmentation errors by morphological processes. According to the clinical application, erosion with a suitable size of the structuring element was applied to reduce the segmented liver region to avoid uncertainty around the edge of the liver. The segmentation results were evaluated by comparing with manual segmentation performed by a board-certified radiologist.

Results

The proposed method was able to produce a good grade output in approximately 81% of all data. Approximately 11% of all data required an easy modification step. The rest of the output, approximately 8%, was an unsuccessful grade and required manual intervention by a user. For the evaluation matrices, percent dice similarity coefficient (%DSC) was in the range 86–92, percent Jaccard index (%JC) was 78–86, and Hausdorff distance (H) was 14–28 mm, respectively. In this study, percent false positive (%FP) and percent false negative (%FN) were applied to evaluate under- and over-segmentation that other evaluation matrices could not handle. The average of operation times could be reduced from 10 s per case using traditional method, to 1.5 s per case using our proposed method.

Conclusion

The experimental results showed that the proposed method provided an effective automated liver segmentation technique, which can be applied clinically for automated median LIC calculation in thalassemia major patients.

Improved R2* liver iron concentration assessment using a novel fuzzy c-mean clustering scheme

BACKGROUND

In thalassemia patients, R2* liver iron concentration (LIC) measurement is a common clinical tool for assessing iron overload and for determining necessary chelator dose and evaluating its efficacy. Despite the importance of accurate LIC measurement, existing methods suffer from LIC variability, especially at the severe iron overload range due to inclusion of vessel parts in LIC calculation. In this study, we build upon previous Fuzzy C-Mean (FCM) clustering work to formulate a scheme with superior performance in segmenting vessel pixels from the parenchyma. Our method (MIX-FCM) combines our novel 2D-FCM with the existing 1D-FCM algorithm. This study further assessed possible optimal clustering parameters (OP scheme) and proposed a semi-automatic (SA) scheme for routine clinical application.

METHODS

Segmentation of liver parenchyma and vessels was performed on T2* images and their LIC maps in 196 studies from 147 thalassemia major patients. We used manual segmentation as the reference. 1D-FCM clustering was performed on the acquired image alone and 2D-FCM used both the acquired image and its LIC data. To execute the MIX-FCM method, the best outcome (OP-MIX-FCM) was selected from the aforementioned methods and was compared to the SA-MIX-FCM scheme. We used the percent value of the normalized interquartile range (nIQR) to its median to evaluate the variability of all methods.

RESULTS

2D-FCM clustering is more effective than 1D-FCM clustering at the severe overload range only, but inferior for other ranges (where 1D-FCM provides suitable results). This complementary performance between the two methods allows MIX-FCM to improve results for all ranges. OP-MIX-FCM clustering error was 2.1 ± 2.3%, compared with 10.3 ± 9.9% and 7.0 ± 11.9% from 1D- and 2D-FCM clustering, respectively. SA-MIX-FCM result was comparable to OP-MIX-FCM result, with both schemes showing ability to decrease overall nIQR by approximately 30%.

CONCLUSION

Our proposed 2D-FCM algorithm is not as superior to 1D-FCM as hypothesized. In contrast, our MIX-FCM method benefits from the best of both methods to obtain the highest segmentation accuracy at all ranges. Moreover, segmentation accuracy of the practical scheme (SA-MIX-FCM) is comparable to segmentation accuracy of the reference scheme (OP-MIX-FCM). Finally, we confirmed that segmentation is crucial to improving LIC assessments, especially at the severe iron overload range.

The influence of the analysis technique on estimating liver iron overload using magnetic resonance imaging T2* quantification

Iron toxicity is the major cause of tissue damage in patients with iron overload. Iron deposits mainly in the liver, where its concentration closely correlates with whole body iron overload. Different techniques have been proposed for estimating iron content, with liver biopsy being the gold standard despite its invasiveness and influence by sampling error. Recently, magnetic resonance imaging (MRI) has been established as an effective technique for evaluating iron overload by measuring T2(*) in the liver. However, various factors associated with the adopted analysis technique, mainly the exponential fitting model and signal averaging method, affect the resulting measurements. In this study, we evaluate the influences of these factors on T2(*) measurement in numerical phantom, calibrated phantoms, and nine patients with different degrees of iron overload. The results show different performances among the fitting models and signal averaging methods, which are affected by SNR, image quality and signal homogeneity inside the selected ROI for analysis.