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