A highway pile number latitude and longitude conversion method, system and storage medium
Patent Information
- Application Number
- CN202611293422.1
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2026-08-25
- Publication Date
- 2026-10-13
AI Technical Summary
但高速公路实际路线可能存在弯曲、转向、上下行分离以及道路路径与直线距离不一致等情况,若仅依据相邻整公里桩号之间的直线距离进行计算,容易导致目标桩号位置偏离真实道路路径,影响桩号位置信息的准确性
[0049]本发明提供的高速公路桩号经纬度转换方法、系统和存储介质通过获取高速公路道路网络数据,并结合起始桩号和终止桩号的经纬度坐标确定目标道路路径,能够使桩号经纬度转换过程更加贴合真实道路走向,通过将已知桩号投影至目标道路路径,并根据投影点之间的实际道路距离计算目标桩号位置,能够降低已知桩号与道路网络节点不完全重合对转换结果的影响,提高非整公里桩号经纬度计算的准确性,通过对目标道路路径进行细分并计算沿道路方向的累计距离,能够根据目标桩号的桩号比例关系自动插值得到对应经纬度坐标,提高高速公路桩号位置信息转换的自动化程度和计算效率,满足高速公路桩号位置信息精确化表达的需求。
Smart Images

Figure CN122903522A_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of transportation geographic information processing technology, specifically to a method, system, and storage medium for converting highway station numbers to latitude and longitude. Background Technology
[0002] Highway mileage markers are important markers indicating locations along the road, typically reflecting the relative distance of a target location within the highway's route. In practical road management and field applications, mileage markers are generally placed at full kilometer marks to indicate road mileage nodes.
[0003] With the development of digital management, refined maintenance, accident location, and road data analysis on highways, relying solely on whole-kilometer markers is no longer sufficient to meet the demand for accurate representation of road location information. In practical applications, it is often necessary to determine the spatial coordinates corresponding to a non-whole-kilometer marker, such as a target marker located between two whole-kilometer markers, or to accurately locate the point based on marker information in electronic maps, navigation systems, and geographic information systems. Existing methods typically rely on manual searching, map comparison, or simple interpolation based on the latitude and longitude of adjacent whole-kilometer markers. However, actual highway routes may have curves, turns, separation of upstream and downstream traffic, and discrepancies between the road path and the straight-line distance. If calculations are based solely on the straight-line distance between adjacent whole-kilometer markers, the target marker location may deviate from the actual road path, affecting the accuracy of the marker location information. Summary of the Invention
[0004] The purpose of this invention is to overcome the shortcomings of the prior art and provide a method, system and storage medium for converting highway station latitude and longitude, which can improve the automation and calculation efficiency of highway station location information conversion and meet the needs of accurate expression of highway station location information.
[0005] To achieve the above objectives, the present invention is implemented using the following technical solution:
[0006] In a first aspect, the present invention provides a method for converting highway station numbers to latitude and longitude, comprising the following steps:
[0007] Obtain road network data of the highway segment where the target station is located and two known station numbers and their latitude and longitude coordinates. The two known station numbers include a starting station number and an ending station number. The starting station number is the known station number that is the nearest whole kilometer ahead of the target station number along the road direction. The ending station number is the known station number that is the nearest whole kilometer ahead of the target station number along the road direction.
[0008] Based on the latitude and longitude coordinates of two known station numbers, candidate road nodes corresponding to the two known station numbers are determined in the road network data respectively;
[0009] Determine the target road path based on two known station numbers and their corresponding candidate road nodes;
[0010] Projecting two known station numbers onto the target road path yields the target road path projection point information for the two known station numbers.
[0011] The target road path is subdivided according to a preset interval to obtain several subdivided coordinate points distributed along the target road path, and the cumulative distance of each subdivided coordinate point along the road direction is calculated.
[0012] Based on the cumulative distance, the ratio between the target station and two known station, and the actual road distance between the projection points of the target road path of the two known station, the position of the target station on the target road path is determined, and the latitude and longitude coordinates of the target station are obtained by interpolation.
[0013] Furthermore, the road network data includes road nodes, road edges, latitude and longitude coordinates of road nodes, road edge length information, and road type information;
[0014] The road network data is obtained by expanding the latitude and longitude coordinates of the two known kilometer markers outward by 0.01° to 0.03° respectively, constructing a rectangular bounding box, using the rectangular bounding box as the acquisition range of the road network data, and filtering highway-type roads from the acquisition range to obtain the road network data.
[0015] Furthermore, based on the latitude and longitude coordinates of the two known station numbers, candidate road nodes corresponding to the two known station numbers are determined in the road network data, including the following steps:
[0016] Traverse each road node in the road network data, calculate the spherical distance between each road node and the latitude and longitude coordinates of any known station, and select a preset number of road nodes as candidate road nodes corresponding to this known station in ascending order of the spherical distance.
[0017] Furthermore, based on two known station numbers and their corresponding candidate road nodes, the target road path is determined, including the following steps:
[0018] The candidate road nodes corresponding to the starting chainage and the candidate road nodes corresponding to the ending chainage are combined in pairs to form several initial node pairs;
[0019] For any initial node pair, determine whether there is a walkable path between the initial node pair. If there is, obtain the shortest path between the initial node pair; otherwise, discard the initial node pair.
[0020] Project the two known station numbers onto the shortest path respectively to obtain the two shortest path projection points corresponding to the two known station numbers. Determine whether the shortest path contains the two shortest path projection points. If it does not contain them, discard the initial node pair; otherwise, retain them as candidate node pairs to obtain a set of candidate node pairs.
[0021] The evaluation function is used to comprehensively evaluate all candidate node pairs in the candidate node pair set to obtain the optimal candidate node pair. The shortest path between the optimal candidate node pairs is then determined as the target road path corresponding to the known station point.
[0022] The evaluation function includes a path length matching term and an endpoint location matching term. The path length matching term is used to characterize the difference between the shortest path length between the candidate node pairs to be evaluated and the target road segment length. The target road segment length is calculated based on the starting station number and the ending station number.
[0023] The endpoint position matching item is used to characterize the sum of the starting distance and the ending distance, wherein the starting distance is the distance between the shortest path start point and the shortest path projection point of the starting station, and the ending distance is the distance between the shortest path end point and the shortest path projection point of the ending station.
[0024] Based on the length difference and the total distance, calculate the comprehensive score of each candidate road path, and select the candidate road path with the smallest comprehensive score as the target road path;
[0025] The comprehensive score is expressed as follows: .
[0026] Furthermore, the two known station numbers are projected onto the target road path to obtain the target road path projection point information for the two known station numbers, including:
[0027] Project the starting station number to the starting segment of the target road path, and project the ending station number to the ending segment of the target road path.
[0028] Based on the projection positions of the starting station number on the starting line segment and the ending station number on the ending line segment, the latitude and longitude coordinates, starting offset distance, and ending offset distance of the target road path projection points for the two known station numbers are obtained respectively. The starting offset distance is the distance between the target road path projection point of the starting station number and the starting node of the target road path, and the ending offset distance is the distance between the target road path projection point of the ending station number and the ending node of the target road path.
[0029] Furthermore, the target road path is subdivided according to a preset interval to obtain several subdivided coordinate points distributed along the target road path, and the cumulative distance of each subdivided coordinate point along the road direction is calculated, including:
[0030] Extract the latitude and longitude coordinates of each road node in the target road path, and construct a path line object based on the latitude and longitude coordinates of each road node;
[0031] Interpolate along the path line object according to a preset interval to generate multiple subdivided coordinate points arranged sequentially according to the road extension direction;
[0032] According to the arrangement order of the subdivided coordinate points, the spherical distance between two adjacent subdivided coordinate points is calculated sequentially, and the spherical distance is accumulated segment by segment to obtain the cumulative distance along the road corresponding to each subdivided coordinate point.
[0033] Furthermore, the location of the target station on the target road path is determined and its latitude and longitude coordinates are obtained by interpolation, including:
[0034] Based on the ratio between the target station and two known station, and the actual road distance between the target road path projection points of the two known station, calculate the path distance from the target station's target road path projection point to the target station.
[0035] Add the distance along the path to the initial offset distance to obtain the target cumulative distance of the target station on the target road path;
[0036] Based on the cumulative distance of each subdivided coordinate point along the road direction, determine the two adjacent subdivided coordinate points where the target cumulative distance is located;
[0037] Calculate the proportion of the target's cumulative distance between two adjacent subdivided coordinate points;
[0038] Interpolation calculations are performed based on the distance ratio to obtain the latitude and longitude coordinates corresponding to the target station number.
[0039] Furthermore, the actual road distance between the known station projection points is calculated as follows:
[0040] Calculate the cumulative distance of the target road path from the starting node to the ending node;
[0041] Subtract the starting offset distance and the ending offset distance from the cumulative distance to obtain the actual road distance between the two known station number target road path projection points along the target road path direction.
[0042] Secondly, the present invention provides a highway station latitude and longitude conversion system for implementing the aforementioned highway station latitude and longitude conversion method, comprising:
[0043] The road network extraction module is used to obtain the road network data of the highway section where the target station is located, as well as two known station numbers and their latitude and longitude coordinates.
[0044] The node matching module is used to determine candidate road nodes corresponding to the two known station numbers in the road network data based on the latitude and longitude coordinates of the two known station numbers, and to determine the target road path based on the two known station numbers and their corresponding candidate road nodes.
[0045] The station projection and road segment subdivision module is used to project two known station numbers onto the target road path to obtain the projection point information of the target road path for the two known station numbers, and to subdivide the target road path according to a preset interval to obtain a number of subdivided coordinate points distributed along the target road path, and to calculate the cumulative distance of each subdivided coordinate point along the road direction.
[0046] The coordinate transformation module is used to determine the position of the target station on the target road path based on the cumulative distance, the station ratio between the target station and two known station, and the actual road distance between the projection points of the target road path of the two known station, and to interpolate to obtain the latitude and longitude coordinates of the target station.
[0047] Thirdly, the present invention provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the highway station latitude and longitude conversion method.
[0048] Compared with the prior art, the beneficial effects achieved by the present invention are as follows:
[0049] The highway station latitude and longitude conversion method, system, and storage medium provided by this invention obtain highway road network data and determine the target road path by combining the latitude and longitude coordinates of the starting and ending station numbers. This makes the station latitude and longitude conversion process more closely resemble the actual road direction. By projecting known station numbers onto the target road path and calculating the target station number position based on the actual road distance between the projection points, the impact of incomplete overlap between known station numbers and road network nodes on the conversion result can be reduced, improving the accuracy of latitude and longitude calculation for non-whole kilometer station numbers. By subdividing the target road path and calculating the cumulative distance along the road direction, the corresponding latitude and longitude coordinates can be automatically interpolated based on the station number ratio of the target station numbers, improving the automation and calculation efficiency of highway station number position information conversion and meeting the need for accurate expression of highway station number position information. Attached Figure Description
[0050] Figure 1 A flowchart illustrating the highway station latitude and longitude conversion method provided in this embodiment of the invention;
[0051] Figure 2 This is a schematic diagram of the road path projected from a known station number in an embodiment of the present invention;
[0052] Figure 3 This is a schematic diagram illustrating the process of calculating the latitude and longitude of the target station in an embodiment of the present invention;
[0053] Figure 4 This is a schematic diagram of the structure of the highway station latitude and longitude conversion system provided in an embodiment of the present invention. Detailed Implementation
[0054] The present invention will be further described below with reference to the accompanying drawings. The following embodiments are only used to more clearly illustrate the technical solution of the present invention, and should not be used to limit the scope of protection of the present invention.
[0055] Example 1
[0056] This embodiment provides a method for converting highway station numbers to latitude and longitude, such as... Figure 1 As shown, it includes the following steps:
[0057] S1. Obtain the road network data of the highway section where the target station is located and two known station numbers and their latitude and longitude coordinates. The two known station numbers include a starting station number and an ending station number. The starting station number is the known station number that is the nearest whole kilometer ahead of the target station number along the road direction. The ending station number is the known station number that is the nearest whole kilometer ahead of the target station number along the road direction.
[0058] S2. Based on the latitude and longitude coordinates of the two known station numbers, determine the candidate road nodes corresponding to the two known station numbers in the road network data respectively;
[0059] S3. Determine the target road path based on two known station numbers and their corresponding candidate road nodes;
[0060] S4. Project the two known station numbers onto the target road path to obtain the target road path projection point information for the two known station numbers;
[0061] S5. Subdivide the target road path according to a preset interval to obtain several subdivided coordinate points distributed along the target road path, and calculate the cumulative distance of each subdivided coordinate point along the road direction.
[0062] S6. Based on the cumulative distance, the ratio between the target station and the two known station, and the actual road distance between the projection points of the target road path of the two known station, determine the position of the target station on the target road path, and interpolate to obtain the latitude and longitude coordinates of the target station.
[0063] First, obtain the road network data of the area where the target highway is located, as well as the station numbers and latitude and longitude coordinates of at least two known station numbers. The road network data includes road nodes, road edges, node latitude and longitude coordinates, and road edge length information. In this embodiment, the known station numbers are the starting station number K760 and the ending station number K762, but can also be other adjacent or non-adjacent known station numbers.
[0064] After acquiring the road network data, based on the latitude and longitude coordinates of each known station, multiple candidate road nodes that are relatively close to that known station are selected from the road network data. Specifically, each road node in the road network data is traversed, and the spherical distance between each road node and the latitude and longitude coordinates of any known station is calculated. A preset number of road nodes are selected as candidate road nodes corresponding to this known station, in ascending order of the spherical distance. The preset number can be set according to the accuracy requirements, such as 5, 10, or other numbers.
[0065] Subsequently, the candidate road nodes corresponding to the starting chainage and the candidate road nodes corresponding to the ending chainage are combined in pairs to form several initial node pairs.
[0066] For any initial node pair, determine whether there is a walkable path between the initial node pair. If there is, obtain the shortest path between the initial node pair; otherwise, discard the initial node pair.
[0067] Project the two known station numbers onto the shortest path respectively to obtain the two shortest path projection points corresponding to the two known station numbers. Determine whether the shortest path contains the two shortest path projection points. If it does not contain them, discard the initial node pair; otherwise, retain them as candidate node pairs to obtain a set of candidate node pairs.
[0068] The evaluation function is used to comprehensively evaluate all candidate node pairs in the candidate node pair set to obtain the optimal candidate node pair. The shortest path between the optimal candidate node pairs is then determined as the target road path corresponding to the known station point.
[0069] The evaluation function includes a path length matching term and an endpoint location matching term. The path length matching term is used to characterize the difference between the shortest path length between the candidate node pairs to be evaluated and the target road segment length. The target road segment length is calculated based on the starting station number and the ending station number.
[0070] The endpoint position matching item is used to characterize the sum of the starting distance and the ending distance, wherein the starting distance is the distance between the shortest path start point and the shortest path projection point of the starting station, and the ending distance is the distance between the shortest path end point and the shortest path projection point of the ending station.
[0071] Based on the length difference and the total distance, calculate the comprehensive score of each candidate road path, and select the candidate road path with the smallest comprehensive score as the target road path.
[0072] The optimal candidate node pair is used to determine the starting and ending nodes of the target road path. In this embodiment, the comprehensive score can be expressed as:
[0073]
[0074] like Figure 2 As shown, a known station number does not necessarily fall exactly on a road network node; therefore, it is necessary to project the known station number onto the target road path. Projecting two known station numbers onto the target road path yields the target road path projection point information for the two known station numbers, including:
[0075] Project the starting station number to the starting segment of the target road path, and project the ending station number to the ending segment of the target road path.
[0076] Based on the projection positions of the starting station number on the starting line segment and the ending station number on the ending line segment, the latitude and longitude coordinates, starting offset distance, and ending offset distance of the target road path projection points for the two known station numbers are obtained respectively. The starting offset distance is the distance between the target road path projection point of the starting station number and the starting node of the target road path, and the ending offset distance is the distance between the target road path projection point of the ending station number and the ending node of the target road path.
[0077] After determining the target road path and the projection points of the target road path at two known station numbers, the target road path is subdivided according to a preset interval to obtain several subdivided coordinate points distributed along the target road path. The preset interval can be set according to the conversion accuracy requirements, such as 5 meters, 10 meters, or other distances.
[0078] For each road subdivision coordinate point, calculate its cumulative distance relative to the starting road node along the target road path direction, including:
[0079] Extract the latitude and longitude coordinates of each road node in the target road path, and construct a path line object based on the latitude and longitude coordinates of each road node;
[0080] Interpolate along the path line object according to a preset interval to generate multiple subdivided coordinate points arranged sequentially according to the road extension direction;
[0081] According to the arrangement order of the subdivided coordinate points, the spherical distance between two adjacent subdivided coordinate points is calculated sequentially, and the spherical distance is accumulated segment by segment to obtain the cumulative distance along the road corresponding to each subdivided coordinate point.
[0082] like Figure 3 As shown, when it is necessary to determine the location of the target station on the target road path and interpolate to obtain the latitude and longitude coordinates of the target station, the following steps are included:
[0083] Based on the ratio between the target station and two known station, and the actual road distance between the target road path projection points of the two known station, calculate the path distance from the target station's target road path projection point to the target station.
[0084] Add the distance along the path to the initial offset distance to obtain the target cumulative distance of the target station on the target road path;
[0085] Based on the cumulative distance of each subdivided coordinate point along the road direction, determine the two adjacent subdivided coordinate points where the target cumulative distance is located;
[0086] Calculate the proportion of the target's cumulative distance between two adjacent subdivided coordinate points;
[0087] Interpolation calculations are performed based on the distance ratio to obtain the latitude and longitude coordinates corresponding to the target station number.
[0088] In this embodiment, the positional ratio of the target station to the two known stations is first calculated based on the station ratio between the target station and the two known stations. The station ratio can be expressed as:
[0089]
[0090] Then, based on the station ratio and the actual road length between the two known station projection points, the distance along the road path from the starting station projection point to the target station is calculated. Since the cumulative distance of the road subdivision coordinate points is calculated from the starting node of the target road path, it is also necessary to add the distance between the starting road node and the starting station projection point to the distance between the starting station projection point and the target station to obtain the target cumulative distance of the target station relative to the starting road node.
[0091] After obtaining the target cumulative distance, the location of the target cumulative distance is found in the cumulative distance sequence corresponding to the road subdivision coordinate points, determining which two adjacent road subdivision coordinate points the target station falls between. If the target cumulative distance is located at the... The road subdivision coordinates and the first Between each road subdivision coordinate point, the proportion of the cumulative distance to the target within that adjacent interval is calculated. Based on this proportion, the latitude and longitude coordinates of two adjacent road subdivision coordinate points are interpolated to obtain the latitude and longitude coordinates corresponding to the target station number. The positional order of the coordinate points for road subdivision.
[0092] Specifically, if the The latitude and longitude coordinates of each road subdivision coordinate point are ( , ), No. The latitude and longitude coordinates of each road subdivision coordinate point are ( , The ratio of the cumulative distance of the target to the two is: Then the latitude and longitude of the target station can be calculated as follows:
[0093]
[0094]
[0095] The above calculations yield the latitude and longitude coordinates corresponding to the target station. This calculation method does not directly interpolate a straight line between two known station coordinates; instead, it interpolates the distance along the road based on the actual road path, road nodes, station projection points, and subdivided road coordinate points. Therefore, it is better suited to actual road scenarios involving curved sections, bends, and numerous nodes on highways.
[0096] In one specific embodiment, given the latitude and longitude coordinates of the starting station K760 and the ending station K762, it is necessary to calculate the latitude and longitude coordinates of the target station K760.5, which includes the following steps:
[0097] The road network acquisition range is determined based on the latitude and longitude coordinates of K760 and K762, and the highway road network data within this range is obtained.
[0098] Candidate road nodes near K760 and K762 are filtered separately, and target road paths that can include the projection points of K760 and K762 are determined.
[0099] Project K760 and K762 onto the target road path respectively, and calculate the actual road length between their projection points;
[0100] Based on the station ratio between K760.5 and K760 and K762, the road distance B of K760.5 relative to the projection point of K760 is calculated. Combined with the distance A from the starting road node to the projection point of K760, the cumulative distance C of K760.5 relative to the starting road node is determined. Then, the latitude and longitude coordinates of K760.5 are obtained by interpolation of the road subdivision coordinate points.
[0101] Example 2
[0102] This embodiment provides a highway station latitude and longitude conversion system to implement the highway station latitude and longitude conversion method described in Embodiment 1, such as... Figure 4As shown, the system includes a road network extraction module, a node matching module, a stationing projection and road segment subdivision module, and a coordinate transformation module.
[0103] Among them, the road network extraction module is used to obtain the road network data of the highway section where the target station is located and two known station numbers and their latitude and longitude coordinates;
[0104] The node matching module is used to determine candidate road nodes corresponding to the two known station numbers in the road network data based on the latitude and longitude coordinates of the two known station numbers, and to determine the target road path based on the two known station numbers and their corresponding candidate road nodes.
[0105] The station projection and road segment subdivision module is used to project two known station numbers onto the target road path to obtain the projection point information of the target road path for the two known station numbers, and to subdivide the target road path according to a preset interval to obtain a number of subdivided coordinate points distributed along the target road path, and to calculate the cumulative distance of each subdivided coordinate point along the road direction.
[0106] The coordinate transformation module is used to determine the position of the target station on the target road path based on the cumulative distance, the station ratio between the target station and two known station, and the actual road distance between the projection points of the target road path of the two known station, and to interpolate to obtain the latitude and longitude coordinates of the target station.
[0107] The methods and systems described in this invention can be implemented using Python programs, or other computer programming languages or geographic information processing platforms. Any system capable of performing the technical processes of road network acquisition, node matching, station projection, path subdivision, cumulative distance calculation, and latitude / longitude interpolation conversion falls within the scope of protection of this invention.
[0108] Example 3
[0109] This embodiment provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the highway station latitude and longitude conversion method as described in Embodiment 1.
[0110] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.
[0111] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.
[0112] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.
[0113] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.
[0114] The embodiments of the present invention have been described above with reference to the accompanying drawings. However, the present invention is not limited to the specific embodiments described above. The specific embodiments described above are merely illustrative and not restrictive. Those skilled in the art can make many other forms under the guidance of the present invention without departing from the spirit and scope of the claims. All of these forms are within the protection scope of the present invention.
Claims
1. A method for converting highway station numbers to latitude and longitude, characterized in that, Includes the following steps: Obtain road network data of the highway segment where the target station is located and two known station numbers and their latitude and longitude coordinates. The two known station numbers include a starting station number and an ending station number. The starting station number is the known station number that is the nearest whole kilometer ahead of the target station number along the road direction. The ending station number is the known station number that is the nearest whole kilometer ahead of the target station number along the road direction. Based on the latitude and longitude coordinates of two known station numbers, candidate road nodes corresponding to the two known station numbers are determined in the road network data respectively; Determine the target road path based on two known station numbers and their corresponding candidate road nodes; Projecting two known station numbers onto the target road path yields the target road path projection point information for the two known station numbers. The target road path is subdivided according to a preset interval to obtain several subdivided coordinate points distributed along the target road path, and the cumulative distance of each subdivided coordinate point along the road direction is calculated. Based on the cumulative distance, the ratio between the target station and two known station, and the actual road distance between the projection points of the target road path of the two known station, the position of the target station on the target road path is determined, and the latitude and longitude coordinates of the target station are obtained by interpolation.
2. The method for converting highway station latitude and longitude according to claim 1, characterized in that: The road network data includes road nodes, road edges, latitude and longitude coordinates of road nodes, road edge length information, and road type information; The road network data is obtained by expanding the latitude and longitude coordinates of the two known kilometer markers outward by 0.01° to 0.03° respectively, constructing a rectangular bounding box, using the rectangular bounding box as the acquisition range of the road network data, and filtering highway-type roads from the acquisition range to obtain the road network data.
3. The method for converting highway station latitude and longitude according to claim 2, characterized in that, Based on the latitude and longitude coordinates of two known station numbers, candidate road nodes corresponding to the two known station numbers are determined in the road network data, including the following steps: Traverse each road node in the road network data, calculate the spherical distance between each road node and the latitude and longitude coordinates of any known station, and select a preset number of road nodes as candidate road nodes corresponding to this known station in ascending order of the spherical distance.
4. The method for converting highway station latitude and longitude according to claim 1, characterized in that, Determining the target road path based on two known station numbers and their corresponding candidate road nodes includes the following steps: The candidate road nodes corresponding to the starting chainage and the candidate road nodes corresponding to the ending chainage are combined in pairs to form several initial node pairs; For any initial node pair, determine whether there is a walkable path between the initial node pair. If there is, obtain the shortest path between the initial node pair; otherwise, discard the initial node pair. Project the two known station numbers onto the shortest path respectively to obtain the two shortest path projection points corresponding to the two known station numbers. Determine whether the shortest path contains the two shortest path projection points. If it does not contain them, discard the initial node pair; otherwise, retain them as candidate node pairs to obtain a set of candidate node pairs. The evaluation function is used to comprehensively evaluate all candidate node pairs in the candidate node pair set to obtain the optimal candidate node pair. The shortest path between the optimal candidate node pairs is then determined as the target road path corresponding to the known station point. The evaluation function includes a path length matching term and an endpoint location matching term. The path length matching term is used to characterize the difference between the shortest path length between the candidate node pairs to be evaluated and the target road segment length. The target road segment length is calculated based on the starting station number and the ending station number. The endpoint position matching item is used to characterize the sum of the starting distance and the ending distance, wherein the starting distance is the distance between the shortest path start point and the shortest path projection point of the starting station, and the ending distance is the distance between the shortest path end point and the shortest path projection point of the ending station. Based on the length difference and the total distance, calculate the comprehensive score of each candidate road path, and select the candidate road path with the smallest comprehensive score as the target road path; The comprehensive score is expressed as follows: .
5. The method for converting highway station latitude and longitude according to claim 1, characterized in that, Projecting two known station numbers onto the target road path yields projection point information for the target road path with the two known station numbers, including: Project the starting station number to the starting segment of the target road path, and project the ending station number to the ending segment of the target road path. Based on the projection positions of the starting station number on the starting line segment and the ending station number on the ending line segment, the latitude and longitude coordinates, starting offset distance, and ending offset distance of the target road path projection points for the two known station numbers are obtained respectively. The starting offset distance is the distance between the target road path projection point of the starting station number and the starting node of the target road path, and the ending offset distance is the distance between the target road path projection point of the ending station number and the ending node of the target road path.
6. The method for converting highway station latitude and longitude according to claim 1, characterized in that, The target road path is subdivided according to a preset interval to obtain several subdivided coordinate points distributed along the target road path, and the cumulative distance of each subdivided coordinate point along the road direction is calculated, including: Extract the latitude and longitude coordinates of each road node in the target road path, and construct a path line object based on the latitude and longitude coordinates of each road node; Interpolate along the path line object according to a preset interval to generate multiple subdivided coordinate points arranged sequentially according to the road extension direction; According to the arrangement order of the subdivided coordinate points, the spherical distance between two adjacent subdivided coordinate points is calculated sequentially, and the spherical distance is accumulated segment by segment to obtain the cumulative distance along the road corresponding to each subdivided coordinate point.
7. The method for converting highway station latitude and longitude according to claim 5, characterized in that, Determine the location of the target station on the target road path and interpolate to obtain the latitude and longitude coordinates of the target station, including: Based on the ratio between the target station and two known station, and the actual road distance between the target road path projection points of the two known station, calculate the path distance from the target station's target road path projection point to the target station. Add the distance along the path to the initial offset distance to obtain the target cumulative distance of the target station on the target road path; Based on the cumulative distance of each subdivided coordinate point along the road direction, determine the two adjacent subdivided coordinate points where the target cumulative distance is located; Calculate the proportion of the target's cumulative distance between two adjacent subdivided coordinate points; Interpolation calculations are performed based on the distance ratio to obtain the latitude and longitude coordinates corresponding to the target station number.
8. The method for converting highway station latitude and longitude according to claim 7, characterized in that, The actual road distance between the known station projection points is calculated as follows: Calculate the cumulative distance of the target road path from the starting node to the ending node; Subtract the starting offset distance and the ending offset distance from the cumulative distance to obtain the actual road distance between the two known station number target road path projection points along the target road path direction.
9. A highway station latitude and longitude conversion system, used to implement the highway station latitude and longitude conversion method as described in any one of claims 1-8, characterized in that, include: The road network extraction module is used to obtain the road network data of the highway section where the target station is located, as well as two known station numbers and their latitude and longitude coordinates. The node matching module is used to determine candidate road nodes corresponding to the two known station numbers in the road network data based on the latitude and longitude coordinates of the two known station numbers, and to determine the target road path based on the two known station numbers and their corresponding candidate road nodes. The station projection and road segment subdivision module is used to project two known station numbers onto the target road path to obtain the projection point information of the target road path for the two known station numbers, and to subdivide the target road path according to a preset interval to obtain a number of subdivided coordinate points distributed along the target road path, and to calculate the cumulative distance of each subdivided coordinate point along the road direction. The coordinate transformation module is used to determine the position of the target station on the target road path based on the cumulative distance, the station ratio between the target station and two known station, and the actual road distance between the projection points of the target road path of the two known station, and to interpolate to obtain the latitude and longitude coordinates of the target station.
10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, implements the highway station latitude and longitude conversion method as described in any one of claims 1-8.