(57)【要約】【課題】被写体を斜め方向から撮像した一枚の写真画像１０１を用いて被写体の寸法を正確に計測する。【解決手段】姿勢基準軸設定部１２０は、写真画像１０１に映っている被写体の姿勢を表す二次元の姿勢基準軸データ１２９を生成する。基準軸ベクトル推定部１４０は、姿勢基準軸データ１２９と撮影パラメータ１３９とに基づいて、被写体の姿勢を表す三次元の基準軸ベクトル１４９を算出する。リファレンス設定部１５０は、基準軸ベクトル１４９を用いて、被写体のリファレンス線分を表すリファレンス線分ベクトル１５９ａを算出する。寸法推定部１６０は、基準軸ベクトル１４９を用いて、被写体の計測線分を表す計測線分ベクトルを算出する。寸法推定部１６０は、既知であるリファレンス線分長１５９ｂと、リファレンス線分ベクトル１５９ａと、計測線分ベクトルとを用いて、計測線分の長さを算出する。【選択図】図１
2014109440.tif 000002
【特許請求の範囲】
【請求項１】
対象物を斜め方向から撮像した対象画像に映っている前記対象物の幅方向と高さ方向と奥行き方向とに沿う３つの基準線分を前記対象画像の二次元の座標値を用いて表す基準線分情報を生成する基準線分情報生成部と、
前記基準線分情報生成部によって生成された前記基準線分情報と、前記対象画像の撮像パラメータとに基づいて、前記対象物の前記幅方向と前記高さ方向と前記奥行き方向とを三次元のベクトル値を用いて表す３つの基準軸ベクトルを算出する基準軸ベクトル算出部と、
前記基準軸ベクトル算出部によって算出された前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルを用いて、前記対象画像に映っている前記対象物のうち長さが既知の線分であるリファレンス線分のベクトル値を、リファレンス線分ベクトルとして算出するリファレンス線分ベクトル算出部と、
前記基準軸ベクトル算出部によって算出された前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルを用いて、前記対象画像に映っている前記対象物のうち長さが未知の線分である計測線分のベクトル値を、計測線分ベクトルとして算出する計測線分ベクトル算出部と、
前記リファレンス線分の既知の長さと、前記リファレンス線分ベクトル算出部によって算出された前記リファレンス線分ベクトルと、前記計測線分ベクトル算出部によって算出された前記計測線分ベクトルとを用いて、前記計測線分の長さを算出する計測線分長算出部と
を備えることを特徴とする寸法推定装置。
【請求項２】
前記リファレンス線分ベクトル算出部は、
前記対象画像が撮像された位置を示す座標点Ｐｃから前記リファレンス線分の一方の端点Ｐｓへの向きを有する単位ベクトルＶｓを算出し、
前記座標点Ｐｃから前記リファレンス線分の他方の端点Ｐｅへの向きを有する単位ベクトルＶｅを算出し、
前記単位ベクトルＶｓと、前記単位ベクトルＶｅと、前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルとを用いて、前記リファレンス線分ベクトルを算出する
ことを特徴とする請求項１記載の寸法推定装置。
【請求項３】
前記リファレンス線分ベクトル算出部は、
前記座標点Ｐｃから前記リファレンス線分の前記端点ＰｓへのベクトルＶｓｔｓを得るための媒介変数ｔｓを、前記３つの基準軸ベクトルのうちいずれかの基準軸ベクトルを用いて算出し、
前記座標点Ｐｃから前記リファレンス線分の前記端点ＰｅへのベクトルＶｅｔｅを得るための媒介変数ｔｅを、前記いずれかの基準軸ベクトルを用いて算出し、
前記ベクトルＶｓｔｓと前記ベクトルＶｅｔｅとを用いて、前記リファレンス線分ベクトルを算出する
ことを特徴とする請求項２記載の寸法推定装置。
【請求項４】
前記リファレンス線分ベクトル算出部は、以下の式（１０１）を計算して前記リファレンス線分ベクトルＬｒを算出する、ことを特徴とする請求項３記載の寸法推定装置。
【数１０１】
2014109440.tif 000032
【請求項５】
前記リファレンス線分ベクトル算出部は、前記３つの基準軸ベクトルの原点Ｐ０の座標値を用いると共に基準軸ベクトルＯｚを前記いずれかの基準軸ベクトルとして用いて以下の式（１０２）を計算し、前記媒介変数ｔｓと前記媒介変数ｔｅとを算出する、ことを特徴とする請求項４記載の寸法推定装置。
【数１０２】
2014109440.tif 000033
【請求項６】
前記計測線分ベクトル算出部は、
前記対象画像が撮像された位置を示す座標点Ｐｃから前記計測線分の一方の端点Ｐｍｓへの向きを有する単位ベクトルＶｍｓを算出し、
前記座標点Ｐｃから前記リファレンス線分の他方の端点Ｐｍｅへの向きを有する単位ベクトルＶｍｅを算出し、
前記単位ベクトルＶｍｓと、前記単位ベクトルＶｍｅと、前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルとを用いて、前記計測線分ベクトルを算出する
ことを特徴とする請求項１から請求項５のいずれかに記載の寸法推定装置。
【請求項７】
前記計測線分ベクトル算出部は、
前記座標点Ｐｃから前記計測線分の前記端点ＰｍｓへのベクトルＶｍｓｔｍｓを得るための媒介変数ｔｍｓを、前記３つの基準軸ベクトルのうちいずれかの基準軸ベクトルを用いて算出し、
前記座標点Ｐｃから前記計測線分の前記端点ＰｍｅへのベクトルＶｍｅｔｍｅを得るための媒介変数ｔｍｅを、前記いずれかの基準軸ベクトルを用いて算出し、
前記ベクトルＶｍｓｔｍｓと前記ベクトルＶｍｅｔｍｅとを用いて、前記リファレンス線分ベクトルを算出する
ことを特徴とする請求項６記載の寸法推定装置。
【請求項８】
前記計測線分ベクトル算出部は、以下の式（１０３）を計算して前記計測線分ベクトルＬｍを算出する、ことを特徴とする請求項７記載の寸法推定装置。
【数１０３】
2014109440.tif 000034
【請求項９】
前記計測線分ベクトル算出部は、前記３つの基準軸ベクトルの原点Ｐ０の座標値を用いると共に基準軸ベクトルＯｘを前記いずれかの基準軸ベクトルとして用いて以下の式（１０４）を計算し、前記媒介変数ｔｍｓと前記媒介変数ｔｍｅとを算出する、ことを特徴とする請求項８記載の寸法推定装置。
【数１０４】
2014109440.tif 000035
【請求項１０】
前記計測線分長算出部は、前記リファレンス線分ベクトルに対する前記計測線分ベクトルの大きさの比率に前記リファレンス線分の既知の長さを乗じた値を前記計測線分の長さとして算出する
ことを特徴とする請求項１から請求項９のいずれかに記載の寸法推定装置。
【請求項１１】
前記基準軸ベクトル算出部は、
前記基準線分情報と前記撮像パラメータとに基づいて、前記対象物の前記幅方向と前記高さ方向と前記奥行き方向とを表す３つの単位ベクトルであって、ＦＰＭ（Ｆｉｒｓｔ  Ｐｅｒｓｐｅｃｔｉｖｅ  Ｍｏｖｉｎｇ）座標系またはＳＰＭ（Ｓｅｃｏｎｄ  Ｐｅｒｓｐｅｃｔｉｖｅ  Ｍｏｖｉｎｇ）座標系の３つの単位ベクトルを算出し、
ＦＰＭ座標系またはＳＰＭ座標系の前記３つの単位ベクトルを用いて、ワールド座標系の３つのベクトルを前記３つの基準軸ベクトルとして算出する
ことを特徴とする請求項１から請求項１０のいずれかに記載の寸法推定装置。
【請求項１２】
対象物を斜め方向から撮像した対象画像に映っている前記対象物の幅方向と高さ方向と奥行き方向とに沿う３つの基準線分を前記対象画像の二次元の座標値を用いて表す基準線分情報を生成する基準線分情報生成処理と、
前記基準線分情報生成処理によって生成された前記基準線分情報と、前記対象画像の撮像パラメータとに基づいて、前記対象物の前記幅方向と前記高さ方向と前記奥行き方向とを三次元のベクトル値を用いて表す３つの基準軸ベクトルを算出する基準軸ベクトル算出処理と、
前記基準軸ベクトル算出処理によって算出された前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルを用いて、前記対象画像に映っている前記対象物のうち長さが既知の線分であるリファレンス線分のベクトル値を、リファレンス線分ベクトルとして算出するリファレンス線分ベクトル算出処理と、
前記基準軸ベクトル算出処理によって算出された前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルを用いて、前記対象画像に映っている前記対象物のうち長さが未知の線分である計測線分のベクトル値を、計測線分ベクトルとして算出する計測線分ベクトル算出処理と、
前記リファレンス線分の既知の長さと、前記リファレンス線分ベクトル算出処理によって算出された前記リファレンス線分ベクトルと、前記計測線分ベクトル算出処理によって算出された前記計測線分ベクトルとを用いて、前記計測線分の長さを算出する計測線分長算出処理と
をコンピュータに実行させるための寸法推定プログラム。
【請求項１３】
前記リファレンス線分ベクトル算出処理は、
前記対象画像が撮像された位置を示す座標点Ｐｃから前記リファレンス線分の一方の端点Ｐｓへの向きを有する単位ベクトルＶｓを算出する処理と、
前記座標点Ｐｃから前記リファレンス線分の他方の端点Ｐｅへの向きを有する単位ベクトルＶｅを算出する処理と、
前記単位ベクトルＶｓと、前記単位ベクトルＶｅと、前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルとを用いて、前記リファレンス線分ベクトルを算出する処理とを含む
ことを特徴とする請求項１２記載の寸法推定プログラム。
【請求項１４】
前記リファレンス線分ベクトル算出処理は、
前記座標点Ｐｃから前記リファレンス線分の前記端点ＰｓへのベクトルＶｓｔｓを得るための媒介変数ｔｓを、前記３つの基準軸ベクトルのうちいずれかの基準軸ベクトルを用いて算出する処理と、
前記座標点Ｐｃから前記リファレンス線分の前記端点ＰｅへのベクトルＶｅｔｅを得るための媒介変数ｔｅを、前記いずれかの基準軸ベクトルを用いて算出する処理と、
前記ベクトルＶｓｔｓと前記ベクトルＶｅｔｅとを用いて、前記リファレンス線分ベクトルを算出する処理とを含む
ことを特徴とする請求項１３記載の寸法推定プログラム。
【請求項１５】
前記リファレンス線分ベクトル算出処理は、以下の式（１０１）を計算して前記リファレンス線分ベクトルＬｒを算出する処理を含む、ことを特徴とする請求項１４記載の寸法推定プログラム。
【数１０５】
2014109440.tif 000036
【請求項１６】
前記リファレンス線分ベクトル算出処理は、前記３つの基準軸ベクトルの原点Ｐ０の座標値を用いると共に基準軸ベクトルＯｚを前記いずれかの基準軸ベクトルとして用いて以下の式（１０２）を計算し、前記媒介変数ｔｓと前記媒介変数ｔｅとを算出する処理を含む、ことを特徴とする請求項１５記載の寸法推定プログラム。
【数１０６】
2014109440.tif 000037
【請求項１７】
前記計測線分ベクトル算出処理は、
前記対象画像が撮像された位置を示す座標点Ｐｃから前記計測線分の一方の端点Ｐｍｓへの向きを有する単位ベクトルＶｍｓを算出する処理と、
前記座標点Ｐｃから前記リファレンス線分の他方の端点Ｐｍｅへの向きを有する単位ベクトルＶｍｅを算出する処理と、
前記単位ベクトルＶｍｓと、前記単位ベクトルＶｍｅと、前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルとを用いて、前記計測線分ベクトルを算出する処理とを含む
ことを特徴とする請求項１２から請求項１６のいずれかに記載の寸法推定プログラム。
【請求項１８】
前記計測線分ベクトル算出処理は、
前記座標点Ｐｃから前記計測線分の前記端点ＰｍｓへのベクトルＶｍｓｔｍｓを得るための媒介変数ｔｍｓを、前記３つの基準軸ベクトルのうちいずれかの基準軸ベクトルを用いて算出する処理と、
前記座標点Ｐｃから前記計測線分の前記端点ＰｍｅへのベクトルＶｍｅｔｍｅを得るための媒介変数ｔｍｅを、前記いずれかの基準軸ベクトルを用いて算出する処理と、
前記ベクトルＶｍｓｔｍｓと前記ベクトルＶｍｅｔｍｅとを用いて、前記リファレンス線分ベクトルを算出する処理とを含む
ことを特徴とする請求項１７記載の寸法推定プログラム。
【請求項１９】
前記計測線分ベクトル算出処理は、以下の式（１０３）を計算して前記計測線分ベクトルＬｍを算出する処理を含む、ことを特徴とする請求項１８記載の寸法推定プログラム。
【数１０７】
2014109440.tif 000038
【請求項２０】
前記計測線分ベクトル算出処理は、前記３つの基準軸ベクトルの原点Ｐ０の座標値を用いると共に基準軸ベクトルＯｘを前記いずれかの基準軸ベクトルとして用いて以下の式（１０４）を計算し、前記媒介変数ｔｍｓと前記媒介変数ｔｍｅとを算出する処理を含む、ことを特徴とする請求項１９記載の寸法推定プログラム。
【数１０８】
2014109440.tif 000039
【請求項２１】
前記計測線分長算出処理は、前記リファレンス線分ベクトルに対する前記計測線分ベクトルの大きさの比率に前記リファレンス線分の既知の長さを乗じた値を前記計測線分の長さとして算出する処理を含む
ことを特徴とする請求項１２から請求項２０のいずれかに記載の寸法推定プログラム。
【請求項２２】
前記基準軸ベクトル算出処理は、
前記基準線分情報と前記撮像パラメータとに基づいて、前記対象物の前記幅方向と前記高さ方向と前記奥行き方向とを表す３つの単位ベクトルであって、ＦＰＭ（Ｆｉｒｓｔ  Ｐｅｒｓｐｅｃｔｉｖｅ  Ｍｏｖｉｎｇ）座標系またはＳＰＭ（Ｓｅｃｏｎｄ  Ｐｅｒｓｐｅｃｔｉｖｅ  Ｍｏｖｉｎｇ）座標系の３つの単位ベクトルを算出する処理と、
ＦＰＭ座標系またはＳＰＭ座標系の前記３つの単位ベクトルを用いて、ワールド座標系の３つのベクトルを前記３つの基準軸ベクトルとして算出する処理とを含む
ことを特徴とする請求項１２から請求項２１のいずれかに記載の寸法推定プログラム。
【請求項２３】
基準線分情報生成部が、対象物を斜め方向から撮像した対象画像に映っている前記対象物の幅方向と高さ方向と奥行き方向とに沿う３つの基準線分を前記対象画像の二次元の座標値を用いて表す基準線分情報を生成し、
基準軸ベクトル算出部が、前記基準線分情報生成部によって生成された前記基準線分情報と、前記対象画像の撮像パラメータとに基づいて、前記対象物の前記幅方向と前記高さ方向と前記奥行き方向とを三次元のベクトル値を用いて表す３つの基準軸ベクトルを算出し、
リファレンス線分ベクトル算出部が、前記基準軸ベクトル算出部によって算出された前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルを用いて、前記対象画像に映っている前記対象物のうち長さが既知の線分であるリファレンス線分のベクトル値を、リファレンス線分ベクトルとして算出し、
計測線分ベクトル算出部が、前記基準軸ベクトル算出部によって算出された前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルを用いて、前記対象画像に映っている前記対象物のうち長さが未知の線分である計測線分のベクトル値を、計測線分ベクトルとして算出し、
計測線分長算出部が、前記リファレンス線分の既知の長さと、前記リファレンス線分ベクトル算出部によって算出された前記リファレンス線分ベクトルと、前記計測線分ベクトル算出部によって算出された前記計測線分ベクトルとを用いて、前記計測線分の長さを算出する
ことを特徴とする寸法推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、例えば、写真に映っている被写体の寸法を計測するための寸法推定装置、寸法推定プログラムおよび寸法推定方法に関するものである。
【背景技術】
【０００２】
写真に映った被写体の寸法を計測する方法としてステレオ計測が知られている。ステレオ計測とは、被写体を異なる方向から撮影し、撮影時の視差を利用して三角測量の原理で被写体の寸法を計測する方法である。
しかし、海上の不審船などを監視する監視業務のように監視対象を撮影する回数または角度が限定される場合、適当な視差を確保しながら監視対象を異なる方向から撮影することができない。そのため、ステレオ計測によって監視対象の正確な寸法を計測することは困難である。
【０００３】
また、大きさが既知である比較部分を含んだ被写体が写真に映っている場合、被写体の計測部分と比較部分との大きさを比較することによって、被写体の計測部分の寸法を計測することができる。
但し、被写体が斜め方向から撮影されている場合、被写体の手前側の部分と奥側の部分とで写真に映る大きさが異なるため、被写体の計測部分の寸法を正確に計測することができない。
【０００４】
図３６は、従来の寸法計測方法の概要図である。
例えば、被写体の円部分の直径「２メートル」が既知であるものとする（図３６の（ａ））。ここで、被写体は正面から撮像されているものとする（図３６の（ｂ））。
この場合、被写体の円部分の画素数「２０」と被写体の高さ方向の画素数「４０」との比率「２（＝４０／２０）」と、被写体の円部分の直径「２メートル」とを用いて、被写体の高さ「４メートル（＝２×２メートル）」を算出することができる（図３６の（ｃ）（ｄ））。
同様に、被写体の円部分の画素数「２０」と被写体の幅方向の画素数「１２０」との比率「６（＝１２０／２０）」と、被写体の円部分の直径「２メートル」とを用いて、被写体の幅「１２メートル（＝６×２メートル）」を算出することができる（図３６の（ｃ）（ｄ））。
但し、被写体が斜め方向から撮像されている場合（図３６の（ｅ））、被写体の手前側の部分ｄ２と奥側の部分ｄ１とで写真に映る大きさが異なるため、被写体の高さまたは幅を正確に計測することができない。
【０００５】
非特許文献１には、斜め方向から撮影された被写体が映っている１枚の写真から被写体が映っている向きを表す三次元の姿勢ベクトルを算出する方法が開示されている。
しかし、被写体が映っている向きを表す三次元の姿勢ベクトルに基づいて被写体の寸法を計測する方法は開示も示唆もされていない。
【０００６】
特許文献１には、視覚に基づく三次元位置および姿勢の認識方法が開示されている。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開平７−９８２１４号公報
【非特許文献】
【０００８】
【非特許文献１】Ｔ．Ｓｈａｋｕｎａｇａ，Ｈ．Ｋａｎｅｋｏ，  ‘Ｐｅｒｓｐｅｃｔｉｖｅ  ａｎｇｌｅ  ｔｒａｎｓｆｏｒｍ：  Ｐｒｉｎｃｉｐｌｅ  ｆｏｒ  ｓｈａｐｅ  ｆｒｏｍ  ａｎｇｌｅｓ’，  Ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｊｏｕｒｎａｌ  ｏｆ  Ｃｏｍｐｕｔｅｒ  Ｖｉｓｉｏｎ，  ｖｏｌ．３，  ｎｏ．３，  ｐｐ．２３９−２５４，  １９８９．
【発明の概要】
【発明が解決しようとする課題】
【０００９】
本発明は、例えば、対象物を斜め方向から撮像した対象画像を用いて対象物の寸法を正確に計測できるようにすることを目的とする。
【課題を解決するための手段】
【００１０】
本発明の寸法推定装置は、
対象物を斜め方向から撮像した対象画像に映っている前記対象物の幅方向と高さ方向と奥行き方向とに沿う３つの基準線分を前記対象画像の二次元の座標値を用いて表す基準線分情報を生成する基準線分情報生成部と、
前記基準線分情報生成部によって生成された前記基準線分情報と、前記対象画像の撮像パラメータとに基づいて、前記対象物の前記幅方向と前記高さ方向と前記奥行き方向とを三次元のベクトル値を用いて表す３つの基準軸ベクトルを算出する基準軸ベクトル算出部と、
前記基準軸ベクトル算出部によって算出された前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルを用いて、前記対象画像に映っている前記対象物のうち長さが既知の線分であるリファレンス線分のベクトル値を、リファレンス線分ベクトルとして算出するリファレンス線分ベクトル算出部と、
前記基準軸ベクトル算出部によって算出された前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルを用いて、前記対象画像に映っている前記対象物のうち長さが未知の線分である計測線分のベクトル値を、計測線分ベクトルとして算出する計測線分ベクトル算出部と、
前記リファレンス線分の既知の長さと、前記リファレンス線分ベクトル算出部によって算出された前記リファレンス線分ベクトルと、前記計測線分ベクトル算出部によって算出された前記計測線分ベクトルとを用いて、前記計測線分の長さを算出する計測線分長算出部とを備える。
【００１１】
前記リファレンス線分ベクトル算出部は、
前記対象画像が撮像された位置を示す座標点Ｐｃから前記リファレンス線分の一方の端点Ｐｓへの向きを有する単位ベクトルＶｓを算出し、
前記座標点Ｐｃから前記リファレンス線分の他方の端点Ｐｅへの向きを有する単位ベクトルＶｅを算出し、
前記単位ベクトルＶｓと、前記単位ベクトルＶｅと、前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルとを用いて、前記リファレンス線分ベクトルを算出する。
【００１２】
前記リファレンス線分ベクトル算出部は、
前記座標点Ｐｃから前記リファレンス線分の前記端点ＰｓへのベクトルＶｓｔｓを得るための媒介変数ｔｓを、前記３つの基準軸ベクトルのうちいずれかの基準軸ベクトルを用いて算出し、
前記座標点Ｐｃから前記リファレンス線分の前記端点ＰｅへのベクトルＶｅｔｅを得るための媒介変数ｔｅを、前記いずれかの基準軸ベクトルを用いて算出し、
前記ベクトルＶｓｔｓと前記ベクトルＶｅｔｅとを用いて、前記リファレンス線分ベクトルを算出する。
【００１３】
前記リファレンス線分ベクトル算出部は、以下の式（１０１）を計算して前記リファレンス線分ベクトルＬｒを算出する。
【００１４】
【数１】
2014109440.tif 000003
【００１５】
前記リファレンス線分ベクトル算出部は、前記３つの基準軸ベクトルの原点Ｐ０の座標値を用いると共に基準軸ベクトルＯｚを前記いずれかの基準軸ベクトルとして用いて以下の式（１０２）を計算し、前記媒介変数ｔｓと前記媒介変数ｔｅとを算出する。
【００１６】
【数２】
2014109440.tif 000004
【００１７】
前記計測線分ベクトル算出部は、
前記対象画像が撮像された位置を示す座標点Ｐｃから前記計測線分の一方の端点Ｐｍｓへの向きを有する単位ベクトルＶｍｓを算出し、
前記座標点Ｐｃから前記リファレンス線分の他方の端点Ｐｍｅへの向きを有する単位ベクトルＶｍｅを算出し、
前記単位ベクトルＶｍｓと、前記単位ベクトルＶｍｅと、前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルとを用いて、前記計測線分ベクトルを算出する。
【００１８】
前記計測線分ベクトル算出部は、
前記座標点Ｐｃから前記計測線分の前記端点ＰｍｓへのベクトルＶｍｓｔｍｓを得るための媒介変数ｔｍｓを、前記３つの基準軸ベクトルのうちいずれかの基準軸ベクトルを用いて算出し、
前記座標点Ｐｃから前記計測線分の前記端点ＰｍｅへのベクトルＶｍｅｔｍｅを得るための媒介変数ｔｍｅを、前記いずれかの基準軸ベクトルを用いて算出し、
前記ベクトルＶｍｓｔｍｓと前記ベクトルＶｍｅｔｍｅとを用いて、前記リファレンス線分ベクトルを算出する。
【００１９】
前記計測線分ベクトル算出部は、以下の式（１０３）を計算して前記計測線分ベクトルＬｍを算出する。
【００２０】
【数３】
2014109440.tif 000005
【００２１】
前記計測線分ベクトル算出部は、前記３つの基準軸ベクトルの原点Ｐ０の座標値を用いると共に基準軸ベクトルＯｘを前記いずれかの基準軸ベクトルとして用いて以下の式（１０４）を計算し、前記媒介変数ｔｍｓと前記媒介変数ｔｍｅとを算出する。
【００２２】
【数４】
2014109440.tif 000006
【００２３】
前記計測線分長算出部は、前記リファレンス線分ベクトルに対する前記計測線分ベクトルの大きさの比率に前記リファレンス線分の既知の長さを乗じた値を前記計測線分の長さとして算出する。
【００２４】
前記基準軸ベクトル算出部は、
前記基準線分情報と前記撮像パラメータとに基づいて、前記対象物の前記幅方向と前記高さ方向と前記奥行き方向とを表す３つの単位ベクトルであって、ＦＰＭ（Ｆｉｒｓｔ  Ｐｅｒｓｐｅｃｔｉｖｅ  Ｍｏｖｉｎｇ）座標系またはＳＰＭ（Ｓｅｃｏｎｄ  Ｐｅｒｓｐｅｃｔｉｖｅ  Ｍｏｖｉｎｇ）座標系の３つの単位ベクトルを算出し、
ＦＰＭ座標系またはＳＰＭ座標系の前記３つの単位ベクトルを用いて、ワールド座標系の３つのベクトルを前記３つの基準軸ベクトルとして算出する。
【００２５】
本発明の寸法推定プログラムは、
対象物を斜め方向から撮像した対象画像に映っている前記対象物の幅方向と高さ方向と奥行き方向とに沿う３つの基準線分を前記対象画像の二次元の座標値を用いて表す基準線分情報を生成する基準線分情報生成処理と、
前記基準線分情報生成処理によって生成された前記基準線分情報と、前記対象画像の撮像パラメータとに基づいて、前記対象物の前記幅方向と前記高さ方向と前記奥行き方向とを三次元のベクトル値を用いて表す３つの基準軸ベクトルを算出する基準軸ベクトル算出処理と、
前記基準軸ベクトル算出処理によって算出された前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルを用いて、前記対象画像に映っている前記対象物のうち長さが既知の線分であるリファレンス線分のベクトル値を、リファレンス線分ベクトルとして算出するリファレンス線分ベクトル算出処理と、
前記基準軸ベクトル算出処理によって算出された前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルを用いて、前記対象画像に映っている前記対象物のうち長さが未知の線分である計測線分のベクトル値を、計測線分ベクトルとして算出する計測線分ベクトル算出処理と、
前記リファレンス線分の既知の長さと、前記リファレンス線分ベクトル算出処理によって算出された前記リファレンス線分ベクトルと、前記計測線分ベクトル算出処理によって算出された前記計測線分ベクトルとを用いて、前記計測線分の長さを算出する計測線分長算出処理と
をコンピュータに実行させる。
【００２６】
前記リファレンス線分ベクトル算出処理は、
前記対象画像が撮像された位置を示す座標点Ｐｃから前記リファレンス線分の一方の端点Ｐｓへの向きを有する単位ベクトルＶｓを算出する処理と、
前記座標点Ｐｃから前記リファレンス線分の他方の端点Ｐｅへの向きを有する単位ベクトルＶｅを算出する処理と、
前記単位ベクトルＶｓと、前記単位ベクトルＶｅと、前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルとを用いて、前記リファレンス線分ベクトルを算出する処理とを含む。
【００２７】
前記リファレンス線分ベクトル算出処理は、
前記座標点Ｐｃから前記リファレンス線分の前記端点ＰｓへのベクトルＶｓｔｓを得るための媒介変数ｔｓを、前記３つの基準軸ベクトルのうちいずれかの基準軸ベクトルを用いて算出する処理と、
前記座標点Ｐｃから前記リファレンス線分の前記端点ＰｅへのベクトルＶｅｔｅを得るための媒介変数ｔｅを、前記いずれかの基準軸ベクトルを用いて算出する処理と、
前記ベクトルＶｓｔｓと前記ベクトルＶｅｔｅとを用いて、前記リファレンス線分ベクトルを算出する処理とを含む。
【００２８】
前記リファレンス線分ベクトル算出処理は、以下の式（１０１）を計算して前記リファレンス線分ベクトルＬｒを算出する処理を含む。
【００２９】
【数５】
2014109440.tif 000007
【００３０】
前記リファレンス線分ベクトル算出処理は、前記３つの基準軸ベクトルの原点Ｐ０の座標値を用いると共に基準軸ベクトルＯｚを前記いずれかの基準軸ベクトルとして用いて以下の式（１０２）を計算し、前記媒介変数ｔｓと前記媒介変数ｔｅとを算出する処理を含む。
【００３１】
【数６】
2014109440.tif 000008
【００３２】
前記計測線分ベクトル算出処理は、
前記対象画像が撮像された位置を示す座標点Ｐｃから前記計測線分の一方の端点Ｐｍｓへの向きを有する単位ベクトルＶｍｓを算出する処理と、
前記座標点Ｐｃから前記リファレンス線分の他方の端点Ｐｍｅへの向きを有する単位ベクトルＶｍｅを算出する処理と、
前記単位ベクトルＶｍｓと、前記単位ベクトルＶｍｅと、前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルとを用いて、前記計測線分ベクトルを算出する処理とを含む。
【００３３】
前記計測線分ベクトル算出処理は、
前記座標点Ｐｃから前記計測線分の前記端点ＰｍｓへのベクトルＶｍｓｔｍｓを得るための媒介変数ｔｍｓを、前記３つの基準軸ベクトルのうちいずれかの基準軸ベクトルを用いて算出する処理と、
前記座標点Ｐｃから前記計測線分の前記端点ＰｍｅへのベクトルＶｍｅｔｍｅを得るための媒介変数ｔｍｅを、前記いずれかの基準軸ベクトルを用いて算出する処理と、
前記ベクトルＶｍｓｔｍｓと前記ベクトルＶｍｅｔｍｅとを用いて、前記リファレンス線分ベクトルを算出する処理とを含む。
【００３４】
前記計測線分ベクトル算出処理は、以下の式（１０３）を計算して前記計測線分ベクトルＬｍを算出する処理を含む。
【００３５】
【数７】
2014109440.tif 000009
【００３６】
前記計測線分ベクトル算出処理は、前記３つの基準軸ベクトルの原点Ｐ０の座標値を用いると共に基準軸ベクトルＯｘを前記いずれかの基準軸ベクトルとして用いて以下の式（１０４）を計算し、前記媒介変数ｔｍｓと前記媒介変数ｔｍｅとを算出する処理を含む。
【００３７】
【数８】
2014109440.tif 000010
【００３８】
前記計測線分長算出処理は、前記リファレンス線分ベクトルに対する前記計測線分ベクトルの大きさの比率に前記リファレンス線分の既知の長さを乗じた値を前記計測線分の長さとして算出する処理を含む。
【００３９】
前記基準軸ベクトル算出処理は、
前記基準線分情報と前記撮像パラメータとに基づいて、前記対象物の前記幅方向と前記高さ方向と前記奥行き方向とを表す３つの単位ベクトルであって、ＦＰＭ（Ｆｉｒｓｔ  Ｐｅｒｓｐｅｃｔｉｖｅ  Ｍｏｖｉｎｇ）座標系またはＳＰＭ（Ｓｅｃｏｎｄ  Ｐｅｒｓｐｅｃｔｉｖｅ  Ｍｏｖｉｎｇ）座標系の３つの単位ベクトルを算出する処理と、
ＦＰＭ座標系またはＳＰＭ座標系の前記３つの単位ベクトルを用いて、ワールド座標系の３つのベクトルを前記３つの基準軸ベクトルとして算出する処理とを含む。
【００４０】
本発明の寸法推定方法は、
基準線分情報生成部が、対象物を斜め方向から撮像した対象画像に映っている前記対象物の幅方向と高さ方向と奥行き方向とに沿う３つの基準線分を前記対象画像の二次元の座標値を用いて表す基準線分情報を生成し、
基準軸ベクトル算出部が、前記基準線分情報生成部によって生成された前記基準線分情報と、前記対象画像の撮像パラメータとに基づいて、前記対象物の前記幅方向と前記高さ方向と前記奥行き方向とを三次元のベクトル値を用いて表す３つの基準軸ベクトルを算出し、
リファレンス線分ベクトル算出部が、前記基準軸ベクトル算出部によって算出された前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルを用いて、前記対象画像に映っている前記対象物のうち長さが既知の線分であるリファレンス線分のベクトル値を、リファレンス線分ベクトルとして算出し、
計測線分ベクトル算出部が、前記基準軸ベクトル算出部によって算出された前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルを用いて、前記対象画像に映っている前記対象物のうち長さが未知の線分である計測線分のベクトル値を、計測線分ベクトルとして算出し、
計測線分長算出部が、前記リファレンス線分の既知の長さと、前記リファレンス線分ベクトル算出部によって算出された前記リファレンス線分ベクトルと、前記計測線分ベクトル算出部によって算出された前記計測線分ベクトルとを用いて、前記計測線分の長さを算出する。
【発明の効果】
【００４１】
本発明によれば、例えば、対象物を斜め方向から撮像した対象画像を用いて対象物の寸法を正確に計測することができる。
【図面の簡単な説明】
【００４２】
【図１】実施の形態１における画像寸法推定装置１００の機能構成図である。
【図２】実施の形態１における画像寸法推定方法を示すフローチャートである。
【図３】実施の形態１における写真画像１０１の一例を示す図である。
【図４】実施の形態１における姿勢基準軸［ｌｘ，ｌｙ，ｌｚ］の一例を示す図である。
【図５】実施の形態１における姿勢基準軸［ｌｘ，ｌｙ，ｌｚ］の一例を示す図である。
【図６】実施の形態１における撮影パラメータ１３９の関係図である。
【図７】実施の形態１における撮影パラメータ設定処理（Ｓ１３０）の一例を示すフローチャートである。
【図８】実施の形態１におけるリファレンス設定部１５０の構成図である。
【図９】実施の形態１におけるリファレンス設定処理（Ｓ１５０）を示すフローチャートである。
【図１０】実施の形態１における基準軸ベクトル投影処理（Ｓ１５１）の概要図である。
【図１１】実施の形態１における基準軸ベクトル移動処理（Ｓ１５３）の概要図である。
【図１２】実施の形態１におけるリファレンス線分設定処理（Ｓ１５４）の概要図である。
【図１３】実施の形態１におけるリファレンス線分設定処理（Ｓ１５５）の概要図である。
【図１４】実施の形態１におけるリファレンス線分ベクトル算出処理（Ｓ１５６）の概要図である。
【図１５】実施の形態１における寸法推定部１６０の構成図である。
【図１６】実施の形態１における寸法推定処理（Ｓ１６０）を示すフローチャートである。
【図１７】実施の形態１におけるＳ１６１からＳ１６５の概要図である。
【図１８】実施の形態１におけるＳ１６６の概要図である。
【図１９】実施の形態１における基準軸ベクトル推定処理（Ｓ１４０）を示すフローチャートである。
【図２０】実施の形態１における基準軸ベクトル算出処理［ＦＰＭ座標系］（Ｓ１４２）を示すフローチャートである。
【図２１】実施の形態１における端点単位ベクトル算出処理（Ｓ１４２Ａ）の概要図である。
【図２２】実施の形態１における交点単位ベクトル算出処理（Ｓ１４２Ｂ）の概要図である。
【図２３】実施の形態１におけるＦＰＭ座標面算出処理（Ｓ１４２Ｃ）の概要図である。
【図２４】実施の形態１におけるＦＰＭ座標面ベクトル算出処理（Ｓ１４２Ｄ）の概要図である。
【図２５】実施の形態１におけるＦＰＭ座標系の単位ベクトル算出処理（Ｓ１４２Ｅ）の概要図である。
【図２６】実施の形態１における被写体座標系の基準軸ベクトルの関係図である。
【図２７】実施の形態１における基準軸ベクトル算出処理［ＳＰＭ座標系］（Ｓ１４３）を示すフローチャートである。
【図２８】実施の形態１における交点単位ベクトル算出処理（Ｓ１４３Ａ）の概要図である。
【図２９】実施の形態１におけるＳＰＭ座標面算出処理（Ｓ１４３Ｂ）の概要図である。
【図３０】実施の形態１における姿勢基準軸投影処理（Ｓ１４３Ｃ）の概要図である。
【図３１】実施の形態１におけるＳＰＭ座標系の単位ベクトル算出処理（Ｓ１４３Ｄ）の概要図である。
【図３２】実施の形態１におけるＳＰＭ座標面内の姿勢基準軸の関係図である。
【図３３】実施の形態１における基準原点算出処理（Ｓ１４３Ｆ）の概要図である。
【図３４】実施の形態１における画像寸法推定装置１００のハードウェア資源の一例を示す図である。
【図３５】実施の形態１における被写体１０２の三次元モデルを表す図である。
【図３６】従来の寸法計測方法の概要図である。
【発明を実施するための形態】
【００４３】
実施の形態１．
被写体を斜め方向から撮像した一枚の画像を用いて被写体の寸法を正確に計測する形態について説明する。
【００４４】
図１は、実施の形態１における画像寸法推定装置１００の機能構成図である。
実施の形態１における画像寸法推定装置１００の機能構成について、図１に基づいて説明する。
【００４５】
画像寸法推定装置１００（寸法推定装置の一例）は、写真画像表示部１１０、姿勢基準軸設定部１２０（基準線分情報生成部の一例）、撮影パラメータ設定部１３０、基準軸ベクトル推定部１４０（基準軸ベクトル算出部の一例）、リファレンス設定部１５０（リファレンス線分ベクトル算出部の一例）、寸法推定部１６０（計測線分ベクトル算出部、計測線分長算出部の一例）および推定結果出力部１７０を備える。
また、画像寸法推定装置１００は、写真画像データベース１８０、装置記憶部１９０および入力機器１０９を備える。
【００４６】
写真画像表示部１１０は、入力機器１０９を用いて入力される写真画像１０１をディスプレイに表示する。
姿勢基準軸設定部１２０は、写真画像１０１に映っている被写体の姿勢基準軸を写真画像１０１の二次元の座標値を用いて示す姿勢基準軸データ１２９を生成する。
撮影パラメータ設定部１３０は、写真画像１０１を撮影した際にカメラに設定された撮影パラメータ１３９を取得する。
基準軸ベクトル推定部１４０は、姿勢基準軸データ１２９と撮影パラメータ１３９とを用いて、写真画像１０１に映っている被写体の姿勢を表す三次元の基準軸ベクトル１４９を算出する。
リファレンス設定部１５０は、基準軸ベクトル１４９を用いて、写真画像１０１に映っている被写体のうち長さが既知の線分であるリファレンス線分を表す三次元のリファレンス線分ベクトル１５９ａを算出する。また、リファレンス設定部１５０は、リファレンス線分の既知の長さとしてリファレンス線分長１５９ｂを取得する。
寸法推定部１６０は、基準軸ベクトル１４９とリファレンス線分ベクトル１５９ａとリファレンス線分長１５９ｂとを用いて、写真画像１０１に映っている被写体のうち長さが未知の線分である計測線分の長さを示す計測線分長１６９を算出する。
推定結果出力部１７０は、計測線分長１６９を出力する。
【００４７】
写真画像データベース１８０は、写真画像１０１と計測線分長１６９とを対応付けて記憶するデータベースである。
【００４８】
装置記憶部１９０は、画像寸法推定装置１００で使用するデータを記憶する。
例えば、写真画像１０１、姿勢基準軸データ１２９、撮影パラメータ１３９、基準軸ベクトル１４９、リファレンス線分ベクトル１５９ａ、リファレンス線分長１５９ｂおよび計測線分長１６９は、装置記憶部１９０に記憶されるデータの一例である。
【００４９】
図２は、実施の形態１における画像寸法推定方法を示すフローチャートである。
実施の形態１における画像寸法推定方法について、図２に基づいて説明する。
【００５０】
Ｓ１１０において、利用者は、キーボードまたはマウスなどの入力機器１０９を用いて写真画像１０１を指定する。
写真画像表示部１１０は、利用者の指定に従って、ＣＤドライブまたはスキャナなどの入力機器１０９から写真画像１０１を取得する。例えば、写真画像表示部１１０は、ＣＤドライブに挿入されたＣＤ（コンパクトディスク）から写真画像１０１を読み出す。または、写真画像表示部１１０は、写真画像１０１を読み取ったスキャナから写真画像１０１を受け取る。但し、写真画像表示部１１０は、ＣＤドライブまたはスキャナ以外の入力機器１０９から写真画像１０１を取得しても構わない。
写真画像表示部１１０は、取得した写真画像１０１をディスプレイに表示する。
図３は、実施の形態１における写真画像１０１の一例を示す図である。
図３に示すように、写真画像１０１は、被写体１０２を斜めの方向から撮像して得られた画像である。
Ｓ１１０の後、処理はＳ１２０に進む。
【００５１】
Ｓ１２０において、利用者は、入力機器１０９を用いて、写真画像１０１に映っている被写体１０２の姿勢を表す３つの線分を被写体１０２の姿勢基準軸として指定する。
被写体１０２の姿勢を表す３つの線分とは、被写体１０２の幅方向を表す線分と、被写体１０２の奥行き方向を表す線分と、被写体１０２の高さ方向を表す線分との３つの線分のことである。
姿勢基準軸設定部１２０は、指定された姿勢基準軸を写真画像１０１の二次元の座標値を用いて特定し、姿勢基準軸を表す姿勢基準軸データ１２９を生成する。写真画像１０１の二次元の座標値とは、写真画像１０１の横方向の座標値および写真画像１０１の縦方向の座標値である。
図４、図５は、実施の形態１における姿勢基準軸［ｌｘ，ｌｙ，ｌｚ］の一例を示す図である。
図４に示すように、利用者は、マウスを操作することによってマウスカーソル１０９ａを移動し、被写体１０２の高さ方向の辺に沿った線分ｌｘと、被写体１０２の奥行き方向の辺に沿った線分ｌｙと、被写体１０２の幅方向の辺に沿った線分ｌｚとを指定する。但し、利用者は、マウス以外の入力機器１０９（例えば、キーボード）を用いて３つの線分ｌｘ、ｌｙ、ｌｚを指定しても構わない。
例えば、姿勢基準軸設定部１２０は、線分ｌｘの両端の二次元の座標値と、線分ｌｙの両端の二次元の座標値と、線分ｌｚの両端の二次元の座標値とを特定し、特定した各座標値を含んだデータを姿勢基準軸データ１２９として生成する。
利用者は、図４に示すように１点で交わっている３つの線分ｌｘ、ｌｙ、ｌｚを指定してもよいし、図５に示すように１点で交わっていない３つの線分ｌｘ、ｌｙ、ｌｚを指定してもよい。
Ｓ１２０の後、処理はＳ１３０に進む。
但し、Ｓ１２０はＳ１３０の後に実行してもよいし、Ｓ１３０と並行して実行しても構わない。
【００５２】
Ｓ１３０において、撮影パラメータ設定部１３０は、写真画像１０１を撮像したカメラの焦点距離、ＣＣＤ（Ｃｈａｒｇｅ  Ｃｏｕｐｌｅｄ  Ｄｅｖｉｃｅ）の縦方向のサイズまたはレンズの縦方向の視野角などを撮影パラメータ１３９として取得する。ＣＣＤはイメージセンサの一例である。
以下、ＣＣＤの縦方向のサイズを縦方向ＣＣＤサイズという。また、レンズの縦方向の視野角を縦方向視野角という。
図６は、実施の形態１における撮影パラメータ１３９の関係図である。焦点距離、縦方向ＣＣＤサイズおよび縦方向視野角は、図６に示すような関係にある。
撮影パラメータ設定部１３０は、焦点距離、縦方向ＣＣＤサイズまたは縦方向視野角に相当する他の撮影パラメータ１３９を取得しても構わない。
【００５３】
図７は、実施の形態１における撮影パラメータ設定処理（Ｓ１３０）の一例を示すフローチャートである。
実施の形態１における撮影パラメータ設定処理（Ｓ１３０）について、図７に基づいて説明する。
【００５４】
Ｓ１３１において、撮影パラメータ設定部１３０は、撮影パラメータ１３９に関する付加情報が写真画像１０１に付加されているか否かを判定する。
例えば、Ｅｘｉｆ（Ｅｘｃｈａｎｇｅａｂｌｅ  ｉｍａｇｅ  ｆｉｌｅ  ｆｏｒｍａｔ）は、写真画像１０１に付加される付加情報の一例である。
写真画像１０１に付加情報が付加されている場合（ＹＥＳ）、処理はＳ１３２に進む。
写真画像１０１に付加情報が付加されていない場合（ＮＯ）、処理はＳ１３６ａに進む。
【００５５】
Ｓ１３２において、撮影パラメータ設定部１３０は、写真画像１０１の付加情報に焦点距離が含まれているか否かを判定する。
写真画像１０１の付加情報に焦点距離が含まれている場合（ＹＥＳ）、処理はＳ１３３ａに進む。
写真画像１０１の付加情報に焦点距離が含まれていない場合（ＮＯ）、処理はＳ１３３ｂに進む。
【００５６】
Ｓ１３３ａにおいて、撮影パラメータ設定部１３０は、写真画像１０１の付加情報から焦点距離を取得する。
Ｓ１３３ａの後、処理はＳ１３４に進む。
【００５７】
Ｓ１３３ｂにおいて、撮影パラメータ設定部１３０は、焦点距離の入力を要求するメッセージをディスプレイに表示する。
そして、利用者は入力機器１０９を用いて焦点距離を入力し、撮影パラメータ設定部１３０は入力された焦点距離を取得する。
Ｓ１３３ｂの後、処理はＳ１３４に進む。
【００５８】
Ｓ１３４において、撮影パラメータ設定部１３０は、写真画像１０１の付加情報に縦方向ＣＣＤサイズが含まれているか否かを判定する。
写真画像１０１の付加情報に縦方向ＣＣＤサイズが含まれている場合（ＹＥＳ）、処理はＳ１３５ａに進む。
写真画像１０１の付加情報に縦方向ＣＣＤサイズが含まれていない場合（ＮＯ）、処理はＳ１３５ｂに進む。
【００５９】
Ｓ１３５ａにおいて、撮影パラメータ設定部１３０は、写真画像１０１の付加情報から縦方向ＣＣＤサイズを取得する。
Ｓ１３５ａの後、処理はＳ１３７に進む。
【００６０】
Ｓ１３５ｂにおいて、撮影パラメータ設定部１３０は、縦方向ＣＣＤサイズの入力を要求するメッセージをディスプレイに表示する。
そして、利用者は入力機器１０９を用いて縦方向ＣＣＤサイズを入力し、撮影パラメータ設定部１３０は入力された縦方向ＣＣＤサイズを取得する。
但し、写真画像１０１の付加情報に焦点距離が含まれる場合、撮影パラメータ設定部１３０は、焦点距離に基づいて縦方向ＣＣＤサイズを算出してもよい。例えば、焦点距離が３５ｍｍ換算の値である場合、縦方向ＣＣＤサイズは３５ｍｍである（ｍｍ：ミリメートル）。
Ｓ１３５ｂの後、処理はＳ１３７に進む。
【００６１】
Ｓ１３６ａにおいて、撮影パラメータ設定部１３０は、Ｓ１３５ｂと同様に、利用者によって入力された縦方向ＣＣＤサイズを取得する。
Ｓ１３６ａの後、処理はＳ１３６ｂに進む。
【００６２】
Ｓ１３６ｂにおいて、撮影パラメータ設定部１３０は、Ｓ１３３ｂと同様に、利用者によって入力された焦点距離を取得する。
Ｓ１３６ｂの後、処理はＳ１３７に進む。
【００６３】
Ｓ１３７において、撮影パラメータ設定部１３０は、縦方向ＣＣＤサイズと焦点距離とを用いて、縦方向視野角を算出する。
例えば、撮影パラメータ設定部１３０は、以下の式（１）を計算することによって縦方向視野角を算出する。
【００６４】
【数９】
2014109440.tif 000011
【００６５】
但し、写真画像１０１の付加情報に縦方向視野角が含まれる場合、撮影パラメータ設定部１３０は写真画像１０１の付加情報から縦方向視野角を取得してもよい。また、撮影パラメータ設定部１３０は、縦方向視野角の入力を要求するメッセージを表示し、利用者によって入力された縦方向視野角を取得してもよい。
Ｓ１３７により、撮影パラメータ設定処理（Ｓ１３０）は終了する。
【００６６】
図２に戻り、画像寸法推定方法の説明を続ける。
Ｓ１３０の後、処理はＳ１４０に進む。
【００６７】
Ｓ１４０において、基準軸ベクトル推定部１４０は、Ｓ１２０で生成された姿勢基準軸データ１２９と、Ｓ１３０で取得された撮影パラメータ１３９とを用いて、写真画像１０１に映っている被写体１０２の姿勢を表す基準軸ベクトル１４９を算出する。基準軸ベクトル１４９とは、被写体１０２の姿勢を表す三次元のベクトル値である。
例えば、基準軸ベクトル推定部１４０は、非特許文献１に開示されている技術を適用して基準軸ベクトル１４９を算出する。
非特許文献１の技術を利用する基準軸ベクトル推定処理（Ｓ１４０）については、別途説明する。
Ｓ１４０の後、処理はＳ１５０に進む。
【００６８】
Ｓ１５０において、リファレンス設定部１５０は、基準軸ベクトル１４９を用いて、写真画像１０１に映っている被写体１０２のリファレンス線分を表すリファレンス線分ベクトル１５９ａを算出する。リファレンス線分とは、被写体１０２のうち長さが既知である線分である。例えば、リファレンス線分は利用者によって指定される。
さらに、リファレンス設定部１５０は、リファレンス線分の長さ（リファレンス線分長１５９ｂ）を取得する。
【００６９】
リファレンス設定部１５０およびリファレンス設定処理（Ｓ１５０）の詳細について、以下に説明する。
【００７０】
図８は、実施の形態１におけるリファレンス設定部１５０の構成図である。
実施の形態１におけるリファレンス設定部１５０の構成について、図８に基づいて説明する。
【００７１】
リファレンス設定部１５０は、基準軸ベクトル投影部１５１、リファレンス線分設定部１５２、リファレンス線分ベクトル算出部１５３（リファレンス線分ベクトル算出部の一例）およびリファレンス線分長取得部１５４を備える。
基準軸ベクトル投影部１５１は、写真画像１０１に映っている被写体１０２に基準軸ベクトル１４９を投影して表示する。
リファレンス線分設定部１５２は、リファレンス線分の情報を取得し、写真画像１０１に映っている被写体１０２にリファレンス線分を描画する。
リファレンス線分ベクトル算出部１５３は、被写体１０２に投影された基準軸ベクトル１４９を用いて、リファレンス線分ベクトル１５９ａを算出する。
リファレンス線分長取得部１５４は、既知であるリファレンス線分長１５９ｂを取得する。
【００７２】
図９は、実施の形態１におけるリファレンス設定処理（Ｓ１５０）を示すフローチャートである。
実施の形態１におけるリファレンス設定処理（Ｓ１５０）について、図９に基づいて説明する。
【００７３】
Ｓ１５１において、基準軸ベクトル投影部１５１は、写真画像１０１に映っている被写体１０２に、基準軸ベクトル１４９を投影して表示する（図１０参照）。
写真画像１０１への投影とは、三次元の基準軸ベクトル１４９を写真画像面の二次元のベクトルに変換することである。写真画像面とは、カメラから焦点距離だけ離れた座標点を含んでカメラの視線方向と直交する平面である。また、写真画像面に被写体１０２を投影した画像が写真画像１０１に相当する。
例えば、基準軸ベクトル投影部１５１は、被写体１０２の３つの線分（姿勢基準軸）の交点、またはいずれかの線分の端点に基準軸ベクトル１４９を投影する。
以下、基準軸ベクトル１４９を表す３つのベクトル［Ｏｘ，Ｏｙ，Ｏｚ］のうち２つのベクトル（例えば、ＯｙとＯｚ）が成す平面を「基準平面」という。
また、被写体１０２のうち大きさが分かっている部分を「リファレンス」という。
図１０は、実施の形態１における基準軸ベクトル投影処理（Ｓ１５１）の概要図である。図１０の「Ｐ０」は基準軸ベクトル１４９の原点を示している。
Ｓ１５１の後、処理はＳ１５２に進む。
【００７４】
Ｓ１５２において、利用者は、被写体１０２のうち大きさが分かっているリファレンスに基準軸ベクトル１４９が投影されているか否かを判定する。
基準軸ベクトル１４９がリファレンスに投影されている場合（ＹＥＳ）、処理はＳ１５６に進む。
基準軸ベクトル１４９がリファレンスに投影されていない場合、つまり、基準軸ベクトル１４９がリファレンス以外の部分に投影されている場合（ＮＯ）、処理はＳ１５３に進む。
【００７５】
Ｓ１５３において、利用者は、マウスを操作し、マウスカーソル１０９ａを用いて基準軸ベクトル１４９を被写体１０２のリファレンスに移動させる。但し、利用者は、マウス以外の入力機器１０９（例えば、キーボード）を用いて、基準軸ベクトル１４９を移動させても構わない。
基準軸ベクトル投影部１５１は、利用者による基準軸ベクトル１４９の移動に従って、基準軸ベクトル１４９を被写体１０２のリファレンスに投影して表示する。
【００７６】
図１１は、実施の形態１における基準軸ベクトル移動処理（Ｓ１５３）の概要図である。
例えば、利用者は、図１１の（１）（２）に示す手順で、基準軸ベクトル１４９を被写体１０２のリファレンスに移動させる。移動前の基準軸ベクトル１４９を点線で記し、移動後の基準軸ベクトル１４９を実線で記す。
（１）利用者は、基準平面（Ｏｙ−Ｏｚ平面）内で基準軸ベクトル１４９を移動させる。
基準軸ベクトル投影部１５１は、基準平面内の移動先に基準軸ベクトル１４９を投影する。
（２）利用者は、基準平面と直交する法線ベクトル（Ｏｘ）の方向に基準軸ベクトル１４９を移動させる。
基準軸ベクトル投影部１５１は、移動先のリファレンスに基準軸ベクトル１４９を投影する。
以下、リファレンスに移動後の基準軸ベクトル１４９が表す基準平面（Ｏｙ−Ｏｚ平面）を「リファレンス平面」という。
Ｓ１５３（図９）の後、処理はＳ１５４に進む。
【００７７】
Ｓ１５４において、利用者は、マウスを操作し、マウスカーソル１０９ａを用いて被写体１０２のリファレンスにリファレンス線分を指定する。リファレンス線分とは、リファレンスのうち長さが既知である線分である。利用者は、マウス以外の入力機器１０９を利用してリファレンス線分を指定しても構わない。
リファレンス線分設定部１５２は、利用者の指定に従ってリファレンス線分を描画する（図１２参照）。
図１２は、実施の形態１におけるリファレンス線分設定処理（Ｓ１５４）の概要図である。
Ｓ１５４の後、処理はＳ１５６に進む。
【００７８】
Ｓ１５５において、利用者は、Ｓ１５４と同様に、被写体１０２のリファレンスにリファレンス線分を指定する。
リファレンス線分設定部１５２は、利用者の指定に従ってリファレンス線分を描画する（図１３参照）。
図１３は、実施の形態１におけるリファレンス線分設定処理（Ｓ１５５）の概要図である。
Ｓ１５５の後、処理はＳ１５６に進む。
【００７９】
Ｓ１５６において、リファレンス線分ベクトル算出部１５３は、基準軸ベクトル１４９を用いてリファレンス線分ベクトル１５９ａを算出する。
リファレンス線分ベクトル１５９ａは、リファレンス線分を表す三次元のベクトル値である。
【００８０】
図１４は、実施の形態１におけるリファレンス線分ベクトル算出処理（Ｓ１５６）の概要図である。
例えば、リファレンス線分ベクトル算出部１５３は、図１４の（１）（２）に示す手順で、リファレンス線分ベクトルＬｒを算出する。
（１）リファレンス線分ベクトル算出部１５３は、写真画像面ｌｖ（Ｅｚ）内のリファレンス線分の端点Ｐｓおよび端点Ｐｅの座標値を取得する。
そして、リファレンス線分ベクトル算出部１５３は、カメラの中心点Ｐｃから端点Ｐｓまたは端点Ｐｅへの向きを示す単位ベクトルＶｓおよび単位ベクトルＶｅを算出する。カメラの中心点Ｐｃは、写真画像１０１を撮像したときのカメラの位置を示す座標点、つまり、写真画像１０１が撮像された位置を示す座標点である。
（２）リファレンス線分ベクトル算出部１５３は、単位ベクトルＶｓおよび単位ベクトルＶｅを用いてリファレンス線分ベクトルＬｒを算出する。
リファレンス線分ベクトルＬｒは、以下の各式を用いて算出することができる。
【００８１】
リファレンスが成すリファレンス平面Γｒは、以下の式（２）で表すことができる。
【００８２】
【数１０】
2014109440.tif 000012
【００８３】
カメラの座標点Ｐｃとリファレンス線分の端点Ｐｓとを通る直線ｌｓは、単位ベクトルＶｓを用いて、以下の式（３ｓ）で表すことができる。
同様に、カメラの座標点Ｐｃとリファレンス線分の端点Ｐｅとを通る直線ｌｅは、単位ベクトルＶｅを用いて、以下の式（３ｅ）で表すことができる。
【００８４】
【数１１】
2014109440.tif 000013
【００８５】
上記の直線ｌｓの式（３ｓ）をｘｙｚ成分に分解すると以下の式（４ｓ）が得られる。
同様に、上記の直線ｌｅの式（３ｅ）をｘｙｚ成分に分解すると以下の式（４ｅ）が得られる。
【００８６】
【数１２】
2014109440.tif 000014
【００８７】
上記の直線ｌｓの分解式（４ｓ）を上記のリファレンス平面Γｒの式（２）に代入して展開すると、以下の式（５ｓ）が得られる。
同様に、直線ｌｅの分解式（４ｅ）を上記のリファレンス平面Γｒの式（２）に代入して展開すると、以下の式（５ｅ）が得られる。
【００８８】
【数１３】
2014109440.tif 000015
【００８９】
但し、上記の式（５ｓ）（５ｅ）において、リファレンス平面ΓｒがＯｘ−Ｏｙ平面でなく、Ｏｘ−Ｏｚ平面である場合、基準軸ベクトルＯｚ（ｘｏｚ，ｙｏｚ，ｚｏｚ）は基準軸ベクトルＯｙ（ｘｏｙ，ｙｏｙ，ｚｏｙ）に置き換わる。
同様に、リファレンス平面がＯｙ−Ｏｚ平面である場合、基準軸ベクトルＯｚ（ｘｏｚ，ｙｏｚ，ｚｏｚ）は基準軸ベクトルＯｘ（ｘｏｘ，ｙｏｘ，ｚｏｘ）に置き換わる。
【００９０】
リファレンス線分ベクトル算出部１５３は、上記の式（５ｓ）を計算して媒介変数ｔｓを算出し、上記の式（５ｅ）を計算して媒介変数ｔｅを算出する。
【００９１】
媒介変数ｔｓは端点Ｐｓを得るための変数値である。つまり、媒介変数ｔｓはカメラの中心点Ｐｃから端点ＰｓへのベクトルＶｓｔｓを得るための変数値である。
同様に、媒介変数ｔｅは端点Ｐｅを得るための変数値である。つまり、媒介変数ｔｅはカメラの中心点Ｐｃから端点ＰｅへのベクトルＶｅｔｅを得るための変数値である。
【００９２】
リファレンス線分ベクトル算出部１５３は、媒介変数値ｔｓと媒介変数値ｔｅとを以下の式（６）に代入し、リファレンス線分ベクトルＬｒを算出する。
【００９３】
【数１４】
2014109440.tif 000016
【００９４】
上記のＳ１５６により、リファレンス線分ベクトル１５９ａ（Ｌｒ）が算出される。
Ｓ１５６（図９）の後、処理はＳ１５７に進む。
【００９５】
Ｓ１５７において、利用者は、入力機器１０９を用いて、既知であるリファレンス線分長１５９ｂを入力する。
リファレンス線分長取得部１５４は、入力されたリファレンス線分長１５９ｂを取得する。
Ｓ１５７の後、リファレンス設定処理（Ｓ１５０）は終了する。
【００９６】
図２に戻り、画像寸法推定方法の説明を続ける。
リファレンス設定処理（Ｓ１５０）の後、処理はＳ１６０に進む。
【００９７】
Ｓ１６０において、寸法推定部１６０は、リファレンス設定処理（Ｓ１５０）と同様に、基準軸ベクトル１４９を用いて計測線分ベクトルを算出する。
そして、寸法推定部１６０は、計測線分ベクトルと、リファレンス線分ベクトル１５９ａと、リファレンス線分長１５９ｂとを用いて、計測線分長１６９を算出する。
計測線分ベクトルとは、計測線分を表す三次元のベクトル値である。計測線分とは、被写体１０２のうち長さが未知である線分、つまり、計測する対象部分の線分である。例えば、計測線分は利用者によって指定される。
【００９８】
寸法推定部１６０および寸法推定処理（Ｓ１６０）の詳細について、以下に説明する。
【００９９】
図１５は、実施の形態１における寸法推定部１６０の構成図である。
実施の形態１における寸法推定部１６０の構成について、図１５に基づいて説明する。
【０１００】
寸法推定部１６０は、基準軸ベクトル投影部１６１、計測線分設定部１６２、計測線分ベクトル算出部１６３（計測線分ベクトル算出部の一例）および計測線分長算出部１６４（計測線分長算出部の一例）を備える。
基準軸ベクトル投影部１６１は、写真画像１０１に映っている被写体１０２に基準軸ベクトル１４９を投影して表示する。
計測線分設定部１６２は、計測線分の情報を取得し、写真画像１０１に映っている被写体１０２に計測線分を描画する。
計測線分ベクトル算出部１６３は、被写体１０２に投影された基準軸ベクトル１４９を用いて、計測線分ベクトルを算出する。
計測線分長算出部１６４は、計測線分ベクトルと、リファレンス線分ベクトル１５９ａと、リファレンス線分長１５９ｂとを用いて、計測線分長１６９を算出する。
【０１０１】
図１６は、実施の形態１における寸法推定処理（Ｓ１６０）を示すフローチャートである。
実施の形態１における寸法推定処理（Ｓ１６０）について、図１６に基づいて説明する。
【０１０２】
寸法推定処理（Ｓ１６０）のＳ１６１からＳ１６６は、リファレンス設定処理（Ｓ１５０）（図９参照）のＳ１５１からＳ１５６と同様である。
つまり、Ｓ１５１からＳ１５６の「リファレンス」を「計測部分」に置き換え、「リファレンス平面」を「計測平面」に置き換え、「リファレンス線分」を「計測線分」に置き換えた処理がＳ１６１からＳ１６６である。
【０１０３】
図１７は、実施の形態１におけるＳ１６１からＳ１６５の概要図である。
Ｓ１６１からＳ１６５の処理概要を図１７に示す。図１７はＳ１５１からＳ１５５の概要図（図１１から図１３）に相当する。
基準軸ベクトル投影部１６１は、写真画像１０１に映っている被写体１０２に基準軸ベクトル１４９を投影して表示する（Ｓ１６１）。
基準軸ベクトル１４９が被写体１０２の計測部分に投影されていない場合（Ｓ１６２：ＹＥＳ）、基準軸ベクトル投影部１６１は、基準軸ベクトル１４９を被写体１０２の計測部分に移動する（Ｓ１６３）。
計測線分設定部１６２は、利用者の指定に従って計測線分を描画する（Ｓ１６４、Ｓ１６５）。
【０１０４】
図１８は、実施の形態１におけるＳ１６６の概要図である。
Ｓ１６６の処理概要を図１８に示す。図１８はＳ１５６の概要図（図１４）に相当する。
計測線分ベクトル算出部１６３は、計測部分に投影された基準軸ベクトル１４９を用いて、計測線分ベクトルＬｍを算出する（Ｓ１６６）。
【０１０５】
つまり、計測線分ベクトル算出部１６３は、カメラの中心点Ｐｃから計測線分の端点Ｐｍｓへの単位ベクトルＶｍｓと、カメラの中心点Ｐｃから計測線分の端点Ｐｍｅへの単位ベクトルＶｍｅとを算出する。
また、計測線分ベクトル算出部１６３は、以下の式（７ｓ）を計算して媒介変数ｔｍｓを算出し、以下の式（７ｅ）を計算して媒介変数ｔｍｅを算出する。
そして、計測線分ベクトル算出部１６３は、媒介変数ｔｍｓと媒介変数ｔｍｅとを用いて以下の式（７ｍ）を計算し、計測線分ベクトルＬｍを算出する。
【０１０６】
【数１５】
2014109440.tif 000017
【０１０７】
但し、上記の式（７ｓ）（７ｅ）において、計測平面ΓｍがＯｙ−Ｏｚ平面でなく、Ｏｘ−Ｏｚ平面である場合、基準軸ベクトルＯｘ（ｘｏｘ，ｙｏｘ，ｚｏｘ）は基準軸ベクトルＯｙ（ｘｏｙ，ｙｏｙ，ｚｏｙ）に置き換わる。
同様に、リファレンス平面がＯｘ−Ｏｙ平面である場合、基準軸ベクトルＯｘ（ｘｏｘ，ｙｏｘ，ｚｏｘ）は基準軸ベクトルＯｚ（ｘｏｚ，ｙｏｚ，ｚｏｚ）に置き換わる。
【０１０８】
媒介変数ｔｍｓは端点Ｐｍｓを得るための変数値である。つまり、媒介変数ｔｍｓはカメラの中心点Ｐｃから端点ＰｍｓへのベクトルＶｍｓｔｍｓを得るための変数値である。
同様に、媒介変数ｔｍｅは端点Ｐｍｅを得るための変数値である。つまり、媒介変数ｔｍｅはカメラの中心点Ｐｃから端点ＰｍｅへのベクトルＶｍｅｔｍｅを得るための変数値である。
【０１０９】
Ｓ１６７（図１６参照）において、計測線分長算出部１６４は、Ｓ１６６で算出した計測線分ベクトルと、リファレンス設定部１５０によって算出されたリファレンス線分ベクトル１５９ａと、リファレンス設定部１５０によって取得されたリファレンス線分長１５９ｂとを用いて、計測線分長１６９を算出する。
【０１１０】
以下に、計測線分長ｌｍ＿ｔｒｕｅを算出するための式（８）を示す。
以下の式（８）において、「Ｌｍ」は計測線分ベクトルであり、「Ｌｒ」はリファレンス線分ベクトル１５９ａであり、「ｌｒ＿ｔｒｕｅ」はリファレンス線分長１５９ｂである。
【０１１１】
【数１６】
2014109440.tif 000018
【０１１２】
つまり、計測線分長算出部１６４は、リファレンス線分ベクトルに対する計測線分ベクトルの大きさの比率をリファレンス線分長に乗じた値を計測線分長１６９として算出する。
Ｓ１６７の後、寸法推定処理（Ｓ１６０）は終了する。
【０１１３】
図２に戻り、画像寸法推定方法の説明を続ける。
寸法推定処理（Ｓ１６０）の後、処理はＳ１７０に進む。
【０１１４】
Ｓ１７０において、推定結果出力部１７０は、Ｓ１６０で算出された計測線分長１６９を出力する。
例えば、推定結果出力部１７０は、計測線分長１６９をディスプレイ９１１に表示する。
また、推定結果出力部１７０は、計測線分長１６９と写真画像１０１（計測線分または計測線分ベクトルなどの情報を含む）とを対応付けて写真画像データベース１８０に記憶する。
計測部分は新たなリファレンスとして利用することができ、計測線分長１６９は新たな既知のリファレンス線分長として利用することができる。
【０１１５】
次に、非特許文献１の技術を利用する基準軸ベクトル推定処理（Ｓ１４０）について説明する。
【０１１６】
図１９は、実施の形態１における基準軸ベクトル推定処理（Ｓ１４０）を示すフローチャートである。
実施の形態１における基準軸ベクトル推定処理（Ｓ１４０）について、図１９に基づいて説明する。
【０１１７】
Ｓ１４１において、基準軸ベクトル推定部１４０は、姿勢基準軸データ１２９に基づいて、被写体１０２の姿勢基準軸［ｌｘ，ｌｙ，ｌｚ］が図４のように一点で交わっているか、または図５のように一点で交わっていないか判定する。
姿勢基準軸［ｌｘ，ｌｙ，ｌｚ］が一点で交わる場合（ＹＥＳ）、処理はＳ１４２に進む。
姿勢基準軸［ｌｘ，ｌｙ，ｌｚ］が一点で交わらない場合（ＮＯ）、処理はＳ１４３に進む。
【０１１８】
Ｓ１４２において、基準軸ベクトル推定部１４０は、非特許文献１で定義されたＦＰＭ座標系を用いて基準軸ベクトル１４９を算出する。
ＦＰＭはＦｉｒｓｔ  Ｐｅｒｓｐｅｃｔｉｖｅ  Ｍｏｖｉｎｇの略である。
基準軸ベクトル算出処理［ＦＰＭ座標系］（Ｓ１４２）の詳細については別途説明する。
Ｓ１４２の後、基準軸ベクトル推定処理（Ｓ１４０）は終了する。
【０１１９】
Ｓ１４３において、基準軸ベクトル推定部１４０は、非特許文献１で定義されたＳＰＭ座標系を用いて基準軸ベクトル１４９を算出する。
ＳＰＭはＳｅｃｏｎｄ  Ｐｅｒｓｐｅｃｔｉｖｅ  Ｍｏｖｉｎｇの略である。
基準軸ベクトル推定処理［ＳＰＭ座標系］（Ｓ１４３）の詳細については別途説明する。
Ｓ１４３の後、基準軸ベクトル推定処理（Ｓ１４０）は終了する。
【０１２０】
図２０は、実施の形態１における基準軸ベクトル算出処理［ＦＰＭ座標系］（Ｓ１４２）を示すフローチャートである。
ＦＰＭ座標系を用いる基準軸ベクトル算出処理（Ｓ１４２）について、図２０に基づいて説明する。
【０１２１】
Ｓ１４２Ａにおいて、基準軸ベクトル推定部１４０は、カメラの中心点Ｐｃから写真画像面ｌｖ（Ｅｚ）内の姿勢基準軸［ｌｘ，ｌｙ，ｌｚ］の端点への単位ベクトル［Ｅｘ，Ｅｙ，Ｅｚ］を算出する（図２１参照）。
以下、単位ベクトル［Ｅｘ，Ｅｙ，Ｅｚ］を「端点単位ベクトル」という。
図２１は、実施の形態１における端点単位ベクトル算出処理（Ｓ１４２Ａ）の概要図である。
Ｓ１４２Ａの後、処理はＳ１４２Ｂに進む。
【０１２２】
Ｓ１４２Ｂにおいて、基準軸ベクトル推定部１４０は、カメラの中心点Ｐｃから写真画像面ｌｖ（Ｅｚ）内の姿勢基準軸の交点ｌ０への単位ベクトルＥｃを算出する（図２２参照）。
以下、単位ベクトルＥｃを「交点単位ベクトル」という。
図２２は、実施の形態１における交点単位ベクトル算出処理（Ｓ１４２Ｂ）の概要図である。
Ｓ１４２Ｂの後、処理はＳ１４２Ｃに進む。
【０１２３】
Ｓ１４２Ｃにおいて、基準軸ベクトル推定部１４０は、カメラの中心点Ｐｃから交点単位ベクトルＥｃの方向に焦点距離ｆだけ離れた座標点Ｐｆを算出する（図２３の（１）参照）。
基準軸ベクトル推定部１４０は、座標点Ｐｆを含んで交点単位ベクトルＥｃと直交する平面ｌｖ（Ｅｃ）を算出する（図２３の（２）参照）。つまり、交点単位ベクトルＥｃは平面ｌｖ（Ｅｃ）の法線ベクトルになる。
以下、平面ｌｖ（Ｅｃ）を「ＦＰＭ座標面」という。また、座標点Ｐｆを「ＦＰＭ座標原点」という。
図２３は、実施の形態１におけるＦＰＭ座標面算出処理（Ｓ１４２Ｃ）の概要図である。
Ｓ１４２Ｃの後、処理はＳ１４２Ｄに進む。
【０１２４】
Ｓ１４２Ｄにおいて、基準軸ベクトル推定部１４０は、端点単位ベクトル［Ｅｘ，Ｅｙ，Ｅｚ］を含んだ直線とＦＰＭ座標面ｌｖ（Ｅｃ）との交点［Ｐｘ’，Ｐｙ’，Ｐｚ’］を算出する（図２４の（１）参照）。
基準軸ベクトル推定部１４０は、ＦＰＭ座標原点Ｐｆから交点［Ｐｘ’，Ｐｙ’，Ｐｚ’］へのベクトル［Ｇｘ，Ｇｙ，Ｇｚ］を算出する（図２４の（２）参照）。
以下、ベクトル［Ｇｘ，Ｇｙ，Ｇｚ］を「ＦＰＭ座標面ベクトル」という。
図２４は、実施の形態１におけるＦＰＭ座標面ベクトル算出処理（Ｓ１４２Ｄ）の概要図である。
Ｓ１４２Ｄの後、処理はＳ１４２Ｅに進む。
【０１２５】
Ｓ１４２Ｅにおいて、基準軸ベクトル推定部１４０は、ＦＰＭ座標面ベクトル［Ｇｘ，Ｇｙ，Ｇｚ］のいずれかのベクトル（例えば、Ｇｘ）を選択する。
基準軸ベクトル推定部１４０は、選択したＦＰＭ座標面ベクトル（Ｇｘ）の単位ベクトルＥａ（＝Ｇｘ／｜Ｇｘ｜）を算出する（図２５参照）。
基準軸ベクトル推定部１４０は、交点単位ベクトルＥｃと単位ベクトルＥａとの外積ベクトルＥｂ（＝Ｅｃ×Ｅａ）を算出する（図２５参照）。
単位べクトルＥａはＦＰＭ座標系のＡ軸の単位ベクトルであり、外積ベクトルＥｂはＦＰＭ座標系のＢ軸の単位ベクトルであり、交点単位ベクトルＥｃはＦＰＭ座標系のＣ軸の単位ベクトルである。単位ベクトル［Ｅａ，Ｅｂ，Ｅｃ］は「正規直交基底」ともいう。
図２５は、実施の形態１におけるＦＰＭ座標系の単位ベクトル算出処理（Ｓ１４２Ｅ）の概要図である。
Ｓ１４２Ｅの後、処理はＳ１４２Ｆに進む。
【０１２６】
Ｓ１４２Ｆにおいて、基準軸ベクトル推定部１４０は、ＦＰＭ座標面ベクトル［Ｇｘ，Ｇｙ，Ｇｚ］を用いて、被写体座標系の基準軸ベクトル［Ｅ１，Ｅ２，Ｅ３］を算出する。
被写体座標系とは被写体１０２に設定したローカル座標系である。また、基準軸ベクトル［Ｅ１，Ｅ２，Ｅ３］は、被写体座標系の単位ベクトル（正規直交基底）である。
図２６は、実施の形態１における被写体座標系の基準軸ベクトルの関係図である。
図２６において、ＧｚとＧｙとが成す見かけの角度Δβ、ＧｘとＧｙとが成す見かけの角度β２およびＧｘとＧｚとが成す見かけの角度β１は、以下の式（１１）で表すことができる。
【０１２７】
【数１７】
2014109440.tif 000019
【０１２８】
さらに、Ｅ１とＥ２とが成す実際の角度ΔαのｔａｎΔαは、Ｅ３とＦＰＭ座標系のＣ軸とが成す角度ψとを用いて、以下の式（１２）で表すことができる。
【０１２９】
【数１８】
2014109440.tif 000020
【０１３０】
ここで、被写体１０２の基準軸ベクトル［Ｅ１，Ｅ２］が成す角度Δαは９０度であるから、上記の式（１２）は、以下の式（１３）に展開することができる。
【０１３１】
【数１９】
2014109440.tif 000021
【０１３２】
そして、被写体座標系の基準軸ベクトル［Ｅ１，Ｅ２，Ｅ３］は、上記の式（１１）と式（１３）との関係により、以下の式（１４）で表すことができる。
以下の式（１４）において、α１はＥ３とＥ１とが成す実際の角度であり、α２はＥ３とＥ２とが成す実際の角度である。
【０１３３】
【数２０】
2014109440.tif 000022
【０１３４】
Ｓ１４２Ｆ（図２０）の後、処理はＳ１４２Ｇに進む。
【０１３５】
Ｓ１４２Ｇにおいて、基準軸ベクトル推定部１４０は、被写体座標系の基準軸ベクトル［Ｅ１，Ｅ２，Ｅ３］をワールド座標系の基準軸ベクトル［Ｏｘ，Ｏｙ，Ｏｚ］に変換するための変換行列Ｍｈを算出する。変換行列Ｍｈは以下の式（１５）を用いて表すことができる。
【０１３６】
【数２１】
2014109440.tif 000023
【０１３７】
上記の式（１５）で用いている符号の意味は以下の通りである。
（ｘＥａ，ｙＥａ，ｚＥａ）は、ＦＰＭ座標系のＡ軸の単位ベクトルＥａ（図２５参照）をワールド座標系で表したベクトル値である。
（ｘＥｂ，ｙＥｂ，ｚＥｂ）は、ＦＰＭ座標系のＢ軸の単位ベクトルＥｂ（図２５参照）をワールド座標系で表したベクトル値である。
（ｘＥｃ，ｙＥｃ，ｚＥｃ）は、ＦＰＭ座標系のＣ軸の単位ベクトルＥｃ（図２５参照）をワールド座標系で表したベクトル値である。
（ａＥ１，ｂＥ１，ｃＥ１）は、被写体座標系の基準軸ベクトルＥ１（図２６参照）をＦＰＭ座標系で表したベクトル値である。
（ａＥ２，ｂＥ２，ｃＥ２）は、被写体座標系の基準軸ベクトルＥ２（図２６参照）をＦＰＭ座標系で表したベクトル値である。
（ａＥ３，ｂＥ３，ｃＥ３）は、被写体座標系の基準軸ベクトルＥ３（図２６参照）をＦＰＭ座標系で表したベクトル値である。
【０１３８】
そして、基準軸ベクトル推定部１４０は、上記の変換行列Ｈｈを用いて、被写体座標系の基準軸ベクトル［Ｅ１，Ｅ２，Ｅ３］をワールド座標系の基準軸ベクトル［Ｏｘ，Ｏｙ，Ｏｚ］に変換し、ワールド座標系の基準軸ベクトル［Ｏｘ，Ｏｙ，Ｏｚ］を被写体１０２の基準軸ベクトル１４９として出力する。
以下に、被写体座標系の基準軸ベクトル［Ｅ１，Ｅ２，Ｅ３］をワールド座標系の基準軸ベクトル［Ｏｘ，Ｏｙ，Ｏｚ］に変換するための式（１６）を示す。
【０１３９】
【数２２】
2014109440.tif 000024
【０１４０】
Ｓ１４２Ｇにより、基準軸ベクトル算出処理［ＦＰＭ座標系］（Ｓ１４２）は終了する。
【０１４１】
図２７は、実施の形態１における基準軸ベクトル算出処理［ＳＰＭ座標系］（Ｓ１４３）を示すフローチャートである。
ＳＰＭ座標系を用いる基準軸ベクトル算出処理（Ｓ１４３）について、図２７に基づいて説明する。
【０１４２】
Ｓ１４３Ａにおいて、基準軸ベクトル推定部１４０は、写真画像面ｌｖ（Ｅｚ）内の姿勢基準軸［ｌｘ，ｌｙ，ｌｚ］から２つの姿勢基準軸（例えば、ｌｘとｌｙ）を選択する。
基準軸ベクトル推定部１４０は、選択した２つの姿勢基準軸（ｌｘ，ｌｙ）の延長上の交点Ｐｘｙを算出する（図２８の（１）参照）。
基準軸ベクトル推定部１４０は、カメラの中心点Ｐｃから交点Ｐｘｙへの単位ベクトルＥｘｙを算出する（図２８の（２）参照）。以下、単位ベクトルＥｘｙを「交点単位ベクトル」という。
図２８は、実施の形態１における交点単位ベクトル算出処理（Ｓ１４３Ａ）の概要図である。
Ｓ１４３Ａの後、処理はＳ１４３Ｂに進む。
【０１４３】
Ｓ１４３Ｂにおいて、基準軸ベクトル推定部１４０は、カメラの中心点Ｐｃから交点単位ベクトルＥｘｙの方向に焦点距離ｆだけ離れた座標点Ｐｆを算出する（図２９の（１）参照）。
基準軸ベクトル推定部１４０は、座標点Ｐｆを含んで交点単位ベクトルＥｘｙと直交する平面ｌｖ（Ｅｃ）を算出する（図２９の（２）参照）。つまり、交点単位ベクトルＥｘｙは平面ｌｖ（Ｅｃ）の法線ベクトルになる。
以下、平面ｌｖ（Ｅｃ）を「ＳＰＭ座標面」という。また、座標点Ｐｆを「ＳＰＭ座標原点」という。
図２９は、実施の形態１におけるＳＰＭ座標面算出処理（Ｓ１４３Ｂ）の概要図である。
Ｓ１４３Ｂの後、処理はＳ１４３Ｃに進む。
【０１４４】
Ｓ１４３Ｃにおいて、基準軸ベクトル推定部１４０は、姿勢基準軸［ｌｘ，ｌｙ，ｌｚ］を写真画像面ｌｖ（Ｅｚ）からＳＰＭ座標面ｌｖ（Ｅｃ）に投影する（図３０参照）。
以下、ＳＰＭ座標面ｌｖ（Ｅｃ）に投影した姿勢基準軸を［ｌｘ’，ｌｙ’，ｌｚ’］と記す。
図３０は、実施の形態１における姿勢基準軸投影処理（Ｓ１４３Ｃ）の概要図である。
Ｓ１４３Ｃの後、処理はＳ１４３Ｄに進む。
【０１４５】
Ｓ１４３Ｄにおいて、基準軸ベクトル推定部１４０は、ＳＰＭ座標面ｌｖ（Ｅｃ）に投影した姿勢基準軸［ｌｘ’，ｌｙ’，ｌｚ’］のうち、Ｓ１４３Ａで選択した姿勢基準軸（ｌｘ，ｌｙ）を投影後の姿勢基準軸（ｌｘ’，ｌｙ’）を選択する。
基準軸ベクトル推定部１４０は、選択した姿勢基準軸（ｌｘ’，ｌｙ’）が成す角の二等分線をＳＰＭ座標系のＳ軸に決定する（図３１の（１）参照）。
さらに、基準軸ベクトル推定部１４０は、ＳＰＭ座標面ｌｖ（Ｅｃ）の法線ベクトルＥｘｙをＳＰＭ座標面のＣ軸に決定し、Ｓ軸とＣ軸とに直交する軸をＳＰＭ座標系のＴ軸に決定する（図３１の（１）参照）。
そして、基準軸ベクトル推定部１４０は、ＳＰＭ座標系のＳ軸の単位ベクトルＥｓと、Ｔ軸の単位ベクトルＥｔと、Ｃ軸の単位ベクトルＥｃとを算出する（図３１の（２）参照）。
図３１は、実施の形態１におけるＳＰＭ座標系の単位ベクトル算出処理（Ｓ１４３Ｄ）の概要図である。
【０１４６】
以下に、ＳＰＭ座標系の単位ベクトル［Ｅｓ，Ｅｔ，Ｅｃ］を算出する式（２１）を示す。
以下の式（２１）で用いている「Ｇｘ’」は姿勢基準軸ｌｘ’と同じ向きを有する単位ベクトルである。また、「Ｇｙ’」は姿勢基準軸ｌｙ’と同じ向きを有する単位ベクトルである。
【０１４７】
【数２３】
2014109440.tif 000025
【０１４８】
Ｓ１４３Ｄの後、処理はＳ１４３Ｅに進む。
【０１４９】
Ｓ１４３Ｅにおいて、基準軸ベクトル推定部１４０は、以下のように、ＳＰＭ座標面内の姿勢基準軸［ｌｘ’，ｌｙ’，ｌｚ’］を用いて被写体座標系の基準軸ベクトル［Ｅ１，Ｅ２，Ｅ３］を算出する。
被写体座標系とは被写体１０２に設定したローカル座標系である。また、基準軸ベクトル［Ｅ１，Ｅ２，Ｅ３］は、被写体座標系の単位ベクトル（正規直交基底）である。
【０１５０】
図３２は、実施の形態１におけるＳＰＭ座標面内の姿勢基準軸の関係図である。
図３２に示すＰＡＴ曲線は非特許文献１で定義されている曲線である。ＰＡＴはＰｅｒｓｐｅｃｔｉｖｅ  Ａｎｇｌｅ  Ｔｒａｎｓｆｏｒｍの略である。
ＰＡＴ曲線は、ｌｘ’とｌｙ’とが成す角度Δβと焦点距離ｆとを用いて、以下の式（２２）で表すことができる。
【０１５１】
【数２４】
2014109440.tif 000026
【０１５２】
また、ｌｚ’を含んだ直線の式は、変数γおよび変数δを用いて以下の式（２３）で表すことができる。
【０１５３】
【数２５】
2014109440.tif 000027
【０１５４】
ＰＡＴ曲線と直線ｌｚ’との交点Ｐ１の極座標値は、上記の式（２２）と上記の式（２３）とを用いて、以下の式（２４）で算出することができる。
【０１５５】
【数２６】
2014109440.tif 000028
【０１５６】
また、ＰＡＴ曲線と直線ｌｚ’との交点Ｐ１の極座標値は、ＰｆからＰ１までの直線の長さａと、ＰｆからＰ１までの直線とＳ軸とが成す角度ηとを用いて、（ａｃｏｓη，ａｓｉｎη）で表すことができる。
【０１５７】
ここで、被写体座標系の基準軸ベクトル［Ｅ１，Ｅ２，Ｅ３］は、以下の式（２５）で表すことができる。
【０１５８】
【数２７】
2014109440.tif 000029
【０１５９】
基準軸ベクトル推定部１４０は、上記の式（２５）を計算することにより、被写体座標系の基準軸ベクトル［Ｅ１，Ｅ２，Ｅ３］を算出する。
Ｓ１４３Ｅ（図２７）の後、処理はＳ１４３Ｆに進む。
【０１６０】
Ｓ１４３Ｆにおいて、基準軸ベクトル推定部１４０は、写真画像面ｌｖ（Ｅｚ）から基準点ｌ０を選択する。例えば、基準軸ベクトル推定部１４０は、姿勢基準軸ｌｘの端点を基準点ｌ０として選択する（図３３の（１）参照）。
基準軸ベクトル推定部１４０は、カメラの中心点Ｐｃから基準点ｌ０への単位ベクトルＥ０を算出する（図３３の（２）参照）。
基準軸ベクトル推定部１４０は、単位ベクトルＥ０を用いて、カメラの中心点Ｐｃから基準点ｌ０への方向に焦点距離ｆだけ離れた座標点Ｐ０を算出する（図３３の（２）参照）。
以下、座標点Ｐ０を「基準原点」という。
図３３は、実施の形態１における基準原点算出処理（Ｓ１４３Ｆ）の概要図である。
Ｓ１４３Ｆの後、処理はＳ１４３Ｇに進む。
【０１６１】
Ｓ１４３Ｇにおいて、基準軸ベクトル推定部１４０は、被写体座標系の基準軸ベクトル［Ｅ１，Ｅ２，Ｅ３］をワールド座標系の基準軸ベクトル［Ｏｘ，Ｏｙ，Ｏｚ］に変換するための変換行列Ｍｈを算出する。変換行列Ｍｈは以下の式（２６）で表すことができる。
【０１６２】
【数２８】
2014109440.tif 000030
【０１６３】
上記の式（２６）で用いている符号の意味は以下の通りである。
（ｘｐ０，ｙｐ０，ｚｐ０）は、基準原点Ｐ０（図３３の（２）参照）をワールド座標系（ｘ，ｙ，ｚ）で表した座標値である。
（ｘＥｓ，ｙＥｓ，ｚＥｓ）は、ＳＰＭ座標系のＳ軸の単位ベクトルＥｓ（図３１の（２）参照）をワールド座標系で表したベクトル値である。
（ｘＥｔ，ｙＥｔ，ｚＥｔ）は、ＦＰＭ座標系のＴ軸の単位ベクトルＥｔ（図３１の（２）参照）をワールド座標系で表したベクトル値である。
（ｘＥｃ，ｙＥｃ，ｚＥｃ）は、ＦＰＭ座標系のＣ軸の単位ベクトルＥｃ（＝Ｅｘｙ）をワールド座標系で表したベクトル値である。
（ａＥ１，ｂＥ１，ｃＥ１）は、被写体座標系の基準軸ベクトルＥ１（式（２５）参照）をＦＰＭ座標系で表したベクトル値である。
（ａＥ２，ｂＥ２，ｃＥ２）は、被写体座標系の基準軸ベクトルＥ２（式（２５）参照）をＦＰＭ座標系で表したベクトル値である。
（ａＥ３，ｂＥ３，ｃＥ３）は、被写体座標系の基準軸ベクトルＥ３（式（２５）参照）をＦＰＭ座標系で表したベクトル値である。
【０１６４】
そして、基準軸ベクトル推定部１４０は、上記の変換行列Ｈｈを用いて、被写体座標系の基準軸ベクトル［Ｅ１，Ｅ２，Ｅ３］をワールド座標系の基準軸ベクトル［Ｏｘ，Ｏｙ，Ｏｚ］に変換し、ワールド座標系の基準軸ベクトル［Ｏｘ，Ｏｙ，Ｏｚ］を被写体１０２の基準軸ベクトル１４９として出力する。
以下に、被写体座標系の基準軸ベクトル［Ｅ１，Ｅ２，Ｅ３］をワールド座標系の基準軸ベクトル［Ｏｘ，Ｏｙ，Ｏｚ］に変換するための式（２７）を示す。
【０１６５】
【数２９】
2014109440.tif 000031
【０１６６】
Ｓ１４３Ｇにより、基準軸ベクトル算出処理［ＳＰＭ座標系］（Ｓ１４３）は終了する。
【０１６７】
図３４は、実施の形態１における画像寸法推定装置１００のハードウェア資源の一例を示す図である。
図３４において、画像寸法推定装置１００（コンピュータの一例）は、ＣＰＵ９０１（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）を備えている。ＣＰＵ９０１は、バス９０２を介してＲＯＭ９０３、ＲＡＭ９０４、通信ボード９０５（通信装置）、ディスプレイ９１１（表示装置）、キーボード９１２、マウス９１３、ドライブ９１４、磁気ディスク装置９２０などのハードウェアデバイスと接続され、これらのハードウェアデバイスを制御する。ドライブ９１４は、ＦＤ（Ｆｌｅｘｉｂｌｅ  Ｄｉｓｋ）、ＣＤ（Ｃｏｍｐａｃｔ  Ｄｉｓｃ）、ＤＶＤ（Ｄｉｇｉｔａｌ  Ｖｅｒｓａｔｉｌｅ  Ｄｉｓｃ）などの記憶媒体を読み書きする装置である。
ＲＯＭ９０３、ＲＡＭ９０４、磁気ディスク装置９２０およびドライブ９１４は記憶装置の一例である。キーボード９１２、マウス９１３および通信ボード９０５は入力装置の一例である。ディスプレイ９１１および通信ボード９０５は出力装置の一例である。
【０１６８】
通信ボード９０５は、有線または無線で、ＬＡＮ（Ｌｏｃａｌ  Ａｒｅａ  Ｎｅｔｗｏｒｋ）、インターネット、電話回線などの通信網に接続している。
【０１６９】
磁気ディスク装置９２０には、ＯＳ９２１（オペレーティングシステム）、プログラム群９２２、ファイル群９２３が記憶されている。
【０１７０】
プログラム群９２２には、実施の形態において「〜部」として説明する機能を実行するプログラムが含まれる。プログラム（例えば、寸法推定プログラム）は、ＣＰＵ９０１により読み出され実行される。すなわち、プログラムは、「〜部」としてコンピュータを機能させるものであり、また「〜部」の手順や方法をコンピュータに実行させるものである。
【０１７１】
ファイル群９２３には、実施の形態において説明する「〜部」で使用される各種データ（入力、出力、判定結果、計算結果、処理結果など）が含まれる。
【０１７２】
実施の形態において構成図およびフローチャートに含まれている矢印は主としてデータや信号の入出力を示す。
フローチャートなどに基づいて説明する実施の形態の処理はＣＰＵ９０１、記憶装置、入力装置、出力装置などのハードウェアを用いて実行される。
【０１７３】
実施の形態において「〜部」として説明するものは「〜回路」、「〜装置」、「〜機器」であってもよく、また「〜ステップ」、「〜手順」、「〜処理」であってもよい。すなわち、「〜部」として説明するものは、ファームウェア、ソフトウェア、ハードウェアまたはこれらの組み合わせのいずれで実装されても構わない。
【０１７４】
実施の形態１により、被写体１０２を斜めの方向から映している一枚の写真画像１０１を用いて、被写体１０２の正しい寸法を計測することができる。
例えば、海上を撮像した写真画像１０１に被写体として不審船が映っている場合、実施の形態１により不審船の寸法を計測することができる。この場合、不審船に取り付けられている浮き輪をリファレンスとして使用することができる。船に搭載される浮き輪は一般的に規格化され、大きさが決まっているためである。
さらに、各国で使用されている船の寸法を含んだデータベースを用意しておき、計測した不審船の寸法とデータベースに記録された寸法とを比較することにより、不審船がどの国の船であるかを特定することもできる。
但し、船は被写体１０２の一例であり、浮き輪はリファレンスの一例である。つまり、実施の形態１によって船以外の被写体１０２の寸法を計測してもよいし、浮き輪以外のリファレンスを利用してもよい。
【０１７５】
また、被写体１０２（例えば、本）の各部分の寸法を算出することにより、被写体１０２の三次元モデルを生成することができる（図３５の（１）参照）。
さらに、被写体１０２の三次元モデルを用いて、異なる方向から見た被写体１０２を表示することができる（図３５の（２）から（５）参照）。
なお、被写体１０２の三次元モデルに写真画像１０１の画素をテクスチャとして張り付けることにより、被写体１０２の三次元モデルをより鮮明に描くことができる。
図３５は、実施の形態１における被写体１０２の三次元モデルを表す図である。
【０１７６】
実施の形態１において、例えば、以下のような寸法推定装置（１００）について説明した。括弧内に対応する構成の符号または名称を記す。
前記寸法推定装置は、基準線分情報生成部（１２０）と、基準軸ベクトル算出部（１４０）と、リファレンス線分ベクトル算出部（１５３）と、計測線分ベクトル算出部（１６３）と、計測線分長算出部（１６４）とを備える。
前記基準線分情報生成部は、対象物（１０２）を斜め方向から撮像した対象画像（１０１）に映っている前記対象物の幅方向と高さ方向と奥行き方向とに沿う３つの基準線分を前記対象画像の二次元の座標値を用いて表す基準線分情報（１２９）を生成する。
前記基準軸ベクトル算出部は、前記基準線分情報生成部によって生成された前記基準線分情報と、前記対象画像の撮像パラメータ（撮影パラメータ）とに基づいて、前記対象物の前記幅方向と前記高さ方向と前記奥行き方向とを三次元のベクトル値を用いて表す３つの基準軸ベクトル（１４９）を算出する。
前記リファレンス線分ベクトル算出部は、前記基準軸ベクトル算出部によって算出された前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルを用いて、前記対象画像に映っている前記対象物のうち長さが既知の線分であるリファレンス線分のベクトル値を、リファレンス線分ベクトル（１５９ａ）として算出する。
前記計測線分ベクトル算出部は、前記基準軸ベクトル算出部によって算出された前記３つの基準軸ベクトルのうち少なくともいずれかの基準軸ベクトルを用いて、前記対象画像に映っている前記対象物のうち長さが未知の線分である計測線分のベクトル値を、計測線分ベクトルとして算出する。
前記計測線分長算出部は、前記リファレンス線分の既知の長さ（１５９ｂ）と、前記リファレンス線分ベクトル算出部によって算出された前記リファレンス線分ベクトルと、前記計測線分ベクトル算出部によって算出された前記計測線分ベクトルとを用いて、前記計測線分の長さ（１６９）を算出する。
【符号の説明】
【０１７７】
１００  画像寸法推定装置、１０１  写真画像、１０２  被写体、１０９  入力機器、１０９ａ  マウスカーソル、１１０  写真画像表示部、１２０  姿勢基準軸設定部、１２９  姿勢基準軸データ、１３０  撮影パラメータ設定部、１３９  撮影パラメータ、１４０  基準軸ベクトル推定部、１４９  基準軸ベクトル、１５０  リファレンス設定部、１５１  基準軸ベクトル投影部、１５２  リファレンス線分設定部、１５３  リファレンス線分ベクトル算出部、１５４  リファレンス線分長取得部、１５９ａ  リファレンス線分ベクトル、１５９ｂ  リファレンス線分長、１６０  寸法推定部、１６１  基準軸ベクトル投影部、１６２  計測線分設定部、１６３  計測線分ベクトル算出部、１６４  計測線分長算出部、１６９  計測線分長、１７０  推定結果出力部、１８０  写真画像データベース、１９０  装置記憶部、９０１  ＣＰＵ、９０２  バス、９０３  ＲＯＭ、９０４  ＲＡＭ、９０５  通信ボード、９１１  ディスプレイ、９１２  キーボード、９１３  マウス、９１４  ドライブ、９２０  磁気ディスク装置、９２１  ＯＳ、９２２  プログラム群、９２３  ファイル群。
【図１】
2014109440.tif 000040
【図２】
2014109440.tif 000041
【図３】
2014109440.tif 000042
【図４】
2014109440.tif 000043
【図５】
2014109440.tif 000044
【図６】
2014109440.tif 000045
【図７】
2014109440.tif 000046
【図８】
2014109440.tif 000047
【図９】
2014109440.tif 000048
【図１０】
2014109440.tif 000049
【図１１】
2014109440.tif 000050
【図１２】
2014109440.tif 000051
【図１３】
2014109440.tif 000052
【図１４】
2014109440.tif 000053
【図１５】
2014109440.tif 000054
【図１６】
2014109440.tif 000055
【図１７】
2014109440.tif 000056
【図１８】
2014109440.tif 000057
【図１９】
2014109440.tif 000058
【図２０】
2014109440.tif 000059
【図２１】
2014109440.tif 000060
【図２２】
2014109440.tif 000061
【図２３】
2014109440.tif 000062
【図２４】
2014109440.tif 000063
【図２５】
2014109440.tif 000064
【図２６】
2014109440.tif 000065
【図２７】
2014109440.tif 000066
【図２８】
2014109440.tif 000067
【図２９】
2014109440.tif 000068
【図３０】
2014109440.tif 000069
【図３１】
2014109440.tif 000070
【図３２】
2014109440.tif 000071
【図３３】
2014109440.tif 000072
【図３４】
2014109440.tif 000073
【図３５】
2014109440.tif 000074
【図３６】
2014109440.tif 000075