2-D weighted least-squares phase unwrapping

• US 5,424,743 A
• Filed: 06/01/1994
• Issued: 06/13/1995
• Est. Priority Date: 06/01/1994
• Status: Expired due to Fees

• Current Assignee: United States Of America As Represented By The Department Of Energy
• Original Assignee: US Department of Energy (Government of the United States of America)
• Inventors: Ghiglia, Dennis C., Romero, Louis A.
• Primary Examiner(s): Sotomayor, John B.
• Application Number:
  US08/251,988
• Time in Patent Office: 377 Days
• Field of Search: 342/25, 342/159, 342/162, 342/191, 342/194, 342/195
• US Class Current: 342/25.00C
• CPC Class Codes: G01S 13/9023 combined with interferometric techniques

Abstract

Weighted values of interferometric signals are unwrapped by determining the least squares solution of phase unwrapping for unweighted values of the interferometric signals; and then determining the least squares solution of phase unwrapping for weighted values of the interferometric signals by preconditioned conjugate gradient methods using the unweighted solutions as preconditioning values. An output is provided that is representative of the least squares solution of phase unwrapping for weighted values of the interferometric signals.

Claims

1. A system for phase unwrapping weighted values of interferometric signals comprising:

means for determining the least squares solution of phase unwrapping for unweighted values of the interferometric signals; and

means for determining the least squares solution of phase unwrapping for weighted values of the interferometric phase signals by preconditioned conjugate gradient methods using the unweighted solutions as preconditioning values.

2. The system of claim 1 wherein said interferometric phase signals are generated by a pair of spaced synthetic aperture radars.

3. A method for phase unwrapping weighted values of interferometric signals comprising:

generating interferometric signals;

determining the least squares solution of phase unwrapping for unweighted values of the interferometric signals;

determining the least squares solution of phase unwrapping for weighted values of the interferometric signals by preconditioned conjugate gradient methods using the unweighted solutions as preconditioning values; and

generating an output representative of the least squares solution of phase unwrapping for weighted values of the interferometric signals.

4. The method of claim 3 wherein said interferometric signals are generated by SAR.

5. The method of claim 3 wherein the step of determining the least squares solution of phase unwrapping for unweighted values of the interferometric signals comprises:

determining the 2D forward direct cosine transform values ρ

_i,j, of the array of values, ρ

_i,j, where ##EQU16## where the x and y superscripts refer to differences in the i and j indices, respectively; and

performing the 2D inverse direct cosine transform of Φ

_i,j to obtain the least-squares unwrapped phase values Φ

_i,j, where ##EQU17##

6. The method of claim 5 wherein said interferometric signals are generated by SAR.

7. The method of claim 5 wherein said step of determining the least squares solution of phase unwrapping for weighted values of the interferometric signals comprises:

(a) determine the modified discrete Laplacian of the orignal wrapped phases ##EQU18## (b) set ρ

_i,j =c_i,j (c) determine the least squares solution as set forth in claim 5;

(d) perform scalar and vector updates, ##EQU19## where ##EQU20## (e) repeat (c) using values of r_k from (c) until ∥

r_k ∥

converges to an acceptable value.

8. The method of claim 7 wherein said interferometric signals are generated by SAR.