Synthetic Aperture Radar Interferometry: Principles and Algorithms
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Interferometric synthetic aperture radar (INSAR) is a new way of performing topography mapping. Among the factors critical to mapping accuracy is the registration of the complex SAR images from repeated orbits. A new algorithm for registering interferometric SAR images is presented. A new figure of merit, the average fluctuation function of the phase difference image, is proposed to evaluate the fringe pattern quality. The process of adjusting the registration parameters according to the fringe pattern quality is optimized through a downhill simplex minimization algorithm. The results of applying the proposed algorithm to register two pairs of Seasat SAR images with a short baseline (75 m) and a long baseline (500 m) are shown. It is found that the average fluctuation function is a very stable measure of fringe pattern quality allowing very accurate registration.
The first topographic and ice-motion maps of the northwestern flank of Hielo Patagonico Norte (HPN, northern Patagonia Icefield), in Chile, were produced using satellite synthetic-aperture interferometric radar data acquired by NASA's Spaceborne Imaging Radar C instrument in October 1994. The topographic map has a IO m vertical precision with a 30 m horizontal spacing, which should be sufficient to serve as a reference for monitoring future mass changes of the icefield. The ice-motion map is accurate to within 4 mm/ d (or 1/ ma). The radar-derived surface topography and ice velocity are used to estimate the ice discharge from the accumulation area of four outlet glaciers, and the calving flux and mass balance of Glaciar San Rafael. The results demonstrate the use of SAR interferometry for monitoring glaciological parameters on a spatial and temporal scale unattainable by any other means.
Explore the source record for details and available documents.
We present interferometric observations of ice-sheet motion in western Greenland based on pairs of ERS-1 synthetic aperture radar (SAR) images. These observations provide the first detailed regional view of ice motion associated with dynamically supported topography near the margin of an ice sheet. The interferograms of this area are much more complicated than other interferograms of ice sheets presented to date. We devote the largest part of this paper to explaining the source of the complexity in these interferograms. A synthetic interferogram based on a simple model helps to illustrate the effects of different components of the ice velocity field in interferometric data and suggests a method for estimating the large-scale ice velocity field from such data.
Methods are presented for using Synthetic Aperture Radar (SAR) interferometry data to estimate surface topography. An expression is given to relate the elevation of a ground point to the phase difference of SAR images received from two spatially separated antennas. An iterative algorithm which solves for the position and elevation of each point in the image simultaneously is developed. One of the critical issues that determines the accuracy of the terrain mapping is the phase unwrapping. An approach to the problem by fringe line detection is proposed. The algorithms are tested with two Seasat SAR images of terrain near Yellowstone National Park. The resultant elevation map is compared with a USGS terrain elevation model. The error of the SAR elevation with respect to the digital terrain map is about 8.2 percent of the total terrain variation.
This paper provides an overview of some of the recent developments in active microwave systems at JPL. The first section summarizes the engineering and scientific results that were obtained from the two recent flights of the Spaceborne Imaging Radar-C/X-band Synthetic Aperture Radar (SIR-C/X-SAR) on the Shuttle Endeavor. Specific examples of the applications of these data to various scientific disciplines will be discussed. The second section of the paper describes the recent progress in the development of an interferometric SAR system for high-resolution digital topography mapping. Examples of measurements from this technique over areas with earthquakes will be shown. The final section of the paper addresses development of meteorological radars for rain and cloud measurements. The design of an airborne rain mapping radar that was developed to support the Tropical Rain Measuring Mission is summarized. The technical design challenges for this high-frequency radar system will be described.
Interferometric Synthetic Aperture Radar (IFSAR) is a technique that allows an automated way to carry out terrain mapping. IFSAR is carried out by first generating a SAR image pair from two antennas that are spatially separated. The phase difference between the SAR image pair is proportional to the topography. After registering the SAR images, the difference in phase in each pixel is extracted to generate an interferogram. Since the phase can only be measured within 2pi radians, phase unwrapping is carried out to extract the absolute phase for each pixel that will be proportional to the local height. While IFSAR algorithm is typically applied to measurement data, it is useful to develop an IFSAR simulator to develop a better understanding of the IFSAR technique. The IFSAR simulator can be used in choosing system parameters, experimenting with processing procedures and mission planning. In this paper we will present an IFSAR simulation methodology to simulate the interferogram based on the shooting and bouncing ray (SBR) technique. SBR is a standard ray-tracing technique used to simulate scattering from large, complex targets. SBR is carried out by shooting rays at the target or scene. At the exit point of each ray, a ray-tube integration is done to find its contribution to the total field. A fast algorithm has been developed for the SBR for simulating SAR images of complex targets. In the IFSAR simulation, we build upon the fast SAR simulation technique. Given the antenna pair configuration, radar system parameters and the geometrical description of the scene, we first simulate two SAR images from each antenna. After post processing the two SAR images, we generate an interferogram. Phase unwrapping is then performed on the interferogram to arrive at the desired terrain map. We will present results from the SBR-based IFSAR simulator. The results will include terrain map reconstruction of urban environments. The reconstruction will be compared to the ground truth to examine the fidelity of the simulation. We will also investigate the effect of multi-bounce scattering in urban environments on phase unwrapping and reconstruction.
A large (M&subw;=7.0) earthquake on May 27, 1995 completely destroyed the town of Neftegorsk in the northern part of Sakhalin Island and caused more than 2000 human deaths. The shallow, right-lateral, strick-slip earthquake resulted in extensive surface ruptures and up to 7 m of horizontal displacement as reported by field workers. The sourthern part of the mainshock epicenter zone was imaged by the JERS-1 SAR (synthetic aperature radar) one month (April 28) before and two weeks after (June 11) the mainshock. Despite drastically changed surface conditions in the 44 days between the two images, due primarily to spring thaw, we obtained reasonably good interferometric correlation with the L-band (24 cm) SAR pair. The interoferogram records the distribution of deformation reflecting displacement during both the mainshock and aftershocks. The ability to map the deformation pattern can aid the assessment and mitigation of damage.
A report discusses Ka-band (35-GHz) radar for mapping the surface topography of glaciers and ice sheets at high spatial resolution and high vertical accuracy, independent of cloud cover, with a swath-width of 70 km. The system is a single- pass, single-platform interferometric synthetic aperture radar (InSAR) with an 8-mm wavelength, which minimizes snow penetration while remaining relatively impervious to atmospheric attenuation. As exhibited by the lower frequency SRTM (Shuttle Radar Topography Mission) AirSAR and GeoSAR systems, an InSAR measures topography using two antennas separated by a baseline in the cross-track direction, to view the same region on the ground. The interferometric combination of data received allows the system to resolve the pathlength difference from the illuminated area to the antennas to a fraction of a wavelength. From the interferometric phase, the height of the target area can be estimated. This means an InSAR system is capable of providing not only the position of each image point in along-track and slant range as with a traditional SAR but also the height of that point through interferometry. Although the evolution of InSAR to a millimeter-wave center frequency maximizes the interferometric accuracy from a given baseline length, the high frequency also creates a fundamental problem of swath coverage versus signal-to-noise ratio. While the length of SAR antennas is typically fixed by mass and stowage or deployment constraints, the width is constrained by the desired illuminated swath width. As the across-track beam width which sets the swath size is proportional to the wavelength, a fixed swath size equates to a smaller antenna as the frequency is increased. This loss of antenna size reduces the two-way antenna gain to the second power, drastically reducing the signal-to-noise ratio of the SAR system. This fundamental constraint of high-frequency SAR systems is addressed by applying digital beam-forming (DBF) techniques to synthesize multiple simultaneous receive beams in elevation while maintaining a broad transmit illumination. Through this technique, a high antenna gain on receive is preserved, thereby reducing the required transmit power and thus enabling high-frequency SARs and high-precision InSAR from a single spacecraft.
The field of synthetic aperture radar changed dramatically over the past decade with the operational introduction of advance radar techniques such as polarimetry and interferometry. Radar polarimetry became an operational research tool with the introduction of the NASA/JPL AIRSAR system in the early 1980's, and reached a climax with the two SIR-C/X-SAR flights on board the space shuttle Endeavour in April and October 1994. Radar interferometry received a tremendous boost when the airborne TOPSAR system was introduced in 1991 by NASA/JPL, and further when data from the European Space Agency ERS-1 radar satellite became routinely available in 1991. Several airborne interferometric SAR systems are either currently operational, or are about to be introduced. Radar interferometry is a technique that allows one to map the topography of an area automatically under all weather conditions, day or night. The real power of radar interferometry is that the images and digital elevation models are automatically geometrically resampled, and could be imported into GIS systems directly after suitable reformatting. When combined with polarimetry, a technique that uses polarization diversity to gather more information about the geophysical properties of the terrain, a very rich multi-layer data set is available to the remote sensing scientist. This talk will discuss the principles of radar interferometry and polarimetry with specific application to the automatic categorization of land cover. Examples will include images acquired with the NASA/JPL AIRSAR/TOPSAR system in Australia and elsewhere.
Under suitable conditions the bottom topography of shallow seas is visible in remote sensing radar imagery. Two experiments were performed to establish which remote sensing technique or combination yields optimal imaging of bottom topography and which hydro-meteorological conditions are favorable. A further goal is to gain experience with these techniques. Two experiments were performed over an area in the North Sea near the measuring platform Meetpost Noordwijk (MPN). The bottom topography in the test area is dominated by sand waves. The crests of the sand waves are perpendicular to the coast line and the dominating (tidal-)current direction. A 4x4 sq km wide section of the test area was studied in more detail. The first experiment was undertaken on 16 Aug. 1989. During the experiment the following remote sensing instruments were used: Landsat-Thematic Mapper, and NASA/JPL Airborne Imaging Radar (AIR). The hydro-meteorological conditions; current, wind, wave, and air and water temperature were monitored by MPN, a ship of Rijkswaterstaat (the OCTANS), and a pitch-and-roll WAVEC-buoy. The second experiment took place on 12 July 1992. During this experiment data were collected with the NASA/JPL polarimetric synthetic aperture radar (SAR), and a five-band helicopter-borne scatterometer. Again the hydro-meteorological conditions were monitored at MPN and the OCTANS. Furthermore, interferometric radar data were collected.
It has been demonstrated and recognized that radar interferometry is a promising method for the determination of digital elevation information and terrain slope from Synthetic Aperture Radar (SAR) data. An important application of Interferometric SAR (InSAR) data in areas with topographic variations is that the derived elevation and slope can be directly used for the absolute radiometric calibration of the amplitude SAR data as well as for scattering mechanisms analysis. On the other hand polarimetric SAR data has long been recognized as permitting a more complete inference of natural surfaces than a single channel radar system. In fact, imaging polarimetry provides the measurement of the amplitude and relative phase of all transmit and receive polarizations. On board the NASA DC-8 aircraft, NASA/JPL operates the multifrequency (P, L and C bands) multipolarimetric radar AIRSAR. The TOPSAR, a special mode of the AIRSAR system, is able to collect single-pass interferometric C- and/or L-band VV polarized data. A possible configuration of the AIRSAR/TOPSAR system is to acquire single-pass interferometric data at C-band VV polarization and polarimetric radar data at the two other lower frequencies. The advantage of this system configuration is to get digital topography information at the same time the radar data is collected. The digital elevation information can therefore be used to correctly calibrate the SAR data. This step is directly included in the new AIRSAR Integrated Processor. This processor uses a modification of the full motion compensation algorithm described by Madsen et al. (1993). However, the Digital Elevation Model (DEM) with the additional products such as local incidence angle map, and the SAR data are in a geometry which is not convenient, since especially DEMs must be referred to a specific cartographic reference system. Furthermore, geocoding of SAR data is important for multisensor and/or multitemporal purposes. In this paper, a procedure to geocode the new AIRSAR/TOPSAR data is presented. As an example an AIRSAR/TOPSAR image acquired in 1994 is geocoded and evaluated in terms of geometric accuracy.
Jurassicprok is a computer program that generates strip-map digital elevation models and other data products from raw data acquired by an airborne synthetic-aperture radar (SAR) system. This software can process data from a variety of airborne SAR systems but is designed especially for the GeoSAR system, which is a dual-frequency (P- and X-band), single-pass interferometric SAR system for measuring elevation both at the bare ground surface and top of the vegetation canopy. Jurassicprok is a modified version of software developed previously for airborne-interferometric- SAR applications. The modifications were made to accommodate P-band interferometric processing, remove approximations that are not generally valid, and reduce processor-induced mapping errors to the centimeter level. Major additions and other improvements over the prior software include the following: a) A new, highly efficient multi-stage-modified wave-domain processing algorithm for accurately motion compensating ultra-wideband data; b) Adaptive regridding algorithms based on estimated noise and actual measured topography to reduce noise while maintaining spatial resolution; c) Exact expressions for height determination from interferogram data; d) Fully calibrated volumetric correlation data based on rigorous removal of geometric and signal-to-noise decorrelation terms; e) Strip range-Doppler image output in user-specified Doppler coordinates; f) An improved phase-unwrapping and absolute-phase-determination algorithm; g) A more flexible user interface with many additional processing options; h) Increased interferogram filtering options; and i) Ability to use disk space instead of random- access memory for some processing steps.
Atmospheric errors due to the troposphere are a limiting error source for spaceborne interferometric synthetic aperture radar (InSAR) imaging. This software generates tropospheric delay maps that can be used to correct atmospheric artifacts in InSAR data. The software automatically acquires all needed GPS (Global Positioning System), weather, and Digital Elevation Map data, and generates a tropospheric correction map using a novel algorithm for combining GPS and weather information while accounting for terrain. Existing JPL software was prototypical in nature, required a MATLAB license, required additional steps to acquire and ingest needed GPS and weather data, and did not account for topography in interpolation. Previous software did not achieve a level of automation suitable for integration in a Web portal. This software overcomes these issues. GPS estimates of tropospheric delay are a source of corrections that can be used to form correction maps to be applied to InSAR data, but the spacing of GPS stations is insufficient to remove short-wavelength tropospheric artifacts. This software combines interpolated GPS delay with weather model precipitable water vapor (PWV) and a digital elevation model to account for terrain, increasing the spatial resolution of the tropospheric correction maps and thus removing short wavelength tropospheric artifacts to a greater extent. It will be integrated into a Web portal request system, allowing use in a future L-band SAR Earth radar mission data system. This will be a significant contribution to its technology readiness, building on existing investments in in situ space geodetic networks, and improving timeliness, quality, and science value of the collected data
Interferometric synthetic aperture radar (IFSAR) measurements at X-, C-, L-, and P-band are used to derive ground topography at meter level resolution. Interpretation of the derived topography requires attention due to the complex interaction of the radar signal with ground cover. The presence of penetrable surfaces such as vegetation, and tree canopies poses a challenge since the depth of penetration depends on a number of parameters such as the operating radar frequency, polarization, incident angle, as well as terrain structure. The dependence of the reconstructed topography on polarization may lead to the characterization of the ground cover. Simulation of interferometric measurements is useful for interpretation of the derived topography (B. Houshmand, Proceedings of URSI, 314, 1997). In this talk , time domain simulations for interferometric measurement for TE- and TM- polarization are presented. Time domain simulation includes the effects of the surface material property as well geometry comparable the radar signal wavelength (B. Houshmand, Proceedings of the URSI, 25, 1998). The IFSAR simulation is carried out in two steps. First, the forward scattering data is generated based on full wave analysis. Next, the electromagnetic information is inverted to generate surface topography. This inversion is based on the well known IFSAR processing technique which is composed of signal compression, and formation of an interferogram. The full wave forward scattering data is generated by the scattered-field formulation of the FDTD algorithm. The simulation is carried out by exciting the computational domain by a radar signal. The scattered field is then computed and translated to the receiving interferometric antennas using the time-domain Huygen's principle. The inversion process starts by compressing the time-domain data. The range compressed data from both receivers are then coregistered to form an interferogram. The resulting interferogram is then related to the ground topography using the radar imaging geometry. In this talk, the simulation results are compared with the C-band TM IFSAR derived topography, and the TE/TM SAR images at L-Band.