(57)【特許請求の範囲】
【請求項１】
外界の三次元形状を取り入れて自己位置を同定する自己位置同定方法であって、
新たな計測位置において、三次元形状上の座標値をコンピュータに入力するデータ入力ステップと、
前記三次元形状の存在する空間領域を、境界表面が互いに直交する直方体からなる複数のボクセルに分割し、各ボクセル位置を記憶する環境モデルを構築するモデル構築ステップと、
前記座標値に対応するボクセルの内部に代表点とその誤差分布を設定し記憶するマッチングステップとを実施し、
先の計測位置おける前記三次元形状のデータが存在しない場合には、新たな計測位置を自己位置と同定し、
先の計測位置おける前記三次元形状のデータが存在する場合には、
先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、近接する誤差分布間の距離の総和が最小になるように位置合わせする精密合わせステップと、
前記精密合わせステップにおける回転量及び並進量から自己位置を同定する自己位置同定ステップとを実施し、
更に前記自己位置を出力装置に出力する出力ステップを実施する、ことを特徴とする自己位置同定方法。
【請求項２】
前記マッチングステップにおいて、ボクセルの内部に代表点とその誤差分布に加え、ボクセル内に物体の存在確率を表す確率値を設定し、記憶する、ことを特徴とする請求項１に記載の自己位置同定方法。
【請求項３】
前記精密合わせステップの前に、先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、計測データおよび誤差分布と近接する代表点を有するボクセル間との距離の総和が最小になるように、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、代表点を有するボクセル間の距離の総和が最小になるように、位置合わせする粗合わせステップを有する、ことを特徴とする請求項１に記載の自己位置同定方法。
【請求項４】
前記精密合わせステップの前に、先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、計測データおよび誤差分布と近接する代表点を有するボクセルの確率値の総和が最大になるように、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、近接するボクセルが持つ確率値の差の総和が最小になるように、位置合わせする粗合わせステップを有する、ことを特徴とする請求項１に記載の自己位置同定方法。
【請求項５】
前記データ入力ステップの後に、現在の計測位置を過去の計測位置の変化から推定し、もしくは現在の計測位置を取得可能なセンサにより取得し、もしくは計測データの距離値だけでなく反射強度値を利用し、照合する範囲を限定する探索範囲限定ステップを有する、ことを特徴とする請求項１に記載の自己位置同定方法。
【請求項６】
前記自己位置同定ステップにおいて、先の計測位置における位置及び姿勢から新たな計測位置の６自由度位置を同定する、ことを特徴とする請求項１に記載の自己位置同定方法。
【請求項７】
前記精密合わせステップにおいて、誤差分布が交差する場合を同一計測点とし、その場合の距離値に分布の一致度から求めた重みを掛け合わせて誤差分布間の距離を算出する、ことを特徴とする請求項１に記載の自己位置同定方法。
【請求項８】
前記モデル構築ステップにおいて、最大のボクセルを必要最小限の分解能に相当する大きさに設定し、かつ単一のボクセル内に複数の被計測点が存在する場合に、単一のボクセル内に単一の被計測点のみが存在するように、該ボクセルを更に分割して階層的に複数のボクセルに分割する、ことを特徴とする請求項１に記載の自己位置同定方法。
【請求項９】
前記自己位置同定ステップの後に、前記環境モデルを更新するモデル更新ステップを有し、該モデル更新ステップにおいて、新たに入力された被計測点の座標値に対応するボクセルを探索し、
原点と被計測点の間に物体が存在しないものとして、その間に位置するボクセル内の代表点と誤差分布を再設定もしくは消去する、ことを特徴とする請求項１に記載の自己位置同定方法。
【請求項１０】
前記自己位置同定ステップの後に、前記環境モデルを更新するモデル更新ステップを有し、該モデル更新ステップにおいて、新たに入力された被計測点の座標値に対応するボクセルを探索し、
該ボクセル内に代表点がない場合に、前記座標値と誤差分布を代表点の座標値と誤差分布として設定する、ことを特徴とする請求項１に記載の自己位置同定方法。
【請求項１１】
前記自己位置同定ステップの後に、前記環境モデルを更新するモデル更新ステップを有し、該モデル更新ステップにおいて、新たに入力された被計測点の座標値に対応するボクセルを探索し、
前記ボクセル内に既に設定した代表点がある場合に、新たに取得した誤差分布と既に設定したボクセル内の誤差分布とを比較し、
誤差分布が互いに重複する場合に、両誤差分布から新たな誤差分布と新たな代表点を再設定し、
誤差分布が互いに重複しない場合に、単一のボクセル内に単一の代表点のみが存在するように、該ボクセルを更に分割して階層的に複数のボクセルに分割する、ことを特徴とする請求項１に記載の自己位置同定方法。
【請求項１２】
前記自己位置同定ステップにおいて、自己位置の同定とともに、自己位置の誤差分布を同定し、前記出力ステップの前に、現在の自己位置と誤差分布および同定した自己位置と誤差分布とからカルマンフィルタによって、自己位置と誤差分布を補正し、出力ステップにおいて、自己位置と誤差分布を出力する、ことを特徴とする請求項１に記載の自己位置同定方法。
【請求項１３】
前記モデル更新ステップにおいて、前記新たに取得した誤差分布と前記既に設定したボクセル内の誤差分布とを比較し、誤差分布が互いに重複する場合に、両誤差分布から新たな誤差分布と新たな代表点を再設定した結果、新たな代表点が他のボクセル内へ移動したとき、
該他のボクセル内に代表点がない場合に、該新たな誤差分布と該新たな代表点を該他のボクセルの内部に設定し、
該他のボクセル内に既に設定した代表点がある場合に、該新たな誤差分布と既に設定した該他のボクセル内の誤差分布とを比較し、（Ａ）誤差分布が互いに重複する場合に、両誤差分布から、または、両誤差分布とボクセル内に既に設定した代表点と新たに入力された被計測点の座標値から、新たな誤差分布と新たな代表点を再設定し、（Ｂ）誤差分布が互いに重複しない場合に、単一のボクセル内に単一の代表点のみが存在するように、該ボクセルを更に分割して階層的に複数のボクセルに分割する、ことを特徴とする請求項１に記載の自己位置同定方法。
【請求項１４】
前記自己位置同定ステップの後に、前記環境モデルを更新するモデル更新ステップを有し、
該モデル更新ステップにおいて、新たに入力された被計測点の座標値およびその誤差分布と、既に設定したボクセル内の代表点およびその誤差分布とから、カルマンフィルタにより新たな代表点と誤差分布を取得して再設定する、ことを特徴とする請求項１に記載の自己位置同定方法。
【請求項１５】
前記精密合わせステップにおいて、前記近接する誤差分布間の距離に関する評価値が最小になるように位置合わせする代わりに、前記近接する誤差分布間に基づく最尤推定値によって定めた一致度に関する評価値が最大となるように、先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、もしくは新たな計測位置に対する環境モデルを回転及び並進させて位置合わせする、ことを特徴とする請求項１に記載の自己位置同定方法。
【請求項１６】
前記一致度に関する評価値の算出式が次の［数１６］で表わされ、
【数１６】
0004650752.tif 000017
この式において、計測点ｊと環境モデル上の代表点ｉとが対応付けられているとし、当該計測点ｊなる計測データが得られる確率をＥＭ（ｉ、ｊ）としており、ω（ｊ）は、環境モデルの中に計測点ｊと対応付けられる代表点が存在する場合は１、それ以外の場合は０としている、ことを特徴とする請求項１５に記載の自己位置同定方法。
【請求項１７】
外界の三次元形状を取り入れて自己位置を同定する自己位置同定方法であって、
新たな計測位置において、三次元形状上の座標値をコンピュータに入力するデータ入力ステップと、
前記三次元形状の存在する空間領域を、境界表面が互いに直交する直方体からなる複数のボクセルに分割し、各ボクセル位置を記憶する環境モデルを構築するモデル構築ステップと、
前記座標値に対応するボクセルの内部に代表点とその誤差分布を設定し記憶するマッチングステップとを実施し、
先の計測位置おける前記三次元形状のデータが存在しない場合には、新たな計測位置を自己位置と同定し、
先の計測位置おける前記三次元形状のデータが存在する場合には、
先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、近接する誤差分布間の距離に関する評価値が最小になるように位置合わせする精密合わせステップと、
前記精密合わせステップにおける回転量及び並進量から自己位置を同定する自己位置同定ステップとを実施する、ことを特徴とする自己位置同定方法。
【請求項１８】
前記自己位置を出力装置に出力する出力ステップを実施する、ことを特徴とする請求項１７に記載の自己位置同定方法。
【請求項１９】
前記マッチングステップにおいて、ボクセルの内部に代表点とその誤差分布に加え、ボクセル内に物体の存在確率を表す確率値を設定し、記憶する、ことを特徴とする請求項１７に記載の自己位置同定方法。
【請求項２０】
前記精密合わせステップの前に、先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、計測データおよび誤差分布と近接する代表点を有するボクセル間との距離に関する評価値が最小になるように、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、代表点を有するボクセル間の距離に関する評価値が最小になるように、位置合わせする粗合わせステップを有する、ことを特徴とする請求項１７に記載の自己位置同定方法。
【請求項２１】
前記精密合わせステップの前に、先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、計測データおよび誤差分布と近接する代表点を有するボクセルの確率値に関する評価値が最大になるように、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、近接するボクセルが持つ確率値の差に関する評価値が最小になるように、位置合わせする粗合わせステップを有する、ことを特徴とする請求項１７に記載の自己位置同定方法。
【請求項２２】
外界の三次元形状の被計測点の座標値から三次元形状を復元し、三次元形状データを出力する三次元形状の計測方法であって、
新たな計測位置において、三次元形状上の座標値をコンピュータに入力するデータ入力ステップと、
前記三次元形状の存在する空間領域を、境界表面が互いに直交する直方体からなる複数のボクセルに分割し、各ボクセル位置を記憶する環境モデルを構築するモデル構築ステップと、
前記座標値に対応するボクセルの内部に代表点とその誤差分布を設定し記憶するマッチングステップとを実施し、
先の計測位置における前記三次元形状のデータが存在しない場合には、新たな計測位置を自己位置と同定し、
先の計測位置おける前記三次元形状のデータが存在する場合には、
先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、近接する誤差分布間の距離に関する評価値が最小になるように位置合わせする精密合わせステップと、
前記精密合わせステップにおける回転量及び並進量から自己位置を同定する自己位置同定ステップとを実施し、
更に前記自己位置、前記自己位置に基づいた前記ボクセル位置、代表点、および、誤差分布の少なくともいずれかを出力装置に出力する出力ステップを有する、ことを特徴とする三次元形状の計測方法。
【請求項２３】
前記出力ステップにおいて、前記ボクセルの位置、代表点の位置および誤差分布の位置の少なくともいずれかを三次元形状の計測値として出力装置に出力するとともに、該計測値の信頼性または精度を示す指標を、該ボクセルの内部の誤差分布の大きさに基づいて、出力装置に出力する、ことを特徴とする請求項２２に記載の三次元形状の計測方法。
【請求項２４】
前記出力ステップにおいて、前記ボクセルの位置、代表点の位置および誤差分布の位置の少なくともいずれかを三次元形状の計測値として出力装置に出力するときに、該ボクセルの内部の誤差分布の大きさが所定の基準値よりも大きい場合に、該計測値の信頼性または精度が所定の基準よりも低いとして、該ボクセルの前記計測値を出力装置に出力しない、ことを特徴とする請求項２２に記載の三次元形状の計測方法。
【請求項２５】
前記マッチングステップの後に、前記環境モデルを更新するモデル更新ステップを有し、該モデル更新ステップにおいて、新たに入力された被計測点の座標値に対応するボクセルを探索し、
該ボクセル内に代表点がない場合に、前記座標値と誤差分布を代表点の座標値と誤差分布として設定する、ことを特徴とする請求項２２に記載の三次元形状の計測方法。
【請求項２６】
前記マッチングステップの後に、前記環境モデルを更新するモデル更新ステップを有し、該モデル更新ステップにおいて、新たに入力された被計測点の座標値に対応するボクセルを探索し、
前記ボクセル内に既に設定した代表点がある場合に、新たに取得した誤差分布と既に設定したボクセル内の誤差分布とを比較し、
誤差分布が互いに重複する場合に、両誤差分布から、または、両誤差分布とボクセル内に既に設定した代表点と新たに入力された被計測点の座標値から、新たな誤差分布と新たな代表点を再設定し、
誤差分布が互いに重複しない場合に、単一のボクセル内に単一の代表点のみが存在するように、該ボクセルを更に分割して階層的に複数のボクセルに分割する、ことを特徴とする請求項２２に記載の三次元形状の計測方法。
【請求項２７】
前記マッチングステップの後に、前記環境モデルを更新するモデル更新ステップを有し、該モデル更新ステップにおいて、新たに入力された被計測点の座標値に対応するボクセルを探索し、前記ボクセル内の代表点および誤差分布の少なくともいずれかが新たに設定される、または再設定される、または該ボクセルを更に分割して階層的に複数のボクセルに分割される場合、前記出力ステップにおいて、当該ボクセルの代表点の位置を三次元形状の計測値として出力装置に出力する、ことを特徴とする請求項２２に記載の三次元形状の計測方法。
【請求項２８】
前記出力ステップにおいて、距離センサの位置から距離センサが位置計測可能な範囲の環境モデル内のボクセルの代表点の位置を三次元形状の計測値として出力装置に出力する、ことを特徴とする請求項２２に記載の三次元形状の計測方法。
【請求項２９】
外界の三次元形状を取り入れて自己位置を同定する自己位置同定装置であって、
三次元形状上の座標値をコンピュータに入力するデータ入力装置と、
前記三次元形状の存在する空間領域を、境界表面が互いに直交する直方体からなる複数のボクセルに分割し、各ボクセル位置を記憶する環境モデルを構築するモデル構築装置と、
前記座標値に対応するボクセルの内部に代表点とその誤差分布を設定し記憶するマッチング装置と、
自己位置を出力装置に出力する出力装置とを備え、
新たな計測位置において、先の計測位置おける前記三次元形状のデータが存在しない場合には、新たな計測位置を自己位置と同定し、
新たな計測位置において、先の計測位置おける前記三次元形状のデータが存在する場合には、
先の計測位置に対する環境モデルに対し、新たな計測位置に対する環境モデルを回転及び並進させ、近接する誤差分布間の距離の総和が最小になるように位置合わせし、
前記位置合わせにおける回転量及び並進量から自己位置を同定する、ことを特徴とする自己位置同定装置。
【請求項３０】
外界の三次元形状を取り入れて自己位置を同定する自己位置同定装置であって、
三次元形状上の座標値をコンピュータに入力するデータ入力装置と、
前記三次元形状の存在する空間領域を、境界表面が互いに直交する直方体からなる複数のボクセルに分割し、各ボクセル位置を記憶する環境モデルを構築するモデル構築装置と、
前記座標値に対応するボクセルの内部に代表点とその誤差分布を設定し記憶するマッチング装置と、
自己位置を出力装置に出力する出力装置とを備え、
新たな計測位置において、先の計測位置おける前記三次元形状のデータが存在しない場合には、新たな計測位置を自己位置と同定し、
新たな計測位置において、先の計測位置おける前記三次元形状のデータが存在する場合には、
先の計測位置に対する環境モデルに対し、新たな計測位置に対する環境モデルを回転及び並進させ、近接する誤差分布間の距離に関する評価値が最小になるように位置合わせし、
前記位置合わせにおける回転量及び並進量から自己位置を同定する、ことを特徴とする自己位置同定装置。
【請求項３１】
三次元形状上の被計測点の座標値から三次元形状を復元し、三次元形状データを出力する三次元形状の計測装置であって、
三次元形状上の座標値をコンピュータに入力するデータ入力装置と、
前記三次元形状の存在する空間領域を、境界表面が互いに直交する直方体からなる複数のボクセルに分割し、各ボクセル位置を記憶する環境モデルを構築するモデル構築装置と、
前記座標値に対応するボクセルの内部に代表点とその誤差分布を設定し記憶するマッチング装置と、
新たな計測位置において、先の計測位置おける前記三次元形状のデータが存在しない場合には、新たな計測位置を自己位置と同定し、
新たな計測位置において、先の計測位置おける前記三次元形状のデータが存在する場合には、
先の計測位置に対する環境モデルに対し、新たな計測位置に対する環境モデルを回転及び並進させ、近接する誤差分布間の距離に関する評価値が最小になるように位置合わせし、
前記位置合わせにおける回転量及び並進量から自己位置を同定し、
更に前記自己位置、前記自己位置に基づいた前記ボクセル位置、代表点、および、誤差分布の少なくともいずれかを出力装置に出力するデータ伝達装置とを備える、ことを特徴とする三次元形状の計測装置。
【発明の詳細な説明】
【発明の背景】
【０００１】
発明の技術分野
本発明は、外界の三次元情報を取り入れて、自己位置を同定する自己位置同定方法と装置および三次元形状の計測方法と装置に関する。
【関連技術の説明】
【０００２】
自己位置を直接同定する手段として、ＧＰＳ、ＤＧＰＳ、磁気マーカ、磁気方式、無線方式、画像マーカ方式、光学式・電磁式マーカ方式等が知られている。しかしこれらの手段は、安定した位置同定が可能であるというメリットはあるが、インフラが必要であり、位置検出が可能な範囲が限られ、設置のための手間を必要とする等の問題点がある。
【０００３】
これに対して、インフラを用いずにセンサにより周囲環境を計測して獲得した環境情報と、新たに獲得した計測データを比較しながら自己位置を算出する自己位置同定手段が種々提案されている。
この目的に適した三次元センサとして、超音波センサ、画像センサ、距離センサが挙がられる。しかし、超音波センサは、計測範囲が狭く精度が低い問題点があり、画像センサは照明条件に影響を受けやすい、という問題点がある。
【０００４】
これらに対し、距離センサは超音波センサに比べて広範囲の距離データを精度良く計測でき、かつ、画像処理のように照明の変化に影響を受けることがないため、外乱光の影響を受けずに安定した計測を行うことができる利点がある。
【０００５】
距離センサを用いた位置同定の手段は、非特許文献１、特許文献１、等に提案されている。
【０００６】
非特許文献１の手段は、三次元スキャン式レーザセンサを用いて周囲環境計測により獲得した環境モデルと、逐次得られる計測データをＩＣＰアルゴリズムにより比較しながら位置合わせを行うものである。
ＩＣＰアルゴリズムは、対応点が既知でない場合の位置合わせ手段の１つであり、先の計測データに対するその後の計測データの最も近傍の点を求め、回転・並進を行いながらその距離の総和が最小になる状態を一致状態となるように解を求めるものである（非特許文献２）。
【０００７】
特許文献１の「自己位置同定装置」は、特に無軌道でかつ外界にランドマーク等の設備を設けないものであり、あらかじめ移動体を目的とする軌道に沿って手動で動作させる（ティーチング）ことで、外環境の地図を自動的に記憶させるものである。
この特許文献１では、図１に示すように、平面内を移動する移動体１０１の、該平面内における周囲外環境物体１０２との位置関係を、離散化したデータとしてメモリ上に保持する外環境地図を用いる。前記データは前記移動体上に固定された点を原点１０３とし、所定の間隔に離散化した前記移動体からの方向データ１０４と、前記方向における所定の間隔に離散化した前記周囲外環境物体までの距離データ１０５との組からなる。
【０００８】
また、本発明に関連する技術として、非特許文献３が開示されている。
【０００９】
【非特許文献１】Ａｎｄｒｅａｓ  Ｎｕｃｈｔｅｒ，ｅｔ  ａｌ．，“６Ｄ  ＳＬＡＭ  ｗｉｔｈ  ａｎ  Ａｐｐｌｉｃａｔｉｏｎ  ｉｎ  Ａｕｔｏｎｏｍｏｕｓ  Ｍｉｎｅ  Ｍａｐｐｉｎｇ”，ＩＥＥＥ  ２００４  Ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｃｏｎｆｅｒｅｎｃｅ  Ｒｏｂｏｔｉｃｓ  ａｎｄ  Ａｕｔｏｍａｔｉｏｎ  （ＩＣＲＡ  ’０４），  Ａｐｒｉｌ  ２００４
【非特許文献２】Ｐａｕｌ  Ｊ．Ｂｅｓｌ，“Ａ  Ｍｅｔｈｏｄ  ｆｏｒ  Ｒｅｇｉｓｔｒａｔｉｏｎ  ｏｆ  ３−Ｄ  Ｓｈａｐｅｓ”，ＩＥＥＥ  Ｔｒａｎｓａｃｔｉｏｎｓ  ｏｆ  Ｐａｔｔｅｒｎ  Ａｎａｌｙｓｉｓ  ａｎｄ  Ｍｅｃｈａｎｉｃａｌ  Ｉｎｔｅｌｌｉｇｅｎｃｅ，  Ｖｏｌ．１４，  Ｎｏ．２，  Ｆｅｂｒｕａｒｙ  １９９２
【非特許文献３】関本清英、他「三次元レーザレーダの開発」、石川島播磨技報Ｖｏｌ．４３  Ｎｏ．４（２００３−７）
【００１０】
【特許文献１】特開２００３?１５７３９号公報、「外環境地図、並びに自己位置同定装置および誘導制御装置」
【発明の開示】
【発明が解決しようとする課題】
【００１１】
三次元レーザレーダのような距離センサを用いる場合、計測される三次元形状上の被計測点は、横方向及び垂直方向に離散した点群となる。この点群の間隔は、計測点からの距離が例えば５０ｍの場合、被計測点の間隔は、例えば横方向で約３１５ｍｍ、垂直方向で約５２５ｍｍに達する。【００１２】
また、複数の計測位置から静止している三次元形状を計測する場合、三次元レーザレーダのような距離センサでは、計測位置毎に、被計測点の位置は通常異なる。
さらに、このような距離センサは、一般に測定距離に例えば約２０ｃｍ前後の誤差を有する。
従って、三次元レーザレーダのような距離センサを用いる場合、以下の制約条件Ａ〜Ｃがある。
条件Ａ：計測データに点数が少ない（例えば、１フレーム＝１６６×５０点）
条件Ｂ：計測データに誤差を含む（例えば、測定距離に約２０ｃｍ前後）
条件Ｃ：計測データは同じ計測点を計るとは限らない。
すなわち、得られる距離データは横方向及び垂直方向に離散した点群であり、計測毎に位置が相違するため対応点がなく、測定距離に比較的大きな誤差を含んでいる。【００１３】
このような距離データを上述した「ＩＣＰアルゴリズム」で処理する場合、以下の問題点がある。
（１）誤差の蓄積
ＩＣＰアルゴリズムは、２つの距離データの重ね合わせ手段であり、先のデータとその後のデータとの比較を繰り返し、その差分を積分しても、対応点がほとんどないため、誤差が蓄積してしまう。
（２）計算量が多い
ＩＣＰアルゴリズムは、繰り返し計算であるため、計算量が膨大となる。すなわち、ＩＣＰアルゴリズムは、計測データの各データ点に対応するモデルデータを探索する必要があるため、モデルデータ点数および計測データ点数が増加すると計算量が増大する。具体的には、モデルデータの点数をＭ、計測データの点数をＮとした場合、例えば全探索時の計算オーダは、Ｏ（Ｍ×Ｎ）となる。
（３）計測点が少ない場合を扱えない
ＩＣＰアルゴリズムは、密な距離データを対象としているため、離散した点群であり、空間的に疎である場合は、誤った結果に収束してしまう。【００１４】
そのため、このような距離センサを用いた位置同定には、以下の要件を満たす必要がある。
（１）メモリ効率の良い周囲環境のためのデータ構造
逐次得られる計測データを全て保存する方法では無尽蔵にメモリを必要とする。よって、周囲環境を計測した結果を効率良く格納するデータ構造が必要である。
（２）計測データの点数不足および誤差に対する安定性
計測データの点数が少ない場合、および、誤差を含む場合であっても、できるだけ位置同定精度が低下しない必要がある。
（３）位置同定のための計算の効率化
周囲環境を計測して獲得した環境情報と計測データの比較により、自己位置同定を行うが、この比較処理に多くの計算を必要とする。
【発明の要約】
【００１５】
本発明は上述した問題点を解決するために創案されたものである。すなわち、本発明の目的は、メモリ効率の良い周囲環境のためのデータ構造を有し、計測データの点数不足および誤差に対する安定性があり、位置同定のための計算が効率的にでき、これにより、外界の三次元情報が空間的に離散した点群であり、計測毎に位置が相違して対応点がなく、比較的大きな誤差を含んでいる場合でも、誤差の蓄積がなく、少ない計算量で、正確な位置合わせができる自己位置同定方法と装置および三次元形状の計測方法と装置を提供することにある。
【００１６】
本発明によれば、外界の三次元形状を取り入れて自己位置を同定する自己位置同定方法であって、
新たな計測位置において、三次元形状上の座標値をコンピュータに入力するデータ入力ステップと、
前記三次元形状の存在する空間領域を、境界表面が互いに直交する直方体からなる複数のボクセルに分割し、各ボクセル位置を記憶する環境モデルを構築するモデル構築ステップと、
前記座標値に対応するボクセルの内部に代表点とその誤差分布を設定し記憶するマッチングステップとを実施し、
先の計測位置おける前記三次元形状のデータが存在しない場合には、新たな計測位置を自己位置と同定し、
先の計測位置おける前記三次元形状のデータが存在する場合には、
先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、近接する誤差分布間の距離に関する評価値（例えば、当該距離の総和）が最小になるように位置合わせする精密合わせステップと、
前記精密合わせステップにおける回転量及び並進量から自己位置を同定する自己位置同定ステップを実施する、ことを特徴とする自己位置同定方法が提供される。なお、更に前記自己位置を出力装置に出力する出力ステップを実施してもよい。
また、本発明によれば、外界の三次元形状の被計測点の座標値から三次元形状を復元し、三次元形状データを出力する三次元形状の計測方法であって、
新たな計測位置において、三次元形状上の座標値をコンピュータに入力するデータ入力ステップと、
前記三次元形状の存在する空間領域を、境界表面が互いに直交する直方体からなる複数のボクセルに分割し、各ボクセル位置を記憶する環境モデルを構築するモデル構築ステップと、
前記座標値に対応するボクセルの内部に代表点とその誤差分布を設定し記憶するマッチングステップとを実施し、
先の計測位置における前記三次元形状のデータが存在しない場合には、新たな計測位置を自己位置と同定し、
先の計測位置おける前記三次元形状のデータが存在する場合には、
先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、近接する誤差分布間の距離に関する評価値が最小になるように位置合わせする精密合わせステップと、
前記精密合わせステップにおける回転量及び並進量から自己位置を同定する自己位置同定ステップとを実施し、
更に前記自己位置、前記自己位置に基づいた前記ボクセル位置、代表点、および、誤差分布の少なくともいずれかを出力装置に出力する出力ステップを有する、ことを特徴とする三次元形状の計測方法が提供される。
なお、上記自己位置同定方法または三次元形状の計測方法において、上記距離に関する評価値は、当該距離の総和の代わりに、当該距離の平均値、当該距離の２乗の総和、または、当該距離の最大値であってもよく、他の適切な評価値であってもよい。
【００１７】
上記自己位置同定方法または三次元形状の計測方法は、以下の内容を含んでもよい。
【００１８】
本発明の好ましい実施形態によれば、前記自己位置同定ステップにおいて、先の計測位置における位置及び姿勢から新たな計測位置の６自由度位置を同定する。
【００１９】
また、前記マッチングステップにおいて、ボクセルの内部に代表点とその誤差分布に加え、ボクセル内に物体の存在確率を表す確率値を設定し、記憶する。
【００２０】
また、前記精密合わせステップの前に、先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、計測データおよび誤差分布と近接する代表点を有するボクセル間との距離に関する評価値（例えば、当該距離の総和）が最小になるように、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、代表点を有するボクセル間の距離に関する評価値（例えば、当該距離の総和）が最小になるように、位置合わせする粗合わせステップを有する。
なお、計測データおよび誤差分布と近接する代表点を有するボクセル間との距離に関する上記評価値は、当該距離の総和の代わりに、当該距離の平均値、当該距離の２乗の総和、または、当該距離の最大値であってもよく、他の適切な評価値であってもよい。
また、代表点を有するボクセル間の距離に関する上記評価値は、当該距離の総和の代わりに、当該距離の総和の代わりに、当該距離の平均値、当該距離の２乗の総和、または、当該距離の最大値であってもよく、他の適切な評価値であってもよい。
【００２１】
または、前記マッチングステップにおいて、ボクセルの内部に物体の存在確率を表す確率値を設定し、記憶する手段を有する場合、前記精密合わせステップの前に、先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、計測データおよび誤差分布と近接する代表点を有するボクセルの確率値に関する評価値（例えば、当該確率値の総和）が最大になるように、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、近接するボクセルが持つ確率値の差に関する評価値（例えば、当該確率値の差の総和）が最小になるように、位置合わせする粗合わせステップを有する。
なお、計測データおよび誤差分布と近接する代表点を有するボクセルの確率値に関する上記評価値は、当該確率値の総和の代わりに、当該確率値の平均値、または、当該確率値の２乗の総和、または、当該確率値の最小値であってもよく、他の適切な評価値であってもよい。
また、近接するボクセルが持つ確率値の差に関する上記評価値は、当該確率値の差の総和の代わりに、当該確率値の差の平均値、当該確率値の差の２乗の総和、または、当該確率値の差の最大値であってもよく、他の適切な評価値であってもよい。
【００２２】
また、前記データ入力ステップの後に、現在の計測位置を過去の計測位置の変化から推定し、もしくは現在の計測位置を取得可能なセンサにより取得し、もしくは計測データの距離値だけでなく反射強度値を利用し、照合する範囲を限定する探索範囲限定ステップを有する。
【００２３】
また、前記精密合わせステップにおいて、誤差分布が交差する場合を同一計測点とし、その場合の距離値に分布の一致度から求めた重みを掛け合わせて誤差分布間の距離を算出する。
【００２４】
また、前記モデル構築ステップにおいて、最大のボクセルを必要最小限の分解能に相当する大きさに設定し、かつ単一のボクセル内に複数の被計測点が存在する場合に、単一のボクセル内に単一の被計測点のみが存在するように、該ボクセルを更に分割して階層的に複数のボクセルに分割する。
【００２５】
また、前記自己位置同定ステップの後に、前記環境モデルを更新するモデル更新ステップを有し、該モデル更新ステップにおいて、新たに入力された被計測点の座標値に対応するボクセルを探索し、
原点と被計測点の間に物体が存在しないものとして、その間に位置するボクセル内の代表点と誤差分布を再設定もしくは消去する。
【００２６】
また、前記自己位置同定ステップの後に、前記環境モデルを更新するモデル更新ステップを有し、該モデル更新ステップにおいて、新たに入力された被計測点の座標値に対応するボクセルを探索し、
該ボクセル内に代表点がない場合に、前記座標値と誤差分布を代表点の座標値と誤差分布として設定する。
【００２７】
また、自己位置同定ステップの後に、前記環境モデルを更新するモデル更新ステップを有し、該モデル更新ステップにおいて、新たに入力された被計測点の座標値に対応するボクセルを探索し、
前記ボクセル内に既に設定した代表点がある場合に、新たに取得した誤差分布と既に設定したボクセル内の誤差分布とを比較し、
誤差分布が互いに重複する場合に、両誤差分布から、または、両誤差分布とボクセル内に既に設定した代表点と新たに入力された被計測点の座標値から、新たな誤差分布と新たな代表点を再設定し、
誤差分布が互いに重複しない場合に、単一のボクセル内に単一の代表点のみが存在するように、該ボクセルを更に分割して階層的に複数のボクセルに分割する。ボクセルの分割には、例えば八分木やＫ−Ｄ木を用いる。
また、自己位置の同定とともに、自己位置の誤差分布を同定し、現在の自己位置と誤差分布および同定した自己位置と誤差分布とからカルマンフィルタによって、自己位置と誤差分布を補正し、自己位置と誤差分布を出力する。
【００２８】
前記モデル更新ステップにおいて、前記新たに取得した誤差分布と前記既に設定したボクセル内の誤差分布とを比較し、誤差分布が互いに重複する場合に、両誤差分布から新たな誤差分布と新たな代表点を再設定した結果、新たな代表点が他のボクセル内へ移動したとき、
該他のボクセル内に代表点がない場合に、該新たな誤差分布と該新たな代表点を該他のボクセルの内部に設定し、
該他のボクセル内に既に設定した代表点がある場合に、該新たな誤差分布と既に設定した該他のボクセル内の誤差分布とを比較し、（Ａ）誤差分布が互いに重複する場合に、両誤差分布から、または、両誤差分布とボクセル内に既に設定した代表点と新たに入力された被計測点の座標値から、新たな誤差分布と新たな代表点を再設定し、（Ｂ）誤差分布が互いに重複しない場合に、単一のボクセル内に単一の代表点のみが存在するように、該ボクセルを更に分割して階層的に複数のボクセルに分割する。
【００２９】
また、前記自己位置同定ステップの後に、前記環境モデルを更新するモデル更新ステップを有し、
該モデル更新ステップにおいて、新たに入力された被計測点の座標値およびその誤差分布と、既に設定したボクセル内の代表点およびその誤差分布とから、カルマンフィルタにより新たな代表点と誤差分布を取得して再設定する。
【００３０】
本発明の別の実施形態によると、前記精密合わせステップにおいて、前記近接する誤差分布間の距離に基づく評価値が最小になるように位置合わせする代わりに、前記近接する誤差分布間に基づく最尤推定値によって定めた一致度に関する評価値（例えば、当該一致度の総乗）が最大となるように、先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、もしくは新たな計測位置に対する環境モデルを回転及び並進させて位置合わせする。
【００３１】
前記一致度に関する評価値の算出式が次の［数１］で表わされ、
【数１】
0004650752.tif 000002
【００３２】
この式において、計測点ｊと環境モデル上の代表点ｉとが対応付けられているとし、当該計測点ｊなる計測データが得られる確率をＥＭ（ｉ、ｊ）としており、ω（ｊ）は、環境モデルの中に計測点ｊと対応付けられる代表点が存在する場合は１、それ以外の場合は０としている。
【００３３】
前記出力ステップにおいて、前記ボクセルの位置、代表点の位置および誤差分布の位置の少なくともいずれかを三次元形状の計測値として出力装置に出力するとともに、該計測値の信頼性または精度を示す指標を、該ボクセルの内部の誤差分布の大きさに基づいて、出力装置に出力する。
【００３４】
前記出力ステップにおいて、前記ボクセルの位置、代表点の位置および誤差分布の位置の少なくともいずれかを三次元形状の計測値として出力装置に出力するときに、該ボクセルの内部の誤差分布の大きさが所定の基準値よりも大きい場合に、該計測値の信頼性または精度が所定の基準よりも低いとして、該ボクセルの前記計測値を出力装置に出力しない。
【００３５】
前記マッチングステップの後に、前記環境モデルを更新するモデル更新ステップを有し、該モデル更新ステップにおいて、新たに入力された被計測点の座標値に対応するボクセルを探索し、前記ボクセル内の代表点および誤差分布の少なくともいずれかが新たに設定される、または再設定される、または該ボクセルを更に分割して階層的に複数のボクセルに分割される場合、前記出力ステップにおいて、当該ボクセルの代表点の位置を三次元形状の計測値として出力装置に出力する。
【００３６】
前記出力ステップにおいて、距離センサの位置から距離センサが位置計測可能な範囲の環境モデル内のボクセルの代表点の位置を三次元形状の計測値として出力装置に出力する。
【００３７】
また、本発明によれば、外界の三次元形状を取り入れて自己位置を同定する自己位置同定装置であって、
三次元形状上の座標値をコンピュータに入力するデータ入力装置と、
前記三次元形状の存在する空間領域を、境界表面が互いに直交する直方体からなる複数のボクセルに分割し、各ボクセル位置を記憶する環境モデルを構築するモデル構築装置と、
前記座標値に対応するボクセルの内部に代表点とその誤差分布を設定し記憶するマッチング装置と、
自己位置を出力装置に出力するデータ伝達装置とを備え、
新たな計測位置において、先の計測位置おける前記三次元形状のデータが存在しない場合には、新たな計測位置を自己位置と同定し、
新たな計測位置において、先の計測位置おける前記三次元形状のデータが存在する場合には、
先の計測位置に対する環境モデルに対し、新たな計測位置に対する環境モデルを回転及び並進させ、近接する誤差分布間の距離に関する評価値（例えば、当該距離の総和）が最小になるように位置合わせし、
前記位置合わせにおける回転量及び並進量から自己位置を同定する、ことを特徴とする自己位置同定装置が提供される。
また、本発明によれば、三次元形状上の被計測点の座標値から三次元形状を復元し、三次元形状データを出力する三次元形状の計測装置であって、
三次元形状上の座標値をコンピュータに入力するデータ入力装置と、
前記三次元形状の存在する空間領域を、境界表面が互いに直交する直方体からなる複数のボクセルに分割し、各ボクセル位置を記憶する環境モデルを構築するモデル構築装置と、
前記座標値に対応するボクセルの内部に代表点とその誤差分布を設定し記憶するマッチング装置と、
新たな計測位置において、先の計測位置おける前記三次元形状のデータが存在しない場合には、新たな計測位置を自己位置と同定し、
新たな計測位置において、先の計測位置おける前記三次元形状のデータが存在する場合には、
先の計測位置に対する環境モデルに対し、新たな計測位置に対する環境モデルを回転及び並進させ、近接する誤差分布間の距離に関する評価値が最小になるように位置合わせし、
前記位置合わせにおける回転量及び並進量から自己位置を同定し、
更に前記自己位置、前記自己位置に基づいた前記ボクセル位置、代表点、および、誤差分布の少なくともいずれかを出力装置に出力するデータ伝達装置とを備える、ことを特徴とする三次元形状の計測装置が提供される。
【００３８】
前記自己位置同定装置または三次元形状の計測装置は、以下の内容を含んでもよい。
【００３９】
前記マッチング装置は、ボクセルの内部に代表点とその誤差分布に加え、ボクセル内に物体の存在確率を表す確率値を設定し、記憶する。
【００４０】
前記自己位置同定装置または三次元形状の計測装置は、前記位置合わせ（精密合わせステップ）の前に、先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、計測データおよび誤差分布と近接する代表点を有するボクセル間との距離に関する評価値（例えば、当該距離の総和）が最小になるように、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、代表点を有するボクセル間の距離に関する評価値（例えば、当該距離の総和）が最小になるように、位置合わせをする粗合わせステップを行う。
【００４１】
前記自己位置同定装置または三次元形状の計測装置は、前記精密合わせステップの前に、先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、計測データおよび誤差分布と近接する代表点を有するボクセルの確率値に関する評価値（例えば、当該確率値の総和）が最大になるように、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、近接するボクセルが持つ確率値の差に関する評価値（例えば、当該確率値の差の総和）が最小になるように、位置合わせする粗合わせステップを行う。
【００４２】
前記自己位置同定装置または三次元形状の計測装置は、前記データ入力ステップの後に、現在の計測位置を過去の計測位置の変化から推定し、もしくは現在の計測位置を取得可能なセンサにより取得し、もしくは計測データの距離値だけでなく反射強度値を利用し、照合する範囲を限定する。
【００４３】
前記自己位置同定装置または三次元形状の計測装置は、先の計測位置における位置及び姿勢から新たな計測位置の６自由度位置を同定する。
【００４４】
前記自己位置同定装置または三次元形状の計測装置は、誤差分布が交差する場合を同一計測点とし、その場合の距離値に分布の一致度から求めた重みを掛け合わせて誤差分布間の距離を算出する。
【００４５】
前記モデル構築装置は、最大のボクセルを必要最小限の分解能に相当する大きさに設定し、かつ単一のボクセル内に複数の被計測点が存在する場合に、単一のボクセル内に単一の被計測点のみが存在するように、該ボクセルを更に分割して階層的に複数のボクセルに分割する。
【００４６】
前記自己位置同定装置または三次元形状の計測装置は、前記自己位置の同定（自己位置同定ステップ）の後に、前記環境モデルを更新するモデル更新装置を有し、
該モデル更新装置は、
新たに入力された被計測点の座標値に対応するボクセルを探索し、
原点と被計測点の間に物体が存在しないものとして、その間に位置するボクセル内の代表点と誤差分布を再設定もしくは消去する。
【００４７】
前記自己位置同定装置または三次元形状の計測装置は、前記自己位置同定ステップの後に、前記環境モデルを更新するモデル更新装置を有し、
該モデル更新装置は、新たに入力された被計測点の座標値に対応するボクセルを探索し、
該ボクセル内に代表点がない場合に、前記座標値と誤差分布を代表点の座標値と誤差分布として設定する。
【００４８】
前記自己位置同定装置または三次元形状の計測装置は、前記自己位置同定ステップの後に、前記環境モデルを更新するモデル更新ステップを有し、該モデル更新ステップにおいて、新たに入力された被計測点の座標値に対応するボクセルを探索し、
前記ボクセル内に既に設定した代表点がある場合に、新たに取得した誤差分布と既に設定したボクセル内の誤差分布とを比較し、
誤差分布が互いに重複する場合に、両誤差分布から、または、両誤差分布とボクセル内に既に設定した代表点と新たに入力された被計測点の座標値から、新たな誤差分布と新たな代表点を再設定し、
誤差分布が互いに重複しない場合に、単一のボクセル内に単一の代表点のみが存在するように、該ボクセルを更に分割して階層的に複数のボクセルに分割する。
【００４９】
前記自己位置同定装置または三次元形状の計測装置は、前記自己位置同定ステップにおいて、自己位置の同定とともに、自己位置の誤差分布を同定し、前記出力ステップの前に、現在の自己位置と誤差分布および同定した自己位置と誤差分布とからカルマンフィルタによって、自己位置と誤差分布を補正し、前記データ伝達装置は、自己位置と誤差分布を出力装置に出力する。
【００５０】
前記モデル更新装置は、前記新たに取得した誤差分布と前記既に設定したボクセル内の誤差分布とを比較し、誤差分布が互いに重複する場合に、両誤差分布から新たな誤差分布と新たな代表点を再設定した結果、新たな代表点が他のボクセル内へ移動したとき、
該他のボクセル内に代表点がない場合に、該新たな誤差分布と該新たな代表点を該他のボクセルの内部に設定し、
該他のボクセル内に既に設定した代表点がある場合に、該新たな誤差分布と既に設定した該他のボクセル内の誤差分布とを比較し、（Ａ）誤差分布が互いに重複する場合に、両誤差分布から、または、両誤差分布とボクセル内に既に設定した代表点と新たに入力された被計測点の座標値から、新たな誤差分布と新たな代表点を再設定し、（Ｂ）誤差分布が互いに重複しない場合に、単一のボクセル内に単一の代表点のみが存在するように、該ボクセルを更に分割して階層的に複数のボクセルに分割する。
【００５１】
前記自己位置同定装置または三次元形状の計測装置は、前記自己位置同定ステップの後に、前記環境モデルを更新するモデル更新装置を有し、
該モデル更新装置は、新たに入力された被計測点の座標値と、既に設定したボクセル内の代表点およびその誤差分布とから、カルマンフィルタにより新たな代表点と誤差分布を取得して再設定する。
【００５２】
前記自己位置同定装置または三次元形状の計測装置は、前記精密合わせステップにおいて、前記近接する誤差分布間の距離に関する評価値（例えば、当該距離の総和）が最小になるように位置合わせする代わりに、前記近接する誤差分布間に基づく最尤推定値によって定めた一致度に関する評価値（例えば、当該一致度の総乗）が最大となるように、先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、もしくは新たな計測位置に対する環境モデルを回転及び並進させて位置合わせする。この場合、前記一致度に関する評価値の算出式が上記［数１］の式で表わされる。
【００５３】
前記データ伝達装置は、前記ボクセルの位置、代表点の位置および誤差分布の位置の少なくともいずれかを三次元形状の計測値として出力装置に出力するとともに、該計測値の信頼性または精度を示す指標を、該ボクセルの内部の誤差分布の大きさに基づいて、出力装置に出力する。
【００５４】
前記データ伝達装置は、前記ボクセルの位置、代表点の位置および誤差分布の位置の少なくともいずれかを三次元形状の計測値として出力装置に出力するときに、該ボクセルの内部の誤差分布の大きさが所定の基準値よりも大きい場合に、該計測値の信頼性または精度が所定の基準よりも低いとして、該ボクセルの前記計測値を出力装置に出力しない。
【００５５】
前記データ入力装置は、距離センサを用いて、三次元形状上の座標値を任意の計測位置を原点とする距離データとして原点を移動しながら順次取得する。
【００５６】
前記マッチング装置の後に、前記環境モデルを更新するモデル更新装置を有し、該モデル更新装置において、新たに入力された被計測点の座標値に対応するボクセルを探索し、前記ボクセル内の代表点および誤差分布の少なくともいずれかが新たに設定される、または再設定される、または該ボクセルを更に分割して階層的に複数のボクセルに分割される場合、前記出力装置は、当該ボクセルの代表点の位置を三次元形状の計測値として出力装置に出力する。
【００５７】
前記データ伝達装置は、距離センサの位置から距離センサが位置計測可能な範囲の環境モデル内のボクセルの代表点の位置を三次元形状の計測値として出力装置に出力する。
【００５８】
発明の効果
上記本発明の方法と装置によれば、三次元形状の存在する空間領域を、複数のボクセルに分割し、各ボクセル位置を記憶するので、計測対象物が大きい場合であっても、データ量をボクセル数に比例する小さいデータサイズに抑えることができる。
また、座標値に対応するボクセルの内部に代表点とその誤差分布を設定し記憶するので、ボクセルの分解能以上の情報を表現することができる。
従って、本発明のデータ構造により、複数視点の計測データを固定サイズに統合できる。
【００５９】
また、ボクセルの内部に物体の存在確率を表す確率値を設定し、保存することによって、誤差分布が代表点の属するボクセルよりも広がっている場合においても、各ボクセルにおける物体の存在有無を代表点が属するボクセルを見つけ、その誤差分布から再計算させることなく、当該ボクセルの確率値だけで容易に判断できるので、検索時間を抑えることができる。
また、現在の計測位置を過去の計測位置の変化から推定し、もしくは現在の計測位置を取得可能なセンサにより取得し、もしくは計測データの距離値だけでなく反射強度を利用し、照合する範囲を制限するので、検索時間を抑えることができる。
【００６０】
また、粗合わせステップにおいて、先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、計測データおよび誤差分布と近接する代表点を有するボクセル間との距離に関する評価値（例えば、当該距離の総和）が最小になるように、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、代表点を有するボクセル間の距離に関する評価値（例えば、当該距離の総和）が最小になるように、
もしくは新たな計測データおよび誤差分布を回転及び並進させ、計測データおよび誤差分布と近接する代表点を有するボクセルの確率値に関する評価値（例えば、当該確率値の総和）が最大になるように、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、近接するボクセルが持つ確率値の差に関する評価値（例えば、当該確率値の差の総和）が最小になるように、位置合わせするので、誤差の蓄積を防ぎながら、短時間に代表点を有するボクセル同士の位置合わせができる。
【００６１】
また、先の計測位置に対する環境モデルに対し、新たな計測位置に対する環境モデルを回転及び並進させ、近接するボクセルが持つ確率値の差に関する評価値（例えば、当該確率値の差の総和）が最小になるように、位置合わせする場合には、物体が存在しない情報も加味して位置合わせするので、精度の向上が図れる。
【００６２】
次いで、精密合わせステップにおいて、先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、近接する誤差分布間の距離に関する評価値（例えば、当該距離の総和）が最小になるように位置合わせするので、誤差分布同士の精密な位置合わせを短時間にできる。
従って、本発明による複数視点からのデータの統合処理により、誤差の蓄積を防ぎながら、高精度な形状取得を可能にできる。
【００６３】
また、本発明のデータ構造は、ボクセル構造を拡張した構造をしているため、データサイズは点群に比べて小さく済ませることができる。すなわち、本発明で提案するデータ構造はボクセル内に１つの代表点を保存するため、計測点に対応するモデル点を探索する計算オーダを１とすることができるため、全体での計算オーダをＯ（Ｎ）に減少することができる。
従って、周囲環境と計測データの位置合わせ（ＩＣＰアルゴリズム）を行う際に探索対象となるデータ量が減少することから、計算を効率化できる。
【００６４】
さらに、従来のＩＣＰアルゴリズムは疎データに対して誤った結果を出力するが、本発明の環境モデルは、ボクセル内に代表点と誤差分布を持っているため、疎データに対応した位置合わせが可能である。
【００６５】
また、現在の自己位置と誤差分布および同定した自己位置と誤差分布とからカルマンフィルタによって、自己位置と誤差分布を補正することによって、より自己位置の精度を向上させることができる。
【００６６】
従って、本発明の方法と装置によれば、誤差を含むデータを正確な情報に補正する機能を有すると共に、長時間の計測に対して高精度に収束する。なおかつ、位置同定処理は、ボクセルの各ボクセルに対応する代表点とその誤差分布を計測点で更新する処理であるため計算量が小さく、演算は周囲ボクセルへの影響を与えずボクセル内で閉じている。よって、高速な処理が可能である。また、計測データはボクセルに逐次統合可能であり、その結果得られた周囲環境情報のメモリサイズは固定サイズを上回ることはない。
【００６７】
さらに、前記モデル更新ステップにおいて、新たに入力された被計測点の座標値およびその誤差分布と、既に設定したボクセル内の代表点およびその誤差分布とから、カルマンフィルタにより新たな代表点と誤差分布を取得して再設定するので、より真値に近い形状を得ることができる。
特に、カルマンフィルタを用いたモデル更新ステップを繰り返すことで、誤差を含むデータであってもカルマンフィルタの効果により真値に収束した高精度な形状が得られる。
【００６８】
また、前記精密合わせステップにおいて、前記近接する誤差分布間の距離に関する評価値（例えば、当該距離の総和）が最小になるように位置合わせする代わりに、前記近接する誤差分布間に基づく最尤推定値によって定めた一致度に関する評価値（例えば、当該一致度の総乗）が最大となるように、先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、もしくは新たな計測位置に対する環境モデルを回転及び並進させて位置合わせすることで、環境モデル、計測データの双方の誤差を考慮して位置合わせすることができる。
【００６９】
また、前記出力ステップにおいて、前記ボクセルの代表点の位置を三次元形状の計測値として出力装置に出力する際に、該計測値の信頼性または精度を示す指標を、該ボクセルの内部の誤差分布の大きさ出力装置に出力するため、計測装置を使用する際に、信頼性の低い計測値をアプリケーションの内容に応じて使用者が取捨選択することができるようになる。
【００７０】
また、前記出力ステップにおいて、前記ボクセルの代表点の位置を三次元形状の計測値として出力装置に出力するときに、該ボクセルの内部の誤差分布の大きさが所定の基準値よりも大きい場合に、該計測値の信頼性または精度が所定の基準よりも低いとして、該ボクセルの前記計測値を出力装置に出力しないようにすることによって、計測装置を使用する際に、そもそも信頼性の高い計測値のみを扱うことができるようになるため、扱うデータ量を削減することや、信頼性の向上につながる。
【００７１】
また、前記データ入力ステップにおいて、距離センサを用いて、三次元形状上の座標値を任意の計測位置を原点とする距離データとして原点を移動しながら順次取得し、特に異なる方向から距離データを得ることによって、分布形状が異なる誤差分布の距離データを統合することが可能となり、精度を向上することができる。
【００７２】
また、前記マッチングステップの後に、前記環境モデルを更新するモデル更新ステップを有し、該モデル更新ステップにおいて、新たに入力された被計測点の座標値に対応するボクセルを探索し、前記ボクセル内の代表点および誤差分布の少なくともいずれかが新たに設定される、または再設定される、または該ボクセルを更に分割して階層的に複数のボクセルに分割される場合のみ、前記出力ステップにおいて、当該ボクセルの代表点の位置、誤差分布およびボクセルの位置の少なくともいずれかを三次元形状の計測値として出力装置に出力することによって、距離センサから新たに得られた被計測点によって影響を受けたボクセル内の代表点等の値が出力される。このため、使用者は従来と同様の運用を想定しながら、距離センサで得られる元々の計測値があたかもより精度の高い計測値に置き換わったように、利用することができるようになる。このように、より精度の高い三次元形状計測が可能になる。
【００７３】
また、前記出力ステップにおいて、距離センサの位置から距離センサが位置計測可能な範囲の環境モデル内のボクセルの代表点の位置を三次元形状の計測値として出力装置に出力することによって、従来の距離センサの計測値の分解能が粗い場合でも、あたかも精度良くかつ分解能が高い距離センサのように利用することが可能となる。このように、より精度の高い三次元形状計測が可能になる。  
本発明のその他の目的及び有利な特徴は、添付図面を参照した以下の説明から明らかになろう。
【図面の簡単な説明】
【００７４】
【図１】特許文献１の「自己位置同定装置」の模式図である。
【図２】本発明の自己位置同定装置の全体構成図である。
【図３】非特許文献２に開示された三次元レーザレーダの構成図である。
【図４Ａ】距離センサで計測された極座標データと誤差の関係を示す図である。
【図４Ｂ】誤差分布を直方体に包含される楕円体として近似する場合を示している。
【図５】本発明の方法を実行するための装置構成図である。
【図６】本発明の方法を示すフローチャートである。
【図７】モデル構築ステップの模式図である。
【図８】構築された環境モデルの模式図である。
【図９Ａ】本発明におけるボクセルデータのデータ構造を示す図であり、各ボクセルデータのメモリレイアウト例を示している。
【図９Ｂ】本発明におけるボクセルデータのデータ構造を示す図であり、レベル２（１，１，０）のボクセルが代表点を持つ場合の例を示している。
【図１０】粗合わせステップＳ６と精密合わせステップＳ７のデータ処理フロー図である。
【図１１】粗合わせステップＳ６の模式図である。
【図１２】精密合わせステップＳ７の模式図である。
【図１３】モデル更新ステップにおけるデータ処理フロー図である。
【図１４】該当するボクセル内に既に設定した代表点がある場合の模式図である。
【図１５】誤差分布が互いに重複する場合に、両誤差分布から新たな誤差分布と新たな誤差分布の中心を再設定した結果、新たな代表点が他のボクセル内へ移動する場合を示している。
【図１６】誤差分布が互いに重複する場合の模式図である。
【図１７】複数の計測位置を原点とする距離データを統合することによって、代表点の誤差分布が縮小し、代表点の精度が向上する様子を示す模式図である。
【図１８】カルマンフィルタを用いたモデル更新ステップにより得られた結果を示す。
【図１９】図１８の一部拡大図である。
【図２０】誤差を考慮した対応付けを示している。
【好ましい実施例の説明】
【００７５】
以下本発明の好ましい実施形態について、図面を参照して説明する。なお、各図において、共通する部分には同一の符号を付し、重複した説明を省略する。
【００７６】
図２は、本発明の自己位置同定装置の全体構成図である。この図において、本発明の自己位置同定装置は、６自由度同定手段であり、距離センサ、オドメータ、カメラ、ＧＰＳ、姿勢センサから外界の三次元情報を取り入れて、自己位置を同定する装置である。なお、本発明において、自己位置とは、自己位置同定装置自体の外界における６自由度の位置と姿勢を意味する。
なお、本発明において、距離センサ以外のオドメータ、カメラ、ＧＰＳ、姿勢センサは不可欠ではなく、必要に応じてオプションとして用いる。以下、距離センサを用いた例を説明する。
【００７７】
図３は、距離センサの一例としての三次元レーザレーダの構成図である。三次元レーザレーダは、例えば非特許文献３に開示されている。
この図に示すように、三次元レーザレーダ１０は、レーダヘッド１２と制御器２０から構成される。レーザダイオード１３から発振されたパルスレーザ光１は、投光レンズ１４で平行光２に整形され、ミラー１８ａ，１８ｂと回転・揺動するポリゴンミラー１５で二次元方向に走査され、測定対象物に照射される。測定対象物から反射されたパルスレーザ光３は、ポリゴンミラー１５を介して受光レンズ１６で集光され、光検出器１７で電気信号に変換される。
【００７８】
制御器２０内の時間間隔カウンタ２１は、レーザダイオード１３のパルス発振タイミングと同期したスタートパルス４と、光検出器１７から出力されたストップパルス５の時間間隔を計測する。信号処理ボード２２は、反射光が検出された時点の時間間隔ｔ、ポリゴンミラーの回転角度θ、揺動角度φを極座標データ（ｒ，θ，φ）として出力する。
ｒは計測位置（レーダヘッド設置位置）を原点とする距離であり、ｒ＝ｃ×ｔ／２  の式で求められる。ここでｃは光速である。
判定処理ユニット２３は、信号処理ボードからの極座標データを、レーダヘッド設置位置を原点とした三次元空間データ（ｘ，ｙ，ｚ）へ変換して、検出処理を行うようになっている。なおこの図で２４はドライブユニットである。
【００７９】
上述した三次元レーザレーダ１０の計測範囲は、例えば、水平画角６０°、垂直画角３０°、最大測定距離５０ｍである。また、位置検出精度は、例えば約２０ｃｍである。
また、計測データを各画素に対して奥行き方向の距離値を持った距離画像で表示する場合、１フレームの計測点数を、横方向１６６点、スキャン方向５０点とすると、１フレームに１６６×５０＝８３００点が表示される。この場合にフレームレートは、例えば約２フレーム／秒である。
【００８０】
この三次元レーザレーダ１０で計測される三次元形状上の被計測点は、横方向にΔθ×ｒ、垂直方向にΔφ×ｒ、互いに離散した点群となる。例えば、Δθ＝６０／１６６×π／１８０＝６．３×１０−３ラジアン、Δφ＝３０／５０×π／１８０＝１０．５×１０−３ラジアン、ｒ＝５０ｍの場合、最も近接する場合でも、被計測点の間隔は、横方向で約３１５ｍｍ、垂直方向で約５２５ｍｍとなる。
【００８１】
本発明では、距離センサとして、例えば、上述した三次元レーザレーダ１０を用いる。しかし、距離センサはこれに限定されず、視差を利用した距離センサ、その他の周知の距離センサを用いることができる。
【００８２】
図４Ａ，図４Ｂは、距離センサで計測された極座標データと誤差の関係を示す図である。
図４Ａに示すように、任意の計測位置を原点とする極座標値（ｒ，θ，φ）を計測結果として計測する。距離センサによる計測結果には、図に示すような誤差分布が通常存在する。
この誤差分布は、誤差分布のｒｓ，θｓ，φｓでの存在確率をＰ（ｒｓ，θｓ，φｓ）とした場合、誤差分布は計測の軸ｒ，θ，φ方向に正規分布しているとし、例えば式（１）で表すことができる。ここで、ｒ，θ，φはセンサからの計測値、σｒ，σθ，σφは  は標準偏差、Ａは規格化定数である。
図４Ｂに示すように、誤差分布は、通常ｒ方向に長い切頭円錐形（左図）に内包される分布であるが、遠方においてａとｂの差は小さい。従って、この誤差分布を直方体に包含される楕円体として安全サイドに近似することができる。
【００８３】
【数２】
0004650752.tif 000003
【００８４】
図５は、本発明の方法を実行するための装置構成図である。この図に示すように、この装置は、データ入力装置３２、外部記憶装置３３、内部記憶装置３４、中央処理装置３５および出力装置３６を備える。
【００８５】
データ入力装置３２は、上述した距離センサを有し、三次元形状上の座標値を記憶装置に入力する。また、例えばゴニオメータ、オドメータ等を併用して、距離センサの位置・姿勢や移動距離も入力するのがよい。なお、データ入力装置３２は、キーボード等の通常の入力手段も有するのがよい。
外部記憶装置３３は、ハードディスク、フロッピー（登録商標）ディスク、磁気テープ、コンパクトディスク等である。外部記憶装置３３は、環境モデルのサイズが大きく後述する内部記憶装置３４に入力された三次元形状上の座標値、ボクセル位置、及び代表点とその誤差分布の全体を保持できない場合には、環境モデルの一部範囲または全体範囲に対する、入力された三次元形状上の座標値、ボクセル位置、及び代表点とその誤差分布の一部または全体を記憶し、かつ本発明の方法を実行するためのプログラムを記憶する。
内部記憶装置３４は、例えばＲＡＭ，ＲＯＭ等であり、環境モデルの一部範囲または全体範囲に対する、入力された三次元形状上の座標値、ボクセル位置、及び代表点とその誤差分布の一部または全体を保管し、かつ演算情報を保管する。
中央処理装置３５（ＣＰＵ）は、モデル構築装置、マッチング装置、粗合わせと精密合わせの位置合わせ装置、モデル更新装置、データ伝達装置として機能し、演算や入出力等を集中的に処理し、内部記憶装置３４と共に、プログラムを実行する。モデル構築装置は、後述のモデル構築ステップを行う装置であり、マッチング装置は、後述のマッチングステップを行う装置であり、位置合わせ装置は、後述の粗合わステップと精密合わせステップを行う装置であり、モデル更新装置は、後述のモデル更新ステップを行う装置であり、データ伝達装置は、出力装置３６へデータを出力する装置である。
出力装置３６は、例えば表示装置、プリンタ、外部装置等であり、内部記憶装置３４および外部記憶装置３３の少なくともいずれかに記憶したデータ及びプログラムの実行結果を出力するようになっている。外部装置とのインターフェイスは、ＬＡＮ、ＵＳＢ、ＩＥＥＥ１３９４等であり、入力された三次元形状上の座標値に対して該当するボクセル内の代表点、誤差分布、ボクセル位置などを付加した結果や、環境モデル全体または環境モデル一部を要求に応じて出力する。
【００８６】
上述した本発明の装置は、上述した距離センサと通常のＰＣ（コンピュータ）を組み合わせたものでもよく、或いは、全体を一体にした装置であってもよい。また、自走可能な装置内に一体的に組み込んでもよい。
【００８７】
図６は、本発明の方法を示すフローチャートである。
本発明の方法は、外界の三次元形状を取り入れて自己位置を同定する自己位置同定方法、および三次元形状の計測方法であり、データ入力ステップＳ１、データ補正ステップＳ２、探索範囲限定ステップＳ３、モデル構築ステップＳ４、マッチングステップＳ５、自己位置同定ステップＳ７、Ｓ１０、粗合わせステップＳ８、精密合わせステップＳ９、モデル更新ステップＳ１１及び出力ステップＳ１２を有する。
なお、これら一連の処理のうち、Ｓ１、Ｓ２、Ｓ３，Ｓ５〜Ｓ１２は、計測データが得られる毎に実施し、Ｓ４は初めて計測データが得られたときにだけ実施する。
【００８８】
データ入力ステップＳ１では、新たな計測位置において、距離センサを用いて、三次元形状上の座標値をコンピュータの記憶装置に入力する。また、例えばゴニオメータ、オドメータ等を併用して、距離センサの位置・姿勢や移動距離も入力するのがよい。
なおこのデータ入力ステップＳ１において、三次元レーザレーダ１０を用いて、三次元形状上の座標値を任意の計測位置を原点とする距離データとして原点を移動しながら順次取得するのがよい。
【００８９】
距離センサとして三次元レーザレーダ１０を用いた場合、三次元形状上の座標値は、任意の計測位置を原点とする距離データであり、極座標値（ｒ，θ，φ）で表される。また、各座標値の誤差分布は、極座標値（ｒ，θ，φ）から演算で求めるか、予め別の入力手段（例えばキーボード）で入力する。
【００９０】
データ補正ステップＳ２では、距離データの補正処理を行い、距離データの精度を向上させる。また、極座標データとオドメータのデータから、任意の固定位置を原点とした三次元空間データ（ｘ，ｙ，ｚ）へ変換してもよい。
距離データの補正処理では、孤立点の除去、統計的処理、等を行う。孤立点は、周囲の点から孤立して存在する点であり、計測データは複数の近接する点で構成されることから、孤立点は誤計測と仮定して除去することができる。統計的処理は、計測データが含む誤差分布を考慮して、複数回の計測を統計処理（例えば平均値等）することで、距離の補正を行う。
さらに、対象とする三次元形状が、直線近似又は平面近似できる場合にはこれらを行うのがよい。
【００９１】
探索範囲限定ステップＳ３では、距離センサの探索範囲を限定する。
探索範囲を限定せずに環境モデルに対する計測データのマッチング処理を行うと、複数の解（被計測点）が得られる可能性がある。そこで、（１）現在のセンサ位置を過去のセンサ位置の変化から推定し、センサ位置推定結果の近傍を探索する、（２）オドメータを用いてセンサ位置を推定し、探索範囲を限定する、（３）距離データのうち、距離値だけでなく、反射強度値も利用して探索結果を絞り込む、等を実施する。
【００９２】
図７は、ボクセルの分割に八分木を用いた場合のモデル構築ステップの模式図である。
モデル構築ステップＳ４では、この図に示すように、三次元形状の存在する空間領域を、境界表面が互いに直交する直方体からなる複数のボクセル６に分割し、各ボクセル位置を記憶する環境モデルを構築する。
ボクセル６の形状は、各辺の長さが等しい立方体でも、各辺の長さが異なる直方体でもよい。
また、ボクセル６の各辺の長さは、最大のボクセル６を必要最小限の分解能に相当する大きさに設定するのがよい。以下、最大のボクセル６をレベル１のボクセルと呼ぶ。
また、単一のボクセル内に複数の被計測点が存在する場合には、単一のボクセル内に単一の被計測点のみが存在するように、例えば八分木を選んだ場合は、ボクセルを更に八分割して階層的に複数のボクセルに分割する。以下、最大のボクセル６の八分割を１回実施した空間領域をレベル２のボクセル、ｋ回実施した空間領域をレベルｋ＋１のボクセルと呼ぶ。
【００９３】
図８は、構築された環境モデルの模式図である。
マッチングステップＳ５では、この図に示すように、三次元形状上の座標値に対応するボクセル６の内部に代表点７とその誤差分布８を設定し記憶する。末端のボクセルは計測値の代表点を１つだけ持つことができる。各ボクセルが計測値の代表点とその誤差分布を持つことで、物体の形状を表す。また、物体の存在確率を表す確率値をボクセルに持たせることもできる。
【００９４】
マッチングステップＳ５において、代表点の絶対位置は次の［数３］の式（２）で与えられる。ここで、（ｘ，ｙ，ｚ）は代表点のボクセルでの相対座標、Ｓｘ，Ｓｙ，Ｓｚはレベル１でのボクセルの一辺の大きさ、ｎｘ（ｋ），ｎｙ（ｋ），ｎｚ（ｋ）はレベルｋでのボクセルの番地、Ｌは求める代表点が存在するレベルである。
【００９５】
【数３】
0004650752.tif 000004
【００９６】
図９Ａ，図９Ｂは、本発明におけるボクセルデータのデータ構造を示す図である。
この図において、図９Ａは、各ボクセルデータのメモリレイアウト例である。この図において、矢印はデータへのリンクを表し、値としてはデータへのポインタを保持する。
図９Ｂは、レベル２（１，１，０）のボクセルが代表点を持つ場合の例を示している。なおこの図において、ｎｕｌｌは空集合を表す。
【００９７】
上述したデータ構造の環境モデルは、以下の特徴を有する。
（１）内容：空間を小直方体で分割して各ボクセルに計測点の代表点と誤差分布を保持する。
（２）精度：ボクセル毎に持つ計測点の代表値相当である。
（３）存在：物体の存在の有無を表現できる。
（４）データ量：ボクセルの個数に比例してメモリを必要とするが、サイズ固定である。
（５）点群からの変換：適しており、計算量は少ない。
（６）アクセス速度：シンプルな構造をしているため、要素へのアクセスが高速である。
【００９８】
またこの特徴から、上述した環境モデルは、以下の効果Ａ〜Ｃをすべて満たしている。
効果Ａ：誤差を考慮した表現が可能である。
効果Ｂ：必要なメモリ量と計算量が一定量以下である。
効果Ｃ：物体の存在だけでなく、物体が存在しないことを表せる。
【００９９】
図６において、マッチングステップＳ５の後に、ステップＳ６で先の計測位置おける同一の三次元形状のデータの有無をチェックする。このチェックで、先の計測位置おける同一の三次元形状のデータが存在しない場合には、新たな計測位置を自己位置と同定する。
この同定は、逐次移動する移動体の初期位置において、既知の６自由度位置（例えばグローバル座標系の原点）で実施するのが好ましい。また、この同定において、新たな計測位置の位置（３自由度）と姿勢（３自由度）を含む６自由度位置を同定するのがよい。
【０１００】
また、ステップＳ６のチェックで、先の計測位置おける同一の三次元形状のデータが存在する場合には、図６において、粗合わせステップＳ８と精密合わせステップＳ９を、マッチングステップＳ５の後に実施する。
図１０は、粗合わせステップＳ８と精密合わせステップＳ９のデータ処理フロー図であり、図１１は粗合わせステップＳ８の模式図、図１２は精密合わせステップＳ９の模式図である。
【０１０１】
図１０において、粗合わせステップＳ８では、図１１に示すように、先の計測位置に対する環境モデルに対し、
新たな計測データおよび誤差分布を回転及び並進させ、計測データおよび誤差分布と近接する代表点を有するボクセル間との距離に関する評価値（例えば、当該距離の総和）が最小になるように、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、代表点を有するボクセル間の距離に関する評価値（例えば、当該距離の総和）が最小になるように、
もしくは、新たな計測データおよび誤差分布を回転及び並進させ、計測データおよび誤差分布と近接する代表点を有するボクセルの確率値に関する評価値（例えば、当該確率値の総和）が最大になるように、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、近接するボクセルが持つ確率値の差に関する評価値（例えば、当該確率値の差の総和）が最小になるように、位置合わせする。
【０１０２】
粗合わせステップＳ６における位置合わせは、環境モデルと計測データをボクセル空間上で表現すること、もしくは環境モデルはボクセル空間上で、計測データは代表点と誤差分布表現することで実施する。現在の計測データが位置（ｘ，ｙ，ｚ）、姿勢（θ，φ，ψ）での計測であったとして、計測データをワールド座標に変換して環境モデルとの一致度を算出する。
一致度の算出には、例えば最短距離法を用いることができる。最短距離法を用いた場合のボクセル間の距離は、２つのボクセル空間をｘ（１）、ｘ（２）、ボクセルの総数Ｉを、ボクセルの値ｘｉ（ｎ）とすると、［数４］の式（３）で定義できる。
計測データの最適な位置・姿勢は、位置（ｘ，ｙ，ｚ）、姿勢（θ，φ，ψ）を変化させることによってεを最小にする最小自乗法によって算出できる。
【０１０３】
また、一致度として、例えば環境モデルと計測データの両ボクセルにおいて、近接する両ボクセルの持つ確率値の差に関する評価値（例えば、当該確率値の差の総和）を用いることができる。この場合は一致度を最小にするように、計測データの最適な位置・姿勢を変化させる。
また、環境モデルはボクセル空間上で、計測データは代表値と誤差分布表現した場合には、計測データの代表値、および誤差分布が近接する環境モデルのボクセルの確率値に関する評価値（例えば、当該確率値の総和）を用いることができる。この場合は一致度を最大にするように、計測データの最適な位置・姿勢を変化させる。
【０１０４】
【数４】
0004650752.tif 000005
【０１０５】
図１０において、精密合わせステップＳ９では、図１２に示すように、先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、近接する誤差分布間の距離に関する評価値（例えば、当該距離の総和）が最小になるように位置合わせする。
精密合わせステップＳ７における環境モデルと計測データの精密合わせ位置合わせには、点群と点群の位置合わせが可能なＩＣＰアルゴリズムに誤差分布を考慮した手法を利用する。位置合わせの初期値には、粗い位置合わせにより得られた位置・姿勢を利用する。
ＩＣＰアルゴリズムに利用する誤差分布間の距離の算出には、例えば誤差分布が交差する場合を同一計測点と考え、その場合の距離値に分布の一致度から求めた重みを掛け合わせて算出する。分布の一致には例えばマハラノビス距離のような距離尺度を利用できる。
この場合の環境モデルと計測データの距離は、環境モデルデータをｐＭｉ、環境モデルデータの誤差分布をΣＭｉ、計測データをＰＤｉ、計測データの誤差分布をΣＤｉ、誤差分布の合成関数をｗ、計測データに対応する環境モデルデータの個数をＮとすると、［数５］の式（４）で定義できる。ここで、Ｔは転置を表す。
計測データの最適な位置・姿勢は、計測データを計測した位置（ｘ，ｙ，ｚ）  、姿勢（θ，φ，ψ）を変化させてＰＤｉを移動することによりεを最小にする最小自乗法によって算出できる。
また、自己位置・姿勢の同定とともに、自己位置の誤差分布を同定し、現在の自己位置と誤差分布および同定した自己位置と誤差分布とからカルマンフィルタによって、自己位置と誤差分布を補正する。
【０１０６】
【数５】
0004650752.tif 000006
【０１０７】
さらに図６において、モデル更新ステップＳ１１は、自己位置同定ステップＳ１０の後に実施し、モデル構築ステップＳ４で構築した環境モデルを更新する。
【０１０８】
図１３は、モデル更新ステップＳ１１におけるデータ処理フロー図である。この図に示すように、ステップＳＴ１で新たに入力された被計測点の座標値に対応するボクセルを探索し、ステップＳＴ２で該当するボクセル内に代表点がない（ボクセルが空である）場合には、ステップＳＴ３で新たに入力された被計測点の座標値と誤差分布を代表点の座標値と誤差分布として設定（新規に登録）する。
また、このステップＳＴ３において、新しい計測位置（原点）と被計測点の間には、原理的に物体が存在しないはずである。従って新しい計測位置（原点）と被計測点の間に位置するボクセル内の代表点と誤差分布を再設定、もしくは消去する。
【０１０９】
図１４は、該当するボクセル内に既に設定した代表点がある場合の模式図である。
図１３のステップＳＴ２で該当するボクセル内に既に設定した代表点がある場合には、ステップＳＴ４で新たに取得した誤差分布と既に設定したボクセル内の誤差分布とを比較する（すなわち異なる点か同一点かを判断する）。
この比較で、誤差分布が互いに重複する場合（図１４の（Ａ））には、ステップＳＴ５で両誤差分布から、または、両誤差分布とボクセル内に既に設定した代表点と新たに入力された被計測点の座標値から、新たな誤差分布と新たな代表点を再設定する（すなわち誤差分布を合成する）。
またこの比較で、誤差分布が互いに重複しない場合（図１４の（Ｂ））には、ステップＳＴ６、ＳＴ７で単一のボクセル内に単一の代表点のみが存在するように、該ボクセルを更に八分割して階層的に複数のボクセルに分割し新規に登録する。
分割と合成の基準は、例えば誤差分布の一致度から判断する。誤差分布の一致度には例えば、マハラノビス距離のような距離尺度を利用できる。また、２つの誤差分布に基づき、両者が同一点を表しているかを統計的検定によって判定してもよい。
【０１１０】
ステップＳＴ５で両誤差分布から新たな誤差分布と新たな誤差分布の中心を再設定した結果、新たな代表点が他のボクセル内へ移動したとき（即ち、ステップＳＴ８で、Ｙｅｓ）、ステップＳＴ２へ戻り、上述の処理を繰り返す。
なお、図１５は、ステップＳＴ５で両誤差分布から、または、両誤差分布とボクセル内に既に設定した代表点と新たに入力された被計測点の座標値から、新たな誤差分布と新たな誤差分布の中心を再設定した結果、新たな代表点が他のボクセル内へ移動する場合を示している。
【０１１１】
ボクセル内に物体の存在確率を表す確率値を設定する場合には、モデル更新ステップＳ１１において、ボクセル内の代表点と誤差分布の新規登録、もしくは再設定、もしくは消去、もしくは分割後に新規登録の処理に応じて、ボクセル内の確率値も統計処理によって、新規登録、もしくは再設定、もしくは消去、もしくは分割後に新規登録を行う。
【０１１２】
図１６は、誤差分布が互いに重複する場合（図１４の（Ａ））の別の模式図である。ステップＳＴ５において、２つの代表点と誤差分布を合成して新たな代表点を誤差分布を設定する手段として、カルマンフィルタを用いることができる。例えば、二次元の場合に、この図に示すように、２つの代表点をそれぞれｘ（１），ｘ’（２）、２つの誤差分布をΣ（１）、Σ’（２）とし、これを合成した代表点をｘ（２）、誤差分布をΣ（２）とすると、代表点ｘ（２）と誤差分布Σ（２）を算出する模式図は図１６のようになる。
【０１１３】
図６において、出力ステップＳ１２では、前記自己位置、前記自己位置に基づいた前記ボクセル位置、代表点、および、誤差分布の少なくともいずれかを出力装置３６に出力する。出力装置３６が表示装置（例えばＣＲＴ）の場合、三次元画像上に自己位置を立体表示するのが好ましい。また、自己位置のデータを別の装置（例えば制御装置、コンピュータ）に転送してもよく、プリンタで出力してもよい。
また、出力ステップＳ１２において、ボクセルの代表点の位置を三次元形状の計測値として出力装置３６に出力するとともに、該計測値の信頼性または精度を示す指標（例えば、数値）を、該ボクセルの内部の誤差分布の大きさに基づいて、出力装置３６に出力してもよい。さらに、出力ステップＳ１２において、ボクセルの代表点の位置を三次元形状の計測値として出力装置３６に出力するときに、該ボクセルの内部の誤差分布の大きさ（広がり）が所定の基準値よりも大きい場合に、該計測値の信頼性または精度が所定の基準よりも低いとして、該ボクセルの前記計測値（即ち、このボクセルの代表点の位置）を出力装置３６に出力しないようにしてもよい。
【０１１４】
また、前記マッチングステップの後に、前記環境モデルを更新するモデル更新ステップを有し、該モデル更新ステップにおいて、新たに入力された被計測点の座標値に対応するボクセルを探索し、前記ボクセル内の代表点および誤差分布の少なくともいずれかが新たに設定される、または再設定される、または該ボクセルを更に分割して階層的に複数のボクセルに分割される場合、前記出力ステップにおいて、当該ボクセルの代表点の位置を三次元形状の計測値として出力装置に出力しても良い。
【０１１５】
また、出力ステップＳ６において、距離センサの位置・姿勢が得られる場合は、その位置から見える範囲の環境モデル内のボクセルの代表点の位置を三次元形状の計測値として出力装置に出力してもよい。距離センサの位置から見える範囲とは、距離センサの位置から距離センサが位置計測可能な範囲であり、距離センサの位置から距離センサが位置計測可能な角度範囲（視野）と、距離センサの位置から距離センサが位置計測可能な距離範囲とを含んでよい。
【０１１６】
図６に示した処理の手順は、新たな計測位置において、新しい計測データが得られる度に、処理を繰り返し行い、内部記憶装置３４および外部記憶装置３３の少なくともいずれかに結果を格納する。処理を高速化するためには、内部記憶装置３４に容量が許す限り結果を格納することが好ましい。
【０１１７】
上述した本発明の方法と装置によれば、三次元形状の存在する空間領域を、複数のボクセル６に分割し、各ボクセル位置を内部記憶装置３４および外部記憶装置３３の少なくともいずれかに記憶するので、計測対象物が大きい場合であっても、データ量をボクセル数に比例する小さいデータサイズに抑えることができる。
また、座標値に対応するボクセル６の内部に代表点７とその誤差分布８を設定し記憶するので、ボクセルの分解能以上の情報を表現することができる。
従って、本発明のデータ構造により、複数視点の計測データを固定サイズに統合できる。
【０１１８】
また、ボクセルの内部に物体の存在確率を表す確率値を設定し、保存することによって、誤差分布が代表点の属するボクセルよりも広がっている場合においても、各ボクセルにおける物体の存在有無を代表点が属するボクセルを見つけ、その誤差分布から再計算させることなく、当該ボクセルの確率値だけで容易に判断できるので、検索時間を抑えることができる。
また、現在の計測位置を過去の計測位置の変化から推定し、もしくは現在の計測位置を取得可能なセンサにより取得し、もしくは計測データの距離値だけでなく反射強度を利用し、照合する範囲を制限するので、検索時間を抑えることができる。
【０１１９】
また、粗合わせステップＳ６において、先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、計測データおよび誤差分布と近接する代表点を有するボクセル間との距離に関する評価値（例えば、当該距離の総和）が最小になるように、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、代表点を有するボクセル間の距離に関する評価値（例えば、当該距離の総和）が最小になるように、
もしくは新たな計測データおよび誤差分布を回転及び並進させ、計測データおよび誤差分布と近接する代表点を有するボクセルの確率値に関する評価値（例えば、当該確率値の総和）が最大になるように、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、近接するボクセルが持つ確率値の差に関する評価値（例えば、当該確率値の差の総和）が最小になるように、位置合わせするので、誤差の蓄積を防ぎながら、短時間に代表点を有するボクセル同士の位置合わせができる。
【０１２０】
また、先の計測位置に対する環境モデルに対し、新たな計測位置に対する環境モデルを回転及び並進させ、近接するボクセルが持つ確率値の差に関する評価値（例えば、当該確率値の差の総和）が最小になるように、位置合わせする場合には、物体が存在しない情報も加味して位置合わせするので、精度の向上が図れる。
【０１２１】
次いで、精密合わせステップＳ７において、先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、もしくは新たな計測位置に対する環境モデルを回転及び並進させ、近接する誤差分布間の距離に関する評価値（例えば、当該距離の総和）が最小になるように位置合わせするので、誤差分布同士の精密な位置合わせを短時間にできる。
従って、本発明による複数視点からのデータの統合処理により、誤差の蓄積を防ぎながら、高精度な形状取得を可能にできる。
【０１２２】
また、本発明のデータ構造は、ボクセル構造を拡張した構造をしているため、データサイズは点群に比べて小さく済ませることができる。従って、周囲環境と計測データの位置合わせ（ＩＣＰアルゴリズム）を行う際に探索対象となるデータ量が減少することから、計算を効率化できる。
【０１２３】
さらに、従来のＩＣＰアルゴリズムは疎データに対して誤った結果を出力するが、本発明の環境モデルは、ボクセル内に代表点と誤差分布を持っているため、疎データに対応した位置合わせが可能である。
【０１２４】
また、モデル構築ステップＳ４において、最大のボクセル９を必要最小限の分解能に相当する大きさに設定し、かつ単一のボクセル９内に複数の被計測点が存在する場合に、単一のボクセル内に単一の被計測点のみが存在するように、該ボクセルを更に八分割して階層的に複数のボクセルに分割するので、データ量を小さいデータサイズに抑えると同時に、分割後のボクセルと代表点を用いて解像度を更に高めることができる。
【０１２５】
特に、三次元形状上の複数の座標値を複数の計測位置を原点とする距離データとして取得し、該距離データの座標値を、前記代表点の座標値とし、距離データの座標値の計測誤差を代表点の誤差分布とすることにより、正確な座標値と誤差分布を用いて複数回の計測を統計的に統合することができ、一層の精度向上が可能となる。
【０１２６】
図１７は、複数の計測位置を原点とする距離データを統合することによって、代表点の誤差分布が縮小し、代表点の精度が向上する様子を示している。このように異なる計測位置（即ち、距離センサである三次元計測器の位置）を原点として得られた距離データは誤差分布の向きも異なるので、これらの距離データを環境モデルを介して逐次統合することによって、代表点の誤差分布が縮小し、代表点の位置精度が向上する。なお、図１７において、三次元計測後の図はコップの２次元断面を表わした模式図であり、三次元計測後の図の破線はコップの実際表面を表わしている。
【０１２７】
また、自己位置の同定とともに、自己位置の誤差分布を同定し、現在の自己位置と誤差分布および同定した自己位置と誤差分布とからカルマンフィルタによって、自己位置と誤差分布を補正することによって、一層の精度向上が可能となる。
【０１２８】
また、原点と被計測点の間に物体が存在しないものとして、その間に位置するボクセル内の代表点と誤差分布を再設定、もしくは消去することにより、誤った計測データの影響を除去することができる。
【０１２９】
また、新たに入力された被計測点の座標値に対応するボクセルを探索し、該ボクセル内に代表点がない場合に、前記座標値と誤差分布を代表点の座標値と誤差分布として設定することにより、代表点の座標値と誤差分布を容易に設定できる。
【０１３０】
更に、前記ボクセル内に既に設定した代表点がある場合に、新たに取得した誤差分布と既に設定したボクセル内の誤差分布とを比較し、
誤差分布が互いに重複する場合に、両誤差分布から、または、両誤差分布とボクセル内に既に設定した代表点と新たに入力された被計測点の座標値から、新たな誤差分布と新たな代表点を再設定し、
誤差分布が互いに重複しない場合に、単一のボクセル内に単一の代表点のみが存在するように、該ボクセルを更に八分割して階層的に複数のボクセルに分割する、ことにより、誤差の蓄積を回避しながら高精度な形状に収束させることができる。
【０１３１】
従って、本発明の方法と装置によれば、誤差を含むデータを正確な情報に補正する機能を有すると共に、長時間の計測に対して高精度に収束する。なおかつ、位置同定処理は、ボクセルの各ボクセルに対応する代表点とその誤差分布を計測点で更新する処理であるため計算量が小さく、演算は周囲ボクセルへの影響を与えずボクセル内で閉じている。よって、高速な処理が可能である。また、計測データはボクセルに逐次統合可能であり、その結果得られた周囲環境情報のメモリサイズは固定サイズを上回ることはない。
【０１３２】
カルマンフィルタを用いたモデル更新ステップについて、詳しく説明する。
【０１３３】
カルマンフィルタを用いたモデル更新ステップの場合には、新たに入力された被計測点の座標値およびその誤差分布と、既に設定したボクセル内の代表点およびその誤差分布とから、カルマンフィルタにより新たな代表点と誤差分布を取得して再設定する。
【０１３４】
各モデル点群の位置m（ｉ）を状態量とし、距離センサの計測点の位置を基に、モデルを次の［数６］で表現する。なお、本実施例では、m（ｉ）は、ボクセル内部の代表点である（以下、同様）。
【０１３５】
【数６】
0004650752.tif 000007
【０１３６】
［数６］において、
Ｌ（ｊ）は、距離センサによる計測位置である。例えば、Ｌ（ｊ）は、距離センサのセンサ座標系において三次元ＬＲＦ（レーザレンジファインダ）の計測点ｊ（ｊ＝１，．．．，Ｎ）の位置Ｌ（ｊ）＝（ｘＬ（ｊ），ｙＬ（ｊ），ｚＬ（ｊ））ｔである。ここで、ｔは転置行列を示す（以下、同様）。
ｈｍ（Ｒｒ，ｔｒ，ｍ（ｉ））は、Ｌ（ｊ）に対する観測系モデルである。
Ｒｒは、距離センサを搭載した移動体（例えば移動ロボット）のワールド座標系に対する姿勢を表す回転行列Ｒｒ＝Ｒ（θｘ，θｙ，θｚ）である。なお、θｘ，θｙ，θｚは、それぞれｘ軸、ｙ軸、ｚ軸周りの回転角を示す（以下、同様）。
ｔｒは、上記移動体のワールド座標系に対する位置を表す並進ベクトルｔｒ＝（ｘ，ｙ，ｚ）である。
ｖＬ（ｉ）は、距離センサの計測値Ｌ（ｊ）に加わる観測ノイズである。
Ｒｓは、センサ座標系の移動体座標系に対する回転行列Ｒｓ＝Ｒ（θｘ，θｙ，θｚ）である。
ｔｓは、センサ座標系の移動体座標系に対する位置を表す併進ベクトルｔｓ＝（ｘ，ｙ，ｚ）である。
【０１３７】
測定対象物は静止しているものであり、測定対象物の位置ｔｒ、姿勢Ｒｒを環境モデルに対して固定する。
【０１３８】
距離センサによる計測点群と、環境モデル点群上の点ｉ（即ち、代表点）を対応づける。この対応付けが行われたモデル点群上の点iに対して次式（４）により更新を行う。なお、距離センサによる計測点群と対応付けが行われたモデル点群上の代表点ｍ（i）に対してのみ次の数［７］により更新を行ってよい。
【０１３９】
【数７】
0004650752.tif 000008
【０１４０】
［数７］において、
添え字ｋは、離散時刻ｋでの値であることを表す。
ｍｋ（ｉ）について、ｍ’ｋ（ｉ）はｍｋ（ｉ）の更新値（事後推定値）を示し、ｍｋ，ｋ−１（ｉ）はｍ’ｋ−１（ｉ）に基づいたｍｋ（ｉ）の予測値（事前推定値）を示す。なお、環境（測定対象物）は静止しているので、ｍｋ，ｋ-１（ｉ）＝ｍ’ｋ-１（ｉ）である。
Σmk（ｉ）は、ボクセル内部の代表点ｍｋ（ｉ）の誤差共分散行列(即ち、上述の誤差分布)である。また、Σｍｋ（ｉ）について、  Σ’ｍｋ（ｉ）はΣｍｋ（ｉ）の更新値（事後推定値）を示し、Σｍｋ，ｋ−１（ｉ）はΣ’ｍｋ−１（ｉ）に基づいたΣｍｋ（ｉ）の予測値（事前推定値）を示す。センサ座標系において三次元LRFの計測点j（j=１,…,Ｎ）の位置をＬ（ｊ）で表し、その誤差共分散行列をΣＬ（ｊ）で表す。ここでＮは、三次元LRFで得られた計測点の総数である。三次元LRFの誤差モデルとして計測距離に関係ない一定の正規分布を仮定する。センサ座標系のｘ軸方向にレーザを照射する場合の誤差共分散行列をΣＳとする。レーザの照射方向に応じて誤差分布も姿勢を変える。ΣＬ（ｊ）は、基準の方向に対するレーザ照射方向を回転行列ＲＬ（ｊ）を用いてΣＬ（ｊ）＝ＲＬ（ｊ）ΣＳＲＬｔ（ｊ）と表される。計測点jのワールド座標系における位置ｚ（ｊ）、およびその誤差共分散行列Σｚ（ｊ）は、それぞれｚ（ｊ）＝Ｒｒ（ＲｓＬ（ｊ）＋ｔｓ）＋ｔｒ、Σｚ（ｊ）＝ＲｒＲｓΣＬ（ｊ）ＲｓｔＲｒｔと表すことができる。
Ｋｍｋ（ｉ） は、 ｍｋ（ｉ）に対するカルマンゲインである。
ｈｍｋ（Ｒｒｋ，ｔｒｋ，ｍｋ，ｋ−１（ｉ））は、Ｌｋ（ｊ）、ｉ＝ｐｋ（ｊ）に対する観測系モデルである。ｉ＝ｐｋ（ｊ）は、計測点ｊに対応付けられた環境地図（即ち、環境モデル）上の点である。
Ｈｍｋは、Ｌｋ（ｊ）、ｉ＝ｐｋ（ｊ）に対する観測系モデルのヤコビアン行列であり、次の［数８］で表わされる。
【０１４１】
【数８】
0004650752.tif 000009
【０１４２】
カルマンフィルタの更新過程によって、環境地図のモデル点群の各点（ボクセルの代表点）の位置と誤差共分散行列の更新値ｍ’ｋ（ｉ）、Σ’mk（ｉ）が得られた段階で、環境モデルの更新を以下の手順で行う。
（１）これら更新値ｍ’ｋ（ｉ）、Σ’mk（ｉ）を、新たな代表点、誤差分布として再設定する。
（２）上述（１）の結果、代表点の位置が別のボクセル内に移動した場合、移動先のボクセルが代表点を保持していないときは、移動後の代表点とその誤差共分散行列を移動先のボクセルに保持させ、移動元のボクセルからは代表点等を取り除く。移動先のボクセルが既に代表点を保持しているときには、2つの代表点において、これらの両誤差分布が重複するかを判断する（上述のＳＴ４における判断と同様）。その後の処理は、図１３のＳＴ４以降の処理と同じであってよい。
（３）モデル点群上の代表点ｍ（i）と対応付けが行われなかった距離センサによる計測点について、当該計測点が含まれるボクセルが代表点を持たない場合は、計測点とその誤差分布をそのボクセルの代表点と誤差分布として追加し保持する。もし、ボクセル内に既に代表点が存在する場合には、ボクセル内にある対応付けが行われなかった他の複数の計測点を含め、既存の代表点と各計測点とが全て異なるボクセルに含まれるように、ボクセルを分割した上で分割後のボクセルに代表点等を継承させる。
【０１４３】
上述のカルマンフィルタを用いたモデル更新ステップを繰り返すことで、序々にボクセル内の誤差共分散行列（即ち、誤差分布）の範囲が小さくなるとともに、ボクセルが分割され易くなる。ボクセルが分割されることによって、初期ボクセルのサイズ以下の変化も表現することが可能となる。
図１８は、カルマンフィルタを用いたモデル更新ステップにより得られた結果を示す。図１９は図１８の一部拡大図である。これら図において、初期のボクセルの1辺の長さを100cmとし、再分割数を6分割まで許している。対象が存在している領域では、ボクセルの再分割を繰り返した結果、計測対象を精度良く表現している。対象が存在しない領域ではボクセルの再分割は行われず、必要十分なデータ量で環境を表現できることがわかる。また、各ボクセル内の代表点の誤差分布も小さく、環境地図を高精度で表現できている。このように、誤差を含むデータであってもカルマンフィルタの効果により、真値に収束した結果が得られる。さらに、この方法では計測データ数を増加させることによって標準偏差が小さくなり、精度のさらなる向上が期待できる。
【０１４４】
また、測定対象物の位置・姿勢は固定しているため、更新を測定対象物の位置・姿勢と独立して行うことができる。なお、距離センサによる計測点群と対応付けが行われたモデル点群上の代表点ｍ（i）に対してのみ、上述のカルマンフィルタによる更新を行うことで、大幅な計算コストの削減が可能になる。
【０１４５】
前記精密合わせステップにおいて、前記近接する誤差分布間の距離の総和が最小になるように位置合わせする代わりに、前記近接する誤差分布間に基づく最尤推定値によって定めた一致度に関する評価値（例えば、当該一致度の総乗）が最大となるように、先の計測位置に対する環境モデルに対し、新たな計測データおよび誤差分布を回転及び並進させ、もしくは新たな計測位置に対する環境モデルを回転及び並進させて位置合わせしてもよい。
この場合について詳しく説明する。
【０１４６】
環境地図（環境モデル）であるモデル点群、およびセンサの計測点群双方に誤差モデルを考慮しているため、前記一致度に関する評価値（例えば、当該一致度の総乗）の算出式（以下、評価関数と言う）にも双方の誤差を取り入れることが可能である。本実施例の場合、単に最も近い点を対応付けするのではなく、尤度の概念を取り入れることによって、評価関数に環境地図上の誤差分布、および計測点の誤差分布を加味し、現時点での環境地図において、各計測点の出現確率が最大となる場合に評価関数も最大となるように定める。
具体的には、計測点jに対応付けられた環境地図上の点ｉ＝ｐ（ｊ）の位置が、平均値（代表点）ｍ（ｉ）、誤差共分散行列Σｍ（ｉ）の正規分布に従うと仮定した上で、三次元LRFで計測した結果、Ｌ（ｊ）なる計測データが得られる確率値Ｐｒ（Ｌ（ｊ）｜ｍ（ｉ），Σｍ（ｉ））を点iと点jとの評価関数ＥＭ（ｉ，ｊ）とし、その総乗が最大となるように次の［数９］で評価関数を定める。
【０１４７】
【数９】
0004650752.tif 000010
【０１４８】
ただし、ω（ｊ）は、モデル点群の中に計測点ｊと対応付けられる点が存在する場合は１、それ以外の場合は０とする。
ここで、Ｐｒ（Ｌ（ｊ）｜ｑ）を環境地図の点がｑの位置にある場合にＬ（ｊ）なる計測データが得られる確率値を表すものとし、Ｐｒ（ｑ｜ｍ（ｉ），Σｍ（ｉ））を平均値ｍ（ｉ）、誤差共分散行列Σｍ（ｉ）の正規分布に従うと仮定した上で環境地図の点がｑの位置にある確率値を表すものとすると、数［１０］が成り立つ。
【０１４９】
【数１０】
0004650752.tif 000011
【０１５０】
Ｐｒ（ｑ｜ｍ（ｉ），Σｍ（ｉ））は正規分布を仮定すると、次の［数１１］となる。
【０１５１】
【数１１】
0004650752.tif 000012
【０１５２】
一方、Ｐｒ（Ｌ（ｊ）｜ｑ）は、Ｌ（ｊ）をｚ（ｊ）で置き換えて、次の［数１２］で近似することができる。
【数１２】
0004650752.tif 000013
【０１５３】
ここで、ｚｋ（ｊ）は、距離センサを搭載した移動体、例えば移動ロボットの位置ｔｒ、姿勢Ｒｒに依存している。実際は、三次元LRFのセンサ座標系の中心から見たｑの向きと計測点Ｌ（ｊ）の向きとは図２０に示すように異なるため、誤差共分散行列Σｚ（ｊ）もｑの向きに合わせて回転変換する必要があるが、対応付けられた環境地図上の点iから大きく離れたところにあるｑの存在確率は低いため、十分な精度で近似できると考えられる。よって、Ｐｒ（Ｌ（ｊ）｜ｍ（ｉ），Σｍ（ｉ））は、次の［数１３］で表すことができる。
【０１５４】
【数１３】
0004650752.tif 000014
【０１５５】
簡単な計算によって、次の［数１４］を得る。
【０１５６】
【数１４】
0004650752.tif 000015
【０１５７】
ただし、α（ｊ），β（ｊ）は次の［数１５］で表すことができる。
【０１５８】
【数１５】
0004650752.tif 000016
【０１５９】
従って、モデル点群の点iと計測点群の点jとの対応付けの一致度を表す評価関数ＥＭ（ｉ，ｊ）は、平均値ｍ（ｉ）、誤差共分散行列Σｍ（ｉ）＋Σｚ（ｊ）の正規分布において、ｚ（ｊ）が得られる確率値に近似できる。この評価関数を用いることによって、環境地図、計測データ双方の誤差を考慮した対応付けが可能となる。
【０１６０】
計測点と環境地図（即ち、環境モデル）との対応付けについて補足説明をする。上記実施例では、誤差分布を考慮した統計的な評価関数を用いるため、評価関数の値を求めなければ対応点を定めることができない。そこで、環境地図上のモデル点群の中で対応付けする候補を予め絞り込み、その候補の中から評価関数の値を基に対応点を求める。具体的には、以下のように定めることができる。
（１）対象とする計測点jの誤差共分散行列ΣＬ（ｊ）の範囲（例えば標準偏差の3倍の範囲）と交わる最上位のボクセルとそのボクセルに隣接している最上位のボクセルを求め、下層のボクセルも含めこれらのボクセル内に存在する代表点を対応点の候補とする。ボクセルが階層構造となっているため、この候補点の探索には計算コストはほとんどかからない。このとき、候補となる代表点がない場合には、対応点がないものとみなす。隣接するボクセルも候補に加える理由は、ボクセル内の代表点の位置によっては、誤差共分散行列の範囲が隣接するボクセルまではみ出すことがあるからである。
（２）候補となるボクセルの代表点iと誤差共分散行列を用いて、評価関数ＥＭ（ｉ，ｊ）の値を求める。
（３）評価関数ＥＭ（ｉ，ｊ）の値が最も大きい代表点iを対応点とする。ただし、評価関数の値がある閾値未満の場合には、対応点がないものとみなす。
本実施例では、対応付けの評価関数ＥＭ（ｉ，ｊ）として、尤度に基づいた式を採用しており、対応点の有無に関して統計的に明確な判断尺度があるため、対応点が存在しないと考えられる場合においても無理に対応付けを行うようなことはない。なお、対応点がない場合には、対象となる計測点はこれまで未計測の部分に相当する点であると解釈し、環境地図に追加する。
【０１６１】
実施形態として、外界の三次元情報を取り入れて、自己の三次元位置を同定する自己位置同定方法と装置、および三次元形状の計測方法と装置について説明したが、二次元形状を三次元形状の特別な場合として見ることにより、外界の二次元情報を取り入れて、自己の二次元位置を同定する自己位置同定方法と装置、および二次元形状の計測方法と装置の形態としても実施できる。
また、上記出力ステップにおいて、前記ボクセル位置、代表点、および誤差分布のすべてを出力しなくてもよく、例えば、これらすべてが無くても三次元形状が把握できる場合や、これらのうち１つ又は２つが必要な場合などにおいては、前記ボクセル位置、代表点、および誤差分布の少なくともいずれかを出力装置に出力してもよい。
【０１６２】
なお、本発明は、上述した実施形態に限定されず、本発明の要旨を逸脱しない範囲で種々に変更することができることは勿論である。
【図１】
0004650752.tif 000018
【図２】
0004650752.tif 000019
【図３】
0004650752.tif 000020
【図４Ａ】
0004650752.tif 000021
【図４Ｂ】
0004650752.tif 000022
【図５】
0004650752.tif 000023
【図６】
0004650752.tif 000024
【図７】
0004650752.tif 000025
【図８】
0004650752.tif 000026
【図９Ａ】
0004650752.tif 000027
【図９Ｂ】
0004650752.tif 000028
【図１０】
0004650752.tif 000029
【図１１】
0004650752.tif 000030
【図１２】
0004650752.tif 000031
【図１３】
0004650752.tif 000032
【図１４】
0004650752.tif 000033
【図１５】
0004650752.tif 000034
【図１６】
0004650752.tif 000035
【図１７】
0004650752.tif 000036
【図１８】
0004650752.tif 000037
【図１９】
0004650752.tif 000038
【図２０】
0004650752.tif 000039