JP7442072B2 - 三次元変位計測方法及び三次元変位計測装置 - Google Patents

Description

本開示は、三次元変位計測方法及び三次元変位計測装置に関する。

コンピュータビジョンの分野における三次元再構成技術では、複数の二次元画像間で特徴点の対応付けを行い、カメラの位置、向き、及び被写体の三次元位置を推定する。また、三次元再構成技術では、カメラキャリブレーション及び三次元点群再構成が行われる。

例えば、特許文献１に記載の装置は、３台以上の複数カメラ間でキャリブレーションを行い、取得したカメラパラメータにより各カメラ座標系を任意視点の仮想カメラ座標系へ変換する。当該装置は、その仮想カメラ座標系において、座標変換後の画像間のブロックマッチングによる対応付けを行い、距離情報を推定する。

特開２０１０－２５０４５２号公報

このような、上記従来技術では、複数の視点から被写体を撮像することにより得られた複数の画像に基づいて、三次元空間における被写体の経時的な変位を効果的に計測することができない。

そこで、本開示は、複数の視点から被写体を撮像することにより得られた複数の画像に基づいて、三次元空間における被写体の経時的な変位を効果的に計測することができる技術を提供する。

本開示の一態様に係る三次元変位計測方法は、Ｎ（Ｎは１より大きい整数）台のカメラを含む複数のカメラによって互いに異なるＭ（ＭはＮより大きい整数）個の視点において撮像されたＭ枚の校正用画像に基づいて、前記Ｎ台のカメラのカメラパラメータを算出するカメラ校正ステップと、（１）異なる時刻毎に前記Ｎ台のカメラのそれぞれによって被写体を撮像することで生成されるＮ枚の画像、及び、（２）前記カメラパラメータ、に基づいて、三次元空間における前記被写体の経時的な変位を計測する変位計測ステップと、を含む。

本開示によれば、複数の視点から被写体を撮像することにより得られた複数の画像に基づいて、三次元空間における被写体の経時的な変位を効果的に計測することができる。

図１は、実施の形態における三次元変位計測が行われる際の概略図である。 図２は、実施の形態における三次元モデリングの原理を示す図である。 図３は、実施の形態における三次元変位計測の原理を示す図である。 図４Ａは、実施の形態におけるカメラの撮影周期と露光時間の一例を示す図である。 図４Ｂは、実施の形態におけるカメラの撮影周期と露光時間の別の一例を示す図である。 図５は、実施の形態における三次元変位計測装置のブロック図である。 図６は、実施の形態における三次元再構成システムのフローチャートである。 図７は、実施の形態における制御部のブロック図である。 図８は、実施の形態における制御部の動作を表すフローチャートである。 図９は、実施の形態における三次元変位計測の際に用いる基線長、撮像対象とカメラとの距離、及び、輻輳角を示す図である。 図１０は、実施の形態におけるカメラ校正に用いる視点位置の決定方法を示す図である。 図１１は、実施の形態における多視点フレームセットの例を示す図である。 図１２は、実施の形態における三次元画像処理部のブロック図である。 図１３は、実施の形態における三次元画像処理部の動作を示すフローチャートである。 図１４は、実施の形態に係る三次元変位計測装置が自由視点映像生成を行う際の動作を示すフローチャートである。 図１５は、自由視点映像生成を行う際の三次元画像処理部の動作を示すフローチャートである。 図１６は、三次元画像処理部の変形例のブロック図である。 図１７は、三次元画像処理部の変形例の動作を示すフローチャートである。 図１８は、三次元空間認識の概要を示す図である。

（本開示の基礎となった知見）
従来技術では、被写体の三次元形状を再構成（モデリング）することで、または再構成結果を用いることで三次元空間認識が実現される。被写体の三次元形状を再構成する三次元空間再構成装置は、同一の場面の映像を撮影する複数のカメラを含む撮像システムから提供される映像データと、校正によって得られた各カメラの位置、姿勢等（以下、「位置姿勢」という。）を示すカメラパラメータと、を用いてモデリングを行う。カメラの位置は視点に対応する。

また、構造物（例えばインフラ用の構造物）、機械部品等といった計測対象の変位を三次元空間において把握することが求められている。しかしながら、従来技術では、計測対象である被写体の三次元空間における経時的な変位を効果的に計測することができなかった。例えば、従来技術では、被写体の三次元空間における経時的な変位を計測するための処理負荷の低減は不十分であった。

一般的に、カメラ校正処理、及び、三次元変位計測処理では、視点数が多い、つまり、画像が多いほど処理負荷が大きくなる一方で精度が向上するというトレードオフの関係がある。このトレードオフの関係に対し、発明者らは、三次元変位計測の精度よりもカメラ校正の精度を優先すると、処理負荷を抑えつつ２つの処理全体としての精度が向上することを見出した。

そこで、本開示の一態様に係る三次元変位計測方法は、Ｎ（Ｎは１より大きい整数）台のカメラを含む複数のカメラによって互いに異なるＭ（ＭはＮより大きい整数）個の視点において撮像されたＭ枚の校正用画像に基づいて、前記Ｎ台のカメラのカメラパラメータを算出するカメラ校正ステップと、（１）異なる時刻毎に前記Ｎ台のカメラのそれぞれによって被写体を撮像することで生成されるＮ枚の画像、及び、（２）前記カメラパラメータ、に基づいて、三次元空間における前記被写体の経時的な変位を計測する変位計測ステップと、を含む。

これにより、三次元変位計測方法では、カメラパラメータの精度が向上するように三次元変位計測処理における視点数Ｎよりも多い視点数Ｍを、カメラ校正処理に用いることで、三次元変位計測処理における精度を向上させることができる。

また、例えば、本開示の一態様に係る三次元変位計測方法において、前記変位計測ステップは、前記Ｎ枚の画像に含まれる第１画像に含まれる、前記被写体の第１対応点と、前記Ｎ枚の画像に含まれる第２画像に含まれ且つ前記第１対応点に対応する第２対応点と、を求める対応点算出ステップと、前記カメラパラメータと、前記第１対応点と、前記第２対応点とに基づいて、前記三次元空間における前記被写体の三次元点を求める三次元点算出ステップと、前記異なる時刻に対応する複数の前記三次元点から、前記三次元点の経時的な変位を求める変位算出ステップと、を含んでもよい。

これにより、三次元変位計測方法では、例えば画像マッチングの手法及び三角測量の原理を用いて、カメラパラメータ、第１対応点及び第２対応点に基づいて、三次元点の経時的な変位を求めることで、被写体の三次元空間における変位を求めることができる。

また、例えば、本開示の一態様に係る三次元変位計測方法は、前記Ｎは、前記カメラパラメータと、前記変位の計測の精度と、に基づいて決定されてもよい。

これにより、三次元変位計測方法では、変位計測ステップで行われる計測の精度を満たすように、カメラの視点数Ｎを定めることができる。よって、三次元変位計測方法は、精度よく被写体の三次元空間における変位を求めることができる。

また、本願発明者らは、カメラ校正に用いる視点を均等に配置するよりも、以下のように配置すれば、カメラ校正の精度が向上することを見出した。例えば、本開示の一態様に係る三次元変位計測方法は、前記Ｎ台のカメラは、互いに隣り合う第１カメラ及び第２カメラを含み、前記Ｍ個の視点は、前記第１カメラが配置される第１位置と、前記第２カメラが配置される第２位置との間の第１中点を含み、前記第１位置と前記第１中点との間の中点を第２中点とすると、前記Ｍ個の視点に含まれ且つ前記第２中点と前記第１位置との間に設けられる視点は、前記Ｍ個の視点に含まれ且つ前記第１中点と前記第２中点との間に設けられる視点よりも多くてもよい。

すなわち、この態様では、カメラ校正に使用する視点を、三次元変位計測に使用する視点の位置を基準として第２中点よりも近い範囲に、第２中点よりも遠い範囲よりも多く設定する。このような不均等な配置により、カメラ校正の精度が向上する。

例えば、カメラ校正のための視点の不均等な配置では、前記カメラ校正ステップでは、前記第１中点と前記第２中点との間には、視点を設けなくてもよい。

また、例えば、本開示の一態様に係る三次元変位計測装置は、Ｎ（Ｎは１より大きい整数）台のカメラを含む複数のカメラによって互いに異なるＭ（ＭはＮより大きい整数）個の視点において撮像されたＭ枚の校正用画像に基づいて前記Ｎ台のカメラのカメラパラメータを算出するカメラ校正部と、（１）異なる時刻毎に前記Ｎ台のカメラのそれぞれによって被写体を撮像することで生成されるＮ枚の画像、及び、（２）前記カメラパラメータ、に基づいて、三次元空間における前記被写体の経時的な変位を計測する変位計測部と、を備える。

これにより、三次元変位計測装置は、本開示の一態様に係る三次元変位計測方法と同様の効果を奏することができる。

以下、実施の形態について図面を参照しながら具体的に説明する。

なお、以下で説明する実施の形態は、いずれも包括的または具体的な例を示すものである。以下の実施の形態で示される数値、形状、材料、構成要素、構成要素の配置位置及び接続形態、ステップ、ステップの順序などは、一例であり、請求の範囲を限定する主旨ではない。また、以下の実施の形態における構成要素のうち、最上位概念を示す独立請求項に記載されていない構成要素については、任意の構成要素として説明される。また、各図は、必ずしも厳密に図示したものではない。各図において、実質的に同一の構成については同一の符号を付し、重複する説明は省略又は簡略化する。

（実施の形態）
本実施の形態に係る三次元変位計測装置は、時刻間で座標軸の一致した時系列三次元モデルを再構成できる。ここで、本開示における三次元再構成を定義する。実空間上に存在する被写体を複数のカメラにより異なる視点で撮影した映像又は画像を多視点映像又は多視点画像と呼ぶ。つまり、多視点画像は、同一の被写体を異なる視点から撮影した複数の二次元画像を含む。また、時系列に撮影された多視点画像を多視点映像と呼ぶ。この多視点画像を用いて被写体を三次元空間に再構成することを三次元再構成と呼ぶ。

具体的には、まず、三次元変位計測装置は、時刻毎に独立して三次元再構成を行うことで、各時刻の三次元モデルを取得する。次に、三次元変位計測装置は、静止カメラ及び静止物体（静止三次元点）を検出し、検出した静止カメラ及び静止物体を用いて、時刻間で三次元モデルの座標合せを行い、座標軸の一致した時系列三次元モデルを生成する。

これにより、三次元変位計測装置は、カメラの固定／非固定又は被写体の移動／静止に関わらず、各時刻の被写体及びカメラの相対位置関係が高精度であり、かつ時間方向の推移情報を利用可能な時系列三次元モデルを生成できる。

図１は、実施の形態における三次元変位計測が行われる際の概略図である。実施の形態における三次元変位計測方法において、被写体１００１を、カメラ１００－１及びカメラ１００－２が撮影している。ここで、被写体１００１は、所定の位置に固定されているコンプレッサ等でもよい。また、カメラ１００－１及びカメラ１００－２は、互いに異なる位置において被写体１００１を撮影範囲に含むような向きで固定されている固定カメラでもよい。例えば、校正済みの１つまたは複数のカメラ（例えば固定カメラ）によって同一空間または同一の被写体１００１が、多視点から撮影されることにより、撮影される空間または被写体１００１が三次元再構成されることができる。すなわち、多視点から撮影された複数の二次元画像の間で、特徴点等の対応付けが行われることで、カメラの位置、カメラの向き、並びに、被写体１００１の三次元変位または三次元位置が推定されることが可能である。つまり、三次元変位計測装置では、多視点から撮影された複数の二次元画像を用いて、カメラ校正及び被写体１００１の三次元変位計測等が行われる。三次元変位計測装置では、さらに、被写体１００１の三次元再構成が行われてもよい。

三次元変位計測は、例えば、インフラ構造物点検システム等に有用である。例えば、計測装置が新幹線等の防音壁等の壁の振動を、壁の奥行方向に向かって計測するというニーズがある。また、計測装置がコンプレッサ等の動作において振動が発生する機器の変位を計測するというニーズがある。

図２は、実施の形態における三次元モデリングの原理を示す図である。まず、本実施の形態における三次元モデリングを定義する。本実施の形態では、三次元変位計測装置が、実空間上に存在する被写体を複数のカメラによって異なる視点から撮影した二次元の多視点映像を用いて、被写体を三次元空間内に再構成することを三次元モデリングという。なお、以下では、三次元モデリングを三次元再構成とも言うことがある。

図２に示されるように、少なくとも２つの視点のカメラ座標及びカメラパラメータを用いることで、点線で示された三角形が形成される。実施の形態における三次元変位計測方法では、各カメラで撮影された画像の画像面上の座標系で表された画像面上の各点は、この三角形を利用した三角測量の原理に基づいて、世界座標系を用いた三次元モデルとして再構成される。画像面上の座標系は、例えば、図２に示される画像面座標１または画像面座標２である。三次元モデルは、被写体が、世界座標系を用いて表される三次元空間内に再構成（つまり、三次元再構成）されたものである。三次元モデルは、多視点から撮影された複数の二次元画像に映る被写体上の複数の点それぞれの三次元位置を表している。ここで、三次元位置は、例えば、三次元空間における直交座標系（ＸＹＺ軸）で示される、Ｘ成分、Ｙ成分及びＺ成分からなる三値の情報で表される。なお、三次元モデル上の各点に対応する情報は、上記の各点の三次元位置を表す情報だけでなく、各点の色、並びに、各点及びその周辺の表面形状を表す情報を含んでもよい。

三次元モデリングに際して、各カメラのカメラパラメータは、予め取得されていてもよいし、三次元モデルの作成と同時に推定されてもよい。ここで、カメラパラメータは、内部パラメータ及び外部パラメータからなる。内部パラメータは、カメラの焦点距離、画像中心等を含み、外部パラメータは、カメラの三次元位置、カメラの三次元における向き等を含む。

図２では、代表的なピンホールカメラモデルを例として用いている。この例では、カメラのレンズ歪みは考慮されていない。三次元モデリングに際して、レンズ歪みが考慮される場合は、三次元モデルに用いられる各点の位置として、各カメラで撮影された画像の画像面上の座標系における点の位置を歪みモデルにより正規化した補正位置が用いられる。

図３は、実施の形態における三次元変位計測の原理を示す図である。本開示の三次元変位計測方法について定義を行う。多視点から撮影された二次元画像を用いて被写体１００１の三次元空間における変位を算出することを三次元変位計測という。図３では、三次元変位計測の仕組みが示される。少なくとも２つの視点におけるそれぞれのカメラ座標及びカメラパラメータを用いることで、図３において点線で示されたような三角形が形成される。この三角形が利用される三角測量の原理を元に、複数の視点において撮影され、画像面上の座標系で表された画像上の点が、世界座標系において再構成される。画像面上の座標系は、例えば、図３に示される画像面座標１または画像面座標２である。

被写体１００１の三次元モデルは、多視点から撮影された複数の二次元画像に映る被写体１００１上の複数の点それぞれの三次元位置により示される。被写体１００１が移動または振動する場合、時刻ｔ及び時刻ｔ＋ｉのそれぞれにおいて、被写体１００１上の三次元点の位置（Ｘ_ｔ，Ｙ_ｔ，Ｚ_ｔ）及び（Ｘ_ｔ＋ｉ，Ｙ_ｔ＋ｉ，Ｚ_ｔ＋ｉ）を再構成すると、時間ｉの間の変位（Ｘ_ｔ＋ｉ－Ｘ_ｔ，Ｙ_ｔ＋ｉ－Ｙ_ｔ，Ｚ_ｔ＋ｉ－Ｚ_ｔ）が算出される。三次元変位計測において、例えば、奥行き方向（Ｚ成分）のみを計測することとしてもよい。すなわち、三次元変位計測において、いわゆるデプスマップを作成することとしてもよい。また、本開示の三次元変位計測は、Ｘ成分、Ｙ成分、及びＺ成分のうち任意の２成分における変位を計測してもよい。

なお、例えば、被写体１００１がコンプレッサ等の振動体である場合、時刻差ｉは、振動体の共振周波数をサンプリングできる周波数の逆数である。しかし、時刻差ｉは、これに限定されない。例えば、被写体１００１が構造物であり、構造物のクラックの成長を計測する場合、時刻差ｉは、構造物に求められる耐久年数に応じて決定される。

また、上述の例では、時刻ｔ及び時刻ｔ＋ｉにおいて三次元点の位置が再構成された。しかし、三次元点位置の再構成回数は、２回に限らず、時刻差ｉ毎に３回以上であってもよい。さらに、三次元点位置の再構成回数が時刻差ｉ毎に３回以上の場合、本開示の三次元変位計測は、三次元点の変位だけでなく、変位の速度及び加速度を算出してもよい。

三次元変位計測における具体的な計算は以下の通りである。まず、カメラ座標系１において、カメラの内部パラメータを表す行列をＫ_１とし、外部パラメータのうち平行移動を表す並進ベクトルをＴ_１とし、回転を表す回転行列をＲ_１とする。また、カメラ座標系２において、カメラの内部パラメータを表す行列をＫ_２とし、外部パラメータのうち平行移動を表す並進ベクトルをＴ_２とし、回転を表す回転行列をＲ_２とする。

次に、画像マッチングにより、被写体１００１上のある点について、画像面座標１と画像面座標２との間の対応点を算出する。ここで、画像マッチングでは、画像の中の矩形領域であるブロックの差分値が最小となる位置を探索するブロックマッチングが行われてもよい。具体的には、ブロックマッチングとして、それぞれの画像に含まれる対応付けられた矩形領域の一致度を、輝度差の総和（ＳＡＤ：Ｓｕｍ ｏｆ Ａｂｓｏｌｕｔｅ Ｄｉｆｆｅｒｅｎｃｅ）、または、輝度差の２乗和（ＳＳＤ：Ｓｕｍ ｏｆ Ｓｑｕａｒｅｄ Ｄｉｆｆｅｒｅｎｃｅ）といった一致度評価関数によって評価し、一致度評価値が最大または最小になる位置を探索する。一致度評価関数として、相互相関（ＣＣ：Ｃｒｏｓｓ Ｃｏｒｒｅｌａｔｉｏｎ）または正規化相互相関（ＮＣＣ：Ｎｏｒｍａｌｉｚｅｄ Ｃｒｏｓｓ Ｃｏｒｒｅｌａｔｉｏｎ）が用いられてもよい。また、画像マッチングでは、ＳＩＦＴ（Ｓｃａｌｅ－Ｉｎｖａｒｉａｎｔ Ｆｅａｔｕｒｅ Ｔｒａｎｓｆｏｒｍ）、ＳＵＲＦ（Ｓｐｅｅｄｅｄ Ｕｐ Ｒｏｂｕｓｔ Ｆｅａｔｕｒｅｓ）等の、画像の小領域から特量ベクトルを算出し、その差分値が最小となる位置を探索する特徴量マッチングが行われてもよい。

画像マッチングの際には、時刻ｔにおける、画像面座標１における対応点をａ_１，ｔ＝（ｘ_１，ｔ，ｙ_１，ｔ）、画像面座標２における対応点をｂ_２，ｔ＝（ｘ_２，ｔ，ｙ_２，ｔ）とする。画像面座標１における対応点ａ_１，ｔは、第１対応点の具体例であり、画像面座標２における対応点ｂ_２，ｔは、第２対応点の具体例である。また、時刻ｔ＋ｉにおける、画像面座標１における対応点をａ_{１，ｔ＋ｉ}＝（ｘ_{１，ｔ＋ｉ}，ｙ_{１，ｔ＋ｉ}）、画像面座標２における対応点をｂ_{２，ｔ＋ｉ}＝（ｘ_{２，ｔ＋ｉ}，ｙ_{２，ｔ＋ｉ}）とする。なお、ここで説明された処理は、対応点算出ステップに対応する。

そして、カメラの内部パラメータ及び外部パラメータを用いて、ａ_１，ｔ、ｂ_２，ｔ、ａ_{１，ｔ＋ｉ}、ｂ_{２，ｔ＋ｉ}のそれぞれに対応する三次元位置を表す三次元点を算出する。時刻ｔのａ_１，ｔに相当する三次元点をＡ_１，ｔ、時刻ｔのｂ_２，ｔに相当する三次元点をＢ_２，ｔ、時刻ｔ＋ｉのａ_{１，ｔ＋ｉ}に相当する三次元点をＡ_{１，ｔ＋ｉ}、時刻ｔ＋ｉのｂ_{２，ｔ＋ｉ}に相当する三次元点をＢ_{２，ｔ＋ｉ}とする。

二次元画像上でのそれぞれの対応点と、それぞれの対応点に相当する三次元点との変換は、以下の式（１）～（４）によって計算される。

ａ_１，ｔ＝Ｋ_１［Ｒ_１｜Ｔ_１］Ａ_１，ｔ 式（１）
ｂ_２，ｔ＝Ｋ_２［Ｒ_２｜Ｔ_２］Ｂ_２，ｔ 式（２）
ａ_{１，ｔ＋ｉ}＝Ｋ_１［Ｒ_１｜Ｔ_１］Ａ_{１，ｔ＋ｉ} 式（３）
ｂ_{２，ｔ＋ｉ}＝Ｋ_２［Ｒ_２｜Ｔ_２］Ｂ_{２，ｔ＋ｉ} 式（４）

このとき、二次元画像上での対応点に相当する三次元点Ａ_１，ｔ及びＢ_２，ｔは、世界座標系で同一の点であるが、式（１）及び式（２）を用いて算出された三次元点Ａ_１，ｔ及びＢ_２，ｔは必ずしも一致しない。同様に、三次元点Ａ_{１，ｔ＋ｉ}及びＢ_{２，ｔ＋ｉ}は、世界座標系で同一の点であるが、式（３）及び式（４）を用いて算出された三次元点Ａ_{１，ｔ＋ｉ}及びＢ_{２，ｔ＋ｉ}は必ずしも一致しない。その場合、式（５）及び式（６）に示されるように、三次元点Ａ_１，ｔ及びＢ_２，ｔの間の中点、または、三次元点Ａ_{１，ｔ＋ｉ}及びＢ_{２，ｔ＋ｉ}の間の中点を算出することで、二次元画像上での対応点に相当する三次元点を算出してもよい。

（Ｘ_ｔ，Ｙ_ｔ，Ｚ_ｔ）＝（Ａ_１，ｔ＋Ｂ_２，ｔ）／２ 式（５）
（Ｘ_ｔ＋ｉ，Ｙ_ｔ＋ｉ，Ｚ_ｔ＋ｉ）＝（Ａ_{１，ｔ＋ｉ}＋Ｂ_{２，ｔ＋ｉ}）／２ 式（６）

なお、二次元画像上での対応点に相当する三次元点が一致しない場合の、三次元点の決定の仕方は、これに限定されない。三次元点Ａ_１，ｔ及びＢ_２，ｔのうちの一方、または、三次元点Ａ_{１，ｔ＋ｉ}及びＢ_{２，ｔ＋ｉ}のうちの一方のみを、それぞれの時刻に対応する二次元画像上での対応点に相当する三次元点として利用することとしてもよい。また、式（５）及び式（６）では、２つの三次元点の間の中点を算出する方法を示したが、他の算出方法であってもよい。例えば、三次元点Ａ_１，ｔ及びＢ_２，ｔを加重平均することで得られた時刻ｔにおける三次元点、または、三次元点Ａ_{１，ｔ＋ｉ}及びＢ_{２，ｔ＋ｉ}を加重平均することで得られた時刻ｔ＋ｉにおける三次元点を、それぞれの時刻に対応する三次元点として用いてもよい。なお、ここで説明された処理は、三次元点算出ステップに対応する。

なお、算出される三次元点は、被写体の特定の点に限らない。三次元点は、対応点から求められる限り、被写体の任意の点である。

以上より、時刻ｔから時刻ｔ＋ｉまでの間の、三次元変位Ｖ_{ｔ，ｔ＋ｉ}は、式（７）により計算される。

Ｖ_{ｔ，ｔ＋ｉ}＝（Ｘ_ｔ＋ｉ－Ｘ_ｔ，Ｙ_ｔ＋ｉ－Ｙ_ｔ、Ｚ_ｔ＋ｉ－Ｚ_ｔ） 式（７）

なお、ここで説明された処理は、変位算出ステップに対応する。

また、上記の三次元変位計測の方法は、一例であり、これに限定されない。

図４Ａは、実施の形態におけるカメラの撮影周期と露光時間の一例を示す図である。また、図４Ｂは、実施の形態におけるカメラの撮影周期と露光時間の別の一例を示す図である。図４Ａ及び図４Ｂの横方向は時間を示し、矩形信号が立っている時間はカメラ１００－１またはカメラ１００－２が露光していることを示す。カメラ１００－１またはカメラ１００－２により画像を取得する際、シャッタが開放されている時間を露光時間と呼ぶ。

露光時間中、レンズを通して撮像素子にさらされたシーンが画像として得られる。図４Ａでは、視点の異なる２台のカメラ１００－１、１００－２で撮影されたフレームでは、露光時間が重複している。これにより２台のカメラ１００－１、１００－２により取得したフレームは、同一時刻のシーンを含んでいる同期フレームと判定される。

一方、図４Ｂでは、２台のカメラ１００－１、１００－２で露光時間の重複が無いため、２台のカメラ１００－１、１００－２により取得したフレームは、同一時刻のシーンを含まない非同期フレームと判定される。図４Ａのように、同期フレームを複数のカメラで撮影することを同期撮影と呼ぶ。

次に、本実施の形態に係る三次元変位計測装置の構成を説明する。図５は、実施の形態における三次元変位計測装置のブロック図である。図５に示す三次元再構成システム３は、複数のカメラ１００－１～１００－ｎと、三次元変位計測装置２０００と、制御部３００とを備える。また、三次元変位計測装置２０００は、受信部２１０と、記憶部２２０と、取得部２３０と、三次元画像処理部２４０と、送信部２５０とを備える。

複数のカメラ１００－１～１００－ｎは撮影空間１０００にある被写体１００１を撮影し、撮影された複数の映像である多視点映像を受信部２１０へ送信する。多視点映像の送信は、インターネットなどの公衆通信網、又は専用通信網のいずれを介してもよい。あるいは、多視点映像は、一度ハードディスクドライブ（ＨＤＤ）、ソリッドステートドライブ（ＳＳＤ）などの外部記憶装置に記憶され、必要な時に三次元変位計測装置２０００に入力されてもよい。あるいは、多視点映像は、クラウドサーバ等の外部記憶装置に一旦ネットワークを介して送信され、記憶される。そして、必要な時に三次元変位計測装置２０００へ送信されてもよい。

また、ｎ台のカメラ１００－１～１００－ｎの各々は、監視カメラなどの固定カメラ、または、固定されていない非固定カメラである。つまり、ｎ台のカメラ１００－１～１００－ｎは、例えば、それぞれ、互いに異なる位置において、互いに異なる姿勢で固定されている固定カメラであってもよく、または、ビデオカメラ、スマートフォンもしくはウェアラブルカメラなどのモバイルカメラであってもよく、撮影機能付きドローンなどの移動カメラであってもよい。なお、ｎは、２以上の整数である。

また、多視点映像には、映像又はフレームのヘッダ情報として、撮影したカメラを特定するカメラＩＤなどのカメラ特定情報が付加されてもよい。

複数のカメラ１００－１～１００－ｎを用いて、毎フレームで同じ時刻の被写体を撮影する同期撮影が行われてもよい。あるいは、複数のカメラ１００－１～１００－ｎに内蔵された時計の時刻を合せ、同期撮影せずに、映像又はフレーム毎に、撮影された画像に撮影時刻情報が付加されてもよいし、撮影順序を示すインデックス番号が付加されてもよい。

多視点映像の映像セット毎、映像毎、又はフレーム毎に、同期撮影されたか、非同期撮影されたかを示す情報がヘッダ情報として付加されてもよい。

また、制御部３００は、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）を用いて実現され、記憶部２２０に記憶されるプログラムを読み出して実行する。制御部３００は、三次元画像処理部２４０におけるカメラ校正、三次元モデリング及び三次元変位計測のうちの１つ以上の処理を実施するために、三次元再構成システム３によって用いられる、視点数及び視点位置を決定する。ここで、視点数とは、カメラ校正、三次元モデリング及び三次元変位計測のうちの１つ以上の処理で用いられる二次元画像を撮影する位置の数のことであり、具体的には、カメラの台数のことであってもよいし、固定カメラ及び移動カメラによって撮影が行われる地点の数のことであってもよい。また、視点位置とは、カメラ校正、三次元モデリング及び三次元変位計測のうちの１つ以上の処理で用いられる二次元画像を撮影する位置のことであり、具体的にはカメラの位置のことである。

そして、受信部２１０は、カメラ１００－１～１００－ｎが撮影した多視点映像を受信し、受信した多視点映像を記憶部２２０に入力する。

記憶部２２０は、入力された多視点映像を記憶する。また、記憶部２２０には、撮影空間１０００が撮影された映像データ、この映像データに付されるタイムスタンプなどのメタ情報、カメラ１００－１～１００－ｎのカメラパラメータ、及び、適用中のフレームレートまたは解像度等の撮影設定が一時的に記憶されている。また、記憶部２２０には、制御部３００が読み出して実行するプログラムが記憶されている。記憶部２２０は、ＲＡＭ、ＤＲＡＭ、ＳＲＡＭ等の揮発性メモリでもよい。また、記憶部２２０は、ＲＯＭ、フラッシュメモリ等の不揮発性メモリで実現されてもよい。

取得部２３０は、記憶部２２０に記憶された複数の多視点映像からフレームを選択し、選択したフレームを多視点フレームセットとして三次元画像処理部２４０に入力する。

三次元画像処理部２４０は、入力された多視点フレームセットを用いて、カメラ校正を行い、カメラ校正の結果を用いて、被写体の三次元モデリングまたは被写体の三次元変位計測を行う。三次元画像処理部２４０は、カメラ校正を行った結果としてカメラパラメータを算出し、送信部２５０に出力する。

送信部２５０は、三次元画像処理部２４０が算出したカメラパラメータ、被写体の三次元モデル及び被写体の三次元変位のうちの少なくとも１つを、三次元変位計測装置２０００の外部へ送信する。

受信部２１０、取得部２３０、三次元画像処理部２４０及び送信部２５０は、それぞれプロセッサ及びメモリによって実現されてもよい。また、受信部２１０、取得部２３０、三次元画像処理部２４０及び送信部２５０は、専用の回路によって実現されてもよい。

図６は、実施の形態における三次元再構成システムのフローチャートである。

まず、制御部３００は、カメラ校正、三次元モデリングまたは三次元変位計測のうちの１つ以上の処理を実施するために、三次元再構成システム３によって用いられる、視点数及び視点位置を決定する（ステップＳ２１００）。視点数及び視点位置の決定の詳細は後述される。

次に、カメラ１００－１～１００－ｎは、制御部３００が決定した視点数及び視点位置に基づき、被写体の多視点映像を撮影する（ステップＳ２２００）。

そして、受信部２１０は、カメラ１００－１～１００－ｎのうちの複数のカメラで撮影した多視点映像を受信し、記憶部２２０に入力する（ステップＳ２３００）。

続いて、記憶部２２０は、入力された多視点映像を記憶する（ステップＳ２４００）。

次に、取得部２３０は、制御部３００が決定した視点数及び視点位置に基づいて、記憶された多視点映像からフレームを選択し、選択したフレームを多視点フレームセットとして三次元画像処理部２４０に入力する（ステップＳ２５００）。

ここで、多視点フレームセットは、カメラ１００－１～１００－ｎが撮影したすべての映像から、それぞれ１フレームずつ選択された複数のフレームによって構成されてもよい。また、多視点フレームセットは、カメラ１００－１～１００－ｎが撮影したすべての映像から、少なくとも１フレーム選択された複数のフレームにより構成されてもよい。また、多視点フレームセットは、多視点映像のうちから２つ以上の映像を選択し、選択された各映像から１フレームずつ選択された複数のフレームによって構成されてもよい。また、多視点フレームセットは、多視点映像のうちから２つ以上の映像を選択し、選択された各映像から少なくとも１フレーム選択された複数のフレームによって構成されてもよい。

なお、多視点フレームセットの各フレームにカメラ特定情報が付加されていない場合は、各フレームのヘッダ情報に個別にカメラ特定情報を付加してもよいし、多視点フレームセットのヘッダ情報に一括して、カメラ特定情報を付加してもよい。

また、多視点フレームセットの各フレームに撮影時刻または撮影順を示すインデックス番号が付加されていない場合、各フレームのヘッダ情報に個別に撮影時刻または撮影順を示すインデックス番号を付加してもよい。また、同様の場合、多視点フレームセットのヘッダ情報に一括して撮影時刻または撮影順を示すインデックス番号を付加してもよい。

そして、三次元画像処理部２４０は、入力された多視点フレームセットを用いて、カメラ校正を行い、カメラ校正の結果を用いて三次元モデリングまたは三次元変位計測を行う（ステップＳ２６００）。三次元モデリングまたは三次元変位計測の詳細は、後述される。

また、ステップＳ２５００及びＳ２６００の処理は、多視点フレームセット毎に繰り返し行われる。ここで、繰り返し行われるステップＳ２５００及びステップＳ２６００において、カメラ校正が一回だけ行われることとしてもよい。その場合、最初に一回実施されたカメラ校正の結果を、繰り返し行われるステップＳ２５００及びステップＳ２６００で、繰り返し使用してもよい。

続いて、送信部２５０は、カメラパラメータ、被写体の三次元モデルまたは被写体の三次元変位を、三次元変位計測装置２０００の外部へ送信する（ステップＳ２７００）。

ここで、三次元再構成システム３は、処理を終了する。

なお、ステップＳ２６００で算出されたカメラパラメータ、被写体の三次元モデルまたは被写体の三次元変位は、ステップＳ２６００より後の処理で利用されてもよい。

なお、ステップＳ２７００で出力されたカメラパラメータ、被写体の三次元モデルまたは被写体の三次元変位は、記憶部２２０に記憶されてもよい。

図７は、実施の形態における制御部のブロック図である。制御部３００は、計測用決定部３１０及びカメラ校正用決定部３２０を備える。計測用決定部３１０及びカメラ校正用決定部３２０は、プロセッサ及びメモリで実現されてもよいし、専用の論理回路で実現されてもよい。

計測用決定部３１０は、三次元再構成システム３が被写体１００１の三次元変位計測を行う際に用いる二次元画像を撮影する視点の視点数、視点位置等を決定する。

また、カメラ校正用決定部３２０は、三次元再構成システム３がカメラ校正を行う際に用いる二次元画像を撮影する視点の視点数、視点位置等を決定する。

なお、三次元再構成システム３が、カメラ校正、三次元モデリング及び三次元変位計測を行うために用いる視点数は、それぞれ異なってもよい。また、制御部３００は、視点位置のみを決定してもよい。すなわち、計測用決定部３１０は、三次元再構成システム３が被写体１００１の三次元変位計測を行う際に用いる二次元画像を撮影する視点の視点位置のみを決定してもよい。また、カメラ校正用決定部３２０は、三次元再構成システム３がカメラ校正を行う際に用いる二次元画像を撮影する視点の視点位置のみを決定してもよい。

なお、制御部３００は、三次元変位計測装置２０００にカメラが導入される前に、カメラの視点数及び視点位置を決定してもよい。また、カメラは、制御部３００が決定したカメラの視点数及び視点位置に応じて、導入されてもよい。

なお、制御部３００は、三次元変位計測に用いられる視点数Ｎと、カメラ校正に用いられる視点数Ｍとを、Ｍ＞Ｎの関係を満たすように決定する。また、さらに、制御部３００は、三次元変位計測に用いられる視点が、カメラ校正に用いられる視点に含まれる（すなわち２個の視点が一致する）ように、三次元変位計測に用いられる視点とカメラ校正に用いられる視点とを決定する。

カメラ１００－１～１００－ｎ、または、取得部２３０は、三次元変位計測に用いられるＮ個の視点として、Ｎ台の固定カメラを使用する。換言すれば、三次元変位計測のために、Ｍ－Ｎ個の視点においてカメラによる撮像はされない。また、カメラ校正に用いられるＭ個の視点のうち、三次元変位計測に用いられるＮ個の視点以外の（Ｍ－Ｎ）個の視点として、（Ｍ－Ｎ）台の固定カメラを使用してもよく、または、（Ｍ－Ｎ）台以下の移動カメラを使用してもよい。（Ｍ－Ｎ）台以下の移動カメラを使用する場合は、（Ｍ－Ｎ）台以下の移動カメラは、移動しながら撮影を行うことで、（Ｍ－Ｎ）個以上の視点からの撮影を実現する。また、カメラ校正に用いられるＭ個の視点のうち、三次元変位計測に用いられるＮ個の視点以外の（Ｍ－Ｎ）個の視点は、固定カメラと移動カメラとを組み合わせて実現されてもよい。

次に、制御部３００の動作を説明する。図８は、実施の形態における制御部の動作を表すフローチャートである。

まず、計測用決定部３１０は、三次元変位計測の視点数及び視点位置を決定する（ステップＳ３１００）。

次に、カメラ校正用決定部３２０は、カメラ校正の視点数及び視点位置を決定する（ステップＳ３２００）。このとき、カメラ校正用決定部３２０は、三次元変位計測の視点数Ｎより、カメラ校正の視点数Ｍが多くなるように設定し、かつ、三次元変位計測の視点が、カメラ校正の視点に含まれるように設定する。

ここで、制御部３００は動作を終了する。なお、ステップＳ３１００において、計測用決定部３１０は、三次元モデリングのための視点数及び視点位置の決定を行ってもよい。

次に、制御部３００における視点の決定方法の一例を説明する。まず、図９を用いて、三次元変位計測のための視点数Ｎの決定について説明する。図９は、実施の形態における三次元変位計測の際に用いる基線長、撮像対象とカメラとの距離、及び、輻輳角を示す図である。図９では、撮像対象である被写体１００１と、ｎ台のカメラ１００－１～１００－ｎのうちの２台のカメラ１００－１、１００－２との位置関係が示されている。

制御部３００は、三次元変位計測に用いられる視点数をＮ、被写体１００１とカメラ１００－１、１００－２との距離をｄ、三次元モデリングまたは三次元変位計測の精度をΔｚと設定する。これらの値は、任意に決定される。視点数Ｎは、三次元変位計測に使用されることが可能であるカメラの台数でもよいし、設定可能な範囲で最大の数でもよい。以下の説明では、視点数Ｎがカメラ１００の台数ｎに等しいものとする。ただし、視点数Ｎはカメラ１００の台数ｎと等しくなくてもよい。また、被写体１００１とカメラ１００－１、１００－２との距離ｄは、撮影空間１０００の大きさ、カメラ１００－１、１００－２によって撮影される二次元画像への被写体１００１の映り方等のうちの少なくとも１つから決定されてもよい。ここで、距離ｄは、図９に示されるように、２台のカメラ１００－１、１００－２の位置を通る線から被写体１００１までの最短距離である。ただし、距離ｄは、図９に示す例に限らない。例えば、距離ｄは、被写体１００１から２台のカメラ１００－１、１００－２までの距離の平均値であってもよい。制御部３００は、三次元モデリングまたは三次元変位計測の精度Δｚ、カメラの焦点距離ｆ、カメラの素子間隔ΔＨ、被写体１００１と２台のカメラ１００－１、１００－２との距離ｄを用いて、２台のカメラ１００－１、１００－２間の間隔である基線長Ｌを、式（８）で算出する。

Ｌ＝（ｄ×ｄ×ΔＨ）／ｆ×Δｚ 式（８）

例えば、三次元変位計測に用いられる視点数Ｎについて、Ｎ＝２の場合は、三次元変位計測に使用されるカメラの台数は２台である。このとき、制御部３００は、基線長Ｌ及び被写体とカメラとの距離ｄを、２台のカメラの視点位置を決定するための情報とする。

三次元変位計測に用いられる視点数Ｎが２より大きいとき、制御部３００は以下の計算を行う。まず、制御部３００は、２台のカメラそれぞれの光軸の輻輳角をθとし、θを、基線長Ｌ及び距離ｄを用いて、式（９）で算出する。

θ＝２ａｒｃｔａｎ（Ｌ／２ｄ） 式（９）

また、制御部３００は、値Ｐを式（１０）で算出する。ここで、値Ｐは、輻輳角θに対応して定まる、三次元再構成システム３で利用されるカメラの台数を示す値である。

Ｐ＝ｆｌｏｏｒ（２π／θ） 式（１０）

ＰがＮより大きい場合、制御部３００は、Ｎを式（１１）で算出する。

Ｎ＝Ｐ 式（１１）

以上のように、制御部３００は、カメラパラメータ及び求められる精度Δｚに基づいて、視点数Ｎを決定する。

また、基線長Ｌ、距離ｄ、及び、輻輳角θを、カメラの視点位置を決定するための情報としてもよい。

なお、上記で説明した被写体１００１及び２台のカメラ１００－１、１００－２の位置関係は、２台のカメラ１００－１、１００－２に限るものではなく、ｎ台のカメラ１００－１～１００－ｎのうちの他の任意の２台の組み合わせについても同様である。

図１０は、実施の形態におけるカメラ校正に用いる視点位置の決定方法を示す図である。なお、図１０は、鉛直方向に沿って複数のカメラを見た図である。まず、制御部３００は、カメラ校正に用いられる視点の視点間の輻輳角φを任意に設定する。ここで、θ＞２π―（Ｎ－１）×θの場合、制御部３００は、カメラ校正に用いられる視点数Ｍを、式（１２）で算出する。

Ｍ＝Ｎ＋ｆｌｏｏｒ（θ／φ）×（Ｎ－１） 式（１２）

ここで、三次元変位計測で使用するＮ個の視点のそれぞれの間に、ｆｌｏｏｒ（θ／φ）個ずつ視点が追加される。追加されるｆｌｏｏｒ（θ／φ）個の視点は、図１０に示されるように、三次元変位計測で使用するＮ個の視点のそれぞれの間に追加される。例えば、制御部３００は、視点を、視点位置ｐ１と視点位置ｐ２との中間位置ｃに１個または２個追加する。ここで、制御部３００は、視点位置ｐ１と中間位置ｃとの間を２等分した位置を位置ｐ３とし、視点位置ｐ２と中間位置ｃとの間を２等分した位置を位置ｐ４とする。そして、制御部３００は、視点位置ｐ１と位置ｐ３との間に、中間位置ｃに配置された視点数をｆｌｏｏｒ（θ／φ）個の視点から引いた数の視点を配置する。同様に、制御部３００は、視点位置ｐ２と位置ｐ４との間に、中間位置ｃに配置された視点数を、ｆｌｏｏｒ（θ／φ）個の視点から引いた数の視点を配置する。このような配置方法によれば、例えば、追加する視点のうち大多数の視点が、三次元変位計測で使用する視点に対応する位置に配置されたカメラ同士の中点からみて、中点よりも三次元変位計測で使用する視点に対応する位置に配置されたカメラの方に近い位置に配置される。すなわち、位置ｐ３（ｐ４）と視点位置ｐ１（ｐ２）との間に配置される校正用の視点は、中間位置ｃと位置ｐ３（ｐ４）との間に配置される校正用の視点よりも多い。また、他の配置例では、位置ｐ３（ｐ４）と中間位置ｃとの間に校正用の視点が設けられない。

このため、三次元変位計測方法では、カメラ校正に使用する視点を、三次元変位計測に使用する視点の位置を基準として位置ｐ３（ｐ４）よりも近い範囲に、位置ｐ３（ｐ４）よりも遠い範囲よりも多く設定する。このような不均等な配置により、カメラ校正の精度が向上する。

次に、多視点フレームセットの詳細について説明する。図１１は、実施の形態における多視点フレームセットの例を示す図である。ここでは、取得部２３０が、５台のカメラ１００－１～１００－５から１フレームずつを選択することで多視点フレームセットを決定する例を説明する。

また、複数のカメラが同期撮影することを仮定している。各フレームのヘッダ情報には、撮影されたカメラを特定するカメラＩＤがそれぞれ１００－１～１００－５として付与されている。また、各フレームのヘッダ情報には、各カメラ内での撮影順序を示すフレーム番号００１～Ｑが付与されており、カメラ間で同じフレーム番号を持つフレームは同時刻の被写体１００１が撮影されたことを示す。

取得部２３０は、多視点フレームセット２００－１～２００－Ｑを三次元画像処理部２４０へ順次出力する。三次元画像処理部２４０は、繰り返し処理により多視点フレームセット２００－１～２００－Ｑを用いて、順次三次元モデリングまたは三次元変位計測を行う。

多視点フレームセット２００－１は、カメラ１００－１のフレーム番号００１、カメラ１００－２のフレーム番号００１、カメラ１００－３のフレーム番号００１、カメラ１００－４のフレーム番号００１、カメラ１００－５のフレーム番号００１の５枚のフレームから構成される。多視点映像の最初のフレームの集合として、これら５枚のフレームが、三次元画像処理部２４０の繰り返し処理１で使用される。すなわち、繰り返し処理１では、フレーム番号００１を撮影した時刻の三次元モデルが再構成され、または、三次元変位が計測される。

図１１に示す例では、多視点フレームセット２００－２では、全てのカメラでフレーム番号が更新される。多視点フレームセット２００－２は、カメラ１００－１のフレーム番号００２、カメラ１００－２のフレーム番号００２、カメラ１００－３のフレーム番号００２、カメラ１００－４のフレーム番号００２、カメラ１００－５のフレーム番号００２の５枚のフレームから構成される。これら５枚のフレームが、三次元画像処理部２４０の繰り返し処理２で使用される。すなわち、繰り返し処理２では、フレーム番号００２を撮影した時刻の三次元モデルが再構成され、または、三次元変位が計測される。

以下、繰り返し処理３以降でも、繰り返し処理２と同様に、全てのカメラでフレーム番号が更新される。これにより、三次元再構成システム３は、各時刻における被写体１００１の三次元モデルを再構成することができ、また、各時刻における被写体１００１の三次元変位を計測することができる。

ただし、各時刻で独立して三次元再構成を行うため、再構成された複数の三次元モデルの座標軸とスケールが一致しているとは限らない。つまり、動く被写体１００１の三次元モデルを取得するためには、各時刻の座標軸及びスケールを合せる必要がある。

その場合、各フレームには撮影時刻が付与されており、その撮影時刻を基に取得部２３０は、同期フレームと非同期フレームを組み合わせた多視点フレームセットを作成する。以下、２台のカメラ間での撮影時刻を用いた同期フレームと非同期フレームの判定方法を説明する。

カメラ１００－１から選択したフレームの撮影時刻をＴ１とし、カメラ１００－２から選択したフレームの撮影時刻をＴ２とし、カメラ１００－１の露光時間をＴＥ１とし、カメラ１００－２の露光時間をＴＥ２とする。ここで、撮影時刻Ｔ１、Ｔ２は、図４Ａ及び図４Ｂの例で露光が開始された時刻、つまり矩形信号の立ち上がりの時刻を指している。

この場合、カメラ１００－１の露光終了時刻はＴ１＋ＴＥ１である。この時、式（１３）又は式（１４）が成立していれば、２台のカメラ１００－１、１００－２は、同じ時刻の被写体１００１を撮影していることになり、２枚のフレームは同期フレームと判定される。

Ｔ１≦Ｔ２≦Ｔ１＋ＴＥ１ 式（１３）
Ｔ１≦Ｔ２＋ＴＥ２≦Ｔ１＋ＴＥ１ 式（１４）

なお、制御部３００が決定した視点数及び視点位置に基づいて、カメラ１００－１～１００－ｎが多視点映像を撮影した場合、取得部２３０は、入力された多視点映像の本数と同数のフレームを多視点フレームセットとしてもよい。一方、カメラ１００－１～１００－ｎが、制御部３００が決定した視点数または視点位置と関係なく、多視点映像を撮影した場合、取得部２３０が、制御部３００が決定した視点数または視点位置に対応した映像を選択し、多視点フレームセットを作成することとしてもよい。

つまり、三次元変位計測装置２０００は、異なる時刻毎に撮像されたＮ枚の画像を１セットとする。三次元変位計測装置２０００は、複数セットのセット毎について、当該セットに含まれるＮ枚の画像とカメラパラメータとを用いて、当該セットの時刻における被写体１００１の三次元位置を計測する。

図１２は、実施の形態における三次元画像処理部のブロック図である。三次元画像処理部２４０は、カメラ校正部４１０及び三次元処理部４２０を備える。なお、三次元処理部４２０は、変位計測部の具体例である。

カメラ校正部４１０は、カメラ校正を行う。具体的には、カメラ校正部４１０は、複数のカメラ１００－１～１００－ｎのそれぞれの内部パラメータ、外部パラメータ及びレンズ歪み係数をカメラパラメータとして算出する。内部パラメータとは、カメラの焦点距離、収差、画像中心等の光学系の特性を示す。外部パラメータとは、三次元空間におけるカメラの位置及び姿勢を示す。

カメラ校正部４１０は、複数のカメラ１００－１～１００－ｎがチェッカボードの白黒の交点を撮影することにより得られたＭ枚のフレームであるＭ枚の校正用画像を用いて内部パラメータ、外部パラメータ及びレンズ歪み係数を別々に算出してもよいし、Ｓｔｒｕｃｔｕｒｅ ｆｒｏｍ ＭｏｔｉｏｎのようにＭ枚のフレーム間の対応点を用いて内部パラメータ、外部パラメータ及びレンズ歪み係数を一括して算出し、全体最適化を行ってもよい。後者の場合のＭ枚のフレームは、チェッカボードが撮像された画像でなくてもよい。

三次元処理部４２０は、多視点フレームセットとカメラパラメータとを用いて、被写体１００１の三次元モデリング、被写体１００１の三次元変位計測等の三次元処理を行う。具体的には、三次元処理部４２０は、多視点フレームセット及びカメラ校正の結果を用いて、三次元モデリングとして、被写体１００１のモデルを三次元点から生成する。三次元モデリングの代わりに、三次元処理部４２０は、三次元変位計測として、多視点フレームセット及びカメラ校正の結果を用いて、被写体１００１の三次元変位を計測してもよい。三次元処理部４２０は、三次元変位を計測する際に、三次元モデリングで生成された被写体１００１の三次元モデルの三次元点のデータを用いてもよい。ここで、三次元変位計測で利用される三次元点は、被写体１００１の任意の点でもよい。

三次元画像処理部２４０、並びに、三次元画像処理部２４０を構成する、カメラ校正部４１０及び三次元処理部４２０は、プロセッサとメモリとによって実現される。また、これらは、専用の論理回路によって実現されてもよい。

図１３は、実施の形態における三次元画像処理部の動作を示すフローチャートである。なお、図１３に示す処理では、制御部３００において決定された視点数の多視点フレームセットが用いられる。

まず、カメラ校正部４１０は、制御部３００が決定した視点数Ｍ及び視点位置に対応する多視点フレームセットを用いて、撮影が行われた各視点に存在するカメラのカメラパラメータを算出する（ステップＳ４１００）。ここで、制御部３００は、少なくともＮ台のカメラについてのみ、カメラパラメータを算出してもよい。なお、ステップＳ４１００は、カメラ校正ステップの一例である。

次に、三次元処理部４２０は、制御部３００が決定した視点数Ｎ及び視点位置に対応する多視点フレームセット、及び、カメラ校正部４１０が算出したカメラパラメータを用いて、被写体１００１の三次元モデルを生成、または、被写体１００１の三次元変位を計測する（ステップＳ４２００）。ここで、三次元処理部４２０は、変位計測部の具体例である。なお、ステップＳ４２００は、変位計測ステップの一例である。

なお、本実施の形態で、基準として用いる視点数は、三次元モデリング又は三次元変位計測に用いられる視点数Ｎである。例えば、視点数Ｎは、撮影空間１０００に常時設置されているカメラの数に等しい。

また、ステップＳ４１００で算出されたカメラパラメータの精度は、ステップＳ４２００で行われる被写体１００１の三次元モデルの生成または被写体１００１の三次元変位の計測の精度に大きな影響を及ぼす。よって、カメラパラメータの精度を高めるため、ステップＳ４１００では、ステップＳ４２００の視点数Ｎより多い視点数Ｍを使用する。このとき、カメラの台数は必ずしもＭ台でなくてもよい。Ｍ台未満の移動カメラによって、複数視点での撮影を行い、撮影された画像内の静止領域が使用されて、カメラ校正が行われてもよい。静止領域は、画像中の、静止物体が映っている領域であり、動物体が映っている動領域を除く領域である。

なお、三次元変位計測装置２０００は、上述しているように三次元画像処理部２４０において、三次元変位計測の他に、三次元モデリングを行ってもよく、さらに自由視点映像生成を行ってもよい。

以下、実施の形態における三次元変位計測装置２０００で行われる自由視点映像生成について説明する。

三次元変位計測装置２０００は、時刻間で座標軸の一致した時系列三次元モデルを再構成できる。具体的には、まず、三次元変位計測装置２０００は、時刻毎に独立して三次元再構成を行うことで、各時刻の三次元モデルを取得する。次に、三次元変位計測装置２０００は、静止カメラ及び静止物体（静止三次元点）を検出し、検出した静止カメラ及び静止物体を用いて、時刻間で三次元モデルの座標合せを行い、座標軸の一致した時系列三次元モデルを生成する。

次に、三次元変位計測装置２０００が自由視点映像生成を行う際の動作を説明する。図１４は、本実施の形態に係る三次元変位計測装置２０００が自由視点映像生成を行う際の動作を示すフローチャートである。

まず、受信部２１０は、複数のカメラ１００－１～１００－ｎで撮影された多視点映像を受信する（Ｓ１０１）。記憶部２２０は、受信された多視点映像を記憶する（Ｓ１０２）。

次に、取得部２３０は、多視点映像からフレームを選択し、多視点フレームセットとして三次元画像処理部２４０へ出力する（Ｓ１０３）。なお、ここで選択される多視点フレームセットは、図６のステップＳ２５００で選択される多視点フレームセットと同様の構成を有していてもよい。

次に、三次元画像処理部２４０は、多視点フレームセットを用いて、カメラ校正処理、三次元モデリング処理及び自由視点映像合成処理を実行することで自由視点映像を生成する（Ｓ１０４）。

また、ステップＳ１０３及びＳ１０４の処理は、多視点フレームセット毎に繰り返し行われる。

最後に、送信部２５０は、カメラパラメータ、被写体の三次元モデル及び自由視点映像の少なくとも１つを外部装置へ送信する（Ｓ１０５）。

次に、三次元画像処理部２４０による自由視点映像生成の詳細について説明する。

制御部３００は、図７～図１０を用いて説明したように、カメラ校正部４１０及び三次元処理部４２０における各処理で最適な視点数を決定する。ここで決定する視点数とは、互いに異なる視点の数を示す。

カメラ校正処理において算出するカメラパラメータの精度は、三次元モデリング処理及び自由視点映像合成処理における精度に大きな影響を及ぼす。よって、制御部３００は、三次元モデリング処理及び自由視点映像合成処理における精度を低下させないために、カメラパラメータの精度が向上するように三次元モデリング処理における視点数Ｎよりも多い視点数Ｍを、カメラ校正処理において用いる多視点フレームセットの視点数として決定する。つまり、制御部３００は、ｎ台のカメラ１００－１～１００－ｎにより撮像されたＮ枚のフレームに、ａ台の非固定カメラにより撮像されたｋ（ｋはａ以上の整数）枚のフレームを加えたＭ枚のフレームを用いてカメラ校正部４１０にカメラ校正処理を実行させる。なお、ａ台の非固定カメラは、必ずしもｋ台でなくてもよく、ａ台の非固定カメラを移動させることによりｋ視点で撮像を行った結果得られた、ｋ枚のフレーム（画像）であってもよい。

また、自由視点映像合成処理において、実カメラによって得られた画像と、仮想視点の画像との対応位置の算出には、実カメラ台数が多いほど大きい処理負荷がかかるため、多くの処理時間を要する。一方で、ｎ台のカメラ１００－１～１００－ｎのうち配置されている位置が近い複数のカメラにおいて得られた複数の画像間において、当該複数の画像から得られるテクスチャ情報が互いに似ている。このため、自由視点映像合成処理に、当該複数の画像の全てを用いても、当該複数の画像のうちの１つの画像を用いても、自由視点映像合成処理により得られる結果への精度はあまり変わらない。よって、制御部３００は、三次元モデリング処理における視点数Ｎよりも少ない視点数ｌを、自由視点映像合成処理において用いる多視点フレームセットの視点数として決定する。

図１５は、自由視点映像生成を行う際の三次元画像処理部２４０の動作を示すフローチャートである。

まず、カメラ校正部４１０は、ステップＳ４１００と同様に、互いに異なる位置に配置されるｎ台のカメラ１００－１～１００－ｎによって異なるＭ視点において撮像されたＭ枚の校正用画像を用いて複数のカメラ１００－１～１００－ｎのカメラパラメータを算出する（Ｓ３１０）。なお、ここでのＭ視点は、制御部３００において決定された視点数に基づく。

なお、カメラ校正部４１０は、固定カメラ及び非固定カメラを含むｎ台のカメラ１００－１～１００－ｎによって得られたＭ枚の校正用画像を用いてカメラ校正処理を行う。カメラ校正処理では、カメラの数が多いほどカメラ間の距離が近くなり、距離が近い複数のカメラの視野が近くなるため、距離が近い複数のカメラから得られる複数の画像の対応付けが容易になる。よって、カメラ校正部４１０は、カメラ校正を行う場合、撮影空間１０００に常時設置されている固定カメラに加えて、非固定カメラを用いて視点数を増やす。

非固定カメラは、少なくとも１台の移動カメラでもよく、非固定カメラとして移動カメラを使用する場合、撮像するタイミングが異なる画像が含まれることとなる。つまり、カメラ校正処理において用いられるＭ枚の校正用画像は、異なるタイミングで撮像された画像を含むことになる。言い換えると、Ｍ枚の校正用画像が構成するＭ視点の多視点フレームセットは、非同期撮影により得られたフレームを含む。このため、カメラ校正部４１０は、Ｍ枚の校正用画像のうちの静止物体が映っている領域である静止領域から得られる特徴点の画像間同士の対応点を利用してカメラ校正処理を行う。よって、カメラ校正部４１０は、静止領域に対応したカメラパラメータを算出する。静止領域は、Ｍ枚の校正用画像のうちの動物体が映っている動領域を除く領域である。フレームに映り込む動領域は、例えば、過去のフレームとの差分を計算したり、背景映像との差分を計算したり、機械学習により動物体の領域を自動検知するなどで検出される。

なお、カメラ校正部４１０は、三次元画像処理部２４０における自由視点映像生成処理において、ステップＳ３１０のカメラ校正処理を常に行わなくてもよく、所定の回数毎に１回行ってもよい。

次に、三次元処理部４２０は、ｎ台のカメラ１００－１～１００－ｎのそれぞれによって撮像されたＮ枚の三次元処理用画像、及び、カメラ校正処理において得られたカメラパラメータ、を用いて三次元モデルを再構成する（Ｓ３２０）。つまり、三次元処理部４２０は、制御部３００において決定された視点数Ｎに基づいて、Ｎ視点において撮像されたＮ枚の三次元処理用画像を用いて三次元モデルを再構成する。これにより、三次元処理部４２０は、Ｎ枚の三次元処理用画像における被写体を三次元点として再構成する。三次元モデリング処理において用いられるＮ枚の三次元処理用画像は、任意のタイミングでｎ台のカメラ１００－１～１００－ｎのそれぞれによって撮像された画像である。つまり、Ｎ枚の三次元処理用画像が構成するＮ視点の多視点フレームセットは、同期撮影により得られた多視点フレームセットである。このため、三次元処理部４２０は、Ｎ枚の三次元処理用画像のうち静止物体及び動物体を含む領域（つまり、全ての領域）を用いて三次元モデリング処理を行う。なお、三次元処理部４２０は、レーザスキャンを用いて被写体の三次元空間上の位置の計測結果を用いてもよいし、多視点ステレオ法のように複数のステレオ画像の対応点を用いて、被写体の三次元空間上の位置を算出してもよい。

次に、三次元処理部４２０は、ｎ台のカメラ１００－１～１００－ｎのうちの、ｌ台のカメラのそれぞれによって撮像されたｌ枚の第３画像、カメラ校正処理において算出されたカメラパラメータ、及び、三次元モデリング処理において再構成された三次元モデル、を用いて自由視点映像を合成する（Ｓ３３０）。つまり、三次元処理部４２０は、制御部３００において決定された視点数ｌに基づいて、ｌ視点において撮像されたｌ枚の第３画像を用いて自由視点映像を合成する。具体的には、三次元処理部４２０は、カメラパラメータ及び三次元モデルにより求めた、実カメラの画像と仮想視点の画像との対応位置を基に、実カメラのテクスチャ情報を用いて仮想視点のテクスチャ情報を算出することで、自由視点映像を合成する。

（効果など）
本実施の形態に係る三次元変位計測方法によれば、カメラパラメータの精度が向上するように三次元変位計測処理における視点数Ｎよりも多い視点数Ｍを、カメラ校正処理に用いることで、三次元変位計測処理における精度を向上させることができる。

また、本実施の形態に係る三次元変位計測方法によれば、精度が高いカメラパラメータ、第１対応点及び第２対応点に基づいて、画像マッチングの手法及び三角測量の原理を用いて、三次元点の経時的な変位を求める。このため、被写体の三次元空間における変位を効果的に求めることができる。

また、本実施の形態に係る三次元変位計測方法によれば、変位計測ステップで行われる計測の精度を満たすように、カメラの視点数Ｎを定めることができる。よって、三次元変位計測方法では、精度よく被写体の三次元空間における変位を求めることができる。

本実施の形態に係る三次元変位計測装置２０００によれば、カメラ校正処理において算出するカメラパラメータの精度が、三次元モデリング処理及び自由視点映像合成処理における精度に大きな影響を及ぼすことを考慮して、カメラパラメータの精度が向上するように三次元モデリング処理における視点数Ｎよりも多い視点数Ｍを、カメラ校正処理において用いる多視点フレームセットの視点数として決定する。このため、三次元モデリング処理及び自由視点映像合成処理における精度を向上させることができる。

また、本実施の形態に係る三次元変位計測装置２０００によれば、三次元モデリング処理における視点数Ｎよりも少ない視点数ｌを、自由視点映像合成処理において用いる多視点フレームセットの視点数として決定することで、自由視点映像を生成するのに要する処理負荷を低減することができる。

（変形例）
図１６は、三次元画像処理部の変形例のブロック図である。本実施の形態における変形例では、三次元画像処理部２４０Ａは、カメラ校正部５１０及び三次元処理部５２０で構成されている。

カメラ校正部５１０は、カメラ校正部４１０と機能、構成等は同一であるが、以下の点において異なる。カメラ校正部５１０は、カメラ校正部５１０の出力をフィードバックされた入力として受け付け、その入力に基づいて再度処理を行う。三次元処理部５２０は、三次元処理部４２０と機能及び構成は同一である。三次元処理部５２０は、カメラ校正部５１０から最終的に出力されたデータに基づいて、処理を行う。なお、カメラ校正部からの最終的な出力とは、例えば、カメラ校正部からの最初の出力を除いた出力である。

図１７は、三次元画像処理部の変形例の動作を示すフローチャートである。

まず、カメラ校正部５１０は、制御部３００が決定した視点数及び視点位置に対応する多視点フレームセットを用いて、撮影が行われた各視点に存在するカメラのカメラパラメータを算出する（ステップＳ５１００）。

次に、カメラ校正部５１０は、ステップＳ５１００で行われたカメラ校正の結果の値を初期値として、制御部３００が決定した視点数及び視点位置に対応する多視点フレームセットを用いて、撮影が行われた各視点に存在するカメラのカメラパラメータを算出する（ステップＳ５２００）。

そして、三次元処理部５２０は、制御部３００が決定した視点数及び視点位置に対応する多視点フレームセット、及び、カメラ校正部５１０が算出したカメラパラメータを用いて、被写体１００１の三次元モデリング、または、被写体の三次元変位計測を行う（ステップＳ５３００）。

カメラ校正処理において、カメラ校正に使用されるカメラの台数が多いほど、カメラ間の距離が近くなる。距離が近いカメラ間では、カメラの視野に入る撮影空間１０００の領域同士が近くなるため、カメラごとに撮影された画像同士の対応付けが容易になる。よって、カメラ校正を行う場合、撮影空間１０００に常時設置されている固定カメラに加えて、非固定カメラを用いて、カメラ校正に用いられる視点の数を、固定カメラで実現される視点の数よりも増加させる。ここで、非固定カメラは、少なくとも１台の移動カメラでもよい。しかしながら、非固定カメラとして移動カメラを使用する場合、異なる視点において撮影された画像の撮影時刻が異なるため、撮影された画像の静止領域のみ対応したカメラパラメータが算出される。

このとき、算出されたカメラパラメータは、固定カメラで撮影された画像における動領域に対応しているとは限らない。ここで、動領域とは、画像中の、動物体が映っている領域である。また、Ｓｔｒｕｃｔｕｒｅ ｆｒｏｍ Ｍｏｔｉｏｎのような方式では、カメラパラメータの全体最適化が実施される。よって、この方式が用いられた場合は、固定カメラのみに着目すると、カメラパラメータが最適化されているとは限らない。

よって、本変形例では、ステップＳ５１００とステップＳ５２００との２つのステップにわたって、２段階のカメラ校正を実施する。

まず、ステップＳ５１００では、カメラ校正部５１０は、撮影空間１０００に常時設置されているＡ台のカメラ以外の少なくとも１台のカメラで撮影された画像を含んだ、Ｂ個の視点が用いられる多視点フレームセットを用いて、大まかなカメラ校正を行う。

次に、ステップＳ５２００では、カメラ校正部５１０は、撮影空間１０００に常時設置されているＡ台のカメラで撮影された画像のみを用いて、カメラ校正部５１０が算出したカメラパラメータをＡ台のカメラ環境で最適化する。ここで、最適化とは、カメラパラメータ算出の際に副次的に得られた三次元点を、撮影された画像上に再投影し、その再投影点と画像上で検出された特徴点との誤差（再投影誤差）を評価値として、評価値を最小化する処理を行うことを指してもよい。

なお、２段階目のカメラ校正で用いる視点数は、撮影空間１０００に常時設置されるカメラ台数Ａより大きくてもよい。また、２段階目のカメラ校正で用いられる視点数は、三次元画像処理（例えば三次元変位計測）で用いられる視点数よりも大きくてもよい。ただし、２段階目のカメラ校正で用いられる視点数は、三次元画像処理（例えば三次元変位計測）で用いられる視点数と等しくてもよいし、それよりも小さくてもよい。

上記の三次元変位計測装置において、ｎ台のカメラ１００－１～１００－ｎは、２つのカメラを有するステレオカメラを含んでいてもよい。ステレオカメラは、互いに略同じ方向を撮像する２つのカメラ、つまり、第１カメラ及び第２カメラを有し、２つのカメラの間の距離が所定距離以下である構成であればよい。このように、ｎ台のカメラ１００－１～１００－ｎがステレオカメラにより構成された場合、ｎ／２台の第１カメラと、ｎ／２台の第２カメラとにより構成される。なお、ステレオカメラが有する２つのカメラは、一体化されていてもよいし、別体であってもよい。

また、ステレオカメラを構成する第１カメラ及び第２カメラは、互いに異なる感度で撮像してもよい。第１カメラは、第１の感度で撮像するカメラである。第２カメラは、第１の感度とは異なる第２の感度で撮像するカメラである。第１カメラと第２カメラとは、色感度が互いに異なるカメラである。

三次元処理部は、ｎ台のカメラ１００－１～１００－ｎの全てによって撮像されることにより得られたｎ枚の三次元処理用画像を用いて三次元モデルを再構成する。三次元処理部は、三次元モデリング処理において、輝度情報を使用するため、色感度の相違に関わらずｎ台のカメラ全てを使用して三次元モデルを高精度に算出することができる。

三次元処理部は、ｎ／２台の第１カメラまたはｎ／２台の第２カメラによって撮像されることにより得られた複数の画像であるｎ／２枚の第３画像、カメラ校正部により算出されたカメラパラメータ、及び、三次元処理部により再構成された三次元モデル、を用いて自由視点映像を合成する。三次元処理部は、自由視点映像生成処理において、ｎ／２台の第１カメラ、及び、ｎ／２台の第２カメラのどちらか一方によるｎ／２枚の画像を使用しても精度に及ぼす影響は小さい。そこで、三次元処理部は、撮影空間１０００の状況に応じて、第１カメラと第２カメラとの一方で撮像されたｎ／２枚の画像を用いて、自由視点合成を実施する。例えば、ｎ／２台の第１カメラは、赤系統の色感度が高いカメラであり、ｎ／２台の第２カメラは、青系統の色感度が高いカメラであるとする。この場合、三次元処理部は、被写体が赤系統の色であれば、赤の色感度が高い第１カメラにより撮像された画像を用い、被写体が青系統の色であれば、青の色感度が高い第２カメラにより撮像された画像を用いて自由視点映像合成処理を実行するように、用いる画像を切り替える。

三次元変位計測装置によれば、撮影空間の状況に応じて感度の異なる２種類のカメラから得られる２種類の画像のうちの一方の画像を用いて自由視点映像合成を行う。このため、精度よく自由視点映像を生成することができる。

なお、第１カメラと第２カメラとは、色感度が互いに異なることに限らず、輝度感度が互いに異なるカメラであってもよい。この場合、三次元処理部は、昼間と夜間、晴天と曇天などの状況に応じて、カメラを切替えることができる。

なお、ステレオカメラを用いるとしたが、必ずしもステレオカメラを用いなくてもよい。よって、ｎ台のカメラは、ｎ／２台の第１カメラ、及び、ｎ／２台の第２カメラにより構成されていることに限らずに、ｉ台の第１カメラと、ｊ台の第２カメラとにより構成されていてもよい。

複数のカメラ１００－１～１００－ｎは、固定カメラ及び非固定カメラにより構成されるとしたが、これに限らずに、全ての複数のカメラが固定カメラにより構成されてもよい。また、三次元モデリングで用いられるｎ枚の画像は、固定カメラにより撮像された画像であるとしたが、非固定カメラにより撮像された画像を含んでいてもよい。

図１８は、三次元空間認識の概要を示す図である。上述のように、校正済みのカメラ（例えば固定カメラ）を用いて同一空間を多視点から撮影することにより撮影する空間を三次元再構成できる（三次元空間再構成）。この三次元再構成されたデータを用いて、三次元変位計測等を行うことができる。これにより、次世代広域監視システム、及び三次元変位計測システムを実現できる。具体的には、車道における交差点等において、通行する車両１０等を固定カメラ１００－ａまたは非固定カメラ１００－ｂを用いて撮影し、車両１０等の三次元変位を測定することが出来る。

上記の三次元変位計測装置によれば、被写体の三次元変位計測を行い、さらに、被写体の三次元モデリングを行ってもよいとしたが、これに限らずに、被写体の三次元変位計測を行わずに被写体の三次元モデリングを行ってもよい。

本開示は、三次元変位計測方法及び三次元変位計測装置に適用でき、例えば、インフラ点検システム、産業用機械モニタリングシステム、三次元空間認識システム及び次世代監視システム等に適用できる。

３ 三次元再構成システム
１０ 車両
１００－１～１００－ｎ カメラ
１００－ａ 固定カメラ
１００－ｂ 非固定カメラ
２００－１～２００－Ｑ 多視点フレームセット
２１０ 受信部
２２０ 記憶部
２３０ 取得部
２４０、２４０Ａ 三次元画像処理部
２５０ 送信部
３００ 制御部
３１０ 計測用決定部
３２０ カメラ校正用決定部
４１０、５１０ カメラ校正部
４２０、５２０ 三次元処理部
１０００ 撮影空間
１００１ 被写体
２０００ 三次元変位計測装置

Claims (9)

1. Ｎ（Ｎは１より大きい整数）台のカメラを含む複数のカメラによって互いに異なるＭ（ＭはＮより大きい整数）個の視点において撮像されたＭ枚の校正用画像に基づいて、前記Ｎ台のカメラのカメラパラメータを算出するカメラ校正ステップと、
   （１）異なる時刻毎に前記Ｎ台のカメラのそれぞれによって被写体を撮像することで生成されるＮ枚の画像、及び、（２）前記カメラパラメータ、に基づいて、三次元空間における前記被写体の振動による前記被写体の経時的な変位を計測する変位計測ステップと、を含み、
   前記異なる時刻は、前記被写体の振動の周波数に応じて設定される
   三次元変位計測方法。
2. 前記振動の周波数は、前記被写体の共振周波数である
   請求項１に記載の三次元変位計測方法。
3. 前記異なる時刻は、前記共振周波数の逆数で表される時間毎における３以上の異なる時刻を含む
   請求項２に記載の三次元変位計測方法。
4. 前記変位計測ステップは、
   前記Ｎ枚の画像に含まれる第１画像に含まれる、前記被写体の第１対応点と、前記Ｎ枚の画像に含まれる第２画像に含まれ且つ前記第１対応点に対応する第２対応点と、を求める対応点算出ステップと、
   前記カメラパラメータと、前記第１対応点と、前記第２対応点とに基づいて、前記三次元空間における前記被写体の三次元点を求める三次元点算出ステップと、
   前記異なる時刻に対応する複数の前記三次元点から、前記三次元点の経時的な変位を求める変位算出ステップと、を含む、
   請求項１～３のいずれか１項に記載の三次元変位計測方法。
5. 前記Ｎは、前記カメラパラメータと、前記変位の計測の精度と、に基づいて決定される、
   請求項１～４のいずれか１項に記載の三次元変位計測方法。
6. 前記Ｎ台のカメラは、互いに隣り合う第１カメラ及び第２カメラを含み、
   前記Ｍ個の視点は、前記第１カメラが配置される第１位置と、前記第２カメラが配置される第２位置との間の第１中点を含み、
   前記第１位置と前記第１中点との間の中点を第２中点とすると、前記Ｍ個の視点に含まれ且つ前記第２中点と前記第１位置との間に設けられる視点は、前記Ｍ個の視点に含まれ且つ前記第１中点と前記第２中点との間に設けられる視点よりも多い
   請求項１～５のいずれか１項に記載の三次元変位計測方法。
7. 前記カメラ校正ステップでは、前記第１中点と前記第２中点との間には、視点を設けない
   請求項６に記載の三次元変位計測方法。
8. 前記Ｎ台のカメラは、非固定カメラを含む
   請求項１～７のいずれか１項に記載の三次元変位計測方法。
9. Ｎ（Ｎは１より大きい整数）台のカメラを含む複数のカメラによって互いに異なるＭ（ＭはＮより大きい整数）個の視点において撮像されたＭ枚の校正用画像に基づいて、前記Ｎ台のカメラのカメラパラメータを算出するカメラ校正部と、
   （１）異なる時刻毎に前記Ｎ台のカメラのそれぞれによって被写体を撮像することで生成されるＮ枚の画像、及び、（２）前記カメラパラメータ、に基づいて、三次元空間における前記被写体の振動による前記被写体の経時的な変位を計測する変位計測部と、を備え、
   前記異なる時刻は、前記被写体の振動の周波数に応じて設定される
   三次元変位計測装置。

Images