7,012 research outputs found
Boundary length of reconstructions in discrete tomography
We consider possible reconstructions of a binary image of which the row and
column sums are given. For any reconstruction we can define the length of the
boundary of the image. In this paper we prove a new lower bound on the length
of this boundary. In contrast to simple bounds that have been derived
previously, in this new lower bound the information of both row and column sums
is combined
Geometric reconstruction methods for electron tomography
Electron tomography is becoming an increasingly important tool in materials
science for studying the three-dimensional morphologies and chemical
compositions of nanostructures. The image quality obtained by many current
algorithms is seriously affected by the problems of missing wedge artefacts and
nonlinear projection intensities due to diffraction effects. The former refers
to the fact that data cannot be acquired over the full tilt range;
the latter implies that for some orientations, crystalline structures can show
strong contrast changes. To overcome these problems we introduce and discuss
several algorithms from the mathematical fields of geometric and discrete
tomography. The algorithms incorporate geometric prior knowledge (mainly
convexity and homogeneity), which also in principle considerably reduces the
number of tilt angles required. Results are discussed for the reconstruction of
an InAs nanowire
Edge-promoting reconstruction of absorption and diffusivity in optical tomography
In optical tomography a physical body is illuminated with near-infrared light
and the resulting outward photon flux is measured at the object boundary. The
goal is to reconstruct internal optical properties of the body, such as
absorption and diffusivity. In this work, it is assumed that the imaged object
is composed of an approximately homogeneous background with clearly
distinguishable embedded inhomogeneities. An algorithm for finding the maximum
a posteriori estimate for the absorption and diffusion coefficients is
introduced assuming an edge-preferring prior and an additive Gaussian
measurement noise model. The method is based on iteratively combining a lagged
diffusivity step and a linearization of the measurement model of diffuse
optical tomography with priorconditioned LSQR. The performance of the
reconstruction technique is tested via three-dimensional numerical experiments
with simulated measurement data.Comment: 18 pages, 6 figure
A parametric level-set method for partially discrete tomography
This paper introduces a parametric level-set method for tomographic
reconstruction of partially discrete images. Such images consist of a
continuously varying background and an anomaly with a constant (known)
grey-value. We represent the geometry of the anomaly using a level-set
function, which we represent using radial basis functions. We pose the
reconstruction problem as a bi-level optimization problem in terms of the
background and coefficients for the level-set function. To constrain the
background reconstruction we impose smoothness through Tikhonov regularization.
The bi-level optimization problem is solved in an alternating fashion; in each
iteration we first reconstruct the background and consequently update the
level-set function. We test our method on numerical phantoms and show that we
can successfully reconstruct the geometry of the anomaly, even from limited
data. On these phantoms, our method outperforms Total Variation reconstruction,
DART and P-DART.Comment: Paper submitted to 20th International Conference on Discrete Geometry
for Computer Imager
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 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
Study of noise effects in electrical impedance tomography with resistor networks
We present a study of the numerical solution of the two dimensional
electrical impedance tomography problem, with noisy measurements of the
Dirichlet to Neumann map. The inversion uses parametrizations of the
conductivity on optimal grids. The grids are optimal in the sense that finite
volume discretizations on them give spectrally accurate approximations of the
Dirichlet to Neumann map. The approximations are Dirichlet to Neumann maps of
special resistor networks, that are uniquely recoverable from the measurements.
Inversion on optimal grids has been proposed and analyzed recently, but the
study of noise effects on the inversion has not been carried out. In this paper
we present a numerical study of both the linearized and the nonlinear inverse
problem. We take three different parametrizations of the unknown conductivity,
with the same number of degrees of freedom. We obtain that the parametrization
induced by the inversion on optimal grids is the most efficient of the three,
because it gives the smallest standard deviation of the maximum a posteriori
estimates of the conductivity, uniformly in the domain. For the nonlinear
problem we compute the mean and variance of the maximum a posteriori estimates
of the conductivity, on optimal grids. For small noise, we obtain that the
estimates are unbiased and their variance is very close to the optimal one,
given by the Cramer-Rao bound. For larger noise we use regularization and
quantify the trade-off between reducing the variance and introducing bias in
the solution. Both the full and partial measurement setups are considered.Comment: submitted to Inverse Problems and Imagin
Inverse scattering for reflection intensity phase microscopy
Reflection phase imaging provides label-free, high-resolution characterization of biological samples, typically using interferometric-based techniques. Here, we investigate reflection phase microscopy from intensity-only measurements under diverse illumination. We evaluate the forward and inverse scattering model based on the first Born approximation for imaging scattering objects above a glass slide. Under this design, the measured field combines linear forward-scattering and height-dependent nonlinear back-scattering from the object that complicates object phase recovery. Using only the forward-scattering, we derive a linear inverse scattering model and evaluate this model's validity range in simulation and experiment using a standard reflection microscope modified with a programmable light source. Our method provides enhanced contrast of thin, weakly scattering samples that complement transmission techniques. This model provides a promising development for creating simplified intensity-based reflection quantitative phase imaging systems easily adoptable for biological research.https://arxiv.org/abs/1912.07709Accepted manuscrip
- …