Unwrapping ground displacement signals in satellite radar interferograms with aid of GPS data and MRF regularization

Synthetic aperture radar (SAR) images acquired by radar satellites at different times can be combined into interferograms that reveal information about the change in range from ground to satellite, wrapped into phase fringes corresponding to half the radar wavelength. We describe a methodology that uses Markov random field (MRF) regularization and simulated annealing optimization to unwrap such differential interferograms. Often, repeated Global Positioning System (GPS) geodetic measurements are available in an area covered by interferograms. Here, such repeated GPS observations are used to provide a complementary measurement of the unwrapped change in range at sparse locations. The process of unwrapping interferograms can be initialized and guided with such sparsely located "correct" values. Both interferograms and GPS observations may include several error factors, which are reduced before combining the two observations. GPS-measured range change is used to eliminate residual orbital error. In the unwrapping procedure, a vectorized lowpass filter is used to gain temporarily increased smoother variation of the phase. For the purposes of initializing the unwrapping process, virtual unwrapped interferograms are created by ordinary kriging of GPS-measured range change. The initial interferograms are then optimized further by using MRF regularization that incorporate the assumption of a smoothly varying displacement field and the relationship of the unwrapped images to the GPS observations. A simulated annealing optimization algorithm is used to find an optimal solution of the MRF regularization. The smoothed unwrapped interferograms are then used to construct unwrapped versions of the unfiltered input interferogram. Several additional image analyses methods are used in the optimization process to make the unwrapping more efficient and faster. The unwrapping technique is applied to unwrap interferograms from the Reykjanes Peninsula, in southwest Iceland.

[1]  Didier Massonnet,et al.  Reduction of the need for phase unwrapping in radar interferometry , 1996, IEEE Trans. Geosci. Remote. Sens..

[2]  J. Besag Spatial Interaction and the Statistical Analysis of Lattice Systems , 1974 .

[3]  A. Leick GPS satellite surveying , 1990 .

[4]  Freysteinn Sigmundsson,et al.  Crustal deformation at the oblique spreading Reykjanes Peninsula, SW Iceland: GPS measurements from 1993 to 1998 , 2001 .

[5]  Sheldon M. Ross,et al.  Introduction to Probability and Statistics for Engineers and Scientists , 1987 .

[6]  Freysteinn Sigmundsson,et al.  Crustal Deformation from 1992 to 1995 at the Mid-Atlantic Ridge, Southwest Iceland, Mapped by Satellite Radar Interferometry , 1997, Science.

[7]  J. Carstensen,et al.  Three‐dimensional surface motion maps estimated from combined interferometric synthetic aperture radar and GPS data , 2002 .

[8]  Freysteinn Sigmundsson,et al.  Three-dimensional surface motion maps estimated from combined insar and gps data , 2002 .

[9]  Stan Z. Li,et al.  Markov Random Field Modeling in Image Analysis , 2001, Computer Science Workbench.

[10]  K. Feigl,et al.  Radar interferometry and its application to changes in the Earth's surface , 1998 .

[11]  Allan Aasbjerg Nielsen,et al.  Geostatistik og analyse af spatielle data , 2004 .

[12]  Fuk K. Li,et al.  Synthetic aperture radar interferometry , 2000, Proceedings of the IEEE.

[13]  R. Hanssen Radar Interferometry: Data Interpretation and Error Analysis , 2001 .

[14]  Donald Geman,et al.  Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images , 1984, IEEE Transactions on Pattern Analysis and Machine Intelligence.