Skip to main content Accessibility help

An Efficient Numerical Method for Mean Curvature-Based Image Registration Model

  • Jin Zhang (a1), Ke Chen (a2), Fang Chen (a3) and Bo Yu (a1)


Mean curvature-based image registration model firstly proposed by Chumchob-Chen-Brito (2011) offered a better regularizer technique for both smooth and nonsmooth deformation fields. However, it is extremely challenging to solve efficiently this model and the existing methods are slow or become efficient only with strong assumptions on the smoothing parameter β. In this paper, we take a different solution approach. Firstly, we discretize the joint energy functional, following an idea of relaxed fixed point is implemented and combine with Gauss-Newton scheme with Armijo's Linear Search for solving the discretized mean curvature model and further to combine with a multilevel method to achieve fast convergence. Numerical experiments not only confirm that our proposed method is efficient and stable, but also it can give more satisfying registration results according to image quality.


Corresponding author

*Corresponding author. Email addresses: (J. Zhang), (K. Chen), (F. Chen), (B. Yu)


Hide All
[1] Acar, R. and Vogel, C.R., Analysis of bounded variation penalty methods for ill-posed problems, Inverse Problems, 10(6), pp. 12171229, 1994.
[2] Bai, Z.-Z., Huang, Y.-M. and Ng, M.K., On preconditioned iterative methods for Burgers equations, SIAM Journal on Scientific Computing, 29(1), pp. 415439, 2007.
[3] Bai, Z.-Z., Huang, Y.-M. and Ng, M.K., On preconditioned iterative methods for certain time-dependent partial differential equations, SIAM Journal on Numerical Analysis, 47(2), pp. 10191037, 2009.
[4] Bai, Z.-Z. and Ng, M.K., On inexact preconditioners for nonsymmetric matrices, SIAM Journal on Scientific Computing, 26(5), pp. 17101724, 2005.
[5] Bai, Z.-Z., Ng, M.K. and Wang, Z.-Q., Constraint preconditioners for symmetric indefinite matrices, SIAM Journal on Matrix Analysis and Applications, 31(2), pp. 410433, 2009.
[6] Brown, L.G., A survey of image registration techniques, ACM Computing Surveys, 24(4), pp. 325376, 1992.
[7] Chan, T.F., Chen, K. and Carter, J.L., Iterative methods for solving the dual formulation arising fromimage restoration, Electronic Transactions on Numerical Analysis, 26, pp. 299311, 2007.
[8] Chan, T.F. and Shen, J.-H., Image Processing and Analysis-Variational, PDE, Wavelet, and Stochastic Methods, SIAM, Philadelphia, 2005.
[9] Chen, H.-M., Arora, M.K. and Varshney, P.K., Mutual information based image registration for remote sensing data, International Journal of Remote Sensing, 24(18), pp. 37013706, 2003.
[10] Chen, K., Matrix Preconditioning Techniques and Applications, Cambridge University Press, UK, 2005.
[11] Chen, K. and Tai, X.-C., A nonlinear multigrid method for total variation minimization from image restoration, Journal of Scientific Computing, 32(2), pp. 115138, 2007.
[12] Chumchob, N. and Chen, K., A robust multigrid approach for variational image registration models, Journal of Computational and Applied Mathematics, 236(5), pp. 653674, 2011.
[13] Chumchob, N., Chen, K. and Brito-Loeza, C., A fourth order variational image registration model and its fast multigrid algorithm, Multiscale Modeling and Simulation, 9(1), pp. 89128, 2011.
[14] Fischer, B. and Modersitzki, J., Curvature based image registration, Journal of Mathematical Imaging and Vision, 18(1), pp. 8185, 2003.
[15] Frohn-Schauf, C., Henn, S., Hömke, L. and Witsch, K., Total variation based image registration, in Proceedings of the International Conference on PDE-Based Image Processing and Related Inverse Problems Series: Mathematics and Visualization, Springer-Verlag, pp. 305323, 2006.
[16] Frohn-Schauf, C., Henn, S. and Witsch, K., Multigrid based total variation image registration, Computing and Visualization in Science, 11(2), pp. 101113, 2008.
[17] Gill, P.E., Murray, W. and Wright, M.H., Practical Optimization, Academic Press, London, 1981.
[18] Gratton, S., Lawless, A.S. and Nichols, N.K., Approximate Gauss-Newton methods for nonlinear least squares problems, SIAM Journal on Optimization, 18(1), pp. 106132, 2007.
[19] Haber, E., Heldmann, S. and Modersitzki, J., Adaptive mesh refinement for non parametric image registration, SIAM Journal on Scientific Computing, 30(6), pp. 30123027, 2008.
[20] Haber, E., Heldmann, S. and Modersitzki, J., A computational framework for image-based constrained registration, Linear Algebra and its Applications, 431(3–4), pp. 459470, 2009.
[21] Haber, E., Horesh, R. and Modersitzki, J., Numerical optimization for constrained image registration, Numerical Linear Algebra with Applications, 17(2-3), pp. 343359, 2010.
[22] Haber, E. and Modersitzki, J., Numerical methods for volume preserving image registration, Inverse Problems, 20, pp. 16211638, 2004.
[23] Haber, E. and Modersitzki, J., A multilevel method for image registration, SIAM Journal on Scientific Computing, 27(5), pp. 15941607, 2006.
[24] Hajnal, J.V., Hawkes, D.J. and Hill, D.L.G., Medical Image Registration, The Biomedical Engineering Series, CRC Press, 2001.
[25] Henn, S., A multigrid method for a fourth-order diffusion equation with application to image processing, SIAM Journal on Scientific Computing, 27(3), pp. 831849, 2005.
[26] Köstler, H., Ruhnau, K. and Wienands, R., Multigrid solution of the optical flow system using a combined diffusion- and curvature-based regularizer, Numerical Linear Algebra with Applications, 15(2-3), pp. 201218, 2008.
[27] Maintz, J.B.A. and Viergever, M.A., A survey of medical image registration, Medical Image Analysis, 2(1), pp. 136, 1998.
[28] Modersitzki, J., Numerical Methods for Image Registration, Oxford University Press, New York, 2004.
[29] Modersitzki, J., FAIR: Flexible Algorithms for Image Registration, SIAM, Philadelphia, 2009.
[30] Nocedal, J. and Wright, S.J., Numerical Optimization, Springer-Verlag, New York, 1999.
[31] Rudin, L., Osher, S. and Fatemi, E., Nonlinear total variation based noise removal algorithms, Physica D: Nonlinear Phenomena, 60(1-4), pp. 259268, 1992.
[32] Savage, J. and Chen, K., An improved and accelerated nonlinear multigrid method for total-variation denoising, International Journal of Computer Mathematics, 82(8), pp. 10011015, 2005.
[33] Sorzano, C., Thévenaz, P. and Unser, M., Elastic registration of biological, images using vector-spline regularization, IEEE Transactions On Biomedical Engineering, 52(4), pp. 652663, 2005.
[34] Stürmer, M., Köstler, H. and Rüde, U., A fast full mulitigrid solver for applications in image processing, Numerical Linear Algebra with Applications, 15(2-3), pp. 187200, 2008.
[35] Thévenaz, P. and Unser, M., Optimization of mutual information for multiresolution image registration, IEEE Transactions on Image Processing, 9(12), pp. 20832099, 2000.
[36] Vogel, C.R., Computational Methods for Inverse Problems, SIAM, Philadelphia, 2002.
[37] Vogel, C.R. and Oman, M.E., Iterative methods for total variation denoising, SIAM Journal on Scientific Computing, 17(1), pp. 227238, 1996.


MSC classification


Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed