Fast numerical propagation in high-NA imaging using the resampling angular spectrum method

Abstract

Numerical propagation calculation is a fundamental research topic in optical engineering. The standard angular spectrum method (ASM) is accurate but time- and memory-consuming, especially for high-NA systems. In this work, we propose a fast and simple numerical propagation method, the resampling ASM (RS-ASM). Numerical propagation can be accelerated by combining a resampling technique with interpolation methods in the angular spectrum domain of a constrained object at the focal plane. RS-ASM has three main advantages: simple implementation, faster calculation than the standard ASM, and SNR enhancement. Here we validate RS-ASM using theory, simulation and experiment. Using the “bilinear” ASM with a proper resampling factor can result in a speed-up factor of up to 20x (for a transformation from the angular spectrum to the E field) and 4x (for a transformation from E field to the angular spectrum), together with a SNR improvement of approximately 2x. For an application example of Gerchberg-Saxton phase reconstruction, the “bilinear” RS-ASM can converge 2.6x faster than the standard ASM.