Open AccessPacked domain Rayleigh-Sommerfeld wavefield propagation for large targets
Author(s) -
Andreas Wuttig,
Mario Kanka,
H. J. Kreuzer,
Rainer Riesenberg
Publication year - 2010
Publication title -
optics express
Language(s) - Uncategorized
Resource type - Journals
SCImago Journal Rank - 1.394
H-Index - 271
ISSN - 1094-4087
DOI - 10.1364/oe.18.027036
Subject(s) - fourier transform , convolution (computer science) , optics , fast fourier transform , kernel (algebra) , pointwise , scalar (mathematics) , frequency domain , algorithm , discrete time fourier transform , computer science , fourier analysis , physics , mathematical analysis , mathematics , short time fourier transform , geometry , artificial intelligence , combinatorics , artificial neural network
For applications in the domain of digital holographic microscopy, we present a fast algorithm to propagate scalar wave fields from a small source area to an extended, parallel target area of coarser sampling pitch, using the first Rayleigh-Sommerfeld diffraction formula. Our algorithm can take full advantage of the fast Fourier transform by decomposing the convolution kernel of the propagation into several convolution kernel patches. Using partial overlapping of the patches together with a soft blending function, the Fourier spectrum of these patches can be reduced to a low number of significant components, which can be stored in a compact sparse array structure. This allows for rapid evaluation of the partial convolution results by skipping over negligible components through the Fourier domain pointwise multiplication and direct mapping of the remaining multiplication results into a Fourier domain representation of the coarsly sampled target patch. The algorithm has been verified experimentally at a numerical aperture of 0.62, not showing any significant resolution limitations.