(57)【特許請求の範囲】
【請求項１】
第１電子機器と、任意の時間に移動可能な複数の第２電子機器と、を有する第１の電子機器群を備える、少なくとも時間の経過に伴って形状が変化する物体を含む被写体の表面の各点までの深度情報を画素単位に有する点群データを取得する点群データ取得システムであって、
前記第２電子機器のそれぞれは、
固有の第２座標系を有し、
複数の、視覚的特徴を有する所定の標識物であって、前記複数の標識物はそれぞれ第２座標系において線形独立となる位置にそれぞれ設けられる、前記複数の標識物と、
１秒あたりに予め設定した複数の頻度となるように設定された所定の周期で前記被写体の点群データを前記第２座標系に基づいて連続して計測し、計測した時刻を示すタイムスタンプを付した前記点群データを取得する３次元計測部と、
を備え、
前記第１電子機器は、
第１座標系を有し、
被写体を前記第１座標系に基づいて計測するとともに、計測した時刻を取得する位置情報計測部を備え、
前記点群データ取得システムは、さらに
基準座標系を有し、
前記第２電子機器それぞれに対して、前記位置情報計測部により計測されたデータ及び計測時刻に基づいて、前記第２電子機器それぞれの備える前記複数の標識物の前記第１座標系における座標値である第２電子機器位置情報を前記計測時刻に対応付けて算出するマーカー位置情報算出部と、
前記計測時刻に対応付けて算出された前記第２電子機器位置情報に基づいて、前記第２電子機器それぞれの備える前記３次元計測部により前記所定の周期で連続して計測された前記点群データにおける各点のデータの前記第２座標系における座標値を、前記タイムスタンプごとに前記基準座標系における座標値に変換する点群データ座標値変換部と、
前記第２電子機器それぞれの備える前記３次元計測部により前記所定の周期で連続して計測された複数の点群データを、前記点群データ座標値変換部により変換された前記基準座標系における座標値に基づいて、前記所定の周期ごとに１つの合成点群データを作成する合成点群データ作成部と、
を備える、点群データ取得システム。
【請求項２】
前記第２電子機器の備える視覚的特徴を有する所定の標識物は、それぞれ表面が予め設定された色で着色された、所定の半径の球であって、前記球の中心がそれぞれ予め設定された三角形の頂点を形成するとともに、前記第２電子機器を一意的に識別できるように、前記複数の標識物の視覚的特徴を前記第２電子機器ごとに異なるように構成されたことを特徴とする、請求項１に記載の点群データ取得システム。
【請求項３】
前記第１の電子機器群は、２つ以上の第１電子機器を備え、
１つの第１電子機器の備える位置情報計測部、及び別の第１電子機器の備える位置情報計測部が、それぞれの備える固有の第１座標系に基づいて、前記第２電子機器の備える前記複数の標識物を計測することを特徴とする請求項１又は請求項２に記載の点群データ取得システム。
【請求項４】
前記点群データ取得システムは、さらに、複数の、視覚的特徴を有する所定の標識物であって、前記複数の標識物はそれぞれ基準座標系において線形独立となる座標位置にそれぞれ設けられる、前記複数の標識物を備え、
前記点群データ座標値変換部は、前記マーカー位置情報算出部により算出された、前記複数の、視覚的特徴を有する所定の標識物の前記第１座標系における座標値、及び前記複数の、視覚的特徴を有する所定の標識物の前記基準座標系における座標値に基づいて、予め前記第１座標系における座標値を前記基準座標系における座標値に変換する基準座標変換行列を生成し、
前記点群データ座標値変換部は、前記第２電子機器位置情報、前記点群データ、及び前記基準座標変換行列に基づいて、前記点群データにおける各点のデータの前記第２座標系における座標値を、前記基準座標系における座標値に変換する、
請求項１乃至請求項３のいずれか１項に記載の点群データ取得システム。
【請求項５】
前記点群データ取得システムは、さらに、
前記第１の電子機器群に属さない、１つの第１電子機器と１つ以上の第２電子機器を備える第２の電子機器群を有し、
前記第２の電子機器群に属する第１電子機器は、さらに、
複数の、視覚的特徴を有する所定の標識物であって、前記複数の標識物はそれぞれ第１座標系において線形独立となる位置にそれぞれ設けられる、前記複数の標識物を備え、
前記第１の電子機器群に属する前記第１電子機器の前記位置情報計測部は、第１座標系に基づいて前記第２の電子機器群に属する前記第１電子機器の備える前記複数の標識物を計測し、
前記マーカー位置情報算出部は、前記第２の電子機器群に属する前記第１電子機器の備える前記複数の標識物の前記第１座標系における座標値である、前記第２の電子機器群に属する前記第１電子機器の位置情報を算出する、請求項１乃至請求項４のいずれか１項に記載の点群データ取得システム。
【請求項６】
第１電子機器と、任意の時間に移動可能な複数の第２電子機器と、を有する第１の電子機器群を備える、少なくとも時間の経過に伴って形状が変化する物体を含む被写体の表面の各点までの深度情報を画素単位に有する点群データを取得する点群データ取得システムであって、
前記第１電子機器は、
第１座標系を有し、
前記第２電子機器それぞれに対して、電磁波を発信する、複数の電磁波発信部であって、前記複数の電磁波発信部は、それぞれ第１座標系の線形独立である任意の位置に設けられるとともに、計測した時刻を取得する前記複数の電磁波発信部と、
を備え、
前記第２電子機器は、
それぞれが固有の第２座標系を有し、
電磁波を受信又は反射する、複数の電磁波受信部であって、前記複数の電磁波受信部は、それぞれ前記第２座標系の線形独立である任意の位置に設けられる、前記複数の電磁波受信部と、
１秒あたりに予め設定した複数の頻度となるように設定された所定の周期で前記被写体の点群データを前記第２座標系に基づいて連続して計測し、計測した時刻を示すタイムスタンプを付した前記点群データを取得する３次元計測部と、
を備え、
前記点群データ取得システムは、さらに
基準座標系を有し、
前記複数の電磁波受信部のそれぞれが、計測時刻に対応付けて前記複数の電磁波発信部のそれぞれから発信された電磁波を受信するのに要した電磁波伝播時間、又は前記複数の電磁波発信部のそれぞれが、前記複数の電磁波受信部のそれぞれから反射された電磁波を受信するのに要した電磁波伝播時間に基づいて、前記複数の電磁波発信部のそれぞれと前記複数の電磁波受信部それぞれとの距離を測定する距離計測部と、
前記距離計測部により計測時刻に対応付けて測定された距離に基づいて、前記複数の電磁波受信部の前記第１座標系における座標値である第２電子機器位置情報をそれぞれ計測時刻に対応付けて算出する第２電子機器位置情報算出部と、
前記計測時刻に対応付けて算出された前記第２電子機器位置情報に基づいて、前記第２電子機器それぞれの備える前記３次元計測部により前記所定の周期で連続して計測された前記点群データにおける各点のデータの前記第２座標系における座標値を、前記タイムスタンプごとに前記基準座標系における座標値に変換する点群データ座標値変換部と、
前記第２電子機器それぞれの備える前記３次元計測部により前記所定の周期で連続して計測された複数の点群データを、前記点群データ座標値変換部により変換された前記基準座標系における座標値に基づいて、前記所定の周期ごとに１つの合成点群データを作成する合成点群データ作成部と、
を備える、点群データ取得システム。
【請求項７】
前記点群データ取得システムは、さらに、
前記第１の電子機器群に属さない、１つの第１電子機器と１つ以上の第２電子機器を備える第２の電子機器群を有し、
前記第２の電子機器群に属する第１電子機器は、さらに、
電磁波を受信又は反射する、複数の電磁波受信部を備え、
前記距離計測部は、さらに前記第１の電子機器群に属する前記第１電子機器の前記複数の電磁波発信部のそれぞれから発信された電磁波を、前記複数の電磁波受信部のそれぞれが、前記複数の電磁波発信部のそれぞれから発信された電磁波を受信するのに要した電磁波伝播時間、又は前記複数の電磁波発信部のそれぞれが、前記複数の電磁波受信部のそれぞれから反射された電磁波を受信するのに要した電磁波伝播時間に基づいて、前記複数の電磁波発信部のそれぞれと前記第２の電子機器群に属する前記第１電子機器の備える前記複数の電磁波受信部それぞれとの距離を測定し、
前記第２電子機器位置情報算出部は、さらに前記距離計測部により測定された前記距離に基づいて、前記第２の電子機器群に属する第１電子機器の複数の電磁波受信部の前記第１の電子機器群に属する第１電子機器における第１座標系における座標値をそれぞれ算出する、請求項６に記載の点群データ取得システム。
【請求項８】
前記合成点群データ作成部は、さらに、
複数の点群データを、前記点群データ座標値変換部により変換された前記基準座標系における座標値に基づいて、重なる部分の位置合わせを行うことで、１つの合成点群データを作成する、請求項１乃至請求項７のいずれか１項に記載の点群データ取得システム。
【請求項９】
前記合成点群データ作成部は、さらに、
複数の点群データを、前記点群データ座標値変換部により変換された前記基準座標系における座標値に基づいて、重なる部分の位置合わせを行い、さらに前記基準座標系における同一座標値を有する部分を有する異なる複数の点群データを用いて、前記複数の点群データをマッチングする点群データマッチング部を備える請求項１乃至請求項８のいずれか１項に記載の点群データ取得システム。
【請求項１０】
少なくとも時間の経過に伴って形状が変化する物体を含む被写体の表面の各点までの深度情報を画素単位に有する点群データを取得する点群データ取得方法であって、
第１座標系を有する第１電子機器の位置情報計測部により、第２座標系をそれぞれ有する任意の時間に移動可能な複数の第２電子機器がそれぞれ備える複数の標識物を計測するとともに、計測した時刻を取得する位置情報計測ステップと、
前記第２電子機器それぞれに対して、前記位置情報計測ステップにおいて計測されたデータ及び計測時刻に基づいて、前記第２電子機器それぞれの備える前記複数の標識物の前記第１座標系における座標値である第２電子機器位置情報を前記計測時刻に対応付けて算出する位置情報算出ステップと、
前記複数の第２電子機器それぞれの３次元計測部により、１秒あたりに予め設定した複数の頻度となるように設定された所定の周期で、前記被写体の点群データを前記第２座標系に基づいて連続して計測し、計測した時刻を示すタイムスタンプを付した前記点群データを取得する３次元計測ステップと、
前記計測時刻に対応付けて算出された前記第２電子機器位置情報に基づいて、前記３次元計測ステップにおいて前記第２電子機器それぞれの備える前記３次元計測部により前記所定の周期で連続して計測された前記点群データにおける各点のデータの前記第２座標系における座標値を、前記タイムスタンプごとに基準座標系における座標値に変換する点群データ座標値変換ステップと、
前記３次元計測ステップにおいて前記第２電子機器それぞれの備える前記３次元計測部により前記所定の周期で連続して計測された複数の点群データを、前記点群データ座標値変換ステップにおいて変換された前記基準座標系における座標値に基づいて、前記所定の周期ごとに１つの合成点群データを作成する合成点群データ作成ステップと、
を備える点群データ取得方法。
【請求項１１】
少なくとも時間の経過に伴って形状が変化する物体を含む被写体の表面の各点までの深度情報を画素単位に有する点群データを取得する点群データ取得方法であって、
第１座標系を有する第１電子機器の複数の電磁波発信部により複数の電磁波を発信するとともに、発信した計測時刻を取得する電磁波発信ステップと、
第２座標系をそれぞれ有する任意の時間に移動可能な複数の第２電子機器それぞれの備える複数の電磁波受信部により、前記電磁波発信ステップにおいて計測時刻に対応付けて発信された前記複数の電磁波のそれぞれを受信する電磁波受信ステップと、
前記複数の電磁波受信部のそれぞれが、計測時刻に対応付けて前記複数の電磁波発信部のそれぞれから発信された電磁波を受信するのに要した電磁波伝播時間、又は前記複数の電磁波発信部のそれぞれが、前記複数の電磁波受信部のそれぞれから反射された電磁波を受信するのに要した電磁波伝播時間に基づいて、前記複数の電磁波発信部のそれぞれと前記複数の電磁波受信部それぞれとの距離を測定する距離計測ステップと、
前記距離計測ステップにおいて計測時刻に対応付けて測定された距離に基づいて、前記複数の電磁波受信部の前記第１座標系における座標値である第２電子機器位置情報をそれぞれ計測時刻に対応付けて算出する第２電子機器位置情報算出ステップと、
１秒あたりに予め設定した複数の頻度となるように設定された所定の周期で前記第２電子機器それぞれの備える３次元計測部により、前記被写体の点群データを前記第２座標系に基づいて連続して計測し、計測した時刻を示すタイムスタンプを付した前記点群データを取得する３次元計測ステップと、
前記計測時刻に対応付けて算出された前記第２電子機器位置情報に基づいて、前記３次元計測ステップにおいて前記第２電子機器それぞれの備える３次元計測部により、前記所定の周期で連続して計測された前記点群データにおける各点のデータの前記第２座標系における座標値を、前記タイムスタンプごとに基準座標系における座標値に変換する点群データ座標値変換ステップと、
前記３次元計測ステップにおいて前記第２電子機器それぞれの備える３次元計測部により、前記所定の周期で連続して計測された複数の点群データを、前記点群データ座標値変換ステップにおいて変換された前記基準座標系における座標値に基づいて、前記所定の周期ごとに１つの合成点群データを作成する合成点群データ作成ステップと、
を備える点群データ取得方法。
【請求項１２】
１つのコンピュータ又は複数のコンピュータに、請求項１０又は請求項１１に記載の方法の各ステップを実行させるためのコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の計測位置からの非接触計測による対象物の外郭を表す点群データの取得に関するシステム及びその方法に関する。
【背景技術】
【０００２】
レーザースキャナ等の三次元計測装置を使用して、計測対象物の表面の位置及び形状を含む外郭を表す点群データを計測し、この点群データから計測対象物のパノラマ画像を作成する三次元データ処理装置が知られている。この三次元データ処理装置は、点群データに基づいて、例えば三次元計測装置の計測位置を原点とする仮想的な球面に対するパノラマ画像を作成することで、点群データの表示を容易化し操作性を向上することが知られている（例えば、特許文献１参照）。
【０００３】
また、レーザースキャナ等の三次元計測装置を使用して、既設の建造物（例えば原子力発電所等のプラントや文化財等）の凹凸を有する非平面形状の計測対象物の輪郭を表す点群データを計測する場合、計測位置を様々に変えて複数の視点から当該計測対象物の輪郭を表す点群データを計測することが知られている。
複数の異なる計測位置から計測された各点群データの位置情報は、当該点群データを計測した各計測位置の位置及び方位に基づく座標系（以下「固有座標系」ともいう）により表される。
このため、これらの異なる計測位置から計測された各点群データの各固有座標系により表される位置情報を同一の座標系（以下、「基準座標系」ともいう）により表される位置情報に変換する必要がある。座標系の変換により、すべての点群データの位置情報を基準座標系により表し、同一座標系により表された各点群データを部分的に重なる部分のマッチングを行い、合成することが知られている。こうすることで、死角が減り、計測対象物の輪郭を把握する三次元データが得られる（例えば、特許文献２参照）。
ここで、マッチングとは、異なる計測位置から得られた複数の点群データにおいて、部分的に重なる部分の位置合わせを行うことで、整合性を保持しつつ統合させることを言う。
【０００４】
上記従来技術において、複数の計測位置から当該計測対象物の外郭を表す点群データを計測する場合に、各三次元計測装置の位置及び方位を予め設定しておく必要がある。このため、例えばレーザースキャナ等の三次元計測装置を予め設定した計測位置から別の計測位置に変更して設置する場合に、再度、変更先の計測位置における三次元計測装置の位置、方位等を正確に測定することが必要となる。３次元計測装置の位置の変更作業は、時間がかかり、簡単なものではない。さらには、例えば、当該計測対象物が、例えば、人の動作、物体の移動のように時間の経過に伴って、動的にその外郭が変化する非静止物体であれば、ひとつの計測装置で多方面からとらえることは、不可能である。
【０００５】
また、マッチングのために、対象構造物やその周辺に複数のマーカーを設置（例えば接着又は塗布）して、その設置した複数のマーカーを基準にして計測データの合成を行うことが知られている（例えば、特許文献３）。しかしながら、複数のマーカーを基準にするためには、マーカーを設置する作業が必要となり、そのマーカーを適切な位置に合わせる等の作業に時間がかかり、計測作業のための前準備に時間が長くなるという課題があった。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２０１０−０９７４１９号公報
【特許文献２】特開２０１３−１６１１３３号公報
【特許文献３】特開２０１０−０６６１６９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
上述した従来技術では、計測対象物の外郭を表す点群データを複数の計測位置から計測する場合に、３次元計測装置の位置、方位等を予め正確に測定しておくことが必要となっていた。
このため、３次元計測装置の位置の設定作業に時間がかかるという課題があった。
また、各３次元計測装置を予め正確に測定された位置に設置することが必要となるため、計測対象物の状況に合わせて、計測位置を即座に変更できないという課題があった。
例えば、動く人間、動く動物、動く物等、時間の経過に伴って形状が変化する物体を計測対象とした場合には、複数の計測装置を使用していても、物体の変形に伴って生じるオクルージョン（手前にある物体が背後の物体をかくす状態）による死角が発生する場合に従来技術では、計測位置を即座に変更して、欠損部分の少ない３次元画像を取得することができないという課題があった。
【０００８】
また、対象物やその周辺に複数のマーカーを設置する場合には、当該マーカーを設置する作業が必要になり、またマーカーを適切な位置に合わせる等の作業に時間がかかるという課題があった。
なお、マーカー等を対象物に設置（例えば接着又は塗布）する作業において、計測対象物に損傷を与える恐れがある場合もあるという課題もあった。
【０００９】
本発明は、複数の計測位置から計測対象物の外郭を表す点群データを取得するシステム又は方法において、計測対象物を計測するために、３次元計測装置を移動させたときに３次元計測装置の位置、方位をリアルタイムに取得することで、当該計測対象物の輪郭を表す複数の点群データを３６０度の角度から取得し、欠損部分の少ない３次元画像を取得するシステムを提供することを目的とする。
【課題を解決するための手段】
【００１０】
（１）本願発明は、第１電子機器と１つ以上の移動可能な第２電子機器とを有する第１の電子機器群を備える、被写体の表面の各点までの深度情報を画素単位に有する点群データを取得する点群データ取得システムであって、前記第２電子機器は、それぞれが、第２座標系（例えば、後述の「第２座標系」）を有し、第２通信部と、複数の、視覚的特徴を有する所定の標識物であって、前記複数の標識物はそれぞれ第２座標系において線形独立となる位置にそれぞれ設けられる、前記複数の標識物（例えば、後述の「マーカー２３」）と、被写体の点群データを前記第２座標系に基づいて計測する３次元計測部と、を備え、前記第１電子機器は、第１座標系（例えば、後述の「第１計測座標系」）を有し、被写体を前記第１座標系に基づいて計測する位置情報計測部（例えば、後述の「位置情報計測部１２」）を備え、前記点群データ取得システムは、さらに基準座標系（例えば、後述の「基準座標系」）を有し、前記位置情報計測部により計測されたデータに基づいて、前記第２電子機器の備える前記複数の標識物の前記第１座標系における座標値である第２電子機器位置情報を算出するマーカー位置情報算出部（例えば、後述の「マーカー位置情報算出部３１」）と、前記第２電子機器位置情報に基づいて、前記第２電子機器の備える前記３次元計測部により計測された前記点群データにおける各点のデータの前記第２座標系における座標値を、前記基準座標系における座標値に変換する点群データ座標値変換部と、前記第２電子機器の備える前記３次元計測部により計測された複数の点群データを、前記点群データ座標値変換部により変換された前記基準座標系における座標値に基づいて、１つの合成点群データを作成する合成点群データ作成部と、を備える、点群データ取得システムに関する。
【００１１】
（２）前記第２電子機器の備える視覚的特徴を有する所定の標識物は、それぞれ表面が予め設定された色で着色された、所定の半径の球であって、前記球の中心がそれぞれ予め設定された三角形の頂点を形成するとともに、前記第２電子機器を一意的に識別できるように、前記複数の標識物の視覚的特徴を前記第２電子機器ごとに異なるように構成されたことを特徴とする、（１）に記載の点群データ取得システム。
【００１２】
（３）前記第１の電子機器群は、２つ以上の第１電子機器を備え、１つの第１電子機器の備える位置情報計測部、及び別の第１電子機器の備える位置情報計測部が、それぞれの備える固有の第１座標系に基づいて、前記第２電子機器の備える前記複数の標識物を計測することを特徴とする（１）又は（２）に記載の点群データ取得システム。
【００１３】
（４）前記点群データ取得システムは、さらに、複数の、視覚的特徴を有する所定の標識物（例えば、後述の「基準座標系マーカー１０Ａ」）であって、前記複数の標識物はそれぞれ基準座標系において線形独立となる座標位置にそれぞれ設けられる、前記複数の標識物を備え、前記点群データ座標値変換部は、前記マーカー位置情報算出部により算出された、前記複数の、視覚的特徴を有する所定の標識物の前記第１座標系における座標値、及び前記複数の、視覚的特徴を有する所定の標識物の前記基準座標系における座標値に基づいて、予め前記第１座標系における座標値を前記基準座標系における座標値に変換する基準座標変換行列を生成し、前記点群データ座標値変換部は、前記第２電子機器位置情報、前記点群データ、及び前記基準座標変換行列に基づいて、前記点群データにおける各点のデータの前記第２座標系における座標値を、前記基準座標系における座標値に変換する、（１）〜（３）のいずれかに記載の点群データ取得システム。
【００１４】
（５）前記点群データ取得システムは、さらに、前記第１の電子機器群に属さない、１つの第１電子機器と１つ以上の第２電子機器を備える第２の電子機器群を有し、前記第２の電子機器群に属する第１電子機器は、さらに、複数の、視覚的特徴を有する所定の標識物であって、前記複数の標識物はそれぞれ第１座標系において線形独立となる位置にそれぞれ設けられる、前記複数の標識物を備え、前記第１の電子機器群に属する前記第１電子機器の前記位置情報計測部は、第１座標系に基づいて前記第２の電子機器群に属する前記第１電子機器の備える前記複数の標識物を計測し、前記マーカー位置情報算出部は、前記第２の電子機器群に属する前記第１電子機器の備える前記複数の標識物の前記第１座標系における座標値である、前記第２の電子機器群に属する前記第１電子機器の位置情報を算出する、（１）〜（４）のいずれかに記載の点群データ取得システム。
【００１５】
（６）本願発明は、第１電子機器と１つ以上の移動可能な第２電子機器とを有する第１の電子機器群を備える、被写体の表面の各点までの深度情報を画素単位に有する点群データを取得する点群データ取得システムであって、前記第１電子機器は、第１座標系（例えば、後述の「第１座標系」）を有し、電磁波を発信する、複数の電磁波発信部であって、前記複数の電磁波発信部は、それぞれ第１座標系の線形独立である任意の位置に設けられる、前記複数の電磁波発信部と、を備え、前記第２電子機器は、それぞれが固有の第２座標系（例えば、後述の「第２座標系」）を有し、電磁波を受信又は反射する、複数の電磁波受信部であって、前記複数の電磁波受信部は、それぞれ前記第２座標系の線形独立である任意の位置に設けられる、前記複数の電磁波受信部と、前記点群データを前記第２座標系に基づいて計測する３次元計測部と、を備え、前記点群データ取得システムは、さらに基準座標系（例えば、後述の「基準座標系」）を有し、前記複数の電磁波受信部のそれぞれが、前記複数の電磁波発信部のそれぞれから発信された電磁波を受信するのに要した電磁波伝播時間、又は前記複数の電磁波発信部のそれぞれが、前記複数の電磁波受信部のそれぞれから反射された電磁波を受信するのに要した電磁波伝播時間に基づいて、前記複数の電磁波発信部のそれぞれと前記複数の電磁波受信部それぞれとの距離を測定する距離計測部と、前記距離計測部により測定された距離に基づいて、前記複数の電磁波受信部の前記第１座標系における座標値である第２電子機器位置情報をそれぞれ算出する第２電子機器位置情報算出部と、前記第２電子機器位置情報に基づいて、前記第２電子機器の備える前記３次元計測部により計測された前記点群データにおける各点のデータの前記第２座標系における座標値を、前記基準座標系における座標値に変換する点群データ座標値変換部と、前記第２電子機器の備える前記３次元計測部により計測された複数の点群データを、前記点群データ座標値変換部により変換された前記基準座標系における座標値に基づいて、１つの合成点群データを作成する合成点群データ作成部と、を備える、点群データ取得システムに関する。
【００１６】
（７）前記点群データ取得システムは、さらに、前記第１の電子機器群に属さない、１つの第１電子機器と１つ以上の第２電子機器を備える第２の電子機器群を有し、前記第２の電子機器群に属する第１電子機器は、さらに、電磁波を受信又は反射する、複数の電磁波受信部を備え、前記距離計測部は、さらに前記第１の電子機器群に属する前記第１電子機器の前記複数の電磁波発信部のそれぞれから発信された電磁波を、前記複数の電磁波受信部のそれぞれが、前記複数の電磁波発信部のそれぞれから発信された電磁波を受信するのに要した電磁波伝播時間、又は前記複数の電磁波発信部のそれぞれが、前記複数の電磁波受信部のそれぞれから反射された電磁波を受信するのに要した電磁波伝播時間に基づいて、前記複数の電磁波発信部のそれぞれと前記第２の電子機器群に属する前記第１電子機器の備える前記複数の電磁波受信部それぞれとの距離を測定し、前記第２電子機器位置情報算出部は、さらに前記距離計測部により測定された前記距離に基づいて、前記第２の電子機器群に属する第１電子機器の複数の電磁波受信部の前記第１の電子機器群に属する第１電子機器における第１座標系における座標値をそれぞれ算出する、（６）に記載の点群データ取得システム。
【００１７】
（８）前記合成点群データ作成部は、さらに、複数の点群データを、前記点群データ座標値変換部により変換された前記基準座標系における座標値に基づいて、重なる部分の位置合わせを行うことで、１つの合成点群データを作成する、（１）〜（７）のいずれかに記載の点群データ取得システム。
【００１８】
（９）前記合成点群データ作成部は、さらに、複数の点群データを、前記点群データ座標値変換部により変換された前記基準座標系における座標値に基づいて、重なる部分の位置合わせを行い、さらに前記基準座標系における同一座標値を有する部分を有する異なる複数の点群データを用いて、前記複数の点群データをマッチングする点群データマッチング部を備える、（１）〜（８）のいずれかに記載の点群データ取得システム。
【００１９】
（１０）本願発明は、被写体の表面の各点までの深度情報を画素単位に有する点群データを取得する点群データ取得方法であって、第１座標系を有する第１電子機器の位置情報計測部により、第２座標系をそれぞれ有する１つ以上の第２電子機器がそれぞれ備える複数の標識物を計測する位置情報計測ステップと、前記位置情報計測ステップにより計測されたデータに基づいて、前記第２電子機器の備える前記複数の標識物の前記第１座標系における座標値である第２電子機器位置情報を算出する位置情報算出ステップと、前記１つ以上の第２電子機器の３次元計測部により、点群データを前記第２座標系に基づいて計測する３次元計測ステップと、前記第２電子機器位置情報に基づいて、前記３次元計測ステップにより計測された前記点群データにおける各点のデータの前記第２座標系における座標値を、基準座標系における座標値に変換する点群データ座標値変換ステップと、前記３次元計測ステップにより計測された複数の点群データを、前記点群データ座標値変換ステップにより変換された前記基準座標系における座標値に基づいて、１つの合成点群データを作成する合成点群データ作成ステップと、を備える点群データ取得方法に関する。
【００２０】
（１１）本願発明は、被写体の表面の各点までの深度情報を画素単位に有する点群データを取得する点群データ取得方法であって、第１座標系を有する第１電子機器の複数の電磁波発信部により複数の電磁波を発信する電磁波発信ステップと、第２座標系をそれぞれ有する１つ以上の第２電子機器の複数の電磁波受信部により、前記電磁波発信ステップにより発信された前記複数の電磁波のそれぞれを受信する超音波受信ステップと、前記複数の電磁波受信部のそれぞれが、前記複数の電磁波発信部のそれぞれから発信された電磁波を受信するのに要した電磁波伝播時間、又は前記複数の電磁波発信部のそれぞれが、前記複数の電磁波受信部のそれぞれから反射された電磁波を受信するのに要した電磁波伝播時間に基づいて、前記複数の電磁波発信部のそれぞれと前記複数の電磁波受信部それぞれとの距離を測定する距離計測ステップと、前記距離計測ステップにより測定された距離に基づいて、前記複数の電磁波受信部の前記第１座標系における座標値である第２電子機器位置情報をそれぞれ算出する第２電子機器位置情報算出ステップと、前記１つ以上の第２電子機器の３次元計測部により、点群データを前記第２座標系に基づいて計測する３次元計測ステップと、前記第２電子機器位置情報に基づいて、前記３次元計測ステップにより計測された前記点群データにおける各点のデータの前記第２座標系における座標値を、基準座標系における座標値に変換する点群データ座標値変換ステップと、前記３次元計測ステップにより計測された複数の点群データを、前記点群データ座標値変換ステップにより変換された前記基準座標系における座標値に基づいて、１つの合成点群データを作成する合成点群データ作成ステップと、を備える点群データ取得方法に関する。
【００２１】
（１２）本願発明は、コンピュータに、（１０）又は（１１）に記載の各ステップを実行させるためのコンピュータプログラムに関する。
【発明の効果】
【００２２】
本発明の点群データ取得システムによれば、対象物を計測するために、３次元計測装置の位置、方位を予め設定した場所に設置することなく、３次元計測装置の計測位置を即座に変更することが必要になった場合に、３次元計測装置を移動させるとともに、３次元計測装置の位置及び方位をリアルタイムに取得することができる。
そうすることで、例えば、動く人間、動く動物、動く物等、時間の経過に伴って形状が変化する物体を計測対象とした場合には、物体の変形に伴って生じるオクルージョン（手前にある物体が背後の物体をかくす状態）による死角が発生するときに、欠損部分を極力少なくするために、３次元計測装置の計測位置を即座に変更して対応することができ、３６０度の角度から計測対象物の輪郭を表す複数の点群データを取得することができる。
【図面の簡単な説明】
【００２３】
【図１Ａ】本発明の第１実施形態に係る点群データ取得システムの一例を示す概略図である。概略的に示すブロック図である。
【図１Ｂ】本発明の第１実施形態に係る第１電子機器の概略ブロック図である。
【図１Ｃ】本発明の第１実施形態に係る第２電子機器の概略ブロック図である。
【図１Ｄ】本発明の第１実施形態に係る座標系算出部の概略ブロック図である。
【図１Ｅ】本発明の第１実施形態に係る３次元データ処理部の概略ブロック図である。
【図１Ｆ】本発明の第１実施形態に係る３次元データ表示制御部の概略ブロック図である。
【図１Ｇ】本発明の第１実施形態に係る３次元データ保存部の概略ブロック図である。
【図２】本発明の第１実施形態に係るマーカーの一例を示す概略図である。
【図３】本発明の第１実施形態に係るキャリブレーションを説明する各座標系間の変換行列の関連を示す概略図である。
【図４】本発明の一実施形態における一連の処理の流れを示す図である。
【図５Ａ】本発明の第２実施形態に係る第１電子機器の概略ブロック図である。
【図５Ｂ】本発明の第２実施形態に係る第２電子機器の概略ブロック図である。
【図５Ｃ】本発明の第２実施形態に係る座標系算出部の概略ブロック図である。
【図６】本発明の第２実施形態に係る電波発信部と電波受信部間における電波伝播時間とそれに基づく直線距離の関係を示す概略図である。
【図７】第１電子機器及び１つ以上の第２電子機器から構成される電子機器群を階層的に配置して、階層化を示す概略図である。
【発明を実施するための形態】
【００２４】
以下、本発明の点群データ取得システムの好ましい実施形態について、図面を参照しながら説明する。ここで、下記の実施形態はいずれも本発明の一例であり、これに限定されるものではない。
【００２５】
［第１実施形態］
図１Ａに示すように、点群データ取得システム１００は、第１電子機器１０と、１つ以上の移動可能な第２電子機器２０と、座標系算出部３０と、３次元データ処理部４０と、３次元データ表示制御部５０と、３次元データ保存部６０と、を含んで構成される。これらは、有線及び／又は無線ネットワーク等を介して接続されている。
後述するように、点群データ取得システム１００は、システムに固有の座標系（以下、「基準座標系」又は「グローバル座標系」ともいう）を有するように構成される。例えば、後述する「基準座標系マーカー」を用いて、キャリブレーションを行うことで、基準座標系を設定することができる。詳細については後述する。
【００２６】
＜第１電子機器１０＞
第１電子機器１０は、有線及び／又は無線通信機能を備える、例えばノートパソコンを含むパソコン、あるいはＣＰＵ、ＤＳＰ等の演算チップを搭載した演算ボードで構成することができる。
【００２７】
第１電子機器１０は、第１通信部１１と、被写体の表面の各点を画素単位で例えばＲＧＢの色情報とともに位置情報を計測する位置情報計測部１２と、を備える。第１電子機器１０の備える位置情報計測部１２は、当該位置情報計測部１２に固有の座標系（以下、「第１計測座標系」という）を有する。
また、第１電子機器１０は、さらに時刻同期部１５を備えることができる。以下、第１電子機器１０の備える上記各機能部を総括して「第１機能部」ともいう。
【００２８】
第１電子機器１０には、当該第１電子機器１０を、前述した第１通信部１１、位置情報計測部１２等、第１機能部として機能させるためのプログラムをインストールし、第１通信部１１、位置情報計測部１２等の各機能部として機能させるようにすることができる。なお、当該プログラムは、第１電子機器１０に予めインストールされていてもよい。また、当該プログラムは、コンピュータ読み取り可能な記録媒体から必要に応じてインストール又は予め設定されたサーバから必要に応じて適宜ダウンロードされてもよい。以下、「コンピュータで読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。
【００２９】
＜第２電子機器２０＞
第２電子機器２０は、有線及び／又は無線通信機能と、を備える、例えばノートパソコンを含むパソコン、タブレット端末等の移動可能な電子機器、あるいはＣＰＵ、ＤＳＰ等の演算チップを搭載した演算ボードで構成することができる。
【００３０】
第２電子機器２０は、第２通信部２１と、３次元計測部２２と、マーカー２３と、を備える。３次元計測部２２は、当該３次元計測部２２に固有の座標系（以下、「第２座標系」という）を有する。同様に、マーカー２３は、当該マーカー２３に固有の座標系（以下、「マーカー座標系」という）を有する。詳細については後述する。
なお、マーカー２３を「位置情報ユニット」と称することができる。
また、第２電子機器２０は、さらに、時刻同期部２５を備えることができる。
以下、第２電子機器２０の備える上記各機能部を総括して「第２機能部」ともいう。
【００３１】
第２電子機器２０には、当該第２電子機器２０を、第２通信部２１、３次元計測部２２等、第２機能部として機能させるためのプログラムをインストールすることで、第２電子機器２０を第２通信部２１、３次元計測部２２等の各機能部として機能させるようにすることができる。なお、当該プログラムは、予めインストールされていてもよい。また、当該プログラムは、コンピュータ読み取り可能な記録媒体から必要に応じてインストール又は予め設定されたサーバから必要に応じて適宜ダウンロードされてもよい。
【００３２】
＜座標系算出部３０＞
座標系算出部３０は、マーカー位置情報算出部３１と、キャリブレーション処理部３２と、を備える。
【００３３】
＜３次元データ処理部４０＞
３次元データ処理部４０は、点群データ座標値変換部４１と、合成点群データ作成部４２と、を備える。
合成点群データ作成部４２は、さらに、第１合成点群データ作成部４２１と、点群データマッチング部４２２と、を備える。
【００３４】
＜３次元データ表示制御部５０＞
３次元データ表示制御部５０は、視点操作部５１と、データ表示方式選択部５２と、データ再生・保存指示部５３と、３次元データ表示処理部５４と、を備える。
【００３５】
＜３次元データ保存部６０＞
３次元データ保存部６０は、源点群データ記憶部６１と、第２電子機器位置情報記憶部６２と、を備える。
【００３６】
座標系算出部３０、３次元データ処理部４０、３次元データ表示制御部５０、及び３次元データ保存部６０は、有線及び／又は無線通信機能を備える、例えば、１台又は複数のパソコン、タブレット端末、スマートフォン、サーバ等の電子機器、あるいはＣＰＵ、ＤＳＰ等の演算チップを搭載した演算ボードで構成することができる。
また、座標系算出部３０、３次元データ処理部４０、３次元データ表示制御部５０、及び３次元データ保存部６０をそれぞれ別のパソコン、タブレット端末、スマートフォン、サーバ等の電子機器により構成された、分散システムとすることができる。
また、座標系算出部３０、３次元データ処理部４０、３次元データ表示制御部５０、及び３次元データ保存部６０の一部又は全てを第１電子機器１０、又は第２電子機器２０により構成することもできる。
また、座標系算出部３０、３次元データ処理部４０、３次元データ表示制御部５０、及び３次元データ保存部６０をクラウド上の仮想サーバ上においてもよい。
【００３７】
座標系算出部３０、３次元データ処理部４０、３次元データ表示制御部５０、及び３次元データ保存部６０を構成する電子機器には、当該電子機器を、座標系算出部３０、３次元データ処理部４０、３次元データ表示制御部５０、及び３次元データ保存部６０として機能させるためのプログラムをインストールすることで、当該電子機器を座標系算出部３０、３次元データ処理部４０、３次元データ表示制御部５０、及び３次元データ保存部６０として機能させるようにすることができる。なお、当該プログラムは、予め当該電子機器にインストールされていてもよい。また、当該プログラムは、コンピュータ読み取り可能な記録媒体から必要に応じてインストール又は予め設定されたサーバから必要に応じて適宜ダウンロードされてもよい。
【００３８】
次に、点群データ取得システム１００の備える各電子機器等について、詳細に説明する。
最初に、第２電子機器２０について説明する。
【００３９】
＜第２電子機器２０について＞
図１Ｃに示すように、第２電子機器２０は、第２通信部２１と、３次元計測部２２、マーカー２３と、を備える。また、第２電子機器２０は、さらに、時刻同期部２５を備えることができる。
【００４０】
[第２通信部２１]
有線、又は無線回線による通信モジュール。第１電子機器１０の第１通信部１１と対向し、第１電子機器１０、第２電子機器２０間の制御信号、各種データをやり取りする。
【００４１】
［３次元計測部２２］
３次元デジタイザ、レーザースキャナ等と同様の３次元計測手段として、「深度センサ」がある。「深度センサ」には主に、（１）Ｔｉｍｅ−ｏｆ−ｆｌｉｇｈｔ方式、（２）三角測量方式、（３）Ｄｅｐｔｈ  ｏｆ  Ｄｅｆｏｒｃｅ方式等があるが、（１）のＴｉｍｅ−ｏｆ−ｆｌｉｇｈｔ方式は、レーザや赤外線等の光を計測対象物に対して照射して反射させ、センサに届くまでの飛行時間をもとに計算し、計測対象物表面の各点までの深度センサのからの距離等（以下、「深度」ともいう）を瞬時に得ることができるセンサである。深度センサは計測対象物各点の深度を画素単位で把握することができる。
ここで、計測対象物とは、例えば、既設の建造物（例えば原子力発電所等のプラントや文化財等）の凹凸を有する非平面形状の計測対象物のみならず、動く人間、動く動物、動く物等、時間の経過に伴って形状が変化する物体を含むことができる。
第２電子機器２０の有する３次元計測部２２は、深度センサを、例えば通常の画像カメラと組み合わせて構成することができる。深度センサを通常の画像カメラと組み合わせて構成される３次元計測部２２を、以下「深度カメラ」ともいう。
３次元計測部２２は、例えば深度カメラを用いて計測対象物を計測することで、計測対象物表面の各点を画素単位で例えばＲＧＢの色情報とともに深度カメラの備える３次元の座標系における座標値（以下「深度情報」ともいう）を併せ持つ、計測対象物の画像データを生成することができる。
【００４２】
３次元計測部２２は、固有の基準位置に基づく３次元の固有座標系（以下、「第２座標系」ともいう）を備えている。３次元計測部２２により計測された計測対象物表面の各点（画素単位）の位置は、第２座標系により、３次元座標（ｘ、ｙ、ｚ）により表される。
したがって、３次元計測部２２により計測する場合には、計測対象物表面の各点（画素単位）の位置を表す３次元座標（ｘ、ｙ、ｚ）と、その３次元座標位置における計測対象物表面のＲＧＢの色情報と、からなる点データの集合を計測データとして取得することができる。
【００４３】
３次元計測部２２は、第２電子機器２０を移動させることで、３次元空間内を移動することができる。ユーザ（例えば、計測作業者、カメラマン）は、計測対象物の状況に応じて、第２電子機器２０を移動させることで、３次元計測部２２の計測位置を変更して、計測対象物の計測を行うことができる。なお、第２電子機器２０の移動については、例えば、第２電子機器２０を搭載した移動装置を自走式又は遠隔操作により、移動させることもできる。例えば、第２電子機器２０を搭載したクアッドコプターを遠隔操作することにより、移動させることもできる。
【００４４】
３次元計測部２２は、計測対象物表面の各点の情報を画素単位で取得することができる。また、３次元計測部２２は、計測対象物の表面の各点の情報を画素単位で、予め設定した時間間隔で取得することができる。例えば、時間間隔として１／３０秒を設定することで、３次元計測部２２は、１秒間に３０コマの画像データ（３０ｆｐｓ）を取得することができる。３次元計測部２２により計測された各計測データには、それぞれ計測時刻となるタイムスタンプを付すことができ、当該タイムスタンプに基づいて計測データ間で関連付けが可能となっている。
すなわち、Ｎ個（１≦Ｎ）の第２電子機器２０及び３次元計測部２２をインデックス１、２、３、・・・Ｎにより識別する場合、ｔを計測時刻とすると、インデックス（識別番号）ｎの３次元計測部２２（以下「３次元計測部２２ｎ」ともいう）により取得した点群データをＱｎ（ｔ）と表現することができる。なお、インデックス（識別番号）ｎの第２電子機器２０は、第２電子機器２０ｎという。
【００４５】
＜３次元計測部２２による測定時刻の同期制御＞
例えば、第２電子機器２０が複数個ある場合（すなわち３次元計測部２２が複数個ある場合）、例えば、第１電子機器１０又は座標系算出部３０のトリガーによる同期制御により、複数個あるすべての３次元計測部２２に対して計測動作を同期させることにより、同時刻に計測対象物を計測させるように構成することができる。トリガーとしては、例えば、第１電子機器１０又は座標系算出部３０から３次元計測部２２ｎに対して計測指示信号を出すことができる。
【００４６】
以下、Ｎは１以上（すなわち、第２電子機器２０は１以上）とする。
ここで、各３次元計測部２２ｎにより計測対象物の画像データを予め設定した時間間隔（例えば、１秒間に３０コマ）で取得する場合、当該点群データは、第２電子機器２０ｎの備える３次元計測部２２ｎ（１≦ｎ≦Ｎ）毎に、１／３０秒単位の計測時刻（ｔ）毎に生成され、通信ネットワークを介して、３次元データ保存部６０の源点群データ記憶部６１にリアルタイムに記憶される。
すなわち、源点群データ記憶部６１には、すべての３次元計測部２２ｎにより計測される、点群データの集合｛Ｑｎ（ｔ）｝（１≦ｎ≦Ｎ）（ｔ：タイムスタンプ）が記憶される。
点群データＱｎ（ｔ）は、計測対象物の表面の各点を、例えばデータ形式｛３次元計測部２２の識別番号ｎ、各点の第２座標系における座標値（ｘ、ｙ、ｚ）、点群取得時刻ｔ、Ｒ（赤成分）、Ｇ（緑成分）、Ｂ（青成分）等｝で表した点データの集合であるといえる。
【００４７】
［マーカー２３］
第２電子機器２０は、３次元計測部２２の備える第２座標系において、線形独立となる座標値が既知である位置に設けられる複数の視覚的特徴を有する所定の標識物（以下、「マーカー２３」という）を備える。なお、インデックス（識別番号）ｎの第２電子機器２０ｎの備えるマーカー２３をマーカー２３（ｎ）という。
第１実施形態は、後述するように、座標系算出部３０の備えるマーカー位置情報算出部３１が位置情報計測部１２による計測情報に基づいて、第２電子機器２０ｎの備えるマーカー２３（ｎ）を検出し、検出したマーカー２３（ｎ）の構成要素（マーカー２３（ｎ）の頂点等）の第１計測座標系における座標値（以下「マーカー位置情報（ｎ）」という）を算出する。
そうすることで、３次元データ処理部４０の備える点群データ座標値変換部４１は、第２電子機器２０ｎの備えるマーカー座標系における座標値を第１電子機器１０の備える位置情報計測部１２における第１計測座標系又は第１計測座標系とは独立した基準座標系（グローバル座標系）における座標値に変換することができる。
【００４８】
なお、測定時刻ｔにおける位置情報計測部１２による計測情報に基づいて算出されるマーカー位置情報（ｎ）を算出する際に用いた位置情報計測部１２による測定時刻ｔと、点群データＱｎ（ｔ）の測定時刻ｔと、を突合せすることで、点群データＱｎ（ｔ）が計測されたときの第２電子機器２０ｎのマーカー位置情報（ｎ）を特定することができる。
測定時刻が必ずしも一致しない場合には、例えば当該時刻の近傍のマーカー位置情報（ｎ）により、線形補完することができる。
【００４９】
図２に、マーカー２３の一例を示す。図２に示すように、マーカー２３は、例えば、各頂点を同一の半径の球体とした３辺を有する三角形の形状をした構造体とすることができる。より具体的には、三角形の各頂点は、各球体の中心と一致するように構成されることが好ましい。また、三角形は直角三角形とすることが好ましいが、直角三角形に限定されるものではなく、任意の三角形でもよい。
球体とすることでどの角度から計測されても計測画像データが円となることから、マーカー位置情報算出部３１は、位置情報計測部１２による計測結果に基づいて、球体を検出することが容易となる。なお、球体の半径、球体の表面の色、及び三角形の辺の長さは、適宜設定することができる。
なお、第１実施形態においては、マーカー２３を所定の色で着色された所定の半径の球の形状をしたものと仮定するが、視覚的特徴を有する所定の標識物であるマーカーの形状は、球体に限定されるものではない。位置情報計測部１２により、マーカーを計測することができ、マーカー位置情報算出部３１により当該マーカーを検出することができる任意の形状であればよい。例えば、色、柄、模様、特別な形状等により識別できるようにしてもよい。
【００５０】
例えば、マーカー２３を構成する三角形を全て相似形となる直角三角形とし、辺の長さを各第２電子機器２０ｎ毎に異なる長さにすることで、辺の長さにより、当該マーカー２３を有する第２電子機器２０を識別するように構成することができる。また、球体の半径を各第２電子機器２０毎に異なる長さにすることで、球体の半径により、当該マーカー２３を有する第２電子機器２０ｎを識別するように構成することができる。また、球体の表面の色を各第２電子機器２０毎に異なる色にすることで、球体の表面の色により、当該マーカー２３を有する第２電子機器２０ｎを識別するように構成することができる。
なお、直角三角形を作成するに際して、３辺の長さの比を３：４：５とすることが、実用上、また直角三角形の作成上、好ましい。マーカー２３の一例として、例えば、球の半径を５ｃｍ、球の表面の色をＲＥＤ、３辺の長さをそれぞれ３ｉｃｍ、４ｉｃｍ、５ｉｃｍ（ｉは、１以上の任意の整数）とすることができる。
【００５１】
［マーカー座標系］
ここで、マーカー２３によるマーカー座標系とは、一例として、例えば、マーカー２３を直角三角形とした場合、当該直角三角形の直角に対応する頂点を原点とする座標系を意味する。
より具体的には、マーカー２３は、前述したように各頂点を同一の半径の球体とした３辺を有する三角形の形状をした構造体と仮定する。例えば、マーカー２３の３辺の長さをそれぞれ３ｉｃｍ、４ｉｃｍ、５ｉｃｍ（ｉは、１以上の任意の整数）となる直角三角形を仮定する。
マーカー２３における直角三角形の直角に対応する頂点をＣ´０、辺の長さが４ｉｃｍとなる辺のＣ´０以外の頂点をＣ´１、辺の長さが３ｉｃｍとなる辺のＣ´０以外の頂点をＣ´２とする。Ｃ´０、Ｃ´１、Ｃ´２は、それぞれ、各球体の中心と一致する。
そして、頂点Ｃ´０を始点として頂点Ｃ´１を終点とするベクトル（以下「ベクトルＣ´０Ｃ´１」という）と、頂点Ｃ´０を始点として頂点Ｃ´２を終点とするベクトル（以下「ベクトルＣ´０Ｃ´２」という）との外積を頂点Ｃ´０を始点として表現した場合における終点をＣ´３とする。
【００５２】
ここで、Ｃ´０を原点、ベクトルＣ´０Ｃ´１をＸ軸、ベクトルＣ´０Ｃ´２をＹ軸、ベクトルＣ´０Ｃ´３をＺ軸とする、マーカー座標系を規定することができる。
具体的には、各点Ｃ´０、Ｃ´１、Ｃ´２、Ｃ´３のマーカー座標系における座標値は、次のとおりである。
Ｃ´０＝（０，０，０）
Ｃ´１＝（４ｉ，０，０）
Ｃ´２＝（０，３ｉ，０）
Ｃ´３＝（０，０，１２ｉ２）
【００５３】
なお、マーカー２３が直角三角形でない場合においても、例えば、Ｃ´０を原点として、ベクトルＣ´０Ｃ´１をＸ軸とし、三角形をＸＹ平面とすることで、マーカー２３によりマーカー座標系を規定することができる。
なお、上記座標値は、一例であって、マーカー２３は任意の三角形とすることができ、その場合の座標値は、当該三角形に基づいて算出される。
【００５４】
［時刻同期部２５］
時刻同期部２５は、例えば、ＩＥＥＥ１５８８に規定されたプロトコルを用いて、第２電子機器２０の内部時刻を基準時刻に同期させる。後述するように、第１電子機器１０、及び／又は座標系算出部３０も、同様に時刻同期部１５を有することができる。そうすることで、例えば、ＩＥＥＥ１５８８に規定されたプロトコルを用いて、第１電子機器１０及び／又は座標系算出部３０の内部時刻を基準時刻に同期させる。
そうすることで、第１電子機器１０及び／又は座標系算出部３０の内部時刻並びに１つ以上の第２電子機器２０の各内部時刻を、同期させることができる。
【００５５】
＜第１電子機器１０について＞
図１Ｂに示すように、第１電子機器１０は、第１通信部１１と、被写体の表面の各点を画素単位で例えばＲＧＢの色情報とともに位置情報を計測する位置情報計測部１２と、を備える。また、第１電子機器１０は、さらに時刻同期部１５を備えるようにしてもよい。
【００５６】
[第１通信部１１]
有線、又は無線回線による通信モジュール。第２電子機器２０の第２通信部２１と対向し、第１電子機器１０、第２電子機器２０間の制御信号、各種データをやり取りする。
【００５７】
［位置情報計測部１２］
位置情報計測部１２は、第２電子機器２０の備える３次元計測部２２と同様に被写体の表面の各点を画素単位で例えばＲＧＢの色情報とともに位置情報を計測する。位置情報計測部１２は、固有の座標系（第１計測座標系）を有する。
位置情報計測部１２は、第２電子機器２０の備えるマーカー２３を検出するための計測情報を、後述するマーカー位置情報算出部３１に提供するものである。
したがって、位置情報計測部１２は、前述したように第１計測座標系を有することから、例えば、位置情報計測部１２が１台しかない場合、当該位置情報計測部１２の有する第１計測座標系を基準座標系とすることができる。しかしながら、位置情報計測部１２が１台しかない場合、何らかの物体が、位置情報計測部１２と第２電子機器２０ｎのマーカー２３との間に介在した場合、位置情報計測部１２から第２電子機器２０ｎのマーカー２３の座標値を計測できなくなる可能性がある。
第１実施形態においては、このようなオクルージョンによる死角を生じないように、複数の位置情報計測部１２を備えることが好ましい。
その場合、位置情報計測部１２を備える第１電子機器１０を複数備えるようにシステムを構成してもよい。また、第１電子機器１０が複数の位置情報計測部１２を備えるように構成してもよい。
ここで、ｉ（１≦ｉ≦Ｍ）（Ｍは位置情報計測部１２の個数）を位置情報計測部１２の識別番号とした場合、インデックス（識別番号）ｉの位置情報計測部１２は、位置情報計測部１２（ｉ）という。同様に、位置情報計測部１２（ｉ）の備える固有の座標系（第１計測座標系）を第１計測座標系（ｉ）という。
なお、第１実施形態では、位置情報計測部１２の個数Ｍは、１以上（１≦Ｍ）として、１個の場合を含むものとする。
【００５８】
［時刻同期部１５］
時刻同期部１５は、例えば、ＩＥＥＥ１５８８に規定されたプロトコルを用いて、第１電子機器１０の内部時刻を基準時刻に同期させる。前述したように、第２電子機器２０も、同様に時刻同期部２５を有し、例えば、ＩＥＥＥ１５８８に規定されたプロトコルを用いて、第１電子機器１０の内部時刻を基準時刻に同期させる。
そうすることで、第１電子機器１０の内部時刻及び座標系算出部３０、複数の第２電子機器２０の各内部時刻の同期をとることができる。
【００５９】
＜座標系算出部３０＞
図１Ｄに示すように、座標系算出部３０は、マーカー位置情報算出部３１と、キャリブレーション処理部３２と、を備える。また、座標系算出部３０は、時刻同期部３３を備えることができる。
【００６０】
［マーカー位置情報算出部３１］
マーカー位置情報算出部３１は、位置情報計測部１２により計測された計測情報に基づいて、第２電子機器２０の備えるマーカー２３を検出し、検出したマーカー２３の構成要素（マーカー２３の頂点等）の第１計測座標系における座標値（マーカー位置情報）を算出する。
【００６１】
ここで、ｉ（１≦ｉ≦Ｍ）を位置情報計測部１２の識別番号、ｎ（１≦ｎ≦Ｎ）を第２電子機器２０（３次元計測部２２）の識別番号、ｔを計測時刻、Ｌを第２電子機器２０のマーカー位置情報を表す場合、時刻ｔにおいて位置情報計測部１２（ｉ）により計測される第２電子機器２０ｎのマーカー位置情報をｉＬｎ（ｔ）と表現することができる。
したがって、ｉＬｎ（ｔ）は、時刻ｔにおける、第１計測座標系（ｉ）による、第２電子機器２０ｎのマーカー位置情報を意味する。
【００６２】
マーカー位置情報算出部３１は、第２電子機器２０ｎ毎に、第１計測座標系（ｉ）による、第２電子機器２０ｎのマーカー位置情報ｉＬｎ（ｔ）を、ネットワークを介して３次元データ保存部６０の第２電子機器位置情報記憶部６２にリアルタイムに記憶するようにしてもよい。
そうすることで、第２電子機器位置情報記憶部６２には、すべてのマーカー構成要素位置情報の集合｛ｉＬｎ（ｔ）｝（１≦ｎ≦Ｎ，１≦ｉ≦Ｍ，ｔ：タイムスタンプ）がリアルタイムに記憶される。
【００６３】
［キャリブレーション処理部３２］
キャリブレーション処理部３２は、計測対象物の点群データ取得に際して、予めキャリブレーションを行うことで、第１計測座標系における座標値を基準座標系における座標値に変換するための基準座標変換行列ＴＷ、及び第２電子機器２０の３次元計測部２２の備える第２座標系における座標値を第２電子機器２０の備えるマーカー座標系における座標値に変換するためのマーカー座標変換行列Ｔｈを算出する。詳細については、後述する。
【００６４】
［時刻同期部３３］
時刻同期部３３は、例えば、ＩＥＥＥ１５８８に規定されたプロトコルを用いて、第１電子機器１０の内部時刻を基準時刻に同期させる。座標系算出部３０は、例えば、ＩＥＥＥ１５８８に規定されたプロトコルを用いて、座標系算出部３０の内部時刻を基準時刻に同期させることができる。
【００６５】
［位置情報計測部１２の基準座標系のキャリブレーションについて］
位置情報計測部１２は、前述したように第１計測座標系を有することから、第１計測座標系を基準座標系とすることができる。しかしながら、位置情報計測部１２が１台しかない場合、何らかの物体が、位置情報計測部１２と第２電子機器２０ｎのマーカー２３（ｎ）との間に介在した場合、位置情報計測部１２から第２電子機器２０ｎのマーカー２３（ｎ）の座標値を計測できなくなる可能性がある。
第１実施形態においては、このようなオクルージョンによる死角を生じないように、複数の位置情報計測部１２（ｉ）（１≦ｉ≦Ｍ）を備えることが好ましい。ただし、Ｍ＝１の場合を含む。
その場合、位置情報計測部１２を備える第１電子機器１０を複数備えるようにシステムを構成してもよい。また、第１電子機器１０が複数の位置情報計測部１２（ｉ）を備えるように構成してもよい。
第１実施形態においては、複数の第１計測座標系（ｉ）が存在する場合に対応できるように、第１計測座標系（ｉ）とは独立した基準座標系を備えるように構成する。
【００６６】
この場合、各位置情報計測部１２（ｉ）を基準座標系における座標値が既知となる場所に設置することで、各第１計測座標系（ｉ）における座標値を基準座標系における座標値に変換するための基準座標変換行列（以下、「ｉＴＷ」ともいう）を予め算出することができる。
【００６７】
これに対して、基準座標系を設定するための基準座標系マーカーを用いて、各第１計測座標系（ｉ）における座標値を基準座標系における座標値に変換するための基準座標変換行列ｉＴＷを予め算出するように、キャリブレーションを行うようにしてもよい。
以下、キャリブレーションについて説明する。
【００６８】
基準座標系を設定するために、基準座標系マーカーを現実空間中に設置する。なお、基準座標系マーカー１０Ａは、キャリブレーションを行った後に、現実空間から撤去される。
【００６９】
基準座標系マーカー１０Ａは、例えば第２電子機器２０の備えるマーカー２３と同様に、各頂点を同一の半径の球体とした３辺を有する三角形の形状をした構造体とすることができる。球体の半径、球体の表面の色、及び三角形の辺の長さは、適宜設定することができる。ここでは、計算式を簡単にするために、例えば、基準座標系マーカー１０Ａは、３辺の長さがそれぞれ３ｊｃｍ、４ｊｃｍ、５ｊｃｍ（ｊは、１以上の任意の整数）となる直角三角形と仮定するが、三角形を任意の形状とした場合も同様に、キャリブレーションをすることができる。
【００７０】
基準座標系マーカー１０Ａにおける直角三角形の直角に対応する頂点をＷ０、辺の長さが４ｊｃｍとなる辺のＷ０以外の頂点をＷ１、辺の長さが３ｊｃｍとなる辺のＷ０以外の頂点をＷ２とする。Ｗ０、Ｗ１、Ｗ２は、それぞれ、各球体の中心と一致する。
そして、頂点Ｗ０を始点として頂点Ｗ１を終点とするベクトル（以下「ベクトルＷ０Ｗ１」という）と、頂点Ｗ０を始点として頂点Ｗ２を終点とするベクトル（以下「ベクトルＷ０Ｗ２」という）との外積ベクトルを頂点Ｗ０を始点として表現した場合の終点をＷ３とする。
【００７１】
このようにして、一例として、基準座標系マーカー１０Ａにより、Ｗ０を原点、ベクトルＷｏＷ１をＸ軸、ベクトルＷ０Ｗ２をＹ軸、ベクトルＷ０Ｗ３をＺ軸とする、基準座標系を規定することができる。
具体的には、各点Ｗ０、Ｗ１、Ｗ２、Ｗ３の基準座標系における座標値は、次のとおりである。
Ｗ０＝（０，０，０）
Ｗ１＝（４ｊ，０，０）
Ｗ２＝（０，３ｊ，０）
Ｗ３＝（０，０，１２ｊ２）
なお、上記座標値は、一例であって、マーカー２３が直角三角形でない任意の三角形とした場合、一例として、例えば、Ｗ０を原点として、ベクトルＷ０Ｗ１をＸ軸とし、三角形Ｗ０Ｗ１Ｗ２をＸＹ平面とすることで、基準座標系を規定することができる。その場合の座標値は、当該三角形に基づいて算出される。
基準座標系は、基準座標系マーカー１０Ａにより任意に設定することができる。
【００７２】
位置情報計測部１２（ｉ）のキャリブレーションのために、マーカー位置情報算出部３１は、位置情報計測部１２（ｉ）により計測された計測情報に基づいて、基準座標系マーカー１０Ａを検出し、検出した基準座標系マーカー１０Ａの構成要素（基準座標系マーカー１０Ａの頂点Ｗ０、Ｗ１、Ｗ２）の第１計測座標系（ｉ）における座標値を算出する。
算出された基準座標系マーカー１０Ａの頂点Ｗ０、Ｗ１、Ｗ２の第１計測座標系（ｉ）における座標値をそれぞれ
Ｗ０＝（ａ０，ｂ０，ｃ０）
Ｗ１＝（ａ１，ｂ１，ｃ１）
Ｗ２＝（ａ２，ｂ２，ｃ２）
Ｗ３＝（ａ３，ｂ３，ｃ３）とする。
ここで、ベクトルＷ０Ｗ３は、ベクトルＷ０Ｗ１とベクトルＷ０Ｗ２との外積であることから、Ｗ３の座標値は、次のように、Ｗ０、Ｗ１、Ｗ２の座標値で表すことができる。
ａ３＝（ｂ１−ｂ０）・（ｃ２−ｃ０）−（ｃ１−ｃ０）（ｂ２−ｂ０）＋ａ０
ｂ３＝（ｃ１−ｃ０）・（ａ２−ａ０）−（ａ１−ａ０）・（ｃ２−ｃ０）＋ｂ０
ｃ３＝（ａ１−ａ０）・（ｂ２−ｂ０）−（ｂ１−ｂ０）・（ａ２−ａ０）＋ｃ０
【００７３】
第１電子機器１０の位置情報計測部１２（ｉ）の備える第１計測座標系（ｉ）における座標値を基準座標系における座標値に変換するための基準座標変換行列をｉＴＷとすると、以下の式１が成立する。
【数１】
0006293110.tif 000002
【００７４】
キャリブレーション処理部３２は、基準座標系マーカー１０Ａの頂点Ｗ０、Ｗ１、Ｗ２の位置を表す第１計測座標系における座標値に基づいて、式２によって、位置情報計測部１２（ｉ）における第１計測座標系（ｉ）における座標値を基準座標系における座標値に変換するための基準座標変換行列ｉＴＷを算出することができる。
【数２】
0006293110.tif 000003
第１電子機器１０は、キャリブレーション処理部３２により算出した基準座標変換行列ｉＴＷを、例えば３次元データ処理部４０の記憶部（図示せず）に記憶するようにしてもよい。
なお、三角形を任意の形状とした場合に、各点Ｗ０、Ｗ１、Ｗ２、Ｗ３の基準座標系における座標値を以下の値、
Ｗ０＝（０，０，０）
Ｗ１＝（ｘ１，０，０）
Ｗ２＝（ｘ２，ｙ２，０）
Ｗ３＝（ｘ３，ｙ３，ｚ３）
とすると、
式２´によって、位置情報計測部１２における第１計測座標系における座標値を基準座標系における座標値に変換するための基準座標変換行列ｉＴＷを算出することができる。
【数３】
0006293110.tif 000004
【００７５】
［基準座標系マーカー１０Ａの検出方法について］
マーカー位置情報算出部３１による基準座標系マーカー１０Ａの検出方法の一例について説明する。なお、ここで説明する検出方法は、マーカー位置情報算出部３１が第２電子機器２０ｎの備えるマーカー２３（ｎ）を検出する際にも適用される。
【００７６】
ここでは、説明を簡単にするために、基準座標系マーカー１０Ａに対応する球体の表面の色をＲＥＤ、球体の半径を５ｃｍ、直角三角形の各辺の長さを３０ｃｍ、４０ｃｍ、５０ｃｍと仮定する。
また、キャリブレーションに際しては、計測対象範囲内には、基準座標系マーカー１０Ａ以外に（特に、半径５ｃｍ、色がＲＥＤとなる）球体を設置しないようにすることが好ましい。
【００７７】
マーカー位置情報算出部３１は、位置情報計測部１２（ｉ）により計測された計測対象範囲内において色がＲＥＤの部分を抽出して、当該抽出したＲＥＤの部分の深度画像データ、すなわち当該抽出したＲＥＤの部分表面の各点のＲＧＢの色情報である（Ｒ）と、各点（画素単位）の位置を表す第１計測座標系（ｉ）における３次元座標（ｘ、ｙ、ｚ）と、からなる点データの集合を計測画像データとして取得し、記憶部（図示せず）に記憶する。
【００７８】
マーカー位置情報算出部３１は、位置情報計測部１２（ｉ）により計測された計測画像データ及び記憶部に格納された色がＲＥＤである計測画像データの位置情報に基づいて、基準座標系マーカー１０Ａを構成する球体を示す円を抽出する。
具体的には、マーカー位置情報算出部３１は、計測画像データから例えば公知のラプラシアン変換、Ｓｏｂｅｌ法又はＣａｎｎｙ法等により、物体の輪郭部分の抽出を行い、エッジ画像を生成する。この際、マーカー位置情報算出部３１は、記憶部に格納された色がＲＥＤである計測画像データとその周辺を含む画像データを対象として、エッジ画像を生成するようにしてもよい。マーカー位置情報算出部３１は、エッジ画像から公知のＨｏｕｇｈ変換により円を抽出する。マーカー位置情報算出部３１は、抽出された円の座標が記憶部に格納された色がＲＥＤである計測画像データに含まれるものを抽出する。マーカー位置情報算出部３１は、さらに、抽出された円の半径を計算することにより、半径が５ｃｍとなる円のみを抽出するように構成される。こうすることで、マーカー位置情報算出部３１は、基準座標系マーカー１０Ａの直角三角形の各頂点に対応する各球体に対応する円を抽出する。
【００７９】
［基準座標系マーカー１０Ａの頂点Ｗ０、Ｗ１、Ｗ２の第１計測座標系における座標値］
マーカー位置情報算出部３１は、抽出された円の深度データの中から最も小さな値となる点、すなわち最も視平面、すなわち位置情報計測部１２（ｉ）（深度カメラ）の平面に近い点の座標を（ａ，ｂ，ｃ）とした場合に、Ｚ軸方向に垂直に（球体の半径となる）５ｃｍプラスした点の座標（ａ，ｂ，ｃ＋５）を球の中心として求めることができる。
また、円上の直径を形成する任意の２点の中心となる座標値を求めることで、各球体の中心の座標値を算出するようにしてもよい。
例えば、抽出された円の直径を形成する円上の任意の２点Ａ１及びＡ２の座標値をそれぞれ、（ｘ１、ｙ１、ｚ１）及び（ｘ２、ｙ２、ｚ２）とした場合、球体の中心の座標値（ａ´、ｂ´、ｃ´）は、式３のように算出することができる。
ａ´＝（ｘ１＋ｘ２）／２
ｂ´＝（ｙ１＋ｙ２）／２
ｃ´＝（ｚ１＋ｚ２）／２                                    （式３）
【００８０】
マーカー位置情報算出部３１は、このようにして基準座標系マーカー１０Ａの各頂点の第１計測座標系（ｉ）における座標値を算出する。その後、マーカー位置情報算出部３１は、算出した基準座標系マーカー１０Ａの各頂点の第１計測座標系（ｉ）における座標値に基づいて、各三角形の辺の長さを算出し、算出した辺の長さに基づいて、各頂点が基準座標系マーカー１０Ａの頂点Ｗ０、Ｗ１、Ｗ２のいずれに対応するか、を判定することができる。
【００８１】
このように、マーカー位置情報算出部３１は、基準座標系マーカー１０Ａの頂点Ｗ０、Ｗ１、Ｗ２の位置を表す第１計測座標系（ｉ）における座標値を算出することができる。
【００８２】
上記の例においては、基準座標系マーカー１０Ａに対応する球体の表面の色をＲＥＤ、球体の半径を５ｃｍと仮定したが、球体の色をＲＥＤ以外の色（例えば、色Ａ）、球体の半径を５ｃｍ以外の長さ（例えば、半径Ｒｃｍ）とすることができる。
この場合、マーカー位置情報算出部３１は、計測対象範囲内の色が色Ａの部分を抽出して、抽出した色Ａの部分表面の各点のＲＧＢの色情報と、各点（画素単位）の位置を表す３次元座標（ｘ、ｙ、ｚ）と、からなる点データの集合を計測画像データとして抽出すればよい。
そうした上で、マーカー位置情報算出部３１は、エッジ画像を生成し、生成したエッジ画像からＨｏｕｇｈ変換により円を抽出し、円の半径を計算することにより、色Ａで半径Ｒの球体の中心を求めるように構成することができる。
以上、マーカー位置情報算出部３１による、基準座標系マーカー１０Ａの検出方法について説明した。
【００８３】
［マーカー２３の検出］
基準座標系マーカー１０Ａの検出方法は、後述する第２電子機器２０ｎの備えるマーカー２３（ｎ）の検出にも適用される。
なお、マーカー２３（ｎ）を構成する三角形の辺の長さを各第２電子機器２０ｎ毎に異なる長さにすることで、位置情報計測部１２の計測対象範囲に複数の第２電子機器２０ｎが計測される場合において、辺の長さにより、当該マーカー２３（ｎ）を有する第２電子機器２０を識別するように構成することができる。
同様に、マーカー２３（ｎ）を構成する球体の半径を各第２電子機器２０ｎ毎に異なる長さに予め設定しておくことで、位置情報計測部１２（ｉ）の計測対象範囲に複数の第２電子機器２０ｎが計測される場合において、球体の半径により、当該マーカー２３（ｎ）を有する第２電子機器２０ｎを識別するように構成することができる。
同様に、マーカー２３（ｎ）を構成する球体の表面の色を各第２電子機器２０ｎ毎に異なる色に予め設定しておくことで、位置情報計測部１２（ｉ）の計測対象範囲に複数の第２電子機器２０ｎが計測される場合において、球体の色により、当該マーカー２３（ｎ）を有する第２電子機器２０ｎを識別するように構成することができる。
【００８４】
［マーカー座標系と第２座標系との間の座標変換］
第１実施形態は、座標系算出部３０の備えるマーカー位置情報算出部３１が、位置情報計測部１２（ｉ）により計測された計測情報に基づいて、第２電子機器２０ｎの備えるマーカー２３（ｎ）を検出し、検出したマーカー２３（ｎ）の構成要素（マーカー２３の頂点等）の第１計測座標系（ｉ）における座標値（マーカー位置情報）を算出し、マーカー位置情報に基づいて、第２電子機器２０ｎの備える第２座標系（ｎ）における座標値を基準座標系における座標値に変換するための座標変換行列を算出するというものである。
このためには、第２電子機器２０ｎの３次元計測部２２ｎの備える第２座標系（ｎ）における座標値を第２電子機器２０ｎの備えるマーカー２３（ｎ）によるマーカー座標系による座標値に変換するためのマーカー座標変換行列（以下、「ｎＴＭ」という）を予め算出しておく必要がある。
【００８５】
この場合、マーカー２３（ｎ）における三角形の頂点（Ｃ´０、Ｃ´１、Ｃ´２）を第２座標系（ｎ）における座標値が既知となる位置に取り付けることで、第２座標系（ｎ）における座標値をマーカー座標系（ｎ）における座標値に変換するためのマーカー座標変換行列ｎＴＭを予め算出することができる。
【００８６】
これに対して、マーカー２３（ｎ）における三角形の頂点（Ｃ´０、Ｃ´１、Ｃ´２）の第２座標系（ｎ）における座標値が必ずしも既知でない場合又は取り付け時に誤差が生じる場合、第２電子機器２０ｎを予め固定した位置においた状態で、キャリブレーションを行うことで、マーカー座標変換行列ｎＴＭを算出することができる。
以下、マーカー２３（ｎ）によるキャリブレーションについて説明する。
【００８７】
図３に、第１計測座標系（ｉ）における座標値を基準座標系における座標値に変換するための基準座標変換行列ｉＴＷ、３次元計測部２２ｎの備える第２座標系（ｎ）における座標値を基準座標系における座標値に変換するための座標変換行列（以下、「ｎＴＷ」という）、第２座標系（ｎ）における座標値をマーカー座標系（ｎ）における座標値に変換するためのマーカー座標変換行列ｎＴＭ、マーカー座標系（ｎ）における座標値を第１計測座標系（ｉ）における座標値に変換するための第１計測座標変換行列（以下、「ＭｎＴｉ」という）間の関連を示す。
【００８８】
マーカー座標変換行列ｎＴＭを算出するために、第２電子機器２０ｎを予め固定した位置においた状態で、３次元計測部２２ｎの備える第２座標系（ｎ）における座標値を基準座標系における座標値に変換するための座標変換行列ｎＴＷを算出する。
これと同時に、第２電子機器２０ｎを予め固定した位置においた状態で、第２電子機器２０ｎの備えるマーカー２３（ｎ）によるマーカー座標系（ｎ）による座標値を第１電子機器１０の位置情報計測部１２（ｉ）の備える第１計測座標系（ｉ）における座標値に変換するための第１計測座標変換行列ＭｎＴｉを算出する。
【００８９】
［第２座標系を基準座標系に変換するための座標変換行列ｎＴＷ］
まず、第２電子機器２０ｎを予め固定した位置においた状態で、３次元計測部２２ｎの備える第２座標系（ｎ）における座標値を基準座標系における座標値に変換するための座標変換行列ｎＴＷを算出するための構成について説明する。
ここで、第２電子機器２０ｎの備える３次元計測部２２ｎによる、基準座標系マーカー１０Ａの検出については、第１電子機器１０の位置情報計測部１２（ｉ）による基準座標系マーカー１０Ａの検出方法と同様の方法が適用できる。
すなわち、マーカー位置情報算出部３１は、第２電子機器２０ｎの３次元計測部２２ｎにより計測された計測情報に基づいて、基準座標系マーカー１０Ａを検出し、その頂点Ｗ０、Ｗ１、Ｗ２の位置を表す第２座標系（ｎ）における座標値を算出する。そうすることで、キャリブレーション処理部３２は、３次元計測部２２ｎの備える第２座標系（ｎ）における座標値を基準座標系における座標値に変換するための座標変換行列ｎＴＷを算出することができる。
【００９０】
なお、上記の例では、座標系算出部３０の備えるマーカー位置情報算出部３１が、３次元計測部２２ｎにより計測した計測画像データに基づいて、基準座標系マーカー１０Ａを検出し、検出した基準座標系マーカー１０Ａの構成要素（基準座標系マーカー１０Ａの頂点Ｗ０、Ｗ１、Ｗ２）の第２座標系（ｎ）における座標値を算出した。
しかしながら、第２電子機器２０ｎにそれぞれマーカー位置情報算出部（図示せず）を設けて、第２電子機器２０ｎのマーカー位置情報算出部により、基準座標系マーカー１０Ａの構成要素（基準座標系マーカー１０Ａの頂点Ｗ０、Ｗ１、Ｗ２）の第２座標系（ｎ）における座標値を算出するように構成してもよい。
【００９１】
［マーカー座標系を第１計測座標系に変換するための第１計測座標変換行列Ｔ０］
次に、第２電子機器２０ｎを予め固定した位置においた状態で、マーカー座標系（ｎ）における座標値を第１計測座標系（ｉ）における座標値に変換するための第１計測座標変換行列ＭｎＴｉを算出するための構成について説明する。
マーカー２３（ｎ）によるマーカー座標系（ｎ）とは、前述したとおり、マーカー２３（ｎ）における三角形の頂点をそれぞれ、原点、Ｘ軸上の点、Ｙ軸上の点とする座標系を意味する。
なお、説明を簡単にするため、マーカー２３（ｎ）における三角形を３辺の比が３：４：５となる直角三角形とする。
具体的には、各点Ｃ´０、Ｃ´１、Ｃ´２、Ｃ´３のマーカー座標系（ｎ）における座標値は、次のとおりである。
Ｃ´０＝（０，０，０）
Ｃ´１＝（４ｉ，０，０）
Ｃ´２＝（０，３ｉ，０）
Ｃ´３＝（０，０，１２ｉ２）
【００９２】
前述したように、マーカー位置情報算出部３１は、位置情報計測部１２（ｉ）による計測情報に基づいて、マーカー２３（ｎ）を検出し、検出したマーカー２３（ｎ）の構成要素（マーカー２３（ｎ）の頂点Ｃ´０、Ｃ´１、Ｃ´２）の第１計測座標系（ｉ）における座標値を算出する。
【００９３】
このようにして、算出されたマーカー２３（ｎ）の頂点Ｃ´０、Ｃ´１、Ｃ´２の第１計測座標系（ｉ）における座標値を
Ｃ´０＝（ｐ０，ｑ０，ｒ０）
Ｃ´１＝（ｐ１，ｑ１，ｒ１）
Ｃ´２＝（ｐ２，ｑ２，ｒ２）
Ｃ´３＝（ｐ３，ｑ３，ｒ３）とする。
ここで、ベクトルＣ´０Ｃ´３は、ベクトルＣ´０Ｃ´１とベクトルＣ´０Ｃ´２との外積であることから、Ｃ´３の座標値は、次のように、Ｃ´０、Ｃ´１、Ｃ´２の座標値で表すことができる。
ｐ３＝（ｑ１−ｑ０）・（ｒ２−ｒ０）−（ｒ１−ｒ０）（ｑ２−ｑ０）＋ｐ０
ｑ３＝（ｒ１−ｒ０）・（ｐ２−ｐ０）−（ｐ１−ｐ０）・（ｒ２−ｒ０）＋ｑ０
ｒ３＝（ｐ１−ｐ０）・（ｑ２−ｑ０）−（ｑ１−ｑ０）・（ｐ２−ｐ０）＋ｒ０
【００９４】
こうすると、以下の式４が成立する。
【数４】
0006293110.tif 000005
【００９５】
このように、キャリブレーション処理部３２は、マーカー２３（ｎ）の頂点Ｃ´０、Ｃ´１、Ｃ´２の位置を表す第１計測座標系（ｉ）における座標値に基づいて、式５によって、マーカー座標系（ｎ）を第１計測座標系（ｉ）に変換するための第１計測座標変換行列ＭｎＴｉを算出することができる。
【数５】
0006293110.tif 000006
なお、上記式は、マーカー２３（ｎ）における三角形を３辺の比が３：４：５となる直角三角形とした場合の一例であって、マーカー２３（ｎ）を任意の三角形とした場合、各点Ｃ´０、Ｃ´１、Ｃ´２、Ｃ´３のマーカー座標系（ｎ）における座標値は、当該三角形に基づいて、算出される。
【００９６】
［第２座標系をマーカー座標系に変換するためのマーカー座標変換行列ｎＴＭ］
第２電子機器２０ｎの３次元計測部２２ｎの備える第２座標系（ｎ）における座標値を第２電子機器２０ｎの備えるマーカー２３（ｎ）によるマーカー座標系（ｎ）による座標値に変換するためのマーカー座標変換行列ｎＴＭとする。
図３を参照すると、第１計測座標系（ｉ）における座標値を基準座標系における座標値に変換するための基準座標変換行列ｉＴＷ、基準座標系における座標値を第２座標系（ｎ）における座標値に変換するための第２座標変換行列ｎＴＷ−１、第２座標系（ｎ）における座標値をマーカー座標系（ｎ）における座標値に変換するためのマーカー座標変換行列ｎＴＭ、マーカー座標系（ｎ）における座標値を第１計測座標系（ｉ）における座標値に変換するための第１計測座標変換行列ＭｎＴｉの順番に掛け合わせることにより、式６に示すように、第１計測座標系における座標値を第１計測座標系における座標値に変換する単位行列Ｅとなる。
ＭｎＴｉ・ｎＴＭ・ｎＴＷ−１・ｉＴＷ＝Ｅ        （式６）
以上から、キャリブレーション処理部３２は、式７によって、第２電子機器２０ｎの３次元計測部２２ｎの備える第２座標系（ｎ）における座標値を第２電子機器２０ｎの備えるマーカー２３（ｎ）によるマーカー座標系（ｎ）における座標値に変換するためのマーカー座標変換行列ｎＴＭを算出することができる。
ｎＴＭ＝ＭｎＴｉ−１・ｉＴＷ−１・ｎＴＷ        （式７）
【００９７】
このように、キャリブレーション処理部３２は、位置情報計測部１２（ｉ）における第１計測座標系（ｉ）における座標値を基準座標系における座標値に変換するための基準座標変換行列ｉＴＷ、及び第２電子機器２０ｎの３次元計測部２２ｎの備える第２座標系（ｎ）における座標値を第２電子機器２０ｎの備えるマーカー座標系（ｎ）における座標値に変換するためのマーカー座標変換行列ｎＴＭを算出する。
ここで、基準座標変換行列ｉＴＷ、及びマーカー座標変換行列ｎＴＭは、それぞれ固定の変換行列である。
以上のように、各位置情報計測部１２（ｉ）（１≦ｉ≦Ｍ）、及び各第２電子機器２０ｎ（１≦ｎ≦Ｎ）のキャリブレーションを行い、第１計測座標系（ｉ）における座標値を基準座標系における座標値に変換するための基準座標変換行列ｉＴＷ、及び第２座標系（ｎ）における座標値をマーカー座標系（ｎ）における座標値に変換するためのマーカー座標変換行列ｎＴＭを算出した後、第１電子機器１０（キャリブレーション処理部３２）は、例えば３次元データ処理部４０の記憶部に記憶するようにしてもよい。その後、基準座標系マーカー１０Ａを現実空間から撤去することができる。
【００９８】
［基準座標変換行列の補正について］
なお、位置情報計測部１２（ｉ）（１≦ｉ≦Ｍ）の備える第１計測座標系（ｉ）における座標値を基準座標系における座標値に変換するための基準座標変換行列ｉＴＷとすると、現実空間上に例えば、基準座標系マーカー１０Ａを適当に置いた場合、当該基準座標系マーカー１０Ａの各頂点Ｗ０、Ｗ１、Ｗ２及びベクトルＷ０Ｗ１とベクトルＷ０Ｗ２との外積ベクトルをＷ０Ｗ３とした場合の、各点Ｗ０、Ｗ１、Ｗ２、Ｗ３の位置情報計測部１２（ｉ）（１≦ｉ≦Ｍ）の第１計測座標系（ｉ）における座標値をそれぞれ、Ｗ０（ｉ）、Ｗ１（ｉ）、Ｗ２（ｉ）、Ｗ３（ｉ）とする。
そうすると、これらの座標値を第１計測座標系（ｉ）における座標値を基準座標系における座標値に変換するための基準座標変換行列ｉＴＷで変換して得られる基準座標系における座標値ｉＴＷ（Ｗ０（ｉ））、ｉＴＷ（Ｗ１（ｉ））、ｉＴＷ（Ｗ２（ｉ））、ｉＴＷ（Ｗ３（ｉ））は、全て同一の座標値となるはずである。
しかしながら、仮に、位置情報計測部１２（ｉ）（１≦ｉ≦Ｍ）による測定誤差により、座標値がずれる場合、例えば、ｉが２以上の第１計測座標系（ｉ）における座標値を基準座標系における座標値を基準座標系に変換した座標値｛ｉＴＷ（Ｗ０（ｉ））、ｉＴＷ（Ｗ１（ｉ））、ｉＴＷ（Ｗ２（ｉ））、ｉＴＷ（Ｗ３（ｉ））｝（ｉ≧２）の値が、１番目の第１計測座標系（１）における座標値を基準座標系における座標値を基準座標系に変換した座標値｛１ＴＷ（Ｗ０（ｉ））、１ＴＷ（Ｗ１（ｉ））、１ＴＷ（Ｗ２（ｉ））、１ＴＷ（Ｗ３（ｉ））｝に一致するように、第１計測座標系（ｉ）（ｉ≧２）における座標値を基準座標系における座標値に変換するための基準座標変換行列ｉＴＷ（ｉ≧２）を補正するようにしてもよい。
そうすることで、補正後の第１計測座標系（ｉ）（ｉ≧２）における座標値を基準座標系における座標値を基準座標系に変換した座標値｛ｉＴＷ（Ｗ０（ｉ））、ｉＴＷ（Ｗ１（ｉ））、ｉＴＷ（Ｗ２（ｉ））、ｉＴＷ（Ｗ３（ｉ））｝（ｉ≧２）の値は、１番目の第１計測座標系（１）における座標値を基準座標系における座標値を基準座標系に変換した座標値｛１ＴＷ（Ｗ０（ｉ））、１ＴＷ（Ｗ１（ｉ））、１ＴＷ（Ｗ２（ｉ））、１ＴＷ（Ｗ３（ｉ））｝に一致するようにできる。
【００９９】
以上、第１電子機器１０の備える第１機能部、第２電子機器２０の備える第２機能部、及び座標系算出部３０の備える座標系算出機能部について説明した。
【０１００】
このように、第１電子機器１０の位置情報計測部１２（ｉ）及び座標系算出部３０は、第２電子機器１０ｎの備えるマーカー２３（ｎ）のマーカー位置情報（ｎ）を正確に算出することができる。そうすることで、３次元計測部２２ｎの視線方向及び視線周りの回転角が決まり、例えば従来のように例えば加速度センサやジャイロセンサ等により視線方向及び視線周りの回転角を算出するよりも、正確な位置情報を得ることができる。
次に、３次元データ処理部４０について説明する。
【０１０１】
＜３次元データ処理部４０について＞
図１Ｅに示すように、３次元データ処理部４０は、点群データ座標値変換部４１と、合成点群データ作成部４２と、を含む。また、合成点群データ作成部４２は、さらに、第１合成点群データ作成部４２１と、点群データマッチング部４２２と、を備える。
【０１０２】
［点群データ座標値変換部４１］
ｎを第２電子機器２０（３次元計測部２２）のインデックス（識別番号）、計測時刻をｔとしたとき、点群データ座標値変換部４１は、点群データ｛Ｑｎ（ｔ）｝（１≦ｎ≦Ｎ）（ｔ：タイムスタンプ）を、同一座標系（基準座標系）の座標に変換するための座標変換行列｛Ｆｎ（ｔ）｝（１≦ｉ≦Ｍ）（１≦ｎ≦Ｎ）（ｔ：タイムスタンプ）を算出し、各点群データＱｎ（ｔ）の各点を表す第２座標系における座標値から基準座標系における座標値への座標変換を行う。
【０１０３】
以下に、点群データ座標値変換部４１が、時刻ｔにおいて位置情報計測部１２（ｉ）により計測される第２電子機器２０ｎのマーカー位置情報ｉＬｎ（ｔ）（ここで、ｉは位置情報計測部１２の識別番号（１≦ｉ≦Ｍ）、ｎは第２電子機器２０の識別番号、ｔはタイムスタンプ）に基づいて、各点群データＱｎ（ｔ）の各点を示す第２座標系における座標値を基準座標系における座標値に変換するための変換行列Ｆｎ（ｔ）を算出する処理について、説明する。
【０１０４】
前述したように、キャリブレーションにより、位置情報計測部１２（ｉ）の備える第１計測座標系（ｉ）（１≦ｉ≦Ｍ）における座標値を基準座標系における座標値に変換するための基準座標変換行列ｉＴＷが、３次元データ処理部４０の記憶部に記憶されている。
同様に、キャリブレーションにより、第２電子機器２０ｎ（１≦ｎ≦Ｎ）の３次元計測部２２ｎの備える第２座標系（ｎ）における座標値を第２電子機器２０ｎの備えるマーカー２３（ｎ）によるマーカー座標系（ｎ）における座標値に変換するためのマーカー座標変換行列ｎＴＭが、３次元データ処理部４０の記憶部に記憶されている。
【０１０５】
まず、点群データ座標値変換部４１は、マーカー位置情報の集合｛ｉＬｎ（ｔ）｝（１≦ｎ≦Ｎ，１≦ｉ≦Ｍ，ｔ：タイムスタンプ）と３次元計測部２２ｎにより計測された点群データの集合｛Ｑｎ（ｔ）｝（１≦ｎ≦Ｎ）（ｔ：タイムスタンプ）のタイムスタンプとを突合せすることで、点群データの計測位置を特定することができる。
なお、点群データＱｎ（ｔ）を測定したときの測定時刻とマーカー位置情報ｉＬｎ（ｔ）を測定したときの測定時刻とが必ずしも一致しない場合には、例えば点群データＱｎ（ｔ）を測定したときの測定時刻の近傍の時刻におけるマーカー２３（ｎ）のマーカー位置情報を用いて、線形補完することができる。
【０１０６】
点群データ座標値変換部４１は、時刻ｔ及び第２電子機器２０の識別番号ｎに基づいて、各点群データＱｎ（ｔ）とマーカー位置情報ｉＬｎ（ｔ）とを対応づけることができる。
点群データ座標値変換部４１は、Ｑｎ（ｔ）に対して、ｉＬｎ（ｔ）を対応付ける。
なお、位置情報計測部１２（ｉ）を複数備える場合、Ｑｎ（ｔ）に対して、異なる位置情報計測部１２（ｉ）及び位置情報計測部１２（ｊ）（ｉ≠ｊ）によりそれぞれ算出された位置情報ｉＬｎ（ｔ）及びｊＬｎ（ｔ）（ｉ≠ｊ）が対応付けられる可能性がある。
そこで、点群データ座標値変換部４１は、Ｑｎ（ｔ）に対して、複数のマーカー位置情報から、所定の条件を満たす１つのマーカー構成要素位置情報ｉＬｎ（ｔ）を選択する。
【０１０７】
例えば、選択するための条件としては、（ｉ、ｔ）に対して、集合｛ｉＬｎ（ｔ）｝の個数が最大となる（ｉ）を選択するようにしてもよい。すなわち、位置情報計測部１２（ｉ）により、マーカー位置情報（ｎ）を取得することができた第２電子機器２０の個数が最大となるｉを選択する。
こうすることで、後述する座標値変換の効率をあげることが可能となる。
【０１０８】
そこで、Ｑｎ（ｔ）に対してｉＬｎ（ｔ）が選択されたとする。
【０１０９】
前述したように、点群データ座標値変換部４１は、マーカー位置情報ｉＬｎ（ｔ）に基づいて、マーカー座標系（ｎ）における座標値を第１計測座標系（ｉ）における座標値に変換するための第１計測座標変換行列ＭｎＴｉ（ｔ）を算出することができる（式５参照）。
【０１１０】
他方、予めキャリブレーションによって算出された第２座標系（ｎ）における座標値をマーカー座標系（ｎ）における座標値に変換するためのマーカー座標変換行列ｎＴＭ、第１計測座標系（ｉ）における座標値を基準座標系における座標値に変換するための基準座標変換行列をｉＴＷとすると、点群データ座標値変換部４１は、タイムスタンプｔにおいて、点群データの点Ｑｎ（ｔ）の第２座標系（ｎ）での座標が（ｘ，ｙ，ｚ）で表される点に対して、マーカー座標変換行列ｎＴＭ、第１計測座標変換行列ＭｎＴｉ（ｔ）、及び基準座標変換行列ｉＴＷによる座標変換を行うことで、基準座標系での座標（ｘｇ，ｙｇ，ｚｇ）を算出する（式８参照）。
すなわち、Ｆｎ（ｔ）＝ｉＴＷ・ＭｎＴｉ（ｔ）・ｎＴＭとなる。
【数６】
0006293110.tif 000007
このようにして、点群データ座標値変換部４１は、点群データＱｎ（ｔ）の各点を示す第２座標系における座標値を基準座標系における座標値にリアルタイムに変換することができる。
【０１１１】
マーカー位置情報算出部３１は、マーカー位置情報ｉＬｎ（ｔ）をリアルタイムに算出して、第２電子機器２０ｎ毎に、ネットワークを介して、３次元データ処理部４０（点群データ座標値変換部４１）に、例えば通信等のインターフェースを介して引き渡すように構成することができる。
なお、マーカー位置情報算出部３１は、マーカー位置情報ｉＬｎ（ｔ）をリアルタイムに算出して、第２電子機器２０ｎ毎に、ネットワークを介して３次元データ保存部６０の第２電子機器位置情報記憶部６２にリアルタイムに記憶し、第２電子機器位置情報記憶部６２を介して３次元データ処理部４０（点群データ座標値変換部４１）に引き渡すようにしてもよい。
このようにして、本発明の点群データ取得システム１００は、点群データＱｎ（ｔ）の各点を基準座標系での座標値で構成される点群データ｛Ｑ´ｎ（ｔ）｝にリアルタイムに変換することができる。
【０１１２】
［合成点群データ作成部４２］
合成点群データ作成部４２は、第１合成点群データ作成部４２１と、点群データマッチング部４２２と、を備える。
【０１１３】
［第１合成点群データ作成部４２１］
第１合成点群データ作成部４２１は、点群データ座標値変換部４１により、時刻ｔにおける各点群データＱｎ（ｔ）の各点の第２座標系での座標値を、基準座標系での座標値に変換した各点群データ｛Ｑ´ｎ（ｔ）｝  （１≦ｎ≦Ｎ）間の基準座標系での座標値が同じ点については、同一点として、重ね合わせることで、合成点群データＱ´（ｔ）を作成することができる。
【０１１４】
［点群データマッチング部４２２］
しかし、各点群データ｛Ｑ´ｎ（ｔ）｝（１≦ｎ≦Ｎ）間の基準座標系での座標値が同じ点であっても、測定誤差等により、位置ずれが発生する可能性がある。
このため、点群データマッチング部４２２は、点群データ座標値変換部４１により、基準座標系での座標に変換した各点群データ｛Ｑ´ｎ（ｔ）｝の各点の座標値に基づいて、それぞれの点群データ｛Ｑ´ｎ（ｔ）｝（１≦ｎ≦Ｎ）から抽出した３次元形状に係る特徴パラメータ（例えば、曲面の曲率（ガウス曲率））に基づいて重なる部分の位置合わせ（マッチング）を、行うことで、位置ずれを補正することができる。
そうすることで、同一の基準座標系での座標で表される複数の点群データ｛Ｑ´ｎ（ｔ）｝（１≦ｎ≦Ｎ）を用いて、より正確な欠損部分の少ない３次元立体となる合成点群データ（以下、「Ｑ（ｔ）」と表す）を取得することができる。
【０１１５】
なお、隣り合う計測時刻（例えば、ｔｉ、ｔｉ＋１）に、同じ３次元計測部２２ｎにより計測された点群データＱｎ（ｔｉ）及びＱｎ（ｔｉ＋１）について、位置ずれがある場合にも、必要に応じて、重なる部分の位置合わせ（マッチング）を行うことで、位置ずれを補正することができる。
なお、位置ずれの補正については、バックグラウンド処理として、例えば、１秒１回程度のサンプリングとしてもよい。
【０１１６】
点群データの位置合わせ（マッチング）は、例えば周知の技術であるＩＣＰ法と呼ばれる位置合わせ手法を用いることもできる。ＩＣＰ法は、位置合わせを行うそれぞれの点群データについて、最近傍頂点間距離の二乗和を反復計算により最小化（収束）させることで位置合わせを行うものである。
しかしながら、ＩＣＰ法は総当たり方式であるため、データ量が多くなると、膨大な時間を要していた。
また、位置合わせのための判定指標として、色情報を用いる方法もある。しかしながら、色は、同一個所であっても、見る角度や光の当たり具合によって違って見えることもあり、したがって、異なる位置から測定された点群データＱｎ（ｔ）とＱｍ（ｔ）（ｎ≠ｍ）とを、色を判定指標とするマッチングは不安定であった。
【０１１７】
点群データマッチング部４２２では、判定指標として、各点の形状面での特性を示す「曲面の曲率」を用いるようにしてもよい。点群データ間において、頂点や尖った点等の特徴点ではなく、「ある程度の広がりを持った面」を考え、曲面の曲率を用いることで、それらの「凹凸具合の一致度」を見ることで、判定指標とする。
曲面の曲率は、座標系に依存しない物理的数学的特性であって、同じ曲面の同一箇所の凸凹度合いは、どの位置（すなわち、どの３次元計測部２２）から見ても同じになることから、マッチングの判定指標として使用することができる。「曲面の曲率」を用いることで、どの３次元計測部２２から測定されても、「曲面」の曲率は同一の値をとることが想定されることから、位置合わせ時のエラー発生を減少させることができる。また、マーカー設置の手間が無く、色情報と違って光が不要なため、光源（太陽等）の具合に左右されることがない。
なお、曲面の曲率としては、例えばガウス曲率、平均曲率等があるが、第１実施形態では、ガウス曲率を適用することができる。
【０１１８】
また、マッチング処理の前処理もまた重要である。
具体的には、点群データマッチング部４２２は、まず複数の点群データを基準座標系での座標値を突き合わせ、その積集合を得ることで、マッチング処理対象となる点群データの範囲を限定する。このことにより、マッチング処理に要する計算量を大幅に減らすことができ、システム全体の速度向上に寄与することができる。また、積集合が空集合の場合は、そもそもマッチング処理をしなくてよい。
次に、点群データマッチング部４２２は、部分的に重なる、点群データＡ及び点群データＢの重なる部分において、「両曲面」上の「ある程度の広がりをもった面」上の各点の曲率を求め、それらを比較して、例えば曲率の差の合計が最小となるような最適な変換を、これら一連の処理を繰り返しながら求める。
【０１１９】
ここで、点群データマッチング部４２２は、曲率の計算を、例えば次のように行うことができる。まず、曲面の曲率を求める１つの点群データＱｎ（ｔ）の点（ｘ０、ｙ０、ｚ０）を原点とし、その点を含む、例えば基準座標系の平面Ｘ及びＹに平行となる、平面Ｘ’Ｙ’を含む局所座標系（ｘ´、ｙ´、ｚ´）を設定する。次に、その点（ｘ０、ｙ０、ｚ０）の近傍点における「曲面」の局所座標を求める。この「曲面」の局所座標値は、点（ｘ０、ｙ０、ｚ０）の近傍にある点群データの点（ｘ、ｙ、ｚ）により算出できる。そうすることで、「曲面」上のｚ値と（ｘ、ｙ）との関係を得ることができる。
例えば近傍点は、ｘ０の隣り合う値を例えばｘ−２、ｘ−１、ｘ０、ｘ＋１、ｘ＋２等として表すと、
（ｘ＋１、ｙ０）、（ｘ＋２、ｙ０）、・・・
（ｘ−１、ｙ０）、（ｘ−２、ｙ０）、・・・
（ｘ０、ｙ＋１）、（ｘ０、ｙ＋２）、・・・
（ｘ０、ｙ−１）、（ｘ０、ｙ−２）、・・・
等として求めることができる。
【０１２０】
次に２次曲面パッチｐ（ｕ、ｖ）＝（ｕ、ｖ、ａｕ２＋ｂｖ２＋ｃｕｖ＋ｄｕ＋ｅｖ＋ｆ）を計測対象物表面に最小二乗法によって局所的に当てはめて、２次元曲面パッチの偏導関数から、曲面の曲率を計算する。
最後に、対象物の曲面が滑らかであることを制約条件として用いて曲率計算の精度を反復計算によって向上させて、算出する。
しかし、この曲率を用いる方法は、正確ではあるが計算量が多くなる。あるいは、物体の場所によっては（エッジ部分等）、うまく適用できない場合がある。その場合は、従来のＩＣＰ法、色情報を用いる方法等を組み合わせることで、対応することができる。
【０１２１】
［直線のＨｏｕｇｈ変換を用いた、マッチング方法］
以上、曲率を用いる方法について、説明したが、直線のＨｏｕｇｈ変換を用いた点群データのマッチング方法について説明する。
【０１２２】
２つの深度画像をＡ，Ｂとし、基準座標系におけるその点群データをそれぞれＡｗ，Ｂｗとする。
点群データマッチング部４２２は、深度画像Ａ及び深度画像Ｂのそれぞれにおいて、エッジ画像を生成し、生成した各エッジ画像からＨｏｇｈ変換を用いて、それぞれ直線を抽出し、それぞれから抽出した直線を第２座標系上の直線（３次元）に変換する。
点群データマッチング部４２２は、深度画像Ａ及び深度画像Ｂそれぞれから抽出した直線を基準座標系上の直線に変換する。
次に、深度画像Ａ及び深度画像Ｂから抽出された基準座標系上の直線（３次元直線）をそれぞれ、４次元上の点として表現する。
より具体的には、点群データマッチング部４２２は、深度画像Ａから得られた３次元直線をＡｗｉ（０≦ｉ≦ｎ１）、深度画像Ｂから得られた３次元直線をＢｗｊ（０≦ｊ≦ｎ２）とした場合、３次元直線Ａｗｉ（０≦ｉ≦ｎ１）及び３次元直線Ｂｗｊ（０≦ｊ≦ｎ２）をそれぞれ４次元空間上の点として以下のように表わすことができる。
Ａｗｉ＝（αＡｉ，βＡｉ，ｌＡｉ，ｍＡｉ）  （０≦ｉ≦ｎ１）
Ｂｗｊ＝（αＢｊ，βＢｊ，ｌＢｊ，ｍＢｊ）  （０≦ｊ≦ｎ２）
【０１２３】
一般に、３次元直線Ｌを４次元上の点として表す変換としては、例えば次の例が挙げられるが、この例に限定されるものではない。
直線Ｌの方向ベクトル（原点Ｏを始点とする方向ベクトル）の終点Ｐから、ＸＺ平面上に垂直に引いた直線のＸＺ平面との交点Ｑとすると、ＯＱとＸ軸との角度α、辺ＱＯと辺ＰＯとの角度をβとする。こうすることで、直線Ｌの方向ベクトルが決定される。
次に、方向ベクトルの法平面上と直線Ｌとの交点の座標を（ｍ、ｌ）とすると、直線Ｌは、４次元上の点（α，β，ｌ，ｍ）で表すことができる。
【０１２４】
そうすることで、点群データマッチング部４２２は、全てのｉ（０≦ｉ≦ｎ１）とｊ（０≦ｊ≦ｎ２）について、４次元上の点ＡｗｉとＢｗｊとの４次元上の距離Ｄｉｊを算出し、その距離の値が０に近い値となる組み合わせを４個を抽出する。
こうすることで、深度画像Ａ及び深度画像Ｂにおいてマッチングする直線を４個抽出することができる。
このようにして、点群データマッチング部４２２は、４個のペア（Ａｗｉ，Ｂｗｊ）を抽出した後、各点Ａｗｉ及びＢｗｊを３次元空間上の直線に変換して、直線Ｂｗｊを直線Ａｗｉに変換するための行列を算出することで、点群データＢｗを点群データＡｗとにマッチングするための変換行列を得ることができる。
【０１２５】
このように、点群データマッチング部４２２は、これら一連の作業を、各点群データ毎に行い、それぞれの一致度合いを見ながら、複数の点群データを合成していくことで、マッチングを行うことができる。そして、その結果は、後述する３次元データ表示制御部５０に、例えば通信等のインターフェースを介して、引き渡すように構成することができる。なお、後述する３次元データ表示制御部５０に引き渡すに際して、例えば、３次元データ保存部６０の記憶部（図示せず）に格納し、当該記憶部を介して引き渡すようにしてもよい。
なお、上述した点群データマッチング部４２２のマッチング処理として、ＩＣＰマッチング処理、曲率を用いたマッチング処理、及び直線を用いたマッチング処理について説明したが、これらのマッチング処理に限定されるものではない。これらのマッチング処理を組み合わせて使用してもよい。また、別のマッチングアルゴリズムを使用してもよい。
【０１２６】
このようにして、点群データマッチング部４２２は、点群データ座標値変換部４１により、同時刻に計測された複数の点群データ（例えば、｛Ｑｎ（ｔ）｝（１≦ｎ≦Ｎ））を、点群データ座標値変換部４１により変換された基準座標系（第１座標系）における座標値に基づいて重なる部分の位置合わせを行った後に、必要に応じて、測定誤差等により、各点群データ｛Ｑｎ（ｔ）｝（１≦ｎ≦Ｎ）間の基準座標系での座標値が同じであっても、位置ずれの発生した場合の補正を行うことで、より正確な、欠損部分の少ない１つの合成点群データを取得することができる。
こうすることで、例えば、形状の変化する計測対象物の多視点同時刻の点群データの合成を行うことができる。
【０１２７】
なお、形状の変化しない静止した計測対象物の場合には、第２電子機器２０ｎを移動させないで、すなわち、基準座標系における座標値が同じ値となる、同一位置において計測された複数の点群データを、点群データ座標値変換部４１により変換された基準座標系（第１座標系）における座標値に基づいて、重なる部分の位置合わせを行うことで、１つの合成点群データを作成することができる。
この場合は、タイムスタンプを省略することができる。
【０１２８】
また、形状の変化しない静止した計測対象物の場合には、１つの第２電子機器２０を当該計測対象物の周りを移動しながら計測することで得られる複数の点群データを、点群データ座標値変換部４１により変換された基準座標系における座標値に基づいて、重なる部分の位置合わせを行うことで、１つの合成点群データを作成することができる。
【０１２９】
＜３次元データ表示制御部５０＞
次に、３次元データ表示制御部５０について説明する。
図１Ｆに示すように、３次元データ表示制御部５０は、視点操作部５１と、データ表示方式選択部５２と、データ再生・保存指示部５３と、３次元データ表示処理部５４と、を備える。
【０１３０】
点群データは、そのままの形式でレンダリングすることは可能であるが、そのデータ形式のままでは、各種の３次元処理には適さないことが多い。
このため、点群データを面形式に変換し、例えばポリゴン、不整三角網のメッシュ、及びＣＡＤモデルとして扱うことができる。
また、点群データは、各点データが三次元座標値を持つため、自由視点からのデータ閲覧及び空間参照を実現できる。具体的には、点群データ自体が有する３次元座標値を用いて、点群データを３次元空間から、例えばパノラマ空間への投影が可能となる。
【０１３１】
［視点操作部５１］
視点操作部５１は、利用者の視点方向、視点座標等の入力を利用者から受け取ることができる。このことにより、表示画像のズーム等のカメラワークと同様の操作が可能となる。
【０１３２】
［データ表示方式選択部５２］
データ表示方式選択部５２は、表示画像の表示方式に関する指示を利用者から受け取ることができる。例えば、３次元空間上の点群データをどのような投影面に投影するか、を選択することができる。例えばパノラマモデルは、球面、円柱、立方体等の投影面を選択できる。
【０１３３】
［データ再生・保存指示部５３］
データ再生・保存指示部５３は、点群データから生成された３次元モデルデータの再生又は保存に関する指示を利用者から受け取ることができる。
【０１３４】
［３次元データ表示処理部５４］
３次元データ表示処理部５４は、合成点群データ作成部４２により作成された合成点群データを、視点操作部５１により入力された利用者の仮想視点に基づいて、マッピング加工を行い、その結果の表示画像をディスプレイ画面に３次元表示する。
また、３次元データ表示処理部５４は、動的な３次元モデルを用いて、視点操作部５１から指定された視点パラメータで３次元空間を眺めた動画像を描画することができる。
【０１３５】
次に、３次元データ保存部６０について説明する。
図１Ｇに示すように、３次元データ保存部６０は、源点群データ記憶部６１と、第２電子機器位置情報記憶部６２と、を備える。
【０１３６】
［源点群データ記憶部６１］
源点群データ記憶部６１は、各３次元計測部２２により計測された点群データに、当該計測時刻となるタイムスタンプを関連付けた点群データ情報を、第２電子機器２０ｎ毎に記憶する。
前述したように、例えば、各３次元計測部２２ｎにより計測対象物表面の各点を画素単位で、予め設定した時間間隔（例えば、１秒間に３０コマ）で取得する場合、源点群データ記憶部６１は、３次元計測部２２ｎ毎に、例えば、１／３０秒単位の計測時刻（ｔ）毎に生成される点群データ｛Ｑｎ（ｔ）｝を記憶する（１≦ｎ≦Ｎ）（ｔ：タイムスタンプ）。
【０１３７】
［第２電子機器位置情報記憶部６２］
第２電子機器位置情報記憶部６２は、第２電子機器２０の備えるマーカー２３のマーカー位置情報にタイムスタンプを関連付けたマーカー位置情報を、第２電子機器２０ｎ毎、位置情報計測部１２（ｉ）毎に記憶する。
前述したように、第２電子機器位置情報記憶部６２は、すべての第２電子機器２０ｎの、時刻ｔに測定される第２電子機器位置情報の集合｛ｉＬｎ（ｔ）｝（１≦ｎ≦Ｎ）（１≦ｉ≦Ｍ）（ｔ：タイムスタンプ）を記憶する。
【０１３８】
以上、これまで説明した本発明の実施形態における、一連の処理の流れを簡単に図示したものが図４である。
【０１３９】
図４を参照すると、ＳＴ１において、基準座標系マーカー１０Ａを所定の場所に設定される。各位置情報計測部１２（ｉ）（１≦ｉ≦Ｍ）が所定の場所に設定される。同様に、第２電子機器２０ｎ（１≦ｎ≦Ｎ）が所定の場所に設定される。
【０１４０】
ＳＴ２において、キャリブレーションを行うことで、第１計測座標系（ｉ）における座標値を基準座標系における座標値に変換するための基準座標変換行列ｉＴＷ、及び第２座標系（ｎ）における座標値をマーカー座標系（ｎ）における座標値に変換するためのマーカー座標変換行列ｎＴＭを算出し、３次元データ処理部４０の記憶部に記憶する。
【０１４１】
ＳＴ３において、基準座標系マーカー１０Ａを現実空間から撤去する。
【０１４２】
ＳＴ４において、例えば、計測作業者が、移動しながら、各第２電子機器２０ｎ（３次元計測部２２ｎ）により、計測対象物の計測を開始する。
【０１４３】
ＳＴ５において、測定時刻ｔ毎に、各第２電子機器２０ｎ（３次元計測部２２ｎ）は点群データＱｎ（ｔ）を取得する。
【０１４４】
ＳＴ６において、測定時刻ｔ毎に、マーカー位置情報算出部３１は、各位置情報計測部１２（ｉ）（１≦ｉ≦Ｍ）の計測した計測情報に基づいて、各第２電子機器２０ｎのマーカー位置情報ｉＬｎ（ｔ）を算出する。
【０１４５】
ＳＴ７において、３次元データ処理部４０により、３次元計測部２２ｎの取得した点群データＱｎ（ｔ）と、各位置情報計測部１２（ｉ）の算出した各第２電子機器２０ｎのマーカー位置情報ｉＬｎ（ｔ）との対応付けがなされる。
【０１４６】
ＳＴ８において、３次元データ処理部４０（点群データ座標値変換部４１）により、点群データＱｎ（ｔ）の各点の第２座標系（ｎ）における座標値を基準座標系における座標値に変換する。
【０１４７】
ＳＴ９において、３次元データ処理部４０（合成点群データ作成部４２）により、同時刻ｔにおける各点群データＱｎ（ｔ）のマッチング・結合（合成）がなされる。
【０１４８】
ＳＴ１０において、３次元データ表示制御部５０により、ユーザの仮想的な視点の指定にしたがい、当該仮想的な視点から見た３次元動画が表示される。
【０１４９】
以上のように、第１実施形態によれば、各３次元計測部２２を自由かつダイナミックに移動させることが可能となり、計測対象物の変形に伴って生じるオクルージョンによる死角が発生する場合においても、３次元計測装置の計測位置を即座に変更して対応することで。欠損部分の少ない３次元画像を取得することができる。
【０１５０】
また、本発明は、このような一連の処理により、合成点群データを視点操作部５１によって指定入力された視点の位置から計測対象物を見たものとしてクリッピングし、欠損部分の少ない三次元物体の画像を、任意の視点からリアルタイムで見ることができる。
【０１５１】
例えば、医療分野において、資料映像として残すために、手術中の患部の状況を画像撮影することがある。その際、カメラが手術の邪魔になるからと、他の位置に移動させたとしても、見る側は、実際のカメラの位置とは無関係に、常に同じ視点からの映像を見ることが可能になる。それゆえ、この方式で資料化された映像は、きわめて有用性の高いものとなる。
【０１５２】
また、手術者の体で隠れる視点からの映像として表示することができる。また、腹腔鏡手術においては、内臓の状況を患者の皮膚の外側からの視点で見ることもできる。
【０１５３】
また、本発明を撮影カメラに応用することで、スポーツ実況中継等において、その都度最適な位置にカメラを移動することで、死角の少ないで欠損部分を少なくし、ありえない視点から選手の動きを、各々の画像装置で表示することができるようになること等、様々なシーンにおいて、「撮影位置とは異なった視点からの画像表示対象物体を見ることができる」ができるようになる。
【０１５４】
また、合成点群データをサーバ上におき、これをネットワーク経由で例えば、スマートフォン、タブレット端末、ノートパソコン等を用いて、複数のユーザが同一の物体を任意の視点から閲覧することができる。
【０１５５】
＜第２実施形態＞
＜電磁波伝播時間を用いた第２電子機器２０の位置情報の算出＞
第１実施形態では、第１電子機器１０の備える位置情報計測部１２により、第２電子機器２０の備えるマーカー２３の位置情報を計測することで、第２電子機器２０の備える３次元計測部２２の備える第２座標系から第１電子機器の備える第１計測座標系への座標変換行列を算出した。
これに対して、第２実施形態では、電磁波（電波）が第１電子機器１０の備える電波発信部１２Ｐｊ（１≦ｊ≦４）から発信されてから、第２電子機器２０の備える電波受信部２３Ｃｋ（１≦ｋ≦３）により受信されるまでの１２個の電波伝播時間測定値（Ｔｊｋ）を測定し、各電波発信部１２Ｐｊと各電波受信部２３Ｃｋ間の直線距離Ｒｊｋ（１≦ｊ≦４，１≦ｋ≦３）を算出し、直線距離Ｒｊｋ（１≦ｊ≦４，１≦ｋ≦３）を用いて三角測量の原理に基づいて、第２電子機器２０の備える電波受信部２３Ｃｋ（１≦ｋ≦３）の位置情報を算出する。
なお、以下、簡単のため、電波発信部１２Ｐｊ（１≦ｊ≦４）及び電波受信部２３Ｃｋ（１≦ｋ≦３）を、それぞれ単にＰｊ（１≦ｊ≦４）とＣｋ（１≦ｋ≦３）と略して説明することもある。
【０１５６】
図５Ａ、図５Ｂ、図５Ｃは、それぞれ、第２実施形態における点群データ取得システムの第１電子機器１０、第２電子機器２０、及び座標系算出部３０の機能ブロックの一例を示す図である。
【０１５７】
図５Ａに示すように、第１電子機器１０は、位置情報計測部１２に換えて、複数（例えば４つ）の電波発信部１２Ｐｊ（１≦ｊ≦４）及び距離計測部１６を備える。
【０１５８】
図５Ｂに示すように、第２電子機器２０は、マーカー２３に換えて、複数（例えば３つ）の電波受信部２３Ｃｋ（１≦ｋ≦３）を備える。複数（例えば３つ）の電波受信部２３Ｃｋ（１≦ｋ≦３）を「位置情報ユニット」ともいう。
ここで、第１電子機器１０は基準座標系（第１座標系）を備えており、電波発信部Ｐｊ（１≦ｊ≦４）は、例えば、基準座標系（第１座標系）のＸ軸、Ｙ軸、Ｚ軸上、線形独立となる位置に設定されるものとする。同様に、電波受信部２３Ｃｋ（１≦ｋ≦３）は、第２電子機器２０側の３次元計測部２２の備える第２座標系の固有座標系のＸ軸、Ｙ軸、Ｚ軸上、線形独立となる位置に設定されるものとする。
【０１５９】
図５Ｃに示すように、座標系算出部３０は、第２電子機器位置情報算出部３１Ａを備える。
電波発信部１２Ｐｊ（１≦ｊ≦４）、距離計測部１６、及び第２電子機器位置情報算出部３１Ａを総称して、「位置情報算出部」という。
【０１６０】
第１電子機器１０は、距離計測部１６を備える。距離計測部１６は、各電波発信部１２Ｐｊ（１≦ｊ≦４）から送信された電波を各電波受信部２３Ｃｋ（１≦ｋ≦３）が受信することにより、送信、受信の時間間隔（伝播時間）を算出し、その距離を計測する。
【０１６１】
より、具体的には、第２電子機器２０の備える電波受信部２３Ｃｋ（１≦ｋ≦３）は、第１電子機器１０の電波発信部１２Ｐｊ（１≦ｊ≦４）から信号を受信すると、例えば、電波受信部２３Ｃｋ（１≦ｋ≦３）の各識別信号（電波受信部２３Ｃｋ（１≦ｋ≦３）の各ＩＤ）を付加して反射させる。
こうすることで、第１電子機器１０が、距離計測部１６を備えることで、電波発信部１２Ｐｊ（１≦ｊ≦４）から送信された電波を各電波受信部２３Ｃｋ（１≦ｋ≦３）が反射した反射波を受信することにより、送信、受信の時間間隔（伝播時間）が算出され、その距離を計測することができる。
【０１６２】
距離計測部１６により、各電波発信部１２Ｐｊ（１≦ｊ≦４）と各電波受信部２３Ｃｋ（１≦ｋ≦３）間の直線距離Ｒｊｋ（１≦ｊ≦４，１≦ｋ≦３）を算出した後は、第２電子機器位置情報算出部３１Ａは、三角測量の原理に基づいて、第２電子機器２０の備える電波受信部２３Ｃｋ（１≦ｋ≦３）の第１電子機器の備える第１座標系（基準座標系）における位置情報を算出することができる。
【０１６３】
なお、タイムスタンプについては、各電波発信部１２Ｐｊ（１≦ｊ≦４）の送信時刻となるタイムスタンプを付して、第２電子機器２０ｎ毎に、第２電子機器位置情報Ｌｎ（ｔ）（ここで、ｎは第２電子機器２０の識別番号、ｔはタイムスタンプ）を３次元データ保存部６０の第２電子機器位置情報記憶部６２にリアルタイムに記憶することができる。
【０１６４】
［変形例］
前述のように、電波受信部２３Ｃｋ（１≦ｋ≦３）が受信した電波を反射させる方式に換えて、各電波発信部１２Ｐｊ（１≦ｊ≦４）から送信された電波を各電波受信部２３Ｃｋ（１≦ｋ≦３）が受信することにより、送信、受信の時間間隔（伝播時間）を算出し、その距離を計測するようにしてもよい。
この場合、第１電子機器１０の備える電波発信部１２Ｐｊ（１≦ｊ≦４）は、電波発信部１２Ｐｊ（１≦ｊ≦４）の各識別信号（電波発信部１２Ｐｊ（１≦ｊ≦４）の各ＩＤ）を付加して送信させるようにしてもよい。
このようにして、距離計測部１６は、第１電子機器１０側の基準座標系（例えば、Ｘ軸、Ｙ軸、Ｚ軸上の線形独立となる４点Ｐｊ（１≦ｊ≦４）から、第２電子機器２０側の第２座標系（例えば、Ｘ軸、Ｙ軸、Ｚ軸上の線形独立となる３点Ｃｋ（１≦ｋ≦３）へのそれぞれの距離Ｒｊｋ（１≦ｊ≦４，１≦ｋ≦３）を計測することができる。
【０１６５】
距離計測部１６により、各電波発信部１２Ｐｊ（１≦ｊ≦４）と各電波受信部２３Ｃｋ（１≦ｋ≦３）間の直線距離Ｒｊｋ（１≦ｊ≦４，１≦ｋ≦３）を算出した後は、三角測量の原理に基づいて、第２電子機器位置情報算出部３１Ａは、第２電子機器２０の有する各３次元計測部２２の第１座標系における位置情報を算出することができる。
【０１６６】
次に、各第２電子機器２０の第２電子機器位置情報の算出方法について説明する。
まず、図６を参照して、距離計測部１６により、電波が第１電子機器１０の各電波発信部１２Ｐｊ（１≦ｊ≦４）から発信されてから、各第２電子機器２０の各電波受信部２３Ｃｋ（１≦ｋ≦３）により受信されるまでの時間を計測することにより２点間の距離を算出する処理について簡単に説明する。
【０１６７】
ここでは、計算式を簡単にするため、各電波受信部２３Ｃｋ（１≦ｋ≦３）を３次元計測部２２の備える第２座標軸に一致して、第２座標軸の原点Ｃ０から１０ｃｍの距離に設置していると仮定する。
なお、図６を参照すると、例えば、レンズの焦点（あるいは、カメラの重心）を原点Ｃ０とし、光軸方向をＺ’軸、水平方向をＸ’軸、垂直方向をＹ’軸とすると、Ｃ１、Ｃ２、Ｃ３はそれぞれ、Ｚ’軸、Ｘ’軸、Ｙ’軸に一致して、原点Ｃ０から一定の距離例えば１０ｃｍ離れた点に設置する。
すなわち、Ｃ０，Ｃ１，Ｃ２，Ｃ３の第２座標系における座標値は、以下のとおりとする（単位ｃｍ）。
Ｃ０の座標値＝（０，０，０）、
Ｃ１の座標値＝（０，０，１０）、
Ｃ２の座標値＝（１０，０，０）
Ｃ３の座標値＝（０，１０，０）
【０１６８】
同様に、計算式を簡単にするため、４つの電波発信部１２Ｐｊ（１≦ｊ≦４）のうち、Ｐ１を基準座標系（第１座標系）の原点に設置し、Ｐ２、Ｐ３、Ｐ４をそれぞれ基準座標軸（第１座標軸）に一致して、基準座標軸（第１座標軸）の原点Ｐ１からそれぞれ１０ｃｍの距離に設置していると仮定する。
図６を参照すると、例えば、４つの電波発信部１２Ｐｊ（１≦ｊ≦４）は、それぞれ基準座標系（第１座標系）の原点、Ｘ軸、Ｙ軸、及びＺ軸に一致して設置され、各Ｐ１〜Ｐ３は、原点Ｐ１から一定の距離例えば１０ｃｍ離れた点に設置する。
すなわち、Ｐ１、Ｐ２、Ｐ３、Ｐ４の基準座標系（第１座標系）における座標値は、以下のとおりとする（単位ｃｍ）。
Ｐ１の座標値＝（０，０，０）、
Ｐ２の座標値＝（１０，０，０）、
Ｐ３の座標値＝（０，１０，０）、
Ｐ４の座標値＝（０，０，１０）
【０１６９】
図６には、各電波発信部１２Ｐｊ（１≦ｊ≦４）、及び各電波受信部２３Ｃｋ（１≦ｋ≦３）の位置関係を示している。
図６に示すように、点Ｃ０、Ｃ１、Ｃ２、及びＣ３の基準座標系（第１座標系）における座標値をそれぞれ（ｘ０，ｙ０，ｚ０）、（ｘ１，ｙ１，ｚ１）、（ｘ２，ｙ２，ｚ２）、及び（ｘ３，ｙ３，ｚ３）と表す。
【０１７０】
図６に示すように、各電波発信部１２Ｐｊ（１≦ｊ≦４）と各電波受信部２３Ｃｋ（１≦ｋ≦３）間の直線距離は全部で１２本あり、これらをＲｊｋ（１≦ｊ≦４、１≦ｋ≦３）で表す。これらの距離Ｒｊｋは、その二点間の電波の伝播時間Ｔｊｋ（１≦ｊ≦４、１≦ｋ≦３）及び電波の速度から算出することができる。
【０１７１】
第１電子機器１０は、各電波発信部１２Ｐｊ（１≦ｊ≦４）から順次、予め設定された時間間隔（例えば、６ミリ秒間隔）で、電波を発信する。各電波受信部２３Ｃｋ（１≦ｋ≦３）は、各電波発信部１２Ｐｊ（１≦ｊ≦４）の発信した電波を反射させることで、各電波発信部１２Ｐｊ（１≦ｊ≦４）は、反射された電波を受信する。
このようにして、距離計測部１６は、各電波発信部１２Ｐｊ（１≦ｊ≦４）から発信されて、各電波受信部２３Ｃｋ（１≦ｋ≦３）により受信されるまでのそれぞれの電波伝播時間測定値Ｔｊｋ（１≦ｊ≦４、１≦ｋ≦３）に基づいて、各電波発信部１２Ｐｊ（１≦ｊ≦４）と各電波受信部２３Ｃｋ（１≦ｋ≦３）間の直線距離Ｒｊｋ（１≦ｊ≦４、１≦ｋ≦３）を算出することができる。
この際、前述したように、各電波発信部１２Ｐｊ（１≦ｊ≦４）の電波送信時刻となるタイムスタンプを直線距離｛Ｒｊｋ｝（１≦ｊ≦4，１≦ｋ≦３）に付すことができ、当該タイムスタンプに基づいて第２電子機器位置計測結果間で関連付けが可能となる。
【０１７２】
＜３次元計測部２２の位置情報の算出について＞
次に、第２電子機器位置情報算出部３１Ａが、各電波発信部１２Ｐｊ（１≦ｊ≦４）と各電波受信部２３Ｃｋ（１≦ｋ≦３）間の各直線距離Ｒｊｋ（１≦ｊ≦４、１≦ｋ≦３）に基づいて、各点Ｃ０、Ｃ１、Ｃ２、Ｃ３の基準座標系（第１座標系）における座標値を算出する処理について、前述の例に基づいて説明する。
ここで、各点Ｃ０、Ｃ１、Ｃ２、Ｃ３の基準座標系（第１座標系）における座標値を（ｘ０，ｙ０，ｚ０）、（ｘ１，ｙ１，ｚ１）、（ｘ２，ｙ２，ｚ２）、及び（ｘ３，ｙ３，ｚ３）とする。
【０１７３】
そうすると、三平方の定理に基づき、Ｐ１とＣ１間の距離、Ｐ２とＣ１間の距離、Ｐ３とＣ１間の距離、Ｐ４とＣ１間の距離は、以下の式９〜式１２が成り立つ。
ｘ１2＋ｙ１2＋ｚ１2  ＝  Ｒ１１２                          （式９）
（ｘ１−１０）2＋ｙ１2＋ｚ１2  ＝  Ｒ２１２                （式１０）
ｘ１2＋（ｙ１−１０）2＋ｚ１2  ＝  Ｒ３１２                （式１１）
ｘ１2＋ｙ１2＋（ｚ１−１０）2  ＝  Ｒ４１２                （式１２）
【０１７４】
第２電子機器位置情報算出部３１Ａは、式９〜式１２に基づいて、Ｃ１の基準座標系（第１座標系）における座標値（ｘ１、ｙ１、ｚ１）を式１３により算出する。
ｘ１＝（Ｒ２１２  −  Ｒ１１２  ＋１００）／２０          
ｙ１＝（Ｒ３１２  −  Ｒ１１２  ＋１００）／２０          
ｚ１＝（Ｒ４１２  −  Ｒ１１２  ＋１００）／２０          
・・・（式１３）
【０１７５】
同様に、第２電子機器位置情報算出部３１Ａは、Ｃ２の基準座標系（第１座標系）における座標値（ｘ２、ｙ２、ｚ２）を、Ｐ１とＣ２間の距離、Ｐ２とＣ２間の距離、Ｐ３とＣ２間の距離、Ｐ４とＣ２間の距離に基づいて、Ｃ２の基準座標系（第１座標系）における座標値（ｘ２、ｙ２、ｚ２）を式１４により算出する。
ｘ２＝（Ｒ２２２  −  Ｒ１２２  ＋１００）／２０          
ｙ２＝（Ｒ３２２  −  Ｒ１２２  ＋１００）／２０          
ｚ２＝（Ｒ４２２  −  Ｒ１２２  ＋１００）／２０          
・・・（式１４）
【０１７６】
同様に、第２電子機器位置情報算出部３１Ａは、Ｃ３の基準座標系（第１座標系）における座標値（ｘ３、ｙ３、ｚ３）を、Ｐ１とＣ３間の距離、Ｐ２とＣ３間の距離、Ｐ３とＣ３間の距離、Ｐ４とＣ３間の距離に基づいて、Ｃ３の基準座標系（第１座標系）における座標値（ｘ３、ｙ３、ｚ３）を式１５により算出する。
ｘ３＝（Ｒ２３２  −  Ｒ１３２  ＋１００）／２０          
ｙ３＝（Ｒ３３２  −  Ｒ１３２  ＋１００）／２０          
ｚ３＝（Ｒ４３２  −  Ｒ１３２  ＋１００）／２０          
・・・（式１５）
【０１７７】
次に、第２電子機器位置情報算出部３１Ａは、Ｃ０の基準座標系（第１座標系）における座標値（ｘ０，ｙ０，ｚ０）を（式１６）〜（式１８）に基づいて算出する。
（ｘ０−ｘ１）2＋（ｙ０−ｙ１）2＋（ｚ０−ｚ１）2  ＝  １０２    （式１６）
（ｘ０−ｘ２）2＋（ｙ０−ｙ２）2＋（ｚ０−ｚ２）2  ＝  １０２    （式１７）
（ｘ０−ｘ３）2＋（ｙ０−ｙ３）2＋（ｚ０−ｚ３）2  ＝  １０２    （式１８）
【０１７８】
こうすることで、第２電子機器位置情報算出部３１Ａは、第２電子機器２０の３次元計測部２２の固有座標系（第２座標系）における原点Ｃ０（０，０，０）、Ｚ´軸上の点Ｃ１（０，０，１０）、Ｘ´軸上の点Ｃ２（１０，０，０）、及びＹ´軸上の点Ｃ３（０，１０，０）の基準座標系（第１座標系）における座標値を各々算出することができる。
なお、Ｃ１、Ｃ２、Ｃ３の位置を計測するためには、第１電子機器１０側の超音波発信部１２は３つあればよいが、本実施形態ではＰ１、Ｐ２、Ｐ３、Ｐ４の４つとした。そうすることで、計算を簡単にするとともに、測定精度を向上させることができる。
【０１７９】
以上のようにして、第２電子機器位置情報算出部３１Ａは、三角測量の原理に基づいて、第２電子機器２０の有する各３次元計測部２２の基準座標系（第１座標系）における位置情報を算出することができる。
【０１８０】
なお、４つの電波発信部１２Ｐｊ（１≦ｊ≦４）の位置は、線形独立である任意の位置とすることができる。また、３つの電波受信部２３Ｃｋ（１≦ｋ≦３）の位置についても、第２座標系の原点Ｃ０とした場合に、始点をＣ０、終点をそれぞれＣ１、Ｃ２、Ｃ３とする３つのベクトルが、線形独立であれば、Ｃ１、Ｃ２、Ｃ３の位置は、任意の位置とすることができる。
【０１８１】
この場合においても、第２電子機器位置情報算出部３１Ａは、上記と同様に、三角測量の原理により、各Ｃｊ（１≦ｊ≦３）の基準座標系（第１座標系）における座標値を、Ｐ１とＣｊ（１≦ｊ≦３）間の距離、Ｐ２とＣｊ（１≦ｊ≦３）間の距離、Ｐ３とＣｊ（１≦ｊ≦３）間の距離、Ｐ４とＣｊ（１≦ｊ≦３）間の距離と、Ｐ１、Ｐ２、Ｐ３、Ｐ４の基準座標系（第１座標系）における座標値に基づいて、求めることができる。
【０１８２】
そして、第２電子機器位置情報算出部３１Ａは、上記と同様に、第２座標系の原点Ｃ０の基準座標系（第１座標系）における座標値を、各Ｃｊ（１≦ｊ≦３）の基準座標系（第１座標系）における座標値と各Ｃｊ（１≦ｊ≦３）の第２座標系における座標値とに基づいて、算出することができる。
【０１８３】
以上のように、第２電子機器位置情報算出部３１Ａは、第２電子機器２０（３次元計測部２２）の位置を表す各電波受信部２３Ｃｋ（１≦ｋ≦３）の基準座標系（第１座標系）における座標値（第２電子機器位置計測結果）を算出することができる。
【０１８４】
なお、前述したように、第２電子機器位置情報算出部３１Ａにより算出された各電波受信部２３Ｃｋ（１≦ｋ≦３）の基準座標系（第１座標系）における座標値は、各電波発信部１２Ｐｊ（１≦ｊ≦４）の送信時刻となるタイムスタンプを付して、第２電子機器２０ｎ毎に、第２電子機器位置情報Ｌｎ（ｔ）（ここで、ｎは第２電子機器２０の識別番号、ｔはタイムスタンプ）を３次元データ保存部６０の第２電子機器位置情報記憶部６２にリアルタイムに記憶することができる。
距離計測部１６及び第２電子機器位置情報算出部３１Ａの一連の処理はリアルタイムに処理することができる。
【０１８５】
［第２実施形態におけるキャリブレーションの適用］
第２実施形態においては、第１電子機器１０の電波発信部１２Ｐｊ（１≦ｊ≦４）に基づく座標系を基準座標系（第１座標系）とし、第２電子機器２０ｎの３次元計測部２２の備える第２座標系とした場合に、第２電子機器２０ｎの電波受信部２３Ｃｋ（１≦ｋ≦３）を第２座標系における座標値が既知となる位置に取り付けることを前提とした。
これに対して、第２電子機器２０ｎの電波受信部２３Ｃｋ（１≦ｋ≦３）の第２座標系における座標値が必ずしも既知でない場合又は取り付け時に誤差が生じる場合、第２電子機器２０ｎを予め固定した位置においた状態で、第１実施形態と同様に、キャリブレーションを行うことができる。
【０１８６】
説明を簡単にするため、第２電子機器２０ｎの電波受信部２３Ｃｋ（１≦ｋ≦３）に基づく座標系を電波受信部座標系という。
ここで、説明を簡単にするため、電波受信部２３Ｃｋ（１≦ｋ≦３）の電波受信部座標系における座標値を、例えば、次のようにする（ｉは、所定の自然数）。
Ｃ１＝（０，０，０）
Ｃ２＝（４ｉ，０，０）
Ｃ３＝（０，３ｉ，０）
そして、ベクトルＣ１Ｃ２とベクトルＣ１Ｃ３との外積を頂点Ｃ１を始点として表現した場合における終点をＣ４とする。
Ｃ４＝（０，０，１２ｉ２）
また、第２電子機器２０ｎの第２座標系における座標値を電波受信部座標系における座標値に変換するための変換行列を電波受信部座標変換行列を、「ｎＴＭ」とする。
【０１８７】
第１実施形態の場合と同様に、基準座標系マーカー１０Ａを用意することで、第１座標系における座標値を基準座標系における座標値に変換するための基準座標変換行列１Ｔｗ及び第２電子機器２０ｎの第２座標系における座標値を電波受信部座標系における座標値に変換するための変換行列（電波受信部座標変換行列ｎＴＭ）を算出することができる。
例えば、基準座標系マーカー１０Ａの球体の中心に電波受信部（Ｗ０、Ｗ１、Ｗ２）を設置する。そして、基準系マーカーの各頂点Ｗ０、Ｗ１、Ｗ２の座標を簡単のため、例えば次のようにする（ｊは、所定の自然数）。
Ｗ０＝（０，０，０）
Ｗ１＝（４ｊ，０，０）
Ｗ２＝（０，３ｊ，０）
そして、ベクトルＷ０Ｗ１とベクトルＷ０Ｗ２との外積を頂点Ｗ０を始点として表現した場合における終点をＷ３とする。
Ｗ３＝（０，０，１２ｊ２）
を設ける。
【０１８８】
［第１電子機器１０のキャリブレーション］
第１実施形態と同様に、第１電子機器１０及び座標系算出部３０は、電磁波に基づいて、基準座標系マーカー１０Ａの有する電波受信部（Ｗ０、Ｗ１、Ｗ２）及びＷ３の第１座標系における座標値を算出する。
そうすることで、第１実施形態と同様に、キャリブレーション処理部３２は、式２と同様に、第１座標系における座標値を基準座標系における座標値に変換するための基準座標変換行列１Ｔｗを算出することができる。
【０１８９】
次に、第１実施形態と同様に、第２電子機器２０ｎを予め固定した位置においた状態で、３次元計測部２２ｎの備える第２座標系における座標値を基準座標系における座標値に変換するための座標変換行列（以下、「ｎＴＷ」という）を算出する。
このため、座標系算出部３０は、例えば、マーカー位置情報算出部３１を備えるように構成する。
そうすることで、第２電子機器２０ｎの３次元計測部２２及びマーカー位置情報算出部３１により基準座標系マーカー１０Ａの頂点Ｗ０、Ｗ１、Ｗ２の位置を表す第２座標系における座標値を算出することにより、キャリブレーション処理部３２は、３次元計測部２２ｎの備える第２座標系における座標値を基準座標系における座標値に変換するための座標変換行列ｎＴＷを算出することができる。
【０１９０】
次に、第２電子機器２０ｎを固定した位置においた状態のままで、第２電子機器２０ｎの備える電波受信部座標系による座標値を第１電子機器１０の第１座標系における座標値に変換するための第１座標変換行列ＭｎＴ１を算出する。
具体的には、第１実施形態と同様に、第１電子機器１０は、電波に基づいて、第２電子機器２０ｎの電波受信部２３Ｃｋ（１≦ｋ≦３）及びＣ４の第１座標系における座標値を算出する。
そうすることで、第１実施形態と同様に、キャリブレーション処理部３２は、第２電子機器２０ｎの電波受信部座標系による座標値を第１電子機器１０の第１座標系における座標値に変換するための第１座標変換行列ＭｎＴ１を算出することができる。
【０１９１】
以上から、第１実施形態と同様に、キャリブレーション処理部３２は、基準座標変換行列１ＴＷ、第１座標変換行列ＭｎＴ１、及び第２座標系における座標値を基準座標系における座標値に変換するための座標変換行列ｎＴＷに基づいて、式７と同様に、第２電子機器２０ｎの３次元計測部２２の備える第２座標系における座標値を第２電子機器２０ｎの備える電波受信部２３Ｃｋ（１≦ｋ≦３）による電波受信部座標系における座標値に変換するための電波受信部座標変換行列ｎＴＭを算出することができる。
第１実施形態の場合と同様に、基準座標変換行列１ＴＷ、及び電波受信部座標変換行列ｎＴＭは、それぞれ固定の変換行列である。
【０１９２】
以上のようにキャリブレーションを行い、基準座標変換行列１ＴＷ、及び電波受信部座標変換行列ｎＴＭを算出した後、例えば３次元データ処理部４０の記憶部（図示せず）に記憶するように構成する。その後、基準座標系マーカー１０Ａを現実空間から撤去することができる。
【０１９３】
その後、点群データ座標値変換部４１は、タイムスタンプｔにおいて、点群データの点Ｑｎ（ｔ）の第２座標系（ｎ）での座標が（ｘ，ｙ，ｚ）で表される点に対して、電波受信部座標変換行列ｎＴＭ、第２電子機器２０ｎの電波受信部座標系による座標値を第１電子機器１０の第１座標系における座標値に変換するための第１座標変換行列ＭｎＴ１（ｔ）、及び第１座標系における座標値を基準座標系における座標値に変換するための基準座標変換行列１Ｔｗによる座標変換を行うことで、基準座標系での座標（ｘｇ，ｙｇ，ｚｇ）を算出する（式１９参照）。
【数７】
0006293110.tif 000008
ここで、Ｆｎ（ｔ）＝１ＴＷ・ＭｎＴ１（ｔ）・ｎＴＭである。
このようにして、点群データ座標値変換部４１は、タイムスタンプｔにおける、点群データＱｎ（ｔ）の各点を基準座標系での座標値に変換した点群データ｛Ｑ´ｎ（ｔ）｝を算出することができる。
【０１９４】
［第３実施形態］
第１実施形態又は第２実施形態では、点群データ取得システム１００は、第１電子機器１０と、当該第１電子機器１０の備える位置情報計測部１２が、その計測対象とする、第２電子機器２０とからなる１つのグループ（以下、「電子機器群」という）から構成される。
これに対して、第１電子機器１０と第２電子機器２０とを有する電子機器群が、階層構造を構成するようにすることができる。
【０１９５】
図７に、電子機器群を階層構造とする場合の一例を示す。図７を参照すると、第（ｎ＋１）層（ｎ≧１）に第１電子機器１０を備え、第（ｎ＋１）層の第１電子機器１０は、第ｎ層の第１電子機器１０からみると、第２電子機器２０のように作用する。
【０１９６】
例えば、第１実施形態において、第（ｎ＋１）層（ｎ≧１）の第１電子機器１０は、マーカー２３を備えることにより、第ｎ層の第１電子機器１０が、第（ｎ＋１）層（ｎ≧１）の第１電子機器１０の位置情報を取得することができる。
そうすることで、点群データ座標値変換部４１は、第（ｎ＋１）層の第１電子機器１０の第１計測座標系における座標値を第ｎ層の第１電子機器１０の第１計測座標系における座標値に変換することができる。
このように反復処理することにより、第（ｎ＋１）層以上（ｎ≧１）にある第１電子機器１０の第１計測座標系における座標値及び第２電子機器２０の第２座標系における座標値を、第１層にある第１電子機器１０の基準座標系における座標値に変換することができる。
【０１９７】
第２実施形態において、電磁波伝播時間を用いて第２電子機器２０の位置情報の算出を測定することにより、第１電子機器１０が第２電子機器２０の位置情報を取得する場合、第（ｎ＋１）層（ｎ≧１）の第１電子機器１０は、電波発信部１２Ｐｊ（１≦ｊ≦４）に加えて、電波受信部２３Ｃｋ（１≦ｋ≦３）を備えることで、第ｎ層の第１電子機器１０の電波発信部１２Ｐｊ（１≦ｊ≦４）から送信された電波を、第（ｎ＋１）層（ｎ≧１）の第１電子機器１０の電波受信部２３Ｃｋ（１≦ｋ≦３）が受信することにより、第ｎ層の第１電子機器１０から、第（ｎ＋１）層の第１電子機器１０の位置情報を算出するように構成することができる。
【０１９８】
こうすることで、点群データ座標値変換部４１は、第（ｎ＋１）層の第１電子機器１０の第１座標系における座標値を第ｎ層の第１電子機器１０の第１座標系における座標値に変換することができる。
このように反復処理することにより、第（ｎ＋１）層以上（ｎ≧１）にある第１電子機器１０の第１計測座標系における座標値及び第２電子機器２０の第２座標系における座標値を、第１層にある第１電子機器１０の基準座標系における座標値に変換することができる。
【０１９９】
以上のように、点群データ取得システム１００において、複数の電子機器群を階層的に構成することにより、３次元計測部２２の個数及び計測位置を柔軟に拡張することができる。すなわち、電子機器群の階層を増やすことにより、３次元計測部２２による測定可能な領域を広げることができる。
【０２００】
以上、本発明の点群データ取得システム１００の好ましい実施形態について説明したが、本発明は上述の実施形態に制限されるものではなく、適宜変更が可能である。
【０２０１】
［システム構成の自由度］
本発明の点群データ取得システム１００において、第１電子機器１０、第２電子機器２０、座標系算出部３０、３次元データ処理部４０、及び３次元データ表示制御部５０、及び３次元データ保存部６０、の有する各機能部をそれぞれ、適宜特定のコンピュータに集中させるか、又は分散させることは、ユーザにとって適宜成しえる設計事項である。
例えば、第１実施形態において、座標系算出部３０、３次元データ処理部４０、３次元データ表示制御部５０、３次元データ保存部６０の機能の一部、一部の組み合わせ、又は全部を、例えば、第１電子機器１０の機能としてもよい。
また、座標系算出部３０、３次元データ処理部４０、３次元データ表示制御部５０、３次元データ保存部６０の機能を１つのコンピュータ又は複数のコンピュータで実装することも当業者が適宜設計できる事項である。
【０２０２】
［位置情報算出部と位置情報ユニット］
前述したように、第１実施形態の第１電子機器１０の備える位置情報計測部１０、座標系算出部３０の備えるマーカー位置情報算出部３１を総称して、「位置情報算出部」といい、第２実施形態の第１電子機器１０の備える電波発信部１２Ｐｊ（１≦ｊ≦４）、第２電子機器位置情報算出部３１Ａを総称して、「位置情報算出部」ということとする。
第１実施形態の第２電子機器２０の備えるマーカー２３を総称して、「位置情報ユニット」といい、第２実施形態の第２電子機器２０の備える電波受信部２３Ｃｋ（１≦ｋ≦３）を総称して、「位置情報ユニット」ということとする。
そうすると、本願発明は、「位置情報算出部」が、第２電子機器２０の備える「位置情報ユニット」の位置情報を算出することにより、第２電子機器２０の備える３次元計測部１２による計測した各点群データを共通の基準座標系における座標に変換することを可能とするものであるといえる。
【０２０３】
［自律計測］
上記実施形態においては、第２電子機器２０（３次元計測部２２）は、例えば、カメラマン等により自由に移動可能とされた。しかし、第２電子機器２０（３次元計測部２２）が、例えば、予め設定された直線レールの上を動く場合等所定の直線の上に限定して移動する場合、特に、レールが直線かつ起伏のない場合に走行距離計によりある地点からの走行距離を計測できる場合、電磁波伝播時間、又は位置情報計測部１２等に基づいて位置を算出する替わりに、自律的に位置を算出するように構成してもよい。
同様に、第２電子機器２０（３次元計測部２２）が、例えば、回転検出器を内蔵したリールにワイヤを巻き付け、ワイヤの巻き取り回数と巻き取り角度を計測して位置を算出することができる場合、自律的に位置を算出するように構成してもよい。
同様に、例えば、第２電子機器２０（３次元計測部２２）をコロにより移動させる場合において、回転計をコロに係合し、当該コロをレール等に押圧してコロの回転数及び角度を計測することで、第２電子機器２０（３次元計測部２２）の位置を算出することができる場合、自律的に位置を算出するように構成してもよい。
【０２０４】
上記実施形態では、深度カメラを３次元計測部２２として適用することを前提としたが、本発明はこれに限定されない。レーザースキャナやデジタイザや電子カメラ等であっても、同様に適用できる。
【０２０５】
また、第１実施形態においては、第１電子機器１０又は座標系算出部３０のトリガーによる同期制御により、複数個あるすべての３次元計測部２２に対して計測動作を同期させることにより、同時刻に計測対象物を計測させるように構成するものとしているが、各３次元計測部２２は、各々別の時刻に計測対象物を計測させるように構成してもよい。
【０２０６】
複数の３次元計測部２２は、各々自分のタイミングで対象物を計測して点群データを取り込み、その結果をリアルタイムに、源点群データ記憶部６１内の割り当てられた領域に格納することができる。その際、点群データにタイムスタンプは付与されるが、その時刻は３次元計測部２２毎に付与される。３次元計測部２２によって性能が違うかもしれないし、同じ性能の３次元計測部２２で同じ対象物を見ても、見る角度によって取得する点群データの量が異なれば、データ取り込みのタイミングもずれることもあるからである。
しかし、複数の点群データを使ってデータ合成、マッチング処理を行う段階になると、データ同士の時刻合わせが必要になる。具体的には、測定時刻が必ずしも一致しない場合には、例えば当該時刻の近傍の３次元計測部２２ｎの計測位置により、線形補完することができる。また、複数の点群データ間において、時刻の近いもの同士を取り出して合成、マッチング処理することもできる。
このように構成することで、３次元計測部２２は、第１電子機器１０のトリガーによる同期制御なしに、各々自分のタイミングで計測対象物を計測することができる。
【０２０７】
また、本システムの自然な拡張として、深度センサ以外の特殊カメラの追加が考えられる。サーモグラフィーのように温度データを取得できるもの、あるいは分光センシングカメラ等、スペクトルグラフを解析して、水分量、糖度、癌細胞と正常細胞の分離、血管の強調等、様々な画像情報を取得できるカメラが存在している。そのような特殊なカメラと第２電子機器２０とを組み合わせることによって、点群データに特殊な情報を付加することができる。これは、コンピュータによる外界の認識に必要な「セグメント化（画像の有意味な分離）」を可能にする。
【０２０８】
これらにより、例えば、機械と人体の分離、正常細胞と癌細胞の分離等の認識処理が、リアルな３次元画像を用いて行うことができるようになる。これらは、ロボットによる手術を行う場合等に有用な技術である。
【符号の説明】
【０２０９】
１０  第１電子機器
１０Ａ  基準座標系マーカー
１１  第１通信部
１２  位置情報計測部
１２Ｐ  電波発信部
１５  時刻同期部
１６  距離計測部
２０  第２電子機器
２１  第２通信部
２２  ３次元計測部
２３  マーカー
２３Ｃ  電波受信部
２５  時刻同期部
３０  座標系算出部
３１  マーカー位置情報算出部
３１Ａ  第２電子機器位置情報算出部
３２  キャリブレーション処理部
３３  時刻同期部
４０  ３次元データ処理部
４１  点群データ座標値変換部
４２  合成点群データ作成部
４２１  第１合成点群データ作成部
４２２  点群データマッチング部
５０  ３次元データ表示制御部
５１  視点操作部
５２  データ表示方式選択部
５３  データ再生・保存指示部
５４  ３次元データ表示処理部
６０  ３次元データ保存部
６１  源点群データ記憶部
６２  第２電子機器位置情報記憶部
１００  点群データ取得システム
【図１Ａ】
0006293110.tif 000009
【図１Ｂ】
0006293110.tif 000010
【図１Ｃ】
0006293110.tif 000011
【図１Ｄ】
0006293110.tif 000012
【図１Ｅ】
0006293110.tif 000013
【図１Ｆ】
0006293110.tif 000014
【図１Ｇ】
0006293110.tif 000015
【図２】
0006293110.tif 000016
【図３】
0006293110.tif 000017
【図４】
0006293110.tif 000018
【図５Ａ】
0006293110.tif 000019
【図５Ｂ】
0006293110.tif 000020
【図５Ｃ】
0006293110.tif 000021
【図６】
0006293110.tif 000022
【図７】
0006293110.tif 000023