Macroscopic modeling and simulations of room evacuation

We analyze numerically two macroscopic models of crowd dynamics: the classical Hughes model and the second order model being an extension to pedestrian motion of the Payne-Whitham vehicular traffic model. The desired direction of motion is determined by solving an eikonal equation with density dependent running cost, which results in minimization of the travel time and avoidance of congested areas. We apply a mixed finite volume-finite element method to solve the problems and present error analysis for the eikonal solver, gradient computation and the second order model yielding a first order convergence. We show that Hughes' model is incapable of reproducing complex crowd dynamics such as stop-and-go waves and clogging at bottlenecks. Finally, using the second order model, we study numerically the evacuation of pedestrians from a room through a narrow exit.Comment: 22 page

Multi-valued geodesic tractography for diffusion weighted imaging

Diffusion-Weighted Imaging (DWI) is a Magnetic Resonance(MR) technique that measures water diffusion characteristics in tissue for a given direction. The diffusion profile in a specific location can be obtained by combining the DWI measurements of different directions. The diffusion profile gives information about the underlying fibrous structure, e.g., in human brain white matter, based on the assumption that water molecules are moving less freely perpendicularly to the fibrous structure. From the DW-MRI measurements often a positive definite second-order tensor is defined, the so-called diffusion tensor (DT). Neuroscientists have begun using diffusion tensor images (DTI) to study a host of various disorders and neurodegenerative diseases including Parkinson, Alzheimer and Huntington. The techniques for reconstructing the fiber tracts based on diffusion profiles are known as tractography or fiber tracking. There are several ways to extract fibers from the raw diffusion data. In this thesis, we explain and apply geodesic-based tractography techniques specifically, where the assumption is that fibers follow the most efficient diffusion propagation paths. A Riemannian manifold is defined using as metric the inverse of the diffusion tensor. A shortest path in this manifold is one with the strongest diffusion along this path. Therefore geodesics (i.e., shortest paths) on this manifold follow the most efficient diffusion paths. The geodesics are often computed from the stationary Hamilton-Jacobi equation (HJ). One characteristic of solving the HJ equation is that it gives only the single-valued viscosity solution corresponding to the minimizer of the length functional. It is also well known that the solution of the HJ equation can develop discontinuities in the gradient space, i.e., cusps. Cusps occur when the correct solution should become multi-valued. HJ methods are not able to handle this situation. To solve this, we developed a multi-valued solution algorithm for geodesic tractography in a metric space defined by given by diffusion tensor imaging data. The algorithm can capture all possible geodesics arriving at a single voxel instead of only computing the first arrival. Our algorithm gives the possibility of applying different cost functions in a fast post-processing. Moreover, the algorithm can be used for capturing possible multi-path connections between two points. In this thesis, we first focus on the mathematical and numerical model for analytic and synthetic fields in twodimensional domains. Later, we present the algorithm in three-dimensions with examples of synthetic and brain data. Despite the simplicity of the DTI model, the tractography techniques using DT are shown to be very promising to reveal the structure of brain white matter. However, DTI assumes that each voxel contains fibers with only one main orientation and it is known that brain white matter has multiple fiber orientations, which can be arbitrary many in arbitrary directions. Recently, High Angular Resolution Diffusion Imaging (HARDI) acquisition and its modeling techniques have been developed to overcome this limitation. As a next contribution we propose an extension of the multi-valued geodesic algorithm to HARDI data. First we introduce the mathematical model for more complex geometries using Finsler geometry. Next, we propose, justify and exploit the numerical methods for computing the multi-valued solution of these equations

A unified framework for document image restoration

Ph.DDOCTOR OF PHILOSOPH

Fast Algorithms for Surface Reconstruction from Point Cloud

We consider constructing a surface from a given set of point cloud data. We explore two fast algorithms to minimize the weighted minimum surface energy in [Zhao, Osher, Merriman and Kang, Comp.Vision and Image Under., 80(3):295-319, 2000]. An approach using Semi-Implicit Method (SIM) improves the computational efficiency through relaxation on the time-step constraint. An approach based on Augmented Lagrangian Method (ALM) reduces the run-time via an Alternating Direction Method of Multipliers-type algorithm, where each sub-problem is solved efficiently. We analyze the effects of the parameters on the level-set evolution and explore the connection between these two approaches. We present numerical examples to validate our algorithms in terms of their accuracy and efficiency