(57)【要約】【課題】画素値の変動が激しい画像や、コントラストが低い画像であっても、精度よく画素ずれ量を検出することができる画像処理装置を得ることを目的とする。【解決手段】相関値演算処理部９が当該フレームの変位量Ｋh（ｘ），Ｋv（ｙ）及び曲率量Ｌh（ｘ），Ｌv（ｙ）と、１フレーム前の変位量Ｋ’h（ｘ），Ｋ’v（ｙ）及び曲率量Ｌ’h（ｘ），Ｌ’v（ｙ）とを用いて、フレーム間の相関値Ｒh（ｘd），Ｒv（ｙd）を演算し、その相関値Ｒh（ｘd），Ｒv（ｙd）からフレーム間の画素ずれ量を検出する。【選択図】図１
2008079007.tif 000002
【特許請求の範囲】
【請求項１】
画像の画素値を示す画像信号を入力する画像信号入力手段と、上記画像信号入力手段により入力された画像信号のうち、水平方向の座標が一致する画像信号同士の画素値を積算して、その画素値の積算結果である水平積算値を出力するとともに、垂直方向の座標が一致する画像信号同士の画素値を積算して、その画素値の積算結果である垂直積算値を出力する画素値積算手段と、上記画素値積算手段から出力された複数の水平積算値を比較して、水平積算値の変位量及び曲率量を検出するとともに、上記画素値積算手段から出力された複数の垂直積算値を比較して、垂直積算値の変位量及び曲率量を検出する検出処理手段と、上記検出処理手段により検出された変位量及び曲率量から２枚の画像間の相関値を演算し、上記相関値から２枚の画像間の画素ずれ量を検出する画素ずれ量検出手段とを備えた画像処理装置。
【請求項２】
検出処理手段は、画素値積算手段から出力された複数の水平積算値のうち、水平方向の間隔が一定の関係にある２つの水平積算値の大小関係を比較し、大小関係の比較結果から水平積算値の変位量を検出する第１の変位量検出手段と、上記画素値積算手段から出力された複数の水平積算値のうち、水平方向の間隔が一定の関係にある３つの水平積算値の大小関係を比較し、大小関係の比較結果から水平積算値の曲率量を検出する第１の曲率量検出手段と、上記画素値積算手段から出力された複数の垂直積算値のうち、垂直方向の間隔が一定の関係にある２つの垂直積算値の大小関係を比較し、大小関係の比較結果から垂直積算値の変位量を検出する第２の変位量検出手段と、上記画素値積算手段から出力された複数の垂直積算値のうち、垂直方向の間隔が一定の関係にある３つの垂直積算値の大小関係を比較し、大小関係の比較結果から垂直積算値の曲率量を検出する第２の曲率量検出手段とから構成されていることを特徴とする請求項１記載の画像処理装置。
【請求項３】
検出処理手段は、大小関係を比較する２つの水平積算値の水平方向の間隔が相互に異なる第１の変位量検出手段を複数個備え、大小関係を比較する３つの水平積算値の水平方向の間隔が相互に異なる第１の曲率量検出手段を複数個備え、大小関係を比較する２つの垂直積算値の垂直方向の間隔が相互に異なる第２の変位量検出手段を複数個備え、大小関係を比較する３つの垂直積算値の垂直方向の間隔が相互に異なる第２の曲率量検出手段を複数個備えていることを特徴とする請求項２記載の画像処理装置。
【請求項４】
画素ずれ量検出手段は、検出処理手段により今回の検出処理で検出された変位量及び曲率量と、上記検出処理手段により前回の検出処理で検出された変位量及び曲率量とを用いて、２枚の画像間の相関値を演算し、上記相関値が最小になる画素ずれ量を２枚の画像間の画素ずれ量として出力することを特徴とする請求項１から請求項３のうちのいずれか１項記載の画像処理装置。
【請求項５】
画素ずれ量検出手段は、２枚の画像間の相関値を演算する際、検出処理手段により今回の検出処理で検出された変位量及び曲率量と、上記検出処理手段により前回の検出処理で検出された変位量及び曲率量との間に重み付けを行うことを特徴とする請求項４記載の画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、２枚の画像間の画素ずれ量を検出して、画素ずれのない安定した画像の出力を可能にする画像処理装置に関するものである。
【背景技術】
【０００２】
従来の画像処理装置は、次のような手法を用いて、２枚の画像間の画素ずれ量を検出するようにしている。
即ち、画像処理装置は、画像センサにより撮影された２枚の画像の射影データを入力すると、２つの射影データの相関値を演算して、２枚のフレーム画像間の動きベクトルを算出し、その動きベクトルを画素ずれ量として検出する（例えば、特許文献１を参照）。
【０００３】
最近では、上記の画像処理装置の他に、相関値を演算する際の演算量を削減することができる画像処理装置が開発されている。
当該画像処理装置は、２枚のフレーム画像間の相関値を演算する際、射影データの傾き情報（例えば、画像の傾斜部の重心座標、傾斜開始座標、傾斜終了座標など）を抽出し、複数の傾き情報を時系列的に比較して、その比較結果から画素ずれ量を検出するようにしている（例えば、特許文献２を参照）。
【０００４】
【特許文献１】特許３７７０２７１号公報（第２図）
【特許文献２】特開平６−３８０９１号公報（段落番号［００１６］から［００３８］、図１）
【発明の開示】
【発明が解決しようとする課題】
【０００５】
従来の画像処理装置は以上のように構成されているので、射影データの傾き情報を抽出することができれば、その傾き情報から画素ずれ量を検出することができる。しかし、射影データの傾き情報を抽出する度に、その傾き情報として、画像の傾斜部の重心座標、傾斜開始座標、傾斜終了座標などのデータをメモリに格納する必要があり、そのメモリに格納するデータ量は常に一定ではないため、画素値変動が激しい被写体においては、傾き情報が狭いエリアのデータとなってしまって、十分な画素ずれ量の検出精度が得られないことがある課題があった。
また、コントラストが低く、画像の傾斜部として認識することが困難な微小傾斜が含まれている画像に関しては、画素ずれ量を検出することができない課題があった。
【０００６】
この発明は上記のような課題を解決するためになされたもので、画素値の変動が激しい画像や、コントラストが低い画像であっても、精度よく画素ずれ量を検出することができる画像処理装置を得ることを目的とする。
【課題を解決するための手段】
【０００７】
この発明に係る画像処理装置は、画像信号入力手段により入力された画像信号のうち、水平方向の座標が一致する画像信号同士の画素値を積算して、その画素値の積算結果である水平積算値を出力するとともに、垂直方向の座標が一致する画像信号同士の画素値を積算して、その画素値の積算結果である垂直積算値を出力する画素値積算手段と、画素値積算手段から出力された複数の水平積算値を比較して、水平積算値の変位量及び曲率量を検出するとともに、画素値積算手段から出力された複数の垂直積算値を比較して、垂直積算値の変位量及び曲率量を検出する検出処理手段とを設け、画素ずれ量検出手段が検出処理手段により検出された変位量及び曲率量から２枚の画像間の相関値を演算し、その相関値から２枚の画像間の画素ずれ量を検出するようにしたものである。
【発明の効果】
【０００８】
この発明によれば、画像信号入力手段により入力された画像信号のうち、水平方向の座標が一致する画像信号同士の画素値を積算して、その画素値の積算結果である水平積算値を出力するとともに、垂直方向の座標が一致する画像信号同士の画素値を積算して、その画素値の積算結果である垂直積算値を出力する画素値積算手段と、画素値積算手段から出力された複数の水平積算値を比較して、水平積算値の変位量及び曲率量を検出するとともに、画素値積算手段から出力された複数の垂直積算値を比較して、垂直積算値の変位量及び曲率量を検出する検出処理手段とを設け、画素ずれ量検出手段が検出処理手段により検出された変位量及び曲率量から２枚の画像間の相関値を演算し、その相関値から２枚の画像間の画素ずれ量を検出するように構成したので、画素値の変動が激しい画像や、コントラストが低い画像であっても、精度よく画素ずれ量を検出することができる効果がある。
【発明を実施するための最良の形態】
【０００９】
実施の形態１．
図１はこの発明の実施の形態１による画像処理装置を示す構成図であり、図において、画像入力部１は画像センサにより撮影された画像の射影データ（画像の輝度値（画素値）を示す画像信号）を入力する処理を実施する。なお、画像入力部１は画像信号入力手段を構成している。
水平データ積算処理部２は画像入力部１により入力された射影データのうち、水平方向の座標ｘが一致する射影データ同士の輝度値Ｐ（ｘ，ｙ）を積算して、その輝度値Ｐ（ｘ，ｙ）の積算結果である水平積算値Ｓh（ｘ）を出力する処理を実施する。
垂直データ積算処理部３は画像入力部１により入力された射影データのうち、垂直方向の座標ｙが一致する射影データ同士の輝度値Ｐ（ｘ，ｙ）を積算して、その輝度値Ｐ（ｘ，ｙ）の積算結果である垂直積算値Ｓv（ｙ）を出力する処理を実施する。
なお、水平データ積算処理部２及び垂直データ積算処理部３から画素値積算手段が構成されている。
【００１０】
変位検出処理部４は水平データ積算処理部２から出力された複数の水平積算値Ｓh（ｘ）のうち、水平方向の間隔が一定の関係にある２つの水平積算値の大小関係（例えば、Ｓh（ｘ）とＳh（ｘ＋ａ）の大小関係）を比較し、大小関係の比較結果から水平積算値の変位量Ｋh（ｘ）を検出する処理を実施する。なお、変位検出処理部４は第１の変位量検出手段を構成している。
曲率検出処理部５は水平データ積算処理部２から出力された複数の水平積算値Ｓh（ｘ）のうち、水平方向の間隔が一定の関係にある３つの水平積算値の大小関係（例えば、Ｓh（ｘ）×２と、Ｓh（ｘ−ｂ）＋Ｓh（ｘ＋ｂ）の大小関係）を比較し、大小関係の比較結果から中央の位置にある水平積算値Ｓh（ｘ）の曲率量Ｌh（ｘ）を検出する処理を実施する。なお、曲率検出処理部５は第１の曲率量検出手段を構成している。
【００１１】
変位検出処理部６は垂直データ積算処理部３から出力された複数の垂直積算値Ｓv（ｙ）のうち、垂直方向の間隔が一定の関係にある２つの垂直積算値の大小関係（例えば、Ｓv（ｙ）とＳv（ｙ＋ｃ）の大小関係）を比較し、大小関係の比較結果から垂直積算値の変位量Ｋv（ｙ）を検出する処理を実施する。なお、変位検出処理部６は第２の変位量検出手段を構成している。
曲率検出処理部７は垂直データ積算処理部３から出力された複数の垂直積算値Ｓv（ｙ）のうち、垂直方向の間隔が一定の関係にある３つの垂直積算値の大小関係（例えば、Ｓv（ｙ）×２と、Ｓv（ｙ−ｄ）＋Ｓv（ｙ＋ｄ）の大小関係）を比較し、大小関係の比較結果から中央の位置にある垂直積算値Ｓv（ｙ）の曲率量Ｌv（ｙ）を検出する処理を実施する。なお、曲率検出処理部７は第２の曲率量検出手段を構成しており、変位検出処理部４、曲率検出処理部５、変位検出処理部６及び曲率検出処理部７から検出処理手段が構成されている。
【００１２】
特徴量バッファ８は変位検出処理部４，６により検出された変位量Ｋh（ｘ），Ｋv（ｙ）を記録するとともに、曲率検出処理部５，７により検出された曲率量Ｌh（ｘ），Ｌv（ｙ）を記録するメモリである。
相関値演算処理部９は変位検出処理部４，６により今回の処理で検出された変位量Ｋh（ｘ），Ｋv（ｙ）と、曲率検出処理部５，７により今回の処理で検出された曲率量Ｌh（ｘ），Ｌv（ｙ）と、特徴量バッファ８に記録されている変位量Ｋ’h（ｘ），Ｋ’v（ｙ）（変位検出処理部４，６により前回の処理で検出された変位量Ｋh（ｘ），Ｋv（ｙ））と、特徴量バッファ８に記録されている曲率量Ｌ’h（ｘ），Ｌ’v（ｙ）（曲率検出処理部５，７により前回の処理で検出された曲率量Ｌh（ｘ），Ｌv（ｙ））とを用いて、２枚の画像間の相関値Ｒh（ｘd），Ｒv（ｙd）を演算し、相関値Ｒh（ｘd），Ｒv（ｙd）が最小になる画素ずれ量ｘd，ｙdを２枚の画像間の画素ずれ量として出力する処理を実施する。
なお、特徴量バッファ８及び相関値演算処理部９から画素ずれ量検出手段が構成されている。
【００１３】
画像出力部１０は画像入力部１から画像の射影データを受けると、相関値演算処理部９から出力された２枚の画像間の画素ずれ量を打ち消す方向に画像の切り出し位置を変更して、画素ずれのない安定した画像を出力する処理を実施する。
【００１４】
次に動作について説明する。
画像入力部１は、画像センサにより撮影された画像の射影データを入力する毎に、その射影データを水平データ積算処理部２、垂直データ積算処理部３及び画像出力部１０に出力する。
【００１５】
水平データ積算処理部２は、画像入力部１から画像の射影データを受けると、画像の射影データのうち、水平方向の座標ｘが一致する射影データ同士の輝度値Ｐ（ｘ，ｙ）を積算して、その輝度値Ｐ（ｘ，ｙ）の積算結果である水平積算値Ｓh（ｘ）を変位検出処理部４及び曲率検出処理部５に出力する。
即ち、水平データ積算処理部２は、下記の式（１）に示すような積算演算を実施して、水平積算値Ｓh（ｘ）を出力する。
【数１】
2008079007.tif 000003
ただし、ｘは水平方向の座標、ｙは垂直方向の座標、ｙmaxは射影データにおける垂直方向の座標の最大値、Ｐ（ｘ，ｙ）は座標（ｘ，ｙ）に位置する画素の輝度値、Ｓh（ｘ）は輝度値Ｐ（ｘ，ｙ）の積算結果である。
【００１６】
垂直データ積算処理部３は、画像入力部１から画像の射影データを受けると、画像の射影データのうち、垂直方向の座標ｙが一致する射影データ同士の輝度値Ｐ（ｘ，ｙ）を積算して、その輝度値Ｐ（ｘ，ｙ）の積算結果である垂直積算値Ｓv（ｙ）を変位検出処理部６及び曲率検出処理部７に出力する。
即ち、垂直データ積算処理部３は、下記の式（２）に示すような積算演算を実施して、水平積算値Ｓv（ｙ）を出力する。
【数２】
2008079007.tif 000004
ただし、ｘmaxは射影データにおける水平方向の座標の最大値、Ｓv（ｙ）は輝度値Ｐ（ｘ，ｙ）の積算結果である。
【００１７】
変位検出処理部４は、水平データ積算処理部２が複数の水平積算値Ｓh（ｘ）を出力すると、複数の水平積算値Ｓh（ｘ）のうち、水平方向の間隔が一定の関係にある２つの水平積算値の大小関係を比較する。
即ち、変位検出処理部４は、水平積算値Ｓh（ｘ）と水平積算値Ｓh（ｘ＋ａ）の大小関係を比較する。
そして、変位検出処理部４は、下記の式（３）に示すように、水平積算値Ｓh（ｘ）が水平積算値Ｓh（ｘ＋ａ）以上であれば、水平積算値の変位量Ｋh（ｘ）として“１”を出力する。
一方、水平積算値Ｓh（ｘ）が水平積算値Ｓh（ｘ＋ａ）より小さければ、水平積算値の変位量Ｋh（ｘ）として“０”を出力する。
【００１８】
【数３】
2008079007.tif 000005
ただし、ａは水平積算値Ｓh（ｘ）と水平積算値Ｓh（ｘ＋ａ）の水平方向の間隔を示す定数である。
なお、式（３）の演算においては、対象座標値の水平積算値Ｓh（ｘ）と一定の距離がある座標値の水平積算値Ｓh（ｘ＋ａ）を参照するため、対象座標値によっては、参照する水平積算値Ｓh（ｘ＋ａ）が存在しない場合がある。そのような場合には、存在する最も近い座標値の水平積算値を参照するようにしてもよいし、水平積算値の変位量Ｋh（ｘ）を“０”の固定値としてもよい。
【００１９】
ここで、図２は定数ａが“１”である場合の変位検出処理部４の動作例を示す説明図である。
例えば、Ｓh（０）→Ｓh（１）や、Ｓh（３）→Ｓh（４）の局面では、水平積算値の変位量Ｋh（０），Ｋh（３）として“０”が出力され、Ｓh（１）→Ｓh（２）や、Ｓh（２）→Ｓh（３）の局面では、水平積算値の変位量Ｋh（１），Ｋh（２）として“１”が出力されている。
これにより、変位検出処理部４から垂直画像サイズ分のデータ列として、水平積算値の変位量Ｋh（ｘ）が出力されることになる。
【００２０】
曲率検出処理部５は、水平データ積算処理部２が複数の水平積算値Ｓh（ｘ）を出力すると、複数の水平積算値Ｓh（ｘ）のうち、水平方向の間隔が一定の関係にある３つの水平積算値の大小関係を比較する。
即ち、曲率検出処理部５は、Ｓh（ｘ）×２と、Ｓh（ｘ−ｂ）＋Ｓh（ｘ＋ｂ）の大小関係を比較する。
そして、曲率検出処理部５は、下記の式（４）に示すように、中央の位置にある水平積算値の２倍であるＳh（ｘ）×２が、水平積算値Ｓh（ｘ−ｂ）と水平積算値Ｓh（ｘ＋ｂ）の和であるＳh（ｘ−ｂ）＋Ｓh（ｘ＋ｂ）以上であれば、水平積算値Ｓh（ｘ）の曲率量Ｌh（ｘ）として“１”を出力する。
一方、Ｓh（ｘ）×２が、Ｓh（ｘ−ｂ）＋Ｓh（ｘ＋ｂ）より小さければ、水平積算値Ｓh（ｘ）の曲率量Ｌh（ｘ）として“０”を出力する。
【００２１】
【数４】
2008079007.tif 000006
ただし、ｂは水平積算値Ｓh（ｘ）と水平積算値Ｓh（ｘ−ｂ），Ｓh（ｘ＋ｂ）の水平方向の間隔を示す定数である。
なお、式（４）の演算においては、対象座標値の水平積算値Ｓh（ｘ）と一定の距離がある座標値の水平積算値Ｓh（ｘ−ｂ），Ｓh（ｘ＋ｂ）を参照するため、対象座標値によっては、参照する水平積算値Ｓh（ｘ−ｂ），Ｓh（ｘ＋ｂ）が存在しない場合がある。そのような場合には、存在する最も近い座標値の水平積算値を参照するようにしてもよいし、水平積算値Ｓh（ｘ）の曲率量Ｌh（ｘ）を“０”の固定値としてもよい。
【００２２】
変位検出処理部６は、垂直データ積算処理部３が複数の垂直積算値Ｓv（ｙ）を出力すると、複数の垂直積算値Ｓv（ｙ）のうち、垂直方向の間隔が一定の関係にある２つの垂直積算値の大小関係を比較する。
即ち、変位検出処理部６は、垂直積算値Ｓv（ｙ）と垂直積算値Ｓv（ｙ＋ｃ）の大小関係を比較する。
そして、変位検出処理部６は、下記の式（５）に示すように、垂直積算値Ｓv（ｙ）が垂直積算値Ｓv（ｙ＋ｃ）以上であれば、垂直積算値の変位量Ｋv（ｙ）として“１”を出力する。
一方、垂直積算値Ｓv（ｙ）が垂直積算値Ｓv（ｙ＋ｃ）より小さければ、垂直積算値の変位量Ｋv（ｙ）として“０”を出力する。
【００２３】
【数５】
2008079007.tif 000007
ただし、ｃは垂直積算値Ｓv（ｙ）と垂直積算値Ｓv（ｙ＋ｃ）の垂直方向の間隔を示す定数である。
なお、式（５）の演算においては、対象座標値の垂直積算値Ｓv（ｙ）と一定の距離がある座標値の垂直積算値Ｓv（ｙ＋ｃ）を参照するため、対象座標値によっては、参照する垂直積算値Ｓv（ｙ＋ｃ）が存在しない場合がある。そのような場合には、存在する最も近い座標値の垂直積算値を参照するようにしてもよいし、垂直積算値の変位量Ｋv（ｙ）を“０”の固定値としてもよい。
【００２４】
曲率検出処理部７は、垂直データ積算処理部３が複数の垂直積算値Ｓv（ｙ）を出力すると、複数の垂直積算値Ｓv（ｙ）のうち、垂直方向の間隔が一定の関係にある３つの垂直積算値の大小関係を比較する。
即ち、曲率検出処理部７は、Ｓv（ｙ）×２と、Ｓv（ｙ−ｄ）＋Ｓv（ｙ＋ｄ）の大小関係を比較する。
そして、曲率検出処理部７は、下記の式（６）に示すように、中央の位置にある垂直積算値の２倍であるＳv（ｙ）×２が、垂直積算値Ｓv（ｙ−ｄ）と垂直積算値Ｓv（ｙ＋ｄ）の和であるＳv（ｙ−ｄ）＋Ｓv（ｙ＋ｄ）以上であれば、垂直積算値Ｓv（ｙ）の曲率量Ｌv（ｙ）として“１”を出力する。
一方、Ｓv（ｙ）×２が、Ｓv（ｙ−ｄ）＋Ｓv（ｙ＋ｄ）より小さければ、垂直積算値Ｓv（ｙ）の曲率量Ｌv（ｙ）として“０”を出力する。
【００２５】
【数６】
2008079007.tif 000008
ただし、ｄは垂直積算値Ｓv（ｙ）と垂直積算値Ｓv（ｙ−ｄ），Ｓv（ｙ＋ｄ）の垂直方向の間隔を示す定数である。
なお、式（６）の演算においては、対象座標値の垂直積算値Ｓv（ｙ）と一定の距離がある座標値の垂直積算値Ｓv（ｙ−ｄ），Ｓv（ｙ＋ｄ）を参照するため、対象座標値によっては、参照する垂直積算値Ｓv（ｙ−ｄ），Ｓv（ｙ＋ｄ）が存在しない場合がある。そのような場合には、存在する最も近い座標値の垂直積算値を参照するようにしてもよいし、垂直積算値Ｓv（ｙ）の曲率量Ｌv（ｙ）を“０”の固定値としてもよい。
【００２６】
上記のようにして、変位検出処理部４，６及び曲率検出処理部５，７により検出された変位量Ｋh（ｘ），Ｋv（ｙ）と曲率量Ｌh（ｘ），Ｌv（ｙ）は、相関値演算処理部９に出力されるほか、特徴量バッファ８に出力されて、変位量Ｋh（ｘ），Ｋv（ｙ）と曲率量Ｌh（ｘ），Ｌv（ｙ）が特徴量バッファ８に記録される。
なお、特徴量バッファ８に記録される変位量Ｋh（ｘ），Ｋv（ｙ）と曲率量Ｌh（ｘ），Ｌv（ｙ）は、相関値演算処理部９における次回の相関演算処理に備えるものである。
【００２７】
相関値演算処理部９は、変位検出処理部４，６から今回の処理で検出された変位量Ｋh（ｘ），Ｋv（ｙ）を受け、曲率検出処理部５，７から今回の処理で検出された曲率量Ｌh（ｘ），Ｌv（ｙ）を受けると、特徴量バッファ８から変位検出処理部４，６により前回の処理で検出された変位量Ｋh（ｘ），Ｋv（ｙ）を一つ前のフレームの画像に基づく変位量Ｋ’h（ｘ），Ｋ’v（ｙ）として読み込み、また、特徴量バッファ８から曲率検出処理部５，７により前回の処理で検出された曲率量Ｌh（ｘ），Ｌv（ｙ）を一つ前のフレームの画像に基づく曲率量Ｌ’h（ｘ），Ｌ’v（ｙ）として読み込みを行う。
そして、相関値演算処理部９は、変位量Ｋh（ｘ），Ｋv（ｙ），Ｋ’h（ｘ），Ｋ’v（ｙ）と、曲率量Ｌh（ｘ），Ｌv（ｙ），Ｌ’h（ｘ），Ｌ’v（ｙ）とを用いて、２枚の画像間（フレーム間）の相関値Ｒh（ｘd），Ｒv（ｙd）を演算する。
【００２８】
即ち、相関値演算処理部９は、下記の式（７）に示すような相関値演算を実施することにより、比較対象の画素を１画素ずつずらしながら、水平方向の画素ずれ量がｘdであるとする場合のフレーム間の相関値Ｒh（ｘd）を演算する。
【数７】
2008079007.tif 000009
ただし、ｅは画素ずれ量の最大検出値を示す定数である。
【００２９】
また、相関値演算処理部９は、下記の式（８）に示すような相関値演算を実施することにより、比較対象の画素を１画素ずつずらしながら、垂直方向の画素ずれ量がｙdであるとする場合のフレーム間の相関値Ｒv（ｙd）を演算する。
【数８】
2008079007.tif 000010
ただし、ｆは画素ずれ量の最大検出値を示す定数である。
【００３０】
相関値演算処理部９は、上記の相関演算を−ｅ≦ｘd≦ｅ、−ｆ≦ｙd≦ｆのそれぞれの範囲で実施して、フレーム間の相関値Ｒh（ｘd），Ｒv（ｙd）を演算すると、フレーム間の相関値Ｒh（ｘd），Ｒv（ｙd）が最小になる画素ずれ量ｘd，ｙdをフレーム間の画素ずれ量として出力する。
ここで、図３はｘdの値を−１、０、１とする場合の相関値演算処理部１２の動作例を示す説明図である。
図３の例では、ｘd＝−１  →  Ｒh（ｘd）＝１０、ｘd＝０  →  Ｒh（ｘd）＝９、ｘd＝１  →  Ｒh（ｘd）＝２であり、ｘdの値が１のとき、フレーム間の相関値Ｒh（ｘd）が最小になるので、水平方向の画素ずれ量ｘdとして“１”が出力される。
【００３１】
画像出力部１０は、画像入力部１から画像の射影データを受け、相関値演算処理部９からフレーム間の画素ずれ量を受けると、その画素ずれ量を打ち消す方向に画像の切り出し位置を変更して、画素ずれのない安定した画像を出力する。
例えば、水平方向の画素ずれ量ｘdが“１”であれば、画像の切り出し位置を前のフレームより１画素だけ左の位置に変更する。
また、垂直方向の画素ずれ量ｙdが“−２”であれば、画像の切り出し位置を前のフレームより２画素だけ下の位置に変更する。
【００３２】
以上で明らかなように、この実施の形態１によれば、画像入力部１により入力された射影データのうち、水平方向の座標ｘが一致する射影データ同士の輝度値Ｐ（ｘ，ｙ）を積算して、その輝度値の積算結果である水平積算値Ｓh（ｘ）を出力する水平データ積算処理部２と、垂直方向の座標ｙが一致する射影データ同士の輝度値Ｐ（ｘ，ｙ）を積算して、その輝度値の積算結果である垂直積算値Ｓv（ｙ）を出力する垂直データ積算処理部３と、水平データ積算処理部２又は垂直データ積算処理部３から出力された複数の水平積算値Ｓh（ｘ）又は複数の垂直積算値Ｓv（ｙ）を比較して、水平積算値の変位量Ｋh（ｘ）又は垂直積算値の変位量Ｋv（ｙ）を検出する変位検出処理部４，６と、水平データ積算処理部２又は垂直データ積算処理部３から出力された複数の水平積算値Ｓh（ｘ）又は複数の垂直積算値Ｓv（ｙ）を比較して、水平積算値Ｓh（ｘ）の曲率量Ｌh（ｘ）又は垂直積算値Ｓv（ｙ）の曲率量Ｌv（ｙ）を検出する曲率検出処理部５，７とを設け、相関値演算処理部９が当該フレームの変位量Ｋh（ｘ），Ｋv（ｙ）及び曲率量Ｌh（ｘ），Ｌv（ｙ）と、１フレーム前の変位量Ｋ’h（ｘ），Ｋ’v（ｙ）及び曲率量Ｌ’h（ｘ），Ｌ’v（ｙ）とを用いて、フレーム間の相関値Ｒh（ｘd），Ｒv（ｙd）を演算し、その相関値Ｒh（ｘd），Ｒv（ｙd）からフレーム間の画素ずれ量を検出するように構成したので、画素値の変動が激しい画像や、コントラストが低い画像であっても、精度よく画素ずれ量を検出することができる効果を奏する。
【００３３】
即ち、この実施の形態１では、水平積算値Ｓh（ｘ）及び垂直積算値Ｓv（ｙ）から変位量Ｋh（ｘ），Ｋv（ｙ）と曲率量Ｌh（ｘ），Ｌv（ｙ）を検出し、射影データの各ライン当り１ビットの情報として、その変位量Ｋh（ｘ），Ｋv（ｙ）と曲率量Ｌh（ｘ），Ｌv（ｙ）を特徴量バッファ８に記録するようにしているので、特徴量として保持するデータ量を大きく削減することが可能である。
また、フレーム間で比較するデータ量が削減されることによって、相関値演算処理部９の相関値演算処理では、各ライン毎に１ビット値の一致もしくは不一致を判定するだけでよくなり、相関値演算の演算量自体も大きく削減することができる。
また、特徴量として互いに性質の異なる２種類の値を検出し、その組み合わせによって相関値の演算を行うことにより、個々の値のみで相関値を算出する場合よりも高精度に画素ずれ量を算出することが可能となる。
【００３４】
実施の形態２．
図４はこの発明の実施の形態２による画像処理装置を示す構成図であり、図において、図１と同一符号は同一または相当部分を示すので説明を省略する。
変位検出処理部４ａは水平データ積算処理部２から出力された複数の水平積算値Ｓh（ｘ）のうち、水平方向の間隔が一定の関係にある２つの水平積算値の大小関係（例えば、Ｓh（ｘ）とＳh（ｘ＋ａ1）の大小関係）を比較し、大小関係の比較結果から水平積算値の変位量Ｋh1（ｘ）を検出する処理を実施する。
変位検出処理部４ｂは水平データ積算処理部２から出力された複数の水平積算値Ｓh（ｘ）のうち、水平方向の間隔が一定の関係にある２つの水平積算値の大小関係（例えば、Ｓh（ｘ）とＳh（ｘ＋ａ2）の大小関係）を比較し、大小関係の比較結果から水平積算値の変位量Ｋh2（ｘ）を検出する処理を実施する。
なお、変位検出処理部４ａ，４ｂは第１の変位量検出手段を構成している。
【００３５】
曲率検出処理部５ａは水平データ積算処理部２から出力された複数の水平積算値Ｓh（ｘ）のうち、水平方向の間隔が一定の関係にある３つの水平積算値の大小関係（例えば、Ｓh（ｘ）×２と、Ｓh（ｘ−ｂ1）＋Ｓh（ｘ＋ｂ1）の大小関係）を比較し、大小関係の比較結果から中央の位置にある水平積算値Ｓh（ｘ）の曲率量Ｌh1（ｘ）を検出する処理を実施する。
曲率検出処理部５ｂは水平データ積算処理部２から出力された複数の水平積算値Ｓh（ｘ）のうち、水平方向の間隔が一定の関係にある３つの水平積算値の大小関係（例えば、Ｓh（ｘ）×２と、Ｓh（ｘ−ｂ2）＋Ｓh（ｘ＋ｂ2）の大小関係）を比較し、大小関係の比較結果から中央の位置にある水平積算値Ｓh（ｘ）の曲率量Ｌh2（ｘ）を検出する処理を実施する。
なお、曲率検出処理部５ａ，５ｂは第１の曲率量検出手段を構成している。
【００３６】
変位検出処理部６ａは垂直データ積算処理部３から出力された複数の垂直積算値Ｓv（ｙ）のうち、垂直方向の間隔が一定の関係にある２つの垂直積算値の大小関係（例えば、Ｓv（ｙ）とＳv（ｙ＋ｃ1）の大小関係）を比較し、大小関係の比較結果から垂直積算値の変位量Ｋv1（ｙ）を検出する処理を実施する。
変位検出処理部６ｂは垂直データ積算処理部３から出力された複数の垂直積算値Ｓv（ｙ）のうち、垂直方向の間隔が一定の関係にある２つの垂直積算値の大小関係（例えば、Ｓv（ｙ）とＳv（ｙ＋ｃ2）の大小関係）を比較し、大小関係の比較結果から垂直積算値の変位量Ｋv2（ｙ）を検出する処理を実施する。
なお、変位検出処理部６ａ，６ｂは第２の変位量検出手段を構成している。
【００３７】
曲率検出処理部７ａは垂直データ積算処理部３から出力された複数の垂直積算値Ｓv（ｙ）のうち、垂直方向の間隔が一定の関係にある３つの垂直積算値の大小関係（例えば、Ｓv（ｙ）×２と、Ｓv（ｙ−ｄ1）＋Ｓv（ｙ＋ｄ1）の大小関係）を比較し、大小関係の比較結果から中央の位置にある垂直積算値Ｓv（ｙ）の曲率量Ｌv1（ｙ）を検出する処理を実施する。
曲率検出処理部７ｂは垂直データ積算処理部３から出力された複数の垂直積算値Ｓv（ｙ）のうち、垂直方向の間隔が一定の関係にある３つの垂直積算値の大小関係（例えば、Ｓv（ｙ）×２と、Ｓv（ｙ−ｄ2）＋Ｓv（ｙ＋ｄ2）の大小関係）を比較し、大小関係の比較結果から中央の位置にある垂直積算値Ｓv（ｙ）の曲率量Ｌv2（ｙ）を検出する処理を実施する。
なお、曲率検出処理部７ａ，７ｂは第２の曲率量検出手段を構成しており、変位検出処理部４ａ，４ｂ、曲率検出処理部５ａ，５ｂ、変位検出処理部６ａ，６ｂ及び曲率検出処理部７ａ，７ｂから検出処理手段が構成されている。
【００３８】
特徴量バッファ１１は変位検出処理部４ａ，４ｂ，６ａ，６ｂにより検出された変位量Ｋh1（ｘ），Ｋh2（ｘ），Ｋv1（ｙ），Ｋv2（ｙ）を記録するとともに、曲率検出処理部５ａ，５ｂ，７ａ，７ｂにより検出された曲率量Ｌh1（ｘ），Ｌh2（ｘ），Ｌv1（ｙ），Ｌv2（ｙ）を記録するメモリである。
相関値演算処理部１２は変位検出処理部４ａ，４ｂ，６ａ，６ｂにより今回の処理で検出された変位量Ｋh1（ｘ），Ｋh2（ｘ），Ｋv1（ｙ），Ｋv2（ｙ）と、曲率検出処理部５ａ，５ｂ，７ａ，７ｂにより今回の処理で検出された曲率量Ｌh1（ｘ），Ｌh2（ｘ），Ｌv1（ｙ），Ｌv2（ｙ）と、特徴量バッファ８に記録されている変位量Ｋ’h1（ｘ），Ｋ’h2（ｘ），Ｋ’v1（ｙ），Ｋ’v2（ｙ）（変位検出処理部４ａ，４ｂ，６ａ，６ｂにより前回の処理で検出された変位量Ｋh1（ｘ），Ｋh2（ｘ），Ｋv1（ｙ），Ｋv2（ｙ））と、特徴量バッファ８に記録されている曲率量Ｌ’h1（ｘ），Ｌ’h2（ｘ），Ｌ’v1（ｙ），Ｌ’v2（ｙ）（曲率検出処理部５ａ，５ｂ，７ａ，７ｂにより前回の処理で検出された曲率量Ｌh1（ｘ），Ｌh2（ｘ），Ｌv1（ｙ），Ｌv2（ｙ））とを用いて、フレーム間の相関値Ｒh（ｘd），Ｒv（ｙd）を演算し、相関値Ｒh（ｘd），Ｒv（ｙd）が最小になる画素ずれ量ｘd，ｙdをフレーム間の画素ずれ量として出力する処理を実施する。
なお、特徴量バッファ１１及び相関値演算処理部１２から画素ずれ量検出手段が構成されている。
【００３９】
次に動作について説明する。
この実施の形態２では、図１の変位検出処理部４，６に相当する検出処理部と、曲率検出処理部５，７に相当する検出処理部をそれぞれ２組ずつ用意している点で上記実施の形態１と相違している。
即ち、この実施の形態２では、図１の変位検出処理部４に相当する検出処理部として変位検出処理部４ａ，４ｂ、図１の曲率検出処理部５に相当する検出処理部として曲率検出処理部５ａ，５ｂ、図１の変位検出処理部６に相当する検出処理部として変位検出処理部６ａ，６ｂ、図１の曲率検出処理部７に相当する検出処理部として曲率検出処理部７ａ，７ｂを用意している。
【００４０】
画像入力部１は、上記実施の形態１と同様に、画像センサにより撮影された画像の射影データを入力する毎に、その射影データを水平データ積算処理部２、垂直データ積算処理部３及び画像出力部１０に出力する。
水平データ積算処理部２は、画像入力部１から画像の射影データを受けると、上記実施の形態１と同様に、画像の射影データのうち、水平方向の座標ｘが一致する射影データ同士の輝度値Ｐ（ｘ，ｙ）を積算して、その輝度値Ｐ（ｘ，ｙ）の積算結果である水平積算値Ｓh（ｘ）を変位検出処理部４ａ，４ｂ及び曲率検出処理部５ａ，５ｂに出力する。
垂直データ積算処理部３は、画像入力部１から画像の射影データを受けると、上記実施の形態１と同様に、画像の射影データのうち、垂直方向の座標ｙが一致する射影データ同士の輝度値Ｐ（ｘ，ｙ）を積算して、その輝度値Ｐ（ｘ，ｙ）の積算結果である垂直積算値Ｓv（ｙ）を変位検出処理部６ａ，６ｂ及び曲率検出処理部７ａ，７ｂに出力する。
【００４１】
変位検出処理部４ａ，４ｂは、水平データ積算処理部２から複数の水平積算値Ｓh（ｘ）を受けると、図１の変位検出処理部４と同様にして、水平積算値の変位量Ｋh1（ｘ），Ｋh2（ｘ）を検出する。
ただし、変位検出処理部４ａ，４ｂは、大小関係を比較する２つの水平積算値Ｓh（ｘ）の水平方向の間隔が相互に異なっている。
即ち、変位検出処理部４ａでは間隔がａ1であり、変位検出処理部４ｂでは間隔がａ2である。
【００４２】
曲率検出処理部５ａ，５ｂは、水平データ積算処理部２から複数の水平積算値Ｓh（ｘ）を受けると、図１の曲率検出処理部５と同様にして、水平積算値Ｓh（ｘ）の曲率量Ｌh1（ｘ），Ｌh2（ｘ）を検出する。
ただし、曲率検出処理部５ａ，５ｂは、大小関係を比較する３つの水平積算値Ｓh（ｘ）の水平方向の間隔が相互に異なっている。
即ち、曲率検出処理部５ａでは間隔がｂ1であり、曲率検出処理部５ｂでは間隔がｂ2である。
【００４３】
変位検出処理部６ａ，６ｂは、垂直データ積算処理部３から複数の垂直積算値Ｓv（ｙ）を受けると、図１の変位検出処理部６と同様にして、垂直積算値の変位量Ｋv1（ｙ），Ｋv2（ｙ）を検出する。
ただし、変位検出処理部６ａ，６ｂは、大小関係を比較する２つの垂直積算値Ｓv（ｙ）の垂直方向の間隔が相互に異なっている。
即ち、変位検出処理部６ａでは間隔がｃ1であり、変位検出処理部６ｂでは間隔がｃ2である。
【００４４】
曲率検出処理部７ａ，７ｂは、垂直データ積算処理部３から複数の垂直積算値Ｓv（ｙ）を受けると、図１の曲率検出処理部７と同様にして、垂直積算値Ｓv（ｙ）の曲率量Ｌv1（ｙ），Ｌv2（ｙ）を検出する。
ただし、曲率検出処理部７ａ，７ｂは、大小関係を比較する３つの垂直積算値Ｓv（ｙ）の垂直方向の間隔が相互に異なっている。
即ち、曲率検出処理部７ａでは間隔がｄ1であり、曲率検出処理部７ｂでは間隔がｄ2である。
【００４５】
上記のようにして、変位検出処理部４ａ，４ｂ，６ａ，６ｂ及び曲率検出処理部５ａ，５ｂ，７ａ，７ｂにより検出された変位量Ｋh1（ｘ），Ｋh2（ｘ），Ｋv1（ｙ），Ｋv2（ｙ）と曲率量Ｌh1（ｘ），Ｌh2（ｘ），Ｌv1（ｙ），Ｌv2（ｙ）は、相関値演算処理部１２に出力されるほか、特徴量バッファ１１に出力されて、変位量Ｋh1（ｘ），Ｋh2（ｘ），Ｋv1（ｙ），Ｋv2（ｙ）と曲率量Ｌh1（ｘ），Ｌh2（ｘ），Ｌv1（ｙ），Ｌv2（ｙ）が特徴量バッファ１１に記録される。
なお、特徴量バッファ１１に記録される変位量Ｋh1（ｘ），Ｋh2（ｘ），Ｋv1（ｙ），Ｋv2（ｙ）と曲率量Ｌh1（ｘ），Ｌh2（ｘ），Ｌv1（ｙ），Ｌv2（ｙ）は、相関値演算処理部１２における次回の相関演算処理に備えるものである。
【００４６】
相関値演算処理部１２は、変位検出処理部４ａ，４ｂ，６ａ，６ｂから今回の処理で検出された変位量Ｋh1（ｘ），Ｋh2（ｘ），Ｋv1（ｙ），Ｋv2（ｙ）を受け、曲率検出処理部５ａ，５ｂ，７ａ，７ｂから今回の処理で検出された曲率量Ｌh1（ｘ），Ｌh2（ｘ），Ｌv1（ｙ），Ｌv2（ｙ）を受けると、特徴量バッファ１１から変位検出処理部４ａ，４ｂ，６ａ，６ｂにより前回の処理で検出された変位量Ｋh1（ｘ），Ｋh2（ｘ），Ｋv1（ｙ），Ｋv2（ｙ）を一つ前のフレームの画像に基づく変位量Ｋ’h1（ｘ），Ｋ’h2（ｘ），Ｋ’v1（ｙ），Ｋ’v2（ｙ）として読み込み、また、特徴量バッファ１１から曲率検出処理部５ａ，５ｂ，７ａ，７ｂにより前回の処理で検出された曲率量Ｌh1（ｘ），Ｌh2（ｘ），Ｌv1（ｙ），Ｌv2（ｙ）を一つ前のフレームの画像に基づく曲率量Ｌ’h1（ｘ），Ｌ’h2（ｘ），Ｌ’v1（ｙ），Ｌ’v2（ｙ）として読み込みを行う。
そして、相関値演算処理部１２は、変位量Ｋh1（ｘ），Ｋh2（ｘ），Ｋv1（ｙ），Ｋv2（ｙ），Ｋ’h1（ｘ），Ｋ’h2（ｘ），Ｋ’v1（ｙ），Ｋ’v2（ｙ）と、曲率量Ｌh1（ｘ），Ｌh2（ｘ），Ｌv1（ｙ），Ｌv2（ｙ），Ｌ’h1（ｘ），Ｌ’h2（ｘ），Ｌ’v1（ｙ），Ｌ’v2（ｙ）とを用いて、フレーム間の相関値Ｒh（ｘd），Ｒv（ｙd）を演算する。
【００４７】
即ち、相関値演算処理部１２は、下記の式（９）に示すような相関値演算を実施することにより、比較対象の画素を１画素ずつずらしながら、水平方向の画素ずれ量がｘdであるとする場合のフレーム間の相関値Ｒh（ｘd）を演算する。
【数９】
2008079007.tif 000011
【００４８】
また、相関値演算処理部１２は、下記の式（１０）に示すような相関値演算を実施することにより、比較対象の画素を１画素ずつずらしながら、垂直方向の画素ずれ量がｙdであるとする場合のフレーム間の相関値Ｒv（ｙd）を演算する。
【数１０】
2008079007.tif 000012
【００４９】
画像出力部１０は、画像入力部１から画像の射影データを受け、相関値演算処理部１２からフレーム間の画素ずれ量を受けると、上記実施の形態１と同様に、その画素ずれ量を打ち消す方向に画像の切り出し位置を変更して、画素ずれのない安定した画像を出力する。
【００５０】
以上で明らかなように、この実施の形態２によれば、大小関係を比較する２つの水平積算値の水平方向の間隔が相互に異なる変位検出処理部４を２組備え、大小関係を比較する３つの水平積算値の水平方向の間隔が相互に異なる曲率検出処理部６を２組備え、大小関係を比較する２つの垂直積算値の垂直方向の間隔が相互に異なる変位検出処理部５を２組備え、大小関係を比較する３つの垂直積算値の垂直方向の間隔が相互に異なる曲率検出処理部７を２組備えるように構成したので、異なる特徴量（変位量、曲率量）の抽出が可能になり、これらの特徴量を組み合わせて使用することで、単独の特徴量では検出精度の低い画像に対しても検出精度を高めることができる効果を奏する。
【００５１】
なお、この実施の形態２では、変位検出処理部及び曲率検出処理部を２組ずつ備えるものについて示したが、要求される検出の精度に応じて、変位検出処理部又は曲率検出処理部のいずれか一方は１組のみ、もしくは、３組以上備えるようにしてもよい。
【００５２】
実施の形態３．
上記実施の形態１では、相関値演算処理部９が当該フレームの変位量Ｋh（ｘ），Ｋv（ｙ）及び曲率量Ｌh（ｘ），Ｌv（ｙ）と、１フレーム前の変位量Ｋ’h（ｘ），Ｋ’v（ｙ）及び曲率量Ｌ’h（ｘ），Ｌ’v（ｙ）とを用いて、フレーム間の相関値Ｒh（ｘd），Ｒv（ｙd）を演算するものについて示したが、相関値演算処理部９がフレーム間の相関値Ｒh（ｘd），Ｒv（ｙd）を演算する際、当該フレームの変位量Ｋh（ｘ），Ｋv（ｙ）及び曲率量Ｌh（ｘ），Ｌv（ｙ）と、１フレーム前の変位量Ｋ’h（ｘ），Ｋ’v（ｙ）及び曲率量Ｌ’h（ｘ），Ｌ’v（ｙ）との間に重み付けを行うようにしてもよい。
【００５３】
具体的には、次のようにして、フレーム間の相関値Ｒh（ｘd），Ｒv（ｙd）を演算する。
即ち、相関値演算処理部９は、下記の式（１１）に示すように、加重倍率係数α，βを乗算することによって重み付けを実施して水平方向の画素ずれ量がｘdであるとする場合のフレーム間の相関値Ｒh（ｘd）を演算する。
【数１１】
2008079007.tif 000013
【００５４】
また、相関値演算処理部９は、下記の式（１２）に示すように、加重倍率係数γ，δを乗算することによって重み付けを実施して垂直方向の画素ずれ量がｙdであるとする場合のフレーム間のＲv（ｙd）を演算する。
【数１２】
2008079007.tif 000014
【００５５】
以上で明らかなように、この実施の形態３によれば、相関値演算処理部９がフレーム間の相関値Ｒh（ｘd），Ｒv（ｙd）を演算する際、当該フレームの変位量Ｋh（ｘ），Ｋv（ｙ）及び曲率量Ｌh（ｘ），Ｌv（ｙ）と、１フレーム前の変位量Ｋ’h（ｘ），Ｋ’v（ｙ）及び曲率量Ｌ’h（ｘ），Ｌ’v（ｙ）との間に重み付けを行うように構成したので、変位検出処理部４，６及び曲率検出処理部５，７における検出精度の違いを加味して、全体の検出精度を高めることができる効果を奏する。
即ち、より検出精度の高い検出処理部の比重を高くして、検出精度の低い検出処理部の比重を低くすることにより、各検出処理部の優先度を設定し、処理精度の向上を図りつつ、より細かく検出精度の調整を行うことが可能になる。
【００５６】
なお、この実施の形態３では、上記実施の形態１における式（７）（８）に適用して、重み付けを実施するものについて示したが、上記実施の形態２における式（９）（１０）に適用して、重み付けを実施するようにしてもよく、同様の効果を奏することができる。
【図面の簡単な説明】
【００５７】
【図１】この発明の実施の形態１による画像処理装置を示す構成図である。
【図２】定数ａが“１”である場合の変位検出処理部４の動作例を示す説明図である。
【図３】ｘdの値を−１、０、１とする場合の相関値演算処理部１２の動作例を示す説明図である。
【図４】この発明の実施の形態２による画像処理装置を示す構成図である。
【符号の説明】
【００５８】
１  画像入力部（画像信号入力手段）、２  水平データ積算処理部（画素値積算手段）、３  垂直データ積算処理部（画素値積算手段）、４，４ａ，４ｂ  変位検出処理部（第１の変位量検出手段、検出処理手段）、５，５ａ，５ｂ  曲率検出処理部（第１の曲率量検出手段、検出処理手段）、６，６ａ，６ｂ  変位検出処理部（第２の変位量検出手段、検出処理手段）、７，７ａ，７ｂ  曲率検出処理部（第２の曲率量検出手段、検出処理手段）、８，１１  特徴量バッファ（画素ずれ量検出手段）、９，１２  相関値演算処理部（画素ずれ量検出手段）、１０  画像出力部。
【図１】
2008079007.tif 000015
【図２】
2008079007.tif 000016
【図３】
2008079007.tif 000017
【図４】
2008079007.tif 000018