Multigrid algorithm based on hybrid smoothers for variational and selective segmentation models

ABSTRACT

Automatic segmentation of an image to identify all meaningful parts is one of the most challenging as well as useful tasks in a number of application areas. This is widely studied. Selective segmentation, less studied, aims to use limited user-specified information to extract one or more interesting objects (instead of all objects). Constructing a fast solver remains a challenge for both classes of model. However, our primary concern is on selective segmentation. In this work, we develop an effective multigrid algorithm, based on a new non-standard smoother to deal with non-smooth coefficients, to solve the underlying partial differential equations of a class of variational segmentation models in the level-set formulation. For such models, non-smoothness (or jumps) is typical as segmentation is only possible if edges (jumps) are present. In comparison with previous multigrid methods which were shown to produce an acceptable mean smoothing rate for related models, the new algorithm can ensure a small and global smoothing rate that is a sufficient condition for convergence. Our rate analysis is by local Fourier analysis and, with it, we design the corresponding iterative solver, improving on an ineffective line smoother. Numerical tests show that the new algorithm outperforms multigrid methods based on competing smoothers.

Keywords:

• Partial differential equations
• multigrid
• fast solvers
• local fourier analysis
• image segmentation
• jump coefficients

2010 AMS SUBJECT CLASSIFICATIONS:

• 65N55
• 65F05
• 65T50
• 62H35
• 93C20

Disclosure statement

No potential conflict of interest was reported by the authors.

Additional information

Funding

The first author wishes to thank the UK EPSRC, the Smith Institute and the Liverpool Heart and Chest Hospital for supporting the work through an Industrial CASE award. The second author is grateful to UK EPSRC for grants EP/K036939/1 and EP/N014499/1.