CN105045274A - Intelligent tower connected graph construction method for unmanned aerial vehicle inspection track planning - Google Patents

Intelligent tower connected graph construction method for unmanned aerial vehicle inspection track planning Download PDF

Info

Publication number
CN105045274A
CN105045274A CN201510218379.8A CN201510218379A CN105045274A CN 105045274 A CN105045274 A CN 105045274A CN 201510218379 A CN201510218379 A CN 201510218379A CN 105045274 A CN105045274 A CN 105045274A
Authority
CN
China
Prior art keywords
tower
msub
mrow
line
sigma
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
CN201510218379.8A
Other languages
Chinese (zh)
Other versions
CN105045274B (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.)
China South Power Grid International Co ltd
Tianjin Aerospace Zhongwei Date Systems Technology Co Ltd
Original Assignee
China South Power Grid International Co ltd
Tianjin Aerospace Zhongwei Date Systems Technology 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 China South Power Grid International Co ltd, Tianjin Aerospace Zhongwei Date Systems Technology Co Ltd filed Critical China South Power Grid International Co ltd
Priority to CN201510218379.8A priority Critical patent/CN105045274B/en
Publication of CN105045274A publication Critical patent/CN105045274A/en
Application granted granted Critical
Publication of CN105045274B publication Critical patent/CN105045274B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides a method for constructing an intelligent tower connectivity graph for planning an unmanned aerial vehicle inspection track, which comprises the following steps of: A. predicting the distribution of the power transmission line towers by adopting a unitary nonlinear regression prediction method to generate a plurality of tower lines, wherein each tower line covers a plurality of towers; B. connecting the plurality of tower lines to form a line connection diagram by using the critical condition of tower distribution; wherein the critical conditions include: the method comprises the following steps of (1) crossing conditions, parallel distribution conditions of a plurality of lines at close distances, tower steering conditions and tower branch conditions; C. and constructing a storage structure of the line connection graph. The invention can utilize the algorithm to carry out intelligent planning on the power transmission line, thereby improving the planning efficiency; the invention provides a whole-course inspection scheme, which not only inspects the route, but also inspects the route, thereby greatly improving the efficiency of line inspection; according to the invention, through an intelligent traversal algorithm, the optimal patrol route at the planned position in a line network is met, and traversal of all towers in the area under the condition of the shortest patrol distance is met.

Description

Intelligent tower connected graph construction method for unmanned aerial vehicle inspection track planning
Technical Field
The invention relates to a method for constructing an intelligent tower connectivity graph for planning an unmanned aerial vehicle inspection track.
Background
The inspection of the overhead line is the basic work of a power grid company on the routine maintenance of the power transmission line, and the inspection modes can be divided into manual inspection, inspection by a manned helicopter and inspection by an unmanned aerial vehicle. Although manual inspection is the most commonly used inspection method, the defects of slow efficiency, restriction by climatic geographic environment and the like exist all the time, the mechanical inspection method is widely researched and applied, especially the safety and high efficiency characteristic of unmanned aerial vehicle inspection, and an overhead power line inspection Flying Robot (FROPI) has greater and greater application value.
The unmanned aerial vehicle system comprises an unmanned aerial vehicle body platform, a task load and a data wireless transmission part. The existing unmanned aerial vehicle platform applied to power transmission line inspection mainly comprises a micro multi-rotor unmanned aerial vehicle, an unmanned aerial vehicle helicopter and the like; the task load mainly comprises the charge of polling imaging equipment such as a visible light camera, a thermal infrared imager and the like; the wireless transmission mainly comprises the transmission (data transmission) of remote control instructions and the transmission (image transmission) of high-bandwidth image data, which are collectively called as a measurement and control data chain.
The unmanned aerial vehicle system is a main body of power transmission line inspection work, an inspection area is firstly determined in the autonomous line inspection process of the unmanned aerial vehicle, the flight line of the unmanned aerial vehicle is planned according to a shooting angle and a flight safety distance on the basis of obtaining accurate line geographic coordinates, and the unmanned aerial vehicle is sent by a starting/landing point to inspect and fly the line according to the flight line uploaded in advance. The unmanned aerial vehicle acquires ideal patrol data as much as possible in the patrol process, and ensures absolute safety and reliability of flight, the flight area and flight route of the unmanned aerial vehicle directly influence patrol quality, patrol efficiency and even the safety of the unmanned aerial vehicle flight and overhead lines, so that the mission planning of the unmanned aerial vehicle for the power transmission line plays a crucial role. Meanwhile, overhead transmission lines are distributed in mountainous areas in many ways, cross-over areas exist, and transmission lines in a large area are complicated, so that the selection of the optimal routing inspection line in a transmission line network in a large area plays a key role in improving the operation efficiency of the unmanned aerial vehicle, improving the energy efficiency ratio and flying safely.
At present, a manual planning mode is mostly adopted for task planning of power transmission line unmanned aerial vehicle inspection, the mode ensures the flight safety of the unmanned aerial vehicle, but the efficiency is low, the requirement of large-scale unmanned aerial vehicle inspection cannot be met, and meanwhile, the manual planning mode is difficult to realize optimal planning in a large area.
The intelligent unmanned aerial vehicle line patrol task planning firstly needs to construct a data structure for all tower positions in the whole overhead line area so as to carry out intelligent algorithm planning.
Chinese patent application 201410283175.8 discloses a high-precision three-dimensional reconstruction method for power transmission lines and corridors. The method comprises the steps of collecting three-dimensional point cloud data of a power transmission line and a corridor through an airborne laser radar, detecting the three-dimensional point cloud data, eliminating wrong points and points with abnormal elevations, automatically classifying the three-dimensional point cloud data by using a voxel classification based method, and finally automatically reconstructing the power transmission line and the corridor through an algorithm.
Chinese patent application 201310676886.7 discloses an unmanned aerial vehicle flight path planning algorithm based on Dubins path and sparse a search. Combining the Dubins path with the sparse A search algorithm, adopting the Dubins path length as a heuristic function of the sparse A search algorithm, and searching nodes in a space by using the heuristic function to realize the flight path planning of the unmanned aerial vehicle.
The defects of the technical scheme are as follows: 1. the existing autonomous task planning means of the overhead transmission line mainly depends on a manual planning mode, the efficiency is low, and the requirement of patrol of large-area complex lines cannot be met; 2. the flight path planning is generally only one-way task planning, namely, the inspection is only carried out on a one-way route, and the task inspection is not carried out in the return process, so that the energy consumption is wasted; 3. in a power transmission line network with cross spanning, the optimal inspection line is difficult to find out from the overall consideration by a manual planning method, so that the flight distance is shortest on the premise of finishing inspection of all lines.
Disclosure of Invention
Aiming at the defects of the prior art, the invention aims to provide a method for constructing a connection diagram of an intelligent tower for planning an inspection flight path of an unmanned aerial vehicle.
In order to achieve the purpose, the invention provides a method for constructing an intelligent tower connectivity graph for unmanned aerial vehicle inspection track planning, which comprises the following steps:
A. predicting the distribution of the power transmission line towers by adopting a unitary nonlinear regression prediction method to generate a plurality of tower lines, wherein each tower line covers a plurality of towers;
in the step, firstly, a three-dimensional space model is simplified into a two-dimensional space model, and the distribution of the power transmission line towers is predicted by adopting a unitary nonlinear regression prediction method; and the prediction mode adopts a confidence interval mode, and new input data (explanatory variables) are predicted and judged according to historical tower data (response variables). Due to the particularity of tower distribution, the special case types are summarized, and the existence condition of each special case is determined so as to be classified in the algorithm processing process.
B. Connecting a plurality of tower lines to form a line connection diagram by using the critical condition of tower distribution; wherein the critical conditions include: the method comprises the following steps of (1) crossing conditions, parallel distribution conditions of a plurality of lines at close distances, tower steering conditions and tower branch conditions;
in the step, if the pole tower distribution reaches a critical condition due to the occurrence of special conditions, critical type judgment is carried out, and processing is carried out; the critical types of tower distribution are divided into the following types:
cross-over situation: predicting that a plurality of points are generated in the interval, the number of the points exceeds the set number of the towers, and intersection points exist;
the parallel distribution condition that the distance of a plurality of lines is short: predicting that a plurality of points appear in the interval, exceed the set number of towers and have no intersection points;
the tower turns to the condition: compared with a prediction equation, a unique inflection point appears in a prediction interval;
the tower branch situation: in contrast to the prediction equation, a plurality of inflection points occur within the prediction interval.
C. And constructing a storage structure of the line connection diagram.
According to the method, the tower connected graph is intelligently constructed under the condition that only the geographical coordinates of the tower exist through a linear regression prediction algorithm, so that intelligent task planning is facilitated.
According to another embodiment of the present invention, step a specifically includes the following steps:
a1, dimension reduction treatment: carrying out dimension reduction processing on the geographic coordinates of the three-dimensional tower to convert the geographic coordinates into two-dimensional coordinates; setting the three-dimensional coordinate of the original A pole tower as (x)t,yt,zt),xtRepresenting the three-dimensional spatial longitude coordinate, y, of the towertRepresenting tower latitude coordinate, ztThe altitude of the tower is shown, and the coordinates of the tower A after the dimensionality reduction are (x)t,yt);
A2, establishment of regression equation: the unitary linear regression prediction model formula applied to the power transmission line mission planning is as follows:
Y ^ t = a + b x t - - - ( 1 )
in the formula xtThe longitude coordinates of the tower at the moment t are shown,representing the estimated latitude coordinate at the time t;
a3, retrieving the prediction step size as N, and obtaining the solving equation of the parameters a and b in the regression equation as follows:
<math><mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mi>a</mi> <mo>=</mo> <mfrac> <mrow> <mi>&Sigma;</mi> <msub> <mi>Y</mi> <mi>i</mi> </msub> </mrow> <mi>N</mi> </mfrac> <mo>-</mo> <mi>b</mi> <mfrac> <mrow> <mi>&Sigma;</mi> <msub> <mi>X</mi> <mi>i</mi> </msub> </mrow> <mi>N</mi> </mfrac> </mtd> </mtr> <mtr> <mtd> <mi>b</mi> <mo>=</mo> <mfrac> <mrow> <mi>N&Sigma;</mi> <msub> <mi>Y</mi> <mi>i</mi> </msub> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <mi>&Sigma;</mi> <msub> <mi>Y</mi> <mi>i</mi> </msub> <mi>&Sigma;</mi> <msub> <mi>X</mi> <mi>i</mi> </msub> </mrow> <mrow> <mi>N&Sigma;</mi> <msup> <msub> <mi>X</mi> <mi>i</mi> </msub> <mn>2</mn> </msup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mi>&Sigma;</mi> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow></math>
wherein,n is a predicted moving step length; because the distance between two base pole towers is different from dozens of meters to hundreds of meters, most of the cases are that the continuous multiple base poles and towers are composed and can approximate to a straight line segment.
According to another embodiment of the present invention, step N of step A3 is 5 m to 10 m.
According to another embodiment of the present invention, step a further comprises step a 4: and (4) performing optimal prediction fitting, and determining a curve determined by the tower coordinates through the optimal fitting curve. Since the curve can be approximated to a straight line within a certain distance except for special cases, a unary linear regression method is adopted in the invention. The curve fitting adopts a 'least square method', namely 'residual square sum minimum' to determine the position of a straight line, and in terms of formula (1), the estimation method of the 'least square method-OLS' should meet the condition that the smaller the difference ei between an actual observed value Yi and a predicted value, the better, the ei represents the formula as follows:
<math><mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mo>&Sigma;</mo> <msup> <msub> <mi>e</mi> <mi>i</mi> </msub> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mo>=</mo> <mi>m</mi> <mi>i</mi> <mi>n</mi> <msup> <mrow> <mo>(</mo> <msub> <mi>Y</mi> <mi>i</mi> </msub> <mo>-</mo> <mover> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mo>^</mo> </mover> <mo>-</mo> <mover> <msub> <mi>&beta;</mi> <mn>2</mn> </msub> <mo>^</mo> </mover> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow></math>
solving by using the claime rule to obtain an OLS estimation formula in the form of an observed value as follows:
<math><mrow> <mfenced open = '{' close = ''> <mtable> <mtr> <mtd> <mrow> <mover> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mo>^</mo> </mover> <mo>=</mo> <mfrac> <mrow> <mo>&Sigma;</mo> <msup> <msub> <mi>X</mi> <mi>i</mi> </msub> <mn>2</mn> </msup> <mo>&Sigma;</mo> <msub> <mi>Y</mi> <mi>i</mi> </msub> <mo>-</mo> <mo>&Sigma;</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>&Sigma;</mo> <msub> <mi>Y</mi> <mi>i</mi> </msub> <msub> <mi>X</mi> <mi>i</mi> </msub> </mrow> <mrow> <mi>N</mi> <mo>&Sigma;</mo> <msup> <msub> <mi>X</mi> <mi>i</mi> </msub> <mn>2</mn> </msup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mrow> <mo>&Sigma;</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mover> <msub> <mi>&beta;</mi> <mn>2</mn> </msub> <mo>^</mo> </mover> <mo>=</mo> <mfrac> <mrow> <mi>N</mi> <mo>&Sigma;</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <msub> <mi>Y</mi> <mi>i</mi> </msub> <mo>-</mo> <mo>&Sigma;</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>&Sigma;</mo> <msub> <mi>Y</mi> <mi>i</mi> </msub> </mrow> <mrow> <mi>N</mi> <mo>&Sigma;</mo> <msup> <msub> <mi>X</mi> <mi>i</mi> </msub> <mn>2</mn> </msup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mrow> <mo>&Sigma;</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow></math>
fitting a final linear equation according to the processing result as follows:
<math><mrow> <mover> <mi>Y</mi> <mo>&OverBar;</mo> </mover> <mo>=</mo> <mover> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mo>^</mo> </mover> <mo>+</mo> <mover> <msub> <mi>&beta;</mi> <mn>2</mn> </msub> <mo>^</mo> </mover> <mover> <mi>X</mi> <mo>&OverBar;</mo> </mover> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow></math>
according to another embodiment of the present invention, step a further comprises step a 5: and constructing a prediction interval. Because the curve equation determined by the actual tower and the prediction equation have a certain deviation, the Y value is subjected to interval prediction, namely a prediction interval of the average value is constructed. And according to the distribution condition of the towers, setting a significance level a, and calculating a prediction interval with the confidence coefficient of the Y average value being 1-a.
According to another embodiment of the present invention, step C specifically includes the following steps:
c1, establishing a tower matrix, wherein the row and column coordinates represent the tower number, and the matrix data are tower geographical position information; the towers are divided into two types according to the prediction result: a node tower and a non-node tower; the nodes comprise a line start-stop tower and a cross tower; the non-node tower is an internal tower only belonging to a single line;
c2, constructing a node tower adjacency list for the node tower;
and C3, for the non-node tower, storing a matrix structure, wherein the matrix row coordinates represent the line to which the non-node tower belongs, and the column coordinates represent the tower number.
In step C2, for the node tower, a chain storage structure-adjacency list is used in consideration of the algorithm storage space and the algorithm efficiency due to the large tower data size.
The construction process of the adjacency list for the node tower is as follows:
c21, sorting according to the geographic position coordinates of the towers and the reference direction;
c22, the data structure of each tower vertex has the following attributes:
a. subscripts: the method is used for marking the tower line;
b. front adjacent point: when the node tower is a cross tower, the previous adjacent line of the line where the front adjacent point is located;
c. rear adjacent point: when the node tower is a cross tower, the next adjacent line of the line where the rear adjacent point is located is adjacent to the line.
In step C3, matrix structure storage is performed on the non-node tower, the matrix row coordinates represent the line to which the non-node tower belongs, and the column coordinates represent the tower number. Note the following:
I. whether the node is a head node, wherein the head node represents one of the end points of one line;
II. Whether it is a cross node: the intersection point value is '1', and the common node is '0';
III, the circuit of the device: the common node has only one line, and the cross node has a plurality of lines.
Compared with the prior art, the invention has the following beneficial effects:
1. the invention can utilize the algorithm to carry out intelligent planning on the power transmission line, thereby improving the planning efficiency;
2. the invention provides a whole-course inspection scheme, which not only inspects the route, but also inspects the route, thereby greatly improving the efficiency of line inspection;
3. aiming at the task planning of a complex power transmission line network, the optimal patrol route is planned in the line network through an intelligent traversal algorithm, and the traversal of all towers in an area under the condition of the shortest patrol distance is met.
The present invention will be described in further detail with reference to the accompanying drawings.
Drawings
Fig. 1 is an overall flowchart of the intelligent tower connectivity graph construction method of embodiment 1;
FIG. 2 is a flow chart of step A in example 1.
Detailed Description
Example 1
The embodiment provides a method for constructing an intelligent tower connectivity graph for planning an unmanned aerial vehicle inspection track, which comprises the following steps of:
A. and predicting the distribution of the power transmission line towers by adopting a unitary nonlinear regression prediction method to generate a plurality of tower lines, wherein each tower line covers a plurality of towers. As shown in fig. 2, step a specifically includes the following steps:
a1, dimension reduction treatment: carrying out dimension reduction processing on the geographic coordinates of the three-dimensional tower to convert the geographic coordinates into two-dimensional coordinates; setting the three-dimensional coordinate of the original A pole tower as (x)t,yt,zt),xtRepresenting the three-dimensional spatial longitude coordinate, y, of the towertRepresenting tower latitude coordinate, ztIndicating the altitude of the towerAnd the coordinates of the A pole tower after dimension reduction are (x)t,yt) (ii) a But the storage structure still keeps the tower elevation information for planning the flight path.
A2, establishment of regression equation: the unitary linear regression prediction model formula applied to the power transmission line mission planning is as follows:
Y ^ t = a + b x t - - - ( 1 )
in the formula xtThe longitude coordinates of the tower at the moment t are shown,representing the estimated latitude coordinate at the time t;
a3, retrieving the prediction step size as N, and obtaining the solving equation of the parameters a and b in the regression equation as follows:
<math><mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mi>a</mi> <mo>=</mo> <mfrac> <mrow> <mi>&Sigma;</mi> <msub> <mi>Y</mi> <mi>i</mi> </msub> </mrow> <mi>N</mi> </mfrac> <mo>-</mo> <mi>b</mi> <mfrac> <mrow> <mi>&Sigma;</mi> <msub> <mi>X</mi> <mi>i</mi> </msub> </mrow> <mi>N</mi> </mfrac> </mtd> </mtr> <mtr> <mtd> <mi>b</mi> <mo>=</mo> <mfrac> <mrow> <mi>N&Sigma;</mi> <msub> <mi>Y</mi> <mi>i</mi> </msub> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <mi>&Sigma;</mi> <msub> <mi>Y</mi> <mi>i</mi> </msub> <mi>&Sigma;</mi> <msub> <mi>X</mi> <mi>i</mi> </msub> </mrow> <mrow> <mi>N&Sigma;</mi> <msup> <msub> <mi>X</mi> <mi>i</mi> </msub> <mn>2</mn> </msup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mi>&Sigma;</mi> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow></math>
wherein,n is a predicted moving step length; because the distance between two base pole towers is different from dozens of meters to hundreds of meters, most of the cases are that the continuous multiple base poles and towers are composed and can approximate to a straight line segment. According to the tower distribution condition of a specific area and the accuracy and rapidity of prediction are guaranteed, the step length N can be set to be 5-10 m.
And A4, best prediction fitting, and determining a curve determined by the tower coordinates through the best fitting curve. Since the curve can be approximated to a straight line within a certain distance except for special cases, a unary linear regression method is adopted in the invention. The curve fitting adopts a 'least square method', namely 'residual square sum minimum' to determine the position of a straight line, and in terms of formula (1), the estimation method of the 'least square method-OLS' should meet the condition that the smaller the difference ei between an actual observed value Yi and a predicted value, the better, the ei represents the formula as follows:
<math><mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <msup> <msub> <mi>&Sigma;e</mi> <mi>i</mi> </msub> <mn>2</mn> </msup> <mo>)</mo> </mrow> <mo>=</mo> <mi>m</mi> <mi>i</mi> <mi>n</mi> <msup> <mrow> <mo>(</mo> <msub> <mi>Y</mi> <mi>i</mi> </msub> <mo>-</mo> <mover> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mo>^</mo> </mover> <mo>-</mo> <mover> <msub> <mi>&beta;</mi> <mn>2</mn> </msub> <mo>^</mo> </mover> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow></math>
solving by using the claime rule to obtain an OLS estimation formula in the form of an observed value as follows:
<math><mrow> <mfenced open = '{' close = ''> <mtable> <mtr> <mtd> <mrow> <mover> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mo>^</mo> </mover> <mo>=</mo> <mfrac> <mrow> <mo>&Sigma;</mo> <msup> <msub> <mi>X</mi> <mi>i</mi> </msub> <mn>2</mn> </msup> <mo>&Sigma;</mo> <msub> <mi>Y</mi> <mi>i</mi> </msub> <mo>-</mo> <mo>&Sigma;</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>&Sigma;</mo> <msub> <mi>Y</mi> <mi>i</mi> </msub> <msub> <mi>X</mi> <mi>i</mi> </msub> </mrow> <mrow> <mi>N</mi> <mo>&Sigma;</mo> <msup> <msub> <mi>X</mi> <mi>i</mi> </msub> <mn>2</mn> </msup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mrow> <mo>&Sigma;</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mover> <msub> <mi>&beta;</mi> <mn>2</mn> </msub> <mo>^</mo> </mover> <mo>=</mo> <mfrac> <mrow> <mi>N</mi> <mo>&Sigma;</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <msub> <mi>Y</mi> <mi>i</mi> </msub> <mo>-</mo> <mo>&Sigma;</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>&Sigma;</mo> <msub> <mi>Y</mi> <mi>i</mi> </msub> </mrow> <mrow> <mi>N</mi> <mo>&Sigma;</mo> <msup> <msub> <mi>X</mi> <mi>i</mi> </msub> <mn>2</mn> </msup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mrow> <mo>&Sigma;</mo> <msub> <mi>X</mi> <mi>i</mi> </msub> </mrow> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow></math>
fitting a final linear equation according to the processing result as follows:
<math><mrow> <mover> <mi>Y</mi> <mo>&OverBar;</mo> </mover> <mo>=</mo> <mover> <msub> <mi>&beta;</mi> <mn>1</mn> </msub> <mo>^</mo> </mover> <mo>+</mo> <mover> <msub> <mi>&beta;</mi> <mn>2</mn> </msub> <mo>^</mo> </mover> <mover> <mi>X</mi> <mo>&OverBar;</mo> </mover> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow></math>
and A5, constructing a prediction interval. Because the curve equation determined by the actual tower and the prediction equation have a certain deviation, the Y value is subjected to interval prediction, namely a prediction interval of the average value is constructed. And according to the distribution condition of the towers, setting a significance level a, and calculating a prediction interval with the confidence coefficient of the Y average value being 1-a.
B. Connecting a plurality of tower lines to form a line connection diagram by using the critical condition of tower distribution; wherein the critical conditions include: the method comprises the following steps of (1) crossing conditions, parallel distribution conditions of a plurality of lines at close distances, tower steering conditions and tower branch conditions;
in the step, if the pole tower distribution reaches a critical condition due to the occurrence of special conditions, critical type judgment is carried out, and processing is carried out; the critical types of tower distribution are divided into the following types:
cross-over situation: predicting that a plurality of points are generated in the interval, the number of the points exceeds the set number of the towers, and intersection points exist;
the parallel distribution condition that the distance of a plurality of lines is short: predicting that a plurality of points appear in the interval, exceed the set number of towers and have no intersection points;
the tower turns to the condition: compared with a prediction equation, a unique inflection point appears in a prediction interval;
the tower branch situation: in contrast to the prediction equation, a plurality of inflection points occur within the prediction interval.
C. And constructing a storage structure of the line connection diagram. The step C specifically comprises the following steps:
c1, establishing a tower matrix, wherein the row and column coordinates represent the tower number, and the matrix data are tower geographical position information; the towers are divided into two types according to the prediction result: a node tower and a non-node tower; the nodes comprise a line start-stop tower and a cross tower; the non-node tower is an internal tower only belonging to a single line;
c2, constructing a node tower adjacency list for the node tower;
and C3, for the non-node tower, storing a matrix structure, wherein the matrix row coordinates represent the line to which the non-node tower belongs, and the column coordinates represent the tower number.
In step C2, for the node tower, a chain storage structure-adjacency list is used in consideration of the algorithm storage space and the algorithm efficiency due to the large tower data size.
The construction process of the adjacency list for the node tower is as follows:
c21, sorting according to the geographic position coordinates of the towers and the reference direction;
c22, the data structure of each tower vertex has the following attributes:
a. subscripts: the method is used for marking the tower line;
b. front adjacent point: when the node tower is a cross tower, the previous adjacent line of the line where the front adjacent point is located;
c. rear adjacent point: when the node tower is a cross tower, the next adjacent line of the line where the rear adjacent point is located is adjacent to the line.
In step C3, matrix structure storage is performed on the non-node tower, the matrix row coordinates represent the line to which the non-node tower belongs, and the column coordinates represent the tower number. Note the following:
I. whether the node is a head node, wherein the head node represents one of the end points of one line;
II. Whether it is a cross node: the intersection point value is '1', and the common node is '0';
III, the circuit of the device: the common node has only one line, and the cross node has a plurality of lines.
Although the present invention has been described with reference to the preferred embodiments, it is not intended to limit the scope of the invention. It will be appreciated by those skilled in the art that changes may be made without departing from the scope of the invention, and it is intended that all matter contained in the above description or shown in the accompanying drawings shall be interpreted as illustrative and not in a limiting sense.

Claims (5)

1. An intelligent tower connected graph construction method for planning an unmanned aerial vehicle inspection track is characterized by comprising the following steps:
A. predicting the distribution of the power transmission line towers by adopting a unitary nonlinear regression prediction method to generate a plurality of tower lines, wherein each tower line covers a plurality of towers;
B. connecting the plurality of tower lines to form a line connection diagram by using the critical condition of tower distribution; wherein the critical conditions include: the method comprises the following steps of (1) crossing conditions, parallel distribution conditions of a plurality of lines at close distances, tower steering conditions and tower branch conditions;
C. and constructing a storage structure of the line connection graph.
2. The connectivity graph construction method according to claim 1, wherein the step a specifically includes the steps of:
a1, dimension reduction treatment: carrying out dimension reduction processing on the geographic coordinates of the three-dimensional tower to convert the geographic coordinates into two-dimensional coordinates;
a2, establishment of regression equation: the unitary linear regression prediction model formula applied to the power transmission line mission planning is as follows:
Y ^ t = a + bx t
in the formula xtThe longitude coordinates of the tower at the moment t are shown,representing the estimated latitude coordinate at the time t;
a3, retrieving the prediction step size as N, and obtaining the solving equation of the parameters a and b in the regression equation as follows:
<math> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mi>a</mi> <mo>=</mo> <mfrac> <msub> <mi>&Sigma;Y</mi> <mi>i</mi> </msub> <mi>N</mi> </mfrac> <mo>-</mo> <mi>b</mi> <mfrac> <msub> <mi>&Sigma;X</mi> <mi>i</mi> </msub> <mi>N</mi> </mfrac> </mtd> </mtr> <mtr> <mtd> <mi>b</mi> <mo>=</mo> <mfrac> <mrow> <mi>N&Sigma;</mi> <msub> <mi>Y</mi> <mi>i</mi> </msub> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>-</mo> <mi>&Sigma;</mi> <msub> <mi>Y</mi> <mi>i</mi> </msub> <mi>&Sigma;</mi> <msub> <mi>X</mi> <mi>i</mi> </msub> </mrow> <mrow> <mi>N&Sigma;</mi> <msup> <msub> <mi>X</mi> <mi>i</mi> </msub> <mn>2</mn> </msup> <mo>-</mo> <msup> <mrow> <mo>(</mo> <mi>&Sigma;</mi> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> </mtd> </mtr> </mtable> </mfenced> </math>
wherein, <math> <mrow> <mi>&Sigma;</mi> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mo>,</mo> </mrow> </math> and N is the predicted moving step.
3. The method for constructing a connectivity graph according to claim 2, wherein in the step a3, the step length N is 5 m to 10 m.
4. The connectivity graph building method according to claim 2, wherein the step a further comprises a step a 4: and (4) performing optimal prediction fitting, and determining a curve determined by the tower coordinates through the optimal fitting curve.
5. The connectivity graph constructing method according to claim 1, wherein the step C specifically includes the steps of:
c1, establishing a tower matrix, wherein the row and column coordinates represent the tower number, and the matrix data are tower geographical position information; the towers are divided into two types according to the prediction result: a node tower and a non-node tower; the nodes comprise a line start-stop tower and a cross tower; the non-node tower is an internal tower only belonging to a single line;
c2, constructing a node tower adjacency list for the node tower;
and C3, for the non-node tower, storing a matrix structure, wherein the matrix row coordinates represent the line to which the non-node tower belongs, and the column coordinates represent the tower number.
CN201510218379.8A 2015-04-30 2015-04-30 Intelligent tower connected graph construction method for unmanned aerial vehicle inspection track planning Active CN105045274B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510218379.8A CN105045274B (en) 2015-04-30 2015-04-30 Intelligent tower connected graph construction method for unmanned aerial vehicle inspection track planning

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510218379.8A CN105045274B (en) 2015-04-30 2015-04-30 Intelligent tower connected graph construction method for unmanned aerial vehicle inspection track planning

Publications (2)

Publication Number Publication Date
CN105045274A true CN105045274A (en) 2015-11-11
CN105045274B CN105045274B (en) 2019-01-04

Family

ID=54451877

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510218379.8A Active CN105045274B (en) 2015-04-30 2015-04-30 Intelligent tower connected graph construction method for unmanned aerial vehicle inspection track planning

Country Status (1)

Country Link
CN (1) CN105045274B (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107527395A (en) * 2016-06-21 2017-12-29 中兴通讯股份有限公司 A kind of unmanned plane method for inspecting, device and unmanned plane
CN108364369A (en) * 2018-01-17 2018-08-03 上海振华重工(集团)股份有限公司 Unmanned plane inspection point determines method, apparatus, medium, equipment and system
CN110825110A (en) * 2019-11-13 2020-02-21 昆明能讯科技有限责任公司 Acquisition flight method for power line visible light point cloud resolving photo
CN111443730A (en) * 2020-05-12 2020-07-24 江苏方天电力技术有限公司 Unmanned aerial vehicle track automatic generation method and device for power transmission line inspection
CN111609855A (en) * 2019-12-27 2020-09-01 北京数字绿土科技有限公司 Method for generating refined routing inspection routes in batch based on tower shapes
CN112414409A (en) * 2020-11-16 2021-02-26 天津航天中为数据系统科技有限公司 Autonomous inspection path planning method based on string structure and aircraft
CN113009929A (en) * 2021-03-01 2021-06-22 中国南方电网有限责任公司超高压输电公司贵阳局 Unmanned aerial vehicle power grid line patrol control method
CN113052339A (en) * 2021-03-31 2021-06-29 广东电网有限责任公司 Unmanned aerial vehicle inspection method and device based on line crossing and storage medium
CN113110577A (en) * 2021-04-15 2021-07-13 中国南方电网有限责任公司超高压输电公司柳州局 Unmanned aerial vehicle flight route planning management system is patrolled and examined to electric wire netting
CN113910227A (en) * 2021-10-11 2022-01-11 国网智能科技股份有限公司 Distribution network vehicle-mounted intelligent inspection robot task-free inspection method and system
CN115686070A (en) * 2022-11-18 2023-02-03 国网江苏省电力有限公司苏州供电分公司 Unmanned aerial vehicle autonomous line patrol path planning method for power distribution/power transmission line

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102353945A (en) * 2011-03-31 2012-02-15 北京航空航天大学 Three-dimensional position reconstructing method based on ISAR (inverse synthetic aperture radar) image sequence for scattering point
US8311686B2 (en) * 2009-08-20 2012-11-13 The Boeing Company Onboard airplane community noise and emissions monitor
CN102819663A (en) * 2012-07-17 2012-12-12 哈尔滨工程大学 Method for forecasting ship wake based on optimized support vector regression parameter
CN102116625B (en) * 2009-12-31 2013-01-23 武汉大学 GIS (geographic information system)-GPS (global position system) navigation method of inspection robot
CN102968554A (en) * 2012-11-01 2013-03-13 武汉大学 Tower pole icing disaster risk prediction method based on safety margin
CN103020733A (en) * 2012-11-27 2013-04-03 南京航空航天大学 Method and system for predicting single flight noise of airport based on weight
CN103915810A (en) * 2014-04-17 2014-07-09 国家电网公司 Lightning protection method for 500 kV circuit JG type pole and tower
CN104049639A (en) * 2014-06-24 2014-09-17 上海大学 Unmanned surface vehicle anti-surge control device and method based on support vector regression

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8311686B2 (en) * 2009-08-20 2012-11-13 The Boeing Company Onboard airplane community noise and emissions monitor
CN102116625B (en) * 2009-12-31 2013-01-23 武汉大学 GIS (geographic information system)-GPS (global position system) navigation method of inspection robot
CN102353945A (en) * 2011-03-31 2012-02-15 北京航空航天大学 Three-dimensional position reconstructing method based on ISAR (inverse synthetic aperture radar) image sequence for scattering point
CN102819663A (en) * 2012-07-17 2012-12-12 哈尔滨工程大学 Method for forecasting ship wake based on optimized support vector regression parameter
CN102968554A (en) * 2012-11-01 2013-03-13 武汉大学 Tower pole icing disaster risk prediction method based on safety margin
CN103020733A (en) * 2012-11-27 2013-04-03 南京航空航天大学 Method and system for predicting single flight noise of airport based on weight
CN103915810A (en) * 2014-04-17 2014-07-09 国家电网公司 Lightning protection method for 500 kV circuit JG type pole and tower
CN104049639A (en) * 2014-06-24 2014-09-17 上海大学 Unmanned surface vehicle anti-surge control device and method based on support vector regression

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
于陈平等: "基于回归分析的航迹关联算法研究", 《舰船电子对抗》 *
郭志军: "应用Excel对一元线性回归模型的分析", 《宁波职业技术学院学报》 *

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107527395A (en) * 2016-06-21 2017-12-29 中兴通讯股份有限公司 A kind of unmanned plane method for inspecting, device and unmanned plane
CN108364369B (en) * 2018-01-17 2023-09-26 上海振华重工(集团)股份有限公司 Unmanned aerial vehicle inspection point determining method, unmanned aerial vehicle inspection point determining device, unmanned aerial vehicle inspection point determining medium, unmanned aerial vehicle inspection point determining equipment and unmanned aerial vehicle inspection point determining system
CN108364369A (en) * 2018-01-17 2018-08-03 上海振华重工(集团)股份有限公司 Unmanned plane inspection point determines method, apparatus, medium, equipment and system
CN110825110A (en) * 2019-11-13 2020-02-21 昆明能讯科技有限责任公司 Acquisition flight method for power line visible light point cloud resolving photo
CN111609855A (en) * 2019-12-27 2020-09-01 北京数字绿土科技有限公司 Method for generating refined routing inspection routes in batch based on tower shapes
CN111609855B (en) * 2019-12-27 2021-05-18 北京数字绿土科技有限公司 Method for generating refined routing inspection routes in batch based on tower shapes
CN111443730A (en) * 2020-05-12 2020-07-24 江苏方天电力技术有限公司 Unmanned aerial vehicle track automatic generation method and device for power transmission line inspection
CN111443730B (en) * 2020-05-12 2021-02-12 江苏方天电力技术有限公司 Unmanned aerial vehicle track automatic generation method and device for power transmission line inspection
CN112414409A (en) * 2020-11-16 2021-02-26 天津航天中为数据系统科技有限公司 Autonomous inspection path planning method based on string structure and aircraft
CN113009929A (en) * 2021-03-01 2021-06-22 中国南方电网有限责任公司超高压输电公司贵阳局 Unmanned aerial vehicle power grid line patrol control method
CN113052339A (en) * 2021-03-31 2021-06-29 广东电网有限责任公司 Unmanned aerial vehicle inspection method and device based on line crossing and storage medium
CN113110577A (en) * 2021-04-15 2021-07-13 中国南方电网有限责任公司超高压输电公司柳州局 Unmanned aerial vehicle flight route planning management system is patrolled and examined to electric wire netting
CN113910227A (en) * 2021-10-11 2022-01-11 国网智能科技股份有限公司 Distribution network vehicle-mounted intelligent inspection robot task-free inspection method and system
CN115686070A (en) * 2022-11-18 2023-02-03 国网江苏省电力有限公司苏州供电分公司 Unmanned aerial vehicle autonomous line patrol path planning method for power distribution/power transmission line
CN115686070B (en) * 2022-11-18 2024-01-23 国网江苏省电力有限公司苏州供电分公司 Autonomous line patrol path planning method for unmanned aerial vehicle of power distribution/transmission line

Also Published As

Publication number Publication date
CN105045274B (en) 2019-01-04

Similar Documents

Publication Publication Date Title
CN105045274B (en) Intelligent tower connected graph construction method for unmanned aerial vehicle inspection track planning
He et al. Research of multi-rotor UAVs detailed autonomous inspection technology of transmission lines based on route planning
CN106371456B (en) A kind of unmanned plane patrolling method and system
CN108109437B (en) Unmanned aerial vehicle autonomous route extraction and generation method based on map features
CN110244750B (en) Unmanned aerial vehicle tour path planning method and device
Salvo et al. Urban traffic analysis through an UAV
WO2017219780A1 (en) Tour inspection method and device for unmanned aerial vehicle, unmanned aerial vehicle, and computer storage medium
Li The UAV intelligent inspection of transmission lines
Lawrance et al. Path planning for autonomous soaring flight in dynamic wind fields
CN110888453A (en) Unmanned aerial vehicle autonomous flight method for constructing three-dimensional real scene based on LiDAR data
CN113342046A (en) Power transmission line unmanned aerial vehicle routing inspection path optimization method based on ant colony algorithm
CN114153224B (en) Unmanned aerial vehicle flight path intelligent planning control system and method based on microclimate
JP2018146350A (en) Flight control system of unmanned aircraft and schedule route setting method
CN112684806A (en) Electric power inspection unmanned aerial vehicle system based on dual obstacle avoidance and intelligent identification
CN110687919A (en) Unmanned aerial vehicle autonomous flight control method for hydropower station water pipeline inspection
CN112215416B (en) Intelligent planning inspection route system and method
CN115275870A (en) Inspection system based on high-altitude line maintenance
CN116483118A (en) Unmanned aerial vehicle autonomous inspection route automatic generation method suitable for ground lead sag change
Xu et al. Iterative construction of low-altitude UAV air route network in urban areas: Case planning and assessment
CN115935610A (en) Method and system for optimizing routing inspection strategy of unmanned aerial vehicle on overhead line
CN115185303A (en) Unmanned aerial vehicle patrol path planning method for national parks and natural protected areas
Chehri et al. Accelerating power grid monitoring with flying robots and artificial intelligence
Yang et al. Multi-base multi-UAV cooperative patrol route planning novel method
CN109445462B (en) Unmanned aerial vehicle robust path planning method in uncertain environment
Battulwar et al. High-resolution modeling of open-pit slopes using UAV and photogrammetry

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant