CN110146883A - The same period accelerated based on minimum and more satellite platform MT-InSAR three-dimensional deformation decomposition methods of different phase - Google Patents
The same period accelerated based on minimum and more satellite platform MT-InSAR three-dimensional deformation decomposition methods of different phase Download PDFInfo
- Publication number
- CN110146883A CN110146883A CN201910411920.5A CN201910411920A CN110146883A CN 110146883 A CN110146883 A CN 110146883A CN 201910411920 A CN201910411920 A CN 201910411920A CN 110146883 A CN110146883 A CN 110146883A
- Authority
- CN
- China
- Prior art keywords
- deformation
- data set
- date
- data
- satellite
- 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
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/904—SAR modes
- G01S13/9058—Bistatic or multistatic SAR
-
- 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
-
- 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/904—SAR modes
- G01S13/9064—Inverse SAR [ISAR]
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 present invention provides a kind of same period accelerated based on minimum and more satellite platform MT-InSAR three-dimensional deformation decomposition methods of different phase, satellite sight can be decomposed thing, north and south and vertical direction to MT-InSAR deformation timing and rate of deformation.The following steps are included: step 1, using the sight of MT-InSAR difference inverting ascending, descending rail satellite image collection to deformation data;Step 2, the obtained sight of ascending, descending rail satellite platform inverting is extracted to public deformation timing;Step 3, spatial registration and the high coherent point of the same name of ascending, descending rail satellite platform image set is extracted;Step 4, it for the same period, different issue according to the feature mixedly appeared, by ascending, descending rail time arrow mixing and duplicate removal processing, establishes the deformation accelerated based on minimum and decomposes equation group, solution obtains three-dimensional deformation sequence.The present invention, according to the new feature mixedly appeared, develops the multi-platform MT-InSAR data three-dimensional decomposition method for this feature for the new Synthetic Aperture Radar satellite constellation same period, different issue, obtains reliable three-dimensional Ground Deformation timing.
Description
Technical field
The invention belongs to the technical fields of remote sensing and geodesic survey, are related to a kind of flat for the ascending, descending rail of the same period and different phase
Platform MT-InSAR sight, can be by satellite direction of visual lines MT-InSAR deformation timing to the three-dimensional exploded method of Ground Deformation timing
Thing, north and south and vertical direction are decomposed with rate of deformation.This method can for the new Synthetic Aperture Radar satellite constellation same period,
Different issue realizes that multi-platform MT-InSAR data three-dimensional decomposes according to the new feature mixedly appeared, obtains reliable three-dimensional Ground Deformation
Timing.
Background technique
Timing synthetic aperture radar interferometry method (MT-InSAR) can obtain reliable high-precision satellite sight side
To (LOS) deformation timing and rate of deformation, it is widely used in the fields such as landslide, earthquake, bridge, city deformation monitoring.
However it is limited to the imaging geometry of single satellite platform, atural object can only be reflected by the deformation quantity of single platform MT-InSAR inverting
In the deformation of satellite direction of visual lines, and generally require to convert the deformation quantity of satellite direction of visual lines in practice east-west direction,
The deflection of North and South direction and vertical direction needs accurately to obtain vertical if City Building and infrastructure Gernral Check-up
With the deformation result of horizontal direction.
Current multi-platform MT-InSAR deformation timing three-dimensional exploded method, is mostly divided from the angle of mathematic(al) manipulation
Solution lacks and introduces necessary physical condition constraint, and obtained mathematical solution may there is no physical significances, and its steady degree also has
It is to be modified.On the other hand, in orbit along with New Generation Radar satellite, especially equipped with same synthetic aperture radar
Multi-satellite operates in the satellite constellation of same track, the same area can repeatedly be imaged interior on the same day, to the same area
When observe for a long time, it may appear that rail lift and drop rail satellite data, the same period or different issue evidence and the same period and different issue are according to mixed
It is miscellaneous situation occur.This is that Synthetic Aperture Radar satellite constellation data of new generation obtains the new feature occurred, existing three-dimensional deformation
Decomposition method does not have consideration.
Summary of the invention
The new spy that the purpose of the present invention is mixedly appear for the new Synthetic Aperture Radar satellite constellation same period, different issue evidence
Sign develops the multi-platform MT-InSAR data three-dimensional decomposition method for this feature, obtains reliable three-dimensional Ground Deformation timing.
Realizing the specific technical solution of the object of the invention is:
A kind of same period and more satellite platform timing SAR interferometry methods (MT-InSAR) of different phase based on minimum acceleration
Three-dimensional deformation decomposition method, this method is for the new Synthetic Aperture Radar satellite constellation same period, different issue according to the new spy mixedly appeared
Sign can obtain reliable three-dimensional Ground Deformation timing;Specifically includes the following steps:
Step 1: using MT-InSAR difference inverting rail lift and dropping the sight of rail satellite platform to deformation data
Using MT-InSAR difference inverting difference rail lift and the sight of rail Synthetic Aperture Radar satellite platform drops to when deformation
Between sequence and rate of deformation;
Step 2: extracting the obtained sight of ascending, descending rail satellite platform inverting to public deformation timing
The MT-InSAR sight that extraction rail lift and drop rail satellite platform inverting obtain respectively is to deformation timing in common time
Common deformation timing in range, and then rail lift and drop rail satellite platform MT-InSAR sight are extracted to deformation timing and rate field
Public space range, and extract common range accordingly;
Step 3: spatial registration and the high coherent point of the same name for extracting ascending, descending rail satellite platform image set
Spatial registration is carried out to rail lift and drop rail satellite platform the MT-InSAR high coherent point respectively extracted, lays equal stress on and samples
Equal resolution, so that the high coherent point that the MT-InSAR of rail lift and drop rail satellite platform is respectively extracted corresponds;Statistics and
It records rail lift and drops quantity and the position of rail satellite platform high coherent point of the same name;
Step 4: for the same period, different issue according to the feature mixedly appeared, at ascending, descending rail time arrow mixing and duplicate removal
Reason is established the deformation accelerated based on minimum and decomposes equation group, and solution obtains three-dimensional deformation sequence;It specifically includes:
A. after step 1-3 processing, the visual angle of N number of difference Synthetic Aperture Radar satellite data set is different, both comprising rising
Rail data, also comprising drop rail data;Each data set has Qi(i=1,2 ..., N) a time series, each data set direction of visual lines
Rate of deformation described with following formula:
Every a line all represents deformation of the Synthetic Aperture Radar satellite platform data collection on some common point in above formula
The initial timing issue of time series, each platform is different,It is the direction of visual lines shape of first data set
Become time series vector to n-th data set direction of visual lines deformation data vector,It is first
Each phase of the first phase to a last phase accumulates deformation quantity in data set,It is that second data are concentrated from the
Each phase of one phase to a last phase accumulates deformation quantity,Be in n-th data set from the first phase to the end
Each phase of one phase accumulates deformation quantity, Q1It is the time series number of first data set, Q2..., QNIt is second data set to
Time series number in N number of data set;
The acquisition date corresponding to the time series of same each data set indicates are as follows:
T in above formulai(i=1,2 ..., N) is the data acquisition date vector of first data set to n-th data set
Data acquisition date vector,It is the data acquisition day that first data concentrates the first phase to a last phase
Phase,It is the data acquisition date that second data concentrates the first phase to a last phase,It is data acquisition date of the first phase to a last phase, Q in n-th data set1It is first data set
Time series number, Q2..., QNIt is second data set to the time series number in n-th data set;
For Synthetic Aperture Radar satellite constellation, different Synthetic Aperture Radar satellites can be interior more on the same day to the same area
Secondary imaging, i.e. contemporaneous data;The single plateau time sequence of the difference corresponding acquisition date may be identical in above formula, i.e. the same period,
May be different, i.e., the different phase;MT-InSAR deformation timing multi-platform for the different phase, will acquire the time sorts according to sequencing;
MT-InSAR deformation timing multi-platform for the same period will carry out date duplicate removal processing to it, guarantee the total date sequence finally obtained
Column are unique;
B. the same period, different issue first do duplicate checking before being ranked up to all acquisition dates, under this process is used according to appearance is mixed
Formula describes:
T=Sort (Unique (T1, T2..., TN))=[t0, t1, t2..., tQ-1] (3)
T is the not repetition time vector ultimately generated in above formula, and Sort is ranking functions, is used to institute's having time according to elder generation
Sequence ordered arrangement afterwards, Unique is duplicate removal function, for the date repeated is only retained one, T1, T2..., TNIt is
Data acquisition date vector from first data set to n-th data set, t0, t1, t2..., tQ-1It is by sequence and duplicate removal
All data set date sequences retained afterwards, Q are unduplicated date sum in N number of data set;
Obtain N number of data set data acquisition date not repeat date sequence T after, by date list dislocation subtract each other, obtain
Q-1 difference result, is expressed as follows with mathematic(al) representation:
Δ T=[Δ t1, Δ t2..., Δ tQ-1] (4)
Δ T is the difference date intervals that will successively subtract each other before and after the date sequence in T, Δ t in above formula1, Δ
t2..., Δ tQ-1It is first to the Q-1 difference date intervals obtained after subtracting each other date sequence disruption in T, Δ t1It is T
Middle t1-t0As a result, Δ t2It is t in T2-t1As a result, Δ tQ-1It is t in TQ-1-tQ-2As a result, Q be N number of data set in do not repeat
Date sum;
It is several according to satellite imagery since the earth's surface deformation result that MT-InSAR method obtains is the deformation of satellite direction of visual lines
What relationship, decomposes thing, north and south and vertical direction for the deformation of satellite direction of visual lines, describes this decomposable process with following formula:
θ is satellite look angle in above formula,For satellite flight direction course angle, dLOSSatellite direction of visual lines deformation quantity is represented,
dEast-WestFor the deformation quantity of east-west direction, dNorth-southFor North and South direction deformation quantity, dUp-DownFor the deformation quantity of vertical direction;
In addition, the nearly north-south in satellite flight direction, the deformation results of inverting are insensitive to north-south, if not considering the shape of North and South direction
Variation amount is reduced to the expression of two-dimensional deformation decomposition method are as follows: by course angleIt is set as 0, it will be by satellite direction of visual lines shape in above formula
Variable is only decomposed into thing and vertical direction;
C. the different date intervals of Q-1 are shared, by the east-west direction in each date intervals, North and South direction and vertical
Direction average deformation speed is indicated with matrix, obtains following formula:
V in above formulaEast-West, VNorth-south, VUp-Down, represent east-west direction, North and South direction and vertical direction average deformation
Velocity vector, υE1, υE2..., υEQ-1It is average deformation velocity component, υ in each phase time interval in east-west direction Δ TN1,
υN2..., υNQ-1, average deformation velocity component, υ in each phase time interval in North and South direction Δ TU1, υU2..., υUQ-1It is vertical
Average deformation velocity component in each phase time interval in the Δ T of direction, Q are unduplicated date sum in N number of data set;
There is the expression of above formula, can be merged with the averaged deformation rate of each date intervals in conjunction with formula (5)
At total deformation quantity sequence, rate-deformation quantity transfer equation is expressed as following formula:
B is that rate turns deformation quantity reconstruct synthetical matrix, B in above formulai(i=1,2 ..., N) is that rate turns deformation quantity reconstruct
Matrix embodies in formula below (10), VEast-West, VNorth-south, VUp-Down, it is east-west direction, the South and the North respectively
To with vertical direction average deformation velocity vector,It is from when the direction of visual lines deformation of first data set
Between sequence vector to the direction of visual lines deformation data vector of n-th data set, θ is satellite look angle,For satellite flight direction
Course angle;
In order to construct Bi(i=1,2 ..., N), by the time series of each data set corresponding acquisition date after sequence
Time interval sequence Δ T in position be set asIt indicates are as follows:
In above formulaIt is the date vector T of first data set1After sequence and duplicate removal
All data set time vector Ts in position number,It arrivesIt is the date vector T of second data set2All data sets after sequence and duplicate removal
The date vector T of position number in time arrow T to n-th data setNAll data set times after sequence and duplicate removal
Position number in vector T, Qi(i=1,2 ..., N) represents the 1 time series quantity for arriving each data set of N;
There are the help of above formula, BiIt can be constructed with following formula:
B in above formulaiTurn deformation quantity restructuring matrix for rate,It is the position-order in formula (9)
Number, Q is unduplicated date sum in N number of data set;
D. by minimum accelerate assume is introduced into the solution of formula (7), by restriction all directions on deformation component when
There is minimum acceleration, i.e., the speed difference between continuous time interval is minimum, so that solving the equation that deformation is decomposed has in sequence
Physical significance, and obtain stable solution;Foregoing description specifically may be expressed as:
In above formula, C is minimum acceleration defining equation group, and α is regularization factors, for balancing the opposite power of minimum acceleration
Weight, to limit the remaining Norm minimum in formula,Respectively east-west direction,
Average deformation speed difference between North and South direction and vertically adjacent date intervals is as a result, Q is unduplicated in N number of data set
Date sum;
Aggregative formula (7)-(11) obtain formula (12), average for solving east-west direction, North and South direction and vertical direction
Deformation velocity.
In above formula shown in B such as formula (8), shown in C such as formula (11), VEast-West, VNorth-south, VUp-Down, respectively eastern
West to, North and South direction and vertical direction average deformation velocity vector,It is the view from first data set
Direction of visual lines deformation data vector of the line direction deformation data vector to n-th data set;
E. the Vertical Square that above-mentioned system of linear equations obtains each high coherent point is solved with crackle Burger-Ma Kuaer special formula method
To, north and south and east-west direction deformation timing.
F. step a-e is applied in all high coherent points of the same name to obtain in public domain based on ascending, descending rail platform MT-
The three-dimensional Surface Deformation Field and rate of deformation of InSAR.
A kind of same period and more satellite platform MT-InSAR three-dimensional shaped variations of different phase based on minimum acceleration that the present invention develops
Satellite direction of visual lines MT-InSAR deformation timing and rate of deformation can be decomposed thing, north and south and vertical direction by solution method.
This method can be realized multi-platform for the new Synthetic Aperture Radar satellite constellation same period, different issue according to the new feature mixedly appeared
MT-InSAR deformation timing three-dimensional exploded obtains reliable three-dimensional Ground Deformation timing.
The present invention has the advantages that (1) for the new Synthetic Aperture Radar satellite constellation ascending, descending rail satellite data same period,
The different phase mixes susceptible shape with the different phase, can obtain reliable three-dimensional Ground Deformation timing and three-dimensional Surface Deformation Field.(2) algorithm
It is steady and efficient.Algorithm has fully considered this minimum physical constraint of the acceleration of adjacent timing deformation, ensure that side to be solved
Journey group positive definite, and algorithm directly solves deformation timing, efficient stable with matrix operation.
The present invention is suitable for the situation that multi-source radar data obtains direction complementation, i.e., at least to have a pair of of ascending, descending rail data
Collection is to guarantee that three-dimensional deformation decomposition result is reliable and stable.
Detailed description of the invention
Fig. 1 flow chart of the present invention;
Fig. 2 is the S1A ascending, descending rail same period, different phase data obtaining time distribution map;
Fig. 3 is ascending, descending rail Synthetic Aperture Radar satellite platform imaging geometry schematic diagram;
Thing and vertical direction averaged deformation rate diagram after the three-dimensional deformation decomposition of Fig. 4 embodiment of the present invention;
Thing that four characteristic points of Fig. 5 embodiment of the present invention are decomposed and vertically to deformation data figure.
Specific embodiment
Embodiment
Below by taking certain studies area as an example, illustrate the specific implementation step of the present invention in embodiment:
Step 1: using MT-InSAR difference inverting rail lift and dropping the sight of rail satellite platform to deformation data
The rail lift and drop rail data of synthetic aperture radar in a certain research area of search covering, are defended using MT-InSAR inverting difference
Star platform (including rail lift and drop rail) sight is to deformation data.This sentences certain research area Sentinel-1A lifting of acquisition
The present invention is described in detail for orbital data.
The information such as following table of area's fluctuating orbit data is studied, data obtaining time distribution such as attached drawing 2:
Table 1.Sentinel-1A orbital data explanation
Step 2: extracting the obtained sight of ascending, descending rail satellite platform inverting to public deformation timing
The MT-InSAR sight that extraction S1A rail lift and drop rail satellite platform inverting obtain respectively is to deformation timing when public
Between deformation timing common in range, and then extract rail lift and drop rail satellite platform MT-InSAR sight to deformation timing and rate
The public space range of field, and common range is extracted accordingly;Ascending, descending rail Synthetic Aperture Radar satellite platform imaging geometry is such as
Shown in attached drawing 3.
Step 3: spatial registration and the high coherent point of the same name for extracting ascending, descending rail satellite platform image set
Spatial registration is carried out to the MT-InSAR of the S1A rail lift and drop rail satellite platform high coherent point respectively extracted, is laid equal stress on
Equal resolution is sampled, so that the high coherent point that the MT-InSAR of rail lift and drop rail satellite platform is respectively extracted corresponds.
Statistics and record rail lift and quantity and the position for dropping rail satellite platform high coherent point of the same name;
Step 4: establishing the three-dimensional deformation accelerated based on minimum and decompose equation group generation three-dimensional deformation sequence
A. after the pretreatment of first three step, the data of 2 different perspectivess with common time and spatial dimension are obtained
Collection (rail lift and drop rail), the two data sets have 23 and 22 deformation datas respectively, two datasets direction of visual lines
Rate of deformation can be described with following formula:
Every a line respectively represents rail lift and drops deformation data of the rail data set on some common point in above formula,It represents from the direction of visual lines deformation data vector of first data set to the sight side of second data set
To deformation data vector,It represents first data and concentrates and tire out from the first phase to single phase of a last phase
Product deformation quantity,It represents second data and concentrates and accumulate deformation quantity from the first phase to single phase of a last phase;
The acquisition date corresponding to same lift rail time series can indicate are as follows:
T in above formulai(i=1,2) represents the number of data acquisition date second data set of vector sum of first data set
According to obtain date vector,First data is represented to concentrate from the first phase to the data acquisition day of a last phase
Phase,Second data is represented to concentrate from the first phase to the data acquisition date of a last phase;
The single plateau time sequence of the difference corresponding acquisition date is possible identical (same period) in above formula, it is also possible to different (different
Phase).MT-InSAR deformation timing multi-platform for the different phase, it is only necessary to which will acquire the time sorts according to sequencing;For the same period
Multi-platform MT-InSAR deformation timing needs to carry out it date duplicate removal processing, guarantees that the total date sequence finally obtained is unique.
B. in view of the influence that the same period different phase mixes, the present invention has first done before being ranked up to all acquisition dates and has once looked into
Weight, this process can be described with following formula:
T=Sort (Unique (T1, T2))=[t0, t1, t2..., t22] (3)
T represents the not repetition time vector ultimately generated in above formula, and Sort represents ranking functions, for pressing institute's having time
According to sequencing ordered arrangement, Unique represents duplicate removal function, for the date repeated is only retained one, T1, T2It represents
The data acquisition date vector of first data set and second data set, t0, t1, t2..., t22It represents by sequence and duplicate removal
All data set date sequences remained afterwards;
It obtains after not repeating date sequence T, date list dislocation is subtracted each other, 22 difference results can be obtained, use mathematics
Expression formula is expressed as follows:
Δ T=[Δ t1, Δ t2..., Δ t22] (4)
Δ T represents the difference date intervals that will successively subtract each other before and after the date sequence in T, Δ t in above formula1, Δ
t2..., Δ t22First to the 22nd difference date intervals obtained after subtracting each other date sequence disruption in T are represented, such as Δ t1
Represent t in T1-t0As a result, Δ t2Represent t in T2-t1As a result, Δ t22Represent t in T22-t21Result;
Since the earth's surface deformation result that MT-InSAR method obtains only reflects the component of satellite direction of visual lines (LOS), root
According to satellite imagery geometrical relationship attached drawing Fig. 1, the deformation in the direction LOS is decomposed into thing and vertical direction, ignores North and South direction point
Amount, using simplified two-dimensional deformation isolation, this decomposable process can be described with following formula:
dLOS=sin θ * dEast-West+cosθ*dUp-Down (5)
θ is satellite look angle, d in above formulaLOSRepresent satellite direction of visual lines deformation quantity, dEast-WestFor the deformation of east-west direction
Amount, dUp-DownFor the deformation quantity of vertical direction.
C. in view of sharing 22 different date intervals, by the east-west direction and vertical direction in each date intervals
Average deformation speed is indicated with matrix, can obtain following formula:
V in above formulaEast-West, VNorth-south, VUp-Down, represent east-west direction, North and South direction and vertical direction average deformation
Velocity vector, υE1, υE2..., υE22, υU1, υU2..., υU22Represent each phase time interval in east-west direction and vertical direction Δ T
Interior average deformation velocity component;
There is the expression of above formula, can be synthesized with the averaged deformation rate of each date intervals in conjunction with formula (5)
Total deformation quantity sequence, rate-deformation quantity transfer equation can be expressed as following formula:
B is that rate turns deformation quantity reconstruct synthetical matrix, B in above formulai(i=1,2) is that rate turns deformation quantity restructuring matrix,
It embodies in formula below (10), VEast-West, VUp-Down, represent east-west direction and vertical direction average deformation speed to
Amount,It represents from the direction of visual lines deformation data vector of first data set to the view of second data set
Line direction deformation data vector, θ are satellite look angle,For satellite flight direction course angle.
In order to construct Bi, need to find the corresponding acquisition date T of time series of lift rail data set respectivelyi(i=1,2)
The position in time interval sequence Δ T after sequence, is set asIt can indicate are as follows:
In above formulaRepresent the date vector T of first data set1After sequence and duplicate removal
All data set time vector Ts in position number,Represent the date of second data set
Vector T2The position number in all data set time vector Ts after sequence and duplicate removal;
There are the help of above formula, BiIt can be constructed with following formula:
B in above formulaiTurn deformation quantity restructuring matrix for rate,Represent the position number in formula (9);
D. arrive system of linear equations most basic here and have been built up and complete, as long as solve formula (7) as a result,
The Deformation partition of east-west direction and vertical direction is obtained as a result, the present invention assumes minimum acceleration to introduce the solution of above-mentioned matrix
In, there is minimum acceleration, i.e. speed between continuous time interval in timing by limiting the deformation component in all directions
Degree difference is minimum, so that solving the equation that deformation is decomposed has physical significance, and can obtain stable solution.Foregoing description specifically can table
It is shown as:
In above formula, C is minimum acceleration defining equation group, and α is regularization factors, for balancing the opposite power of minimum acceleration
Weight, to limit the remaining Norm minimum in formula,Represent east-west direction and vertically adjacent
Average deformation speed difference result between date intervals;
To sum up the complete same period and the multi-platform MT-InSAR three-dimensional deformation decomposition technique of different phase ascending, descending rail can indicate are as follows:
In above formula shown in B such as formula (8), shown in C such as formula (11), VEast-West, VUp-Down, represent east-west direction and hang down
Histogram to average deformation velocity vector,Represent from the direction of visual lines deformation data of first data set to
Measure the direction of visual lines deformation data vector to second data set;
E. the Vertical Square that above-mentioned system of linear equations obtains each high coherent point is solved with crackle Burger-Ma Kuaer special formula method
To, north and south and east-west direction deformation timing.
F. step a-e is applied in all high coherent points of the same name to obtain in public domain based on ascending, descending rail platform MT-
The three-dimensional Surface Deformation Field and rate of deformation of InSAR.Three-dimensional deformation in the present embodiment decompose after thing and vertical direction it is flat
Equal rate of deformation figure is as shown in Fig. 4, and obtained thing will be decomposed and vertically draw to deformation data by choosing four characteristic points
It produces and, as shown in Fig. 5.
The content that this specification is not described in detail belongs to the prior art well known to professional and technical personnel in the field.
Claims (1)
1. a kind of same period accelerated based on minimum and more satellite platform MT-InSAR three-dimensional deformation decomposition methods of different phase, feature are existed
In, this method comprising the following specific steps
Step 1: using MT-InSAR difference inverting rail lift and dropping the sight of rail satellite platform to deformation data
Using MT-InSAR difference inverting difference rail lift and the sight of rail Synthetic Aperture Radar satellite platform drops to deformation time sequence
Column and rate of deformation;
Step 2: extracting the obtained sight of ascending, descending rail satellite platform inverting to public deformation timing
The MT-InSAR sight that extraction rail lift and drop rail satellite platform inverting obtain respectively is to deformation timing in common time range
Interior common deformation timing, and then extract rail lift and drop public affairs of the rail satellite platform MT-InSAR sight to deformation timing and rate field
Cospace range, and common range is extracted accordingly;
Step 3: spatial registration and the high coherent point of the same name for extracting ascending, descending rail satellite platform image set
Spatial registration is carried out to rail lift and drop rail satellite platform the MT-InSAR high coherent point respectively extracted, laying equal stress on, it is identical to sample
Resolution ratio, so that the high coherent point that the MT-InSAR of rail lift and drop rail satellite platform is respectively extracted corresponds;Statistics and record
The quantity and position of rail lift and drop rail satellite platform high coherent point of the same name;
Step 4: building ascending, descending rail time arrow mixing and duplicate removal processing according to the feature mixedly appeared for the same period, different issue
Equation group is decomposed in the minimum deformation accelerated that is based on, and solution obtains three-dimensional deformation sequence;It specifically includes:
A. after step 1-3 processing, it both includes rail lift number that the visual angle of N number of difference Synthetic Aperture Radar satellite data set is different
According to also comprising drop rail data;Each data set has Qi(i=1,2 ..., N) a time series, the shape of each data set direction of visual lines
Variable Rate is described with following formula:
Every a line all represents deformation time of the Synthetic Aperture Radar satellite platform data collection on some common point in above formula
The initial timing issue of sequence, each platform is different,When being the direction of visual lines deformation of first data set
Between sequence vector to n-th data set direction of visual lines deformation data vector,It is first number
Deformation quantity is accumulated according to each phase of the concentration first phase to a last phase,It is that second data is concentrated from first
Each phase of phase to a last phase accumulates deformation quantity,It is in n-th data set from the first phase to last one
Each phase of phase accumulates deformation quantity, Q1It is the time series number of first data set, Q2..., QNIt is second data set to N
The time series number that a data are concentrated;
The acquisition date corresponding to the time series of same each data set indicates are as follows:
T in above formulai(i=1,2 ..., N) is the data acquisition date vector of first data set to the data of n-th data set
Date vector is obtained,It is the data acquisition date that first data concentrates the first phase to a last phase,It is the data acquisition date that second data concentrates the first phase to a last phase,
It is data acquisition date of the first phase to a last phase, Q in n-th data set1It is the time series number of first data set,
Q2..., QNIt is second data set to the time series number in n-th data set;
For Synthetic Aperture Radar satellite constellation, different Synthetic Aperture Radar satellites to the same area can on the same day repeatedly at
Picture, i.e. contemporaneous data;The single plateau time sequence of the difference corresponding acquisition date may be identical in above formula, i.e. the same period, it is also possible to
Difference, i.e., different phase;MT-InSAR deformation timing multi-platform for the different phase, will acquire the time sorts according to sequencing;For
The same period multi-platform MT-InSAR deformation timing will carry out date duplicate removal processing to it, guarantee the total date sequence finally obtained only
One;
B. the same period, different issue be according to mixing appearances, to all acquisition dates be ranked up it is preceding first do duplicate checking, this process with following formula come
Description:
T=Sort (Unique (T1, T2..., TN))=[t0, t1, t2..., tQ-1] (3)
T is the not repetition time vector ultimately generated in above formula, and Sort is ranking functions, is used to institute's having time according to successively suitable
Sequence ordered arrangement, Unique is duplicate removal function, for the date repeated is only retained one, T1, T2..., TNIt is from
Data acquisition date vector of one data set to n-th data set, t0, t1, t2..., tQ-1It is by sorting and being protected after duplicate removal
All data set date sequences stayed, Q are unduplicated date sum in N number of data set;
Obtain N number of data set data acquisition date do not repeat date sequence T after, date list dislocation is subtracted each other, Q-1 is obtained
Difference result is expressed as follows with mathematic(al) representation:
Δ T=[Δ t1, Δ t2..., Δ tQ-1] (4)
Δ T is the difference date intervals that will successively subtract each other before and after the date sequence in T, Δ t in above formula1, Δ t2..., Δ
tQ-1It is first to the Q-1 difference date intervals obtained after subtracting each other date sequence disruption in T, Δ t1It is t in T1-t0's
As a result, Δ t2It is t in T2-t1As a result, Δ tQ-1It is t in TQ-1-tQ-2As a result, Q be N number of data set in the unduplicated date it is total
Number;
Since the earth's surface deformation result that MT-InSAR method obtains is the deformation of satellite direction of visual lines, closed according to satellite imagery geometry
System, decomposes thing, north and south and vertical direction for the deformation of satellite direction of visual lines, describes this decomposable process with following formula:
θ is satellite look angle in above formula,For satellite flight direction course angle, dLOsRepresent satellite direction of visual lines deformation quantity, dEast-West
For the deformation quantity of east-west direction, dNorth-SouthFor North and South direction deformation quantity, dUp-DownFor the deformation quantity of vertical direction;In addition, defending
The nearly north-south of star heading, the deformation results of inverting are insensitive to north-south, if not considering the deformation component of North and South direction,
It is reduced to the expression of two-dimensional deformation decomposition method are as follows: by course angleIt is set as 0, it will be by satellite direction of visual lines deformation quantity only in above formula
It is decomposed into thing and vertical direction;
C. Q-1 different date intervals are shared, by east-west direction, North and South direction and the vertical direction in each date intervals
Average deformation speed is indicated with matrix, obtains following formula:
V in above formulaEast-West,VUp-Down, represent east-west direction, North and South direction and vertical direction average deformation speed
Spend vector, vE1, vE2..., vEQ-1It is average deformation velocity component, v in each phase time interval in east-west direction Δ TN1,
vN2..., vNQ-1, average deformation velocity component, v in each phase time interval in North and South direction Δ TU1, vU2..., vUQ-1It is vertical
Average deformation velocity component in each phase time interval in the Δ T of direction, Q are unduplicated date sum in N number of data set;
There is the expression of above formula, can be merged into always with the averaged deformation rate of each date intervals in conjunction with formula (5)
Deformation quantity sequence, rate-deformation quantity transfer equation is expressed as following formula:
B is that rate turns deformation quantity reconstruct synthetical matrix, B in above formulai(i=1,2 ..., N) is that rate turns deformation quantity restructuring matrix,
It is embodied in formula below (10), VEast-West,VUp-Down, it is east-west direction, North and South direction respectively
With vertical direction average deformation velocity vector,It is from when the direction of visual lines deformation of first data set
Between sequence vector to the direction of visual lines deformation data vector of n-th data set, θ is satellite look angle,For satellite flight direction
Course angle;
In order to construct Bi(i=1,2 ..., N), by the time series of each data set it is corresponding obtain the date after sequence when
Between position in intervening sequence Δ T be set asIt indicates are as follows:
In above formulaIt is the date vector T of first data set1Institute after sequence and duplicate removal
There is the position number in data set time vector T,It arrivesIt is the date vector T of second data set2All data sets after sequence and duplicate removal
The date vector T of position number in time arrow T to n-th data setNAll data set times after sequence and duplicate removal
Position number in vector T, Qi(i=1,2 ..., N) represents the 1 time series quantity for arriving each data set of N;
There are the help of above formula, BiIt is constructed with following formula:
B in above formulaiTurn deformation quantity restructuring matrix for rate,It is the position number in formula (9), Q
For date sum unduplicated in N number of data set;
D. minimum is accelerated to assume to be introduced into the solution of formula (7), by limiting the deformation component in all directions in timing
With minimum acceleration, i.e., the speed difference between continuous time interval is minimum, so that solving the equation that deformation is decomposed has physics
Meaning, and obtain stable solution;Foregoing description specifically may be expressed as:
In above formula, C is minimum acceleration defining equation group, and α is regularization factors, for balancing the relative weighting of minimum acceleration,
To limit the remaining Norm minimum in formula,Respectively east-west direction, north and south
Average deformation speed difference between direction and vertically adjacent date intervals is as a result, Q is the unduplicated date in N number of data set
Sum;
Aggregative formula (7)-(11) obtain formula (12), for solving east-west direction, North and South direction and vertical direction average deformation
Speed;
In above formula shown in B such as formula (8), shown in C such as formula (11), VEast-West,VUp-Down, respectively thing
Direction, North and South direction and vertical direction average deformation velocity vector,It is the view from first data set
Direction of visual lines deformation data vector of the line direction deformation data vector to n-th data set;
E. above-mentioned system of linear equations is solved with crackle Burger-Ma Kuaer special formula method obtain the vertical direction of each high coherent point, south
North and east-west direction deformation timing;
F. step a-e is applied in all high coherent points of the same name to obtain in public domain based on ascending, descending rail platform MT-
The three-dimensional Surface Deformation Field and rate of deformation of InSAR.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910411920.5A CN110146883B (en) | 2019-05-17 | 2019-05-17 | Simultaneous and asynchronous multi-satellite platform MT-InSAR three-dimensional deformation decomposition method based on minimum acceleration |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910411920.5A CN110146883B (en) | 2019-05-17 | 2019-05-17 | Simultaneous and asynchronous multi-satellite platform MT-InSAR three-dimensional deformation decomposition method based on minimum acceleration |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110146883A true CN110146883A (en) | 2019-08-20 |
CN110146883B CN110146883B (en) | 2022-04-05 |
Family
ID=67595465
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910411920.5A Active CN110146883B (en) | 2019-05-17 | 2019-05-17 | Simultaneous and asynchronous multi-satellite platform MT-InSAR three-dimensional deformation decomposition method based on minimum acceleration |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110146883B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111458709A (en) * | 2020-06-08 | 2020-07-28 | 河南大学 | Satellite-borne radar wide-area earth surface two-dimensional deformation field monitoring method and device |
CN113091598A (en) * | 2021-04-06 | 2021-07-09 | 中国矿业大学 | Method for defining stability grade range of goaf building site by InSAR |
CN113281748A (en) * | 2021-05-24 | 2021-08-20 | 西南石油大学 | Surface deformation monitoring method |
EP4361678A1 (en) * | 2022-10-27 | 2024-05-01 | Planetek Italia S.r.l. | System for determining the components of the displacement of the persistent scatterers and method thereof |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150042510A1 (en) * | 2012-11-07 | 2015-02-12 | Neva Ridge Technologies | SAR Point Cloud Generation System |
CN108983232A (en) * | 2018-06-07 | 2018-12-11 | 中南大学 | A kind of InSAR two dimension earth's surface deformation monitoring method based on adjacent rail data |
CN109540159A (en) * | 2018-10-11 | 2019-03-29 | 同济大学 | A kind of quick complete automatic Pilot method for planning track |
-
2019
- 2019-05-17 CN CN201910411920.5A patent/CN110146883B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150042510A1 (en) * | 2012-11-07 | 2015-02-12 | Neva Ridge Technologies | SAR Point Cloud Generation System |
CN108983232A (en) * | 2018-06-07 | 2018-12-11 | 中南大学 | A kind of InSAR two dimension earth's surface deformation monitoring method based on adjacent rail data |
CN109540159A (en) * | 2018-10-11 | 2019-03-29 | 同济大学 | A kind of quick complete automatic Pilot method for planning track |
Non-Patent Citations (5)
Title |
---|
GUANYU MA ET AL.: "On the Effects of InSAR Temporal Decorrelation and Its Implications for Land Cover Classification:The Case of the Ocean-Reclaimed Lands of the Shanghai Megacity", 《SENSORS》 * |
QIANG WANG ET AL.: "Investigation of the ground displacement in Saint Petersburg, Russia, using multiple-track differential synthetic aperture radar interferometry", 《INT J APPL EARTH OBS GEOINFORMATION》 * |
刘国祥等: "基于多卫星平台永久散射体雷达干涉提取三维地表形变速度场", 《地球物理学报》 * |
杨红磊等: "基于DInSAR和MAI技术揭示地震三维形变场", 《地球物理学进展》 * |
禹雷: "基于多级联合策略的多平台 MT-InSAR 形变时间序列联合分析—以上海市新成陆区为例", 《中国知网》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111458709A (en) * | 2020-06-08 | 2020-07-28 | 河南大学 | Satellite-borne radar wide-area earth surface two-dimensional deformation field monitoring method and device |
CN111458709B (en) * | 2020-06-08 | 2023-12-22 | 河南大学 | Method and device for monitoring wide-area earth surface two-dimensional deformation field of spaceborne radar |
CN113091598A (en) * | 2021-04-06 | 2021-07-09 | 中国矿业大学 | Method for defining stability grade range of goaf building site by InSAR |
CN113281748A (en) * | 2021-05-24 | 2021-08-20 | 西南石油大学 | Surface deformation monitoring method |
CN113281748B (en) * | 2021-05-24 | 2022-02-11 | 西南石油大学 | Surface deformation monitoring method |
EP4361678A1 (en) * | 2022-10-27 | 2024-05-01 | Planetek Italia S.r.l. | System for determining the components of the displacement of the persistent scatterers and method thereof |
Also Published As
Publication number | Publication date |
---|---|
CN110146883B (en) | 2022-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110146883A (en) | The same period accelerated based on minimum and more satellite platform MT-InSAR three-dimensional deformation decomposition methods of different phase | |
Hiesinger et al. | Ages and stratigraphy of lunar mare basalts in Mare Frigoris and other nearside maria based on crater size‐frequency distribution measurements | |
Haxby et al. | Digital images of combined oceanic and continental data sets and their use in tectonic studies | |
Diener et al. | Proto-groups at 1.8< z< 3 in the zCOSMOS-deep sample | |
Bleem et al. | A New Reduction of the Blanco Cosmology Survey: An Optically Selected Galaxy Cluster Catalog and a Public Release of Optical Data Products | |
Bickel et al. | Automated detection of lunar rockfalls using a convolutional neural network | |
Conn et al. | Slicing the monoceros overdensity with SUPRIME-CAM | |
Gladstone et al. | The Lyman‐α Sky Background as Observed by New Horizons | |
CN112505696B (en) | Strip mine side slope deformation dynamic monitoring method based on space-borne SAR | |
Bickel et al. | A labeled image dataset for deep learning-driven rockfall detection on the Moon and Mars | |
Bisi et al. | Three‐dimensional reconstructions of the early November 2004 Coordinated Data Analysis Workshop geomagnetic storms: Analyses of STELab IPS speed and SMEI density data | |
Chabrillat et al. | EnMAP pre-launch and start phase: mission update | |
Neukum | Meteorite bombardment and dating of planetary surfaces | |
Raynal et al. | Thomas Quarry I, Morocco | |
Riedel et al. | The Solar Neighborhood. XLIII. Discovery of New Nearby Stars with μ< 0.″18 yr− 1 (TINYMO Sample) | |
Trilling et al. | The DECam Ecliptic Exploration Project (DEEP). I. Survey Description, Science Questions, and Technical Demonstration | |
McNeill et al. | Asteroid Photometry from the Transiting Exoplanet Survey Satellite: A Pilot Study | |
Rutkofske et al. | Applications of modern digital mapping systems to Assist inclusion of persons with disabilities in geoscience education and research | |
Mathews et al. | Examining Urban Built‐up Volume: Three‐Dimensional Analyses with Lidar and Radar Data | |
Reeves | Aeromagnetic and gravity features of Gondwana and theirrelation to continental break-up: More pieces, less puzzle | |
Hood | Bulk magnetization properties of the Fra Mauro and Reiner Gamma formations | |
Zhang et al. | Surface Subsidence Monitoring of Mining Areas in Hunan Province Based on Sentinel-1A and DS-InSAR | |
Florio et al. | Magnetic anomalies of the tyrrhenian sea revisited: a processing workflow for enhancing the resolution of aeromagnetic data | |
CN109828272A (en) | A kind of SAR difference chromatography method and device | |
Hatch et al. | Application of a new method for calculation of low-frequency wave vectors |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 |