2006
DOI: 10.1364/ao.45.001102
|View full text |Cite
|
Sign up to set email alerts
|

Fast-Fourier-transform based numerical integration method for the Rayleigh-Sommerfeld diffraction formula

Abstract: The numerical calculation of the Rayleigh-Sommerfeld diffraction integral is investigated. The implementation of a fast-Fourier-transform (FFT) based direct integration (FFT-DI) method is presented, and Simpson's rule is used to improve the calculation accuracy. The sampling interval, the size of the computation window, and their influence on numerical accuracy and on computational complexity are discussed for the FFT-DI and the FFT-based angular spectrum (FFT-AS) methods. The performance of the FFT-DI method … Show more

Help me understand this report

Search citation statements

Order By: Relevance

Paper Sections

Select...
1
1
1
1

Citation Types

1
158
0
1

Year Published

2006
2006
2024
2024

Publication Types

Select...
7
2

Relationship

0
9

Authors

Journals

citations
Cited by 272 publications
(160 citation statements)
references
References 18 publications
(28 reference statements)
1
158
0
1
Order By: Relevance
“…The Nyquist sampling criterion states that if f max is the maximum (spatial) frequency with which the field changes in the sampling space, the aliasing is avoided when the sampling step d in the sampling window satisfies the condition 2d f À1 max . However, since the aperture is finite in its extent, the Fourier transform of the field at the aperture has infinite extent in the spatial-frequency domain [33]. Therefore, the Nyquist criterion for sampling the U 0 x 0 ; y 0 ð Þ, given in Eq.…”
Section: Samplingmentioning
confidence: 99%
See 1 more Smart Citation
“…The Nyquist sampling criterion states that if f max is the maximum (spatial) frequency with which the field changes in the sampling space, the aliasing is avoided when the sampling step d in the sampling window satisfies the condition 2d f À1 max . However, since the aperture is finite in its extent, the Fourier transform of the field at the aperture has infinite extent in the spatial-frequency domain [33]. Therefore, the Nyquist criterion for sampling the U 0 x 0 ; y 0 ð Þ, given in Eq.…”
Section: Samplingmentioning
confidence: 99%
“…The Rayleigh-Sommerfeld integral has an analytical solution for the former and the analytical solution for the latter can be derived using Fresnel approximation. When a circular aperture with radius a is illuminated with a plane wave, the diffracted field on the z axis given by the Rayleigh-Sommerfeld diffraction integral is [33].…”
Section: Appendix Amentioning
confidence: 99%
“…(4) but we can use also a numerical approach. We have used a fast Fourier Transform based direct-integration method, [18], which uses the Rayleigh-Sommerfeld approach to determine the intensity distribution at the observation plane placed at a distance z 2 from the grating. Using this numerical approach, the optical intensity after the mask is shown, Fig.…”
Section: Numerical Simulationsmentioning
confidence: 99%
“…Remind that these sort of double size calculations are typical for accurate convolutional techniques (e.g. [7]). Second, the extra area of images appeared due to zero padding is essential for the recursive inverse performed in the frequency domain [1].…”
Section: Observation Modelsmentioning
confidence: 99%