CN111815973A - Signalized intersection analysis method and related equipment - Google Patents

Signalized intersection analysis method and related equipment Download PDF

Info

Publication number
CN111815973A
CN111815973A CN202010622684.4A CN202010622684A CN111815973A CN 111815973 A CN111815973 A CN 111815973A CN 202010622684 A CN202010622684 A CN 202010622684A CN 111815973 A CN111815973 A CN 111815973A
Authority
CN
China
Prior art keywords
sequence
preset
headway
signalized intersection
time period
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
CN202010622684.4A
Other languages
Chinese (zh)
Other versions
CN111815973B (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.)
Ping An International Smart City Technology Co Ltd
Original Assignee
Ping An International Smart City Technology Co Ltd
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 Ping An International Smart City Technology Co Ltd filed Critical Ping An International Smart City Technology Co Ltd
Priority to CN202010622684.4A priority Critical patent/CN111815973B/en
Publication of CN111815973A publication Critical patent/CN111815973A/en
Application granted granted Critical
Publication of CN111815973B publication Critical patent/CN111815973B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/07Controlling traffic signals
    • G08G1/081Plural intersections under common control
    • G08G1/082Controlling the time between beginning of the same phase of a cycle at adjacent intersections
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/07Controlling traffic signals
    • G08G1/081Plural intersections under common control
    • G08G1/083Controlling the allocation of time between phases of a cycle

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Traffic Control Systems (AREA)

Abstract

The invention relates to artificial intelligence and provides a signalized intersection analysis method and related equipment. The method comprises the following steps: acquiring a passing time sequence of each preset lane in each preset time period at the signalized intersection; calculating the head time distance sequence of each preset lane in each preset time period of the signalized intersection according to the passing time sequence; acquiring signal timing data of each preset lane in each preset time period at the signalized intersection; correcting and matching the signal timing data according to the passing time sequence; identifying a saturated headway sequence of each preset lane in each preset time period of the signalized intersection according to the headway sequence; removing abnormal data in the saturated headway sequence according to the corrected and matched signal timing data to obtain a target saturated headway sequence; and calculating the traffic capacity of the signalized intersection according to the target saturated headway sequence. The invention can improve the accuracy of the traffic capacity analysis of the signalized intersection.

Description

Signalized intersection analysis method and related equipment
Technical Field
The invention relates to the field of artificial intelligence, in particular to a signalized intersection analysis method and related equipment.
Background
Signal intersections serve as main components of urban road networks, and the overall operation efficiency of an urban road traffic system is greatly restricted. The traffic capacity is one of the most core indexes of the signalized intersection, and directly influences the bearing capacity of the whole traffic network. Therefore, how to correctly analyze and determine the traffic capacity of the signalized intersection is important, and accordingly, the signalized timing scheme is optimized, and the signalized intersection is designed or reconstructed, so that the traffic capacity of the signalized intersection is improved, and the traffic condition of the signalized intersection is improved.
The traffic capacity of a signal intersection is influenced by various aspects such as traffic flow composition, vehicle performance, driving environment, driving behavior, road traffic management conditions and the like, and accurate calculation of the traffic capacity is a very complicated problem. The existing signalized intersection traffic capacity analysis method mainly comprises a regression analysis method, a large amount of off-line data are adopted to carry out linear regression fitting, the real-time requirement cannot be met, the estimation precision is low, and an accurate analysis result cannot be obtained.
Disclosure of Invention
In view of the above, there is a need for a signalized intersection analysis method, device, computer device and storage medium, which can improve the accuracy of the signalized intersection traffic capacity analysis.
A first aspect of the present application provides a signalized intersection analysis method, including:
acquiring a passing time sequence of each preset lane in each preset time period at the signalized intersection;
calculating a vehicle headway sequence of each preset lane in each preset time period of the signalized intersection according to the vehicle passing time sequence;
acquiring signal timing data of each preset lane in each preset time period at the signalized intersection;
correcting and matching the signal timing data according to the passing time sequence;
identifying a saturated headway sequence of each preset lane in each preset time period of the signalized intersection according to the headway sequence;
removing abnormal data in the saturated headway sequence according to the corrected and matched signal timing data to obtain a target saturated headway sequence;
and calculating the traffic capacity of the signalized intersection according to the target saturated headway sequence.
In another possible implementation manner, the performing a correction matching on the signal timing data according to the passing time sequence includes:
mapping the passing time sequence and the signal timing data of each preset lane in each preset time period, calculating the red light running rate of each signal stage of the preset lane in the preset time period, and calculating the total red light running rate of all the signal stages of the preset lane in the preset time period according to the red light running rate of each signal stage of the preset lane in the preset time period;
and adjusting the green light turn-on time of each preset lane in each preset time period to enable the total red light running rate of the preset lane in the preset time period to obtain the minimum value.
In another possible implementation manner, the identifying, according to the headway sequence, a saturated headway sequence of each preset lane in each preset time period at the signalized intersection includes:
calculating the difference value delta e (n) of two adjacent headway distances in the headway distance sequence;
judging whether the expectation of the delta e (n) is 0 or not through hypothesis testing;
and if the expectation of delta e (n) is 0, judging that the headway sequence is a saturated headway sequence.
In another possible implementation, the determining whether the expectation of Δ e (n) is 0 by hypothesis testing includes:
whether the expectation of Δ e (n) is 0 is judged by DF hypothesis test.
In another possible implementation manner, the removing abnormal data in the saturated headway sequence according to the signal timing data includes:
judging whether each headway in the saturated headway sequence is greater than or equal to the red light duration;
and if the headway is greater than or equal to the red light duration, removing the headway from the saturated headway sequence to obtain a first filtered saturated headway sequence.
In another possible implementation manner, the removing abnormal data in the saturated headway sequence according to the signal timing data further includes:
calculating a preset quantile of the first filtered saturated headway sequence;
judging whether each headway in the first filtered saturated headway sequence is larger than or equal to the preset quantile or not;
if the headway is larger than or equal to the preset quantile, removing the headway from the first filtered saturated headway sequence to obtain a second filtered saturated headway sequence;
performing DF hypothesis test on the second filtered saturated headway sequence, and if the DF hypothesis test is passed, taking the second filtered saturated headway sequence as the target saturated headway sequence;
and if the DF hypothesis test is not passed, updating the first filtered saturated headway sequence into the second filtered saturated headway sequence and returning to the calculation of the preset quantiles of the first filtered saturated headway sequence.
In another possible implementation manner, the obtaining of the passing time sequence of each preset lane in each preset time period at the signalized intersection includes:
and acquiring the passing time sequence through an electric alarm detector of the electric alarm monitoring system.
A second aspect of the present application provides a signalized intersection analysis device, the device comprising:
the first acquisition module is used for acquiring a passing time sequence of each preset lane in each preset time period at the signalized intersection;
the first calculation module is used for calculating a vehicle headway sequence of each preset lane in each preset time period of the signalized intersection according to the vehicle passing time sequence;
the second acquisition module is used for acquiring signal timing data of each preset lane in each preset time period at the signalized intersection;
the correction module is used for correcting and matching the signal timing data according to the passing time sequence;
the identification module is used for identifying a saturated headway sequence of each preset lane in each preset time period of the signalized intersection according to the headway sequence;
the filtering module is used for removing abnormal data in the saturated headway sequence according to the corrected and matched signal timing data to obtain a target saturated headway sequence;
and the second calculation module is used for calculating the traffic capacity of the signalized intersection according to the target saturated headway sequence.
A third aspect of the application provides a computer device comprising a processor for implementing the signalized intersection analysis method when executing a computer program stored in a memory.
A fourth aspect of the present application provides a storage medium having stored thereon a computer program which, when executed by a processor, implements the signalized intersection analysis method.
The method and the device solve the problem that the existing signalized intersection communication capacity analysis method cannot obtain an accurate analysis result. The invention corrects and matches the signal timing data, corrects the time deviation of the signal timing data and improves the accuracy of the analysis of the communication capacity of the signalized intersection. The abnormal data in the saturated headway is identified and removed under the conditions of factors such as the peak preset time period, the traffic light period alternation and the like, so that the accuracy of analyzing the traffic capacity of the signalized intersection is further improved. The invention can obtain accurate signalized intersection analysis results.
Drawings
Fig. 1 is a flowchart of a signalized intersection analysis method according to an embodiment of the present invention.
Fig. 2 is a structural diagram of a signalized intersection analysis apparatus according to an embodiment of the present invention.
Fig. 3 is a schematic diagram of a computer device provided by an embodiment of the present invention.
Detailed Description
In order that the above objects, features and advantages of the present invention can be more clearly understood, a detailed description of the present invention will be given below with reference to the accompanying drawings and specific embodiments. It should be noted that the embodiments and features of the embodiments of the present application may be combined with each other without conflict.
In the following description, numerous specific details are set forth to provide a thorough understanding of the present invention, and the described embodiments are merely a subset of the embodiments of the present invention, rather than a complete embodiment. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. The terminology used in the description of the invention herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention.
Preferably, the signalized intersection analysis method of the present invention is applied to one or more computer devices. The computer device is a device capable of automatically performing numerical calculation and/or information processing according to a preset or stored instruction, and the hardware includes, but is not limited to, a microprocessor, an Application Specific Integrated Circuit (ASIC), a Programmable Gate Array (FPGA), a Digital Signal Processor (DSP), an embedded device, and the like.
The computer device can be a desktop computer, a notebook, a palm computer, a cloud server and other computing devices. The computer equipment can carry out man-machine interaction with a user through a keyboard, a mouse, a remote controller, a touch panel or voice control equipment and the like.
Example one
Fig. 1 is a flowchart of a signalized intersection analysis method according to an embodiment of the present invention. The signalized intersection analysis method is applied to computer equipment. The signalized intersection analysis method analyzes the traffic capacity of the signalized intersection.
The signalized intersection analysis method can be applied to intelligent traffic, is used for accurately analyzing the traffic capacity of signalized intersections, can optimize a signal timing scheme, design or reconstruct the signalized intersections according to the analysis result, improves the traffic capacity of the signalized intersections, and improves the traffic conditions of the signalized intersections, thereby promoting the construction of smart cities.
As shown in fig. 1, the signalized intersection analysis method includes:
101, acquiring a passing time sequence of each preset lane in each preset time period at the signalized intersection.
The preset time period and the preset lane can be set as required. In one embodiment, the predetermined time period is a predetermined time period of a traffic peak period (e.g., 7:00 to 9:00 or 17:00 to 19: 00). For example, the traffic capacity of signalized intersections at 7:00-9:00 am is analyzed, and 24 preset time periods are set, wherein each preset time period is 5 minutes. The preset lanes may be one or more lanes of a signalized intersection.
The passing time sequence of each preset lane at each preset time period at the signalized intersection is a time sequence formed by the time when each vehicle of the preset lane passes through the stop line of the signalized intersection at the preset time period. For example, if a preset time period is 7:00-7:05 and a preset lane is marked as lane 1, the time when each vehicle on the 7:00-7:05 lane 1 passes through the stop line constitutes a passing time sequence of the 7:00-7:05 lane 1.
In an embodiment, the obtaining of the passing time sequence of each preset lane at each preset time period at the signalized intersection includes: and acquiring the passing time sequence through an electric alarm detector of the electric alarm monitoring system.
The electric alarm detector is arranged above a road behind a stop line of the signalized intersection, the electric alarm detector is used for collecting the time of the vehicle in each preset lane in each preset time period passing through the stop line, and the time of the vehicle in each preset lane in each preset time period passing through the stop line is formed into the vehicle passing time sequence.
In one embodiment, a detection area is arranged on the road surface close to the stop line behind the stop line, the time when the vehicle passes through the detection area is collected through an electric alarm detector, and the time when the vehicle passes through the detection area is taken as the time when the vehicle passes through the stop line.
The electric alarm detector can be arranged above a road surface 20-30 meters behind a stop line of a signalized intersection, the electric alarm detector captures the time when a vehicle passes through the detection area, and the time when the vehicle passes through the detection area is taken as the time when the vehicle passes through the stop line.
The passing time sequence of each preset lane for each preset time period can be recorded as
Figure BDA0002563582010000071
Figure BDA0002563582010000072
Where M (M ═ 1, … M) is the slot number, N (N ═ 1, …, N +1) is the vehicle number, L (L ═ 1,2, …, L) is the lane number,
Figure BDA0002563582010000073
and the time when the nth vehicle of the ith preset lane passes through the stop line is the mth preset time period.
And 102, calculating a headway sequence of each preset lane in each preset time period of the signalized intersection according to the passing time sequence.
The sequence of the headway of each preset lane in each preset time period at the signalized intersection is a sequence formed by the headways of all vehicles passing through the stop line in the preset time period in the preset lane.
The headway refers to the time difference between adjacent vehicles passing through the stop line. The headway sequence of each preset lane in each preset time period can be recorded as
Figure BDA0002563582010000074
Where M (M-1, … M) is a preset time period, N-1, …, N is a vehicle number, L (L-1, 2, …, L) is a lane number,
Figure BDA0002563582010000075
the time headway of the nth vehicle of the ith preset lane in the mth preset time period.
The head time distance of the nth vehicle of the ith preset lane in the mth preset time period is calculated according to the following formula:
Figure BDA0002563582010000076
wherein the content of the first and second substances,
Figure BDA0002563582010000077
the time when the (n +1) th vehicle of the ith preset lane passes through the stop line for the mth preset time period,
Figure BDA0002563582010000078
and the time when the nth vehicle of the ith preset lane passes through the stop line is the mth preset time period.
And 103, acquiring signal timing data of each preset lane in each preset time period at the signalized intersection.
The signal timing data indicates the turn-on time of signal lights (red and green lights) at the signalized intersection. The signal timing data may be obtained from a signal control system.
In order to avoid mutual conflict among traffic flows in different directions, the signal control system distributes the passing time in each direction according to the traffic demands in each direction. The signal timing data of the ith preset lane for the mth preset time period may be represented as:
Figure BDA0002563582010000081
wherein the content of the first and second substances,
Figure BDA0002563582010000082
the red light turn-on time of the t signal period of the ith preset lane in the mth preset time period,
Figure BDA0002563582010000083
the green light turn-on time of the T signal cycle of the ith preset lane in the mth preset time period is T (T is 1, … T) which is a cycle number.
And 104, performing correction matching on the signal timing data according to the passing time sequence.
The acquired vehicle passing time sequence and signal timing data usually have certain time deviation, and the deviation can influence the accuracy of the analysis result of the signalized intersection.
For example, in one embodiment, the passing time sequence is obtained from the electric police monitoring system, the signal timing data is obtained from the signal control system, and the two systems respectively have respective timing mechanisms, and the clocks of the two systems have certain deviation, so that when the passing time sequence is mapped with the corresponding signal timing data, more red light running phenomena exist, that is, more time of the vehicle passing through a stop line is concentrated in a red light time range.
In order to improve the accuracy of analysis of the signalized intersection, the signal timing data of each preset lane in each preset time period is corrected and matched.
In an embodiment, the performing a modified match on the signal timing data according to the passing time sequence includes:
mapping the passing time sequence and the signal timing data of each preset lane in each preset time period, calculating the red light running rate of each signal stage of the preset lane in the preset time period, and calculating the total red light running rate of all the signal stages of the preset lane in the preset time period according to the red light running rate of each signal stage of the preset lane in the preset time period;
and adjusting the green light turn-on time of each preset lane in each preset time period to enable the total red light running rate of the preset lane in the preset time period to obtain the minimum value.
The signal stage is divided according to the changing times of the right of way at the signalized intersection in a signal period, and the right of way has several changing times in the signal period.
The total red light running rate of all signal phases is the sum of the red light running rates of all signal phases. For example, if the signal intersection includes four signal phases, the total red light running rate is the sum of the red light running rates of the first signal phase (i.e., phase 1), the second signal phase (i.e., phase 2), the third signal phase (i.e., phase 3), and the fourth signal phase (i.e., phase 4).
Adjusting the green light turn-on time of each preset lane in each preset time period to make the total red light running rate of the preset lane in the preset time period obtain the minimum value, namely performing target optimization, and can be expressed as:
Figure BDA0002563582010000091
s.t.
0≤os≤C(m)
Figure BDA0002563582010000092
wherein red _ occupancys(M) is the red light running rate of the mth preset time period signal stage S, M is 1, … M, S is 1, … S is the number of signal stages, osAnd C (m) represents the signal period duration of the mth preset time period. O is not less than 0sC (m) represents that the correction quantity of the turn-on time of the green light of each signal stage needs to be adjusted between 0 and C (m),
Figure BDA0002563582010000093
Figure BDA0002563582010000094
the sum of the correction amounts of the turn-on time of the green lamp indicating the respective signal stages is equal to 0.
And 105, identifying a saturated headway sequence of each preset lane in each preset time period of the signalized intersection according to the headway sequence.
Due to the periodic replacement of traffic light signals at the signalized intersection, the vehicles stop and go through the signalized intersection. When the red light is turned on, the vehicles close to the stop line are parked and queued, the subsequently arriving vehicles are influenced by the front queued vehicles, the speed is gradually reduced to zero, and the vehicles are added at the tail of the queuing queue. When the green light is turned on, vehicles close to the stop line start to drive away from the signalized intersection, and subsequent vehicles in line (4 th vehicle in line starts) drive away from the signalized intersection at the saturated flow rate, and the line is gradually dissipated. Therefore, from the 4 th queued vehicle to the last queued vehicle, the headway should fluctuate within the saturated headway range, and the fluctuation (error) is irregular, i.e., the errors of the headway should be independently and equally distributed, and the average value is 0.
In one embodiment, headway may be represented as:
Figure BDA0002563582010000095
wherein the content of the first and second substances,
Figure BDA0002563582010000096
the time headway of the nth vehicle of the ith preset lane in the mth preset time period is shown, mu is the saturated time headway, and e (n) is a random error term. Further, assuming that the random error term e (n) satisfies a mean of 0 and a variance of
Figure BDA0002563582010000101
The variance of the saturated headway mu and the random error term e (n)
Figure BDA0002563582010000102
The time headway of the nth vehicle of the ith preset lane in the mth preset time period is the amount to be estimated
Figure BDA0002563582010000103
Is expected to
Figure BDA0002563582010000104
Sum variance
Figure BDA0002563582010000105
Can be expressed as:
Figure BDA0002563582010000106
Figure BDA0002563582010000107
namely the acquired saturated headway sequence
Figure BDA0002563582010000108
Obedience mean is mu and variance is
Figure BDA0002563582010000109
Is normally distributed. The headway of the (n-1) th vehicle may be represented as:
Figure BDA00025635820100001010
according to
Figure BDA00025635820100001011
And
Figure BDA00025635820100001012
the relationship between two adjacent headway in the headway sequence can be established:
Figure BDA00025635820100001013
setting Δ e (n) ═ e (n) -e (n-1),
Figure BDA00025635820100001014
the expectation and variance of Δ e (n) can be expressed as:
e (Δ E (n)) ═ E (E) (n)) -E (n-1)) ═ 0 formula 8
Figure BDA00025635820100001015
According to E (Δ E (n)) ═ E (E) (n)) ═ 0 and
Figure BDA00025635820100001016
Figure BDA00025635820100001017
the head time distance error term difference delta e (n) of adjacent vehicles can be obtained, the mean value is 0, and the variance is
Figure BDA00025635820100001018
Is normally distributed.
Based on this, in an embodiment, the identifying, according to the headway sequence, the saturated headway sequence of each preset lane for each preset time period at the signalized intersection includes:
calculating the difference value delta e (n) of two adjacent headway distances in the headway distance sequence;
judging whether the expectation of the delta e (n) is 0 or not through hypothesis testing;
and if the expectation of delta e (n) is 0, judging that the headway sequence is a saturated headway sequence.
In one embodiment, the expectation of Δ e (n) is determined to be 0 by DF hypothesis testing.
For better checking the time headway obtained by observation
Figure BDA00025635820100001019
If the vehicle headway is saturated, k (n) ═ Δ e (n) is set,
Figure BDA00025635820100001020
the following expression can be obtained:
Figure BDA0002563582010000111
wherein the content of the first and second substances,
Figure BDA0002563582010000112
and p is a parameter to be estimated, and is a vehicle headway sequence of the ith preset lane in the mth preset time period. Based on the formula
Figure BDA0002563582010000113
And judging whether the detected headway is the saturated headway or not by carrying out hypothesis test on whether rho is equal to 1 or not. If passing hypothesis testing, i.e.
Figure BDA0002563582010000114
Then, E (Δ E), (n) ═ E (k) (n) ═ 0 can be obtained, and further, E (Δ E), (n) ═ E (k) (n) — 0 can be obtained
Figure BDA0002563582010000115
The observed headway sequence is the saturated headway.
Further, the parameters ρ and k (n), the estimators of which are estimated, are estimated using a least squares method
Figure BDA0002563582010000116
And
Figure BDA0002563582010000117
can be expressed as:
Figure BDA0002563582010000118
Figure BDA0002563582010000119
according to the theory of regression analysis,
Figure BDA00025635820100001110
obeying a t distribution with a degree of freedom of (N-2), wherein
Figure BDA00025635820100001111
Is composed of
Figure BDA00025635820100001112
The smaller the standard error is, the closer the acquired headway sequence is to the real saturated headway sequence is, and the closer the estimated parameters are to the real parameters is, namely, the higher the accuracy of the estimated saturated headway is.
Figure BDA00025635820100001113
The specific expression is as follows:
Figure BDA00025635820100001114
Figure BDA00025635820100001115
Figure BDA00025635820100001116
further, according to Dickey and Fuller theory, due to the estimators
Figure BDA00025635820100001117
The statistical headway of the nth and the (n-1) th vehicles is included at the same time, the hypothesis test aiming at the rho being 1 does not satisfy the t test any more, but satisfies the DF test, and the DF test statistic is constructed:
Figure BDA00025635820100001118
the DF statistic constructed above obeys DF distribution, and the critical values of different confidence degrees of the DF distribution are met during hypothesis test. In the invention, when the sample amount of the headway sequence is more than 50 and the DF test statistic is [ -2.25,1.66], the test rho is assumed to be not rejected as 1, namely, the observed headway sequence mean value can be regarded as saturated headway data.
And 106, removing abnormal data in the saturated headway sequence according to the corrected and matched signal timing data to obtain a target saturated headway sequence.
In the above steps, the time headway sequence checked by the DF hypothesis has a potential abnormal time headway due to factors such as a preset time period of receiving a peak, a period alternation of traffic lights, and the like. Therefore, the abnormal headway needs to be filtered, and the estimation precision of the saturated headway is ensured. Common abnormal headway includes the following situations: (1) due to the red light time, the time interval between the vehicle passing through the stop line in the first period and the vehicle passing through the stop line in the last period is more than or equal to the red light time; (2) the green light is turned on, the vehicle is accelerated to start and leave a stop line, and the time interval of the vehicle head is larger due to the fact that the vehicle is still in an accelerated state; (3) during green light, there may also be a large headway due to random arrival of vehicles, etc.
In an embodiment, the removing abnormal data in the saturated headway sequence according to the modified and matched signal timing data includes:
judging whether each headway in the saturated headway sequence is greater than or equal to the red light duration;
and if the headway is greater than or equal to the red light duration, removing the headway from the saturated headway sequence to obtain a first filtered saturated headway sequence.
In an embodiment, removing abnormal data in the saturated headway sequence according to the modified and matched signal timing data further includes:
calculating a preset quantile (e.g., 85% quantile) of the first filtered saturated headway sequence;
judging whether each headway in the first filtered saturated headway sequence is larger than or equal to the preset quantile or not;
if the headway is larger than or equal to the preset quantile, removing the headway from the first filtered saturated headway sequence to obtain a second filtered saturated headway sequence;
performing DF hypothesis test on the second filtered saturated headway sequence, and if the DF hypothesis test is passed, taking the second filtered saturated headway sequence as a target saturated headway sequence;
and if the DF hypothesis test is not passed, updating the first filtered saturated headway sequence into the second filtered saturated headway sequence and returning to the calculation of the preset quantiles of the first filtered saturated headway sequence.
And 107, calculating the traffic capacity of the signalized intersection according to the target saturated headway.
In an embodiment, an average value of the target saturated headway sequence of each preset lane in each preset time period may be calculated, and the average value is used as the target saturated headway sequence of the preset lane in the preset time period
Figure BDA0002563582010000131
According to
Figure BDA0002563582010000132
The traffic capacity of the preset lane can be obtained:
Figure BDA0002563582010000133
wherein the content of the first and second substances,
Figure BDA0002563582010000134
the traffic capacity of a first preset lane m in a preset time period, S is the saturation flow, tegThe effective green duration.
The unit of traffic capacity is pcu/h,
Figure BDA0002563582010000135
the larger the value of (a), the stronger the trafficability of the signalized intersection. On the contrary, the method can be used for carrying out the following steps,
Figure BDA0002563582010000136
the smaller the value of (a), the weaker the traffic capacity of the signalized intersection.
In other embodimentsIn other ways, the target saturated headway sequence can be calculated according to the target saturated headway sequence
Figure BDA0002563582010000137
For example, a median of the target saturated headway sequence of each preset lane in each preset time period may be calculated, and the median is used as the target saturated headway of the preset lane in the preset time period
Figure BDA0002563582010000138
According to
Figure BDA0002563582010000139
And calculating the traffic capacity of the preset lane.
The signalized intersection analysis method solves the problem that an accurate analysis result cannot be obtained by an existing signalized intersection communication capacity analysis method. The signal intersection analysis method is used for correcting and matching the signal timing data, correcting the time deviation of the signal timing data and improving the accuracy of the analysis of the communication capacity of the signal intersection. And abnormal data exist in the saturated head time interval under the factors of the peak preset time period, the traffic light period alternation and the like, and the signalized intersection analysis method identifies the abnormal data in the saturated head time interval and eliminates the abnormal data, so that the accuracy of the signalized intersection traffic capacity analysis is further improved. The signalized intersection analysis method can obtain an accurate signalized intersection analysis result.
Example two
Fig. 2 is a structural diagram of a signalized intersection analysis apparatus according to a second embodiment of the present invention. The signalized intersection analysis device 20 is applied to a computer device. The signalized intersection analysis device 20 analyzes the traffic capacity of the signalized intersection.
The signalized intersection analysis device 20 can be applied to intelligent traffic, is used for accurately analyzing the traffic capacity of signalized intersections, can optimize a signal timing scheme, design or reconstruct signalized intersections according to analysis results, improves the traffic capacity of signalized intersections, improves the traffic conditions of signalized intersections, and accordingly promotes the construction of smart cities.
As shown in fig. 2, the signalized intersection analyzing apparatus 20 may include a first obtaining module 201, a first calculating module 202, a second obtaining module 203, a correcting module 204, an identifying module 205, a filtering module 206, and a second calculating module 207.
The first obtaining module 201 is configured to obtain a passing time sequence of each preset lane in each preset time period at the signalized intersection.
The preset time period and the preset lane can be set as required. In one embodiment, the predetermined time period is a predetermined time period of a traffic peak period (e.g., 7:00 to 9:00 or 17:00 to 19: 00). For example, the traffic capacity of signalized intersections at 7:00-9:00 am is analyzed, and 24 preset time periods are set, wherein each preset time period is 5 minutes. The preset lanes may be one or more lanes of a signalized intersection.
The passing time sequence of each preset lane at each preset time period at the signalized intersection is a time sequence formed by the time when each vehicle of the preset lane passes through the stop line of the signalized intersection at the preset time period. For example, if a preset time period is 7:00-7:05 and a preset lane is marked as lane 1, the time when each vehicle on the 7:00-7:05 lane 1 passes through the stop line constitutes a passing time sequence of the 7:00-7:05 lane 1.
In an embodiment, the obtaining of the passing time sequence of each preset lane at each preset time period at the signalized intersection includes: and acquiring the passing time sequence through an electric alarm detector of the electric alarm monitoring system.
The electric alarm detector is arranged above a road behind a stop line of the signalized intersection, the electric alarm detector is used for collecting the time of the vehicle in each preset lane in each preset time period passing through the stop line, and the time of the vehicle in each preset lane in each preset time period passing through the stop line is formed into the vehicle passing time sequence.
In one embodiment, a detection area is arranged on the road surface close to the stop line behind the stop line, the time when the vehicle passes through the detection area is collected through an electric alarm detector, and the time when the vehicle passes through the detection area is taken as the time when the vehicle passes through the stop line.
The electric alarm detector can be arranged above a road surface 20-30 meters behind a stop line of a signalized intersection, the electric alarm detector captures the time when a vehicle passes through the detection area, and the time when the vehicle passes through the detection area is taken as the time when the vehicle passes through the stop line.
The passing time sequence of each preset lane for each preset time period can be recorded as
Figure BDA0002563582010000151
Figure BDA0002563582010000152
Where M (M ═ 1, … M) is the slot number, N (N ═ 1, …, N +1) is the vehicle number, L (L ═ 1,2, …, L) is the lane number,
Figure BDA0002563582010000153
and the time when the nth vehicle of the ith preset lane passes through the stop line is the mth preset time period.
The first calculating module 202 is configured to calculate a headway sequence of each preset lane in each preset time period of the signalized intersection according to the passing time sequence.
The sequence of the headway of each preset lane in each preset time period at the signalized intersection is a sequence formed by the headways of all vehicles passing through the stop line in the preset time period in the preset lane.
The headway refers to the time difference between adjacent vehicles passing through the stop line. The headway sequence of each preset lane in each preset time period can be recorded as
Figure BDA0002563582010000154
Where M (M-1, … M) is a preset time period, N-1, …, N is a vehicle number, L (L-1, 2, …, L) is a lane number,
Figure BDA0002563582010000155
the time headway of the nth vehicle of the ith preset lane in the mth preset time period.
The head time distance of the nth vehicle of the ith preset lane in the mth preset time period is calculated according to the following formula:
Figure BDA0002563582010000156
wherein the content of the first and second substances,
Figure BDA0002563582010000157
the time when the (n +1) th vehicle of the ith preset lane passes through the stop line for the mth preset time period,
Figure BDA0002563582010000158
and the time when the nth vehicle of the ith preset lane passes through the stop line is the mth preset time period.
The second obtaining module 203 is configured to obtain signal timing data of each preset lane at each preset time period at the signalized intersection.
The signal timing data indicates the turn-on time of signal lights (red and green lights) at the signalized intersection. The signal timing data may be obtained from a signal control system.
In order to avoid mutual conflict among traffic flows in different directions, the signal control system distributes the passing time in each direction according to the traffic demands in each direction. The signal timing data of the ith preset lane for the mth preset time period may be represented as:
Figure BDA0002563582010000161
wherein the content of the first and second substances,
Figure BDA0002563582010000162
the red light turn-on time of the t signal period of the ith preset lane in the mth preset time period,
Figure BDA0002563582010000163
is as followsThe green light turn-on time of the T signal period of the ith preset lane in the m preset time periods is the period number T (T is 1, … T).
And the correcting module 204 is configured to perform correction matching on the signal timing data according to the passing time sequence.
The acquired vehicle passing time sequence and signal timing data usually have certain time deviation, and the deviation can influence the accuracy of the analysis result of the signalized intersection.
For example, in one embodiment, the passing time sequence is obtained from the electric police monitoring system, the signal timing data is obtained from the signal control system, and the two systems respectively have respective timing mechanisms, and the clocks of the two systems have certain deviation, so that when the passing time sequence is mapped with the corresponding signal timing data, more red light running phenomena exist, that is, more time of the vehicle passing through a stop line is concentrated in a red light time range.
In order to improve the accuracy of analysis of the signalized intersection, the signal timing data of each preset lane in each preset time period is corrected and matched.
In an embodiment, the performing a modified match on the signal timing data according to the passing time sequence includes:
mapping the passing time sequence and the signal timing data of each preset lane in each preset time period, calculating the red light running rate of each signal stage of the preset lane in the preset time period, and calculating the total red light running rate of all the signal stages of the preset lane in the preset time period according to the red light running rate of each signal stage of the preset lane in the preset time period;
and adjusting the green light turn-on time of each preset lane in each preset time period to enable the total red light running rate of the preset lane in the preset time period to obtain the minimum value.
The signal stage is divided according to the changing times of the right of way at the signalized intersection in a signal period, and the right of way has several changing times in the signal period.
The total red light running rate of all signal phases is the sum of the red light running rates of all signal phases. For example, if the signal intersection includes four signal phases, the total red light running rate is the sum of the red light running rates of the first signal phase (i.e., phase 1), the second signal phase (i.e., phase 2), the third signal phase (i.e., phase 3), and the fourth signal phase (i.e., phase 4).
Adjusting the green light turn-on time of each preset lane in each preset time period to make the total red light running rate of the preset lane in the preset time period obtain the minimum value, namely performing target optimization, and can be expressed as:
Figure BDA0002563582010000171
s.t.
0≤os≤C(m)
Figure BDA0002563582010000172
wherein red _ occupancys(M) is the red light running rate of the mth preset time period signal stage S, M is 1, … M, S is 1, … S is the number of signal stages, osAnd C (m) represents the signal period duration of the mth preset time period. O is not less than 0sC (m) represents that the correction quantity of the turn-on time of the green light of each signal stage needs to be adjusted between 0 and C (m),
Figure BDA0002563582010000173
Figure BDA0002563582010000174
the sum of the correction amounts of the turn-on time of the green lamp indicating the respective signal stages is equal to 0.
The identifying module 205 is configured to identify a saturated headway sequence of each preset lane at each preset time period of the signalized intersection according to the headway sequence.
Due to the periodic replacement of traffic light signals at the signalized intersection, the vehicles stop and go through the signalized intersection. When the red light is turned on, the vehicles close to the stop line are parked and queued, the subsequently arriving vehicles are influenced by the front queued vehicles, the speed is gradually reduced to zero, and the vehicles are added at the tail of the queuing queue. When the green light is turned on, vehicles close to the stop line start to drive away from the signalized intersection, and subsequent vehicles in line (4 th vehicle in line starts) drive away from the signalized intersection at the saturated flow rate, and the line is gradually dissipated. Therefore, from the 4 th queued vehicle to the last queued vehicle, the headway should fluctuate within the saturated headway range, and the fluctuation (error) is irregular, i.e., the errors of the headway should be independently and equally distributed, and the average value is 0.
In one embodiment, headway may be represented as:
Figure BDA0002563582010000181
wherein the content of the first and second substances,
Figure BDA0002563582010000182
the time headway of the nth vehicle of the ith preset lane in the mth preset time period is shown, mu is the saturated time headway, and e (n) is a random error term. Further, assuming that the random error term e (n) satisfies a mean of 0 and a variance of
Figure BDA0002563582010000183
The variance of the saturated headway mu and the random error term e (n)
Figure BDA0002563582010000184
The time headway of the nth vehicle of the ith preset lane in the mth preset time period is the amount to be estimated
Figure BDA0002563582010000185
Is expected to
Figure BDA0002563582010000186
Sum variance
Figure BDA0002563582010000187
Can be expressed as:
Figure BDA0002563582010000188
Figure BDA0002563582010000189
namely the acquired saturated headway sequence
Figure BDA00025635820100001810
Obedience mean is mu and variance is
Figure BDA00025635820100001811
Is normally distributed. The headway of the (n-1) th vehicle may be represented as:
Figure BDA00025635820100001812
according to
Figure BDA00025635820100001813
And
Figure BDA00025635820100001814
the relationship between two adjacent headway in the headway sequence can be established:
Figure BDA00025635820100001815
setting Δ e (n) ═ e (n) -e (n-1),
Figure BDA00025635820100001816
the expectation and variance of Δ e (n) can be expressed as:
e (Δ E (n)) ═ E (E) (n)) -E (n-1)) ═ 0 formula 8
Figure BDA00025635820100001817
According to E (Δ E (n)) ═ E (E) (n)) ═ 0 and
Figure BDA00025635820100001818
Figure BDA00025635820100001819
the head time distance error term difference delta e (n) of adjacent vehicles can be obtained, the mean value is 0, and the variance is
Figure BDA0002563582010000191
Is normally distributed.
Based on this, in an embodiment, the identifying, according to the headway sequence, the saturated headway sequence of each preset lane for each preset time period at the signalized intersection includes:
calculating the difference value delta e (n) of two adjacent headway distances in the headway distance sequence;
judging whether the expectation of the delta e (n) is 0 or not through hypothesis testing;
and if the expectation of delta e (n) is 0, judging that the headway sequence is a saturated headway sequence.
In one embodiment, the expectation of Δ e (n) is determined to be 0 by DF hypothesis testing.
For better checking the time headway obtained by observation
Figure BDA0002563582010000192
If the vehicle headway is saturated, k (n) ═ Δ e (n) is set,
Figure BDA0002563582010000193
the following expression can be obtained:
Figure BDA0002563582010000194
wherein the content of the first and second substances,
Figure BDA0002563582010000195
and p is a parameter to be estimated, and is a vehicle headway sequence of the ith preset lane in the mth preset time period. Based on the formula
Figure BDA0002563582010000196
Judging whether the detected headway is 1 by carrying out hypothesis test on whether rho is equal to 1Or else, the head time interval is saturated. If passing hypothesis testing, i.e.
Figure BDA0002563582010000197
Then, E (Δ E), (n) ═ E (k) (n) ═ 0 can be obtained, and further, E (Δ E), (n) ═ E (k) (n) — 0 can be obtained
Figure BDA0002563582010000198
The observed headway sequence is the saturated headway.
Further, the parameters ρ and k (n), the estimators of which are estimated, are estimated using a least squares method
Figure BDA0002563582010000199
And
Figure BDA00025635820100001910
can be expressed as:
Figure BDA00025635820100001911
Figure BDA00025635820100001912
according to the theory of regression analysis,
Figure BDA00025635820100001913
obeying a t distribution with a degree of freedom of (N-2), wherein
Figure BDA00025635820100001914
Is composed of
Figure BDA00025635820100001915
The smaller the standard error is, the closer the acquired headway sequence is to the real saturated headway sequence is, and the closer the estimated parameters are to the real parameters is, namely, the higher the accuracy of the estimated saturated headway is.
Figure BDA00025635820100001916
The specific expression is as follows:
Figure BDA00025635820100001917
Figure BDA00025635820100001918
Figure BDA0002563582010000201
further, according to Dickey and Fuller theory, due to the estimators
Figure BDA0002563582010000202
The statistical headway of the nth and the (n-1) th vehicles is included at the same time, the hypothesis test aiming at the rho being 1 does not satisfy the t test any more, but satisfies the DF test, and the DF test statistic is constructed:
Figure BDA0002563582010000203
the DF statistic constructed above obeys DF distribution, and the critical values of different confidence degrees of the DF distribution are met during hypothesis test. In the invention, when the sample amount of the headway sequence is more than 50 and the DF test statistic is [ -2.25,1.66], the test rho is assumed to be not rejected as 1, namely, the observed headway sequence mean value can be regarded as saturated headway data.
And the filtering module 206 is configured to remove abnormal data in the saturated headway sequence according to the modified and matched signal timing data to obtain a target saturated headway sequence.
In the above steps, the time headway sequence checked by the DF hypothesis has a potential abnormal time headway due to factors such as a preset time period of receiving a peak, a period alternation of traffic lights, and the like. Therefore, the abnormal headway needs to be filtered, and the estimation precision of the saturated headway is ensured. Common abnormal headway includes the following situations: (1) due to the red light time, the time interval between the vehicle passing through the stop line in the first period and the vehicle passing through the stop line in the last period is more than or equal to the red light time; (2) the green light is turned on, the vehicle is accelerated to start and leave a stop line, and the time interval of the vehicle head is larger due to the fact that the vehicle is still in an accelerated state; (3) during green light, there may also be a large headway due to random arrival of vehicles, etc.
In an embodiment, the removing abnormal data in the saturated headway sequence according to the modified and matched signal timing data includes:
judging whether each headway in the saturated headway sequence is greater than or equal to the red light duration;
and if the headway is greater than or equal to the red light duration, removing the headway from the saturated headway sequence to obtain a first filtered saturated headway sequence.
In an embodiment, removing abnormal data in the saturated headway sequence according to the modified and matched signal timing data further includes:
calculating a preset quantile (e.g., 85% quantile) of the first filtered saturated headway sequence;
judging whether each headway in the first filtered saturated headway sequence is larger than or equal to the preset quantile or not;
if the headway is larger than or equal to the preset quantile, removing the headway from the first filtered saturated headway sequence to obtain a second filtered saturated headway sequence;
performing DF hypothesis test on the second filtered saturated headway sequence, and if the DF hypothesis test is passed, taking the second filtered saturated headway sequence as a target saturated headway sequence;
and if the DF hypothesis test is not passed, updating the first filtering saturated headway sequence into the second filtering saturated headway sequence by an iteration factor i +1, and returning to the step of calculating the preset quantile of the first filtering saturated headway sequence.
And the second calculating module 207 is used for calculating the traffic capacity of the signalized intersection according to the target saturated headway.
In an embodiment, an average value of the target saturated headway sequence of each preset lane in each preset time period may be calculated, and the average value is used as the target saturated headway sequence of the preset lane in the preset time period
Figure BDA0002563582010000211
According to
Figure BDA0002563582010000212
The traffic capacity of the preset lane can be obtained:
Figure BDA0002563582010000213
wherein the content of the first and second substances,
Figure BDA0002563582010000214
the traffic capacity of a first preset lane m in a preset time period, S is the saturation flow, tegThe effective green duration.
The unit of traffic capacity is pcu/h,
Figure BDA0002563582010000215
the larger the value of (a), the stronger the trafficability of the signalized intersection. On the contrary, the method can be used for carrying out the following steps,
Figure BDA0002563582010000216
the smaller the value of (a), the weaker the traffic capacity of the signalized intersection.
In other embodiments, other ways of calculating the target saturated headway from the target saturated headway sequence may be employed
Figure BDA0002563582010000217
For example, a median of the target saturated headway sequence of each preset lane in each preset time period may be calculated, and the median is used as the target saturated headway of the preset lane in the preset time period
Figure BDA0002563582010000218
According to
Figure BDA0002563582010000219
And calculating the traffic capacity of the preset lane.
The signalized intersection analysis device 20 of the second embodiment solves the problem that the existing signalized intersection communication capacity analysis method cannot obtain an accurate analysis result. The passing time sequence and the signal timing data usually have certain time deviation, the deviation can influence the accuracy of the analysis result of the signalized intersection, and the signalized intersection analysis device 20 corrects and matches the signal timing data, corrects the time deviation of the signal timing data and improves the accuracy of the analysis of the communication capacity of the signalized intersection. Due to factors such as a peak preset time period and traffic light period alternation, abnormal data can exist in the saturated head time interval, the signalized intersection analysis device 20 identifies the abnormal data in the saturated head time interval and eliminates the abnormal data, and the accuracy of analyzing the traffic capacity of the signalized intersection is further improved. The signalized intersection analysis device 20 can obtain an accurate signalized intersection analysis result.
EXAMPLE III
The present embodiment provides a storage medium, where a computer program is stored on the storage medium, and when the computer program is executed by a processor, the computer program implements the steps in the above-mentioned signalized intersection analysis method embodiment, for example, 101-107 shown in fig. 1:
101, acquiring a passing time sequence of each preset lane in each preset time period at the signalized intersection;
102, calculating a vehicle headway sequence of each preset lane in each preset time period of the signalized intersection according to the vehicle passing time sequence;
103, acquiring signal timing data of each preset lane in each preset time period at the signalized intersection;
104, correcting and matching the signal timing data according to the passing time sequence;
105, identifying a saturated headway sequence of each preset lane in each preset time period of the signalized intersection according to the headway sequence;
106, removing abnormal data in the saturated headway sequence according to the corrected and matched signal timing data to obtain a target saturated headway sequence;
and 107, calculating the traffic capacity of the signalized intersection according to the target saturated headway sequence.
Alternatively, the computer program, when executed by the processor, implements the functions of the modules in the above device embodiments, for example, the module 201 and 207 in fig. 2:
the first obtaining module 201 is configured to obtain a passing time sequence of each preset lane in each preset time period at the signalized intersection;
the first calculating module 202 is configured to calculate a headway sequence of each preset lane in each preset time period of the signalized intersection according to the passing time sequence;
the second obtaining module 203 is configured to obtain signal timing data of each preset lane at each preset time period at the signalized intersection;
a correction module 204, configured to perform correction matching on the signal timing data according to the passing time sequence;
the identification module 205 is configured to identify a saturated headway sequence of each preset lane in each preset time period of the signalized intersection according to the headway sequence;
the filtering module 206 is configured to remove abnormal data in the saturated headway sequence according to the modified and matched signal timing data to obtain a target saturated headway sequence;
and the second calculating module 207 is used for calculating the traffic capacity of the signalized intersection according to the target saturated headway sequence.
Example four
Fig. 3 is a schematic diagram of a computer device according to a fourth embodiment of the present invention. The computer device 30 comprises a memory 301, a processor 302 and a computer program 303, such as a signalized intersection analysis program, stored in the memory 301 and executable on the processor 302. The processor 302 executes the computer program 303 to implement the steps in the intersection analysis method embodiment, such as 101-107 shown in fig. 1. Alternatively, the computer program, when executed by the processor, implements the functions of the modules in the above-described device embodiments, such as the module 201 and 207 in fig. 2.
Illustratively, the computer program 303 may be partitioned into one or more modules that are stored in the memory 301 and executed by the processor 302 to perform the present method. The one or more modules may be a series of computer program instruction segments capable of performing specific functions, which are used to describe the execution of the computer program 303 in the computer device 30. For example, the computer program 303 may be divided into modules in fig. 2.
The computer device 30 may be a desktop computer, a notebook, a palm computer, a cloud server, or other computing devices. Those skilled in the art will appreciate that the schematic diagram 3 is merely an example of the computer device 30 and does not constitute a limitation of the computer device 30, and may include more or less components than those shown, or combine certain components, or different components, for example, the computer device 30 may also include input and output devices, network access devices, buses, etc.
The Processor 302 may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other Programmable logic device, discrete Gate or transistor logic, discrete hardware components, etc. A general purpose processor may be a microprocessor or the processor 302 may be any conventional processor or the like, the processor 302 being the control center for the computer device 30 and connecting the various parts of the overall computer device 30 using various interfaces and lines.
The memory 301 may be used to store the computer program 303, and the processor 302 may implement various functions of the computer device 30 by running or executing the computer program or module stored in the memory 301 and calling data stored in the memory 301. The memory 301 may mainly include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required by at least one function (such as a sound playing function, an image playing function, etc.), and the like; the storage data area may store data created according to the use of the computer device 30. Further, the memory 301 may include a non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), at least one magnetic disk storage device, a Flash memory device, or other non-volatile solid state storage device.
The modules integrated by the computer device 30 may be stored in a storage medium if they are implemented in the form of software functional modules and sold or used as separate products. Based on such understanding, all or part of the flow of the method according to the embodiments of the present invention may also be implemented by a computer program, which may be stored in a storage medium and executed by a processor, to instruct related hardware to implement the steps of the embodiments of the method. Wherein the computer program comprises computer program code, which may be in the form of source code, object code, an executable file or some intermediate form, etc. The computer-readable medium may include: any entity or device capable of carrying said computer program code, recording medium, U-disk, removable hard disk, magnetic disk, optical disk, computer Memory, Read-Only Memory (ROM).
Further, the computer usable storage medium may mainly include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required for at least one function, and the like; the storage data area may store data created according to the use of the blockchain node, and the like.
The block chain is a novel application mode of computer technologies such as distributed data storage, point-to-point transmission, a consensus mechanism, an encryption algorithm and the like. A block chain (Blockchain), which is essentially a decentralized database, is a series of data blocks associated by using a cryptographic method, and each data block contains information of a batch of network transactions, so as to verify the validity (anti-counterfeiting) of the information and generate a next block. The blockchain may include a blockchain underlying platform, a platform product service layer, an application service layer, and the like.
In the embodiments provided in the present invention, it should be understood that the disclosed system, apparatus and method may be implemented in other ways. For example, the above-described apparatus embodiments are merely illustrative, and for example, the division of the modules is only one logical functional division, and other divisions may be realized in practice.
The modules described as separate parts may or may not be physically separate, and parts displayed as modules may or may not be physical modules, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment.
In addition, functional modules in the embodiments of the present invention may be integrated into one processing module, or each of the modules may exist alone physically, or two or more modules are integrated into one module. The integrated module can be realized in a hardware form, and can also be realized in a form of hardware and a software functional module.
The integrated module implemented in the form of a software functional module may be stored in a storage medium. The software functional module is stored in a storage medium and includes several instructions to enable a computer device (which may be a personal computer, a server, or a network device) or a processor (processor) to execute some steps of the methods according to the embodiments of the present invention.
It will be evident to those skilled in the art that the invention is not limited to the details of the foregoing illustrative embodiments, and that the present invention may be embodied in other specific forms without departing from the spirit or essential attributes thereof. The present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein. Any reference signs in the claims shall not be construed as limiting the claim concerned. Furthermore, it is to be understood that the word "comprising" does not exclude other modules or steps, and the singular does not exclude the plural. A plurality of modules or means recited in the system claims may also be implemented by one module or means in software or hardware. The terms first, second, etc. are used to denote names, but not any particular order.
Finally, it should be noted that the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting, and although the present invention is described in detail with reference to the preferred embodiments, it should be understood by those skilled in the art that modifications or equivalent substitutions may be made on the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention.

Claims (10)

1. A signalized intersection analysis method, the method comprising:
acquiring a passing time sequence of each preset lane in each preset time period at the signalized intersection;
calculating a vehicle headway sequence of each preset lane in each preset time period of the signalized intersection according to the vehicle passing time sequence;
acquiring signal timing data of each preset lane in each preset time period at the signalized intersection;
correcting and matching the signal timing data according to the passing time sequence;
identifying a saturated headway sequence of each preset lane in each preset time period of the signalized intersection according to the headway sequence;
removing abnormal data in the saturated headway sequence according to the corrected and matched signal timing data to obtain a target saturated headway sequence;
and calculating the traffic capacity of the signalized intersection according to the target saturated headway sequence.
2. The signalized intersection analysis method according to claim 1, wherein the modifying and matching of the signalized timing data according to the passing time sequence comprises:
mapping the passing time sequence and the signal timing data of each preset lane in each preset time period, calculating the red light running rate of each signal stage of the preset lane in the preset time period, and calculating the total red light running rate of all the signal stages of the preset lane in the preset time period according to the red light running rate of each signal stage of the preset lane in the preset time period;
and adjusting the green light turn-on time of each preset lane in each preset time period to enable the total red light running rate of the preset lane in the preset time period to obtain the minimum value.
3. The signalized intersection analysis method according to claim 1, wherein the identifying the saturated headway sequence for each preset lane for each preset time period at the signalized intersection according to the headway sequence comprises:
calculating the difference value delta e (n) of two adjacent headway distances in the headway distance sequence;
judging whether the expectation of the delta e (n) is 0 or not through hypothesis testing;
and if the expectation of delta e (n) is 0, judging that the headway sequence is a saturated headway sequence.
4. The signalized intersection analysis method of claim 3, wherein the determining whether the expectation of Δ e (n) is 0 by hypothesis testing comprises:
whether the expectation of Δ e (n) is 0 is judged by DF hypothesis test.
5. The signalized intersection analysis method according to claim 1, wherein the removing of the abnormal data in the saturated headway sequence according to the signalized timing data comprises:
judging whether each headway in the saturated headway sequence is greater than or equal to the red light duration;
and if the headway is greater than or equal to the red light duration, removing the headway from the saturated headway sequence to obtain a first filtered saturated headway sequence.
6. The signalized intersection analysis method according to claim 5, wherein the removing of the abnormal data in the saturated headway sequence according to the signalized timing data further comprises:
calculating a preset quantile of the first filtered saturated headway sequence;
judging whether each headway in the first filtered saturated headway sequence is larger than or equal to the preset quantile or not;
if the headway is larger than or equal to the preset quantile, removing the headway from the first filtered saturated headway sequence to obtain a second filtered saturated headway sequence;
performing DF hypothesis test on the second filtered saturated headway sequence, and if the DF hypothesis test is passed, taking the second filtered saturated headway sequence as the target saturated headway sequence;
and if the DF hypothesis test is not passed, updating the first filtered saturated headway sequence into the second filtered saturated headway sequence and returning to the step of calculating the preset quantiles of the first filtered saturated headway sequence.
7. The signalized intersection analysis method according to claim 1, wherein the obtaining of the passing time sequence of each preset lane for each preset time period at the signalized intersection comprises:
and acquiring the passing time sequence through an electric alarm detector of the electric alarm monitoring system.
8. A signalized intersection analysis device, the device comprising:
the first acquisition module is used for acquiring a passing time sequence of each preset lane in each preset time period at the signalized intersection;
the first calculation module is used for calculating a vehicle headway sequence of each preset lane in each preset time period of the signalized intersection according to the vehicle passing time sequence;
the second acquisition module is used for acquiring signal timing data of each preset lane in each preset time period at the signalized intersection;
the correction module is used for correcting and matching the signal timing data according to the passing time sequence;
the identification module is used for identifying a saturated headway sequence of each preset lane in each preset time period of the signalized intersection according to the headway sequence;
the filtering module is used for removing abnormal data in the saturated headway sequence according to the corrected and matched signal timing data to obtain a target saturated headway sequence;
and the second calculation module is used for calculating the traffic capacity of the signalized intersection according to the target saturated headway sequence.
9. A computer device comprising a processor for executing a computer program stored in a memory to implement a signalized intersection analysis method according to any one of claims 1 to 7.
10. A storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the signalized intersection analysis method according to any one of claims 1 to 7.
CN202010622684.4A 2020-06-30 2020-06-30 Signalized intersection analysis method and related equipment Active CN111815973B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010622684.4A CN111815973B (en) 2020-06-30 2020-06-30 Signalized intersection analysis method and related equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010622684.4A CN111815973B (en) 2020-06-30 2020-06-30 Signalized intersection analysis method and related equipment

Publications (2)

Publication Number Publication Date
CN111815973A true CN111815973A (en) 2020-10-23
CN111815973B CN111815973B (en) 2022-12-16

Family

ID=72855922

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010622684.4A Active CN111815973B (en) 2020-06-30 2020-06-30 Signalized intersection analysis method and related equipment

Country Status (1)

Country Link
CN (1) CN111815973B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113129594A (en) * 2021-04-23 2021-07-16 武汉科技大学 Intersection saturation flow rate calculation method based on bayonet data
CN113643531A (en) * 2021-07-20 2021-11-12 东北大学 Intersection lane saturation flow rate calculation method based on small time zone division statistics
CN113936453A (en) * 2021-09-09 2022-01-14 上海宝康电子控制工程有限公司 Information identification method and system based on headway

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102509456A (en) * 2011-11-21 2012-06-20 青岛海信网络科技股份有限公司 Saturation flow determination method and device
CN107564290A (en) * 2017-10-13 2018-01-09 公安部交通管理科学研究所 A kind of urban road intersection saturation volume rate computational methods
CN110070734A (en) * 2019-05-14 2019-07-30 东南大学 Signalized intersections saturation headway estimation method based on gauss hybrid models
CN110100271A (en) * 2016-11-29 2019-08-06 赫尔环球有限公司 For using the method, apparatus and computer program product of data traffic signal estimating road traffic situation

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102509456A (en) * 2011-11-21 2012-06-20 青岛海信网络科技股份有限公司 Saturation flow determination method and device
CN110100271A (en) * 2016-11-29 2019-08-06 赫尔环球有限公司 For using the method, apparatus and computer program product of data traffic signal estimating road traffic situation
CN107564290A (en) * 2017-10-13 2018-01-09 公安部交通管理科学研究所 A kind of urban road intersection saturation volume rate computational methods
CN110070734A (en) * 2019-05-14 2019-07-30 东南大学 Signalized intersections saturation headway estimation method based on gauss hybrid models

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
顾叶华: "环形信号交叉口饱和流率研究", 《科技视界》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113129594A (en) * 2021-04-23 2021-07-16 武汉科技大学 Intersection saturation flow rate calculation method based on bayonet data
CN113643531A (en) * 2021-07-20 2021-11-12 东北大学 Intersection lane saturation flow rate calculation method based on small time zone division statistics
CN113643531B (en) * 2021-07-20 2022-09-20 东北大学 Intersection lane saturation flow rate calculation method based on small time zone division statistics
CN113936453A (en) * 2021-09-09 2022-01-14 上海宝康电子控制工程有限公司 Information identification method and system based on headway

Also Published As

Publication number Publication date
CN111815973B (en) 2022-12-16

Similar Documents

Publication Publication Date Title
CN111815973B (en) Signalized intersection analysis method and related equipment
WO2021126243A1 (en) Systems and methods for detecting and responding to anomalous traffic conditions
CN111105628A (en) Parking lot portrait construction method and device
CN105448092A (en) Analysis method and apparatus of associated vehicles
CN115146865A (en) Task optimization method based on artificial intelligence and related equipment
CN112559831A (en) Link monitoring method and device, computer equipment and medium
CN108615359B (en) Vehicle foothold analysis method and device
CN110110339B (en) Japanese hydrologic forecast error correction method and system
CN114647684A (en) Traffic prediction method and device based on stacking algorithm and related equipment
CN114201482A (en) Dynamic population distribution statistical method and device, electronic equipment and readable storage medium
CN110867077A (en) Road optimization method based on traffic video structured data quality control
CN112652164B (en) Traffic time interval dividing method, device and equipment
CN110956826B (en) Method and device for generating traffic signal timing scheme and storage medium
CN111369790B (en) Vehicle passing record correction method, device, equipment and storage medium
CN111679959A (en) Computer performance data determination method and device, computer equipment and storage medium
CN115562934A (en) Service flow switching method based on artificial intelligence and related equipment
CN114627642A (en) Traffic jam identification method and device
CN114490262A (en) Database monitoring method, device, equipment and storage medium
CN114419889A (en) Time interval dividing method and device, electronic equipment and storage medium
CN114066288A (en) Intelligent data center-based emergency detection method and system for operation road
CN112185117B (en) Optimized evaluation method and device based on electric alarm data
CN115457766B (en) Method and system for predicting road congestion state
CN115022014B (en) Login risk identification method, device, equipment and storage medium
CN116630014A (en) Risk identification method, apparatus, device and storage medium
CN116821206A (en) Road experience speed determination method, device, vehicle, medium and equipment

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
GR01 Patent grant