(57)【特許請求の範囲】
【請求項１】
３次元カラー変換テーブルを含むディジタル画像形成システムにおいて、輝度、レッドとグリーンとの第１の色差、及びイエローとブルーとの第２の色差の３つの成分からなる入力カラー値を出力カラー値へ変換するカラー変換方法であって、前記３次元カラー変換テーブルは、各次元が入力カラー値の成分に対応し、各次元に複数のノードが設けられ、各ノードが出力カラー値を定義付けるノードの配列を含み、
前記カラー変換方法が、
前記入力カラー値を受け取るステップと、
前記入力カラー値の２つの成分に対してマルチレベルハーフトーン処理を施して、該２つの成分に対応する２次元の各々について１つのノードを選択するステップと、
前記入力カラー成分のうち残った他の１つの入力カラー成分を用いて選択された前記複数のノード同士を補間して出力カラー値を求めるステップと、
を備える前記カラー変換方法。
【請求項２】
入力カラーを定義付けると共に輝度、レッドとグリーンとの第１の色差、及びイエローとブルーとの第２の色差の３つの入力カラー成分を受け取る受取器と、前記入力カラーを出力カラーへ変換する変換プロセッサと、前記出力カラーを出力用媒体上へ生成する画像出力装置と、を有する電子画像形成システムにおいて、前記変換プロセッサが、
各次元が入力カラー値の成分に対応し、各次元に複数のノードが設けられ、各ノードが出力カラー値を定義付けるノードの配列を含む３次元テーブルと、
前記入力カラー成分の２つの成分を受け取ると共に、該２つの成分に対応する２次元について１つのノードを識別する第一の量子化器と、
前記入力カラー成分のうちの残った１つの入力カラー成分に基づいて識別されたノードに含まれる予め定義付けられたカラーデータから出力カラー値を計算する補間器と、
を備える前記変換プロセッサ。
【発明の詳細な説明】
【０００１】
本発明は、多次元カラー変換テーブルを含むディジタル画像形成システムにおいて入力カラー値を出力カラー値へ変換するカラー変換方法及び変換プロセッサに関する。
【発明の属する技術分野】
【０００２】
コンピュータ及び他の電子装置は一般的にＲＧＢ等の３−Ｄ（三次元）座標で色を表示する。一方、多くのプリンタは、入力値、例えば、Ｒ、Ｇ、Ｂに対応するシアン、マゼンタ、イエロー（ＣＭＹ）を三次元着色剤空間で又はシアン、マゼンタ、イエロー、及びブラック（ＣＭＹＫ）を四次元着色剤空間で印刷する。デバイス非依存色空間（カラースペース）はしばしば画像の中間的表示として使用される。このようなデバイス非依存色空間として一般的に選択されるのはＬＣ1Ｃ2で総称的に表されるルミナンス（輝度）−クロミナンス（色差）空間である。Ｌ成分はルミナンス（輝度）又は明度を表し、Ｃ1 及びＣ2はそれぞれレッド−グリーン及びイエロー−ブルーの変化を表す二つのクロミナンス座標である。このような空間の例としてはＣＩＥＬＡＢ（国際照明委員会（ＣＩＥ）が１９７６年に定めた色空間の一つ）が挙げられる。従って、変換は入力ＲＧＢからＬＣ1Ｃ2へ及びＬＣ1Ｃ2からプリンタ着色剤空間へ導き出される。
【０００３】
【従来の技術】
本発明は広範囲の変換に適用可能であるが、論議はルミナンス（輝度）−クロミナンス（色差）空間ＬＣ1Ｃ2からプリンタ着色剤空間ＣＹＭＫへの変換に集中する。このような変換は一般的に３−Ｄルックアップテーブル（ＬＵＴ）によって実行され、この３−ＤＬＵＴは各ディジタルＬＣ1Ｃ2入力値を、これらの入力値がプリンタによって受信される前に、対応する出力ＣＭＹＫ値へ変換する。
【０００４】
望ましい着色剤作用を有するプリンタは、シアン対レッド、マゼンタ対グリーン、及びイエロー対ブルーからなる１対1の対応を有する。これは、印刷時に、シアン着色剤はレッド光、マゼンタ着色剤はグリーン光、イエロー着色剤はブルー光のみを吸収することを意味する。しかしながら、一般的なプリンタ着色剤はこの理想的な作用を逸脱し、実際は、意図される吸収バンド以外の電磁スペクトルのバンド内で光を吸収してしまう。これらの所謂望ましくない吸収によって着色剤同士の相互作用が発生し、これにより、プリンタを駆動するディジタル値と得られる測色応答との間に複雑な非線型の関係が生じる。「測色」とラベル付けされた応答又は他の値は平均的な人間の観察者が見た時の印刷されたカラーの測定を称し、ＣＩＥＬＡＢ等のデバイス非依存色座標で表される。ＣＭＹＫ値の全体範囲に及ぶ測色応答のモデリングは従って通常単純な関数によっては達成することはできず、実際、多くのパラメータ及び測定回数を要する。プリンタを十分に特徴付けるのに要する測定回数は優に１０００回を超える。このような複雑な関数を割安な演算コストで表すため、通常、測色空間とＣＭＹＫ値とのマッピングを近似するカラー補正ＬＵＴが作成される。特に、このカラー補正ＬＵＴは、非線型性及び着色剤の望ましくない吸収を補正する。ＬＵＴは、ルミナンス−クロミナンス空間ＬＣ1Ｃ2において特定されたあらゆる入力カラー毎に、印刷及び測定時に要求されるＬＣ1Ｃ2カラーを出力する対応するＣＭＹＫを検索するが、この出力カラーがデバイスの再生可能な全領域内にあることが前提とされる。
【０００５】
ＬＵＴを作成するためにＣＭＹＫディジタル値の予め定義付けられたセットがプリンタへ送られる。プリンタは対応するセットのカラーパッチを印刷する。キャリブレーション（較正）用カラーパッチが測定されると共に、測色用ＬＣ1Ｃ2座標が各パッチ毎に即ち予め定義付けられた各ＣＭＹＫ毎に求められる。測定されたＬＣ1Ｃ2座標の各々は次に三次元空間内での三次元ベクトル位置を識別する。各ＬＣ1Ｃ2座標は一般的にＬ、Ｃ1、及びＣ2の各成分毎に８−ビット値で表される。このような２４ビットのＬＣ1Ｃ2座標は２５６3個のロケーション（記憶場所）にアドレスすることが可能であるが、ルックアップテーブル（ＬＵＴ）は一般に各ノードがＣＭＹＫ値を記憶する１６ｘ１６ｘ１６（４０９６）テーブルロケーションのようなより小さなサイズへ小分けされる。ＬＣ1Ｃ2の中間地点でのＣＹＭＫ値は、ＬＵＴノード同士の補間の形態によって決まる。ルックアップテーブルのサイズは、ルックアップテーブル（即ち、出力上の忠実度）の所望される精度と多数の値を記憶するための費用との折合による。従って、キャリブレーション用パッチが生成された後は、各測定されたＬＣ1Ｃ2座標は対応する既知のＣＭＹＫ値を有することになる。あいにく、これらのＬＣ1Ｃ2座標は概してルックアプテーブルのノードロケーション（即ち、三次元交差地点）と完全には一致しないものである。従って、ノード上に配置されるＣＭＹＫ値は、シェパード(Shepard)のアルゴリズム等の多次元データフィッティング（当てはめ）技術などによって推定される。
【０００６】
【発明が解決しようとする課題】
ルックアップテーブルの動作は図示することによって理解される。図１によれば、入力ＬＣ1Ｃ2値１０はプリンタ特有のＣＭＹＫ値への変換の為にテーブル１２へ送られる。この変換は入力ＬＣ1Ｃ2座標ロケーション１０に最も近接したノード１４に対応する既知のＣＭＹＫ値を補間することによって達成される。カラーは三次元で定義付けられるので、補間も同様に三次元で行われる。３−Ｄ補間技術の一般的な例としては、三線、四面体、及びプリズム補間が挙げられる。これらの補間のうち、四面体補間は処理が最も速い方法であり、四つのノード同士の補間で済む。（三線及びプリズム技術はそれぞれ八つ及び六つのノードを使用する）。これらの技術はすべて各画素の各出力信号毎にいくつかの乗算、加算、桁送り、及び比較の演算が行われることを要求すると共にしばしば専用のハードウェアと共に実行される。次に、この補間されたＣＭＹＫ値はプリンタによって出力される。残念ながら、三次元補間は多くのアプリケーションに相当な計算のための負担をかけるという欠点がある。
【０００７】
【課題を解決するための手段】
本発明は上述した問題等を克服する新しく改良されたカラー補正方法を提供する。
【０００８】
本発明によれば、マルチレベルハーフトーン処理と補間を組み合わせるルックアップテーブルベースのカラー変換のための改良された方法が提供される。特に、ルックアップテーブルへ入力される色空間がルミナンス−クロミナンス座標内にある場合は、クロミナンス座標の各々へマルチレベルハーフトーン処理が施され、これら二つの座標に沿ってルックアップテーブルのノードとの符合が達成される。このマルチレベルハーフトーン処理のステップはクロミナンスチャネルの高い空間周波数において誤差を引き起こす。しかしながら、これらの誤差は人間の視覚系によって容易には知覚されない。最終的には、出力カラー値を得るために１−Ｄ補間がルミナンス成分上で実行される。
【０００９】
本発明の第１の態様は、３次元カラー変換テーブルを含むディジタル画像形成システムにおいて、輝度、レッドとグリーンとの第１の色差、及びイエローとブルーとの第２の色差の３つの成分からなる入力カラー値を出力カラー値へ変換するカラー変換方法であって、前記カラー変換テーブルは、各次元が入力カラー値の成分に対応し、各次元に複数のノードが設けられ、各ノードが出力カラー値を定義付けるノードの配列を含み、前記カラー変換方法が、前記入力カラー値を受け取るステップと、前記入力カラー値の２つの成分に対してマルチレベルハーフトーン処理を施して、該２つの成分に対応する２次元の各々について１つのノードを選択するステップと、前記入力カラー成分のうち残った他の１つの入力カラー成分を用いて選択された前記複数のノード同士を補間して出力カラー値を求めるステップと、を備えるカラー変換方法である。
【００１０】
本発明の第１の態様によれば、入力カラー値を出力カラー値へ変換するカラー変換方法は、ノードの配列を有する３次元カラー変換テーブルを含む。３次元カラー変換テーブルは、各次元が入力カラー値の成分に対応し、各次元に複数のノードが設けられ、各ノードは出力カラー値を定義付ける。このカラー変換方法は輝度、レッドとグリーンとの第１の色差、及びイエローとブルーとの第２の色差の３つの成分からなる入力カラー成分を備えたこの入力カラー値を受け取ることと、この入力カラー値の２つの成分に対してマルチレベルハーフトーン処理を施して、該２つの成分に対応する２次元の各々について１つのノードを選択することを含む。入力カラー成分のうち残った他の１つの入力カラー成分を用いて、選択された複数のノードを補間して出力カラー値を求める。
【００１１】
本発明の第２の態様は、入力カラーを定義付けると共に輝度、レッドとグリーンとの第１の色差、及びイエローとブルーとの第２の色差の３つの入力カラー成分を受け取る受取器と、前記入力カラーを出力カラーへ変換する変換プロセッサと、前記出力カラーを出力用媒体上へ生成する画像出力装置と、を有する電子画像形成システムにおいて、前記変換プロセッサが、各次元が入力カラー値の成分に対応し、各次元に複数のノードが設けられ、各ノードが出力カラー値を定義付けるノードの配列を含む３次元テーブルと、前記入力カラー成分の２つの成分を受け取ると共に、該２つの成分に対応する２次元について１つのノードを識別する第一の量子化器と、前記入力カラー成分のうちの残った１つの入力カラー成分に基づいて識別されたノードに含まれる予め定義付けられたカラーデータから出力カラー値を計算する補間器と、を備える変換プロセッサである。
【００１２】
本発明の第２の態様によれば、受取器と変換プロセッサを含む電子画像形成システムが提供されている。この変換プロセッサは各次元が入力カラー値の成分に対応し、各次元に複数のノードが設けられ、各ノードが出力カラー値を定義付けるノードの配列を含む３次元テーブルを含む。第一の量子化器は、入力カラー成分の２つの成分を受け取ると共に、該２つの成分に対応する２次元について１つのノードを識別する。補間器は、入力カラー成分のうちの残った１つの入力カラー成分に基づいて識別されたノードに含まれる予め定義付けられたカラーデータから出力カラー値を計算する。
【００１３】
本発明の利点は高価な３−Ｄ補間がはるかに単純な処理であるハーフトーン処理及び１−Ｄ補間へ還元され、これによって視覚的品質を殆ど損わずに計算のための実質的なコストの節約が可能になることである。本発明の他の利点は、ハーフトーン処理及び１−Ｄ補間のステップによって、ルックアップテーブル内で予め演算される量のための記憶装置及びメモリの量が少なくて済み、これによって他のグリッド構造をルックアップテーブルに使用できることである。
【００１４】
本発明は使用されている部品及びこれら部品の配列をあらゆる工程及び構成において変形が可能である。従って、添付されている図面は図示だけを目的としており、本発明の範囲を限定するものではない。
【００１５】
【発明の実施の形態】
３−Ｄ（３次元）ルックアップテーブル（LUT）ベースのカラー補正に要する補間演算のコストを削減するための方法が提供されている。入力カラー値はルミナンス（輝度）−クロミナンス（色差）（ＬＣ1Ｃ2）座標で表されるが、本発明はその機能性を失わずに他の入力フォーマットにも同様に適用することが可能である。マルチレベルハーフトーン処理は、Ｃ1及びＣ2を限定されたセット(集合)のレベル即ちＣ1及びＣ2に沿ったＬＵＴノードロケーションに符合するレベルに制限するように、使用される。ＬＵＴ演算によってＬＣ1Ｃ2座標での３−Ｄ補間からＬ座標だけに沿って行われる１−Ｄ補間へ削減される。従来の技術の３−Ｄ補間に比べ、本発明が計算のためのコスト面において十分な利得を得つつ受容可能な画像が得られることが試験によって示されている。
【００１６】
本発明は、カラー補正を３−Ｄ補間の問題から１−Ｄ補間の問題へ還元させることによってカラー補正のコストと処理速度が有効に削減されることを認識する。一般に、デバイス非依存の画像値は、例えば、ＣＩＥＬＡＢカラースペース等のルミナンス−クロミナンス座標として、記憶される。概して、人間の視覚系は高い空間周波数ではルミナンス上の歪みよりもクロミナンス上の誤差に対して感応しない。即ち、観察者は高い空間周波数においてはクロミナンスよりルミナンスにおける小さな誤差を検知しやすい。この認識が本発明の基礎となっている。
【００１７】
具体的には、図２にはＣ1−Ｃ2単一平面２０が示されている。入力カラー２２は入力カラー成分Ｌin、Ｃ1in、及びＣ2inを備える。後続のカラー補正演算において計算のためのコストを削減するためにＣ1及びＣ2チャネルでいくらかの歪みがあってもよい。従って、Ｃ1成分は、各々が（Ｃ2及びL軸に沿った）複数のノードに符合する有限セット（集合）の次元又はレベル２４、２５、２６、２７、及び２８に限定される。同様に、Ｃ2成分はレベル３１、３２、３３、３４、及び３５に合致するように限定される。入力カラー値がＬＵＴ内のレベル又は予め定義付けられた次元に常に強制的に符合させられる場合、その次元に沿って補間を実行する必要がないことは当業者に理解されよう。つまり、三つの入力成分のうち、二つ（Ｃ1及びＣ2）がマルチレベルハーフトーン処理によって求められる。これにより、カラー補正は以下にも示すように残った軸（例えば、ルミナンス又はＬ）に沿った１−Ｄ補間へと還元される。
【００１８】
信号の精度を小さなセットのレベル又は次元へ還元することは量子化の問題である。一般的なカラー補正ＬＵＴは、Ｌ、Ｃ1、及びＣ2の各次元に沿って１６個のノードロケーションを割当てることによって（図１参照）、結果的に、合計１６3＝４０９６ノードのエントリとなる。従って、本発明は、一般的な８−ビット入力に対しては、各クロミナンスチャネルごとに２５６の可能な入力レベルを１６の出力ノードロケーション（記憶場所）へ量子化マッピングすることを要求する。直接的量子化には各入力値を最も近接したノードへマッピングすることを伴う。しかしながら、実験によれば、この方法によって輪郭上に望ましくないアーティファクト(人工的発生物)が生じてしまうことがわかっている。
【００１９】
量子化された画像に生じるこのようなアーティファクトを最小にする方法は、ディザリング又はより一般的にはマルチレベルハーフトーン処理を介して行われる。これによって、量子化による誤差を敢えて高い空間周波数へ効果的に変換し、人間の視覚系における他の無感応性を利用する。そして、好ましくは、ハーフトーン処理を高い空間周波数におけるクロミナンスチャネルC1及びC2.に施すことによって大幅に減少した人間の視覚系における感応性を高い周波数のクロミナンスの誤差として利用する。マルチレベルハーフトーン処理には二つの公知の方法即ちスクリーニング（遮蔽）と誤差拡散がある。本発明では、あるゆるスクリーニング技術の計算上の容易さを享受できることからブルーノイズの確率論的スクリーンを用い、誤差拡散からはより望ましい質的作用を少々取り入れることが好ましい。同等に活用可能なマルチレベルハーフトーン処理における他の方法は拡散ドットスクリーンを用いることである。さらに図２を参照すると、このハーフトーン処理はＬＵＴに入力される入力カラー２２として視覚化される。一般的な例としては、入力カラー２２はC1及びC2.両軸上の予め定義付けられたレベル２４−２８と３１−３５との間に置かれている。先ず、C1軸に対して入力カラー２２をやや下位のレベル２５又はやや上位のレベル２６のどちらの上へ置くべきかが決定される。補間が実行されず、入力カラーはこれらのレベルの一つへ「強制的に」置かれ、これによって誤差が引き起こされる。同様の決定がＣ2軸についても行われ、入力カラー２２を下位のレベル３３又は上位のレベル３４へ強制的に置く。仮想上の図示の為に閾値処理により入力カラーがノード４０に置かれるものとする。
【００２０】
図３に示されるように、変換プロセッサ５０は、三つの成分即ちＬ２２L、Ｃ1２２A、及びＣ2２２Bを備える入力カラー２２を受け取る。図示されている実施の形態において、Ｃ1２２A及びＣ2２２Bは別個に量子化される。マルチレベルハーフトーン処理における最初のステップは、Ｃ1軸に沿ったＬＵＴにおける二つの隣接ノード（例えば、図２のレベル２５及び２６）が形成する間隔内の画素位置（x, y）での座標Ｃ1（ｘ、ｙ）２２Aを正規化することである。正規化された座標Ｃ1’は以下の等式によって求められる。
【００２１】
C1' (x, y)=[C1(x, y)- C1L(x, y)]/[C1G(x, y)- C1L(x, y)]   等式（１）
【００２２】
この数式において、Ｃ1（ｘ、ｙ）は元の入力座標であり、Ｃ1L（ｘ、ｙ）はＣ1（ｘ、ｙ）未満で最も近接したＬＵＴノードであり、Ｃ1G（ｘ、ｙ）はＣ1（ｘ、ｙ）より上で最も近接したＬＵＴノードである。従って、正規化された座標Ｃ1’（ｘ、ｙ）は常に０と１の間にある。この正規化演算は前もって交互に実行され、１−Ｄルックアップテーブルアレイ（配列）内に記憶されてもよい。正規化された値Ｃ1’（ｘ、ｙ）は、次にハーフトーンスクリーン閾値、さもなければ、ディザ信号Ｄ（ｘ、ｙ）５２と呼ばれ、これも０と１の間で正規化される値と比較される。Ｃ1’（ｘ、ｙ）がＤ（ｘ、ｙ）以下である場合、Ｃ1L（ｘ、ｙ）は出力ハーフトーンレベルとして選択され、さもなければ、Ｃ1G（ｘ、ｙ）が選択される。この出力レベルはＣ1qで表される。同一の正規化演算がＣ2（ｘ、ｙ）２２Bに対しても行われる。
【００２３】
この時点で、図４を相互に参照すると、これらの量子化された値又は次元２５及び３３が２−Ｄルックアップテーブル５４へ入力され、量子化されたＣ1のクロミナンスレベル又は平面である例えば２５と、Ｃ2のクロミナンスレベル又は平面３３と、の交差によって定義付けられるルミナンスライン６０が決定される。この結果は１−Ｄルックアップ及び補間器５６を通過し、これによって、入力カラーのＬ成分２２Lは入力カラー値２２を囲僥する二つのノード６４及び６６の間での残りの次元（即ちＬ）を基に（一次元）補間することが可能となり、これによって出力カラーが求められる。
【００２４】
図５には他の実施の形態が示されており、この実施の形態では、３−ＤＬＵＴにおける２−Ｄクロミナンス平面に対して、図２に示されたシ−ケンシャルグリッドの代替として使用されるシ−ケンシャルグリッド７０が用いられる。図２に示されるような厳密に矩形のグリッド配置よりも、このシ−ケンシャルグリッド７０によってＣ1−Ｃ2平面上でのよりフレキシブルなノード配置が可能となる。また、このシ−ケンシャルグリッド７０によってＬＵＴ変換が最大の変化を遂げるか又は視覚的に重要とされる位置にノードが配置されることが可能とされ、その結果、一般的に、ＬＵＴの精度とＬＵＴのサイズの間で優れたトレードオフが行われる。この実施の形態において、マルチレベルハーフトーン処理は入力カラー７２の第１のクロミナンス座標へ適用される。これはＣ1の値又は成分をレベル７４又はレベル７５へマッピングする。所与の例においては、ハーフトーン処理は選択される可能性のあるノードをレベル７４によって定義付けられた二次元空間へ限定する。次のステップにおいては、マルチレベルハーフトーン処理はノード７６とノード７７との間のレベル７４に沿ったＣ2次元で施される。所与の例においては、これにより結果として入力カラーはノード７７にマッピングされる。最終的に、本実施の形態と同様、Ｌ軸に沿った１−Ｄ補間が後に続いて行われる２−Ｄルックアップが出力カラーを得るために施される。
【００２５】
この他の実施の形態における処理による変換プロセッサ８０のダイヤグラムが図６に示されている。この変換プロセッサ８０の処理フローは、一次元に沿ってハーフトーン処理された出力Ｃ1qが第２の軸Ｃ2に沿ってＬＵＴノードのセット（Ｃ1に沿ったノード値が異なれば当然異なるが）を求めるために使用される他は、あらゆる本質的な点で、（図３に示した）第１の実施の形態におけるプロセッサ５０と同一である。従って、ハーフトーン処理は二つのクロミナンス軸に沿って別個に行われるというよりも寧ろ順番に実行される。
【００２６】
表１において、三つの方法、即ち三線補間、四面体補間、及び本発明による本実施の形態が、Ｎ個の出力カラーに対して各画素毎にカラー補正演算を実行する為に必要とされる計算量に関して、それぞれ比較されている。この演算は、１）入力カラーを囲僥するサブセルのノードの検索、及び２）実際の補間を含む。本発明においては、マルチレベルハーフトーン処理における比較ステップがコスト分析に含まれていたことに注目されたい。三つのケース全部において、入力カラーに依存しない任意の量は予め計算され記憶されるものと仮定される。これらの量には隣接ノード間の出力値同士の差が含まれると共に等式（１）によって与えられる入力値の正規化が伴う。表１の最後の行は四面体補間に対して本発明を適用することによって達成されたコストの節約度を示す。
【００２７】
この節約度は三線補間と比較するとかなり大きい。
【表１】
0004429439.tif 000002
【００２８】
【発明の効果】
本発明は、高価な３−Ｄ補間がはるかに単純な処理であるハーフトーン処理及び１−Ｄ補間へ還元され、これによって視覚的品質を殆ど損わずに計算のための実質的なコストの節約が可能となるという利点を有する。本発明はさらに、ハーフトーン処理及び１−Ｄ補間のステップによりルックアップテーブル内で予め演算される量のための記憶装置及びメモリの量が少なくて済み、これによって他のグリッド構造をルックアップテーブルに使用可能とする利点を有する。
【図面の簡単な説明】
【図１】従来の技術の三次元補間ルックアップテーブルを示す線図である。
【図２】  C1C2平面と、入力カラーのオフノード入力を示す図である。
【図３】本発明の実施の形態における論理的流れを示す線図である。
【図４】他の次元でのマルチレベルのハーフトーン処理の後に続いて行われる残りの次元に沿った補間を示す図である。
【図５】本発明によるシーケンシャルグリッドを示す線図である。
【図６】本発明の他の実施の形態における論理的流れを示す線図である。
【符号の説明】
Ｌin、Ｃ1in、及びＣ2in：入力カラー成分
【図１】
0004429439.tif 000003
【図２】
0004429439.tif 000004
【図３】
0004429439.tif 000005
【図４】
0004429439.tif 000006
【図５】
0004429439.tif 000007
【図６】
0004429439.tif 000008