Exploiting Information Geometry to Improve the Convergence Properties of Variational Active Contours

In this paper we seek to exploit information geometry in order to define the Riemannian structure of the statistical manifold associated with the Chan-Vese active contour model. This Riemannian structure is obtained through a relationship between the contour's Mumford-Shah energy functional and the likelihood of the categorical latent variables of a Gaussian mixture model. Accordingly the natural metric of the statistical manifold formed by the contours is determined by their Fisher information matrix. Mathematical developments show that this matrix has a closed-form analytic expression and is diagonal. Based on this, we subsequently develop a natural gradient algorithm for the Chan-Vese active contour, with application to image segmentation. Because the proposed method performs optimization on the parameters natural manifold it attains dramatically faster convergence rates than the Euclidean gradient descent algorithm commonly used in the literature. Experiments performed on standard test images from the active contour literature are presented and confirm that the proposed natural gradient algorithm delivers accurate segmentation results in few iterations. Comparisons with methods from the state of the art show that the proposed method converges extremely fast, and could improve significantly the speed of many existing image segmentation methods based on the Chan-Vese active contour as well as enable its application to new problems. A MATLAB implementation of the proposed method is available at http://www.stats.bris.ac.uk/~p12320/code/SmoothNaturalGradient4ChanVese.rar.

[1]  T. Steihaug The Conjugate Gradient Method and Trust Regions in Large Scale Optimization , 1983 .

[2]  Guillermo Sapiro,et al.  New Possibilities with Sobolev Active Contours , 2007, International Journal of Computer Vision.

[3]  Ke Chen,et al.  On Two Multigrid Algorithms for Modeling Variational Multiphase Image Segmentation , 2009, IEEE Transactions on Image Processing.

[4]  Mohamed-Jalal Fadili,et al.  Region-Based Active Contours with Exponential Family Observations , 2009, Journal of Mathematical Imaging and Vision.

[5]  Shun-ichi Amari,et al.  Differential-geometrical methods in statistics , 1985 .

[6]  Michael Hintermüller,et al.  An Inexact Newton-CG-Type Active Contour Approach for the Minimization of the Mumford-Shah Functional , 2004, Journal of Mathematical Imaging and Vision.

[7]  J. Sethian,et al.  FRONTS PROPAGATING WITH CURVATURE DEPENDENT SPEED: ALGORITHMS BASED ON HAMILTON-JACOB1 FORMULATIONS , 2003 .

[8]  Nicholas Ayache,et al.  Iconic feature based nonrigid registration: the PASHA algorithm , 2003, Comput. Vis. Image Underst..

[9]  Shun-ichi Amari,et al.  Why natural gradient? , 1998, Proceedings of the 1998 IEEE International Conference on Acoustics, Speech and Signal Processing, ICASSP '98 (Cat. No.98CH36181).

[10]  R. Deriche Recursively Implementing the Gaussian and its Derivatives , 1993 .

[11]  Bodo Rosenhahn,et al.  Analysis of Numerical Methods for Level Set Based Image Segmentation , 2009, ISVC.

[12]  Robert J. Renka,et al.  Image segmentation with a Sobolev gradient method , 2009 .

[13]  Jean-Philippe Pons,et al.  Generalized Gradients: Priors on Minimization Flows , 2007, International Journal of Computer Vision.

[14]  Guillermo Sapiro,et al.  Generalized Newton-Type Methods for Energy Formulations in Image Processing , 2009, SIAM J. Imaging Sci..

[15]  Xavier Bresson,et al.  Fast Global Minimization of the Active Contour/Snake Model , 2007, Journal of Mathematical Imaging and Vision.

[16]  Rick Chartrand,et al.  A Faster-converging Algorithm for Image Segmentation with a Modified Chan-Vese Model , 2007, IPCV.

[17]  Anthony J. Yezzi,et al.  Conformal metrics and true "gradient flows" for curves , 2005, Tenth IEEE International Conference on Computer Vision (ICCV'05) Volume 1.

[18]  Tony F. Chan,et al.  Active contours without edges , 2001, IEEE Trans. Image Process..

[19]  C. R. Rao,et al.  Information and the Accuracy Attainable in the Estimation of Statistical Parameters , 1992 .

[20]  Engin Mendi,et al.  Quasi-Newton Minimization for Active Contours with Chan-Vese Model , 2009, IPCV.