"GPS衛星から送信された擬似距離及びドップラー周波数を含むGPS情報を受信するGPS受信機から該GPS情報をエポック毎に取得すると共に、自装置が搭載された移動体の移動量を含むINS情報を検出するINS装置から該INS情報をエポック毎に取得する取得手段と、前記取得手段により取得された複数エポック分のGPS情報とINS情報とを統合して、前記移動体の位置を測位する測位手段と、前記INS装置の観測誤差を、各エポックにおける前記移動体から見た前記GPS衛星の方向を用いて前記GPS受信機の観測誤差に換算して該GPS受信機の観測誤差と統合し、統合した観測誤差に基づいて、前記測位手段による測位の信頼度を計算する計算手段と、を含む測位装置。前記測位手段は、前記取得手段により取得された複数エポック分のINS情報から前記移動体の移動軌跡を算出し、該移動軌跡を前記擬似距離を用いて位置あわせすることにより、前記移動体の位置を測位する請求項1記載の測位装置。前記測位手段により測位された移動体の位置との残差が予め定めた閾値以上となる擬似距離を含むGPS情報、または該GPS情報を送信したGPS衛星から送信された全てのGPS情報を、前記取得手段により取得したGPS情報から排除して、前記測位手段による測位、及び前記計算手段による信頼度の計算が再度行われるように制御する制御手段を含む請求項1または請求項2記載の測位装置。コンピュータを、GPS衛星から送信された擬似距離及びドップラー周波数を含むGPS情報を受信するGPS受信機から該GPS情報をエポック毎に取得すると共に、自装置が搭載された移動体の移動量を含むINS情報を検出するINS装置から該INS情報をエポック毎に取得する取得手段、前記取得手段により取得された複数エポック分のGPS情報とINS情報とを統合して、前記移動体の位置を測位する測位手段、及び前記INS装置の観測誤差を、各エポックにおける前記移動体から見た前記GPS衛星の方向を用いて前記GPS受信機の観測誤差に換算して該GPS受信機の観測誤差と統合し、統合した観測誤差に基づいて、前記測位手段による測位の信頼度を計算する計算手段として機能させるための測位プログラム。しかし、特許文献1に記載の技術では、機器の誤差特性から計算される信頼度の範囲を単純な円で表現しているが、その具体的な計算方法は明記されていない。また、特許文献2においても、具体的な誤差範囲の計算方法は記載されていない。また、特許文献3に記載の技術では、GPS情報とINS情報とをカルマンフィルタを用いて統合し、カルマンフィルタの共分散を信頼度としている。しかし、カルマンフィルタを適用した場合において、観測値に非ガウスノイズが存在する場合、信頼度の精度が低くなる可能性がある。また、特許文献1~3に記載の技術はいずれも、複数エポックの観測値を用いる方式ではないため、複数エポックのGPS情報とINS情報とを統合する場合の信頼度については記載されていない。本発明は、上記の問題点を解決するためになされたもので、GPS情報とINS情報とを統合して測位する場合に、測位解の信頼度を精度良く計算することができる測位装置及びプログラムを提供することを目的とする。上記目的を達成するために、本発明の測位装置は、GPS衛星から送信された擬似距離及びドップラー周波数を含むGPS情報を受信するGPS受信機から該GPS情報をエポック毎に取得すると共に、自装置が搭載された移動体の移動量を含むINS情報を検出するINS装置から該INS情報をエポック毎に取得する取得手段と、前記取得手段により取得された複数エポック分のGPS情報とINS情報とを統合して、前記移動体の位置を測位する測位手段と、前記INS装置の観測誤差を、各エポックにおける前記移動体から見た前記GPS衛星の方向を用いて前記GPS受信機の観測誤差に換算して該GPS受信機の観測誤差と統合し、統合した観測誤差に基づいて、前記測位手段による測位の信頼度を計算する計算手段と、を含んで構成されている。本発明の測位装置によれば、取得手段が、GPS衛星から送信された擬似距離及びドップラー周波数を含むGPS情報を受信するGPS受信機からGPS情報をエポック毎に取得すると共に、自装置が搭載された移動体の移動量を含むINS情報を検出するINS装置からINS情報をエポック毎に取得する。そして、測位手段が、取得手段により取得された複数エポック分のGPS情報とINS情報とを統合して、移動体の位置を測位し、計算手段が、INS装置の観測誤差を、各エポックにおける移動体から見たGPS衛星の方向を用いてGPS受信機の観測誤差に換算してGPS受信機の観測誤差と統合し、統合した観測誤差に基づいて、測位手段による測位の信頼度を計算する。このように、複数エポック分のGPS情報とINS情報とを統合して、移動体の位置を測位する際に、INS装置の観測誤差をGPS受信機の観測誤差に換算して、INS装置の観測誤差とGPS受信機の観測誤差とを統合するため、INS装置の観測誤差も加味して、測位解の信頼度を精度良く計算することができる。また、前記測位手段は、前記取得手段により取得された複数エポック分のINS情報から前記移動体の移動軌跡を算出し、該移動軌跡を前記擬似距離を用いて位置あわせすることにより、前記移動体の位置を測位することができる。このような測位は、各エポックの位置ではなく、移動軌跡の位置のみを計算するだけでよいため、演算量を削減することができる。GPS情報のみの観測方程式においても、上記のように、INS装置の観測誤差をGPS受信機の観測誤差に換算して、INS装置の観測誤差とGPS受信機の観測誤差とを統合するため、INS装置の観測誤差も加味して、測位解の信頼度を精度良く計算することができる。また、本発明の測位装置は、前記測位手段により測位された移動体の位置との残差が予め定めた閾値以上となる擬似距離を含むGPS情報、または該GPS情報を送信したGPS衛星から送信された全てのGPS情報を、前記取得手段により取得したGPS情報から排除して、前記測位手段による測位、及び前記計算手段による信頼度の計算が再度行われるように制御する制御手段を含んで構成することができる。これにより、マルチパス等の影響を受けたつじつまの合わないGPS情報を排除することができるため、測位精度がより向上する。また、本発明の測位プログラムは、コンピュータを、GPS衛星から送信された擬似距離及びドップラー周波数を含むGPS情報を受信するGPS受信機から該GPS情報をエポック毎に取得すると共に、自装置が搭載された移動体の移動量を含むINS情報を検出するINS装置から該INS情報をエポック毎に取得する取得手段、前記取得手段により取得された複数エポック分のGPS情報とINS情報とを統合して、前記移動体の位置を測位する測位手段、及び前記INS装置の観測誤差を、各エポックにおける前記移動体から見た前記GPS衛星の方向を用いて前記GPS受信機の観測誤差に換算して該GPS受信機の観測誤差と統合し、統合した観測誤差に基づいて、前記測位手段による測位の信頼度を計算する計算手段として機能させるためのプログラムである。なお、本発明のプログラムを記憶する記憶媒体は、特に限定されず、ハードディスクであってもよいし、ROMであってもよい。また、CD-ROMやDVDディスク、光磁気ディスクやICカードであってもよい。更にまた、該プログラムを、ネットワークに接続されたサーバ等からダウンロードするようにしてもよい。以上説明したように、本発明の測位装置及びプログラムによれば、複数エポック分のGPS情報とINS情報とを統合して、移動体の位置を測位する際に、INS装置の観測誤差をGPS受信機の観測誤差に換算して、INS装置の観測誤差とGPS受信機の観測誤差とを統合するため、INS装置の観測誤差も加味して、測位解の信頼度を精度良く計算することができる、という効果が得られる。以下、図面を参照して本発明の実施の形態を詳細に説明する。本実施の形態では、車両に搭載され、自車両の位置を測位する測位装置に本発明を適用した例について説明する。図1に示すように、本実施の形態の測位装置10は、GPS衛星から送信されたGPS情報を受信するGPS受信機12と、自車両の移動量を含むINS情報を取得するINS装置14と、GPS情報及びINS情報に基づいて測位解及びその信頼度を算出するコンピュータ16とを含んで構成されている。GPS受信機12で受信されるGPS情報は、擬似距離、ドップラー周波数、及びエフェメリスを含む。また、INS装置14は、速度センサ、ヨーレイトセンサ、及びジャイロセンサを備え、各センサにより検出された検出値に基づいて、自車両の移動量を算出する。コンピュータ16は、CPU、後述する測位処理ルーチンを実現するためのプログラムを記憶したROM、データを一時的に記憶するRAM、及びHDD等の記憶装置を含んで構成されている。このコンピュータ16を以下で説明する測位処理ルーチンに従って機能ブロックで表すと、図1に示すように、Mエポック分のGPS情報及びINS情報に基づいて自車両の走行軌跡を算出する走行軌跡算出部20と、エポック毎にGPS受信機12で受信したGPS情報及びINS装置14で取得したINS情報を記憶するMエポックデータ記憶部22と、算出した走行軌跡に基づいて、Mエポック分の自車両の位置の初期値を算出する初期値算出部24と、INS装置14の観測誤差以下、INS誤差というを擬似距離誤差に換算すると共に、初期値に基づくGPS情報の観測方程式を導出する誤差換算方程式導出部26と、観測方程式から測位解を計算する測位解計算部28と、残差に基づいて整合性のないGPS衛星を選択する衛星選択部30とを含んだ構成で表すことができる。ここで、本実施の形態の原理について説明する。GPS情報とINS情報とを統合した測位の方式として、図2に示すようなタイトカップリング方式がある。この方式では、GPS受信機及びINS装置の観測を元に、Mエポックの拘束条件を求める。すなわち、下記1式に示すように、擬似距離ρGPS情報と走行軌跡ΔxINS情報との連立方程式を解くことにより、各エポックでの自車両の位置x1,x2,,xMを求める。タイトカップリング方式では、エポック数Mが増加すると、1式内の行列Gのサイズも増大し、演算コストが大きくなるという問題がある。また、GPS情報とINS情報とを統合した測位の別の方式として、図3に示すようなバンドルアジャストメント方式がある。この方式では、予め方位を定めて、INS情報から走行軌跡を算出しておき、その走行軌跡をGPS情報で位置合わせすることにより、走行軌跡の始点の位置を求める。すなわち、下記2式に示すように、擬似距離ρGPS情報と走行軌跡の位置あわせ量dxとの連立方程式を解くことにより、各エポックでの自車両の位置x1,x2,,xMを求める。バンドルアジャストメント方式では、INS情報から算出した走行軌跡が正しいと仮定しているため、2式に示すようにパラメータを減らして演算量を削減できる。しかし、方程式にINS情報の変数が含まれていないため、INS誤差を信頼度に反映できないという問題がある。そこで、本実施の形態では、図4に示すように、バンドルアジャストメント方式において、INS誤差を擬似距離誤差に換算することで、演算量を増加させることなく、INS誤差も含んだ信頼度を算出する。以下、iをGPS衛星番号1~n、jをエポック番号1~Mとし、推定する状態[xM,yM,zM,CbM]Mエポック目の位置及びクロックバイアスとして、各部について詳述する。走行軌跡算出部20は、INS装置14から自車両の移動量を取得し、GPS受信機12からドップラー周波数を取得する。そして、自車両の走行軌跡として、1エポック目からjエポック目までの移動量ベクトルΔxj、1エポック目からjエポック目までの時間変化Δtj、クロックドリフトCd、及び方位角θの値を求める。これらの値は、例えば、参考文献GPSドップラーと慣性センサの統合による車両軌跡推定手法第43回情報処理学会ITS研究会の手法を用いることができる。初期値算出部24は、1エポック目の自車両の位置の初期値を定め、走行軌跡算出部20で算出された各値に基づいて、下記3式により、各エポックでの自車両の位置の初期値を算出し、下記4式により、各エポックでのINS誤差を算出し、下記5式により、各エポックのクロックバイアスを算出し、下記6式により、各エポックでのクロックバイアス誤差を算出する。ここで、Rは方位角θの回転行列、σvはINS装置14で検出される速度のばらつき、及びσcdはクロックドリフトのばらつきである。なお、σv及びσcdは固定値を用いることができる。誤差換算方程式導出部26は、下記7式により、INS誤差を擬似距離誤差に換算する。ここで、σρは擬似距離誤差、gi,jはGPS受信機12からjエポック目のGPS衛星iへの単位方向ベクトルである。従って、7式は、INS誤差を、GPS衛星の方向に投影して、GPS受信機12の観測誤差以下、GPS誤差という。ここでは、擬似距離誤差σρ及びクロックバイアス誤差σcbに加算することで、GPS誤差とINS誤差とを統合した観測誤差を計算していることを表している。なお、7式では、擬似距離誤差と換算したINS誤差との二乗和の平方根により、統合した観測誤差を計算する例を示しているが、各誤差の絶対値の和などを用いて計算してもよく、観測誤差の統合方法は限定されない。また、誤差換算方程式導出部26は、下記8式及び9式に示す初期値算出部24で算出された初期値と真の値との関係に基づいて、下記10式及び11式の観測方程式を導出する。なお、10式は、1つのGPS衛星に関する観測方程式であり、11式は、Mエポック分の全GPS衛星に関する観測方程式である。また、このときの観測値の共分散行列は、下記12式となる。測位解計算部28は、誤差換算方程式導出部26で導出された観測方程式を解いて、下記13式により自車両の位置xの測位解を計算する。また、その共分散Qxを下記14により計算する。このQxが測位解の信頼度である。衛星選択部30は、計算された測位解とのつじつまが合わないGPS衛星からのGPS情報を排除する。つじつまが合わないとは、図5に示すように、各エポックでの測位解が示す位置とそのエポックでの衛星位置との距離と、そのエポックでの擬似距離との差が大きいことをいう。具体的には、11式の観測方程式で推定される状態と、実際の観測値との残差を、下記15式により計算する。残差rが所定の閾値以上となるGPS衛星をつじつまの合わないGPS衛星として排除し、その衛星番号を誤差換算方程式導出部26に出力する。誤差換算方程式導出部26では、衛星選択部30から出力された衛星番号のGPS衛星から受信したGPS情報を排除して、再度、観測方程式を導出し、測位解計算部28で、再度、測位解を計算する。次に、図6を参照して、本実施の形態の測位装置10において実行される測位処理ルーチンについて説明する。ステップ100で、エポック毎にGPS受信機12で受信したGPS情報及びINS装置14で取得したINS情報を取得して、Mエポックデータ記憶部22記憶する。次に、ステップ102で、INS装置14から自車両の移動量を取得し、GPS受信機12からドップラー周波数を取得する。そして、自車両の走行軌跡として、1エポック目からjエポック目までの移動量ベクトルΔxj、1エポック目からjエポック目までの時間変化Δtj、クロックドリフトCd、及び方位角θの値を求める。次に、ステップ104で、1エポック目の自車両の位置の初期値を定め、上記ステップ102で求めた各値に基づいて、3式により、各エポックでの自車両の位置の初期値を算出し、4式により、各エポックでのINS装置14の観測誤差を算出し、5式により、各エポックのクロックバイアスを算出し、6式により、各エポックでのクロックバイアス誤差を算出する。次に、ステップ106で、7式により、INS誤差を擬似距離誤差に換算する。次に、ステップ108で、上記ステップ104で求めた初期値と真の値との関係に基づいて、11式の観測方程式を導出する。次に、ステップ110で、上記ステップ108で導出された観測方程式を解いて、13式により自車両の位置xの測位解を計算し、14式により、その信頼度共分散Qxを計算する。次に、ステップ112で、15式により残差rを計算し、残差rが所定の閾値以上となる擬似距離を送信したGPS衛星、すなわち、つじつまの合わないGPS衛星の有無を判定する。次に、ステップ114で、上記ステップ112での判定結果に基づいて、上記ステップ110で計算された測位解に整合性があるか否かを判定する。上記ステップ112で、つじつまの合わないGPS衛星が存在すると判定された場合には、本ステップで整合性なしと判定されて、ステップ116へ移行する。ステップ116では、取得したMエポック分のGPS情報のうち、上記ステップ112でつじつまが合わない判定されたGPS衛星から受信したGPS情報を排除して、ステップ106へ戻る。これにより、つじつまの合わないGPS衛星からのGPS情報が排除された状態で、再度、観測方程式が導出され、測位解が計算される。上記ステップ114で、測位解に整合性あり、すなわち、つじつまの合わないGPS衛星が存在しないと判定された場合には、ステップ118へ移行し、上記ステップ112で計算された測位解及び信頼度を出力して、測位処理ルーチンを終了する。以上説明したように、本実施の形態の測位装置によれば、複数エポック分のGPS情報及びINS情報を統合して測位する際に、INS誤差をGPS衛星の方向に投影してGPS誤差に換算することで、測位解の信頼度として、INS誤差も加味した精度の高い信頼度を得ることができる。また、GPS情報とINS情報との統合にバンドルアジャストメント方式を用いることで、方程式のパラメータを少なくすることができ、演算量を削減することができ、このような方式でも、誤差換算を行うことにより、INS誤差も加味した信頼度を得ることができる。なお、上記実施の形態では、衛星選択部でつじつまの合わないGPS衛星を選択し、そのGPS衛星からのGPS情報を排除する場合について説明したが、そのGPS衛星からの全てのGPS情報を排除するのではなく、つじつまの合わないGPS情報のみを排除するようにしてもよい。"