WO2000054006A2 - Single-pass interferometric synthetic aperture radar - Google Patents

Single-pass interferometric synthetic aperture radar Download PDF

Info

Publication number
WO2000054006A2
WO2000054006A2 PCT/US2000/005747 US0005747W WO0054006A2 WO 2000054006 A2 WO2000054006 A2 WO 2000054006A2 US 0005747 W US0005747 W US 0005747W WO 0054006 A2 WO0054006 A2 WO 0054006A2
Authority
WO
WIPO (PCT)
Prior art keywords
images
image
focal plane
values
recited
Prior art date
Application number
PCT/US2000/005747
Other languages
French (fr)
Other versions
WO2000054006A3 (en
Inventor
Carter M. Glass
Lori K. Harrison
Original Assignee
Lockheed Martin Corporation
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Lockheed Martin Corporation filed Critical Lockheed Martin Corporation
Priority to AU57219/00A priority Critical patent/AU5721900A/en
Publication of WO2000054006A2 publication Critical patent/WO2000054006A2/en
Publication of WO2000054006A3 publication Critical patent/WO2000054006A3/en

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021SAR image post-processing techniques
    • G01S13/9023SAR image post-processing techniques combined with interferometric techniques

Definitions

  • the present invention relates to the application of interferometric principles to determine terrain height from data collected by a synthetic aperture radar (SAR). It discloses a method of performing interferometric height estimation from data that may be collected on a single orbital pass of a radar-carrying satellite with a single transmit- receive antenna.
  • the invention may be implemented utilizing the techniques disclosed herein together with aspects of known imaging techniques, such as disclosed in C. Jakowatz, Jr., D. Wahl, P. Eichel, D. Ghiaglia and P. Thompson, SPOTLIGHT-MODE SYNTHETIC APERTURE RADAR: A Signal Process Approach (1996).
  • IFS AR is the primary technique which holds the promise of accurate, large area topographic mapping from space
  • a practical, cost-effective collection system has yet to be proposed.
  • the present invention provides a solution to the above-noted problems, and is referred to herein as "SPINS AR" imaging.
  • the SPINS AR imaging approach uses two complex images constructed from overlapping synthetic apertures.
  • the images may be formed on a single orbital satellite pass using a single antenna.
  • an interferometric baseline is created by the orbital curvature associated with the along-track separation of the synthetic aperture centers.
  • the height estimates are formed using coherent data collected during a contiguous time interval.
  • the SPINSAR invention offers additional benefits. For example, no phase unwrapping and no terrain tie points are needed. Further, height estimation accuracy and elevation post spacing are comparable or superior to those obtainable with other orbital interferometric schemes. Moreover, computational burden is comparable to that of other interferometric techniques.
  • the present invention provides a method for determining height information for a planetary surface region of interest, and includes the step of forming at least two complex images of the planetary surface region in different slant planes from synthetic aperture radar SAR image data, wherein the two images comprise partially overlapping and partially non-overlapping data content (i.e., due to the partially overlapping synthetic apertures corresponding with the two images). More particularly, the two SAR images are formed so as to define a convergence angle ⁇ c between the mid- aperture slant plane normal vectors of the two images (e.g., ⁇ c may be established between about .05° to 1° in typical applications).
  • each of the images should be obtained in spotlight-mode and utilize a common ground reference point (GRP) and focal-plane normal.
  • GRP ground reference point
  • the overlapping data content should comprise at least about three percent (3%) of the image data content of each, and more preferably, the overlapping data content should be from about five percent (5%) to fifty percent (50%) of the total image data content of each.
  • the inventive method further includes the steps of utilizing the at least two SAR images to obtain spatially variant differential layover values and employing the at least two SAR images to obtain interferometric phase values.
  • the spatially variant differential layover values and interferometric phase values may be utilized in an estimating step to yield height values for the planetary surface region of interest.
  • the inventive method may include the step of acquiring the image data employed to form the images utilizing a single synthetic aperture radar antenna.
  • a single synthetic aperture radar antenna may be located on a satellite traveling along an orbital path relative to the planetary surface region of interest.
  • the image data may be acquired during a single orbital pass of the satellite over the planetary surface region of interest.
  • the image data used to form the SAR images may be acquired during a continuous time interval, e.g., less than about 60 seconds, wherein temporal decorrelation is substantially avoided.
  • both overlapping and non-overlapping image data content portions of the at least two SAR images should be utilized. Further in this regard, it has been found that substantially all of the overlapping and non-overlapping data content portions may be advantageously employed. In this regard, and as will be further described, it has been recognized that the restriction or total avoidance of aperture trimming in the present invention allows for the obtainment of differential layover values from which height can be estimated and sufficient phase correlation for purposes of interferometric processing.
  • the inventive method may comprise the step of projecting the at least two slant plane SAR images into corresponding focal plane images.
  • sampled pixel values for each slant plane SAR image may be utilized to obtain corresponding projected focal plane image pixel values, wherein projected focal plane image pixel values may then be obtained at each pixel location via interpolation.
  • the method may further include the step of rotating a first of the focal plane images (i.e., a "slave" image), to match a focal plane grid of a second of the focal plane images (i.e., a "master” image). Thereafter, the method may include the step of removing positional phase differences between the rotated slave image and the master image at each pixel location within the master image to obtain an adjusted master plane image.
  • the rotated slave image and adjusted master image may be processed to produce corresponding amplitude-only images.
  • the amplitude-only images may be utilized for registration of the rotated slave image and adjusted master image.
  • the method may further include the step of decomposing the image data of each of the amplitude-only images into a corresponding plurality of sub-images, and correlating the sub-images corresponding with the rotated slave image and adjusted master image to obtain corresponding spatially variant differential layover values (i.e., resulting from the surface height of the region of interest).
  • the sub-image, spatially variant differential layover values may be interpolated to obtain values for each pixel location.
  • each of the rotated slave and adjusted master images may be decomposed into a corresponding plurality of image blocks. Then, each complex pixel in each image block corresponding with the adjusted master image may be multiplied by a complex conjugate of the corresponding pixel of the corresponding image block of the related slave image, wherein a plurality of complex-difference image blocks are obtained. In turn, such plurality of complex-difference image blocks may be utilized to form an interferogram.
  • interferogram formation may be achieved by first removing interferometric-phase equivalents of the previously determined spatially variant differential layover values for each pixel location of the complex difference image blocks . Then, a phase value at each pixel location may be extracted via conventional arctangent operation.
  • the interferogram may be conditioned (e.g., via filtering techniques) and transformed to height values. More particularly, a phase-to- height conversion factor may be applied to the interferogram to obtain the height values . Similarly, layover-to-height conversion factors may be applied to the above-noted differential layover values noted hereinabove. Then, statistical processing techniques may be utilized to obtain best fit planes for the interferometric height estimates and for the SAR height estimates. The best fit plane for the interferometric height estimates may then be employed with the interferometric height values to obtain residual interferometric heights which are adjusted to fit the SAR height estimate plane. The resulting adjusted values reflect the final height estimates for the planetary surface region of interest.
  • Both IFSAR and its relative SAR stereo are based on the differential layover of scattering centers which lie out of the focal plane.
  • the differential layover is induced by differences in the intersections of the iso-range spheres and iso-Doppler cones with the focal planes of an image pair. The intersections are different, thus generating differential layover proportional to out-of-plane height, provided the scattering center is "viewed" from suitably different aspects by the radar.
  • the curvature of the orbit produces differential layover proportional to out-of-plane height.
  • the SAR stereo principle estimates height by transforming differential layover directly to height. Interferometry, on the other hand, measures the phase equivalent of the slant-range difference between laid-over scattering centers on the two images and then compares this to the phase difference one would have had if the scattering centers had been in the focal planes.
  • the height-measurement sensitivity improves.
  • the scattering-center signatures begin to decorrelate in amplitude, making it increasingly difficult to register them for purposes of measuring interferometric phase.
  • the phase between corresponding scattering centers also decorrelates with increased viewing-angle difference, introducing an interferometric error source called "baseline decorrelation". So, as the along-orbit distance between aperture centers increases, estimate sensitivity improves but at the expense of degraded pixel registration and phase correlation. The result is a "bathtub" performance curve.
  • Baseline decorrelation is somewhat tractable in the case of sidelooking aircraft systems.
  • the angle diversity needed to generate differential layover is achieved by the grazing-angle difference between the boresight directions of the two antennas.
  • the line-of-sight wedges formed by the center of the imaged scene and the antenna positions at the beginning and end of each data collection interval are coincident in the along-track direction.
  • baseline decorrelation is caused only by the small grazing-angle difference in viewing the scene. Contrast this situation with that of two-pass satellite interferometry where the collection wedges may not be coincident in the along-track direction. They may even be disjoint.
  • baseline decorrelation is problematic for single-pass orbital interferometry for the following reasons: (1) if one forms two images with different along-track aperture centers and the collection regions are disj oint, then baseline decorrelation is total and the attempt fails; (2) if, on the other hand, the apertures overlap and are trimmed to eliminate the data which does not correlate, then the overlapping regions have the same grazing angle so that no differential layover is generated, and again interferometry fails. Fortunately, the consequences of baseline decorrelation are not so dire. To be sure, correlation does decade as one increases the along-track distance between aperture centers, and decorrelation is more or less total, depending on terrain content, when the apertures become disjoint.
  • the present inventors have determined that if one uses partially overlapped apertures to form two images from data collected on a single orbital pass, then sufficient phase correlation can be achieved for interferometry and the orbital curvature induces differential layover from which height can be estimated by interferometry.
  • the inventors have further recognized that in order to utilize the differential layover, aperture trimming should be restricted or omitted so that sufficient differential layover associated with the non-overlapping image regions is obtained, even at the price of reduced correlation.
  • Fig.2 illustrates sub-steps corresponding with process step 40 of the embodiment of Figs. lA - lC.
  • Fig.3 illustrates sub-steps corresponding with process step 50 of the embodiment of Figs. 1A - IC.
  • the illustrated SPINSAR embodiment of Figs. 1 -3 begins with the formation of at least two SAR images, in different slant planes, for a planetary surface region of interest (step 10).
  • the two images are formed so as to define a convergence angle ⁇ c between the mid-aperture slant plane normal vectors of the two images.
  • Each of the images should be obtained in spotlight-mode and utilize a common ground reference point (GRP) and focal-plane normal.
  • GRP ground reference point
  • the two SAR images are formed so as to have partially overlapping image data content (i.e., via partially overlapping apertures).
  • the overlapping data content will be at least about three percent (3%) and most preferably between about five percent (5%) and fifty percent (50%) of the data content of each.
  • the two SAR images may be formed using image data acquired via a single transmit/receive SAR antenna.
  • the image data is acquired by transmitting microwave energy pulses at a predetermined frequency toward the surface region of interest, and receiving resultant microwave energy reflected from the region.
  • an image-formation processor may be employed to complete the processing steps contemplated by the embodiment shown in Fig. 1. Such processing may be completed either on-board an imaging spacecraft, at a ground-based location, or at some other location remote from the imaging spacecraft.
  • the antenna may be located on a satellite orbiting the earth, wherein the two images are obtained during a single-pass over the region of interest.
  • the two complex images may each comprise a matrix of data samples, or pixels, of at least about 1000 x 1000.
  • image data acquired by the transmit/receive antenna may be readily processed to define the two images, wherein data obtained during a first time segment is used to create the first image and data obtained during a second time segment is used to form the second image, and wherein the first and second time segments partially overlap so as to define the degree of overlapping image data content.
  • the first and second time segments may comprise a single continuous image acquisition time frame during which imaging data is received by the transmit/receive antenna. By way of example, such acquisition time frame may be less than about 60 seconds.
  • the corresponding image data may be used to yield two corresponding, projected focal plane images (step 20).
  • the slant plane image data for each of the SAR images may be sampled and the sample values may be employed to obtain corresponding, projected focal plane values.
  • the projected focal plane values may be utilized to obtain image values at each pixel location.
  • one of the projected focal plane images may be referred to as a "master” while the other projected focal plane image may be referred to as "slave”.
  • the slave image data may then be rotated so as to match the focal plane grid of the master image data (step 22).
  • the rotated slave image and the master image are coincident except for differential layover induced by the height of objects within the imaged region of interest that lie outside of the focal plane.
  • differential layover can be extracted in further processing and used for both SAR height estimation and in interferometric processing.
  • residual position phase differences between the master and updated slave images may be determined and removed from the master image (step 30).
  • an image pixel at any given location in the images will generally have a residual position phase, i.e., the differential phase which would have been present if an object present at such location had been in the corresponding focal plane.
  • the position phase is the phase equivalent of the slant-range difference between a given pixel location and the common ground reference point (GRP).
  • GRP common ground reference point
  • the described process embodiment also provides for image registration and spatially variant differential layover determination (step 40).
  • the adjusted master image data from step 30 and rotated slave image data may be detected, or processed, to produce amplitude-only pixel values for each image, which data may be processed for correlation purposes and to obtain spatially variant differential layover values (step 40).
  • each of the adjusted master and rotated slave images may be processed to produce amplitude-only pixel values (step 41), and the amplitude-only images may be decomposed into "sub-images" of a predetermined size (step 43).
  • each image may be decomposed into a plurality of sub-images each having dimensions in the order of 100 x 100 pixels.
  • Such sub-images corresponding with each of the master and slave images may then be correlated to determine a spatially variant differential layover values for each sub-image occasioned by the height of the imaged regions of interest (step 45).
  • three steps may be completed in accordance with basic SAR stereo image processing principles.
  • the layover values for each of the sub- images may be median filtered to remove outliers and smoothed with a low-pass filter. Then, the layover values may be interpolated so as to provide a differential layover value for each of the pixel locations within each of the sub-images (step 47).
  • steps 45 and 47 may be repeated after further decomposing the master and slave images into corresponding, refined sub-images of an even lesser dimension than noted above (e.g., sub-images comprising 40 x 40 pixels), wherein the initially-determined differential layover values are employed in obtaining refined layover values.
  • step 40 further processing may be carried out with respect to the adjusted master and rotated slave images.
  • Such complex images may be decomposed into smaller corresponding "squares", or image clocks, wherein complex image values may be utilized with the spatially variant differential layover values obtained in step 40 to obtain complex-difference squares (step 50).
  • the adjusted master image and rotated slave image may be decomposed into squares (step 51 ) having dimensions which in effect determine the spacing of the independent elevation posts that may be employed.
  • image squares in the order of 7 x 7 pixels and 15 x 15 pixels are typical.
  • each complex pixel value may be multiplied by the complex conjugate of the associated pixel of the corresponding rotated slave image square, taking into account the corresponding spatially variant differential layover values for the given pixel obtained in step 40.
  • the complex-different squares obtained in step 50 may then be smoothed by averaging the in-phase and quadrature components (step 60).
  • an interferogram may be formed from the smoothed complex-difference squares (step 70). More particularly, such formation may be achieved by first removing the interferometric-phase equivalent of the spatially variant differential layover at each pixel location (i.e., as determined in step 40), followed by phase extraction via the conventional arctangent operation.
  • phase equivalent of the differential layover values avoids a need for phase unwrapping since the interferogram contains only the
  • the interferogram may be median filtered to remove outliers and smoothed with a low-pass filter (step 80). Then, the conditioned interferogram may be transformed to height values utilizing a phase-to-height conversion factor (step 80).
  • the phase-to-height conversion factor may be constructed via simulation using collection parameters determined for each imaging opportunity.
  • the above-described SAR stereo differential layover values obtained in step 40 may also be converted to height for each pixel location utilizing a layover-to-height conversion factor (step 82).
  • the layover-to-height conversion factor may be developed via simulation of the relevant imaging geometry.
  • the interferometric height estimates (i.e., from step 80) and SAR stereo height estimates (i.e., from step 82) may each be statistically processed to determine "best fit" planes for each (step 90). Such processing may utilize least-squares processing techniques. Then, the best fit plane for the interferometric estimates may be utilized with each of the interferometric estimated height values to obtain residual interferometric height values (step 92). In turn, such residual values may be adjusted to the best fit plane for the SAR estimated height values determined in step 90, thereby yielding final height estimates.
  • the SAR images referred to in step 10 above should be obtained via partially overlapped apertures in order to provide for the phase correlation needed for interferometric processing.
  • the imaged terrain consists of uniformly distributed, uncorrelated point scattering centers. This intuitively satisfying model has been used successfully in other SAR applications and is a worst case with respect to complex correlation.
  • the terrain is viewed from along-track aperture centers which differ by the slant-plane angle do. As it has been shown, see e.g., Zebker et al. supra, that the complex correlation coefficient for such terrain is given by: 2 - Wa
  • Wa is the azimuth resolution of the SAR and ⁇ is the radar wavelength.
  • the correlation coefficient is zero when: ⁇
  • the interferometric baseline for the described SPINSAR embodiment is generated by the angle ⁇ c between mid-aperture slant plane normal vectors of the two images formed in step 10 above.
  • ⁇ c is called the convergence angle.
  • the relationship between interferometric phase difference, ⁇ int, and out-of-plane height, h, has the form:
  • ⁇ h f( ⁇ c, x. y) - ⁇ ⁇ i nt .
  • ⁇ dif (xm, ym) M(xm, ym) • S(xs,
  • ⁇ pos( ⁇ m > y m ) - - • (RM ( ⁇ m ' y m ) - Rm G ) - ( R Q ( ⁇ q 5 y ) - R sG J
  • RmG slant range from mid-aperture of master image to the GRP
  • RsG slant range from mid-aperture of slave image to the GRP
  • RM(xm,Ym) slant range from mid-aperture of master image to point
  • RQ(xq, yq) slant range from mid-aperture of slave image to point Q(xq, yq), where Q(xq,yq) is the location of M(xm, ym) in the coordinate system of the slave image, i.e. the location of the scattering center at
  • RG is the range from mid-aperture to the GRP.
  • the associated slant-plane position phase is:
  • slant range, r, from the GRP is related to focal plane distance, y, from the GRP by:
  • ⁇ g is the mid-aperture grazing angle between slant and focal planes.
  • ⁇ os (x m s y ⁇ ) 4, ⁇ .ly m .cosl ⁇ g ⁇ J -y q .co ⁇ g s .
  • ⁇ rot be the angle between master and slave focal-plane y-axes. Then:
  • ⁇ pos(Xm,ym) ⁇ -- ym -(c ⁇ S( ⁇ gm) -COS( ⁇ rot) - COS( ⁇ g S )) -Xm -Sin( ⁇ rot)
  • Conversion Factors As indicated above, the described embodiment requires a conversion factor from interferometric phase difference to out-of-plane height (step 80). In addition, since we compensate the complex-difference image (i.e., obtained in step 50) with the phase equivalent of differential layover to avoid phase unwrapping (step 70), we also need the conversion from differential layover to equivalent interferometric phase difference. Neither of these conversion factors can be expressed as simple formulas for arbitrary radar parameters and imaging geometries. Instead they are developed by simulation for each imaging opportunity. The method is as follows:
  • Spotlight-mode imaging is assumed, with image formation based on the polar algorithm.
  • Tsx m [UR, x Ns m Ur m Ns m ]
  • Tsx [UR S x Ns s Ur s Ns s ]
  • Tsx Tsx " '
  • Tsx Tsx s " '
  • Focal-plane x-axis (cross-range) vectors
  • Tfx m [ufx m NfxUfx m Nf]
  • Tfx [ufx s Nf ⁇ Ufx s Nf]
  • Tfs wear Txs refuse # Tfx society
  • Tfs s Txs s » Tfx s
  • (xLm,yLm,0) be the layover point associated with the elevated point. Both of-these points project to the same (a,r) point in the slant plane; thus:
  • Equating (a,r) in these expressions gives the layover point in terms of the elevated point as follows:
  • the Interferometric Phase Difference and Conversion Factor (i.e., for step 80).
  • the slant- range differences from aperture center to the elevated and layover point for the master and slave images are as follows:
  • the interferometric phase difference is then the master-slave difference for the elevated point and the master-slave difference for the layover point.
  • the latter difference is the position phase which would have been observed if the scattering center had truly been at the layover point; thus,
  • phase-to-height conversion factor Kph is given by:
  • the conversion factor has a planar variation with (xm, ym); hence, it must be evaluated at three or more points fit to a plane, and interpolated at every interferogram point.
  • Tms Txf ⁇ • Tfx m

Abstract

A new process for interferometric height estimation using orbital synthetic aperture radar (SAR) is disclosed. The process uses coherent SAR data collected on a single pass (10) and may be collected by a single-antenna radar mounted on an orbiting satellite. The data is used to construct two overlapping synthetic apertures (10). The along-orbit separation of the synthetic aperture centers combines with orbital curvature to form a suitable interferometric baseline. The overlap provides adequate correlation to support image registration and phase-difference estimation (30). Single-pass collection eliminates temporal decorrelation between the interferometric image pair. It also eliminates the differential satellite ephemeris errors which plague two-pass and two-vehicle schemes. Both phase unwrapping and the need for elevation tie points are eliminated by pre-compensating the interferometric data with height estimates from a SAR-stereo pair formed from the same collected data (40).

Description

SINGLE-PASS INTERFEROMETRIC SYNTHETIC APERTURE RADAR
RELATED APPLICATIONS This application claims priority under 35 U.S.C. §119 to prior U.S. Provisional Patent Application Serial Number 60/123,269, filed March 8, 1999, the entirety of which is hereby incorporated by reference.
FIELD OF THE INVENTION The present invention relates to the application of interferometric principles to determine terrain height from data collected by a synthetic aperture radar (SAR). It discloses a method of performing interferometric height estimation from data that may be collected on a single orbital pass of a radar-carrying satellite with a single transmit- receive antenna. The invention may be implemented utilizing the techniques disclosed herein together with aspects of known imaging techniques, such as disclosed in C. Jakowatz, Jr., D. Wahl, P. Eichel, D. Ghiaglia and P. Thompson, SPOTLIGHT-MODE SYNTHETIC APERTURE RADAR: A Signal Process Approach (1996).
BACKGROUND OF THE INVENTION
The principles of interferometry are currently being used to extract terrain-height information from synthetic aperture radar images. This is done primarily using aircraft, but demonstrations with orbital SAR have also been somewhat successful. See e.g., Li, Fuk and R.M. Goldstein, "Studies of Multibaseline Spaceborne Interferometric Synthetic Aperture Radars", IEEE Transactions on Geoscience and Remote Sensing. Vol. 28, No. 1, January, 1990, pp. 888-97, and Goldstein, R.M., Zebker, A., and C.L. Werner, "Satellite Radar Interferometry: Two-Dimensional Phase Unwrapping", Radio Science, Vol. 23, No. 4, July-August, 1988, pp. 713-720. Aircraft practice implements the interferometric baseline by mounting two or more spatially separated radars, or one radar and extra receive antennas, on the wing tips and in the fuselage, for example.
To obtain terrain height for large, remote areas, aircraft interferometric synthetic aperture radar (IFSAR) systems have been proven to be expensive and potentially vulnerable. As a result, attention has turned to spaceborne platforms, but here difficulties arise. Optimal baselines from space are quite long, see e.g., Li et al. supra, so that physically separated antennas on a single satellite are impractical if height accuracies better than a few meters are needed. Multiple imaging passes from one or more radar- carrying vehicles are feasible, but because of uncertainties in vehicle ephemeris, numerous known elevation tie points are needed. A more severe problem is temporal terrain decorrelation between imaging passes. Ideally, one would like no more than a few seconds between collections. However, this is impractical with multiple passes, and often with multiple satellites.
As a consequence of these problems, even though IFS AR is the primary technique which holds the promise of accurate, large area topographic mapping from space, a practical, cost-effective collection system has yet to be proposed.
SUMMARY OF THE INVENTION Overview
The present invention provides a solution to the above-noted problems, and is referred to herein as " SPINS AR" imaging. The SPINS AR imaging approach uses two complex images constructed from overlapping synthetic apertures. Of note, the images may be formed on a single orbital satellite pass using a single antenna. In a primary application, an interferometric baseline is created by the orbital curvature associated with the along-track separation of the synthetic aperture centers. The height estimates are formed using coherent data collected during a contiguous time interval. In addition to the important advantages associated with a single pass and single radar approach (e.g., the temporal decorrelation interval is a maximum of a few seconds and height errors caused by differential vehicle ephemeris errors are substantially avoided.), the SPINSAR invention offers additional benefits. For example, no phase unwrapping and no terrain tie points are needed. Further, height estimation accuracy and elevation post spacing are comparable or superior to those obtainable with other orbital interferometric schemes. Moreover, computational burden is comparable to that of other interferometric techniques.
Characterization As indicated, the present invention provides a method for determining height information for a planetary surface region of interest, and includes the step of forming at least two complex images of the planetary surface region in different slant planes from synthetic aperture radar SAR image data, wherein the two images comprise partially overlapping and partially non-overlapping data content (i.e., due to the partially overlapping synthetic apertures corresponding with the two images). More particularly, the two SAR images are formed so as to define a convergence angle θc between the mid- aperture slant plane normal vectors of the two images (e.g., θc may be established between about .05° to 1° in typical applications). Further, each of the images should be obtained in spotlight-mode and utilize a common ground reference point (GRP) and focal-plane normal. Preferably, the overlapping data content should comprise at least about three percent (3%) of the image data content of each, and more preferably, the overlapping data content should be from about five percent (5%) to fifty percent (50%) of the total image data content of each.
The inventive method further includes the steps of utilizing the at least two SAR images to obtain spatially variant differential layover values and employing the at least two SAR images to obtain interferometric phase values. In turn, the spatially variant differential layover values and interferometric phase values may be utilized in an estimating step to yield height values for the planetary surface region of interest.
In conjunction with forming the at least two SAR images, and as noted above, the inventive method may include the step of acquiring the image data employed to form the images utilizing a single synthetic aperture radar antenna. Such antenna may be located on a satellite traveling along an orbital path relative to the planetary surface region of interest. In turn, the image data may be acquired during a single orbital pass of the satellite over the planetary surface region of interest. Further, the image data used to form the SAR images may be acquired during a continuous time interval, e.g., less than about 60 seconds, wherein temporal decorrelation is substantially avoided.
It should be noted that, in both the utilizing and employing steps identified above (i.e., to obtain spatially variant differential layover values and interferometric phase values, respectively), both overlapping and non-overlapping image data content portions of the at least two SAR images should be utilized. Further in this regard, it has been found that substantially all of the overlapping and non-overlapping data content portions may be advantageously employed. In this regard, and as will be further described, it has been recognized that the restriction or total avoidance of aperture trimming in the present invention allows for the obtainment of differential layover values from which height can be estimated and sufficient phase correlation for purposes of interferometric processing. In a further aspect of the present invention, the inventive method may comprise the step of projecting the at least two slant plane SAR images into corresponding focal plane images. By way of example, sampled pixel values for each slant plane SAR image may be utilized to obtain corresponding projected focal plane image pixel values, wherein projected focal plane image pixel values may then be obtained at each pixel location via interpolation. After forming the at least two focal plane images (i.e., corresponding with the at least two slant plane SAR images, the method may further include the step of rotating a first of the focal plane images (i.e., a "slave" image), to match a focal plane grid of a second of the focal plane images (i.e., a "master" image). Thereafter, the method may include the step of removing positional phase differences between the rotated slave image and the master image at each pixel location within the master image to obtain an adjusted master plane image.
The rotated slave image and adjusted master image may be processed to produce corresponding amplitude-only images. In turn, the amplitude-only images may be utilized for registration of the rotated slave image and adjusted master image. In this regard, the method may further include the step of decomposing the image data of each of the amplitude-only images into a corresponding plurality of sub-images, and correlating the sub-images corresponding with the rotated slave image and adjusted master image to obtain corresponding spatially variant differential layover values (i.e., resulting from the surface height of the region of interest). The sub-image, spatially variant differential layover values may be interpolated to obtain values for each pixel location.
After registration, and in another aspect of the inventive method, each of the rotated slave and adjusted master images may be decomposed into a corresponding plurality of image blocks. Then, each complex pixel in each image block corresponding with the adjusted master image may be multiplied by a complex conjugate of the corresponding pixel of the corresponding image block of the related slave image, wherein a plurality of complex-difference image blocks are obtained. In turn, such plurality of complex-difference image blocks may be utilized to form an interferogram. In conjunction with this aspect, interferogram formation may be achieved by first removing interferometric-phase equivalents of the previously determined spatially variant differential layover values for each pixel location of the complex difference image blocks . Then, a phase value at each pixel location may be extracted via conventional arctangent operation.
Following interferogram formation, the interferogram may be conditioned (e.g., via filtering techniques) and transformed to height values. More particularly, a phase-to- height conversion factor may be applied to the interferogram to obtain the height values . Similarly, layover-to-height conversion factors may be applied to the above-noted differential layover values noted hereinabove. Then, statistical processing techniques may be utilized to obtain best fit planes for the interferometric height estimates and for the SAR height estimates. The best fit plane for the interferometric height estimates may then be employed with the interferometric height values to obtain residual interferometric heights which are adjusted to fit the SAR height estimate plane. The resulting adjusted values reflect the final height estimates for the planetary surface region of interest.
Additional aspects and extensions of the present invention will be apparent to those skilled in the art upon further consideration of the descriptions that follow.
Theoretical Basis
Both IFSAR and its relative SAR stereo are based on the differential layover of scattering centers which lie out of the focal plane. The differential layover is induced by differences in the intersections of the iso-range spheres and iso-Doppler cones with the focal planes of an image pair. The intersections are different, thus generating differential layover proportional to out-of-plane height, provided the scattering center is "viewed" from suitably different aspects by the radar. When the images are formed with aperture centers separated along a single orbital pass, as per the SPINSAR embodiment herein, the curvature of the orbit produces differential layover proportional to out-of-plane height. The SAR stereo principle estimates height by transforming differential layover directly to height. Interferometry, on the other hand, measures the phase equivalent of the slant-range difference between laid-over scattering centers on the two images and then compares this to the phase difference one would have had if the scattering centers had been in the focal planes.
As one increases the along-orbit distance between the two synthetic aperture centers, the height-measurement sensitivity improves. Unfortunately, as the distance increases, the scattering-center signatures begin to decorrelate in amplitude, making it increasingly difficult to register them for purposes of measuring interferometric phase. Moreover, the phase between corresponding scattering centers also decorrelates with increased viewing-angle difference, introducing an interferometric error source called "baseline decorrelation". So, as the along-orbit distance between aperture centers increases, estimate sensitivity improves but at the expense of degraded pixel registration and phase correlation. The result is a "bathtub" performance curve.
Baseline decorrelation is somewhat tractable in the case of sidelooking aircraft systems. Here the angle diversity needed to generate differential layover is achieved by the grazing-angle difference between the boresight directions of the two antennas. In this situation, the line-of-sight wedges formed by the center of the imaged scene and the antenna positions at the beginning and end of each data collection interval are coincident in the along-track direction. Hence, baseline decorrelation is caused only by the small grazing-angle difference in viewing the scene. Contrast this situation with that of two-pass satellite interferometry where the collection wedges may not be coincident in the along-track direction. They may even be disjoint. According to current theory, and depending on how one models "diffuse" terrain, only the overlapping portion of the collection wedges is useful for interferometry. See e.g., Jakowatz, C.V., Wahl, D.E., Ghliga, D.C., and P.A. Thompson, Spotlight-Mode Svnthetic Aperture Radar: A Signal Processing Approach, Kluwer Academic Publishers,
1996, Boston, pp. 283,285, and Mareshal, N., "Tomographic Formulation of Interferometric SAR for Terrain Elevation Mapping", IEEE Transactions on Geoscience and Remote Sensing, Vol. 33, No. 3, May 1995, pp. 726-739. The non-overlapping data does not correlate and is believed to contribute only noise. Accordingly, standard interferometric processing practice is to discard the non-overlapping data, a process called "aperture trimming" or "wavenumber filtering". In the extreme case where the collection wedges are disjoint, i.e. the synthetic apertures do not overlap, none of the data correlates, and interferometry presumably cannot be done.
In view of the foregoing, it may appear that baseline decorrelation is problematic for single-pass orbital interferometry for the following reasons: (1) if one forms two images with different along-track aperture centers and the collection regions are disj oint, then baseline decorrelation is total and the attempt fails; (2) if, on the other hand, the apertures overlap and are trimmed to eliminate the data which does not correlate, then the overlapping regions have the same grazing angle so that no differential layover is generated, and again interferometry fails. Fortunately, the consequences of baseline decorrelation are not so dire. To be sure, correlation does decade as one increases the along-track distance between aperture centers, and decorrelation is more or less total, depending on terrain content, when the apertures become disjoint. See e.g., Zebker, H.A., "Decorrelation in Interferometric Radar Echoes", IEEE Transactions on Geoscience and Remote Sensing, Vol. 30, No. 5, September 1992, pp. 950-959. Nevertheless, enough correlation remains to support excellent IFSAR imaging results from single-pass collections.
For example, successful single-image height estimation algorithms based on the proportionality between residual dispersed-domain quadratic phase and out-of-plane object height often use the so-called map-drift algorithm which depends on the residual amplitude correlation of disjoint sub-aperture images. Single-pass orbital SAR stereo also depends on amplitude correlation of disjoint collections with even larger angle diversity, and height accuracies of a few meters have been achieved. On the other hand, both theory and empirical evidence suggest that phase correlation is more sensitive than amplitude correlation to viewing-angle difference. As a result, the present inventors have recognized that SAR stereo imaging apertures must have at least a few percent of overlap to support single-pass IFSAR.
More particularly, the present inventors have determined that if one uses partially overlapped apertures to form two images from data collected on a single orbital pass, then sufficient phase correlation can be achieved for interferometry and the orbital curvature induces differential layover from which height can be estimated by interferometry. In this regard, the inventors have further recognized that in order to utilize the differential layover, aperture trimming should be restricted or omitted so that sufficient differential layover associated with the non-overlapping image regions is obtained, even at the price of reduced correlation. By using this approach SPINSAR has achieved height accuracies and elevation post spacings of a couple of meters from authentic SAR data.
DESCRIPTION OF THE DRAWINGS Figs. 1 A, IB and IC illustrate process steps of one embodiment of the inventive method of the present invention.
Fig.2 illustrates sub-steps corresponding with process step 40 of the embodiment of Figs. lA - lC.
Fig.3 illustrates sub-steps corresponding with process step 50 of the embodiment of Figs. 1A - IC.
DETAILED DESCRIPTION Referring to Figs. 1 - 3, one embodiment of the present invention will be described. That description will be followed by a review of underlying algorithms which support the described embodiment. Additional embodiments and adaptations will be apparent to those skilled in the art. Process Description The illustrated SPINSAR embodiment of Figs. 1 -3 begins with the formation of at least two SAR images, in different slant planes, for a planetary surface region of interest (step 10). In this regard, the two images are formed so as to define a convergence angle θc between the mid-aperture slant plane normal vectors of the two images. Each of the images should be obtained in spotlight-mode and utilize a common ground reference point (GRP) and focal-plane normal. Of note, the two SAR images are formed so as to have partially overlapping image data content (i.e., via partially overlapping apertures). Preferably, the overlapping data content will be at least about three percent (3%) and most preferably between about five percent (5%) and fifty percent (50%) of the data content of each. In a primary application, the two SAR images may be formed using image data acquired via a single transmit/receive SAR antenna. The image data is acquired by transmitting microwave energy pulses at a predetermined frequency toward the surface region of interest, and receiving resultant microwave energy reflected from the region. After image data acquisition and storage, an image-formation processor may be employed to complete the processing steps contemplated by the embodiment shown in Fig. 1. Such processing may be completed either on-board an imaging spacecraft, at a ground-based location, or at some other location remote from the imaging spacecraft.
In earth imaging applications, the antenna may be located on a satellite orbiting the earth, wherein the two images are obtained during a single-pass over the region of interest. By way of example, the two complex images may each comprise a matrix of data samples, or pixels, of at least about 1000 x 1000. As will be appreciated, image data acquired by the transmit/receive antenna may be readily processed to define the two images, wherein data obtained during a first time segment is used to create the first image and data obtained during a second time segment is used to form the second image, and wherein the first and second time segments partially overlap so as to define the degree of overlapping image data content. The first and second time segments may comprise a single continuous image acquisition time frame during which imaging data is received by the transmit/receive antenna. By way of example, such acquisition time frame may be less than about 60 seconds.
After the two slant plane images are formed (step 10), the corresponding image data may be used to yield two corresponding, projected focal plane images (step 20). By way of example, the slant plane image data for each of the SAR images may be sampled and the sample values may be employed to obtain corresponding, projected focal plane values. In turn, the projected focal plane values may be utilized to obtain image values at each pixel location. For purposes of further description, one of the projected focal plane images may be referred to as a "master" while the other projected focal plane image may be referred to as "slave". The slave image data may then be rotated so as to match the focal plane grid of the master image data (step 22). At this point, the rotated slave image and the master image are coincident except for differential layover induced by the height of objects within the imaged region of interest that lie outside of the focal plane. As will be described, such differential layover can be extracted in further processing and used for both SAR height estimation and in interferometric processing. In the later regard, and in order to facilitate interferometric phase extraction, residual position phase differences between the master and updated slave images may be determined and removed from the master image (step 30). In this regard, it will be understood that an image pixel at any given location in the images will generally have a residual position phase, i.e., the differential phase which would have been present if an object present at such location had been in the corresponding focal plane. For spotlight- mode imaging using a polar algorithm, the position phase is the phase equivalent of the slant-range difference between a given pixel location and the common ground reference point (GRP). As such, removal of the position phase differences between the master and related slave images may be achieved via algorithmic processing.
As indicated in Fig. 1A, the described process embodiment also provides for image registration and spatially variant differential layover determination (step 40). In this regard, the adjusted master image data from step 30 and rotated slave image data may be detected, or processed, to produce amplitude-only pixel values for each image, which data may be processed for correlation purposes and to obtain spatially variant differential layover values (step 40). More particularly, and referring to Fig. 2, each of the adjusted master and rotated slave images may be processed to produce amplitude-only pixel values (step 41), and the amplitude-only images may be decomposed into "sub-images" of a predetermined size (step 43). By way of example, each image may be decomposed into a plurality of sub-images each having dimensions in the order of 100 x 100 pixels. Such sub-images corresponding with each of the master and slave images may then be correlated to determine a spatially variant differential layover values for each sub-image occasioned by the height of the imaged regions of interest (step 45). As will be appreciated by those skilled in the art, three steps may be completed in accordance with basic SAR stereo image processing principles. The layover values for each of the sub- images may be median filtered to remove outliers and smoothed with a low-pass filter. Then, the layover values may be interpolated so as to provide a differential layover value for each of the pixel locations within each of the sub-images (step 47). Additionally, in order to enhance precision, steps 45 and 47 may be repeated after further decomposing the master and slave images into corresponding, refined sub-images of an even lesser dimension than noted above (e.g., sub-images comprising 40 x 40 pixels), wherein the initially-determined differential layover values are employed in obtaining refined layover values.
Referring now to Fig. IB, further steps relating to the interferometric processing aspects of the described embodiment will be reviewed. Specifically, after registration step 40, further processing may be carried out with respect to the adjusted master and rotated slave images. Such complex images may be decomposed into smaller corresponding "squares", or image clocks, wherein complex image values may be utilized with the spatially variant differential layover values obtained in step 40 to obtain complex-difference squares (step 50). More particularly, and with reference now to Fig. 3, the adjusted master image and rotated slave image may be decomposed into squares (step 51 ) having dimensions which in effect determine the spacing of the independent elevation posts that may be employed. By way of example, image squares in the order of 7 x 7 pixels and 15 x 15 pixels are typical. Then, for each adjusted master image square, each complex pixel value may be multiplied by the complex conjugate of the associated pixel of the corresponding rotated slave image square, taking into account the corresponding spatially variant differential layover values for the given pixel obtained in step 40. The complex-different squares obtained in step 50 may then be smoothed by averaging the in-phase and quadrature components (step 60).
Next, an interferogram may be formed from the smoothed complex-difference squares (step 70). More particularly, such formation may be achieved by first removing the interferometric-phase equivalent of the spatially variant differential layover at each pixel location (i.e., as determined in step 40), followed by phase extraction via the conventional arctangent operation.
As will be appreciated, removal of the phase equivalent of the differential layover values avoids a need for phase unwrapping since the interferogram contains only the
"delta" between the true height and the coarse height values determinable by the SAR stereo algorithm referred to in step 40 hereinabove. Further, it should be noted that, as a consequence of the short baseline utilized in the described embodiment, the phase equivalent of the deltas will be less than 360°. Further algorithmic explanation of the described phase equivalent removal is provided hereinbelow. Following formation, the interferogram may be median filtered to remove outliers and smoothed with a low-pass filter (step 80). Then, the conditioned interferogram may be transformed to height values utilizing a phase-to-height conversion factor (step 80).
The phase-to-height conversion factor may be constructed via simulation using collection parameters determined for each imaging opportunity.
In conjunction with the described embodiment, the above-described SAR stereo differential layover values obtained in step 40 may also be converted to height for each pixel location utilizing a layover-to-height conversion factor (step 82). Again, the layover-to-height conversion factor may be developed via simulation of the relevant imaging geometry.
Next, the interferometric height estimates (i.e., from step 80) and SAR stereo height estimates (i.e., from step 82) may each be statistically processed to determine "best fit" planes for each (step 90). Such processing may utilize least-squares processing techniques. Then, the best fit plane for the interferometric estimates may be utilized with each of the interferometric estimated height values to obtain residual interferometric height values (step 92). In turn, such residual values may be adjusted to the best fit plane for the SAR estimated height values determined in step 90, thereby yielding final height estimates.
Algorithmic Support
1. Aperture Overlap and Baseline Decorrelation
As indicated above, the SAR images referred to in step 10 above should be obtained via partially overlapped apertures in order to provide for the phase correlation needed for interferometric processing. In this regard, assume the imaged terrain consists of uniformly distributed, uncorrelated point scattering centers. This intuitively satisfying model has been used successfully in other SAR applications and is a worst case with respect to complex correlation. Assume further that the terrain is viewed from along-track aperture centers which differ by the slant-plane angle do. As it has been shown, see e.g., Zebker et al. supra, that the complex correlation coefficient for such terrain is given by: 2 - Wa
0) p = ι - dφ λ
where Wa is the azimuth resolution of the SAR and λ is the radar wavelength. The correlation coefficient is zero when: λ
(2) dφ
2 - Wa
Recall, however, that in order to obtain clutter cells with uniformly weighted impulse response width Wa, the synthetic aperture must subtend a slant-plane integration angle θi given by:
Λ . 0.886 - λ
(3) θι = .
2 • Wa
Substituting Eq. (3) into Eq. (2) yields:
(4) dφ = θ i .
Ψ 0.886
This expression tells us that the correlation coefficient diminishes to zero when the synthetic aperture centers are separated approximately by the integration angle. Thus, the correlation becomes zero at the point where the apertures become disjoint. Accordingly, to preserve the correlation necessary for interferometry the coherent data collection interval must be decomposed such that the two apertures have some overlap.
The interferometric baseline for the described SPINSAR embodiment is generated by the angle θc between mid-aperture slant plane normal vectors of the two images formed in step 10 above. In the SAR-stereo community θc is called the convergence angle. The relationship between interferometric phase difference, φint, and out-of-plane height, h, has the form:
(5) h = f(θc, x, y) - φ int ,
where f(.) decreases monotonically with θc, which in turn increases with dφ. The error standard deviations are related by:
(6) σ h = f(θc, x. y) - σ φ int . To minimize σh we need a large convergence angle, and hence large dφ; but this conflicts with the need for correlation. The result is a bathtub response for σh versus dφ. The only parameter which improves both convergence and correlation simultaneously is the azimuth resolution Wa, the finer the better. These comments are as far as one can go with theory in the general case. This is because p is unpredictable for authentic terrain (although the model used here is thought to be worst case). Moreover, f(θc, x, y) is a complicated function of radar parameters and imaging geometry which must be computed by simulation for each imaging opportunity.
2. Position Phase
As indicated above, residual phase differences between the master and slave images are removed from the master image in step 30 prior to the interferometric processing steps 50-70 of the described embodiment.
In this regard, suppose a scattering center lies at location M(xm,ym) in the master image and at S(xs,ys) in the slave. Then the phase difference is given by:
θ dif (xm, ym) = M(xm, ym) S(xs,
Figure imgf000016_0001
Figure imgf000016_0002
where Arg(.) is the argument and the overbar denotes the complex conjugate. To get the interferometric phase, one must subtract from θdif the phase difference one would have observed at M(xm,ym) if it corresponded to an in-plane scattering center. This is called the position phase θpos(xm,ym). It is given by:
θ pos(χm> ym) = - - (RM (χm' ym) - RmG ) - (RQ (χq5 y ) - RsG J
where:
RmG = slant range from mid-aperture of master image to the GRP, RsG = slant range from mid-aperture of slave image to the GRP, RM(xm,Ym) = slant range from mid-aperture of master image to point
M(xm,ym), and, RQ(xq, yq) = slant range from mid-aperture of slave image to point Q(xq, yq), where Q(xq,yq) is the location of M(xm, ym) in the coordinate system of the slave image, i.e. the location of the scattering center at
M(xm,ym) if the differential layover had been zero.
Let (a,r) be azimuth and range measured from the GRP in the range-Doppler slant plane. Then the first two terms of a binomial expansion of the Pythagorean Theorem give a good approximation for the range difference between (1) the mid-aperture vehicle location and P(a,r) and (2) the mid-aperture vehicle location and the GRP. The result is: a 2 + r 2 dR = r +
2 RG
where RG is the range from mid-aperture to the GRP. The associated slant-plane position phase is:
Figure imgf000017_0001
For typical imaging geometry RG is much larger than the dimensions of the image, in which case we can further approximate: θPos(a,r) = 4_π -r
But slant range, r, from the GRP is related to focal plane distance, y, from the GRP by:
r = y-cos(φg) ,
where Ψg is the mid-aperture grazing angle between slant and focal planes. Hence:
4-π- us g) θ pos(χ,y)= λ φ -y
Thus, we can write the expression for position-phase difference between master and slave images as:
θ os(xm syπι) = 4,π .lym.coslφgπιJ -yq .co^φgs. Let θrot be the angle between master and slave focal-plane y-axes. Then:
yq = xm -sin(θrot) - ym • cos(θrot)
and:
θpos(Xm,ym) = ^-- ym -(cθS(φgm) -COS(θrot) - COS(φgS)) -Xm -Sin(θrot)
This expression shows that the position phase in the complex-difference image, which is in the focal plane of the master image, is a plane passing through zero phase at the GRP where xm =ym = 0 and having slopes:
- cos(φgs) -cos(θrot)) .
Figure imgf000018_0001
3. Conversion Factors As indicated above, the described embodiment requires a conversion factor from interferometric phase difference to out-of-plane height (step 80). In addition, since we compensate the complex-difference image (i.e., obtained in step 50) with the phase equivalent of differential layover to avoid phase unwrapping (step 70), we also need the conversion from differential layover to equivalent interferometric phase difference. Neither of these conversion factors can be expressed as simple formulas for arbitrary radar parameters and imaging geometries. Instead they are developed by simulation for each imaging opportunity. The method is as follows:
Except where noted, all vector quantities are given in Earth-centered fixed (ECF) axes in which the rotation of the Earth has been referred to the radar-carrying satellite.
Spotlight-mode imaging is assumed, with image formation based on the polar algorithm.
Simulation inputs:
Pm,Ps Mid-aperture Earth-center-to-satellite position vectors for the master and slave images.
G Earth-center-to-GRP vector.
Vm,Vs Mid-aperture master and slave satellite velocity vectors. λ Transmit wavelength.
M(xm,Ym) Master image focal-plane locations at which conversion factors are desired.
Coordinate transformations and unit vectors:
Mid-aperture slant ranges and slant-plane unit vectors:
Figure imgf000019_0001
Slant-plane normal vectors: Nsm = ,Xm x R"1, Nss = ,XS x Rs,
|Vm X Rm| |Vs X Rsj
Slant-plane to EC1: transformations:
Tsxm= [UR, x Nsm Urm Nsm ] Tsx= [URS x Nss Urs Nss ]
Tsx =Tsx "' Tsx=Tsxs "'
Focal-plane normal:
Figure imgf000019_0002
Focal-plane x-axis (cross-range) vectors:
UfXm = Rm Nf ufχs = ^l f
|RmxNfj |RsxNf|
Focal-plane to ECF transformations: Tfxm=[ufxm NfxUfxm Nf] Tfx=[ufxs NfχUfxs Nf]
Focal -to-slant plane transformations:
Tfs„=Txs„#Tfx„ Tfss=Txss »Tfxs Assume a scattering center elevated by h units at location (xm,ym,h) in the focal plane of the master image. Let (xLm,yLm,0) be the layover point associated with the elevated point. Both of-these points project to the same (a,r) point in the slant plane; thus:
Figure imgf000020_0002
Equating (a,r) in these expressions gives the layover point in terms of the elevated point as follows:
xLm Tfsm >0,0 yLm Tfs ml,0
Figure imgf000020_0001
where we have used projections rather than iso-range spheres and iso-Doppler cones with negligible error for realistic imaging geometries.
The Interferometric Phase Difference and Conversion Factor (i.e., for step 80). The slant- range differences from aperture center to the elevated and layover point for the master and slave images are as follows:
Figure imgf000021_0001
Figure imgf000021_0003
dR m =
Figure imgf000021_0004
where the vertical bars denote the magnitude of the enclosed vector quantity. The interferometric phase difference is then the master-slave difference for the elevated point and the master-slave difference for the layover point. The latter difference is the position phase which would have been observed if the scattering center had truly been at the layover point; thus,
4 - π
Φmt- -.f' dRhm - dRh H dRLm - dRLc
The phase-to-height conversion factor Kph is given by:
Figure imgf000021_0002
For typical image sizes and geometries the conversion factor has a planar variation with (xm, ym); hence, it must be evaluated at three or more points fit to a plane, and interpolated at every interferogram point.
Conversion From Differential Layover to Interferometric Phase Difference (i.e., for step 70). The differential layover found by a spatially variant correlation of small regions in the master and slave images (i.e., via step 40) is transformed to equivalent interferometric phase difference and used to compensate the complex difference image before the phase is extracted (step 70). The focal plane location of the elevated point in the slave image is denoted as (xs,ys,hs). The transformation from master to slave focal plane is:
Tms = Txf ς • Tfx m
Then:
Figure imgf000022_0002
The out-of-plane elevations are identical in the master and slave images since the focal planes have the same normal vector and differ by only a rotation. The layover point (xLs,yLs) in the focal plane of the slave image is then:
1 xs xLc TfsSQ,0 TfSs0,l Tfs Tfs Tfs s0,0 s0,l s0,2 ys yLs Tfs Tfs si,0 SU Tfssi,0 Tfssi,l Tfssi,2 h
Before we can compute the differential layover, this point must be referred back to the master image. Denote this referred point by (xLsm,yLsm); then,
Figure imgf000022_0001
The differential layover is therefore:
Figure imgf000023_0001
Three conversions from differential layover to height are available:
dLx m dLy m dLx m dLy m
KxLh KyLh KepiLh =
Φint
These conversions have negligle variation with (xm,ym). In practice, it is computationally efficient to rotate both master and slave images by the epipolar angle— the angle which makes ICx or Ky zero—in which case the layover is all in the epipolar direction. In this situation the correlation needed to determine the differential layover is one-dimensional rather than two-dimensional. When epipolar correlation is used, one obtains the layover associated with the numerator of Kepi and uses this expression for the conversion. (The correct sign must be applied to the square root.) Finally, the conversion from layover to equivalent interferometric phase is obtained by combining the conversion factors:
Kepi h
KLφ = κφh(χm > ym )
* * *
The above-described embodiment is not intended to limit the scope of the present invention. Other embodiments, extensions and modifications will be apparent to those skilled in the art. For example, the invention may be employed in conjunction with 3 or more SAR images to yield height estimates for a region of interest. All such extensions are intended to be within the scope of the present invention as defined by the claims that follow.

Claims

What is Claimed is:
1. A method for determining height information for a planetary surface region comprising: forming at least two images of a planetary surface region in different slant planes from synthetic aperture radar image data, wherein said at least two images comprise partially overlapping and partially non-overlapping data content; utilizing said at least two images to obtain spatially variant differential layover values; employing said at least two images to obtain interferometric phase values; and estimating height values for said planetary surface region utilizing said spatially variant differential layover values and said interferometric phase values.
2. A method as recited in Claim 1, wherein said forming step comprises: acquiring said image data utilizing a single synthetic aperture radar antenna.
3. A method as recited in Claim 2, wherein said single synthetic aperture radar antenna is located on a satellite on an orbital path relative to said planetary surface region, and wherein said image data is acquired during a single orbital pass of said satellite over said planetary surface region.
4. A method as recited in Claim 2, wherein said image data used to form said at least two images is acquired during a continuous time interval.
5. A method as recited in Claim 4, wherein said continuous time interval is less than about 60 seconds.
6. A method as recited in Claim 1, wherein said partially overlapping data content comprises at least about three percent (3%) of the image data corresponding with each of said at least two images.
7. A method as recited in Claim 6, wherein said overlapping data content comprises between about five percent (5%) and fifty percent (50%) of the image data corresponding with each of said at least two images.
8. A method as recited in Claim 1 , where said utilizing and employing steps each comprise: using said overlapping and non-overlapping data content portions of said at least two images.
9. A method as recited in Claim 1, wherein substantially all of said overlapping and non-overlapping data content portions of said at least two images are used in each of said utilizing and employing steps.
10. A method as recited in Claim 1 , further comprising: projecting each of said at least two slant plane images to obtain corresponding focal plane images; rotating a first of said focal plane images to match a focal plane grid of a second of said focal plane images; and removing residual position phase differences between said rotated first focal plane image and said second focal plane image at each pixel location within said second focal plane image to obtain an adjusted second focal plane image.
1 1. A method as recited in Claim 10, wherein said utilizing step comprises: processing said rotated first focal plane image and said adjusted second focal plane image to produce corresponding amplitude-only images; decomposing each of said amplitude-only images into a corresponding plurality of sub-images; correlating said corresponding plurality of sub-images to determine said spatially variant differential layover values.
12. A method as recited in Claim 11 , wherein said employing step comprises: decomposing each of said first and adjusted second focal plane images into a corresponding plurality of image blocks; multiplying each complex pixel in each image block corresponding with said adjusted second focal plane image by a complex conjugate of a corresponding pixel of a corresponding image block of the first focal plane image, wherein a plurality of complex-difference image blocks are obtained; and using said plurality of complex-difference image squares to form an interferogram.
13. A method as recited in Claim 12, wherein said using step comprises: removing interferometric-phase equivalents of said spatially variant differential layover values for each pixel location of said complex-image squares; and extracting said interferometric phase values for each pixel location of said complex image squares.
14. A method as recited in Claim 13, wherein said estimating step comprises: transforming said interferogram to obtain first estimated height values; and utilizing said spatially variant differential layover values to obtain second estimated height values.
15. A method as recited in Claim 14, wherein said estimating step further comprises: statistically determining a first best fit plane for said first estimated height values and a second best fit plane for said second estimated height values; using said first best fit plane and said first estimated height values to obtain residual height values; and adjusting said residual height values to said second best fit plane to obtain estimated height values for said planetary surface region.
16. A method for determining height information for a planetary surface region, comprising: acquiring synthetic aperture radar image data with a single synthetic aperture radar antenna located on a satellite traveling along an orbital path relative to the planetary surface region, wherein said image data is acquired in a spotlight-mode during a continuous time interval of a single orbital pass of said satellite over said planetary surface region; forming at least two images of said planetary surface region in different slant planes, wherein said at least two images comprise partially overlapping and partially non- overlapping data content and have a common ground reference point; registering said at least two images and obtaining spatially variant differential layover values; employing said at least two images to form an interferogram; extracting interferometric phase values from said interferogram; and estimating height values for said planetary surface region utilizing said spatially variant differential layover values and said interferometric phase values.
17. A method as recited in Claim 16, wherein said partially overlapping data content comprises between about five percent (5%) and fifty percent (50%) of the image data corresponding with each of said at least two images.
18. A method as recited in Claim 17, wherein substantially all of said overlapping and non-overlapping data content portions of said at least two images are used in each of said registering and employing steps.
19. A method as recited in Claim 16, further comprising: projecting each of said at least two slant plane images to obtain corresponding focal plane images; rotating a first of said focal plane images to match a focal plane grid of a second of said focal plane images; and removing residual position phase differences between said rotated first focal plane image and said second focal plane image at each pixel location within said second focal plane image to obtain an adjusted second focal plane image, hwerein said first and second adjusted focal plane images are employed in said registering step.
20. A method as recited in Claim 19, further comprising: processing said rotated first focal plane image and said adjusted second focal plane image to produce corresponding amplitude-only images; decomposing each of said amplitude-only images into a corresponding plurality of sub-images; correlating said corresponding plurality of sub-images to determine said spatially variant differential layover values.
21. A method as recited in Claim 20, wherein said employing step comprises: decomposing each of said first and adjusted second focal plane images into a corresponding plurality of image blocks; multiplying each complex pixel in each image block corresponding with said adjusted second focal plane image by a complex conjugate of a corresponding pixel of a corresponding image block of the first focal plane image, wherein a plurality of complex-difference image blocks are obtained; and using said plurality of complex-difference image squares to form an interferogram.
22. A method as recited in Claim 16, wherein said estimating step comprises: employing said interferometric phase values to obtain first estimated height values; and utilizing said spatially variant differential layover values to obtain second estimated height values.
23. A method as recited in Claim 22, wherein said estimating step further comprises: statistically determining a first best fit plane for said first estimated height values and a second best fit plane for said second estimated height values; using said first best fit plane and said first estimated height values to obtain residual height values; and adjusting said residual height values to said second best fit plane to obtain estimated height values for said planetary surface region.
PCT/US2000/005747 1999-03-08 2000-03-03 Single-pass interferometric synthetic aperture radar WO2000054006A2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
AU57219/00A AU5721900A (en) 1999-03-08 2000-03-03 Single-pass interferometric synthetic aperture radar

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US12326999P 1999-03-08 1999-03-08
US60/123,269 1999-03-08

Publications (2)

Publication Number Publication Date
WO2000054006A2 true WO2000054006A2 (en) 2000-09-14
WO2000054006A3 WO2000054006A3 (en) 2001-01-18

Family

ID=22407677

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2000/005747 WO2000054006A2 (en) 1999-03-08 2000-03-03 Single-pass interferometric synthetic aperture radar

Country Status (2)

Country Link
AU (1) AU5721900A (en)
WO (1) WO2000054006A2 (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007071113A1 (en) * 2005-12-22 2007-06-28 Qifeng Yu A three-images isoline correlation interferometric method for generating an interferogram of radar
US8154435B2 (en) 2008-08-22 2012-04-10 Microsoft Corporation Stability monitoring using synthetic aperture radar
CN102495407A (en) * 2011-11-14 2012-06-13 电子科技大学 Characterization method of similarity of polarized synthetic aperture radar image
CN102955157A (en) * 2011-08-26 2013-03-06 中国科学院空间科学与应用研究中心 Fast correlation coefficient method for interferometric synthetic aperture radar image precise registration
CN104297751A (en) * 2014-10-20 2015-01-21 西安电子科技大学 Main and auxiliary SAR image fusion method of distributed satellite SAR system
WO2015151134A1 (en) * 2014-04-04 2015-10-08 三菱電機株式会社 Radar signal processing device
WO2016201050A1 (en) * 2015-06-12 2016-12-15 Delta ID Inc. Apparatuses and methods for iris based biometric recognition
CN110618409A (en) * 2019-09-09 2019-12-27 长沙理工大学 Multi-channel InSAR interferogram simulation method and system considering overlapping and shading
CN110865372A (en) * 2018-08-27 2020-03-06 中国人民解放军61646部队 Target height information extraction method based on synthetic aperture radar multi-azimuth observation
CN111239736A (en) * 2020-03-19 2020-06-05 中南大学 Single-baseline-based surface elevation correction method, device, equipment and storage medium
EP3757610A4 (en) * 2018-03-30 2021-04-21 Mitsubishi Electric Corporation Radar image processing device and radar image processing method

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5332999A (en) * 1993-02-05 1994-07-26 Agence Spatiale Europeenne Process for generating synthetic aperture radar interferograms
US5552787A (en) * 1995-10-10 1996-09-03 The United States Of America As Represented By The Secretary Of The Navy Measurement of topography using polarimetric synthetic aperture radar (SAR)
US5608405A (en) * 1995-10-06 1997-03-04 Lockheed Martin Corporation Method of generating visual representation of terrain height from SAR data employing multigrid analysis
US5620708A (en) * 1993-02-22 1997-04-15 Alza Corporation Compositions and methods for the oral delivery of active agents
US5659318A (en) * 1996-05-31 1997-08-19 California Institute Of Technology Interferometric SAR processor for elevation
US5680138A (en) * 1995-12-15 1997-10-21 Lockheed Martin Corporation Synthetic aperture radar simulation
US5835055A (en) * 1996-03-20 1998-11-10 Atlantis Scientific Inc. Method for iterative disk masking and automatic error repair for phase unwrapping
US6011625A (en) * 1998-07-08 2000-01-04 Lockheed Martin Corporation Method for phase unwrapping in imaging systems

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5332999A (en) * 1993-02-05 1994-07-26 Agence Spatiale Europeenne Process for generating synthetic aperture radar interferograms
US5620708A (en) * 1993-02-22 1997-04-15 Alza Corporation Compositions and methods for the oral delivery of active agents
US5608405A (en) * 1995-10-06 1997-03-04 Lockheed Martin Corporation Method of generating visual representation of terrain height from SAR data employing multigrid analysis
US5552787A (en) * 1995-10-10 1996-09-03 The United States Of America As Represented By The Secretary Of The Navy Measurement of topography using polarimetric synthetic aperture radar (SAR)
US5680138A (en) * 1995-12-15 1997-10-21 Lockheed Martin Corporation Synthetic aperture radar simulation
US5835055A (en) * 1996-03-20 1998-11-10 Atlantis Scientific Inc. Method for iterative disk masking and automatic error repair for phase unwrapping
US5659318A (en) * 1996-05-31 1997-08-19 California Institute Of Technology Interferometric SAR processor for elevation
US6011625A (en) * 1998-07-08 2000-01-04 Lockheed Martin Corporation Method for phase unwrapping in imaging systems

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007071113A1 (en) * 2005-12-22 2007-06-28 Qifeng Yu A three-images isoline correlation interferometric method for generating an interferogram of radar
US8154435B2 (en) 2008-08-22 2012-04-10 Microsoft Corporation Stability monitoring using synthetic aperture radar
CN102955157A (en) * 2011-08-26 2013-03-06 中国科学院空间科学与应用研究中心 Fast correlation coefficient method for interferometric synthetic aperture radar image precise registration
CN102955157B (en) * 2011-08-26 2014-04-02 中国科学院空间科学与应用研究中心 Fast correlation coefficient method for interferometric synthetic aperture radar image precise registration
CN102495407A (en) * 2011-11-14 2012-06-13 电子科技大学 Characterization method of similarity of polarized synthetic aperture radar image
CN102495407B (en) * 2011-11-14 2013-03-13 电子科技大学 Characterization method of similarity of polarized synthetic aperture radar image
WO2015151134A1 (en) * 2014-04-04 2015-10-08 三菱電機株式会社 Radar signal processing device
CN104297751A (en) * 2014-10-20 2015-01-21 西安电子科技大学 Main and auxiliary SAR image fusion method of distributed satellite SAR system
CN104297751B (en) * 2014-10-20 2017-03-08 西安电子科技大学 A kind of major-minor SAR image fusion method of Distributed Spaceborne SAR System
WO2016201050A1 (en) * 2015-06-12 2016-12-15 Delta ID Inc. Apparatuses and methods for iris based biometric recognition
EP3757610A4 (en) * 2018-03-30 2021-04-21 Mitsubishi Electric Corporation Radar image processing device and radar image processing method
CN110865372A (en) * 2018-08-27 2020-03-06 中国人民解放军61646部队 Target height information extraction method based on synthetic aperture radar multi-azimuth observation
CN110618409A (en) * 2019-09-09 2019-12-27 长沙理工大学 Multi-channel InSAR interferogram simulation method and system considering overlapping and shading
CN111239736A (en) * 2020-03-19 2020-06-05 中南大学 Single-baseline-based surface elevation correction method, device, equipment and storage medium

Also Published As

Publication number Publication date
AU5721900A (en) 2000-09-28
WO2000054006A3 (en) 2001-01-18

Similar Documents

Publication Publication Date Title
US9417323B2 (en) SAR point cloud generation system
Brenner et al. Radar imaging of urban areas by means of very high-resolution SAR and interferometric SAR
CA2248420C (en) Improved method for phase unwrapping in imaging systems
Ignatenko et al. ICEYE Microsatellite SAR Constellation Status Update: Evaluation of first commercial imaging modes
Scheiber et al. Coregistration of interferometric SAR images using spectral diversity
Sansosti et al. Geometrical SAR image registration
EP0654681B1 (en) Hyper-precision SAR interferometry using a dual-antenna multi-pass SAR system
US5923278A (en) Global phase unwrapping of interferograms
JP6320577B2 (en) Synthetic aperture radar signal processor
Ulander et al. Ultra-wideband SAR interferometry
US20040004569A1 (en) Method and system for mutual coherent synthetic aperture radiometry
Magnard et al. Analysis of a maximum likelihood phase estimation method for airborne multibaseline SAR interferometry
CN110018474A (en) Three-D imaging method based on geostationary orbit synthetic aperture radar chromatographic technique
Mancon et al. On the phase calibration by multisquint analysis in TOPSAR and stripmap interferometry
WO2000054006A2 (en) Single-pass interferometric synthetic aperture radar
Margot et al. Digital elevation models of the Moon from Earth-based radar interferometry
Rigling et al. Three-dimensional surface reconstruction from multistatic SAR images
EP3757610B1 (en) Radar image processing device and radar image processing method
Nitti et al. Evaluation of DEM-assisted SAR coregistration
Madsen et al. EMISAR single pass topographic SAR interferometer modes
Eineder et al. SAR interferometry with TerraSAR-X
Sommer et al. Image reconstruction without an inertial navigation system using backprojection autofocus for synthetic aperture radar
Moccia Fundamentals of bistatic synthetic aperture radar
Li et al. An approach for spaceborne InSAR DEM inversion integrated with stereo-SAR method
Ding et al. Deramp range migration processing for space-borne spotlight synthetic aperture radar

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A2

Designated state(s): AE AL AM AT AU AZ BA BB BG BR BY CA CH CN CR CU CZ DE DK DM EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX NO NZ PL PT RO RU SD SE SG SI SK SL TJ TM TR TT TZ UA UG US UZ VN YU ZA ZW

AL Designated countries for regional patents

Kind code of ref document: A2

Designated state(s): GH GM KE LS MW SD SL SZ TZ UG ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE BF BJ CF CG CI CM GA GN GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
AK Designated states

Kind code of ref document: A3

Designated state(s): AE AL AM AT AU AZ BA BB BG BR BY CA CH CN CR CU CZ DE DK DM EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX NO NZ PL PT RO RU SD SE SG SI SK SL TJ TM TR TT TZ UA UG US UZ VN YU ZA ZW

AL Designated countries for regional patents

Kind code of ref document: A3

Designated state(s): GH GM KE LS MW SD SL SZ TZ UG ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE BF BJ CF CG CI CM GA GN GW ML MR NE SN TD TG

REG Reference to national code

Ref country code: DE

Ref legal event code: 8642

122 Ep: pct application non-entry in european phase