Disclosure of Invention
Aiming at the problems, the invention provides a high-voltage line detection method of a high-resolution millimeter wave radar image so as to improve the low-altitude flight and safe obstacle avoidance capability of a helicopter/large and medium unmanned aerial vehicle.
The invention aims to provide a high-voltage line detection method of a millimeter wave radar image, which comprises the following steps:
s1: performing Constant False Alarm Rate (CFAR) detection on radar echoes in one antenna scanning period to obtain a target measuring point;
s2: performing iterative agglomeration on the target point;
s3: eliminating clutter area data;
s4: dividing a high-voltage tower point and a high-voltage line point;
s5: converting the coordinates of the condensation points into coordinates in a three-dimensional rectangular coordinate system with the radar as an origin;
s6: extracting a straight line as a candidate high-voltage line based on Hough transform (Hough transform); a
S7: setting a characteristic segment of a candidate high-voltage wire;
s8: extracting a plurality of suspected high-voltage wires according to Bragg effect (Bragg effect) characteristics;
s9: associating a plurality of suspected high-voltage wires to identify the high-voltage wires;
s10: and combining the high-voltage wires at the overlapped part to obtain a final high-voltage wire detection result.
The high-voltage line detection method of the millimeter wave radar image provided by the invention is further characterized in that iterative agglomeration is carried out according to the threshold values of the distance and the orientation in the step S2 so as to remove the false alarm in the target measuring point obtained in the step S1.
The method for detecting the high-voltage line of the millimeter wave radar image further has the characteristic that clutter region data are removed according to the radial distance and the azimuth threshold of the clutter region in the step S3.
The method for detecting the high-voltage line of the millimeter wave radar image, provided by the invention, is further characterized in that the step S4 comprises the following steps:
s4.1: setting a radial distance threshold range and an azimuth threshold range of a high-voltage tower point, wherein the radial distance threshold of the high-voltage tower point is smaller than a radial distance threshold of a clutter area, and the azimuth threshold of the high-voltage tower point is smaller than the azimuth threshold of the clutter area;
s4.2: and selecting candidate high-voltage tower points from the residual data after the clutter area data are removed in the S3, wherein the distance length and the direction bit length of the candidate high-voltage tower points are both in the range set by the S4.1, and the rest points are candidate high-voltage line points.
The method for detecting the high-voltage line of the millimeter wave radar image further has the characteristic that in S7, the incidence angle of the current radar wave is recorded as alpha, and the incidence angle range is selected as alpha in the candidate high-voltage linedown≤α≤αupThe straight line segment of the candidate high-voltage line is taken as a characteristic segment of the candidate high-voltage line.
The method for detecting the high-voltage line of the millimeter wave radar image, provided by the invention, is further characterized in that the step S8 comprises the following steps:
s8.1: extracting an angle-power relation curve of a candidate high-voltage line characteristic section from radar echoes, and searching a maximum power peak value angle and marking as beta based on an incidence angle alpha of a current radar wave;
s8.2: over a range of incident angles [ alpha ]downBeta) searching the maximum peak point in turn, and recording the incidence angle as alphad1With a power of Pd1,Δα1=β-αd1
S8.3: setting the incidence angle range of the symmetrical peak point search as delta and a power judgment threshold PthresFrom the range of incident angles [ beta + delta alpha ]1-δ/2,β+Δα1+δ/2]The maximum peak point is searched, wherein,
β+Δα1+δ/2≤αupif the maximum peak point is searched, the incidence angle is recorded as alphau1And its power Pu1Satisfy constraint | Pu1-Pd1|≤PthresCalculating a recognition confidence coefficient rho:
s8.4: if rho is greater than 0.8, the candidate high-voltage line is considered to have the Bragg effect, and the candidate high-voltage line is identified as a suspected high-voltage line; otherwise, the next step is carried out:
s8.5: in the range [ alpha ]downBeta) searching the peak point with the second highest power, and recording the corresponding incidence angle as alphad2With a power of Pd2,Δα2=β-αd2At its corresponding incident angle range [ beta + delta alpha ]2-δ/2,β+Δα2+δ/2]Searching for the second largest peak point, where β + Δ α2+δ/2≤αupAnd judging whether the peak point power meets the constraint condition | Pu2-Pd2|≤PthresIf yes, calculating confidence coefficient, otherwise, analogizing, and calculating the confidence coefficient at [ alpha ]downBeta) at most before co-completion of NpeakSearching for the maximum peak point.
The method for detecting the high-voltage line of the millimeter wave radar image, provided by the invention, is further characterized in that the step S9 comprises the following steps:
s9.1: judging whether the slopes of the suspected high-voltage lines are the same or not, performing S9.2 on the same suspected high-voltage line data, and performing S9.3 on different suspected high-voltage line data:
s9.2: a calculation of the confidence of the association recognition is performed,
wherein n is the number of suspected high-voltage lines with the same slope, and S9.4 is carried out after calculation;
s9.3: the associated identification confidence coefficient corresponding to the suspected high-voltage line with the slope different from the data of the rest suspected high-voltage lines is zetai=ρiThen S9.4 is carried out;
s9.4: setting a threshold, respectively judging the associated identification confidence degrees of the suspected high-voltage wires acquired in S9.2 and S9.3, if the associated identification confidence degrees are not greater than the threshold, judging the suspected high-voltage wires as false high-voltage wires, otherwise, judging the suspected high-voltage wires as real high-voltage wires.
The method for detecting a high-voltage line of a millimeter wave radar image provided by the invention is further characterized in that the step S10 comprises the following steps:
s10.1: setting the radial distance threshold and the azimuth threshold of the high-voltage wire combination as R respectivelytower_thresAnd thetatower_thres;
S10.2: selecting an end point of a certain high-voltage line, and calculating the absolute values of the radial distance difference and the azimuth difference between the end point and a line point and a tower point on the other high-voltage line to be delta RtowerAnd Δ θtowerWhen Δ R istower<Rtower_thresAnd Δ θtower<θtower_thresIf so, judging that the end point falls on the other high-voltage line, judging whether the other end point on the high-voltage line has a coincidence point on the other high-voltage line, and if so, judging that the high-voltage line is coincident with the other high-voltage line;
s10.3: if the superposed high-voltage wires exist, merging is carried out, the strip containing the largest number of the tower wire points is reserved, and if the superposed high-voltage wires do not exist, data are directly output.
Compared with the prior art, the invention has the beneficial effects that:
the high-voltage wire detection method of the millimeter wave radar image provided by the invention aims at the high-voltage wire identification problem of the high-resolution millimeter wave radar image, comprehensively utilizes the special Bragg effect characteristics of the high-voltage wire, reduces the degree of dependence on the high-voltage tower detection, has high-voltage wire extraction accuracy, can detect and identify a plurality of high-voltage wires in the distribution of the same group of high-voltage towers, and can be applied to the field of millimeter wave radar collision avoidance.
Detailed Description
In order to make the technical means, the creation features, the achievement purposes and the effects of the invention easy to understand, the following embodiments are specifically described in the detection method provided by the invention with reference to the drawings.
In the description of the embodiments of the present invention, it should be understood that the terms "central", "longitudinal", "lateral", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", etc. indicate orientations or positional relationships based on those shown in the drawings, and are only used for convenience in describing and simplifying the description of the present invention, but do not indicate or imply that the referred device or element must have a specific orientation, be constructed and operated in a specific orientation, and thus, should not be construed as limiting the present invention.
Furthermore, the terms "first," "second," "third," and the like are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicit to a number of indicated technical features. Thus, a feature defined as "first," "second," etc. may explicitly or implicitly include one or more of that feature. In the description of the invention, the meaning of "a plurality" is two or more unless otherwise specified.
The terms "mounted," "connected," and "coupled" are to be construed broadly and may, for example, be fixedly coupled, detachably coupled, or integrally coupled; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meaning of the above terms in the creation of the present invention can be understood by those of ordinary skill in the art through specific situations.
As shown in fig. 1, a method for detecting a high-voltage line of a millimeter wave radar image includes the following steps:
s1: performing Constant False Alarm Rate (CFAR) detection on radar echoes in one antenna scanning period to obtain a target measuring point;
s2: iterative agglomeration of target points:
setting neighborhood radial distance threshold RthresAnd an orientation threshold thetathresA radial distance threshold RthresAnd an orientation threshold thetathresThe target detection points in the neighborhood are condensed into a target point, and the radial distance and the azimuth value of the target point are condensedFocusing the mean value of the radial distance and the azimuth of the target point, and similarly, performing iterative condensation;
s3: eliminating clutter area data:
setting radial distance threshold of clutter area
And azimuth threshold
The elimination distance in the condensation point is longer than
Or an azimuthal length greater than
A clutter zone of;
s4: dividing a high-voltage tower point and a high-voltage line point;
s5: converting the radial distance, azimuth angle and pitch angle of the condensation point into a three-dimensional rectangular coordinate system with radar as the origin, wherein the height direction is the Z axis, and projecting the Z axis into an X-Y rectangular coordinate system
S6: extracting straight lines as candidate high-voltage lines based on Hough transformation (Hough transformation) in the condensation point X-Y rectangular coordinate system distribution diagram obtained at S5;
s7: recording the incidence angle of the current radar wave as alpha, and selecting the incidence angle range as alpha in the candidate high-voltage linedown≤α≤αupThe straight line segment of the candidate high-voltage line is taken as a characteristic segment of the candidate high-voltage line;
s8: extracting a plurality of suspected high-voltage wires according to Bragg effect (Bragg effect) characteristics;
s9: associating a plurality of suspected high-voltage wires to identify the high-voltage wires;
s10: and combining the high-voltage wires at the overlapped part to obtain a final high-voltage wire detection result.
In some embodiments, iterative clustering is performed according to the distance and orientation thresholds in S2 to remove false alarms in the target station obtained in S1.
In some embodiments, clutter region data is removed based on the clutter region radial distance and the bearing threshold in S3.
In some embodiments, the S4 includes the following steps:
s4.1: setting a radial distance threshold range and an azimuth threshold range of a high-voltage tower point, wherein the radial distance threshold of the high-voltage tower point is smaller than a radial distance threshold of a clutter area, and the azimuth threshold of the high-voltage tower point is smaller than the azimuth threshold of the clutter area;
s4.2: and selecting candidate high-voltage tower points from the residual data after the clutter area data are removed in the S3, wherein the distance length and the direction bit length of the candidate high-voltage tower points are both in the range set by the S4.1, and the rest points are candidate high-voltage line points.
In some embodiments, in S7, the incident angle of the current radar wave is recorded as α, and the incident angle range is selected as α in the candidate high-voltage linedown≤α≤αupThe straight line segment of the candidate high-voltage line is taken as a characteristic segment of the candidate high-voltage line.
In some embodiments, the S8 includes the following steps:
s8.1: extracting an angle-power relation curve of a candidate high-voltage line characteristic section from radar echoes, and searching a maximum power peak value angle and marking as beta based on an incidence angle alpha of a current radar wave;
s8.2: over a range of incident angles [ alpha ]downBeta) searching the maximum peak point in turn, and recording the incidence angle as alphad1With a power of Pd1,Δα1=β-αd1
S8.3: setting the incidence angle range of the symmetrical peak point search as delta and a power judgment threshold PthresFrom the range of incident angles [ beta + delta alpha ]1-δ/2,β+Δα1+δ/2]The maximum peak point is searched, wherein,
β+Δα1+δ/2≤αupif the maximum peak point is searched, the incidence angle is recorded as alphau1And its power Pu1Satisfy constraint | Pu1-Pd1|≤PthresCalculating a recognition confidence coefficient rho:
s8.4: if rho is greater than 0.8, the candidate high-voltage line is considered to have the Bragg effect, and the candidate high-voltage line is identified as a suspected high-voltage line; otherwise, the next step is carried out:
s8.5: in the range [ alpha ]downBeta) searching the peak point with the second highest power, and recording the corresponding incidence angle as alphad2With a power of Pd2,Δα2=β-αd2At its corresponding incident angle range [ beta + delta alpha ]2-δ/2,β+Δα2+δ/2]Searching for the second largest peak point, where β + Δ α2+δ/2≤αupAnd judging whether the peak point power meets the constraint condition | Pu2-Pd2|≤PthresIf yes, calculating confidence coefficient, otherwise, analogizing, and calculating the confidence coefficient at [ alpha ]downBeta) at most before co-completion of NpeakSearching for the maximum peak point.
In some embodiments, the S9 includes the following steps:
s9.1: judging whether the slopes of the suspected high-voltage lines are the same or not, performing S9.2 on the same suspected high-voltage line data, and performing S9.3 on different suspected high-voltage line data:
s9.2: a calculation of the confidence of the association recognition is performed,
wherein n is the number of suspected high-voltage lines with the same slope, and S9.4 is carried out after calculation;
s9.3: the associated identification confidence coefficient corresponding to the suspected high-voltage line with the slope different from the data of the rest suspected high-voltage lines is zetai=ρiThen S9.4 is carried out;
s9.4: setting a threshold, respectively judging the associated identification confidence degrees of the suspected high-voltage wires acquired in S9.2 and S9.3, if the associated identification confidence degrees are not greater than the threshold, judging the suspected high-voltage wires as false high-voltage wires, otherwise, judging the suspected high-voltage wires as real high-voltage wires.
In some embodiments, the S10 includes:
s10.1: setting the radial distance threshold and the azimuth threshold of the high-voltage wire combination as R respectivelytower_thresAnd thetatower_thres;
S10.2: selecting an end point of a certain high-voltage line, and calculating the absolute values of the radial distance difference and the azimuth difference between the end point and a line point and a tower point on the other high-voltage line to be delta RtowerAnd Δ θtowerWhen Δ R istower<Rtower_thresAnd Δ θtower<θtower_thresIf so, judging that the end point falls on the other high-voltage line, judging whether the other end point on the high-voltage line has a coincidence point on the other high-voltage line, and if so, judging that the high-voltage line is coincident with the other high-voltage line;
s10.3: if the superposed high-voltage wires exist, merging is carried out, the strip containing the largest number of the tower wire points is reserved, and if the superposed high-voltage wires do not exist, data are directly output.
In some embodiments, the test is performed using high voltage line echo data collected by certain W-band FMCW millimeter wave radar radiation with a range resolution of 0.15 meters, an antenna scan cycle millimeter wave radar echo pattern is shown in FIG. 2, where the antenna scans ranging from-50 to 50 in azimuth and 100 to 600 in radial azimuth, FIG. 2(b) is an enlarged view of the incident angle range [ -525 ] and radial distance range [ 450510 ] of FIG. 2(a), and the optical image of the high voltage tower/line in the radiation scene is shown in FIG. 3,
step 1: performing CFAR detection on radar echoes in one antenna scanning period to obtain target detection points, wherein the length of a protection window half window of the CFAR detection is 5, the length of a reference window half window is 20, a detection threshold is 80 decibels, and fig. 4 is a target detection point distribution diagram of a corresponding area in fig. 2;
step 2: taking neighborhood distance threshold Rthres0.5 and an orientation threshold θthresCarrying out iterative aggregation on the candidate target detection points as 0.4;
and step 3: the radial distance of the clutter zone is set to the threshold
And an azimuth threshold of
Removing impuritiesWave zone processing;
and 4, step 4: respectively selecting the distance direction threshold of the condensation high-pressure tower point as Rt′hres1=1.5,Rt′hres210 and an orientation threshold θt′hres1=0.5,θt′hres2Determining candidate high-voltage towers and high-voltage line points as 2;
and 5: converting the condensation points into a three-dimensional rectangular coordinate system with the radar as the origin, and projecting the three-dimensional rectangular coordinate system into an X-Y rectangular coordinate system, wherein the distribution of the XY rectangular coordinate system of the regional condensation points shown in FIG. 4 is shown in FIG. 5;
step 6: extracting a straight line as a candidate high-voltage line based on Hough transformation in the rectangular coordinate system distribution diagram of the condensation points X-Y obtained by the processing of thestep 5;
and 7: take alpha respectivelydownAnd alphaupThe feature segments of the 4 selected candidate high-voltage lines are shown in fig. 7 and 30 degrees, and are respectively represented by different lines, the "+" signs of the corresponding lines respectively represent candidate high-voltage tower points and high-voltage line points on the candidate high-voltage lines, and fig. 6 is a schematic diagram of the candidate high-voltage lines extracted in fig. 5(b), and the schematic diagram is sequentially named as a first candidate high-voltage line, a second candidate high-voltage line, a third candidate high-voltage line and a fourth candidate high-voltage line from top to bottom.
And 8: the distribution of the Bragg scattering points of the characteristic segments of the 4 candidate high-voltage wires is shown in FIG. 7, wherein N is takenpeak=5,δ=1,PthresTaking a third candidate high-voltage line characteristic segment as an example, where the maximum peak point corresponds to an incident angle β of 11.4 °, and the maximum peak points are searched within a range of [ -5 °,11.4 °, and the incident angles are α, respectivelyd12.5 DEG, power Pd1135.1 db,. DELTA.alpha1At [19.8 °,20.8 ° ], 8.9 °]Searching the corresponding peak point with the incidence angle alphau120.1 deg. and power Pd1The power difference is 3.1 db, 132 db, and the confidence in the identification of the third candidate high voltage line is calculated to be 0.85. Similarly, the confidence of the identification of the second, first and fourth candidate high voltage lines is 0.89,0.86 and 0.82, respectively.
And step 9: calculating the associated recognition confidence coefficients of the 4 high-voltage line groups with the same slope as follows: 0.97,1,0.98 and 0.94, the group of suspected high voltage lines are all identified as high voltage lines.
Step 10: selecting the distance difference threshold of the merging of the high-voltage lines as Rtower_thresThe azimuth difference threshold is theta 3tower_thresAnd (5) combining and outputting high-voltage wires which possibly have an overlapping part or are close to each other in tower points. Fig. 8 is a result of the detection, identification and enhanced display of the high voltage lines in fig. 2(b), wherein the high voltage lines are all represented by distinct lines.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and is not to be construed as limiting the invention, and any modifications, equivalents and improvements made within the spirit and principle of the present invention are intended to be included within the scope of the present invention. The above description is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, several modifications and variations can be made without departing from the technical principle of the present invention, and these modifications and variations should also be regarded as the protection scope of the present invention.