(57)【要約】【課題】離散的なγ補正データを多次元関数により補間して連続的なγ補正テーブルを作成するときに、傾きの逆転を生じず、かつ滑らかなγ補正テーブルを得る。【解決手段】単調増加もしくは単調減少の離散的なγ補正データを元に、その入力階調値の所定の区間毎に多次元関数を割り当て、離散的なγ補正データを補間することにより連続的なγ補正データを算出する（Ｓ11〜14）。算出された連続的なγ補正データの前記区間毎に傾きの逆転の有無を判定し(Ｓ15,16)、逆転ありと判断された部分を含む区間について、逆転を生じない補間処理により連続的なγ補正データを算出する（Ｓ18）。逆転を生じない補間処理によりγ補正データを算出した区間と、その区間に隣接する区間とが接続する点においてγ補正データが滑らかに接続されるよう、前者以外の区間のγ補正データを、離散的なγ補正データの区間毎に多次元関数を割り当てて補間処理によって連続的なγ補正データを算出し直す（Ｓ19,13,14）。【選択図】  図３
2010219743.tif 000002
【特許請求の範囲】
【請求項１】
離散的なγ補正データを補間することで、連続的なγ補正テーブルを作成するγ補正テーブル作成手段と、そのγ補正テーブルを用いて、画像データの階調特性を変換するγ変換手段とを備えた画像形成装置であって、
前記γ補正テーブル作成手段は、
前記画像データの階調値が取り得る値の範囲で単調増加もしくは単調減少のいずれかである離散的なγ補正データを元に、その離散的なγ補正データの入力階調値の所定の区間毎に多次元関数を割り当て、前記離散的なγ補正データを補間することにより連続的なγ補正データを算出する第１のγ補正データ算出手段と、
その算出された連続的なγ補正データの前記区間毎に傾きの逆転の有無を判定する逆転判定手段と、
傾きの逆転があると判定された部分を含む区間について、前記多次元関数による補間処理に代えて傾きの逆転を生じない補間処理によって連続的なγ補正データを算出する第２のγ補正データ算出手段と、
その第２のγ補正データ算出手段がγ補正データを算出した区間と、その区間に隣接する区間とが接続する点においてγ補正データが滑らかに接続されるよう、前者以外の区間のγ補正データを、離散的なγ補正データの区間毎に多次元関数を割り当てて補間処理によって連続的なγ補正データを算出し直す第３のγ補正データ算出手段と
を有することを特徴とする画像形成装置。
【請求項２】
請求項１に記載された画像形成装置において、
前記第１のγ補正データ算出手段および第３のγ補正データ算出手段は、三次スプライン補間により多次元関数を割り当てることを特徴とする画像形成装置。
【請求項３】
請求項１または２に記載された画像形成装置において、
前記逆転を生じない補間処理は、直線補間処理であることを特徴とする画像形成装置。
【請求項４】
請求項１から３のいずれかに記載された画像形成装置において、
前記第３のγ補正データ算出手段は、前記第２のγ補正データ算出手段がγ補正データを算出した区間と、その区間に隣接する区間とが接続する点において傾きが一致するよう、γ補正データを算出し直すことを特徴とする画像形成装置。
【請求項５】
請求項４に記載された画像形成装置において、
前記接続する点において、前記第２のγ補正データ算出手段がγ補正データを算出した区間の傾きと、その区間に隣接する区間の傾きとの差が所定値より小さいか否かを判定する傾き差判定手段を有し、前記第３のγ補正データ算出手段は、前記傾きの差が所定値以上のときのみγ補正データを算出し直すことを特徴とする画像形成装置。
【請求項６】
画像データの階調特性を変換するためのγ変換に用いるγ補正テーブルを作成する方法であって、
前記画像データの階調値が取り得る値の範囲で単調増加もしくは単調減少のいずれかである離散的なγ補正データを元に、その離散的なγ補正データの入力階調値の所定の区間毎に多次元関数を割り当て、前記離散的なγ補正データを補間することにより連続的なγ補正データを算出する第１のγ補正データ算出ステップと、
その算出された連続的なγ補正データの前記区間毎に傾きの逆転の有無を判定する逆転判定ステップと、
傾きの逆転があると判断された部分を含む区間について、前記多次元関数による補間処理に代えて傾きの逆転を生じない補間処理によって連続的なγ補正データを算出する第２のγ補正データ算出ステップと、
その第２のγ補正データ算出ステップでγ補正データを算出した区間と、その区間に隣接する区間とが接続する点においてγ補正データが滑らかに接続されるよう、前者以外の区間のγ補正データを、離散的なγ補正データの区間毎に多次元関数を割り当てて補間処理によって連続的なγ補正データを算出し直す第３のγ補正データ算出ステップと
を有することを特徴とするγ補正テーブル作成方法。
【請求項７】
画像形成装置のコンピュータに、請求項６に記載されたγ補正テーブル作成方法の各ステップを実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、デジタル複写機、ファクシミリ、プリンタ等の画像形成装置、そのγ補正テーブル作成方法およびプログラムに関する。
【背景技術】
【０００２】
デジタル複写機等の画像形成装置は、原稿をＣＣＤ（電荷結合素子）等のイメージセンサにより読み取り、電気信号に変換した後、レーザープリンタ、インクジェットプリンタ、熱転写プリンタ等の出力装置を用いて印字記録することで複写画像を得るものである。
【０００３】
この場合、読み取った原稿画像に忠実な画像を複製する必要がある。しかし、特に出力装置においては、その記録特性がロットや環境によって異なり、また経時的に変動するため、原稿画像に対して複写画像の忠実性が損なわれるという問題がある。この忠実性は、写真画像のような階調性が要求される画像において特に厳しく、さらにカラー画像を形成する場合にはモノクロ画像よりも高い品質が要求される。
【０００４】
このため、画像形成装置は、γ補正機能を備え、適宜γ補正データを調整できるように構成されている。まず、図６を用いてγ補正の原理について説明する。
【０００５】
図６において、第１象限（１）は、印刷データの出力レベル（ＶOUT）と画像形成装置が形成する画像濃度との関係を表し、プリンタの現状のγ特性を示す。第３象限（３）は、所定の出力レベル（ＶIN）に対して期待される画像濃度を表す。即ち目標γ特性である。この現状のγ特性から目標のγ特性へと変換するのが、第４象限（４）に示すγ補正カーブである。第２象限（２）には恒等変換テーブルが示されている。
【０００６】
ここでγ補正カーブデータは、出力レベルＶIN，ＶOUT が８ビットデータの場合、256段の入力値に対して256値の出力値の対応を規定するテーブルであるから、256バイトのデータ量となる。このデータは、次にγ補正データが調整されるまで、不揮発性のメモリに保存しておく必要がある。このときに必要なメモリ容量を減らすために、図６の第４象限（４）に示すように、丸印でプロットした点（８点ある）のみ記憶すれば、入出力の座標データが１バイトずつなので、合計16バイトのメモリで済むことになる。
【０００７】
ここで、離散的なデータを補間して連続的なテーブルを作成するため、３次スプライン関数のように区分的に３次関数の係数を決めることで、多次元関数の次数を制限して計算の負荷を軽減しつつ、柔軟なカーブに対応する方法が知られている。
【０００８】
しかしながら、このような多次元関数を用いた補間方法においては、滑らかなγ補正カーブを作成できるものの、そのカーブに予期しない振動が生じる場合がある。例えば図７に示すような入、出力階調特性を持つ離散的なγ補正データを補間して連続的なγ補正テーブルを作成した場合、図８に示すようになる。図９は図８のγ補正テーブルの傾き（一次微分）であり、入力階調値129から140の間で負の値を取っている、すなわち逆転していることを示している。つまり、元となる離散点データは、単調に増加しているにもかかわらず、γ補正テーブルでは入力階調値が129から140の間で減少してしまうため、このγ補正テーブルを用いてγ補正処理を行った場合、画像上においても濃度の逆転が生じるという不具合が発生し、所望の出力画像を得ることができないという問題がある。
【０００９】
このような問題に対処した画像形成装置として、多次元関数によって補間データを区間毎に作成し、傾き逆転が発生している区間では、作成されたγ補正データを直線補間データで置き換えることにより、画像上においても濃度の逆転がなく、かつ滑らかなγ補正テーブルを作成し、良好な再生画像を得ることのできる画像形成装置がある（特許文献１参照）。図１０は、このように作成されたγ補正テーブルであり、図１１はその傾きを示すグラフである。
【００１０】
ところが、このγ補正テーブルの場合、逆転が発生している区間では直線補間により作成されたγ補正データと置き換えることで逆転のないγ補正データを得ることが出来るが、多次元関数によって作成された区間と、直線補間により作成された区間との境界部でγ補正テーブルの傾きが一致しないことにより、急激な階調変化をもたらす可能性がある。図１０からも明らかなように、図９のγ補正テーブルでは、特に入力階調値159付近で急激に出力階調値が変化している。
【発明の概要】
【発明が解決しようとする課題】
【００１１】
本発明は、このような問題を解決するためになされたもので、その目的は、離散的なγ補正データを多次元関数により補間して連続的なγ補正テーブルを作成するときに、傾きの逆転を生じず、かつ滑らかなγ補正テーブルが得られるようにすることである。
【課題を解決するための手段】
【００１２】
本発明の画像形成装置は、離散的なγ補正データを補間することで、連続的なγ補正テーブルを作成するγ補正テーブル作成手段と、そのγ補正テーブルを用いて、画像データの階調特性を変換するγ変換手段とを備えた画像形成装置であって、前記γ補正テーブル作成手段は、前記画像データの階調値が取り得る値の範囲で単調増加もしくは単調減少のいずれかである離散的なγ補正データを元に、その離散的なγ補正データの入力階調値の所定の区間毎に多次元関数を割り当て、前記離散的なγ補正データを補間することにより連続的なγ補正データを算出する第１のγ補正データ算出手段と、その算出された連続的なγ補正データの前記区間毎に傾きの逆転の有無を判定する逆転判定手段と、傾きの逆転があると判定された部分を含む区間について、前記多次元関数による補間処理に代えて傾きの逆転を生じない補間処理によって連続的なγ補正データを算出する第２のγ補正データ算出手段と、その第２のγ補正データ算出手段がγ補正データを算出した区間と、その区間に隣接する区間とが接続する点においてγ補正データが滑らかに接続されるよう、前者以外の区間のγ補正データを、離散的なγ補正データの区間毎に多次元関数を割り当てて補間処理によって連続的なγ補正データを算出し直す第３のγ補正データ算出手段とを有することを特徴とする画像形成装置である。
本発明のγ補正テーブル作成方法は、画像データの階調特性を変換するためのγ変換に用いるγ補正テーブルを作成する方法であって、前記画像データの階調値が取り得る値の範囲で単調増加もしくは単調減少のいずれかである離散的なγ補正データを元に、その離散的なγ補正データの入力階調値の所定の区間毎に多次元関数を割り当て、前記離散的なγ補正データを補間することにより連続的なγ補正データを算出する第１のγ補正データ算出ステップと、その算出された連続的なγ補正データの前記区間毎に傾きの逆転の有無を判定する逆転判定ステップと、傾きの逆転があると判断された部分を含む区間について、前記多次元関数による補間処理に代えて傾きの逆転を生じない補間処理によって連続的なγ補正データを算出する第２のγ補正データ算出ステップと、その第２のγ補正データ算出ステップでγ補正データを算出した区間と、その区間に隣接する区間とが接続する点においてγ補正データが滑らかに接続されるよう、前者以外の区間のγ補正データを、離散的なγ補正データの区間毎に多次元関数を割り当てて補間処理によって連続的なγ補正データを算出し直す第３のγ補正データ算出ステップとを有することを特徴とするγ補正テーブル作成方法である。
【００１３】
［作用］
本発明によれば、離散的なγ補正データを多次元関数により補間して連続的なγ補正テーブルを作成するときに、画像データの階調値が取り得る値の範囲で単調増加もしくは単調減少のいずれかである離散的なγ補正データを元に、その離散的なγ補正データの入力階調値の所定の区間毎に多次元関数を割り当て、前記離散的なγ補正データを補間することにより連続的なγ補正データを算出し、その算出された連続的なγ補正データの前記区間毎に傾きの逆転の有無を判定し、傾きの逆転があると判断された部分を含む区間について、前記多次元関数による補間処理に代えて傾きの逆転を生じない補間処理によって連続的なγ補正データを算出し、その第２のγ補正データ算出ステップでγ補正データを算出した区間と、その区間に隣接する区間とが接続する点においてγ補正データが滑らかに接続されるよう、前者以外の区間のγ補正データを、離散的なγ補正データの区間毎に多次元関数を割り当てて補間処理によって連続的なγ補正データを算出し直すことにより、γ補正テーブルを作成する。
【発明の効果】
【００１４】
本発明によれば、離散的なγ補正データを多次元関数により補間して連続的なγ補正テーブルを作成するときに、傾きの逆転を生じず、かつ滑らかなγ補正テーブルが得られる。
【図面の簡単な説明】
【００１５】
【図１】本発明の第１の実施形態の画像形成装置の概略構成を示すブロック図である。
【図２】本発明の第１の実施形態の画像形成装置における画像処理の概略を示すフローチャートである。
【図３】本発明の第１の実施形態の画像形成装置におけるγ補正テーブル算出手段の動作のフローチャートである。
【図４】本発明の第１の実施形態の画像形成装置における三次元スプライン関数の係数及び区間両端の一次微分値、二次微分値を示す図である。
【図５】本発明の第２の実施形態の画像形成装置におけるγ補正テーブル算出手段の動作のフローチャートである。
【図６】γ補正の原理を示す図である。
【図７】離散的なγ補正データを示す図である。
【図８】図７に示すγ補正データを補間して作成したγ補正テーブルを示す図である。
【図９】図８のγ補正テーブルの傾きを示す図である。
【図１０】傾き逆転が発生している区間のγ補正データを直線補間データで置き換えたγ補正テーブルを示す図である。
【図１１】図１０のγ補正テーブルの傾きを示す図である。
【発明を実施するための形態】
【００１６】
以下、本発明の実施形態について図面を参照しながら説明する。
［第１の実施形態］
図１は、本発明の第１の実施形態の画像形成装置の概略構成を示すブロック図である。
この図に示すように、本実施形態の画像形成装置は、スキャナ１と、ＣＰＵ２と、スキャナγ変換回路３と、遅延メモリ４と、拡張ＩＰＵ５と、空間フィルタ処理、色補正処理、ＵＣＲ／ＵＣＡ処理、テストパターン発生処理６と、変倍回路７と、プリンタγ変換回路８と、階調処理回路９と、特定原稿検出回路10と、像域自動分離回路／カラー原稿自動検知回路11と、ＦＩＦＯ12と、Ｉ／Ｆ（インタフェース）回路13と、プリンタ14とを有する。
【００１７】
スキャナ１は、例えば、原稿載置台に載置された原稿を読み取るための機構であり、ＡＤＦ（自動原稿送り装置）が設けられている場合には、１枚ずつ搬出されてくる原稿を読み取る。ＣＰＵ２は、画像形成装置全体の動作制御を司るため、図示されていないＲＯＭに記憶されている制御プログラムを読み出し、当該読み出されたプログラムに基づいて、各回路における動作制御命令を送出する。スキャナγ変換回路３は、スキャナ１により読み取った原稿の画像データのγ特性を変換する。
【００１８】
遅延メモリ４は、入力された画像データと像域信号とのタイミングを調整するために、ＣＰＵ２からの動作指示に基づいて入力された画像データを所定時間遅延させる。拡張ＩＰＵ５は、色変換等の加工処理を行う。画像処理部６は、平滑化処理やエッジ強調処理等を行うため、空間フィルタ、色補正処理、ＵＣＲ／ＵＣＡ処理、テストパターン発生処理等を行う。変倍回路７は、主走査方向における変倍を行う。プリンタγ変換回路８は、プリンタ14におけるγ特性を変換する。階調処理回路９は、離散的な階調変換処理を行う。特定原稿検出回路10は、紙幣や株券などの通常法的に複写することを禁じられている原稿が載置されているか否かを認識する。像域自動分離回路・カラー原稿自動検知回路11は、スキャナ１により読み取られた原稿の画像データから像形成されている領域を分離したり、カラー原稿であるか否かの判定を行ったりする。
【００１９】
ＦＩＦＯ12は、階調処理回路９にて階調変換処理を行われた画像データの出力タイミングを調整するために設けられている。Ｉ／Ｆ回路13は、プリンタ14と接続するためのインタフェースユニットである。
【００２０】
上述した各部は、パラメータを設定して処理条件の変更を行うため、必要に応じてＣＰＵバスと接続されている。ここで、プリンタγ補正カーブデータをプリンタγ変換回路８に設定することで、図６の第３象限（３）に示される目標特性を所望のプリンタγ特性に調整することができる。
【００２１】
スキャナ１で読み取られた原稿の画像データは、スキャナγ変換回路３、遅延メモリ４、拡張ＩＰＵ５、画像処理部６、変倍回路７と、プリンタγ変換回路８、階調処理回路９、特定原稿検出回路10、および像域自動分離回路・カラー原稿自動検知回路11とにより適宜処理を施され、ＦＩＦＯ12、Ｉ／Ｆ回路13を介してプリンタ14に送られて画像が再生される。
【００２２】
図２は、本実施形態における画像処理の概略を示すフローチャートである。この図に示すように、本実施形態の画像処理は、画像入力工程（ステップＳ１）、γ補正工程（ステップＳ２）、画像出力工程（ステップＳ３）からなる。
【００２３】
画像入力工程（ステップＳ１）は画像を一画素ずつ入力する。ここで画像データは、０以上２５５以下の整数値の階調値を一画素毎に持つ。γ補正工程（ステップＳ２）は後述のγ補正テーブル算出手段が算出したγ補正テーブルに基づき、画像入力工程（ステップＳ１）で生成された画像データの階調値を一画素ずつ目標階調値に変換する。画像出力工程（ステップＳ３）はγ補正工程で変換された階調値を一画素ずつ出力する。
【００２４】
γ補正テーブル算出手段について説明する。この手段はＣＰＵ２により実行される演算の機能により実現される手段であり、その結果得られたγ補正テーブルは、プリンタγ変換回路８に設定される。図３は、γ補正テーブル算出手段の動作のフローチャートである。ここでは、多次元関数による補間には三次関数を用い、傾きの逆転しない補間関数として直線補間を用いる。
【００２５】
ステップＳ11では、測定などにより予め算出し、単調増加するよう整形した離散γテーブルを入力する。ここでは、離散γテーブルとして図７に示すテーブルを入力したものとする。ステップＳ12では、三次関数を求めるために両端、即ち、入力値０と入力値255における条件を指定する。ここでは、両端における二次微分値が０である、自然スプラインを用いたものとする。
【００２６】
ステップＳ13では、与えられた離散γテーブルと両端条件を用いて、区間０から区間７に対して、Ｃ++言語プログラムを用いて、区間毎の三次関数を算出する。ここで、区間０は入力値０から31、区間1は入力値32から63、区間２は入力値64から95、区間３は入力値96から126、区間４は入力値127から158、区間５は入力値159から190、区間６は入力値191から222、区間７は入力値223から255とする。
【００２７】
区間毎の三次関数は、各区間の最小値をx0として、
f(x)=a・(x-x0)3+b・(x-x0)2+c・(x-x0)+d  …式［１］
と表現した場合に、xにおける一次微分値f'(x)は、
f'(x)=3a・(x-x0)2+2b・(x-x0)+c  …式［２］
で示され、二次微分値f''(x)は、
f''(x)=6a・(x-x0)+2b  …式［３］
で示される。
【００２８】
そこで、図７の入力値と出力値の関係と、入力値0と入力値255における二次微分値0を条件として与えることで、各区間における各係数が図４のように求まる。また、図４には左右端における一次微分値、二次微分値も併せて示す。
【００２９】
この係数は上記のように条件を与えて連立方程式を解くことで求まる。詳細な方法については例えば、「Ｃによるスプライン関数」、菅野敬祐ほか、東京電機大学出版局（１９９３年）の４４ページから５２ページで解説されている。
【００３０】
ステップＳ14では、ステップＳ13で算出した区間毎の三次関数に基づいて、階調値０から255まで1刻みで入力値と出力値を対応付けたγ補正テーブルを算出する。これにより図８のγ補正テーブルが求まる。
【００３１】
ステップＳ15では、γ補正テーブルの傾きに逆転がないか否か、即ち、入力値iに対応する出力値をf(i)で表すとして、入力値aに対して、f(a)>f(a+1)となる部分がないか否か、区間毎に検査する。なお、この検査はγ補正テーブルで求めても良いし、得られた区間毎の三次関数の一次微分値が、区間内で負となることがあるか否かを計算で求めても良い。
【００３２】
ステップＳ16では、ステップＳ15で逆転した区間があるか否かにより分岐する。逆転した区間がない場合は、ステップＳ17へ進み、γ補正テーブルを出力する。逆転した区間がある場合は、ステップＳ18に進む。このステップでは、逆転した区間を直線補間で置き換える。図８では、入力値127から158の区間４において逆転しているため、この区間を直線補間とする。
【００３３】
次のステップＳ19では、ステップＳ18にて直線補間で置き換えた区間の傾きを、境界条件に利用する。ここでは、直線補間で置き換えた区間４の傾きは、左端が(入力,出力)＝(127,119)、右端が(159,132)であるから、(132-119)÷(159-127)=0.40625である。
【００３４】
そこで、区間０から区間３に対しては左端で二次微分値が０、右端で一次微分値が0.40625、また、区間５から区間７に対しては左端で一次微分値が0.40625、右端で二次微分値が0となるよう条件を与える。
【００３５】
次いで、ステップＳ13で、区間０から区間３に対して区間毎にＣ++言語プログラムを用いて三次関数を求め、また、区間５から区間７に対して区間毎にＣ++言語プログラムを用いて三次関数を求める。以下は同様に、ステップＳ14でγ補正テーブルを算出し、ステップＳ15で逆転を検知して、逆転区間がある場合は再度逆転区間に対して直線補間を行い、三次関数を求める。そして、ステップＳ16で、全ての区間で逆転区間がなくなったと判定した場合、γ補正テーブルを出力する。
【００３６】
なお、三次関数が一意に定まるためには4点が与えられている必要があるが、4点未満の点しか与えられていない以下のような条件では、区間毎にそれぞれ直線補間を行う。例えば、図５のステップＳ16において、区間0から区間7のうち区間2のみで逆転が生じていた場合、区間2に対して直線補間を行い、区間3から区間7までに対しては、それぞれ区間毎に三次関数による補間を行う。そして、区間0から区間1に対しては、三次関数による補間を用いずに区間毎に直線補間を行う。
【００３７】
以上のように、本発明の第1の実施形態の画像形成装置によれば、傾きの逆転を生じず、かつ、滑らかなγ補正テーブルを得ることができる。
【００３８】
［第２の実施形態］
図５は、本発明の第２の実施形態の画像形成装置におけるγ補正テーブル算出手段の動作のフローチャートである。この図において、図３（第１の実施形態）と同じステップには、図３と同じ参照符号を付した。また、本実施形態の画像形成装置の概略構成は第１の実施形態（図１）と同じである。
【００３９】
本実施形態は、図３のステップＳ19の次にステップＳ20を追加したものである。このステップでは、直線補間することとなった区間と、その区間に隣接する区間の傾きを比較し、差が閾値未満である場合は区間毎の三次関数を改めて算出しないことを決定する。
【００４０】
前記閾値を0.5とすると、区間4の傾きは直線補間化されたことで0.40625となり、区間3の右端での一次微分値は図４に示すとおり0.038950である。直線補間する区間4の一次微分値と、それに隣接する三次関数を用いて補間する区間3の右端での一次微分値との差の絶対値0.3673が閾値0.5より小さいため区間0から区間3に対しては三次関数を改めて算出しない。一方、区間4の一次微分値と、区間5の左端での一次微分値との差の絶対値1.38582は閾値0.5より大きいため、区間5から区間7に対しては実施例1と同様に三次関数を改めて算出する。
【００４１】
本発明の第２の実施形態の画像形成装置によれば、階調の急激な変化を感じさせずに、高速にγ補正テーブルを得ることが出来る。
【符号の説明】
【００４２】
２・・・ＣＰＵ、８・・・プリンタγ変換手段、Ｓ１・・・画像入力工程、Ｓ２・・・γ補正工程、Ｓ３・・・画像出力工程。
【先行技術文献】
【特許文献】
【００４３】
【特許文献１】特開２００１−９４７８８号公報
【図１】
2010219743.tif 000003
【図２】
2010219743.tif 000004
【図３】
2010219743.tif 000005
【図４】
2010219743.tif 000006
【図５】
2010219743.tif 000007
【図６】
2010219743.tif 000008
【図７】
2010219743.tif 000009
【図８】
2010219743.tif 000010
【図９】
2010219743.tif 000011
【図１０】
2010219743.tif 000012
【図１１】
2010219743.tif 000013