(57)【要約】【課題】撮影環境の並進及び回転のスケールに依存することなく、世界座標系の距離との線形性を維持して、撮影画像からカメラの位置姿勢の類似度を算出する装置等を提供する。【解決手段】装置は、参照対象が映る撮影画像の画像座標値と当該参照対象の参照画像の参照世界座標値との座標変換に基づいて、当該撮影画像におけるカメラと参照対象との相対的な撮影位置姿勢ベクトル（並進t及び回転θ）を推定する撮影位置姿勢推定手段と、撮影位置姿勢ベクトルを用いて、カメラの世界座標値を表すカメラ視点Ｐcam(x,y,z)から、参照対象の世界座標値Ｐtargetに対するカメラ注視点Ｐlookat(u,v,w)に向けた撮影視線分ベクトルを算出する撮影視線分算出手段と、異なる撮影画像における撮影視線分ベクトルの類似度を、カメラの位置姿勢の類似度として算出する類似度算出手段とを有する。【選択図】図５
2019075018.tif 000002
【特許請求の範囲】
【請求項１】
世界座標値Ｐtargetを既知とする参照対象が映る撮影画像から、異なるカメラの位置姿勢の類似度を算出する装置であって、
参照対象が映る参照画像と、当該参照対象の参照世界座標値とを対応付けて予め記憶した参照対象記憶手段と、
参照対象が映る撮影画像の画像座標値と当該参照対象の参照画像の参照世界座標値との座標変換に基づいて、当該撮影画像におけるカメラと参照対象との相対的な撮影位置姿勢ベクトル（並進t及び回転θ）を推定する撮影位置姿勢推定手段と、
前記撮影位置姿勢ベクトルを用いて、前記カメラの世界座標値を表すカメラ視点Ｐcam(x,y,z)から、参照対象の世界座標値Ｐtargetに対するカメラ注視点Ｐlookat(u,v,w)に向けた撮影視線分ベクトルを算出する撮影視線分算出手段と、
異なる撮影画像における撮影視線分ベクトルの類似度を、カメラの位置姿勢の類似度として算出する類似度算出手段と
を有することを特徴とする装置。
【請求項２】
提示情報と、当該提示情報の提示視線分ベクトルとを対応付けて予め記憶した提示視線分記憶手段と、
前記類似度算出手段で検索された提示視線分ベクトルに対応する提示情報を、前記撮影画像に重畳的に表示する提示情報表示手段と
を更に有し、
前記類似度算出手段は、前記撮影位置姿勢推定手段で用いられた参照対象について、前記撮影視線分ベクトルに最も類似する提示視線分ベクトルを、前記提示視線分記憶手段から検索する
ことを特徴とする請求項１に記載の装置。
【請求項３】
前記提示視線分記憶手段は、同一の参照対象について異なる複数の提示情報それぞれに提示視線分ベクトルを予め記憶している
ことを特徴とする請求項２に記載の装置。
【請求項４】
前記撮影視線分算出手段は、カメラ視点Ｐcamと、当該カメラ視点Ｐcamからの撮影方向に伸びる半直線に対して前記参照対象の世界座標値Ｐtargetから下した垂線との交点であるカメラ注視点Ｐlookatとを結ぶ撮影視線分ベクトルを算出する
ことを特徴とする請求項１から３のいずれか１項に記載の装置。
【請求項５】
前記参照画像の世界座標値Ｐtargetは、対称面Ｓtargetに基づくものであり、
前記撮影視線分算出手段は、カメラ視点Ｐcamと、当該カメラ視点Ｐcamからの撮影方向に伸びる半直線と前記対称面Ｓtargetとの交点であるカメラ注視点Ｐlookatとを結ぶ撮影視線分ベクトルを算出する
ことを特徴とする請求項１から３のいずれか１項に記載の装置。
【請求項６】
前記類似度算出手段における前記撮影視線分ベクトルの類似度は、視線分ベクトルにおけるカメラ視点(x,y,z)間の距離とカメラ注視点(u,v,w)間の距離との和である
ことを特徴とする請求項１から５のいずれか１項に記載の装置。
【請求項７】
前記類似度算出手段における前記撮影視線分ベクトルの類似度は、視線分ベクトル間のユークリッド距離である
ことを特徴とする請求項１から５のいずれか１項に記載の装置。
【請求項８】
前記類似度算出手段における前記撮影視線分ベクトルの類似度は、異なるカメラ視点Ｐcam(x,y,z)間のユークリッド距離と、異なるカメラ注視点Ｐlookat(u,v,w)間のユークリッド距離との和である
ことを特徴とする請求項１から５のいずれか１項に記載の装置。
【請求項９】
前記撮影画像は、時間軸上の所定フレームレートにおける連続的なものであり、
前記類似度算出手段は、時刻tの撮影画像の撮影視線分ベクトルと、類似度算出対象となる提示視線分ベクトルとの類似度に、時刻tの撮影画像の撮影視線分ベクトルと、時刻t-1の撮影画像によって検索された提示視線分ベクトルとの類似度を重み付けて算出する
ことを特徴とする請求項１から８のいずれか１項に記載の装置。
【請求項１０】
前記撮影画像は、時間軸上の所定フレームレートにおける連続的なものであり、
前記類似度算出手段は、時刻tの撮影画像の撮影視線分ベクトルと、類似度算出対象となる提示視線分ベクトルとの類似度に、時刻t-1の撮影画像の撮影視線分ベクトルと、時刻tの撮影画像の撮影視線分ベクトルとの移動量を、時刻t-1の撮影画像によって検索された提示視線分ベクトルに加算して得られる時刻tの提示視線分ベクトルの予想値と、類似度算出対象となる提示視線分ベクトルとの類似度を重み付けて算出する
ことを特徴とする請求項１から８のいずれか１項に記載の装置。
【請求項１１】
世界座標値Ｐtargetを既知とする参照対象が映る撮影画像から、異なるカメラの位置姿勢の類似度を算出するようにコンピュータを機能させるプログラムであって、
参照対象が映る参照画像と、当該参照対象の参照世界座標値とを対応付けて予め記憶した参照対象記憶手段と、
参照対象が映る撮影画像の画像座標値と当該参照対象の参照画像の参照世界座標値との座標変換に基づいて、当該撮影画像におけるカメラと参照対象との相対的な撮影位置姿勢ベクトルを推定する撮影位置姿勢推定手段と、
前記撮影位置姿勢ベクトルを用いて、前記カメラの世界座標値を表すカメラ視点Ｐcam(x,y,z)から、参照対象の世界座標値Ｐtargetに対するカメラ注視点Ｐlookat(u,v,w)に向けた撮影視線分ベクトルを算出する撮影視線分算出手段と、
異なる撮影画像における撮影視線分ベクトルの類似度を、カメラの位置姿勢の類似度として算出する類似度算出手段と
してコンピュータを機能させることを特徴とするプログラム。
【請求項１２】
世界座標値Ｐtargetを既知とする参照対象が映る撮影画像から、異なるカメラの位置姿勢の類似度を算出する装置の位置姿勢類似度算出方法であって、
前記装置は、
参照対象が映る参照画像と、当該参照対象の参照世界座標値とを対応付けて予め記憶した参照対象記憶部を有し、
参照対象が映る撮影画像の画像座標値と当該参照対象の参照画像の参照世界座標値との座標変換に基づいて、当該撮影画像におけるカメラと参照対象との相対的な撮影位置姿勢ベクトルを推定する第１のステップと、
前記撮影位置姿勢ベクトルを用いて、前記カメラの世界座標値を表すカメラ視点Ｐcam(x,y,z)から、参照対象の世界座標値Ｐtargetに対するカメラ注視点Ｐlookat(u,v,w)に向けた撮影視線分ベクトルを算出する第２のステップと、
異なる撮影画像における撮影視線分ベクトルの類似度を、カメラの位置姿勢の類似度として算出する第３のステップと
を実行することを特徴とする装置の位置姿勢類似度算出方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数のカメラの位置姿勢における類似度を算出する技術に関する。
【背景技術】
【０００２】
従来、撮影対象物とカメラとの相対的な位置姿勢に応じて、撮影画像に別のコンテンツをリアルタイムに表示する技術がある（例えば特許文献１参照）。
この技術によれば、撮影対象物に、参照対象としての「マーカ」が貼り付けられており、そのマーカが映る撮影画像を用いて、カメラの位置姿勢を推定する。また、参照対象が映る参照画像毎に、そのカメラの位置姿勢に応じたコンテンツ（提示情報）を予め記憶している。ここでのコンテンツとは、拡張現実(Augmented Reality)画像として、撮影画像に重畳的に表示可能なものである。
そして、参照対象が映り込む撮影画像から、そのカメラの位置姿勢と、予め記憶された複数の参照対象が映る参照画像の位置姿勢とを、「距離関数」に基づいて比較する。距離が短い（位置姿勢が類似する）と推定された参照対象に紐づく提示情報を、撮影画像に重畳的に表示する。
【０００３】
尚、カメラと参照対象との相対的な位置姿勢（並進ベクトルｔ及び回転ベクトルθ）を推定する技術もある（例えば非特許文献１参照）。この技術によれば、参照対象となるマーカの参照画像を予め記憶しており、撮影画像からその参照画像をマッチングによって位置姿勢を推定する。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１６−７１４９６号公報
【非特許文献】
【０００５】
【非特許文献１】加藤, M. Billinghurst, 浅野, 橘, "マーカー追跡に基づく拡張現実感システムとそのキャリブレーション", 日本バーチャルリアリティ学会論文誌, vol.4, 20 no.4, pp.607-617, 1999.、[online]、［平成２９年１０月７日検索］、インターネット＜URL:http://intron.kz.tsukuba.ac.jp/tvrsj/4.4/kato/p-99_VRSJ4_4.pdf＞
【非特許文献２】D.G.Lowe, Distinctive image features from scale-invariant key points, Proc. of Int. Journal of Computer Vision (IJCV), 60(2) pp.91-110 (2004)
【非特許文献３】H.Bay, T.Tuytelaars, and L.V.Gool, SURF: Speed Up Robust Features, Proc. of Int. Conf. of ECCV, (2006)
【非特許文献４】Ethan Rublee, Vincent Rabaud, Kurt Konolige, Gary R. Bradski: ORB: An efficient alternative to SIFT or SURF. ICCV 2011: 2564-2571.
【発明の開示】
【発明が解決しようとする課題】
【０００６】
図１は、従来技術における位置姿勢の評価を表す説明図である。
【０００７】
例えば以下のような距離関数を用いて、この距離が短いほど、カメラの位置姿勢は類似する（近い）と判定される。
位置姿勢ベクトルに基づく第１の距離関数
画像特徴量に基づく第２の距離関数
【０００８】
＜位置姿勢ベクトルに基づく第１の距離関数＞
第１の距離関数として、カメラと参照対象との相対的な位置姿勢を表す並進ベクトルｔ及び回転ベクトルθのユークリッド距離を用いる。
【０００９】
図１（ａ）によれば、世界座標系上で、カメラCAMに対して、カメラcam1又はcam2のいずれの位置姿勢が近いか、を判定している。
CAMとcam1との並進の差：Δｔ1＝100mm
CAMとcam1との回転の差：Δθ1＝45°
CAMとcam2との並進の差：Δｔ2＝120mm
CAMとcam1との回転の差：Δθ2＝10°
重みｗ(＝1）を付けて、回転ｔと並進θとを加算すると、例えば以下のようになる。
Ｄ1＝Δｔ1＋ｗ・Δθ1＝100＋1*45＝145
Ｄ2＝Δｔ2＋ｗｖΔθ2＝120＋1*10＝130
この場合、Ｄ1＞Ｄ2となり、カメラCAMの位置姿勢には、カメラcam2の位置姿勢の方が近いと推定される。
一方で、撮影環境のスケールが変わって、例えば並進ｔが10000mmオーダとなった場合、以下のようになる。
Ｄ1＝Δｔ1＋ｗ・Δθ1＝10000＋45＝10045
Ｄ2＝Δｔ2＋ｗ・Δθ2＝12000＋10＝12010
この場合、Ｄ1＜Ｄ2となり、カメラCAMの位置姿勢には、カメラcam1の位置姿勢の方が近いと推定される。
【００１０】
このように、撮影環境の並進ｔ及び回転θのスケールに応じて、パラメータに与える相対的な重みｗを調整する必要がある。また、そもそも、重みｗをどのように適切に決定するかも不明となる。
【００１１】
＜画像特徴量に基づく第２の距離関数＞
第２の距離関数として、局所領域（画像特徴量）同士の画像座標のユークリッド距離を用いる。
【００１２】
画像座標でのユークリッド距離と、世界座標でのユークリッド距離との関係は、非線形となる。画像特徴量の画像座標での距離が小さいにも関わらず、世界座標での位置姿勢が大きく異なる場合がある。即ち、撮影画像における横ずれや拡大縮小の変化があると、距離が遠い（類似度が低い）と推定される。
【００１３】
図１（ｂ）によれば、カメラCAMの位置姿勢に最も近い（類似度が高い）カメラを推定しようとしている。
カメラcam1は、画像内での拡大による影響で、画像特徴量間の距離が長い（類似度が低い）と推定される。
カメラcam2は、画像内での平行移動による影響で、画像特徴量間の距離が長い（類似度が低い）と推定される。
一方で、カメラcam3は、画像内で画像特徴量が幾何的に大きく歪んでいるにも関わらず、画像特徴量間の距離が短い（類似度が高い）と推定される。
結果的に、図１（ｂ）によれば、カメラCAMには、カメラcam3が、位置姿勢としては最も近いと推定される。このとき、実際の位置姿勢から乖離した形で、提示情報（コンテンツ）を撮影画像に重畳させた場合、ユーザに対して不自然な印象を与えてしまう。
【００１４】
これに対し、本願の発明者らは、カメラの位置姿勢の類似度評価について、以下の２点の問題点を除去することができないか、と考えた。
（１）撮影対象物との相対的なカメラの位置姿勢を表すベクトルのユークリッド距離を用いる場合に、撮影環境における並進及び回転のスケールへの依存性の問題
（２）画像特徴量の画像座標でのユークリッド距離を用いる場合に、世界座標系における距離との非線形の問題
【００１５】
そこで、本発明は、撮影環境における並進及び回転のスケールに依存することなく、世界座標系における距離との線形性を維持して、撮影画像からカメラの位置姿勢の類似度を算出する装置、プログラム及び方法を提供することを目的とする。
【課題を解決するための手段】
【００１６】
本発明によれば、世界座標値Ｐtargetを既知とする参照対象が映る撮影画像から、異なるカメラの位置姿勢の類似度を算出する装置であって、
参照対象が映る参照画像と、当該参照対象の参照世界座標値とを対応付けて予め記憶した参照対象記憶手段と、
参照対象が映る撮影画像の画像座標値と当該参照対象の参照画像の参照世界座標値との座標変換に基づいて、当該撮影画像におけるカメラと参照対象との相対的な撮影位置姿勢ベクトル（並進t及び回転θ）を推定する撮影位置姿勢推定手段と、
撮影位置姿勢ベクトルを用いて、カメラの世界座標値を表すカメラ視点Ｐcam(x,y,z)から、参照対象の世界座標値Ｐtargetに対するカメラ注視点Ｐlookat(u,v,w)に向けた撮影視線分ベクトルを算出する撮影視線分算出手段と、
異なる撮影画像における撮影視線分ベクトルの類似度を、カメラの位置姿勢の類似度として算出する類似度算出手段と
を有することを特徴とする。
【００１７】
本発明の装置における他の実施形態によれば、
提示情報と、当該提示情報の提示視線分ベクトルとを対応付けて予め記憶した提示視線分記憶手段と、
類似度算出手段で検索された提示視線分ベクトルに対応する提示情報を、撮影画像に重畳的に表示する提示情報表示手段と
を更に有し、
類似度算出手段は、撮影位置姿勢推定手段で用いられた参照対象について、撮影視線分ベクトルに最も類似する提示視線分ベクトルを、提示視線分記憶手段から検索することも好ましい。
【００１８】
本発明の装置における他の実施形態によれば、
提示視線分記憶手段は、同一の参照対象について異なる複数の提示情報それぞれに提示視線分ベクトルを予め記憶していることも好ましい。
【００１９】
本発明の装置における他の実施形態によれば、
撮影視線分算出手段は、カメラ視点Ｐcamと、当該カメラ視点Ｐcamからの撮影方向に伸びる半直線に対して参照対象の世界座標値Ｐtargetから下した垂線との交点であるカメラ注視点Ｐlookatとを結ぶ撮影視線分ベクトルを算出することも好ましい。
【００２０】
本発明の装置における他の実施形態によれば、
参照画像の世界座標値Ｐtargetは、対称面Ｓtargetに基づくものであり、
撮影視線分算出手段は、カメラ視点Ｐcamと、当該カメラ視点Ｐcamからの撮影方向に伸びる半直線と対称面Ｓtargetとの交点であるカメラ注視点Ｐlookatとを結ぶ撮影視線分ベクトルを算出することも好ましい。
【００２１】
本発明の装置における他の実施形態によれば、
類似度算出手段における撮影視線分ベクトルの類似度は、視線分ベクトルにおけるカメラ視点(x,y,z)間の距離とカメラ注視点(u,v,w)間の距離との和であることも好ましい。
【００２２】
本発明の装置における他の実施形態によれば、
類似度算出手段における撮影視線分ベクトルの類似度は、視線分ベクトル間のユークリッド距離であることも好ましい。
【００２３】
本発明の装置における他の実施形態によれば、
類似度算出手段における撮影視線分ベクトルの類似度は、異なるカメラ視点Ｐcam(x,y,z)間のユークリッド距離と、異なるカメラ注視点Ｐlookat(u,v,w)間のユークリッド距離との和であることも好ましい。
【００２４】
本発明の装置における他の実施形態によれば、
撮影画像は、時間軸上の所定フレームレートにおける連続的なものであり、
類似度算出手段は、時刻tの撮影画像の撮影視線分ベクトルと、類似度算出対象となる提示視線分ベクトルとの類似度に、時刻tの撮影画像の撮影視線分ベクトルと、時刻t-1の撮影画像によって検索された提示視線分ベクトルとの類似度を重み付けて算出することも好ましい。
【００２５】
本発明の装置における他の実施形態によれば、
撮影画像は、時間軸上の所定フレームレートにおける連続的なものであり、
類似度算出手段は、時刻tの撮影画像の撮影視線分ベクトルと、類似度算出対象となる提示視線分ベクトルとの類似度に、時刻t-1の撮影画像の撮影視線分ベクトルと、時刻tの撮影画像の撮影視線分ベクトルとの移動量を、時刻t-1の撮影画像によって検索された提示視線分ベクトルに加算して得られる時刻tの提示視線分ベクトルの予想値と、類似度算出対象となる提示視線分ベクトルとの類似度を重み付けて算出することも好ましい。
【００２６】
本発明によれば、世界座標値Ｐtargetを既知とする参照対象が映る撮影画像から、異なるカメラの位置姿勢の類似度を算出するようにコンピュータを機能させるプログラムであって、
参照対象が映る参照画像と、当該参照対象の参照世界座標値とを対応付けて予め記憶した参照対象記憶手段と、
参照対象が映る撮影画像の画像座標値と当該参照対象の参照画像の参照世界座標値との座標変換に基づいて、当該撮影画像におけるカメラと参照対象との相対的な撮影位置姿勢ベクトルを推定する撮影位置姿勢推定手段と、
撮影位置姿勢ベクトルを用いて、カメラの世界座標値を表すカメラ視点Ｐcam(x,y,z)から、参照対象の世界座標値Ｐtargetに対するカメラ注視点Ｐlookat(u,v,w)に向けた撮影視線分ベクトルを算出する撮影視線分算出手段と、
異なる撮影画像における撮影視線分ベクトルの類似度を、カメラの位置姿勢の類似度として算出する類似度算出手段と
してコンピュータを機能させることを特徴とする。
【００２７】
本発明によれば、世界座標値Ｐtargetを既知とする参照対象が映る撮影画像から、異なるカメラの位置姿勢の類似度を算出する装置の位置姿勢類似度算出方法であって、
装置は、
参照対象が映る参照画像と、当該参照対象の参照世界座標値とを対応付けて予め記憶した参照対象記憶部を有し、
参照対象が映る撮影画像の画像座標値と当該参照対象の参照画像の参照世界座標値との座標変換に基づいて、当該撮影画像におけるカメラと参照対象との相対的な撮影位置姿勢ベクトルを推定する第１のステップと、
撮影位置姿勢ベクトルを用いて、カメラの世界座標値を表すカメラ視点Ｐcam(x,y,z)から、参照対象の世界座標値Ｐtargetに対するカメラ注視点Ｐlookat(u,v,w)に向けた撮影視線分ベクトルを算出する第２のステップと、
異なる撮影画像における撮影視線分ベクトルの類似度を、カメラの位置姿勢の類似度として算出する第３のステップと
を実行することを特徴とする。
【発明の効果】
【００２８】
本発明の装置、プログラム及び方法によれば、撮影環境における並進及び回転のスケールに依存することなく、世界座標系における距離との線形性を維持して、撮影画像からカメラの位置姿勢の類似度を算出することができる。
【図面の簡単な説明】
【００２９】
【図１】従来技術における位置姿勢の評価を表す説明図である。
【図２】本発明における仮想現実の提示情報を表示する説明図である。
【図３】本発明における装置の機能構成図である。
【図４】本発明におけるフローチャートである。
【図５】本発明におけるカメラ注視点の決定方法を表す説明図である。
【図６】図１に対応した本発明の位置姿勢の類似度を表す説明図である。
【図７】斜方撮影における本発明の位置姿勢の類似度を表す説明図である。
【発明を実施するための形態】
【００３０】
以下では、図面を用いて、本発明を実施するための形態について詳細に説明する。
【００３１】
図２は、本発明における仮想現実の提示情報を表示する説明図である。
【００３２】
図２（ａ）によれば、撮影対象物に、マーカ（参照対象）が貼り付けられている。マーカは、例えば２次元的なＱＲコード（登録商標）であってもよい。その撮影対象物は、位置姿勢が異なる４台のカメラcam1、cam2、cam3、cam4によって撮影されている。
図２（ｂ）によれば、各カメラの撮影画像が表されている。撮影画像にマーカが映り込む場合、そのマーカは、カメラ毎に異なる角度から映る。本発明によれば、マーカが映り込む撮影画像から、そのカメラの位置姿勢が抽出される。
図２（ｃ）によれば、カメラの位置姿勢に応じた提示情報（コンテンツ）が表されている。コンテンツが立体形状の場合、カメラの位置姿勢に応じて、異なる形状の提示情報が選択される。
図２（ｄ）によれば、仮想現実画像には、撮影画像に映るマーカの上に、例えば直方体のコンテンツが重畳的に表示される。
図２（ｅ）によれば、ユーザは、カメラcam1,cam2,cam4によって撮影された撮影画像について、撮影対象物にコンテンツが重畳されているように視認することができる。
【００３３】
図３は、本発明における装置の機能構成図である。
図４は、本発明におけるフローチャートである。
【００３４】
図３（ａ）によれば、本発明の装置１は、例えばスマートフォンのような、カメラとディスプレイとを搭載した端末として表されている。勿論、スマートフォンに限られず、タブレットやパーソナルコンピュータであってもよい。
図３（ｂ）によれば、本発明の装置１は、サーバとして表されている。サーバは、通信インタフェースを介して、ユーザ端末のカメラによって撮影された撮影画像を受信し、提示情報（コンテンツ）を重畳した仮想現実画像をそのユーザ端末へ返信する。
尚、本発明を１台の装置のみで実現することなく、本発明の機能を区分して、ネットワークを介したサーバ及び端末のシステムとして構成することもできる。
【００３５】
装置１は、参照対象記憶部１０１と、提示視線分記憶部１０２と、撮影位置姿勢推定部１１と、撮影視線分算出部１２と、類似度算出部１３と、提示情報表示部１４とを有する。これら機能構成部は、装置に搭載されたコンピュータを機能させるプログラムを実行することによって実現される。また、これら機能構成部の処理の流れは、カメラの位置姿勢の類似度を算出する方法としても理解できる。
【００３６】
［参照対象記憶部１０１］
参照対象記憶部１０１は、参照対象が映る「参照画像」と、当該参照対象の「参照世界座標値」とを対応付けて予め記憶したものである。
図４によれば、参照対象記憶部１０１のテーブルには、参照画像と参照世界座標とが対応付けて記憶されている。
【００３７】
参照画像は、回転やスケールの変化にロバストな、例えばＳＩＦＴ(Scale-Invariant Feature Transform)やＳＵＲＦ(Speeded Up Robust Features)の特徴量であってもよいしコンパクトな局所特徴を高速に抽出するＯＲＢ(Oriented FAST and Rotated BRIEF)の特徴量であってもよい（例えば非特許文献２、３、４参照）。
例えばＳＩＦＴの場合、１枚の画像からは１２８次元の局所特徴の集合が抽出される。ＳＩＦＴとは、スケールスペースを用いて特徴的な局所領域を解析し、そのスケール変化及び回転に不変となる局所特徴を記述する技術である。
また、ＳＵＲＦの場合、ＳＩＦＴよりも高速処理が可能であって、１枚の画像から６４次元の局所特徴の集合が抽出される。ＳＩＦＴは、処理コストが高く且つリアルタイムマッチングが困難であるのに対し、ＳＵＲＦは、積分画像を利用することによって処理を高速化している。
更に、ＯＲＢの場合、スマートフォンやタブレットのようなモバイル端末が普及する中で、拡張現実感の用途について、省メモリ化及び高速マッチング化に適する。
【００３８】
［撮影位置姿勢推定部１１］
撮影位置姿勢推定部１１は、最初に、カメラで撮影された撮影画像から、局所特徴量抽出アルゴリズムによって、参照対象記憶部１０１に記憶された参照画像を抽出する。
局所特徴量抽出アルゴリズムは、前述したＳＩＦＴ、ＳＵＲＦ、ＯＲＢのような周知の既存技術である。具体的には、撮影画像の中から、相対的な輝度勾配に基づく局所特徴量を用いて、参照画像の特徴量に類似する局所領域を抽出する。抽出される局所領域は、回転、拡大縮小、射影変化（射影変換による歪み）のいずれか又はその任意の組み合わせに対して、不変な性質を有する。
但し、参照画像を抽出できない場合、その撮影画像に重畳的に表示される提示情報（コンテンツ）はなく、提示情報表示部１４からその撮影画像のままで出力される。
【００３９】
そして、撮影位置姿勢推定部１１は、参照対象が映る「撮影画像の画像座標値」と、当該参照対象の「参照世界座標値」との座標変換に基づいて、当該撮影画像におけるカメラと参照対象との相対的な「撮影位置姿勢ベクトル」を推定する。
「撮影位置姿勢ベクトル」は、撮影対象物に貼り付けられたマーカ（参照対象）を基準として、カメラの相対的な位置姿勢（並進ｔ及び回転θ）である。
そして、撮影位置姿勢ベクトルは、撮影視線分算出部１２へ出力される。
【００４０】
［撮影視線分算出部１２］
撮影視線分算出部１２は、マーカ（参照対象）の世界座標値Ｐtargetを、ユーザによって又は別途の機能によって予め登録している。
そして、撮影視線分算出部１２は、撮影位置姿勢ベクトルを用いて、カメラの世界座標値を表すカメラ視点Ｐcam(x,y,z)から、参照対象の世界座標値Ｐtargetに対するカメラ注視点Ｐlookat(u,v,w)に向けた「撮影視線分ベクトル」を算出する。視線分ベクトルは、(x,y,z,u,v,w)の６要素からなる。
そして、撮影視線分算出部１２は、算出した視線分ベクトル(x,y,z,u,v,w)を、類似度算出部１３へ出力する。
【００４１】
図５は、本発明におけるカメラ注視点の決定方法を表す説明図である。
【００４２】
＜カメラ注視点の第１の決定方法＞
図５（ａ）によれば、撮影視線分算出部１２は、カメラ視点Ｐcam(x,y,z)と、当該カメラ視点Ｐcam(x,y,z)からの撮影方向に伸びる半直線に対して参照対象の世界座標値Ｐtarget(u,v,w)から下した垂線との交点であるカメラ注視点Ｐlookat(u,v,w)とを結ぶ「撮影視線分ベクトル(x,y,z,u,v,w)」を算出する。
＜カメラ注視点の第２の決定方法＞
図５（ｂ）によれば、参照画像の世界座標値Ｐtargetは、対称面Ｓtarget（平面）に基づくものであるとする。ここで、撮影視線分算出部１２は、カメラ視点Ｐcam(x,y,z)と、当該カメラ視点Ｐcam(x,y,z)からの撮影方向に伸びる半直線と対称面Ｓtargetとの交点であるカメラ注視点Ｐlookat(u,v,w)とを結ぶ「撮影視線分ベクトル(x,y,z,u,v,w)」を算出する。
【００４３】
［提示視線分記憶部１０２］
提示視線分記憶部１０２は、「提示視線分ベクトル」と、「提示情報」とを予め記憶したものである。
図４によれば、提示視線分記憶部１０２は、同一の参照対象について異なる複数の提示情報それぞれに提示視線分ベクトルを予め記憶している
「提示視線分ベクトル」は、提示情報（コンテンツ）に対する、カメラの視線分ベクトルである。これは、提示情報について、前述した撮影視線分算出部１２と同様に算出されたものである。
【００４４】
［類似度算出部１３］
類似度算出部１３は、異なる撮影画像における撮影視線分ベクトルの類似度を、カメラの位置姿勢の類似度として算出する。
類似度算出部１３は、撮影位置姿勢推定部１１で用いられた参照対象について、撮影視線分ベクトルに最も類似する提示視線分ベクトルを、提示視線分記憶部１０２から検索する。このとき、類似度算出部１３は、撮影視線分ベクトルCと提示視線分ベクトルCi'との類似度ｄ(Ci',C)を算出する。
i_min＝argimin ｄ(Ci',C)
C：撮影視線分ベクトル
Ci'：検索対象となる提示視線分ベクトル(i∈N、検索対象となる総数)
Ci_min'：撮影視線分ベクトルに最も類似する提示視線分ベクトル
argimin ｄ  ：ｄを最小化するiの集合
そして、類似度算出部１３は、検索された提示視線分ベクトルCi_min'に対応する提示情報を、提示情報表示部１４へ出力する。
【００４５】
撮影画像は、時間軸上の所定フレームレートにおける連続的なものであってもよい。この場合、時刻範囲Δt毎に、カメラの撮影画像I(t)が入力され、その都度、リアルタイムに、撮影位置姿勢推定部１１、撮影視線分算出部１２、類似度算出部１３及び提示情報表示部１４が実行される。
i_min(t)＝argimin ｄ(Ci',C(t))
【００４６】
ここで、類似度算出における異なる実施形態について説明する。
＜類似度算出の第１の実施形態＞
類似度算出部１３における撮影視線分ベクトルの類似度は、視線分ベクトルにおけるカメラ視点(x,y,z)間の距離とカメラ注視点(u,v,w)間の距離との和である。
ｄ(C',C)＝ｒ＋ｇ
＝｜Ｐ'cam(x,y,z)−Ｐcam(x,y,z)｜
＋｜Ｐ'lookat(u,v,w)−Ｐlookat(u,v,w)｜
ｒ：カメラ視点(x,y,z)間の距離｜Ｐ'cam(x,y,z)−Ｐcam(x,y,z)｜
ｇ：カメラ注視点(u,v,w)間の距離｜Ｐ'lookat(u,v,w)−Ｐlookat(u,v,w)｜
【００４７】
＜類似度算出の第２の実施形態＞
類似度算出部１３における撮影視線分ベクトルの類似度は、視線分ベクトル間のユークリッド距離である。即ち、６次元空間における距離として算出することができる。
ｄ(C',C)＝√((x'-x)2＋(y'-y)2＋(z'-z)2＋(u'-u)2＋(v'-v)2＋(w'-w)2)
【００４８】
＜類似度算出の第３の実施形態＞
類似度算出部１３における撮影視線分ベクトルの類似度は、異なるカメラ視点Ｐcam(x,y,z)間のユークリッド距離と、異なるカメラ注視点Ｐlookat(u,v,w)間のユークリッド距離との和である。
ｄ(C',C)＝√((x'-x)2＋(y'-y)2＋(z'-z)2)＋√((u'-u)2＋(v'-v)2＋(w'-w)2)
【００４９】
＜類似度算出の第４の実施形態＞
撮影画像が、時間軸上の所定フレームレートにおける連続的なものである場合、類似度算出部１３は、時刻tの撮影画像の撮影視線分ベクトルと、類似度算出対象となる提示視線分ベクトルとの類似度に、時刻tの撮影画像の撮影視線分ベクトルと、時刻t-1によって検索された提示視線分ベクトルとの類似度を重み付けて算出する。
【００５０】
時系列の撮影画像をリアルタイムに処理する場合、時刻t毎に、撮影画像の撮影視線分ベクトルC(t)に最も類似する提示視線分ベクトルCi_min(t)'を検索する。一方で、提示視線分ベクトルCi'は、不規則かつ離散的な位置姿勢に基づいて予め記憶されたものである。そのために、最も類似すると判定された提示視線分ベクトルCi_min(t)'は、直近の過去時刻t-1における既に判定された提示視線分ベクトルCi_min(t-1)'から大きく変動したものとなる可能性がある。そのような突発的な変動は、時刻tにおける過去の撮影視線分ベクトルC(t-1)から、突発的に乖離した形で提示視線分ベクトルCi_min(t)'を選択することとなる。このとき、時系列上で連続してリアルタイムに提示視線分ベクトルを選択する場合、鋸歯状に変動する選択結果となってしまうこともある。
【００５１】
そのために、時刻t毎に判定される提示視線分ベクトルCi_min(t)'が、時間軸上で滑らかに変動するように、直近の過去時刻t-1における既に判定された提示視線分ベクトルCi_min(t-1)'から、時刻tにおける提示視線分ベクトルCi_min(t)'が極端に変動することを抑制するべく、コスト項を追加した以下のような式とする。
i_min(t)＝argimin｛αｄ(Ci',C(t))＋(1-α)ｄ(Ci',Ci_min(t-1)')｝
α：0＜α＜1
【００５２】
＜類似度算出の第５の実施形態＞
類似度算出部１３は、時刻t毎に判定される提示視線分ベクトルCi_min(t)'が、時間軸上で滑らかに変動するように、直近の過去時刻t-1における既に判定された提示視線分ベクトルCi_min(t-1)'から、時刻tにおける提示視線分ベクトルCi_min(t)'が極端に変動することを抑制するべく、時刻t-1の撮影画像の撮影視線分ベクトルC(t-1)と、時刻tの撮影画像の撮影視線分ベクトルC(t)との移動量(C(t)-C(t-1))を、Ci_min(t-1)'に加算して得られるCi_min(t)'の予想値と、類似度算出対象となる提示視線分ベクトルとの類似度を、コスト項として追加した以下のような式とする。
i_min(t)＝argimin｛αd(Ci',C(t))＋
(1−α)ｄ(Ci',Ci_min(t-1)'+{C(t)−C(t-1)})｝
α：0<α<1
これによって、時刻tの提示視線分ベクトルCi_min(t)'の変動を平滑化することができる。この場合、提示視線分記憶部１０２の中で、平滑化された提示視線分ベクトルCi_min(t)'に対応する提示情報が選択される。
【００５３】
［提示情報表示部１４］
提示情報表示部１４は、類似度算出部１３で検索された提示視線分ベクトルに対応する提示情報を、撮影画像に重畳的に表示する。
これによって、撮影画像にコンピュータグラフィックでコンテンツを描画することによって、高品位な仮想現実画像をユーザに提供することができる。
【００５４】
図６は、図１に対応した本発明の位置姿勢の類似度を表す説明図である。
【００５５】
図６（ａ）は、図１（ａ）と比較して、回転を、カメラ注視点の座標として扱うことによって、回転と並進との間の重み付けが不要となる。
カメラCAMの視線分ベクトルとカメラcam1の視線分ベクトルとの類似度は、例えば以下のように算出される。
Ｄ1＝ｒ1＋ｇ1＝100mm＋300mm＝400mm
カメラCAMの視線分ベクトルとカメラcam2の視線分ベクトルとの類似度は、例えば以下のように算出される。
Ｄ2＝ｒ2＋ｇ2＝120mm＋100mm＝220mm
この場合も、Ｄ1＞Ｄ2となり、カメラcam2の位置姿勢の方が近いと推定される。
一方で、撮影環境のスケールが変わって、例えば並進ｔが10000mmオーダとなっても、以下のようになる。
Ｄ1＝ｒ1＋ｇ1＝10000＋30000＝40000
Ｄ2＝ｒ2＋ｇ2＝12000＋10000＝22000
この場合も、Ｄ1＞Ｄ2となり、カメラcam2の位置姿勢の方が近いと推定され、位置姿勢の比較に変化がなく、スケールを考慮しなくても、カメラcam2の位置姿勢の方が近いと推定される。
【００５６】
図６（ｂ）は、図１（ｂ）と比較して、カメラCAMに対して、cam1及びcam2が同程度の類似度であり、cam3は類似度が低い（距離が長い）と判定される。本発明によれば、拡大縮小、横ずれ、視線軸に対する回転の影響によって、位置姿勢の類似度が低くなることがない。
【００５７】
図７は、斜方撮影における本発明の位置姿勢の類似度を表す説明図である。
【００５８】
図７によれば、カメラCAMとカメラcam1、cam2、cam3の類似度の変化は、回転角の変化に対して、線形に近い関係となる。カメラのアングルによって回転の類似度への寄与率が変わりづらく、正対していても回転変化への感度が悪くならない。
【００５９】
以上、詳細に説明したように、本発明の装置、プログラム及び方法によれば、撮影環境における並進及び回転のスケールに依存することなく、世界座標系における距離との線形性を維持して、撮影画像からカメラの位置姿勢の類似度を算出することができる。
【００６０】
前述した本発明における種々の実施形態によれば、当業者は、本発明の技術思想及び見地の範囲における種々の変更、修正及び省略を容易に行うことができる。前述の説明はあくまで例であって、何ら制約しようとするものではない。本発明は、特許請求の範囲及びその均等物として限定するものにのみ制約される。
【符号の説明】
【００６１】
１  装置
１０１  参照対象記憶部
１０２  提示視線分記憶部
１１  撮影位置姿勢推定部
１２  撮影視線分算出部
１３  類似度算出部
１４  提示情報表示部
【図１】
2019075018.tif 000003
【図２】
2019075018.tif 000004
【図３】
2019075018.tif 000005
【図４】
2019075018.tif 000006
【図５】
2019075018.tif 000007
【図６】
2019075018.tif 000008
【図７】
2019075018.tif 000009