## Abstract

In the present work, the geometry and basic parameters of interferometric synthetic aperture radar (InSAR) geophysics system are addressed. Equations of pixel height and displacement evaluation are derived. Synthetic aperture radar (SAR) signal model based on linear frequency modulation (LFM) waveform and image reconstruction procedure are suggested. The concept of pseudo InSAR measurements, interferogram, and differential interferogram generation is considered. Interferogram and differential interferogram are generated based on a surface model and InSAR measurements. Results of numerical experiments are provided.

### Keywords

- InSAR
- geometry
- signal modeling
- SAR interferogram
- SAR differential interferograms

## 1. Introduction

Synthetic aperture radar (SAR) is a coherent microwave imaging instrument capable to provide for data all weather, day and night, guaranteeing global coverage surveillance. SAR interferometry is based on processing two or more complex valued SAR images obtained from different SAR positions [1, 2, 3, 4]. The InSAR is a system intends for geophysical measurements and evaluation of topography, slopes, surface deformations (volcanoes, earthquakes, ice fields), glacier studies, vegetation growth, etc. The estimation of topographic height with essential accuracy is performed by the interferometric distance difference measured based on two SAR echoes from the same surface. Changes in topography (displacement), precise to a fraction of a radar wavelength, can be evaluated by differential interferogram generated by three or more successive complex SAR images [5, 6]. Demonstration of time series InSAR processing in Beijing using a small stack of Gaofen-3 differential interferograms is discussed in [7].

A general overview of the InSAR principles and the recent development of the advanced multi-track InSAR combination methodologies, which allow to discriminate the 3-D components of deformation processes and to follow their temporal evolution, are presented in [8]. The combination of global navigation satellite system (GNSS) and InSAR for future Australian datums is discussed in [9].

A high-precision DEM extraction method based on InSAR data and quality assessment of InSAR DEMs is suggested in [10, 11]. InSAR digital surface model (DSM) and time series analysis based on C-band Sentinel-1 TOPS data are presented in [12, 13]. DEM registration, alignment, and evaluation for SAR interferometry, deformation monitoring by ground-based SAR interferometry (GB-InSAR), a field test in dam, and an improved approach to estimate large-gradient deformation using high-resolution TerraSAR-X data are discussed in [14, 15, 16]. *InSAR Time-Series Estimation of the Ionospheric Phase Delay: An Extension of the Split Range-Spectrum Technique* and InSAR data coherence estimation using 2D fast Fourier transform are performed in [17, 18].

In comparison with the results described in the aforementioned publications, the main goal of the present work is to suggest an analytical model of multi-pass InSAR geometry and derive analytical expressions of current distances between SAR’s positions and individual pixels on the surface and to describe principal InSAR parameters: topographic height and topographic displacement from the position of InSAR modelling. The focus is on the two modelling approaches: first, by the definition of real scenario, geometry, and kinematics and SAR signal models and corresponding complex image reconstruction and interferogram and differential interferogram generation and, second, the process of pseudo SAR measurements and interferogram generation that is analytically described. Results of numerical experiments with real data are provided.

The rest of the chapter is organized as follows. In Section 2, 3D InSAR geometry and kinematics are analytically described. In Section 3 and Section 4, analytical expressions of InSAR relief measurements and relief displacement measurements are presented. In Section 5 and Section 6, SAR waveform, deterministic signal model, and image reconstruction algorithm are described. In Section 7, numerical results of InSAR modelling based on the geometry, kinematics, and signal models are provided. In Section 8 and Section 9, a pseudo InSAR modelling of geophysical measurements and numerical results are presented, respectively. Conclusion remarks are made in Section 10.

## 2. InSAR geometry and kinematics

Assume a three-pass SAR system viewing three-dimensional (3-D) surface presented by discrete resolution elements, pixels. Each pixel is defined by the third coordinate

The basic geometric SAR characteristic is the time-dependent distance vector from SAR to each pixel on the surface in the *n*-th SAR pass at the *p*-th moment defined by

where *n* = 1–3 is the number of SAR passes and *n*-th SAR pass at the *p-*th moment, *n*-th SAR pass,

where *n*-th pass, defined by the following equation.

where *n*-th pass, measured at the initial moment;

Eq. (4) can be used to model a SAR signal from the *ij*-th pixel in the *n*-th SAR pass by calculation of the respective time delay and phase of the signal.

## 3. InSAR relief measurements

The distances to *m*-th and *n*-th pass (*m* ≠ *n*) at the moment of imaging can be defined by the cosine’s theorem, i.e.,

where *Oxy*. The look angle *ij*-th pixel on the surface with respect to *m*-th SAR position in the moment of imaging can be written as

The distance difference,

## 4. InSAR measurements of relief displacement

Consider a three-pass SAR interferometry (Figure 1). Let A and B be the two positions of imaging which can be defined by two passes of the same spaceborne SAR in different time (two pass interferometry). The third position C is defined by the third pass of the spaceborne SAR. The surface displacement,

where

Given the SAR wavelength λ, the phase differences proportional to range differences related to a particular pixel before and after displacement in the moment of imaging can be written as.

Neglecting the term

The displacement

For surface displacement

## 5. SAR waveform and deterministic signal model

The SAR transmits a series of electromagnetic waveforms to the surface, which are described analytically by the sequence of linear frequency modulation (chirp) pulses as follows

where *A* is the amplitude of the transmitted pulses, *M* is an emitted pulse number for synthesis of the aperture,

The SAR signal, reflected by *n*-th pass, can be expressed as

where

The parameter

SAR signal reflected from the entire illuminated surface is an interference of elementary signals of scattering pixels and can be written as

The time dwell *t* of the SAR signal return for each transmitted pulse *n*-th pass, *n*-th pass. Hence, in discrete form SAR signal can be rewritten as

The expressions derived in Section 2 and Section 5 can be used for modeling the SAR signal return in case the satellites are moving rectilinearly in 3-D coordinate system.

## 6. SAR image reconstruction

The complex image reconstruction includes the following operations: frequency demodulation, range compression, coarse range alignment, precise phase correction, and azimuth compression. The frequency demodulation is performed by multiplication of Eq. (20) with a complex conjugated function

Thus, the range distributed frequency demodulated SAR return in *n*-th pass for *p*-th pulse can be written as

The range compression of the LFM demodulated SAR signal is performed by cross correlation with a reference function,

where *K* is the full number of LFM samples, the range bins where SAR signal is registered, and by Fourier transform

for each

The range alignment and higher-order phase correction are beyond of the scope of the present work. The azimuth compression is accomplished by Fourier transform of the range compressed signal, *n*-th pass data can be expressed as

for each

The complex SAR image extracted from the *n*-th pass data preserves phases defined by distances from the satellite to each pixel at the moment of imaging. Based on pixel phases and image co-registration, a complex interferograms and differential interferograms can be created.

## 7. InSAR modeling: numerical results

The SAR signal model and imaging algorithm are illustrated by results of numerical experiments. Consider three pass satellite SAR system with position coordinates at the moment of imaging as follows.

Coordinates of vector-velocity of the satellite are

where *I* = 128 pixels; *J* = 128 pixels;

Normalized amplitude of reflected signals from every pixel ^{9} Hz. Frequency bandwidth is *K* = 512. Emitted pulse number is *M* = 512. Digital geometry description and SAR signal modeling are performed based on the theory in Sections 3 and 4. The complex images are extracted from the SAR signal by applying correlation range compression and FFT azimuth compression. Based on a priori-known kinematical parameters of satellites and coordinates of reference point from the surface autofocusing phase correction of the SAR signals registered in the both passes can be implemented.

The real and imaginary components of the SAR complex signal measured in the first SAR pass are depicted in Figure 2.

The complex SAR image’s amplitude and phase obtained in the first SAR pass are depicted in Figure 3. The orientation of the surface’s image (Figure 3a) in the frame is defined by the position of the SAR at the moment of imaging.

The real and imaginary components of the SAR complex signal measured in the second SAR pass are depicted in Figure 4.

The complex SAR image’s amplitude and phase obtained in the second SAR pass are depicted in Figure 5. It can be seen that the shape of the surface (the amplitude of the complex image) is similar to the shape of the surface obtained by the first SAR pass. In contrast, the phase structures of both complex images are different based on the different SAR positions in respect of the surface in the first and second pass at the moment of imaging.

By co-registration of the first and third SAR complex images, a complex SAR interferogram can be created with components in a coherent map and interferometric phase depicted in Figure 6.

The real and imaginary components of the SAR complex signal obtained in the third SAR pass is depicted in Figure 7.

The complex SAR image’s amplitude and phase obtained in the third SAR pass are depicted in Figure 8. The shape of the surface obtained in the third SAR pass is similar to the shape of the surface obtained by the first and second SAR passes. Comparing phase structures of the three complex SAR images, it can be noticed that they are different based on the different SAR’s positions in respect to the surface at the moment of imaging.

Under pixel co-registration of the first and third SAR complex images, a complex SAR interferogram can be created with components in a coherent map and interferometric phase depicted in Figure 9.

Due to precise under pixel co-registrations of the first and second and the first and third SAR complex images, the phase interferograms depicted in Figures 6b and 9b, respectively, are characterized with the similar structures.

## 8. Pseudo InSAR modeling of geophysical measurements

Consider three-pass InSAR geometry (Figure 1). The vector distances from the SAR positions to each

After distance measurements from the master SAR position A and slave SAR positions B and C, respectively, to each

without pixel displacement

with pixel displacement

In order to unwrap the interferometric phases, standard algorithms, MATLAB *unwrap* function, 2-D Costantini phase unwrapping based on network programming, and 2D Goldstein branch cut phase unwrapping, can be applied.

## 9. Pseudo InSAR geophysical measurements: numerical results

Consider a GeoTIFF file of Dilijan region in Caucasus, Armenia, located at the geographical coordinates 40° 44′ 27″ north and 44° 51′ 47″ east longitude. Consider 2-pass InSAR scenario. Coordinates of SAR positions in the moment of imaging are the following: master SAR position A, **,** b. Interferogram wrapped phases and unwrapped phases are presented in Figure 10c**,** d, respectively.

Consider a three-pass InSAR scenario and a surface before (Figure 11a) and after (Figure 11b) displacement described by MATLAB function *peaks*. Coordinates of SAR positions in the moment of imaging are the following: master SAR position A,

Distances to the surface at the moment of imaging as pseudo collar maps measured from SAR positions A, B, and C are presented in Figure 12a–c, respectively. AB interferogram without surface displacement and AC interferogram with surface displacement are presented in Figure 12d**,** e, respectively. Differential interferogram AB-AC is presented in Figure 12f.

The differential interferogram obtained by pixel subtraction of interferograms in Figure 12d**,** e is presented in Figure 12f. It illustrates the displacement of the surface. Only deformed part of the surface as differential fringes is depicted. The pseudo InSAR modeling can be applied to generate interferograms and differential interferograms based on real geophysical measurements and Geo TIFF maps of the observed surface.

## 10. Conclusions

A multi-pass InSAR system has been theoretically analyzed and numerically experimented. Geometry and kinematics of multi-pass InSAR scenario have been analytically described. Mathematical expressions for definition of current distance vectors between SAR system and surface’s pixels are derived. The basic InSAR parameters are defined. Analytical expressions to calculate pixel heights and pixel displacement have been derived. A model of linear frequency modulated SAR signal, reflected from the topographic surface, has been developed. An image reconstruction algorithm has been described. Numerical results verifying InSAR geometry, kinematics, and signal models are provided. Based on geometrical, kinematical, and signal models, numerical interferograms of a topographic surface have been created.

A pseudo InSAR approach has been applied to model processes of interferograms and differential interferogram generation using GeoTIFF files and measurements of distances from SAR positions to each pixels of the observed surface at the moment of imaging. Based on distance vector description of the InSAR scenario, the interferometric phase and interferometric differential phase have been analytically described. Pseudo InSAR geophysical measurements and interferograms and differential interferogram generation have been illustrated by results of numerical experiments.

In conclusion, the results in the present work can be applied for analysis and modeling of SAR interferometric processes in scenarios with different geometric, kinematics, and geological structures as well as for generating pseudo SAR interferograms based on the geophysical measurements and topographic maps.

## Conflict of interest

The author declares no conflict of interest.