(57)【要約】【課題】印刷機において面積変調階調や濃度変調階調により印刷した印刷物の色に、簡易な印刷デバイスで印刷する印刷物の色が一致するように、印刷機の面積変調階調や濃度変調階調によるインキ調整に対応した印刷デバイスの制御値を調整する印刷色調整システムを提供する。【解決手段】本発明の印刷色調整システムは、面積変調階調や濃度変調階調により印刷した印刷物に用いたインキの各々のベタ印刷の測色値から、面積変調階調や濃度変調階調の指令値から印刷される印刷物の色の予測値を求める印刷色予測モデルを用い、印刷物を印刷する第１印刷機のカラープロファイルを生成するカラープロファイル生成部と、カラープロファイルで、印刷物の画像データのドットの指令値に対応する色の予測値を求め、この予測値を印刷物の色の色情報とし、色情報を第２印刷機のカラープロファイルにより制御値に変換して、第２印刷機に対して出力する印刷制御部とを備える。【選択図】図１
2019016924.tif 000002
【特許請求の範囲】
【請求項１】
面積変調階調や濃度変調階調により印刷した印刷物に用いたインキの各々のベタ印刷の測色値から、面積変調階調や濃度変調階調の指令値から印刷される前記印刷物の色の予測値を求める印刷色予測モデルを用い、前記印刷物を印刷する第１印刷機のカラープロファイルを生成するカラープロファイル生成部と、
前記カラープロファイルにより、前記印刷物の画像データにおける各ドットの前記指令値に対応する色の予測値を求め、当該予測値を印刷物の色の色情報とし、当該色情報を第２印刷機のカラープロファイルにより制御値に変換して、前記第２印刷機に対して出力する印刷制御部と
を備えることを特徴とする印刷色調整システム。
【請求項２】
前記印刷色予測モデルにおける予測色の予測の計算に用いるパラメータの調整を行い、複数の段階に調整された前記パラメータに対応し、複数の前記印刷色予測モデルを生成して、複数の前記カラープロファイルを生成させる調整制御部をさらに備え、
前記印刷制御部が、前記カラープロファイルの各々から得られる制御値それぞれを、前記第２印刷機に出力する
ことを特徴とする請求項１に記載の印刷色調整システム。
【請求項３】
前記印刷色予測モデルが、前記インキを重ね刷りした際の前記予測値の予測に、インキの各々のベタ印刷の吸収係数及び散乱係数を用いて、前記印刷物の各ドットにおける分光反射率をクベルカ・ムンクの式により求める処理を含んでいる
ことを特徴とする請求項２に記載の印刷色調整システム。
【請求項４】
前記クベルカ・ムンクの式におけるベタ印刷の前記インキの厚さを前記パラメータとして用いる
ことを特徴とする請求項３に記載の印刷色調整システム。
【請求項５】
前記クベルカ・ムンクの式におけるベタ印刷の前記インキの色彩値を前記パラメータとして用いる
ことを特徴とする請求項３または請求項４に記載の印刷色調整システム。
【請求項６】
前記指令値の中間調における前記予測色の階調特性を前記パラメータとして用いる
ことを特徴とする請求項３から請求項５のいずれか一項に記載の印刷色調整システム。
【請求項７】
前記インキを重ね刷りした際の前記予測色を変更する際、他のインキに重ねる前記インキの分光反射率を求める前記クベルカ・ムンクの式におけるベタ印刷の前記インキの厚さを前記パラメータとして用いる
ことを特徴とする請求項３から請求項６のいずれか一項に記載の印刷色調整システム。
【請求項８】
前記調整制御部が、前記前記印刷物の版ずれを調整する際、前記印刷物の前記インキの版の各々の位置をずらし、
前記印刷制御部が、位置がずらされた前記版の重なりにおいて、前記印刷物の各ドットの色の予測値を求める
ことを特徴とする請求項２から請求項７のいずれか一項に記載の印刷色調整システム。
【請求項９】
前記調整制御部が、前記画像データの前記予測値に対してノイズを重畳させる
ことを特徴とする請求項２から請求項８のいずれか一項に記載の印刷色調整システム。
【請求項１０】
面積変調階調や濃度変調階調により印刷した印刷物に用いたインキの各々のベタ印刷の測色値を測定する測色過程と、
前記測色値から、面積変調階調や濃度変調階調の指令値から印刷される前記印刷物の色の予測値を求める印刷色予測モデルを用い、前記印刷物を印刷する第１印刷機のカラープロファイルを生成する過程と、
前記カラープロファイルにより、前記印刷物の画像データにおける各ドットの前記指令値に対応する色の予測値を求め、当該予測値を印刷物の色の色情報とし、当該色情報を第２印刷機のカラープロファイルにより制御値に変換して、前記第２印刷機に対して出力する過程と、
前記第１印刷機で印刷された前記印刷物の色と、前記第２印刷機で印刷された印刷物との色を比較し、前記印刷色予測モデルにおける調整する必要のあるパラメータを選択する過程と、
前記パラメータを多段階に変更し、複数の前記印刷色予測モデルを生成し、当該印刷色予測モデルの各々に対応する、前記カラープロファイルそれぞれを生成する過程と、
前記カラープロファイルの各々から得られる色情報に対応する制御値それぞれにより、前記第２印刷機に前記画像データを印刷させる過程と
を含むことを特徴とする印刷色調整方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、インキを網点で印刷した印刷物の色に、インクジェットプリンタやデジタル印刷機で出力する印刷物の色を調整する印刷色調整システムおよび印刷色調整方法に関する。
【背景技術】
【０００２】
グラビア印刷、スクリーン印刷及びオフセット印刷などにおいて、面積変調（面積変調階調表現）や濃度変調（濃度変調階調表現）により再現色を調整し、印刷媒体に対して所定の画像の印刷が行われる。
これらの印刷は、多量の枚数を印刷する大量印刷の場合に適しており、少量の枚数の印刷にはコストの問題で適していない。
このため、印刷物の枚数が少量足りなくなった場合に、上記印刷で増刷を行うとコストがかかるため、所定の枚数を余分に印刷し、在庫として保管する必要がある。
【０００３】
しかし、在庫として保管する場合、余分に印刷した在庫が使われなくなった場合、在庫のためのコストが無駄になる。
また、在庫を保管する場所を確保する必要となり、複数種類の印刷物の在庫を保管するためには、広い保管場所を用意しなければならず、多くの無駄なスペースが使われる。
【０００４】
そのため、少量の増刷が必要になった場合、印刷物の出力機として、インクジェットプリンタやデジタル印刷機などの簡易な印刷デバイスが用いられる。
このとき、大量印刷された印刷物と、簡易な印刷デバイスの印刷物との色が一致している必要があり、一致するためのカラーマネージメントが行われる。すなわち、大量印刷を行う印刷機のカラープロファイル（面積変調階調表現を所定の基準の色情報（たとえば、ＣＩＥＬＡＢ値）に変換）と、簡易な印刷デバイスのカラープロファイル（基準の色情報を印刷デバイスの制御値に変換）との各々を作成し、大量印刷された印刷物と、簡易な印刷デバイスの印刷物との色を一致させる（例えば、特許文献１及び特許文献２を参照）。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特表２００７−５１６６６３号公報
【特許文献２】特開２００６−３２９７５３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、印刷機における印刷のばらつきにより印刷物の色が上記カラープロファイルからずれていたり、あるいは意図的に色が変えられていた場合、カラープロファイルを用いて印刷した簡易な印刷デバイスの印刷物と、印刷機にて印刷した印刷物との色が一致しない。
また、特色に関しては、カラープロファイルが作成されていない場合が多く、あったとしても上述した理由により、簡易な印刷デバイスの印刷物と、印刷機にて印刷した印刷物との色が一致させることは困難である。
【０００７】
また、簡易な印刷デバイスが印刷する色を、印刷デバイスの制御値（たとえば、ＣＭＹＫ値）を調整し、印刷機の印刷物の色に合わせることが考えられる。
しかし、複数の特色で印刷されている場合、いずれかの特色の色を調整しようとし、例えばＣＭＹＫ値を調整することにより、すべての特色の色が変化してしまい、いずれかの特色の色のみを変更することができず、全体の色が印刷物の色とずれてしまう。
【０００８】
本発明は、このような状況に鑑みてなされたもので、印刷機において面積変調階調表現や濃度変調階調により印刷した印刷物の色に、簡易な印刷デバイスにより印刷する印刷物の色が一致するように、印刷機の調整方法に対応する、印刷デバイスの制御値を調整する印刷色調整システムおよび印刷色調整方法を提供する。
【課題を解決するための手段】
【０００９】
本発明の印刷色調整システムは、面積変調階調や濃度変調階調により印刷した印刷物に用いたインキの各々のベタ印刷の測色値から、面積変調階調や濃度変調階調の指令値から印刷される前記印刷物の色の予測値を求める印刷色予測モデルを用い、前記印刷物を印刷する第１印刷機のカラープロファイルを生成するカラープロファイル生成部と、前記カラープロファイルにより、前記印刷物の画像データにおける各ドットの前記指令値に対応する色の予測値を求め、当該予測値を印刷物の色の色情報とし、当該色情報を第２印刷機のカラープロファイルにより制御値に変換して、前記第２印刷機に対して出力する印刷制御部とを備えることを特徴とする。
【００１０】
本発明の印刷色調整システムは、前記印刷色予測モデルにおける前記予測色の予測の計算に用いるパラメータの調整を行い、複数の段階に調整された前記パラメータに対応し、複数の前記印刷色予測モデルを生成して、複数の前記カラープロファイルを生成させる調整制御部をさらに備え、前記印刷制御部が、前記カラープロファイルの各々から得られる制御値それぞれを、前記第２印刷機に出力することを特徴とする。
【００１１】
本発明の印刷色調整システムは、前記印刷色予測モデルが、前記インキを重ね刷りした際の前記予測値の予測に、インキの各々のベタ印刷の吸収係数及び散乱係数を用いて、前記印刷物の各ドットにおける分光反射率をクベルカ・ムンクの式により求める処理を含んでいることを特徴とする。
【００１２】
本発明の印刷色調整システムは、前記クベルカ・ムンクの式におけるベタ印刷の前記インキの厚さを前記パラメータとして用いることを特徴とする。
【００１３】
本発明の印刷色調整システムは、前記クベルカ・ムンクの式におけるベタ印刷の前記インキの色彩値を前記パラメータとして用いることを特徴とする。
【００１４】
本発明の印刷色調整システムは、前記指令値の中間調における前記予測色の階調特性を前記パラメータとして用いることを特徴とする。
【００１５】
本発明の印刷色調整システムは、前記インキを重ね刷りした際の前記予測色を変更する際、他のインキに重ねる前記インキの分光反射率を求める前記クベルカ・ムンクの式におけるベタ印刷の前記インキの厚さを前記パラメータとして用いる
ことを特徴とする。
【００１６】
本発明の印刷色調整システムは、前記調整制御部が、前記前記印刷物の版ずれを調整する際、前記印刷物の前記インキの版の各々の位置をずらし、前記印刷制御部が、位置がずらされた前記版の重なりにおいて、前記印刷物の各ドットの色の予測値を求めることを特徴とする。
【００１７】
本発明の印刷色調整システムは、前記調整制御部が、前記画像データの前記予測値に対してノイズを重畳させることを特徴とする。
【００１８】
本発明の印刷色調整方法は、面積変調階調や濃度変調階調により印刷した印刷物に用いたインキの各々のベタ印刷の測色値を測定する測色過程と、前記測色値から、面積変調階調や濃度変調階調の指令値から印刷される前記印刷物の色の予測値を求める印刷色予測モデルを用い、前記印刷物を印刷する第１印刷機のカラープロファイルを生成する過程と、前記カラープロファイルにより、前記印刷物の画像データにおける各ドットの前記指令値に対応する色の予測値を求め、当該予測値を印刷物の色の色情報とし、当該色情報を第２印刷機のカラープロファイルにより制御値に変換して、前記第２印刷機に対して出力する過程と、前記第１印刷機で印刷された前記印刷物の色と、前記第２印刷機で印刷された印刷物との色を比較し、前記印刷色予測モデルにおける調整する必要のあるパラメータを選択する過程と、前記パラメータを多段階に変更し、複数の前記印刷色予測モデルを生成し、当該印刷色予測モデルの各々に対応する、前記カラープロファイルそれぞれを生成する過程と、前記カラープロファイルの各々から得られる色情報に対応する制御値それぞれにより、前記第２印刷機に前記画像データを印刷させる過程とを含むことを特徴とする。
【発明の効果】
【００１９】
以上説明したように、本発明によれば、印刷機において面積変調階調や濃度変調階調により印刷した印刷物の色に、簡易な印刷デバイスにより印刷する印刷物の色が一致するように、印刷機の面積変調階調や濃度変調階調によるインキ調整に対応する、印刷デバイスの制御値を調整する印刷色調整システムおよび印刷色調整方法を提供することが可能となる。
【図面の簡単な説明】
【００２０】
【図１】本発明の第１の実施形態による印刷色調整システムの構成例を示すブロック図である。
【図２】図１における色予測モデル生成部１２の構成例を示すブロック図である。
【図３】指令網点面積率と、指令網点面積率で形成される網点における濃度階調領域の出現率との対応関係を説明する図である。
【図４】特色インキを構成する原色インキの配合比を決定し、特色インキの吸収係数Ｋｔ（λ）、散乱係数Ｓｔ（λ）を算出するフローチャートである。
【図５】下地としてのインキの上に重ねて印刷されるインキの濃度階調分光反射率Ｒｉｍ（λ）の算出を説明する図である。
【図６】特色インキの重ね合わせによりノイゲバウア原色を算出する処理を示すフローチャートである。
【図７】下地のインキの濃度階調領域と、下地のインキ（特色インキまたは原色インキ）に重ねて印刷されるインキの濃度階調領域との重なりの出現率の割合の算出を説明する図である。
【図８】図７における領域Ｑ１から領域Ｑ９のそれぞれの出現率の算出結果を示すテーブルの図である。
【図９】本実施形態における色予測モデル生成部１２及びカラープロファイル生成部１３によるカラープロファイルの生成処理の動作例を説明するフローチャートである。
【図１０】本実施形態における印刷色調整システムによる印刷機の印刷物の色に対し、印刷デバイスの印刷する色を同様とする調整処理の動作例を説明するフローチャートである。
【発明を実施するための形態】
【００２１】
本発明においては、印刷機（商業印刷機）でグラビア印刷、スクリーン印刷及びオフセット印刷などの手法で印刷媒体に対してインキを網点形状に印刷して画像の色を再現する網点印刷物と、同様な色味で、同一の画像データをインクジェットプリンタやデジタル印刷機など簡易な印刷機（印刷デバイス）で印刷する際、網点形状をモデル化した色予測モデルのパラメータを調整することにより、インクジェットプリンタの色調整を行う。
以下、本実施形態においては、グラビア印刷を対象として説明を行う。
グラビア印刷では、印刷媒体に対してインキを印刷する際、階調度を示す指令網点面積率に対応した網点が印刷媒体表面に形成される。網点は、インキが印刷される印刷媒体表面の面積（面積変調階調表現）と、印刷されるインキの厚さ（濃度変調階調表現）とが、指令網点面積率に応じて変化する。例えば、印刷において形成される網点は、山の構造に相似しているということもでき、大きな山ほど裾野が広く高さも高く、この大きな山に比較して小さな山ほど裾野が狭く高さも低い。すなわち、印刷における網点は、インキが印刷される面積だけでなく、印刷されるインキの厚さも、指令網点面積率に応じて変化する。
【００２２】
このため、本発明においては、印刷媒体（例えば、用紙）に対して印刷された網点印刷物における各ドットの単色あるいは特色の各階調度の色を、上記予測モデルによって色予測を行い、その予測された結果（例えばＣＩＥＬＡＢ値等）を、印刷するインクジェットプリンタのＩＣＣ（International Color Consortium）プロファイルなどにより、ＣＭＹＫ形式などの色データに変換してインクジェットプリンタに対して出力する。そして、すでに述べたように、インクジェットプリンタの印刷における色を調整する際、インクジェットプリンタの制御値を制御するのではなく、色予測モデルのパラメータを調整することにより、インクジェットプリンタの色調整を行う。
【００２３】
以下、図１における印刷色調整システムの構成例について、図面を参照して説明する。
図１は、本発明の第１の実施形態による印刷色調整システムの構成例を示すブロック図である。
図１において、印刷色調整システム１は、入力制御部１１、色予測モデル生成部１２、カラープロファイル生成部１３、調整制御部１４、制御値生成部１５、表示部１６、印刷データ出力部１７及び記憶部１８の各々を備えている。
入力制御部１１は、外部装置から供給される印刷物の印刷に用いるインキの版毎の画像データを、記憶部１８に書き込んで記憶させる。
また、ユーザは、印刷物の印刷に用いたインキのベタ印刷の測色値の測定を行い、不図示の入力手段（例えば、キーボードなど）により、各インキの測色値と刷り順とを印刷色調整システム１に対して入力する。これにより、入力制御部１１は、入力される各インキのベタ印刷の測色値と、各インキの刷り順とを、上記印刷物に対応させて記憶部１８に書き込んで記憶させる。
【００２４】
色予測モデル生成部１２は、印刷機に与えられる指令網点面積率と印刷部分の予測分光反射率Ｒ（λ）との対応を示す色予測モデルを生成する（後に詳細に説明）。ここで、指令網点面積率は、網点印刷物である参照印刷物を印刷する制御値として、グラビア印刷などを行う商業印刷機に与えられる数値である。すなわち、色予測モデル生成部１２は、各インキの刷り順に対応し、インキそれぞれの指令網点面積率の組合わせ毎の予測分光反射率Ｒ（λ）を求め、この予測分光反射率Ｒ（λ）から測色値を算出し、指令網点面積率と測色値との対応を示す色予測モデル（色予測テーブル）を生成する。
カラープロファイル生成部１３は、上記色予測モデル（色予測テーブル）に基づき、公知のＩＣＣプロファイルフォーマットに変換し、変換後のカラープロファイル（色予測プロファイルとも示す）を生成する。
【００２５】
調整制御部１４は、色予測モデルにおける予測分光反射率Ｒ（λ）を求める際に用いる式で用いられるパラメータを、複数の段階に調整する。そして、調整制御部１４は、調整した複数の段階のパラメータを色予測モデル生成部１２に対して出力し、それぞれのパラメータに対応する色予測モデルを生成させる。ここで、複数種類のパラメータを多段階で調整した場合、複数種類のパラメータの各段階の組合わせの数の色予測モデルが生成される。例えば、５種類のパラメータの各々を、それぞれ５段階に調整した場合、５５＝３１２５の組合わせとなり、３１２５個の色予測モデルが生成される。また、これにより、カラープロファイル生成部１３は、３１２５個の色予測モデルに基づき、３１２５個のカラープロファイルを生成する。
【００２６】
制御値生成部１５は、カラープロファイル生成部１３の生成したカラープロファイルを用いて、画像データの画素毎に、インキの各々に対する指令網点面積率により印刷される色の予測の測色値を求める。そして、制御値生成部１５は、予測された測色値に対応する制御値を、少量印刷をしようとしている印刷デバイスのデバイスカラープロファイルから抽出して求める。
【００２７】
表示部１６は、例えば液晶ディスプレイなどであり、処理に必要なデータの入力画面、また上記印刷デバイスに対する制御値を生成する処理に対する操作を行う操作画面などが表示される。
印刷データ出力部１７は、画像データの各々の画素に対する制御値を画像データ制御値として、印刷デバイスに対して出力して画像の印刷を実行させる。
【００２８】
＜印刷機のカラープロファイルの生成処理＞
以下に説明する色予測モデルは、印刷における網点の構造を、指令網点面積率に対応して、インキが印刷された網点の印刷媒体表面の面積と、印刷される網点におけるインキの厚さとを表す計算モデル（後述するコアフリンジモデル）を用いる。すなわち、印刷媒体（例えば、用紙）に対して印刷された単色あるいは特色の各階調度の色予測を行う際、指令網点面積率で生成される各階調度の網点が、複数の濃度階調領域（等高線と同様な構成）から生成されている網点形状をモデル化した計算モデルを用いる。
【００２９】
以下、図１における色予測モデル生成部１２について、図面を参照して説明する。
図２は、図１における色予測モデル生成部１２の構成例を示すブロック図である。この図２において、色予測モデル生成部１２は、入力部１０１、濃度階調分光反射率算出部１０２、濃度階調出現率算出部１０３、濃度階調出現率テーブルデータベース１０４、測定分光反射率データベース１０５、吸収係数・散乱係数データベース１０６、吸収係数・散乱係数算出部１０７、分光反射率予測部１０８、分光光学濃度予測部１１１、混合予測部１１４、特色インキ分光反射率算出部１２１、特色インキ配合比決定部１２２、特色インキ濃度階調出現率算出部１２３、一時記憶部１２５、予測パラメータデータベース１２６、近似色データベース１２８、特色インキ再現色算出部１２９、濃度階調分光光学濃度算出部２０２、濃度階調出現率算出部２０３及び濃度階調出現率テーブルデータベース２０４を備えている。
【００３０】
入力部１０１は、例えば、外部のコンピュータと入力制御部１１を介して接続されており、ユーザが設定する原色インキや特色インキの各原色の指令網点面積率の指令値などのデータを入力する。
また、入力部１０１は、入力制御部１１を介して入力されるユーザが設定するインキの各原色の指令網点面積率の指令値などのデータを、色予測モデル生成部１２内部の各部に出力する構成としても良い。
【００３１】
吸収係数・散乱係数データベース１０６には、原色インキ毎に、指令網点面積率が１００％、すなわちベタで印刷媒体（例えば、コート紙などの用紙）に印刷された印刷部分から求めたインキの着色層の散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）が、上記外部のコンピュータなどにより予め書き込まれて記憶されている。ここで、散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）を求める際、白地及び黒地の各々の印刷媒体に対して、それぞれ原色インキをベタで印刷して上記印刷部分を作成する。
そして、白地及び黒地の各々の印刷媒体表面にベタで印刷されたインキの着色層の分光反射率をそれぞれ測定する。ここで、白地の印刷媒体表面における印刷部分の分光反射率を白地測定分光反射率とし、黒地の印刷媒体表面における印刷部分の分光反射率を黒地測定分光反射率とする。
この求めた白地測定分光反射率及び黒地測定分光反射率から、原色インキをベタで印刷媒体に印刷した印刷部分における散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）を求める。また、散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）は、所定の波長の範囲で複数の波長λで求められている。  
【００３２】
近似色データベース１２８には、配合された原色インキの組合せが異なる複数種類の特色インキ（参照特色インキ）について、この特色インキの測色値と、原色の配合比が予め書き込まれて記憶されている。
【００３３】
測定分光反射率データベース１０５には、原色インキ毎の複数の指令網点面積率で印刷媒体に印刷された印刷部分の測定分光反射率Ｒｓ（λ）が上記外部のコンピュータなどにより予め書き込まれて記憶されている。例えば、この印刷部分はｍ段階の指令網点面積率で印刷されたステップチャートとなっている。測定分光反射率Ｒｓ（λ）を求める際、原色インキを用いて、複数の指令網点面積率の網点を、実際の印刷に用いる印刷媒体に対してそれぞれ印刷する。ここで、ｓは指令網点面積率である。そして、印刷媒体に印刷された印刷部分の分光反射率を、指令網点面積率の網点を有する印刷部分毎に測定する。また、測定分光反射率データベース１０５には、実際の印刷に用いられる印刷媒体の下地分光反射率Ｒ０（λ）が、上述した測定分光反射率Ｒｓ（λ）と同様に、上記外部のコンピュータなどにより予め書き込まれて記憶されている。
【００３４】
一時記憶部１２５は、本実施形態における各部の途中計算結果が記憶される記憶部であり、内部に散乱吸収係数テーブル、特色インキ配合比テーブル、特色インキ濃度階調出現率テーブル、ノイゲバウア原色テーブル、分光反射率テーブルなどが書き込まれて記憶されている。
予測パラメータデータベース１２６には、本実施形態においては混合予測部１１４で用いる後述する重み係数ｗが、記憶されている。この重み係数ｗは、教師となるデータを印刷して、重み係数により作成したモデルから得られた分光反射率と、印刷された教師データとしての分光反射率との誤差が最小となるようにして求められている。また、印刷色の分光反射率予測モデルとしてユール・ニールセン修正ノイゲバウアモデルを用いる場合には、ｎ値が設定されていても良い。
【００３５】
・分光反射率による濃度階調出現率の算出
濃度階調分光反射率算出部１０２は、入力部１０１より供給される原色インキの種類に応じて、吸収係数・散乱係数データベース１０６から、対応する原色インキの吸収係数Ｋ（λ）及び散乱係数Ｓ（λ）を読み出す。また、濃度階調分光反射率算出部１０２は、測定分光反射率データベース１０５から、入力部１０１より供給される原色インキの種類及び印刷媒体の種類に応じて、指令網点面積率毎の測定分光反射率Ｒｓ（λ）と、印刷媒体の下地分光反射率Ｒ０（λ）との各々を読み出す。
【００３６】
そして、濃度階調分光反射率算出部１０２は、吸収係数・散乱係数データベース１０６測定分光反射率データベース１０５から読み出した原色インキの吸収係数Ｋ（λ）及び散乱係数Ｓ（λ）と、印刷媒体の下地分光反射率Ｒ０（λ）と濃度階調の膜厚係数Ｘｍ（後述するように、膜厚係数Ｘｍは濃度階調領域毎に変更する）との各々を、下記の（１）式（クベルカ・ムンクの式）に代入し、濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）の各々を算出する。
以下の（１）式において、ａ（λ）は、散乱係数Ｓ（λ）と吸収係数Ｋ（λ）とを加算し、加算結果を散乱係数Ｓ（λ）により除算した数値である。ｂ（λ）は、ａ（λ）を二乗した数値から１を減算し、減算結果の平方根を計算した数値である。
【００３７】
【数１】
2019016924.tif 000003
【００３８】
本実施形態において、上述した（１）式に示すクベルカ・ムンク（Kubelka-Munk）の式、すなわち（１）式における印刷されたインキの膜厚係数Ｘｍは、印刷媒体に原色インキがベタで印刷された印刷部分を基にし、印刷部分の濃度階調を示す数値として用いる。すなわち、膜厚係数Ｘｍは任意に設定されており、例えば、最もインキの膜厚が厚いベタの印刷部分を膜厚が１００％として、上記膜厚係数を１とし、この１を濃度階調領域の膜厚の段階数ｍに対応させてｍ分割する。例えば、指令網点面積率の示す濃度階調領域の膜厚の段階が５段階であれば、ｍ＝１、２、３、４、５であり、膜厚係数Ｘｍは、それぞれの濃度階調領域の膜厚の段階に対応して、Ｘ１＝１．０、Ｘ２＝０．８、Ｘ３＝０．６、Ｘ４＝０．４、Ｘ５＝０．２とする。
【００３９】
膜厚係数Ｘｍを、すでに述べたように、下地の印刷媒体の下地分光反射率Ｒ０（λ）と、吸収係数Ｋ（λ）と、散乱係数Ｓ（λ）とともにクベルカ・ムンクの式である（１）式に代入し、それぞれの指令網点面積率の網点が含む濃度階調領域の分光反射率として、濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）の各々を算出する。この濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）が、後述する計算モデルで用いる、網点を構成する複数の濃度階調領域の各々の分光反射率として用いられる。
【００４０】
濃度階調出現率算出部１０３は、濃度階調分光反射率算出部１０２から、濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）の各々を読み込む。また、濃度階調出現率算出部１０３は、指令網点面積率毎の測定分光反射率Ｒｓ（λ）の各々を、測定分光反射率データベース１０５から読み込む。
そして、濃度階調出現率算出部１０３は、以下の（２）式（計算モデル）に対して、濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）の各々を代入し、後述する処理により算出分光反射率Ｒ’（ｓ，λ）を算出する。
【００４１】
【数２】
2019016924.tif 000004
【００４２】
ここで、濃度階調出現率算出部１０３は、以下の（３）式において、出現率（用紙の印刷部分における網点を構成する各濃度階調の濃度階調領域の面積の比率）の数値を変更しつつ、算出分光反射率Ｒ’（ｓ，λ）を算出する。そして、濃度階調出現率算出部１０３は、算出された算出分光反射率Ｒ’（ｓ，λ）の各々と、測定分光反射率Ｒｓ（λ）との平均二乗誤差ＲＭＳＥを、指令網点面積率毎に所定の波長範囲において求める。
濃度階調出現率算出部１０３は、算出分光反射率Ｒ’（ｓ，λ）と測定分光反射率Ｒｓ（λ）との平均二乗誤差が最も小さくなる濃度階調領域の各々の出現率を求める。ここで、ｓは指令網点面積率である。
【００４３】
そして、濃度階調出現率算出部１０３は、濃度階調領域の各々の出現率により、各濃度階調領域の出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）を求める。ここで、濃度階調出現率算出部１０３は、濃度階調領域毎に、得られた出現率を用いて、この出現率を指令網点面積率ｓの２次関数などをフィッティングし、出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）を求めてもよい。濃度階調出現率算出部１０３は、求めた各濃度階調領域の出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）を、濃度階調出現率テーブルデータベース１０４に書き込んで記憶させる。
【００４４】
【数３】
2019016924.tif 000005
【００４５】
上記（３）式において、濃度階調出現率算出部１０３は、波長λが３８０ｎｍから７３０ｎｍをｎ分割した刻み幅により、各波長λにおける誤差の二乗を加算した平均二乗誤差ＲＭＳＥを、指令網点面積率毎に求めている。
上述したように、濃度階調出現率算出部１０３は、（２）式の計算モデルを用いて、指令網点面積率によるインキが印刷された印刷部分において、網点が含む濃度階調毎に、濃度階調の濃度階調分光反射率Ｒｉｍ（λ）と濃度階調の出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）とを乗算した結果を加算して、指令網点面積率の印刷部分の分分光光反射率Ｒ’（ｓ，λ）（算出分光反射率）を算出する。
【００４６】
図３は、指令網点面積率と、指令網点面積率で形成される網点における濃度階調領域の出現率との対応関係を説明する図である。図３（ａ）において、横軸が指令網点面積率ｓを示し、縦軸が出現率ａを示している。出現率関数ａｍ（ｓ）により、指令網点面積率における各濃度階調領域の出現率を求めることができる。濃度階調出現率算出部１０３は、すでに説明したように濃度階調領域毎の出現率関数ａｍ（ｓ）の各々を、（３）式で求めた指令網点面積率に対応した濃度階調領域毎の出現率を２次関数の近似式として関数化することにより求める。この図３（ａ）は、ｍ＝４の場合、すなわち４階調の濃度階調の濃度階調領域の各々の出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、ａ４（ｓ）の各々を示している。
【００４７】
また、図３（ｂ）は、図３（ａ）における指令網点面積率ｓ１、ｓ２、ｓ３、ｓ４の各々における網点の平面視における形状を示す図である。指令網点面積率ｓ１においては、濃度階調領域Ｐ１のみが形成されている。また、指令網点面積率ｓ２においては、濃度階調領域Ｐ１の内部に濃度階調領域Ｐ２が形成される。指令網点面積率ｓ３においては、濃度階調領域Ｐ２のみが形成されている。また、指令網点面積率ｓ４においては、濃度階調領域Ｐ３の内部に濃度階調領域Ｐ４が形成される。上述したように、本実施形態においては、（２）式の計算モデルを用いて、この図２（ｂ）に示す網点構造として、グラビア印刷における網点の構造をモデル化している。
【００４８】
図２に戻り、濃度階調出現率算出部１０３は、求めた濃度階調領域の指令網点面積率における各濃度階調領域の出現率を示す出現率関数ａｍ（ｓ）を、濃度階調出現率テーブルデータベース１０４に対し、濃度階調領域の濃度階調分光反射率と対応させて書き込んで記憶させる。同様に、濃度階調出現率算出部１０３は、他の原色インキの各々に対しても、濃度階調領域の指令網点面積率における各濃度階調領域の出現率を示す出現率関数ａｍ（ｓ）を、濃度階調出現率テーブルデータベース１０４に対して、書き込んで記憶させる。
【００４９】
特色インキ分光反射率算出部１２１は、近似色データベース１２８における原色インキの組み合わせが異なる特色インキから、測色値が印刷物に用いたインキのベタ印刷の測色値と最も近い特色インキを読み出す。この近似色データベース１２８には特色インキの測色値と原色の配合比が記憶されており、読み出した特色インキの原色の配合比を、ベタ印刷の測色値を再現する特色インキの原色の配合比として設定する。
また、特色インキ分光反射率算出部１２１は、測定分光反射率データベース１０５から、印刷媒体の下地分光反射率Ｒ０（λ）を読み出す。
吸収係数・散乱係数算出部１０７は、吸収係数・散乱係数データベース１０６から、設定した配合比に含まれる原色インキを構成する原色インキ、例えば原色インキ＃１及び原色インキ＃２の各々の吸収係数Ｋ１（λ）及びＫ２（λ）と、散乱係数Ｓ１（λ）及びＳ２（λ）との各々を読み出す。
【００５０】
そして、吸収係数・散乱係数算出部１０７は、以下の（４）式により、原色インキを混色した特色インキの吸収係数Ｋｔ（λ）及び散乱係数Ｓｔ（λ）を算出する。この特色インキの場合、混ぜ合わせる原色インキの比率に応じて、以下の（４）式において、特色インキの散乱係数Ｓｔ（λ）及び吸収係数Ｋｔ（λ）を算出する。また、特色インキ分光反射率算出部１２１は、求めた特色インキの散乱係数Ｓｔ（λ）及び吸収係数Ｋｔ（λ）の各々を、一時記憶部１２５の散乱吸収係数テーブルに書き込んで記憶させる。
【００５１】
【数４】
2019016924.tif 000006
【００５２】
上記（４）式において、係数α及び係数βの各々は、原色インキ＃１及び原色インキ＃２それぞれの混ぜ合わせる比率を示す係数である。原色インキ＃１の吸収係数Ｋ１（λ）に対して係数αを乗算し、原色インキ＃２の吸収係数Ｋ２に対して係数βを乗算し、加算した数値を、特色インキの吸収係数Ｋｔ（λ）としている。同様に、原色インキ＃１の散乱係数Ｓ１（λ）に対して係数αを乗算し、原色インキ＃２の散乱係数Ｓ２（λ）に対して係数βを乗算し、加算した数値を、特色インキの散乱係数Ｓｔ（λ）としている。
【００５３】
これにより、特色インキ分光反射率算出部１２１は、（１）式に対して、吸収係数Ｋｔ（λ）及び散乱係数Ｓｔ（λ）と、印刷媒体の下地分光反射率Ｒ０（λ）と濃度階調の膜厚係数Ｘｍとの各々を代入し、特色インキの分光反射率ＲＫＭ（λ）を求める。
特色インキ再現色算出部１２９は、特色インキ分光反射率算出部１２１が算出した特色インキの分光反射率ＲＫＭ（λ）に対して観察環境の光源の分光分布と、標準観測者を設定し、測色値（例えばＬ*ａ*ｂ*値）に変換し、特色インキ配合比決定部１２２に対して出力する。
【００５４】
特色インキ配合比決定部１２２は、特色インキ再現色算出部１２９が求めた測色値と、ベタ印刷の測色値との色差を確認する。そして、特色インキ配合比決定部１２２は、この色差が予め設定された許容範囲内である場合、特色インキ分光反射率算出部１２１が求めた色における原色インキ各々の配合比を、見本の特色インキを構成する原色インキの種類と、その配合比として、一時記憶部１２５の特色インキ配合比テーブルに対し、特色インキの識別情報である特色インキ識別情報とともに書き込んで記憶させる。
【００５５】
図４は、特色インキを構成する原色インキの配合比を決定し、特色インキの吸収係数Ｋｔ（λ）、散乱係数Ｓｔ（λ）を算出するフローチャートである。
ステップＳ１０１：
ユーザは再現に用いるベタ印刷の測色値を測定し、色予測モデル生成部１２に対して入力する。
【００５６】
ステップＳ１０２：
特色インキ分光反射率算出部１２１は、印刷媒体の分光反射率である下地分光反射率Ｒ０（λ）を、測定分光反射率データベース１０５から読み出す。
【００５７】
ステップＳ１０３：
次に、特色インキ分光反射率算出部１２１は、特色インキを構成する原色インキの組み合わせを、近似色データベース１２８から読み出す。例えば、このとき、特色インキ分光反射率算出部１２１は、ベタ印刷から取得した測色値と最も近い測色値を持つ特色インキを、近似色データベース１２８から抽出して選択する。
【００５８】
ステップＳ１０４：
吸収係数・散乱係数算出部１０７は、ステップＳ１０３において選択された特色インキを構成する原色インキの散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）を、吸収係数・散乱係数データベース１０６から読み出す。
【００５９】
ステップＳ１０５：
特色インキ配合比決定部１２２は、特色インキを構成する原色インキ各々の配合比に初期値を設定する。初期値には、ステップＳ１０３で近似色データベース１２８から抽出した特色の配合比を設定する。
【００６０】
ステップＳ１０６：
次に、吸収係数・散乱係数算出部１０７は、（４）式に示すように、原色インキの各々の配合比を、それぞれ原色インキの散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）に乗じる。
そして、吸収係数・散乱係数算出部１０７は、乗算結果を加算することにより、その配合比（例えば、（４）式におけるα、β）における特色インキの散乱係数Ｓｔ（λ）及び吸収係数Ｋｔ（λ）を算出する。
【００６１】
特色インキ分光反射率算出部１２１は、（１）式に対し、印刷媒体の下地分光反射率Ｒ０（λ）と、散乱係数Ｓｔ（λ）及び吸収係数Ｋｔ（λ）とを代入して、特色インキの分光反射率を算出する。さらに、特色インキ再現色算出部１２９は、特色インキ分光反射率算出部１２１が算出した特色インキの分光反射率から、測色値を算出する。
特色インキ配合比決定部１２２は、ベタ印刷の測色値Ｌ*ａ*ｂ*値と、特色インキ再現色算出部１２９の求めた測色値Ｌ*ａ*ｂ*値との色差を求める。
【００６２】
ステップＳ１０７：
そして、特色インキ配合比決定部１２２は、上記色差が予め設定した許容範囲内であるか否かの判定を行う。
このとき、特色インキ配合比決定部１２２は、色差が予め設定した許容範囲内である場合、処理をステップＳ１０８へ進める。一方、特色インキ配合比決定部１２２は、色差が予め設定した許容範囲内にない場合、処理をステップＳ１１０へ進める。
【００６３】
ステップＳ１０８：
次に、特色インキ配合比決定部１２２は、ベタ印刷に対応する特色インキにおける原色インキの配合比を現在の配合比とする。
そして、特色インキ配合比決定部１２２は、特色インキにおける原色インキの種類と、原色インキ各々の配合比を、一時記憶部１２５の配合比テーブルに書き込んで記憶させる。
【００６４】
ステップＳ１０９：
次に、特色インキ配合比決定部１２２は、ベタ印刷に対応する特色インキの散乱係数Ｓｔ（λ）及び吸収係数Ｋｔ（λ）を、現在の散乱係数Ｓｔ（λ）及び吸収係数Ｋｔ（λ）とする。
そして、特色インキ配合比決定部１２２は、特色インキの散乱係数Ｓｔ（λ）及び吸収係数Ｋｔ（λ）を、一時記憶部１２５の散乱吸収係数テーブルに書き込んで記憶させる。
【００６５】
ステップＳ１１０：
特色インキ配合比決定部１２２は、配合比の変更回数が予め設定された規定以内か否かの判定を行う。
このとき、特色インキ配合比決定部１２２は、配合比の変更回数が予め設定された規定回数以内である場合、処理をステップＳ１１１へ進める。一方、特色インキ配合比決定部１２２は、配合比の変更回数が予め設定された規定回数を超えた場合、処理をステップＳ１１２へ進める。
【００６６】
ステップＳ１１１：
特色インキ配合比決定部１２２は、特色インキを構成する原色インキの各々の配合比を変更する。
そして、特色インキ配合比決定部１２２は、配合比の変更回数を計数するカウンタをインクリメント（１つカウント値を上げる）し、処理をステップＳ１０６へ進める。
【００６７】
ステップＳ１１２：
特色インキ配合比決定部１２２は、原色インキの組み合わせを変更するか否かの選択画面を、色予測モデル生成部１２の図示しない表示部に表示する。
特色インキ配合比決定部１２２は、ユーザが原色インキの組み合わせの変更を選択した場合、処理をステップＳ１１３へ進める。一方、特色インキ配合比決定部１２２は、ユーザが原色インキの組み合わせを変更しない選択をした場合、処理をステップＳ１０８へ進める。
【００６８】
ステップＳ１１３：
特色インキ配合比決定部１２２は、配合比の変更回数を計数するカウンタをリセットし、計数値、すなわち変更回数を「０」とする。
【００６９】
ステップＳ１１４：
特色インキ配合比決定部１２２は、 特色インキ分光反射率算出部１２１に対し、原色インキの組み合わせを変更する制御信号を出力する。
特色インキ分光反射率算出部１２１は、特色インキを構成する原色インキの組み合わせを、近似色データベース１２８から新たに読み出す。
【００７０】
図２に戻り、特色インキ濃度階調出現率算出部１２３は、原色インキを所定の比率で混ぜ合わせて生成した特色インキの濃度階調領域毎の出現率関数ａｍ（ｓ）を求める。そして、特色インキ濃度階調出現率算出部１２３は、求めた出現率関数ａｍ（ｓ）を一時記憶部１２５の特色インキ濃度階調出現率テーブルに対して書き込んで記憶させる。ここで、特色インキ濃度階調出現率算出部１２３は、特色インキを構成する原色インキのいずれかの出現率関数ａｍ（ｓ）を、濃度階調出現率テーブルデータベース１０４から読み出し、この特色インキの出現率関数ａｍ（ｓ）とする。また、特色インキ濃度階調出現率算出部１２３は、特色インキを構成する原色インキの各々の出現率関数ａｍ（ｓ）を、それぞれの原色インキの配合比により組合わせ、特色インキの出現率関数ａｍ（ｓ）としても良い。
【００７１】
・分光反射率予測部１０８の動作
次に、分光反射率予測部１０８は、拡張ノイゲバウア原色算出部１０８１、拡張ノイゲバウア原色出現率算出部１０８２及び分光反射率算出部１０８３を備えている。
拡張ノイゲバウア原色算出部１０８１は、重ね合わせる色の順番に従い、下地になるインキ（原色インキあるいは特色インキ）と、下地の表面に印刷されるインキとを決定する。また、拡張ノイゲバウア原色算出部１０８１は、重ね合わせるインキの各々の指令網点面積率を、入力部１０１から読み込む。
【００７２】
拡張ノイゲバウア原色算出部１０８１は、下地のインキに対して重なる原色インキの散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）を、吸収係数・散乱係数データベース１０６から読み込む。また、拡張ノイゲバウア原色算出部１０８１は、下地のインキに対して重なる特色インキの散乱係数Ｓｔ（λ）及び吸収係数Ｋｔ（λ）を、一時記憶部１２５の散乱吸収係数テーブルから読み込む。
【００７３】
次に、拡張ノイゲバウア原色算出部１０８１は、（１）式に下地となるインキ（特色インキまたは原色インキ）の網点における濃度階調領域の濃度階調分光反射率Ｒｉｍ（λ）を下地分光反射率Ｒ０（λ）として代入し、下地のインキに対して重なるインキの散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）と、濃度階調領域の膜厚係数Ｘｍとの各々を代入することにより、下地のインキの上に重なって印刷されるインキの網点を構成する濃度階調領域における濃度階調分光反射率Ｒｉｍ（λ）を算出する。
【００７４】
ここで、拡張ノイゲバウア原色算出部１０８１は、上述した下地のインキの網点の濃度階調領域と、この下地のインキの上に印刷されるインキの網点を構成する濃度階調領域の各々の重なり部分の濃度階調分光反射率Ｒｉｍ（λ）の計算を、下地のインキの網点における濃度階調領域と、下地のインキの網点上に印刷されるインキ（特色インキまたは原色インキ）の網点の濃度階調領域との組合わせの全てにおいて行う。
【００７５】
図５は、下地としてのインキの上に重ねて印刷されるインキの濃度階調分光反射率Ｒｉｍ（λ）の算出を説明する図である。
濃度階調分光反射率算出部１０２は、印刷媒体の分光反射率を下地の下地分光反射率Ｒ０（λ）とし、インキ１０００の領域の吸収係数Ｋ（λ）と散乱係数Ｓ（λ）とを用いて、（１）式から印刷媒体に対して印刷されたインキ１０００（特色インキまたは原色インキ）の領域の網点の濃度階調領域の濃度階調分光反射率Ｒｉｍ（λ）を求める。
【００７６】
そして、拡張ノイゲバウア原色算出部１０８１は、領域１００２に示す分光反射率のインキ１０００の上に対し、インキの領域１００４の上部に印刷された領域１００５に示す分光反射率ＲＫＭ（λ）を算出する。ここで、拡張ノイゲバウア原色算出部１０８１は、インキ１０００の領域における網点の濃度階調領域の分光反射率ＲＫＭ（λ）を下地の下地分光反射率Ｒ０（λ）とし、領域１００３に示すインキ１００１の吸収係数Ｋ（λ）と散乱係数Ｓ（λ）と、濃度階調領域の膜厚係数Ｘｍを、（１）式に代入して、領域１００４におけるインキ１０００の領域に重ねて印刷されたインキ１００１（特色インキまたは原色インキ）の領域１００４の網点における濃度階調領域の分光反射率ＲＫＭ（λ）を求める。
【００７７】
これにより、拡張ノイゲバウア原色算出部１０８１は、下地のインキ（特色インキまたは原色インキ）の網点における濃度階調領域と、下地のインキの網点上に印刷されるインキの網点の濃度階調領域との組合わせにおいて、下地のインキの網点上に印刷される上部のインキ（特色インキまたは原色インキ）の網点の濃度階調領域の重なり領域における分光反射率ＲＫＭ（λ）を、後述するように、下地のインキの網点と上部のインキ網点とにおける濃度階調領域との組合わせの全てについて算出することにより、インキが重ねて印刷された印刷部分における分光反射率ＲＫＭ（λ）を求める。拡張ノイゲバウア原色算出部１０８１は、求めた分光反射率ＲＫＭ（λ）を一時記憶部１２５のノイゲバウア原色テーブルに書き込んで記憶させる。
【００７８】
図６は、特色インキの重ね合わせによりノイゲバウア原色を算出する処理を示すフローチャートである。また、重ね合わせるインキとしては特色インキだけでなく、原色インキと特色インキとを組み合わせても良い。
ステップＳ２０１：
ユーザは、特色インキの組み合わせに対応し、それぞれの組み合わせにおける特色インキの刷り順を入力し、色予測システムに対して設定する。
【００７９】
ステップＳ２０２：
拡張ノイゲバウア原色算出部１０８１は、印刷媒体の分光反射率である下地分光反射率Ｒ０（λ）を、測定分光反射率データベース１０５から読み出す。
【００８０】
ステップＳ２０３：
拡張ノイゲバウア原色算出部１０８１は、重ね合わせる特色インキの散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）の各々を、一時記憶部１２５の散乱吸収係数テーブルから読み込む。
【００８１】
ステップＳ２０４：
拡張ノイゲバウア原色算出部１０８１は、ユーザが入力する膜厚補正値である膜厚係数を分光反射率に反映させるため、すなわち（１）式において用いる膜厚係数を内部の記憶部に記憶させる。
【００８２】
ステップＳ２０５：
拡張ノイゲバウア原色算出部１０８１は、特色インキあるいは原色インキを下地として印刷した印刷物の分光反射率ＲＫＭ（λ）を算出する。このとき、拡張ノイゲバウア原色算出部１０８１は、重ね合わせる特色インキにおけるそれぞれの膜厚係数の組み合わせの全てのノイゲバウア原色を算出する。すなわち、拡張ノイゲバウア原色算出部１０８１は、下地の印刷媒体に印刷した特色インキの分光反射率ＲＫＭ（λ）を求め、この求めた分光反射率ＲＫＭ（λ）を新たな下地の分光反射率とする。そして、拡張ノイゲバウア原色算出部１０８１は、求めた分光反射率ＲＫＭ（λ）の特色インキに対して、新たな特色インキを重ね合わせた際の分光反射率を算出し、これをノイゲバウア原色の分光反射率ＲＫＭ（λ）とする。
ここで、下地がプロセスインキのベタなど、予め分光反射率の実測値が既知である場合、媒体に印刷した際の分光反射率を算出するのではなく、実測値である分光反射率を、重ね合わせる下地の分光反射率として用いても良い。
そして、拡張ノイゲバウア原色算出部１０８１は、算出した分光反射率ＲＫＭ（λ）を、一時記憶部１２５のノイゲバウア原色テーブルに書き込んで記憶させる。
【００８３】
ステップＳ２０６：
拡張ノイゲバウア原色算出部１０８１は、特色インキの組み合わせの全てのノイゲバウア原色の算出が終了したか否かの判定を行う。
このとき、拡張ノイゲバウア原色算出部１０８１は、全てのノイゲバウア原色の算出が終了した場合、処理を終了する。一方、拡張ノイゲバウア原色算出部１０８１は、全てのノイゲバウア原色の算出が終了していない場合、処理をステップＳ２０７へ進める。
【００８４】
ステップＳ２０７：
拡張ノイゲバウア原色算出部１０８１は、特色インキの組み合わせを変更し、重ね合わせる特色インキの順番を、次の特色インキの組み合わせの順番とし、処理をステップＳ２０５に進める。
【００８５】
図２に戻り、拡張ノイゲバウア原色出現率算出部１０８２は、濃度階調出現率テーブルデータベース１０４から、出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）を読み出す。
そして、拡張ノイゲバウア原色出現率算出部１０８２は、下地のインキの指令網点面積率のおける網点で出現する濃度領域と、この下地のインキに重ねるインキの指令網点面積率における網点で出現する濃度改良領域との各々の出現確率を、濃度階調出現率テーブルデータベース１０４から読み出した出現確率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）から求める。
また、拡張ノイゲバウア原色出現率算出部１０８２は、下地のインキの濃度階調領域と、下地のインキに重ねて印刷されるインキの濃度階調領域との重なりの出現率の割合を算出する。
【００８６】
図７は、下地のインキの濃度階調領域と、下地のインキ（特色インキまたは原色インキ）に重ねて印刷されるインキの濃度階調領域との重なりの出現率の割合の算出を説明する図である。
図７においては、簡単のため重ね合わせるインキを２種類とし、濃度階調領域の種類も２種類としている。しかしながら、重ね合わせるインキと濃度階調領域の種類との各々が３以上の複数でも、下地のインキの濃度階調領域と、下地のインキに重ねて印刷されるインキの濃度階調領域との重なりの出現率の割合の算出は、以下の説明と同様に行うことができる。
【００８７】
図７（ａ）は、濃度階調領域が２種類ある場合を示しており、インキ＃１の指令網点面積率（設定網％）と、指令網点面積率における網点に出現する濃度階調領域であるコア及びフリンジの各々の出現率との対応関係を示している。また、インキ＃２もインキ＃１と同様の対応関係を有している。インキ＃１及びインキ＃２の各々は、原色インキあるいは特色インキのいずれかである。また、インキ＃１及びインキ＃２の各々において、濃度階調領域が濃度階調コア領域１及び濃度階調フリンジ領域２のそれぞれ２種類があり、濃度階調コア領域１は膜厚１００％であり、濃度階調フリンジ領域２は膜厚５０％である。また、他インキも同様の対応関係を有している。
【００８８】
図７（ｂ）は、濃度階調コア領域１及び濃度階調フリンジ領域２の各々の重なりを示している。図７（ｂ）において、例えばインキ＃１がシアン（Ｃ）であり、インキ＃２がマジェンダ（Ｍ）である。濃度階調コア領域１及び濃度階調フリンジ領域２の各々の重なりの組み合わせとしては、領域Ｑ１から領域Ｑ９の９種類がある。領域Ｑ１は、シアンの濃度階調コア領域１のみの領域である。領域Ｑ２は、シアンの濃度階調フリンジ領域２のみの領域である。領域Ｑ３は、マジェンダの濃度階調コア領域１のみの領域である。領域Ｑ４は、マジェンダの濃度階調フリンジ領域２のみの領域である。領域Ｑ５は、シアン及びマジェンダ各々の濃度階調コア領域１が重なった領域である。領域Ｑ６は、マジェンダの濃度階調フリンジ領域２とシアンの濃度階調コア領域１が重なった領域である。領域Ｑ７は、マジェンダの濃度階調コア領域１とシアンの濃度階調フリンジ領域２が重なった領域である。領域Ｑ８は、シアン及びマジェンダの各々の濃度階調フリンジ領域２が重なった領域である。領域Ｑ９は、シアン及びマジェンダのいずれのインキも存在しない領域である。
【００８９】
図８は、図７における領域Ｑ１から領域Ｑ９のそれぞれの出現率の算出結果を示すテーブルの図である。図８においては、原色インキＣ（シアン）及びＭ（マジェンダ）で説明しているが、特色インキでも同様に出現率が算出される。このテーブルにおいて、Ｃはシアンの濃度階調領域であり、Ｍはマジェンダの濃度階調領域であり、ＣＭはシアン及びマジェンダの濃度階調領域が重なった領域であり、Ｗはシアン及びマジェンダのいずれのインキも存在しない領域を示している。図８において、出現率α１はシアンの濃度階調コア領域１の出現率であり、出現率α２はシアンの濃度階調フリンジ領域２の出現率である。出現率α０は、出現率α１及び出現率α２の各々を加算したものである（α０＝α１＋α２）。出現率β１はマジェンダの濃度階調コア領域１の出現率であり、出現率β２はマジェンダの濃度階調フリンジ領域２の出現率である。出現率β０は、出現率β１及び出現率β２の各々を加算したものである（β０＝β１＋β２）。
【００９０】
領域Ｑ１及び領域Ｑ２の単次の出現率は、シアンのインキが出現する出現率α０に対し、マジェンダのインキが出現しない率である（１−β０）を乗算したα０＊（１−β０）となり、シアンのインキのみの領域の出現率を示す。本実施形態において、＊は乗算を示している。
領域Ｑ３及び領域Ｑ４の単次の出現率は、マジェンダのインキが出現する出現率β０に対し、シアンのインキが出現しない率である（１−α０）を乗算したβ０＊（１−α０）となり、マジェンダのインキのみの領域の出現率を示す。
【００９１】
領域Ｑ５から領域Ｑ８の単次の出現率は、シアンのインキが出現する出現率α０に対し、マジェンダのインキが出現する出現率β０を乗算したα０＊β０となり、シアンのインキとマジェンダのインキとが重なっている領域の出現率を示す。
領域Ｑ９の単次の出現率は、シアンのインキが出現しない率である（１−α０）に対し、マジェンダのインキが出現しない率である（１−β０）を乗算した（１−α０）＊（１−β０）となり、シアンのインキとマジェンダのインキとの双方が存在しない領域の出現率を示す。
【００９２】
領域Ｑ１の副次の出現率は、シアンのインキの網点の濃度階調コア領域１のみの出現率を示しており、濃度階調コア領域１の出現率α１を、濃度階調コア領域１の出現率α１及び濃度階調フリンジ領域２の出現率α２の加算結果により除算した率である。
領域Ｑ２の副次の出現率は、シアンのインキの網点の濃度階調フリンジ領域２のみの出現率を示しており、濃度階調コア領域１の出現率α１を、濃度階調コア領域１の出現率α１及び濃度階調フリンジ領域２の出現率α２の加算結果により除算し、この除算結果を１から減算した率である。
【００９３】
領域Ｑ３の副次の出現率は、マジェンダのインキの網点の濃度階調コア領域１のみの出現率を示しており、濃度階調コア領域１の出現率β１を、濃度階調コア領域１の出現率β１及び濃度階調フリンジ領域２の出現率β２の加算結果により除算した率である。
領域Ｑ４の副次の出現率は、マジェンダのインキの網点の濃度階調フリンジ領域２のみの出現率を示しており、濃度階調コア領域１の出現率β１を、濃度階調コア領域１の出現率β１及び濃度階調フリンジ領域２の出現率β２の加算結果により除算し、この除算結果を１から減算した率である。
【００９４】
領域Ｑ５の副次の出現率は、シアン及びマジェンダの各々のインキの網点の濃度階調コア領域１の重なり部分の出現率を示しており、シアンのインキの網点の濃度階調コア領域１の出現率α１を、濃度階調コア領域１の出現率α１及び濃度階調フリンジ領域２の出現率α２の加算結果により除算した率と、マジェンダのインキの網点の濃度階調コア領域１の出現率β１を、濃度階調コア領域１の出現率β１及び濃度階調フリンジ領域２の出現率β２の加算結果により除算した率とを乗算した結果の率である。
【００９５】
領域Ｑ６の副次の出現率は、シアンのインキの網点の濃度階調コア領域１とマジェンダのインキの網点の濃度階調フリンジ領域２の重なり部分の出現率を示しており、シアンのインキの網点の濃度階調コア領域１の出現率α１を、濃度階調コア領域１の出現率α１及び濃度階調フリンジ領域２の出現率α２の加算結果により除算した率と、マジェンダのインキの網点の濃度階調コア領域１の出現率β１を、濃度階調コア領域１の出現率β１及び濃度階調フリンジ領域２の出現率β２の加算結果により除算した率を１から減算した率とを乗算した結果の率である。
【００９６】
領域Ｑ７の副次の出現率は、シアンのインキの網点の濃度階調フリンジ領域２とマジェンダのインキの網点の濃度階調コア領域１の重なり部分の出現率を示しており、シアンのインキの網点の濃度階調コア領域１の出現率α１を、濃度階調コア領域１の出現率α１及び濃度階調フリンジ領域２の出現率α２の加算結果により除算した率を１から減算した率と、マジェンダのインキの網点の濃度階調コア領域１の出現率β１を、濃度階調コア領域１の出現率β１及び濃度階調フリンジ領域２の出現率β２の加算結果により除算した率とを乗算した結果の率である。
【００９７】
領域Ｑ８の副次の出現率は、シアンのインキの網点の濃度階調フリンジ領域２とマジェンダのインキの網点の濃度階調フリンジ領域２の重なり部分の出現率を示しており、シアンのインキの網点の濃度階調コア領域１の出現率α１を、濃度階調コア領域１の出現率α１及び濃度階調フリンジ領域２の出現率α２の加算結果により除算した率を１から減算した率と、マジェンダのインキの網点の濃度階調コア領域１の出現率β１を、濃度階調コア領域１の出現率β１及び濃度階調フリンジ領域２の出現率β２の加算結果により除算した率を１から減算した率とを乗算した結果の率である。
【００９８】
領域Ｑ９の副次の出現率は、シアンのインキとマジェンダのインキとの双方が存在しない領域の出現率を示し、「１」である。
【００９９】
上述したように、使用するインキの重ね合わせに用いる出現率の式は、予め設定され図８に示すテーブルとして、濃度階調出現率テーブルデータベース１０４に書き込んで記憶されている。
図２に戻り、拡張ノイゲバウア原色出現率算出部１０８２は、重ね合わせるインキの種類と、重ね合わせるインキの網点を示す指令網点面積率との各々の組み合わせにより、濃度階調出現率テーブルデータベース１０４の図８のテーブルの式を読み込む。
また、拡張ノイゲバウア原色出現率算出部１０８２は、濃度階調出現率テーブルデータベース１０４から、シアンとマジェンダとの各々の濃度階調コア領域１及び濃度階調フリンジ領域２それぞれの出現率α１、α２、β１、β２を読み込む。
【０１００】
拡張ノイゲバウア原色算出部１０８１は、シアンが下地で、マジェンダがシアンに重ね合わされる場合、シアンの分光反射率を下地の下地分光反射率Ｒ０（λ）として、指令網点面積率で出現するマジェンダの濃度階調領域の膜厚を用い、（１）式により、重なり部分における濃度階調領域の分光反射率ＲＫＭ（λ）を算出する。例えば、拡張ノイゲバウア原色算出部１０８１は、領域Ｑ５における分光反射率ＲＫＭ（λ）を求める際、印刷媒体上のシアンのインキの膜厚１００％における濃度階調分光反射率Ｒｉｍ（λ）を下地の下地分光反射率Ｒ０（λ）として、マジェンダのインキの膜厚１００％を重ねた際の重なり部分における印刷部分の分光反射率ＲＫＭ（λ）を算出する。
【０１０１】
同様に、拡張ノイゲバウア原色算出部１０８１は、領域Ｑ７における分光反射率ＲＫＭ（λ）を求める際、印刷媒体上のシアンのインキの膜厚５０％における濃度階調分光反射率Ｒｉｍ（λ）を下地の下地分光反射率Ｒ０（λ）とし、マジェンダのインキの膜厚１００％を重ねた際の分光反射率ＲＫＭ（λ）を算出する。同様に、拡張ノイゲバウア原色の全ての組み合わせについて分光反射率ＲＫＭ（λ）を算出する。
次に、拡張ノイゲバウア原色出現率算出部１０８２は、領域Ｑ１の拡張ノイゲバウア原色出現率を算出する場合、単次の出現率α０＊（１−β０）と、副次の出現率α１／（α１＋α２）を乗算することで求める。また、拡張ノイゲバウア原色出現率算出部１０８２は、領域Ｑ２の拡張ノイゲバウア原色出現率を算出する場合、単次の出現率α０＊（１−β０）と、副次の出現率｛１−α１／（α１＋α２）｝を乗算することで求める。領域Ｑ３及び領域Ｑ４の各々は、シアンがマジェンダに変更になっただけで、上述した領域Ｑ１及び領域Ｑ２と同様に求めることができる。
拡張ノイゲバウア原色出現率算出部１０８２は、領域Ｑ５の拡張ノイゲバウア原色出現率を算出する場合、単次の出現率α０＊β０と、副次の出現率｛α１／（α１＋α２）｝＊｛β１／（β１＋β２）｝を乗算することで求める。
拡張ノイゲバウア原色出現率算出部１０８２は、領域Ｑ６の拡張ノイゲバウア原色出現率を算出する場合、単次の出現率α０＊β０と、副次の出現率｛α１／（α１＋α２）｝＊｛１−β１／（β１＋β２）｝を乗算することで求める。
拡張ノイゲバウア原色出現率算出部１０８２は、領域Ｑ７の拡張ノイゲバウア原色出現率を算出する場合、単次の出現率α０＊β０と、副次の出現率｛１−α１／（α１＋α２）｝＊｛β１／（β１＋β２）｝を乗算することで求める。
拡張ノイゲバウア原色出現率算出部１０８２は、領域Ｑ８の拡張ノイゲバウア原色出現率を算出する場合、単次の出現率α０＊β０と、副次の出現率｛１−α１／（α１＋α２）｝＊｛１−β１／（β１＋β２）｝を乗算することで求める。
拡張ノイゲバウア原色出現率算出部１０８２は、領域Ｑ９の拡張ノイゲバウア原色出現率を算出する場合、単次の出現率（１−α０）＊（１−β０）と、副次の出現率１を乗算することで求める。
【０１０２】
そして、分光反射率算出部１０８３は、上述した領域Ｑ１から領域Ｑ９の拡張ノイゲバウア原色の各々の分光反射率に拡張ノイゲバウア原色出現率を乗算し、波長毎に加算する。
分光反射率算出部１０８３は、シアンのインキの網点が印刷された印刷媒体に対し、マジェンダのインキの網点を重ねて印刷した際における印刷部分の第１予測分光反射率ＲＤ１（λ）を算出する。
【０１０３】
・分光光学濃度による濃度階調出現率の算出
また、濃度階調分光光学濃度算出部２０２は、入力部１０１より供給される原色インキの種類に応じて、吸収係数・散乱係数データベース１０６から、原色インキがベタで印刷媒体に印刷された印刷部分から求めた着色層の吸収係数Ｋ（λ）及び散乱係数Ｓ（λ）を読み出す。また、濃度階調分光光学濃度算出部２０２は、測定分光反射率データベース１０５から、入力部１０１より供給される原色インキの種類及び印刷媒体の種類に応じて、指令網点面積率毎の測定分光反射率Ｒｓ（λ）と、印刷媒体の下地分光反射率Ｒ０（λ）との各々を読み出す。
【０１０４】
そして、濃度階調分光光学濃度算出部２０２は、読み出した原色インキの吸収係数Ｋ（λ）及び散乱係数Ｓ（λ）と、印刷媒体の下地分光反射率Ｒ０（λ）と、濃度階調の膜厚係数Ｘｍとの各々を、すでに説明した（１）式（クベルカ・ムンクの式）に代入し、濃度階調領域の分光反射率として、濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）を算出する。
また、濃度階調分光光学濃度算出部２０２は、求めた濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）の各々を、以下の（５）式により、濃度階調分光光学濃度ＯＤｉ１（λ）、ＯＤｉ２（λ）、ＯＤｉ３（λ）、…、ＯＤｉｍ（λ）にそれぞれ変換する。
【０１０５】
【数５】
2019016924.tif 000007
【０１０６】
濃度階調分光光学濃度算出部２０２においても、濃度階調分光反射率算出部１０２と同様に、（１）式（クベルカ・ムンクの式）における印刷されたインキの膜厚係数Ｘｍは、印刷媒体に原色インキがベタで印刷された印刷部分を基にし、印刷部分の濃度階調を示す数値として用いる。すなわち、膜厚係数Ｘｍは任意に設定されており、例えば、最もインキの膜厚が厚いベタの印刷部分を１００％として、上記膜厚係数を１とし、この１を濃度階調領域の膜厚の段階数ｍに対応させてｍ分割する。例えば、指令網点面積率の示す濃度階調領域の膜厚の段階が５段階であれば、ｍ＝１、２、３、４、５であり、膜厚係数Ｘｍは、それぞれの濃度階調領域の膜厚の段階に対応して、Ｘ１＝１．０、Ｘ２＝０．８、Ｘ３＝０．６、Ｘ４＝０．４、Ｘ５＝０．２とする。
【０１０７】
膜厚係数Ｘｍを、すでに述べたように、印刷媒体の下地分光反射率Ｒ０（λ）と、吸収係数Ｋ（λ）と、散乱係数Ｓ（λ）とともに、クベルカ・ムンクの式である（１）式に代入し、それぞれの指令網点面積率の網点が含む濃度階調領域の分光反射率として、濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）の各々を算出する。そして、濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）の各々を、（５）式により、濃度階調分光光学濃度ＯＤｉ１（λ）、ＯＤｉ２（λ）、ＯＤｉ３（λ）、…、ＯＤｉｍ（λ）にそれぞれ変換する。この濃度階調分光光学濃度ＯＤｉ１（λ）、ＯＤｉ２（λ）、ＯＤｉ３（λ）、…、ＯＤｉｍ（λ）が、後述する計算モデルで用いる、網点を構成する複数の濃度階調領域の各々の分光光学濃度として用いられる。
【０１０８】
濃度階調出現率算出部２０３は、濃度階調分光光学濃度算出部２０２から、濃度階調分光光学濃度ＯＤｉ１（λ）、ＯＤｉ２（λ）、ＯＤｉ３（λ）、…、ＯＤｉｍ（λ）の各々を読み込む。また、濃度階調出現率算出部２０３は、指令網点面積率毎の測定分光反射率Ｒｓ（λ）の各々を、測定分光反射率データベース１０５から読み込む。
そして、濃度階調出現率算出部２０３は、（５）式を用いることにより、測定分光反射率Ｒｓ（λ）の各々を、測定分光光学濃度ＯＤｓ（λ）にそれぞれ変換する。
次に、濃度階調出現率算出部２０３は、以下の（６）式（計算モデル）に対して、濃度階調分光光学濃度ＯＤｉ１（λ）、ＯＤｉ２（λ）、ＯＤｉ３（λ）、…、ＯＤｉｍ（λ）の各々を代入し、後述する処理により算出分光光学濃度ＯＤ’（ｓ，λ）を算出する。
【０１０９】
【数６】
2019016924.tif 000008
【０１１０】
ここで、濃度階調出現率算出部２０３は、（６）式において、出現率（用紙の印刷部分における網点を構成する各濃度階調の濃度階調領域の面積の比率）の数値を変更しつつ、算出分光光学濃度ＯＤ’（ｓ，λ）を算出する。そして、濃度階調出現率算出部２０３は、以下の（７）式において算出された算出分光光学濃度ＯＤ’（ｓ，λ）の各々と、測定分光光学濃度ＯＤｓ（λ）との平均二乗誤差ＲＭＳＥを、指令網点面積率毎に所定の波長範囲において求める。濃度階調出現率算出部２０３は、算出分光光学濃度ＯＤ’（ｓ，λ）と測定分光光学濃度ＯＤｓ（λ）との平均二乗誤差が最も小さくなる濃度階調領域の各々の出現率を求める。ここで、ｓは指令網点面積率である。
【０１１１】
【数７】
2019016924.tif 000009
【０１１２】
上記（７）式において、濃度階調出現率算出部２０３は、波長λが３８０ｎｍから７３０ｎｍをｎ分割した刻み幅により、算出分光光学濃度ＯＤ’（ｓ，λ）と測定分光光学濃度ＯＤｓ（λ）とにおける各波長λに誤差の二乗を加算した平均二乗誤差ＲＭＳＥを、指令網点面積率毎に求めている。
【０１１３】
そして、濃度階調出現率算出部２０３は、濃度階調領域の各々の出現率により、各濃度階調領域の出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）を求める。ここで、濃度階調出現率算出部２０３は、濃度階調領域毎に、得られた出現率を用いて、この出現率を指令網点面積率ｓの２次関数などをフィッティングし、出現率関数を求めてもよい。濃度階調出現率算出部２０３は、求めた各濃度階調領域の出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）を、濃度階調出現率テーブルデータベース２０４に書き込んで記憶させる。出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）については、分光反射率予測部１０８における説明と同様に、（６）式の計算モデルで用いる濃度階調領域の出現率を求める関数であり、図３の説明におけるように、グラビア印刷による網点の構造をモデル化するために用いている。
【０１１４】
上述したように、濃度階調出現率算出部２０３は、（６）式の計算モデルを用いて、指令網点面積率によるインキ（原色インキ）が印刷された印刷部分において、網点が含む濃度階調毎に、濃度階調の濃度階調分光光学濃度ＯＤｉｍ（λ）と濃度階調の出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）とを乗算した結果を加算して、指令網点面積率の印刷部分の分光光学濃度を算出する。
【０１１５】
また、濃度階調出現率算出部２０３は、求めた濃度階調領域の指令網点面積率における各濃度階調領域の出現率を示す出現率関数ａｍ（ｓ）を、濃度階調出現率テーブルデータベース２０４に対し、濃度階調領域の濃度階調分光光学濃度と対応させて書き込んで記憶させる。同様に、濃度階調出現率算出部２０３は、他の原色インキの各々に対しても、濃度階調領域の指令網点面積率における各濃度階調領域の出現率を示す出現率関数ａｍ（ｓ）を、濃度階調出現率テーブルデータベース２０４に対して、書き込んで記憶させる。
【０１１６】
また、原色インキを所定の比率で混ぜ合わせて生成した特色インキの場合、上述した原色インキの濃度階調領域の出現率関数ａｍ（ｓ）を濃度階調出現率テーブルデータベース２０４から読み出して用いる。ここで、特色インキで混ぜ合わせる原色インキのいずれかの出現率関数ａｍ（ｓ）を用いても良いし、また混ぜ合わせる原色インキそれぞれの出現率関数ａｍ（ｓ）を配合率に応じて組合わせて用いても良い。
この特色インキの場合、分光反射率予測部１０８における説明と同様に、混ぜ合わせる原色インキの比率に応じて、すでに説明した（４）式において、特色インキの散乱係数Ｓｔ（λ）及び吸収係数Ｋｔ（λ）を算出する。
【０１１７】
上記（４）式において、係数α及び係数βの各々は、原色インキ＃１及び原色インキ＃２それぞれの混ぜ合わせる比率を示す係数である。原色インキ＃１の吸収係数Ｋ１（λ）に対して係数αを乗算し、原色インキ＃２の吸収係数Ｋ２に対して係数βを乗算し、加算した数値を、特色インキの吸収係数Ｋｔ（λ）としている。同様に、原色インキ＃１の散乱係数Ｓ１（λ）に対して係数αを乗算し、原色インキ＃２の散乱係数Ｓ２（λ）に対して係数βを乗算し、加算した数値を、特色インキの散乱係数Ｓｔ（λ）としている。
濃度階調分光光学濃度算出部２０２は、吸収係数・散乱係数データベース１０６から、特色インキを構成する原色のインキ＃１及び原色のインキ＃２の各々の散乱係数Ｓ１（λ）及びＳ２（λ）と、吸収係数Ｋ１（λ）及びＫ２（λ）との各々を読み出す。また、濃度階調分光光学濃度算出部１０２Ａは、測定分光反射率データベース１０５から、印刷媒体の下地分光反射率Ｒ０（λ）を読み出す。
次に、濃度階調分光光学濃度算出部２０２は、算出した特色インキの散乱係数Ｓｔ（λ）及び吸収係数Ｋｔ（λ）と、印刷媒体の下地分光反射率Ｒ０（λ）と、濃度階調の膜厚係数Ｘｍとの各々を、上記（１）式に代入して、濃度階調領域の分光反射率として、濃度階調分光反射率Ｒｉ１（λ）、Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）を算出する。
濃度階調分光光学濃度算出部２０２は、算出した濃度階調分光反射率Ｒｉ１（λ）、
Ｒｉ２（λ）、Ｒｉ３（λ）、…、Ｒｉｍ（λ）の各々を、（５）式を用い、濃度階調領
域毎の濃度階調光学濃度ＯＤｉ１（λ）、ＯＤｉ２（λ）、ＯＤｉ３（λ）、ノ、ＯＤｉ
ｍ（λ）にそれぞれ変換する。
【０１１８】
・分光光学濃度予測部１１１の動作
次に、分光光学濃度予測部１１１は、拡張ノイゲバウア原色算出部１１１１、拡張ノイゲバウア原色出現率算出部１１１２、分光光学濃度算出部１１１３及び分光反射率算出部１１１４を備えている。
拡張ノイゲバウア原色算出部１１１１は、重ね合わせる色の順番に従い、下地になるインキ（原色インキあるいは特色インキ）と、下地の表面に印刷されるインキとを決定する。また、拡張ノイゲバウア原色算出部１１１１は、重ね合わせるインキの各々の指令網点面積率を、入力部１０１から読み込む。
【０１１９】
また、拡張ノイゲバウア原色算出部１１１１は、下地のインキの上に印刷されるインキ（原色インキあるいは特色インキ）の散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）を読み込む。
次に、拡張ノイゲバウア原色算出部１１１１は、すでに図５で説明されたように、（１）式に下地となるインキの濃度階調分光反射率Ｒｉｍ（λ）と、下地のインキの上に印刷されるインキの散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）と、濃度階調領域の膜厚係数Ｘｍとの各々を代入し、下地のインキの上に印刷される網点の濃度階調分光反射率Ｒｉｍ（λ）を算出する。
【０１２０】
ここで、拡張ノイゲバウア原色算出部１１１１は、上述した下地のインキの網点上に印刷されるインキの網点の濃度階調分光反射率Ｒｉｍ（λ）の計算を、下地のインキの網点における濃度階調領域と、下地のインキの網点上に印刷されるインキの網点の濃度階調領域との重ね合う部分の組合わせの全てにおいて行う。そして、拡張ノイゲバウア原色算出部１１１１は、得られた下地のインキの網点上に印刷されるインキの濃度階調分光反射率Ｒｉｍ（λ）を、（５）式により、濃度階調分光光学濃度ＯＤｉｍ（λ）に変換する。拡張ノイゲバウア原色算出部１１１１は、求めた濃度階調分光光学濃度ＯＤｉｍ（λ）を一時記憶部１２５のノイゲバウア原色テーブルに書き込んで記憶させる。
【０１２１】
図２に戻り、拡張ノイゲバウア原色出現率算出部１１１２は、図７の説明と同様に、下地のインキの濃度階調領域と、下地のインキに重ねて印刷されるインキの濃度階調領域との重なりの出現率の割合を算出する。
【０１２２】
また、使用するインキ（原色インキあるいは特色インキ）の重ね合わせに用いる出現率の式は、分光反射率予測部１０８において説明したように、図８に示す出現率のテーブルとして濃度階調出現率テーブルデータベース２０４に書き込んで記憶されている。
拡張ノイゲバウア原色出現率算出部１１１２は、重ね合わせるインキの種類と、重ね合わせるインキの網点を示す指令網点面積率との各々の組み合わせにより、濃度階調出現率テーブルデータベース２０４の図８のテーブルの式を読み込む。
また、拡張ノイゲバウア原色出現率算出部１１１２は、濃度階調出現率テーブルデータベース２０４から、シアンとマジェンダとの各々の濃度階調コア領域１及び濃度階調フリンジ領域２それぞれの出現率α１、α２、β１、β２を読み込む。
【０１２３】
また、拡張ノイゲバウア原色算出部１１１１は、シアンが下地で、マジェンダがシアンに重ね合わされる場合、シアンの分光反射率を下地の下地分光反射率Ｒ０（λ）として、指令網点面積率で出現するマジェンダの濃度階調領域の膜厚を用い、（１）式により、重なり部分における濃度階調領域の分光反射率ＲＫＭ（λ）を算出する。例えば、拡張ノイゲバウア原色算出部１１１１は、領域Ｑ５における分光反射率ＲＫＭ（λ）を求める際、印刷媒体上のシアンのインキの膜厚１００％における濃度階調分光反射率Ｒｉｍ（λ）を下地の下地分光反射率Ｒ０（λ）として、マジェンダのインキの膜厚１００％を下地に重ねた際の重なり部分における印刷部分の分光反射率ＲＫＭ（λ）を算出する。そして、拡張ノイゲバウア原色算出部１１１１は、算出した分光反射率ＲＫＭ（λ）を、（５）式により、分光光学濃度ＯＤＫＭ（λ）に変換する。
【０１２４】
同様に、拡張ノイゲバウア原色算出部１１１１は、領域Ｑ７における分光反射率ＲＫＭ（λ）を求める際、印刷媒体上のシアンのインキの膜厚５０％における濃度階調分光反射率Ｒｉｍ（λ）を下地の下地分光反射率Ｒ０（λ）とし、マジェンダのインキの膜厚１００％を重ねた際の分光反射率ＲＫＭ（λ）を算出する。そして、拡張ノイゲバウア原色算出部１１１１は、算出した分光反射率ＲＫＭ（λ）を、（５）式により、分光光学濃度ＯＤＫＭ（λ）に変換する。同様に、拡張ノイゲバウア原色の全ての組み合わせについて分光光学濃度ＯＤＫＭ（λ）を算出する。
拡張ノイゲバウア原色出現率算出部１１１２は、領域Ｑ１の拡張ノイゲバウア原色出現率を算出する場合、単次の出現率α０＊（１−β０）と、副次の出現率α１／（α１＋α２）を乗算することで求める。また、拡張ノイゲバウア原色出現率算出部１１１２は、領域Ｑ２の拡張ノイゲバウア原色出現率を算出する場合、単次の出現率α０＊（１−β０）と、副次の出現率｛１−α１／（α１＋α２）｝を乗算することで求める。領域Ｑ３及び領域Ｑ４の各々は、シアンがマジェンダに変更になっただけで、上述した領域Ｑ１及び領域Ｑ２と同様に求めることができる。
拡張ノイゲバウア原色出現率算出部１１１２は、領域Ｑ５の拡張ノイゲバウア原色出現率を算出する場合、単次の出現率α０＊β０と、副次の出現率｛α１／（α１＋α２）｝＊｛β１／（β１＋β２）｝を乗算することで求める。
拡張ノイゲバウア原色出現率算出部１１１２は、領域Ｑ６の拡張ノイゲバウア原色出現率を算出する場合、単次の出現率α０＊β０と、副次の出現率｛α１／（α１＋α２）｝＊｛１−β１／（β１＋β２）｝を乗算することで求める。
拡張ノイゲバウア原色出現率算出部１１１２は、領域Ｑ７の拡張ノイゲバウア原色出現率を算出する場合、単次の出現率α０＊β０と、副次の出現率｛１−α１／（α１＋α２）｝＊｛β１／（β１＋β２）｝を乗算することで求める。
拡張ノイゲバウア原色出現率算出部１１１２は、領域Ｑ８の拡張ノイゲバウア原色出現率を算出する場合、単次の出現率α０＊β０と、副次の出現率｛１−α１／（α１＋α２）｝＊｛１−β１／（β１＋β２）｝を乗算することで求める。
拡張ノイゲバウア原色出現率算出部１１１２は、領域Ｑ９の拡張ノイゲバウア原色出現率を算出する場合、単次の出現率（１−α０）＊（１−β０）と、副次の出現率１を乗算することで求める。
そして、分光光学濃度算出部１１１３は、上述した領域Ｑ１から領域Ｑ９の拡張ノイゲバウア原色の各々の光学濃度に拡張ノイゲバウア原色出現率を乗算し加算する。
分光光学濃度算出部１１１３は、シアンのインキの網点が印刷された印刷媒体に対し、マジェンダのインキの網点を重ねて印刷した際における印刷部分の予測光学濃度ＯＤ（λ）を算出する。
【０１２５】
そして、分光光学濃度算出部１１１３は、上述した領域Ｑ１から領域Ｑ９の拡張ノイゲバウア原色の各々の分光光学濃度に拡張ノイゲバウア原色出現率を乗算し、波長毎に加算する。
分光光学濃度算出部１１１３は、シアンのインキの網点が印刷された印刷媒体に対し、マジェンダのインキの網点を重ねて印刷した際における印刷部分の予測分光光学濃度ＯＤＤ（λ）を算出する。
【０１２６】
分光反射率算出部１１１４は、分光光学濃度算出部１１１３が算出した予測分光光学濃度ＯＤＤ（λ）を、以下の（８）式を用い、第２予測分光反射率ＲＤ２（λ）に変換する。
【０１２７】
【数８】
2019016924.tif 000010
【０１２８】
・混合予測部１１４の動作
混合予測部１１４は、予測パラメータデータベース１２６から、第１予測分光反射率ＲＤ１（λ）及び第２予測分光反射率ＲＤ２（λ）に乗算する重み係数ｗを読み出す。例えば、混合予測部１１４は、ｗＲＤ１（λ）と（１−ｗ）ＲＤ２（λ）とを加算して、統合予測分光反射率ＲＤ（λ）を求める。そして、混合予測部１１４は、統合予測分光反射率ＲＤ（λ）を予測分光反射率Ｒ（λ）としてカラープロファイル生成部１３に対して出力する。
【０１２９】
＜カラープロファイルの生成＞
そして、カラープロファイル生成部１３は、この予測分光反射率Ｒ（λ）に基づき、観測光源を定めて、三刺激値ＸＹＺやＣＩＥＬＡＢ値等を算出し、再現色の予測を行う色予測テーブルを作成する。つまり、このカラープロファイル生成部１３は、色予測処理の対象である入力データ（本実施例ではＣＭＹ値）の情報と、予測分光反射率Ｒ（λ）に基づき予測される再現色を表現する再現色情報とを対応付けるカラープロファイルとしての色予測テーブルを作成する。なお、本実施例においては、カラープロファイル生成部１３がこの色予測テーブルを公知のＩＣＣ（International Color Consortium）プロファイルフォーマットで作成し、記憶部１８に書き込んで記憶させる。
【０１３０】
図９は、本実施形態における色予測モデル生成部１２及びカラープロファイル生成部１３によるカラープロファイルの生成処理の動作例を説明するフローチャートである。
【０１３１】
ステップＳ４０１：
拡張ノイゲバウア原色出現率算出部１０８２は、一時記憶部１２５のノイゲバウア原色テーブルから、重なり合う領域の特色インキのノイゲバウア原色の分光反射率ＲＫＭ（λ）を読み込む。
同様に、拡張ノイゲバウア原色出現率算出部１１１２は、一時記憶部１２５のノイゲバウア原色テーブルから、重なり合う領域の特色インキのノイゲバウア原色の分光光学濃度ＯＤＫＭ（λ）を読み込む。
【０１３２】
ステップＳ４０２：
混合予測部１１４は、分光反射率予測部１０８が予測した第１予測分光反射率ＲＤ１（λ）と、分光光学濃度予測部１１１が予測した第２予測分光反射率ＲＤ２（λ）を統合する際、第１予測分光反射率ＲＤ１（λ）及び第２予測分光反射率ＲＤ２（λ）の各々に乗算する重み係数ｗを、予測パラメータデータベース１２６から読み出す。
【０１３３】
ステップＳ４０３：
次に、拡張ノイゲバウア原色出現率算出部１０８２は、一時記憶部１２５における色予測テーブルにおける重ね合わせるインキの指令網点面積率のマトリクスにおいて、未計算の指令網点面積率を抽出する。
そして、拡張ノイゲバウア原色出現率算出部１０８２は、抽出した指令網点面積率の組み合わせを、計算対象として設定する。上述した処理は、拡張ノイゲバウア原色出現率算出部１１１２が行う構成としても良い。
【０１３４】
ステップＳ４０４：
分光反射率予測部１０８は、上記指令網点面積率の組み合わせにより、第１予測分光反射率ＲＤ１（λ）を算出する。
分光光学濃度予測部１１１は、上記指令網点面積率の組み合わせにより、第２予測分光反射率ＲＤ２（λ）を算出する。
【０１３５】
また、混合予測部１１４は、第１予測分光反射率ＲＤ１（λ）及び第２予測分光反射率ＲＤ２（λ）にステップＳ４０２で読み出した重み係数ｗと（１−ｗ）を積算し、ｗＲＤ１（λ）と（１−ｗ）ＲＤ２（λ）とを加算して、統合予測分光反射率ＲＤ（λ）を求める。さらに、算出した統合予測分光反射率ＲＤ（λ）に対して、観測光源の分光分布と標準観測者を設定することにより、測色値を算出する。
【０１３６】
ステップＳ４０５：
そして、混合予測部１１４は、算出した測色値を、重ね合わせるインキ（特色インキあるいは原色インキ）の指令網点面積率の組み合わせに対応させ、記憶部１８の色予測テーブルに書き込んで記憶させる。
【０１３７】
ステップＳ４０６：
拡張ノイゲバウア原色出現率算出部１０８２は、記憶部１８の色予測テーブルにおける全ての指令網点面積率の組み合わせの統合予測分光反射率ＲＤ（λ）の算出が終了したか否かの判定を行う。
そして、拡張ノイゲバウア原色出現率算出部１０８２は、全ての指令網点面積率の組み合わせの統合予測分光反射率ＲＤ（λ）の算出が終了した場合、処理をステップＳ４０７へ進める。一方、拡張ノイゲバウア原色出現率算出部１０８２は、全ての指令網点面積率の組み合わせの統合予測分光反射率ＲＤ（λ）の算出が終了していない場合、処理をステップＳ４０３へ進める。上述した処理は、拡張ノイゲバウア原色出現率算出部１１１２が行う構成としても良い。
【０１３８】
ステップＳ４０７：
カラープロファイル生成部１３は、記憶部１８の色予測テーブルのデータを読み出し、公知のＩＣＣプロファイルフォーマットに変換し、変換後の色予測プロファイルのデータを記憶部１８に書き込んで記憶させる。
【０１３９】
上述したように、本実施形態によれば、原色インキから特色インキの吸収特性及び散乱特性を推定し、特色インキの分光反射率を算出し、特色インキの重ね刷りされた際の重なり領域の分光反射率を求めることができる。
このため、本実施形態によれば、面積変調階調表現と、濃度変調階調表現との双方の色表現を有するグラビア印刷などの印刷において特色インキを重ね刷りした印刷物における再現色の色予測を高い精度で容易に行うことができる。
【０１４０】
また、本実施形態によれば、拡張ノイゲバウア原色の分光反射率と、その出現率から算出した第１予測分光反射率ＲＤ１（λ）と、拡張ノイゲバウア原色の分光光学濃度と、その出現率から算出した第２予測分光反射率ＲＤ２（λ）とを、予め求めた重み係数ｗにより混合するため、より実測値の分光反射率に近い予測分光反射率として混合予測分光反射率ＲＷ（λ）を、予測分光反射率Ｒ（λ）として求めることができる。
【０１４１】
＜印刷色調整システムの動作説明＞
図１０は、本実施形態における印刷色調整システムによる印刷機の印刷物の色に対し、印刷デバイス（例えば、インクジェットプリンタ）の印刷する色を同様とする調整処理の動作例を説明するフローチャートである。
ステップＳ５０１：
ユーザは、印刷機が印刷物を印刷する際に用いたインキのベタ印刷のカラーパッチの色の測定を行い、それぞれのインキの測色値を得る。そして、ユーザは、入力手段（キーボード、タッチパネルなど）から、測定したインキの各々の測色値を、印刷色調整システム１に対して入力する。
また、ユーザは、上記印刷物の版の各々の画像データ（各画素の指令網点面積率のデータ）を、外部装置から印刷色調整システム１に対して入力する。
【０１４２】
ステップＳ５０２：
入力制御部１１は、入力手段から供給された測色値を、色予測モデル生成部１２に対して出力する。また、入力制御部１１は、入力手段から供給された印刷物の画像データを記憶部１８に対して書き込んで記憶させる。
そして、色予測モデル生成部１２は、すでに説明したように、色予測モデル（色予測テーブル）を生成し、生成した色予測モデルを記憶部１８に書き込んで記憶させる。
【０１４３】
ステップＳ５０３：
カラープロファイル生成部１３は、記憶部１８の色予測モデル（色予測テーブル）を参照し、インキの各々の指令網点面積率の組合わせに対応した予測分光反射率Ｒ（λ）の各々を、公知のＩＣＣプロファイルフォーマットに変換し、変換後のカラープロファイル（色予測プロファイル）のデータを記憶部１８に書き込んで記憶させる。
【０１４４】
ステップＳ５０４：
調整制御部１４は、画像データにおける画像から、印刷デバイスの印刷する印刷物の色を評価する参照画像とする部分画像を抽出するための画像を、表示部１６の表示画面に表示する。この参照画像は、特に画像データの画像において表現したい色が含まれる部分の部分画像である。
そして、調整制御部１４は、表示部１６の表示画面に表示された画像データの画像において、ユーザが選択した部分画像の領域を参照画像として画像データから切り出す。
【０１４５】
ステップＳ５０５：
制御値生成部１５は、調整制御部１４が印刷物の画像データから切り出した参照画像の画素の各々におけるインキそれぞれの指令網点面積率の組合わせを、記憶部１８のカラープロファイルを参照して予想測色値に変換する。
また、制御値生成部１５は、記憶部１８の印刷デバイスのデバイスカラープロファイルを参照し、参照画像の画素の各々の予想測色値を印刷デバイスの制御値に変換し、記憶部１８に書き込んで記憶させる。
そして、印刷データ出力部１７は、記憶部１８から参照画像の各画素の制御値を読み出し、読み出した参照画像の印刷物を印刷するための制御値を、印刷デバイスに対して出力する。
【０１４６】
ステップＳ５０６：
これにより、ユーザは、印刷デバイスの印刷した参照画像の印刷物と、印刷機の印刷物における参照画像に対応する領域との各々の色を観察して比較する。そして、ユーザは、印刷デバイスの印刷した参照画像の印刷物の色が、印刷機の印刷物における参照画像に対応する領域の色に対して、自身の許容範囲の違いに収まっている（すなわち、一致している）か否かの判定を行う。
【０１４７】
ステップＳ５０７：
ユーザは、印刷デバイスの印刷した参照画像の印刷物の色と、印刷機の印刷物における参照画像に対応する領域の色とが、一致していると判定した場合、処理をステップＳ５０８へ進める。
一方、印刷デバイスの印刷した参照画像の印刷物の色と、印刷機の印刷物における参照画像に対応する領域の色とが、一致していないと判定した場合、処理をステップＳ５０９へ進める。
【０１４８】
ステップＳ５０８：
ユーザは、所定の枚数を印刷するため、印刷色調整システム１に対して、入力手段を介して印刷枚数を入力する。
これにより、制御値生成部１５は、印刷機の印刷物の画像データにおける全ての画素の各々におけるインキそれぞれの指令網点面積率の組合わせを、記憶部１８のカラープロファイルを参照して予想測色値に変換し、記憶部１８に記憶させるとともに、印刷データ出力部１７に対して出力する。
そして、印刷データ出力部１７は、制御値生成部１５から供給される画像データの制御値を、印刷デバイスに対して出力し、印刷枚数の数の印刷物を印刷する。
【０１４９】
ステップＳ５０９：
ユーザは、印刷デバイスの参照画像の印刷物と、印刷機の印刷物との色の違いなどが許容範囲を超えている場合、色予測モデルにおけるパラメータの変更処理を行う。この変更処理については、後に詳述する。
そして、調整制御部１４は、ユーザが入力手段から入力する調整対象のパラメータの種類と、調整割合とに基づき、調整対象のパラメータの各々に対して多段階の調整パラメータを生成し、色予測モデル生成部１２に対して出力し、新たな色予測モデルの生成を行わせる。
【０１５０】
これにより、色予測モデル生成部１２は、多段階の複数の調整パラメータの各々に対応する色予測モデルを生成する。例えば、調整対象のパラメータの種類が２個であり、これらのパラメータを５段階に調整する場合、色予測モデル生成部１２は、５×５＝２５の組合わせに対応した２５個の色予測モデルを生成する。
そして、色予測モデル生成部１２は、生成した色予測モデルを記憶部１８に書き込んで記憶させる。
【０１５１】
ステップＳ５１０：
カラープロファイル生成部１３は、記憶部１８における調整された色予測モデルを参照し、調整対象のパラメータの組合わせの色予測モデルの各々に対応するカラープロファイルを生成する。
そして、カラープロファイル生成部１３は、生成したカラープロファイルの各々を、記憶部１８に書き込んで記憶させる。
【０１５２】
ステップＳ５１１：
制御値生成部１５は、調整制御部１４が印刷物の画像データから切り出した参照画像の画素の各々におけるインキそれぞれの指令網点面積率の組合わせを、記憶部１８の調整後のカラープロファイルの各々を順次参照し、カラープロファイル毎の予想測色値に変換する。すなわち、制御値生成部１５は、カラープロファイル毎の予想測色値に変換された参照画像のデータを生成する。
また、制御値生成部１５は、記憶部１８の印刷デバイスのデバイスカラープロファイルを参照し、各参照画像のデータにおける画素の各々の予想測色値を印刷デバイスの制御値に変換し、記憶部１８に書き込んで記憶させる。これにより、制御値生成部１５は、カラープロファイルの数の異なる参照画像の制御値のデータを生成する。
そして、印刷データ出力部１７は、記憶部１８から参照画像の各々の各画素の制御値を読み出し、読み出した各参照画像の制御値を、印刷デバイスに対して出力する。
【０１５３】
これにより、ステップＳ５０６に処理が移行し、ユーザは、印刷デバイスの印刷した複数の参照画像の印刷物の各々と、印刷機の印刷物における参照画像に対応する領域とのそれぞれの色を観察して比較する。そして、ユーザは、印刷デバイスの印刷した参照画像の印刷物のいずれかの色が、印刷機の印刷物における参照画像に対応する領域の色に対して、自身の許容範囲の違いに収まっている（すなわち、一致している）か否かの判定を行う。
そして、ステップＳ５０７において、印刷機の印刷物と色が一致している参照画像があれば処理をＳ５０８へ進め、一致している参照画像がなければ、処理をステップＳ５０９へ進める。
【０１５４】
＜色予測モデルのパラメータ調整＞
ユーザが印刷機の印刷物と印刷デバイスの印刷物との各々の色を比較し、一致していないと判定した場合、すでに述べたように色予測モデルのパラメータの調整を行う。
調整対象のパラメータとしては、調整したい項目に対応して以下に示す種類がある。
【０１５５】
ａ．印刷の濃度を調整する場合
調整するパラメータは、（１）式に示すクベルカ・ムンクの式における膜厚係数Ｘｍの数値の変更である。調整対象のインキの色の濃度を濃くしたいのであれば、膜厚係数Ｘｍを大きくし、一方、調整対象のインキの色の濃度を薄くしたいのであれば、膜厚係数Ｘｍを小さくする。調整制御部１４は、多段階に変更する場合、所定の比率で変化させ、例えば、濃くする変化の場合、現在の膜厚係数を５段階に調整するのであれば、膜厚係数Ｘｍに対して１．０５、１．１０、１．２５、１．２０、１．２５などを乗算し、微少変化させた調整用のパラメータを生成する。また、調整制御部１４は、薄くする変化の場合、現在の膜厚係数を５段階に調整するのであれば、膜厚係数Ｘｍに対して０．９５、０．９０、０．８５、０．８０、０．７５などを乗算し、微少変化させた調整用のパラメータとして用いる。色予測モデル生成部１２は、調整対象のインキの色の濃度を多段階に調整した色予測モデルを生成する。
上記比率は、ユーザが入力手段で入力しても良いし、予め設定しておいてもよい。
【０１５６】
ｂ．印刷の色彩値を調整する場合
印刷媒体のいずれかのインキの色彩値を変更したい場合、印刷物に用いたインキのベタ印刷の測色値により選択したインキ自体を他のインキに変更する。ユーザがこのパラメータの変更の指示を入力した場合、特色インキ分光反射率算出部１２１はベタ印刷と近い測色値のインキ（特色インキ）を、すでに選択されている以外に、測色値が近い順に複数個を近似色データベース１２８から抽出し、抽出した複数個を表示部１６の表示画面に表示する。インキの色彩値を多段階に変更するため、表示された特色インキから、ユーザが多段階の数の特色インキを選択する。
【０１５７】
これより、調整制御部１４は、選択された特色インキの各々に対して、（４）式において特色インキにおける配合されたインキの配合比（α及びβの値）を変化させて、この特色インキの吸収係数Ｋ（λ）及び散乱係数Ｓ（λ）を調整し、配合比毎に得られる分光反射率を算出し、この分光反射率から求めた予測測色値が、印刷物の印刷に用いたインキのベタ印刷の測色値と同一となる配合比を選択する。色予測モデル生成部１２は、この選択された特色インキと、このインキの配合比による特色インキの吸収係数Ｋ（λ）及び散乱係数Ｓ（λ）を用いて、色予測モデルを生成する。この結果、色予測モデル生成部１２は、ユーザが選択した複数の特色インキと、この特色インキの吸収係数Ｋ（λ）及び散乱係数Ｓ（λ）の組合わせに対応した色予測モデルを生成する。
【０１５８】
ｃ．印刷の階調特性を調整する場合
以下の（９）式における出現率関数ａ１（ｓ）、ａ２（ｓ）、ａ３（ｓ）、…、ａｍ（ｓ）の各々における関数内のパラメータを変更し、図３（ａ）に示す出現率関数の各々の形状及び位置などを調整する。
これにより、各インキの濃度の段階の階調特性を任意に変化させることができる。この関数内のパラメータを複数段階に変更し、複数個の階調特性の図７に対応するテーブルを作成し、それぞれのテーブルに対応する色予測モデルを生成する。
【０１５９】
ｄ．印刷物のトラッピング状態に合わせる調整を印刷デバイスに対して行う場合
インキを重ね合わせた構成において、上側に重ねるインキの濃度階調分光反射率ＲＫＭを（１）式のクベルカ・ムンクの式で求める際、膜厚係数Ｘｍを変更する。この膜厚係数Ｘｍの変更については、すでに説明した「ａ．印刷の濃度を調整する場合」の場合と同様である。そして、色予測モデル生成部１２は、調整対象の重なりの上側のインキの色の濃度を多段階に調整し、トラッピングの状態を変更させた色予測モデルを生成する。
【０１６０】
ｅ．印刷物における各版の重ね合わせのずれの調整を印刷デバイスに対して行う場合
調整制御部１４は、記憶部１８に記憶されている版の画像データを読み出し、それぞれの版を合成して印刷物の画像データを生成させる際、各版を完全に一致させて重ねるのではなく、所定のずれ（微少なずれ）量を持たせて重ね合わせる。これにより、各版の重なりにより、各画素の指令網点面積率の組合わせが異なる。複数個の異なるずれ量により、この個数に対応する数の画像データが生成される。
そして、制御値生成部１５は、この複数の印刷物の画像データの各々の画素におけるインキそれぞれの指令網点面積率の組合わせを、記憶部１８のカラープロファイルを参照して予想測色値に変換する。制御値生成部１５は、記憶部１８の印刷デバイスのデバイスカラープロファイルを参照し、画素の各々の予想測色値を印刷デバイスの制御値に変換し、記憶部１８に書き込んで記憶させる。
【０１６１】
ｆ．印刷物の粒状感に合わせる調整を印刷デバイスに対して行う場合
調整制御部１４は、画像データに対してノイズを重畳するフィルタ処理、例えばガウシアンノイズを重畳させるフィルタ処理を行うことにより、画像データにおける画素の各々の指令網点面積率を変化させ、画像に粒状のノイズを重畳させる。この重畳させるノイズの形状を多段階に変更し、調整された複数の画像データを得ることができる。
そして、制御値生成部１５は、この複数の印刷物の画像データの各々の画素におけるインキそれぞれの指令網点面積率の組合わせを、記憶部１８のカラープロファイルを参照して予想測色値に変換する。制御値生成部１５は、記憶部１８の印刷デバイスのデバイスカラープロファイルを参照し、画素の各々の予想測色値を印刷デバイスの制御値に変換し、記憶部１８に書き込んで記憶させる。
【０１６２】
本実施形態によれば、上述したように、印刷機でグラビア印刷、スクリーン印刷及びオフセット印刷などの手法で印刷媒体に対してインキを網点形状に印刷して画像の色を再現する網点印刷物と、同様な色で同一の画像データをインクジェットプリンタやデジタル印刷機など簡易な印刷機（印刷デバイス）で印刷する際、網点形状をモデル化した色予測モデルのパラメータを調整することにより、印刷デバイスの印刷物の色調整を行うことにより、印刷デバイスの制御値を変更することで、印刷デバイスの印刷物の色調整を、網点印刷物を印刷したインキ単位に対応して行うことが可能となり、網点印刷物の所定のインキに対応した印刷の調整を行う際、この所定のインキ以外の他のインキの色に対して影響を与えないため、印刷デバイスの印刷物における調整対象の色の各々の調整を容易に行うことができる。
【０１６３】
また、本実施形態によれば、大量に印刷が可能な商業用の印刷機で印刷した印刷物を、少量追加したい際、少量印刷ではコストがかかる商業用の印刷機でなく、インクジェットプリンタやデジタル印刷機などの簡易な印刷機を用い、商業用の印刷機で印刷した印刷物と同様な印刷物を印刷する可能となり、低コストかつリアルタイムに、少量の印刷物の追加を行うことができる。
また、本実施形態によれば、大量に印刷が可能な商業用の印刷機で印刷した印刷物が足りなくなった場合、商業用の印刷機で印刷した印刷物と同様な印刷物を、低コストかつリアルタイムに簡易な印刷機を用いて印刷することが可能となり、従来のように、足りなくなることを予測して余分に印刷して在庫を持つ必要がなくなる。
【０１６４】
また、本実施形態においては、特色インキの吸収特性及び散乱特性を原色インキから推定し、特色インキの分光反射率を算出し、印刷物の特色インキの重ね刷りされた領域の分光反射率を求め、印刷物の各画素の測色値を推定する色予測モデルを用いているが、指令網点面積率の組合わせから測色値を推定できる色予測モデルであれば、いずれの構成を用いても良い。
【０１６５】
なお、本発明における図１の印刷色調整システムの機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより印刷デバイスの印刷する印刷物の色調整の制御を行ってもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。
また、「コンピュータシステム」は、ホームページ提供環境（あるいは表示環境）を備えたＷＷＷシステムも含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムが送信された場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリ（ＲＡＭ）のように、一定時間プログラムを保持しているものも含むものとする。
【０１６６】
また、上記プログラムは、このプログラムを記憶装置等に格納したコンピュータシステムから、伝送媒体を介して、あるいは、伝送媒体中の伝送波により他のコンピュータシステムに伝送されてもよい。ここで、プログラムを伝送する「伝送媒体」は、インターネット等のネットワーク（通信網）や電話回線等の通信回線（通信線）のように情報を伝送する機能を有する媒体のことをいう。また、上記プログラムは、前述した機能の一部を実現するためのものであっても良い。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であっても良い。
【符号の説明】
【０１６７】
１…印刷色調整システム
１１…入力制御部
１２…色予測モデル生成部
１３…カラープロファイル生成部
１４…調整制御部
１５…制御値生成部
１６…表示部
１７…印刷データ出力部
１８…記憶部
１０１…入力部
１０２…濃度階調分光反射率算出部
１０３，２０３…濃度階調出現率算出部
１０４，２０４…濃度階調出現率テーブルデータベース
１０５…測定分光反射率データベース
１０６…吸収係数・散乱係数データベース
１０７…吸収係数・散乱係数算出部
１０８…分光反射率予測部
１１０…出力部
１１１…分光光学濃度予測部
１１２…重み係数算出部
１１３…重み係数データベース
１１４…混合予測部
１２１…特色インキ分光反射率算出部
１２２…特色インキ配合比決定部
１２３…特色インキ濃度階調出現率算出部
１２５…一時記憶部
１２６…予測パラメータデータベース
１２８…近似色データベース
２０２…濃度階調分光光学濃度算出部
１０８１，１１１１…拡張ノイゲバウア原色算出部
１０８２，１１１２…拡張ノイゲバウア原色出現率算出部
１０８３，１１１４…分光反射率算出部
１１１３…分光光学濃度算出部
【図１】
2019016924.tif 000011
【図２】
2019016924.tif 000012
【図３】
2019016924.tif 000013
【図４】
2019016924.tif 000014
【図５】
2019016924.tif 000015
【図６】
2019016924.tif 000016
【図７】
2019016924.tif 000017
【図８】
2019016924.tif 000018
【図９】
2019016924.tif 000019
【図１０】
2019016924.tif 000020