CN116311916A - Method for estimating discrete point isochronous characteristic value of land surface generated by facing traffic isochronous line - Google Patents

Method for estimating discrete point isochronous characteristic value of land surface generated by facing traffic isochronous line Download PDF

Info

Publication number
CN116311916A
CN116311916A CN202310149481.1A CN202310149481A CN116311916A CN 116311916 A CN116311916 A CN 116311916A CN 202310149481 A CN202310149481 A CN 202310149481A CN 116311916 A CN116311916 A CN 116311916A
Authority
CN
China
Prior art keywords
discrete points
discrete
isochronous
class
points
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202310149481.1A
Other languages
Chinese (zh)
Other versions
CN116311916B (en
Inventor
马莹莹
陈曦
许明朗
谢文艺
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN202310149481.1A priority Critical patent/CN116311916B/en
Publication of CN116311916A publication Critical patent/CN116311916A/en
Application granted granted Critical
Publication of CN116311916B publication Critical patent/CN116311916B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G08—SIGNALLING
    • G08G—TRAFFIC CONTROL SYSTEMS
    • G08G1/00—Traffic control systems for road vehicles
    • G08G1/01—Detecting movement of traffic to be counted or controlled
    • G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0125—Traffic data processing
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00—Road transport of goods or passengers
    • Y02T10/10—Internal combustion engine [ICE] based vehicles
    • Y02T10/40—Engine management systems

Landscapes

  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Traffic Control Systems (AREA)

Abstract

The invention discloses a land surface discrete point isochronous characteristic value estimation method for traffic isochronous line generation, which comprises the following steps: 1) Rasterizing a research area of the urban road network, wherein a grid central point is a land surface discrete point; 2) The discrete points of the land surface are divided into the following steps according to the space connection relation between the discrete points of the land surface and the discrete points of the road line and the road section: 1. 2,3 class surface discrete points, and assigning a maximum value for all land block surface discrete points as an initial value; 3) Calculating the characteristic values of discrete points of the class 1 and class 2 surfaces; 4) Calculating the isochronous characteristic values of the 3-class surface discrete points, checking whether the isochronous characteristic values of all the land surface discrete points are smaller than the initial value, and repeating the step 4) until the condition is met if the isochronous characteristic values of all the land surface discrete points are not met; 5) And deleting the 1-class surface discrete points, and outputting the isochronous characteristic values of the 2-class surface discrete points and the 3-class surface discrete points until the estimation is completed. The invention breaks through the defect of the existing traffic isochrone generation method in consideration of discrete points of the land surface, and refines the generation process of the traffic isochrone.

Description

面向交通等时线生成的地块面离散点等时特征值估计方法Estimation Method of Isochronous Eigenvalues of Discrete Points of Land Parcel for Traffic Isochrone Line Generation

技术领域technical field

本发明涉及城市道路路网时效性能评估与分析的技术领域，尤其是指一种面向交通等时线生成的地块面离散点等时特征值估计方法。The invention relates to the technical field of time-effect performance evaluation and analysis of urban road networks, in particular to a method for estimating isochronous eigenvalues of discrete points on a plot surface oriented to generation of traffic isochrones.

背景技术Background technique

路网的时效性能是路网在时间维度承载交通流能力的体现，是路网规划设计和管理的重要评价参数。交通等时线是从路网中某点出发，行程时间相等的所有点相连形成的闭合曲线，是交通网络上行程时间在空间分布状态的反映，是研究路网空间时效性能的重要手段。而路网离散点的等时特征值估计是交通等时线生成过程中的重要环节。The timeliness performance of the road network is the embodiment of the road network's ability to carry traffic flow in the time dimension, and is an important evaluation parameter for road network planning, design and management. The traffic isochrone is a closed curve formed by connecting all points with equal travel time starting from a certain point in the road network. It is a reflection of the spatial distribution of travel time on the traffic network and an important means to study the spatial time performance of the road network. The estimation of isochronous eigenvalues of discrete points in road network is an important link in the process of generating traffic isochrones.

路网离散点根据其所在位置可分为道路线离散点和地块面离散点。其中道路线离散点的等时特征值可以通过在线地图服务直接获取或是通过最短路算法等计算得到，而地块面离散点脱离道路而存在，其等时特征值无法像道路线离散点一样简单计算得到。已有的等时线生成方法中，要么采用均匀离散，没有区分道路线离散点和地块面离散点；要么仅考虑道路线离散点，没有考虑地块面离散点的等时特征值估计。前者无法生成准确的交通等时线，且耗时较长；后者可能会因为道路的离散空间分布无法得到理想的插值结果，从而生成带有“空洞”、“气泡”的交通等时线。The road network discrete points can be divided into road line discrete points and plot surface discrete points according to their location. Among them, the isochronous eigenvalues of the discrete points of the road line can be obtained directly through the online map service or calculated by the shortest path algorithm, while the discrete points of the plot surface exist away from the road, and their isochronous eigenvalues cannot be the same as the discrete points of the road line. It can be calculated simply. In the existing isochrone generation methods, either uniform discretization is adopted without distinguishing the discrete points of the road line and the discrete points of the plot surface; or only the discrete points of the road line are considered, and the isochronous characteristic value estimation of the discrete points of the plot surface is not considered. The former cannot generate accurate traffic isochrones and takes a long time; the latter may not be able to obtain ideal interpolation results due to the discrete spatial distribution of roads, thus generating traffic isochrones with "holes" and "bubbles".

本发明提出的面向交通等时线生成的地块面离散点等时特征值估计方法细化了交通等时线的生成流程，为高效且准确的生成交通等时线打好基础。本发明的方法对基于路网不均匀离散方法得到的地块面离散点进行分类，根据地块面离散点与道路线离散点的空间关系将地块面离散点分为三类，不同类别的地块面离散点的等时特征值分别根据道路线离散点等时特征值的估计结果计算得到。The method for estimating the isochronous eigenvalues of the discrete points of the land plot surface oriented to the generation of the traffic isochrone proposed by the present invention refines the generation process of the traffic isochrone, and lays a solid foundation for efficient and accurate generation of the traffic isochrone. The method of the present invention classifies the discrete points on the plot surface obtained based on the uneven discrete method of the road network, and divides the discrete points on the plot surface into three categories according to the spatial relationship between the discrete points on the plot surface and the discrete points on the road line. The isochronous eigenvalues of the discrete points on the plot surface are calculated according to the estimated results of the isochronous eigenvalues of the discrete points on the road line.

发明内容Contents of the invention

本发明的目的在于克服现有技术的缺点和不足，提供一种面向交通等时线生成的地块面离散点等时特征值估计方法，突破已有的交通等时线生成方法在对地块面离散点的考虑缺失，精细化了交通等时线的生成过程。The purpose of the present invention is to overcome the shortcomings and deficiencies of the prior art, to provide a method for estimating the isochronous eigenvalues of the discrete points of the plot surface for traffic isochrone generation, and to break through the existing traffic isochrone generation method in the land plot The lack of consideration of surface discrete points refines the generation process of traffic isochrones.

为实现上述目的，本发明所提供的技术方案为：面向交通等时线生成的地块面离散点等时特征值估计方法，包括以下步骤：In order to achieve the above object, the technical solution provided by the present invention is: a method for estimating the isochronous eigenvalues of the discrete points of the plot surface generated for traffic isochrones, comprising the following steps:

1)将城市路网的研究区域栅格化，即分成多个边长为s的方格，方格中心点即为地块面离散点Pp(i)，i＝1,2,3,…,n，n代表地块面离散点数量；1) Grid the research area of the urban road network, that is, divide it into multiple grids with side length s, and the center point of the grid is the discrete point P p (i) on the plot surface, i=1,2,3, ..., n, n represents the number of discrete points on the plot surface;

2)通过地块面离散点Pp(i)与道路线离散点PL(j)和路段ek的空间连接关系，j＝1,2,3,…,m，m代表道路线离散点数量，将地块面离散点分为：1类面离散点

Figure BDA0004090280680000021
2类面离散点/>
Figure BDA0004090280680000022
和3类面离散点/>
Figure BDA0004090280680000023
并为所有地块面离散点赋一个极大值M作为初始值；2) Through the spatial connection relationship between the discrete point P p (i) of the plot surface and the discrete point P L (j) of the road line and the road section e k , j=1,2,3,...,m, m represents the discrete point of the road line Quantity, divide the discrete points on the plot surface into: Type 1 discrete points on the surface
Figure BDA0004090280680000021
2 types of surface discrete points/>
Figure BDA0004090280680000022
and 3 types of surface discrete points/>
Figure BDA0004090280680000023
And assign a maximum value M as the initial value for all discrete points on the plot surface;

3)估计不同类别的地块面离散点的等时特征值，由已知的道路线离散点等时特征值

Figure BDA0004090280680000024
分别计算1类面离散点等时特征值/>
Figure BDA0004090280680000025
和2类面离散点等时特征值/>
Figure BDA0004090280680000026
3) Estimate the isochronous eigenvalues of the discrete points of different types of plots, from the known isochronous eigenvalues of the discrete points of the road line
Figure BDA0004090280680000024
Calculating the isochronous eigenvalues of discrete points on a class 1 surface respectively/>
Figure BDA0004090280680000025
and isochronous eigenvalues of discrete points on two types of surfaces/>
Figure BDA0004090280680000026

4)基于3类面离散点

Figure BDA0004090280680000027
近邻的其它地块面离散点，计算3类面离散点等时特征值/>
Figure BDA0004090280680000028
检查是否所有地块面离散点的等时特征值均小于初始值M，若不满足则重复步骤4)直到满足条件；4) Discrete points based on three types of surfaces
Figure BDA0004090280680000027
Calculate the isochronous eigenvalues of the three types of surface discrete points on other adjacent plot surface discrete points/>
Figure BDA0004090280680000028
Check whether the isochronous eigenvalues of all discrete points on the plot surface are less than the initial value M, if not, repeat step 4) until the condition is met;

5)删去1类面离散点，输出2类面离散点和3类面离散点的等时特征值，至此，所有地块面离散点的等时特征值均估计完毕。5) Delete the discrete points of type 1 surface, and output the isochronous eigenvalues of discrete points of type 2 and type 3 surfaces. So far, the isochronous eigenvalues of the discrete points of all plot surfaces have been estimated.

进一步，在步骤1)中，通过栅格化后得到的方格的边长s来确定离散精度，以方格中心确定均匀分布的地块面离散点。Further, in step 1), the discrete precision is determined by the side length s of the grid obtained after rasterization, and the evenly distributed discrete points on the plot surface are determined by the center of the grid.

进一步，在步骤2)中，根据地块面离散点Pp(i)所在的方格范围内，与道路线离散点PL(j)和路段ek在空间上的关联，能够将地块面离散点分为以下三类离散点，判断标准如下：Further, in step 2), according to the spatial association between the discrete point P p (i) on the plot surface and the spatial association with the discrete point P L (j) of the road line and the road section e k within the grid area, the plot can be Surface discrete points are divided into the following three types of discrete points, and the judgment criteria are as follows:

1类面离散点

Figure BDA0004090280680000031
地块面离散点Pp(i)所属方格同时与道路线离散点PL(j)和路段ek相交；Type 1 Surface Discrete Points
Figure BDA0004090280680000031
The grid of the discrete point P p (i) on the plot surface intersects the discrete point P L (j) of the road line and the road section e k at the same time;

2类面离散点

Figure BDA0004090280680000032
地块面离散点Pp(i)所属方格不与道路线离散点PL(j)相交，但与路段ek相交；2 types of surface discrete points
Figure BDA0004090280680000032
The square to which the discrete point P p (i) of the plot surface belongs does not intersect with the discrete point P L (j) of the road line, but intersects with the road segment e k ;

3类面离散点

Figure BDA0004090280680000033
地块面离散点Pp(i)所属方格与道路线离散点PL(j)和路段ek均不相交；3 types of surface discrete points
Figure BDA0004090280680000033
The grid of the discrete point P p (i) on the plot surface does not intersect with the discrete point P L (j) of the road line and the road section e k ;

根据上述的判断标准，确定每个地块面离散点的归属，同时，将一个极大值M作为初始值，赋值给所有地块面离散点，并且，该初始值用于后续步骤中地块面离散点的等时特征值的计算，作为判断一个地块面离散点是否被估计赋值的基准。According to the above judgment criteria, determine the attribution of discrete points on each plot surface, and at the same time, assign a maximum value M as an initial value to all discrete points on the plot surface, and this initial value is used for the plots in subsequent steps The calculation of isochronous eigenvalues of surface discrete points is used as a benchmark for judging whether a plot surface discrete point is estimated or not.

进一步，在步骤3)中，道路线离散点等时特征值

Figure BDA0004090280680000034
作为已知条件输入，用于估计地块面离散点等时特征值；Further, in step 3), the isochronous eigenvalues of the road line discrete points
Figure BDA0004090280680000034
As a known condition input, it is used to estimate the isochronous eigenvalues of discrete points on the plot surface;

1类面离散点等时特征值

Figure BDA0004090280680000035
等于与该1类面离散点/>
Figure BDA0004090280680000036
所属方格相交的道路线离散点等时特征值的平均值，如下所示：Isochronous Eigenvalues of Discrete Points on Type 1 Surface
Figure BDA0004090280680000035
is equal to the discrete point of the class 1 surface />
Figure BDA0004090280680000036
The average value of the isochronous eigenvalues of the road line discrete points intersected by the grid is as follows:

Figure BDA0004090280680000037
Figure BDA0004090280680000037

式中，

Figure BDA0004090280680000038
表示对应1类面离散点/>
Figure BDA0004090280680000039
所属方格内相交的所有道路线离散点的集合；/>
Figure BDA00040902806800000310
表示集合/>
Figure BDA00040902806800000311
内某一道路线离散点/>
Figure BDA00040902806800000312
的等时特征值，其中c表示集合/>
Figure BDA00040902806800000313
中的索引编号；In the formula,
Figure BDA0004090280680000038
Indicates the discrete point corresponding to a type 1 surface />
Figure BDA0004090280680000039
A collection of discrete points of all road lines intersecting in the grid;/>
Figure BDA00040902806800000310
represents a collection />
Figure BDA00040902806800000311
Discrete point within a certain road line />
Figure BDA00040902806800000312
isochronous eigenvalues of , where c represents the set />
Figure BDA00040902806800000313
the index number in the ;

2类面离散点

Figure BDA00040902806800000314
与路段相邻，所以2类面离散点等时特征值/>
Figure BDA00040902806800000315
为与其直线距离最小的道路线离散点的等时特征值加上误差时间，如下所示：2 types of surface discrete points
Figure BDA00040902806800000314
Adjacent to the road section, so the isochronous eigenvalues of the discrete points of the 2 types of surfaces />
Figure BDA00040902806800000315
Add the error time to the isochronous characteristic value of the discrete point of the road line with the smallest straight-line distance, as follows:

Figure BDA00040902806800000316
Figure BDA00040902806800000316

式中，PL(N)表示与对应2类面离散点

Figure BDA0004090280680000041
直线距离最近的道路线离散点，/>
Figure BDA0004090280680000042
表示PL(N)的等时特征值；/>
Figure BDA0004090280680000043
和vp分别表示/>
Figure BDA0004090280680000044
和PL(N)两点间的欧式距离和估计步行速度，用于估计误差时间。In the formula, P L (N) represents the discrete points corresponding to the two types of surfaces
Figure BDA0004090280680000041
Straight line distance to the nearest discrete point on the road line, />
Figure BDA0004090280680000042
Indicates the isochronous eigenvalue of PL (N); />
Figure BDA0004090280680000043
and v p represent /> respectively
Figure BDA0004090280680000044
The Euclidean distance and estimated walking speed between two points and PL (N) are used to estimate the error time.

进一步，在步骤4)中，因为3类面离散点

Figure BDA0004090280680000045
不与路段相邻，则3类面离散点等时特征值/>
Figure BDA0004090280680000046
为临近的地块面离散点的等时特征值的最小值加上误差时间，如下所示：Further, in step 4), because the three types of surface discrete points
Figure BDA0004090280680000045
is not adjacent to the road section, then the isochronous eigenvalues of the three types of surface discrete points/>
Figure BDA0004090280680000046
Add the error time to the minimum value of the isochronous eigenvalues of the discrete points of the adjacent plot surface, as follows:

Figure BDA0004090280680000047
Figure BDA0004090280680000047

Figure BDA0004090280680000048
Figure BDA0004090280680000048

式中，Pp near表示与3类面离散点

Figure BDA0004090280680000049
相邻的地块面离散点的集合，/>
Figure BDA00040902806800000410
表示集合Pp near内某一地块面离散点Pp near(b)的等时特征值，其中b表示集合Pp near中的索引编号；s和vp分别表示地块面离散点所属方格的边长和估计步行速度，Δt0表示相邻地块面离散点之间的时间误差。In the formula, P p near represents the discrete point with the three types of surfaces
Figure BDA0004090280680000049
A collection of discrete points on adjacent parcel faces, />
Figure BDA00040902806800000410
Indicates the isochronous eigenvalue of a discrete point P p near (b) on a block surface in the set P p near , where b represents the index number in the set P p near ; s and v p represent the location of the discrete point on the block surface grid side length and estimated walking speed, Δt 0 represents the time error between discrete points on adjacent plot faces.

进一步，在步骤5)中，因为1类面离散点空间位置与道路线离散点接近，所以在辅助计算3类面离散点等时特征值之后，在生成交通等时线的后续过程中删去，以避免数据冗余。Furthermore, in step 5), since the spatial position of the discrete points of type 1 surface is close to the discrete points of the road line, after the auxiliary calculation of the isochronous eigenvalues of the discrete points of type 3 surface, it is deleted in the subsequent process of generating the traffic isochrone , to avoid data redundancy.

本发明与现有技术相比，具有如下优点与有益效果：Compared with the prior art, the present invention has the following advantages and beneficial effects:

1)细化了交通等时线的生成过程，在考虑道路线离散点的基础上，还额外考虑了地块面离散点等时特征值的估计。1) The generation process of the traffic isochrone is refined. On the basis of considering the discrete points of the road line, the estimation of the isochronous eigenvalues of the discrete points of the plot surface is also considered.

2)根据空间连接关系，将地块面离散点分为了三类，根据不同的位置关系采用不同的等时特征值估计方法，提高了地块面离散点等时特征值估计的准确性。2) According to the spatial connection relationship, the discrete points on the plot surface are divided into three categories, and different isochronous eigenvalue estimation methods are used according to different positional relationships, which improves the accuracy of the isochronous eigenvalue estimation of the discrete points on the plot surface.

3)在计算2类和3类地块面离散点的等时特征值时，均考虑了不同的时间误差，即考虑了交通影响因素，使得地块面等时特征值的估计更贴合实际，计算结果也更为准确。3) When calculating the isochronous eigenvalues of the discrete points of the 2nd and 3rd types of plots, different time errors are considered, that is, traffic factors are considered, so that the estimation of the isochronous eigenvalues of the plots is more realistic , and the calculation result is more accurate.

附图说明Description of drawings

图1为本发明方法的流程图。Fig. 1 is the flowchart of the method of the present invention.

图2为实施案例路网示意图。Figure 2 is a schematic diagram of the road network of the implementation case.

图3为三种面离散点示意图。Figure 3 is a schematic diagram of discrete points on three surfaces.

图4为实施案例地块面离散点等时特征值计算结果示意图。Figure 4 is a schematic diagram of the calculation results of the isochronous eigenvalues of the discrete points on the land surface of the implementation case.

具体实施方式Detailed ways

下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。The present invention will be further described in detail below in conjunction with the embodiments and the accompanying drawings, but the embodiments of the present invention are not limited thereto.

如图1所示，本实施例公开了一种面向交通等时线生成的地块面离散点等时特征值估计方法，具体步骤如下：As shown in Figure 1, this embodiment discloses a method for estimating isochronous eigenvalues of discrete points on a plot face generated for traffic isochrones, and the specific steps are as follows:

1)将城市路网的研究区域栅格化，即分成多个边长为s的方格，方格中心点即为地块面离散点Pp(i)，i＝1,2,3,…,n，n代表地块面离散点数量；其中，通过栅格化后得到的方格的边长s来确定离散精度，以方格中心确定均匀分布的地块面离散点。1) Grid the research area of the urban road network, that is, divide it into multiple grids with side length s, and the center point of the grid is the discrete point P p (i) on the plot surface, i=1,2,3, ..., n, n represents the number of discrete points on the plot surface; among them, the discretization accuracy is determined by the side length s of the grid obtained after rasterization, and the evenly distributed discrete points on the plot surface are determined by the center of the grid.

2)通过地块面离散点Pp(i)与道路线离散点PL(j)和路段ek的空间连接关系，j＝1,2,3,…,m，m代表道路线离散点数量，将地块面离散点分为：1类面离散点

Figure BDA0004090280680000051
2类面离散点/>
Figure BDA0004090280680000052
和3类面离散点/>
Figure BDA0004090280680000053
并为所有地块面离散点赋一个极大值M作为初始值，具体如下：2) Through the spatial connection relationship between the discrete point P p (i) of the plot surface and the discrete point P L (j) of the road line and the road section e k , j=1,2,3,...,m, m represents the discrete point of the road line Quantity, divide the discrete points on the plot surface into: Type 1 discrete points on the surface
Figure BDA0004090280680000051
2 types of surface discrete points/>
Figure BDA0004090280680000052
and 3 types of surface discrete points/>
Figure BDA0004090280680000053
And assign a maximum value M as the initial value for all discrete points on the plot surface, as follows:

根据地块面离散点Pp(i)所在的方格范围内，与道路线离散点PL(j)和路段ek在空间上的关联，能够将地块面离散点分为以下三类离散点，判断标准如下：According to the spatial relationship between the discrete point P p (i) of the plot surface and the discrete point P L (j) of the road line and the road segment e k within the grid range, the discrete points of the plot surface can be divided into the following three types Discrete points, the judgment criteria are as follows:

1类面离散点

Figure BDA0004090280680000054
地块面离散点Pp(i)所属方格同时与道路线离散点PL(j)和路段ek相交；Type 1 Surface Discrete Points
Figure BDA0004090280680000054
The grid of the discrete point P p (i) on the plot surface intersects the discrete point P L (j) of the road line and the road section e k at the same time;

2类面离散点

Figure BDA0004090280680000061
地块面离散点Pp(i)所属方格不与道路线离散点PL(j)相交，但与路段ek相交；2 types of surface discrete points
Figure BDA0004090280680000061
The square to which the discrete point P p (i) of the plot surface belongs does not intersect with the discrete point P L (j) of the road line, but intersects with the road segment e k ;

3类面离散点

Figure BDA0004090280680000062
地块面离散点Pp(i)所属方格与道路线离散点PL(j)和路段ek均不相交；3 types of surface discrete points
Figure BDA0004090280680000062
The grid of the discrete point P p (i) on the plot surface does not intersect with the discrete point P L (j) of the road line and the road section e k ;

根据上述的判断标准，确定每个地块面离散点的归属，同时，将一个极大值M作为初始值，赋值给所有地块面离散点，并且，该初始值用于后续步骤中地块面离散点的等时特征值的计算，作为判断一个地块面离散点是否被估计赋值的基准。According to the above judgment criteria, determine the attribution of discrete points on each plot surface, and at the same time, assign a maximum value M as an initial value to all discrete points on the plot surface, and this initial value is used for the plots in subsequent steps The calculation of isochronous eigenvalues of surface discrete points is used as a benchmark for judging whether a plot surface discrete point is estimated or not.

3)估计不同类别的地块面离散点的等时特征值，由已知的道路线离散点等时特征值tPL(j)，分别计算1类面离散点等时特征值

Figure BDA0004090280680000063
和2类面离散点等时特征值/>
Figure BDA0004090280680000064
具体如下：3) Estimate the isochronous eigenvalues of different types of discrete points on the plot surface, and calculate the isochronous eigenvalues of the discrete points of the first type of surface from the known isochronous eigenvalues t PL(j) of the road line discrete points
Figure BDA0004090280680000063
and isochronous eigenvalues of discrete points on two types of surfaces/>
Figure BDA0004090280680000064
details as follows:

道路线离散点等时特征值

Figure BDA0004090280680000065
作为已知条件输入，用于估计地块面离散点等时特征值；Isochronous Eigenvalues of Discrete Points on Road Lines
Figure BDA0004090280680000065
As a known condition input, it is used to estimate the isochronous eigenvalues of discrete points on the plot surface;

1类面离散点等时特征值

Figure BDA0004090280680000066
等于与该1类面离散点/>
Figure BDA0004090280680000067
所属方格相交的道路线离散点等时特征值的平均值，如下所示：Isochronous Eigenvalues of Discrete Points on Type 1 Surface
Figure BDA0004090280680000066
is equal to the discrete point of the class 1 surface />
Figure BDA0004090280680000067
The average value of the isochronous eigenvalues of the road line discrete points intersected by the grid is as follows:

Figure BDA0004090280680000068
Figure BDA0004090280680000068

式中，

Figure BDA0004090280680000069
表示对应1类面离散点/>
Figure BDA00040902806800000610
所属方格内相交的所有道路线离散点的集合；/>
Figure BDA00040902806800000611
表示集合/>
Figure BDA00040902806800000612
内某一道路线离散点/>
Figure BDA00040902806800000613
的等时特征值，其中c表示集合/>
Figure BDA00040902806800000614
中的索引编号；In the formula,
Figure BDA0004090280680000069
Indicates the discrete point corresponding to a type 1 surface />
Figure BDA00040902806800000610
A collection of discrete points of all road lines intersecting in the grid;/>
Figure BDA00040902806800000611
represents a collection />
Figure BDA00040902806800000612
Discrete point within a certain road line />
Figure BDA00040902806800000613
isochronous eigenvalues of , where c represents the set />
Figure BDA00040902806800000614
the index number in the ;

2类面离散点

Figure BDA00040902806800000615
与路段相邻，所以2类面离散点等时特征值/>
Figure BDA00040902806800000616
为与其直线距离最小的道路线离散点的等时特征值加上误差时间，如下所示：2 types of surface discrete points
Figure BDA00040902806800000615
Adjacent to the road section, so the isochronous eigenvalues of the discrete points of the 2 types of surfaces />
Figure BDA00040902806800000616
Add the error time to the isochronous characteristic value of the discrete point of the road line with the smallest straight-line distance, as follows:

Figure BDA00040902806800000617
Figure BDA00040902806800000617

式中，PL(N)表示与对应2类面离散点

Figure BDA0004090280680000071
直线距离最近的道路线离散点，/>
Figure BDA0004090280680000072
表示PL(N)的等时特征值；/>
Figure BDA0004090280680000073
和vp分别表示/>
Figure BDA0004090280680000074
和PL(N)两点间的欧式距离和估计步行速度，用于估计误差时间。In the formula, P L (N) represents the discrete points corresponding to the two types of surfaces
Figure BDA0004090280680000071
Straight line distance to the nearest discrete point on the road line, />
Figure BDA0004090280680000072
Indicates the isochronous eigenvalue of PL (N); />
Figure BDA0004090280680000073
and v p represent /> respectively
Figure BDA0004090280680000074
The Euclidean distance and estimated walking speed between two points and PL (N) are used to estimate the error time.

4)基于3类面离散点

Figure BDA0004090280680000075
近邻的其它地块面离散点，计算3类面离散点等时特征值/>
Figure BDA0004090280680000076
检查是否所有地块面离散点的等时特征值均小于初始值M，若不满足则重复步骤4)直到满足条件；4) Discrete points based on three types of surfaces
Figure BDA0004090280680000075
Calculate the isochronous eigenvalues of the three types of surface discrete points on other adjacent plot surface discrete points/>
Figure BDA0004090280680000076
Check whether the isochronous eigenvalues of all discrete points on the plot surface are less than the initial value M, if not, repeat step 4) until the condition is met;

因为3类面离散点

Figure BDA0004090280680000077
不与路段相邻，则3类面离散点等时特征值/>
Figure BDA0004090280680000078
为临近的地块面离散点的等时特征值的最小值加上误差时间，如下所示：Because the three types of surface discrete points
Figure BDA0004090280680000077
is not adjacent to the road section, then the isochronous eigenvalues of the three types of surface discrete points/>
Figure BDA0004090280680000078
Add the error time to the minimum value of the isochronous eigenvalues of the discrete points of the adjacent plot surface, as follows:

Figure BDA0004090280680000079
Figure BDA0004090280680000079

Figure BDA00040902806800000710
Figure BDA00040902806800000710

式中，Pp near表示与3类面离散点

Figure BDA00040902806800000711
相邻的地块面离散点的集合，/>
Figure BDA00040902806800000712
表示集合Pp near内某一地块面离散点Pp near(b)的等时特征值，其中b表示集合Pp near中的索引编号；s和vp分别表示地块面离散点所属方格的边长和估计步行速度，Δt0表示相邻地块面离散点之间的时间误差。In the formula, P p near represents the discrete point with the three types of surfaces
Figure BDA00040902806800000711
A collection of discrete points on adjacent parcel faces, />
Figure BDA00040902806800000712
Indicates the isochronous eigenvalue of a discrete point P p near ( b) on a block surface in the set P p near , where b represents the index number in the set P p near ; s and v p represent the location of the discrete point on the block surface grid side length and estimated walking speed, Δt 0 represents the time error between discrete points on adjacent plot faces.

5)因为1类面离散点空间位置与道路线离散点接近，所以在辅助计算3类面离散点等时特征值之后，在生成交通等时线的后续过程中删去，以避免数据冗余；删去1类面离散点，输出2类面离散点和3类面离散点的等时特征值，至此，所有地块面离散点的等时特征值均估计完毕。5) Because the spatial position of the discrete points of type 1 surface is close to the discrete points of the road line, after the auxiliary calculation of the isochronous eigenvalues of the discrete points of type 3 surface, it is deleted in the subsequent process of generating the traffic isochrone to avoid data redundancy ; Delete the discrete points of type 1 surface, and output the isochronous eigenvalues of the discrete points of type 2 and type 3 surfaces. So far, the isochronous eigenvalues of all the discrete points of the plot surface have been estimated.

下面我们选取广州市天河区中心商务区路网作为案例，路网面积为8.56km2。路网重点研究区域北至天河路，南至临江大道，西至广州大道中，东至猎德大道，具体如图2所示。该路网内共有140条路段，其中快速路、主干路、次干路和支路路段占比分别为27.2％、5.7％、35.7％和31.4％。对研究区域的道路条件现状、交叉口现状和信号配时现状情况进行调研，获取并整理基础路网数据与信号配时数据。Below we take the road network of the central business district of Tianhe District, Guangzhou as a case, and the area of the road network is 8.56km 2 . The key research area of the road network extends to Tianhe Road in the north, Linjiang Avenue in the south, Guangzhou Avenue Middle in the west, and Liede Avenue in the east, as shown in Figure 2. There are 140 road sections in this road network, among which expressway, trunk road, secondary trunk road and branch road sections account for 27.2%, 5.7%, 35.7% and 31.4% respectively. Investigate the status quo of road conditions, intersections and signal timing in the study area, and obtain and organize basic road network data and signal timing data.

道路先根据一定规则离散成点。而路网目标范围作栅格化，即划分得到排列整齐的多个小方格，以方格的中心点为地块面离散点。再根据空间连接确定地块面离散点的分类，各分类情况如图3所示。根据不同的计算方法，分别估计不同类别地块面离散点的等时特征值。最后即由图2所示路网得到了图4所示的带有等时特征值的地块面离散点(图中包含了道路线离散点，去除了1类面离散点)。本发明方法生成的带有等时特征值的地块面离散点，考虑了空间位置关系和交通影响因素，可有效提高后续生成的交通等时线的精度。The road is first discretized into points according to certain rules. The target range of the road network is gridded, that is, it is divided into a number of neatly arranged small squares, and the center point of the square is the discrete point of the plot surface. Then, according to the spatial connection, the classification of discrete points on the plot surface is determined, and the classification of each classification is shown in Figure 3. According to different calculation methods, the isochronous eigenvalues of discrete points on different types of plots are estimated respectively. Finally, from the road network shown in Figure 2, the discrete points on the plot surface with isochronous eigenvalues shown in Figure 4 are obtained (the discrete points of the road line are included in the figure, and the discrete points of the first type of surface are removed). The discrete points on the plot surface with isochronous characteristic values generated by the method of the present invention take into account the spatial position relationship and traffic influence factors, and can effectively improve the accuracy of the subsequently generated traffic isochrones.

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。The above-mentioned embodiment is a preferred embodiment of the present invention, but the embodiment of the present invention is not limited by the above-mentioned embodiment, and any other changes, modifications, substitutions, combinations, Simplifications should be equivalent replacement methods, and all are included in the protection scope of the present invention.

Claims (6)

1. The method for estimating the land surface discrete point isochronal characteristic value generated by facing the traffic isochrone is characterized by comprising the following steps of:
1) The investigation region of the urban road network is rasterized,namely, dividing the land into a plurality of square grids with side length s, wherein the center point of each square grid is the discrete point P of the land surface p (i) I=1, 2,3,..n, n represents the number of discrete points of the plot;
2) Through discrete point P of land surface p (i) Discrete point P from road line L (j) And road segment e k J=1, 2,3,..m, m represents the number of road line discrete points, dividing the land surface discrete points into: class 1 surface discrete points
Figure FDA0004090280670000011
Class 2 surface discrete points
Figure FDA0004090280670000012
And class 3 surface discrete points->
Figure FDA0004090280670000013
A maximum value M is given to all land surface discrete points as an initial value;
3) Estimating isochronous characteristic values of discrete points of land areas of different categories, and using known isochronous characteristic values of discrete points of road lines
Figure FDA0004090280670000014
Respectively calculating the characteristic value +.1 class surface discrete points and the like>
Figure FDA0004090280670000015
And class 2 surface discrete point isochronal eigenvalue +.>
Figure FDA0004090280670000016
4) Based on 3 kinds of surface discrete points
Figure FDA0004090280670000017
Calculating the characteristic values of 3 kinds of discrete points and the like of other neighboring land block surface discrete points
Figure FDA0004090280670000018
Check whether or not to all ofThe isochronous characteristic values of the discrete points of the block surface are smaller than the initial value M, and if the isochronous characteristic values are not met, the step 4) is repeated until the condition is met;
5) And deleting the 1-class surface discrete points, outputting the isochronous characteristic values of the 2-class surface discrete points and the 3-class surface discrete points, and finishing the estimation of the isochronous characteristic values of all the land block surface discrete points.
2. The method for estimating the land surface discrete point isochronous characteristic value generated for the traffic isochronous line according to claim 1, wherein the method comprises the following steps: in step 1), discrete accuracy is determined by the side length s of the grid obtained after rasterization, and uniformly distributed land discrete points are determined by the grid center.
3. The method for estimating the land surface discrete point isochronous characteristic value for traffic isochronous line generation according to claim 2, wherein the method comprises the steps of: in step 2), according to the land surface discrete point P p (i) Within the square area, the points P are separated from the road line L (j) And road segment e k The space association can divide the discrete points of the land surface into the following three types of discrete points, and the judgment standard is as follows:
class 1 surface discrete points
Figure FDA0004090280670000021
Discrete point P of land surface p (i) The square grid and the discrete point P of the road line are simultaneously adopted L (j) And road segment e k Intersecting;
class 2 surface discrete points
Figure FDA0004090280670000022
Discrete point P of land surface p (i) The square does not have a discrete point P with the road line L (j) Intersecting but with road section e k Intersecting;
class 3 surface discrete points
Figure FDA0004090280670000023
Discrete point P of land surface p (i) Discrete points P of the square and road line L (j) And road segment e k Are all disjoint;
and determining the attribution of each land surface discrete point according to the judging standard, and simultaneously, assigning a maximum value M to all the land surface discrete points as an initial value, wherein the initial value is used for calculating the isochronal characteristic value of the land surface discrete point in the subsequent step and is used as a reference for judging whether the land surface discrete point is estimated and assigned.
4. The traffic isochrone-oriented land surface discrete point isochronal characteristic value estimation method of claim 3, wherein the method is characterized by: in step 3), the characteristic values of discrete points of the road line and the like
Figure FDA0004090280670000024
As known condition input, the method is used for estimating characteristic values such as discrete points of the land surface;
class 1 surface discrete point isochronous feature value
Figure FDA0004090280670000025
Equal to +.1 class of surface discrete points>
Figure FDA0004090280670000026
The average value of the characteristic values of the discrete points of the road line intersected by the square is as follows:
Figure FDA0004090280670000027
in the method, in the process of the invention,
Figure FDA0004090280670000028
representing the discrete point of the corresponding class 1 surface>
Figure FDA0004090280670000029
A collection of all road line discrete points intersected in the square;
Figure FDA00040902806700000210
representation set->
Figure FDA00040902806700000211
Discrete points of a certain road line in the interior->
Figure FDA00040902806700000212
Wherein c represents the set +.>
Figure FDA00040902806700000213
Index number in (a);
class 2 surface discrete points
Figure FDA00040902806700000214
Adjacent to road section, so 2 kinds of surface discrete points are equal in characteristic value +.>
Figure FDA00040902806700000215
Error time is added to the isochronous characteristic value of the road line discrete point with the smallest straight line distance, and the error time is as follows:
Figure FDA00040902806700000216
wherein P is L (N) represents and corresponds to the discrete points of the class 2 surface
Figure FDA00040902806700000217
Discrete points of the road line with the nearest straight line, +.>
Figure FDA0004090280670000031
Representing P L An isochronous characteristic value of (N); />
Figure FDA0004090280670000032
And v p Respectively indicate->
Figure FDA0004090280670000033
And P L (N) Euclidean distance between two points and estimated walking speed, for estimating error time.
5. The method for estimating the land surface discrete point isochronous characteristic value generated for the traffic isochronous line according to claim 4, wherein the method comprises the following steps: in step 4), because of the 3-class surface discrete points
Figure FDA0004090280670000034
If the characteristic value is not adjacent to the road section, the characteristic value of 3 kinds of surface discrete points and the like is->
Figure FDA0004090280670000035
Error time is added to the minimum value of the isochronous feature values of adjacent land surface discrete points as follows:
Figure FDA0004090280670000036
Figure FDA0004090280670000037
wherein P is p near Representation and class 3 surface discrete points
Figure FDA0004090280670000038
A set of discrete points of adjacent plots, +.>
Figure FDA0004090280670000039
Representing the set P p near Discrete point P of certain land surface p near (b) Of (2), wherein b Representing the set P p near Index number in (a); s and v p Respectively representing the side length of the square to which the discrete points of the land surface belong and the estimated walking speed, delta t 0 Representing the time error between discrete points of adjacent plots.
6. The method for estimating the land surface discrete point isochronous characteristic value based on the traffic isochronous line generation according to claim 5, wherein the method comprises the following steps: in step 5), since the 1-class surface discrete point spatial position is close to the road line discrete point, after the auxiliary calculation of the 3-class surface discrete point or the like characteristic value, the characteristic value is deleted in the subsequent process of generating the traffic isochrone so as to avoid data redundancy.
CN202310149481.1A 2023-02-21 2023-02-21 Estimation method of isochronal eigenvalues of discrete points on land surface for traffic isochrone generation Active CN116311916B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310149481.1A CN116311916B (en) 2023-02-21 2023-02-21 Estimation method of isochronal eigenvalues of discrete points on land surface for traffic isochrone generation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310149481.1A CN116311916B (en) 2023-02-21 2023-02-21 Estimation method of isochronal eigenvalues of discrete points on land surface for traffic isochrone generation

Publications (2)

Publication Number Publication Date
CN116311916A true CN116311916A (en) 2023-06-23
CN116311916B CN116311916B (en) 2024-06-07

Family

ID=86800578

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310149481.1A Active CN116311916B (en) 2023-02-21 2023-02-21 Estimation method of isochronal eigenvalues of discrete points on land surface for traffic isochrone generation

Country Status (1)

Country Link
CN (1) CN116311916B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6266703B1 (en) * 1992-12-29 2001-07-24 International Business Machines Corporation Method and apparatus for providing confirmation notification for isochronous data
JP2006038907A (en) * 2004-07-22 2006-02-09 Nissan Motor Co Ltd Navigation system, and map display method in navigation system
CN108022005A (en) * 2017-11-24 2018-05-11 武汉大学 A kind of efficient approximate shceme optimal transmission method
US20220164913A1 (en) * 2020-11-24 2022-05-26 Verizon Connect Ireland Limited Systems and methods for utilizing models to determine real time estimated times of arrival for scheduled appointments
CN114898013A (en) * 2022-07-15 2022-08-12 深圳市城市交通规划设计研究中心股份有限公司 Traffic isochronal generation method, electronic device and storage medium

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6266703B1 (en) * 1992-12-29 2001-07-24 International Business Machines Corporation Method and apparatus for providing confirmation notification for isochronous data
JP2006038907A (en) * 2004-07-22 2006-02-09 Nissan Motor Co Ltd Navigation system, and map display method in navigation system
CN108022005A (en) * 2017-11-24 2018-05-11 武汉大学 A kind of efficient approximate shceme optimal transmission method
US20220164913A1 (en) * 2020-11-24 2022-05-26 Verizon Connect Ireland Limited Systems and methods for utilizing models to determine real time estimated times of arrival for scheduled appointments
CN114898013A (en) * 2022-07-15 2022-08-12 深圳市城市交通规划设计研究中心股份有限公司 Traffic isochronal generation method, electronic device and storage medium

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
CHRYSOVALANTIS ANASTASIOU 等: "Time-Dependent Reachability Analysis: A Data-Driven Approach", 《2019 20TH IEEE INTERNATIONAL CONFERENCE ON MOBILE DATA MANAGEMENT (MDM)》, 8 August 2019 (2019-08-08) *
胡继华;程智锋;钟广鹏;黄泽;: "一种公交等时线的计算方法及其应用", 交通运输系统工程与信息, no. 03, 15 June 2013 (2013-06-15) *
陈曦: "路网单中心交通等时线模型生成算法及特征参数研究", 《万方数据知识服务平台》, 15 July 2023 (2023-07-15) *

Also Published As

Publication number Publication date
CN116311916B (en) 2024-06-07

Similar Documents

Publication Publication Date Title
CN108399748B (en) A Road Travel Time Prediction Method Based on Random Forest and Clustering Algorithm
CN106023589B (en) A kind of track of vehicle reconstructing method based on bayonet socket data
CN108109382A (en) A kind of congestion points based on composite network, congestion line, the discovery method of congestion regions
CN101986102B (en) Method for matching electronic map in urban geographic information system
CN108152789B (en) Utilize the passive track-corelation data correlation and localization method of RSS information
CN107563566A (en) A kind of run time interval prediction method between bus station based on SVMs
CN107248007A (en) The rural residential area intelligent optimization collocation method that a kind of double-deck domain is oriented to
CN107730113A (en) A kind of quantitative evaluation method of the urban road network planning based on function
CN116611586B (en) Newly built road network flow prediction method and system based on double-layer heterogeneous network
CN107704531A (en) The automatic merging method of road network intersection node in a kind of OSM map datums
CN107679653B (en) An OD Allocation Method Based on Predominant Travel Distance
CN110120154B (en) Traffic road condition prediction method using detector data under large-scale road network
CN115662124B (en) GPS track data road section flow matching method based on network coding
Kianfar et al. Optimizing freeway traffic sensor locations by clustering global-positioning-system-derived speed patterns
CN111462492A (en) A key road segment detection method based on Ritchie flow
CN108388540A (en) Road network Algorithms of Selecting based on fuzzy overall evaluation
CN116311916B (en) Estimation method of isochronal eigenvalues of discrete points on land surface for traffic isochrone generation
CN110580377A (en) An optimization method for design parameters of grass ditch based on response surface method
CN111008730B (en) Method and device for constructing crowd aggregation prediction model based on urban spatial structure
CN103646537A (en) Public traffic route adjustment influence set generation technique
CN103996284A (en) Urban traffic network optimization method and system
CN107121143A (en) A kind of road choosing method of collaboration POI data
CN113284030B (en) Urban traffic network community division method
CN117037492A (en) Urban traffic model construction method integrating mobile phone signaling and travel investigation data
CN114169580B (en) Traffic equal-time-circle calculation method for regional hub

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant