Track vertical segmentation and missing time calculation method and device
By classifying and computing aircraft flight data, the flight time and missing time of aircraft at different altitude levels are evaluated, and the problem of insufficient rigor of calculating missing time in the prior art is solved, and effective positioning and solving aircraft communication interruptions are achieved. .
Patent Information
- Application Number
- CN202410960568.1
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-07-17
- Publication Date
- 2024-12-10
AI Technical Summary
The prior art lacks rigor in calculating the flight time and missing time of aircraft in different altitude layers, and cannot effectively evaluate the continuity of surveillance data, making it difficult to locate and solve the communication interruption problem.
By acquiring the aircraft's flight data, using 24-bit address codes to classify and calculate the time for flights to enter and leave the airspace, calculate the loss time of the track point based on the monitoring data update frequency and air pressure altitude, and evaluate the continuity of different altitude layers.
Accurate calculation of the flight time and missing time of the aircraft at different altitude levels is achieved, the continuity of surveillance data is evaluated, and the communication interruption problem is helped to ensure the normal processing of aircraft signals.
Smart Images
Figure CN119107849A_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to the technical field of civil aviation data processing, and in particular to a method for calculating vertical segmentation of a track and missing time, and also to a device for calculating vertical segmentation of a track and missing time, which uses aircraft track information to calculate actual flight time and missing time at different high levels within a given airspace. Background Art
[0002] In order to ensure the safety of aircraft during flight, it is necessary to monitor the aircraft's track information in real time, and the monitoring information can reflect the aircraft's flight route and status. By analyzing the track information, we can better understand the aircraft status and air traffic conditions.
[0003] At present, aviation surveillance technology in China has been relatively perfected. Civil aircraft rely on S-mode radar and ADS-B as the main means of surveillance, sending track information to the ground system, which receives and processes it. Because the current surveillance system has the characteristic of a small message time interval, it is suitable for real-time surveillance of aircraft to obtain flight status and reflect air traffic conditions, so as to better carry out air traffic control work.
[0004] Civil aviation controllers will assign different flight altitudes to aircraft based on different aircraft performance, weather and terrain characteristics of the area, flight status, etc. According to the flight altitude layer allocation standard determined by the Civil Aviation Administration of China in 2007, an altitude layer will be set every 300 meters for flights between 600 meters and 8,400 meters, an altitude layer will be set every 8,400 meters to 8,900 meters, an altitude layer will be set every 300 meters for flights between 8,900 meters and 12,500 meters, and an altitude layer will be set every 600 meters above 12,500 meters.
[0005] In actual operation, since the airborne and ground equipment are interfered by terrain, software and hardware during communication, the signal may be lost at a specified altitude in a certain area. If this happens, it is necessary to locate the area and altitude where the missing time occurs to troubleshoot the problem, determine the cause of the missing time problem, and ensure the normal communication and processing of aircraft signals by the airborne equipment and ground system.
[0006] The existing methods mainly determine the continuity of the track from a horizontal perspective, and do not perform rigorous calculations and analyses on the time of missing signals at different altitudes. Summary of the invention
[0007] In order to overcome the defects of the prior art, the technical problem to be solved by the present invention is to provide a method for calculating the vertical segmentation of the track and the missing time, which can accurately calculate the flight time and missing time of each altitude layer passed by the flight, and evaluate the continuity of the monitoring data in a given area at different altitude layers.
[0008] The technical solution of the present invention is: the method for calculating the vertical segmentation of the track and the missing time comprises the following steps: (I) Obtain all flight data and classify the flight data according to the 24-bit address code; (II) obtain flight data within a specified airspace area; (III) calculate the time each flight enters and leaves the airspace; (IV) Calculate, based on the existing monitoring data update frequency, whether the time interval between two adjacent messages of each flight meets the established frequency requirements; (V) The track points that do not meet the established update frequency requirements are sorted out, and the time lost by the adjacent track points through these altitude layers is calculated based on the pressure altitudes of the adjacent track points.
[0009] The present invention obtains the track data received by the ground monitoring system, obtains the track points of each flight in the specified airspace, can accurately calculate the flight time and missing time of each altitude layer passed by the flight, and evaluate the continuity of the monitoring data at different altitude layers in a given area.
[0010] Another method for calculating the vertical segmentation of the track and the missing time is also provided, comprising the following steps: (1) In the specified rectangular airspace, the 24-bit address code of the aircraft is used as the unique identifier of the aircraft to extract the track. The total number of tracks is m, and j = 1; (2) If j <= m, execute step (3); otherwise, execute step (12); (3) Take the jth track and count the number of track points as n; (4) If n<13, execute step (11); otherwise, execute step (5); (5) Calculate the start time t of the jth track entering the airspace in ; (6) Calculate the deadline t for the jth track to leave the airspace out ; (7) Get {P1,…P i ,…P n The time of each point {t1,…t i ,…t n} and pressure altitude {fl1,…fl i ,…fl n}; (8) Calculate the FLT of each altitude layer between the time when the jth track enters the airspace and the time corresponding to the first track. in ; (9) Calculate the FLT of each altitude layer between the time corresponding to the last track point in the airspace and the time of leaving the airspace for the jth track out ; (10) Calculate the loss time vector of each altitude layer between the first track point and the last track point of the jth track in the airspace, denoted as MT j ; (11) Let j = j + 1 and execute step (2). (12) End.
[0011] A device for calculating vertical segmentation of a track and missing time is also provided, comprising: a first acquisition module configured to acquire all flight data and classify the flight data according to a 24-bit address code; A second acquisition module configured to acquire flight data within a specified airspace range; A first calculation module configured to calculate the time when each flight enters and leaves the airspace; A second calculation module is configured to calculate whether the time interval between two adjacent messages of each flight meets the established frequency requirement according to the existing monitoring data update frequency; The third calculation module is configured to sort out the track points that do not meet the established update frequency requirements, and calculate the time when the adjacent tracks pass through these altitude layers based on the pressure altitudes of the adjacent track points.
[0012] Another device for calculating the vertical segmentation and missing time of a track is also provided, comprising: A track extraction module is configured to extract tracks within a specified rectangular airspace using the aircraft's 24-bit address code as the aircraft's unique identifier; The entry and exit airspace time calculation module is configured to calculate the start time of a certain track entering the airspace and the end time of leaving the airspace; The lost time vector calculation module is configured to obtain the time and pressure altitude of each point in the track point, calculate the lost time of each altitude layer from the time of entering the airspace to the time corresponding to the first track point, and the lost time of each altitude layer from the time corresponding to the last track point in the airspace to the time of leaving the airspace, and calculate the lost time vector of each altitude layer between the first track point and the last track point of the track in the airspace. BRIEF DESCRIPTION OF THE DRAWINGS
[0013] Figure 1Shown is an overall flow chart of a method for vertical segmentation of a track and calculation of missing time according to the present invention.
[0014] Figure 2 Shown is an overall flow chart of another method for vertical track segmentation and missing time calculation according to the present invention.
[0015] Figure 3 Shown Figure 2 Flowchart of step (5).
[0016] Figure 4 Shown Figure 2 Flowchart of step (6).
[0017] Figure 5 Shown Figure 2 Flow chart of step (8).
[0018] Figure 6 Shown Figure 2 Flow chart of step (9).
[0019] Figure 7 Shown Figure 2 Flow chart of step (10).
[0020] Figure 8 A schematic diagram showing the situation where a track point is missing when an aircraft enters the airspace. DETAILED DESCRIPTION
[0021] In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and specific embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.
[0022] In order to make the description of the present disclosure more detailed and complete, the following provides an illustrative description of the implementation and specific embodiments of the present invention; however, this is not the only form of implementing or using the specific embodiments of the present invention. The implementation covers the features of multiple specific embodiments and the method steps and sequences used to construct and operate these specific embodiments. However, other specific embodiments can also be used to achieve the same or equivalent functions and step sequences.
[0023] First, the following are the parameter symbols and their meanings used in the present invention: FRQ is the time interval constant for sending messages specified by the ADS-B system; P i , P j are two consecutive track points, the i-th point and the j-th point; t i , t jis the time between two consecutive track points, the i-th point and the j-th point; t in is the time at which the aircraft is inferred to enter the airspace; t out is the time at which the aircraft is inferred to have left the airspace; fl i , fl j is the pressure altitude value of the i-th and j-th two consecutive track targets; fl in It is the pressure altitude value inferred when the aircraft enters the airspace; fl out It is the pressure altitude value inferred when the aircraft leaves the airspace; ta i is the steering angle value of the i-th point track target; gs i is the ground speed of the ith track point; dis i is the minimum distance from the ith track point to the boundary; FRQFL i is the altitude layer value corresponding to the pressure altitude of the ith trajectory point; FL in It is the altitude layer value corresponding to the pressure altitude when entering the airspace; FL out It is the altitude layer value corresponding to the pressure altitude when leaving the airspace; FLT in Yes in The time lost corresponding to the altitude layers passed between t1 and t2; FLT out Yes n to out The time lost corresponding to the altitude layers passed between them; FL ijn is the adjacent trajectory point P i , P j The height value of the nth altitude layer between; FLT ijn is the adjacent trajectory point P i , P j There is a loss time of n altitude layers between them; MT j is the vector of the sum of the missing time of the jth aircraft at each altitude level.
[0024] like Figure 1 As shown in FIG. 1 , the method for calculating the vertical segmentation of the track and the missing time includes the following steps: (I) Obtain all flight data and classify the flight data according to the 24-bit address code; (II) obtain flight data within a specified airspace area; (III) calculate the time each flight enters and leaves the airspace; (IV) Calculate, based on the existing monitoring data update frequency, whether the time interval between two adjacent messages of each flight meets the established frequency requirements; (V) The track points that do not meet the established update frequency requirements are sorted out, and the time lost in these altitude layers passed by the adjacent track points is calculated based on the pressure altitudes of the adjacent track points.
[0025] The present invention obtains the track data received by the ground monitoring system, obtains the track points of each flight in the specified airspace, can accurately calculate the flight time and missing time of each altitude layer passed by the flight, and evaluate the continuity of the monitoring data at different altitude layers in a given area.
[0026] like Figure 2 As shown, another method for calculating the vertical segmentation of the track and the missing time is provided, comprising the following steps: (1) In the specified rectangular airspace, the 24-bit address code of the aircraft is used as the unique identifier of the aircraft to extract the track. The total number of tracks is m, and j = 1; (2) If j <= m, execute step (3); otherwise, execute step (12); (3) Take the jth track and count the number of track points as n; (4) If n<13, execute step (11); otherwise, execute step (5); (5) Calculate the start time t of the jth track entering the airspace in ; (6) Calculate the deadline t for the jth track to leave the airspace out ; (7) Get {P1,…P i ,…P n The time of each point {t1,…t i ,…t n} and pressure altitude {fl1,…fl i ,…fl n}; (8) Calculate the FLT of each altitude layer between the time when the jth track enters the airspace and the time corresponding to the first track. in ; (9) Calculate the FLT of each altitude layer between the time corresponding to the last track point in the airspace and the time of leaving the airspace for the jth track out ; (10) Calculate the loss time vector of each altitude layer between the first track point and the last track point of the jth track in the airspace, denoted as MT j ; (11) Let j = j + 1 and execute step (2). (12) End.
[0027] Preferably, if Figure 3 As shown, the step (5) includes the following sub-steps: (5.1) Obtain the speed gs1 of the first track point. If the speed gs1 < 350 km / h, the aircraft did not lose any track points when entering the airspace, and execute step (5.17). Otherwise, execute step (5.2). (5.2) Calculate the minimum distance dis from the first track point P1 to the airspace boundary line in , if dis in If it is less than 3×FRQ×gs1, execute step (5.3); otherwise, execute step (5.7); (5.3) Take the 2nd to 13th points among the n track points as {P2,P 3, P i …,P 13}, let i=2; (5.4) If i=13, execute step (5.7), otherwise execute step (5.5); (5.5) Calculate P i Minimum distance to the airspace boundary line dis i , if dis i If it is less than 3×FRQ×gs1, execute step (5.6); otherwise, execute step (5.17); (5.6) Let i = i + 1 and execute step (5.4); (5.7) When the aircraft enters the airspace, there is a missing track point. The vertical distance between the P1 track point and the upper boundary of the rectangle is d1, the vertical distance between the right boundary is d2, the vertical distance between the lower boundary is d3, and the vertical distance between the left boundary is d4. Figure 8 As shown; (5.8) Get the heading angle ta1 of the P1 track point. If ta1∈(0,π / 2), execute step (5.9); if ta1∈(π / 2,π), execute step (5.10); if ta1∈(π, 3π / 2), execute step (5.11); if ta1∈(3π / 2, 2π), execute step (5.12); if ta1= π / 2, execute step (5.13); if ta1=π, execute step (5.14); if =3π / 2, execute step (5.15); if ta1=2π, execute step (5.16); (5.9) The angle between P1 and E and due north is α: if ta1<α, then t in is t1-(d1 / |gs1cosθ|), if ta1>α, then t in is t1-(d2 / |gs1sinθ|), execute step (5.18); (5.10) The angle between P1 and F and due north is α: if ta1<α, then t in is t1-(d2 / |gs1sinθ|), if ta1>α then t in is t1-(d4 / |gs1sinθ|), execute step (5.18); (5.11) The angle between P1 and G and due north is α: if ta1<α, then t in is t1-(d4 / |gs1sinθ|), if ta1>α, then t in For t1-(d3 / |gs1sinθ|), execute step (5.18); (5.12) The angle between P1 and H and due north is α: if ta1<α, then t in is t1-(d4 / |gs1sinθ|), if ta1>α then t in is t1-(d1 / |gs1cosθ|), execute step (5.18); (5.13) If ta1 = π / 2, then t in =t1-(d1 / |gs1|); (5.14) If ta1 = π, then t in =t1-(d2 / |gs1|), execute step (5.18); (5.15) If ta1 = 3π / 2, then t in =t1-(d3 / |gs1|), execute step (5.18); (5.16) If ta1 = 2π, then t in =t1-(d4 / |gs1|), execute step (5.18); (5.17)t in =t1; (5.18) End.
[0028] Preferably, if Figure 4 As shown, the step (6) includes the following sub-steps: (6.1) Get the speed gs of the last track point n , if the speed gs n <350, execute step (6.7), otherwise execute step (6.2); (6.2) Calculate the minimum distance dis from the last track point to the airspace boundary line out , if dis out Less than 3×FRQ×gs n , execute step (6.3), otherwise execute step (6.7); (6.3) In a traversal manner, take {P n-13 ,…,P i ,…,P n}, let i=n-1; (6.4) If i=n-13, execute step (6.7), otherwise execute step (6.5); (6.5) Calculate the minimum distance dis from the ith track point to the airspace boundary line i , if dis i Less than 3×FRQ×gs n , execute step (6.6), otherwise execute step (6.17); (6.6) Let i = i-1 and execute step (6.4); (6.7) Get the vertical distance between track point P1 and the upper boundary of the rectangle as d1, the vertical distance between the right boundary as d2, the vertical distance between the lower boundary as d3, and the vertical distance between the left boundary as d4, as follows: Figure 8 As shown; (6.8) In case of missing time for an aircraft, obtain P n The turning angle ta of the track point n If n ∈(0,π / 2), execute step (6.9), if ta n ∈(π / 2,π), execute step (6.10), if ta n ∈(π, 3π / 2), execute step (6.11), if ta n ∈(3π / 2, 2π), execute step (6.12), if ta n = π / 2, execute step (6.13), if ta1=π, execute step (6.14), if = 3π / 2, execute step (6.15), if ta1=2π, execute step (6.16); (6.9) Connection P n The angle between point E and true north is α: if ta n <α, then t out t n +(d1 / |gs n cosθ|), if ta n >α, then t out t n +(d2 / |gs nsinθ|), execute step (6.18); (6.10) Connection P n The angle between point F and true north is α: if ta n <α, then t out t n +(d2 / |gs n sinθ|), if ta n >α then t out t n +(d3 / |gs n cosθ|), execute step (6.18); (6.11) Connection P n The angle between point G and true north is α: if ta n <α, then t out t n +(d3 / |gs n cosθ|), if ta n >α, then t out t n +(d4 / |gs n sinθ|), execute step (6.18); (6.12) Connection P n The angle between point H and true north is α: if ta n <α, then t out t n +(d4 / |gs n sinθ|); if ta n >α, then t out t n +(d1 / |gs n cosθ|), execute step (6.18); (6.13) If n = π / 2, then time t out =t n +(d1 / |gs n |), execute step (6.18); (6.14) If n =π, then time t out =t n +(d2 / |gs n |), execute step (6.18); (6.15) If n =3π / 2, then time t out =t n +(d3 / |gs n |), execute step (6.18); (6.16) Ifn = 2π, then the time t out = t n +(d4 / |gs n |), perform step (6.18); (6.17)t out = t n ; (6.18)Calculate to obtain t out , and end.
[0029] Preferably, as Figure 5 shown, step (8) includes the following sub-steps: (8.1)Obtain the barometric altitudes fl1 and fl2 of the first two track points P1 and P2. If fl1 > fl2, perform step (8.2); if fl1 < fl2, perform step (8.5); if fl1 = fl2, perform step (8.3); (8.2)Calculate the barometric altitude fl when entering the airspace in = [(t1 - t in )(fl1 - fl2)] / (t2 - t1) + fl1; According to fl in and fl1, find the corresponding altitude layer values FL in and FL1. When FL in = FL1, perform step (8.3); otherwise, perform step (8.4); (8.3)All the lost time is at the altitude layer where the P1 track point is located, and its lost time is FL in = (t1 - t in ), perform step (8.7); (8.4)When the track is in an ascending state, FL in is 1 FL altitude layer more than FL1, and there is lost time at the altitude adjacent to the P1 point. Therefore, fl in corresponds to FLT of the FL1 + 1 altitude layer in = |(t1 - t n ) × [1 - (FL1 - fl1) / (fl in - fl2)]|, FLT1 of the FL1 altitude layer = |(t1 - t n ) × [(FL1 - fl1) / (fl in - fl2)]|, perform step (8.7); (8.5)Calculate the barometric altitude fl when entering the airspace in = fl1 - (t1 - t in )(fl2 - fl1) / (t2 - t1). When fl in = fl1, perform step (8.3); otherwise, perform step (8.6); (8.6) When the track is in a descending state, fl in The corresponding altitude layer is one FL layer less than the FL layer corresponding to fl1, and the FLT of FL1-1 altitude layer in =|(t1-t n )×[1-(FL1-fl1) / (fl in -fl1)]|, FL1 altitude layer FLT1=|(t1-t n )×[(FL1-fl1) / (fl in -fl1)]|, execute step (8.7); (8.7) End.
[0030] Preferably, if Figure 6 As shown, the step (9) includes the following sub-steps: (9.1) Get the first two track points P n-1 and P n The pressure altitude fl n-1 and fl n , if fl n-1 >fl n , execute step (9.2), if fl n-1 <fl n , execute step (9.5), otherwise execute step (9.3); (9.2) Calculate the pressure altitude fl when entering the airspace out =fl n -|(t n -t out )(fl n-1 -fl n ) / (t n -t n -1)|, according to fl out and fl n Find the corresponding altitude layer value FL out and FL n , when FL out =FL n When , execute step (9.3), otherwise execute step (9.4); (9.3) The lost time is all in P n The altitude layer where the track point is located, and its loss time is FLT out =( t out -t n ), execute step (9.7); (9.4) When the track is in an ascending state, fl out The corresponding altitude layer ratio fl nThe corresponding FL has one more FL altitude layer, and the altitude adjacent to the point has lost time, fl out Corresponding to FL n +1 level FLT out =|(t out -t n )×[1-(FL n -fl n ) / (fl out -fl n )]|, FL n FLT n =|(t out -t n )×[(FL n -fl n ) / (fl out -fl n )]|, execute step (9.7); (9.5) Calculate the pressure altitude fl when entering the airspace in =fl1-(t1-t in )(fl2-fl1) / (t2-t1), when fl in = fl1, execute step (9.3), otherwise execute step (9.6); (9.6) When the track is in a descending state, fl out The corresponding altitude layer ratio fl n The corresponding FL is one less FL level, FL n -1 level FLT out =|(t out -t n )×[1-(FL n -fl n ) / (fl out -fl n )]|, FL n FLT n =|(t out -t n )×[(FL n -fl n ) / (fl out -fl n )]|, execute step (9.7); (9.7) End.
[0031] Preferably, if Figure 7 As shown, the step (10) includes the following sub-steps: (10.1) Take the jth track and count the number of track points as n; (10.2) Take two adjacent track points Pi ,P j Perform the calculation. Let \(i = 1\) and \(j = i + 1\). (10.3)If \(j>n\), then execute step (10.14); otherwise, execute step (10.4). (10.4)If \(t\) i -t j >FRQ, then execute step (10.5); otherwise, execute step (10.13). (10.5)Find the corresponding altitude level values FL i and FL j according to fl i and FL j . If FL i =FL j , then execute step (10.6); otherwise, execute step (10.7). (10.6)FL i =t i -t j -FRQ, and execute step (10.13). (10.7)Take all the altitude levels FL i between fl j . The altitude level where fl i-jn is located is \(m\), and the altitude level where fl i is located is \(n\). Let \(q = m\). j 所处高度层为n,令q=m; (10.8)If \(q = m\), then execute step (10.9); if \(q = n\), then execute step (10.10); if \(q < n\), then execute step (10.11); otherwise, execute step (10.13). (10.9)FLT i-jm =|(t i -t j -FRQ)(FL i -fl i ) / (fl i -fl j )|, and execute step (10.12). (10.10)FLT i-jn =|(t i -t j -FRQ)(FL i -fl i ) / (fl i -fl j )|, and execute step (10.12). (10.11)FLT i-jq =|(t i -t j-FRQ)(FL i-jq-1 -FL i-jq ) / (fl i -fl j )| , execute step (10.12); (10.12) Let q = q + 1 and execute step (10.8); (10.13) Let i = i + 1, j = j + 1, and execute step (10.3); (10.14) MT j =FLT in +FLT out +Σ n n=0 FLT i-jn , execute step (10.15); (10.15)End.
[0032] Those skilled in the art can understand that all or part of the steps in the above-mentioned embodiment method can be completed by instructing the relevant hardware through a program, and the program can be stored in a computer-readable storage medium. When the program is executed, it includes the steps of the above-mentioned embodiment method, and the storage medium can be: ROM / RAM, disk, CD, memory card, etc. Therefore, corresponding to the method of the present invention, the present invention also includes a device for calculating the vertical segmentation of the track and the missing time, which is usually represented in the form of functional modules corresponding to the steps of the method. The device includes: a first acquisition module configured to acquire all flight data and classify the flight data according to a 24-bit address code; A second acquisition module configured to acquire flight data within a specified airspace range; A first calculation module configured to calculate the time when each flight enters and leaves the airspace; A second calculation module is configured to calculate whether the time interval between two adjacent messages of each flight meets the established frequency requirement according to the existing monitoring data update frequency; The third calculation module is configured to sort out the track points that do not meet the established update frequency requirements, and calculate the time when the adjacent tracks pass through these altitude layers based on the pressure altitudes of the adjacent track points.
[0033] Another device for calculating the vertical segmentation and missing time of a track is also provided, comprising: A track extraction module is configured to extract tracks within a specified rectangular airspace using the aircraft's 24-bit address code as the aircraft's unique identifier; The entry and exit airspace time calculation module is configured to calculate the start time of a certain track entering the airspace and the end time of leaving the airspace; The lost time vector calculation module is configured to obtain the time and pressure altitude of each point in the track point, calculate the lost time of each altitude layer from the time of entering the airspace to the time corresponding to the first track point, and the lost time of each altitude layer from the time corresponding to the last track point in the airspace to the time of leaving the airspace, and calculate the lost time vector of each altitude layer between the first track point and the last track point of the track in the airspace.
[0034] The above description is only a preferred embodiment of the present invention and does not limit the present invention in any form. Any simple modification, equivalent changes and modifications made to the above embodiments based on the technical essence of the present invention are still within the protection scope of the technical solution of the present invention.
Claims
1. A method for calculating track vertical segmentation and missing time, characterized in that: The method comprises the following steps: (I) Obtain all flight data and classify the flight data according to the 24-bit address code; (II) obtain flight data within a specified airspace area; (III) calculate the time each flight enters and leaves the airspace; (IV) Calculate, based on the existing monitoring data update frequency, whether the time interval between two adjacent messages of each flight meets the established frequency requirements; (V) The track points that do not meet the established update frequency requirements are sorted out, and the time lost in these altitude layers passed by the adjacent track points is calculated based on the pressure altitudes of the adjacent track points.
2. A method for calculating track vertical segmentation and missing time, characterized in that: The method comprises the following steps: (1) In the specified rectangular airspace, the 24-bit address code of the aircraft is used as the unique identifier of the aircraft to extract the track. The total number of tracks is m, and j = 1; (2) If j <= m, execute step (3); otherwise, execute step (12); (3) Take the jth track and count the number of track points as n; (4) If n<13, execute step (11); otherwise, execute step (5); (5) Calculate the start time t of the jth track entering the airspace in ; (6) Calculate the deadline t for the jth track to leave the airspace out ; (7) Get {P1,…P i ,…P n The time of each point {t1,…t i ,…t n } and pressure altitude {fl1,…fl i ,…fl n }; (8) Calculate the FLT of each altitude layer between the time of entry into the airspace and the time corresponding to the first track point of the jth track. in ; (9) Calculate the FLT of each altitude layer between the time corresponding to the last track point of the jth track and the time of leaving the airspace. out ; (10) Calculate the loss time vector of each altitude layer between the first track point and the last track point of the jth track in the airspace, denoted as MT j ; (11) Let j = j + 1 and execute step (2). (12) End.
3. The method for calculating track vertical segmentation and missing time according to claim 2, characterized in that: The step (5) comprises the following sub-steps: (5.1) Obtain the speed gs1 of the first track point. If the speed gs1 < 350 km / h, the aircraft did not lose any track points when entering the airspace, and execute step (5.17). Otherwise, execute step (5.2). (5.2) Calculate the minimum distance dis from the first track point P1 to the airspace boundary line in , if dis in If it is less than 3×FRQ×gs1, execute step (5.3); otherwise, execute step (5.7); (5.3) Take the 2nd to 13th points among the n track points as {P2, P3…, P 13 }, let i=2; (5.4) If i=13, execute step (5.7), otherwise execute step (5.5); (5.5) Calculate P i Minimum distance to the airspace boundary line dis i , if dis i If it is less than 3×FRQ×gs1, execute step (5.6); otherwise, execute step (5.17); (5.6) Let i = i + 1 and execute step (5.4); (5.7) When the aircraft enters the airspace, there is a missing track point. The vertical distance between the P1 track point and the upper boundary of the rectangle is d1, the vertical distance between the right boundary is d2, the vertical distance between the lower boundary is d3, and the vertical distance between the left boundary is d4; (5.8) Get the heading angle ta1 of the P1 track point. If ta1∈(0,π / 2), execute step (5.9); if ta1∈(π / 2,π), execute step (5.10); if ta1∈(π, 3π / 2), execute step (5.11); if ta1∈(3π / 2, 2π), execute step (5.12); if ta1= π / 2, execute step (5.13); if ta1=π, execute step (5.14); if ta1=3π / 2, execute step (5.15); if ta1= 2π, execute step (5.16); (5.9) The angle between P1 and E and due north is α: if ta1<α, then t in is t1-(d1 / |gs1cosθ|), if ta1>α, then t in is t1-(d2 / |gs1sinθ|), execute step (5.18); (5.10) The angle between P1 and F and due north is α: if ta1<α, then t in is t1-(d2 / |gs1sinθ|), if ta1>α then t in is t1-(d4 / |gs1sinθ|), execute step (5.18); (5.11) The angle between P1 and G and due north is α: if ta1<α, then t in is t1-(d4 / |gs1sinθ|), if ta1>α, then t in For t1-(d3 / |gs1sinθ|), execute step (5.18); (5.12) The angle between P1 and H and due north is α: if ta1<α, then t in is t1-(d4 / |gs1sinθ|), if ta1>α then t in is t1-(d1 / |gs1cosθ|), execute step (5.18); (5.13) If ta1 = π / 2, then t in =t1-(d1 / |gs1|), execute step (5.18); (5.14) If ta1 = π, then t in =t1-(d2 / |gs1|), execute step (5.18); (5.15) If ta1 = 3π / 2, then t in =t1-(d3 / |gs1|), execute step (5.18); (5.16) If ta1 = 2π, then t in =t1-(d4 / |gs1|), execute step (5.18); (5.17)t in =t1; (5.18) End.
4. The method for calculating track vertical segmentation and missing time according to claim 3, characterized in that: The step (6) comprises the following sub-steps: (6.1) Get the speed gs of the last track point n , if the speed gs n <350, execute step (6.7), otherwise execute step (6.2); (6.2) Calculate the minimum distance dis from the last track point to the airspace boundary line out , if dis out Less than 3×FRQ×gs n , execute step (6.3), otherwise execute step (6.7); (6.3) In a traversal manner, take {P n-13 ,…,P i ,…,P n }, let i=n-1; (6.4) If i=n-13, execute step (6.7), otherwise execute step (6.5); (6.5) Calculate the minimum distance dis from the ith track point to the airspace boundary line i , if dis i Less than 3×FRQ×gs n , execute step (6.6), otherwise execute step (6.17); (6.6) Let i = i-1 and execute step (6.4); (6.7) Obtain the vertical distance between track point P1 and the upper boundary of the rectangle as d1, the vertical distance between the right boundary as d2, the vertical distance between the lower boundary as d3, and the vertical distance between the left boundary as d4; (6.8) If the aircraft is missing time, obtain P n The turning angle ta of the track point n If n ∈(0,π / 2), execute step (6.9), if ta n ∈(π / 2,π), execute step (6.10), if ta n ∈(π, 3π / 2), execute step (6.11), if ta n ∈(3π / 2, 2π), execute step (6.12), if ta n = π / 2, execute step (6.13), if ta n =π, execute step (6.14), if ta n =3π / 2, execute step (6.15), if ta n =2π, execute step (6.16); (6.9) Connection P n The angle between point E and true north is α: if ta n <α, then t out t n +(d1 / |gs n cosθ|), if ta n >α, then t out t n +(d2 / |gs n sinθ|), execute step (6.18); (6.10) Connection P n The angle between point F and true north is α: if ta n <α, then t out t n +(d2 / |gs n sinθ|), if ta n >α then t out t n +(d3 / |gs n cosθ|), execute step (6.18); (6.11) Connection P n The angle between point G and true north is α: if ta n <α, then t out t n +(d3 / |gs n cosθ|), if ta n >α, then t out t n +(d4 / |gs n sinθ|), execute step (6.18); (6.12) Connection P n The angle between point H and true north is α: if ta n <α, then t out t n +(d4 / |gs n sinθ|); if ta n >α, then t out t n +(d1 / |gs n cosθ|), execute step (6.18); (6.13) If n= π / 2, then time t out =t n +(d1 / |gs n |)], execute step (6.18); (6.14) If n =π, then time t out =t n +(d2 / |gs n |), execute step (6.18); (6.15) If n =3π / 2, then time t out =t n +(d3 / |gs n |), execute step (6.18); (6.16) If n =2π, then time t out =t n +(d4 / |gs n |), execute step (6.18); (6.17)t out =t n ; (6.18) Calculate t out ,Finish.
5. The method for calculating track vertical segmentation and missing time according to claim 4, characterized in that: The step (8) comprises the following sub-steps: (8.1) Obtain the barometric altitudes fl1 and fl2 of the first two track points P1 and P2. If fl1 > fl2, execute step (8.2); if fl1 < fl2, execute step (8.5); if fl1 = fl2, execute step (8.3). (8.2) Calculate the pressure altitude fl when entering the airspace in =[(t1-t in )(fl1-fl2)] / (t2-t1)+fl1; according to fl in and fl1 to find the corresponding altitude layer value FL in and FL1, when FL in =FL1, execute step (8.3), otherwise execute step (8.4); (8.3) The loss time is all at the altitude layer where the P1 track point is located, and the loss time is FL in =(t1-t in ), execute step (8.7); (8.4) When the track is in the ascending state, FL in There is one more FL level than FL1. There is a loss of time at the altitude adjacent to P1. in FLT corresponding to FL1+1 in =|(t1-t n )× [1-(FL1-fl1) / (fl in -fl2)]|, FL1 altitude layer FLT1=|(t1-t n )×[(FL1-fl1) / (fl in -fl2)]|, execute step (8.7); (8.5) Calculate the pressure altitude fl when entering the airspace in =fl1-(t1-t in )(fl2-fl1) / (t2-t1), when fl in = fl1, execute step (8.3), otherwise execute step (8.6); (8.6) When the track is in a descending state, fl in The corresponding altitude layer is one FL layer less than the FL layer corresponding to fl1, and the FLT of FL1-1 altitude layer in =|(t1-t n )×[1-(FL1-fl1) / (fl in -fl1)]|, FL1 altitude layer FLT1=|(t1-t n )×[(FL1-fl1) / (fl in -fl1)]|, execute step (8.7); (8.7) End.
6. The method for calculating track vertical segmentation and missing time according to claim 5, characterized in that: The said step (9) includes the following sub-steps: (9.1) Get the first two track points P n-1 and P n The pressure altitude fl n-1 and fl n , if fl n-1 >fl n , execute step (9.2), if fl n-1 <fl n , execute step (9.5), otherwise execute step (9.3); (9.2) Calculate the pressure altitude fl when entering the airspace out =fl n -|(t n -t out )(fl n-1 -fl n ) / (t n -t n -1)|, according to fl out and fl n Find the corresponding altitude layer value FL out and FL n , when FL out =FL n When , execute step (9.3), otherwise execute step (9.4); (9.3) The lost time is all in P n The altitude layer where the track point is located, and its loss time is FLT out =( t out -t n ), execute step (9.7); (9.4) When the track is in an ascending state, fl out The corresponding altitude layer ratio fl n The corresponding FL has one more FL altitude layer, and the altitude adjacent to the point has lost time, fl out Corresponding to FL n +1 level FLT out =|(t out -t n )×[1-(FL n -fl n ) / (fl out -fl n )]|, FL n FLT n =|(t out -t n )×[(FL n -fl n ) / (fl out -fl n )]|, execute step (9.7); (9.5) Calculate the pressure altitude fl when entering the airspace in =fl1-(t1-t in )(fl2-fl1) / (t2-t1), when fl in = fl1, execute step (9.3), otherwise execute step (9.6); (9.6) When the track is in a descending state, fl out The corresponding altitude layer ratio fl n The corresponding FL is one less FL level, FL n -1 level FLT out =|(t out -t n )×[1-(FL n -fl n ) / (fl out -fl n )]|, FL n FLT n =|(t out -t n )×[(FL n -fl n ) / (fl out -fl n )]|, execute step (9.7); (9.7) End.
7. The method for calculating track vertical segmentation and missing time according to claim 6, characterized in that: The said step (10) includes the following sub-steps: (10.1) Select the j-th track and count the number of track points as n. (10.2) Take two adjacent track points P i ,P j Perform calculations, let i=1, j=i+1; (10.3) If j > n, execute step (10.14); otherwise, execute step (10.4). (10.4) If t i -t j >FRQ, then execute step (10.5), otherwise execute step (10.13); (10.5) According to fl i and fl j Find the corresponding altitude layer value FL i and FL j , if FL i =FL j , then execute step (10.6), otherwise execute step (10.7); (10.6) FL i =t i -t j -FRQ, execute step (10.13); (10.7) Take fl i and fl j All levels between FL i-jn , fl i The altitude layer is m, fl j The altitude layer is n, let q = m; (10.8) If q = m, execute step (10.9); if q = n, execute step (10.10); if q < n, execute step (10.11); otherwise, execute step (10.13). (10.9) FLT i-jm =|(t i -t j -FRQ)(FL i -fl i ) / (fl i -fl j )| , execute step (10.12); (10.10) FLT i-jn =|(t i -t j -FRQ)(FL i -fl i ) / (fl i -fl j )| , execute step (10.12); (10.11) FLT i-jq =|(t i -t j -FRQ)(FL i-jq-1 -FL i-jq ) / (fl i -fl j )| , execute step (10.12); (10.12) Let q = q + 1 and execute step (10.8). (10.13) Let i = i + 1, j = j + 1, and execute step (10.3). (10.14) MT j =FLT in +FLT out +Σ n n=0 FLT i-jn , execute step (10.15); (10.15) End.
8. A device for calculating vertical segmentation and missing time of a track, characterized in that: It includes: The first acquisition module, which is configured to acquire all flight data and classify the flight data according to the 24-bit address code. The second acquisition module, which is configured to acquire the flight data within the specified airspace range. The first calculation module, which is configured to calculate the time when each flight enters and leaves the airspace. The second calculation module, which is configured to calculate whether the time interval between two adjacent messages of each flight meets the established frequency requirement according to the existing surveillance data update frequency. The third calculation module, which is configured to sort out the track points that do not meet the established update frequency requirement and calculate the time lost at these altitude layers passed by adjacent tracks according to the barometric altitudes of adjacent track points.
9. A device for calculating vertical segmentation and missing time of a track, characterized in that: It includes: The track extraction module, which is configured to extract tracks within the specified rectangular airspace range using the 24-bit address code of the aircraft as the unique identifier of the aircraft. The entering and leaving airspace time calculation module, which is configured to calculate the start time when a certain track enters the airspace and the cut-off time when it leaves the airspace. The lost time vector calculation module, which is configured to obtain the time and barometric altitude of each point in the track points, calculate the lost time at each altitude layer from the time when a certain track enters the airspace to the time corresponding to the first track point, and the lost time at each altitude layer from the time corresponding to the last track point in the airspace to the time when it leaves the airspace, and calculate the lost time vector at each altitude layer between the first track point and the last track point of this track within the airspace.