Graph- and finite element-based total variation models for the inverse problem in diffuse optical tomography

Total variation (TV) is a powerful regularization method that has been widely applied in different imaging applications, but is difficult to apply to diffuse optical tomography (DOT) image reconstruction (inverse problem) due to complex and unstructured geometries, non-linearity of the data fitting and regularization terms, and non-differentiability of the regularization term. We develop several approaches to overcome these difficulties by: i) defining discrete differential operators for unstructured geometries using both finite element and graph representations; ii) developing an optimization algorithm based on the alternating direction method of multipliers (ADMM) for the non-differentiable and non-linear minimization problem; iii) investigating isotropic and anisotropic variants of TV regularization, and comparing their finite element- and graph-based implementations. These approaches are evaluated on experiments on simulated data and real data acquired from a tissue phantom. Our results show that both FEM and graph-based TV regularization is able to accurately reconstruct both sparse and non-sparse distributions without the over-smoothing effect of Tikhonov regularization and the over-sparsifying effect of L$_1$ regularization. The graph representation was found to out-perform the FEM method for low-resolution meshes, and the FEM method was found to be more accurate for high-resolution meshes.Comment: 24 pages, 11 figures. Reviced version includes revised figures and improved clarit

Magnetic susceptibility anisotropy of myocardium imaged by cardiovascular magnetic resonance reflects the anisotropy of myocardial filament α-helix polypeptide bonds.

BackgroundA key component of evaluating myocardial tissue function is the assessment of myofiber organization and structure. Studies suggest that striated muscle fibers are magnetically anisotropic, which, if measurable in the heart, may provide a tool to assess myocardial microstructure and function.MethodsTo determine whether this weak anisotropy is observable and spatially quantifiable with cardiovascular magnetic resonance (CMR), both gradient-echo and diffusion-weighted data were collected from intact mouse heart specimens at 9.4 Tesla. Susceptibility anisotropy was experimentally calculated using a voxelwise analysis of myocardial tissue susceptibility as a function of myofiber angle. A myocardial tissue simulation was developed to evaluate the role of the known diamagnetic anisotropy of the peptide bond in the observed susceptibility contrast.ResultsThe CMR data revealed that myocardial tissue fibers that were parallel and perpendicular to the magnetic field direction appeared relatively paramagnetic and diamagnetic, respectively. A linear relationship was found between the magnetic susceptibility of the myocardial tissue and the squared sine of the myofiber angle with respect to the field direction. The multi-filament model simulation yielded susceptibility anisotropy values that reflected those found in the experimental data, and were consistent that this anisotropy decreased as the echo time increased.ConclusionsThough other sources of susceptibility anisotropy in myocardium may exist, the arrangement of peptide bonds in the myofilaments is a significant, and likely the most dominant source of susceptibility anisotropy. This anisotropy can be further exploited to probe the integrity and organization of myofibers in both healthy and diseased heart tissue

Weighted Mean Curvature

In image processing tasks, spatial priors are essential for robust computations, regularization, algorithmic design and Bayesian inference. In this paper, we introduce weighted mean curvature (WMC) as a novel image prior and present an efficient computation scheme for its discretization in practical image processing applications. We first demonstrate the favorable properties of WMC, such as sampling invariance, scale invariance, and contrast invariance with Gaussian noise model; and we show the relation of WMC to area regularization. We further propose an efficient computation scheme for discretized WMC, which is demonstrated herein to process over 33.2 giga-pixels/second on GPU. This scheme yields itself to a convolutional neural network representation. Finally, WMC is evaluated on synthetic and real images, showing its superiority quantitatively to total-variation and mean curvature.Comment: 12 page