A Method and System for Calculating Route Congestion Index Based on Route Segment Spatial Analysis

Abstract

This invention discloses a method and system for calculating airway congestion index based on airway segment spatial analysis. The method includes: S1, dividing the airway into three-dimensional segment calculation airspaces based on airway segments in the airway GIS spatial model; S2, collecting historical ADS-B data by hourly time series in the three-dimensional segment calculation airspace, and collecting ADS-B data from the past N2 days of the study day by hour in the three-dimensional segment calculation airspace to obtain the segment baseline traffic flow; S3, collecting ADS-B data of the study area and calculating the number of aircraft in the three-dimensional segment calculation airspace, and then calculating the segment congestion index; S4, using the proportion of the study aircraft's segment distance to the total airway distance as the weight of each segment, filtering the segment congestion index of all segments of the study aircraft's entire airway, and calculating the overall airway congestion index of the study aircraft. This invention can accurately quantify and obtain the segment congestion index and the airway congestion index, and can measure and characterize the airway congestion situation.

Description

Technical Field

[0001] This invention relates to the field of airway congestion measurement, and in particular to a method and system for measuring airway congestion index based on airway segment spatial analysis. Background Technology

[0002] Just as vehicles experience traffic congestion on roads, civil aviation aircraft also face similar challenges in airspace. With the rapid development of the civil aviation industry, air traffic volume is constantly increasing, making airway congestion a more prominent issue. This is especially true during peak hours. Current research on airway congestion indices in the civil aviation field employs two main approaches: one is based on flight plans for traffic forecasting. However, an aircraft's flight plan represents a single airway, and air traffic involves multiple aircraft flights. This method of congestion prediction based on flight plans is heavily influenced by human factors and lacks flexibility. The other approach, with the widespread application of ADS-B (Automatic Dependent Surveillance-Broadcast) technology, estimates airway traffic through ADS-B data simulation and statistical analysis. This typically involves manual calculations, which lack consideration for airway width and altitude, have significant statistical errors, and cannot quantify the congestion status of airways and their different locations. Summary of the Invention

[0003] The purpose of this invention is to provide a method and system for calculating the airway congestion index based on airway segment spatial analysis. This method divides the airway into segments using each intersection point as a dividing point, and calculates the corresponding three-dimensional airspace based on these segments in the airway GIS spatial model. It selects the segment baseline flow rate for each three-dimensional airspace calculation using segment flow data datasets from N1 years and N2 days of historical flight data. It also collects real-time flight data from ADS-B data to classify segment congestion indices, and obtains the airway congestion index through segment weighting.

[0004] The objective of this invention is achieved through the following technical solution:

[0005] A method for calculating the airway congestion index based on airway segment spatial analysis, the method comprising:

[0006] S1. Construct a route GIS spatial model of the study area. Based on the intersection of the flight plan route projections of all aircraft in the study area as the dividing points, divide all flight plan routes into several segments. In the route GIS spatial model, calculate the airspace i corresponding to the three-dimensional segments of the segments based on the segment division.

[0007] S2. Collect historical ADS-B data of all aircraft flights in the study area for N1 years. The ADS-B data includes time, latitude and longitude, and elevation information. The historical ADS-B data is counted hourly according to time series and aggregated into the three-dimensional flight segment calculation airspace to obtain the flight segment flow count set A for N1 years, aggregated by hour. Take the ADS-B data of historical flights over the past N2 days of the study date and count them hourly, aggregating them into the three-dimensional flight segment calculation airspace to obtain the flight segment flow count set B for N2 days, aggregated by hour. Merge the flight segment flow count set A and the flight segment flow count set B, and extract the median or average of each hour as the flight segment baseline flow R for the three-dimensional flight segment calculation airspace i and hour m. im ;

[0008] S3. Collect ADS-B data of the study area and calculate the number of aircraft Z in the three-dimensional flight segment calculation airspace i one hour before the study time. im The segment congestion index M for the three-dimensional flight segment calculation airspace i and hour m is calculated according to the following formula. im ,

[0009] S4. The weight of each flight segment is determined by the proportion of the aircraft's flight segment distance to the total route distance. im The segment congestion index of all segments of the entire flight route of the aircraft under study is screened, and the predicted route congestion index C of the entire flight route of the aircraft under study at the study time is calculated according to the following formula. im , P represents the sequence number of the route segment.

[0010] This invention provides a second method for calculating the airway congestion index. The difference between this method and the first method lies in the following: In method S1, the three-dimensional airway calculation airspace i is subdivided into several altitude layers; in method S2, the airway calculation airspace i is subdivided into altitude layers to obtain the airway baseline flow R for the three-dimensional airway calculation airspace i, altitude layer j, and hour m. i，j，m In method S3, the reference flow rate M for the three-dimensional flight segment calculation airspace i, altitude layer j, and hour m is obtained at different altitude layers. i，j，m , Z i，j，m To collect ADS-B data for the study area and calculate the number of aircraft in the three-dimensional flight segment computational airspace i, altitude layer j, and one hour before the study time m; in method S4, based on the flight plan route of the study aircraft, the segment congestion index of all segments of the entire route at the corresponding altitude layer is screened, and the predicted route congestion index C of the entire route of the study aircraft at the study time is calculated according to the following formula. im ,

[0011] Preferably, the altitude layer of the three-dimensional flight segment calculation airspace i is divided according to the true flight path angle.

[0012] Preferably, the width of the three-dimensional flight segment calculation airspace i is N3 kilometers extending to the left and right of the flight segment as the center.

[0013] Preferably, the method for calculating the airway congestion index of the present invention further includes method S5:

[0014] S5. Obtain the airway congestion index for all aircraft in the study area according to methods S1 to S4, and visualize it in the airway GIS spatial model.

[0015] Preferably, the projection plane of the route GIS spatial model is a latitude and longitude plane, and the route congestion index of all aircraft is projected onto the projection plane of the route GIS spatial model for visualization.

[0016] Preferably, all aircraft routes in the study area are visualized according to segment congestion index, with the segment congestion index represented by colors from green to red from low to high.

[0017] Preferably, in method S2, N1 is taken as three years, and the collected historical ADS-B data is subjected to data denoising, cleaning and invalid data removal processing, including the removal of outliers and duplicates.

[0018] A system for calculating airway congestion index based on segmented spatial analysis of airways includes an ADS-B data acquisition module, an airway GIS spatial model, a segment airspace division module, a segment baseline flow calculation module, and an airway congestion index calculation module. The ADS-B data acquisition module collects historical ADS-B data from all aircraft in the study area over N1 years. The ADS-B data includes time, latitude and longitude, and elevation information. The segment airspace division module divides all planned flight routes into segments based on the intersection points of the flight plan route projections of all aircraft in the study area. In the inter-segment model, the three-dimensional segment calculation airspace i is divided based on the segment division. The segment baseline flow calculation module collects historical ADS-B data by hourly counts according to the time series into the three-dimensional segment calculation airspace, obtaining the segment flow count set A collected by hour for N1 years. The ADS-B data of historical flights in the past N2 days are taken and collected by hourly counts into the three-dimensional segment calculation airspace, obtaining the segment flow count set B collected by hour for N2 days. The segment flow count set A and the segment flow count set B are merged, and the median or average of each hour is taken as the segment baseline flow R of the three-dimensional segment calculation airspace i and hour m. imThe ADS-B data acquisition module also acquires ADS-B data of the study area and inputs it into the route congestion index calculation module. The route congestion index calculation module is used to calculate the number of aircraft Z in the three-dimensional flight segment calculation airspace i one hour m before the study time. im The segment congestion index M for the three-dimensional flight segment calculation airspace i and hour m is calculated according to the following formula. im , The route congestion index calculation module uses the proportion of the aircraft's segment distance to the total route distance as the weight W for each segment. im The segment congestion index of all segments of the entire flight route of the aircraft under study is screened, and the predicted route congestion index C of the entire flight route of the aircraft under study at the study time is calculated according to the following formula. im ,

[0019] Preferably, the airway congestion index calculation system of the present invention further includes a visualization module, which is used to visualize the airway GIS spatial model.

[0020] Compared with the prior art, the present invention has the following advantages and beneficial effects:

[0021] (1) This invention divides each intersection of the route into segments and calculates the three-dimensional segment airspace corresponding to the segment in the route GIS spatial model. The segment reference flow is selected from the segment flow dataset of N1-year historical flight data and N2-day historical flight data. Real-time flight data of ADS-B data is collected and classified to calculate the segment congestion index of the three-dimensional segment calculation airspace. The route congestion index is obtained by segment weighting.

[0022] (2) The present invention can obtain traffic congestion index data of air route network and visualize it, and can accurately and timely reflect the changes in air traffic congestion in air route network. Attached Figure Description

[0023] Figure 1 This is a schematic diagram of the method flow for calculating the airway congestion index of the present invention;

[0024] Figure 2 This is a schematic diagram illustrating the height layer division as shown in the embodiment;

[0025] Figure 3 This is a schematic diagram illustrating the vectorization and transformation of flight plan and route data for all aircraft in a certain region, as an example of a study area, and its construction within a route GIS spatial model.

[0026] Figure 4 This is a schematic diagram illustrating the segmentation of a planned flight route as an example in the embodiment;

[0027] Figure 5 This is a schematic diagram illustrating the projection of the route congestion index onto the projection plane of the route GIS spatial model, as exemplified in this embodiment. Detailed Implementation

[0028] The present invention will be further described in detail below with reference to embodiments:

[0029] Example 1

[0030] like Figure 1 As shown, a method for calculating the airway congestion index based on airway segment spatial analysis includes the following steps:

[0031] S1. Construct a route GIS spatial model of the study area. Based on the intersection points of the flight plan route projections of all aircraft in the study area as dividing points, divide all flight plan routes into several segments. (Preferably, the flight plan routes of all aircraft in the study area are constructed in the route GIS spatial model using vectorized data. The vectorized flight plan route data is as follows...) Figure 3 (As shown). Taking the flight plan route B208 as an example, as... Figure 4 The diagram illustrates the division of flight segments using route intersections as dividing points. In the route GIS spatial model, the airspace i corresponding to the three-dimensional flight segment is calculated based on the segment division. The study area in this embodiment is selected according to the actual research situation; for example, it could be a specific civil aviation control area or several civil aviation control areas, or even a specific country region.

[0032] In some embodiments, the altitude layer of the three-dimensional flight segment calculation airspace i is divided according to the true flight path angle. Specific examples are as follows: Figure 2 As shown, the true flight path angles (in this embodiment, the flight path angles are divided according to the flight plan routes of all aircraft in the study area, measured from the start and turning points) are in the range of 0 to 179 degrees. At altitudes from 900 meters to 8100 meters, each altitude layer is 600 meters apart; at altitudes from 8900 meters to 12500 meters, each altitude layer is 600 meters apart; and at altitudes above 12500 meters, each altitude layer is 1200 meters apart. The true flight path angles are in the range of 180 to 359 degrees. At altitudes from 600 meters to 8400 meters, each altitude layer is 600 meters apart; at altitudes from 9200 meters to 12200 meters, each altitude layer is 600 meters apart; and at altitudes above 13100 meters, each altitude layer is 1200 meters apart.

[0033] In some embodiments, the width of the three-dimensional segment calculation airspace i is N3 kilometers extending left and right from the center of the segment (for example, N3 kilometers is 10 kilometers).

[0034] S2. Collect historical ADS-B data of all aircraft flights in the study area for N1 years (in some embodiments, N1 is taken as three years; the collected historical ADS-B data undergoes data denoising, cleaning, and invalid data removal, including removing outliers and duplicates, to ensure data accuracy and reliability). The ADS-B data includes time, latitude and longitude, and elevation information. The historical ADS-B data is divided into hourly segments according to the time series (24 hours a day, collected hourly, with each hourly segment grouped together for calculation; due to large differences in traffic flow at different times, this embodiment uses detailed segmentation, dividing the day into 24 hours, and calculating the traffic flow value separately for each hour, so that each flight segment has a baseline traffic flow for each hour). The counts are collected in the three-dimensional flight segment calculation airspace to obtain the flight segment traffic count set A for N1 years, collected by hour. For example: If we define s as a flight segment dataset, y as the year, m as the month, d as the day, h as the hour, H as the altitude layer, and D as the flight segment direction, we establish a flight segment traffic dataset: Q(s, y, m, d, h, H, D). Q(s, y, m, d, h, H, D) represents the traffic flow value of flight segment s at altitude layer H and direction D, on day d and hour h in year y. We construct the flight segment traffic function: R(s, h, H, D). R(s, h, H, D) represents the flight segment traffic flow of flight segment s at hour h. The data correspondence between the flight segment traffic function and the flight segment traffic dataset is: R(s, h, H, D) = {Q(s, h, d, m, y, H, D)}.

[0035] ADS-B data from historical flights over the past N2 days are collected and counted hourly in the three-dimensional segment calculation airspace to obtain segment flow count set B for N2 days, collected hourly. Segment flow count set A and segment flow count set B are merged, and the median or average of each hour is extracted (the median is used in this example) as the segment baseline flow R for three-dimensional segment calculation airspace i and hour m. im For example: N1 is a three-year period, N2 is a ten-day period, and the traffic data for the same month (m), the same day (d), and the same hour from the historical data of y1, y2, and y3 over the past three years are as follows:

[0036] q1=Q(s,h,d,m,y1,H,D)

[0037] q1=Q(s,h,d,m,y1,H,D)

[0038] q3=Q(s,h,d,m,y3,H,D)

[0039] The following is the traffic data for the same hour over the past 10 days:

[0040] p n ={Q(s, h, d)} n ,m,y,H,D)},

[0041] Then R(s, h, H, D) = {Q (s, h, d, m, y, H, D)} = [q1, q2, q3, p1, p2, p3...p 10 In this embodiment, the median is taken as the reference flow rate for flight segment s at altitude H, direction D, when the flight segment is at altitude h.

[0042] Baseline(s,h,H,D)=Median(R(s,h,H,D)).

[0043] S3. Collect ADS-B data of the study area and calculate the number of aircraft Z in the three-dimensional flight segment calculation airspace i one hour before the study time. im The segment congestion index M for the three-dimensional flight segment calculation airspace i and hour m is calculated according to the following formula. im , For example: The real-time traffic flow of a flight segment is at altitude H in the direction of flight segment D, and the number of aircraft passing through in the hour preceding the current time h is z, so A(s, h, H, D) = z. The segment congestion index M is the ratio of the real-time traffic flow to the baseline traffic flow of the segment. The congestion index of flight segment s at altitude H in the direction of flight segment D is:

[0044] The congestion index for each flight segment is calculated by weighting the traffic flow percentage for each altitude level and direction as follows:

[0045]

[0046] The congestion index for flight segment s is as follows:

[0047]

[0048] S4. The weight of each flight segment is determined by the proportion of the aircraft's flight segment distance to the total route distance. im The segment congestion index of all segments of the entire flight route of the aircraft under study is screened, and the predicted route congestion index C of the entire flight route of the aircraft under study at the study time is calculated according to the following formula. im , P represents the sequence number of the route segment.

[0049] S5. Obtain the route congestion index for all aircraft in the study area according to methods S1-S4, and visualize it in the route GIS spatial model. In some embodiments, the projection plane of the route GIS spatial model (the projection plane is the latitude and longitude plane of the route GIS spatial model excluding elevation data) is the latitude and longitude plane, and the route congestion index of all aircraft is projected onto the projection plane of the route GIS spatial model for visualization, such as... Figure 5The diagram illustrates the projection of the airway congestion index of a study area onto the projection plane of the airway GIS spatial model. In some embodiments, all aircraft airways in the study area are visualized according to segment congestion indices, with the segment congestion indices represented by colors from green to red, from low to high.

[0050] Example 2

[0051] A method for calculating the airway congestion index based on airway segment spatial analysis, the method comprising:

[0052] S1. Construct a route GIS spatial model of the study area. Based on the intersection of the flight plan route projections of all aircraft in the study area as the dividing points, divide all flight plan routes into several segments. In the route GIS spatial model, divide the corresponding three-dimensional segment calculation airspace i based on the segment. The three-dimensional segment calculation airspace i is subdivided into several altitude layers. For example, the three-dimensional segment calculation airspace i includes J altitude layers, which are ordered from bottom to top. Altitude layer j is the j-th altitude layer.

[0053] S2. Collect historical ADS-B data of all aircraft flights in the study area for N1 years. The ADS-B data includes time, latitude and longitude, and elevation information. The historical ADS-B data is counted hourly according to the time series and aggregated into the corresponding altitude layer of the three-dimensional flight segment calculation airspace to obtain the flight segment flow count set A for N1 years, aggregated by hour. Take the ADS-B data of historical flights over the past N2 days of the study date and count them hourly, aggregating them into the corresponding altitude layer of the three-dimensional flight segment calculation airspace to obtain the flight segment flow count set B for N2 days, aggregated by hour. Merge the flight segment flow count set A and the flight segment flow count set B, and extract the median or average of each hour (the median is used in this example) as the flight segment baseline flow R for the three-dimensional flight segment calculation airspace i and hour m. im The baseline flow rate R for the three-dimensional flight segment calculation airspace i, altitude layer j, and hour m is obtained at each of the three-dimensional flight segment calculation airspace i, altitude layer j, and hour m. i，j，m .

[0054] S3. Collect ADS-B data of the study area and calculate the number of aircraft Z in the three-dimensional flight segment calculation airspace i one hour before the study time. i，j，m Z i，j，m To collect ADS-B data for the study area and calculate the number of aircraft in the three-dimensional flight segment computational airspace i, altitude layer j, and hour m prior to the study time, the flight segment baseline flow M in the three-dimensional flight segment computational airspace i, altitude layer j, and hour m is calculated using the following formula. i，j，m ,

[0055] S4. The weight of each flight segment is determined by the proportion of the aircraft's flight segment distance to the total route distance.im The segment congestion index of all segments of the entire flight route of the aircraft under study is screened, and the predicted route congestion index C of the entire flight route of the aircraft under study at the study time is calculated according to the following formula. im , P represents the sequence number of the route segment.

[0056] S5. Obtain the airway congestion index for all aircraft in the study area according to methods S1 to S5, and visualize it in the airway GIS spatial model.

[0057] Example 3

[0058] A system for calculating airway congestion index based on airway segment spatial analysis includes an ADS-B data acquisition module, an airway GIS spatial model, a segment airspace division module, a segment baseline flow calculation module, and an airway congestion index calculation module. The ADS-B data acquisition module collects historical ADS-B data for all aircraft's flights in the study area over a period of N1 years. This ADS-B data includes time, latitude, longitude, and elevation information. The segment airspace division module divides all planned flight routes into segments based on the intersection points of the flight plan route projections of all aircraft in the study area. In the airway GIS spatial model, the corresponding three-dimensional segment calculation airspace i is calculated based on the segment division. The segment baseline flow calculation module aggregates the historical ADS-B data by hourly counts according to the time series into the three-dimensional segment calculation airspace, obtaining the segment flow count set A for year N1, aggregated by hour. ADS-B data from the past N2 days of historical flights are collected and counted hourly in the three-dimensional flight segment calculation airspace, resulting in a flight segment flow count set B for N2 days, collected hourly. The flight segment flow count set A and the flight segment flow count set B are merged, and the median or average of each hour is extracted as the flight segment baseline flow R for the three-dimensional flight segment calculation airspace i and hour m. im The ADS-B data acquisition module also collects ADS-B data of the study area and inputs it into the route congestion index calculation module. The route congestion index calculation module is used to calculate the number of aircraft Z in the three-dimensional flight segment calculation airspace i one hour before the study time. im The segment congestion index M for the three-dimensional flight segment calculation airspace i and hour m is calculated according to the following formula. im , The route congestion index calculation module uses the proportion of the aircraft's segment distance to the total route distance as the weight W for each segment. im The segment congestion index of all segments of the entire flight route of the aircraft under study is screened, and the predicted route congestion index C of the entire flight route of the aircraft under study at the study time is calculated according to the following formula. im ,

[0059] The airway congestion index calculation system of the present invention also includes a visualization module, which is used to visualize the airway GIS spatial model.

[0060] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A method for calculating the airway congestion index based on airway segment spatial analysis, characterized in that: The methods include: S1. Construct a route GIS spatial model of the study area. Based on the intersection of the flight plan route projections of all aircraft in the study area as the dividing points, divide all flight plan routes into several segments. In the route GIS spatial model, calculate the airspace i corresponding to the three-dimensional segments of the segments based on the segment division. S2. Collect historical ADS-B data of all aircraft flights in the study area for N1 years. The ADS-B data includes time, latitude and longitude, and elevation information. The historical ADS-B data is counted hourly according to time series and aggregated into the three-dimensional flight segment calculation airspace to obtain flight segment flow count set A for N1 years, aggregated by hour. Take the ADS-B data of historical flights for the past N2 days of the study date and count them hourly, aggregating them into the three-dimensional flight segment calculation airspace to obtain flight segment flow count set B for N2 days, aggregated by hour. Merge flight segment flow count set A and flight segment flow count set B, and extract the median or average for each hour as the flight segment baseline flow for the three-dimensional flight segment calculation airspace i and hour m. ; S3. Collect ADS-B data of the study area and calculate the number of aircraft in the three-dimensional flight segment calculation airspace i one hour before the study time m. The segment congestion index for the three-dimensional flight segment calculation airspace i and hour m is calculated according to the following formula. , ; S4. The weight of each flight segment is determined by the proportion of the aircraft's flight segment distance to the total route distance. The segment congestion index of all segments of the entire flight route of the aircraft under study is screened, and the predicted route congestion index of the entire flight route of the aircraft under study at the study time is calculated according to the following formula. , , This refers to the sequence number of a segment within the air route.

2. The method for calculating the airway congestion index based on airway segment spatial analysis according to claim 1, characterized in that: In method S1, the three-dimensional flight segment calculation airspace i is subdivided into several altitude layers; in method S2, the flight segment reference flow is obtained by subdividing the three-dimensional flight segment calculation airspace i into altitude layers, and obtaining the flight segment reference flow for three-dimensional flight segment calculation airspace i, altitude layer j, and hour m respectively. In method S3, the baseline flow rates for the three-dimensional flight segment calculation airspace i, altitude layer j, and hour m are obtained at different altitude layers. , ,in To collect ADS-B data for the study area and calculate the number of aircraft in the three-dimensional flight segment computational airspace i, altitude layer j, and one hour before the study time m; in method S4, based on the flight plan route of the study aircraft, the segment congestion index of all segments of the entire route at the corresponding altitude layer is screened, and the predicted route congestion index of the entire route of the study aircraft at the study time is calculated according to the following formula. , .

3. The method for calculating the airway congestion index based on airway segment spatial analysis according to claim 2, characterized in that: The altitude layer of the three-dimensional flight segment calculation airspace i is divided according to the true flight path angle.

4. The method for calculating the airway congestion index based on airway segment spatial analysis according to claim 1 or 2, characterized in that: The width of the calculated airspace i for the three-dimensional flight segment is N3 kilometers extending to the left and right of the flight segment as the center.

5. The method for calculating the airway congestion index based on airway segment spatial analysis according to claim 1, characterized in that: It also includes method S5: S5. Obtain the airway congestion index for all aircraft in the study area according to methods S1 to S4, and visualize it in the airway GIS spatial model.

6. The method for calculating the airway congestion index based on airway segment spatial analysis according to claim 5, characterized in that: The projection plane of the route GIS spatial model is a latitude and longitude plane, and the route congestion index of all aircraft is projected onto the projection plane of the route GIS spatial model for visualization.

7. The method for calculating the airway congestion index based on airway segment spatial analysis according to claim 5, characterized in that: The study area visualizes the congestion index of all aircraft routes according to flight segments, with the congestion index ranging from low to high using a color scheme from green to red.

8. The method for calculating the airway congestion index based on airway segment spatial analysis according to claim 1, characterized in that: In method S2, N1 is taken as three years. The collected historical ADS-B data is subjected to data denoising, cleaning and invalid data removal. Invalid data removal includes removing outliers and duplicates.

9. A system for calculating airway congestion index based on airway segment spatial analysis, characterized in that: The system includes an ADS-B data acquisition module, a route GIS spatial model, a segment airspace division module, a segment baseline flow calculation module, and a route congestion index calculation module. The ADS-B data acquisition module collects historical ADS-B data for all aircraft's N1-year flight history in the study area. This ADS-B data includes time, latitude, longitude, and elevation information. The segment airspace division module divides all planned flight routes into segments based on the intersection points of the flight plan route projections of all aircraft in the study area. These segments are then defined in the route GIS spatial model. The corresponding three-dimensional flight segment calculation airspace i; the flight segment baseline flow calculation module collects historical ADS-B data by hourly count according to time series in the three-dimensional flight segment calculation airspace, obtaining the flight segment flow count set A collected by hour for N1 years; takes the historical ADS-B data of the past N2 days of the study date and collects it by hourly count in the three-dimensional flight segment calculation airspace, obtaining the flight segment flow count set B collected by hour for N2 days; merges the flight segment flow count set A and the flight segment flow count set B, and takes the median or average of each hour as the flight segment baseline flow for the three-dimensional flight segment calculation airspace i and hour m. ; The ADS-B data acquisition module also collects ADS-B data of the study area and inputs it into the route congestion index calculation module. The route congestion index calculation module is used to calculate the number of aircraft in the three-dimensional flight segment calculation airspace i one hour m before the study time. The segment congestion index for the three-dimensional flight segment calculation airspace i and hour m is calculated according to the following formula. , ; The route congestion index calculation module uses the proportion of the aircraft's segment distance to the total route distance as the weight of each segment. The segment congestion index of all segments of the entire flight route of the aircraft under study is screened, and the predicted route congestion index of the entire flight route of the aircraft under study at the study time is calculated according to the following formula. , .

10. The route congestion index calculation system based on route segment spatial analysis according to claim 9, characterized in that: It also includes a visualization module, which is used to visualize the airway GIS spatial model.

Images