"ステレオ画像間の誤差関数と奥行き方向への距離又は視差である距離指標の連続性を示すスムーズ関数とで規定される確率に関するエネルギーを最小化する手法により、前記ステレオ画像の一方である基準画像の画素ごとの前記距離指標を示す距離指標情報を推定すると共に、前記距離指標情報の信頼度を示す評価値を計算する距離指標情報推定装置であって、前記ステレオ画像が入力されると共に、予め設定された距離指標候補ごとに、予め設定された誤差関数閾値以下となるように、入力された前記ステレオ画像の類似度である前記誤差関数を計算する誤差関数計算部と、前記距離指標候補ごとに、予め設定されたスムーズ関数閾値以下となるように、マルコフランダムフィールドにおける隣接画素間の前記距離指標候補の差分絶対値である前記スムーズ関数を計算するスムーズ関数計算部と、前記誤差関数計算部で計算した誤差関数と前記スムーズ関数計算部で計算したスムーズ関数とを含む前記エネルギーが集合したエネルギー集合を生成し、当該エネルギー集合内における隣接画素間で、所定の漸化式によって前記エネルギーを生成して伝搬させることで、前記距離指標候補及び前記画素ごとの前記エネルギーを更新するエネルギー伝搬部と、前記エネルギー伝搬部で更新されたエネルギーを示す評価関数を計算する評価関数計算部と、前記距離指標候補のうち、前記評価関数計算部で計算された評価関数が最小となる距離指標を、前記距離指標情報として推定する距離指標推定部と、前記評価関数におけるエネルギー分布の尖度を、前記評価値として求める評価値計算部と、を備えることを特徴とする距離指標情報推定装置。前記評価値計算部は、前記エネルギー分布の尖度を、テクスチャに起因する推定エラーを判定するための前記評価値であるテクスチャ評価値として計算するテクスチャ評価値計算部、をさらに備えることを特徴とする請求項1に記載の距離指標情報推定装置。前記テクスチャ評価値計算部は、nfが前記距離指標候補fqの数、bqfqが前記エネルギー、bq ̄が前記エネルギーの平均値、及び、σbqが前記エネルギーの標準偏差と定義された以下の式9を用いて、前記テクスチャ評価値Kqを計算することを特徴とする請求項2に記載の距離指標情報推定装置。前記テクスチャ評価値計算部で計算されたテクスチャ評価値が予め設定されたテクスチャ閾値以下であるか否かを判定し、当該テクスチャ評価値が当該テクスチャ閾値以下の場合、前記テクスチャに起因する推定エラーが生じたと判定するテクスチャ推定エラー判定部、をさらに備えることを特徴とする請求項2又は請求項3に記載の距離指標情報推定装置。前記ステレオ画像を撮影したカメラの光軸交点をカメラキャリブレーションにより計算し、前記ステレオ画像内の被写体が収まるように予め設定された半径の仮想球が当該光軸交点を中心に生成され、前記中央カメラから前記仮想球の奥側表面までの前記距離指標を計算し、計算した当該距離指標を、前記中央カメラの撮影画像で前記テクスチャに起因する推定エラーが生じた領域内の各画素の前記距離指標情報として割り当てるテクスチャ推定エラー抑制部、をさらに備えることを特徴とする請求項2から請求項4の何れか一項に記載の距離指標情報推定装置。前記誤差関数計算部は、前記基準画像における処理対象画素と、前記ステレオ画像の他方である参照画像で前記処理対象画素に対応する対応画素との色情報の差分絶対値を求め、当該差分絶対値を予め設定した原色数で除算して重み付ける前記誤差関数を計算することを特徴とする請求項1から請求項5の何れか一項に記載の距離指標情報推定装置。ステレオ画像間の誤差関数と奥行き方向への距離又は視差である距離指標の連続性を示すスムーズ関数とで規定される確率に関するエネルギーを最小化する手法により、前記ステレオ画像の一方である基準画像の画素ごとの前記距離指標を示す距離指標情報を推定すると共に、前記距離指標情報の信頼度を示す評価値を計算するために、コンピュータを、前記ステレオ画像が入力されると共に、予め設定された距離指標候補ごとに、予め設定された誤差関数閾値以下となるように、入力された前記ステレオ画像の類似度である前記誤差関数を計算する誤差関数計算部、前記距離指標候補ごとに、予め設定されたスムーズ関数閾値以下となるように、マルコフランダムフィールドにおける隣接画素間での前記距離指標候補の差分絶対値である前記スムーズ関数を計算するスムーズ関数計算部、前記誤差関数計算部で計算した誤差関数と前記スムーズ関数計算部で計算したスムーズ関数とを含む前記エネルギーが集合したエネルギー集合を生成し、当該エネルギー集合内における隣接画素間で、所定の漸化式によって前記エネルギーを生成して伝搬させることで、前記距離指標候補及び前記画素ごとの前記エネルギーを更新するエネルギー伝搬部、前記エネルギー伝搬部で更新されたエネルギーを示す評価関数を計算する評価関数計算部、前記距離指標候補のうち、前記評価関数計算部で計算された評価関数が最小となる距離指標を、前記距離指標情報として推定する距離指標推定部、前記評価関数で表されるエネルギー分布の尖度を、前記評価値として求める評価値計算部、として機能させるための距離情報推定プログラム。しかし、従来の距離推定技術では、ステレオ画像にオクルージョンが生じた領域やテクスチャが無い領域が含まれる場合、これらに起因して、距離情報の推定結果にエラーが生じることが多い。この場合、精度が高い距離情報を得るために、これら推定エラーが生じた領域に対して、これら推定エラーの抑制処理を施す必要がある。この推定エラーの抑制処理の可否を判定するために、推定された距離情報の客観的な信頼度を提示することが強く要望されている。そこで、本発明は、前記した問題を解決し、距離指標情報の客観的な信頼度である評価値を提示できる距離指標情報推定装置及びそのプログラムを提供することを課題とする。前記した課題に鑑みて、本願第1発明に係る距離指標情報推定装置は、ステレオ画像間の誤差関数と奥行き方向への距離又は視差である距離指標の連続性を示すスムーズ関数とで規定される確率に関するエネルギーを最小化する手法により、ステレオ画像の一方である基準画像の画素ごとの距離指標を示す距離指標情報を推定すると共に、距離指標情報の信頼度を示す評価値を計算する距離指標情報推定装置であって、誤差関数計算部と、スムーズ関数計算部と、エネルギー伝搬部と、評価関数計算部と、距離指標推定部と、評価値計算部とを備えることを特徴とする。かかる構成によれば、距離指標情報推定装置は、誤差関数計算部によって、ステレオ画像が入力されると共に、予め設定された距離指標候補ごとに、予め設定された誤差関数閾値以下となるように、入力されたステレオ画像の類似度である誤差関数を計算する。つまり、誤差関数計算部は、誤差関数の値が誤差関数閾値を超えないようにトランケーションを行う。また、距離指標情報推定装置は、スムーズ関数計算部によって、距離指標候補ごとに、予め設定されたスムーズ関数閾値以下となるように、マルコフランダムフィールドにおける隣接画素間の距離指標候補の差分絶対値であるスムーズ関数を計算する。つまり、スムーズ関数計算部は、スムーズ関数の値がスムーズ関数閾値を超えないようにトランケーションを行う。また、距離指標情報推定装置は、エネルギー伝搬部によって、誤差関数計算部で計算した誤差関数とスムーズ関数計算部で計算したスムーズ関数とを含むエネルギーが集合したエネルギー集合を生成し、当該エネルギー集合内における隣接画素間で、所定の漸化式によってエネルギーを生成して伝搬させることで、距離指標候補及び画素ごとのエネルギーを更新する。また、距離指標情報推定装置は、評価関数計算部によって、エネルギー伝搬部で更新されたエネルギーを示す評価関数を計算する。そして、距離指標情報推定装置は、距離指標推定部によって、距離指標候補のうち、評価関数計算部で計算された評価関数が最小となる距離指標を、距離指標情報として推定する。ここで、前記したトランケーションでエネルギーに上限が設けられるため、評価関数のばらつきが一定範囲に収まって、距離指標情報の精度が向上する。また、距離指標情報推定装置は、評価値計算部によって、評価関数で表されるエネルギー分布の尖度を、評価値として求める。ここで、評価値は、エネルギー分布の尖り具合である尖度に基づいて計算されるため、その精度が向上し、距離情報の客観的な信頼度を示す情報となる。なお、本発明において、距離指標とは、奥行き方向への距離奥行き、又は、この距離の指標になる視差のことである。つまり、距離指標情報推定装置は、距離指標情報として、画素ごとの距離を示す距離情報、又は、画素ごとの視差を示す視差情報を推定する。また、本願第2発明に係る距離指標情報推定装置は、評価値計算部が、エネルギー分布の尖度を、テクスチャに起因する推定エラーを判定するためのテクスチャ評価値として計算する。このとき、テクスチャ評価値計算部は、nfが距離指標候補fqの数、bqfqがエネルギー、bq ̄がエネルギーの平均値、及び、σbqがエネルギーの標準偏差と定義された式9を用いて、テクスチャ評価値Kqを計算することが好ましい本願第3発明。かかる構成によれば、距離指標情報推定装置は、エネルギー分布の尖り具合である尖度を考慮してテクスチャ評価値を計算するため、テクスチャ評価値の精度を向上させることができる。また、本願第4発明に係る距離指標情報推定装置は、テクスチャ評価値計算部で計算されたテクスチャ評価値が予め設定されたテクスチャ閾値以下であるか否かを判定し、テクスチャ評価値がテクスチャ閾値以下の場合、テクスチャに起因する推定エラーが生じたと判定するテクスチャ推定エラー判定部、をさらに備えることを特徴とする。かかる構成によれば、距離指標情報推定装置は、テクスチャ評価値がテクスチャに起因する推定エラーが生じた画素で小さくなる性質を利用して、テクスチャに起因する推定エラーを判定する。また、本願第5発明に係る距離指標情報推定装置は、ステレオ画像を撮影したカメラの光軸交点をカメラキャリブレーションにより計算し、ステレオ画像内の被写体が収まるように予め設定された半径の仮想球が当該光軸交点を中心に生成され、中央カメラから仮想球の奥側表面までの距離指標を計算し、計算した距離指標を、中央カメラの撮影画像でテクスチャに起因する推定エラーが生じた領域内の各画素の距離指標情報として割り当てるテクスチャ推定エラー抑制部、をさらに備えることを特徴とする。かかる構成によれば、距離指標情報推定装置は、テクスチャに起因する推定エラーが生じたと判定された領域に対して、その推定エラーを抑制する処理を行う。また、本願第6発明に係る距離指標情報推定装置は、誤差関数計算部が、基準画像における処理対象画素と、ステレオ画像の他方である参照画像で処理対象画素に対応する対応画素との色情報の差分絶対値を求め、当該差分絶対値を予め設定した原色数で除算して重み付ける誤差関数を計算することを特徴とする。かかる構成によれば、距離指標情報推定装置は、被写体の輪郭を鮮鋭化して、精度が高い距離情報を推定することができる。なお、本願第1発明に係る距離指標情報推定装置は、一般的なコンピュータを、誤差関数計算部、スムーズ関数計算部、エネルギー伝搬部、評価関数計算部、距離指標推定部及び評価値計算部として動作させる距離指標情報推定プログラムによって実現することもできる。このプログラムは、通信回線を介して配布しても良く、CD-ROMやフラッシュメモリ等の記録媒体に書き込んで配布しても良い。本発明は、以下のような優れた効果を奏する。本願第1発明によれば、前記したトランケーションでエネルギーに上限が設けられるため、評価関数のばらつきが一定範囲に収まって、距離指標情報の精度を向上させることができる。さらに、本願第1発明によれば、エネルギー分布の尖り具合である尖度を考慮して評価値を計算するため、評価値の精度を向上させ、距離指標情報の客観的な信頼度として評価値を提示することができる。本願第2,3発明によれば、エネルギー分布の尖度を考慮してテクスチャ評価値を計算するため、テクスチャ評価値の精度を向上させることができる。本願第4発明によれば、テクスチャ評価値がテクスチャに起因する推定エラーが生じた画素で小さくなる性質を利用して、テクスチャに起因する推定エラーを正確に判定することができる。本願第5発明によれば、テクスチャに起因する推定エラーが生じたと判定された領域に対して、その推定エラーを抑制する処理を正確に行うことができる。本願第6発明によれば、被写体の輪郭を鮮鋭化して、精度が高い距離指標情報を推定することができる。本発明の概略:信頼度伝搬法による距離推定本発明の各実施形態では、ステレオ画像間の誤差関数と距離指標の連続性を示すスムーズ関数とで規定されるエネルギーを最小化する手法として、信頼度伝搬法を用いている。そこで、最初に、本発明の実施形態における信頼度伝搬法の概略を説明する。なお、以後の各実施形態では、請求項に記載の距離指標が、奥行き方向への距離奥行きであることとして説明する。つまり、請求項に記載の距離指標情報として、画素ごとの距離を示す距離情報を推定することになる。この信頼度伝搬法BP法:BeliefPropagation法は、距離情報にマルコフランダムフィールドMRF:MarkovRandomFieldモデルを仮定し、ステレオ画像間の誤差関数とスムーズ関数とで構成規定されるエネルギーを最小化することにより、ノイズが少なく滑らかな距離情報を生成するものである。ここで、信頼度伝搬法は、図1に示すように、処理対象となる画素pの隣接画素sから、メッセージと呼ばれる各画素にどの距離が割り当てられるかという確率に関するエネルギーmを受け取り、それらからメッセージを更新して隣接画素qに伝搬させる。そして、信頼度伝搬法は、この処理を画像内の全ての画素に対して繰り返し行うことで各画素の距離情報を求める。なお、請求項に記載のエネルギーが、信頼度伝搬法ではメッセージとなる。ここで、図1では、枠白塗四角形が基準画像の画素に対応し、矢印mが画素間のメッセージエネルギーの伝搬を示している。また、図1に示すように、このメッセージを伝搬する画素がノードとして縦横に配列されたグラフをエネルギー集合メッセージ集合と呼ぶ。さらに、文字sが付された画素sは、メッセージの伝搬元となる隣接画素を示し、文字pが付された画素pは、処理対象の画素を示し、文字qが付された画素qは、メッセージの伝搬先となる隣接画素を示す。なお、図1では、説明を簡易にするために、一部符号のみを図示した。図2を参照し、信頼度伝搬法による距離推定の手順を説明する。図2に示すように、信頼度伝搬法では、メッセージの生成に必要な画素間の誤差関数Dpfp以後、誤差関数Dpを計算するステップS1。本発明において、誤差関数に、SSD、SAD、ZNCC等の類似度が利用可能であり、ここでは、ステレオ画像の対応画素間で色情報の差分絶対値により類似度を示す式1を利用する。この式1によれば、画素単位で類似度を計算するため、被写体の輪郭を鮮鋭化することができる。ここで、pは処理対象の画素、fpは画素pの距離候補、cは色情報、r,g,bはRGB値、dpは距離候補fpに対応する視差、Icはステレオ画像の一方である基準画像での画素値、Ic ́はステレオ画像の他方である参照画像での画素値、λdataは重み係数を表わす。なお、距離候補距離指標候補とは、任意の範囲内で予め設定された距離の候補のことである図4等参照。すなわち、式1は、基準画像Icにおける処理対象の画素pの色情報Icpと、参照画像Ic ́における対応画素pdpの色情報Ic ́pdpとの差分絶対値を予め設定された原色数例えば、3で除算して、重み係数λdataで重み付けることを示す。この対応画素pdpとは、参照画像Ic ́において、処理対象の画素pに対応する画素のことである。また、式1において、視差の逆数にステレオカメラの間隔ベースラインを乗算して距離が計算できることを利用して、距離候補fpから視差dpを逆算できる。このベースラインは、例えば、OpenCV等のカメラキャリブレーションで求めることができるURLhttpopencvjp。続いて、信頼度伝搬法では、下記の式2に示すように、距離候補fpと、画素qの距離候補fqとの差分絶対値|fp-fq|で定義されるスムーズ関数Vfp-fq以後、スムーズ関数Vの計算を行うステップS2。このスムーズ関数Vは、距離情報の滑らかさ例えば、ノイズの少なさ、言い換えるなら、距離の連続性を示す。ここで、qは画素pのメッセージ伝搬先となる隣接画素である。そして、信頼度伝搬法では、下記の式3で規定される生成式によって、メッセージを生成して伝搬するステップS3。mはメッセージ、tは反復回数、Np\qは画素pへメッセージを渡す画素q以外の4近傍の画素集合、sはその画素集合の要素画素、minは最小値を返す関数である。この式3は、漸化式になっており、受け取ったメッセージをもとにメッセージを更新する処理を繰り返し行うことを示す。例えば、信頼度伝搬法では、ラスタスキャンを行うように処理対象画素pを移動させながら、その処理対象画素pでメッセージを更新することで、メッセージの伝搬を全画素で行う。全画素についてメッセージの更新処理を、予め設定したT回例えば、10回まで反復したとする。この場合、画素qの距離候補fqに関する評価関数bqfq以後、評価関数bqは、下記の式4で表すことができる。そして、信頼度伝搬法では、この評価関数bqが最小になる距離候補fqを、各画素の最終的な距離として推定するステップS4。第1実施形態[距離情報推定装置の構成]以下、本発明の各実施形態について、適宜図面を参照しながら詳細に説明する。なお、各実施形態において、同一の機能を有する手段には同一の符号を付し、説明を省略した。図3を参照し、本発明の第1実施形態に係る距離情報推定装置距離指標情報推定装置1の構成について説明する。距離情報推定装置1は、信頼度伝搬法により距離情報を推定して、この距離情報の信頼度を示す評価値を計算すると共に、この評価値によって、距離情報の推定エラーが生じたか否かを判定するものである。図3に示すように、距離情報推定装置1は、誤差関数計算部10と、スムーズ関数計算部11と、メッセージ生成伝搬部エネルギー伝搬部12と、評価関数計算部13と、距離推定部距離指標推定部14と、評価値計算部15と、オクルージョン領域判定部オクルージョン推定エラー判定部17と、無テクスチャ領域判定部テクスチャ推定エラー判定部19とを備える。本実施形態では、距離情報推定装置1は、評価値として、オクルージョンに起因する推定エラーを判定するための評価値Oqオクルージョン評価値、以後評価値Oと、テクスチャに起因する推定エラーを判定するための評価値Kqテクスチャ評価値、以後評価値Kとを計算する。そして、距離情報推定装置1は、評価値O及び評価値Kを用いて、距離情報の推定エラーとして、オクルージョンに起因する推定エラーと、テクスチャに起因する推定エラーとを判定することとする。テクスチャに起因する推定エラーとは、例えば、ステレオ画像でのテクスチャが存在しない領域テクスチャレス領域において、マッチング精度が低下して、距離情報の推定エラーが生じることである。距離情報推定装置1は、左カメラCL及び右カメラCRステレオカメラから、被写体不図示を撮影したステレオ画像基準画像及び参照画像が入力される。この左カメラCL及び右カメラCRは、静止画や動画を撮影可能な撮影カメラであり、図3に示すように、一定間隔ベースラインだけ離して配置される。ここで、例えば、左カメラCL及び右カメラCRのうち左カメラCLを基準カメラとし、左カメラCLの撮影画像を基準画像とし、右カメラCRの撮影画像を参照画像とする。誤差関数計算部10は、ステレオ画像が入力されると共に、距離候補ごとに、予め設定された誤差関数閾値Tdata以下となるように、入力されたステレオ画像の類似度である誤差関数Dpを計算するものである。この誤差関数計算部10は、下記の式5及び式6を用いて、図4に示すように画素s,p,qについて、距離候補ごとの誤差関数Dpを計算する。ここで、ifは、後段の条件式を満たす場合、前段に記述した値を出力する関数である。本実施形態では、距離候補fs,fp,fqがそれぞれ、図4に示すように、0,1,,K,,N-1の範囲内で予め設定されたこととして説明する但し、K,Nは、K<N-1を満たす整数。なお、図4では、文字s,p,qを付した四辺形が画素s,p,qである。このとき、誤差関数計算部10は、式6に示すように、誤差関数Dpの値が誤差関数閾値Tdataを超えないように、トランケーションを行っている。具体的には、誤差関数計算部10は、誤差関数Dpの値が誤差関数閾値Tdataを超える場合、誤差関数閾値Tdataを誤差関数Dpの値としてメッセージ生成伝搬部12に出力する。その一方、誤差関数計算部10は、誤差関数Dpの値が誤差関数閾値Tdata以下の場合、メッセージ生成伝搬部12に誤差関数Dpの値をそのまま出力する。スムーズ関数計算部11は、距離候補ごとに、予め設定されたスムーズ関数閾値Tsmooth以下となるように、スムーズ関数Vを計算するものである。このスムーズ関数計算部11は、図4に示すように、画素pの距離候補fp0~N-1と、画素qの距離候補fq0~N-1との全組み合わせについて、下記の式7を用いて、スムーズ関数Vを計算する。このとき、スムーズ関数計算部11は、式7に示すように、スムーズ関数Vの値がスムーズ関数閾値Tsmoothを超えないように、トランケーションを行っている。具体的には、スムーズ関数計算部11は、スムーズ関数Vの値がスムーズ関数閾値Tsmoothを超える場合、スムーズ関数閾値Tsmoothをスムーズ関数Vの値としてメッセージ生成伝搬部12に出力する。その一方、スムーズ関数計算部11は、スムーズ関数Vの値がスムーズ関数閾値Tsmooth以下の場合、メッセージ生成伝搬部12にスムーズ関数Vの値をそのまま出力する。ここで、マルコフランダムフィールドおよびスムーズ関数Vについて補足する。このマルコフランダムフィールドは、基準画像で互いに隣接する画素の色、輝度や距離が似ていることを示すモデルである。そして、スムーズ関数Vは、このマルコフランダムフィールドを踏まえた距離の連続性に関する関数であり、基準画像の各画素に距離候補を割り当てた際、互いに隣接する画素の距離候補が同じであればエネルギーが小さく、異なればエネルギーが大きくなる。例えば、マルコフランダムフィールドでは、ある画素の距離候補が1で、その隣接画素の距離候補が1である場合、エネルギーは0となる。また、ある画素の距離候補が1で、その画素に隣接する画素の距離候補が2である場合、例えば、エネルギーは1となる。さらに、ある画素の距離候補が1で、その隣接画素の距離候補が10である場合、例えば、エネルギーは9となる。つまり、信頼度伝搬法において、エネルギーが最小となる距離候補が最終的な距離となることを考慮すれば、このマルコフランダムフィールドにより、互いに隣接する画素は、同じような距離になり易くなる。メッセージ生成伝搬部エネルギー伝搬部12は、誤差関数計算部10で計算された誤差関数Dpと、スムーズ関数計算部11で計算されたスムーズ関数Vとが入力される。そして、メッセージ生成伝搬部12は、図1と同様に、これら誤差関数Dpおよびスムーズ関数Vを含むメッセージが集合したエネルギー集合を生成するものである。具体的には、メッセージ生成伝搬部12は、基準画像の画素位置に対応した位置に画素ごとのメッセージを配置することで、エネルギー集合をモデル化生成する。次に、メッセージ生成伝搬部12は、前記した式3を用いて、エネルギー集合内における隣接画素間で、メッセージを生成して伝搬することで、画素ごとのメッセージを更新する。まず、メッセージ生成伝搬部12は、図5,図6に示すように、画素sにおいて、距離候補fs0~N-1のメッセージms0~msN-1を生成して、画素sから画素pに伝搬させる。次に、メッセージ生成伝搬部12は、画素pにおいて、距離候補fp0~N-1のメッセージmp→q0~mp→qN-1を生成して、画素pから画素qに伝搬させる。評価関数計算部13は、前記した式4を用いて、メッセージ生成伝搬部12で更新されたメッセージを示す評価関数bqを計算するものである。そして、評価関数計算部13は、計算した評価関数bqを、距離推定部14と、評価値O計算部16と、評価値K計算部18とに出力する。距離推定部距離指標推定部14は、評価関数計算部13で計算された評価関数bqが入力される。そして、距離推定部14は、距離候補fqのうち、この評価関数bqが最小となる距離候補fqを、距離情報として推定するものである。すなわち、距離推定部14は、前記した式4の評価関数bqが最小となる距離候補fqを、基準画像の画素qの最終的な距離情報として推定する。ここでは、距離推定部14は、距離情報として、画素の距離が手前側である程、その画素の輝度が高く、画素の距離が奥側である程、その画素の輝度が低くなる距離画像を生成して出力する。この距離画像は、例えば、被写体の3次元立体モデルの生成に利用できる。評価値計算部15は、評価関数bqにおけるエネルギー分布の形状により、評価値を求めるものであり、評価値O計算部オクルージョン評価値計算部16と、評価値K計算部テクスチャ評価値計算部18とを備える。評価値O計算部オクルージョン評価値計算部16は、評価関数計算部13で計算された評価関数bqが入力される。この評価関数bqは、基準画像の各画素について、距離候補fqごとのメッセージを示す情報、つまり、距離候補fqとメッセージとを対応付けたエネルギー分布メッセージ分布と言える。具体的には、評価値O計算部16は、メッセージが最も高くなる最大値とメッセージが最も低くなる最小値との差分をこの最大値で除算することで、評価値Oを計算するものである。具体的には、評価値O計算部16は、下記の式8を用いて評価値Oを計算して、オクルージョン領域判定部17に出力する。ここで、maxは、最大値を返す関数である。従って、maxbqfqは、画素qにおけるメッセージの最大値を示す。また、minbqfqは、画素qにおけるメッセージの最小値を示す。オクルージョン領域判定部オクルージョン推定エラー判定部17は、評価値O計算部16から評価値Oが入力される。そして、オクルージョン領域判定部17は、この評価値Oが予め設定されたオクルージョン閾値To以下であるか否かを判定し、評価値Oがオクルージョン閾値To以下の場合、オクルージョンに起因する推定エラーが生じたと判定するものである。ここで、オクルージョンに起因する推定エラーが生じていない画素では、図7aのエネルギー分布のように、メッセージの最大値maxと最小値minとの差分ΔEOが、大きくなる。一方、オクルージョンに起因する推定エラーが生じた画素では、図7bのエネルギー分布のように、差分ΔEOが図7aよりも小さくなる。また、図7の差分ΔEOは、式8の分子であることから、評価値Oに比例する。つまり、オクルージョン領域判定部17は、評価値Oが、オクルージョンに起因する推定エラーが生じた画素で、オクルージョンに起因する推定エラーが生じていない画素よりも小さくなる性質を利用して、オクルージョンに起因する推定エラーを判定している。そして、オクルージョン領域判定部17は、評価値Oがオクルージョン閾値To以下と判定された画素領域をオクルージョン領域として出力する。このオクルージョン領域は、例えば、被写体の3次元立体モデルを生成する際、オクルージョンに起因する推定エラーが生じた領域を把握するために利用できる。なお、オクルージョン領域判定部17は、評価値Oが閾値Toを超える場合、オクルージョンに起因する推定エラーが生じていないと判定することは言うまでもない。評価値K計算部テクスチャ評価値計算部18は、評価関数計算部13で計算された評価関数bqが入力される。そして、評価値K計算部18は、基準画像の各画素について、このエネルギー分布の尖度を評価値Kとして計算するものである。具体的には、評価値K計算部18は、下記の式9に示すように、距離候補fqの個数nf例えば、0~N-1までの合計N個と、メッセージと、メッセージの平均値bq ̄と,メッセージの標準偏差σbqとに基づいて、評価値Kを計算する。そして、評価値K計算部18は、計算した評価値Kを、無テクスチャ領域判定部19に出力する。無テクスチャ領域判定部テクスチャ推定エラー判定部19は、評価値K計算部18から、評価値Kが入力される。そして、無テクスチャ領域判定部19は、この評価値Kが予め設定されたテクスチャ閾値TK以下であるか否かを判定し、評価値Kがテクスチャ閾値TK以下の場合、テクスチャに起因する推定エラーが生じたと判定するものである。ここで、テクスチャに起因する推定エラーが生じていない画素では、図8aに示すように、エネルギー分布が尖っており、その尖度を示す評価値Kが大きくなる。一方、テクスチャに起因する推定エラーが生じた画素では、図8bに示すように、エネルギー分布が尖っておらず、評価値Kが図8aよりも小さくなる。この尖度Kurtosisは、エネルギー分布の尖り具合を示すものであり、正規分布と比べて、以下のような特徴がある。すなわち、尖度が大きくなるほど、鋭いピークと長い尾とを有するエネルギー分布になる図8a参照。その一方、尖度が小さくなるほど、より丸みがかったピークと短い尾とを有するエネルギー分布になる図8b参照。なお、尾とは、図8において、トランケーションによって形成される、メッセージの上限となる平坦部である。つまり、無テクスチャ領域判定部19は、評価値Kが、テクスチャに起因する推定エラーが生じた画素で、テクスチャに起因する推定エラーが生じていない画素よりも小さくなる性質を利用して、テクスチャに起因する推定エラーを判定している。そして、無テクスチャ領域判定部19は、評価値Kがテクスチャ閾値TK以下と判定された画素領域を無テクスチャ領域として出力する。このオクルージョン領域は、例えば、被写体の3次元立体モデルを生成する際、テクスチャに起因する推定エラーが生じた領域を把握するために利用できる。なお、無テクスチャ領域判定部19は、評価値Kがテクスチャ閾値TKを超える場合、テクスチャに起因する推定エラーが生じていないと判定することは言うまでもない。[距離情報推定装置の動作]図9を参照し、距離情報推定装置1の動作について説明する適宜図3参照。距離情報推定装置1は、誤差関数計算部10によって、距離候補ごとに、誤差関数閾値Tdata以下となるように誤差関数Dpを計算する。つまり、誤差関数計算部10は、前記した式5及び式6を用いて、誤差関数Dpを計算して、トランケーションを行うステップS10。距離情報推定装置1は、スムーズ関数計算部11によって、距離候補ごとに、スムーズ関数閾値Tsmooth以下となるように、スムーズ関数Vを計算する。つまり、スムーズ関数計算部11は、前記した式7を用いて、スムーズ関数Vを計算して、トランケーションを行うステップS11。距離情報推定装置1は、メッセージ生成伝搬部12によって、図1と同様にエネルギー集合を生成する。そして、メッセージ生成伝搬部12は、前記した式3を用いて、エネルギー集合内における隣接画素間で、メッセージを生成して伝搬するステップS12。距離情報推定装置1は、評価関数計算部13によって、前記した式4を用いて、距離候補ごとに評価関数bqを計算するステップS13。そして、距離情報推定装置1は、距離推定部14によって、距離候補のうち、この評価関数bqが最小となる距離を、距離情報として推定するステップS14。距離情報推定装置1は、評価値O計算部16によって、評価値Oを計算する。つまり、評価値O計算部16は、前記した式8を用いて、メッセージの最大値とメッセージの最小値との差分をこの最大値で除算することで、評価値Oを計算するステップS15。距離情報推定装置1は、オクルージョン領域判定部17によって、評価値Oがオクルージョン閾値To以下であるか否かを判定し、評価値Oがオクルージョン閾値To以下の場合、オクルージョンに起因する推定エラーが生じたと判定する。ステップS16。距離情報推定装置1は、評価値K計算部18によって、エネルギー分布の尖度である評価値Kを計算する。つまり、評価値K計算部18は、前記した式9に示すように、距離候補の個数nfと、メッセージと、メッセージの平均値bq ̄と、メッセージの標準偏差σbqとに基づいて、評価値Kを計算するステップS17。距離情報推定装置1は、無テクスチャ領域判定部19によって、評価値Kがテクスチャ閾値TK以下であるか否かを判定し、評価値Kがテクスチャ閾値TK以下の場合、テクスチャに起因する推定エラーが生じたと判定するステップS18。以上のように、本発明の第1実施形態に係る距離情報推定装置1は、メッセージの最小値だけでなくメッセージの最大値も考慮して評価値Oを計算すると共に、エネルギー分布の尖度である評価値Kを計算する。そして、距離情報推定装置1は、オクルージョンやテクスチャに起因する推定エラーが生じた画素で、評価値O及び評価値Kが小さくなる性質を利用して、これら推定エラーを判定する。これによって、距離情報推定装置1は、オクルージョン及びテクスチャに起因する推定エラーの判定精度を向上させて、これら推定エラーが生じた領域を正確に提示することができる。さらに、距離情報推定装置1は、図7,図8に示すように、トランケーションによって、メッセージの上限が設けられるため、評価関数bqのばらつきが一定範囲に収まって、距離情報の精度を向上させることができる。第2実施形態[距離情報推定装置の構成]図10を参照し、本発明の第2実施形態に係る距離情報推定装置1Bの構成について、第1実施形態と異なる点を説明する。この距離情報推定装置1Bは、オクルージョン及びテクスチャに起因する推定エラーを抑制する点が、第1実施形態と異なる。このため、距離情報推定装置1Bは、テクスチャ推定エラー抑制部20をさらに備える。図10に示すように、距離情報推定装置1Bは、左カメラCL及び中央カメラCCの撮影画像である左ステレオ画像と、右カメラCR及び中央カメラCCの撮影画像である右ステレオ画像とが入力される。この中央カメラCCは、左カメラCL及び右カメラCRと同様、静止画や動画を撮影可能な撮影カメラである。また、左カメラCL及び右カメラCRは、それぞれ、中央カメラCCから一定間隔離して左右に配置される。ここでは、中央カメラCCを基準カメラとし、中央カメラCCの撮影画像を基準画像とし、左カメラCL及び右カメラCRの撮影画像を参照画像とする。誤差関数計算部10と、スムーズ関数計算部11と、メッセージ生成伝搬部12と、評価関数計算部13と、距離推定部14とは、左ステレオ画像及び右ステレオ画像それぞれの処理を行う以外、図3の各手段と同様であるため、詳細な説明を省略する。ここで、後記するオクルージョンに起因する推定エラーの抑制には、左ステレオ画像及び右ステレオ画像の両方以下、両ステレオ画像から求めた評価関数bqが必要になる。このため、評価関数計算部13は、両ステレオ画像から求めた評価関数bqを評価値O計算部16Bに出力する。さらに、距離推定部14は、両ステレオ画像から求めた距離情報距離画像を、評価値O計算部16Bに出力する。一方、後記するテクスチャに起因する推定エラーの抑制には、左ステレオ画像及び右ステレオ画像の何れか一方から求めた評価関数bqがあればよい。従って、評価関数計算部13は、左ステレオ画像から求めた評価関数bqを、評価値K計算部18に出力する。評価値O計算部16Bは、両ステレオ画像から評価値Oを計算して、オクルージョン領域判定部17に出力する。また、評価値O計算部16Bは、両ステレオ画像から計算した評価値Oと、両ステレオ画像の距離情報とを用いて、オクルージョンに起因する推定エラーを抑制する。なお、評価値Oの計算は、図3の評価値O計算部16と同様のため、説明を省略する。具体的には、評価値O計算部16Bは、基準画像の各画素について、左ステレオ画像から計算した評価値Oと、右ステレオ画像から計算した評価値Oとを比較する。この比較結果に基づいて、評価値O計算部16Bは、両ステレオ画像の距離情報のうち、評価値Oが大きい側の距離情報を、その画素の距離情報として、テクスチャ推定エラー抑制部20に出力する。すなわち、評価値O計算部16Bは、オクルージョンに起因する推定エラーが生じた画素で評価値Oが小さくなることを利用して、距離推定部14で生成された両ステレオ画像の距離情報を、この推定エラーを抑制した1つの距離情報に統合する。オクルージョン領域判定部17は、図3と同様、両ステレオ画像のオクルージョン領域を出力するため、詳細な説明を省略する。なお、このオクルージョン領域は、オクルージョンに起因する推定エラーの抑制に必要ではなく、参考情報として提示される。評価値K計算部18は、図3と同様、左ステレオ画像についての評価値Kを計算するため、詳細な説明を省略する。無テクスチャ領域判定部19は、図3と同様、左ステレオ画像についての無テクスチャ領域を生成して、テクスチャ推定エラー抑制部20に出力する。テクスチャ推定エラー抑制部20は、評価値O計算部16Bから距離情報が入力され、無テクスチャ領域判定部19から無テクスチャ領域が入力される。そして、テクスチャ推定エラー抑制部20は、この無テクスチャ領域を用いて、この距離情報から、テクスチャに起因する推定エラーを抑制するものである。図11,図12を参照して、テクスチャに起因する推定エラーの抑制について説明する適宜図10参照。説明を簡易にするため、図11に示すように、土俵Aの上で二人の被写体力士O1,O2が取組を行っている例をあげて説明する。この場合、図10の中央カメラCCは、土俵Aの正面中央に位置し、取組を行っている被写体O1,O2の側面方向図11の矢印方向から基準画像Pを撮影したこととする。なお、図11では、基準画像Pに含まれない土俵Aと被写体O1,O2との一部を破線で図示した。また、図11では、基準画像Pのうち、土俵Aの奥側地面となる背景Bの領域が無テクスチャ領域TXであり、ハッチングで図示した。この例では、テクスチャ推定エラー抑制部20は、図11のハッチング部分に対して、テクスチャに起因する推定エラーの抑制処理を行うことになる。ここで、中央カメラCCと、土俵Aと、被写体O1,O2とを上側より見ると、これらの位置関係は、図12のようになる。図12では、土俵A及び被写体O1,O2を破線で図示すると共に、後記する距離Dを一転鎖線で図示した。また、図11のハッチング部分無テクスチャ領域TXは、球表面の一部を切り取った曲面形状になる図12の吹き出し。まず、テクスチャ推定エラー抑制部20は、左カメラCLと、中央カメラCCと、右カメラCRとの光軸LAが交わる光軸交点CEを、カメラキャリブレーションにより計算する。ここで、テクスチャ推定エラー抑制部20は、全てのカメラで光軸交点CEを求める必要がなく、基準画像を撮影する中央カメラCCが含まれていればよい。つまり、テクスチャ推定エラー抑制部20は、中央カメラCCと左カメラCLとの組、又は、中央カメラCCと右カメラCRとの一方の組から、光軸交点CEを求めてもよい。次に、テクスチャ推定エラー抑制部20は、予め設定された半径rの仮想球VCを、光軸交点CEを中心に生成する。ここで、テクスチャ推定エラー抑制部20は、全ての被写体例えば、土俵A、被写体O1,O2が仮想球VCに含まれるように、半径rが予め設定されている。そして、テクスチャ推定エラー抑制部20は、中央カメラCCから仮想球VCの奥側表面までの距離Dを計算し、この距離Dを、無テクスチャ領域TX内の各画素の距離情報として割り当てる。すなわち、無テクスチャ領域TX内の各画素の距離情報は、下記の計算式を用いて割り当てることができる。まず、テクスチャ推定エラー抑制部20は、下記の10及び式11の連立方程式により、仮想球VCと中央カメラCCの光軸LAとの交点X,Y,Zを計算する。前記した式10は、仮想球VCを表した方程式である。ここで、a,b,cが仮想球VCの中心の世界座標光軸交点CEの世界座標であり、rが仮想球VCの半径である。前記した式11は、中央カメラCCの光学主点と無テクスチャ領域TX内の各画素とを通過する直線を表した方程式である。ここで、e,g,iが無テクスチャ領域TX内の各画の世界座標と光学主点の世界座標とにより計算される傾き既知であり、f,h,jが光学主点の世界座標既知である。そして、テクスチャ推定エラー抑制部20は、下記の12により、中央カメラCCから仮想球VCの奥側表面までの距離Dを計算する。ここで、zは仮想球VCの奥側表面までの距離Dであり、i,jは基準画像の座標既知であり、Aはカメラキャリブレーションにより計算される内部パラメータ既知であり、Rはカメラキャリブレーションにより計算される回転行列既知であり、Uはカメラキャリブレーションにより計算される並進行列既知である。[距離情報推定装置の動作]図13を参照し、距離情報推定装置1Bの動作について説明する適宜図10参照。ステップS20~S24の処理は、図9のステップS10~S14と同様のため、説明を省略する。距離情報推定装置1Bは、評価値O計算部16Bによって、無テクスチャ領域及び距離情報を用いて、オクルージョンに起因する推定エラーを抑制する。つまり、評価値O計算部16Bは、両ステレオ画像から評価値Oを計算し、これら評価値Oを比較して、評価値Oが大きい側のステレオ画像から生成した距離情報を、基準画像での各画素の距離情報として出力するステップS25。ステップS26~S28の処理は、図9のステップS16~S18と同様のため、説明を省略する。距離情報推定装置1Bは、テクスチャ推定エラー抑制部20によって、無テクスチャ領域を用いて、距離情報から、テクスチャに起因する推定エラーを抑制する。つまり、テクスチャ推定エラー抑制部20は、左カメラCLと、中央カメラCCと、右カメラCRとの光軸交点CEをカメラキャリブレーションにより計算する。そして、テクスチャ推定エラー抑制部20は、半径rの仮想球VCを、光軸交点CEを中心に生成する。さらに、テクスチャ推定エラー抑制部20は、中央カメラCCから仮想球VCの奥側表面までの距離Dを計算して、無テクスチャ領域TX内の各画素の距離情報として割り当てるステップS29。以上のように、本発明の第2実施形態に係る距離情報推定装置1Bは、オクルージョン及びテクスチャに起因する推定エラーを抑制して、精度が高い距離情報を生成することができる。この距離情報を用いれば、例えば、奥行き感が高い、インテグラルフォトグラフィIP方式の立体映像を生成することができる。第3実施形態[距離情報推定装置の構成]図14を参照し、本発明の第3実施形態に係る距離情報推定装置1Cの構成について、第1実施形態と異なる点を説明する。この距離情報推定装置1Cは、メッセージ伝搬の際、隣接画素間で色が大きく異なるときにその伝搬を制限する点が、第1実施形態と異なる。このため、距離情報推定装置1Cは、メッセージ伝搬制限判定部エネルギー伝搬判定部21をさらに備える。ここで、誤差関数計算部10と、スムーズ関数計算部11と、評価値計算部15と、オクルージョン領域判定部17と、無テクスチャ領域判定部19とは、図3の各手段と同様であるため、詳細な説明を省略する。メッセージ伝搬制限判定部エネルギー伝搬判定部21は、基準画像における隣接画素間の色情報の差分絶対値が予め設定された色情報閾値Tc以下であるか否かによって、隣接画素間でメッセージを伝搬できるか否かを判定するものである。このメッセージ伝搬制限判定部21は、メッセージの伝搬の制限を判定するために、隣接画素間の色情報の差分絶対値に関する色情報閾値Tcが予め設定される。このとき、色情報閾値Tcの値が小さくなってメッセージ伝搬の制限が強くなると、重み係数λdataの値も小さくなってスムーズ関数の影響が大きくなるように、色情報閾値Tc及び重み係数λdataを予め設定することが好ましい。また、メッセージ伝搬制限判定部21は、メッセージ伝搬元の画素sと処理対象の画素pとの色情報の差分絶対値が色情報閾値Tc以下であるか否かによって、メッセージ伝搬元の画素sと処理対象の画素pとの間でメッセージを伝搬できるか否かを判定する。具体的には、メッセージ伝搬制限判定部21は、下記の式13に示すように、画素s,pの色情報の差分絶対値と色情報閾値Tcとを比較することにより、メッセージ伝搬元の重み係数λs,pを計算する。ここで、メッセージ伝搬制限判定部21は、式13に示すように、画素s,pの色情報の差分絶対値が色情報閾値Tc以下の場合、画素s,p間でメッセージを伝搬できると判定し、重み係数λs,pを1として計算する。一方、メッセージ伝搬制限判定部21は、画素s,pの色情報の差分絶対値が色情報閾値Tcを超える場合、画素s,p間でメッセージを伝搬できないと判定し、重み係数λs,pを0として計算する。そして、メッセージ伝搬制限判定部21は、計算した重み係数λs,pをメッセージ生成伝搬部12Cに出力する。また、メッセージ伝搬制限判定部21は、メッセージ伝搬先の画素qと処理対象の画素pとの色情報の差分絶対値が色情報閾値Tc以下であるか否かによって、メッセージ伝搬先の画素qと処理対象の画素pとの間でメッセージを伝搬できるか否かを判定する。具体的には、メッセージ伝搬制限判定部21は、下記の式14に示すように、画素p,qの色情報の差分絶対値と色情報閾値Tcとを比較することにより、メッセージ伝搬先の重み係数λp,qを計算する。ここで、メッセージ伝搬制限判定部21は、式14に示すように、画素p,qの色情報の差分絶対値が色情報閾値Tc以下の場合、画素p,q間でメッセージを伝搬できると判定し、重み係数λp,qを1として計算する。一方、メッセージ伝搬制限判定部21は、画素p,qの色情報の差分絶対値が色情報閾値Tcを超える場合、画素p,q間でメッセージを伝搬できないと判定し、重み係数λp,qを0として計算する。そして、メッセージ伝搬制限判定部21は、計算した重み係数λp,qを評価関数計算部13Cに出力する。メッセージ生成伝搬部12Cは、メッセージ伝搬制限判定部21からメッセージ伝搬元の重み係数λs,pが入力される。そして、メッセージ生成伝搬部12Cは、下記の式15を用いて、メッセージを生成して伝搬する。このとき、メッセージ生成伝搬部12Cは、式15に示すように、画素p,qの色情報の差分絶対値と色情報閾値Tcとを比較して、メッセージを伝搬するか又は制限するかを判定する。具体的には、メッセージ生成伝搬部12Cは、画素p,qの色情報の差分絶対値が色情報閾値Tc以下の場合、画素pから画素qへのメッセージを伝搬させる。一方、メッセージ生成伝搬部12Cは、色情報の差分絶対値が色情報閾値Tcを超える場合、画素pから画素qへメッセージを伝搬させない。図15,図16を参照して、メッセージの伝搬の制限について詳細に説明する適宜図14参照。この図15,図16では、図1と同様のエネルギー集合において、枠白塗四角形が明るい色の画素であり、ハッチングされた枠が暗い色の画素である。ここでは、明るい色の画素と暗い色の画素との色情報の差分絶対値が色情報閾値Tcを超えることとして説明する。例えば、被写体と背景との境界といったように、隣接画素間で色が大きく異なる場合を考える。図15aに示すように、画素s,pで色が大きく異なる場合、距離情報推定装置1Cは、前記した式13よりメッセージ伝搬元の重み係数λs,pを0として、これら2画素s、pでのメッセージ伝搬を制限する。また、図15bに示すように、画素p,qで色が大きく異なる場合、距離情報推定装置1Cは、前記した式14よりメッセージ伝搬先の重み係数λp,qを0として、これら2画素p、qでのメッセージ伝搬を制限する。一方、図16aに示すように画素s,p,qの全てが明るい色、又は、図16bに示すように画素s,p,qの全てが暗い色の場合、距離情報推定装置1Cは、前記した式13及び式14より重み係数λs,p,λp,qを1として、画素s,p,q間でメッセージを伝搬する。つまり、距離情報推定装置1Cでは、メッセージ伝搬元の重み係数λs,pの値によって、画素sから画素pへのメッセージ伝搬を制御し、メッセージ伝搬先の重み係数λp,qの値によって、画素pから画素qへのメッセージ伝搬を制御する。評価関数計算部13Cは、メッセージ伝搬制限判定部21で計算された重み係数λp,qが入力され、下記の式16で定義される評価関数bqを計算する。そして、評価関数計算部13Cは、計算した評価関数bqを、距離推定部14Cと、評価値O計算部16と、評価値K計算部18とに出力する。距離推定部14Cは、図3の距離推定部14と同様、距離情報を推定する。また、距離推定部14Cは、ノイズ除去フィルタ14aを備える。このノイズ除去フィルタ14aは、距離推定部14Cが推定した距離情報距離画像にメディアンフィルタ等のノイズ除去処理を施すものである。このようにして、ノイズの影響が少ない滑らかな距離情報を推定することができる。[距離情報推定装置の動作]図17を参照し、距離情報推定装置1Cの動作について説明する適宜図14参照。ステップS30,S31の処理は、図9のステップS10~S11と同様のため、説明を省略する。距離情報推定装置1Cは、メッセージ伝搬制限判定部21によって、隣接画素間の色情報の差分絶対値が色情報閾値Tc以下であるか否かによって、隣接画素間でメッセージを伝搬できるか否かを判定する。具体的には、メッセージ伝搬制限判定部21は、前記した式13に示すように、画素s,pの色情報の差分絶対値と色情報閾値Tcとを比較することにより、画素s,p間におけるメッセージ伝搬の判定結果として、メッセージ伝搬元の重み係数λs,pを計算する。さらに、メッセージ伝搬制限判定部21は、前記した式14に示すように、画素p,qの色情報の差分絶対値と色情報閾値Tcとを比較することにより、画素p,q間におけるメッセージ伝搬の判定結果として、メッセージ伝搬先の重み係数λp,qを計算するステップS32。距離情報推定装置1Cは、メッセージ生成伝搬部12Cによって、前記した式15を用いて、メッセージを生成して伝搬するステップS33。距離情報推定装置1Cは、評価関数計算部13Cによって、前記した式16で定義される評価関数bqを計算するステップS34。なお、ステップS35~S39の処理は、図9のステップS14~S28と同様のため、説明を省略する。以上のように、本発明の第3実施形態に係る距離情報推定装置1Cは、隣接画素間の色情報の差分絶対値が大きいときにメッセージの伝搬を制限する。これによって、距離情報推定装置1Cは、被写体の膨張を抑制して被写体の輪郭を鮮鋭化し、精度が高い距離情報を生成することができる。なお、各実施形態では、誤差関数とスムーズ関数とで規定される確率に関するエネルギーを最小化する手法の一例として信頼度伝搬法を説明したが、本発明は、これに限定されない。本発明では、この手法として、例えば、ビタビアルゴリズムを用いることもできる。なお、各実施形態では、距離指標情報として距離情報を推定することとして説明したが、本発明は、これに限定されない。つまり、本発明は、前記した距離候補を視差候補に置き換えれば、視差情報を推定する視差情報推定装置を実現することができる。なお、各実施形態では、評価値として評価値K,Oを計算することとして説明したが、本発明は、これに限定されない。また、各実施形態では、評価値K,Oの両方を計算することとして説明したが、本発明は、評価値K,Oの何れか一方のみを計算してもよい。この場合、本発明は、テクスチャ又はオクルージョンに起因する推定エラーの何れか一方だけを抑制してもよい。変形例1以下、図18を参照して、本発明の変形例1に係る距離情報推定装置1Dについて説明する。図18に示すように、距離情報推定装置1Dは、図10の距離情報推定装置1Bに、図14のメッセージ伝搬制限判定部21を付加したものである。従って、距離情報推定装置1Dは、オクルージョン及びテクスチャに起因する推定エラーを抑制すると共に被写体の膨張を抑制し、精度が高い距離情報を生成することができる。なお、図18の各手段は、図10,図14と同様であるため、詳細な説明を省略する。以下、本発明の実施例として、推定エラーの判定精度を実験した結果について説明する。本実施例では、図3の距離情報推定装置1を用いて実験を行った。また、この距離情報推定装置1に入力するステレオ画像として、力士被写体の取組を撮影した画像を準備した。左カメラCLが撮影した基準画像を図19aに示し、右カメラCRが撮影した参照画像を図19bに示す。このとき、左カメラCL及び右カメラCRは、互いに2メートルだけ離した。また、左カメラCLから被写体までの距離が25メートルであり、距離候補を2センチメートルごとに設定した。さらに、λdataを0.07、誤差関数閾値Tdataを15、スムーズ関数閾値Tsmoothを1.7に設定した。図19のステレオ画像から、距離情報推定装置1が生成した距離画像を図20aに示す。図20aの距離画像には、オクルージョンに起因する推定エラーを生じた領域を赤色で図示した。また、テクスチャに起因する推定エラーを生じた領域を青色で図示した。図20aの距離画像から、距離情報推定装置1が、評価値Oを用いてオクルージョン領域を出力した結果を、図20bに示す。図20bの画像において、輝度値の低い領域がオクルージョンに起因する推定エラー領域が生じた領域オクルージョン領域であり、図20aの赤色領域とほとんど一致している。このことより、距離情報推定装置1は、オクルージョンに起因する推定エラーの判定精度が高いことが分かる。また、図20aの距離画像に対し、距離情報推定装置1が、オクルージョンに起因する推定エラーの抑制処理を行った。この推定エラーを抑制した後の距離画像を図21aに示す。ここで、図21aの距離画像にも、図20aと同様に赤色と、青色とを図示した。そして、図21aの距離画像から、距離情報推定装置1が、評価値Kを用いて無テクスチャ領域を出力した結果を、図21bに示す。図21bの画像において、輝度値の低い領域がテクスチャに起因する推定エラー領域が生じた領域例えば、土俵奥の地面であり、図21aの青色領域とほとんど一致している。このことより、距離情報推定装置1は、テクスチャに起因する推定エラーの判定精度が高いことが分かる。"