(57)【要約】      （修正有）【課題】画像の再現色を効率的に算出する。【解決手段】プロファイル生成部1103は、出力デバイスが媒体上に形成したプライマリ色の色票に可視光域の光とより波長が短い可視光域外の光を入射した場合の可視光域の分光反射率の測定結果から、出力デバイスが媒体上に形成する、入力画像信号の色空間の格子点に対応する色票それぞれに可視光域と可視光域外の光を入射した場合の可視光域の分光反射率を演算し、格子点に対応する色票それぞれの分光反射率を記述した分光反射率プロファイルを生成して、格子点の画像信号値と分光反射率プロファイルの関係をプロファイルテーブルに記録する。観察光データ入力部1105は、画像を観察する観察光の分光輝度値を取得する。分光演算部1106は、分光輝度値、プロファイルテーブル、および、分光反射率プロファイルを使用して、出力デバイスが形成する画像を観察光の下で観察する場合の画像の再現色を算出する。【選択図】図４
2011155519.tif 000002
【特許請求の範囲】
【請求項１】
出力デバイスが媒体上に形成したプライマリ色の色票に可視光域の光、および、前記可視光域より波長が短い可視光域外の光を入射した場合の前記可視光域の分光反射率の測定結果から、前記出力デバイスが前記媒体上に形成する、入力画像信号の色空間の格子点に対応する色票それぞれに前記可視光域および前記可視光域外の光を入射した場合の前記可視光域の分光反射率を演算する演算手段と、
前記格子点に対応する色票それぞれの前記分光反射率を記述した分光反射率プロファイルを生成して、前記格子点の画像信号値と前記分光反射率プロファイルの関係をプロファイルテーブルに記録する生成手段と、
画像を観察する観察光の分光輝度値を取得する取得手段と、
前記分光輝度値、前記プロファイルテーブル、および、前記分光反射率プロファイルを使用して、前記出力デバイスが形成する画像を前記観察光の下で観察する場合の前記画像の再現色を算出する算出手段とを有することを特徴とする色処理装置。
【請求項２】
前記演算手段は、前記出力デバイスの色分解テーブルを参照して注目格子点の画像信号値に対応する色材量値を取得し、前記色材量値から前記プライマリ色が前記注目格子点に対応する色票上で占める面積の比率を示す面積率を取得し、前記面積率と前記分光反射率プロファイルから前記注目格子点に対応する色票の分光反射率を演算することを特徴とする請求項1に記載された色処理装置。
【請求項３】
前記演算手段は、前記プライマリ色ごとに演算した前記分光反射率を波長ごとに積算して前記注目格子点に対応する色票の分光反射率にすることを特徴とする請求項2に記載された色処理装置。
【請求項４】
出力デバイスが媒体上に形成したプライマリ色の色票に可視光域の光、および、前記可視光域より波長が短い可視光域外の光を入射した場合の前記可視光域の分光反射率の測定結果から、前記出力デバイスが前記媒体上に形成する、入力画像信号の色空間の格子点に対応する色票それぞれに前記可視光域および前記可視光域外の光を入射した場合の前記可視光域の分光反射率を演算し、
前記格子点に対応する色票それぞれの前記分光反射率を記述した分光反射率プロファイルを生成して、前記格子点の画像信号値と前記分光反射率プロファイルの関係をプロファイルテーブルに記録し、
画像を観察する観察光の分光輝度値を取得し、
前記分光輝度値、前記プロファイルテーブル、および、前記分光反射率プロファイルを使用して、前記出力デバイスが形成する画像を前記観察光の下で観察する場合の前記画像の再現色を算出することを特徴とする色処理方法。
【請求項５】
コンピュータ装置を制御して、請求項1から請求項3の何れか一項に記載された色処理装置の各手段として機能させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、蛍光材を含む媒体や色材を使用して形成される画像の再現色を算出する色処理に関する。
【背景技術】
【０００２】
画像形成に使用する媒体および色材には、媒体の白色度や色材の彩度の向上（以下、発色性の向上）のために、蛍光材が付与される場合がある。蛍光材は、紫外光を代表とする可視光域より波長が短い領域（以下、可視光域外）の光によって、可視光域内の波長をもつ光を励起する特性がある。蛍光材を含む媒体や色材によって形成された画像を観察すると、観察光（照明光）が含む可視光域外の光の波長の違いにより色再現に違いが生じる。従って、蛍光材を含む媒体や色材を使用して形成された出力画像の色再現を推定するには、観察光により励起された蛍光材の分光反射率と、色材の分光反射率の両者を考慮する必要がある。
【０００３】
特許文献1は、媒体の蛍光成分（増白剤）の影響を考慮して、出力画像の色再現を推定する方法を開示する。つまり、タングステンランプのような基準光源の下で媒体の分光反射率Pref(λ)を測定する。また、蛍光灯のような観察光源の下で媒体の分光反射率Pflu(λ)を測定し、差分Pflu(λ)-Pref(λ)を媒体の蛍光成分に起因する差成分ΔP(λ)にする。次に、基準光源の下で色材によって形成した試料の分光反射率Rref(λ)を測定する。そして、基準光源の下の試料の分光反射率Rref(λ)と差成分ΔP(λ)から、観察光源の下の試料の分光反射率を推定することで、試料の再現色を推定する。つまり、任意の観察光の下における試料の再現色を推定するには、媒体について二回の測定と、各色の試料について一回の測定が必要である。
【０００４】
特許文献2は、媒体および色材の蛍光成分を考慮して、出力画像の色再現を推定する方法を開示する。つまり、媒体の種類と754組の色材量値（CMYK値）の組み合わせそれぞれについて、300〜780nmの入射光波長に対する380〜780nmの反射光波長の反射率を測定し記録した二次元励起テーブルを用意する。そして、二次元励起テーブルに記述された分光反射率と、観察光源の完全拡散面の三刺激値XnYnZnから表色値（Lab値）を計算し、色材量値に対する再現色を推定する。なお、任意の色材量値に対する再現色は、複数の色材量値と表色値の対応関係に基づき線形補間により推定する。
【０００５】
特許文献2の二次元励起テーブルを作成するには、可視光域外の300nmから可視光域の780nmまで10nm間隔の入射光について測定するとして、色材量値一組当り49回の測定が必要になる。その上、特許文献2の技術は、二次元励起テーブルに記述された分光反射率から計算した観察光の下の表色値と色材量値の関係から任意の色材量値に対する再現色を推定する。色材量値と表色値は非線形の関係にあることが多い。非線形な関係から再現色を高精度に推定するには、二次元励起テーブルには多数のエントリが必要で、特許文献2は754組の色材量値をエントリとする二次元励起テーブルを記載する。つまり、特許文献2のように出力色を754色にする場合、媒体一種類当り49×754=36,652回の測定を必要とし、複数の媒体に対応するための測定は膨大な数になる。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特許3862997号公報
【特許文献２】特開2003-110867号公報
【非特許文献】
【０００７】
【非特許文献１】三宅洋一「分光画像処理入門」（財）東京大学出版会、2006年2月、ISBN 9784130621410
【発明の概要】
【発明が解決しようとする課題】
【０００８】
本発明は、蛍光材を含む媒体や色材を使用して形成される画像の再現色を効率的に算出することを目的とする。
【課題を解決するための手段】
【０００９】
本発明は、前記の目的を達成する一手段として、以下の構成を備える。
【００１０】
本発明にかかる色処理は、出力デバイスが媒体上に形成したプライマリ色の色票に可視光域の光、および、前記可視光域より波長が短い可視光域外の光を入射した場合の前記可視光域の分光反射率の測定結果から、前記出力デバイスが前記媒体上に形成する、入力画像信号の色空間の格子点に対応する色票それぞれに前記可視光域および前記可視光域外の光を入射した場合の前記可視光域の分光反射率を演算し、前記格子点に対応する色票それぞれの前記分光反射率を記述した分光反射率プロファイルを生成して、前記格子点の画像信号値と前記分光反射率プロファイルの関係をプロファイルテーブルに記録し、画像を観察する観察光の分光輝度値を取得し、前記分光輝度値、前記プロファイルテーブル、および、前記分光反射率プロファイルを使用して、前記出力デバイスが形成する画像を前記観察光の下で観察する場合の前記画像の再現色を算出することを特徴とする。
【発明の効果】
【００１１】
本発明によれば、蛍光材を含む媒体や色材を使用して形成される画像の再現色を効率的に算出することができる。
【図面の簡単な説明】
【００１２】
【図１】実施例1の色処理装置の構成例を説明するブロック図。
【図２】プロファイルテーブルの一例を説明する図。
【図３】分光反射率プロファイルの一例を説明する図。
【図４】任意の観察光下における形成色の色再現を推定する機能構成例を説明するブロック図。
【図５】任意の観察光下における形成色の色再現を推定する処理を説明するフローチャート。
【図６】任意の観察光下における入力画像の色再現を推定する処理を説明するフローチャート。
【図７】入力信号値を色材量値に色分解する色分解テーブルの一例を示す図。
【図８】各色材の色材量値に対する有効面積率を示すテーブルを説明する図。
【図９】分光ノイゲバウア方程式による分光反射率プロファイルを生成する機能構成例を説明するブロック図。
【図１０】分光ノイゲバウア方程式による分光反射率プロファイルを生成する処理を説明するフローチャート。
【図１１】プロファイルテーブルの一例を説明する図。
【図１２】プロファイル生成部による分光反射プロファイルPiの生成の詳細を説明するフローチャート。
【図１３】分光反射率の算出の詳細を説明するフローチャート。
【発明を実施するための形態】
【００１３】
以下、本発明にかかる実施例の色処理を図面を参照して詳細に説明する。
【実施例１】
【００１４】
［装置の構成］
図1のブロック図により実施例1の色処理装置の構成例を説明する。
【００１５】
CPU101は、RAM103をワークメモリとして、ROM102やハードディスクドライブ(HDD)105に格納されたオペレーティングシステム(OS)や各種プログラムを実行し、システムバス109を介して、後述する構成を制御する。各種プログラムには、後述する色処理（画像処理）を実現するアプリケーションプログラム（以下、色処理AP）が含まれる。また、ROM102やHDD105には、後述する色処理に利用される各種データが格納されている。
【００１６】
CPU101によって実行される色処理APは、モニタ106にユーザインタフェイス(UI)を表示する。ユーザは、操作部104のキーボードやマウスを操作して、UIに対して指示やデータを入力する。色処理APは、ユーザ指示に従い、汎用インタフェイス(I/F)108を介して、プリンタ111にカラーチャートを印刷させる、測定器112にカラーチャートの各パッチの分光反射率を測定させる、などの処理を実行する。また、取得した分光反射率の測定結果に基づき、後述する分光反射率の推定など、様々な色処理を行う。
【００１７】
［分光反射率プロファイル］
HDD105には、色処理APが利用する、出力デバイス（本実施例ではプリンタ111）が入力する画像信号に従い形成する色（形成色）の分光反射率を推定するためのプロファイルテーブルが格納されている。プロファイルテーブルは、出力デバイスの入力画像信号の値を格子点の入力値とし、各格子点に対応する形成色の分光反射率プロファイルの番号（ポインタ）を記述した多次元ルックアップテーブル(LUT)である。例えば出力デバイスがRGB色空間の画像信号値を入力する場合、格子点はRGB値の離散的な組み合わせに対応する。
【００１８】
図2によりプロファイルテーブルの一例を説明する。図2の例は、各格子点に対して離散的なRGB値の組み合わせを入力値として、n番目の格子点に対応するプロファイルの番号Pn（またはアドレスポインタ）を記録したLUTの例を示している。このLUTは、9×9×9の格子点を有するから、0〜728番の各格子点に対してP0〜P728番のプロファイルが定義されている。なお、P0番の分光反射率プロファイルは、媒体上に画像を形成しない媒体自体の分光反射率を示す。
【００１９】
図3により分光反射率プロファイルの一例を説明する。分光反射率プロファイルは、ある形成色に対して300〜780nmの光を入射した場合の380〜780nmの光の反射率を10nmごとに記述した二次元テーブルである。つまり、分光反射率プロファイルを参照すれば、ある形成色について、色材が含む蛍光成分により可視光域外の波長λinの光が入射された場合に励起される可視光域の波長λoの光の反射率[%]を知ることができる。また、可視光域の波長λinの光が入射した場合の可視光域の波長λoの光の反射率[%]を知ることができる。
【００２０】
離散的な信号値の間の信号値に対する分光反射率プロファイルを得るには、当該信号値を囲む格子点の信号値に対応する分光反射率プロファイルを取得する。そして、それら分光反射率プロファイルの各波長ごとの反射率から当該信号値の各波長ごとの反射率を線形補間して、分光反射率プロファイルを生成する。この線形補間により、離散的な信号値の間の信号値に対する分光反射率プロファイルの精度を向上するには、プロファイルテーブルの格子点数を充分な数にする必要がある。図2には9×9×9=729の格子点を有するプロファイルテーブルを示したが、例えば17×17×17=4913の格子点をもつプロファイルテーブルを利用すればより高精度な推定が可能になる。
【００２１】
［観察光下の色再現］
画像は任意の照明下で観察される。従って、形成色の再現色を推定するには、形成色の分光反射率を推定するとともに、観察光の下における色再現を推定する必要がある。
【００２２】
図4のブロック図により任意の観察光下における形成色の色再現を推定する機能構成例を説明する。この機能構成は、色処理APによって実現される。また、図5のフローチャートにより任意の観察光下における形成色の色再現を推定する処理を説明する。
【００２３】
信号値入力部1102は、例えばユーザが操作部104を用いて指示したRGB値を入力する(S101)。プロファイル生成部1103は、例えばHDD105の所定領域に格納されたプロファイルテーブルと分光反射率プロファイルを参照して、入力RGB値に対応する形成色の分光反射率プロファイルを生成する(S102)。形成色の分光反射率プロファイルは、入力信号値および形成色に対応付けて例えばHDD105の所定領域に格納される。
【００２４】
次に、観察光データ入力部1105は、測定器112（例えば分光光度計）から形成色の観察光の分光輝度値の測定結果を入力し(S103)、測定結果を観察光データとして観察光の名称や記号に対応付けて例えばHDD105の所定領域に格納する。分光演算部1106は、観察光データと、入力RGB値に対応する形成色の分光反射率プロファイルを積算する(S104)。そして、演算結果の分光輝度プロファイルを入力信号値、形成色および観察光源の名称や記号に関連付けて例えばHDD105の所定領域に格納する(S105)。
【００２５】
分光反射率プロファイルは、入射光と反射光の波長の組み合わせごとに反射率を示す二次元テーブルである。従って、入射光の波長ごとに観察光データが示す分光輝度と反射率を積算し、反射光の波長ごとに積を総和して、反射光の波長と輝度の関係を示す一次元のテーブルを形成色の分光輝度プロファイルとする。
【００２６】
［画像の色再現の推定］
図6のフローチャートにより任意の観察光下における入力画像の色再現を推定する処理を説明する。
【００２７】
色処理APは、例えばユーザが操作部104を用いて指示する、観察光源および画像データの情報を入力する(S201)。そして、指示された画像データをHDD105の所定領域に格納し(S202)、画像データからラスタ順に注目画素の例えばRGB値（以下、注目RGB値）を取り出す(S203)。そして、注目RGB値と観察光源に対応するCIEXYZ値（以下、対応XYZ値）がHDD105の所定領域に格納されているか否かを判定し(S204)、対応XYZ値が格納されている場合は処理をステップS210に進める。
【００２８】
対応するXYZ値が格納されていない場合、色処理APは、注目RGB値と観察光源に対応する分光輝度プロファイル（以下、対応分光輝度プロファイル）が例えばHDD105に存在するか否かを判定する(S205)。対応分光輝度プロファイルが存在する場合は処理をステップS209に進める。
【００２９】
また、対応分光輝度プロファイルが存在しない場合、色処理APは、注目RGB値に対応する分光反射率プロファイル（以下、対応分光反射率プロファイル）が例えばHDD105に存在するか否かを判定する(S206)。対応分光反射率プロファイルが存在する場合は、図5に示すステップS103からS105の処理によって対応分光輝度プロファイルを生成し(S207)、処理をステップS209に進める。
【００３０】
また、対応分光反射率プロファイルが存在しない場合、色処理APは、図5に示すステップS102からS105の処理によって対応分光輝度プロファイルを生成し(S208)、処理をステップS209に進める。
【００３１】
色処理APは、対応分光輝度プロファイルを取得または生成すると、次式によって観察光下の注目RGB値に対する形成色のCIEXYZ値（対応XYZ値）を計算する(S209)。
X = K∫S(λ)x(λ)dλ
Y = K∫S(λ)y(λ)dλ
Z = K∫S(λ)z(λ)dλ        …(1)
K = 100/∫S(λ)y(λ)dλ
ここで、積分範囲は可視光域(380〜780nm)、
S(λ)は分光輝度、
x(λ)y(λ)z(λ)は等色関数。
【００３２】
色処理APは、対応XYZ値を取得または計算すると、当該XYZ値を観察光源、RGB値に関連付けてHDD105の所定領域にラスタ順に格納する(S210)。そして、ステップS211の判定により、入力した画像データの全画素についてステップS203からS210の処理を行ったか否かを判定し、未了であれば処理をステップS203に戻し、次の画素を処理する。
【００３３】
入力した画像データの全画素についてステップS203からS210の処理が終了した場合、色処理APは、HDD105の所定領域にラスタ順に格納したXYZ値を読み出して、XYZ値が表す画像を例えばモニタ106に表示する(S212)。その際、モニタ106のデバイスプロファイルを参照して、XYZ値をRGB値に変換することは言うまでもない。例えばモニタ106に表示される画像は、ユーザが指示した画像データに基づきプリンタ111が形成する画像をユーザが指示する観察光源の下で観察したプレビュー画像に相当する。
【００３４】
このように、可視光域外の入射光に対する可視光域の反射率を考慮した分光反射率プロファイルの利用により、媒体や色材に含まれる蛍光成分と観察光の関係を考慮した高精度のプレビュー画像が得られる。
【００３５】
［測定の効率化」
上記では、一種類の媒体に対して、729組または4913組の入力信号値に対応する形成色に対して300〜780nmの10nmステップの入射光に対する380〜780nmの反射率を測定した分光反射率プロファイルを使用する例を説明した。つまり、一種類の媒体について、49×729=35,721回または49×4,913=240,737回の測定が必要になる。そこで、本実施例においては、格子点数より少ない分光反射率プロファイルから、分光ノイゲバウア方程式を利用して、各格子点の分光放射率プロファイルを算出する。なお、分光ノイゲバウア方程式の詳細は、非特許文献1に記載されている。
【００３６】
●分光ノイゲバウア方程式の概要
本実施例においては、四色プリンタを前提にするので、ノイゲバウアプライマリとして、例えば、次のプライマリ色の色票（以下、プライマリ色票）の分光反射率を測定したデータを用意する。
濃度0%（媒体の色）、
出力デバイスのべた印刷色：一次色C、M、Y、K（四色）、
二次色CM、MY、YK、MK、CK（五色）、
三次色CMY、CMK、MYK、CYK(四色）、
四次色CMYK（一色）。
【００３７】
上記15色のプライマリ色票の測定結果から任意の形成色の色を推定する。なお、単色の中間調として網点印刷による単色の中間調（25%、50%、75%）や、二次色以上の中間調を追加してもよい。分光ノイゲバウア方程式は次式で示される。
R(λ) = {ΣpFp・Rp(λ)1/n}n      …(2)
ここで、Rp(λ)はノイゲバウアプライマリpの分光反射率、
Fpは各ノイゲバウアプライマリpが色票上で占める面積の比率（部分面積率）、
nは色材における光の透過・散乱の非線形性を考慮するための修正パラメータ。
【００３８】
出力デバイスの色材数がLの場合、ノイゲバウアプライマリの総数Pは2L個になる。CMYK四色印刷におけるノイゲバウアプライマリの総数Pは16であり、W、K、Y、YK、M、MK、MY(=R)、MYK(=RK)、C、CK、CY(=G)、CYK(=GK)、CM(=B)、CMK(=BK)、CMY、CMYKに対応する。
【００３９】
部分面積率Fpは、単色の色材が印刷物（例えば色票）上で占める面積の比率を示す有効面積率から、デミチェルの方程式によって求められる。CMYK四色印刷における下地（媒体自体）、単色（単一色材）、混色（色材の混在）の部分面積率Fpを次式に示す。
FW = (1 - ac)(1 - am)(1 - ay)(1 - ak)
FK = (1 - ac)(1 - am)(1 - ay)ak
FY = (1 - ac)(1 - am)ay(1 - ak)
FYK = (1 - ac)(1 - am)ayak
FM = (1 - ac)am(1 - ay)(1 - ak)
FMK = (1 - ac)am(1 - ay)ak
FC = ac(1 - am)(1 - ay)(1 - ak)
FCKw = ac(1 - am)(1 - ay)ak
FR = (1 - ac)amay(1 - ak)
FRK = (1 - ac)amayak
FG = ac(1 - am)ay(1 - ak)
FGK = ac(1 - am)ayak
FB = acam(1 - ay)(1 - ak)
FBK = acam(1 - ay)ak
FCMY = acamay(1 - ak)
FCMYK = acamayak        …(3)
ここで、ac、am、ay、akは各色材の有効面積率（図8参照）。
【００４０】
図7により入力信号値を色材量値に色分解する色分解テーブルの一例を示す。色分解テーブルは、図2に示す分光反射率プロファイルテーブルと同様の離散的なRGB値（格子点）に対する、CMYK各色材の使用量に相当する色材量値を示す。なお、格子点の間のRGB値に対応する色材量値は、当該RGB値を囲む格子点の色材量値から補間演算する。
【００４１】
プリンタのような出力デバイスにおいて、二次色、三次色、四次色のべた印刷が重なることはほとんどない。従って、色分解テーブルを参照して、ノイゲバウアプライマリに対応する出力デバイスの入力信号値に対する色材量値を取得する。
【００４２】
図8により各色材の色材量値に対する有効面積率を示すテーブルを説明する。テーブル601〜604はそれぞれシアンC色材、マゼンタM色材、イエローY色材、ブラックK色材の有効面積率ac、am、ay、akを示す。なお、テーブル601〜604において色材量値は離散的な値を示す。従って、テーブルに存在しない色材量値に対応する有効面積率は、当該色材量値を挟むテーブルに存在する色材量値から補間演算する。また、各色材の有効面積率の取得方法として様々な手法が知られている。それらの手法を利用して有効面積率を取得してもよい。
【００４３】
●分光反射率プロファイルの生成
図9のブロック図により分光ノイゲバウア方程式による分光反射率プロファイルを生成する機能構成例を説明する。この機能構成は、色処理APによって実現される。また、図10のフローチャートにより分光ノイゲバウア方程式による分光反射率プロファイルを生成する処理を説明する。
【００４４】
分光反射率入力部1001は、例えば300nmから780nmを10nm間隔で、プライマリ色票それぞれの分光反射率を測定器112（例えば分光光度計）から入力する(S121)。プロファイル生成部1002は、入力された各プライマリ色票の分光反射率を図3に示す二次元テーブルの形式にしてHDD105の所定領域に格納する(S122)。
【００４５】
図11によりプロファイルテーブルの一例を説明する。プロファイルテーブルは、色票番号で示される各プライマリ色票の色材量値と、二次元テーブル（分光反射率プロファイル）の番号Niの対応関係を示す。つまり、プロファイル生成部1002は、分光反射率プロファイルNiを生成するとともに、プライマリ色票と二次元テーブルの関係を記述したプロファイルテーブルを生成する。なお、図11には色材量値の最大を255として示す。しかし、CMYKの各色材量値の合計によって媒体に供給される色材量が、媒体が許容可能な色材量（最大色材量）を超える場合、媒体に供給する色材量は最大色材量に制限する必要がある。例えば、色材量値255に対応する色材量が100%、出力デバイスの最大色材量が200%とすると、色票番号15の色材量値は(255×2/4=127, 127, 127, 127)になる。また、色票番号14の色材量値は(255×2/3=170, 0, 170, 170)になる。
【００４６】
面積率算出部1003は、HDD105の所定領域に格納された有効面積率テーブル（図8）を参照して、プロファイルテーブル（図2）の各入力信号値（各格子点）に対する有効面積率を算出する(S123)。プロファイル生成部1002は、プロファイルテーブル（図2）の各入力信号値に対応する分光反射率プロファイルPi（図3）を生成する(S124)。その際、プロファイルテーブル（図2）、色分解テーブル（図7）、面積率算出部1003が算出した有効面積率、プロファイルテーブル（図11）、および、分光反射率プロファイルNiを参照して、分光ノイゲバウア方程式により分光反射率を演算する。そして、生成した分光反射率プロファイルPiをHDD105の所定領域に格納する(S125)。
【００４７】
図12のフローチャートによりプロファイル生成部1002による分光反射プロファイルPiの生成の詳細を説明する。
【００４８】
プロファイル生成部1002は、カウンタjを0に初期化し(S801)、プロファイルテーブル（図2）の格子点数Gを取得し(S802)、プロファイルテーブル（図2）のj番目の格子点の入力信号値を取得する(S803)。そして、色分解テーブル（図7）、有効面積率、プロファイルテーブル（図11）、および、分光反射率プロファイルNiを参照し、分光ノイゲバウア方程式を用いてj番目の格子点の入力信号値に対する各波長の分光反射率を算出する(S804)。
【００４９】
プロファイル生成部1002は、算出した各波長の分光反射率を分光反射率プロファイルPjとしてHDD105の所定領域に格納し(S805)、カウンタjとG-1を比較する(S806)。j＜G-1の場合はカウンタjをインクリメントして(S807)、処理をステップS803に戻し、分光反射率プロファイルPjが未算出の格子点がなくなるまでステップS803からS805の処理を繰り返す。
【００５０】
図13のフローチャートにより分光反射率の算出(S804)の詳細を説明する。
【００５１】
プロファイル生成部1002は、出力デバイスの色材数Lを取得し(S901)、注目格子点の入力信号値を取得する(S902)。そして、色分解テーブルを参照して、入力信号値に対する色材量値を取得し(S903)、面積率算出部1003から、入力信号値に対する各色材の有効面積率を取得する(S904)。そして、カウンタpを0に初期化し(S905)、入力信号値に対する部分面積率Fp（例えば、式(3)に示すFWからFCMYKの何れか）を算出する(S906)。さらに、色票番号pのプライマリ色票pの分光反射率プロファイルNpを取得する(S907)。
【００５２】
次に、プロファイル生成部1002は、分光反射率プロファイルNpを使用して、λin=300nmから380nm未満に対する可視光域の分光反射率Rp(λ)に部分面積率Fpを積算する。つまり、分光ノイゲバウア方程式により、蛍光による分光反射率Rp(λ)を算出する(S908)。そして、λin=380nmから780nmに対する分光反射率Rp(λ)に部分面積率Fpを積算する。つまり、分光ノイゲバウア方程式により、可視光の分光反射率Rp(λ)を算出する(S909)。
【００５３】
次に、プロファイル生成部1002は、カウンタpと2Lを比べて(S910)、p＜2Lならば、カウンタpをインクリメントし(S911)、処理をステップS906へ戻す。つまり、入力信号値に対する部分面積率Fpそれぞれの蛍光による分光反射率と可視光の分光反射率を算出する。
【００５４】
p個のノイゲバウアプライマリについて上記の分光反射率の算出が終了すると、プロファイル生成部1002は、ステップS908、S909で算出した分光放射率を波長ごとに積算（式(2)参照）する。そして、積算結果から二次元テーブル（分光反射率プロファイルPi）を生成する(S912)。
【００５５】
離散的な入力信号値を格子点とするプロファイルテーブル（図2）の各格子点に対応する、蛍光を含む可視光域の分光反射率プロファイルPiを利用して、蛍光成分を含む媒体、色材を使用する場合の再現色を高精度に演算することが可能である。その際、分光反射率プロファイルPiは、ノイゲバウアプライマリに対応する少数のプライマリ色票に可視光域外から可視光域の光を入射した場合の分光反射率の測定結果から生成する。例えば、色材数がLの場合は、少なくとも2L個のプライマリ色票の分光反射率を測色結果から分光反射率プロファイルPiを生成することができ、分光反射率プロファイルPiを生成するための測定の効率化を図ることができる。
【００５６】
［変形例］
上記実施例においては、多次元ルックアップテーブルと補間処理を用いて形成色の分光反射率（再現色）を算出する例を説明した。しかし、線形補間を用いると、色を推定する精度が劣化する場合がある。より精度が高い色（分光反射率）を推定するには、ノイゲバウアプライマリの分光反射率プロファイルNiを用いて、分光ノイゲバウア方程式により、形成色の分光反射率（再現色）を演算すればよい。
【００５７】
［その他の実施例］
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステムあるいは装置のコンピュータ（又はCPUやMPU等）がプログラムを読み出して実行する処理である。
【図１】
2011155519.tif 000003
【図２】
2011155519.tif 000004
【図３】
2011155519.tif 000005
【図４】
2011155519.tif 000006
【図５】
2011155519.tif 000007
【図６】
2011155519.tif 000008
【図７】
2011155519.tif 000009
【図８】
2011155519.tif 000010
【図９】
2011155519.tif 000011
【図１０】
2011155519.tif 000012
【図１１】
2011155519.tif 000013
【図１２】
2011155519.tif 000014
【図１３】
2011155519.tif 000015