Phase Wrapping Method

• Phase wrapping is the process of reducing a continuous phase measurement to a principal interval, making values separated by 2π indistinguishable.
• Techniques include phase unwrapping, fringe-order determination, and frequency shifting to reconstruct absolute phase in applications like interferometry and profilometry.
• Advanced methods leverage derivative integration, residue manipulation, and deep learning for optimized, real-time phase recovery and denoising.

Phase wrapping denotes the reduction of a continuous phase, or of a phase-difference quantity, to a principal interval such as $(-\pi,\pi]$ or $[-\pi,\pi)$, so that values differing by integer multiples of $2\pi$ become observationally equivalent. In adaptive-optics notation, the wrapped phase is

$$\Phi_w=\mathrm{Arg}(e^{i\Phi})=\tan^{-1}\!\left(\frac{\mathrm{Im}(e^{i\Phi})}{\mathrm{Re}(e^{i\Phi})}\right)=\Phi-2\pi\,\mathrm{floor}\!\left[\frac{\Phi+\pi}{2\pi}\right],$$

while in interferometry and fringe projection the same phenomenon appears after trigonometric demodulation, in direction-of-arrival estimation it appears in wrapped phase-difference observations, and in speech analysis it appears in wrapped short-time Fourier phase (Huerta et al., 19 Aug 2025, Qi et al., 2022, Chen et al., 2021, Ai et al., 2022). Phase wrapping methods therefore form a family of techniques that either reconstruct a continuous or absolute phase, reduce the number of wraps before reconstruction, or exploit the geometry of the wrapped representation directly.

1. Formal structure of wrapped phase

Across the cited literatures, the essential ambiguity is identical: a measured phase value is only known modulo $2\pi$. In speech STFT analysis, the complex spectrogram is written as

$S(t,k)=A(t,k)e^{j\phi(t,k)},$

with wrapped phase restricted to $(-\pi,\pi]$. In interferometry and fringe projection, the normalized interferogram function $F$ or the arctangent-demodulated phase yields the principal value of a cosine or sine model, again forcing discontinuities at branch boundaries. In array processing, the measured phase difference of a sensor pair is observed through

$\psi_{ij}=\mathrm{mod}(\phi_{ij}+\pi,2\pi)-\pi,$

so different physical angles can produce the same wrapped observation when the inter-sensor spacing exceeds $\lambda/2$ (Ai et al., 2022, Berejnov et al., 2024, Chen et al., 2021).

The standard reconstruction target is an absolute or continuous phase. In structured-light profilometry this is written explicitly as

$[-\pi,\pi)$0

where $[-\pi,\pi)$1 is the fringe order. In polynomial phase estimation, the same ambiguity is handled by integer offsets $[-\pi,\pi)$2 added to wrapped observations before least-squares regression. In adaptive optics, continuous-surface deformable mirrors require an unwrapped phase even when the sensor itself returns only the principal value (Qi et al., 2022, McKilliam et al., 2012, Huerta et al., 19 Aug 2025).

Several works also replace naive Euclidean phase differences by circular discrepancies. In neural speech phase prediction, the anti-wrapping function

$[-\pi,\pi)$3

measures the shortest angular distance on the circle and prevents error expansion near the $[-\pi,\pi)$4 boundary. This same circular viewpoint underlies wrapped-distance formulations in other domains, even when the optimization problem is expressed differently (Ai et al., 2022, Ai et al., 2024).

2. Fringe projection and structured-light strategies

In structured-light and fringe-projection profilometry, phase wrapping methods are commonly organized around fringe-order determination, carrier removal, and residue management. The bidirectional coding method for absolute phase measurement combines four phase-shifting patterns with two bidirectional coded patterns, so that the wrapped phase is obtained by four-step demodulation and the fringe order is determined locally from a local stair-phase code and a partition code. Its defining rule is

$[-\pi,\pi)$5

and the absolute phase is then $[-\pi,\pi)$6. The paper reports six images total, namely 4 phase-shift patterns plus 2 bidirectional stair-coded patterns, and gives RMSE values of approximately $[-\pi,\pi)$7 rad in Experiment 1, $[-\pi,\pi)$8 rad in Experiment 2, and $[-\pi,\pi)$9 rad in Experiment 3 relative to a three-frequency reference (Qi et al., 2022).

A second line of work reduces wraps before any explicit unwrapping. “Precise and Fast Phase Wraps Reduction in Fringe Projection Profilometry” constructs the complex wrapped phase $2\pi$0, estimates the carrier peak through iterative local discrete Fourier transform upsampling, and then applies a non-integer spatial-domain frequency shift. In the reported simulation, the RMSE values were $2\pi$1 rad for conventional integer shifting, $2\pi$2 rad for zero-padding, and $2\pi$3 rad for the proposed method, with corresponding times of $2\pi$4 s, $2\pi$5 s, and $2\pi$6 s (Wang et al., 2016). The related “Improved method for phase wraps reduction in profilometry” used 1D zero-padding to estimate non-integer carrier frequencies and reported $2\pi$7 s for Gdeisat’s method versus $2\pi$8 s for the proposed method; in a second experiment Herráez unwrapping times were $2\pi$9 s without carrier removal, $$\Phi_w=\mathrm{Arg}(e^{i\Phi})=\tan^{-1}\!\left(\frac{\mathrm{Im}(e^{i\Phi})}{\mathrm{Re}(e^{i\Phi})}\right)=\Phi-2\pi\,\mathrm{floor}\!\left[\frac{\Phi+\pi}{2\pi}\right],$$0 s for Gdeisat, and $$\Phi_w=\mathrm{Arg}(e^{i\Phi})=\tan^{-1}\!\left(\frac{\mathrm{Im}(e^{i\Phi})}{\mathrm{Re}(e^{i\Phi})}\right)=\Phi-2\pi\,\mathrm{floor}\!\left[\frac{\Phi+\pi}{2\pi}\right],$$1 s for the proposed method (Du et al., 2016).

Residue-based methods manipulate the wrapped map in a different way. “An algorithm to increase the residues of wrapped-phase in spatial domain” applies a spatial-domain modulation to $$\Phi_w=\mathrm{Arg}(e^{i\Phi})=\tan^{-1}\!\left(\frac{\mathrm{Im}(e^{i\Phi})}{\mathrm{Re}(e^{i\Phi})}\right)=\Phi-2\pi\,\mathrm{floor}\!\left[\frac{\Phi+\pi}{2\pi}\right],$$2, which is equivalent to a spectral shift by the Fourier shift property, in order to increase the density of residues used by branch-cut unwrappers. In the reported hand experiment, the residue count increased from $$\Phi_w=\mathrm{Arg}(e^{i\Phi})=\tan^{-1}\!\left(\frac{\mathrm{Im}(e^{i\Phi})}{\mathrm{Re}(e^{i\Phi})}\right)=\Phi-2\pi\,\mathrm{floor}\!\left[\frac{\Phi+\pi}{2\pi}\right],$$3 to $$\Phi_w=\mathrm{Arg}(e^{i\Phi})=\tan^{-1}\!\left(\frac{\mathrm{Im}(e^{i\Phi})}{\mathrm{Re}(e^{i\Phi})}\right)=\Phi-2\pi\,\mathrm{floor}\!\left[\frac{\Phi+\pi}{2\pi}\right],$$4, and runtime decreased from $$\Phi_w=\mathrm{Arg}(e^{i\Phi})=\tan^{-1}\!\left(\frac{\mathrm{Im}(e^{i\Phi})}{\mathrm{Re}(e^{i\Phi})}\right)=\Phi-2\pi\,\mathrm{floor}\!\left[\frac{\Phi+\pi}{2\pi}\right],$$5 s to $$\Phi_w=\mathrm{Arg}(e^{i\Phi})=\tan^{-1}\!\left(\frac{\mathrm{Im}(e^{i\Phi})}{\mathrm{Re}(e^{i\Phi})}\right)=\Phi-2\pi\,\mathrm{floor}\!\left[\frac{\Phi+\pi}{2\pi}\right],$$6 s relative to the earlier FFT-based procedure (Du et al., 2016). This suggests a distinction within profilometry between methods that suppress wraps and methods that deliberately reshape the wrap pattern to expose inconsistent regions more clearly.

3. Geometry-, projection-, and wrap-counting methods

Some phase wrapping methods do not treat wrapping as a nuisance to be removed immediately, but instead as a geometric structure from which the unknown quantity can be inferred. In non-uniform linear arrays, wrapped phase differences across sensor pairs form a wrapped phase-difference pattern (WPDP) that evolves along piecewise linear segments in an $$\Phi_w=\mathrm{Arg}(e^{i\Phi})=\tan^{-1}\!\left(\frac{\mathrm{Im}(e^{i\Phi})}{\mathrm{Re}(e^{i\Phi})}\right)=\Phi-2\pi\,\mathrm{floor}\!\left[\frac{\Phi+\pi}{2\pi}\right],$$7-dimensional space. The phase-difference projection (PDP) algorithm exploits the fact that these WPDP lines are parallel to the spacing vector $$\Phi_w=\mathrm{Arg}(e^{i\Phi})=\tan^{-1}\!\left(\frac{\mathrm{Im}(e^{i\Phi})}{\mathrm{Re}(e^{i\Phi})}\right)=\Phi-2\pi\,\mathrm{floor}\!\left[\frac{\Phi+\pi}{2\pi}\right],$$8 and perpendicular to the hyperplane $$\Phi_w=\mathrm{Arg}(e^{i\Phi})=\tan^{-1}\!\left(\frac{\mathrm{Im}(e^{i\Phi})}{\mathrm{Re}(e^{i\Phi})}\right)=\Phi-2\pi\,\mathrm{floor}\!\left[\frac{\Phi+\pi}{2\pi}\right],$$9. A measured wrapped phase-difference vector $2\pi$0 is projected onto that hyperplane,

$2\pi$1

matched to the nearest precomputed projection point, and then unwrapped without any angular grid search. The paper reports online multiplication count $2\pi$2, and for a 3-sensor array gives $2\pi$3 multiplications for PDP versus $2\pi$4 for MLE with a fine grid. In simulation, PDP closely matched MLE with a $2\pi$5 grid and reached the CRLB at high SNR across the tested NULA configurations (Chen et al., 2021).

An analogous strategy appears in millimetre-VLBI frequency phase transfer. Conventional FPT fails for non-integer frequency ratios because each wrap at the reference frequency produces a fractional-cycle jump at the target frequency. The phase-wrap counting (PWC) method writes the residual discontinuity as

$2\pi$6

where $2\pi$7 is the fractional part of the frequency ratio and $2\pi$8 is the reference-band wrap count, and then compensates it explicitly. In the reported Event Horizon Telescope simulation, PWC calibrated 230-GHz phases using 86-GHz phase solutions, achieved fractional peak flux recovery $2\pi$9 versus $S(t,k)=A(t,k)e^{j\phi(t,k)},$0 for conventional fringe-fitting plus self-calibration, and ran in under $S(t,k)=A(t,k)e^{j\phi(t,k)},$1 s using two logical cores and $S(t,k)=A(t,k)e^{j\phi(t,k)},$2 GB RAM (Simelane et al., 25 Jun 2026).

Both examples treat wrap indices as latent geometric or combinatorial variables rather than purely local discontinuities. A plausible implication is that when the wrapped observation has a rigid low-dimensional structure, line identity or wrap-count estimation can replace generic unwrapping.

4. Continuous recovery without explicit unwrapping

A distinct family of methods bypasses modular reconstruction by deriving differential equations directly from the interferogram. In the continuous phase unwrapping method, one begins from

$S(t,k)=A(t,k)e^{j\phi(t,k)},$3

differentiates, and eliminates the trigonometric functions to obtain

$S(t,k)=A(t,k)e^{j\phi(t,k)},$4

The phase is then reconstructed by piecewise integration with sign alternation at extrema. In the synthetic parabola example $S(t,k)=A(t,k)e^{j\phi(t,k)},$5, the method reproduced the exact continuous phase over the full interval (Berejnov et al., 2024).

“Continuous Recovery of Phase from Single Interferogram” extends this construction to straight paths in two-dimensional two-beam interferograms and makes the directional derivative explicit: $S(t,k)=A(t,k)e^{j\phi(t,k)},$6 The same paper also gives a multiple-beam thin-film variant with a factor $S(t,k)=A(t,k)e^{j\phi(t,k)},$7 in the denominator and discusses Newton-type versus Fizeau-type interferograms, root classification, and straight-path integration with sign changes only at true extrema (Berejnov et al., 8 Sep 2025).

A related wrapping-free line replaces explicit branch management by weighted least-squares integration of phase derivatives. “Robust Wrapping-free Phase Retrieval Method Based on Weighted Least-square Method” computes carrier-free quadrature signals, extracts $S(t,k)=A(t,k)e^{j\phi(t,k)},$8 and $S(t,k)=A(t,k)e^{j\phi(t,k)},$9, detects singular regions through the derivative variance correlation map, and solves a weighted Poisson-type equation. In the reported simulations, Perciante’s method gave RMSE $(-\pi,\pi]$0 rad in a shadowed case and $(-\pi,\pi]$1 rad in a shear/discontinuity case, while the proposed binary-mask weighted method reduced these to $(-\pi,\pi]$2 rad and $(-\pi,\pi]$3 rad respectively (Wang et al., 2017).

These methods are unified by a shared principle: they recover continuous phase by integrating derivative information or by solving a global field equation, rather than by piecing together locally wrapped samples.

5. Optimization, denoising, and real-time reconstruction

Adaptive optics emphasizes the computational side of phase wrapping. Complex-field wavefront sensors estimate a wrapped phase that must be unwrapped for continuous-surface deformable mirrors, and the reported real-time target is $(-\pi,\pi]$4 millisecond or faster. In a $(-\pi,\pi]$5 Kolmogorov-screen simulation with $(-\pi,\pi]$6 trials, four mature methods were benchmarked: Fast2D, Zernike Gradient, DFT, and LSPV (Huerta et al., 19 Aug 2025).

Method	RMS WFE without / with boundary	Latency (ms)
Fast2D	$(-\pi,\pi]$7 / $(-\pi,\pi]$8	$(-\pi,\pi]$9
LSPV	$F$0 / $F$1	$F$2
Zernike	$F$3 / $F$4	$F$5
DFT	$F$6 / $F$7	$F$8

None of these implementations reached the $F$9 ms CPU target. Fast2D remained near the numerical floor even with a binary circular aperture, LSPV degraded at the boundary but stayed diffraction-limited, Zernike was essentially unchanged by the circular mask, and DFT improved with masking because boundary discontinuities were reduced (Huerta et al., 19 Aug 2025).

Wrapped phase denoising introduces another optimization layer. “A variational model for wrapped phase denoising” does not unwrap at all; it denoises $\psi_{ij}=\mathrm{mod}(\phi_{ij}+\pi,2\pi)-\pi,$0 and $\psi_{ij}=\mathrm{mod}(\phi_{ij}+\pi,2\pi)-\pi,$1 jointly with isotropic total variation and a soft unit-circle penalty $\psi_{ij}=\mathrm{mod}(\phi_{ij}+\pi,2\pi)-\pi,$2, then reconstructs $\psi_{ij}=\mathrm{mod}(\phi_{ij}+\pi,2\pi)-\pi,$3. The paper reports fixed-point iteration counts and runtimes of $\psi_{ij}=\mathrm{mod}(\phi_{ij}+\pi,2\pi)-\pi,$4 versus $\psi_{ij}=\mathrm{mod}(\phi_{ij}+\pi,2\pi)-\pi,$5 iterations and $\psi_{ij}=\mathrm{mod}(\phi_{ij}+\pi,2\pi)-\pi,$6 s versus $\psi_{ij}=\mathrm{mod}(\phi_{ij}+\pi,2\pi)-\pi,$7 s in one synthetic case, and $\psi_{ij}=\mathrm{mod}(\phi_{ij}+\pi,2\pi)-\pi,$8 versus $\psi_{ij}=\mathrm{mod}(\phi_{ij}+\pi,2\pi)-\pi,$9 iterations and $\lambda/2$0 s versus $\lambda/2$1 s in another, when compared with gradient descent (May-Cen et al., 2023).

Digital wavefront sensors provide still another route. “On Phase Unwrapping via Digital Wavefront Sensors” observes that common wavefront-sensor forward models depend only on $\lambda/2$2, so they are insensitive to whether the incoming phase is wrapped. A wrapped phase can therefore be digitally propagated through a Shack–Hartmann or Fourier-type wavefront sensor, and the smooth wavefront returned by CuReD, PCuReD, or NOPE functions as an unwrapped solution. On the OCTOPUS dataset, digital PWFS with NOPE, linear start, and smoothness index $\lambda/2$3 achieved relative error $\lambda/2$4, compared with $\lambda/2$5 for digital SH-WFS at $\lambda/2$6 subapertures and $\lambda/2$7 for MATLAB unwrap (Hubmer et al., 2024).

6. Learning-based prediction and broader uses

In neural speech generation, several recent works predict wrapped phase directly instead of unwrapping it. “Neural Speech Phase Prediction based on Parallel Estimation Architecture and Anti-Wrapping Losses” uses a residual convolutional network followed by a parallel estimation architecture that outputs pseudo real and imaginary parts and computes phase with a formula equivalent to $\lambda/2$8, guaranteeing values in $\lambda/2$9. Training uses the anti-wrapping loss

$[-\pi,\pi)$00

together with group-delay and instantaneous-angular-frequency losses. On analysis-synthesis from amplitude, the reported objective scores were SNR $[-\pi,\pi)$01 dB, F0-RMSE $[-\pi,\pi)$02 cent, and RTF $[-\pi,\pi)$03, while the MOS was $[-\pi,\pi)$04, statistically indistinguishable from ground truth with $[-\pi,\pi)$05 (Ai et al., 2022). The low-latency extension retained the same anti-wrapping construction and, through causal convolutions and knowledge distillation, reduced latency to $[-\pi,\pi)$06 ms while keeping the wrapped-phase prediction paradigm (Ai et al., 2024).

In ptychography, phase wrapping is tied to frequency-domain asymptotes rather than to a conventional spatial unwrapping stage. “Phase offset method of ptychographic contrast reversal correction” shows that phase-wrap asymptotes in the frequency response of single-side-band ptychography can cause contrast reversals without requiring dynamical scattering. The proposed correction applies a rigid offset

$[-\pi,\pi)$07

to the nonzero spatial-frequency components, and the paper states that this phase offset method is generally superior to post-collection defocus, with greater reliability and generally stronger contrast across thin and thick samples (Hofer et al., 2023).

The term also appears in a physically distinct sense outside modulo-$[-\pi,\pi)$08 inversion. In galactic dynamics, “Phase Wrapping of Epicyclic Perturbations in the Wobbly Galaxy” uses phase wrapping to describe the winding of stellar epicyclic phases after a perturbing impulse. In that usage, phase wrapping means time-dependent shearing of $[-\pi,\pi)$09 and $[-\pi,\pi)$10, not principal-value folding. The paper reports vertical-gradient ranges $[-\pi,\pi)$11 km s$[-\pi,\pi)$12 kpc$[-\pi,\pi)$13 for the E1 run and $[-\pi,\pi)$14 km s$[-\pi,\pi)$15 kpc$[-\pi,\pi)$16 for E2, emphasizing that “phase wrapping” is not semantically identical across all disciplines (Vega et al., 2015).

Taken together, these works suggest that phase wrapping methods are best understood not as a single algorithm but as a technical family organized around one invariant: the observable phase is periodic, whereas the sought quantity is usually continuous, absolute, or structurally interpretable. Different fields therefore either reconstruct the hidden $[-\pi,\pi)$17-equivalence class, suppress its visible discontinuities, or exploit its geometric consequences directly.