(57)【特許請求の範囲】
【請求項１】
画素に対応する多値の画像データの値に応じた数のドットの配置を決定するための２値化処理に用いる２値化用パターンを画素の位置に対応付けて配置したパターンテーブルを用い、前記画素の位置と前記画素に対応する多値の画像データの値とに応じて２値化用パターンを選択することで２値化処理を行うデータ処理装置であって、
前記多値の画像データが取り得る値のうちのある１つの値に対応するパターンテーブルは、異なるドット配置を決定するための複数種類の２値化用パターンを含み、
前記ある１つの値に対応するパターンテーブルにおける前記複数種類の２値化用パターンそれぞれの配置は、非周期で、かつ周波数成分が存在する空間周波数領域のうちの半分よりも低周波側の領域に存在する低周波数成分の積分値が前記空間周波数領域のうちの半分よりも高周波側の領域に存在する高周波数成分の積分値よりも少ない特性を有することを特徴とするデータ処理装置。
【請求項２】
前記パターンテーブルを格納したメモリを更に具え、前記メモリに格納されたパターンテーブルを参照することにより前記２値化処理を行うことを特徴とする請求項１に記載のデータ処理装置。
【請求項３】
画素に対応する多値の画像データの値に応じた数のドットの配置を決定するための２値化処理に用いる２値化用パターンを画素の位置に対応付けて配置したパターンテーブルを用い、前記画素の位置と前記画素に対応する多値の画像データの値とに応じて２値化用パターンを選択することで２値化処理を行うデータ処理方法であって、
前記多値の画像データが取り得る値のうちのある１つの値に対応するパターンテーブルは、異なるドット配置を決定するための複数種類の２値化用パターンを含み、
前記ある１つの値に対応するパターンテーブルにおける前記複数種類の２値化用パターンそれぞれの配置は、非周期で、かつ周波数成分が存在する空間周波数領域のうちの半分よりも低周波側の領域に存在する低周波数成分の積分値が前記空間周波数領域のうちの半分よりも高周波側の領域に存在する高周波数成分の積分値よりも少ない特性を有することを特徴とするデータ処理方法。
【請求項４】
情報処理装置に読み取られることにより、該情報処理装置を、画素に対応する多値の画像データの値に応じた数のドットの配置を決定するための２値化処理に用いる２値化用パターンを画素の位置に対応付けて配置したパターンテーブルを用い、前記画素の位置と前記画素に対応する多値の画像データの値とに応じて２値化用パターンを選択することで２値化処理を行うデータ処理装置として機能させるプログラムであって、
前記多値の画像データが取り得る値のうちのある１つの値に対応するパターンテーブルは、異なるドット配置を決定するための複数種類の２値化用パターンを含み、
前記ある１つの値に対応するパターンテーブルにおける前記複数種類の２値化用パターンそれぞれの配置は、非周期で、かつ周波数成分が存在する空間周波数領域のうちの半分よりも低周波側の領域に存在する低周波数成分の積分値が前記空間周波数領域のうちの半分よりも高周波側の領域に存在する高周波数成分の積分値よりも少ない特性を有することを特徴とするプログラム。
【請求項５】
コンピュータを、請求項１または２に記載のデータ処理装置として機能させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、データ処理装置およびデータ処理方法に関し、詳しくは、量子化した画像データを２値化する際に用いるドット配置パターンの展開パターンに関するものである。
【背景技術】
【０００２】
昨今のパーソナルコンピュータ等、情報処理機器の普及にともない、画像を印刷する記録装置もその種類、内容など様々な形態で急速に発展、普及して来ている。種々の記録装置の中でも、特に吐出口からインクを吐出して紙、布、プラスチックシート、ＯＨＰ用シートなどの記録媒体に記録を行うインクジェット記録装置が主流の一つとなりつつある。この種の装置は、ノンインパクト型で低騒音の記録方式であること、高密度かつ高速な記録動作が可能であること、カラー記録にも容易に対応できること、低廉であることなど、極めて優れた特長を有したものである。
【０００３】
このような利点を有したインクジェット記録装置では、インク滴の小液滴化による画像の高解像度化を行い、それによって記録画像のさらなる高画質化などが図られている。しかし、このような高解像度化によって、装置内で処理すべきデータの量は膨大なものとなる。その結果、コンピュータなどのホスト装置とインクジェット記録装置とからなるプリントシステムにおいて、ホスト装置での処理時間やホスト装置からインクジェット記録装置へのデータ転送時間などが長くなるといった問題が生じる。さらに、インクジェット記録装置本体で必要とされるメモリー容量も大きくしなければならず、コストアップをもたらすという問題もある。
【０００４】
ドット配置パターンを用いて２値化を行う方法は、このような問題に対し有効なものである。この方法は、例えば、ホスト装置において画像データを比較的低解像度で量子化するものである。そして、その量子化された多値の画像データをインクジェット記録装置に転送し、記録装置本体で受信した画像データをドット配置パターンによって２値（ドット）データに展開する（特許文献１、特許文献２参照）。このように、ドット配置パターンを用いる２値化は、低解像度のデータ処理、転送を行うものであることから、データ処理、転送の時間が長大なものとなることを防ぐことができ、また、記録装置におけるメモリー容量も小さなものとすることができる。
【０００５】
【特許文献１】特開平０９−０４６５２２号公報
【特許文献２】特許第３６１０２５０号明細書
【発明の開示】
【発明が解決しようとする課題】
【０００６】
しかしながら、従来のドット配置パターンを用いた２値化処理には、特に、そのパターンの展開の仕方において問題を生じることがある。
【０００７】
図１（ａ）および（ｂ）は、従来のドット配置パターンの展開の仕方を説明する図である。
【０００８】
図１（ａ）は、５値に量子化されたデータを２値化するためのドット配置パターンで、５値（レベル０〜４）のうちレベル１に対応したドット配置パターンを示している。同図において、黒塗りの部分がドットを形成する（“１”）２値データ、白の部分がドットを形成しない（“０”）２値データを示している。このように、図に示す例のドット配置パターンは２小画素×２小画素の単位で構成されたものであり、このパターンを用いることにより、最終的に展開される２値画像データに対して縦横とも１／２の解像度で、データ処理、転送などを行うことが可能となる。
【０００９】
以上のようなドット配置パターンの量子化された多値データに対する展開は次のように行われる。例えば、レベル１の多値データが入力すると、その多値データの画素位置に応じて、図１（ａ）に示すパターンＡ、Ｂ、Ｃ、Ｄのいずれかが選択される。そして、その多値データに対応した画素（２小画素×２小画素のサイズ）に対して縦横とも２倍の解像度の位置にドットデータが設定される。
【００１０】
このドット配置パターンを選択すべく、画素位置ごとにパターンを配置ないし展開する方法は、図１（ｂ）に示すように周期的に配置ないし展開する方法（以下、「タイル展開」ともいう）が知られている。その他に、パターンを表すＡ、Ｂ、Ｃ、Ｄについて擬似乱数を発生させ、その乱数に従ってＡ、Ｂ、Ｃ、Ｄを配置ないし展開する方法（以下、「擬似ランダム展開」ともいう）も知られている。そして、２値化処理では、上述のように配置が定められたパターンＡ、Ｂ、Ｃ、Ｄの展開パターンに従い、多値データの画素位置に応じて１つのドット配置パターンを選択する。すなわち、パターンＡ、Ｂ、Ｃ、Ｄは、多値データが５値の場合、レベル０〜４のそれぞれについて定められており、入力する多値データが示すレベルに応じたパターンＡ、Ｂ、Ｃ、Ｄのいずれかが上記のようにして選択される。
【００１１】
しかしながら、上述した従来の展開パターンに従ってドット配置パターンの選択を行う場合、展開パターンにおけるドット配置パターンの選択方法に起因して、形成されるドットに偏りを生じることがある。例えば、図１（ｂ）に示すように周期的な展開パターンを用いた場合、記録装置で実際に印刷される画像にテクスチャーが現れることがある。すなわち、記録ヘッドの吐出精度や装置精度の誤差によってインクの着弾位置がずれると、そのずれがドットの周期的な偏りとして現れ、それがテクスチャーとして視認されることがある。また、擬似乱数に基づく展開パターンを用いる場合は、いくつかのドットが偏って大きなドットの塊を形成し、それが粒状感として視認されることがある。
【００１２】
以上のように、従来のどのような展開パターンでも、パターンＡ、Ｂ、Ｃ、Ｄのような、レベルごとの複数種類のドット配置パターンが適切に分散して配置されていないため、形成されるドットに偏りを生じ記録画像において画質を劣化させることになる。
【００１３】
本発明の目的は、ドット配置パターンを用いた２値化処理において、形成されるドットの偏りを生じないようにすることを可能とするデータ処理装置およびデータ処理方法を提供することである。
【課題を解決するための手段】
【００１４】
そのために本発明では、画素に対応する多値の画像データの値に応じた数のドットの配置を決定するための２値化処理に用いる２値化用パターンを画素の位置に対応付けて配置したパターンテーブルを用い、前記画素の位置と前記画素に対応する多値の画像データの値とに応じて２値化用パターンを選択することで２値化処理を行うデータ処理装置であって、前記多値の画像データが取り得る値のうちのある１つの値に対応するパターンテーブルは、異なるドット配置を決定するための複数種類の２値化用パターンを含み、前記ある１つの値に対応するパターンテーブルにおける前記複数種類の２値化用パターンそれぞれの配置は、非周期で、かつ周波数成分が存在する空間周波数領域のうちの半分よりも低周波側の領域に存在する低周波数成分の積分値が前記空間周波数領域のうちの半分よりも高周波側の領域に存在する高周波数成分の積分値よりも少ない特性を有することを特徴とする。
【００１５】
また、画素に対応する多値の画像データの値に応じた数のドットの配置を決定するための２値化処理に用いる２値化用パターンを画素の位置に対応付けて配置したパターンテーブルを用い、前記画素の位置と前記画素に対応する多値の画像データの値とに応じて２値化用パターンを選択することで２値化処理を行うデータ処理方法であって、前記多値の画像データが取り得る値のうちのある１つの値に対応するパターンテーブルは、異なるドット配置を決定するための複数種類の２値化用パターンを含み、前記ある１つの値に対応するパターンテーブルにおける前記複数種類の２値化用パターンそれぞれの配置は、非周期で、かつ周波数成分が存在する空間周波数領域のうちの半分よりも低周波側の領域に存在する低周波数成分の積分値が前記空間周波数領域のうちの半分よりも高周波側の領域に存在する高周波数成分の積分値よりも少ない特性を有することを特徴とする。
【発明の効果】
【００１６】
以上の構成によれば、画素に対応する多値画像データが取り得る値のうちのある１つの値（ある階調レベル）に対応したパターンテーブル（展開パターン）における複数種類のドット配置パターン（２値化用パターン）それぞれの配置が、非周期で、かつ低周波数成分が高周波数成分より少ない特性を有するようになる。これにより、同じドット配置パターンの配置が非周期となるため、テクスチャーを低減することができ、また、同じドット配置パターンの配置の偏りが少なく高分散となるため、粒状感を低減することができる。
【発明を実施するための最良の形態】
【００１７】
以下、図面を参照して本発明の実施形態を詳細に説明する。
【００１８】
図２は、本発明の一実施形態に係る、主にインクジェット記録装置におけるドット配置パターン（以下、「インデックスパターン」ともいう）を用いた２値化処理の構成を示すブロック図である。
【００１９】
同図に示すように、本実施形態の記録装置２０２は、ホスト装置２０１からＹ、Ｍ、Ｃ、Ｋの各信号について５値の量子化された画像データを受信し、このデータに基づいて記録ヘッドからＹ、Ｍ、Ｃ、Ｋそれぞれのインクを吐出して記録媒体に記録を行う。詳細には、受信バッファ２０２２は、ホスト装置２０１から５値の画像データを受信する。そして、インデックスパターン割り付けモジュール２０２３を実行することによって、受信バッファ２０２２の５値画像データにインデックスパターンを割り付ける。この処理は、以下の各実施形態で後述される展開パターンのテーブルとインデックスパターンを用いて行われる。展開パターンテーブル格納ユニット２０２４はこの展開パターンテーブルを格納し、また、インデックスパターン格納ユニット２０２６はインデックスパターンを格納する。展開バッファ（プリントバッファ）２０２５には、インデックスパターン割付モジュール２０２３によって割り付けられたインデックスパターン（ドット配置パターン）により２値化された結果であるドットデータが展開される。そして、このようにして得られるＹ、Ｍ、Ｃ、Ｋそれぞれのドットデータに基づいて記録ヘッドからインクを吐出して記録を行う。
【００２０】
ＣＰＵ２０２１は、以下で後述される各実施形態のインデックスパターン（ドット配置パターン）の展開処理の実行を制御する。このようにＣＰＵ２０２１を有した記録装置の制御構成はインデックスパターンの展開処理に関してデータ処理装置を構成する。ＣＰＵ２０２１はこれらの処理制御とともに、記録装置２０２における記録ヘッドを用いた記録動作などの公知の制御を実行する。なお、ホスト装置２０１は、本実施形態では、５値の量子化データを作成し、それを記録装置に渡すものであるが、この形態に限られない。ホスト装置２０１において、以下の各実施形態で後述されるインデックスパターンの展開処理まで実行し、その結果の２値データ（ドットデータ）を記録装置２０２に渡すようにしてもよい。この場合、ホスト装置２０１は、インデックスパターンの展開処理に関してデータ処理装置を構成する。
【００２１】
本実施形態では、インデックスパターン（ドット配置パターン）格納ユニット２０２６に、例えば図３に示すような、５値の画像データのレベル０〜レベル４に対応したインデックスパターン（ドット配置パターン）が、４つのグループごとに格納されている。すなわち、本実施形態では、各レベルについて４通りのパターンがあり（レベル０および４は実質的に１通り）、それらが後述される各実施形態の展開パターン（例えば、図１４のような展開パターン）に従って選択される。なお、レベル２のインデックスパターンは、図４（ａ）や図４（ｂ）に示すパターンもあり、それらを用いることもできる。また、グループの数は４つに限られるものでなく、量子化によるレベル数とそれに対応するインデックスパターンの数に応じて定めることができる。さらに、図３に示す例では、各レベルのドット配置パターンは前のレベルを保存した関係になっている。例えば、グループ１において、レベル２のパターンのドット位置は左上と右下であるが、左上のドットは同じグループ１のレベル１の左上のドットをそのまま維持した配置であり、前のレベルを保存している。しかし、本発明を適用する上で、このようなレベル保存の関係が保たれる必要がないことはもちろんである。
【００２２】
以下では、上述の記録システムで用いられるインデックスパターンを選択するべくその配置を定めた展開パターンまたはその製造方法のいくつかの実施形態を説明するが、その前に、展開パターン製造の基本的な方法やそこで用いる斥力計算などの概念を説明する。
【００２３】
（展開パターンの製法）
以下の基本的な展開パターン製造方法の説明では、説明の簡略化のため、インデックスパターン（ドット配置パターン）の配置の仕方を定めた展開パターンを、「プレーン」と称することもある。
【００２４】
本発明の実施形態に係る展開パターンの製造方法は、先ず、レベルごとにインデックスパターンの複数のグループそれぞれの展開パターンのプレーンを規定する。図３に示した例では、レベル０〜４のそれぞれについて、グループ１〜４の４プレーンが規定される。説明の簡略化のため、このプレーンが３つの場合について説明する。先ず、図５（ａ）〜（ｄ）に示すように、プレーンＡ１、Ａ２、Ａ３の３つのプレーンを規定する。そして、同じプレーン内のインデックスパターン間や異なるプレーン間に斥力を作用させるとともに、異なるプレーンのインデックスパターンの重なりを認める。そして、そのような重なり同士の間にも斥力を作用させて、それぞれのプレーン内のインデックスパターンの配置を定める。なお、図５（ａ）〜（ｄ）において、黒または斜線で示す最小の正方形は、インデックスパターンの内容に関わりなく、プレーンにおけるその配置を示すものである。
【００２５】
プレーンにおけるインデックスパターンの配置を定める方法は、大別して、複数のプレーンを同時に定める方法（同時生成）と、プレーンごとに順次定める方法（プレーンごとの生成）の２つの方法のいずれかである。さらに、上記２つの生成方法それぞれについて、具体的にインデックスパターンの配置を定める仕方として、展開パターンにおいて総てのインデックスパターンの配置を予め定めこれらを移動させながら、展開パターン全体で分散性を上げて行く方法がある。以下、この方法を「配置移動法」という。また、その他に、生成される展開パターン全体で分散性を上げながらインデックスパターンを１つずつ配置して行く方法（以下、「順次配置法」）も実行することができる。
【００２６】
配置移動法
配置移動法によるインデックスパターンの配置決定処理の概略は次のとおりである。
【００２７】
例えば配置率が５０％プレーンのインデックスパターンの配置を定める場合、プレーンＡ１、Ａ２、Ａ３それぞれについて、１ビットのデータが“１”であるデータが配置可能位置の５０％に配された初期配置を例えば、誤差拡散法などの２値化処理によって得る。そして、このデータが配置された位置にインデックスパターンを配置することにより、インデックスパターンの初期配置を得る。なお、この２値化の手法を用いてインデックスパターンの初期配置を得るのは、用いる２値化の手法に応じてある程度、初期状態で分散性のよい配置を得ることができるからである。そして、それによって、その後の最終的な配置決定までの演算時間ないし収束時間を短くできるからである。換言すれば、本発明を適用する上で初期配置を得る方法は本質ではなく、例えば、展開パターンのプレーンにおいて、１ビットのデータが“１”であるデータをランダムに配置した初期配置であってもよい。
【００２８】
次に、上記のようにして得たそれぞれのプレーンＡ１、Ａ２、Ａ３の総てのインデックスパターンについて斥力ポテンシャルを計算する。具体的には、
（ｉ）同一プレーン内のインデックスパターン間に距離に応じた斥力を与える。
（ｉｉ）さらに、異なるプレーン間のインデックスパターンにも斥力を与える。
（ｉｉｉ）同一プレーンと異なるプレーン間に異なる斥力を与える。
（ｉｖ）異なるプレーンのインデックスパターンの重なりを認め、インデックスパターンの重なり（２つのインデックスパターンの重なり、３つのインデックスパターンの重なり、…）同士も組み合わせに応じた斥力を与える。
【００２９】
図６は、本実施形態に係る基本斥力ポテンシャルＥ（ｒ）の関数を模式的に示す図である。
【００３０】
同図に示すように、本実施形態で規定する斥力関数は、その斥力が及ぶ範囲をｒ＝１６（１６個分のインデックスパターンが配置される位置）までとする。このような距離とともに減衰するポテンシャルを用いることにより、基本的に、インデックスパターンが接近して配置されるとエネルギーが高い状態、すなわち不安定な状態となり、収束計算の結果、接近した配置はできるだけ選択されないようにすることができる。なお、この斥力の形状は、配置可能位置全体に対するインデックスパターンの配置の割合により決定することがより望ましい。
【００３１】
また、複数のインデックスパターンが重なるインデックスパターン配置を考えるとき、インデックスパターンを配置できる位置（解像度１２００ｄｐｉの場合、１インチ四方に１２００×１２００個の可能位置がある）以上に重ねてインデックスパターンを配置する。このため、各インデックスパターンについて斥力ポテンシャルを計算する際には、インデックスパターンの上にインデックスパターンが重なることを考慮する。このため、ｒ＝０において有限の斥力ポテンシャルを持つように関数を定義する。これにより、インデックスパターンの重なりをも考慮した分散が可能となる。
【００３２】
本実施形態では、同一プレーンのインデックスパターン同士に関してαＥ（ｒ）、異なるプレーン間のインデックスパターン同士に関してβＥ（ｒ）、重なるインデックスパターン同士に関してγｓ（ｎ）Ｅ（ｒ）の斥力ポテンシャルを与えて計算を行う。つまり、あるインデックスパターンが存在することによるポテンシャルは、距離ｒ以内の範囲にある、同プレーンのインデックスパターン、異なるプレーンのインデックスパターン、さらには異なるプレーンの重なるインデックスパターンについての斥力ポテンシャルが加算されたものである。
【００３３】
上記の斥力ポテンシャルにおいて、係数α、β、γは重み付け係数であり、例えば、α＝３、β＝１、γ＝３の値を用いる。このα、β、γの値によってインデックスパターンの分散性が影響を受ける。このα、β、γの値は、例えば、実際には実験を行い、記録される記録画像を参照した最適化により求めることができる。
【００３４】
また、係数ｓ（ｎ）は、重なるインデックスパターンを分散させるためにγに加えてさらに積算する係数である。この係数ｓ（ｎ）は、重なりが多いほどそれらのインデックスパターンをより分散させるべく重なりの数に応じた値とするものである。本願発明者等の実験によれば、次の２つの式いずれかによって求められるｓ（ｎ）を用いることにより、分散に関してよい結果を得ることができる。
【００３５】
【数１】
0004519875.tif 000002
【００３６】
すなわち、ｎを重なりの数とするとき、組合せの数の和をｓ（ｎ）とするものである。詳細には、斥力を計算する注目インデックスパターンに対して重なる（同じプレーンまたは異なるプレーンにおける同じ位置の）インデックスパターンを調べるとともに、注目インデックスパターンから距離ｒに位置するインデックスパターンを調べる。この場合に、注目インデックスパターンおよびそれと同じ位置で重なる他のプレーンのインデックスパターンと、距離ｒにある各プレーンのその位置で同じように重なるインデックスパターンの共通する重なりの数をｎとする。そして、これら２つの位置間の重なったインデックスパターン同士による斥力を考える。
【００３７】
この場合、例えば、ある２つの位置間でプレーンＡ１、プレーンＡ２およびプレーンＡ３にそれぞれ共通にインデックスパターンが存在する例を考えると、ｎ＝３となる。そして、それらの位置間には３つのインデックスパターンの重なりに起因する斥力を作用させる。ここで、３つのインデックスパターンの重なりによる斥力を考えるとき、３つのインデックスパターンの重なりとともに、２つのインデックスパターンの重なり同士や１つのインデックスパターン同士の斥力が多重的に作用すると考える。換言すれば、プレーンＡ３を考えなければ、プレーンＡ１とプレーンＡ２の２つのインデックスパターンの重なりと考えることができ、また、プレーンＡ２を考えなければプレーンＡ１とプレーンＡ３の２つのインデックスパターンの重なりとも考えられる。プレーンＡ１を考えなければプレーンＡ２とプレーンＡ３の重なりと考えられる。このようなインデックスパターンが重なることの多重的な効果を計算するために、重なりの組合せによる斥力を定義し上記のようなｓ（ｎ）を用いる。これによれば、分散性のよいインデックスパターンの配置を得ることができることが実験上確認されている。
【００３８】
なお、本発明の実施形態は、インデックスパターンの配置を定める展開パターンを定める処理である。そして、各レベルつき複数のプレーン（グループ１〜４）のインデックスパターンの配置は相補的である。このため、プレーン間のインデックスパターンの重なりは排除される必要がある。この点から、斥力ポテンシャルの計算では、γを実質的に重なりを排除できる程度の大きな値とする。あるいは、インデックスパターンの移動において処理的に重なりを排除するようにする。
【００３９】
上記のようにして、総てのインデックスパターンの斥力ポテンシャルを合計した総エネルギーが求まると、この総エネルギーを減衰させる処理を行う。
【００４０】
この処理では、総てのインデックスパターンについて、順に距離ｒが４以内のインデックスパターン配置可能位置の中で斥力ポテンシャルが最も下がる位置にインデックスパターンを移す。このような処理を繰り返していくことによって、総てのインデックスパターンの斥力ポテンシャルの合計値である総エネルギーを低下させて行く。すなわち、この総エネルギーが順次減少して行く過程は、インデックスパターンの配置が順次分散性を高める過程、つまりインデックスパターンの配置の低周波数成分が順次少なくなって行く過程である。
【００４１】
そして、総エネルギーの低下率を計算し、それが所定値以下であると判断すると、エネルギー減衰処理を終了する。なお、この所定値は、例えば、実際に印刷を行った結果をもとに、低周波数成分が適切に抑えられた画像を記録できる低下率として求めることができる。最後に、上記のように総エネルギーの低下率が所定値以下となった各プレーンの状態を最終的なインデックスパターンの配置として設定する。
【００４２】
図５（ａ）〜（ｄ）は、上述した斥力ポテンシャルの計算と総エネルギーの減衰処理を模式的に示している。詳しくは、３プレーンＡ１、Ａ２、Ａ３を斜視図で示し、また、特にインデックスパターンの移動を平面図で示す図である。ここで、最小の正方形はインデックスパターンの配置可能位置を示し、３プレーンの重なりにおいて重なる位置がプレーン間で同じ配置可能位置に対応する。
【００４３】
図５（ａ）は、同一プレーンにインデックスパターンが存在する場合にそれらインデックスパターン間の斥力によってポテンシャルが加えられる（増す）ことを説明する図である。図に示す例は、プレーンＡ１の注目位置のインデックスパターンＤｏと同じプレーンで距離ｒ離れた位置にインデックスパターンが１個存在する例である。この場合、α＝３が適用され、インデックスパターンＤｏのポテンシャルとして１×αＥ（ｒ）のポテンシャルが加えられる。
【００４４】
図５（ｂ）は、注目インデックスパターンＤｏとは異なるプレーン（プレーンＡ２、Ａ３）にインデックスパターンが存在する場合に、それら２個のインデックスパターンとの関係で加えられる斥力ポテンシャルを説明する図である。異なるプレーン間のインデックスパターンとの関係であるから、β＝１が適用されインデックスパターンＤｏのポテンシャルとしてインデックスパターン２個分の２×βＥ（ｒ）のポテンシャルが加えられる。
【００４５】
図５（ｃ）は、上記の２つの場合に加え、異なるプレーンの同一位置にインデックスパターンが存在してインデックスパターンの重なりが存在する場合に、それらのインデックスパターンとの関係で加えられる斥力ポテンシャルを説明する図である。この場合は、図５（ａ）および（ｂ）の場合に加え、注目インデックスパターンＤｏのプレーンＡ１と異なるプレーンＡ３の同じ位置にインデックスパターンが存在する。これによって、同プレーンの斥力ポテンシャル１×αＥ（ｒ）と、同じ位置の異なるプレーンの１個のインデックスパターンによる斥力ポテンシャル１×βＥ（０）と、異なるプレーンの２個のインデックスパターンによる斥力ポテンシャル２×βＥ（ｒ）と、重なる数ｎ＝２でγ＝３が適用される、重なりによる斥力ポテンシャルγｓ（２）×Ｅ（ｒ）のポテンシャルが加えられる。この結果、図５（ｃ）に示すインデックスパターンの配置において注目インデックスパターンＤｏが存在することによる斥力ポテンシャルの合計は、１×βＥ（０）＋１×αＥ（ｒ）＋２×βＥ（ｒ）＋γｓ（２）×Ｅ（ｒ）となる。
【００４６】
図５（ｄ）は、図５（ｃ）に示すインデックスパターンの配置において、インデックスパターンＤｏを移動させることにより、そのインデックスパターンの斥力ポテンシャルの合計が変化することを説明する図である。図５（ｄ）に示すように、インデックスパターンＤｏ（プレーンＡ１のインデックスパターン）が同じプレーンの隣の位置に移ると、そのインデックスパターンＤｏが存在することによる斥力ポテンシャルの合計は次のようになる。すなわち、距離がｒ２、重なり同士の数ｎが０となることなどにより、斥力ポテンシャルの合計はβＥ（１）＋１×αＥ（ｒ２）＋２×βＥ（ｒ２）に変化する。そして、図５（ｃ）に示すインデックスパターンの配置の場合の斥力ポテンシャルの合計１×βＥ（０）＋２×βＥ（ｒ）＋１×αＥ（ｒ）＋γｓ（２）×Ｅ（ｒ）と、図５（ｄ）のインデックスパターンＤｏが移動したことによる斥力の合計とを比較し、この移動前後の斥力ポテンシャルの合計の変化を知ることができる。
【００４７】
なお、この斥力ポテンシャルの合計は、上記の説明では、２つの位置またはインデックスパターンを移動させたときは３つの位置のインデックスパターンによるエネルギーの合計を求めるものとしている。しかし、これは説明を簡易にするためであり、実際は、これらのインデックスパターン以外に存在し得る他の位置のインデックスパターンを含めたインデックスパターンとの関係に基づく斥力ポテンシャルの積分として求められるものであることはもちろんである。
【００４８】
図５（ａ）〜（ｃ）に示したように斥力ポテンシャルの合計が計算される各インデックスパターンの中で、例えばインデックスパターンＤｏが斥力ポテンシャルの合計が最も大きい場合、図５（ｄ）で説明したようにその移動前後の斥力ポテンシャルの変化を求める。そして、移動の前後で最も斥力ポテンシャルの合計が低くなる位置にインデックスパターンＤｏを移動させる。このような処理を繰り返すことによって３つのプレーン全体の総エネルギーを下げることができる。すなわち、３つのプレーンの重なりにおいてインデックスパターンの分布が、低周波数成分が少なく良好に分散された配置となる。
【００４９】
順次配置法
この方法は、上述したように、プレーンにおいてインデックスパターンが未だ配置されていない部分に順次インデックスパターンを配置して行く方法である。例えば、図５（ａ）〜（ｃ）に示した３つのプレーンに対して、順次１つずつインデックスパターンを配置し、それを繰り返すことによってそれぞれのプレーンで、そのプレーンの配置率に応じたインデックスパターンの配置を行う。この場合、先ず、インデックスパターンを配置しようとするときに、その位置のインデックスパターンとプレーンＡ１、Ａ２、Ａ３において既に配置されているインデックスパターンとの間に発生する斥力ポテンシャルを計算する。斥力ポテンシャルの計算自体は、上述の配置移動法で説明したものと同じである。異なるのは、インデックスパターンＤｏを新たに置くと仮定したときに既に配置されている同じプレーンＡ１や異なるプレーンＡ２、Ａ３のインデックスパターンとの間で斥力ポテンシャルを計算する点である。すなわち、図５（ａ）〜（ｃ）に示す例のように、インデックスパターンＤｏが同図の位置に既に置いてあるのではなく、インデックスパターンＤｏを新たに置くと仮定したときに既に配置されている他のインデックスパターンとの間で斥力ポテンシャルを計算する。以上からも明らかなように、未だインデックスパターンが１つも配置されていない最初の段階では、インデックスパターンをどこにおいても斥力ポテンシャルは同じ値となる。
【００５０】
次に、それぞれのプレーン位置に置いたとしたときに計算される斥力ポテンシャルの中で、最小のポテンシャルエネルギーとなる位置を決定する。そして、その最小のエネルギーとなる位置が複数ある場合は、例えば、乱数を用いてその複数の位置の中から１つの位置を決定する。なお、本実施形態では、同じプレーンでは既にインデックスパターンが配置されている位置には重ねて配置しないという条件の下で、最小エネルギーの位置を決定する。これは、重み付け係数や斥力ポテンシャル関数などのパラメータによっては、斥力ポテンシャルの計算において同じプレーンで重ねた場合の方が他のプレーンのインデックスパターンとの関係などでエネルギーが最小となることがある。そして、その場合に、プレーンは１つの位置に１つのインデックスパターンのみが許されるので重なりを禁ずるようにするためである。決定した最小ポテンシャルエネルギーの位置にインデックスパターンを配置する。すなわち、その位置のデータを“１”とする。そして、３つのプレーンＡ１、Ａ２、Ａ３について各１つずつインデックスパターンが配置されたか否かを判定する。配置されていない場合には、上記の処理を繰返す。
【００５１】
プレーンＡ１、Ａ２、Ａ３とこの順で１つずつインデックスパターンを配置すると、３つのプレーンそれぞれの全配置可能位置に対して５０％までインデックスパターンが配置されるまで上記処理を繰り返し、５０％まで配置されると本処理を終了する。
【００５２】
以上説明した順次配置法によっても上述した配置移動法と同様の特性を持つプレーンを得ることができる。すなわち、順次配置法による３プレーンは、それらの重なりにおいてもインデックスパターンが良好に分散されたものとなる。
【００５３】
なお、上述したプレーン製法の他の特徴として、インデックスパターンの配置が規則的に繰り返されるような周期パターンが生成されることはない、ということがある。例えば、千鳥パターンやベイヤー型の配置が繰り返されるような周期性を持ったパターンは生成されない。万が一生成されたとしても、斥力ポテンシャルのパラメータを設定し直すことで周期パターンを避ける状態に収束させることができる。このように本実施形態の展開パターン製法によって生成されるプレーンにおけるインデックスパターンの配置は非周期となる。
【００５４】
以上説明した、基本的な展開パターンの製法を用いた、本発明のいくつかの実施形態にかかる展開パターン製造方法を次に説明する。
【００５５】
なお、上記展開パターンの製法では、プレーンの数を“３”とし、３つのプレーンＡ１〜Ａ３夫々についてのインデックスパターンの配置率を“５０％”としている。一方、以下の本実施形態では、プレーンの数を“４”とし、４つのプレーンＰ１〜Ｐ４夫々についてのインデックスパターンの配置率を“２５％”としている。
【００５６】
（実施形態１）
図７は、本発明の第１の実施形態に係る、ある１つのレベルの４つのインデックスパターンのグループに対応した４つのプレーンＰ１〜Ｐ４それぞれのインデックスパターンの配置に関する斥力を計算するための概念を示す図である。これら４つのプレーンにおけるインデックスパターンの配置を定める処理では、それぞれのプレーンにおけるインデックスパターン同士の斥力ポテンシャルや、１つの、例えばプレーンＰ１におけるインデックスパターンと他のプレーンＰ２〜Ｐ４におけるインデックスパターンとの間の斥力ポテンシャルを計算する。そして、前述した斥力ポテンシャル計算の結果に基づき、プレーンＰ１〜Ｐ４のインデックスパターンの配置を定める。
【００５７】
図８は、本実施形態に係わる４つのプレーンＰ１〜Ｐ４におけるインデックスパターン配置を定める処理を示すフローチャートである。同図に示す処理は、配置移動法によるインデックスパターンの配置方法を示している。
【００５８】
先ず、ステップ４０１で、２５％濃度のグレーの画像を得る。次に、ステップ４０２で、誤差拡散法などによって２値化を行う。これにより、それぞれのプレーンについて、全体の配置可能位置の２５％にインデックスパターンが配置された初期配置を得ることができる。なお、この配置はプレーンごとに同じ配置であっても異なる配置であってもよい。このように２５％ずつの配置とするのは、各レベルについてインデックスパターンの４つのグループがある場合、最終的に製造される展開パターン全体において、それぞれのグループ（図３に示す例ではグループ１〜４）のインデックスパターンを同じ割合で配置するためである。ここで、各プレーンのサイズは、本実施形態では１２８画素×１２８画素であり、その１画素は、縦横とも２値画像データの１／２の解像度である。
【００５９】
次に、ステップ４０３で、上記のようにして得たグループ１〜４のプレーンの総てのインデックスパターンについて、前述したようにして斥力ポテンシャルを計算する。これにより、総てのインデックスパターンの斥力ポテンシャルを合計した総エネルギーが求まる。さらに、ステップ４０４では、総てのインデックスパターンについて順に所定の範囲内で斥力ポテンシャルが小さくなる位置にインデックスパターンを移す。このような処理を繰り返していくことによって、総エネルギーを低下させていく。すなわち、この総エネルギーが徐々に順次減少していく過程は、インデックスパターンの配置が順次分散性を高める過程、つまりインデックスパターンの配置の低周波数成分が少なくなっていく過程である。
【００６０】
ステップ４０５で、Ｓ４０４における総エネルギーの低下率を計算しそれが所定以下であると判断するとステップ４０５でＳ４０４の処理を終了させ、各プレーンの最終的な配置、すなわち各グループの展開パターンのテーブルを決定する。
【００６１】
以上の処理を実行することによって、レベルごとに４つのグループ１〜４（プレーンＰ１〜Ｐ４）の展開パターンを作成することができる。そして、レベルごとの展開パターンは、図１４のようにテーブル化されて、展開パターンテーブル格納ユニット２０２４（図２）に格納される。なお、図１４に示す例では、レベル０および４については、グループ間でインデックスパターンが同じである。このような場合、これらのレベルの展開パターンは、特に上述した本実施形態の処理を用いる必要はない。従って、上述の処理によらず予め定めるようにしてもよい。
【００６２】
図１４は、図８に示す処理によって作成された４つのグループ１〜４の展開パターンをテーブル化したものであって、ここでは、展開パターンに含まれる４種類のインデックスパターンＡ〜Ｄの配置を模式的に示している。この図１４に示されるように、本実施形態では、展開パターンにおける複数種類（ここでは４種類）のインデックスパターンそれぞれの配置が、非周期で、かつ低周波数成分が高周波数成分より少ない特性を有するようになっている。なお、「低周波数成分」並びに「高周波数成分」の定義については後述する通りである。
【００６３】
図１４に示されるように展開パターンは、多値レベル（レベル０〜４）夫々について画素位置（Ｘ，Ｙ）に対応付けられたインデックスパターンＡ〜Ｄを有している。そして、本実施形態では、多値レベルと画素位置に応じてインデックスパターンＡ〜Ｄを選択する。例えば、画素位置（１、１）の多値レベルが“１”であれば、画素位置（１、１）についてレベル１のインデックスパターンＤを選択する。同様に、画素位置（２、０）の多値レベルが“２”であれば、画素位置（２、０）についてレベル２のインデックスパターンＡが選択する。そして、こうして選択されたインデックスパターンにより、各画素のドットデータ（２値画像データ）が決定される。
【００６４】
図９は、図１４に示される展開パターンの中から一つのプレーン（グループ１のインデックスパターンＡ）の展開パターンを抜き出し、当該抜き出した箇所をドットで示した図である。同図に示すように、１つのプレーンのインデックスパターンの配置は非周期であり、かつ分散性が良好であることがわかる。これは他のプレーン２〜４についても同様である。そして、インデックスパターンの配置の分散が良い本実施形態では、後述の図１１に示すように、低周波数の成分が高周波数の成分より少ないものとなっている。また、以上のように、それぞれのプレーン１〜４のインデックスパターンの配置が非周期でかつ分散が良好であるということは、インデックスパターンの配置が相補的であることから、複数（４つ）のプレーン相互のインデックスパターンの配置も非周期で、かつ良好に分散していることを意味する。これにより、同じインデックスパターンは偏りなく配置される。
【００６５】
次に、本実施形態の展開パターンにおけるインデックスパターンの配置を、その周波数特性を示すパワースペクトルによって評価する。以下で説明するパワースペクトルは、１２８画素×１２８画素のサイズを有する展開パターンのうち、ある一つのグループのインデックスパターン（ここでは、インデックスパターンＡ）の配置についてパワースペクトルを求めたものである。ここで、パワースペクトルは、２次元空間周波数を１次元として扱える、「T. Mitsa and K. J. Parker, “Digital Halftoning using a Blue Noise Mask”, Proc. SPIE 1452, pp.47-56(1991)」に記載のradially averaged power spectrum に基づくものである。
【００６６】
なお、本明細書において「低周波数領域」とは、周波数成分が存在する空間周波数領域のうち、半分より低い側の周波数領域を指す。そして、この低周波数領域に存在する周波数成分を「低周波数成分」という。一方、「高周波数領域」とは、周波数成分が存在する空間周波数領域のうち、半分より高い側の周波数領域を指す。そして、この高周波数領域に存在する周波数成分を「高周波数成分」という。
【００６７】
また、「低周波数成分が高周波数成分よりも少ない」とは、低周波数領域に存在する周波数成分（低周波数成分）の積分値が高周波数領域に存在する周波数成分（高周波数成分）の積分値よりも小さい状態を指す。
【００６８】
図１１は、図９（あるいは図１４）に示すインデックスパターン（インデックスパターンＡ）の配置ついて周波数解析を行った結果を示す図である。「自相互分散展開」と記されているのが、本実施形態の展開パターンが示す周波数特性である。同図には、比較のため、従来の周期的にインデックスパターンを展開する「タイル展開」と擬似乱数の配列に従ってインデックスパターンを展開する「擬似ランダム展開」による周波数解析の結果も併せて示している。
【００６９】
図１１に示すように、従来のタイル展開は、インデックスパターンが複数組み合わされたより大きなサイズのパターンが周期的に繰り返されるため、その周波数のところにピークが存在する（図中の矢印）。そのため、記録ヘッドの吐出精度誤差や記録装置の機械的誤差に弱い。すなわち、このような誤差よりインクの着弾位置ズレが生じると、それがテクスチャーとして視認されやすい。一方、擬似ランダム展開は、インデックスパターンの配置が非周期的であるため、インクの着弾位置ズレが生じても、タイル展開のようなテクチャーは視認されにくい。しかし、低周波数領域に比較的大きな成分を持っているため、ノイズが目立ってザラツキ（粒状感）があるように感じられてしまう。
【００７０】
従来の二つの展開方法と比較して、本実施形態の展開パターンは、タイル展開のような周期性を示すようなピークがない。つまり、非周期の特性を示すため、インクの着弾位置ズレが生じてもテクスチャーが発生しにくい。これに加え、本実施形態の展開パターンは、低周波数領域の成分が高周波数領域の成分より少ない特性を有している。このように低周波数成分が少ないことから、インデックスパターンの配置の偏りが少なく分散性が良いといえる。そのため、擬似ランダム展開のような粒状感が少ない。
【００７１】
なお、以上説明した図１１の周波数特性は、展開パターンに含まれるインデックスパターンＡの配置に関するものであったが、その他のインデックスパターンＢ、Ｃ，Ｄにおいても同様の結果が得られる。すなわち、インデックスパターンＢ、Ｃ、Ｄの夫々の配置に関しても、非周期で、且つ低周波数成分が高周波数成分よりも少なくなっている。
【００７２】
図１０は、レベル１のベタ画像（総ての画素がレベル１のデータである画像）を、図８に示す処理によって作成された展開パターン（図１４の展開パターン）に従って２値化して得られるドットデータ（２値画像データ）を示す図である。具体的には、図１４の展開パターンに従って、レベル１を示す画素の位置に対応したインデックスパターンＡ〜Ｄを選択し、選択したインデックスパターンをドットデータとしたものである。
【００７３】
図１０から明らかなように、本実施形態の展開パターンに従ったインデックスパターンによって２値化されたドットデータ（２値画像データ）は、ドットあるいはそれらが集まった大きな塊が分散していることがわかる。つまり、インデックスパターンの配置と同様、ドットデータの配置も、非周期且つ高分散となる。因みに、ドットデータの周波数特性も、上記図１１に示した周波数特性と同様、低周波数成分が高周波数成分よりも少なくなっている。
【００７４】
図１２は、図１０に示すドットデータを偶数カラムと奇数カラムに分割した場合のドットの配置を示す図である。図１２に示すように、ある方向に分割した場合でもその方向によらずドットの偏りがない配置を得ることができる。
【００７５】
本実施形態によれば、各グループ（各プレーン）におけるインデックスパターンの配置について、プレーン内およびプレーン間に斥力ポテンシャルの計算が適用される。これにより、非周期的かつ低周波数成分が少ない展開パターンを得ることができ、その結果、分散性が高く、インクの着弾位置ズレによるテクスチャーの発生がない２値化画像を得ることができる。
【００７６】
（実施形態２）
本発明の第二の実施形態は、メモリーにおいて展開パターンのテーブルが占める容量を小さくするために小さいサイズの展開パターンを作成し、それを繰り返し用いる形態に関する。例えば、比較的小さいサイズの横６４画素×縦３２画素の大きさを有する展開パターンを繰り返し使用して２値化処理を行う。すなわち、展開パターンにおける画素配列の１つの方向にその展開パターンを所定量シフトするとともに、他の方向では展開パターンの一端部が他の端部に隣接するように移動させていることにより展開パターンを繰り返して用いる。本実施形態のインデックスパターンは、上述した第一の実施形態と同じであり、レベル０〜４のそれぞれについて、図３に示す４つのグループ１〜４のパターンである。また、それらの展開パターンも第一の実施形態と略同様に作成される。
【００７７】
異なる点は、それぞれのプレーンの初期画像を２５％濃度の、横６４画素×縦３２画素のサイズで作成する。また、展開パターンを繰り返し使用するために、斥力ポテンシャル計算を行うときに左端は右端と隣接し、上端は下端と隣接しているという境界条件を定める。例えば、斥力ポテンシャルの計算において、右端近傍の位置にインデックスパターンを移動させあるいは配置しようとするときに、左端近傍の位置に既に配置されたインデックスパターンとの距離は、約１プレーン分の幅の距離でなく、左端と右端が隣接したときの近い距離とする境界条件を定める。
【００７８】
以上の形態において、図１３に示すように展開パターンを記録ヘッドのノズル配列方向にシフトさせて用いる場合は、シフト量に応じた境界条件を入れる。ここで、図１３ではシフトする方向をノズル列方向としているが、それに限らずノズル列と垂直な方向あるいはその両方であっても良い。
【００７９】
以上のように、展開パターンのサイズを小さくしても、斥力ポテンシャルを計算するときにそれに応じた境界条件を入れることにより第一の実施形態と同様、非周期的かつ低周波数成分が少ない特性をもつ展開パターンを得ることができる。その結果、分散性が高く、インクの着弾位置ズレによるテクスチャーの発生がない２値画像データを得ることができる。
【００８０】
（その他の実施形態）
上述の各実施形態では、インクジェット記録装置が本発明のデータ処理装置として機能し、ドット配置パターン（インデックスパターン）を用いた２値化処理およびそれに関連した処理を行うものとしが、本発明の適用が、この構成に限られないことはもちろんである。例えば、パーソナルコンピュータなどが本発明のデータ処理装置として機能し、上記各実施形態で説明したドット配置パターン（インデックスパターン）を用いた２値化処理を行う構成であってもよい。
【００８１】
さらに、本発明は、上述した各実施形態の機能を実現する、図８に示したフローチャートの手順を実現するプログラムコード、またはそれを記憶した記憶媒体によっても実現することができる。また、システムあるいは装置のコンピュータ（またはＣＰＵやＭＰＵ）が記憶媒体に格納されたプログラムコードを読出し実行することによっても達成される。この場合、記憶媒体から読み出されたプログラムコード自体が前述した実施形態の機能を実現することになり、そのプログラムコードを記憶した記憶媒体は本発明を構成することになる。
【００８２】
プログラムコードを供給するための記憶媒体としては、例えば、フロッピー（登録商標）ディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＣＤ−Ｒ、磁気テープ、不揮発性のメモリカード、ＲＯＭなどを用いることができる。
【００８３】
また、コンピュータが読み出したプログラムコードを実行することにより、前述した実施形態の機能が実現されるだけでなく、そのプログラムコードの指示に基づき、コンピュータ上で稼動しているＯＳが実際の処理の一部または全部を行うものであってもよい。
【図面の簡単な説明】
【００８４】
【図１】（ａ）および（ｂ）は、従来のインデックスパターンの展開の仕方を説明する図である。
【図２】本発明の一実施形態に係る、主にインクジェット記録装置におけるインデックスパターンを用いた２値化処理の構成を示すブロック図である。
【図３】本発明の一実施形態に係る、５値の画像データのレベル０〜レベル４ごとのインデックスパターンの４つのグループを示す図である。
【図４】（ａ）および（ｂ）は、本発明の一実施形態に係るインデックスパターンの他の例を示す図である。
【図５】（ａ）〜（ｄ）は、本発明の一実施形態に係る斥力ポテンシャルの計算と総エネルギーの減衰処理を模式的に示す図である。
【図６】本発明の一実施形態に係る基本斥力ポテンシャルＥ（ｒ）の関数を模式的に示す図である。
【図７】本発明の第一の実施形態に係る、ある１つのレベルの４つのインデックスパターンのグループに対応した４つのプレーンそれぞれにおけるインデックスパターンの配置に関する斥力を計算するための概念を示す図である。
【図８】第一の実施形態に係わる４つのプレーンにおけるインデックスパターン配置を定める処理を示すフローチャートである。
【図９】図８に示す処理によって得られた一つのプレーンの展開パターンを示す図である。
【図１０】レベル１のベタ画像を、図８に示す処理によって作成された展開テーブルに従って２値化して得られるドットデータ（２値画像データ）を示す図である。
【図１１】図９に示すインデックスパターンの配置ついて周波数解析を行った結果を示す図である。
【図１２】図１０に示すドットデータを偶数カラムと奇数カラムに分割した場合のドットの配置を示す図である。
【図１３】本発明の第二の実施形態に係わり、展開パターンを記録ヘッドのノズル配列方向にシフトさせて用いるため例を説明する図である。
【図１４】図８に示す処理によって得られた４つのグループの展開パターンをテーブル化したものであって、展開パターンに含まれる４種類のインデックスパターンＡ〜Ｄの配置を模式的に示した図である。
【符号の説明】
【００８５】
２０１    ホスト装置
２０２    インクジェット記録装置
２０２１    ＣＰＵ
２０２２    受信バッファ
２０２３    インデックスパターン割付モジュール
２０２４    展開パターンテーブル格納ユニット
２０２５    展開バッファ
２０２６    ドット配置パターン格納ユニット
【図１】
0004519875.tif 000003
【図２】
0004519875.tif 000004
【図３】
0004519875.tif 000005
【図４】
0004519875.tif 000006
【図５】
0004519875.tif 000007
【図６】
0004519875.tif 000008
【図７】
0004519875.tif 000009
【図８】
0004519875.tif 000010
【図９】
0004519875.tif 000011
【図１０】
0004519875.tif 000012
【図１１】
0004519875.tif 000013
【図１２】
0004519875.tif 000014
【図１３】
0004519875.tif 000015
【図１４】
0004519875.tif 000016