Method and system for automatic calculation of five-sided slant range between aircraft touch-downs

Abstract

The application discloses an aircraft landing interval five-edge slant distance automatic calculation method and system, relates to the field of aviation, and comprises the following steps: acquiring a data set of incoming flights, and sorting the flights in the data set of incoming flights to obtain an ordered flight sequence; for each flight in the ordered flight sequence, determining whether the flight can be used as a preceding aircraft; if yes, acquiring an original track point set of the flight; screening track points meeting preset requirements from the original track point set of the flight, and taking the track points as an effective track point set of the flight; obtaining a runway and an actual landing site of the flight according to the effective track point set of the flight; according to the actual landing site of the flight, finding out a following aircraft of the flight, and determining a spatial distance between the flight and the following aircraft and an interval type between the flight and the following aircraft. The application can accurately calculate the spatial distance between the preceding aircraft and the following aircraft, automatically determine the interval type, and does not need manual intervention.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of aircraft landing interval calculation, and particularly relates to a method and system for automatically calculating five-edge slant distance of aircraft landing interval. BACKGROUND

[0002] Airport approach control is an important part of air traffic management, which needs to accurately calculate the landing interval between the front and rear aircrafts to ensure flight safety. The landing interval includes the same runway and slant distance: the same runway refers to the front and rear aircrafts using the same side and same column of the runway (such as 02L and 02R), and the slant distance refers to the front and rear aircrafts using different sides or different columns of the runway (such as 02L and 20L). The five-edge slant distance refers to the actual spatial distance between the front and rear aircrafts in the approach phase (five-edge approach), which is an important indicator for evaluating whether the control interval meets the safety standard.

[0003] Currently, the main method for calculating the landing interval of aircrafts is the interval calculation method based on the database timestamp, which specifically includes: querying the landing timestamp (ALDT field) of the arrival flight from the flight database; sorting the flights according to the landing timestamp; for each flight, finding the first flight after the landing time as the rear aircraft; calculating the landing time difference between the front and rear aircrafts as the interval; and manually judging whether it is the same runway or slant distance according to the runway number of the front and rear aircrafts.

[0004] However, this method only uses the landing timestamp (ALDT) in the database, which may be the planned time or the system recorded time, and not the actual touch-down time of the aircraft. Due to the possibility of re-approach, circling, etc. during the approach of the aircraft, there may be a large deviation between the actual landing point and the database record, resulting in inaccurate interval calculation. Moreover, this method can only calculate the time interval and cannot calculate the actual spatial distance between the front and rear aircrafts. In the approach control, the spatial distance is more accurate than the time interval in reflecting the safety interval, especially in the case of different runways and different approach paths, two pairs of aircrafts with the same time interval may have a large difference in spatial distance. Moreover, this method simply finds the rear aircraft according to the time sequence, and when multiple front aircrafts land almost simultaneously (different runways), the same rear aircraft may be identified as the rear aircraft of multiple front aircrafts, resulting in data duplication and statistical errors. For example, A aircraft and B aircraft land almost simultaneously, and they may both identify C aircraft as their rear aircraft.

[0005] Moreover, this method needs manual judgment of the same runway and slant distance, and cannot automatically identify according to the runway configuration. The judgment of the same runway and slant distance needs to consider the side (N side or S side) and column (same column or different column) of the runway, and manual judgment is prone to errors and low efficiency. In addition, when the flight path data is interrupted or missing, this method cannot accurately determine the runway and landing point, resulting in calculation failure or inaccurate results. SUMMARY

[0006] To solve the technical problems in the background art, the present application provides an aircraft landing interval five-edge slant distance automatic calculation method and system.

[0007] The aircraft landing interval five-edge slant distance automatic calculation method based on flight path data provided by the present application comprises:

[0008] Obtain a set of arrival flight data, and sort the flights in ascending order according to the landing time stamp of the flights in the set of arrival flight data to obtain an ordered flight sequence;

[0009] For each flight in the ordered flight sequence, determine whether the flight can be used as a preceding aircraft; if so, obtain the original flight path point set of the flight; select the flight path points that meet the preset requirements from the original flight path point set of the flight, and combine the flight path points that meet the preset requirements as the effective flight path point set of the flight; perform runway automatic identification on the flight according to the effective flight path point set of the flight to obtain the runway of the flight; perform landing site identification on the runway of the flight to obtain the actual landing site of the flight; find the following aircraft of the flight from the ordered flight sequence according to the actual landing site of the flight, and calculate the spatial distance between the flight and the following aircraft;

[0010] Determine the interval type of the spatial distance between the flight and the following aircraft according to the configuration information of the runways of the flight and the following aircraft.

[0011] Preferably, the original flight path point set of the flight is obtained, specifically comprising:

[0012] Obtain all flight path points of the flight from the flight path database according to the flight number and flight date of the flight, and combine them into the original flight path point set.

[0013] Preferably, the runway of the flight is obtained by performing runway automatic identification on the flight according to the effective flight path point set of the flight, specifically comprising:

[0014] Obtain the configuration information of all preset runways of the flight; initialize the count of flight path points in each preset runway;

[0015] For each flight path point in the effective flight path point set of the flight, use the point-in-polygon algorithm to determine the preset runway on which each flight path point falls;

[0016] When traversing the effective flight path point set of the flight, select the preset runway with the most flight path points as the runway of the flight.

[0017] Preferably, the actual landing site of the flight is obtained by performing landing site identification on the runway of the flight, specifically comprising:

[0018] acquire the coordinate of the theoretical landing point of the preset runway configuration of the flight;

[0019] for each track point in the set of valid track points of the flight, calculate the distance between each track point and the theoretical landing point using the Haversine formula:

[0020] determine whether the distance between each track point and the theoretical landing point is less than the first preset distance; if yes, take the track point as the nearest point; if no, end the process;

[0021] when the set of valid track points of the flight is traversed, determine whether there is a nearest point; if yes, select the track point corresponding to the minimum value from the nearest point as the actual landing point of the flight.

[0022] Preferably, according to the actual landing point of the flight, the aircraft behind the flight is found from the ordered flight sequence, specifically comprising:

[0023] According to the timestamp of the actual landing point of the flight, find the aircraft with a landing timestamp greater than the timestamp of the actual landing point of the flight from the ordered flight sequence as the candidate aircraft behind the flight;

[0024] construct a candidate aircraft list for the flight and initialize it;

[0025] for each candidate aircraft of the flight, determine whether the candidate aircraft has been used; if yes, end the process; if no, find the coordinate point of the candidate aircraft at the landing point timestamp of the flight;

[0026] According to the coordinate of the actual landing point of the flight, determine whether the coordinate point of the candidate aircraft is within the approach route area: if no, end the process; if yes, calculate the spatial distance between the candidate aircraft and the aircraft in front and behind the flight;

[0027] determine whether the spatial distance between the candidate aircraft and the aircraft in front and behind the flight is less than or equal to the second preset distance; if no, end the process; if yes, add the candidate aircraft and the spatial distance between the aircraft in front and behind to the candidate list;

[0028] When the candidate aircraft of the flight is traversed, determine whether the candidate aircraft list is empty; if yes, determine that the flight has no aircraft behind; if no, find the candidate aircraft corresponding to the minimum value of the spatial distance between the aircraft in front and behind from the candidate aircraft list as the final aircraft behind the flight, and mark the aircraft behind the flight as used.

[0029] Preferably, determining whether the coordinate point of the candidate aircraft is within the approach route area specifically comprises:

[0030] acquire the configuration information of the runway of the candidate aircraft;

[0031] According to the configuration information of the runway of the candidate follower aircraft, a corresponding approach route area polygon is obtained; and a point-in-polygon judgment algorithm is used to judge whether the coordinate point of the candidate follower aircraft is in the approach route area.

[0032] Preferably, when each flight in the ordered flight sequence is traversed, a five-edge slant range result table is generated; wherein the five-edge slant range result table includes leader aircraft information, follower aircraft information, interval information, and a validity flag.

[0033] Preferably, before judging whether each flight in the ordered flight sequence can serve as a leader aircraft, the method further comprises: constructing a used follower aircraft set and initializing for the ordered flight sequence;

[0034] After judging the interval type of the spatial distance between the flight and its follower aircraft, the method further comprises: adding the flight, the follower aircraft of the flight, the spatial distance therebetween, and the interval type to the used follower aircraft set.

[0035] Preferably, before marking the follower aircraft of the flight as used, the method further comprises: double-checking the follower aircraft of the flight;

[0036] The double-checking comprises: a first check and a second check.

[0037] The first check is to check whether the follower aircraft of the flight is already in the used follower aircraft set; if yes, the check fails; if no, the check passes.

[0038] The second check is to check whether the follower aircraft of the flight is already used by other flights in a preset database; if yes, the check fails; if no, the check passes.

[0039] If either of the checks fails, the double-checking fails; if both the first check and the second check pass, the double-checking succeeds.

[0040] Preferably, the interval type comprises a same runway and a slant range.

[0041] The judgment of the interval type of the spatial distance between the flight and its follower aircraft according to the configuration information of the runways of the flight and its follower aircraft comprises:

[0042] According to the configuration information of the runway of the flight, the side and the column to which the runway threshold of the flight belongs are obtained; wherein the side is divided into an N side and an S side.

[0043] According to the configuration information of the runway of the follower aircraft, the side and the column to which the runway threshold of the follower aircraft belongs are obtained.

[0044] determine the interval type of the spatial distance between the flight and the subsequent flight according to the side and column to which the runway head of the flight and the subsequent flight belongs;

[0045] When the runway heads of the flight and the subsequent flight are on the same side and belong to the same column, it is determined that the interval type of the spatial distance between the flight and the subsequent flight is the same runway.

[0046] When the runway heads of the flight and the subsequent flight are not on the same side, or are on the same side but belong to different columns, it is determined that the interval type of the spatial distance between the flight and the subsequent flight is the oblique distance.

[0047] In a second aspect, the present application further provides an aircraft landing interval five-edge oblique distance automatic calculation system based on flight path data, comprising:

[0048] An acquisition unit is configured to acquire a set of inbound flight data and a set of original flight path points of each flight in the set of inbound flight data.

[0049] A processing unit is configured to sort the flights in the set of inbound flight data in ascending order according to the landing time stamp of each flight, and obtain an ordered flight sequence.

[0050] For each flight in the ordered flight sequence, it is determined whether the flight can be used as a preceding flight. If yes, the flight path points meeting the preset requirements are selected from the set of original flight path points of the flight, and the flight path points meeting the preset requirements are combined as the set of effective flight path points of the flight. The runway of the flight is automatically identified according to the set of effective flight path points of the flight, and the actual landing site of the flight is identified according to the runway of the flight. The subsequent flight of the flight is found from the ordered flight sequence according to the actual landing site of the flight, and the spatial distance between the flight and the subsequent flight is calculated. The interval type of the spatial distance between the flight and the subsequent flight is determined according to the configuration information of the runway between the flight and the subsequent flight.

[0051] In the present application, the aircraft landing interval five-edge oblique distance automatic calculation method and system proposed in the present application can select the flight path points meeting the preset requirements from the set of original flight path points of the flight that can be used as a preceding flight in the ordered flight sequence, and combine the flight path points meeting the preset requirements as the set of effective flight path points of the flight. The runway of the flight is automatically identified according to the set of effective flight path points of the flight, and the actual landing site of the flight is identified according to the runway of the flight. The subsequent flight of the flight is found from the ordered flight sequence according to the actual landing site of the flight, and the spatial distance between the flight and the subsequent flight is calculated. The interval type of the spatial distance between the flight and the subsequent flight is determined according to the configuration information of the runway between the flight and the subsequent flight. BRIEF DESCRIPTION OF DRAWINGS

[0052] Figure 1 A flowchart of an aircraft landing interval five-side slant distance automatic calculation method based on track data in an embodiment of the present application. DETAILED DESCRIPTION

[0053] It should be noted that the embodiments in the present application and the features in the embodiments can be combined with each other without conflict. The present application will be described in detail below with reference to the accompanying drawings and in combination with the embodiments.

[0054] In a first aspect, with reference to Figure 1 The present application provides an aircraft landing interval five-side slant distance automatic calculation method based on track data, which comprises the following steps:

[0055] Obtaining a set of arrival flight data;

[0056] According to the landing time stamp of the flight in the set of arrival flight data, the flights are sorted in ascending order to obtain an ordered flight sequence;

[0057] For each flight in the ordered flight sequence, it is determined whether the flight can be used as a preceding aircraft. If not, the process is ended. If yes, the original track point set of the flight is obtained;

[0058] From the original track point set of the flight, track points meeting the preset requirements are selected, and the selected track points are combined as the effective track point set of the flight;

[0059] According to the effective track point set of the flight, runway automatic identification is performed on the flight to obtain the runway of the flight;

[0060] The landing point of the runway of the flight is identified to obtain the actual landing point of the flight;

[0061] According to the actual landing point of the flight, the following aircraft of the flight is found from the ordered flight sequence, and the spatial distance between the flight and the following aircraft is calculated;

[0062] According to the configuration information of the runway between the flight and the following aircraft, the interval type of the spatial distance between the flight and the following aircraft is determined.

[0063] The application can screen the track points meeting the preset requirements from the original track point set of the flight as the preceding aircraft in the ordered flight sequence, combine the track points meeting the preset requirements as the effective track point set of the flight, perform runway automatic identification on the flight according to the effective track point set of the flight, obtain the runway of the flight, perform landing site identification on the runway of the flight, automatically identify the actual landing site of the flight, accurately calculate the spatial distance between the preceding and following aircrafts based on the effective track point sets of the preceding and following aircrafts, and provide a more accurate interval evaluation index than the time interval.

[0064] In the embodiment, the inbound flight data set F includes a plurality of flights, and each flight includes a flight number, a date, a landing timestamp, an airline, a runway, a speed, an aircraft type, and a bitstream level.

[0065] In the embodiment, before determining whether each flight in the ordered flight sequence can serve as the preceding aircraft, the ordered flight sequence is constructed with a used following aircraft set and initialized to avoid repeated identification.

[0066] After determining the interval type of the spatial distance between the flight and the following aircraft, the flight, the following aircraft, the spatial distance therebetween, and the interval type are added to the used following aircraft set.

[0067] In the embodiment, the used following aircraft set is maintained, and each pair of preceding and following aircrafts determined is recorded in the used following aircraft set. The used following aircraft set is used to check whether the following aircraft has been used when a subsequent following aircraft is determined, so that the following aircraft can only be used by one preceding aircraft, and the situation that multiple preceding aircrafts simultaneously identify the same following aircraft is avoided.

[0068] In the embodiment, the ordered flight sequence is represented as ;

[0069] In the embodiment, after the track data of each flight is stored in the track database, the track data of each flight is subjected to data integrity verification, and the verification is repeated at most 5 times with an interval of 1 second each time, so as to ensure that the track data has been successfully written into the track database. If the verification fails, a warning is recorded but the subsequent analysis is continued, so as to avoid infinite waiting.

[0070] In the embodiment, the used following aircraft set of the ordered flight sequence is initialized, and the initialization specifically includes:

[0071] The used following aircraft set of the ordered flight sequence is subjected to idempotency processing to delete the existing landing interval record corresponding to the ordered flight sequence currently processed, so as to ensure the consistency of data when the data is reanalyzed.

[0072] In the embodiment, the original track point set of the flight is acquired, specifically including:

[0073] According to the flight number and flight date of the flight, all track points of the flight are acquired from the track database, and are combined into the original track point set.

[0074] It should be understood that the original coordinates of the track points are WGS84 coordinates, which need to be converted into GCJ02 coordinates (Mars coordinates used by Gaode Map) to match the map system.

[0075] The preset requirement in the embodiment is that the height is less than or equal to 600 meters. It should be understood that the track points with the height less than or equal to 600 meters represent the trajectory of the aircraft in the approach phase (five-edge approach).

[0076] In a further embodiment, after the track points meeting the preset requirement are filtered out from the original track point set of the flight, the method further includes:

[0077] Attribute data of each track point meeting the preset requirement is extracted; wherein the attribute data includes longitude and latitude coordinates (lon, lat), a timestamp ts, a height alt_m, a speed speed_kmh and a wake turbulence category wtc.

[0078] It should be understood that the height alt_m is converted from the MEASURED FLIGHT LEVEL field (unit: feet to meters, conversion factor: 1 foot = 0.3048 meters). The timestamp ts is calculated based on the TOD field (offset seconds of UTC time of the day): the date is extracted from the adc_insert_time or adc_insert_utc_time field, the timestamp of UTC 00:00:00 of the day is calculated, the TOD seconds are added, and 8 hours (offset of East Eight Zone) are added to obtain the timestamp of East Eight Zone. The speed is acquired by using an intelligent searching mechanism: the INDICATED AIRSPEED (indicated airspeed, unit: knots) of the current point is used preferentially, and if it is empty, the AIRSPEED (unit: value > 500 is considered as km / h, otherwise it is considered as knots and converted into km / h) is used; if the speed of the current point is empty, at most 3 subsequent points are searched backward to try to acquire the speed value, ensuring the integrity of the speed data.

[0079] In the embodiment, the runway of the flight is obtained by performing runway automatic identification on the flight according to the valid track point set of the flight, specifically including:

[0080] Configuration information of all preset runways of the flight is acquired, and the count of the track points in each preset runway is initialized;

[0081] For each track point in the valid track point set of the flight, the point-in-polygon algorithm is used to determine the preset runway where each track point falls;

[0082] When traversing the valid track point set of the flight, the preset runway with the most track points is selected as the runway of the flight.

[0083] In this embodiment, by counting the number of track points in the polygon range of each preset runway, the preset runway with the most points is identified as the runway used by the aircraft, which can handle abnormal track data and improve the robustness of the calculation.

[0084] The counting of track points in each runway is initialized, specifically including:

[0085] An empty counter dictionary is established to record the number of track points matched by each runway. Each preset runway in the input preset runway configuration is traversed, and an entry is created in the counter dictionary for each preset runway, and the initial count is set to 0.

[0086] For each track point in the valid track point set of the flight, the point-in-polygon algorithm is used to determine the preset runway where each track point falls, specifically including:

[0087] Read each point in the valid track point set in order , get the latitude and longitude coordinates of the point and form a coordinate point;

[0088] For the current coordinate point, each preset runway is checked in turn:

[0089] For each preset runway, get the boundary polygon coordinates of the runway, and use the "point-in-polygon" algorithm to determine whether the current coordinate point is within the polygon range of the runway. If so, increment the count of the preset runway in the counter by 1, and stop checking the remaining runways for the current point and start processing the next track point (because a point theoretically belongs to only one runway). If not, continue to check the next runway until the current coordinate point is determined to be within the polygon range of a preset runway.

[0090] The preset runway with the most track points is selected as the runway of the flight, specifically including:

[0091] After traversing and judging all track points in the valid track point set, the runway counter is checked. Find the preset runway with the largest count in the counter dictionary. The preset runway is used as the runway of the flight.

[0092] In this embodiment, the runway of the flight is identified for the landing site, and the actual landing site of the flight is obtained, specifically including:

[0093] obtain the coordinate of the theoretical landing point of the preset runway configuration of the flight;

[0094] For each track point in the set of valid track points of the flight, the distance between each track point and the theoretical landing point is calculated using the Haversine formula:

[0095] determine whether the distance between each track point and the theoretical landing point is less than the first preset distance; if yes, the track point is taken as the nearest point; if no, the process ends and the calculation and determination of the next track point begins;

[0096] When the set of valid track points of the flight is traversed, it is determined whether there is a nearest point; if yes, the track point corresponding to the minimum value is selected from the nearest point, i.e., the track point closest to the runway landing point, as the actual landing point of the flight, and the attribute information of the actual landing point is recorded; if no, the information of the track point closest to the theoretical landing point is recorded as diagnostic information, and an empty value is returned for subsequent analysis of the case where the landing point is not found.

[0097] The present embodiment can reliably identify the position of the actual contact of the aircraft with the runway through reasonable radius setting and coordinate system processing.

[0098] In one specific embodiment, the first preset distance is 180 meters to avoid missing the actual landing point due to GPS error.

[0099] It should be noted that the coordinate of the theoretical landing point of the runway configuration is in the WGS84 coordinate system, and the track points in the present embodiment have been converted to the GCJ02 coordinate system, so the coordinate systems need to be unified for distance calculation in the same coordinate system.

[0100] wherein the distance between each track point and the theoretical landing point is specifically:

[0101] distance = Haversine Distance([lon, lat], runway_landing_point); wherein distance is the distance between the track point and the theoretical landing point, runway_landing_point is the theoretical landing point, (lon, lat) is the latitude and longitude coordinate of the flight, and Haversine Distance is the Haversine algorithm.

[0102] It should be noted that the Haversine formula is used to calculate the great circle distance (spherical distance) between two points on the surface of the Earth, which takes into account the curvature of the Earth and is suitable for calculating the actual spatial distance between aircraft.

[0103] The attribute information of the actual landing point of the flight as the front aircraft in the embodiment includes coordinates (lon_front, lat_front), a timestamp ts_front, a speed speed_kmh_front, and a wake turbulence category wtc_front.

[0104] The embodiment automatically identifies the closest track point as the actual landing point within a range of 180 meters by analyzing track points with an altitude of less than or equal to 600 meters and combining the landing point coordinates of the runway configuration, and obtains an accurate timestamp and coordinates of the actual landing point.

[0105] In the embodiment, the rear aircraft of the flight is found from the ordered flight sequence according to the actual landing point of the flight, specifically including:

[0106] According to the timestamp of the actual landing point of the flight, a flight with a landing timestamp greater than the timestamp of the actual landing point of the flight is found from the ordered flight sequence as a candidate rear aircraft of the flight.

[0107] A candidate rear aircraft list is constructed for the flight and initialized.

[0108] For each candidate rear aircraft of the flight, it is determined whether the candidate rear aircraft has been used, that is, whether it has been in the used rear aircraft set; if yes, the process ends and the candidate rear aircraft is skipped; if not, the coordinate point of the candidate rear aircraft at the landing point timestamp of the flight is found.

[0109] According to the coordinates of the actual landing point of the flight, it is determined whether the coordinate point of the candidate rear aircraft is within the approach route area: if not, the process ends and the candidate rear aircraft is skipped; if yes, the front-rear aircraft spatial distance of the candidate rear aircraft is calculated.

[0110] It is determined whether the front-rear aircraft spatial distance of the candidate rear aircraft is less than or equal to a second preset distance; if not, the process ends and the candidate rear aircraft is skipped; if yes, the candidate rear aircraft and the front-rear aircraft spatial distance are added to the candidate list.

[0111] When the candidate rear aircraft of the flight is traversed, it is determined whether the candidate rear aircraft list is empty; if yes, it is determined that the flight has no rear aircraft, the process ends, and the next flight begins; if not, the candidate rear aircraft corresponding to the minimum value of the front-rear aircraft spatial distance is found from the candidate rear aircraft list as the final rear aircraft of the flight, and the rear aircraft of the flight is marked as used.

[0112] It needs to be understood that in the process of finding the coordinate point of the candidate post-machine at the time stamp of the landing point of the flight, first, the track point with the closest time stamp to the actual landing point of the flight (pre-machine) is found from the set of valid track points of the candidate post-machine; if no exact matching track point is found, the earliest available track point of the post-machine after the time stamp of the actual landing point (the time stamp must be greater than or equal to the time stamp of the actual landing point) is found, ensuring that even if there is a time interval in the track data, a suitable coordinate point can be found; record its coordinates, speed and wake level. The speed acquisition also uses an intelligent search mechanism: preferentially use the indicated airspeed I, if empty, use the ground speed or other speed, and if the current point speed is empty, search back for at most 3 subsequent points.

[0113] The embodiment is thus arranged, which can handle abnormal situations such as interruption and absence of track data, and effectively improves the calculation robustness.

[0114] The candidate post-machine list in the embodiment is arranged in ascending order according to the landing time stamp of the candidate post-machine.

[0115] The second preset distance in the embodiment is 20 km.

[0116] In a further embodiment, the judgment of whether the coordinate point of the candidate post-machine is in the approach route area specifically includes:

[0117] The configuration information of the runway of the candidate post-machine is obtained.

[0118] According to the configuration information of the runway of the candidate post-machine, the corresponding approach route area polygon is obtained; the point-in-polygon judgment algorithm is used to judge whether the coordinate point of the candidate post-machine is in the approach route area.

[0119] The embodiment judges the approach route area by the track point, which can handle abnormal situations such as interruption and absence of track data, and improves the calculation robustness.

[0120] It needs to be understood that the configuration information of the runway of the candidate post-machine in the embodiment is preferentially determined from the track data of the candidate post-machine, and if it cannot be determined, the ARWY (runway) field in the preset database is used as a backup.

[0121] In the embodiment, the spatial distance between the flight and the post-machine thereof is calculated, specifically including:

[0122] The Haversine formula is used to calculate the distance between the actual landing point of the flight and the coordinate point of the post-machine of the flight at the time stamp of the actual landing point of the flight.

[0123] wherein distance_km = Haversine Distance([lon_front, lat_front], [lon_rear_at_time, lat_rear_at_time]);

[0124] wherein distance_km represents the distance between the actual landing point of the flight and the coordinate point of the rear aircraft at the time stamp of the actual landing point of the flight, [lon_front, lat_front] is the coordinate of the actual landing point of the flight, and [lon_rear_at_time, lat_rear_at_time] is the coordinate point of the rear aircraft of the flight at the time stamp of the actual landing point of the flight.

[0125] In a further embodiment, before marking the rear aircraft of the flight as used, the rear aircraft of the flight is double-checked.

[0126] The double checking in the embodiment includes a first checking and a second checking.

[0127] The first checking is to check whether the rear aircraft of the flight is an existing rear aircraft in the used rear aircraft set; if yes, it is determined that the checking fails; if no, it is determined that the checking passes.

[0128] The second checking is to check whether the rear aircraft of the flight is a rear aircraft used by other flights in the preset database; if yes, it is determined that the checking fails; if no, it is determined that the checking passes.

[0129] If any of the checkings fails, it is determined that the double checking fails; if both the first checking and the second checking pass, it is determined that the double checking succeeds.

[0130] When the double checking succeeds, the rear aircraft of the flight is marked as used and added to the used rear aircraft set.

[0131] When multiple front aircrafts land almost at the same time, the first front aircraft finds the rear aircraft and marks it as used first, and the subsequent front aircrafts are automatically excluded from the rear aircraft when searching (whether in the memory set or in the record in the database), so as to avoid repeated identification. This mechanism can handle the repeated identification problem across batches and time periods, ensure that a rear aircraft can only be used by one front aircraft, and ensure data accuracy.

[0132] The interval types in the embodiment include the same runway and the slant distance.

[0133] In the embodiment, according to the configuration information of the runway of the flight and the rear aircraft thereof, the interval type of the spatial distance between the flight and the rear aircraft thereof is determined, specifically including:

[0134] Based on the runway configuration information for this flight, the side and column to which the runway head belongs are determined; the side is divided into the N side and the S side.

[0135] Based on the runway configuration information of the aircraft following the flight, the side and column to which the runway head of the aircraft following the flight belongs can be obtained;

[0136] Determine the type of spatial interval between the flight and the aircraft following it based on the side and column of the runway to which the flight and the aircraft following it belong.

[0137] Among them, if the runways of the flight and the following aircraft are on the same side and of the same type, the spatial distance interval between the flight and the following aircraft is determined to be the same runway.

[0138] If the runways of the flight and the following aircraft are not on the same side, or on the same side but in different columns, the spatial distance between the flight and the following aircraft is determined to be a slant distance.

[0139] In this embodiment, when traversing each flight in the ordered flight sequence, a pentagonal slant distance result table is generated; wherein, the pentagonal slant distance result table includes preceding aircraft information, following aircraft information, interval information, and valid markers.

[0140] The preceding aircraft information includes flight number, airline, runway, actual landing timestamp, speed, aircraft type, and wake turbulence level.

[0141] The information for the following aircraft includes flight number, airline, runway, actual landing timestamp, speed (using the speed of the following aircraft at the landing point timestamp of the preceding aircraft), aircraft type, and wake turbulence level.

[0142] The interval information includes spatial distance (unit: kilometers) and interval type (1=same track, 2=slant distance).

[0143] The validity flag includes the valid field, where 1 = valid, indicating that there is a follow-up mechanism and the judgment is normal; 0 = invalid, indicating that there is no follow-up mechanism or the judgment failed.

[0144] The pentagonal slope distance result table in this embodiment also includes notes information. These notes information includes a "notes" field. The "notes" field is in JSON format and is used to record debugging information such as the reasons for errors and failures during the judgment process.

[0145] In this case, the runway to which the following aircraft belongs is determined by querying the flight path data of the following aircraft and counting the number of flight path points falling in each runway area. If it cannot be determined by the flight path data, the runway field in the database is used as a backup.

[0146] This embodiment retrieves the inbound flight data set from a preset database and stores the generated pentagonal slant distance result table into the database.

[0147] In a second aspect, the present application provides an aircraft landing interval five-side slant distance automatic calculation system based on flight path data, comprising:

[0148] An acquisition unit is configured to acquire a set of inbound flight data and a set of original flight path points of each flight in the set of inbound flight data;

[0149] A processing unit is configured to sort the flights in the set of inbound flight data in ascending order according to the landing time stamp of each flight, to obtain an ordered flight sequence;

[0150] For each flight in the ordered flight sequence, it is determined whether the flight can be used as a preceding aircraft. If yes, the flight path points meeting the preset requirements are filtered from the set of original flight path points of the flight, and the flight path points meeting the preset requirements are combined as the set of effective flight path points of the flight. The runway of the flight is automatically identified according to the set of effective flight path points of the flight, to obtain the runway of the flight. The landing point of the runway of the flight is identified, to obtain the actual landing point of the flight. The flight following the flight in the ordered flight sequence is found according to the actual landing point of the flight, and the spatial distance between the flight and the flight following the flight is calculated. The spatial distance between the flight and the flight following the flight is determined to be of an interval type according to the configuration information of the runway between the flight and the flight following the flight.

[0151] In this embodiment, a set of inbound flight data is acquired.

[0152] In this embodiment, the set of original flight path points of the flight is acquired, specifically including:

[0153] The flight number and flight date of the flight are used to acquire all flight path points of the flight from a flight path database, and the flight path points are combined as the set of original flight path points.

[0154] In this embodiment, the runway of the flight is automatically identified according to the set of effective flight path points of the flight, specifically including:

[0155] The configuration information of all preset runways of the flight is acquired, and the count of flight path points in each preset runway is initialized.

[0156] For each flight path point in the set of effective flight path points of the flight, the point-in-polygon algorithm is used to determine the preset runway on which each flight path point falls.

[0157] When the set of effective flight path points of the flight is traversed, the preset runway with the largest number of flight path points is selected as the runway of the flight.

[0158] In this embodiment, the actual landing point of the flight is identified according to the runway of the flight, specifically including:

[0159] The coordinates of the theoretical landing point configured by the preset runway of the flight are acquired.

[0160] For each track point in the valid track point set of the flight, the distance between each track point and the theoretical landing site is calculated using the Haversine formula:

[0161] It is judged whether the distance between each track point and the theoretical landing site is less than the first preset distance; if yes, the track point is taken as the nearest point; if no, the process ends;

[0162] When the valid track point set of the flight is traversed, it is judged whether there is a nearest point; if yes, the track point corresponding to the minimum value is selected from the nearest point as the actual landing site of the flight.

[0163] Among them, according to the actual landing site of the flight, the aircraft behind the flight is found from the ordered flight sequence, which specifically includes:

[0164] According to the timestamp of the actual landing site of the flight, the flight whose landing timestamp is greater than the timestamp of the actual landing site of the flight is found from the ordered flight sequence as the candidate aircraft behind the flight;

[0165] A candidate aircraft list is constructed for the flight and initialized;

[0166] For each candidate aircraft of the flight, it is judged whether the candidate aircraft has been used; if yes, the process ends; if no, the coordinate point of the candidate aircraft at the landing site timestamp of the flight is found;

[0167] According to the coordinates of the actual landing site of the flight, it is judged whether the coordinate point of the candidate aircraft is within the approach route area; if no, the process ends; if yes, the front and rear aircraft spatial distance of the candidate aircraft and the flight is calculated;

[0168] It is judged whether the front and rear aircraft spatial distance of the candidate aircraft and the flight is less than or equal to the second preset distance; if no, the process ends; if yes, the candidate aircraft and the front and rear aircraft spatial distance are added to the candidate list;

[0169] When the candidate aircraft of the flight is traversed, it is judged whether the candidate aircraft list is empty; if yes, it is determined that the flight has no aircraft, and the process ends; if no, the candidate aircraft corresponding to the minimum value of the front and rear aircraft spatial distance is found from the candidate aircraft list as the final aircraft behind the flight, and the aircraft behind the flight is marked as used.

[0170] The process of judging whether the coordinate point of the candidate aircraft is within the approach route area in the embodiment specifically includes:

[0171] The configuration information of the runway of the candidate aircraft is obtained;

[0172] According to the configuration information of the runway of the candidate follower aircraft, a corresponding approach route area polygon is obtained; and a point-in-polygon judgment algorithm is used to judge whether the coordinate point of the candidate follower aircraft is in the approach route area.

[0173] The calculation process of the spatial distance between the flight and the follower aircraft in the embodiment specifically includes:

[0174] The Haversine formula is used to calculate the distance between the actual landing site of the flight and the coordinate point of the follower aircraft at the timestamp of the actual landing site of the flight.

[0175] The interval types in the embodiment include the same runway and the diagonal distance;

[0176] The process of judging the interval type of the spatial distance between the flight and the follower aircraft in the embodiment specifically includes:

[0177] According to the configuration information of the runway of the flight, the side and the column to which the runway threshold of the flight belongs are obtained; wherein, the side is divided into N side and S side;

[0178] According to the configuration information of the runway of the follower aircraft of the flight, the side and the column to which the runway threshold of the follower aircraft of the flight belongs are obtained;

[0179] According to the side and the column to which the runway threshold of the flight and the follower aircraft belongs, the interval type of the spatial distance between the flight and the follower aircraft is judged;

[0180] When the runway thresholds of the flight and the follower aircraft are on the same side and belong to the same column, it is determined that the interval type of the spatial distance between the flight and the follower aircraft is the same runway;

[0181] When the runway thresholds of the flight and the follower aircraft are not on the same side, or are on the same side but belong to different columns, it is determined that the interval type of the spatial distance between the flight and the follower aircraft is the diagonal distance.

[0182] The processing module in the embodiment is further used to build a used follower aircraft set and initialize for the ordered flight sequence before judging whether the flight in the ordered flight sequence can be used as a leader aircraft;

[0183] After judging the interval type of the spatial distance between the flight and the follower aircraft, the flight, the follower aircraft of the flight, the spatial distance therebetween and the interval type are added to the used follower aircraft set.

[0184] The processing module in the embodiment is further used to double-check the follower aircraft of the flight before marking the follower aircraft of the flight as used;

[0185] The double-checking includes a first double-checking and a second double-checking.

[0186] The first re-inspection is to check whether the back aircraft of the flight is a back aircraft already existing in the used back aircraft set; if yes, it is determined that the re-inspection fails; if no, it is determined that the re-inspection passes;

[0187] The second re-inspection is to check whether the back aircraft of the flight is a back aircraft already used by other flights in the database; if yes, it is determined that the re-inspection fails; if no, it is determined that the re-inspection passes;

[0188] If any of the inspections fails, it is determined that the double re-inspection fails; if the first re-inspection and the second re-inspection both pass, it is determined that the double re-inspection succeeds.

[0189] The above merely describes the preferred embodiments of the present application, but the protection scope of the present application is not limited to this. Any person skilled in the art, according to the technical solution and the inventive concept of the present application, makes equivalent replacement or change within the technical range disclosed by the present application, which should be covered in the protection scope of the present application.

Claims

1. A method for automatic calculation of aircraft landing interval five-sided slant range based on track data, characterized in that, include: Obtain the inbound flight data set, and sort the flights in ascending order according to their arrival timestamps to obtain an ordered flight sequence; for each flight in the ordered flight sequence, determine whether the flight can be used as the preceding flight; If so, obtain the original set of waypoints for the flight; select waypoints that meet the preset requirements from the original set of waypoints for the flight, and combine the waypoints that meet the preset requirements as the set of valid waypoints for the flight. Obtain the configuration information of all preset runways for the flight; initialize the count of track points in each preset runway; for each track point in the set of valid track points for the flight, use the point-in-polygon judgment algorithm to determine the preset runway where each track point lands; when traversing the set of valid track points for the flight, select the preset runway with the most track points as the runway for the flight. Obtain the coordinates of the theoretical landing point of the flight's preset runway configuration; for each track point in the set of valid track points for the flight, use the Haversine formula to calculate the distance between each track point and the theoretical landing point: determine whether the distance between each track point and the theoretical landing point is less than a first preset distance; if so, take the track point as the closest point; if not, the process ends; after traversing the set of valid track points for the flight, determine whether there is a closest point; If so, then select the track point corresponding to the minimum value from the nearest points as the actual landing point of the flight; Based on the actual landing timestamp of the flight, find flights with landing timestamps greater than the actual landing timestamp of the flight from the ordered flight sequence as candidate successors; construct and initialize a candidate successor list for the flight; for each candidate successor, determine whether the candidate successor has been used; if yes, the process ends; if no, find the coordinates of the candidate successor at the landing timestamp of the flight; based on the coordinates of the actual landing point of the flight, determine whether the coordinates of the candidate successor are within the approach route area; if no, the process ends; if yes, calculate the spatial distance between the candidate successor and the preceding and following aircraft of the flight; determine whether the spatial distance between the candidate successor and the preceding and following aircraft of the flight is less than or equal to a second preset distance; if no, the process ends. If so, add the candidate rear vehicle and the spatial distance between the front and rear vehicles to the candidate list; When iterating through the candidate aircraft following a flight, check if the list of candidate aircraft following a flight is empty; if it is, determine that the flight has no aircraft following a flight. If not, the candidate aircraft corresponding to the minimum spatial distance between the preceding and following aircraft is found from the candidate aircraft list and used as the final aircraft for the flight. The aircraft for the flight is marked as used, and the spatial distance between the flight and its aircraft is calculated. Based on the runway configuration information of the flight and its aircraft, the interval type of the spatial distance between the flight and its aircraft is determined.

2. The automatic calculation method for the pentagonal slant distance of aircraft landing interval based on flight track data according to claim 1, characterized in that, Obtain the original set of waypoints for this flight, specifically including: Based on the flight number and flight date, retrieve all track points for the flight from the track database and combine them into an original track point set.

3. The automatic calculation method for the pentagonal slant distance of aircraft landing interval based on flight track data according to claim 1, characterized in that, Determining whether the coordinates of the candidate aircraft are within the approach route area specifically includes: Obtain the runway configuration information for the candidate aircraft; Based on the runway configuration information of the candidate aircraft, the corresponding approach route area polygon is obtained; using a point-within-polygon judgment algorithm, it is determined whether the coordinates of the candidate aircraft are within the approach route area.

4. The automatic calculation method for the pentagonal slant distance of aircraft landing interval based on flight track data according to claim 1, characterized in that, Interval types include same-track and slant-distance; Among these, based on the runway configuration information of the flight and its following aircraft, the type of spatial interval between the flight and its following aircraft is determined, specifically including: Based on the runway configuration information for this flight, the side and column to which the runway head belongs are determined; the side is divided into the N side and the S side. Based on the runway configuration information of the flight and its following aircraft, the side and column to which the runway head of the flight and its following aircraft belong can be obtained; Determine the type of spatial interval between the flight and the aircraft following it based on the side and column of the runway head. When the runway noses of the flight and the following aircraft are on the same side and of the same type, the spatial distance interval between the flight and the following aircraft is determined to be of the same runway type. If the runway noses of the flight and the following aircraft are not on the same side, or are on the same side but belong to different columns, then the spatial distance between the flight and the following aircraft is determined to be slant distance.

5. The automatic calculation method for the pentagonal slant distance of aircraft landing interval based on flight track data according to claim 1, characterized in that, Before determining whether a flight can be used as a preceding flight for each flight in an ordered flight sequence, the process also includes: constructing and initializing a set of used successor flights for the ordered flight sequence; After determining the interval type of the spatial distance between the flight and its successor aircraft, the process also includes: adding the flight, the successor aircraft of the flight, the spatial distance between them, and the interval type to the set of used successor aircraft.

6. The automatic calculation method for the pentagonal slant distance of aircraft landing interval based on flight track data according to claim 5, characterized in that, Before marking the rear aircraft of the flight as used, the process also includes: double-checking the rear aircraft of the flight; The double check includes: a first check and a second check; The first check is to check whether the following aircraft of the flight is an existing following aircraft in the used following aircraft group; if yes, the check is deemed to have failed; if no, the check is deemed to have passed. The second check is to check whether the aircraft following the flight is a pre-set aircraft in the database that has already been used by other flights; if so, the check is deemed to have failed; if not, the check is deemed to have passed. If either of the two tests fails, the double test is considered to have failed; if both the first and second tests pass, the double test is considered to have succeeded.

7. The automatic calculation method for the pentagonal slant distance of aircraft landing interval based on flight track data according to claim 1, characterized in that, When traversing each flight in the ordered flight sequence, a pentagonal slant distance result table is generated; the pentagonal slant distance result table includes preceding flight information, following flight information, interval information, and valid markers.

8. An automatic calculation system for the pentagonal slant distance of aircraft landing interval based on flight track data, applied to the automatic calculation method for the pentagonal slant distance of aircraft landing interval based on flight track data as described in any one of claims 1-7, characterized in that, include: The acquisition unit is used to acquire the inbound flight data set and the original track point set of each flight in the inbound flight data set; The processing unit is used to sort the flights in ascending order based on the arrival timestamps of the flights in the inbound flight data set to obtain an ordered flight sequence; For each flight in an ordered flight sequence, determine whether the flight can be used as the preceding flight; If so, select the set of set points that meet the preset requirements from the original set of set points for the flight, and combine the set points that meet the preset requirements as the set of valid set points for the flight; based on the set of valid set points for the flight, automatically identify the runway for the flight to obtain the runway; identify the landing point for the flight to obtain the actual landing point; based on the actual landing point, find the aircraft following the flight from the ordered flight sequence and calculate the spatial distance between the flight and the aircraft following it; based on the runway configuration information between the flight and the aircraft following it, determine the interval type of the spatial distance between the flight and the aircraft following it.

Images