(57)【要約】【課題】特殊な環境下でも、計測対象を高精度で効率良く計測することができる３次元位置計測システム及び方法を提供する。【解決手段】互いに異なる方向から撮影する第１及び第２カメラ１１Ａ、１１Ｂと、第１及び第２カメラ１１Ａ、１１Ｂで撮影した画像の処理により、計測対象点の３次元位置を計測するデータ処理装置２０とを有し、データ処理装置２０は、基準格子画像と計測画像を取得する画像取得部２１と、基準格子画像と計測画像とを用いて算出した複数の直線について、直線同士を最短で結ぶ線分の距離に基づいて、計測対象点の算出に最適な２つの直線を設定し、設定した２つの直線同士を最短で結ぶ線分を求め、当該線分の中点の３次元座標を算出し、当該３次元座標を計測対象点の３次元位置とする画像解析部２２Ａとを有する。【選択図】図１
2018146363.tif 000002
【特許請求の範囲】
【請求項１】
対象物を異なる方向から撮影する複数の撮影手段と、
前記複数の撮影手段で撮影した画像を取得し、前記画像の処理により、計測対象点の３次元位置を計測する処理手段とを有し、
前記処理手段は、
前記計測の前に、前記複数の撮影手段で３次元座標が既知である格子の画像を各々撮影し、撮影した前記画像を各々複数の基準格子画像として取得する計測前画像取得部と、
前記計測時に、前記複数の撮影手段で前記計測対象点を各々撮影し、撮影した前記画像を各々複数の計測画像として取得する計測時画像取得部と、
前記複数の基準格子画像の１つと、それに対応する前記複数の計測画像の１つにより、前記撮影手段から前記計測対象点への視線となる線を算出する線算出部と、
前記複数の線から前記計測対象点の３次元位置を算出する３次元位置算出部とを有する
ことを特徴とする３次元位置計測システム。
【請求項２】
対象物を互いに異なる方向から撮影する第１及び第２の撮影手段と、
前記第１及び第２の撮影手段で撮影した画像を取得し、前記画像の処理により、計測対象点の３次元位置を計測する処理手段とを有し、
前記処理手段は、
前記計測の前に、目盛となる基準格子を表面に有する格子定規を、３次元座標が既知である複数の設置位置に設置しながら、一端側から他端側へ順に１〜ｎ番目とする前記設置位置の各々において、前記第１及び第２の撮影手段で前記格子定規を各々撮影し、撮影した前記画像を各々第１及び第２の基準格子画像として取得する計測前画像取得部と、
前記計測時に、前記第１及び第２の撮影手段で前記計測対象点を各々撮影し、撮影した前記画像を各々第１及び第２の計測画像として取得する計測時画像取得部と、
ｉ番目とｊ番目の前記第１の基準格子画像の各々に前記第１の計測画像を重ね合わせ、各々の前記基準格子画像中の前記基準格子に対する前記第１の計測画像中の前記計測対象点の位置に基づいて、各々の前記基準格子上に前記計測対象点があると仮定した場合の前記計測対象点の３次元座標を各々ｉ番目とｊ番目の第１の対象点座標として算出すると共に、ｉ番目の前記第１の対象点座標とｊ番目の前記第１の対象点座標とを通る第１の直線を算出する第１直線算出部と、
ｉ番目とｊ番目の前記第２の基準格子画像の各々に前記第２の計測画像を重ね合わせ、各々の前記基準格子画像中の前記基準格子に対する前記第２の計測画像中の前記計測対象点の位置に基づいて、各々の前記基準格子上に前記計測対象点があると仮定した場合の前記計測対象点の３次元座標を各々ｉ番目とｊ番目の第２の対象点座標として算出すると共に、ｉ番目の前記第２の対象点座標とｊ番目の前記第２の対象点座標とを通る第２の直線を算出する第２直線算出部と、
前記第１の直線と前記第２の直線とを最短で結ぶ線分の距離又は前記線分の中点の３次元座標に基づいて、前記計測対象点の算出に最適な前記第１の直線と前記第２の直線を設定する直線設定部と、
最適な前記第１の直線と最適な前記第２の直線とを最短で結ぶ線分を求め、当該線分の中点の３次元座標を算出し、当該３次元座標を前記計測対象点の３次元位置とする３次元位置算出部とを有する
ことを特徴とする３次元位置計測システム。
但し、ｎは２以上の自然数、ｉ＝１〜ｎの自然数、ｊ＝１〜ｎの自然数、ｉ≠ｊである。
【請求項３】
請求項２に記載の３次元位置計測システムにおいて、
前記処理手段は、
前記第１の直線と前記第２の直線とを最短で結ぶ線分の距離を算出する距離算出部と、
複数の前記距離の中から最小値を算出する最小値算出部とを有し、
前記直線設定部は、
前記第１直線算出部と前記第２直線算出部とを用いて、１〜ｎ−１の自然数であるｋの全てにおいて、ｉ＝ｋ及びｊ＝ｋ＋１として、ｋ番目の前記第１の直線とｋ番目の前記第２の直線とを各々算出し、
前記距離算出部を用いて、１〜ｎ−１の自然数であるｋの全てにおいて、ｋ番目の前記第１の直線とｋ番目の前記第２の直線とを最短で結ぶ線分のｋ番目の前記距離を各々算出し、
前記最小値算出部を用いて、１〜ｎ−１番目の前記距離の中から最小値となるｍ番目の前記距離を算出し、
最適な前記第１の直線と最適な前記第２の直線として、ｍ番目の前記第１の直線とｍ番目の前記第２の直線を設定する
ことを特徴とする３次元位置計測システム。
但し、ｍ＝１〜ｎ−１の自然数である。
【請求項４】
請求項２に記載の３次元位置計測システムにおいて、
前記処理手段は、
前記第１の直線と前記第２の直線とを最短で結ぶ線分を求め、当該線分の中点の３次元座標を算出する中点座標算出部を有し、
前記直線設定部は、
前記第１直線算出部と前記第２直線算出部とを用いて、ｉ＝１及びｊ＝ｎとして、初期の前記第１の直線と初期の前記第２の直線とを各々算出し、
前記中点座標算出部を用いて、初期の前記第１の直線と初期の前記第２の直線とを最短で結ぶ初期の前記線分を求め、初期の前記線分の中点の３次元座標を算出し、
前記計測対象点の推定位置となる初期の前記線分の中点の３次元座標に基づいて、前記計測対象点を挟んで前記計測対象点に最も近い２つの前記設置位置で撮影したｍ番目とｍ＋１番目の前記第１及び第２の基準格子画像を選択し、
前記第１直線算出部と前記第２直線算出部とを用いて、ｉ＝ｍ及びｊ＝ｍ＋１として、最適な前記第１の直線と最適な前記第２の直線とを各々算出する
ことを特徴とする３次元位置計測システム。
【請求項５】
請求項１から請求項４のいずれか１つに記載の３次元位置計測システムにおいて、
前記対象物に所定のパターン映像を投影する投影手段を有し、
前記処理手段は、前記パターン映像に基づいて、同一の前記計測対象点を前記第１の計測画像中と前記第２の計測画像中で特定する対象点特定部を有する
ことを特徴とする３次元位置計測システム。
【請求項６】
請求項２から請求項５のいずれか１つに記載の３次元位置計測システムにおいて、
前記格子定規は、前記基準格子を表裏両方の表面に有し、
前記第１の撮影手段は前記格子定規の一方の表面を撮影する位置に設置され、前記第２の撮影手段は前記格子定規の他方の表面を撮影する位置に設置されている
ことを特徴とする３次元位置計測システム。
【請求項７】
請求項１から請求項６のいずれか１つに記載の３次元位置計測システムにおいて、
前記対象物は、航空機の翼の風洞模型であり、
前記計測対象点は、前記風洞模型に対する着氷風洞試験において、前記風洞模型の表面に付着する氷の頂点である
ことを特徴とする３次元位置計測システム。
【請求項８】
複数の撮影手段を用いて、対象物を異なる方向から撮影し、撮影した画像の処理により、計測対象点の３次元位置を計測する３次元位置計測方法において、
前記計測の前に、前記複数の撮影手段で３次元座標が既知である格子の画像を各々撮影し、撮影した前記画像を各々複数の基準格子画像として取得する計測前画像取得工程と、
前記計測時に、前記複数の撮影手段で前記計測対象点を各々撮影し、撮影した前記画像を各々複数の計測画像として取得する計測時画像取得工程と、
前記複数の基準格子画像の１つと、それに対応する前記複数の計測画像の１つにより、前記撮影手段から前記計測対象点への視線となる線を算出する線算出工程と、
前記複数の線から前記計測対象点の３次元位置を算出する３次元位置算出工程とを有する
ことを特徴とする３次元位置計測方法。
【請求項９】
第１及び第２の撮影手段を用いて、対象物を互いに異なる方向から撮影し、撮影した画像の処理により、計測対象点の３次元位置を計測する３次元位置計測方法において、
前記計測の前に、目盛となる基準格子を表面に有する格子定規を、３次元座標が既知である複数の設置位置に設置しながら、一端側から他端側へ順に１〜ｎ番目とする前記設置位置の各々において、前記第１及び第２の撮影手段で前記格子定規を各々撮影し、撮影した前記画像を各々第１及び第２の基準格子画像として取得する計測前画像取得工程と、
前記計測時に、前記第１及び第２の撮影手段で前記計測対象点を各々撮影し、撮影した前記画像を各々第１及び第２の計測画像として取得する計測時画像取得工程と、
ｉ番目とｊ番目の前記第１の基準格子画像の各々に前記第１の計測画像を重ね合わせ、各々の前記基準格子画像中の前記基準格子に対する前記第１の計測画像中の前記計測対象点の位置に基づいて、各々の前記基準格子上に前記計測対象点があると仮定した場合の前記計測対象点の３次元座標を各々ｉ番目とｊ番目の第１の対象点座標として算出すると共に、ｉ番目の前記第１の対象点座標とｊ番目の前記第１の対象点座標とを通る第１の直線を算出する第１直線算出工程と、
ｉ番目とｊ番目の前記第２の基準格子画像の各々に前記第２の計測画像を重ね合わせ、各々の前記基準格子画像中の前記基準格子に対する前記第２の計測画像中の前記計測対象点の位置に基づいて、各々の前記基準格子上に前記計測対象点があると仮定した場合の前記計測対象点の３次元座標を各々ｉ番目とｊ番目の第２の対象点座標として算出すると共に、ｉ番目の前記第２の対象点座標とｊ番目の前記第２の対象点座標とを通る第２の直線を算出する第２直線算出工程と、
前記第１の直線と前記第２の直線とを最短で結ぶ線分の距離又は前記線分の中点の３次元座標に基づいて、前記計測対象点の算出に最適な前記第１の直線と前記第２の直線を設定する直線設定工程と、
最適な前記第１の直線と最適な前記第２の直線とを最短で結ぶ線分を求め、当該線分の中点の３次元座標を算出し、当該３次元座標を前記計測対象点の３次元位置とする３次元位置算出工程とを有する
ことを特徴とする３次元位置計測方法。
但し、ｎは２以上の自然数、ｉ＝１〜ｎの自然数、ｊ＝１〜ｎの自然数、ｉ≠ｊである。
【請求項１０】
請求項９に記載の３次元位置計測方法において、
前記第１の直線と前記第２の直線とを最短で結ぶ線分の距離を算出する距離算出工程と、
複数の前記距離の中から最小値を算出する最小値算出工程とを有し、
前記直線設定工程は、
前記第１直線算出工程と前記第２直線算出工程とを用いて、１〜ｎ−１の自然数であるｋの全てにおいて、ｉ＝ｋ及びｊ＝ｋ＋１として、ｋ番目の前記第１の直線とｋ番目の前記第２の直線とを各々算出し、
前記距離算出工程を用いて、１〜ｎ−１の自然数であるｋの全てにおいて、ｋ番目の前記第１の直線とｋ番目の前記第２の直線とを最短で結ぶ線分のｋ番目の前記距離を各々算出し、
前記最小値算出工程を用いて、１〜ｎ−１番目の前記距離の中から最小値となるｍ番目の前記距離を算出し、
最適な前記第１の直線と最適な前記第２の直線として、ｍ番目の前記第１の直線とｍ番目の前記第２の直線を設定する
ことを特徴とする３次元位置計測方法。
但し、ｍ＝１〜ｎ−１の自然数である。
【請求項１１】
請求項９に記載の３次元位置計測方法において、
前記第１の直線と前記第２の直線とを最短で結ぶ線分を求め、当該線分の中点の３次元座標を算出する中点座標算出工程を有し、
前記直線設定工程は、
前記第１直線算出工程と前記第２直線算出工程とを用い、ｉ＝１及びｊ＝ｎとして、初期の前記第１の直線と初期の前記第２の直線とを各々算出し、
前記中点座標算出工程を用いて、初期の前記第１の直線と初期の前記第２の直線とを最短で結ぶ初期の前記線分を求め、初期の前記線分の中点の３次元座標を算出し、
前記計測対象点の推定位置となる初期の前記線分の中点の３次元座標に基づいて、前記計測対象点を挟んで前記計測対象点に最も近い２つの前記設置位置で撮影したｍ番目とｍ＋１番目の前記第１及び第２の基準格子画像を選択し、
前記第１直線算出工程と前記第２直線算出工程とを用いて、ｉ＝ｍ及びｊ＝ｍ＋１として、最適な前記第１の直線と最適な前記第２の直線とを各々算出する
ことを特徴とする３次元位置計測方法。
【請求項１２】
請求項８から請求項１１のいずれか１つに記載の３次元位置計測方法において、
前記対象物に所定のパターン映像を投影する投影手段を用い、前記パターン映像に基づいて、同一の前記計測対象点を前記第１の計測画像中と前記第２の計測画像中で特定する対象点特定工程を有する
ことを特徴とする３次元位置計測方法。
【請求項１３】
請求項９から請求項１２のいずれか１つに記載の３次元位置計測方法において、
前記格子定規は、前記基準格子を表裏両方の表面に有し、
前記第１の撮影手段は前記格子定規の一方の表面を撮影する位置に設置され、前記第２の撮影手段は前記格子定規の他方の表面を撮影する位置に設置されている
ことを特徴とする３次元位置計測方法。
【請求項１４】
請求項８から請求項１３のいずれか１つに記載の３次元位置計測方法において、
前記対象物は、航空機の翼の風洞模型であり、
前記計測対象点は、前記風洞模型に対する着氷風洞試験において、前記風洞模型の表面に付着する氷の頂点である
ことを特徴とする３次元位置計測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、計測対象の３次元位置を計測する３次元位置計測システム及び方法に関する。
【背景技術】
【０００２】
航空機では、その安全のため、「着氷」を防止しており、そのため、主翼では、エンジンからの暖気（ブリード）を利用した防氷システムが一般に用いられている。しかしながら、エンジンから分岐する暖気供給量はエンジン燃費に直結する。航空機の安全性能と燃費性能を両立するためには、暖気供給量を適正に設定することが必要となるが、低減させすぎると、蒸発しきれず残る水分が固定翼に付着する氷、いわゆる、ランバックアイスが発生する。これは、航空機の安全性を考慮すると、注意が必要であり、防氷システムの暖機供給量の最適化には、ランバックアイスの高さ、形状を適切に評価することが必要となる。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００９−２７０９１５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
上記着氷を評価するための性能評価試験の１つとして、着氷風洞試験がある。当該試験は、通常の風洞試験と同様に閉鎖された風路に水滴を噴射し、試験断面として設置された供試体に積雪する氷の大きさ、形状を評価する試験である。
【０００５】
氷の形状を評価する場合、着氷の途中で通風を停止し、残った氷を切断して、型取する方法がある。この方法は、単純明快な方法であるが、用意された断面のみしか計測できない、氷が解け始める前に計測する必要がある、通風中の形状は計測できないなどの問題がある。また、近年では、通風停止後に残った氷を３次元計測器（例えば、３次元レーザトラッカなど）により実測する方法が通例である。この方法は、精度が高く、また、計測場所を限ることもないが、計測に非常に時間がかかる、氷が解け始める前に計測する必要がある、通風中の形状は計測できないなどの問題がある。このように、３次元計測器により実測する方法であっても、試験効率の低下や計測中の氷の熔解を伴う。
【０００６】
なお、脱落が発生する前の氷の高さを計測し、脱落直前の高さを予測（線形に外挿）することで、評価すべき高さを推定するなどの方法が考えられるが、この手法では、脱落直前の瞬間を見誤ることにより、航空機の安全性能に影響がある氷を見逃す可能性がある。また、見逃しを恐れて、設計上、過度な高さを推定することにより、航空機性能を過小評価する可能性があるため、これも適切ではない。また、そもそも、通風停止時の温度上昇や風速変動程度で脱落してしまうような条件の氷の形状は、上述した方法では計測できない。
【０００７】
本発明は上記課題に鑑みなされたもので、特殊な環境下でも、計測対象を高精度で効率良く計測することができる３次元位置計測システム及び方法を提供することを目的とする。
【課題を解決するための手段】
【０００８】
上記課題を解決する第１の発明に係る３次元位置計測システムは、
対象物を異なる方向から撮影する複数の撮影手段と、
前記複数の撮影手段で撮影した画像を取得し、前記画像の処理により、計測対象点の３次元位置を計測する処理手段とを有し、
前記処理手段は、
前記計測の前に、前記複数の撮影手段で３次元座標が既知である格子の画像を各々撮影し、撮影した前記画像を各々複数の基準格子画像として取得する計測前画像取得部と、
前記計測時に、前記複数の撮影手段で前記計測対象点を各々撮影し、撮影した前記画像を各々複数の計測画像として取得する計測時画像取得部と、
前記複数の基準格子画像の１つと、それに対応する前記複数の計測画像の１つにより、前記撮影手段から前記計測対象点への視線となる線を算出する線算出部と、
前記複数の線から前記計測対象点の３次元位置を算出する３次元位置算出部とを有する
ことを特徴とする。
【０００９】
上記課題を解決する第２の発明に係る３次元位置計測システムは、
対象物を互いに異なる方向から撮影する第１及び第２の撮影手段と、
前記第１及び第２の撮影手段で撮影した画像を取得し、前記画像の処理により、計測対象点の３次元位置を計測する処理手段とを有し、
前記処理手段は、
前記計測の前に、目盛となる基準格子を表面に有する格子定規を、３次元座標が既知である複数の設置位置に設置しながら、一端側から他端側へ順に１〜ｎ番目とする前記設置位置の各々において、前記第１及び第２の撮影手段で前記格子定規を各々撮影し、撮影した前記画像を各々第１及び第２の基準格子画像として取得する計測前画像取得部と、
前記計測時に、前記第１及び第２の撮影手段で前記計測対象点を各々撮影し、撮影した前記画像を各々第１及び第２の計測画像として取得する計測時画像取得部と、
ｉ番目とｊ番目の前記第１の基準格子画像の各々に前記第１の計測画像を重ね合わせ、各々の前記基準格子画像中の前記基準格子に対する前記第１の計測画像中の前記計測対象点の位置に基づいて、各々の前記基準格子上に前記計測対象点があると仮定した場合の前記計測対象点の３次元座標を各々ｉ番目とｊ番目の第１の対象点座標として算出すると共に、ｉ番目の前記第１の対象点座標とｊ番目の前記第１の対象点座標とを通る第１の直線を算出する第１直線算出部と、
ｉ番目とｊ番目の前記第２の基準格子画像の各々に前記第２の計測画像を重ね合わせ、各々の前記基準格子画像中の前記基準格子に対する前記第２の計測画像中の前記計測対象点の位置に基づいて、各々の前記基準格子上に前記計測対象点があると仮定した場合の前記計測対象点の３次元座標を各々ｉ番目とｊ番目の第２の対象点座標として算出すると共に、ｉ番目の前記第２の対象点座標とｊ番目の前記第２の対象点座標とを通る第２の直線を算出する第２直線算出部と、
前記第１の直線と前記第２の直線とを最短で結ぶ線分の距離又は前記線分の中点の３次元座標に基づいて、前記計測対象点の算出に最適な前記第１の直線と前記第２の直線を設定する直線設定部と、
最適な前記第１の直線と最適な前記第２の直線とを最短で結ぶ線分を求め、当該線分の中点の３次元座標を算出し、当該３次元座標を前記計測対象点の３次元位置とする３次元位置算出部とを有する
ことを特徴とする。
但し、ｎは２以上の自然数、ｉ＝１〜ｎの自然数、ｊ＝１〜ｎの自然数、ｉ≠ｊである。
【００１０】
上記課題を解決する第３の発明に係る３次元位置計測システムは、
上記第２の発明に記載の３次元位置計測システムにおいて、
前記処理手段は、
前記第１の直線と前記第２の直線とを最短で結ぶ線分の距離を算出する距離算出部と、
複数の前記距離の中から最小値を算出する最小値算出部とを有し、
前記直線設定部は、
前記第１直線算出部と前記第２直線算出部とを用いて、１〜ｎ−１の自然数であるｋの全てにおいて、ｉ＝ｋ及びｊ＝ｋ＋１として、ｋ番目の前記第１の直線とｋ番目の前記第２の直線とを各々算出し、
前記距離算出部を用いて、１〜ｎ−１の自然数であるｋの全てにおいて、ｋ番目の前記第１の直線とｋ番目の前記第２の直線とを最短で結ぶ線分のｋ番目の前記距離を各々算出し、
前記最小値算出部を用いて、１〜ｎ−１番目の前記距離の中から最小値となるｍ番目の前記距離を算出し、
最適な前記第１の直線と最適な前記第２の直線として、ｍ番目の前記第１の直線とｍ番目の前記第２の直線を設定する
ことを特徴とする。
但し、ｍ＝１〜ｎ−１の自然数である。
【００１１】
上記課題を解決する第４の発明に係る３次元位置計測システムは、
上記第２の発明に記載の３次元位置計測システムにおいて、
前記処理手段は、
前記第１の直線と前記第２の直線とを最短で結ぶ線分を求め、当該線分の中点の３次元座標を算出する中点座標算出部を有し、
前記直線設定部は、
前記第１直線算出部と前記第２直線算出部とを用いて、ｉ＝１及びｊ＝ｎとして、初期の前記第１の直線と初期の前記第２の直線とを各々算出し、
前記中点座標算出部を用いて、初期の前記第１の直線と初期の前記第２の直線とを最短で結ぶ初期の前記線分を求め、初期の前記線分の中点の３次元座標を算出し、
前記計測対象点の推定位置となる初期の前記線分の中点の３次元座標に基づいて、前記計測対象点を挟んで前記計測対象点に最も近い２つの前記設置位置で撮影したｍ番目とｍ＋１番目の前記第１及び第２の基準格子画像を選択し、
前記第１直線算出部と前記第２直線算出部とを用いて、ｉ＝ｍ及びｊ＝ｍ＋１として、最適な前記第１の直線と最適な前記第２の直線とを各々算出する
ことを特徴とする。
【００１２】
上記課題を解決する第５の発明に係る３次元位置計測システムは、
上記第１〜第４のいずれか１つの発明に記載の３次元位置計測システムにおいて、
前記対象物に所定のパターン映像を投影する投影手段を有し、
前記処理手段は、前記パターン映像に基づいて、同一の前記計測対象点を前記第１の計測画像中と前記第２の計測画像中で特定する対象点特定部を有する
ことを特徴とする。
【００１３】
上記課題を解決する第６の発明に係る３次元位置計測システムは、
上記第２〜第５のいずれか１つの発明に記載の３次元位置計測システムにおいて、
前記格子定規は、前記基準格子を表裏両方の表面に有し、
前記第１の撮影手段は前記格子定規の一方の表面を撮影する位置に設置され、前記第２の撮影手段は前記格子定規の他方の表面を撮影する位置に設置されている
ことを特徴とする。
【００１４】
上記課題を解決する第７の発明に係る３次元位置計測システムは、
上記第１〜第６のいずれか１つの発明に記載の３次元位置計測システムにおいて、
前記対象物は、航空機の翼の風洞模型であり、
前記計測対象点は、前記風洞模型に対する着氷風洞試験において、前記風洞模型の表面に付着する氷の頂点である
ことを特徴とする。
【００１５】
上記課題を解決する第８の発明に係る３次元位置計測方法は、
複数の撮影手段を用いて、対象物を異なる方向から撮影し、撮影した画像の処理により、計測対象点の３次元位置を計測する３次元位置計測方法において、
前記計測の前に、前記複数の撮影手段で３次元座標が既知である格子の画像を各々撮影し、撮影した前記画像を各々複数の基準格子画像として取得する計測前画像取得工程と、
前記計測時に、前記複数の撮影手段で前記計測対象点を各々撮影し、撮影した前記画像を各々複数の計測画像として取得する計測時画像取得工程と、
前記複数の基準格子画像の１つと、それに対応する前記複数の計測画像の１つにより、前記撮影手段から前記計測対象点への視線となる線を算出する線算出工程と、
前記複数の線から前記計測対象点の３次元位置を算出する３次元位置算出工程とを有する
ことを特徴とする。
【００１６】
上記課題を解決する第９の発明に係る３次元位置計測方法は、
第１及び第２の撮影手段を用いて、対象物を互いに異なる方向から撮影し、撮影した画像の処理により、計測対象点の３次元位置を計測する３次元位置計測方法において、
前記計測の前に、目盛となる基準格子を表面に有する格子定規を、３次元座標が既知である複数の設置位置に設置しながら、一端側から他端側へ順に１〜ｎ番目とする前記設置位置の各々において、前記第１及び第２の撮影手段で前記格子定規を各々撮影し、撮影した前記画像を各々第１及び第２の基準格子画像として取得する計測前画像取得工程と、
前記計測時に、前記第１及び第２の撮影手段で前記計測対象点を各々撮影し、撮影した前記画像を各々第１及び第２の計測画像として取得する計測時画像取得工程と、
ｉ番目とｊ番目の前記第１の基準格子画像の各々に前記第１の計測画像を重ね合わせ、各々の前記基準格子画像中の前記基準格子に対する前記第１の計測画像中の前記計測対象点の位置に基づいて、各々の前記基準格子上に前記計測対象点があると仮定した場合の前記計測対象点の３次元座標を各々ｉ番目とｊ番目の第１の対象点座標として算出すると共に、ｉ番目の前記第１の対象点座標とｊ番目の前記第１の対象点座標とを通る第１の直線を算出する第１直線算出工程と、
ｉ番目とｊ番目の前記第２の基準格子画像の各々に前記第２の計測画像を重ね合わせ、各々の前記基準格子画像中の前記基準格子に対する前記第２の計測画像中の前記計測対象点の位置に基づいて、各々の前記基準格子上に前記計測対象点があると仮定した場合の前記計測対象点の３次元座標を各々ｉ番目とｊ番目の第２の対象点座標として算出すると共に、ｉ番目の前記第２の対象点座標とｊ番目の前記第２の対象点座標とを通る第２の直線を算出する第２直線算出工程と、
前記第１の直線と前記第２の直線とを最短で結ぶ線分の距離又は前記線分の中点の３次元座標に基づいて、前記計測対象点の算出に最適な前記第１の直線と前記第２の直線を設定する直線設定工程と、
最適な前記第１の直線と最適な前記第２の直線とを最短で結ぶ線分を求め、当該線分の中点の３次元座標を算出し、当該３次元座標を前記計測対象点の３次元位置とする３次元位置算出工程とを有する
ことを特徴とする。
但し、ｎは２以上の自然数、ｉ＝１〜ｎの自然数、ｊ＝１〜ｎの自然数、ｉ≠ｊである。
【００１７】
上記課題を解決する第１０の発明に係る３次元位置計測方法は、
上記第９の発明に記載の３次元位置計測方法において、
前記第１の直線と前記第２の直線とを最短で結ぶ線分の距離を算出する距離算出工程と、
複数の前記距離の中から最小値を算出する最小値算出工程とを有し、
前記直線設定工程は、
前記第１直線算出工程と前記第２直線算出工程とを用いて、１〜ｎ−１の自然数であるｋの全てにおいて、ｉ＝ｋ及びｊ＝ｋ＋１として、ｋ番目の前記第１の直線とｋ番目の前記第２の直線とを各々算出し、
前記距離算出工程を用いて、１〜ｎ−１の自然数であるｋの全てにおいて、ｋ番目の前記第１の直線とｋ番目の前記第２の直線とを最短で結ぶ線分のｋ番目の前記距離を各々算出し、
前記最小値算出工程を用いて、１〜ｎ−１番目の前記距離の中から最小値となるｍ番目の前記距離を算出し、
最適な前記第１の直線と最適な前記第２の直線として、ｍ番目の前記第１の直線とｍ番目の前記第２の直線を設定する
ことを特徴とする。
但し、ｍ＝１〜ｎ−１の自然数である。
【００１８】
上記課題を解決する第１１の発明に係る３次元位置計測方法は、
上記第９の発明に記載の３次元位置計測方法において、
前記第１の直線と前記第２の直線とを最短で結ぶ線分を求め、当該線分の中点の３次元座標を算出する中点座標算出工程を有し、
前記直線設定工程は、
前記第１直線算出工程と前記第２直線算出工程とを用い、ｉ＝１及びｊ＝ｎとして、初期の前記第１の直線と初期の前記第２の直線とを各々算出し、
前記中点座標算出工程を用いて、初期の前記第１の直線と初期の前記第２の直線とを最短で結ぶ初期の前記線分を求め、初期の前記線分の中点の３次元座標を算出し、
前記計測対象点の推定位置となる初期の前記線分の中点の３次元座標に基づいて、前記計測対象点を挟んで前記計測対象点に最も近い２つの前記設置位置で撮影したｍ番目とｍ＋１番目の前記第１及び第２の基準格子画像を選択し、
前記第１直線算出工程と前記第２直線算出工程とを用いて、ｉ＝ｍ及びｊ＝ｍ＋１として、最適な前記第１の直線と最適な前記第２の直線とを各々算出する
ことを特徴とする。
【００１９】
上記課題を解決する第１２の発明に係る３次元位置計測方法は、
上記第８〜第１１のいずれか１つの発明に記載の３次元位置計測方法において、
前記対象物に所定のパターン映像を投影する投影手段を用い、前記パターン映像に基づいて、同一の前記計測対象点を前記第１の計測画像中と前記第２の計測画像中で特定する対象点特定工程を有する
ことを特徴とする。
【００２０】
上記課題を解決する第１３の発明に係る３次元位置計測方法は、
上記第９〜第１２のいずれか１つの発明に記載の３次元位置計測方法において、
前記格子定規は、前記基準格子を表裏両方の表面に有し、
前記第１の撮影手段は前記格子定規の一方の表面を撮影する位置に設置され、前記第２の撮影手段は前記格子定規の他方の表面を撮影する位置に設置されている
ことを特徴とする。
【００２１】
上記課題を解決する第１４の発明に係る３次元位置計測方法は、
上記第８〜第１３のいずれか１つの発明に記載の３次元位置計測方法において、
前記対象物は、航空機の翼の風洞模型であり、
前記計測対象点は、前記風洞模型に対する着氷風洞試験において、前記風洞模型の表面に付着する氷の頂点である
ことを特徴とする。
【発明の効果】
【００２２】
本発明によれば、特殊な環境下でも、計測対象を高精度で効率良く計測することができる。例えば、着氷風洞試験で計測を行う場合には、通風を停止することなく、着氷する氷の高さや形状をリアルタイムかつ高精度で計測することができ、また、任意の場所の氷を計測することができる。その結果、従来と比較して、計測効率が大幅に向上する。
【図面の簡単な説明】
【００２３】
【図１】（ａ）は、本発明に係る３次元位置計測システムの実施形態の一例（実施例１）を示すブロック図であり、（ｂ）は、その画像取得部を示すブロック図であり、（ｃ）は、その画像解析部を示すブロック図である。
【図２】図１に示した３次元位置計測システムで実施する基準格子画像の取得方法を示すフローチャートである。
【図３】（ａ）は、図２に示した基準格子画像の取得方法で使用する格子定規を示す正面図であり、（ｂ）は、図２に示した基準格子画像の取得方法において、撮影される格子定規の設置手順を説明する斜視図である。
【図４】取得した複数の基準格子画像の基準格子から構成される３次元格子を示す斜視図である。
【図５】図１に示した３次元位置計測システムで実施する３次元位置計測方法を示すフローチャートである。
【図６】図５に示した３次元位置計測方法を説明する斜視図である。
【図７】（ａ）は、図６に示した風洞模型の表面に着氷した氷の模式図であり、（ｂ）は、（ａ）に示した氷にグリッドパターンを投影した場合の模式図である。
【図８】図５に示した３次元位置計測方法において、計測対象点の撮影を説明する図である。
【図９】図５に示した３次元位置計測方法において、２つの基準格子への計測対象点のトレースを説明する図である。
【図１０】図５に示した３次元位置計測方法において、２つのカメラから計測対象点への２つの直線と１〜ｎ番目の基準格子との関係を説明する図である。
【図１１】図５に示した３次元位置計測方法において、２つの基準格子から求めた２つの直線の交点（線分の中点）を説明する図である。
【図１２】（ａ）は、本発明に係る３次元位置計測システムの実施形態の他の一例（実施例２）を示すブロック図であり、（ｂ）は、その画像取得部を示すブロック図であり、（ｃ）は、その画像解析部を示すブロック図である。
【図１３】図１２に示した３次元位置計測システムで実施する３次元位置計測方法を示すフローチャートの前半である。
【図１４】図１２に示した３次元位置計測システムで実施する３次元位置計測方法を示すフローチャートの後半である。
【発明を実施するための形態】
【００２４】
以降、本発明に係る３次元位置計測システム及び方法の実施形態について、図１〜図１４を参照して、説明を行う。
【００２５】
［実施例１］
まず、図１を参照して、本実施例の３次元位置計測システムの構成を説明する。ここで、図１（ａ）は、本実施例の３次元位置計測システムを示すブロック図であり、図１（ｂ）は、その画像取得部を示すブロック図であり、図１（ｃ）は、その画像解析部を示すブロック図である。
【００２６】
本実施例の３次元位置計測システムは、第１カメラ１１Ａ（第１の撮影手段）、第２カメラ１１Ｂ（第２の撮影手段）、パターン投影装置１２（投影手段）、データ処理装置２０（処理手段）、出力装置３０を有している。
【００２７】
第１カメラ１１Ａ及び第２カメラ１１Ｂは、計測対象をステレオ視するカメラであり、例えば、後述の図３（ｂ）に示すように、計測対象を互いに異なる方向から撮影している。第１カメラ１１Ａ及び第２カメラ１１Ｂで撮影した画像は、データ処理装置２０に送信される。
【００２８】
パターン投影装置１２は、プロジェクタなどの装置であり、第１カメラ１１Ａ及び第２カメラ１１Ｂによるステレオ視において、同一の計測対象点を特定するために、所定のパターンを計測対象に投影している。例えば、後述の図６に示すように、格子状のグリッドパターンＶを計測対象に投影している。グリッドパターンＶを投影する場合には、隣り合う線同士の色を互いに変えれば、同一の計測対象点の特定がより容易になる。なお、投影するパターンは、同一の計測対象点の特定ができれば、必ずしも、グリッドパターンＶである必要はなく、他のパターンや模様でも良い。
【００２９】
データ処理装置２０は、例えば、コンピュータなどの装置であり、装置構成としては、演算装置、記憶装置、入出力装置などからなり（図示省略）、機能構成としては、画像取得部２１、画像解析部２２Ａ、記憶部２３を有している。
【００３０】
データ処理装置２０において、画像取得部２１は、計測前画像取得部２１１（計測前画像取得工程）、計測時画像取得部２１２（計測時画像取得工程）を有している。画像取得部２１の各部については、後述の図２及び図５に示すフローチャートと共に説明するが、第１カメラ１１Ａ及び第２カメラ１１Ｂで撮影した画像を取得することで、後述する基準格子画像や計測画像の取得を行っている。
【００３１】
また、画像解析部２２Ａは、対象点特定部２２１（対象点特定工程）、第１直線算出部２２２（第１直線算出工程）、第２直線算出部２２３（第２直線算出工程）、直線設定部２２４（直線設定工程）、距離算出部２２５（距離算出工程）、最小値算出部２２６（最小値算出工程）、３次元位置算出部２２８（３次元位置算出工程）を有している。画像解析部２２Ａの各部についても、後述の図５に示すフローチャートと共に説明するが、第１カメラ１１Ａ及び第２カメラ１１Ｂから取得した基準格子画像や計測画像を処理することにより、計測対象点の３次元位置の計測を実施している。なお、特許請求の範囲の記載において、線算出部は、上記第１直線算出部２２２及び第２直線算出部２２３と対応し、線算出工程は、上記第１直線算出工程及び第２直線算出工程と対応している。
【００３２】
また、記憶部２３は、画像解析部２２Ａで実施する計測対象点の３次元位置の計測のための画像解析プログラムや後述する基準格子画像のデータを記憶している。
【００３３】
出力装置３０は、例えば、画像や情報を表示するディスプレイなどの装置であり、第１カメラ１１Ａ及び第２カメラ１１Ｂで撮影した画像や画像解析により得られた計測対象点の３次元位置を表示している。更に、出力装置３０として、画像や情報を印刷するプリンタなどの装置があっても良い。
【００３４】
次に、第１カメラ１１Ａ、第２カメラ１１Ｂ及びパターン投影装置１２の配置について、後述する図３及び図６を参照して説明する。第１カメラ１１Ａ、第２カメラ１１Ｂ及びパターン投影装置１２の配置は、後述する基準格子画像の取得時も、計測対象点の３次元位置の計測時も同じである。
【００３５】
風洞試験室（図示省略）の内側には、航空機の主翼の一部を試験断面として模型化した風洞模型Ｍが配置されている。第１カメラ１１Ａ、第２カメラ１１Ｂ及びパターン投影装置１２は、風洞試験室の外側に配置されており、風洞試験室の窓（窓ガラス）を介して、風洞模型Ｍの計測対象の領域を撮影したり、パターンを投影したりしている。
【００３６】
上述したように、第１カメラ１１Ａ及び第２カメラ１１Ｂは、互いに異なる方向から撮影しており、例えば、第１カメラ１１Ａは風洞試験室の側壁の窓から、第２カメラ１１Ｂは風洞試験室の天井の窓から撮影している。パターン投影装置１２は、計測対象の領域にパターンを投影できれば、第１カメラ１１Ａ又は第２カメラ１１Ｂと同じ窓でも良いし、これらとは別の窓でも良い。このような配置により、後述する基準格子画像の取得や計測対象点の３次元位置の計測が行われている。
【００３７】
［基準格子画像の取得］
計測対象点の３次元位置の計測の前に、基準格子画像の取得が行われる。
【００３８】
本実施例の３次元位置計測システムにおける基準格子画像の取得方法について、図１と共に、図２〜図４も参照して説明する。ここで、図２は、図１に示した３次元位置計測システムで実施する基準格子画像の取得方法を示すフローチャートである。また、図３（ａ）は、図２に示した基準格子画像の取得方法で使用する格子定規を示す正面図であり、図３（ｂ）は、図２に示した基準格子画像の取得方法において、撮影される格子定規の設置手順を説明する斜視図である。また、図４は、取得した複数の基準格子画像の基準格子から構成される３次元格子を示す斜視図である。
【００３９】
基準格子画像の取得には、図３（ａ）に示す格子定規４１を使用する。格子定規４１の表面には、目盛となる基準格子４１ａが形成されており、基準格子４１ａは、予め定めた一定間隔の縦線及び横線から構成されている。この基準格子４１ａは、格子定規４１の一方の表面のみにあっても良いが、より広角での撮影のためには、表裏両方の表面にあった方が良い。この場合、第１カメラ１１Ａを格子定規４１の一方の表面を撮影する位置に設置し、第２カメラ１１Ｂを格子定規４１の他方の表面を撮影する位置に設置すれば良い。
【００４０】
格子定規４１の表裏両方の表面に基準格子４１ａを形成する場合には、例えば、表側において、各縦線及び各横線を形成し、各縦線及び各横線の交点４１ｂの位置に貫通孔を形成し、裏側において、形成した貫通孔を基準にして、各縦線及び各横線を形成することで、格子定規４１の表裏で基準格子４１ａの各縦線及び各横線の位置を同じにしている。形成した貫通孔には、撮影時に交点４１ｂが認識され易いように、例えば、石膏などの目立つ物を詰めるようにしても良い。
【００４１】
また、風洞模型Ｍの表面には、設置位置を変える際に格子定規４１を設置し易いように、各々の設置位置を示す設置位置線４２が形成されている。ここでは、設置位置線４２の一例として、格子状の線が形成されているが、格子定規４１は、図３（ｂ）に示すように、一方向（例えば、図中に示す横方向）に順に移動させるので、設置位置線４２は縦方向の線のみで形成されていても良い。
【００４２】
そして、設置位置線４２の各々の設置位置には、後述する図１０に示すように、一端側から他端側へ順に１〜ｎ番目と順番が付されている。また、各々の設置位置は、風洞実験室の原点からの３次元座標が既知である。なお、ｎは２以上の自然数であり、計測対象の領域の大きさなどに応じて、予め規定されている。
【００４３】
上述したように、格子定規４１の基準格子４１ａは目盛の役割を果たし、設置位置線４２の各々の設置位置は３次元座標が既知であるので、各々の設置位置における格子定規４１の基準格子４１ａの各縦線、各横線及び各交点４１ｂは、その３次元座標が既知又は容易に算出可能である。このことを利用して、後述する計測対象点の３次元位置の計測が行われている。
【００４４】
また、風洞模型Ｍの表面には、格子定規４１の上端及び下端を支えるガイド４３が設けられている。このガイド４３は、基準格子画像の取得時のみ取り付け、計測対象点の３次元位置の計測時には取り外す。
【００４５】
図３（ａ）に示した格子定規４１を用い、図３（ｂ）に示すように、風洞模型Ｍの表面において、格子定規４１の設置位置を変えながら、各々の設置位置において、第１カメラ１１Ａ及び第２カメラ１１Ｂで撮影を行い、基準格子画像の取得が行われる。基準格子画像の取得の具体的な手順について、以下、図２に示すフローチャートを参照して説明する。なお、基準格子画像の取得の際には、通風前の通常の環境下で、つまり、通風も水滴の噴射も行われない環境下で行われる。
【００４６】
（ステップＳ１）
格子定規４１を設置位置線４２の１番目の設置位置に設置する。例えば、格子定規４１は、その表面が風洞実験室の３次元座標のＹＺ平面に平行になるように配置されている。この場合、設置位置線４２の２〜ｎ番目の設置位置に設置する場合には、格子定規４１をＸ軸方向に平行に移動すれば良い。例えば、図３（ｂ）では、格子定規４１が移動する横方向がＸ軸方向となる。
【００４７】
（ステップＳ２）
第１カメラ１１Ａで、設置位置線４２の１番目の設置位置の格子定規４１を撮影する。そして、撮影した格子定規４１の画像を１番目の第１の基準格子画像として計測前画像取得部２１１で取得し、記憶部２３へ記憶する。この１番目の第１の基準格子画像には、設置位置線４２の１番目の設置位置における格子定規４１の基準格子４１ａが撮影されている。
【００４８】
（ステップＳ３）
第１カメラ１１Ａと同様に、第２カメラ１１Ｂで、設置位置線４２の１番目の設置位置の格子定規４１を撮影する。そして、撮影した格子定規４１の画像を１番目の第２の基準格子画像として計測前画像取得部２１１で取得し、記憶部２３へ記憶する。この１番目の第２の基準格子画像には、設置位置線４２の１番目の設置位置における格子定規４１の基準格子４１ａが撮影されている。
【００４９】
なお、格子定規４１の表裏両方の表面に基準格子４１ａが形成されており、第１カメラ１１Ａが格子定規４１の一方の表面を撮影し、第２カメラ１１Ｂが格子定規４１の他方の表面を撮影する場合には、取得した第１の基準格子画像と第２の基準格子画像と間において、格子定規４１の厚さ分の座標を補正するようにすれば、計測誤差を抑制することができる。
【００５０】
（ステップＳ４）
上述したステップＳ１〜Ｓ３と同様に、格子定規４１を設置位置線４２の２〜ｎ番目の設置位置に順に設置し、２〜ｎ番目の設置位置の各々において、第１カメラ１１Ａ及び第２カメラ１１Ｂで格子定規４１を各々撮影し、撮影した格子定規４１の画像を２〜ｎ番目の第１及び第２の基準格子画像として計測前画像取得部２１１で各々取得し、記憶部２３へ各々記憶する。
【００５１】
上述した手順により、格子定規４１を設置位置線４２の１番目からｎ番目までの設置位置に順次設置しながら、各々の設置位置において、第１カメラ１１Ａ及び第２カメラ１１Ｂで撮影し、撮影した画像を１〜ｎ番目の第１及び第２の基準格子画像として取得し、記憶部２３へ記憶することになる。
【００５２】
この結果、取得した１〜ｎ番目の第１及び第２の基準格子画像により、これらの画像中に撮影された基準格子４１ａを基準格子ＢＧ1〜ＢＧnとすると、基準格子ＢＧ1〜ＢＧnからなる３次元格子群Ｇを風洞模型Ｍの表面の空間上に擬似的に作成することになる（図４及び図１０参照）。３次元格子群Ｇを構成する基準格子ＢＧ1〜ＢＧnは、３次元座標を示す３次元の定規の役割を果たすことになり、２つの基準格子ＢＧi及びＢＧjを用いて、計測対象点を内挿補間することにより、計測対象点の３次元位置を求めることができる。
【００５３】
なお、ここでは、上述した構成の格子定規４１を用いたが、３次元座標が既知であれば、他の格子状の物も使用可能であり、例えば、格子その物も使用可能である。
【００５４】
［計測対象の３次元位置計測］
基準格子画像の取得後、着氷風洞試験が行われ、計測対象の３次元位置の計測が行われることになる。
【００５５】
本実施例の３次元位置計測システムにおける３次元位置計測方法について、図１と共に、図５〜図１１も参照して説明する。ここで、図５は、図１に示した３次元位置計測システムで実施する３次元位置計測方法を示すフローチャートである。また、図６は、図５に示した３次元位置計測方法を説明する斜視図である。また、図７（ａ）は、図６に示した風洞模型の表面に着氷した氷の模式図であり、図７（ｂ）は、図７（ａ）に示した氷にグリッドパターンを投影した場合の模式図である。
【００５６】
また、図８〜図１１は、図５に示した３次元位置計測方法を説明する図であり、図８は、計測対象点の撮影を説明する図、図９は、２つの基準格子への計測対象点のトレースを説明する図、図１０は、各カメラから計測対象点への２つの直線と１〜ｎ番目の基準格子との関係を説明する図、図１１は、２つの基準格子から求めた２つの直線の交点（線分の中点）を説明する図である。
【００５７】
着氷風洞試験では、通風が開始され、その風路に水滴が噴射される。そのような環境下において、計測対象の３次元位置を計測する際には、図６に示すように、２つの第１カメラ１１Ａ及び第２カメラ１１Ｂで着氷の様子を捉える。この際、風洞模型Ｍの計測対象の領域にパターン投影装置１２からグリッドパターンＶを投影しており、第１カメラ１１Ａ及び第２カメラ１１Ｂは、投影されたグリッドパターンＶと共に、風洞模型Ｍの表面に付着する氷、即ち、ランバックアイスＲを撮影する。なお、図中の黒く塗りつぶした部分がランバックアイスである。
【００５８】
氷は半透明、無地であり、また、氷の表層は平坦であるので、そのままでは、氷上の計測対象点、例えば、着氷した氷の頂点を特定することは難しい。特に、ステレオ視の場合、同一の計測対象点を第１カメラ１１Ａ及び第２カメラ１１Ｂで各々特定する必要があるが、上述した理由により、各々のカメラで同一の計測対象点を特定することが難しい。そこで、本実施例では、パターン投影装置１２からグリッドパターンＶを投影することで、いわゆる、プロジェクションマッピングを行うことで、氷上の計測対象点の特定を可能としている。
【００５９】
例えば、図７（ａ）に示すようなランバックアイスＲが風洞模型Ｍの表面に付着したとする。氷は半透明、無地であり、その表層は平坦であるので、計測対象点Ｔを第１カメラ１１Ａ及び第２カメラ１１Ｂで各々特定することは難しい。本実施例では、上述したように、パターン投影装置１２からグリッドパターンＶを投影しているので、図７（ｂ）に示すように、グリッドパターンＶの縦線及び横線を目安にして、計測対象点Ｔを第１カメラ１１Ａ及び第２カメラ１１Ｂで各々特定するが可能になる。
【００６０】
このように、通風及び水滴の噴射が行われる着氷風洞試験では、図６に示すように、風洞模型Ｍの計測対象の領域にグリッドパターンＶを投影するパターン投影装置１２も用いて、計測対象の３次元位置の計測が行われる。計測対象の３次元位置の計測の具体的な手順について、以下、図５に示すフローチャートと共に、図８〜図１１も参照して説明する。
【００６１】
（ステップＳ１１）
第１カメラ１１Ａで、風洞模型Ｍの計測対象点Ｔを含む領域を撮影し、撮影した画像を第１の計測画像として計測時画像取得部２１２で取得し、取得した第１の計測画像中の計測対象点Ｔを対象点特定部２２１で特定する。
【００６２】
図８を参照して、ステップＳ１１を説明する。第１カメラ１１Ａで撮影した画像は、第１の計測画像ＩＡとして取得される。風洞模型Ｍの表面に付着したランバックアイスＲの計測対象点Ｔは、投影されたグリッドパターンＶにより、その位置が特定可能である。例えば、ランバックアイスＲの高さを計測したい場合には、ランバックアイスＲの頂点を計測対象点Ｔとすれば良い。そして、第１の計測画像ＩＡにおいて、特定した計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）が求まる。なお、第１の計測画像ＩＡにおいて、原点ＯA-pix、横軸ＸA-pix、縦軸ＹA-pixである。
【００６３】
（ステップＳ１２）
１番目の第１の基準格子画像中に撮影された格子定規４１の基準格子４１ａである基準格子ＢＧ1上に、計測対象点Ｔが存在すると仮定する。その場合、ステップＳ１１で求めた計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）を１番目の第１の基準格子画像中にトレースすることで、当該画像中の基準格子ＢＧ1に対する計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）に基づいて、仮定された計測対象点Ｔの３次元座標を１番目の第１の対象点座標ＰA1（ＸA1，ＹA1，ＺA1）として算出することができる。
【００６４】
図９を参照して、ステップＳ１２を説明する。１番目の第１の基準格子画像とステップＳ１１で撮影した第１の計測画像ＩＡとを重ね合わせることにより、１番目の基準格子ＢＧ1上に計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）をトレースすることができる（トレース画像ＴＡ１参照）。上述したように、格子定規４１の基準格子４１ａは目盛の役割を果たし、格子定規４１の設置位置は３次元座標が既知であるので、１番目の基準格子ＢＧ1の各線、各交点の３次元座標は既知又は容易に算出できる。この１番目の基準格子ＢＧ1に基づいて、仮定された計測対象点Ｔの対象点座標ＰA1（ＸA1，ＹA1，ＺA1）を算出することができる。
【００６５】
（ステップＳ１３）
２番目の第１の基準格子画像中に撮影された格子定規４１の基準格子４１ａである基準格子ＢＧ2上に、計測対象点Ｔが存在すると仮定する。その場合も、ステップＳ１１で求めた計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）を２番目の第１の基準格子画像中にトレースすることで、当該画像中の基準格子ＢＧ2に対する計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）に基づいて、仮定された計測対象点Ｔの３次元座標を２番目の第１の対象点座標ＰA2（ＸA2，ＹA2，ＺA2）として算出することができる。
【００６６】
再び、図９を参照して、ステップＳ１３を説明する。２番目の第１の基準格子画像とステップＳ１１で撮影した第１の計測画像ＩＡとを重ね合わせることにより、２番目の基準格子ＢＧ2上に計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）をトレースすることができる（トレース画像ＴＡ２参照）。上述したように、格子定規４１の基準格子４１ａは目盛の役割を果たし、格子定規４１の設置位置は３次元座標が既知であるので、２番目の基準格子ＢＧ2の各線、各交点の３次元座標は既知又は容易に算出できる。この２番目の基準格子ＢＧ2に基づいて、仮定された計測対象点Ｔの対象点座標ＰA2（ＸA2，ＹA2，ＺA2）を算出することができる。
【００６７】
（ステップＳ１４）
上述したステップＳ１２及びＳ１３により、３次元空間上の２点となる対象点座標ＰA1と対象点座標ＰA2が求められるので、これらを用いて、対象点座標ＰA1と対象点座標ＰA2を通る１番目の第１の直線ＬA1,2を算出する（図１１参照）。この第１の直線ＬA1,2は、第１カメラ１１Ａから計測対象点Ｔへの視線となる直線であり、言い換えれば、視線となる直線と１番目の基準格子ＢＧ1の交点が対象点座標ＰA1となり、視線となる直線と２番目の基準格子ＢＧ2の交点が対象点座標ＰA2となる（図１０参照）。このように、互いに隣接する位置となる１番目の基準格子ＢＧ1と２番目の基準格子ＢＧ2とを用いて、第１の直線ＬA1,2を算出している。
【００６８】
上述したステップＳ１２〜Ｓ１４が第１直線算出部２２２で実施される。この第１直線算出部２２２は、サブルーチンとして機能し、その場合、ｉ番目とｊ番目の第１の基準格子画像と第１の計測画像ＩＡとに基づいて、ｉ番目とｊ番目の第１の対象点座標を算出し、ｉ番目の第１の対象点座標とｊ番目の第１の対象点座標とを通る第１の直線を算出する。ここで、ｉ＝１〜ｎの自然数、ｊ＝１〜ｎの自然数、ｉ≠ｊである。従って、上述したステップＳ１２〜Ｓ１４では、ｉ＝１、ｊ＝２として、第１の直線ＬA1,2を算出している。
【００６９】
（ステップＳ１５）
第１カメラ１１Ａと同様に、第２カメラ１１Ｂで、風洞模型Ｍの計測対象点Ｔを含む領域を撮影し、撮影した画像を第２の計測画像として計測時画像取得部２１２で取得し、取得した第２の計測画像中の計測対象点Ｔを対象点特定部２２１で特定する。
【００７０】
再び、図８を参照して、ステップＳ１５を説明する。第２カメラ１１Ｂで撮影した画像は、第２の計測画像ＩＢとして取得される。風洞模型Ｍの表面に付着したランバックアイスＲの計測対象点Ｔは、投影されたグリッドパターンＶにより、その位置が特定可能である。つまり、第１カメラ１１Ａで特定した位置と同じ位置が、第２カメラ１１Ｂで特定可能である。そして、第２の計測画像ＩＢにおいて、特定した計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）が求まる。なお、第２の計測画像ＩＢにおいて、原点ＯB-pix、横軸ＸB-pix、縦軸ＹB-pixである。
【００７１】
（ステップＳ１６）
１番目の第２の基準格子画像中に撮影された格子定規４１の基準格子４１ａである基準格子ＢＧ1上に、計測対象点Ｔが存在すると仮定する。その場合、ステップＳ１５で求めた計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）を１番目の第２の基準格子画像中にトレースすることで、当該画像中の基準格子ＢＧ1に対する計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）に基づいて、仮定された計測対象点Ｔの３次元座標を１番目の第２の対象点座標ＰB1（ＸB1，ＹB1，ＺB1）として算出することができる。
【００７２】
再び、図９を参照して、ステップＳ１６を説明する。１番目の第２の基準格子画像とステップＳ１５で撮影した第２の計測画像ＩＢとを重ね合わせることにより、１番目の基準格子ＢＧ1上に計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）をトレースすることができる（トレース画像ＴＢ１参照）。上述したように、１番目の基準格子ＢＧ1の各線、各交点の３次元座標は既知又は容易に算出できるので、この１番目の基準格子ＢＧ1に基づいて、仮定された計測対象点Ｔの対象点座標ＰB1（ＸB1，ＹB1，ＺB1）を算出することができる。
【００７３】
（ステップＳ１７）
２番目の第２の基準格子画像中に撮影された格子定規４１の基準格子４１ａである基準格子ＢＧ2上に、計測対象点Ｔが存在すると仮定する。その場合も、ステップＳ１５で求めた計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）を２番目の第２の基準格子画像中にトレースすることで、当該画像中の基準格子ＢＧ2に対する計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）に基づいて、仮定された計測対象点Ｔの３次元座標を２番目の第２の対象点座標ＰB2（ＸB2，ＹB2，ＺB2）として算出することができる。
【００７４】
再び、図９を参照して、ステップＳ１７を説明する。２番目の第２の基準格子画像とステップＳ１５で撮影した第２の計測画像ＩＢとを重ね合わせることにより、２番目の基準格子ＢＧ2上に計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）をトレースすることができる（トレース画像ＴＢ２参照）。上述したように、２番目の基準格子ＢＧ2の各線、各交点の３次元座標は既知又は容易に算出できるので、この２番目の基準格子ＢＧ2に基づいて、仮定された計測対象点Ｔの対象点座標ＰB2（ＸB2，ＹB2，ＺB2）を算出することができる。
【００７５】
（ステップＳ１８）
上述したステップＳ１６及びＳ１７により、３次元空間上の２点となる対象点座標ＰB1と対象点座標ＰB2が求められるので、これらを用いて、対象点座標ＰB1と対象点座標ＰB2を通る１番目の第２の直線ＬB1,2を算出する（図１１参照）。この第２の直線ＬB1,2は、第２カメラ１１Ｂから計測対象点Ｔへの視線となる直線であり、言い換えれば、視線となる直線と１番目の基準格子ＢＧ1の交点が対象点座標ＰB1となり、視線となる直線と２番目の基準格子ＢＧ2の交点が対象点座標ＰB2となる（図１０参照）。このように、互いに隣接する位置となる１番目の基準格子ＢＧ1と２番目の基準格子ＢＧ2とを用いて、第２の直線ＬB1,2を算出している。
【００７６】
上述したステップＳ１６〜Ｓ１８が第２直線算出部２２３で実施される。この第２直線算出部２２３も、サブルーチンとして機能し、その場合、ｉ番目とｊ番目の第２の基準格子画像と第２の計測画像ＩＢとに基づいて、ｉ番目とｊ番目の第２の対象点座標を算出し、ｉ番目の第２の対象点座標とｊ番目の第２の対象点座標とを通る第２の直線を算出する。ここでも、ｉ＝１〜ｎの自然数、ｊ＝１〜ｎの自然数、ｉ≠ｊである。従って、上述したステップＳ１６〜Ｓ１８では、ｉ＝１、ｊ＝２として、第２の直線ＬB1,2を算出している。
【００７７】
そして、上述したステップＳ１１〜Ｓ１８により、１番目の基準格子ＢＧ1と２番目の基準格子ＢＧ2の間において、１番目の第１の直線ＬA1,2と第２の直線ＬB1,2が算出されることになる。
【００７８】
（ステップＳ１９）
１番目の基準格子ＢＧ1と２番目の基準格子ＢＧ2の間において、１番目の第１の直線ＬA1,2と第２の直線ＬB1,2とを結ぶ最短の線分の距離である最小距離ｈ1,2を距離算出部２２５で算出する。
【００７９】
（ステップＳ２０）
２〜ｎ番目の第１及び第２の基準格子画像についても、上述したステップＳ１１〜Ｓ１９と同様の手順を行って、２〜ｎ−１番目の第１の直線ＬAk,k+1と第２の直線ＬBk,k+1を各々算出する。そして、２番目の基準格子ＢＧ2と３番目の基準格子ＢＧ3の間、３番目の基準格子ＢＧ3と４番目の基準格子ＢＧ4の間、・・・、ｎ−１番目の基準格子ＢＧn-1とｎ番目の基準格子ＢＧnの間においても、各々、２〜ｎ−１番目の最小距離ｈk,k+1を算出する。なお、ここでのｋは、２〜ｎ−１の自然数である。
【００８０】
つまり、ステップＳ１１〜Ｓ２０では、１〜ｎ−１の自然数ｋの全てにおいて、第１直線算出部２２２と第２直線算出部２２３とを用いて、ｉ＝ｋ及びｊ＝ｋ＋１として、ｋ番目の第１の直線ＬAk,k+1とｋ番目の第２の直線ＬBk,k+1とを各々算出し、そして、距離算出部２２５を用いて、ｋ番目の第１の直線ＬAk,k+1とｋ番目の第２の直線ＬBk,k+1とを最短で結ぶ線分のｋ番目の最小距離ｈk,k+1を各々算出している。
【００８１】
（ステップＳ２１）
１〜ｎ−１番目の最小距離ｈ1,2 〜ｈn-1,nの中から、これらの最小値となるｍ番目の最小距離ｈm,m+1を最小値算出部２２６で算出する。これは、計測対象点Ｔがｍ番目の基準格子ＢＧmとｍ＋１番目の基準格子ＢＧm+1と間に存在することを意味する。ここで、ｍ＝１〜ｎ−１の自然数である。
【００８２】
上述したステップＳ１１〜Ｓ２１が直線設定部２２４で実施される。つまり、直線設定部２２４では、ｋ＝１〜ｎ−１とすると、第１直線算出部２２２と第２直線算出部２２３とで算出した１〜ｎ−１番目の第１の直線ＬAk,k+1と第２の直線ＬBk,k+1の中から、距離算出部２２５及び最小値算出部２２６を経て算出した最小値となる最小距離ｈm,m+1に基づいて、計測対象点Ｔの算出に最適な第１の直線ＬAm,m+1と第２の直線ＬBm,m+1とを設定している。
【００８３】
（ステップＳ２２）
計測対象点Ｔの算出に最適な第１の直線ＬAm,m+1と第２の直線ＬBm,m+1に基づいて、第１の直線ＬAm,m+1と第２の直線ＬBm,m+1を結ぶ最短の線分（最小距離ｈm,m+1）の中点の３次元座標（ＸT，ＹT，ＺT）を３次元位置算出部２２８で算出する。この３次元座標（ＸT，ＹT，ＺT）が、求める計測対象点Ｔの３次元位置となる。
【００８４】
なお、理想的には、最小距離ｈm,m+1は「０」となり、その場合、第１の直線ＬAm,m+1と第２の直線ＬBm,m+1は交差し、この交点が求める計測対象点Ｔとなる。しかしながら、誤差やレンズの歪みなどのカメラパラーメータの影響により、第１の直線ＬAm,m+1と第２の直線ＬBm,m+1は、厳密には、直線ではなく、曲線となり、最小距離ｈm,m+1が「０」となることは少ない。
【００８５】
そのため、第１の直線ＬAm,m+1と第２の直線ＬBm,m+1を結ぶ最短の線分の中点を、求める計測対象点Ｔとしている。この場合、最小となる最小距離ｈm,m+1を有する第１の直線ＬAm,m+1と第２の直線ＬBm,m+1を選択することは、計測対象点Ｔを挟んで計測対象点Ｔに最も近い２つの基準格子ＢＧmと基準格子ＢＧm+1を選択することになる。つまり、計測対象点Ｔは、計測対象点Ｔを挟んで計測対象点Ｔに最も近い２つの基準格子ＢＧmと基準格子ＢＧm+1とを用いて、内挿補間により求められた内挿点となる。そのため、最も精度よく位置を算出することができ、位置精度の向上を図ることができる。
【００８６】
従来、風洞模型に対する着氷風洞試験では、風洞通風中の着氷形状を計測することはできなかった。着氷風洞試験で着氷形状を計測するためには、風洞通風中に風洞試験室の外からステレオ視を行えば良いのだが、風洞試験室の窓を介して、つまり、密度媒体の異なる窓ガラスを透過して行うため、一般的なステレオ視では、カメラパラメータ（レンズ歪曲など）の影響を見積もらなければならず、計測精度を担保することが困難であった。
【００８７】
これに対して、本実施例の３次元位置計測システムは、一般的なステレオ視による方法ではなく、計測前に作成された複数の基準格子ＢＧ1〜ＢＧnを用いて、計測対象点Ｔを内挿補間することにより求めている。従って、本実施例の３次元位置計測システムを用いることにより、特殊な環境下、例えば、風洞試験室での着氷風洞試験で着氷を計測する場合でも、着氷の高さ、形状を高精度で計測することが可能となる。また、閉鎖された風洞内に設置された風洞模型において、任意の場所かつ時々刻々と変化する氷の形状を高精度でリアルタイムに計測することが可能となる。また、従来のように、通風停止作業や形状取得のための計測時間がなくなり、試験効率が大幅に向上する。
【００８８】
［実施例２］
図１２（ａ）は、本実施例の３次元位置計測システムを示すブロック図であり、図１２（ｂ）は、その画像取得部を示すブロック図であり、図１２（ｃ）は、その画像解析部を示すブロック図である。また、図１３〜図１４は、図１２に示した３次元位置計測システムで実施する３次元位置計測方法を示すフローチャートであり、図１３は、その前半、図１４は、その後半である。
【００８９】
本実施例の３次元位置計測システムは、画像解析部２２Ｂを除き、図１に示した実施例１の３次元位置計測システムと同等の構成で良い。従って、同等の構成には同じ符号を付し、重複する説明は省略する。
【００９０】
本実施例の３次元位置計測システムにおいて、画像解析部２２Ｂは、対象点特定部２２１、第１直線算出部２２２、第２直線算出部２２３、直線設定部２２４、中点座標算出部２２７（中点座標算出工程）、３次元位置算出部２２８を有している。つまり、実施例１の３次元位置計測システムの画像解析部２２Ａにおける距離算出部２２５及び最小値算出部２２６に代えて、中点座標算出部２２７を備えている。この画像解析部２２Ｂの各部については、後述の図１３及び図１４に示すフローチャートと共に説明する。
【００９１】
また、本実施例の３次元位置計測方法で用いる基準格子画像の取得方法は、実施例１で説明した方法で良いので、基準格子画像の取得方法についても、その説明を省略する。
【００９２】
［計測対象の３次元位置計測］
実施例１と同様に、基準格子画像の取得後、着氷風洞試験が行われ、計測対象の３次元位置の計測が行われる。
【００９３】
本実施例の３次元位置計測システムにおける３次元位置計測方法について、図１２と共に、図１３及び図１４、そして、上述した図８〜図１１も参照して説明を行う。なお、本実施例の３次元位置計測方法において、着氷風洞試験については、実施例１の図６及び図７で説明した通りであるので、着氷風洞試験についても、その説明を省略する。
【００９４】
（ステップＳ３１）
第１カメラ１１Ａで、風洞模型Ｍの計測対象点Ｔを含む領域を撮影し、撮影した画像を第１の計測画像として計測時画像取得部２１２で取得し、取得した第１の計測画像中の計測対象点Ｔを対象点特定部２２１で特定する。
【００９５】
上述した図８を参照して、ステップＳ３１を説明する。第１カメラ１１Ａで撮影した画像は、第１の計測画像ＩＡとして取得される。風洞模型Ｍの表面に付着したランバックアイスＲの計測対象点Ｔは、投影されたグリッドパターンＶにより、その位置が特定可能である。そして、第１の計測画像ＩＡにおいて、特定した計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）が求まる。
【００９６】
（ステップＳ３２）
１番目の第１の基準格子画像中に撮影された格子定規４１の基準格子４１ａである基準格子ＢＧ1上に、計測対象点Ｔが存在すると仮定する。その場合、ステップＳ３１で求めた計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）を１番目の第１の基準格子画像中にトレースすることで、当該画像中の基準格子ＢＧ1に対する計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）に基づいて、仮定された計測対象点Ｔの３次元座標を第１の対象点座標ＰA1（ＸA1，ＹA1，ＺA1）として算出することができる。
【００９７】
上述した図９を参照して、ステップＳ３２を説明する。１番目の第１の基準格子画像とステップＳ３１で撮影した第１の計測画像ＩＡとを重ね合わせることにより、１番目の基準格子ＢＧ1上に計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）をトレースすることができる（トレース画像ＴＡ１参照）。上述したように、１番目の基準格子ＢＧ1の各線、各交点の３次元座標は既知又は容易に算出できるので、この１番目の基準格子ＢＧ1に基づいて、仮定された計測対象点Ｔの対象点座標ＰA1（ＸA1，ＹA1，ＺA1）を算出することができる。
【００９８】
（ステップＳ３３）
ｎ番目の第１の基準格子画像中に撮影された格子定規４１の基準格子４１ａである基準格子ＢＧn上に、計測対象点Ｔが存在すると仮定する。その場合も、ステップＳ３１で求めた計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）をｎ番目の第１の基準格子画像中にトレースすることで、当該画像中の基準格子ＢＧnに対する計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）に基づいて、仮定された計測対象点Ｔの３次元座標を第１の対象点座標ＰAn（ＸA2，ＹAn，ＺAn）として算出することができる。
【００９９】
再び、図９を参照して、ステップＳ３３を説明する。但し、ここでは、図９中の「２」を「ｎ」に読み替える。ｎ番目の第１の基準格子画像とステップＳ３１で撮影した第１の計測画像ＩＡとを重ね合わせることにより、ｎ番目の基準格子ＢＧn上に計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）をトレースすることができる。上述したように、格子定規４１の基準格子４１ａは目盛の役割を果たし、格子定規４１の設置位置は３次元座標が既知であるので、ｎ番目の基準格子ＢＧnの各線、各交点の３次元座標は既知又は容易に算出できる。このｎ番目の基準格子ＢＧnに基づいて、仮定された計測対象点Ｔの対象点座標ＰAn（ＸA2，ＹAn，ＺAn）を算出することができる。
【０１００】
（ステップＳ３４）
上述したステップＳ３２及びＳ３３により、３次元空間上の２点となる対象点座標ＰA1と対象点座標ＰAnが求められるので、これらを用いて、対象点座標ＰA1と対象点座標ＰAnを通る第１の直線ＬA1,nを算出する（図１１参照）。但し、ここでは、図１１中の「２」を「ｎ」に読み替える。第１の直線ＬA1,nは、第１カメラ１１Ａから計測対象点Ｔへの視線となる直線であり、言い換えれば、視線となる直線と１番目の基準格子ＢＧ1の交点が対象点座標ＰA1となり、視線となる直線とｎ番目の基準格子ＢＧnの交点が対象点座標ＰAnとなる（図１０参照）。このように、両端の位置となる１番目の基準格子ＢＧ1とｎ番目の基準格子ＢＧnとを用いて、第１の直線ＬA1,nを算出している。
【０１０１】
上述したステップＳ３２〜Ｓ３４が第１直線算出部２２２で実施される。サブルーチンとして機能する第１直線算出部２２２は、上述したステップＳ３２〜Ｓ３４では、ｉ＝１、ｊ＝ｎとして、第１の直線ＬA1,nを算出している。
【０１０２】
（ステップＳ３５）
第１カメラ１１Ａと同様に、第２カメラ１１Ｂで、風洞模型Ｍの計測対象点Ｔを含む領域を撮影し、撮影した画像を第２の計測画像として計測時画像取得部２１２で取得し、取得した第２の計測画像中の計測対象点Ｔを対象点特定部２２１で特定する。
【０１０３】
再び、図８を参照して、ステップＳ３５を説明する。第２カメラ１１Ｂで撮影した画像は、第２の計測画像ＩＢとして取得される。風洞模型Ｍの表面に付着したランバックアイスＲの計測対象点Ｔは、投影されたグリッドパターンＶにより、その位置が特定可能である。つまり、第１カメラ１１Ａで特定した位置と同じ位置が、第２カメラ１１Ｂで特定可能である。そして、第２の計測画像ＩＢにおいて、特定した計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）が求まる。
【０１０４】
（ステップＳ３６）
１番目の第２の基準格子画像中に撮影された格子定規４１の基準格子４１ａである基準格子ＢＧ1上に、計測対象点Ｔが存在すると仮定する。その場合、ステップＳ３５で求めた計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）を１番目の第２の基準格子画像中にトレースすることで、当該画像中の基準格子ＢＧ1に対する計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）に基づいて、仮定された計測対象点Ｔの３次元座標を第２の対象点座標ＰB1（ＸB1，ＹB1，ＺB1）として算出することができる。
【０１０５】
再び、図９を参照して、ステップＳ３６を説明する。１番目の第２の基準格子画像とステップＳ３５で撮影した第２の計測画像ＩＢとを重ね合わせることにより、１番目の基準格子ＢＧ1上に計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）をトレースすることができる（トレース画像ＴＢ１参照）。上述したように、１番目の基準格子ＢＧ1の各線、各交点の３次元座標は既知又は容易に算出できるので、この１番目の基準格子ＢＧ1に基づいて、仮定された計測対象点Ｔの対象点座標ＰB1（ＸB1，ＹB1，ＺB1）を算出することができる。
【０１０６】
（ステップＳ３７）
ｎ番目の第２の基準格子画像中に撮影された格子定規４１の基準格子４１ａである基準格子ＢＧn上に、計測対象点Ｔが存在すると仮定する。その場合も、ステップＳ３５で求めた計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）をｎ番目の第２の基準格子画像中にトレースすることで、当該画像中の基準格子ＢＧnに対する計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）に基づいて、仮定された計測対象点Ｔの３次元座標を２番目の第２の対象点座標ＰBn（ＸBn，ＹBn，ＺBn）として算出することができる。
【０１０７】
再び、図９を参照して、ステップＳ３７を説明する。但し、ここでは、図９中の「２」を「ｎ」に読み替える。ｎ番目の第２の基準格子画像とステップＳ３５で撮影した第２の計測画像ＩＢとを重ね合わせることにより、ｎ番目の基準格子ＢＧn上に計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）をトレースすることができる。上述したように、ｎ番目の基準格子ＢＧnの各線、各交点の３次元座標は既知又は容易に算出できるので、このｎ番目の基準格子ＢＧnに基づいて、仮定された計測対象点Ｔの対象点座標ＰBn（ＸBn，ＹBn，ＺBn）を算出することができる。
【０１０８】
（ステップＳ３８）
上述したステップＳ３６及びＳ３７により、３次元空間上の２点となる対象点座標ＰB1と対象点座標ＰBnが求められるので、これらを用いて、対象点座標ＰB1と対象点座標ＰBnを通る第２の直線ＬB1,nを算出する（図１１参照）。但し、ここでは、図１１中の「２」を「ｎ」に読み替える。第２の直線ＬB1,nは、第２カメラ１１Ｂから計測対象点Ｔへの視線となる直線であり、言い換えれば、視線となる直線と１番目の基準格子ＢＧ1の交点が対象点座標ＰB1となり、視線となる直線とｎ番目の基準格子ＢＧnの交点が対象点座標ＰBnとなる（図１０参照）。このように、両端の位置となる１番目の基準格子ＢＧ1とｎ番目の基準格子ＢＧnとを用いて、第２の直線ＬB1,nを算出している。
【０１０９】
上述したステップＳ３６〜Ｓ３８が第２直線算出部２２３で実施される。サブルーチンとして機能する第２直線算出部２２３は、上述したステップＳ３６〜Ｓ３８では、ｉ＝１、ｊ＝ｎとして、第２の直線ＬB1,nを算出している。
【０１１０】
そして、上述したステップＳ３１〜Ｓ３８により、１番目の基準格子ＢＧ1とｎ番目の基準格子ＢＧnの間において、第１の直線ＬA1,nと第２の直線ＬB1,nが算出されることになる。
【０１１１】
（ステップＳ３９）
１番目の基準格子ＢＧ1とｎ番目の基準格子ＢＧnの間において、第１の直線ＬA1,nと第２の直線ＬB1,nとを結ぶ最短の線分を求め、求めた線分の中点の３次元座標を中点座標算出部２２７で算出する。このようにして、１番目の基準格子ＢＧ1とｎ番目の基準格子ＢＧnから、つまり、両端の２つの基準格子ＢＧ1と基準格子ＢＧnから、初期の２つの直線となる第１の直線ＬA1,n及び第２の直線ＬB1,nを求め、第１の直線ＬA1,nと第２の直線ＬB1,nとを結ぶ最短の線分の中点の３次元座標を求めており、これが、計測対象点Ｔの初期の推定位置となる。
【０１１２】
（ステップＳ４０）
ステップＳ３９で求めた線分の中点の３次元座標に基づいて、この中点、即ち、初期の推定位置となる計測対象点Ｔを挟んで計測対象点Ｔに最も近い位置となる２つの設置位置で撮影した格子定規４１の画像として、ｍ番目とｍ＋１番目の第１及び第２の基準格子画像を求めている。つまり、ｍ番目の基準格子ＢＧmとｍ＋１番目の基準格子ＢＧm+1を求めている。
【０１１３】
（ステップＳ４１）
ｍ番目の第１の基準格子画像中に撮影された格子定規４１の基準格子４１ａである基準格子ＢＧm上に、計測対象点Ｔが存在すると仮定する。その場合、ステップＳ３１で求めた計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）をｍ番目の第１の基準格子画像中にトレースすることで、当該画像中の基準格子ＢＧmに対する計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）に基づいて、仮定された計測対象点Ｔの３次元座標を第１の対象点座標ＰAm（ＸAm，ＹAm，ＺAm）として算出することができる。
【０１１４】
再び、図９を参照して、ステップＳ４１を説明する。但し、ここでは、図９中の「１」を「ｍ」に読み替える。ｍ番目の第１の基準格子画像とステップＳ３１で撮影した第１の計測画像ＩＡとを重ね合わせることにより、ｍ番目の基準格子ＢＧm上に計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）をトレースすることができる。上述したように、格子定規４１の基準格子４１ａは目盛の役割を果たし、格子定規４１の設置位置は３次元座標が既知であるので、ｍ番目の基準格子ＢＧmの各線、各交点の３次元座標は既知又は容易に算出できる。このｍ番目の基準格子ＢＧmに基づいて、仮定された計測対象点Ｔの対象点座標ＰAm（ＸAm，ＹAm，ＺAm）を算出することができる。
【０１１５】
（ステップＳ４２）
ｍ＋１番目の第１の基準格子画像中に撮影された格子定規４１の基準格子４１ａである基準格子ＢＧm+1上に、計測対象点Ｔが存在すると仮定する。その場合も、ステップＳ３１で求めた計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）をｍ＋１番目の第１の基準格子画像中にトレースすることで、当該画像中の基準格子ＢＧm+1に対する計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）に基づいて、仮定された計測対象点Ｔの３次元座標を第１の対象点座標ＰAm+1（ＸAm+1，ＹAm+1，ＺAm+1）として算出することができる。
【０１１６】
再び、図９を参照して、ステップＳ４２を説明する。但し、ここでは、図９中の「２」を「ｍ＋１」に読み替える。ｍ＋１番目の第１の基準格子画像とステップＳ３１で撮影した第１の計測画像ＩＡとを重ね合わせることにより、ｍ＋１番目の基準格子ＢＧm+1上に計測対象点Ｔのピクセル座標（ｘpix,A,t、ｙpix,A,t）をトレースすることができる。上述したように、格子定規４１の基準格子４１ａは目盛の役割を果たし、格子定規４１の設置位置は３次元座標が既知であるので、ｍ＋１番目の基準格子ＢＧm+1の各線、各交点の３次元座標は既知又は容易に算出できる。このｍ＋１番目の基準格子ＢＧm+1に基づいて、仮定された計測対象点Ｔの対象点座標ＰAm+1（ＸAm+1，ＹAm+1，ＺAm+1）を算出することができる。
【０１１７】
（ステップＳ４３）
上述したステップＳ４１及びＳ４２により、３次元空間上の２点となる対象点座標ＰAmと対象点座標ＰAm+1が求められるので、これらを用いて、対象点座標ＰAmと対象点座標ＰAm+1を通る第１の直線ＬAm,m+1を算出する（図１１参照）。但し、ここでは、図１１中の「１」を「ｍ」に、「２」を「ｍ＋１」に読み替える。第１の直線ＬAm,m+1は、第１カメラ１１Ａから計測対象点Ｔへの視線となる直線であり、言い換えれば、視線となる直線とｍ番目の基準格子ＢＧmの交点が対象点座標ＰAmとなり、視線となる直線とｍ＋１番目の基準格子ＢＧm+1の交点が対象点座標ＰAm+1となる（図１０参照）。このように、計測対象点Ｔを挟んで計測対象点Ｔに最も近いｍ番目の基準格子ＢＧmとｍ＋１番目の基準格子ＢＧm+1とを用いて、計測対象点Ｔの算出に最適な第１の直線ＬAm,m+1を算出している。
【０１１８】
上述したステップＳ４１〜Ｓ４３が第１直線算出部２２２で実施される。サブルーチンとして機能する第１直線算出部２２２は、上述したステップＳ４１〜Ｓ４３では、ｉ＝ｍ、ｊ＝ｍ＋１として、第１の直線ＬAm,m+1を算出している。
【０１１９】
（ステップＳ４４）
ｍ番目の第２の基準格子画像中に撮影された格子定規４１の基準格子４１ａである基準格子ＢＧm上に、計測対象点Ｔが存在すると仮定する。その場合も、ステップＳ３５で求めた計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）をｍ番目の第２の基準格子画像中にトレースすることで、当該画像中の基準格子ＢＧmに対する計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）に基づいて、仮定された計測対象点Ｔの３次元座標を第２の対象点座標ＰBm（ＸBm，ＹBm，ＺBm）として算出することができる。
【０１２０】
再び、図９を参照して、ステップＳ４４を説明する。但し、ここでは、図９中の「１」を「ｍ」に読み替える。ｍ番目の第２の基準格子画像とステップＳ３５で撮影した第２の計測画像ＩＢとを重ね合わせることにより、ｍ番目の基準格子ＢＧm上に計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）をトレースすることができる。上述したように、ｍ番目の基準格子ＢＧmの各線、各交点の３次元座標は既知又は容易に算出できるので、このｍ番目の基準格子ＢＧmに基づいて、仮定された計測対象点Ｔの対象点座標ＰBm（ＸBm，ＹBm，ＺBm）を算出することができる。
【０１２１】
（ステップＳ４５）
ｍ＋１番目の第２の基準格子画像中に撮影された格子定規４１の基準格子４１ａである基準格子ＢＧm+1上に、計測対象点Ｔが存在すると仮定する。その場合も、ステップＳ３５で求めた計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）をｍ＋１番目の第２の基準格子画像中にトレースすることで、当該画像中の基準格子ＢＧm+1に対する計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）に基づいて、仮定された計測対象点Ｔの３次元座標を第２の対象点座標ＰBm+1（ＸBm+1，ＹBm+1，ＺBm+1）として算出することができる。
【０１２２】
再び、図９を参照して、ステップＳ４５を説明する。但し、ここでは、図９中の「２」を「ｍ＋１」に読み替える。ｍ＋１番目の第２の基準格子画像とステップＳ３５で撮影した第２の計測画像ＩＢとを重ね合わせることにより、ｍ＋１番目の基準格子ＢＧm+1上に計測対象点Ｔのピクセル座標（ｘpix,B,t、ｙpix,B,t）をトレースすることができる。上述したように、ｍ＋１番目の基準格子ＢＧm+1の各線、各交点の３次元座標は既知又は容易に算出できるので、このｍ＋１番目の基準格子ＢＧm+1に基づいて、仮定された計測対象点Ｔの対象点座標ＰBm+1（ＸBm+1，ＹBm+1，ＺBm+1）を算出することができる。
【０１２３】
（ステップＳ４６）
上述したステップＳ４４及びＳ４５により、３次元空間上の２点となる対象点座標ＰBmと対象点座標ＰBm+1が求められるので、これらを用いて、対象点座標ＰBmと対象点座標ＰBm+1を通る第２の直線ＬBm,m+1を算出する（図１１参照）。但し、ここでは、図１１中の「１」を「ｍ」に、「２」を「ｍ＋１」に読み替える。第２の直線ＬBm,m+1は、第２カメラ１１Ｂから計測対象点Ｔへの視線となる直線であり、言い換えれば、視線となる直線とｍ番目の基準格子ＢＧmの交点が対象点座標ＰBmとなり、視線となる直線とｍ＋１番目の基準格子ＢＧm+1の交点が対象点座標ＰBm+1となる（図１０参照）。このように、計測対象点Ｔを挟んで計測対象点Ｔに最も近いｍ番目の基準格子ＢＧmとｍ＋１番目の基準格子ＢＧm+1とを用いて、計測対象点Ｔの算出に最適な第２の直線ＬBm,m+1を算出している。
【０１２４】
上述したステップＳ４４〜Ｓ４６が第２直線算出部２２３で実施される。サブルーチンとして機能する第２直線算出部２２３は、上述したステップＳ４４〜Ｓ４６では、ｉ＝ｍ、ｊ＝ｍ＋１として、第２の直線ＬBm,m+1を算出している。
【０１２５】
そして、上述したステップＳ３１〜Ｓ４６が直線設定部２２４で実施される。つまり、直線設定部２２４では、第１直線算出部２２２と第２直線算出部２２３とを用いて、ｉ＝１及びｊ＝ｎとして、初期の第１の直線ＬA1,nと初期の第２の直線ＬB1,nとを各々算出し、中点座標算出部２２７を用いて、初期の第１の直線ＬA1,nと初期の第２の直線ＬB1,nとを最短で結ぶ初期の線分を求め、初期の線分の中点の３次元座標を算出している。そして、計測対象点Ｔの推定位置となる初期の線分の中点の３次元座標に基づいて、計測対象点をＴ挟んで計測対象点Ｔに最も近い２つの設置位置で撮影したｍ番目とｍ＋１番目の第１及び第２の基準格子画像を選択し、再び、第１直線算出部２２２と第２直線算出部２２３とを用いて、ｉ＝ｍ及びｊ＝ｍ＋１として、計測対象点Ｔの算出に最適な第１の直線ＬAm,m+1と第２の直線ＬBm,m+1とを各々算出している。
【０１２６】
（ステップＳ４７）
ステップＳ４３で算出した第１の直線ＬAm,m+1とステップＳ４６で算出した第２のＬBm,m+1を結ぶ最短の線分の中点の３次元座標（ＸT，ＹT，ＺT）を３次元位置算出部２２８で算出する。この３次元座標（ＸT，ＹT，ＺT）が、求める計測対象点Ｔの３次元位置となる。
【０１２７】
本実施例でも、計測前に作成された複数の基準格子ＢＧ1〜ＢＧnを用いており、計測対象点Ｔは、計測対象点Ｔを挟んで最近傍の位置となる２つの基準格子ＢＧmと基準格子ＢＧm+1とを用いて、内挿補間により求められた内挿点となる。そのため、最も精度よく位置を算出することができ、位置精度の向上を図ることができる。従って、本実施例の３次元位置計測システムを用いることにより、実施例１と同様に、特殊な環境下、例えば、風洞試験室での着氷風洞試験で着氷を計測する場合でも、計測対象を高精度で効率良くリアルタイムで計測することができる。
【０１２８】
なお、上記各実施例では、２つのカメラで撮影した画像により算出した２つの直線から計測対象点の３次元位置を算出する例を示したが、本発明は、これに限らず、３つ以上のカメラで撮影して計測対象点の３次元位置を算出しても良い。
【産業上の利用可能性】
【０１２９】
本発明は、計測対象の３次元計測に適用されるものであり、特に、特殊な環境下における計測に好適なものであり、例えば、航空機などの着氷風洞試験の際の着氷位置の計測に適用可能である。
【符号の説明】
【０１３０】
１１Ａ  第１カメラ
１１Ｂ  第２カメラ
１２  パターン投影装置
２０  データ処理装置
２１  画像取得部
２１１  計測前画像取得部
２１２  計測時画像取得部
２２Ａ、２２Ｂ  画像解析部
２２１  対象点特定部
２２２  第１直線算出部
２２３  第２直線算出部
２２４  直線設定部
２２５  距離算出部
２２６  最小値算出部
２２７  中点座標算出部
２２８  ３次元位置算出部
４１  格子定規
４１ａ  基準格子
【図１】
2018146363.tif 000003
【図２】
2018146363.tif 000004
【図３】
2018146363.tif 000005
【図４】
2018146363.tif 000006
【図５】
2018146363.tif 000007
【図６】
2018146363.tif 000008
【図７】
2018146363.tif 000009
【図８】
2018146363.tif 000010
【図９】
2018146363.tif 000011
【図１０】
2018146363.tif 000012
【図１１】
2018146363.tif 000013
【図１２】
2018146363.tif 000014
【図１３】
2018146363.tif 000015
【図１４】
2018146363.tif 000016