Two Fast Alternating Direction Optimization Methods for
Multiphase Segmentation
Lu Tan
1
,WeiboWei
2
, and Zhenkuan Pan
2
1,2
College of Information Engineering, Qingdao University, Qingdao, Shandong, China
luna086@163.com, njustwwb@163.com
Abstract - In this paper, two new methods associated with
alternating direction optimization, fast alternating direction
method of multipliers(FastADMM) and fast alternating
minimization algorithm(FastAMA), are proposed for image
segmentation using the multiphase Chan-Vese model, which is
on the basis of piecewise constant optimal approximations.
For these methods, we incorporate the variable splitting
approach and a ‘reset’ condition in order to update the
Lagrange multiplier and make sure the value of energy
functional is always positive. The Osher and Stethian level set
method, binary level set functions, thresholding method and
projection formula are applied in the implementation. Finally,
numerical results with rapid convergence are obtained by our
methods, which are also compared with those of some other
fast variational methods to demonstrate better effectiveness of
our methods.
Keywords: Multiphase segmentation, fast alternating
direction method of multipliers, fast alternating minimization
algorithm, active contours, level sets.
1 Introduction
In image segmentation, major advances were made in
two-phase image segmentation[1-3] in the early days.
Mumford-Shah model proposed by Mumford D and Shah J [4]
is regarded as the most significant region-based model. It has
been extended to a great deal of applications. In 2001, Tony F.
Chan and Luminita A. Vese proposed the Chan-Vese model
[5] for active contours to detect objects in a given image. It is
one of the simplified variants of Mumford-Shah model.
Nevertheless, as the complexity of the images increases, 2-
phase image segmentation is not able to meet the actual needs.
Therefore, multiphase segmentation is applied to satisfy the
demands. On the basis of Potts model [6][7] from statistical
mechanics, Zhao et al. [8] started to study multiphase motion
segmentation by using the level set method and proposed a
model which can represent n different regions by n level set
functions. In order to reduce the number of level set functions,
Chan et al. continued their work and proposed multiphase
segmentation model [9] which is a generalization of Chan-
Vese model. Their scheme can naturally avoid "overlap" and
"leakage" problem. But there are still some problems about
solving the global optimization, accuracy, stability and speed.
Alternating direction method of multipliers(ADMM) was
first described by Glowinski and Marocco [10] and alternating
minimization algorithm(AMA) was presented by Tseng [11].
These techniques are commonly known as the Split Bregman
Method [12], and are known to be an efficient solver for
problems involving the total-variation norm [13]. These
methods can be accelerated using optimal first order methods,
of the type first proposed by Nesterov [14]. And accelerated
variants of ADMM and AMA can be called FastADMM and
FastAMA.
In our paper, we design methods (FastADMM and
FastAMA) that can achieve computational efficiency and own
even faster convergence to solve the functional of multiphase
Chan-Vese model based on binary level sets framework
[15,16]. The gradient descent method (GDM) [17],
Chambolle’s dual method (DM) [18], alternating direction
method of multipliers(ADMM) [19] i.e. the augmented
Lagrangian method (ALM) [20], and alternating minimization
algorithm(AMA) [11] are used to compare with our methods.
But results of these methods which are used in multiphase
segmentation model will be obtained in a slow convergence.
Tom Goldstein et al. [19] introduced FastADMM and applied
it to solve the TV model as an example and some other
strongly convex problems. A predictor-corrector type
acceleration step is used in this method.
The remaining of this paper is organized as follows. In
Section 2, the binary level set formulation of the functional of
multiphase Chan-Vese model used in our paper is reviewed
along with its four traditional solution methods. Our proposed
methods are discussed in Section 3 and its iterative discrete
formulas for implementation will be presented in detail. In
Section 4, some numerical experiments are given to illustrate
the effectiveness of our method by comparing with other
methods. Finally a conclusion is given in section 5.
2 The multiphase Chan-Vese model and
its four traditional methods
2.1 The binary level set based formulation
In order to separate an image domain
into n
subdomains with
1
n
i
i
and
i
j
ij
∃
.Veseand
Chan defined up to
2
m
n
phases and
m
level set functions.
This way makes sure that each pixel
! ∀
,x
y #
will belong
to one, and only one phase.
Int'l Conf. IP, Comp. Vision, and Pattern Reco
nition
IPCV'15
113