CN106526590A - Method for monitoring and resolving three-dimensional ground surface deformation of industrial and mining area by means of multi-source SAR image - Google Patents

Method for monitoring and resolving three-dimensional ground surface deformation of industrial and mining area by means of multi-source SAR image Download PDF

Info

Publication number
CN106526590A
CN106526590A CN201610961794.7A CN201610961794A CN106526590A CN 106526590 A CN106526590 A CN 106526590A CN 201610961794 A CN201610961794 A CN 201610961794A CN 106526590 A CN106526590 A CN 106526590A
Authority
CN
China
Prior art keywords
deformation
radar
monitoring
dimensional
coherence
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201610961794.7A
Other languages
Chinese (zh)
Other versions
CN106526590B (en
Inventor
王志伟
于胜文
陶秋香
刘国林
高腾飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University of Science and Technology
Original Assignee
Shandong University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shandong University of Science and Technology filed Critical Shandong University of Science and Technology
Priority to CN201610961794.7A priority Critical patent/CN106526590B/en
Publication of CN106526590A publication Critical patent/CN106526590A/en
Application granted granted Critical
Publication of CN106526590B publication Critical patent/CN106526590B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

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

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a method for monitoring and resolving three-dimensional ground surface deformation of an industrial and mining area by means of a multi-source SAR image. The method is based on SAR image data of different rails before and after deformation. Deformation information in different radar satellite visual line directions are obtained through combining D-InSAR technology and Offset-tracking technology. Deformation information in different radar satellite azimuth directions is acquired through combining MAI technology and Offset-tracking technology. Based on this fact, different deformation information monitoring precisions in the viewing line direction and the azimuth direction are considered. Weights are determined in a Helmert variance component estimation manner. A three-dimensional land surface deformation conversion model is established and resolved, thereby obtaining a three-dimensional deformation field of the industrial and mining area. The method can settle a problem of easy restriction by factors such as decorrelation and small monitorable deformation gradient. Furthermore, the true ground surface deformation condition is more comprehensively monitored through the multi-source SAR image.

Description

A kind of fusion multi-source SAR images industrial and mining area three-dimensional earth's surface deformation monitoring and calculation method
Technical field
The present invention relates to a kind of fusion multi-source SAR images industrial and mining area three-dimensional earth's surface deformation monitoring and calculation method.
Background technology
Industrial and mining area Ground Deformation serious harm landing ground Architectural Equipment and natural environment, affect the mankind living environment, The security of the lives and property and local economic development.Carry out comprehensive monitoring to mining area surface deformation, further investigation deformation formation mechenism with Changing Pattern, to reasonable recovery of subterranean mineral resources, control industrial and mining area sustainable development is significant.
What traditional geodetic level(l)ing, GPS monitoring mining area surface deformation were obtained is all the deformation data of a point, with survey The development of means and further expanding for Ground Subsidence Monitoring range of application are painted, its shortcoming is also increasingly projected.
Spaceborne InSAR (Interferometric Synthetic Aperture Radar, InSAR) technology is earth observation from space new technique developed in recent years, its principle be based on deformation before and after two scape SAR The differential interferometry phase place of image obtains Ground Deformation information, such as D-InSAR, MAI technology.In the application of mining area surface deformation monitoring In obtain some successful cases both at home and abroad, monitoring accuracy has reached grade.But receive SAR image wavelength, angle of incidence, ground distributor The isoparametric impact of resolution, InSAR technologies only have monitoring capability to the deformation in the range of certain deformation gradient.Due to industrial and mining area There is different degrees of Ground Deformation with the carrying out of exploitation in earth's surface, is likely to be breached meter level, seriously beyond with difference in the short time The monitoring capability of the conventional InSAR technologies based on point interferometric phase information, causes SAR to affect dephasing to do.Secondly, InSAR skills Art can only obtain one-dimensional or two-dimentional earth's surface deformation data, rather than reflect the three-dimensional deformation of Ground Deformation feature directly, comprehensively .
Side-play amount tracking technique (Offset-Tracking) is by two width SAR images based on amplitude information Carry out precision registration, obtain between registration point side-play amount to calculate deformation, it is not necessary to carry out phase unwrapping, the phase to SAR image Dryness is insensitive, can make up that above-mentioned routine is dry by dephasing based on interferometric phase information technology and can to monitor deformation gradient little etc. The deficiency that factor is limited.In addition, the in-orbit SAR satellites of increasing different platform can provide multi-direction for areal Deformation monitoring result so that easily realized based on the monitoring of multi-source SAR image three-dimensionals Ground Deformation.At the same time different tracks SAR images are monitored to areal, can obtain more fully Ground Deformation information.
The content of the invention
It is an object of the invention to propose a kind of fusion multi-source SAR images industrial and mining area three-dimensional earth's surface deformation monitoring and resolving side Method, to overcome conventional InSAR technology to be done and single direction Ground Deformation information can only be provided and limited by big gradient deformation, dephasing System, beneficial to the deformation characteristics for really reflecting industrial and mining area Deformation Field comprehensively.
To achieve these goals, the present invention is adopted the following technical scheme that：
A kind of fusion multi-source SAR images industrial and mining area three-dimensional earth's surface deformation monitoring and calculation method, comprise the steps：
A obtain multi-source SAR images interfere to coherence map；
B is smoothed to coherence map；
C combines above-mentioned coherence map and merges D-InSAR technologies and Offset-tracking technical limit spacing radar line of sights to shape Change information；
D combines above-mentioned coherence map and merges MAI technologies and Offset-tracking technical limit spacings radar bearing and believes to deformation Breath；
The radar line of sight that e is obtained using step c and step d is adopted to deformation data and radar bearing to deformation data Helmert variance components estimates are set up three-dimensional earth's surface deformation model and are resolved.
Preferably, step a is specially：
The covering of collection research area's multi-source SAR images are carried out being imaged, after multiple look processing, is carried out based on Precise Orbit data Thick configuration is processed, and calculates initial offset；Again using the method for registering based on correlation coefficient, side-play amount multinomial is fitted, Processed by resampling after multinomial coefficient is calculated under criterion of least squares and complete precision registration；After each self registration SAR images carry out conjugate multiplication respectively, and the coherence map of SAR images is calculated while obtaining respective differential interferometry phase diagram.
Preferably, step c is specially：
C.1 according to InSAR technologies, maximum can monitor deformation gradient theory in actual applications, determine that D-InSAR can be monitored Coherence's threshold value；
C.2 the coherence map in the coherence's threshold value and step b that are c.1 obtained using step, is more than to studying coherent value in area Step c.1 middle threshold value region using D-InSAR technologies carry out Ground Deformation monitoring obtain radar line of sight to deformation data；
C.3 in research area coherent value less than c.1 middle threshold value region using Offset-tracking technologies at Reason, obtains distance and radar line of sight is further converted into after deformation data to deformation data；
C.4 by step c.2 with step c.3 in the radar line of sight that obtains respectively melted to deformation data according to pixel position Close, so as to obtain whole research area's radar line of sight to Ground Deformation information；
C.5 repeat the above steps c.1- c.4, so as to be respectively processed to different track SAR image datas, realize by step Monitoring of the different radar line of sights to deformation data.
Preferably, the step c.1 middle D-InSAR can monitor coherence's threshold value determination process it is as follows：
The computing formula that largest deformation gradient former can be monitored in theory is as follows：
Wherein, dmaxLargest deformation gradient, λ are wavelength, and u is Pixel size；
The relation of largest deformation gradient and coherence can actually be monitored：Dmax=dmax+0.002(γ-1)；
Wherein, DmaxCan actually to monitor largest deformation gradient, γ is coherent value；
Found out, with the reduction of coherent value, there is a less coherent value and cause D by formulamaxIt is changed into 0, i.e., prison is not measured Deformation data；Each sensor is calculated based on this can actually monitor the critical coherent value of largest deformation gradient as threshold value.
Preferably, step d is specially：
D.1 according to MAI technologies and Offset-tracking technical monitorings precision and the relation of coherence, two kinds of skills are chosen When art monitoring accuracy height changes, critical coherent value is used as coherence's threshold value；
D.2 the coherence map in the coherence's threshold value and step b d.1 determined according to step is big to studying coherent value in area Earth's surface orientation deformation monitoring is carried out using MAI technologies in the region of d.1 middle threshold value radar bearing is obtained to deformation data；
D.3 in research area coherent value less than d.1 middle threshold value region using Offset-tracking technologies at Reason, in acquisition research area, remaining area radar bearing is to deformation data；
D.4 the radar bearing that d.3 d.2 step obtained respectively with step is merged to deformation data, it is whole so as to obtain Individual research area radar bearing is to Ground Deformation information；
D.5 repeat the above steps d.1- d.4, so as to process to different track SAR image datas, realize different by step Monitoring of the radar bearing to deformation data.
Preferably, the step d.1 middle coherence's threshold value size be 0.8.
Preferably, step e is specially：
Multi-source SAR images radar line of sight is obtained to deformation data and radar bearing to after deformation data, by different resolution Monitoring result is resampled to equal resolution；Process is iterated by pixel using the method for Helmert variance components estimates, is counted The weight of each pixel is calculated, three-dimensional Ground Deformation is set up and is resolved model and resolved.
Preferably, the process three-dimensional Ground Deformation resolving model set up in step e and is resolved is as follows：
Three-dimensional Ground Deformation resolves model and is expressed as：R=Bd；
Wherein, d=(du,de,dn)TRepresent three-dimensional earth's surface deformation data, duRepresent vertical direction deformation data, deRepresent east West is to deformation data, dnRepresent North and South direction deformation data；Represent difference in radar fix system Direction deformation data,1 sight line of track is represented to deformation data,2 sight line of track is represented to deformation data,Represent track 1 orientation deformation data,Represent 2 orientation deformation data of track；B is three-dimensional deformation model coefficient matrix；
It is expressed as according to SAR satellite earth observation principle schematics：
Wherein, α1Represent 1 radar bearing angle of track, α2Represent 2 radar bearing angle of track, θ1Represent that 1 radar of track is incident Angle, θ2Represent 2 radar angle of incidence of track；
According to the principle of least square, earth's surface three-dimensional deformation information is calculated：D=(BTPB)-1 BTPR；
Wherein, P be two different track sight lines to the power battle array corresponding to orientation deformation quantity；According to Helmert variances Component estimation formulas, the relational expression between component of variance and observation residual error is expressed as：S θ=Wθ；
In formula,Represent error in the weight unit of two class observed quantities；And WθIt is respectively with S：
In formula,Ni=Bi TPiBi, i=1,2；P1And P2The power battle array of two class observed quantities is represented respectively；tr Representing matrix asks mark computing；V1And V2The correction of two class observed quantities is represented respectively；
Given initial power battle array P, estimates to θ according to Helmert Variance Components Estimations after first time adjustment, then Weight is recalculated according to the following formula：
Wherein,WithRepresent that the valuation of battle array is weighed in two class observed quantities respectively；OrderMeter is iterated again Calculate, untilAgain weigh surely, and then calculate three-dimensional deformation information.
The invention has the advantages that：
The inventive method is monitored to industrial and mining area Ground Deformation based on multi-source SAR image data, can be more comprehensively true Ground obtains the space-time characterisation of industrial and mining area Ground Deformation.For industrial and mining area Ground Deformation feature, can by merging various InSAR technologies So that effectively solving dephasing is dry, the factor such as deformation gradient is big restriction, InSAR technologies have been widened in industrial and mining area large scale Ground Deformation The application prospect of monitoring, the cost and technical limitations that reduce change monitoring.On this basis, mould is resolved by setting up three-dimensional deformation Type, more comprehensively can effectively obtain industrial and mining area three-dimensional earth's surface deformation data, more directly understand and grasp research area's earth's surface shape Change situation, effectively pre- Geological disaster prevention measure.Also the research to carry out Ground Deformation mechanism and Deformation Law simultaneously is provided more Direct basic data, this is significant to mining area suitable development.
Description of the drawings
Fig. 1 is the stream of a kind of fusion multi-source SAR images industrial and mining area three-dimensional earth's surface deformation monitoring and calculation method in the present invention Cheng Tu；
Fig. 2 is to deformation monitoring in the present invention using the rail lift ALOS PALSAR data sight lines after D-InSAR technical finesses Result schematic diagram；
Fig. 3 is to merge D-InSAR technologies and the rail lift ALOS after Offseet-tracking technical finesses in the present invention PALSAR data sight line is to deformation monitoring result schematic diagram；
Fig. 4 is to be supervised to deformation using the drop rail ENVISAT ASAR data sight lines after D-InSAR technical finesses in the present invention Survey result schematic diagram；
Fig. 5 is to merge D-InSAR technologies and the drop rail ENVISAT after Offseet-tracking technical finesses in the present invention ASAR data sight line is to deformation monitoring result schematic diagram；
Fig. 6 is using the rail lift ALOS PALSAR data orientation deformation monitoring results after MAI technical finesses in the present invention Schematic diagram；
Fig. 7 is to merge MAI technologies and the rail lift ALOS PALSAR after Offseet-tracking technical finesses in the present invention Data orientation deformation monitoring result schematic diagram；
Fig. 8 is using the drop rail ENVISAT ASAR data orientations deformation monitoring knot after MAI technical finesses in the present invention Fruit schematic diagram；
Fig. 9 is to merge MAI technologies and the drop rail ENVISAT of the result after Offseet-tracking technical finesses in the present invention ASAR data orientation deformation monitoring result schematic diagrams；
Figure 10 is the vertical direction result schematic diagram of three-dimensional earth's surface deformation monitoring in the present invention；
Figure 11 is the east-west direction result schematic diagram of three-dimensional earth's surface deformation monitoring in the present invention；
Figure 12 is the North and South direction result schematic diagram of three-dimensional earth's surface deformation monitoring in the present invention；
Figure 13 is SAR satellite earth observation principle schematics in the present invention.
Specific embodiment
The present invention ultimate principle be：
By being processed to covering multi-source SAR images in same research area using various InSAR technologies, will be various InSAR technologies carry out fusion and obtain the multi-direction Ground Deformation information in the region, further merge multi-source monitoring result and resolve three-dimensional Ground Deformation.
Below in conjunction with the accompanying drawings and specific embodiment is described in further detail to the present invention：
It is with reference to shown in Fig. 1, a kind of to merge multi-source SAR images industrial and mining area three-dimensional earth's surface deformation monitoring and calculation method, including Following several steps：
A obtain multi-source SAR images interfere to coherence map.Specifically acquisition methods are：
The covering of collection research area's multi-source SAR images are carried out being imaged, after multiple look processing, is carried out based on Precise Orbit data Thick configuration is processed, and calculates initial offset；Again using the method for registering based on correlation coefficient, side-play amount multinomial is fitted, Processed by resampling after multinomial coefficient is calculated under criterion of least squares and complete precision registration；After each self registration SAR images carry out conjugate multiplication respectively, and the coherence map of SAR images is calculated while obtaining respective differential interferometry phase diagram.
Wherein, coherence map refers to the foundation for evaluating two width SAR image similarity degrees, generates in differential interferometry process.
B is smoothed to coherence map.
Coherence map is the basis that follow-up data is processed, but as the impact of the factors such as various noises is led in actual process In causing the coherence map for obtaining, each pixel coherent value spatial continuity is poor, i.e., between adjacent picture elements, coherent value is changed greatly, it is impossible to The difference accessible successive range of InSAR technologies in subsequent processes is determined based on this.Therefore, subsequently located First the coherence map of coherent value spatial spreading is smoothed before reason.
In order to the verity for ensureing coherence map as far as possible can be using the method for Kriging regression in processing procedure.
C merges D-InSAR technologies and Offset-tracking technical limit spacing radar line of sights to deformation data.
To obtain industrial and mining area high accuracy earth's surface deformation monitoring result, easily occur what dephasing was done for industrial and mining area Ground Deformation Phenomenon, realizes full basin Ground Deformation monitoring by merging D-InSAR technologies and Offset-tracking technical approach.
Step c is specially：
C.1 largest deformation gradient theory can be monitored in actual applications according to InSAR technologies, determines that D-InSAR can be monitored Coherence's threshold value.
Industrial and mining area Ground Deformation has the features such as coverage is little, deformation magnitude is big, big gradient deformation easily occurs, poor Divide in interferogram and show as contrast of fringes comparatively dense, cause solution to twine mistake and then reliable deformation data cannot be obtained；And deformation Gradient is excessive to be easily caused deformation area and dephasing occurs and do, and reduces the precision of monitoring result.
InSAR technologies maximum can monitor deformation gradient theory to be proposed by Massonnett and Feigl earliest, and gives Theoretical model is gone out.On this basis, Baran et al. considers that the factor such as coherence affects, and in finding practical application, InSAR is most Deformation gradient value can be monitored greatly less than theoretical value, maximum can monitor deformation gradient in actual applications to give InSAR technologies Model.
Step of the present invention c.1 middle D-InSAR can monitor coherence's threshold value determination process it is as follows：
Massonnett and Feigl advanced a theory in 1998 on can monitor largest deformation gradient former computing formula such as Under：
Wherein, dmaxLargest deformation gradient, λ are wavelength, and u is Pixel size.
But by as factor is affected by dephasing is dry etc., largest deformation gradient theory can be monitored in actual application Model is less than theoretic.Affect in view of factors such as coherences, maximum can be monitored deformation gradient and be wanted than theoretical value in practice Little, Baran et al. gave the relation that can actually monitor largest deformation gradient and coherence in 2005：
Dmax=dmax+0.002(γ-1)；Wherein, DmaxCan actually to monitor largest deformation gradient, γ is coherent value.
Found out, with the reduction of coherent value, there is a less coherent value and cause D by formulamaxIt is changed into 0, i.e., prison is not measured Deformation data；Each sensor is calculated based on this can actually monitor the critical coherent value of largest deformation gradient as threshold value.
Wherein, coherent value can carry out Ground Deformation monitoring using D-InSAR technologies more than the region of threshold value, and coherent value is little Ground Deformation monitoring can be carried out using Offset-tracking technologies in the region of threshold value.
C.2 the coherence map in the coherence's threshold value and step b that are c.1 obtained using step, is more than to studying coherent value in area Step c.1 middle threshold value region using D-InSAR technologies carry out Ground Deformation monitoring obtain radar line of sight to deformation data.
C.3 in research area coherent value less than c.1 middle threshold value region using Offset-tracking technologies at Reason, obtains distance and radar line of sight is further converted into after deformation data to deformation data.
Offset_tracking technologies are to calculate registering side-play amount based on two scape image amplitude informations, will further be matched somebody with somebody Quasi- side-play amount is converted to deformation data, calculates orientation deformation on this basis.
The technology maximum feature is insensitive to coherence, it is not necessary to which phase unwrapping just can obtain deformation data, Ke Yiyou Effect overcomes the application of D-InSAR technologies to limit to, and the area that dephasing closes weight tight can be caused to large deformation gradient to provide deformation details.
Fig. 2 and Fig. 4 represent respectively two groups of different SAR images of ALOS and ASAR using D-InSAR technical limit spacings sight line to Deformation monitoring result, in figure, white space as loses coherent area movement.For the region is not using requiring to coherence Offset-tracking technologies are processed, and Fig. 3 and Fig. 5 is to merge final monitoring result after two kinds of technologies, by comparison diagram 2 and Fig. 3 and Fig. 4 and Fig. 5 can be found that effect is obvious, realize the monitoring of Low coherence area Ground Deformation.In addition, It can be found that different track SAR images are observed to earth's surface with different view by Fig. 2 to Fig. 5, obtains and be included in deformation The different monitoring results of aspect such as scope, deformation distribution, can obtain more fully Ground Deformation information.
C.4 by step c.2 with step c.3 in the radar line of sight that obtains respectively melted to deformation data according to pixel position Close, so as to obtain whole research area's radar line of sight to Ground Deformation information；
C.5 in order to obtain research area three-dimensional earth's surface deformation data, at least three different directions deformation datas are needed.In repetition C.1- c.4 step, so as to be respectively processed to different track SAR image datas, realizes different radar line of sights to shape to state step The monitoring of change information, can obtain more comprehensively real surface deformation data.
D merges MAI technologies and Offset-tracking technical limit spacing radar bearings to deformation data.
MAI technologies are that, based on interferometric phase information earth observation technology, different from InSAR technologies, the technology is used divides The InSAR Processing Algorithm for splitting wave beam obtains the orientation deformation quantity of degree of precision：
Two pairs of principal and subordinate's images are regarded first with two width haplopias plural number image (SLC) division generation is forward and backward, by forward and backward regarding master Interference treatment is carried out respectively from image obtains forward sight and backsight interferogram, it is then forward and backward regarding interferogram conjugate multiplication by this two Multiple aperture differential interferometry figure of the width comprising orientation deformation data is produced, after filtering, after the process of step such as phase unwrapping Obtain the deformation quantity of orientation.Similar with the mechanism that conventional InSAR monitors deformation, the coherence of SAR images is to govern the skill The principal element of art application, the quality of coherence directly affect the precision that MAI measures orientation deformation.Relative to based on amplitude The Offset-tracking technologies of information, its monitoring accuracy of high coherent area movement are high, and its monitoring accuracy of Low coherence area compared with It is low.
Image coherence affects larger to the precision of MAI technologies, and the general mulch farming thing of industrial and mining area earth's surface, coherence compared with Difference, in addition, in the short time, to have deepened image dephasing dry for big magnitude Ground Deformation.In order to obtain the employing of orientation deformation data Fusion MAI technologies and the Offset-tracking technology insensitive to coherence carry out Ground Deformation monitoring.
It is 1/50 pixel that Offset-tracking technologies calculate the precision of side-play amount, for SAR data, equivalent to Precision in orientation can reach 7.5cm and sight line to for 14cm.
MAI phase standards difference can be expressed as
With the raising of coherence, the phase standard difference for calculating is less and less.Show in pertinent literature when coherent value is big When 0.8, the standard deviation of MAI is better than Offset-tracking technologies.Therefore, the inventive method is big to studying coherent value in area Processed using MAI technologies in 0.8 region, Offset- is adopted to studying region of the coherent value less than 0.8 in area Tracking technologies are processed, and finally, two monitoring results are carried out fusion and realizes full basin Ground Deformation monitoring.
Step d is specially：
D.1 according to MAI technologies and Offset-tracking technical monitorings precision and the relation of coherence, two kinds of skills are chosen When art monitoring accuracy height changes, used as coherence's threshold value, its size is 0.8 to critical coherent value.
D.2 the coherence map in the coherence's threshold value and step b d.1 determined according to step is big to studying coherent value in area Earth's surface orientation deformation monitoring is carried out using MAI technologies in the region of d.1 middle threshold value radar bearing is obtained to deformation data.
D.3 in research area coherent value less than d.1 middle threshold value region using Offset-tracking technologies at Reason, in acquisition research area, remaining area radar bearing is to deformation data.
D.4 the radar bearing that d.3 d.2 step obtained respectively with step is merged to deformation data, it is whole so as to obtain Individual research area radar bearing is to Ground Deformation information.
Fig. 6 represents orientation deformation of the two groups of different SAR images of ALOS and ASAR using MAI technical limit spacings respectively with Fig. 8 Monitoring result, in figure, white space as loses coherent area movement.Lose coherent area movement and adopt Offset-tracking technical limit spacing deformation Information, Fig. 7 and Fig. 9 are to merge final monitoring result after two kinds of technologies.Can be with by contrasting Fig. 6 and Fig. 7 and Fig. 8 and Fig. 9 It was found that effect is obvious, by merge two kinds of technical limit spacings including setting basin center Low coherence zone aspect to deformation Information, realizes the monitoring of Low coherence area Ground Deformation.
D.5 repeat the above steps d.1- d.4, so as to process to different track SAR image datas, realize different by step Monitoring of the radar bearing to deformation data, can obtain more comprehensively real surface deformation data.
So far, multi-source SAR images are processed by merging various InSAR technologies, obtaining industrial and mining area includes two Different sight to two different azimuths in four interior different directions deformation datas.
E obtain multi-source SAR images (at least two different tracks) sight line to orientation deformation data after, by different resolutions Rate monitoring result is resampled to equal resolution.
Monitoring result is divided into by the difference for obtaining monitoring result precision in view of different technologies according to the difference of acquisition modes Two classes, are iterated process by pixel using the method for Helmert variance components estimates, calculate the weight of each pixel, build Vertical three-dimensional Ground Deformation resolves model and is resolved.
Figure 10, Figure 11 and Figure 12 respectively study three-dimensional the vertical of earth's surface deformation monitoring, thing and North and South direction result figure.
Wherein, Helmert variance components estimates are estimating using the weighted sum of squares of all kinds of corrections after each adjustment The variance of unit weight of all kinds of observations, and then realize the weight of all kinds of observations is determined according to precision.
Step e is specially：Three-dimensional Ground Deformation resolves model and is expressed as：R=Bd；
Wherein, d=(du,de,dn)TRepresent three-dimensional earth's surface deformation data, duRepresent vertical direction deformation data, deRepresent east West is to deformation data, dnRepresent North and South direction deformation data；Represent difference in radar fix system Direction deformation data,1 sight line of track is represented to deformation data,2 sight line of track is represented to deformation data,Represent track 1 orientation deformation data,Represent 2 orientation deformation data of track；B is three-dimensional deformation model coefficient matrix；
As shown in figure 13, SAR satellite earth observations principle schematic is represented by：
Wherein, α1Represent 1 radar bearing angle of track, α2Represent 2 radar bearing angle of track, θ1Represent that 1 radar of track is incident Angle, θ2Represent 2 radar angle of incidence of track；
According to the principle of least square, earth's surface three-dimensional deformation information is calculated：D=(BTPB)-1 BTPR；
Wherein, P be two different track sight lines to the power battle array corresponding to orientation deformation quantity；According to Helmert variances Component estimation formulas, the relational expression between component of variance and observation residual error is expressed as：S θ=Wθ；
In formula,Represent error in the weight unit of two class observed quantities；And WθIt is respectively with S：
In formula,Ni=Bi TPiBi, i=1,2；P1And P2The power battle array of two class observed quantities is represented respectively；tr Representing matrix asks mark computing；V1And V2The correction of two class observed quantities is represented respectively；
Given initial power battle array P, estimates to θ according to Helmert Variance Components Estimations after first time adjustment, then Weight is recalculated according to the following formula：
Wherein,WithRepresent that the valuation of battle array is weighed in two class observed quantities respectively；OrderMeter is iterated again Calculate, untilAgain weigh surely, and then calculate three-dimensional deformation information.
Certainly, only presently preferred embodiments of the present invention described above, the present invention are not limited to enumerate above-described embodiment, should When explanation, any those of ordinary skill in the art are under the teaching of this specification, all equivalent substitutes for being made, bright Aobvious variant, all falls within the essential scope of this specification, ought to be protected by the present invention.

Claims (8)

1. it is a kind of to merge multi-source SAR images industrial and mining area three-dimensional earth's surface deformation monitoring and calculation method, it is characterised in that including as follows Step：
A obtain multi-source SAR images interfere to coherence map；
B is smoothed to coherence map；
C combines above-mentioned coherence map and merges D-InSAR technologies and Offset-tracking technical limit spacings radar line of sight and believes to deformation Breath；
D combines above-mentioned coherence map and merges MAI technologies and Offset-tracking technical limit spacing radar bearings to deformation data；
The radar line of sight that e is obtained using step c and step d to deformation data and radar bearing to deformation data, using Helmert Variance components estimate is set up three-dimensional earth's surface deformation model and is resolved.
2. one kind according to claim 1 merges multi-source SAR images industrial and mining area three-dimensional earth's surface deformation monitoring and calculation method, Characterized in that, step a is specially：
The covering of collection research area's multi-source SAR images are carried out being imaged, after multiple look processing, is slightly matched somebody with somebody based on Precise Orbit data Process is put, initial offset is calculated；Again using the method for registering based on correlation coefficient, side-play amount multinomial is fitted, most A young waiter in a wineshop or an inn is taken advantage of to be processed by resampling after calculating multinomial coefficient under criterion and completes precision registration；By the SAR shadows after each self registration As carrying out conjugate multiplication respectively, the coherence map of SAR images while obtaining respective differential interferometry phase diagram, is calculated.
3. one kind according to claim 1 merges multi-source SAR images industrial and mining area three-dimensional earth's surface deformation monitoring and calculation method, Characterized in that, step c is specially：
C.1 according to InSAR technologies, maximum can monitor deformation gradient theory in actual applications, determine that D-InSAR can be monitored relevant Property threshold value；
C.2 the coherence map in the coherence's threshold value and step b that are c.1 obtained using step, is more than step to studying coherent value in area C.1 the region of middle threshold value carries out Ground Deformation monitoring using D-InSAR technologies and obtains radar line of sight to deformation data；
C.3 in research area coherent value processed using Offset-tracking technologies less than the region of c.1 middle threshold value, Obtain distance and radar line of sight is further converted into after deformation data to deformation data；
C.4 by step c.2 with step c.3 in the radar line of sight that obtains respectively merged to deformation data according to pixel position, So as to obtain whole research area's radar line of sight to Ground Deformation information；
C.5 repeat the above steps c.1- c.4, so as to be respectively processed to different track SAR image datas, realize different by step Monitoring of the radar line of sight to deformation data.
4. one kind according to claim 3 merges multi-source SAR images industrial and mining area three-dimensional earth's surface deformation monitoring and calculation method, Characterized in that, the step c.1 middle D-InSAR can monitor coherence's threshold value determination process it is as follows：
The computing formula that largest deformation gradient former can be monitored in theory is as follows：
Wherein, dmaxLargest deformation gradient, λ are wavelength, and u is Pixel size；
The relation of largest deformation gradient and coherence can actually be monitored：Dmax=dmax+0.002(γ-1)；
Wherein, DmaxCan actually to monitor largest deformation gradient, γ is coherent value；
Found out, with the reduction of coherent value, there is a less coherent value and cause D by formulamaxIt is changed into 0, i.e., prison does not measure deformation Information；Each sensor is calculated based on this can actually monitor the critical coherent value of largest deformation gradient as threshold value.
5. one kind according to claim 1 merges multi-source SAR images industrial and mining area three-dimensional earth's surface deformation monitoring and calculation method, Characterized in that, step d is specially：
D.1 according to MAI technologies and Offset-tracking technical monitorings precision and the relation of coherence, two kinds of technology prisons are chosen When surveying precision height change, critical coherent value is used as coherence's threshold value；
D.2 the coherence map in the coherence's threshold value and step b d.1 determined according to step, to studying in area coherent value more than d.1 The region of middle threshold value carries out earth's surface orientation deformation monitoring using MAI technologies and obtains radar bearing to deformation data；
D.3 in research area coherent value processed using Offset-tracking technologies less than the region of d.1 middle threshold value, In acquisition research area, remaining area radar bearing is to deformation data；
D.4 the radar bearing that d.3 d.2 step obtained respectively with step is merged to deformation data, is entirely ground so as to obtain Area's radar bearing is studied carefully to Ground Deformation information；
D.5 d.1- d.4 step, so as to process to different track SAR image datas, realizes different radars to repeat the above steps The monitoring of orientation deformation data.
6. one kind according to claim 5 merges multi-source SAR images industrial and mining area three-dimensional earth's surface deformation monitoring and calculation method, Characterized in that, the step d.1 middle coherence's threshold value size be 0.8.
7. one kind according to claim 1 merges multi-source SAR images industrial and mining area three-dimensional earth's surface deformation monitoring and calculation method, Characterized in that, step e is specially：
Acquisition multi-source SAR images radar line of sight, to after deformation data different resolution is monitored to deformation data and radar bearing As a result it is resampled to equal resolution；Process is iterated by pixel using the method for Helmert variance components estimates, is calculated The weight of each pixel, sets up three-dimensional Ground Deformation and resolves model and resolved.
8. one kind according to claim 7 merges multi-source SAR images industrial and mining area three-dimensional earth's surface deformation monitoring and calculation method, Characterized in that, the process three-dimensional Ground Deformation resolving model set up in step e and is resolved is as follows：
Three-dimensional Ground Deformation resolves model and is expressed as：R=Bd；
Wherein, d=(du,de,dn)TRepresent three-dimensional earth's surface deformation data, duRepresent vertical direction deformation data, deRepresent between east and west To deformation data, dnRepresent North and South direction deformation data；Represent different directions in radar fix system Deformation data,1 sight line of track is represented to deformation data,2 sight line of track is represented to deformation data,Represent 1 side of track Position to deformation data,Represent 2 orientation deformation data of track；B is three-dimensional deformation model coefficient matrix；
It is expressed as according to SAR satellite earth observation principle schematics：
B = cosθ 1 - cos ( α 1 - 3 π 2 ) sinθ 1 - sin ( α 1 - 3 π 2 ) sinθ 1 cosθ 2 - cos ( α 2 - 3 π 2 ) sinθ 2 - sin ( α 2 - 3 π 2 ) sinθ 2 0 sin ( α 1 - 3 π 2 ) - cos ( α 1 - 3 π 2 ) 0 sin ( α 2 - 3 π 2 ) - cos ( α 2 - 3 π 2 ) ;
Wherein, α1Represent 1 radar bearing angle of track, α2Represent 2 radar bearing angle of track, θ1Represent 1 radar angle of incidence of track, θ2Table Show 2 radar angle of incidence of track；
According to the principle of least square, earth's surface three-dimensional deformation information is calculated：D=(BTPB)-1BTPR；
Wherein, P be two different track sight lines to the power battle array corresponding to orientation deformation quantity；According to Helmert components of variance Estimation formulas, the relational expression between component of variance and observation residual error is expressed as：S θ=Wθ；
In formula,Represent error in the weight unit of two class observed quantities；And WθIt is respectively with S：
W θ = V 1 T P 1 V 1 V 2 T P 2 V 2 ; S = 2 - 2 t r ( N 1 N - 1 ) + t r ( N 1 N - 1 ) 2 t r ( N 1 N - 1 N 2 N - 1 ) t r ( N 2 N - 1 N 1 N - 1 ) 2 - 2 t r ( N 2 N - 1 ) + t r ( N 2 N - 1 ) 2 ;
In formula,Ni=Bi TPiBi, i=1,2；P1And P2The power battle array of two class observed quantities is represented respectively；Tr is represented Matrix Calculating mark computing；V1And V2The correction of two class observed quantities is represented respectively；
Given initial power battle array P, estimates to θ according to Helmert Variance Components Estimations after first time adjustment, then according to Following formula recalculates weight：
Wherein,WithRepresent that the valuation of battle array is weighed in two class observed quantities respectively；OrderCalculating is iterated again, UntilAgain weigh surely, and then calculate three-dimensional deformation information.
CN201610961794.7A 2016-11-04 2016-11-04 A kind of fusion multi-source SAR image industrial and mining area three-dimensional earth's surface deformation monitorings and calculation method Active CN106526590B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610961794.7A CN106526590B (en) 2016-11-04 2016-11-04 A kind of fusion multi-source SAR image industrial and mining area three-dimensional earth's surface deformation monitorings and calculation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610961794.7A CN106526590B (en) 2016-11-04 2016-11-04 A kind of fusion multi-source SAR image industrial and mining area three-dimensional earth's surface deformation monitorings and calculation method

Publications (2)

Publication Number Publication Date
CN106526590A true CN106526590A (en) 2017-03-22
CN106526590B CN106526590B (en) 2018-09-25

Family

ID=58326298

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610961794.7A Active CN106526590B (en) 2016-11-04 2016-11-04 A kind of fusion multi-source SAR image industrial and mining area three-dimensional earth's surface deformation monitorings and calculation method

Country Status (1)

Country Link
CN (1) CN106526590B (en)

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106842199A (en) * 2017-01-11 2017-06-13 湖南科技大学 It is a kind of to merge the method that different resolution SAR data monitors Ground Deformation
CN107037428A (en) * 2017-03-27 2017-08-11 中国科学院遥感与数字地球研究所 It is a kind of to improve the method that spaceborne dual station difference InSAR extracts deformation precision
CN107102332A (en) * 2017-05-11 2017-08-29 中南大学 The three-dimensional earth's surface deformation monitoring methods of InSAR based on variance components estimate and strees strain model
CN108921174A (en) * 2018-06-04 2018-11-30 中国矿业大学 Based on the big gradient deformation extracting method in mining area under the constraint of multistage match window
CN110017762A (en) * 2019-03-01 2019-07-16 四川省地质工程勘察院 A method of utilizing Offset-Tracking technical monitoring Wriggle-changing Body deformation
CN110058234A (en) * 2019-05-20 2019-07-26 太原理工大学 A method of it resolving mining area surface and settles three-dimensional deformation
CN110702050A (en) * 2019-11-12 2020-01-17 山东科技大学 Surface three-dimensional deformation monitoring method based on InSAR
CN110765934A (en) * 2019-10-22 2020-02-07 成都理工大学 Geological disaster identification method based on multi-source data fusion
CN111650579A (en) * 2020-06-12 2020-09-11 中南大学 InSAR mining area three-dimensional deformation estimation method and device for rock migration parameter adaptive acquisition and medium
CN111780660A (en) * 2020-07-13 2020-10-16 内蒙古工业大学 Mining area three-dimensional multi-level deformation optimization method and optimization device
CN111965646A (en) * 2020-10-20 2020-11-20 航天宏图信息技术股份有限公司 Satellite-borne radar data processing method, device, equipment and storage medium
CN112147610A (en) * 2020-09-28 2020-12-29 武汉大学 InSAR interferogram down-sampling method fusing topology analysis technology
CN112184902A (en) * 2020-09-21 2021-01-05 东华理工大学 Underground mining face inversion method for boundary crossing mining identification
CN112284332A (en) * 2020-08-31 2021-01-29 北京四象爱数科技有限公司 High-rise building settlement monitoring result three-dimensional positioning method based on high-resolution INSAR
CN112540370A (en) * 2020-12-08 2021-03-23 鞍钢集团矿业有限公司 InSAR and GNSS fused strip mine slope deformation measurement method
CN112667970A (en) * 2021-02-04 2021-04-16 淮阴师范学院 Deformation monitoring method for super-huge structure combined by multi-source observation technology
CN112965065A (en) * 2021-02-03 2021-06-15 北京观微科技有限公司 Mining area detection method cooperating with SAR image phase, intensity and coherence information
CN113091598A (en) * 2021-04-06 2021-07-09 中国矿业大学 Method for defining stability grade range of goaf building site by InSAR
CN114236541A (en) * 2021-12-08 2022-03-25 电子科技大学 Large-area earth surface three-dimensional deformation calculation method based on Sentinel-1 satellite SAR image
CN115512239A (en) * 2022-10-18 2022-12-23 闽江学院 Satellite platform multi-star sensor attitude self-adaptive fusion method
CN116047519A (en) * 2023-03-30 2023-05-02 山东建筑大学 Point selection method based on synthetic aperture radar interferometry technology
CN117213443A (en) * 2023-11-07 2023-12-12 江苏省地质调查研究院 Construction and updating method of ground settlement monitoring network with integration of heaves, earth and depth
CN117970331A (en) * 2024-04-02 2024-05-03 中国电建集团华东勘测设计研究院有限公司 InSAR earth surface deformation monitoring method and system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6911931B2 (en) * 2002-10-24 2005-06-28 The Regents Of The University Of California Using dynamic interferometric synthetic aperature radar (InSAR) to image fast-moving surface waves
CN102645650A (en) * 2012-03-06 2012-08-22 北京北科安地科技发展有限公司 Landslide dynamic identifying and monitoring technology based on synthetic aperture radar differential interferometry (D-InSAR)
CN103528489A (en) * 2013-10-22 2014-01-22 安徽理工大学 Monitoring method and device for subsidence of rock strata in overlying bending and sinking zone of goaf
CN105158760A (en) * 2015-08-10 2015-12-16 中南大学 Method for inverting underground fluid volume change and three dimension surface deformation using InSAR

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6911931B2 (en) * 2002-10-24 2005-06-28 The Regents Of The University Of California Using dynamic interferometric synthetic aperature radar (InSAR) to image fast-moving surface waves
CN102645650A (en) * 2012-03-06 2012-08-22 北京北科安地科技发展有限公司 Landslide dynamic identifying and monitoring technology based on synthetic aperture radar differential interferometry (D-InSAR)
CN103528489A (en) * 2013-10-22 2014-01-22 安徽理工大学 Monitoring method and device for subsidence of rock strata in overlying bending and sinking zone of goaf
CN105158760A (en) * 2015-08-10 2015-12-16 中南大学 Method for inverting underground fluid volume change and three dimension surface deformation using InSAR

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HU JUN等: "Inferring three-dimensional surface displacement field by combining SAR interferometric phase and amplitude information of ascending and descending orbits", 《SCIENCE CHINA》 *
王洪友: "基于D-InSAR技术和MAI技术加权获取巴姆三维地震同震形变场", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *

Cited By (36)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106842199A (en) * 2017-01-11 2017-06-13 湖南科技大学 It is a kind of to merge the method that different resolution SAR data monitors Ground Deformation
CN107037428B (en) * 2017-03-27 2019-11-12 中国科学院遥感与数字地球研究所 A method of it improving spaceborne dual station difference InSAR and extracts deformation precision
CN107037428A (en) * 2017-03-27 2017-08-11 中国科学院遥感与数字地球研究所 It is a kind of to improve the method that spaceborne dual station difference InSAR extracts deformation precision
CN107102332A (en) * 2017-05-11 2017-08-29 中南大学 The three-dimensional earth's surface deformation monitoring methods of InSAR based on variance components estimate and strees strain model
CN108921174A (en) * 2018-06-04 2018-11-30 中国矿业大学 Based on the big gradient deformation extracting method in mining area under the constraint of multistage match window
CN108921174B (en) * 2018-06-04 2022-07-05 中国矿业大学 Mining area large gradient deformation extraction method based on multistage matching window constraint
CN110017762A (en) * 2019-03-01 2019-07-16 四川省地质工程勘察院 A method of utilizing Offset-Tracking technical monitoring Wriggle-changing Body deformation
CN110058234A (en) * 2019-05-20 2019-07-26 太原理工大学 A method of it resolving mining area surface and settles three-dimensional deformation
CN110765934A (en) * 2019-10-22 2020-02-07 成都理工大学 Geological disaster identification method based on multi-source data fusion
CN110765934B (en) * 2019-10-22 2021-02-19 成都理工大学 Geological disaster identification method based on multi-source data fusion
CN110702050A (en) * 2019-11-12 2020-01-17 山东科技大学 Surface three-dimensional deformation monitoring method based on InSAR
CN111650579A (en) * 2020-06-12 2020-09-11 中南大学 InSAR mining area three-dimensional deformation estimation method and device for rock migration parameter adaptive acquisition and medium
CN111650579B (en) * 2020-06-12 2022-09-30 中南大学 InSAR mining area three-dimensional deformation estimation method and device for rock migration parameter adaptive acquisition and medium
CN111780660A (en) * 2020-07-13 2020-10-16 内蒙古工业大学 Mining area three-dimensional multi-level deformation optimization method and optimization device
CN111780660B (en) * 2020-07-13 2022-03-15 内蒙古工业大学 Mining area three-dimensional multi-level deformation optimization method and optimization device
CN112284332B (en) * 2020-08-31 2021-10-08 北京四象爱数科技有限公司 High-rise building settlement monitoring result three-dimensional positioning method based on high-resolution INSAR
CN112284332A (en) * 2020-08-31 2021-01-29 北京四象爱数科技有限公司 High-rise building settlement monitoring result three-dimensional positioning method based on high-resolution INSAR
CN112184902A (en) * 2020-09-21 2021-01-05 东华理工大学 Underground mining face inversion method for boundary crossing mining identification
CN112147610B (en) * 2020-09-28 2022-08-05 武汉大学 InSAR interferogram down-sampling method fusing topology analysis technology
CN112147610A (en) * 2020-09-28 2020-12-29 武汉大学 InSAR interferogram down-sampling method fusing topology analysis technology
CN111965646B (en) * 2020-10-20 2021-01-26 航天宏图信息技术股份有限公司 Satellite-borne radar data processing method, device, equipment and storage medium
CN111965646A (en) * 2020-10-20 2020-11-20 航天宏图信息技术股份有限公司 Satellite-borne radar data processing method, device, equipment and storage medium
CN112540370A (en) * 2020-12-08 2021-03-23 鞍钢集团矿业有限公司 InSAR and GNSS fused strip mine slope deformation measurement method
CN112965065A (en) * 2021-02-03 2021-06-15 北京观微科技有限公司 Mining area detection method cooperating with SAR image phase, intensity and coherence information
CN112667970A (en) * 2021-02-04 2021-04-16 淮阴师范学院 Deformation monitoring method for super-huge structure combined by multi-source observation technology
CN112667970B (en) * 2021-02-04 2024-04-12 淮阴师范学院 Deformation monitoring method for oversized structure combined by multi-source observation technology
CN113091598A (en) * 2021-04-06 2021-07-09 中国矿业大学 Method for defining stability grade range of goaf building site by InSAR
CN114236541A (en) * 2021-12-08 2022-03-25 电子科技大学 Large-area earth surface three-dimensional deformation calculation method based on Sentinel-1 satellite SAR image
CN114236541B (en) * 2021-12-08 2023-05-16 电子科技大学 Large-area earth surface three-dimensional deformation calculation method based on Sentinel-1 satellite SAR image
CN115512239A (en) * 2022-10-18 2022-12-23 闽江学院 Satellite platform multi-star sensor attitude self-adaptive fusion method
CN115512239B (en) * 2022-10-18 2023-10-03 闽江学院 Multi-star sensor gesture self-adaptive fusion method for satellite platform
CN116047519A (en) * 2023-03-30 2023-05-02 山东建筑大学 Point selection method based on synthetic aperture radar interferometry technology
CN117213443A (en) * 2023-11-07 2023-12-12 江苏省地质调查研究院 Construction and updating method of ground settlement monitoring network with integration of heaves, earth and depth
CN117213443B (en) * 2023-11-07 2024-03-19 江苏省地质调查研究院 Construction and updating method of ground settlement monitoring network with integration of heaves, earth and depth
CN117970331A (en) * 2024-04-02 2024-05-03 中国电建集团华东勘测设计研究院有限公司 InSAR earth surface deformation monitoring method and system
CN117970331B (en) * 2024-04-02 2024-05-31 中国电建集团华东勘测设计研究院有限公司 InSAR earth surface deformation monitoring method and system

Also Published As

Publication number Publication date
CN106526590B (en) 2018-09-25

Similar Documents

Publication Publication Date Title
CN106526590B (en) A kind of fusion multi-source SAR image industrial and mining area three-dimensional earth's surface deformation monitorings and calculation method
Sansosti et al. Geometrical SAR image registration
Wang et al. Improved SAR amplitude image offset measurements for deriving three-dimensional coseismic displacements
He et al. Complete three-dimensional near-field surface displacements from imaging geodesy techniques applied to the 2016 Kumamoto earthquake
Montazeri et al. Three-dimensional deformation monitoring of urban infrastructure by tomographic SAR using multitrack TerraSAR-X data stacks
CN106226767B (en) Mining area three-D sequential deformation monitoring method based on single radar imagery geometry SAR images
Tiwari et al. Monitoring of landslide activity at the Sirobagarh landslide, Uttarakhand, India, using LiDAR, SAR interferometry and geodetic surveys
CN113238228B (en) Three-dimensional earth surface deformation obtaining method, system and device based on level constraint
Raucoules et al. Combined use of space-borne SAR interferometric techniques and ground-based measurements on a 0.3 km2 subsidence phenomenon
Yang et al. Monitoring urban subsidence with multi-master radar interferometry based on coherent targets
Xing et al. Time series ground subsidence inversion in mining area based on CRInSAR and PSInSAR integration
Spore et al. Collection, processing, and accuracy of mobile terrestrial lidar survey data in the coastal environment
Artese et al. Monitoring of land subsidence in Ravenna municipality using integrated SAR-GPS techniques: Description and first results
Goel Advanced stacking techniques and applications in high resolution SAR interferometry
Yuwono et al. Time series of land subsidence rate on coastal demak using gnss cors udip and dinsar
Ge et al. Integrating corner reflectors and PSInSAR technique to monitor regional land subsidence
Zhu et al. Monitoring of Surface Subsidence of the Mining Area Based on SBAS.
Zamiri Aghdam et al. Monitoring of Urmia Lake Bridge Subsidence during 2014-2021 Using DInSAR-SBAS Method and GPS Data
Du Three dimensional temporal-spatial model based SAR algorithms for monitoring large gradient deformations in mining area
Liu et al. A New Deformation Enhancement Method based on Multi-temporal InSAR for Landslide Surface Stability Assessment
Wu et al. Ground subsidece monitoring along the levee of Huangpu River by multi-platform high-resolution SAR images
Li et al. Land Subsidence Monitoring with Multi-Track SAR Data
Karlsen Measuring velocities of a surge type glacier with SAR interferometry using ALOS-2 data
Zhang et al. Multi-platform persistent scatterer SAR interferometry time series analyzing
Vrankić et al. Praćenje geofizičkih procesa na području zapadne Hercegovine satelitskim opažanjima

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant