[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
Role of SAM Chain Length in Enhancing the Sensitivity of Nanopillar Modified Electrodes for Glucose Detection
Next Article in Special Issue
Using a LIDAR Vegetation Model to Predict UHF SAR Attenuation in Coniferous Forests
Previous Article in Journal
A New Method for Node Fault Detection in Wireless Sensor Networks
Previous Article in Special Issue
A Fast Level Set Method for Synthetic Aperture Radar Ocean Image Segmentation
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Sensors
Volume 9
Issue 2
10.3390/s90201259
sensors-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Improvement of the Accuracy of InSAR Image Co-Registration Based On Tie Points – A Review

by
Weibao Zou
1,*,
Yan Li
2,
Zhilin Li
3 and
Xiaoli Ding
3
1
Shenzhen Institute of Advanced Technology, P.R. China
2
The University of Southern Queensland, Australia
3
The Hong Kong Polytechnic University, Hong Kong, P.R. China
*
Author to whom correspondence should be addressed.
Sensors 2009, 9(2), 1259-1281; https://doi.org/10.3390/s90201259
Submission received: 12 January 2009 / Revised: 11 February 2009 / Accepted: 17 February 2009 / Published: 24 February 2009
(This article belongs to the Special Issue Synthetic Aperture Radar (SAR))
Browse Figures

<p>Imaging geometry of SAR interferometry.</p> ">
<p>Imaging geometry of repeat-pass SAR interferometry.</p> ">
<p>Imaging geometry of across-track SAR interferometry.</p> ">
<p>Imaging geometry of along-track SAR interferometry.</p> ">
<p>SAR image processing procedures.</p> ">
<p>Procedures of InSAR image co-registration.</p> ">
<p>A pair of SAR images (a) and (b) and its co-registered slave image. (a) a master image; (b) a slave image; (c) a co-registered slave image.</p> ">
<p>Feature points extracted by different ways. (a) feature points extracted by Interest Operators; (b) feature points extracted by wavelet transform.</p> ">
<p>Feature points extracted by different ways. (a) feature points extracted by Interest Operators; (b) feature points extracted by wavelet transform.</p> ">
<p>A pair of SAR images (Tai Lam in Hong Kong). (a) the master image; (b) the slave image.</p> ">
Versions Notes

Abstract

:
Interferometric Synthetic Aperture Radar (InSAR) is a new measurement technology, making use of the phase information contained in the Synthetic Aperture Radar (SAR) images. InSAR has been recognized as a potential tool for the generation of digital elevation models (DEMs) and the measurement of ground surface deformations. However, many critical factors affect the quality of InSAR data and limit its applications. One of the factors is InSAR data processing, which consists of image co-registration, interferogram generation, phase unwrapping and geocoding. The co-registration of InSAR images is the first step and dramatically influences the accuracy of InSAR products. In this paper, the principle and processing procedures of InSAR techniques are reviewed. One of important factors, tie points, to be considered in the improvement of the accuracy of InSAR image co-registration are emphatically reviewed, such as interval of tie points, extraction of feature points, window size for tie point matching and the measurement for the quality of an interferogram.

1. Introduction

Interferometric Synthetic Aperture Radar (InSAR) has undergone rapid development since its first proposal by Graham in 1974. Within more than three decades, InSAR has gained increasing attention from researchers in a variety of areas due to its advantages of all weather condition, rapid and accurate topographic data collection [1]. As a result, it has matured to be a widely applied technique in many fields [2], such as surface deformation monitoring [3-5]; forest management and classification [6-9]; ocean current and glacier movement monitoring [10-12]; hydrologic studies [13] and polar research [14,15].
InSAR is a potential technique for generating digital elevation models (DEM) by using the phase component of the complex radar signal [16]. A DEM is one of the most demanded products in the remote sensing community [16]. It is widely used within the geoscientific community, e.g. for mapping purposes, geomorphologic studies based on slope and aspect maps, and as a layer in geographical information systems (GIS) for combining relief data with thematic information. A properly equipped spaceborne InSAR system can be used to produce a highly accurate global DEM with its advantages in significantly less time and at significantly lower cost than other systems [17]. This greatly attracts researchers' attention in military field. In the navigation field, a highly accurate DEM is one of the important information for improving the accuracy of the navigational system [17-21].
InSAR system is used on the Shuttle Radar Topography Mission (SRTM) in February 2000 which is sponsored by the NASA and National Geospatial-Intelligence Agency (NGA) to acquire spatially-continuous elevation information over 80% of the Earth's land mass in a single 10-day Space Shuttle flight [22-26]. The produced DTMs are with 30 m (1 arc-second, SRTM1) and 90 m (3 arc-seconds, SRTM3) grid size and its vertical accuracy ranges from a few meters to decameters [27-28]. It is the first spaceborne single-pass interferometric SAR and will produce the first near-global, high-resolution, digital elevation map. Such a global map will be constructed significantly sooner than with other systems by taking advantage of the unique opportunity offered through augmentation of the previously flown NASA SIR-C and the X-SAR [22].
Since then, a variety of research has being carried out on SRTM data. One kind of research is on the assessment of the accuracy of SRTM DEM. The vertical accuracy of the SRTM Digital Terrain Elevation Data (DTED) level 1 in Southern Greece was assessed by a reference DEM produced from contour lines digitization [23]. Elevations generated by NASA's Laser Vegetation Imaging Sensor (LVIS) were used to assess the validation of SRTM elevations at five study sites of variable relief and land-cover [24]. The results show that there is a linear relationship between the SRTM-LVIS elevation differences and canopy vertical extent [24]. The research carried out by Rodriguez et al. documents the results of the SRTM validation effort using global data set generated by SRTM [29]. An accuracy assessment of the C-band SRTM DEM was performed through comparisons with independent elevation data obtained from ICESat and terrestrial height information. It is found that the comparisons can be greatly influenced by dynamic changes of the Earth's surface or instrument-based specifications [30]. Walker et al. presented a comprehensive application-specific assessment of the quality of SRTM C-and X-band DEMs in terms of vertical accuracy. The testing results indicated that the quality of SRTM data may be higher than previously thought [25]. It was implemented to evaluate and compare the utility of SRTM and Advanced Spaceborne Thermal Emission and Reflection Radiometer DEMs with baseline DTED-1 topography for mapping lahar flooding hazards from the volcano of Citlaltepetl, Mexico [31]. A comparison of water stages derived from LiDAR, topographic contours and SRTM was implemented. A surprisingly good performance by the SRTM suggests that this is a potentially valuable source for initial flood information extraction in large, topographically homogeneous floodplains [32].
Another is on the comparison between SRTM DEM and other resource derived DEMs. The elevation differences between SRTM C-band 1 and 3 arcsecond resolution DEMs and ICESat 1,064 nm altimeter channel elevation data generated in areas of variable topography and vegetable cover were studied [33]. Geomorphometric comparison between SRTM and the National Elevation Dataset (NED) was implemented to demonstrate how the two sources represent the same surface morphology [34]. The image differences between DEMs are sensitive to even small amounts of misregistration between the two DEMs [35]. The various levels of misregiatration have different impacts on image differences between SRTM and DEMs [35].
Moreover, some research has focused on the applications of SRTM DEMs in many fields. The SRTM C- and X-band DEM data for Ohio and the Amazon were used to analyze the capability of measuring water elevations. The results suggest a great potential for space-based, laterally-spatial measurements of water surface elevations [36]. The accuracy of the Consultative Group for International Agriculture Research (CGIAR) DEM was assessed by comparing elevation values from processed Consultative Group for International Agriculture Research Consortium for Spatial Information (CGIAR-CSI) SRTM data [37]. The SRTM DEM was also used as a reference to provide an adjustment of other elevation data with a precision of five meters [38]. The SRTM data was used to validate the quality of CORONA DEMs [39]. The evaluation of SRTM data was implemented for lahar modeling [40]. An archaeological survey based on the SRTM terrain model was conducted to detect the ancient settlement mounds in the Near East [41]. For many data voids contained in SRTM DEM, a method to fill these voids, based on Triangular Irregular Networks, using a fill surface derived from Russian military maps is presented in [42]. Based on experiments with SRTM data and a few of other data sources, Schetselaar et al. presented a data processing and integration methodology of geological field mapping of the LOT 2/LOT 3 project areas and their extensions [43]. A new quantitative topographic algorithm, called HAND (Height Above the Nearest Drainage), based on SRTM DEM was developed. The HAND terrain descriptor produces a normalized digital elevation model that can be applied to classify terrain in a manner that is related to local soil water conditions. This increases usability of the SRTM DEM and provides a new quantitative view on the steady state landscape, one that was missing in the repertoire of terrain descriptors [26].
In short, InSAR is enjoying widespread applications in many fields [3-21]. It is a promising technique to solve some difficult problems in a study area accurately, economically, efficiently and rapidly. There are, however, numerous critical factors affecting the quality of InSAR data. One such factor is the SAR image co-registration that determines both the robustness and the accuracy of InSAR technology and limits its applications and developments. There is no doubt that improving the accuracy of image co-registration in InSAR would greatly benefit the development of InSAR technology.
After this section, the principle of SAR interferommetric processing is introduced. Then, the techniques for improving the accuracy of InSAR image co-registration are described in Section 3, which includes the basic schedule of InSAR image co-registration, determination of tie point interval, extraction of feature points, strategy of tie point matching, determination of window size, selection of transformation model and measurement of the interferogram quality. Finally, a summary is drawn.

2. Synthetic Aperture Radar Interferometric Processing

2.1. Acquisition of SAR Images

Normally, there are three ways to acquire Synthetic Aperture Radar (SAR) images, according to the number of antenna mounted on the platform and the orientation of the formed baseline [3, 16]: repeat-pass, across-track and along-track. There is no difference on the procedure of InSAR image co-registration regardless of the interferometric way the images are acquired.
In order to understand the principle of SAR image acquisition, the general InSAR imaging geometry is illustrated in Figure 1. Two antennas s1 and s2 on ideally parallel flight paths are separated by a baseline B. The slant range from antenna s1 to an illuminated point on the ground is r1; r2 is the slant range from antenna s2 to the same point. With the looking angle θ, the angle of baseline with respect to horizontal line α and the flying height H, the imaging geometry of InSAR is fixed.
The height can be calculated by the following equations:
sin ( θ − α ) = r 1 2 − r 2 2 + B 2 2 r 1 B ≈ − λ φ unwrapped 2 π B
θ = α − arcsin [ λ φ unwrapped 4 π B ]
h = H − r 1 cos θ
where H is the antenna s1 is the height above the reference plane and h is the target height. r1 and r2 are the slant ranges from antennae to the point on the ground, respectively. θ is the looking angle. α is the angle of baseline with respect to horizontal line. φunwrapped is the unwrapped interferometric phase. B is the baseline. λ is the wavelength of radar signal. Therefore, the DEM data can be derived from radar interferogram.

2.1.1. Acquiring SAR Images by Repeat-pass

The term “repeat-pass” in this paper means to pass the same area repeatedly. In this approach, only one antenna is used and the SAR system acquires SAR data sets by passing the same area twice, covering it with a slightly different viewing geometry. The imaging geometry of repeat-pass SAR interferometry is illustrated in Figure 2. From the figure, it is noticed that, on the first pass, the radar wave is transmitted from antenna s1, and after interaction with the terrain, the backscattered return is also recorded by antenna s1. The signal is then processed to a complex SAR image, so does s2 on the second pass. The signal is processed to another complex SAR image. These SAR images are called the master image and the slave image, respectively, as shown later.
This approach is most suited to spaceborne systems. This is because the precise location of the flight path is required, and satellites typically have much more precise and stable orbital paths in the absence of the atmosphere than aircrafts. Normally, the repeat-pass approach is adopted to obtain SAR images to generate a DEM [44-47]. It is also used to do surface deformation measurements [3-5].

2.1.2. Acquiring SAR Images by Across-track

Another approach to acquire SAR images is by across-track. Its imaging geometry is illustrated in Figure 3. In this approach, two SAR antennae, S1 and S2, are mounted on the same platform. The SAR data are acquired at the same time and the baseline, B, is perpendicular to the flight path. The system, therefore, records the master image and the slave image at the same time but from two slightly different locations.
The main problem with the geometry of the across-track approach in the airborne configuration due to errors caused by the aircraft roll cannot be distinguished from the influence of the terrain slope. This problem is less critical in the spaceborne case, as a satellite track is more stable than an airborne flight path.
The across-track interferometer is also used for DEM generation [48]. In addition, it can also be available for sea surface mapping [49]. A system combining along-track and across track InSAR is carried out for measurement of ocean waves [50].

2.1.3. Acquiring SAR Images by Along-track

The along-track means to acquire SAR images by two SAR antennae, s1 and s2. The two antennae follow each other at a short distance on the same orbit and possibly on the same platform for SAR data acquisition. The baseline, B, is parallel to the flight direction. Its imaging geometry is illustrated in Figure 4. At present, the along-track approach is only applicable to airborne SAR systems, as it requires two antennae on the same platform. The along-track approach is sensitive to movement of the scatters in the range direction. It can be used to measure Earth-surface velocities. As such, this technique holds promise for the detection of slowly moving ground targets. It was developed in mapping ocean currents. Above sea, it produces a map of the current velocity component parallel to the radar look direction. Now the master image and the slave image are recorded from the same position but at slightly different times. Therefore, the phase differences between the corresponding signals are caused by the movement of the measured object, e.g. water current. The moving surface leads to a Doppler shift according to the phase velocity of the water waves. All stationary targets are not visible whereas the moving ones can be seen in the radar imagery. The along-track mode can be used to measure the velocity of ocean currents and on-ground objects [51,52].
The quality of interferometric measurements depends on the relative distance of the antenna(s), the so-called base line. Across-track and along-track are only possible using airborne systems with two antennae. The yaw and pitch cause baseline components in y- and z-direction produce additional phase differences. The base line can, therefore, be readily found to the required precision. In principle, repeat-pass way is also possible for aircraft, but the flight track and the aircraft attitude should be known with sufficient accuracy.

2.2. Interferometric Processing of SAR Images

There are five steps in the procedures for processing InSAR images. They can be summarized as shown in Figure 5 as follows:
  • Acquisition of SAR images;
  • Co-registration of two SAR images;
  • Generation of interferogram;
  • Phase unwrapping of the wrapped interferometric phase;
  • Geocoding of DEM.
  • DEM is the final product of InSAR image processing.

3. Procedures of InSAR Image Co-registration

3.1. Co-registration of InSAR Image

A lot of research has been carried out to reduce the errors that are caused by interferometric processing. The InSAR processing consists of co-registration, interferogram generation, phase unwrapping and geocoding. Research has been carried out on co-registration [44,45,53-56], generation of interferograms [44], phase unwrapping [58-61], and geocoding [62].
The co-registration of SAR complex images is the first step in interferometric processing and it is one of the most important processing procedures involved in InSAR. For the quality of the final products, accurate co-registration of the two input images is a prerequisite. Based on a good image registration, ideally up to 1/10 of a pixel, interferometric products with reliable quality could be achieved.
The SAR image co-registration chain consists of two steps: i.e., coarse co-registration and fine co-registration. The fine co-registration is an important step, which includes selection of tiepoints, tiepoint matching, transformation models and resampling method shown in Figure 6. A pair of SAR images and its co-registered slave image are shown in Figure 7, respectively. The distribution and the window size of tie points for image co-registration will greatly affect the final generation of SAR image co-registration. Some researches have been implemented on these aspects [44,45,53-56].
An automatic approach with multi-step image matching algorithm was proposed by Liao [63]. In the first step, some interesting tie points have been selected by an operator in a grid form in the master image. The grid size is 9 × 9. The tie point matching was carried out based on a normalized correction coefficient. A set of window sizes during tie point matching is adopted to do experiments. Good results of InSAR image processing are obtained based on the window size, 63 × 63. Therefore, it was selected as the optimum window size for tie point matching. The next step based on relaxation technique is used to ensure the reliability. At last, least-squares adjustment is carried out to improve the accuracy of co-registration. The multi-step strategy is recommended as the first selection for image co-registration.
Another method for image co-registration in the Fourier domain, using Fast Fourier Transform (FFT), is performed [64]. Indeed, this method works very fast. However, it is just suitable for the stationary regions. For non-stationary regions, the result is unreliable.
There are three kinds of issues in image co-registration. The first is the selection of tie point that includes two factors: the distribution of tie points; and the number of tie point to be selected. For the distribution of tie points, theoretically, more tie points will result in a more reliable registration. It has been stated in many literatures that tie points with features should be selected and evenly distributed over the whole image. However, too many tie points would result in a dramatic increase in calculation but it does not necessarily result in a more reliable co-registration. The second is the tie point matching that includes two factors: the window size, and the method for matching. The third is the mathematical models for transformations. They will be described in the following sub-sections, respectively.

3.2. Determination of Interval of Tie Points for Co-registration

Conventionally, tie points are usually selected in a grid form, i.e. at grid nodes. Tie points grid sizes of 2×10, 9×9, 6×6 and 8×8 have been adopted by Zebker et al. [1], Wan [65], Yang and Wang [66] and Liu et al. [67], respectively. At present, the selection of tie points is determined mainly by experience.
The interval of tie points greatly influences the quality of image co-registration and the accuracy of final product, DEM. The optimum interval of tie points for image co-registration is systematically conducted [44, 68]. In the study, only the interval of tie points is considered while other factors are kept unchanged. Four pairs of SAR images of 1,760×400 pixel size of the Hong Kong area were used for testing. The effects of interval of tie points are assessed by a relative measure for the quality of the resultant interferogram and an absolute measure for the accuracy of the resultant DEM. Results show that the effect of tie point interval on the accuracy of the final DEM is not linear. When the tie point interval is smaller than 273×44 pixels (273 pixels in row and 44 pixels in column), the variation in the resultant DEM accuracy is not significant. It is also noticeable that an interval of 205×34 pixels always results in the best or very good results. An interval of around 200×30 pixels (200 pixels in row and 30 pixels in column) is an appropriate choice for the selection of tie points for image co-registration in hilly areas.

3.3. Extraction of Feature Points for Tie Point Matching

Generally, the points with special features, such as intersection and airports, are possibly selected as tie points and if feature points evenly distributes in the image will improve the reliability of co-registration. However, such feature points can not always be found in an image. Therefore, the connotation of feature points should be enlarged. Some methods have been used to extract feature points from images. One of the methods is Interest Operators [69], which is based on the maximum gradient to extract feature points as shown in Figure 8(a). However, each pixel gradient is calculated based on the original whole image. Therefore, too much calculation is needed.
Wavelets can also be used to extract feature points from images [70]. Therefore, an enlarged connotation of feature points in wavelet domain for InSAR image co-registration is developed by Zou [71]. A method to extract feature points with wavelet transform to serve as tie points for co-registration is implemented shown in Figure 8(b). The extraction of feature points is based on the wavelet gradient modulus maximum. Compared with Interest Operators, with wavelets the wavelet gradient is implemented not in the original image but at its highest level so that it reduces the complexity of computation and makes the processing fast. In order to make the co-registration more reliable, the feature points extracted by wavelet and the points selected in grid form are combined together used as tie points [71]. Good experimental results have been achieved. Therefore, this can be considered as a recommended operation for feature extraction.

3.4. Strategy for Tie Point Matching

The second issues in image co-registration should be considered is the tie point matching that includes two factors: the window size, and the strategy for matching. Generally, there are two kinds of methods for tie point matching: one is based on spectral maximization; another is the image cross-correlation maximization [72].
For the first method, a small window centered around the tie point in the master image is selected, and it is moved on the slave image. When the peak of the Fourier Transform of the product of the master image multiplied by the complex conjugate of the slave image is maximized, the matching is regarded as the best one. A large amount of calculation is carried out in the method.
About the second method, the cross-correlation coefficient is used as a parameter for evaluating the matching results. The principle is that given a window of pixels on the master image, there will be many possible matching windows of the same size, but the one matching with the highest cross-correlation coefficient is regarded as the best matching. The central pixels of the two windows are then regarded as the corresponding image points found.
Let z1(m, n) and z2(m, n) represent the two SAR complex images:
z 1 ( m , n ) = a 1 ( m , n ) e j φ 1 ( m , n )
z 2 ( m , n ) = a 2 ( m , n ) e j φ 2 ( m , n )
where m = 0, 1, 2, Λ, M-1, n = 0, 1, 2, Λ, N-1, and the size of the image is M×N. The cross-correlation coefficient γ is expressed as:
γ = | E ( z 1 z 2 ∗ ) | 2 E ( | z 1 | 2 ) E ( | z 2 | 2 )
The expanded form of equation (5) can be written as:
γ ( m , n ) = | ∑ l = − ( L − 1 ) / 2 ( L − 1 ) / 2 ∑ k = − ( K − 1 ) / 2 ( K − 1 ) / 2 z 1 ( m + l , n + k ) z 2 ∗ ( m + l , n + k ) e − j φ ( m + l , n + k ) | ∑ l = − ( L − 1 ) / 2 ( L − 1 ) / 2 ∑ k = − ( K − 1 ) / 2 ( K − 1 ) / 2 | z 1 ( m + l , n + k ) | 2 ∑ l = − ( L − 1 ) / 2 ( L − 1 ) / 2 ∑ k = − ( K − 1 ) / 2 ( K − 1 ) / 2 | z 2 ( m + l , n + k ) | 2
where L×K is the image patch size, φ(m, n) is the phase difference between two images and * denotes the complex conjugate. It can be imagined that the computation will be heavy due to the involvement of φ. In order to reduce the computation, an alternative solution can be carried out, making use of the power correlation coefficient ρ̂ [56]. The computation of ρ̂ is as follows:
ρ ^ ( m , n ) = | ∑ l = − ( L − 1 ) / 2 ( L − 1 ) / 2 ∑ k = − ( K − 1 ) / 2 ( K − 1 ) / 2 | z 1 ( m + l , n + k ) | 2 | z 2 ( m + l , n + k ) | 2 | ∑ l = − ( L − 1 ) / 2 ( L − 1 ) / 2 ∑ k = − ( K − 1 ) / 2 ( K − 1 ) / 2 | z 1 ( m + l , n + k ) | 4 ∑ l = − ( L − 1 ) / 2 ( L − 1 ) / 2 ∑ k = − ( K − 1 ) / 2 ( K − 1 ) / 2 | z 2 ( m + l , n + k ) | 4
The relation between γ̂ and ρ̂ is as follows:
γ ^ = { 2 ρ ^ − 1 ρ ^ > 0.5 0 ρ ^ ≤ 0.5
However, because of the influence of noise, correlation coefficient array maybe has minimax values and maximum correlation coefficient is out of searching window sometimes [73]. So a new correlation coefficient is defined to settle the problem, in which the key point is to form a section of 4×4 pixels. The equation of correlation coefficient is as follows:
γ ( m , n ) new = ∑ m = 1 M / 4 ∑ n = 1 N / 4 | ∑ m = 1 4 ∑ n = 1 4 Z 1 ( m , n ) Z 2 ∗ ( m , n ) | ∑ m = 1 4 ∑ n = 1 4 | Z 1 ( m , n ) | 2 ∑ m = 1 4 ∑ n = 1 4 | Z 2 ( m , n ) | 2 M 4 N 4

3.5. Determination of Window Size for Tie Point Matching

The window size is another important factor to be considered for tie point matching, which can greatly affect the result of tie point matching. Too large a window size would result in not only a dramatic calculation but also the appearance of pseudo-tie point. However, too small a window size would make the co-registration unreliable. A multi-window sized 21×21 pixels based on maximum cross-correlation is developed by Li and Fan [74]. The window size, 33×33 pixels, 9×9 pixels and 43×43 pixels have been adopted by Zebker et al. [55], Weydahl [53] and Liao [63]. However, no theory has been developed to guide the selection of window size.
Recently a method for automated determination of an optimum window size for tie point matching was presented [45], which aims to find a mathematical method to determinate an optimum window size. The method is based on the decomposition of auto-correlation of the SAR image with wavelet transform to select an optimum window size for tie point matching. In this method, the auto-correlation of the SAR images, which reflects the similarity between image pixels, is used as a basis for analysis. To measure such a similarity, auto-correlation coefficient is the most widely used measure. It is defined as:
R ( d ) = Cov ( d ) V
where R(d) is the auto-correlation coefficient of all pixels with a distance of d part; Cov(d) is the covariance of all pixels with a distance of d part, and V is the variance of all pixels. Respectively, they are as follows:
V = ∑ i = 1 m ∑ j = 1 n ( | Z i , j | − K ) 2 m n − 1
Cov ( d ) = ∑ i = 1 m ∑ j = 1 n [ ( | Z i , j | − K ) ( | Q i , j + d | − K ) ] + ∑ i = 1 n ∑ j = 1 m [ ( | Z i , j | − K ) ( | Q i + d , j | − K ) ] 2 m n − 1
and
K = ∑ i = 1 M ∑ j = 1 N | Z i , j | M N
where m = 0, 1, 2, Λ, M-1, n = 1, 2, Λ, N-1, and the size of the image is M×N; Zi,j is one pixel with a coordinate (i, j) (expressed in pixel) in the image, its grey value is |Zi,j|; Qi,j+d is another pixel with a distance d apart from the pixel Zi,j in column direction, |Qi,j+d| is its gray value; Qi+d,j is another pixel with a distance d apart from the pixel Zi,j in row direction, |Qi+d,j| is its gray value; K is the average grey value of all pixels. A pair of SAR images and its curve of auto-correlation coefficient are shown in Figures 9 and 10, respectively.
Then, the auto-correlation function is decomposed into waves of various frequencies by wavelet transform as shown in Figure 11. By a combined analysis of the variations of wave amplitudes with frequency, the optimum window size for tie point matching can be determined. When wavelet db1 is used to decompose the auto-correlation coefficient, the jumping point (on the lowest-frequency component), at which the percentage amplitude is smaller than 15% of the amplitude at d = 0 and the difference of amplitude in percentage is smaller than 10%, can be selected as the optimum window size. In this way, an optimum window size could be determined for any pair of SAR images.

3.6. Selection of Transformation Models

Another issue for image co-registration are the mathematical models used for transformations. The transformation models can also affect the reliability of co-registration. In InSAR, two complex images are used, which cover the same scene but are taken from slightly different positions. Geometrically, they are in different coordinate systems. Their orientations could also be quite different. Therefore, an operation is required to bring them down to an identical coordinate system so as to have an identical orientation. There are two possible solutions: either to bring both down to an absolute ground coordinate or to fit one image into the coordinate system of the other. Normally, co-registration refers to the latter. In order to bring one image to fit into the coordinate system of the other, the relationship between these two images needs to be established. In InSAR practice, some sorts of polynomials are used as approximate models for the transformation between them. The models include geometric transformation models, and radiometric transformation models for re-sampling.
In order to enhance the accuracy of the results, different transformation models could be adopted in different situations. In the case of a photograph, one frame of an image is taken at one time. For each frame, there are six orientation elements (three rotations and three translations). The elements are the same for all image pixels on the same frame. If the image coordinate system is that the origin being the center of the image and the x-axis being along the flight direction, an arbitrary Cartesian coordinate system could possibly be used in the image plane. It is as follows [75]:
x = A 1 X + A 2 Y + A 3 Z + A 4 A 9 X + A 10 Y + A 11 Z + 1
y = A 5 X + A 6 Y + A 7 Z + A 8 A 9 X + A 10 Y + A 11 Z + 1
where, A1∼A11 are the coefficients; x and y are the coordinates of an image point; X, Y and Z are the coordinates of its position on the ground.
In the case of scanner images, image pixels are acquired scan by scan. A frame of an image may consist of thousands of scan lines. In this case, there are six orientation elements and there is a need of at least three control points for the computation of coefficients for each scan line. This is not feasible. The problem could be simplified if the relationship among the scan lines is known. However, this is normally not the situation. In this case, a single polynomial could be used to approximate the geometric model transformation for all scan lines [75]:
x = a 0 + a 1 X + a 2 Y + a 3 X 2 + a 4 Y 2 + a 5 X Y + ⋯
y = b 0 + b 1 X + b 2 Y + b 3 X 2 + b 4 Y 2 + b 5 X Y + ⋯
where, a0∼a5 and b0∼b5 are the coefficients; x and y are the coordinates of an image point; X and Y are the coordinates of its position on the ground. Indeed, in this case, the errors caused by distortions are re-distributed rather than removed by the polynomial function.
The latter case is suitable for InSAR practice. There are several functions used as the polynomial. They are known as nearest-neighbor interpolation, bilinear interpolation, and bicubic interpolation [76]. The bicubic function is adopted as the polynomial for the transformation between the master and slave images [44, 45, 68]. In InSAR image co-registration, (x, y) is the pixel's coordinate in the slave image while (X, Y) is the pixel's coordinate in the master image.
Zou et al. presented a method for an accurate co-registration of spacebore repeat-pass InSAR image based on Matrix Transformation, which considers the offsets errors generated by sensor trajectory excursion in repeat-pass style [77]. The coordinate transformation between the master image and the slave image is implemented by a matrix in which both shifting and rotating procedures are carried out.

3.7. Measurement for Quality of an Interferogram

The quality of an interferogram is an indicator for the reliability of SAR image co-registration. So far, the root mean square error (RMSE) of the residuals at tie points after least-squares adjustment was adopted as a measure for the quality of an interferogram. It implies that the smaller the RMSE, the better the interferogram. However, it is not always the case. In some cases, when the RMS is small, the resultant interferogram is not so good [68].
From the literature, it can be found that visual appearance of the fringes is widely used to indicate the quality of interferograms [68,78], that is, clear fringes imply good quality of an interferogram. In other words, if the interferometric fringes can be observed clearly, the quality of interferogram must be good. On the other hand, the blurry fringes will be found in bad interferogram. Visual inspection is a kind of manual interpretation.
The coherence of interferometric SAR is another indicator for the quality of interferogram [16,56]. The coherence between two SAR images can be theoretically defined on a pixel basis as the correlation coefficient and varies in the range of 0 to 1. It significantly influences the accuracy of phase in interferogram. So the degree of coherence can be used as a quality measure for the interferogram. The stronger the coherence, the better the quality. The coherence of an image varies from area to area or even from pixel to pixel. Therefore coherence image is produced to show the quality variation of an interferogram over the entire area. Figure 12 is a pair of SAR image and its corresponding coherence image is shown in Figure 13 [79]. In this image, some areas are brighter while some areas are darker. The brighter areas correspond to good coherence and the dark areas correspond to bad coherence.
To make a coherence measure quantitative, histogram of the coherence image is used. More precisely, the mean of the distribution and the standard deviation (STD) from the mean are the two quantitative parameters for the coherence. However, it is found even though the coherence images are quite different and the distributions of the pixel number in different interval of coherence are quite different, the actual values for the means and STD can be quite similar [79]. Therefore, such statistical parameters are not sufficient to differentiate the quality of interferograms. Based on the above, Li et al. [79] has proposed a more reliable quantitative measure for the InSAR interferogram based on the “sum of phase differences” (SPD).
The SPD method is based on the characteristics of the interferogram. Because an interferogram is generated by means of the product of the master image and the complex conjugate of the registered (slave) image, it contains an interferometric phase. Different colour fringes express different phase values in the interferogram. In an ideal case, the phase value changes gradually from 0 to 2π radians. The value of the phase difference between two pixels that are in neighboring fringes should be infinitesimal. If the two pixels are in the same colour fringe, the phase difference between them should be zero. Therefore, the smaller the SPD, the better the quality of the interferogram.
The SPDlocal of one pixel is calculated as follows [79]:
SPD local ( x , y ) = 1 8 ∑ l = − 1 1 ∑ k = − 1 1 | φ flattened ( x , y ) − φ flattened ( x + l , y + k ) |
where x = 0,1,2, ⋯, p − 1; y = 0,1,2, ⋯, q − 1 and the interferogram image size is p × q.
The SPDwhole of the whole interferogram is obtained by summarizing the SPDlocal of all pixels as follows:
SPD whole = ∑ x = 0 p − 1 ∑ y = 0 q − 1 SPD local ( x , y )
The indicator is adopted to measure the quality of the interferogram by Zou [71,80].

4. Summary

This paper aims to provide researchers with more information about the achievements in the study of InSAR image co-registration. Co-registration is the first and one of the most important steps in the interferometric InSAR data processing. In this step, the distribution of tie points, the size of window for tie point matching, extraction of feature points and the measurement for the quality of an interferogram can affect the accuracy of the co-registration, which in turn determines the quality of any InSAR products, and thus the quality of the final InSAR product, e.g. DEM. Therefore, the enhanced techniques on these critical factors for improving the accuracy of InSAR image co-registration with tie points are reviewed in this paper. Specially, the most recent achievements on extraction of feature points, effects of interval of tie points, determination of an optimum window size and an indicator of an interferogram quality are emphatically reviewed.

Acknowledgments

This work was supported by the Shenzhen Bureau of Science Technology & Information under the Research Project “Selection of Optimum Window Size for InSAR Image Co-registration”. The authors greatly appreciate the valuable comments and suggestions received from the reviewers of this paper.

References

  1. Zebker, H.A.; Rosen, P.A.; Goldstein, R.M.; Gabriel, A.; Werner, C.L. On the derivation of coseismic displacement fields using differential radar interferometry: The Landers earthquake. J. Geophys. Res. 1994, 99, 19617–19634. [Google Scholar]
  2. Hanssen, R.F. Radar Interferometry Data Interpretation and Error Analysis; Kluwer Acadamic Publishers Press: Dordrecht, Netherlands, 2000; p. 308. [Google Scholar]
  3. Rosen, P.A.; Hensley, S.; Zebker, H.A.; Webb, F.H.; Fielding, E.J. Surface deformation and coherence measurements of Kilauea Volcano, Hawaii, from SIR-C radar interferometry. J. Geophys. Res. 1996, 101, 23109–23125. [Google Scholar]
  4. Tobita, M.; Fujiwara, S.; Ozawa, S.; Rosen, P.A. Deformation of the 1995 North Sakhalin earthquake detected by JERS-1/SAR interferometry. Earth Planet. Space 1998, 50, 313–335. [Google Scholar]
  5. Liu, G.X. Mapping of earth deformations with satellite SAR interferometry: a study of its accuracy and reliability performances. Ph.D Dissertation, The Hong Kong Polytechnic University, Hong Kong, P.R. China, 2003. [Google Scholar]
  6. Askne, J.; Hagberg, J.O. Potential of interferometric SAR for classification of land surfaces. Proceedings of the International Geoscience and Remote Sensing Symposium, Tokyo, Japan, August 18-21, 1993; 3, pp. 985–987.
  7. Varkekamp, C.; Hoekman, D.H. Segmentation of high-resolution InSAR data of a tropic forest using Fourier parameterized deformable models. Int. J. Remote Sens. 2001, 22, 2339–2350. [Google Scholar]
  8. Wegmuller, U.; Werner, C. Retrieval of vegetation parameters with SAR interferometry. IEEE Trans. Geosci. Remote Sens. 1997, 35, 18–24. [Google Scholar]
  9. Weydahl, D.J. Analysis of ERS SAR coherence images acquired over vegetated areas and urban features. Int. J. Remote Sen. 2001, 22, 2811–2830. [Google Scholar]
  10. Mattar, K.E.; Vachon, P.W.; Geudtner, D.; Gray, A.L.; Cumming, I.G.; Brubman, M. Validation of Alpine glacier velocity measurements using ERS tandem-mission SAR data. IEEE Trans. Geosci. Remote Sens. 1998, 36, 974–984. [Google Scholar]
  11. Joughin, I.R.; Kwok, R.; Fahnestock, A. Interferometric estimation of three-dimensional ice-flow using ascending and descending passes. IEEE Trans. Geosci. Remote Sens. 1998, 36, 25–37. [Google Scholar]
  12. Wangensteen, B.; Weydahl, D.J.; Hagen, J.O. Mapping glacier velocities at Spitsbergen using ERS tandem SAR data. Proceedings of International Geoscience and Remote Sensing Symposium, Hamburg, Germany, June 28 – July 2, 1998; 4, pp. 1954–1956.
  13. Ludwig, R.; Lampart, G.; Mauser, W. The determination of hydrological model parameters from airborne interferometric SAR-data. Int. J. Remote Sens. 1998, 2609–2611. [Google Scholar]
  14. Madsen, S.N.; Zebker, H.A.; Martin, J. Topographic mapping using radar interferometry: processing techniques. IEEE Trans. Geosci. Remote Sens. 1993, 31, 246–256. [Google Scholar]
  15. Shi, J.C.; Dozier, J.; Rott, H. Snow mapping in Alpine regions with synthetic aperture radar. IEEE Trans. Geosci. Remote Sens. 1994, 32, 152–158. [Google Scholar]
  16. Gens, R. Quality assessment of SAR interferometric data. Ph.D. Dissertation, Hanover University, Hanover, Germany, 1998. [Google Scholar]
  17. Hilland, J.E.; Stuhr, F.V.; Anthony, F.; Imel, D.; Shen, Y.; Jordan, R.; Caro, E.R. Future NASA spaceborne SAR missions. IEEE AES Syst. Mag. 1998, 11, 9–16. [Google Scholar]
  18. Burgett, S. Target location in WGS-84 coordinates using synthetic aperture radar images. Proceedings of the Institute of Navigation's 49th Annual Meeting; 1993; pp. 57–65. [Google Scholar]
  19. Guarino, R.; Ibsen, P. Integrated GPS/INS/SAR/GMTI radar precision targeting flight test results. Proceedings of the Institute of Navigation's 51th Annual Meeting; 1995; pp. 373–379. [Google Scholar]
  20. Li, J.; Zelnio, E. Target detection with Synthetic Aperture Radar. IEEE Trans. Aeros. Electron. Syst. 1996, 32, 613–627. [Google Scholar]
  21. Fotnaro, G. Trajectory deviations in airborne SAR: analysis and compensation. IEEE Trans. Aeros. Electron. Syst. 1996, 35, 997–1009. [Google Scholar]
  22. Zyl, J.J.V. The shuttle radar topography mission (SRTM): a breakthrough in remote sensing of topography. Acta Astron. 2001, 48, 559–565. [Google Scholar]
  23. Miliaresis, G.C.; Paraschou, C.V.E. Vertical accuracy of the SRTM DTED level 1 of Crete. Int. J. Appl. Earth Observ. Geoinf. 2005, 7, 49–59. [Google Scholar]
  24. Hofton, M.A.; Dubayah, R.; Blair, J.B.; Rabine, D. Validation of SRTM Elevations over vegetated and non-vegetated terrain using medium footprint Lidar. Photogramm. Eng. Remote Sens. 2006, 72, 279–285. [Google Scholar]
  25. Walker, W.S.; Kellndorfer, J.M.; Pierce, L.E. Quality assessment of SRTM C-and X-band interferometric data: implications for the retrieval of vegetation canopy height. Remote Sens. Environ. 2007, 106, 428–448. [Google Scholar]
  26. Renno, C.D.; Nobre, A.D.; Cuartas, L.A.; Soares, J.V.; Hodnett, M.G.; Tomasella, J.; Waterloo, M.J. HAND, a new terrain descriptor using SRTM-DEM: mapping terra-firme rainforest environments in Amazonia. Remote Sens. Environ. 2008, 112, 3469–3481. [Google Scholar]
  27. Berthier, E.; Arnaud, Y.; Vincent, C.; Remy, F. Biases of SRTM in high-mountain areas: implications for the monitoring of glacier volume changes. Geophys. Res. Lett. 2006, 33, L08502. [Google Scholar]
  28. Kaab, A. Remote sensing of permafrost-related problems and hazards. Permafr. Periglac. Proc. 2008, 19, 107–136. [Google Scholar]
  29. Rodriguez, E.; Morris, C.S.; Belz, J.E. A global assessment of the SRTM performance. Photogramm. Eng. Remote Sens. 2006, 72, 249–260. [Google Scholar]
  30. Braun, A.; Fotopoulos, G. Assessment of SRTM, ICESat, and survey control monument elevations in Canada. Photogramm. Eng. Remote Sens. 2007, 73, 1333–1342. [Google Scholar]
  31. Hubbard, B.E.; Sheridan, M.F.; Carrasco-Nunez, G.; Diaz-Castellon, R.; Rodriguez, S.R. Comparative lahar hazard mapping at volcan Citlaltepetl, Mexico using SRTM, ASTER and DTED-1 digital topographic data. J. Volcanol. Geothermal. Res. 2007, 160, 99–124. [Google Scholar]
  32. Schumann, G.; Matgen, P.; Cutler, M.E.J.; Black, A.; Hoffmann, L.; Pfister, L. Comparison of remotely sensed water stages from LiDAR, topographic contours and SRTM. ISPRS J. Photogramm. Remote Sens. 2008, 63, 283–296. [Google Scholar]
  33. Carabajal, C.C.; Harding, D.J. SRTM C-band and ICESat laser altimetry elevation comparisons as a function of tree cover and relief. Photogramm. Eng. Remote Sens. 2006, 72, 287–298. [Google Scholar]
  34. Guth, P.L. Geomorphometry from SRTM: Comparison to NED. Photogramm. Eng. Remote Sens. 2006, 72, 269–277. [Google Scholar]
  35. Van, Niel T.G.; McVicar, T.R.; Li, L.T.; Gallant, J.C.; Yang, Q.K. The impact of misregistration on SRTM and DEM image differences. Remote Sens. Environ. 2008, 112, 2430–2442. [Google Scholar]
  36. Kiel, B.; Alsdorf, D.; LeFavour, G. Capability of SRTM C- and X-band DEM data to measure water elevations in Ohio and the Amazon. Photogramm. Eng. Remote Sens. 2006, 72, 313–320. [Google Scholar]
  37. Gorokhovich, Y.; Voustianiouk, A. Accuracy assessment of the processed SRTM-based elevation data by CGIAR using field data from USA and Thailand and its relation to the terrain characteristics. Remote Sens. Environ. 2006, 104, 409–415. [Google Scholar]
  38. Goncalves, J.A.; Korgado, A.M. Use of the SRTM DEM as a geo-referencing tool by elevation matching. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2008, XXXVII, 879–883. [Google Scholar]
  39. Galiatsatos, N.; Donoghue, D.N.M.; Philip, G. High resolution elevation data derived from stereoscopic CORONA imagery with minimal ground control: an approach using Ikonos and SRTM data. Photogramm. Eng. Remote Sens. 2007, 74, 1093–1106. [Google Scholar]
  40. Huggel, C.; Schneider, D.; Miranda, P.J.; Grannados, H.D.; Kaab, A. Evaluation of ASTER and SRTM DEM data for lahar modeling: a case study on lahars from Popocatepetl volcano, Mexico. J. Volcanol. Geothermal. Res. 2008, 170, 99–110. [Google Scholar]
  41. Menze, B.H.; Ur, J.A.; Sherrat, A.G. Detection of ancient settlement mounds: archaeological survey based on the SRTM terrain model. Photogramm. Eng. Remote Sens. 2006, 72, 321–327. [Google Scholar]
  42. Luedeling, E.; Siebert, S.; Buerkert, A. Filling the voids in the SRTM elevation model – a TIN-based delta surface approach. ISPRS J. Photogramm. Remote Sens. 2007, 62, 283–294. [Google Scholar]
  43. Schetselaar, E.M.; Tiainen, M.; Woldai, T. Integrated geological interpretation of remotely sensed data to support geological mapping in Mozambique. Geol. Surv. Finland 2008, 48, 35–63. [Google Scholar]
  44. Zou, W.; Li, Z.; Ding, X. Effects of the intervals of tiepoints used in co-registration on the accuracy of digital elevation models (DEM) generated by InSAR. Photogramm. Rec. 2006, 21, 232–254. [Google Scholar]
  45. Zou, W.; Li, Z.; Ding, X. Determination of optimum window size for SAR image co-registration with decomposition of auto-correlation. Photogramm. Rec. 2007, 22, 237–255. [Google Scholar]
  46. Masoud, A.; Raghavan, V.; Masumoto, S.; Shiono, K. Repeat-pass JERS-1 INSAR for DEM Generation in Safaga Area, Red Sea Coast of Egypt. Geoinformatics 2002, 13, 84–85. [Google Scholar]
  47. Zhou, C.; Ge, L.; E, D.; Chang, H. A case study of using external DEM in InSAR DEM generation. Geo-Spat. Inf. Sci. 2005, 8, 14–18. [Google Scholar]
  48. Bamler, R.; Eineder, M.; Kampes, B.; Runge, H.; Adam, N. SRTM and beyond: current situation and new developments in spaceborne InSAR. Proceedings of the ISPRS Workshop on High Resolution Mapping from Space, Istanbul, Turkey, October 6-8, 2003; University Hannover Germany, International Society for Photogrammetry and Remote Sensing.
  49. Schulz-stellenfleth, J.; Horstmann, J.; Lehner, S.; Rosenthal, W. Sea surface imaging with an across-track interferometric synthetic aperture radar: The SINEWAVE experiment. IEEE Trans. Geosci. Remote Sens. 2001, 39, 2017–2028. [Google Scholar]
  50. Lehner, S.; Gunther, H.; Horstmann, J.; Bao, M.; Schulz-Stellenfleth, J. Joint along-across track interferometry of ocean waves. Geoscience and Remote Sensing Symposium; 2001; 1, pp. 581–583. [Google Scholar]
  51. Gray, A.L.; Van der Kooij, M.W.A.; Mattar, K.E.; Farris-Manning, P.J. Progress in the development of the CCRS along-track interferometer. Geoscience and Remote Sensing Symposium; 1994; 4, pp. 2285–2287. [Google Scholar]
  52. Eberhard, G.; Hartmut, R. Tight formation flying for an along-track SAR interferometer. Acta Astron. 2004, 55, 473–485. [Google Scholar]
  53. Weydahl, D.J. Change detection in SAR images. In IEEE IGRSS; Espoo, Finland, June 3-6 1991; pp. 1421–1424. [Google Scholar]
  54. Lin, Q. New Approaches in Interferometric SAR Data Processing. IEEE Trans. Geosci. Remote Sens. 1992, 30, 560–567. [Google Scholar]
  55. Zebker, H.A.; Werner, C.L.; Rosen, P.A.; Hensley, S. Accuracy of topographic maps derived from ERS-1 interferometric radar. IEEE Trans. Geosci. Remote Sens. 1994, 32, 823–836. [Google Scholar]
  56. Guarnieri, A.M.; Prati, C. SAR Interferometry: A “Quick and Dirty” Coherence Estimatior for DATA Browsing. IEEE Trans. Geosci. Remote Sens. 1997, 35, 660–669. [Google Scholar]
  57. Bamler, R.; Just, D. phase statistics and decorrelation in SAR interferograms. IEEE International Geoscience and Remote Sensing Symposium, Tokoy, Japan, Aug 18-21, 1993; pp. 980–984.
  58. Lim, I.; Yeo, T.S.; Ng, C.S.; Lu, Y.H.; Zhang, B.C. Phase noise filter for interferometric SAR. IEEE Int. Geosci. Remote Sens. 1997, 1, 445–447. [Google Scholar]
  59. Ainsworth, T.L.; Lee, J.S. A joint adaptive interferometric phase unwrapping and filtering algorithm. IEEE IGRSS. 1998, 1, 71–73. [Google Scholar]
  60. Bamler, R.; Adam, N.; Davidson, G.W.; Just, D. Noise-induced slope distortion in 2-D phase unwrapping by linear estimators with application to SAR interferometry. IEEE Trans. Geosci. Remote Sens. 1998, 36, 913–921. [Google Scholar]
  61. Bao, M.; Kwoh, L. K.; Singh, K. An improved least-squares method for INSAR phase unwrapping. IEEE IGRSS. 1998, 1, 62–64. [Google Scholar]
  62. Wivell, C.E.; Steinwand, D.R.; Kelly, G.G.; Meyer, D.J. Evaluation of terrain models for the geocoding and terrain correction of synthetic aperture radar (SAR) image. IEEE Trans. Geosci. Remote Sens. 1992, 30, 1137–1143. [Google Scholar]
  63. Liao, M.; Lin, H.; Zhang, Z.X. Automatic registration of InSAR data based on Least-square matching and multi-step strategy. Photogramm. Eng. Remote Sens. 2004, 70, 1139–1144. [Google Scholar]
  64. Abdelfattah, R.; Nicolas, J.-M. Sub-pixelic Image Registration for SAR Interferometry Coherence Optimization. Proceedings of the XXth ISPRS Congress, Istanbul, Turkey, July 12-23, 2004; XXXV, p. 3.
  65. Wan, F.G. Digital Processing of Remotely Sensed Images; Huazhong Polytechnical University, 1990. [Google Scholar]
  66. Yang, Q. Y.; Wang, C. Registration of INSAR complex images and interferogram enhancement. J. Remote Sens. 1999, 3, 122–131. (In Chinese) [Google Scholar]
  67. Liu, G.X.; Ding, X.L.; Li, Z.L.; Chen, Y.Q.; Zhang, G.B. Co-registration of satellite SAR complex images. Acta Geod. Cartogr. Sinica. 2001, 30, 60–66. (in Chinese). [Google Scholar]
  68. Zou, W.; Li, Z.; Ding, X.; Chen, Y.; Liu, G. Effects of the interval of tiepoints on the reliability of SAR image co-registration. International Archives of Photogrammetry, Remote Sensing and Spatial Information Sciences, Xi'an, China, Aug 20-23, 2002; XXXIV, pp. 639–644.
  69. Urban, M. Harris Interest Operator. http://cmp.felk.cvut.cz/cmp/courses/dzo/resources/lecture_harris_urban.pdf, 2003. (accessed 2 January 2008).
  70. Moigne, J.L.; Campbell, W.J.; Cromp, R.F. An Automated Parallel Image Registration Technique Based on the Correlation of Wavelet Features. IEEE Trans. Geosci. Remote Sens. 2002, 40, 1849–1864. [Google Scholar]
  71. Zou, W.B. Improving the Accuracy of Image Co-registration in InSAR. Ph.D. Dissertation, Hong Kong Polytechnic University, 2005. [Google Scholar]
  72. Franceschetti, G.; Lanari, R. Synthetic Aperture Radar Processing; CRC Press: New York, NY, USA, 1999. [Google Scholar]
  73. Shi, X.; Zhang, Y.; Jiang, J. InSAR image registration using modified correlation coefficient algorithm. ISAPE '06 7th International Symposium on Antennas, Propagation & EM Theory; 2006; pp. 1–4. [Google Scholar]
  74. Li, Z.; Fan, X.T. Registration between remote sensing images based on multi-window cross-correlation. J. Image Graph. A 2001, 6, 129–132. (in Chinese). [Google Scholar]
  75. Li, Z. L. Orthoimage Generation and Measurement from Single Images. In Geographical data acquisition; Chen, Y.Q., Lee, Y.CH., Eds.; Springer: Wien, NewYork, 2001; pp. 173–183. [Google Scholar]
  76. Wolf, P.R.; Dewitt, B.A. Elements of Photogrammetry with Applications in GIS, 3rd Ed. ed; McGraw-Hill Science Engineering: New York, NY, USA, 2000. [Google Scholar]
  77. Zou, B.; Hao, L.; Bao, X. An accurate co-registration method of spaceborne repeat-pass InSAR based on matrix transformation. IEEE IGRSS 2005, 7, 4560–4563. [Google Scholar]
  78. Eldhuset, K.; Aanvik, F.; Aksnes, K. First results from ERS tandem InSAR processing on Svalbard. http://www.geo.unizh.ch/rsl/fringe96/papers/eldhuset-et-al/, 1996. (accessed 8 January 2008).
  79. Li, Z.; Zou, W.B.; Ding, X.L.; Chen, Y.Q.; Liu, G.X. A quantitative measure for the quality of InSAR interferogram based on phase differences. Photogramm. Eng. Remote Sens. 2004, 70, 1131–1137. [Google Scholar]
  80. Li, Z.W.; Ding, X.L.; Huang, C.; Zheng, D.W.; Zou, W.B. Filtering method for SAR interferograms with strong noise. Int. J. Remote Sens. 2006, 27, 2991–3000. [Google Scholar]
Sensors 09 01259f1 1024
Figure 1. Imaging geometry of SAR interferometry.
Figure 1. Imaging geometry of SAR interferometry.
Sensors 09 01259f1
Sensors 09 01259f2 1024
Figure 2. Imaging geometry of repeat-pass SAR interferometry.
Figure 2. Imaging geometry of repeat-pass SAR interferometry.
Sensors 09 01259f2
Sensors 09 01259f3 1024
Figure 3. Imaging geometry of across-track SAR interferometry.
Figure 3. Imaging geometry of across-track SAR interferometry.
Sensors 09 01259f3
Sensors 09 01259f4 1024
Figure 4. Imaging geometry of along-track SAR interferometry.
Figure 4. Imaging geometry of along-track SAR interferometry.
Sensors 09 01259f4
Sensors 09 01259f5 1024
Figure 5. SAR image processing procedures.
Figure 5. SAR image processing procedures.
Sensors 09 01259f5
Sensors 09 01259f6 1024
Figure 6. Procedures of InSAR image co-registration.
Figure 6. Procedures of InSAR image co-registration.
Sensors 09 01259f6
Sensors 09 01259f7 1024
Figure 7. A pair of SAR images (a) and (b) and its co-registered slave image. (a) a master image; (b) a slave image; (c) a co-registered slave image.
Figure 7. A pair of SAR images (a) and (b) and its co-registered slave image. (a) a master image; (b) a slave image; (c) a co-registered slave image.
Sensors 09 01259f7
Sensors 09 01259f8a 1024Sensors 09 01259f8b 1024
Figure 8. Feature points extracted by different ways. (a) feature points extracted by Interest Operators; (b) feature points extracted by wavelet transform.
Figure 8. Feature points extracted by different ways. (a) feature points extracted by Interest Operators; (b) feature points extracted by wavelet transform.
Sensors 09 01259f8aSensors 09 01259f8b
Sensors 09 01259f9 1024
Figure 9. A pair of SAR images (Tai Lam in Hong Kong). (a) the master image; (b) the slave image.
Figure 9. A pair of SAR images (Tai Lam in Hong Kong). (a) the master image; (b) the slave image.
Sensors 09 01259f9
Sensors 09 01259f10 1024
Figure 10. Auto-correlation coefficient of the Tai Lam image.
Figure 10. Auto-correlation coefficient of the Tai Lam image.
Sensors 09 01259f10
Sensors 09 01259f11 1024
Figure 11. Components of the auto-correlation of the Tai Lam image decomposed by wave db1. d1 ∼ d4 are the high-frequency components and a4 is the lowest-frequency component.
Figure 11. Components of the auto-correlation of the Tai Lam image decomposed by wave db1. d1 ∼ d4 are the high-frequency components and a4 is the lowest-frequency component.
Sensors 09 01259f11
Sensors 09 01259f12 1024
Figure 12. A pair of SAR images. (a) the master image; (b) the slave image.
Figure 12. A pair of SAR images. (a) the master image; (b) the slave image.
Sensors 09 01259f12
Sensors 09 01259f13 1024
Figure 13. The corresponding coherence image.
Figure 13. The corresponding coherence image.
Sensors 09 01259f13

Share and Cite

MDPI and ACS Style

Zou, W.; Li, Y.; Li, Z.; Ding, X. Improvement of the Accuracy of InSAR Image Co-Registration Based On Tie Points – A Review. Sensors 2009, 9, 1259-1281. https://doi.org/10.3390/s90201259

AMA Style

Zou W, Li Y, Li Z, Ding X. Improvement of the Accuracy of InSAR Image Co-Registration Based On Tie Points – A Review. Sensors. 2009; 9(2):1259-1281. https://doi.org/10.3390/s90201259

Chicago/Turabian Style

Zou, Weibao, Yan Li, Zhilin Li, and Xiaoli Ding. 2009. "Improvement of the Accuracy of InSAR Image Co-Registration Based On Tie Points – A Review" Sensors 9, no. 2: 1259-1281. https://doi.org/10.3390/s90201259

APA Style

Zou, W., Li, Y., Li, Z., & Ding, X. (2009). Improvement of the Accuracy of InSAR Image Co-Registration Based On Tie Points – A Review. Sensors, 9(2), 1259-1281. https://doi.org/10.3390/s90201259

Article Metrics

Back to TopTop