Open access peer-reviewed chapter

InSAR Modeling of Geophysics Measurements

Written By

Andon Lazarov, Dimitar Minchev and Chavdar Minchev

Submitted: 13 June 2019 Reviewed: 22 August 2019 Published: 06 November 2019

DOI: 10.5772/intechopen.89293

From the Edited Volume

Geographic Information Systems in Geospatial Intelligence

Edited by Rustam B. Rustamov

Chapter metrics overview

1,099 Chapter Downloads

View Full Metrics

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.

Advertisement

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 zijxijyij in 3-D coordinate system Oxyz. Let A, B, and C, be the SAR positions of imaging. Between every SAR position, C32=3 InSAR baselines can be drawn.

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

Rijnp=RnpRij=xijnpyijnpzijnpT,E1

where n = 1–3 is the number of SAR passes and Rnp=R0n+V.p.Tp is the distance vector in the n-th SAR pass at the p-th moment, R0n is the initial distance vector in the n-th SAR pass, Rij is the constant distance vector of the ij- th pixel on the surface, and xijnp, yijnp, and zijnp are the current coordinates of Rijnp written by the expression.

xijnp=xnpxij,yijnp=ynpyij,zijnp=znpzijE2

where xij=iΔX, yij=jΔY, and zij=zijxijyij is the pixel’s discrete coordinates and xnp, ynp, and znp are the SAR current coordinates in the n-th pass, defined by the following equation.

xnp=x0nVxpTp,ynp=y0nVypTp,znp=z0nVzpTp,E3

where x0n, y0n, and z0n are the SAR initial coordinates in the n-th pass, measured at the initial moment; Tp is the time repetition period; p is the number of the emitted pulse; V=VxVyVzT is the SAR vector velocity; Vx=Vcosα, Vy=Vcosβ, and Vz=Vcosδ are the components of vector velocity; cosα, cosβ, and cosδ=1cos2αcos2β are the guiding cosines; and V is the module of the vector velocity V. Modulus of the current distance vector Rijnp is defined by

Rijnp=xijnp2+yijnp2+zijnp212.E4

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.

Advertisement

3. InSAR relief measurements

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

Rijn=Rijm2+Bmn22BmnRijmcosπ2θmijαmn]12,E5

where Bmn is the modulus of the baseline vector, θmij is the look angle, and αmn is a priory known tilt angle, the angle between the baseline vector and plane Oxy. The look angle θmij and height hm of an ij-th pixel on the surface with respect to m-th SAR position in the moment of imaging can be written as

θmij=αmn+arcsinRijm2+Bmn2Rijn22BmnRijm,E6
zij=hmRijmcosθmij.E7

The distance difference, ΔRijmn=RijnRijm, can be expressed by the interferometric phase difference ΔRijmn=λ2πΔϕijmn. In case Rijm can be measured, i.e., Rijn=Rijm+ΔRijmn, then

θmij=αmn+arcsinBmn2Rijmλ2πBmnΔϕijmn1+λ4πRijmΔϕijmn,E8
zij=hmRijm.cosαmn+arcsinBmn2Rijmλ2πBmnΔϕijmn.1+λ4πRijmΔϕijmn.E9
Advertisement

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, Δzij, due, for instance, to an earthquake could derive from two SAR interferograms built before and after the seismic impact. The temporal baseline, the time scale over which the displacement is measured, must follow the dynamics of the geophysical phenomenon. Short-time baseline is applied for monitoring fast surface changes. Long temporal baseline is used for monitoring slow geophysics phenomena (subsidence). The interferometry phase before event is derived from complex images acquired by A and B SAR positions in the moment of imaging, while the interferometry phase after event is derived from complex images acquired by A and C SAR positions in the moment of imaging. The distances Rij1,Rij2, Rij3, and Rijd3 after standard manipulations are written as.

Figure 1.

InSAR geometry and kinematics.

Rij2Rij1B1sinθijα1+B122Rij1,Rij3Rij1B2sinθijα2+B222Rij1,E10
Rijd3Rij3Δzcosθij+B2Rij1sinα2+Δzij22Rij1,

where Rij1,Rij2, and Rij3 are the slant ranges from A, B, and C positions of SAR system to the observed pixel in the moment of imaging before the surface displacement and Rijd3 is the slant range to from C SAR position to the observed pixel after Δzij surface displacement.

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.

ϕAB=4πλRij1Rij2,ϕAC=4πλRij1Rij3,ϕdAC=4πλRij1Rijd3E11

Neglecting the term Δz2/2Rij1 in Eq. (10) can be rewritten as.

ϕAB=4πλB1sinθijα1B122Rij1;ϕAC=4πλB2sinθijα2B222Rij1E12
ϕdAC=4πλB2sinθijα2B222Rij1+Δzcosθij+B2Rij1sinα2ϕAC+Δzcosθij+B2Rij1sinα2E13

The displacement Δzij is extracted from the differential interferometric phase difference ΔΦd=ϕdACϕAB. Considering B2/Rij1<<1, then ΔΦd=ΔΦ+4πλΔzijcosθij, where

ΔΦ=4πλB2sinθijα2B1sinθijα1B22B122Rij1.E14

For surface displacement zij can be written as

Δzij=λ4πΔΦdΔΦcosθij.E15
Advertisement

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

St=p=1MAexpjωtpTp+btpTp2,E16

where A is the amplitude of the transmitted pulses, Tp is the pulse repetition period, ω=2π.c/λ is the angular frequency, p=1,M¯ is the index of LFM emitted pulse, M is an emitted pulse number for synthesis of the aperture, c=3×108 m/s is the light speed in vacuum, ΔF is the LFM pulse bandwidth, b=π.ΔF/T is the chirp rate, and T is the time LFM pulse width.

The SAR signal, reflected by ij-th pixel and registered in the n-th pass, can be expressed as

Sijnt=aijzijrectttijnTexpjωttijn+bttijn2E17
rectttijnpT=1,0<ttijnpT1,E18

where aijzij is the reflection coefficient of the pixel from the surface.

The parameter aijzij is a function of surface geometry; tijnp=Rij1p+Rijnpc is the time propagation of the reflected signal from the ij-th scattering pixel registered in the n-th pass.

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

Snt=ijaijzijrectttijnTexpjωttijn+bttijn2.E19

The time dwell t of the SAR signal return for each transmitted pulse p can be expressed as t=tijminnp+kΔT, where k=kijminnp,kijmaxnp¯ is the sample number of the SAR return measured on range direction in n-th pass, kijminn=inttijminnp/ΔT, kijmaxn=inttijmaxnp/ΔT, ΔT=1/2ΔF is the sample time width, and kmaxnp is the number of the furthest range bin where SAR signal is registered in n-th pass. Hence, in discrete form SAR signal can be rewritten as

Ṡnkp=ijaijzijrectttijnTexpjωk1ΔTtijnp+bk1ΔTtijnp2.E20

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.

Advertisement

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 expjωk1ΔT+bk1ΔT2.

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

Ŝ̇nkp=ijaijzijrectk1ΔTtijnT.expjωtijnp+bk1ΔTtijnp2.E21

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

ṠRnk̂p=k=1KŜ̇nkpexpjbkk̂1ΔT2E22

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

ṠRnk̂p=k=1KŜ̇nkp.expj2πkk̂Kmaxn,E23

for each p=1,M¯ and k̂=1,K¯.

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, ṠRnk̂p. The complex image extracted from the n-th pass data can be expressed as

İnk̂p̂=p=1MṠRnpk̂expj2πpp̂M,E24

for each p̂=1,M¯, k̂=1,K¯.

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.

Advertisement

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.

x01=0m;y01=10.103m,z01=100.103m,x02=0m,y02=10,1.103m,z02=100.103m,
x03=0m,y03=10,2.103m,z03=100.103m.

Coordinates of vector-velocity of the satellite are vx=0 m/s, vy=600 m/s, and vz=0 m/s. The surface observed by the SAR system is modeled by the following equation

zij=31xij2expxij2yij+1210xij5xij3yij5expxij2yij213expxij+12yij2,E25

where xij=iΔX,yij=jΔY, i=1,I¯, j=1,J¯, I = 128 pixels; J = 128 pixels; ΔX; ΔY- the spatial resolution of the pixels.

Normalized amplitude of reflected signals from every pixel aij=0.001. The spatial resolution of the pixel are ΔX=ΔY=2 m. Wavelength is 0.03 m. Carrier frequency is 3.109 Hz. Frequency bandwidth is ΔF=250 MHz. Pulse repetition period is Tp=25.103 s. LFM pulse duration is T=5.106 s. Sample time duration is ΔT=1,95.108 s. LFM sample number 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.

Figure 2.

The real (a) and imaginary (b) component of the SAR complex signal measured in the first SAR pass.

Figure 3.

The amplitude (a) and phase (b) component of the SAR complex image obtained in the first pass.

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

Figure 4.

The real (a) and imaginary (b) component of the SAR complex signal measured in the second SAR pass.

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.

Figure 5.

The amplitude (a) and phase (b) components of the SAR complex image obtained in the second pass.

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.

Figure 6.

The coherent map (a) and interferometric phase (b) of the complex SAR interferogram created by the first and second SAR complex images.

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

Figure 7.

The real (a) and imaginary (b) component of the SAR complex signal measured in the third SAR pass.

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.

Figure 8.

The amplitude (a) and phase (b) component of the SAR complex image obtained in the third SAR pass.

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.

Figure 9.

The coherent map (a) and interferometric phase (b) of the complex SAR interferogram created by the first and third SAR complex images.

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.

Advertisement

8. Pseudo InSAR modeling of geophysical measurements

Consider three-pass InSAR geometry (Figure 1). The vector distances from the SAR positions to each ij-th pixel from the region of interest are RijS=RSRij, where S=A,B,C denotes the SAR position at the moment of imaging, RS=xSySzSTdenotes the SAR vector position, and Rij=xijyijzijTdenotes the ij-th pixel vector position. Coordinates of SAR positions in the moment of imaging are as follows: for a master SAR position A, xA, yA, zA; for a slave SAR position B, xB,yB, zB; and for a slave SAR position C, xC, yC, zC.

After distance measurements from the master SAR position A and slave SAR positions B and C, respectively, to each ij-th pixel on the surface and co-registration of so obtained master image and slave images, the instrumental interferometric phase differences are calculated as follows

  • without pixel displacement

ϕijAB=4πλRijARijB2π.max2λRijARij,ϕijAC=4πλRijARijC2π.max2λRijARijC

  • with pixel displacement

ϕij,dAC=4πλRijARij,dC2π.max2λRijARij,dC.

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.

Advertisement

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, xA=0 m, yA=300.3×103 m, zA=3×105 m and slave SAR position B xB=0 m, yB=300×103 m, zB=3×105 m. Wavelength is 0.05 m. Distances at the moment of imaging from the SAR position A and SAR position B to each pixel on the surface are illustrated in Figure 10a,b. Interferogram wrapped phases and unwrapped phases are presented in Figure 10c,d, respectively.

Figure 10.

Distances from SAR position A (a) and SAR position B (b) to each pixel on the surface in pseudo color map, interferogram wrapped phases (c) and unwrapped phases (d).

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, xA=350 km, yA=350 km, zA=800 km; slave SAR position B xB=351.5 km, yB=350 km, zB=800 km; and slave SAR position C xC=350 km, yC=351.2 km, zC=800 m. Wavelength is 0.03 m.

Figure 11.

Surface (peaks) before (a) and after (b) displacement.

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 12ac, 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.

Figure 12.

Distances to the surface measured from SAR positions A (a), B (b), and C (c). AB interferogram (d), AC interferogram (e) with surface displacement, and differential interferogram AB-AC (f).

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.

Advertisement

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.

References

  1. 1. Pepe A, Bonano M, Zhao Q, Yang T, Wang H. The use of C-/X-band time-gapped SAR data and geotechnical models for the study of Shanghai’s ocean-reclaimed lands through the SBAS-DInSAR technique. Remote Sensing. 2016;8:911. DOI: 10.3390/rs8110911
  2. 2. Hu J, Ding X, Li Z, et al. Kalman-filter based approach for multi-sensor, multitrack and multitemporal InSAR. IEEE Transactions on Geoscience and Remote Sensing. 2013;51:4226-4239. DOI: 10.1109/TGRS.2012.2227759
  3. 3. Casu F, Manconi A, Pepe A, Lanari R. Deformation time-series generation in areas characterized by large displacement dynamics: The SAR amplitude pixel-offset SBAS technique. IEEE Transactions on Geoscience and Remote Sensing. 2011;49:2752-2763. DOI: 10.1109/TGRS.2010.2104325
  4. 4. Samsonov S, d’Oreye N. Multidimensional time-series analysis of ground deformation from multiple InSAR data sets applied to Virunga Volcanic Province. Geophysical Journal International. 2012;191:1095-1108. DOI: 10.1111/j.1365-246X.2012.05669.x
  5. 5. Sansosti E, Berardino P, et al. How second-generation SAR systems are impacting the analysis of ground deformation. International Journal of Applied Earth Observation and Geoinformation. 2014;28:1-11. DOI: 10.1016/j.jag.2013.10.007. Available from: https://www.sciencedirect.com/science/article/pii/S0303243413001323 [Accessed: 11 August 2019]
  6. 6. Vettore A, Ponte S, et al. Space-based surface change detection with differential Synthetic Aperture Radar (SAR) Interferometry: Potentialities and preliminary investigations. In: Symposium on Geospatial Theory, Processing and Applications; Ottawa; 2002. Available from: https://www.semanticscholar.org/paper/Space-based-surface-change-detection-with-Synthetic-Ponte-Crocetto/10591de3ab1c236a30f2eeaff3ab1f3e8da1af17 [Accessed: 11 August 2019]
  7. 7. Wang J, Yu W, Deng Y, Wang R, et al. Demonstration of time-series InSAR processing in Beijing using a small stack of Gaofen-3 differential interferograms. Journal of Sensors. 2019;2019:4204580. DOI: 10.1155/2019/4204580
  8. 8. Pepe A, Calò F. Review of interferometric synthetic aperture radar (InSAR) multi-track approaches for the retrieval of Earth’s surface displacements. Applied Sciences. 2017;7:1264. DOI: 10.3390/app7121264
  9. 9. Fuhrmann T, Garthwaite M, Lawrie S, Brown N. Combination of GNSS and InSAR for future Australian datums. In: International Global Navigation Satellite Systems Association IGNSS Symposium 2018; Colombo Theatres, Kensington Campus, UNSW Australia; 7–9 February 2018. Available from: http://www.ignss2018.unsw.edu.au/sites/ignss2018/files/u80/Papers/IGNSS2018_paper_17.pdf [Accessed: 11 August 2019]
  10. 10. Xinshuang W, Lingling L, Xiaoliang Sh, Xitao H, Weil G. A high precision DEM extraction method based on InSAR data. In: ISPRS Annals of the Photogrammetry, Remote Sensing and Spatial Information Sciences, Volume IV-3, 2018 ISPRS TC III Mid-Term Symposium “Developments, Technologies and Applications in Remote Sensing”, 7–10 May, Beijing, China; 2018. Available from: https://pdfs.semanticscholar.org/7f92/aca4273b2afb1b9347724d72eb4c13bcd1cc.pdf [Accessed: 10 August 2019]
  11. 11. Crosetto M, Crippa B. Quality assessment of interferometric SAR DEMs. In: International Archives of Photogrammetry and Remote Sensing. Vol. XXXIII, Part B1; Amsterdam; 2000. Available from: https://pdfs.semanticscholar.org/376f/4f7f12602827801b954c7fb79aae2585b9c3.pdf [Accessed: 28 July 2019]
  12. 12. Letsios V, Faraslis I, Stathakis D. InSAR digital surface model (DSM) using Sentinel 1 and spatial data creation. In: AGILE 2019—Limassol; June 17–20, 2019. Available from: https://agile-online.org/images/conference_2019/documents/short_papers/50_Upload_your_PDF_file.pdf [Accessed: 28 July 2019]
  13. 13. Liang C, Agram P, Simons M, Fielding EJ. Ionospheric correction of InSAR time series analysis of C-band Sentinel-1 TOPS data. This paper is a non-peer reviewed preprint published at EarthArXiv and submitted for consideration by IEEE TGRS. DOI: 10.1109/TGRS.2019.2908494. Available from: https://eartharxiv.org/atxr7/download
  14. 14. Li Z, Bethel J. DEM registration, alignment and evaluation for SAR interferometry. In: The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences. Vol. XXXVII. Part B1. Beijing; 2008. pp. 11-116. Available from: https://www.isprs.org/proceedings/XXXVII/congress/1_pdf/19.pdf [Accessed: 07 August 2019]
  15. 15. Mao W-J, Chang W-L. Deformation monitoring by ground-based SAR interferometry (GB-InSAR): A field test in dam. Advances in Information Sciences and Service Sciences (AISS). 2015;7(2):133-140. Available from: http://www.globalcis.org/aiss/ppl/AISS3808PPL.pdf [Accessed: 05 August 2019]
  16. 16. Liu D, Chen R, Riedel B, Niemeier W. An improved approach to estimate large-gradient deformation using high resolution TerraSAR-X data. Solid Earth Discussions. 2014;6:2759-2778. DOI: 10.5194/sed-6-2759-2014. Available from: https://www.solid-earth-discuss.net/6/2759/2014/sed-6-2759-2014.pdf [Accessed: 01 August 2019]
  17. 17. Fattahi H, Simons M, Agram P. InSAR time-series estimation of the ionospheric phase delay: An Extension of the split Range-Spectrum technique. IEEE Transactions on Geoscience and Remote Sensing. 2017;55(2):777-786. Available from: http://web.gps.caltech.edu/∼simons/publications/pdfs/Fattahi_2017.pdf [Accessed: 01 August 2019]
  18. 18. Sosnovsky A, Kobernichenko V, Vinogradova N, Tsogtbaatar Ol. InSAR data coherence estimation using 2D fast Fourier transform. CEUR Workshop Proceedings (CEUR-WS.org). 2017;1814:98-105. Available from: http://ceur-ws.org/Vol-1814/paper-12.pdf [Accessed: 10 August 2019]

Written By

Andon Lazarov, Dimitar Minchev and Chavdar Minchev

Submitted: 13 June 2019 Reviewed: 22 August 2019 Published: 06 November 2019