Statistics for Spatial Data.

[1] The Shuttle Radar Topography Mission produced the most complete, highest-resolution digital elevation model of the Earth. The project was a joint endeavor of NASA, the National Geospatial-Intelligence Agency, and the German and Italian Space Agencies and flew in February 2000. It used dual radar antennas to acquire interferometric radar data, processed to digital topographic data at 1 arc sec resolution. Details of the development, flight operations, data processing, and products are provided for users of this revolutionary data set.

/pdf/the-shuttle-radar-topography-mission-4x49lt5xnj.pdf

The Shuttle Radar Topography Mission

Temporal and geometrical decorrelation often prevents SAR interferometry from being an operational tool for surface deformation monitoring and topographic profile reconstruction. Moreover, atmospheric disturbances can strongly compromise the accuracy of the results. The authors present a complete procedure for the identification and exploitation of stable natural reflectors or permanent scatterers (PSs) starting from long temporal series of interferometric SAR images. When, as it often happens, the dimension of the PS is smaller than the resolution cell, the coherence is good even for interferograms with baselines larger than the decorrelation one, and all the available images of the ESA ERS data set can be successfully exploited. On these pixels, submeter DEM accuracy and millimetric terrain motion detection can be achieved, since atmospheric phase screen (APS) contributions can be estimated and removed. Examples are then shown of small motion measurements, DEM refinement, and APS estimation and removal in the case of a sliding area in Ancona, Italy. ERS data have been used.

Permanent scatterers in SAR interferometry

We present a new differential synthetic aperture radar (SAR) interferometry algorithm for monitoring the temporal evolution of surface deformations. The presented technique is based on an appropriate combination of differential interferograms produced by data pairs characterized by a small orbital separation (baseline) in order to limit the spatial decorrelation phenomena. The application of the singular value decomposition method allows us to easily "link" independent SAR acquisition datasets, separated by large baselines, thus increasing the observation temporal sampling rate. The availability of both spatial and temporal information in the processed data is used to identify and filter out atmospheric phase artifacts. We present results obtained on the data acquired from 1992 to 2000 by the European Remote Sensing satellites and relative to the Campi Flegrei caldera and to the city of Naples, Italy, that demonstrate the capability of the proposed approach to follow the dynamics of the detected deformations.

/pdf/a-new-algorithm-for-surface-deformation-monitoring-based-on-ot8ep5chnm.pdf

A new algorithm for surface deformation monitoring based on small baseline differential SAR interferograms

On February 22, 2000 Space Shuttle Endeavour landed at Kennedy Space Center, completing the highly successful 11-day flight of the Shuttle Radar Topography Mission (SRTM). Onboard were over 300 high-density tapes containing data for the highest resolution, most complete digital topographic map of Earth ever made. SRTM is a cooperative project between NASA and the National Imagery and Mapping Agency (NIMA) of the U.S. Department of Defense. The mission was designed to use a single-pass radar interferometer to produce a digital elevation model (DEM) of the Earth's land surface between about 60 deg north and 56 deg south latitude. When completed, the DEM will have 30 m pixel spacing and about 15 m vertical accuracy. Two orthorectified image mosaics (one from the ascending passes with illumination from the southeast and one from descending passes with illumination from the southwest) will also be produced.

Interferometric synthetic aperture radar observations provide a means for obtaining high-resolution digital topographic maps from measurements of amplitude and phase of two complex radar images. The phase of the radar echoes may only be measured modulo 2 pi; however, the whole phase at each point in the image is needed to obtain elevations. An approach to 'unwrapping' the 2 pi ambiguities in the two-dimensional data set is presented. It is found that noise and geometrical radar layover corrupt measurements locally, and these local errors can propagate to form global phase errors that affect the entire image. It is shown that the local errors, or residues, can be readily identified and avoided in the global phase estimation. A rectified digital topographic map derived from the unwrapped phase values is presented.

http://paristech.institutoptique.fr/site.php?id=704&fileid=7037

Satellite radar interferometry: Two-dimensional phase unwrapping

A radar interferometric technique for topographic mapping of surfaces, implemented utilizing a single synthetic aperture radar (SAR) system in a nearly repeating orbit, is discussed. The authors characterize the various sources contributing to the echo correlation statistics, and isolate the term which most closely describes surficial change. They then examine the application of this approach to topographic mapping of vegetated surfaces which may be expected to possess varying backscatter over time. It is found that there is decorrelation increasing with time but that digital terrain model generation remains feasible. The authors present such a map of a forested area in Oregon which also includes some nearly unvegetated lava flows. Such a technique could provide a global digital terrain map. >

Decorrelation in interferometric radar echoes

A technique is described, based on synthetic aperture radar (SAR) interferometry, which uses SAR images for measuring very small (1 cm or less) surface motions with good resolution (10 m) over swaths of up to 50 km. The method was applied to a Seasat data set of an imaging site in Imperial Valley, California, where motion effects were observed that were identified with movements due to the expansion of water-absorbing clays. The technique can be used for accurate measurements of many geophysical phenomena, including swelling and buckling in fault zones, residual displacements from seismic events, and prevolcanic swelling.

Mapping small elevation changes over large areas: Differential radar interferometry

[1] We present here a new InSAR persistent scatterer (PS) method for analyzing episodic crustal deformation in non-urban environments, with application to volcanic settings. Our method for identifying PS pixels in a series of interferograms is based primarily on phase characteristics and finds low-amplitude pixels with phase stability that are not identified by the existing amplitude-based algorithm. Our method also uses the spatial correlation of the phases rather than a well-defined phase history so that we can observe temporally-variable processes, e.g., volcanic deformation. The algorithm involves removing the residual topographic component of flattened interferogram phase for each PS, then unwrapping the PS phases both spatially and temporally. Our method finds scatterers with stable phase characteristics independent of amplitudes associated with man-made objects, and is applicable to areas where conventional InSAR fails due to complete decorrelation of the majority of scatterers, yet a few stable scatterers are present.

/pdf/a-new-method-for-measuring-deformation-on-volcanoes-and-3x99a7e12s.pdf

A new method for measuring deformation on volcanoes and other natural terrains using InSAR persistent scatterers

The production of high-resolution topographic maps derived from interferometric synthetic aperture radar observations of the earth is reported. Topographic maps are typically determined from stereo-pair optical photographs. Vertical relief causes the same terrain to appear in a slightly different projection for differing look angles, and this shift in appearance is interpreted in terms of the height of the terrain. The radar interferometric approach is related to the stereo technique in that the terrain is viewed at two different angles; however, in this case, the angular separation of the antennas is extremely small, on the order of a milliradian or less, as compared to tens of degrees for the optical case. Thus, the geometrical distortion and subsequent rectification correction algorithms are much less severe in the reduction of interferometric data.

Howard A. Zebker

Papers

Satellite radar interferometry: Two-dimensional phase unwrapping

Decorrelation in interferometric radar echoes

Mapping small elevation changes over large areas: Differential radar interferometry

A new method for measuring deformation on volcanoes and other natural terrains using InSAR persistent scatterers

Topographic mapping from interferometric synthetic aperture radar observations