JP2016042298A - Accident information calculation apparatus and program - Google Patents

Accident information calculation apparatus and program Download PDF

Info

Publication number
JP2016042298A
JP2016042298A JP2014166068A JP2014166068A JP2016042298A JP 2016042298 A JP2016042298 A JP 2016042298A JP 2014166068 A JP2014166068 A JP 2014166068A JP 2014166068 A JP2014166068 A JP 2014166068A JP 2016042298 A JP2016042298 A JP 2016042298A
Authority
JP
Japan
Prior art keywords
accident
value
probability
feature
calculated
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
JP2014166068A
Other languages
Japanese (ja)
Other versions
JP6365107B2 (en
Inventor
優司 伊藤
Yuji Ito
優司 伊藤
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.)
Toyota Central R&D Labs Inc
Original Assignee
Toyota Central R&D Labs Inc
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 Toyota Central R&D Labs Inc filed Critical Toyota Central R&D Labs Inc
Priority to JP2014166068A priority Critical patent/JP6365107B2/en
Publication of JP2016042298A publication Critical patent/JP2016042298A/en
Application granted granted Critical
Publication of JP6365107B2 publication Critical patent/JP6365107B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

PROBLEM TO BE SOLVED: To make it possible to calculate the degree of occurrence of accidents using an obtainable amount of statistics.SOLUTION: A statistical-information acquisition section 30 acquires driving information and accident statistical information. A calculation section 32 calculates a parameter used in a predefined function for approximating an accident occurrence probability for a value of each feature quantity, on the basis of a constraint condition in relation to the probability that the value of each feature quantity is obtained under an accident occurrence condition as well as the driving information and the accident statistical information, and calculates a function using the calculated parameter. A feature-quantity acquisition section 28 acquires the value of each of the plurality of feature quantities related to the driving behavior of a vehicle driven by a driver of interest. A driving-behavior evaluation section 34 evaluates, as an accident-free probability, the product of the accident-free probabilities which are calculated to correspond to a plurality of types of accidents, respectively for the driving behavior of the vehicle driven by the driver of interest on the basis of the acquired value of each of the plurality of feature quantities and the function calculated for each of the plurality of types of accidents calculated by the calculation section 32.SELECTED DRAWING: Figure 2

Description

本発明は、事故情報算出装置、及びプログラムに係り、特に、事故が発生する度合いを算出する事故情報算出装置、及びプログラムに関する。   The present invention relates to an accident information calculation apparatus and program, and more particularly, to an accident information calculation apparatus and program for calculating the degree of occurrence of an accident.

従来、交通事故に結びつきやすい要因を考慮した上で運転を診断するため、道路上の交通規則を取得し、その規則と運転挙動に関する情報を用いてドライバの運転特性を診断する方法が知られている(特許文献1)。   Conventionally, in order to diagnose driving considering factors that are likely to lead to traffic accidents, there is a known method for obtaining traffic rules on the road and diagnosing driver driving characteristics using information on the rules and driving behavior. (Patent Document 1).

また、運転者の過去の運転履歴(車両の速度、ブレーキ及びアクセル操作の頻度、車間距離、車線切り替えの頻度、追い越しの頻度、クラクション鳴動の頻度等)、運転状況(天気、時間帯、気温、湿度、カーブ走行か直線走行かなどの区別、高速道路か一般道路かなどの道路種別等)、運転者の体調の情報(注意力、体温、血圧等)に基づき、危険度を推定する方法が知られている(特許文献2)。   In addition, the driver's past driving history (vehicle speed, brake and accelerator operation frequency, inter-vehicle distance, lane switching frequency, overtaking frequency, horn ringing frequency, etc.), driving status (weather, time zone, temperature, There is a method to estimate the degree of risk based on humidity, distinction between curve driving and straight driving, road type such as highway or ordinary road, etc.) and physical condition information of the driver (attention, body temperature, blood pressure, etc.) Known (Patent Document 2).

また、各走行状況において、基準となる運転挙動データと自車の運転挙動データとの乖離度に基づくことで、運転の危険性を診断する方法が知られている(特許文献3)。   Also, a method for diagnosing driving risk based on the degree of deviation between the driving behavior data serving as a reference and the driving behavior data of the host vehicle in each driving situation is known (Patent Document 3).

また、過去の事故統計情報と、各事故の発生に対する事故要因の確率テーブルを用いて、ベイジアンネットワークにより自車の運転状態に対する事故の発生確率を算出する方法が知られている(特許文献4)。   Also, a method of calculating the probability of occurrence of an accident with respect to the driving state of the vehicle by a Bayesian network using past accident statistical information and an accident factor probability table for the occurrence of each accident is known (Patent Document 4). .

また、評価対象者を含む複数運転者の運転操作情報から得た統計的な運転情報と、評価対象者の運転情報とを用いて、他の運転者との相対的な評価を行うことで、運転を診断する方法が知られている(特許文献5)。   In addition, by using the statistical driving information obtained from the driving operation information of a plurality of drivers including the evaluation subject and the driving information of the evaluation subject, by performing a relative evaluation with other drivers, A method for diagnosing driving is known (Patent Document 5).

特開2006−209455号公報JP 2006-209455 A 特開2006−277058号公報JP 2006-277058 A 特開2006−343904号公報JP 2006-343904 A 特開2009−282873号公報JP 2009-282873 A 特開2012−247854号公報JP 2012-247854 A

しかし、上記特許文献4の手法では、過去に発生した事故に対する事故要因の確率分布を算出する際に困難が生じる。例えば、交通事故発生時の天候等、観測可能な情報は確率分布化が可能であるが、事故発生直前の運転挙動は事故要因に該当し得るにも関わらず観測困難である。さらに、速度超過の有無等の離散情報ではなく、加速や脇見時間等の連続値の場合は大量の事例データが必要となり、確率分布の取得が更に困難であるという問題がある。   However, the method of Patent Document 4 has difficulty in calculating the probability distribution of accident factors for accidents that have occurred in the past. For example, observable information such as the weather at the time of a traffic accident can be probability-distributed, but the driving behavior immediately before the occurrence of the accident is difficult to observe although it may correspond to the accident factor. Furthermore, in the case of continuous values such as acceleration and aside time rather than discrete information such as the presence or absence of excess speed, a large amount of case data is required, which makes it more difficult to obtain a probability distribution.

本発明では、上記問題点を解決するために成されたものであり、取得可能な統計量を用いて、事故が発生する度合いを算出することができる事故情報算出装置、及びプログラムを提供することを目的とする。   In the present invention, there is provided an accident information calculation device and a program which are made to solve the above-described problems and which can calculate the degree of occurrence of an accident using an acquirable statistic. With the goal.

上記目的を達成するために、本発明の事故情報算出装置は、車両の運転行動に関する特徴量の各値が得られる度合いを表す運転情報と、前記特徴量に起因して生じる事故が発生する度合いを表す事故統計情報とを取得する統計情報取得手段と、前記統計情報取得手段によって取得された前記運転情報及び前記事故統計情報と、前記事故が発生した条件下での前記特徴量の各値が得られる度合いに関する制約条件とに基づいて、前記特徴量の各値に対して前記事故が発生する度合いを算出する算出手段と、を含んで構成されている。   In order to achieve the above object, the accident information calculation apparatus according to the present invention provides driving information indicating the degree to which each value of the characteristic amount related to the driving behavior of the vehicle is obtained, and the degree of occurrence of an accident caused by the characteristic amount. Statistical information acquisition means for acquiring accident statistical information representing, the driving information and the accident statistical information acquired by the statistical information acquisition means, and each value of the feature value under the condition where the accident occurred Calculation means for calculating the degree of occurrence of the accident for each value of the feature amount based on the constraint condition relating to the degree to be obtained.

また、本発明のプログラムは、コンピュータを、車両の運転行動に関する特徴量の各値が得られる度合いを表す運転情報と、前記特徴量に起因して生じる事故が発生する度合いを表す事故統計情報とを取得する統計情報取得手段、及び前記統計情報取得手段によって取得された前記運転情報及び前記事故統計情報と、前記事故が発生した条件下での前記特徴量の各値が得られる度合いに関する制約条件とに基づいて、前記特徴量の各値に対して前記事故が発生する度合いを算出する算出手段として機能させるためのプログラムである。   In addition, the program of the present invention, the computer, driving information representing the degree to which each value of the feature value related to the driving behavior of the vehicle is obtained, accident statistics information representing the degree of occurrence of the accident caused by the feature value, A statistical information acquiring means for acquiring the driving information and the accident statistical information acquired by the statistical information acquiring means, and a constraint on the degree to which each value of the feature value is obtained under the condition where the accident occurred And a program for functioning as a calculation means for calculating the degree of occurrence of the accident for each value of the feature amount.

本発明によれば、統計情報取得手段によって、車両の運転行動に関する特徴量の各値が得られる度合いを表す運転情報と、特徴量に起因して生じる事故が発生する度合いを表す事故統計情報とを取得する。   According to the present invention, the driving information representing the degree to which each value of the characteristic amount related to the driving behavior of the vehicle is obtained by the statistical information acquisition means, and the accident statistical information representing the degree of occurrence of an accident caused by the characteristic amount, To get.

そして、算出手段によって、統計情報取得手段によって取得された運転情報及び事故統計情報と、事故が発生した条件下での特徴量の各値が得られる度合いに関する制約条件とに基づいて、特徴量の各値に対して事故が発生する度合いを算出する。   Then, based on the driving information and accident statistical information acquired by the statistical information acquisition unit and the constraint condition regarding the degree of each characteristic value obtained under the condition where the accident occurred, the calculation unit calculates the feature amount. The degree of accident occurrence is calculated for each value.

このように、運転情報及び事故統計情報と、事故が発生した条件下での特徴量の各値が得られる度合いに関する制約条件とに基づいて、特徴量の各値に対して事故が発生する度合いを算出することにより、取得可能な統計量を用いて、事故が発生する度合いを算出することができる。   As described above, the degree of occurrence of an accident for each value of the feature amount based on the driving information and the accident statistical information, and the constraint condition regarding the degree to which each value of the feature amount is obtained under the condition where the accident occurred. By calculating, the degree of occurrence of an accident can be calculated using the statistics that can be acquired.

また、本発明前記算出手段は、前記統計情報取得手段によって取得された前記運転情報及び前記事故統計情報と、前記事故が発生した条件下での前記特徴量の各値が得られる度合いに関する制約条件とに基づいて、前記特徴量の各値に対して前記事故が発生する度合いを近似した予め定められた関数で用いるパラメータを算出し、前記算出されたパラメータを用いた前記関数を、前記特徴量の各値に対して前記事故が発生する度合いとして算出するようにすることができる。   Further, the calculation means according to the present invention includes a restriction condition regarding the driving information and the accident statistical information acquired by the statistical information acquisition means, and a degree to which each value of the feature value under the condition where the accident occurs is obtained. Based on the above, a parameter used in a predetermined function approximating the degree of occurrence of the accident for each value of the feature value is calculated, and the function using the calculated parameter is calculated as the feature value. It is possible to calculate the degree of occurrence of the accident for each value.

また、前記算出手段は、前記統計情報取得手段によって取得された前記運転情報及び前記事故統計情報と、前記事故が発生した条件下での前記特徴量の各値が得られる度合いに関する制約条件とに基づいて、ベイズの定理を用いて、前記パラメータを算出するようにすることができる。   Further, the calculation means includes the driving information and the accident statistical information acquired by the statistical information acquisition means, and a constraint condition regarding a degree to which each value of the feature value is obtained under the condition where the accident occurs. Based on this, the parameter can be calculated using Bayes' theorem.

また、本発明は、対象ドライバの車両の運転行動に関する特徴量の値を取得する特徴量取得手段と、前記特徴量取得手段によって取得された前記特徴量の値、及び前記算出手段によって算出された前記特徴量の各値に対して前記事故が発生する度合いに基づいて、前記対象ドライバの車両の運転行動について、前記事故が発生する度合い又は前記事故が発生しない度合いを評価する運転行動評価手段と、を更に含むことができる。   Further, the present invention provides a feature amount acquisition unit that acquires a feature amount value related to the driving behavior of the vehicle of the target driver, the feature amount value acquired by the feature amount acquisition unit, and the calculation unit. Driving behavior evaluation means for evaluating the degree of occurrence of the accident or the degree of no occurrence of the accident with respect to the driving behavior of the vehicle of the target driver based on the degree of occurrence of the accident for each value of the feature amount; Can further be included.

また、前記統計情報取得手段は、車両の運転行動に関する特徴量の各値が得られる確率を表す前記運転情報と、前記特徴量に起因して生じる事故が発生する確率を表す前記事故統計情報とを取得し、前記算出手段は、前記統計情報取得手段によって取得された前記運転情報及び前記事故統計情報と、前記事故が発生した条件下での前記特徴量の各値が得られる確率に関する前記制約条件とに基づいて、前記特徴量の各値に対して前記事故が発生する確率を近似した予め定められた関数で用いるパラメータを算出し、算出された前記パラメータを用いた前記関数を算出するようにすることができる。これによって、取得可能な統計量を用いて、事故が発生する確率を算出することができる。   Further, the statistical information acquisition means includes the driving information representing a probability that each value of a feature amount related to driving behavior of a vehicle is obtained, and the accident statistical information representing a probability that an accident caused by the feature amount occurs. And the calculation means is the constraint on the driving information and the accident statistical information acquired by the statistical information acquisition means, and the probability that each value of the feature value is obtained under the condition where the accident occurred. A parameter used in a predetermined function approximating the probability of occurrence of the accident with respect to each value of the feature amount based on a condition, and calculating the function using the calculated parameter Can be. Thereby, the probability that an accident will occur can be calculated using the statistics that can be acquired.

また、前記統計情報取得手段は、車両の運転行動に関する複数次元の特徴量の各次元についての各値が得られる確率を表す前記運転情報と、前記複数次元の特徴量に起因して生じる事故が発生する確率を表す前記事故統計情報とを取得し、前記算出手段は、前記統計情報取得手段によって取得された前記運転情報及び前記事故統計情報と、前記複数次元の特徴量の各次元についての前記事故が発生した条件下での前記次元の各値が得られる確率に関する前記制約条件とに基づいて、前記特徴量の各次元の各値に対して前記事故が発生する確率を近似した予め定められた関数で用いるパラメータを算出し、算出された前記パラメータを用いた前記関数を算出するようにすることができる。これによって、取得可能な統計量を用いて、事故が発生する確率を算出することができる。   Further, the statistical information acquisition means includes the driving information indicating the probability that each value of the multi-dimensional feature value related to the driving behavior of the vehicle is obtained, and an accident caused by the multi-dimensional feature value. The accident statistical information representing the probability of occurrence is acquired, and the calculation means is the driving information and accident statistical information acquired by the statistical information acquisition means, and the dimensions of each of the dimensions of the plurality of dimensions. Based on the constraint condition regarding the probability that each value of the dimension is obtained under the condition where the accident occurred, a predetermined value approximating the probability that the accident will occur for each value of each dimension of the feature quantity is determined in advance. The parameter used in the function can be calculated, and the function using the calculated parameter can be calculated. Thereby, the probability that an accident will occur can be calculated using the statistics that can be acquired.

また、前記統計情報取得手段は、車両の運転行動に関する複数の特徴量についての各値が得られる確率を表す前記運転情報と、前記複数の特徴量の各々に起因して生じる複数種類の事故の各々が発生する確率を表す前記事故統計情報とを取得し、前記算出手段は、前記複数種類の事故の各々について、前記統計情報取得手段によって取得された前記運転情報及び前記事故統計情報と、前記事故が発生した条件下での前記事故に対応する前記特徴量の各値が得られる確率に関する前記制約条件とに基づいて、前記特徴量の各値に対して前記事故が発生する確率を近似した予め定められた関数で用いるパラメータを算出し、算出された前記パラメータを用いた前記関数を算出するようにすることができる。これによって、取得可能な統計量を用いて、事故が発生する確率を算出することができる。   Further, the statistical information acquisition means includes the driving information indicating the probability of obtaining each value for a plurality of feature values related to driving behavior of the vehicle, and a plurality of types of accidents caused by each of the plurality of feature values. The accident statistical information representing the probability of each occurrence, the calculation means, for each of the plurality of types of accidents, the driving information and accident statistical information acquired by the statistical information acquisition means, Based on the constraint condition regarding the probability that each value of the feature value corresponding to the accident under the condition where the accident occurred is obtained, the probability that the accident will occur is approximated for each value of the feature value A parameter used in a predetermined function can be calculated, and the function using the calculated parameter can be calculated. Thereby, the probability that an accident will occur can be calculated using the statistics that can be acquired.

また、本発明は、対象ドライバの車両の運転行動に関する前記複数の特徴量の各々の値を取得する特徴量取得手段と、前記特徴量取得手段によって取得された前記複数の特徴量の各々の値、及び前記算出手段によって前記複数種類の事故の各々について算出された前記関数に基づいて、前記複数種類の事故の各々について前記事故が発生しない確率を算出し、前記対象ドライバの車両の運転行動について、前記複数種類の事故の各々に対して算出される前記事故が発生しない確率の積を、事故が起こらない確率を表す無事故確率として評価する運転行動評価手段とを更に含むようにすることができる。これによって、取得可能な統計量を用いて、ドライバの運転行動を評価することができる。   Further, the present invention provides a feature amount acquisition unit that acquires each value of the plurality of feature amounts related to the driving behavior of the vehicle of the target driver, and each value of the plurality of feature amounts acquired by the feature amount acquisition unit. And a probability that the accident does not occur for each of the plurality of types of accidents based on the function calculated for each of the plurality of types of accidents by the calculation means, and the driving behavior of the vehicle of the target driver And a driving behavior evaluation means for evaluating a product of the probability that the accident does not occur calculated for each of the plurality of types of accidents as an accident-free probability representing a probability that the accident does not occur. . Thus, the driving behavior of the driver can be evaluated using the statistics that can be acquired.

また、前記特徴量取得手段は、対象ドライバの車両の運転行動に関する前記特徴量の値、及び前記特徴量の値の継続時間を取得し、前記運転行動評価手段は、前記特徴量取得手段によって取得された前記特徴量の値及び前記継続時間と、前記算出手段によって前記複数種類の事故の各々について算出された、予め定められた継続時間までに前記特徴量の各値に対して前記事故が発生する確率を近似した前記関数とに基づいて、前記対象ドライバの車両の運転行動について、前記複数種類の事故の各々に対して算出される前記事故が発生しない確率の積を、前記無事故確率として評価するようにすることができる。これによって、時々刻々と変化する特徴量に対して無事故確率の算出をすることができる。   In addition, the feature amount acquisition unit acquires the value of the feature amount related to the driving behavior of the vehicle of the target driver and the duration of the value of the feature amount, and the driving behavior evaluation unit acquires the feature amount acquisition unit. The accident occurs for each value of the feature amount by a predetermined duration calculated for each of the plurality of types of accidents by the calculating means and the duration value of the feature amount. On the basis of the function approximating the probability of performing, the product of the probability that the accident does not occur calculated for each of the plurality of types of accidents is evaluated as the accident-free probability for the driving behavior of the target driver's vehicle To be able to. As a result, it is possible to calculate the accident-free probability for the feature quantity that changes from moment to moment.

また、前記特徴量取得手段は、対象ドライバの車両の運転行動に関する前記特徴量の値、及び前記特徴量の値が得られているときの走行距離を取得し、前記運転行動評価手段は、前記特徴量取得手段によって取得された前記特徴量の値及び前記走行距離と、前記算出手段によって前記複数種類の事故の各々について算出された、予め定められた走行距離までに前記特徴量の各値に対して前記事故が発生する確率を近似した前記関数とに基づいて、前記対象ドライバの車両の運転行動について、前記複数種類の事故の各々に対して算出される前記事故が発生しない確率の積を、前記無事故確率として評価するようにすることができる。これによって、時々刻々と変化する特徴量に対して無事故確率の算出をすることができる。   Further, the feature amount acquisition means acquires the value of the feature amount related to the driving behavior of the vehicle of the target driver, and the travel distance when the value of the feature amount is obtained, and the driving behavior evaluation means The value of the feature amount and the travel distance acquired by the feature amount acquisition means, and each value of the feature amount by the predetermined travel distance calculated for each of the plurality of types of accidents by the calculation means. On the basis of the function that approximates the probability that the accident will occur, the product of the probability that the accident will not occur calculated for each of the plurality of types of accidents for the driving behavior of the target driver's vehicle. , And can be evaluated as the accident-free probability. As a result, it is possible to calculate the accident-free probability for the feature quantity that changes from moment to moment.

また、前記関数は、前記関数の値が最大となるときの前記特徴量の予め定められた値、又は前記関数の値が最小となるときの前記特徴量の予め定められた値を含むようにすることができる。   Further, the function includes a predetermined value of the feature amount when the value of the function is maximum, or a predetermined value of the feature amount when the value of the function is minimum. can do.

また、前記算出手段は、前記統計情報取得手段によって取得された前記運転情報及び前記事故統計情報と、前記複数次元の特徴量の各次元についての前記事故が発生した条件下での前記次元の各値が得られる確率に関する前記制約条件とに基づいて、ベイズの定理を用い、かつ前記複数次元の特徴量の各次元が互いに独立又は従属であると仮定して、前記パラメータを算出するようにすることができる。   In addition, the calculation unit may include the driving information and the accident statistical information acquired by the statistical information acquisition unit, and each dimension of the dimensions under the condition where the accident has occurred for each dimension of the multi-dimensional feature amount. Based on the constraint on the probability of obtaining a value, the parameter is calculated using Bayes' theorem and assuming that the dimensions of the feature quantities of the plurality of dimensions are independent or dependent on each other. be able to.

また、前記関数は、前記特徴量の各値に対して前記事故が発生する度合いを近似したガウス関数とすることができる。   The function may be a Gaussian function that approximates the degree of occurrence of the accident for each value of the feature amount.

また、前記関数は、前記特徴量の各値に対して前記事故が発生する度合いを近似した線形関数とすることができる。   The function may be a linear function that approximates the degree of occurrence of the accident with respect to each value of the feature amount.

また、前記関数は、前記特徴量の各値に対して前記事故が発生する度合いを近似したシグモイド関数とすることができる。   The function may be a sigmoid function that approximates the degree of occurrence of the accident for each value of the feature amount.

以上説明したように、本発明の事故情報算出装置、及びプログラムによれば、運転情報及び事故統計情報と、事故が発生した条件下での特徴量の各値が得られる度合いに関する制約条件とに基づいて、特徴量の各値に対して事故が発生する度合いを算出することにより、取得可能な統計量を用いて、事故が発生する度合いを算出することができる、という効果が得られる。   As described above, according to the accident information calculation apparatus and program of the present invention, the driving information and the accident statistical information, and the constraint condition relating to the degree to which each value of the feature value under the condition where the accident occurred can be obtained. On the basis of this, by calculating the degree of occurrence of an accident for each value of the feature amount, an effect that the degree of occurrence of the accident can be calculated by using an acquirable statistic.

無事故確率を算出する概念を示す図である。It is a figure which shows the concept which calculates an accident-free probability. 本発明の実施の形態の事故情報算出装置の構成を示すブロック図である。It is a block diagram which shows the structure of the accident information calculation apparatus of embodiment of this invention. ベイズの定理を適用させて統計量を算出する概念を示す図である。It is a figure which shows the concept which calculates a statistic applying Bayes' theorem. 本発明の実施の形態の事故情報算出装置におけるパラメータ算出処理ルーチンの内容を示すフローチャートである。It is a flowchart which shows the content of the parameter calculation process routine in the accident information calculation apparatus of embodiment of this invention. 本発明の実施の形態の事故情報算出装置における運転診断処理ルーチンの内容を示すフローチャートである。It is a flowchart which shows the content of the driving | operation diagnosis processing routine in the accident information calculation apparatus of embodiment of this invention.

<概要>
ドライバの運転能力を維持・向上させるためには、継続的且つ的確な措置が非常に重要となる。しかし、現状では、ドライバの運転行動を診断するための評価指標が確立していない。本発明の実施の形態では、ドライバの運転行動を診断する評価指標の一つとして、「無事故確率(事故を起こさない確率)」を定義し、無事故確率の算出方法を説明する。
<Overview>
In order to maintain and improve the driving ability of the driver, continuous and appropriate measures are very important. However, at present, an evaluation index for diagnosing the driver's driving behavior has not been established. In the embodiment of the present invention, “accident-free probability (probability of not causing an accident)” is defined as one of evaluation indexes for diagnosing the driving behavior of the driver, and a method for calculating the accident-free probability will be described.

図1に本発明の実施の形態の概念図を示す。図1に示すように、本発明の実施の形態では、取得された特徴量の値に対して事故が発生する確率に基づいて、事故が起こらない確率を表す無事故確率を算出する。   FIG. 1 shows a conceptual diagram of an embodiment of the present invention. As shown in FIG. 1, in the embodiment of the present invention, an accident-free probability representing the probability that an accident will not occur is calculated based on the probability that an accident will occur with respect to the acquired feature value.

以下、図面を参照して本発明の実施の形態を詳細に説明する。本発明の実施の形態では、車両に搭載され、取得された特徴量の値に対して事故が発生する確率を算出し、算出された事故が発生する確率に基づいて、事故が起こらない確率を表す無事故確率を算出する事故情報算出装置に本発明を適用した場合を例に説明する。   Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings. In the embodiment of the present invention, the probability that an accident will occur is calculated for the value of the feature quantity that is mounted on the vehicle and acquired, and the probability that no accident will occur is calculated based on the calculated probability that the accident will occur. A case where the present invention is applied to an accident information calculation apparatus that calculates an accident-free probability is described as an example.

[第1の実施の形態]
<事故情報算出装置の構成>
本発明の第1の実施の形態に係る事故情報算出装置について説明する。図2に示すように、本発明の実施の形態に係る事故情報算出装置10は、車両の横方向加速度を検出する横方向加速度センサ12と、車両の前後方向加速度を検出する前後方向加速度センサ14と、車両の速度を検出する車速センサ16と、自車両と先行車両との間の車間距離を検出する車間距離センサ18と、道路領域を含む自車両前方を逐次撮像して前方画像を出力する前方カメラ20と、車室内のドライバの目領域を含む顔画像を逐次撮像して出力するドライバカメラ22と、事故統計情報と運転情報とが格納された統計情報データベース24と、無事故確率を算出するコンピュータ26と、コンピュータ26によって算出された無事故確率を出力する出力装置40とを備えている。
[First Embodiment]
<Configuration of accident information calculation device>
An accident information calculation apparatus according to the first embodiment of the present invention will be described. As shown in FIG. 2, the accident information calculation apparatus 10 according to the embodiment of the present invention includes a lateral acceleration sensor 12 that detects lateral acceleration of the vehicle, and a longitudinal acceleration sensor 14 that detects longitudinal acceleration of the vehicle. A vehicle speed sensor 16 that detects the speed of the vehicle, an inter-vehicle distance sensor 18 that detects an inter-vehicle distance between the host vehicle and a preceding vehicle, and sequentially images the front of the host vehicle including the road area and outputs a front image. The front camera 20, a driver camera 22 that sequentially captures and outputs a face image including the driver's eye area in the vehicle interior, a statistical information database 24 that stores accident statistical information and driving information, and an accident-free probability are calculated. A computer 26 and an output device 40 that outputs the accident-free probability calculated by the computer 26 are provided.

車間距離センサ18は、車間距離として自車両と先行車両との相対距離を検出してコンピュータ26に出力する。車間距離センサ18としては、例えば、レーザレーダやミリ波レーダ等を用いることができる。   The inter-vehicle distance sensor 18 detects the relative distance between the host vehicle and the preceding vehicle as the inter-vehicle distance and outputs it to the computer 26. As the inter-vehicle distance sensor 18, for example, a laser radar, a millimeter wave radar, or the like can be used.

前方カメラ20は、車両に搭載され、道路領域を含む自車両前方を撮像して画像を出力する。前方カメラ20は、自車両前方を撮像し、画像の画像信号を生成する撮像部(図示省略)と、撮像部で生成された画像信号をA/D変換するA/D変換部(図示省略)と、A/D変換された画像信号を一時的に格納するための画像メモリ(図示省略)とを備えている。   The front camera 20 is mounted on the vehicle, images the front of the host vehicle including the road area, and outputs an image. The front camera 20 images the front of the host vehicle and generates an image signal of the image (not shown), and an A / D converter (not shown) that A / D converts the image signal generated by the imaging unit. And an image memory (not shown) for temporarily storing the A / D converted image signal.

ドライバカメラ22は、ドライバの顔画像を撮像して出力する。ドライバカメラ22は、ドライバの顔を撮像し、画像の画像信号を生成する撮像部(図示省略)と、撮像部で生成された画像信号をA/D変換するA/D変換部(図示省略)と、A/D変換された画像信号を一時的に格納するための画像メモリ(図示省略)とを備えている。   The driver camera 22 captures and outputs a driver's face image. The driver camera 22 images a driver's face and generates an image signal of the image (not shown), and an A / D conversion unit (not shown) that A / D converts the image signal generated by the imaging unit. And an image memory (not shown) for temporarily storing the A / D converted image signal.

統計情報データベース24には、車両の運転行動に関する複数の特徴量Aについての各値が得られる確率p(A)を表す運転情報と、複数の特徴量Aの各々に起因して生じる複数種類の事故Xの各々が発生する確率p(X)を表す事故統計情報とが格納されている。 The statistical information database 24 is generated due to driving information representing the probability p (A i ) that each value of the plurality of feature values A i relating to the driving behavior of the vehicle is obtained, and each of the plurality of feature values A i. Accident statistical information indicating the probability p (X i ) of occurrence of each of a plurality of types of accidents X i is stored.

ここで、Aはドライバの運転行動に関する特徴量(速度、加速度、車間距離、脇見時間、等)を表す。また、Xは、特徴量Aの値に起因して生じる事故(スリップ事故、追突事故、脇見事故、等)を表す。また、p(X)は、事故Xが発生する確率を表し、p(A)は、特徴量Aの値が得られる確率を表す。 Here, A i represents a feature amount (speed, acceleration, inter-vehicle distance, side-viewing time, etc.) related to the driving behavior of the driver. X i represents an accident (slip accident, rear-end collision, side-view accident, etc.) caused by the value of the feature value A i . Further, p (X i ) represents the probability that the accident X i will occur, and p (A i ) represents the probability that the value of the feature amount A i will be obtained.

事故統計情報は、「l[km]走行当たりの事故率」「t[sec]走行当たりの事故率」等のように、任意の量当たりの事故率として正規化され、統計情報データベース24に格納される。例えば、事故Xが追突事故の場合、全国の1年間の追突事故の件数N、全国の1年間の平均走行距離L[km]、全国の自動車所有台数Cの情報を用いて、l[km]走行当たりの追突事故が起こる確率p(X)は、以下の式(1)のように予め算出され、統計情報データベース24に格納される。確率p(X)は、複数種類の事故の各々について算出される。 The accident statistical information is normalized as an accident rate per arbitrary amount, such as “accident rate per l [km] travel”, “accident rate per t [sec] travel”, and stored in the statistical information database 24. Is done. For example, in the case of an accident X i is rear-end collisions, with the number N i of 1 year rear-end collisions national average mileage of 1 year nationwide L [miles], the information of car ownership number C nationwide, l [ km] The probability p (X i ) of occurrence of a rear-end collision per driving is calculated in advance as in the following equation (1) and stored in the statistical information database 24. The probability p (X i ) is calculated for each of a plurality of types of accidents.

また、運転情報の確率p(A)は、事故統計情報が算出されたときと同じ母集団で計算されることが好ましい。例えば、事故統計情報を全国の事故件数から求めた場合、運転情報の確率p(A)は全国の車両の走行情報を統計した際の分布でなければならない。現状では、運転情報は取得困難であるため、実際には少数のドライバの運転データから計算せざるを得ない事もあるが、多数のドライバの走行データを蓄積していき、蓄積された走行データを用いて運転情報の確率p(A)を求めることで、運転情報の確率p(A)は全国の平均分布に漸近していくと考えられる。 In addition, the driving information probability p (A i ) is preferably calculated in the same population as when the accident statistical information is calculated. For example, when the accident statistical information is obtained from the number of accidents in the whole country, the probability p (A i ) of the driving information must be a distribution obtained when the traveling information of the vehicles in the whole country is statistics. Currently, it is difficult to obtain driving information, so it may actually be necessary to calculate from the driving data of a small number of drivers, but the driving data of a large number of drivers will be accumulated, by obtaining the probability p (a i) of the operation information using the probability p (a i) of the operation information is considered to gradually approaches the average distribution across the country.

コンピュータ26は、事故情報算出装置10全体の制御を司るCPU、後述する運転診断処理ルーチンのプログラム等を記憶した記憶媒体としてのROM、ワークエリアとしてデータを一時格納するRAM、及びこれらを接続するバスを含んで構成されている。このような構成の場合には、各構成要素の機能を実現するためのプログラムをROMやHDD等の記憶媒体に記憶しておき、これをCPUが実行することによって、各機能が実現されるようにする。   The computer 26 includes a CPU that controls the entire accident information calculation apparatus 10, a ROM as a storage medium that stores a program of a driving diagnosis processing routine that will be described later, a RAM that temporarily stores data as a work area, and a bus that connects these It is comprised including. In the case of such a configuration, a program for realizing the function of each component is stored in a storage medium such as a ROM or HDD, and each function is realized by executing the program by the CPU. To.

このコンピュータ26をハードウエアとソフトウエアとに基づいて定まる機能実現手段毎に分割した機能ブロックで説明すると、図2に示すように、横方向加速度センサ12によって検出された車両の横方向加速度、前後方向加速度センサ14によって検出された車両の前後方向加速度、車速センサ16によって検出された車両の速度、車間距離センサ18によって検出された自車両と先行車両との間の車間距離、前方カメラ20によって撮像された前方画像、及びドライバカメラ22によって撮像されたドライバの顔画像に基づいて、複数の特徴量を取得する特徴量取得部28と、統計情報データベース24に格納された事故統計情報と運転情報とを取得する統計情報取得部30と、特徴量取得部28によって取得された複数の特徴量と、統計情報取得部30によって取得された事故統計情報及び運転情報とに基づいて、事故が発生する確率を近似した予め定められた関数で用いるパラメータを算出する算出部32と、算出部32によって算出されたパラメータを用いた関数に基づいて、対象ドライバの運転を評価する運転行動評価部34とを備えている。   When the computer 26 is described with function blocks divided for each function realizing means determined based on hardware and software, as shown in FIG. 2, the lateral acceleration of the vehicle detected by the lateral acceleration sensor 12, the longitudinal The longitudinal acceleration of the vehicle detected by the directional acceleration sensor 14, the vehicle speed detected by the vehicle speed sensor 16, the inter-vehicle distance between the host vehicle and the preceding vehicle detected by the inter-vehicle distance sensor 18, and imaged by the front camera 20. A feature amount acquisition unit 28 that acquires a plurality of feature amounts based on the forward image that has been captured and the driver's face image captured by the driver camera 22; accident statistical information and driving information stored in the statistical information database 24; A plurality of feature amounts acquired by the statistical information acquisition unit 30 and the feature amount acquisition unit 28; Based on accident statistics information and driving information acquired by the meter information acquisition unit 30, a calculation unit 32 that calculates parameters used in a predetermined function that approximates the probability of occurrence of an accident, and a calculation unit 32 And a driving behavior evaluation unit 34 that evaluates the driving of the target driver based on a function using the parameters.

特徴量取得部28は、運転診断の対象となるドライバの車両の運転行動に関する複数の特徴量の各々の値を取得する。   The feature amount acquisition unit 28 acquires each value of a plurality of feature amounts related to the driving behavior of the vehicle of the driver that is the target of the driving diagnosis.

例えば、特徴量取得部28は、前方カメラ20によって撮像された前方画像から、画像処理によって一時停止標識を検出し、一時停止標識が検出されたときに車速センサ16によって検出された車両の速度を、一時停止時の通過速度として取得する。   For example, the feature amount acquisition unit 28 detects a stop sign by image processing from the front image captured by the front camera 20, and detects the vehicle speed detected by the vehicle speed sensor 16 when the stop sign is detected. , Get as the passing speed at the time of pause.

また、特徴量取得部28は、ドライバカメラ22によって撮像されたドライバの顔画像に基づいて、例えばパターンマッチング処理によりドライバの顔を検出する。次に、特徴量取得部28は、検出されたドライバの顔に基づいて、複数の異なる顔向きの顔画像に基づいて学習されたパターンマッチング処理により、ドライバの顔向きを検出する。そして特徴量取得部28は、ドライバの顔向きに基づいて、脇見運転が継続している時間を算出し、脇見運転が継続している時間と、脇見運転が継続しているときに車速センサ16によって検出された車両の速度とを取得する。そして、特徴量取得部28は、取得した脇見運転が継続している時間と、車速とに基づいて、脇見運転が継続した脇見運転距離を取得する。   In addition, the feature amount acquisition unit 28 detects the driver's face by, for example, pattern matching processing based on the driver's face image captured by the driver camera 22. Next, the feature amount acquisition unit 28 detects the face direction of the driver by pattern matching processing learned based on a plurality of face images having different face directions based on the detected face of the driver. Then, the feature amount acquisition unit 28 calculates the time during which the aside operation is continued based on the driver's face direction, and the vehicle speed sensor 16 when the aside operation is continued. And the vehicle speed detected by. And the feature-value acquisition part 28 acquires the armpit driving distance which the armpit driving continued based on the acquired time and the vehicle speed.

また、特徴量取得部28は、前方カメラ20によって撮像された前方画像から、画像処理によって信号標識の赤信号を検出し、赤信号が検出されたときの交差点の車両の通過時間を取得する。   The feature amount acquisition unit 28 detects a red signal of a signal sign from the front image captured by the front camera 20 by image processing, and acquires the passing time of the vehicle at the intersection when the red signal is detected.

そして、特徴量取得部28は、横方向加速度センサ12によって検出された横方向加速度と、前後方向加速度センサ14によって検出された前後方向加速度と、車間距離センサ18によって検出された車間距離と、一時停止時の通過速度と、脇見運転距離と、赤信号時の交差点通過時間とを、ドライバの運転行動に関する複数の特徴量として取得する。   The feature amount acquisition unit 28 then detects the lateral acceleration detected by the lateral acceleration sensor 12, the longitudinal acceleration detected by the longitudinal acceleration sensor 14, the inter-vehicle distance detected by the inter-vehicle distance sensor 18, The passing speed at the time of stoppage, the aside driving distance, and the intersection passing time at the time of a red signal are acquired as a plurality of feature amounts related to the driving behavior of the driver.

統計情報取得部30は、統計情報データベース24に格納された事故統計情報と運転情報とを取得する。   The statistical information acquisition unit 30 acquires accident statistical information and driving information stored in the statistical information database 24.

以下では、まず、本実施の形態で算出される無事故確率の定義と、無事故確率を推定するための定式化について説明する。   In the following, first, the definition of the accident-free probability calculated in the present embodiment and the formulation for estimating the accident-free probability will be described.

車両のドライバの運転行動(運転走行データ)から無事故確率(事故が起こらない確率)を推定するためには、事故に起因する特徴量の取得及び取得された特徴量を用いた定式化が課題となるが、無事故確率の真値が取得出来ない限り、一般に回帰分析等による定式化は困難である。そこで、本発明の実施の形態では、実際の事故件数に基づく事故統計情報や、運転データの分布を表す運転情報等の統計情報に着目する。本発明の実施の形態では、以下の近似(A)(B)を用いる。   In order to estimate the accident-free probability (probability that no accident will occur) from the driving behavior (driving data) of the driver of the vehicle, the acquisition of the feature value resulting from the accident and the formulation using the acquired feature value are issues. However, unless a true value of accident-free probability can be obtained, it is generally difficult to formulate by regression analysis or the like. Therefore, in the embodiment of the present invention, attention is focused on accident statistical information based on the actual number of accidents and statistical information such as driving information representing the distribution of driving data. In the embodiment of the present invention, the following approximations (A) and (B) are used.

(A)全ての事故Xが生じる確率は互いに独立である。
(B)特徴量Aのみに起因して事故Xが生じる。
(A) The probabilities that all accidents X i occur are independent of each other.
(B) The accident X i occurs only due to the feature value A i .

上記の近似(A)に関しては、一台の車両(一人のドライバの運転行動)に着目している以上、複数種類の事故が同時に発生する事は殆ど起こりえないため、妥当と考えられる。また、上記の近似(B)に関しては、様々な要因が重なって事故が生じる場合もあるが、例えば「脇見による事故」は「脇見」が主原因というように、多くの事故の要因は単一もしくは少数の特徴量が支配的になると考えられる。そこで、第1の実施の形態では、Aを一次元特徴量として定義した上で無事故確率の定式化を行う。上記の近似(A)(B)を用いることで、無事故確率(事故を起こらない確率)psafeは以下の式(2)で与えられる。 Regarding the above approximation (A), since attention is paid to one vehicle (driving behavior of one driver), it is considered that a plurality of types of accidents hardly occur at the same time, so that it is considered appropriate. In addition, regarding the above approximation (B), an accident may occur due to the overlap of various factors. For example, the “accident due to wakimi” is mainly caused by “wakimi”. Alternatively, it is considered that a small number of features become dominant. Therefore, in the first embodiment, the accident-free probability is formulated after defining A i as a one-dimensional feature value. By using the above approximations (A) and (B), the accident-free probability (probability of not causing an accident) p safe is given by the following equation (2).

なお、p(X|A)は、特徴量Aの値に対して事故Xが発生する確率を表す。 Note that p (X i | A i ) represents the probability that an accident X i will occur with respect to the value of the feature value A i .

次に、特徴量Aの値に対して事故Xが発生する確率p(X|A)を定義する。本発明の実施の形態では、特徴量Aの各値に対して事故Xが発生する確率p(X|A)の各々を、関数を用いて近似する。p(X|A)を、以下の式(3)又は式(4)に示すように、ガウス関数を用いて近似し、定義する。 Next, the probability p (X i | A i ) that an accident X i will occur with respect to the value of the feature value A i is defined. In the embodiment of the present invention, each probability p (X i | A i ) of occurrence of an accident X i is approximated using a function for each value of the feature value A i . p (X i | A i ) is approximated and defined using a Gaussian function as shown in the following formula (3) or formula (4).

上記式(3)におけるAmaxは、1(100%)の確率で事故Xが生じる特徴量Aの値を表し、予め定められた値である。また、上記式(4)におけるAminは、1(100%)の確率で事故Xが生じない特徴量Aの値を表し、予め定められた値である。また、σはガウス関数のパラメータを表す。 A max in the above equation (3) represents a value of the feature amount A i at which the accident X i occurs with a probability of 1 (100%), and is a predetermined value. Further, A min in the above formula (4) represents a value of the feature value A i at which the accident X i does not occur with a probability of 1 (100%), and is a predetermined value. Σ i represents a parameter of a Gaussian function.

算出部32は、複数種類の事故の各々について、統計情報取得部30によって取得された運転情報及び事故統計情報と、事故Xが発生した条件下での事故Xに対応する特徴量Aの各値が得られる確率に関する制約条件とに基づいて、ベイズの定理を用いて、特徴量Aの各値に対して事故Xが発生する確率を近似したガウス関数で用いるパラメータを算出し、算出されたパラメータを用いたガウス関数を、特徴量の各値に対して事故が発生する確率p(X|A)として算出する。 Calculation unit 32, for each of a plurality of types of accidents, and operation information and accident statistics that are acquired by the statistical information acquiring unit 30, a feature quantity corresponding to accident X i under conditions accident X i occurs A i The parameters used in the Gaussian function that approximates the probability of occurrence of the accident X i for each value of the feature value A i are calculated using the Bayes' theorem based on the constraint on the probability of obtaining each value of Then, a Gaussian function using the calculated parameters is calculated as the probability p (X i | A i ) that an accident will occur for each value of the feature quantity.

具体的には、算出部32は、特徴量Aと事故Xとに関する知見を用いて、上記式(3)及び式(4)に含まれるパラメータσ(未知数)を推定する。例えば、特徴量Aと事故Xとに関する知見から、近似したガウス関数の値が最大となる特徴量、及び近似したガウス関数の値が最小となる特徴量の少なくとも一方を予め設定する。 Specifically, the calculation unit 32 estimates the parameter σ i (unknown number) included in the above formulas (3) and (4) by using knowledge about the feature amount A i and the accident X i . For example, the knowledge and the feature quantity A i and accidents X i, feature quantity value of the Gaussian function approximating becomes maximum, and the value of the Gaussian function approximating to preset at least one of smallest feature amount.

例えば、「A:加速度、X:スリップ事故」の場合、加速度が静止摩擦係数と重力加速度との積を超えた際にはスリップ事故が生じるものとし、「Amax=静止摩擦係数×重力加速度」と近似する。(あるいは、「加速していないときの事故確率は0%」を前提に「Amin=0[m/s]」としてもよい。)また、「A:車間距離(又はTTC:Time to Collision)、X:車間距離不足による追突事故」の場合、車間距離が0mのときには必ず追突事故が生じるため、「Amax=0[m]」とする事が出来る。以上の近似により、p(X|A)に残された未知数はパラメータσのみとなる。表1に、本発明の実施の形態で用いる特徴量A、事故X、Amax、Aminの関係を示す。 For example, in the case of “A i : acceleration, X i : slip accident”, it is assumed that a slip accident occurs when the acceleration exceeds the product of the static friction coefficient and the gravitational acceleration, and “A max = static friction coefficient × gravity”. Approximate acceleration. (Alternatively, “A min = 0 [m / s 2 ]” may be set on the assumption that the accident probability when not accelerating is 0%.) Also, “A i : inter-vehicle distance (or TTC: Time to Collision), X i : rear-end collision due to lack of inter-vehicle distance ”, a rear-end collision always occurs when the inter-vehicle distance is 0 m. Therefore,“ A max = 0 [m] ”can be obtained. As a result of the above approximation, the only unknown remaining in p (X i | A i ) is the parameter σ i . Table 1 shows the relationship between the feature value A i , the accident X i , A max , and A min used in the embodiment of the present invention.

上記式(3)又は式(4)のパラメータσは、特徴量Aの値がどれだけ安全なのかを決定するパラメータであり、パラメータの値を上手く設定することで、実際の安全性と矛盾の少ない定式化が可能になると考えられる。従って、近似したガウス関数は、ガウス関数の値が最大となるときの特徴量Aの予め定められた値、又はガウス関数の値が最小となるときの特徴量Aの予め定められた値を含む。 The parameter σ i in the above formula (3) or (4) is a parameter that determines how safe the value of the feature value A i is. By setting the parameter value well, It is thought that formulation with few contradictions is possible. Therefore, the approximated Gaussian function is a predetermined value of the feature value A i when the value of the Gauss function is maximum, or a predetermined value of the feature value A i when the value of the Gauss function is minimum. including.

算出部32は、前述した通り、実際の統計情報である事故統計情報及び運転情報を用いて、ベイズの定理に従って、近似したガウス関数のパラメータσを算出する。そして、算出部32は、算出されたパラメータσを用いたガウス関数を算出する。 As described above, the calculation unit 32 calculates the parameter σ i of the approximated Gaussian function according to the Bayes' theorem using accident statistical information and driving information that are actual statistical information. Then, the calculation unit 32 calculates a Gaussian function using the calculated parameter σ i .

まず、ベイズの定理を適用することで、p(X|A)を以下の式(5)のように変形する。 First, p (X i | A i ) is transformed into the following equation (5) by applying Bayes' theorem.

ここで、p(A |X)は、事故Xが発生した条件下での特徴量Aの値が得られる確率を表し、未知の値である。また、確率p(A)に関しては、以下の式(6)に示すように、p(A)の両側無限積分の値は1となる。 Here, p (A i | X i ) represents the probability of obtaining the value of the feature value A i under the condition where the accident X i has occurred, and is an unknown value. Regarding the probability p (A i ), the value of the two-sided infinite integral of p (A i ) is 1 as shown in the following formula (6).

また、事故Xが発生した条件下での特徴量Aの値が得られる確率p(A|X)に関する制約条件を以下の式(7)に示す。確率p(A|X)に関しては、その分布形状の取得は困難であるが、事故Xは特徴量Aのみに起因して生じるため、以下の式(9)に示すように、p(A|X)の両側無限積分の値は1となる。 Further, a constraint condition regarding the probability p (A i | X i ) at which the value of the feature value A i under the condition where the accident X i occurs is shown in the following formula (7). Regarding the probability p (A i | X i ), it is difficult to obtain the distribution shape, but since the accident X i is caused only by the feature amount A i , as shown in the following equation (9), The value of the two-sided infinite integral of p (A i | X i ) is 1.

算出部32は、上記の式(1)、式(3)、式(4)、式(5)、式(6)、及び式(7)を連立することで、図3に示すように、パラメータσを算出することが可能となる。パラメータσの算出を全ての確率p(A|X)に対して行えば、上記式(2)の無事故確率psafeの算出式が適用可能となる。 As shown in FIG. 3, the calculation unit 32 combines the above equations (1), (3), (4), (5), (6), and (7). The parameter σ i can be calculated. If the calculation of the parameter σ i is performed for all probabilities p (A i | X i ), the calculation formula for the accident-free probability p safe in the above formula (2) can be applied.

以下では、p(X|A)の未知のパラメータであるσの導出を行う際の過程を、上記式(3)、式(5)、及び式(7)を例に説明する。 In the following, the process for deriving σ i , which is an unknown parameter of p (X i | A i ), will be described using the above formulas (3), (5), and (7) as an example.

特徴量Aの値が連続値として定義される場合、上記式(3)を上記式(5)に代入し変形することで、以下の式(8)に示すように、p(A|X)が与えられる。 When the value of the feature value A i is defined as a continuous value, the above equation (3) is substituted into the above equation (5) and deformed, and as shown in the following equation (8), p (A i | X i ) is given.

また、上記式(8)のp(A|X)を、上記式(7)に代入することで、以下の式(9)が得られる。 Further, the above equation p of (8) | a (A i X i), by substituting the above equation (7), the formula (9) below is obtained.

上記式(9)において、値が未知であるp(A|X)が消えて、未知数はパラメータσのみとなるため、上記式(9)を解くことでパラメータσを得ることができる。ここで、事前に取得しておいたp(A)がAの関数として与えられており、上記式(9)の積分が解析的に可能である場合には、実際に積分してパラメータσを求めれば良い。上記式(9)の積分が困難な場合、またはp(A)がAの関数として定義困難な場合には、A=ai,k(k=1,2,・・・,Nai)としてp(A)を離散化すればよく、上記式(7)及び式(9)はそれぞれ以下の式(10)及び式(11)として近似可能である。 In the formula (9), the value p is unknown | disappears (A i X i), since the unknown is only parameter sigma i, to obtain the parameter sigma i by solving the above expression (9) it can. Here, when p (A i ) acquired in advance is given as a function of A i and the integration of the above equation (9) is analytically possible, it is actually integrated and parameter What is necessary is just to obtain | require (sigma) i . If the integral of the equation (9) is difficult, or p if (A i) are difficult defined as a function of A i is, A i = a i, k (k = 1,2, ···, N It is only necessary to discretize p (A i ) as a i ), and the above equations (7) and (9) can be approximated as the following equations (10) and (11), respectively.

したがって、上記式(11)を解くことで、未知のパラメータσが得られる。特徴量Aの値が離散値A=ai,k(k=1,2,・・・,Nai)として定義されている場合も同様に、式(11)を解くことで、未知のパラメータσが得られる。また、上記では式(3)のp(X|A)を例に未知のパラメータを算出したが、式(4)を用いる場合にも、式(9)もしくは式(11)のp(X|A)、またはp(X|ai,k)に代入することで、同様にして未知のパラメータの算出が可能である。 Therefore, the unknown parameter σ i can be obtained by solving the above equation (11). Similarly, when the value of the feature value A i is defined as a discrete value A i = a i, k (k = 1, 2,..., N ai ), Parameter σ i is obtained. In the above description, unknown parameters are calculated by taking p (X i | A i ) of Equation (3) as an example. However, when Equation (4) is used, p (X) of Equation (9) or Equation (11) is also used. By substituting into X i | A i ) or p (X i | a i, k ), unknown parameters can be calculated in the same manner.

以上のように本発明の実施の形態では、特徴量A、事故X、Amax、Amin、p(X|A)を知見や仮説を基に定義していき、推測困難な未知数に関しては事故や運転データの統計情報を用いて算出することで、無事故確率の定式化を行った。 As described above, in the embodiment of the present invention, the feature amount A i , the accident X i , A max , A min , and p (X i | A i ) are defined based on knowledge and hypothesis, and are difficult to estimate. For the unknowns, the accident-free probability was formulated by using statistical information on accidents and driving data.

運転行動評価部34は、特徴量取得部28によって取得された複数の特徴量の各々の値、及び算出部32によって複数種類の事故の各々について算出された、ガウス関数を用いて近似された確率p(X|A)に基づいて、複数種類の事故の各々について当該事故Xが発生しない確率(1−p(X|A))を算出する。そして、運転行動評価部34は、対象ドライバの車両の運転行動について、複数種類の事故の各々に対して算出される事故が発生しない確率(1−p(X|A))の積を、事故が起こらない確率を表す無事故確率psafeとして算出する。 The driving behavior evaluation unit 34 is a probability approximated using a Gaussian function calculated for each of a plurality of types of accidents by the calculation unit 32 and each value of the plurality of feature amounts acquired by the feature amount acquisition unit 28. Based on p (X i | A i ), the probability (1-p (X i | A i )) that the accident X i does not occur is calculated for each of a plurality of types of accidents. Then, the driving behavior evaluation unit 34 calculates the product of the probability of occurrence of an accident (1-p (X i | A i )) calculated for each of a plurality of types of accidents regarding the driving behavior of the target driver's vehicle. , The accident-free probability psafe representing the probability that no accident will occur is calculated.

出力装置40は、運転行動評価部34によって算出された無事故確率psafeを結果として出力する。 The output device 40 outputs the accident-free probability p safe calculated by the driving behavior evaluation unit 34 as a result.

<事故情報算出装置10の作用>
次に、第1の実施の形態に係る事故情報算出装置10の作用について説明する。まず、コンピュータ26によって、図4に示すパラメータ算出処理ルーチンが実行される。
<Operation of the accident information calculation device 10>
Next, the operation of the accident information calculation apparatus 10 according to the first embodiment will be described. First, the computer 26 executes a parameter calculation processing routine shown in FIG.

ステップS100において、統計情報取得部30によって、統計情報データベース24に格納された事故統計情報を取得する。   In step S <b> 100, the statistical information acquisition unit 30 acquires accident statistical information stored in the statistical information database 24.

ステップS102において、統計情報取得部30によって、統計情報データベース24に格納された運転情報を取得する。   In step S <b> 102, the driving information stored in the statistical information database 24 is acquired by the statistical information acquisition unit 30.

ステップS104において、算出部32によって、複数種類の事故の各々について、上記ステップS102で取得された運転情報及び上記ステップS100で取得された事故統計情報と、事故Xが発生した条件下での事故Xに対応する特徴量Aの各値が得られる確率に関する制約条件とに基づいて、ベイズの定理を用いて、上記式(1)、式(3)、式(4)、式(5)、式(6)、及び式(7)を連立して、特徴量Aの各値に対して事故Xが発生する確率をガウス関数を用いて近似した上記式(3)又は式(4)に示す関数で用いるパラメータを算出し、算出されたパラメータを用いた関数を、特徴量Aの各値に対して事故が発生する確率p(X|A)として算出する。 In step S104, the calculation unit 32, an accident for each of the plurality of types of accidents, and accident statistics that are acquired by the operation information and the step S100 which is obtained in step S102, under conditions of an accident X i occurs Based on the constraint condition relating to the probability that each value of the feature value A i corresponding to X i is obtained, the above formulas (1), (3), (4), (5) ), Equation (6), and Equation (7) are combined to approximate the probability of occurrence of accident X i for each value of feature value A i using a Gaussian function. The parameter used in the function shown in 4) is calculated, and the function using the calculated parameter is calculated as the probability p (X i | A i ) that an accident will occur for each value of the feature value A i .

次に、診断対象のドライバの運転操作により車両が走行し、横方向加速度センサ12によって車両の横方向加速度が検出され、前後方向加速度センサ14によって車両の前後方向加速度が検出され、車速センサ16によって車両の速度が検出され、車間距離センサ18によって自車両と先行車両との間の車間距離が検出され、前方カメラ20によって前方画像が撮像され、ドライバカメラ22によってドライバの顔画像が撮像されているときに、コンピュータ26によって、図5に示す運転診断処理ルーチンが実行される。   Next, the vehicle travels by the driving operation of the driver to be diagnosed, the lateral acceleration of the vehicle is detected by the lateral acceleration sensor 12, the longitudinal acceleration of the vehicle is detected by the longitudinal acceleration sensor 14, and the vehicle speed sensor 16 The vehicle speed is detected, the inter-vehicle distance sensor 18 detects the inter-vehicle distance between the host vehicle and the preceding vehicle, the front camera 20 captures the front image, and the driver camera 22 captures the driver's face image. Sometimes, the computer 26 executes the driving diagnosis processing routine shown in FIG.

ステップS200において、特徴量取得部28によって、横方向加速度センサ12によって検出された横方向加速度と、前後方向加速度センサ14によって検出された前後方向加速度と、一時停止時の通過速度と、車間距離センサ18によって検出された車間距離と、一時停止時の通過速度と、脇見運転距離と、赤信号時の交差点通過時間とを、ドライバの運転行動に関する複数の特徴量として取得する。   In step S200, the lateral acceleration detected by the lateral acceleration sensor 12, the longitudinal acceleration detected by the longitudinal acceleration sensor 14, the passing speed at the time of pause, the inter-vehicle distance sensor, by the feature amount acquisition unit 28. The inter-vehicle distance detected by 18, the passing speed at the time of the temporary stop, the armpit driving distance, and the intersection passing time at the time of the red light are acquired as a plurality of feature amounts related to the driving behavior of the driver.

ステップS202において、運転行動評価部34によって、上記ステップS200で取得された複数の特徴量Aの各々の値、及びパラメータ算出処理ルーチンの上記ステップS104で算出された複数種類の事故の各々について算出されたガウス関数に基づいて、複数種類の事故の各々について当該事故Xが発生しない確率(1−p(X|A))を算出する。 In step S202, the driving action evaluation unit 34, each value of the plurality of feature quantity A i obtained in step S200, and for each of the plurality of types of accidents that have been calculated in step S104 of parameter calculation routine calculates Based on the Gaussian function, a probability (1-p (X i | A i )) that the accident X i does not occur is calculated for each of a plurality of types of accidents.

ステップS204において、運転行動評価部34によって、対象ドライバの車両の運転行動について、上記ステップS202で算出された複数種類の事故の各々に対して算出される事故が発生しない確率(1−p(X|A))の積を、事故が起こらない確率を表す無事故確率psafeとして算出する。 In step S204, the probability (1−p (X) that the accident calculated for each of the plurality of types of accidents calculated in step S202 regarding the driving behavior of the target driver's vehicle by the driving behavior evaluation unit 34 does not occur. The product of i | A i )) is calculated as the accident-free probability p safe representing the probability that no accident will occur.

ステップS112において、出力装置40によって、上記ステップS204で算出された無事故確率psafeを結果として出力する。 In step S112, the output unit 40, and outputs the accident probability p safe calculated in step S204 as a result.

このように、第1の実施の形態の事故情報算出装置によれば、複数種類の事故の各々について、運転情報及び事故統計情報と、事故が発生した条件下での事故に対応する特徴量の各値が得られる確率に関する制約条件とに基づいて、特徴量の各値に対して事故が発生する確率を近似した予め定められた関数で用いるパラメータを算出し、算出されたパラメータを用いた当該関数を算出し、対象ドライバの車両の運転行動に関する複数の特徴量の各々の値を取得し、複数の特徴量の各々の値、及び算出された複数種類の事故の各々について算出された関数に基づいて、複数種類の事故の各々について事故が発生しない確率を算出し、対象ドライバの車両の運転行動について、無事故確率として評価することにより、取得可能な統計量を用いてドライバの運転を診断することができる。   As described above, according to the accident information calculation apparatus of the first embodiment, for each of a plurality of types of accidents, the driving information and the accident statistical information, and the feature amount corresponding to the accident under the condition in which the accident has occurred. Based on the constraint on the probability that each value is obtained, the parameter used in a predetermined function that approximates the probability of occurrence of an accident for each value of the feature value is calculated, and the parameter using the calculated parameter The function is calculated, each value of the plurality of feature values related to the driving behavior of the vehicle of the target driver is obtained, and each of the plurality of feature values and the calculated function for each of the plurality of types of accidents are obtained. Based on this, the probability that no accident will occur for each of the multiple types of accidents is calculated, and the driving behavior of the target driver's vehicle is evaluated as an accident-free probability. It is possible to diagnose the operation of the driver.

また、取得が困難な統計量を用いずに、ドライバの運転を診断することができる。また、ドライバの運転を精度よく診断することができる。   In addition, the driver's driving can be diagnosed without using statistics that are difficult to obtain. In addition, the operation of the driver can be diagnosed with high accuracy.

また、各事故と要因間の大規模な情報を必要とせず容易に観測可能な統計量を用いて、経時変化し得る実際の交通事情に則した客観的な運転の診断及び無事故確率の算出が可能となる。   In addition, by using statistics that can be easily observed without requiring large-scale information between accidents and factors, objective driving diagnosis and accident-free probability calculation can be performed in accordance with actual traffic conditions that can change over time. It becomes possible.

また、加速度、脇見運転時間等の特徴量と、スリップ事故、脇見事故等の事故の種類を関連付けたうえで、条件付確率に関する定理に基づくことで、多数ドライバから収集した運転情報の統計量及び事故の統計量から、各運転状態に対する無事故確率が推定可能となり、これに基づき運転の良し悪しを診断することができる。容易に観測可能な統計量のみで実際の交通事情に則した客観的な運転の診断が期待され、且つ統計量のみを逐次更新することで交通事情の経時変化にも対応が可能となる。   In addition, by associating features such as acceleration and armpit driving time with the types of accidents such as slip accidents and armpit accidents, and based on theorems on conditional probabilities, statistics of driving information collected from many drivers and The accident-free probability for each driving state can be estimated from the accident statistics, and based on this, it is possible to diagnose whether driving is good or bad. Objective driving diagnosis according to actual traffic conditions is expected only with easily observable statistics, and it is possible to cope with changes in traffic conditions over time by sequentially updating only statistics.

[第2の実施の形態]
次に、第2の実施の形態について説明する。なお、第2の実施の形態に係る事故情報算出装置の構成は、第1の実施の形態と同様であるため、同一符号を付して説明を省略する。
[Second Embodiment]
Next, a second embodiment will be described. In addition, since the structure of the accident information calculation apparatus which concerns on 2nd Embodiment is the same as that of 1st Embodiment, it attaches | subjects the same code | symbol and abbreviate | omits description.

第2の実施の形態では、特徴量の各値に対して事故が発生する確率を近似した関数が、複数のパラメータを含む点が、第1の実施の形態と異なっている。   The second embodiment is different from the first embodiment in that a function that approximates the probability that an accident will occur for each value of a feature value includes a plurality of parameters.

上記第1の実施の形態で示した確率p(X|A)を近似する関数は、任意に選択することができる。例えば、以下の式(12)又は式(13)に示すような、線形関数やシグモイド関数等によって、特徴量の各値に対して事故が発生する確率p(X|A)を近似することができる。 The function that approximates the probability p (X i | A i ) shown in the first embodiment can be arbitrarily selected. For example, the probability p (X i | A i ) of occurrence of an accident is approximated for each value of the feature quantity by a linear function, a sigmoid function, or the like as shown in the following formula (12) or formula (13). be able to.

なお、上記式(12)及び式(13)におけるu,vはパラメータを表す。上記第1の実施の形態においては、上記式(1)、式(3)、式(4)、式(5)、式(6)、及び式(7)を連立して算出可能な未知数は1つであるため、未知数が2つ以上になる第2の実施の形態では、p(X|A)を定義した場合には別の拘束条件が必要となる。 Note that u i and v i in the above equations (12) and (13) represent parameters. In the first embodiment, the unknowns that can be calculated by combining the equations (1), (3), (4), (5), (6), and (7) are as follows: Since there is one, in the second embodiment in which the number of unknowns is two or more, another constraint condition is required when p (X i | A i ) is defined.

従って、第1の実施の形態の上記式(3)や上記式(4)の代わりに、上記式(12)や上記式(13)のp(X|A)を採用した場合、未知のパラメータであるu、vのいずれか一方は経験的に定めたり、別の情報に基づいて定めたりしておく必要がある。 Therefore, when p (X i | A i ) of the above formula (12) or the above formula (13) is adopted instead of the above formula (3) or the above formula (4) of the first embodiment, the unknown one of u i, v i is a parameter, it is necessary to or determined based on or determined empirically, other information.

あるいは、上記式(12)や上記式(13)に関して、例えばAmax、Aminのいずれかが既知である場合、以下の式のいずれかが成立する。 Alternatively, for example, when any one of A max and A min is known with respect to the above formula (12) and the above formula (13), one of the following formulas is established.

すなわち、上記式(14)及び式(15)を用いることで、上記式(12)及び式(13)は以下の式(16)及び式(17)に示すように、1つの未知のパラメータで定義可能である。 That is, by using the above formulas (14) and (15), the above formulas (12) and (13) are expressed by one unknown parameter as shown in the following formulas (16) and (17). It can be defined.

ここで、Zは十分大きな正の値、Z−∞は十分小さな負の値として事前に定義しておく。以上のようにして、p(X|A)の複数の未知のパラメータであるu、vの導出が可能となる。 Here, Z is sufficiently large positive value, Z -∞ is predefine as a sufficiently small negative value. As described above, it is possible to derive u i and v i which are a plurality of unknown parameters of p (X i | A i ).

なお、第2の実施の形態に係る事故情報算出装置の他の構成及び作用については、第1の実施の形態と同様であるため、説明を省略する。   In addition, about the other structure and effect | action of the accident information calculation apparatus which concern on 2nd Embodiment, since it is the same as that of 1st Embodiment, description is abbreviate | omitted.

このように、第2の実施の形態の事故情報算出装置によれば、複数種類の事故の各々について、事故統計情報及び運転情報と、事故が発生した条件下での事故に対応する特徴量の各値が得られる確率に関する制約条件とに基づいて、特徴量の各値に対して事故が発生する確率を近似した予め定められた関数で用いる複数のパラメータを算出し、算出された複数のパラメータを用いた当該関数を算出し、対象ドライバの車両の運転行動に関する複数の特徴量の各々の値を取得し、複数の特徴量の各々の値、及び算出された複数種類の事故の各々について算出された関数に基づいて、複数種類の事故の各々について事故が発生しない確率を算出し、対象ドライバの車両の運転行動について、無事故確率として評価することにより、取得可能な統計量を用いてドライバの運転を診断することができる。   As described above, according to the accident information calculation apparatus of the second embodiment, for each of a plurality of types of accidents, the accident statistics information and the driving information, and the feature amount corresponding to the accident under the condition in which the accident has occurred. Based on the constraint on the probability of obtaining each value, calculate a plurality of parameters to be used in a predetermined function that approximates the probability that an accident will occur for each value of the feature value, and calculate the plurality of parameters To calculate each function of the plurality of feature values related to the driving behavior of the vehicle of the target driver, and to calculate each value of the plurality of feature values and each of the calculated types of accidents. Statistics that can be obtained by calculating the probability that an accident will not occur for each of multiple types of accidents based on the calculated function and evaluating the driving behavior of the target driver's vehicle as an accident-free probability It is possible to diagnose operation of the driver used.

[第3の実施の形態]
次に、第3の実施の形態について説明する。なお、第3の実施の形態に係る事故情報算出装置の構成は、第1の実施の形態と同様であるため、同一符号を付して説明を省略する。
[Third Embodiment]
Next, a third embodiment will be described. In addition, since the structure of the accident information calculation apparatus which concerns on 3rd Embodiment is the same as that of 1st Embodiment, it attaches | subjects the same code | symbol and abbreviate | omits description.

第3の実施の形態では、特徴量Aを多次元特徴量へ拡張し、特徴量Aをn次元ベクトルとして取得する点が、第1及び第2の実施の形態と異なっている。 The third embodiment is different from the first and second embodiments in that the feature quantity A i is expanded to a multidimensional feature quantity and the feature quantity A i is acquired as an n-dimensional vector.

第3の実施の形態の統計情報データベース24には、車両の運転行動に関する複数次元の特徴量Aの各次元についての各値が得られる確率p(A)を表す運転情報と、複数次元の特徴量Aに起因して生じる事故Xが発生する確率を表す事故統計情報とが格納されている。 The statistical information database 24 according to the third embodiment includes driving information representing the probability p (A i ) that each value of the multi-dimensional feature value A i relating to the driving behavior of the vehicle is obtained, and multi-dimensional And accident statistical information representing the probability of occurrence of an accident X i caused by the feature amount A i of the.

また、第3の実施の形態に係る事故情報算出装置のコンピュータでは、特徴量取得部28は、複数種類の事故Xの各々について、対象ドライバの運転行動に関する複数次元の特徴量Aの値を取得する。
例えば、特徴量取得部28は、上記表1の脇見による事故に対して、脇見運転距離だけでなく、車速センサ16によって検出された車両の速度、前方カメラ20によって撮像された前方画像から得られる歩行者の有無等を、対象ドライバの運転行動に関する複数次元の特徴量A(=(Ai,1,Ai,2,・・・,Ai,n))の各次元の値として取得する。
また、上記表1の追突事故に対して、車間距離だけでなく、ドライバカメラ22から取得したドライバ画像から検出したドライバの閉眼度合い等を、対象ドライバの運転行動に関する複数次元の特徴量A(=(Ai,1,Ai,2,・・・,Ai,n))の各次元の値として取得してもよい。
Further, in the computer of the accident information calculation apparatus according to the third embodiment, the feature amount acquisition unit 28 has a value of the multi-dimensional feature amount A i regarding the driving behavior of the target driver for each of the plurality of types of accidents X i. To get.
For example, the feature amount acquisition unit 28 is obtained from the vehicle speed detected by the vehicle speed sensor 16 and the front image captured by the front camera 20 in addition to the sideward driving distance for the accident caused by the side look in Table 1 above. The presence / absence of a pedestrian and the like is acquired as a value of each dimension of a multidimensional feature amount A i (= (A i, 1 , A i, 2 ,..., A i, n )) related to the driving behavior of the target driver. To do.
In addition, for the rear-end collision in Table 1, not only the inter-vehicle distance but also the degree of driver's eye closed detected from the driver image acquired from the driver camera 22, the multi-dimensional feature amount A i ( = (A i, 1 , A i, 2 ,..., A i, n )) may be acquired as values of each dimension.

算出部32は、統計情報取得部30によって取得された運転情報及び事故統計情報と、複数次元の特徴量Aの各次元についての事故Xが発生した条件下での当該次元の各値が得られる確率に関する制約条件とに基づいて、特徴量Aの各次元の各値に対して事故Xが発生する確率を近似した予め定められた関数で用いるパラメータを算出し、算出されたパラメータを用いた関数を算出する。 The calculation unit 32 calculates the driving information and the accident statistical information acquired by the statistical information acquisition unit 30, and each value of the dimension under the condition in which the accident X i occurs for each dimension of the multi-dimensional feature value A i. Parameters used in a predetermined function approximating the probability of occurrence of an accident X i for each value of each dimension of the feature value A i are calculated on the basis of the constraint conditions relating to the obtained probability, and the calculated parameters The function using is calculated.

算出部32は、ベイズの定理を用い、かつ複数次元の特徴量Aの各次元が互いに独立又は従属であると仮定して、パラメータを算出する。
具体的には、複数次元の特徴量Aの各次元の値は、互いに独立又は互いに従属して得られると仮定する。特徴量取得部28によって取得される複数次元の特徴量Aが、A=(Ai,1,Ai,2,・・・,Ai,n)というようにn次元である場合には、各特徴Ai,1, Ai,2,・・・, Ai,nを独立または従属と仮定することで、上記式(7)を、以下の式(18)又は式(19)に示すように近似し、さらに以下の式(20)を用いることで、n個の未知数が算出可能となる。以下の式(20)は、複数次元の特徴量Aの各次元についての事故Xが発生した条件下での当該次元の各値が得られる確率に関する制約条件の一例である。
The calculation unit 32 calculates the parameters using Bayes' theorem and assuming that the dimensions of the multi-dimensional feature quantity A i are independent or dependent from each other.
Specifically, it is assumed that the values of each dimension of the multi-dimensional feature quantity A i are obtained independently or subordinate to each other. When the multi-dimensional feature quantity A i acquired by the feature quantity acquisition unit 28 is n-dimensional such as A i = (A i, 1 , A i, 2 ,..., A i, n ). , each feature a i, 1, a i, 2, ···, by assuming a i, independently or dependent of n, the above formula (7), the following equation (18) or (19) By using the following approximation and further using the following equation (20), n unknowns can be calculated. The following equation (20) is an example of a constraint condition regarding the probability that each value of the dimension is obtained under the condition in which the accident X i has occurred for each dimension of the multi-dimensional feature value A i .

したがって、p(X|A)の未知のパラメータがn個以下になるようにp(X|A)を定義すればよく、例えばp(X|A)を以下の式(21)の多変量正規分布の関数で近似する方法が考えられる。 Thus, p may be defined to | | (A i X i) , for example p unknown parameters (X i A i) is p such that n or less (X i | A i) the following equation ( A method of approximating with the function of the multivariate normal distribution of 21) is conceivable.

はパラメータを表し、Σはパラメータ行列を表し、Arefは予め設定したベクトルを表す。なお、第3の実施の形態で求めることができるパラメータの数はn個であるため、上記式(21)のパラメータの数がn個以下になるように、一部のパラメータを経験的に定めたり、別の情報に基づいて予め定めておく。以上のようにして、多次元特徴量に起因する事故に対しても対応可能になると考えられる。 v i represents a parameter, Σ i represents a parameter matrix, and A ref represents a preset vector. Since the number of parameters that can be obtained in the third embodiment is n, some parameters are determined empirically so that the number of parameters in the above equation (21) is n or less. Or predetermined based on other information. As described above, it is considered that it is possible to cope with an accident caused by a multidimensional feature amount.

運転行動評価部34は、特徴量取得部28によって取得された特徴量の値、及び算出部32によって算出された複数種類の事故の各々について算出されたガウス関数に基づいて、複数種類の事故の各々について当該事故Xが発生しない確率(1−p(X|A))を算出する。そして、運転行動評価部34は、対象ドライバの車両の運転行動について、複数種類の事故の各々に対して算出される事故が発生しない確率(1−p(X|A))の積を、事故が起こらない確率を表す無事故確率psafeとして算出する。 The driving behavior evaluation unit 34 calculates the plurality of types of accidents based on the feature value acquired by the feature amount acquisition unit 28 and the Gaussian function calculated for each of the plurality of types of accidents calculated by the calculation unit 32. The probability (1-p (X i | A i )) that the accident X i does not occur is calculated for each. Then, the driving behavior evaluation unit 34 calculates the product of the probability of occurrence of an accident (1-p (X i | A i )) calculated for each of a plurality of types of accidents regarding the driving behavior of the target driver's vehicle. , The accident-free probability psafe representing the probability that no accident will occur is calculated.

なお、第3の実施の形態に係る事故情報算出装置の他の構成及び作用については、第1の実施の形態と同様であるため、説明を省略する。   In addition, about the other structure and effect | action of the accident information calculation apparatus which concern on 3rd Embodiment, since it is the same as that of 1st Embodiment, description is abbreviate | omitted.

このように、第3の実施の形態の事故情報算出装置によれば、運転情報及び事故統計情報と、複数次元の特徴量の各次元についての事故が発生した条件下での当該次元の各値が得られる確率に関する制約条件とに基づいて、特徴量の各次元の各値に対して事故が発生する確率を近似した予め定められた関数で用いるパラメータを算出し、算出されたパラメータを用いた関数を算出し、対象ドライバの車両の運転行動に関する複数の特徴量の各々の値を取得し、複数の特徴量の各々の値、及び算出された複数種類の事故の各々について算出された関数に基づいて、複数種類の事故の各々について事故が発生しない確率を算出し、対象ドライバの車両の運転行動について、無事故確率として評価することにより、複数次元の特徴量に対してドライバの運転を診断することができる。   As described above, according to the accident information calculation apparatus of the third embodiment, each value of the dimension under the condition in which the accident occurs with respect to each dimension of the driving information, the accident statistical information, and the multi-dimensional feature amount. Based on the constraint on the probability of obtaining a parameter, a parameter used in a predetermined function approximating the probability of occurrence of an accident is calculated for each value of each dimension of the feature value, and the calculated parameter is used. The function is calculated, each value of the plurality of feature values related to the driving behavior of the vehicle of the target driver is obtained, and each of the plurality of feature values and the calculated function for each of the plurality of types of accidents are obtained. Based on this, the probability that no accident will occur for each of the multiple types of accidents is calculated, and the driving behavior of the target driver's vehicle is evaluated as an accident-free probability. It is possible to diagnose the operation.

なお、第3の実施の形態では、上記式(21)の多変量正規分布の関数に代えて、以下の式(22)の線形関数、式(23)のシグモイド関数で近似することもできる。   In the third embodiment, instead of the function of the multivariate normal distribution of the above equation (21), approximation can be performed by a linear function of the following equation (22) and a sigmoid function of the equation (23).

及びwi,jはパラメータを表す。なお、上記式(21)と同様に、第3の実施の形態で求めることができるパラメータの数はn個であるため、上記式(22)〜式(23)のパラメータの数がn個以下になるように、一部のパラメータを経験的に定めたり、別の情報に基づいて予め定めておく。 v i and w i, j represent parameters. Since the number of parameters that can be obtained in the third embodiment is n as in the above equation (21), the number of parameters in the above equations (22) to (23) is n or less. Therefore, some parameters are determined empirically or in advance based on other information.

[第4の実施の形態]
次に、第4の実施の形態について説明する。なお、第4の実施の形態に係る事故情報算出装置の構成は、第1の実施の形態と同様であるため、同一符号を付して説明を省略する。
[Fourth Embodiment]
Next, a fourth embodiment will be described. In addition, since the structure of the accident information calculation apparatus which concerns on 4th Embodiment is the same as that of 1st Embodiment, it attaches | subjects the same code | symbol and abbreviate | omits description.

第4の実施の形態では、時々刻々と変化する特徴量Aに対して無事故確率の算出を行う点が、第1〜第3の実施の形態と異なっている。 The fourth embodiment is different from the first to third embodiments in that the accident-free probability is calculated for the feature value A i that changes every moment.

上記第1〜第3の実施の形態で算出したpsafe(A,A,・・・)は、ある瞬間の状態(A,A,・・・)が一定期間持続した際の事故を起こさない確率を意味している。例えば、p(X)を「l[km]走行当たりの事故確率」とした場合、psafe(A,A,・・・)は「状態(A,A,・・・)を維持してl[km]走行した場合の事故確率」を意味する。
しかし、実際には同じ運転状態を持続する事は稀であるため、第4の実施の形態では、時々刻々と変化する運転状態に対しての無事故確率の算出を行う。
The p safe (A 1 , A 2 ,...) Calculated in the first to third embodiments is obtained when the state (A 1 , A 2 ,. It means the probability of not causing an accident. For example, when p (X i ) is “1 [km] accident probability per driving”, p safe (A 1 , A 2 ,...) Is “state (A 1 , A 2 ,...)”. This means the accident probability when driving 1 [km] while maintaining
However, since it is rare that the same driving state is actually maintained, in the fourth embodiment, the accident-free probability is calculated for the driving state that changes every moment.

はじめに、任意の物理量Lを定義する。ただし同時に、特徴量Aに対する条件付き事故確率p(X|A)は「特徴量Aの値が一定のまま対象の運転がLを満たすまでに事故が発生する確率」として定義される。例えば、Lは走行距離や走行時間等が該当し、Lを走行距離とするとp(X|A)は「特徴量Aの値が一定のままL[km]の走行を満たすまでに事故が発生する確率」となる。ここで、物理量Lに対して特徴量Aの状態である区間の総和をl(A)とすると、対象の運転がl(A)を満たすまでに事故Xが発生しない確率(事故Xに対する無事故確率)qsafe(X|A,l(A))を定義したとき、以下の式(24)が成立する。 First, an arbitrary physical quantity L i is defined. But simultaneously, conditional accident probability p for the feature amount A i (X i | A i) is defined as the "probability that remain with the operational values of certain characteristic amount A i is an accident occurs until satisfy L i" Is done. For example, L i is applicable such as travel distance and travel time, when the travel distance L i p (X i | A i) is the running of the "while the value of the feature amount A i is constant L i [km] The probability that an accident will occur before it is satisfied. Here, if the sum of the physical quantity L i is the state of the feature quantity A i relative section and l i (A i), accident X i is not generated until the operation of the subject satisfies the l i (A i) When the probability (accident-free probability for accident X i ) q safe (X i | A i , l i (A i )) is defined, the following equation (24) is established.

したがって、特徴量Aの値が時々刻々と変化する中、対象の運転がLを満たすまでの無事故確率qsafeは、以下の式(25)のように与えられる。 Accordingly, while the value of the feature value A i changes every moment, the accident-free probability q safe until the target driving satisfies L i is given by the following equation (25).

ここで、SはAの非可算集合(Aの取り得る全ての値の集合)である。 Here, S i is a non-countable set of A i (a set of all possible values of A i).

例えば、任意の物理量Lとして走行距離を用いる場合には、特徴量取得部28は、対象ドライバの車両の運転行動に関する特徴量Aの値、及び特徴量Aの値が得られているときの走行距離l(A)を取得する。
そして、運転行動評価部34は、特徴量取得部28によって取得された特徴量Aの値及び走行距離l(A)と、算出部32によって複数種類の事故の各々について算出された、予め定められた走行距離Lまでに特徴量Aの各値に対して事故Xが発生する確率p(X|A)を近似した関数とに基づいて、対象ドライバの車両の運転行動について、複数種類の事故の各々に対して算出される事故が発生しない確率(1−p(X|A))の積を、無事故確率qsafe(A,A,・・・)を算出する。
For example, when the travel distance is used as the arbitrary physical quantity L i , the feature quantity acquisition unit 28 obtains the value of the feature quantity A i and the value of the feature quantity A i regarding the driving behavior of the target driver's vehicle. The travel distance l i (A i ) is acquired.
The driving behavior evaluation unit 34 calculates the value of the feature amount A i and the travel distance l i (A i ) acquired by the feature amount acquisition unit 28 and each of the plurality of types of accidents by the calculation unit 32. probability p accidents X i is generated by the travel distance L i which is predetermined for each value of the characteristic quantity a i |, based on the function obtained by approximating (X i a i), the operation of the vehicle target driver For the behavior, the product of the probability of occurrence of an accident (1-p (X i | A i )) calculated for each of a plurality of types of accidents is used as the accident-free probability q safe (A 1 , A 2 ,... ) Is calculated.

なお、第3の実施の形態に係る事故情報算出装置の他の構成及び作用については、第1〜第3の実施の形態と同様であるため、説明を省略する。   In addition, about the other structure and effect | action of the accident information calculation apparatus which concern on 3rd Embodiment, since it is the same as that of 1st-3rd Embodiment, description is abbreviate | omitted.

第4の実施の形態の事故情報算出装置によれば、取得された特徴量の値、及び継続時間又は走行時間と、予め定められた継続時間又は走行時間までに特徴量の各値に対して事故が発生する確率を近似した関数とに基づいて、対象ドライバの車両の運転行動について、複数種類の事故の各々に対して算出される事故が発生しない確率の積を、無事故確率として評価することにより、時々刻々と変化する特徴量に対して無事故確率の算出をすることができる。   According to the accident information calculation apparatus of the fourth embodiment, the acquired feature value, the duration or the running time, and each value of the feature amount up to a predetermined duration or the running time. Based on a function that approximates the probability that an accident will occur, the product of the probability that no accident will occur for each type of accident will be evaluated as the accident-free probability for the driving behavior of the target driver's vehicle Thus, it is possible to calculate the accident-free probability for the feature quantity that changes every moment.

なお、上記の第4の実施の形態では、任意の物理量Lとして走行距離を用いた場合を例に説明したが、任意の物理量Lとして走行時間を用いるようにしてもよい。この場合には、p(X|A)は「特徴量Aの値が一定のまま継続した走行時間がL時間の走行を満たすまでに事故が発生する確率」となる。特徴量取得部28は、対象ドライバの車両の運転行動に関する特徴量Aの値、及び特徴量Aの値の継続時間l(A)を取得する。そして、運転行動評価部34は、特徴量取得部28によって取得された特徴量Aの値及び継続時間l(A)と、算出部32によって複数種類の事故の各々について算出された、予め定められた継続時間Lまでに特徴量Aの各値に対して事故Xが発生する確率p(X|A)を近似した関数とに基づいて、対象ドライバの車両の運転行動について、複数種類の事故の各々に対して算出される事故が発生しない確率(1−p(X|A))の積を、無事故確率qsafe(A,A,・・・)を算出する。 In the fourth embodiment described above, the case of using the travel distance as any physical quantity L i has been described as an example, may be used travel time as any physical quantity L i. In this case, p (X i | A i ) is "the probability that the travel time value of the feature amount A i is continued remained constant accident occurs until meeting the travel of L i time". The feature amount acquisition unit 28 acquires the value of the feature amount A i related to the driving behavior of the vehicle of the target driver and the duration l i (A i ) of the value of the feature amount A i . Then, the driving behavior evaluation unit 34 calculates the value of the feature amount A i acquired by the feature amount acquisition unit 28 and the duration l i (A i ), and each of a plurality of types of accidents by the calculation unit 32. Based on a function that approximates the probability p (X i | A i ) of occurrence of an accident X i for each value of the feature value A i by a predetermined duration L i , the driving of the target driver's vehicle For the behavior, the product of the probability of occurrence of an accident (1-p (X i | A i )) calculated for each of a plurality of types of accidents is used as the accident-free probability q safe (A 1 , A 2 ,... ) Is calculated.

また、上記の実施の形態では、事故情報算出装置が車両に搭載されている場合を例に説明したが、これに限定されるものではなく、例えばサーバ等によって、事故情報算出装置が実現されてもよい。この場合には、複数のドライバの運転行動に関する特徴量の値を収集して、車両の運転行動に関する複数の特徴量Aについての各値が得られる確率p(A)を表す運転情報を求めるようにしてもよい。また、リアルタイムに逐次特徴量を取得し、対象ドライバの運転行動を評価してもよい。 In the above embodiment, the case where the accident information calculation device is mounted on the vehicle has been described as an example. However, the present invention is not limited to this. For example, the accident information calculation device is realized by a server or the like. Also good. In this case, characteristic values relating to driving behaviors of a plurality of drivers are collected, and driving information representing the probability p (A i ) at which each value for the plurality of characteristic values A i relating to the driving behavior of the vehicle is obtained. You may make it ask. Further, the feature amount may be sequentially acquired in real time to evaluate the driving behavior of the target driver.

また、上記の実施の形態では、取得された特徴量の値に対して事故が発生する確率を算出し、算出された事故が発生する確率に基づいて、事故が起こらない確率を表す無事故確率を算出する事故情報算出装置に本発明を適用させた場合を例に説明したが、これに限定されるものではない。
例えば、統計情報取得部30は、車両の運転行動に関する特徴量の各値が得られる度合いを表す運転情報と、特徴量に起因して生じる事故が発生する度合いを表す事故統計情報とを取得し、算出部32は、統計情報取得部30によって取得された運転情報及び事故統計情報と、事故が発生した条件下での特徴量の各値が得られる度合いに関する制約条件とに基づいて、特徴量の各値に対して事故が発生する度合いを算出してもよい。この場合には、対象ドライバの無事故確率を算出するわけではなく、特徴量の各値に対して事故が発生する度合いを算出することができる。
また、特徴量の各値に対して事故が発生する度合いを算出する場合には、算出部32は、統計情報取得部30によって取得された運転情報及び事故統計情報と、事故が発生した条件下での特徴量の各値が得られる度合いに関する制約条件とに基づいて、特徴量の各値に対して事故が発生する度合いを近似した予め定められた関数で用いるパラメータを算出し、算出されたパラメータを用いた関数を、特徴量の各値に対して前記事故が発生する度合いとして算出してもよい。
Further, in the above embodiment, the probability of occurrence of an accident is calculated for the acquired feature value, and the accident-free probability representing the probability that no accident will occur is calculated based on the calculated probability of occurrence of the accident. Although the case where the present invention is applied to the accident information calculation apparatus to be calculated has been described as an example, the present invention is not limited to this.
For example, the statistical information acquisition unit 30 acquires driving information indicating the degree to which each value of the characteristic amount related to the driving behavior of the vehicle is obtained, and accident statistical information indicating the degree of occurrence of an accident caused by the characteristic amount. The calculation unit 32 is based on the driving information and accident statistical information acquired by the statistical information acquisition unit 30 and the restriction condition regarding the degree to which each value of the characteristic value under the condition where the accident occurred is obtained. The degree of occurrence of an accident may be calculated for each value. In this case, the accident-free probability of the target driver is not calculated, but the degree of occurrence of an accident can be calculated for each feature value.
Further, when calculating the degree of occurrence of an accident for each value of the feature amount, the calculation unit 32 calculates the driving information and accident statistical information acquired by the statistical information acquisition unit 30 and the conditions under which the accident occurred. The parameters used in a predetermined function that approximates the degree of occurrence of an accident with respect to each value of the feature amount are calculated based on the restriction condition regarding the degree that each value of the feature amount is obtained at A function using parameters may be calculated as the degree of occurrence of the accident for each value of the feature value.

または、特徴量取得部28は、対象ドライバの車両の運転行動に関する特徴量の値を取得し、運転行動評価部34は、特徴量取得部28によって取得された特徴量の値、及び算出部32によって算出された特徴量の各値に対して事故が発生する度合いに基づいて、対象ドライバの車両の運転行動について、事故が発生する度合い又は事故が発生しない度合いを評価してもよい。この場合には、確率を算出するわけではなく、事故が発生する度合い又は事故が発生しない度合いを算出することができる。   Alternatively, the feature amount acquisition unit 28 acquires the value of the feature amount related to the driving behavior of the vehicle of the target driver, and the driving behavior evaluation unit 34 calculates the value of the feature amount acquired by the feature amount acquisition unit 28 and the calculation unit 32. The degree of occurrence of an accident or the degree of no occurrence of an accident may be evaluated for the driving behavior of the target driver's vehicle based on the degree of occurrence of an accident for each value of the feature amount calculated by the above. In this case, the probability is not calculated, but the degree at which an accident occurs or the degree at which no accident occurs can be calculated.

上記のように、特徴量の各値に対して事故が発生する度合い又は事故が発生しない度合いを算出する場合には、確率を算出しないため、上記式(7)に示したベイズの定理は、確率法則として用いずに、非線形関数同士の関係として用いる。   As described above, when calculating the degree of occurrence of an accident or the degree of no occurrence of an accident for each value of the feature value, the probability is not calculated, so the Bayes' theorem shown in the above equation (7) is It is not used as a probability law but as a relationship between nonlinear functions.

また、上記の実施の形態では、無事故確率を算出する場合を例に説明したが、これに限定されるものではなく、例えば、統計情報取得部30は、車両の運転行動に関する特徴量の各値が得られる確率を表す運転情報と、特徴量に起因して生じる事故が発生する確率を表す事故統計情報とを取得し、算出部32は、統計情報取得部30によって取得された運転情報及び事故統計情報と、事故が発生した条件下での特徴量の各値が得られる確率に関する制約条件とに基づいて、特徴量の各値に対して事故が発生する確率を近似した予め定められた関数で用いるパラメータを算出し、算出されたパラメータを用いた関数のみを算出してもよい。   Further, in the above-described embodiment, the case where the accident-free probability is calculated has been described as an example. However, the present invention is not limited to this. For example, the statistical information acquisition unit 30 includes each value of the feature amount related to the driving behavior of the vehicle. Driving information indicating the probability of obtaining an accident and accident statistical information indicating the probability of occurrence of an accident caused by the feature amount, the calculation unit 32 acquires the driving information and the accident acquired by the statistical information acquisition unit 30 Predetermined function that approximates the probability that an accident will occur for each value of the feature based on statistical information and the constraint on the probability that each value of the feature will be obtained under the condition where the accident occurred It is also possible to calculate the parameters used in step 1 and calculate only the function using the calculated parameters.

また、上記の実施の形態では、特徴量取得部28によって取得された特徴量の各値に対して事故が発生する確率を算出する場合を例に説明したが、特徴量取得部28によって取得された特徴量の各値に対して事故が発生しない確率を算出してもよい。   In the above embodiment, the case where the probability of occurrence of an accident is calculated for each value of the feature amount acquired by the feature amount acquisition unit 28 has been described as an example. The probability that an accident will not occur may be calculated for each value of the feature amount.

また、上記の実施の形態では、横方向加速度センサ12によって検出された横方向加速度と、前後方向加速度センサ14によって検出された前後方向加速度と、一時停止時の通過速度と、車間距離センサ18によって検出された車間距離と、一時停止時の通過速度と、脇見運転距離と、赤信号時の交差点通過時間とを、ドライバの運転行動に関する複数の特徴量として取得する場合を例に説明したが、これに限定されるものではなく、他の特徴量を取得してもよい。   In the above embodiment, the lateral acceleration detected by the lateral acceleration sensor 12, the longitudinal acceleration detected by the longitudinal acceleration sensor 14, the passing speed at the time of pause, and the inter-vehicle distance sensor 18 are used. Although the case where the detected inter-vehicle distance, the passing speed at the time of the temporary stop, the armpit driving distance, and the intersection passing time at the time of the red light are acquired as a plurality of feature amounts is described as an example, However, the present invention is not limited to this, and other feature amounts may be acquired.

なお、本発明のプログラムは、記録媒体に格納して提供することができる。   The program of the present invention can be provided by being stored in a recording medium.

10 事故情報算出装置
12 横方向加速度センサ
14 前後方向加速度センサ
16 車速センサ
18 車間距離センサ
20 前方カメラ
22 ドライバカメラ
24 統計情報データベース
26 コンピュータ
28 特徴量取得部
30 統計情報取得部
32 算出部
34 運転行動評価部
40 出力装置
DESCRIPTION OF SYMBOLS 10 Accident information calculation apparatus 12 Lateral direction acceleration sensor 14 Longitudinal direction acceleration sensor 16 Vehicle speed sensor 18 Inter-vehicle distance sensor 20 Front camera 22 Driver camera 24 Statistical information database 26 Computer 28 Feature-value acquisition part 30 Statistical information acquisition part 32 Calculation part 34 Driving behavior Evaluation unit 40 output device

Claims (16)

車両の運転行動に関する特徴量の各値が得られる度合いを表す運転情報と、前記特徴量に起因して生じる事故が発生する度合いを表す事故統計情報とを取得する統計情報取得手段と、
前記統計情報取得手段によって取得された前記運転情報及び前記事故統計情報と、前記事故が発生した条件下での前記特徴量の各値が得られる度合いに関する制約条件とに基づいて、前記特徴量の各値に対して前記事故が発生する度合いを算出する算出手段と、
を含む事故情報算出装置。
Statistical information acquisition means for acquiring driving information indicating the degree to which each value of the characteristic amount related to the driving behavior of the vehicle is obtained, and accident statistical information indicating the degree of occurrence of an accident caused by the characteristic amount;
Based on the driving information and the accident statistical information acquired by the statistical information acquisition means, and a constraint condition regarding the degree to which each value of the characteristic value under the condition in which the accident has occurred is obtained. Calculating means for calculating the degree of occurrence of the accident for each value;
Accident information calculation device.
前記算出手段は、前記統計情報取得手段によって取得された前記運転情報及び前記事故統計情報と、前記事故が発生した条件下での前記特徴量の各値が得られる度合いに関する制約条件とに基づいて、前記特徴量の各値に対して前記事故が発生する度合いを近似した予め定められた関数で用いるパラメータを算出し、前記算出されたパラメータを用いた前記関数を、前記特徴量の各値に対して前記事故が発生する度合いとして算出する請求項1記載の事故情報算出装置。   The calculation means is based on the driving information and the accident statistical information acquired by the statistical information acquisition means, and a constraint condition relating to the degree to which each value of the feature value is obtained under the condition where the accident has occurred. Calculating a parameter used in a predetermined function approximating the degree of occurrence of the accident with respect to each value of the feature value, and converting the function using the calculated parameter to each value of the feature value. The accident information calculation apparatus according to claim 1, wherein the accident information calculation apparatus calculates the degree of occurrence of the accident. 前記算出手段は、前記統計情報取得手段によって取得された前記運転情報及び前記事故統計情報と、前記事故が発生した条件下での前記特徴量の各値が得られる度合いに関する制約条件とに基づいて、ベイズの定理を用いて、前記パラメータを算出する請求項2記載の事故情報算出装置。   The calculation means is based on the driving information and the accident statistical information acquired by the statistical information acquisition means, and a constraint condition relating to the degree to which each value of the feature value is obtained under the condition where the accident has occurred. The accident information calculation apparatus according to claim 2, wherein the parameter is calculated using Bayes' theorem. 対象ドライバの車両の運転行動に関する特徴量の値を取得する特徴量取得手段と、
前記特徴量取得手段によって取得された前記特徴量の値、及び前記算出手段によって算出された前記特徴量の各値に対して前記事故が発生する度合いに基づいて、前記対象ドライバの車両の運転行動について、前記事故が発生する度合い又は前記事故が発生しない度合いを評価する運転行動評価手段と、を更に含む請求項1〜3の何れか1項記載の事故情報算出装置。
Feature amount acquisition means for acquiring a value of a feature amount related to the driving behavior of the vehicle of the target driver;
Based on the value of the feature quantity acquired by the feature quantity acquisition means and the degree of occurrence of the accident with respect to each value of the feature quantity calculated by the calculation means, the driving behavior of the vehicle of the target driver The accident information calculation device according to any one of claims 1 to 3, further comprising: a driving behavior evaluation unit that evaluates the degree of occurrence of the accident or the degree of occurrence of the accident.
前記統計情報取得手段は、車両の運転行動に関する特徴量の各値が得られる確率を表す前記運転情報と、前記特徴量に起因して生じる事故が発生する確率を表す前記事故統計情報とを取得し、
前記算出手段は、前記統計情報取得手段によって取得された前記運転情報及び前記事故統計情報と、前記事故が発生した条件下での前記特徴量の各値が得られる確率に関する前記制約条件とに基づいて、前記特徴量の各値に対して前記事故が発生する確率を近似した予め定められた関数で用いるパラメータを算出し、算出された前記パラメータを用いた前記関数を算出する請求項2〜請求項4の何れか1項記載の事故情報算出装置。
The statistical information acquisition means acquires the driving information indicating a probability that each value of a feature amount related to driving behavior of a vehicle is obtained, and the accident statistical information indicating a probability that an accident caused by the feature amount occurs. And
The calculation means is based on the driving information and the accident statistical information acquired by the statistical information acquisition means, and the constraint condition regarding the probability that each value of the feature value is obtained under the condition where the accident has occurred. And calculating a parameter to be used for a predetermined function approximating the probability of occurrence of the accident for each value of the feature amount, and calculating the function using the calculated parameter. 5. The accident information calculation apparatus according to any one of items 4.
前記統計情報取得手段は、車両の運転行動に関する複数次元の特徴量の各次元についての各値が得られる確率を表す前記運転情報と、前記複数次元の特徴量に起因して生じる事故が発生する確率を表す前記事故統計情報とを取得し、
前記算出手段は、前記統計情報取得手段によって取得された前記運転情報及び前記事故統計情報と、前記複数次元の特徴量の各次元についての前記事故が発生した条件下での前記次元の各値が得られる確率に関する前記制約条件とに基づいて、前記特徴量の各次元の各値に対して前記事故が発生する確率を近似した予め定められた関数で用いるパラメータを算出し、算出された前記パラメータを用いた前記関数を算出する請求項2〜請求項4の何れか1項記載の事故情報算出装置。
The statistical information acquisition means generates the driving information indicating the probability of obtaining each value of each dimension of the multi-dimensional feature quantity related to the driving behavior of the vehicle and an accident caused by the multi-dimensional feature quantity. The accident statistics information representing the probability,
The calculation means includes the driving information and the accident statistical information acquired by the statistical information acquisition means, and each value of the dimension under the condition in which the accident has occurred for each dimension of the multi-dimensional feature amount. Based on the constraint condition regarding the obtained probability, a parameter used in a predetermined function approximating the probability of occurrence of the accident is calculated for each value of each dimension of the feature value, and the calculated parameter The accident information calculation apparatus according to any one of claims 2 to 4, wherein the function is calculated using the function.
前記統計情報取得手段は、車両の運転行動に関する複数の特徴量についての各値が得られる確率を表す前記運転情報と、前記複数の特徴量の各々に起因して生じる複数種類の事故の各々が発生する確率を表す前記事故統計情報とを取得し、
前記算出手段は、前記複数種類の事故の各々について、前記統計情報取得手段によって取得された前記運転情報及び前記事故統計情報と、前記事故が発生した条件下での前記事故に対応する前記特徴量の各値が得られる確率に関する前記制約条件とに基づいて、前記特徴量の各値に対して前記事故が発生する確率を近似した予め定められた関数で用いるパラメータを算出し、算出された前記パラメータを用いた前記関数を算出する請求項5又は請求項6記載の事故情報算出装置。
The statistical information acquisition means includes the driving information indicating the probability that each value for a plurality of feature values related to driving behavior of a vehicle is obtained, and each of a plurality of types of accidents caused by each of the plurality of feature values. The accident statistics information representing the probability of occurrence,
The calculation means includes, for each of the plurality of types of accidents, the driving information and the accident statistical information acquired by the statistical information acquisition means, and the feature amount corresponding to the accident under the condition where the accident has occurred. Based on the constraint condition relating to the probability that each value is obtained, a parameter used in a predetermined function that approximates the probability that the accident will occur for each value of the feature amount is calculated, and the calculated The accident information calculation apparatus according to claim 5 or 6, wherein the function using a parameter is calculated.
対象ドライバの車両の運転行動に関する前記複数の特徴量の各々の値を取得する特徴量取得手段と、
前記特徴量取得手段によって取得された前記複数の特徴量の各々の値、及び前記算出手段によって前記複数種類の事故の各々について算出された前記関数に基づいて、前記複数種類の事故の各々について前記事故が発生しない確率を算出し、前記対象ドライバの車両の運転行動について、前記複数種類の事故の各々に対して算出される前記事故が発生しない確率の積を、事故が起こらない確率を表す無事故確率として評価する運転行動評価手段とを更に含む請求項7記載の事故情報算出装置。
Feature quantity acquisition means for acquiring each value of the plurality of feature quantities relating to the driving behavior of the vehicle of the target driver;
Based on the value of each of the plurality of feature amounts acquired by the feature amount acquisition unit and the function calculated for each of the plurality of types of accidents by the calculation unit, for each of the plurality of types of accidents Calculate the probability that an accident will not occur, and for the driving behavior of the target driver's vehicle, the product of the probability that the accident will not occur calculated for each of the multiple types of accidents represents the probability that no accident will occur The accident information calculation device according to claim 7, further comprising driving behavior evaluation means for evaluating as a probability.
前記特徴量取得手段は、対象ドライバの車両の運転行動に関する前記特徴量の値、及び前記特徴量の値の継続時間を取得し、
前記運転行動評価手段は、前記特徴量取得手段によって取得された前記特徴量の値及び前記継続時間と、前記算出手段によって前記複数種類の事故の各々について算出された、予め定められた継続時間までに前記特徴量の各値に対して前記事故が発生する確率を近似した前記関数とに基づいて、前記対象ドライバの車両の運転行動について、前記複数種類の事故の各々に対して算出される前記事故が発生しない確率の積を、前記無事故確率として評価する
請求項8記載の事故情報算出装置。
The feature amount acquisition means acquires the value of the feature amount related to the driving behavior of the vehicle of the target driver, and the duration of the value of the feature amount,
The driving behavior evaluation means includes a value of the feature amount acquired by the feature amount acquisition means and the duration, and a predetermined duration calculated for each of the plurality of types of accidents by the calculation means. Based on the function that approximates the probability of occurrence of the accident for each value of the feature amount, the driving behavior of the vehicle of the target driver is calculated for each of the plurality of types of accidents The accident information calculation apparatus according to claim 8, wherein a product of probabilities that an accident does not occur is evaluated as the accident-free probability.
前記特徴量取得手段は、対象ドライバの車両の運転行動に関する前記特徴量の値、及び前記特徴量の値が得られているときの走行距離を取得し、
前記運転行動評価手段は、前記特徴量取得手段によって取得された前記特徴量の値及び前記走行距離と、前記算出手段によって前記複数種類の事故の各々について算出された、予め定められた走行距離までに前記特徴量の各値に対して前記事故が発生する確率を近似した前記関数とに基づいて、前記対象ドライバの車両の運転行動について、前記複数種類の事故の各々に対して算出される前記事故が発生しない確率の積を、前記無事故確率として評価する
請求項8記載の事故情報算出装置。
The feature amount acquisition means acquires the value of the feature amount related to the driving behavior of the vehicle of the target driver, and the travel distance when the value of the feature amount is obtained,
The driving behavior evaluation unit is configured to obtain a predetermined travel distance calculated for each of the plurality of types of accidents by the value of the feature amount and the travel distance acquired by the feature amount acquisition unit and the calculation unit. Based on the function that approximates the probability of occurrence of the accident for each value of the feature amount, the driving behavior of the vehicle of the target driver is calculated for each of the plurality of types of accidents The accident information calculation apparatus according to claim 8, wherein a product of probabilities that an accident does not occur is evaluated as the accident-free probability.
前記関数は、前記関数の値が最大となるときの前記特徴量の予め定められた値、又は前記関数の値が最小となるときの前記特徴量の予め定められた値を含む
請求項5〜請求項10の何れか1項記載の事故情報算出装置。
The function includes a predetermined value of the feature amount when the value of the function is maximum, or a predetermined value of the feature amount when the value of the function is minimum. The accident information calculation device according to claim 10.
前記算出手段は、前記統計情報取得手段によって取得された前記運転情報及び前記事故統計情報と、前記複数次元の特徴量の各次元についての前記事故が発生した条件下での前記次元の各値が得られる確率に関する前記制約条件とに基づいて、ベイズの定理を用い、かつ前記複数次元の特徴量の各次元が互いに独立又は従属であると仮定して、前記パラメータを算出する請求項6記載の事故情報算出装置。   The calculation means includes the driving information and the accident statistical information acquired by the statistical information acquisition means, and each value of the dimension under the condition in which the accident has occurred for each dimension of the multi-dimensional feature amount. 7. The parameter according to claim 6, wherein the parameter is calculated on the basis of the constraint on the obtained probability, using Bayes' theorem, and assuming that the dimensions of the feature quantities of the plurality of dimensions are independent or dependent from each other. Accident information calculation device. 前記関数は、前記特徴量の各値に対して前記事故が発生する度合いを近似したガウス関数である請求項2又は請求項3記載の事故情報算出装置。   4. The accident information calculation apparatus according to claim 2, wherein the function is a Gaussian function that approximates the degree of occurrence of the accident for each value of the feature amount. 前記関数は、前記特徴量の各値に対して前記事故が発生する度合いを近似した線形関数である請求項2又は請求項3記載の事故情報算出装置。   4. The accident information calculation apparatus according to claim 2, wherein the function is a linear function that approximates the degree of occurrence of the accident for each value of the feature amount. 前記関数は、前記特徴量の各値に対して前記事故が発生する度合いを近似したシグモイド関数である請求項2又は請求項3記載の事故情報算出装置。   4. The accident information calculation apparatus according to claim 2, wherein the function is a sigmoid function that approximates the degree of occurrence of the accident for each value of the feature amount. コンピュータを、
車両の運転行動に関する特徴量の各値が得られる度合いを表す運転情報と、前記特徴量に起因して生じる事故が発生する度合いを表す事故統計情報とを取得する統計情報取得手段、及び
前記統計情報取得手段によって取得された前記運転情報及び前記事故統計情報と、前記事故が発生した条件下での前記特徴量の各値が得られる度合いに関する制約条件とに基づいて、前記特徴量の各値に対して前記事故が発生する度合いを算出する算出手段
として機能させるためのプログラム。
Computer
Statistical information acquisition means for acquiring driving information indicating the degree to which each value of the characteristic amount related to driving behavior of the vehicle is obtained, and accident statistical information indicating the degree of occurrence of an accident caused by the characteristic amount, and the statistics Each value of the feature value based on the driving information and the accident statistical information acquired by the information acquisition means, and a constraint on the degree to which each value of the feature value is obtained under the condition where the accident occurred A program for functioning as a calculating means for calculating the degree of occurrence of the accident.
JP2014166068A 2014-08-18 2014-08-18 Accident information calculation device and program Expired - Fee Related JP6365107B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2014166068A JP6365107B2 (en) 2014-08-18 2014-08-18 Accident information calculation device and program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2014166068A JP6365107B2 (en) 2014-08-18 2014-08-18 Accident information calculation device and program

Publications (2)

Publication Number Publication Date
JP2016042298A true JP2016042298A (en) 2016-03-31
JP6365107B2 JP6365107B2 (en) 2018-08-01

Family

ID=55592012

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014166068A Expired - Fee Related JP6365107B2 (en) 2014-08-18 2014-08-18 Accident information calculation device and program

Country Status (1)

Country Link
JP (1) JP6365107B2 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106355883A (en) * 2016-10-20 2017-01-25 同济大学 Risk evaluation model-based traffic accident happening probability acquiring method and system
DE102017104519A1 (en) 2016-03-04 2017-09-07 Denso Corporation Rotary electrical machine and method of making the same
JP2019128744A (en) * 2018-01-23 2019-08-01 株式会社ユピテル System and program and the like
JP2019532372A (en) * 2017-04-18 2019-11-07 ベイジン ディディ インフィニティ テクノロジー アンド ディベロップメント カンパニー リミティッド System and method for determining a driver's safety score
WO2020136840A1 (en) 2018-12-27 2020-07-02 パイオニア株式会社 Accident index calculation device, information providing device, content selection device, insurance premium setting device, accident index calculation method, and program

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1195795A (en) * 1997-09-16 1999-04-09 Nippon Telegr & Teleph Corp <Ntt> Voice quality evaluating method and recording medium
JP2006243856A (en) * 2005-03-01 2006-09-14 Hitachi Ltd Operation diagnosis method and its device
JP2006350737A (en) * 2005-06-16 2006-12-28 Yafoo Japan Corp Computer processing method and program
JP2007072393A (en) * 2005-09-09 2007-03-22 Ntt Data Corp Speech recognition device and program
JP2007128169A (en) * 2005-11-01 2007-05-24 Japan Association For Safety Of Hazardous Materials Risk evaluation device and computer program
JP2008165412A (en) * 2006-12-27 2008-07-17 Toshiba Corp Performance-computing device
WO2009113265A1 (en) * 2008-03-11 2009-09-17 パナソニック株式会社 Tag sensor system and sensor device, and object position estimating device and object position estimating method
JP2009282873A (en) * 2008-05-26 2009-12-03 Hitachi Ltd Accident warning reporting system using bayesian network
JP2010097392A (en) * 2008-10-16 2010-04-30 Chugoku Electric Power Co Inc:The Facility deterioration prediction system and facility deterioration prediction method
JP2010134745A (en) * 2008-12-05 2010-06-17 Fujitsu Ten Ltd Drive recorder and video recording method for the same
JP2013190909A (en) * 2012-03-13 2013-09-26 Toshiba Corp System fatigue degree evaluation device and method for the same

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1195795A (en) * 1997-09-16 1999-04-09 Nippon Telegr & Teleph Corp <Ntt> Voice quality evaluating method and recording medium
JP2006243856A (en) * 2005-03-01 2006-09-14 Hitachi Ltd Operation diagnosis method and its device
JP2006350737A (en) * 2005-06-16 2006-12-28 Yafoo Japan Corp Computer processing method and program
JP2007072393A (en) * 2005-09-09 2007-03-22 Ntt Data Corp Speech recognition device and program
JP2007128169A (en) * 2005-11-01 2007-05-24 Japan Association For Safety Of Hazardous Materials Risk evaluation device and computer program
JP2008165412A (en) * 2006-12-27 2008-07-17 Toshiba Corp Performance-computing device
WO2009113265A1 (en) * 2008-03-11 2009-09-17 パナソニック株式会社 Tag sensor system and sensor device, and object position estimating device and object position estimating method
JP2009282873A (en) * 2008-05-26 2009-12-03 Hitachi Ltd Accident warning reporting system using bayesian network
JP2010097392A (en) * 2008-10-16 2010-04-30 Chugoku Electric Power Co Inc:The Facility deterioration prediction system and facility deterioration prediction method
JP2010134745A (en) * 2008-12-05 2010-06-17 Fujitsu Ten Ltd Drive recorder and video recording method for the same
JP2013190909A (en) * 2012-03-13 2013-09-26 Toshiba Corp System fatigue degree evaluation device and method for the same

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102017104519A1 (en) 2016-03-04 2017-09-07 Denso Corporation Rotary electrical machine and method of making the same
CN106355883A (en) * 2016-10-20 2017-01-25 同济大学 Risk evaluation model-based traffic accident happening probability acquiring method and system
CN106355883B (en) * 2016-10-20 2019-11-29 同济大学 A kind of method, apparatus and system obtaining traffic accident probability of happening
JP2019532372A (en) * 2017-04-18 2019-11-07 ベイジン ディディ インフィニティ テクノロジー アンド ディベロップメント カンパニー リミティッド System and method for determining a driver's safety score
JP2019128744A (en) * 2018-01-23 2019-08-01 株式会社ユピテル System and program and the like
WO2020136840A1 (en) 2018-12-27 2020-07-02 パイオニア株式会社 Accident index calculation device, information providing device, content selection device, insurance premium setting device, accident index calculation method, and program

Also Published As

Publication number Publication date
JP6365107B2 (en) 2018-08-01

Similar Documents

Publication Publication Date Title
JP6365107B2 (en) Accident information calculation device and program
Wang et al. A safety-based approaching behavioural model with various driving characteristics
CN111489588B (en) Vehicle driving risk early warning method and device, equipment and storage medium
US10220841B2 (en) Method and system for assisting a driver of a vehicle in driving the vehicle, vehicle and computer program
CN108062095B (en) Object tracking using sensor fusion within a probabilistic framework
WO2018103313A1 (en) Traffic accident occurrence risk prediction method, device and system
Wang et al. What is the level of volatility in instantaneous driving decisions?
Klingelschmitt et al. Combining behavior and situation information for reliably estimating multiple intentions
WO2020042984A1 (en) Vehicle behavior detection method and apparatus
CN112487617A (en) Collision model-based risk prevention method, device, equipment and storage medium
JP7413503B2 (en) Evaluating vehicle safety performance
Noh et al. Risk assessment for automatic lane change maneuvers on highways
EP3940487B1 (en) Estimation of probability of collision with increasing severity level for autonomous vehicles
KR20200078274A (en) Method and system for evaluating safety operation index using vehicle driving information collection device
US11511737B2 (en) Apparatus and method for processing vehicle signals to compute a behavioral hazard measure
Martinsson et al. Clustering vehicle maneuver trajectories using mixtures of hidden markov models
Alpar et al. Intelligent collision warning using license plate segmentation
Ucar et al. Abnormal driving behavior detection system
Weng et al. A finite-sampling, operational domain specific, and provably unbiased connected and automated vehicle safety metric
Chang et al. Driving safety monitoring and warning for connected and automated vehicles via edge computing
US20220126818A1 (en) Systems and methods for identifying high-risk driving situations from driving data
CN114446042B (en) Method, device, equipment and storage medium for early warning traffic accidents
JP6358051B2 (en) Transition prediction data generation device and transition prediction device
Lv et al. A novel rear-end collision warning algorithm in VANET
DE112020001776T5 (en) VEHICLE CONTROL SYSTEM

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20170306

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20171220

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180109

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180312

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20180618

R150 Certificate of patent or registration of utility model

Ref document number: 6365107

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees