Document not found! Please try again

Sparse reconstruction in digital holography

7 downloads 262 Views 3MB Size Report
Sparse reconstruction in digital holography. Loïc Denis. Observatoire de Lyon joint work with Corinne Fournier, Ferréol Soulez and Eric Thiébaut. Loïc Denis.
Sparse reconstruction in digital holography Loïc Denis Observatoire de Lyon

joint work with Corinne Fournier, Ferréol Soulez and Eric Thiébaut

Loïc Denis

Sparse reconstruction in digital holography

Context Digital holography: 2 possible meanings 1

Generate numerically a hologram for 3D visualization:

Loïc Denis

Sparse reconstruction in digital holography

Context Digital holography: 2 possible meanings 1

Generate numerically a hologram for 3D visualization:

Loïc Denis

Sparse reconstruction in digital holography

Context Digital holography: 2 possible meanings 1

Generate numerically a hologram for 3D visualization:

2

Record a hologram on a digital camera and reconstruct it numerically:

Applications: - metrology (particles, surface deformation) - (fast) imaging of semi-transparent objects Loïc Denis

Sparse reconstruction in digital holography

Context (cont) Digital holography: 3D+t particle imaging Particle tracking is an essential technique in experimental fluid dynamics. Holography has long been used to instantaneously measure particle sizes and 3D locations:

Loïc Denis

Sparse reconstruction in digital holography

Context (cont) Digital holography: 3D+t particle imaging Particle tracking is an essential technique in experimental fluid dynamics. Holography has long been used to instantaneously measure particle sizes and 3D locations. With digital holography time-resolved measurements can be done:

Loïc Denis

Sparse reconstruction in digital holography

What is measured in particle holography Direct model Diffraction can be modeled as a convolution with a (complexed valued) chirp kernel (Fresnel’s approximation). A spherical opaque particle creates concentric fringes on the hologram:

Loïc Denis

Sparse reconstruction in digital holography

What is measured in particle holography Direct model Diffraction can be modeled as a convolution with a (complexed valued) chirp kernel (Fresnel’s approximation). A spherical opaque particle creates concentric fringes on the hologram. Inter-particle interferences can be neglected → linear direct model: Direct model: d = c 1 − H · ϑ + , with d the data, c an offset, H the dictionnary of diffraction patterns, ϑ the unknown particle distribution. Loïc Denis

Sparse reconstruction in digital holography

What is measured in particle holography Direct model Diffraction can be modeled as a convolution with a (complexed valued) chirp kernel (Fresnel’s approximation). A spherical opaque particle creates concentric fringes on the hologram . Inter-particle interferences can be neglected → linear direct model: Direct model: ¯ · ϑ + , d¯ = H with proper centering.

Loïc Denis

Sparse reconstruction in digital holography

The detection/reconstruction problems 1- Detection problem From a particle hologram, detect all particles (3D location + size)

2- Reconstruction problem From a hologram, reconstruct the 2D or 3D transmittance distrib.

Loïc Denis

Sparse reconstruction in digital holography

Detection: a CLEAN/matching-pursuit approach

Loïc Denis

Sparse reconstruction in digital holography

Detection: a CLEAN/matching-pursuit approach → Ferréol Soulez’s PhD thesis 1- Global detection The model is shift-invariant → correlationbased detection of the “best matching” diffraction pattern. 2- Local optimization The underlying “dictionary” should ideally be continuous → local (LevenbergMarquardt) optimization. 3- Residuals update Subtract the fitted diffraction pattern, and re-iterate steps 1-3 to detect another particle. Loïc Denis

Sparse reconstruction in digital holography

Detection: a CLEAN/matching-pursuit approach (cont) Results: 1- 3D+t reconstruction. 2- Out-of-field particle detection:

→ accuracy estimation with Cramèr-Rao bounds (cf. Corinne’s talk) Loïc Denis

Sparse reconstruction in digital holography

Reconstruction: minimization of sparsity-promoting `1 norm

Principle: `1 norm minimization is widely used to get a solution that contains many zeros (i.e., a sparse solution). The reconstruction problem is recast into a (non-smooth) minimization problem: ˆ sparse = arg min ϑ ϑ

1 kHϑ ¯ 2

− d¯ k22 + τ kϑk1

Using a weighted `1 norm is useful to account for non-normalized dictionaries.

Loïc Denis

Sparse reconstruction in digital holography

Reconstruction: minimization of sparsity-promoting `1 norm Results:

Loïc Denis

Sparse reconstruction in digital holography

Conclusion

1- Hologram reconstruction as a linear inverse problem 2- CLEAN/matching-pursuit approach to particle detection gave unprecedented results 3- Global optimization with sparsity-promoting `1 norm improves on classical (diffraction-based) reconstruction

Loïc Denis

Sparse reconstruction in digital holography

Ferreol Soulez, Loic Denis, Corinne Fournier, Eric Thiebaut, and Charles Goepfert. Inverse-problem approach for particle digital holography: accurate location based on local optimization. Journal of the Optical Society of America A, 24(4):1164–1171, April 2007. Ferreol Soulez, Loic Denis, Eric Thiebaut, Corinne Fournier, and Charles Goepfert. Inverse problem approach in particle digital holography: out-of-field particle detection made possible. Journal of the Optical Society of America A, 24(12):3708–3716, December 2007. L Denis, D A Lorenz, and D Trede. Greedy solution of ill-posed problems: error bounds and exact inversion. Inverse Problems, 25(11):115017, 2009. Loic Denis, Dirk Lorenz, Eric Thiebaut, Corinne Fournier, and Dennis Trede. Inline hologram reconstruction with sparsity constraints. Optics Letters, 34(22):3475–3477, November 2009.

Loïc Denis

Sparse reconstruction in digital holography