Abstract
Finding a meaningful 1–1 correspondence between different data, such as images or surface data, has important applications in various fields. It involves the optimization of certain energy functionals over the space of all diffeomorphisms. This type of optimization problems (called the diffeomorphism optimization problems, DOPs) is especially challenging, since the bijectivity of the mapping has to be ensured. Recently, a method, called the Beltrami holomorphic flow (BHF), has been proposed to solve the DOP using quasi-conformal theories (Lui et al. in J Sci Comput 50(3):557–585, 2012). The optimization problem is formulated over the space of Beltrami coefficients (BCs), instead of the space of all diffeomorphisms. BHF iteratively finds a sequence of BCs associated with a sequence of diffeomorphisms, using the gradient descent method, to minimize the energy functional. The use of BCs effectively controls the smoothness and bijectivity of the mapping, and hence makes it easier to handle the constrained optimization problem. However, the algorithm is computationally expensive. In this paper, we propose an efficient splitting algorithm, based on the classical alternating direction method of multiplier (ADMM), to solve the DOP. The basic idea is to split the energy functional into two energy terms: one involves the BC whereas the other involves the quasi-conformal map. Alternating minimization scheme can then be applied to minimize the energy functional. The proposed method significantly speeds up the previous BHF approach. It also extends the previous BHF algorithm to Riemann surfaces of arbitrary topologies, such as multiply-connected shapes. Experiments have been carried out on synthetic together with real surface data, which demonstrate the efficiency and efficacy of the proposed algorithm to solve the DOP.
Similar content being viewed by others
References
Lui, L.M., Wong, T.W., Zeng, W., Gu, X.F., Thompson, P.M., Chan, T.F., Yau, S.T.: Optimization of surface registrations using Beltrami holomorphic flow. J. Sci. Comput. 50(3), 557–585 (2012)
Angenent, S., Haker, S., Tannenbaum, A., Kikinis, R.: On the Laplace–Beltrami operator and brain surface flattening. IEEE Trans. Med. Imaging 18(8), 700–711 (1999)
Durrleman, S., Pennec, X., Trouve, A., Thompson, P., Ayache, N.: Measuring brain variability via sulcal lines registration: a diffeomorphic approach. In: Medical Image Computing and Computer-Assisted Intervention (MICCAI 2007) Lecture Notes in Computer Science 4791, pp. 675–682 Springer, Berlin, Heidelberg (2007)
Durrleman, S., Pennec, X., Trouve, A., Thompson, P., Ayache, N.: Inferring brain variability from diffeomorphic deformations of currents: an integrative approach. Med. Image Anal. 12, 626–637 (2008)
Fischl, B., Sereno, M.I., Tootell, R.B., Dale, A.M.: High-resolution intersubject averaging and a coordinate system for the cortical surface. Hum. Brain Mapp. 8, 272–284 (1999)
Gardiner, F.: Quasiconformal Teichmüller Theory. American Mathematical Society, Providence (2000)
Glaunès, J., Vaillant, M., Miller, M.I.: Landmark matching via large deformation diffeomorphisms on the sphere. J. Math. Imaging Vis. 20, 179–200 (2004)
Gu, X., Wang, Y., Chan, T., Thompson, P., Yau, S.-T.: Genus zero surface conformal mapping and its application to brain surface mapping. IEEE Trans. Med. Imaging 23(8), 949–958 (2004)
Haker, S., Angenent, S., Tannenbaum, A., Kikinis, R., Sapiro, G., Halle, M.: Conformal surface parameterization for texture mapping. IEEE Trans. Vis. Comput. Graph. 6(8), 181189 (2000)
Hurdal, M., Stephenson, K.: Cortical cartography using the discrete conformal approach of circle packings. NeuroImage 23, S119S128 (2004)
Hurdal, M.K., Stephenson, K.: Discrete conformal methods for cortical brain flattening. Neuroimage 45, 86–98 (2009)
Joshi, S., Miller, M.: Landmark matching via large deformation diffeomorphisms. IEEE Trans. Image Process. 9(8), 13571370 (2000)
Ju, L., Stern, J., Rehm, K., Schaper, K., Hurdal, M., Rottenberg, D.: Cortical surface flattening using least square conformal mapping with minimal metric distortion. In: IEEE International Symposium on Biomedical, Imaging pp. 77–80 (2004)
Lehto, O., Virtanen, K.: Quasiconformal Conformal Mappings in the Plane. Springer, New York (1973)
Gardiner, F., Lakic, N.: Quasiconformal Teichmuller Theory. American Mathematical Society, Providence (2000)
Leow, A., Yu, C., Lee, S., Huang, S., Protas, H., Nicolson, R., Hayashi, K., Toga, A., Thompson, P.: Brain structural mapping using a novel hybrid implicit/explicit framework based on the level-set method. NeuroImage 24(3), 910–927 (2005)
Lepore, N., Leow, P.T.A.D. : Landmark matching on the sphere using distance functions. In: IEEE International Symposium on Biomedical Imaging (ISBI2006), April 6–9 2006.
Lord, N.A., Ho, J., Vemuri, B., Eisenschenk, S.: Simultaneous registration and parcellation of bilateral hippocampal surface pairs for local asymmetry quantification. IEEE Trans. Med. Imaging 26(4), 471–478 (2007)
Lui, L., Wang, Y., Chan, T., Thompson, P.: Brain anatomical feature detection by solving partial differential equations on general manifolds. Discrete Contin. Dyn. Syst. B 7(3), 605–618 (2007)
Lui, L., Wang, Y., Chan, T., Thompson, P.: Landmark constrained genus zero surface conformal mapping and its application to brain mapping research. Appl. Numer. Math. 5, 847–858 (2007)
Lui, L.M., Thiruvenkadam, S., Wang, Y., Chan, T., Thompson, P.: Optimized conformal parameterization of cortical surfaces using shape based matching of landmark curves. SIAM J. Imaging Sci. 3(1), 52–78 (2010)
Lui, L.M., Wong, T.W., Gu, X.F., Thompson, P.M., Chan, T.F., Yau, S.T.: Hippocampal Shape Registration using Beltrami Holomorphic flow. In: Medical Image Computing and Computer Assisted Intervention(MICCAI). Part II, LNCS, vol. 6362, pp. 323–330 (2010)
Zeng, W., Lui, L.M., Luo, F., Chan, T.F., Yau, S.T., Gu, X.F.: Computing quasiconformal maps using an auxiliary metric and discrete curvature flow. Numer. Math. 121(4), 671–703 (2012)
Lui, L.M., Lam, K.C., Wong, T.W., Gu, X.F.: Texture map and video compression using Beltrami representation. SIAM J. Imaging Sci. 6(4), 1880–1902 (2013)
Lvy, B., Petitjean, S., Ray, N., Maillot, J.: Least squares conformal maps for automatic texture atlas generation. ACM Trans. Graph. (TOG) 21(3), 362–371 (2002)
Zeng, W., Gu, X.: Registration for 3D surfaces with large deformations using quasi-conformal curvature flow. In: IEEE Conference on Computer Vision and Pattern Recognition (CVPR11), Jun 20–25, Colorado Springs. Colorado, USA (2011)
Lyttelton, O., Bouchera, M., Robbinsa, S., Evans, A.: An unbiased iterative group registration template for cortical surface analysis. NeuroImage 34, 1535–1544 (2007)
Glowinski, R., Le Tallee, P.: Augmented Lagrangian and Operator-splitting Methods in Nonlinear Mechanics. Studies in Applied and Numerical Mathematics. Society for Industrial and Applied Mathematics (1989)
Boyd, S., Parikh, N., Chu, E., Peleato, B., Eckstein, J.: Distributed optimization and statistical learning via the alternating direction method of multipliers. Found. Trends Mach. Learn. 3(1), 1–122 (2010)
Glowinski, R., Marrocco, A.: Sur lapproximation par elements Lnis dordre un, et la resolution par penalisation-dualite dune classe de problemes de dirichlet nonlineaires. Rev. Francaise dAut. Inf. Rech. Oper. R–2, 41–76 (1975)
Shi, Y., Thompson, P.M., Dinov, I., Osher, S., Toga, A.W.: Direct cortical mapping via solving partial differential equations on implicit surfaces. Med. Image Anal. 11(3), 207–223 (2007)
Tosun, D., Rettmann, M., Prince, J.: Mapping techniques for aligning sulci across multiple brains. Med. Image Anal. 8, 295309 (2004)
Wang, Y., Chiang, M.-C., Thompson, P.M.: Automated surface matching using mutual information applied to Riemann surface structures. Proc. Med. Image Comput. Comput. Assist. Interv. MICCAI 2005(2), 666–674 (2005)
Jacobson, A., Sorkine, O.: A cotangent Laplacian for images as surfaces. ACM Trans. Graph. 25(3), 646–653 (2012)
Weber, O., Myles, A., Zorin, D.: Computing extremal quasiconformal maps. Comput. Graph. Forum 31(5), 16791689 (2012)
Wong, T.W., Zhao, H.: Computation of quasiconformal surface maps using discrete Beltrami flow. UCLA CAM report 12-85 (2013)
Lui, L.M., Lam, K.C., Yau, S.T., Gu, X.F. : Teichmüller extremal mapping and its applications to landmark matching registration. SIAM J. Imaging Sci. (2013). arXiv:1211.2569 ( arXiv:1210.8025)
Ng, T.C., Gu, X.F., Lui, L.M. : Computing Teichmüller extremal map of multiply-connected domains using Beltrami holomorphic flow. J. Sci. Comput. (2013). doi:10.1007/s10915-013-9791-z
Wang, Y., Lui, L., Gu, X., Hayashi, K., Chan, T., Toga, A., Thompson, P., Yau, S.: Brain surface conformal parameterization using riemann surface structure. IEEE Trans. Med. Imaging 26(6), 853–865 (2007)
Wang, Y., Lui, L.M., Chan, T.F., Thompson, P.M.: Optimization of brain conformal mapping with landmarks. Proc. Med. Image Comput. Comput. Assist. Interv. MICCAI 2005, 675–683 (2005)
Zeng, W., Luo, F., Yau, S.-T., Gu, X. : Surface quasi-conformal mapping by solving Beltrami equations. In: IMA Conference on the Mathematics of Surfaces XIII (2009)
Gu, X.F., Yau, S.T.: Computational Conformal Geometry. International Press, Boston (2007)
Acknowledgments
Lok Ming Lui’s research was supported in part by HKRGC GRF grant (CUHK project ID: 404612) and CUHK FIS grant 1902036.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Lui, L.M., Ng, T.C. A Splitting Method for Diffeomorphism Optimization Problem Using Beltrami Coefficients. J Sci Comput 63, 573–611 (2015). https://doi.org/10.1007/s10915-014-9903-4
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10915-014-9903-4