Total variation optimization for imaging through turbid media with transmission matrix

Abstract. With the transmission matrix (TM) of the whole optical system measured, the image of the object behind a turbid medium can be recovered from its speckle field by means of an image reconstruction algorithm. Instead of Tikhonov regularization algorithm (TRA), the total variation minimization by augmented Lagrangian and alternating direction algorithms (TVAL3) is introduced to recover object images. As a total variation (TV)-based approach, TVAL3 allows to effectively damp more noise and preserve more edges compared with TRA, thus providing more outstanding image quality. Different levels of detector noise and TM-measurement noise are successively added to analyze the antinoise performance of these two algorithms. Simulation results show that TVAL3 is able to recover more details and suppress more noise than TRA under different noise levels, thus providing much more excellent image quality. Furthermore, whether it be detector noise or TM-measurement noise, the reconstruction images obtained by TVAL3 at SNR=15  dB are far superior to those by TRA at SNR=50  dB.

[1]  C. Beenakker.  Random-matrix theory of quantum transport , 1996, cond-mat/9612179.

[2]  E. G. van Putten,et al.  Scattering lens resolves sub-100 nm structures with visible light. , 2011, Physical review letters.

[3]  G. Pedrini,et al.  Phase retrieval of arbitrary complex-valued fields through aperture-plane modulation , 2007 .

[4]  Sylvain Gigan,et al.  Image transmission through an opaque material. , 2010, Nature communications.

[5]  J R Fienup,et al.  Phase retrieval algorithms: a comparison. , 1982, Applied optics.

[6]  Florent Krzakala,et al.  Reference-less measurement of the transmission matrix of a highly scattering material using a DMD and phase retrieval techniques. , 2015, Optics express.

[7]  Yin Zhang,et al.  Compressive sensing for 3d data processing tasks: applications, models and algorithms , 2011 .

[8]  H. H. Rachford,et al.  The Numerical Solution of Parabolic and Elliptic Differential Equations , 1955 .

[9]  Emmanuel J. Candès,et al.  NESTA: A Fast and Accurate First-Order Method for Sparse Recovery , 2009, SIAM J. Imaging Sci..

[10]  J. Borwein,et al.  Two-Point Step Size Gradient Methods , 1988 .

[11]  I. Yamaguchi,et al.  Three-dimensional microscopy with phase-shifting digital holography. , 1998, Optics letters.

[12]  J. Bertolotti,et al.  Speckle correlation resolution enhancement of wide-field fluorescence imaging , 2015 .

[13]  Junfeng Yang,et al.  A New Alternating Minimization Algorithm for Total Variation Image Reconstruction , 2008, SIAM J. Imaging Sci..

[14]  Junfeng Yang,et al.  A Fast Alternating Direction Method for TVL1-L2 Signal Reconstruction From Partial Fourier Data , 2010, IEEE Journal of Selected Topics in Signal Processing.

[15]  A Tikhonov,et al.  Solution of Incorrectly Formulated Problems and the Regularization Method , 1963 .

[16]  S. Popoff,et al.  Controlling light through optical disordered media: transmission matrix approach , 2011, 1107.5285.

[17]  Yefeng Guan,et al.  Image restoration through thin turbid layers by correlation with a known object. , 2013, Optics express.

[18]  Feng,et al.  Memory effects in propagation of optical waves through disordered media. , 1988, Physical review letters.

[19]  José M. Bioucas-Dias,et al.  A New TwIST: Two-Step Iterative Shrinkage/Thresholding Algorithms for Image Restoration , 2007, IEEE Transactions on Image Processing.

[20]  Laurent Daudet,et al.  Imaging With Nature: Compressive Imaging Using a Multiply Scattering Medium , 2013, Scientific Reports.

[21]  Michael S Feld,et al.  Overcoming the diffraction limit using multiple light scattering in a highly disordered medium. , 2011, Physical review letters.

[22]  Yin Zhang,et al.  An efficient augmented Lagrangian method with applications to total variation minimization , 2013, Computational Optimization and Applications.

[23]  Tengfei Wu,et al.  Simulation on light refocusing through a highly scattering turbid medium using circular Gaussian distribution model , 2013 .

[24]  Tengfei Wu,et al.  An adaptive total variation image reconstruction method for speckles through disordered media , 2013, Optics & Photonics - Optical Engineering + Applications.

[25]  D. Conkey,et al.  Genetic algorithm optimization for focusing through turbid media in noisy environments. , 2012, Optics express.

[26]  Vicente Durán,et al.  Image transmission through dynamic scattering media by single-pixel photodetection. , 2014, Optics express.

[27]  M. Hestenes.  Multiplier and gradient methods , 1969 .

[28]  A. Mosk,et al.  Focusing coherent light through opaque strongly scattering media. , 2007, Optics letters.

[29]  M. Fink,et al.  Non-invasive single-shot imaging through scattering layers and around corners via speckle correlations , 2014, Nature Photonics.

[30]  J. Bertolotti,et al.  Non-invasive imaging through opaque scattering layers , 2012, Nature.

[31]  M. J. D. Powell,et al.  A method for nonlinear constraints in minimization problems , 1969 .

[32]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[33]  Chengbo Li.  An efficient algorithm for total variation regularization with applications to the single pixel camera and compressive sensing , 2010 .

[34]  A. Mosk,et al.  Phase control algorithms for focusing light through turbid media , 2007, 0710.3295.

[35]  Junfeng Yang,et al.  A Fast Algorithm for Edge-Preserving Variational Multichannel Image Restoration , 2009, SIAM J. Imaging Sci..

[36]  Feng,et al.  Correlations and fluctuations of coherent wave transmission through disordered media. , 1988, Physical review letters.

[37]  O. Katz,et al.  Looking around corners and through thin turbid layers in real time with scattered incoherent light , 2012, Nature Photonics.

[38]  S. Popoff,et al.  Measuring the transmission matrix in optics: an approach to the study and control of light propagation in disordered media. , 2009, Physical review letters.

[39]  Eero P. Simoncelli,et al.  Image quality assessment: from error visibility to structural similarity , 2004, IEEE Transactions on Image Processing.

[40]  William W. Hager,et al.  A Nonmonotone Line Search Technique and Its Application to Unconstrained Optimization , 2004, SIAM J. Optim..

[41]  Jonghee Yoon,et al.  Measuring optical transmission matrices by wavefront shaping. , 2015, Optics express.