Damien P. Kelly

Numerical calculation of the Fresnel transform

In this paper, we address the problem of calculating Fresnel diffraction integrals using a finite number of uniformly spaced samples. General and simple sampling rules of thumb are derived that allow the user to calculate the distribution for any propagation distance. It is shown how these rules can be extended to fast-Fourier-transform-based algorithms to increase calculation efficiency. A comparison with other theoretical approaches is made.

Three-dimensional static speckle fields.Part I. Theory and numerical investigation

When monochromatic light is scattered from an optically rough surface a complicated three-dimensional (3D) field is generated. These fields are often described by reference to the 3D volume (extent) of their speckles, leading to the definition of lateral (x,y) and longitudinal speckle sizes (z). For reasons of mathematical simplicity the longitudinal speckle size is often derived by examining the decorrelation of the speckle field for a single point lying on axis, i.e., x=y=0, and this size is generally assumed to be representative for other speckles that lie further off-axis. Some recent theoretical results, however, indicate that in fact longitudinal speckle size gets smaller as the observation position moves to off-axis spatial locations. In this paper (Part I), we review the physical argument leading to this conclusion and support this analysis with a series of robust numerical simulations. We discuss, in some detail, computational issues that arise when simulating the propagation of speckle fields numerically, showing that the spectral method is not a suitable propagation algorithm when the autocorrelation of the scattering surface is assumed to be delta correlated. In Part II [J. Opt. Soc. Am. A28, 1904 (2011)] of this paper, experimental results are provided that exhibit the predicted variation of longitudinal speckle size as a function of position in x and y. The results are not only of theoretical interest but have practical implications, and in Part II a method for locating the optical system axis is proposed and experimentally demonstrated.

Numerical solution of nonparaxial scalar diffraction integrals for focused fields

In this paper, we present sampling conditions for fast-Fourier-transform-based field propagations. The input field and the propagation kernel are analyzed in a combined manner to derive sampling criteria that guarantee accurate calculation results in the output plane. These sampling criteria are also applicable to the propagation of general fields. For focal field calculations, geometrical optics is used to obtain a priori knowledge about the input and output fields. This a priori knowledge is used to determine an optimum balance between computational load and calculation accuracy. In a numerical example, correct results are obtained even though both the input field and the propagation kernel are sampled below the Nyquist rate. Finally, we show how chirp z-transform-based zoom-algorithms may be analyzed using the same techniques.

Speckle orientation in paraxial optical systems

The statistical properties of speckles in paraxial optical systems depend on the system parameters. In particular, the speckle orientation and the lateral dependence (x and y) of the longitudinal speckle size can vary significantly. For example, the off-axis longitudinal correlation length remains equal to the on-axis size for speckles in a Fourier transform system, while it decreases dramatically as the observation position moves off axis in a Fresnel system. In this paper, we review the speckle correlation function in general linear canonical transform (LCT) systems, clearly demonstrating that speckle properties can be controlled by introducing different optical components, i.e., lenses and sections of free space. Using a series of numerical simulations, we examine how the correlation function changes for some typical LCT systems. The integrating effect of the camera pixel and the impact this has on the measured first- and second-order statistics of the speckle intensities is also examined theoretically. A series of experimental results are then presented to confirm several of these predictions. First, the effect the pixel size has on the measured first-order speckle statistics is demonstrated, and second, the orientation of speckles in a Fourier transform system is measured, showing that the speckles lie parallel to the optical axis.

Fast nonparaxial scalar focal field calculations

An efficient algorithm for calculating nonparaxial scalar field distributions in the focal region of a lens is discussed. The algorithm is based on fast Fourier transform implementations of the first Rayleigh-Sommerfeld diffraction integral and assumes that the input field at the pupil plane has a larger extent than the field in the focal region. A sampling grid is defined over a finite region in the output plane and referred to as a tile. The input field is divided into multiple separate spatial regions of the size of the output tile. Finally, the input tiles are added coherently to form a summed tile, which is propagated to the output plane. Since only a single tile is propagated, there are significant reductions of computational load and memory requirements. This method is combined either with a subpixel sampling technique or with a chirp z-transform to realize smaller sampling intervals in the output plane than in the input plane. For a given example the resulting methods enable a speedup of approximately 800× in comparison to the normal angular spectrum method, while the memory requirements are reduced by more than 99%.

Speckle suppression by doubly scattering systems

Speckle suppression in a two-diffuser system is examined. An analytical expression for the speckle space-time correlation function is derived, so that the speckle suppression mechanism can be investigated statistically. The grain size of the speckle field illuminating the second diffuser has a major impact on the speckle contrast after temporal averaging. It is shown that, when both the diffusers are rotating, the one with the lower rotating speed determines the period of the speckle correlation function. The coherent length of the averaged speckle intensity is shown to equal the mean speckle size of the individual speckle pattern before averaging. Numerical and experimental results are presented to verify our analysis in the context of speckle reduction.

K speckle: space–time correlation function of doubly scattered light in an imaging system

The scattering of coherent monochromatic light at an optically rough surface, such as a diffuser, produces a speckle field, which is usually described by reference to its statistical properties. For example, the real and imaginary parts of a fully developed speckle field can be modeled as a random circular Gaussian process. When such a speckle field is used to illuminate a second diffuser, the statistics of the resulting doubly scattered field are in general no longer Gaussian, but rather follow a K distribution. In this paper we determine the space-time correlation function of such a doubly scattered speckle field that has been imaged by a single lens system. A space-time correlation function is derived that contains four separate terms; similar to the Gaussian case it contains an average DC term and a fluctuating AC term. However, in addition there are two terms that are related to contributions from each of the diffusers independently. We examine how our space-time correlation function varies as the diffusers are rotated at different speeds and as the point spread function of the imaging system is changed. A series of numerical simulations are used to confirm different aspects of the theoretical analysis. We then finish with a discussion of our results and some potential applications, including controlling spatial coherence and speckle reduction.

Experimental evaluation of inline free-space holography systems

It is important to be able to quantify, theoretically and experimentally, the performance of coherent digital systems, so that their suitability for a given metrology application can be assessed. Here, a free-space inline digital holographic system is investigated. To isolate the scattered object field, phase-shifting interferometry (PSI) techniques are used. Several sequential holographic measurements are made, where the phase of the reference field is stepped by a known amount relative to the scattered object field between captures. Under ideal conditions such as noise-free electronics, vibration-free environments, and perfect reference and illuminating object waves, this system will be diffraction limited. However, real systems suffer from experimental error and noise effects. In this paper, we examine a PSI digital holographic imaging system considering all prominent error sources. An experimental metric is defined that quantifies how far from the theoretical ideal a real system is performing. By carefully optimizing our system, following our recommended guidelines, we approach diffraction limited imaging, surpassing the Nyquist sampling rate of the CCD/CMOS device.

Two-step phase-shift interferometry with known but arbitrary reference waves: a graphical interpretation

There are many applications in biology and metrology where it is important to be able to measure both the amplitude and phase of an optical wave field. There are several different techniques for making this type of measurement, including digital holography and phase retrieval methods. In this paper we propose an analytical generalization of this two-step phase-shifting algorithm. We investigate how to reconstruct the object signal if both reference waves are different in phase and amplitude. The resulting equations produce two different solutions and hence an ambiguity remains as to the correct solution. Because of the complexity of the generalized analytical expressions we propose a graphical-vectorial method for solution of this ambiguity problem. Combining our graphical method with a constraint on the amplitude of the object field we can unambiguously determine the correct result. The results of the simulation are presented and discussed.

Iterative reconstruction of digital holograms from three intensity measurements

Abstract. A practical technique is presented based on DH, for the reconstruction of a wavefront from three recorded intensity images. Combining the off-axis Fourier spatial filtering (OFSF) technique with iterative phase retrieval algorithms, it is shown how the twin image can be eliminated. The proposed method overcomes system geometry constraints and improves both the flexibility and resolution associated with OFSF-based DH. It also overcomes the cost problem associated with phase-shifting interferometry-based DH. In order to demonstrate the performance of the proposed DH method, both simulation and experiment results for objects having smooth and rough surfaces are presented.