CN113740852A - Method for monitoring deformation of power transmission tower by ground-based SAR (synthetic aperture radar) of artificial corner reflector - Google Patents

Method for monitoring deformation of power transmission tower by ground-based SAR (synthetic aperture radar) of artificial corner reflector Download PDF

Info

Publication number
CN113740852A
CN113740852A CN202111298360.0A CN202111298360A CN113740852A CN 113740852 A CN113740852 A CN 113740852A CN 202111298360 A CN202111298360 A CN 202111298360A CN 113740852 A CN113740852 A CN 113740852A
Authority
CN
China
Prior art keywords
phase
iron tower
deformation
sar
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.)
Granted
Application number
CN202111298360.0A
Other languages
Chinese (zh)
Other versions
CN113740852B (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.)
Hubei Chaoneng Electric Power Co ltd
Wuhan University WHU
Maintenance Branch of State Grid Hubei Electric Power Co Ltd
Original Assignee
Hubei Chaoneng Electric Power Co ltd
Wuhan University WHU
Maintenance Branch of State Grid Hubei Electric Power 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 Hubei Chaoneng Electric Power Co ltd, Wuhan University WHU, Maintenance Branch of State Grid Hubei Electric Power Co Ltd filed Critical Hubei Chaoneng Electric Power Co ltd
Priority to CN202111298360.0A priority Critical patent/CN113740852B/en
Publication of CN113740852A publication Critical patent/CN113740852A/en
Application granted granted Critical
Publication of CN113740852B publication Critical patent/CN113740852B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9094Theoretical aspects
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B15/00Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons
    • G01B15/06Measuring arrangements characterised by the use of electromagnetic waves or particle radiation, e.g. by the use of microwaves, X-rays, gamma rays or electrons for measuring the deformation in a solid

Landscapes

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

Abstract

The invention provides a method for monitoring deformation of a power transmission tower by using a ground-based SAR of an artificial corner reflector, which comprises the steps of obtaining ground-based SAR time sequence image data of a monitoring area comprising the power transmission tower and a reference point corner reflector, measuring geographic coordinates of a steel tower CR and a reference point CR, and inputting the geographic coordinates as initial data into a phase time sequence analysis software module of the steel tower CR point; the method comprises the steps of utilizing a foundation SAR to aim at an artificial corner reflector fixed on a power transmission iron tower for scanning to obtain a time sequence SAR sequence image, utilizing a geometric positioning method to extract a phase sequence of a pixel where a CR target is located from the time sequence image, and then utilizing a three-dimensional space-time phase unwrapping method to obtain a time sequence phase result of the corner reflector for iron tower deformation monitoring. A daily cycle thermal expansion model, an atmospheric error model, a breeze vibration model, a phase jump detection algorithm and a coherence iterative analysis technology are introduced into CR data processing to realize accurate processing of the deformation phase of the CR point of the iron tower.

Description

Method for monitoring deformation of power transmission tower by ground-based SAR (synthetic aperture radar) of artificial corner reflector
Technical Field
The invention relates to the technical field of aperture radar remote sensing, in particular to a method for monitoring deformation of a power transmission tower by using a ground SAR (synthetic aperture radar) of an artificial corner reflector.
Background
In recent years, the construction of transmission lines in China is accelerated, particularly the construction scale of ultra-high voltage transmission lines is continuously increased, and the total length of the transmission lines of 110 kilovolts and more is 103.34 ten thousand kilometers as the end of 2018. The geological environment of the passing area of the high-voltage transmission line is complex and various, and the safe operation of the transmission line is influenced by the hidden danger of disaster weather and geological disaster. At present, the monitoring of the running state of a power grid mainly adopts various monitoring modes such as sensor on-line monitoring, manual inspection, unmanned aerial vehicle inspection, helicopter inspection, robot inspection, infrared imaging temperature measuring instrument and the like, and the monitoring technologies are not enough in the aspects of monitoring slow micro deformation and iron tower inclination change caused by the influence of geological disasters on a power transmission iron tower, and are difficult to meet the requirements of high precision, non-contact and remote measurement.
In recent years, the ground-based SAR technology has been successfully applied to real-time monitoring and early warning of landslide deformation, but due to the problems of image resolution and scattering characteristics of iron tower targets, most power transmission iron towers are difficult to remotely identify and monitor by ground-based SAR equipment. By arranging artificial corner reflectors on some large buildings, related scientific researchers can also utilize the ground SAR technology to carry out deformation monitoring work of large bridges and super high-rise buildings, but compared with a high-voltage transmission iron tower erected in a vegetation high-coverage mountain area, the technology is difficult to be transplanted for deformation monitoring, and the main factors are that accurate identification and positioning of an iron tower target from a ground SAR image are difficult to achieve under the condition that the iron tower target is far away, and error phases caused by physical characteristics (such as thermal expansion and the influence of breeze vibration) of the iron tower target are required to be additionally eliminated.
Therefore, a monitoring technology for monitoring the target deformation of the iron tower by the ground-based SAR, which can overcome the above problems, is needed.
Disclosure of Invention
The invention aims to provide a method for monitoring the deformation of a power transmission iron tower by using a foundation SAR of an artificial corner reflector, aiming at overcoming the defects of the prior art, and the method is used for accurately identifying the target of the corner reflector on the iron tower and realizing the high-precision deformation monitoring of the power transmission iron tower.
In order to achieve the purpose, the invention adopts the following technical scheme:
the invention provides a method for monitoring deformation of a power transmission tower by using a foundation SAR (synthetic aperture radar) based on an artificial corner reflector, which comprises the following steps of:
s1, acquiring ground SAR time sequence image data of a monitoring area foundation including a power transmission iron tower and a reference point corner reflector, measuring geographic coordinates of an iron tower CR and a reference point CR, inputting the geographic coordinates as initial data into the iron tower CR point phase time sequence analysis software module, and calculating the coordinate point position of the CR in the SAR image by utilizing an SAR constellation equation;
determining the row and column number of the CR in the SAR image through the coordinate points, wherein the formula is as follows:
Figure 209955DEST_PATH_IMAGE001
(1)
wherein the content of the first and second substances,
Figure 74005DEST_PATH_IMAGE002
and
Figure 967006DEST_PATH_IMAGE003
position vectors of the sensor and the target point, respectively;
Figure 523889DEST_PATH_IMAGE004
and
Figure 384398DEST_PATH_IMAGE005
velocity vectors of the sensor and the target point, respectively;
Figure 52140DEST_PATH_IMAGE006
is a radar wavelength, obtained by imaging parameters;
Figure 799647DEST_PATH_IMAGE007
is the Doppler frequency;
Figure 527432DEST_PATH_IMAGE008
is the distance between the sensor to the target point;
s2, adopting GMT local time to carry out data diversity, generating a data set every 24 hours, wherein the reference point CR is selected as a stable point around the iron tower and used as a starting point of phase unwrapping of the interferogram, and CR point phases are obtained through the phase unwrapping, and the CR point phases comprise the following parts:
Figure 78499DEST_PATH_IMAGE009
(2)
wherein the content of the first and second substances,
Figure 156789DEST_PATH_IMAGE010
the phase of the CR point after the t interference pattern is unwrapped;
Figure 414595DEST_PATH_IMAGE011
a deformation phase of CR;
Figure 172335DEST_PATH_IMAGE012
for atmospheric phase, solving by model and interpolation;
Figure 617223DEST_PATH_IMAGE013
calculating the phase caused by thermal expansion according to the temperature and the thermal expansion coefficient;
Figure 502133DEST_PATH_IMAGE014
for phase errors caused by noise points, the noise phase is ignored due to the good ground-based SAR coherence.
Further, still include: s3, establishing a monitoring area by using an unmanned aerial vehicle LiDAR, and establishing an iron tower stress model by combining the iron tower trend and the power transmission line extension direction; the iron tower stress model takes the extending direction of a power transmission conductor as the longitudinal direction, and takes the direction which is supposed to be vertical to the extending direction of the conductor on the iron tower as the lateral direction, so as to carry out the local decomposition of the stress of the power transmission iron tower;
when the tower type and the line are tangent towers, the stressed longitudinal direction and the stressed lateral direction of the iron towers on all the straight lines are in the same coordinate system;
when the iron tower is a corner tower or a corner exists in the extension direction of the power transmission line, the longitudinal direction and the lateral direction of the iron tower are respectively set according to the front and rear iron tower directions and the wire extension direction.
Further, based on the stress model of the iron tower and the 24-hour air temperature change characteristics, carrying out the daily period modeling analysis of the thermal expansion of the iron tower, extracting the CR phase periodic change data caused by the thermal expansion of the iron tower, and establishing a corresponding deformation model as follows:
Figure 614446DEST_PATH_IMAGE015
(3)
wherein the content of the first and second substances,
Figure 543088DEST_PATH_IMAGE016
is the length of the time base of the interferogram;
Figure 740851DEST_PATH_IMAGE017
is the thermal expansion coefficient of the iron tower.
Further, SAR time series image data based on the SAR time;
when an SAR observation sampling interval is short or the atmosphere in a scene is homogeneous, a distance-phase function model is constructed by using corresponding information of a stable point outside an iron tower target area in the SAR image interferogram sequence to simulate an atmosphere phase screen and is used for eliminating the atmospheric error in the SAR image interferogram, wherein the distance-phase function model is as follows:
Figure 163873DEST_PATH_IMAGE018
(4)
wherein the content of the first and second substances,
Figure 927430DEST_PATH_IMAGE019
the unwrapped phase value of the ith stable pixel element;
Figure 902339DEST_PATH_IMAGE020
the slope distance of the ith pixel element;
Figure 712032DEST_PATH_IMAGE021
the atmospheric delay coefficient to be solved;
when the SAR observation sampling interval is long or the atmosphere in the scene is non-homogeneous, the iron tower and the CR target are firstly masked by using the strength information of the iron tower and the CR target, and the atmospheric phase of the iron tower and the CR target is interpolated by using the filtering phase of the stable point outside the iron tower target area, so that the atmospheric error of CR is eliminated.
Further, in S2, the phase unwrapping is that when the deformation phase of the iron tower exceeds pi/2 in one data set, a phase jump is generated, so that the deformation result is distorted, and a phase accumulation integral and a phase jump estimation method of the SAR time series image data are used to find out the phase jump caused by the deformation magnitude, accurately accumulate phase ambiguity integers, and accurately extract the deformation of the iron tower in a single data set.
Further, in S3, when the tower shape and the line are tangent towers, the stressed longitudinal direction and the stressed lateral direction of the iron tower on all the straight lines are in the same coordinate system, and when the iron tower is stressed by settlement, the longitudinal deformation is caused
Figure 407587DEST_PATH_IMAGE022
Decomposition of deformation from the line of sight of radar
Figure 423385DEST_PATH_IMAGE023
The conversion formula is as follows:
Figure 834775DEST_PATH_IMAGE024
(5)
wherein the content of the first and second substances,
Figure 616917DEST_PATH_IMAGE022
longitudinal deformation;
Figure 303114DEST_PATH_IMAGE025
elevation angle for ground based SAR beams;
Figure 244525DEST_PATH_IMAGE023
for solving the CR deformation, the calculation formula needs to be combined with the radar wavelength
Figure 685871DEST_PATH_IMAGE026
Comprises the following steps:
Figure 79943DEST_PATH_IMAGE027
(6)
when the iron tower is subjected to transverse traction force, the iron tower deforms laterally
Figure 648458DEST_PATH_IMAGE028
Can be decomposed by deformation of radar line of sight
Figure 444376DEST_PATH_IMAGE023
The conversion formula of (1) is as follows:
Figure 322202DEST_PATH_IMAGE029
(7)
wherein the content of the first and second substances,
Figure 203571DEST_PATH_IMAGE028
is laterally deformed.
The invention has the beneficial effects that: the method comprises the steps of utilizing a foundation SAR to aim at an artificial Corner Reflector (CR) fixed on a power transmission iron tower for scanning to obtain a time sequence SAR sequence image, utilizing a constellation equation to extract a phase sequence of a pixel where a CR target is located from the time sequence image, and then utilizing a three-dimensional space-time phase unwrapping method to obtain a deformation result of the corner reflector for iron tower deformation monitoring. A daily cycle thermal expansion model, an atmospheric error model, a breeze vibration model, a phase jump detection algorithm and a coherence iterative analysis technology are introduced into CR data processing to realize accurate acquisition of the deformation phase of the CR point of the iron tower. And finally, establishing an iron tower stress model based on the extending direction of the transmission conductor, decomposing the deformation quantity acquired by the ground SAR into a two-dimensional space, and realizing accurate calculation of the deformation and the inclination of the iron tower.
Drawings
FIG. 1 is a flow chart of a method for monitoring deformation of a transmission tower by ground-based SAR monitoring of an artificial corner reflector according to the present invention;
FIG. 2 is an intensity plot of an SAR image of a monitored area containing a pylon and a reference point corner reflector;
FIG. 3 is a schematic diagram of a three-dimensional LiDAR point cloud for a pylon.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to the accompanying drawings. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
Referring to fig. 1, a method for monitoring deformation of a power transmission tower by using a ground-based SAR of an artificial corner reflector includes the following steps:
s1, acquiring ground SAR time sequence image data of a monitoring area foundation including a power transmission iron tower and a reference point corner reflector, measuring geographic coordinates of an iron tower CR (artificial corner reflector) and the reference point CR, inputting the data into the iron tower CR point phase time sequence analysis software module as initial data, and calculating coordinate point positions of all CRs in the SAR image by utilizing an SAR constellation equation;
wherein, the iron tower CR is a CR placed on the iron tower; the reference point CR is CR placed at the reference point.
Determining the row and column number of the CR in the SAR image through the coordinate points, wherein the formula is as follows:
Figure 44619DEST_PATH_IMAGE001
(1)
wherein the content of the first and second substances,
Figure 960622DEST_PATH_IMAGE002
and
Figure 743770DEST_PATH_IMAGE003
position vectors of the sensor and the target point, respectively;
Figure 112435DEST_PATH_IMAGE004
and
Figure 19824DEST_PATH_IMAGE005
velocity vectors of the sensor and the target point, respectively;
Figure 524754DEST_PATH_IMAGE006
is a radar wavelength, obtained by imaging parameters;
Figure 744383DEST_PATH_IMAGE007
is the Doppler frequency;
Figure 600344DEST_PATH_IMAGE008
is the distance between the sensor to the target point;
that is, the combination parameter file can be
Figure 314353DEST_PATH_IMAGE030
And
Figure 939369DEST_PATH_IMAGE031
the row and column numbers of the CR dots are converted.
The ground-based SAR device comprises a transmitting antenna (a part of linear scanning device comprising an antenna orbit);
the CR is fixed on the power transmission iron tower and is properly adjusted according to the scanning direction of the ground SAR so that the axis direction of the CR is aligned to the direction of the ground SAR wave beam;
s2, performing data diversity by adopting GMT local time, and generating a data set every 24 hours, wherein the data for performing the data diversity is acquired by ground-based SAR scanning;
referring to fig. 2, the reference point CR should be selected as a stable point around the iron tower, and used as a starting point for phase unwrapping of the interferogram, and a CR point phase is obtained through phase unwrapping, where the CR point phase includes the following parts:
Figure 329899DEST_PATH_IMAGE009
(2)
wherein the content of the first and second substances,
Figure 407577DEST_PATH_IMAGE010
the phase of the CR point after the t interference pattern is unwrapped;
Figure 659698DEST_PATH_IMAGE011
a deformation phase of CR;
Figure 139220DEST_PATH_IMAGE012
for atmospheric phase, solving by model and interpolation;
Figure 700652DEST_PATH_IMAGE013
calculating the phase caused by thermal expansion according to the temperature and the thermal expansion coefficient;
Figure 265625DEST_PATH_IMAGE014
for phase errors caused by noise points, the noise phase is ignored due to the good ground-based SAR coherence.
Wherein the function of the formula (2) is to obtain
Figure 321437DEST_PATH_IMAGE032
By the obtained
Figure 921046DEST_PATH_IMAGE032
And judging the deformation degree of the iron tower according to the value.
Further comprising: s3, establishing a monitoring area by using an unmanned aerial vehicle LiDAR (laser radar technology), and establishing an iron tower stress model by combining the iron tower trend and the power transmission line extension direction; the iron tower stress model takes the extending direction of a power transmission conductor as the longitudinal direction, and takes the direction which is supposed to be vertical to the extending direction of the conductor on the iron tower as the lateral direction, so as to carry out the local decomposition of the stress of the power transmission iron tower;
the monitoring area comprises a three-dimensional point cloud model of the power transmission tower;
when the tower type and the line are tangent towers, the stressed longitudinal direction and the stressed lateral direction of the iron towers on all the straight lines are in the same coordinate system;
when the iron tower is a corner tower or a corner exists in the extension direction of the power transmission line, the longitudinal direction and the lateral direction of the iron tower are respectively set according to the front and rear iron tower directions and the wire extension direction.
Further comprising: based on the stress model of the iron tower and the 24-hour air temperature change characteristics, carrying out modeling analysis on the daily period of thermal expansion of the iron tower, extracting CR phase periodic change data caused by the thermal expansion of the iron tower, and establishing a corresponding deformation model as follows:
Figure 653378DEST_PATH_IMAGE015
(3)
wherein the content of the first and second substances,
Figure 705648DEST_PATH_IMAGE016
for the interferenceThe time base length of the graph;
Figure 302501DEST_PATH_IMAGE017
is the thermal expansion coefficient of the iron tower.
Thus, in the formula (2) is obtained,
Figure 756616DEST_PATH_IMAGE033
the value of (c).
Based on the SAR time sequence image data, the situation that when the SAR observation sampling interval is short or the atmosphere in the scene is homogeneous and the situation that when the SAR observation sampling interval is long or the atmosphere in the scene is inhomogeneous can be obtained;
when an SAR observation sampling interval is short or the atmosphere in a scene is homogeneous, a distance-phase function model is constructed by using corresponding information of a stable point outside an iron tower target area in the SAR image interferogram sequence to simulate an atmosphere phase screen and is used for eliminating the atmospheric error in the SAR image interferogram, wherein the distance-phase function model is as follows:
Figure 659850DEST_PATH_IMAGE018
(4)
wherein the content of the first and second substances,
Figure 199416DEST_PATH_IMAGE019
the unwrapped phase value of the ith stable pixel element;
Figure 597030DEST_PATH_IMAGE020
the slope distance of the ith pixel element;
Figure 905652DEST_PATH_IMAGE021
the atmospheric delay coefficient to be solved;
thus, in the formula (2) obtained by solving,
Figure 979787DEST_PATH_IMAGE034
an atmospheric phase value;
when the SAR observation sampling interval is long or the atmosphere in the scene is non-homogeneous, the iron tower and the CR target are firstly masked by using the strength information of the iron tower and the CR target, and the atmospheric phase of the iron tower and the CR target is interpolated by using the filtering phase of the stable point outside the iron tower target area, so that the atmospheric error of CR is eliminated.
In the embodiment of the present invention, since in the formula (2),
Figure 741070DEST_PATH_IMAGE035
and the CR point phase after the interference pattern unwrapping is obtained by phase unwrapping the original interference pattern.
Atmospheric phase value
Figure 676796DEST_PATH_IMAGE036
Can be obtained by the formula (4);
phase value caused by thermal expansion
Figure 839924DEST_PATH_IMAGE037
Can be obtained by the formula (3);
phase error caused by noise point
Figure 84961DEST_PATH_IMAGE038
Can be omitted;
the deformed phase value of CR can be obtained by the formula (2)
Figure 599119DEST_PATH_IMAGE039
In S2, the phase unwrapping process includes: and when the deformation phase of the iron tower exceeds pi/2 in one data set, phase jump is generated to distort the deformation result, and the phase jump caused by the deformation magnitude is found out by adopting the phase accumulated integral and phase jump estimation method of the SAR time sequence image data, the phase ambiguity integer is accurately accumulated, and the accurate extraction of the deformation of the iron tower in a single data set is realized.
Because the stress of the iron tower mainly comprises the tower shape and the line of the iron tower being a tangent tower and the iron tower being a corner tower or the extension direction of the power transmission line having a corner,
referring to fig. 3, in particular,when the tower type and the line are tangent towers, the stressed iron towers on all the straight lines have longitudinal and lateral directions in the same coordinate system, and when the iron towers are subjected to settlement stress, the iron towers deform longitudinally
Figure 338536DEST_PATH_IMAGE022
Decomposition of deformation from the line of sight of radar
Figure 356170DEST_PATH_IMAGE023
The conversion formula is as follows:
Figure 37687DEST_PATH_IMAGE024
(5)
wherein the content of the first and second substances,
Figure 773562DEST_PATH_IMAGE022
longitudinal deformation;
Figure 972462DEST_PATH_IMAGE025
elevation angle for ground based SAR beams;
Figure 982619DEST_PATH_IMAGE023
for solving the CR deformation, the calculation formula needs to be combined with the radar wavelength
Figure 38300DEST_PATH_IMAGE026
Comprises the following steps:
Figure 792629DEST_PATH_IMAGE027
(6)
when the iron tower is subjected to transverse traction force, the iron tower deforms laterally
Figure 342690DEST_PATH_IMAGE028
Can be decomposed by deformation of radar line of sight
Figure 131655DEST_PATH_IMAGE023
The conversion formula of (1) is as follows:
Figure 764761DEST_PATH_IMAGE029
(7)
wherein the content of the first and second substances,
Figure 599862DEST_PATH_IMAGE028
is laterally deformed.
When the iron tower is a corner tower or the extension direction of the power transmission line has a corner, the deformation does not belong to the technical content of the invention and is not further elaborated.
The above-mentioned embodiments only express the embodiments of the present invention, and the description thereof is more specific and detailed, but not construed as limiting the scope of the present invention. It should be noted that, for a person skilled in the art, several variations and modifications can be made without departing from the inventive concept, which falls within the scope of the present invention. Therefore, the protection scope of the present patent shall be subject to the appended claims.

Claims (6)

1. The method for monitoring the deformation of the power transmission tower by the foundation SAR of the artificial corner reflector is characterized by comprising the following steps of:
s1, acquiring ground SAR time sequence image data of a monitoring area foundation including a power transmission iron tower and a reference point corner reflector, measuring geographic coordinates of an iron tower CR and a reference point CR, inputting the geographic coordinates as initial data into the iron tower CR point phase time sequence analysis software module, and calculating the coordinate point position of the CR in the SAR image by utilizing an SAR constellation equation;
determining the row and column number of the CR in the SAR image through the coordinate points, wherein the formula is as follows:
Figure 449731DEST_PATH_IMAGE001
(1)
wherein the content of the first and second substances,
Figure 924575DEST_PATH_IMAGE002
and
Figure 562361DEST_PATH_IMAGE003
position vectors of the sensor and the target point, respectively;
Figure 9523DEST_PATH_IMAGE004
and
Figure 543404DEST_PATH_IMAGE005
velocity vectors of the sensor and the target point, respectively;
Figure 912068DEST_PATH_IMAGE006
is a radar wavelength, obtained by imaging parameters;
Figure 71654DEST_PATH_IMAGE007
is the Doppler frequency;
Figure 248688DEST_PATH_IMAGE008
is the distance between the sensor to the target point;
s2, adopting GMT local time to carry out data diversity, generating a data set every 24 hours, wherein the reference point CR is selected as a stable point around the iron tower and used as a starting point of phase unwrapping of the interferogram, and CR point phases are obtained through the phase unwrapping, and the CR point phases comprise the following parts:
Figure 78104DEST_PATH_IMAGE009
(2)
wherein the content of the first and second substances,
Figure 806501DEST_PATH_IMAGE010
the phase of the CR point after the t interference pattern is unwrapped;
Figure 645144DEST_PATH_IMAGE011
a deformation phase of CR;
Figure 129215DEST_PATH_IMAGE012
for atmospheric phase, by means of a dieType and interpolation to solve;
Figure 395111DEST_PATH_IMAGE013
calculating the phase caused by thermal expansion according to the temperature and the thermal expansion coefficient;
Figure 269527DEST_PATH_IMAGE014
for phase errors caused by noise points, the noise phase is ignored due to the good ground-based SAR coherence.
2. The method for monitoring deformation of power transmission towers through ground-based SAR monitoring of artificial corner reflectors according to claim 1, further comprising:
s3, establishing a monitoring area by using an unmanned aerial vehicle LiDAR, and establishing an iron tower stress model by combining the iron tower trend and the power transmission line extension direction; the iron tower stress model takes the extending direction of a power transmission conductor as the longitudinal direction, and takes the direction which is supposed to be vertical to the extending direction of the conductor on the iron tower as the lateral direction, so as to carry out the local decomposition of the stress of the power transmission iron tower;
when the tower type and the line are tangent towers, the stressed longitudinal direction and the stressed lateral direction of the iron towers on all the straight lines are in the same coordinate system;
when the iron tower is a corner tower or a corner exists in the extension direction of the power transmission line, the longitudinal direction and the lateral direction of the iron tower are respectively set according to the front and rear iron tower directions and the wire extension direction.
3. The method for monitoring deformation of power transmission towers through ground-based SAR monitoring of artificial corner reflectors according to claim 2, wherein the method comprises the following steps:
based on the stress model of the iron tower and the 24-hour air temperature change characteristics, carrying out modeling analysis on the daily period of thermal expansion of the iron tower, extracting CR phase periodic change data caused by the thermal expansion of the iron tower, and establishing a corresponding deformation model as follows:
Figure 787227DEST_PATH_IMAGE015
(3)
wherein the content of the first and second substances,
Figure 266750DEST_PATH_IMAGE016
is the length of the time base of the interferogram;
Figure 31443DEST_PATH_IMAGE017
is the thermal expansion coefficient of the iron tower.
4. The method for monitoring deformation of power transmission towers through ground-based SAR monitoring of artificial corner reflectors according to claim 1, wherein the method comprises the following steps:
based on the SAR time sequence image data, when an SAR observation sampling interval is short or the atmosphere in a scene is homogeneous, a distance-phase function model is constructed by using corresponding information of a stable point outside an iron tower target region in the SAR image interferogram sequence to simulate an atmosphere phase screen so as to eliminate the atmosphere error in the SAR image interferogram, wherein the distance-phase function model is as follows:
Figure 393154DEST_PATH_IMAGE018
(4)
wherein the content of the first and second substances,
Figure 448966DEST_PATH_IMAGE019
the unwrapped phase value of the ith stable pixel element;
Figure 782996DEST_PATH_IMAGE020
the slope distance of the ith pixel element;
Figure 718591DEST_PATH_IMAGE021
the atmospheric delay coefficient to be solved;
when the SAR observation sampling interval is long or the atmosphere in the scene is non-homogeneous, the iron tower and the CR target are firstly masked by using the strength information of the iron tower and the CR target, and the atmospheric phase of the iron tower and the CR target is interpolated by using the filtering phase of the stable point outside the iron tower target area, so that the atmospheric error of CR is eliminated.
5. The method for monitoring deformation of power transmission towers through ground-based SAR monitoring of artificial corner reflectors according to claim 1, wherein the method comprises the following steps: in the step S2, the phase unwrapping is that when the deformation phase of the iron tower exceeds pi/2 in one data set, phase jump is generated, so that the deformation result is distorted, and a phase accumulated integral and phase jump estimation method of the SAR time series image data is adopted to find out the phase jump caused by the deformation magnitude, accurately accumulate the phase ambiguity integer, and accurately extract the deformation of the iron tower in the single data set.
6. The method for monitoring deformation of a pylon according to claim 2, wherein in S3:
when the tower type and the line are tangent towers, the stressed iron towers on all the straight lines have longitudinal and lateral directions in the same coordinate system, and when the iron towers are subjected to settlement stress, the iron towers deform longitudinally
Figure 567598DEST_PATH_IMAGE022
Decomposition of deformation from the line of sight of radar
Figure 161521DEST_PATH_IMAGE023
The conversion formula is as follows:
Figure 615637DEST_PATH_IMAGE024
(5)
wherein the content of the first and second substances,
Figure 518870DEST_PATH_IMAGE022
longitudinal deformation;
Figure 58436DEST_PATH_IMAGE025
elevation angle for ground based SAR beams;
Figure 283199DEST_PATH_IMAGE023
for solving the CR deformation, the calculation formula needs to be combined with the radar wavelength
Figure 591821DEST_PATH_IMAGE026
Comprises the following steps:
Figure 665956DEST_PATH_IMAGE027
(6)
when the iron tower is subjected to transverse traction force, the iron tower deforms laterally
Figure 427239DEST_PATH_IMAGE028
Can be decomposed by deformation of radar line of sight
Figure 628544DEST_PATH_IMAGE023
The conversion formula of (1) is as follows:
Figure 791672DEST_PATH_IMAGE029
(7)
wherein the content of the first and second substances,
Figure 36708DEST_PATH_IMAGE028
is laterally deformed.
CN202111298360.0A 2021-11-04 2021-11-04 Method for monitoring deformation of power transmission tower by ground-based SAR (synthetic aperture radar) of artificial corner reflector Active CN113740852B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111298360.0A CN113740852B (en) 2021-11-04 2021-11-04 Method for monitoring deformation of power transmission tower by ground-based SAR (synthetic aperture radar) of artificial corner reflector

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111298360.0A CN113740852B (en) 2021-11-04 2021-11-04 Method for monitoring deformation of power transmission tower by ground-based SAR (synthetic aperture radar) of artificial corner reflector

Publications (2)

Publication Number Publication Date
CN113740852A true CN113740852A (en) 2021-12-03
CN113740852B CN113740852B (en) 2022-04-12

Family

ID=78727348

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111298360.0A Active CN113740852B (en) 2021-11-04 2021-11-04 Method for monitoring deformation of power transmission tower by ground-based SAR (synthetic aperture radar) of artificial corner reflector

Country Status (1)

Country Link
CN (1) CN113740852B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114966685A (en) * 2022-05-24 2022-08-30 中国水利水电科学研究院 Dam deformation monitoring and predicting method based on InSAR and deep learning

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0815426A (en) * 1994-06-30 1996-01-19 Nec Corp Terrain alteration extracting device using interference type synthetic aperture radar
CN101770027A (en) * 2010-02-05 2010-07-07 河海大学 Ground surface three-dimensional deformation monitoring method based on InSAR and GPS data fusion
CN110456346A (en) * 2019-06-28 2019-11-15 深圳市水务规划设计院股份有限公司 A kind of electric power pylon inclination monitoring method based on InSAR technology
CN110596702A (en) * 2019-07-29 2019-12-20 武汉大学 River basin reservoir bank deformation monitoring method based on ground stable platform SAR

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0815426A (en) * 1994-06-30 1996-01-19 Nec Corp Terrain alteration extracting device using interference type synthetic aperture radar
CN101770027A (en) * 2010-02-05 2010-07-07 河海大学 Ground surface three-dimensional deformation monitoring method based on InSAR and GPS data fusion
CN110456346A (en) * 2019-06-28 2019-11-15 深圳市水务规划设计院股份有限公司 A kind of electric power pylon inclination monitoring method based on InSAR technology
CN110596702A (en) * 2019-07-29 2019-12-20 武汉大学 River basin reservoir bank deformation monitoring method based on ground stable platform SAR

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
陈怡曲: "基于InSAR的形变监测技术研究", 《中国优秀博硕士学位论文全文数据库(硕士)•基础科学辑》 *
陶李 等: "基于InSAR技术在输变电工程形变监测中的应用", 《电力勘测设计》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114966685A (en) * 2022-05-24 2022-08-30 中国水利水电科学研究院 Dam deformation monitoring and predicting method based on InSAR and deep learning
CN114966685B (en) * 2022-05-24 2023-04-07 中国水利水电科学研究院 Dam deformation monitoring and predicting method based on InSAR and deep learning

Also Published As

Publication number Publication date
CN113740852B (en) 2022-04-12

Similar Documents

Publication Publication Date Title
Tarchi et al. Remote monitoring of buildings using a ground-based SAR: application to cultural heritage survey
US20200103530A1 (en) Method for extracting elevation control point with assistance of satellite laser altimetry data
Pieraccini et al. Terrain mapping by ground-based interferometric radar
CN110456346B (en) Power transmission tower inclination monitoring method based on InSAR technology
CN112034460B (en) Circular arc aperture radar imaging method and radar based on antenna phase pattern compensation
US20100045513A1 (en) Stability monitoring using synthetic aperture radar
CN109100723B (en) High-altitude wind inversion method based on Doppler weather radar data
KR102151362B1 (en) Image decoding apparatus based on airborn using polar coordinates transformation and method of decoding image using the same
CN110441770A (en) Three-dimensional deformation measurement method based on multi-section MIMO-SAR joint observation
CN109613583A (en) The passive object localization method of the time difference is surveyed based on the direction finding of Dan Xingyu earth station and joint
Noferini et al. DEM by ground-based SAR interferometry
CN110940966B (en) Laser footprint plane positioning method based on laser height measurement satellite footprint image
CN111398959A (en) InSAR time sequence earth surface deformation monitoring method based on earth surface stress strain model
Nico et al. Generation of digital terrain models with a ground-based SAR system
CN113740852B (en) Method for monitoring deformation of power transmission tower by ground-based SAR (synthetic aperture radar) of artificial corner reflector
CN112444188B (en) Multi-view InSAR sea wall high-precision three-dimensional deformation measurement method
CN112904337A (en) Side slope deformation time sequence monitoring method based on Offset Tracking technology
CN115407289A (en) Satellite-borne SAR geometric calibration method based on control points provided by wind driven generator
CN111650570A (en) Three-dimensional atmospheric correction method and system for ground-based interference radar
CN112379342A (en) Echo simulation and echo characteristic parameter precision estimation method for satellite-borne cloud detection radar
CN111487621A (en) Sea surface flow field inversion method based on radar image and electronic equipment
Hosseiny et al. Structural displacement monitoring using ground-based synthetic aperture radar
CN106483511B (en) A kind of Electromagnetic Simulation method of transmission line of electricity in satellite-borne synthetic aperture radar imaging
CN113064161B (en) Wave spectrometer cross spectrum calculation method based on double sub-pulse reconstruction
CN114239379A (en) Transmission line geological disaster analysis method and system based on deformation detection

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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 430000, No. 60-65, Erhe village, wulidun, Hanyang District, Wuhan City, Hubei Province

Applicant after: Super high voltage company of State Grid Hubei Electric Power Co.,Ltd.

Applicant after: HUBEI CHAONENG ELECTRIC POWER Co.,Ltd.

Applicant after: WUHAN University

Address before: 430000, No. 60-65, Erhe village, wulidun, Hanyang District, Wuhan City, Hubei Province

Applicant before: STATE GRID HUBEI ELECTRIC POWER CO., LTD. MAINTENANCE Co.

Applicant before: HUBEI CHAONENG ELECTRIC POWER Co.,Ltd.

Applicant before: WUHAN University

GR01 Patent grant
GR01 Patent grant