(57)【要約】【課題】検査画像列と参照画像列との間の比較に際し、画素値の誤差をある程度許容しながらオフセットの検出を行うこと。【解決手段】検査画像列及び参照画像列を取得する画像列取得部と、各参照フレームについてフレーム全体又はフレームを分割した各部に相当するブロックごとに参照画像変動量を算出する参照画像解析部と、前記参照画像変動量に基づいて前記参照画像列における画像の変動が特徴的に現れる参照フレーム番号、ブロックレベル、及びブロック番号を含むデータセットを選択する特徴位置選択部と、前記データセットに含まれるブロックレベル及びブロック番号における検査画像変動量が最も大きい検査フレーム番号を検出する検査画像解析部と、前記検査フレーム番号と前記参照フレーム番号との差に相当する前記検査画像列と前記参照画像列との間のオフセットを算出するオフセット算出部と、を備える画像処理装置を提供する。【選択図】図２
2010092234.tif 000002
【特許請求の範囲】
【請求項１】
検査の対象とする検査画像列及び当該検査画像列と比較される参照画像列を取得する画像列取得部と；
前記参照画像列に含まれる各参照フレームについて、フレーム全体又はフレームを分割した各部に相当するブロックごとに、前参照フレームに対する画像の変動の程度を表す参照画像変動量を算出する参照画像解析部と；
各参照フレームについて算出された前記ブロックごとの前記参照画像変動量に基づいて、前記参照画像列における画像の変動が特徴的に現れる参照フレーム番号、ブロックレベル、及びブロック番号を含むデータセットを選択する特徴位置選択部と；
前記検査画像列に含まれる各検査フレームについて、前記データセットに含まれる前記ブロックレベル及び前記ブロック番号における前検査フレームに対する検査画像変動量を計算し、当該検査画像変動量が最も大きい検査フレーム番号を検出する検査画像解析部と；
前記検査画像解析部により検出された前記検査フレーム番号と、当該検査フレーム番号の検出に用いられた前記データセットに含まれる前記参照フレーム番号とを用いて、前記検査画像列と前記参照画像列との間のオフセットを算出するオフセット算出部と；
を備える画像処理装置。
【請求項２】
前記特徴位置選択部は、前記ブロックレベルごとに前記参照画像変動量が大きい順に所定の数の前記データセットを選択する、請求項１に記載の画像処理装置。
【請求項３】
前記特徴位置選択部は、前記参照画像変動量が予め与えられる所定の閾値を上回っている複数の前記データセットが存在する場合には、当該複数の前記データセットを選択する、請求項１に記載の画像処理装置。
【請求項４】
前記特徴位置選択部は、さらに、各データセットに含まれる前記ブロックレベルと当該データセットに対応する前記参照画像変動量とに基づいて、各データセットの優先度を計算する、請求項２又は請求項３のいずれかに記載の画像処理装置。
【請求項５】
前記検査画像解析部は、前記特徴位置選択部により複数の前記データセットが選択された場合には、前記優先度の最も大きいデータセットについて前記検査画像変動量が最も大きい検査フレーム番号の検出を行う、請求項４に記載の画像処理装置。
【請求項６】
前記検査画像解析部は、さらに、検出した前記検査フレーム番号に対応する前記検査画像変動量と当該検査フレーム番号の検出に用いたデータセットに対応する前記参照画像変動量との差を比較する、請求項５に記載の画像処理装置。
【請求項７】
前記検査画像解析部は、前記検査画像変動量と前記参照画像変動量との差が所定の閾値よりも大きい場合には、さらに、前記優先度の次点のデータセットについて前記検査画像変動量が最も大きい検査フレーム番号の検出を行う、請求項６に記載の画像処理装置。
【請求項８】
前記優先度は、前記ブロックレベルに応じた１つの参照フレーム内のブロック数と前記参照画像変動量との積から所定の調整量を除いた値として計算される、請求項４に記載の画像処理装置。
【請求項９】
前記画像処理装置は、前記オフセット算出部により算出された前記オフセットに応じて、前記検査画像列又は前記参照画像列のフレーム位置を補正する画像列比較部、をさらに備える、請求項１に記載の画像処理装置。
【請求項１０】
前記画像列比較部は、さらに、フレーム位置を補正した後の前記検査画像列と前記参照画像列との間で画像の内容が一致するか否かを判定する、請求項９に記載の画像処理装置。
【請求項１１】
検査の対象とする検査画像列及び当該検査画像列と比較される参照画像列を取得するステップと；
前記参照画像列に含まれる各参照フレームについて、フレーム全体又はフレームを分割した各部に相当するブロックごとに、前参照フレームに対する画像の変動の程度を表す参照画像変動量を算出するステップと；
各参照フレームについて算出された前記ブロックごとの前記参照画像変動量に基づいて、前記参照画像列における画像の変動が特徴的に現れる参照フレーム番号、ブロックレベル、及びブロック番号を含むデータセットを選択するステップと；
前記検査画像列に含まれる各検査フレームについて、前記データセットに含まれる前記ブロックレベル及び前記ブロック番号における前検査フレームに対する検査画像変動量を計算し、当該検査画像変動量が最も大きい検査フレーム番号を検出するステップと；
検出された前記検査フレーム番号と、当該検査フレーム番号の検出に用いられた前記データセットに含まれる前記参照フレーム番号とを用いて、前記検査画像列と前記参照画像列との間のオフセットを算出するステップと；
を含む、画像検査方法。
【請求項１２】
画像処装置を制御するコンピュータを：
検査の対象とする検査画像列及び当該検査画像列と比較される参照画像列を取得する画像列取得部と；
前記参照画像列に含まれる各参照フレームについて、フレーム全体又はフレームを分割した各部に相当するブロックごとに、前参照フレームに対する画像の変動の程度を表す参照画像変動量を算出する参照画像解析部と；
各参照フレームについて算出された前記ブロックごとの前記参照画像変動量に基づいて、前記参照画像列における画像の変動が特徴的に現れる参照フレーム番号、ブロックレベル、及びブロック番号を含むデータセットを選択する特徴位置選択部と；
前記検査画像列に含まれる各検査フレームについて、前記データセットに含まれる前記ブロックレベル及び前記ブロック番号における前検査フレームに対する検査画像変動量を計算し、当該検査画像変動量が最も大きい検査フレーム番号を検出する検査画像解析部と；
前記検査画像解析部により検出された前記検査フレーム番号と、当該検査フレーム番号の検出に用いられた前記データセットに含まれる前記参照フレーム番号とを用いて、前記検査画像列と前記参照画像列との間のオフセットを算出するオフセット算出部と；
として機能させるための、プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像処理装置、画像検査方法、及びプログラムに関する。
【背景技術】
【０００２】
昨今のデジタル技術の発達により、デジタルカメラや携帯電話端末など様々な機器に動画を処理する画像処理機能が組み込まれている。かかる画像処理機能の開発に際しては、通常、動作保証やデグレードの防止などの観点から、開発された機能により出力された検査画像列と予め想定される参照画像列とを比較することにより、画像列に対して予期しない影響が生じていないかが検査される。
【０００３】
こうした検査画像列と参照画像列との比較は、従来、目視により行われることが多かった。これに対し、画像列の比較を自動化し正確に行うためには、検査画像列と参照画像列との比較の起点とするフレームが時間軸上で一致していることが求められる。即ち、検査画像列と参照画像列との間に時間軸上のずれが存在している場合には、時間軸上のずれ（以下、オフセットという）を検出して補正（位置合わせ）した後に、画像列の比較を行わなければならない。
【０００４】
検査画像列と参照画像列との間のオフセットを検出する方法としては、例えば、画像間の動きを検出する手法を応用することが考えられる。例えば、下記特許文献１及び下記特許文献２には、複数のフレーム間の画素値の差分を用いてフレーム間の動きを検出する手法が開示されている。
【０００５】
【特許文献１】特開平１０−１３６３７３号公報
【特許文献２】特開２００１−２５０２１号公報
【０００６】
即ち、このような動き検出の考え方を応用し、例えば、検査画像列と参照画像列のうち一方を時間軸方向に順次移動させながら、双方の画像列間の画素値の差分が最も少ない位置に位置合わせをすることも考えられる。
【発明の開示】
【発明が解決しようとする課題】
【０００７】
しかしながら、例えば過去にキャプチャされた参照画像列を使用する場合などでは、検査画像列と参照画像列にそれぞれ含まれるフレーム内の画素値は完全には一致せず、ある程度の範囲の誤差が存在する。そのため、フレーム間の画素値の差分を単純に用いると、画像の内容の特徴的な変動がその誤差に埋もれてしまい、オフセットを正確に検出して画像列を比較することは困難となる。
【０００８】
そこで、本発明は、上記問題に鑑みてなされたものであり、本発明の目的とするところは、検査画像列と参照画像列との間の比較に際し、画素値の誤差をある程度許容しながらオフセットの検出を行うことのできる、新規かつ改良された画像処理装置、画像検査方法、及びプログラムを提供することにある。
【課題を解決するための手段】
【０００９】
上記課題を解決するために、本発明のある観点によれば、検査の対象とする検査画像列及び当該検査画像列と比較される参照画像列を取得する画像列取得部と、前記参照画像列に含まれる各参照フレームについて、フレーム全体又はフレームを分割した各部に相当するブロックごとに、前参照フレームに対する画像の変動の程度を表す参照画像変動量を算出する参照画像解析部と、各参照フレームについて算出された前記ブロックごとの前記参照画像変動量に基づいて、前記参照画像列における画像の変動が特徴的に現れる参照フレーム番号、ブロックレベル、及びブロック番号を含むデータセットを選択する特徴位置選択部と、前記検査画像列に含まれる各検査フレームについて、前記データセットに含まれる前記ブロックレベル及び前記ブロック番号における前検査フレームに対する検査画像変動量を計算し、当該検査画像変動量が最も大きい検査フレーム番号を検出する検査画像解析部と、前記検査画像解析部により検出された前記検査フレーム番号と、当該検査フレーム番号の検出に用いられた前記データセットに含まれる前記参照フレーム番号とを用いて、前記検査画像列と前記参照画像列との間のオフセットを算出するオフセット算出部と、を備える画像処理装置が提供される。
【００１０】
また、前記特徴位置選択部は、前記ブロックレベルごとに前記参照画像変動量が大きい順に所定の数の前記データセットを選択してもよい。
【００１１】
また、前記特徴位置選択部は、前記参照画像変動量が予め与えられる所定の閾値を上回っている複数の前記データセットが存在する場合には、当該複数の前記データセットを選択してもよい。
【００１２】
また、前記特徴位置選択部は、さらに、各データセットに含まれる前記ブロックレベルと当該データセットに対応する前記参照画像変動量とに基づいて、各データセットの優先度を計算してもよい。
【００１３】
また、前記検査画像解析部は、前記特徴位置選択部により複数の前記データセットが選択された場合には、前記優先度の最も大きいデータセットについて前記検査画像変動量が最も大きい検査フレーム番号の検出を行ってもよい。
【００１４】
また、前記検査画像解析部は、さらに、検出した前記検査フレーム番号に対応する前記検査画像変動量と当該検査フレーム番号の検出に用いたデータセットに対応する前記参照画像変動量との差を比較してもよい。
【００１５】
また、前記検査画像解析部は、前記検査画像変動量と前記参照画像変動量との差が所定の閾値よりも大きい場合には、さらに、前記優先度の次点のデータセットについて前記検査画像変動量が最も大きい検査フレーム番号の検出を行ってもよい。
【００１６】
また、前記優先度は、前記ブロックレベルに応じた１つの参照フレーム内のブロック数と前記参照画像変動量との積から所定の調整量を除いた値として計算されてもよい。
【００１７】
また、前記画像処理装置は、前記オフセット算出部により算出された前記オフセットに応じて、前記検査画像列又は前記参照画像列のフレーム位置を補正する画像列比較部、をさらに備えてもよい。
【００１８】
また、前記画像列比較部は、さらに、フレーム位置を補正した後の前記検査画像列と前記参照画像列との間で画像の内容が一致するか否かを判定してもよい。
【００１９】
上記課題を解決するために、本発明の別の観点によれば、検査の対象とする検査画像列及び当該検査画像列と比較される参照画像列を取得するステップと、前記参照画像列に含まれる各参照フレームについて、フレーム全体又はフレームを分割した各部に相当するブロックごとに、前参照フレームに対する画像の変動の程度を表す参照画像変動量を算出するステップと、各参照フレームについて算出された前記ブロックごとの前記参照画像変動量に基づいて、前記参照画像列における画像の変動が特徴的に現れる参照フレーム番号、ブロックレベル、及びブロック番号を含むデータセットを選択するステップと、前記検査画像列に含まれる各検査フレームについて、前記データセットに含まれる前記ブロックレベル及び前記ブロック番号における前検査フレームに対する検査画像変動量を計算し、当該検査画像変動量が最も大きい検査フレーム番号を検出するステップと、検出された前記検査フレーム番号と、当該検査フレーム番号の検出に用いられた前記データセットに含まれる前記参照フレーム番号とを用いて、前記検査画像列と前記参照画像列との間のオフセットを算出するステップと、を含む画像検査方法が提供される。
【００２０】
上記課題を解決するために、本発明の別の観点によれば、画像処装置を制御するコンピュータを、検査の対象とする検査画像列及び当該検査画像列と比較される参照画像列を取得する画像列取得部と、前記参照画像列に含まれる各参照フレームについて、フレーム全体又はフレームを分割した各部に相当するブロックごとに、前参照フレームに対する画像の変動の程度を表す参照画像変動量を算出する参照画像解析部と、各参照フレームについて算出された前記ブロックごとの前記参照画像変動量に基づいて、前記参照画像列における画像の変動が特徴的に現れる参照フレーム番号、ブロックレベル、及びブロック番号を含むデータセットを選択する特徴位置選択部と、前記検査画像列に含まれる各検査フレームについて、前記データセットに含まれる前記ブロックレベル及び前記ブロック番号における前検査フレームに対する検査画像変動量を計算し、当該検査画像変動量が最も大きい検査フレーム番号を検出する検査画像解析部と、前記検査画像解析部により検出された前記検査フレーム番号と、当該検査フレーム番号の検出に用いられた前記データセットに含まれる前記参照フレーム番号とを用いて、前記検査画像列と前記参照画像列との間のオフセットを算出するオフセット算出部と、として機能させるためのプログラムが提供される。
【発明の効果】
【００２１】
以上説明したように、本発明に係る画像処理装置、画像検査方法、及びプログラムによれば、検査画像列と参照画像列との間の比較に際し、画素値の誤差をある程度許容しながらオフセットの検出を行うことができる。
【発明を実施するための最良の形態】
【００２２】
以下に添付図面を参照しながら、本発明の好適な実施の形態について詳細に説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。
【００２３】
また、以下の順序にしたがって当該「発明を実施するための最良の形態」を説明する。
１．検査画像列と参照画像列の概要
２．一実施形態に係る画像処理装置の構成
３．一実施形態に係る画像処理の流れ
４．まとめ
【００２４】
＜１．検査画像列と参照画像列の概要＞
まず、図１を用いて、本発明の一実施形態において使用される検査画像列と参照画像列の概要について説明する。図１には、本発明の一実施形態に係る画像処理装置１０のほか、被検査装置２０及び参照画像生成装置３０が示されている。
【００２５】
被検査装置２０には、検査の対象となる検査画像列をソース画像列から生成するソフトウェアが搭載される。被検査装置２０に搭載されるソフトウェアは、例えば、開発の目的とする画像処理機能を実現するためにアップデートされたソフトウェアであって、入力されるソース画像列に対し所定の画像処理を施した後、検査画像列を出力する。当該ソフトウェアは、例えば、被検査装置２０に設けられたＲＯＭ（Read  Only  Memory）などに格納される。なお、被検査装置２０に搭載されるソフトウェアにより実現される画像処理機能は、例えば、符号化や復号化、画像の加工、高画質化など、任意の機能であってよい。
【００２６】
参照画像生成装置３０は、被検査装置２０に搭載されるソフトウェアの動作保証などのために検査画像列と比較される参照画像列を、ソース画像列から生成する。参照画像生成装置３０による参照画像列の生成処理は、典型的には、既に動作保証のための検査を通過した画像処理機能によって行われる。例えば、参照画像生成装置３０による参照画像列の生成処理は、被検査装置２０においてアップデートされる前の過去のモデル（旧バージョンのソフトウェア）に基づく画像処理であってもよい。
【００２７】
被検査装置２０及び参照画像生成装置３０に入力されるソース画像列は、被検査装置２０に搭載されるソフトウェアの動作を保証するために十分なフレーム数又は内容を有する任意の画像列であってよい。
【００２８】
画像処理装置１０は、被検査装置２０から出力される検査画像列と参照画像生成装置３０から出力される参照画像列とを比較し、例えば、画像列の内容の一致又は不一致を判定してその判定結果を出力する。それにより、例えば、被検査装置２０に搭載されるソフトウェアによりデグレードが生じていないかなどを検査することができる。
【００２９】
ここで、検査画像列と参照画像列との間には、例えば、被検査装置２０又は参照画像生成装置３０から出力される画像列のキャプチャタイミングの違いなどを原因として、時間軸上のずれ、即ちオフセットが生じている場合がある。検査画像列と参照画像列との間にオフセットが生じている場合、画像処理装置１０は、例えば、オフセットの程度（例えば、フレーム数何枚分など）を算出して位置合わせをした上で、画像列の比較を行う必要がある。そこで、本実施形態において、画像処理装置１０は、次節より詳しく述べる方法により、検査画像列と参照画像列との間のオフセットを算出する。
【００３０】
＜２．一実施形態に係る画像処理装置の構成＞
次に、図２〜図６を用いて、一実施形態に係る画像処理装置１０の構成についてより詳しく説明する。
【００３１】
図２は、画像処理装置１０の論理的な構成を示すブロック図である。図２を参照すると、画像処理装置１０は、画像列取得部１１０、参照画像解析部１２０、特徴位置選択部１３０、検査画像解析部１４０、オフセット算出部１５０、及び画像列比較部１６０を備える。
【００３２】
［画像列取得部］
画像列取得部１１０は、例えば、図１に示した被検査装置２０から出力される検査画像列をキャプチャすることにより取得し、検査画像解析部１４０及び画像列比較部１６０へ出力する。また、画像列取得部１１０は、例えば、図１に示した参照画像生成装置３０から出力される参照画像列をキャプチャすることにより取得し、参照画像解析部１２０及び画像列比較部１６０へ出力する。その代わりに、画像列取得部１１０は、例えば、事前にキャプチャされた検査画像列又は参照画像列を画像処理装置１０の内部又は外部の記憶装置から読み込んで取得し、前述の各処理部へ出力してもよい。
【００３３】
［参照画像解析部］
参照画像解析部１２０は、画像列取得部１１０から入力された参照画像列に含まれる各参照フレームについて、フレーム全体又はフレームを分割した各部に相当するブロックごとに、参照画像変動量Ｅｒ（Ｎｒ，Ｌｖ，ｂ）を算出する。ここで、Ｎｒは参照フレーム番号、Ｌｖはブロックの分割の程度を表すブロックレベル、ｂは分割された各ブロックに付されるブロック番号である。
【００３４】
図３は、本実施形態で扱うブロックの例を示す説明図である。図３には、ブロックレベルＬｖ＝０、１、２の３通りのブロックの集合が示されており、各ブロックにはブロック番号ｂが付されている。
【００３５】
このうち、図３（Ａ）は、ブロックレベルＬｖ＝０として、フレーム全体を１つのブロックとした例を示している。この場合の１フレーム内のブロック数は１個である。また、図３（Ｂ）は、ブロックレベルＬｖ＝１として、フレーム全体を縦横２×２個のブロックに分割した例を示している。この場合の１フレーム内のブロック総数は４個である。また、図３（Ｃ）は、ブロックレベルＬｖ＝２として、フレーム全体を縦横４×４個のブロックに分割した例を示している。この場合の１フレーム内のブロックの数は１６個である。なお、ブロックの分割方法やブロックレベルの定義はかかる例に限定されない。例えば、縦横３×３個のブロックをブロックレベルＬｖ＝２として定義してもよい。
【００３６】
参照画像変動量Ｅｒ（Ｎｒ，Ｌｖ，ｂ）は、参照フレーム番号Ｎｒ内のブロックレベルＬｖ、ブロック番号ｂにより特定されるブロックにおける、前の参照フレームに対する画像の変動の程度を表す指標である。参照画像変動量Ｅｒ（Ｎｒ，Ｌｖ，ｂ）は、例えば、ブロックレベルＬｖ及びブロック番号ｂにより特定されるブロックに含まれる画素位置（ｘ，ｙ）を用いて、式（１）により定義され得る。
【００３７】
【数１】
2010092234.tif 000003
【００３８】
ここで、ｅ（Ｎｒ，ｘ，ｙ）は画素単位の画像変動量を表す。例えば、画素位置（ｘ，ｙ）におけるＮｒ番目のフレームの画素値をｐ（Ｎｒ，ｘ，ｙ）とすると、画素単位での画像変動量ｅ（Ｎｒ，ｘ，ｙ）は、式（２）のように定義され得る。
【００３９】
【数２】
2010092234.tif 000004
【００４０】
ここで、ＣｍｐＦ（ｐ１，ｐ２）は、２つの画素ｐ１及びｐ２の画素値の比較を行う比較関数である。ＣｍｐＦ（ｐ１，ｐ２）は、例えば、式（３）に示すように、画素ｐ１及びｐ２のコンポーネントごとの差分絶対値をベクトル値として出力する関数であってもよい。
【００４１】
【数３】
2010092234.tif 000005
【００４２】
また、ＣｍｐＦ（ｐ１，ｐ２）は、例えば、式（４）に示すように、画素ｐ１及びｐ２のコンポーネントごとの誤差割合をベクトル値として出力する関数であってもよい。
【００４３】
【数４】
2010092234.tif 000006
【００４４】
なお、画素のコンポーネント（ｃ０，ｃ１，ｃ２）とは、例えばＲＧＢなどの三次元色空間において画素値を構成する個々の座標値を意味する。
【００４５】
また、ＣｍｐＦ（ｐ１，ｐ２）は、例えば、式（５）に示すように、画素ｐ１及びｐ２の類似度をスカラー値として出力する関数であってもよい。なお、式（５）において、ｐ１・ｐ２は、ベクトル値ｐ１及びｐ２の内積を表す。
【００４６】
【数５】
2010092234.tif 000007
【００４７】
さらに、ＣｍｐＦ（ｐ１，ｐ２）は、例えば、式（６）に示すように、画素ｐ１及びｐ２のコンポーネントごとの差分絶対値のうちの最大値を算出する関数であってもよい。
【００４８】
【数６】
2010092234.tif 000008
【００４９】
その代わりに、ＣｍｐＦ（ｐ１，ｐ２）は、式（７）〜（９）のいずれかに示すように、画素ｐ１及びｐ２のコンポーネントごとの差分絶対値の平均値や自乗和の算出、又は許容誤差の判定を行う関数であってもよい。
【００５０】
【数７】
2010092234.tif 000009
【００５１】
なお、式（９）におけるＴ１、Ｔ２、Ｔ３は、各コンポーネントについての許容誤差を表す閾値である。この場合、差分絶対値が閾値を超えるコンポーネントが１つでもあれば、比較関数ＣｍｐＦは許容誤差の判定結果として真値“１”を出力し、そうでなければ比較関数ＣｍｐＦは“０”を出力する。
【００５２】
また、式（６）〜（９）において、画素ｐ１及びｐ２のコンポーネントごとの差分絶対値を用いる代わりに、式（４）に示したような画素ｐ１及びｐ２のコンポーネントごとの
誤差割合を用いてもよい。
【００５３】
なお、比較関数ＣｍｐＦがベクトル値を出力する場合は、参照画像変動量Ｅｒもまたベクトル値であってもよい。その場合、後述する特徴位置選択部１３０における参照画像変動量Ｅｒの大きさの判定は、例えばベクトルの長さを基準として行われてもよく、又はベクトルの個々のコンポーネントごとに行われてもよい。
【００５４】
図４は、参照画像解析部１２０により参照画像変動量Ｅｒ（Ｎｒ，Ｌｖ，ｂ）が算出される様子を一例として描いた説明図である。
【００５５】
図４を参照すると、参照フレーム番号Ｎｒ−１からＮｒにかけて、フレーム内の下部中央左側に２名の人物が登場している。また、フレーム内のそれ以外の部分には変動は見られない。この場合、ブロックレベルＬｖ＝０のブロック番号ｂ＝１、ブロックレベルＬｖ＝１のブロック番号ｂ＝３、並びにブロックレベルＬｖ＝２のブロック番号ｂ＝１０及び１４において、それぞれ所定の参照画像変動量Ｅｒ（Ｎｒ，Ｌｖ，ｂ）＞０が算出される。それ以外のブロックでは、参照画像変動量Ｅｒ（Ｎｒ，Ｌｖ，ｂ）＝０となる。
【００５６】
参照画像解析部１２０は、参照画像列に含まれる各参照フレームについて、このように複数のブロックレベルＬｖにわたってブロック番号ｂごとに参照画像変動量Ｅｒ（Ｎｒ，Ｌｖ，ｂ）を順次計算し、計算結果を特徴位置選択部１３０へ出力する。
【００５７】
［特徴位置選択部］
特徴位置選択部１３０は、参照画像解析部１２０から入力された前述の参照画像変動量Ｅｒ（Ｎｒ，Ｌｖ，ｂ）に基づいて、参照画像列における画像の変動が特徴的に現れる参照フレーム番号、ブロックレベル、及びブロック番号を含むデータセットを選択する。
【００５８】
図５は、特徴位置選択部１３０による前述のデータセットの選択処理の一例を説明するための説明図である。
【００５９】
図５を参照すると、４枚の参照フレーム（Ｎｒ＝１〜４）について、３段階のブロックレベル（Ｌｖ＝０，１，２）ごとに参照画像変動量Ｅｒ（Ｎｒ，Ｌｖ，ｂ）が網掛けの濃淡により示されている。特徴位置選択部１３０は、例えば、このような入力データの中から、ブロックレベルＬｖごとに、参照画像変動量Ｅｒが大きい順に上限数Ｍ（Ｌｖ）個の参照フレーム番号、ブロックレベル、及びブロック番号を含むデータセットを選択する。
【００６０】
図５の例では、ブロックレベルＬｖ＝０の場合の上限数Ｍ（０）＝１である。そして、ブロックレベルＬｖ＝０については、参照フレーム番号Ｎｒ＝３、ブロック番号ｂ＝１の参照画像変動量Ｅｒが最も大きい。よって、選択されるデータセット（Ｎｒ，Ｌｖ，ｂ）＝（３，０，１）となる。なお、選択されるデータセットには図５において矢印記号を付している。
【００６１】
また、ブロックレベルＬｖ＝１の場合の上限数Ｍ（１）＝２である。図５に濃淡で示された参照画像変動量Ｅｒに基づき、ブロックレベルＬｖ＝１について選択されるデータセット（Ｎｒ，Ｌｖ，ｂ）＝（３，１，２）、（４，１，２）となる。
【００６２】
また、ブロックレベルＬｖ＝２の場合の上限数Ｍ（２）＝４である。図５に濃淡で示された参照画像変動量Ｅｒに基づき、ブロックレベルＬｖ＝２について選択されるデータセット（Ｎｒ，Ｌｖ，ｂ）＝（１，２，１４）、（３，２，３）、（３，２，４）、（４，２，３）となる。
【００６３】
即ち、図５の例の場合、特徴位置選択部１３０は、（Ｎｒ，Ｌｖ，ｂ）＝（３，０，１）、（３，１，２）、（４，１，２）、（１，２，１４）、（３，２，３）、（３，２，４）、（４，２，３）の計７個のデータセットを選択する。
【００６４】
なお、特徴位置選択部１３０は、ブロックレベルＬｖごとに予め与えられる所定の閾値を参照画像変動量Ｅｒが上回っているか否かを各参照フレームのブロックごとに判定し、当該閾値を上回っているデータセットのみを選択してもよい。また、次に述べるデータセットごとの優先度を計算した後に、予め与えられる所定の閾値を上回る優先度を有するデータセットのみを選択してもよい。
【００６５】
そして、特徴位置選択部１３０は、図５を用いて例示した手法により選択されたデータセット（Ｎｒ、Ｌｖ、ｂ）及び対応する参照画像変動量Ｅｒに基づいて、検査画像列の検査に使用する順序を決定するためのデータセットごとの優先度を計算する。
【００６６】
例えば、データセットごとの優先度Ｐは、次式により定義され得る。
【００６７】
【数８】
2010092234.tif 000010
【００６８】
ここで、Ｂ（Ｌｖ）は、ブロックレベルＬｖに応じた１フレーム内のブロック数である。例えば、図３に示したブロックレベルを用いる場合には、Ｂ（０）＝１、Ｂ（１）＝４、Ｂ（２）＝１６となる。このように参照画像変動量Ｅｒにブロック総数を乗算することにより、フレーム全体の緩やかな画像の変動に１ブロック内の急激な画像の変動が埋もれることなく、参照画像列内の画像の変動の特徴位置が識別される。
【００６９】
また、Ｍｇｎ（Ｌｖ）は、ブロックレベルＬｖに応じて予め与えられる調整量である。Ｍｇｎ（Ｌｖ）の値は、典型的には、チューニングによって決定される。式（１０）に調整量Ｍｇｎ（Ｌｖ）を導入することで、フレーム全体の急激な画像の変動を１ブロック内の同等の画像の変動よりも優先して扱うことが可能となる。
【００７０】
特徴位置選択部１３０は、このようにデータセットごとの優先度Ｐを計算し、選択したデータセット（Ｎｒ，Ｌｖ，ｂ）と対応する参照画像変動量Ｅｒと共に検査画像解析部１４０へ出力する。
【００７１】
［検査画像解析部］
検査画像解析部１４０は、画像列取得部１１０から入力される検査画像列において画像の変動が特徴的に現れる検査フレーム番号を検出する。
【００７２】
より具体的には、検査画像解析部１４０は、まず、最も優先度の大きい前述のデータセットに含まれるブロックレベルＬｖ及びブロック番号ｂを所与の値とし、各検査フレームについて前検査フレームに対する検査画像変動量Ｅｔ（Ｎｒ，Ｌｖ，ｂ）を順次計算する。そして、検査画像解析部１４０は、計算された検査画像変動量Ｅｔ（Ｎｒ，Ｌｖ，ｂ）が最も大きい検査フレーム番号Ｎｔを、検査画像列において画像の変動が特徴的に現れる検査フレーム番号Ｎｔｍａｘとして検出する。なお、検査画像変動量Ｅｔ（Ｎｒ，Ｌｖ，ｂ）は、式（１）〜（９）を用いて説明した参照画像変動量Ｅｒ（Ｎｒ，Ｌｖ，ｂ）の算出式と同様にして計算され得る。
【００７３】
その結果、検査フレーム番号Ｎｔｍａｘの検出に使用されたデータセットに含まれる参照フレーム番号ＮｒをＮｒｍａｘとすると、検査フレーム番号Ｎｔｍａｘと参照フレーム番号Ｎｒｍａｘとの差を、画像列間のオフセットとして扱うことができる。
【００７４】
ここで、検査画像解析部１４０は、さらに、検出した検査フレーム番号Ｎｔｍａｘに対応する検査画像変動量Ｅｔとその検出に用いたデータセットに対応する参照画像変動量Ｅｒとを比較することにより、検出結果の妥当性を検証してもよい。
【００７５】
例えば、検査画像解析部１４０は、まず、検査フレーム番号Ｎｔｍａｘに対応する検査画像変動量Ｅｔと対応する参照画像変動量Ｅｒとの差を計算する。そして、検査画像解析部１４０は、計算した差が予め定義された所定の閾値よりも大きい場合には、検出した検査フレーム番号Ｎｔｍａｘが参照画像列における特徴位置と正しく対応していないと判定する。その場合、検査画像解析部１４０は、例えば、特徴位置選択部１３０から入力されたデータセットのうち優先度が次点のデータセットについてあらためて検査画像変動量Ｅｔが最も大きい検査フレーム番号の検出を行う。検査画像解析部１４０は、このような処理を優先度の大きいデータセットから順に行い、検出された検査フレーム番号が妥当であると検証された時点で、その検査フレーム番号Ｎｔｍａｘ及び対応する参照フレーム番号Ｎｒｍａｘをオフセット算出部１５０へ出力する。
【００７６】
このような検出結果の妥当性の検証により、画像処理装置１０によるオフセットの算出の信頼性が高められる。
【００７７】
［オフセット算出部］
オフセット算出部１５０は、前述した検査フレーム番号Ｎｔｍａｘ、及び当該検査フレーム番号Ｎｔｍａｘに対応する参照フレーム番号Ｎｒｍａｘを用いて、検査画像列と参照画像列との間のオフセットを算出する。
【００７８】
検査画像列と参照画像列との間のオフセットは、例えば、次式で与えられる。
【００７９】
【数９】
2010092234.tif 000011
【００８０】
そして、オフセット算出部１５０は、このように算出したオフセットを画像列比較部１６０へ出力する。
【００８１】
［画像列比較部］
画像列比較部１６０は、オフセット算出部１５０から入力されるオフセットに応じて、画像列取得部１１０から入力される検査画像列又は参照画像列のフレーム位置を補正する。そして、画像列比較部１６０は、フレーム位置を補正した後の検査画像列と参照画像列との間で画像の内容が一致するか否かを判定する。
【００８２】
より具体的には、画像列比較部１６０は、例えばオフセットが正の値であれば、検査画像列の検査フレーム番号からオフセットを減算した後、検査画像列と参照画像列の１番目のフレームから順に画像の内容を比較する。即ち、この場合には、検査画像列の最初のオフセット分の検査フレームはスキップされる。また、画像列比較部１６０は、例えばオフセットが負の値であれば、参照画像列の参照フレーム番号からオフセットを減算した後、検査画像列と参照画像列の１番目のフレームから順に画像の内容を比較する。即ち、この場合には、参照画像列の最初のオフセット分の参照フレームはスキップされる。
【００８３】
画像列比較部１６０による検査画像列と参照画像列との間の比較は、例えば、式（１）を応用し、フレーム全体に含まれる画素位置（Ｘ，Ｙ）についての総画像変動量Ｅａを求めることにより行われる。フレーム位置の補正後のＮ番目の検査フレーム及び参照フレームについての検査画像列と参照画像列との間の総画像変動量Ｅａ（Ｎ）は、次式で与えられる。
【００８４】
【数１０】
2010092234.tif 000012
【００８５】
ここで、ｅ´（Ｎ，Ｘ，Ｙ）は、検査画像列と参照画像列との間の画素単位の画像変動量を表す。例えば、画素位置（Ｘ，Ｙ）におけるＮ番目の検査フレームの画素値をｐｔ（Ｎ，Ｘ，Ｙ）、参照フレームの画素値をｐｒ（Ｎ，Ｘ，Ｙ）とすると、画素単位での画像変動量ｅ´（Ｎ，Ｘ，Ｙ）は、式（１３）のように定義される。
【００８６】
【数１１】
2010092234.tif 000013
【００８７】
そして、画像列比較部１６０は、例えば、式（１２）及び（１３）に従って計算した総画像変動量Ｅａ（Ｎ）が予め与えられた閾値を超える場合には、当該フレーム番号Ｎにおいて検査画像列と参照画像列とは一致しないと判定する。なお、Ｅａ（Ｎ）がベクトル値となる場合には、閾値はベクトルの要素の数だけ与えられる。
【００８８】
ここで、画像列比較部１６０は、画像の内容が一致しないフレームが存在した場合には、一致しない箇所を特定可能な任意の情報を出力する。例えば、画像列比較部１６０は、不一致となった検査フレーム番号若しくは参照フレーム番号、総画像変動量Ｅａ、又は不一致の原因となった画素位置や検査フレーム若しくは参照フレームの画素値などを出力することができる。
【００８９】
ここまで、図２〜図５を用いて、本実施形態に係る画像処理装置１０の構成について説明した。次に、本実施形態に係る画像処理装置１０による画像処理の流れについて、フローチャートを用いて説明する。
【００９０】
＜３．一実施形態に係る画像処理の流れ＞
図６は、画像処理装置１０による処理の全体的な流れを示すフローチャートである。
【００９１】
図６において、まず、画像列取得部１１０は、検査の対象とする検査画像列及び当該検査画像列と比較される参照画像列を取得する（Ｓ２０２）。そして、画像列取得部１１０は、取得した参照画像列を参照画像解析部１２０及び画像列比較部１６０へ、検査画像列を検査画像解析部１４０及び画像列比較部１６０へ出力する。
【００９２】
次に、参照画像解析部１２０は、参照画像列に含まれる各参照フレームについて、ブロックレベルＬｖ及びブロック番号ｂで特定されるブロックごとに、前参照フレームに対する画像の変動の程度を表す参照画像変動量Ｅｒ（Ｎｒ，Ｌｖ，ｂ）を算出する（Ｓ２０４）。なお、ここでの処理の流れは、図７を用いて後にさらに説明する。
【００９３】
その後、特徴位置選択部１３０は、算出された前述の参照画像変動量Ｅｒに基づいて、参照画像列における画像の変動が特徴的に現れる参照フレーム番号Ｎｒ、ブロックレベルＬｖ、及びブロック番号ｂを含むデータセットを選択する（Ｓ２０６）。なお、ここでの処理の流れは、図８を用いて後にさらに説明する。
【００９４】
そして、検査画像解析部１４０は、各検査フレームについて前述のデータセットに含まれるブロックレベルＬｖ及びブロック番号ｂにおける検査画像変動量Ｅｔを計算し、当該検査画像変動量Ｅｔが最も大きい検査フレーム番号Ｎｔｍａｘを検出する（Ｓ２０８）。なお、ここでの処理の流れは、図９を用いて後にさらに説明する。
【００９５】
さらに、オフセット算出部１５０は、検査画像解析部１４０により検出された検査フレーム番号Ｎｔｍａｘと対応する参照フレーム番号Ｎｒｍａｘとの差に相当する、検査画像列と参照画像列との間のオフセットを算出する（Ｓ２１０）。
【００９６】
そして、画像列比較部１６０は、オフセット算出部１５０により算出されたオフセットに応じて検査画像列又は参照画像列のフレーム位置を補正した後、当該検査画像列と参照画像列との間で画像の内容が一致するか否かを判定する（Ｓ２１２）。なお、ここでの処理の流れは、図１０を用いて後にさらに説明する。
【００９７】
図７は、参照画像解析部１２０による参照画像解析処理のより具体的な流れの一例を示すフローチャートである。
【００９８】
図７において、まず、参照画像解析部１２０は、画像列取得部１１０から入力された参照画像列に含まれる参照フレームに対し、注目する参照フレームを参照フレーム番号の順に順次設定する（Ｓ３０２）。
【００９９】
そして、参照画像解析部１２０は、注目した参照フレームについて、複数のブロックレベルＬｖにわたってブロック番号ｂごとに参照画像変動量Ｅｒを算出する（Ｓ３０４）。
【０１００】
その後、参照画像解析部１２０は、全ての参照フレームについて参照画像変動量Ｅｒの算出が終了したか否かを判定する（Ｓ３０６）。ここで、参照画像変動量Ｅｒの算出が終了していない参照フレームが残っていれば、処理はＳ３０２へ戻り、次の参照フレームについて前述した処理が繰り返される。一方、全ての参照フレームについて参照画像変動量Ｅｒの算出が終了していれば、参照画像解析部１２０による参照画像解析処理は終了する。
【０１０１】
図８は、特徴位置選択部１３０による特徴位置選択処理のより具体的な流れの一例を示すフローチャートである。
【０１０２】
図８において、まず、特徴位置選択部１３０は、例えば、ブロックレベルＬｖごとに参照画像変動量Ｅｒの大きい順に所定の数のデータセット（Ｎｒ，Ｌｖ，ｂ）を選択する（Ｓ４０２）。その代わりに、特徴位置選択部１３０は、参照画像変動量Ｅｒが予め与えられる所定の閾値を上回っているデータセット（Ｎｒ，Ｌｖ，ｂ）をここで選択してもよい。
【０１０３】
次に、特徴位置選択部１３０は、各データセットに含まれるブロックレベルＬｖと対応する参照画像変動量Ｅｒとに基づいて、各データセットの優先度Ｐを計算する（Ｓ４０４）。このとき、特徴位置選択部１３０は、優先度Ｐを、ブロックレベルＬｖに応じた参照フレーム内のブロック数と参照画像変動量Ｅｒとの積から所定の調整量Ｍｇｎ（Ｌｖ）を除いた値として計算することができる。
【０１０４】
そして、特徴位置選択部１３０は、選択したデータセット（Ｎｒ，Ｌｖ，ｂ）、参照画像変動量Ｅｒ、及び計算した優先度Ｐを検査画像解析部１４０へ出力し（Ｓ４０６）、特徴位置選択部１３０による特徴位置選択処理は終了する。
【０１０５】
図９は、検査画像解析部１４０による検査画像解析処理のより具体的な流れの一例を示すフローチャートである。
【０１０６】
図９において、まず、検査画像解析部１４０は、特徴位置選択部１３０から入力されたデータセットのうち、例えば優先度Ｐの最も大きいデータセットを取得する（Ｓ５０２）。
【０１０７】
次に、検査画像解析部１４０は、取得したデータセットに含まれるブロックレベルＬｖ及びブロック番号ｂにおいて検査画像列内で検査画像変動量Ｅｔが最も大きくなる検査フレーム番号Ｎｔｍａｘを検出する（Ｓ５０４）。
【０１０８】
そして、検査画像解析部１４０は、検出した検査フレーム番号Ｎｔｍａｘに対応する検査画像変動量Ｅｔと検出に用いたデータセットに対応する参照画像変動量Ｅｒとの差が許容範囲内であるかを、所定の閾値を用いて判定する（Ｓ５０６）。
【０１０９】
ここで、検査画像変動量Ｅｔと参照画像変動量Ｅｒとの差が許容範囲内でなければ、検査画像解析部１４０は、優先度Ｐが次点のデータセットを取得し（Ｓ５０８）、検査画像変動量Ｅｔが最も大きい検査フレーム番号Ｎｔｍａｘの検出をあらためて行う。
【０１１０】
一方、検査画像変動量Ｅｔと参照画像変動量Ｅｒとの差が許容範囲内であれば、検査画像解析部１４０は、検出した検査フレーム番号Ｎｔｍａｘと対応する参照フレーム番号Ｎｒｍａｘをオフセット算出部１５０へ出力する（Ｓ５１０）。そして、検査画像解析部１４０による検査画像解析処理は終了する。
【０１１１】
図１０は、画像列比較部１６０による画像列比較処理のより具体的な流れの一例を示すフローチャートである。
【０１１２】
図１０において、まず、画像列比較部１６０は、オフセット算出部１５０により算出されたオフセットに応じて、検査画像列又は参照画像列のフレーム番号をシフトさせる（Ｓ６０２）。
【０１１３】
次に、画像列比較部１６０は、フレーム番号をシフトさせた後の検査画像列の検査フレームと参照画像列の参照フレームとを順次比較する（Ｓ６０４）。そして、両フレームが一致しているか否かが判定される（Ｓ６０６）。フレームの一致の判定は、例えば、前述した式（１２）及び式（１３）に従って計算した総画像変動量Ｅａが予め与えられる閾値を超えるか否かを判定することなどにより行うことができる。
【０１１４】
Ｓ６０６において、検査フレームと参照フレームとが一致していない場合には、一致しない箇所を特定可能な情報として、例えば、検査フレーム番号、参照フレーム番号、総画像変動量Ｅａ、又は不一致の原因となった画素位置や画素値などが出力される（Ｓ６０８）。
【０１１５】
その後、全ての検査フレームと参照フレームとの比較が終了したか否かが判定される（Ｓ６１０）。ここで、比較が終了していないフレームが残っていれば、処理はＳ６０４へ戻り、次のフレームの比較が行われる。一方、全てのフレームの比較が終了していれば、画像列比較部１６０による画像列比較処理は終了する。
【０１１６】
なお、本節で説明した一連の処理又はその一部をハードウェアで実現するかソフトウェアで実現するかは問わない。いずれかの処理をソフトウェアで実行させる場合には、ソフトウェアを構成するプログラムが、例えば図１１に示した汎用コンピュータなどを用いて実行される。
【０１１７】
図１１において、ＣＰＵ（Central  Processing  Unit）９０２は、汎用コンピュータの動作全般を制御する。ＲＯＭ９０４には、一連の処理の一部又は全部を記述したプログラム又はデータが格納される。ＲＡＭ（Random  Access  Memory）９０６には、処理の実行時にＣＰＵ９０２により用いられるプログラムやデータなどが一時的に記憶される。
【０１１８】
ＣＰＵ９０２、ＲＯＭ９０４、及びＲＡＭ９０６は、バス９１０を介して相互に接続される。バス９１０にはさらに、入出力インタフェース９１２が接続される。
【０１１９】
入出力インタフェース９１２は、ＣＰＵ９０２、ＲＯＭ９０４、及びＲＡＭ９０６と、入力装置９２０、出力装置９２２、記憶装置９２４、通信装置９２６、及びドライブ９３０とを接続するためのインタフェースである。
【０１２０】
入力装置９２０は、例えばマウスやキーボードなどの入力装置を介して、ユーザからの指示や情報入力を受け付ける。出力装置９２２は、例えばＣＲＴ（Cathode  Ray  Tube）、液晶ディスプレイ、ＯＬＥＤ（Organic  Light  Emitting  Diode）などの表示装置、又はスピーカなどの音声出力装置を介してユーザに判定結果などの情報を出力する。
【０１２１】
記憶装置９２４は、例えばハードディスクドライブ又はフラッシュメモリなどにより構成され、プログラムやプログラムデータ、画像データなどを記憶する。通信装置９２６は、ＬＡＮ又はインターネットなどのネットワークを介する通信処理を行う。ドライブ９３０は、必要に応じて汎用コンピュータに設けられ、例えばドライブ９３０にはリムーバブルメディア９３２が装着される。
【０１２２】
＜４．まとめ＞
ここまで、図１〜図１１を用いて、本発明の一実施形態に係る画像処理装置１０の構成と処理の流れについて説明を行った。本実施形態によれば、検査画像列と参照画像列との間にある程度の誤差が存在する場合にも、検査画像列と参照画像列との間のオフセットを算出して時間方向の位置合わせを自動的に行うことが可能となる。それにより、位置合わせ後に自動的に検査画像列と参照画像列とを比較して画像変動量を算出し、画像変動量が所定の範囲内かに応じて検査画像列が合格基準を満たすか否かを判定することが可能となる。
【０１２３】
そして、こうした合格判定を用いることで、画像処理を行うソフトウェアの開発中に発生した副作用的な不具合の検知を自動的に行う回帰試験環境を整えることが可能となり、製品開発の評価効率の向上と品質の向上が期待される。
【０１２４】
また、本実施形態によれば、検査画像列と参照画像列との間のオフセットの算出の基準となる画像の変動の特徴位置は、フレーム全体だけでなくフレームの一部分であるブロックの単位で選択され得る。それにより、例えばフレーム全体のシーンチェンジが存在しない画像列であっても、ブロックの単位で特徴的な画像の変動が存在すれば、その画像の変動が誤差に埋もれることなく検出される。
【０１２５】
以上、添付図面を参照しながら本発明の好適な実施形態について説明したが、本発明は係る例に限定されないことは言うまでもない。当業者であれば、特許請求の範囲に記載された範疇内において、各種の変更例又は修正例に想到し得ることは明らかであり、それらについても当然に本発明の技術的範囲に属するものと了解される。
【０１２６】
例えば、本実施形態に係る一連の処理又はその一部を、必ずしも図６〜図１０のフローチャートに記載された順序に沿って実行しなくてもよい。各処理ステップは、並列的あるいは個別に独立して実行される処理を含んでもよい。
【図面の簡単な説明】
【０１２７】
【図１】一実施形態に係る検査画像列と参照画像列の概要を示す模式図である。
【図２】一実施形態に係る画像処理装置の論理的な構成を示すブロック図である。
【図３】一実施形態に係るブロックレベルとブロック番号について説明するための説明図である。
【図４】一実施形態に係る参照画像解析部の処理を説明するための説明図である。
【図５】一実施形態に係る特徴位置選択部の処理を説明するための説明図である。
【図６】一実施形態に係る画像処理装置による処理の流れを示すフローチャートである。
【図７】一実施形態に係る参照画像解析処理の流れを示すフローチャートである。
【図８】一実施形態に係る特徴位置選択処理の流れを示すフローチャートである。
【図９】一実施形態に係る検査画像解析処理の流れを示すフローチャートである。
【図１０】一実施形態に係る画像列比較処理の流れを示すフローチャートである。
【図１１】汎用コンピュータの構成例を示すブロック図である。
【符号の説明】
【０１２８】
１０          画像処理装置
１１０        画像列取得部
１２０        参照画像解析部
１３０        特徴位置選択部
１４０        検査画像解析部
１５０        オフセット算出部
１６０        画像列比較部
【図１】
2010092234.tif 000014
【図２】
2010092234.tif 000015
【図３】
2010092234.tif 000016
【図４】
2010092234.tif 000017
【図５】
2010092234.tif 000018
【図６】
2010092234.tif 000019
【図７】
2010092234.tif 000020
【図８】
2010092234.tif 000021
【図９】
2010092234.tif 000022
【図１０】
2010092234.tif 000023
【図１１】
2010092234.tif 000024