(57)【要約】【課題】面積変調階調と、濃度変調階調との双方の色表現を有するグラビア印刷等における印刷物の色予測を精度よく行うことができる色予測装置、色予測方法およびプログラムを提供する。【解決手段】印刷機が印刷する再現色を予測する色予測装置であって、指令網点面積率に対応する濃度階調領域及び濃度階調領域の出現率とを求める濃度階調出現率算出部と、網点による再現色を示す分光反射率を、出現する濃度階調領域及び濃度階調領域の出現率に基づいて算出する分光データ算出部と、デバイス値と測色値の関係が記述されている色特性記述ファイルを読み込み、測色値を計算する測色値算出部とを備える。【選択図】図１
2015161523.tif 000002
【特許請求の範囲】
【請求項１】
指令網点面積率と、前記指令網点面積率でインキを印刷することにより形成される網点が含む濃度階調領域と、前記濃度階調領域の各々の出現率との関係に基づき、供給される前記指令網点面積率に対応する濃度階調領域及び前記濃度階調領域の出現率とを求める濃度階調出現率算出部と、
濃度階調領域及び前記濃度階調領域の各々の出現率と、分光反射率との対応関係を示す第１の計算モデルにより、前記指令網点面積率により形成された網点による再現色を示す分光反射率を、出現する前記濃度階調領域及び前記濃度階調領域の前記出現率に基づいて算出する分光データ算出部と、
デバイス値と測色値の関係が記述されている色特性記述ファイルを読み込み、デバイス値に対する測色値を計算する測色値算出部とを備える、色予測装置。
【請求項２】
前記分光データ算出部が、
前記濃度階調領域及び前記濃度階調領域の各々の出現率と、前記濃度階調領域の分光反射率から求めた分光光学濃度との対応関係を示す第２の計算モデルにより、前記指令網点面積率により形成された網点による再現色を示す分光光学濃度を算出する光学濃度算出部と、
前記分光光学濃度から、前記再現色の前記分光反射率を算出する分光反射率算出部とを有する、請求項１に記載の色予測装置。
【請求項３】
前記第１の計算モデルが、指令網点面積率によるインキが印刷された印刷部分の網点が含む前記濃度階調領域毎に、前記濃度階調領域の分光反射率と前記濃度階調領域の前記出現率とを乗算した結果を加算して、前記指令網点面積率の印刷部分の分光反射率を算出する式であり、
前記指令網点面積率と、前記濃度階調領域の出現率との対応関係を示す対応関数が予め書き込まれて記憶されている出現率テーブルデータベースをさらに、含み、
前記分光データ算出部が、前記出現率テーブルデータベースから、前記指令網点面積率に対応する前記濃度階調領域の分光反射率と、前記濃度階調領域の出現率を読み出し、前記第１の計算モデルにより、前記指令網点面積率に対応する前記印刷部分の分光反射率を求める、請求項１に記載の色予測装置。
【請求項４】
前記第２の計算モデルが、指令網点面積率によるインキが印刷された印刷部分の網点が含む前記濃度階調領域毎に、前記濃度階調領域の分光光学濃度と前記濃度階調の前記出現率とを乗算した結果を加算して、前記指令網点面積率の印刷部分の分光光学濃度を算出する式であり、
前記指令網点面積率と、前記濃度階調領域の出現率との対応関係を示す対応関数が予め書き込まれて記憶されている出現率テーブルデータベースをさらに、含み、
前記分光データ算出部が、前記出現率テーブルデータベースから、前記指令網点面積率に対応する前記濃度階調領域の分光光学濃度と、前記濃度階調領域の出現率を読み出し、前記第２の計算モデルにより、前記指令網点面積率に対応する前記印刷部分の分光光学濃度を求める、請求項２に記載の色予測装置。
【請求項５】
前記インキをベタで印刷した印刷部分を測定した分光反射率から吸収係数及び散乱係数を求め、求めた前記吸収係数及び前記散乱係数により、前記濃度階調領域のインキの厚さに対応する前記濃度階調領域毎の濃度階調分光反射率をクベルカ・ムンクの式により算出し、前記印刷部分の実測の分光反射率と、前記濃度階調領域の分光反射率を用いて前記第１の計算モデルにより算出する分光反射率との差が最小となるように、前記第１の計算モデルにおける前記出現率をフィッティングし、得られた前記濃度階調領域毎の出現率を用い、前記濃度階調毎の指令網点面積率と、前記濃度階調領域の出現率との対応関係を示す前記対応関数が求められている、請求項３に記載の色予測装置。
【請求項６】
前記インキをベタで印刷した印刷部分を測定した分光反射率から吸収係数及び散乱係数を求め、求めた前記吸収係数及び前記散乱係数により、前記濃度階調領域のインキの厚さに対応する前記濃度階調領域毎の濃度階調分光光学濃度をクベルカ・ムンクの式により算出し、前記印刷部分の実測の分光反射率から求めた分光光学濃度と、前記濃度階調分光光学濃度を用いて前記第２の計算モデルにより算出する分光光学濃度との差が最小となるように、前記計算モデルにおける前記出現率をフィッティングし、得られた前記濃度階調領域毎の出現率を用い、前記濃度階調毎の指令網点面積率と、前記濃度階調領域の出現率との対応関係を示す前記対応関数が求められている、請求項４に記載の色予測装置。
【請求項７】
前記インキが複数の原色のインキを複数混合した特色インキの場合、特色インキに含まれる原色インキの比率に応じて吸収係数及び散乱係数の各々を加算し、特色インキの吸収係数及び散乱係数を求める、請求項３から６の何れかに記載の色予測装置。
【請求項８】
前記インキが複数の原色インキを複数混合した特色インキの場合、特色インキに含まれる原色インキの比率に応じて、原色インキの前記指令網点面積率と、前記濃度階調領域の出現率との対応関係を示す前記対応関数の各々を加算し、特色インキの前記指令網点面積率と、前記濃度階調領域の出現率との前記対応関係を示す対応関数を求める、請求項３から７の何れかに記載の色予測装置。
【請求項９】
プロセッサを備えるコンピュータが実行する色予測方法であって、
前記コンピュータが、
指令網点面積率と、前記指令網点面積率でインキを印刷することにより形成される網点が含む濃度階調領域と、前記濃度階調領域の各々の出現率との関係に基づき、供給される前記指令網点面積率に対応する濃度階調領域及び前記濃度階調領域の出現率とを求める濃度階調出現率算出過程と、
濃度階調領域及び前記濃度階調領域の各々の出現率と、分光反射率との対応関係を示す計算モデルにより、前記指令網点面積率により形成された網点による再現色を示す分光反射率を、出現する前記濃度階調領域及び前記濃度階調領域の前記出現率に基づいて算出する分光データ算出過程と、
デバイス値と測色値の関係が記述されている色特性記述ファイルを読み込み、デバイス値に対する測色値を計算する測色値算過程とを含む処理を実行する、色予測方法。
【請求項１０】
コンピュータに実行させるプログラムであって、
前記コンピュータを、
指令網点面積率と、前記指令網点面積率でインキを印刷することにより形成される網点が含む濃度階調領域と、前記濃度階調領域の各々の出現率との関係に基づき、供給される前記指令網点面積率に対応する前記濃度階調領域及び前記濃度階調領域の出現率とを求める濃度階調出現率算出手段と、
前記濃度階調領域及び前記濃度階調領域の各々の出現率と、分光反射率との対応関係を示す計算モデルにより、前記指令網点面積率により形成された網点による再現色を示す分光反射率を、出現する前記濃度階調領域及び前記濃度階調領域の前記出現率に基づいて算出する分光データ算出手段と、
デバイス値と測色値の関係が記述されている色特性記述ファイルを読み込み、デバイス値に対する測色値を計算する測色値算手段として機能させるための、プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、印刷機が印刷する再現色のカラーマッチングにおいて、この再現色を予測する色予測装置、色予測方法およびプログラムに関する。
【背景技術】
【０００２】
例えば、カラーマッチングにおいて面積変調により再現色を調整する場合、与えられた網点面積率で各色が印刷される部分の色を予測する方法として、例えば、マーレイ・デービス式を用いる方法（非特許文献１）が知られている。
【０００３】
上述したモデルは、網点の大きさが異なっていても、網点が原色単体ベタの均一な厚さのインキにより印刷媒体表面に印刷されていることを前提としている。このため、網点におけるドットのエッジにおける光の散乱による構成分を正確に表現することができていない。
【０００４】
この対応として、修正ユールニールセン式を用いる方法（非特許文献１）や、網点におけるドットのエッジにおける光の散乱構成分を、より正確に考慮するモデルもある（例えば、特許文献１参照）。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００６−３２９７５３号公報
【非特許文献】
【０００６】
【非特許文献１】日本色彩学会編、「カラーイメージング（色彩科学講座）」、朝倉書店、２００４年、Ｐ．１４８
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、従来技術は、オフセット印刷などの面積変調階調表現、すなわち網点の大きさが異なっていても、網点が原色単体ベタの均一な厚さとなる印刷に対応した色予測方法である。
したがって、網点の大きさが異なり、かつ網点の大きさにより網点を形成するインキの厚さが不均一となるグラビア印刷などのような面積変調階調表現と濃度変調階調表現が混在した印刷方式に対応させることはできない。
すなわち、従来の色予測は、オフセット印刷に代表される面積変調階調表現の印刷方式に対応しており、印刷媒体表面の網点の１つ１つが網点面積率に係わらず、同じ色、濃度であることを前提とした物理モデルにより行われている。
【０００８】
一方、上述したグラビア印刷のように、面積変調階調表現に対し、濃度変調階調表現が加わった場合、網点面積率が変わると、網点の色及び濃度も変わることになる。
このため、特許文献１などの従来例においては、網点面積によらずに色、濃度が一律である場合の色予測を前提としているため、面積変調階調表現に対し、濃度変調階調表現が加わったグラビア印刷等の色予測を精度よく行うことが困難である。
【課題を解決するための手段】
【０００９】
上述した課題を解決するための本発明の一局面は、指令網点面積率と、指令網点面積率でインキを印刷することにより形成される網点が含む濃度階調領域と、濃度階調領域の各々の出現率との関係に基づき、供給される指令網点面積率に対応する濃度階調領域及び濃度階調領域の出現率とを求める濃度階調出現率算出部と、濃度階調領域及び濃度階調領域の各々の出現率と、分光反射率との対応関係を示す第１の計算モデルにより、指令網点面積率により形成された網点による再現色を示す分光反射率を、出現する濃度階調領域及び濃度階調領域の出現率に基づいて算出する分光データ算出部と、デバイス値と測色値の関係が記述されている色特性記述ファイルを読み込み、デバイス値に対する測色値を計算する測色値算出部とを備える、色予測装置である。
【００１０】
また、分光データ算出部が、濃度階調領域及び濃度階調領域の各々の出現率と、濃度階調領域の分光反射率から求めた分光光学濃度との対応関係を示す第２の計算モデルにより、指令網点面積率により形成された網点による再現色を示す分光光学濃度を算出する光学濃度算出部と、分光光学濃度から、再現色の分光反射率を算出する分光反射率算出部とを有してもよい。
【００１１】
また、第１の計算モデルが、指令網点面積率によるインキが印刷された印刷部分の網点が含む濃度階調領域毎に、濃度階調領域の分光反射率と濃度階調領域の出現率とを乗算した結果を加算して、指令網点面積率の印刷部分の分光反射率を算出する式であり、指令網点面積率と、濃度階調領域の出現率との対応関係を示す対応関数が予め書き込まれて記憶されている出現率テーブルデータベースをさらに、含み、分光データ算出部が、出現率テーブルデータベースから、指令網点面積率に対応する濃度階調領域の分光反射率と、濃度階調領域の出現率を読み出し、第１の計算モデルにより、指令網点面積率に対応する印刷部分の分光反射率を求めてもよい。
【００１２】
また、第２の計算モデルが、指令網点面積率によるインキが印刷された印刷部分の網点が含む濃度階調領域毎に、濃度階調領域の分光光学濃度と濃度階調の出現率とを乗算した結果を加算して、指令網点面積率の印刷部分の分光光学濃度を算出する式であり、指令網点面積率と、濃度階調領域の出現率との対応関係を示す対応関数が予め書き込まれて記憶されている出現率テーブルデータベースをさらに、含み、分光データ算出部が、出現率テーブルデータベースから、指令網点面積率に対応する濃度階調領域の分光光学濃度と、濃度階調領域の出現率を読み出し、第２の計算モデルにより、指令網点面積率に対応する印刷部分の分光光学濃度を求めてもよい。
【００１３】
また、インキをベタで印刷した印刷部分を測定した分光反射率から吸収係数及び散乱係数を求め、求めた吸収係数及び散乱係数により、濃度階調領域のインキの厚さに対応する濃度階調領域毎の濃度階調分光反射率をクベルカ・ムンクの式により算出し、印刷部分の実測の分光反射率と、濃度階調領域の分光反射率を用いて第１の計算モデルにより算出する分光反射率との差が最小となるように、第１の計算モデルにおける出現率をフィッティングし、得られた濃度階調領域毎の出現率を用い、濃度階調毎の指令網点面積率と、濃度階調領域の出現率との対応関係を示す対応関数が求められてもよい。
【００１４】
また、インキをベタで印刷した印刷部分を測定した分光反射率から吸収係数及び散乱係数を求め、求めた吸収係数及び散乱係数により、濃度階調領域のインキの厚さに対応する濃度階調領域毎の濃度階調分光光学濃度をクベルカ・ムンクの式により算出し、印刷部分の実測の分光反射率から求めた分光光学濃度と、濃度階調分光光学濃度を用いて第２の計算モデルにより算出する分光光学濃度との差が最小となるように、計算モデルにおける出現率をフィッティングし、得られた濃度階調領域毎の出現率を用い、濃度階調毎の指令網点面積率と、濃度階調領域の出現率との対応関係を示す対応関数が求められてもよい。
【００１５】
また、インキが複数の原色のインキを複数混合した特色インキの場合、特色インキに含まれる原色インキの比率に応じて吸収係数及び散乱係数の各々を加算し、特色インキの吸収係数及び散乱係数を求めてもよい。
【００１６】
また、インキが複数の原色インキを複数混合した特色インキの場合、特色インキに含まれる原色インキの比率に応じて、原色インキの指令網点面積率と、濃度階調領域の出現率との対応関係を示す対応関数の各々を加算し、特色インキの指令網点面積率と、濃度階調領域の出現率との対応関係を示す対応関数を求めてもよい。
【００１７】
また、本発明の他の局面は、プロセッサを備えるコンピュータが実行する色予測方法であって、コンピュータが、指令網点面積率と、指令網点面積率でインキを印刷することにより形成される網点濃度階調領域と、濃度階調領域の各々の出現率との関係に基づき、供給される指令網点面積率に対応する濃度階調領域及び濃度階調領域の出現率とを求める濃度階調出現率算出過程と、濃度階調領域及び濃度階調領域の各々の出現率と、分光反射率との対応関係を示す計算モデルにより、指令網点面積率により形成された網点による再現色を示す分光反射率を、出現する濃度階調領域及び濃度階調領域の出現率に基づいて算出する分光データ算出過程と、デバイス値と測色値の関係が記述されている色特性記述ファイルを読み込み、デバイス値に対する測色値を計算する測色値算過程とを含む処理を実行する、色予測方法である。
【００１８】
また、本発明の他の局面は、コンピュータに実行させるプログラムであって、コンピュータを、指令網点面積率と、指令網点面積率でインキを印刷することにより形成される網点濃度階調領域と、濃度階調領域の各々の出現率との関係に基づき、供給される指令網点面積率に対応する濃度階調領域及び濃度階調領域の出現率とを求める濃度階調出現率算出手段と、濃度階調領域及び濃度階調領域の各々の出現率と、分光反射率との対応関係を示す計算モデルにより、指令網点面積率により形成された網点による再現色を示す分光反射率を、出現する濃度階調領域及び濃度階調領域の出現率に基づいて算出する分光データ算出手段と、デバイス値と測色値の関係が記述されている色特性記述ファイルを読み込み、デバイス値に対する測色値を計算する測色値算手段として機能させるためのプログラムである。
【発明の効果】
【００１９】
以上説明したように、本発明によれば、網点を形成する濃度階調領域の分光反射率と濃度階調領域の各々の出現率との対応関係により、指令網点面積率により形成された網点における分光反射率を求めているため、面積変調階調と、濃度変調階調との双方の色表現を有するグラビア印刷等における印刷物の色予測を行うことができる。また、上述した方法により、特色の任意の指令網点面積率に対する分光反射率が求められるのと同時に、デバイス値と測色値の関係が色特性記述ファイルとして取得できる場合には、色特性記述ファイルを参照してその再現色を予測することができるため、例えば、プロセス色で印刷される画像部と、特色で印刷される平網やグラデーション等がレイアウトされた画像に対する色予測を、簡便に実現することができる。
【図面の簡単な説明】
【００２０】
【図１】本発明の第１の実施形態による色予測装置１の構成例を示すブロック図である。
【図２】指令網点面積率と、指令網点面積率で形成される網点における濃度階調領域の出現率との対応関係を説明する図である。
【図３】本発明の第１の実施形態による色予測装置１による、濃度階調出現率テーブルを算出する際の動作例を示すフローチャートである。
【図４】本発明の第１の実施形態による色予測装置１による、特色インキの任意の指令網点面積率に対する分光反射率を予測する処理の動作例を示すフローチャートである。
【図５】本発明の第１の実施形態による色予測装置１による、画像の色を予測する処理の動作例を示すフローチャートである。
【図６】プロセスインキで印刷される部分と、特色インキで印刷される部分から構成される画像の例を示す図である。
【図７】本発明の第２の実施形態による色予測装置１Ａの構成例を示すブロック図である。
【図８】本発明の第２の実施形態による色予測装置１Ａによる、濃度階調出現率テーブルを算出する際の動作例を示すフローチャートである。
【図９】本発明の第２の実施形態による色予測装置１Ａによる、特色インキの任意の指令網点面積率に対する分光反射率を予測する処理の動作例を示すフローチャートである。
【図１０】本発明の第２の実施形態による色予測装置１Ａによる、画像の色を予測する処理の動作例を示すフローチャートである。
【発明を実施するための形態】
【００２１】
＜第１の実施形態＞
本発明においては、グラビア印刷により印刷媒体に対してインキを印刷する際、階調度を示す指令網点面積率に対応した網点が印刷媒体表面に形成される。このグラビア印刷における網点は、インキが印刷される印刷媒体表面の面積（面積変調階調表現）と、印刷されるインキの厚さ（濃度変調階調表現）とが、指令網点面積率に応じて変化する。例えば、グラビア印刷において形成される網点は、山の構造に相似しているということもでき、大きな山ほど裾野が広く高さも高く、この大きな山に比較して小さな山ほど裾野が狭く高さも低い。すなわち、グラビア印刷の網点は、インキが印刷される面積だけでなく、印刷されるインキの厚さも、指令網点面積率に応じて変化する。
【００２２】
このため、本実施形態においては、グラビア印刷における網点の構造を、指令網点面積率に対応して、インキが印刷された網点の印刷媒体表面の面積と、印刷される網点におけるインキの厚さとを表す計算モデルを用いる。本実施形態においては、印刷媒体（例えば、用紙）に対してグラビア印刷により印刷した単色の各階調度の色予測を行う際、指令網点面積率で生成される各階調度の網点が、複数の濃度階調領域（等高線と同様な構成）から生成されている網点形状をモデル化した計算モデルを用いる。
【００２３】
以下、本発明の第１の実施形態について、図面を参照して説明する。
図１は、本発明の第１の実施形態による色予測装置１の構成例を示すブロック図である。この図１において、色予測装置１は、入力部１０１、濃度階調分光反射率算出部１０２、吸収係数・散乱係数データベース１０３、濃度階調出現率テーブル算出部１０４、測定分光反射率データベース１０５、濃度階調出現率テーブルデータベース１０６、分光反射率予測部１０７（分光データ算出部）（濃度階調出現率算出部１０７１、分光反射率算出部１０７２を含む）、測色値算出部１０８、プロセス色測色値算出部１０９、色特性記述ファイル群１１０及び出力部１１１を備えている。なお、図１に示すブロックおよびその接続関係は例示であって、接続関係を図示しないブロック間でデータの受け渡しが行われてもよい。例えば、後述するように吸収係数・散乱係数データベース１０３、測定分光反射率データベース１０５、濃度階調出現率テーブルデータベース１０６が記憶する各種データは、入力部１０１が受け付けたデータに基づいて設定、更新されてもよい。また、各ブロックが表す機能単位も例示であって、これらの機能の一部または全部は、例えば記憶手段、プロセッサ等を備えたコンピュータによって実現可能である。
【００２４】
入力部１０１は、例えば、外部のコンピュータと接続されており、ユーザが設定するインキの各原色の指令網点面積率の指定値、特色ベタの分光反射率、特色に含まれる原色インキの比率などのデータを入力する。また、本実施形態による色予測装置１で画像の色を予測する場合には、入力部１０１は色を予測する画像のレイアウトデータを受け付け、適切な処理部に、画像の画素値を受け渡す機能を提供する。
また、入力部１０１は、キーボードあるいはタッチパネルなどの入力手段を有し、この入力手段から入力されるユーザが設定するインキの各原色の指令網点面積率の指定値などのデータを、色予測装置１内部の各部に出力する構成としても良い。また、入力部１０１は、反射物体の分光反射率を測定できる入力手段を有し、特色ベタの分光反射率など、必要な測定データを直接測定し、色予測装置１内部の各部に出力する構成としても良い。
【００２５】
吸収係数・散乱係数データベース１０３には、原色のインキ毎に、指令網点面積率が１００％、すなわちベタで印刷媒体（例えば、コート紙などの用紙）に印刷された印刷部分から求めたインキの着色層の散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）が、上記外部のコンピュータなどにより予め書き込まれて記憶されている。
【００２６】
散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）を求める際には、白地及び黒地の各々の印刷媒体に対して、それぞれ原色のインキをベタでグラビア印刷して上記印刷部分を作成する。そして、白地及び黒地の各々の印刷媒体表面にベタで印刷されたインキの着色層の分光反射率をそれぞれ測定する。ここで、白地の印刷媒体表面における印刷部分の分光反射率を白地測定分光反射率とし、黒地の印刷媒体表面における印刷部分の分光反射率を黒地測定分光反射率とする。
この求めた白地測定分光反射率及び黒地測定分光反射率から、原色のインキをベタで印刷媒体に印刷した印刷部分における散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）を求める。また、散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）は、所定の波長の範囲で複数の波長λで求められている。
【００２７】
測定分光反射率データベース１０５には、複数の指令網点面積率で印刷媒体に印刷された印刷部分の測定分光反射率Ｒｓ（λ）が上記外部のコンピュータなどにより予め書き込まれて記憶されている。例えば、この印刷部分はｍ段階の指令網点面積率で印刷されたステップチャートとなっている。測定分光反射率Ｒｓ（λ）を求める際、原色のインキを用いて、複数の指令網点面積率の網点を、実際の印刷に用いる印刷媒体に対してそれぞれグラビア印刷する。ここで、ｓは指令網点面積率である。そして、印刷媒体に印刷された印刷部分の分光反射率を、指令網点面積率の網点を有する印刷部分毎に測定する。また、測定分光反射率データベース１０５には、実際の印刷に用いられる印刷媒体の下地分光反射率Ｒ０（λ）が、上述した測定分光反射率Ｒｓ（λ）と同様に、上記外部のコンピュータなどにより予め書き込まれて記憶されている。
【００２８】
濃度階調分光反射率算出部１０２は、入力部１０１より供給される原色のインキの種類に応じて、吸収係数・散乱係数データベース１０３から、対応する原色のインキの吸収係数Ｋ（λ）及び散乱係数Ｓ（λ）を読み出す。また、濃度階調分光反射率算出部１０２は、測定分光反射率データベース１０５から、入力部１０１より供給される印刷媒体の種類に応じて、印刷媒体の下地分光反射率Ｒ０（λ）を読み出す。
【００２９】
そして、濃度階調分光反射率算出部１０２は、吸収係数・散乱係数データベース１０３から読み出した原色のインキの吸収係数Ｋ（λ）及び散乱係数Ｓ（λ）と、測定分光反射率データベース１０５から読み出した印刷媒体の下地分光反射率Ｒ０（λ）と、所定の濃度階調の膜厚係数Ｘｍ（後述するように、膜厚係数Ｘｍは濃度階調領域毎に変更する）との各々を、下記の（１）式（クベルカ・ムンクの式）に代入し、濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）の各々を算出する。
以下の（１）式において、ａ（λ）は、散乱係数Ｓ（λ）と吸収係数Ｋ（λ）とを加算し、加算結果を散乱係数Ｓ（λ）により除算した数値である。ｂ（λ）は、ａ（λ）を二乗した数値から１を減算し、減算結果の平方根を計算した数値である。
【００３０】
【数１】
2015161523.tif 000003
【００３１】
本実施形態において、上述したクベルカ・ムンクの式、すなわち（１）式における膜厚係数Ｘｍは、印刷媒体に原色のインキがベタで印刷された印刷部分を基に、印刷部分の濃度階調を示している。膜厚係数Ｘｍは任意に設定することができ、例えば、最もインキの膜厚が厚いベタの印刷部分を１００％として膜厚係数を１とし、これを基準にして、濃度階調領域の膜厚の段階ごとに膜厚係数Ｘｍを設定する。例えば、指令網点面積率の示す濃度階調領域の膜厚の段階が５段階であれば、ｍ＝１、２、３、４、５であり、膜厚係数Ｘｍは、それぞれの濃度階調領域の膜厚の段階に対応して、Ｘ１＝１．０、Ｘ２＝０．８、Ｘ３＝０．６、Ｘ４＝０．４、Ｘ５＝０．２とする。
【００３２】
すでに述べたように、この膜厚係数Ｘｍと、下地の印刷媒体の下地分光反射率Ｒ０（λ）と、吸収係数Ｋ（λ）と、散乱係数Ｓ（λ）とをクベルカ・ムンクの式である（１）式に代入し、指令網点面積率の網点が含む濃度階調領域の分光反射率として、濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）の各々を算出する。この濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）が、後述する計算モデルで用いる、網点を構成する複数の濃度階調領域の各々の分光反射率として用いられる。
【００３３】
また、濃度階調分光反射率算出部１０２は、入力部１０１より供給される特色インキに含まれる原色インキの種類に応じて、吸収係数・散乱係数データベース１０３から、対応する原色のインキの吸収係数Ｋ（λ）及び散乱係数Ｓ（λ）を読み出す。濃度階調分光反射率算出部１０２は、読み出したそれぞれの原色の吸収係数Ｋ（λ）、散乱係数Ｓ（λ）に対して、入力部１０１より供給される特色インキに含まれる原色インキの比率を乗じて加算することにより、特色インキの吸収係数Ｋｔ（λ）及び散乱係数Ｓｔ（λ）を算出する。例えば、特色インキが２つの原色インキを混ぜ合わせたものである場合、以下の（２）式により、特色インキの吸収係数Ｋｔ（λ）及び散乱係数Ｓｔ（λ）を算出する。
【００３４】
【数２】
2015161523.tif 000004
【００３５】
上記（２）式において、係数α及び係数βの各々は、原色のインキ＃１及び原色のインキ＃２それぞれの混ぜ合わせる比率を示す係数である。原色のインキ＃１の吸収係数Ｋ１（λ）に対して係数αを乗算し、原色のインキ＃２の吸収係数Ｋ２（λ）に対して係数βを乗算し、各乗算結果を加算した数値を、特色インキの吸収係数Ｋｔ（λ）としている。同様に、原色のインキ＃１の散乱係数Ｓ１（λ）に対して係数αを乗算し、原色のインキ＃２の散乱係数Ｓ２（λ）に対して係数βを乗算し、各乗算結果を加算した数値を、特色インキの散乱係数Ｓｔ（λ）としている。また、濃度階調分光反射率算出部１０２は、測定分光反射率データベース１０５から、入力部１０１より供給される印刷媒体の種類に応じて、印刷媒体の下地分光反射率Ｒ０（λ）を読み出す。
【００３６】
そして、濃度階調分光反射率算出部１０２は、特色の吸収係数Ｋｔ（λ）、散乱係数Ｓｔ（λ）、下地の印刷媒体の下地分光反射率Ｒ０（λ）を（１）式に代入して得られるＲＫＭ（λ）と、入力部１０１より供給される特色インキベタの分光反射率Ｒｔ（λ）との平均二乗誤差ＲＭＳＥが最も小さくなるように、特色ベタの膜厚係数Ｘｔを求める。平均二乗誤差ＲＭＳＥは以下の（３）式で求める。
【００３７】
【数３】
2015161523.tif 000005
【００３８】
特色インキの濃度階調毎の膜厚係数Ｘｍは、このＸｔを基準にして、濃度階調領域の膜厚の段階ごとに膜厚係数Ｘｍを設定する。例えば、濃度階調領域の膜厚の段階が５段階であれば、ｍ＝１、２、３、４、５であり、膜厚係数Ｘｍは、それぞれの濃度階調領域の膜厚の段階に対応して、Ｘ１＝１．０Ｘｔ、Ｘ２＝０．８Ｘｔ、Ｘ３＝０．６Ｘｔ、Ｘ４＝０．４Ｘｔ、Ｘ５＝０．２Ｘｔとする。
【００３９】
すでに述べたように、この膜厚係数Ｘｍと、下地の印刷媒体の下地分光反射率Ｒ０（λ）と、吸収係数Ｋｔ（λ）と、散乱係数Ｓｔ（λ）とをクベルカ・ムンクの式である（１）式に代入し、指令網点面積率の網点が含む濃度階調領域の分光反射率として、濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）の各々を算出する。この濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）が、後述する計算モデルで用いる、網点を構成する複数の濃度階調領域の各々の分光反射率として用いられる。
【００４０】
濃度階調出現率テーブル算出部１０４は、濃度階調分光反射率算出部１０２から、原色の濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）の各々を読み込む。また、濃度階調出現率テーブル算出部１０４は、指令網点面積率毎の測定分光反射率Ｒｓ（λ）の各々を、測定分光反射率データベース１０５から読み込む。
そして、濃度階調出現率テーブル算出部１０４は、以下の（４）式（第１の計算モデル）に対して、濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）の各々を代入し、後述する処理により算出分光反射率
【００４１】
【数４】
2015161523.tif 000006
【００４２】
を算出する。
【００４３】
【数５】
2015161523.tif 000007
【００４４】
ここで、濃度階調出現率テーブル算出部１０４は、（４）式において、出現率（用紙の印刷部分における網点を構成する各濃度階調の濃度階調領域の面積の比率）の数値を変更しつつ、算出分光反射率
【００４５】
【数６】
2015161523.tif 000008
【００４６】
を算出する。そして、濃度階調出現率テーブル算出部１０４は、算出された算出分光反射率
【００４７】
【数７】
2015161523.tif 000009
【００４８】
の各々と、測定分光反射率Ｒｓ（λ）との平均二乗誤差ＲＭＳＥを、指令網点面積率ｓ毎に所定の波長範囲において求める。濃度階調出現率テーブル算出部１０４は、算出分光反射率
【００４９】
【数８】
2015161523.tif 000010
【００５０】
と測定分光反射率Ｒｓ（λ）との平均二乗誤差が最も小さくなる濃度階調領域の各々の出現率を求める。ここで、ｓは指令網点面積率である。平均二乗誤差ＲＭＳＥは（５）式で求める。
【００５１】
そして、濃度階調出現率テーブル算出部１０４は、濃度階調領域の各々の出現率により、指令網点面積率と各濃度階調領域の出現率との対応関係を示す対応関数である出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）を求める。ここで、濃度階調出現率テーブル算出部１０４は、濃度階調領域毎に、得られた出現率を用いて、この出現率を指定網点面積率ｓの２次関数などをフィッティング（近似）し、出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）を求めてもよい。濃度階調出現率テーブル算出部１０４は、求めた指令網点面積率と各濃度階調領域の出現率との対応関係を示す対応関数である出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）を、濃度階調出現率テーブルデータベース１０６に書き込んで記憶させる。
【００５２】
【数９】
2015161523.tif 000011
【００５３】
上述したように、濃度階調出現率テーブル算出部１０４は、（４）式の計算モデルを用いて、指令網点面積率ｓによるインキが印刷された印刷部分において、網点が含む濃度階調毎に、濃度階調の濃度階調分光反射率Ｒｉｍ（λ）と濃度階調の出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）とを乗算した結果を加算し、さらにいずれの濃度階調領域も出現しない下地部分の出現率（１−ａ１（ｓ）−ａ２（ｓ）−…−ａｍ（ｓ））に下地分光反射率Ｒ０（λ）を乗じた結果を加算して、指令網点面積率ｓの印刷部分の分光反射率
【００５４】
【数１０】
2015161523.tif 000012
【００５５】
（算出分光反射率）を算出する。
【００５６】
図２は、指令網点面積率ｓと、指令網点面積率ｓで形成される網点における濃度階調領域の出現率ａとの対応関係を説明する図である。図２の（ａ）において、横軸が指令網点面積率ｓを示し、縦軸が出現率ａを示している。出現率関数ａｍ（ｓ）により、指令網点面積率ｓにおける各濃度階調領域の出現率ａｍ（ｓ）を求めることができる。濃度階調出現率テーブル算出部１０４は、すでに説明したように濃度階調領域毎の出現率関数ａｍ（ｓ）の各々を、（５）式で求めた指令網点面積率ｓに対応した濃度階調領域毎の出現率を２次関数の近似式として関数化することにより求める。この図２の（ａ）は、ｍ＝４の場合、すなわち４階調の濃度階調の濃度階調領域の各々の出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、ａ４（ｓ）の各々を示している。
【００５７】
また、図２の（ｂ）は、図２の（ａ）における指令網点面積率ｓ１、ｓ２、ｓ３、ｓ４の各々における網点の平面視における形状を示す図である。指令網点面積率ｓ１においては、濃度階調領域Ｐ１のみが形成されている。また、指令網点面積率ｓ２においては、濃度階調領域Ｐ１の内部に濃度階調領域Ｐ２が形成される。指令網点面積率ｓ３においては、濃度階調領域Ｐ２のみが形成されている。また、指令網点面積率ｓ４においては、濃度階調領域Ｐ３の内部に濃度階調領域Ｐ４が形成される。上述したように、本実施形態においては、（４）式の計算モデルを用いて、この図２の（ｂ）に示す網点構造として、グラビア印刷における網点の構造をモデル化している。
【００５８】
図１に戻り、濃度階調出現率テーブル算出部１０４は、求めた濃度階調領域の指令網点面積率ｓにおける各濃度階調領域の出現率を示す出現率関数ａｍ（ｓ）を、濃度階調出現率テーブルデータベース１０６に対し、濃度階調領域の濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）と対応させて書き込んで記憶させる。同様に、濃度階調出現率テーブル算出部１０４は、他の原色のインキの各々に対しても、濃度階調領域の指令網点面積率ｓにおける各濃度階調領域の出現率を示す出現率関数ａｍ（ｓ）を、濃度階調出現率テーブルデータベース１０６に対して、書き込んで記憶させる。
【００５９】
次に、分光反射率予測部１０７は、濃度階調出現率算出部１０７１と、分光反射率算出部１０７２を備えている。
濃度階調出現率算出部１０７１は、入力部１０１より供給される特色インキに含まれる原色インキの種類に応じて、濃度階調出現率テーブルデータベース１０６から、それぞれの原色に対応する濃度階調出現率テーブルを読み出す。そして、入力部より供給される特色インキに含まれる原色インキの比率を、それぞれの原色インキに対応する濃度階調出現率テーブルに乗じて加算することにより、特色インキの濃度階調出現率テーブルを算出する。濃度階調出現率算出部１０７１は、上述のように求めた特色の濃度階調出現率テーブルに、入力部１０１より供給される特色の指令網点面積率ｓを代入することにより、指令網点面積率に対する濃度階調出現率ａ１、ａ２、…、ａｍを算出する。
【００６０】
分光反射率算出部１０７２は、濃度階調分光反射率算出部１０２から、特色の濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）を読み出す。さらに、濃度階調出現率算出部１０７１から、指令網点面積率に対する濃度階調出現率ａ１、ａ２、…、ａｍを読み出す。また、測定分光反射率データベース１０５から、下地の印刷媒体の下地分光反射率Ｒ０（λ）を読み出す。分光反射率算出部１０７２は、これら読み出した値を、（４）式に代入し、指令網点面積率ｓにより形成された網点による再現色を示す特色の予測分光反射率
【００６１】
【数１１】
2015161523.tif 000013
【００６２】
を算出する。
【００６３】
測色値算出部１０８は、この予測分光反射率
【００６４】
【数１２】
2015161523.tif 000014
【００６５】
に基づき、観測光源を定めて、三刺激値ＸＹＺやＣＩＥＬＡＢ値等を算出し、結果を出力部１１１に出力する。
【００６６】
色特性記述ファイル群１１０には、デバイス値であるプロセスインキの掛け合わせと測色値の関係が記述されたファイルが、上記外部のコンピュータなどにより予め書き込まれて記憶されている。このファイルを記述するフォーマットの好適な例としては、公知のＩＣＣプロファイルフォーマットが挙げられる。本実施例では、色特性記述ファイルとしてＩＣＣプロファイルフォーマットで記述されたファイルを使用する。
【００６７】
プロセス色測色値算出部１０９は、入力部１０１より供給される印刷条件（印刷機、用紙）に対応する色特性記述ファイルを、色特性記述ファイル群１１０から読み出す。そして、入力部１０１から供給されるプロセスインキの掛け合わせに対して、色特性記述ファイルを参照することにより、再現される色の測色値を算出し、結果を出力部１１１に出力する。
【００６８】
次に図３は、第１の実施形態による色予測装置１による、濃度階調出現率テーブルを算出する際の動作例を示すフローチャートである。
ステップＳ１１：
入力部１０１は、ユーザが指定した原色インキの種類を、濃度階調分光反射率算出部１０２及び、濃度階調出現率テーブル算出部１０４に対して出力する。
【００６９】
ステップＳ１２：
濃度階調分光反射率算出部１０２は、使用する原色インキの吸収係数Ｋ（λ）、散乱係数Ｓ（λ）を、吸収係数・散乱係数データベース１０３から読み出す。
【００７０】
ステップＳ１３：
濃度階調分光反射率算出部１０２は、下地となる印刷媒体の下地分光反射率Ｒ０（λ）と、吸収係数Ｋ（λ）と、散乱係数Ｓ（λ）と、濃度階調の膜厚係数Ｘｍとの各々を（１）式に対して代入し、濃度階調分光反射率Ｒｉｍ（λ）を求める。
【００７１】
ステップＳ１４：
濃度階調出現率テーブル算出部１０４は、使用する原色インキのステップチャートの測定分光反射率Ｒｓ（λ）を、測定分光反射率データベース１０５から読み込む。
【００７２】
ステップＳ１５：
濃度階調出現率テーブル算出部１０４は、濃度階調分光反射率算出部１０２から濃度階調分光反射率Ｒｉｍ（λ）を読み込み、（４）式に代入して予測分光反射率
【００７３】
【数１３】
2015161523.tif 000015
【００７４】
を算出する。そして、ステップチャートの測定分光反射率Ｒｓ（λ）に対して、（５）式に示したＲＭＳＥが最小となるように、最適化により、濃度階調出現率ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）を求める。
【００７５】
ステップＳ１６：
濃度階調出現率テーブル算出部１０４は、ステップＳ１５で求めた、指令網点面積率ｓと濃度階調出現率ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）の関係をテーブルとして、濃度階調出現率テーブルデータベース１０６に書き出す。なお、書き込むテーブルの形式としては、離散的な指令網点面積率ｓと濃度階調出現率ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）の関係を一次元ルックアップテーブルの形式で書き出しても良いし、あるいは、濃度階調出現率を指令網点面積率の関数として多項式、例えば二次式で近似し、その係数を書き出すようにしても良い。
なお、図３に示したフローチャートの処理は、特色に使用される原色毎に１度行えば良く、色予測の度に行う必要はない。
【００７６】
次に、図４は、第１の実施形態による色予測装置１による、特色インキの任意の指令網点面積率に対する分光反射率を予測する処理の動作例を示すフローチャートである。
【００７７】
ステップＳ２１：
入力部１０１は、ユーザが指定した特色インキの種類、特色インキに含まれる原色インキの種類と比率、特色インキベタの測定分光反射率を、濃度階調分光反射率算出部１０２及び、分光反射率予測部１０７に対して出力する。
【００７８】
ステップＳ２２：
濃度階調分光反射率算出部１０２は、使用する特色インキに含まれる原色インキの吸収係数Ｋ（λ）及び散乱係数Ｓ（λ）を、吸収係数・散乱係数データベース１０３から読み出す。そして、特色の吸収係数Ｋｔ（λ）と、散乱係数Ｓｔ（λ）とを、読み出した原色の吸収係数、散乱係数に、それぞれの原色インキが混ぜ合わされる割合をかけて足し合わせることにより算出する。
例えば、使用する特色インキに含まれる原色インキが２色の場合、原色インキ＃１の吸収係数Ｋ１（λ）、散乱係数Ｓ１（λ）と、原色インキ＃２の吸収係数Ｋ２（λ）、散乱係数Ｓ２（λ）を読み出す。そして、原色インキ＃１がα、原色インキ＃２がβの割合で混合されている場合、原色のインキ＃１の吸収係数Ｋ１（λ）に対して係数αを乗算し、原色のインキ＃２の吸収係数Ｋ２（λ）に対して係数βを乗算し、各乗算結果を加算した数値を、特色インキの吸収係数Ｋｔ（λ）とし、同様に、原色のインキ＃１の散乱係数Ｓ１（λ）に対して係数αを乗算し、原色のインキ＃２の散乱係数Ｓ２（λ）に対して係数βを乗算し、各乗算結果を加算した数値を、特色インキの散乱係数Ｓｔ（λ）とする。
【００７９】
ステップＳ２３：
濃度階調分光反射率算出部１０２は、（１）式にステップＳ２２で求めた特色の吸収係数Ｋｔ（λ）と、散乱係数Ｓｔ（λ）と、下地となる印刷媒体の下地分光反射率Ｒ０（λ）を代入して、特色ベタの予測分光反射率ＲＫＭ（λ）を算出する。そして、（３）式により、予測分光反射率ＲＫＭ（λ）と、特色ベタの測定分光反射率Ｒｔ（λ）とのＲＭＳＥが最小となるように、最適化により特色インキベタの膜厚係数Ｘｔを算出する。
なお、本実施例では、特色ベタの分光反射率を入力し、予測分光反射率ＲＫＭ（λ）とのＲＭＳＥが最小となるように最適化により膜厚係数Ｘｔを求めているが、特色ベタの測色値を入力として受け取り、予測分光反射率ＲＫＭ（λ）から算出される測色値との色差が最小となるように最適化により膜厚係数Ｘｔを求める変形例も好適である。
【００８０】
ステップＳ２４：
濃度階調分光反射率算出部１０２は、下地となる印刷媒体の下地分光反射率Ｒ０（λ）と、特色の吸収係数Ｋｔ（λ）と、散乱係数Ｓｔ（λ）と、濃度階調の膜厚係数Ｘｍとの各々を（１）式に対して代入し、濃度階調分光反射率Ｒｉｍ（λ）を求める。なお、このとき膜厚係数Ｘｍは、ステップＳ２３で求めたベタの膜厚係数Ｘｔを基準として、濃度階調ごとに求めた値を用いる。
【００８１】
ステップＳ２５：
濃度階調出現率算出部１０７１は、使用する特色インキに含まれる原色インキの種類に応じて、対応する濃度階調出現率テーブルを、濃度階調出現率テーブルデータベース１０６から読み出す。そして、特色インキに含まれる原色インキの割合に応じて、これらテーブルを足し合わせることにより、特色インキの濃度階調出現率テーブルを算出する。
【００８２】
ステップＳ２６：
入力部１０１は、ユーザが指定した特色の指令網点面積率ｓを、分光反射率予測部１０７に対して出力する。
【００８３】
ステップＳ２７：
濃度階調出現率算出部１０７１は、ステップＳ２５で算出した特色インキの濃度階調出現率テーブルを、ステップＳ２６で指定された指令網点面積率ｓで参照し、各濃度階調の出現率を算出する。
【００８４】
ステップＳ２８：
分光反射率算出部１０７２は、濃度階調分光反射率算出部１０２から、濃度階調分光反射率Ｒｉｍ（λ）を読み出す。また、濃度階調出現率算出部１０７１から、各濃度階調の出現率を読み出す。さらに、測定分光反射率データベース１０５から、下地の印刷媒体の下地分光反射率Ｒ０（λ）を読み出す。これら読み出した値を、（４）式に代入して、指令網点面積率により形成された網点による再現色を示す予測分光反射率
【００８５】
【数１４】
2015161523.tif 000016
【００８６】
を算出する。
【００８７】
なお、同一の特色について、異なる指令網点面積率に対する分光反射率を算出する際は、図４に示すフローチャートのステップＳ２１からＳ２５は１度だけ行えば良く、異なる指令網点面積率の指定に対して、ステップＳ２６からＳ２８を繰り返せば良い。
【００８８】
次に、図５は、第１の実施形態による色予測装置１による、画像の色を予測する処理の動作例を示すフローチャートである。
【００８９】
ステップＳ３１：
入力部１０１は、ユーザが指定する画像データと、画像で使用される特色の情報（特色の種類、特色に含まれる原色の種類と割合、特色ベタの測定分光反射率）、プロセス部の印刷の情報（印刷機、用紙）を受け付ける。画像データは、図６の例に示すように、プロセスインキで印刷される画像部６０１と、特色インキで印刷される平網や文字６０２や、特色のグラデーション６０３で構成される画像の各画素の情報が記述されたデータである。
【００９０】
ステップＳ３２：
入力部１０１は、画像データの画素のデータを順次読み出す。そして、画素に指定された色が、プロセスインキで印刷される色か、特色で印刷される色かを判断する。画素に指定された色が、プロセスインキで印刷される色（プロセス色）であればステップＳ３３へ、特色で印刷される色であればステップＳ３５に分岐する。
【００９１】
ステップＳ３３：
プロセス色測色値算出部１０９は、入力部１０１から供給されるプロセス印刷の情報（印刷機、用紙）に基づいて、デバイス値であるプロセスの掛け合わせとそれに対する測色値が記述された色特性記述ファイルを、色特性記述ファイル群１１０から読み出す。
【００９２】
ステップＳ３４：
プロセス色測色値算出部１０９は、入力部１０１から供給される画素値（ＣＭＹＫ）に基づいて、色特性記述ファイルを参照し、対応するプロセス色の測色値を算出する。
【００９３】
ステップＳ３５：
画素が特色で印刷される色である場合、入力部１０１から供給される特色の情報（特色の種類、特色に含まれる原色の種類と割合、特色ベタの測定分光反射率）と、画素値により示される特色の指令網点面積率に従って、図４に示したフローチャートの手順により、特色の分光反射率を算出する。
【００９４】
ステップＳ３６：
ステップＳ３５で算出した分光反射率に基づいて、観察光源を定めることにより、ＸＹＺ値や、ＣＩＥＬＡＢ値などの測色値を算出する。
【００９５】
ステップＳ３７：
全ての画素について上記処理を繰り返すことにより、全ての画素が測色値の情報を持つ画像を取得することができる。
【００９６】
＜第２の実施形態＞
以下、本発明の第２の実施形態について、図面を参照して説明する。
図７は、本発明の第２の実施形態による色予測装置１Ａの構成例を示すブロック図である。この図７において、色予測装置１Ａは、入力部１０１、濃度階調分光光学濃度算出部１０２Ａ、吸収係数・散乱係数データベース１０３、濃度階調出現率テーブル算出部１０４Ａ、測定分光反射率データベース１０５、濃度階調出現率テーブルデータベース１０６Ａ、分光反射率予測部１１２（分光データ算出部）（濃度階調出現率算出部１１２１、分光光度濃度算出部１１２２、分光反射率算出部１１２３を含む）、測色値算出部１０８、プロセス色測色値算出部１０９、色特性記述ファイル群１１０及び出力部１１１を備えている。
【００９７】
濃度階調分光光学濃度算出部１０２Ａは、入力部１０１より供給される原色のインキの種類に応じて、吸収係数・散乱係数データベース１０３から、対応する原色のインキの吸収係数Ｋ（λ）及び散乱係数Ｓ（λ）を読み出す。また、濃度階調分光反射率算出部１０２Ａは、測定分光反射率データベース１０５から、入力部１０１より供給される印刷媒体の種類に応じて、印刷媒体の下地分光反射率Ｒ０（λ）を読み出す。
【００９８】
そして、濃度階調分光反射率算出部１０２Ａは、吸収係数・散乱係数データベース１０３から読み出した原色のインキの吸収係数Ｋ（λ）及び散乱係数Ｓ（λ）と、測定分光反射率データベース１０５から読み出した印刷媒体の下地分光反射率Ｒ０（λ）と、所定の濃度階調の膜厚係数Ｘｍ（後述するように、膜厚係数Ｘｍは濃度階調領域毎に変更する）との各々を、（１）式（クベルカ・ムンクの式）に代入し、濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）の各々を算出する。
【００９９】
また、濃度階調分光光学濃度算出部１０２Ａは、求めた濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）の各々を、（６）式により、濃度階調分光光学濃度Ｄｉ１（λ）、Ｄｉ２（λ）、Ｄｉ３（λ）、…、Ｄｉｍ（λ）にそれぞれ変換する。
【０１００】
【数１５】
2015161523.tif 000017
【０１０１】
本実施形態においても、第１の実施形態と同様に、（１）式（クベルカ・ムンクの式）における膜厚係数Ｘｍは、印刷媒体に原色のインキがベタで印刷された印刷部分を基に、印刷部分の濃度階調を示している。膜厚係数Ｘｍは任意に設定することができ、第１の実施形態と同様に、例えば、最もインキの膜厚が厚いベタの印刷部分を１００％として、の膜厚係数Ｘ１を１とし、これを基準にして、濃度階調領域の膜厚の段階ごとに膜厚係数Ｘｍを設定する。例えば、指令網点面積率の示す濃度階調領域の膜厚の段階が５段階であれば、ｍ＝１、２、３、４、５であり、膜厚係数Ｘｍは、それぞれの濃度階調領域の膜厚の段階に対応して、Ｘ１＝１．０、Ｘ２＝０．８、Ｘ３＝０．６、Ｘ４＝０．４、Ｘ５＝０．２とする。
【０１０２】
すでに述べたように、この膜厚係数Ｘｍと、下地の印刷媒体の下地分光反射率Ｒ０（λ）と、吸収係数Ｋ（λ）と、散乱係数Ｓ（λ）とを（１）式に代入し、指令網点面積率の網点が含む濃度階調領域の分光反射率として、濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）の各々を算出する。そして、濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）の各々を、（６）式により、濃度階調分光光学濃度Ｄｉ１（λ）、Ｄｉ２（λ）、Ｄｉ３（λ）、…、Ｄｉｍ（λ）にそれぞれ変換する。この濃度階調分光光学濃度Ｄｉ１（λ）、Ｄｉ２（λ）、Ｄｉ３（λ）、…、Ｄｉｍ（λ）が、後述する計算モデルで用いる、網点を構成する複数の濃度階調領域の各々の分光光学濃度として用いられる。
【０１０３】
また、濃度階調分光光学濃度算出部１０２Ａは、入力部１０１より供給される特色インキに含まれる原色インキの種類に応じて、吸収係数・散乱係数データベース１０３から、対応する原色のインキの吸収係数Ｋ（λ）及び散乱係数Ｓ（λ）を読み出す。濃度階調分光光学濃度算出部１０２Ａは、読み出したそれぞれの原色の吸収係数Ｋ（λ）、散乱係数Ｓ（λ）に対して、入力部１０１より供給される特色インキに含まれる原色インキの比率を乗じて加算することにより、特色インキの吸収係数Ｋｔ（λ）及び散乱係数Ｓｔ（λ）を算出する。例えば、特色インキが２つの原色インキを混ぜ合わせたものである場合、（２）式により、特色インキの吸収係数Ｋｔ（λ）及び散乱係数Ｓｔ（λ）を算出する。
【０１０４】
また、濃度階調分光光学濃度算出部１０２Ａは、測定分光反射率データベース１０５から、入力部１０１より供給される印刷媒体の種類に応じて、印刷媒体の下地分光反射率Ｒ０（λ）を読み出す。
【０１０５】
そして、濃度階調分光光学濃度算出部１０２Ａは、特色の吸収係数Ｋｔ（λ）、散乱係数Ｓｔ（λ）、下地の印刷媒体の下地分光反射率Ｒ０（λ）を（１）式に代入して得られるＲＫＭ（λ）と、入力部１０１より供給される特色インキベタの分光反射率Ｒｔ（λ）との平均二乗誤差ＲＭＳＥが最も小さくなるように、特色ベタの膜厚係数Ｘｔを求める。平均二乗誤差ＲＭＳＥは（３）式で求める。
【０１０６】
特色インキの濃度階調毎の膜厚係数Ｘｍは、このＸｔを基準にして、濃度階調領域の膜厚の段階ごとに膜厚係数Ｘｍを設定する。例えば、濃度階調領域の膜厚の段階が５段階であれば、ｍ＝１、２、３、４、５であり、膜厚係数Ｘｍは、それぞれの濃度階調領域の膜厚の段階に対応して、Ｘ１＝１．０Ｘｔ、Ｘ２＝０．８Ｘｔ、Ｘ３＝０．６Ｘｔ、Ｘ４＝０．４Ｘｔ、Ｘ５＝０．２Ｘｔとする。
【０１０７】
すでに述べたように、この膜厚係数Ｘｍと、下地の印刷媒体の下地分光反射率Ｒ０（λ）と、吸収係数Ｋｔ（λ）と、散乱係数Ｓｔ（λ）とをクベルカ・ムンクの式である（１）式に代入し、指令網点面積率の網点が含む濃度階調領域の分光反射率として、濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）の各々を算出する。そして、濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）の各々を、（６）式により、濃度階調分光光学濃度Ｄｉ１（λ）、Ｄｉ２（λ）、Ｄｉ３（λ）、…、Ｄｉｍ（λ）にそれぞれ変換する。この濃度階調分光光学濃度Ｄｉ１（λ）、Ｄｉ２（λ）、Ｄｉ３（λ）、…、Ｄｉｍ（λ）が、後述する計算モデルで用いる、網点を構成する複数の濃度階調領域の各々の分光光学濃度として用いられる。
【０１０８】
濃度階調出現率テーブル算出部１０４Ａは、濃度階調分光光学濃度算出部１０２Ａから、濃度階調分光光学濃度Ｄｉ１（λ）、Ｄｉ２（λ）、Ｄｉ３（λ）、…、Ｄｉｍ（λ）の各々を読み込む。また、濃度階調出現率テーブル算出部１０４Ａは、指令網点面積率毎の測定分光反射率Ｒｓ（λ）の各々を、測定分光反射率データベース１０５から読み込む。
そして、濃度階調出現率テーブル算出部１０４Ａは、（６）式を用いることにより、測定分光反射率Ｒｓ（λ）の各々を、測定光学濃度Ｄｓ（λ）にそれぞれ変換する。
次に、濃度階調出現率テーブル算出部１０４Ａは、以下の（７）式（第２の計算モデル）に対して、濃度階調光学濃度Ｄｉ１（λ）、Ｄｉ２（λ）、Ｄｉ３（λ）、…、Ｄｉｍ（λ）の各々を代入し、後述する処理により算出光学濃度
【０１０９】
【数１６】
2015161523.tif 000018
【０１１０】
を算出する。
【０１１１】
【数１７】
2015161523.tif 000019
【０１１２】
ここで、濃度階調出現率テーブル算出部１０４Ａは、すでに説明した（７）式において、出現率（用紙の印刷部分における網点を構成する各濃度階調の濃度階調領域の面積の比率）の数値を変更しつつ、算出光学濃度
【０１１３】
【数１８】
2015161523.tif 000020
【０１１４】
を算出する。そして、濃度階調出現率テーブル算出部１０４Ａは、算出された算出光学濃度
【０１１５】
【数１９】
2015161523.tif 000021
【０１１６】
の各々と、測定光学濃度Ｄｓ（λ）との平均二乗誤差、すなわち（８）式に示すＲＭＳＥを、指令網点面積率ｓ毎に所定の波長範囲において求める。濃度階調出現率テーブル算出部１０４Ａは、算出光学濃度
【０１１７】
【数２０】
2015161523.tif 000022
【０１１８】
と測定光学濃度Ｄｓ（λ）との平均二乗誤差が最も小さくなる濃度階調領域の各々の出現率を求める。ここで、ｓは指令網点面積率である。
【０１１９】
【数２１】
2015161523.tif 000023
【０１２０】
そして、濃度階調出現率テーブル算出部１０４Ａは、濃度階調領域の各々の出現率により、指令網点面積率と各濃度階調領域の出現率との対応関係を示す対応関数である出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）を求める。ここで、濃度階調出現率テーブル算出部１０４Ａは、濃度階調領域毎に、得られた出現率を用いて、この出現率を指定網点面積率ｓの２次関数などをフィッティング（近似）し、出現率関数を求めてもよい。濃度階調出現率テーブル算出部１０４Ａは、求めた指令網点面積率と各濃度階調領域の出現率との対応関係を示す対応関数である出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）を、濃度階調出現率テーブルデータベース１０６Ａに書き込んで記憶させる。出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）については、第１の実施形態と同様に、（７）式の計算モデルで用いる濃度階調領域の出現率を求める関数であり、図２の説明におけるように、グラビア印刷による網点の構造をモデル化するために用いている。
【０１２１】
上述したように、濃度階調出現率テーブル算出部１０４Ａは、（７）式の計算モデルを用いて、指令網点面積率によるインキが印刷された印刷部分において、網点が含む濃度階調毎に、濃度階調の濃度階調分光光学濃度Ｄｉｍ（λ）と濃度階調の出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）とを乗算した結果を加算して、指令網点面積率の印刷部分の分光光学濃度を算出する。
【０１２２】
また、濃度階調出現率テーブル算出部１０４Ａは、求めた濃度階調領域の指令網点面積率における各濃度階調領域の出現率を示す出現率関数ａｍ（ｓ）を、濃度階調出現率テーブルデータベース１０６Ａに対し、濃度階調領域の濃度階調分光光学濃度と対応させて書き込んで記憶させる。同様に、濃度階調出現率テーブル算出部１０４Ａは、他の原色のインキの各々に対しても、濃度階調領域の指令網点面積率における各濃度階調領域の出現率を示す出現率関数ａｍ（ｓ）を、濃度階調出現率テーブルデータベース１０６Ａに対して、書き込んで記憶させる。
【０１２３】
次に、分光反射率予測部１１２は、濃度階調出現率算出部１１２１と、分光光学濃度算出部１１２２、分光反射率算出部１１２３を備えている。
濃度階調出現率算出部１１２１は、入力部１０１より供給される特色インキに含まれる原色インキの種類に応じて、濃度階調出現率テーブルデータベース１０６Ａから、それぞれの原色に対応する濃度階調出現率テーブルを読み出す。そして、入力部１０１より供給される特色インキに含まれる原色インキの比率を、それぞれの原色インキに対応する濃度階調出現率テーブルに乗じて加算することにより、特色インキの濃度階調出現率テーブルを算出する。濃度階調出現率算出部１１２１は、上述のように求めた特色の濃度階調出現率テーブルに、入力部１０１より供給される特色の指令網点面積率を代入することにより、指令網点面積率に対する濃度階調出現率ａ１、ａ２、…、ａｍを算出する。
【０１２４】
分光光学濃度算出部１１２２は、濃度階調分光光学濃度算出部１０２Ａから、特色の濃度階調分光光学濃度Ｄｉ１（λ）、Ｄｉ２（λ）、Ｄｉ３（λ）、…、Ｄｉｍ（λ）を読み出す。さらに、濃度階調出現率算出部１１２１から、指令網点面積率に対する濃度階調出現率ａ１、ａ２、…、ａｍを読み出す。また、測定分光反射率データベース１０５から、下地の印刷媒体の下地分光反射率Ｒ０（λ）を読み出す。分光光学濃度算出部１１２２は、これら読み出した値を、（７）式に代入し、指令網点面積率により形成された網点による再現色を示す特色の予測分光光学濃度
【０１２５】
【数２２】
2015161523.tif 000024
【０１２６】
を算出する。
【０１２７】
分光反射率算出部１１２３は、分光光学濃度算出部１１２２が算出した予測分光光学濃度
【０１２８】
【数２３】
2015161523.tif 000025
【０１２９】
を、以下の（９）式を用い、予測分光反射率
【０１３０】
【数２４】
2015161523.tif 000026
【０１３１】
に変換する。
【０１３２】
【数２５】
2015161523.tif 000027
【０１３３】
次に、図８は、第２の実施形態による色予測装置１Ａによる、濃度階調出現率テーブルを算出する際の動作例を示すフローチャートである。
【０１３４】
ステップＳ４１：
入力部１０１は、ユーザが指定した原色インキの種類を、濃度階調分光光学濃度算出部１０２Ａ及び、濃度階調出現率テーブル算出部１０４Ａに対して出力する。
【０１３５】
ステップＳ４２：
濃度階調分光光学濃度算出部１０２Ａは、使用する原色インキの吸収係数Ｋ（λ）、散乱係数Ｓ（λ）を、吸収係数・散乱係数データベース１０３から読み出す。
【０１３６】
ステップＳ４３：
濃度階調分光光学濃度算出部１０２Ａは、下地となる印刷媒体の下地分光反射率Ｒ０（λ）と、吸収係数Ｋ（λ）と、散乱係数Ｓ（λ）と、濃度階調の膜厚係数Ｘｍとの各々を（１）式に対して代入し、濃度階調分光反射率Ｒｉｍ（λ）を求める。そして、求めた濃度階調分光反射率Ｒｉｍ（λ）から、（６）式により濃度階調分光光学濃度Ｄｉｍ（λ）を算出する。
【０１３７】
ステップＳ４４：
濃度階調出現率テーブル算出部１０４Ａは、使用する原色インキのステップチャートの測定分光反射率Ｒｓ（λ）を、測定分光反射率データベース１０５から読み込む。
【０１３８】
ステップＳ４５：
濃度階調出現率テーブル算出部１０４Ａは、濃度階調分光光学濃度算出部１０２Ａから濃度階調分光光学濃度Ｄｉｍ（λ）を読み込み、（７）式に代入して予測分光光学濃度
【０１３９】
【数２６】
2015161523.tif 000028
【０１４０】
を算出する。そして、ステップチャートの測定分光反射率Ｒｓ（λ）を、（６）式に代入して測定分光光学濃度Ｄｔ（λ）を求め、（８）式に示すＲＭＳＥが最小となるように、最適化により濃度階調出現率ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）を求める。
【０１４１】
ステップＳ４６：
濃度階調出現率テーブル算出部１０４Ａは、ステップＳ４５で求めた、指令網点面積率ｓと濃度階調出現率ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）の関係をテーブルとして、濃度階調出現率テーブルデータベース１０６Ａに書き出す。なお、書き込むテーブルの形式としては、離散的な指令網点面積率ｓと濃度階調出現率ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）の関係を一次元ルックアップテーブルの形式で書き出しても良いし、あるいは、濃度階調出現率を指令網点面積率の関数として多項式、例えば二次式で近似し、その係数を書き出すようにしても良い。
【０１４２】
なお、図８に示したフローチャートの処理は、特色に使用される原色毎に１度行えば良く、色予測の度に行う必要はない。
【０１４３】
次に、図９は、第２の実施形態による色予測装置１Ａによる、特色インキの任意の指令網点面積率に対する分光反射率を予測する処理の動作例を示すフローチャートである。
【０１４４】
ステップＳ５１：
入力部１０１は、ユーザが指定した特色インキの種類、特色インキに含まれる原色インキの種類と比率、特色インキベタの測定分光反射率を、濃度階調分光光学濃度算出部１０２Ａ及び、分光反射率予測部１１２に対して出力する。
【０１４５】
ステップＳ５２：
濃度階調分光光学濃度算出部１０２Ａは、使用する特色インキに含まれる原色インキの吸収係数Ｋ（λ）、散乱係数Ｓ（λ）を、吸収係数・散乱係数データベース１０３から読み出す。そして、特色の吸収係数Ｋｔ（λ）と、散乱係数Ｓｔ（λ）とを、読み出した原色の吸収係数、散乱係数に、それぞれの原色インキが混ぜ合わされる割合をかけて足し合わせることにより算出する。
例えば、使用する特色インキに含まれる原色インキが２色の場合、原色インキ＃１の吸収係数Ｋ１（λ）、散乱係数Ｓ１（λ）と、原色インキ＃２の吸収係数Ｋ２（λ）、散乱係数Ｓ２（λ）を読み出す。そして、原色インキ＃１がα、原色インキ＃２がβの割合で混合されている場合、原色のインキ＃１の吸収係数Ｋ１（λ）に対して係数αを乗算し、原色のインキ＃２の吸収係数Ｋ２（λ）に対して係数βを乗算し、各乗算結果を加算した数値を、特色インキの吸収係数Ｋｔ（λ）とし、同様に、原色のインキ＃１の散乱係数Ｓ１（λ）に対して係数αを乗算し、原色のインキ＃２の散乱係数Ｓ２（λ）に対して係数βを乗算し、各乗算結果を加算した数値を、特色インキの散乱係数Ｓｔ（λ）とする。
【０１４６】
ステップＳ５３：
濃度階調分光光学濃度算出部１０２Ａは、（１）式にステップＳ５２で求めた特色の吸収係数Ｋｔ（λ）と、散乱係数Ｓｔ（λ）と、下地となる印刷媒体の下地分光反射率Ｒ０（λ）を代入して、特色ベタの予測分光反射率ＲＫＭ（λ）を算出する。そして、（３）式により、予測分光反射率ＲＫＭ（λ）と、特色ベタの測定分光反射率Ｒｔ（λ）とのＲＭＳＥが最小となるように、最適化により特色インキベタの膜厚係数Ｘｔを算出する。
なお、本実施例では、特色ベタの分光反射率を入力し、予測分光反射率ＲＫＭ（λ）とのＲＭＳＥが最小となるように最適化により膜厚係数Ｘｔを求めているが、特色ベタの測色値を入力として受け取り、予測分光反射率ＲＫＭ（λ）から算出される測色値との色差が最小となるように最適化により膜厚係数を求める変形例も好適である。
【０１４７】
ステップＳ５４：
濃度階調分光光学濃度算出部１０２Ａは、下地となる印刷媒体の下地分光反射率Ｒ０（λ）と、特色の吸収係数Ｋｔ（λ）と、散乱係数Ｓｔ（λ）と、濃度階調の膜厚係数Ｘｍとの各々を（１）式に対して代入し、濃度階調分光反射率Ｒｉｍ（λ）を求める。なお、このとき膜厚係数Ｘｍは、ステップＳ５３で求めたベタの膜厚係数Ｘｔを基準として、濃度階調ごとに求めた値を用いる。求めた濃度階調分光反射率Ｒｉｍ（λ）を（６）式に代入し、濃度階調分光光学濃度Ｄｉｍ（λ）を求める。
【０１４８】
ステップＳ５５：
濃度階調出現率算出部１１２１は、使用する特色インキに含まれる原色インキの種類に応じて、対応する濃度階調出現率テーブルを、濃度階調出現率テーブルデータベース１０６Ａから読み出す。そして、特色インキに含まれる原色インキの割合に応じて、これらテーブルを足し合わせることにより、特色インキの濃度階調出現率テーブルを算出する。
【０１４９】
ステップＳ５６：
入力部１０１は、ユーザが指定した特色の指令網点面積率を、分光反射率予測部１１２に対して出力する。
【０１５０】
ステップＳ５７：
濃度階調出現率算出部１１２１は、ステップＳ５５で算出した特色インキの濃度階調出現率テーブルを、ステップＳ５６で指定された指令網点面積率で参照し、各濃度階調の出現率を算出する。
【０１５１】
ステップＳ５８：
分光光学濃度算出部１１２２は、濃度階調分光光学濃度算出部１０２Ａから、濃度階調分光反射率Ｄｉｍ（λ）を読み出す。また、濃度階調出現率算出部１１２１から、各濃度階調の出現率を読み出す。さらに、測定分光反射率データベース１０５から、下地の印刷媒体の下地分光反射率Ｒ０（λ）を読み出す。これら読み出した値を、（７）式に代入して、指令網点面積率により形成された網点による再現色を示す予測分光光学濃度
【０１５２】
【数２７】
2015161523.tif 000029
【０１５３】
を算出する。
【０１５４】
ステップＳ５９：
分光反射率算出部１１２３は、ステップＳ５８で算出した予測分光光学濃度
【０１５５】
【数２８】
2015161523.tif 000030
【０１５６】
を、以下の（８）式を用い、予測分光反射率
【０１５７】
【数２９】
2015161523.tif 000031
【０１５８】
に変換する。
【０１５９】
なお、同一の特色について、異なる指令網点面積率に対する分光反射率を算出する際は、図９に示すフローチャートのステップＳ５１からＳ５５は１度だけ行えば良く、異なる指令網点面積率の指定に対して、ステップＳ５６からＳ５９を繰り返せば良い。
【０１６０】
次に、図１０は、第２の実施形態による色予測装置１Ａによる、画像の色を予測する処理の動作例を示すフローチャートである。
【０１６１】
ステップＳ６１：
入力部１０１は、ユーザが指定する画像データと、画像で使用される特色の情報（特色の種類、特色に含まれる原色の種類と割合、特色ベタの測定分光反射率）、プロセス部の印刷の情報（印刷機、用紙）を受け付ける。画像データは、図６の例に示すように、プロセスインキで印刷される画像部６０１と、特色インキで印刷される平網や文字６０２や、特色のグラデーション６０３で構成される画像の各画素の情報が記述されたデータである。
【０１６２】
ステップＳ６２：
入力部１０１は、画像データの画素のデータを順次読み出す。そして、画素に指定された色が、プロセスインキで印刷される色か、特色で印刷される色かを判断する。画素に指定された色が、プロセスインキで印刷される色（プロセス色）であればステップＳ６３へ、特色で印刷される色であればＳ６５に分岐する。
【０１６３】
ステップＳ６３：
プロセス色測色値算出部１０９は、入力部１０１から供給されるプロセス印刷の情報（印刷機、用紙）に基づいて、デバイス値であるプロセスの掛け合わせとそれに対する測色値が記述された色特性記述ファイルを、色特性記述ファイル群１１０から読み出す。
【０１６４】
ステップＳ６４：
プロセス色測色値算出部１０９は、入力部１０１から供給される画素値（ＣＭＹＫ）に基づいて、色特性記述ファイルを参照し、対応するプロセス色の測色値を算出する。
【０１６５】
ステップＳ６５：
画素が特色で印刷される色である場合、入力部１０１から供給される特色の情報（特色の種類、特色に含まれる原色の種類と割合、特色ベタの測定分光反射率）と、画素値により示される特色の指令網点面積率に従って、図９に示したフローチャートの手順により、特色の分光反射率を算出する。
【０１６６】
ステップＳ６６：
ステップＳ６５で算出した分光反射率に基づいて、観察光源を定めることにより、ＸＹＺ値や、ＣＩＥＬＡＢ値などの測色値を算出する。
【０１６７】
ステップＳ６７：
全ての画素について上記処理を繰り返すことにより、全ての画素が測色値の情報を持つ画像を取得することができる。
【０１６８】
なお、本実施例では、デバイス値と測色値の関係を色特性記述ファイルとして取得可能な例として、プロセスインキによる印刷を取り上げたが、適用範囲はこれに限定されない。例えば、使用頻度の高い数色（２〜３色程度など）の特色について、色特性記述ファイルを作成しておき、このファイルに記述された特色に関しては、色特性記述ファイルを参照して色を予測する方法をとることも可能である。
【０１６９】
なお、本実施例では、印刷としてグラビア印刷を対象として説明したが、対象とする印刷方式はこれに限定されない。本発明で採用した色予測モデルは、面積変調階調と、濃度変調階調との双方の色表現を有するグラビア印刷等をモデル化するために好適であるのと同時に、基本的に面積変調階調による色表現を有するオフセット印刷でも、網点の周辺部で発生する光学的ドットゲインを精度良くモデル化することができるため、オフセット印刷に適用しても精度良く再現色を予測することができる。
【０１７０】
また、本実施例では、画像として、プロセスインキで印刷される箇所と、特色インキで印刷される箇所に重なりのない画像を例として扱ったが、適用範囲はこれに限定されない。プロセスインキで印刷される箇所と、特色インキで印刷される箇所に重なりがある画像の場合、例えば、本実施例のごとく、プロセス部分に対する測色値と特色部分に対する測色値をそれぞれ算出し、重なり合う部分では、それら測色値を合成することにより、重ね合わせの色を予測する。例えば、測色値がＸＹＺ値で与えられるとき、それらの値を用紙のＸＹＺ値で割った値の常用対数をＸＹＺ濃度として算出し、プロセス部分、特色部分それぞれのＸＹＺ濃度を加算して求めた合成のＸＹＺ濃度を、再度ＸＹＺに戻すことで、合成後のＸＹＺ値を算出する。ここでは単純な合成方法を例としてあげたが、より複雑な合成方法を採用しても良い。
【０１７１】
なお、本発明における図１の色予測装置１及び、図７の色予測装置１Ａの機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより色予測の制御を行ってもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。
また、「コンピュータシステム」は、ホームページ提供環境（あるいは表示環境）を備えたＷＷＷシステムも含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムが送信された場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリ（ＲＡＭ）のように、一定時間プログラムを保持しているものも含むものとする。
【０１７２】
また、上記プログラムは、このプログラムを記憶装置等に格納したコンピュータシステムから、伝送媒体を介して、あるいは、伝送媒体中の伝送波により他のコンピュータシステムに伝送されてもよい。ここで、プログラムを伝送する「伝送媒体」は、インターネット等のネットワーク（通信網）や電話回線等の通信回線（通信線）のように情報を伝送する機能を有する媒体のことをいう。また、上記プログラムは、前述した機能の一部を実現するためのものであっても良い。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であっても良い。また、本発明は、プロセッサを備えるコンピュータが実行する色予測方法としてとらえることができる。
【産業上の利用可能性】
【０１７３】
本発明は、印刷機等が印刷する再現色のカラーマッチングにおける再現色の予測に有用である。
【符号の説明】
【０１７４】
１、１Ａ            色予測装置
１０１              入力部
１０２              濃度階調分光反射率算出部
１０２Ａ            濃度階調分光光学濃度算出部
１０３              吸収係数・散乱係数データベース
１０４、１０４Ａ    濃度階調出現率テーブル算出部
１０５              測定分光反射率データベース
１０６、１０６Ａ    濃度階調出現率テーブルデータベース
１０７、１１２      分光反射率予測部
１０８              測色値算出部
１０９              プロセス色測色値算出部
１１０              色特性記述ファイル群
１１１              出力部
１０７１、１１２１  濃度階調出現率算出部
１０７２            分光反射率算出部
１１２２            分光光学濃度算出部
１１２３            分光反射率算出部
【図１】
2015161523.tif 000032
【図２】
2015161523.tif 000033
【図３】
2015161523.tif 000034
【図４】
2015161523.tif 000035
【図５】
2015161523.tif 000036
【図６】
2015161523.tif 000037
【図７】
2015161523.tif 000038
【図８】
2015161523.tif 000039
【図９】
2015161523.tif 000040
【図１０】
2015161523.tif 000041