Comparison of Otsu and an adapted Chan–Vese method to determine thyroid active volume using Monte Carlo generated SPECT images

Abstract

Background

The Otsu method and the Chan–Vese model are two methods proven to perform well in determining volumes of different organs and specific tissue fractions. This study aimed to compare the performance of the two methods regarding segmentation of active thyroid gland volumes, reflecting different clinical settings by varying the parameters: gland size, gland activity concentration, background activity concentration and gland activity concentration heterogeneity.

Methods

A computed tomography was performed on three playdough thyroid phantoms with volumes 20, 35 and 50 ml. The image data were separated into playdough and water based on Hounsfield values. Sixty single photon emission computed tomography (SPECT) projections were simulated by Monte Carlo method with isotope Technetium-99 m ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{\text {99m}}$$\end{document} Tc). Linear combinations of SPECT images were made, generating 12 different combinations of volume and background: each with both homogeneous thyroid activity concentration and three hotspots of different relative activity concentrations (48 SPECT images in total). The relative background levels chosen were 5 %, 10 %, 15 % and 20 % of the phantom activity concentration and the hotspot activities were 100 % (homogeneous case) 150 %, 200 % and 250 %. Poisson noise, (coefficient of variation of 0.8 at a 20 % background level, scattering excluded), was added before reconstruction was done with the Monte Carlo-based SPECT reconstruction algorithm Sahlgrenska Academy reconstruction code (SARec). Two different segmentation algorithms were applied: Otsu’s threshold selection method and an adaptation of the Chan–Vese model for active contours without edges; the results were evaluated concerning relative volume, mean absolute error and standard deviation per thyroid volume, as well as dice similarity coefficient.

Results

Both methods segment the images well and deviate similarly from the true volumes. They seem to slightly overestimate small volumes and underestimate large ones. Different background levels affect the two methods similarly as well. However, the Chan–Vese model deviates less and paired t-testing showed significant difference between distributions of dice similarity coefficients ( p-value  \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$<0.01$$\end{document} ).

Conclusions

The investigations indicate that the Chan–Vese model performs better and is slightly more robust, while being more challenging to implement and use clinically. There is a trade-off between performance and user-friendliness.