Shape Analysis with Hyperbolic Wasserstein Distance

Topological Embedding of Human Brain Networks with Applications to Dynamics of Temporal Lobe Epilepsy

We introduce a novel, data-driven topological data analysis (TDA) approach for embedding brain networks into a lower-dimensional space in quantifying the dynamics of temporal lobe epilepsy (TLE) obtained from resting-state functional magnetic resonance imaging (rs-fMRI). This embedding facilitates the orthogonal projection of 0D and 1D topological features, allowing for the visualization and modeling of the dynamics of functional human brain networks in a resting state. We then quantify the topological disparities between networks to determine the coordinates for embedding. This framework enables us to conduct a coherent statistical inference within the embedded space. Our results indicate that brain network topology in TLE patients exhibits increased rigidity in 0D topology but more rapid flections compared to that of normal controls in 1D topology.

Topological state-space estimation of functional human brain networks

We introduce an innovative, data-driven topological data analysis (TDA) technique for estimating the state spaces of dynamically changing functional human brain networks at rest. Our method utilizes the Wasserstein distance to measure topological differences, enabling the clustering of brain networks into distinct topological states. This technique outperforms the commonly used k-means clustering in identifying brain network state spaces by effectively incorporating the temporal dynamics of the data without the need for explicit model specification. We further investigate the genetic underpinnings of these topological features using a twin study design, examining the heritability of such state changes. Our findings suggest that the topology of brain networks, particularly in their dynamic state changes, may hold significant hidden genetic information.

Unified topological inference for brain networks in temporal lobe epilepsy using the Wasserstein distance

Persistent homology offers a powerful tool for extracting hidden topological signals from brain networks. It captures the evolution of topological structures across multiple scales, known as filtrations, thereby revealing topological features that persist over these scales. These features are summarized in persistence diagrams, and their dissimilarity is quantified using the Wasserstein distance. However, the Wasserstein distance does not follow a known distribution, posing challenges for the application of existing parametric statistical models. To tackle this issue, we introduce a unified topological inference framework centered on the Wasserstein distance. Our approach has no explicit model and distributional assumptions. The inference is performed in a completely data driven fashion. We apply this method to resting-state functional magnetic resonance images (rs-fMRI) of temporal lobe epilepsy patients collected from two different sites: the University of Wisconsin-Madison and the Medical College of Wisconsin. Importantly, our topological method is robust to variations due to sex and image acquisition, obviating the need to account for these variables as nuisance covariates. We successfully localize the brain regions that contribute the most to topological differences. A MATLAB package used for all analyses in this study is available at https://github.com/laplcebeltrami/PH-STAT.

Unified Topological Inference for Brain Networks in Temporal Lobe Epilepsy Using the Wasserstein Distance

Persistent homology offers a powerful tool for extracting hidden topological signals from brain networks. It captures the evolution of topological structures across multiple scales, known as filtrations, thereby revealing topological features that persist over these scales. These features are summarized in persistence diagrams, and their dissimilarity is quantified using the Wasserstein distance. However, the Wasserstein distance does not follow a known distribution, posing challenges for the application of existing parametric statistical models. To tackle this issue, we introduce a unified topological inference framework centered on the Wasserstein distance. Our approach has no explicit model and distributional assumptions. The inference is performed in a completely data driven fashion. We apply this method to resting-state functional magnetic resonance images (rs-fMRI) of temporal lobe epilepsy patients collected from two different sites: the University of Wisconsin-Madison and the Medical College of Wisconsin. Importantly, our topological method is robust to variations due to sex and image acquisition, obviating the need to account for these variables as nuisance covariates. We successfully localize the brain regions that contribute the most to topological differences. A MATLAB package used for all analyses in this study is available at https://github.com/laplcebeltrami/PH-STAT.

Solid waste shape description and generation based on spherical harmonics and probability density function

Transport and separation processes of solid waste can only be modelled successfully with discrete element methods in case the shape of the particles can be described accurately. The existing techniques for morphological data acquisition, such as computed tomography, laser scanning technique, optical interferometer, stereo photography and structured light technique, are laborious and require a large amount of realistic solid waste samples. Therefore, there is a pressing need for an alternative method to describe the shape of solid waste particles and to generate multiple variations of particles with almost similar shapes. In this paper, a new method to describe solid waste particles is proposed that is frequency-based and uses spherical harmonics (SHs). Additionally, a new shape generation method is introduced that uses the shape description of a single particle to generate an array of related shapes based on a probability density function with a dimensionless control factor η. The newly proposed methods were successfully applied to describe the complex shapes of pieces of metal and plastic scrap. The shapes of these pieces of scrap can be described adequately with 15° of SH expansion and the overall divergence is within 0.1 mm. Five different values for η were tested, which generated shapes with the same distribution as the original particle. Rising levels of η cause the morphological variation of the generated particles to increase. These new methods improve the modelling of transportation and separation processes.

Computing Univariate Neurodegenerative Biomarkers with Volumetric Optimal Transportation: A Pilot Study

Changes in cognitive performance due to neurodegenerative diseases such as Alzheimer's disease (AD) are closely correlated to the brain structure alteration. A univariate and personalized neurodegenerative biomarker with strong statistical power based on magnetic resonance imaging (MRI) will benefit clinical diagnosis and prognosis of neurodegenerative diseases. However, few biomarkers of this type have been developed, especially those that are robust to image noise and applicable to clinical analyses. In this paper, we introduce a variational framework to compute optimal transportation (OT) on brain structural MRI volumes and develop a univariate neuroimaging index based on OT to quantify neurodegenerative alterations. Specifically, we compute the OT from each image to a template and measure the Wasserstein distance between them. The obtained Wasserstein distance, Wasserstein Index (WI) for short to specify the distance to a template, is concise, informative and robust to random noise. Comparing to the popular linear programming-based OT computation method, our framework makes use of Newton's method, which makes it possible to compute WI in large-scale datasets. Experimental results, on 314 subjects (140 Aβ + AD and 174 Aβ- normal controls) from the Alzheimer's Disease Neuroimaging Initiative (ADNI) baseline dataset, provide preliminary evidence that the proposed WI is correlated with a clinical cognitive measure (the Mini-Mental State Examination (MMSE) score), and it is able to identify group difference and achieve a good classification accuracy, outperforming two other popular univariate indices including hippocampal volume and entorhinal cortex thickness. The current pilot work suggests the application of WI as a potential univariate neurodegenerative biomarker.

Hyperbolic Wasserstein Distance for Shape Indexing

Shape space is an active research topic in computer vision and medical imaging fields. The distance defined in a shape space may provide a simple and refined index to represent a unique shape. This work studies the Wasserstein space and proposes a novel framework to compute the Wasserstein distance between general topological surfaces by integrating hyperbolic Ricci flow, hyperbolic harmonic map, and hyperbolic power Voronoi diagram algorithms. The resulting hyperbolic Wasserstein distance can intrinsically measure the similarity between general topological surfaces. Our proposed algorithms are theoretically rigorous and practically efficient. It has the potential to be a powerful tool for 3D shape indexing research. We tested our algorithm with human face classification and Alzheimer's disease (AD) progression tracking studies. Experimental results demonstrated that our work may provide a succinct and effective shape index.

Enhancing Diffusion MRI Measures By Integrating Grey and White Matter Morphometry With Hyperbolic Wasserstein Distance

In order to improve the preclinical diagnose of Alzheimer's disease (AD), there is a great deal of interest in analyzing the AD related brain structural changes with magnetic resonance image (MRI) analyses. As the major features, variation of the structural connectivity and the cortical surface morphometry provide different views of structural changes to determine whether AD is present on presymptomatic patients. However, the large scale tensor-valued information and relatively low imaging resolution in diffusion MRI (dMRI) have created huge challenges for analysis. In this paper, we propose a novel framework that improves dMRI analysis power by fusing cortical surface morphometry features from structural MRI (sMRI). We first compute the hyperbolic harmonic maps between cortical surfaces with the landmark constraints thus to precisely evaluate surface tensor-based morphometry. Meanwhile, the graph-based analysis of structural connectivity derived from dMRI is conducted. Next, we fuse these two features via the optimal mass transportation (OMT) and eventually the Wasserstein distance (WD) based single image index is computed as a potential clinical multimodality imaging score. We apply our framework to brain images of 20 AD patients and 20 matched healthy controls, randomly chosen from the Alzheimer's Disease Neuroimaging Initiative (AD-NI2) dataset. Our preliminary experimental results of group classification outperformed those of some other single dMRI-based features, such as regional hippocampal volume, mean scores of fractional anisotropy (FA) and mean axial (MD). The novel image fusion pipeline and simple imaging score of structural changes may benefit the preclinical AD and AD prevention research.

An Optimal Transportation based Univariate Neuroimaging Index

The alterations of brain structures and functions have been considered closely correlated to the change of cognitive performance due to neurodegenerative diseases such as Alzheimer's disease. In this paper, we introduce a variational framework to compute the optimal transformation (OT) in 3D space and propose a univariate neuroimaging index based on OT to measure such alterations. We compute the OT from each image to a template and measure the Wasserstein distance between them. By comparing the distances from all the images to the common template, we obtain a concise and informative index for each image. Our framework makes use of the Newton's method, which reduces the computational cost and enables itself to be applicable to large-scale datasets. The proposed work is a generic approach and thus may be applicable to various volumetric brain images, including structural magnetic resonance (sMR) and fluorodeoxyglucose positron emission tomography (FDG-PET) images. In the classification between Alzheimer's disease patients and healthy controls, our method achieves an accuracy of 82.30% on the Alzheimers Disease Neuroimaging Initiative (ADNI) baseline sMRI dataset and outperforms several other indices. On FDG-PET dataset, we boost the accuracy to 88.37% by leveraging pairwise Wasserstein distances. In a longitudinal study, we obtain a 5% significance with p-value = 1.13×10⁵ in a t-test on FDG-PET. The results demonstrate a great potential of the proposed index for neuroimage analysis and the precision medicine research.

Conformal invariants for multiply connected surfaces: Application to landmark curve-based brain morphometry analysis

Landmark curves were widely adopted in neuroimaging research for surface correspondence computation and quantified morphometry analysis. However, most of the landmark based morphometry studies only focused on landmark curve shape difference. Here we propose to compute a set of conformal invariant-based shape indices, which are associated with the landmark curve induced boundary lengths in the hyperbolic parameter domain. Such shape indices may be used to identify which surfaces are conformally equivalent and further quantitatively measure surface deformation. With the surface Ricci flow method, we can conformally map a multiply connected surface to the Poincaré disk. Our algorithm provides a stable method to compute the shape index values in the 2D (Poincaré Disk) parameter domain. The proposed shape indices are succinct, intrinsic and informative. Experimental results with synthetic data and 3D MRI data demonstrate that our method is invariant under isometric transformations and able to detect brain surface abnormalities. We also applied the new shape indices to analyze brain morphometry abnormalities associated with Alzheimer' s disease (AD). We studied the baseline MRI scans of a set of healthy control and AD patients from the Alzheimer' s Disease Neuroimaging Initiative (ADNI: 30 healthy control subjects vs. 30 AD patients). Although the lengths of the landmarks in Euclidean space, cortical surface area, and volume features did not differ between the two groups, our conformal invariant based shape indices revealed significant differences by Hotelling' s T² test. The novel conformal invariant shape indices may offer a new sensitive biomarker and enrich our brain imaging analysis toolset for studying diagnosis and prognosis of AD.