Interprovincial road section automatic identification method and system based on geometric space calculation

Abstract

The invention discloses an inter-provincial road section automatic identification method and system based on geometric space calculation. The method comprises the following steps: acquiring space geometric data of a road network and a provincial administrative division; constructing a first geometric object corresponding to the road section and a second geometric object corresponding to the provincial boundary based on the spatial geometric data; calculating a spatial relationship between the first geometric object and the second geometric object, and screening candidate road sections intersecting with the provincial boundary based on the spatial relationship; performing expansion processing in the road network based on the candidate road section to obtain a target road section; and generating planar inter-provincial road section area data based on the target road section. According to the method, automatic identification of the inter-provincial road sections is realized, the technical problem that the public security and traffic management department cannot automatically calibrate cross-provincial traffic events due to lack of an inter-provincial road section space information base is solved, and the informatization and intelligence level of traffic scheduling is effectively improved.

Description

Technical Field

[0001] This invention belongs to the interdisciplinary field of intelligent traffic management and spatial information technology, and in particular relates to an automatic identification method and system for inter-provincial road sections based on geometric spatial calculation. Background Technology

[0002] Inter-provincial road segment incident management is a core business of public security traffic management departments. However, at present, public security organs lack a nationwide database of inter-provincial road segment information. They are unable to automatically identify inter-provincial areas from the data level of currently used maps such as Baidu and Gaode. They are also unable to automatically identify whether a specific event such as congestion, severe weather, or abnormal parking belongs to an inter-provincial area based on specific event information. This affects the national dispatch informatization and intelligentization capabilities of public security traffic management departments. Summary of the Invention

[0003] To address the aforementioned technical problems, this invention proposes an automatic identification method and system for inter-provincial road sections based on geometric space calculation, thereby resolving the issues present in the prior art.

[0004] To achieve the above objectives, the present invention provides an automatic identification method for inter-provincial road sections based on geometric space calculation, comprising:

[0005] Acquire spatial geometric data of road networks and provincial administrative divisions;

[0006] Based on the spatial geometric data, a first geometric object corresponding to the road segment and a second geometric object corresponding to the provincial boundary are constructed.

[0007] Calculate the spatial relationship between the first geometric object and the second geometric object, and filter candidate road segments that intersect with the provincial boundary based on the spatial relationship;

[0008] Based on the candidate road segments, the road network is expanded to obtain the target road segments.

[0009] Based on the target road segment, area data of inter-provincial road segments is generated in a planar format.

[0010] Optionally, the process of acquiring spatial geometric data of road networks and provincial administrative divisions includes:

[0011] Obtain the original road network trajectory coordinate data and provincial administrative division boundary coordinate data;

[0012] The road network trajectory coordinate data and provincial administrative division boundary coordinate data are batch converted into unified target coordinate system data.

[0013] Optionally, the process of constructing the first geometric object corresponding to the road segment and the second geometric object corresponding to the provincial boundary includes:

[0014] The road network trajectory coordinate data after unifying the coordinate system is deserialized using WKT and constructed as a linear geometric object, which serves as the first geometric object.

[0015] The coordinate data of provincial administrative division boundaries after unifying the coordinate system are deserialized using WKT and constructed into polygonal geometric objects, which serve as the second geometric object.

[0016] Optionally, calculating the spatial relationship between the first geometric object and the second geometric object includes:

[0017] Extract the extreme values ​​of all coordinate points of the first geometric object and construct the minimum bounding rectangle; determine the spatial overlap relationship between the minimum bounding rectangle and the second geometric object; when the spatial overlap relationship exists, perform a precise spatial relationship calculation between the first geometric object and the second geometric object;

[0018] When the result of the precise spatial relationship calculation is an intersection, the road segment represented by the corresponding first geometric object is identified as the candidate road segment.

[0019] Optionally, the precise intersection relationship between the first geometric object and the second geometric object can be calculated using a geometric object intersection function.

[0020] Optionally, the process of expanding the road network based on the candidate road segments includes:

[0021] Based on the topological connection relationship of the road network, the direct upstream adjacent road segment and the direct downstream adjacent road segment of the candidate road segment are obtained; the candidate road segment is connected with the direct upstream adjacent road segment and the direct downstream adjacent road segment to obtain a continuous extended road segment, which is used as the target road segment.

[0022] Optionally, based on the target road segment, area-based inter-provincial road segment data is generated, including:

[0023] The geometric objects of the target road segment are converted into a projected coordinate system suitable for linear distance calculation;

[0024] A buffer analysis operation is performed on the geometric object in the projected coordinate system to generate planar inter-provincial road segment area data;

[0025] Convert the area data of inter-provincial road sections in planar form back to the target coordinate system.

[0026] This invention also provides an automatic identification system for inter-provincial road sections based on geometric space calculation, used to implement the above method, comprising:

[0027] The data acquisition module is used to acquire spatial geometric data of road networks and provincial administrative divisions;

[0028] A geometric object construction module, connected to the data acquisition module, is used to construct a first geometric object corresponding to a road segment and a second geometric object corresponding to a provincial boundary based on the spatial geometric data.

[0029] A spatial relationship calculation module, connected to the geometric object construction module, is used to calculate the spatial relationship between the first geometric object and the second geometric object;

[0030] The candidate road segment screening module is connected to the spatial relationship calculation module and is used to screen candidate road segments that intersect with the provincial boundary based on the spatial relationship.

[0031] A road segment expansion module, connected to the candidate road segment filtering module, is used to perform expansion processing on the road network based on the candidate road segments to obtain the target road segment;

[0032] The region generation module, connected to the road segment expansion module, is used to generate area-based inter-provincial road segment region data based on the target road segment.

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

[0034] This invention, based on the geographical coordinate trajectory information of the national road network and publicly available national administrative division information, employs WKT deserialization processing on the road network trajectory information to construct spatial data containing geometric objects from each coordinate of the road network. It also organizes publicly available administrative division information, and, based on a unified coordinate system standard, determines the spatial intersection relationship between road segment information and provincial boundaries item by item. For intersecting road segments, it expands the selection of certain upstream and downstream areas, thereby expanding short-distance road segments into longer-distance inter-provincial road segment information usable by public security traffic management departments. Simultaneously, to ensure the reliability of matching various inter-provincial events, it employs geometric expansion on the linear data, expanding road segments into regions by setting custom buffers on both sides of the linear data. This data constructs an automated national inter-provincial road segment database to support public security traffic management work. Attached Figure Description

[0035] The accompanying drawings, which form part of this application, are used to provide a further understanding of this application. The illustrative embodiments and descriptions of this application are used to explain this application and do not constitute an undue limitation of this application. In the drawings:

[0036] Figure 1 This is a flowchart of a method according to an embodiment of the present invention;

[0037] Figure 2 This is a flowchart of inter-provincial road segment spatial identification based on geometric spatial calculation, according to an embodiment of the present invention. Detailed Implementation

[0038] It should be noted that, unless otherwise specified, the embodiments and features described in this application can be combined with each other. This application will now be described in detail with reference to the accompanying drawings and embodiments.

[0039] It should be noted that the steps shown in the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although a logical order is shown in the flowchart, in some cases the steps shown or described may be executed in a different order than that shown here.

[0040] Example 1

[0041] like Figure 1 As shown, this embodiment provides an automatic identification method for inter-provincial road segments based on geometric space calculation, including:

[0042] Using a dedicated program, administrative division information provided by public interfaces of map providers such as Baidu and Gaode is downloaded and compiled into a list of information for each province, which can be stored in CSV format. The boundary information of all provinces is then compiled and summarized for later use, thus obtaining a set of left-hand boundary points for all provinces in China.

[0043] The obtained national expressway network data includes the coordinates of points along each kilometer of the highway. To subsequently calculate inter-provincial road sections, a unified coordinate system is required.

[0044] The latitude and longitude trajectories of each road link (road segment connection) in the road network data are processed using WKT deserialization to construct linear geometric objects. Currently, in various types of raw road network data, although each road is divided into different road links, and each link reflects a corresponding latitude and longitude trajectory, this data is only string information containing latitude and longitude data. To facilitate subsequent spatial calculations, the raw latitude and longitude information must first be converted into corresponding geometric objects at the data level. Therefore, the following processing is performed:

[0045] (1) Standardize the coordinate system. It is necessary to verify the coordinate system type of the current road network information. Since road network data can be downloaded from the public internet, but the relevant latitude and longitude types are not labeled, it is necessary to manually verify the coordinate system type of the road network information first. This is done by selecting points on the map and projecting them one by one to verify whether the coordinate system is WGS84, GCJ02, or BD09. After determining the coordinate system type, based on the map to be used, determine the coordinate system type to be converted, and convert all coordinate points to the corresponding coordinate system type.

[0046] (2) Extract the latitude and longitude information of each link and convert it into a geometric object. After the coordinate system transformation is completed, the corresponding new coordinates can be stored in a CSV or database file, but the relevant coordinate strings are all of string type and cannot be used in subsequent spatial calculations. Therefore, the WKT deserialization method is adopted to construct the relevant data into a linear geometric object through software program.

[0047] Since highway and inter-provincial data are discrete points at the data level, making it difficult to judge solely from the data itself, this embodiment uses WKT deserialization to construct each data point into a line segment. The WKT deserialization method uses the shapely tool, calling common methods developed during program development.

[0048] The boundary data of each provincial administrative division is processed using WKT deserialization. Provincial administrative divisions can be obtained from the internet through the external interface of the corresponding map company, or they can be compiled manually. However, it is essential to ensure that the coordinate system type is consistent with the road network. After obtaining the provincial administrative divisions, the boundary coordinate string data of each province is deserialized using WKT to obtain the corresponding polygonal geometric object of that provincial division.

[0049] The spatial topological relationships between road segments and provincial-level geometric objects are calculated one by one. The intersection of lines and polygons is used as a characteristic; if they intersect, it proves that the link of that road segment crosses a provincial boundary, and therefore must be an inter-provincial road link. This allows for further expansion. For example... Figure 2 As shown, the specific steps include:

[0050] (1) A rough calculation is performed using an outer rectangle filtering method. Before performing the intersection judgment between geometric objects, the four maximum values ​​of minx, miny, maxx, and maxy of the link geometric object of each road segment are obtained from all coordinate points to construct a minimum bounding rectangle that can contain the geometric object. By comparing the boundary of the rectangle with the polygonal regions of different provinces, if there is no overlap, the precise geometric calculation is skipped, thereby improving the overall processing efficiency; if there is overlap, the precise calculation is further performed. This improves the overall system computational efficiency.

[0051] (2) Determine the positional relationship between linear and polygonal geometric objects based on whether they intersect. Calculate the relative positional relationship between the link geometric objects of the road segment and the polygonal geometric objects of the provincial boundary, and calculate whether a line-plane intersection relationship exists for each.

[0052] Due to the massive amount of data in the national road network and inter-provincial data, this embodiment uses inter-provincial regions as the basis for comparison to improve efficiency. Each highway alignment data point is compared with inter-provincial road segments to determine if they intersect. The intersection algorithm can be called using the `intersection` function of the `shapely` utility.

[0053] The `intersection` function is typically used; it outputs True as long as the two lines intersect.

[0054] To address the challenge of excessive data volume, especially since much of the highway data does not intersect with provincial boundaries, this embodiment considers first performing spatial polygon intersection calculations. This involves determining whether the current highway alignment data spatially intersects with the provincial boundary of the province being compared. If they do not intersect, it indicates the current segment is outside the provincial boundary, and the process is skipped, moving on to another segment. If they intersect, it proves the segment is on or inside the current provincial boundary, allowing for further verification of whether the road segment's alignment intersects with inter-provincial alignment data. This significantly reduces the resource consumption for spatial computation.

[0055] The specific method involves constructing the current line segment into a minimum rectangle. The extreme coordinates of the line segment are extracted, and the corresponding rectangle is constructed. The polygon of this rectangle and the intersecting polygons are then used to calculate whether they intersect.

[0056] Finally, by iterating through the data, all results where the intersection of the road segment alignment and the inter-provincial alignment is True can be stored in the database for later use.

[0057] Extend the road segment link both upstream and downstream. Based on the upstream and downstream data of the link, extend it to obtain a continuous road segment area of ​​a certain length. Specifically, based on the basic road network information, find the adjacent upstream and downstream data of the link (adjacent upstream and downstream data are the upstream and downstream mileage data arranged in order of road name and kilometer marker in the highway information table database). Obtain the road segment information of the required extension length by calculating distance or field annotation. If kilometer marker data is available, use the kilometer marker data for judgment; otherwise, calculate based on latitude and longitude.

[0058] Furthermore, based on the obtained linear coordinates, a buffer value (floating data, which can be set according to actual needs) is set, and the coordinate line is shifted up and down by a custom distance to obtain the expanded inter-provincial road segment area coordinate sequence.

[0059] Geometric expansion calculations are performed on road segment information to obtain inter-provincial road segment areas. First, the coordinates of the road segments are uniformly converted to the WGS84 coordinate system, and then further converted to Web Mercator coordinates using projection to support accurate distance calculations. Based on this, a buffer expansion operation is performed on the line-type object, expanding the geometric object of the line type to both sides (horizontally or vertically by the buffer value) to generate a surface-type geometric object. Finally, the coordinates of the expanded object are converted back to the WGS84 coordinate system and ultimately to the required coordinate system type. Then, using this line segment as the center, the buffer value is expanded horizontally or vertically by the specified distance.

[0060] This embodiment also provides an automatic inter-provincial road segment identification system based on geometric space calculation, used to implement the above method, including:

[0061] The data acquisition module is used to acquire spatial geometric data of road networks and provincial administrative divisions;

[0062] A geometric object construction module, connected to the data acquisition module, is used to construct a first geometric object corresponding to a road segment and a second geometric object corresponding to a provincial boundary based on the spatial geometric data.

[0063] A spatial relationship calculation module, connected to the geometric object construction module, is used to calculate the spatial relationship between the first geometric object and the second geometric object;

[0064] The candidate road segment screening module is connected to the spatial relationship calculation module and is used to screen candidate road segments that intersect with the provincial boundary based on the spatial relationship.

[0065] A road segment expansion module, connected to the candidate road segment filtering module, is used to perform expansion processing on the road network based on the candidate road segments to obtain the target road segment;

[0066] The region generation module, connected to the road segment expansion module, is used to generate area-based inter-provincial road segment region data based on the target road segment.

[0067] The above are merely preferred embodiments of this application, but the scope of protection of this application is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the scope of the technology disclosed in this application should be included within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A method for automatic identification of inter-provincial road sections based on geometric spatial calculation, characterized in that, Includes the following steps: Acquire spatial geometric data of road networks and provincial administrative divisions; Based on the spatial geometric data, a first geometric object corresponding to the road segment and a second geometric object corresponding to the provincial boundary are constructed. Calculate the spatial relationship between the first geometric object and the second geometric object, and filter candidate road segments that intersect with the provincial boundary based on the spatial relationship; Based on the candidate road segments, the road network is expanded to obtain the target road segments. Based on the target road segment, area data of inter-provincial road segments is generated in a planar format.

2. The method for automatic identification of inter-provincial road sections based on geometric space calculation according to claim 1, characterized in that, The process of obtaining spatial geometric data of road networks and provincial administrative divisions includes: Obtain the original road network trajectory coordinate data and provincial administrative division boundary coordinate data; The road network trajectory coordinate data and provincial administrative division boundary coordinate data are batch converted into unified target coordinate system data.

3. The method for automatic identification of inter-provincial road sections based on geometric space calculation according to claim 1, characterized in that, The process of constructing the first geometric object corresponding to the road segment and the second geometric object corresponding to the provincial boundary includes: The road network trajectory coordinate data after unifying the coordinate system is deserialized using WKT and constructed as a linear geometric object, which serves as the first geometric object. The coordinate data of provincial administrative division boundaries after unifying the coordinate system are deserialized using WKT and constructed into polygonal geometric objects, which serve as the second geometric object.

4. The method for automatic identification of inter-provincial road sections based on geometric space calculation according to claim 1, characterized in that, Calculating the spatial relationship between the first geometric object and the second geometric object includes: Extract the extreme values ​​of all coordinate points of the first geometric object and construct the minimum bounding rectangle; determine the spatial overlap relationship between the minimum bounding rectangle and the second geometric object; when the spatial overlap relationship exists, perform a precise spatial relationship calculation between the first geometric object and the second geometric object; When the result of the precise spatial relationship calculation is an intersection, the road segment represented by the corresponding first geometric object is identified as the candidate road segment.

5. The method for automatic identification of inter-provincial road sections based on geometric space calculation according to claim 4, characterized in that, The precise intersection relationship between the first geometric object and the second geometric object is calculated using the geometric object intersection function.

6. The method for automatic identification of inter-provincial road sections based on geometric space calculation according to claim 1, characterized in that, The process of expanding the road network based on the candidate road segments includes: Based on the topological connection relationship of the road network, the direct upstream adjacent road segment and the direct downstream adjacent road segment of the candidate road segment are obtained; the candidate road segment is connected with the direct upstream adjacent road segment and the direct downstream adjacent road segment to obtain a continuous extended road segment, which is used as the target road segment.

7. The method for automatic identification of inter-provincial road sections based on geometric space calculation according to claim 3, characterized in that, Based on the target road segment, area-based inter-provincial road segment data is generated, including: The geometric objects of the target road segment are converted into a projected coordinate system suitable for linear distance calculation; A buffer analysis operation is performed on the geometric object in the projected coordinate system to generate planar inter-provincial road segment area data; Convert the area data of inter-provincial road sections in planar form back to the target coordinate system.

8. An automatic identification system for inter-provincial road sections based on geometric space calculation, used to implement the method described in any one of claims 1-7, characterized in that, include: The data acquisition module is used to acquire spatial geometric data of road networks and provincial administrative divisions; A geometric object construction module, connected to the data acquisition module, is used to construct a first geometric object corresponding to a road segment and a second geometric object corresponding to a provincial boundary based on the spatial geometric data. A spatial relationship calculation module, connected to the geometric object construction module, is used to calculate the spatial relationship between the first geometric object and the second geometric object; The candidate road segment screening module is connected to the spatial relationship calculation module and is used to screen candidate road segments that intersect with the provincial boundary based on the spatial relationship. A road segment expansion module, connected to the candidate road segment filtering module, is used to perform expansion processing on the road network based on the candidate road segments to obtain the target road segment; The region generation module, connected to the road segment expansion module, is used to generate area-based inter-provincial road segment region data based on the target road segment.

Images