## 1.

## Introduction

With the successful integration of synthetic aperture radar imaging principle and interferometric measuring technique, synthetic aperture radar interferometry^{1}^{,}^{2} (InSAR) can be used to precisely measure the three-dimensional position and subtle change of a certain spot on the ground. In the past decade, it has been developed with great breakthroughs and has become an important branch in the field of radar remote sensing. Currently, its application has been extended to many fields, such as disaster monitoring, resources detecting, etc. The potential for its application is great.

The single look complex image used in the InSAR data processing includes four parts: two real parts and two imaginary parts. In the current methods, all these parts will be used to precisely register the complex image pair and generate interferometric phase images, which, however, are loaded with high-level decorrelation noise. The noise has a great impact on the phase unwrapping and the recovery of the digital elevation model (DEM) with high precision. This has become the bottleneck in InSAR data processing.^{3} Many scholars present filtering methods to effectively reduce the noise.^{4}5.^{–}^{6} Based on the contoured correlation fringe method used in electronic speckle pattern interferometry (ESPI),^{7} the author further proposes a new InSAR data processing method, which is named contoured correlation interferometry (CCI) and includes an image pair precise registration method^{8} and generation of interferometric phase images.^{9} In this method, the interferometric phase image, can be obtained by any three of the four real and imaginary parts of the two InSAR complex images with almost no phase noises and blurring. The principle of this new method has been published in *Applied Physics Letters*.^{9} This paper will introduce CCI with great details and do some makeups and extensions, so a complete InSAR data processing method is formed. The CCI formula will be induced and discussed from different aspects, the algorithms and procedures will be illustrated systematically, the algorithm of the fringe orientation map will be improved, the selection of the contoured window’s parameter will be discussed, and the results from different methods will be compared. Also some new data processing results about Etna volcano are given.

## 2.

## Defining of the Contoured Window of the InSAR Phase Image

The interferometric fringe pattern or the phase image has the following features: the fringe field is formed by the flow field in the gray level orientation, with the largest phase gradient in the fringe normal orientation and the smallest change of the gray level (phase) in the fringe tangent orientation, namely the phase being a constant. The fringe contour, on which the phase is keeping a constant, is the contoured gray level line in the fringe pattern. The various low-pass filtering on the fringe contour will effectively filter the noise without the phase being damaged. The authors have proposed a series of fringe pattern processing methods employing the fringe orientation information, such as spin filter,^{10} contoured-window filter,^{11} etc., which have been applied successfully in the processing of various interferometric fringe patterns, including holography, moiré patterns, ESPI fringe patterns, etc., and the latest InSAR phase images.^{12} The particularities of the fringe contour make many novel fringe processing methods possible.^{7}8.^{–}^{9}

As one of the key steps, the generation of the contoured window will be introduced first. In this paper, the contoured window is secured based on the tracking along the fringe orientation, and therefore the fringe orientation precision will determine the precision of the contoured window. Here we propose to replace the traditional plane-fit method^{8} with gradient method to reach an orientation result with better adaptation and greater accuracy.

## 2.1.

### Calculation of Fringe Orientation Map

In Refs. 7 through 9, the plane-fit method is employed to calculate the fringe orientation, but the plane-fit method is susceptible to the calculation window and achieves the best result only when the size of the calculation window approaches half of the local fringe’s width. Here we propose that the gradient method^{13} often be used in the processing of fringe orientation maps to calculate the InSAR fringe orientation. In an ideal condition, the fringe tangential orientation can be derived from

## (1)

$$\theta (x,y)=t{g}^{-1}[\frac{\partial I(x,y)}{\partial x}/\frac{\partial I(x,y)}{\partial y}],$$## (2)

$$\theta =\frac{1}{2}\mathrm{arc}tg\frac{\sum \mathrm{sin}\text{\hspace{0.17em}}2\theta}{\sum \mathrm{cos}\text{\hspace{0.17em}}2\theta}=\frac{1}{2}\mathrm{arc}tg\frac{\sum (2\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}\theta \text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\theta )}{\sum ({\mathrm{cos}}^{2}\text{\hspace{0.17em}}\theta -{\mathrm{sin}}^{2}\text{\hspace{0.17em}}\theta )}.$$## (3)

$$\theta =\frac{1}{2}\mathrm{arc}tg\frac{\sum (2{I}_{x}{I}_{y})}{\sum ({I}_{y}^{2}-{I}_{x}^{2})}.$$## 2.2.

### Determination of the Fringe Contoured Window

The interferometric fringe has the character of obvious fringe orientation flow, so we can track it along the fringe tangential direction to obtain the fringe contour. Suppose the coordinate of the current pixel is ${P}_{0}({x}_{0},{y}_{0})$, and its fringe orientation is ${\theta}_{0}$. We track it along the positive and negative directions of the fringe respectively and obtain the pixels ${P}_{i}({x}_{i},{y}_{i})$ and ${P}_{-i}({x}_{-i},{y}_{-i})$ [as shown in Fig. 3(a)] by

## (5)

$$\{\begin{array}{l}{x}_{\pm i}={x}_{\pm (i-1)}\pm \mathrm{cos}\text{\hspace{0.17em}}{\theta}_{\pm (i-1)}\\ {y}_{\pm i}={y}_{\pm (i-1)}\pm \mathrm{sin}\text{\hspace{0.17em}}{\theta}_{\pm (i-1)}\end{array},$$## 2.3.

### Determination of the Size of Fringe-Contoured Window

In the InSAR phase image, the fringe density and the fringe orientation sometimes vary greatly.^{14} For the dense fringe area, smaller windows should be used to preserve the fringe details; for the sparse fringe area, a bigger window should be used to filter the noise better and get a reliable result. Therefore the window’s size has a great impact on the calculation.

The phase value of the InSAR phase image being relatively consistent, as shown in the left part of Fig. 4, we can binarize the phase image with the phase value level threshold $\pi $ to obtain the better binary fringe image as shown by the right part of Fig. 4. Then the fringe width of every pixel can be calculated based on the space between the fringes in the binary fringe image, as shown in Fig. 5. To reduce the impact of the noise, we can first do the periodic pivoting filter in small windows to the InSAR phase image. The gray level values in Fig. 5 correspond with the density of the corresponding fringe in Fig. 4, the sparser the fringe is, the higher the gray level value is, the bigger the filter window is, and vice versa.

## 3.

## Contoured Correlation Interferometry to Generate InSAR Phase Images

As illustrated above, the interferometric phase image generation method is important to the InSAR data processing. The author has proposed CCI^{6} to generate interferometric phase images of high qualities. In this section, we will illustrate CCI in great details. As a comparison, the widely applied complex conjugate multiplicative method will be first introduced.

## 3.1.

### Current Complex Conjugate Multiplicative Method

The precisely registered complex image pair ${V}_{1}(r,x)$ and ${V}_{2}(r,x)$ are represented by Eqs. (6) and (7):

## (6)

$${V}_{1}(r,x)={A}_{1}{e}^{i{\varphi}_{1}}={A}_{1}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}{\varphi}_{1}+i{A}_{1}\text{\hspace{0.17em}}\mathrm{sin}\text{\hspace{0.17em}}{\varphi}_{1}\phantom{\rule{0ex}{0ex}}={A}_{1}\text{\hspace{0.17em}}\mathrm{cos}({\varphi}_{1c}+{u}_{1})+i{A}_{1}\text{\hspace{0.17em}}\mathrm{sin}({\varphi}_{1c}+{u}_{1})\phantom{\rule{0ex}{0ex}}={a}_{1}+{b}_{1}i,$$## (7)

$${V}_{2}(r,x)={A}_{2}{e}^{i{\phi}_{2}}={A}_{2}\text{\hspace{0.17em}}\mathrm{cos}({\phi}_{2c}+{u}_{2})+i{A}_{2}\text{\hspace{0.17em}}\mathrm{sin}({\phi}_{2c}+{u}_{2})\phantom{\rule{0ex}{0ex}}={a}_{2}+{b}_{2}i,$$## (8)

$${\varphi}_{1}-{\varphi}_{2}=({\varphi}_{1c}-{\varphi}_{2c})+({u}_{1}-{u}_{2})=\mathrm{\Delta}{\varphi}_{c}+\mathrm{\Delta}u\phantom{\rule{0ex}{0ex}}=\mathrm{arc}\mathrm{tan}\left\{\frac{\mathrm{Im}[{V}_{1}(r,x)\xb7{V}_{2}^{*}(r,x)]}{\mathrm{Re}[{V}_{1}(r,x)\xb7{V}_{2}^{*}(r,x)]}\right\}\phantom{\rule{0ex}{0ex}}=\mathrm{arctan}\left[\frac{{a}_{2}{b}_{1}-{a}_{1}{b}_{2}}{{a}_{1}{a}_{2}+{b}_{1}{b}_{2}}\right].$$After being registered, the two complex images correspond with the echo signal from the same area. Accordingly, most phase noises are correlative,^{15} the position and density of the speckle in the two images are basically consistent. Therefore the difference $\mathrm{\Delta}u$ can filter most correlative noises. However, due to different look angles of the two images, registration errors, and the time difference of the earth surface, etc., there still exists the decorrelation noise $\mathrm{\Delta}u$, which is usually much lower than the correlation noise and whose density depends on the degree of decorrelation.

## 3.2.

### CCI Formulae Induction

The correlation formula involved in CCI takes various forms: direct correlation, standardized correlation, and standardized covariance correlation, etc. Now take direct correlation formula, for example.

The direct correlation formula is as follows:

where ${\langle \xb7\rangle}_{m\times n}$ means calculating the mean value of some variable within the $m\times n$ pixels.As for the two complex images, ${V}_{1}(r,x)$ and ${V}_{2}(r,x)$, we take their real parts ${a}_{1}$ and ${a}_{2}$, respectively, and have

## (10)

$${\langle {a}_{1}{a}_{2}\rangle}_{m\times n}={\langle {A}_{1}\text{\hspace{0.17em}}\mathrm{cos}({\varphi}_{1c}+{u}_{1})\xb7{A}_{2}\text{\hspace{0.17em}}\mathrm{cos}({\varphi}_{2c}+{u}_{2})\rangle}_{m\times n}\phantom{\rule{0ex}{0ex}}={\langle \frac{1}{2}{A}_{1}{A}_{2}\xb7[\mathrm{cos}({\varphi}_{1c}+{\varphi}_{2c}+{u}_{1}+{u}_{2})+\mathrm{cos}({\varphi}_{1c}-{\varphi}_{2c}+{u}_{1}-{u}_{2})]\rangle}_{m\times n},$$^{16}Eq. (11) can be obtained on the window with the size $m\times n$.

## (11)

$${\langle \mathrm{cos}\text{\hspace{0.17em}}{u}_{i}\rangle}_{m\times n}={\langle \mathrm{sin}\text{\hspace{0.17em}}{u}_{i}\rangle}_{m\times n}=0\phantom{\rule[-0.0ex]{2em}{0.0ex}}(i=1,2).$$## (12)

$${\langle \mathrm{cos}({\varphi}_{1c}+{\varphi}_{2c}+{u}_{1}+{u}_{2})\rangle}_{m\times n}={\langle \mathrm{sin}({\varphi}_{1c}+{\varphi}_{2c}+{u}_{1}+{u}_{2})\rangle}_{m\times n}=0.$$## (13)

$${\langle {a}_{1}{a}_{2}\rangle}_{m\times n}={\langle \frac{1}{2}{A}_{1}{A}_{2}\xb7\mathrm{cos}(\mathrm{\Delta}{\varphi}_{c}+\mathrm{\Delta}u)\rangle}_{m\times n}\phantom{\rule{0ex}{0ex}}=\frac{1}{2}{\langle {A}_{1}{A}_{2}\rangle}_{m\times n}\xb7{\langle \mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}u-\mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}\text{\hspace{0.17em}}\mathrm{sin}\mathrm{\Delta}u\rangle}_{m\times n}.$$## (14)

$${\langle \mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}u\rangle}_{m\times n}=0,\phantom{\rule[-0.0ex]{2em}{0.0ex}}{\langle \mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}u\rangle}_{m\times n}\ne 0.$$## (16)

$${\langle {a}_{1}{a}_{2}\rangle}_{m\times n}=\frac{1}{2}{\langle {A}_{1}{A}_{2}\rangle}_{m\times n}\xb7[\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}{\langle \mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}u\rangle}_{m\times n}-\mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}{\langle \mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}u\rangle}_{m\times n}]\phantom{\rule{0ex}{0ex}}=\frac{1}{2}{\langle {A}_{1}{A}_{2}\rangle}_{m\times n}\xb7\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}{\langle \mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}u\rangle}_{m\times n}.$$## (17)

$${\langle {a}_{1}{b}_{2}\rangle}_{m\times n}=-\frac{1}{2}{\langle {A}_{1}{A}_{2}\rangle}_{m\times n}\xb7\mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}{\langle \mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}u\rangle}_{m,n}.$$## (18)

$$\frac{{\langle -{a}_{1}{b}_{2}\rangle}_{m\times n}}{{\langle {a}_{1}{a}_{2}\rangle}_{m\times n}}=\frac{\mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}{\langle \mathrm{cos}\mathrm{\Delta}u\rangle}_{m\times n}}{\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}{\langle \mathrm{cos}\mathrm{\Delta}u\rangle}_{m\times n}}=\frac{\mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}}{\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}}.$$## (19)

$$\mathrm{\Delta}{\varphi}_{c}=\mathrm{arctan}\left(\frac{{\langle -{a}_{1}{b}_{2}\rangle}_{m\times n}}{{\langle {a}_{1}{a}_{2}\rangle}_{m\times n}}\right)=\mathrm{arctan}\left(\frac{\mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}}{\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}}\right).$$The above conclusion is strictly right only on the fringe-contoured window that satisfies Eq. (15). With CCI, the decorrelation noise $\mathrm{\Delta}u$ on the contour can be theoretically removed completely with the phase value $\mathrm{\Delta}{\varphi}_{c}$ undamaged. Therefore the best window for the correlation calculation is the fringe-contoured window, through which we can get an interferometric phase image with best quality and highest precision. Next we will point out that when the condition of contour can’t be satisfied, we can use the rectangular window. When $\mathrm{\Delta}{\varphi}_{c}$ does not vary greatly within the rectangular window, we can still arrive at the similar result based on Eqs. (13) and (19).

The nature of the correlation coefficient is the resemblance degree of two functions. The physical significance is greater with CCI generating the phase image by measuring the resemblance degree of the phases of the two complex images through the correlation calculation. For the conventional complex conjugate multiplicative method, CCI is an alternative in the InSAR data processing and is supposed to be widely used.

## 3.3.

### Formulae Induction of the Correlation Interferometry with Rectangular Window

The best window for the generation of InSAR phase images with the correlation interferometry is the fringe-contoured window, which can be obtained from the phase image generated with the complex conjugate multiplicative method. But with the complex conjugate multiplicative method, all the four parts of the two complex images are needed. If we want to use only three of the four parts to finish the tasks of registering the InSAR complex images, generating the interferometric phase image and processing the data, we need to use the correlation interferometry with rectangular windows, which needs only three parts. We can first generate the suboptimal interferometric phase images, based on which we can obtain the fringe-contoured window. Finally we get the best interferometric phase image with CCI. Here the correlation interferometry with rectangular windows is discussed with three different sizes of windows:

1. The rectangular window with the correlation calculation window is $1\times 1$ ($m=1$; $n=1$) (single pixel calculation). When the calculation window is $1\times 1$, i.e., a single pixel, the data of the single pixel does not have any smoothing and suppressing effect on the random noise, and the conditions of Eqs. (11) and (12) are not satisfied. Then Eq. (10) can be changed into

We can have the expression of ${\langle {a}_{1}{b}_{2}\rangle}_{m\times n}$ in the same way. Thus## (20)

$${\langle {a}_{1}{a}_{2}\rangle}_{m\times n}=\frac{1}{2}{A}_{1}{A}_{2}\xb7[\mathrm{cos}({\varphi}_{1c}+{\varphi}_{2c}+{u}_{1}+{u}_{2})+\mathrm{cos}(\mathrm{\Delta}{\varphi}_{c}+\mathrm{\Delta}u)].$$Under this condition we cannot work out the phase $\mathrm{\Delta}{\varphi}_{c}$ with the correlation interferometry, and what’s more, the components of the original random speckle fields ${u}_{1}$, ${u}_{2}$, and $\mathrm{\Delta}u$ are preserved.## (21)

$$\mathrm{arctan}\left(\frac{{\langle -{a}_{1}{b}_{2}\rangle}_{m\times n}}{{\langle {a}_{1}{a}_{2}\rangle}_{m\times n}}\right)=\mathrm{arctan}\left[\frac{-\mathrm{sin}({\varphi}_{1c}+{\varphi}_{2c}+{u}_{1}+{u}_{2})+\mathrm{sin}(\mathrm{\Delta}{\varphi}_{c}+\mathrm{\Delta}u)}{\mathrm{cos}({\varphi}_{1c}+{\varphi}_{2c}+{u}_{1}+{u}_{2})+\mathrm{cos}(\mathrm{\Delta}{\varphi}_{c}+\mathrm{\Delta}u)}\right].$$2. The correlation calculation window is a small rectangular window (for example, $m=3$ or 5; $n=3$ or 5).

When the calculation window is small, the correlation calculation has a certain smoothing effect but still does not satisfy Eqs. (11) and (12).

The key induction condition being only approximately satisfied, there exists great error in the above calculation, resulting in the generated interferometric phase image with serious decorrelation noise.## (22)

$$\mathrm{\Delta}{\varphi}_{c}\approx \mathrm{arctan}\left(\frac{{\langle -{a}_{1}{b}_{2}\rangle}_{m\times n}}{{\langle {a}_{1}{a}_{2}\rangle}_{m\times n}}\right)\approx \mathrm{arctan}\left(\frac{\mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}}{\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}}\right).$$3. The correlation calculation window is large rectangular window (for example, $m\ge 7$, $n\ge 7$). When the rectangular window is large enough to satisfy Eqs. (11) and (12), the phase distribution in the rectangular window is not a constant, namely $\mathrm{\Delta}{\varphi}_{c}\ne \mathrm{const}$. Equation (10) is changed into

Equation (17) is changed into## (23)

$${\langle {a}_{1}{a}_{2}\rangle}_{m\times n}=\frac{1}{2}{\langle {A}_{1}{A}_{2}\rangle}_{m\times n}\xb7{\langle \mathrm{cos}(\mathrm{\Delta}{\varphi}_{c}+\mathrm{\Delta}u)\rangle}_{m\times n}\phantom{\rule{0ex}{0ex}}=\frac{1}{2}{\langle {A}_{1}{A}_{2}\rangle}_{m\times n}\xb7[{\langle \mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}\rangle}_{m\times n}{\langle \mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}u\rangle}_{m\times n}-{\langle \mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}\rangle}_{m\times n}{\langle \mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}u\rangle}_{m\times n}]\phantom{\rule{0ex}{0ex}}=\frac{1}{2}{\langle {A}_{1}{A}_{2}\rangle}_{m\times n}\xb7{\langle \mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}\rangle}_{m\times n}{\langle \mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}u\rangle}_{m\times n}.$$Therefore,## (24)

$${\langle {a}_{1}{b}_{2}\rangle}_{m\times n}=\frac{1}{2}{\langle {A}_{1}{A}_{2}\rangle}_{m\times n}\xb7{\langle \mathrm{sin}(\mathrm{\Delta}{\varphi}_{c}+\mathrm{\Delta}u)\rangle}_{m\times n}\phantom{\rule{0ex}{0ex}}=-\frac{1}{2}{\langle {A}_{1}{A}_{2}\rangle}_{m\times n}\xb7{\langle \mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}\rangle}_{m\times n}{\langle \mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}u\rangle}_{m\times n}.$$## (25)

$$\mathrm{arctan}\left(\frac{{\langle {a}_{1}{b}_{2}\rangle}_{m\times n}}{{\langle {a}_{1}{a}_{2}\rangle}_{m\times n}}\right)=\mathrm{arctan}\left[\frac{\mathrm{sin}(\mathrm{\Delta}{\varphi}_{c}+\mathrm{\Delta}u)}{\mathrm{cos}(\mathrm{\Delta}{\varphi}_{c}+\mathrm{\Delta}u)}\right]\phantom{\rule{0ex}{0ex}}=\mathrm{arctan}\left[\frac{{\langle \mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}\rangle}_{m\times n}{\langle \mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}u\rangle}_{m\times n}}{{\langle \mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}\rangle}_{m\times n}{\langle \mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}u\rangle}_{m\times n}}\right]\phantom{\rule{0ex}{0ex}}=\mathrm{arctan}\left(\frac{{\langle \mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}\rangle}_{m\times n}}{{\langle \mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}\rangle}_{m\times n}}\right)\approx \mathrm{\Delta}{\varphi}_{c}.$$In the above induction, due to the large calculation window, the random speckle fields ${u}_{1}$, ${u}_{2}$ and the decorrelation noise are all suppressed with the wanted phase $\mathrm{\Delta}{\varphi}_{c}$ being kept. But we have to use approximate expressions, for the rectangular window does not satisfy the condition that $\mathrm{\Delta}{\varphi}_{c}$ is a constant:

Accordingly, the correlation interferometric phase image generated in the large rectangular window is a mean phase image that has been averaged, has the blurring effect, and is approximate. Only when the correlation calculation window is the fringe-contoured window, the induction condition will be strictly satisfied, and we will obtain the phase field from Eq. (19).## (26)

$${\langle \mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}\rangle}_{m\times n}\approx \mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}\varphi ,\phantom{\rule[-0.0ex]{2em}{0.0ex}}{\langle \mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}\rangle}_{m\times n}\approx \mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}\varphi .$$

## 3.4.

### Principle Steps to Generate InSAR Interferometric Phase Image with CCI Method

To sum up, the principle steps to generate InSAR interferometric phase image with CCI:

1. Generate the interferometric phase image with phase noises by the four parts conventional complex conjugate multiplicative method, or generate the interferometric phase image with blurring effect by the three parts correlation interferometry of the rectangular window, from which we work out the fringe-orientation map and establish the fringe-contoured window.

2. For each pair of corresponding pixels in the two complex images, we take the current pixel as the center and establish a contoured window of size $m\times n$, in which we do the correlation calculation to the data of the two complex images with Eq. (9), and thus we get the correlation coefficient ${C}_{1}$.

3. In the same way, we do the correlation calculation to the real parts (imaginary parts) of the main image data and the imaginary parts (real parts) of the complementary image data with Eq. (12) to work out the correlation coefficient ${C}_{2}(-{C}_{2})$.

4. Work out the inverse tangent of the ratio between ${C}_{2}$ and ${C}_{1}$, and we will obtain the interferometric phase image.

## 4.

## Co-Registration Based on Three Parts of Two Complex Images and Contoured Windows for Synthetic Aperture Radar Interferometry

In the InSAR data processing, the precise registration technique of the two single look complex images is one of the key factors in improving the measurement accuracy. The registration error is one main source of the noises in the interferometric phase image. The improvement of the registration accuracy and reduction of registration error is of great significance in improving the quality of the interferometric phase image, reducing the difficulty of phase unwrapping, and improving the unwrapping precision.

With the illumination of the above CCI method to generate InSAR phase images, we propose the new InSAR registration criterion based on the correlation coefficient of three parts, according to which only any three parts of the complex images are needed, and the registration result is better than that using the conventional relevant coefficient registration criterion.

The formula induction of the new registration criterion $g$ is similar to that in CCI and thus it will only be briefly introduced. As for the detailed induction, please refer to Ref. 10. The correlation formulae used can be such mathematic expression as standardized correlation, standardized covariance correlation, etc. Here we take the standardized correlation as an example to illustrate the induction:

## (27)

$$C(r,x)=\frac{{\langle {f}_{1}\xb7{f}_{2}\rangle}_{m\times n}}{\sqrt{{\langle {f}_{1}^{2}\rangle}_{m\times n}\xb7{\langle {f}_{2}^{2}\rangle}_{m\times n}}}.$$## (28)

$${\langle {a}_{1}{a}_{2}\rangle}_{m\times n}={\langle {A}_{1}\text{\hspace{0.17em}}\mathrm{cos}({\varphi}_{1c}+{u}_{1})\xb7{A}_{2}\text{\hspace{0.17em}}\mathrm{cos}({\varphi}_{2c}+{u}_{2})\rangle}_{m\times n}\phantom{\rule{0ex}{0ex}}={\langle \frac{1}{2}{A}_{1}{A}_{2}\xb7[\mathrm{cos}({\varphi}_{1c}+{\varphi}_{2c}+{u}_{1}+{u}_{2})+\mathrm{cos}({\varphi}_{1c}-{\varphi}_{2c}+{u}_{1}-{u}_{2})]\rangle}_{m\times n}.$$## (29)

$${\langle {a}_{1}{a}_{2}\rangle}_{m\times n}={\langle {A}_{1}\text{\hspace{0.17em}}\mathrm{cos}({\varphi}_{1c}+{u}_{1})\xb7{A}_{2}\text{\hspace{0.17em}}\mathrm{cos}({\varphi}_{2c}+{u}_{2})\rangle}_{m\times n}\phantom{\rule{0ex}{0ex}}={\langle \frac{1}{2}{A}_{1}{A}_{2}\xb7[\mathrm{cos}({\varphi}_{1c}+{\varphi}_{2c}+{u}_{1}+{u}_{2})+\mathrm{cos}({\varphi}_{1c}-{\varphi}_{2c}+{u}_{1}-{u}_{2})]\rangle}_{m\times n}\phantom{\rule{0ex}{0ex}}=\frac{1}{2}{\langle {A}_{1}{A}_{2}\rangle}_{m\times n}\xb7{\langle 1+\mathrm{cos}(2{\varphi}_{1c}+2{u}_{1})\rangle}_{m\times n}$$## (30)

$${\langle {a}_{1}{a}_{2}\rangle}_{m\times n}=\frac{1}{2}{\langle {A}_{1}{A}_{2}\rangle}_{m\times n}\xb7{\langle 1+\mathrm{cos}(2{\varphi}_{2c}+2{u}_{2})\rangle}_{m\times n},$$## (31)

$${C}_{1}=\frac{{\langle {A}_{1}{A}_{2}\rangle}_{m\times n}\xb7{\langle \mathrm{cos}({\varphi}_{1c}+{\varphi}_{2c}+{u}_{1}+{u}_{2})+\mathrm{cos}({\varphi}_{1c}+{\varphi}_{2c}+{u}_{1}-{u}_{2})\rangle}_{m\times n}}{\sqrt{{\langle {A}_{1}^{2}\rangle}_{m\times n}\xb7{\langle 1+\mathrm{cos}(2{\varphi}_{1c}+2{u}_{1})\rangle}_{m\times n}\xb7{\langle {A}_{2}^{2}\rangle}_{m\times n}\xb7{\langle 1+\mathrm{cos}(2{\varphi}_{2c}+2{u}_{2})\rangle}_{m\times n}}}\phantom{\rule{0ex}{0ex}}=\frac{{\langle {A}_{1}{A}_{2}\rangle}_{m\times n}}{\sqrt{{\langle {A}_{1}^{2}\rangle}_{m\times n}\xb7{\langle {A}_{2}^{2}\rangle}_{m\times n}}}\xb7{\langle \mathrm{cos}(\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}+\mathrm{\Delta}u)\rangle}_{m\times n}.$$Based on ${C}_{1}$ and ${C}_{2}$, we define the following registration criterion:

## (32)

$$g=\sqrt{{C}_{1}^{2}+{C}_{2}^{2}}\phantom{\rule{0ex}{0ex}}=\frac{|{\langle {A}_{1}{A}_{2}\rangle}_{m\times n}|}{\sqrt{{\langle {A}_{1}^{2}\rangle}_{m\times n}\xb7{\langle {A}_{2}^{2}\rangle}_{m\times n}}}\xb7\sqrt{{\langle \mathrm{cos}(\mathrm{\Delta}{\varphi}_{c}+\mathrm{\Delta}u\rangle}_{m\times n}^{2}+{\langle \mathrm{sin}(\mathrm{\Delta}{\varphi}_{c}+\mathrm{\Delta}u\rangle}_{m\times n}^{2}}.$$$\mathrm{\Delta}{\varphi}_{c}$ is a constant within the contoured window. Put it into Eq. (32) and we have

## (33)

$${\langle \mathrm{cos}(\mathrm{\Delta}{\varphi}_{c}+\mathrm{\Delta}u)\rangle}_{m\times n}={\langle \mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}u-\mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}\mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}u\rangle}_{m\times n}\phantom{\rule{0ex}{0ex}}=\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}{\langle \mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}u\rangle}_{m\times n}$$## (34)

$${\langle \mathrm{sin}(\mathrm{\Delta}{\varphi}_{c}+\mathrm{\Delta}u)\rangle}_{m\times n}={\langle \mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}u-\mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}\mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}u\rangle}_{m\times n}\phantom{\rule{0ex}{0ex}}=\mathrm{sin}\text{\hspace{0.17em}}\mathrm{\Delta}{\varphi}_{c}{\langle \mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}u\rangle}_{m\times n}.$$## (35)

$$g=\sqrt{{C}_{1}^{2}+{C}_{2}^{2}}=\frac{|{\langle {A}_{1}{A}_{2}\rangle}_{m\times n}|}{|{\langle {A}_{1}^{2}\rangle}_{m\times n}\xb7{\langle {A}_{2}^{2}\rangle}_{m\times n}|}\xb7{\langle \mathrm{cos}\text{\hspace{0.17em}}\mathrm{\Delta}u\text{\hspace{0.17em}}\rangle}_{m\times n}.$$Figure 6 is the flowchart of the three parts registration method.

## 5.

## Experimental Results of CCI in InSAR Data Processing

To prove the effectiveness of CCI proposed in this paper, we would like to use CCI to process real InSAR data, which are from ERS1/2. Figure 7 is the result comparison of three parts registration method and the four parts correlation coefficient registration method. Compared in the figure are the correlation coefficients of the corresponding master image and the slave image after they have been registered with the two methods. The $x$-axis is the correlation coefficient, and the $y$-axis is the number of pixels. The figure shows that the correlation coefficients with the method proposed in this paper are all higher than those with the conventional method. This means that the method proposed in this paper is more accurate, though only three parts are used.

Figure 8(a) is the interferometric phase image generated with the conventional complex conjugate multiplicative method after the registration with three parts registration method; Fig. 8(b) is the phase image generated with CCI. Figure 8 shows that the three parts registration method is right and effective, and the phase image generated with CCI three parts is obviously better than that generated with the conventional four parts. With CCI, we can get the smooth phase image with almost no phase noises, while also the fringe phase is well kept without the blurring effect.

As illustrated previously, in the CCI formula induction, we can choose direct correlation, standardized correlation, and standardized covariance correlation, etc. The original data in Fig. 9 is the same as that in Fig. 8. Figure 9(a) is the interferometric phase image generated by calculation of each pixel (the window is $1\times 1$) with the covariance correlation formula CCI method, which is different from the direct correlation formula induced in this paper. From the figure, we can see that the phase can’t be obtained by the covariance correlation formula. Figure 9(b) is the result of direct correlation method proposed in this paper. Since we, in fact, did not do the correlative average calculation, the conditions of Eqs. (11) and (12) are not satisfied. Therefore, the original phase ${\varphi}_{1c}$, ${\varphi}_{2c}$, and the random noises ${u}_{1}$ and ${u}_{2}$ are all preserved in the generated image. The noise is larger than the decorrelation noise with the conventional complex conjugate multiplicative method.

Figure 9(c) and 9(d) are the interferometric phase images generated with the covariance correlation formula CCI method after being correlated in the rectangular windows $5\times 5$ and $9\times 9$, respectively. With the window increasing, the random phase noises decreases quickly. Figure 9 shows that after the correlation calculation in a window of certain size, the phase noises will be quickly removed, for the phase noises conform to the statistical rules in Eqs. (11), (12), and (14). Figure 9 also shows that the result is worse when the covariance correlation formula CCI method with which the mean value is deducted is adopted in a small window, compared with the CCI method proposed in this paper. Accordingly, the proposed direct correlation formula should be adopted.

The comparison also shows that we can obtain a better phase image with the rectangular window three parts generation method, but it is worse than that generated with the contoured window, for it has some blurring effect. Although only an approximate fringe pattern is generated with the rectangular window, this is one of the important steps in the three parts registration and the interferometric phase image generation, for it provides the fringe orientation with the three parts, making the need of the method for the three parts a closed loop and the method perfect. Without such a step, the fringe orientation would be obtained only through the fringe pattern generated with the conventional complex conjugate multiplicative method and thus the four parts would have to be used, but not the three parts in the last stage.

Figure 10 is the result comparison of the real airborne InSAR data processing with CCI and the conventional complex conjugate multiplicative method. Due to the small decorrelation noise with the airborne InSAR data, the complex conjugate multiplicative method can also arrive at a good result. But the noise in the phase image generated with CCI is smaller. Figure 11 is the result comparison of the CCI processing of different real and imaginary parts of the two complex images in different areas with the same set of data. Take any three parts from the four parts (${a}_{1}$, ${a}_{2}$, ${b}_{1}$, ${b}_{2}$), and we have four different combinations (${a}_{1}$, ${b}_{1}$, ${a}_{2}$), (${a}_{1}$, ${b}_{2}$, ${a}_{2}$), (${b}_{1}$, ${a}_{1}$, ${b}_{2}$), and (${b}_{1}$, ${a}_{2}$, ${b}_{2}$). From the image, we can see that the processing results of the four sets with CCI are almost the same, proving that the selection of parts with CCI is arbitrary.

Figure 12 is the processing result of the same area with serious decorrelation noise with the conventional complex conjugate multiplicative method (left half) and CCI (right half). It shows that for the area with lower correlation, it is difficult to process by the conventional method, but CCI can still result in good phase message. The theoretical analysis and many experimental results show that as long as the approximately right fringe-orientation map (namely, the fringe-contoured window), is obtained, we can to a great extent reduce or even remove the decorrelation effect by doing the correlation calculation with CCI on the approximate contour. For the area with serious decorrelation, we can still obtain the phase undisturbed and the phase noises removed.

Figure 13(a) and 13(b) are the results of InSAR data over Mount Etna volcano in Italy whose fringe pattern is with fast variations. Figures 12 and 13 show that CCI can still result in a good phase message in the area with poor correlation, which is difficult to process by the conventional method. Both theoretical analysis and many experimental results show that as long as the approximately right fringe-orientation map (namely, the fringe-contoured window) is obtained, we can to a great extent reduce or even remove the decorrelation effect by doing the correlation calculation with CCI on the approximate contoured windows. For the area with serious decorrelation, we can still obtain the undisturbed phase with the phase noises being removed.

Due to the slower change of the fringe orientation than the phase (namely, the fringe), the recovery of the approximate fringe orientation is much easier and more reliable than that of the right phase, especially for the area with obvious decorrelation. Therefore CCI can effectively improve the adaptation and the applicability of the InSAR data. Besides, the experiment shows that CCI is not quite sensitive to the precision of the contoured window; the contoured window with some error still being able to get better results than the conventional rectangular window.

## 6.

## Conclusion

In this paper, the proposed registration method and phase image-generation method based on three parts of the InSAR complex image pair are systematically illustrated and expanded. The adaptation and the precision of CCI method, the improved key steps, and the determination of the contoured window are also analyzed. The processing result of the real data has proved that the three parts registration method is more accurate than the conventional correlation coefficient-based registration method and is a brand-new concept of InSAR data-processing method. The phase image with the method proposed in this paper contains almost no phase noises, preserves the fringe at the same time, and is free of the blurring effect. As for the conventional method, while the phase noises are reduced, the signal is blurred. Although CCI is more complicated than the conventional method, the generated phase image is void of the phase noises, which saves the trouble of filtering in post-processing.

What’s more important, the proposed registration method and the interferometric phase image-generation method contribute to forming a complete InSAR data-processing method that involves only three parts. For the space-borne InSAR system, if the imaging is to be completed on the satellite, only any three of the four parts need to be transmitted back to the earth to finish the InSAR data processing, greatly reducing the transmission load, which is of great significance for the transmitting and processing of the space-borne InSAR data. One of the critical steps with CCI is the obtaining of the contoured windows, for which we will do some further research concerning how to determine the contoured windows with better self-adaptation, higher efficiency, and greater precision.

## Acknowledgments

This research is funded by the National Nature Science Foundation of China (Grant No. 40901215 & 11002156) and the Foundation of National University of Defense Technology.

## References

P. A. RosenS. HensleyI. R. Joughin, “Synthetic aperture radar interferometry,” Proc. IEEE 88(3), 333–382 (2000), http://dx.doi.org/10.1109/5.838084.IEEPAD0018-9219Google Scholar

R. BamlerP. Hartl, “Synthetic aperture radar interferometry,” Inverse Probl. 14(4), R1–R54 (1998), http://dx.doi.org/10.1088/0266-5611/14/4/001.INPEEY0266-5611Google Scholar

J.-S. Leeet al., “A new technique for noise filtering of SAR interferometric phase images,” IEEE Trans. Geosci. Rem. Sens. 36(5), 1456–1465 (1998), http://dx.doi.org/10.1109/36.701024.IGRSD20196-2892Google Scholar

H. LiG. Liao, “An estimation method for InSAR interferometric phase based on MMSE criterion,” IEEE Trans. Geosci. Rem. Sens. 48(3), 1457–1469 (2010), http://dx.doi.org/10.1109/TGRS.2009.2031100.IGRSD20196-2892Google Scholar

B. YongB. Mercer, “Interferometric SAR phase filtering in the wavelet domain using simultaneous detection and estimation,” IEEE Trans. Geosci. Rem. Sens. 49(4), 1396–1416 (2011), http://dx.doi.org/10.1109/TGRS.2010.2076286.IGRSD20196-2892Google Scholar

Q. Wanget al., “An efficient and adaptive approach for noise filtering of SAR interferometric phase images,” IEEE Geosci. Rem. Sens. Lett. 8(6), 1140–1144 (2011), http://dx.doi.org/10.1109/LGRS.2011.2158289.IGRSBY1545-598XGoogle Scholar

Q. Yuet al., “Single-phase-step method with contoured correlation fringe patterns for ESPI,” Opt. Express 12(20), 4980–4985 (2004), http://dx.doi.org/10.1364/OPEX.12.004980.OPEXFF1094-4087Google Scholar

Q. Yuet al., “Co-registration based on three parts of two complex images and contoured windows for synthetic aperture radar interferometry,” IEEE Geosci. Rem. Sens. Lett. 4(2), 288–292 (2007), http://dx.doi.org/10.1109/LGRS.2007.894146.IGRSBY1545-598XGoogle Scholar

Q. YuS. FuH. Mayer, “Generation of speckle-reduced phase images from three complex parts for SAR interferometry,” Appl. Phys. Lett. 88(11), 114106 (2006), http://dx.doi.org/10.1063/1.2185250.APPLAB0003-6951Google Scholar

Q. Yu, “Spin filtering process and automatic extraction of fringe center-lines from interferometric patterns,” Appl. Opt. 27(18), 3782–3784 (1988), http://dx.doi.org/10.1364/AO.27.003782.APOPAI0003-6935Google Scholar

Q. YuX. SunX. Liu, “Spin filtering with curve windows for interferometric fringes,” Appl. Opt. 41(14), 2650–2654 (2002), http://dx.doi.org/10.1364/AO.41.002650.APOPAI0003-6935Google Scholar

Q. Yuet al., “An adaptive contoured window filter for interferometric synthetic aperture radar,” IEEE Geosci. Rem. Sens. Lett. 4(1), 23–26 (2007), http://dx.doi.org/10.1109/LGRS.2006.883527.IGRSBY1545-598XGoogle Scholar

M. KassA. Witkin, “Analyzing oriented patterns,” Comp. Vis. Graph. Image Proc. 37(3), 362–385 (1987).CVGPDB0734-189Xhttp://dx.doi.org/10.1016/0734-189X(87)90043-0Google Scholar

S. Fuet al., “Directionally adaptive filter for synthetic aperture radar interferometric phase images,” IEEE Trans. Geosci. Rem. Sens. 51(1), 552–559 (2012), http://dx.doi.org/10.1109/TGRS.2012.2202911.IGRSD20196-2892Google Scholar

R. BamlerP. Hartl, “Synthetic aperture radar interferometry,” Inverse Probl. 14(4), R1–R54 (1998), http://dx.doi.org/10.1088/0266-5611/14/4/001.INPEEY0266-5611Google Scholar

J. W. Goodman, Laser Speckle and Related Phenomena, Springer-Verlag, Berlin (1975).Google Scholar

## Biography

**Jianhua Shi** is an associate professor at the National University of Defense Technology in China. She received her PhD degrees in optical engineering from the National University of Defense Technology in 2004. She is the author of about 20 journal papers and has written two book chapters. Her current research interests include interferogram processing and applied laser technology.

**Sihua Fu** received his BS, MS, and PhD degrees in aeronautical and astronautical science and technology from the National University of Defense Technology, P. R. China, in 1999, 2002, and 2006, respectively. He was with the National University of Defense Technology as an assistant professor in 2006. From 2010, he is an associate professor with the College of Opto-Electronic Science and Engineering. His research topics include image measurement technologies, InSAR data processing and optical interferometic fringes processing.