Quarterly of Applied Mathematics

Quarterly of Applied Mathematics

Online ISSN 1552-4485; Print ISSN 0033-569X



Estimating diffeomorphic mappings between templates and noisy data: Variance bounds on the estimated canonical volume form

Authors: Daniel J. Tward, Partha P. Mitra and Michael I. Miller
Journal: Quart. Appl. Math. 77 (2019), 467-488
MSC (2010): Primary 92C55, 49K40
DOI: https://doi.org/10.1090/qam/1527
Published electronically: November 20, 2018
Full-text PDF

Abstract | References | Similar Articles | Additional Information

Abstract: Anatomy is undergoing a renaissance driven by the availability of large digital data sets generated by light microscopy. A central computational task is to map individual data volumes to standardized templates. This is accomplished by regularized estimation of a diffeomorphic transformation between the coordinate systems of the individual data and the template, building the transformation incrementally by integrating a smooth flow field. The canonical volume form of this transformation is used to quantify local growth, atrophy, or cell density. While multiple implementations exist for this estimation, less attention has been paid to the variance of the estimated diffeomorphism for noisy data. Notably, there is an infinite dimensional unobservable space defined by those diffeomorphisms which leave the template invariant. These form the stabilizer subgroup of the diffeomorphic group acting on the template. The corresponding flat directions in the energy landscape are expected to lead to increased estimation variance. Here we show that a least-action principle used to generate geodesics in the space of diffeomorphisms connecting the subject brain to the template removes the stabilizer. This provides reduced-variance estimates of the volume form. Using simulations we demonstrate that the asymmetric large deformation diffeomorphic mapping methods (LDDMM), which explicitly incorporate the asymmetry between idealized template images and noisy empirical images, provide lower variance estimators than their symmetrized counterparts (cf. ANTs). We derive Cramer-Rao bounds for the variances in the limit of small deformations. Analytical results are shown for the Jacobian in terms of perturbations of the vector fields and divergence of the vector field.

References [Enhancements On Off] (What's this?)

  • [1] Rick A. Adams, Klaas Enno Stephan, Harriet R. Brown, Christopher D. Frith, and Karl J. Friston, The computational anatomy of psychosis, Frontiers in Psychiatry 4 (2013), 47.
  • [2] Y. Amit, U. Grenander, and M. Piccioni, Structural image restoration through deformable templates, Journal of the American Statistical Association 86 (1991), no. 414, 376-381.
  • [3] J. Ashburner, Computational anatomy with the spm software, Magnetic Resonance Imaging 27 (2009), 1163-1174.
  • [4] J. Ashburner and K. J. Friston, Computational anatomy, Statistical Parametric Mapping The Analysis of Functional Brain Images (K. J. Friston, J. Ashburner, S. J. Kiebel, T. E. Nichols, and W. D. Penny, eds.), Academic Press, 2007, pp. 49-100.
  • [5] John Ashburner and Karl J. Friston, Voxel-based morphometry-the methods, Neuroimage 11 (2000), no. 6, 805-821.
  • [6] Brian Avants and James C. Gee, Geodesic estimation for large deformation anatomical shape averaging and interpolation, Neuroimage 23 (2004), S139-S150.
  • [7] Brian B. Avants, Charles L. Epstein, Murray Grossman, and James C. Gee, Symmetric diffeomorphic image registration with cross-correlation: evaluating automated labeling of elderly and neurodegenerative brain, Medical Image Analysis 12 (2008), no. 1, 26-41.
  • [8] Brian B. Avants, Murray Grossman, and James C. Gee, Symmetric diffeomorphic image registration: Evaluating automated labeling of elderly and neurodegenerative cortex and frontal lobe, Proceedings of the Third International Conference on Biomedical Image Registration (Berlin, Heidelberg), WBIR'06, Springer-Verlag, 2006, pp. 50-57.
  • [9] R. Bajcsy and S. Kovačič, Multiresolution elastic matching, Computer Vision, Graphics, and Image Processing 46 (1989), 1-21.
  • [10] R. Bajcsy, R. Lieberson, and M. Reivich, A computerized system for the elastic matching of deformed radiographic images to idealized atlas images., J Comput Assist Tomogr (1983), no. 7, 618-625.
  • [11] M. Faisal Beg, Michael I. Miller, Alain Trouvé, and Laurent Younes, Computing Large Deformation Metric Mappings via Geodesic Flows of Diffeomorphisms, International Journal of Computer Vision 61 (2005), no. 2, 139-157.
  • [12] Mirza Faisal Beg, Patrick A. Helm, Elliot McVeigh, Michael I. Miller, and Raimond L. Winslow, Computational cardiac anatomy using mri, Magnetic Resonance in Medicine: An Official Journal of the International Society for Magnetic Resonance in Medicine 52 (2004), no. 5, 1167-1174.
  • [13] Mirza Faisal Beg and Ali Khan, Symmetric data attachment terms for large deformation image registration, IEEE Transactions on Medical Imaging 26 (2007), no. 9, 1179-1189.
  • [14] F. L. Bookstein, Principal warps: Thin-plate splines and the decomposition of deformations, IEEE Trans. Pattern Anal. Mach. Intell. 11 (1989), no. 6, 567-585.
  • [15] F. L. Bookstein, Biometrics, biomathematics and the morphometric synthesis, Bulletin of Mathematical Biology 58 (1996), no. 2, 313-365.
  • [16] Fred L. Bookstein, Thin-plate splines and the atlas problem for biomedical images, Biennial International Conference on Information Processing in Medical Imaging, Springer, 1991, pp. 326-342.
  • [17] Vincent Camion and Laurent Younes, Geodesic interpolating splines, Energy Minimization Methods in Computer Vision and Pattern Recognition, Springer Berlin/Heidelberg, 2001, pp. 513-527.
  • [18] Can Ceritoglu, Lei Wang, Lynn D. Selemon, John G. Csernansky, Michael I. Miller, and J. Tilak Ratnanather, Large deformation diffeomorphic metric mapping registration of reconstructed 3d histological section images and in vivo mr images, Frontiers in Human Neuroscience 4 (2010), 43.
  • [19] Gary Christensen, Michael I. Miller, and Richard D. Rabbit, Deformable templates using large deformation kinematics, IEEE Transactions of Medical Imaging 5 (1995), 1435-1447.
  • [20] Gary E. Christensen, Sarang C. Joshi, and Michael I. Miller, Volumetric transformation of brain anatomy, IEEE Transactions on Medical Imaging 16 (1997), no. 6, 864-877.
  • [21] G. E. Christensen and H. J. Johnson, Consistent estimation, IEEE Tran. on Med. Imaging 20 (2001), no. 7.
  • [22] Jia Du, Alvina Goh, and Anqi Qiu, Diffeomorphic metric mapping of high angular resolution diffusion imaging based on riemannian structure of orientation distribution functions, IEEE Transactions on Medical Imaging 31 (2012), no. 5, 1021-1033.
  • [23] Jia Du, A. Pasha Hosseinbor, Moo K. Chung, Barbara B. Bendlin, Gaurav Suryawanshi, Andrew L. Alexander, and Anqi Qiu, Diffeomorphic metric mapping and probabilistic atlas generation of hybrid diffusion imaging based on bfor signal basis, Medical image Analysis 18 (2014), no. 7, 1002-1014.
  • [24] Jia Du, L. Younes, and A. Qiu, Whole brain diffeomorphic metric mapping via integration of sulcan and gyral curves, cortical surfaces, and images, Neuroimage 56 (2011), no. 1, 162-173.
  • [25] Paul Dupuis, Ulf Grenander, and Michael I. Miller, Variational problems on flows of diffeomorphisms for image matching, Quart. Appl. Math. 56 (1998), no. 3, 587–600. MR 1632326, https://doi.org/10.1090/qam/1632326
  • [26] Stanley Durrleman, Xavier Pennec, Alain Trouvé, Paul Thompson, and Nicholas Ayache, Inferring brain variability from diffeomorphic deformations of currents: an integrative approach., Medical Image Analysis 12 (2008), no. 5, 626-637.
  • [27] James C. Gee, Martin Reivich, and Ruzena Bajcsy, Elastically deforming a three-dimensional atlas to match anatomical brain images, University of Pennsylvania Institute for Research in Cognitive Science Technical Report No. IRCS-93-37 (1993).
  • [28] Joan Glaunès, Anqi Qiu, Michael I. Miller, and Laurent Younes, Large deformation diffeomorphic metric curve matching, International Journal of Computer Vision 80 (2008), no. 3, 317-336.
  • [29] Joan Glaunes, Alain Trouvé, and Laurent Younes, Diffeomorphic matching of distributions: A new approach for unlabelled point-sets and sub-manifolds matching, Computer Vision and Pattern Recognition, 2004. CVPR 2004. Proceedings of the 2004 IEEE Computer Society Conference on, vol. 2, Ieee, 2004, pp. II-712.
  • [30] Ulf Grenander and Michael I. Miller, Pattern theory: from representation to inference, Oxford University Press, Oxford, 2007. MR 2285439
  • [31] Ulf Grenander and Michael I. Miller, Computational anatomy: an emerging discipline, Quart. Appl. Math. 56 (1998), no. 4, 617–694. Current and future challenges in the applications of mathematics (Providence, RI, 1997). MR 1668732, https://doi.org/10.1090/qam/1668732
  • [32] Hákon Gudbjartsson and Samuel Patz, The rician distribution of noisy mri data, Magnetic Resonance in Medicine 34 (1995), no. 6, 910-914.
  • [33] J. W. Haller, A. Banerjee, G. E. Christensen, M. Gado, S. Joshi, M. I. Miller, Y. Sheline, M. W. Vannier, and J. G. Csernansky, Three-dimensional hippocampal mr morphometry with high-dimensional transformation of neuroanatomic atlas, Radiology 202 (1997), no. 2, 504-10.
  • [34] Allan R. Jones, Caroline C. Overly, and Susan M. Sunkin, The allen brain atlas: 5 years and beyond, Nature Reviews Neuroscience 10 (2009), no. 11, 821.
  • [35] Sarang Joshi, Brad Davis, Matthieu Jomier, and Guido Gerig, Unbiased diffeomorphic atlas construction for computational anatomy, NeuroImage 23 (2004), S151-S160.
  • [36] Sarang C. Joshi and Michael I. Miller, Landmark matching via large deformation diffeomorphisms, IEEE Trans. Image Process. 9 (2000), no. 8, 1357–1370. MR 1808275, https://doi.org/10.1109/83.855431
  • [37] Bilge Karaçali and Christos Davatzikos, Estimating topology preserving and smooth displacement fields, IEEE Transactions on Medical Imaging 23 (2004), no. 7, 868-880.
  • [38] David G. Kendall, A survey of the statistical theory of shape, Statist. Sci. 4 (1989), no. 2, 87–120. With discussion and a reply by the author. MR 1007558
  • [39] Michael D. Ketcha, Tharindu De Silva, Runze Han, Ali Uneri, Joseph Goerres, Matthew W. Jacobson, Sebastian Vogt, Gerhard Kleinszig, and Jeffrey H. Siewerdsen, Effects of image quality on the fundamental limits of image registration accuracy, IEEE Transactions on Medical Imaging 36 (2017), no. 10, 1997-2009.
  • [40] Ali R Khan, Lei Wang, and Mirza Faisal Beg, Multistructure large deformation diffeomorphic brain registration, IEEE Transactions on Biomedical Engineering 60 (2013), no. 2, 544-553.
  • [41] Jan Kybic, Fast no ground truth image registration accuracy evaluation: Comparison of bootstrap and hessian approaches, Biomedical Imaging: From Nano to Macro, 2008. ISBI 2008. 5th IEEE International Symposium on, IEEE, 2008, pp. 792-795.
  • [42] Jan Kybic, Bootstrap resampling for image registration uncertainty estimation without ground truth, IEEE Trans. Image Process. 19 (2010), no. 1, 64–73. MR 2729958, https://doi.org/10.1109/TIP.2009.2030955
  • [43] M. I. Miller, Computational anatomy: shape, growth, and atrophy comparison via diffeomorphisms, Neuroimage 29 (2004), no. Suppl 1:S, 19-33.
  • [44] M. I. Miller and L. Younes, Group Actions, Homeomorphisms, and Matching: A General Framework, International Journal of Computer Vision 41 (2001), no. 4, 61-84.
  • [45] Michael I. Miller, Alain Trouvé, and Laurent Younes, On the metrics and euler-lagrange equations of computational anatomy, Annual Review of Biomedical Engineering 4 (2002), no. 1, 375-405.
  • [46] Michael I. Miller, Alain Trouvé, and Laurent Younes, Geodesic shooting for computational anatomy, J. Math. Imaging Vision 24 (2006), no. 2, 209–228. MR 2227097, https://doi.org/10.1007/s10851-005-3624-0
  • [47] Michael I. Miller, Alain Trouvé, and Laurent Younes, Diffeomorphometry and geodesic positioning systems for human anatomy, Technology 1:36 (2014).
  • [48] Michael I. Miller, Alain Trouvé, and Laurent Younes, Hamiltonian systems and optimal control in computational anatomy: 100 years since d'arcy thompson., Annual Review of Biomed Engineering (2015), no. 17, 447-509.
  • [49] Susumu Mori, Kenichi Oishi, Hangyi Jiang, Li Jiang, Xin Li, Kazi Akhter, Kegang Hua, Andreia V. Faria, Asif Mahmood, Roger Woods, et al., Stereotaxic white matter atlas based on diffusion tensor imaging in an icbm template, Neuroimage 40 (2008), no. 2, 570-582.
  • [50] X. Pennec, From Riemannian Geometry to Computational Anatomy, Elements (2011).
  • [51] Petter Risholm, Steve Pieper, Eigil Samset, and William M. Wells, Summarizing and visualizing uncertainty in non-rigid registration, International Conference on Medical Image Computing and Computer-Assisted Intervention, Springer, 2010, pp. 554-561.
  • [52] Laurent Risser, François-Xavier Vialard, Robin Wolz, Maria Murgasova, Darryl D. Holm, and Daniel Rueckert, Simultaneous Multi-scale Registration Using Large Deformation Diffeomorphic Metric Mapping., IEEE Transactions on Medical Imaging 30 (2011), no. 10, 1746-59.
  • [53] Laurent Risser, François-Xavier Vialard, Habib Y. Baluwala, and Julia A. Schnabel, Piecewise-diffeomorphic image registration: Application to the motion estimation between 3d ct lung images with sliding conditions, Medical image analysis 17 (2013), no. 2, 182-193.
  • [54] J. H. Siewerdsen, L. E. Antonuk, Y. El-Mohri, J. Yorkston, W. Huang, J. M. Boudry, and I. A. Cunningham, Empirical and theoretical investigation of the noise performance of indirect detection, active matrix flat-panel imagers (amfpis) for diagnostic radiology, Medical Physics 24 (1997), no. 1, 71-89.
  • [55] Amber L. Simpson, Burton Ma, Elvis C. S. Chen, Randy E. Ellis, and A. James Stewart, Using registration uncertainty visualization in a user study of a simple surgical task, International Conference on Medical Image Computing and Computer-Assisted Intervention, Springer, 2006, pp. 397-404.
  • [56] Ivor J. A. Simpson, Julia A. Schnabel, Adrian R. Groves, Jesper L. R. Andersson, and Mark W. Woolrich, Probabilistic inference of regularisation in non-rigid registration, NeuroImage 59 (2012), no. 3, 2438-2451.
  • [57] Ivor J. A. Simpson, Mark W. Woolrich, Adrian R. Groves, and Julia A. Schnabel, Longitudinal brain mri analysis with uncertain registration, International Conference on Medical Image Computing and Computer-Assisted Intervention, Springer, 2011, pp. 647-654.
  • [58] Donald L. Snyder, Random point processes, Wiley-Interscience [John Wiley & Sons], New York-London-Sydney, 1975. MR 0501325
  • [59] Stefan Sommer, Mads Nielsen, François Lauze, and Xavier Pennec, A multi-scale kernel bundle for LDDMM: towards sparse deformation description across space and scales., Information Processing in Medical Imaging Proceedings of the Conference 22 (2011), no. 17, 624-635.
  • [60] Aristeidis Sotiras, Christos Davatzikos, and Nikos Paragios, Deformable medical image registration: A survey, IEEE Transactions on Medical Imaging 32 (2013), no. 7, 1153-1190.
  • [61] X. Tang, S. Mori, and M. I. Miller, Segmentation via multi-atlas lddmm, MICCAI 2012 Workshop on Multi-Atlas Labeling (2012), 123-133.
  • [62] Paul M. Thompson and Arthur W. Toga, A framework for computational anatomy, Computing and Visualization in Science 5 (2002), no. 1, 13-34.
  • [63] A. Toga and P. M. Thompson, Maps of the brain, Anat Rec. 265 (2001), no. 2, 37-53.
  • [64] Alain Trouvé, An approach of pattern recognition through infinite dimensional group action, Rapport de recherche du LMENS (1995).
  • [65] Mauro Piccioni, Sergio Scarlatti, and Alain Trouvé, A variational problem arising from speech recognition, SIAM J. Appl. Math. 58 (1998), no. 3, 753–771. MR 1616603, https://doi.org/10.1137/S0036139995292537
  • [66] Alain Trouvé and Laurent Younes, Shape spaces, Handbook of mathematical methods in imaging. Vol. 1, 2, 3, Springer, New York, 2015, pp. 1759–1817. MR 3560100
  • [67] D. Tward, M. Miller, A. Trouve, and L. Younes, Parametric surface diffeomorphometry for low dimensional embeddings of dense segmentations and imagery, IEEE Transactions on Pattern Analysis and Machine Intelligence PP (2016), no. 99, 1-1.
  • [68] Daniel J. Tward and Jeffrey H. Siewerdsen, Cascaded systems analysis of the 3d noise transfer characteristics of flat-panel cone-beam ct, Medical Physics 35 (2008), no. 12, 5510-5529.
  • [69] D. J. Tward, C. Ceritoglu, A. Kolasny, G. M. Sturgeon, W. P. Segars, M. I. Miller, and J. T. Ratnanather, Patient specific dosimetry phantoms using multichannel lddmm of the whole body, International Journal of Biomedical Imaging 2011.
  • [70] D. J. Tward, J. Ma, M. I. Miller, and L. Younes, Robust diffeomorphic mapping via geodesically controlled active shapes, International Journal of Biomedical Imaging (2013), no. Article ID 205494, 19 pages.
  • [71] F. Vadakkumpadan, H. Arevalo, C. Ceritoglu, M. Miller, and N. Trayanova, Image-based estimation of ventricular fiber orientations for personalized modeling of cardiac electrophysiology, IEEE Transactions on Medical Imaging 31 (2012), no. 5, 1051 - 1060.
  • [72] Marc Vaillant and Joan Glaunès, Surface matching via currents, Proceedings of Information Processing in Medical Imaging (IPMI 2005), number 3565 in Lecture Notes in Computer Science, 2005, pp. 381-392.
  • [73] Tom Vercauteren, Xavier Pennec, Aymeric Perchant, and Nicholas Ayache, Symmetric log-domain diffeomorphic registration: A demons-based approach, International Conference on Medical Image Computing and Computer-Assisted Intervention, Springer, 2008, pp. 754-761.
  • [74] Tom Vercauteren, Xavier Pennec, Aymeric Perchant, and Nicholas Ayache, Diffeomorphic demons: Efficient non-parametric image registration, NeuroImage 45 1 Suppl (2009), S61-72.
  • [75] François-Xavier Vialard, Laurent Risser, Daniel Rueckert, and Colin J. Cotter, Diffeomorphic 3D image registration via geodesic shooting using an efficient adjoint calculation, Int. J. Comput. Vis. 97 (2012), no. 2, 229–241. MR 2897670, https://doi.org/10.1007/s11263-011-0481-8
  • [76] Laurent Younes, Shapes and diffeomorphisms, Applied Mathematical Sciences, vol. 171, Springer-Verlag, Berlin, 2010. MR 2656312
  • [77] Laurent Younes, Felipe Arrate, and Michael I. Miller, Evolutions equations in computational anatomy, NeuroImage 45 (2009), no. 1, S40-S50.

Similar Articles

Retrieve articles in Quarterly of Applied Mathematics with MSC (2010): 92C55, 49K40

Retrieve articles in all journals with MSC (2010): 92C55, 49K40

Additional Information

Daniel J. Tward
Affiliation: Center for Imaging Science, Johns Hopkins University, Baltimore, Maryland 21218
Email: dtward@cis.jhu.edu

Partha P. Mitra
Affiliation: Cold Spring Harbor Laboratory, Cold Spring Harbor, New York 11724
Email: mitra@cshl.edu

Michael I. Miller
Affiliation: Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland 21218
Email: mim@cis.jhu.edu

DOI: https://doi.org/10.1090/qam/1527
Keywords: Computational anatomy, morphometry, cell density, Hamiltonian dynamics
Received by editor(s): September 26, 2018
Received by editor(s) in revised form: October 8, 2018
Published electronically: November 20, 2018
Additional Notes: This work was supported by the National Institutes of Health [P41-EB015909, R01-EB020062, R01-NS102670, R01-MH105660, U19-AG033655, U19-MH114821, U01-MH114824]; National Science Foundation 16-569 NeuroNex contract 1707298; the Kavli Neuroscience Discovery Institute; the Crick-Clay Professorship, CSHL; and the H N Mahabala Chair, IIT Madras.
Dedicated: Twenty years after Computational Anatomy was hatched by Ulf Grenander and Michael Miller at the Division of Applied Mathematics at Brown University, we revisit the core geodesic equations and nuisance non-identifiable stabilizing subgroup of the deformable template model examining the variance in estimating the fundamental form and associated Cramer-Rao Bound in the small deformation limit. This picks up on a tradition started by Ulf’s thesis advisor, which we would like to think Ulf would have appreciated.
Article copyright: © Copyright 2018 Brown University