
技术领域technical field
本发明涉及道路交通安全技术领域,具体而言涉及一种基于组合地理自回归匹配的道路安全分析方法。The invention relates to the technical field of road traffic safety, in particular to a road safety analysis method based on combined geographic autoregressive matching.
背景技术Background technique
随着社会经济发展,社会汽车保有量持续上升,与此同时,道路安全等问题也随之产生。为了提高道路安全,交通安全管理政策在许多国家与地区开始实施。与此同时,交通安全管理政策实施后的效果也成为了相关领域的研究热点。其中倾向得分方法成为了研究道路安全管理政策与道路安全因果关系的有效技术手段,该方法的有效性得到了很多研究的认证。但是道路交通系统是一个空间网络结构,相邻交通小区之间是相互制约相互影响的,即在没有任何外界条件的影响下,相邻小区的属性也会导致临近小区产生相应的交通事故,该现象称为地理空间相关性。然而倾向得分方法并没有考虑到该问题的限制在进行因果分析中,从而导致分析结果产生一定的偏差,影响政策的有效性判断。With the development of society and economy, the number of cars in the society continues to rise, and at the same time, problems such as road safety also arise. In order to improve road safety, traffic safety management policies have been implemented in many countries and regions. At the same time, the effect of traffic safety management policy implementation has also become a research hotspot in related fields. Among them, the propensity score method has become an effective technical means to study the causal relationship between road safety management policies and road safety, and the effectiveness of this method has been verified by many studies. However, the road traffic system is a spatial network structure. Adjacent traffic cells interact with each other, that is, without the influence of any external conditions, the attributes of adjacent cells will also lead to corresponding traffic accidents in adjacent cells. The phenomenon is called geospatial correlation. However, the propensity score method does not take into account the limitations of this problem in the causal analysis, which leads to certain deviations in the analysis results and affects the judgment of the effectiveness of the policy.
发明内容SUMMARY OF THE INVENTION
本发明的目的在于提供一种基于组合地理自回归匹配的道路安全分析方法,以解决现有技术中的问题。The purpose of the present invention is to provide a road safety analysis method based on combined geographic autoregressive matching to solve the problems in the prior art.
为实现上述目的,本发明提供如下技术方案:To achieve the above object, the present invention provides the following technical solutions:
一种基于组合地理自回归匹配的道路安全分析方法,基于预设时间周期内所采集的目标区域内各个目标道路的道路信息,对目标区域内各个待优化目标道路进行道路安全分析,获得各个待优化目标道路的安全优化条件,利用安全优化条件,对待优化目标道路进行安全优化,包括以下步骤:A road safety analysis method based on combined geographic autoregressive matching, based on the road information of each target road in the target area collected within a preset time period, carries out road safety analysis on each target road to be optimized in the target area, and obtains each target road to be optimized. Optimize the safety optimization conditions of the target road, and use the safety optimization conditions to perform safety optimization of the target road to be optimized, including the following steps:
步骤A、在预设时间周期内获取目标区域中各个目标道路的道路信息、以及各个目标道路交通事故发生数量,随后进入步骤B;Step A, obtaining the road information of each target road in the target area and the number of traffic accidents on each target road within a preset time period, and then entering step B;
步骤B、分别针对各个目标道路,基于其所对应的道路信息,提取该目标道路对应的包含交通事故影响因素特征的预设各类型特征的特征值,基于所提取的交通事故影响因素特征对该目标道路进行预处理,预测得到该目标道路对应当前时间周期向未来时间方向的预设时间周期内的交通事故预测数量,随后进入步骤C;Step B: For each target road, based on the corresponding road information, extract the characteristic values of the preset various types of features corresponding to the target road including the characteristics of the influencing factors of traffic accidents, and based on the extracted characteristics of the influencing factors of traffic accidents. The target road is preprocessed, and the predicted number of traffic accidents in the preset time period corresponding to the current time period to the future time direction of the target road is predicted to be obtained, and then step C is entered;
步骤C、针对各个目标道路,基于各个目标道路的道路信息,计算并得到各个目标道路的安全系数,基于各目标道路分别所对应的安全系数,将各目标道路划分为优化目标道路和待优化目标道路,对该目标区域内的各个优化目标道路和各个待优化目标道路依次进行安全系数匹配,获得安全匹配组,即获得该目标区域对应的各个安全匹配组,随后进入步骤D;Step C. For each target road, based on the road information of each target road, calculate and obtain the safety factor of each target road, and divide each target road into optimization target roads and targets to be optimized based on the corresponding safety factors of each target road. road, perform safety factor matching on each optimization target road in the target area and each target road to be optimized in turn to obtain a safety matching group, that is, obtain each safety matching group corresponding to the target area, and then enter step D;
步骤D、分别针对各个安全匹配组,计算该安全匹配组所对应的待优化目标道路与优化目标道路之间的安全差距,进一步根据安全系数,获取影响各目标道路安全系数的影响因素,基于该安全匹配组内对应优化目标道路的道路信息、以及该优化目标道路所对应安全系数的影响因素,对待优化目标道路进行安全优化。Step D: For each safety matching group, calculate the safety gap between the target road to be optimized and the optimized target road corresponding to the safety matching group, and further obtain the influencing factors affecting the safety coefficient of each target road according to the safety factor, based on the safety factor. The road information corresponding to the optimization target road in the safety matching group and the influencing factors of the safety factor corresponding to the optimization target road are used for safety optimization of the target road to be optimized.
进一步地,前述的道路信息包括目标道路的人流密度D,目标道路的经济水平GDP,目标道路的日交通平均交通量Q、目标道路的道路网密度L,目标道路的公交站点密度B,目标道路的轨道站点密度R,目标道路的交通节点密度S。Further, the aforementioned road information includes the population density D of the target road, the economic level GDP of the target road, the average daily traffic volume Q of the target road, the road network density L of the target road, the bus stop density B of the target road, and the target road. The density of track stations R, the density of traffic nodes S of the target road.
进一步地,前述的步骤B中,对目标道路进行预处理,根据以下公式:Further, in the aforementioned step B, the target road is preprocessed according to the following formula:
计算并获得目标道路i的交通事故预测数量,xi为目标道路i对应的交通事故影响因素,βi为该目标道路i的交通事故预测数量对应的回归系数,εi为该目标道路i的交通事故预测数量对应的误差项,表示回归系数βi服从正态分布,基于该目标道路i的交通事故发生数量对所获交通事故预测数量进行修正,获得交通事故修正数量N′i,N′i=N″i-Ni,其中N″i为交通事故发生数量。Calculate and obtain the predicted number of traffic accidents on the target road i,xi is the traffic accident influencing factor corresponding to the target road i, βi is the regression coefficient corresponding to the predicted number of traffic accidents on the target road i, and εi is the target road i. The error term corresponding to the predicted number of traffic accidents, Represents that the regression coefficient βi obeys a normal distribution, and the predicted number of traffic accidents is corrected based on the number of traffic accidents on the target road i to obtain the corrected number of traffic accidents N′i , N′i =N″i -Ni , where N″i is the number of traffic accidents.
进一步地,前述的步骤C中,根据以下公式:Further, in the aforementioned step C, according to the following formula:
计算目标区域内目标道路i的安全系数pn,其中,为常数系数,μ1-μ7为各类型道路信息分别所对应的回归系数。Calculate the safety factor pn of the target road i in the target area, where, are constant coefficients, and μ1 -μ7 are regression coefficients corresponding to each type of road information.
进一步地,前述的步骤C中,基于目标区域内各个目标道路分别所对应的安全系数,对个优化目标道路和各待优化目标道路依次进行安全系数匹配,即当优化目标道路的安全系数与待优化目标道路的安全系数差最小时,该优化目标道路与该待优化目标道路为安全匹配组。Further, in the aforementioned step C, based on the safety factors corresponding to each target road in the target area, the safety factor matching is performed on each optimization target road and each target road to be optimized in turn, that is, when the safety factor of the optimized target road is the same as that of the target road to be optimized. When the safety factor difference of the optimization target road is the smallest, the optimization target road and the target road to be optimized are a safety matching group.
进一步地,前述的步骤D中,根据以下公式:Further, in the aforementioned step D, according to the following formula:
计算并获得目标区域内安全匹配组I中优化目标道路与优化目标道路之间的安全差距ATT,其中,n为该目标区域内所含目标道路的数量,N′iS为该安全匹配组中优化目标道路对应的交通事故修正系数,N′iK为该安全匹配组中待优化目标道路的交通事故修正系数。Calculate and obtain the safety gap ATT between the optimized target road and the optimized target road in the safety matching group I in the target area, where n is the number of target roads contained in the target area, and N′iS is the optimized target road in the safety matching group. The traffic accident correction coefficient corresponding to the target road, N′iK is the traffic accident correction coefficient of the target road to be optimized in the safety matching group.
本发明所述一种基于组合地理自回归匹配的道路安全分析方法,采用以上技术方案与现有技术相比,具有以下技术效果:根据本发明技术方案通过对交通事故数量进行预测,并对预测结果进行修正,使用地理自回归模型消除地理空间影响,能够较为准确的判断交通安全对道路安全事故的影响因素与影响效果,根据所获修正结果,对目标道路进行安全优化,消除因为目标区域内不同道路之前因为自身地理因素互相影响道路安全的问题,有效提升目标道路的安全性。A road safety analysis method based on combined geographic autoregressive matching according to the present invention, compared with the prior art by using the above technical solution, has the following technical effects: according to the technical solution of the present invention, by predicting the number of traffic accidents, and predicting The results are corrected, and the geographic autoregressive model is used to eliminate the influence of geographic space, which can more accurately judge the influencing factors and effects of traffic safety on road safety accidents. Previously, different roads affected each other's road safety due to their own geographical factors, which effectively improved the safety of the target road.
附图说明Description of drawings
图1为本发明示例性实施例的一种道路安全分析方法的流程示意图。FIG. 1 is a schematic flowchart of a road safety analysis method according to an exemplary embodiment of the present invention.
具体实施方式Detailed ways
为了更了解本发明的技术内容,特举具体实施例并配合所附图式说明如下。In order to better understand the technical content of the present invention, specific embodiments are given and described below in conjunction with the accompanying drawings.
在本公开中参照附图来描述本发明的各方面,附图中示出了许多说明的实施例。本公开的实施例不必定义在包括本发明的所有方面。应当理解,上面介绍的多种构思和实施例,以及下面更加详细地描述的那些构思和实施方式可以以很多方式中任意一种来实施,这是因为本发明所公开的构思和实施例并不限于任何实施方式。另外,本发明公开的一些方面可以单独使用,或者与本发明公开的其他方面的任何适当组合来使用。Aspects of the invention are described in this disclosure with reference to the accompanying drawings, in which a number of illustrative embodiments are shown. Embodiments of the present disclosure are not necessarily defined to include all aspects of the invention. It should be understood that the various concepts and embodiments described above, as well as those described in greater detail below, can be implemented in any of a number of ways, as the concepts and embodiments disclosed herein do not limited to any implementation. Additionally, some aspects of the present disclosure may be used alone or in any suitable combination with other aspects of the present disclosure.
结合图1,本发明示例性提出一种基于组合地理自回归匹配的道路安全分析方法,基于预设时间周期内所采集的目标区域上对应各个目标道路的道路信息,道路信息包括目标道路的人流密度D,目标道路的经济水平GDP,目标道路的日交通平均交通量Q、目标道路的道路网密度L,目标道路的公交站点密度B,目标道路的轨道站点密度R,目标道路的交通节点密度S,对各个目标道路进行道路安全分析,获得各个待优化目标道路的安全优化条件,对待优化目标道路进行安全优化,包括以下步骤:1, the present invention exemplarily proposes a road safety analysis method based on combined geographic autoregressive matching, based on the road information corresponding to each target road in the target area collected within a preset time period, the road information includes the flow of people on the target road. Density D, economic level GDP of the target road, average daily traffic volume Q of the target road, road network density L of the target road, bus station density B of the target road, rail station density R of the target road, traffic node density of the target road S, carry out road safety analysis on each target road, obtain the safety optimization conditions of each target road to be optimized, and perform safety optimization on the target road to be optimized, including the following steps:
步骤A、获取目标区域在预设时间周期内对应各个目标道路的道路信息、以及各个目标道路交通事故发生数量,随后进入步骤B。Step A: Obtain road information of each target road corresponding to each target road in the target area within a preset time period, and the number of traffic accidents on each target road, and then proceed to Step B.
步骤B、分别针对各个目标道路,基于目标道路所对应的道路信息,提取该目标道路对应的包含交通事故影响因素特征的预设各类型特征的特征值,基于交通事故影响因素对目标道路进行预处理,根据以下公式:Step B: For each target road, based on the road information corresponding to the target road, extract the characteristic values of the preset various types of characteristics corresponding to the target road including the characteristics of the traffic accident influencing factors, and predict the target road based on the traffic accident influencing factors. processing, according to the following formula:
计算并获得目标道路i对应当前时间点向未来时间方向的预设时间周期内的交通事故预测数量,xi为目标道路i对应的交通事故影响因素,βi为该目标道路i的交通事故预测数量对应的回归系数,εi为该目标道路i的交通事故预测数量对应的误差项,为地理回归系数遵循正态分布计算公式如下,Calculate and obtain the predicted number of traffic accidents in the preset time period from the current time point to the future time direction corresponding to the target road i, xi is the traffic accident influencing factor corresponding to the target road i, βi is the traffic accident prediction of the target road i The regression coefficient corresponding to the number, εi is the error term corresponding to the predicted number of traffic accidents on the target road i, for geographic regression coefficients to follow a normal distribution Calculated as follows,
其中wki取值为1时为该目标区域内目标道路k与目标道路i相邻,反之wki的取值为0,基于该目标道路i的交通事故发生数量对所获交通事故预测数量进行修正,获得交通事故修正数量N′i,N′i=N″i-Ni,其中N″i为交通事故发生数量,随后进入步骤C。Among them, when wki is 1, the target road k is adjacent to the target road i in the target area. Otherwise, the value of wki is 0. Based on the number of traffic accidents on the target road i, the predicted number of traffic accidents obtained is calculated. Correction, obtain the traffic accident correction number N′i , N′i =N″i -Ni , where N″i is the number of traffic accidents, and then proceed to step C.
步骤C、基于各个目标道路的道路信息,根据以下公式:Step C, based on the road information of each target road, according to the following formula:
计算目标区域内目标道路i的安全系数pn,其中,为常数系数,μ1-μ7为各类型道路信息分别所对应的回归系数,并将该目标区域内的安全系数大于预设阈值的目标道路划分为优化目标道路,将安全系数小于预设阈值的目标道路划分为待优化目标道路,Calculate the safety factor pn of the target road i in the target area, where, is a constant coefficient, μ1 - μ7 are the regression coefficients corresponding to each type of road information respectively, and the target road in the target area whose safety factor is greater than the preset threshold is divided into optimized target roads, and the safety factor is smaller than the preset threshold. The target road is divided into the target road to be optimized,
根据优化目标道路和待优化目标道路分别所对应的安全系数,提取影响目标道路安全系数的影响因素,对该目标区域内的优化目标道路和待优化目标道路进行安全系数匹配,当优化目标道路的安全系数与待优化目标道路的安全系数差最小时,该优化目标道路与该待优化目标道路为安全匹配组,即获得该目标区域对应的各个安全匹配组,随后进入步骤D。According to the corresponding safety factors of the optimized target road and the target road to be optimized, the influencing factors affecting the safety factor of the target road are extracted, and the safety factors of the optimized target road and the target road to be optimized in the target area are matched. When the difference between the safety factor and the safety factor of the target road to be optimized is the smallest, the optimization target road and the target road to be optimized are a safety matching group, that is, each safety matching group corresponding to the target area is obtained, and then step D is entered.
步骤D、分别针对各个安全匹配组,据以下公式:Step D, respectively for each security matching group, according to the following formula:
计算并获得目标区域内安全匹配组I中优化目标道路与优化目标道路之间的安全差距ATT,其中,n为该目标区域内所含目标道路的数量,N′iS为该安全匹配组中优化目标道路对应的交通事故修正系数,N′iK为该安全匹配组中待优化目标道路的交通事故修正系数,基于该安全匹配组内对应优化目标道路的道路信息、以及安全系数的影响因素,对待优化目标道路进行安全优化。Calculate and obtain the safety gap ATT between the optimized target road and the optimized target road in the safety matching group I in the target area, where n is the number of target roads contained in the target area, and N′iS is the optimized target road in the safety matching group. The traffic accident correction coefficient corresponding to the target road, N′iK is the traffic accident correction coefficient of the target road to be optimized in the safety matching group, based on the road information corresponding to the optimized target road in the safety matching group and the influencing factors of the safety coefficient, Optimize the target road for safety optimization.
实施例:Example:
采集目标区域内各个目标道路的道路信息,如表1-1所示:Collect road information of each target road in the target area, as shown in Table 1-1:
表1-1样本数据采集统计表Table 1-1 Sample data collection statistics table
基于步骤B中过程,通过以下公式:Based on the process in step B, by the following formula:
N′A1=N″A1-NA1N′A1 =N″A1 -NA1
依次获得交通事故预测数量、交通事故修正数量消除地理空间影响条件下的事故数量,该地理空间所带来的影响如区域1与区域2相邻互为匹配,在没有安全优化下,区域2的存在会使区域1本身就产生一定数量的事故,假设A1与C1互为安全匹配组,则满足pA1-pC1=min(pA1-pCi)的条件。The number of traffic accident predictions and the number of traffic accident corrections are obtained in turn. The number of accidents under the condition of eliminating the influence of geographic space, the influence of the geographic space, such as the adjacent areas of area 1 and area 2, match each other. Without safety optimization, the area of area 2 There is a certain number of accidents that will cause the area 1 itself. Assuming that A1 and C1 are mutually safe matching groups, the condition of pA1 -pC1 =min(pA1 -pCi ) is satisfied.
基于步骤C中过程,以A1为例,计算并获得道路A1的安全系数:Based on the process in step C, taking A1 as an example, calculate and obtain the safety factor of road A1:
根据步骤D中所描述过程,According to the procedure described in step D,
根据所获安全差距,利用优化目标道路的安全影响因素对待优化目标道路进行安全优化。According to the obtained safety gap, the safety optimization of the target road to be optimized is carried out by using the safety influencing factors of the optimized target road.
虽然本发明已以较佳实施例揭露如上,然其并非用以限定本发明。本发明所属技术领域中具有通常知识者,在不脱离本发明的精神和范围内,当可作各种的更动与润饰。因此,本发明的保护范围当视权利要求书所界定者为准。Although the present invention has been disclosed above with preferred embodiments, it is not intended to limit the present invention. Those skilled in the art to which the present invention pertains can make various changes and modifications without departing from the spirit and scope of the present invention. Therefore, the protection scope of the present invention should be determined according to the claims.
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210100208.5ACN114495501B (en) | 2022-01-27 | 2022-01-27 | Road safety analysis method based on combined geographic autoregressive matching |
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210100208.5ACN114495501B (en) | 2022-01-27 | 2022-01-27 | Road safety analysis method based on combined geographic autoregressive matching |
| Publication Number | Publication Date |
|---|---|
| CN114495501Atrue CN114495501A (en) | 2022-05-13 |
| CN114495501B CN114495501B (en) | 2022-11-15 |
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210100208.5AActiveCN114495501B (en) | 2022-01-27 | 2022-01-27 | Road safety analysis method based on combined geographic autoregressive matching |
| Country | Link |
|---|---|
| CN (1) | CN114495501B (en) |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140074402A1 (en)* | 2012-09-12 | 2014-03-13 | Lexisnexis Risk Solutions Fl Inc. | Systems and methods for determining risks associated with driving routes |
| CN107909247A (en)* | 2017-10-26 | 2018-04-13 | 东南大学 | A kind of city macroscopic view Road Traffic Safety Factors analysis method based on spatial level Bayesian model |
| CN109448369A (en)* | 2018-10-26 | 2019-03-08 | 中交第公路勘察设计研究院有限公司 | Highway real time execution Risk Calculation method |
| CN112767697A (en)* | 2021-01-19 | 2021-05-07 | 东南大学 | Traffic safety evaluation method, system and device based on traffic conflict prediction |
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140074402A1 (en)* | 2012-09-12 | 2014-03-13 | Lexisnexis Risk Solutions Fl Inc. | Systems and methods for determining risks associated with driving routes |
| CN107909247A (en)* | 2017-10-26 | 2018-04-13 | 东南大学 | A kind of city macroscopic view Road Traffic Safety Factors analysis method based on spatial level Bayesian model |
| CN109448369A (en)* | 2018-10-26 | 2019-03-08 | 中交第公路勘察设计研究院有限公司 | Highway real time execution Risk Calculation method |
| CN112767697A (en)* | 2021-01-19 | 2021-05-07 | 东南大学 | Traffic safety evaluation method, system and device based on traffic conflict prediction |
| Title |
|---|
| 牛建峰等: "区域公路网交通安全态势评估关键技术的探讨", 《上海公路》* |
| Publication number | Publication date |
|---|---|
| CN114495501B (en) | 2022-11-15 |
| Publication | Publication Date | Title |
|---|---|---|
| CN107610469B (en) | Day-dimension area traffic index prediction method considering multi-factor influence | |
| CN117523299B (en) | Image recognition method, system and storage medium based on computer network | |
| CN105427594B (en) | A kind of public transport section volume of the flow of passengers acquisition methods and system based on two-way passenger flow of getting on the bus | |
| CN109872535A (en) | A kind of current prediction technique of wisdom traffic, device and server | |
| CN111931978A (en) | Urban rail transit passenger flow state prediction method based on space-time characteristics | |
| CN103984994B (en) | Method for predicting urban rail transit passenger flow peak duration | |
| CN111105065A (en) | Rural water supply system and method based on machine learning | |
| CN108597227A (en) | Road traffic flow forecasting method under freeway toll station | |
| Wang et al. | Measure vulnerability of metro network under cascading failure | |
| CN114881381B (en) | Method and system for predicting urban water level based on improved convolutional neural network | |
| CN116504050B (en) | Special event city partitioning method based on multi-mode public transport trip data | |
| CN111881961A (en) | Power distribution network fault risk grade prediction method based on data mining | |
| CN110288139A (en) | A kind of Regional Road Network method for predicting under the event of subway in short-term | |
| CN109377755B (en) | Lane traffic efficiency optimizing grouping device and method | |
| CN106327867B (en) | Bus punctuation prediction method based on GPS data | |
| CN118350858A (en) | Method and device for predicting airline passenger quantity, electronic equipment and storage medium | |
| CN114495501B (en) | Road safety analysis method based on combined geographic autoregressive matching | |
| CN113269367B (en) | A Metro Passenger Flow Prediction Method Based on Edge Collaborative Alternating Direction Multiplier Method | |
| CN113112795B (en) | Road condition prediction method, device and equipment | |
| CN115204755B (en) | Method, device, electronic equipment, and readable storage medium for measuring and calculating service area entry rate | |
| CN117522561A (en) | Enterprise credit risk determination method and device | |
| Liu | A spatial data statistical model of urban road traffic accidents. | |
| Zhang et al. | Constrained multinomial Probit route choice modeling for passengers in large-scaled metro networks in China | |
| Liang et al. | Optimizing bus stop spacing using the simulated annealing algorithm with spatial interaction coverage model | |
| CN114757702A (en) | Method and system for constructing business demand index of virtual power plant and adapting communication mode |
| 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 | ||
| GR01 | Patent grant |