Unified Heat Kernel Regression for Diffusion, Kernel Smoothing and Wavelets on Manifolds and Its Application to Mandible Growth Modeling in CT Images

We present a novel kernel regression framework for smoothing scalar surface data using the Laplace-Beltrami eigenfunctions. Starting with the heat kernel constructed from the eigenfunctions, we formulate a new bivariate kernel regression framework as a weighted eigenfunction expansion with the heat kernel as the weights. The new kernel regression is mathematically equivalent to isotropic heat diffusion, kernel smoothing and recently popular diffusion wavelets. Unlike many previous partial differential equation based approaches involving diffusion, our approach represents the solution of diffusion analytically, reducing numerical inaccuracy and slow convergence. The numerical implementation is validated on a unit sphere using spherical harmonics. As an illustration, we have applied the method in characterizing the localized growth pattern of mandible surfaces obtained in CT images from subjects between ages 0 and 20 years by regressing the length of displacement vectors with respect to the template surface.Comment: Accepted in Medical Image Analysi

Validation of Transcranial Electrical Stimulation (TES) Finite Element Modeling Against MREIT Current Density Imaging in Human Subjects

abstract: Transcranial electrical stimulation (tES) is a non-invasive brain stimulation therapy that has shown potential in improving motor, physiological and cognitive functions in healthy and diseased population. Typical tES procedures involve application of weak current (< 2 mA) to the brain via a pair of large electrodes placed on the scalp. While the therapeutic benefits of tES are promising, the efficacy of tES treatments is limited by the knowledge of how current travels in the brain. It has been assumed that the current density and electric fields are the largest, and thus have the most effect, in brain structures nearby the electrodes. Recent studies using finite element modeling (FEM) have suggested that current patterns in the brain are diffuse and not concentrated in any particular brain structure. Although current flow modeling is useful means of informing tES target optimization, few studies have validated tES FEM models against experimental measurements. MREIT-CDI can be used to recover magnetic flux density caused by current flow in a conducting object. This dissertation reports the first comparisons between experimental data from in-vivo human MREIT-CDI during tES and results from tES FEM using head models derived from the same subjects. First, tES FEM pipelines were verified by confirming FEM predictions agreed with analytic results at the mesh sizes used and that a sufficiently large head extent was modeled to approximate results on human subjects. Second, models were used to predict magnetic flux density, and predicted and MREIT-CDI results were compared to validate and refine modeling outcomes. Finally, models were used to investigate inter-subject variability and biological side effects reported by tES subjects. The study demonstrated good agreements in patterns between magnetic flux distributions from experimental and simulation data. However, the discrepancy in scales between simulation and experimental data suggested that tissue conductivities typically used in tES FEM might be incorrect, and thus performing in-vivo conductivity measurements in humans is desirable. Overall, in-vivo MREIT-CDI in human heads has been established as a validation tool for tES predictions and to study the underlying mechanisms of tES therapies.Dissertation/ThesisDoctoral Dissertation Biomedical Engineering 201

Unfolding the hippocampus: An intrinsic coordinate system for subfield segmentations and quantitative mapping

The hippocampus, like the neocortex, has a morphological structure that is complex and variable in its folding pattern, especially in the hippocampal head. The current study presents a computational method to unfold hippocampal grey matter, with a particular focus on the hippocampal head where complexity is highest due to medial curving of the structure and the variable presence of digitations. This unfolding was performed on segmentations from high-resolution, T2-weighted 7T MRI data from 12 healthy participants and one surgical patient with epilepsy whose resected hippocampal tissue was used for histological validation. We traced a critical image feature composed of the hippocampal sulcus and stratum radiatum lacunosum-moleculare, (SRLM) in these images, then employed user-guided semi-automated techniques to detect and subsequently unfold the surrounding hippocampal grey matter. This unfolding was performed by solving Laplace\u27s equation in three dimensions of interest (long-axis, proximal-distal, and laminar). The resulting ‘unfolded coordinate space’ provides an intuitive way of mapping the hippocampal subfields in 2D space (long-axis and proximal-distal), such that similar borders can be applied in the head, body, and tail of the hippocampus independently of variability in folding. This unfolded coordinate space was employed to map intracortical myelin and thickness in relation to subfield borders, which revealed intracortical myelin differences that closely follow the subfield borders used here. Examination of a histological resected tissue sample from a patient with epilepsy reveals that our unfolded coordinate system has biological validity, and that subfield segmentations applied in this space are able to capture features not seen in manual tracing protocols

An automated, geometry-based method for hippocampal shape and thickness analysis

The hippocampus is one of the most studied neuroanatomical structures due to its involvement in attention, learning, and memory as well as its atrophy in ageing, neurological, and psychiatric diseases. Hippocampal shape changes, however, are complex and cannot be fully characterized by a single summary metric such as hippocampal volume as determined from MR images. In this work, we propose an automated, geometry-based approach for the unfolding, point-wise correspondence, and local analysis of hippocampal shape features such as thickness and curvature. Starting from an automated segmentation of hippocampal subfields, we create a 3D tetrahedral mesh model as well as a 3D intrinsic coordinate system of the hippocampal body. From this coordinate system, we derive local curvature and thickness estimates as well as a 2D sheet for hippocampal unfolding. We evaluate the performance of our algorithm with a series of experiments to quantify neurodegenerative changes in Mild Cognitive Impairment and Alzheimer's disease dementia. We find that hippocampal thickness estimates detect known differences between clinical groups and can determine the location of these effects on the hippocampal sheet. Further, thickness estimates improve classification of clinical groups and cognitively unimpaired controls when added as an additional predictor. Comparable results are obtained with different datasets and segmentation algorithms. Taken together, we replicate canonical findings on hippocampal volume/shape changes in dementia, extend them by gaining insight into their spatial localization on the hippocampal sheet, and provide additional, complementary information beyond traditional measures. We provide a new set of sensitive processing and analysis tools for the analysis of hippocampal geometry that allows comparisons across studies without relying on image registration or requiring manual intervention

Proceedings of the Fourth International Workshop on Mathematical Foundations of Computational Anatomy - Geometrical and Statistical Methods for Biological Shape Variability Modeling (MFCA 2013), Nagoya, Japan

International audienceComputational anatomy is an emerging discipline at the interface of geometry, statistics and image analysis which aims at modeling and analyzing the biological shape of tissues and organs. The goal is to estimate representative organ anatomies across diseases, populations, species or ages, to model the organ development across time (growth or aging), to establish their variability, and to correlate this variability information with other functional, genetic or structural information. The Mathematical Foundations of Computational Anatomy (MFCA) workshop aims at fostering the interactions between the mathematical community around shapes and the MICCAI community in view of computational anatomy applications. It targets more particularly researchers investigating the combination of statistical and geometrical aspects in the modeling of the variability of biological shapes. The workshop is a forum for the exchange of the theoretical ideas and aims at being a source of inspiration for new methodological developments in computational anatomy. A special emphasis is put on theoretical developments, applications and results being welcomed as illustrations. Following the first edition of this workshop in 2006, second edition in New-York in 2008, the third edition in Toronto in 2011, the forth edition was held in Nagoya Japan on September 22 2013