CN114994633A - Method for secondary radar point track correlation processing - Google Patents
Method for secondary radar point track correlation processing Download PDFInfo
- Publication number
- CN114994633A CN114994633A CN202210930747.1A CN202210930747A CN114994633A CN 114994633 A CN114994633 A CN 114994633A CN 202210930747 A CN202210930747 A CN 202210930747A CN 114994633 A CN114994633 A CN 114994633A
- Authority
- CN
- China
- Prior art keywords
- track
- point
- correlation
- score
- matching
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details 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
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a method for secondary radar point track correlation processing, which comprises the following steps: s1, receiving trace point data according to the sector and changing the coordinate system of the position information; s2, traversing with effective track nodes in the track sector, calculating the identity matching marks of the current track and the current track nodes, and adaptively adjusting the relevant wave gate threshold; s3, calculating the position difference between the current detection point track and the current track node, and determining a point track correlation score according to the position difference and a correlation gate threshold; s4, counting the point track correlation scores of the effective track nodes of the current track sector, and finding the track node with the maximum point track correlation score as the optimal point track matching target, so as to facilitate the subsequent track updating processing.
Description
Technical Field
The invention relates to the technical field of secondary radar point track processing, in particular to a method for secondary radar point track correlation processing.
Background
The secondary radar data processing system consists of a point track processing part and a track processing part. And the trace point processing is to finish the aggregation of multiple responses to the target into trace point data, and comprises the steps of aggregating multiple response positions of a single target, extracting identity codes, height information, information amplitude and other response information, calculating the azimuth quality and the response times and eliminating basic false data. And the track processing is to complete point track correlation, track establishment, track filtering and prediction and track management. The point track correlation is to realize the correct interconnection between the point track data and the track nodes. And the track establishment is to establish a corresponding target track node according to the target track data. And in the track filtering and predicting, the optimal estimation method is adopted to process the track position data by adopting a filtering algorithm, and the estimated value is used for updating the track position and predicting the position which is most likely to appear in the next period of the target, so that the influence of observation noise is reduced, and more accurate track position information is obtained. And the track management is to extrapolate or clear overtime track nodes to ensure the timeliness of the track nodes. The function of the point track related processing is known, the point track related processing is a very critical link in the track processing, plays a very important role in the track processing, and is very helpful to the performance of secondary radar equipment.
The traditional secondary radar point track correlation processing method adopts a position nearest neighbor method and has the same identity code, and when targets are sparse in a clean airspace, the method is small in calculation amount, low in operation complexity and good in target correlation effect. With the increase of air targets, finer and more standardized navigation management area division and complicated geographic and electromagnetic environment influence, the following disadvantages can occur in the traditional point track related processing.
(1) When the air space is crossed, the identity code of the target is switched, and the phenomena of track fracture and batch change can occur due to unsuccessful track point correlation.
(2) When multiple targets fly at a short distance, target related errors occur, so that errors occur in subsequent track updating.
(3) Because the traditional method only has 1 gate number aiming at the target correlation and the gate threshold is fixed, under the interference of external factors, the detection effect of the point track direction is deteriorated, and the phenomena of heating track extrapolation, even track fracture, track batch change and track loss can occur due to unsuccessful point track correlation.
Disclosure of Invention
Aiming at the defects in the prior art, the method for the secondary radar point track correlation processing solves the problems that the identity code switching of the target occurs when the airspace is crossed, and the track is broken and batched; error correlation occurs during multi-target close-range flight, and then the track filtering effect is poor, or track extrapolation or track fracture batch change is performed; under the condition that the azimuth detection is not ideal, a large amount of track extrapolation problems occur due to unsuccessful point track correlation.
In order to achieve the purpose of the invention, the invention adopts the technical scheme that: a method for secondary radar point track correlation processing, comprising the steps of:
s1, receiving trace data according to the sector and changing the coordinate system of the position information;
s2, traversing with effective track nodes in the track sector, calculating the identity matching marks of the current track and the current track nodes, and adaptively adjusting the relevant wave gate threshold;
s3, calculating the position difference between the current detection point track and the current track node, and determining a point track correlation score according to the position difference and a correlation gate threshold;
s4, counting the point track correlation scores of the effective track nodes of the current track sector, and finding the track node with the maximum point track correlation score as the optimal point track matching target, so as to facilitate the subsequent track updating processing.
Further: the specific steps of step S1 are:
s11, receiving trace point data sent by the trace point processing system according to the sector;
and S12, transforming the space coordinate system of the received trace point data.
Further, the method comprises the following steps: and (3) finishing the elimination of basic false targets in a trace point processing system, wherein trace point data without the same identity code within the distance difference of 10km with the same sector are carried with the calculated trace point azimuth quality after trace point processing, and counting whether other targets exist within the range of 3km and 10km of the trace point.
Further: the specific steps of step S2 are:
s21, traversing with the effective track node in the corresponding track sector after receiving the trace point data, calculating the time difference, exiting the traversal if the time difference does not meet the track requirement, continuing to calculate whether the identity code carried by the trace point data is matched with the identity code of the current track node if the time difference meets the requirement, setting the identity matching identifier to be 1 if the matching is successful, and clearing the identity matching identifier to be 0 if the matching is not successful;
s22, aiming at the immature track and the mature track, self-adaptively adjusting a related wave gate threshold according to the identity matching identifier and the point track azimuth quality; and aiming at the extrapolated tracks, adaptively adjusting the threshold of the related wave gate according to the extrapolation times, the identity matching identifier and the point track azimuth quality, and dividing the wave gate into a small wave gate, a medium wave gate, a large wave gate and an extra-large wave gate through the threshold of the related wave gate.
Further: the specific steps of step S3 are:
s31, calculating a distance difference, an azimuth difference, an amplitude difference and a height difference between the track point position and the historical position or the predicted position of the track node; a height difference exists on the premise that the height value is effective;
s32, matching and comparing the calculated distance difference and the direction difference with a small wave gate, a middle wave gate, a big wave gate and a big wave gate in sequence, if the matching is successful, setting different point track relevant basic scores and setting a wave gate type identifier, if the identity matching identifier is 1 and the point track relevant basic Score is greater than zero, adding an identity code auxiliary Score code to the point track relevant basic Score to obtain a point track relevant Score, otherwise, entering the step S33;
s33, performing signal amplitude correlation matching on the premise that the point track correlation Score is larger than zero, adding the point track correlation Score to the signal amplitude auxiliary Score Score1 amplitude if the matching result is within 1dB, adding the point track correlation Score to the Score2 amplitude if the matching result is within 3dB, and otherwise, entering the step S34;
s34, on the premise that the point track correlation Score is nonzero and the height value is effective, if the height difference value is within the height layer interval of civil aviation planning, accumulating the point track correlation Score and the height of a height correlation auxiliary Score, otherwise, entering the step S35;
s35, if the identity matching identifier is 1, completing the matching scoring mechanism, and exiting the traversal; if the identity matching identifier is 0, continuously judging whether other target point track data exist in the trace point data 10km or not, if not, exiting the traversal, if the trace point data exist in the trace point data 10km, carrying out identity code consistency matching on the current track node and other trace point data which are not the current trace point in the sector, if identity codes are consistent, resetting the score related to the current point track, if the identity codes are different, exiting the traversal, and continuously executing the next circulation processing.
Further: the specific steps of step S4 are:
s41, when the track points and the track nodes in the sector are traversed, counting the track point related scores of the track nodes in the sector, and finding the track node with the maximum track point related score as the best track point matching result;
s42, clearing the point track related scores of other track nodes in the track sector, setting the wave gate type as invalid, only keeping the track node with the maximum point track related score, and giving the position index of the best matching track node.
The invention has the beneficial effects that: the invention uses the position as the basis for correlation, and uses the identity code, the signal amplitude and the height as the auxiliary correlation, thereby greatly improving the probability of successful target correlation, so that the identity code switching of the target occurs when the airspace is crossed, and the phenomena of track breakage and track batch number switching of the target can not occur.
The method is based on the correlation based on the position, takes the identity code, the signal amplitude and the height as auxiliary correlation, and simultaneously, under the premise that multi-target correlation exists, not only is the flight path matched based on the point path, but also the point path is reversely matched by the flight path node, so that the error correlation is greatly reduced, and the accuracy of the point flight path correlation is improved.
According to the track node attribute, the track point azimuth quality, the identity code and the extrapolation frequency, the size of the relevant wave gate is adjusted in a self-adaptive mode, the relevant correct probability is increased to the maximum extent, the relevant wave gate indication is assigned, the phenomena of track extrapolation and track restart batch are greatly reduced, and in addition, the relevant wave gate indication can also provide reference for a track filtering coefficient.
With the application of the invention on the interrogator equipment, the detection observation effect of the secondary radar is improved, the performance of the equipment is improved, the high recognition of the user is obtained, the workload of the operator on duty is reduced, and the powerful guarantee is provided for ensuring the flight safety of the airplane.
Drawings
FIG. 1 is a flow chart of a secondary radar point track correlation process of the present invention;
FIG. 2 is a flow chart of the design related to the spot track of the immature track in step S3;
FIG. 3 is a flow chart of the point track correlation design of the matured track in step S3 of the present invention;
FIG. 4 is a flow chart of the design related to the point track of the extrapolated track in step S3;
FIG. 5 is a comparison of a number of statistics using the method of the present invention.
Detailed Description
The following description of the embodiments of the present invention is provided to facilitate the understanding of the present invention by those skilled in the art, but it should be understood that the present invention is not limited to the scope of the embodiments, and it will be apparent to those skilled in the art that various changes may be made without departing from the spirit and scope of the invention as defined and defined in the appended claims, and all matters produced by the invention using the inventive concept are protected.
As shown in fig. 1, a method for secondary radar point track correlation processing includes the following steps:
s1, receiving trace data according to the sector and changing the coordinate system of the position information; the method comprises the following specific steps:
step 11: receiving trace point data sent by a trace point processing system according to a sector, wherein the trace point processing system finishes the elimination of a basic false target, and the trace point data without the same identity code has the same distance difference of 10km with the sector; after the trace point azimuth processing, the calculated trace point azimuth quality is carried; counting trace data of other targets in the range of 3km and 10km of the trace;
step 12: because the secondary radar detection is carried out under a polar coordinate system, a rectangular coordinate system is often adopted in the track processing, and the received point track data is transformed into a space coordinate system: assuming that the coordinate values in the polar coordinate system are (r, Ɵ, β) and the coordinate values in the rectangular coordinate system are (x, y, z), the formula is converted:
s2, traversing with effective track nodes in the track sector, calculating the identity matching marks of the current track and the current track nodes, and self-adaptively adjusting the relevant wave gate threshold; the specific process is as follows:
step 21: after receiving the trace point data, traversing and matching the trace point data with effective track nodes in corresponding track sectors, calculating time difference, and if the time difference does not meet the track processing requirement, exiting the traversal; if the time difference meets the requirement, whether the identity code carried by the trace point data is matched with the identity code of the current track node or not is continuously calculated, a mark is set to be 1 if the matching is successful, and the mark is clear to be 0 if the matching is unsuccessful;
step 22: aiming at the immature track and the mature track, when the identity matching identifier is 0, directly reducing the threshold of the relevant wave gate by half; when the identity matching identifier is 1 and the trace point quality is high, keeping the current relevant gate threshold; when the identity matching identifier is 1 and the trace point quality is medium quality, the amplification of the direction threshold related to the current relevant wave gate is kept by 1.5 times, and the distance wave gate is kept unchanged; when the identity matching identifier is 1 and the trace point quality is low, amplifying the direction threshold related to the current relevant wave gate by 2 times, and keeping the distance wave gate unchanged;
step 23: for the extrapolated track, if the identity matching identifier is 0, halving the threshold of the relevant wave gate; if the identity matching identifier is 0, amplifying the relevant gate threshold by (1.0+0.25n) times according to the extrapolation times (n), and if the direction quality of the trace point is high quality, keeping the current relevant gate threshold; if the point trace azimuth quality is the medium quality, amplifying the azimuth-related wave gate threshold by 1.25 times, and keeping the distance from the azimuth-related wave gate unchanged; and if the azimuth quality of the point trace is low, amplifying the threshold of the azimuth-related wave gate by 1.5 times, and keeping the distance-related wave gate unchanged.
S3, calculating the position difference between the current detection point track and the current track node, and determining a point track correlation score according to the position difference and a correlation gate threshold; the method comprises the following specific steps:
step 31: calculating a distance difference value, an azimuth difference value, an amplitude difference value and a height difference value (on the premise that the height is effective) between the track point position and the historical position or the predicted position of the track node;
step 32: matching and comparing the distance difference and the azimuth difference with a wavelet gate, a middle wave gate, a big wave gate and a super-big wave gate in sequence according to the calculated distance difference and azimuth difference, and if the matching is successful, setting different point flight path related basic scores (Score in Score is small, Score is large, and Score is super-big) and setting the types of the related wave gates; if the identity code mark is 1 and the point track correlation Score is greater than zero, adding the point track correlation Score and the identity code auxiliary Score code, and if not, continuing to execute the next step;
step 33: on the premise that the point track correlation Score is larger than zero, performing correlation matching of signal amplitude, if the point track correlation Score is within 1dB, adding the signal amplitude to an auxiliary Score Score1 amplitude, if the point track correlation Score is within 3dB, adding a Score2 amplitude (the Score1 amplitude is larger than the Score2 amplitude), and if the point track correlation Score is not larger than zero, continuing the next step;
step 34: on the premise that the point track correlation Score is not zero and the height value is effective, if the height difference is within the height layer interval of civil aviation planning, accumulating the point track correlation Score and the height of the height correlation auxiliary Score, and if not, continuing the next step;
step 35: if the identity matching identifier is 1, completing the matching scoring mechanism, and exiting the traversal; if the identity matching identifier is 0, continuously judging whether other target trace data exist in the trace data 10km or not, and if not, exiting the traversal; if other targets exist within 10km, identity code consistency matching is carried out on the current track node and other trace point data which are not current trace point in the sector, and if identity codes are consistent, the current point track related fraction is reset; if the identity codes are different, exiting the traversal, and continuing to execute the next loop processing.
Point track correlation processing of immature tracks, as shown in fig. 2:
in order to reduce the false rate of the target, the initial conditions of the flight path are stricter, and when the identity codes are inconsistent, the loop exits; when the identity codes are consistent, continuously judging whether the direction and the distance of the detection point are in a wavelet gate (taking the updated position of the previous circle as the center), if so, setting the point track correlation score as s1 and setting the correlation type as the wavelet gate; if the detected point is not in the wavelet gate, continuously judging whether the azimuth and the distance of the detected point are in the middle wave gate (the updated position of the previous circle is taken as the center), if so, setting the point track related fraction as s2 and setting the type of the related wave gate as the middle wave gate; if the detected point is not in the medium wave gate, continuously judging whether the azimuth and the distance of the detected point are in the large wave gate (the updated position of the previous circle is the center), if so, setting the point track related fraction as s3(s1> s2> s3) and setting the type of the related wave gate as the large wave gate; if the target position is not in the big wave gate, otherwise, the traversal is regarded as that the target position is not related, and the traversal can be quitted;
under the condition that the point track related fraction is larger than 0, judging whether the amplitude difference is in the range of 0 dB-1 dB, and if so, adding the point track related fraction to an auxiliary score s 4; if the amplitude difference is within the range of 1 dB-3 dB, if so, adding the auxiliary score s5 to the point track related score; if the amplitude difference is larger than 3dB, the point track related fraction is not processed, and the next step is continuously executed;
and under the condition that the point track related score is larger than 0 and the height is effective, judging whether the altitude difference is in the altitude layer interval range, if so, adding the auxiliary score s6 to the point track related score, otherwise, not processing the point track related score, and exiting the traversal.
Point track correlation processing of mature track, as shown in fig. 3:
when the identity codes are inconsistent, automatically reducing the threshold value of the relevant wave gate by half; then, under the same condition when the identity codes are consistent, continuously judging whether the azimuth and the distance of the detection point are in a wavelet gate (the predicted position of the previous circle is taken as the center), if the azimuth and the distance are in the wavelet gate and the identity codes are consistent, setting the point track correlation score as y1 and setting the correlation type as the wavelet gate; if the point flight path correlation score is within the wavelet gate and the identity codes are inconsistent, setting the point flight path correlation score as y2 and setting the correlation type as the wavelet gate; if not, continuously judging whether the azimuth and the distance of the detection point are in a medium wave gate (the predicted position of the previous circle is taken as the center); if the navigation path is in the medium wave gate and the identity codes are consistent, setting the point track correlation score as y3 and setting the correlation wave gate type as the medium wave gate; if the point flight path is in the medium wave gate and the identity codes are not consistent, setting the point flight path related score as y4 and setting the related wave gate type as the medium wave gate; if the detected point is not in the medium wave gate, continuously judging whether the azimuth and the distance of the detected point are in a large wave gate (the predicted position of the previous circle is taken as the center), if the detected point is in the large wave gate and the identity codes are consistent, setting the point track correlation score as y5 and setting the type of the relevant wave gate as the large wave gate; if the point track correlation score is within the big wave gate and the identity codes are not consistent, setting the point track correlation score as y6 and setting the correlation wave gate type as the big wave gate; if the detected point is not in the big wave gate, whether the azimuth and the distance of the detected point are in the big wave gate (the updated position of the previous circle is taken as the center and is calculated according to the speed, the time, the acceleration, the turning coefficient and the deceleration coefficient) is continuously judged, if the detected point is in the big wave gate and the identity codes are consistent, the point track related fraction is y7(y1> y3> y5> y2> y4> y6> y7), and the related wave gate type is the big wave gate; otherwise, the traversal is regarded as impossible to be related to the target position, and the traversal can be quitted.
Under the condition that the point track related fraction is larger than 0, judging whether the amplitude difference is in the range of 0 dB-1 dB, and if so, adding the point track related fraction to an auxiliary score y 8; if the amplitude difference is within the range of 1 dB-3 dB, if so, adding the auxiliary score y9 to the point track related score (y8> y 9); if the amplitude difference is larger than 3dB, the relevant fraction is not processed, and the next step is continuously executed;
and under the condition that the point track related score is larger than 0 and the height is effective, judging whether the height difference is in the range of the interval between the altitudes, if so, adding the point track related score to an auxiliary score y10, otherwise, not processing the point track related score, and continuing to execute the next step.
When the detected point track is within the range of 10km, if other target data exist, traversing other point track data except the detected point track in the point track sector in a reverse direction by taking a current track node as a base point, matching identity codes, resetting the point track related fraction of the current track node if the identity codes are successfully matched, and if the identity codes are unsuccessfully matched, continuing to execute identity code matching processing until all other point tracks in the point track sector are traversed and finished;
and if no other target data exists in the range of 10km of the detected point track, the point track related score is not processed, and the current traversal exits.
Point track correlation processing of extrapolated tracks, as shown in fig. 4:
when the identity codes are consistent, continuously judging whether the azimuth and the distance of the detection point are in a wavelet gate (the predicted position of the previous circle is taken as the center), if so, setting the point track related fraction as m1 and setting the related wavelet gate type as the wavelet gate; if the detected point is not in the wavelet gate, continuously judging whether the azimuth and the distance of the detected point are in the middle wave gate (the predicted position of the previous circle is taken as the center), if so, setting the point track related fraction as m2 and setting the type of the related wave gate as the middle wave gate; if the detected point is not in the medium wave gate, continuously judging whether the azimuth and the distance of the detected point are in a large wave gate (the predicted position of the previous circle is taken as the center), if so, setting the point track related fraction as m3 and setting the type of the related wave gate as the large wave gate; if the point trajectory correlation score is not in the grand gate, continuously judging whether the azimuth and the distance of the detection point are in the grand gate (the updated position of the previous circle is taken as the center), if so, setting the point trajectory correlation score as m4(m1> m2> m3> m4) and setting the correlation gate type as the grand gate; if the target position is not within the extra-large wave gate, the traversal is regarded as impossible to be related to the target position, and the traversal can be quitted. When the identity codes are inconsistent and other trace point data exist in the range of 3km of the target, the loop is exited; when the identity codes are inconsistent and no other trace data exists within the range of 3km of the target, the same processing as the identity code consistency is carried out.
Under the condition that the point track related fraction is larger than 0, judging whether the amplitude difference is in the range of 0 dB-1 dB, and if so, adding an auxiliary score m5 to the point track related fraction; if the amplitude difference is within the range of 1 dB-3 dB, if so, adding the auxiliary score m6 to the point track related score (m5> m 6); if the amplitude difference is larger than 3dB, the relevant fraction is not processed, and the next step is continuously executed;
and under the condition that the point track correlation score is larger than 0 and the height is effective, judging whether the altitude difference is in the altitude layer interval range, if so, adding an auxiliary score m7 to the point track correlation score, otherwise, not processing the point track correlation score, and exiting the traversal.
S4, counting the point track correlation scores of the effective track nodes of the current track sector, and finding the track node with the maximum point track correlation score as the optimal point track matching target, so as to facilitate the subsequent track updating processing; the specific process is as follows:
step 41: when the track points and track nodes in the sector are traversed, counting track point related scores of the track nodes in the sector, and if a maximum score value is found, the track point is the best point track matching;
step 42: and (3) clearing the point track related scores of other track nodes in the track sector and setting the types of gates as invalid, only keeping the track node with the maximum point track point related score, and giving a position index value of the best matching track node so as to be applied when the subsequent track is updated.
As shown in fig. 5, in the dual-machine hot backup machine, the method of the present invention is applied to the processing part related to the point track in the track processing system of the machine a, the traditional method is adopted to the processing part related to the point track in the track processing system of the machine B, and on the premise that other track processing systems are not changed, a large amount of data is collected to compare the results obtained by statistics, as can be seen from fig. 5, when the identity code is switched during the cross-air traffic control airspace, the success rate of target matching is 100%, no track change occurs, the multi-target error matching rate is 36.78% less, and the track extrapolation rate is 41.32% less.
The invention provides the method which uses the position as the basis for correlation and uses the identity code, the signal amplitude and the height as the auxiliary correlation, thereby greatly improving the probability of successful target correlation, so that the identity code switching of the target occurs when the airspace is crossed with the navigation management, and the phenomena of track breakage and track batch number switching of the target can not occur.
The invention provides a method for matching a flight path based on position-based correlation and identity codes, signal amplitude and height as auxiliary correlation, and on the premise that multi-target correlation exists, the flight path is matched based on a point path, and meanwhile, the point path is reversely matched by a flight path node, so that error correlation is greatly reduced, and the accuracy of point-flight path correlation is improved.
According to the design of the invention, the size of the related wave gate is adjusted in a self-adaptive manner according to the track node attribute, the track point azimuth quality, the identity code and the extrapolation frequency, the related correct probability is increased to the maximum extent, the related wave gate indication is assigned, the phenomena of track extrapolation and track restart batch are greatly reduced, and the related wave gate indication can also provide reference for a track filtering coefficient.
With the application of the invention on the interrogator equipment, the detection observation effect of the secondary radar is improved, the performance of the equipment is improved, the high recognition of the user is obtained, the workload of the operator on duty is reduced, and the powerful guarantee is provided for ensuring the flight safety of the airplane.
Claims (6)
1. A method for secondary radar point track correlation processing, comprising the steps of:
s1, receiving trace data according to the sector and changing the coordinate system of the position information;
s2, traversing effective track nodes in the track sector, calculating identity matching marks of the current track and the current track nodes, and adaptively adjusting a relevant wave gate threshold;
s3, calculating the position difference between the current detection point track and the current track node, and determining a point track correlation score according to the position difference and a correlation gate threshold;
s4, counting the point track correlation scores of the effective track nodes of the current track sector, and finding the track node with the maximum point track correlation score as the optimal point track matching target, so as to facilitate the subsequent track updating processing.
2. The method for secondary radar spot track correlation processing according to claim 1, wherein the specific steps of S1 are as follows:
s11, receiving trace point data sent by the trace point processing system according to the sector;
and S12, transforming the space coordinate system of the received trace point data.
3. The method as claimed in claim 2, wherein the elimination of false targets is completed in the trace processing system, the difference between the distances of the false targets and the sector is 10km, the trace data without the same identity code is obtained, the calculated azimuth quality of the trace is carried after the trace processing, and the trace data of other targets in the range of 3km and 10km of the trace is counted.
4. The method for secondary radar spot track correlation processing according to claim 1, wherein the specific steps of S2 are as follows:
s21, traversing with the effective track node in the corresponding track sector after receiving the trace point data, calculating the time difference, exiting the traversal if the time difference does not meet the track requirement, continuing to calculate whether the identity code carried by the trace point data is matched with the identity code of the current track node if the time difference meets the requirement, setting the identity matching identifier to be 1 if the matching is successful, and clearing the identity matching identifier to be 0 if the matching is not successful;
s22, aiming at the immature track and the mature track, self-adaptively adjusting a related wave gate threshold according to the identity matching identifier and the point track azimuth quality; and aiming at the extrapolated tracks, adaptively adjusting the threshold of the related wave gate according to the extrapolation times, the identity matching identifier and the point track azimuth quality, and dividing the wave gate into a small wave gate, a medium wave gate, a large wave gate and an extra-large wave gate through the threshold of the related wave gate.
5. The method for secondary radar spot track correlation processing according to claim 1, wherein the specific steps of S3 are as follows:
s31, calculating a distance difference value, an orientation difference value, an amplitude difference value and a height difference value between the track position and the historical position or the predicted position of the track node; a height difference exists on the premise that the height value is effective;
s32, matching and comparing the calculated distance difference and the direction difference with a small wave gate, a middle wave gate, a big wave gate and a big wave gate in sequence, if the matching is successful, setting different point track relevant basic scores and setting a wave gate type identifier, if the identity matching identifier is 1 and the point track relevant basic Score is greater than zero, adding an identity code auxiliary Score code to the point track relevant basic Score to obtain a point track relevant Score, otherwise, entering the step S33;
s33, performing signal amplitude correlation matching on the premise that the point track correlation Score is larger than zero, adding the point track correlation Score to the signal amplitude auxiliary Score Score1 amplitude if the matching result is within 1dB, adding the point track correlation Score to the Score2 amplitude if the matching result is within 3dB, and otherwise, entering the step S34;
s34, on the premise that the point track correlation Score is not zero and the height value is effective, if the height difference value is in the height layer interval of civil aviation planning, accumulating the point track correlation Score and the height of the height correlation auxiliary Score, otherwise, entering the step S35;
s35, if the identity matching identifier is 1, completing the matching scoring mechanism, and exiting the traversal; if the identity matching identifier is 0, continuously judging whether other target point track data exist in the trace point data 10km or not, if not, exiting the traversal, if the trace point data exist in the trace point data 10km, carrying out identity code consistency matching on the current track node and other trace point data which are not the current trace point in the sector, if identity codes are consistent, resetting the score related to the current point track, if the identity codes are different, exiting the traversal, and continuously executing the next circulation processing.
6. The method for secondary radar spot track correlation processing according to claim 1, wherein the specific steps of S4 are as follows:
s41, when the track points and the track nodes in the sector are traversed, counting the track point related scores of the track nodes in the sector, and finding the track node with the largest track point related score as the best track point matching result;
s42, clearing the point track related scores of other track nodes in the track sector, setting the wave gate type as invalid, only keeping the track node with the maximum point track related score, and giving the position index of the best matching track node.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210930747.1A CN114994633B (en) | 2022-08-04 | 2022-08-04 | Method for secondary radar point track correlation processing |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210930747.1A CN114994633B (en) | 2022-08-04 | 2022-08-04 | Method for secondary radar point track correlation processing |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114994633A true CN114994633A (en) | 2022-09-02 |
| CN114994633B CN114994633B (en) | 2022-11-08 |
Family
ID=83023247
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210930747.1A Active CN114994633B (en) | 2022-08-04 | 2022-08-04 | Method for secondary radar point track correlation processing |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114994633B (en) |
Citations (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5001650A (en) * | 1989-04-10 | 1991-03-19 | Hughes Aircraft Company | Method and apparatus for search and tracking |
| US20060268111A1 (en) * | 2005-05-31 | 2006-11-30 | Objectvideo, Inc. | Multi-state target tracking |
| CA2585023A1 (en) * | 2006-11-09 | 2008-05-09 | Raytheon Canada Limited | Track quality based multi-target tracker |
| CN101655561A (en) * | 2009-09-14 | 2010-02-24 | 南京莱斯信息技术股份有限公司 | Federated Kalman filtering-based method for fusing multilateration data and radar data |
| CN103076605A (en) * | 2013-01-06 | 2013-05-01 | 中国电子科技集团公司第十研究所 | Secondary surveillance radar track extraction method for multimode polling and S-mold roll-calling interrogation |
| CN104749576A (en) * | 2015-04-08 | 2015-07-01 | 四川九洲空管科技有限责任公司 | Multi-radar track association and fusion method |
| CN104808197A (en) * | 2015-05-06 | 2015-07-29 | 四川九洲空管科技有限责任公司 | Multi-surveillance-source flying target parallel track processing method |
| CN104991235A (en) * | 2015-06-15 | 2015-10-21 | 南京航空航天大学 | Method for rapid tracking target based on radar trace points |
| CN105510896A (en) * | 2015-12-03 | 2016-04-20 | 成都能通科技有限公司 | Weighted nearest-neighbor data association method for centralized multi-radar data processing process |
| CN105807280A (en) * | 2016-04-26 | 2016-07-27 | 南京鹏力系统工程研究所 | Echo fused target track association method based on track state estimation |
| CN106934324A (en) * | 2015-12-30 | 2017-07-07 | 南京理工大学 | Based on the radar data correlating methods for simplifying many hypothesis algorithms |
| CN107063259A (en) * | 2017-03-08 | 2017-08-18 | 四川九洲电器集团有限责任公司 | A kind of Data Association and electronic equipment |
| CN107656266A (en) * | 2017-08-29 | 2018-02-02 | 安徽四创电子股份有限公司 | S mode secondary radar information fusion system and its fusion method based on ADS B |
| CN107977987A (en) * | 2017-11-20 | 2018-05-01 | 北京理工大学 | A kind of UAV system multiple target detection tracking, instruction system and method |
| CN109143181A (en) * | 2018-07-25 | 2019-01-04 | 南京理工大学 | A kind of radar plot Trace Association method convenient for programming |
| CN110824467A (en) * | 2019-11-15 | 2020-02-21 | 中山大学 | Multi-target tracking data association method and system |
| CN112346096A (en) * | 2020-11-10 | 2021-02-09 | 中国人民解放军海军航空大学 | High-low orbit remote sensing satellite ship target track and point track correlation method and system |
| CN112946626A (en) * | 2021-03-11 | 2021-06-11 | 中国电子科技集团公司第三十八研究所 | Airborne phased array radar track correlation method |
| CN113064155A (en) * | 2021-03-18 | 2021-07-02 | 沈阳理工大学 | Optimization method of track association under multi-target tracking of aerial radar |
| CN114670851A (en) * | 2022-04-20 | 2022-06-28 | 北京理工大学重庆创新中心 | Driving assistance system, method, terminal and medium based on optimizing tracking algorithm |
| CN114740466A (en) * | 2022-03-30 | 2022-07-12 | 南通大学 | Radar target trace filtering method based on scoring mechanism |
-
2022
- 2022-08-04 CN CN202210930747.1A patent/CN114994633B/en active Active
Patent Citations (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5001650A (en) * | 1989-04-10 | 1991-03-19 | Hughes Aircraft Company | Method and apparatus for search and tracking |
| US20060268111A1 (en) * | 2005-05-31 | 2006-11-30 | Objectvideo, Inc. | Multi-state target tracking |
| CA2585023A1 (en) * | 2006-11-09 | 2008-05-09 | Raytheon Canada Limited | Track quality based multi-target tracker |
| CN101655561A (en) * | 2009-09-14 | 2010-02-24 | 南京莱斯信息技术股份有限公司 | Federated Kalman filtering-based method for fusing multilateration data and radar data |
| CN103076605A (en) * | 2013-01-06 | 2013-05-01 | 中国电子科技集团公司第十研究所 | Secondary surveillance radar track extraction method for multimode polling and S-mold roll-calling interrogation |
| CN104749576A (en) * | 2015-04-08 | 2015-07-01 | 四川九洲空管科技有限责任公司 | Multi-radar track association and fusion method |
| CN104808197A (en) * | 2015-05-06 | 2015-07-29 | 四川九洲空管科技有限责任公司 | Multi-surveillance-source flying target parallel track processing method |
| CN104991235A (en) * | 2015-06-15 | 2015-10-21 | 南京航空航天大学 | Method for rapid tracking target based on radar trace points |
| CN105510896A (en) * | 2015-12-03 | 2016-04-20 | 成都能通科技有限公司 | Weighted nearest-neighbor data association method for centralized multi-radar data processing process |
| CN106934324A (en) * | 2015-12-30 | 2017-07-07 | 南京理工大学 | Based on the radar data correlating methods for simplifying many hypothesis algorithms |
| CN105807280A (en) * | 2016-04-26 | 2016-07-27 | 南京鹏力系统工程研究所 | Echo fused target track association method based on track state estimation |
| CN107063259A (en) * | 2017-03-08 | 2017-08-18 | 四川九洲电器集团有限责任公司 | A kind of Data Association and electronic equipment |
| CN107656266A (en) * | 2017-08-29 | 2018-02-02 | 安徽四创电子股份有限公司 | S mode secondary radar information fusion system and its fusion method based on ADS B |
| CN107977987A (en) * | 2017-11-20 | 2018-05-01 | 北京理工大学 | A kind of UAV system multiple target detection tracking, instruction system and method |
| CN109143181A (en) * | 2018-07-25 | 2019-01-04 | 南京理工大学 | A kind of radar plot Trace Association method convenient for programming |
| CN110824467A (en) * | 2019-11-15 | 2020-02-21 | 中山大学 | Multi-target tracking data association method and system |
| CN112346096A (en) * | 2020-11-10 | 2021-02-09 | 中国人民解放军海军航空大学 | High-low orbit remote sensing satellite ship target track and point track correlation method and system |
| CN112946626A (en) * | 2021-03-11 | 2021-06-11 | 中国电子科技集团公司第三十八研究所 | Airborne phased array radar track correlation method |
| CN113064155A (en) * | 2021-03-18 | 2021-07-02 | 沈阳理工大学 | Optimization method of track association under multi-target tracking of aerial radar |
| CN114740466A (en) * | 2022-03-30 | 2022-07-12 | 南通大学 | Radar target trace filtering method based on scoring mechanism |
| CN114670851A (en) * | 2022-04-20 | 2022-06-28 | 北京理工大学重庆创新中心 | Driving assistance system, method, terminal and medium based on optimizing tracking algorithm |
Non-Patent Citations (10)
| Title |
|---|
| WANG, MH: "A tracking gate algorithm based on optimal performance of data association", 《2000 5TH INTERNATIONAL CONFERENCE ON SIGNAL PROCESSING PROCEEDINGS》 * |
| WEI S H: "A Performance Optimized Tracking Gate Algorithm Based on Data Association", 《MECHANICAL ENGINEERING & AUTOMATION》 * |
| ZHANG CHENG: "Association Gate Algorithm for Multi-target Tracking Based on Association Performance EvaluationTracking Associat gate", 《JOURNAL OF EAST CHINA UNIVERSITY OF SCIENCE AND TECHNOLOGY (NATURAL SCIENCE EDITION)》 * |
| 何桂萍: "ADS-B与雷达组合监视数据融合方法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
| 冯涛: "空管一次雷达抗风电场干扰目标检测方法", 《雷达科学与技术》 * |
| 刘高峰: "多目标航迹处理中相关波门的确定", 《火力与指挥控制》 * |
| 周敏: "多目标跟踪中航迹关联技术研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
| 李春雨: "基于无源时差定位的目标跟踪算法研究与实现", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
| 莫然等: "基于小波滤波的JPDA算法多目标跟踪", 《沈阳化工大学学报》 * |
| 高又松: "基于DSP的雷达多目标相关处理", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114994633B (en) | 2022-11-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107255803B (en) | Ship track logic initiation algorithm | |
| CN112612862B (en) | Grid map positioning method based on point cloud registration | |
| CN114859339A (en) | Multi-target tracking method based on millimeter wave radar | |
| CN111832462B (en) | Frequency hopping signal detection and parameter estimation method based on deep neural network | |
| CN116153078B (en) | Road safety assessment method and device based on millimeter wave radar and storage medium | |
| CN112505681B (en) | Four-side two-dimensional phased array radar multi-target tracking processing method based on TAS | |
| CN113993172B (en) | Ultra-dense network switching method based on user movement behavior prediction | |
| CN110806575A (en) | Cooperative and non-cooperative unmanned aerial vehicle identification method and system based on multi-source information | |
| CN109239702B (en) | Airport low-altitude flying bird number statistical method based on target state set | |
| CN109143192A (en) | Multiple target effectively originates confirmation management method, system and computer readable storage medium | |
| CN114994633B (en) | Method for secondary radar point track correlation processing | |
| CN109270520B (en) | Processing method for acquiring secondary radar response target identity code based on amplitude information | |
| CN116520310A (en) | Maneuvering multi-target tracking method under Doppler blind area | |
| CN109828253B (en) | Multi-station radar quantitative fusion target detection method | |
| CN109656271B (en) | Track soft association method based on data association idea | |
| CN115542307B (en) | High-speed scene multi-radar track fusion method based on high-precision map | |
| CN115119141B (en) | UWB positioning method based on Kalman filtering and applicable to complex indoor environment | |
| CN112083388B (en) | Trace point quality evaluation method based on secondary radar response signal | |
| CN115508787A (en) | Sliding window water-flowing type point trace condensation processing method and system based on three-coordinate radar | |
| CN113359100A (en) | Attribute data association method for tracking target radar radiation source | |
| CN115097437A (en) | Underwater target tracking track approaching and crossing solution method based on label multi-Bernoulli pre-detection tracking algorithm | |
| CN109978170B (en) | Mobile equipment identification method based on multiple elements | |
| CN110286363B (en) | Target long-time tangential flight track speed difference re-correction method | |
| RU2716495C1 (en) | Method and system for multi-purpose tracking in two-position radar systems | |
| CN112491393B (en) | Linear filtering method based on unknown covariance matrix of observed noise |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |