JP5214355B2

JP5214355B2 - 車両走行軌跡観測システム、車両走行軌跡観測方法、およびそのプログラム

Info

Publication number: JP5214355B2
Application number: JP2008183228A
Authority: JP
Inventors: 弘和赤羽; 粋史小宮; 良太堀口
Original assignee: I-TRANSPORT LAB. CO., LTD.
Current assignee: I-TRANSPORT LAB. CO., LTD.
Priority date: 2008-07-14
Filing date: 2008-07-14
Publication date: 2013-06-19
Anticipated expiration: 2028-07-14
Also published as: JP2010020729A

Description

本発明は、車両走行軌跡観測システム、車両走行軌跡観測方法、およびそのプログラムに関する。

従来の車両走行軌跡観測装置は、複数のビデオカメラを観測区間に配置し、ビデオカメラ各々の撮像結果から得られた走行車両の軌跡を合成することで、長い観測区間の観測を可能としている（例えば、特許文献１参照）。
また、ビデオカメラと、ＧＰＳ（ＧｌｏｂａｌＰｏｓｉｔｉｏｎｉｎｇＳｙｓｔｅｍ；全地球測位システム）、加速度計、３軸ジャイロなどの測定機器と搭載し、走行しながら周囲の走行車両の撮影と、自車両の位置および姿勢角を測定とを行い、周囲の走行車両の観測をしているものもある（例えば、特許文献２参照）。
特開２００３−０８５６８５号公報特開２００７−１４８６１５号公報

しかしながら、特許文献１に示す車両走行軌跡観測装置にあっては、ビデオカメラが撮影した撮影画像から走行車両の位置を算出するための座標変換式を生成するために、撮影画像中の少なくとも４つの標定点の実空間での位置を必要とするが、路面に適当な目標物が無く標定点が決められない、あるいは、標定点の実空間での位置を得るために測量を行おうとしても高速道路など人が立ち入ることができず測量ができず得られない、航空測量を行おうとしてもトンネル内である、標定点の上に橋梁があるなど観測区間の上に構造物があるために標定点の実空間での位置を得るための航空測量が行えず、座標変換式を生成できないという問題がある。

本発明は、このような事情に鑑みてなされたもので、その目的は、撮影画像から走行車両の位置を算出するための座標変換式を、容易に生成できる車両走行軌跡観測装置、車両走行軌跡観測方法、およびそのプログラムを提供することにある。

（１）この発明は上述した課題を解決するためになされたもので、本発明の車両走行軌跡観測システムは、観測区間を撮影するように設置された路側カメラと、前記路側カメラにて撮影された前記観測区間を走行する走行車両の走行軌跡を測定する車両走行軌跡観測装置とを具備する車両走行軌跡観測システムにおいて、自車両の位置を表す情報を検出しながら前記観測区間を走行する計測車両を備え、前記路側カメラは、前記観測区間を走行する前記計測車両を撮影し、前記車両走行軌跡観測装置は、前記路側カメラの位置と、前記路側カメラが撮影した少なくとも４つの撮影画像各々における前記計測車両の座標と、前記計測車両が検出した自車両の位置を表す情報であって、前記少なくとも４つの撮影画像各々の撮影時刻の位置を表す情報とを用いて、前記撮影画像中の座標から、前記観測区間の路面を近似した平面上の位置を算出する座標変換式を生成する座標変換式生成部と、前記座標変換式生成部が生成した座標変換式を用いて、前記路側カメラが前記観測区間を走行する走行車両を撮影した複数の撮影画像各々における前記走行車両の位置を算出し、前記走行車両の走行軌跡を生成する走行軌跡生成部とを具備することを特徴とする。

（２）また、本発明の車両走行軌跡観測システムは、上述の車両走行軌跡観測システムであって、前記計測車両は、自車両の位置を表す情報に加えて、自車両の姿勢角を表す情報の検出と、前記路側カメラの撮影と、前記路側カメラまでの距離の測定とを行いながら走行し、前記車両走行軌跡観測装置は、前記計測車両の位置を表す情報および姿勢角を表す情報と、前記路側カメラを撮影した画像と、前記路側カメラまでの距離とから、前記路側カメラの位置を算出する路側カメラ位置算出部を具備し、前記座標変換式生成部が前記座標変換式を生成する際に用いる前記路側カメラの位置は、前記路側カメラ位置算出部が算出した前記路側カメラの位置であることを特徴とする。

（３）また、本発明の車両走行軌跡観測システムは、上述のいずれかの車両走行軌跡観測システムであって、前記座標変換式生成部は、前記観測区間を複数の領域に分割し、該領域毎に前記座標変換式を生成し、前記走行軌跡生成部は、前記走行車両の位置を算出する際に、前記領域のうち、前記撮影画像中の前記走行車両の座標に応じた領域、または算出済の前記走行車両の位置に応じた領域を選択し、選択した前記領域に対応する前記座標変換式を用いることを特徴とする。

（４）また、本発明の車両走行軌跡観測方法は、観測区間を撮影するように設置された路側カメラと、前記路側カメラにて撮影された前記観測区間を走行する走行車両の挙動を測定する車両走行軌跡観測装置と、自車両の位置を検出しながら走行する計測車両とを具備する車両走行軌跡観測システムにおける車両走行軌跡観測方法であって、前記計測車両が、自車両の位置を検出しながら前記観測区間を走行する第１の過程と、前記路側カメラが、前記観測区間を走行する前記計測車両を撮影する第２の過程と、前記走行車両挙動測定装置が、前記路側カメラの位置と、前記路側カメラが撮影した少なくとも４つの撮影画像各々における前記計測車両の座標と、前記計測車両が検出した自車両の位置を表す情報であって、前記少なくとも４つの撮影画像各々の撮影時刻の位置を表す情報とを用いて、前記該撮影画像中の座標から、前記観測区間の路面を近似した平面上の位置を算出する座標変換式を生成する第３の過程と、前記走行車両挙動測定装置が、前記が第３の過程にて生成した座標変換式を用いて、前記路側カメラが前記観測区間を走行する走行車両を撮影した複数の撮影画像各々における前記走行車両の位置を算出し、前記走行車両の走行軌跡を生成する第４の過程とを備えることを特徴とする。

（５）また、本発明のプログラムは、観測区間を撮影するように設置された路側カメラと、前記路側カメラにて撮影された前記観測区間を走行する走行車両の挙動を測定する車両走行軌跡観測装置と、自車両の位置を検出しながら前記観測区間を走行する計測車両とを具備する車両走行軌跡観測システムにおける前記車両走行軌跡観測装置のコンピュータを、前記路側カメラの位置と、前記路側カメラが撮影した少なくとも４つの撮影画像各々における前記計測車両の座標と、前記計測車両が検出した自車両の位置を表す情報であって、前記少なくとも４つの撮影画像各々の撮影時刻の位置を表す情報とを用いて、前記撮影画像中の座標から、前記観測区間の路面を近似した平面上の位置を算出する座標変換式を生成する座標変換式生成部、前記走行車両挙動測定装置が、前記座標変換式生成部が生成した座標変換式を用いて、前記路側カメラが撮影した複数の撮影画像各々における前記走行車両の位置を算出し、前記走行車両の走行軌跡を生成する走行軌跡生成部として機能させるプログラム。

この発明によれば、計測車両が、自車の位置を測定しながら観測区間を走行し、該走行を路側カメラが撮影することで、路側カメラが撮影した複数の撮影画像各々における計測車両の座標と、計測車両が検出した自車両の位置を表す情報であって、複数の撮影画像各々の撮影時刻の位置を表す情報とに基づき、路側カメラの撮影画像中の座標から位置を算出する座標変換式を、容易に生成することができる。

以下、図面を参照して、本発明の実施の形態について説明する。図１は、この発明の一実施形態による車両走行軌跡観測システムの構成を示す概略ブロック図である。本実施形態による車両走行軌跡観測システムは、車両走行軌跡観測装置１０、計測車両１１、ＧＰＳ（ＧｌｏｂａｌＰｏｓｉｔｉｏｎｉｎｇＳｙｓｔｅｍ；全地球測位システム）時計１２、路側カメラ１３、表示装置１４を具備する。

ＧＰＳ時計１２は、全地球測位システムの人工衛星が送信するＧＰＳ信号に含まれる時刻情報により内部時計を校正し、ビデオカメラのフレームレートである１／２９．９７秒間隔の時刻を表すタイムコード（ＳＭＰＴＥ（ＳｏｃｉｅｔｙｏｆＭｏｔｉｏｎＰｉｃｔｕｒｅａｎｄＴｅｌｅｖｉｓｉｏｎＥｎｇｉｎｅｅｒｓ）タイムコード）と同期信号（ＢｌａｃｋＢｕｒｓｔ信号）を出力する。

路側カメラ１３は、車両走行軌跡観測システムが走行車両の軌跡を観測する区間である観測区間を撮影するように固定した画角で路側に設置されたビデオカメラであり、そのフレームタイミングは、ＧＰＳ時計１２の出力する同期信号に従う。また、路側カメラ１３は、撮影した各フレームに、ＧＰＳ時計１２の出力するタイムコードを記録する。本実施形態では、路側カメラ１３が撮影した撮影画像は、路側カメラ１３に挿入されたビデオテープにいったん記録される。

計測車両１１は、ＧＰＳ時計１１１と自車位置測定器１１２を備え、自車両の位置（実空間における位置）を表す情報を検出しながら、観測区間を走行する。好ましくは、自車両の位置を表す情報に加えて、自車両の姿勢角を表す情報を検出し、固定した画角で設置されたビデオカメラである車載カメラ１１０を備え、走行中に路側カメラ１３を撮影する。本実施形態では、計測車両１１は、自車両の位置を表す情報と姿勢角を表す情報とを検出する自車位置測定器１１２として、進行方向の速度計、３軸の加速度計、ＧＰＳ測位器、３軸回り角速度を測定する３軸ジャイロを備える。

ＧＰＳ測位器は、全地球測位システムの人工衛星が送信するＧＰＳ信号を受信し、ＧＰＳ信号を受信したアンテナの実空間における位置を検出する。なお、ＧＰＳ時計１１１は、ＧＰＳ時計１２と同様に、ＧＰＳ信号に含まれる時刻情報により内部時計を校正し、ビデオカメラのフレームレートである１／２９．９７秒間隔のタイムコードと同期信号を出力する。車載カメラ１１０のフレームタイミングは、ＧＰＳ時計１１１の出力する同期信号に従い、各フレームに、ＧＰＳ時計１１１の出力するタイムコードを記録する。また、車載カメラ１１０は、車載カメラ１１０と路側カメラ１３との距離を検出するレーザー測距装置を備える。自車位置測定器１１２は、ＧＰＳ時計１１１の出力する同期信号に同期させて各値を検出し、検出した各値をＧＰＳ時計１１１の出力するタイムコードとともに記録する。

車両走行軌跡観測装置１０は、車両検出部１００、計測車両位置算出部１０１、標定点選択部１０２、路側カメラ位置算出部１０３、座標変換式生成部１０４、座標変換式記憶部１０５、座標変換式選択部１０６、座標変換部１０７、車両位置記憶部１０８、軌跡表示部１０９を具備する。車両検出部１００は、路側カメラ１３が撮影した撮影画像中の計測車両１１および走行車両を検出し、検出した計測車両１１および走行車両の撮影画像中の座標を、検出したフレームのタイムコードおよび車両を識別する情報とともに出力する。車両検出部１００は、計測車両１１について検出した座標とタイムコードとを標定点選択部１０２に出力する。また車両検出部１００は、走行車両について検出した座標とタイムコードとを座標変換式選択部１０６に出力する。なお、車両検出部１００は、複数の走行車両について観測を行っているときは、座標とタイムコードに加えて、走行車両を識別する情報も座標変換式選択部１０６に出力する。

なお、撮影画像中の計測車両１１および走行車両の検出は、車両の大きさ、輝度、ＲＧＢ値などの外観に関する特徴情報を用いた画像認識処理により、計測車両１１および走行車両を検出し、そのナンバープレートなどの追跡点の座標を検出することで行なってもよいし、撮影画像を表示装置１４などに表示させ、表示された撮影画像中の計測車両１１および走行車両をマウスなどの入力デバイスを用いてオペレータが指定した追跡点の座標を取得することで、検出するようにしてもよい。このとき、車両検出部１００が検出する計測車両１１（追跡点）の座標は、路面を近似する平面を規定する標定点（詳細は後述）の撮影画像上の座標として用いられる。この計測車両１１の追跡点が路面から離れている場合は、追跡点の高さに仮想路面があるとみなし、標定点は、この仮想路面上の点とする。

計測車両位置算出部１０１は、計測車両１１の自車位置測定器１１２が測定した各タイムコードにおける進行方向の速度、３軸の加速度、３軸の角速度、ＧＰＳ測位器が観測した測位結果の実空間での位置とを受けて、計測車両１１の状態遷移方程式と観測方程式とに対して、拡張カルマンスムーザを適用して、各タイムコードにおける計測車両１１の実空間での位置（世界座標系での座標値）と姿勢角とを推定する。なお、本実施形態における計測車両位置算出部１０１が用いる計測車両１１の状態遷移方程式と観測方程式の詳細については、後述する。ここで、実空間での位置を表す世界座標系とは、例えば、南北方向にｘ軸、東西方向にｙ軸をとり、鉛直方向にｚ軸をとった座標系である。

標定点選択部１０２は、観測区間を複数の領域に分割し、各々の領域について路側カメラ１３の撮影画像から実空間の位置を算出するための座標変換式の対象となる路面を規定する標定点を、計測車両位置算出部１０１が算出した計測車両１１の実空間での位置の中から、各々の領域について少なくとも３つ選択する。この選択は、マウスなどの入力デバイスを用いて、オペレータが指定したものを選択するようにしてもよいし、一定時間毎の位置を選択するようにしてもよい。標定点選択部１０２は、計測車両位置算出部１０１が算出した位置の中から少なくとも３つを標定点として選択すると、車両検出部１００から受けた計測車両１１の撮影画像中の座標の中から、選択した標定点のタイムコードと、タイムコードが一致する座標を選択し、標定点（実空間での位置）と対応付けて、座標変換式生成部１０４に出力する。

このとき、座標変換式は、観測区間の路面を平面で近似し、該平面から高さｈを持った追跡点（走行車両のナンバープレートなど）の撮影画像上の座標から実空間における位置を表す座標への座標変換を表す式である。すなわち、撮影画像中の追跡点が、実空間において路面から高さｈを持った位置に存在すると仮定して、座標変換することで、この追跡点の実空間における位置を算出する式である。標定点は、この路面を規定するための点であるので、３つの標定点は、実空間において正三角形に近い形状を成していることが望ましい。計測車両１１を１回だけ走行させた場合、その軌跡は、ほぼ直線であるため、正三角形に近い形状を成す標定点を選択することが困難になる。これを避けるため、計測車両１１を、複数回あるいは複数台、異なる車線などを走行させ、これらの計測結果から標定点を選択するとよい。なお、本実施形態では、観測区間を３つの標定点を頂点とする三角形の領域に分割し、各々の領域について座標変換式を求める。

また、計測車両位置算出部１０１が算出する計測車両の位置は、自車位置測定器１１２のＧＰＳ測位器のアンテナの位置であるため、路側カメラ１３の撮影画像路側カメラ１３の撮影画像から車両検出部１００が検出する追跡点とは一致していないという問題がある。この問題を解決するには、計測車両位置算出部１０１が算出した計測車両１１の姿勢角と、既知の計測車両１１の形状とに基づき、ＧＰＳ測位器のアンテナが設置された位置から、標定点にするのに適した位置（例えば、ナンバープレートなど、路側カメラ１３の撮影画像に含まれ、車両検出部１００が検出可能な点）の値を算出し、これを標定点にすればよい。

路側カメラ位置算出部１０３は、計測車両１１の位置を表す情報および姿勢角を表す情報と、車載カメラ１１０が路側カメラ１３を撮影した画像と、車載カメラ１１０から路側カメラ１３までの距離とから、路側カメラ１３の位置を算出する。より具体的には、路側カメラ位置算出部１０３は、車載カメラ１１０が路側カメラ１３を撮影した画像から、路側カメラ１３を検出し、路側カメラ１３の画像中の座標を取得する。この路側カメラ１３の検出は、車両検出部１００における計測車両１１および走行車両の検出と同様に、路側カメラ１３の大きさ、輝度、ＲＧＢ値などの外観に関する特徴情報を用いた画像認識処理により、路側カメラ１３を検出することで行なってもよいし、撮影画像を表示装置１４などに表示させ、表示された撮影画像中の路側カメラ１３をマウスなどの入力デバイスを用いてオペレータが指定した座標を取得することで、検出するようにしてもよい。

次に、路側カメラ位置算出部１０３は、検出した路側カメラ１３の路側カメラ１３の状態遷移方程式と観測方程式とに対して、拡張カルマンスムーザを適用して、各タイムコードに対応する計測車両１１から路側カメラ１３への相対座標を推定する。そして、路側カメラ位置算出部１０３は、推定した路側カメラ１３への相対座標を、該相対座標のタイムコードに対応する計測車両１１の位置と姿勢角とを用いて、路側カメラ１３の実空間における位置に変換する。なお、本実施形態における路側カメラ位置算出部１０３が用いる路側カメラ１３の状態遷移方程式と観測方程式の詳細については、後述する。

座標変換式生成部１０４は、路側カメラ１３の位置と、路側カメラ１３が撮影した少なくとも４つの撮影画像各々における計測車両１１の座標と、計測車両１１が検出した自車両の位置を表す情報であって、これら少なくとも４つの撮影画像各々の撮影時刻の位置を表す情報とを用いて、すなわち路側カメラ位置算出部１０３が算出した路側カメラ１３の実空間における位置と、標定点選択部１０２が選択した３つの標定点に少なくとも１つを加えた標定点の撮影画像における座標と、該標定点の実空間における位置とを用いて、撮影画像中の座標から、前記観測区間を近似した平面上の位置を算出する座標変換式を生成する。この際に追加する標定点は、標定点選択部１０２が選択した３つの標定点のそれぞれが形成する隣接する他の三角形領域の標定点のうちで、標定点選択部１０２が選択した３つの標定点が成す３次元平面からの高低差が最小のものを選択する。

このとき、路側カメラ１３の位置は、本実施形態では、路側カメラ位置算出部１０３が算出した路側カメラ１３の実空間における位置を用いるが、測量により求めた値を用いてもよい。しかし、高速道路など人の立ち入りに制限がある場所などの理由で測量の実施が困難な場合であっても、本実施形態のように車載カメラ１１０の撮像結果と、計測車両位置算出部１０１が算出した計測車両１１の位置および姿勢角を用いて、路側カメラ位置算出部１０３が路側カメラ１３の位置を算出することで、路側カメラ１３の実空間における位置を得ることができる。

ここで、座標変換式を生成するとは、観測区間を分割した領域各々について、各領域の標定点の情報に基づき座標変換式である式（１）、（１）−１のパラメータα、β、γ、δ、ｂ１〜ｂ８を算出することである。本実施形態における座標変換式の詳細については後述する。図２は、本実施形態における座標変換式のうちの式（１）による、追跡点の路面からの高さの補正方法を説明する図である。図２に示すように、Ｘｃ、Ｙｃは、路側カメラ１３の実空間での位置を表す座標のうちのＸ座標とＹ座標の値であり、路側カメラ位置算出部１０３から得られる。ｈは、路側カメラ１３で撮影した撮影画像から車両検出部１００が走行車両を検出する追跡点Ｔ（例えば、ナンバープレート）の路面からの高さであり、初期値として予め固定値を決めておき、本実施形態では、後述する座標変換部１０７が推定する。なお、高さｈの初期値は、車種毎に予め固定値を決めておき、走行車両の車種に応じて選択した固定値を用いてもよい。図２には表れないが、ｘ、ｙは、路側カメラ１３による撮影画像上の走行車両などの追跡点Ｔの座標であり、車両検出部１００から得られる。また、路側カメラ１３と追跡点Ｔを通る直線と路面との交点である投影点ＰのＸＹ座標（Ｘｐ、Ｙｐ）は、走行車両などの追跡点Ｔの座標路の撮影画像上の座標（ｘ、ｙ）を用いて、２次元射影変換式である式（１）−１により表される。

式（１）のパラメータα、β、γ、δは、標定点選択部１０２が選択した標定点Ｈ１〜３を含む平面、すなわち路面を近似した平面を表す正規化された式（２）を満たす値である。座標変換式生成部１０４は、これらα、β、γ、δの値を、観測区間を分割した領域各々に対して標定点選択部１０２が選択した３つの標定点の実空間での位置を表す座標から求める。

さらに、撮影画像上の座標（ｘ、ｙ）と実空間における位置のＸＹ座標値（Ｘｐ、Ｙｐ）との関係を表す式である式（１）−１を、ｘ、ｙについて解くと式（３）が得られる。座標変換式生成部１０４は、この式（３）に標定点の値（撮影画像上の座標（ｘ、ｙ）と実空間における位置のＸＹ座標値（Ｘｐ、Ｙｐ））を代入すると得られるｂ１〜ｂ８を未知数とした方程式を、最小二乗法などを用いて解き、ｂ１〜ｂ８の値を得る。このとき、未知数は８つあるので、２つの式からなる式（３）を解くには、最少でも４つの標定点の値が必要となる。そこで、本実施形態では、座標変換式生成部１０４は、観測区間を分割した領域各々に対して標定点選択部１０２が選択した３つの標定点に加えて、隣接する領域の標定点のうち、前述の３つの標定点と一致せず、かつ、標定点選択部１０２が選択した３つの標定点が成す３次元平面（式（２）の平面）からのＺ軸方向（鉛直方向）の差が最小の標定点の値を用いて、式（３）を解く。

座標変換式生成部１０４は、観測区間を分割した路面の領域各々の路側カメラ１３による撮影画像中の範囲を表す情報（本実施形態では標定点選択部１０２が選択した３つの標定点の座標）とともに、該領域に対応する座標変換式のパラメータα、β、γ、δ、ｂ１〜ｂ８の値を座標変換式記憶部１０５に格納する。

座標変換式記憶部１０５は、観測区間を分割した領域の範囲を表す情報と、各々の領域に対応する座標変換式を表す情報とを対応付けて記憶する。本実施形態では、観測区間を分割した領域は、３つの標定点を頂点とする三角形であるので、観測区間を分割した領域の範囲を表す情報として、標定点選択部１０２が選択した３つの標定点の路側カメラ１３で撮影した画像中の座標を記憶する。また、座標変換式を表す情報として、式（１）、（１）−１のパラメータα、β、γ、δ、ｂ１〜ｂ８の値を記憶する。

座標変換式選択部１０６は、車両検出部１００から走行車両の撮影画像中の座標とタイムコードを受け、対象の走行車両の位置をいずれの撮影画像についても算出していないときは、該座標に応じた領域を選択し、選択した該領域に対応する座標変換式のパラメータを座標変換式記憶部１０５から読み出し、車両検出部１００から受けた走行車両の座標とタイムコードとともに、読み出した座標変換式のパラメータを座標変換部１０７に出力する。対象の走行車両の位置をいずれかの撮影画像について算出しているときは、座標変換式選択部１０６は、車両位置記憶部１０８が記憶する算出済の走行車両の位置（例えば、１フレーム前の撮影画像から算出した実空間における位置など、最も時間的に近い撮影画像から算出した位置が好ましい）に応じた領域を選択し、選択した該領域に対応する座標変換式のパラメータを座標変換式記憶部１０５から読み出し、車両検出部１００から受けた走行車両の座標とタイムコードとともに、読み出した座標変換式のパラメータを座標変換部１０７に出力する。なお、撮影画像中の走行車両の座標に応じた領域の選択は、該座標を含む領域を選択する、あるいは、該座標に最も重心が近い領域を選択することで行う。また、算出済の走行車両の位置に応じた領域の選択も、該位置を含む領域を選択する、あるいは、該位置に最も重心が近い領域を選択することで行う。

座標変換部１０７は、座標変換式生成部１０４が生成した座標変換式を用いて、路側カメラ１３が観測区間を走行する走行車両を撮影した複数の撮影画像各々に対応する走行車両の位置を算出し、走行車両の走行軌跡を生成する。ここで、座標変換部１０７は、座標変換式生成部１０４が生成した座標変換式として、座標変換式選択部１０６から受けた座標変換式のパラメータを前述の式（１）に代入して作成した座標変換式を用いる。

複数の撮影画像各々における走行車両の位置（実空間における位置のＸＹ座標値）は、このパラメータを用いて作成した座標変換式に、座標変換式選択部１０６から受けた走行車両の座標を代入して演算することでも得られる。しかし、このように座標変換式により算出した結果には観測誤差などが含まれているため、座標変換式による算出結果である走行車両の位置を繋いだ軌跡は、滑らかな線にはならず、実際にはあり得ないジグザグな線になることがある。
そこで、本実施形態では、座標変換部１０７は、上述の座標変換式に誤差分散の項を加えた式と、運動方程式に基づく走行車両の状態遷移方程式（式（３）−１）とに対して、拡張カルマンスムーザを適用して、走行車両の位置を含む状態量を最小二乗推定により逐次更新しながら平滑化する。これにより、座標変換部１０７は、路側カメラ１３が撮影した複数の撮影画像各々に対応する走行車両の位置（実空間における位置のＸＹ座標値と路面からの高さｈ）を算出し、座標変換式選択部１０６から受けたタイムコードとともに、走行車両毎に車両位置記憶部１０８に記憶させる。

状態遷移方程式（式（３）−１）のベクトルＸ（ｋ）、ベクトルＷ、行列Ｆを式（３）−２に示す。状態遷移方程式（式（３）−１）は運動方程式であり、走行車両の位置、速度、加速度、加加速度、追跡点の高さを状態量Ｘ（ｋ）としている。なお、３次項以下は無視している。加加速度とは加速度の１次微分でアクセルペダルまたはブレーキペダルの踏み込み量に相当する。また、追跡点の高さｈを状態遷移方程式に含める事により座標変換において生じた誤差の補正に活用している。なお、ｈの時間変動は無いものとした。以上の状態量ベクトルの誤差分散には時間変動は無いものとして値を設定した。

なお、本実施形態では、座標変換式選択部１０６と座標変換部１０７とで、走行軌跡生成部として機能する。
車両位置記憶部１０８は、走行車両毎に、座標変換部１０７が算出した走行車両の実空間における位置の座標値を、対応するタイムコードとともに記憶する。
軌跡表示部１０９は、車両位置記憶部１０８が記憶する各走行車両の実空間における位置の座標値を読み出し、各走行車両の走行軌跡を表示装置１４に表示させる。

図３は、本実施形態における座標変換式記憶部１０５が記憶内容例を示す図である。図３に示すように、本実施形態における座標変換式記憶部１０５は、観測区間を分割した三角形の全ての領域について、各領域に対応する撮影画像中の領域の３つの頂点、すなわち標定点選択部１０２が選択した３つの標定点である頂点１、頂点２、頂点３の座標（例えば、「ｘ１１、ｙ１１」「ｘ１２、ｙ１２」、「ｘ１３、ｙ１３」）と、座標変換式のパラメータ（例えば、α＝α１、β＝β１、γ＝γ１、δ＝δ１、ｂ１＝ｃ１１、ｂ２＝ｃ１２、ｂ３＝ｃ１３、ｂ４＝ｃ１４、ｂ５＝ｃ１５、ｂ６＝ｃ１６、ｂ７＝ｃ１７、ｂ８＝ｃ１８）とを対応付けて記憶する。（図３のα１〜α３、β１〜β３、γ１〜γ３、δ１〜δ３、ｘ１１〜ｘ３３、ｙ１１〜ｙ３３、ｃ１１〜ｃ３８は、各々数値を表す）

図４は、本実施形態における車両位置記憶部１０８が記憶する走行車両の位置情報の例を示す図である。図４に示す例は、１台の走行車両に対する位置情報である。車両位置記憶部１０８は、車両走行軌跡観測装置１０が観測する走行車両台数分の位置情報を記憶する。図４に示すように、車両位置記憶部１０８は、各タイムコードの時刻（例えば、Ｙ１年Ｍ１月Ｄ１日ｈ１時ｍ１分ｓ１秒ｔ１ミリ秒を表す「Ｙ１／Ｍ１／Ｄ１／ｈ１：ｍ１：ｓ１：ｔ１」）各々に対応付けて、該時刻における走行車両の実空間における位置を表す世界座標系の値（例えば、Ｘ座標「ｘ１」、Ｙ座標「ｙ１」）を記憶する。（図４のＹ１、Ｍ１、Ｄ１、ｈ１〜ｈ３、ｍ１〜ｍ３、ｓ１〜ｓ３、ｔ１〜ｔ３、ｘ１〜ｘ３、ｙ１〜ｙ３は、各々数値を表す）

図５は、本実施形態における車両走行軌跡観測システムの動作を説明するフローチャートである。まず、計測車両１１が、ＧＰＳ時計１１１が受信するＧＰＳ信号により同期された車載カメラ１１０による撮影と自車位置測定器１１２による測定とを行いながら観測区間を走行し、ＧＰＳ時計１２が受信するＧＰＳ信号により同期された路側カメラ１３が、その計測車両１１を撮影する（Ｓ１）。ここでは、このステップＳ１を観測区間にある全ての車線について行う。

次に、路側カメラ１３の撮影結果を車両検出部１００に入力し、自車位置測定器１１２の測定結果を計測車両位置算出部１０１に入力する。自車位置測定器１１２の測定結果を受けて、計測車両位置算出部１０１は、計測車両１１の状態遷移方程式と観測方程式とに対して、拡張カルマンスムーザを適用して、各タイムコードにおける計測車両１１の実空間での位置（世界座標系での座標値）と姿勢角とを算出する（Ｓ２）。

次に、車載カメラ１１０の撮影結果を路側カメラ位置算出部１０３に入力する。路側カメラ位置算出部１０３は、車載カメラ１１０の撮影結果と、ステップＳ２で算出した計測車両１１の位置・姿勢角とから、路側カメラ１３の実空間での位置を算出する（Ｓ３）。次に、標定点選択部１０２は、計測車両位置算出部１０１が算出した各タイムコードにおける計測車両１１の位置、すなわち計測車両１１の走行軌跡から、標定点を選択する（Ｓ４）。このとき、標定点を頂点とする正三角形にできるだけ近い形の領域で観測区間を分割できるように標定点を選択する。

次に、座標変換式生成部１０４は、路側カメラ１３の位置と、観測区間を分割した領域毎に、標定点選択部１０２が選択した該領域の標定点の撮影画像における座標と、実空間における位置とに基づき、撮影画像中の座標から位置を算出する座標変換式を生成し、該座標変換式のパラメータを座標変換式記憶部１０５に記憶させる（Ｓ５）。次に、走行車両を、観測区間を走行させ、路側カメラ１３が、これを撮影する（Ｓ６）。

路側カメラ１３の撮影結果を車両検出部１００に入力する。車両検出部１００は、撮影画像中の走行車両の座標を検出する（Ｓ７）。座標変換式選択部１０６は、車両検出部１００から撮影画像中の走行車両の座標を受けて、該座標を含む領域の座標変換式のパラメータを座標変換式記憶部１０５が記憶するパラメータから選択する。続いて、車両検出部１００は、読み出したパラメータの座標変換式で走行車両の座標を変換し、走行車両の実空間での位置を算出する（Ｓ８）。この走行車両の実空間での位置を繋いで、走行車両の走行軌跡を表示する。

このように、自車位置測定器１１２を備えた計測車両１１が、自車の位置を測定しながら観測区間を走行し、該走行を路側カメラ１３が撮影することで、路側カメラ１３が撮影した複数の撮影画像各々における計測車両１１の座標と、計測車両１１が検出した自車両の位置を表す情報であって、複数の撮影画像各々の撮影時刻の位置を表す情報とに基づき、すなわち計測車両１１の走行軌跡から選択した標定点に基づき、路側カメラ１３の撮影画像中の座標から位置を算出する座標変換式を、容易に生成することができる。
また、計測車両１１に車載カメラ１１０を搭載し、計測車両１１が、自車の位置と姿勢角とを測定しながら観測区間を走行し、該走行中に車載カメラ１１０で路側カメラ１３を撮影しつつ、路側カメラ１３までの距離を測定することで、路側カメラ１３の位置を算出し、路側カメラ１３の撮影画像から位置を算出するための座標変換式を容易に生成することができる。

また、観測区間を分割した領域各々について、各領域の標定点の情報に基づき座標変換式を生成するため、路側カメラ１３のレンズを通して得られた撮影画像の各領域ごとに座標変換式を最適化できることになる。したがって、路側カメラ１３として、写真測量用の高価なレンズを用いたカメラを使用しなくても、路側カメラ１３の汎用レンズの歪みを補正し、所要の座標変換精度を実現することができる。

次に、計測車両位置算出部１０１にて使用する、状態遷移方程式および観測方程式について説明する。状態遷移方程式は、状態変数（ここではベクトルｘ）の時間経過による遷移関係を記述した式であり、一般的に式（４）のように、ある時刻ｔからΔｔ経過したときの状態変数を、時刻ｔにおける状態変数に状態遷移ベクトル関数ｆを施したベクトルと白色雑音ベクトルｗの和で表される。

計測車両位置算出部１０１では、式（５）および式（６）に示すように、実空間の絶対座標系である世界座標系ＸＹＺにおける自車両の位置、速度、加速度、加加速度（加速度の時間微分）、自車両姿勢方向の加速度バイアス、自車両の姿勢角、自車両の姿勢角の角速度、自車両の姿勢角の角加速度、３軸ジャイロの移動座標系ｘｙｚに対する取付不整角、移動座標系ｘｙｚに３軸ジャイロの取付不整角を加えた方向の角速度のバイアスを状態変数とする。これらの状態変数は、それぞれ３要素からなるため、状態変数ベクトルｘは、合計で３０の要素からなるベクトルとなる。なお、ここで、＜ψ、θ、φ＞で表される姿勢角は、ヨー角（φ）、ピッチ角（θ）、ロール角（ψ）の順に回転変換して得られるＮａｖｉｇａｔｉｏｎ角と呼ばれる姿勢角である（ベクトルの要素の並び順と回転変換の順は逆である）。

本実施形態では、式（５）、（６）に示した状態変数に関する状態遷移方程式の状態遷移ベクトル関数ｆを式（７）、白色雑音ベクトルｗを式（８）、（９）とする。

状態遷移ベクトル関数ｆ（式（７））を参照すると分かるように、状態遷移方程式は線形方程式となっているので、状態遷移方程式は式（１０）のように行列表現できる。このとき、状態遷移行列Ｆは式（１１）のようになる。

観測方程式は、観測機器により観測した値と、状態変数の関係を表す式であり、一般的に、式（１２）のように、観測変数ベクトルｘ^ｍを、状態変数ベクトルｘに観測ベクトル関数ｈを施したベクトルと白色雑音ベクトルｅの和で表される。

計測車両位置算出部１０１は、自車位置測定器１１２であるＧＰＳ時計、速度計、３軸ジャイロの出力である世界座標系ＸＹＺにおける自車の位置、進行方向の速度、３軸ジャイロ取付方向に対する加速度と角速度が、観測変数ベクトルｘ^ｍの要素となるので、観測変数ベクトルｘ^ｍは、式（１３）のように表され、観測ベクトル関数ｈは式（１４）となる。

観測ベクトル関数ｈ（式（１４））は、非線形であるため、状態遷移ベクトル関数ｆのように行列表現することはできないが、観測ベクトル関数ｈを状態変数ベクトルｘで偏微分した、観測行列Ｈを用いて、近似的に表すことができる。このとき、観測方程式は、式（１５）のように表され、観測行列Ｈは、式（１６）で表される。

また、計測車両位置算出部１０１における観測方程式の白色雑音ベクトルｅは、式（２５）である。

次に、路側カメラ位置算出部１０３における状態遷移方程式は、計測車両位置算出部１０１と同様に、式（４）で表される。本実施形態では、式（２６）に示すように、計測車両１１に固定された移動座標系ｘｙｚにおける車載カメラ１１０からの路側カメラ１３の相対位置、加速度、加加速度のｘｙ成分、車載カメラ１１０の取付高さと路側カメラ１３の取付高さの差、車載カメラ１１０の取付角と焦点距離、路側カメラ１３までの距離を測定するレーザー測距装置の測距点と車載カメラ１１０によるトラッキング点との較差を状態変数ベクトルｘ_ｉとする。また、状態遷移ベクトル関数ｆは、式（２７）とする。

状態遷移ベクトル関数ｆ（式（２６））を参照すると分かるように、路側カメラ１０３位置算出部１０３においても、状態遷移方程式は線形方程式となっているので、状態遷移方程式は式（１０）のように行列表現できる。このとき、状態遷移行列Ｆは式（２８）のようになる。

また、路側カメラ１０３位置算出部１０３における状態遷移方程式の白色雑音ベクトルｗは、式（２９）とする。

次に、路側カメラ位置算出部１０３における観測方程式は、計測車両位置算出部１０１と同様に、式（３０）で表される。本実施形態における路側カメラ位置算出部１０３では、車載カメラ１１０の出力画像より抽出した路側カメラ１３の座標および車載カメラ１１０のレーザー測距装置の出力である路側カメラ１３までの距離が、観測変数ベクトルｙ_ｉの要素となるので、観測変数ベクトルｙ_ｉは、式（３１）のように表され、観測ベクトル関数ｈは式（３２）となる。

計測車両位置算出部１０１と同様に、路側カメラ位置算出部１０３においても、観測ベクトル関数ｈ式（３２）は、非線形であるため、状態遷移ベクトル関数ｆのように行列表現することはできないが、観測ベクトル関数ｈを状態変数ベクトルｘ_ｉで偏微分した、観測行列Ｈを用いて、近似的に表すことができる。このとき、観測方程式は、式（３３）のように表され、観測行列Ｈは、式（３４）で表される。

また、路側カメラ位置算出部１０３における観測方程式の白色雑音ベクトルｅは、式（３７）である。

次に、本実施形態における座標変換式（式（１））について説明する。一般に、平面（Ｘｐ、Ｙｐ）から平面（ｘ、ｙ）への２次元射影変換は、式（３８）で表される。これが、式（１）−１である。

図６と図７とは、本実施形態における座標変換式である式（１）を説明する図である。図６は、路側カメラ１３と追跡点Ｐとを結ぶ線に垂直かつ、路面Ｒに平行な方向から見たときの、路側カメラ１３と走行車両と路面Ｒとの関係を示す図である。図７は、世界座標系のＸ軸、Ｙ軸を水平面にとり、Ｚ軸方向から見たときの路側カメラ１３と走行車両と投影点Ｐとの関係を示す図である。これらの図において、路側カメラ１３から投影点Ｐまでの距離をＬｐ、路側カメラ１３から追跡点Ｔまでの距離をＬ、路面Ｒから路側カメラ１３までの高さをｈｃ、路面Ｒから追跡点Ｔまでの高さをｈとする。

図６に示すように、追跡点Ｔの路面Ｒへの垂線と路面Ｒとの交点、投影点Ｐ、追跡点Ｔの３点を頂点とする三角形と、路側カメラ１３の路面Ｒへの垂線と路面Ｒとの交点、投影点Ｐ、路側カメラ１３の３点を頂点とする三角形とは相似である。従って、式（３９）が得られる。式（３９）を変形すると式（４０）が得られる。
ｈｃ：ｈ＝Ｌｐ：Ｌｐ−Ｌ …（３９）
ｈｃ：ｈｃ−ｈ＝Ｌｐ：Ｌ …（４０）

また、路面Ｒを表す正規化された式は、αＸ＋βＹ＋δＺ＋γ＝０であることから、路側カメラ１３（Ｘｃ、Ｙｃ、Ｚｃ）から路面Ｒまでの距離ｈｃは、式（４１）で表される。
ｈｃ＝αＸｃ＋βＹｃ＋δＺｃ＋γ …（４１）
また、図７に示すように、路側カメラ１３の水平面への射影と追跡点Ｔの水平面への射影と点（Ｘ，Ｙｃ）とからなる三角形と、路側カメラ１３の水平面への射影と投影点Ｐと点（Ｘｐ，Ｙｃ）とからなる三角形とは相似である。従って、Ｘｐ−Ｘｃ：Ｘ−Ｘｃ＝Ｌ’ｐ：Ｌ’である。また、Ｌ’とＬ’ｐは、図６のＬとＬｐの水平面への射影であることから、Ｌ’ｐ：Ｌ’＝Ｌｐ：Ｌである。従って、式（４２）が得られる。
Ｘｐ−Ｘｃ：Ｘ−Ｘｃ＝Ｌｐ：Ｌ …（４２）
式（４０）と式（４２）とから、Ｘｐ−Ｘｃ：Ｘ−Ｘｃ＝ｈｃ：ｈｃ−ｈであり、この式と、式（４１）とから、式（４３）が得られる。この式（４３）を変形すると、式（１）のＸｐに関する式と一致する。また、式（１）のＹｐに関する式も同様にして得られる。このように、本実施形態における座標変換式である式（１）は、実空間における位置が座標（Ｘ、Ｙ）で路面からの高さｈの追跡点に関する撮影画像上の座標（ｘ、ｙ）を求める式となっている。

なお、本実施形態において路側カメラ１３は一つであるとして説明したが、複数の路側カメラ１３を備え、各路側カメラ１３について観測区間を分割した領域毎に座標変換式を生成し、これらの座標変換式を用いて算出した走行車両の走行軌跡を合成するようにしてもよい。
また、本実施形態において、車載カメラ１１０は、レーザー測距装置を備え、路側カメラ１３までの距離を測定するとして説明したが、レーザー測距装置を備えなくてもよい。この場合、路側カメラ位置算出部１０３における状態変数ベクトル（式（２６））における要素ｓ_ｉ’（レーザー測距装置の測距点とトラッキング点との較差）、状態遷移行列Ｆ（式（２８））の最終列、白色雑音ベクトルｗ（式（２９））の１４番目（最後）の要素ｗ_ｓ、観測変数ベクトル（式（３１））の３番目の要素ｇ_ｉ（ｔ）、観測行列Ｈ（式（３４））の３行目、白色雑音ベクトルｅ（式（３７））の３番目の要素が不要となる。

なお、座標変換式記憶部１０５および車両位置記憶部１０８は、ハードディスク装置や光磁気ディスク装置、フラッシュメモリ等の不揮発性のメモリや、ＲＡＭ（ＲａｎｄｏｍＡｃｃｅｓｓＭｅｍｏｒｙ）のような揮発性のメモリ、あるいはこれらの組み合わせにより構成されるものとする。

また、この車両走行軌跡観測装置１０には、表示装置１４に加えて、周辺機器として入力装置等（図示せず）が接続されるものとする。ここで、入力装置とはキーボード、マウス等の入力デバイスのことをいう。表示装置１４とはＣＲＴ（ＣａｔｈｏｄｅＲａｙＴｕｂｅ）や液晶表示装置等のことをいう。

また、図１における車両検出部１００、計測車両位置算出部１０１、標定点選択部１０２、路側カメラ位置算出部１０３、座標変換式生成部１０４、座標変換式選択部１０６、座標変換部１０７、軌跡表示部１０８の機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより各部の処理を行ってもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。

また、「コンピュータシステム」は、ＷＷＷシステムを利用している場合であれば、ホームページ提供環境（あるいは表示環境）も含むものとする。
また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含むものとする。また上記プログラムは、前述した機能の一部を実現するためのものであっても良く、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであっても良い。

以上、この発明の実施形態を図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計変更等も含まれる。

本発明は、二輪自動車および自動車の走行軌跡を観測する車両走行軌跡観測システムに用いて好適であるが、これに限定されない。

この発明の一実施形態による車両走行軌跡観測システムの構成を示すブロック図である。同実施形態における座標変換式のうちの式（１）による、追跡点の路面からの高さの補正方法を説明する図である。同実施形態における座標変換式記憶部１０５が記憶内容例を示す図である。同実施形態における車両位置記憶部１０８が記憶する走行車両の位置情報の例を示す図である。同実施形態における車両走行軌跡観測システムの動作を説明するフローチャートである。同実施形態における座標変換式である式（１）を説明する図である。同実施形態における座標変換式である式（１）を説明する図である。

符号の説明

１０…車両走行軌跡観測装置
１００…車両検出部
１０１…計測車両位置算出部
１０２…標定点選択部
１０３…路側カメラ位置算出部
１０４…座標変換式生成部
１０５…座標変換式記憶部
１０６…座標変換式選択部
１０７…座標変換部
１０８…車両位置記憶部
１０９…軌跡表示部
１１…計測車両
１１０…車載カメラ
１１１…ＧＰＳ時計
１１２…自車位置測定器
１２…ＧＰＳ時計
１３…路側カメラ
１４…表示装置

Claims

観測区間を撮影するように設置された路側カメラと、前記路側カメラにて撮影された前記観測区間を走行する走行車両の走行軌跡を測定する車両走行軌跡観測装置とを具備する車両走行軌跡観測システムにおいて、
自車両の位置を表す情報を検出しながら前記観測区間を走行する計測車両を備え、
前記路側カメラは、前記観測区間を走行する前記計測車両を撮影し、
前記車両走行軌跡観測装置は、
前記路側カメラの位置と、前記路側カメラが撮影した少なくとも４つの撮影画像各々における前記計測車両の座標と、前記計測車両が検出した自車両の位置を表す情報であって、前記少なくとも４つの撮影画像各々の撮影時刻の位置を表す情報とを用いて、前記撮影画像中の座標から、前記観測区間の路面を近似した平面上の位置を算出する座標変換式を生成する座標変換式生成部と、
前記座標変換式生成部が生成した座標変換式を用いて、前記路側カメラが前記観測区間を走行する走行車両を撮影した複数の撮影画像各々に対応する前記走行車両の位置を算出し、前記走行車両の走行軌跡を生成する走行軌跡生成部と
を具備し、
前記計測車両は、自車両の位置を表す情報に加えて、自車両の姿勢角を表す情報の検出と、前記路側カメラの撮影と、前記路側カメラまでの距離の測定とを行いながら走行し、
前記車両走行軌跡観測装置は、前記計測車両の位置を表す情報および姿勢角を表す情報と、前記路側カメラを撮影した画像と、前記路側カメラまでの距離とから、前記路側カメラの位置を算出する路側カメラ位置算出部を具備し、
前記座標変換式生成部が用いる前記路側カメラの位置は、前記路側カメラ位置算出部が算出した前記路側カメラの位置であること
を特徴とする車両走行軌跡観測システム。
観測区間を撮影するように設置された路側カメラと、前記路側カメラにて撮影された前記観測区間を走行する走行車両の挙動を測定する車両走行軌跡観測装置と、自車両の位置を検出しながら走行する計測車両とを具備する車両走行軌跡観測システムにおける車両走行軌跡観測方法であって、
前記計測車両が、自車両の位置と、自車両の姿勢角を表す情報の検出と、前記路側カメラの撮影と、前記路側カメラまでの距離の測定を検出しながら前記観測区間を走行する第１の過程と、
前記路側カメラが、前記観測区間を走行する前記計測車両を撮影する第２の過程と、
前記車両走行軌跡観測装置が、前記計測車両の位置を表す情報および姿勢角を表す情報と、前記路側カメラを撮影した画像と、前記路側カメラまでの距離とから、前記路側カメラの位置を算出する第３の過程と、
前記車両走行軌跡観測装置が、前記路側カメラの位置と、前記路側カメラが撮影した少なくとも４つの撮影画像各々における前記計測車両の座標と、前記計測車両が検出した自車両の位置を表す情報であって、前記少なくとも４つの撮影画像各々の撮影時刻の位置を表す情報とを用いて、前記該撮影画像中の座標から、前記観測区間の路面を近似した平面上の位置を算出する座標変換式を生成する第４の過程と、
前記車両走行軌跡観測装置が、前記第４の過程にて生成した座標変換式を用いて、前記路側カメラが前記観測区間を走行する走行車両を撮影した複数の撮影画像各々における前記走行車両の位置を算出し、前記走行車両の走行軌跡を測定する第５の過程と
を備えることを特徴とする車両走行軌跡観測方法。
観測区間を撮影するように設置された路側カメラと、前記路側カメラにて撮影された前記観測区間を走行する走行車両の挙動を測定する車両走行軌跡観測装置と、自車両の位置を検出しながら前記観測区間を走行する計測車両とを具備する車両走行軌跡観測システムにおける前記車両走行軌跡観測装置のコンピュータを、
前記計測車両の位置を表す情報および姿勢角を表す情報と、前記路側カメラを撮影した画像と、前記路側カメラまでの距離とから、前記路側カメラの位置を算出する路側カメラ位置算出部、
前記路側カメラの位置と、前記路側カメラが撮影した少なくとも４つの撮影画像各々における前記計測車両の座標と、前記計測車両が検出した自車両の位置を表す情報であって、前記少なくとも４つの撮影画像各々の撮影時刻の位置を表す情報とを用いて、前記撮影画像中の座標から、前記観測区間の路面を近似した平面上の位置を算出する座標変換式を生成する座標変換式生成部、
前記走行車両挙動測定装置が、前記座標変換式生成部が生成した座標変換式を用いて、前記路側カメラが前記観測区間を走行する走行車両を撮影した複数の撮影画像各々における前記走行車両の位置を算出し、前記走行車両の走行軌跡を生成する走行軌跡生成部
として機能させるプログラム。