Unfolding wrapped phase

Phase unwrapping is the final step in phase extraction methods, which consists of recovering the correct phase from the wrapped phase by removing 2π discontinuities. The difference between the correct phase and the wrapped phase is the phase wrapping map. A new method for phase unwrapping is presented by identifying the phase wrapping map as a sequence of binary valued intermediate wrapping maps and iteratively removing them producing the correct phase by phase-wrapped unfolding. A path-following algorithm is presented to exemplify the phase wrapped unfolding method. © The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI. [DOI: 10.1117/1.


Introduction
Phase extraction methods are used in different coherent signal processing areas, such as optical and microwave interferometry, 1 magnetic resonance imaging, 2 synthetic aperture radar, 3 and Fourier transform optical coherent tomography, 4 among others.In these areas, the signal modulates a harmonic function and phase extraction algorithms are used to recover the signal.The recovered signal within the principal branch of the harmonic function exhibits a series of discontinuities, which need to be removed to obtain the real signal.This process, known as phase unwrapping, is the final phase extraction step.We collect these discontinuities in a wrapping map, normally overlooked based on the principle that the proposed methods aim to recover the real phase map knowing its expected continuity properties.Since phase unwrapping was first showed in optical interferometry phase extraction, 5,6 there have been extensive studies to recover the correct phase map.Because of the presence of noise, signal discontinuities, and undersampling, phase unwrapping is considered to be one of the most difficult problems in mathematics and engineering. 7][10] Independent of the approach, a definite solution for phase unwrapping remains an open problem.
On this note, an alternative solution to the phase unwrapping problem focusing on the properties of the phase wrapping map instead of the real phase signal properties is presented.

Phase Unfolding
Phase extraction methods retrieve a wrapped phase from a harmonic function confined in the main branch of the function, which is normally a 2π band in either [0;2π) or [−π; π).Because phase is a relative quantity, we can focus only on the [0;2π) band without loss of generality.In order to obtain an insight on the phase unwrapping problem instead of looking at the problem as the final step in a phase extraction methods as is customary, let us look at the problem from the solution point of view.In other words, starting from the correct phase map U [Fig. 1 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 1 ; 3 2 6 ; 5 5 3 The aim of the phase unwrapping methods aforementioned is to add a correct multiple of 2π to the wrapped phase such that a continuous unwrapped phase map is reconstructed.These methods focus on the properties of the expected phase map U and the challenge is to restore a function free of the 2π jumps.
As mentioned, instead of focusing in the properties of U, here we investigate the properties of the phase wrapping map M. To be more precise, focus on an alternative construction of M. To do such and bearing in mind that if the phase map U is sampled with scaled harmonic functions with a multiple of 2π periodicity, in principle, we can obtain the correct solution after an adequate phase unwrapping procedure.Considering using powers of two multiples of 2π (2 n π) is of particular interest.By using this approach, an insight on the phase wrapping map is obtained and it is found that it can be properly decomposed by simpler phase wrapping maps.
The first task is to obtain the minimum power of two multiples of 2π, which contains U.This is, look for N, which is the minimum integer such that E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 2 ; 3 2 6 ; 2 9 1 jUj ≤ 2 N π: (2) After determining this number, consider wrapping planes in subpowers of two multiples of 2π∶2 N−1 π, 2 N−2 π; : : : ; 4π, 2π (Fig. 2) in such a way that iterative wrapping maps M 1 ; M 2 ; : : :

and 3(h)] are obtained by wrapping the immediate previous wrapped function (
The wrapping map down to the main branch (2π) is therefore the sum of the binary intermediate wrapping maps.
E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 4 ; 3 2 6 ; 3 2 4 Using this approach instead of finding the correct multiple of 2π for phase unwrapping, the solution is reduced to reconstruct each step backward.In other words, instead of unwrapping the phase, we iteratively unfold the phase.
E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 5 ; 3 2 6 ; 2 3 9 where φ n ¼ 2 N−n π.Unfolding the phase consists on finding and removing iteratively M n .The benefit is that the partial wrapping map M n consists of only zeroes and ones.
In spite of the apparent advantage using this binary approach, the unfolding problem is as complex as the unwrapping approach.This complexity is made evident by manipulating an intermediate wrapped function U w;n placed on top of itself, thus mimicking one unfolding step.Using this construction, it is found that each point has two possible solutions (Fig. 4): Equation ( 6) is the correct solution [Fig.4(b) solid line], whereas Eq. ( 7) is a false solution, which increases the amount of wrapping [Fig.4(b) dotted line].Nevertheless, the correct solution is found among only two possible solutions and phase unwrapping has been reduced to find the correct one among these two.
3 Path-Following Phase Unfolding Algorithm Different approaches for unfolding the phase have been explored.Here, a straightforward path following method is presented.In order to implement the basic path following phase unfolding process, notice that if the intermediate phase wrapped function U n is shifted by φ n ∕2, the following folding step produces the same wrapped function U nþ1 .This property is used to simplify the unfolding process.Thus, the path following phase unfolding process can be performed by splitting the task down into the following steps: This is a binary function, which by construction contains the position of the unknown phase jumps of M n .It also contains intercalated the phase jumps of O n 3. Initialize an unfolding variable S w to zero. 4. Starting from the first sample of the left, calculate the difference in the function resulting from Eq. ( 8) between the current sample (j) and its direct adjacent right-hand neighbor (j þ 1).Only pay attention to nonzero differences.
a.For the first nonzero value, which can be either positive or negative, assign 1 or −1 to an unfolding sign variable sgn w , correspondingly.
5. For the nonzero difference points identified in step 4, if the difference between adjacent values of U n added to S w is larger than a. if S w is zero, assign sgn w ϕ n to it.b. if S w is nonzero, assign zero to it.
This step selects between the two possible paths M n or O n .6.In parallel to steps 4 and 5, starting from the first sample on the left, add S w to U w;n (U w;n−1 ¼ U w;n þ S w ). 7. Continue steps 4, 5, and 6 to the right-hand samples until all the samples have been processed.This concludes an unfolding step and U w;n−1 has been obtained.8. Double the value of φ n (φ n−1 ¼ 2φ n ) and repeat steps 1 through 7 using the obtained U w;n−1 .9. Continue until the difference in Eq. ( 8) is constant and/ or smaller than a given threshold (φ∕4 for instance) and exit in step 2.
Using this procedure after N − 1 iterations, in which the final value of φ is 2 N π, the correct phase map U is recovered.The algorithm performs the unfolding procedure by locating phase jumps on a binary function M nþ1 governed by Eq. ( 8).Therefore, in contrast to regular phase unwrapping techniques instead of searching for discontinuities on a stepwise continuous function U n at each point to execute the unwrapping decision like ice walking the unfolding decision is made only on abrupt changes in M nþ1 in a cliff hopping manner.Phase unfolding produces N − 1 times more operations than phase unwrapping methods but has the advantage that the number of discontinuity correcting decisions is reduced to the phase jumps clearly identified on M nþ1 [Eq.( 8)], which decrease on each iteration.Phase unfolding also has the same restrictions than conventional phase unwrapping: is limited by signal noise, under sampling beyond Shannon-Nyquist limit and discontinuities in the wrapped signal.Furthermore, the first iteration in the unfolding algorithm requires an additional downfolding the wrapped phase to half of the main branch [Fig.3(h)], thus increasing the sampling limitations to fulfill the Shannon-Nyquist limit.On the other hand, errors due to false wrapping do not propagate as it does in the path following phase unwrapping because of the local phase jump detection in Eq. (7).
As an example of the phase unfolding algorithm, a unitary circle function closely related to third order spherical aberration with 15π amplitude is analyzed (Fig. 6). Figure 6(a) shows the expected result of a phase extraction method confined in the [0;2π) band.In this case, N is 4 (2 4 π > 15π > 2 3 π), therefore unfolding phase unwrapping is achieved in three steps as shown.On each consecutive step [Figs.6(b)-6(d)], the phase unfolds as a folding cup opening until the correct phase map U is recovered in Fig. 6(d).
For the simple path following the method demonstrated, phase unfolding presents some advantages with respect to phase unwrapping, namely the correct identification of phase jumps, a decrease number of discontinuity correcting decisions and splitting the problem to binary maps [Eq.( 5)].It also presents disadvantages such as tighter limitation in the sampling period and the number of computational steps.Both approaches are similarly sensitive to signal noise or discontinuities in the signal.Comparable computational time has been observed using either approach.Further research is required to incorporate simultaneously the properties of both the real phase map U and the phase wrapping map M in phase unwrapping or unfolding algorithms.Additional studies also need to be conducted to fully determine the effect of the wrapping map decomposition [Eqs.( 4) and ( 5)] on either path following or minimum norm unwrapping methods.

Conclusions
In conclusion, an alternative method for phase unwrapping, named phase unfolding, has been presented by exploring the properties of the phase wrapping map decomposed in a series of binary intermediate wrapping maps.An effective iterative algorithm for path following phase unfolding was also presented, identifying some advantages and disadvantages with respect to standard phase unwrapping.
(a)], the wrapped phase U w is obtained [Fig.1(b)] by adding multiples of 2π [Fig.1(c)] to contain U w in the [0;2π) band.The form of the integer-valued correcting field [Fig.1(c)], or phase wrapping map M, is dependent on the properties of U.
by the next subpower of two multiples of 2π.By doing so the n'th wrapping map has binary values: zero or 2 N−nþ1 π.Notice that by construction, the position, but not necessarily the sign, of the phase jumps for each wrapping map is contained in the next one.Thus, the following intermediate wrapped functions are constructed [Figs.3(a), 3(c), 3(e), and 3(g)]: *Address all correspondence to: Carlos Gerardo Treviño-Palacios, E-mail: carlost@inaoep.mxOptical Engineering 110503-1 November 2015 • Vol.54(11) OE Letters Downloaded From: https://www.spiedigitallibrary.org/journals/Optical-Engineering on 12/7/2018 Terms of Use: https://www.spiedigitallibrary.org/terms-of-use E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 3 ; 6 3 ; 1 9 2

Fig. 2
Fig. 2 Alternative wrapping planes separated in a 2 N−n π sequence for the function on Fig. 1.In this case, N ¼ 4 and planes for M 1 ; M 2 , and M 3 are considered.

1 . 2 .
Start with an intermediate wrapped function U w;n ¼ U w;n−1 þ M n φ n and roll the phase within its [0; φ n ] domain by φ n ∕2.This is done by adding φ n ∕2 to values lower than half the step amplitude (φ n ) and subtracting φ n ∕2 to the upper half values [Figs.5(a) and 5 (b)].The resulting intermediate wrapped phase is U w;n−1 þ O n φ n , where O n is the correspondent phase wrapping map for the φ n ∕2 shifted U w;n function.a. for the first iteration, use φ N−1 ¼ 2π and U w;n ¼ U w .Subtract the previous two functions [Fig.5(c)]: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 8 ; 6 3 ; 1 6 3

Fig. 4
Fig. 4 Graphical principle of the phase unfold method.(a) Original U w;n intermediate wrapped phase and (b) U w;n duplicated and placed on top of itself showing the correct solution U w;n−1 (solid line) and false solution in Eq. (7) (dotted line).

Fig. 5
Fig. 5 (a) Intermediate wrapped function U w;n ¼ U w;n−1 − M n φ n , the arrows show the required displacements to roll the phase by φ n ∕2, as described in the text, (b) complementary intermediate wrapped function after a φ n ∕2 roll, U w;n − O n φ n and (c) difference between M n φ n and O n φ n signifying M nþ1 φ n .

Fig. 6
Fig. 6 Phase unfolding algorithm operating on the unitary circle function −10π {6ðx 2 þ y 2 Þ 2 − 6ðx 2 þ y 2 Þ}.Only half of the function is shown to better appreciate the effect of the phase unfolding algorithm.(a) Original wrapped function U w ¼ U w;3 in the [0;2π) domain, (b) first unfolding iteration U w;2 in the [0;4π) domain, (c) second iteration U w;1 in the [0; 8π) domain, and (d) third and final iteration recovering the correct phase map U w;0 ¼ U with a 15π amplitude.