Space partitioning method of road network under influence of urban facility complexity

Abstract

The invention discloses a spatial partitioning method of a road network under the influence of urban facility complexity, which comprises the following steps of: firstly, dividing road network data in a city into different road sections, and generating corresponding planar buffer areas for the different road sections; secondly, spatial superposition is carried out on POI data in a road section buffer area and a city, and POI quantity characteristics in different road section buffer area ranges are calculated; thirdly, spatial superposition is carried out on the height data of the road section buffer areas and the urban internal buildings, and the height characteristics of the buildings in different road section buffer areas are calculated; fourthly, based on POI quantity features and building height features corresponding to different road sections, calculating urban facility complexity indexes in the range of the road section buffer area; and 5, generating different road network partitions according to the sequence of the facility complexity indexes from large to small and the principle that the buffer space is not superposed. The invention provides a space partitioning method of a traffic road network based on POI data and building height data for a space distribution rule of a road network under the influence of different facilities in a city. The demonstration shows that the method can represent the diversity of commercial activities and building distribution near different road sections by using the complexity index of urban facilities, and further generates the space subareas of the traffic road network based on the diversity of the commercial activities and building distribution. The invention provides a more intuitive and accurate traffic road network space partitioning method for the aspects of urban traffic management, urban planning, infrastructure operation and maintenance, urban sustainable development and the like.

Description

Technical Field

[0001] The present invention relates to the fields of geographic information technology, urban traffic management, and in particular to a spatial partitioning method for a road network under the influence of the complexity of urban facilities. Background Art

[0002] The urban traffic network is the most important and active part of the urban system operation. The effective management of the urban traffic network is of great significance for facilitating the lives of urban residents and promoting urban development. The spatial distribution of urban internal infrastructure is highly complex. Different urban facilities have different demands and impacts on traffic. Therefore, spatial zoning of the urban road network based on the differences in urban facilities where the roads are located will help achieve targeted and refined management of urban traffic.

[0003] With the widespread use of electronic maps, the spatial distribution of facilities within cities can be obtained in batches, forming a type of geographic big data represented by POI (Point of Internet) data. In addition, lidar technology, oblique photogrammetry technology, and remote sensing technology can realize the calculation of building heights within cities. The above two types of technologies provide data support for characterizing the complexity of urban facilities. At the same time, the popularization of open source geographic information technology has made the road network distribution data within the city publicly accessible, and the more commonly used is OSM (Open Street Map) data. Based on this, combining POI data, building height data, and OSM road data within the city is expected to construct a spatial partitioning method for urban transportation road networks that takes into account the complexity of urban facilities. Summary of the invention

[0004] The present invention aims at efficient management of spatial zoning of urban traffic, and provides a method for spatial zoning of road networks under the influence of urban facility complexity. The technical problems solved include the measurement of urban facility complexity and the spatial zoning of traffic road networks based on urban facility complexity.

[0005] To achieve the above object, the present invention adopts the following technical solutions:

[0006] A spatial partitioning method for road networks under the influence of urban facility complexity includes five parts: road segment unit division and buffer zone generation, calculation of POI quantity characteristics corresponding to road segments, calculation of building height characteristics corresponding to road segments, calculation of urban facility complexity index corresponding to road segments, and road network spatial partitioning based on facility complexity order.

[0007] The steps of road segment unit division and buffer zone generation are as follows:

[0008] S1: Based on the OSM road network, the road network is divided into different sections;

[0009] S2: For the vector data of the divided road segments, add fields to implement the numbering of different road segments, and then add fields to calculate the geometric length of the road segments;

[0010] S3: Generate a spatial buffer for each road section, the distance of the buffer is set to d meters, and the buffer data contains the corresponding road section number;

[0011] The calculation steps of the POI quantity characteristics corresponding to the road section are as follows:

[0012] S4: spatially intersect the planar buffer data of all road sections with the point POI data, extract the POI data corresponding to different buffers according to the number field, and calculate the number of POIs in the buffer;

[0013] S5: Calculate the quantity characteristics of the POIs in the buffer zone according to the length information of different road sections and the quantity information of the POI points in the corresponding buffer zone;

[0014] Calculation of building height features corresponding to the road section. The steps are as follows:

[0015] S6: superimpose the planar buffers of all road sections with the building height data in grid form, and calculate the average building heights within different buffer ranges;

[0016] S7: Calculate building height characteristics of the buffer zone according to the length information of different road sections and the average height information of the buildings in the corresponding buffer zone;

[0017] The calculation steps of the urban facility complexity index corresponding to the road section are as follows:

[0018] S8: Based on the POI quantity characteristics and building height characteristics corresponding to different road sections, the two are added and coupled to calculate the urban facility complexity index;

[0019] The steps of road network spatial partitioning based on the order of facility complexity are as follows:

[0020] S9: Summarize the urban facility complexity corresponding to all road segments and rearrange them in order from largest to smallest;

[0021] S10: Based on the principle that the buffer spaces of the road sections in descending order do not overlap, a partition number is assigned to each road section and the road sections around it that are less than d meters away, and finally a spatial partition is generated. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only schematic diagrams of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying creative work.

[0023] Figure 1 Invention technology flow chart;

[0024] Figure 2 A road section division result diagram provided by the present invention;

[0025] Figure 3 The POI quantity feature histogram provided by the present invention;

[0026] Figure 4 A characteristic histogram of building heights provided by the present invention;

[0027] Figure 5 A distribution map of the complexity of urban facilities provided by the present invention;

[0028] Figure 6 The urban road network spatial zoning map provided by the present invention; DETAILED DESCRIPTION

[0029] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0030] The following is the road segment unit division and buffer zone generation part:

[0031] Step 1: Collect and pre-process the open source OSM road network data in the target area, and divide the road network into different sections using location points such as road forks and traffic lights.

[0032] Step 2: Add fields to the road segment data in vector form and number them uniformly without duplication.

[0033] Step 3: Generate a buffer for each road segment, with the distance set to d meters. The number field of the buffer is consistent with the number field of the road segment.

[0034] The following is the calculation part of the POI quantity characteristics corresponding to the road section:

[0035] Step 1: Spatially intersect the road segment buffer data of the vector surface element with the vector point element POI data so that each POI point has a corresponding buffer number.

[0036] Step 2: Based on the intersection result attribute table, count the total number of POI points corresponding to each buffer zone number.

[0037] Step 3: Based on the length of each road segment and the total number of POI points in the road segment buffer, calculate the POI quantity feature corresponding to the road segment. The calculation method is as follows:

[0038]

[0039] in, is the number of POIs corresponding to road section r, N r is the number of POI points within the buffer zone of road segment r, L r is the length of road segment r.

[0040] During the calculation process, the number of POIs corresponding to a certain road section may be 0. In order not to affect the calculation of the POI quantity feature, the number of POIs is set to 0.1 in this case, which can ensure that the number of POIs is still at the minimum and the calculation process can proceed normally.

[0041] The following is the calculation part of the building height characteristics corresponding to the road section:

[0042] Step 1: Spatially overlay and clip the road section buffer data of the vector surface element and the building height data in the form of raster data so that each buffer area has corresponding building height data.

[0043] Step 2: Based on the building height data in each road section buffer, calculate the average building height in each buffer.

[0044] Step 3: Based on the length of each road segment and the average building height within the road segment buffer, calculate the building height feature corresponding to the road segment. The calculation method is as follows:

[0045]

[0046] in, is the building height feature corresponding to road section r, MH r is the average building height within the buffer zone of road section r, L r is the length of road segment r.

[0047] During the calculation process, the average building height corresponding to a certain road section may be 0. In order not to affect the calculation of the building height characteristics, the average building height is set to 0.1 in this case, that is, the average building height is still at the minimum at this time, and the calculation process can proceed normally.

[0048] The following is the calculation part of the urban facility complexity index corresponding to the road section:

[0049] Step 1: Based on the POI quantity characteristics and building height characteristics corresponding to the road section, calculate the urban facility complexity index corresponding to the road section. The calculation method is as follows:

[0050]

[0051] Among them, USCI r is the urban scene corresponding to road section r, is the POI quantity feature corresponding to road section r, is the building height feature corresponding to road section r.

[0052] The following is a section of the road network spatial partitioning based on the order of facility complexity:

[0053] Step 1: Summarize the complexity of urban facilities corresponding to all road sections and rearrange them in order from large to small.

[0054] Step 2: Set z = 1, take out the first ranked road section, that is, the road section with the largest urban scene complexity index, and set the partition label z for this road section and all road sections within the buffer range of d meters near this road section.

[0055] Step 3: Set z=z+1, take out the next road segment in the facility complexity ranking. If the taken road segment already has a partition label, continue to take the next road segment in the ranking. If the taken road segment does not have a partition label, further filter all road segments within the buffer range of d meters near this road segment. If the filtered adjacent road segments have partition labels, delete these road segments with partition labels. Set the partition labels of the current road segment and the adjacent road segments after deleting the partition labels to z.

[0056] Step 4: Repeat the previous step until all the sorted road segments are taken in sequence and all road segments have partition labels.

[0057] Step 5: Generate the spatial partition of the urban road network based on the partition labels of different road sections.

[0058] Finally, a part of the core urban area of ​​Xi'an was selected as the case area. The traffic network was derived from Open StreetMap, and the spatial resolution of the building height data was 10 meters × 10 meters. The buffer distance d = 200 meters was set, and the method of this patent was used to perform spatial partitioning of the traffic network in the case area. The results of the road network segment division in the area are shown in the figure below. Figure 2 As shown in the figure, the histogram of the calculation results of the POI quantity characteristics of each road section is as follows Figure 3 As shown in the figure, the histogram of the calculation results of the building height characteristics of each road section is as follows Figure 4 As shown in Figure 2, the spatial distribution of the urban facility complexity index results for each road section is as follows: Figure 5 As shown in the figure, the spatial distribution of the final road network partition results is as follows Figure 6 shown.

[0059] The above description of the disclosed embodiments enables those skilled in the art to implement or use the present invention. Various modifications to these embodiments will be apparent to those skilled in the art, and the general principles defined by the scheme may be implemented in other embodiments without departing from the spirit or scope of the present invention. Therefore, the present invention will not be limited to the embodiments shown in the scheme, but will conform to the widest scope consistent with the principles and novel features disclosed in the present invention.

Claims

1. A method for spatial partitioning of road networks under the influence of urban facility complexity, including road section unit division and buffer zone generation, calculation of POI quantity characteristics corresponding to road sections, calculation of building height characteristics corresponding to road sections, calculation of urban facility complexity index corresponding to road sections, and road network spatial partitioning based on facility complexity order. The steps of road segment unit division and buffer zone generation are as follows: S1: Based on the OSM road network, the road network is divided into different sections; S2: For the vector data of the divided road segments, add fields to implement different road segment numbers, and then add fields to calculate the geometric length of the road segment; S3: Generate a spatial buffer for each road section, the distance of the buffer is set to d meters, and the buffer data contains the corresponding road section number; The calculation steps of the POI quantity characteristics corresponding to the road section are as follows: S4: spatially intersecting the planar buffer data of all road sections with the point POI data, and extracting the POI data corresponding to different buffers according to the number field; S5: Calculate the quantity characteristics of the POIs in the buffer zone according to the length information of different road sections and the number of POIs in the corresponding buffer zone; Calculation of building height features corresponding to the road section. The steps are as follows: S6: superimpose the planar buffers of all road sections with the building height data in grid form, and calculate the average building heights within different buffer ranges; S7: Calculate building height characteristics of the buffer zone according to the length information of different road sections and the average height information of the buildings in the corresponding buffer zone; The calculation steps of the urban facility complexity index corresponding to the road section are as follows: S8: Based on the POI quantity characteristics and building height characteristics corresponding to different road sections, the two are coupled to calculate the urban facility complexity index; The steps of road network spatial partitioning based on the order of facility complexity are as follows: S9: Summarize the urban facility complexity corresponding to all road segments and rearrange them in order from largest to smallest; S10: Based on the principle that the buffer spaces of the road sections in descending order do not overlap, a partition number is assigned to each road section and the road sections around it that are less than d meters away, and finally a spatial partition is generated.

2. The urban facility complexity calculation according to claim 1 is characterized in that, in steps S2 to S8, the complexity of urban POI facilities and building facilities in road section distribution is considered in turn, and then the two complexities are integrated to realize the measurement calculation of urban facility complexity.

3. The road network space partitioning generation method according to claim 1 has the advantage that, in steps S9 to S10, the corresponding road network partition is determined for each road section in turn according to the order of urban facility complexity from large to small and based on the idea that the buffer space of the road section does not overlap.

Images