1,100 research outputs found
Matrix probing: a randomized preconditioner for the wave-equation Hessian
This paper considers the problem of approximating the inverse of the
wave-equation Hessian, also called normal operator, in seismology and other
types of wave-based imaging. An expansion scheme for the pseudodifferential
symbol of the inverse Hessian is set up. The coefficients in this expansion are
found via least-squares fitting from a certain number of applications of the
normal operator on adequate randomized trial functions built in curvelet space.
It is found that the number of parameters that can be fitted increases with the
amount of information present in the trial functions, with high probability.
Once an approximate inverse Hessian is available, application to an image of
the model can be done in very low complexity. Numerical experiments show that
randomized operator fitting offers a compelling preconditioner for the
linearized seismic inversion problem.Comment: 21 pages, 6 figure
Regularized Newton Methods for X-ray Phase Contrast and General Imaging Problems
Like many other advanced imaging methods, x-ray phase contrast imaging and
tomography require mathematical inversion of the observed data to obtain
real-space information. While an accurate forward model describing the
generally nonlinear image formation from a given object to the observations is
often available, explicit inversion formulas are typically not known. Moreover,
the measured data might be insufficient for stable image reconstruction, in
which case it has to be complemented by suitable a priori information. In this
work, regularized Newton methods are presented as a general framework for the
solution of such ill-posed nonlinear imaging problems. For a proof of
principle, the approach is applied to x-ray phase contrast imaging in the
near-field propagation regime. Simultaneous recovery of the phase- and
amplitude from a single near-field diffraction pattern without homogeneity
constraints is demonstrated for the first time. The presented methods further
permit all-at-once phase contrast tomography, i.e. simultaneous phase retrieval
and tomographic inversion. We demonstrate the potential of this approach by
three-dimensional imaging of a colloidal crystal at 95 nm isotropic resolution.Comment: (C)2016 Optical Society of America. One print or electronic copy may
be made for personal use only. Systematic reproduction and distribution,
duplication of any material in this paper for a fee or for commercial
purposes, or modifications of the content of this paper are prohibite
- …