(57)【要約】【課題】  従来以上に低コストでシンプルな手段で、光学系の収差の影響を軽減することを実現した画像読取装置の提供。【解決手段】  アレイセンサのラスタ走査によって画像を読み取る際に、光学系の収差の影響をデジタル変倍処理によって軽減する。ラスタ走査の順に画素信号を入力し、入力の順に１つのラスタ走査の全画素よりも少ない限られた数の入力画素の画素値をメモリに記憶する。そして、収差の特性を反映して作成されたデータテーブルを用いて、次の出力画素の座標を求める。求めた座標に最も近い座標に対応する入力画素の画素値を、前記メモリから読み出し、少なくともこの画素値を用いて次の出力画素の画素値を求める。これを前記入力の順に出力する。【選択図】  図４
2010068124.tif 000002
【特許請求の範囲】
【請求項１】
アレイセンサと光学系を有する読取手段を用い、前記アレイセンサのラスタ走査によって画像を読み取って、画像データに変換して出力する画像読取装置であって、
前記光学系の収差の影響を、デジタル変倍処理によって軽減する信号処理手段を有し、
該信号処理手段は、
前記ラスタ走査の順に各画素の画素値を入力する入力手段；
１つのラスタ走査の全画素よりも少ない限られた数の入力画素の画素値を一時記憶するメモリ；
前記収差の特性を反映して作成されたデータテーブルを用いて、次の出力画素の座標位置を求める第１処理手段；および
前記第１処理手段で求めた座標位置に対応する座標位置の入力画素の画素値を、前記メモリから読み出し、少なくともこの画素値を用いて次の出力画素の画素値を求める第２処理手段
を有することを特徴とする画像読取装置。
【請求項２】
前記収差は、前記光学系の歪曲収差および色収差である、請求項１記載の画像読取装置。
【請求項３】
前記アレイセンサは、色成分毎に個別に読み取りを行なう複数のラインセンサを有し、前記信号処理手段は各色成分毎に前記処理を行なう、請求項１又は２記載の画像読取装置。
【請求項４】
前記メモリは、複数の色成分毎に、前記限られた数の画素値を格納するセルを有するＦＩＦＯメモリであり、
前記第２処理手段は、色成分毎に、前記第１処理手段で求めた座標位置に対応する最も近い座標位置の入力画素の画素値を保持する前記ＦＩＦＯメモリのセルを指定する手段を有する、請求項３記載の画像読取装置。
【請求項５】
前記第１処理手段は、１つのラスタ走査の全画素を複数の区間に区切って、前記収差の特性を反映して各区間における変倍率と画素数を設定して処理を行なう、請求項１乃至４のいずれか記載の画像読取装置。
【請求項６】
前記第１処理手段は、前記データテーブルを参照して、初期の座標位置に対して、１画素を処理する毎に所定の加算値を積算していくと共に、１つのラスタの処理の中で前記所定の加算値を変更する手段を有する、請求項１乃至５のいずれか記載の画像読取装置。
【請求項７】
前記データテーブルは、［加算値］と［処理画素数］のデータ組を複数有するものであり、ある［加算値］で［処理画素数］だけ画素を処理する毎に、前記データテーブルの次のデータ組を参照して、前記加算値を変更する、請求項６記載の画像読取装置。
【請求項８】
前記第１処理手段は、１つのラスタに対して前半と後半で対称的な処理を行なう、請求項１乃至７のいずれか記載の画像読取装置。
【請求項９】
前記第２処理手段は、前記第１処理手段で求めた座標位置に最も近傍した座標位置とその次に近傍した座標位置の２つの入力画素の画素値から、補間処理によって次の出力画素の画素値を求める手段を有する、請求項１乃至８のいずれか記載の画像読取装置。
【請求項１０】
前記第２処理手段は、前記第１処理手段で求めた座標位置に最も近傍した座標位置の入力画素の画素値を、次の出力画素の画素値とする、請求項１乃至８のいずれか記載の画像読取装置。
【請求項１１】
アレイセンサと光学系を有する読取系を用い、前記アレイセンサのラスタ走査によって画像を読み取って、画像データに変換して出力する画像読取方法であって、
前記光学系の収差の影響を、デジタル変倍処理によって軽減する信号処理を有し、
該信号処理は、
前記ラスタ走査の順に各画素の画素値を入力するステップ；
１つのラスタ走査の全画素を複数の区間に区切って、前記収差の特性を反映して各区間の変倍率を設定し、これに基づいて出力画素の座標位置およびその画素値を演算処理によって求めるステップ；および
前記求めた画素値を、前記演算処理の順に出力するステップ
を有することを特徴とする画像読取方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、光学スキャナ等の画像読取装置に関するものである。
【背景技術】
【０００２】
画像読取装置である光学スキャナには、等倍光学系を用いたいわゆるＣＩＳタイプと、縮小光学系を採用したいわゆるＣＣＤタイプがある。
【０００３】
ＣＣＤタイプのスキャナは、ＣＣＤラインセンサと縮小光学系を有する読取手段を用い、原稿画像をセンサのラスタ走査によって読み取って、各画素の信号から画像データに変換して出力する。ラインセンサによる主走査方向のラスタ走査と、原稿とラインセンサとの副走査方向の相対走査の組み合わせによって、二次元の画像を読み取る。これは例えば、特許文献１に開示されている。
【特許文献１】特開２０００−２４４７１７号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
ＣＣＤタイプのスキャナの小型化を進めるためには、内蔵する縮小光学系のコンパクト化が有効である。しかし、コンパクトにするために光学系の光路長を短くするほどに、光学設計の余裕がなくなって、光学系の各種の収差、たとえば歪曲収差や色収差が大きくなる。
【０００５】
収差が大きくなると、以下のような問題が生じる。
【０００６】
図１１に示すように、ラスタ走査の１ライン内において、読み取った画像の解像度が局所的に異なってしまう。さらに、原稿上の同じ位置に対応するセンサ画素の座標位置が、ＲＧＢの各色成分で異なってしまうため、読み取って得た画像データ上では色ズレとなってしまう現象が生じる。
【０００７】
図１２は、光路長を短くする代償として大きな収差を持つ縮小光学系を用いたスキャナ装置における、収差による影響を示す図である。解像度ｄ（ｄｐｉ）で読み取って得た画像の、位置ズレと色ズレが発生した状況を示すものである。
【０００８】
原稿の中心をＣＣＤセンサの主走査方向中心に合わせた場合を考える。図１２（ｃ）は、原稿上での距離がｉ×Ａインチ（ｉは整数、Ａは一定値）となるサンプル点Ｐ［ｉ］について、収差ゼロの理想の光学系で読み取った場合である。得られる画像は、Ｐ［０］の読み取り画像とＰ［ｉ］の読み取り画像間の距離（＝ピクセル数）はｉＡｄであり、ｉに比例する。この値をＱ［ｉ］とする。
【０００９】
しかし、歪曲収差のある光学系で読み取った画像ではこのような理想形とはならない。また、色収差によってＲＧＢの各色成分の間でも差異を生じる。図１２（ｄ）はＲＧＢの各色成分毎の画像データを示す。
【００１０】
この値をそれぞれＱｒ［ｉ］、Ｑｇ［ｉ］、Ｑｂ［ｉ］として、Ｄｒ［ｉ］＝Ｑ［ｉ］−Ｑｒ［ｉ］、Ｄｇ［ｉ］＝Ｑ［ｉ］−Ｑｇ［ｉ］、Ｄｂ［ｉ］＝Ｑ［ｉ］−Ｑｂ［ｉ］とする。すると、Ｄｒ［ｉ］、Ｄｇ［ｉ］、Ｄｂ［ｉ］（単位はピクセル）が、原稿上の点Ｐ［ｉ］における画素ズレ量であり、この３つの値の差が同じく点Ｐ［ｉ］における色ズレ量である。
【００１１】
図１２（ａ）のグラフは、Ｄｒ［ｉ］、Ｄｇ［ｉ］、Ｄｂ［ｉ］を、実際に歪みのある光学系で測定しプロットしたものである。このグラフから判るように、センサの中心では画素ズレ・色ズレともに存在しない。しかし、センサ端に近づくにつれて、センサ中心からの距離の比例倍よりも大きくなる、いわゆる画素ズレが発生する。また、ＲＧＢ間の色ズレについては、どの色が最も大きくずれるかは画素の座標位置によって異なっている。これらの現象はセンサ中心（＝光学系中心）を基準に左右対称に起きる。
【００１２】
この例のような大きな収差を持った光学系を用いた場合、得られた画像に何も補正をしなければ、主走査方向で中心から画像端に近づくほど大きく縮んだ画像となってしまう。加えて、座標位置によって異なる色ズレが発生してしまう。
【００１３】
上記の特許文献１は、この課題を全てではないが部分的に解決する手法を提案している。ＣＣＤラインセンサの各画素に対し、その画素を間引くか否かを指定する間引きビット情報を、予めメモリに保持しておく。そして、センサからラスタ走査順に入力されてくる画像の各画素データに対して、その間引きビット情報を参照し、間引く画素であればその画素データは出力しないようにする。これにより、縮小光学系の歪みに起因する画素の位置ズレを補正するものである。
【００１４】
しかし、特許文献１の装置は、色収差による色ズレについての課題認識は無く、その解決についての開示も示唆もない。また、入力画素を間引いて解像度を調整するため、センサが本来持っている画素数の一部を捨てることになり、捨てる画素を見込んで高解像度なセンサが必要となる。このようなセンサは一般に高価である。また、間引きビット情報を保有するメモリは、センサの画素数に応じた容量が必要である。センサが高解像度化するほどに大容量のメモリが必要になるため、これも低コスト化の妨げる要因となる。
【００１５】
本発明は、このような課題の認識に鑑みてなされたものである。
【００１６】
本発明の目的は、低コストでシンプルな手段で、光学系の収差の影響を軽減することを実現した画像読取装置を提供することである。本発明の別の目的は、従来以上に小型の画像読取装置を低コストに実現することである。
【課題を解決するための手段】
【００１７】
上記課題を解決する本発明は、アレイセンサと光学系を有する読取手段を用い前記アレイセンサのラスタ走査によって画像を読み取って、画像データに変換して出力する画像読取装置であって、前記光学系の収差の影響を、デジタル変倍処理によって軽減する信号処理手段を有し、該信号処理手段は、前記ラスタ走査の順に各画素の画素値を入力する入力手段；１つのラスタ走査の全画素よりも少ない限られた数の入力画素の画素値を一時記憶するメモリ；前記収差の特性を反映して作成されたデータテーブルを用いて、次の出力画素の座標位置を求める第１処理手段；および、前記第１処理手段で求めた座標位置に対応する、最も近い入力画素の画素値を、前記メモリから読み出し、少なくともこの画素値を用いて次の出力画素の画素値を求める第２処理手段を有することを特徴とするものである。
【００１８】
本発明の別の形態は、アレイセンサと光学系を有する読取系を用い、前記アレイセンサのラスタ走査によって画像を読み取って、画像データに変換して出力する画像読取方法であって、前記光学系の収差の影響を、デジタル変倍処理によって軽減する信号処理を有し、該信号処理は、前記ラスタ走査の順に各画素の画素値を入力するステップ；１つのラスタ走査の全画素を複数の区間に区切って、前記収差の特性を反映して各区間の変倍率を設定し、これに基づいて出力画素の座標位置およびその画素値を演算処理によって求めるステップ；および前記求めた画素値を、前記演算処理の順に出力するステップを有することを特徴とするものである。
【発明の効果】
【００１９】
本発明によれば、低コストでシンプルな手段で、光学系の収差の影響を軽減することを実現した画像読取装置を提供することができる。また、従来以上に小型の画像読取装置を低コストに実現することができる。
【発明を実施するための最良の形態】
【００２０】
本発明の実施の形態の説明に先立ち、その基本概念について説明する。
【００２１】
入力画像に対して、主走査方向のラスタ走査内の座標位置に応じて倍率の異なる変倍処理を、リアルタイムに、ＲＧＢの各色成分毎に別々に行なうものである。
【００２２】
１つのラスタ（主走査方向の１ライン）を複数区間に分けて、異なる区間に対しては異なる変倍率を適用する。１つの区間の中での変倍率は一定値（Ｑ［ｉ＋１］−Ｑ［ｉ］）／（Ｐｒ［ｉ＋１］−Ｐｒ［ｉ］）（Ｒ各画素についての式）とする。これを、ラスタ順に画素値の列がシーケンシャルに入力されてくる画像データに対して、以下の手順でリアルタイムにデジタル変倍処理を行い、光学系の収差（歪曲収差、色収差等）の影響を軽減する。これにより、図１２で説明したような画素ズレを補正すると共に、ＲＧＢ別々に処理を行なうことで、色ズレの補正も同時に行なうことができる。
【００２３】
画像をＲＧＢの色成分毎に個別に読み取る複数のラインセンサの出力をＡ／Ｄ変換して、センサのラスタ走査の順にシーケンシャルに入力されてくる画像データに対し、以下のようなリアルタイム変倍処理を行なう。なお、必要に応じて、Ａ／Ｄ変換した信号にシェーディング処理やガンマ変換処理等の前処理を行なってから、その画像をラスタ走査の順に入力する。
（１）入力ＲＧＢ画素値それぞれを、最大画素ズレ量に対応した限られた画素数分だけ、小容量のＦＩＦＯメモリに一時記憶する。ＦＩＦＯメモリは適宜更新される。
（２）処理開始の初期画素の座標位置および初期位相を、ＲＧＢそれぞれ設定する。
（３）出力すべき画素の座標位置に対応する入力画素の画素値をＦＩＦＯメモリから読み出し、その値と位相とから出力画素値をＲＧＢそれぞれ計算する。計算した結果の画素値は、ＲＧＢ同時に出力する。
（４）画素を出力する毎に、次に出力すべき画素の座標位置に対応する入力画素の座標位置および位相をＲＧＢそれぞれ再計算し、その画素がＦＩＦＯメモリに入力されるのを待つ。
（５）上記（３）（４）の処理を、１ライン分の画素についてラスタ順に繰り返す。
【００２４】
ここで、（４）の処理においては、１つのラスタ走査の全画素を複数の区間に区切って、光学系の収差の特性を反映して各区間の変倍率と画素数を設定する。
【００２５】
出力画素の座標位置は、初期の座標位置からの画素間距離の積算で得る。ここで言う画素間距離とは、リアルタイム変倍処理における、入力画素レートに対する出力画素レートの比であり、別の言い方をすれば、ある区間での局所的な変倍率の逆数である。そして、１ラインの処理の最中に、区間に応じて画素間距離を変更することで、変倍率の異なる変倍処理を行ない、光学系の収差の影響を軽減する。
【００２６】
画素間距離の変更を実現する具体的な手法は次のとおりである。予めデータテーブルとして、［加算値］（ＲＧＢ別に、ある区間の中での局所的な画素間距離）と、［処理画素数］（その区間で処理して出力する画素の数）のデータ組を複数（ｎ組）用意する。ある区間で［処理画素数］を出力して、次の区間に移行する際には、そのデータテーブルにアクセスして、次の区間に対応した画素間距離を変更する。
【００２７】
以上のようなリアルタイム変倍処理とすることで、大容量の一時記憶メモリが不要で、ごく僅かな容量のメモリで高速に処理を行なうことができる。また、入力解像度が大きくなったとしても、［処理画素数］の値を解像度に応じて大きくすれば済むため、必要なメモリが解像度に比例して増大することはない。
【００２８】
＜第一の実施形態＞
以下、本発明を適用した画像読取装置の実施の形態を説明する。なお、本発明を適用する画像読取装置は、単機能のスキャナ装置であってもよいし、スキャナ機能部を組み込んだ所謂マルチファンクションプリンタであってもよい。
【００２９】
図１は、スキャナ装置の全体構成を示すブロック図である。
【００３０】
制御チップ１０は、装置全体を制御するもので、四角線で囲んだ多くの機能ブロックがワンチップ化されている。各機能ブロックは以下のとおりである。
【００３１】
スキャナ画像処理部Ａ１００は、ＡＦＥ（アナログフロントエンド）２４００でＡ／Ｄ変換されたライン単位の画像データを、制御チップ１０の外部から受信する。そして、シェーディング処理など１ラインで行える（副走査方向の複数ライン分の画像データを必要としない）所定の画像処理を行なう。その結果は、ＤＭＡＣ−Ａ（ＤＭＡコントローラ）１１０を用いて画像バス１１００およびメモリ制御部５００を介してメインメモリ２０００へ書き込む。
【００３２】
スキャナ画像処理部Ｂ２００は、ＤＭＡＣ−Ｂ（ＤＭＡコントローラ）２１０を用いてメインメモリ２０００から複数ラインの画像データを読み出し、フィルタ処理など複数ライン分の画像データを必要とする所定の画像処理を行なう。そして再びＤＭＡＣ−Ｂ２１０を用いて画像バス１１００およびメモリ制御部５００を介してメインメモリ２０００へ書き出す。メモリ制御部５００は、画像バス１１００、制御バス１２００とメインメモリ２０００との間の調停やインターフェース制御を行なう。メインメモリ２０００は、ＣＰＵ６００や各処理ブロックが作業用メモリとして使用する。
【００３３】
ＣＰＵ６００は装置全体の制御を司る。モータ制御部７００は、アレイセンサを移動させるモータの制御を行なう。パネルＩ／Ｆ８００は、操作部２５００とのインターフェース制御を行なう。ＳＩＯ（シリアルＩ／Ｏ）９００は、ＲＯＭ２６００とのインターフェース制御を行なう。ＵＳＢデバイスＩ／Ｆ１０００は、スキャナ装置外部のＰＣ１３００との通信処理を行なう。これらの各ブロックは、共通の制御バス１２００に接続されている。
【００３４】
画像読み取り部２３００は、縮小光学系とアレイセンサから構成される読取手段と、これを原稿に対して副走査方向に移動走査させる副走査機構を有する。アレンセンサによる主走査方向のラスタ走査と、副走査機構による副走査との組み合わせによって二次元画像を読み取る。アレイセンサは、ＲＧＢ各色にそれぞれ対応した３つのＣＣＤラインセンサを有している。なお、色成分ＲＧＢ３色をさらに多色化して、例えばシアン、マゼンタ、エメラルドなど、他の色成分を検出するラインセンサを付加するようにしてもよい。
【００３５】
ＡＦＥ（アナログフロントエンド）２４００は、画像読み取り部２３００からラスタ走査の順にアナログ出力される画像データをＡ／Ｄ変換し、ライン単位の画像データとして出力する。この信号は、スキャナ画像処理部Ａ１００に入力される。
【００３６】
操作部２５００は、装置のユーザへ各種情報を出力するＬＣＤと、ユーザが直接操作するキーなどから構成され、パネルＩ／Ｆ８００を介して制御チップ１０との情報の入出力を行なう。ＲＯＭ２６００は、ＣＰＵ６００の動作プログラムを格納している。
【００３７】
図２は、スキャナ画像処理部Ａ１００の処理ブロックの構成を示す。
【００３８】
シェーディング処理ブロック１０１は、センサの画素毎の感度ばらつきを正規化する。シェーディングメモリ１０１ａは、シェーディング処理で必要となる、各画素の明暗シェーディングデータを格納するものである。ガンマ変換処理部１０２は、シェーディング処理を施されて正規化された各画素値の輝度特性を補正する。ガンマテーブルメモリ１０２ａは、ガンマ変換で必要となる、入出力の輝度特性の対応データを格納するものである。
【００３９】
歪み補正処理部１０３は、光学系の収差の影響を、デジタル変倍処理によって軽減するための信号処理手段であり、本発明の特徴となる処理を行なう。詳細については後述する。
【００４０】
図３は、スキャナ画像処理部Ｂ２００の処理ブロックの構成を示す。
【００４１】
文字判定処理部２０３は、ある画素を中心としたｎ×ｎ画素の画素値を調べ、その画素が文字画像の一部であるか否かを判定する。この処理のためには、副走査方向ｎ−１ライン分の画像データが必要なので、その分のデータをラインバッファメモリ２０３ａに格納し、随時そこから画像データを読み出すことにより処理を行なう。
【００４２】
フィルタ処理部２０４は、ある画素を中心としたｍ×ｍ画素の画素値から、フィルタ演算によりエッジ強調処理やモアレ抑制処理を行なう。文字判定処理部２０３と同様に副走査方向ｍ−１ライン分の画像データが必要なので、その分のデータをラインバッファメモリ２０４ａに格納し、随時そこから画像データを読み出すことにより処理を行なう。
【００４３】
図４は、図２の歪み補正処理部１０３の処理ブロックの構成を示す。
【００４４】
歪み補正処理部１０３にはγ変換処理部１０２から、ＲＧＢ別の信号線（入力手段に対応）を経てＲＧＢの画像信号がシーケンシャルに入力される。
【００４５】
ＦＩＦＯメモリ１０３１は、ラスタ走査の順で入力されてくるＲＧＢ画像データを数画素分、一時的に記憶する。これは入力の順に、１つのラスタ走査の全画素よりも少ない限られた数の入力画素の画素値を一時記憶するメモリである。
【００４６】
部分倍率テーブルメモリ１０３２は、入力画素に対する出力画素のＲＧＢそれぞれの部分倍率を示す情報（後述する［加算値］）と、その倍率で処理して出力する画素数（後述する［処理画素数］）のデータ組をｎ組分保持する。このデータテーブルは、縮小光学系が有する収差の特性を反映して作成されたものである。部分倍率テーブルメモリ１０３２は、ＣＰＵバスとつながっており、テーブルの内容はＣＰＵ６００から読み書きすることができる。
【００４７】
画素位置積算器１０３３は、部分倍率テーブルメモリ１０３２から読み込んだデータをもとに、出力すべきＲＧＢ画素の入力に対する座標位置をサブピクセル単位で積算する。
【００４８】
ＦＩＦＯアクセス制御部１０３４は、画素位置積算器１０３３の結果に従って、ＦＩＦＯメモリ１０３１の更新、および、ＦＩＦＯメモリ１０３１からどの座標位置のＲＧＢ画素値データを読み出すかを制御する。
【００４９】
補間演算処理部１０３５は、ＦＩＦＯアクセス制御部１０３４の制御に従ってＦＩＦＯメモリ１０３１から読み出されたＲＧＢ各２画素ずつの画素値と、画素位置積算器１０３３で得られる座標位置とから、出力すべき画素値をＲＧＢそれぞれ計算する。計算した結果は、ＲＧＢ別の信号線（出力手段に対応）を経て、画像信号としてＤＭＡ−Ａ１１０にシーケンシャルに出力する。
【００５０】
以下、光学系の収差の影響をリアルタイムのデジタル変倍処理によって軽減する信号処理の手法について詳細に説明する。
【００５１】
図５（ａ）は、主走査方向にＡインチ刻みに点Ｐ［ｉ］（ｉ＝−１０．．１０、ｉは整数）の描かれた原稿である。図５（ｂ）は、その原稿を収差ゼロの歪のない理想光学系で解像度ｄ（ｄｐｉ）で読み取ったときの読み取り画像である。図５（ｃ）は、歪みある現実的な光学系で解像度ｄ（ｄｐｉ）で読み取ったときの読み取り画像である。
【００５２】
Ｐ［０］がセンサの中心に位置するように原稿位置を調整するものとする。収差による光学歪みはセンサ中心を軸に左右対称に現れるため、以降はｉ≧０の半分について考える。実際には１ラインに対して前半と後半で対称的な処理を行なう。
【００５３】
図５（ｂ）において、Ｐ［ｉ］の各点の読み取り画像がｑ［ｉ］になるものとし、ｑ［ｉ］とｑ［ｉ＋１］の間の距離（ピクセル数）をＱ［ｉ］とする。図５（ｂ）は歪みがない理想形なので、Ｑ［ｉ］はｉによらず一定値ｄＡである。
【００５４】
同様に、図５（ｃ）において、Ｐ［ｉ］の各点の読み取り画像がＲ（赤）画像においてｑｒ［ｉ］になるとし、ｑｒ［ｉ］とｑｒ［ｉ＋１］との間の距離（ピクセル数）をＱｒ［ｉ］とする。
【００５５】
ｑ［０］とｑ［ｉ］の間の距離はｉ＊Ｑ［０］ピクセルとなり、また先の図１２のグラフの画素ズレ量データからｑｒ［０］とｑｒ［ｉ］の間の距離はｉ＊Ｑ［０］−Ｄｒ［ｉ］ピクセルとなる。よって、Ｑｒ［ｉ］＝（ｉ＋１）Ｑ［０］−Ｄｒ［ｉ＋１］−（ｉＱ［０］−Ｄｒ［ｉ］）＝Ｑ［０］＋Ｄｒ［ｉ］−Ｄｒ［ｉ＋１］となる。
【００５６】
光学系の歪みにより、ｑ［ｉ］とｑ［ｉ＋１］の間の区間の距離（ピクセル数）は、歪みのない場合のＱ［０］に対して、Ｑｒ［ｉ］／Ｑ［０］＝１＋（Ｄｒ［ｉ］−Ｄｒ［ｉ＋１］）／Ｑ［０］倍に縮小される。この値を１／Ｔｒ［ｉ］とする。Ｇ（緑）・Ｂ（青）についても同様に１／Ｔｇ［ｉ］、１／Ｔｂ［ｉ］とする。
【００５７】
歪みのある光学系でのＰｌの読み取り画像のＲＧＢについての画素ズレ量（初期画素ズレ量）は、次のように近似できる。ここで、原稿面上でセンサ左端に対応する点をＰｌとし、Ｐ［０］とＰｌの間の距離Ｌ（インチ）がＬ＝（Ｋ＋ｃ）Ａ（Ｋは整数、０≦ｃ＜１）を満たすとする。
ＤＬｒ＝（１−ｃ）Ｄｒ［Ｋ］＋ｃＤｒ［Ｋ＋１］
ＤＬｇ＝（１−ｃ）Ｄｇ［Ｋ］＋ｃＤｇ［Ｋ＋１］
ＤＬｂ＝（１−ｃ）Ｄｂ［Ｋ］＋ｃＤｂ［Ｋ＋１］
歪みのある光学系で読み取った画像データでは、センサ端からこの数だけ無効画素（原稿の読み取り範囲外に相当する）が入力されることになる。
【００５８】
以上をまとめると、このような読み取り画像を補正するためには、Ｒ（赤）については図６に示すフローで処理を行なう。詳細は次のとおりである。
ｓ１．  入力開始から、初期画素ズレ量分の入力画素ＤＬｒ個（ピクセル）を読み飛ばす。
ｓ２．  入力順に入力画素に対し、Ｔｒ［Ｋ＋１］倍の変倍処理を行ない、ｃ＊ｄＡ（一定値）画素分をくり返し出力する。
ｓ３．  ｊ＝Ｋ
ｓ４．  入力順に入力画素に対し、Ｔｒ［ｊ］倍の変倍処理を行ない、ｄＡ（一定値）画素分をくり返し出力する。
ｓ５．  ｊ＝０ならばｓ７へ
ｓ６．  ｊ＝ｊ−１、ｓ４へ
ｓ７．  入力順に入力画素に対し、Ｔｒ［ｊ］倍の変倍処理を行ない、ｄＡ（一定値）画素分をくり返し出力する。
ｓ８．  ｊ＝ｊ＋１
ｓ９．  ｊ＝Ｋ＋１ならばｓ１０へ、でなければｓ７へ
ｓ１０．  入力順に入力画素に対し、Ｔｒ［Ｋ＋１］倍の変倍処理を行ない、ｃ＊ｄＡ（一定値）画素分をくり返し出力する。
【００５９】
以上のＲ（赤）の処理と並行して、Ｇ（緑）・Ｂ（青）の入力画素についても、初期画素ズレ量ＤＬｇ・ＤＬｂ、変倍率Ｔｇ［ｉ］・Ｔｂ［ｉ］での変倍処理を行なう。こうして、図５（ｃ）の画像データから、図５（ｄ）のような歪が補正された画像データを生成することができる。
【００６０】
以上の手順をまとめると、まず、ラスタ順に入力されてくる画像に対し、初期画素ズレ量分の画素数だけ読み飛ばす。その後、ＲＧＢそれぞれ異なる倍率で変倍処理を行ない、且つその変倍率を所定の画素数出力する毎に変更する。こうして、画素ズレおよび色ズレを補正した画像を得ることができる。
【００６１】
縮小光学系の光学歪みは、レンズ中心を軸として左右対称に現れる。そこで、図６に示すフローチャートの処理では、ｄＡ（一定値）画素出力する毎に変倍率が、Ｔｒ［Ｋ＋１］→Ｔｒ［Ｋ］→Ｔｒ［Ｋ−１］→…→Ｔｒ［１］→Ｔｒ［０］→Ｔｒ［０］→Ｔｒ［１］→Ｔｒ［２］→…→Ｔｒ［Ｋ］→Ｔｒ［Ｋ＋１］と変化させる。すなわち、１ライン（１つのラスタ出力）に対して前半と後半で対称的な処理を行なう。
【００６２】
図７は、デジタル変倍処理の補間処理の概念を示す図である。ここでは、左端の初期画素とその次の画素を出力する処理を説明する。
【００６３】
Ｒ（赤）の初期画素ズレ量ＤＬｒ（ピクセル）の計算法は、先に説明したとおりである。ここで、この値の整数部をｓｋｉｐ＿ｐｉｘ＿ｒ、小数部をｆｉｒｓｔ＿ｐｈａｓｅ＿ｒとする。ｆｉｒｓｔ＿ｐｈａｓｅ＿ｒは位相を意味する。
【００６４】
歪みのある光学系では、原稿左端Ｐｌに対応する画素は、入力開始からｓｋｉｐ＿ｐｉｘ＿ｒ番目のＲ入力画素と、ｓｋｉｐ＿ｐｉｘ＿ｒ＋１番目のＲ入力画素の間の座標位置にある。正確には、それら画素の間をｆｉｒｓｔ＿ｐｈａｓｅ＿ｒ：１−ｆｉｒｓｔ＿ｐｈａｓｅ＿ｒの比率に分割した座標位置にある。
【００６５】
従って、Ｐｌ（小文字のエル）に対応する出力画素値Ｒｌは、ｓｋｉｐ＿ｐｉｘ＿ｒ番目の入力画素のＲ画素値Ｒ０と、ｓｋｉｐ＿ｐｉｘ＿ｒ＋１番目の入力画素のＲ画素値Ｒ１の２つの入力画素の値から、次のような補間演算処理で求めることができる。
Ｒｌ（小文字のエル）＝Ｒ０＊（１−ｆｉｒｓｔ＿ｐｈａｓｅ＿ｒ）＋Ｒ１＊ｆｉｒｓｔ＿ｐｈａｓｅ＿ｒ
次に、Ｐｌの次に出力すべき画素値を求める。Ｐｌ付近でのＲ変倍率はＴｒ［Ｋ＋１］である。原稿上でＰｌの１／ｄインチ（スキャナ解像度の逆数＝光学歪みがない場合の読み取り画像における１画素幅に相当）右にある点に対応する入力画素の座標位置は、入力開始からＤＬｒ＋（１／Ｔｒ［Ｋ＋１］）（ピクセル）の場所にある。ここでは、ＤＬｒに加算する値をＲ画素間距離と呼ぶ。この入力画素の座標位置について、整数部と小数部に分離し、整数部の値から出力画素の座標位置に対応する入力画素の座標位置を、小数部の値から位相を得ることができるので、上記Ｒｌと同様の補間演算処理で出力画素値を計算することができる。これはＰｌ付近のＲ画素について、Ｔｒ［Ｋ＋１］の倍率で変倍処理を行なうのと等価である。これを繰り返してｃ＊ｄＡ（一定値）画素分だけ順次出力する。
【００６６】
その後、変倍率をＴｒ［Ｋ＋１］からＴｒ［Ｋ］に変更して、ｄＡ（一定値）分の画素を出力するまで処理を繰り返す（図６参照）。以上の手順により、１ラスタ分の入力画素列から、歪みが補正された１ラスタ分の出力画素列を得ることができる。別の色成分Ｇ・Ｂについても同様の処理を並行して行ない、ＲＧＢすべての色成分について出力画像を得る。
【００６７】
次に、歪み補正処理部１０３での詳細な処理手順について説明する。
【００６８】
予め、初期画素ズレ量ＤＬｒ、ＤＬｇ、ＤＬｂと、サンプル点Ｐ［ｉ］とＰ［ｉ＋１］の間の区間の部分変倍率Ｔｒ［ｉ］、Ｔｇ［ｉ］、Ｔｂ［ｉ］  （ｉ＝０．．Ｋ）の値から、以下（１）（２）（３）の値を計算しておく。
（１）  ｓｋｉｐ＿ｐｉｘ＿ｒ、ｓｋｉｐ＿ｐｉｘ＿ｇ、ｓｋｉｐ＿ｐｉｘ＿ｂ：それぞれ  ＤＬｒ、ＤＬｇ、ＤＬｂの整数部
（２）  ｆｉｒｓｔ＿ｐｈａｓｅ＿ｒ、ｆｉｒｓｔ＿ｐｈａｓｅ＿ｇ、ｆｉｒｓｔ＿ｐｈａｓｅ＿ｂ：それぞれＤＬｒ、ＤＬｇ、ＤＬｂの小数部
（３）  ＲＧＢの、区間ごとの画素間距離ｘｒ［ｊ］、ｘｇ［ｊ］、ｘｂ［ｊ］  （０≦ｊ≦２Ｋ−１）
ただしｊ≦Ｋのときｘｒ［ｊ］＝＝１／Ｔｒ［Ｋ−ｊ］、ｘｇ［ｊ］＝１／Ｔｇ［Ｋ−ｊ］、ｘｂ［ｊ］＝１／Ｔｂ［Ｋ−ｊ］
Ｋ＜ｊ＜２Ｋのときｘｒ［ｊ］＝１／Ｔｒ［Ｋ＋１＋ｊ］、ｘｇ［ｊ］＝１／Ｔｇ［Ｋ＋１＋ｊ］、ｘｂ［ｊ］＝１／Ｔｂ［Ｋ＋１＋ｊ］
ここで、（３）のｉ→ｊの変換は、Ｔｒ［ｉ］がセンサ中心で左右対称であることを利用したものである。画像入力は原稿左端から行われるため、各ｉに対応する区間をセンサ左端からの順番に変換する。
【００６９】
（１）（２）については、画素位置積算器１０３３内のレジスタに設定しておく。
【００７０】
（３）については、［処理画素数］と［Ｒ加算値］、［Ｇ加算値］、［Ｂ加算値］のデータ組を（２Ｋ＋１組）複数有するデータテーブルを作成して、下記のような形式で、部分倍率テーブルメモリ１０３２に書き込んでおく。このデータテーブルは、光学系が有する収差の特性を反映したものである。
【００７１】
【表１】
2010068124.tif 000003
【００７２】
ここで、区間出力画素数（＝［処理画素数］）とは、各区間ごとの出力画素数ｏｕｔｐｉｘ［ｊ］であり、ある区間ｊで処理して出力する画素数を意味する。Ｒ画素間距離（＝［Ｒ加算値］）、Ｇ画素間距離（＝［Ｇ加算値］）、Ｂ画素間距離（＝［Ｂ加算値］）は、各色における、ある区間ｊで処理する際の所定の加算値を意味する。
【００７３】
ｏｕｔｐｉｘ［ｊ］は、各区間における出力画素数で、図５の例ではｊ＝０  ｏｒ  ２Ｋ−１のときｏｕｔｐｉｘ［ｊ］＝ｃ＊ｄＡ、１≦ｊ≦２Ｋ−２のときｏｕｔｐｉｘ［ｊ］＝ｄＡとなる。また、ｊ＝２Ｋでｏｕｔｐｉｘ［ｊ］＝０と設定することで、１ラスタ分の処理終了を判断するフラグとする。
【００７４】
なお、より一層の精度向上を求めるのであれば、１≦ｊ≦２Ｋ−２のときｏｕｔｐｉｘ［ｊ］＝ｄＡ（一定値）とするのではなく、ｄＡの値をｊに応じた非一定値として設定するようにしてもよい。こうすることで、より正確に光学系の収差を反映したデータテーブルとすることができる。
【００７５】
図８は、歪み補正処理部１０３の処理手順を示すフローチャートである。また、図９は、ＦＩＦＯメモリを用いた処理を概念を示す図である。以下に処理手順を説明する。
ｓ１．  ラスタ画像の入力開始。シェーディング処理およびγ変換処理を施されたラスタ画像データを、１画素ずつＲＧＢ同時にＦＩＦＯメモリ１０３１へ入力する。なお、ＦＩＦＯの深さ（ＦＩＦＯ内に画素値を格納できる画素数）をｄｅｐｔｈとする。図９の例では５である。
ｓ２．  画素位置積算器１０３３の内部変数を初期化する。ｐｌａｃｅ＿ｒ（ｇ，ｂ）は入力に対する出力画素の座標位置の積算値、ｐｈａｓｅ＿ｒ（ｇ，ｂ）は入力画素に対する出力画素の位相である。ａｒｅａ＿ｃｏｕｎｔｅｒはデータテーブルにアクセスする際のテーブルインデックス値で、倍率が変わる際に増加する。
ｐｌａｃｅ＿ｒ＝ｓｋｉｐ＿ｐｉｘ＿ｒ、ｐｈａｓｅ＿ｒ＝ｆｉｒｓｔ＿ｐｈａｓｅ＿ｒ、
ｐｌａｃｅ＿ｇ＝ｓｋｉｐ＿ｐｉｘ＿ｇ、ｐｈａｓｅ＿ｇ＝ｆｉｒｓｔ＿ｐｈａｓｅ＿ｇ、
ｐｌａｃｅ＿ｂ＝ｓｋｉｐ＿ｐｉｘ＿ｂ、ｐｈａｓｅ＿ｂ＝ｆｉｒｓｔ＿ｐｈａｓｅ＿ｂ、
ａｒｅａ＿ｃｏｕｎｔｅｒ＝０
ｓ３．  画素位置積算器１０３３は、最初に、部分倍率テーブルメモリ１０３２が保持するデータテーブルのエントリ０（＝ａｒｅａ＿ｃｏｕｎｔｅｒ）をアクセスする。そして、ＲＧＢ画素間距離ｘｒ［０］、ｘｇ［０］、ｘｂ［０］と、区間０の出力画素数ｏｕｔｐｉｘ［０］を画素位置積算器１０３３の内部変数ｘｒ、ｘｇ、ｘｂ、ｏｕｔ＿ｃｏｕｎｔｅｒにそれぞれセットする。ａｒｅａ＿ｃｏｕｎｔｅｒ＝ａｒｅａ＿ｃｏｕｎｔｅｒ＋１
ｓ４．  ｓｋｉｐ＿ｐｉｘ＿ｒ、ｓｋｉｐ＿ｐｉｘ＿ｇ、ｓｋｉｐ＿ｐｉｘ＿ｂのうち最大値をｓｋｉｐ＿ｐｉｘとする。
ｓ５．  ｓｋｉｐ＿ｐｉｘ＋１画素目までの入力画素データをＦＩＦＯメモリ１０３１へ入力する。このとき、ＦＩＦＯ内に存在する最新画素が、ラスタ最初からの何番目の入力画素であるかを表す内部変数ｆｉｆｏ＿ｔｏｐ＝ｓｋｉｐ＿ｐｉｘ＋１となる。
ｓ６．  画素位置積算器１０３３内で、次の内部変数を更新する。
ｆｉｆｏ＿ｒ０＝ｄｅｐｔｈ−（ｆｉｆｏ＿ｔｏｐ−ｐｌａｃｅ＿ｒ）、ｆｉｆｏ＿ｒ１＝ｆｉｆｏ＿ｒ０＋１
ｆｉｆｏ＿ｇ０＝ｄｅｐｔｈ−（ｆｉｆｏ＿ｔｏｐ−ｐｌａｃｅ＿ｇ）、ｆｉｆｏ＿ｇ１＝ｆｉｆｏ＿ｇ０＋１
ｆｉｆｏ＿ｂ０＝ｄｅｐｔｈ−（ｆｉｆｏ＿ｔｏｐ−ｐｌａｃｅ＿ｂ）、ｆｉｆｏ＿ｂ１＝ｆｉｆｏ＿ｂ０＋１  の各値を計算する。
【００７６】
これらの値は、ＲＧＢそれぞれにおいて、出力すべき画素に対応する２つの入力画素が格納されているＦＩＦＯのセル位置を示す。
ｓ７．  ＦＩＦＯアクセス制御部１０３４は、ＦＩＦＯメモリ１０３１内の、ｆｉｆｏ＿ｒ０とｆｉｆｏ＿ｒ１で指定されるセルにアクセスし、該当するＲ画素値を読み出してそれぞれ変数ｉｎｒ０、ｉｎｒ１に格納する。
【００７７】
図９の例では、ｄｅｐｔｈ＝５、ｐｌａｃｅ＿ｒ＝ｓｋｉｐ＿ｐｉｘ＿ｒ＝ｓｋｉｐ＿ｐｉｘ＝９、ｆｉｆｏ＿ｂｏｔｔｏｍ＝１０、ｆｉｆｏ＿ｒ０＝４、ｆｉｆｏ＿ｒ１＝５である。すなわち、ＦＩＦＯメモリ１０３１の５つのセルの内、セル番号（４）と（５）に格納されているＲ画素値をそれぞれ読み出して、ｉｎｒ０、ｉｎｒ１に格納する。Ｇ（ｆｉｆｏ＿ｇ０，１）・Ｂ（ｆｉｆｏ＿ｂ０，１）についても同様に、ｉｎｇ０、ｉｎｇ１、ｉｎｂ０、ｉｎｂ１に格納する。
ｓ８．  補間演算処理部１０３５は、画素位置積算器１０３３から出力画素の位相ｐｈａｓｅ＿ｒ、ｐｈａｓｅ＿ｇ、ｐｈａｓｅ＿ｂの値を受け取り、次のように補間演算する。そしてそれら出力値ｏｕｔｒ、ｏｕｔｇ、ｏｕｔｂは、ＤＭＡＣ−Ａ１１０へＲＧＢ同時に出力する。
ｏｕｔｒ＝ｉｎｒ０＊（１−ｐｈａｓｅ＿ｒ）＋ｉｎｒ１＊ｐｈａｓｅ＿ｒ
ｏｕｔｇ＝ｉｎｇ０＊（１−ｐｈａｓｅ＿ｇ）＋ｉｎｇ１＊ｐｈａｓｅ＿ｇ
ｏｕｔｂ＝ｉｎｂ０＊（１−ｐｈａｓｅ＿ｂ）＋ｉｎｂ１＊ｐｈａｓｅ＿ｂ
ｓ９．  ｏｕｔ＿ｃｏｕｎｔｅｒ＝ｏｕｔ＿ｃｏｕｎｔｅｒ−１
ｓ１０．  ｏｕｔ＿ｃｏｕｎｔｅｒ＝０ならばｓ１１へ、でなければｓ１３へ
ｓ１１．  画素位置積算器１０３３は、部分倍率テーブルメモリ１０３２のａｒｅａ＿ｃｏｕｎｔｅｒエントリをアクセスして、次の値を読み出し、各内部変数にセットする。
ｘｒ＝ｘｒ［ａｒｅａ＿ｃｏｕｎｔｅｒ］、ｘｇ＝ｘｇ［ａｒｅａ＿ｃｏｕｎｔｅｒ］、ｘｂ＝ｘｂ［ａｒｅａ＿ｃｏｕｎｔｅｒ］、
ｏｕｔ＿ｃｏｕｎｔｅｒ＝ｏｕｔｐｉｘ［ａｒｅａ＿ｃｏｕｎｔｅｒ］
ｓ１２．  ｏｕｔ＿ｃｏｕｎｔｅｒ＝０なら、このラスタの処理は終了。でなければｓ１３へ
ｓ１３．  ａｒｅａ＿ｃｏｕｎｔｅｒ＝ａｒｅａ＿ｃｏｕｎｔｅｒ＋１
ｓ１４．  画素位置積算器１０３３は、次の出力画素のための下記各内部変数の値を更新する。
ｐｌａｃｅ＿ｒ＝ｐｌａｃｅ＿ｒ＋ｐｈａｓｅ＿ｒ＋ｘｒの整数部、ｐｈａｓｅ＿ｒ＝ｐｌａｃｅ＿ｒ＋ｐｈａｓｅ＿ｒ＋ｘｒの小数部
ｐｌａｃｅ＿ｇ＝ｐｌａｃｅ＿ｇ＋ｐｈａｓｅ＿ｇ＋ｘｇの整数部、ｐｈａｓｅ＿ｇ＝ｐｌａｃｅ＿ｇ＋ｐｈａｓｅ＿ｇ＋ｘｇの小数部
ｐｌａｃｅ＿ｂ＝ｐｌａｃｅ＿ｂ＋ｐｈａｓｅ＿ｂ＋ｘｂの整数部、ｐｈａｓｅ＿ｂ＝ｐｌａｃｅ＿ｂ＋ｐｈａｓｅ＿ｂ＋ｘｂの小数部
ｓ１５．  「ｐｌａｃｅ＿ｒ＜ｆｉｆｏ＿ｔｏｐかつｐｌａｃｅ＿ｇ＜ｆｉｆｏ＿ｔｏｐかつｐｌａｃｅ＿ｂ＜ｆｉｆｏ＿ｔｏｐ」ならばｓ６へ、でなければｓ１６へ
ｓ１６．  ＦＩＦＯメモリ１０３１に次の画素を入力する。ｆｉｆｏ＿ｔｏｐ＝ｆｉｆｏ＿ｔｏｐ＋１    ｓ６へ
図１０は、上記処理フローでの処理途中の状態を示す図である。
【００７８】
ｓｋｉｐ＿ｐｉｘ＿ｒ（ｇ，ｂ）とｆｉｒｓｔ＿ｐｈａｓｅ＿ｒ（ｇ，ｂ）で表される、ＲＧＢそれぞれの出力画素の初期座標位置０に対し、入力画素間距離に対する出力ＲＧＢ画素間距離ｘｒ、ｘｇ、ｘｂの値を積算していく。これにより、出力画素の座標位置に対応する入力画素の座標位置（ＦＩＦＯメモリのセル番号）を取得する。そして、その入力画素の画素値（ＦＩＦＯメモリの記憶値）に基づいて補間演算で出力画素の画素値を求める。
【００７９】
ここで、予め設定した画素数（ｏｕｔｐｉｘ［０，１，…］）を出力する毎に、ｘｒ、ｘｇ、ｘｂの値を変更することで、図５に示す各区間について倍率の異なる変倍処理を行っている。
【００８０】
さらに、入力画素をＲＧＢそれぞれに対応したＦＩＦＯメモリに一時蓄え、積算したＲＧＢ画素の座標位置に応じて、ＦＩＦＯメモリの読み出しセル位置を制御して、並行処理を行なっている。これにより、位置ズレと同時に色ズレの補正も行っている。
【００８１】
なお、上記では、次の出力画素の座標位置に最も近傍した座標位置とその次に近傍した座標位置の２つの入力画素の画素値から、位相による重み付けを行った補間処理によって次の出力画素の画素値を求めている。別の方法として、位相を見ることで２画素の内どちらが近いか、すなわち最近傍画素を選択して、その画素の値をそのまま出力するようにしてもよい。すなわち、次の出力画素の座標位置に最も近傍した座標位置の入力画素（最近傍画素）の画素値を、次の出力画素の画素値とするようにしてもよい。これらの方法に共通して言えることは、次の出力画素の座標位置に対応した、最近傍の座標位置の入力画素の画素値を少なくとも用いて、次の出力画素の画素値を求めるということである。
【００８２】
以上説明した実施の形態の形態の画像読取装置によれば、極めてシンプルな構成で、大容量のメモリを用いずとも、光学系の収差の影響である位置ズレおよび色ズレを軽減することができる。また、従来以上に小型の画像読取装置を低コストに実現することができる。
【図面の簡単な説明】
【００８３】
【図１】スキャナ装置の全体ブロック図
【図２】スキャナ画像処理部Ａのブロック図
【図３】スキャナ画像処理部Ｂのブロック図
【図４】歪み補正処理部の内部構成を示すブロック図
【図５】歪み補正変倍処理の内容を表す図
【図６】歪み補正変倍処理の手順を示すフローチャート図
【図７】歪み補正変倍処理において、１番目と２番目の画素を出力する際の処理内容を示す図
【図８】歪み補正変倍処理の詳細な手順を示すフローチャート
【図９】ＦＩＦＯメモリを用いた処理を概念を示す図
【図１０】区間毎に異なる変倍率の処理を概念示す図
【図１１】縮小光学系の収差が原因で、画素ズレおよび色ズレが生じることを示す図
【図１２】画素ズレおよび色ズレの程度を示すグラフ図
【符号の説明】
【００８４】
１０３  歪み補正処理部（信号処理手段に対応）
１０４、１０３１  ＦＩＦＯメモリ（メモリに対応）
１０３２  部分倍率テーブルメモリ（データテーブルに対応）
１０３３  画素位置積算器（第１処理手段に対応）
１０３４  ＦＩＦＯアクセス制御部（第２処理手段に対応）
１０３５  補間演算処理部（第２処理手段に対応）
２３００  画像読み取り部（読取手段、読取系に対応）
【図１】
2010068124.tif 000004
【図２】
2010068124.tif 000005
【図３】
2010068124.tif 000006
【図４】
2010068124.tif 000007
【図５】
2010068124.tif 000008
【図６】
2010068124.tif 000009
【図７】
2010068124.tif 000010
【図８】
2010068124.tif 000011
【図９】
2010068124.tif 000012
【図１０】
2010068124.tif 000013
【図１１】
2010068124.tif 000014
【図１２】
2010068124.tif 000015