Hostname: page-component-76fb5796d-wq484 Total loading time: 0 Render date: 2024-04-30T05:57:08.863Z Has data issue: false hasContentIssue false
  • English
  • Français

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

Part of: Partial differential equations, initial value and time-dependent initial-boundary value problems Numerical linear algebra

Published online by Cambridge University Press:  31 January 2017

Jin Zhang ,
Ke Chen ,
Fang Chen  and
Bo Yu
Jin Zhang*
Affiliation:
School of Mathematical Sciences, Dalian University of Technology, Liaoning 116024, P R China
Ke Chen*
Affiliation:
Centre for Mathematical Imaging Techniques and Department of Mathematical Sciences, University of Liverpool, United Kingdom
Fang Chen*
Affiliation:
School of Applied Science, Beijing Information Science and Technology University, Beijing 100192, P R China
Bo Yu*
Affiliation:
School of Mathematical Sciences, Dalian University of Technology, Liaoning 116024, P R China
*
*Corresponding author. Email addresses:zhangjinsunny321@163.com (J. Zhang), k.chen@liv.ac.uk (K. Chen), chenfreesky@126.com (F. Chen), yubo@dlut.edu.cn (B. Yu)
*Corresponding author. Email addresses:zhangjinsunny321@163.com (J. Zhang), k.chen@liv.ac.uk (K. Chen), chenfreesky@126.com (F. Chen), yubo@dlut.edu.cn (B. Yu)
*Corresponding author. Email addresses:zhangjinsunny321@163.com (J. Zhang), k.chen@liv.ac.uk (K. Chen), chenfreesky@126.com (F. Chen), yubo@dlut.edu.cn (B. Yu)
*Corresponding author. Email addresses:zhangjinsunny321@163.com (J. Zhang), k.chen@liv.ac.uk (K. Chen), chenfreesky@126.com (F. Chen), yubo@dlut.edu.cn (B. Yu)
Get access

Abstract

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.

Keywords

Deformable image registrationregularizationmultilevelmean curvature

MSC classification

Secondary: 65F10: Iterative methods for linear systems 65M55: Multigrid methods; domain decomposition 68U10: Image processing
Type
Research Article
Information
East Asian Journal on Applied Mathematics , Volume 7 , Issue 1 , February 2017 , pp. 125 - 142
Copyright
Copyright © Global-Science Press 2017 

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1] Acar, R. and Vogel, C.R., Analysis of bounded variation penalty methods for ill-posed problems, Inverse Problems, 10(6), pp. 1217–1229, 1994.CrossRefGoogle Scholar
[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. 415–439, 2007.Google Scholar
[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. 1019–1037, 2009.Google Scholar
[4] Bai, Z.-Z. and Ng, M.K., On inexact preconditioners for nonsymmetric matrices, SIAM Journal on Scientific Computing, 26(5), pp. 1710–1724, 2005.CrossRefGoogle Scholar
[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. 410–433, 2009.Google Scholar
[6] Brown, L.G., A survey of image registration techniques, ACM Computing Surveys, 24(4), pp. 325–376, 1992.Google Scholar
[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. 299–311, 2007.Google Scholar
[8] Chan, T.F. and Shen, J.-H., Image Processing and Analysis-Variational, PDE, Wavelet, and Stochastic Methods, SIAM, Philadelphia, 2005.Google Scholar
[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. 3701–3706, 2003.Google Scholar
[10] Chen, K., Matrix Preconditioning Techniques and Applications, Cambridge University Press, UK, 2005.Google Scholar
[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. 115–138, 2007.CrossRefGoogle Scholar
[12] Chumchob, N. and Chen, K., A robust multigrid approach for variational image registration models, Journal of Computational and Applied Mathematics, 236(5), pp. 653–674, 2011.CrossRefGoogle Scholar
[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. 89–128, 2011.Google Scholar
[14] Fischer, B. and Modersitzki, J., Curvature based image registration, Journal of Mathematical Imaging and Vision, 18(1), pp. 81–85, 2003.Google Scholar
[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. 305–323, 2006.Google Scholar
[16] Frohn-Schauf, C., Henn, S. and Witsch, K., Multigrid based total variation image registration, Computing and Visualization in Science, 11(2), pp. 101–113, 2008.CrossRefGoogle Scholar
[17] Gill, P.E., Murray, W. and Wright, M.H., Practical Optimization, Academic Press, London, 1981.Google Scholar
[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. 106–132, 2007.Google Scholar
[19] Haber, E., Heldmann, S. and Modersitzki, J., Adaptive mesh refinement for non parametric image registration, SIAM Journal on Scientific Computing, 30(6), pp. 3012–3027, 2008.Google Scholar
[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. 459–470, 2009.Google Scholar
[21] Haber, E., Horesh, R. and Modersitzki, J., Numerical optimization for constrained image registration, Numerical Linear Algebra with Applications, 17(2-3), pp. 343–359, 2010.Google Scholar
[22] Haber, E. and Modersitzki, J., Numerical methods for volume preserving image registration, Inverse Problems, 20, pp. 1621–1638, 2004.Google Scholar
[23] Haber, E. and Modersitzki, J., A multilevel method for image registration, SIAM Journal on Scientific Computing, 27(5), pp. 1594–1607, 2006.Google Scholar
[24] Hajnal, J.V., Hawkes, D.J. and Hill, D.L.G., Medical Image Registration, The Biomedical Engineering Series, CRC Press, 2001.CrossRefGoogle Scholar
[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. 831–849, 2005.Google Scholar
[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. 201–218, 2008.Google Scholar
[27] Maintz, J.B.A. and Viergever, M.A., A survey of medical image registration, Medical Image Analysis, 2(1), pp. 1–36, 1998.Google Scholar
[28] Modersitzki, J., Numerical Methods for Image Registration, Oxford University Press, New York, 2004.Google Scholar
[29] Modersitzki, J., FAIR: Flexible Algorithms for Image Registration, SIAM, Philadelphia, 2009.Google Scholar
[30] Nocedal, J. and Wright, S.J., Numerical Optimization, Springer-Verlag, New York, 1999.Google Scholar
[31] Rudin, L., Osher, S. and Fatemi, E., Nonlinear total variation based noise removal algorithms, Physica D: Nonlinear Phenomena, 60(1-4), pp. 259–268, 1992.Google Scholar
[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. 1001–1015, 2005.Google Scholar
[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. 652–663, 2005.CrossRefGoogle ScholarPubMed
[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. 187–200, 2008.Google Scholar
[35] Thévenaz, P. and Unser, M., Optimization of mutual information for multiresolution image registration, IEEE Transactions on Image Processing, 9(12), pp. 2083–2099, 2000.Google Scholar
[36] Vogel, C.R., Computational Methods for Inverse Problems, SIAM, Philadelphia, 2002.Google Scholar
[37] Vogel, C.R. and Oman, M.E., Iterative methods for total variation denoising, SIAM Journal on Scientific Computing, 17(1), pp. 227–238, 1996.Google Scholar