GMTSAR: An InSAR Processing System Based on Generic Mapping Tools

GMTSAR is an open source (GNU General Public License) InSAR processing system designed for users familiar with Generic Mapping Tools (GMT). The code is written in C and will compile on any computer where GMT and NETCDF are installed. The system has three main components: 1) a preprocessor for each satellite data type (e.g., ERS, Envisat, and ALOS) to convert the native format and orbital information into a generic format; 2) an InSAR processor to focus and align stacks of images, map topography into phase, and form the complex interferogram; 3) a postprocessor, mostly based on GMT, to filter the interferogram and construct interferometric products of phase, coherence, phase gradient, and line-of-sight displacement in both radar and geographic coordinates. GMT is used to display all the products as postscript files and kml-images for Google Earth. A set of C-shell scripts has been developed for standard 2- pass processing as well as image alignment for stacking and time series. ScanSAR processing is also possible but requires a knowledgeable user. Users are welcome to contribute to this effort.

[1]  Josef Mittermayer,et al.  TOPS Sentinel-1 and TerraSAR-X Processor Comparison based on Simulated Data , 2010 .

[2]  Thomas A. Hennig,et al.  The Shuttle Radar Topography Mission , 2001, Digital Earth Moving.

[3]  Howard A. Zebker,et al.  ScanSAR-to-Stripmap Mode Interferometry Processing Using ENVISAT/ASAR Data , 2007, IEEE Transactions on Geoscience and Remote Sensing.

[4]  D. Sandwell,et al.  Phase gradient approach to stacking interferograms , 1998 .

[5]  Ronald N. Bracewell,et al.  The Fourier Transform and Its Applications , 1966 .

[6]  R. Harger.  Synthetic aperture radar systems , 1970 .

[7]  F. Webb,et al.  Surface deformation and coherence measurements of Kilauea Volcano, Hawaii, from SIR C radar interferometry , 1996 .

[8]  Michele Manunta,et al.  Geometrical SAR image registration , 2006, IEEE Transactions on Geoscience and Remote Sensing.

[9]  Ron Kwok,et al.  Measurement of ice-sheet topography using satellite-radar interferometry , 1996 .

[10]  Ian G. Cumming,et al.  Application of efficient linear FM matched filtering algorithms to synthetic aperture radar processing , 1985 .

[11]  C. W. Chen,et al.  Network approaches to two-dimensional phase unwrapping: intractability and two new algorithms. , 2000, Journal of the Optical Society of America. A, Optics, image science, and vision.

[12]  Paul N. Swarztrauber,et al.  Vectorizing the FFTs , 1982 .

[13]  Michael Eineder,et al.  ScanSAR processing using standard high precision SAR algorithms , 1996, IEEE Trans. Geosci. Remote. Sens..

[14]  William H. Press,et al.  Numerical recipes in C , 2002 .

[15]  Robert J. Mellors,et al.  Land subsidence in the Cerro Prieto Geothermal Field, 1 Baja California, Mexico, from 1994 to 2005. An integrated analysis of DInSAR, levelingand geological data. , 2011 .

[16]  Josef Mittermayer,et al.  TOPS Imaging With TerraSAR-X: Mode Design and Performance Analysis , 2010, IEEE Transactions on Geoscience and Remote Sensing.

[17]  Richard M. Goldstein,et al.  Radar ice motion interferometry , 1997 .

[18]  David T. Sandwell,et al.  Coseismic slip model of the 2008 Wenchuan earthquake derived from joint inversion of interferometric synthetic aperture radar, GPS, and field data , 2010 .

[19]  John C. Curlander,et al.  Synthetic Aperture Radar: Systems and Signal Processing , 1991 .