(57)【特許請求の範囲】
【請求項１】
光学系を通してスキャナ中にある複数の角度における輝度の遷移を持つ光学校正パターンを走査するステップと、
走査した前記光学校正パターンから点拡がり関数を計算するステップと、
画像画素の走査線を得る為に画像を走査するステップと、
前記点拡がり関数を用いて前記走査線の対応区分を修正し、前記光学系に起因する収差の幾分かを除去するステップと、を含む画像スキャナを用いた画像走査方法。
【請求項２】
前記複数の角度における輝度の遷移を持つ光学校正パターンは、白黒交互のくさび型が放射状に並んだ半円形が配列されていることを特徴とする請求項１に記載の画像スキャナを用いた画像走査方法。
【請求項３】
前記複数の角度における輝度の遷移を持つ光学校正パターンは、同心円が配列されていることを特徴とする請求項１に記載の画像スキャナを用いた画像走査方法。
【請求項４】
画像の走査前に光学系の第１の点拡がり関数を決めるステップと、
前記光学系を通して、画像画素を得るために前記画像を走査するステップと、
前記画像の走査中に、スキャナ内にある複数の角度における輝度の遷移を持つ複数の光学校正パターンの少なくとも１つを前記光学系を通して走査するステップと、
前記画像の走査中に走査された前記光学校正パターンの１つから第２の点拡がり関数を計算するステップと、
前記光学系が変化したかどうかを決定するため、前記第１の点拡がり関数と前記第２の点拡がり関数を比較するステップとを含む画像スキャナにおける光学系をモニタする方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、一般的にはイメージスキャナ、デジタルコピー機、ファクシミリ及びその類似機器に用いられる光学系に関し、より具体的には、スキャナ内部にてリアルタイムで光学系の点拡がり関数を計測した後に光学収差を補償することに関する。
【０００２】
【従来の技術】
イメージスキャナは原稿又は写真上の可視画像、或は透明媒体中の画像を、コンピュータにより複写、記憶又は処理するに適した電子形式へと変換する。画像スキャナは独立した装置、或はコピー機やファクシミリの一部、また或は多目的装置の一部である場合もある。反射式イメージスキャナは代表的には制御された光源を有し、光は原稿の表面で反射し、光学系を通って感光性素子のアレイ上に至る。感光性素子は受けた光の強度を電子信号へと変換する。透過式イメージスキャナでは、光は例えばポジ写真スライドのような透光性の画像を透過し、光学系を通って感光性素子のアレイ上へと至る。
【０００３】
一般的にスキャナは、光学レンズ系又は光学導波管系を利用して画像の焦点を光センサーアレイ上に合わせる。レンズ系及び光学導波管系では例えば球面収差やコマ収差、非点収差、視野曲率収差、色収差、動きぶれ、迷光反射（ゴースト）等といった様々な光学収差問題が起こり得る。収差は光学系をより精密かつ複雑にすることにより低減することが可能であるが、そうするとコストや大きさ、重量を増す結果となる。比較的単純で低コストの光学系の利用を可能にすることにより、全体的なシステムコストやサイズを低減する必要がある。
【０００４】
周知の固定収差はある程度、幾何学的校正により排除することが出来る。例えば魚眼レンズにより作成された画像の空間画素ずらし値を決定する為に、ドット格子の画像を利用することが出来る。何種類かの収差に対する他の補償値は、レンズ及び光学デザインソフトウエアの数学的モデルを用いて決定することが出来る。
【０００５】
光学収差の補償を実行する為の他の方法では、線形系理論及び系のインパルス応答が使われる。光学系において、インパルス応答、即ち暗い背景に対する高輝度点オブジェクトの画像を点拡がり関数と呼ぶ。空間領域点拡がり関数の２次元フーリエ変換を光学伝達関数と呼ぶ。周波数領域光学伝達関数は大きさと相を含む。２次元、大きさのみ、変調伝達関数の縦割りは変調伝達関数と呼ばれる。画像が光学系によりぶれた場合、ぶれた画像の変換値をその同じ不完全光学部品の光学伝達関数で割ることでぶれの一部を周波数領域中で除去することが出来る。空間領域におけるコンボリューション（たたみ込み積分）は周波数領域における乗算に等しい。よって、かわりに点拡がり関数から計算した適切な２次元核でコンボリューションすることにより、空間領域におけるぶれを幾分取り去ることも可能である。エイリアシングの為、そしてＳ／Ｎ比と画像シャープニングとの間に可能な所望のトレードオフの為に、所望の変調伝達関数は、実際のインパルスに対する応答に基づく理想変調伝達関数よりも低い周波数応答を有する場合もある。
【０００６】
通常、レンズの数学的モデルを利用した収差の判定、又は点拡がり関数の判定は、軸上画像形成の為の１度きりの校正として実行される。しかしながら、レンズの点拡がり関数は製造工程により、経時により、温度により、湿度により、そして画像焦点面と対象物焦点面の位置によって変化する。光学系が実質的にコンピュータモデルと異なる場合、又は点拡がり関数が判定された後にレンズ系が変化した場合、収差を補償する為のフィルタリング処理は、画像の劣化を抑えるのではなく、むしろその画像を更に悪化させてしまう結果となる。
【０００７】
具体的には、低コストのプラスチックレンズには比較的に温度に敏感なものもある。画像スキャナは、一般的にかなりの熱を放射するランプを含んでいる。コピー機や多目的装置もまた、トナーを融解して用紙上に塗布する為の、又はインクを乾燥させる為のヒーターを含んでいるものがある。走査の実行要求以前にランプ及びヒーターをＯＦＦ、又は低電力スタンバイ状態にすることは出来る。走査の間の画像スキャナの内部温度は、ある時間の間（単一画像の走査必要な時間に対して考えると長い）に著しく変化する場合もある。更に小型のものにおいては、画像走査レンズ系は比較的レンズに近い原稿に焦点を合わせる為、広い視野にわたって正確でなければならない。収差のいくつか、例えば非点収差、コマ収差、像面湾曲等は、空間的に変化する。従って軸上での投影に基づいて一度判定された単一の核は画像スキャナ中のレンズ系については適切ではない場合も有り得る。
【０００８】
【発明が解決しようとする課題】
画像スキャナ製品において、単純性、低コスト性、そして完成された製品内部における収差の特性付け及び収差補償の迅速性が要求されている。加えて、スキャナ製品において、レンズの視野にわたって空間的に変化する収差を補償することが必要である。更に、温度、湿度及びその他の変化要因を調整する為に要する補償値を各走査処理前に決定することが必要とされている。最後に、走査中の収差をモニタし、必要であれば再走査することが必要である。
【０００９】
本発明は上記事情に鑑みて成されたもので、スキャナ内において、ほぼリアルタイムで測定を行い、デジタル画像処理法を利用してその不完全性を部分的に補償することを目的とする。
【００１０】
【課題を解決するための手段】
上記目的を達成するために、１実施例においては、一連の２次元光学ターゲットが原稿走査領域外に設けられている。ターゲットは走査線の全体長にわたって伸び、レンズ系の視野全体にわたる補償を可能にしている。他の実施例においては、ターゲットが原稿の側辺に沿って置かれており、走査中の補償をモニタすることが出来るようになっている。
【００１１】
個別のターゲットは、それぞれ走査線の小区分に対する光学伝達関数の推測値を得るに適している。ターゲットの一例には、複数の角度における輝度（白から黒、又は黒から白）の段階関数を持つものもある。各段階関数は１つの変調伝達関数の推測値を得るに適している。１つの個別ターゲットに対する複数の変調伝達関数の集合がその個別ターゲットの空間位置に対応する視角の光学伝達関数を推測する為に利用される。
【００１２】
他の実施例においては、２次元の既知のランダムパターンがターゲットとして設けられている。既知のランダムパターンの一部分を既知のランダムパターンの同じ部分の走査画像と相互相関すると、その既知のランダムパターンの部分に対応する走査線の部分の点拡がり関数を推定することが出来る。
【００１３】
次に個々の光学伝達関数をそれぞれ用いて光学収差を補償する。１実施例においては、光学伝達関数は走査線の対応区分のコンボリューション核を計算する為に利用される。走査線の各区分内の走査画像画素はその区分に適切な核でコンボリューションされる。他の実施例においては、補償画像の計算に反復（イタレーション）的ソリューションが用いられている。
【００１４】
一連のターゲット、又は走査線の幅にわたって連続ランダムターゲットを設けることにより、スキャナ内で光学伝達関数を走査に適切な湿度及び温度下での組立レンズの位置の関数として判定することが出来るようになる。更なるターゲットを原稿の側辺に沿って設けると、走査進行中に確認を行なうことが可能である。結果として、より小型で低コストのレンズを利用することが可能となり、最終走査画像から画像劣化をある程度排除することが出来るようになる。
【００１５】
【発明の実施の形態】
図１は一例としてあげる画像スキャナ（１００）の一部を描いたものである。スキャナ（１００）は透明プラテン（１０２）を含んでいる。画像作成の対象である原稿（１０４）は透明プラテン上に下向きに置かれている。光線（１０６）はランプ（１０８）から発され、原稿の表面で反射し、その後ミラー（１１０、１１２、１１４）によって反射され、レンズ系（１１６）を通過して最終的に光センサアレイ（１１８）へと至る。その他の構成も一般的である。通常、レンズ系（１１６）は複数の要素から成るもので良い。本発明に基づくスキャナにおいては、レンズ系（１１６）の光学性能はスキャナ内部で特性付けられ、スキャナはデジタル画像処理を利用してレンズ収差を部分的に補償する。特に重要なのは、レンズ校正ターゲット（１２０）である。ターゲット（１２０）は、出来ればターゲット上の校正パターンが原稿面と同じ焦点面に来るように、また原稿（１０４）の画像形成を妨害しないように置かれることが望ましい。
【００１６】
図２はターゲット（１２０）の１実施例である。ターゲット（１２０）はスキャナ（１００）により走査可能の最高原稿幅と少なくとも同寸の幅を持つものが望ましい。図２の実施例においては、６個の同一パターン（２００）が示されており、各々は白黒交互のくさび型が放射状に並んだ半円形の配列である。６個というのは単に説明の便宜上そうしているだけであって、実際のターゲットはそのレンズ系が必要とする数に応じて６個よりも少なくても多くても良い。
【００１７】
ターゲット（１２０）上のパターンの各々は、エイリアシングの為に黒から白（及び白から黒）への輝度の遷移は真の段階関数ではなく、かわりに既知の変化率で迅速なもので、理想的にはターゲットの変調伝達関数が既知であるべきである。輝度遷移端の各々は、ノイズを低減する為にその放射の伸びに沿って平均化されている。端拡がり関数を各端について計算することが出来る。結果として得られた端拡がり関数を、特定の１角における点拡がり関数の１次元成分を算出する為に微分することが出来る。例えば、Ｋｅｎｎｅｔｈ  Ｒ．Ｃａｓｔｌｅｍａｎによる「Ｄｉｇｉｔａｌ  Ｉｍａｇｅ  Ｐｒｏｃｅｓｓｉｎｇ」（Ｐｒｅｎｔｉｃｅ  Ｈａｌｌ、１９９６）の第１６章を参照されたい。点拡がり関数が円対称であった場合、単一の端から計算された１次元点拡がり関数を回転させて２次元点拡がり関数を得ることが出来る。しかしながら、一般的には正確に２次元点拡がり関数を決定するには様々な方向を向いた数多くの端が必要となる。
【００１８】
図３は他のターゲットパターンを描いたものである。図３においては、ターゲットは同心円である。複数の遷移でノイズを平均化することが出来、同心円により任意のいずれの角における端拡がり関数も測定可能となる。
【００１９】
図２及び図３に描かれた一連の不連続校正パターンのかわりに、連続したランダムの２次元パターンを持つ校正ターゲットを用いても良い。空間領域において適したパターンを、例えば一定の振幅とランダムな相を持つ２次元スペクトルの２次元逆フーリエ変換により得ることが出来る。既知のランダムパターン及びそのパターンの走査画像との相互相関値は、その系の点拡がり関数である。走査線の異なる区分の点拡がり関数を決定する上で必要に応じてランダムパターンを論理的に様々なサイズのセクションへと分割することが出来る。
【００２０】
点拡がり関数を補償の為に利用する場合、レンズ系は線形系であると想定される。原画像ｆ（ｘ，ｙ）及びレンズ系の空間領域点拡がり関数ｈ（ｘ，ｙ）が与えられると、結果として得られる走査画像ｇ（ｘ，ｙ）は、
ｇ（ｘ，ｙ）＝ｈ（ｘ，ｙ）＊ｆ（ｘ，ｙ）
であり、ここで「＊」は２次元コンボリューションを示す。離散形式では以下の通りである。
【００２１】
【数１】
0004021594.tif 000002
【００２２】
走査画像ｇ（ｘ，ｙ）が与えられ、ｈ（ｘ，ｙ）が与えられると、原画像はデコンボリューションにより復元（ぶれは減少している）される。
ｆ（ｘ，ｙ）＝ｓ（ｘ，ｙ）＊ｇ（ｘ，ｙ）
ここでｓ（ｘ，ｙ）は空間領域核（デコンボリューションマトリクス、又は逆フィルタ）であり、そして場合によっては所望の点拡がり関数もまた組み込む為に変更される。周波数領域における系の点拡がり関数ｈ（ｕ，ｖ）が与えられると、ｓ（ｘ，ｙ）は１／（ｈ（ｕ，ｖ））の逆フーリエ変換を計算することにより得ることができる。かわりに空間領域の点拡がり関数ｈ（ｘ，ｙ）が与えられた場合は、ｓ（ｘ，ｙ）は一連の連立方程式を解くことにより計算することが出来る。
【００２３】
かわりに、ｓ（ｘ，ｙ）を特に計算するのではなく、イタレーションによってｆ（ｘ，ｙ）を直接的に算出することが出来る。例えば、Ｍｉｃｈａｅｌ  Ｅｌａｄ及びＡｒｉｅ  Ｆｅｕｅｒによる「Ｒｅｓｔｏｒａｔｉｏｎ  ｏｆ  ａ  Ｓｉｎｇｌｅ  Ｓｕｐｅｒｒｅｓｏｌｕｔｉｏｎ  Ｉｍａｇｅ  ｆｒｏｍ  Ｓｅｖｅｒａｌ  Ｂｌｕｒｒｅｄ，  Ｎｏｉｓｙ、ａｎｄ  ＵｎｄｅｒｓａｍｐｌｅｄＭｅａｓｕｒｅｄ  Ｉｍａｇｅｓ」（画像処理に関するＩＥＥＥ発表論文Ｖｏｌ．６、Ｎｏ．１２、ｐｐ１６４６−１６５８、１９９７年１２月）を参照されたい。Ｓｔｅｅｐｅｓｔ  Ｄｅｓｃｅｎｔアルゴリズムを用いて、次の反復等式を使うことが出来る。
【００２４】
【数２】
0004021594.tif 000003
【００２５】
ここでｈは上記に説明した空間領域点拡がり関数であって、初期化ベクトルｆ−０はいかなるベクトルでも良く、そしてμ＞０は小さい。
【００２６】
反復的手法の利点は、点拡がり関数補償を２次元再サンプリングと組み合わせることが出来る点にある。しばしば、所望の走査画像が、可能な最高光学サンプリングレートよりも低いサンプリングレートを有する場合もある。高解像度のカラー画像には大きい記憶容量が必要とされ、そしてその画像分析又は通信チャンネル上での伝送には多大な時間が必要である。よってオペレーターが容量や処理時間又は伝送時間を削減する為に低いサンプリングレートを選択する場合もある。画像を文書又はディスプレイに入るように調整する為に、プリンタやディスプレイ装置の解像度で決まるサンプリングレートでスケーリングしなければならない場合もある。結果として、所望のサンプリングレートとは、いずれかの任意の数字となる。最後に幾つかの走査コンフィギュレーションにおいては、サンプリングレートが下げられると走査時間を短くすることが出来るものがある。走査時間は、特定のＳ／Ｎ比を得る為に必要な画素毎の露出時間によってしばしば限定される。隣接するＣＣＤ素子の電荷を併せることにより、より高速ではあるが低サンプリングレートでの走査が可能である。
【００２７】
上記のいずれかの理由から、スキャナのオペレーターが最高光学サンプリングレートの何分の１かのサンプリングレートを要求することがしばしばある。代表的には、単一の走査線のサンプルがスキャナにより最高光学サンプリングレートで読み取られ、中間の補間サンプルが算出され、そして補間サンプルを再サンプリングすることにより所望のサンプリングレートが得られる。例えば、スキャナは１インチあたり１２００画素の光学サンプリングレートを持ち、１インチあたり４８００画素にまで補間することができ、そして４８００画素／インチのサンプルを再サンプリングしてインチあたり９３８画素（９３８は任意の数の説明上の一例にしか過ぎない）の画像を作ることが出来るといったようなことである。移動式リニアセンサアレイを有するスキャナにおいては、補間及び再サンプリングにより走査線に沿った任意の解像度の提供が可能であり、走査方向の解像度はセンサアレイの変換速度により制御することが出来る。かわりに、４個の隣接センサの強度値を用いてある画素の値を計算する為に、双線形補間法を利用しても良い。しかしながら、上記に触れたＥｌａｄの論文でも説明されているように、（点拡がり関数を用いることによる）ぶれの補間及び２次元再サンプリングを反復的手法に組み合わせることも可能である。
【００２８】
上記にも説明したように、反復的手法には利点はあるものの、多くのスキャナは空間領域核処理用に専用ハードウエアサポートを有している。一般的に、点拡がり関数は核を計算する為に用いられる唯一の情報である。
【００２９】
端拡がりデータからの核の計算は、スキャナ内部で実施することが可能である。周辺装置である場合はかわりに端拡がりデータをホストコンピュータにアップロードすることも出来る。そうすればホストコンピュータは核を計算し、その核をスキャナへとダウンロードすることが出来る。また、画像スキャナは、例えば画像シャープニングの為等、画像を他の核でコンボリューションすることが必要な場合もある。様々な核を合わせてコンボリューションすることが出来、また、結果として得られる組み合わせられた核を走査画像データに適用することが出来る。幾つかの核情報はホストコンピュータにより判定することが出来る。装置の独立性の為に、スキャナが内部スキャナ核情報をホストコンピュータへとアップロードしなければならない場合もあり得る。するとホストコンピュータはスキャナに特定の核をホストコンピュータにより決められた核と共にコンボリューションし、結果の組み合わせられた核をスキャナへとダウンロードすることが出来る。例えば、米国特許出願番号第０８／７７５、０６１号を参照されたい。本発明によれば、スキャナは端拡がりデータを他の核情報に加えてホストコンピュータへとアップロードすることが出来る。
【００３０】
走査線がＮ個のパターンでＮ個の領域に分割されており、点拡がり関数は各領域内においては一定であると仮定する。各領域についてコンボリューション核が対応する点拡がり関数及び所望の変調伝達関数から計算される。核のサイズは出来ればいずれのぶれ範囲よりも大きいことが望まれる。例えば、１画素のぶれ範囲がいずれの方向に置いても４画素分を越えてはいない場合、９×９値の核を用いれば良い。空間領域において、コンボリューション核は各領域についてその領域の各画素に作用させる為に用いられる。原稿の端に近い画素はフィルタリングせずに残すことも、或は空白データを使って端を伸ばし、フィルタリングの範囲を端まで伸ばすようにすることも出来る。光学補償核が他の核と共にコンボリューションされると、他の核と共にコンボリューションされたその光学補償核は走査されている領域に依存することになる。
【００３１】
図４は、原稿の側辺に沿った更なる校正パターンが描かれてたスキャナの上面図である。図４でがパターン（４００）が原稿の一方の側にしか描かれていないが、原稿の両側に設けても良い。走査中、温度変化が比較的急激に起きた場合、走査前に計算された点拡がり関数は不適切となってしまう可能性がある。原稿の側辺に沿った更なるパターンにより、走査中に少なくとも１つの点拡がり関数が変化しているかどうかをチェックすることが可能となる。変化をモニターする目的であれば、点拡がり関数全体を計算する必要が無い場合もある。例えば、１つの角度における１つの端拡がり関数をモニタすれば充分な場合があるということである。測定した関数がリアルタイムで変化している場合、そしてレンズ系の変化が一貫している場合、図４のようなモニタ用パターンは点拡がり関数のリアルタイムでの変更を実現するに充分な情報を提供することが出来る。かわりに、単にスキャナが開始点に戻り、再校正及び再走査を行なうものであっても良い。
【００３２】
校正ターゲットをスキャナに組み込むと、解像度及び動きぶれを含むオペレーター制御変数の補償も容易になる。センサアレイと走査原稿との間に連続した相対的な移動がある場合、各画素の露光時間中の移動に起因していくらかのぶれが生じる。移動量が既知である場合、動きぶれを補償核に組み込むことが出来る。例えば上記で触れたＥｌａｄによる論文を参照されたい。多くのスキャナにおいては、走査線に直交する方向における光学サンプリングレートは相対的走査運動の速度により決定する。例えば、オペレーターがより高い解像度を望む場合、相対的走査運動は比較的に遅く、オペレータが低解像度を望む場合は相対的走査運動は比較的に高速となる。校正パターンがスキャナの内部にある場合（図２、図３）、パターンはオペレーターにより指定される（Ｘ、Ｙ）解像度で走査されなければならない。すると画像に特定の解像度及び動きぶれは、その校正プロセスに固有のものとなる。
【００３３】
つまり、複数の校正ターゲットパターンをスキャナ中に設ければ、光学特性がレンズ毎に異なったり温度や経時で変化する可能性がある、より低コストのレンズ系を利用出来るということである。加えて、校正ターゲットパターンをスキャナに設けることにより、画像に特異的な補償［例えば特定の（Ｘ、Ｙ）解像度及び動作ぶれ］の計算が容易になる。更に、事前に決められた点拡がりデータの大容量の表を記憶する為のＲＯＭを必要としない。
【００３４】
本発明の上記の説明は描写及び説明の目的の為に示したものである。それは本発明を開示の特定の形態に限定する為に網羅することを意図したものではなく、上記の説明に照らして他の改変及び変更が可能である。実施例は本発明の原理及びその実用的用途を最もわかりやすく説明する為に、また、そうすることにより当業者が企図する特定の利用法に適した様々な実施例及び様々な改変形態において本発明をうまく利用出来るように選択及び記述したものである。添付請求項は、従来技術により制約されるものを除き、他の実施例をも含むと解釈されることを意図するものである。
【００３５】
以下、本発明の実施の形態を要約して挙げる。
【００３６】
１．  光学系（１１６）を通してスキャナ中にある光学校正パターン（１２０）を走査するステップと、走査した光学校正パターンから点拡がり関数を計算するステップと、画像画素を得る為に画像（１０４）を走査するステップと、点拡がり関数を用いて画像画素を修正し、光学系に起因する収差の幾分かを除去するステップと、を含む画像スキャナを用いた画像走査方法。
【００３７】
２．  前記点拡がり関数を計算するステップが、更に点拡がり関数から空間領域デジタル核を計算するステップを含み、前記点拡がり関数を用いるステップが、更に画像画素を空間領域デジタル核でコンボリューションするステップを含む上記１に記載の画像スキャナを用いた画像走査方法。
【００３８】
３．  核の各々が走査線の１部分に対応する複数の異なるデジタル核を計算するステップと、デジタル核の各々をその核に対応する走査線の部分に入る画像画素でコンボリューションするステップと、を更に含む上記２に記載の画像スキャナを用いた画像走査方法。
【００３９】
４．  前記点拡がり関数を用いるステップが、光学系に起因する収差が幾分除去された新たな画像を形成する為に画像画素を修正する等式を反復するステップを更に含む上記１に記載の画像スキャナを用いた画像走査方法。
【００４０】
５．  光学系（１１６）と、該光学系の光学伝達関数を決定するのに適した少なくとも１つの光学ターゲット（１２０）と、を含む画像スキャナ。
【００４１】
６．  走査中、前記光学伝達関数を確認する為に利用可能なように少なくとも１つの追加の光学ターゲット（４００）が配置されている上記５に記載の画像スキャナ。
【００４２】
【発明の効果】
以上説明した本発明によれば、画像スキャンに先立ちスキャナ内の光学校正パターンをスキャンして点拡がり関数を得、この関数に従いスキャンされた画像の光学収差を補償することができる。
【図面の簡単な説明】
【図１】本発明に基づく校正ターゲットを有する画像スキャナの一例の側面の断面図である。
【図２】校正ターゲットの一例の平面図である。
【図３】他の校正ターゲットの平面図である。
【図４】図１のスキャナの上面図であり、原稿の側辺に沿った更なるターゲットを描いたものである。
【符号の説明】
１０４  画像
１１６  光学系
１２０  光学校正パターン（ターゲット）
４００  追加の光学校正パターン（ターゲット）
【図１】
0004021594.tif 000004
【図２】
0004021594.tif 000005
【図３】
0004021594.tif 000006
【図４】
0004021594.tif 000007