EP1749288A1 - Method for determination of turning rates in a traffic network - Google Patents

Method for determination of turning rates in a traffic network

Info

Publication number
EP1749288A1
EP1749288A1 EP06763261A EP06763261A EP1749288A1 EP 1749288 A1 EP1749288 A1 EP 1749288A1 EP 06763261 A EP06763261 A EP 06763261A EP 06763261 A EP06763261 A EP 06763261A EP 1749288 A1 EP1749288 A1 EP 1749288A1
Authority
EP
European Patent Office
Prior art keywords
traffic
determined
subnetwork
exit
turning
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
EP06763261A
Other languages
German (de)
French (fr)
Other versions
EP1749288B1 (en
Inventor
Jürgen Mück
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.)
Siemens AG
Original Assignee
Siemens AG
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 Siemens AG filed Critical Siemens AG
Publication of EP1749288A1 publication Critical patent/EP1749288A1/en
Application granted granted Critical
Publication of EP1749288B1 publication Critical patent/EP1749288B1/en
Not-in-force legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G08—SIGNALLING
    • G08G—TRAFFIC CONTROL SYSTEMS
    • G08G1/00—Traffic control systems for road vehicles
    • G08G1/01—Detecting movement of traffic to be counted or controlled
    • G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions

Definitions

  • the invention relates to a method for determining turning rates in a road network according to the preamble of claim 1 and applications thereof in various road traffic engineering methods.
  • a basic task of traffic control systems in cities is the online determination of the traffic situation in the road network in order to inform about the traffic situation and to optimally control the connected subsystems.
  • systems are used for the large-area determination of the traffic situation, but also for the precise determination of the traffic status in subnetworks and for the optimization of associated traffic signal systems.
  • An essential task of these methods is the determination of the traffic flows in the road network, whereby the determination of the turn currents at junctions is a central algorithmic question.
  • This known method is characterized high demands on the data supply - such as an aggregation of measured data in intervals of two to three seconds - and a complex network modeling from. Partly also very special positions for the measuring cross sections are needed.
  • the model equation for the traffic intensities used as the basis for the estimation method either the same measurement interval for the left and right side is used as the time reference or the exit traffic volume for the measurement interval k is calculated from the access traffic intensities of the preceding measurement interval k-1.
  • This method suffers from the disadvantage that the estimation result depends strongly on the travel times between the measurement cross sections of the considered subnetwork.
  • the invention is therefore based on the object to provide a method of the type mentioned, which is robust in terms of travel times between the measuring cross sections and still works fast and accurate.
  • the object is achieved according to the invention by a generic method in which the access traffic volumes of a plurality of preceding measurement intervals are taken into account in the weighted sum for the exit traffic volume of a given measurement interval, with a forward-related turning rate to be determined being the sum of the corresponding turning rates gives the measurement intervals considered in the model equation.
  • the generalized time reference for the model equation has shown that the method according to the invention for determining the turning rates is robust with regard to travel times between the measuring cross sections and thus robust with respect to the size of the considered subnetworks. It is fast and has a previously unknown in practice accuracy. Finally, in contrast to previously used methods, the method according to the invention does not require any calibration.
  • At least one backward-related partial Network of the road network in which measuring cross sections are taken into account at a driveway and at exits of the subnet, set up a model equation in which the driveway traffic intensity is calculated as the weighted sum of the exit traffic volumes and the weighting factors correspond to the backward-related turning rates indicate in each case the proportion of an exit traffic volume which has flowed through the considered access, the turning rates being calculated from the model equation by means of a mathematical estimation method, wherein in the weighted sum for the inbound traffic strength of a given measurement interval the exit traffic strengths of a plurality of subsequent measurement intervals are taken into account, and wherein a backward-related turning rate to be determined results as the sum of the corresponding turning rates of the measuring intervals considered in the model equation.
  • the mathematical estimation method uses an extended, in particular non-linear Kalman filter, since it is a stochastic system with noise effects.
  • the stochastic parameters of the extended Kalman filter can be estimated in advance from the statistical analyzes of the data.
  • the filter is robust with regard to the parameterization and only requires the current measured values.
  • the non-linear Kalman filter is significantly more accurate, requiring less computation time and fewer data series than correlation analysis.
  • the proposed filter requires less calibration effort than heuristic methods of Operations Research.
  • the estimation method is interrupted when a traffic congestion is detected at a measuring cross section. This ensures the retention of the billing rates estimated before the onset of overload so that misjudgments be prevented due to the dammed vehicles. These are in fact a buffer that destroys the correlation between incoming and outgoing traffic flows.
  • forward-related and / or backward-related subnetworks are considered around the node where measurement cross sections in the approach and exit points of the node are taken into account, the turning rates being upwards be determined.
  • the turn rates can be estimated with advantage.
  • the turning rates for the entrances and exits of the subnetwork are determined according to the above-mentioned method, whereby measuring cross sections are considered only at the edge of the subnetwork but not in its interior, so that the source-destination traffic flows for this subnetwork are calculated from the determined turning rates and the recorded traffic volumes. In doing so, it is ensured that all relevant access measuring sections are included in the model equation for each one exit measuring section; Similarly, in a backward-related subnetwork to an access measuring cross section all relevant exit measuring cross sections are included. In this way, in subnetworks of limited size, a direct dynamic estimation of source-destination streams can be performed.
  • the number of measurement intervals considered is increased with increasing size of the considered subnetwork. If the measurement cross sections are close to each other, it is sufficient to consider a smaller number of preceding or subsequent measurement intervals. If, as the size of the considered subnetwork increases, the travel times between If the access and exit measurement cross sections grow, a larger number of measurement intervals must also be taken into account.
  • the measurement intervals considered are extended as the size of the considered subnetwork increases. Increasing the aggregation intervals to, for example, five minutes reduces noise disturbances in the estimation process.
  • a method for determining the traffic intensity at a road cross section of a road network are for a subnetwork of the road network, an entrance or exit the road cross section and its other entrances and / or exits measuring cross sections, according to the above-mentioned method provided determined turning rates and calculated from the provided bends and the detected at the measurement cross sections traffic levels of the other entrances and / or exits the traffic volume at the road cross section of an entrance or exit.
  • known turning rates for a subnetwork traffic volumes in its entrances or exits can be determined therefrom where no measured value is present.
  • the traffic intensity determined for the roadway cross section is used as a substitute value for a faulty or failed measuring cross section.
  • the traffic intensity is detected at a measuring cross section and at the second end point no measuring cross section is arranged, the traffic intensity at the second end point according to the above-described method and from this, the number of vehicles in the carriageway section is determined by temporal integration of the difference between the traffic volume flowing into the roadway section and the traffic volume flowing therefrom.
  • a homogeneous system of equations is first set up from the vehicle preservation of the actual forward and reverse-related turning rates for the correction factors to be determined, then off the optimization of the homogeneous system of equations, together with a constraint which excludes the trivial solution, whereby the correction factors result as a solution to the optimization problem.
  • the determined correction factors are divided by their median value. It is assumed that in a subnet less than half of all measurement cross sections count too many vehicles and fewer than half count too few vehicles, so that the median value of the list of specific correction factors can be used as the reference value. Due to the mentioned correction division this is then at the value one.
  • a considered subnetwork of the road network is divided into island networks, each island network comprising only measuring cross sections at its network edge, and correction factors for the island grids are determined.
  • a parameter is calculated in the solution of the optimization problem whose value is used as a measure of the estimated quality of the turning rates.
  • this parameter is ideally close to zero, if the turning rates are accurately estimated and between the measuring cross sections neither vehicle losses nor vehicle increases occur and all measurement errors are proportional.
  • an error message is output if the value of the parameter exceeds a predefinable limit.
  • This is an indicator of inaccurately estimated turning rates, unrecorded access or exit traffic levels, or measurement errors of a non-proportional nature, which may arise, for example, if not all relevant access and exit traffic levels in the subnetwork are measured.
  • two correction factors are determined for each measuring cross section, which is shared by two adjoining island networks, the correction factors of each island network being scaled such that the correction factors of common measuring cross sections are equal to one another.
  • the ratio of the correction factors within each island network remains unchanged, but common estimation errors between the island networks are adjusted.
  • the traffic intensities recorded at the measuring cross sections and the turning rates determined by estimation are calibrated by means of the determined correction factors.
  • FIG. 1 subnetworks of a network section with a node
  • FIG. 2 the forward-related subnetwork of FIG. 1 with turn-off connections
  • FIG. 1 subnetworks of a network section with a node
  • FIG. 3 shows the backward-related subnetwork from FIG. 1 with turn-related connections, FIG.
  • FIG. 4 shows an island network around a node
  • FIG. 5 shows the time course of estimated turning rates of the island network according to FIG. 4
  • FIG. 1 shows a network section of a road network having a node, for example a city road network in which the turning rates of the traffic flows are to be determined for the purpose of traffic control.
  • the roadway from driveway 1 to exit 3 does not include any measuring cross-sections in the exemplary embodiment described. At all other entrances and exits of the junction are measuring cross sections with detectors for detecting entrance traffic volumes q "(n) and exit traffic volumes q ° ut (n) in predetermined measuring intervals n.
  • Basic element of the inventive method for dynamic estimation of turning rates 1 shows a first subnetwork fw whose network edge is shown as a dashed-dotted line and which measuring cross-sections in exit 1 as well as in the relevant access roads 3 and 4.
  • the access roads 3 and 4 are relevant, because on them partial traffic flows into the subnetwork fw, which flow through the exit 1 from the subnetwork fw.
  • FIG. 2 shows the subnetwork fw from FIG. 1 with the associated turn-off relationships.
  • the turning rate m ⁇ "(k) indicates the proportion of traffic volume q j n (k) measured on access 3, which flows out of subnetwork fw through exit 1 and thus contributes to the measured traffic q o ut (k)
  • a - naloges applies to the turning rate mf "(k) in relation to the entrance traffic volume q ⁇ fk).
  • the subnetwork fw (forward) models the turn-related relationships forward in time.
  • the exit traffic volume q ° ut (k) can be modeled as a weighted sum of the entrance traffic volumes q j n (k) and q ⁇ fk) for a given measurement interval k, the weighting factors corresponding to the corresponding turning rate m ⁇ "(k) and mf" (k).
  • r can be used at r for a j exit accesses relevant:
  • the turning rates m i "(k) can be determined by means of a mathematical estimation method.
  • the estimation method used here a generalized time reference.
  • n k-1 , k-2, ...
  • the inventive method is robust in terms of travel times between the measuring cross sections, with high accuracy and sufficient speed.
  • FIG. 1 shows such a subnetwork bw whose network edge is shown as a dashed line and which comprises measuring cross sections in the access road 3 and in the relevant exits 1, 2 and 4.
  • the exits 1, 2 and 4 are relevant because flow on them partial traffic flows from the sub-network bw, which have flowed through the driveway 3 in the sub-network bw.
  • FIG. 3 shows this subnetwork bw with the associated turn-off relationships.
  • the turning rate m ⁇ fk) indicates the proportion of traffic volume q ⁇ Ck measured in exit 2, which has flowed into the subnetwork bw via driveway 3 and thus contributes to the traffic strength q j n (k) measured there.
  • the subnetwork bw thus models the turn-related relationships backwards in time. In general, at s for departures j relevant exits i taking into account z measuring intervals n again the weighted sum can be applied:
  • ⁇ M ⁇ ⁇ ⁇ k + i - D • q ° ut ( k + 1-D
  • an extended Kalman filter for estimating the turning rates is used as the mathematical estimation method.
  • the measurement cross sections are close to each other, for example one to two intervening traffic signal systems, a smaller number of measurement intervals to be taken into account is sufficient, according to experience three or four.
  • the model equation can also be applied in larger subnetworks whose detectors are evaluated only at the edge of the network and not inside the network. It merely has to be ensured that in the forward-related case to an exit measuring cross-section of a sub-network all relevant access measuring cross-sections are included, or that in the backward-related case to an access measuring cross section all relevant exit measuring cross sections are included.
  • the estimation process is interrupted if an overload is detected on one of the measuring cross sections used for a subnetwork.
  • an accounting approach is used to determine traffic jams in access roads to traffic lights. This determines the number of vehicles, ie the traffic jam length, in one of the driveways by integrating the traffic volumes at the end points of the roadway section over time. If, as is usually the case in practice, only one measuring cross-section is located at one of the two end points, the traffic volume at the other end point can be estimated via turning rates determined according to the invention.
  • a further advantage of an application of the method according to the invention is described below for the case where a measuring cross section is located at the entrance to a roadway section, while the traffic volume at the exit from the roadway section is estimated via turning rates, since there is no measuring cross section at this roadway cross section.
  • the above model equations are characterized by the fact that the turn rates of the entrances and exits are estimated in relation to the exit or access traffic volume. If the measured values of the traffic intensities have proportional errors, these are compensated in the turn rates. The calculated access traffic volumes are consistent with the measured exit traffic. In this way, the quality of balances is significantly increased without the need for special calibrations for each balanced road section.
  • q out M fw • q in
  • q in and q out are column vectors whose components represent the access traffic volumes q ⁇ n and exit traffic volumes q ° ut of all measurement cross sections in the entrances and exits of the node arms i, while M fw represents a (nxn)
  • Matrix means whose elements are the turn rates mf ".
  • the column vectors q in and q out comprise in this case all the measuring cross sections in the same sequence, even if there is no driving relationship between components q ⁇ the right side and q, the left side of the equation.
  • ItI 1 - 0.
  • detectors of the measuring cross-sections capture the actual traffic volumes only with a certain accuracy. This may result from their location, e.g. Passing through vehicles of two lanes, but also through a faulty internal calibration process, in which the measured values drift, as it often occurs in practice over time.
  • n in access and out n exit measurement cross sections results in an over-determined homogeneous system of equations for the U 1 with n in + n out ⁇ n equations.
  • the simplest case is that of a single traffic signal system in which all relevant entrances and exits are recorded.
  • Wi and W2 represent selectable weights.
  • the formulated optimization problem can be formally applied to a network as a whole.
  • this has two disadvantages: On the one hand, the computational effort for the optimization increases disproportionately with the number of measurement cross sections. On the other hand, the fact that many measuring cross sections in the network have further measuring cross sections in both directions leads to leveling effects of the estimated correction factor
  • the method can be refined by properly decomposing the network.
  • the application to subnetworks has proven to be favorable, which are defined such that from any point within such a network only those measuring points are part of this subnetwork which can be reached directly via network edges.
  • Each such subnet constitutes, as it were, an island network with entrances and exits into the It summarizes all forward or backward related subnets that have the same entry and exit measurement cross sections.
  • FIG. 4 shows such a stand-alone grid, which generally arises around the traffic signal system at the usual detector equipment at nodes. It has measuring cross-sections only at the network edge, which is illustrated in FIG. 4 with a double-dashed line.
  • the parameter P resulting from the solution to the optimization problem for an island grid provides further information about possible detector disturbances: the value P of the solution to the optimization problem for an island grid ideally approaches zero if the turning rates are exactly estimated and between the measuring cross sections no vehicle losses (sinks) or increases (sources) occur, and all measurement errors are proportional. If P results in a value significantly greater than zero, e.g. Two, this is an indication of inaccurately estimated turning rates, missed entrances and exits, or measurement errors other than proportional. Typically, such an error arises when not all traffic flows in the relevant entrances and exits of such a subnetwork are detected or e.g. due to incorrect supply / wiring, measuring sections were assigned incorrectly to the sections.
  • substitute values if detectors or entire measuring sections have failed. Such failures are either already detected in the hardware of the devices and reported further or they can be detected via simple plausibility checks. In this case, substitute values can simply be formed based on the measured and calibrated traffic intensities and turning rates of the surrounding measuring cross sections. This represents a considerable leap in quality over known methods, which in the simple case simply replace missing measured values with previously-dependent time-dependent values or-very costly-determine these via pre-recorded, day-type-specific hydrographs and the current time.
  • the turn rates estimated by the described method are accurate enough for both traffic estimation and reuse in adaptive network control techniques and far surpassing the estimation quality of conventional metering based methods.
  • the estimation of the correction factors U 1 for all measuring cross sections on the one hand enables an online correction of the turning rates and traffic volume counts.
  • the thus possible method for network-related replacement value formation works without hydrographs or pre-supplied default count values and also works with several closely spaced detector failures, since in this case substitutions can also be made recursively.
  • the consistency check provides valuable information for the maintenance of the detector network.
  • Corresponding notification mechanisms allow a city's maintenance service to efficiently and quickly repair defective detectors, with significant cost savings.

Landscapes

  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Traffic Control Systems (AREA)
  • Data Exchanges In Wide-Area Networks (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Road Repair (AREA)

Abstract

Methods for determining turning rates in a road network are provided. Traffic volumes are recorded at measurement cross-sections at predefinable measuring intervals. For at least one forward-related subnetwork of the road network in which measurement cross-sections are taken into account at an exit and at entries of the subnetwork, a model equation is formulated in which the exit traffic volume is set as the weighted sum of the entry traffic volumes and the weighting factors correspond to the forward-related turning rates which specify in each case the portion of an entry traffic volume flowing out through the exit taken into account, and wherein the forward-related turning rates are calculated on the basis of the model equation using a mathematical estimation method.

Description

Beschreibungdescription
Verfahren zur Ermittlung von Abbiegeraten in einem StraßennetzMethod for determining turning rates in a road network
Die Erfindung betrifft ein Verfahren zur Ermittlung von Abbiegeraten in einem Straßennetz gemäß dem Oberbegriff des Patentanspruches 1 sowie Anwendungen dessen in verschiedenen straßenverkehrstechnischen Verfahren .The invention relates to a method for determining turning rates in a road network according to the preamble of claim 1 and applications thereof in various road traffic engineering methods.
Eine Basisaufgabe von Verkehrssteuerungssystemen in Städten ist die Online-Ermittlung der Verkehrslage im Straßennetz, um über die Verkehrssituation zu informieren und die angeschlossenen Subsysteme optimal zu steuern. Hierbei kommen Systeme zur großflächigen Bestimmung der Verkehrslage, aber auch zur genauen Ermittlung des Verkehrszustandes in Teilnetzen und zur Optimierung zugeordneter Lichtsignalanlagen zum Einsatz. Eine wesentliche Aufgabe dieser Verfahren ist die Bestimmung der Verkehrsströme im Straßennetz, wobei die Ermittlung der Abbiegeströme an Knotenpunkten eine zentrale algorithmische Fragestellung ist.A basic task of traffic control systems in cities is the online determination of the traffic situation in the road network in order to inform about the traffic situation and to optimally control the connected subsystems. In this case, systems are used for the large-area determination of the traffic situation, but also for the precise determination of the traffic status in subnetworks and for the optimization of associated traffic signal systems. An essential task of these methods is the determination of the traffic flows in the road network, whereby the determination of the turn currents at junctions is a central algorithmic question.
Aus der Dissertation "Ein Verfahren zur gekoppelten Schätzung von Kantenbelastungen, Abbiegequoten und Störungen in Stadtstraßennetzen", veröffentlicht in der Schriftenreihe der Arbeitsgruppe Automatisierungstechnik, Technische Universität Hamburg-Harburg, Heft 20, Mai 2001, ist ein Verfahren der eingangs genannten Art bekannt. Der Schätzer wird im Zeittakt von zwei Sekunden mit sämtlichen gemessenen knotenzu- und - wegfließenden Verkehrsstärken versorgt. Zur Abbiegequotenschätzung kommt eine rein dynamische Methode zum Einsatz, die auch der Klasse der rekursiven Verfahren zuzurechnen ist. Durch eine zeitlich hoch aufgelöste Berechnungsweise ergibt sich die Unterteilung in phasengruppenorientierte Subsysteme automatisch. Bei der Berechnung der Abbiegequoten werden die Veränderungen der Ströme von Aus- und Einfahrt in einem Zeitintervall k gegenüber dem vorangehenden Zeitintervall k-1 berücksichtigt. Dieses bekannte Verfahren zeichnet sich durch hohe Anforderungen an die Datenbereitstellung - etwa eine Aggregation von Messdaten in Intervallen von zwei bis drei Sekunden - und eine aufwendige Netzmodellierung aus. Teilweise werden auch ganz spezielle Positionen für die Messquerschnitte benötigt. Bei der als Grundlage für das Schätzverfahren herangezogenen Modellgleichung für die Verkehrsstärken wird als Zeitbezug entweder das gleiche Messintervall für die linke und rechte Seite verwendet oder die Ausfahrt- Verkehrsstärke zum Messintervall k aus den Zufahrt- Verkehrsstärken des vorangegangenen Messintervalls k-1 berechnet. Dieses Verfahren leidet unter dem Nachteil, dass das Schätzergebnis stark von den Reisezeiten zwischen den Messquerschnitten des betrachteten Teilnetzes abhängt.From the dissertation "A method for the coupled estimation of edge loads, turn rates and disturbances in city road networks", published in the series of the working group automation technology, Hamburg University of Technology, Issue 20, May 2001, a method of the type mentioned is known. The estimator is supplied with all measured node-to-and-out traffic volumes at a time interval of two seconds. For the turn rate estimation, a purely dynamic method is used, which is also attributable to the class of recursive methods. By a high-resolution calculation method, the subdivision into phase-group-oriented subsystems automatically results. When calculating the turn rates, the changes in the flows of exit and exit are taken into account in a time interval k compared to the preceding time interval k-1. This known method is characterized high demands on the data supply - such as an aggregation of measured data in intervals of two to three seconds - and a complex network modeling from. Partly also very special positions for the measuring cross sections are needed. In the model equation for the traffic intensities used as the basis for the estimation method, either the same measurement interval for the left and right side is used as the time reference or the exit traffic volume for the measurement interval k is calculated from the access traffic intensities of the preceding measurement interval k-1. This method suffers from the disadvantage that the estimation result depends strongly on the travel times between the measurement cross sections of the considered subnetwork.
Der Erfindung liegt daher die Aufgabe zugrunde, ein Verfahren der eingangs genannten Art bereitzustellen, welches robust hinsichtlich der Reisezeiten zwischen den Messquerschnitten ist und trotzdem schnell und genau arbeitet.The invention is therefore based on the object to provide a method of the type mentioned, which is robust in terms of travel times between the measuring cross sections and still works fast and accurate.
Die Aufgabe wird erfindungsgemäß gelöst durch ein gattungsgemäßes Verfahren, bei dem in der gewichteten Summe für die Ausfahrt-Verkehrsstärke eines gegebenen Messintervalls die Zufahrt-Verkehrsstärken einer Mehrzahl vorangehender Messintervalle berücksichtigt werden, wobei sich eine zu ermittelnde vorwärts-bezogene Abbiegerate als Summe aus den entsprechenden Abbiegeraten der in der Modellgleichung berücksichtigten Messintervalle ergibt. Durch den verallgemeinerten Zeitbezug für die Modellgleichung hat sich gezeigt, dass das erfindungsgemäße Verfahren zur Ermittlung der Abbiegeraten robust hinsichtlich der Reisezeiten zwischen den Messquerschnitten und damit robust hinsichtlich der Größe der betrachteten Teilnetze ist. Es ist schnell und weist eine bisher in der Praxis nicht bekannte Genauigkeit auf. Schließlich erfordert das erfindungsgemäße Verfahren im Gegensatz zu bisher eingesetzten Verfahren keine Kalibrierung.The object is achieved according to the invention by a generic method in which the access traffic volumes of a plurality of preceding measurement intervals are taken into account in the weighted sum for the exit traffic volume of a given measurement interval, with a forward-related turning rate to be determined being the sum of the corresponding turning rates gives the measurement intervals considered in the model equation. The generalized time reference for the model equation has shown that the method according to the invention for determining the turning rates is robust with regard to travel times between the measuring cross sections and thus robust with respect to the size of the considered subnetworks. It is fast and has a previously unknown in practice accuracy. Finally, in contrast to previously used methods, the method according to the invention does not require any calibration.
In einer vorteilhaften Ausgestaltung des erfindungsgemäßen Verfahrens wird für wenigstens ein rückwärts-bezogenes Teil- netz des Straßennetzes, bei welchen Messquerschnitte an einer Zufahrt und an Ausfahrten des Teilnetzes berücksichtigt werden, eine Modellgleichung aufgestellt, in der die Zufahrt- Verkehrsstärke als gewichtete Summe der Ausfahrt-Verkehrsstärken angesetzt wird und die Gewichtungsfaktoren den rück- wärts-bezogenen Abbiegeraten entsprechen, welche jeweils den Anteil einer Ausfahrt-Verkehrsstärke angeben, der durch die berücksichtigte Zufahrt eingeflossen ist, wobei ausgehend von der Modellgleichung mittels eines mathematischen Schätzverfahrens die Abbiegeraten berechnet werden, wobei in der ge- wichteten Summe für die Zufahrtverkehrsstärke eines gegebenen Messintervalls die Ausfahrtverkehrsstärken einer Mehrzahl nachfolgender Messintervalle berücksichtigt werden, und wobei sich eine zu ermittelnde rückwärts-bezogene Abbiegerate als Summe aus den entsprechenden Abbiegeraten der in der Modellgleichung berücksichtigten Messintervalle ergibt. Durch die Schätzung sowohl zeitlich vorwärts- als auch rückwärtsbezogener Abbiegeraten gewinnt das erfindungsgemäße Verfahren noch weiter an Robustheit und Genauigkeit.In an advantageous embodiment of the method according to the invention, at least one backward-related partial Network of the road network, in which measuring cross sections are taken into account at a driveway and at exits of the subnet, set up a model equation in which the driveway traffic intensity is calculated as the weighted sum of the exit traffic volumes and the weighting factors correspond to the backward-related turning rates indicate in each case the proportion of an exit traffic volume which has flowed through the considered access, the turning rates being calculated from the model equation by means of a mathematical estimation method, wherein in the weighted sum for the inbound traffic strength of a given measurement interval the exit traffic strengths of a plurality of subsequent measurement intervals are taken into account, and wherein a backward-related turning rate to be determined results as the sum of the corresponding turning rates of the measuring intervals considered in the model equation. By estimating both forward and backward related turn rates, the inventive method gains even more robustness and accuracy.
Vorzugsweise wird bei dem mathematischen Schätzverfahren ein erweitertes, insbesondere nichtlineares Kaiman-Filter angewandt, da es sich um ein stochastisches System mit Rauscheffekten handelt. Die stochastischen Parameter des erweiterten Kaiman-Filters können aus den statistischen Analysen der Daten vorab abgeschätzt werden. Gleichzeitig ist das Filter robust hinsichtlich der Parametrierung und benötigt nur die aktuellen Messwerte. Das nichtlineare Kaiman-Filter ist erheblich genauer, benötigt weniger Rechenzeit und weniger Datenreihen als Korrelationsanalysen. Daneben erfordert das vorgeschlagene Filter einen geringeren Kalibrieraufwand als heuristische Methoden des Operations Research.Preferably, the mathematical estimation method uses an extended, in particular non-linear Kalman filter, since it is a stochastic system with noise effects. The stochastic parameters of the extended Kalman filter can be estimated in advance from the statistical analyzes of the data. At the same time, the filter is robust with regard to the parameterization and only requires the current measured values. The non-linear Kalman filter is significantly more accurate, requiring less computation time and fewer data series than correlation analysis. In addition, the proposed filter requires less calibration effort than heuristic methods of Operations Research.
In einer besonderen Ausführungsform des erfindungsgemäßen Verfahrens wird das Schätzverfahren unterbrochen, wenn an einem Messquerschnitt eine Verkehrsüberlastung erfasst wird. Dies sichert die Beibehaltung der vor Auftreten der Überlastung geschätzten Abbiegeraten, so dass Fehleinschätzungen aufgrund der aufgestauten Fahrzeuge verhindert werden. Diese stellen nämlich einen Puffer dar, der die Korrelation zwischen zu- und abfließenden Verkehrsströmen zerstört.In a particular embodiment of the method according to the invention, the estimation method is interrupted when a traffic congestion is detected at a measuring cross section. This ensures the retention of the billing rates estimated before the onset of overload so that misjudgments be prevented due to the dammed vehicles. These are in fact a buffer that destroys the correlation between incoming and outgoing traffic flows.
In einer bevorzugten Anwendung der Erfindung für ein Verfahren zur Ermittlung von Abbiegeraten an einem Knotenpunkt des Straßennetzes werden vorwärts-bezogene und/oder rückwärtsbezogene Teilnetze um den Knotenpunkt betrachtet, bei welchen Messquerschnitte in den Zufahrten und Ausfahrten des Knotenpunktes berücksichtigt werden, wobei die Abbiegeraten nach oben beschriebenem Verfahren ermittelt werden. Durch eine geeignete Wahl der Teilnetze um eine - gegebenenfalls Lichtsignal-gesteuerte - Kreuzung lassen sich die Abbiegequoten mit Vorteil schätzen.In a preferred application of the invention for a method of determining turning rates at a node of the road network, forward-related and / or backward-related subnetworks are considered around the node where measurement cross sections in the approach and exit points of the node are taken into account, the turning rates being upwards be determined. By a suitable choice of subnets to a - optionally light signal-controlled - intersection, the turn rates can be estimated with advantage.
In einer weiteren Anwendung der Erfindung bei einem Verfahren zur Ermittlung von Quelle-Ziel-Verkehrsströmen eines Teilnetzes werden für die Zufahrten und Ausfahrten des Teilnetzes die Abbiegeraten gemäß oben genanntem Verfahren ermittelt, wobei Messquerschnitte nur am Rand des Teilnetzes aber nicht in seinem Inneren berücksichtigt werden, so dass aus den ermittelten Abbiegeraten und den erfassten Verkehrsstärken die Quelle-Ziel-Verkehrsströme für dieses Teilnetz berechnet werden. Dabei wird sichergestellt, dass zu je einem Ausfahrt- Messquerschnitt alle relevanten Zufahrt-Messquerschnitte in die Modellgleichung einbezogen werden; analog sind bei einem rückwärts-bezogenen Teilnetz zu einem Zufahrt-Messquerschnitt alle relevanten Ausfahrt-Messquerschnitte einzubeziehen. Auf diese Weise kann in Teilnetzen begrenzter Größe eine direkte dynamische Schätzung von Herkunfts-Ziel-Strömen durchgeführt werden.In a further application of the invention in a method for determining source-destination traffic flows of a subnetwork, the turning rates for the entrances and exits of the subnetwork are determined according to the above-mentioned method, whereby measuring cross sections are considered only at the edge of the subnetwork but not in its interior, so that the source-destination traffic flows for this subnetwork are calculated from the determined turning rates and the recorded traffic volumes. In doing so, it is ensured that all relevant access measuring sections are included in the model equation for each one exit measuring section; Similarly, in a backward-related subnetwork to an access measuring cross section all relevant exit measuring cross sections are included. In this way, in subnetworks of limited size, a direct dynamic estimation of source-destination streams can be performed.
Vorzugsweise wird die Anzahl der berücksichtigten Messintervalle mit zunehmender Größe des betrachteten Teilnetzes vergrößert. Falls die Messquerschnitte nahe beieinander liegen, genügt es, eine kleinere Anzahl von vorangehenden bzw. nachfolgenden Messintervallen zu berücksichtigen. Wenn mit zunehmender Größe des betrachteten Teilnetzes die Reisezeiten zwi- sehen den Zufahrt- und Ausfahrt-Messquerschnitten wachsen, muss auch eine größere Anzahl von Messintervallen berücksichtigt werden.Preferably, the number of measurement intervals considered is increased with increasing size of the considered subnetwork. If the measurement cross sections are close to each other, it is sufficient to consider a smaller number of preceding or subsequent measurement intervals. If, as the size of the considered subnetwork increases, the travel times between If the access and exit measurement cross sections grow, a larger number of measurement intervals must also be taken into account.
Bevorzugt werden die berücksichtigten Messintervalle mit zunehmender Größe des betrachteten Teilnetzes verlängert. Eine Vergrößerung der Aggregationsintervalle auf beispielsweise fünf Minuten reduziert die Störungen des Schätzprozesses durch Rauschen.Preferably, the measurement intervals considered are extended as the size of the considered subnetwork increases. Increasing the aggregation intervals to, for example, five minutes reduces noise disturbances in the estimation process.
In einer anderen vorteilhaften Anwendung der Erfindung auf ein Verfahren zur Ermittlung der Verkehrsstärke an einem Fahrbahnquerschnitt eines Straßennetzes werden für ein Teilnetz des Straßennetzes, dessen eine Zu- oder Ausfahrt den Fahrbahnquerschnitt und dessen andere Zu- und/oder Ausfahrten Messquerschnitte aufweisen, gemäß oben genanntem Verfahren ermittelte Abbiegeraten bereitgestellt und aus den bereitgestellten Abbiegeraten und den an den Messquerschnitten er- fassten Verkehrsstärken der anderen Zu- und/oder Ausfahrten die Verkehrsstärke am Fahrbahnquerschnitt der einen Zu- bzw. Ausfahrt berechnet. Bei bekannten Abbiegeraten für ein Teilnetz lassen sich daraus Verkehrsstärken in dessen Zu- oder Ausfahrten bestimmen, wo kein Messwert vorliegt.In another advantageous application of the invention to a method for determining the traffic intensity at a road cross section of a road network are for a subnetwork of the road network, an entrance or exit the road cross section and its other entrances and / or exits measuring cross sections, according to the above-mentioned method provided determined turning rates and calculated from the provided bends and the detected at the measurement cross sections traffic levels of the other entrances and / or exits the traffic volume at the road cross section of an entrance or exit. With known turning rates for a subnetwork, traffic volumes in its entrances or exits can be determined therefrom where no measured value is present.
Vorzugsweise wird die für den Fahrbahnquerschnitt ermittelte Verkehrsstärke als Ersatzwert für einen fehlerhaften oder ausgefallenen Messquerschnitt herangezogen.Preferably, the traffic intensity determined for the roadway cross section is used as a substitute value for a faulty or failed measuring cross section.
In einer anderen bevorzugten Anwendung der Erfindung auf ein Verfahren zur Ermittlung der Fahrzeuganzahl innerhalb eines Fahrbahnabschnitts, an dessen ersten Endpunkt die Verkehrsstärke an einem Messquerschnitt erfasst wird und an dessen zweiten Endpunkt kein Messquerschnitt angeordnet ist, wird die Verkehrsstärke am zweiten Endpunkt gemäß vorstehend beschriebenem Verfahren und daraus die Fahrzeuganzahl im Fahrbahnabschnitt durch zeitliche Integration der Differenz aus der in den Fahrbahnabschnitt zufließenden Verkehrsstärke und der aus diesem abfließenden Verkehrsstärke ermittelt. Durch diesen Bilanzierungsansatz können beispielsweise Staulängen in Zufahrten von Lichtsignal-gesteuerten Knotenpunkten ermittelt werden, auch wenn nur ein Messquerschnitt an einem der beiden Endpunkte des betrachteten Fahrbahnabschnittes in der Zufahrt vorhanden ist.In another preferred application of the invention to a method for determining the number of vehicles within a lane section, at the first end point, the traffic intensity is detected at a measuring cross section and at the second end point no measuring cross section is arranged, the traffic intensity at the second end point according to the above-described method and from this, the number of vehicles in the carriageway section is determined by temporal integration of the difference between the traffic volume flowing into the roadway section and the traffic volume flowing therefrom. By This balancing approach, for example, jam lengths can be determined in access roads of light signal-controlled nodes, even if only a measuring cross-section is present at one of the two endpoints of the considered lane section in the driveway.
In ebenfalls vorteilhafter Anwendung der Erfindung auf ein Verfahren zur Bestimmung von Korrekturfaktoren für Abbiegeraten, welche nach oben beschriebenem Verfahren ermittelt werden, wird zunächst aus der Fahrzeugerhaltung der tatsächlichen vorwärts- und rückwärts-bezogenen Abbiegeraten für die zu bestimmenden Korrekturfaktoren ein homogenes Gleichungssystem aufgestellt, dann aus dem homogenen Gleichungssystem zusammen mit einer die triviale Lösung ausschließenden Nebenbedingung ein Optimierungsproblem gewonnen, wobei die Korrekturfaktoren sich als Lösung des Optimierungsproblems ergeben. Hierdurch lassen sich beispielsweise konstante prozentuale Fehler bei der Erfassung der Verkehrsstärken kompensieren, die sich aus der speziellen Lage, aber auch durch einen defekten internen Kalibriervorgang von Detektoren ergeben können.In a likewise advantageous application of the invention to a method for determining correction factors for turning rates, which are determined by the method described above, a homogeneous system of equations is first set up from the vehicle preservation of the actual forward and reverse-related turning rates for the correction factors to be determined, then off the optimization of the homogeneous system of equations, together with a constraint which excludes the trivial solution, whereby the correction factors result as a solution to the optimization problem. This makes it possible, for example, to compensate for constant percentage errors in the detection of the traffic volumes, which can result from the special position, but also from a defective internal calibration process of detectors.
Vorzugsweise werden die bestimmten Korrekturfaktoren durch ihren Medianwert dividiert. Dabei geht man davon aus, dass in einem Teilnetz weniger als die Hälfte aller Messquerschnitte zu viele Fahrzeuge und weniger als die Hälfte zu wenige Fahrzeuge zählen, so dass als Referenzwert der Medianwert der Liste der bestimmten Korrekturfaktoren verwendet werden kann. Durch die genannte Korrekturdivision liegt dieser dann beim Wert eins .Preferably, the determined correction factors are divided by their median value. It is assumed that in a subnet less than half of all measurement cross sections count too many vehicles and fewer than half count too few vehicles, so that the median value of the list of specific correction factors can be used as the reference value. Due to the mentioned correction division this is then at the value one.
In einer bevorzugten Ausgestaltung dieser Anwendung wird ein betrachtetes Teilnetz des Straßennetzes in Inselnetze zerlegt, wobei jedes Inselnetz nur Messquerschnitte an seinem Netzrand umfasst, und es werden Korrekturfaktoren für die Inselnetze bestimmt. Durch diese geeignete Zerlegung in Inselnetze wird einerseits der Rechenaufwand für die Optimierung reduziert und andererseits Nivellierungseffekte der geschätz- ten Korrekturfaktoren vermieden, die in Netzen auftreten, bei denen viele Messstellen in beide Richtungen weitere Messquerschnitte haben.In a preferred embodiment of this application, a considered subnetwork of the road network is divided into island networks, each island network comprising only measuring cross sections at its network edge, and correction factors for the island grids are determined. By means of this suitable decomposition in island grids, on the one hand, the computational outlay for the optimization is reduced, and on the other hand, the leveling effects of the estimated correction factors that occur in networks in which many measuring points have further measuring cross sections in both directions.
Mit Vorteil wird hierbei überprüft, ob die bestimmten Korrekturfaktoren innerhalb eines vorgegebenen Wertebereiches liegen. Bei Verlassen des vorgebbaren Wertebereiches liegt eine derart große Abweichung zwischen geschätzten und gemessenen Größen vor, dass auf diese Weise eine Fehlermeldung ausgegeben werden kann.It is advantageously checked here whether the determined correction factors lie within a predetermined value range. When leaving the predefinable value range, there is such a large deviation between estimated and measured variables that an error message can be output in this way.
In vorteilhafter Ausgestaltung dieser Anwendung wird bei der Lösung des Optimierungsproblems eine Kenngröße berechnet, deren Wert als Maß für die Schätzqualität der Abbiegeraten verwendet wird. Für ein Inselnetz geht diese Kenngröße im Idealfall gegen Null, wenn die Abbiegeraten exakt geschätzt sind und zwischen den Messquerschnitten weder Fahrzeugverluste noch Fahrzeugzunahmen auftreten und alle Messfehler proportionaler Art sind.In an advantageous embodiment of this application, a parameter is calculated in the solution of the optimization problem whose value is used as a measure of the estimated quality of the turning rates. For an island grid, this parameter is ideally close to zero, if the turning rates are accurately estimated and between the measuring cross sections neither vehicle losses nor vehicle increases occur and all measurement errors are proportional.
Vorzugsweise wird eine Fehlermeldung ausgegeben, wenn der Wert der Kenngröße eine vorgebbare Schranke übersteigt. Dies ist ein Indikator für ungenau geschätzte Abbiegeraten, nicht erfasste Zufahrt- oder Ausfahrt-Verkehrsstärken oder Messfehler nicht proportionaler Art, die beispielsweise entstehen, wenn nicht alle relevanten Zufahrt- und Ausfahrt-Verkehrsstärken im Teilnetz gemessen werden.Preferably, an error message is output if the value of the parameter exceeds a predefinable limit. This is an indicator of inaccurately estimated turning rates, unrecorded access or exit traffic levels, or measurement errors of a non-proportional nature, which may arise, for example, if not all relevant access and exit traffic levels in the subnetwork are measured.
In einer weiteren bevorzugten Ausgestaltung der beschriebenen Anwendung der Erfindung werden für jeden Messquerschnitt, den sich zwei aneinander grenzende Inselnetze teilen, zwei Korrekturfaktoren bestimmt, wobei die Korrekturfaktoren je eines Inselnetzes derart skaliert werden, dass sich die Korrekturfaktoren gemeinsamer Messquerschnitte einander angleichen. Durch diesen weiteren Optimierungsschritt bleibt das Verhältnis der Korrekturfaktoren innerhalb jedes Inselnetzes unverändert, wobei aber gemeinsame Schätzfehler zwischen den Inselnetzen abgeglichen werden. Vorzugsweise werden die an den Messquerschnitten erfassten Verkehrsstärken und die durch Schätzung ermittelten Abbiegeraten mittels der bestimmten Korrekturfaktoren kalibriert.In a further preferred embodiment of the described application of the invention, two correction factors are determined for each measuring cross section, which is shared by two adjoining island networks, the correction factors of each island network being scaled such that the correction factors of common measuring cross sections are equal to one another. By means of this further optimization step, the ratio of the correction factors within each island network remains unchanged, but common estimation errors between the island networks are adjusted. Preferably, the traffic intensities recorded at the measuring cross sections and the turning rates determined by estimation are calibrated by means of the determined correction factors.
Weitere Vorteile des erfindungsgemäßen Verfahrens sowie seiner bevorzugten Anwendungen ergeben sich aus einem konkreten Ausführungsbeispiel, welches nachfolgend anhand der Zeichnungen näher beschrieben wird, in derenFurther advantages of the method according to the invention and its preferred applications will become apparent from a concrete embodiment, which is described below with reference to the drawings, in whose
FIG 1 Teilnetze eines Netzausschnittes mit Knotenpunkt, FIG 2 das vorwärts-bezogene Teilnetz aus FIG 1 mit Abbiegebeziehungen,1 subnetworks of a network section with a node, FIG. 2 the forward-related subnetwork of FIG. 1 with turn-off connections, FIG.
FIG 3 das rückwärts-bezogene Teilnetz aus FIG 1 mit Abbiegebeziehungen,3 shows the backward-related subnetwork from FIG. 1 with turn-related connections, FIG.
FIG 4 ein Inselnetz um einen Knotenpunkt und FIG 5 den zeitlichen Verlauf erfindungsgemäß geschätzter Abbiegeraten des Inselnetzes aus FIG 44 shows an island network around a node, and FIG. 5 shows the time course of estimated turning rates of the island network according to FIG. 4
schematisch veranschaulicht sind.are illustrated schematically.
FIG 1 zeigt einen einen Knotenpunkt aufweisenden Netzausschnitt eines Straßennetzes, etwa eines Stadtstraßennetzes, in welchem die Abbiegeraten der Verkehrsströme zum Zwecke der Verkehrssteuerung ermittelt werden sollen. Der Knotenpunkt weist vier Knotenarme i (i = 1, ..., 4) mit Zu- und Ausfahrten zum bzw. vom Knotenpunkt auf, wobei Knotenarm 2 im dargestellten Ausführungsbeispiel nur eine Ausfahrt umfasst. Die Fahrbahn von Zufahrt 1 nach Ausfahrt 3 umfasst im beschriebenen Ausführungsbeispiel keine Messquerschnitte. An allen anderen Zu- und Ausfahrten des Knotenpunktes befinden sich Messquerschnitte mit Detektoren zum Erfassen von Einfahrt- Verkehrsstärken q"(n) und Ausfahrt-Verkehrsstärken q°ut(n) in vorgebbaren Messintervallen n. Grundelement des erfindungsgemäßen Verfahrens zur dynamischen Schätzung von Abbiegeraten ist eine geeignete Zerlegung des Netzausschnittes in Teilnetze. FIG 1 zeigt ein erstes Teilnetz fw, dessen Netzrand als Strichpunktlinie dargestellt ist und welches Messquerschnitte in der Ausfahrt 1 sowie in den relevanten Zufahrten 3 und 4 umfasst. Die Zufahrten 3 und 4 sind relevant, da auf ihnen Teilverkehrströme in das Teilnetz fw zufließen, welche durch die Ausfahrt 1 aus dem Teilnetz fw abfließen.FIG. 1 shows a network section of a road network having a node, for example a city road network in which the turning rates of the traffic flows are to be determined for the purpose of traffic control. The node has four node arms i (i = 1, ..., 4) with entrances and exits to or from the node, with node arm 2 in the illustrated embodiment includes only one exit. The roadway from driveway 1 to exit 3 does not include any measuring cross-sections in the exemplary embodiment described. At all other entrances and exits of the junction are measuring cross sections with detectors for detecting entrance traffic volumes q "(n) and exit traffic volumes q ° ut (n) in predetermined measuring intervals n. Basic element of the inventive method for dynamic estimation of turning rates 1 shows a first subnetwork fw whose network edge is shown as a dashed-dotted line and which measuring cross-sections in exit 1 as well as in the relevant access roads 3 and 4. The access roads 3 and 4 are relevant, because on them partial traffic flows into the subnetwork fw, which flow through the exit 1 from the subnetwork fw.
In FIG 2 ist das Teilnetz fw aus FIG 1 mit den zugehörigen Abbiegebeziehungen dargestellt. Die Abbiegerate m§"(k) gibt den Anteil der auf Zufahrt 3 gemessenen Verkehrsstärke qj n(k) an, der durch Ausfahrt 1 aus dem Teilnetz fw abfließt und damit zu der dort gemessenen Verkehrstärke q°ut(k) beiträgt. A- naloges gilt für die Abbiegerate mf"(k) in Bezug auf die Einfahrt-Verkehrsstärke q^fk) . Das Teilnetz fw (forward) modelliert die Abbiegebeziehungen also zeitlich vorwärts-bezogen. Bei sich langsam ändernden Verkehrsgeschehen kann für ein vorgegebenes Messintervall k die Ausfahrt-Verkehrsstärke q°ut(k) als gewichtete Summe der Einfahrt-Verkehrsstärken qj n(k) und q^fk) modelliert werden, wobei die Gewichtungsfaktoren den entsprechenden Abbiegeraten m§"(k) und mf"(k) entsprechen.FIG. 2 shows the subnetwork fw from FIG. 1 with the associated turn-off relationships. The turning rate m§ "(k) indicates the proportion of traffic volume q j n (k) measured on access 3, which flows out of subnetwork fw through exit 1 and thus contributes to the measured traffic q o ut (k) A - naloges applies to the turning rate mf "(k) in relation to the entrance traffic volume q ^ fk). The subnetwork fw (forward) models the turn-related relationships forward in time. For slowly changing traffic events, the exit traffic volume q ° ut (k) can be modeled as a weighted sum of the entrance traffic volumes q j n (k) and q ^ fk) for a given measurement interval k, the weighting factors corresponding to the corresponding turning rate m§ "(k) and mf" (k).
Allgemein kann bei r für eine Ausfahrt j relevanten Zufahrten i angesetzt werden:In general, r can be used at r for a j exit accesses relevant:
Von dieser Modellgleichung ausgehend können die Abbiegeraten mi"(k) mittels eines mathematischen Schätzverfahrens ermittelt werden. Erfindungsgemäß wird dem Schätzverfahren jedoch ein verallgemeinerter Zeitbezug zugrunde gelegt. Es werden außer dem aktuell betrachteten Messintervall k noch eine Mehrzahl von vorangehenden Messintervallen n = k-1, k-2, ... berücksichtigt. Für die vorwärts-bezogene Modellgleichung, in der insgesamt z Messintervalle einbezogen werden, ergibt sich daher:Starting from this model equation, the turning rates m i "(k) can be determined by means of a mathematical estimation method. According to the invention the estimation method used here, however, a generalized time reference. There are other than the currently considered measurement interval k or a plurality of preceding measuring intervals n = k-1 , k-2, ... For the forward-related model equation, in which a total of z measurement intervals are included, the result is:
q°ut(k) = ∑ ∑ <(k - l + D • qr(k -i + i) Die zu schätzenden Abbiegeraten mf"(k) ergeben sich als Summen der Abbiegeraten mf"(k - 1 + 1) über die berücksichtigten Messintervalle 1 = 1, ..., z:q ° ut (k) = Σ Σ <(k-1 + D • qr (k -i + i) The turning rates mf "(k) to be estimated are the sum of the turning rates mf" (k-1 + 1) over the measuring intervals 1 = 1,..., Z considered:
Durch diesen Ansatz wird das erfindungsgemäße Verfahren robust hinsichtlich der Reisezeiten zwischen den Messquerschnitten, bei hoher Genauigkeit und hinreichender Schnelligkeit .By this approach, the inventive method is robust in terms of travel times between the measuring cross sections, with high accuracy and sufficient speed.
Mit Vorteil wird dieser Ansatz erfindungsgemäß auch für zeitlich rückwärts-bezogene Teilnetze bw (backward) angewandt. FIG 1 zeigt ein solches Teilnetz bw, dessen Netzrand als Strichlinie dargestellt ist und welches Messquerschnitte in der Zufahrt 3 sowie in den relevanten Ausfahrten 1, 2 und 4 umfasst. Die Ausfahrten 1, 2 und 4 sind relevant, da auf ihnen Teilverkehrströme aus dem Teilnetz bw abfließen, welche durch die Zufahrt 3 in das Teilnetz bw eingeflossen sind. InAdvantageously, this approach is also applied according to the invention for time-backward-related subnetwork bw (backward). FIG. 1 shows such a subnetwork bw whose network edge is shown as a dashed line and which comprises measuring cross sections in the access road 3 and in the relevant exits 1, 2 and 4. The exits 1, 2 and 4 are relevant because flow on them partial traffic flows from the sub-network bw, which have flowed through the driveway 3 in the sub-network bw. In
FIG 3 ist dieses Teilnetz bw mit den zugehörigen Abbiegebeziehungen dargestellt. Die Abbiegerate m^fk) gibt den Anteil der in Ausfahrt 2 gemessenen Verkehrsstärke q^Ck) an, der über Zufahrt 3 in das Teilnetz bw eingeflossen ist und damit zu der dort gemessenen Verkehrstärke qj n(k) beiträgt. Analoges gilt für die Abbiegeraten m^fk) bzw. m^fk) in Bezug auf die Ausfahrt-Verkehrsstärke q°ut(k) bzw. q°ut(k) . Das Teilnetz bw modelliert die Abbiegebeziehungen also zeitlich rückwärtsbezogen. Allgemein kann bei s für eine Zufahrt j relevanten Ausfahrten i unter Berücksichtigung von z Messintervallen n wieder die gewichtete Summe angesetzt werden:FIG. 3 shows this subnetwork bw with the associated turn-off relationships. The turning rate m ^ fk) indicates the proportion of traffic volume q ^ Ck measured in exit 2, which has flowed into the subnetwork bw via driveway 3 and thus contributes to the traffic strength q j n (k) measured there. The same applies to the turning rates m ^ fk) or m ^ fk) in relation to the exit traffic q o ut (k) or q o ut (k). The subnetwork bw thus models the turn-related relationships backwards in time. In general, at s for departures j relevant exits i taking into account z measuring intervals n again the weighted sum can be applied:
<M = Σ Σ <<k + i - D • q°ut(k + 1-D<M = Σ Σ << k + i - D • q ° ut ( k + 1-D
Die zu schätzenden Abbiegeraten mf"(k) ergeben sich analog als Summen der Abbiegeraten m^"(k + 1 - 1) über die berücksichtigten Messintervalle 1 = 1, ..., z: m?(k) = ∑ 1O + 1 - D , j = 1 , . . . , sThe turning rates mf "(k) to be estimated are obtained analogously as sums of the turning rates m 1" (k + 1 - 1) over the measuring intervals 1 = 1,... m ? (k) = Σ 1 O + 1 - D, j = 1,. , , , s
1 = 11 = 1
Erfindungsgemäß wird als mathematisches Schätzverfahren ein erweitertes Kaiman-Filter zur Schätzung der Abbiegeraten eingesetzt .According to the invention, an extended Kalman filter for estimating the turning rates is used as the mathematical estimation method.
Falls die Messquerschnitte nahe beieinander liegen, etwa ein bis zwei dazwischen liegende Lichtsignalanlagen, genügt eine geringere Anzahl an zu berücksichtigenden Messintervallen, erfahrungsgemäß drei oder vier.If the measurement cross sections are close to each other, for example one to two intervening traffic signal systems, a smaller number of measurement intervals to be taken into account is sufficient, according to experience three or four.
Wird der allgemeine Ansatz mit z > 3 angewendet, kann die Modellgleichung auch in größeren Teilnetzen angewendet werden, deren Detektoren nur am Netzrand und nicht im Netzinneren ausgewertet werden. Es muss lediglich sichergestellt werden, dass im vorwärts-bezogenen Fall zu einem Ausfahrt- Messquerschnitt eines Teilnetzes alle relevanten Zufahrt- Messquerschnitte einbezogen werden, bzw. dass im rückwärtsbezogenen Fall zu einem Zufahrt-Messquerschnitt alle relevanten Ausfahrt-Messquerschnitte einbezogen werden.If the general approach with z> 3 is used, the model equation can also be applied in larger subnetworks whose detectors are evaluated only at the edge of the network and not inside the network. It merely has to be ensured that in the forward-related case to an exit measuring cross-section of a sub-network all relevant access measuring cross-sections are included, or that in the backward-related case to an access measuring cross section all relevant exit measuring cross sections are included.
Auf diese Weise kann in Teilnetzen begrenzter Größe eine direkte dynamische Schätzung von Quelle-Ziel-Strömen durchgeführt werden. Die Parameter des Kaiman-Filters, etwa Fehlervarianzen, müssen allerdings entsprechend angepasst werden und die Schätzqualität ist nicht so hoch wie bei nahe beieinander liegenden Messquerschnitten. Sinnvoll ist hier die Vergrößerung der Messintervalle, d.h. der Aggregationszeiträume, auf beispielsweise fünf Minuten, um die Störung des Schätzprozesses durch Rauschen zu reduzieren.In this way, in subnetworks of limited size, a direct dynamic estimation of source-destination streams can be performed. However, the parameters of the Kalman filter, such as error variances, need to be adjusted accordingly and the estimation quality is not as high as for closely spaced measurement cross sections. It makes sense here to increase the measurement intervals, i. the aggregation periods, for example, five minutes to reduce the noise of the estimation process.
Mit Vorteil wird der Schätzvorgang unterbrochen, wenn auf einem der für ein Teilnetz verwendeten Messquerschnitte eine Überlastung detektiert wird. Verschiedentlich wird für die Ermittlung von Staulängen in Zufahrten zu Lichtsignalanlagen ein Bilanzierungsansatz verwendet. Dieser bestimmt die Anzahl der Fahrzeuge, d.h. die Staulänge, in einer der Zufahrten durch Integration der Verkehrsstärken an den Endpunkten des Fahrbahnabschnitts über die Zeit. Wenn sich - was in der Praxis meist der Fall ist - nur ein Messquerschnitt an einem der beiden Endpunkte befindet, kann die Verkehrsstärke am anderen Endpunkt über erfindungsgemäß ermittelte Abbiegeraten geschätzt werden.Advantageously, the estimation process is interrupted if an overload is detected on one of the measuring cross sections used for a subnetwork. Occasionally, an accounting approach is used to determine traffic jams in access roads to traffic lights. This determines the number of vehicles, ie the traffic jam length, in one of the driveways by integrating the traffic volumes at the end points of the roadway section over time. If, as is usually the case in practice, only one measuring cross-section is located at one of the two end points, the traffic volume at the other end point can be estimated via turning rates determined according to the invention.
Ein weiterer Vorteil einer Anwendung des erfindungsgemäßen Verfahrens wird nachfolgend für den Fall beschrieben, dass sich ein Messquerschnitt an der Zufahrt zu einem Fahrbahnabschnitt befindet, während die Verkehrsstärke an der Ausfahrt aus dem Fahrbahnabschnitt über Abbiegeraten geschätzt wird, da sich an diesem Fahrbahnquerschnitt kein Messquerschnitt befindet. Die oben genannten Modellgleichungen zeichnen sich dadurch aus, dass die Abbiegeraten der Zu- bzw. Ausfahrten in Bezug zur Ausfahrt- bzw. Zufahrt-Verkehrsstärke geschätzt werden. Weisen die Messwerte der Verkehrsstärken proportionale Fehler auf, werden diese in den Abbiegeraten kompensiert. Die daraus errechneten Zufahrt-Verkehrsstärken sind konsistent zur gemessenen Ausfahrt-Verkehrsstärke. Auf diese Weise wird die Qualität von Bilanzierungen erheblich gesteigert, ohne dass spezielle Kalibrierungen für jeden bilanzierten Fahrbahnabschnitt erforderlich sind.A further advantage of an application of the method according to the invention is described below for the case where a measuring cross section is located at the entrance to a roadway section, while the traffic volume at the exit from the roadway section is estimated via turning rates, since there is no measuring cross section at this roadway cross section. The above model equations are characterized by the fact that the turn rates of the entrances and exits are estimated in relation to the exit or access traffic volume. If the measured values of the traffic intensities have proportional errors, these are compensated in the turn rates. The calculated access traffic volumes are consistent with the measured exit traffic. In this way, the quality of balances is significantly increased without the need for special calibrations for each balanced road section.
Für die weitere Beschreibung wird eine Matrix-Schreibweise der Schätzergebnisse eingeführt. Betrachtet man ein Messintervall k - nachfolgend nicht mehr dargestellt -, lassen sich die vorwärts-bezogenen Schätzungen in Matrix-Schreibweise zusammenfassen. Die Formel beschreibt, wie man anhand der vorwärts-bezogenen Abbiegeraten von den Zufahrt-Verkehrsstärken auf die Ausfahrt-Verkehrsstärken schließen kann:For the further description, a matrix notation of the estimation results is introduced. Considering a measuring interval k - not shown in the following -, the forward-related estimates can be summarized in matrix notation. The formula describes how to deduce the exit traffic volumes from the access traffic volumes on the basis of the forward-related turning rates:
qout = Mfw • qin In dieser Vektorgleichung sind qin bzw. qout Spaltenvektoren, deren Komponenten die Zufahrt-Verkehrsstärken q^n bzw. Ausfahrt-Verkehrsstärken q°ut aller Messquerschnitte in den Zu- bzw. Ausfahrten der Knotenarme i darstellen, während Mfw eine (n x n) -Matrix bedeutet, deren Elemente die Abbiegeraten mf" sind.q out = M fw • q in In this vector equation q in and q out are column vectors whose components represent the access traffic volumes q ^ n and exit traffic volumes q ° ut of all measurement cross sections in the entrances and exits of the node arms i, while M fw represents a (nxn) Matrix means whose elements are the turn rates mf ".
Entsprechend lässt sich für die rückwärts-bezogene Propagati- on formulieren:Correspondingly, it can be stated for the backward-related propagation:
Die Spaltenvektoren qin und qout umfassen dabei jeweils alle Messquerschnitte in der gleichen Reihenfolge, auch wenn es zwischen Komponenten q± der rechten Seite und q-, der linken Seite der Gleichung keine Fahrbeziehung gibt. Für das entsprechende Element der Matrix M gilt in diesem Fall ItI1-, = 0.The column vectors q in and q out comprise in this case all the measuring cross sections in the same sequence, even if there is no driving relationship between components q ± the right side and q, the left side of the equation. For the corresponding element of the matrix M, in this case ItI 1 - = 0.
Im Allgemeinen kann man davon ausgehen, dass Detektoren der Messquerschnitte die tatsächlichen Verkehrsstärken nur mit einer gewissen Genauigkeit erfassen. Dies kann sich aus ihrer Lage ergeben, z.B. Überfahrt durch Fahrzeuge zweier Spuren, aber auch durch einen defekten internen Kalibriervorgang, bei dem die Messwerte einen Drift erfahren, wie er sich im Laufe der Zeit in der Praxis häufig einstellt.In general, it can be assumed that detectors of the measuring cross-sections capture the actual traffic volumes only with a certain accuracy. This may result from their location, e.g. Passing through vehicles of two lanes, but also through a faulty internal calibration process, in which the measured values drift, as it often occurs in practice over time.
Geht man in erster Näherung von einer konstanten prozentualen Abweichung f1 = 1 + Δ aus, kann der Zusammenhang zwischen tatsächlicher Verkehrsstärke q^ und gemessener Verkehrsstärke Ig1 wie folgt formuliert werden:Assuming a first approximation of a constant percentage deviation f 1 = 1 + Δ, the relationship between the actual traffic volume q ^ and the measured traffic intensity Ig 1 can be formulated as follows:
q, = f, • q, r i = i r ... , nq, = f, • q, ri = ir ..., n
Wenn mf" bzw. m^" die geschätzten Abbiegeraten sind, gilt:If mf "or m ^" are the estimated turning rates, then:
q°ut = M fw • qin bzw . qin = Mbw • qout Die Korrekturfaktoren fx lassen sich in eine Diagonalmatrix F zusammenfassen :q ° ut = M fw • q in resp. q in = M bw • q out The correction factors f x can be summarized in a diagonal matrix F:
F = diag (fi, ... , flf ... , fn)F = diag (fi, ..., f lf ..., f n )
Unter der physikalisch sinnvollen Annahme, dass die Diagonalelemente fi ungleich Null sind, existiert die Inverse F"1 von F, die ebenfalls Diagonalenform hatUnder the physically reasonable assumption that the diagonal elements fi are not equal to zero, there exists the inverse F "1 of F, which also has diagonal form
F"1 = diag^"1, ... , f;1, ... , f;1) ,F "1 = diag ^ " 1 , ..., f; 1 , ..., f; 1 ),
so dass folgt :so that follows:
qout = (F • Mfw • F"1) • qin bzw . qin = (F • Mbw • F"1 ) • qout q out = (F • M fw • F '1 ) • q in or q in = (F • M bw • F ' 1 ) • q out
Aus dem Vergleich der Beziehungen für die tatsächlichen und geschätzten Verkehrsstärken lässt sich die Beziehung zwischen den realen und den geschätzten Matrizen der Abbiegeraten ableiten:From the comparison of the relationships for the actual and estimated traffic levels, the relationship between the real and the estimated matrices of the turn rates can be deduced:
Mfw = F Mfw ,-i bzw. Mbw = F Mbw -,-1M fw = F Mfw, -i and M bw = F Mbw -, - 1
Für die Elemente der geschätzten Matrizen ergibt sich darausFor the elements of the estimated matrices results from it
π£ = ^-m* bzw. m£ = ^ . A£π £ = ^ -m * or m £ = ^. A £
Für die realen Abbiegebeziehungen gilt unter der Annahme der Fahrzeugerhaltung vorwärts-bezogen betrachtet am Zufahrt- Messquerschnitt i bzw. rückwärts-bezogen betrachtet am Ausfahrt-Messquerschnitt iWith regard to the real turning relationships, assuming vehicle preservation, when viewed in the forward direction, at the access measuring cross section i or, viewed in the reverse direction, at the exit measuring cross section i
was unmittelbar zu den Bestimmungsgleichungen für die Elemente der Matrix F mit den Korrekturfaktoren führt: ∑ f: - KW 3 " f, = 0 bzw. ± f3 • π£ - f, = 0which leads directly to the equations of determination for the elements of the matrix F with the correction factors: Σ f : - K W 3 "f, = 0 or ± f 3 • π £ - f, = 0
D=I D=ID = I D = I
für alle Spalten i in Mfw bzw. Mbw , die nicht nur aus Nullen bestehen.for all columns i in M fw or M bw , which do not consist only of zeros.
Aus nin Zufahrt- und nout Ausfahrt-Messquerschnitten ergibt sich ein überbestimmtes homogenes Gleichungssystem für die U1 mit nin + nout ≥ n Gleichungen. Den Grenzfall nin + nout = n erhält man, wenn jeder einzelne Messquerschnitt nur in jeweils einer Richtung einen benachbarten Messquerschnitt aufweist. Dieser Fall tritt in der Praxis dann auf, wenn Messquerschnitte nur am Netzrand vorliegen. Der einfachste Fall ist der einer einzelnen Lichtsignalanlage, in der alle relevanten Zu- und Ausfahrten erfasst werden.From n in access and out n exit measurement cross sections results in an over-determined homogeneous system of equations for the U 1 with n in + n out ≥ n equations. The limit case n in + n out = n is obtained if each individual measuring cross-section only has an adjacent measuring cross-section in one direction. This case occurs in practice when measuring cross-sections are present only at the edge of the network. The simplest case is that of a single traffic signal system in which all relevant entrances and exits are recorded.
Um die physikalisch nicht sinnvolle Triviallösung F = O auszuschließen, wird die folgende Anforderung an die Lösung ergänzt :In order to exclude the physically irrational trivial solution F = O, the following requirement is added to the solution:
fe - I)2 = 0fe - I) 2 = 0
Trifft die Annahme zu, dass Abweichungen lediglich zu proportional verfälschten Messwerten führen, entsteht hierdurch keine Einschränkung, da für jede Lösung F des homogenen Gleichungssystems auch F = X-F' eine Lösung darstellt.If the assumption applies that deviations only lead to proportionally falsified measured values, this results in no restriction, since for each solution F of the homogeneous system of equations F = X-F 'also represents a solution.
Aus letzten Gleichungen ergibt sich letztendlich die Formulierung eines geeigneten nichtlinearen Optimierungsproblems, wonach die KenngrößeUltimately, the formulation of a suitable nonlinear optimization problem results from the last equations, according to which the parameter
W1 • ∑ (sf f + W1 • ∑ (sf)2 + W2 • (fl - I)2 p _ i≡l i≡l nW 1 • Σ (sf f + W 1 • Σ (sf) 2 + W 2 • ( fl - I) 2 p _ i≡li≡ln
mit sfwWith s fw
= Σ D = I f3 • K ~ fi fÜr Σ D = I K ≠ 0 ' S°nSt S"W = 0 = Σ D = I f 3 • K ~ f i for Σ D = IK ≠ 0 ' S ° nSt S " W = 0
sf = ∑ f, • K - f, für £ m£ ≠ 0 , sonst sf = 0sf = Σ f, • K - f, for £ m £ ≠ 0, otherwise sf = 0
D =I D=ID = I D = I
zu minimieren ist, wobei Wi und W2 wählbare Gewichtungen darstellen.to minimize, where Wi and W2 represent selectable weights.
Da der erste Faktor fi willkürlich gewählt wurde, obwohl er möglicherweise tatsächlich ungleich 1 ist, kann abschließend noch eine Korrektur der ermittelten Lösung F' erfolgen. Geht man davon aus, dass in einem Netz weniger als die Hälfte aller Messquerschnitte zu viele Fahrzeuge (fi > 1) und weniger als die Hälfte aller Messquerschnitte zu wenige Fahrzeuge (fx < 1) schätzen, kann als Referenz der Medianwert von F verwendet werden, da dieser bei 1 liegen muss. Als Lösung des Problems kann dann verwendet werden:Since the first factor fi has been arbitrarily chosen, even though it may not actually be equal to 1, a final correction of the determined solution F 'can take place. Assuming that in a net less than half of all measuring sections estimate too many vehicles (fi> 1) and less than half of all measuring sections too few vehicles (f x <1), the median value of F can be used as a reference because it has to be 1. As a solution to the problem can then be used:
11
F = λ • F' mit λ = median(F)F = λ • F 'with λ = median (F)
Das formulierte Optimierungsproblem kann formal auf ein Netz als Ganzes angewendet werden. Dies hat allerdings zwei Nachteile: Einerseits nimmt der Rechenaufwand für die Optimierung mit der Zahl der Messquerschnitte überproportional zu. Andererseits führt die Tatsache, dass im Netz viele Messquerschnitte in beide Richtungen weitere Messquerschnitte haben, zu Nivellierungseffekten der geschätzten Korrekturfakto¬The formulated optimization problem can be formally applied to a network as a whole. However, this has two disadvantages: On the one hand, the computational effort for the optimization increases disproportionately with the number of measurement cross sections. On the other hand, the fact that many measuring cross sections in the network have further measuring cross sections in both directions leads to leveling effects of the estimated correction factor
Um diesen Nachteil zu umgehen, kann das Verfahren verfeinert werden, indem das Netz geeignet zerlegt wird. Als günstig hat sich die Anwendung auf Teilnetze erwiesen, die so definiert sind, dass von jedem Punkt innerhalb eines solchen Netzes nur solche Messstellen Teil dieses Teilnetzes sind, die über Netzkanten direkt erreichbar sind. Jedes solche Teilnetz stellt quasi ein Inselnetz mit Zu- und Ausfahrten in das üb- rige Netz dar. Es fasst alle vorwärts- bzw. rückwärtsbezogenen Teilnetze zusammen, die gleiche Zu- bzw. Ausfahrt- Messquerschnitte haben. FIG 4 zeigt ein solches Inselnetz, das bei üblicher Detektorausstattung an Knotenpunkten in der Regel um die Lichtsignalanlage herum entsteht. Es weist Messquerschnitte nur am Netzrand auf, der in FIG 4 mit einer Strichdoppelpunktlinie veranschaulicht ist.To circumvent this disadvantage, the method can be refined by properly decomposing the network. The application to subnetworks has proven to be favorable, which are defined such that from any point within such a network only those measuring points are part of this subnetwork which can be reached directly via network edges. Each such subnet constitutes, as it were, an island network with entrances and exits into the It summarizes all forward or backward related subnets that have the same entry and exit measurement cross sections. FIG. 4 shows such a stand-alone grid, which generally arises around the traffic signal system at the usual detector equipment at nodes. It has measuring cross-sections only at the network edge, which is illustrated in FIG. 4 with a double-dashed line.
Die eben beschriebenen Nachteile sind damit aufgehoben. Gleichzeitig liefert die Kenngröße P, die sich aus der Lösung des Optimierungsproblems für ein Inselnetz ergibt, weitergehende Hinweise über eventuelle Detektorstörungen: Der Wert P der Lösung des Optimierungsproblems für ein Inselnetz geht im Idealfall gegen Null, wenn die Abbiegeraten exakt geschätzt sind und zwischen den Messquerschnitten keine Fahrzeugverluste (Senken) oder -zunahmen (Quellen) auftreten, sowie alle Messfehler proportionaler Art sind. Ergibt sich für P ein Wert erheblich größer als Null, z.B. Zwei, ist dies ein Indiz für ungenau geschätzte Abbiegeraten, nicht erfasste Zu- und Ausfahrten oder Messfehler anderer als proportionaler Art. Typischerweise entsteht ein solcher Fehler dann, wenn nicht alle Verkehrsflüsse in den relevanten Zu- und Ausfahrten eines solchen Teilnetzes erfasst sind oder z.B. durch Fehlversorgungen / -Verdrahtungen Messquerschnitte den Strecken falsch zugewiesen wurden.The disadvantages just described are thus repealed. At the same time, the parameter P resulting from the solution to the optimization problem for an island grid provides further information about possible detector disturbances: the value P of the solution to the optimization problem for an island grid ideally approaches zero if the turning rates are exactly estimated and between the measuring cross sections no vehicle losses (sinks) or increases (sources) occur, and all measurement errors are proportional. If P results in a value significantly greater than zero, e.g. Two, this is an indication of inaccurately estimated turning rates, missed entrances and exits, or measurement errors other than proportional. Typically, such an error arises when not all traffic flows in the relevant entrances and exits of such a subnetwork are detected or e.g. due to incorrect supply / wiring, measuring sections were assigned incorrectly to the sections.
FIG 5 zeigt den Verlauf der erfindungsgemäß ermittelten rück- wärts-bezogenen Abbiegeraten ItI1 (k) für die Verkehrsströme zwischen Knotenarm 3 und Knotenarmen 4, 1 und 2. Dazu wird eine Liste mit Korrekturfaktoren f1 sowie der Wert der Kenngröße P ausgegeben, was nicht dargestellt ist.5 shows the course of the backward-related turning rates ItI 1 (k) determined according to the invention for the traffic flows between node arm 3 and node arms 4, 1 and 2. For this purpose, a list with correction factors f 1 and the value of the parameter P are output not shown.
An allen Messquerschnitten zwischen Inselnetzen werden mit diesem Ansatz jeweils zwei Werte für den Messfehler geschätzt. In einem abschließenden Ausgleichsverfahren können diese über einen weiteren Optimierungsschritt paarweise aneinander angeglichen werden, indem alle Messfehler U1 aller Teilnetze mit Teilnetz-spezifischen Korrekturfaktoren multip- liziert werden. Diese Korrekturfaktoren lassen das Verhältnis der fx innerhalb jedes Teilnetzes unverändert, führen aber zu einem Abgleich gemeinsamer Schätzfehler zwischen den Teilnetzen.At each measuring cross section between isolated networks, two values for the measurement error are estimated using this approach. In a final compensation process, these can be matched to one another in pairs by means of a further optimization step by multiplying all measurement errors U 1 of all subnetworks with subnetwork-specific correction factors. be published. These correction factors leave the ratio of the f x within each subnet unchanged, but result in a balancing of common estimation errors between the subnets.
Schließlich lassen sich alle gemessenen Verkehrsstärken überFinally, all measured traffic volumes can be over
qf' = f. • q, q f '= f . • q ,
und alle Abbiegeraten überand all turning rates over
kalibrieren, was die Schätzung der Verkehrslage qualitativ verbessert .calibrate, which qualitatively improves the estimate of the traffic situation.
Eine weitere Anwendung der Ergebnisse ist die Bildung von Ersatzwerten, falls Detektoren oder ganze Messquerschnitte ausgefallen sind. Solche Ausfälle werden entweder schon in der Hardware der Geräte erkannt und weitergemeldet oder sie lassen sich über einfache Plausibilitätskontrollen erfassen. In diesem Fall können Ersatzwerte einfach anhand der gemessenen und kalibrierten Verkehrsstärken und Abbiegeraten der umgebenden Messquerschnitte gebildet werden. Dies stellt einen erheblichen Qualitätssprung gegenüber bekannten Verfahren dar, die fehlende Messwerte im einfachen Fall lediglich durch vorab fest versorgte zeitabhängige Werte ersetzen oder - sehr aufwendig - diese über vorab erfasste, Tagestyp-spezifische Ganglinien und den aktuellen Zeitpunkt bestimmen.Another application of the results is the formation of substitute values if detectors or entire measuring sections have failed. Such failures are either already detected in the hardware of the devices and reported further or they can be detected via simple plausibility checks. In this case, substitute values can simply be formed based on the measured and calibrated traffic intensities and turning rates of the surrounding measuring cross sections. This represents a considerable leap in quality over known methods, which in the simple case simply replace missing measured values with previously-dependent time-dependent values or-very costly-determine these via pre-recorded, day-type-specific hydrographs and the current time.
Die mit dem beschriebenen Verfahren geschätzten Abbiegeraten weisen eine Genauigkeit auf, die sowohl für die Verkehrslageschätzung als auch für die Weiterverwendung in adaptiven Netzsteuerungsverfahren ausreicht, und die die Schätzqualität konventioneller, umlegungsbasierter Verfahren weit übertreffen. Die Schätzung der Korrekturfaktoren U1 für alle Messquerschnitte ermöglicht einerseits eine online-Korrektur der Abbiegeraten und Verkehrsstärke-Zählwerte. Das damit mögliche Verfahren zur netzbezogenen Ersatzwertbildung kommt ohne Ganglinien oder vorab versorgte Default-Zählwerte aus und funktioniert auch bei mehreren dicht beieinander liegenden Detektorausfällen, da in diesem Fall Ersetzungen auch rekursiv vorgenommen werden können.The turn rates estimated by the described method are accurate enough for both traffic estimation and reuse in adaptive network control techniques and far surpassing the estimation quality of conventional metering based methods. The estimation of the correction factors U 1 for all measuring cross sections on the one hand enables an online correction of the turning rates and traffic volume counts. The thus possible method for network-related replacement value formation works without hydrographs or pre-supplied default count values and also works with several closely spaced detector failures, since in this case substitutions can also be made recursively.
Andererseits entstehen aus der Konsistenzprüfung wertvolle Hinweise für die Wartung des Detektornetzes. Entsprechende Benachrichtigungsmechanismen ermöglichen dem Wartungsdienst einer Stadt eine effiziente und beschleunigte Reparatur defekter Detektoren, verbunden mit nennenswerten Kosteneinsparungen. On the other hand, the consistency check provides valuable information for the maintenance of the detector network. Corresponding notification mechanisms allow a city's maintenance service to efficiently and quickly repair defective detectors, with significant cost savings.

Claims

Patentansprüche claims
1. Verfahren zur Ermittlung von Abbiegeraten (mf*(k)) in einem1. A method for determining turning rates (mf * (k)) in one
Straßennetz,Road network,
- wobei an Messquerschnitten in vorgebbaren Messintervallen (k) Verkehrsstärken (q(k)) erfasst werden,where traffic intensities (q (k)) are detected at measurement cross sections in predeterminable measurement intervals (k),
- wobei für wenigstens ein vorwärts-bezogenes Teilnetz (fw) des Straßennetzes, bei welchem Messquerschnitte an einer Ausfahrt (j) und an Zufahrten (i) des Teilnetzes (fw) berücksichtigt werden, eine Modellgleichung aufgestellt wird, in der die Ausfahrt-Verkehrsstärke (q°ut(k)) als gewichtete Summe der Zufahrt-Verkehrsstärken (q"(k)) angesetzt wird und diefor at least one forward-related subnetwork (fw) of the road network, in which measuring cross sections are taken into account at an exit (j) and at accesses (i) of the subnetwork (fw), a model equation is set up in which the exit traffic volume ( q ° ut (k)) is the weighted sum of the access traffic volumes (q "(k)) and the
Gewichtungsfaktoren den vorwärts-bezogenen Abbiegeraten (mf"(k)) entsprechen, welche jeweils den Anteil einer Zufahrt- Verkehrsstärke (q"(k)) angeben, der durch die berücksichtigte Ausfahrt (j) abfließt,Weighting factors correspond to the forward-related turning rate (mf "(k)), which respectively indicate the proportion of an access traffic volume (q" (k)) that flows through the considered exit (j),
- und wobei ausgehend von der Modellgleichung mittels eines mathematischen Schätzverfahrens die vorwärts-bezogenen Abbiegeraten (mf"(k)) berechnet werden, d a d u r c h g e k e n n z e i c h n e t,- and where, based on the model equation, the forward-related turning rates (mf "(k)) are calculated by means of a mathematical estimation method, d a d u c h e c e n e c e n e,
- dass in der gewichteten Summe für die Ausfahrt- Verkehrsstärke (q°ut(k)) eines gegebenen Messintervalls (k) die Zufahrt-Verkehrsstärken (q"(n)) einer Mehrzahl (z) vorangehender Messintervalle (n = k, k-1, k-2, ..., k-z+1) berücksichtigt werden,- that in the weighted sum for the exit traffic volume (q ° ut (k)) of a given measurement interval (k), the access traffic intensities (q "(n)) of a plurality of (z) preceding measurement intervals (n = k, k) 1, k-2, ..., k-z + 1),
- wobei sich eine zu ermittelnde vorwärts-bezogene Abbiegerate (mf"(k)) als Summe aus den entsprechenden Abbiegeraten- Where to be determined forward-related turning rate (mf "(k)) as the sum of the corresponding turning rates
(mf"(n)) der in der Modellgleichung berücksichtigten Messintervalle (n = k, k-1, k-2,..., k-z+1) ergibt.(mf "(n)) gives the measurement intervals considered in the model equation (n = k, k-1, k-2, ..., k-z + 1).
2. Verfahren nach Anspruch 1, d a d u r c h g e k e n n z e i c h n e t, dass für wenigstens ein rückwärts-bezogenes Teilnetz (bw) des Straßennetzes, bei welchem Messquerschnitte an einer Zufahrt (j) und an Ausfahrten (i) des Teilnetzes (bw) berücksichtigt werden, eine Modellgleichung aufgestellt wird, in der die Zufahrt- Verkehrsstärke (q"(k)) als gewichtete Summe der Ausfahrt- Verkehrsstärken (q°ut(k)) angesetzt wird und die Gewichtungsfaktoren den rückwärts-bezogenen Abbiegeraten (m^"(k)) entsprechen, welche jeweils den Anteil einer Ausfahrt- Verkehrsstärke (q°ut(k)) angeben, der durch die berücksichtigte Zufahrt (j) eingeflossen ist,2. The method according to claim 1, characterized in that for at least one backward-related subnetwork (bw) of the road network, in which measuring cross sections at an access (j) and at exits (i) of the subnetwork (bw) are taken into account, a model equation is set up in which the access volume (q "(k)) is the weighted sum of the Traffic intensities (q ° ut (k)) and the weighting factors correspond to the reverse-related turning rates (m ^ "(k)) which respectively indicate the proportion of an exit traffic volume (q ° ut (k)) determined by the considered access (j),
- wobei ausgehend von der Modellgleichung mittels eines mathematischen Schätzverfahrens die Abbiegeraten (m^"(k)) berechnet werden,in which, based on the model equation by means of a mathematical estimation method, the turning rates (m ^ "(k)) are calculated,
- wobei in der gewichteten Summe für die Zufahrt- Verkehrsstärke (q"(k)) eines gegebenen Messintervalls (k) die- wherein in the weighted sum for the access traffic volume (q "(k)) of a given measurement interval (k) the
Ausfahrt-Verkehrsstärken (q°ut(n)) einer Mehrzahl (z) nachfolgender Messintervalle (n = k, k+1, k+2, ..., k+z-1) berücksichtigt werden,Exit traffic intensities (q ° ut (n)) of a plurality of (z) subsequent measurement intervals (n = k, k + 1, k + 2,..., K + z-1) are taken into account,
- und wobei sich eine zu ermittelnde rückwärts-bezogene Abbiegerate (m^"(k)) als Summe aus den entsprechenden Abbiegeraten (m^"(n)) der in der Modellgleichung berücksichtigten Messintervalle (n = k, k+1, k+2,..., k+z-1) ergibt.- and where a backward-related turning rate (m ^ "(k)) to be determined as the sum of the corresponding turning rates (m ^" (n)) of the measuring intervals considered in the model equation (n = k, k + 1, k + 2, ..., k + z-1).
3. Verfahren nach Anspruch 1 oder 2, d a d u r c h g e k e n n z e i c h n e t, dass bei dem mathematischen Schätzverfahren ein erweitertes Kaiman-Filter angewandt wird.3. Method according to claim 1 or 2, characterized in that in the mathematical estimation method an extended Kalman filter is used.
4. Verfahren nach einem der Ansprüche 1 bis 3, d a d u r c h g e k e n n z e i c h n e t, dass das Schätzverfahren unterbrochen wird, wenn an einem Messquerschnitt eine Verkehrsüberlastung erfasst wird.4. Method according to one of claims 1 to 3, characterized in that the estimation method is interrupted if a traffic congestion is detected at a measuring cross section.
5. Verfahren zur Ermittlung von Abbiegeraten an einem Knotenpunkt des Straßennetzes, d a d u r c h g e k e n n z e i c h n e t, dass vorwärtsbezogene und/oder rückwärts-bezogene Teilnetze (fw, bw) um den Knotenpunkt betrachtet werden, bei welchen Messquerschnitte in den Zufahrten (j, i) und Ausfahrten (i, j) des Knotenpunktes berücksichtigt werden, und dass die Abbiegeraten ( mf"(k) , mf"(k) ) gemäß einem Verfahren nach einem der Ansprüche 1 bis 4 ermittelt werden. 5. A method for determining turning rates at a node of the road network, characterized in that forward-related and / or backward-related subnets (fw, bw) are considered around the node, in which measuring cross sections in the driveways (j, i) and exits (i , j) of the node, and that the turning rates (mf "(k), mf" (k)) are determined according to a method according to any one of claims 1 to 4.
6. Verfahren zur Ermittlung von Quelle-Ziel-Verkehrsströmen eines Teilnetzes, d a d u r c h g e k e n n z e i c h n e t, dass für die Zufahrten und Ausfahrten des Teilnetzes die Abbiegeraten gemäß einem Verfahren nach einem der Ansprüche 1 bis 4 ermittelt werden, wobei Messquerschnitte nur am Rand des Teilnetzes und nicht in seinem Inneren berücksichtigt werden, und dass aus den ermittelten Abbiegeraten und den erfassten Verkehrsstärken die Quelle-Ziel-Verkehrsströme für dieses Teilnetz berechnet werden.6. A method for determining source-destination traffic flows of a subnetwork, characterized in that for the accesses and exits of the subnetwork, the turning rates are determined according to a method of one of claims 1 to 4, wherein measuring cross sections only at the edge of the subnet and not in his Be considered inside, and that are calculated from the calculated turning rates and the detected traffic strengths, the source-destination traffic flows for this subnet.
7. Verfahren nach Anspruch 6, d a d u r c h g e k e n n z e i c h n e t, dass die Anzahl (z) der berücksichtigten Messintervalle (n) mit zunehmender Größe des betrachteten Teilnetzes vergrößert wird.7. The method as claimed in claim 6, wherein the number (z) of the measuring intervals (n) taken into account is increased with increasing size of the considered subnetwork.
8. Verfahren nach Anspruch 6 oder 7, d a d u r c h g e k e n n z e i c h n e t, dass die berücksichtigten Messintervalle (n) mit zunehmender Größe des betrachteten Teilnetzes verlängert werden.8. Method according to claim 6 or 7, characterized in that the measuring intervals (n) taken into account are lengthened with increasing size of the considered subnetwork.
9. Verfahren zur Ermittlung der Verkehrsstärke an einem Fahrbahnquerschnitt eines Straßennetzes, d a d u r c h g e k e n n z e i c h n e t, dass für ein Teilnetz des Straßennetzes, dessen eine Zu- oder Ausfahrt den Fahrbahnquerschnitt und dessen andere Zu- und/oder Ausfahrten Messquerschnitte aufweisen, gemäß einem Verfahren nach einem der Ansprüche 1 bis 4 ermittelte Abbiegeraten bereitgestellt werden, und dass aus den bereitgestellten Abbiegeraten und den an den Messquerschnitten erfassten Verkehrsstärken der anderen Zu- und/oder Ausfahrten die Verkehrstärke am Fahrbahnquerschnitt der einen Zu- bzw. Ausfahrt berechnet wird.9. A method for determining the traffic volume on a road cross section of a road network, characterized in that for a subnetwork of the road network whose one entrance or exit the road cross section and whose other entrances and / or exits measuring cross-sections, according to a method according to one of claims 1 are provided to 4 determined turning rates, and that the traffic strength is calculated on the road cross section of an entrance or exit from the provided turning data and recorded at the measuring cross sections traffic levels of the other entrances and / or exits.
10. Verfahren nach Anspruch 9, d a d u r c h g e k e n n z e i c h n e t, dass die für den Fahrbahnquerschnitt ermittelte Verkehrsstärke als Ersatzwert für einen fehlerhaften oder ausgefallenen Messquerschnitt herangezogen wird.10. The method according to claim 9, characterized in that the determined for the roadway traffic volume as a substitute value is used for a faulty or failed measuring cross-section.
11. Verfahren zur Ermittlung der Fahrzeuganzahl innerhalb eines Fahrbahnabschnitts, an dessen ersten Endpunkt die Verkehrsstärke an einem Messquerschnitt erfasst wird und an dessen zweiten Endpunkt kein Messquerschnitt angeordnet ist, d a d u r c h g e k e n n z e i c h n e t, dass die Verkehrsstärke am zweiten Endpunkt gemäß einem Verfahren nach Anspruch 9 oder 10 und daraus die Fahrzeuganzahl im Fahrbahnabschnitt durch zeitliche Integration der Differenz aus der in den Fahrbahnabschnitt zufließenden Verkehrsstärke und der aus diesem abfließenden Verkehrsstärke ermittelt wird.11. A method for determining the number of vehicles within a lane section, at the first end point, the traffic intensity is detected at a measuring cross section and at the second end point no measuring cross section is arranged, characterized in that the traffic strength at the second end point according to a method according to claim 9 or 10 and therefrom the number of vehicles in the lane section is determined by temporal integration of the difference between the traffic volume flowing into the lane section and the traffic volume flowing therefrom.
12. Verfahren zur Bestimmung von Korrekturfaktoren für gemäß einem Verfahren nach einem der Ansprüche 1 bis 4 ermittelte Abbiegeraten, d a d u r c h g e k e n n z e i c h n e t, dass zunächst aus der Fahrzeugerhaltung der tatsächlichen vorwärts- und rück- wärts-bezogenen Abbiegeraten für die zu bestimmenden Korrekturfaktoren (fi) ein homogenes Gleichungssystem aufgestellt wird, dass dann aus dem homogenen Gleichungssystem zusammen mit einer die triviale Lösung ausschließenden Nebenbedingung ein Optimierungsproblem gewonnen wird, und dass die Korrekturfaktoren (fi) sich als Lösung des Optimierungsproblems ergeben.12. A method for determining correction factors for determined according to a method according to one of claims 1 to 4 turning rates, characterized in that first from the vehicle preservation of the actual forward and backward-related turning rates for the correction factors to be determined (fi) a homogeneous system of equations is set up, that then an optimization problem is obtained from the homogeneous system of equations together with a secondary condition excluding the trivial solution, and that the correction factors (fi) result as a solution to the optimization problem.
13. Verfahren nach Anspruch 12, d a d u r c h g e k e n n z e i c h n e t, dass die bestimmten Korrekturfaktoren (fx) durch ihren Medianwert dividiert werden.13. The method according to claim 12, characterized in that the determined correction factors (f x ) are divided by their median value.
14. Verfahren nach Anspruch 12 oder 13, d a d u r c h g e k e n n z e i c h n e t, dass ein betrachtetes Teilnetz des Straßennetzes in Inselnetze zerlegt wird, wobei jedes Inselnetz nur Messquerschnitte an seinem Netzrand umfasst, und dass Korrekturfaktoren (f±) für die Inselnetze bestimmt werden. 14. The method according to claim 12 or 13, characterized in that a considered subnetwork of the road network is decomposed into island networks, each island network comprising only measuring cross sections at its network edge, and that correction factors (f ± ) are determined for the island grids.
15. Verfahren nach einem der Ansprüche 12 bis 14, d a d u r c h g e k e n n z e i c h n e t, dass überprüft wird, ob die bestimmten Korrekturfaktoren (f±) innerhalb eines vorgebbaren Wertebereiches liegen.15. The method according to any one of claims 12 to 14, characterized in that it is checked whether the determined correction factors (f ± ) are within a predefinable value range.
16. Verfahren nach einem der Ansprüche 12 bis 15, d a d u r c h g e k e n n z e i c h n e t, dass bei der Lösung des Optimierungsproblems eine Kenngröße (P) berechnet wird, deren Wert als Maß für die Schätzqualität der Abbiegeraten (m*"(k), mf"(k) ) verwendet wird.16. The method according to any one of claims 12 to 15, characterized in that in the solution of the optimization problem, a parameter (P) is calculated whose value as a measure of the estimated quality of the turning rates (m * "(k), mf" (k)) is used.
17. Verfahren nach Anspruch 16, d a d u r c h g e k e n n z e i c h n e t, dass eine Fehlermeldung ausgegeben wird, wenn der Wert der Kenngröße (P) eine vorgebbare Schranke übersteigt.17. Method according to claim 16, wherein a fault message is output if the value of the parameter exceeds a predeterminable limit.
18. Verfahren nach einem der Ansprüche 12 bis 17, d a d u r c h g e k e n n z e i c h n e t, dass für jeden Messquerschnitt, den sich zwei aneinander grenzende Inselnetze teilen, zwei Korrekturfaktoren bestimmt werden und dass die Korrekturfaktoren je eines Inselnetzes derart skaliert werden, dass sich die Korrekturfaktoren gemeinsamer Messquerschnitte einander angleichen.18. Method according to claim 12, wherein two correction factors are determined for each measuring cross section shared by two adjoining island networks, and that the correction factors of each island network are scaled such that the correction factors of common measuring cross sections are equal to one another.
19. Verfahren nach einem der Ansprüche 12 bis 18, d a d u r c h g e k e n n z e i c h n e t, dass die an den Messquerschnitten erfassten Verkehrsstärken (q"(k), q°ut(k)) und die durch Schätzung ermittelten Abbiegeraten (mf"(k), mf"(k) ) mittels der Korrekturfaktoren (f±) kalibriert werden. 19. The method according to any one of claims 12 to 18, characterized in that the detected at the measuring cross sections traffic volumes (q "(k), q ° ut (k)) and the determined by estimation turning rates (mf" (k), mf "( k)) are calibrated by means of the correction factors (f ± ).
EP06763261A 2005-05-31 2006-05-24 Method for determination of turning rates in a traffic network Not-in-force EP1749288B1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE102005024953A DE102005024953A1 (en) 2005-05-31 2005-05-31 Method for determining turning rates in a road network
PCT/EP2006/062571 WO2006128819A1 (en) 2005-05-31 2006-05-24 Method for determination of turning rates in a traffic network

Publications (2)

Publication Number Publication Date
EP1749288A1 true EP1749288A1 (en) 2007-02-07
EP1749288B1 EP1749288B1 (en) 2008-02-27

Family

ID=36939069

Family Applications (1)

Application Number Title Priority Date Filing Date
EP06763261A Not-in-force EP1749288B1 (en) 2005-05-31 2006-05-24 Method for determination of turning rates in a traffic network

Country Status (8)

Country Link
US (1) US7894979B2 (en)
EP (1) EP1749288B1 (en)
CN (1) CN101006481B (en)
AT (1) ATE387695T1 (en)
DE (2) DE102005024953A1 (en)
DK (1) DK1749288T3 (en)
NO (1) NO20065797L (en)
WO (1) WO2006128819A1 (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4770858B2 (en) * 2008-03-28 2011-09-14 アイシン・エィ・ダブリュ株式会社 Signalized intersection information acquisition apparatus, signalized intersection information acquisition method, and signalized intersection information acquisition program
DE102008022349A1 (en) * 2008-05-02 2009-11-12 Deutsches Zentrum für Luft- und Raumfahrt e.V. Method and device for determining tailback lengths at traffic lights
US8798897B2 (en) 2010-11-01 2014-08-05 International Business Machines Corporation Real-time traffic analysis through integration of road traffic prediction and traffic microsimulation models
CA2955961A1 (en) 2014-07-28 2016-02-04 Econolite Group, Inc. Self-configuring traffic signal controller
DE102017208854A1 (en) * 2017-05-24 2018-11-29 Volkswagen Aktiengesellschaft A method, apparatus and computer readable storage medium having instructions for determining applicable traffic rules for a motor vehicle
KR20210144411A (en) * 2020-05-22 2021-11-30 현대자동차주식회사 Apparatus and method for recognizing driving lane of vehicle

Family Cites Families (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3663937A (en) * 1970-06-08 1972-05-16 Thiokol Chemical Corp Intersection ingress-egress automatic electronic traffic monitoring equipment
US3920967A (en) * 1974-02-22 1975-11-18 Trw Inc Computerized traffic control apparatus
DE3870709D1 (en) 1987-05-27 1992-06-11 Siemens Ag METHOD FOR MEASURING THE INTENSITY OF ROAD TRAFFIC.
JP2712844B2 (en) * 1990-04-27 1998-02-16 株式会社日立製作所 Traffic flow measurement device and traffic flow measurement control device
JP3414843B2 (en) * 1993-06-22 2003-06-09 三菱電機株式会社 Transportation control device
DE4408547A1 (en) * 1994-03-14 1995-10-12 Siemens Ag Process for traffic detection and traffic situation detection on highways, preferably motorways
DE4415737A1 (en) * 1994-05-04 1995-11-23 Siemens Ag Procedure for the ongoing measurement of the current traffic flows at a traffic junction
JP3435623B2 (en) * 1996-05-15 2003-08-11 株式会社日立製作所 Traffic flow monitoring device
DE19647127C2 (en) * 1996-11-14 2000-04-20 Daimler Chrysler Ag Process for automatic traffic monitoring with dynamic analysis
US5921701A (en) * 1997-06-25 1999-07-13 Clayton; Robert F. Traffic interchange
DE19729914A1 (en) 1997-07-04 1999-01-07 Mannesmann Ag Process for the analysis of a traffic network, traffic analysis, traffic forecast as well as creation of a historical traffic database and traffic analysis and forecasting center
DE19835979B4 (en) * 1998-08-08 2005-01-05 Daimlerchrysler Ag Method for monitoring traffic conditions and vehicle inflow control in a road network
KR20000025578A (en) * 1998-10-13 2000-05-06 윤종용 Method for modeling traffic road and method for recognizing traffic lane using the same
DE19944075C2 (en) * 1999-09-14 2002-01-31 Daimler Chrysler Ag Traffic condition monitoring method for a traffic network with effective bottlenecks
DE10036789A1 (en) * 2000-07-28 2002-02-07 Daimler Chrysler Ag Method for determining the traffic condition in a traffic network with effective bottlenecks
DE10136646A1 (en) 2001-07-20 2003-02-06 Ddg Ges Fuer Verkehrsdaten Mbh Method for balancing of vehicle flows from input traffic data relating to a particular road stretch, whereby an instantaneous flow can be determined and from it useful traffic data that can be quickly transmitted to drivers
US20040088392A1 (en) * 2002-03-18 2004-05-06 The Regents Of The University Of California Population mobility generator and simulator
US7135989B2 (en) * 2004-09-27 2006-11-14 Gregory Fife Parsons Parallel flow vehicle turn system for traffic intersections

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See references of WO2006128819A1 *

Also Published As

Publication number Publication date
DK1749288T3 (en) 2008-06-23
US7894979B2 (en) 2011-02-22
NO20065797L (en) 2007-03-05
EP1749288B1 (en) 2008-02-27
CN101006481B (en) 2010-10-06
CN101006481A (en) 2007-07-25
ATE387695T1 (en) 2008-03-15
US20080069000A1 (en) 2008-03-20
DE502006000385D1 (en) 2008-04-10
DE102005024953A1 (en) 2006-12-07
WO2006128819A1 (en) 2006-12-07

Similar Documents

Publication Publication Date Title
DE69329119T2 (en) PREDICTION METHOD FOR ROAD TRAFFIC PARAMETERS
EP1749288B1 (en) Method for determination of turning rates in a traffic network
EP2280383B1 (en) Method for determining traffic information for a section of a road network and traffic calculator to implement the method
EP3611709A1 (en) Traffic flow simulator
EP0740280A2 (en) Disturbance detection method for road traffic
DE102012204306A1 (en) A method of controlling provision of traffic information data for updating traffic information
DE1574054B1 (en) DEVICE FOR DETERMINING CAR TRAFFIC FLOW
WO2013120765A1 (en) Method for model construction for a travel-time database
DE102017104274A1 (en) Method for determining a fault location in a ring-structured electrical network
EP2161698B1 (en) Method for coordinating light signal-controlled nodes in a street network
EP0501193B1 (en) Method for the automatic coordination of an independent intersection traffic light control device with one or more neighbouring intersections
DE102005056709A1 (en) Method and device for coupling at least two independent bus systems
DE102017008667A1 (en) Method for dynamic route planning
DE102016220561A1 (en) Method for determining driving routes for land vehicles
WO2012025188A1 (en) Method for predicting the electrical power an electrochemical energy store can output to a consumer
EP2945032B1 (en) Method for determining the switching times and/or the heating characteristic curve of a heating system
AT502459A1 (en) METHOD FOR DETERMINING TRANSPORT CONDITION
EP2887332B1 (en) Method and system for detection of a traffic situation on a stretch of road
DE102014203366A1 (en) Method and apparatus for determining the direction of current flowing through a circuit breaker
DE102016010028A1 (en) A technique for selecting a route in a traffic network by a car navigation device
DE102011087803B4 (en) Method and apparatus for merging classified traffic situation information
EP1056063A1 (en) Method and device for the determination of the traffic conditions of a roadsegment
DE102008012569B4 (en) Method and device for detecting events by means of at least two counters
DE102008045238A1 (en) Method for starting up a system for determining the energy flows in an energy distribution device
WO2004025836A1 (en) Method and device for detecting period length fluctuations of periodic signals

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20061213

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IS IT LI LT LU LV MC NL PL PT RO SE SI SK TR

AX Request for extension of the european patent

Extension state: AL BA HR MK YU

RIC1 Information provided on ipc code assigned before grant

Ipc: G08G 1/01 20060101AFI20070119BHEP

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

DAX Request for extension of the european patent (deleted)
GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IS IT LI LT LU LV MC NL PL PT RO SE SI SK TR

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

Free format text: NOT ENGLISH

GBT Gb: translation of ep patent filed (gb section 77(6)(a)/1977)

Effective date: 20080227

REG Reference to a national code

Ref country code: CH

Ref legal event code: NV

Representative=s name: SIEMENS SCHWEIZ AG

Ref country code: CH

Ref legal event code: EP

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

Free format text: LANGUAGE OF EP DOCUMENT: GERMAN

REF Corresponds to:

Ref document number: 502006000385

Country of ref document: DE

Date of ref document: 20080410

Kind code of ref document: P

REG Reference to a national code

Ref country code: GR

Ref legal event code: EP

Ref document number: 20080401067

Country of ref document: GR

REG Reference to a national code

Ref country code: SE

Ref legal event code: TRGR

REG Reference to a national code

Ref country code: DK

Ref legal event code: T3

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080627

Ref country code: ES

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080607

Ref country code: FI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080227

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: PL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080227

Ref country code: LV

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080227

Ref country code: SI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080227

REG Reference to a national code

Ref country code: IE

Ref legal event code: FD4D

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080227

Ref country code: SK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080227

Ref country code: PT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080721

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: RO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080227

EN Fr: translation not filed
PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MC

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20080531

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080227

26N No opposition filed

Effective date: 20081128

REG Reference to a national code

Ref country code: CH

Ref legal event code: PCAR

Free format text: SIEMENS SCHWEIZ AG;INTELLECTUAL PROPERTY FREILAGERSTRASSE 40;8047 ZUERICH (CH)

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080527

Ref country code: FR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20081219

Ref country code: EE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080227

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CY

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080227

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LU

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20080524

Ref country code: HU

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080828

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: TR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080227

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: BE

Payment date: 20100603

Year of fee payment: 5

GBPC Gb: european patent ceased through non-payment of renewal fee

Effective date: 20100524

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GB

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20100524

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: SE

Payment date: 20110510

Year of fee payment: 6

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: DK

Payment date: 20110512

Year of fee payment: 6

BERE Be: lapsed

Owner name: SIEMENS A.G.

Effective date: 20110531

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: BE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20110531

REG Reference to a national code

Ref country code: SE

Ref legal event code: EUG

REG Reference to a national code

Ref country code: DK

Ref legal event code: EBP

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20120525

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: DK

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20120531

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GR

Payment date: 20150511

Year of fee payment: 10

Ref country code: IT

Payment date: 20150514

Year of fee payment: 10

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: NL

Payment date: 20160519

Year of fee payment: 11

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: CZ

Payment date: 20160524

Year of fee payment: 11

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: AT

Payment date: 20160407

Year of fee payment: 11

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: DE

Payment date: 20160720

Year of fee payment: 11

Ref country code: CH

Payment date: 20160802

Year of fee payment: 11

REG Reference to a national code

Ref country code: GR

Ref legal event code: ML

Ref document number: 20080401067

Country of ref document: GR

Effective date: 20161207

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GR

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20161207

Ref country code: IT

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20160524

REG Reference to a national code

Ref country code: CH

Ref legal event code: PCOW

Free format text: NEW ADDRESS: WERNER-VON-SIEMENS-STRASSE 1, 80333 MUENCHEN (DE)

REG Reference to a national code

Ref country code: DE

Ref legal event code: R119

Ref document number: 502006000385

Country of ref document: DE

REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

REG Reference to a national code

Ref country code: NL

Ref legal event code: MM

Effective date: 20170601

REG Reference to a national code

Ref country code: AT

Ref legal event code: MM01

Ref document number: 387695

Country of ref document: AT

Kind code of ref document: T

Effective date: 20170524

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CZ

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20170524

Ref country code: AT

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20170524

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CH

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20170531

Ref country code: LI

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20170531

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: NL

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20170601

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: DE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20171201