Backward discrete wave field propagation modeling as an inverse problem: toward perfect reconstruction of wave field distributions.

We consider reconstruction of a wave field distribution in an input/object plane from data in an output/diffraction (sensor) plane. We provide digital modeling both for the forward and backward wave field propagation. A novel algebraic matrix form of the discrete diffraction transform (DDT) originated in Katkovnik et al. [Appl. Opt. 47, 3481 (2008)] is proposed for the forward modeling that is aliasing free and precise for pixelwise invariant object and sensor plane distributions. This "matrix DDT" is a base for formalization of the object wave field reconstruction (backward propagation) as an inverse problem. The transfer matrices of the matrix DDT are used for calculations as well as for the analysis of conditions when the perfect reconstruction of the object wave field distribution is possible. We show by simulation that the developed inverse propagation algorithm demonstrates an improved accuracy as compared with the standard convolutional and discrete Fresnel transform algorithms.

[1]  G. Sherman.  Integral-Transform Formulation of Diffraction Theory* , 1967 .

[2]  J. Goodman.  Introduction to Fourier optics , 1969 .

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

[4]  Mario Bertero,et al.  Introduction to Inverse Problems in Imaging , 1998 .

[5]  Alan V. Oppenheim,et al.  Discrete-time signal processing (2nd ed.) , 1999 .

[6]  I. Yamaguchi,et al.  Image formation in phase-shifting digital holography and applications to microscopy. , 2001, Applied optics.

[7]  Thierry Blu,et al.  Fresnelets: new multiresolution wavelet bases for digital holography , 2003, IEEE Trans. Image Process..

[8]  Fabrice Labeau,et al.  Discrete Time Signal Processing , 2004 .

[9]  T. Kreis.  Handbook of Holographic Interferometry: Optical and Digital Methods , 2004 .

[10]  L. Onural.  Some mathematical properties of the uniformly sampled quadratic phase function and associated issues in digital Fresnel diffraction simulations , 2004 .

[11]  John T. Sheridan,et al.  Generalizing, optimizing, and inventing numerical algorithms for the fractional Fourier, Fresnel, and linear canonical transforms. , 2005, Journal of the Optical Society of America. A, Optics, image science, and vision.

[12]  Emmanuel J. Candès,et al.  Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information , 2004, IEEE Transactions on Information Theory.

[13]  Anbo Wang,et al.  Fast-Fourier-transform based numerical integration method for the Rayleigh-Sommerfeld diffraction formula. , 2006, Applied optics.

[14]  I. Aizenberg,et al.  Discrete Generalized Fresnel Functions and Transforms in an Arbitrary Discrete Basis , 2006, IEEE Transactions on Signal Processing.

[15]  Jaakko Astola,et al.  Local Approximation Techniques in Signal and Image Processing (SPIE Press Monograph Vol. PM157) , 2006 .

[16]  L. Onural.  Exact analysis of the effects of sampling of the scalar diffraction field. , 2007, Journal of the Optical Society of America. A, Optics, image science, and vision.

[17]  J. Astola,et al.  Discrete diffraction transform for propagation, reconstruction, and design of wavefield distributions. , 2008, Applied optics.

[18]  Don H. Johnson,et al.  Statistical Signal Processing , 2009, Encyclopedia of Biometrics.