CN110456345B - Building inclination monitoring method based on InSAR technology - Google Patents

Building inclination monitoring method based on InSAR technology Download PDF

Info

Publication number
CN110456345B
CN110456345B CN201910571244.8A CN201910571244A CN110456345B CN 110456345 B CN110456345 B CN 110456345B CN 201910571244 A CN201910571244 A CN 201910571244A CN 110456345 B CN110456345 B CN 110456345B
Authority
CN
China
Prior art keywords
building
inclination
monitoring
insar
monitoring points
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.)
Active
Application number
CN201910571244.8A
Other languages
Chinese (zh)
Other versions
CN110456345A (en
Inventor
平扬
王明洲
龚春龙
熊寻安
陈凯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Water Planning And Design Institute Co ltd
Original Assignee
Shenzhen Water Planning And Design Institute Co 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 Shenzhen Water Planning And Design Institute Co ltd filed Critical Shenzhen Water Planning And Design Institute Co ltd
Priority to CN201910571244.8A priority Critical patent/CN110456345B/en
Publication of CN110456345A publication Critical patent/CN110456345A/en
Application granted granted Critical
Publication of CN110456345B publication Critical patent/CN110456345B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021—SAR image post-processing techniques
    • G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C15/00—Surveying instruments or accessories not provided for in groups G01C1/00 - G01C13/00
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C9/00—Measuring inclination, e.g. by clinometers, by levels
    • 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

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

Abstract

The invention discloses a building inclination monitoring method based on an InSAR technology, relates to a synthetic aperture radar interferometry (InSAR) technology, and particularly relates to an InSAR inclination monitoring technology facing buildings and high-rise structures. The method comprises the following steps: the first step is as follows: acquiring a differential interferogram set by using a differential interferometry technique (DInSAR); the second step is that: acquiring the spatial position and deformation information of a building surface monitoring point set by using a time series InSAR technology; the third step: and calculating the inclination direction and inclination of the building and the high-rise structure by using the InSAR monitoring result and the satellite parameter information. The method comprises the steps of firstly, acquiring the spatial position and deformation information of a building surface monitoring point based on an InSAR technology, and acquiring the inclination direction and inclination of a building (structure) by combining satellite image parameters and orbit data. The method can provide monitoring results in a large-range wide area regularly, greatly improves the working efficiency, does not need to arrange monitoring equipment on site and manually operate, and greatly reduces the manual operation cost.

Description

Building inclination monitoring method based on InSAR technology
Technical Field
The invention relates to the technical field of Synthetic Aperture Radar (InSAR) interferometry, buildings, high-rise structures, tilt monitoring and the like, in particular to an InSAR tilt monitoring technology for buildings and high-rise structures.
Background
The construction (structure) may be inclined and deformed due to the influence of the deformation of the foundation or the structure during the construction operation and use. The main situations are as follows: the peripheral load of the building foundation is changed greatly, such as a large amount of piled soil; the foundation of the building per se is changed greatly, such as the foundation is soaked; the bearing structure of the building is changed or damaged due to the impact of strong external force; and the earthquake, landslide, flood and the like can be encountered. The building inclination deformation has great harm to the building and has direct influence on the service life of the building. Therefore, during the building construction operation and use period, the continuous and effective inclination monitoring is carried out on the building, the operation state of the building is reliably acquired, and the method has important significance for guaranteeing the safety of the building and the lives and properties of people.
For the inclination monitoring of a building (structure), common measurement methods include a total station observation method, a vertical alignment method, a relative settlement observation method and the like. The measurement method usually needs to set a measurement mark on site, and needs to consume a large amount of manpower and material resources for regular manual field operation. Due to the advantages of large-range, high-precision and all-weather measurement all day long, the InSAR technology is gradually maturely applied to the fields of urban ground settlement observation, landslide deformation monitoring and the like. However, this technique is generally used to acquire the deformation of the target in the radar sight line direction, and in special cases, the movement of the target in the vertical direction or the horizontal direction can be obtained through projective transformation, and is not used for building (structure) tilt monitoring.
The Synthetic Aperture Radar interferometry (Interferometric Synthetic Aperture Radar, inssar) refers to a technique for obtaining surface elevation and deformation information by using a satellite-borne or ground-based Radar image for interference processing. The method utilizes a radar to emit microwaves to a target area, then receives echoes reflected by the target to obtain SAR complex image pairs imaged by the same target area, and the SAR complex image pairs are subjected to conjugate multiplication to obtain an interferogram. And calculating the tiny change of the surface elevation of the target area according to the phase value of the interferogram, and using the tiny change for deformation monitoring.
Buildings include single-storey buildings, multi-storey buildings, high-rise buildings and super high-rise buildings. High-rise structures generally refer to structures with large height and relatively small cross section, such as signal towers, power transmission towers, chimneys, piers and the like.
The inclination monitoring generally refers to acquiring the inclination direction and inclination of a building by using observation methods such as a total station observation method, a vertical alignment method, a relative settlement observation method and a photogrammetry method.
Disclosure of Invention
In order to achieve the purpose of the invention, the invention provides a building inclination monitoring method based on InSAR technology.
The specific technical scheme of the invention is as follows: a building inclination monitoring method based on InSAR technology is characterized by comprising the following steps:
the first step is as follows: acquiring a differential interference pattern set by using a differential interference measurement technology;
the second step is that: acquiring the spatial position and deformation information of a building surface monitoring point set by using a time series InSAR technology;
the third step: and calculating the inclination direction and inclination of the building and the high-rise structure by using the InSAR monitoring result and the satellite parameter information.
Further, in the first step, a differential interferogram set is obtained through the steps of image registration, interference relative combination mode selection, interferogram generation, terrain phase simulation, terrain phase removal and the like;
the external DEM data used in the simulated terrain phase step may comprise low resolution SRTM data or high resolution terrain data generated based on oblique photography techniques.
Furthermore, in the second step, the spatial position and deformation information of the building surface monitoring point set are obtained through the steps of coherent target point selection, phase unwrapping, parameter estimation, orbit error phase estimation, atmospheric error phase estimation, seasonal temperature model phase estimation, deformation sequence acquisition, geocoding and the like. The seasonal temperature model phase estimation aims to eliminate the phase caused by the expansion and contraction effect of the building.
Further, the seasonal temperature model phase estimation method is based on a relation between the phase difference of the monitoring points and the thermal expansion factor of the building, and the relation between the seasonal temperature model phase and the thermal expansion factor of the building is expressed as follows:
Figure BDA0002110953580000031
in the formula (I), the compound is shown in the specification,
Figure BDA0002110953580000032
the phase difference value of monitoring points i and j in the kth unwrapped differential interferogram is the phase from which the elevation error and the linear deformation rate are removed; λ is the wavelength of the electromagnetic wave emitted by the SAR satellite; delta TempkRepresenting the atmospheric temperature difference of the shooting time of the two SAR images corresponding to the kth scene difference interference pattern; Δ Hi,jThe elevation difference value of the monitoring points i and j is obtained; tep is the thermal expansion coefficient of the monitored target building, which is to be estimated;
Figure BDA0002110953580000033
and obtaining phase residuals corresponding to monitoring points i and j in the kth interferogram.
Further, in the third step, the following steps are specifically applied to calculate the inclination direction and the inclination of the building and the towering structure:
selecting any two monitoring points on the three-dimensional surface of the building, wherein the positions of the monitoring points are represented by X, Y, H, X is an east coordinate value, Y is a north coordinate value, H represents the elevation of the monitoring points, the inclination direction of the building is set to be theta, the inclination is set to be q, and delta H is set to be delta Hi,jIs the elevation difference value, delta D, of two monitoring points i and j of the vertical face of the buildingi,jFor the distance between the projections of two monitoring points i and j of the building facade on the horizontal plane, the deformation difference value between any two monitoring points of the building solid surface can be expressed as follows:
ΔdXi,j＝ΔHi,j*cosθ*q
ΔdYi,j＝ΔHi,j*sinθ*q
ΔdHi,j＝ΔDi,j*cos(θ-αi,j)*q
in the formula,. DELTA.dXi,j、ΔdYi,j、ΔdHi,jRespectively representing the deformation difference values alpha of the two monitoring points i and j in the east direction, the north direction and the elevation directioni,jThe included angle between the connecting line of the projection points of the two monitoring points i and j on the horizontal plane and the east direction can be specifically expressed as follows:
Figure BDA0002110953580000041
Figure BDA0002110953580000042
in the formula,. DELTA.Xi,j、ΔYi,jThe coordinate difference of the two monitoring points i and j in the east direction and the north direction is obtained;
the InSAR result is the projection of the deformation of the monitoring points in the radar sight line direction, the image parameters and the orbit data of the radar satellite are considered, and the projection of the deformation difference value of the two monitoring points i and j in the radar sight line direction (LOS)
Figure BDA0002110953580000043
Expressed as:
Figure BDA0002110953580000044
in the formula, betai,jThe incidence angle of the radar satellite image is usually between 20 and 50 degrees, and the radar satellite image incidence angle and the radar satellite orbit inclination angle parameter can be obtained by utilizing an SAR image parameter file.
Further, InSAR data processing is carried out by utilizing the high-resolution radar satellite image, a plurality of effective monitoring points can be obtained on the vertical surface of the building, at most N x (N-1)/2 equations are formed for the N monitoring points according to a relational expression of the deformation of the InSAR monitoring points and the inclination value of the building, and the inclination direction and the inclination of the building can be obtained by adopting parameter estimation methods such as a least square method, a maximum likelihood method and the like.
Further, in the third step, a functional relation between the difference value of the deformation of any two InSAR monitoring points on the three-dimensional surface of the building in the radar sight line direction and the inclination of the building is formed;
the InSAR result is the projection of the deformation of the monitoring points in the radar sight line direction, the image parameters and the orbit data of the radar satellite are considered, and the projection of the deformation difference value of the two monitoring points i and j in the radar sight line direction (LOS)
Figure BDA0002110953580000051
Expressed as:
Figure BDA0002110953580000052
in the formula, betai,jThe incidence angle of the radar satellite image is usually between 20 and 50 degrees, and the radar satellite image incidence angle and the radar satellite orbit inclination angle parameter can be obtained by utilizing an SAR image parameter file.
Further, in the third step, calculating the inclination direction and inclination of the building according to a relational expression of the deformation of the InSAR monitoring points and the inclination value of the building;
the method comprises the steps of utilizing a high-resolution radar satellite image to process InSAR data, obtaining a plurality of effective monitoring points on a building facade, forming at most N x (N-1)/2 equations for N monitoring points according to a relational expression of deformation of the InSAR monitoring points and a building inclination value, and obtaining the inclination direction and the inclination of the building by adopting parameter estimation methods such as a least square method, a maximum likelihood method and the like.
The invention has the beneficial effects that: the invention provides an InSAR inclination monitoring method facing buildings and high-rise structures. The method can provide monitoring results in a large-range wide area regularly, greatly improves the working efficiency, does not need to arrange monitoring equipment on site and manually operate, and greatly reduces the manual operation cost.
Drawings
FIG. 1 is a flow chart of the working principle of the first embodiment of the present invention;
FIG. 2a is a schematic structural diagram of a building InSAR inclination monitoring test;
FIG. 2b is a schematic view of the exterior profile structure of a high-rise building;
fig. 3 is a schematic diagram of the accumulated deformation of the InSAR tilt monitoring test of the building.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
On the contrary, the invention is intended to cover alternatives, modifications, equivalents and alternatives which may be included within the spirit and scope of the invention as defined by the appended claims. Furthermore, in the following detailed description of the present invention, certain specific details are set forth in order to provide a better understanding of the present invention. It will be apparent to one skilled in the art that the present invention may be practiced without these specific details.
As shown in fig. 1, which is a flow chart of the operation principle of the first embodiment of the present invention, the method for monitoring InSAR tilt for buildings and high-rise structures according to the present invention can realize continuous and effective periodic monitoring of buildings (structures). The key implementation strategy of the method mainly comprises differential interference processing, InSAR time sequence analysis, and inclination direction and inclination estimation, and the specific implementation steps are shown in FIG. 1.
(1) Differential interference processing
The method mainly comprises the steps of carrying out interference processing on an SAR image sequence according to an interference pattern combination mode to obtain a differential interference pattern set, wherein the differential interference pattern set comprises image registration, interference relative combination mode selection, interference pattern generation, terrain phase simulation, terrain phase removal and the like.
And A1, SAR image registration, namely registering the SAR image sequence with the main image, wherein the registering comprises coarse registration and fine registration. The rough registration mainly utilizes satellite orbit parameters or manually selects a small number of feature points, and calculates the offset of the image to be registered (auxiliary image) relative to the main image in the azimuth direction (row direction) and the distance direction (column direction).
A2 interference relative combination mode selection: common combination modes include a single main image combination mode, a small baseline set combination mode, and a free combination mode. The basic idea of the single-master image combination mode is to select one image as a common master image, and all the other images are slave images. The small baseline set combination mode is to select a short baseline interference pair by using a spatial baseline threshold value method so as to weaken the image of spatial incoherent. The free combination mode is that all images in the SAR image set can be subjected to interference pairing in an arbitrary combination mode.
A3 interferogram generation: and according to the interference relative combination relation, carrying out complex conjugate multiplication on the main and auxiliary images after the fine registration to generate a corresponding interference image. Distance-wise spectral filtering may be performed relatively for interference with a longer spatial baseline to eliminate distance-wise spectral shifts.
A4 simulated terrain phase: and establishing a mapping relation between the image coordinate and the geodetic coordinate according to the SAR imaging parameters by using external DEM data to obtain elevation data of an SAR image coordinate system, and calculating the simulated terrain phase of each pixel by using the image orbit data and the elevation data of the SAR image. The external DEM data used includes, but is not limited to, low resolution SRTM data or high resolution topographic data generated based on oblique photography techniques.
A5 go to terrain phase: and removing the simulated terrain phase from the interference pattern to obtain a differential interference pattern. The differential interference phase comprises a surface deformation phase, an elevation error phase, an atmospheric error phase, a track error phase, a noise phase and the like.
(2) InSAR timing analysis
The method mainly comprises the steps of obtaining deformation sequences and elevation information of effective monitoring points on the three-dimensional surface of a building (structure) through time sequence analysis, and comprises the contents of coherent target point selection, phase unwrapping, parameter estimation, orbit error phase estimation, atmospheric error phase estimation, seasonal temperature model phase estimation, deformation sequence obtaining, geocoding and the like.
B1 coherent target point selection step: the coherent target point is selected as a target point which utilizes the intensity or phase information of the SAR image to identify the surface of the iron tower with higher signal-to-noise ratio. Common coherent point selection methods include amplitude dispersion, coherence threshold, intensity threshold, signal-to-noise threshold, and the like. Taking the coherence threshold method as an example, firstly, a coherence coefficient is calculated based on the differential interference pattern, and a pixel point with the coherence coefficient larger than a set threshold is a candidate coherence target point.
B2 phase unwrapping: the phase unwrapping mainly solves the problem of 2 pi ambiguity in an interference pattern, and restores the absolute phase value of a pixel point, and can be divided into space phase unwrapping and three-dimensional phase unwrapping. Common spatial phase unwrapping methods include the twitter method, the least cost stream method, and the least squares method. Three-dimensional phase unwrapping is a method that combines both time and space dimensional phase unwrapping.
B3 parameter estimation: the estimated parameters include elevation values and deformation rates of the monitoring points, etc. Firstly, establishing a functional relation among phase difference values, elevation difference values and deformation rate difference values among InSAR monitoring points:
Figure BDA0002110953580000091
in the formula (I), the compound is shown in the specification,
Figure BDA0002110953580000092
and λ is the wavelength of the electromagnetic wave transmitted by the SAR satellite, and is the phase difference value of monitoring points i and j in the unwrapped differential interferogram of the kth scene. Delta TkRepresenting the difference value of the shooting time of two SAR images adopted by the kth scene difference interferogram, delta Vi,jAs a monitoring pointThe difference in the deformation rates of i and j,
Figure BDA0002110953580000093
is the vertical baseline length, R, of the kth view differential interferogrami,j、βi,jRespectively, the slant range and the incidence angle of the radar satellite, Δ Hi,jFor the elevation difference of monitoring points i and j,
Figure BDA0002110953580000094
and obtaining phase residuals corresponding to monitoring points i and j in the kth interferogram.
The differential interference image set generated by the SAR image sequence is utilized, and based on a functional relation of a phase difference value, an elevation difference value and a deformation rate difference value among InSAR monitoring points, estimation of the elevation difference value and the deformation rate difference value among the building surface monitoring points can be obtained by adopting parameter estimation methods such as a least square method, a maximum likelihood method and the like. 1 monitoring point with known ground elevation is selected in a stable area near a building, and then the elevation and the deformation rate of the building three-dimensional surface monitoring point can be obtained.
B4 orbit error phase estimation: the track error phase is the phase error caused by the inaccuracy of the baseline estimate. Common orbit error correction methods include baseline accurate estimation and interferometric phase error correction. The baseline accurate estimation is a method for refining the InSAR interference baseline by using a ground control point. The interference phase error correction mostly adopts a mode of fitting a polynomial model. Finally, the estimated track error phase is removed from the differential interferogram.
B5 atmospheric error phase estimation: the atmospheric delay correction method mainly comprises an empirical method and a prediction method. The empirical method mainly utilizes the time-independent characteristic of the atmospheric delay phase, realizes the estimation of the atmospheric delay phase through time-dimensional high-pass filtering and spatial low-pass filtering, and removes the atmospheric delay phase from an interferogram. The prediction method is to directly calculate the atmospheric delay amount when the SAR image is imaged by using external meteorological data (temperature, air pressure, humidity or water vapor content information). Finally, the estimated atmospheric error phase is removed from the differential interferogram.
B6 seasonal temperature model phase estimation step: the seasonal temperature model phase removal is mainly to remove deformation phases generated by the influence of seasonal thermal expansion and cold contraction effects on the building.
The seasonal temperature model phase estimation method is based on a relational expression of a phase difference of monitoring points and a building thermal expansion factor, the building thermal expansion factor is estimated by adopting a least square method and a maximum likelihood method, so that the seasonal temperature model phase is obtained, and finally the estimated seasonal temperature model phase is removed from a differential interference map. The seasonal temperature model phase versus building thermal expansion factor may be expressed as:
Figure BDA0002110953580000101
in the formula (I), the compound is shown in the specification,
Figure BDA0002110953580000102
the phase difference values of monitoring points i and j in the k th unwrapped differential interferogram are the phases with the elevation error and the linear deformation rate removed. λ is the wavelength of the electromagnetic wave emitted by the SAR satellite. Delta TempkAnd representing the atmospheric temperature difference of the shooting time of the two SAR images corresponding to the kth scene difference interference pattern. Δ Hi,jIs the elevation difference between monitoring points i and j. tep is the thermal expansion coefficient of the monitored target building, which is to be evaluated.
Figure BDA0002110953580000103
And obtaining phase residuals corresponding to monitoring points i and j in the kth interferogram.
The InSAR data processing is carried out by utilizing the high-resolution radar satellite image, a plurality of effective monitoring points can be obtained on the vertical face of the building, based on the relational expression, the thermal expansion coefficient of the building can be obtained by adopting parameter estimation methods such as a least square method, a maximum likelihood method and the like, and the seasonal temperature model phase estimation is realized.
B7 deformation sequence acquisition step: the elevation error phase, the orbit error phase, the atmospheric error phase and the seasonal temperature model phase are removed from the differential interference phase, and then the phase sequence is converted into a deformation sequence based on satellite parameter information.
B8 geocoding step: and geocoding is to establish a mapping relation between SAR image coordinates and geodetic coordinates, and geocoding is carried out on monitoring point deformation monitoring results so that the monitoring points can be superposed and analyzed with other basic geographic information. In general, the geodetic coordinates of the monitoring points are acquired as longitude and latitude information under the ellipsoids of the CGCS2000 or the WGS 84. In consideration of the convenience of subsequent analysis, the invention can carry out projection transformation on the longitude and latitude coordinates of the monitoring points.
(3) Tilt direction and tilt estimation
The main content of the part is to estimate the inclination direction and the inclination of a building (structure) by using the resolving parameters and the deformation information of InSAR monitoring points. The specific estimation procedure is as follows:
selecting any two monitoring points on the three-dimensional surface of the building, in the embodiment, calculating the inclination direction and the inclination of the building by using all the monitoring points with different heights on the three-dimensional surface of the building, wherein the positions of the monitoring points are represented by X, Y, H, X is an east coordinate value, Y is a north coordinate value, H is an elevation of the monitoring points, the inclination direction of the building is set to be theta, the inclination is set to be q, and delta H is set to be delta Hi,jIs the elevation difference value, delta D, of two monitoring points i and j of the vertical face of the buildingi,jFor the distance between the projections of two monitoring points i and j of the building facade on the horizontal plane, the deformation difference value between any two monitoring points of the building solid surface can be expressed as follows:
ΔdXi,j＝ΔHi,j*cosθ*q
ΔdYi,j＝ΔHi,j*sinθ*q
ΔdHi,j＝ΔDi,j*cos(θ-αi,j)*q
in the formula,. DELTA.dXi,j、ΔdYi,j、ΔdHi,jRespectively representing the deformation difference values alpha of the two monitoring points i and j in the east direction, the north direction and the elevation directioni,jThe included angle between the connecting line of the projection points of the two monitoring points i and j on the horizontal plane and the east direction can be specifically expressed as follows:
Figure BDA0002110953580000121
Figure BDA0002110953580000122
in the formula,. DELTA.Xi,j、ΔYi,jThe coordinate difference of the two monitoring points i and j in the east direction and the north direction is obtained;
the InSAR result is the projection of the deformation of the monitoring points in the radar sight line direction, the image parameters and the orbit data of the radar satellite are considered, and the projection of the deformation difference value of the two monitoring points i and j in the radar sight line direction (LOS)
Figure BDA0002110953580000123
Expressed as:
Figure BDA0002110953580000124
in the formula, betai,jThe incident angle of the radar satellite image is usually between 20 ° and 50 °. For the inclination angle of the orbit of the radar satellite, in this embodiment, taking the radar satellite (terrasaar-X, COSMO-SkyMed, Sentinel-1, etc.) currently in service as an example, and taking the local latitude of shenzhen as a reference, the inclination angle of the orbit of the ascending radar satellite is about-11 ° 00 ', and the inclination angle of the orbit of the descending radar satellite is about 11 ° 40'. The radar satellite image incidence angle and the radar satellite orbit inclination angle parameters can be obtained by utilizing an SAR image parameter file.
The method comprises the steps of utilizing a high-resolution radar satellite image to process InSAR data, obtaining a plurality of effective monitoring points on a building facade, forming at most N x (N-1)/2 equations for N monitoring points according to a relational expression of deformation of the InSAR monitoring points and a building inclination value, and obtaining the inclination direction and the inclination of the building by adopting parameter estimation methods such as a least square method, a maximum likelihood method and the like.
Further, in this embodiment, as shown in fig. 2a and 2b, an elevated rail TerraSAR-X image with a resolution of 3 meters is selected to perform the building InSAR tilt monitoring test. In the experiment, 44 scenes of SAR images are collected together, the starting time is 11 months and 27 days in 2013, the cut-off time is 2016 years and 10 months and 29 days, the image operation waveband is an X waveband, and the single scene image coverage range is 30km multiplied by 30 km. Fig. 2a shows the relative geometry of the elevated SAR satellite to the building in this example, where the angle of incidence phi of the SAR satellite is 37.3 deg.. Fig. 2b is an average intensity graph of the SAR image, and it can be seen that the reflected signal of the high-rise building is strong and the external contour is clear.
In the test, a common main image interference combination mode is adopted to carry out interference processing on the SAR image sequence, 43 differential interferograms are obtained in total, and then an InSAR time sequence analysis technology is utilized to obtain a deformation sequence and elevation information of the effective monitoring points on the three-dimensional surface of the building.
As shown in fig. 3, the position information of the effective monitoring point of the three-dimensional surface of the building in the test and the accumulated deformation amount during the monitoring period. In the figure, the X axis points to the north direction, the Y axis points to the east direction, the H axis is the elevation, and the color depth of the points represents the magnitude of the accumulated deformation of the monitoring points in the radar sight line direction. Analysis shows that the InSAR technology can obtain more effective monitoring points on the surface of the building, the accumulated deformation of the monitoring points at different positions has larger difference, the monitoring points at the bottom area of the building are shown to sink in the radar sight direction, the maximum deformation is-13 mm, the monitoring points at the top area of the building are shown to lift in the radar sight direction, and the maximum deformation is 38 mm.
Based on the relation between the deformation of InSAR monitoring points and the inclination value (the inclination direction and the inclination rate) of the building, the inclination direction and the inclination rate of the building in the case are obtained by a project group by adopting a least square method. The building has an inclination direction of 36 degrees in the southwest and an inclination of 0.08 percent (the ratio of the settlement difference of the two monitoring points in the inclination direction to the distance).
The method for monitoring the InSAR inclination facing the building and the high-rise structure is applied, firstly, the method acquires the inclination direction and the inclination of the building (structure) based on the space position and the deformation information of the monitoring point on the surface of the building, which are acquired by the InSAR technology, and combines the satellite image parameters and the orbit data. The method can provide monitoring results in a large-range wide area regularly, greatly improves the working efficiency, does not need to arrange monitoring equipment on site and manually operate, and greatly reduces the manual operation cost.

Claims (7)

1. A building inclination monitoring method based on InSAR technology is characterized by comprising the following steps:
the first step is as follows: acquiring a differential interference pattern set by using a differential interference measurement technology;
the second step is that: acquiring the spatial position and deformation information of a building surface monitoring point set by using a time series InSAR technology;
the third step: calculating the inclination direction and inclination of the building and the high-rise structure by using the InSAR monitoring result and the satellite parameter information;
in the third step, the following steps are specifically applied to calculate the inclination direction and the inclination of the building and the high-rise structure:
selecting any two monitoring points on the three-dimensional surface of the building, wherein the positions of the monitoring points are represented by X, Y, H, X is an east coordinate value, Y is a north coordinate value, H represents the elevation of the monitoring points, the inclination direction of the building is set to be theta, the inclination is set to be q, and delta H is set to be delta Hi,jIs the elevation difference value, delta D, of two monitoring points i and j of the vertical face of the buildingi,jFor the distance between the projections of two monitoring points i and j of the building facade on the horizontal plane, the deformation difference value between any two monitoring points of the building solid surface can be expressed as follows:
ΔdXi,j＝ΔHi,j*cosθ*q
ΔdYi,j＝ΔHi,j*sinθ*q
ΔdHi,j＝ΔDi,j*cos(θ-αi,j)*q
in the formula,. DELTA.dXi,j、ΔdYi,j、ΔdHi,jRespectively representing the deformation difference values alpha of the two monitoring points i and j in the east direction, the north direction and the elevation directioni,jThe included angle between the connecting line of the projection points of the two monitoring points i and j on the horizontal plane and the east direction can be specifically expressed as follows:
Figure FDA0002617184360000021
Figure FDA0002617184360000022
in the formula,. DELTA.Xi,j、ΔYi,jThe coordinate difference of the two monitoring points i and j in the east direction and the north direction is obtained;
the InSAR result is the projection of the deformation of the monitoring points in the radar sight line direction, the image parameters and the orbit data of the radar satellite are considered, and the projection of the deformation difference value of the two monitoring points i and j in the radar sight line direction (LOS)
Figure FDA0002617184360000023
Expressed as:
Figure FDA0002617184360000024
in the formula, betai,jThe incidence angle of the radar satellite image is usually between 20 and 50 degrees, and the radar satellite image incidence angle and the radar satellite orbit inclination angle parameter can be obtained by utilizing an SAR image parameter file.
2. The method for monitoring the inclination of the building based on the InSAR technology as claimed in claim 1, characterized in that in the first step, a differential interferogram set is obtained through the steps of image registration, selection of an interference relative combination mode, generation of interferograms, terrain phase simulation, terrain phase removal and the like;
the external DEM data used in the simulated terrain phase step may comprise low resolution SRTM data or high resolution terrain data generated based on oblique photography techniques.
3. The method for monitoring the inclination of the building based on the InSAR technology as claimed in claim 1, wherein in the second step, the spatial position and the deformation information of the monitoring point set on the surface of the building are obtained through the steps of coherent target point selection, phase unwrapping, parameter estimation, orbit error phase estimation, atmospheric error phase estimation, seasonal temperature model phase estimation, deformation sequence acquisition, geocoding and the like, wherein the seasonal temperature model phase estimation aims at eliminating the phase caused by the expansion and contraction effect of the building.
4. The InSAR technology-based building tilt monitoring method according to claim 3, wherein the seasonal temperature model phase estimation method is based on a relationship between a monitoring point phase difference and a building thermal expansion factor, and the relationship between the seasonal temperature model phase and the building thermal expansion factor is expressed as:
Figure FDA0002617184360000031
in the formula (I), the compound is shown in the specification,
Figure FDA0002617184360000032
the phase difference value of monitoring points i and j in the kth unwrapped differential interferogram is the phase from which the elevation error and the linear deformation rate are removed; λ is the wavelength of the electromagnetic wave emitted by the SAR satellite; delta TempkRepresenting the atmospheric temperature difference of the shooting time of the two SAR images corresponding to the kth scene difference interference pattern; Δ Hi,jThe elevation difference value of the monitoring points i and j is obtained; tep is the thermal expansion coefficient of the monitored target building, which is to be estimated;
Figure FDA0002617184360000033
and obtaining phase residuals corresponding to monitoring points i and j in the kth interferogram.
5. The method for monitoring the inclination of the building based on the InSAR technology as claimed in claim 1, wherein the InSAR data processing is performed by using a high-resolution radar satellite image, a plurality of effective monitoring points can be obtained on the facade of the building, at most N (N-1)/2 equations are formed for N monitoring points according to the relational expression of the deformation of the InSAR monitoring points and the inclination value of the building, and the inclination direction and the inclination of the building can be obtained by adopting a parameter estimation method such as a least square method, a maximum likelihood method and the like.
6. The method for monitoring the inclination of the building based on the InSAR technology as claimed in claim 1, wherein in the third step, the functional relation between the difference of the deformation of any two InSAR monitoring points on the stereoscopic surface of the building in the direction of the radar sight line and the inclination direction and inclination of the building is expressed;
the InSAR result is the projection of the deformation of the monitoring points in the radar sight line direction, the image parameters and the orbit data of the radar satellite are considered, and the projection of the deformation difference value of the two monitoring points i and j in the radar sight line direction (LOS)
Figure FDA0002617184360000041
Expressed as:
Figure FDA0002617184360000042
in the formula, betai,jThe incidence angle of the radar satellite image is usually between 20 and 50 degrees, and the radar satellite image incidence angle and the radar satellite orbit inclination angle parameter can be obtained by utilizing an SAR image parameter file.
7. The method for monitoring the inclination of the building based on the InSAR technology as claimed in claim 1, wherein in the third step, the method for calculating the inclination direction and inclination of the building is based on the relational expression between the deformation of the InSAR monitoring points and the inclination value of the building;
the method comprises the steps of utilizing a high-resolution radar satellite image to process InSAR data, obtaining a plurality of effective monitoring points on a building facade, forming at most N x (N-1)/2 equations for N monitoring points according to a relational expression of deformation of the InSAR monitoring points and a building inclination value, and obtaining the inclination direction and the inclination of the building by adopting parameter estimation methods such as a least square method, a maximum likelihood method and the like.
CN201910571244.8A 2019-06-28 2019-06-28 Building inclination monitoring method based on InSAR technology Active CN110456345B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910571244.8A CN110456345B (en) 2019-06-28 2019-06-28 Building inclination monitoring method based on InSAR technology

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910571244.8A CN110456345B (en) 2019-06-28 2019-06-28 Building inclination monitoring method based on InSAR technology

Publications (2)

Publication Number Publication Date
CN110456345A CN110456345A (en) 2019-11-15
CN110456345B true CN110456345B (en) 2020-11-10

Family

ID=68481796

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910571244.8A Active CN110456345B (en) 2019-06-28 2019-06-28 Building inclination monitoring method based on InSAR technology

Country Status (1)

Country Link
CN (1) CN110456345B (en)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111323768B (en) * 2020-02-26 2023-12-12 北京佳格天地科技有限公司 Building change identification method, device, terminal and storage medium
CN111398959B (en) * 2020-04-07 2023-07-04 中南大学 InSAR time sequence earth surface deformation monitoring method based on earth surface stress strain model
CN111610519A (en) * 2020-06-11 2020-09-01 石家庄铁道大学 Non-contact type measuring method for dynamic deformation of small bridge structure
CN112284332B (en) * 2020-08-31 2021-10-08 北京四象爱数科技有限公司 High-rise building settlement monitoring result three-dimensional positioning method based on high-resolution INSAR
CN112268517B (en) * 2020-10-13 2023-04-28 内蒙古电力(集团)有限责任公司乌海超高压供电局 Method for monitoring deformation of power transmission tower equipment by using PSInSAR
CN112485790B (en) * 2020-11-23 2023-11-24 中大智能科技股份有限公司 Track non-contact deformation high-precision measurement method based on K-band radar
CN113589226B (en) * 2021-02-05 2023-09-22 北京深蓝长盛科技有限公司 Method and system for calculating body inclination of athlete
CN114491770A (en) * 2022-02-21 2022-05-13 江苏建研建设工程质量安全鉴定有限公司 House inclination detection and calculation method based on three-dimensional laser scanning point cloud
CN115790360A (en) * 2023-02-03 2023-03-14 中大智能科技股份有限公司 Three-dimensional deformation measurement method
CN116973877B (en) * 2023-09-22 2023-12-12 南京楚航科技有限公司 Millimeter wave radar deformation measurement method, system and measurement truth value calibration method

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102927934A (en) * 2012-11-07 2013-02-13 中南大学 Method for obtaining mining area earth surface three-dimensional deformation fields through single interferometric synthetic aperture radar (InSAR) interference pair
CN103675790A (en) * 2013-12-23 2014-03-26 中国国土资源航空物探遥感中心 Method for improving earth surface shape change monitoring precision of InSAR (Interferometric Synthetic Aperture Radar) technology based on high-precision DEM (Digital Elevation Model)
CN105954747A (en) * 2016-06-20 2016-09-21 中国电力工程顾问集团中南电力设计院有限公司 Tower foundation stability analyzing method based on three-dimensional deformation monitoring of unfavorable geologic body of power grid
CN106772377A (en) * 2017-01-18 2017-05-31 深圳市路桥建设集团有限公司 A kind of building deformation monitoring method based on InSAR
WO2018027332A1 (en) * 2016-08-08 2018-02-15 Comercial E Industrial Gesecology Limitada Method and system for the analysis and generation of early or predictive alerts concerning the stability of slopes in open-pit mines

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102927934A (en) * 2012-11-07 2013-02-13 中南大学 Method for obtaining mining area earth surface three-dimensional deformation fields through single interferometric synthetic aperture radar (InSAR) interference pair
CN103675790A (en) * 2013-12-23 2014-03-26 中国国土资源航空物探遥感中心 Method for improving earth surface shape change monitoring precision of InSAR (Interferometric Synthetic Aperture Radar) technology based on high-precision DEM (Digital Elevation Model)
CN105954747A (en) * 2016-06-20 2016-09-21 中国电力工程顾问集团中南电力设计院有限公司 Tower foundation stability analyzing method based on three-dimensional deformation monitoring of unfavorable geologic body of power grid
WO2018027332A1 (en) * 2016-08-08 2018-02-15 Comercial E Industrial Gesecology Limitada Method and system for the analysis and generation of early or predictive alerts concerning the stability of slopes in open-pit mines
CN106772377A (en) * 2017-01-18 2017-05-31 深圳市路桥建设集团有限公司 A kind of building deformation monitoring method based on InSAR

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Revised Improved DINSAR Algorithm for Monitoring the Inclination Displacement of Top Position of Electric Power Transmission Tower;Baolong Wu等;《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》;20180630;第15卷(第6期);第877-881页 *
The Thermal Expansion Component of Persistent Scatterer Interferometry Observations;Oriol Monserrat等;《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》;20110930;第8卷(第5期);第865页 *
建筑物倾斜与沉降监测结果综合分析;花向红等;《地理空间信息》;20090228;第7卷(第1期);第120-122页 *
监测城市基础设施健康的星载ＭＴ－ＩｎＳＡＲ方法介绍;林珲等;《测绘学报》;20171031;第46卷(第10期);第1421-1433页 *

Also Published As

Publication number Publication date
CN110456345A (en) 2019-11-15

Similar Documents

Publication Publication Date Title
CN110456345B (en) Building inclination monitoring method based on InSAR technology
CN110456346B (en) Power transmission tower inclination monitoring method based on InSAR technology
Liu et al. Detecting subsidence in coastal areas by ultrashort-baseline TCPInSAR on the time series of high-resolution TerraSAR-X images
CN109031301A (en) Alpine terrain deformation extracting method based on PSInSAR technology
CN110017762B (en) Method for monitoring deformation of creep body by using Offset-Tracking technology
CN108663678B (en) Multi-baseline InSAR phase unwrapping algorithm based on mixed integer optimization model
CN112051571A (en) LOS (line of sight) deformation variable estimation method of novel differential InSAR (interferometric synthetic Aperture Radar)
CN112711021B (en) Multi-resolution InSAR (interferometric synthetic Aperture Radar) interactive interference time sequence analysis method
CN107918127A (en) A kind of road slope deformation detecting system and method based on vehicle-mounted InSAR
CN112444188B (en) Multi-view InSAR sea wall high-precision three-dimensional deformation measurement method
Nikolakopoulos et al. Updating the 1: 50.000 topographic maps using ASTER and SRTM DEM: the case of Athens, Greece
Feng et al. An improved geometric calibration model for spaceborne SAR systems with a case study of large-scale Gaofen-3 images
Yu et al. Digital Elevation Model generation using ascending and Descending multi-baseline ALOS/PALSAR radar images
CN115343709B (en) Terrain inversion method suitable for distributed spaceborne SAR system
Selvakumaran et al. Waterloo bridge monitoring: Comparing measurements from earth and space
CN113219414B (en) Novel method for eliminating satellite interference radar surface deformation direction blurring
CN114488144A (en) SAR offset three-dimensional deformation estimation method and system with additional DEM constraint
Nonaka et al. The relationships between errors of DEM and the height of ambiguity of Sentinel-1
Kelevitz et al. Novel corner-reflector array application in essential infrastructure monitoring
CN114706074A (en) Surface three-dimensional deformation InSAR measurement method 2T3D-InSAR based on two-dimensional observation
Rao et al. Comparison of DEMs derived from INSAR and optical stereo techniques
Zhang et al. Repeat-pass SAR interferometry experiments with Gaofen-3: a case study of Ningbo and Nanjing area
Du et al. Insar Reveals the Long Term Subsidence and Potential Landdegradation in Mexico City from 2004 to 2018 with Five Sar Sensors
Sarago et al. Operational use of SAR interferometry for surface and infrastructures movement monitoring
CN117572378B (en) Mountain settlement analysis method and device based on InSAR and Beidou data

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