知財求人 - 知財ポータルサイト「IP Force」
▶ コミサリヤ・ア・レネルジ・アトミク・エ・オ・エネルジ・アルテルナテイブの特許一覧
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】6869230
(24)【登録日】2021年4月15日
(45)【発行日】2021年5月12日
(54)【発明の名称】物体を認識する方法およびシステム
(51)【国際特許分類】
G06N 7/00 20060101AFI20210426BHJP
B60W 40/04 20060101ALI20210426BHJP
【ＦＩ】
G06N7/00 150
B60W40/04
【請求項の数】17
【全頁数】24
(21)【出願番号】特願2018-512239(P2018-512239)
(86)(22)【出願日】2016年9月22日
(65)【公表番号】特表2018-537745(P2018-537745A)
(43)【公表日】2018年12月20日
(86)【国際出願番号】EP2016072530
(87)【国際公開番号】WO2017050890
(87)【国際公開日】20170330
【審査請求日】2019年7月17日
(31)【優先権主張番号】1558919
(32)【優先日】2015年9月22日
(33)【優先権主張国】FR
(73)【特許権者】
【識別番号】311015001
【氏名又は名称】コミサリヤ・ア・レネルジ・アトミク・エ・オ・エネルジ・アルテルナテイブ
(74)【代理人】
【識別番号】110001173
【氏名又は名称】特許業務法人川口國際特許事務所
(72)【発明者】
【氏名】モッタン，ジュリアン
(72)【発明者】
【氏名】プスチニ・パスクアル，ディエゴ
(72)【発明者】
【氏名】ラコトバオ・アンドリアマエファ，ティアナ
【審査官】
金田 孝之
(56)【参考文献】
【文献】
特開２０１２−０４８６４２（ＪＰ，Ａ）
【文献】
特表２００９−５０８１９９（ＪＰ，Ａ）
【文献】
Tiana Rakotovao et al.，Real-time power-efficient integration of multi-sensor occupancy grid on many-core，2015 IEEE International Workshop on Advanced Robotics and its Social Impacts (ARSO) [online]，２０１５年 ６月３０日，[2020年8月12日検索]，ＵＲＬ，https://ieeexplore.ieee.org/document/7428211
(58)【調査した分野】（Int.Cl.，ＤＢ名）
Ｇ０６Ｎ ３／００−９９／００
Ｂ６０Ｗ ４０／０４
(57)【特許請求の範囲】
【請求項1】
コンピュータまたは専用デジタル電子回路（ＭＴＤ１，ＭＴＤ２）により実装された物体（Ｏ）認識方法であって以下のステップ、すなわち
ａ）１個以上のセンサ（Ｃ１．．．ＣＮＣ）から得られる、前記物体の複数の距離測定値（ｚ１．．．ｚＮＣ）を取得するステップと、
ｂ）前記距離測定値の各々に対し、前記センサの環境の離散化された空間表現を与える占有格子（ＧＯ）上の対応センサの逆モデルを適用することにより、物体による前記占有格子のセルの集合の占有確率を決定するステップと、
ｃ）各セルがステップｂで推定された占有確率の融合により計算された占有確率を示す統合された占有格子を構築するステップとを含み、
前記逆センサモデルの各々が、対応する占有格子の各セルおよび各距離測定値に同一の有限集合内で選択された確率クラスを関連付ける離散的なモデルであり、前記確率クラスの各々が整数添え字により識別されること、および前記ステップｃ）において、統合された占有格子の各セルの占有確率が前記ステップｂ）で決定された確率クラスの添え字に対して実行された整数計算により決定されることを特徴とする方法。
ａ）１個以上のセンサ（Ｃ１．．．ＣＮＣ）から得られる、前記物体の複数の距離測定値（ｚ１．．．ｚＮＣ）を取得するステップと、
ｂ）前記距離測定値の各々に対し、前記センサの環境の離散化された空間表現を与える占有格子（ＧＯ）上の対応センサの逆モデルを適用することにより、物体による前記占有格子のセルの集合の占有確率を決定するステップと、
ｃ）各セルがステップｂで推定された占有確率の融合により計算された占有確率を示す統合された占有格子を構築するステップとを含み、
前記逆センサモデルの各々が、対応する占有格子の各セルおよび各距離測定値に同一の有限集合内で選択された確率クラスを関連付ける離散的なモデルであり、前記確率クラスの各々が整数添え字により識別されること、および前記ステップｃ）において、統合された占有格子の各セルの占有確率が前記ステップｂ）で決定された確率クラスの添え字に対して実行された整数計算により決定されることを特徴とする方法。
【請求項2】
確率クラスの前記有限個の集合が１個以上の部分集合の和集合により形成されていて、前記ステップｃ）において、同一の部分集合に属する２個の確率クラスを融合した結果得られるクラスも同様に前記部分集合に属する、請求項１に記載の方法。
【請求項3】
確率クラスの前記有限個の集合が、確率区間［０、１］において、離散化間隔が０〜０．５の間で増大し、次いで０．５〜１の間で減少するように非一様に離散化されている、請求項１または２に記載の方法。
【請求項4】
有限個の確率クラスのＳｋで示す前記集合が、次式により定義される２個の部分集合
の和集合
により形成されていて、添え字ｎが相対整数値を取り、パラメータｋが非ゼロ相対整数であり、且つ前記ステップｃ）において、２個の確率クラスｐｉ，ｐｊ∈Ｓｋの間の融合Ｆ（ｐｉ、ｐｊ）が、添え字「ｉ」、「ｊ」が同一符号を有する場合に次式
・ｉ，ｊ＜０ならばＦ（ｐｉ，ｐｊ）＝ｐｉ＋ｊ−ｋ・ｉ・ｊ
・ｉ，ｊ＞０ならばＦ（ｐｉ，ｐｊ）＝ｐｉ＋ｊ＋ｋ・ｉ・ｊ
を適用することにより計算される、請求項２および３に記載の方法。
【数1】
【数2】
・ｉ，ｊ＜０ならばＦ（ｐｉ，ｐｊ）＝ｐｉ＋ｊ−ｋ・ｉ・ｊ
・ｉ，ｊ＞０ならばＦ（ｐｉ，ｐｊ）＝ｐｉ＋ｊ＋ｋ・ｉ・ｊ
を適用することにより計算される、請求項２および３に記載の方法。
【請求項5】
前記ステップｃ）において、２個の確率クラスｐｉ，ｐｊ∈Ｓｋの間の融合Ｆ（ｐｉ、ｐｊ）が、添え字「ｉ」、「ｊ」が逆符号を有する場合に次式
・ｉ＜０，ｊ＞０且つ｜ｉ｜≧｜ｊ｜ならば、Ｆ（ｐｉ，ｐｊ）＝ｐ（ｉ＋ｊ）÷（１＋ｋ・ｊ）
・ｉ＜０，ｊ＞０且つ｜ｉ｜≦｜ｊ｜ならば、Ｆ（ｐｉ，ｐｊ）＝ｐ（ｉ＋ｊ）÷（１−ｋ・ｉ）
を適用することにより計算され、記号÷は整数除算を表す、請求項４に記載の方法。
・ｉ＜０，ｊ＞０且つ｜ｉ｜≧｜ｊ｜ならば、Ｆ（ｐｉ，ｐｊ）＝ｐ（ｉ＋ｊ）÷（１＋ｋ・ｊ）
・ｉ＜０，ｊ＞０且つ｜ｉ｜≦｜ｊ｜ならば、Ｆ（ｐｉ，ｐｊ）＝ｐ（ｉ＋ｊ）÷（１−ｋ・ｉ）
を適用することにより計算され、記号÷は整数除算を表す、請求項４に記載の方法。
【請求項6】
ｋ＝１である、請求項４または５に記載の方法。
【請求項7】
有限個の確率クラスのＧｐで示す前記集合が、次式で定義される２個の部分集合
の和集合
により形成され、添え字ｎが相対整数値を取り、ｐが厳密に０．５〜１の間に値を有するパラメータであり、確率ｐｎが次式
ｐ０＝０．５
ｐ１＝ｐ
ｐｎ＋１＝Ｆ（ｐｎ，ｐ）∀ｎ＞１
ｐ−１＝１−ｐ
ｐｎ−１＝Ｆ（ｐｎ，ｐ−１）∀ｎ＜−１
のように再帰的に定義され、且つ前記ステップｃ）において、２個の確率クラスｐｉ，ｐｊ∈Ｇｐの間の融合Ｆ（ｐｉ、ｐｊ）が次式
Ｆ（ｐｉ，ｐｊ）＝ｐｉ＋ｊ
を適用することにより計算される、請求項３に記載の方法。
【数3】
【数4】
ｐ０＝０．５
ｐ１＝ｐ
ｐｎ＋１＝Ｆ（ｐｎ，ｐ）∀ｎ＞１
ｐ−１＝１−ｐ
ｐｎ−１＝Ｆ（ｐｎ，ｐ−１）∀ｎ＜−１
のように再帰的に定義され、且つ前記ステップｃ）において、２個の確率クラスｐｉ，ｐｊ∈Ｇｐの間の融合Ｆ（ｐｉ、ｐｊ）が次式
Ｆ（ｐｉ，ｐｊ）＝ｐｉ＋ｊ
を適用することにより計算される、請求項３に記載の方法。
【請求項8】
前記ステップｃ）が、同一位置に配置されていないセンサに関連付けられた占有格子のセルの占有確率に基づいて前記統合された占有格子を構築すべくフレームの変更を実行するステップを含んでいる、請求項１〜７のいずれか１項に記載の方法。
【請求項9】
次式
を適用することにより実装される少なくとも１個の前記センサの逆モデルを対応する占有格子上に構築する事前ステップをも含み、
ここに、
・Ｐ（ｏｉ│ｚ）が占有格子の添え字「ｉ」のセルの占有確率を表し、前記セルが前記センサからの距離の昇順に並んでおり、
・Ｎが前記占有格子のセルの個数であり、
・ｚが前記センサから得られる距離測定値を表し、
・ｐ（ｚ│ｘｋ）が、前記センサにより近い他の全てのセル以外の占有格子のセルｋの位置ｘｋが物体により占有されている場合に前記センサが与える測定値「ｚ」の確率密度を表す前記センサの直接モデルを表している、請求項１〜８のいずれか１項に記載の方法。
【数5】
ここに、
・Ｐ（ｏｉ│ｚ）が占有格子の添え字「ｉ」のセルの占有確率を表し、前記セルが前記センサからの距離の昇順に並んでおり、
・Ｎが前記占有格子のセルの個数であり、
・ｚが前記センサから得られる距離測定値を表し、
・ｐ（ｚ│ｘｋ）が、前記センサにより近い他の全てのセル以外の占有格子のセルｋの位置ｘｋが物体により占有されている場合に前記センサが与える測定値「ｚ」の確率密度を表す前記センサの直接モデルを表している、請求項１〜８のいずれか１項に記載の方法。
【請求項10】
少なくとも１個の前記占有格子の空間解像度を決定する事前ステップをも含んでいて、前記空間解像度が、ｍａｘ［Ｐ（ｏｉ│ｚｋ）］≧ｐｍｉｎであるように選択されていて、Ｐ（ｏｉ│ｚｋ）が、前記センサの逆モデルであり、ｐｍｉｎが前記集合に属し、且つ０．５よりも大きい最小確率クラスである、請求項１〜９のいずれか１項に記載の方法。
【請求項11】
少なくとも１個の前記占有格子の空間解像度を決定する事前ステップをも含んでいて、前記空間解像度が、ｍａｘ［Ｐ（ｏｉ│ｚｋ）］＝ｐｍｉｎであるように選択されていて、Ｐ（ｏｉ│ｚｋ）が、前記センサの逆モデルであり、ｐｍｉｎが前記集合に属し、且つ０．５よりも大きい最小確率クラスである、請求項１〜９のいずれか１項に記載の方法。
【請求項12】
物体を認識するシステムであって、
−１個以上のセンサから得られる、前記物体の距離測定値を表す複数の信号（ｚ１，．．．，ｚＮＣ）を受信する少なくとも１個の入力ポートと、
−前記信号を入力として受信し、それらを用いて請求項１〜８に記載の方法を適用することにより、統合された占有格子を構築すべく構成されたデータ処理モジュール（ＭＴＤ１，ＭＴＤ２）と、
−前記統合された占有格子を表す信号（ｇｆｕｓ）の少なくとも１個の出力ポートを含むシステム。
−１個以上のセンサから得られる、前記物体の距離測定値を表す複数の信号（ｚ１，．．．，ｚＮＣ）を受信する少なくとも１個の入力ポートと、
−前記信号を入力として受信し、それらを用いて請求項１〜８に記載の方法を適用することにより、統合された占有格子を構築すべく構成されたデータ処理モジュール（ＭＴＤ１，ＭＴＤ２）と、
−前記統合された占有格子を表す信号（ｇｆｕｓ）の少なくとも１個の出力ポートを含むシステム。
【請求項13】
前記物体の複数の距離測定値を表す信号を生成すべく適合されていて、前記入力ポートまたはポート群に接続された１個以上の距離センサ（Ｃ１，．．．，ＣＮＣ）をも含んでいる、請求項１２に記載のシステム。
【請求項14】
前記データ処理モジュールが、占有格子の各セルに関連付けられた確率クラスの添え字を表す各距離測定値に整数ベクトルを関連付けるセンサの逆モデルを対応関係テーブルの形式で保存しているメモリを含んでいて占有確率（ＣＯ１．．．ＣＯｎＮＣ）を計算する少なくとも１個のハードウェアブロックを含んでいる、請求項１２または１３に記載のシステム。
【請求項15】
前記データ処理モジュールが、各々の占有格子のセルに関連付けられた確率クラスの添え字を表す複数の整数を入力として受信し、前記統合された占有格子のセルに関連付けられた確率クラスの添え字を計算すべく構成された、いわゆる統合（Ｆ）用の整数計算ハードウェアブロックを含んでいる、請求項１２〜１４のいずれか１項に記載のシステム。
【請求項16】
請求項１４に依存する場合に、前記データ処理モジュールが、
−占有確率を計算する対応ハードウェアブロックから得られ、且つ占有格子の各セルに関連付けられた確率クラスの添え字を表す、少なくとも１個の第１の整数ベクトルを入力として受信し、
−前記第１の整数ベクトルを、前記統合された占有格子と空間的に一致する別の占有格子の各セルに関連付けられた確率クラスの添え字を表す第２の整数ベクトルに変換して、
−前記第２の整数ベクトルを前記統合ハードウェアブロックへの入力として提供すべく構成された、フレーム（Ｒ１，．．．，ＲＮＣ）のいわゆる変更を行う少なくとも１個の計算ブロックをも含んでいる、請求項１５に記載のシステム。
−占有確率を計算する対応ハードウェアブロックから得られ、且つ占有格子の各セルに関連付けられた確率クラスの添え字を表す、少なくとも１個の第１の整数ベクトルを入力として受信し、
−前記第１の整数ベクトルを、前記統合された占有格子と空間的に一致する別の占有格子の各セルに関連付けられた確率クラスの添え字を表す第２の整数ベクトルに変換して、
−前記第２の整数ベクトルを前記統合ハードウェアブロックへの入力として提供すべく構成された、フレーム（Ｒ１，．．．，ＲＮＣ）のいわゆる変更を行う少なくとも１個の計算ブロックをも含んでいる、請求項１５に記載のシステム。
【請求項17】
前記データ処理モジュールが、同一のセンサまたは同一位置に配置されたセンサにより連続的に取得された距離測定値を表す信号（ｚ１、．．．、ｚＮＣ）を入力として受信し、連続的な取得時点に対応する複数の前記信号に基づいて統合された占有格子を構築すべく構成されている、請求項１２または１３に記載のシステム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、計算能力およびエネルギー消費の観点から効率的な仕方でマルチセンサ融合を実行して、物体の位置を認識および推定する方法およびシステムに関係する。
【背景技術】
【0002】
「物体」とは、個別性を示し、且つ適当なセンサにより検出および識別可能な任意の物理的対象または物質を指す。従って、自然または人工の無生物、植物、動物、人間だけでなく、雲のように空気中に漂う液体または固体粒子、または実際に液体または気体の物質もまた物体であると考えられる。
【0003】
本発明は特に、ロボット、ドローン、自律走行車両等のナビゲーション分野、より一般的には認識の分野に適用できる。
【0004】
ロボットに組み込み可能な計算手段の発展に伴い、ロボット工学の応用が、工業生産からホームオートメーションまで、さらには宇宙および深海探査から一般消費市場向けの遊具ドローンまでに近年拡大している。ロボット工学用途で実行される作業は次第に複雑になっており、ロボットは未知の環境で動き回れることが求められるケースが増えている。このため、認識の手段および技術の開発、すなわち周囲の空間を発見および解釈を可能にすることが益々重要になっている。ロボット工学において認識を用いる重要な用途としてナビゲーションがあり、これはロボットの目標とする行き先を指定して、未知且つ潜在的に移動する障害を回避すべく誘導しながら、到着させることである。ロボットは従って、自身で自身の軌道を設定する責任がある。典型的な例が、注力されている研究主題となっている自律走行車である。
【0005】
最大限に死角をなくしながら環境全体を把握し、且つセンサが故障するリスクを低減させるには一般に複数のセンサを統合する必要がある。恐らくは異なる種類の複数のセンサが同一空間を占める場合、その各々から抽出された情報を組み合わせることができなければならない。従ってマルチセンサ融合に言及する。
【0006】
認識技術には主として二つの系統が存在する。表面形状手順は周囲の空間の対象の表面形状を識別することを目的とし、占有格子に基づく手順は特定の位置が障害物（より一般的には物体）により占有されているか否かを判定することを目的としている。本発明は占有格子に基づく技術に関する。
【0007】
占有格子に基づくマルチセンサ認識および融合手順の理論的基礎がＡ．Ｅｌｆｅｓによる論文「Ｏｃｃｕｐａｎｃｙ ｇｒｉｄｓ：ａ ｓｔｏｃｈａｓｔｉｃ ｓｐａｔｉａｌ ｒｅｐｒｅｓｅｎｔａｔｉｏｎ ｆｏｒ ａｃｔｉｖｅ ｒｏｂｏｔ ｐｅｒｃｅｐｔｉｏｎ」（Ｓｉｘｔｈ Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｕｎｃｅｒｔａｉｎｔｙ ｉｎ ＡＩ，１９９０）に記述されている。本文献は手順の実用的な実施には触れておらず、その直接的適用には複雑な浮動小数点計算を必要とする。
【0008】
Ｋ．Ｋｏｎｏｌｉｇｅによる論文「Ｉｍｐｒｏｖｅｄ ｏｃｃｕｐａｎｃｙ ｇｒｉｄｓ ｆｏｒ ｍａｐ ｂｕｉｌｄｉｎｇ」（Ａｕｔｏｎｏｍｏｕｓ Ｒｏｂｏｔｓ，４，３５１−３６７，１９９７）、Ｊ．Ａｄａｒｖｅ他による論文「Ｃｏｍｐｕｔｉｎｇ ｏｃｃｕｐａｎｃｙ ｇｒｉｄｓ ｆｒｏｍ ｍｕｌｔｉｐｌｅ ｓｅｎｓｏｒｓ ｕｓｉｎｇ ｌｉｎｅａｒ ｏｐｉｎｉｏｎ ｐｏｏｌｓ」（Ｐｒｏｃｅｅｄｉｎｇｓ，ＩＥＥＥ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｒｏｂｏｔｉｃｓ ａｎｄ Ａｕｔｏｍａｔｉｏｎ，２０１２）、およびＴ．Ｒａｋｏｔｏｖａｏ他による論文「Ｒｅａｌ−ｔｉｍｅ ｐｏｗｅｒ−ｅｆｆｉｃｉｅｎｔ ｉｎｔｅｇｒａｔｉｏｎ ｏｆ ｍｕｌｔｉ−ｓｅｎｓｏｒ ｏｃｃｕｐａｎｃｙ ｇｒｉｄ ｏｎ ｍａｎｙ ｃｏｒｅ」（２０１５ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｗｏｒｋｓｈｏｐ ｏｎ Ａｄｖａｎｃｅｄ Ｒｏｂｏｔｉｃｓ ａｎｄ ｉｔｓ Ｓｏｃｉａｌ Ｉｍｐａｃｔ，Ｊｕｎｅ ３０，２０１５）に占有格子に基づく技術の発展が記述されている。この場合もまた、これらの技術の実装に多くの浮動小数点計算が必要になる。
【0009】
文献米国特許第２０１４／０３５７７５号明細書、仏国特許第２００６／０５０８６０号明細書、および独国特許第１０２００９００７３９５号明細書に、地上車両の自律走行に応用された占有格子に基づくマルチセンサ認識および融合の手順とシステムが記述されている。これらの手順は全て、実装のために浮動小数点計算を必要とする。
【先行技術文献】
【特許文献】
【0010】
【特許文献1】米国特許出願公開第２０１４／０３５７７５号明細書
【特許文献2】仏国特許第２００６／０５０８６０号明細書
【特許文献3】独国特許発明第１０２００９００７３９５号明細書
【非特許文献】
【0011】
【非特許文献1】Ａ．Ｅｌｆｅｓ、「Ｏｃｃｕｐａｎｃｙ ｇｒｉｄｓ：ａ ｓｔｏｃｈａｓｔｉｃ ｓｐａｔｉａｌ ｒｅｐｒｅｓｅｎｔａｔｉｏｎ ｆｏｒ ａｃｔｉｖｅ ｒｏｂｏｔ ｐｅｒｃｅｐｔｉｏｎ」（Ｓｉｘｔｈ Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｕｎｃｅｒｔａｉｎｔｙ ｉｎ ＡＩ，１９９０）
【非特許文献2】Ｋ．Ｋｏｎｏｌｉｇｅ、「Ｉｍｐｒｏｖｅｄ ｏｃｃｕｐａｎｃｙ ｇｒｉｄｓ ｆｏｒ ｍａｐ ｂｕｉｌｄｉｎｇ」（Ａｕｔｏｎｏｍｏｕｓ Ｒｏｂｏｔｓ，４，３５１−３６７，１９９７）
【非特許文献3】Ｊ．Ａｄａｒｖｅ他、「Ｃｏｍｐｕｔｉｎｇ ｏｃｃｕｐａｎｃｙ ｇｒｉｄｓ ｆｒｏｍ ｍｕｌｔｉｐｌｅ ｓｅｎｓｏｒｓ ｕｓｉｎｇ ｌｉｎｅａｒ ｏｐｉｎｉｏｎ ｐｏｏｌｓ」（Ｐｒｏｃｅｅｄｉｎｇｓ，ＩＥＥＥ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｃｏｎｆｅｒｅｎｃｅ ｏｎ Ｒｏｂｏｔｉｃｓ ａｎｄ Ａｕｔｏｍａｔｉｏｎ，２０１２）
【非特許文献4】Ｔ．Ｒａｋｏｔｏｖａｏ他、「Ｒｅａｌ−ｔｉｍｅ ｐｏｗｅｒ−ｅｆｆｉｃｉｅｎｔ ｉｎｔｅｇｒａｔｉｏｎ ｏｆ ｍｕｌｔｉ−ｓｅｎｓｏｒ ｏｃｃｕｐａｎｃｙ ｇｒｉｄ ｏｎ ｍａｎｙ ｃｏｒｅ」（２０１５ Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｗｏｒｋｓｈｏｐ ｏｎ Ａｄｖａｎｃｅｄ Ｒｏｂｏｔｉｃｓ ａｎｄ ｉｔｓ Ｓｏｃｉａｌ Ｉｍｐａｃｔ，Ｊｕｎｅ ３０，２０１５）
【発明の概要】
【発明が解決しようとする課題】
【0012】
しかし、浮動小数点計算の実行は計算能力の観点から相当なリソースを必要とし、これは組み込みシステムに固有な制約と両立し難い。正確には、ＩＥＥＥ標準７５４により定義された浮動小数点形式では、３個の要素、すなわち符号（１ビット）、仮数（２３または５２ビット）、および指数（８または１１ビット）で数を表す。浮動小数点で表された数を用いる計算の実行は、整数計算の実行よりもはるかに複雑である（すなわち、はるかに多くの基本演算を必要とする）。これは従って、より高速なプロセッサおよび／または専用ハードウエアアクセラレータの使用を必要とし、コスト、大きさ、および電力消費の観点から好ましくない影響を伴う。
【0013】
１０ｃｍ×１０ｃｍの１０００００個のセルに離散化された２０ｍ×５０ｍの環境、および２５Ｈｚで動作する２個のセンサを一例として考える。浮動小数点計算によりＡ．Ｅｌｆｅｓの理論を実装することにより、従来技術に従いデータ融合を実行する必要がある場合、必要な計算能力は５〜５０ＧＦｌｏｐｓ（毎秒数十億回の浮動小数点演算）のオーダーである。自動車分野での用途により現実的な仮定である５０Ｈｚで動作する８個のセンサを有するセル１５００００個の格子を考えた場合、６０〜６００ＧＦｌｏｐｓの計算能力が必要となる。
【0014】
本発明は、より少ない計算リソースで済み、従って組み込み型のソリューションにより適合されたマルチセンサ融合により物体を認識する手順を提供することを目的とする。本発明は特に、−一方で、必ずしも浮動小数点演算に対応していない簡単な組み込み計算装置の使用を可能とし、
−他方で、浮動小数点計算対応している場合であっても、そのような計算の実行を回避または大幅に制限することにより、計算装置のエネルギー消費を減らすことを目的としている。そのようなエネルギー消費の低減自体が有利であり、更に、発熱を抑えることを可能になるため、放熱の観点からの制約を緩和しつつ計算装置の堅牢性および寿命が向上する。
【0015】
一例として、発明者らは、本発明により、自動車分野で用いられるＩＳＯ２６２６２ＡＳＩＬ−Ｄ認定プラットフォームで利用できるものと同等の計算性能をもたらすデジタルシステム上での融合が実行可能になる点に気付くことができた。これは、従来技術で知られている融合技術を用いても不可能である。
【0016】
以下では障害物の認識への応用を特に考慮するが、本発明はこの典型的なケースに一切限定されない。他の可能な用途は、雲や降雨、または実際に汚染物質の濃度、あるいは公共空間における人々の分布および動きの検出に関する。これら全てのケースにおいて、センサを利用して物体を検出し、その距離を測定して、測定値を用いて占有格子を計算することが可能である。
【課題を解決するための手段】
【0017】
上述の目的を実現可能にする本発明の主題は、コンピュータまたは専用デジタル電子回路により実装された物体認識方法であって以下のステップ、すなわちａ）１個以上のセンサから得られる、前記物体の複数の距離測定値を取得するステップと、
ｂ）前記距離測定値の各々に対し、前記センサの環境の離散化された空間表現を与える占有格子上の対応センサの逆モデルを適用することにより、物体による前記占有格子のセルの集合の占有確率を決定するステップと、
ｃ）各セルがステップｂで推定された占有確率の融合により計算された占有確率を示す統合された占有格子を構築するステップを含み、
前記逆センサモデルの各々が、対応する占有格子の各セルおよび各距離測定値に同一の有限集合内で選択された確率クラスを関連付ける離散的なモデルであり、前記確率クラスの各々が整数添え字により識別されること、および前記ステップｃ）において、統合された占有格子の各セルの占有確率が前記ステップｂ）で決定された確率クラスの添え字に対して実行された整数計算により決定されることを特徴とする方法である。
【0018】
本発明の別の主題は、物体を認識するシステムであって、−１個以上のセンサから得られる、前記物体の距離測定値を表す複数の信号を受信する少なくとも１個の入力ポートと、
−前記信号を入力として受信し、それらを用いて上で定義した方法を適用することにより、統合された占有格子を構築すべく構成されたデータ処理モジュールと、
−前記統合された占有格子を表す信号の少なくとも１個の出力ポートを含んでいる。
【0019】
本発明の他の特徴、詳細、および利点は、例示的に示す添付図面を参照しながら以下の記述を精査することにより明らかになろう。
【図面の簡単な説明】
【0020】
【図1】距離センサの「直接」モデルを示す。
【図2】占有格子を示す。
【図3】距離センサの「逆」モデルを示す。
【図4】占有格子上の逆モデルの空間離散化を占めす。
【図5A】占有格子上のセンサの逆モデルを計算する手順を示す。
【図5B】占有格子上のセンサの逆モデルを計算する手順を示す。
【図5C】占有格子上のセンサの逆モデルを計算する手順を示す。
【図5D】占有格子上のセンサの逆モデルを計算する手順を示す。
【図6A】占有格子の最適空間解像度を選択する手順を示す。
【図6B】占有格子の最適空間解像度を選択する手順を示す。
【図6C】占有格子の最適空間解像度を選択する手順を示す。
【図6D】占有格子の最適空間解像度を選択する手順を示す。
【図7A】様々な確率クラスのシステムを示す。
【図7B】様々な確率クラスのシステムを示す。
【図8】占有格子上のセンサの逆モデルを量子化する２個の手順を示す。
【図9A】本発明の第１の実施形態による障害物認識システムを示す。
【図9B】本発明の第１の実施形態による障害物認識システムを示す。
【図10A】本発明の第２の実施形態による障害物認識システムを示す。
【図10B】本発明の第２の実施形態による障害物認識システムを示す。
【図11】本発明の第３の実施形態による障害物認識システムを示す。
【発明を実施するための形態】
【0021】
以下の詳細な記述において、障害物認識のケースに言及している。しかし、全ての記述内容はより一般的に任意の種類の物体の認識に適用できる。
【0022】
通常、ナビゲーションに用いられるセンサは、周囲の障害物との距離に関して助言するものであり、以下では距離センサについて言及する。センサの精度、生じ得る誤差、または解像度を考察すべく確率モデルを導入する。その趣旨は、センサが出力した測定値が必ずしも障害物とセンサとの正確な距離を示す訳ではなく、従ってセンサの反応を把握して障害物が所与の距離にある確率の観点から推論することが適切であるというものである。
【0023】
ｄを障害物とセンサとの実際の距離、ｚをセンサの出力とすれば、障害物の実際の位置と、センサが検知したその推定値との関係をモデル化する（「直接モデル」）条件付き確率密度関数ｐ（ｚ｜ｄ）に関心が持たれる。図１にセンサの例示的な直接モデルを示す。長さ５０ｍの線形空間は考え、センサからｄ＝２５ｍに障害物が位置していると仮定する。ガウス関数によりモデル化可能な誤差を有するセンサの場合、最尤反応ｚは２５ｍ付近であるが、曲線により定義される確率密度により他の値も可能である。理想的なセンサの場合、ｐ（ｚ｜ｄ）＝δ（ｚ−ｄ）が得られ（δはディラックのデルタ）、且つ測定値は常に真の距離に等しい。センサの直接モデルは実験的に決定することができ、典型的には、構築子が与えるデータ（ガウス関数の場合、モデルの特徴付けに標準偏差の値で充分である）に基づいて構築することができる。
【0024】
以下において、Ωで１、２、または３次元の空間ベンチマークを表すものとする。占有格子ＧＯは、Ωの連続的且つ有界な部分集合を、セルと呼ばれ、添え字ｉ∈［０，Ｎ−１］で示すＮ個の部分に分割したものである。添え字ｉのセルをｃｉで示す。一般性を失うことなく、以下において、１個の距離センサＣ（または同一位置に配置された複数のセンサ）により観察される１次元の占有格子を考える。遠いセンサほど大きい添え字ｉが付される（従って、ｃ０はセンサに最も近いセル、ｃＮ−１は最も遠いセンサである）。本構成を図２に示す。
【0025】
障害物ＡはΩの有界な連続部分集合である。Ａ∩ｃｉ≠φの場合セルｃｉは障害物Ａにより占有されていると言い、Ａ∩ｃｉ＝φの場合Ａにより占有されていないと言う。換言すれば、障害物がセルの一部でも覆うならば、当該セルは占有されているとみなす。他の方式も可能であるが、いずれの場合もセルは空いているか占有されているかのいずれかでなければならない。
【0026】
格子の各セルに対して、当該セルが障害物を含んでいるか否かを把握する、２通りの結果｛占有、空｝の一方を取り得る２値ランダム実験「状態」を考える。セルｃｉの状態をｅｉで示し、ｏｉは結果ｅｉ＝占有を示し、ｖｉは結果ｅｉ＝空を示す。１個の格子において全てのセルが独立していると考えられるため、次式が成り立つ。∀ｉ，ｊ∈［０，Ｎ−１］，Ｐ（ｏｉ∧ｏｊ）＝Ｐ（ｏｉ）・Ｐ（ｏｊ） （１）
ここに、∧は論理演算子「ａｎｄ」であり、Ｐ（．）は事象の確率を示す（小文字「ｐ」で表す確率密度と混同しないように）。
【0027】
また、障害物の位置は、不確実な距離センサを利用することでしか分からないものと考えられ、より一般的に【数1】
【数2】
【0028】
センサから得られる測定値ｚにより、セルｃｉの占有確率Ｐ（ｏｉ｜ｚ）を決定することが可能になる。所与の測定値ｚに対して、確率の集合Ｐ（ｏｉ｜ｚ）∀ｉ∈ ［０，Ｎ−１］は格子上のセンサの逆モデルを構成する。センサの直接モデルはセンサの反応を現実世界の関数と見なして示唆するのに対し、逆モデルは、採用した現実世界のモデルである占有格子に対する測定の影響を表すため、逆モデルという用語が正当化される。
【0029】
図３に、ｚ＝２５ｍである場合の距離センサの逆モデルの典型的な例を示す。センサから２４．２５ｍ未満の距離にあるセルに対して占有確率が準ゼロであって、距離２５ｍ（センサが与える測定値に対応する）でピークに達することが確認できよう。２５ｍを越えた場合、占有確率は、障害物よりも遠方にあって、当該障害物により遮蔽されているためセンサに届かないセルの占有状態が全く分からないこと示す値０．５で安定するまで低下し続ける。
【0030】
本明細書での主な適用例によれば、図３は平滑化された曲線を用いる逆モデルを表す。より正確な表現は、格子のセルの限度に対応する位置だけを示すものであろう。実際に、「部分的に」占有されたセルを「完全に」占有された別のセルから区別するのは不可能であり、全てのケースにおいて障害物までの距離を、対応するセルまでの距離として推定する。これは格子によりもたらされる空間誤差である。
【0032】
「占有」および「障害物距離」の概念が完全に同等である訳ではない点を強調すべきである。実際、センサ距離ｚに障害物にあると言った場合、単に特定のセルが占有されていることだけなく、当該センサにより近い格子の他のセルは空いていることも示す（さもなければ第１の障害物はｚ未満の距離にあった筈である）。上述の図２において、障害物Ａは添え字ｉのセル（黒）に位置しており、添え字がｉよりも小さいセルは空いていることを示すため白く描かれ、添え字がｉよりも大きいセルは占有状態が未知であることを示すためにグレーに描かれている。
【0033】
自身の（直接）モデル【数3】
【数4】
【数5】
【0034】
式（２）は、格子（ｘｉ）のセルの境界上の点で評価されるセンサのモデルが、対応する格子構成、すなわち、セルｉよりも近いセルが空であり、セルｉは占有されていて、セルｉよりも遠方のセルの占有状態が決定されていない格子に対するセンサの反応の確率密度に等しいことを示す。Ａ．Ｅｌｆｅｓが上述の論文でセンサから逆モデルを構築する方法を提案したのはこの情報を利用したことによる。本方法について以下で説明する。
【0035】
ベイズの定理により、センサの逆モデルＰ（ｏｉ｜ｚ）を次式【数6】
【0036】
従って次式【数7】
【0037】
項ｐ（ｚ｜ｏｉ）に対して可能な格子構成は形式【数8】
【数9】
【0038】
従って、次式【数10】
【0039】
項【数11】
【数12】
【数13】
【数14】
【数15】
【数16】
【0040】
上述の手順の主な制約は、式（６）の和の計算を現実的に不可能にする、和の項数の指数関数的爆発から生じる。実際、従来技術では、逆モデルの解析的式を、直接ｐ（ｚ｜ｏｉ）の関数として表すことなく一般に用いている。これは、実験により直接得られる唯一の関係である直接モデルとの関係、従ってモデリングにより生じる誤差を測定する能力が失われることを意味する。
【0041】
本発明の第１の目的は、式（６）、（４）の使用に関して一切近似値を導入することなく、Ｎの指数関数ではなく線形な逆モデルを計算する簡素化された手順である。以下に、一般性を失うことなく本開示を簡素化する目的で、１次元占有格子のケースを考えることにし、障害物の位置をスカラー変数ｘで表し、この場合、式（６）を次式【数17】
【0042】
Ｎ個のセルを有する１次元格子を考えているため、ｘｋはセルｃｊ、ｊ∈［０〜Ｎ−１］のうち１個の位置の値しか取り得ない。従って、和（７）の項はＮ個の異なる値しか取り得ない。ｐ（ｚ｜ｘｋ）の同じ値を与える全ての格子を因数分解することにより、式（７）を２Ｎ−１ではなくＮ個の項の和に簡約化することができる。計算の複雑度は指数関数から線形関数となる。
【0043】
まずｅｉ＝ｏｉの場合を考えることから始める。ｃｉを占有された状態に固定し、第１の障害物の距離をセルｃｋの距離あるように、従って当該距離が実際にｘｋであるように固定する。従って、セルｃｉが占有されていることを知った上で第１の障害物が位置ｘｋに見られるように格子の個数を求める。３個の典型的なケースが考えられる。
【0044】
１．ｋ＜ｉ。セルｋは占有され（１個のセル）、ｊ＜ｋのセルｃは空であり（ｋ個のセル）、セルｃｉは占有されていて（１個のセル）、他は占有されているかまたは空である。これを図５Ａに示し、空きセルを白、占有されているセルを黒、および不確定状態のセルをグレーで表している。従って、固定状態のｋ＋２個のセルおよび不確定状態のＮ−ｋ−２個のセルがあり、そのような格子が２Ｎ−ｋ−２個ある。２．ｋ＝ｉ（図５Ｂ参照）。セルｋ＝ｉは占有され（１個のセル）、ｊ＜ｉのセルｃｊは空であり（ｉ個のセル）、他は占有されているかまたは空である。この結果、状態が固定された１＋ｉ個のセルおよび状態が不確定のＮ−ｉ−１個のセルが生じる。結局そのような格子が２Ｎ−ｉ−１個存在する。
３．ｋ＞ｉ。この構成は不可能である。その理由は、ｋが厳密にｉよりも大きい場合、最も近い障害物はｘｋではなくｘｉに見られるであろう。
【0045】
従ってｅｉ＝ｏｉのケースで式（７）は次式【数18】
【0046】
ｅｉ＝ｖｉのケースに同じ推論を繰り返すことができる。この場合もまた３通りの可能性を区別することができる。１．ｋ＜ｉ、すなわち図５Ｃに示す状況。この条件を満たす２Ｎ−ｋ−２個のセルがある。
２．ｋ＝ｉ、すなわち不可能な状況。
３．ｋ＞ｉ、すなわち図５Ｄに示す状況。この条件を満たす２Ｎ−ｋ−２個のセルがある。
【0047】
従ってｅｉ＝ｖｉのケースの式（７）は次式【数19】
【0048】
式（８）および（９）を（４）に代入することにより、格子のサイズＮに関して、線形の複雑度を有する直接モデルに基づいて占有格子上の１次元逆センサモデルを構築することができる。【数20】
【0049】
このような簡素化はまた、特に極または球面座標を用いてより大きい次元でも可能である。
【0050】
逆モデルの構造は格子の定義に大いに依存する。従って、逆モデルに対する空間解像度の変化の影響を調べることが興味深い。図６Ａ〜６Ｄに、４個の異なる空間解像度格子１ｍ（６Ａ）、５０ｃｍ（６Ｂ）、２５ｃｍ（６Ｃ）、１２．５ｃｍ（６Ｄ）に対する同一のセンサの逆モデルを示す。解像度が増大する（格子の間隔が狭まる）につれて、逆モデルの最大値が減少し、０．５に収束する傾向がある点に注意されたい。実際、センサの精度よりも高い精度で障害物の位置を知ることができると期待すべきではない。逆に、センサの精度よりも極めて低い精度で占有を把握することで満足できるならば、障害物の有無を高い確度（逆モデルの最大が０．９９９９９４に等しい図６Ａのケース）で決定することができる。これらの考察により、格子の空間解像度を最適化することが可能になる。実際、逆モデルの最大値が「重要」と考えられる閾値よりも大きいままである（０．５よりも厳密に大きく、且つ厳密に１未満）格子の最大解像度を決定可能にする探索を行うことが可能である。
【0051】
センサの精度と格子の解像度との間に存在する関係が明らかなのは逆モデルが直接モデルに基づいて計算される（７〜１０）場合に限る点に注目すると興味深い。この関係は、従来技術のように、逆モデルの近似的な解析的式で満足したならば失われる。これは本発明が提案する方法の追加的な利点となる。
【0052】
同一の占有格子上の２個のセンサの逆モデルに基づいて、２個のセンサのデータの融合を次式【数21】
【0053】
この浮動小数点計算は、格子の各セルに対して、少なくともセンサの取得周波数と同程度の周波数で実行されなければならず、従って相当の計算能力を必要とする。
【0054】
本発明の第２の目的は、占有格子における複数のセンサから得られるデータを、浮動小数点計算を行わずにベイズ的に融合する方法であり、従って占有格子の構築に必要な計算の複雑度を大幅に低下させ、従ってより広範な応用分野、特に組み込み計算容量が極めて制限される分野に適用可能にする。
【0055】
本発明の上述の態様は、区間［０，１］における確率を、整数添え字により識別される確率クラスにより離散的に表すものである。
【0056】
以下において、［０，１］の可算部分集合の要素ｐｎは従って、相対整数添え字「ｎ」により特徴付けることができ、「確率クラスの系」【数22】
【数23】
【0057】
特に興味深いケースは、系の２個の確率クラスの融合の結果もまた当該系に帰属するようなクラスの系である。形式的には∀ｐｉ，ｐｊ∈Ｓ、Ｆ（ｐ１、ｐ２）∈Ｓである。次いで、融合が一切誤差または近似をもたらさないため、クラスの「無誤差」系に言及する。従って、対応クラスの添え字により確率値にラベル付けすることができ、融合の結果もまた、添え字により識別される。従って融合問題は、２個の整数添え字に別の整数添え字を関連付ける適当な機能Ｆｄを決定することである。形式的に、【数24】
【0058】
Ｆｄ（ｋ，ｌ）の計算は、添え字ｋとｌ、および整数添え字演算の知識だけを必要とする。情報ｐｋ、ｐｌの融合の計算に浮動小数点計算は不要である。更に、クラスの系を考えた場合、Ｆｄ（ｋ，ｌ）を利用して得られた添え字は、式（１１）を適用することにより浮動小数点数を用いて得られた値と厳密に同一の確率値を示す。当該手順は従って、浮動小数点計算の観点から無誤差の確率クラスの融合を可能にする。
【0059】
系Ｓが個々に無誤差であるクラスの複数の部分系の和集合として構築されるケースを考えることにより上述の方式を一般化することが可能である。この場合、全体としての系Ｓは無誤差とは限らない。従って、融合関数Ｆｄの定義に近似ステップを導入する必要がある。【数25】
【0060】
従って、近似値モデルに応じて、Ｆｄ（ｋ，ｌ）＝ｉまたはＦｄ（ｋ，ｌ）＝ｉ＋１を切り捨て、切り上げ、または四捨五入することにより選択することができる。しかし、誤差は依然としてクラスの系が対応する幅により抑えられる。
【0061】
無誤差系の自明な例がＳ＝｛１／２，１｝である。１／２，１、および０と異なる確率を含むクラスのあらゆる系が必然的に無限個の要素を含んでいる。実際、明らかな実装上の理由のために、有限個の確率クラスの系だけを考える。しかし、センサの個数が有限であり、それらの出力（量子化およびデジタル化された）が有限個の値しか取れないことを前提とすれば、有限個の確率クラスの系に基づく場合であっても「無誤差」融合が実行可能であろうことが証明されている。
【0062】
第１の例示的な応用上の利点は、以下のクラスの２個の無誤差部分系の和集合【数26】
【数27】
【数28】
−０．５の確率は不確実な占有を示し、従ってこの値の付近で精度が極めて高いことは不要である一方、極値０および１の付近では精度は有用である。
−｜ｎ｜のある値を越えたならば、確率クラスｐｎの様々な値は互いに極めて近くなり、従ってクラスの系の切り捨てにより生じる誤差は無視できる。
【0063】
０または１に極めて近い確率だけが精密にサンプリングされるのを防ぐべく、ｋの小さい値（例えば５を越えない）を選択することが好適である。実際、センサの精度が高いほどｋの値を大きくすることができる。
【0064】
図４の逆モデルは、空間的に離散化されているが、各セルの占有確率は、区間［０，１］内で任意の値を取ることができる。現時点では、クラスＳｋの系（より一般的にはクラスＳの系）に属するように、確率値自体が離散化または量子化されている。従って、本発明によれば、図４の逆モデルの確率値をクラスＳの系の要素で近似する必要がある。
【0065】
第１の可能性は、量子化誤差を最小化すべく、図８の曲線ＭＩで表される逆モデルの値をＳの最も近い要素で代替することである。その結果は、確率クラスの系がＳ１（Ｓｋでｋ＝１）である場合、同じ図８の曲線ＭＱＰで表されている。この方式ではセルの占有確率が過小に推定される恐れがあることが分かり、これは障害物検出用途において受容できない恐れがある。代替方式は、理論上の逆モデルの値を、クラスＳの系（同じく系Ｓｋの場合に、図８の曲線ＭＱＥ）の最小上界で近似する。従って、占有確率が過小に推定されることは絶対になく、これは障害物の検出にとって利点であろう。人数を数える等、他の用途では、この種の近似は一方で偽陽性を引き起こす恐れがある。
【0066】
考える近似値の種類が何であれ、空間的に離散化された逆系（図８の曲線ＭＩ）は極めて少ない個数の値（図８の例では７個、センサの誤差によるが格子の最大解像度を取れば最大１８個）を取る。従って、逆モデルでの近似に必要な系Ｓの要素の個数もまた極めて少ない。従って考察する対象を、原理的に無限且つ可算クラスＳｋの系の小さいサイズの有限部分集合に限定してもよい。
【0067】
図６Ａ〜６Ｄを参照しながら上で述べたように、所与のセンサに対して、空間占有格子の解像度が高いほど、逆モデルの最大値は０．５に近づく。しかし、クラスＳの系に属する確率値に限定すれば、０．５よりも大きい最小値に対応する確率クラスｐｍｉｎが存在する。Ｓ１のケースではｐｍｉｎ＝ｐ１＝２／３であり、より一般的には、クラスＳｋの系ではｐｍｉｎ＝ｐ１＝（ｋ＋１）／（ｋ＋２）である。占有格子の最適解像度は従って、逆モデルの最大値がｐｍｉｎに等しいものである。これ以降は格子の間隔を狭くしても何らの情報も得られないまま計算負荷が増大する。
【0068】
上記にて説明したように、タイプＳｋの確率クラスの系を用いる利点は、データ融合に添え字「ｎ」すなわち整数の計算だけを必要とする点である。
【0069】
最初に【数29】
【数30】
Ｆ（ｐｉ，ｐｊ）＝ｐｉ＋ｊ−ｋ・ｉ・ｊ （１５）
が得られる。
【0070】
最初に【数31】
【数32】
Ｆ（ｐｉ，ｐｊ）＝ｐｉ＋ｊ＋ｋ・ｉ・ｊ （１６）
が得られる。
【0071】
ｉ＝０ならば、確率ｐｉは０．５に等しく、従ってＦ（ｐｉ，ｐｊ）＝ｐｊとなり、占有確率が０．５であるセンサは何らの情報も与えない。
【0072】
融合すべき確率の一方が０．５未満であるためクラス【数33】
【数34】
【数35】
【0073】
整数比の計算結果は一般に整数とはならない。しかし、結果が整数、すなわち値が実際の結果に近い確率クラスの添え字となる整数除算「÷」を実行することは可能である。上式は従って次式・ｉ＜０，ｊ＞０且つ｜ｉ｜≧｜ｊ｜ならば、Ｆ（ｐｉ，ｐｊ）＝ｐ（ｉ＋ｊ）÷（１＋ｋ・ｊ） （１９）
・ｉ＜０，ｊ＞０且つ｜ｉ｜≦｜ｊ｜ならば、Ｆ（ｐｉ，ｐｊ）＝ｐ（ｉ＋ｊ）÷（１−ｋ・ｉ） （２０）
となる。
【0074】
式（１９）および（２０）の適用によりもたらされる誤差は、２個の連続するクラス間の最大距離により抑えられる。最大距離の２個のクラスはｐ０およびｐ１であり、従って最大の誤差は次式Ｅ（ｋ）＝ｐ１−ｐ０＝ｋ／（２ｋ＋４） （２１）
により与えられる。
【0075】
ここで考える状況では２個のセンサが矛盾する情報を与える点に注意されたい。特定のケースでは、式（１９）、（２０）ではなく、従来技術で知られた矛盾を管理する規則を適用する方が好ましいであろう。例えば、特に潜在的に危険な衝突を回避する必要がある場合、占有確率の最も高い推定値を取ることが好ましいかもしれない。この場合、単に次式Ｆ（ｐｉ，ｐｊ）＝ｍａｘ（ｐｉ，ｐｊ） （２２）
が得られる。
【0076】
応用上の利点を示すクラスの別の例示的な系は再帰的に定義することができる。
【0077】
ｐを厳密に０．５と１の間に存在する占有確率、すなわち０．５＜ｐ＜１とする。従って数列ｐｎが次式ｐ０＝０．５
ｐ１＝ｐ
ｐ２＝Ｆ（ｐ，ｐ）
ｐ３＝Ｆ（ｐ２，ｐ）
．．．
ｐｎ＋１＝Ｆ（ｐｎ，ｐ） （２３）
のように再帰的に定義される。
【0078】
次いで、ｐｎの定義を次式ｐ０＝０．５
ｐ−１＝１−ｐ
ｐ−２＝Ｆ（ｐ−１，ｐ−１）
ｐ−３＝Ｆ（ｐ−２，ｐ−１）
．．．
ｐｎ−１＝Ｆ（ｐｎ，ｐ−１） （２４）
のようにｎの負の整数値に拡張する。
【0079】
次いで、パラメータｐ∈］０．５，１［を用いて以下のクラスの２個の系【数36】
【数37】
【数38】
【数39】
【0080】
上式から、【数40】
【数41】
【0081】
【数42】
【0082】
式（２７）によりＧｐ内のクラスを融合する次式【数43】
【0083】
更に、パラメータｐにより、量子化により生じた誤差を微細に制御することが可能になる。実際、ｐ＝１／２＋εと置けば次式Ｅ（ｐ）＝ｐ１−ｐ０＝ε （２９）
が得られる。
【0084】
系Ｇｐが極めて興味深いのは、融合全体を可能な最も簡単な方式（整数追加）で実行可能且つ無誤差にし、且つパラメータｐの選択により、量子化により生じた誤差を抑制することが可能になる。
【0085】
図７Ｂにパラメータｐの３個の値すなわち０．５２、０．５５、および０．７について系Ｇｐを示す。
【0086】
図９Ａに、距離センサＣ、例えば機械式走査レーザーテレメータ（ＬＩＤＡＲ）を備えた地上車両Ｖへの本発明の応用を示す。当該センサは車両前面の空間に対し１次元走査を複数回の実行し、各走査が取得「シート」を画定すべく構成されている。好適には、異なる高さで複数のシートＮ１、Ｎ２．．．が生成される。各走査を行う間、センサは測定値ｚのベクトルを生成し、各測定値は、各方向における歩行者、他の車両、道路の端の木等、障害物の存在、およびその距離を表す（例えばｚが最大値を取るとき、センサの走査範囲内で障害物が検出されなかったことを示す）。一変型例として、同一位置に配置された（すなわち、距離測定の同一原点を有する）複数のセンサにより複数の取得シートを同時に生成することが可能になる。
【0087】
図９Ｂに、本出願に適した障害物認識システムを示す。本システムは、前記センサＣ（または同一位置に配置されたセンサの組）、およびセンサの各取得シートに対応する測定値のベクトルを入力として受信し、これらの取得シートのデータを融合して得られた占有格子を表す信号（典型的には整数のベクトル）を自身の出力端から配信するデータ処理モジュールＭＴＤ１受信を含んでいる。
【0088】
図９Ｂの実施形態において、データ処理モジュールＭＴＤ１は、占有確率を計算する複数のハードウェアブロックＣＯ１．．．ＣＯＮＣ、およびいわゆる統合または融合計算用のハードウェアブロックＦを含んでいる。占有確率を計算する各ブロックＣＯｋは、確率クラスの系、例えばＳ１により離散化された、センサＣの添え字ｋのシートの逆モデルを対応する文書テーブルの形式で格納しているメモリを含んでいる。ここで、融合された各種のシートの測定値であるため「シートの逆モデル」に言及する。複数のシートを取得すためにセンサを１個だけ用いる場合、当該１個のセンサは実際、各々がシートを取得する複数のセンサと同等であり、（たとえこれらの逆モデルが全て同一であっても）各々が自身の逆モデルを有している。
【0089】
各処理ブロックＣＯｋは従って各取得シートｚｋ（すなわちｚ１．．．ｚＮＣ）に対応する測定値を入力として受信し、占有格子を出力として格子の各種セルに関連付けられた確率クラスの添え字を表す整数ｇｋのベクトルの形式で配信する。従って、格子ｇｋは、シートｋ単独の測定値、すなわち測定値ｚｋのベクトルを利用して推定された占有情報を含んでいる。
【0090】
統合ハードウェアブロックＦは、各々が式（１５）、（１６）、（１９）および（２０）を実装した４個の組み合わせ論理回路Ｆ１、Ｆ２、Ｆ３およびＦ４を含んでいる。統合ハードウェアブロックＦは自身の入力端で占有格子ｇ１．．．ｇＮＣを受信し、整数すなわち当該統合格子の各種セルに関連付けられた確率クラスの添え字のベクトルで表された「統合」占有格子を自身の出力端から配信する。
【0091】
上記にて説明したように、式（１９）および（２０）は実装されない場合があり、従って、ブロックＦ３およびＦ４は存在しないか、または検出の矛盾を管理する、例えば式２１で表される規則（最大確率の選択）を実装する回路で代替されてもよい。
【0092】
各種の取得シートに関連付けられた逆モデルが同一ならば、ブロックＣＯ１．．．ＣＯＮＣも同一であって、連続的に処理を実行しながら占有確率を計算する単一のハードウェアブロックで代替することができる。
【0093】
データ処理モジュールＭＴＤ１はまた、他の任意の種類の距離センサに関連付けられていてよい。
【0094】
図１０Ａ、１０Ｂは、各種の視点から実行される測定を利用して構築される占有格子を提供すべく協働する、異なる位置に配置された複数のセンサを用いる本発明の別の実施形態に関係している。これらのセンサは、精度、範囲、視野および／または取得速度の観点から技術的に多様であってよい。本実施形態において、第１の障害物の距離は測定を行うセンサに関する情報の項目である。シナリオの模式的な例を様々な位置に配置されていて異なる範囲および視野を有する２個のセンサＣ１、Ｃ２を示す図１０Ａで表す。従って、障害物ＯはＣ１、Ｃ２により完全に異なる距離で観察される。
【0095】
本実施形態において、主な困難さは、占有格子とセンサが各々自身に関連付けられた固有のフレームを有していることである。従って、障害物の位置を評価するためにフレームの変更を実行する必要がある。
【0096】
図１０Ｂに、本発明の上述の実施形態による障害物認識システムを示す。本システムは、一般に、同一位置には配置されておらず、且つ潜在的に異なる種類の「ＮＣ」個のセンサＣ１、Ｃ２．．．ＣＮＣおよびデータ処理モジュールＭＴＤ２を含んでいる。データ処理モジュールＭＴＤ２は、占有確率を計算するハードウェアブロックＣＯ１．．．ＣＯＮＣと、統合ハードウェアブロックＦとの間に挿入されたフレーム変更ブロックＲ１．．．ＲＮＣをも含んでいる点で図９Ｂのデータ処理モジュールＭＴＤ１とは異なる。これらのブロックＲｋの各々は、各センサのフレームから、データ融合が実行されるいわゆる「統合」占有格子のフレームへの変更を実行する一般に浮動小数点計算装置を含んでいる。フレームを変更すべく実行される計算は、センサＣｋのフレーム内で既知の位置の占有（整数ｇｋのベクトルで表す）を統合格子のフレームの対応セルに再び割当てるものである。統合格子のセルの占有を表す整数のベクトルを【数44】
【0097】
一変型例として、フレームを更用する式は、モジュールＲｋに含まれるメモリに保存された変換テーブルに保持することができる。従って、当該ケースであっても、浮動小数点計算を迂回して整数に対する演算だけを実行することが可能である。一方、これらの変換テーブルは相当大量になり得るため、それらの保存にはシリコン表面積の観点から無視できないコストを要する場合がある。
【0098】
図１１に、図示しない単一のセンサＣがスカラーまたはベクトル測定値ｚｔ，ｚｔ＋１，．．．ｚｔ＋ｍ．．．を所定の用途で必要とされるよりもＮ倍速い取得速度で取得する本発明の第３の実施形態を示す。占有確率ＣＯを計算するハードウェアブロックは測定値の各々に対して占有格子ｇｔ，ｇｔ＋１，．．．ｇｔ＋ｍ．．．を生成する。次に、ハードウェア融合ブロックＦが、連続する時点で取得された当該格子のうちＮ個を単一の統合された格子ｇｆｕｓに融合する。統合された占有格子ｇｆｕｓは従って、センサ用途よりもＮ倍遅い速度で生成される。例えば、センサＣが１００Ｈｚの速度動作するのに対し、想定する用途では１０Ｈｚの速度で充分な場合、１０個の連続的に取得された値を融合することができる。
【0099】
図９Ｂ、１０Ｂ、および１１の実施形態において、複数の処理または少なくともそのいくつかをハードウェア計算ブロック、すなわち専用デジタル回路で実行するケースを考えた。しかし、本発明はまた、完全にまたは部分的に、複数の処理または少なくともそのいくつかを適切にプログラミングされた汎用プロセッサにより実行されるソフトウェアにより実装されていてよい。