zbMATH — the first resource for mathematics

A Mumford-Shah level-set approach for the inversion and segmentation of SPECT/CT data. (English) Zbl 1213.94015
Summary: This paper presents a level-set based approach for the simultaneous reconstruction and segmentation of the activity as well as the density distribution from tomography data gathered by an integrated SPECT/CT scanner. Activity and density distributions are modeled as piecewise constant functions. The segmenting contours and the corresponding function values of both the activity and the density distribution are found as minimizers of a Mumford-Shah like functional over the set of admissible contours and – for fixed contours – over the spaces of piecewise constant density and activity distributions which may be discontinuous across their corresponding contours. For the latter step a Newton method is used to solve the nonlinear optimality system. Shape sensitivity calculus is used to find a descent direction for the cost functional with respect to the geometrical variables which leads to an update formula for the contours in the level-set framework. A heuristic approach for the insertion of new components for the activity as well as the density function is used. The method is tested for synthetic data with different noise levels.

94A08 Image processing (compression, reconstruction, etc.) in information and communication theory
34K29 Inverse problems for functional-differential equations
44A12 Radon transform
49M05 Numerical methods based on necessary conditions
65K10 Numerical optimization and variational techniques
92C50 Medical applications (general)
Full Text: DOI