Disclosure of Invention
The invention provides a traffic safety analysis method based on an information attenuation model and driven by license plate data mining, and the traffic safety analysis method has the advantages of real-time and dynamic evaluation of road traffic safety and simple comprehensive evaluation method.
In order to achieve the purpose, the invention provides the following technical scheme:
a traffic safety analysis method based on an information attenuation model and driven by license plate data mining comprises the following steps:
s01, selecting characteristic evaluation factors for evaluating traffic safety;
and S02, quantifying each characteristic evaluation factor, and calculating the influence coefficient of each characteristic evaluation factor. The influence of historical illegal and accident factors on traffic safety is calculated by adopting an information attenuation model.
And S03, comprehensively evaluating the traffic safety level of the urban road through the influence coefficients of the characteristic evaluation factors of the plurality of characteristics.
The method is used for calculating the influence coefficient of the pre-selected characteristic evaluation factor in real time based on the license plate identification data, accident data and the like acquired in real time, and calculating the traffic safety level of the urban road by integrating various influence coefficients acquired through the information attenuation model. The evaluation method is based on the license plate identification data acquired in real time, so that the information source is accurate, the method is convenient to apply to the actual work of evaluating traffic safety, and the method can be applied in engineering.
Preferably, the characteristic evaluation factors selectable in the present invention include weather conditions, traffic flow conditions, traffic events, traffic violations, and historical traffic accidents.
Preferably, the traffic flow state comprises sudden change of flow, speed difference of road sections, proportion of road sections and vehicle types and headway; the traffic events include traffic accident events and road construction events.
According to timeliness of the data flow, the importance of information contained in the data flow is attenuated continuously along with time, so that characteristic evaluation factors distinguish a current event and a historical event, the current event needs to be calculated dynamically in real time, and the historical event needs to be calculated through a historical model. The influence coefficients of characteristic evaluation factors such as sudden change of flow, speed difference of road sections, proportion of vehicle types of road sections, headway and the like contained in the traffic flow state need to be dynamically calculated in real time.
Preferably, the step of dynamically calculating the influence coefficient of the flow rate mutation in real time comprises the following steps:
a1, calculating the rising amplitude of the sudden flow change;
real-time calculation is carried out according to a unit cycle, road section flow is calculated according to license plate identification data, and the difference between the road section flow of the current cycle and the road section flow of the previous unit cycle is larger than a threshold phi1Or the ratio of the traffic flow of the current period to the traffic flow of the previous period is larger than a threshold value phi2If so, generating sudden flow change or the flow deviates from the historical normal rule, and simultaneously calculating the rising amplitude of the sudden flow change;
b1, obtaining influence coefficients of different flow sudden changes corresponding to the rising amplitude according to calculation;
the rising amplitude is more than 0% and less than or equal to 15%, and the corresponding influence coefficient is 0.001;
the rising amplitude is more than 15 percent and less than or equal to 30 percent, and the corresponding influence coefficient is 0.125;
the rising amplitude is more than 30 percent and less than or equal to 45 percent, and the corresponding influence coefficient is 0.375;
the rising amplitude is more than 45% and less than or equal to 60%, and the corresponding influence coefficient is 0.625;
the rising amplitude is more than 60 percent, and the influence coefficient of the corresponding flow sudden change is 0.875. More specifically, the step of dynamically calculating the influence coefficient of the speed difference of the road section in real time comprises the following steps:
a2, calculating the speed difference value of the road section;
the real-time calculation is carried out in a unit cycle, the license plate number, the passing time and the position information of the vehicle are respectively acquired according to the upstream and downstream electric alarm bayonets of the road section with the known length, the same vehicle is determined through the matching of the upstream and downstream license plates, and the time required by the vehicle running is obtained through acquiring the passing time of the vehicle on the upstream and downstream electric alarm bayonets. And determining the distance traveled by the vehicle according to the position information of the vehicle, and obtaining the average speed of the road section as distance/time.
Calculating to obtain the average speed of the road section;
the speed difference value of the road section is the difference between the current period of the same road section and the average speed of the road section in the previous period;
b2, obtaining different influence coefficients corresponding to the speed difference of the road section according to calculation;
the speed difference value of the road section is more than 0km/h and less than or equal to 15km/h, and the corresponding influence coefficient is 0.001;
the speed difference value of the road section is more than 15km/h and less than or equal to 20km/h, and the corresponding influence coefficient is 0.125;
the speed difference value of the road section is more than 20km/h and less than or equal to 25km/h, and the corresponding influence coefficient is 0.375;
the speed difference value of the road section is more than 25km/h and less than or equal to 30km/h, and the corresponding influence coefficient is 0.625;
the speed difference value of the road section is more than 30km/h, and the corresponding influence coefficient is 0.875.
Preferably, the step of dynamically calculating the influence coefficient of the road section vehicle type proportion in real time comprises the following steps:
a3, calculating the ratio of the large vehicle in the daytime to the total vehicle flow in the daytime;
calculating in real time in a unit period, and obtaining the number of large vehicles in the daytime and the total flow in the daytime according to the vehicle type information acquired by the data of the electric police card port, so as to calculate the ratio of the large vehicles in the daytime to the total flow of the vehicles in the daytime;
b3, obtaining different influence coefficients corresponding to the ratio of the white-day large-sized vehicle to the daytime vehicle total flow according to calculation;
the ratio of the white-day large-sized vehicle to the total vehicle flow in the daytime is more than 0% and less than or equal to 20%, and the corresponding influence coefficient is 0.001;
the ratio of the white-day large-sized vehicles to the total vehicle flow in the daytime is more than 20% and less than or equal to 35%, and the corresponding influence coefficient is 0.125;
the ratio of the white-day large-sized vehicles to the total vehicle flow in the daytime is more than 35% and less than or equal to 50%, and the corresponding influence coefficient is 0.375;
the ratio of the white-day large-sized vehicles to the total vehicle flow in the daytime is more than 50% and less than or equal to 65%, and the corresponding influence coefficient is 0.625;
the ratio of the white-day large-sized vehicles to the total vehicle flow in the day is more than 65%, and the corresponding influence coefficient is 0.875.
In determining the impact coefficients of traffic events:
the residual traffic capacity value of the road section with the traffic incident is more than 80 percent of the pre-designed traffic capacity value, and the influence coefficient of the corresponding traffic incident is 0.125;
the residual traffic capacity value of the road section with the traffic incident is more than 50 percent of the pre-designed traffic capacity value and less than or equal to 80 percent of the pre-designed traffic capacity value, and the influence coefficient of the corresponding traffic incident is 0.375;
the residual traffic capacity value of the road section with the traffic incident is less than or equal to 50 percent of the pre-designed traffic capacity value, and the influence coefficient of the corresponding traffic incident is 0.625;
when the remaining traffic capacity value of the traffic incident road section is 0, the influence coefficient of the corresponding traffic incident is 0.875.
In determining the impact coefficients of traffic law violation:
using the mode support degree included in the factor in the information attenuation model to change along with the time t and carrying out attenuation by using an attenuation function POWER (lambda, ti/TLi), wherein, lambda is larger than 0, then:
ai=bi*POWER(0.5,ti/TLi)
wherein, aiIs shown at tiInfluence coefficient of total number of internal violation on traffic safety; TLiHalf-decay time; t is tiA time interval representing the time of occurrence of the most recent event from the current time; biRepresenting an initial coefficient;
at tiWhen the average number of illegal acts per month is less than or equal to 4 and greater than 0, the corresponding biIs 0.125;
at tiWhen the average number of illegal acts per month is less than or equal to 8 and greater than 4, the corresponding bi0.375;
at tiIn the interior of said container body,when the average number of illegal acts per month is less than or equal to 12 and greater than 8, the corresponding biIs 0.625;
at tiWhen the average number of illegal acts per month is greater than 12, the corresponding biIs 0.875.
In determining the impact factors of historical traffic accidents:
using the mode support degree included in the factor in the information attenuation model to change along with the time t and carrying out attenuation by using an attenuation function POWER (lambda, ti/TLi), wherein, lambda is larger than 0, then:
Ai=Bi*POWER(λ,Ti/TLi)
wherein A isiIs shown at TiInfluence coefficient of traffic accident to traffic safety; TLiFor half-decay time, TiA time interval representing the time of occurrence of the most recent event from the current time; b isiThe initial coefficient is represented as a function of time,
at TiWhen the average monthly traffic accident frequency is more than 1 and less than or equal to 2, the corresponding BiIs 0.25;
at TiWhen the average monthly traffic accident frequency is more than 2 and less than or equal to 4, the corresponding BiIs 0.50;
at TiWhen the average monthly traffic accident frequency is more than 4 and less than or equal to 6, the corresponding BiIs 0.75;
at TiWhen the average monthly traffic accident frequency is more than 6, the corresponding BiIs 1.
Among the factors that influence the meteorological conditions:
meteorological conditions are generally divided into: fog, rain, wind, freezing, snow. If severe weather such as heavy rain, strong wind, heavy fog, freezing and the like occurs, a1 is 0.125, otherwise a1 is 0.001;
as a preferred batch, step S03 is specifically:
according to the influence coefficient of the one-way characteristic factor, if a certain item has no score, the influence coefficient is 0.001. When a plurality of feature evaluation factors are combined, the influence of a single feature evaluation factor increases, and as a result, (0, 1):
f=min(max(ai)+avg(ai),0.99)
wherein f represents an influence coefficient of a characteristic evaluation factor of a plurality of characteristics, aiA coefficient value representing each characteristic evaluation factor;
the link safety factor decreases as the influence coefficient of the characteristic evaluation factor increases, but the speed of the decrease decreases, and as a result, (0, 10):
S=10*(1-EXP(-LOG(f,0.05)))
wherein S represents a traffic safety index for the current time slot link.
Dividing the traffic safety index into 4 traffic safety grades according to the traffic safety index, and establishing a corresponding evaluation set which is divided into { unsafe, general safe and safer } - { 0-2.5, 2.5-5, 5-7.5, 7.5-10 }.
The invention has the beneficial effects that:
the method is used for calculating the influence coefficient of the pre-selected characteristic evaluation factor in real time based on the license plate identification data, accident data and the like acquired in real time, and calculating the traffic safety level of the urban road by integrating various influence coefficients acquired through the information attenuation model. The evaluation method is based on the license plate identification data acquired in real time, so that the information source is accurate, the method is convenient to apply to the actual work of evaluating traffic safety, and the method can be applied in engineering.
example 1
A traffic safety analysis method based on an information attenuation model and driven by license plate data mining comprises the following steps:
s01, selecting weather conditions, traffic flow states, traffic incidents, traffic illegal behaviors and historical traffic accidents for evaluating characteristic evaluation factors of traffic safety; the traffic flow state comprises the following characteristic evaluation factors: flow rate sudden change, speed difference of road sections, vehicle type proportion of the road sections and vehicle head time distance; the traffic events include the following characteristic evaluation factors: traffic accident incidents and road construction incidents.
S02, calculating the influence coefficients of all the characteristic evaluation factors selected in the step S01 on the road traffic safety. See below, respectively:
calculating the influence coefficient of the sudden change of the flow on the road traffic safety: firstly, real-time calculation is carried out in a unit period, road section flow is calculated according to license plate identification data, and the difference between the road section flow in the current period and the road section flow in the previous unit period is larger than a threshold value phi1Or the ratio of the traffic flow of the current period to the traffic flow of the previous period is larger than a threshold value phi2If the flow rate suddenly changes or the flow rate deviates from the historical normal rule, the rising amplitude of the sudden flow rate change needs to be calculated. When the rising amplitude is more than 0% and less than or equal to 15%, the corresponding influence coefficient is 0.001, the rising amplitude is more than 15% and less than or equal to 30%, and the corresponding influence coefficient is 0.125; the rising amplitude is more than 30 percent and less than or equal to 45 percent, and the corresponding influence coefficient is 0.375; the rising amplitude is more than 45% and less than or equal to 60%, and the corresponding influence coefficient is 0.625; the rise amplitude is greater than 60%, and the corresponding influence coefficient is 0.875.
Calculating the influence coefficient of the speed difference of the road section on the traffic safety: the real-time calculation is carried out in a unit cycle, the license plate number, the passing time and the position information of the vehicle are respectively acquired according to the upstream and downstream electric alarm bayonets of the road section with the known length, the same vehicle is determined through the matching of the upstream and downstream license plates, and the time required by the vehicle running is obtained through acquiring the passing time of the vehicle on the upstream and downstream electric alarm bayonets. And determining the distance traveled by the vehicle according to the position information of the vehicle, and obtaining the average speed of the road section as distance/time. Thereby calculating the average speed of the road section. When the speed difference value of the road section is more than 0km/h and less than or equal to 15km/h, the corresponding influence coefficient is 0.001; the speed difference value of the road section is more than 15km/h and less than or equal to 20km/h, and the corresponding influence coefficient is 0.125; the speed difference value of the road section is more than 20km/h and less than or equal to 25km/h, and the corresponding influence coefficient is 0.375; the speed difference value of the road section is more than 25km/h and less than or equal to 30km/h, and the corresponding influence coefficient is 0.625; the speed difference value of the road section is more than 30km/h, and the corresponding influence coefficient is 0.875.
Calculating the influence coefficient of the vehicle type proportion on the road traffic safety: the real-time calculation is carried out in a unit cycle, and according to the vehicle type information acquired by the data of the electric police card port, the influence coefficient of the large-sized vehicle which runs in the daytime is calculated as an example: the number of the white-sky large-sized vehicles and the total daytime flow are calculated, and therefore the ratio of the white-sky large-sized vehicles to the total daytime vehicle flow is calculated. When the ratio of the large vehicle to the total vehicle flow in the daytime is more than 0% and less than or equal to 20%, the corresponding influence coefficient is 0.001; the ratio of the white-day large-sized vehicles to the total vehicle flow in the daytime is more than 20% and less than or equal to 35%, and the corresponding influence coefficient is 0.125; the ratio of the white-day large-sized vehicles to the total vehicle flow in the daytime is more than 35% and less than or equal to 50%, and the corresponding influence coefficient is 0.375; the ratio of the white-day large-sized vehicles to the total vehicle flow in the daytime is more than 50% and less than or equal to 65%, and the corresponding influence coefficient is 0.625; the ratio of the white-day large-sized vehicles to the total vehicle flow in the day is more than 65%, and the corresponding influence coefficient is 0.875.
Calculating the influence coefficient of the traffic event on the road safety: the influence of traffic incident factors on road traffic safety is analyzed, and the traffic incident generally comprises traffic accidents, vehicle faults, scattered objects and other incidents influencing road traffic capacity. When the residual traffic capacity value of the road section with the traffic incident is more than 50 percent of the pre-designed traffic capacity value and less than or equal to 80 percent of the pre-designed traffic capacity value, the influence coefficient of the corresponding traffic incident is 0.375; the residual traffic capacity value of the road section with the traffic incident is less than or equal to 50 percent of the pre-designed traffic capacity value, and the influence coefficient of the corresponding traffic incident is 0.625; when the remaining traffic capacity value of the traffic incident road section is 0, the influence coefficient of the corresponding traffic incident is 0.875.
Calculating the influence coefficient of the traffic violation behaviors on road safety: because the traffic violation belongs to the historical event, and the historical event needs to be calculated through the historical model, the mode support degree included by the traffic violation factor in the information attenuation model changes with the time t and is attenuated by the attenuation function POWER (λ, ti/TLi), wherein λ > 0, then:
ai=bi*POWER(0.5,ti/TLi)
wherein, aiIs shown at tiInfluence coefficient of total number of internal violation on traffic safety; TLiHalf-decay time; t is tiA time interval representing the time of occurrence of the most recent event from the current time; biRepresenting the initial coefficients. The initial coefficient values are divided into the following cases: at tiWhen the average number of illegal acts per month is less than or equal to 4 and greater than 0, the corresponding biIs 0.125; at tiWhen the average number of illegal acts per month is less than or equal to 8 and greater than 4, the corresponding bi0.375; at tiWhen the average number of illegal acts per month is less than or equal to 12 and greater than 8, the corresponding biIs 0.625; at tiWhen the average number of illegal acts per month is greater than 12, the corresponding biIs 0.875. The total number of illegal activities can be counted for nearly 3 months, and the average number of illegal activities per month can be obtained.
Calculating the influence coefficient of the traffic accident on the road safety: since the traffic accident belongs to a historical event, and the historical event needs to be calculated through a historical model, the mode support degree included by the factor in the information attenuation model changes along with the time t and is attenuated by an attenuation function POWER (lambda, ti/TLi), wherein lambda is greater than 0, then:
Ai=Bi*POWER(λ,Ti/TLi)
wherein A isiIs shown at TiInfluence coefficient of traffic accident to traffic safety; TLiFor half-decay time, TiA time interval representing the time of occurrence of the most recent event from the current time; b isiRepresenting the initial coefficients. The initial coefficient values are divided into the following cases: at TiWhen the average monthly traffic accident frequency is more than 1 and less than or equal to 2, the corresponding BiIs 0.25; at TiWhen the average monthly traffic accident frequency is more than 2 and less than or equal to 4, the corresponding BiIs 0.5; at TiWhen the average monthly traffic accident frequency is more than 4 and less than or equal to 6, the corresponding BiIs 0.75; at TiAverage number of monthly traffic accidentsWhen the number is greater than 6, corresponding BiIs 1. The total number of illegal activities can be counted for data of nearly 3 months, and the average number of traffic accidents per month can be obtained.
Calculating the influence coefficient of the meteorological condition road safety: if severe weather such as heavy rain, strong wind, heavy fog, freezing and the like occurs, a1 is 0.125, otherwise a1 is 0.001;
and S03, according to the influence coefficient of the one-way characteristic factor, if any item has no score, the influence coefficient is 0.001. When a plurality of feature evaluation factors are combined, the influence of a single feature evaluation factor increases, and as a result, (0, 1):
f=min(max(ai)+avg(ai),0.99);
wherein f represents an influence coefficient of a characteristic evaluation factor of a plurality of characteristics, aiRepresenting the influence coefficient corresponding to each characteristic evaluation factor;
the link safety factor decreases as the influence coefficient of the characteristic evaluation factor increases, but the speed of the decrease decreases, and as a result, (0, 10):
S=10*(1-EXP(-LOG(f,0.05)))
wherein S represents a traffic safety index of a road segment at the current time period, aiA coefficient value representing each characteristic evaluation factor;
dividing the traffic safety index into 4 traffic safety grades according to the traffic safety index, and establishing a corresponding evaluation set which is divided into { unsafe, general safe and safer } - { 0-2.5, 2.5-5, 5-7.5, 7.5-10 }.