Method and system for calculating cross-sectional distance of bus routes

Abstract

The present invention discloses a method and system for calculating cross-sectional distances for bus routes. The method comprises: Step 1: collecting route linear data and station data; Step 2: standardizing the linear data for each route; Step 3: preliminarily calculating the station bearing angles for each station on each route; and Step 4: calculating the cross-sectional distance for each route. The present invention significantly improves the adaptability of calculating cross-sectional distances for bus routes using route linear data and route station data. The present invention provides excellent calculation results for complex situations such as loop-like routes, routes that pass through stations twice, and individual station location data with significant deviations. The present invention can accurately obtain cross-sectional distances, further optimizing network layout and improving operational efficiency.

Description

Technical Field

[0001] The present invention relates to the field of urban intelligent transportation technology, and in particular to a method and system for calculating the cross-sectional distance of a bus route. Background Art

[0002] The cross-sectional distance of a bus route (the distance along the road centerline between two adjacent stops on a route is the cross-sectional distance between the two stops) is a key indicator of urban public transportation. In bus route planning, calculating cross-sectional distance is crucial for optimizing network layout, improving operational efficiency, and improving service quality. This is primarily due to the following factors:

[0003] 1. Assessing capacity needs:

[0004] Calculate the passenger flow per unit length (people / km) through cross-sectional distance, identify high-demand sections, and analyze whether the transport capacity of sections with high passenger flow density is sufficient.

[0005] 2. Optimize operational efficiency:

[0006] Analyze the section speed of each section through section distance and travel time, find out the low-speed section and further analyze the improvement measures.

[0007] 3. Evaluate connection efficiency

[0008] The cross-sectional distance needs to be coordinated with the location of connecting facilities (such as subway stations and bicycle parking spots) to improve the efficiency of multimodal transport.

[0009] 4. Meet planning indicators and specifications:

[0010] Industry standard requirements: Many cities’ bus planning specifications specify the cross-sectional distance of routes (e.g. ≥250m ≤500m) to ensure a reasonable network layout.

[0011] 5. Analyze public transportation indicators:

[0012] Passenger travel distance: Calculate the route distance between the boarding and alighting stops based on the two stops;

[0013] Line operating load rate: actual passenger turnover / rated vehicle turnover.

[0014] In short, cross-sectional distance is fundamental geographic data in bus network planning, directly impacting the practicality, economy, and sustainability of routes. Its calculation requires integration with GIS (Geographic Information System) tools or field measurements to achieve scientific decision-making.

[0015] There are currently two methods for calculating bus section distance:

[0016] The main disadvantage of this method is the high cost, as it requires actual driving.

[0017] The calculation is performed using GIS-related methods based on the line geometry data (data representing the direction of a line using an ordered sequence of longitude and latitude is called line geometry data) and the location data of line stations.

[0018] Linear data is an ordered sequence of longitude and latitude points. Each longitude and latitude point is called a linear point. Linear points are numbered from 1 to N ; There is a straight line segment between two adjacent linear points, which is represented by the numbers of the previous and next linear points. For example Figure 1 The linear shape has 7 linear points, numbered from 1 to 7; there are 6 linear segments: segment 1-2, up to segment 6-7. The idea of ​​this method is to project the station onto the linear data, and then calculate the spherical distance between the projection points. Draw a perpendicular line from the station to the segment, the foot of the perpendicular is the projection point of the station on the linear shape, and the distance from the point to the segment is the projection distance. If the projection distance is less than a certain threshold (such as 30 meters), it is considered possible to project. In addition, if the station is on a cusp and cannot be projected onto the segment, the linear point closest to it (if the distance is less than the threshold) is taken as the projection point, such as Figure 2 shown.

[0019] The main problem of this method is that it has high requirements for the accuracy of site location and linear data, and it cannot handle the situation when a site has multiple close projection points, which makes it impossible to accurately obtain the cross-sectional distance, thus affecting the network layout and operational efficiency.

[0020] The following example illustrates the problem:

[0021] 1. If Figure 3 In the two cases (a) and (b) shown, it is impossible to determine the exact projection point based on the projection distance. Many cities have such bus route data.

[0022] (a) The route follows a loop, heading east and then turning west, with stations on the north and south sides of the road. The station location or alignment data is inaccurate here, as the station is drawn between two roads, making it difficult to determine whether it is on the south or north alignment.

[0023] (b) The route and station locations are the same as in (a), except that this road is a narrow, two-way road with no median, so the drawn lines overlap in this section. In this extreme case, even if the data is completely accurate, it is impossible to accurately calculate the station's projection point on the line based solely on the location data, and therefore it is impossible to accurately calculate the cross-sectional distance.

[0024] 2. If Figure 4 As shown in the figure, the line shape passes through station A twice. Based on the position data alone, it is impossible to know whether the station should be projected onto the line shape that passes through the first time or the line shape that passes through the second time. This situation is rare but it does exist.

[0025] 3. If Figure 5 As shown, the latitude and longitude data of some stations are incorrect. It seems that the stations are on the northern line and far away from the southern line, but in fact the stations are on the southern line. Summary of the Invention

[0026] The technical problem to be solved by the embodiments of the present invention is to provide a method and system for calculating the cross-sectional distance of a bus line, so as to accurately obtain the cross-sectional distance, optimize the line network layout and improve the operational efficiency.

[0027] In order to solve the above technical problems, an embodiment of the present invention proposes a method for calculating the cross-sectional distance of a bus route, comprising:

[0028] Step 1: Collect the line linear data and station data;

[0029] Step 2: Standardize the linear data of each line;

[0030] Step 3: Preliminary calculation of the station direction angle of each station on each line;

[0031] Step 4: Preliminarily calculate the projection points of each station on the line. If there are multiple projection points on the line for a station, simplify the projection points based on the station bearing angle of the station. List the combinations of station projection information sequences on the line, find the most likely combination of line station projection information sequences, and recalculate unreliable projection information. Update the station bearing angle of each station on each line, and then calculate the cross-sectional distance of the line segment. Among them, the most likely meaning is: the station projection sequence contains the largest number of ascending station projection numbers; the smaller of the two linear point numbers of the line segment where the projection point is located is the projection number.

[0032] Accordingly, an embodiment of the present invention further provides a bus route cross-sectional distance calculation system, comprising:

[0033] Data preparation module: collects line linear data and station data;

[0034] Data standardization module: standardizes the linear data of each line;

[0035] Direction angle calculation module: preliminarily calculates the station direction angle of each station on each line;

[0036] Section distance calculation module: Preliminary calculation of the projection points of each station on the line. If there are multiple projection points on the line for a station, the projection points are simplified according to the station direction angle of the station. The combinations of station projection information sequences on the line are listed to find the most likely combination of line station projection information sequences, and the unreliable projection information is recalculated. The station direction angle of each station on each line is updated, and then the section distance of the line segment is calculated. Among them, the most likely meaning is: the station projection sequence contains the largest number of ascending station projection numbers; the smaller of the two linear point numbers of the line segment where the projection point is located is the projection number.

[0037] The beneficial effects of the present invention are as follows: the present invention can greatly improve the adaptability of calculating the cross-sectional distance of bus routes through line route data and line station data. The present invention can have good calculation results for complex situations such as quasi-loop lines, lines passing through stations twice, and individual station position data with large deviations. The present invention can accurately obtain cross-sectional distances, further optimize line network layout and improve operational efficiency. BRIEF DESCRIPTION OF THE DRAWINGS

[0038] Figure 1 It is a schematic diagram of the latitude and longitude points of linear data.

[0039] Figure 2 It is a schematic diagram of the site projected onto the line.

[0040] Figure 3 It is a schematic diagram showing that the projection point of a station on a line cannot be determined based on the projection distance. (a) is a schematic diagram showing that the station is between two line segments, and (b) is a schematic diagram showing that the two line segments overlap.

[0041] Figure 4 It is a schematic diagram showing that the line passes through the site twice.

[0042] Figure 5 This is a schematic diagram when the data deviation of individual stations is large.

[0043] Figure 6 It is a schematic diagram of the projection numbers and direction angles of the sites of the present invention.

[0044] Figure 7 It is a schematic diagram of the standardized line shape data of the present invention.

[0045] Figure 8 It is a schematic diagram of the orientation angle of the preliminary calculation site of the present invention.

[0046] Figure 9 It is a flowchart of a method for calculating cross-sectional distances of bus routes according to an embodiment of the present invention. DETAILED DESCRIPTION

[0047] It should be noted that, unless there is a conflict, the embodiments in this application and the features in the embodiments can be combined with each other. The present invention is further described in detail below with reference to the drawings and specific embodiments.

[0048] In the embodiments of the present invention, if there are directional indications (such as up, down, left, right, front, back, etc.), they are only used to explain the relative position relationship and movement status of the various components under a specific posture (as shown in the accompanying drawings). If the specific posture changes, the directional indication will also change accordingly.

[0049] In addition, the terms "first," "second," and so on, used in this disclosure are for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of technical features being referred to. Therefore, features specified as "first" or "second" may explicitly or implicitly include at least one of these features.

[0050] Please refer to Figure 9 The method for calculating the cross-sectional distance of a bus route according to an embodiment of the present invention includes steps 1 to 4.

[0051] First, two concepts are introduced: the projection number of the site and the site direction angle.

[0052] The smaller linear point number of the line segment where the projection point is located is the projection number. Figure 6 As shown, the projection point of site A is on line segment 2-3, and its projection number is 2; the projection point of site B is on line segment 4-5, and its projection number is 4.

[0053] The angle corresponding to the tangent direction of the projection point of the station on the linear shape is the direction angle of the station. Let the north be 0°, and the clockwise direction is from 0° to 359°. Its actual meaning is the forward driving direction angle of the road where the station is located. Figure 6 As shown, the tangent direction of the projection point of station A is the direction of vector 2->3, which is 105°; the tangent direction of the projection point of station B is the direction of vector 4->5, which is 95°.

[0054] Secondly, let’s introduce the concept: the two angles are in the same direction.

[0055] Because multiple lines may pass through the same station, the calculated bearing angles may differ for different lines. If the data is generally accurate, the angles may not be exactly the same, but the directions will be consistent. For example, if one angle is 85° and the other is 95°, both are strictly facing east.

[0056] The direction angles θ1 and θ2 of two stations are defined as being consistent if their difference is no greater than a given threshold (e.g., 60°). Generally, the two possible direction angles calculated for the same station differ by approximately 0° or 180°, corresponding to the two sides of the road. Note that the direction angle range is [0°, 359°]. The absolute value of the difference cannot be directly considered. The angle difference Δ must be calculated as follows:

[0057] △=|θ1-θ2|, if△>180°then△=360°-△;

[0058] In addition, if higher accuracy is sought during specific calculations, the longitude and latitude of the geographic coordinate system can be converted into meters of the plane rectangular coordinate system. However, generally, directly approximating the longitude and latitude as the values ​​of the horizontal and vertical axes of the rectangular coordinate system can also meet the accuracy requirements.

[0059] Step 1: Data preparation. Line section distance calculation requires complete line basic data:

[0060] 1) Linear data: Linear data is an ordered sequence of longitude and latitude points.

[0061] 2) Route station data: which stations are included in the route, the station numbers and their latitude and longitude locations.

[0062] Step 2: Standardize the linear data of each line (assuming there are N lines in total):

[0063] 1) Linear data standardization, such as Figure 7 As shown in the following example: a) Duplicate linear points or those that are too close together (e.g., 3 meters) are deleted; b) if the distance between two points is far (e.g., greater than 50 meters), linear points are added evenly spaced in between. This is done so that the projection numbers of the stations on the linear will be assigned later, and it is important that the projection numbers of the two stations are unique.

[0064] 2) Ensure that the site locations and linear latitude and longitude points are in the same coordinate system, such as the GCJ02 coordinate system, CGCS2000 coordinate system, etc. Perform coordinate system conversion if there is any inconsistency.

[0065] Step 3: Preliminary calculation of the station's direction angle:

[0066] For each line, project the station onto the line and calculate the station's bearing angle. If a station has multiple projection points on the line, it means that the bearing angle cannot be determined at the moment, and the station is omitted, that is, only the bearing angle of the station that can be determined at this time is retained.

[0067] like Figure 8As shown, site A can be projected onto line segment 2-3 with a bearing angle of 110°. Site A can also be projected onto line segment 12-13 with a bearing angle of 283°. Similarly, site B projects onto line segment 6-7 with a bearing angle of 135°. Site C projects onto line segment 9-10 with a bearing angle of 279°. Site D projects onto line segment 3-4 with a bearing angle of 95°. Site E projects onto line segment 1-2 with a bearing angle of 110° or line segment 13-14 with a bearing angle of 283°.

[0068] The actual bearing angle for a station is only one. We are currently unsure which station to select, A or E. This step will only retain stations with only one bearing angle result. For this example route, retain the bearing angle information for stations B, C, and D.

[0069] In addition, for site E, if there is only one projection for another line and the calculated direction angle is 280°, the direction angle information of site E, 280°, is also retained.

[0070] If different lines give inconsistent results for the same station, the majority rule applies and a single azimuth angle is chosen. For example, if three lines pass through a station and the calculated azimuth angles are 358°, 179°, and 1°, respectively, 358° and 1° are consistent, but not 179°. Therefore, either 358° or 1° is chosen.

[0071] Note: This is a preliminary calculation and may be updated later. For example, for Site D, due to inaccurate site or linear data, the current calculated bearing angle is 95°, indicating east. However, the site is actually west-facing. A more detailed calculation will illustrate this with an example, and the site bearing angle will be updated to 275°.

[0072] Step 4, calculate the cross-sectional distance of each line:

[0073] Sub-step 41: Preliminary calculation of the projection points of each line station on the line shape.

[0074] For each station on the route, the line segments on the line are traversed to see if the station can be projected onto the line segment. If it can be projected and the projection distance is less than a certain threshold (such as 30 meters), the projection information of the station is recorded: the corresponding projection point, projection number, and direction angle.

[0075] like Figure 8As shown, site A can be projected onto line segment 2-3 or line segment 12-13, and its projection information is recorded as [(2,p1,110),(12,p2,283)], where the triple (2,p1,110) represents the projection number, projection point, and direction angle, respectively. Site B is located at a cusp and cannot be projected onto a line segment. In this case, the nearest linear point 6 is used as the projection point, and the projection information is recorded as [(6,p3,135)]. Site C projects onto line segment 9-10: [(9,p4,279)]. Site D projects onto line segment 3-4: [(2,p5,95)]. Site E also has two projection options: [(1,p6,110),(13,p7,283)].

[0076] Sub-step 42: simplifying the projection points.

[0077] According to the station direction angle result in step 3, the situation where the station has multiple projection points is simplified, that is, the projection point information whose station direction angle is consistent with the station direction angle stored in step 3 is selected from the multiple projection point information.

[0078] For example, in the previous step, both sites A and E have two projection points. However, in step 3, the direction angle for site E is 282°. Therefore, the original two projection options for site E, [(1,p6,110),(13,p7,283)], are simplified to [(13,p7,283)], because the directions of 282° and 283° are the same.

[0079] Sub-step 43: List the combinations of projection information sequences of each line station.

[0080] Arrange the projection point information of each station on the line in the order of stations to form a line station projection information sequence. Since each station may have more than one projection point, there will be multiple line station projection information sequences formed in the end. Obviously, the number of sequences = ,in For the i The number of projection points for each site.

[0081] For the example above, there are 5 stations on the line: A, B, C, D, and E. Since station A has 2 projection points and the other stations have 1 projection point each, the final projection point combination of the line stations is indivual:

[0082] Combination 1: Site A (2, p1, 110), Site B (6, p3, 135), Site C (9, p4, 279), Site D (3, p5, 95), Site E (13, p7, 283).

[0083] Combination 2: Site A (12, p2, 283), Site B (6, p3, 135), Site C (9, p4, 270), Site D (3, p5, 95), Site E (13, p7, 283).

[0084] Sub-step 44: Find the most likely combination of route site projection information sequences for each route.

[0085] The most likely meaning is that the site projection sequence contains the maximum number of site projection numbers in ascending order. If a site's projection number is in the maximum ascending projection number sequence, the site's projection information is considered credible. Otherwise, the site's projection information is considered unreliable and requires further processing.

[0086] For example, in the previous example:

[0087] The projection number sequence of combination one is (2, 6, 9, 3, 13), and the maximum number of ascending site projection numbers is "2, 6, 9, 13", which is 4 in total.

[0088] The projection number sequence of combination 2 is (12, 6, 9, 3, 13), and the maximum number of ascending site projection numbers is "6, 9, 13", which is 3 in total.

[0089] Combination 1 contains the largest number of ascending site projection numbers, so combination 1 is selected. The projection information for sites A, B, C, and E in combination 1 is trustworthy. The projection information for site D is untrustworthy because its projection number is not in the largest ascending site projection number sequence.

[0090] If there are multiple projection sequences with the maximum number of ascending stations, this scenario is unlikely to occur in practice, but theoretically possible. This means that existing data alone cannot accurately select the projection point for a specific station. In this case, it is necessary to use bus route station implementation specifications or manually determine the location. Utilizable specifications include: 1) The distance between two adjacent stations is generally greater than 200 meters and less than 600 meters; 2) Stations with the same name on a route (not the initial or final station) generally travel in opposite directions.

[0091] Sub-step 45: recalculate the untrustworthy projection information of each line.

[0092] In the most likely sequence found in the previous step, there may be unreliable site projection information. This step is to recalculate these unreliable site projection information.

[0093] The method involves reducing the projection range of the station onto the linear form and removing the projection distance threshold. Removing the projection threshold ensures that the station will be projected onto the linear form regardless of the projection distance. This assumes that the linear form station data is generally accurate, that is, consistent with general data reality.

[0094] If the k If the projection information of a site is unreliable, search for the projection number of the previous trusted site projection information closest to the site. , and the projection number of the next trusted site projection information closest to the site . Then in the linear [ ] section to re-project. k If there is no reliable site projection information in front of the site, the linear segment starts from the first linear point. k There is no credible site projection information behind the site, and the linear segment ends at the last linear point.

[0095] If the recalculated station direction angle is inconsistent with the station direction angle in the result of step 3, the result of step 3 is updated.

[0096] For example, in the previous example, the projection point of site D in combination 1 needs to be recalculated. D The previous trusted site projection information is the projection information of site C, whose projection number is 9, so the starting point of the re-projected segment is 9+1=10; site D The next trusted site projection information is the projection information of site E, whose projection number is 13, so the end point of the re-projected segment is 12. It can be imagined that site D Finally, the projection is made onto line segment 11-12, and the calculated direction angle is 277°.

[0097] From this example, we can see that although the algorithm did not initially provide correct projection information for site D due to data errors, through subsequent refined steps, the algorithm eventually provided correct projection information.

[0098] Sub-step 46, updating the station direction angle of step 3.

[0099] After the previous step, each station on the line is correctly projected onto the line. However, the calculated station bearing angle may not be consistent with the bearing angle stored in step 3. Since this step is a refined calculation result, if the station bearing angle is inconsistent with the one in step 3, the result in step 3 will be updated.

[0100] For example, for station D, the final result of the previous step is a direction angle of 277°. This is the same as the station angle stored in step 3. D The direction angle of station D is inconsistent with 95°, so the direction angle of station D in step 3 is updated to 277°.

[0101] Sub-step 47, calculating the section distance of each line segment.

[0102] After completing the above steps, all stations are correctly projected onto the line. The spherical distance between the projection points of two adjacent stations is then calculated to determine the distances between each section of the route. This distance is then used for further applications in bus planning. For example, for new route planning, the shortest path through high-demand sections is prioritized when determining route orientation. For existing line optimization, the ratio of cross-sectional passenger flow to distance (passenger flow density) is used to identify inefficient sections and determine whether adjustments are necessary. For competitive analysis, the advantages of bus distance (e.g., "last mile" coverage) compared to other modes of transportation, such as subways and shared bikes, are compared.

[0103] The present invention improves the adaptability of calculating bus line cross-sectional distances using line shape data and station data, that is, solves the problem of how to deal with the situation when there are multiple projection point candidates for a station on the line shape or the data deviation of an individual station is large.

[0104] The main idea of ​​the present invention is to make full use of the sequence number information and the direction angle information of the line stations:

[0105] 1) Traverse each station of each line and calculate the projection point of the station on the line geometry. If there is only one projection point, calculate the direction angle of the station and store the result.

[0106] 2) If there are multiple projection points on the line shape of the previous station, check whether the direction angle information of the station is stored in the result of the previous step. If so, select the projection point with the same direction angle as the result of the previous step.

[0107] 3) Based on the results of the previous two steps, a combination of projection point sequences for each station is generated for each route (since a station may have multiple possible projection points). The most reasonable combination is selected from these combinations. Reasonable means that the projection number of the station with the earlier number is lower than the projection number of the station with the later number.

[0108] 4) Recalculate the unreasonable site projection in the previous step: reduce the projection range and re-project.

[0109] 5) The length of the projection points of two adjacent stations along the line shape is the cross-sectional distance of the line.

[0110] The cross-sectional distance calculation system of the bus route embodiment of the present invention includes a data preparation module, a data standardization module, a direction angle calculation module, and a cross-sectional distance calculation module.

[0111] Data preparation module: collects line linear data and station data.

[0112] Data standardization module: standardizes the linear data of each line.

[0113] Direction angle calculation module: preliminarily calculates the station direction angle of each station on each line.

[0114] Section distance calculation module: Preliminary calculation of the projection points of each station on the line. If there are multiple projection points on the line for a station, the projection points are simplified according to the station direction angle of the station. The combinations of station projection information sequences on the line are listed to find the most likely combination of line station projection information sequences, and the unreliable projection information is recalculated. The station direction angle of each station on each line is updated, and then the section distance of the line segment is calculated. Among them, the most likely meaning is: the station projection sequence contains the largest number of ascending station projection numbers; the smaller of the two linear point numbers of the line segment where the projection point is located is the projection number.

[0115] As an embodiment, the data standardization module further performs coordinate system conversion on the linear data and the station data that are not in the same coordinate system, wherein the standardization includes:

[0116] a) deleting duplicate linear points or linear points whose distance is less than a first preset distance;

[0117] b) If the distance between the two points is greater than a second preset distance, linear points are added between the two points at equal intervals.

[0118] As an implementation method, the direction angle calculation module projects the site onto the line for each line and calculates the site direction angle of the site. If a site has multiple projection points on the line, the site is omitted, that is, only the site direction angle of the site that can be determined at this time is retained.

[0119] As an implementation method, the section distance calculation module calculates the section distance according to the following steps:

[0120] For each station on the line, traverse the line segments on the line to determine whether the station can be projected onto the line segment. If it can be projected and the projection distance is less than a third preset distance, record the projection information of the station: the corresponding projection point, projection number, and station direction angle; if not, take the linear point closest to the station as the projection point;

[0121] If there are multiple projection points on the line of the station, then select the projection point whose station direction angle is consistent with the station direction angle calculated by the direction angle calculation module from the multiple projection points;

[0122] Arrange the projection points of each station on the line in station order to form a line station projection information sequence;

[0123] Find the most likely combination of line station projection information sequences. If the projection number of a station is in the largest ascending projection number sequence, the projection information of the station is considered credible; otherwise, the projection information of the station is considered unreliable.

[0124] For the unreliable site projection information in the most likely line site projection information sequence combination, the projection range of the site on the line is reduced, the projection distance threshold is removed, and the projection point is recalculated.

[0125] If the recalculated station direction angle is inconsistent with the station direction angle calculated by the direction angle calculation module, the station direction angle of the station is updated according to the recalculated result;

[0126] The cross-sectional distances of the line are obtained by calculating the spherical distance between the projection points of two adjacent stations.

[0127] While embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions, and variations may be made to these embodiments without departing from the principles and spirit of the invention, and that the scope of the invention is defined by the appended claims and their equivalents.

Claims

1. A method for calculating the cross-sectional distance of a bus route, characterized in that: include: Step 1: Collect the line linear data and station data; Step 2: Standardize the linear data of each line; Step 3: Preliminary calculation of the station direction angle of each station on each line; Step 4: Preliminarily calculate the projection points of each station on the line. If a station has multiple projection points on the line, simplify the projection points based on the station's station bearing angle. List the combinations of station projection information sequences on the line, find the most likely combination of line station projection information sequences, and recalculate any unreliable projection information. Update the station bearing angle of each station on each line, and then calculate the line section distance. The most likely meaning is: the station projection sequence contains the largest number of ascending station projection numbers; the smaller of the two linear point numbers on the line segment where the projection point is located is the projection number. In step 3, for each line, perform a station-to-line projection and calculate the station direction angle of the station. If a station has multiple projection points on the line, the station is omitted. That is, only the station direction angle of the station that can be determined at this time is retained. Step 4 includes: Sub-step 41: For each station on the route, traverse the line segments on the line to determine whether the station can be projected onto the line segment. If projection is possible and the projection distance is less than a third preset distance, record the projection information of the station: the corresponding projection point, projection number, and station direction angle. If not, take the linear point closest to the station as the projection point. Sub-step 42: If the station has multiple projection points on the line shape, select a projection point whose station direction angle is consistent with the station direction angle calculated in step 3 from the multiple projection points; Sub-step 43: Arrange the projection points of each station on the line in station order to form a line station projection information sequence; Sub-step 44: Find the most likely combination of projection information sequences of route sites, wherein if the projection number of a site is in the maximum ascending projection number sequence, the projection information of the site is considered to be credible; otherwise, the projection information of the site is considered to be unreliable; Sub-step 45: for the unreliable site projection information in the most likely line site projection information sequence combination, narrow the projection range of the site on the line, remove the projection distance threshold, and then recalculate the projection point; Sub-step 46: If the recalculated station direction angle is inconsistent with the station direction angle calculated in step 3, then update the station direction angle of the station according to the recalculated result; Sub-step 47: Obtain the distances of each section of the line by calculating the spherical distance between the projection points of two adjacent stations.

2. The method for calculating the cross-sectional distance of a bus route according to claim 1, wherein: In step 2, coordinate system conversion is also performed on the data in the linear data and station data that are not in the same coordinate system. The standardization includes: a) deleting duplicate linear points or linear points whose distance is less than a first preset distance; b) If the distance between the two points is greater than a second preset distance, linear points are added between the two points at equal intervals.

3. A bus route section distance calculation system, characterized in that: include: Data preparation module: collects line linear data and station data; Data standardization module: standardizes the linear data of each line; Direction angle calculation module: preliminarily calculates the station direction angle of each station on each line; Section distance calculation module: Preliminary calculation of the projection points of each station on the line. If a station has multiple projection points on the line, the projection points are simplified based on the station bearing angle of the station. The combination of station projection information sequences on the line is listed, and the most likely combination of line station projection information sequences is found. Unreliable projection information is recalculated. The station bearing angle of each station on each line is updated, and then the section distance of the line segment is calculated. Among them, the most likely meaning is: the station projection sequence contains the largest number of ascending station projection numbers; the smaller of the two linear point numbers of the line segment where the projection point is located is the projection number. The direction angle calculation module projects the station onto the line for each line and calculates the station direction angle of the station. If a station has multiple projection points on the line, it will be omitted. That is, only the station direction angle of the station that can be determined at this time is retained. The section distance calculation module calculates the section distance according to the following steps: For each station on the line, traverse the line segments on the line to determine whether the station can be projected onto the line segment. If it can be projected and the projection distance is less than a third preset distance, record the projection information of the station: the corresponding projection point, projection number, and station direction angle; if not, take the linear point closest to the station as the projection point; If there are multiple projection points on the line of the station, then select the projection point whose station direction angle is consistent with the station direction angle calculated by the direction angle calculation module from the multiple projection points; Arrange the projection points of each station on the line in station order to form a line station projection information sequence; Find the most likely combination of line station projection information sequences. If the projection number of a station is in the largest ascending projection number sequence, the projection information of the station is considered credible; otherwise, the projection information of the station is considered unreliable. For the unreliable site projection information in the most likely line site projection information sequence combination, the projection range of the site on the line is narrowed, the projection distance threshold is removed, and the projection point is recalculated; If the recalculated station direction angle is inconsistent with the station direction angle calculated by the direction angle calculation module, the station direction angle of the station is updated according to the recalculated result; The cross-sectional distances of the line are obtained by calculating the spherical distance between the projection points of two adjacent stations.

4. The bus route section distance calculation system according to claim 3, characterized in that: The data standardization module also performs coordinate system conversion on linear data and station data that are not in the same coordinate system. Standardization includes: a) deleting duplicate linear points or linear points whose distance is less than a first preset distance; b) If the distance between the two points is greater than a second preset distance, linear points are added between the two points at equal intervals.

Images