Module laplacian 

Laplacian-based phase unwrapping

Uses the Laplacian operator to unwrap phase without path dependence. The wrapped phase Laplacian equals the true Laplacian, so we can recover the true phase by solving a Poisson equation.

Reference: Schofield, M.A., Zhu, Y. (2003). “Fast phase unwrapping algorithm for interferometric applications.” Optics Letters, 28(14):1194-1196. https://doi.org/10.1364/OL.28.001194

Reference implementation: https://github.com/kamesy/QSM.jl

Functions

laplacian_unwrap

Laplacian phase unwrapping

solve_poisson_fft 🔒

Solve Poisson equation using FFT (periodic boundary conditions)

wrap 🔒

Wrap angle to [-π, π]

wrapped_laplacian_periodic 🔒

Compute wrapped Laplacian of phase with periodic boundary conditions