(57)【特許請求の範囲】
【請求項１】
複数の色材を記録媒体に付着させる印刷を印刷装置に実行させるにあたり、前記色材の量の組み合わせである色材量セットに基づいた印刷を前記印刷装置に実行させる印刷制御装置であって、
前記複数の色材は、所定の光源下で略同一の色彩値を有しつつ互いに異なる分光反射率特性を有する等色色材を少なくとも１組含んで構成されており、
前記印刷制御装置は、
印刷対象となる分光反射率又は印刷対象となる色彩値の指示を受付ける受付手段と、
分光反射率を受付けた場合に、前記分光反射率に基づいて決定した前記色材量セットに基づいて前記印刷装置に印刷を実行させる分光反射率印刷制御手段と、
色彩値を受付けた場合に、前記色彩値に基づいて決定した前記色材量セットに基づいて前記印刷装置に印刷を実行させる色彩値印刷制御手段と、
を備えることを特徴とする印刷制御装置。
【請求項２】
前記色彩値印刷制御手段は、略同一の色彩値を有する前記等色色材のそれぞれを一度の印刷にわたって、ほぼ同じ比率で使用した印刷を実行させることを特徴とする請求項１に記載の印刷制御装置。
【請求項３】
前記色彩値印刷制御手段は、略同一の色彩値を有する前記等色色材を、残量に応じて又はほぼ同じ量で、使用する印刷を実行させることを特徴とする請求項２に記載の印刷制御装置。
【請求項４】
前記受付手段は、受信した印刷データを解析することで、該印刷データで指示されている分光反射率又は色彩値を受付けることを特徴とする請求項１に記載の印刷制御装置。
【請求項５】
前記分光反射率印刷制御手段は、波長によって異なる重みを加味しつつ前記ターゲット分光反射率への近似性を評価する評価値に基づいて、指示された分光反射率に近似する分光反射率を前記記録媒体上に示す前記色材量セットを予測し、当該予測したインク量セットに基づいて前記印刷装置に印刷を実行させることを特徴とする請求項１に記載の印刷制御装置。
【請求項６】
複数の色材を記録媒体に付着させる印刷を印刷装置に実行させるにあたり、前記色材の量の組み合わせである色材量セットを指定し、当該色材量セットに基づく印刷を前記印刷装置に実行させる印刷制御装置であって、
前記複数の色材は、所定の光源下で略同一の色彩値を有しつつ互いに異なる分光反射率特性を有する等色色材を少なくとも１組含んで構成されており、
前記印刷制御装置は、
ターゲットの分光反射率をターゲット分光反射率として取得するターゲット分光反射率取得手段と、
波長によって異なる重みを加味しつつ前記ターゲット分光反射率への近似性を評価する評価値に基づいて、当該ターゲット分光反射率に近似する分光反射率を前記記録媒体上に示す前記色材量セットを予測し、当該予測したインク量セットを前記印刷装置に指定する印刷制御手段と、
を備えることを特徴とする印刷制御装置。
【請求項７】
複数の色材を記録媒体に付着させる印刷を実行する印刷装置と、前記色材の量の組み合わせである色材量セットに基づいた印刷を前記印刷装置に指定し、前記色材量セットに基づく印刷を実行させる印刷制御装置とからなる印刷システムであって、
前記複数の色材は、所定の光源下で略同一の色彩値を有しつつ互いに異なる分光反射率特性を有する等色色材を少なくとも１組含んで構成されており、
前記印刷制御装置は、
印刷対象となる分光反射率又は印刷対象となる色彩値の指示を受付ける受付手段と、
分光反射率を受付けた場合に、前記分光反射率に基づいて決定した前記色材量セットに基づいて前記印刷装置に印刷を実行させる分光反射率印刷制御手段と、
色彩値を受付けた場合に、前記色彩値に基づいて決定した前記色材量セットに基づいて前記印刷装置に印刷を実行させる色彩値印刷制御手段と、
を備え、
前記印刷装置は、
前記印刷制御手段が予測した前記色材量セットに基づいて印刷を実行する印刷実効手段を備えることを特徴とする印刷システム。
【請求項８】
複数の色材を記録媒体に付着させる印刷を印刷装置に実行させるにあたり、前記色材の量の組み合わせである色材量セットに基づいた印刷を前記印刷装置に実行させる機能をコンピュータに実現させるための印刷制御プログラムであって、
前記複数の色材は、所定の光源下で略同一の色彩値を有しつつ互いに異なる分光反射率特性を有する等色色材を少なくとも１組含んで構成されており、
前記印刷制御プログラムは、
印刷対象となる分光反射率又は印刷対象となる色彩値の指示を受付ける受付機能と、
分光反射率を受付けた場合に、前記分光反射率に基づいて決定した前記色材量セットに基づいて前記印刷装置に印刷を実行させる分光反射率印刷制御機能と、
色彩値を受付けた場合に、前記色彩値に基づいて決定した前記色材量セットに基づいて前記印刷装置に印刷を実行させる色彩値印刷制御機能と、
をコンピュータに実現させることを特徴とする印刷制御プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、印刷制御装置、印刷システムおよび印刷制御プログラムに関し、特にターゲットを再現させるための印刷制御装置、印刷システムおよび印刷制御プログラムに関する。
【背景技術】
【０００２】
分光的な再現性に注目した印刷方法が提案されている（特許文献１参照。）。この文献においては、ターゲット画像に分光的かつ測色的に一致するような印刷を行うために、プリンティングモデルを使用し、ターゲットの分光反射率（ターゲットスペクトル）にフィッティングするようにプリンタ色（ＣＭＹＫＯＧ）の組み合わせを最適化している。このようにすれば、当該プリンタ色（ＣＭＹＫＯＧ）に基づく印刷を行うことにより、分光的にターゲット画像が再現でき、結果として測色的にも再現性の高い印刷結果を得ることができる。
また特許文献２には、所定の条件下では同じ色彩値を示しつつも、異なる分光反射率を示すインクを、有価証券等の偽造防止等に用いることが記載されている。
【特許文献１】特表２００５−５０８１２５号公報
【特許文献２】特開２０００−３２７９７８号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
前述した特許文献１の技術によれば測色的にも再現性の高い印刷結果を得られるが、インクセット自体が所定の分光反射率分布を有しているため、インクセットをどのように組合せても完全には再現できない波長域が存在する。そのため、既存の色（印刷済みの色、自然界に存在する色、絵画や文化財や文書等に既に着色されている色、絵の具の色等）の分光反射率分布に基づいて印刷を行っても、光源次第では既存の色と異なる色合いに見える可能性がある。但し、特許文献２のように、あえてメタメリズムを発生させて偽造防止に利用したり、元々メタメリズムの存在していた既存の色（印刷済みの色、自然界に存在する色、絵画や文化財や文書等に既に着色されている色、絵の具の色等）を正確に再現したりすることにも意味がある。
【０００４】
本発明は、前記課題にかんがみてなされたもので、メタメリズムの発生をコントロール可能な印刷制御装置、印刷システムおよび印刷制御プログラムを提供する。
【課題を解決するための手段】
【０００５】
前記課題を解決するために、本願の印刷制御装置は、本願の印刷制御装置ににて制御される印刷装置における複数の色材について、所定の光源下で同一の色彩値を有しつつ互いに異なる分光反射率特性を有する等色色材を少なくとも１組含むように構成しつつ、受付手段と分光反射率印刷制御手段と色彩値印刷制御手段とを備えた構成としてある。受付手段は、印刷対象となる分光反射率又は印刷対象となる色彩値の指示を受付けるものである。分光反射率印刷制御手段は、分光反射率を受付けた場合に、前記分光反射率に基づいて決定した前記色材量セットに基づいて前記印刷装置に印刷を実行させるものである。色彩値印刷制御手段は、色彩値を受付けた場合に、前記色彩値に基づいて決定した前記色材量セットに基づいて前記印刷装置に印刷を実行させるものである。このように所定の光源下では同一の色彩値を示しつつも分光反射率特性の異なる等色色材間では、他の所定の光源下では互いに異なる色彩値を示すことになる。従って、これらの等色色材の使用割合をコントロールすることでメタメリズムの発生をコントロール可能になる。すなわち印刷対象を分光反射率で指示された場合には、該分光反射率を再現するように印刷を行ってメタメリズムをコントロールし、所望のメタメリズムを再現したり、メタメリズムの発生しにくい印刷を行ったり、メタメリズムを気にしない印刷を行ったりする。一方、印刷対象を色彩値で指示された場合には、所定の光源下において該色彩値が再現されるように印刷を行う。従って、メタメリズムを気にしない印刷であれば、印刷対象を色彩値で指定することによりメタメリズムを意識しない印刷が可能である。無論、所定の光源下において所望の色彩値が再現される印刷（メタメリズムを気にしない印刷）を、印刷対象を分光反射率で指定する印刷で行っても構わない。
【０００６】
なお、前記受付手段は、実際に前記印刷対象について行われた分光反射率測定に基づいて前記印刷対象の分光反射率や色彩値を指示されてもよいし、ユーザー等から前記印刷対象の分光反射率や色彩値の選択を受付けてもよい。また、前記印刷装置は少なくとも複数の前記色材を前記記録媒体に付着させることができればよく、インクジェットプリンタやレーザープリンタや昇華型プリンタ等の種々の印刷装置に本発明を適用することができる。
【０００７】
また、前記等色色材の色彩値を指定した印刷を実行する際に好適な一例として、前記色彩値印刷制御手段は、略同一の色彩値を有する等色色材のそれぞれを一度の印刷にわたってほぼ同じ比率で使用した印刷を実行させたり、略同一の色彩値を有する等色色材を残量に応じて又はほぼ同じ量で使用する印刷を実行させたりしてもよい。すなわち、メタメリズムが発生しても構わない印刷（分光反射率を一致させる必要の無い印刷）においては、各等色色材を仮想的に１つの色材と見做した印刷を行っても問題ない。したがって、等色色材のすべてを１つの色彩値を有した色材として利用できるようになり、該色彩値による印刷可能枚数が増加することになる。その際、各等色色材の使用頻度に偏りが生じないように印刷を行ったり、色材の残量の多い等色色材から先に使用したりするようにすることで、何れか１つの等色色材が不足する自体の発生を防止できる。特に、略同一の色彩値を有する等色色材のそれぞれを一度の印刷にわたってほぼ同じ比率で使用した印刷を実行させると、印刷物の途中で各等色色材の使用比率が変化しないので、メタメリズムのせいで同じ色彩値であるはずの点同士が光源によっては互いに異なった色に見えてしまったり、濃淡が逆転してしまったりする可能性を防止できる。
【０００８】
また、受付手段が分光反射率と色彩値との何れで印刷指示を受付けているかを判断する際に、前記受付手段は、受信した印刷データを解析することで、該印刷データで指示されている分光反射率又は色彩値を受付けるように構成してもよい。印刷データの解析とは、前記印刷データに対し予め分光反射率と色彩値とを区別可能なフラグを付しておいて該フラグの解析であってもよいし、分光反射率と色彩値とで各データに特有の形式を解析により判別してもよい。
【０００９】
また、色材量セットを予測精度を向上させるために、前記分光反射率印刷制御手段は、波長によって異なる重みを加味しつつ前記ターゲット分光反射率への近似性を評価する評価値に基づいて、指示された分光反射率に近似する分光反射率を前記記録媒体上に示す前記色材量セットを予測し、当該予測したインク量セットに基づいて前記印刷装置に印刷を実行させてもよい。波長ごとに異なる重みを課しつつ前記ターゲット分光反射率への近似性を評価する評価値に基づいて前記色材量セットを予測するため、近似性の重要度が高い波長域を優先させて近似させることができる。
【００１０】
また、前記重みの好適な一例として、人間の目の分光感度特性に基づいて設定されるものを適用してもよい。このようにすることにより、人間の分光感度が敏感な波長について優先的に分光反射率を近似させることができ、視覚的に再現精度のよい印刷結果を得ることができる。より具体的な例として、前記重みが３刺激値に対応する各等色関数の線形結合に基づいて設定されるようにしてもよい。このようにすれば、３刺激値に対応する各等色関数に対応する波長域を総合的に重視することが可能な前記重みを設定することができる。
【００１１】
さらに、前記重みが前記ターゲット分光反射率に基づいて設定されてもよい。例えば、前記ターゲット分光反射率が強いスペクトルを有している波長域については、分光反射率の近似性が最終的に視覚に大きく影響すると考えられるため、当該波長域を優先させて近似させることが望ましい。また、前記重みが所定の光源の分光エネルギー分布に基づいて設定されてもよい。所定の光源の分光エネルギー分布に応じて前記重みを設定することにより、光源が強いスペクトルを有している波長域について優先させて近似させることでき、当該光源における視覚的再現性を向上させることができる。また、単一光源のみならず複数光源での再現性を総合的に考慮して、前記重みを複数光源の分光エネルギーの線形結合に基づいて設定されるようにしてもよい。
【００１２】
さらに、本発明の技術的思想は、具体的な印刷制御装置にて具現化されるのみならず、その方法としても具現化することができる。すなわち、上述した印刷制御装置が行う各手段に対応する工程を有する方法としても本発明を特定することができる。むろん、上述した印刷制御装置がプログラムを読み込んで上述した各手段を実現する場合には、当該各手段に対応する機能を実行させるプログラムや当該プログラムを記録した各種記録媒体においても本発明の技術的思想が具現化できることは言うまでもない。なお、本発明の印刷制御装置は、単一の装置のみならず、複数の装置によって分散して存在可能であることはいうまでもない。例えば、印刷制御装置が有する各手段が、パーソナルコンピュータ上で実行されるプリンタドライバと、プリンタの双方において分散することも可能である。また、プリンタ等の印刷装置に本発明の印刷制御装置の各手段を包含させることも可能である。
【発明を実施するための最良の形態】
【００１３】
以下、下記の順序に従って本発明の実施形態を説明する。
１．印刷制御装置の構成：
２．印刷データ生成処理：
３．印刷制御処理：
３−１．１Ｄ−ＬＵＴ作成処理：
３−２．印刷制御データ生成処理：
４．分光プリンティングモデル：
５．変形例：
５−１：変形例１：
５−２：変形例２：
５−３：変形例３：
５−４：変形例４：
５−５：変形例５：
【００１４】
１．印刷制御装置の構成：
図１は、本発明の一実施形態にかかる印刷制御装置のハードウェア構成を示している。同図において、印刷制御装置は主にコンピュータ１０によって構成されており、コンピュータ１０はＣＰＵ１１とＲＡＭ１２とＲＯＭ１３とハードディスクドライブ（ＨＤＤ）１４と汎用インターフェイス（ＧＩＦ）１５とビデオインターフェイス（ＶＩＦ）１６と入力インターフェイス（ＩＩＦ）１７とバス１８とから構成されている。バス１８は、コンピュータ１０を構成する各要素１１〜１７の間でのデータ通信を実現するものであり、図示しないチップセット等によって通信が制御されている。ＨＤＤ１４には、オペレーティングシステム（ＯＳ）を含む各種プログラムを実行するためのプログラムデータ１４ａが記憶されており、当該プログラムデータ１４ａをＲＡＭ１２に展開しながらＣＰＵ１１が当該プログラムデータ１４ａに準じた演算を実行する。ＧＩＦ １５は、例えばＵＳＢ規格に準じたインターフェイスを提供するものであり、外部のプリンタ２０と分光反射率計３０をコンピュータ１０に接続させている。ＶＩＦ １６はコンピュータ１０を外部のディスプレイ４０に接続し、ディスプレイ４０に画像を表示するためのインターフェイスを提供する。ＩＩＦ １７はコンピュータ１０を外部のキーボード５０ａとマウス５０ｂに接続し、キーボード５０ａとマウス５０ｂからの入力信号をコンピュータ１０が取得するためのインターフェイスを提供する。
【００１５】
図２は、コンピュータ１０にて実行されるプログラムのソフトウェア構成を概略的なデータの流れとともに示している。同図において、コンピュータ１０では、おもにＯＳ Ｐ１と見本印刷アプリケーション（ＡＰＬ）Ｐ２とプリンタドライバ（ＰＤＶ）Ｐ３と分光反射率計ドライバ（ＭＤＶ）Ｐ４とディスプレイドライバ（ＤＤＶ）Ｐ５が実行されている。ＯＳ Ｐ１は、各プログラムが使用可能なＡＰＩのひとつとして画像機器インターフェイス（ＧＤＩ） Ｐ１ａとスプーラＰ１ｂを提供しており、ＡＰＬ Ｐ２の要求に応じてＧＤＩ Ｐ１ａが呼び出され、さらにＧＤＩ Ｐ１ａの要求に応じてＰＤＶ Ｐ３やＤＤＶ Ｐ５が呼び出される。ＧＤＩ Ｐ１ａはコンピュータ１０がプリンタ２０やディスプレイ４０等の画像出力装置における画像出力を制御するための汎用的な仕組みを提供し、一方のＰＤＶ Ｐ３やＤＤＶ Ｐ５はプリンタ２０やディスプレイ４０の機種固有の処理等を提供する。また、スプーラＰ１ｂは、ＡＰＬ Ｐ２やＰＤＶ Ｐ３やプリンタ２０の間に介在し、ジョブのコントロール等を実行する。ＡＰＬ Ｐ２は、見本チャートＳＣを印刷するためのアプリケーションプログラムであり、ＲＧＢビットマップ形式の印刷データＰＤを生成し、ＧＤＩ Ｐ１ａに対して当該印刷データＰＤを出力する。また、印刷データＰＤを生成するにあたっては、ＭＤＶ Ｐ４からターゲットの測色データＭＤを取得する。ＭＤＶ Ｐ４は、ＡＰＬ Ｐ２の要求に応じて分光反射率計３０を制御し、当該制御によって得られた分光反射率データＲＤをＡＰＬ Ｐ２に出力する。
【００１６】
ＡＰＬ Ｐ２が生成した印刷データＰＤはＧＤＩ Ｐ１ａやスプーラＰ１ｂを経由してＰＤＶ Ｐ３に出力され、ＰＤＶ Ｐ３が印刷データＰＤに基づいてプリンタ２０に出力可能な印刷制御データＣＤを生成する処理を実行する。ＰＤＶ Ｐ３が生成した印刷制御データＣＤはＯＳ Ｐ１が提供するスプーラＰ１ｂを介してプリンタ２０に出力され、プリンタ２０が当該印刷制御データＣＤに基づく動作を行うことにより見本チャートＳＣを印刷用紙上に印刷させる。以上においては、全体の処理の流れ概略的に説明したが、以下、フローチャートを用いて各プログラムＰ１〜Ｐ４が実行する処理を詳細に説明する。
【００１７】
２．印刷データ生成処理
図３は、ＡＰＬ Ｐ２が実行する印刷データ生成処理の流れを示している。図２に示すようにＡＰＬ Ｐ２はＵＩ部（ＵＩＭ）Ｐ２ａと測定制御部（ＭＣＭ）Ｐ２ｂと印刷データ生成部（ＰＤＧ）Ｐ２ｃとから構成されており、これらの各モジュールＰ２ａ，Ｐ２ｂ，Ｐ２ｃが図３に示す各ステップを実行する。ステップＳ１００においては、ＵＩＭ Ｐ２ａがＧＤＩ Ｐ１ａおよびＤＤＶ Ｐ５を介して、見本チャートＳＣを印刷させる旨の印刷指示を受付けるためのＵＩ画面を表示させる。前記ＵＩ画面においては、見本チャートＳＣのテンプレートを示す表示が設けられている。
【００１８】
図４は、前記ＵＩ画面の一例を示している。同図において、前記テンプレートＴＰが表示されており、当該テンプレートＴＰにはカラーパッチをレイアウトするための１２個の枠ＦＬ１〜ＦＬ１２が設けられている。前記ＵＩ画面には各枠ＦＬ１〜ＦＬ１２をマウス５０ｂのクリックによって選択することが可能となっており、枠ＦＬ１〜ＦＬ１２をクリックすると分光反射率測定を開始させるか否かを指示するための選択ウィンドウＷが表示される。また、前記ＵＩ画面においては、見本チャートＳＣの印刷を実行させるか否かを指示するためのボタンＢも設けられている。ステップＳ１１０においては、ＵＩＭ Ｐ２ａがマウス５０ｂによる各枠ＦＬ１〜ＦＬ１２のクリックが検出し、検出された場合にはステップＳ１２０にて分光反射率測定を開始させるか否かを指示するための選択ウィンドウＷを表示させる。ステップＳ１３０においては選択ウィンドウＷにおけるマウス５０ｂのクリックを検出し、キャンセルがクリックされた場合にはステップＳ１１０に戻る。一方、分光反射率測定実行がクリックされた場合には、ステップＳ１４０においてＭＣＭ Ｐ２ｂがＭＤＶ Ｐ４を介して分光反射率計３０にターゲットＴＧの分光反射率Ｒ（λ）としてターゲット分光反射率Ｒt（λ）の測定を実行させ、当該ターゲット分光反射率Ｒt（λ）を格納した分光反射率データＲＤを取得する。
【００１９】
ステップＳ１４０におけるターゲット分光反射率Ｒt（λ）の測定が完了すると、最も標準的な光源であるＤ６５光源を照射したときのＣＩＥＬＡＢ色空間における色彩値（Ｌ*ａ*ｂ*値）を算出する。そして、当該Ｌ*ａ*ｂ*値を所定のＲＧＢプロファイルを使用してＲＧＢ値に変換し、当該ＲＧＢ値を表示用ＲＧＢ値として取得する。なお、ＲＧＢプロファイルは絶対色空間としてのＣＩＥＬＡＢ色空間と本実施形態のＲＧＢ色空間との等色関係を規定したプロファイルであり、例えばＩＣＣプロファイルを使用することができる。
【００２０】
図５は、ステップＳ１４０において、分光反射率データＲＤから表示用ＲＧＢ値を算出する様子を模式的に示している。例えばターゲットＴＧについてターゲット分光反射率Ｒt（λ）を測定した結果、図示するようターゲット分光反射率Ｒt（λ）の分布を示す分光反射率データＲＤが得られたとする。なお、ターゲットＴＧは、分光的な再現の目標とする物体表面を意味し、例えば他の印刷装置や塗装装置、絵の具の塗布等で形成した人工的物体表面や、自然物の表面等が該当する。一方、Ｄ６５光源は図示するような可視波長域において不均一な分光エネルギーＰ（λ）の分布を有しており、ターゲットＴＧにＤ６５光源を照射したときの各波長の反射光の分光エネルギーは、ターゲット分光反射率Ｒt（λ）と分光エネルギーＰ（λ）を各波長について掛け合わせた値となる。さらに、反射光の分光エネルギーのスペクトルに対して人間の分光感度特性に応じた等色関数ｘ（λ），ｙ（λ），ｚ（λ）をそれぞれ畳み込み積分し、係数ｋによって正規化することにより、３刺激値Ｘ，Ｙ，Ｚを得る。以上を数式で表すと下記の（１）式となる。
【数１】
0004924414.tif 000002
【００２１】
３刺激値Ｘ，Ｙ，Ｚを所定の変換式によって変換することにより、ターゲットＴＧにＤ６５光源を照射したときの色を示すＬ*ａ*ｂ*値を得ることができ、さらにＲＧＢプロファイルを使用することにより、表示用ＲＧＢ値を得ることができる。ステップＳ１４５においては、テンプレートＴＰにおいてクリックされた枠ＦＬ１〜ＦＬ１２を前記表示用ＲＧＢ値で塗りつぶした表示に更新する。これにより、標準的な光源であるＤ６５光源でのターゲットＴＧの色をＵＩ画面にて感覚的に把握することが可能となる。ステップＳ１４５が完了すると、ステップＳ１５０において固有のインデックスを生成するとともに、当該インデックスと、ステップＳ１１０にてクリックされた枠ＦＬ１〜ＦＬ１２の位置情報と表示用ＲＧＢ値を分光反射率データＲＤに対応付けてＲＡＭ１２に記憶する。ステップＳ１５０が完了すると、ステップＳ１１０に戻り、ステップＳ１２０〜Ｓ１５０を繰り返して実行する。これにより、他の枠ＦＬ１〜ＦＬ１２を選択し、他の枠ＦＬ１〜ＦＬ１２について他の色の物体表面が再現されたターゲットＴＧのターゲット分光反射率Ｒt（λ）の測定を行うことができる。
【００２２】
本実施形態においては、それぞれ異なる１２種類のターゲットＴＧ１〜ＴＧ１２が用意されており、ターゲットＴＧ１〜ＴＧ１２のそれぞれについてのターゲット分光反射率Ｒt（λ）が分光反射率測データＲＤとして取得されるものとする。従って、ステップＳ１５０においては、各枠ＦＬ１〜ＦＬ１２についての分光反射率測データＲＤと固有のインデックスを対応付けたデータが順次ＲＡＭに記憶されていくこととなる。なお、インデックスはそれぞれの値が固有となるように生成されればよく、インクリメントによって生成してもよいし、重複しない乱数によって生成してもよい。
【００２３】
ステップＳ１１０において、各枠ＦＬ１〜ＦＬ１２のクリックが検出されない場合には、ステップＳ１６０にて見本チャートＳＣの印刷を実行させる旨のボタンＢのクリックを検出し、検出されない場合にはステップＳ１１０に戻る。一方、見本チャートＳＣの印刷を実行させる旨のボタンＢのクリックが検出された場合には、ステップＳ１７０にてＰＤＧ Ｐ２ｃが印刷データＰＤを生成する。
【００２４】
図６は、印刷データＰＤの構成を模式的に示している。同図において、印刷データＰＤはドットマトリクス状に配列した多数の画素によって構成されており、各画素が４バイト（８ビット×４）の情報を有している。印刷データＰＤは、図４に示したテンプレートＴＰと同様の画像を示しており、テンプレートＴＰの各枠ＦＬ１〜ＦＬ１２に対応する領域以外の画素は、テンプレートＴＰに対応する色のＲＧＢ値を有している。ＲＧＢ各チャネルの階調値はそれぞれ８ビット（２５６階調）によって表現され、上述した４バイトのうち３バイトがＲＧＢ値を格納するために使用される。例えば、テンプレートＴＰの各枠ＦＬ１〜ＦＬ１２以外の色が（Ｒ，Ｇ，Ｂ）＝（１２８，１２８，１２８）の一様な中間グレーで表される場合、印刷データＰＤにおける各枠ＦＬ１〜ＦＬ１２に対応する領域以外の画素は（Ｒ，Ｇ，Ｂ）＝（１２８，１２８，１２８）の色情報を有することとなる。なお、残りの１バイトは使用されない。
【００２５】
一方、テンプレートＴＰの各枠ＦＬ１〜ＦＬ１２に対応する画素も４バイトの情報を有しており、通常、ＲＧＢ値が格納される３バイトを使用してインデックスを格納する。このインデックスは、ステップＳ１５０にて各枠ＦＬ１〜ＦＬ１２ごとに生成した固有のものであり、ＰＤＧ Ｐ２ｃはインデックスをＲＡＭ１２から取得し、各枠ＦＬ１〜ＦＬ１２に対応する画素に対応するインデックスを格納する。このようにインデックスをＲＧＢ値の代わりに格納した各枠ＦＬ１〜ＦＬ１２に対応する画素については、残りの１バイトを使用してインデックスが格納された旨のフラグを立てる。これにより、各画素がＲＧＢ値を格納しているか、インデックスを格納しているかを判別することができる。本実施形態では、インデックスを格納するために３バイトを使用することができるため、３バイト以下の情報量で表現できるインデックスをステップＳ１５０で生成しておく必要がある。以上のようにしてビットマップ形式の印刷データＰＤが生成できると、ステップＳ１８０において、ＰＤＧ Ｐ２ｃがインデックステーブルＩＤＢを生成する。
【００２６】
図７は、インデックステーブルＩＤＢの一例を示している。同図において、各枠ＦＬ１〜ＦＬ１２に対応して生成された固有のインデックスのそれぞれに対して、測定によって得られたターゲット分光反射率Ｒt（λ）およびＤ６５光源におけるＬ*ａ*ｂ*値に対応する表示用ＲＧＢ値が格納されている。インデックステーブルＩＤＢの生成が完了すると、印刷データＰＤはＧＤＩ Ｐ１ａやスプーラＰ１ｂを経由してＰＤＶ Ｐ３に出力される。印刷データＰＤは、外形上、通常のＲＧＢビットマップ形式と変わらないため、ＯＳ Ｐ１が提供するＧＤＩ Ｐ１ａやスプーラＰ１ｂにおいても通常の印刷ジョブと同様に処理することができる。一方、インデックステーブルＩＤＢは、直接、ＰＤＶ Ｐ３に出力される。なお、本実施形態においては、インデックステーブルＩＤＢを新たに生成するようにしたが、既存のインデックステーブルＩＤＢにインデックスとターゲット分光反射率Ｒt（λ）と表示用ＲＧＢ値の新たな対応関係を追記するようにしてもよい。また、以上の印刷データ生成処理と後述する印刷制御処理は必ずしも同一の装置において連続して実行する必要はなく、印刷データ生成処理と印刷制御処理を例えばＬＡＮやインターネット等の通信回線によって接続された複数のコンピュータ上にて個別に実行してもよい。
【００２７】
３．印刷制御処理
図８は、ＰＤＶ Ｐ３が実行する印刷制御処理の全体的な流れを示している。図２に示すようにＰＤＶ Ｐ３は、１Ｄ−ＬＵＴ生成部（ＬＵＧ）Ｐ３ａと印刷制御データ生成部（ＣＤＧ）Ｐ３ｂとから構成されており、図８に示す１Ｄ−ＬＵＴ生成処理（ステップＳ２００）をＬＵＧ Ｐ３ａが担当し、一方の印刷制御データ生成処理（ステップＳ３００）をＣＤＧ Ｐ３ｂが担当する。１Ｄ−ＬＵＴ生成処理は印刷制御データ生成処理に先行して行われてもよいし、１Ｄ−ＬＵＴ生成処理と印刷制御データ生成処理を並行して行うようにしてもよい。
【００２８】
３−１．１Ｄ−ＬＵＴ生成処理
図９は、１Ｄ−ＬＵＴ生成処理の流れを示している。図２に示すようにＬＵＧ Ｐ３ａは、インク量セット算出モジュール（ＩＣＭ）Ｐ３ａ１と分光反射率予測モジュール（ＲＰＭ）Ｐ３ａ２と評価値算出モジュール（ＥＣＭ）Ｐ３ａ３とＬＵＴ出力モジュール（ＬＯＭ） Ｐ３ａ４とから構成されている。ステップＳ２１０においては、ＩＣＭ Ｐ３ａ１がインデックステーブルＩＤＢを取得する。ステップＳ２２０においては、インデックステーブルＩＤＢから一つのインデックスを選択し、当該インデックスに対応付けられている分光反射率データＲＤを取得する。ステップＳ２３０においては、ＩＣＭ Ｐ３ａ１が前記分光反射率データＲＤが示すターゲット分光反射率Ｒt（λ）と同様の分光反射率Ｒ（λ）が再現可能なインク量セットを算出する処理を行う。その際に、上述したＲＰＭ Ｐ３ａ２とＥＣＭ Ｐ３ａ３を使用する。
【００２９】
ここで本発明のインクセットＩｎｋについて説明する。インクセットＩｎｋは、Ｃ１，Ｃ２，Ｃ３，Ｍ１，Ｍ２，Ｍ３，Ｙ１，Ｙ２，Ｙ３，Ｋ，ｌｃ，ｌｍの各色で構成されている。このインクセットＩｎｋにおいて、Ｃ１とＣ２とＣ３は、所定の標準光源下における色彩値（Ｌ＊ａ＊ｂ＊等、絶対色空間における座標）がいずれも略等しく（何れもシアンと見做せる）、分光反射率分布が互いに異なるものが採用されている。同様に、Ｍ１とＭ２とＭ３は所定の標準光源下において色彩値がＭでありつつ分光反射率分布が互いに異なっており、Ｙ１とＹ２とＹ３は所定の標準光源下において色彩値がＹでありつつ分光反射率分布が互いに異なっている。以下、Ｃ１とＣ２とＣ３の関係、Ｍ１とＭ２とＭ３の関係、Ｙ１とＹ２とＹ３の関係を、等色関係とし、等色関係にあるインクの組を等色インクセットと呼ぶことにする。
【００３０】
図１０は、Ｃにおける等色インクセットの分光反射率分布の一例を示すグラフであり、図１１はＭにおける等色インクセットの分光反射率分布の一例を示すグラフであり、図１２はＹにおける等色インクセットの分光反射率分布の一例を示すグラフである。図１０〜図１２に示すように、等色インクセットを構成する各インクは、各インク色に１つのインクを使用する場合に比べて、各々幅狭のピークを持つ分光反射率分布を示している。また、等色インクセットを構成するインク間で分光反射率分布に現れるピーク位置が互いに異なるものの、Ｌ＊ａ＊ｂ＊値が互いにほぼ等しくなっている。したがって、メタメリズムを考慮しなくてもよい印刷においては、等色インクセットを構成する各インクは同一色として扱うことができる。
【００３１】
以上の説明においては、ＣＭＹが等色インクセットを有するように記載してあるが、無論、ＣＭＹの何れか一色のみが等色インクセットを有するように構成してもよいし、ＣＭＹＫｌｃｌｍの中の任意の組合せが等色インクセットで構成されるようにしてもよい。また、等色インクセットを構成する各インクはＤ６５光源に照射された場合に同一の色彩値を呈する例について記載してあるが、等色インクセットを構成する各インクがどのような光源下で等色関係になるかは適宜選択可能であり、標準昼光系のＤ５０光源，Ｄ５５光源，白熱電球系のＡ光源、蛍光ランプ系のＦ１１光源等、様々な光源を採用可能である。無論、このような標準光源のみならず、等色関係を維持したい光源は任意に選択可能である。
【００３２】
図１３は、前記分光反射率データＲＤが示すターゲット分光反射率Ｒt（λ）と同様の分光反射率Ｒ（λ）が再現可能なインク量セットを算出する処理の流れを模式的に示している。ＲＰＭ Ｐ３ａ２は、ＩＣＭ Ｐ３ａ１からのインク量セットφ（ｄＣ１，ｄＣ２，ｄＣ３，ｄＭ１，ｄＭ２，ｄＭ３，ｄＹ１，ｄＹ２，ｄＹ３，ｄＫ，ｄｌｃ，ｄｌｍ）の入力に応じて、当該インク量セットφに基づいてプリンタ２０が所定の印刷用紙にインクを吐出させたときの分光反射率Ｒ（λ）を予測し、当該分光反射率Ｒ（λ）を予測分光反射率Ｒs（λ）としてＥＣＭ Ｐ３ａ３に出力する。このように等色インクセットを用いることにより、従来のような幅広のピークを持ったインクでターゲット分光反射率に対するフィッティングを行った場合に比べて、高精度なフィッティングが可能となり、よりターゲット分光反射率Ｒｔ（λ）との一致率の高い予測分光反射率Ｒｓ（λ）を再現可能となる。すなわち以下に説明を行う評価において良好な結果を得ることができる。
【００３３】
ＥＣＭ Ｐ３ａ３は、分光反射率データＲＤが示すターゲット分光反射率Ｒt（λ）と予測分光反射率Ｒs（λ）の差分Ｄ（λ）を各波長λについて算出し、各波長λごとに重みが課せられた重み関数ｗ（λ）を当該差分Ｄ（λ）に乗算する。この値の二乗平均の平方根を評価値Ｅ（φ）として算出する。以上の計算を数式で表すと下記の（２）式のように表すことができる。
【数２】
0004924414.tif 000003
【００３４】
前記の（２）式において、Ｎは波長λの有限区分数を意味する。前記の（２）式において、評価値Ｅ（φ）が小さければ小さいほど、各波長λにおけるターゲット分光反射率Ｒt（λ）と予測分光反射率Ｒs（λ）の差が少ないということができる。すなわち、評価値Ｅ（φ）が小さければ小さいほど、入力したインク量セットφによってプリンタ２０が印刷したときに記録媒体上にて再現される分光反射率Ｒ（λ）と、対応するターゲットＴＧから得られたターゲット分光反射率Ｒt（λ）が近似しているということができる。さらに、前記の（１）式によれば、光源の変動に応じてインク量セットφによってプリンタ２０が印刷したときの記録媒体と対応するターゲットＴＧが示す絶対的な色彩値は双方とも変動するものの、分光反射率Ｒ（λ）が近似すれば光源の変動に拘わらず相対的には同じ色に知覚されるということができる。従って、評価値（φ）が小さくなるインク量セットφによれば、あらゆる光源においてターゲットＴＧと同じ色に知覚される印刷結果を得ることができるということができる。
【００３５】
また、本実施形態において、重み関数ｗ（λ）は下記の（３）式のものを使用する。
【数３】
0004924414.tif 000004
前記の（３）式においては、等色関数ｘ（λ），ｙ（λ），ｚ（λ）を加算することにより、重み関数ｗ（λ）が定義されている。なお、前記の（３）式の右辺全体に所定の係数を乗算して、重み関数ｗ（λ）の値の範囲を正規化してもよい。前記の（１）式によれば、等色関数ｘ（λ），ｙ（λ），ｚ（λ）が大きい波長域ほど、色彩値（Ｌ*ａ*ｂ*値）に大きく影響するということができる。従って、等色関数ｘ（λ），ｙ（λ），ｚ（λ）を加算した重み関数ｗ（λ）を使用すれば、色への影響が大きい波長域を重視した二乗誤差が評価可能な評価値Ｅ（φ）を得ることができる。例えば、人間の目に知覚されない近紫外波長域においてはｗ（λ）が０となり、当該波長域における差分Ｄ（λ）は評価値Ｅ（φ）の増大に寄与しないこととなる。
【００３６】
すなわち、必ずしも全可視波長域においてターゲット分光反射率Ｒt（λ）と予測分光反射率Ｒs（λ）との差が小さくなくても、人間の目に特に強く知覚される波長域においてターゲット分光反射率Ｒt（λ）と予測分光反射率Ｒs（λ）とが似ていれば小さい値の評価値Ｅ（φ）を得ることができ、人間の目に知覚に即した分光反射率Ｒ（λ）の近似性の指標として評価値Ｅ（φ）を使用することができる。算出された評価値Ｅ（φ）はＩＣＭ Ｐ３ａ１に返される。すなわち、ＩＣＭＰ ３ａ１が任意のインク量セットφをＲＰＭ Ｐ３ａ２とＥＣＭ Ｐ３ａ３に出力することにより、最終的に評価値Ｅ（φ）がＩＣＭ Ｐ３ａ１に返される構成となっている。ＩＣＭ Ｐ３ａ１は任意のインク量セットφに対応して評価値Ｅ（φ）を得ることを繰り返し実行することにより、目的関数としての評価値Ｅ（φ）が極小化するようなインク量セットφの最適解を算出する。この最適解を算出する手法としては、様々な最適化手法を用いることができるが、例えば勾配法といった非線形最適化手法を用いることが望ましい。
【００３７】
図１４は、ステップＳ２３０においてインク量セットφが最適化されていく様子を模式的に示している。同図において、インク量セットφが最適化されていくにつれて、インク量セットφで印刷を行った場合の予測分光反射率Ｒs（λ）がターゲット分光反射率Ｒt（λ）に近づいていく。また、重み関数ｗ（λ）を使用することにより、等色関数ｘ（λ），ｙ（λ），ｚ（λ）が大きい波長域ほど、予測分光反射率Ｒs（λ）のターゲット分光反射率Ｒt（λ）への拘束が強くなっており、予測分光反射率Ｒs（λ）がターゲット分光反射率Ｒt（λ）の差が小さくなっている。このように、等色関数ｘ（λ），ｙ（λ），ｚ（λ）が大きく、視覚に大きく影響する波長域について優先的に予測分光反射率Ｒs（λ）をターゲットＴＧのターゲット分光反射率Ｒt（λ）に拘束することができるため、任意の光源を照射したときの見た目が近くなるようなインク量セットφを算出することができる。以上により、いずれの光源においてもターゲットＴＧと似たような見た目をプリンタ２０にて再現させることが可能なインク量セットφを算出することができる。なお、最適化の終了条件は、インク量セットφ更新の繰り返し回数としてもよいし、評価値Ｅ（φ）の閾値としてもよい。
【００３８】
以上のようにして、ステップＳ２３０においてＩＣＭ Ｐ３ａ１がターゲットＴＧと同様の分光反射率Ｒ（λ）が再現可能なインク量セットφを算出すると、ステップＳ２４０においてインデックステーブルＩＤＢに記述されたインデックスのすべてがステップＳ２２０にて選択されたか否かを判定し、すべて選択されていない場合にはステップＳ２２０に戻り、次のインデックスを選択する。このようにすることにより、すべてのインデックスについてターゲットＴＧと同様の色が再現可能なインク量セットφを算出することができる。すなわち、印刷データ生成処理（図２）のステップＳ１４０において測色を行ったすべてのターゲットＴＧ１〜ＴＧ１２についてターゲットＴＧ１〜ＴＧ１２と同様の分光反射率Ｒ（λ）が再現可能なインク量セットφを算出することができる。ステップＳ２４０において、すべてのインデックスについて最適なインク量セットφが算出されたことが判定されると、ステップＳ２５０において、ＬＯＭ Ｐ３ａ４が１Ｄ−ＬＵＴを生成し、当該１Ｄ−ＬＵＴをＣＤＧ Ｐ３ｂに出力する。
【００３９】
図１５は、１Ｄ−ＬＵＴの一例を示している。同図において、各インデックスに対応して最適なインク量セットφが格納されている。すなわち、各ターゲットＴＧ１〜ＴＧ１２について、各ターゲットＴＧ１〜ＴＧ１２と似たような見た目をプリンタ２０にて再現させることが可能なインク量セットφ（ｄＣ１，ｄＣ２，ｄＣ３，ｄＭ１，ｄＭ２，ｄＭ３，ｄＹ１，ｄＹ２，ｄＹ３，ｄＫ，ｄｌｃ，ｄｌｍ）を記述した１Ｄ−ＬＵＴを用意することができる。１Ｄ−ＬＵＴをＣＤＧ Ｐ３ｂに出力すると、１Ｄ−ＬＵＴ生成処理が完了し、次の印刷制御データ生成処理（ステップＳ３００）を実行させる。
【００４０】
３−２．印刷制御データ生成処理
図１６は、印刷制御データ生成処理の流れを示している。図２に示すようにＣＤＧ Ｐ３ｂは、モード判別モジュール（ＭＩＭ）Ｐ３ｂ１とインデックス分版モジュール（ＩＳＭ）Ｐ３ｂ２とＲＧＢ分版モジュール（ＣＳＭ）Ｐ３ｂ３とハーフトーンモジュール（ＨＴＭ）Ｐ３ｂ４とラスタ化モジュール（ＲＴＭ）Ｐ３ｂ５とインク残量監視モジュール（ＩＷＭ）Ｐ３ｂ６と等色分版モジュール（ＭＳＭ）Ｐ３ｂ７とから構成されている。ステップＳ３１０においては、モード判別モジュール（ＭＩＭ）Ｐ３ｂ１が印刷データＰＤを取得する。ステップＳ３２０において、ＭＩＭ Ｐ３ｂ１は印刷データＰＤから一の画素を選択する。ステップＳ３３０において、ＭＩＭ Ｐ３ｂ１は当該選択した画素にインデックスが格納された旨のフラグが立っているか否かを判定する。当該フラグが立っていないと判定された場合には、ステップＳ３４０にてＣＳＭ Ｐ３ｂ３が３Ｄ−ＬＵＴを参照して、当該画素についての色変換（分版）を実行する。
【００４１】
図１７は、３Ｄ−ＬＵＴの一例を示している。同図において、３Ｄ−ＬＵＴはＲＧＢ値とインク量セットψ（ｄＣ，ｄＭ，ｄＹ，ｄＫ，ｄｌｃ，ｄｌｍ）との対応関係が色空間における複数の代表的な座標について記述されたテーブルであり、ＣＳＭ Ｐ３ｂ３は３Ｄ−ＬＵＴを参照して当該画素が有するＲＧＢ値に対応するインク量セットψを取得する。その際に、３Ｄ−ＬＵＴに直接記述されていないＲＧＢ値については補間演算を行うことにより、対応するインク量セットψを取得する。なお、３Ｄ−ＬＵＴの作成方法として、特開２００６−８２４６０号公報等を採用することができる。当該公報においては、特定光源における色の再現性や、再現色の階調性や、粒状性や、再現色の光源非依存性や、ガマットや、インクデューティが総合的に良好となる３Ｄ−ＬＵＴが作成される。
【００４２】
３Ｄ−ＬＵＴによる色変換が行われた画素については、ステップＳ３４５においてＭＳＭ Ｐ３ｂ７が等色インクへの分版（以下、等色分版と称する。）を行い、インク量セットψ（ｄＣ，ｄＭ，ｄＹ，ｄＫ，ｄｌｃ，ｄｌｍ）からインク量セットψ（ｄＣ１，ｄＣ２，ｄＣ３，ｄＭ１，ｄＭ２，ｄＭ３，ｄＹ１，ｄＹ２，ｄＹ３，ｄＫ，ｄｌｃ，ｄｌｍ）を得る。なお、３Ｄ−ＬＵＴによって色変換が行われた画素については、メタメリズムを考慮する必要が無いため、ＣＭＹで指定された画素の印刷にはどの等色インクを使用してもよい。したがって、各等色インクの残量に応じて、等色インクをランダムに使用したり、特定の等色インクが優先的に使用されるような分版を行ったりできる。
【００４３】
そこで等色分版を行うにあたり、ＭＳＭ Ｐ３ｂ７はＩＷＭ Ｐ３ｂ６からインク残量情報を取得し、インク残量に応じた各等色インクへの分版割合を行うようにする。例えば、Ｃ１の残量が４０ｍｌ、Ｃ２の残量が３０ｍｌ、Ｃ３の残量が２０ｍｌであれば、インク残量の多いＣ１を優先的に使用するような分版割合に決定する。分版割合の具体的な一例としては、ｄＣ１：ｄＣ２：ｄＣ３＝３：４：６の割合で分版を行うことでＣ１，Ｃ２，Ｃ３が同時に用尽されるようにしたり、各インク残量が等しくなるまで残量の多いインクを優先的に使うようにしたりすることが考えられる。
【００４４】
インク残量情報を出力するＩＷＭ Ｐ３ｂ６は、プリンタ２０のインク吐出状況を監視することにより、各インクの残量を把握している。インク吐出状況を監視するには、例えば、プリンタ２０が備えるセンサによりインク吐出ノズルのインク吐出を監視したり、各印刷ヘッドの使用時間、各インク吐出ノズルの吐出回数、各印刷ヘッドにおける印字枚数、等から算出したりすることができる。プリンタ２０のセンサには、例えば、インクを吐出するノズルを挟んで設けられた発光素子と受光素子により、その間をインク滴が通過することによる光の強弱の変化を検出し、各ノズルの動作を確認するもの等で実現可能である。
【００４５】
また、メタメリズムの防止や濃淡の逆転を防止する意味では、一度の印刷においては各等色インクの使用比率を一定にすると好適である。すなわち、印刷分版の途中で、各等色インクの比率が変化してしまうと、同じ色彩値であるはずのドット同士が光源によっては互いに異なった色に見えてしまったり、色彩の濃淡が逆転してしまったりする可能性があるためである。
【００４６】
以上のように、ステップＳ３４０において３Ｄ−ＬＵＴでＣＭＹまでの分版を行い、続くステップＳ３４５で等色インクへの分版を行うことにより、３Ｄ−ＬＵＴが大きくなりすぎずに済む。すなわち等色インクまで含めた分版を３Ｄ−ＬＵＴで実現しようとすると、各一次色上に５グリッドずつの３Ｄ−ＬＵＴであれば、ＣＭＹＫｌｃｌｍの場合には５６の格子点についての情報が必要であるのに対し、前述のようにＣＭＹ各色がそれぞれ３つの等色インクを備えるようにした場合には５１２の格子点についての情報が必要となる。
【００４７】
一方、ステップＳ３３０において、選択した画素にインデックスが格納された旨のフラグが立っていると判定された場合には、ステップＳ３５０にてＩＳＭ Ｐ３ｂ２が１Ｄ−ＬＵＴを参照して、当該画素についての色変換（分版）を実行する。すなわち、インデックスが格納された旨のフラグが立っている画素から、インデックスを取得するとともに、１Ｄ−ＬＵＴにて当該インデックスに対応付けられているインク量セットφを取得する。ステップＳ３４０またはステップＳ３５０のいずれかにおいて、当該画素についてのインク量セットφが取得できると、ステップＳ３６０においてすべての画素についてインク量セットφが取得できたか否かを判定する。ここでインク量セットφが未取得の画素が残っている場合には、ステップＳ３２０に戻り次の画素を選択する。以上のステップＳ３２０のように、分光反射率分布が互いに異なる等色インクを使い分けてターゲットの分光反射率を表現することにより、ターゲットのメタメリズムをより正確に反映したインク量セットφが作成される。
【００４８】
以上の処理を繰り返して実行することにより、すべての画素についてインク量セットφを取得することができる。すべての画素についてインク量セットφが取得できると、すべての画素がインク量セットφで表現された印刷データＰＤに変換されたこととなる。以上のように各画素について１Ｄ−ＬＵＴと３Ｄ−ＬＵＴのいずれを使用するかを判定することにより、インデックスが格納された枠Ｆ１〜Ｆ１２に対応する画素については、各光源において各ターゲットＴＧ１〜ＴＧ１２に近い色が再現可能なインク量セットφを取得することができるとともに、ＲＧＢ値が格納された画素については３Ｄ−ＬＵＴの作成指針（例えば、粒状性を重視する等。）に基づく色再現が可能なインク量セットφを取得させることができる。
【００４９】
ステップＳ３７０においては、各画素がインク量セットφで表現された印刷データＰＤをＨＴＭ Ｐ３ｂ４が取得し、ハーフトーン処理を実行する。ＨＴＭ Ｐ３ｂ４はハーフトーン処理をするにあたっては公知のディザ法や誤差拡散法等を使用することができる。ハーフトーン処理が完了した印刷データＰＤにおいては、各画素が各インクを吐出させるか否を示す吐出信号を有している。ステップＳ３８０では、ハーフトーン処理が完了した印刷データＰＤをＲＴＭ Ｐ３ｂ５が取得し、当該印刷データＰＤにおける吐出信号をプリンタ２０が有する印刷ヘッドの各走査パスおよび各ノズルに割り振る処理を実行する。以上によりプリンタ２０に出力可能な印刷制御データＣＤが生成でき、プリンタ２０の制御に必要な信号を添付した印刷制御データＣＤをスプーラＰ１ｂおよびプリンタ２０に出力する。これによりプリンタ２０が印刷用紙上にインクを吐出して、見本チャートＳＣを形成する。
【００５０】
以上のようにして印刷用紙上に形成された見本チャートＳＣの枠ＦＬ１〜ＦＬ１２に対応する領域において、各ターゲットＴＧ１〜ＴＧ１２のターゲット分光反射率Ｒt（λ）を再現することができる。すなわち、枠ＦＬ１〜ＦＬ１２に対応する領域は、ターゲットＴＧ１〜ＴＧ１２の複数光源下での色に追従するようなインク量セットφで印刷されているため、各光源下においてターゲットＴＧ１〜ＴＧ１２と似たような色を再現することができる。例えば、見本チャートＳＣを室内で視認したときの各枠ＦＬ１〜ＦＬ１２に対応する領域の色はターゲットＴＧ１〜ＴＧ１２を室内で視認したときの色を再現するし、見本チャートＳＣを室外で視認したときの各枠ＦＬ１〜ＦＬ１２に対応する領域の色もターゲットＴＧ１〜ＴＧ１２を室外で視認したときの色を再現することができる。
【００５１】
なお、究極的には、ターゲットＴＧ１〜ＴＧ１２と完全に同じ分光反射率Ｒ（λ）を有する見本チャートＳＣを再現すれば、いかなる光源においてもターゲットＴＧ１〜ＴＧ１２と同様の色を再現することができる。しかしながら、プリンタ２０が使用可能なインク（色材の種類）はＣＭＹＫｌｃｌｍに限られているため、ターゲットＴＧ１〜ＴＧ１２と完全に同じ分光反射率Ｒ（λ）が再現可能なインク量セットφを求めることは実質的に不可能である。また、知覚色に影響しない波長域についてもターゲットＴＧ１〜ＴＧ１２と同様な分光反射率Ｒ（λ）が再現可能なインク量セットφを求めても、視覚的な再現精度の実現においては無駄となる。これに対して、本発明では、等色関数ｘ（λ），ｙ（λ），ｚ（λ）に基づく重み付けを行った評価値Ｅ（φ）を利用してターゲット分光反射率Ｒt（λ）への近似性を評価しているため、視覚的に十分な精度が達成可能なインク量セットφを求めることができる。
【００５２】
一方、印刷用紙上に形成された見本チャートＳＣの枠ＦＬ１〜ＦＬ１２に対応する領域においては、上述した３Ｄ−ＬＵＴに基づいたインク量セットφによって印刷がなされることとなる。そのため、当該領域についての印刷パフォーマンスは３Ｄ−ＬＵＴに基づくものとなる。上述したとおり本実施形態においては枠ＦＬ１〜ＦＬ１２以外の領域は中間グレーの一様な画像を示すが、当該領域において３Ｄ−ＬＵＴが目標とする印刷パフォーマンスを満足させることができる。すなわち、再現色の階調性や、粒状性や、再現色の光源非依存性や、ガマットや、インクデューティが総合的に良好となる印刷を実現することができる。
【００５３】
４．分光プリンティングモデル
図１８は、本実施形態のプリンタ２０の印刷方式を模式的に示している。なお等色インクを含めた分光プリンティングモデルとするには、以下のモデルに各等色インクを一次色として加えてやればよい。図１８において、プリンタ２０は、ＣＭＹＫｌｃｌｍのインクごとに複数のノズル２１ａ，２１ａ・・・を備えた印刷ヘッド２１を備えており、ノズル２１ａ，２１ａ・・・が吐出するＣＭＹＫｌｃｌｍのインクごとのインク量を上述したインク量セットψ（ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlm）によって指定された量とする制御が印刷制御データＣＤに基づいて行われる。各ノズル２１ａ，２１ａ・・・が吐出したインク滴は印刷用紙上において微細なドットとなり、多数のドットの集まりによってインク量セットψ（ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlm）に応じたインク被覆率の印刷画像が印刷用紙上に形成されることとなる。
【００５４】
ＲＰＭ Ｐ３ａ２が使用する予測モデル（分光プリンティングモデル）は、本実施形態のプリンタ２０で使用され得る任意のインク量セットψ（ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlm）で印刷を行った場合の分光反射率Ｒ（λ）を予測分光反射率Ｒs（λ）として予測するための予測モデルである。分光プリンティングモデルにおいては、インク量空間における複数の代表点について実際にカラーパッチを印刷し、その分光反射率Ｒ（λ）を分光反射率計によって測定することにより得られた分光反射率データベースＲＤＢを用意する。そして、この分光反射率データベースＲＤＢを使用したセル分割ユール・ニールセン分光ノイゲバウアモデル(Cellular Yule-Nielsen Spectral Neugebauer Model)による予測を行うことにより、正確に任意のインク量セットψ（ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlm）で印刷を行った場合の分光反射率Ｒ（λ）を予測する。
【００５５】
図１９は、分光反射率データベースＲＤＢを示している。同図に示すように分光反射率データベースＲＤＢはインク量空間（本実施形態では６次元であるが、図の簡略化のためＣＭ面のみ図示。）における複数の格子点のインク量セットψ（ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlm）について実際に印刷／測定をして得られた分光反射率Ｒ（λ）が記述されたルックアップテーブルとなっている。例えば、各インク量軸を分割する５グリッドの格子点を発生させる。ここでは５13個もの格子点が発生し、膨大な量のカラーパッチの印刷／測定をすることが必要となるが、実際にはプリンタ２０にて同時に搭載可能なインク数や同時に吐出可能なインクデューティの制限があるため、印刷／測定をする格子点の数は絞られることとなる。
【００５６】
さらに、一部の格子点のみ実際に印刷／測定をし、他の格子点については実際に印刷／測定を行った格子点の分光反射率Ｒ（λ）に基づいて分光反射率Ｒ（λ）を予測することにより、実際に印刷／測定を行うカラーパッチの個数を低減させてもよい。分光反射率データベースＲＤＢは、プリンタ２０が印刷可能な印刷用紙ごとに用意されている必要がある。厳密には、分光反射率Ｒ（λ）は印刷用紙上に形成されたインク膜（ドット）による分光透過率と印刷用紙の反射率によって決まるものであり、印刷用紙の表面物性（ドット形状が依存）や反射率の影響を大きく受けるからである。次に、分光反射率データベースＲＤＢを使用したセル分割ユール・ニールセン分光ノイゲバウアモデルによる予測を説明する。
【００５７】
ＲＰＭ Ｐ３ａ２は、ＩＣＭ Ｐ３ａ１の要請に応じて分光反射率データベースＲＤＢを使用したセル分割ユール・ニールセン分光ノイゲバウアモデルによる予測を実行する。この予測にあたっては、ＩＣＭ Ｐ３ａ１から予測条件を取得し、この予測条件を設定する。具体的には、印刷用紙やインク量セットψを印刷条件として設定する。例えば、光沢紙を印刷用紙として予測を行う場合には、光沢紙にカラーパッチを印刷することにより作成した分光反射率データベースＲＤＢが設定される。
【００５８】
分光反射率データベースＲＤＢの設定ができると、ＩＣＭ Ｐ３ａ１から入力されたインク量セットψ（ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlm）を分光プリンティングモデルに適用する。セル分割ユール・ニールセン分光ノイゲバウアモデルは、よく知られた分光ノイゲバウアモデルとユール・ニールセンモデルとに基づいている。なお、以下の説明では、説明の簡略化のためＣＭＹの３種類のインクを用いた場合のモデルについて説明するが、同様のモデルを本実施形態のＣＭＹＫｌｃｌｍを含む任意のインクセットを用いたモデルに拡張することは容易である。また、セル分割ユール・ニールセン分光ノイゲバウアモデルについては、Color Ｒes Appl 25, 4-19, 2000およびＲ Balasubramanian, Optimization of the spectral Neugebauer model for printer characterization, J. Electronic Imaging 8(2), 156-166 (1999)を参照。
【００５９】
図２０は、分光ノイゲバウアモデルを示す図である。分光ノイゲバウアモデルでは、任意のインク量セット（ｄc，ｄm，ｄy）で印刷したときの印刷物の予測分光反射率Ｒs（λ）は、以下の（４）式で与えられる。
【数４】
0004924414.tif 000005
【００６０】
ここで、ａiはｉ番目の領域の面積率であり、Ｒi（λ）はｉ番目の領域の分光反射率である。添え字ｉは、インクの無い領域（ｗ）と、シアンインクのみの領域（ｃ）と、マゼンタインクのみの領域（ｍ）と、イエローインクのみの領域（ｙ）と、マゼンタインクとイエローインクが吐出される領域（ｒ）と、イエローインクとシアンインクが吐出される領域（ｇ）と、シアンインクとマゼンタインクが吐出される領域（ｂ）と、ＣＭＹの３つのインクが吐出される領域（ｋ）をそれぞれ意味している。また、ｆc，ｆm，ｆyは、ＣＭＹ各インクを１種類のみ吐出したときにそのインクで覆われる面積の割合（「インク被覆率(Ink area coverage)」と呼ぶ）である。
【００６１】
インク被覆率ｆc，ｆm，ｆyは、図２０（Ｂ）に示すマーレイ・デービスモデルで与えられる。マーレイ・デービスモデルでは、例えばシアンインクのインク被覆率ｆcは、シアンのインク量ｄcの非線形関数であり、例えば１次元ルックアップテーブルによってインク量ｄcをインク被覆率ｆcに換算することができる。インク被覆率ｆc，ｆm，ｆyがインク量ｄc，ｄm，ｄyの非線形関数となる理由は、単位面積に少量のインクが吐出された場合にはインクが十分に広がるが、多量のインクが吐出された場合にはインクが重なり合うためにインクで覆われる面積があまり増加しないためである。他の種類のＭＹインクについても同様である。
【００６２】
分光反射率に関するユール・ニールセンモデルを適用すると、前記（４）式は以下の（５ａ）式または（５ｂ）式に書き換えられる。
【数５】
0004924414.tif 000006
ここで、ｎは１以上の所定の係数であり、例えばｎ＝１０に設定することができる。前記の（５ａ）式および（５ｂ）式は、ユール・ニールセン分光ノイゲバウアモデル(Yule-Nielsen Spectral Neugebauer Model)を表す式である。
【００６３】
本実施形態で採用するセル分割ユール・ニールセン分光ノイゲバウアモデル(Cellular Yule-Nielsen Spectral Neugebauer Model)は、上述したユール・ニールセン分光ノイゲバウアモデルのインク量空間を複数のセルに分割したものである。
【００６４】
図２１（Ａ）は、セル分割ユール・ニールセン分光ノイゲバウアモデルにおけるセル分割の例を示している。ここでは、説明の簡略化のために、ＣＭインクのインク量ｄc，ｄmの２つの軸を含む２次元インク量空間でのセル分割を描いている。なお、インク被覆率ｆc，ｆmは上述したマーレイ・デービスモデルにてインク量ｄc，ｄmと一意の関係にあるため、インク被覆率ｆc，ｆmを示す軸と考えることもできる。白丸は、セル分割のグリッド点（「格子点」と呼ぶ）であり、２次元のインク量（被覆率）空間が９つのセルＣ１〜Ｃ９に分割されている。各格子点に対応するインク量セット（ｄc，ｄm）は、分光反射率データベースＲＤＢに規定された格子点に対応するインク量セットとされている。すなわち、上述した分光反射率データベースＲＤＢを参照することにより、各格子点の分光反射率Ｒ（λ）を得ることができる。従って、各格子点の分光反射率Ｒ（λ）00，Ｒ（λ）10，Ｒ（λ）20・・・Ｒ（λ）33は、分光反射率データベースＲＤＢから取得することができる。
【００６５】
実際には、本実施形態ではセル分割もＣＭＹＫｌｃｌｍの６次元インク量空間で行うとともに、各格子点の座標も６次元のインク量セットψ（ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlm）によって表される。そして、各格子点のインク量セットψ（ｄc，ｄm，ｄy，ｄk，ｄlc，ｄlm）に対応する格子点の分光反射率Ｒ（λ）が分光反射率データベースＲＤＢ（例えば光沢紙のもの）から取得されることとなる。
【００６６】
図２１（Ｂ）は、セル分割モデルにて使用するインク被覆率ｆcとインク量ｄcとの関係を示している。ここでは、１種類のインクのインク量の範囲０〜ｄcmaxも３つの区間に分割されており、各区間毎に０から１まで単調に増加する非線形の曲線によってセル分割モデルにて使用する仮想的なインク被覆率ｆcが求められる。他のインクについても同様にインク被覆率ｆm，ｆyが求められる。
【００６７】
図２１（Ｃ）は、図２１（Ａ）の中央のセルＣ５内にある任意のインク量セット（ｄc，ｄm）にて印刷を行った場合の予測分光反射率Ｒs（λ）の算出方法を示している。インク量セット（ｄc，ｄm）にて印刷を行った場合の分光反射率Ｒ（λ）は、以下の（６）式で与えられる。
【数６】
0004924414.tif 000007
ここで、（６）式におけるインク被覆率ｆc，ｆmは図２１（Ｂ）のグラフで与えられる値である。また、セルＣ５を囲む４つの格子点に対応する分光反射率Ｒ（λ）11，（λ）12，（λ）21，（λ）22は分光反射率データベースＲＤＢを参照することにより取得することができる。これにより、（６）式の右辺を構成するすべての値を確定することができ、その計算結果として任意のインク量セットψ（ｄc，ｄm）にて印刷を行った場合の予測分光反射率Ｒs（λ）を算出することができる。波長λを可視波長域にて順次シフトさせていくことにより、可視波長域における予測分光反射率Ｒs（λ）を得ることができる。インク量空間を複数のセルに分割すれば、分割しない場合に比べて予測分光反射率Ｒs（λ）をより精度良く算出することができる。以上のようにして、ＲＰＭ Ｐ３ａ２がＩＣＭ Ｐ３ａ１の要請に応じて予測分光反射率Ｒs（λ）を予測することができる。
【００６８】
５．変形例
５−１：変形例１
図２２は、変形例においてＥＣＭ Ｐ３ａ３が設定する重み関数ｗ（λ）を模式的に示している。同図において、ターゲットＴＧから得られたターゲット分光反射率Ｒt（λ）が示されており、当該ターゲット分光反射率Ｒt（λ）と各等色関数ｘ（λ），ｙ（λ），ｚ（λ）との相関係数ｃx，ｃy，ｃzがＥＣＭ Ｐ３ａ３によって算出されている。そして、下記の（７）式によって本変形にかかる重み関数ｗ（λ）を算出する。
【数７】
0004924414.tif 000008
前記の（７）式においては、ターゲットＴＧから得られたターゲット分光反射率Ｒt（λ）との相関が高い等色関数ｘ（λ），ｙ（λ），ｚ（λ）ほど線形結合の際の重みが大きくなるようにされている。以上のようにして得られた重み関数ｗ（λ）においては、ターゲットＴＧのターゲット分光反射率Ｒt（λ）が大きい波長域についての重みを強調することができる。従って、各光源下での反射光の分光エネルギーのスペクトルが強くなりがちな波長域を重視した評価値Ｅ（φ）を得ることができる。すなわち、特にターゲットＴＧのターゲット分光反射率Ｒt（λ）が大きい波長域については、ターゲットＴＧのターゲット分光反射率Ｒt（λ）と予測分光反射率Ｒs（λ）とのずれを許容しないようなインク量セットφの最適解を得ることができる。むろん、重み関数ｗ（λ）は各等色関数ｘ（λ），ｙ（λ），ｚ（λ）に由来しているため、人間の知覚に適合した評価値Ｅ（φ）を得ることができる。
【００６９】
５−２：変形例２
図２３は、別の変形例においてＥＣＭ Ｐ３ａ３が設定する重み関数ｗ（λ）を模式的に示している。同図において、ターゲットＴＧから得られたターゲット分光反射率Ｒt（λ）をそのまま重み関数ｗ（λ）として適用している。このようにすることによっても、特にターゲットＴＧのターゲット分光反射率Ｒt（λ）が大きい波長域についてターゲットＴＧの分光反射率Ｒ（λ）とターゲット分光反射率Ｒt（λ）とのずれを許容しないようなインク量セットφの最適解を得ることができる。
【００７０】
５−３：変形例３
図２４は、別の変形例においてＥＣＭ Ｐ３ａ３が設定する重み関数ｗ（λ）を模式的に示している。同図において、５種類（標準昼光系のＤ５０光源，Ｄ５５光源，Ｄ６５光源、白熱電球系のＡ光源、蛍光ランプ系のＦ１１光源）の各光源の分光エネルギーＰ D50（λ），Ｐ D55（λ），Ｐ D65（λ），ＰA（λ），Ｐ F11（λ）が示されている。本変形例においては、下記の（８）式によって、これらの分光エネルギーＰ D50（λ），Ｐ D55（λ），Ｐ D65（λ），ＰA（λ），Ｐ F11（λ）を線形結合することにより重み関数ｗ（λ）を算出する。
【数８】
0004924414.tif 000009
【００７１】
前記の（８）式において、ｗ1〜ｗ5は各光源についての重みを設定する重み係数である。このように、光源の分光エネルギー分布Ｐ D50（λ），Ｐ D55（λ），Ｐ D65（λ），ＰA（λ），Ｐ F11（λ）に由来する重み関数ｗ（λ）を設定することにより、各光源下での反射光の分光エネルギーのスペクトルが強くなりがちな波長域を重視した評価値Ｅ（φ）を得ることができる。また、重み係数ｗ1〜ｗ5を調整することも可能である。例えば全光源における色の再現性をバランスよく確保したい場合にはｗ1＝ｗ2＝ｗ3＝ｗ4＝ｗ5とすればよいし、人工光源における色の再現性を重視したい場合にはｗ1，ｗ2，ｗ3＜ｗ4，ｗ5とすればよい。
【００７２】
５−４：変形例４
図２５は、変形例においてディスプレイ４０に表示されるＵＩ画面を示している。同図において、ＵＩ画面において複数のターゲット分光反射率Ｒt（λ）のグラフが表示されている。このようなＵＩ画面を表示させることにより、ユーザーがステップＳ１４０においてターゲットＴＧのターゲット分光反射率Ｒt（λ）を測定する代わりに、所望の波形のグラフをターゲットＴＧのターゲット分光反射率Ｒt（λ）として選択することができる。このようにすることにより、実際に分光反射率測定をしなくてもターゲット分光反射率Ｒt（λ）を設定することができる。むろん、グラフの波形をユーザーが直接編集できるようにしてもよい。例えば、新規の物体表面の開発を行う際に目標とするターゲット分光反射率Ｒt（λ）に編集しておけば、実際に物体表面を試作することなく、目標とするターゲット分光反射率Ｒt（λ）を有する見本チャートＳＣをプリンタ２０によって印刷させることができる。
【００７３】
５−５：変形例５
図２６は、変形例にかかる評価値（φ）を模式的に説明している。同図において、ターゲットＴＧのターゲット分光反射率Ｒt（λ）に対して前記５種類の光源を照射したときの色彩値（ターゲット色彩値）を上述した（１）式，図５によって算出する。一方、ＲＰＭ Ｐ３ａ２が予測した予測分光反射率Ｒs（λ）に対して前記５種類の光源を照射したときの色彩値（予測色彩値）も上述した（１）式（Ｒt（λ）をＲs（λ）に置き換えて使用），図５によって算出する。そして、各光源におけるターゲット色彩値と予測色彩値の色差ΔＥ（ΔＥ2000）をＣＩＥ  ＤＥ２０００の色差式に基づいて算出する。そして、各光源についての色差ΔＥをΔＥ D50，ΔＥ D55，ΔＥ D65，ΔＥA，ΔＥ F11とし、下記の（９）式によって評価値Ｅ（φ）を算出する。
【数９】
0004924414.tif 000010
前記の（２）式において、ｗ1〜ｗ5は各光源の重みを設定する重み係数であり、上述した変形例３の重み係数ｗ1〜ｗ5とほぼ同様の性質を有する。ここでも全光源における色の再現性をバランスよく確保したい場合にはｗ1＝ｗ2＝ｗ3＝ｗ4＝ｗ5とすればよいし、人工光源における色の再現性を重視したい場合にはｗ1，ｗ2，ｗ3＜ｗ4，ｗ5とすればよい。
【図面の簡単な説明】
【００７４】
【図１】印刷制御装置のハードウェア構成を示すブロック図である。
【図２】印刷制御装置のソフトウェア構成を示すブロック図である。
【図３】印刷データ生成処理の流れを示すフローチャートである。
【図４】ＵＩ画面の一例を示す図である。
【図５】分光反射率に基づいて色彩値を算出するための計算を説明する図である。
【図６】印刷データを示す図である。
【図７】インデックステーブルを示す図である。
【図８】印刷制御処理の全体的な流れを示すフローチャートである。
【図９】１Ｄ−ＬＵＴ生成処理の流れを示すフローチャートである。
【図１０】シアンの等色インクセットの分光反射率分布のグラフである。
【図１１】マゼンタの等色インクセットの分光反射率分布のグラフである。
【図１２】イエローの等色インクセットの分光反射率分布のグラフである。
【図１３】インク量セットを最適化する処理の流れを示す模式図である。
【図１４】インク量セットが最適化されていく様子を示す模式図である。
【図１５】１Ｄ−ＬＵＴを示す図である。
【図１６】印刷制御データ生成処理の流れを示すフローチャートである。
【図１７】３Ｄ−ＬＵＴを示す図である。
【図１８】プリンタの印刷方式を示す模式図である。
【図１９】分光反射率データベースを示す図である。
【図２０】分光ノイゲバウアモデルを示す図である。
【図２１】セル分割ユール・ニールセン分光ノイゲバウアモデルを示す図である。
【図２２】変形例にかかる重み関数を示す模式図である。
【図２３】変形例にかかる重み関数を示す模式図である。
【図２４】変形例にかかる重み関数を示す模式図である。
【図２５】変形例にかかるＵＩ画面を示す図である。
【図２６】変形例にかかる評価値を示す模式図である。
【符号の説明】
【００７５】
１０…コンピュータ、１１…ＣＰＵ、１２…ＲＡＭ、１３…ＲＯＭ、１４…ＨＤＤ、１５…ＧＩＦ、１６…ＶＩＦ、１７…ＩＩＦ、１８…バス、Ｐ１…ＯＳ、Ｐ１ａ…ＧＤＩ、Ｐ１ｂ…スプーラ、Ｐ２…ＡＰＬ、Ｐ２ａ…ＵＩＭ、Ｐ２ｂ…ＭＣＭ、Ｐ２ｃ…ＰＤＧ、Ｐ３…ＰＤＶ、Ｐ３ａ１…ＩＣＭ、Ｐ３ａ２…ＲＰＭ、Ｐ３ａ３…ＥＣＭ、Ｐ３ａ４…ＬＯＭ、Ｐ４…ＭＤＶ、Ｐ５…ＤＤＶ。
【図１】
0004924414.tif 000011
【図２】
0004924414.tif 000012
【図３】
0004924414.tif 000013
【図４】
0004924414.tif 000014
【図５】
0004924414.tif 000015
【図６】
0004924414.tif 000016
【図７】
0004924414.tif 000017
【図８】
0004924414.tif 000018
【図９】
0004924414.tif 000019
【図１０】
0004924414.tif 000020
【図１１】
0004924414.tif 000021
【図１２】
0004924414.tif 000022
【図１３】
0004924414.tif 000023
【図１４】
0004924414.tif 000024
【図１５】
0004924414.tif 000025
【図１６】
0004924414.tif 000026
【図１７】
0004924414.tif 000027
【図１８】
0004924414.tif 000028
【図１９】
0004924414.tif 000029
【図２０】
0004924414.tif 000030
【図２１】
0004924414.tif 000031
【図２２】
0004924414.tif 000032
【図２３】
0004924414.tif 000033
【図２４】
0004924414.tif 000034
【図２５】
0004924414.tif 000035
【図２６】
0004924414.tif 000036