JP2009262702A - Safe driving support system - Google Patents

Safe driving support system Download PDF

Info

Publication number
JP2009262702A
JP2009262702A JP2008113012A JP2008113012A JP2009262702A JP 2009262702 A JP2009262702 A JP 2009262702A JP 2008113012 A JP2008113012 A JP 2008113012A JP 2008113012 A JP2008113012 A JP 2008113012A JP 2009262702 A JP2009262702 A JP 2009262702A
Authority
JP
Japan
Prior art keywords
state
driver
safe driving
model
support system
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
JP2008113012A
Other languages
Japanese (ja)
Other versions
JP5161643B2 (en
Inventor
Shunsuke Tezuka
俊介 手塚
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.)
Subaru Corp
Original Assignee
Fuji Heavy 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 Fuji Heavy Industries Ltd filed Critical Fuji Heavy Industries Ltd
Priority to JP2008113012A priority Critical patent/JP5161643B2/en
Publication of JP2009262702A publication Critical patent/JP2009262702A/en
Application granted granted Critical
Publication of JP5161643B2 publication Critical patent/JP5161643B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

<P>PROBLEM TO BE SOLVED: To support safe driving by adaptively learning relevance between driving operation of a driver and an external environment corresponding to various state changes in cycle and speed generated to the driver during driving, and by recognizing deviation from an ordinary internal state of the driver. <P>SOLUTION: A correspondence between risk of travel environment recognized in an environment risk recognition part 2 and a characteristic amount of driving operation data quantized by an operation characteristic amount discretization part 5 is learned in a model learning part 6 to obtain a model parameter according to the cycle set by a schedule control part 9. A state estimation part 7 estimates an internal state of the driver from the relevance between the operation characteristic amount and an environment risk level by use of the already learned model parameter by the estimated cycle set by the schedule control part 9. An alarm-support part 8 compares the internal state of the driver and the environment risk level, evaluates whether a vehicle is in an essentially safe or dangerous state, and performs an alarm, the operation support, and the like. <P>COPYRIGHT: (C)2010,JPO&INPIT

Description

本発明は、移動体を運転するドライバの状態を推定して安全運転のための支援を行う安全運転支援システムに関する。   The present invention relates to a safe driving support system that estimates the state of a driver who drives a moving body and provides support for safe driving.

近年、自動車等の移動体の安全な移動を可能とすることを目的として、運転中のドライバの状態を監視して覚醒度の低下や居眠りを検知し、安全を確保する技術が開発されている。ドライバ状態の推定は、主として、生体状態を計測して行うものと、運転操作データから推定するものとがあり、これらの技術に関して従来から各種提案がなされている。   In recent years, with the goal of enabling safe movement of moving objects such as automobiles, techniques have been developed to ensure the safety by monitoring the state of the driver while driving to detect a decrease in wakefulness and falling asleep. . The estimation of the driver state is mainly performed by measuring the living body state and that estimated from the driving operation data, and various proposals have been made regarding these techniques.

生体状態の計測によるドライバ状態の推定に関する技術としては、特許文献1〜5に開示の技術がある。特許文献1の技術は、ドライバの顔面に取り付けた電極や眼球を撮影するカメラ等を用いて眼球運動を計測することでドライバの状態を推定するものであり、特許文献2の技術は、ドライバの瞼の開度を元に覚醒度を推定するものである。   As a technique related to estimation of a driver state by measuring a biological state, there are techniques disclosed in Patent Documents 1 to 5. The technology of Patent Document 1 estimates the driver's state by measuring eye movement using an electrode attached to the driver's face, a camera that captures the eyeball, and the like. The degree of arousal is estimated based on the opening of the heel.

また、特許文献3や特許文献4の技術は、ドライバの心拍信号を計測することでドライバの状態の推定を行うものであり、特許文献5の技術は、脳内電流信号を検知するものである。   Moreover, the technique of patent document 3 and patent document 4 estimates a driver's state by measuring a driver's heartbeat signal, and the technique of patent document 5 detects a brain current signal. .

一方、運転操作データによるドライバ状態の推定に関する技術としては、特許文献6,7に開示の技術がある。特許文献6の技術は、ドライバの運転操作データを予め決められた条件で判別し、車両制御特性の変更に反映させるものであり、特許文献7の技術は、車両状態データからのファジー推論によってドライバの意図や心理状態の推定を行うものである。
特開2003−230552号公報 特開2004−89272号公報 特許第3596198号公報 特開平7−10024号公報 特表2006−524157号公報 特許第3058966号公報 特許第3036183号公報
On the other hand, there are technologies disclosed in Patent Documents 6 and 7 as technologies related to estimation of the driver state based on driving operation data. The technique of Patent Document 6 discriminates driving operation data of a driver under a predetermined condition and reflects it in a change in vehicle control characteristics. The technique of Patent Document 7 is based on fuzzy inference from vehicle state data. Is intended to estimate the intention and psychological state.
JP 2003-230552 A JP 2004-89272 A Japanese Patent No. 3596198 Japanese Patent Laid-Open No. 7-10024 JP-T-2006-524157 Japanese Patent No. 3058966 Japanese Patent No. 3036183

前述したように、自動車の安全な運転を実現するためには、ドライバの状態を推定する技術を用いて警報や制御特性の変更・操作補助などを行うことが有効であり、これまでに各種提案がなされている。   As mentioned above, in order to realize safe driving of automobiles, it is effective to perform warnings, change control characteristics, assist operations, etc. using technology that estimates the state of the driver. Has been made.

しかしながら、文献1〜5に開示されているような、ドライバ状態推定のために生体計測を行う技術では、一般的に、ノイズや個人差の大きい測定量を扱わなくてはならず,自動車の使用される環境、使用する人間の多様性に対応しきれないという問題がある。   However, in the technology for performing living body measurement for estimating the driver state as disclosed in References 1 to 5, generally, a measurement amount having a large noise or individual difference must be dealt with. There is a problem that it cannot cope with the environment and the diversity of people who use it.

すなわち、特許文献1では、眼球運動の計測のため顔面に電極を配置する例を挙げているが、実際の自動車運転時にドライバに身体的拘束を伴う測定装置を装着しなければならないことは、煩わしいばかりでなく現実性に欠ける。   That is, Patent Document 1 gives an example in which electrodes are arranged on the face for measurement of eye movement, but it is troublesome that a driver must be equipped with a measurement device that involves physical restraint during actual driving. Not only lacks reality.

特許文献2では、ドライバの目を撮影するカメラからの情報を基にデータを計測しているが、車両に付加的な装置を付けるコストや、ドライバの姿勢の変化や西日等の外部環境からの光の影響によって撮影が正常に行えない虞がある。   In Patent Document 2, data is measured based on information from a camera that captures the eyes of the driver. However, from the cost of adding additional devices to the vehicle, changes in the driver's attitude, and the external environment such as West Japan There is a possibility that photographing cannot be performed normally due to the influence of light.

特許文献3ではシートならびにステアリングに設置した電極その他の装置、また特許文献4では心理状態検出手段とされるだけの記載であるが、共に心拍を計測する例が挙げられている。しかし、これらの技術も、心電という比較的ノイズに攪乱されやすいデータを計測していることや、車両に付加的な装置が必要であること、さらに各個人によって差異のある心拍の特徴をもって判定を行わなくてはならない点等、自動車の実環境で使用されることを鑑みるに問題があると言わざるを得ない。   Patent Document 3 describes only electrodes and other devices installed on the seat and the steering, and Patent Document 4 merely describes psychological state detection means, but both provide examples of measuring heartbeats. However, these technologies also measure the electrocardiogram, which is relatively easily disturbed by noise, require additional equipment on the vehicle, and determine the characteristics of the heartbeat that vary from individual to individual. It must be said that there is a problem in view of the fact that it is used in the actual environment of an automobile, such as the point that must be performed.

特許文献5では、MRI等の装置を車載することによる非侵襲での計測が説明されているが、文献内で述べられている装置は、直接的或いは間接的にドライバの姿勢その他の条件を拘束するものであり、そのコスト等も含めて、同様に現時点での工業的応用は現実的ではないと言わざるを得ない。   Patent Document 5 describes non-invasive measurement by mounting a device such as MRI on the vehicle, but the device described in the document directly or indirectly restrains the driver's posture and other conditions. Similarly, it must be said that the industrial application at the present time is not realistic, including its cost.

また、特許文献1〜5に共通する事項として、事前のキャリブレーションがあるにせよ、一般的な人間に共通と思われる兆候をもって覚醒度や疲労度を測ろうとしており、各個人の違いや外部環境との関連性については積極的に扱わず、むしろ個人的差異を前処理で平滑化させることに注力しているものも見受けられる。このような手法では、ある程度の効果は見込めるものの、その先のより高い精度でドライバ状態を推定する場合には、個人差や外部環境を積極的に扱う仕組みを持たないために。その性能に一定の限界がある。   In addition, as a matter common to Patent Documents 1 to 5, we are trying to measure the degree of arousal and fatigue with signs that are considered common to general humans even though there is advance calibration. There are some that do not actively deal with environmental relevance, but rather focus on smoothing out personal differences through preprocessing. Although such a method can be expected to some extent, there is no mechanism to actively handle individual differences and the external environment when estimating the driver state with higher accuracy. There is a certain limit to its performance.

一方、ドライバ状態推定のために操作あるいは車両データを用いる特許文献6,7では、ドライバの操作が予め決められた条件を満たす場合に、その意図に対応すると思われる予め決められた方向へ車両特性を変化させている。この方法では、一般的な人間に共通と考えられる操作の特徴をもって車両特性を変えているため、真の意味でドライバ一人ひとりの個性に合わせた状態推定がなされるわけでなく、人によっては誤った状態に推定されてしまう虞がある。   On the other hand, in Patent Documents 6 and 7, which use operation or vehicle data for estimating the driver state, when the driver's operation satisfies a predetermined condition, the vehicle characteristics in a predetermined direction that seems to correspond to the intention. Is changing. In this method, the vehicle characteristics are changed with the characteristics of the operation that are considered common to general humans, so the state estimation according to each driver's individuality is not made in the true sense, and it is wrong for some people There is a risk of being estimated to the state.

以上のように、従来の技術は、ノイズや付加コストについて考慮すべき課題があるばかりでなく、ドライバ状態の推定に関して画一的指標・条件による能力不足があり、更には、運転中にドライバに生じる周期や速度が様々な状態変化を広く推定することが困難であり、自動車の安全な運転を実現するためには不十分である。   As described above, the conventional technology not only has problems to be considered with respect to noise and additional costs, but also has insufficient capability due to uniform indicators and conditions for estimating the driver state. It is difficult to broadly estimate various state changes due to the period and speed generated, which is insufficient to realize safe driving of a car.

本発明は上記事情に鑑みてなされたもので、ドライバの運転操作と外部環境との関連性を、運転中にドライバに生じる周期や速度が様々な状態変化に対応して適応的に学習し、ドライバの普段の内部状態からの逸脱を認識して安全な運転を支援することのできる安全運転支援システムを提供することを目的としている。   The present invention has been made in view of the above circumstances, adaptively learning the relationship between the driving operation of the driver and the external environment in response to various state changes in the period and speed generated in the driver during driving, An object of the present invention is to provide a safe driving support system that can recognize a deviation from the normal internal state of a driver and support safe driving.

上記目的を達成するため、本発明による安全運転支援システムは、移動体を運転するドライバの状態を推定して安全運転のための支援を行う安全運転支援システムであって、上記移動体の外部環境をセンシングして特徴量を抽出し、該特徴量に基づいて外部環境に含まれる環境リスクを認識する環境リスク認識部と、上記環境リスクとドライバの運転操作との対応を確率的状態遷移モデルを用いて学習し、学習モデルを構築するモデル学習部と、上記学習モデルに基づいて、走行中のドライバの運転操作の反映としての車両挙動データを含む運転操作データからドライバのリスク認識状態をドライバの内部状態として推定する状態推定部と、上記環境リスクと上記ドライバの内部状態とを比較し、上記移動体の安全運転に係る支援情報を取得する運転支援部とを備え、更に、上記確率的状態遷移モデルへのデータ入力を複数の異なる周期に設定し、上記確率的状態遷移モデルによる推定周期を適応的に可変するスケジュール制御部を備えることを特徴とする。   In order to achieve the above object, a safe driving support system according to the present invention is a safe driving support system that estimates the state of a driver driving a mobile body and provides support for safe driving, and includes an external environment of the mobile body. To extract the feature quantity and to recognize the environmental risk included in the external environment based on the feature quantity, and the probabilistic state transition model for the correspondence between the environmental risk and the driving operation of the driver. The model learning unit that learns and constructs the learning model, and based on the learning model, the driver's risk recognition state is determined from the driving operation data including vehicle behavior data as a reflection of the driving operation of the driving driver. Compare the environmental risk with the driver's internal state, and obtain support information related to safe driving of the moving body And a schedule control unit that adaptively varies the estimation period based on the stochastic state transition model by setting data input to the stochastic state transition model in a plurality of different periods. It is characterized by.

本発明によれば、ドライバの運転操作と外部環境との関連性を、運転中にドライバに生じる周期や速度が様々な状態変化に対応して適応的に学習することができ、ドライバの普段の内部状態からの逸脱を認識して安全な運転を支援することができる。   According to the present invention, the relationship between the driving operation of the driver and the external environment can be learned adaptively in response to various state changes in the period and speed generated in the driver during driving. Recognize deviations from internal conditions and support safe driving.

以下、図面を参照して本発明の実施の形態を説明する。図1〜図16は本発明の実施の一形態に係り、図1は安全運転支援システムの基本構成図、図2は安全運転支援システムの学習フェーズ及び推定フェーズを示す説明図、図3は操作特徴量と自己組織化マップのネットワークを示す説明図、図4は確率的状態遷移を示す説明図、図5は学習後の状態遷移確率を示すグラフ、図6は学習後の出力確率を示すグラフ、図7は最尤系列による状態推定の説明図、図8は前向きアルゴリズムによる状態推定の説明図、図9はタイムドリブンによる状態推定を示す説明図、図10はイベントドリブンによる状態推定を示す説明図、図11はイベントドリブンによる遷移トリガの発生を示す説明図、図12はタイムドリブン及びイベントドリブンによる状態遷移確率を示すグラフ、図13はタイムドリブン及びイベントドリブンによる内部状態推定例を示すグラフ、図14は適応的サンプリングの説明図、図15はタイムドリブンと適応的サンプリングとを併用した内部状態推定例を示すグラフ、図16は推定試験結果を示す説明図である。   Embodiments of the present invention will be described below with reference to the drawings. 1 to 16 relate to an embodiment of the present invention, FIG. 1 is a basic configuration diagram of a safe driving support system, FIG. 2 is an explanatory diagram showing a learning phase and an estimation phase of the safe driving support system, and FIG. FIG. 4 is an explanatory diagram showing a stochastic state transition, FIG. 5 is a graph showing a state transition probability after learning, and FIG. 6 is a graph showing an output probability after learning. 7 is an explanatory diagram of state estimation by the maximum likelihood sequence, FIG. 8 is an explanatory diagram of state estimation by a forward algorithm, FIG. 9 is an explanatory diagram showing state estimation by time driven, and FIG. 10 is an explanatory diagram showing state estimation by event driven. 11 is an explanatory diagram showing the occurrence of a transition trigger by event driven, FIG. 12 is a graph showing the state transition probability by time driven and event driven, and FIG. 13 is time driven. 14 is a graph showing an example of internal state estimation by event driven, FIG. 14 is an explanatory diagram of adaptive sampling, FIG. 15 is a graph showing an example of internal state estimation using both time driven and adaptive sampling, and FIG. It is explanatory drawing shown.

本発明による安全運転支援システムは、自動車等の移動体を運転する際に、外部環境に含まれるリスクとドライバが意識しているリスクとを比較してドライバの普段の内部状態からの逸脱を認識し、不必要な介入をすることなく適切な運転支援を可能とするものであり、「システムとドライバの不協和」を解消しつつ、自動車を運転する際の安全性を向上することができる。   The safe driving support system according to the present invention recognizes a deviation from the normal internal state of the driver by comparing the risk included in the external environment with the risk recognized by the driver when driving a moving body such as an automobile. Thus, it is possible to provide appropriate driving assistance without unnecessary intervention, and it is possible to improve the safety when driving a car while eliminating the “dissonance between system and driver”.

以下、本実施の形態においては、自動車の走行支援について説明する。本実施の形態の安全運転支援システムは、ドライバ状態(ドライバの内部状態)を推定するための主たるデータとして、ドライバの運転操作の反映としての車両挙動データを含めてハンドル・アクセル・ブレーキ等の操作データを用い、また、ドライバ状態推定の指標として、画一的な条件を用いることなく、ドライバの普段の運転における「車両周囲環境と操作との関連性」を確率的規範モデルとして用いている。   Hereinafter, in the present embodiment, driving assistance for automobiles will be described. The safe driving support system according to the present embodiment includes an operation of a steering wheel, an accelerator, a brake, etc. including vehicle behavior data as a reflection of the driving operation of the driver as main data for estimating the driver state (internal state of the driver). Data is used, and the “relationship between the surrounding environment of the vehicle and the operation” in the normal driving of the driver is used as a probabilistic reference model as a driver state estimation index without using uniform conditions.

その際、本システムでは、継続的に車両を運転することによってモデルを逐次更新し、個人個人に合わせて状態推定精度を向上させていくと共に、ドライバ状態の推定周期を最適に設定し、推定性能を高めるようにしている。   At that time, the system continuously updates the model by continuously driving the vehicle, improves the state estimation accuracy according to the individual, sets the driver state estimation period optimally, and To increase.

例えば、「一瞬の気の緩み」に相当する周期の推定精度を重要視すると、それより周期の長い「ぼんやり」や「覚醒度の低下」、「疲労の蓄積」といった状態に対する推定精度が悪化する虞があり、運転中のある特定の周波数の状態変化に対しては高い推定能力を持つが、それ以外の周波数の状態変化については、相対的に推定精度が悪化する場合がある。このため、本システムでは、状態推定の周期を最適に設定することで多様な内部状態の推定を可能としており、これにより、運転中にドライバに生じる周期や速度が様々な状態変化を広く推定することができる。   For example, if importance is attached to the estimation accuracy of the period corresponding to “moment relaxation of moment”, the estimation accuracy for states such as “blurred”, “decreased arousal”, and “accumulation of fatigue” with a longer period deteriorates. There is a possibility that it has a high estimation capability for a state change at a specific frequency during operation, but the estimation accuracy may be relatively deteriorated for state changes at other frequencies. For this reason, in this system, it is possible to estimate various internal states by optimally setting the state estimation cycle, thereby widely estimating state changes with various cycles and speeds generated during driving. be able to.

具体的には、図1に示すように、本実施の形態における安全運転支援システム1は、単一のコンピュータシステム或いはネットワーク等を介して接続された複数のコンピュータシステムで構成され、カメラやレーダー等のセンシングデバイスによる外界情報から車両周囲の外部環境に含まれるリスク(環境リスク)を認識する環境リスク認識部2、環境リスクを離散化する環境リスク離散化部3、ドライバの運転操作を計測する操作計測部4、計測されたドライバの運転操作データの特徴量を離散化(量子化)する操作特徴量離散化部5、運転操作(操作特徴量)と認識環境(環境リスクレベル)との対応を確率的規範モデルを用いて学習するモデル学習部6、学習されたモデルと現在の操作データ・車外環境との関連性からドライバ状態を推定する状態推定部7、推定されたドライバ状態に応じて車両を安全に運転するための支援情報を取得し、警報や操作支援などを行う警報・支援部8を備えて構成され、更に、ドライバの覚醒度や疲労度等の状態や走行状況の変化に対応して最適な推定を行うため、確率的規範モデルへのデータ入力を複数の異なる周期に設定し、確率的規範モデルによる推定周期を適応的に可変するスケジュール制御部9を備えている。   Specifically, as shown in FIG. 1, the safe driving support system 1 in the present embodiment includes a single computer system or a plurality of computer systems connected via a network or the like, such as a camera or a radar. Environmental risk recognition unit 2 for recognizing risks (environmental risks) included in the external environment around the vehicle from the external information by the sensing device, environmental risk discretization unit 3 for discretizing the environmental risk, and operation for measuring the driving operation of the driver The measurement unit 4, the operation feature quantity discretization unit 5 that discretizes (quantizes) the feature quantity of the measured driving operation data of the driver, the correspondence between the driving operation (operation feature quantity) and the recognition environment (environmental risk level). Model learning unit 6 that learns using a probabilistic normative model, the driver state from the relationship between the learned model and the current operation data / exterior environment A state estimation unit 7 to be determined, and a warning / support unit 8 that obtains support information for safely driving the vehicle according to the estimated driver state, and performs warnings and operation support, and further includes a driver In order to perform optimal estimation in response to changes in conditions such as arousal level and fatigue level and running conditions, data input to the probabilistic normative model is set to multiple different periods, and the estimation period based on the probabilistic normative model is set. An adaptively variable schedule control unit 9 is provided.

尚、本実施の形態においては、外部環境のセンシングデバイスとして車載カメラを用いている。   In this embodiment, an in-vehicle camera is used as a sensing device for the external environment.

この安全運転支援システム1のシステム動作は、便宜上、図2に示すように、ドライバの内部状態のモデルを学習する「学習フェーズ」と、得られたモデルを元にドライバの内部状態を推定する「推定フェーズ」とに分けられる。   As shown in FIG. 2, the system operation of the safe driving support system 1 includes a “learning phase” in which a model of the internal state of the driver is learned, and an internal state of the driver is estimated based on the obtained model. It is divided into “estimation phase”.

学習フェーズでは、車両前方を車載カメラで撮像し、学習によって走行環境のリスクレベルを認識すると、このリスクレベルと、ドライバの運転操作データの特徴量を量子化した操作特徴量データとの対応関係を学習し、学習済みパラメータを取得して学習モデルを構築する。   In the learning phase, when the front of the vehicle is imaged with an in-vehicle camera and the risk level of the driving environment is recognized by learning, the correspondence between this risk level and the operation feature value data obtained by quantizing the feature value of the driver's drive operation data is obtained. Learn and acquire learned parameters to build a learning model.

また、推定フェーズでは、学習済みのモデルパラメータを用いて操作特徴量と環境リスクレベルとの関連性を獲得し、現在のドライバの内部状態を推定する。推定されたドライバの現在の内部状態は、現在の走行環境における環境リスクレベルと比較され、車両が本質的に安全か危険な状態にあるかが評価される。   In the estimation phase, the relationship between the operation feature amount and the environmental risk level is acquired using the learned model parameter, and the current internal state of the driver is estimated. The estimated current internal state of the driver is compared with the environmental risk level in the current driving environment to assess whether the vehicle is inherently safe or dangerous.

この学習フェーズと推定フェーズとは、実際には互いに分離して進行するものではなく、双方が同時に進行する。それにより、ドライバが走行する環境の変化やドライバ自身の特性変化に適応する機能を実現することができ、車両のユーザに対して、走行すればするほど高性能になっていくシステムを提供することができる。また、学習フェーズにおけるモデルパラメータの取得は、走行中のオンライン学習によって全てを取得することも可能であるが、一部、シミュレータ等を用いたオフライン学習を併用することにより、学習速度の向上と認識精度の向上とを期待することができる。   The learning phase and the estimation phase do not actually proceed separately from each other, but both proceed simultaneously. As a result, it is possible to realize a function that can adapt to changes in the environment in which the driver is driving and changes in the characteristics of the driver itself, and to provide a system that increases the performance as the vehicle travels. Can do. In addition, the model parameters in the learning phase can all be acquired by online learning while driving. However, it is recognized that the learning speed can be improved by using offline learning in part with a simulator. An improvement in accuracy can be expected.

次に、安全運転支援システム1の各部の機能について説明する。環境リスク認識部2は、車載カメラから得られる画像の特徴量と、そのときのリスクレベルとの関連性を学習することにより、車外環境の状態を単一のスカラ値(若しくはベクトル)へ縮退変換する。環境リスク離散化部3は、環境リスク認識部2からの環境リスクレベルを閾値を用いたデータ分割等により離散化し、確率的規範モデルの観測ステートとしてモデル学習部6へ伝達する。   Next, the function of each part of the safe driving support system 1 will be described. The environmental risk recognition unit 2 degenerates the state of the environment outside the vehicle into a single scalar value (or vector) by learning the relationship between the feature quantity of the image obtained from the in-vehicle camera and the risk level at that time. To do. The environmental risk discretization unit 3 discretizes the environmental risk level from the environmental risk recognition unit 2 by data division or the like using a threshold value, and transmits it to the model learning unit 6 as an observation state of the probabilistic normative model.

画像特徴量からのリスクレベルの認識は、例えば、本出願人による特願2007−77625号において提案されたオンラインリスク学習システムの技術を採用することができる。この技術は、特願2007−77625号に詳述されているように、アクセルの急激な戻し操作やブレーキ踏込みといったイベントにより、N次元ベクトルの画像特徴量(エッジ情報、動き情報、色情報等)を1次元の状態に変換し、この状態と車両情報(ドライバの操作情報)から作成された教師情報との相関関係から環境に含まれるリスクを学習・認識するものである。   For the recognition of the risk level from the image feature amount, for example, the technology of the online risk learning system proposed in Japanese Patent Application No. 2007-77625 by the present applicant can be adopted. As described in detail in Japanese Patent Application No. 2007-77625, this technique is based on an image feature amount (edge information, motion information, color information, etc.) of an N-dimensional vector depending on an event such as a sudden accelerator return operation or a brake depression. Is converted into a one-dimensional state, and the risk contained in the environment is learned and recognized from the correlation between this state and the teacher information created from the vehicle information (driver operation information).

尚、本実施の形態では、画像特徴量から抽出したリスクレベルを用いる例について説明するが、リスクレベルとしては、これに限定されるものではなく、例えば、車間距離等からリスクレベルを抽出するようにしても良い。   In the present embodiment, an example using the risk level extracted from the image feature amount will be described. However, the risk level is not limited to this, and for example, the risk level is extracted from the inter-vehicle distance or the like. Anyway.

操作特徴量離散化部5は、センサからのデータ或いは車内ネットワーク(図示せず)を介して操作計測部4で取得したアクセル・ブレーキ・ステアリング等のドライバの操作データを、その測定頻度の特徴に応じて量子化し、確率的規範モデルの観測シンボルとしてモデル学習部6及び状態推定部7へ伝達する。   The operation feature quantity discretization unit 5 uses the data from the sensor or the operation data of the driver such as accelerator, brake, steering, etc. acquired by the operation measurement unit 4 via the in-vehicle network (not shown) as the characteristics of the measurement frequency. In response, the signal is quantized and transmitted to the model learning unit 6 and the state estimation unit 7 as an observation symbol of the probabilistic reference model.

すなわち、運転操作データは、そのままでは情報量が非常に多く、リスクとの関連を学習するには扱いが困難である。このため、観測されるデータの分布(出現傾向)を考慮して適切に量子化を行うことで、データに含まれる情報量が失われることを防止しつつ、特徴を学習するための統計的処理を可能とする。   In other words, the driving operation data has a very large amount of information as it is, and it is difficult to handle the relationship with the risk. For this reason, statistical processing to learn features while preventing loss of information contained in the data by appropriately quantizing in consideration of the distribution (appearance tendency) of the observed data Is possible.

観測データの量子化は、閾値を用いたデータ分割やデータ縮退によって行うことができ、本実施の形態では、自己組織化マップ(SOM;Self-Organizing Maps)を用いて運転操作データを量子化する。SOMは、生物の大脳皮質のうち視覚野等をモデル化したニューラルネットワークの一種であり、M次元に並べられたユニットが、それぞれベクトル値(通常入力との結線の重みと呼ばれる)を持ち、入力に対して勝者ユニットをベクトルの距離を基準として決定する。   Quantization of observation data can be performed by data division or data degeneration using thresholds. In this embodiment, driving operation data is quantized using a self-organizing map (SOM). . SOM is a type of neural network that models the visual cortex of the cerebral cortex of living organisms. Each unit arranged in M dimensions has a vector value (usually called the weight of the connection with the input). The winner unit is determined based on the vector distance.

そして、勝者ユニット及びその周辺のユニットの参照ベクトル値を、入力ベクトルに近づくように更新してゆく。これを繰り返すことで、全体が入力データの分布を最適に表現できるように競合学習し、この競合学習に基づいて入力情報の次元を圧縮すると共に、データの特徴に応じてクラスタリングや可視化を行うことができる。   Then, the reference vector values of the winner unit and the surrounding units are updated so as to approach the input vector. By repeating this, competitive learning is performed so that the entire distribution of the input data can be expressed optimally, the dimensions of the input information are compressed based on this competitive learning, and clustering and visualization are performed according to the characteristics of the data Can do.

尚、本実施の形態では、S0Mを用いた教師無し競合学習により入力データを縮退する例について説明するが、教師有り競合学習であるベクトル量子化(LVQ;Learning Vector Quantization)モデルを用いることも可能である。   In this embodiment, an example in which input data is degenerated by unsupervised competitive learning using S0M will be described. However, a vector quantization (LVQ; Learning Vector Quantization) model that is supervised competitive learning can also be used. It is.

以下に、SOMパラメータの例を示す。採用した特徴量は、舵角速度、アクセル開度・ブレーキ圧・アクセル開度変化量・ブレーキ圧変化量・車両速度である。但し、舵角速度については絶対値を用い、また、アクセル開度とブレーキ圧は、それぞれ正規化し、-100(ブレーキ圧最大)から+100(アクセル全開)までの無次元量としている。アクセル並びにブレーキの変化量も、この正規化値の変化量として表している。結果として入力特徴量の次元は4とし、これをSOMによって1次元に圧縮し、ユニット数256に量子化している。   Examples of SOM parameters are shown below. The feature values adopted are the steering angular speed, accelerator opening, brake pressure, accelerator opening change, brake pressure change, and vehicle speed. However, absolute values are used for the rudder angular velocity, and the accelerator opening and the brake pressure are normalized to be dimensionless amounts from -100 (maximum brake pressure) to +100 (accelerator fully open). The change amount of the accelerator and the brake is also expressed as the change amount of the normalized value. As a result, the dimension of the input feature quantity is 4, which is compressed to one dimension by SOM and quantized to 256 units.

[SOMパラメータ]
SOM次元数 1
ユニット数 256
入力特徴量 ADSTR(舵角速度絶対値)
PEDAL(アクセル・ブレーキ操作正規化量)
DPEDAL(PEDAL値変化量)
SPEED(車両速度)
[SOM parameters]
Number of SOM dimensions 1
256 units
Input feature value ADSTR (steering angular velocity absolute value)
PEDAL (accelerator / brake operation normalized amount)
DPEDAL (PEDAL value change amount)
SPEED (vehicle speed)

図3は、SOMで学習されたデータを示しており(但し、入力4次元のうち3次元のみ)、市街地を走行したデータ(5時間走行相当)をSOMで学習したものである。雲のようにプロットされているのが測定データで、ひも状に配置されているのが1次元SOMのネットワークである。図3から測定データの密度と分布に応じて学習が行なわれ、適宜、ユニットが配置されていることが分かる。測定されたデータは、最近傍のSOMのユニットで近似されることで、量子化されると共に、ユニット番号で表される1次元のデータに縮退される。   FIG. 3 shows data learned by SOM (however, only three of the four input dimensions), and data obtained by traveling in an urban area (corresponding to 5-hour driving) is learned by SOM. The measurement data is plotted like a cloud, and the one-dimensional SOM network is arranged in a string shape. It can be seen from FIG. 3 that learning is performed according to the density and distribution of measurement data, and that units are appropriately arranged. The measured data is quantized and reduced to one-dimensional data represented by a unit number by being approximated by the nearest SOM unit.

尚、このSOMのユニット配置は、計測される運転データを学習することで決まり、個人による差や、同一人物でも時間の推移による運転特徴の変化に適応していくことが期待できる。   The unit arrangement of the SOM is determined by learning measured driving data, and it can be expected that the same person will adapt to changes in driving characteristics due to time variations even for the same person.

また、SOM及び環境リスクレベルの離散化の周期は、予め固定した単一の周期ではなく、スケジュール制御部9により状態推定の周期に合わせて設定される。この離散化周期と状態推定の周期については、後述する。   Further, the period for discretizing the SOM and the environmental risk level is not a single fixed period, but is set by the schedule control unit 9 according to the state estimation period. The discretization period and the state estimation period will be described later.

モデル学習部6は、操作特徴量とリスクレベルとの関連性を学習によって獲得し、ドライバの内部状態を推定するための確率的計算を可能とする。人間の行動は、図4(a)に示すように、安心、緊張、不安、焦り、怒りといった心的状態と、その遷移に応じて変化し、必ずしも確定的ではなく、確率的な行動として表現することができる。同様に、ドライバの運転行動は、図4(b)に示すように、先行車への追従、追越、駐車、車線変更、合流といったシーンと、その遷移に対して、確率的な操作出力となって現れる。   The model learning unit 6 acquires the relationship between the operation feature amount and the risk level by learning, and enables probabilistic calculation for estimating the internal state of the driver. As shown in FIG. 4 (a), human behavior changes according to mental states such as relief, tension, anxiety, impatience, and anger, and transitions thereof, and is not necessarily deterministic but expressed as probabilistic behavior. can do. Similarly, as shown in FIG. 4 (b), the driver's driving behavior includes a stochastic operation output for scenes such as following a preceding vehicle, passing, parking, lane change, and merging. It appears.

従って、モデル学習部6では、人間行動の確率的振る舞いを表すための規範モデルとして、確率的状態遷移モデルの一種である隠れマルコフモデル(HMM;Hidden Markov Model)を用いてドライバの内部状態をモデル化する。HMMは、対象の内部状態(ステート)が確率的な条件分岐によって遷移することと、遷移したステートによって異なる確率で外部に信号が出力されることを想定したモデルである。   Therefore, the model learning unit 6 models the internal state of the driver using a hidden Markov model (HMM), which is a kind of probabilistic state transition model, as a reference model for expressing the probabilistic behavior of human behavior. Turn into. The HMM is a model that assumes that a target internal state (state) transitions by a probabilistic conditional branch and that a signal is output to the outside with a different probability depending on the transitioned state.

HMMを用いたモデルでは、ドライバが意識しているリスクレベルを推定するタスクにおいては、図4(c)に示すように、現在意識しているリスクレベルがHMMのステートに相当し、そのときに観測される運転操作データが外部に出力される信号に相当する。図4(c)においては、リスクレベルを5段階として各ステートに1〜5の番号を付与した場合を例示しており、番号1は、ドライバが外部環境のリスクが最も低いと認識している状態、番号5は、ドライバが外部環境のリスクが最も高いと認識している状態を示している。尚、実際に扱うリスクレベルは、0,1,…,9の10段階である。   In the model using the HMM, in the task of estimating the risk level that the driver is aware of, the risk level that is currently conscious corresponds to the state of the HMM, as shown in FIG. Observed driving operation data corresponds to a signal output to the outside. FIG. 4C illustrates an example in which the risk levels are set to 5 levels and numbers 1 to 5 are assigned to the respective states. The number 1 recognizes that the risk of the external environment is the lowest. The state, number 5, indicates a state where the driver recognizes that the risk of the external environment is the highest. In addition, the risk level that is actually handled is 10 levels of 0, 1,.

このように、運転操作データを離散化してデータの出現傾向を求め、ドライバの内部状態を確率的モデルとして近似することで、実際の運転における操作データのように、外部環境から確定的に導出することのできない情報の扱いを適切に行うことが可能となる。但し、推定時に観測される操作データがどのステートから出力されたのかを推定するためには、状態遷移確率と操作出力確率の2つの確率計算を行う必要がある。   In this way, the driving operation data is discretized, the appearance tendency of the data is obtained, and the internal state of the driver is approximated as a probabilistic model, so that it is deterministically derived from the external environment like the operation data in actual driving. It is possible to appropriately handle information that cannot be handled. However, in order to estimate from which state the operation data observed at the time of estimation is output, it is necessary to perform two probability calculations: a state transition probability and an operation output probability.

このため、モデル学習部6では、先ず環境リスク認識部2から環境リスク離散化部3を介して伝達されたスカラ値を基に、統計的手法によって状態遷移確率を学習し、次に、操作特徴量離散化部5からのデータを基に、ステート毎の操作特徴量の観測確率分布(操作出力確率)を学習する。例えば、舵角・アクセル開度・ブレーキ圧力、さらに、操作の反映である速度・ヨーレート・加速度等をそれぞれ適宜離散化し、各離散値毎の観測回数をカウントして統計的に確率を計算する。   For this reason, the model learning unit 6 first learns the state transition probability by a statistical method based on the scalar value transmitted from the environmental risk recognition unit 2 via the environmental risk discretization unit 3, and then operates the operation feature. Based on the data from the quantity discretization unit 5, the observation probability distribution (operation output probability) of the operation feature amount for each state is learned. For example, the steering angle, the accelerator opening, the brake pressure, and the speed, yaw rate, acceleration, and the like that reflect the operation are discretized as appropriate, and the number of observations for each discrete value is counted to statistically calculate the probability.

以下、HMMを用いた状態遷移確率及び操作出力確率の計算処理について説明する。   Hereinafter, calculation processing of the state transition probability and the operation output probability using the HMM will be described.

[HMMの仕様例]
駆動周波数 可変
ステート数 10
出力信号 SOMユニット番号(0〜255)
ネットワーク構造 全ステート相互接続型
[Specification example of HMM]
Drive frequency variable Number of states 10
Output signal SOM unit number (0 to 255)
Network structure All-state interconnection type

[状態遷移確率の計算]
一般に、自動車の運転を想定する場合、車を運転する度に故意ではないにも拘わらず必ず事故を起こしてしまうドライバが存在するとは考えにくい。言い換えれば、巨視的にはドライバは適宜走行環境のリスクレベルに対応した運転操作を行っていると見なすことができる。このような前提の元に、ある程度長い時間範囲のデータを集めて統計的に処理した場合、HMMにおけるドライバの内部状態の遷移は、走行環境のリスクレベルの遷移に依存していると仮定することができる。
[Calculation of state transition probability]
In general, when driving a car, it is unlikely that there will be a driver who always causes an accident even though it is unintentional every time the car is driven. In other words, it can be regarded macroscopically that the driver is appropriately performing a driving operation corresponding to the risk level of the driving environment. Based on these assumptions, when data is collected and statistically processed over a relatively long time range, it is assumed that the transition of the internal state of the driver in the HMM depends on the transition of the risk level of the driving environment. Can do.

従って、遷移確率の計算は、以下の(1−1)〜(1−5)の手順に従って実施し、学習時のリスクレベルの遷移確率を求めることで、ドライバ内部状態の遷移確率を算出する。
(1−1)状態遷移モデルの駆動周波数を定義
(1−2)内部状態のステート数を定義(例えば、10ステート)
(1−3)リスクレベルをステート数分に離散化
(1−4)各ステート間の遷移回数をカウント
(1−5)各ステート間の統計的な遷移確率を計算
Therefore, the transition probability is calculated according to the following procedures (1-1) to (1-5), and the transition probability of the driver internal state is calculated by obtaining the transition probability of the risk level at the time of learning.
(1-1) Define drive frequency of state transition model (1-2) Define number of states in internal state (for example, 10 states)
(1-3) Discrete the risk level by the number of states (1-4) Count the number of transitions between each state (1-5) Calculate the statistical transition probability between each state

図5に、市街地での計測データから学習で獲得したHMMの状態遷移確率のグラフを示す。ステート番号は大きいほど高いリスクレベルであることを示している。このグラフは、所定の駆動周波数でのHMMによる学習結果を示すのであり、学習結果から、自己遷移(同じ番号のステートにとどまる遷移)と、上下1ステート分の遷移確率が大きいこと、高リスクステートでは急激にリスクレベルが下がる遷移確率も比較的大きいこと等が見て取れる。これは、一般的に、中低程度のリスクは連続的に上下することや、高いリスク要因は車両の通過と共に急激に解消すること等に相当し、一般的な運転状況の感覚との乖離は無いものと判断することができる。   FIG. 5 shows a graph of the state transition probability of the HMM obtained by learning from the measurement data in the urban area. The larger the state number, the higher the risk level. This graph shows the learning result by the HMM at a predetermined driving frequency. From the learning result, the self-transition (transition staying in the same numbered state) and the transition probability for one state above and below are large, the high risk state Then, it can be seen that the transition probability that the risk level suddenly decreases is relatively large. This is generally equivalent to the fact that moderate and moderate risks rise and fall continuously, and high risk factors disappear rapidly with the passing of the vehicle. It can be judged that there is nothing.

[出力信号確率の計算]
あるステートに遷移した状態における、操作データ特徴の出力確率の学習には、前述したSOMを用いる。計測された操作データをSOMユニット番号に次元圧縮並びに量子化し、各ユニット番号が観測される回数をカウントすることで、統計的な出力信号確率を算出する。ステート毎に得られる各ユニットの出力確率が異なることで、リスクレベルによって操作傾向が変わることをモデル化することができ、図6に示すようなステート別の出力確率を得ることができる。
[Calculation of output signal probability]
The SOM described above is used for learning the output probability of the operation data feature in the state transitioned to a certain state. The measured operation data is dimensionally compressed and quantized into SOM unit numbers, and the statistical output signal probability is calculated by counting the number of times each unit number is observed. Since the output probability of each unit obtained for each state is different, it can be modeled that the operation tendency changes depending on the risk level, and the output probability for each state as shown in FIG. 6 can be obtained.

以上の状態遷移確率と出力確率が学習によって獲得されると、次に、状態推定部7では、両者の情報を用いてドライバの内部状態を推定する。この内部状態は、最尤系列推定による手法や前向きアルゴリズム(forward algorithm)を用いた手法で推定することができ、システムの条件等に応じて適宜選択される。最尤系列推定による手法では、ドライバの内部状態が離散的な数値で表現されるステート番号で出力され、前向きアルゴリズムでは、ドライバの内部状態がステート番号を連続的な数値で緻密に表現した期待値で出力される。   When the above state transition probability and output probability are acquired by learning, the state estimation unit 7 next estimates the internal state of the driver using both pieces of information. This internal state can be estimated by a technique using maximum likelihood sequence estimation or a technique using a forward algorithm, and is appropriately selected according to system conditions and the like. In the maximum likelihood sequence estimation method, the driver's internal state is output as a discrete number of state numbers, and in the forward-looking algorithm, the driver's internal state is an expected value that expresses the state number precisely as a continuous number. Is output.

[最尤系列推定による内部状態の推定]
最尤系列推定は、現在観測されている操作データがどのステートから出力されているとするのが最も尤もらしいかを推定(最尤推定)する手法であり、遷移系列を計算する際の最も尤もらしい系列とは、発生する確率が最大となる系列を特定することに相当する。ここでは、HMMに対する最尤系列推定手法の一つであるビタビ・アルゴリズム(Viterbi algorithm)を用いて、計測された操作データの時系列データからドライバの内部状態の遷移系列を計算し、得られたステート遷移系列に沿ってドライバの内部状態が遷移していると推定する。
[Internal state estimation by maximum likelihood sequence estimation]
Maximum likelihood sequence estimation is a method that estimates the most likely state from which the currently observed operation data is output (maximum likelihood estimation), and the maximum likelihood when calculating the transition sequence. A probable sequence corresponds to specifying a sequence with the highest probability of occurrence. Here, the transition sequence of the internal state of the driver is calculated from the time series data of the measured operation data using the Viterbi algorithm, which is one of the maximum likelihood sequence estimation methods for the HMM, and obtained. It is estimated that the internal state of the driver is transitioning along the state transition series.

ビタビ・アルゴリズムは、状態遷移確率と出力信号確率を元に、時系列観測データの先頭から、順次、各時間ステップにおいてそれぞれのステートに存在する最大確率を計算する手法である。このアルゴリズムは、動的計画法に類似した計算手法を用いるために計算量が少ないという特徴を持っている。   The Viterbi algorithm is a method of calculating the maximum probability existing in each state at each time step sequentially from the beginning of the time series observation data based on the state transition probability and the output signal probability. This algorithm has a feature that the calculation amount is small because a calculation method similar to dynamic programming is used.

ビタビ・アルゴリズムでは、観測データの最後まで計算を行った結果、最終ステップでの最大確率を持つステートを、そのステップでの推定内部状態として確定する。次に、そこから1ステップずつ過去にさかのぼるバックトラックという操作を行い、最大の状態遷移確率を持つステートを確定していく。最終的には,すべての時間ステップにおいて存在する確率が最も高いステートを特定することで、最尤系列を推定する。尚、確率計算の始めのステップでは、HMMの状態遷移確率から求められる事前確率を、各ステートへの存在確率として用いる。   In the Viterbi algorithm, the state having the maximum probability at the final step is determined as the estimated internal state at that step as a result of calculation up to the end of the observation data. Next, an operation called backtracking goes back one step at a time from there, and the state having the maximum state transition probability is determined. Finally, the maximum likelihood sequence is estimated by specifying the state having the highest probability of being present in all time steps. In the first step of probability calculation, the prior probability obtained from the state transition probability of the HMM is used as the existence probability in each state.

詳細には、以下の(2−1)〜(2−4)のステップに従って逐次計算を行うことで、最尤系列を推定する。但し、π:事前確率、δ:状態存在確率、φ:バックトラック、a:状態遷移確率、b:出力信号確率、P:推定確率、q:推定状態系列、o:零ベクトルとし、各変数の添字1は、初期値であることを示す。   Specifically, the maximum likelihood sequence is estimated by performing sequential calculation according to the following steps (2-1) to (2-4). Where π: prior probability, δ: state existence probability, φ: backtrack, a: state transition probability, b: output signal probability, P: estimated probability, q: estimated state series, o: zero vector, Subscript 1 indicates an initial value.

(2−1)各状態i=1,…,Nに対して、変数の初期化を行う。
δ1(i)=π11(o1)
φ1(i)=0
(2−2)各時刻t=1,…,T−1、各状態j=1,…,Nについて、以下の再帰計算を実行する。
δt+1(j)=maxi[δt(i)aij]bj(ot+1)
φt+1(j)=argmaxi[δt(i)aij
(2−3)再帰計算の終了
P=maxiδT(i)
T=argmaxiδT(i)
(2−4)バックトラックによる最適状態遷移系列の復元(T=T−1,…,1に対して以下を実行)
T=φt+1(qt+1)
(2-1) Variables are initialized for each state i = 1,.
δ 1 (i) = π 1 b 1 (o 1 )
φ 1 (i) = 0
(2-2) The following recursive calculation is executed for each time t = 1,..., T-1 and each state j = 1,.
δ t + 1 (j) = max it (i) a ij ] b j (o t + 1 )
φ t + 1 (j) = argmax it (i) a ij ]
(2-3) End of recursive calculation P = max i δ T (i)
q T = argmax i δ T (i)
(2-4) Restoration of optimal state transition sequence by backtracking (the following is executed for T = T-1,..., 1)
q T = φ t + 1 (q t + 1 )

以上の計算では、ビタビ・アルゴリズムはバックトラックを用いる特性上、時間的に後の事象を確定してから過去方向にステートを確定させることになり、そのままではオンラインでのリアルタイム状態推定に適用することには難がある。そこで、オンライン状態推定を行う場合には、全ての時間ステップをビタビ・アルゴリズムにおける最終ステップとして扱い、各ステップにおける最大存在確率を持つステートを推定ステートとして出力する。   In the above calculation, the Viterbi algorithm has a characteristic of using backtracking, which means that after the event is confirmed in time, the state is determined in the past direction, and it can be applied to online real-time state estimation as it is. There are difficulties. Therefore, when performing online state estimation, all time steps are treated as final steps in the Viterbi algorithm, and a state having the maximum existence probability in each step is output as an estimated state.

例えば、図7(a)に示すように、オンライン状態推定では、時刻t=0における状態推定を、これまでの学習結果から決定される事前確率を用いて行い、時刻t=0で事前確率×出力確率が最大となるステートを出力し、以降の時間ステップでは、遷移確率×出力確率が最大となるステートを遷移系列として出力する。すなわち、各時間ステップにおいては、各ステートは、図7(a)に太線で示すような遷移の確率が高いが、これらの確率のうち、最大確率を持つステートを、オンライン状態推定での遷移ステートとする。その結果、図7(a)に示すオンライン状態推定では、最大存在確率を持つステート1,1,3,2,4が状態遷移系列として出力される。   For example, as shown in FIG. 7A, in the online state estimation, the state estimation at time t = 0 is performed using the prior probability determined from the learning result so far, and the prior probability at time t = 0 × The state having the maximum output probability is output, and in the subsequent time steps, the state having the maximum transition probability × output probability is output as a transition sequence. That is, in each time step, each state has a high transition probability as shown by a bold line in FIG. 7A. Of these probabilities, the state having the maximum probability is the transition state in the online state estimation. And As a result, in the online state estimation shown in FIG. 7A, the states 1, 1, 3, 2, and 4 having the maximum existence probability are output as the state transition series.

このオンライン状態推定に対して、オフライン状態推定では、時刻t4からトラックバックすると、時刻t4で最大存在確率を持つステート4は、図7(a)に太線で示すように、1ステップ前の時刻t3のステップ3からの遷移であり、時刻t3のステップ3は、時刻t2のステート2からの遷移、時刻t2のステート2は時刻t1のステート1からの遷移、時刻t2のステート2は時刻t1のステート1からの遷移であることが分かる。従って、オフライン状態推定では、図7(b)に示すように、オンライン状態推定で求めた系列のうち、時刻t2,t3のステートがトラックバックにより変更され、1→1→2→3→4の最尤系列が推定状態として出力される。   In contrast to the online state estimation, in the offline state estimation, when the track is tracked back from time t4, the state 4 having the maximum existence probability at time t4 is shown at the time t3 one step before, as indicated by the bold line in FIG. Step 3 is a transition from Step 3, Step 3 at time t3 is a transition from State 2 at time t2, State 2 at Time t2 is a transition from State 1 at Time t1, and State 2 at Time t2 is State 1 at Time t1 It turns out that it is a transition from. Therefore, in the offline state estimation, as shown in FIG. 7B, among the sequences obtained by the online state estimation, the states at times t2 and t3 are changed by the trackback, and the state of 1 → 1 → 2 → 3 → 4 is changed. A likelihood sequence is output as an estimated state.

[前向きアルゴリズムによる内部状態の推定]
以上のビタビ・アルゴリズムに対して、オンラインでの適用や推定精度の向上をより強く意図する場合には、前向きアルゴリズムを用いて内部状態を推定する。前向きアルゴリズムは、イベントと各ステートとのトレリス上を前向きに辿る手法であり、以下の(3−1),(3−2)のステップに従って逐次計算を行うことで、前向き確率(状態存在確率)αを計算する。尚、以下の式中において、π:事前確率、a:状態遷移確率、b:出力信号確率、o:零ベクトルとし、各変数の添字1は初期値であることを示す。
[Internal state estimation by forward algorithm]
When the above Viterbi algorithm is strongly intended to be applied online or to improve estimation accuracy, the forward state algorithm is used to estimate the internal state. The forward algorithm is a method of tracing forward on a trellis between an event and each state, and performing forward calculation according to the following steps (3-1) and (3-2), thereby providing a forward probability (state existence probability). α is calculated. In the following equations, π: prior probability, a: state transition probability, b: output signal probability, o: zero vector, and subscript 1 of each variable indicates an initial value.

(3−1)各状態i=1,…,Nに対して、前向き確率の初期化を行う。
α1(i)=π11(o1)
(3−2)各時刻t=1,…,T−1、各状態j=1,…,Nについて、前向き確率を再帰的に計算する。尚、再帰計算におけるΣはj=1〜Nについての総和である。
αt+1(j)=[Σαt(i)aij]bj(ot+1)
(3-1) The forward probability is initialized for each state i = 1,.
α 1 (i) = π 1 b 1 (o 1 )
(3-2) A forward probability is recursively calculated for each time t = 1,..., T−1 and each state j = 1,. Note that Σ in the recursive calculation is the sum for j = 1 to N.
α t + 1 (j) = [Σα t (i) a ij ] b j (o t + 1 )

この前向きアルゴリズムによる状態推定は、図8に示すように、時刻t=0,1,2,3,4,…の各ステップにおいて全ての遷移確率を計算しており、存在確率の低いステートも考慮した状態推定であることから、高精細な内部状態の推定が可能となる。尚、確率計算の始めのステップでは、HMMの状態遷移確率から求められる事前確率を、各ステートへの存在確率として用いる。   In the state estimation by this forward algorithm, as shown in FIG. 8, all transition probabilities are calculated at each step of time t = 0, 1, 2, 3, 4,. Therefore, the internal state can be estimated with high definition. In the first step of probability calculation, the prior probability obtained from the state transition probability of the HMM is used as the existence probability in each state.

逐次計算によって求められた各ステートへの状態存在確率αからは、以下に示すように期待値μが計算される。この期待値μがステート番号を離散的な値でなく連続的な数値で緻密に表現したドライバ内部状態として出力され、ドライバの意識しているリスクレベルとなる。尚、xは確率変数の値(ここではステート番号)であり、Σはi=1〜Nの総和である。
μ=Σxiαi
From the state existence probability α for each state obtained by the sequential calculation, an expected value μ is calculated as shown below. The expected value μ is output as a driver internal state in which the state number is precisely expressed by a continuous numerical value instead of a discrete value, and becomes a risk level that the driver is aware of. Note that x is the value of the random variable (here, the state number), and Σ is the sum of i = 1 to N.
μ = Σx i α i

[推定周期]
以上のHMMの遷移周期は、一つの周期に固定されることなく、スケジュール制御部9において、複数の異なる時間周期や観測シンボルの変化に応じた周期として管理され、多様な内部状態の推定を可能としている。具体的には、SOM及びHMMを階層化して各層を異なる時間周期で駆動する手法(タイムドリブン)、モデル前段(SOM出力)の観測シンボルをイベントとして、このイベントの状態によってHMMを駆動する手法(イベントドリブン)、適応的なサンプリング周期に合わせてHMMを駆動する手法があり、システム条件等によって適宜選択・組わせて用いられる。次に、各手法について説明する。
[Estimation period]
The transition period of the above HMM is not fixed to one period, but is managed as a period according to a plurality of different time periods and changes in observation symbols in the schedule control unit 9, and various internal states can be estimated. It is said. Specifically, the SOM and HMM are hierarchized to drive each layer at different time periods (time-driven), and the observation symbol in the previous stage of the model (SOM output) is used as an event to drive the HMM according to the state of this event ( (Event driven), there is a method of driving the HMM in accordance with an adaptive sampling period, which is appropriately selected and combined depending on system conditions and the like. Next, each method will be described.

<タイムドリブン>
タイムドリブンの特徴は、一定の周期で状態が遷移すると仮定することである。このため、タイムドリブンでは、図9に示すように、環境リスク離散化部3,操作特徴量離散化部5,モデル学習部6、状態推定部7をそれぞれ階層化し、各層を、それぞれの時間周期で並列に動作させる。
<Time driven>
The feature of time-driven is that it is assumed that the state transitions at a constant period. For this reason, in the time driven, as shown in FIG. 9, the environmental risk discretization unit 3, the operation feature quantity discretization unit 5, the model learning unit 6, and the state estimation unit 7 are hierarchized, and each layer is divided into respective time periods. In parallel.

図9においては、10s,3s,300ms,100msの時間周期で各層が並列に動作する例を示しており、それぞれの周期で環境リスクレベルを離散化した観測ステートと操作データをSOMによって離散化(量子化)した観測シンボルとがモデル学習部6に入力され、学習ステップ(学習フェーズ)で更新されたモデルパラメータが推定ステップ(推定フェーズ)で用いられる。推定ステップでは、図9に示すように、対応する階層(周期)で離散化した操作データの特徴量からドライバ内部状態の推定を行い、各層から推定状態が並列に出力される。   FIG. 9 shows an example in which each layer operates in parallel in time periods of 10 s, 3 s, 300 ms, and 100 ms, and the observation state and operation data obtained by discretizing the environmental risk level in each period are discretized by the SOM ( The quantized observation symbols are input to the model learning unit 6 and the model parameters updated in the learning step (learning phase) are used in the estimation step (estimation phase). In the estimation step, as shown in FIG. 9, the internal state of the driver is estimated from the feature quantity of the operation data discretized at the corresponding hierarchy (cycle), and the estimated state is output in parallel from each layer.

この場合、各階層毎にSOMを独立して設けても良く、或いは、複数の層で同一のSOMを共有し、サンプリング周期のみを異なる周期としても良い。また、図9においては、操作データは各層へ共通のデータを入力する例を示しているが、層によって入力される操作データを異なるものとしても良い。更に、駆動周期が同一の複数の階層に、異なる操作データを入力するようにしても良い。   In this case, SOMs may be provided independently for each layer, or the same SOM may be shared by a plurality of layers, and only the sampling periods may be different periods. FIG. 9 shows an example in which common data is input to each layer as the operation data. However, the operation data input by each layer may be different. Further, different operation data may be input to a plurality of layers having the same drive cycle.

この階層化により、運転中にドライバに生じる周期や速度が様々な状態変化を広く推定することが可能になる。例えば、一瞬の気のゆるみであれば、100ms〜1s程度の周期の階層によって学習・推定を行い、覚醒度低下や疲労の蓄積によって徐々にドライバが外部環境との対応が取れなくなっていく場合には、10s或いは更に長周期の階層により、学習・推定を行う。   By this hierarchization, it becomes possible to widely estimate state changes with various periods and speeds generated in the driver during driving. For example, in the case of momentary relaxation, learning / estimation is performed with a hierarchy of about 100 ms to 1 s, and the driver gradually becomes unable to cope with the external environment due to a decrease in alertness or accumulation of fatigue. Performs learning / estimation based on a layer of 10 s or a longer period.

このような階層化により、様々な状態変化に対応した周期で状態遷移モデルを学習させることができ、多様なドライバ内部状態の推定が可能となる。また、得られた複数の層の推定出力を後段で総合的に判断することで、より詳細な内部状態の推定を行うことが可能となる。   By such hierarchization, it is possible to learn a state transition model at a period corresponding to various state changes, and it is possible to estimate various driver internal states. In addition, it is possible to estimate the internal state in more detail by comprehensively determining the estimated outputs of the plurality of layers obtained later.

<イベントドリブン>
イベントドリブンの特徴は、外部環境リスクへのドライバ応答は運転操作として表れる、逆に言えば運転操作に表れないステート遷移は無いと仮定し、イベント発生に応じて状態推定を行うことである。すなわち、図10に示すように、前段で出力される離散化結果(SOMからの観測シンボル)を利用して運転操作やドライバの内部状態の変化を監視し、観測シンボルによるベント発生で遷移トリガを発生させる。
<Event driven>
A feature of event-driven is that a driver response to an external environmental risk appears as a driving operation. In other words, it is assumed that there is no state transition that does not appear in the driving operation, and state estimation is performed according to the occurrence of an event. That is, as shown in FIG. 10, the discretization result (observed symbol from the SOM) output in the previous stage is used to monitor changes in driving operation and the internal state of the driver, and a transition trigger is triggered by the occurrence of a vent by the observed symbol. generate.

例えば、図11に示すように、観測シンボルaへの変化点で観測シンボルaを遷移トリガとして出力した後、観測シンボルがa→bに変化したとき、c→dに変化したとき、d→eに変化したとき、そのときの観測シンボルを遷移トリガとして後段に出力する。尚、観測シンボルの変化は、SOM出力が所定のユニット差(例えば、ユニット差1以上)で変化したときとする。   For example, as shown in FIG. 11, after the observation symbol a is output as a transition trigger at the change point to the observation symbol a, the observation symbol changes from a to b, and when the observation symbol changes from c to d, d to e. When it changes to, the observation symbol at that time is output to the subsequent stage as a transition trigger. The observation symbol changes when the SOM output changes with a predetermined unit difference (for example, a unit difference of 1 or more).

また、同一の離散化結果が一定時間以上連続した場合も一つのイベント発生として捉え、図11に示すように、観測シンボルcの状態が規定以上連続した場合には、そのときの観測シンボルを出力して再トリガをかける。   Also, when the same discretization result continues for a certain time or more, it is regarded as one event occurrence, and as shown in FIG. 11, when the state of the observation symbol c continues for a predetermined time or more, the observation symbol at that time is output. And re-trigger.

このようなイベントドリブン化により、ドライバの運転操作に変化があった場合、その変化のタイミングに合わせて状態遷移モデルを遷移させることができ、遷移周期を固定した場合に比較して、より柔軟な推定器を構築することが可能になる。すなわち、ドライバの内部状態の変化が緩慢な場合は、その運転操作の変化も同様であると考えられ、逆に内部状態が頻繁に変化している場合は、運転操作に頻繁な変化が発生すると考えられる。従って、イベントドリブン化によって状態遷移モデルの遷移周期を可変にすることにより、ドライバの内部状態の変化に対して適応的に推定演算を実行することが可能となる。   With this event-driven approach, if there is a change in the driving operation of the driver, the state transition model can be transitioned according to the timing of the change, which is more flexible than when the transition period is fixed. It becomes possible to build an estimator. That is, when the change in the internal state of the driver is slow, the change in the driving operation is considered to be the same, and conversely, when the internal state is frequently changing, the driving operation is frequently changed. Conceivable. Therefore, by making the transition period of the state transition model variable by event-driven, it is possible to adaptively execute the estimation calculation with respect to the change in the internal state of the driver.

図12は、タイムドリブン及びイベントドリブンによる状態遷移確率を示すグラフであり、図12(a)は、300ms周期のタイムドリブン、図12(b)はSOM出力のユニット差1以上でのイベントドリブン、図12(c)はSOM出力のユニット差5以上でのイベントドリブンによる状態遷移確率を示している。   FIG. 12 is a graph showing the state transition probability by time driven and event driven. FIG. 12A is a time driven with a period of 300 ms, FIG. 12B is an event driven when the unit difference of SOM output is 1 or more, FIG. 12C shows the state transition probability due to event driven when the unit difference of the SOM output is 5 or more.

図12からは、タイムドリブンに比較してイベントドリブンでは自己遷移が相対的に減少し、特に高リスク側の自己遷移が低下していることがわかる。これは、一般的に、高リスクはドライバの運転操作によって回避され、低リスクは運転操作によって自己遷移されていることを示している。   From FIG. 12, it can be seen that self-transition is relatively decreased in event-driven as compared to time-driven, and in particular, self-transition on the high risk side is decreased. This generally indicates that the high risk is avoided by the driving operation of the driver, and the low risk is self-transitioned by the driving operation.

また、図13にイベントドリブンによる内部状態推定の例を示す。図13においては、EDがイベントドリブン、TDが遷移周期を一定時間に固定したタイムドリブン、Rが外部環境リスクを示しており、時間T2〜T3の期間中のH1,H2の箇所、期間T3〜T4の期間中のH3の箇所において、TDの際に見られる推定結果のハンチングがEDでは抑えられている。これは、内部状態推定の遷移周期が適応的に変化しているためと考えられる。   FIG. 13 shows an example of internal state estimation by event driving. In FIG. 13, ED is event-driven, TD is time-driven with the transition period fixed at a fixed time, R indicates external environmental risk, H1 and H2 locations in the period from time T2 to T3, and period T3- In the portion of H3 during the period of T4, the hunting of the estimation result seen during TD is suppressed in the ED. This is presumably because the transition period of the internal state estimation is adaptively changed.

尚、このイベントドリブンは、複数の階層に適用することも可能であり、また、階層によってイベントドリブンとタイムドリブンとを組み合わせることも可能である。   Note that this event driven can be applied to a plurality of hierarchies, and event driven and time driven can be combined depending on the hierarchies.

<適応的サンプリング>
適応的サンプリングは、環境リスクレベルや操作データの離散化の際、特に長周期の階層で問題となる顕著なリスクや顕著な操作のサンプリング漏れを防ぐものであり、前述のタイムドリブンやイベントドリブンと組み合わせることにより、有効性をより向上することができる。
<Adaptive sampling>
Adaptive sampling is used to prevent significant risks and significant sampling omissions that are problematic in the long-period hierarchy when discriminating environmental risk levels and operation data. By combining, the effectiveness can be further improved.

すなわち、図14(a)に示すように、周期が短い階層の場合、サンプリングの間隔が短いためリスクレベルや運転操作の時間的変化を十分に捉えることが可能であるが、一方、図14(b)に示すように、周期が長い階層では、サンプリング間隔が長くなるためにスパイク状のリスクや運転操作はその時間的変化を十分に捉えられないことがある。そこで、図14(c)に示すように、前回のサンプリングから今回のサンプリングの間に顕著な値を示した場合には、その値を代表値として今回のサンプリング値とする適応的サンプリング処理を追加することで、見かけのサンプリング周期を変更する。   That is, as shown in FIG. 14 (a), in the case of a hierarchy with a short cycle, since the sampling interval is short, it is possible to sufficiently capture the risk level and the temporal change of the driving operation. As shown in b), in a hierarchy with a long cycle, since the sampling interval becomes long, the spike-like risk and the driving operation may not be able to sufficiently capture the temporal change. Therefore, as shown in FIG. 14 (c), when a remarkable value is shown between the previous sampling and the current sampling, an adaptive sampling process in which the current sampling value is used as the representative value is added. As a result, the apparent sampling period is changed.

この適応的サンプリングを採用することにより、状態遷移周期が長い推定器のための離散化であっても、必要な情報を見落とすことが無くなり、規定のサンプリング周期をTとした場合、転舵角速度等の絶対値が大きい操作やスパイク状のリスク等を見落とすことなく、最大2Tまでの周期を拡大することが可能となる。適応的サンプリングの代表値としては、例えばリスクの場合には最大値、運転操作の場合には最大絶対値等が考えられるほか、最小値や平均値、最頻値、中央値等でも良い。   By adopting this adaptive sampling, even when discretization is performed for an estimator having a long state transition period, it is not necessary to overlook necessary information. It is possible to extend the period up to 2T without overlooking operations with a large absolute value, spike-like risks, and the like. As representative values of adaptive sampling, for example, a maximum value in the case of risk, a maximum absolute value in the case of driving operation, and the like may be considered, and a minimum value, an average value, a mode value, a median value, and the like may be used.

図15に、階層化と適応的サンプリングとを併用した例を示す。階層は前述の図9に示した4層(第0層〜第3層)で、第0層が100ms、第1層が300ms、第2層が3s、第3層が10sの遷移周期を持ち、時間0〜T1’〜T2’〜T3’〜T4’の期間における内部状態推定結果の変化を示している。100ms,300msの層は、外部環境リスクRの短周期的変化に対してドライバの内部状態が反応しているかを評価する上で重要な結果を与え、ドライバの「よそ見」や「気の緩み」に相当する状態を推定するのに適した出力を示している。一方、3s,10sの層は、外部環境リスクRの長周期的変化に対してドライバの内部状態が追従しているかを評価する上で重要な結果を与え、「ぼんやり」や「勘違い」等の疲労度や覚醒度に近いドライバ内部状態を評価するのに適した出力を示している。   FIG. 15 shows an example in which hierarchization and adaptive sampling are used together. The layers are the four layers (0th to 3rd layers) shown in FIG. 9 described above. The 0th layer has a transition period of 100 ms, the first layer is 300 ms, the second layer is 3 s, and the third layer is 10 s. The change of the internal state estimation result in the period of time 0-T1'-T2'-T3'-T4 'is shown. The 100 ms and 300 ms layers give important results in evaluating whether the internal state of the driver is responding to short-term changes in the external environmental risk R, and the driver's “look away” and “relaxation”. The output suitable for estimating the state corresponding to is shown. On the other hand, the 3s and 10s layers give important results in evaluating whether the internal state of the driver follows the long-term changes in the external environmental risk R, such as “disappointing” and “misunderstanding”. The output is suitable for evaluating the internal state of the driver that is close to the degree of fatigue and arousal.

[ドライバの意識しているリスクレベルの出力]
状態推定部7でドライバの内部状態として各ステートへの存在確率から計算されたステートの番号(或いは期待値)は、一次ローパスフィルタ(例えば、カットオフ周波数0.3Hz)でフィルタリングされた後、警報・支援部8に出力される。この出力値がドライバの意識しているリスクレベルに相当し、ドライバが対処するリスクレベルと見做せるものであり、出力値と環境リスクレベル値とを比較することにより、車両が本質的に安全或いは危険な状態にあるかを評価することが可能となる。
[Output of the driver's conscious risk level]
The state number (or expected value) calculated from the probability of existence in each state as an internal state of the driver by the state estimation unit 7 is filtered by a primary low-pass filter (for example, a cutoff frequency of 0.3 Hz), and then alarmed. Output to the support unit 8. This output value corresponds to the risk level that the driver is aware of and can be regarded as the risk level that the driver deals with. By comparing the output value with the environmental risk level value, the vehicle is inherently safe. Or it becomes possible to evaluate whether it is in a dangerous state.

すなわち、警報・支援部8での比較結果、状態推定部7で推定されたステートが環境リスク認識部2から現在得られているスカラ値の属するステートと同じであれば、ドライバは環境を正常に認識した上で正常な操作を行っていると判定される。一方、状態推定部7で推定されたステートが、環境リスク認識部2から現在得られているスカラ値の属するステートと異なる場合には、ドライバの状態が環境の要求するレベルと食い違っていると判定される。   That is, if the state estimated by the state estimating unit 7 is the same as the state to which the scalar value currently obtained from the environmental risk recognizing unit 2 belongs as a result of the comparison in the alarm / support unit 8, the driver normalizes the environment. It is determined that a normal operation is performed after recognition. On the other hand, when the state estimated by the state estimation unit 7 is different from the state to which the scalar value currently obtained from the environmental risk recognition unit 2 belongs, it is determined that the driver state is inconsistent with the level required by the environment. Is done.

例えば、推定されるステートのリスクレベルの方が環境リスク認識部2から得られるリスクレベルより低い場合には、ドライバの危険認識度が低く、安全性が損なわれる可能性が高いとして、警報・支援部8から警報・支援を行う。逆に、推定されるステートのリスクレベルの方が環境リスク認識部2から得られるリスクレベルより高い場合にも、ドライバが過度に緊張状態にあることが推定され、同様に安全性が損なわれる可能性が高いとして警報を行う。   For example, if the risk level of the estimated state is lower than the risk level obtained from the environmental risk recognition unit 2, the driver's risk recognition level is low, and there is a high possibility that safety will be impaired. Alert / support from Department 8. Conversely, if the estimated risk level of the state is higher than the risk level obtained from the environmental risk recognition unit 2, it is estimated that the driver is overly tensioned, and safety may be similarly impaired. Alarm is given as high probability.

図16に状態推定の走行試験で計測されたデータの一部を示す。S1通り(片側一車線)からA,B,Cの記号で示す間の狭路(センターライン無し)を通ってS2通り(片側一車線)に至り、更に、記号Dで示す位置からS3通り(片側一車線)までを走行したときの環境リスクレベル及び内部状態リスクレベルの変化が操作特徴量の変化と共に示されている。   FIG. 16 shows a part of data measured in the state estimation running test. From S1 street (one side lane) to the S2 street (one lane) through the narrow path (no center line) indicated by the symbols A, B, C, and further from the position indicated by symbol D, S3 street ( Changes in the environmental risk level and the internal state risk level when traveling up to one lane on one side are shown together with changes in the operation feature amount.

図16に示す試験結果では、操作特徴量とリスクレベルとの関係から、ドライバは正常な覚醒度と緊張感を保って運転しており、交差点や狭い路といったリスクレベルの高い状況に対応している様子が、推定された内部状態すなわちドライバが意識しているリスクレベルが高くなっていることから確認できる。特に、グラフ中に示した記号A,B,C,Dの時点は、特徴的な高リスクイベントであり、主として交差点と狭路である。そのほかのリスクレベルが高くなっている時点は、主として歩行者や自転車・対向車等と遭遇した部分であり、この部分でも、同様の推定結果が得られていることが分かる。   In the test results shown in FIG. 16, the driver is driving with a normal arousal level and a sense of tension based on the relationship between the operation feature quantity and the risk level, and corresponds to a high risk level situation such as an intersection or a narrow road. Can be confirmed from the estimated internal state, that is, the risk level perceived by the driver is high. In particular, the points A, B, C, and D shown in the graph are characteristic high-risk events, mainly intersections and narrow streets. When the other risk level is high, it is mainly a portion where a pedestrian, a bicycle, an oncoming vehicle, or the like is encountered, and it can be seen that the same estimation result is obtained also in this portion.

以上のように、本実施の形態における安全運転支援システムは、車両に搭載されているアクセル開度・ブレーキ圧・舵角等のデータを入力として、走行環境のリスクレベルとの対応関係を確率的状態遷移モデルで近似してドライバの運転操作の特徴を学習し、走行環境のリスクレベルとドライバの内部状態との対応関係を評価する。   As described above, the safe driving support system according to the present embodiment uses the data such as the accelerator opening, the brake pressure, and the steering angle that are mounted on the vehicle as inputs, and the correspondence relationship with the risk level of the driving environment is stochastic. The characteristics of the driving operation of the driver are learned by approximation with the state transition model, and the correspondence between the risk level of the driving environment and the internal state of the driver is evaluated.

このシステムは、外部環境のリスクレベルだけで警報や操作アシストその他の介入の要否を判断するのではなく、ドライバが外部環境を正しく把握しているか否かを評価指標とするものであり、これにより、安全システムの介入の要否を走行環境の危険度のみで判断するために起こる不必要な介入による違和感・煩わしさを解消し、必要なときに適切な支援を提供することができ、個々のドライバに適応した高度な安全システムを構築することが可能となる。   This system does not judge whether alarms, operation assistance, or other interventions are necessary based only on the risk level of the external environment, but rather uses the evaluation index as to whether the driver correctly understands the external environment. This eliminates the discomfort and annoyance caused by unnecessary interventions that occur when judging whether or not a safety system intervention is necessary based only on the risk of the driving environment, and provides appropriate support when necessary. It is possible to construct an advanced safety system adapted to the driver.

しかも、本システムでは、状態推定のための主たるデータとして運転操作データを用いるため、各種車両制御用に既に取り付けられているハンドルやアクセル・ブレーキの操作を観測するセンサを流用することができ、生体計測装置のように追加コストの必要が無くなる。これらのセンサは、ノイズに対して比較的頑強であることで、状態推定を安定して行なうことができ、さらには、ドライバの通常の運転操作を妨げたり姿勢を拘束することが無いため、状態推定に悪影響を与える不必要な干渉を排除することができる。   Moreover, in this system, since the driving operation data is used as the main data for state estimation, it is possible to divert the sensors for observing the operation of the steering wheel and accelerator / brake already installed for various vehicle control, There is no need for additional cost as in the case of a measuring device. Since these sensors are relatively robust against noise, it is possible to stably estimate the state, and further, the normal driving operation of the driver is not hindered and the posture is not restrained. Unnecessary interference that adversely affects the estimation can be eliminated.

また、本システムでは、車外環境と運転操作との対応関係を規範モデルとしているため、予め条件やテンプレートを決める必要が無くなり、学習によってモデルを更新することで、ドライバ毎の個人差や、ドライバ毎の走行環境の差に対応することができ、走行を重ねることで状態推定精度を向上させていくことができる。   In addition, in this system, since the correspondence relationship between the environment outside the vehicle and the driving operation is a normative model, there is no need to determine conditions and templates in advance, and by updating the model through learning, individual differences for each driver and each driver Thus, the state estimation accuracy can be improved by repeated traveling.

この場合、規範となる学習モデルを確率的状態遷移モデルとして、バラツキのある操作データを統計的に処理することで、ノイズや外乱に大きく影響されることが無く、環境とドライバとの本質的な関連性を学習することができる。また、学習時に確率的状態遷移モデルのステート遷移を、外部環境から認識されるリスクレベルの遷移と対応させることで、推定された内部状態がどのリスクレベルに対応するのかを推定することができる。   In this case, the learning model that becomes the norm is used as a probabilistic state transition model, and the operation data with variations is statistically processed. Relevance can be learned. In addition, it is possible to estimate which risk level the estimated internal state corresponds to by matching the state transition of the probabilistic state transition model with the transition of the risk level recognized from the external environment at the time of learning.

すなわち、推定された内部状態の対応するリスクレベルが認識された外部環境のリスクレベルより低いときには、たとえ外部環境のリスクレベルが高い場合でなくとも、ドライバが安全に運転できる状態に無いとして、警報・操作支援を行うことができる。一方、推定された内部状態の対応するリスクレベルが認識された外部環境のリスクレベルと同じかより高いときには、たとえ外部環境のリスクレベルが高いとしても、ドライバはその環境に対応できる状態にあると推定できるので、無駄な警報や不必要な操作支援を抑制することができる。   In other words, when the risk level corresponding to the estimated internal state is lower than the recognized external environment risk level, the driver is not in a safe driving state even if the external environment risk level is high.・ Operation support can be performed. On the other hand, when the risk level corresponding to the estimated internal state is the same or higher than the recognized external environment risk level, the driver is in a state that can handle the environment even if the external environment risk level is high. Since it can be estimated, useless warnings and unnecessary operation support can be suppressed.

しかも、ドライバの内部状態を推定する際に、状態遷移モデルによる推定周期を固定するのではなく、異なる時間周期での階層化、イベント発生の有無、適応的サンプリングといった周期で内部状態を推定することにより、運転中にドライバに生じる周期や速度が様々な状態変化を広く推定することが可能になり、高度なシステムを構築することが可能となる。   In addition, when estimating the internal state of the driver, the internal state is estimated at a period such as hierarchization at different time periods, the occurrence of an event, or adaptive sampling, rather than fixing the estimation period based on the state transition model. Thus, it becomes possible to widely estimate various state changes in the period and speed generated in the driver during driving, and it becomes possible to construct an advanced system.

すなわち、階層化による推定では、様々な状態変化を、対応する階層の状態遷移モデルに学習させることができ、多様なドライバ内部状態の推定が可能となる。また、得られた複数の層の推定出力を後段で総合的に判断することで、より詳細な内部状態の推定を行うことが可能となる。   That is, in the estimation by hierarchization, various state changes can be learned by the state transition model of the corresponding hierarchy, and various driver internal states can be estimated. In addition, it is possible to estimate the internal state in more detail by comprehensively determining the estimated outputs of the plurality of layers obtained later.

また、ドライバの内部状態の変化が緩慢な場合には、その運転操作の変化も同様であると考えられ、逆に内部状態が頻繁に変化している場合は運転操作に頻繁な変化が発生すると考えられることから、イベントドリブン化により、ドライバの運転操作に変化があった場合に、その変化のタイミングに合わせて状態遷移モデルを遷移させることができ、遷移周期を固定した場合よりもドライバの内部状態の変化に適応的に推定演算を実行することが可能となる。   In addition, when the change in the internal state of the driver is slow, the change in the driving operation is also considered to be the same, and conversely, if the internal state changes frequently, a frequent change occurs in the driving operation. As a result, if there is a change in the driving operation of the driver due to event driven, the state transition model can be transitioned according to the timing of the change, and the internal state of the driver can be improved compared to when the transition period is fixed. It is possible to execute the estimation calculation adaptively to the change of the state.

更に、適応的サンプリングによって、長い駆動周期を持つ確率的状態遷移モデルの層に入力されるリスクレベルや操作データの離散化の際に、観測すべき特徴的な入力値をサンプリング漏れすることなく学習・推定を行うことが可能になる。   In addition, adaptive sampling enables learning without missing sampling characteristic input values to be observed when discriminating the risk level and operation data input to the layer of the probabilistic state transition model having a long drive cycle.・ Estimation can be performed.

安全運転支援システムの基本構成図Basic configuration of safe driving support system 安全運転支援システムの学習フェーズ及び推定フェーズを示す説明図Explanatory drawing showing the learning phase and estimation phase of the safe driving support system 操作特徴量と自己組織化マップのネットワークを示す説明図Explanatory diagram showing the network of operation features and self-organizing maps 確率的状態遷移を示す説明図Explanatory diagram showing stochastic state transition 学習後の状態遷移確率を示すグラフGraph showing state transition probability after learning 学習後の出力確率を示すグラフGraph showing output probability after learning 最尤系列による状態推定の説明図Illustration of state estimation by maximum likelihood sequence 前向きアルゴリズムによる状態推定の説明図Illustration of state estimation by forward algorithm タイムドリブンによる状態推定を示す説明図Explanatory diagram showing state estimation by time-driven イベントドリブンによる状態推定を示す説明図Explanatory diagram showing state estimation by event driven イベントドリブンによる遷移トリガの発生を示す説明図Explanatory diagram showing the occurrence of transition trigger by event driven タイムドリブン及びイベントドリブンによる状態遷移確率を示すグラフGraph showing state transition probability by time-driven and event-driven タイムドリブン及びイベントドリブンによる内部状態推定例を示すグラフGraph showing examples of internal state estimation by time-driven and event-driven 適応的サンプリングの説明図Illustration of adaptive sampling タイムドリブンと適応的サンプリングとを併用した内部状態推定例を示すグラフA graph showing an example of internal state estimation using both time-driven and adaptive sampling 推定試験結果を示す説明図Explanatory diagram showing estimation test results

符号の説明Explanation of symbols

1 安全運転支援システム
2 車外環境認識部
3 環境リスク離散化部
5 操作特徴量離散化部
6 モデル学習部
7 状態推定部
8 警報・支援部
DESCRIPTION OF SYMBOLS 1 Safe driving assistance system 2 External environment recognition part 3 Environmental risk discretization part 5 Operation feature discretization part 6 Model learning part 7 State estimation part 8 Warning / support part

Claims (9)

移動体を運転するドライバの状態を推定して安全運転のための支援を行う安全運転支援システムであって、
上記移動体の外部環境をセンシングして特徴量を抽出し、該特徴量に基づいて外部環境に含まれる環境リスクを認識する環境リスク認識部と、
上記環境リスクとドライバの運転操作との対応を確率的状態遷移モデルを用いて学習し、学習モデルを構築するモデル学習部と、
上記学習モデルに基づいて、走行中のドライバの運転操作の反映としての車両挙動データを含む運転操作データからドライバのリスク認識状態をドライバの内部状態として推定する状態推定部と、
上記環境リスクと上記ドライバの内部状態とを比較し、上記移動体の安全運転に係る支援情報を取得する運転支援部とを備え、
更に、上記確率的状態遷移モデルへのデータ入力を複数の異なる周期に設定し、上記確率的状態遷移モデルによる推定周期を適応的に可変するスケジュール制御部を備えることを特徴とする安全運転支援システム。
A safe driving support system that estimates the state of a driver driving a moving body and provides support for safe driving,
An environmental risk recognition unit that senses an external environment of the mobile body to extract a feature amount, and recognizes an environmental risk included in the external environment based on the feature amount;
A model learning unit that learns the correspondence between the environmental risk and the driving operation of the driver using a probabilistic state transition model, and constructs a learning model;
Based on the learning model, a state estimation unit that estimates a driver's risk recognition state as a driver's internal state from driving operation data including vehicle behavior data as a reflection of the driving operation of the driver while traveling,
A driving support unit that compares the environmental risk with the internal state of the driver and obtains support information related to safe driving of the mobile body,
Furthermore, a safe driving support system comprising a schedule control unit that sets data input to the stochastic state transition model in a plurality of different periods and adaptively varies an estimation period based on the stochastic state transition model .
上記スケジュール制御部は、
上記確率的状態遷移モデルを時間的に階層化して駆動することを特徴とする請求項1記載の安全運転支援システム。
The schedule control unit
2. The safe driving support system according to claim 1, wherein the stochastic state transition model is driven in a hierarchical manner.
上記階層化を、上記運転操作データの種別に応じて実施することを特徴とする請求項2記載の安全運転支援システム。   The safe driving support system according to claim 2, wherein the hierarchization is performed according to a type of the driving operation data. 上記スケジュール制御部は、
モデル入力側のイベント発生に応じて遷移トリガを出力し、該遷移トリガによって上記確率的状態遷移モデルを駆動することを特徴とする請求項1〜3の何れか一に記載の安全運転支援システム。
The schedule control unit
The safe driving support system according to any one of claims 1 to 3, wherein a transition trigger is output in response to an event occurrence on the model input side, and the stochastic state transition model is driven by the transition trigger.
上記確率的状態遷移モデルへの入力データが変化したとき、イベント発生として上記遷移トリガを出力することを特徴とする請求項4記載の安全運転支援システム。   5. The safe driving support system according to claim 4, wherein when the input data to the stochastic state transition model changes, the transition trigger is output as an event occurrence. 上記確率的状態遷移モデルへの入力データが変化しない状態が一定時間以上継続したとき、イベント発生として上記遷移トリガを出力することを特徴とする請求項4記載の安全運転支援システム。   5. The safe driving support system according to claim 4, wherein when the state in which the input data to the stochastic state transition model does not change continues for a certain time or more, the transition trigger is output as an event occurrence. 上記スケジュール制御部は、
上記確率的状態遷移モデルへのデータ入力周期を、サンプリング間の代表値で可変することを特徴とする請求項1〜6の何れか一に記載の安全運転支援システム。
The schedule control unit
The safe driving support system according to any one of claims 1 to 6, wherein a data input period to the probabilistic state transition model is changed by a representative value between samplings.
上記代表値として、サンプリング間の最大値或いは最小値を用いることを特徴とする請求項7記載の安全運転支援システム。   8. The safe driving support system according to claim 7, wherein a maximum value or a minimum value between samplings is used as the representative value. 上記代表値として、サンプリング間の平均値或いは最頻値或いは中央値を用いることを特徴とする請求項7記載の安全運転支援システム。   8. The safe driving support system according to claim 7, wherein an average value, a mode value, or a median value between samplings is used as the representative value.
JP2008113012A 2008-04-23 2008-04-23 Safe driving support system Expired - Fee Related JP5161643B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2008113012A JP5161643B2 (en) 2008-04-23 2008-04-23 Safe driving support system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008113012A JP5161643B2 (en) 2008-04-23 2008-04-23 Safe driving support system

Publications (2)

Publication Number Publication Date
JP2009262702A true JP2009262702A (en) 2009-11-12
JP5161643B2 JP5161643B2 (en) 2013-03-13

Family

ID=41389117

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008113012A Expired - Fee Related JP5161643B2 (en) 2008-04-23 2008-04-23 Safe driving support system

Country Status (1)

Country Link
JP (1) JP5161643B2 (en)

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013022211A (en) * 2011-07-20 2013-02-04 Nissan Motor Co Ltd Device for estimating driver's fatigue level
JP2015056058A (en) * 2013-09-12 2015-03-23 株式会社デンソー Driving support device
US9132775B2 (en) 2010-11-30 2015-09-15 Toyota Jidosha Kabushiki Kaisha Mobile object target state determination device and program
US9405727B2 (en) 2011-02-23 2016-08-02 Toyota Jidosha Kabushiki Kaisha Driving support device, driving support method, and driving support program
CN107444398A (en) * 2016-04-14 2017-12-08 丰田自动车株式会社 Drive auxiliary device
CN107618512A (en) * 2017-08-23 2018-01-23 清华大学 Driving behavior safe evaluation method based on people's car environment multi-data source
JP2018032229A (en) * 2016-08-25 2018-03-01 三菱電機ビルテクノサービス株式会社 Equipment inspection order setup apparatus and equipment inspection order setup method
CN108423005A (en) * 2017-02-15 2018-08-21 福特全球技术公司 The generation of the Controlling model based on feedback for autonomous vehicle
JP2018169700A (en) * 2017-03-29 2018-11-01 マツダ株式会社 Vehicle driving support system
JP2018169702A (en) * 2017-03-29 2018-11-01 マツダ株式会社 Vehicle driving support system
CN108773373A (en) * 2016-09-14 2018-11-09 北京百度网讯科技有限公司 Method and apparatus for operating automatic driving vehicle
JP2019507443A (en) * 2016-03-01 2019-03-14 ヴァレオ、コンフォート、アンド、ドライビング、アシスタンスValeo Comfort And Driving Assistance Personalization apparatus and method for monitoring motor vehicle drivers
JP2019212276A (en) * 2018-06-06 2019-12-12 緯創資通股▲ふん▼有限公司Wistron Corporation Method, processing device and system for driving prediction
CN111295699A (en) * 2017-08-31 2020-06-16 松下知识产权经营株式会社 Assistance method, assistance system using the assistance method, and assistance device
JP2020140309A (en) * 2019-02-27 2020-09-03 トヨタ自動車株式会社 Evaluation device
JP2020140310A (en) * 2019-02-27 2020-09-03 トヨタ自動車株式会社 Evaluation device
CN111717217A (en) * 2020-06-30 2020-09-29 重庆大学 Driver intention identification method based on probability correction

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06251159A (en) * 1993-03-01 1994-09-09 Nippon Telegr & Teleph Corp <Ntt> Operation recognizing device
JPH0728862A (en) * 1993-07-14 1995-01-31 Toshiba Corp Hierarchized state transition model
JPH0990983A (en) * 1995-09-26 1997-04-04 Nippon Telegr & Teleph Corp <Ntt> Voice recognition method
JPH09185716A (en) * 1995-11-02 1997-07-15 Yamaha Corp Operation judging method and its device
JPH1199848A (en) * 1997-09-30 1999-04-13 Nissan Motor Co Ltd Vehicle action estimating device
JP2004164582A (en) * 2002-09-24 2004-06-10 Ricoh Co Ltd Information processing apparatus and program generation method
JP2004206207A (en) * 2002-12-24 2004-07-22 Toyota Central Res & Dev Lab Inc Danger sensing degree estimation device, safe driving evaluation device, and alarm device
JP2004212641A (en) * 2002-12-27 2004-07-29 Toshiba Corp Voice input system and terminal device equipped with voice input system
JP2004341340A (en) * 2003-05-16 2004-12-02 Toshiba Tec Corp Speaker recognition device

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06251159A (en) * 1993-03-01 1994-09-09 Nippon Telegr & Teleph Corp <Ntt> Operation recognizing device
JPH0728862A (en) * 1993-07-14 1995-01-31 Toshiba Corp Hierarchized state transition model
JPH0990983A (en) * 1995-09-26 1997-04-04 Nippon Telegr & Teleph Corp <Ntt> Voice recognition method
JPH09185716A (en) * 1995-11-02 1997-07-15 Yamaha Corp Operation judging method and its device
JPH1199848A (en) * 1997-09-30 1999-04-13 Nissan Motor Co Ltd Vehicle action estimating device
JP2004164582A (en) * 2002-09-24 2004-06-10 Ricoh Co Ltd Information processing apparatus and program generation method
JP2004206207A (en) * 2002-12-24 2004-07-22 Toyota Central Res & Dev Lab Inc Danger sensing degree estimation device, safe driving evaluation device, and alarm device
JP2004212641A (en) * 2002-12-27 2004-07-29 Toshiba Corp Voice input system and terminal device equipped with voice input system
JP2004341340A (en) * 2003-05-16 2004-12-02 Toshiba Tec Corp Speaker recognition device

Cited By (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9132775B2 (en) 2010-11-30 2015-09-15 Toyota Jidosha Kabushiki Kaisha Mobile object target state determination device and program
US9405727B2 (en) 2011-02-23 2016-08-02 Toyota Jidosha Kabushiki Kaisha Driving support device, driving support method, and driving support program
JP2013022211A (en) * 2011-07-20 2013-02-04 Nissan Motor Co Ltd Device for estimating driver's fatigue level
JP2015056058A (en) * 2013-09-12 2015-03-23 株式会社デンソー Driving support device
JP2019507443A (en) * 2016-03-01 2019-03-14 ヴァレオ、コンフォート、アンド、ドライビング、アシスタンスValeo Comfort And Driving Assistance Personalization apparatus and method for monitoring motor vehicle drivers
CN107444398A (en) * 2016-04-14 2017-12-08 丰田自动车株式会社 Drive auxiliary device
CN107444398B (en) * 2016-04-14 2019-10-22 丰田自动车株式会社 Drive auxiliary device
JP2018032229A (en) * 2016-08-25 2018-03-01 三菱電機ビルテクノサービス株式会社 Equipment inspection order setup apparatus and equipment inspection order setup method
CN108773373A (en) * 2016-09-14 2018-11-09 北京百度网讯科技有限公司 Method and apparatus for operating automatic driving vehicle
CN108773373B (en) * 2016-09-14 2020-04-24 北京百度网讯科技有限公司 Method and device for operating an autonomous vehicle
CN108423005B (en) * 2017-02-15 2022-12-27 福特全球技术公司 Generation of feedback-based control model for autonomous vehicles
CN108423005A (en) * 2017-02-15 2018-08-21 福特全球技术公司 The generation of the Controlling model based on feedback for autonomous vehicle
JP2018169702A (en) * 2017-03-29 2018-11-01 マツダ株式会社 Vehicle driving support system
JP2018169700A (en) * 2017-03-29 2018-11-01 マツダ株式会社 Vehicle driving support system
CN107618512A (en) * 2017-08-23 2018-01-23 清华大学 Driving behavior safe evaluation method based on people's car environment multi-data source
CN111295699A (en) * 2017-08-31 2020-06-16 松下知识产权经营株式会社 Assistance method, assistance system using the assistance method, and assistance device
CN111295699B (en) * 2017-08-31 2023-03-21 松下知识产权经营株式会社 Assistance method, assistance system using the assistance method, and assistance device
JP2019212276A (en) * 2018-06-06 2019-12-12 緯創資通股▲ふん▼有限公司Wistron Corporation Method, processing device and system for driving prediction
JP2020140309A (en) * 2019-02-27 2020-09-03 トヨタ自動車株式会社 Evaluation device
JP7074698B2 (en) 2019-02-27 2022-05-24 トヨタ自動車株式会社 Evaluation device
JP7173702B2 (en) 2019-02-27 2022-11-16 トヨタ自動車株式会社 Evaluation device
CN111626537A (en) * 2019-02-27 2020-09-04 丰田自动车株式会社 Evaluation device
JP2020140310A (en) * 2019-02-27 2020-09-03 トヨタ自動車株式会社 Evaluation device
CN111626537B (en) * 2019-02-27 2023-09-12 丰田自动车株式会社 Evaluation device
CN111717217A (en) * 2020-06-30 2020-09-29 重庆大学 Driver intention identification method based on probability correction
CN111717217B (en) * 2020-06-30 2022-11-08 重庆大学 Driver intention identification method based on probability correction

Also Published As

Publication number Publication date
JP5161643B2 (en) 2013-03-13

Similar Documents

Publication Publication Date Title
JP5161643B2 (en) Safe driving support system
CN111741884B (en) Traffic distress and road rage detection method
US11345234B2 (en) Method and apparatus for detecting status of vehicle occupant
JP2009073465A (en) Safe driving support system
Sigari et al. A driver face monitoring system for fatigue and distraction detection
CN104730949B (en) Affective user interface in autonomous vehicle
JP5005492B2 (en) Safe driving support system
US20210357701A1 (en) Evaluation device, action control device, evaluation method, and evaluation program
US7839292B2 (en) Real-time driving danger level prediction
US20190370580A1 (en) Driver monitoring apparatus, driver monitoring method, learning apparatus, and learning method
JP6219952B2 (en) How to determine the operating status of a driver
US7394393B2 (en) Adaptive driver workload estimator
EP3535646A1 (en) System and method for driver distraction determination
KR101276770B1 (en) Advanced driver assistance system for safety driving using driver adaptive irregular behavior detection
Wu et al. Reasoning-based framework for driving safety monitoring using driving event recognition
JP7212607B2 (en) Machine-learning-based driver abnormality detection method and apparatus using vehicle CAN BUS signal
US9111400B2 (en) System for detecting abnormal driving behavior
US20240029467A1 (en) Modular Predictions For Complex Human Behaviors
Sharma et al. Drowsiness warning system using artificial intelligence
JP2009096365A (en) Risk recognition system
Gite et al. Context–aware assistive driving: an overview of techniques for mitigating the risks of driver in real-time driving environment
Mittal et al. Driver drowsiness detection using machine learning and image processing
Ponomarev et al. Adaptation and personalization in driver assistance systems
Wu et al. Fuzzy logic based driving behavior monitoring using hidden markov models
Amira et al. Driver drowsiness detection and tracking based on YOLO with Haar cascades and ERNN

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20110324

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20121031

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: 20121120

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20121214

R150 Certificate of patent or registration of utility model

Ref document number: 5161643

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

Free format text: PAYMENT UNTIL: 20151221

Year of fee payment: 3

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees