Regularized non-convex image reconstruction in digital holographic microscopy

Inverse problem approaches for image reconstruction can improve resolution recovery over spatial filtering methods while reducing interference artifacts in digital off-axis holography. Prior works implemented explicit regularization operators in the image space and were only able to match intensity measurements approximatively. As a consequence, convergence to a strictly compatible solution was not possible. In this paper, we replace the non-convex image reconstruction problem for a sequence of surrogate convex problems. An iterative numerical solver is designed using a simple projection operator in the data domain and a Nesterov acceleration of the simultaneous Kaczmarz method. For regularization, the complex-valued object wavefield image is represented in the multiresolution CDF 9/7 wavelet domain and an energy-weighted preconditioning promotes minimum-norm solutions. Experiments demonstrate improved resolution recovery and reduced spurious artifacts in reconstructed images. Furthermore, the method is resilient to additive Gaussian noise and subsampling of intensity measurements. c © 2017 Optical Society of America OCIS codes: (090.1995) Digital holography; (100.3190) Inverse problems; (100.3010) Image reconstruction techniques; (100.5070) Phase retrieval; (100.7410) Wavelets; (110.1758) Computational imaging. References and links 1. M. K. Kim, “Principles and techniques of digital holographic microscopy,” SPIE Rev. 1, 8005 (2010). 2. E. Cuche, P. Marquet, and C. Depeursinge, “Simultaneous amplitude-contrast and quantitative phase-contrast microscopy by numerical reconstruction of fresnel off-axis holograms,” Appl. Opt. 38, 6994–7001 (1999). 3. N. Pavillon, C. S. Seelamantula, J. Kühn, M. Unser, and C. Depeursinge, “Suppression of the zero-order term in off-axis digital holography through nonlinear filtering,” Appl. Opt. 48, H186–H195 (2009). 4. C. S. Seelamantula, N. Pavillon, C. Depeursinge, and M. Unser, “Exact complex-wave reconstruction in digital holography,” J. Opt. Soc. Am. A 28, 983–992 (2011). 5. S. Sotthivirat and J. A. Fessler, “Reconstruction from digital holograms by statistical methods,” Conference Record of the Thirty-Seventh Asilomar Conference on Signals, Systems and Computers 2, 1928–1932 (2003). 6. S. Sotthivirat and J. A. Fessler, “Penalized-likelihood image reconstruction for digital holography,” J. Opt. Soc. Am. A 21, 737–750 (2004). 7. J. A. Fessler and S. Sotthivirat, “Simplified digital holographic reconstruction using statistical methods,” 2004 International Conference on Image Processing (ICIP’04) 4, 2435–2438 (2004). 8. A. Bourquard, N. Pavillon, E. Bostan, C. Depeursinge, and M. Unser, “A practical inverse-problem approach to digital holographic reconstruction,” Opt. Express 21, 3417–3433 (2013). 9. L. I. Rudin, S. Osher, and E. Fatemi, “Nonlinear total variation based noise removal algorithms,” Phys. D: Nonlin. Phenom. 60, 259–268 (1992). 10. Y. Rivenson, A. Stern, and B. Javidi, “Compressive Fresnel holography,” J. Displ. Technol. 6, 506–509 (2010). 11. D. Donoho, “Compressed sensing,” IEEE Trans. Inform. Theory 52, 1289–1306 (2006). 12. E. J. Candès, M. B. Wakin, and S. Boyd, “Enhancing sparsity by reweighted l1 minimization,” J. Fourier Anal. Appl. 14, 877–905 (2008). 13. M. M. Marim, M. Atlan, E. Angelini, and J.-C. Olivo-Marin, “Compressed sensing with off-axis frequency-shifting holography,” Opt. Lett. 35, 871–873 (2010). 14. A. Cohen, I. Daubechies, and J.-C. Feauveau, “Biorthogonal bases of compactly supported wavelets,” Commun. Pure Appl. Math. 45, 485–560 (1992). Vol. 25, No. 14 | 10 Jul 2017 | OPTICS EXPRESS 16491 #290649 https://doi.org/10.1364/OE.25.016491 Journal © 2017 Received 14 Mar 2017; revised 10 May 2017; accepted 10 May 2017; published 3 Jul 2017 15. S. Bettens, H. Yan, D. Blinder, H. Ottevaere, C. Schretter, and P. Schelkens, “Studies on the sparsifying operator in compressive digital holography,” Opt. Express (2017, submitted). 16. M. Liebling, T. Blu, and M. Unser, “Fresnelets: New multiresolution wavelet bases for digital holography,” IEEE Trans. Image Process. 12, 29–43 (2003). 17. T.-C. Poon, Digital holography and three-dimensional display: Principles and Applications (Springer Science & Business Media, 2006). 18. J. W. Goodman, Introduction to Fourier Optics (Roberts and Company Publishers, 2004), 3rd ed. 19. S. G. Mallat, “A Theory for Multiresolution Signal Decomposition: The Wavelet Representation,” IEEE Trans. Pattern Anal. Mach. Intel. 11, 674–693 (1989). 20. A. A. Petrosian and F. G. Meyer, Wavelets in signal and image analysis: from theory to practice, vol. 19 (Springer Science & Business Media, 2013). 21. P. Schelkens, A. Skodras, and T. Ebrahimi, The JPEG 2000 Suite (Wiley Publishing, 2009). 22. J. R. Fienup, “Phase retrieval algorithms: a comparison,” Appl. Opt. 21, 2758–2769 (1982). 23. A. H. Andersen and A. C. Kak, “Simultaneous algebraic reconstruction technique (SART): a superior implementation of the ART algorithm,” Ultrason. Imag. 6, 81–94 (1984). 24. S. Kaczmarz, “Angenäherte Auflösung von Systemen linearer Gleichungen,” Bulletin International de l’Académie Polonaise des Sciences et des Lettres 35, 355–357 (1937). 25. B. T. Polyak and A. B. Juditsky, “Acceleration of stochastic approximation by averaging,” SIAM J. Contr. Optim. 30, 838–855 (1992). 26. D. Needell, N. Srebro, and R. Ward, “Stochastic gradient descent and the randomized Kaczmarz algorithm,” arXiv preprint arXiv:1310.5715 (2013). 27. K. Wei, “Phase retrieval via Kaczmarz methods,” arXiv preprint arXiv:1502.01822 (2015). 28. Y. Nesterov, “A method for unconstrained convex minimization problem with the rate of convergence O(1/k2),” Doklady an SSSR 269, 543–547 (1983). 29. Y. Nesterov, “Smooth minimization of non-smooth functions,” Math. Program. 103, 127–152 (2005). 30. D. Kim and J. A. Fessler, “Optimized first-order methods for smooth convex minimization,” Math. Program. 159, 81–107 (2016). 31. C. Schretter and H. Niederreiter, “A direct inversion method for non-uniform quasi-random point sequences,” Monte Carlo Meth. Appl. 19, 1–9 (2013). 32. W. Sweldens, “The lifting scheme: A custom-design construction of biorthogonal wavelets,” Appl. Comput. Harmon. Anal. 3, 186 – 200 (1996). 33. I. Daubechies and W. Sweldens, “Factoring wavelet transforms into lifting steps,” J. Fourier Anal. Appl. 4, 247–269 (1998).