CN109031301A - Alpine terrain deformation extracting method based on PSInSAR technology - Google Patents

Alpine terrain deformation extracting method based on PSInSAR technology Download PDF

Info

Publication number
CN109031301A
CN109031301A CN201811123711.2A CN201811123711A CN109031301A CN 109031301 A CN109031301 A CN 109031301A CN 201811123711 A CN201811123711 A CN 201811123711A CN 109031301 A CN109031301 A CN 109031301A
Authority
CN
China
Prior art keywords
phase
differential
deformation
image
interferogram
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.)
Withdrawn
Application number
CN201811123711.2A
Other languages
Chinese (zh)
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.)
Electric Power Research Institute of Yunnan Power System Ltd
Original Assignee
Electric Power Research Institute of Yunnan Power System Ltd
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 Electric Power Research Institute of Yunnan Power System Ltd filed Critical Electric Power Research Institute of Yunnan Power System Ltd
Priority to CN201811123711.2A priority Critical patent/CN109031301A/en
Publication of CN109031301A publication Critical patent/CN109031301A/en
Withdrawn legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021SAR image post-processing techniques
    • G01S13/9023SAR image post-processing techniques combined with interferometric techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/904SAR modes

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)
  • Image Processing (AREA)

Abstract

本发明公开一种基于PSInSAR技术的山区地形形变提取方法,通过星载卫星获取多幅高分辨率SAR图像,并选取n幅SAR图像,采用综合相干系数法选择其中一幅作主图像,其余作辅图像,再让主图像分别与辅图像配准并进行干涉处理,获得相干系数图和干涉相位图;根据相干系数图,采用相应的振幅离差和相干系数双重阈值法选取PS点并得到PS点三角网;对干涉相位图进行去平滤波,并利用外部DEM数据,对干涉图进行差分干涉处理,得到差分干涉图;分别利用网络平差法、伪3D相位解缠法和Lambda算法三种三维不规则网络相位解缠法对差分干涉图进行相位解缠,得到最后的形变场。

The invention discloses a method for extracting terrain deformation in mountainous areas based on PSInSAR technology. A plurality of high-resolution SAR images are obtained through spaceborne satellites, and n SAR images are selected, and one of them is selected as the main image by a comprehensive coherence coefficient method, and the rest are used as the main image. Then let the main image be registered with the auxiliary image and perform interference processing to obtain the coherence coefficient map and the interferometric phase map; according to the coherence coefficient map, use the corresponding amplitude deviation and coherence coefficient double threshold method to select PS points and get PS Point triangulation; perform flattening and filtering on the interferometric phase map, and use external DEM data to perform differential interferometric processing on the interferogram to obtain a differential interferogram; use network adjustment method, pseudo 3D phase unwrapping method and Lambda algorithm respectively The three-dimensional irregular network phase unwrapping method unwraps the phase of the differential interferogram to obtain the final deformation field.

Description

Alpine terrain deformation extracting method based on PSInSAR technology
Technical field
The present invention relates to synthetic aperture radar image-forming technical field more particularly to a kind of mountain areas based on PSInSAR technology Landform deformation extracting method.
Background technique
Synthetic aperture radar (SAR) has the advantages that round-the-clock, round-the-clock, is widely applied to sees over the ground both at home and abroad Survey field, interfering synthetic aperture radar (InSAR) application that wherein eighties proposes are even more important.
PSInSAR (Permanent scatterers synthetic aperture radar interferometry) detection Ground Deformation is that (difference is dry in DInSAR Relate to measurement) on the basis of develop, the differential interferometry phase of PSInSAR is obtained using DInSAR method.DInSAR is visited Surveying Ground Deformation is to survey on the basis of height to develop from InSAR again, by obtaining the complex pattern pair of ground areal, according to The geometrical relationship of target and aerial position forms interference line figure, interferes and contain in line figure in complex pattern to upper generation phase difference The precise information of the difference of the upward point of oblique distance and two aerial positions, utilizes the orbit parameter, system parameter and base of observation platform Geometrical relationship between line length can measure the three-dimensional position of every bit on image with precise volume, then by redundant observation or outside Portion's terrain data rejects the influence of topography from interferometric phase, to obtain the phase of Ground Deformation, and then is finally inversed by Ground Deformation. For the influence of conventional DInSAR phase misalignment correlation and atmosphere delay, Ferretti, which is proposed, only only keeps track thunder in imaging region The method that those dephasings close serious resolution cell, these target (such as ground are abandoned up to the relatively stable target of scattering properties The corner or roof of building, it is also possible to exposed rock) influence of noise hardly is closed by dephasing, even if in time many years The interference at interval still maintains higher interference correlation in, these stable targets are referred to as Permanent scatterers (PS). Since Permanent scatterers can keep high relevant in some time is spaced, and Space Baseline away from be more than Critical baseline away from feelings Under condition, be also able to maintain high coherence, can make full use of in this way Long baselines away from interference image to improving number to the maximum extent According to utilization rate by carrying out time series analysis to these PS points, eliminate atmosphere to find out the PS point in survey region It influences, just can accurately measure the deformation quantity of PS point, to monitor the movement on ground, and accurately reflect monitored region Relative displacement.
PSInSAR method is similar to control measurement, it obtains the information of whole region by the authentic communication on point, even if Interference fringe cannot be formed in entire research area, can also use the deformation of PSInSAR method detection earth's surface.
Summary of the invention
The present invention provides a kind of alpine terrain deformation extracting method based on PSInSAR technology, with solve the prior art without Method accurately extracts the problem of alpine terrain deformation.
A kind of alpine terrain deformation extracting method based on PSInSAR technology provided by the invention, comprising:
SAR image is chosen, the SAR image includes master image and auxiliary image;
Using the master image and the auxiliary image registration and interference processing is carried out, obtains coherence factor figure and interferometric phase Figure;
PS point is chosen on the coherence factor figure using amplitude deviation and the dual thresholds method of coherence factor, obtains PS point The triangulation network;
The interferometric phase image is subjected to flat filtering processing, using dem data to treated the interferometric phase image Differential interferometry processing is carried out, differential interferometry figure is obtained;
Phase unwrapping processing is carried out to the differential interferometry figure, obtains Deformation Field.
In one embodiment of the invention, SAR image is chosen to specifically include:
N width SAR image is chosen, selects the width in the SAR image as the master map using comprehensive coherence factor method Picture, remaining n-1 width is as the auxiliary image.
In one embodiment of the invention, phase unwrapping processing is carried out to the differential interferometry figure to specifically include:
It is utilized respectively network adjustment method, puppet 3D phase unwrapping method and Lambda algorithm and phase is carried out to the differential interferometry figure Solution twines processing.
In one embodiment of the invention, the network adjustment method specifically includes:
The differential mode pattern three times of PS point is established, the coherence factor K in the differential mode pattern three times is calculated1,K2,B, institute State differential mode pattern three times are as follows:
Plumb line method of loci is introduced, Δ δ ' is chosenHSearch range and step-length with Δ v', by between the adjacent PS point of control The suitable Δ δ ' of multiple correlation coefficient γ search one by oneHWith Δ v', wherein
It regard linear deformation speed difference corresponding to the maximum value of the coherence factor and elevation amendment difference as parameter Estimation Value, establishes adjustment function model, and the adjustment function model is,
Solve linear equation:
In one embodiment of the invention, the puppet 3D phase unwrapping method specifically includes:
Delaunay triangulation network connection is done to the PS point and establishes the triangulation network, time dimension is connected by time relationship;
Residue points calculating is carried out along the side of the triangulation network;
The optimization initial value that space 2D phase unwrapping is done using time dimension disentanglement fruit is carried out using minimum Lp norm method The optimal solution of space dimension.
In one embodiment of the invention, the Lambda algorithm specifically includes:
In n-1 width differential interferometry figure, the phase type of differential interferometry three times of the adjacent PS point of the triangulation network is formed into equation group:
It is denoted as y1=A1·a+B1·b;
In y1=A1·a+B1Imitation observation equation y is added in b2=A2·a+B2B constitutes y=Aa+Bb, wherein A2For 2*m null matrix, B2For 2*2 unit matrix;
It is acquired according to the principle of least square:
C=(A B)
It obtainsComplete cycle float solution, search for obtain the fixed solution of phase complete cycle using fuzziness decorrelation method
By the fixed solutionIt brings solution into and twines phase expression formulaIn reuse minimum two Multiplication obtains:
After completing phase unwrapping, the Deformation Field is obtained by mutually high conversion and geocoding.
A kind of alpine terrain deformation extracting method based on PSInSAR technology provided by the invention, is obtained by piggyback satellite Several High Resolution SAR Images are taken, and choose n width SAR image, wherein a width makees master map using comprehensive coherence factor method selection Picture, remaining makees auxiliary image, then allow master image respectively with auxiliary image registration and carry out interference processing, obtain coherence factor figure and interference Phase diagram;According to coherence factor figure, PS point is chosen using corresponding amplitude deviation and coherence factor dual thresholds method and obtains PS The point triangulation network;Flat filtering is carried out to interferometric phase image, and utilizes external dem data, differential interferometry processing is carried out to interference pattern, Obtain differential interferometry figure;It is utilized respectively three kinds of network adjustment method, puppet 3D phase unwrapping method and Lambda algorithm irregular nets of three-dimensional Network phase unwrapping method carries out phase unwrapping to differential interferometry figure, obtains Deformation Field to the end.
Detailed description of the invention
Fig. 1 is the flow chart of the alpine terrain deformation extracting method provided in an embodiment of the present invention based on PSInSAR technology;
Fig. 2 is the remote sensing image of the Mao County test block in the embodiment of the present invention;
Fig. 3 is the Mao County triangulation diagram in the embodiment of the present invention;
Fig. 4 is the Erlongshan Mountains differential interferometry figure in the embodiment of the present invention;
Fig. 5 is the Chengdu differential interferometry figure in the embodiment of the present invention.
Specific embodiment
Present invention will be further explained below with reference to the attached drawings and examples.
It is the stream of the alpine terrain deformation extracting method provided in an embodiment of the present invention based on PSInSAR technology referring to Fig. 1 Cheng Tu, this method comprises:
S100: SAR image is chosen, the SAR image includes master image and auxiliary image.
SAR image is chosen to specifically include:
N width SAR image is chosen, selects the width in the SAR image as the master map using comprehensive coherence factor method Picture, remaining n-1 width is as the auxiliary image.
S200: using the master image and the auxiliary image registration and carrying out interference processing, obtains coherence factor figure and does Relate to phase diagram.
S300: PS point is chosen on the coherence factor figure using amplitude deviation and the dual thresholds method of coherence factor, is obtained To the PS point triangulation network.
S400: carrying out flat filtering processing for the interferometric phase image, using dem data to treated the interference phase Bitmap carries out differential interferometry processing, obtains differential interferometry figure.
Dem data is digital elevation model (Digital Elevation Model), is by limited landform altitude number The factually existing digitized simulation (i.e. the digital expression of topographical surface form) to ground surface or terrain, is with one group of orderly array of values A kind of actual ground model of form expression ground elevation.In the embodiment of the present invention, using external dem data, to going flat filtering Treated, and interferometric phase image carries out differential interferometry processing, obtains differential interferometry figure.
S500: phase unwrapping processing is carried out to the differential interferometry figure, obtains Deformation Field.
It is utilized respectively network adjustment method, puppet 3D phase unwrapping method and three kinds of Lambda algorithm three-dimensional scramble network phase solutions Quick pushing manipulation carries out phase unwrapping to differential interferometry figure, obtains Deformation Field to the end.
Wherein, network adjustment method specifically includes:
The differential mode pattern three times of PS point is established, the coherence factor K in the differential mode pattern three times is calculated1,K2,B, institute State differential mode pattern three times are as follows:
Plumb line method of loci is introduced, Δ δ ' is chosenHSearch range and step-length with Δ v', by between the adjacent PS point of control The suitable Δ δ ' of multiple correlation coefficient γ search one by oneHWith Δ v', wherein
It regard linear deformation speed difference corresponding to the maximum value of the coherence factor and elevation amendment difference as parameter Estimation Value, establishes adjustment function model, and the adjustment function model is,
Solve linear equation:
Pseudo- 3D phase unwrapping method specifically includes:
Delaunay triangulation network connection is done to the PS point and establishes the triangulation network, time dimension is connected by time relationship.
Residue points calculating is carried out along the side of the triangulation network.Residue points necessarily occur in some triangle, corresponding triangle Certain point in the dual network of net.
The optimization initial value that space 2D phase unwrapping is done using time dimension disentanglement fruit is carried out using minimum Lp norm method The optimal solution of space dimension.
Lambda algorithm is a kind of searching algorithm based on least-squares estimation adjustment, is currently widely applied be based on GPS carrier phase is positioned, and the algorithm of appearance and orientation is surveyed.
Lambda algorithm specifically includes:
In n-1 width differential interferometry figure, the phase type of differential interferometry three times of the adjacent PS point of the triangulation network is formed into equation group:
It is denoted as y1=A1·a+B1·b。
Due to Matrix Properties, equation group intangibility, therefore, in y1=A1·a+B1Imitation observation equation y is added in b2= A2·a+B2B constitutes y=Aa+Bb, whereinA2For 2*m null matrix, B2For 2*2 unit matrix.
It is acquired according to the principle of least square:
C=(A B)
It obtainsComplete cycle float solution, search for obtain the fixed solution of phase complete cycle using fuzziness decorrelation method
By the fixed solutionIt brings solution into and twines phase expression formulaIn reuse minimum two Multiplication obtains:
After completing phase unwrapping, the Deformation Field is obtained by mutually high conversion and geocoding.
The embodiment of the present invention is described in detail below with reference to Fig. 2, Fig. 3, Fig. 4 and Fig. 5.
Survey region is Mao County, and Erlongshan Mountains is located in river blueness block handover region, under the influence of plate motion and rainfall, mudstone The geological disasters such as stream take place frequently, and Chengdu surface configuration changes climate and geological conditions influence is smaller, are affected by culture, The embodiment of the present invention uses the TerraSAR-X image of precision 1m, Mao County, Erlongshan Mountains and Chengdu image information and master map such as 1 institute of table Show.
1 test block image information of table and master map
Place Initial time Terminal time Master map Quantity
Mao County 16/04/23 17/12/19 16/06/17 18
Erlongshan Mountains 14/04/04 15/04/02 14/11/21 15
Chengdu 14/05/20 16/09/02 15/01/17 14
Fig. 2 and Fig. 3 illustrates the PS point triangular network in Mao County trial zone and Mao County area, includes apparent 7 in region Electric power pylon, PS point screening technique use amplitude deviation and coherence factor Integrated Selection, obvious striped in the differential interferometry figure of Mao County Appear in steel tower, rock and ridge region.
Fig. 4 and Fig. 5 illustrates the differential interferometry figure of Erlongshan Mountains and In Chengdu, in the region of Erlongshan Mountains it is affected by noise compared with Greatly, interference fringe is only present in steel tower part, and Chengdu region is typical Plain urban area, in left side culture's the commercial house Area, the middle of the road line take out existing obvious interference fringe.
Network adjustment method, puppet 3D phase unwrapping method and Lambda are used to the differential interferometry figure in Mao County, Erlongshan Mountains and Chengdu Algorithm carries out solution and twines processing, wherein Mao County region linear deformation amount and elevation correction value are as shown in tables 2 and 3.
2 Mao County PS point Linear deformation rate of table is distributed mm/y
v≤-2 -2<v≤-1 -1<v≤0 0<v≤1 1<v≤2 v>2
Network adjustment 0.0313 0.247 0.2615 0.2206 0.2181 0.0306
Pseudo- 3D 0.0218 0.2517 0.2664 0.2254 0.2224 0.0215
Lambda 0 0.1889 0.3510 0.2891 0.1801 0
3 Mao County elevation correction value m of table
h≤-15 -15<h≤-8 -8<h≤0 0<h≤8 8<h≤15 h>15
Network adjustment 0 0.0367 0.4145 0.5073 0.0415 0
Lambda 0 0.1984 0.2653 0.3179 0.2184 0
As seen from the table, Mao County rate of deformation is uniformly distributed, and PS point Ground Deformation performance lifting or downward trend are unified, tower Body region has obvious rise to sink to changing, and elevation correction value deviation trend is unified, and both of which can get trusted elevation knot Fruit.
Erlongshan Mountains regional deformation rate concentrates in -5~+7mm/y range, and 80% or more is distributed within the scope of ± 5mm/y, Ascendant trend is presented in PS point more than half, and the PS point in steel tower region is stablized in 0mm/y or so, has stability, wherein 3D phase Position solution twines algorithm resolving rate and more concentrates relatively.Elevation correction value is distributed within the scope of ± 10m, and Lambda algorithm is more than 80% collection In in ± 5m.In electric power pylon region, there is network adjustment smaller elevation to correct difference, more confidence level.Chengdu region The PS dot variable of rate of deformation 90% concentrates within the scope of ± 2.5mm/y, algorithms of different have certain proportion across 2.5~ 5mm/y has the point greater than 5mm/y deformation.Horizontal road region all show as it is unified it is stage rise or fall, elevation amendment Value 85% concentrates within the scope of ± 10m, and network adjustment method calculates the sparse elevation correction value of about 40m or so, this elevation Amendment may be built by commodity building in region to be formed.
Erlongshan Mountains elevation correction value is minimum, and it is related to be chronically at no developing zone with Erlongshan Mountains.There is artificial trace in Mao County area Mark and crustal movement combined influence, elevation also has to be changed to a certain degree.The elevation of In Chengdu has discrete sharp point, is built by city If work causes.
Table 4, table 5 and table 6 illustrate deformation of the earliest auxiliary image with respect to master image in three regions and are distributed.
4 Mao County deformation quantity of table counts mm
d≤-2.4 -2.4<d≤-1.2 -1.2<d≤0 0<d≤1.2 1.2<d≤2.4 d>2.4
Network adjustment 0.098 0.275 0.164 0.147 0.214 0.099
Pseudo- 3D 0.083 0.225 0.156 0.170 0.283 0.082
Lambda 0.004 0.319 0.220 0.180 0.273 0.002
5 Erlongshan Mountains deformation quantity of table counts mm
d≤-4.5 -4.5<d≤-3 -3<d≤0 0<d≤3 3<d≤4.5 d>4.5
Network adjustment 0 0 0.4845 0.515 0 0
Pseudo- 3D 0 0 0.545 0.455 0 0
Lambda 0.005 0.035 0.416 0.504 0.037 0.003
6 Chengdu deformation quantity of table counts cm
d≤-6 -6<d≤-3 -3<d≤0 0<d≤3 3<d≤6 d>6
Network adjustment 0 0.0865 0.4930 0.3301 0.0904 0
Pseudo- 3D 0 05230 0.5097 0.4553 0 0
Lambda 0 0.001 0.4368 0.5581 0.0040 0
In each research area, PS point in SAR image maximum imaging time interval have total deformation distribution, trend with Linear deformation rate is related to time interval.The deformation of maximum ± 30mm occurs for Mao County, and nearly 60% PS point deformation is in ± 12mm In deformation range.- 30~24mm deformation occurs for Erlongshan Mountains, and 70% deformation concentrates between ± 12mm.Chengdu occurs maximum The deformation of 120mm, most of PS dot changes are distributed within the scope of 30m.By Sichuan Electric Power Network artificial line walking record show Mao County and An iron rake with three to six teeth periphery is transmitted electricity in image imaging time as occurrence of large-area geological disaster in Erlongshan Mountains, does not also have collapse of iron tower case, says Bright deformation, which is extracted, has confidence level, and steel tower is not threatened by deformation safely.
The specific embodiment of invention described above is not intended to limit the scope of the present invention..

Claims (6)

1.一种基于PSInSAR技术的山区地形形变提取方法,其特征在于,包括:1. A mountain terrain deformation extraction method based on PSInSAR technology, is characterized in that, comprising: 选取SAR图像,所述SAR图像包括主图像和辅图像;Selecting a SAR image, the SAR image includes a main image and an auxiliary image; 利用所述主图像与所述辅图像配准并进行干涉处理,得到相干系数图和干涉相位图;Using the main image to register with the auxiliary image and performing interference processing to obtain a coherence coefficient map and an interferometric phase map; 利用振幅离差和相干系数的双重阈值法在所述相干系数图上选取PS点,得到PS点三角网;Utilize the double threshold method of amplitude deviation and coherence coefficient to select PS point on described coherence coefficient figure, obtain PS point triangulation; 将所述干涉相位图进行去平滤波处理,利用DEM数据对处理后的所述干涉相位图进行差分干涉处理,得到差分干涉图;performing flattening and filtering processing on the interferogram, and performing differential interference processing on the processed interferogram by using DEM data to obtain a differential interferogram; 对所述差分干涉图进行相位解缠处理,得到形变场。Phase unwrapping is performed on the differential interferogram to obtain a deformation field. 2.根据权利要求1所述的基于PSInSAR技术的山区地形形变提取方法,其特征在于,选取SAR图像具体包括:2. the mountain terrain deformation extraction method based on PSInSAR technology according to claim 1, is characterized in that, selecting SAR image specifically comprises: 选取n幅SAR图像,采用综合相干系数法选择所述SAR图像中的一幅作为所述主图像,其余n-1幅作为所述辅图像。Selecting n pieces of SAR images, using the comprehensive coherence method to select one of the SAR images as the main image, and the remaining n-1 pieces as the auxiliary images. 3.根据权利要求2所述的基于PSInSAR技术的山区地形形变提取方法,其特征在于,对所述差分干涉图进行相位解缠处理具体包括:3. the mountain terrain deformation extraction method based on PSInSAR technology according to claim 2, is characterized in that, carrying out phase unwrapping process to described differential interferogram specifically comprises: 分别利用网络平差法、伪3D相位解缠法和Lambda算法对所述差分干涉图进行相位解缠处理。The phase unwrapping process of the differential interferogram is carried out by network adjustment method, pseudo 3D phase unwrapping method and Lambda algorithm respectively. 4.根据权利要求3所述的基于PSInSAR技术的山区地形形变提取方法,其特征在于,所述网络平差法具体包括:4. the mountain terrain deformation extraction method based on PSInSAR technology according to claim 3, is characterized in that, described network adjustment method specifically comprises: 建立PS点的三次差分模型式,计算所述三次差分模型式中的相干系数K1,K2,B,所述三次差分模型式为:Establish the three-time difference model formula of the PS point, and calculate the coherence coefficients K 1 , K 2 , B in the three-time difference model formula, and the three-time difference model formula is: 引入铅垂线轨迹法,选取Δδ'H和Δv'的搜索范围以及步长,通过控制相邻PS点间的复相关系数γ逐步搜索合适的Δδ'H和Δv',其中,Introduce the plumb line trajectory method, select the search range and step size of Δδ' H and Δv', and gradually search for the appropriate Δδ' H and Δv' by controlling the complex correlation coefficient γ between adjacent PS points, where, 将所述相干系数的最大值所对应的线性形变速度差和高程修正差作为参数估计值,建立平差函数模型,所述平差函数模型为,Using the linear deformation velocity difference and the elevation correction difference corresponding to the maximum value of the coherence coefficient as parameter estimates, an adjustment function model is established, and the adjustment function model is, 求解线性方程:Solve a linear equation: 5.根据权利要求4所述的基于PSInSAR技术的山区地形形变提取方法,其特征在于,所述伪3D相位解缠法具体包括:5. the mountain terrain deformation extraction method based on PSInSAR technology according to claim 4, is characterized in that, described pseudo 3D phase unwrapping method specifically comprises: 对所述PS点做Delaunay三角网连接建立三角网,时间维通过时间关系连接;Do Delaunay triangulation connection to described PS point and set up triangulation, time dimension is connected by time relation; 沿所述三角网的边进行残差点计算;Carry out residual point calculation along the edge of the triangulation; 利用时间维解缠结果做空间2D相位解缠的最优化初始值,利用最小Lp范数法进行空间维最优求解。The optimal initial value of the space 2D phase unwrapping is made by using the unwrapping results of the time dimension, and the optimal solution of the space dimension is carried out by using the minimum Lp norm method. 6.根据权利要求5所述的基于PSInSAR技术的山区地形形变提取方法,其特征在于,所述Lambda算法具体包括:6. the mountain terrain deformation extraction method based on PSInSAR technology according to claim 5, is characterized in that, described Lambda algorithm specifically comprises: 在n-1幅差分干涉图中,将三角网相邻PS点的三次差分干涉相位式组成方程组:In n-1 differential interferograms, the cubic differential interferometric phase equations of adjacent PS points of the triangular network are composed of equations: 记为y1=A1·a+B1·b;Recorded as y 1 =A 1 ·a+B 1 ·b; 在y1=A1·a+B1·b中加入伪观测方程y2=A2·a+B2·b构成y=A·a+B·b,其中 A2为2*m零矩阵,B2为2*2单位矩阵;Add the pseudo-observation equation y 2 =A 2 ·a+B 2 ·b to y 1 =A 1 ·a+B 1 ·b to form y=A·a+B·b, where A 2 is a 2*m zero matrix, and B 2 is a 2*2 identity matrix; 根据最小二乘原理求得:According to the principle of least squares: 得到的整周浮动解,利用模糊度去相关方法搜索得到相位整周的固定解 get The floating solution of the whole cycle of , using the ambiguity decorrelation method to search for the fixed solution of the phase cycle 将所述固定解带入解缠相位表达式 中再次使用最小二乘法得到: the fixed solution into the unwrapped phase expression Using the method of least squares again to get: 完成相位解缠后,通过相高转换和地理编码获得所述形变场。After phase unwrapping, the deformation field is obtained by phase height transformation and geocoding.
CN201811123711.2A 2018-09-26 2018-09-26 Alpine terrain deformation extracting method based on PSInSAR technology Withdrawn CN109031301A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811123711.2A CN109031301A (en) 2018-09-26 2018-09-26 Alpine terrain deformation extracting method based on PSInSAR technology

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811123711.2A CN109031301A (en) 2018-09-26 2018-09-26 Alpine terrain deformation extracting method based on PSInSAR technology

Publications (1)

Publication Number Publication Date
CN109031301A true CN109031301A (en) 2018-12-18

Family

ID=64618086

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811123711.2A Withdrawn CN109031301A (en) 2018-09-26 2018-09-26 Alpine terrain deformation extracting method based on PSInSAR technology

Country Status (1)

Country Link
CN (1) CN109031301A (en)

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109489625A (en) * 2018-12-20 2019-03-19 国家测绘地理信息局第六地形测量队 A kind of urban area Ground Deformation monitoring method
CN110568440A (en) * 2019-09-10 2019-12-13 四川省地质工程勘察院集团有限公司 method for monitoring deformation of complex mountain area based on DS-InSAR technology
CN111059998A (en) * 2019-12-31 2020-04-24 中国地质大学(北京) A high-resolution time-series InSAR deformation monitoring method and system
CN111308468A (en) * 2019-11-27 2020-06-19 北京东方至远科技股份有限公司 Method for automatically identifying deformation risk area based on In SAR technology
CN111537999A (en) * 2020-03-04 2020-08-14 云南电网有限责任公司电力科学研究院 A Robust and Efficient Decomposed Projection Autofocus Method
CN111598929A (en) * 2020-04-26 2020-08-28 云南电网有限责任公司电力科学研究院 Two-dimensional unwrapping method based on time series differential interferometric synthetic aperture radar data
CN112213723A (en) * 2020-09-11 2021-01-12 速度时空信息科技股份有限公司 Method for monitoring landslide in real time by using SBAS technology
CN112268517A (en) * 2020-10-13 2021-01-26 内蒙古电力(集团)有限责任公司乌海超高压供电局 Method for monitoring deformation of power transmission tower equipment by PSInSAR
CN112363166A (en) * 2020-12-10 2021-02-12 长安大学 InSAR phase unwrapping method and system based on reliable redundant network
CN112698328A (en) * 2020-11-30 2021-04-23 四川大学 Phase unwrapping method and system for monitoring dam and landslide deformation GB-SAR
CN113191374A (en) * 2021-05-19 2021-07-30 甘肃省地震局(中国地震局兰州地震研究所) PolSAR image ridge line extraction method based on pyramid attention network
CN113311433A (en) * 2021-05-28 2021-08-27 北京航空航天大学 InSAR interferometric phase two-step unwrapping method combining quality map and minimum cost flow
CN113866765A (en) * 2021-09-24 2021-12-31 中国科学院精密测量科学与技术创新研究院 PS-InSAR measurement method based on multi-component time coherent model
CN114910907A (en) * 2022-05-07 2022-08-16 苏州深蓝空间遥感技术有限公司 River slope landslide risk monitoring method based on corner reflector and PS-InSAR technology

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102608601A (en) * 2012-03-08 2012-07-25 北京北科安地科技发展有限公司 Method for recognizing artificial corner reflector in SAR (synthetic aperture radar) image
CN105425236A (en) * 2015-11-24 2016-03-23 中清能绿洲科技股份有限公司 Automatic topographic surveying apparatus based on synthetic aperture radar imaging
CN106023157A (en) * 2016-05-10 2016-10-12 电子科技大学 Mountain area surface micro deformation information extraction method based on SAR images
CA3012551A1 (en) * 2016-02-03 2017-08-10 Deutsches Zentrum Fur Luft - Und Raumfahrt E.V. Method for calibrating an active sensor system
CN107132539A (en) * 2017-05-03 2017-09-05 中国地质科学院探矿工艺研究所 Landslide early-stage identification method of time sequence InSAR (interferometric synthetic Aperture Radar) based on small baseline set

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102608601A (en) * 2012-03-08 2012-07-25 北京北科安地科技发展有限公司 Method for recognizing artificial corner reflector in SAR (synthetic aperture radar) image
CN105425236A (en) * 2015-11-24 2016-03-23 中清能绿洲科技股份有限公司 Automatic topographic surveying apparatus based on synthetic aperture radar imaging
CA3012551A1 (en) * 2016-02-03 2017-08-10 Deutsches Zentrum Fur Luft - Und Raumfahrt E.V. Method for calibrating an active sensor system
CN106023157A (en) * 2016-05-10 2016-10-12 电子科技大学 Mountain area surface micro deformation information extraction method based on SAR images
CN107132539A (en) * 2017-05-03 2017-09-05 中国地质科学院探矿工艺研究所 Landslide early-stage identification method of time sequence InSAR (interferometric synthetic Aperture Radar) based on small baseline set

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
唐李: ""基于InSAR的输电通道地形形变提取研究"", 《中国优秀博硕士学位论文全文数据库(硕士) 工程科技Ⅱ辑》 *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109489625A (en) * 2018-12-20 2019-03-19 国家测绘地理信息局第六地形测量队 A kind of urban area Ground Deformation monitoring method
CN110568440A (en) * 2019-09-10 2019-12-13 四川省地质工程勘察院集团有限公司 method for monitoring deformation of complex mountain area based on DS-InSAR technology
CN111308468A (en) * 2019-11-27 2020-06-19 北京东方至远科技股份有限公司 Method for automatically identifying deformation risk area based on In SAR technology
CN111059998A (en) * 2019-12-31 2020-04-24 中国地质大学(北京) A high-resolution time-series InSAR deformation monitoring method and system
CN111537999A (en) * 2020-03-04 2020-08-14 云南电网有限责任公司电力科学研究院 A Robust and Efficient Decomposed Projection Autofocus Method
CN111537999B (en) * 2020-03-04 2023-06-30 云南电网有限责任公司电力科学研究院 A Robust and Efficient Autofocus Method for Decomposed Projection
CN111598929B (en) * 2020-04-26 2023-02-17 云南电网有限责任公司电力科学研究院 Two-Dimensional Unwrapping Method Based on Time-Series Differential Interferometric Synthetic Aperture Radar Data
CN111598929A (en) * 2020-04-26 2020-08-28 云南电网有限责任公司电力科学研究院 Two-dimensional unwrapping method based on time series differential interferometric synthetic aperture radar data
CN112213723A (en) * 2020-09-11 2021-01-12 速度时空信息科技股份有限公司 Method for monitoring landslide in real time by using SBAS technology
CN112268517A (en) * 2020-10-13 2021-01-26 内蒙古电力(集团)有限责任公司乌海超高压供电局 Method for monitoring deformation of power transmission tower equipment by PSInSAR
CN112698328A (en) * 2020-11-30 2021-04-23 四川大学 Phase unwrapping method and system for monitoring dam and landslide deformation GB-SAR
CN112363166A (en) * 2020-12-10 2021-02-12 长安大学 InSAR phase unwrapping method and system based on reliable redundant network
CN112363166B (en) * 2020-12-10 2023-08-22 长安大学 A Method and System for InSAR Phase Unwrapping Based on Reliable Redundant Network
CN113191374A (en) * 2021-05-19 2021-07-30 甘肃省地震局(中国地震局兰州地震研究所) PolSAR image ridge line extraction method based on pyramid attention network
CN113311433B (en) * 2021-05-28 2022-08-02 北京航空航天大学 InSAR interferometric phase two-step unwrapping method combining quality map and minimum cost flow
CN113311433A (en) * 2021-05-28 2021-08-27 北京航空航天大学 InSAR interferometric phase two-step unwrapping method combining quality map and minimum cost flow
CN113866765A (en) * 2021-09-24 2021-12-31 中国科学院精密测量科学与技术创新研究院 PS-InSAR measurement method based on multi-component time coherent model
CN113866765B (en) * 2021-09-24 2022-12-13 中国科学院精密测量科学与技术创新研究院 PS-InSAR measurement method based on multi-component time coherent model
CN114910907A (en) * 2022-05-07 2022-08-16 苏州深蓝空间遥感技术有限公司 River slope landslide risk monitoring method based on corner reflector and PS-InSAR technology

Similar Documents

Publication Publication Date Title
CN109031301A (en) Alpine terrain deformation extracting method based on PSInSAR technology
Dawei et al. Unmanned aerial vehicle (UAV) photogrammetry technology for dynamic mining subsidence monitoring and parameter inversion: A case study in China
CN101770027B (en) Surface three-dimensional deformation monitoring method based on InSAR and GPS data fusion
Catalão et al. Merging GPS and atmospherically corrected InSAR data to map 3-D terrain displacement velocity
CN106066478B (en) Merge the mining area surface deformation calculation method of pixel offset tracking and short baseline set
Neckel et al. Recent mass balance of the Purogangri Ice Cap, central Tibetan Plateau, by means of differential X-band SAR interferometry
CN110888130A (en) A method for monitoring surface deformation in coal mining area based on ascending and descending orbit time series InSAR
CN106772377A (en) A kind of building deformation monitoring method based on InSAR
CN105954747A (en) Tower foundation stability analyzing method based on three-dimensional deformation monitoring of unfavorable geologic body of power grid
CN107389029A (en) A kind of surface subsidence integrated monitor method based on the fusion of multi-source monitoring technology
CN106950556A (en) Heritage area deformation monitoring method based on distributed diffusion body sequential interference SAR technology
CN101551455A (en) 3D terrain imaging system of interferometric synthetic aperture radar and elevation mapping method thereof
CN110456346A (en) A kind of electric power pylon inclination monitoring method based on InSAR technology
Zeybek et al. Accurate determination of the Taşkent (Konya, Turkey) landslide using a long-range terrestrial laser scanner
CN114608527B (en) A fusion method of D-InSAR and UAV photogrammetry based on surface subsidence characteristics of coal mining
CN105824022A (en) Method for monitoring three-dimensional deformation of unfavorable geologic body under power grid
CN109100719A (en) Combine plotting method with the topographic map of optical image based on satellite-borne SAR image
CN113238228B (en) Method, system and device for acquiring 3D surface deformation based on horizontal constraints
CN107180448A (en) A kind of mining area DEM change detecting methods based on earth&#39;s surface invariant features
Guimarães et al. Radargrammetric approaches to the flat relief of the amazon coast using COSMO-SkyMed and TerraSAR-X datasets
Abdulrahman Determination of the local geoid model in Duhok Region, University of Duhok Campus as a Case study
Feng et al. A hierarchical network densification approach for reconstruction of historical ice velocity fields in East Antarctica
CN114674277B (en) Deep mining earth surface subsidence monitoring method by combining full-field survey lines
CN115236667B (en) A method for estimating potential landslide volume by fusing multi-source SAR data
Liu et al. Correction of positional errors and geometric distortions in topographic maps and DEMs using a rigorous SAR simulation technique

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
WW01 Invention patent application withdrawn after publication

Application publication date: 20181218

WW01 Invention patent application withdrawn after publication