This repository contains codes for in-line holography muti-height Gerchberg-Saxton phase retrieval with automatic affine transform preprocessing.
Follow the 'main_GSwithAT.m' code-file for full guide how to use this algorithm.
Firstly, each in-line hologram is reconstructed (propagated to the object plane), then 'AutoAffineTransform.m' function is applied to extract and match features between n-th and last reconstruction. Basing on these features, affine transforms (AT) are estimated and applied to the input holograms to correct the xy shift and magnification mismatch between holograms (gif below)
After this preprocessing step, Gerchberg-Saxton (GS) multi-height algorithm is performed to retrieve object phase with minimized twin image noise comparing to single-frame angular spectrum backpropagation (fig. below)
Affine transform-based twin-image suppression for in-line lensless digital holographic microscopy (M. J. Marzejon, M. Rogalski, M. Trusiak) in Rosen, J., Alford, S., Allan, B., et al, "Roadmap on computational methods in optical imaging and holography [invited]," Appl. Phys. B 130, 166 (2024). https://doi.org/10.1007/s00340-024-08280-3
Mikołaj Rogalski,
mikolaj.rogalski.dokt@pw.edu.pl
Institute of Micromechanics and Photonics,
Warsaw University of Technology, Poland