CN111487619A - Method for detecting ground inclination angle of transmission tower by utilizing SAR-X high-resolution bunching mode - Google Patents

Method for detecting ground inclination angle of transmission tower by utilizing SAR-X high-resolution bunching mode Download PDF

Info

Publication number
CN111487619A
CN111487619A CN202010320737.7A CN202010320737A CN111487619A CN 111487619 A CN111487619 A CN 111487619A CN 202010320737 A CN202010320737 A CN 202010320737A CN 111487619 A CN111487619 A CN 111487619A
Authority
CN
China
Prior art keywords
ground
corridor
tower
coordinates
main image
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.)
Pending
Application number
CN202010320737.7A
Other languages
Chinese (zh)
Inventor
黄然
姜虹云
卢勇
沈志
金晶
周仿荣
文刚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Electric Power Research Institute of Yunnan Power Grid Co Ltd
Original Assignee
Electric Power Research Institute of Yunnan Power Grid 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 Electric Power Research Institute of Yunnan Power Grid Co Ltd filed Critical Electric Power Research Institute of Yunnan Power Grid Co Ltd
Priority to CN202010320737.7A priority Critical patent/CN111487619A/en
Publication of CN111487619A publication Critical patent/CN111487619A/en
Pending legal-status Critical Current

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/9021SAR image post-processing techniques
    • G01S13/9023SAR image post-processing techniques combined with interferometric techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021SAR image post-processing techniques
    • G01S13/9027Pattern recognition for feature extraction
    • 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
    • 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/414Discriminating targets with respect to background clutter
    • 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/418Theoretical aspects

Landscapes

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

Abstract

The application provides a method for detecting the ground inclination angle of a transmission tower by utilizing an SAR-X high-resolution bunching mode. The method comprises the following steps: acquiring a main image intensity map of a power transmission corridor and a power grid tower which take a tower coordinate as a center and have a radius of more than 5 times of a root length of the tower; according to a segmentation threshold value, segmenting a ground corridor part on the main image intensity image; calculating the backscattering coefficient sum of the ground corridor part; extracting the pole tower ground plane part in the main image intensity graph according to the backscattering coefficient sum, and outputting the coordinates of the ground corridor pixel points from the main image intensity graph; converting the coordinates of the pixel points of the ground corridor to the coordinates of the geocentric system; calculating the final coordinates of the ground corridor pixel points according to the coordinates of the ground center system; and calculating the ground inclination angle of the transmission tower according to the final coordinate. The utility model belongs to the technical field of power grid transmission line characteristic quantity monitoring, need not to erect the instrument, use not influenced by topography, have convenient and fast, characteristics that measurement accuracy is high.

Description

Method for detecting ground inclination angle of transmission tower by utilizing SAR-X high-resolution bunching mode
Technical Field
The application relates to the technical field of power grid transmission line characteristic quantity monitoring, in particular to a method for detecting a ground inclination angle of a transmission tower by utilizing an SAR-X high-resolution bunching mode.
Background
The ground inclination angle is an important parameter for protecting the overhead line from lightning strike, the shielding effect of the ground on the side phase conductor during the lightning strike is represented, and according to the shielding failure electrical geometric model principle, the corresponding impact distance is increased along with the increase of the lightning current amplitude, and the shielding failure exposed surface is reduced. When the lightning current value reaches the shielding failure critical lightning current, the shielding failure exposure arc surface is zero, at the moment, the lightning guide cannot develop towards the lead, the probability that the lightning current directly hits the lead is very low, and the possibility of hitting the ground is very high. The smaller the ground inclination angle, the better the shielding effect of the ground on the wire when the lightning stroke occurs, and the weaker the ground inclination angle is, the weaker the ground inclination angle is. The ground inclination angle is different according to the landform where the transmission tower is located and the line trend. In the calculation of the overhead line shielding failure critical lightning current, the ground inclination angle is a parameter which must be acquired first, so that the accurate measurement of the ground inclination angle is particularly important.
The development of Synthetic Aperture Radar (SAR), which is an active microblog sensor, provides a new technology for solving the problem, and the SAR acquires scattering information of a target by transmitting and receiving specific electromagnetic waves. The SAR can be free from the limitation of cloud, fog, sleet, night and other conditions, can monitor the target all weather and all day long, and has obvious superiority under the condition of large-scale natural disasters. With the continuous development of the SAR technology, the application of SAR remote sensing is deepened and expanded continuously, the SAR remote sensing system can be used for identifying and classifying targets such as roads, bridges, buildings, vehicles and the like, can also be used for forest vegetation and ice and snow monitoring, topographic survey, urban overall monitoring, disaster early warning and the like, and has great application potential in the fields of civil use, national defense construction and the like.
The total station used in the prior art is a surveying and mapping tool playing a key role in surveying and mapping, is simple to use, and has the characteristics of high speed measurement and high precision. However, in the practical application process, the test accuracy is affected when vibration occurs, and when the instrument is used under strong sunlight, the instrument needs to be put on an umbrella and a sun cover is arranged, and after the instrument is used, a specially-assigned person needs to be arranged for protection.
Disclosure of Invention
In order to solve the problems, the application provides a method for detecting the ground inclination angle of the transmission tower by utilizing an SAR-X high-resolution bunching mode, so as to solve the problems of more use limitation and difficult maintenance of a surveying and mapping instrument in the prior art.
In order to achieve the purpose, the method is realized by the following technical scheme:
randomly selecting a two-dimensional SAR image of a power transmission corridor at a certain time;
obtaining a main image intensity image by interfering the two-dimensional SAR image of the power transmission corridor;
dividing a ground corridor part on the main image intensity image according to a division threshold value;
calculating the backscattering coefficient sum of the ground corridor part;
extracting a pole tower ground plane part in the main image intensity graph according to the backscattering coefficient sum, and outputting a ground corridor pixel point coordinate on the main image intensity graph;
converting the coordinates of the pixel points of the ground corridor to the coordinates of a geocentric system;
calculating the final coordinate of the ground corridor pixel point according to the geocentric system coordinate;
and calculating the ground inclination angle of the position where the tower is located according to the final coordinates of the ground corridor pixel points.
Further, the two-dimensional SAR image of the power transmission corridor takes a pole tower coordinate as a center, and the length of a pole tower root opening which is more than 5 times of the pole tower root opening is taken as a radius.
Further, the step of segmenting the ground corridor part on the main image intensity map according to the segmentation threshold value comprises the following specific steps:
calculating the intensity average value of each point of the main image intensity map;
marking the point with the intensity smaller than the average value in the main image intensity image as a noise point;
all the noise points are removed to obtain calculation points;
and calculating the intensity average value of all the calculation points, and determining the intensity average value as a segmentation threshold value.
Dividing the main image intensity map according to the segmentation threshold;
and determining the part of the main image intensity image with the intensity larger than the segmentation threshold value as a ground corridor part, and determining the part of the main image intensity image with the intensity smaller than or equal to the segmentation threshold value as a background part.
Further, the step of calculating the backscattering coefficient sum of the ground corridor part comprises the following steps:
obtaining structural parameters of the ground corridor part: the center point coordinates and the radius length of the tower, wherein the radius length is 5 times of the tower root opening length at minimum;
determining the backscattering coefficient of all circular circumferences of the ground corridor part with integral multiple of the pole tower root opening length as the radius to the maximum circular circumference according to the structural parameters;
the sum of the backscatter coefficients for all the terrestrial corridor sections is calculated.
Further, the specific step of outputting the coordinates of the ground corridor pixel points on the main image intensity map includes:
marking the pixel points of the ground corridor part as a first color, and outputting a marking point coordinate on the main image intensity graph;
traversing the coordinates of the mark points according to a preset template;
determining points, the matching degree of which is in accordance with a preset matching threshold value, of the preset template as fine extraction points;
marking the fine extraction point as a second color, and outputting the fine extraction point coordinates on the main image intensity map.
Further, the step of converting the coordinates of the pixel points of the ground corridor to the coordinates of the geocentric system further comprises:
and fitting a satellite orbit equation according to the 12 position vectors and the velocity vector of the required satellite in the geocentric coordinate system:
Figure BDA0002461289110000031
Figure BDA0002461289110000032
Figure BDA0002461289110000033
wherein Xsf, Ysf and XZf are satellite position vectors, Vxsf, VYsf and VZsf are satellite velocity vectors, i is the number of rows of pixel points on an image, PRF is pulse repetition frequency, t is time interval, a0、a1、a3、b1、b2、b3、c1、c2、c3A set of fixed unknown parameters that need to be solved for.
Further, the specific step of establishing the functional relationship between the slant range image coordinate and the earth center system coordinate is as follows;
determining a distance equation according to:
F1=(Xsf-Xgf)2+(Ysf-Ygf)2+(Zsf-Zgf)2-R2
wherein Xgf, Ygf and Zgf are position vectors of ground points in a geocentric coordinate system, and R is a slope distance;
determining the doppler equation according to:
Figure BDA0002461289110000034
wherein V is the wave velocity of incident wave, lambda is the wavelength of incident wave, fDCIs the Doppler frequency shift;
determining an earth ellipsoid equation according to:
Figure BDA0002461289110000035
wherein N is the semiaxis of the ellipse of the earth, N (1-e)2) Is the short semi-axis of the earth ellipse, and h is the elevation of the ground point.
Further, according to the distance equation, the Doppler equation and the earth ellipsoid equation, calculating the position of the pixel point of the ground corridor part in the geocentric coordinate system to obtain the following formula:
Δx=B-1·F(Xk)
wherein, Deltax is the correction amount,
Figure BDA0002461289110000036
b is F (X)k) The value obtained by deriving Xgf, Ygf, Zgf is the B value according to the following formula:
Figure BDA0002461289110000041
calculating an absolute value of the correction amount;
if the absolute value of the correction amount is larger than or equal to the preset correction amount threshold, carrying out correction calculation again until the absolute value of the correction amount is smaller than the preset correction amount threshold;
and determining that the corresponding Xk value when the absolute value of the correction amount is smaller than the preset correction amount threshold is the final coordinate of the pixel point of the ground corridor part.
Compared with the prior art, the beneficial effects of this application are:
the method is a method for detecting the ground inclination angle of a transmission tower by utilizing an SAR-X high-resolution bunching mode, and compared with the prior art, 1. the method does not need to erect related instruments on a surveying and mapping site, is not limited by terrain and climate, and saves labor cost; 2. the method has the advantages of high measurement precision, high accuracy and wide test range.
Drawings
In order to more clearly explain the technical solution of the present application, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious to those skilled in the art that other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a flowchart of a method for detecting a transmission tower ground inclination angle by using an SAR-X high resolution bunching mode.
Detailed Description
Reference will now be made in detail to embodiments, examples of which are illustrated in the accompanying drawings. When the following description refers to the accompanying drawings, like numbers in different drawings represent the same or similar elements unless otherwise indicated. The embodiments described in the following examples do not represent all embodiments consistent with the present application. But merely as exemplifications of systems and methods consistent with certain aspects of the application, as recited in the claims.
Referring to fig. 1, a flowchart of a method for detecting a transmission tower ground inclination angle by using an SAR-X high resolution bunching mode is shown. The technical scheme of the application is as follows:
randomly selecting two power transmission corridor two-dimensional SAR images of a certain time phase with the radius of 1 kilometer and taking the coordinates of a line tower as the center; and obtaining a main image intensity image by interfering the two-dimensional SAR images of the power transmission corridor.
The steps of determining the segmentation threshold are as follows:
calculating the intensity average value of each point of the main image intensity map;
marking the point with the intensity smaller than the average value in the main image intensity image as a noise point;
all the noise points are removed to obtain calculation points;
and calculating the intensity average value of all the calculation points, and determining the intensity average value as a segmentation threshold value.
According to a segmentation threshold value, the specific steps of segmenting the ground corridor part on the main image intensity map comprise:
dividing the main image intensity map according to the segmentation threshold;
and determining the part of the main image intensity image with the intensity larger than the segmentation threshold value as a ground corridor part, and determining the part of the main image intensity image with the intensity smaller than or equal to the segmentation threshold value as a background part.
Calculating the backscattering coefficient sum of the ground corridor part, and the specific steps are as follows:
obtaining structural parameters of the ground corridor part: the center point coordinates and the radius length of the tower, wherein the radius length is 5 times of the tower root opening length at minimum;
determining the backscattering coefficient of all circular circumferences of the ground corridor part with integral multiple of the pole tower root opening length as the radius to the maximum circular circumference according to the structural parameters;
the sum of the backscatter coefficients for all the terrestrial corridor sections is calculated.
According to the backscattering coefficient sum, extracting a pole tower ground plane part in the main image intensity graph, and outputting coordinates of ground corridor pixel points on the main image intensity graph, wherein the method specifically comprises the following steps:
marking the pixel points of the ground corridor part as a first color, and outputting a marking point coordinate on the main image intensity graph;
traversing the coordinates of the mark points according to a preset template;
determining points, the matching degree of which is in accordance with a preset matching threshold value, of the preset template as fine extraction points;
marking the fine extraction point as a second color, and outputting the fine extraction point coordinate, namely the ground corridor pixel point coordinate, on the main image intensity map.
Converting the coordinates to geocentric system coordinates:
and fitting a satellite orbit equation according to the 12 position vectors and the velocity vector of the required satellite in the geocentric coordinate system:
Figure BDA0002461289110000051
Figure BDA0002461289110000052
Figure BDA0002461289110000053
wherein Xsf, Ysf and XZf are satellite position vectors, Vxsf, VYsf and VZsf are satellite velocity vectors, i is the number of rows of pixel points on an image, PRF is pulse repetition frequency, t is time interval, a0、a1、a3、b1、b2、b3、c1、c2、c3A set of fixed unknown parameters that need to be solved for.
The specific steps of establishing the functional relationship between the slant range image coordinates and the coordinates of the earth center system are as follows:
determining a distance equation according to:
F1=(Xsf-Xgf)2+(Ysf-Ygf)2+(Zsf-Zgf)2-R2
wherein Xgf, Ygf and Zgf are position vectors of ground points in a geocentric coordinate system, and R is a slope distance;
determining the doppler equation according to:
Figure BDA0002461289110000061
wherein V is the wave velocity of incident wave, lambda is the wavelength of incident wave, fDCIs the Doppler frequency shift;
determining an earth ellipsoid equation according to:
Figure BDA0002461289110000062
wherein N is the semiaxis of the ellipse of the earth, N (1-e)2) Is the short semi-axis of the earth ellipse, and h is the elevation of the ground point.
Calculating the position of the pixel point of the ground corridor part in the geocentric coordinate system according to the distance equation, the Doppler equation and the earth ellipsoid equation to obtain the following formula:
Δx=B-1·F(Xk)
wherein, Deltax is the correction amount,
Figure BDA0002461289110000063
b is F (X)k) The value obtained by deriving Xgf, Ygf, Zgf is the B value according to the following formula:
Figure BDA0002461289110000064
calculating an absolute value of the correction amount;
if the absolute value of the correction amount is larger than or equal to the preset correction amount threshold, carrying out correction calculation again until the absolute value of the correction amount is smaller than the preset correction amount threshold;
and determining that the corresponding Xk value when the absolute value of the correction amount is smaller than the preset correction amount threshold is the final coordinate of the pixel point of the ground corridor part.
And calculating the ground inclination angle of the position where the tower is located according to the final coordinate.
The embodiments provided in the present application are only a few examples of the general concept of the present application, and do not limit the scope of the present application. Any other embodiments extended according to the scheme of the present application without inventive efforts will be within the scope of protection of the present application for a person skilled in the art.

Claims (9)

1. The method for detecting the ground inclination angle of the transmission tower by utilizing the SAR-X high-resolution bunching mode is characterized by comprising the following steps of:
randomly selecting a two-dimensional SAR image of a power transmission corridor at a certain time phase;
obtaining a main image intensity image by interfering the two-dimensional SAR image of the power transmission corridor;
dividing a ground corridor part on the main image intensity image according to a division threshold value;
calculating the backscattering coefficient sum of the ground corridor part;
extracting a pole tower ground plane part in the main image intensity graph according to the backscattering coefficient sum, and outputting a ground corridor pixel point coordinate on the main image intensity graph;
converting the coordinates of the pixel points of the ground corridor to the coordinates of a geocentric system;
calculating the final coordinate of the ground corridor pixel point according to the geocentric system coordinate;
and calculating the ground inclination angle of the position where the tower is located according to the final coordinates of the ground corridor pixel points.
2. The method for detecting the ground inclination angle of the transmission tower by utilizing the SAR-X high-resolution bunching mode as claimed in claim 1, wherein the two-dimensional SAR image of the transmission corridor takes the coordinates of the tower as the center and the root length of the tower which is more than 5 times as the radius.
3. The method for detecting the ground inclination angle of the transmission tower by using the SAR-X high-resolution bunching mode as claimed in claim 1, wherein the step of determining the segmentation threshold comprises the steps of:
calculating the intensity average value of each point of the main image intensity map;
marking the point with the intensity smaller than the average value in the main image intensity image as a noise point;
all the noise points are removed to obtain calculation points;
and calculating the intensity average value of all the calculation points, and determining the intensity average value as a segmentation threshold value.
4. The method for detecting the ground inclination angle of the transmission tower by using the SAR-X high resolution bunching mode as claimed in claim 1, wherein the specific step of dividing the ground corridor part on the main image intensity map according to the division threshold value is as follows:
dividing the main image intensity map according to the segmentation threshold;
and determining the part of the main image intensity image with the intensity larger than the segmentation threshold value as a ground corridor part, and determining the part of the main image intensity image with the intensity smaller than or equal to the segmentation threshold value as a background part.
5. The method for detecting the ground inclination angle of the transmission tower by using the SAR-X high resolution bunching mode as claimed in claim 1, wherein the step of calculating the backscattering coefficient sum of the ground corridor part comprises the following steps:
obtaining structural parameters of the ground corridor part: the center point coordinates and the radius length of the tower, wherein the radius length is 5 times of the tower root opening length at minimum;
determining the backscattering coefficient of all circular circumferences of the ground corridor part with integral multiple of the pole tower root opening length as the radius to the maximum circular circumference according to the structural parameters;
the sum of the backscatter coefficients for all the terrestrial corridor sections is calculated.
6. The method for detecting the ground inclination angle of the transmission tower by using the SAR-X high resolution bunching mode as claimed in claim 1, wherein the specific step of outputting the coordinates of the pixel points of the ground corridor on the main image intensity map comprises:
marking the pixel points of the ground corridor part as a first color, and outputting a marking point coordinate on the main image intensity graph;
traversing the coordinates of the mark points according to a preset template;
determining points, the matching degree of which is in accordance with a preset matching threshold value, of the preset template as fine extraction points;
marking the fine extraction point as a second color, and outputting the fine extraction point coordinates on the main image intensity map.
7. The method for detecting the ground inclination angle of the transmission tower by using the SAR-X high resolution bunching mode as claimed in claim 1, wherein the step of converting the coordinates of the pixel points of the ground corridor to the coordinates of the geocentric system further comprises:
and fitting a satellite orbit equation according to the 12 position vectors and the velocity vector of the required satellite in the geocentric coordinate system:
Figure FDA0002461289100000021
Figure FDA0002461289100000022
Figure FDA0002461289100000023
wherein Xsf, Ysf and XZf are satellite position vectors, Vxsf, VYsf and VZsf are satellite velocity vectors, i is the number of rows of pixel points on an image, PRF is pulse repetition frequency, t is time interval, a0、a1、a2、a3、b1、b2、b3、c1、c2、c3A set of fixed unknown parameters that need to be solved for.
8. The method for detecting the ground inclination angle of the transmission tower by utilizing the SAR-X high-resolution bunching mode as claimed in claim 1, wherein the specific step of establishing the functional relationship between the slant range image coordinate and the earth center system coordinate is as follows;
determining a distance equation according to:
F1=(Xsf-Xgf)2+(Ysf-Ygf)2+(Zsf-Zgf)2-R2
wherein Xgf, Ygf and Zgf are position vectors of ground points in a geocentric coordinate system, and R is a slope distance;
determining the doppler equation according to:
Figure FDA0002461289100000024
wherein V is the wave velocity of incident wave, lambda is the wavelength of incident wave, fDCIs the Doppler frequency shift;
determining an earth ellipsoid equation according to:
Figure FDA0002461289100000031
wherein N is the semiaxis of the ellipse of the earth, N (1-e)2) Is the short semi-axis of the earth ellipse, and h is the elevation of the ground point.
9. The method for detecting the ground inclination angle of the transmission tower by using the SAR-X high resolution bunching mode as claimed in claim 1, wherein the position of the pixel point of the ground corridor part in the geocentric coordinate system is calculated according to the distance equation, the Doppler equation and the earth ellipsoid equation, and the following formula is obtained:
Δx=B-1·F(Xk)
wherein, Deltax is the correction amount,
Figure FDA0002461289100000032
b is F (X)k) The value obtained by deriving Xgf, Ygf, Zgf is the B value according to the following formula:
Figure FDA0002461289100000033
calculating an absolute value of the correction amount;
if the absolute value of the correction amount is larger than or equal to the preset correction amount threshold, carrying out correction calculation again until the absolute value of the correction amount is smaller than the preset correction amount threshold;
determining X corresponding to the absolute value of the correction amount smaller than the preset correction amount thresholdkAnd the value is the final coordinate of the pixel point of the ground corridor part.
CN202010320737.7A 2020-04-22 2020-04-22 Method for detecting ground inclination angle of transmission tower by utilizing SAR-X high-resolution bunching mode Pending CN111487619A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010320737.7A CN111487619A (en) 2020-04-22 2020-04-22 Method for detecting ground inclination angle of transmission tower by utilizing SAR-X high-resolution bunching mode

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010320737.7A CN111487619A (en) 2020-04-22 2020-04-22 Method for detecting ground inclination angle of transmission tower by utilizing SAR-X high-resolution bunching mode

Publications (1)

Publication Number Publication Date
CN111487619A true CN111487619A (en) 2020-08-04

Family

ID=71813006

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010320737.7A Pending CN111487619A (en) 2020-04-22 2020-04-22 Method for detecting ground inclination angle of transmission tower by utilizing SAR-X high-resolution bunching mode

Country Status (1)

Country Link
CN (1) CN111487619A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109085588A (en) * 2018-09-26 2018-12-25 云南电网有限责任公司电力科学研究院 Based on the inclined monitoring method of Terra SAR-X high-resolution beam bunching mode data power grid iron tower
CN112684449A (en) * 2021-03-22 2021-04-20 北京东方至远科技股份有限公司 Water area power line sag inversion method and device based on SAR technology

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102072992A (en) * 2010-10-26 2011-05-25 国网电力科学研究院武汉南瑞有限责任公司 Detailed terrain data-based transmission line failure-shielding and lightning-protection performance evaluation method
CN102955160A (en) * 2011-08-19 2013-03-06 湖北省电力公司电力科学研究院 Three-dimensional laser radar technology based transmission line tower parameter determination method
CN106557600A (en) * 2015-09-30 2017-04-05 中国电力科学研究院 A kind of power transmission tower unequal leg and the emulation mode on basis
CN107680102A (en) * 2017-08-28 2018-02-09 国网甘肃省电力公司电力科学研究院 A kind of airborne cloud data electric force pole tower extraction method based on space constraint
CN109085588A (en) * 2018-09-26 2018-12-25 云南电网有限责任公司电力科学研究院 Based on the inclined monitoring method of Terra SAR-X high-resolution beam bunching mode data power grid iron tower
CN110441767A (en) * 2019-09-06 2019-11-12 云南电网有限责任公司电力科学研究院 The measurement method and system of power transmission line corridor screen of trees head room
CN110595433A (en) * 2019-08-16 2019-12-20 太原理工大学 Binocular vision-based transmission tower inclination measurement method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102072992A (en) * 2010-10-26 2011-05-25 国网电力科学研究院武汉南瑞有限责任公司 Detailed terrain data-based transmission line failure-shielding and lightning-protection performance evaluation method
CN102955160A (en) * 2011-08-19 2013-03-06 湖北省电力公司电力科学研究院 Three-dimensional laser radar technology based transmission line tower parameter determination method
CN106557600A (en) * 2015-09-30 2017-04-05 中国电力科学研究院 A kind of power transmission tower unequal leg and the emulation mode on basis
CN107680102A (en) * 2017-08-28 2018-02-09 国网甘肃省电力公司电力科学研究院 A kind of airborne cloud data electric force pole tower extraction method based on space constraint
CN109085588A (en) * 2018-09-26 2018-12-25 云南电网有限责任公司电力科学研究院 Based on the inclined monitoring method of Terra SAR-X high-resolution beam bunching mode data power grid iron tower
CN110595433A (en) * 2019-08-16 2019-12-20 太原理工大学 Binocular vision-based transmission tower inclination measurement method
CN110441767A (en) * 2019-09-06 2019-11-12 云南电网有限责任公司电力科学研究院 The measurement method and system of power transmission line corridor screen of trees head room

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
许浩文: "基于分层结构的铁塔散射模型研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》, pages 042 - 2488 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109085588A (en) * 2018-09-26 2018-12-25 云南电网有限责任公司电力科学研究院 Based on the inclined monitoring method of Terra SAR-X high-resolution beam bunching mode data power grid iron tower
CN109085588B (en) * 2018-09-26 2022-06-07 云南电网有限责任公司电力科学研究院 Method for monitoring power grid iron tower inclination based on Terra SAR-X high-resolution bunching mode data
CN112684449A (en) * 2021-03-22 2021-04-20 北京东方至远科技股份有限公司 Water area power line sag inversion method and device based on SAR technology
CN112684449B (en) * 2021-03-22 2021-07-16 北京东方至远科技股份有限公司 Water area power line sag inversion method and device based on SAR technology

Similar Documents

Publication Publication Date Title
Bühler et al. Snow depth mapping in high-alpine catchments using digital photogrammetry
Kubota et al. Verification of high-resolution satellite-based rainfall estimates around Japan using a gauge-calibrated ground-radar dataset
Cummins et al. An overview of lightning locating systems: History, techniques, and data uses, with an in-depth look at the US NLDN
CN110044327B (en) Infrastructure settlement monitoring method and system based on SAR data and GNSS data
LÓPEZ-MORENO et al. Ground-based remote-sensing techniques for diagnosis of the current state and recent evolution of the Monte Perdido Glacier, Spanish Pyrenees
KR101880616B1 (en) Method of sea fog prediction based on Sea surface winds and sea fog information from satellite
CN110488151B (en) Remote sensing technology-based power transmission line vegetation hazard early warning system and method
CN109085588B (en) Method for monitoring power grid iron tower inclination based on Terra SAR-X high-resolution bunching mode data
CN111487619A (en) Method for detecting ground inclination angle of transmission tower by utilizing SAR-X high-resolution bunching mode
Alpers et al. Observation of a nonlinear wave disturbance in the marine atmosphere by the synthetic aperture radar aboard the ERS 1 satellite
Kaasalainen et al. Brief communication" Application of mobile laser scanning in snow cover profiling"
Luce et al. Validation of winds measured by MU radar with GPS radiosondes during the MUTSI campaign
CN112198490B (en) Ground clutter suppression method and device for airborne millimeter wave cloud detection radar and radar
Nag et al. Lightning locating systems: characteristics and validation techniques
CN111540007B (en) Method for estimating oil tank storage capacity by SAR image
Lin et al. An ERS-1 synthetic aperture radar image of a tropical squall line compared with weather radar data
Gu et al. Comparative analysis of microwave brightness temperature data in Northeast China using AMSR-E and MWRI products
Bühler et al. Spatially continuous mapping of snow depth in high alpine catchments using digital photogrammetry
Rouveure et al. Radar-based SLAM without odometric sensor
KR20170104100A (en) System and Method For Ground Clutter Removing of WindProfiler
Rahardjo et al. Estimating additional flood warning time using raindrops intensity data captured by X-Band MP radar
Wang et al. Comparison of Doppler-Derived Sea Ice Radial Surface Velocity Measurement Methods From Sentinel-1A IW Data
OHKI et al. Flood Area Detection Using ALOS-2 PALSAR-2 Data for the 2015 Heavy Rainfall Disaster in the Kanto and Tohoku Area, Japan
KL Continental-scale detection of cloud-to-ground lightning
CN117406778B (en) Unmanned plane laser radar ground-imitating flight method based on geospatial data

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination