Source author record

Allen Tannenbaum

Allen Tannenbaum appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
17topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2022arXiv

Regularized Optimal Mass Transport with Nonlinear Diffusion

In this paper, we combine nonlinear diffusion with the regularized optimal mass transport (rOMT) model. As we will demonstrate, this new approach provides further insights into certain applications of fluid flow analysis in the brain. From the point of view of image processing, the anisotropic diffusion method, based on Perona-Malik, explicitly considers edge information. Applied to rOMT analysis of glymphatic transport based on DCE-MRI data, this new framework appears to capture a larger advection-dominant volume.

preprint2022arXiv

Wasserstein Image Local Analysis: Histogram of Orientations, Smoothing and Edge Detection

The Histogram of Oriented Gradient is a widely used image feature, which describes local image directionality based on numerical differentiation. Due to its ill-posed nature, small noise may lead to large errors. Conventional HOG may fail to produce meaningful directionality results in the presence of noise, which is common in medical radiographic imaging. We approach the directionality problem from a novel perspective by the use of the optimal transport map of a local image patch to a uni-color patch of its mean. We decompose the transport map into sub-work costs in different directions. We evaluated the ability of the optimal transport to quantify tumor heterogeneity from brain MRI images of patients with glioblastoma multiforme from the TCIA. By considering the entropy difference of the extracted local directionality within tumor regions, we found that patients with higher entropy in their images, had statistically significant worse overall survival (p $=0.008$), which indicates that tumors exhibiting flows in many directions may be more malignant, perhaps reflecting high tumor histologic grade, a reflection of histologic disorganization. We also explored the possibility of solving classical image processing problems such as smoothing and edge detection via optimal transport. By looking for a 2-color patch with minimum transport distance to a local patch, we derive a nonlinear shock filter, which preserves edges. Moreover, we found that the color difference of the computed 2-color patch indicates whether there is a large change in color, i.e., an edge in the given patch. In summary, we expand the usefulness of optimal transport as an image local analysis tool, to extract robust measures of imaging tumor heterogeneity for outcomes prediction as well as image pre-processing. Because of its robust nature, we find it offers several advantages over the classical approaches.

preprint2020arXiv

An Analytical Approach for Insulin-like Growth Factor Receptor 1 and Mammalian Target of Rapamycin Blockades in Ewing Sarcoma

We present preliminary results that quantify network robustness and fragility of Ewing sarcoma (ES), a rare pediatric bone cancer that often exhibits de novo or acquired drug resistance. By identifying novel proteins or pathways susceptible to drug targeting, this formalized approach promises to improve preclinical drug development and may lead to better treatment outcomes. Toward that end, our network modeling focused upon the IGF-1R-PI3K-Akt-mTOR pathway, which is of proven importance in ES. The clinical response and proteomic networks of drug-sensitive parental cell lines and their drug-resistant counterparts were assessed using two small molecule inhibitors for IGF-1R (OSI-906 and NVP-ADW-742) and an mTOR inhibitor (mTORi), MK8669, such that protein-to-protein expression networks could be generated for each group. For the first time, mathematical modeling proves that drug resistant ES samples exhibit higher degrees of overall network robustness (e.g., the ability of a system to withstand random perturbations to its network configuration) to that of their untreated or short-term (72-hour) treated samples. This was done by leveraging previous work, which suggests that Ricci curvature, a key geometric feature of a given network, is positively correlated to increased network robustness. More importantly, given that Ricci curvature is a local property of the system, it is capable of resolving pathway fragility. In this note, we offer some encouraging yet limited insights in terms of system-level robustness of ES and lay the foundation for scope of future work in which a complete study will be conducted.

preprint2020arXiv

Fisher-Rao Regularized Transport Analysis of the Glymphatic System and Waste Drainage

In this work, a unified representation of all the time-varying dynamics is accomplished with a Lagrangian framework for analyzing Fisher-Rao regularized dynamical optimal mass transport (OMT) derived flows. While formally equivalent to the Eulerian based Schrödinger bridge OMT regularization scheme, the Fisher-Rao approach allows a simple and interpretable methodology for studying the flows of interest in the present work. The advantage of the proposed Lagrangian technique is that the time-varying particle trajectories and attributes are displayed in a single visualization. This provides a natural capability to identify and distinguish flows under different conditions. The Lagrangian analysis applied to the glymphatic system (brain waste removal pathway associated with Alzheimer's Disease) successfully captures known flows and distinguishes between flow patterns under two different anesthetics, providing deeper insights into altered states of waste drainage.

preprint2020arXiv

Multimarginal Wasserstein Barycenter for Stain Normalization and Augmentation

Variations in hematoxylin and eosin (H&E) stained images (due to clinical lab protocols, scanners, etc) directly impact the quality and accuracy of clinical diagnosis, and hence it is important to control for these variations for a reliable diagnosis. In this work, we present a new approach based on the multimarginal Wasserstein barycenter to normalize and augment H&E stained images given one or more references. Specifically, we provide a mathematically robust way of naturally incorporating additional images as intermediate references to drive stain normalization and augmentation simultaneously. The presented approach showed superior results quantitatively and qualitatively as compared to state-of-the-art methods for stain normalization. We further validated our stain normalization and augmentations in the nuclei segmentation task on a publicly available dataset, achieving state-of-the-art results against competing approaches.

preprint2019arXiv

Kernel Wasserstein Distance

The Wasserstein distance is a powerful metric based on the theory of optimal transport. It gives a natural measure of the distance between two distributions with a wide range of applications. In contrast to a number of the common divergences on distributions such as Kullback-Leibler or Jensen-Shannon, it is (weakly) continuous, and thus ideal for analyzing corrupted data. To date, however, no kernel methods for dealing with nonlinear data have been proposed via the Wasserstein distance. In this work, we develop a novel method to compute the L2-Wasserstein distance in a kernel space implemented using the kernel trick. The latter is a general method in machine learning employed to handle data in a nonlinear manner. We evaluate the proposed approach in identifying computerized tomography (CT) slices with dental artifacts in head and neck cancer, performing unsupervised hierarchical clustering on the resulting Wasserstein distance matrix that is computed on imaging texture features extracted from each CT slice. Our experiments show that the kernel approach outperforms classical non-kernel approaches in identifying CT slices with artifacts.

preprint2016arXiv

Interactive Image Segmentation From A Feedback Control Perspective

Image segmentation is a fundamental problem in computational vision and medical imaging. Designing a generic, automated method that works for various objects and imaging modalities is a formidable task. Instead of proposing a new specific segmentation algorithm, we present a general design principle on how to integrate user interactions from the perspective of feedback control theory. Impulsive control and Lyapunov stability analysis are employed to design and analyze an interactive segmentation system. Then stabilization conditions are derived to guide algorithm design. Finally, the effectiveness and robustness of proposed method are demonstrated.

preprint2016arXiv

Matrix Optimal Mass Transport: A Quantum Mechanical Approach

In this paper, we describe a possible generalization of the Wasserstein 2-metric, originally defined on the space of scalar probability densities, to the space of Hermitian matrices with trace one, and to the space of matrix-valued probability densities. Our approach follows a computational fluid dynamical formulation of the Wasserstein-2 metric and utilizes certain results from the quantum mechanics of open systems, in particular the Lindblad equation. It allows determining the gradient flow for the quantum entropy relative to this matricial Wasserstein metric. This may have implications to some key issues in quantum information theory.

preprint2016arXiv

Regularization and Interpolation of Positive Matrices

We consider certain matricial analogues of optimal mass transport of positive definite matrices of equal trace. The framework is motivated by the need to devise a suitable geometry for interpolating positive definite matrices in ways that allow controlling the apparent tradeoff between "aligning up their eigenstructure" and "scaling the corresponding eigenvalues". Indeed, motivation for this work is provided by power spectral analysis of multivariate time series where, linear interpolation between matrix-valued power spectra generates push-pop artifacts. Push-pop of power distribuion is objectionable as it corresponds to unrealistic response of scatterers.

preprint2016arXiv

Robust transport over networks

We consider transport over a strongly connected, directed graph. The scheduling amounts to selecting transition probabilities for a discrete-time Markov evolution which is designed to be consistent with certain initial and final marginals. The random evolution is selected to be closest to a prior measure on paths in the relative entropy sense, i.e., a Schroedinger bridge between the two marginals. This is an atypical stochastic control problem where the control consists in suitably modifying the transition mechanism. The prior can incorporate cost of traversing edges or allocate equal probability to all paths of equal length connecting any two given nodes, i.e., a uniform measure on paths. This latter choice relies on the so-called Ruelle-Bowen random walk and gives rise to a scheduling that tends to utilize all paths as uniformly as the topology allows. Thus, when the Ruelle-Bowen law is taken as prior, the transportation plan tends to lessen congestion and ensure a level of robustness. We show that the Ruelle-Bowen law is itself a Schroedinger bridge albeit with a prior that is not a probability measure. The paradigm of Schroedinger bridges as a mechanism for scheduling transport on networks can be adapted to graphs that are not strongly connected as well as to weighted graphs. The latter leads to transportation plans that effect a compromise between robustness and transportation cost.

preprint2015arXiv

A Quantitative Analysis of Localized Robustness of MYCN in Neuroblastoma

The amplification of the gene MYCN (V-myc myelocytomatosis viral-valeted oncogene, neuroblastoma derived) has been a well-documented indicator for poor prognosis in neuroblastoma, a childhood cancer. Unfortunately, there has been limited success in understanding MYCN functionality in the landscape of neuroblastoma and more importantly, given that MYCN has been deemed undruggable, the need to potentially illuminate key opportunities that indirectly target MYCN is of great interest. To this end, this work employs an emerging quantitative technique from network science, namely network curvature, to quantify the biological robustness of MYCN and its surrounding neighborhood. In particular, when amplified in Stage IV cancer, MYCN exhibits higher curvature (more robust) than those samples with under expressed MYCN levels. When examining the surrounding neighborhood, the above argument still holds for network curvature, but is lost when only analyzing differential expression - a common technique amongst oncologists and computational/molecular biologists. This finding points to the problem (and possible solution) of drug targeting in the context of complexity and indirect cell signaling affects that have often been obfuscated through traditional techniques.

preprint2015arXiv

Graph Curvature and the Robustness of Cancer Networks

The importance of studying properties of networks is manifest in diverse fields ranging from biology, engineering, physics, chemistry, neuroscience, and medicine. The functionality of networks with regard to performance, throughput, reliability and robustness is strongly linked to the underlying geometric and topological properties of the network and this is the focus of this paper, especially as applied to certain biological networks. The fundamental mathematical abstraction of a network as a weighted graph brings to bear the tools of graph theory--a highly developed subject of mathematical research. But more importantly, recently proposed geometric notions of curvature on very general metric measure spaces allow us to utilize a whole new set of tools and ideas that help quantify functionality and robustness of graphs. In particular, robustness is closely connected to network entropy which, in turn, is very closely related to curvature. We will see that there are a number of alternative notions of discrete curvature that are compatible with the classical Riemannian definition, each having its own advantages and disadvantages, and are relevant to networks of interest. We will concentrate on the role of curvature for certain key cancer networks in order to quantitatively indicate their apparent functional robustness relative to their normal counterparts.

preprint2015arXiv

Market Fragility, Systemic Risk, and Ricci Curvature

Measuring systemic risk or fragility of financial systems is a ubiquitous task of fundamental importance in analyzing market efficiency, portfolio allocation, and containment of financial contagions. Recent attempts have shown that representing such systems as a weighted graph characterizing the complex web of interacting agents over some information flow (e.g., debt, stock returns, shareholder ownership) may provide certain keen insights. Here, we show that fragility, or the ability of system to be prone to failures in the face of random perturbations, is negatively correlated with geometric notion of Ricci curvature. The key ingredient relating fragility and curvature is entropy. As a proof of concept, we examine returns from a set of stocks comprising the S\&P 500 over a 15 year span to show that financial crashes are more robust compared to normal "business as usual" fragile market behavior - i.e., Ricci curvature is a "crash hallmark." Perhaps more importantly, this work lays the foundation of understanding of how to design systems and policy regulations in a manner that can combat financial instabilities exposed during the 2007-2008 crisis.

preprint2013arXiv

Linear models based on noisy data and the Frisch scheme

We address the problem of identifying linear relations among variables based on noisy measurements. This is, of course, a central question in problems involving "Big Data." Often a key assumption is that measurement errors in each variable are independent. This precise formulation has its roots in the work of Charles Spearman in 1904 and of Ragnar Frisch in the 1930's. Various topics such as errors-in-variables, factor analysis, and instrumental variables, all refer to alternative formulations of the problem of how to account for the anticipated way that noise enters in the data. In the present paper we begin by describing the basic theory and provide alternative modern proofs to some key results. We then go on to consider certain generalizations of the theory as well applying certain novel numerical techniques to the problem. A central role is played by the Frisch-Kalman dictum which aims at a noise contribution that allows a maximal set of simultaneous linear relations among the noise-free variables --a rank minimization problem. In the years since Frisch's original formulation, there have been several insights including trace minimization as a convenient heuristic to replace rank minimization. We discuss convex relaxations and certificates guaranteeing global optimality. A complementary point of view to the Frisch-Kalman dictum is introduced in which models lead to a min-max quadratic estimation error for the error-free variables. Points of contact between the two formalisms are discussed and various alternative regularization schemes are indicated.

preprint2013arXiv

Matrix-valued Monge-Kantorovich Optimal Mass Transport

We formulate an optimal transport problem for matrix-valued density functions. This is pertinent in the spectral analysis of multivariable time-series. The "mass" represents energy at various frequencies whereas, in addition to a usual transportation cost across frequencies, a cost of rotation is also taken into account. We show that it is natural to seek the transportation plan in the tensor product of the spaces for the two matrix-valued marginals. In contrast to the classical Monge-Kantorovich setting, the transportation plan is no longer supported on a thin zero-measure set.

preprint2012arXiv

Volumetric Mapping of Genus Zero Objects via Mass Preservation

In this work, we present a technique to map any genus zero solid object onto a hexahedral decomposition of a solid cube. This problem appears in many applications ranging from finite element methods to visual tracking. From this, one can then hopefully utilize the proposed technique for shape analysis, registration, as well as other related computer graphics tasks. More importantly, given that we seek to establish a one-to-one correspondence of an input volume to that of a solid cube, our algorithm can naturally generate a quality hexahedral mesh as an output. In addition, we constrain the mapping itself to be volume preserving allowing for the possibility of further mesh simplification. We demonstrate our method both qualitatively and quantitatively on various 3D solid models