Fast uncertainty quantification of tracer distribution in the brain interstitial fluid with multilevel and quasi Monte Carlo

Efficient uncertainty quantification algorithms are key to understand the propagation of uncertainty -- from uncertain input parameters to uncertain output quantities -- in high resolution mathematical models of brain physiology. Advanced Monte Carlo methods such as quasi Monte Carlo (QMC) and multilevel Monte Carlo (MLMC) have the potential to dramatically improve upon standard Monte Carlo (MC) methods, but their applicability and performance in biomedical applications is underexplored. In this paper, we design and apply QMC and MLMC methods to quantify uncertainty in a convection-diffusion model of tracer transport within the brain. We show that QMC outperforms standard MC simulations when the number of random inputs is small. MLMC considerably outperforms both QMC and standard MC methods and should therefore be preferred for brain transport models.Comment: Multilevel Monte Carlo, quasi Monte Carlo, brain simulation, brain fluids, finite element method, biomedical computing, random fields, diffusion-convectio

Finite Element Modeling Driven by Health Care and Aerospace Applications

This thesis concerns the development, analysis, and computer implementation of mesh generation algorithms encountered in finite element modeling in health care and aerospace. The finite element method can reduce a continuous system to a discrete idealization that can be solved in the same manner as a discrete system, provided the continuum is discretized into a finite number of simple geometric shapes (e.g., triangles in two dimensions or tetrahedrons in three dimensions). In health care, namely anatomic modeling, a discretization of the biological object is essential to compute tissue deformation for physics-based simulations. This thesis proposes an efficient procedure to convert 3-dimensional imaging data into adaptive lattice-based discretizations of well-shaped tetrahedra or mixed elements (i.e., tetrahedra, pentahedra and hexahedra). This method operates directly on segmented images, thus skipping a surface reconstruction that is required by traditional Computer-Aided Design (CAD)-based meshing techniques and is convoluted, especially in complex anatomic geometries. Our approach utilizes proper mesh gradation and tissue-specific multi-resolution, without sacrificing the fidelity and while maintaining a smooth surface to reflect a certain degree of visual reality. Image-to-mesh conversion can facilitate accurate computational modeling for biomechanical registration of Magnetic Resonance Imaging (MRI) in image-guided neurosurgery. Neuronavigation with deformable registration of preoperative MRI to intraoperative MRI allows the surgeon to view the location of surgical tools relative to the preoperative anatomical (MRI) or functional data (DT-MRI, fMRI), thereby avoiding damage to eloquent areas during tumor resection. This thesis presents a deformable registration framework that utilizes multi-tissue mesh adaptation to map preoperative MRI to intraoperative MRI of patients who have undergone a brain tumor resection. Our enhancements with mesh adaptation improve the accuracy of the registration by more than 5 times compared to rigid and traditional physics-based non-rigid registration, and by more than 4 times compared to publicly available B-Spline interpolation methods. The adaptive framework is parallelized for shared memory multiprocessor architectures. Performance analysis shows that this method could be applied, on average, in less than two minutes, achieving desirable speed for use in a clinical setting. The last part of this thesis focuses on finite element modeling of CAD data. This is an integral part of the design and optimization of components and assemblies in industry. We propose a new parallel mesh generator for efficient tetrahedralization of piecewise linear complex domains in aerospace. CAD-based meshing algorithms typically improve the shape of the elements in a post-processing step due to high complexity and cost of the operations involved. On the contrary, our method optimizes the shape of the elements throughout the generation process to obtain a maximum quality and utilizes high performance computing to reduce the overheads and improve end-user productivity. The proposed mesh generation technique is a combination of Advancing Front type point placement, direct point insertion, and parallel multi-threaded connectivity optimization schemes. The mesh optimization is based on a speculative (optimistic) approach that has been proven to perform well on hardware-shared memory. The experimental evaluation indicates that the high quality and performance attributes of this method see substantial improvement over existing state-of-the-art unstructured grid technology currently incorporated in several commercial systems. The proposed mesh generator will be part of an Extreme-Scale Anisotropic Mesh Generation Environment to meet industries expectations and NASA\u27s CFD visio

Impact of within-voxel heterogeneity in fibre geometry on spherical deconvolution

Axons in white matter have been shown to have varying geometries within a bundle using ex vivo imaging techniques, but what does this mean for diffusion MRI (dMRI) based spherical deconvolution (SD)? SD attempts to estimate the fibre orientation distribution function (fODF) by assuming a single dMRI fibre response function (FRF) for all white matter populations and deconvolving this FRF from the dMRI signal at each voxel to estimate the fODF. Variable fibre geometry within a bundle however suggests the FRF might not be constant even within a single voxel. We test what impact realistic fibre geometry has on SD by simulating the dMRI signal in a range of realistic white matter numerical phantoms, including synthetic phantoms and real axons segmented from electron microscopy. We demonstrate that variable fibre geometry leads to a variable FRF across axons and that in general no single FRF is effective to recover the underlying fibre orientation distribution function (fODF). This finding suggests that assuming a single FRF can lead to misestimation of the fODF, causing further downstream errors in techniques such as tractography

Doctor of Philosophy

dissertationImage-based biomechanics, particularly numerical modeling using subject-specific data obtained via imaging, has proven useful for elucidating several biomechanical processes, such as prediction of deformation due to external loads, applicable to both normal function and pathophysiology of various organs. As the field evolves towards applications that stretch the limits of imaging hardware and acquisition time, the information traditionally expected as input for numerical routines often becomes incomplete or ambiguous, and requires specific acquisition and processing strategies to ensure physical accuracy and compatibility with predictive mathematical modeling. These strategies, often derivatives or specializations of traditional mechanics, effectively extend the nominal capability of medical imaging hardware providing subject-specific information coupled with the option of using the results for predictive numerical simulations. This research deals with the development of tools for extracting mechanical measurements from a finite set of imaging data and finite element analysis in the context of constructing structural atlases of the heart, understanding the biomechanics of the venous vasculature, and right ventricular failure. The tools include: (1) application of Hyperelastic Warping image registration to displacement-encoded MRI for reconstructing absolute displacement fields, (2) combination of imaging and a material parameter identification approach to measure morphology, deformation, and mechanical properties of vascular tissue, and (3) extrapolation of diffusion tensor MRI acquired at a single time point for the prediction the structural changes across the cardiac cycle with mechanical simulations. Selected tools were then applied to evaluate structural changes in a reversible animal model for right ventricular failure due to pressure overload

Mapping complex cell morphology in the grey matter with double diffusion encoding MR: a simulation study

This paper investigates the impact of cell body (soma) size and branching of cellular projections on diffusion MR imaging (dMRI) and spectroscopy (dMRS) signals for both standard single diffusion encoding (SDE) and more advanced double diffusion encoding (DDE) measurements using numerical simulations. The aim is to study the ability of dMRI/dMRS to characterize the complex morphology of brain grey matter, focusing on these two distinctive features. To this end, we employ a recently developed framework to create realistic meshes for Monte Carlo simulations, covering a wide range of soma sizes and branching orders of cellular projections, for diffusivities reflecting both water and metabolites. For SDE sequences, we assess the impact of soma size and branching order on the signal b-value dependence as well as the time dependence of the apparent diffusion coefficient (ADC). For DDE sequences, we assess their impact on the mixing time dependence of the signal angular modulation and of the estimated microscopic anisotropy, a promising contrast derived from DDE measurements. The SDE results show that soma size has a measurable impact on both the b-value and diffusion time dependence, for both water and metabolites. On the other hand, branching order has little impact on either, especially for water. In contrast, the DDE results show that soma size has a measurable impact on the signal angular modulation at short mixing times and the branching order significantly impacts the mixing time dependence of the signal angular modulation as well as of the derived microscopic anisotropy, for both water and metabolites. Our results confirm that soma size can be estimated from SDE based techniques, and most importantly, show for the first time that DDE measurements show sensitivity to the branching of cellular projections, paving the way for non-invasive characterization of grey matter morphology

Mapping complex cell morphology in the grey matter with double diffusion encoding MR: A simulation study

This paper investigates the impact of cell body (namely soma) size and branching of cellular projections on diffusion MR imaging (dMRI) and spectroscopy (dMRS) signals for both standard single diffusion encoding (SDE) and more advanced double diffusion encoding (DDE) measurements using numerical simulations. The aim is to investigate the ability of dMRI/dMRS to characterize the complex morphology of brain cells focusing on these two distinctive features of brain grey matter. To this end, we employ a recently developed computational framework to create three dimensional meshes of neuron-like structures for Monte Carlo simulations, using diffusion coefficients typical of water and brain metabolites. Modelling the cellular structure as realistically connected spherical soma and cylindrical cellular projections, we cover a wide range of combinations of sphere radii and branching order of cellular projections, characteristic of various grey matter cells. We assess the impact of spherical soma size and branching order on the b-value dependence of the SDE signal as well as the time dependence of the mean diffusivity (MD) and mean kurtosis (MK). Moreover, we also assess the impact of spherical soma size and branching order on the angular modulation of DDE signal at different mixing times, together with the mixing time dependence of the apparent microscopic anisotropy (μA), a promising contrast derived from DDE measurements. The SDE results show that spherical soma size has a measurable impact on both the b-value dependence of the SDE signal and the MD and MK diffusion time dependence for both water and metabolites. On the other hand, we show that branching order has little impact on either, especially for water. In contrast, the DDE results show that spherical soma size has a measurable impact on the DDE signal's angular modulation at short mixing times and the branching order of cellular projections significantly impacts the mixing time dependence of the DDE signal's angular modulation as well as of the derived μA, for both water and metabolites. Our results confirm that SDE based techniques may be sensitive to spherical soma size, and most importantly, show for the first time that DDE measurements may be more sensitive to the dendritic tree complexity (as parametrized by the branching order of cellular projections), paving the way for new ways of characterizing grey matter morphology, non-invasively using dMRS and potentially dMRI

Geometric models of brain white matter for microstructure imaging with diffusion MRI

The research presented in this thesis models the diffusion-weighted MRI signal within brain white matter tissue. We are interested in deriving descriptive microstructure indices such as white matter axon diameter and density from the observed diffusion MRI signal. The motivation is to obtain non-invasive reliable biomarkers for early diagnosis and prognosis of brain development and disease. We use both analytic and numerical models to investigate which properties of the tissue and aspects of the diffusion process affect the diffusion signal we measure. First we develop a numerical method to approximate the tissue structure as closely as possible. We construct three-dimensional meshes, from a stack of confocal microscopy images using the marching cubes algorithm. The experiment demonstrates the technique using a biological phantom (asparagus). We devise an MRI protocol to acquire data from the sample. We use the mesh models as substrates in Monte-Carlo simulations to generate synthetic MRI measurements. To test the feasibility of the method we compare simulated measurements from the three-dimensional mesh with scanner measurements from the same sample and simulated measurements from an extruded mesh and much simpler parametric models. The results show that the three-dimensional mesh model matches the data better than the extruded mesh and the parametric models revealing the sensitivity of the diffusion signal to the microstructure. The second study constructs a taxonomy of analytic multi-compartment models of white matter by combining intra- and extra-axonal compartments from simple models. We devise an imaging protocol that allows diffusion sensitisation parallel and perpendicular to tissue fibres. We use the protocol to acquire data from two fixed rat brains, which allows us to fit, study and evaluate the models. We conclude that models which incorporate non-zero axon radius describe the measurements most accurately. The key observation is a departure of signals in the parallel direction from the two-compartment models, suggesting restriction, most likely from glial cells or binding of water molecules to the membranes. The addition of the third compartment can capture this departure and explain the data. The final study investigates the estimates using in vivo brain diffusion measurements. We adjust the imaging protocol to allow an in vivo MRI acquisition of a rat brain and compare and assess the taxonomy of models. We then select the models that best explain the in vivo data and compare the estimates with those from the ex vivo measurements to identify any discrepancies. The results support the addition of the third compartment model as per the ex vivo findings, however the ranking of the models favours the zero radius intra-axonal compartments

Multiscale Cohort Modeling of Atrial Electrophysiology : Risk Stratification for Atrial Fibrillation through Machine Learning on Electrocardiograms

Patienten mit Vorhofflimmern sind einem fünffach erhöhten Risiko für einen ischämischen Schlaganfall ausgesetzt. Eine frühzeitige Erkennung und Diagnose der Arrhythmie würde ein rechtzeitiges Eingreifen ermöglichen, um möglicherweise auftretende Begleiterkrankungen zu verhindern. Eine Vergrößerung des linken Vorhofs sowie fibrotisches Vorhofgewebe sind Risikomarker für Vorhofflimmern, da sie die notwendigen Voraussetzungen für die Aufrechterhaltung der chaotischen elektrischen Depolarisation im Vorhof erfüllen. Mithilfe von Techniken des maschinellen Lernens könnten Fibrose und eine Vergrößerung des linken Vorhofs basierend auf P Wellen des 12-Kanal Elektrokardiogramms im Sinusrhythmus automatisiert identifiziert werden. Dies könnte die Basis für eine nicht-invasive Risikostrat- ifizierung neu auftretender Vorhofflimmerepisoden bilden, um anfällige Patienten für ein präventives Screening auszuwählen. Zu diesem Zweck wurde untersucht, ob simulierte Vorhof-Elektrokardiogrammdaten, die dem klinischen Trainingssatz eines maschinellen Lernmodells hinzugefügt wurden, zu einer verbesserten Klassifizierung der oben genannten Krankheiten bei klinischen Daten beitra- gen könnten. Zwei virtuelle Kohorten, die durch anatomische und funktionelle Variabilität gekennzeichnet sind, wurden generiert und dienten als Grundlage für die Simulation großer P Wellen-Datensätze mit genau bestimmbaren Annotationen der zugrunde liegenden Patholo- gie. Auf diese Weise erfüllen die simulierten Daten die notwendigen Voraussetzungen für die Entwicklung eines Algorithmus für maschinelles Lernen, was sie von klinischen Daten unterscheidet, die normalerweise nicht in großer Zahl und in gleichmäßig verteilten Klassen vorliegen und deren Annotationen möglicherweise durch unzureichende Expertenannotierung beeinträchtigt sind. Für die Schätzung des Volumenanteils von linksatrialem fibrotischen Gewebe wurde ein merkmalsbasiertes neuronales Netz entwickelt. Im Vergleich zum Training des Modells mit nur klinischen Daten, führte das Training mit einem hybriden Datensatz zu einer Reduzierung des Fehlers von durchschnittlich 17,5 % fibrotischem Volumen auf 16,5 %, ausgewertet auf einem rein klinischen Testsatz. Ein Long Short-Term Memory Netzwerk, das für die Unterscheidung zwischen gesunden und P Wellen von vergrößerten linken Vorhöfen entwickelt wurde, lieferte eine Genauigkeit von 0,95 wenn es auf einem hybriden Datensatz trainiert wurde, von 0,91 wenn es nur auf klinischen Daten trainiert wurde, die alle mit 100 % Sicherheit annotiert wurden, und von 0,83 wenn es auf einem klinischen Datensatz trainiert wurde, der alle Signale unabhängig von der Sicherheit der Expertenannotation enthielt. In Anbetracht der Ergebnisse dieser Arbeit können Elektrokardiogrammdaten, die aus elektrophysiologischer Modellierung und Simulationen an virtuellen Patientenkohorten resul- tieren und relevante Variabilitätsaspekte abdecken, die mit realen Beobachtungen übereinstim- men, eine wertvolle Datenquelle zur Verbesserung der automatisierten Risikostratifizierung von Vorhofflimmern sein. Auf diese Weise kann den Nachteilen klinischer Datensätze für die Entwicklung von Modellen des maschinellen Lernens entgegengewirkt werden. Dies trägt letztendlich zu einer frühzeitigen Erkennung der Arrhythmie bei, was eine rechtzeitige Auswahl geeigneter Behandlungsstrategien ermöglicht und somit das Schlaganfallrisiko der betroffenen Patienten verringert