CN113884054B - Tower displacement settlement monitoring method and system - Google Patents
Tower displacement settlement monitoring method and system Download PDFInfo
- Publication number
- CN113884054B CN113884054B CN202111150066.5A CN202111150066A CN113884054B CN 113884054 B CN113884054 B CN 113884054B CN 202111150066 A CN202111150066 A CN 202111150066A CN 113884054 B CN113884054 B CN 113884054B
- Authority
- CN
- China
- Prior art keywords
- target point
- image
- phase
- displacement
- tower
- 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
Links
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 87
- 238000012544 monitoring process Methods 0.000 title claims abstract description 61
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000004062 sedimentation Methods 0.000 claims abstract description 49
- 238000012545 processing Methods 0.000 claims abstract description 35
- 238000010183 spectrum analysis Methods 0.000 claims abstract description 27
- 230000005540 biological transmission Effects 0.000 claims abstract description 23
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 18
- 238000007781 pre-processing Methods 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 62
- 239000006185 dispersion Substances 0.000 claims description 15
- 238000001228 spectrum Methods 0.000 claims description 8
- 230000008030 elimination Effects 0.000 claims description 7
- 238000003379 elimination reaction Methods 0.000 claims description 7
- 230000003595 spectral effect Effects 0.000 claims description 7
- 238000012935 Averaging Methods 0.000 claims 2
- 238000005516 engineering process Methods 0.000 description 16
- 238000003860 storage Methods 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 5
- 208000028257 Joubert syndrome with oculorenal defect Diseases 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 238000012876 topography Methods 0.000 description 3
- 238000004590 computer program Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- DMBHHRLKUKUOEG-UHFFFAOYSA-N diphenylamine Chemical compound C=1C=CC=CC=1NC1=CC=CC=C1 DMBHHRLKUKUOEG-UHFFFAOYSA-N 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000802 evaporation-induced self-assembly Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000005305 interferometry Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C5/00—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' 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
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The application provides a method and a system for monitoring displacement and settlement of a tower, wherein the method comprises the following steps: acquiring an SAR image, and preprocessing the SAR image to obtain a registration image set; wherein the registered image set comprises a plurality of registered images; performing spectrum analysis on each registration image, and selecting a target point according to a spectrum analysis result; performing noise processing on the target point through differential interference processing and SVD decomposition to obtain ground deformation information; acquiring the current deformation and the historical deformation of the transmission tower; according to the ground deformation information and the historical deformation quantity, a tower displacement sedimentation monitoring model is established, and a displacement sedimentation judgment threshold value is obtained; and obtaining a tower displacement sedimentation judgment result through the current deformation quantity and the displacement sedimentation judgment threshold value. The application can obtain deformation information with high precision and high time resolution.
Description
Technical Field
The application relates to the technical field of tower displacement sedimentation monitoring, in particular to a tower displacement sedimentation monitoring method and system.
Background
The pole tower is basic equipment in overhead distribution lines, if geological settlement occurs around the pole tower, the power transmission lines can be damaged when serious, the pole tower collapses, the local area is subjected to power failure, and great inconvenience is brought to production and life, so that measures are required to be taken to detect and early warn geology and the pole tower nearby the pole tower, and whether geological disaster potential hazards exist nearby the pole tower or not is judged according to actual conditions.
The more common monitoring methods at present are InSAR, GPS and optical remote sensing technologies. These conventional geological measurement techniques often suffer from the disadvantages of large workload, long measurement time, low degree of automation, etc. Some measuring means can only measure the relative displacement of a geological deformation area, and the measuring range is limited for a site with a larger area, for example, the traditional InSAR technology is easy to be influenced by related factors such as atmospheric delay, space-time interference and the like, the time resolution is low, and the timeliness of monitoring the deformation of geological disasters is not strong; only the deformation during two-scene SAR images, rather than dynamic subsidence sequences, is typically acquired. In addition, in the monitoring and early warning problem of geological disasters, especially geological settlement, the remote sensing image is influenced by factors such as weather and environment, so that InSAR can be combined with the high-time resolution characteristic of the Beidou positioning navigation system in practical application, and high-precision topographic data are utilized to calculate and predict deformation degree, sweep range and other data, but the quantity of variables required to be referred for prediction is relatively large, and the accuracy of the acquired topographic data is still different from the actual demand.
Disclosure of Invention
In order to solve the problems in the prior art, the application provides a tower displacement settlement monitoring method and system, which can monitor the surface deformation by combining Beidou positioning and SBAS-InSAR technology, thereby reducing the influence of atmospheric and space-time interference, obtaining deformation information with high precision and high time resolution, and comparing the deformation information with a displacement settlement judgment threshold value to obtain a transmission line displacement settlement judgment result.
The first aspect of the application provides a method for monitoring displacement and settlement of a tower, which comprises the following steps:
acquiring an SAR image, and preprocessing the SAR image to obtain a registration image set; wherein the registered image set comprises a plurality of registered images;
performing spectrum analysis on each registration image, and selecting a target point according to a spectrum analysis result;
performing noise processing on the target point through differential interference processing and SVD decomposition to obtain ground deformation information;
acquiring the current deformation and the historical deformation of the transmission tower;
according to the ground deformation information and the historical deformation quantity, a tower displacement sedimentation monitoring model is established, and a displacement sedimentation judgment threshold value is obtained;
and obtaining a tower displacement sedimentation judgment result through the current deformation quantity and the displacement sedimentation judgment threshold value.
Further, preprocessing the SAR image to obtain a registered image set, specifically:
selecting one SAR image as a main image, and registering the SAR images at different moments with the main image to obtain registered images corresponding to the SAR images at different moments;
all the registered images are formed into a registered image set.
Further, the spectrum analysis is performed on each of the registered images, and a target point is selected according to the spectrum analysis result, specifically:
selecting a target point by performing spectrum analysis on each registered image, wherein the target point satisfies the following conditions:
wherein N is the total number of registered images, alpha i For the spectral correlation value of pixel i in the image, beta i Is the dispersion of pixel i, α A and βA The average value of the spectrum correlation value and the dispersion in the SAR image is respectively.
Further, the performing noise processing on the target point through differential interference processing and SVD decomposition to obtain ground deformation information includes:
performing error elimination and topographic phase interference removal on the differential interference pattern of the target point, and converting the phase of the differential interference pattern of the target point into a matrix to obtain a coefficient matrix;
and processing the coefficient matrix through SVD decomposition to eliminate the influence of interference and errors.
Further, the phase of the differential interferogram of the target point is expressed by the following formula:
wherein ,Δφi (x, y) is the phase of the ith differential interferogram of the target point with coordinates (x, y), phi (t) z X, y) is the interference phase of the target point with coordinates (x, y) in the main image, phi (t) f X, y) is the interference phase of the target point with coordinates (x, y) in the auxiliary image, d (t) z X, y) is the displacement of the target point with coordinates (x, y) in the main image in the radar direction, d (t) f X, y) is the displacement of the target point with coordinates of (x, y) in the auxiliary image in the radar direction, and lambda is the radar wavelength;
the phase of the differential interferogram of the target point is converted into a matrix, and the matrix is expressed by the following formula:
wherein ,is the phase matrix of the differential interferogram, delta is the coefficient of the phase matrix of the differential interferogram, +.>Is a phase matrix, and H is a coefficient matrix.
A second aspect of the present application provides a tower displacement settlement monitoring system comprising:
the registration module is used for acquiring SAR images and preprocessing the SAR images to obtain a registration image set; wherein the registered image set comprises a plurality of registered images;
the target point selection module is used for carrying out spectrum analysis on each registration image and selecting a target point according to a spectrum analysis result;
the noise processing module is used for carrying out noise processing on the target point through differential interference processing and SVD decomposition to obtain ground deformation information;
the transmission tower deformation acquisition module is used for acquiring the current deformation and the historical deformation of the transmission tower;
the tower displacement sedimentation monitoring model building module is used for building a tower displacement sedimentation monitoring model according to the ground deformation information and the historical deformation quantity to obtain a displacement sedimentation judgment threshold;
and the tower displacement sedimentation judging module is used for obtaining a tower displacement sedimentation judging result through the current deformation and the displacement sedimentation judging threshold value.
Further, the registration module is further configured to:
selecting one SAR image as a main image, and registering the SAR images at different moments with the main image to obtain registered images corresponding to the SAR images at different moments;
all the registered images are formed into a registered image set.
Further, the target point selection module is further configured to:
selecting a target point by performing spectrum analysis on each registered image, wherein the target point satisfies the following conditions:
wherein N is the total number of registered images, alpha i For the spectral correlation value of pixel i in the image, beta i Is the dispersion of pixel i, α A and βA The average value of the spectrum correlation value and the dispersion in the SAR image is respectively.
Further, the noise processing module is further configured to:
performing error elimination and topographic phase interference removal on the differential interference pattern of the target point, and converting the phase of the differential interference pattern of the target point into a matrix to obtain a coefficient matrix;
and processing the coefficient matrix through SVD decomposition to eliminate the influence of interference and errors.
Further, the phase of the differential interferogram of the target point is expressed by the following formula:
wherein ,Δφi (x, y) is the phase of the ith differential interferogram of the target point with coordinates (x, y), phi (t) z X, y) is the interference phase of the target point with coordinates (x, y) in the main image, phi (t) f X, y) is the interference phase of the target point with coordinates (x, y) in the auxiliary image, d (t) z X, y) is the displacement of the target point with coordinates (x, y) in the main image in the radar direction, d (t) f X, y) is the displacement of the target point with coordinates of (x, y) in the auxiliary image in the radar direction, and lambda is the radar wavelength;
the phase of the differential interferogram of the target point is converted into a matrix, and the matrix is expressed by the following formula:
wherein ,phase matrix for differential interferogramsDelta is the coefficient of the phase matrix of the differential interferogram,/->Is a phase matrix, and H is a coefficient matrix.
Compared with the prior art, the embodiment of the application has the beneficial effects that:
the application provides a method and a system for monitoring displacement and settlement of a tower, wherein the method comprises the following steps: acquiring an SAR image, and preprocessing the SAR image to obtain a registration image set; wherein the registered image set comprises a plurality of registered images; performing spectrum analysis on each registration image, and selecting a target point according to a spectrum analysis result; performing noise processing on the target point through differential interference processing and SVD decomposition to obtain ground deformation information; acquiring the current deformation and the historical deformation of the transmission tower; according to the ground deformation information and the historical deformation quantity, a tower displacement sedimentation monitoring model is established, and a displacement sedimentation judgment threshold value is obtained; and obtaining a tower displacement sedimentation judgment result through the current deformation quantity and the displacement sedimentation judgment threshold value. The application combines Beidou positioning and SBAS-InSAR technology, and can monitor the surface deformation, thereby reducing the influence of atmosphere and space-time interference. In addition, because the Beidou positioning system has the advantages of quick positioning, full-automatic monitoring and full-automatic monitoring, deformation quantity of the designated direction coordinates of the tower is monitored by combining Beidou positioning, deformation information with high precision and high time resolution can be obtained, and a displacement sedimentation judgment result of the power transmission line is obtained by comparing the deformation information with a displacement sedimentation judgment threshold. The method can improve the limitation of the deformation monitoring technology under the support of high-time-efficiency Beidou positioning and high-precision data, and plays a great role in monitoring, evaluating and early warning of geological disasters and pole and tower settlement.
Drawings
In order to more clearly illustrate the technical solutions of the present application, the drawings that are needed in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present application, and that other drawings can be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method for monitoring displacement and settlement of a tower according to an embodiment of the present application;
FIG. 2 is a flow chart of a small baseline set method provided by an embodiment of the application;
FIG. 3 is a general flow chart of a Beidou positioning and SBAS-InSAR bi-directional fusion technique provided by an embodiment of the present disclosure;
FIG. 4 is a schematic diagram of a tower displacement settlement monitoring system according to an embodiment of the present application;
fig. 5 is a block diagram of an electronic device according to an embodiment of the present application.
Detailed Description
The following description of the embodiments of the present application will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present application, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application.
It should be understood that the step numbers used herein are for convenience of description only and are not limiting as to the order in which the steps are performed.
It is to be understood that the terminology used in the description of the application is for the purpose of describing particular embodiments only and is not intended to be limiting of the application. As used in this specification and the appended claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise.
The terms "comprises" and "comprising" indicate the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.
The term "and/or" refers to any and all possible combinations of one or more of the associated listed items, and includes such combinations.
With the rapid development of information technology, it is now possible to monitor and evaluate geological settlement disasters by using various technologies. Synthetic aperture radar interferometry (InSAR) is a radar technology with long development time, is commonly used in the fields of remote sensing and mapping and the field of ground subsidence monitoring and analysis, has good sensitivity to regional micro deformation, and can improve the observation precision to a certain extent and enhance the integrity and predictability of a study object compared with the traditional technology. In addition, the Beidou positioning technology can better meet the early warning requirement of geological disasters, has good practicability in the aspect of geological settlement monitoring due to the advantages of quick positioning, high measurement accuracy, uninterrupted throughout the day and the like, and can play an important role in geological disaster early warning work.
Liu Li et al (Liu Li, li Changan, gao Junhua, yu Deqing, xia Le) discuss [ J ]. Geological science and technology information based on geological disaster monitoring key problems of Beidou and InSAR, 2019,38 (06): 141-149.) propose a time sequence InSAR technology based on CR constraint, and a high-time resolution earth surface deformation information of a geological disaster area is obtained by combining a multi-platform lifting rail InSAR cooperative monitoring geological disaster method with the high-time resolution characteristic of Beidou, so that a relatively perfect geological disaster monitoring technology method is formed. Feng Zhijiang et al (Feng Zhijiang, hu Danhui, wang Tao, zhou Xueming, ma Jianguo, zhang Yaodong, zhou Ding, lu Ping) discloses a transmission line displacement sedimentation determination method and system based on Beidou positioning [ P ] Beijing city: CN110749307A,2020-02-04. The transmission line displacement sedimentation determination method and system based on Beidou positioning are disclosed, coordinates and deformation of a transmission tower are calculated by utilizing a differential positioning algorithm, then a displacement sedimentation determination threshold value is determined by using an established prediction model, and compared with the deformation of a currently set direction coordinate to realize the determination of the displacement sedimentation of the transmission tower, so that a more reliable prediction method is required to provide deformation information data for geological disaster prediction, and more reasonable modeling, parameter inversion and monitoring and early warning are sought by analyzing the change rule of the deformation information data.
In a first aspect.
Referring to fig. 1, an embodiment of the present application provides a method for monitoring displacement and settlement of a tower, including:
s10, acquiring SAR images, and preprocessing the SAR images to obtain a registration image set. Wherein the registered image set comprises a plurality of registered images.
S20, carrying out spectrum analysis on each registration image, and selecting a target point according to a spectrum analysis result.
S30, performing noise processing on the target point through differential interference processing and SVD decomposition to obtain ground deformation information.
S40, acquiring the current deformation quantity and the historical deformation quantity of the transmission tower.
And S50, establishing a tower displacement sedimentation monitoring model according to the ground deformation information and the historical deformation quantity to obtain a displacement sedimentation judgment threshold.
And S60, obtaining a tower displacement sedimentation judgment result through the current deformation quantity and the displacement sedimentation judgment threshold value.
In a specific implementation of the embodiment of the present application, the step S10 is specifically:
selecting one SAR image as a main image, and registering the SAR images at different moments with the main image to obtain registered images corresponding to the SAR images at different moments;
all the registered images are formed into a registered image set.
In a specific implementation of the embodiment of the present application, the step S20 specifically includes:
selecting a target point by performing spectrum analysis on each registered image, wherein the target point satisfies the following conditions:
wherein N is the total number of registered images, alpha i For the spectral correlation value of pixel i in the image, beta i Is the dispersion of pixel i, α A and βA The average value of the spectrum correlation value and the dispersion in the SAR image is respectively.
In a specific implementation manner of the embodiment of the present application, the step S30 includes:
performing error elimination and topographic phase interference removal on the differential interference pattern of the target point, and converting the phase of the differential interference pattern of the target point into a matrix to obtain a coefficient matrix;
and processing the coefficient matrix through SVD decomposition to eliminate the influence of interference and errors.
Preferably, the phase of the differential interferogram of the target point is expressed by the following formula:
wherein ,Δφi (x, y) is the phase of the ith differential interferogram of the target point with coordinates (x, y), phi (t) z X, y) is the interference phase of the target point with coordinates (x, y) in the main image, phi (t) f X, y) is the interference phase of the target point with coordinates (x, y) in the auxiliary image, d (t) z X, y) is the displacement of the target point with coordinates (x, y) in the main image in the radar direction, d (t) f X, y) is the displacement of the target point with coordinates of (x, y) in the auxiliary image in the radar direction, and lambda is the radar wavelength;
the phase of the differential interferogram of the target point is converted into a matrix, and the matrix is expressed by the following formula:
wherein ,is the phase matrix of the differential interferogram, delta is the coefficient of the phase matrix of the differential interferogram, +.>Is a phase matrix, and H is a coefficient matrix.
In a specific implementation manner of the embodiment of the present application, the step S40 includes:
and establishing an observation equation through Beidou positioning, and calculating to obtain coordinates of the transmission tower.
And calculating the current deformation and the historical deformation according to the coordinates of the transmission tower.
The method provided by the application combines Beidou positioning and SBAS-InSAR technology, and can monitor the surface deformation, thereby reducing the influence of atmospheric and space-time interference. In addition, because the Beidou positioning system has the advantages of quick positioning, full-automatic monitoring and full-automatic monitoring, deformation quantity of the designated direction coordinates of the tower is monitored by combining Beidou positioning, deformation information with high precision and high time resolution can be obtained, and a displacement sedimentation judgment result of the power transmission line is obtained by comparing the deformation information with a displacement sedimentation judgment threshold. The method can improve the limitation of the deformation monitoring technology under the support of high-time-efficiency Beidou positioning and high-precision data, and plays a great role in monitoring, evaluating and early warning of geological disasters and pole and tower settlement.
In another embodiment of the application, the disaster monitoring and early warning for the geological region near the tower mainly comprises the following parts: data processing, differential interference calculation, time/space domain deformation estimation and deformation calculation, beidou positioning combination, disaster early warning, damage evaluation and other research works. After monitoring and early warning, grading the severity of the damage of the tower by combining the real-time information and the historical information, and setting operation and maintenance measures according to the grades to complete the design of the whole tower displacement settlement monitoring and early warning method.
SBAS-InSAR monitoring surface deformation:
the application adopts a small baseline set method to carry out the inversion of the surface deformation, and the main steps are as shown in figure 2:
(1) Small baseline set registration: firstly, selecting an SAR as a main image, and registering images at different moments with the main image. The space-time baseline distribution is distributed to different interference sets to generate a differential interference pattern and a plurality of interference pairs, so that the interference pairs in the same set have small space baseline distance, and the interference pairs in different sets have large space baseline distance.
(2) Selecting PS points: the application adopts a frequency spectrum characteristic method: the PS points are selected by carrying out spectrum analysis on the registered SAR images, and the conditions to be satisfied are shown as follows:
wherein N is the total number of SAR images, alpha i For the spectral correlation value of pixel i in the image, beta i Is the dispersion of pixel i, α A and βA The average value of the spectrum correlation value and the dispersion in the SAR image is respectively.
After the PS point is extracted, operations such as phase cancellation are continued.
(3) Differential interference processing: the historical topography data, the topography digital elevation model and the like are synthesized, error elimination is carried out on the differential interference diagram, the topography phase interference is removed, and the phase of the ith differential interference diagram of the PS point with coordinates of (x, y) is expressed as follows:
wherein ,φ(tz X, y) and phi (t f X, y) are the interference phases of the PS points of the primary and secondary images, d (t) z X, y) and d (t f X, y) are displacements in the radar direction, respectively, and λ is the radar wavelength.
Convert it into a matrix form and decompose it in the next step:
wherein ,is a differential phase matrix, delta is the coefficient of the phase matrix of the differential interference pattern, and H is the coefficient momentArray (S)>Is a phase matrix.
(4) H is processed by SVD decomposition, and the influence of interference and errors is eliminated. Wherein, let H be the coefficient matrix of the differential phase matrix of SAR image, carry on SVD decomposition to it:
H＝USV T
the main information of the image will then be converted into a series of singular values sigma i S is a diagonal matrix, U and V are sigma i Orthogonal matrix of u i 、v i Each corresponding to U, V. By using the least square inversion idea, the expressions of the phases at different moments can be obtained as follows:
the column vectors of the deformation phases at different moments in the step (4) are respectively filtered to finally separate out atmospheric error phases, the error phases are removed to obtain the needed time change and deformation quantity information, and the required time change and deformation quantity information is combined with the tower displacement information of the next stage to jointly determine the tower displacement settlement quantity.
2. Shaft tower subsides and judges based on big dipper location:
deformation monitoring to ground has been realized in the last stage, and in this stage, subside the judgement to the shaft tower through big dipper location. The improved Beidou CORS settlement monitoring method is mainly adopted, the carrier phase is optimized in data processing by using a satellite positioning principle continuously carried out by CORS, and data such as the carrier phase of a reference station and the carrier phase of a monitoring station are processed simultaneously by a double-difference resolving mode, so that more accurate settlement of the tower is obtained, and the aim of monitoring the settlement of the tower is fulfilled.
(1) Firstly, the coordinates of a transmission tower at the current moment need to be obtained, an observation equation is established through Beidou positioning, beidou CORS geological settlement monitoring is adopted, the method mainly comprises real-time ambiguity resolution and error correction, and a parameter estimation equation is as follows:
in the above formula, lambda is the phase wavelength, N is the ambiguity,d is the geometric distance from the satellite to the receiver, S is the multipath effect error, and gamma is the receiver noise; i. j is a satellite label; and r and u are the marks of the main reference station and the mobile station respectively, and the coordinates of the transmission tower can be obtained through calculation by a differential positioning algorithm.
(2) And (3) setting a fixed point for the tower, calculating the current deformation quantity and the historical deformation quantity of the tower according to the tower coordinates, and performing time sequencing on the historical deformation quantity to generate deformation quantities with a plurality of fixed time intervals, wherein the removal of gross errors is considered in the calculation of the historical deformation quantity, and the detail is shown in (3).
(3) The method for removing the gross error is as follows: and combining the auxiliary judgment of the geological deformation in the previous stage to obtain the geological deformation at the corresponding moment, and arranging the same discrete time series observed quantity according to the sequence from small to large in time.
In terms of error reduction, mean Absolute Error (MAE) metrics are utilized to measure the concentration and degree of dispersion of the data,as the deformation data average value, eliminating the data with overlarge errors:
(4) and (3) establishing a prediction model to predict the deformation of the coordinates in the appointed direction of the fixed point set at the current moment as a displacement sedimentation judgment threshold, and comparing the calculated current deformation of the coordinates in the appointed direction of the fixed point set with the displacement sedimentation judgment threshold to obtain a tower displacement sedimentation judgment result. The general scheme flow is shown in figure 3.
The advantages of this approach are three:
(1) The application fully plays the advantages of the SBAS-InSAR technology and the Beidou positioning method, can realize high-precision and high-time-resolution real-time monitoring, effectively early warn the geological settlement range caused by geological disasters and the disaster receiving degree of the towers, and provide rapid and effective data support for the early warning of the disaster receiving degree of the geological disasters and the post-disaster operation and maintenance of the towers.
(2) The application adopts an SBAS-InSAR small baseline set registration method. Complex terrain deformation information is decomposed by SVD, so that the influence of interference and errors is reduced, and the calculation accuracy is improved.
(3) According to the application, the Beidou positioning technology is combined, the current deformation quantity and the historical deformation quantity of the tower are calculated by using an improved Beidou CORS settlement monitoring method, and the rough differences in the variables are removed, so that real-time tower settlement monitoring information with higher precision and higher time resolution can be obtained.
The second aspect.
Referring to fig. 4, an embodiment of the present application provides a tower displacement settlement monitoring system, which includes:
the registration module 10 is configured to acquire a SAR image, and perform preprocessing on the SAR image to obtain a registered image set. Wherein the registered image set comprises a plurality of registered images.
The target point selection module 20 is configured to perform spectrum analysis on each of the registered images, and select a target point according to a spectrum analysis result.
The noise processing module 30 is configured to perform noise processing on the target point through differential interference processing and SVD decomposition, so as to obtain ground deformation information.
And the transmission tower deformation acquisition module 40 is used for acquiring the current deformation amount and the historical deformation amount of the transmission tower.
The tower displacement sedimentation monitoring model building module 50 is configured to build a tower displacement sedimentation monitoring model according to the ground deformation information and the historical deformation amount, so as to obtain a displacement sedimentation determination threshold.
And the tower displacement sedimentation judging module 60 is configured to obtain a tower displacement sedimentation judging result according to the current deformation amount and the displacement sedimentation judging threshold value.
In a specific implementation of the embodiment of the present application, the registration module 10 is further configured to:
selecting one SAR image as a main image, and registering the SAR images at different moments with the main image to obtain registered images corresponding to the SAR images at different moments;
all the registered images are formed into a registered image set.
In a specific implementation manner of the embodiment of the present application, the target point selection module 20 is further configured to:
selecting a target point by performing spectrum analysis on each registered image, wherein the target point satisfies the following conditions:
wherein N is the total number of registered images, alpha i For the spectral phase of pixel i in the imageOff value, beta i Is the dispersion of pixel i, α A and βA The average value of the spectrum correlation value and the dispersion in the SAR image is respectively.
In a specific implementation of the embodiment of the present application, the noise processing module 30 is further configured to:
performing error elimination and topographic phase interference removal on the differential interference pattern of the target point, and converting the phase of the differential interference pattern of the target point into a matrix to obtain a coefficient matrix;
and processing the coefficient matrix through SVD decomposition to eliminate the influence of interference and errors.
Preferably, the phase of the differential interferogram of the target point is expressed by the following formula:
wherein ,Δφi (x, y) is the phase of the ith differential interferogram of the target point with coordinates (x, y), phi (t) z X, y) is the interference phase of the target point with coordinates (x, y) in the main image, phi (t) f X, y) is the interference phase of the target point with coordinates (x, y) in the auxiliary image, d (t) z X, y) is the displacement of the target point with coordinates (x, y) in the main image in the radar direction, d (t) f X, y) is the displacement of the target point with coordinates of (x, y) in the auxiliary image in the radar direction, and lambda is the radar wavelength;
the phase of the differential interferogram of the target point is converted into a matrix, and the matrix is expressed by the following formula:
wherein ,is the phase matrix of the differential interferogram, delta is the coefficient of the phase matrix of the differential interferogram, +.>Is a phase matrix, and H is a coefficient matrix.
In a specific implementation manner of the embodiment of the present application, the transmission tower deformation obtaining module 40 is further configured to:
and establishing an observation equation through Beidou positioning, and calculating to obtain coordinates of the transmission tower.
And calculating the current deformation and the historical deformation according to the coordinates of the transmission tower.
The application combines Beidou positioning and SBAS-InSAR technology, and can monitor the surface deformation, thereby reducing the influence of atmosphere and space-time interference. In addition, because the Beidou positioning system has the advantages of quick positioning, full-automatic monitoring and full-automatic monitoring, deformation quantity of the designated direction coordinates of the tower is monitored by combining Beidou positioning, deformation information with high precision and high time resolution can be obtained, and a displacement sedimentation judgment result of the power transmission line is obtained by comparing the deformation information with a displacement sedimentation judgment threshold. The method can improve the limitation of the deformation monitoring technology under the support of high-time-efficiency Beidou positioning and high-precision data, and plays a great role in monitoring, evaluating and early warning of geological disasters and pole and tower settlement.
In a third aspect.
The present application provides an electronic device including:
a processor, a memory, and a bus;
the bus is used for connecting the processor and the memory;
the memory is used for storing operation instructions;
the processor is configured to, by invoking the operation instruction, cause the processor to execute an operation corresponding to a tower displacement settlement monitoring method according to the first aspect of the present application.
In an alternative embodiment, an electronic device is provided, as shown in fig. 5, the electronic device 5000 shown in fig. 5 includes: a processor 5001 and a memory 5003. The processor 5001 is coupled to the memory 5003, e.g., via bus 5002. Optionally, the electronic device 5000 may also include a transceiver 5004. It should be noted that, in practical applications, the transceiver 5004 is not limited to one, and the structure of the electronic device 5000 is not limited to the embodiment of the present application.
The processor 5001 may be a CPU, general purpose processor, DSP, ASIC, FPGA or other programmable logic device, transistor logic device, hardware component, or any combination thereof. Which may implement or perform the various exemplary logic blocks, modules and circuits described in connection with this disclosure. The processor 5001 may also be a combination of computing functions, e.g., including one or more microprocessor combinations, a combination of a DSP and a microprocessor, etc.
Bus 5002 may include a path to transfer information between the aforementioned components. Bus 5002 may be a PCI bus or an EISA bus, among others. The bus 5002 may be divided into an address bus, a data bus, a control bus, and the like. For ease of illustration, only one thick line is shown in fig. 5, but not only one bus or one type of bus.
The memory 5003 may be, but is not limited to, ROM or other type of static storage device, RAM or other type of dynamic storage device, which can store static information and instructions, EEPROM, CD-ROM or other optical disk storage, optical disk storage (including compact disk, laser disk, optical disk, digital versatile disk, blu-ray disc, etc.), magnetic disk storage or other magnetic storage devices, or any other medium capable of carrying or storing desired program code in the form of instructions or data structures and capable of being accessed by a computer.
The memory 5003 is used for storing application program codes for implementing the inventive arrangements and is controlled to be executed by the processor 5001. The processor 5001 is operative to execute application code stored in the memory 5003 to implement what has been shown in any of the method embodiments described previously.
Among them, electronic devices include, but are not limited to: mobile terminals such as mobile phones, notebook computers, digital broadcast receivers, PDAs (personal digital assistants), PADs (tablet computers), PMPs (portable multimedia players), in-vehicle terminals (e.g., in-vehicle navigation terminals), and the like, and stationary terminals such as digital TVs, desktop computers, and the like.
A fourth aspect.
The present application provides a computer readable storage medium having a computer program stored thereon, which when executed by a processor implements a method for monitoring displacement and settlement of a shaft tower according to the first aspect of the present application.
Yet another embodiment of the present application provides a computer-readable storage medium having stored thereon a computer program which, when run on a computer, causes the computer to perform the corresponding ones of the foregoing method embodiments.
Claims (6)
1. A method for monitoring displacement and settlement of a tower, comprising the steps of:
acquiring an SAR image, and preprocessing the SAR image to obtain a registration image set; wherein the registered image set comprises a plurality of registered images;
performing spectrum analysis on each registration image, and selecting a target point according to a spectrum analysis result;
performing noise processing on the target point through differential interference processing and SVD decomposition to obtain ground deformation information;
acquiring the current deformation and the historical deformation of the transmission tower;
according to the ground deformation information and the historical deformation quantity, a tower displacement sedimentation monitoring model is established, and a displacement sedimentation judgment threshold value is obtained;
obtaining a tower displacement sedimentation judgment result through the current deformation quantity and the displacement sedimentation judgment threshold value;
the spectrum analysis is performed on each registration image, and a target point is selected according to a spectrum analysis result, specifically:
selecting a target point by performing spectrum analysis on each registered image, wherein the target point satisfies the following conditions:
wherein N is the total number of registered images, alpha i For the spectral correlation value of pixel i in the image, beta i Is the dispersion of pixel i, α A and βA Respectively averaging the spectrum correlation value and the dispersion in the SAR image;
the noise processing is carried out on the target point through differential interference processing and SVD decomposition to obtain ground deformation information, and the method comprises the following steps:
performing error elimination and topographic phase interference removal on the differential interference pattern of the target point, and converting the phase of the differential interference pattern of the target point into a matrix to obtain a coefficient matrix;
and processing the coefficient matrix through SVD decomposition to eliminate the influence of interference and errors.
2. The method for monitoring displacement and settlement of towers according to claim 1, wherein the preprocessing is performed on the SAR image to obtain a registered image set, specifically:
selecting one SAR image as a main image, and registering the SAR images at different moments with the main image to obtain registered images corresponding to the SAR images at different moments;
all the registered images are formed into a registered image set.
3. A tower displacement settlement monitoring method as claimed in claim 1, wherein the phase of the differential interferogram of the target point is expressed by the following formula:
wherein ,Δφi (x, y) is the phase of the ith differential interferogram of the target point with coordinates (x, y), phi (t) z X, y) is the interference phase of the target point with coordinates (x, y) in the main image, phi (t) f X, y) is the interference phase of the target point with coordinates (x, y) in the auxiliary image, d (t) z X, y) is the displacement of the target point with coordinates (x, y) in the main image in the radar direction, d (t) f X, y) is the displacement of the target point with coordinates of (x, y) in the auxiliary image in the radar direction, and lambda is the radar wavelength;
the phase of the differential interferogram of the target point is converted into a matrix, and the matrix is expressed by the following formula:
wherein ,is the phase matrix of the differential interferogram, delta is the coefficient of the phase matrix of the differential interferogram, +.>Is a phase matrix, and H is a coefficient matrix.
4. A tower displacement settlement monitoring system, comprising:
the registration module is used for acquiring SAR images and preprocessing the SAR images to obtain a registration image set; wherein the registered image set comprises a plurality of registered images;
the target point selection module is used for carrying out spectrum analysis on each registration image and selecting a target point according to a spectrum analysis result;
the noise processing module is used for carrying out noise processing on the target point through differential interference processing and SVD decomposition to obtain ground deformation information;
the transmission tower deformation acquisition module is used for acquiring the current deformation and the historical deformation of the transmission tower;
the tower displacement sedimentation monitoring model building module is used for building a tower displacement sedimentation monitoring model according to the ground deformation information and the historical deformation quantity to obtain a displacement sedimentation judgment threshold;
the tower displacement sedimentation judging module is used for obtaining a tower displacement sedimentation judging result through the current deformation and the displacement sedimentation judging threshold value;
wherein, the target point selection module is further configured to:
selecting a target point by performing spectrum analysis on each registered image, wherein the target point satisfies the following conditions:
wherein N is the total number of registered images, alpha i For the spectral correlation value of pixel i in the image, beta i Is the dispersion of pixel i, α A and βA Respectively averaging the spectrum correlation value and the dispersion in the SAR image;
the noise processing module is further configured to:
performing error elimination and topographic phase interference removal on the differential interference pattern of the target point, and converting the phase of the differential interference pattern of the target point into a matrix to obtain a coefficient matrix;
and processing the coefficient matrix through SVD decomposition to eliminate the influence of interference and errors.
5. The tower displacement settlement monitoring system of claim 4, wherein the registration module is further configured to:
selecting one SAR image as a main image, and registering the SAR images at different moments with the main image to obtain registered images corresponding to the SAR images at different moments;
all the registered images are formed into a registered image set.
6. A tower displacement settlement monitoring system as claimed in claim 4, wherein the phase of the differential interferogram of the target point is expressed by the following formula:
wherein ,Δφi (x, y) is the phase of the ith differential interferogram of the target point with coordinates (x, y), phi (t) z X, y) is the interference phase of the target point with coordinates (x, y) in the main image, phi (t) f X, y) is the interference phase of the target point with coordinates (x, y) in the auxiliary image, d (t) z X, y) is the displacement of the target point with coordinates (x, y) in the main image in the radar direction, d (t) f X, y) is the displacement of the target point with coordinates of (x, y) in the auxiliary image in the radar direction, and lambda is the radar wavelength;
the phase of the differential interferogram of the target point is converted into a matrix, and the matrix is expressed by the following formula:
wherein ,is the phase matrix of the differential interferogram, delta is the coefficient of the phase matrix of the differential interferogram, +.>Is a phase matrix, and H is a coefficient matrix.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111150066.5A CN113884054B (en) | 2021-09-29 | 2021-09-29 | Tower displacement settlement monitoring method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111150066.5A CN113884054B (en) | 2021-09-29 | 2021-09-29 | Tower displacement settlement monitoring method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113884054A CN113884054A (en) | 2022-01-04 |
| CN113884054B true CN113884054B (en) | 2023-09-26 |
Family
ID=79008001
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111150066.5A Active CN113884054B (en) | 2021-09-29 | 2021-09-29 | Tower displacement settlement monitoring method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113884054B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117804405B (en) * | 2024-02-28 | 2024-05-07 | 广州煜能电气有限公司 | Transmission line tower deformation monitoring method and system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108663017A (en) * | 2018-08-13 | 2018-10-16 | 伟志股份公司 | A kind of method of ground settlement along monitoring city underground |
| CN111307109A (en) * | 2020-03-18 | 2020-06-19 | 云南电网有限责任公司带电作业分公司 | Power grid tower equipment geological settlement hidden danger monitoring method and system based on SAR satellite data |
| CN112034463A (en) * | 2020-11-06 | 2020-12-04 | 航天宏图信息技术股份有限公司 | High-precision settlement monitoring method and device for SAR data |
| CN112146622A (en) * | 2020-10-23 | 2020-12-29 | 湖南航天智远科技有限公司 | Method for monitoring geological settlement along power transmission line |
-
2021
- 2021-09-29 CN CN202111150066.5A patent/CN113884054B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108663017A (en) * | 2018-08-13 | 2018-10-16 | 伟志股份公司 | A kind of method of ground settlement along monitoring city underground |
| CN111307109A (en) * | 2020-03-18 | 2020-06-19 | 云南电网有限责任公司带电作业分公司 | Power grid tower equipment geological settlement hidden danger monitoring method and system based on SAR satellite data |
| CN112146622A (en) * | 2020-10-23 | 2020-12-29 | 湖南航天智远科技有限公司 | Method for monitoring geological settlement along power transmission line |
| CN112034463A (en) * | 2020-11-06 | 2020-12-04 | 航天宏图信息技术股份有限公司 | High-precision settlement monitoring method and device for SAR data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113884054A (en) | 2022-01-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108981559B (en) | Real-time deformation monitoring method and system based on Beidou foundation enhancement system | |
| CN109029368B (en) | Image space compensation remote sensing image/SAR image high-precision geometric positioning post-processing method | |
| KR100686287B1 (en) | Distorting Modeling method for Transforming the Presize Position of Partial/Positional information | |
| CN108761510B (en) | Method for measuring level elevation difference by using gravity field model based on terrain correction | |
| Zanutta et al. | Monitoring geodynamic activity in the Victoria Land, East Antarctica: Evidence from GNSS measurements | |
| CN108230375B (en) | Registration method of visible light image and SAR image based on structural similarity rapid robustness | |
| Zhou et al. | GPS/terrestrial 3D laser scanner combined monitoring technology for coal mining subsidence: a case study of a coal mining area in Hebei, China | |
| Lovas et al. | Terrestrial laser scanning in deformation measurements of structures | |
| CN107918127A (en) | A kind of road slope deformation detecting system and method based on vehicle-mounted InSAR | |
| CN108919304B (en) | POS error compensation method in mobile measurement system based on reference plane | |
| CN109540095A (en) | Roadbed settlement monitoring method based on satellite navigation and least square | |
| CN112882032B (en) | Method and device for dynamically monitoring geological disaster SAR in key area of gas pipeline | |
| CN113884054B (en) | Tower displacement settlement monitoring method and system | |
| CN112577470B (en) | Method and system for monitoring dynamic subsidence basin of mining area by fusion of UAV and InSAR | |
| CN113138402A (en) | RTK-based ambiguity fixing method and device and storage medium | |
| CN102508279A (en) | Method for processing GNSS (global navigation satellite system) positioning posture measuring value of satellite navigation system and GNSS positioning posture measuring instrument | |
| CN103760586B (en) | A kind of method of quick detection in GPS attitude measurements with cycle slip is repaired | |
| CN115877421A (en) | Deformation detection method and device for geological sensitive area of power transmission channel | |
| Taşçi | Analysis of dam deformation measurements with the robust and non-robust methods | |
| Stępniak et al. | Outliers and uncertainties in GNSS ZTD estimates from double-difference processing and precise point positioning | |
| Raquet et al. | Use of a Covariance Analysis Technique for Predicting Performance of Regional‐Area Differential Code and Carrier‐Phase Networks | |
| CN114463932A (en) | Non-contact construction safety distance active dynamic recognition early warning system and method | |
| CN112363166A (en) | InSAR phase unwrapping method and system based on reliable redundant network | |
| CN109977499B (en) | Beidou three-frequency/static level meter cable-stayed bridge monitoring method based on position constraint | |
| CN111751854A (en) | GNSS single-epoch double-difference integer ambiguity resolution checking method, receiver and tower crane |
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 |