(57)【要約】【課題】  複数の撮像素子を用いて取得された被写体像の複数の画像を、高速かつ高精度に接続する画像処理の技術を提供することにある。【解決手段】  画像処理方法は、同一平面内に相互に間隔をおいて複数の撮像素子を配列した撮像装置を用いて、同一平面に結像する被写体像を撮像して第１の画像群を生成する工程と、各々の撮像素子の撮像範囲が一部重複するように、被写体像と撮像装置との相対位置をシフトさせる工程と、シフト後に撮像装置で被写体像を撮像して第２の画像群を生成する工程と、第１の画像群および第２の画像群における画像の重複領域のうち、同一平面内で相互の距離が最も離れた重複領域の一対を選択し、選択された重複領域の一対を解析して第１の画像群に対する第２の画像群のシフト量を求める工程と、シフト量に基づいて各々の画像間の位置合わせを行い、第１の画像群と第２の画像群とを接続する工程とを備えることを特徴とする。【選択図】    図１
2009031909.tif 000002
【特許請求の範囲】
【請求項１】
同一平面内に相互に間隔を空けて複数の撮像素子を配列した撮像装置を用いて、前記平面に結像する被写体像を撮像して第１の画像群を生成する工程と、
前記各々の撮像素子の撮像範囲が一部重複するように、前記被写体像と前記撮像装置との相対位置をシフトさせる工程と、
前記シフト後に前記撮像装置で前記被写体像を撮像して第２の画像群を生成する工程と、
前記第１の画像群および前記第２の画像群における画像の重複領域のうち、前記平面内で相互の距離が最も離れた前記重複領域の一対を選択し、選択された前記重複領域の一対を解析して前記第１の画像群に対する前記第２の画像群のシフト量を求める工程と、
前記シフト量に基づいて各々の画像間の位置合わせを行い、前記第１の画像群と前記第２の画像群とを接続する工程と、
を備えることを特徴とする画像処理方法。
【請求項２】
同一平面内に相互に間隔を空けて複数の撮像素子を配列した撮像装置を用いて、前記各々の撮像素子の撮像範囲が一部重複するように前記平面に結像する被写体像と前記撮像装置との相対位置をシフトさせて、それぞれ撮像した第１の画像群と第２の画像群とを用いる画像処理方法であって、
前記撮像装置から前記第１の画像群のデータと前記第２の画像群のデータとを取得する工程と、
前記第１の画像群および前記第２の画像群における画像の重複領域のうち、前記平面内で相互の距離が最も離れた前記重複領域の一対を選択し、選択された前記重複領域の一対を解析して前記第１の画像群に対する前記第２の画像群のシフト量を求める工程と、
前記シフト量に基づいて各々の画像間の位置合わせを行い、前記第１の画像群と前記第２の画像群とを接続する工程と、
を備えることを特徴とする画像処理方法。
【請求項３】
請求項１または請求項２に記載の画像処理方法において、
前記シフト量は、前記第１の画像群および前記第２の画像群における前記重複領域の画像のデータまたは前記画像から抽出される特徴点のデータを用いた相関計算によって求めることを特徴とする画像処理方法。
【請求項４】
請求項２または請求項３に記載の画像処理方法の各工程を、コンピュータに実行させることを特徴とするプログラム。
【請求項５】
同一平面内に相互に間隔を空けて配列される複数の撮像素子を有し、前記平面に結像する被写体像を前記各々の撮像素子で撮像する撮像部と、
前記各々の撮像素子の撮像範囲が一部重複するように、前記被写体像と前記撮像部との相対位置をシフトさせるシフト部と、
前記相対位置をシフトさせて前記撮像部でそれぞれ撮像した第１の画像群および第２の画像群における画像の重複領域のうち、前記平面内で相互の距離が最も離れた前記重複領域の一対を選択し、選択された前記重複領域の一対を解析して前記第１の画像群に対する前記第２の画像群のシフト量を求める演算部と、
前記シフト量に基づいて各々の画像間の位置合わせを行い、前記第１の画像群と前記第２の画像群とを接続する画像接続部と、
を備えることを特徴とする画像処理装置。
【請求項６】
同一平面内に相互に間隔を空けて複数の撮像素子を配列した撮像装置を用いて、前記各々の撮像素子の撮像範囲が一部重複するように前記平面に結像する被写体像と前記撮像装置との相対位置をシフトさせて、それぞれ撮像した第１の画像群と第２の画像群とを用いる画像処理装置であって、
前記第１の画像群のデータと前記第２の画像群のデータとを取得する画像入力部と、
前記第１の画像群および前記第２の画像群における画像の重複領域のうち、前記平面内で相互の距離が最も離れた前記重複領域の一対を選択し、選択された前記重複領域の一対を解析して前記第１の画像群に対する前記第２の画像群のシフト量を求める演算部と、
前記シフト量に基づいて各々の画像間の位置合わせを行い、前記第１の画像群と前記第２の画像群とを接続する画像接続部と、
を備えることを特徴とする画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の撮像素子を有する撮像部と被写体像との相対位置を変化させて被写体像の複数の画像を撮像し、これら複数の画像を合成して広視野化を図る画像処理方法、画像処理装置およびプログラムに関する。
【背景技術】
【０００２】
従来からパノラマ画像合成等のように、撮像範囲が一部重複するようにシフトさせながら撮像した複数の画像を繋ぎ合わせることで、撮像装置の視野よりも広い範囲の合成画像を得る手段が公知である。
【０００３】
複数の画像を接続して合成画像を得る方法としては、例えば、重複する領域の画像を用いて相関を求める方法（特許文献１）や、重複する領域の画像に含まれるエッジ成分等の特徴量に着目して相関を求める方法（特許文献２）等が知られている。
【特許文献１】特許第２８８８７１３号
【特許文献２】特許第３４３２７３６号
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかしながら、特許文献１や特許文献２のような画像合成方法では、複数の画像を繋ぎ合わせるために重複する全ての領域で相関処理を行う必要がある。したがって、接続する画像の数に応じて演算量が著しく増加し、また画像処理を行うコンピュータ等に要求される性能も高くなってしまう。
【０００５】
上記従来技術が有する問題に鑑み、本発明の目的は、複数の撮像素子を用いて取得される複数の画像を、より少ない演算量で精度良く接続する画像処理の技術を提供することにある。
【課題を解決するための手段】
【０００６】
第１の発明に係る画像処理方法は、同一平面内に相互に間隔をおいて複数の撮像素子を配列した撮像装置を用いて、同一平面に結像する被写体像を撮像して第１の画像群を生成する工程と、各々の撮像素子の撮像範囲が一部重複するように、被写体像と撮像装置との相対位置をシフトさせる工程と、シフト後に撮像装置で被写体像を撮像して第２の画像群を生成する工程と、第１の画像群および第２の画像群における画像の重複領域のうち、同一平面内で相互の距離が最も離れた重複領域の一対を選択し、選択された重複領域の一対を解析して第１の画像群に対する第２の画像群のシフト量を求める工程と、シフト量に基づいて各々の画像間の位置合わせを行い、第１の画像群と第２の画像群とを接続する工程とを備えることを特徴とする。
【０００７】
また、第２の発明は、同一平面内に相互に間隔をおいて複数の撮像素子を配列した撮像装置を用いて、各々の撮像素子の撮像範囲が一部重複するように同一平面に結像する被写体像と撮像装置との相対位置をシフトさせて、それぞれ撮像した第１の画像群と第２の画像群とを用いる画像処理方法である。この画像処理装置は、撮像装置から第１の画像群のデータと第２の画像群のデータとを取得する工程と、第１の画像群および第２の画像群における画像の重複領域のうち、同一平面内で相互の距離が最も離れた重複領域の一対を選択し、選択された重複領域の一対を解析して第１の画像群に対する第２の画像群のシフト量を求める工程と、シフト量に基づいて各々の画像間の位置合わせを行い、第１の画像群と第２の画像群とを接続する工程とを備えることを特徴とする。
【０００８】
なお、上記発明に関する構成を、画像処理装置や、コンピュータプログラムおよびプログラムを記憶した記憶媒体等に変換して表現したものも本発明の具体的態様として有効である。
【発明の効果】
【０００９】
本発明によれば、画像の重複領域のうち、互いの距離が最も離れた重複領域の一対における相関計算を行うことで、高い接続精度を保ちつつ画像処理の演算時間を大幅に短縮することができる。
【発明を実施するための最良の形態】
【００１０】
図１は、本実施形態の画像処理システムの全体構成を示す概要図である。本実施形態に係る画像処理システムは、顕微鏡によって結像された試料の像を分割して撮像した複数の画像を合成処理することで、広視野かつ高解像度の合成画像を生成する。当該画像処理システムは、画像処理制御ユニット１、顕微鏡ユニット２および撮像ユニット３とから成る。画像処理制御ユニット１は、顕微鏡ユニット２および撮像ユニット３にそれぞれ通信可能に接続される。
【００１１】
画像処理制御ユニット１は、本実施形態の画像処理システムの統括的な制御を行うコンピュータである。画像処理制御ユニット１は、ＣＰＵ１０、記憶装置１１および入出力Ｉ／Ｆ１２とから成る。さらに、入出力Ｉ／Ｆ１２には、ユーザからの入力を受け付ける入力装置１５と、画像処理の途中経過や処理結果を表示する出力装置１６とが接続される。ＣＰＵ１０は、ユーザからの入力に応じて、入出力Ｉ／Ｆ１２を経由して顕微鏡ユニット２および撮像ユニット３を制御する。なお、ＣＰＵ１０は、一般的な中央演算装置を用いることができる。記憶装置１１は、ハードディスク等を適宜用いることができる。一方、入力装置１５としてはキーボードやマウス等を、出力装置１６としてはモニタやプリンタ等をそれぞれ適宜選択して使用できる。
【００１２】
顕微鏡ユニット２は、顕微光学系（２１〜２５）、顕微鏡制御部２６および顕微鏡Ｉ／Ｆ２７とから成る。顕微鏡制御部２６は、画像処理制御ユニット１からの指示を顕微鏡Ｉ／Ｆ２７を通じて受け、それに従って顕微光学系を制御する。一方、顕微光学系は、照明光源２１、コンデンサレンズ２２、試料台２３、対物レンズユニット２４およびズームレンズ群２５とから成る。なお、図１は顕微鏡ユニット２の主要部分のみを示す。例えば、図１において、照明光源２１とコンデンサレンズ２２との間に設置される照度調節部や試料台２３におけるピント調節部等は省略されている。
【００１３】
顕微鏡ユニット２において、照明光源２１からの光は、コンデンサレンズ２２によって集光され、試料台２３の試料を照射する。そして、試料を透過した光束は対物レンズユニット２４およびズームレンズ群２５によって集光され、撮像ユニット３の撮像部３１の受光面に試料の像が結像する。
【００１４】
撮像ユニット３は、顕微鏡ユニット２から出力された光束に基づいて、試料の顕微鏡像を撮像する。撮像ユニット３は、撮像部３１、Ａ／Ｄ変換部３２、シフト機構部３３、メモリ３４、撮像制御部３５および通信Ｉ／Ｆ３６とから成る。なお、撮像部３１、Ａ／Ｄ変換部３２、シフト機構部３３、メモリ３４および通信Ｉ／Ｆ３６は、それぞれ撮像制御部３５と接続される。撮像制御部３５は、通信Ｉ／Ｆ３６を通じて受け取った画像処理制御ユニット１からの指示に従って、撮像部３１およびシフト機構部３３を制御し、その結果を画像処理制御ユニット１に返す。
【００１５】
撮像部３１は、図１のＺ軸方向である顕微鏡ユニット３の光軸に対して垂直なＸＹ平面上に同一性能および同一画素数の４つの撮像素子を有する。撮像素子は、ＣＣＤやＣＭＯＳの半導体のイメージセンサ等を適宜選択して用いることができる。図２は、顕微鏡ユニット２側から見た撮像部３１の配置を示す。各々の撮像素子４０ａ〜４０ｄは、相互に間隔を空けて図２のように配列されている。撮像部３１の各撮像素子４０ａ〜４０ｄは、顕微鏡ユニット２を通過した光束による結像を光電変換によってアナログの画像信号を生成する。
【００１６】
Ａ／Ｄ変換部３２は、撮像部３１の撮像素子４０ａ〜４０ｄからの画像信号の出力に対してＡ／Ｄ変換を行う。
【００１７】
シフト機構部３３は、顕微鏡ユニット２の光軸（図１のＺ軸方向）に対して垂直な図２のＸＹ面で撮像部３１をシフトさせる。シフト機構部３３の動作により、顕微鏡ユニット２による試料の像と撮像部３１の受光面との相対位置は変化する。
【００１８】
メモリ３４は、撮像部３１が取得した試料の顕微鏡像の画像を一時的に格納し保持する。
【００１９】
なお、本実施形態では、撮像部３１の各々の撮像素子の間隔は、シフト機構部３３による撮像部３１の１回のシフト動作によって、各々の撮像素子の撮像範囲が重複するように設定されている。各々の撮像範囲の重複する幅は全て等しく、本実施形態では、上記の重複する幅は、例えば、撮像素子の撮像範囲の約１０％になるように調節されている。サンプリング定理で保証される範囲であれば、重複する幅をこれよりも小さくすることが可能である。
【００２０】
次に、本実施形態に係る画像処理システムにおける処理の手順について、図３のフローチャートおよび図４を参照しながら説明する。なお、図４は、試料像５０に対する撮像部３１のシフト動作に伴って、撮像素子４０ａ〜４０ｄの撮像位置がどのように変化するかを示している。なお、図４は、４つの撮像素子４０ａ〜４０ｄのうち、撮像素子４０ａの撮像位置の変化のみを示している。撮像素子４０ａ〜４０ｄの相互の位置関係は一定に保たれることから、他の撮像素子４０ｂ〜４０ｄの撮像位置は、撮像素子４０ａのものと同様の変化を示す。
【００２１】
具体的には、ユーザが、入力装置１５を使って画像処理の初期設定等とともに、画像処理プログラムのコマンドを入力し開始命令を出す。画像処理制御ユニット１のＣＰＵ１０は、その命令を入出力Ｉ／Ｆ１２を通じて受け、記憶装置１１に格納されている撮像制御プログラムを実行する。その結果、図３のステップＳ１０１からステップＳ１０８の処理が行われる。
【００２２】
ステップＳ１０１では、ＣＰＵ１０が、顕微鏡ユニット２の顕微鏡制御部２６に対して試料台２３上にある試料の像を撮像面３１に結像させるように指示する。同時に、ＣＰＵ１０は、シフト機構部３３に対して、撮像部３１の位置をシフトさせる指令を出す。これにより、撮像素子４０ａ〜４０ｄの位置は、図４の斜線で示された初期位置となる。
【００２３】
ステップＳ１０２では、撮像ユニット３の撮像制御部３５が、撮像部３１に対して第１画像の撮像命令を出す。撮像部３１は、図４に示す初期位置で試料像５０の一部を撮像する。撮像制御部３５は、Ａ／Ｄ変換部３２によってデジタル信号に変換された４つの第１画像をメモリ３４に一時的に格納する。ＣＰＵ１０は、撮像制御部３５に対し、メモリ３４に格納された４つの第１画像を、画像処理制御ユニット１へ転送するよう指示する。その後、４つの第１画像のデータは記憶装置１１に保存される。
【００２４】
ステップＳ１０３では、ＣＰＵ１０が、撮像ユニット３の撮像制御部３５に対して、シフト機構部３３に撮像部３１をシフトさせる指示を行う。これにより、撮像素子４０ａは図４の第２画像の撮像位置に移動する。
【００２５】
ステップＳ１０４では、ステップＳ１０２と同様に、撮像制御部３５が、撮像部３１に対して第２画像の撮像を指示する。ＣＰＵ１０は、撮像制御部３５に対して４つの第２画像を画像処理制御ユニット１へ転送するよう指示し、４つの第２画像を記憶装置１１に保存する。
【００２６】
ステップＳ１０５では、ＣＰＵ１０が、各撮像素子４０ａ〜４０ｄによって取得された第１画像および第２画像において、各々の画像が重複する領域のうちで、互いに最も離れた重複領域を選択する。図５は、撮像素子４０ａ〜４０ｄで撮像された第１画像（実線の四角形）と第２画像（斜線部の四角形）との間で生じる６つの重複領域６０ａ〜６０ｆを示す。ＣＰＵ１０は、この６つの重複領域のうちで、互いに最も距離が離れた重複領域の組み合わせ（重複領域６０ａと重複領域６０ｆ、または重複領域６０ｃと重複領域６０ｄの一方）を選択する。
【００２７】
ステップＳ１０６では、ＣＰＵ１０は、ステップＳ１０５で選択された一対の重複領域を用いて、撮像素子４０ａ〜４０ｄで撮像された第１画像と第２画像との位置関係を求めるための相関計算を行う。本実施形態では、ＣＰＵ１０は、相関計算において、重複領域の全範囲を演算対象とする。また、相関計算は画像で行っても良く、あるいは各々の画像から抽出されるエッジ成分等を画像の特徴量として用いても良い。そして、ＣＰＵ１０は、重複領域同士の位置をずらして各々の位置でパターンマッチングを行う。マッチングの度合いを示す相関値が最大となる位置を基に、ＣＰＵ１０は、第１画像に対する第２画像のシフト量を求め、シフト量データとして記憶装置１１に保持する。
【００２８】
ステップＳ１０７では、ＣＰＵ１０が、各撮像素子による第１画像から第４画像までの全ての画像が取得されているか否かを判定する。全ての画像が取得されていれば、一連の撮像を終了しステップＳ１０８（ＹＥＳ側）へ移行する。一方、全ての画像が取得されていないならば、ステップＳ１０３（ＮＯ側）に戻る。そして、ＣＰＵ１０は、上述したステップＳ１０３からステップＳ１０６までの動作を繰り返して、各撮像素子による第３画像および第４画像のデータ等を取得する。
【００２９】
ステップＳ１０８では、ＣＰＵ１０は、ステップＳ１０６で求めたシフト量データを用いて、撮像部３１の各撮像素子４０ａ〜４０ｄで撮像した第１画像から第４画像全てを繋ぎ合わせ、試料像５０の一枚の合成画像を生成する。ＣＰＵ１０は、生成された合成画像を記憶装置１１に保存するとともに、出力装置１６にその合成画像を出力表示して、画像処理を終了する。
【００３０】
なお、撮像素子４０ａ〜４０ｄによる第３画像取得の際のステップＳ１０５において、撮像素子４０ａ〜４０ｄによって撮像された第２画像（実線の四角形）と第３画像（斜線部の四角形）とによる重複領域は、図６のように７０ａ〜７０ｆで示される領域となる。そして、ＣＰＵ１０は、互いに最も離れた重複領域として、重複領域７０ａと重複領域７０ｆ、または重複領域７０ｂと重複領域７０ｅとの組み合わせにいずれかを選択する。
【００３１】
同様に、撮像素子４０ａ〜４０ｄによる第４画像取得の際のステップＳ１０５において、撮像素子４０ａ〜４０ｄによって撮像された第３画像（実線の四角形）と第４画像（斜線部の四角形）とによる重複領域は、図７のように８０ａ〜８０ｆで示される領域となる。そして、ＣＰＵ１０は、互いに最も離れた重複領域として、重複領域８０ａと重複領域８０ｆ、または重複領域８０ｃと重複領域８０ｄとの組み合わせのいずれかを選択する。なお、本実施形態では、互いに最も離れた重複領域を用いて計算を行ったが、計算不能の領域が含まれている場合には、次に互いの距離が離れた領域を用いて計算を行っても良い。
【００３２】
このように本実施形態は、複数の撮像素子の間隔を機械的に固定された撮像装置によって撮像した画像を用い、互いに最も距離が離れた重複領域における画像のみに着目し、相関計算を行う。これにより、２カ所それぞれの重複領域で求めたシフト量と重複領域間の距離との対比に基づいて各４つ１組の画像相互のシフト量を求めることから、全ての重複領域における相関処理と比較してより高い精度で画像相互の位置合わせが可能となる。
【００３３】
また、本実施形態において、全ての重複領域で相関計算を行わないことから、６つ全ての重複領域での相関計算を行う場合と比較して、処理時間が１／３に短縮される。撮像素子をＮ個用いるとすると、処理時間は従来技術と比較して２／Ｎになる。すなわち、位置合わせの時の演算負荷を大幅に低下して、高速な演算処理が可能となる。また、処理装置の回路規模を小型化することも可能となる。
≪本実施形態の補足事項≫
本実施形態では、４つの撮像素子を有する撮像部３１を用いたが、本発明の構成はこれに限定されない。撮像部３１に設置する撮像素子の数は、被写体像の視野の大きさ、撮像素子の撮像面の大きさ、要求される被写体像の視野をスキャンする時間等の状況に応じて適宜変更することができる。
【００３４】
なお、本実施形態では、撮像部３１の撮像素子４０ａ〜４０ｄの間隔を、撮像素子の撮像面の幅より狭い間隔にしたが、本発明はこれに限定されない。撮像部３１に設置する各撮像素子の間隔は適宜変更することができる。例えば、４つの撮像素子４０ａ〜４０ｄを、撮像部３１による２回のシフト動作で、撮像素子４０ａの撮像範囲が撮像素子４０ｂまたは撮像素子４０ｃの撮像範囲と重複するように設定しても良い。撮像部３１の撮像素子４０ａ〜４０ｄの間隔を変更した場合には、本実施形態で示したシフトパターンを含む様々なシフトパターンのうち最適なものを適宜選択し用いることが好適である。
【００３５】
なお、本実施形態では、シフト量を求める相関計算を、第１画像から第４画像をそれぞれ取得する毎に行ったが、本発明はこれに限定されない。例えば、全ての画像を取得した後に、一括して相関計算から画像合成までの処理を行うにしても良い。
【００３６】
なお、本実施形態では、重複領域の全範囲を対象としてシフト量を求める相関計算を行う例を説明したが、本発明はこれに限定されない。例えば、重複領域のうちの一部に着目して、部分画像のデータまたは部分画像から抽出される特徴量を用いて相関計算を行っても良い。
【００３７】
なお、本実施形態では、試料像５０と撮像部３１との相対位置を変化させるために、撮像部３１をシフトさせたが、本発明はこれに限定されない。例えば、レンズ等の光学部材または試料台２３をシフトさせて上記の相対位置を変化さるようにしても良い。
【００３８】
なお、本実施形態では、顕微鏡ユニット２として、透過照明による位相差顕微鏡を用いて行ったが、本発明はこれに限定されない。例えば、落射照明による顕微鏡に用いることもできる。また、例えば、顕微鏡以外の光学装置（望遠鏡等）によって結像される被写体像を分割撮像する場合にも本発明を適用できる。
【００３９】
なお、本発明は、画像処理システムの画像処理制御ユニット１、顕微鏡ユニット２および撮像ユニット３を一体とした画像処理装置に対しても適用可能である。
【００４０】
なお、本発明は、画像処理システムの画像処理制御ユニット１としてのプログラムを備えたコンピュータに対しても適用可能である。
【００４１】
なお、本発明は、画像処理システムの画像処理制御ユニット１としてのコンピュータプログラムを記憶する記録媒体に対しても適用可能である。
【００４２】
なお、本発明は、その精神またはその主要な特徴から逸脱することなく他の様々な形で実施することができる。そのため、上述した実施形態はあらゆる点で単なる例示にすぎず、限定的に解釈されてはならない。本発明は、特許請求の範囲によって示されるものであって、本発明は明細書本文には何ら拘束されない。さらに、特許請求の範囲の均等範囲に属する変形や変更は、全て本発明の範囲内である。
【産業上の利用可能性】
【００４３】
本発明は、以上説明したように、広視野で高解像度の画像の高速処理が要求される、顕微鏡等の画像の画像処理装置等に応用可能な技術である。
【図面の簡単な説明】
【００４４】
【図１】実施形態の画像処理システムの全体構成を示す概要図
【図２】撮像部３１における４つの撮像素子の配置例を示す図
【図３】本発明に係る画像処理装置における手順のフローチャート
【図４】撮像部３１の受光面に結像する試料像５０に対する４つの撮像素子の初期位置と撮像素子４０ａの撮像位置の変化を示す図
【図５】第１画像および第２画像による重複領域を示す図
【図６】第２画像および第３画像による重複領域を示す図
【図７】第３画像および第４画像による重複領域を示す図
【符号の説明】
【００４５】
１  画像処理制御ユニット、２  顕微鏡ユニット、３  撮像ユニット、１０  ＣＰＵ、１２  入出力Ｉ／Ｆ、１５  入力装置、１６  出力装置、３１  撮像部、３４  シフト機構部、３５  撮像制御部、３６  通信Ｉ／Ｆ
【図１】
2009031909.tif 000003
【図２】
2009031909.tif 000004
【図３】
2009031909.tif 000005
【図４】
2009031909.tif 000006
【図５】
2009031909.tif 000007
【図６】
2009031909.tif 000008
【図７】
2009031909.tif 000009