"複数の視点画像から視点画像間の視差を推定する立体画像処理装置であって、前記視点画像間での画素の類似度及び前記視点画像間の視差の連続性に関する評価値を算出する評価値算出部と、該評価値を最小とする各画素の視差の組合せを、各画素の視差と選択して視差マップを作成する視差マップ作成部とを備え、前記視差の連続性は、注目画素の視差と該注目画素の横及び縦に隣接する画素の視差との連続性であり、前記評価値は、前記評価値算出部において、前記画素の類似度が高いほど小さく算出され、前記視差の連続性が高いほど小さく算出される、ことを特徴とする立体画像処理装置。前記画素の類似度は、一方の視点画像における注目画素を含む所定範囲の周辺領域中の画素と、他方の視点画像における該注目画素に対応する画素を含む所定範囲の周辺領域中の画素との類似度であることを特徴とする請求項1に記載の立体画像処理装置。前記画素の類似度は、少なくとも画素の輝度値を用いて算出されることを特徴とする請求項1又は2に記載の立体画像処理装置。前記視差の連続性は、視差の横成分及び/又は縦成分を用いて算出されることを特徴とする請求項1~3のいずれか1項に記載の立体画像処理装置。前記視差マップ作成部は、前記評価値における前記視差の連続性に基づき、オクルージョン領域を検出するオクルージョン領域検出部を有することを特徴とする請求項1~4のいずれか1項に記載の立体画像処理装置。前記オクルージョン領域検出部は、前記評価値における前記視差の連続性に基づき、横方向のオクルージョン領域と縦方向のオクルージョン領域の両方を検出することを特徴とする請求項5に記載の立体画像処理装置。前記オクルージョン領域検出部は、前記評価値における視差の連続性に基づき、少なくとも横方向のオクルージョン領域を検出し、前記評価値算出部は、前記横方向のオクルージョン領域に隣接する右側と左側の画素間の視差の横成分の差分の大きさが、前記横方向のオクルージョン領域の横幅と等しくなるという条件下で、前記評価値の算出を行うことを特徴とする請求項5又は6に記載の立体画像処理装置。前記オクルージョン領域検出部は、前記評価値における視差の連続性に基づき、少なくとも縦方向のオクルージョン領域を検出し、前記評価値算出部は、前記縦方向のオクルージョン領域に隣接する上側と下側の画素間の視差の縦成分の差分の大きさが、前記縦方向のオクルージョン領域の縦幅と等しくなるという条件下で、前記評価値の算出を行うことを特徴とする請求項5又は6に記載の立体画像処理装置。前記評価値算出部は、前記視点画像間の視差の縦成分についての探索範囲を、注目画素の座標に応じて変化させることを特徴とする請求項1~8のいずれか1項に記載の立体画像処理装置。前記評価値算出部は、前記視点画像間の視差の横成分についての探索範囲を、注目画素の座標に応じて変化させることを特徴とする請求項1~8のいずれか1項に記載の立体画像処理装置。複数の視点画像から視点画像間の視差を推定する立体画像処理方法であって、前記視点画像間での画素の類似度及び前記視点画像間の視差の連続性に関する評価値を算出するステップと、前記評価値を最小とする各画素の視差の組合せを、各画素の視差と選択して視差マップを作成するステップとを有し、前記視差の連続性は、注目画素の視差と該注目画素の横及び縦に隣接する画素の視差との連続性であり、前記評価値は、前記評価値を算出するステップにおいて、前記画素の類似度が高いほど小さく算出され、前記視差の連続性が高いほど小さく算出される、ことを特徴とする立体画像処理方法。コンピュータに、複数の視点画像から視点画像間の視差を推定する立体画像処理を実行させるためのプログラムであって、前記立体画像処理は、前記視点画像間での画素の類似度及び前記視点画像間の視差の連続性に関する評価値を算出するステップと、該評価値を最小とする各画素の視差の組合せを、各画素の視差と選択して視差マップを作成するステップとを含み、前記視差の連続性は、注目画素の視差と該注目画素の横及び縦に隣接する画素の視差との連続性であり、前記評価値は、前記評価値を算出するステップにおいて、前記画素の類似度が高いほど小さく算出され、前記視差の連続性が高いほど小さく算出される、ことを特徴とするプログラム。しかしながら、特許文献1、非特許文献1に記載の技術をはじめとする従来の視差推定方法においては、視点画像間の対応点同士は、縦方向にずれが無いことが想定されており、縦方向にずれが生じている画像に対しては正しい視差が推定できないという問題があった。実際、縦方向のずれは、例えば交差法によって撮影された画像において生じる。交差法では光軸が平行でないため、エピポーラ線の傾きが左目用画像と右目用画像で異なり、領域によって程度の異なる縦方向のずれが生じる。すなわち、視差は横成分と縦成分を持つ。さらに、視差の縦成分が異なる領域が隣接した場合には、縦方向の遮蔽によるオクルージョン領域が生じる。具体例を挙げて説明する。図2A,図2Bは交差法撮影によって縦方向にずれが生じている画像の例であり、図2Aは左目用画像2L、図2Bは右目用画像2Rを表している。なお、図2A,図2Bにおける点線は、左目用画像2Lと右目用画像2Rとで同じ行の位置を示すために描いたものである。右目用画像2R内では左目用画像2Lよりも家が上方に位置し、縦方向にずれが生じていることが分かる。人物の位置も縦方向にずれているが、ずれの量は家よりも小さい。斜線で示した領域は、図1A,図1Bの例と同様、前景である人物に遮蔽されて生じたオクルージョン領域である。人物と家とでは縦方向のずれ量が異なっているため、右目用画像2Rにおいて、人物の頭より上の部分にもオクルージョン領域が存在している。すなわち、この部分は縦方向の遮蔽によって生じたオクルージョンである。図1A,図1Bの例では横方向にのみ遮蔽していたが、図2A,図2Bの例では、横方向と縦方向の両方で遮蔽が生じている。以降では、縦方向の遮蔽により生じるオクルージョンを縦方向オクルージョンと呼び、横方向の遮蔽によって生じる従来のオクルージョンを横方向オクルージョンと呼ぶことにする。そして、このような縦方向オクルージョン領域は、上述した従来の視差推定方法では考慮されておらず、視差推定の誤りの要因となり得る。縦方向のずれや縦方向オクルージョンは、左目用画像と右目用画像でレンズの倍率がずれている場合にも生じる。さらに、左目用画像と右目用画像の撮影位置が縦にずれている場合や、撮影方向が縦にずれている場合にも発生する。また、縦に配置されたカメラで撮影した視点画像間の視差を推定する場合においても、同様に横方向のずれや横方向オクルージョンが発生するという課題がある。さらにより一般的には、斜めに配置されたカメラで撮影した視点画像間の視差を推定する場合においても、同様の課題がある。本発明は、上述のような実情に鑑みてなされたものであり、その目的は、視差に縦成分がある画像の場合にも正しく視差を推定することが可能な立体画像処理装置、立体画像処理方法、及びプログラムを提供することにある。上記課題を解決するために、本発明の第1の技術手段は、複数の視点画像から視点画像間の視差を推定する立体画像処理装置であって、前記視点画像間での画素の類似度及び前記視点画像間の視差の連続性に関する評価値を算出する評価値算出部と、該評価値を最小とする各画素の視差の組合せを、各画素の視差と選択して視差マップを作成する視差マップ作成部とを備え、前記視差の連続性は、注目画素の視差と該注目画素の横及び縦に隣接する画素の視差との連続性であり、前記評価値は、前記評価値算出部において、前記画素の類似度が高いほど小さく算出され、前記視差の連続性が高いほど小さく算出される、ことを特徴としたものである。第2の技術手段は、第1の技術手段において、前記画素の類似度は、一方の視点画像における注目画素を含む所定範囲の周辺領域中の画素と、他方の視点画像における該注目画素に対応する画素を含む所定範囲の周辺領域中の画素との類似度であることを特徴としたものである。第3の技術手段は、第1又は第2の技術手段において、前記画素の類似度は、少なくとも画素の輝度値を用いて算出されることを特徴としたものである。第4の技術手段は、第1~第3のいずれか1の技術手段において、前記視差の連続性は、視差の横成分及び/又は縦成分を用いて算出されることを特徴としたものである。第5の技術手段は、第1~第4のいずれか1の技術手段において、前記視差マップ作成部は、前記評価値における前記視差の連続性に基づき、オクルージョン領域を検出するオクルージョン領域検出部を有することを特徴としたものである。第6の技術手段は、第5の技術手段において、前記オクルージョン領域検出部は、前記評価値における前記視差の連続性に基づき、横方向のオクルージョン領域と縦方向のオクルージョン領域の両方を検出することを特徴としたものである。第7の技術手段は、第5又は第6の技術手段において、前記オクルージョン領域検出部は、前記評価値における視差の連続性に基づき、少なくとも横方向のオクルージョン領域を検出し、前記評価値算出部は、前記横方向のオクルージョン領域に隣接する右側と左側の画素間の視差の横成分の差分の大きさが、前記横方向のオクルージョン領域の横幅と等しくなるという条件下で、前記評価値の算出を行うことを特徴としたものである。第8の技術手段は、第5又は第6の技術手段において、前記オクルージョン領域検出部は、前記評価値における視差の連続性に基づき、少なくとも縦方向のオクルージョン領域を検出し、前記評価値算出部は、前記縦方向のオクルージョン領域に隣接する上側と下側の画素間の視差の縦成分の差分の大きさが、前記縦方向のオクルージョン領域の縦幅と等しくなるという条件下で、前記評価値の算出を行うことを特徴としたものである。第9の技術手段は、第1~第8のいずれか1の技術手段において、前記評価値算出部は、前記視点画像間の視差の縦成分についての探索範囲を、注目画素の座標に応じて変化させることを特徴としたものである。第10の技術手段は、第1~第8のいずれか1の技術手段において、前記評価値算出部は、前記視点画像間の視差の横成分についての探索範囲を、注目画素の座標に応じて変化させることを特徴としたものである。第11の技術手段は、複数の視点画像から視点画像間の視差を推定する立体画像処理方法であって、前記視点画像間での画素の類似度及び前記視点画像間の視差の連続性に関する評価値を算出するステップと、前記評価値を最小とする各画素の視差の組合せを、各画素の視差と選択して視差マップを作成するステップとを有し、前記視差の連続性は、注目画素の視差と該注目画素の横及び縦に隣接する画素の視差との連続性であり、前記評価値は、前記評価値を算出するステップにおいて、前記画素の類似度が高いほど小さく算出され、前記視差の連続性が高いほど小さく算出される、ことを特徴としたものである。第12の技術手段は、コンピュータに、複数の視点画像から視点画像間の視差を推定する立体画像処理を実行させるためのプログラムであって、前記立体画像処理は、前記視点画像間での画素の類似度及び前記視点画像間の視差の連続性に関する評価値を算出するステップと、該評価値を最小とする各画素の視差の組合せを、各画素の視差と選択して視差マップを作成するステップとを含み、前記視差の連続性は、注目画素の視差と該注目画素の横及び縦に隣接する画素の視差との連続性であり、前記評価値は、前記評価値を算出するステップにおいて、前記画素の類似度が高いほど小さく算出され、前記視差の連続性が高いほど小さく算出される、ことを特徴としたものである。本発明によれば、視差に縦成分がある画像の場合にも誤りの少ない視差推定を行うことができる。以下、添付図面を参照しながら本発明の実施の形態について詳細に説明する。図面において同じ機能を有する部分については同じ符号を付し、繰り返しの説明は省略する。第一の実施形態図3は、本発明の第一の実施形態に係る立体画像処理装置の概略構成例を示すブロック図で、図4は、本発明の第一の実施形態に係る評価値の求め方を説明するための図である。図3で示すように、本実施形態の立体画像処理装置10は、複数の視点画像を入力する入力部11と、入力された複数の視点画像について視点画像間での画素値の類似度及び視点画素間の視差の連続性に関する評価値を算出する評価値算出部12と、算出された評価値を用いて視差マップを作成する視差マップ作成部13とを備えている。視差マップ作成部13は、算出された評価値を最小とする各画素の視差の組合せを、各画素の視差と選択して視差マップを作成する。無論、入力部11の代わりに複数の視点画像を記憶した記憶部を備え、評価値算出部12がそこから複数の視点画像を読み出すような構成を採用してもよい。そして、本発明における評価値算出部12では、視差の連続性として、注目画素の視差とその注目画素の横及び縦に隣接する画素の視差との連続性を用いる。さらに、本発明における評価値算出部12では、評価値は、画素の類似度が高いほど小さく算出され、視差の連続性が高いほど小さく算出されるものとする。以下の説明では、理解を容易にするために、まず、評価値算出部12で注目画素の視差とその注目画素の横に隣接する画素の視差との連続性を算出し、その結果の連続性と算出した画素の類似度とを用いて評価値を算出する例を挙げて説明する。注目画素の縦に隣接する画素の視差との連続性の算出、並びにその結果及び算出した画素の類似性を用いた評価値の算出については後述する。本実施形態の立体画像処理装置10は、入力部11で入力された視点画像間について、1行ずつ独立に視差推定を行う。この視差推定は、評価値算出部12及び視差マップ作成部13が実行し、各行の処理は、動的計画法を用いて行う。つまり、評価値算出部12及び視差マップ作成部13は、視差の横成分と縦成分とを、動的計画法を用いて検出する。評価値算出部12では、行内の各画素の視差値の組合せを経路として扱い、各組合せに対する評価値を算出する。但し、動的計画法では全ての組合せについての評価値を求めずに、評価値を最小にする可能性のある組合せだけを選んで経路を探索するため、計算量が削減される。視差マップ作成部13は、オクルージョン領域検出部101、視差横成分検出部102、視差縦成分検出部103を有し、評価値算出部12で算出された評価値を最小とする視差値の組合せを選択しオクルージョンマップ、視差横成分マップ、視差縦成分マップを生成し、出力する。オクルージョンマップ、視差横成分マップ、視差縦成分マップを併せて、視差マップと呼ぶことにする。ここでは、好ましい例として、視差マップ作成部13がオクルージョン領域検出部101を有する例、つまり視差マップ作成部13がオクルージョン領域の検出を行う例を挙げるが、視差マップ作成部13は少なくとも視差の横成分と縦成分とを検出して、視差横成分マップ及び視差縦成分マップでなる視差マップを作成できる構成であればよい。このような構成において、立体画像処理装置10では、まず入力部11で入力された視点画像を予め輝度画像に変換する。ここでは、入力部11を通して入力される複数の視点画像が左目用画像及び右目用画像の合計2枚から構成される例を挙げる。例えば、入力部11が左目用画像と右目用画像とを入力し、輝度画像への変換も施した上で、評価値算出部12に出力する。無論、入力部11が予め輝度画像に変換された左目用画像及び右目用画像を入力してもよい。以降では、左目用画像を基準とした右目用画像の視差マップを推定する場合について説明する。また、ここでは、視差を求めるのに最適な輝度画像つまり画像の輝度値を用いて処理を行うことを前提としているが、輝度値に限らず、左目用画像及び右目用画像の画素値、例えばRGBデータであればそのいずれか1又は複数色の値を用いて処理を行うこともできる。評価値算出部12では、まず、左目用画像の輝度値に対する右目用画像の輝度値の類似度Mx,y,p,qを、次の1式を用いて、0≦x≦w-1,0≦p≦P,-Q≦q≦Qの範囲で計算する。ここで、L、Rはそれぞれ左目用画像と右目用画像の輝度値、x,yは画像中の座標を示し、yは推定を行っている注目行を示す。p,qはそれぞれ視差の横成分,縦成分を表す。P,Qはそれぞれ視差の横成分、縦成分の探索範囲を示し、p,qはそれぞれ0≦p≦P、-Q≦q≦Qの範囲の値のみをとることとする。wは画像の横画素数、2a1は-a≦i≦a及び-a≦j≦aで定義されるブロックサイズである。aの値は任意に決めておけばよいし、探索範囲は矩形に限ったものでもない。なお、1式の計算の際、座標が画像の外の位置を示した場合は、その座標を近傍の画像端の座標で置き換えるなどの画像端処理を行えばよい。1式の類似度Mx,y,p,qは、左目用画像に対する右目用画像の輝度値の差分右目用画像の画素x,yの輝度値に対して、p,qだけずらした位置の左目用画像の輝度値の差分の絶対値和を用いて定義した値であり、小さい程、比較対象の画素値の類似度が高いと言える。但し、この例に限らず、類似度は、差分平方和、正規化相互相関などの類似度を表す他の関数で定義してもよい。また、探索範囲をP,Qで定めた例を挙げて本実施形態を説明しているように、評価値算出部12は、注目画素の所定範囲上記ブロックサイズの範囲の周辺領域に対して、視点画像間での輝度値の類似度を用いて評価値を算出することが望ましい。つまり、類似度として、一方の視点画像における注目画素を含む所定範囲の周辺領域中の画素と、他方の視点画像におけるその注目画素に対応する画素を含む所定範囲の周辺領域中の画素との類似度を用いることが好ましい。但し、探索範囲を特に決めずに画像の全領域について類似度を求め、評価値を算出することもできる。次に、下の2式に従い、評価値Ux,pを0≦x≦w-1,0≦p≦Pを満たす全てのx,pについて算出する。評価値Ux,pは、上述の類似度Mx,y,p,qの最小値に基づいて計算した値である。類似度Mx,y,p,qの値が小さい程、比較対象の画素値が類似しているためである。計算の手順としては、2式を用い、まずx=0の場合を計算し、その後x=1からw-1まで順に計算する。2式は再帰的になっており、x=0から順に累積して最後に得られたUw-1,pが注目行内の視差値の組合せに対する評価値となる。0≦p≦Pであるため、Uw-1,pはP1個得られ、P1通りの経路に対する評価値である。4式のAxは、右目用画像の注目画素x座標がxの視差横成分がpであると仮定した場合の、その左隣の視差横成分の候補を示している。右目用画像の注目画素xは左目用画像の画素xpと対応している。Axの要素を3つに分類して説明する。[I]要素x-1,p左隣の画素x-1と注目画素の視差横成分が同じ場合である。右目用画像の画素x-1と左目用画像の画素xp-1が対応する。[II]要素x-1,p-1,...,x-1,0左隣の画素の視差横成分が注目画素より小さい場合である。例えば、要素x-1,p-2は、右目用画像の画素x-1と左目用画像の画素xp-3が対応することを示す。このとき、左目用画像のxp-1とxp-2の2つの画素は、右目用画像内に対応点が無い。つまり、左目用画像に横方向オクルージョン領域がある場合である。[III]要素x-2,p1,...,x-P-1p,P注目画素の左隣の画素が横方向オクルージョン領域の場合である。この場合は、左目用画像において注目画素と対応した画素xpの左隣の画素xp-1について、対応画素を右目用画像から選ぶ。例えば、要素x-3,p2は、右目用画像の画素x-3と左目用画像の画素xp-1が対応することを示す。このとき、右目用画像のx-1とx-2の2つの画素は、左目用画像内に対応点が無く、横方向オクルージョン領域となる。また、右目用画像の画素x-3と注目画素との視差横成分の差分が2であり、横方向オクルージョン領域の幅と等しい。つまり、横方向オクルージョン領域の横幅がその右側と左側の画素間の視差横成分の差分と等しくなるように制約されている。図4は4式で示す集合Axの範囲を示しており、横軸はx座標、縦軸は視差横成分である。格子模様を付けた座標は注目画素xを示し、その視差横成分がpである。このとき、黒で示したP1個の座標がAxの要素に相当する。また、図4中、斜線で示した座標は、画像端における視差の探索範囲外の領域を示しており、次式の集合Bで表わされる。2式の計算において、この領域では、Ux,pの値を、無限大とみなせる大きな値で置き換える。また、2式の計算において、minqで最小値を与えたqの値、つまりMx,y,p,qが最小となったときのqの値を、Kqx,pとして保持しておく。同様に、minxb,pbで最小値を与えたxbの値、pbの値、つまりUxb,pbDpxp,pbが最小となったときのxbの値、pbの値を、それぞれKxx,p、Kpx,pとして保持しておく。このKxとKpは、評価値Uに対する経路を表している。また、xb、pbは、それぞれ右目用画像の集合Axにおける横方向の座標、視差の横成分の値である。なお、2~4式の技術的な意味については後述する。視差マップ作成部13では、算出した評価値Uw-1,pとKxx,p、Kpx,p、Kqx,pとを用いて、注目行の各画素の視差マップを作成する。まず、動的計画法に基づき、評価値算出部12で得られたP1個の評価値Uw-1,pの中から、最小の評価値を選択する。但し、本例では式5の領域においてUx,pの値を無限大としたため、自動的にUw-1,0が最小となる。そして、その評価値における経路をKxとKpを参照して右端画素x=w-1から順に左へ向かってたどり、各画素の視差横成分と視差縦成分を検出する。併せてそれらの視差横成分からオクルージョン領域の検出も行う。視差マップ作成部13における処理の一例について、その流れを図5を参照しながら説明する。例として、Uw-1,0が最小となった場合について示したものである。まず、ステップS1において、変数x,xb,pbの値を初期化する。xは注目行における右目用画像の横方向の座標であり、この処理例では注目行の右端画素から順に処理を行うために、初期値としてw-1が入力されている。xbは右目用画像の集合Axにおける横方向の座標であり、xと同様の理由からw-1が初期値として入力されている。pbは右目用画像の集合Axにおける視差の横成分の値であり、初期値として0が入力されている。次に、ステップS2において、xbとxの値を比較する。値が異なる場合、すなわち、評価値算出部12においてAxの中から要素[III]が選択されていた場合、ステップS3に進む。ステップS3では、座標x,yが横方向オクルージョン領域であると判定し、オクルージョンマップにオクルージョン領域に該当することを示すフラグを出力する。その後、ステップS8に進む。ステップS2において値が同じである場合、すなわち、評価値算出部12においてAxの中から要素[I]又は[II]が選択されていた場合、ステップS4に進み、pにpbの値を入力する。続いて、ステップS5において、座標x,yの視差横成分の値をpと検出し、視差横成分マップに出力する。さらに、ステップS6において、座標x,yの視差縦成分の値をKqx,pと検出し、視差縦成分マップに出力する。ステップS5,S6の処理後は、ステップS7において、xbの値、pbの値を、それぞれKxx,p、Kpx,pで更新する。また、ステップS3の処理後、並びにステップS7の処理後は、ステップS8において、x=0となったか否かを判定する。そして、x=0である場合は左端画素まで処理が終了しているため、図5の処理を終了し、それ以外の場合は、ステップS9において左隣りの画素に移動し、ステップS2に戻る。上述した処理において、ステップS2,S3はオクルージョン領域の判定を行う処理であり、これらの処理はオクルージョン領域検出部101が実行する。ステップS4,S5は視差横成分を検出する処理であり、これらの処理は視差横成分検出部102が実行する。ステップS6は視差縦成分を検出する処理であり、この処理は視差縦成分検出部103が実行する。このようにして、注目行内の全画素の視差全てのxについての視差が推定され、同時に、行内のオクルージョン領域が検出される。そして、評価値算出部12と視差マップ作成部13が全ての行に対して処理を実行することで、画像全体の視差が推定され、視差マップが作成される。次に、上述した2~4式について、その技術的意味を説明する。2式のminqの項は、視差の縦成分を探索する効果を持つ。各視差横成分pに対して、視差縦成分qの値を探索範囲内で変化させた際に最小となる類似度を用いて評価値に加算する。また、2式のminxb,pbの項は、左隣りの画素における最小評価値を用いて評価値を累積する動的計画法の手法が適用されている。但し、非特許文献1と同様に、4式で定めたAxの範囲を用いることによって、オクルージョン領域の類似度を除いて評価値を累積している。また、4式で示す集合Axは、オクルージョン領域の左端と右端に隣接する画素間の視差横成分の差分の大きさを、オクルージョン領域の横幅と等しくなるように制約するためのものである。つまり、ここでは、評価値算出部12は、横方向オクルージョン領域に隣接する右側と左側の画素間の視差横成分の差分の大きさが横方向オクルージョン領域の横幅と等しくなるという条件下で、評価値の算出を行っている。よって、オクルージョン領域検出部101では、横方向のオクルージョン領域を検出していると言える。このような範囲の集合Axを採用することで、オクルージョン領域検出と視差推定の誤りを生じにくくする効果がある。また、2式の定数項c0は、対応点の個数に応じて評価値に報酬を与え、オクルージョン領域を誤って広く検出することを防ぐ効果がある。また、3式のDpxは、視差の横成分の連続性を示す値であり、視差の横成分が不連続な場合にペナルティ定数項c1を与えてその視差を検出されにくくするものであり、平滑化の効果を持つが、この式に限らず、2つの視差の差分を用いた線形関数などの他の関数によって平滑化の効果を与えてもよい。無論、平滑化の効果を鑑みなければ、2式におけるDpxは常に0、つまりDpxの項を除いた式を採用することもできる。このように、2式の評価値Uw-1,pは、オクルージョン領域と推定した画素を除いた全ての画素における類似度の縦成分探索範囲内での最小値と、視差の横成分の横方向の連続性を示す値との加算値の合計、に基づき算出される。評価値算出部12では、このようにして視点画像間での画素の類似度、並びに視点画像間の視差の横方向の連続性つまり、注目画素の視差とその注目画素の横に隣接する画素の視差との連続性を用いて、評価値を算出する。よって、算出される評価値はこれら類似度及び連続性に関する値と言える。そして、評価値算出部12で算出される評価値は、視点画素間の画素の類似度が高いほど上記Mx,y,p,qが小さいほど小さく算出され、視差の横方向の連続性が高いほど上記Dpxが小さいほど小さく算出されているのが分かる。また、本実施形態では、視差の横成分の範囲を0≦p≦Pとすることでpが正値のみを取る場合について説明したが、pが負値を取り得る場合についても、同様の手法を適用することができる。また、右目用画像の視差マップを求める場合について説明したが、逆に右目用画像に対する左目用画像の視差マップについても同様に求めることができる。上述のとおり、本実施形態の立体画像処理装置によれば、視差推定を行う際に、視点画像間の画素の類似度によって視差の横成分と縦成分の両方を推定し、オクルージョン領域の検出を同時に行うことにより、視差に縦成分のある画像に対しても誤りの少ない視差推定を行うことができる。この効果は、視差の連続性として、注目画素の視差とその注目画素の横に隣接する画素の視差との連続性を用いるだけでも得られるものである。なお、本実施形態では、視差マップ作成部13が、オクルージョン領域の検出と視差の横成分と縦成分の検出処理を行うことを前提に説明したが、オクルージョン領域の検出を行わない場合にも、視差の両方向成分を推定しているため、視差に縦成分のある画像に対しても誤り少なく視差を推定することができる。次に、第一の実施形態の主たる特徴として、上述した横方向の連続性を採用した例を参照しながら、注目画素の視差とその注目画素の横及び縦に隣接する画素の視差との連続性を算出し、その結果及び算出した画素の類似性を用いて評価値を算出する方法について説明する。ここでは、評価値は、画素の類似度が高いほど小さく算出され、視差の横方向及び縦方向の連続性が高いほど小さく算出されることになる。第一の実施形態は、視差の横方向と縦方向の両方に平滑化の効果を与えるように上述した横方向のみの連続性を採用した例を変更したものである。よって、本実施形態における立体画像処理装置の概略構成例は、上述したように図1の立体画像処理装置10で示されるが、評価値算出部12における処理の方法が上述した横方向のみの連続性を採用した例とは異なる。その相違点について説明すると、上述した横方向のみの連続性を採用した例における評価値算出部12は、注目画素の視差の横成分について横方向の連続性を求め、評価値の算出に用いている。これに対し、第一の実施形態における評価値算出部12は、注目画素の視差の横成分と縦成分について、それぞれ横方向と縦方向の連続性を求め、評価値の算出に用いている。より具体的な例を挙げると、第一の実施形態では、上端の行から順に1行ずつ視差推定を行うが、その際、前の行で求めた視差の値を用いる。前行で求めた視差の横成分、縦成分をそれぞれFpx、Fqxとして説明する。評価値算出部12では、まず、上述した横方向のみの連続性を採用した例と同様に類似度Mx,y,p,qを計算する。評価値算出部12では、次に、下の6式に従い、評価値Ux,pを0≦x≦w-1,0≦p≦Pを満たす全てのx,pについて算出する。計算の手順としては、6式を用い、まずx=0の場合を計算し、その後x=1からw-1まで順に計算する。ここで、Dqx,Dpy,Dqyはそれぞれ次式で示す関数であり、c2,c3,c4,c5,c6は定数である。occは、その画素がオクルージョン領域であることを示す定数である。Dpx、Axは、それぞれ3式、4式をそのまま用い、残りの処理は、上述した横方向のみの連続性を採用した例と同様に行う。6式の計算において、minqで最小値を与えたqの値をKqx,pとして保持しておく。同様に、minxb,pbで最小値を与えたxbの値、pbの値をそれぞれKxx,p、Kpx,pとして保持しておく。なお、6~9式の技術的な意味については後述する。但し、注目行が画像の上端の行の場合は、Fpx、Fqxに値が入っていないため、6式の代わりに下の10式を用いる。視差マップ作成部13では、上述した横方向のみの連続性を採用した例と同様の処理を行う。注目行内の全ての画素の視差推定が終了した時点で、行内の全ての視差の横成分と縦成分をFpxとFqxとして保持する。但し、オクルージョン領域と検出された画素については、FpxとFqxの両方にoccを代入する。その後、次の行に移る。以上の処理を全ての行で行うことで、画像全体の視差が推定され、視差マップが作成される。また、本実施形態では、最上行から順に処理を行った例を挙げたが、これに限らず、例えば最下行から順に上に向かって処理を行ってもよい。他の例として、中心行の推定を最初に行い、中心から順に上に向かって最上行まで処理を行い、その後、中心から順に下に向かって最下行まで処理を行ってもよい。これは、前行の誤りが注目行に伝播する場合があることを考慮し、誤りが生じにくい行から推定処理を行うことによって、画像全体での誤りを少なく抑える効果がある。ここで、上述した6~9式について、その技術的意味を説明する。7式のDqxは視差縦成分の横方向の連続性を示す値であり、8式のDpyと9式のDqyはそれぞれ視差の横成分と縦成分の縦方向の連続性を示す値である。Dqx、Dpy、Dqyはそれぞれ不連続な場合にペナルティを与えている。これにより、2次元的な平滑化の効果を与えることができ、特に、画像中の模様が少ない領域、若しくは、繰返し模様のように、局所的な類似度だけでは視差を一意に決められない領域に対して、視差推定の誤りを減らすことができる。そして、6式の評価値Uw-1,pは、これらの連続性を示す値と1式の類似度Mx,y,p,qとに基づいて算出された値である。より具体的には、6式で示す評価値Uw-1,pは、オクルージョン領域と推定した画素を除いた全ての画素における類似度の縦成分探索範囲内での最小値と、視差の横成分と縦成分の横方向と縦方向の連続性を示す値との加算値の合計、に基づき算出される。評価値算出部12では、このようにして視点画像間での画素の類似度、並びに視点画像間の視差の横方向及び縦方向の連続性つまり、注目画素の視差とその注目画素の横及び縦に隣接する画素の視差との連続性を用いて、評価値を算出する。よって、算出される評価値はこれら類似度及び連続性に関する値と言える。そして、評価値算出部12で算出される評価値は、視点画素間の画素の類似度が高いほど上記Mx,y,p,qが小さいほど小さく算出され、視差の横方向及び縦方向の連続性が高いほど上記Dpx及び上記Dqyが小さいほど小さく算出されているのが分かる。なお、10式については、基本的に6式と同様の意味をもっており、上述したように注目行が上端の行であってFpx,Fqxの値が入っていない場合について特化したものである。上述のとおり、本実施形態の立体画像処理装置によれば、視差推定を行う際に、視差に対して横方向と縦方向に平滑化を行うことによって、視差に縦成分があり、かつ模様が曖昧な画像に対しても誤りの少ない視差推定を行うことができる。また、注目画素の視差の縦成分についての連続性を評価値の算出に用いればよく、視差の横成分の連続性を用いなくてもよいが、本実施形態では、視差の性質上横成分の連続性の方がより重要視されることから、横成分の連続性も用いている。さらに、本実施形態では、視差の縦成分の連続性を用いるに際し、縦成分の横方向の連続性も併せて用い、視差の横成分の連続性を用いるに際し、横成分の縦方向の連続性も併せて用いているが、より好ましい例を挙げたに過ぎず、少なくとも縦成分の縦方向の連続性と横成分の横方向の連続性を用いていればよい。なお、上述した横方向のみの連続性を採用した例では、視差の横成分について横方向の連続性を評価値の算出に用いた例を挙げたが、ここでの説明を適用するなどして、視差横成分の縦方向の連続性も用いるように構成することもできる。第二の実施形態第二の実施形態は、第一の実施形態における視差縦成分の探索範囲を修正したものである。一般的に、レンズ歪は画像中心から遠い程大きくなり、また、交差法撮影における光軸の角度によって生じる視差縦成分は、被写体の距離が一定の場合には画像中心から遠い程、大きくなる。このため、例えば、撮影したカメラの光軸やレンズ等の特性が測定によって予め分かっている場合、視差縦成分の分布傾向が予測できる。本実施形態では、評価値算出部12において、このように視差縦成分の画像内での分布傾向が予め分かっている場合、それに応じて視点画像間の視差縦成分の探索範囲つまり上記所定範囲の縦方向の範囲を座標に応じて変動させる。つまり、第二の実施形態における評価値算出部12は、視点画像間の視差の縦成分についての探索範囲を注目画素の座標に応じて変化させる。例えば、画像の中心付近では視差縦成分の大きさが小さいことが分かっている画像の場合、2式においてminを求める範囲-Q≦q≦Qの代わりに、次の11式を用いる。但しw,hはそれぞれ画像の横画素数、縦画素数である。上述のとおり、本実施形態の立体画像処理装置によれば、視差縦成分の画像内での分布傾向が予め分かっている場合、それに応じて視差縦成分の探索範囲を座標に応じて変動させることで、誤りの少ない視差推定を行うことができる。また、第二の実施形態では、より好ましい例として、第一の実施形態における視差横成分の探索範囲も修正することができる。すなわち、同様に、視差横成分の画像内での分布傾向が予め分かっている場合、それに応じて視点画像間の視差縦成分の探索範囲つまり上記所定範囲の横方向の範囲を座標に応じて変動させることで、誤りの少ない視差推定を行うことができる。なお、視差横成分についても、一般的に、レンズ歪は画像中心から遠い程大きくなり、また、交差法撮影における光軸の角度によって生じる視差横成分は、被写体の距離が一定の場合には画像中心から遠い程、大きくなるため、分布傾向の予測は可能である。このような視差横成分の探索範囲を変動させる処理と、視差縦成分の探索範囲を変動させる処理とは、いずれか一方のみを採用することもできる。第三の実施形態第三の実施形態は、横方向オクルージョン領域と縦方向オクルージョン領域の両方を検出できるように第一の実施形態を変更したものである。無論、ここで説明する第三の実施形態は、第二の実施形態と併せて採用することができる。図6は、本発明の第三の実施形態に係る立体画像処理装置の概略構成例を示すブロック図である。図6で示すように、本実施形態の立体画像処理装置20は、入力部21と、評価値算出部22と、視差マップ作成部23とを有している。視差マップ作成部23は、横方向オクルージョン領域検出部201、視差横成分検出部202、縦方向オクルージョン領域検出部203、視差縦成分検出部204を有し、それぞれ横方向オクルージョンマップ、視差横成分マップ、縦方向オクルージョンマップ、視差縦成分マップを生成、出力する。横方向オクルージョンマップ、縦方向オクルージョンマップ、視差横成分マップ、視差縦成分マップを併せて、視差マップと呼ぶことにする。本実施形態では、上端の行から順に1行ずつ視差推定を行うが、前行で求めた視差の縦成分Fqxと前行で求めたオクルージョンフラグFcxを用いる。入力部21は、第一の実施形態の入力部11と同様である。評価値算出部22では、まず、第一の実施形態の1式と同様に、注目行yについて類似度Mx,y,p,qを計算する。次に、下の12式に従い、評価値Ux,pを0≦x≦w-1,0≦p≦Pを満たす全てのx,pについて算出する。計算の手順としては、12式を用い、まず、x=0の場合を計算し、その後x=1から順にw-1まで計算する。Axは、第一の実施形態と同様に式4を満たすx,pの集合である。式12式は再帰的になっており、最後に得られたUw-1,pが注目行内の視差値の組合せに対する評価値となる。0≦p≦Pであるため、Uw-1,pはP1個得られ、P1通りの経路に対する評価値である。14式のAyは、右目用画像の注目行yにおける注目画素座標x,yの視差横成分がpであると仮定した場合の視差縦成分の候補を示している。Fcx=0のとき、注目画素の上隣の画素x,y-1は横方向オクルージョン領域でも縦方向オクルージョン領域でもなく、その視差縦成分はFqxである。また、その視差横成分をFpと示すと、右目用画像の画素x,y-1と左目用画像の画素xFp,y-1Fqxが対応している。この場合のAyの要素を3つに分類して説明する。[i]要素y,Fqx上隣の画素x,y-1と注目画素の視差縦成分が同じ場合である。右目用画像の画素x,yと左目用画像の画素xp,yFqxが対応する。[ii]要素y,Fqx1,...,y,Q注目画素の視差縦成分が、上隣の画素より大きい場合である。例えば、要素y,Fqx2は、右目用画像の画素x,yと左目用画像の画素xp,yFqx2が対応することを示す。このとき、左目用画像のxp,yFqxとxp,yFqx1の2つの画素は右目用画像内に対応点が無いと考える。つまり、左目用画像に縦方向オクルージョンがある場合である。[iii]要素y1,Fqx-1,...,yFqxQ,-Q注目画素が縦方向オクルージョン領域の場合である。この場合は、左目用画像の画素xp′,yFqxの対応画素を右目用画像から選ぶ。例えば、要素y2,Fqx-2は、右目用画像の画素x,y2と左目用画像の画素xp′,yFqxが対応することを示す。このとき、右目用画像のx,yとx,y1の2つの画素は、左目用画像内に対応点が無く、縦方向オクルージョン領域となる。また、右目用画像のx,y-1とx,y2との視差縦成分の差分が2であり、縦方向オクルージョン領域の縦幅と等しい。つまり、縦方向オクルージョン領域の縦幅がその上側と下側の画素間の視差縦成分の差分と等しくなるように制約されている。Fcx≠0のとき、注目画素の上隣の画素は横方向オクルージョン領域又は縦方向オクルージョン領域である。Fcx1は、注目画素からオクルージョン領域の上端に隣接する画素までの距離であり、その画素の視差縦成分がFqxに保持されている。その画素の視差横成分をFpと示すと、右目用画像のx,y-1-Fcxと左目用画像のxFp,y-1-FcxFqxが対応している。この場合のAyの要素を2つに分類して説明する。[iv]要素y,Fqx-Fcx注目画素が縦方向オクルージョン領域でない場合である。右目用画像のx,yと左目用画像のxp,yFqx-Fcxが対応する。注目画素の視差縦成分はFqx-Fcxであり、縦方向オクルージョン領域の上端に隣接する画素との視差縦成分の差分がFcxとなる。これは、縦方向オクルージョン領域の縦幅と等しい。[v]要素y1,Fqx-Fcx-1,...,yFqx-FcxQ,-Q注目画素が縦方向オクルージョン領域の場合である。この場合は、左目用画像xp′,yFqx-Fcxの対応画素を右目用画像から選ぶ。12式の計算において、minyb,qbで最小値を与えたybとqbの値、つまりM′x,yb,p,qbが最小となったときのybの値、qbの値を、それぞれKyx,p,Kqx,pに保持しておく。また、第一の実施形態と同様に、minxb,pbで最小値を与えたxb,pbの値をそれぞれKxx,p、Kpx,pに保持しておく。但し、注目行が画像の上端の行の場合は、Fqx、Fcxに値が入っていないため、12式の代わりに2式を用いる。これらのKx、Kp、Ky、Kqは、評価値Uに対する経路を表している。なお、12~14式の技術的な意味については後述する。視差マップ作成部23では、算出された評価値から視差の横成分及び縦成分を検出するが、算出した評価値Uw-1,pとKxx,p、Kpx,p、Kqx,pとを用いて、注目行の各画素について、視差の横成分縦成分の検出と共にオクルージョン領域か否かの判定も行う。オクルージョン領域は、縦方向オクルージョンと横方向オクルージョンの両方を判定する。視差マップ作成部23における処理の一例について、その流れを図7を参照しながら説明する。例として、Uw-1,0が最小となった場合について示したものである。まず、ステップS11において、変数x,xb,pbの値を初期化する。この初期化については図5のステップS1で説明したとおりである。次に、ステップS12において、xbとxの値を比較する。値が異なる場合、ステップS13に進む。ステップS13では、座標x,yが横方向オクルージョンであると判定し、横方向オクルージョンマップに横方向オクルージョン領域に該当することを示すフラグを出力する。その後、ステップS21でFqx、Fcxの値をそれぞれ-Q-1、0に更新し、ステップS23に進む。ステップS12において値が同じである場合、ステップS14に進み、pにpbの値を入力する。さらに、ステップS15において、座標x,yの視差横成分の値をpと検出し、視差横成分マップに出力する。次に、ステップS16において、座標x,yの視差縦成分Kqx,pと、前行の視差縦成分から前行の縦オクルージョン領域幅を引いた差分、すなわちFqx-Fcxの値を比較する。Kqx,pがFqx-Fcx以上の場合、すなわち評価値算出部22においてAyの中から要素[i]、[ii]、[iv]のいずれかが選択されていた場合、ステップS18において、座標x,yの視差縦成分の値をKqx,pと検出し、視差縦成分マップに出力する。ステップS18の後、ステップS19においてFqx、Fcxの値をそれぞれKqx,p、0で更新し、ステップS22においてxb、pbの値をそれぞれKxx,p、Kpx,pで更新する。ステップS16においてKqx,pがFqx-Fcx未満の場合、すなわち評価値算出部22においてAyの中から要素[iii]又は[v]が選択されていた場合、ステップS17に進む。ステップS17では、座標x,yが縦方向オクルージョンであると判定し、縦方向オクルージョンマップに縦方向オクルージョン領域に該当することを示すフラグを出力する。ステップS17の後、ステップS20においてFcの値を1だけ増やすように更新し、その後、ステップS22においてxb、pbの値をそれぞれKxx,p、Kpx,pで更新する。ここでFqxの値は更新せずそのままとする。最後に、ステップS21の処理後、並びにステップS22の処理後は、ステップS23においてx=0となったか否かを判定する。x=0である場合は左端画素まで処理が終了しているため、図7の処理を終了し、それ以外の場合は、ステップS24において左隣りの画素に移動し、ステップS12に戻る。但し、注目行が画像の上端の行の場合は、予めFqx=-Q-1,Fcx=0で初期化しておく。なお、ステップS21における-Q-1とは、Kqx,pがとり得る値の範囲、すなわち視差の縦成分探索範囲よりも小さい値であり、次にステップS16を通る際に、Kqx,pがFqx-Fcx以上であると判定させるために用いる値である。上述した処理において、ステップS12,S13は、横方向オクルージョン領域の判定を行う処理であり、これをまとめて横方向オクルージョン領域検出部201が実行する。ステップS14,S15は、視差横成分を検出する処理であり、これをまとめて視差横成分検出部202が実行する。ステップS16,S17は、縦方向オクルージョン領域の判定を行う処理であり、これをまとめて縦方向オクルージョン領域検出部203が実行する。ステップS18は、視差縦成分を検出する処理であり、これを視差縦成分検出部204が実行する。このようにして、注目行内の全画素の視差が推定され、同時に、行内の横方向オクルージョン領域と縦方向オクルージョン領域が検出される。評価値算出部22と視差マップ作成部23を全ての行で行うことで、画像全体の視差が推定され、視差マップが作成される。次に、上述した12~14式について、その技術的意味を説明する。12式のminyb,qbの項は、視差の縦成分を探索する効果を持つが、14式で定めたAyの範囲と13式を用いることによって、縦方向のオクルージョン領域の類似度を用いず、代わりにオクルージョン領域の下端に隣接する画素の類似度を用いて評価値を累積している。13式は、yb=y、すなわち縦方向オクルージョン領域でない場合には、類似度としてそのままのMx,yb,p,qbを用い、yb≠y、すなわち縦方向オクルージョン領域の場合には、縦方向オクルージョン領域の下端に隣接する画素の類似度を用いている。この際、オクルージョン領域の下端に隣接する画素の視差横成分はまだ検出されていないため、minp′により視差横成分探索範囲内で最小となる類似度を用いる。14式で示す集合Ayは、縦方向オクルージョン領域の上端と下端に隣接する画素間の視差縦成分の差分の大きさが、縦方向オクルージョン領域の縦幅と等しくなるように制約するためのものである。つまり、ここでは、評価値算出部22は、縦方向のオクルージョン領域に隣接する上側と下側の画素間の視差縦成分の差分の大きさが、縦方向オクルージョン領域の縦幅と等しくなるという条件下で、評価値の算出を行う。これは、オクルージョン領域検出と視差推定の誤りを生じにくくする効果がある。そして、12式の評価値Uw-1,pは、横方向オクルージョン領域と推定した画素を除いた全ての画素における類似度の縦成分探索範囲内での最小値と、視差の横成分の横方向の連続性を示す値との加算値の合計、に基づき算出される。評価値算出部22では、このようにして視点画像間での画素の類似度、並びに視点画像間の視差の横方向の連続性を用いて、評価値を算出する。そして、評価値算出部22で算出される評価値は、視点画素間の画素の類似度が高いほど上記M′x,yb,p,qbが小さいほど小さく算出され、視差の横方向の連続性が高いほど上記Dpxが小さいほど小さく算出されているのが分かる。そして、本実施形態においても、12式を6式と同様にして視点画像間の視差の縦方向の連続性も用いて評価値を算出するように改良すればよい。これにより、評価値算出部22で算出される評価値は、視点画素間の画素の類似度が高いほど上記M′x,yb,p,qbが小さいほど小さく算出され、視差の横方向及び縦方向の連続性が高いほど上記Dpx、Dpy、Dqx及び上記Dqyが小さいほど小さく算出されるようになる。上述のとおり、本実施形態の立体画像処理装置によれば、視差推定を行う際に、縦方向オクルージョン領域の検出を同時に行うことにより、視差に縦成分のある画像に対してもさらに誤りの少ない視差推定を行うことができる。第一~第三の実施形態について以上、本発明の各実施形態においては上記の通り、入力部11を通して入力される複数の視点画像が左目用画像及び右目用画像の2枚の視点画像、すなわち横に配置されたカメラによる視点画像から構成される例を挙げた。但し、視点画像の数は2枚に限らず、3枚以上であっても、各実施形態において説明した処理により同様に視点画像間の視差は推定できる。例えば、ある基準となる視点画像に対する他の視点画像の視差を推定することもできるし、各視点画像間の全ての組合せについて視差を推定することもできる。オクルージョン領域の推定についても同様である。また、縦に配置されたカメラで撮影した複数の視点画像に対しても、同様に視差やオクルージョン領域を推定することも可能である。その際、各実施形態における処理をそのまま適用することも、縦方向と横方向の処理を入れ替えて適用することも可能である。さらには、斜めに配置されたカメラで撮影した複数の視点画像に対しても、同様に視差やオクルージョン領域を推定することも可能である。また、本発明の各実施形態において、動的計画法を適用して視差推定に係る各行の処理を行うことを例に挙げたが、これに限ったものではなく、例えばブロックマッチング、グラフカットなどの他の手法を適用してもよい。本発明は、上記の各実施形態に限定されるものではない。これらの実施形態にそれぞれ開示された技術的手段を適宜組合せて得られる実施形態についても本発明の技術的範囲に含まれる。また、本発明の第一~第三の実施形態に係る立体画像処理装置は、各種レコーダや各種記録メディア再生装置などの映像出力機器に組み込むことができる。さらに、本発明は、このような立体画像処理装置に表示デバイスを外部接続又は内部に搭載したテレビ装置やモニタ装置などの立体画像表示装置としての形態も採り得る。また、図3,図6で例示した立体画像処理装置10,20をはじめ、各実施形態に係る立体画像処理装置は、例えばマイクロプロセッサ又はDSP:DigitalSignalProcessor、メモリ、バス、インターフェイス、周辺装置などのハードウェアと、これらのハードウェア上にて実行可能なソフトウェアとにより実現できる。上記ハードウェアの一部又は全部は集積回路/ICIntegratedCircuitチップセットとして搭載することができ、その場合、上記ソフトウェアは上記メモリに記憶しておければよい。また、本発明の各構成要素の全てをハードウェアで構成してもよく、その場合についても同様に、そのハードウェアの一部又は全部を集積回路/ICチップセットとして搭載することも可能である。また、各実施形態に係る立体画像処理装置は単に、CPUCentralProcessingUnit、作業領域としてのRAMRandomAccessMemory、制御用のプログラムの格納領域としてのROMReadOnlyMemoryやEEPROMElectricallyErasableProgrammableROM等の記憶装置などで構成することもできる。その場合、上記制御用のプログラムは、本発明に係る処理を実行するための後述の立体画像処理プログラムを含むことになる。この立体画像処理プログラムは、PC内に立体画像処理用のアプリケーションソフトとして組み込み、PCを立体画像処理装置として機能させることもできる。また、この立体画像処理プログラムは、クライアントPCから実行可能な状態でWebサーバ等の外部サーバに格納されていてもよい。以上、本発明に係る立体画像処理装置を中心に説明したが、本発明は、この立体画像処理装置における制御の流れを例示したように、立体画像処理方法としての形態も採り得る。この立体画像処理方法は、複数の視点画像から視点画像間の視差を推定する立体画像処理方法であって、評価値算出部が、視点画像間での画素の類似度及び視点画像間の視差の連続性に関する評価値を算出するステップと、視差マップ作成部が、その評価値を最小とする、視点画像間の各画素の組合せを、視差として検出して視差マップを作成するステップとを有するものとする。ここで、視差の連続性は、注目画素の視差とその注目画素の横及び縦に隣接する画素の視差との連続性であり、評価値は、評価値を算出するステップにおいて、画素の類似度が高いほど小さく算出され、視差の連続性が高いほど小さく算出されるものとする。その他の応用例については、立体画像処理装置について説明したとおりである。また、本発明は、その立体画像処理方法をコンピュータにより実行させるための立体画像処理プログラムとしての形態も採り得る。つまり、この立体画像処理プログラムは、コンピュータに、複数の視点画像から視点画像間の視差を推定する立体画像処理を実行させるためのプログラムである。そして、この立体画像処理は、視点画像間での画素の類似度及び視点画像間の視差の連続性に関する評価値を算出するステップと、その評価値を最小とする、視点画像間の各画素の組合せを、視差として検出して視差マップを作成するステップと、を含むものとする。ここで、視差の連続性は、注目画素の視差とその注目画素の横及び縦に隣接する画素の視差との連続性であり、評価値は、評価値を算出するステップにおいて、画素の類似度が高いほど小さく算出され、視差の連続性が高いほど小さく算出されるその他の応用例については、立体画像処理装置について説明したとおりである。また、その立体画像処理プログラムをコンピュータにより読み取り可能な記録媒体に記録したプログラム記録媒体としての形態についても容易に理解することができる。このコンピュータとしては、上述したように、汎用のPCに限らず、マイクロコンピュータやプログラム可能な汎用の集積回路/チップセットなど、様々な形態のコンピュータが適用できる。また、このプログラムは、可搬の記録媒体を介して流通させるに限らず、インターネット等のネットワークを介して、また放送波を介して流通させることもできる。ネットワークを介して受信するとは、外部サーバの記憶装置などに記録されたプログラムを受信することを指す。"