JP2002170190A - Abnormality detection device in traffic stream and method therefor - Google Patents

Abnormality detection device in traffic stream and method therefor

Info

Publication number
JP2002170190A
JP2002170190A JP2001229852A JP2001229852A JP2002170190A JP 2002170190 A JP2002170190 A JP 2002170190A JP 2001229852 A JP2001229852 A JP 2001229852A JP 2001229852 A JP2001229852 A JP 2001229852A JP 2002170190 A JP2002170190 A JP 2002170190A
Authority
JP
Japan
Prior art keywords
vehicle
traffic flow
road
vehicles
abnormality
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
JP2001229852A
Other languages
Japanese (ja)
Other versions
JP3820935B2 (en
Inventor
Osamu Hattori
理 服部
Hiroshi Shimoura
弘 下浦
Kenji Tenmoku
健二 天目
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.)
Sumitomo Electric Industries Ltd
Original Assignee
Sumitomo Electric Industries Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sumitomo Electric Industries Ltd filed Critical Sumitomo Electric Industries Ltd
Priority to JP2001229852A priority Critical patent/JP3820935B2/en
Publication of JP2002170190A publication Critical patent/JP2002170190A/en
Application granted granted Critical
Publication of JP3820935B2 publication Critical patent/JP3820935B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Traffic Control Systems (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)

Abstract

PROBLEM TO BE SOLVED: To examine changes of the arrangement of a group of vehicles in one traffic lane (alignment) to detect the occurrence of an unexpected event on a road because, if an unexpected event such as an accident or a disaster occurs on the road, changes of traffic lanes are frequently performed, and as a result, a vehicle moves onto the other land and looks as if it disappears, another vehicle enters the lane from the other lane and it looks that the unexpected event occurs, and the order of vehicles is changed as a result of passing so that the alignment is changed, although the alignment is not changed unless the unexpected even occurs on the road. SOLUTION: An abnormality detection device in a traffic stream is provided with vehicle sensors 3, 5 installed at a plurality of vehicle observation points to observe a feature amount of passing vehicles, a matching decision part 22 deciding a degree of agreement of the alignment of a group of vehicles by corresponding vehicles passing respective vehicle observation points mutually, and an abnormality decision part 24 detecting the occurrence of the unexpected event on the road based on the results of judgement of the matching decision part 22.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、道路に車両感知器
等を設置して交通計測データを集め、この交通計測デー
タによって、突発事象の発生による交通流の異常を検知
することができる交通流の異常検知装置及び方法に関す
るものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention collects traffic measurement data by installing a vehicle sensor or the like on a road, and uses the traffic measurement data to detect a traffic flow abnormality caused by an unexpected event. And a method and apparatus for detecting an abnormality.

【0002】[0002]

【従来の技術】道路上に交通事故、災害などの突発事象
が発生したとき、この突発事象に基づく交通流の異常を
いち早く検知して、後続の車両に知らせたり、後続の車
両を誘導したりする必要がある。従来、道路にカメラを
設置して、画像処理をして交通流の異常を検知すること
が行われているが(特開平7−21488号公報、特開
平10−40490号公報など参照)、道路の広い範囲
にわたってカメラを設置することは経費がかかり、ま
た、夜間や悪天候時の検知が困難である。
2. Description of the Related Art When an unexpected event such as a traffic accident or disaster occurs on a road, an abnormality in traffic flow based on the incident is detected promptly to notify a following vehicle or guide a following vehicle. There is a need to. 2. Description of the Related Art Conventionally, a camera is installed on a road and image processing is performed to detect an abnormality in a traffic flow (see JP-A-7-21488, JP-A-10-40490, etc.). Installing cameras over a wide area is costly and difficult to detect at night or in bad weather.

【0003】そこで、道路の随所に設置した車両感知器
を使って道路の交通量、車両の速度などを測定し、これ
らの測定値に基づいて交通流の異常を監視することが行
われている。この監視装置によれば、交通量が少ないの
に速度が急激に落ち、その状態が一定時間継続したとき
に事故発生と判定している。
[0003] In view of this, it has been practiced to measure the traffic volume of a road, the speed of a vehicle, and the like using a vehicle sensor installed at various places on the road, and to monitor the abnormality of the traffic flow based on these measured values. . According to this monitoring device, it is determined that an accident has occurred when the speed drops sharply despite a small traffic volume and the state continues for a certain period of time.

【0004】[0004]

【発明が解決しようとする課題】ところが前記の監視装
置では、車両の走行速度に基づいて判定しているため、
自然渋滞時に突発事象が発生した場合に区別がつきにく
く、検知精度が低下するという問題があった。そこで、
発明者は、道路に突発事象が発生したときは、車線変更
が頻繁に行われることに着目した。ある車線に注目すれ
ば、車線変更の結果、車両が他の車線に移って消滅した
ように見えたり、他の車線から入ってきて発生したよう
に見えたり、追い越しの結果車両の順番が替わったりす
ることがある。
However, in the above-described monitoring device, the determination is made based on the traveling speed of the vehicle.
There is a problem that it is difficult to distinguish when an unexpected event occurs during a natural traffic jam, and the detection accuracy is reduced. Therefore,
The inventor has noted that lane changes are frequently performed when an unexpected event occurs on a road. If you focus on one lane, as a result of the lane change, the vehicle may appear to have disappeared by moving to another lane, may appear to be coming from another lane, or the order of vehicles may be changed as a result of passing May be.

【0005】一方、道路の上下流で、車両感知器やカメ
ラを設置することにより、車両を同定する技術が開発さ
れている。本発明は、道路を走行する複数の車両を車両
群としてとらえ、車両同定技術を利用して、車両群の中
の並びの変化に基づいて、道路上の突発事象の発生を検
知することができる交通流の異常検知装置及び方法を実
現することを目的とする。
On the other hand, there has been developed a technology for identifying a vehicle by installing a vehicle detector and a camera on the upper and lower sides of a road. The present invention regards a plurality of vehicles traveling on a road as a vehicle group, and can detect the occurrence of an unexpected event on the road based on a change in the arrangement in the vehicle group using the vehicle identification technology. It is an object of the present invention to implement a traffic flow abnormality detection device and method.

【0006】[0006]

【課題を解決するための手段】本発明の交通流の異常検
知装置は、複数の車両観測地点に設置され、通過車両の
特徴量を観測する車両特徴量観測手段と、それぞれの車
両観測地点を通過した車両同士を対応付けることによ
り、車両群の並び(アライメント)の一致度、すなわち
マッチング率を判定するマッチング判定手段と、マッチ
ング判定手段の判定結果に基づいて、道路上の突発事象
の発生を検知する異常判定手段とを有するものである
(請求項1)。
The traffic flow abnormality detecting device according to the present invention is installed at a plurality of vehicle observation points, and measures vehicle characteristic amount observing means for observing characteristic amounts of passing vehicles. By associating the vehicles that have passed with each other, a matching judgment unit for judging the degree of coincidence of the group of vehicles (alignment), that is, a matching ratio, and the occurrence of a sudden event on the road is detected based on the judgment result of the matching judgment unit. Abnormality determination means for performing the determination (claim 1).

【0007】道路上の突発事象がなければ、車両群のア
ライメントは、あまり変化しないが、道路に突発事象が
発生したときは、車線変更が頻繁に行われ、その結果、
車両が他の車線に移って消滅したように見えたり、他の
車線から入ってきて発生したように見えたり、追い越し
の結果車両の順番が替わったりするので、アライメント
は、変化する。そこで、車両の特徴量を観測し、複数の
車両観測地点を通過した車両同士を対応付けることによ
り、アライメントの一致度、すなわちマッチング率を調
べ、この一致度の低下に基づいて、道路上の突発事象の
発生を検知することができる。
[0007] If there is no sudden event on the road, the alignment of the group of vehicles does not change much. However, when a sudden event occurs on the road, lane changes are frequently performed.
The alignment changes because the vehicle appears to have disappeared as it moves to another lane, appears to have originated from another lane, or as a result of overtaking. Therefore, by observing the feature amount of the vehicle and associating the vehicles that have passed through the plurality of vehicle observation points with each other, the degree of coincidence of the alignment, that is, the matching rate, is examined. Can be detected.

【0008】一致度の低下は、しきい値との比較により
判定することができる(請求項2)。このしきい値は、
実際に道路上の突発事象が発生した時点前後の一致度を
記録しておき、本発明の実施により道路上の突発事象の
発生を最も精度よく検知することができるような値に選
べばよい。車両特徴量として、車長、車高、車両の画像
などを用いることができる(請求項3〜5)。
The decrease in the degree of coincidence can be determined by comparison with a threshold value. This threshold is
The degree of coincidence before and after the actual occurrence of the sudden event on the road is recorded, and the value may be selected so that the occurrence of the sudden event on the road can be detected with the highest accuracy by implementing the present invention. The vehicle length, the vehicle height, the image of the vehicle, and the like can be used as the vehicle feature amount (claims 3 to 5).

【0009】前記しきい値は、時間帯、曜日、催事のあ
る日ない日など、過去のいろいろな条件での検知データ
を蓄積し、時間帯や、曜日、催事の有無などに応じて、
最適なしきい値を設定するようにしてもよい(請求項
6)。前記しきい値は、交通状態に応じて自動的に決定
され、記憶されている値であってもよい(請求項7)。
本発明の交通流の異常検知装置では、道路上の突発事象
が発生していない正常時、道路が渋滞していない時はマ
ッチング率は自然に低下し、道路が渋滞している時はマ
ッチング率は自然に上昇する。したがって、道路が渋滞
していない時には前記しきい値を自動的に低下させて誤
検知率を下げ、道路が渋滞している時には前記しきい値
を自動的に上昇させて検知漏れを減らすこととする。
The threshold value is obtained by accumulating detection data under various conditions in the past, such as a time zone, a day of the week, and a day without an event.
An optimum threshold value may be set (claim 6). The threshold value may be a value that is automatically determined according to a traffic condition and stored.
In the traffic flow abnormality detection device of the present invention, the matching rate naturally decreases when no sudden event on the road occurs, when the road is not congested, and when the road is congested, the matching rate decreases. Rises naturally. Therefore, when the road is not congested, the threshold value is automatically lowered to reduce the false detection rate, and when the road is congested, the threshold value is automatically raised to reduce detection omission. I do.

【0010】また、本発明の交通流の異常検知装置は、
道路上の突発事象が検知された場合に、その突発事象の
発生を外部に知らせる情報提供手段をさらに有していて
もよい(請求項8)。例えば路側ビーコン、移動無線電
話、放送などの通信手段を用いて、車両のドライバに情
報を提供することができる。前記情報提供手段は、当該
地域内若しくは当該区間内を走行する車両、又は当該地
域内若しくは当該区間内への走行が予想される車両に対
して情報を提供してもよい(請求項9)。情報を提供す
る車両を、一定の関係する車両に限定する趣旨である。
[0010] The traffic flow abnormality detecting device of the present invention comprises:
When a sudden event on the road is detected, the information processing device may further include information providing means for notifying the outside of the occurrence of the sudden event (claim 8). For example, information can be provided to the driver of the vehicle using communication means such as a roadside beacon, a mobile radio telephone, and broadcasting. The information providing means may provide information to a vehicle traveling in the area or the section or a vehicle expected to travel in the area or the section (claim 9). The purpose is to limit the vehicles that provide the information to certain related vehicles.

【0011】情報提供手段は、すでに予定されている道
路上の事象に対しては、その事象の発生を検知しても外
部に知らせなくてもよい(請求項10)。例えば、道路
工事等のためある時間から車線が制限されることが分か
っているときは、その時間に当該工事区間の道路上に事
象を検知しても、外部に知らせない。これは外部に知ら
せることによる混乱を防止するためである。前記交通流
監視手段は、判定の基礎となった評価値の大きさ(尤
度)に応じて、段階的な判定を行い、前記情報提供手段
は、この交通流監視手段による段階的な判定の結果によ
って異常情報の内容を変えることが好ましい(請求項1
1)。突発事象発生の尤度(確実性)に応じて、例えば
「この先事故・止まれ」、「前方注意」など情報提供の
内容を変えることにより、ドライバなどに、より適切な
情報を与えることができる。
The information providing means does not have to detect the occurrence of the event on the road already scheduled or notify it to the outside (claim 10). For example, when it is known that the lane is restricted from a certain time due to road construction or the like, even if an event is detected on the road in the construction section at that time, it is not notified to the outside. This is to prevent confusion by notifying the outside. The traffic flow monitoring means makes a stepwise determination according to the magnitude (likelihood) of the evaluation value on which the determination is based, and the information providing means makes a stepwise determination by the traffic flow monitoring means. It is preferable to change the content of the abnormality information depending on the result.
1). By changing the content of information provision such as “accident / stop ahead” or “caution ahead” according to the likelihood (certainty) of the occurrence of the sudden event, more appropriate information can be given to the driver or the like.

【0012】前記交通流監視手段は、複数の道路区間で
交通流の異常を検知した場合、判定の基礎となった評価
値の大きさに応じて、異常発生区間を特定してもよい
(請求項12)。突発事象発生の尤度(確実性)の一番
高い区間を異常発生区間とすることにより、後続のドラ
イバなどに発生区間の情報や回避ルートの情報を知らせ
ることができる。また、本発明の交通流の異常検知方法
(請求項13)は、請求項1記載の交通流の異常検知装
置と同一の発明に係る方法である。
[0012] When the traffic flow monitoring means detects a traffic flow abnormality in a plurality of road sections, the traffic flow monitoring means may specify the abnormality occurrence section in accordance with the magnitude of the evaluation value serving as the basis of the determination. Item 12). By setting the section having the highest likelihood (certainty) of the occurrence of the sudden event as the abnormality occurrence section, it is possible to inform the subsequent driver or the like of the information of the occurrence section and the information of the avoidance route. The traffic flow abnormality detection method according to the present invention (claim 13) is the same method as the traffic flow abnormality detection device according to claim 1.

【0013】[0013]

【発明の実施の形態】以下、高速道路の交通流監視を例
にとって、本発明の実施の形態を、添付図面を参照しな
がら詳細に説明する。図1は、交通流の異常検知をする
ための交通流監視システムを示す概略図である。2車線
の高速道路1に、ループ埋め込み式の車両感知器5が、
間隔をおいて車線ごとに設置されている。また、車両の
上方から車高を測定する超音波式の車両感知器3も、車
線ごとに設置されている。これらの車両感知器3,5
は、一次処理装置4につながれている。一次処理装置4
は、各車両の速度、車高の検知を行う。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS The embodiments of the present invention will be described below in detail with reference to the accompanying drawings, taking monitoring of a traffic flow on a highway as an example. FIG. 1 is a schematic diagram showing a traffic flow monitoring system for detecting a traffic flow abnormality. On a two-lane highway 1, a vehicle sensor 5 embedded in a loop
It is installed in each lane at intervals. An ultrasonic vehicle sensor 3 for measuring the vehicle height from above the vehicle is also provided for each lane. These vehicle sensors 3,5
Are connected to the primary processing unit 4. Primary processing unit 4
Detects the speed and height of each vehicle.

【0014】また、高速道路1には、事故情報や路面情
報などを車両に知らせるための可変表示板6が設けられ
ている。また、車両と双方向通信を行う路側ビーコン7
が設けられている。さらに、高速道路1に接続する一般
道路2には、高速道路1の事故情報や路面情報などを、
高速道路1に入ろうとする車両に知らせるための可変表
示板9が設けられている。
The highway 1 is provided with a variable display panel 6 for notifying the vehicle of accident information and road surface information. In addition, a roadside beacon 7 that performs bidirectional communication with a vehicle.
Is provided. Further, on the general road 2 connected to the highway 1, accident information and road surface information of the highway 1 are provided.
A variable display board 9 is provided to notify a vehicle entering the highway 1.

【0015】交通管理センター10内部のコンピュータ
は、各区間に設置された一次処理装置4、路側ビーコン
7、可変表示板6などと、有線通信網12(無線通信網
であってもよい)を通して接続されている。また、コン
ピュータは、建設省、警察庁、消防庁などの関係機関1
3と通信回線を通して結ばれており、放送局14とも通
信回線を通して結ばれている。図2は、交通管理センタ
ー10のコンピュータ11の機能ブロック図である。
The computer in the traffic management center 10 is connected to the primary processing unit 4, the roadside beacon 7, the variable display panel 6, etc. installed in each section through a wired communication network 12 (which may be a wireless communication network). Have been. Computers are also used by related organizations such as the Ministry of Construction, the National Police Agency, and the Fire Service.
3 via a communication line, and is also connected to the broadcasting station 14 via a communication line. FIG. 2 is a functional block diagram of the computer 11 of the traffic management center 10.

【0016】コンピュータ11の入力処理部21には、
車両感知器3,5の感知信号がマッチング判定部22を
介して入力される。入力処理部21は、車両感知器3の
出力に基づいて各車両の最大車高を算出するとともに、
車両感知器5の2つのループの出力時間差に基づいて車
両の速度を測定し、これと車両の感知時間とに基づいて
車長を算出する。車両が通過するごとに車高、車長が算
出されるので、1又は複数の車高、車長のデータ列が車
線ごとに得られる。
The input processing unit 21 of the computer 11 includes:
The sensing signals of the vehicle sensors 3 and 5 are input via the matching determination unit 22. The input processing unit 21 calculates the maximum vehicle height of each vehicle based on the output of the vehicle sensor 3,
The speed of the vehicle is measured based on the output time difference between the two loops of the vehicle sensor 5, and the vehicle length is calculated based on the measured speed and the vehicle sensing time. Since the vehicle height and vehicle length are calculated each time the vehicle passes, one or a plurality of vehicle height and vehicle length data strings are obtained for each lane.

【0017】コンピュータ11のマッチング判定部22
の行う車両群マッチング処理方法をフローチャート(図
3)を用いて説明する。なお、車両群マッチング処理方
法には、以下に説明する方法以外にも、種々の方法が公
知であるが、以下では、その一例として、車両の対応付
けをp次元アライメント問題に帰着させて行う方法を説
明する(小林他「2地点車両観測情報からの全域的交通
流解析アルゴリズム」(社)情報処理学会 第72回ア
ルゴリズム研究会 平成12年3月21日発表)。
The matching determining section 22 of the computer 11
Will be described with reference to a flowchart (FIG. 3). Various methods are known as the vehicle group matching processing method in addition to the method described below. In the following, as an example, a method of associating vehicles with a p-dimensional alignment problem is performed. (Kobayashi et al., "Global Traffic Flow Analysis Algorithm from Two-Point Vehicle Observation Information" (published by Information Processing Society of Japan 72nd Algorithm Study Group, March 21, 2000)).

【0018】以下、1つの車線のみに着目する。他の車
線から当該車線に進入した車両や、当該車線から他の車
線に退出する車両は、交差点から当該車線に進入した車
両や、当該車線から交差点に退出した車両と同様に扱
う。走行方向上流の車両感知器3,5で感知した車高、
車長データ列を入力し(ステップS1)、走行方向下流の
車両感知器3,5で感知した車高、車長データ列を入力
する(ステップS2)。
Hereinafter, attention will be focused on only one lane. A vehicle that has entered the lane from another lane or a vehicle that has exited from the lane to another lane is treated in the same manner as a vehicle that has entered the lane from an intersection or exited from the lane to an intersection. The vehicle height detected by the vehicle detectors 3, 5 upstream of the traveling direction,
The vehicle length data string is input (step S1), and the vehicle height and vehicle length data strings detected by the vehicle sensors 3, 5 downstream in the traveling direction are input (step S2).

【0019】上流の車両感知器で感知した処理対象車両
数をnとし、車高データ列をhA1,hA2,hA3,‥‥,
hAnで表し(代表するときは「hA」で表す。)、車長
データ列をlA1,lA2,lA3,‥‥,lAn(代表すると
きは「lA」で表す。)で表す。下流の車両感知器で感
知した処理対象車両数をmとし、車高データ列をhB1,
hB2,hB3,‥‥,hBm(代表するときは「hB」で表
す。)で表し、車長データ列をlB1,lB2,lB3,‥
‥,lBm(代表するときは「lB」で表す。)で表す。
The number of vehicles to be processed detected by the upstream vehicle detector is defined as n, and the vehicle height data sequence is defined as hA1, hA2, hA3,.
The vehicle length data sequence is represented by lA1, lA2, lA3,..., lAn (representative is represented by "lA"). The number of vehicles to be processed detected by the downstream vehicle sensor is m, and the vehicle height data sequence is hB1,
hB2, hB3,..., hBm (represented by “hB” when represented), and the vehicle length data sequence is represented by lB1, lB2, lB3,.
‥, IBm (representatively represented by IB).

【0020】上流の車両感知器で車高hAと検出された
車高hの車両が、下流の車両感知器で車高hBと観測さ
れる確率は、車両感知器の観測誤差がガウス分布に従う
と仮定すれば、下記(1)式で表される。
The probability that a vehicle having the vehicle height h detected by the upstream vehicle sensor as the vehicle height hA is observed as the vehicle height hB by the downstream vehicle sensor is based on the observation error of the vehicle sensor following a Gaussian distribution. If it is assumed, it is expressed by the following equation (1).

【0021】[0021]

【数1】 (Equation 1)

【0022】ここでσ2h,Aは上流の車両感知器で観測さ
れる車高の分散、σ2h,Bは上流の車両感知器で観測され
る車高の分散である。同様に、上流の車両感知器で車長
lAと検出された車長lの車両が、下流の車両感知器で
車長lBと観測される確率は、下記(2)式で表される。
Here, σ 2 h, A is the variance of the vehicle height observed by the upstream vehicle sensor, and σ 2 h, B is the variance of the vehicle height observed by the upstream vehicle sensor. Similarly, the probability that the vehicle having the vehicle length l detected by the upstream vehicle detector as the vehicle length lA is observed as the vehicle length IB by the downstream vehicle detector is expressed by the following equation (2).

【0023】[0023]

【数2】 (Equation 2)

【0024】ここでσ2l,Aは上流の車両感知器で観測さ
れる車長の分散、σ2l,Bは上流の車両感知器で観測され
る車長の分散である。 (a)コスト評価基準として車高又は車長の差を用いる場
合 上流の車両感知器で車高hAと検出された車両が、下流
の車両感知器で車高hBと検出される事後確率密度Ph
(ha,hb)は、下記(3)式で表される。
Here, σ 2 l, A is the variance of the vehicle length observed by the upstream vehicle sensor, and σ 2 l, B is the variance of the vehicle length observed by the upstream vehicle sensor. (a) When the difference in vehicle height or vehicle length is used as a cost evaluation criterion The posterior probability density Ph in which a vehicle detected as a vehicle height hA by an upstream vehicle sensor is detected as a vehicle height hB by a downstream vehicle sensor
(ha, hb) is represented by the following equation (3).

【0025】[0025]

【数3】 (Equation 3)

【0026】同様に、上流の車両感知器で車長lAと検
出された車両が、下流の車両感知器で車長lBと検出さ
れる事後確率密度Pl(la,lb)は、下記(4)式で表され
る。
Similarly, the posterior probability density Pl (la, lb) that the vehicle detected as the vehicle length IA by the upstream vehicle sensor is detected as the vehicle length IB by the downstream vehicle sensor is represented by the following (4). It is expressed by an equation.

【0027】[0027]

【数4】 (Equation 4)

【0028】σ2h,A+σ2h,B=σ2hと表し、σ2l,A+σ
2l,B=σ2lと書くと、前記(3)(4)式は、それぞれ、(5)
(6)式のように書き換えられる。
Σ 2 h, A + σ 2 h, B = σ 2 h, σ 2 l, A + σ
2 l, B = Writing and sigma 2 l, wherein (3) (4), respectively, (5)
It can be rewritten as in equation (6).

【0029】[0029]

【数5】 (Equation 5)

【0030】[0030]

【数6】 (Equation 6)

【0031】(b) コスト評価基準として車高又は車長の
比を用いる場合 前記(5)(6)式に代えて、次の(7)(8)式を用いる。
(B) When the ratio of vehicle height or vehicle length is used as a cost evaluation criterion The following expressions (7) and (8) are used instead of the expressions (5) and (6).

【0032】[0032]

【数7】 (Equation 7)

【0033】[0033]

【数8】 (Equation 8)

【0034】以上に掲げた確率密度の自然対数を「マッ
チコスト」という。車高のみを用いる場合の、上流の車
両感知器で車高hAと検出された車両と、下流の車両感
知器で車高hBと検出された車両とのマッチコストd(h
a,hb)は、 d(ha,hb)=ln[Ph(ha,hb)] (9) となる。車長のみを用いる場合は、上流の車両感知器で
車長lAと検出された車両と、下流の車両感知器で車長
lBと検出された車両とのマッチコストd(la,lb)は、 d(la,lb)=ln[Pl(la,lb)] (10) となる。
The natural logarithm of the probability density listed above is called "match cost". When only the vehicle height is used, the match cost d (h) between the vehicle detected as the vehicle height hA by the upstream vehicle sensor and the vehicle detected as the vehicle height hB by the downstream vehicle sensor.
a, hb) becomes d (ha, hb) = ln [Ph (ha, hb)] (9) When only the vehicle length is used, the match cost d (la, lb) between the vehicle detected as the vehicle length IA by the upstream vehicle detector and the vehicle detected as the vehicle length IB by the downstream vehicle detector is: d (la, lb) = ln [Pl (la, lb)] (10)

【0035】車高と車長を併用する場合は、上流の車両
感知器で車高hA、車長lAと検出された車両と、下流の
車両感知器で車高hB、車長lBと検出された車両とのマ
ッチコストd(ha,la;hb,lb)は、 d(ha,la;hb,lb) =ln[Ph(ha,hb)]+ln[Pl(la,lb)] (11) となる。
When both the vehicle height and the vehicle length are used, the vehicle height hA and the vehicle length IA are detected by the upstream vehicle sensor, and the vehicle height hB and the vehicle length IB are detected by the downstream vehicle sensor. The match cost d (ha, la; hb, lb) for the vehicle is d (ha, la; hb, lb) = ln [Ph (ha, hb)] + ln [Pl (la, lb)] (11) Becomes

【0036】マッチング判定部22は、前記(9)(10)(1
1)式のいずれかによって、上流の車両感知器で感知した
車両、下流の車両感知器で感知した車両の全組み合わせ
について、マッチコストを算出する(ステップS3)。次
に、マッチング判定部22は、「ギャップコスト」を取
得する。このギャップコストは、車両感知器の特性に応
じた定数として記憶されているものである。マッチコス
トとギャップコストとを総称して「スコア」という。
The matching judgment section 22 performs the operations (9), (10), (1)
The match cost is calculated for all combinations of the vehicle sensed by the upstream vehicle sensor and the vehicle sensed by the downstream vehicle sensor by one of the formulas (Step S3). Next, the matching determination unit 22 acquires “gap cost”. This gap cost is stored as a constant corresponding to the characteristics of the vehicle sensor. The match cost and the gap cost are collectively referred to as “score”.

【0037】ギャップコストには、インターナルギャッ
プコスト(IGC)とイクスターナルギャップコスト
(EXGC)との2種類がある。IGCは、両地点の間
で追い越しが発生してその順序が入れ替わってしまい、
対応が付かないような場合に設定される。これは、車
高、車長のいずれかが3σ以上離れることはなく、3σ
以上離れるとそれは別の車両とみなすとの仮定に基づ
く。
There are two types of gap costs, an internal gap cost (IGC) and an external gap cost (EXGC). In the IGC, an overtaking occurs between the two points and the order is changed,
This is set when no response can be made. This means that either the vehicle height or the vehicle length will not be more than 3σ apart, and 3σ
It is based on the assumption that if you leave more than that, you consider it another vehicle.

【0038】車高のみを用いる場合、IGCの具体的な
値を示すと、(12)式のようになる。 d(ha,−)=d(−,hb)=(1/2)ln[Ph(0,3σh)] =−(1/2)ln((2π)1/2σh)−9/4 (12) 式中1/2を用いているのは、ギャップコストを用いる
場合は2本の枝を使うことになるので、枝1本分のコス
トにするためである(以下同じ)。車長のみを用いる場
合、IGCの具体的な値を示すと、(13)式のようにな
る。
When only the vehicle height is used, the specific value of IGC is expressed by the following equation (12). d (ha,-) = d (-, hb) = (1/2) ln [Ph (0,3.sigma.h)] =-(1/2) ln ((2.pi.) 1 / 2.sigma.h) -9/4 ( 12) The reason why 1/2 is used in the equation is that the cost is equivalent to one branch because two branches are used when the gap cost is used (the same applies hereinafter). When only the vehicle length is used, the specific value of the IGC is expressed by Expression (13).

【0039】 d(la,−)=d(−,lb)=(1/2)ln[Pl(0,3σl)] =−(1/2)ln((2π)1/2σl)−9/4 (13) 車高と車長を併用する場合は、IGCの具体的な値を示
すと、(14)式のようになる。 d(ha,−;−,−)=d(−,la;−,−) =d(−, −; hb,−)=d(−, −;−,lb) =(1/2){ln[Ph(0,3σh)]+ln[Pl(0,0)]} =(1/2){ln[Ph(0,0)]+ln[Pl(0, 3σl)]} =−(1/2)ln(2πσhσl)−9/4 (14) EXGCは、上流地点を通過した車両がまだ下流地点を
通過せず、車両の対応が付けられない場合に設定する。
このEXGCの値は、実際に、本装置設置後の初期設定
等の時点で、目視などで最適なマッチングが得られてい
ると確認された場合に、アライメントの最大スコアを算
出し、その最大スコアに基づいて決定されるものであ
る。
D (la, −) = d (−, lb) = (1/2) ln [Pl (0,3σl)] = − (1/2) ln ((2π) 1 / 2σl) -9 / 4 (13) When both the vehicle height and the vehicle length are used, the specific value of the IGC is expressed by the equation (14). d (ha, −; −, −) = d (−, la; −, −) = d (−, −; hb, −) = d (−, −; −, lb) = (1/2) { ln [Ph (0,3σh)] + ln [Pl (0,0)]} = (1/2) {ln [Ph (0,0)] + ln [Pl (0,3σl)]} = − (1 / 2) ln (2πσhσl) -9/4 (14) EXGC is set when a vehicle that has passed the upstream point has not yet passed the downstream point and cannot be associated with a vehicle.
The value of this EXGC is calculated by calculating the maximum score of the alignment when it is confirmed by visual observation or the like that the optimal matching is actually obtained at the time of initial setting after installation of the apparatus. Is determined based on

【0040】さらに、現実には対応し得ないと思われる
車両間の対応スコアは−∞にしておくことも考えられる
(例えば旅行時間が負になる、経験上推定される旅行時
間と比べると、あまりにも短いあるいは長い、など)。
以上の各コストが得られると、マッチコスト又はギャッ
プコストの和が最大となるように、上流の車両感知器で
感知したn台の車両と、下流の車両感知器で感知したm
台の車両との対応付けを行う。このため、2次元の文字
列アライメント問題として定式化する。
Furthermore, it is conceivable that the correspondence score between vehicles that cannot be handled in reality is set to -∞ (for example, when the travel time becomes negative, as compared with the travel time estimated from experience, Too short or long, etc.).
When the above costs are obtained, the n vehicles detected by the upstream vehicle detector and the m vehicles detected by the downstream vehicle detector are set so that the sum of the match costs or the gap costs is maximized.
It associates with one vehicle. For this reason, it is formulated as a two-dimensional character string alignment problem.

【0041】上流地点を通過した車両をa1,a2,a
3,‥‥,anで表し、下流地点を通過した車両をb1,
b2,b3,‥‥,bnで表す。行列(aibj)(1<i<n,1<
j<m)をアライメントと呼ぶ。2次元の文字列アライメン
ト問題は、2次元格子状有向グラフ上の最長路問題に帰
着できる(一般的な呼び方は「最短路問題」であるが、
ここではコストの和の最長のパスを求めているので、
「最長路問題」という)。図4は、2次元格子状有向グ
ラフを描いた図である。図4では、車両aiと車両bjと
を対応させる斜めの枝が実線で表されている。この枝長
は、前述したマッチコストd((9)〜(11)式のいずれ
か)に相当する。枠の内側の縦横の破線枝は、車両ai
と車両bjとの対応が付かない場合を表し、その枝長
は、前述したIGCに相当する。枠の外側の一点鎖線の
枝は、対応する車両がない場合を表し、その枝長は、前
述したEXGCに相当する。
The vehicles passing the upstream point are referred to as a1, a2, a
The vehicle passing through the downstream point is denoted by b1,
Expressed as b2, b3,..., bn. Matrix (aibj) (1 <i <n, 1 <
j <m) is called alignment. The two-dimensional character string alignment problem can be reduced to the longest path problem on a two-dimensional grid-like directed graph (the general term is "shortest path problem,
Here we are looking for the longest path of the sum of the costs,
"The longest road problem"). FIG. 4 is a diagram illustrating a two-dimensional grid-like directed graph. In FIG. 4, diagonal branches that correspond to the vehicle ai and the vehicle bj are represented by solid lines. This branch length corresponds to the above-mentioned match cost d (any one of equations (9) to (11)). The vertical and horizontal broken lines inside the frame indicate the vehicle ai
And the vehicle bj do not correspond, and the branch length corresponds to the above-mentioned IGC. The dashed-line branch outside the frame indicates a case where there is no corresponding vehicle, and the branch length corresponds to the above-described EXGC.

【0042】図4の左上の点から、右下の点に至るコス
トの和の最長のパスが、車両のもっともらしい対応付け
を示す解となる。この最長路問題は、動的計画(DP;D
ynamic Programming)法で解くことができる(ステップ
S4)(下記[1][2]参照)。 [1]D. Gusfield. "Aigorithms on Strings, Trees,
and Sequences." Cambridge University Prass, 1997. [2]S.B.Needleman and C.D.Wunsch."A general meth
od applicable to thesearch for similarities in the
amino acid sequences of two proteins." Journal of
Molecular Biology,48, pp.443-453, 1970. 例えば、n=3,m=4とし、最長路問題を解いた結
果、図5に示すような経路(太い実線)が得られたとす
る。この図5から、車両a1はb1に対応し、車両a2と
車両b2は対応せず(インターナルギャップ)、車両a3
はb3に対応し、車両b4に対応するものがない(イクス
ターナルギャップ)、ことが分かる。この原因は、車両
a2とa3は上流地点を通過した後入れ替わった、と考え
られ、車両b4は途中の交差点から進入してきた若しく
は車線間移動により入ってきた、と考えられる。
The longest path of the sum of the costs from the upper left point to the lower right point in FIG. 4 is the solution indicating the plausible association of the vehicle. The longest path problem is based on dynamic programming (DP; D
dynamic programming) (step S4) (see [1] and [2] below). [1] D. Gusfield. "Aigorithms on Strings, Trees,
and Sequences. "Cambridge University Prass, 1997. [2] SBNeedleman and CDWunsch." A general meth
od applicable to thesearch for similarities in the
amino acid sequences of two proteins. "Journal of
Molecular Biology, 48, pp. 443-453, 1970. For example, it is assumed that n = 3 and m = 4, and a route (thick solid line) as shown in FIG. 5 is obtained as a result of solving the longest path problem. From FIG. 5, vehicle a1 corresponds to b1, vehicle a2 does not correspond to vehicle b2 (internal gap), and vehicle a3
Corresponds to b3 and there is no corresponding to vehicle b4 (external gap). It is considered that this is because vehicles a2 and a3 have been switched after passing the upstream point, and vehicle b4 has entered from an intersection on the way or has entered due to movement between lanes.

【0043】以上の解析結果を出力する(ステップS
5)。そして、出力された結果に基づいて、一致する車
両同士に着目して、旅行時間を推定することができる。
ここで、EXGCの値を決定する方法を説明する。EX
GCの値をgとおく。上流でn台車両が観測され、下流
でm台の車両が観測されたとする。ただしn<mとす
る。目視などにより、上流地点で観測された車両a1〜
anは、下流地点において観測された車両b1〜bnに対
応していることが分かっているものとする。しかし、上
流地点で比較的遅い時刻に観測された車両an+1〜am
は、下流地点の観測時点において、まだ下流に到達して
おらず、対応するものがない。
The above analysis results are output (step S
5). Then, based on the output result, the travel time can be estimated by focusing on the matching vehicles.
Here, a method of determining the value of EXGC will be described. EX
Let the value of GC be g. It is assumed that n vehicles are observed upstream and m vehicles are observed downstream. Note that n <m. Vehicles a1-observed at the upstream point by visual inspection
An is assumed to correspond to the vehicles b1 to bn observed at the downstream point. However, vehicles an + 1 to am observed at a relatively late time at the upstream point
Has not reached downstream yet at the time of observation of the downstream point, and there is no corresponding one.

【0044】このときのアライメントの、マッチング部
分のスコアをS′(A)、非マッチング部分のスコアを
g・x(A)と書く。S′(A)は車両感知器の実測値
から算出される値である。gは求めたいEXGCの値で
あり、x(A)は非マッチング台数を示し、x(A)=
n−mである。EXGCの値gを「最適なマッチングが
得られている場合の、マッチング部分のスコアS′
(A)を、n+m−x(A)」で割ったもの、と定義す
る。
In this case, the score of the matching portion of the alignment is expressed as S '(A), and the score of the non-matching portion is expressed as g · x (A). S '(A) is a value calculated from the measured value of the vehicle sensor. g is the value of EXGC to be obtained, x (A) indicates the number of non-matching units, and x (A) =
nm. The value g of the EXGC is changed to “the score S ′ of the matching portion when the optimal matching is obtained.
(A) divided by "n + mx (A)".

【0045】 g=S′(A)/[n+m−x(A)] (15) たとえば、n=4,m=3とした場合の2次元格子状有
向グラフを描くと図6のようになる。図6において、マ
ッチング部分のスコアS′(A)、EXGCの値g・x
(A)を示している。ただし、x(A)=n−m=4−
3=1となるので、図6でg・x(A)と示したもの
は、gそのものを示している。この図6の例では、EX
GCの値gは、前記(15)に基づき、 g=S′(A)/[4+3−1]=S′(A)/6 (16) となる。
G = S ′ (A) / [n + mx (A)] (15) For example, if a two-dimensional lattice-like directed graph where n = 4 and m = 3 is drawn, it becomes as shown in FIG. In FIG. 6, the score S '(A) of the matching part and the value g · x of the EXGC
(A) is shown. Where x (A) = nm = 4-
Since 3 = 1, g.x (A) in FIG. 6 indicates g itself. In the example of FIG.
The value g of GC is g = S '(A) / [4 + 3-1] = S' (A) / 6 based on (15).

【0046】実際には、上下流でマッチングが確認され
ている車両群を観測してEXGCの値gを求める、とい
う処理を複数回行い、求められた複数のgの平均をと
り、この平均値を最終的にEXGCの値gとして決定す
ればよい。前記のIGCとEXGCとを使った実施形態
では、最新時点側(観測時刻の遅いほう)の車両対応付
け部分(mの最も大きいところのデータ)は、EXGC
の値に大きく依存してしまう。EXGCの値は、前述し
たように統計的に求められる値なので、変動する。
In practice, a process of observing a group of vehicles for which matching has been confirmed upstream and downstream to obtain the EXGC value g is performed a plurality of times, and the average of the obtained plurality of g is calculated. May be finally determined as the EXGC value g. In the embodiment using the IGC and the EXGC, the vehicle association portion (data with the largest m) on the latest time point (latest observation time) is EXGC
Greatly depends on the value of. The value of EXGC fluctuates because it is a value obtained statistically as described above.

【0047】そこで、時系列最新時点側のEXGCを必
要としないようにアルゴリズムを拡張する方法を説明す
る。前述したとおり、EXGCの値gは、 g=S′(A)/[n+m−x(A)] (15) で表される。一方、トータルスコアS(A,g)は、 S(A,g)=S′(A)+g・x(A) (17) で表される。(15)式を(17)式に代入すると、 S(A,g) =S′(A)+S′(A)x(A)/[n+m−x(A)] =(n+m)S′(A)/[n+m−x(A)] (18) となる。n+mは一定なので、S′(A)/[n+m−
x(A)]を最大にするようなパスが最適な解となる。
Therefore, a method of extending the algorithm so as not to require the EXGC at the latest time point in the time series will be described. As described above, the value g of EXGC is represented by g = S '(A) / [n + mx (A)] (15). On the other hand, the total score S (A, g) is represented by S (A, g) = S '(A) + g.x (A) (17) By substituting equation (15) into equation (17), S (A, g) = S '(A) + S' (A) x (A) / [n + mx (A)] = (n + m) S '( A) / [n + mx (A)] (18) Since n + m is constant, S '(A) / [n + m-
x (A)] is the optimal solution.

【0048】上に述べたことを図解すると、次のように
なる。図7は、時系列最新時点側のイクスターナルギャ
ップを除いた2次元格子状有向グラフを描いた図であ
る。この図7のグラフにおいて、上流側通過車両a1,
‥‥,anの添え字をi(1≦i≦n)とし、最新時点
側の対応点をviとする。各点viでのアライメントの最
大スコアを求め、そのスコアを(i+m)−x(A)で
割った値が最大になる点を最新時点側の対応点vi0とす
る。
The above is illustrated as follows. FIG. 7 is a diagram illustrating a two-dimensional lattice-like directed graph excluding the external gap at the latest time point in the time series. In the graph of FIG. 7, the upstream passing vehicles a1,
The subscripts of ‥‥ and an are i (1 ≦ i ≦ n), and the corresponding point on the latest time point is vi. The maximum score of the alignment at each point vi is obtained, and the point at which the value obtained by dividing the score by (i + m) -x (A) becomes the maximum is defined as the latest corresponding point vi 0 .

【0049】以上で、本発明の実施の形態を説明した
が、本発明の実施は、前記の形態に限定されるものでは
ない。例えば、本発明において、車長や車高のデータ以
外に、カメラで車両の画像データを取得してアライメン
トに利用することができる。この場合、画像間のスコア
を定める必要があるが、画像検索などの分野で考案され
ている「画像間の距離」、具体的には、マッチディスタ
ンス(下記[3]参照)、EMD(Earth Mover's Dista
nce)(下記[4]参照)を利用することができる。
Although the embodiment of the present invention has been described above, the embodiment of the present invention is not limited to the above embodiment. For example, in the present invention, in addition to vehicle length and vehicle height data, camera image data can be acquired by a camera and used for alignment. In this case, it is necessary to determine a score between images, but the “distance between images” devised in fields such as image search, specifically, match distance (see [3] below), EMD (Earth Mover's Dista
nce) (see [4] below).

【0050】[3]M.Werman, S.Peieg, and A.Rosenfe
ld."A distance metric for multi-dimensional histgr
ams." Computer, Vision, Graphics, and Image Proces
sing, 32, pp.328-336, 1985. [4]Y.Rubner, C Tomasi, and L.J.Guibas."The Eart
h Mover's Distance as a Metric for Image Retrieva
l. "Technical Report STAN-CS-TN-98-86, Department
of Computer Science, Stanford University, Septembe
r 1998. マッチディスタンスは、画像内の画素の色情報を利用し
た距離で、2つの画像ヒストグラムの累積ヒストグラム
間のL1距離として与えられる。具体的には、各車両の
後方より撮影した画像から、輝度のヒストグラムと色相
のヒストグラムを作る。輝度のヒストグラムは、すべて
の画素を輝度より32レベルに分け、度数は各レベルの
画素数とする。色相のヒストグラムは色相を30段階に
分け、度数は、各段階に含まれる画素の彩度の総和とす
る。画像の距離は、輝度と色相のヒストグラムによる距
離の加重和とした。
[3] M. Werman, S. Peieg, and A. Rosenfe
ld. "A distance metric for multi-dimensional histgr
ams. "Computer, Vision, Graphics, and Image Proces
sing, 32, pp.328-336, 1985. [4] Y. Rubner, C Tomasi, and LJ Guibas. "The Eart
h Mover's Distance as a Metric for Image Retrieva
l. "Technical Report STAN-CS-TN-98-86, Department
of Computer Science, Stanford University, Septembe
r 1998. The match distance is a distance using color information of a pixel in an image, and is given as an L1 distance between cumulative histograms of two image histograms. Specifically, a histogram of the luminance and a histogram of the hue are created from an image photographed from behind each vehicle. In the luminance histogram, all pixels are divided into 32 levels from the luminance, and the frequency is the number of pixels at each level. The hue histogram divides the hue into 30 levels, and the frequency is the sum of the saturations of the pixels included in each level. The image distance was a weighted sum of the distances based on the histogram of luminance and hue.

【0051】EMDは、マッチディスタンスを一般化し
た距離で、画素の色情報の他に位置情報も用いる。各画
像の画素を色と位置の近さによってクラスタリングし
(下記[5]参照)、それらの間の最小費用流として与
えられる。今回は、各画像を16クラスタにクラスタリ
ングした。 [5]M.Inaba, H.Imai, and N.Katoh."Application of
weighted Voronoi diagrams and randomization to va
riance-based k-clastering." In Proceedingsof the 1
0th ACM Symposium on Computational Geometry, pp.33
2-339, 1994車両aiと車両bjとの画像間の距離をDij
としたとき、アライメントのスコアDtotalとして、前
記(9)(10)(11)式のいずれかと、画像間の距離Dijとの
線形和を採用する。
The EMD is a generalized distance of the match distance, and uses position information in addition to pixel color information. The pixels of each image are clustered by color and proximity (see [5] below) and given as the least cost stream between them. This time, each image was clustered into 16 clusters. [5] M. Inaba, H. Imai, and N. Katoh. "Application of
weighted Voronoi diagrams and randomization to va
riance-based k-clastering. "In Proceedingsof the 1
0th ACM Symposium on Computational Geometry, pp.33
2-339, 1994 The distance between the image of the vehicle ai and the image of the vehicle bj is Dij
Then, as the alignment score Dtotal, a linear sum of any of the expressions (9), (10), and (11) and the distance Dij between images is adopted.

【0052】 Dtotal=d(ha,la;hb,lb) +λDij (18) ここで定数λは重み付け係数である。また、車長・車高
のギャップコストをGとし、実対応車両間の距離の最大
値をMとしたとき、総ギャップコストGtotalは、 Gtotal=G−λM/2 (19) として与えられる。
Dtotal = d (ha, la; hb, lb) + λDij (18) Here, the constant λ is a weighting coefficient. Further, assuming that the gap cost of the vehicle length / height is G and the maximum value of the distance between the corresponding vehicles is M, the total gap cost Gtotal is given as Gtotal = G−λM / 2 (19).

【0053】また、前記の実施の形態では、車両観測地
点は、2地点としていたが、これを任意の複数地点に拡
張することも可能である。車両観測地点がp(pは2以
上の整数)地点あれば、車両の対応付けをp次元アライ
メント問題に帰着させることが可能である。以上のよう
にして、上流地点を通過したn台の車両と、下流地点を
通過したm台の車両との対応付けが行われた。
In the above-described embodiment, the number of vehicle observation points is two. However, the number of vehicle observation points can be extended to arbitrary plural points. If the vehicle observation point is a point p (p is an integer of 2 or more), the association of the vehicle can be reduced to a p-dimensional alignment problem. As described above, the association between the n vehicles passing the upstream point and the m vehicles passing the downstream point was performed.

【0054】次に、異常判定部24の行う交通流の異常
検知処理をフローチャート(図8)を用いて説明する。
この異常検知処理は、上流地点を通過した車両と、下流
地点を通過した車両との対応付けをするのに同期して、
行う。異常判定部24は、上下流地点での対応付け結果
を入力する(ステップT1)。そして、下流地点を通過し
た車両のうち、何割が上流地点を通過した車両と対応付
けされているのかを調べる。この割合をマッチング率と
いう。例えば、下流地点を通過したm台の車両のうち、
p台が、上流地点を通過した車両と対応するとすれば、
マッチング率は、p/mとなる。このマッチング率をし
きい値と比較する。
Next, the traffic flow abnormality detection processing performed by the abnormality determination unit 24 will be described with reference to a flowchart (FIG. 8).
This abnormality detection processing, in synchronization with the vehicle that has passed the upstream point and the vehicle that has passed the downstream point,
Do. The abnormality determination unit 24 inputs the result of association at the upstream and downstream points (step T1). Then, of the vehicles that have passed the downstream point, it is checked what percentage is associated with the vehicle that has passed the upstream point. This ratio is called a matching ratio. For example, out of m vehicles that passed the downstream point,
Assuming that p vehicles correspond to vehicles that passed the upstream point,
The matching rate is p / m. This matching ratio is compared with a threshold.

【0055】このしきい値は、実際に道路上の突発事象
が発生した時点前後のマッチング率のデータを記録して
おき、本発明の実施により道路上の突発事象の発生を最
も精度よく検知することができるような値に選べばよ
い。さらに、時間帯、曜日、催事のある日ない日など、
過去のいろいろな条件での車線利用率のデータを蓄積
し、時間帯や、曜日、催事の有無などに応じて、最適な
しきい値を設定するようにしてもよい。
As the threshold value, data of the matching rate before and after the time when the sudden event on the road actually occurs is recorded, and the occurrence of the sudden event on the road is detected with the highest accuracy by implementing the present invention. Choose a value that can be used. In addition, such as time of day, day of the week, days with no events,
It is also possible to accumulate data on lane utilization rates under various conditions in the past, and set an optimal threshold value according to the time zone, the day of the week, the presence or absence of an event, and the like.

【0056】さらに、旅行時間を測定し、旅行時間が比
較的短い場合、つまり自由流に近い場合、しきい値を比
較的小さくし、旅行時間が比較的長い場合、つまり渋滞
の場合、しきい値を比較的大きくする(ステップT2〜
T4)。旅行時間に代えて、交通量、占有率に基づいて
渋滞かどうかを判断してもよい。自由流に近い場合、突
発事象が発生していない平常時でも、マッチング率は比
較的低めになる。したがって、しきい値を下げないと誤
検知する可能性が高くなる。一方、渋滞時は、平常時で
もマッチング率は比較的高くなる。したがって、しきい
値を上げないと、検知漏れが多くなる。
Further, the travel time is measured, and when the travel time is relatively short, that is, when it is close to free flow, the threshold value is made relatively small. When the travel time is relatively long, that is, when traffic is congested, the threshold is set. Make the value relatively large (step T2
T4). Instead of the travel time, it may be determined whether or not there is a traffic jam based on the traffic volume and the occupancy. When the flow is close to the free flow, the matching rate is relatively low even in a normal time when no sudden event occurs. Therefore, the possibility of erroneous detection increases unless the threshold is lowered. On the other hand, during a traffic jam, the matching rate is relatively high even in normal times. Therefore, if the threshold value is not increased, detection omission increases.

【0057】このステップT2〜T4におけるしきい値
の増減は、人為的に設定してもよいが、自動的に行うよ
うにすることが、省力のために望ましい。異常判定部2
4は、マッチング率としきい値とを比較し(ステップT
5)、しきい値よりも低ければ、上下流地点の間で突発
事象が発生したと判断する(ステップT7)。しきい値よ
り高ければ、さらに、過去(例えば1回前)の上下流地
点を通過した車両の対応付で得られたマッチング率がし
きい値よりも高かったかどうかを調べ(ステップT6)、
高ければ異常なしと判断する。低ければ、突発事象が続
いていると判断する。
The increase or decrease of the threshold value in steps T2 to T4 may be set artificially, but it is desirable to automatically perform the increase or decrease in order to save labor. Abnormality determination unit 2
4 compares the matching rate with the threshold (step T
5) If it is lower than the threshold value, it is determined that a sudden event has occurred between the upstream and downstream points (step T7). If it is higher than the threshold value, it is further checked whether or not the matching rate obtained by associating the vehicles that have passed the past (for example, one time before) upstream and downstream points was higher than the threshold value (step T6).
If it is higher, it is determined that there is no abnormality. If it is low, it is determined that the sudden event has continued.

【0058】このように、マッチング率が2回以上連続
してしきい値よりも高い場合に、異常なしと判断するの
は、突発事象が起こっても偶然マッチングした場合にお
ける誤った判断を避けるためである。以上に説明した突
発事象の発生が複数の区間で判定される場合がある。こ
の場合は、各区間におけるマッチング率を比較して、も
っともマッチング率の低い区間を突発事象発生区間とし
て特定することができる。
As described above, when the matching rate is higher than the threshold value twice or more consecutively, it is determined that there is no abnormality in order to avoid an erroneous determination in the case of a sudden match even if an unexpected event occurs. It is. The occurrence of the sudden event described above may be determined in a plurality of sections. In this case, by comparing the matching rates in the sections, the section having the lowest matching rate can be specified as the sudden event occurrence section.

【0059】図9は、突発事象発生区間を特定する処理
を説明するためのフローチャートである。まず、いまま
で説明した突発事象発生判定処理を、それぞれの監視対
象道路区間1,2,‥‥,i,‥‥で行う(ステップW
1)。すべての監視対象道路区間1,2,‥‥,iで同
処理が終了すれば(ステップW2のYES)、突発事象発生
と判定された区間があるかどうか調べる(ステップW
3)。そして、各区間で算出されたマッチング率を比較
する(ステップW4)。このマッチング率が最小を示す
区間を、突発事象発生区間と特定する(ステップW
5)。
FIG. 9 is a flowchart for explaining the process of specifying the sudden event occurrence section. First, the sudden event occurrence determination processing described so far is performed in each of the monitored road sections 1, 2,.
1). If the processing is completed in all the monitored road sections 1, 2,..., I (YES in step W2), it is checked whether or not there is a section determined to have an unexpected event (step W).
3). Then, the matching rates calculated in the respective sections are compared (step W4). The section in which the matching rate shows the minimum is specified as the sudden event occurrence section (step W
5).

【0060】以下、マッチング率の逆数を「評価値」と
いうことにする。すなわち、評価値=1/(マッチング
率)。図10は、5つの区間での検知処理結果から得ら
れた評価値の時間推移を示すグラフである。このグラフ
によれば、事故は8時20分に発生し、各区間1〜3で
の評価値が上がっている。特に区間2の評価値が最大で
あるので、区間2が突発事象発生区間と特定することが
できる。
Hereinafter, the reciprocal of the matching ratio is referred to as “evaluation value”. That is, the evaluation value = 1 / (matching rate). FIG. 10 is a graph showing the time transition of the evaluation values obtained from the detection processing results in five sections. According to this graph, the accident occurred at 8:20, and the evaluation values in the sections 1 to 3 have increased. In particular, since the evaluation value of the section 2 is the maximum, the section 2 can be specified as the sudden event occurrence section.

【0061】以上のようにして突発事象の発生及びその
発生区間が決定されると、交通管理センター10は、可
変表示板6,9に、突発事象の発生を表示し、路側ビー
コン7を通して車両に突発事象の発生を通知する。この
通知にあたっては、前記評価値を複数のしきい値と比較
する。1つは「検知しきい値」であり、他の1つは「注
意しきい値」である。「検知しきい値」>「注意しきい
値」の関係がある。
When the occurrence of the sudden event and the section where the sudden event has occurred are determined as described above, the traffic management center 10 displays the occurrence of the sudden event on the variable display boards 6 and 9, and transmits the information to the vehicle through the roadside beacon 7. Notifies the occurrence of a sudden event. In this notification, the evaluation value is compared with a plurality of thresholds. One is a “detection threshold” and the other is a “attention threshold”. There is a relationship of “detection threshold”> “attention threshold”.

【0062】評価値が検知しきい値を超えていれば突発
事象の発生尤度が十分に高く「突発事象発生」と判断す
る。この検知しきい値が高すぎると検知漏れが多くな
り、検知しきい値が低すぎると誤検知が増える。この検
知しきい値は、検知漏れ率や誤検知率の実績に基づき、
自動的に決定されるようにしてもよい。検知しきい値を
超えていなければ、この評価値を注意しきい値と比較す
る。注意しきい値を超えていれば、突発事象の発生尤度
が中程度に高く、「突発事象の発生の可能性が高い注意
状態」と判断する。
If the evaluation value exceeds the detection threshold value, the likelihood of occurrence of a sudden event is sufficiently high and it is determined that "the occurrence of a sudden event". If this detection threshold is too high, detection omissions will increase, and if the detection threshold is too low, false detections will increase. This detection threshold is based on the results of detection omission rate and false detection rate,
It may be determined automatically. If the detection threshold is not exceeded, this evaluation value is compared with the caution threshold. If the caution threshold value is exceeded, the likelihood of the occurrence of the sudden event is moderately high, and it is determined that "the attention state in which the possibility of occurrence of the sudden event is high".

【0063】注意しきい値を超えていなければ、突発事
象の発生尤度が低く、「突発事象の発生なし」と判断す
る。「突発事象発生」と判定されていれば、交通管理セ
ンター10の出力処理部25は、可変表示板6,9に
「この先事故・止まれ」のような運転者の警告を与える
メッセージを表示し、路側ビーコン7を通して車両にも
危険区間である旨を通知する。
If the caution threshold value is not exceeded, the likelihood of occurrence of a sudden event is low, and it is determined that "no sudden event has occurred". If it is determined that “an unexpected event has occurred”, the output processing unit 25 of the traffic management center 10 displays a message giving a driver warning such as “accident / stop ahead” on the variable display plates 6 and 9, The vehicle is also notified of the dangerous section through the roadside beacon 7.

【0064】「突発事象の発生の可能性が高い注意状
態」と判定されていれば、出力処理部25は、可変表示
板6,9に「前方注意」のように運転者の注意を喚起す
るようなメッセージを表示し、路側ビーコン7を通して
車両にも走行注意区間である旨を通知する。そして、こ
れらの情報を通信回線を通して関係機関13や放送局1
4に連絡する。
If it is determined that the caution state is highly likely to cause a sudden event, the output processing unit 25 prompts the driver to pay attention to the variable display plates 6 and 9 as "forward warning". Such a message is displayed, and the vehicle is notified through the roadside beacon 7 to the effect that the section is a traveling caution section. Then, the information is transmitted to the related organization 13 or the broadcasting station 1 through a communication line.
Contact 4

【0065】なお、すでに道路工事などが予定され、交
通流の異常が予想されている場合は、出力処理部25
は、当該時刻に突発事象の発生を判定しても、この判定
に基づいて可変表示板6,9に突発事象の発生を表示す
ることはなく、関係機関13や放送局14に連絡するこ
ともない。以上で本発明の実施の形態を説明したが、本
発明の実施は、前記に限られるものではない。前記の例
では、高速道路を想定していたが、一般道路であっても
よい。車線数が2車線の道路を想定したが、車線数は、
2に限られるものではなく、1車線であっても3以上の
車線であってもよい。
If road construction or the like has already been scheduled and traffic flow abnormalities are expected, the output processing unit 25
Does not display the occurrence of the catastrophic event on the variable display boards 6 and 9 based on this judgment even if the occurrence of the catastrophic event is determined at the time, and may contact the related organization 13 or the broadcasting station 14. Absent. Although the embodiment of the present invention has been described above, the embodiment of the present invention is not limited to the above. In the above example, a highway is assumed, but a general road may be used. Assuming a road with two lanes, the number of lanes is
The number of lanes is not limited to two, and may be one lane or three or more lanes.

【0066】また、複数埋め込み式の車両感知器5に代
えて、道路の脇に設置されるドップラー式の車両感知器
を用いてもよい。また、道路にテレビカメラを設置して
画像処理により車両通過台数、車高、車長、通過速度な
どを検知してもよい。また、マッチング率を算出する方
法として、他の公知の方法を採用てもよい。例えば、車
両の画像を撮影し、画像処理を用いて車両の特徴、例え
ば色やプレートナンバーを認識する方法を採用してもよ
い。
Further, a Doppler type vehicle sensor installed beside a road may be used in place of the plural embedded type vehicle sensor 5. Alternatively, a television camera may be installed on the road to detect the number of vehicles passing, vehicle height, vehicle length, passing speed, and the like by image processing. Further, as a method of calculating the matching ratio, another known method may be employed. For example, a method of taking an image of a vehicle and recognizing a feature of the vehicle, for example, a color or a plate number, using image processing may be adopted.

【0067】[0067]

【発明の効果】以上のように本発明の交通流の異常検知
装置によれば、道路上の突発事象の発生を精度よく検知
することができる。
As described above, according to the traffic flow abnormality detecting device of the present invention, it is possible to accurately detect the occurrence of a sudden event on a road.

【図面の簡単な説明】[Brief description of the drawings]

【図1】交通流の異常検知をするための交通流監視シス
テムを示す概略図である。
FIG. 1 is a schematic diagram showing a traffic flow monitoring system for detecting a traffic flow abnormality.

【図2】交通管理センター10のコンピュータ11の機
能ブロック図である。
FIG. 2 is a functional block diagram of a computer 11 of the traffic management center 10.

【図3】マッチング判定部22の行う車両群マッチング
処理方法を説明するためのフローチャートである。
FIG. 3 is a flowchart for explaining a vehicle group matching processing method performed by a matching determination unit 22;

【図4】2次元格子状有向グラフを描いた図である。FIG. 4 is a diagram illustrating a two-dimensional grid-like directed graph.

【図5】n=3,m=4とし、最長路問題を解いた結果
を示すための、2次元格子状有向グラフの図である。
FIG. 5 is a diagram of a two-dimensional lattice-like directed graph for showing the result of solving the longest path problem with n = 3 and m = 4.

【図6】n=4,m=3とし、車両b4のみマッチング
されていない場合の2次元格子状有向グラフである。
FIG. 6 is a two-dimensional grid-like directed graph when n = 4 and m = 3 and only vehicle b4 is not matched.

【図7】時系列最新時点側のイクスターナルギャップを
必要としないアルゴリズムを説明するための、2次元格
子状有向グラフを描いた図である。
FIG. 7 is a diagram illustrating a two-dimensional grid-like directed graph for explaining an algorithm that does not require an external gap at the latest time point in a time series.

【図8】異常判定部24が行う突発事象発生を監視する
処理を説明するためのフローチャートである。
FIG. 8 is a flowchart illustrating a process performed by the abnormality determination unit for monitoring occurrence of a sudden event.

【図9】突発事象発生区間を特定する処理を説明するた
めのフローチャートである。
FIG. 9 is a flowchart illustrating a process of specifying a sudden event occurrence section.

【図10】5つの区間での検知処理結果から得られた評
価値の時間推移を示すグラフである。
FIG. 10 is a graph showing a time transition of an evaluation value obtained from a detection processing result in five sections.

【符号の説明】[Explanation of symbols]

1 高速道路 2 一般道路 3 車両感知器 4 一次処理装置 5 車両感知器 6 可変表示板 7 路側ビーコン 9 可変表示板 10 交通管理センター 11 コンピュータ 13 関係機関 14 放送局 22 マッチング判定部 24 異常判定部 DESCRIPTION OF SYMBOLS 1 Expressway 2 General road 3 Vehicle detector 4 Primary processing unit 5 Vehicle detector 6 Variable display board 7 Roadside beacon 9 Variable display board 10 Traffic management center 11 Computer 13 Related organization 14 Broadcasting station 22 Matching determination unit 24 Abnormality determination unit

フロントページの続き (51)Int.Cl.7 識別記号 FI テーマコート゛(参考) G08G 1/04 G08G 1/04 D (72)発明者 天目 健二 大阪市此花区島屋一丁目1番3号 住友電 気工業株式会社大阪製作所内 Fターム(参考) 5B057 AA16 AA19 BA11 BA30 CH01 CH11 CH14 DA13 DC03 DC36 5H180 AA01 BB02 BB04 BB05 BB15 CC04 CC11 CC15 CC18 EE03 JJ03 JJ05 JJ27 5L096 BA02 BA04 CA02 FA06 FA70 FA76 GA51 GA53 HA08 HA09 JA03 JA11 Continued on the front page (51) Int.Cl. 7 Identification code FI Theme coat II (Reference) G08G 1/04 G08G 1/04 D (72) Inventor Kenji Tenme 1-3-1 Shimaya, Konohana-ku, Osaka Sumitomo Electric Energy F-term (reference) in Osaka Seisakusho Co., Ltd.

Claims (13)

【特許請求の範囲】[Claims] 【請求項1】複数の車両観測地点に設置され、通過車両
の特徴量を観測する車両特徴量観測手段と、 それぞれの車両観測地点を通過した車両同士を対応付け
ることにより、車両群の並びの一致度を判定するマッチ
ング判定手段と、 マッチング判定手段の判定結果に基づいて、道路上の突
発事象の発生を検知する異常判定手段とを有することを
特徴とする交通流の異常検知装置。
1. A vehicle feature amount observing means installed at a plurality of vehicle observing points and observing a feature amount of a passing vehicle, and the vehicles passing through each vehicle observing point are associated with each other, so that the arrangement of the vehicle groups is matched. A traffic flow abnormality detection device, comprising: a matching determination unit that determines a degree; and an abnormality determination unit that detects occurrence of a sudden event on a road based on a determination result of the matching determination unit.
【請求項2】異常判定手段は、車両群の並びの一致度が
しきい値よりも小さくなったときに、道路上の突発事象
の発生を検知するものであることを特徴とする請求項1
記載の交通流の異常検知装置。
2. The system according to claim 1, wherein said abnormality judging means detects occurrence of a sudden event on the road when the degree of coincidence of the group of vehicles becomes smaller than a threshold value.
Abnormality detection device for traffic flow described.
【請求項3】車両特徴量観測手段は、車長を観測するも
のである請求項1記載の交通流の異常検知装置。
3. The traffic flow abnormality detecting device according to claim 1, wherein the vehicle characteristic amount observing means observes a vehicle length.
【請求項4】車両特徴量観測手段は、車高を観測するも
のである請求項1記載の交通流の異常検知装置。
4. The traffic flow abnormality detecting device according to claim 1, wherein the vehicle characteristic amount observing means observes a vehicle height.
【請求項5】車両特徴量観測手段は、車両の画像を撮影
するものである請求項1記載の交通流の異常検知装置。
5. The traffic flow abnormality detecting device according to claim 1, wherein the vehicle characteristic amount observing means captures an image of the vehicle.
【請求項6】前記しきい値は、時間帯、曜日、催事の有
無などに応じて統計的に求められ、記憶されている値で
ある請求項2記載の交通流の異常検知装置。
6. The traffic flow abnormality detecting device according to claim 2, wherein the threshold value is a value statistically obtained and stored according to a time zone, a day of the week, the presence or absence of an event, and the like.
【請求項7】前記しきい値は、交通状態に応じて自動的
に決定され、記憶されている値である請求項2記載の交
通流の異常検知装置。
7. The traffic flow abnormality detecting device according to claim 2, wherein said threshold value is a value automatically determined and stored according to a traffic condition.
【請求項8】交通流監視手段により道路上の突発事象が
検知された場合に、その突発事象の発生を外部に知らせ
る情報提供手段をさらに有することを特徴とする請求項
1の交通流の異常検知装置。
8. The traffic flow abnormality according to claim 1, further comprising information providing means for notifying the occurrence of the sudden event to the outside when the sudden event on the road is detected by the traffic flow monitoring means. Detection device.
【請求項9】情報提供手段は、当該地域内若しくは当該
区間内を走行する車両、又は当該地域内若しくは当該区
間内への走行が予想される車両に対して情報を提供する
ものであることを特徴とする請求項8記載の交通流の異
常検知装置。
9. The information providing means provides information to a vehicle traveling in the area or the section or a vehicle expected to travel in the area or the section. 9. The traffic flow abnormality detection device according to claim 8, wherein:
【請求項10】情報提供手段は、すでに予定されている
道路上の事象に対しては、その事象の発生を検知しても
外部に知らせないことを特徴とする請求項8記載の交通
流の異常検知装置。
10. The traffic flow according to claim 8, wherein the information providing means does not notify the outside of the event on the road which has already been scheduled even if it detects the occurrence of the event. Anomaly detection device.
【請求項11】前記交通流監視手段は、判定の基礎とな
った評価値の大きさに応じて、段階的な判定を行い、前
記情報提供手段は、この交通流監視手段による段階的な
判定の結果によって異常情報の内容を変えることを特徴
とする請求項8記載の交通流の異常検知装置。
11. The traffic flow monitoring means makes a stepwise determination according to the magnitude of the evaluation value on which the determination is based, and the information providing means makes a stepwise determination by the traffic flow monitoring means. 9. The traffic flow abnormality detection device according to claim 8, wherein the content of the abnormality information is changed according to the result of the above.
【請求項12】前記交通流監視手段は、複数の道路区間
で交通流の異常を検知した場合、判定の基礎となった評
価値の大きさに応じて、異常発生発生区間を特定するこ
とを特徴とする請求項1記載の交通流の異常検知装置。
12. When the traffic flow monitoring means detects a traffic flow abnormality in a plurality of road sections, the traffic flow monitoring means specifies the abnormality occurrence section in accordance with the magnitude of the evaluation value on which the determination is based. The traffic flow abnormality detection device according to claim 1, wherein:
【請求項13】複数の車両観測地点において、道路を走
行する車両を計測し、 この計測結果に基づいて、通過車両の特徴量を観測し、 前記特徴量に基づいて、それぞれの車両観測地点を通過
した車両同士を対応付けることにより、車両群の並びの
一致度を判定し、この判定結果に基づいて、道路上の突
発事象の発生を検知することを特徴とする交通流の異常
検知方法。
13. A vehicle traveling on a road is measured at a plurality of vehicle observation points, a characteristic amount of a passing vehicle is observed based on the measurement result, and each vehicle observation point is identified based on the characteristic amount. A traffic flow abnormality detection method characterized by determining the degree of coincidence of a group of vehicles by associating vehicles that have passed with each other, and detecting occurrence of a sudden event on a road based on the determination result.
JP2001229852A 2000-09-22 2001-07-30 Traffic flow anomaly detection apparatus and method Expired - Fee Related JP3820935B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2001229852A JP3820935B2 (en) 2000-09-22 2001-07-30 Traffic flow anomaly detection apparatus and method

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2000-289296 2000-09-22
JP2000289296 2000-09-22
JP2001229852A JP3820935B2 (en) 2000-09-22 2001-07-30 Traffic flow anomaly detection apparatus and method

Publications (2)

Publication Number Publication Date
JP2002170190A true JP2002170190A (en) 2002-06-14
JP3820935B2 JP3820935B2 (en) 2006-09-13

Family

ID=26600567

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2001229852A Expired - Fee Related JP3820935B2 (en) 2000-09-22 2001-07-30 Traffic flow anomaly detection apparatus and method

Country Status (1)

Country Link
JP (1) JP3820935B2 (en)

Also Published As

Publication number Publication date
JP3820935B2 (en) 2006-09-13

Similar Documents

Publication Publication Date Title
Lee et al. An application of a deep learning algorithm for automatic detection of unexpected accidents under bad CCTV monitoring conditions in tunnels
Ma et al. A wireless accelerometer-based automatic vehicle classification prototype system
US11380105B2 (en) Identification and classification of traffic conflicts
US9154741B2 (en) Apparatus and method for processing data of heterogeneous sensors in integrated manner to classify objects on road and detect locations of objects
JP6180482B2 (en) Methods, systems, products, and computer programs for multi-queue object detection and analysis (multi-queue object detection and analysis)
JP3753014B2 (en) Traffic flow anomaly detection apparatus and method
US20220004818A1 (en) Systems and Methods for Evaluating Perception System Quality
JP5992681B2 (en) Unusual condition detection system for congestion
Sanchez et al. Vehicle re-identification using wireless magnetic sensors: Algorithm revision, modifications and performance analysis
JP2007026300A (en) Traffic flow abnormality detector and traffic flow abnormality detection method
Nguyen et al. Real-time validation of vision-based over-height vehicle detection system
CN116453064B (en) Method for identifying abnormal road conditions of tunnel road section based on monitoring data
Yang et al. Use of sensor data to identify secondary crashes on freeways
WO2022046951A1 (en) Statistical image processing-based anomaly detection for cable cut prevention
CN114495011A (en) Non-motor vehicle and pedestrian illegal intrusion identification method based on target detection, storage medium and computer equipment
CN105405297B (en) A kind of automatic detection method for traffic accident based on monitor video
CN104103079B (en) The intrusion detection method of a kind of moving target and device
JP2020101908A (en) Traffic state prediction device and traffic state prediction method
JP2002083394A (en) Device and method for detecting abnormality in traffic flow
Kezebou et al. A deep neural network approach for detecting wrong-way driving incidents on highway roads
JP2002170190A (en) Abnormality detection device in traffic stream and method therefor
JP2024500979A (en) Road traffic extraction for unknown anomaly detection using distributed optical fiber sensing
Mostafi et al. Real-Time Jaywalking Detection and Notification System using Deep Learning and Multi-Object Tracking
Hubner et al. Audio-video sensor fusion for the detection of security critical events in public spaces
CN112200835A (en) Traffic accident detection method and device, electronic equipment and storage medium

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20050517

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20050607

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20050803

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20060530

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20060612

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090630

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100630

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110630

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110630

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120630

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130630

Year of fee payment: 7

LAPS Cancellation because of no payment of annual fees