JP6326955B2 - 色予測システムおよび色予測方法 - Google Patents

Description

本発明は、印刷機が印刷する再現色のカラーマッチングにおいて、この再現色を予測する色予測システムおよび色予測方法に関する。

例えば、カラーマッチングにおいて面積変調（面積変調階調表現）により再現色を調整する場合、与えられた網点面積率で各色が印刷される部分の色を予測する方法として、例えば、ユール・ニールセン分光ノイゲバウアモデル（ Ｙｕｌｅ−Ｎｉｅｌｓｅｎ Ｓｐｅｃｔｒａｌ Ｎｅｕｇｅｂａｕｅｒ Ｍｏｄｅｌ )を利用する方法（例えば、特許文献１参照）や、分光拡張ノイゲバウアモデルを利用する方法（例えば、非特許文献１参照）などの技術が知られている。
ここでいう「分光拡張ノイゲバウアモデル」とは、ノイゲバウア原色と網点面積率を波長の関数となるように分光に拡張したモデルである。

また、上述したモデルは、網点の大きさが異なっていても、網点が原色単体ベタの均一な厚さのインキにより印刷媒体表面に印刷されていることを前提としている。このため、網点におけるドットのエッジにおける光の散乱による構成分を正確に表現することができていない。

この対応として、この「分光拡張ノイゲバウアモデル」を拡張し、網点におけるドットのエッジにおける光の散乱構成分を、より正確に考慮するためモデルもある（例えば、特許文献２参照）。

特表２００７−５１６６６３号公報 特開２００６−３２９７５３号公報

日本印刷学会誌 第４２巻第５号（２００５）

上述した従来例は原色インキを重ね刷りして行った印刷において、印刷媒体の表面に再現される色の予測を行うものである。
一方、特色インキは使用する都度、複数の原色インキを任意の配合比で配合して生成するものであり、原色インキのように実際に測定された分光反射率が得ることができない。このため、特色インキの階調毎の分光反射率を求められず、上述の「分光拡張ノイゲバウアモデル」を拡張したモデルを用いて色の算出を行うことができない。
また、上述の従来例においては、印刷媒体表面に原色インキを組合わせた特色インキを重ね刷りした印刷によって再現される色の予測を行うことが困難である。

本発明は、このような状況に鑑みてなされたもので、特色インキを印刷媒体に重ね刷りした印刷結果における再現色を予測する色予測システム及び色予測方法を提供する。

本発明の色予測システムは、原色インキ毎の吸収係数及び散乱係数が予め書き込まれて記憶されている吸収係数・散乱係数データベースと、色を再現したい色見本から測定した測色値を再現するのに好適な原色インキの組合せで構成される特色インキを選択し、当該特色インキに含まれる原色インキの吸収係数及び散乱係数を前記吸収係数・散乱係数データベースから読み出し、指定の配合比に応じて特色インキの吸収係数及び散乱係数を求め、求めた吸収係数及び散乱係数から再現したい特色インキの分光反射率を算出する特色インキ分光反射率算出部と、特色インキ分光反射率算出部が算出した分光反射率から再現色を求める特色インキ再現色算出部と、前記色見本から測定した測色値との差分を求め、適宜前記指定の配合比を修正して当該差分が所定の許容範囲以内となる配合比を算出する特色インキ配合比決定部と、前記特色インキを重ね刷りした際のノイゲバウア原色を、前記特色インキ配合比決定部が決定した特色の配合比に基づいて算出した特色インキ各々の吸収係数及び散乱係数から算出する分光データ算出部とを備えることを特徴とする。

本発明の色予測システムは、原色インキ毎に指令網点面積率に応じた濃度階調領域の出現率を示す出現率関数が記憶された濃度階調出現率テーブルデータベースと、前記特色インキに含まれる前記原色インキの前記出現率関数から、前記原色インキの配合情報に基づいて前記特色インキの出現率関数を求める特色濃度階調出現率算出部と、特色インキの濃度階調領域の分光反射率を求める特色濃度階調分光反射率算出部とを備え、前記分光データ算出部が、前記特色インキの濃度階調領域及び当該濃度階調領域の各々の出現率と、分光反射率との対応関係を示す計算モデルにより、前記指令網点面積率により形成された網点による前記特色の再現色を示す分光反射率を算出することを特徴とする。

本発明の色予測システムは、複数の特色インキの網点から形成される再現色を推定する場合、前記分光データ算出部により求めた前記分光反射率や前記出現率関数に基づいて各々の網点における重なり合う領域の組合わせとその出現率を求め、網点が他の網点の上部に印刷された分光反射率を、推定することを特徴とする。

本発明の色予測システムは、前記分光データ算出部が前記計算モデルから求めた前記特色インキの分光反射率と、前記分光データ算出部が分光光学濃度を用いて求めた分光反射率とを混合する混合比を含む色予測パラメータが記憶されている予測パラメータデータベースをさらに有することを特徴とする。

本発明の色予測システムは、前記分光データ算出部に対し、重ね合わせて印刷する前記特色インキ毎の指令網点面積率の組合せ毎に分光反射率を算出させ、色空間における再現色のプロファイルを作成する色予測テーブル作成部をさらに有することを特徴とする。

本発明の色予測システムは、色空間における連続的な座標値毎に、当該座標値の再現色を再現する特色インキ各々の指令網点面積率を求める特色分解テーブル作成部をさらに有することを特徴とする。

本発明の色予測システムは、特色インキを印刷する印刷媒体、及び重ね刷りにおける前記特色インキの印刷順、及び特色インキのインキ膜厚を任意に設定することを特徴とする。

本発明の色予測方法は、特色インキ分光反射率算出部が、色を再現したい色見本から測定した測色値を再現するのに好適な原色インキの組合せで構成される特色インキを選択し、当該特色インキに含まれる原色インキの吸収係数及び散乱係数を、原色インキ毎の吸収係数及び散乱係数が予め書き込まれて記憶されている吸収係数・散乱係数データベースから読み出し、指定の配合比に応じて特色インキの吸収係数及び散乱係数を求め、求めた吸収係数及び散乱係数から再現したい特色インキの分光反射率を算出する特色インキ分光反射率算出過程と、特色インキ再現色算出部が、特色インキ分光反射率算出部が算出した分光反射率から再現色を求める特色インキ再現色算出過程と、特色インキ配合比決定部が、前記色見本から測定した測色値との差分を求め、適宜前記指定の配合比を修正して当該差分が所定の許容範囲以内となる配合比を算出する特色インキ配合比決定過程と、分光データ算出部が、前記特色インキを重ね刷りした際のノイゲバウア原色を、前記特色インキ配合比決定部が決定した特色の配合比に基づいて算出した特色インキ各々の吸収係数及び散乱係数から算出する分光データ算出過程とを含むことを特徴とする。

以上説明したように、本発明によれば、特色インキを配合する原色インキの配合比を求め、特色インキの吸収係数及び散乱係数を原色インキの吸収係数及び散乱係数から求め、色見本の色に近い特色インキの分光反射率を生成し、色見本の色に近い特色インキの単色の再現色を得るため、特色インキを印刷媒体に重ね刷りした印刷結果における再現色を予測することが可能となる。

本発明の第１の実施形態による色予測システム１の構成例を示すブロック図である。 指令網点面積率と、指令網点面積率で形成される網点における濃度階調領域の出現率との対応関係を説明する図である。 特色インキを構成する原色インキの配合比を決定するフローチャートである。 下地としてのインキの上に重ねて印刷されるインキの濃度階調分光反射率Ｒ_ｉｍ（λ）の算出を説明する図である。 特色インキの重ね合わせによりノイゲバウア原色を算出する処理を示すフローチャートである。 下地のインキの濃度階調領域と、下地のインキ（特色インキまたは原色インキ）に重ねて印刷されるインキの濃度階調領域との重なりの出現率の割合の算出を説明する図である。 図６における領域Ｑ１から領域Ｑ９のそれぞれの出現率の算出結果を示すテーブルの図である。 本実施形態における色予測テーブル作成部１０９の色予測テーブルの生成処理の動作例を説明するフローチャートである。 本実施形態における特色分解テーブル作成部１３０の特色分解テーブルの生成処理の動作例を説明するフローチャートである。 本発明の第２の実施形態による色予測システム２の構成例を示すブロック図である。 特色インキの階調毎の分光反射率を算出する処理を示すフローチャートである。

本発明においては、グラビア印刷、スクリーン印刷及びオフセット印刷などにおいて、印刷媒体に対してインキを印刷する際、階調度を示す指令網点面積率に対応した網点が印刷媒体表面に形成される。このグラビア印刷における網点は、インキが印刷される印刷媒体表面の面積（面積変調階調表現）と、印刷されるインキの厚さ（濃度変調階調表現）とが、指令網点面積率に応じて変化する。例えば、印刷において形成される網点は、山の構造に相似しているということもでき、大きな山ほど裾野が広く高さも高く、この大きな山に比較して小さな山ほど裾野が狭く高さも低い。すなわち、印刷における網点は、インキが印刷される面積だけでなく、印刷されるインキの厚さも、指令網点面積率に応じて変化する。

このため、本実施形態においては、印刷における網点の構造を、指令網点面積率に対応して、インキが印刷された網点の印刷媒体表面の面積と、印刷される網点におけるインキの厚さとを表す計算モデル（後述するコアフリンジモデル）を用いる。本実施形態においては、印刷媒体（例えば、用紙）に対して印刷された単色あるいは特色の各階調度の色予測を行う際、指令網点面積率で生成される各階調度の網点が、複数の濃度階調領域（等高線と同様な構成）から生成されている網点形状をモデル化した計算モデルを用いる。

以下、本発明の第１の実施形態について、図面を参照して説明する。
図１は、本発明の第１の実施形態による色予測システム１の構成例を示すブロック図である。この図１において、色予測システム１は、入力部１０１、濃度階調分光反射率算出部１０２、濃度階調出現率算出部１０３、濃度階調出現率テーブルデータベース１０４、測定分光反射率データベース１０５、吸収係数・散乱係数データベース１０６、吸収係数・散乱係数算出部１０７、分光反射率予測部１０８、色予測テーブル作成部１０９、出力部１１０、分光光学濃度予測部１１１、混合予測部１１４、特色インキ分光反射率算出部１２１、特色インキ配合比決定部１２２、特色インキ濃度階調出現率算出部１２３、特色分解部１２４、一時記憶部１２５、予測パラメータデータベース１２６、網点面積率初期値セットデータベース１２７、近似色データベース１２８、特色インキ再現色算出部１２９、特色分解テーブル作成部１３０、濃度階調分光光学濃度算出部２０２、濃度階調出現率算出部２０３及び濃度階調出現率テーブルデータベース２０４を備えている。

入力部１０１は、例えば、外部のコンピュータと接続されており、ユーザが設定する原色インキや特色インキの各原色の指令網点面積率の指定値などのデータを入力する。
また、入力部１０１は、キーボードあるいはタッチパネルなどの入力手段を有し、この入力手段から入力されるユーザが設定するインキの各原色の指令網点面積率の指定値などのデータを、色予測システム１内部の各部に出力する構成としても良い。

吸収係数・散乱係数データベース１０６には、原色インキ毎に、指令網点面積率が１００％、すなわちベタで印刷媒体（例えば、コート紙などの用紙）に印刷された印刷部分から求めたインキの着色層の散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）が、上記外部のコンピュータなどにより予め書き込まれて記憶されている。ここで、散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）を求める際、白地及び黒地の各々の印刷媒体に対して、それぞれ原色インキをベタで印刷して上記印刷部分を作成する。
そして、白地及び黒地の各々の印刷媒体表面にベタで印刷されたインキの着色層の分光反射率をそれぞれ測定する。ここで、白地の印刷媒体表面における印刷部分の分光反射率を白地測定分光反射率とし、黒地の印刷媒体表面における印刷部分の分光反射率を黒地測定分光反射率とする。
この求めた白地測定分光反射率及び黒地測定分光反射率から、原色インキをベタで印刷媒体に印刷した印刷部分における散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）を求める。また、散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）は、所定の波長の範囲で複数の波長λで求められている。

近似色データベース１２８には、配合された原色インキの組合せが異なる複数種類の特色インキ（参照特色インキ）について、この特色インキの測色値と、原色の配合比が予め書き込まれて記憶されている。

測定分光反射率データベース１０５には、原色インキ毎の複数の指令網点面積率で印刷媒体に印刷された印刷部分の測定分光反射率Ｒ_ｓ（λ）が上記外部のコンピュータなどにより予め書き込まれて記憶されている。例えば、この印刷部分はｍ段階の指令網点面積率で印刷されたステップチャートとなっている。測定分光反射率Ｒ_ｓ（λ）を求める際、原色インキを用いて、複数の指令網点面積率の網点を、実際の印刷に用いる印刷媒体に対してそれぞれ印刷する。ここで、ｓは指令網点面積率である。そして、印刷媒体に印刷された印刷部分の分光反射率を、指令網点面積率の網点を有する印刷部分毎に測定する。また、測定分光反射率データベース１０５には、実際の印刷に用いられる印刷媒体の下地分光反射率Ｒ_０（λ）が、上述した測定分光反射率Ｒ_ｓ（λ）と同様に、上記外部のコンピュータなどにより予め書き込まれて記憶されている。

一時記憶部１２５は、本実施形態における各部の途中計算結果が記憶される記憶部であり、内部に散乱吸収係数テーブル、特色インキ配合比テーブル、特色インキ濃度階調出現率テーブル、ノイゲバウア原色テーブル、分光反射率テーブル、色予測テーブル、特色分解テーブルなどが書き込まれて記憶されている。
予測パラメータデータベース１２６には、本実施形態においては混合予測部１１４で用いる後述する重み係数ｗが、記憶されている。この重み係数ｗは、教師となるデータを印刷して、重み係数により作成したモデルから得られた分光反射率と、印刷された教師データとしての分光反射率との誤差が最小となるようにして求められている。また、印刷色の分光反射率予測モデルとしてユール・ニールセン修正ノイゲバウアモデルを用いる場合には、ｎ値が設定されていても良い。

網点面積率初期値セットデータベース１２７には、特色分解部１２４が使用する指令網点面積率の初期値が、使用する特色インキの色数毎に予め書き込まれて記憶されている。

＜分光反射率による濃度階調出現率の算出＞
濃度階調分光反射率算出部１０２は、入力部１０１より供給される原色インキの種類に応じて、吸収係数・散乱係数データベース１０６から、対応する原色インキの吸収係数Ｋ（λ）及び散乱係数Ｓ（λ）を読み出す。また、濃度階調分光反射率算出部１０２は、測定分光反射率データベース１０５から、入力部１０１より供給される原色インキの種類及び印刷媒体の種類に応じて、指令網点面積率毎の測定分光反射率Ｒ_ｓ（λ）と、印刷媒体の下地分光反射率Ｒ_０（λ）との各々を読み出す。

そして、濃度階調分光反射率算出部１０２は、吸収係数・散乱係数データベース１０６測定分光反射率データベース１０５から読み出した原色インキの吸収係数Ｋ（λ）及び散乱係数Ｓ（λ）と、印刷媒体の下地分光反射率Ｒ_０（λ）と濃度階調の膜厚係数Ｘ_ｍ（後述するように、膜厚Ｘ_ｍは濃度階調領域毎に変更する）との各々を、下記の（１）式（クベルカ・ムンクの式）に代入し、濃度階調分光反射率Ｒ_ｉ１（λ）、Ｒ_ｉ２（λ）、Ｒ_ｉ３（λ）、…、Ｒ_ｉｍ（λ）の各々を算出する。
以下の（１）式において、ａ（λ）は、散乱係数Ｓ（λ）と吸収係数Ｋ（λ）とを加算し、加算結果を散乱係数Ｓ（λ）により除算した数値である。ｂ（λ）は、ａ（λ）を二乗した数値から１を減算し、減算結果の平方根を計算した数値である。

本実施形態において、上述したクベルカ・ムンクの式、すなわち（１）式における印刷されたインキの膜厚係数Ｘ_ｍは、印刷媒体に原色インキがベタで印刷された印刷部分を基にし、印刷部分の濃度階調を示す数値として用いる。すなわち、膜厚係数Ｘ_ｍは任意に設定されており、例えば、最もインキの膜厚が厚いベタの印刷部分を膜厚が１００％として、上記膜厚係数を１とし、この１を濃度階調領域の膜厚の段階数ｍに対応させてｍ分割する。例えば、指令網点面積率の示す濃度階調領域の膜厚の段階が５段階であれば、ｍ＝１、２、３、４、５であり、膜厚係数Ｘ_ｍは、それぞれの濃度階調領域の膜厚の段階に対応して、Ｘ_１＝１．０、Ｘ_２＝０．８、Ｘ_３＝０．６、Ｘ_４＝０．４、Ｘ_５＝０．２とする。

膜厚係数Ｘ_ｍを、すでに述べたように、下地の印刷媒体の下地分光反射率Ｒ_０（λ）と、吸収係数Ｋ（λ）と、散乱係数Ｓ（λ）とともにクベルカ・ムンクの式である（１）式に代入し、それぞれの指令網点面積率の網点が含む濃度階調領域の分光反射率として、濃度階調分光反射率Ｒ_ｉ１（λ）、Ｒ_ｉ２（λ）、Ｒ_ｉ３（λ）、…、Ｒ_ｉｍ（λ）の各々を算出する。この濃度階調分光反射率Ｒ_ｉ１（λ）、Ｒ_ｉ２（λ）、Ｒ_ｉ３（λ）、…、Ｒ_ｉｍ（λ）が、後述する計算モデルで用いる、網点を構成する複数の濃度階調領域の各々の分光反射率として用いられる。

濃度階調出現率算出部１０３は、濃度階調分光反射率算出部１０２から、濃度階調分光反射率Ｒ_ｉ１（λ）、Ｒ_ｉ２（λ）、Ｒ_ｉ３（λ）、…、Ｒ_ｉｍ（λ）の各々を読み込む。また、濃度階調出現率算出部１０３は、指令網点面積率毎の測定分光反射率Ｒ_ｓ（λ）の各々を、測定分光反射率データベース１０５から読み込む。
そして、濃度階調出現率算出部１０３は、以下の（２）式（計算モデル）に対して、濃度階調分光反射率Ｒ_ｉ１（λ）、Ｒ_ｉ２（λ）、Ｒ_ｉ３（λ）、…、Ｒ_ｉｍ（λ）の各々を代入し、後述する処理により算出分光反射率Ｒ’（ｓ，λ）を算出する。

ここで、濃度階調出現率算出部１０３は、以下の（３）式において、出現率（用紙の印刷部分における網点を構成する各濃度階調の濃度階調領域の面積の比率）の数値を変更しつつ、算出分光反射率Ｒ’（ｓ，λ）を算出する。そして、濃度階調出現率算出部１０３は、算出された算出分光反射率Ｒ’（ｓ，λ）の各々と、測定分光反射率Ｒ_ｓ（λ）との平均二乗誤差ＲＭＳＥを、指令網点面積率毎に所定の波長範囲において求める。
濃度階調出現率算出部１０３は、算出分光反射率Ｒ’（ｓ，λ）と測定分光反射率Ｒ_ｓ（λ）との平均二乗誤差が最も小さくなる濃度階調領域の各々の出現率を求める。ここで、ｓは指令網点面積率である。

そして、濃度階調出現率算出部１０３は、濃度階調領域の各々の出現率により、各濃度階調領域の出現率関数ａ_１（ｓ）、ａ_２（ｓ）、ａ_３（ｓ）、…、ａ_ｍ（ｓ）を求める。ここで、濃度階調出現率算出部１０３は、濃度階調領域毎に、得られた出現率を用いて、この出現率を指定網点面積率ｓの２次関数などをフィッティングし、出現率関数ａ_１（ｓ）、ａ_２（ｓ）、ａ_３（ｓ）、…、ａ_ｍ（ｓ）を求めてもよい。濃度階調出現率算出部１０３は、求めた各濃度階調領域の出現率関数ａ_１（ｓ）、ａ_２（ｓ）、ａ_３（ｓ）、…、ａ_ｍ（ｓ）を、濃度階調出現率テーブルデータベース１０４に書き込んで記憶させる。

上記（３）式において、濃度階調出現率算出部１０３は、波長λが３８０ｎｍから７３０ｎｍをｎ分割した刻み幅により、各波長λにおける誤差の二乗を加算した平均二乗誤差ＲＭＳＥを、指令網点面積率毎に求めている。
上述したように、濃度階調出現率算出部１０３は、（２）式の計算モデルを用いて、指令網点面積率によるインキが印刷された印刷部分において、網点が含む濃度階調毎に、濃度階調の濃度階調分光反射率Ｒ_ｉｍ（λ）と濃度階調の出現率関数ａ_１（ｓ）、ａ_２（ｓ）、ａ_３（ｓ）、…、ａ_ｍ（ｓ）とを乗算した結果を加算して、指令網点面積率の印刷部分の分光反射率Ｒ’（ｓ，λ）（算出分光反射率）を算出する。

図２は、指令網点面積率と、指令網点面積率で形成される網点における濃度階調領域の出現率との対応関係を説明する図である。図２（ａ）において、横軸が指令網点面積率ｓを示し、縦軸が出現率ａを示している。出現率関数ａ_ｍ（ｓ）により、指令網点面積率における各濃度階調領域の出現率を求めることができる。濃度階調出現率算出部１０３は、すでに説明したように濃度階調領域毎の出現率関数ａ_ｍ（ｓ）の各々を、（３）式で求めた指令網点面積率に対応した濃度階調領域毎の出現率を２次関数の近似式として関数化することにより求める。この図２（ａ）は、ｍ＝４の場合、すなわち４階調の濃度階調の濃度階調領域の各々の出現率関数ａ_１（ｓ）、ａ_２（ｓ）、ａ_３（ｓ）、ａ_４（ｓ）の各々を示している。

また、図２（ｂ）は、図２（ａ）における指令網点面積率ｓ_１、ｓ_２、ｓ_３、ｓ_４の各々における網点の平面視における形状を示す図である。指令網点面積率ｓ_１においては、濃度階調領域Ｐ_１のみが形成されている。また、指令網点面積率ｓ_２においては、濃度階調領域Ｐ_１の内部に濃度階調領域Ｐ_２が形成される。指令網点面積率ｓ_３においては、濃度階調領域Ｐ_２のみが形成されている。また、指令網点面積率ｓ_４においては、濃度階調領域Ｐ_３の内部に濃度階調領域Ｐ_４が形成される。上述したように、本実施形態においては、（２）式の計算モデルを用いて、この図２（ｂ）に示す網点構造として、グラビア印刷における網点の構造をモデル化している。

図１に戻り、濃度階調出現率算出部１０３は、求めた濃度階調領域の指令網点面積率における各濃度階調領域の出現率を示す出現率関数ａ_ｍ（ｓ）を、濃度階調出現率テーブルデータベース１０４に対し、濃度階調領域の濃度階調分光反射率と対応させて書き込んで記憶させる。同様に、濃度階調出現率算出部１０３は、他の原色インキの各々に対しても、濃度階調領域の指令網点面積率における各濃度階調領域の出現率を示す出現率関数ａ_ｍ（ｓ）を、濃度階調出現率テーブルデータベース１０４に対して、書き込んで記憶させる。

特色インキ分光反射率算出部１２１は、近似色データベース１２８における原色インキの組み合わせが異なる特色インキから、測色値が色見本の測色値と最も近い特色インキを読み出す。この近似色データベース１２８には特色インキの測色値と原色の配合比が記憶されており、読み出した特色インキの原色の配合比を、色見本を再現する特色インキの原色の配合比として設定する。
また、特色インキ分光反射率算出部１２１は、測定分光反射率データベース１０５から、印刷媒体の下地分光反射率Ｒ_０（λ）を読み出す。
吸収係数・散乱係数算出部１０７は、吸収係数・散乱係数データベース１０６から、設定した配合比に含まれる原色インキを構成する原色インキ、例えば原色インキ＃１及び原色インキ＃２の各々の吸収係数Ｋ_１（λ）及びＫ_２（λ）と、散乱係数Ｓ_１（λ）及びＳ_２（λ）との各々を読み出す。

そして、吸収係数・散乱係数算出部１０７は、以下の（４）式により、原色インキを混色した特色インキの吸収係数Ｋ_ｔ（λ）及び散乱係数Ｓ_ｔ（λ）を算出する。この特色インキの場合、混ぜ合わせる原色インキの比率に応じて、以下の（４）式において、特色インキの散乱係数Ｓ_ｔ（λ）及び吸収係数Ｋ_ｔ（λ）を算出する。また、特色インキ分光反射率算出部１２１は、求めた特色インキの散乱係数Ｓ_ｔ（λ）及び吸収係数Ｋ_ｔ（λ）の各々を、一時記憶部１２５の散乱吸収係数テーブルに書き込んで記憶させる。

上記（４）式において、係数α及び係数βの各々は、原色インキ＃１及び原色インキ＃２それぞれの混ぜ合わせる比率を示す係数である。原色インキ＃１の吸収係数Ｋ_１（λ）に対して係数αを乗算し、原色インキ＃２の吸収係数Ｋ_２に対して係数βを乗算し、加算した数値を、特色インキの吸収係数Ｋ_ｔ（λ）としている。同様に、原色インキ＃１の散乱係数Ｓ_１（λ）に対して係数αを乗算し、原色インキ＃２の散乱係数Ｓ_２（λ）に対して係数βを乗算し、加算した数値を、特色インキの散乱係数Ｓ_ｔ（λ）としている。

これにより、特色インキ分光反射率算出部１２１は、（１）式に対して、吸収係数Ｋ_ｔ（λ）及び散乱係数Ｓ_ｔ（λ）と、印刷媒体の下地分光反射率Ｒ_０（λ）と濃度階調の膜厚係数Ｘ_ｍとの各々を代入し、特色インキの分光反射率Ｒ_ＫＭ（λ）を求める。
特色インキ再現色算出部１２９は、特色インキ分光反射率算出部１２１が算出した特色インキの分光反射率Ｒ_ＫＭ（λ）に対して観察環境の光源の分光分布と、標準観測者を設定し、測色値（例えばＬ*ａ*ｂ*値）に変換し、特色インキ配合比決定部１２２に対して出力する。

特色インキ配合比決定部１２２は、特色インキ再現色算出部１２９が求めた測色値と、色見本の測色値との色差を確認する。そして、特色インキ配合比決定部１２２は、この色差が予め設定された許容範囲内である場合、特色インキ分光反射率算出部１２１が求めた色における原色インキ各々の配合比を、見本の特色インキを構成する原色インキの種類と、その配合比として、一時記憶部１２５の特色インキ配合比テーブルに対し、特色インキの識別情報である特色インキ識別情報とともに書き込んで記憶させる。

図３は、特色インキを構成する原色インキの配合比を決定し、特色インキの吸収係数Ｋ_ｔ（λ）、散乱係数Ｓ_ｔ（λ）を算出するフローチャートである。
ステップＳ１０１：
ユーザは再現に用いる色見本の測色値を測定し、色予測システム１に対して入力する。

ステップＳ１０２：
特色インキ分光反射率算出部１２１は、印刷媒体の分光反射率である下地分光反射率Ｒ_０（λ）を、測定分光反射率データベース１０５から読み出す。

ステップＳ１０３：
次に、特色インキ分光反射率算出部１２１は、特色インキを構成する原色インキの組み合わせを、近似色データベース１２８から読み出す。例えば、このとき、特色インキ分光反射率算出部１２１は、色見本から取得した測色値と最も近い測色値を持つ特色インキを、近似色データベース１２８から抽出して選択する。

ステップＳ１０４：
吸収係数・散乱係数算出部１０７は、ステップＳ１０３において選択された特色インキを構成する原色インキの散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）を、吸収係数・散乱係数データベース１０６から読み出す。

ステップＳ１０５：
特色インキ配合比決定部１２２は、特色インキを構成する原色インキ各々の配合比に初期値を設定する。初期値には、ステップＳ１０３で近似色データベース１２８から抽出した特色の配合比を設定する。

ステップＳ１０６：
次に、吸収係数・散乱係数算出部１０７は、（４）式に示すように、原色インキの各々の配合比を、それぞれ原色インキの散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）に乗じる。
そして、吸収係数・散乱係数算出部１０７は、乗算結果を加算することにより、その配合比（例えば、（４）式におけるα、β）における特色インキの散乱係数Ｓ_ｔ（λ）及び吸収係数Ｋ_ｔ（λ）を算出する。

特色インキ分光反射率算出部１２１は、（１）式に対し、印刷媒体の下地分光反射率Ｒ_０（λ）と、散乱係数Ｓ_ｔ（λ）及び吸収係数Ｋ_ｔ（λ）とを代入して、特色インキの分光反射率を算出する。さらに、特色インキ再現色算出部１２９は、特色インキ分光反射率算出部１２１が算出した特色インキの分光反射率から、測色値を算出する。
特色インキ配合比決定部１２２は、色見本の測色値Ｌ*ａ*ｂ*値と、特色インキ再現色算出部１２９の求めた測色値Ｌ*ａ*ｂ*値との色差を求める。

ステップＳ１０７：
そして、特色インキ配合比決定部１２２は、上記色差が予め設定した許容範囲内であるか否かの判定を行う。
このとき、特色インキ配合比決定部１２２は、色差が予め設定した許容範囲内である場合、処理をステップＳ１０８へ進める。一方、特色インキ配合比決定部１２２は、色差が予め設定した許容範囲内にない場合、処理をステップＳ１１０へ進める。

ステップＳ１０８：
次に、特色インキ配合比決定部１２２は、色見本に対応する特色インキにおける原色インキの配合比を現在の配合比とする。
そして、特色インキ配合比決定部１２２は、特色インキにおける原色インキの種類と、原色インキ各々の配合比を、一時記憶部１２５の配合比テーブルに書き込んで記憶させる。

ステップＳ１０９：
次に、特色インキ配合比決定部１２２は、色見本に対応する特色インキの散乱係数Ｓ_ｔ（λ）及び吸収係数Ｋ_ｔ（λ）を、現在の散乱係数Ｓ_ｔ（λ）及び吸収係数Ｋ_ｔ（λ）とする。
そして、特色インキ配合比決定部１２２は、特色インキの散乱係数Ｓ_ｔ（λ）及び吸収係数Ｋ_ｔ（λ）を、一時記憶部１２５の散乱吸収係数テーブルに書き込んで記憶させる。

ステップＳ１１０：
特色インキ配合比決定部１２２は、配合比の変更回数が予め設定された規定以内か否かの判定を行う。
このとき、特色インキ配合比決定部１２２は、配合比の変更回数が予め設定された規定回数以内である場合、処理をステップＳ１１１へ進める。一方、特色インキ配合比決定部１２２は、配合比の変更回数が予め設定された規定回数を超えた場合、処理をステップＳ１１２へ進める。

ステップＳ１１１：
特色インキ配合比決定部１２２は、特色インキを構成する原色インキの各々の配合比を変更する。
そして、特色インキ配合比決定部１２２は、配合比の変更回数を計数するカウンタをインクリメント（１つカウント値を上げる）し、処理をステップＳ１０６へ進める。

ステップＳ１１２：
特色インキ配合比決定部１２２は、原色インキの組み合わせを変更するか否かの選択画面を、色予測システム１の図示しない表示部に表示する。
特色インキ配合比決定部１２２は、ユーザが原色インキの組み合わせの変更を選択した場合、処理をステップＳ１１３へ進める。一方、特色インキ配合比決定部１２２は、ユーザが原色インキの組み合わせを変更しない選択をした場合、処理をステップＳ１０８へ進める。

ステップＳ１１３：
特色インキ配合比決定部１２２は、配合比の変更回数を計数するカウンタをリセットし、計数値、すなわち変更回数を「０」とする。

ステップＳ１１４：
特色インキ配合比決定部１２２は、 特色インキ分光反射率算出部１２１に対し、原色インキの組み合わせを変更する制御信号を出力する。
特色インキ分光反射率算出部１２１は、特色インキを構成する原色インキの組み合わせを、近似色データべース１２８から新たに読み出す。

図１に戻り、特色インキ濃度階調出現率算出部１２３は、原色インキを所定の比率で混ぜ合わせて生成した特色インキの濃度階調領域毎の出現率関数ａ_ｍ（ｓ）を求める。そして、特色インキ濃度階調出現率算出部１２３は、求めた出現率関数ａ_ｍ（ｓ）を一時記憶部１２５の特色インキ濃度階調出現率テーブルに対して書き込んで記憶させる。ここで、特色インキ濃度階調出現率算出部１２３は、特色インキを構成する原色インキのいずれかの出現率関数ａ_ｍ（ｓ）を、濃度階調出現率テーブルデータベース１０４から読み出し、この特色インキの出現率関数ａ_ｍ（ｓ）とする。また、特色インキ濃度階調出現率算出部１２３は、特色インキを構成する原色インキの各々の出現率関数ａ_ｍ（ｓ）を、それぞれの原色インキの配合比により組合わせ、特色インキの出現率関数ａ_ｍ（ｓ）としても良い。

＜分光反射率予測部１０８の動作＞
次に、分光反射率予測部１０８は、拡張ノイゲバウア原色算出部１０８１、拡張ノイゲバウア原色出現率算出部１０８２及び分光反射率算出部１０８３を備えている。
拡張ノイゲバウア原色算出部１０８１は、重ね合わせる色の順番に従い、下地になるインキ（原色インキあるいは特色インキ）と、下地の表面に印刷されるインキとを決定する。また、拡張ノイゲバウア原色算出部１０８１は、重ね合わせるインキの各々の指令網点面積率を、入力部１０１から読み込む。

拡張ノイゲバウア原色算出部１０８１は、下地のインキに対して重なる原色インキの散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）を、吸収係数・散乱係数データベース１０６から読み込む。また、拡張ノイゲバウア原色算出部１０８１は、下地のインキに対して重なる特色インキの散乱係数Ｓ_ｔ（λ）及び吸収係数Ｋ_ｔ（λ）を、一時記憶部１２５の散乱吸収係数テーブルから読み込む。

次に、拡張ノイゲバウア原色算出部１０８１は、（１）式に下地となるインキ（特色インキまたは原色インキ）の網点における濃度階調領域の濃度階調分光反射率Ｒ_ｉｍ（λ）を下地分光反射率Ｒ_０（λ）として代入し、下地のインキに対して重なるインキの散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）と、濃度階調領域の膜厚係数Ｘ_ｍとの各々を代入することにより、下地のインキの上に重なって印刷されるインキの網点を構成する濃度階調領域における濃度階調分光反射率Ｒ_ｉｍ（λ）を算出する。

ここで、拡張ノイゲバウア原色算出部１０８１は、上述した下地のインキの網点の濃度階調領域と、この下地のインキの上に印刷されるインキの網点を構成する濃度階調領域の各々の重なり部分の濃度階調分光反射率Ｒ_ｉｍ（λ）の計算を、下地のインキの網点における濃度階調領域と、下地のインキの網点上に印刷されるインキ（特色インキまたは原色インキ）の網点の濃度階調領域との組合わせの全てにおいて行う。

図４は、下地としてのインキの上に重ねて印刷されるインキの濃度階調分光反射率Ｒ_ｉｍ（λ）の算出を説明する図である。
濃度階調分光反射率算出部１０２は、印刷媒体の分光反射率を下地の分光反射率Ｒ_０（λ）とし、インキ１０００の領域の吸収係数Ｋ（λ）と散乱係数Ｓ（λ）とを用いて、（１）式から印刷媒体に対して印刷されたインキ１０００（特色インキまたは原色インキ）の領域の網点の濃度階調領域の分光反射率Ｒ_ｉｍ（λ）を求める。

そして、拡張ノイゲバウア原色算出部１０８１は、領域１００２に示す分光反射率のインキ１０００の上に対し、インキの領域１００４の上部に印刷された領域１００５に示す分光反射率Ｒ_ＫＭ（λ）を算出する。ここで、拡張ノイゲバウア原色算出部１０８１は、インキ１０００の領域における網点の濃度階調領域の分光反射率Ｒ_ＫＭ（λ）を下地の分光反射率Ｒ_０（λ）とし、領域１００３に示すインキ１００１の吸収係数Ｋ（λ）と散乱係数Ｓ（λ）と、濃度階調領域の膜厚係数Ｘ_ｍを、（１）式に代入して、領域１００４におけるインキ１０００の領域に重ねて印刷されたインキ１００１（特色インキまたは原色インキ）の領域１００４の網点における濃度階調領域の分光反射率Ｒ_ＫＭ（λ）を求める。

これにより、拡張ノイゲバウア原色算出部１０８１は、下地のインキ（特色インキまたは原色インキ）の網点における濃度階調領域と、下地のインキの網点上に印刷されるインキの網点の濃度階調領域との組合わせにおいて、下地のインキの網点上に印刷される上部のインキ（特色インキまたは原色インキ）の網点の濃度階調領域の重なり領域における分光反射率Ｒ_ＫＭ（λ）を、後述するように、下地のインキの網点と上部のインキ網点とにおける濃度階調領域との組合わせの全てについて算出することにより、インキが重ねて印刷された印刷部分における分光反射率Ｒ_ＫＭ（λ）を求める。拡張ノイゲバウア原色算出部１０８１は、求めた分光反射率Ｒ_ＫＭ（λ）を一時記憶部１２５のノイゲバウア原色テーブルに書き込んで記憶させる。

図５は、特色インキの重ね合わせによりノイゲバウア原色を算出する処理を示すフローチャートである。また、重ね合わせるインキとしては特色インキだけでなく、原色インキと特色インキとを組み合わせても良い。
ステップＳ２０１：
ユーザは、特色インキの組み合わせに対応し、それぞれの組み合わせにおける特色インキの刷り順を入力し、色予測システムに対して設定する。

ステップＳ２０２：
拡張ノイゲバウア原色算出部１０８１は、印刷媒体の分光反射率である下地分光反射率Ｒ_０（λ）を、測定分光反射率データベース１０５から読み出す。

ステップＳ２０３：
拡張ノイゲバウア原色算出部１０８１は、重ね合わせる特色インキの散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）の各々を、一時記憶部１２５の散乱吸収係数テーブルから読み込む。

ステップＳ２０４：
拡張ノイゲバウア原色算出部１０８１は、ユーザが入力する膜厚補正値である膜厚係数を分光反射率に反映させるため、すなわち（１）式において用いる膜厚係数を内部の記憶部に記憶させる。

ステップＳ２０５：
拡張ノイゲバウア原色算出部１０８１は、特色インキあるいは原色インキを下地として印刷した印刷物の分光反射率Ｒ_ＫＭ（λ）を算出する。このとき、拡張ノイゲバウア原色算出部１０８１は、重ね合わせる特色インキにおけるそれぞれの膜厚係数の組み合わせの全てのノイゲバウア原色を算出する。すなわち、拡張ノイゲバウア原色算出部１０８１は、下地の印刷媒体に印刷した特色インキの分光反射率Ｒ_ＫＭ（λ）を求め、この求めた分光反射率Ｒ_ＫＭ（λ）を新たな下地の分光反射率とする。そして、拡張ノイゲバウア原色算出部１０８１は、求めた分光反射率Ｒ_ＫＭ（λ）の特色インキに対して、新たな特色インキを重ね合わせた際の分光反射率を算出し、これをノイゲバウア原色の分光反射率Ｒ_ＫＭ（λ）とする。
ここで、下地がプロセスインキのベタなど、予め分光反射率の実測値が既知である場合、媒体に印刷した際の分光反射率を算出するのではなく、実測値である分光反射率を、重ね合わせる下地の分光反射率として用いても良い。
そして、拡張ノイゲバウア原色算出部１０８１は、算出した分光反射率Ｒ_ＫＭ（λ）を、一時記憶部１２５のノイゲバウア原色テーブルに書き込んで記憶させる。

ステップＳ２０６：
拡張ノイゲバウア原色算出部１０８１は、特色インキの組み合わせの全てのノイゲバウア原色の算出が終了したか否かの判定を行う。
このとき、拡張ノイゲバウア原色算出部１０８１は、全てのノイゲバウア原色の算出が終了した場合、処理を終了する。一方、拡張ノイゲバウア原色算出部１０８１は、全てのノイゲバウア原色の算出が終了していない場合、処理をステップＳ２０７へ進める。

ステップＳ２０７：
拡張ノイゲバウア原色算出部１０８１は、特色インキの組み合わせを変更し、重ね合わせる特色インキの順番を、次の特色インキの組み合わせの順番とし、処理をステップＳ２０５に進める。

図１に戻り、拡張ノイゲバウア原色出現率算出部１０８２は、下地のインキの濃度階調領域と、下地のインキに重ねて印刷されるインキの濃度階調領域との重なりの出現率の割合を算出する。

図６は、下地のインキの濃度階調領域と、下地のインキ（特色インキまたは原色インキ）に重ねて印刷されるインキの濃度階調領域との重なりの出現率の割合の算出を説明する図である。
図６においては、簡単のため重ね合わせるインキを２種類とし、濃度階調領域の種類も２種類としている。しかしながら、重ね合わせるインキと濃度階調領域の種類との各々が３以上の複数でも、下地のインキの濃度階調領域と、下地のインキに重ねて印刷されるインキの濃度階調領域との重なりの出現率の割合の算出は、以下の説明と同様に行うことができる。

図６（ａ）は、濃度階調領域が２種類ある場合を示しており、インキ＃１の指令網点面積率（設定網％）と、指令網点面積率における網点に出現する濃度階調領域であるコア及びフリンジの各々の出現率との対応関係を示している。また、インキ＃２もインキ＃１と同様の対応関係を有している。インキ＃１及びインキ＃２の各々は、原色インキあるいは特色インキのいずれかである。また、インキ＃１及びインキ＃２の各々において、濃度階調領域が濃度階調コア領域１及び濃度階調フリンジ領域２のそれぞれ２種類があり、濃度階調コア領域１は膜厚１００％であり、濃度階調フリンジ領域２は膜厚５０％である。また、他インキも同様の対応関係を有している。

図６（ｂ）は、濃度階調コア領域１及び濃度階調フリンジ領域２の各々の重なりを示している。図６（ｂ）において、例えばインキ＃１がシアン（Ｃ）であり、インキ＃２がマジェンダ（Ｍ）である。濃度階調コア領域１及び濃度階調フリンジ領域２の各々の重なりの組み合わせとしては、領域Ｑ１から領域Ｑ９の９種類がある。領域Ｑ１は、シアンの濃度階調コア領域１のみの領域である。領域Ｑ２は、シアンの濃度階調フリンジ領域２のみの領域である。領域Ｑ３は、マジェンダの濃度階調コア領域１のみの領域である。領域Ｑ４は、マジェンダの濃度階調フリンジ領域２のみの領域である。領域Ｑ５は、シアン及びマジェンダ各々の濃度階調コア領域１が重なった領域である。領域Ｑ６は、マジェンダの濃度階調フリンジ領域２とシアンの濃度階調コア領域１が重なった領域である。領域Ｑ７は、マジェンダの濃度階調コア領域１とシアンの濃度階調フリンジ領域２が重なった領域である。領域Ｑ８は、シアン及びマジェンダの各々の濃度階調フリンジ領域２が重なった領域である。領域Ｑ９は、シアン及びマジェンダのいずれのインキも存在しない領域である。

図７は、図６における領域Ｑ１から領域Ｑ９のそれぞれの出現率の算出結果を示すテーブルの図である。図７においては、原色インキＣ（シアン）及びＭ（マジェンダ）で説明しているが、特色インキでも同様に出現率が算出される。このテーブルにおいて、Ｃはシアンの濃度階調領域であり、Ｍはマジェンダの濃度階調領域であり、ＣＭはシアン及びマジェンダの濃度階調領域が重なった領域であり、Ｗはシアン及びマジェンダのいずれのインキも存在しない領域を示している。図７において、出現率α１はシアンの濃度階調コア領域１の出現率であり、出現率α２はシアンの濃度階調フリンジ領域２の出現率である。出現率α０は、出現率α１及び出現率α２の各々を加算したものである（α０＝α１＋α２）。出現率β１はマジェンダの濃度階調コア領域１の出現率であり、出現率β２はマジェンダの濃度階調フリンジ領域２の出現率である。出現率β０は、出現率β１及び出現率β２の各々を加算したものである（β０＝β１＋β２）。

領域Ｑ１及び領域Ｑ２の単次の出現率は、シアンのインキが出現する出現率α０に対し、マジェンダのインキが出現しない率である（１−β０）を乗算したα０＊（１−β０）となり、シアンのインキのみの領域の出現率を示す。本実施形態において、＊は乗算を示している。
領域Ｑ３及び領域Ｑ４の単次の出現率は、マジェンダのインキが出現する出現率β０に対し、シアンのインキが出現しない率である（１−α０）を乗算したβ０＊（１−α０）となり、マジェンダのインキのみの領域の出現率を示す。

領域Ｑ５から領域Ｑ８の単次の出現率は、シアンのインキが出現する出現率α０に対し、マジェンダのインキが出現する出現率β０を乗算したα０＊β０となり、シアンのインキとマジェンダのインキとが重なっている領域の出現率を示す。
領域Ｑ９の単次の出現率は、シアンのインキが出現しない率である（１−α０）に対し、マジェンダのインキが出現しない率である（１−β０）を乗算した（１−α０）＊（１−β０）となり、シアンのインキとマジェンダのインキとの双方が存在しない領域の出現率を示す。

領域Ｑ１の副次の出現率は、シアンのインキの網点の濃度階調コア領域１のみの出現率を示しており、濃度階調コア領域１の出現率α１を、濃度階調コア領域１の出現率α１及び濃度階調フリンジ領域２の出現率α２の加算結果により除算した率である。
領域Ｑ２の副次の出現率は、シアンのインキの網点の濃度階調フリンジ領域２のみの出現率を示しており、濃度階調コア領域１の出現率α１を、濃度階調コア領域１の出現率α１及び濃度階調フリンジ領域２の出現率α２の加算結果により除算し、この除算結果を１から減算した率である。

領域Ｑ３の副次の出現率は、マジェンダのインキの網点の濃度階調コア領域１のみの出現率を示しており、濃度階調コア領域１の出現率β１を、濃度階調コア領域１の出現率β１及び濃度階調フリンジ領域２の出現率β２の加算結果により除算した率である。
領域Ｑ４の副次の出現率は、マジェンダのインキの網点の濃度階調フリンジ領域２のみの出現率を示しており、濃度階調コア領域１の出現率β１を、濃度階調コア領域１の出現率β１及び濃度階調フリンジ領域２の出現率β２の加算結果により除算し、この除算結果を１から減算した率である。

領域Ｑ５の副次の出現率は、シアン及びマジェンダの各々のインキの網点の濃度階調コア領域１の重なり部分の出現率を示しており、シアンのインキの網点の濃度階調コア領域１の出現率α１を、濃度階調コア領域１の出現率α１及び濃度階調フリンジ領域２の出現率α２の加算結果により除算した率と、マジェンダのインキの網点の濃度階調コア領域１の出現率β１を、濃度階調コア領域１の出現率β１及び濃度階調フリンジ領域２の出現率β２の加算結果により除算した率とを乗算した結果の率である。

領域Ｑ６の副次の出現率は、シアンのインキの網点の濃度階調コア領域１とマジェンダのインキの網点の濃度階調フリンジ領域２の重なり部分の出現率を示しており、シアンのインキの網点の濃度階調コア領域１の出現率α１を、濃度階調コア領域１の出現率α１及び濃度階調フリンジ領域２の出現率α２の加算結果により除算した率と、マジェンダのインキの網点の濃度階調コア領域１の出現率β１を、濃度階調コア領域１の出現率β１及び濃度階調フリンジ領域２の出現率β２の加算結果により除算した率を１から減算した率とを乗算した結果の率である。

領域Ｑ７の副次の出現率は、シアンのインキの網点の濃度階調フリンジ領域２とマジェンダのインキの網点の濃度階調コア領域１の重なり部分の出現率を示しており、シアンのインキの網点の濃度階調コア領域１の出現率α１を、濃度階調コア領域１の出現率α１及び濃度階調フリンジ領域２の出現率α２の加算結果により除算した率を１から減算した率と、マジェンダのインキの網点の濃度階調コア領域１の出現率β１を、濃度階調コア領域１の出現率β１及び濃度階調フリンジ領域２の出現率β２の加算結果により除算した率とを乗算した結果の率である。

領域Ｑ８の副次の出現率は、シアンのインキの網点の濃度階調フリンジ領域２とマジェンダのインキの網点の濃度階調フリンジ領域２の重なり部分の出現率を示しており、シアンのインキの網点の濃度階調コア領域１の出現率α１を、濃度階調コア領域１の出現率α１及び濃度階調フリンジ領域２の出現率α２の加算結果により除算した率を１から減算した率と、マジェンダのインキの網点の濃度階調コア領域１の出現率β１を、濃度階調コア領域１の出現率β１及び濃度階調フリンジ領域２の出現率β２の加算結果により除算した率を１から減算した率とを乗算した結果の率である。

領域Ｑ９の副次の出現率は、シアンのインキとマジェンダのインキとの双方が存在しない領域の出現率を示し、「１」である。

上述したように、使用するインキの重ね合わせに用いる出現率の式は、予め設定され図７に示すテーブルとして、濃度階調出現率テーブルデータベース１０４に書き込んで記憶されている。
図１に戻り、拡張ノイゲバウア原色出現率算出部１０８２は、重ね合わせるインキの種類と、重ね合わせるインキの網点を示す指令網点面積率との各々の組み合わせにより、濃度階調出現率テーブルデータベース１０４の図７のテーブルの式を読み込む。
また、拡張ノイゲバウア原色出現率算出部１０８２は、濃度階調出現率テーブルデータベース１０４から、シアンとマジェンダとの各々の濃度階調コア領域１及び濃度階調フリンジ領域２それぞれの出現率α１、α２、β１、β２を読み込む。

拡張ノイゲバウア原色算出部１０８１は、シアンが下地で、マジェンダがシアンに重ね合わされる場合、シアンの分光反射率を下地の分光反射率Ｒ_０（λ）として、指令網点面積率で出現するマジェンダの濃度階調領域の膜厚を用い、（１）式により、重なり部分における濃度階調領域の分光反射率Ｒ_ＫＭ（λ）を算出する。例えば、拡張ノイゲバウア原色算出部１０８１は、領域Ｑ５における分光反射率Ｒ_ＫＭ（λ）を求める際、印刷媒体上のシアンのインキの膜厚１００％における濃度階調分光反射率Ｒ_ｉｍ（λ）を下地の分光反射率Ｒ_０（λ）として、マジェンダのインキの膜厚１００％を重ねた際の重なり部分における印刷部分の分光反射率Ｒ_ＫＭ（λ）を算出する。

同様に、拡張ノイゲバウア原色算出部１０８１は、領域Ｑ７における分光反射率Ｒ_ＫＭ（λ）を求める際、印刷媒体上のシアンのインキの膜厚５０％における濃度階調分光反射率Ｒ_ｉｍ（λ）を下地の分光反射率Ｒ_０（λ）とし、マジェンダのインキの膜厚１００％を重ねた際の分光反射率Ｒ_ＫＭ（λ）を算出する。同様に、拡張ノイゲバウア原色の全ての組み合わせについて分光反射率Ｒ_ＫＭ（λ）を算出する。

そして、分光反射率算出部１０８３は、上述した領域Ｑ１から領域Ｑ９の拡張ノイゲバウア原色の各々の分光反射率に拡張ノイゲバウア原色出現率を乗算し、波長毎に加算する。
分光反射率算出部１０８３は、シアンのインキの網点が印刷された印刷媒体に対し、マジェンダのインキの網点を重ねて印刷した際における印刷部分の第１予測分光反射率Ｒ_Ｄ１（λ）を算出する。

＜分光光学濃度による濃度階調出現率の算出＞
また、濃度階調分光光学濃度算出部２０２は、入力部１０１より供給される原色インキの種類に応じて、吸収係数・散乱係数データベース１０６から、原色インキがベタで印刷媒体に印刷された印刷部分から求めた着色層の吸収係数Ｋ（λ）及び散乱係数Ｓ（λ）を読み出す。また、濃度階調分光光学濃度算出部２０２は、測定分光反射率データベース１０５から、入力部１０１より供給される原色インキの種類及び印刷媒体の種類に応じて、指令網点面積率毎の測定分光反射率Ｒ_ｓ（λ）と、印刷媒体の下地分光反射率Ｒ_０（λ）との各々を読み出す。

そして、濃度階調分光光学濃度算出部２０２は、読み出した原色インキの吸収係数Ｋ（λ）及び散乱係数Ｓ（λ）と、印刷媒体の下地分光反射率Ｒ_０（λ）と、濃度階調の膜厚係数Ｘ_ｍとの各々を、すでに説明した（１）式（クベルカ・ムンクの式）に代入し、濃度階調領域の分光反射率として、濃度階調分光反射率Ｒ_ｉ１（λ）、Ｒ_ｉ２（λ）、Ｒ_ｉ３（λ）、…、Ｒ_ｉｍ（λ）を算出する。
また、濃度階調分光光学濃度算出部２０２は、求めた濃度階調分光反射率Ｒ_ｉ１（λ）、Ｒ_ｉ２（λ）、Ｒ_ｉ３（λ）、…、Ｒ_ｉｍ（λ）の各々を、以下の（５）式により、濃度階調分光光学濃度ＯＤ_ｉ１（λ）、ＯＤ_ｉ２（λ）、ＯＤ_ｉ３（λ）、…、ＯＤ_ｉｍ（λ）にそれぞれ変換する。

濃度階調分光光学濃度算出部２０２においても、濃度階調分光反射率算出部１０２と同様に、（１）式（クベルカ・ムンクの式）における印刷されたインキの膜厚係数Ｘ_ｍは、印刷媒体に原色インキがベタで印刷された印刷部分を基にし、印刷部分の濃度階調を示す数値として用いる。すなわち、膜厚係数Ｘ_ｍは任意に設定されており、例えば、最もインキの膜厚が厚いベタの印刷部分を１００％として、上記膜厚係数を１とし、この１を濃度階調領域の膜厚の段階数ｍに対応させてｍ分割する。例えば、指令網点面積率の示す濃度階調領域の膜厚の段階が５段階であれば、ｍ＝１、２、３、４、５であり、膜厚係数Ｘ_ｍは、それぞれの濃度階調領域の膜厚の段階に対応して、Ｘ_１＝１．０、Ｘ_２＝０．８、Ｘ_３＝０．６、Ｘ_４＝０．４、Ｘ_５＝０．２とする。

膜厚係数Ｘ_ｍを、すでに述べたように、印刷媒体の下地分光反射率Ｒ_０（λ）と、吸収係数Ｋ（λ）と、散乱係数Ｓ（λ）とともに、クベルカ・ムンクの式である（１）式に代入し、それぞれの指令網点面積率の網点が含む濃度階調領域の分光反射率として、濃度階調分光反射率Ｒ_ｉ１（λ）、Ｒ_ｉ２（λ）、Ｒ_ｉ３（λ）、…、Ｒ_ｉｍ（λ）の各々を算出する。そして、濃度階調分光反射率Ｒ_ｉ１（λ）、Ｒ_ｉ２（λ）、Ｒ_ｉ３（λ）、…、Ｒ_ｉｍ（λ）の各々を、（５）式により、濃度階調分光光学濃度ＯＤ_ｉ１（λ）、ＯＤ_ｉ２（λ）、ＯＤ_ｉ３（λ）、…、ＯＤ_ｉｍ（λ）にそれぞれ変換する。この濃度階調分光光学濃度ＯＤ_ｉ１（λ）、ＯＤ_ｉ２（λ）、ＯＤ_ｉ３（λ）、…、ＯＤ_ｉｍ（λ）が、後述する計算モデルで用いる、網点を構成する複数の濃度階調領域の各々の分光光学濃度として用いられる。

濃度階調出現率算出部２０３は、濃度階調分光光学濃度算出部２０２から、濃度階調分光光学濃度ＯＤ_ｉ１（λ）、ＯＤ_ｉ２（λ）、ＯＤ_ｉ３（λ）、…、ＯＤ_ｉｍ（λ）の各々を読み込む。また、濃度階調出現率算出部２０３は、指令網点面積率毎の測定分光反射率Ｒ_ｓ（λ）の各々を、測定分光反射率データベース１０５から読み込む。
そして、濃度階調出現率算出部２０３は、（５）式を用いることにより、測定分光反射率Ｒ_ｓ（λ）の各々を、測定分光光学濃度ＯＤ_ｓ（λ）にそれぞれ変換する。
次に、濃度階調出現率算出部２０３は、以下の（６）式（計算モデル）に対して、濃度階調分光光学濃度ＯＤ_ｉ１（λ）、ＯＤ_ｉ２（λ）、ＯＤ_ｉ３（λ）、…、ＯＤ_ｉｍ（λ）の各々を代入し、後述する処理により算出分光光学濃度ＯＤ’（ｓ，λ）を算出する。

ここで、濃度階調出現率算出部２０３は、以下の（７）式において、出現率（用紙の印刷部分における網点を構成する各濃度階調の濃度階調領域の面積の比率）の数値を変更しつつ、算出分光光学濃度ＯＤ’（ｓ，λ）を算出する。そして、濃度階調出現率算出部２０３は、算出された算出分光光学濃度ＯＤ’（ｓ，λ）の各々と、測定分光光学濃度ＯＤ_ｓ（λ）との平均二乗誤差ＲＭＳＥを、指令網点面積率毎に所定の波長範囲において求める。濃度階調出現率算出部２０３は、算出分光光学濃度ＯＤ’（ｓ，λ）と測定分光光学濃度ＯＤ_ｓ（λ）との平均二乗誤差が最も小さくなる濃度階調領域の各々の出現率を求める。ここで、ｓは指令網点面積率である。

そして、濃度階調出現率算出部２０３は、濃度階調領域の各々の出現率により、各濃度階調領域の出現率関数ａ_１（ｓ）、ａ_２（ｓ）、ａ_３（ｓ）、…、ａ_ｍ（ｓ）を求める。ここで、濃度階調出現率算出部２０３は、濃度階調領域毎に、得られた出現率を用いて、この出現率を指定網点面積率ｓの２次関数などをフィッティングし、出現率関数を求めてもよい。濃度階調出現率算出部２０３は、求めた各濃度階調領域の出現率関数ａ_１（ｓ）、ａ_２（ｓ）、ａ_３（ｓ）、…、ａ_ｍ（ｓ）を、濃度階調出現率テーブルデータベース２０４に書き込んで記憶させる。出現率関数ａ_１（ｓ）、ａ_２（ｓ）、ａ_３（ｓ）、…、ａ_ｍ（ｓ）については、分光反射率予測部１０８における説明と同様に、（６）式の計算モデルで用いる濃度階調領域の出現率を求める関数であり、図２の説明におけるように、グラビア印刷による網点の構造をモデル化するために用いている。

上記（７）式において、濃度階調出現率算出部２０３は、波長λが３８０ｎｍから７３０ｎｍをｎ分割した刻み幅により、算出分光光学濃度ＯＤ’（ｓ，λ）と測定分光光学濃度ＯＤ_ｓ（λ）とにおける各波長λに誤差の二乗を加算した平均二乗誤差ＲＭＳＥを、指令網点面積率毎に求めている。
上述したように、濃度階調出現率算出部２０３は、（６）式の計算モデルを用いて、指令網点面積率によるインキ（原色インキ）が印刷された印刷部分において、網点が含む濃度階調毎に、濃度階調の濃度階調分光光学濃度ＯＤ_ｉｍ（λ）と濃度階調の出現率関数ａ_１（ｓ）、ａ_２（ｓ）、ａ_３（ｓ）、…、ａ_ｍ（ｓ）とを乗算した結果を加算して、指令網点面積率の印刷部分の分光光学濃度を算出する。

また、濃度階調出現率算出部２０３は、求めた濃度階調領域の指令網点面積率における各濃度階調領域の出現率を示す出現率関数ａ_ｍ（ｓ）を、濃度階調出現率テーブルデータベース２０４に対し、濃度階調領域の濃度階調分光光学濃度と対応させて書き込んで記憶させる。同様に、濃度階調出現率算出部２０３は、他の原色インキの各々に対しても、濃度階調領域の指令網点面積率における各濃度階調領域の出現率を示す出現率関数ａ_ｍ（ｓ）を、濃度階調出現率テーブルデータベース２０４に対して、書き込んで記憶させる。

また、原色インキを所定の比率で混ぜ合わせて生成した特色インキの場合、上述した原色インキの濃度階調領域の出現率関数ａ_ｍ（ｓ）を濃度階調出現率テーブルデータベース２０４から読み出して用いる。ここで、特色インキで混ぜ合わせる原色インキのいずれかの出現率関数ａ_ｍ（ｓ）を用いても良いし、また混ぜ合わせる原色インキそれぞれの出現率関数ａ_ｍ（ｓ）を配合率に応じて組合わせて用いても良い。
この特色インキの場合、分光反射率予測部１０８における説明と同様に、混ぜ合わせる原色インキの比率に応じて、すでに説明した（４）式において、特色インキの散乱係数Ｓ_ｔ（λ）及び吸収係数Ｋ_ｔ（λ）を算出する。

上記（４）式において、係数α及び係数βの各々は、原色インキ＃１及び原色インキ＃２それぞれの混ぜ合わせる比率を示す係数である。原色インキ＃１の吸収係数Ｋ_１（λ）に対して係数αを乗算し、原色インキ＃２の吸収係数Ｋ_２に対して係数βを乗算し、加算した数値を、特色インキの吸収係数Ｋ_ｔ（λ）としている。同様に、原色インキ＃１の散乱係数Ｓ_１（λ）に対して係数αを乗算し、原色インキ＃２の散乱係数Ｓ_２（λ）に対して係数βを乗算し、加算した数値を、特色インキの散乱係数Ｓ_ｔ（λ）としている。

＜分光光学濃度予測部１１１の動作＞
次に、分光光学濃度予測部１１１は、拡張ノイゲバウア原色算出部１１１１、拡張ノイゲバウア原色出現率算出部１１１２、分光光学濃度算出部１１１３及び分光反射率算出部１１１４を備えている。
拡張ノイゲバウア原色算出部１１１１は、重ね合わせる色の順番に従い、下地になるインキ（原色インキあるいは特色インキ）と、下地の表面に印刷されるインキとを決定する。また、拡張ノイゲバウア原色算出部１１１１は、重ね合わせるインキの各々の指令網点面積率を、入力部１０１から読み込む。

また、拡張ノイゲバウア原色算出部１１１１は、下地のインキの上に印刷されるインキ（原色インキあるいは特色インキ）の散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）を読み込む。
次に、拡張ノイゲバウア原色算出部１１１１は、すでに図４で説明されたように、（１）式に下地となるインキの濃度階調分光反射率Ｒ_ｉｍ（λ）と、下地のインキの上に印刷されるインキの散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）と、濃度階調領域の膜厚係数Ｘ_ｍとの各々を代入し、下地のインキの上に印刷される網点の濃度階調分光反射率Ｒ_ｉｍ（λ）を算出する。

ここで、拡張ノイゲバウア原色算出部１１１１は、上述した下地のインキの網点上に印刷されるインキの網点の濃度階調分光反射率Ｒ_ｉｍ（λ）の計算を、下地のインキの網点における濃度階調領域と、下地のインキの網点上に印刷されるインキの網点の濃度階調領域との重ね合う部分の組合わせの全てにおいて行う。そして、拡張ノイゲバウア原色算出部１１１１は、得られた下地のインキの網点上に印刷されるインキの濃度階調分光反射率Ｒ_ｉｍ（λ）を、（５）式により、濃度階調分光光学濃度ＯＤ_ｉｍ（λ）に変換する。拡張ノイゲバウア原色算出部１１１１は、求めた濃度階調分光光学濃度ＯＤ_ｉｍ（λ）を一時記憶部１２５のノイゲバウア原色テーブルに書き込んで記憶させる。

図１に戻り、拡張ノイゲバウア原色出現率算出部１１１２は、図６の説明と同様に、下地のインキの濃度階調領域と、下地のインキに重ねて印刷されるインキの濃度階調領域との重なりの出現率の割合を算出する。

また、使用するインキ（原色インキあるいは特色インキ）の重ね合わせに用いる出現率の式は、分光反射率予測部１０８において説明したように、図７に示す出現率のテーブルとして濃度階調出現率テーブルデータベース２０４に書き込んで記憶されている。
拡張ノイゲバウア原色出現率算出部１１１２は、重ね合わせるインキの種類と、重ね合わせるインキの網点を示す指令網点面積率との各々の組み合わせにより、濃度階調出現率テーブルデータベース２０４の図７のテーブルの式を読み込む。
また、拡張ノイゲバウア原色出現率算出部１１１２は、濃度階調出現率テーブルデータベース２０４から、シアンとマジェンダとの各々の濃度階調コア領域１及び濃度階調フリンジ領域２それぞれの出現率α１、α２、β１、β２を読み込む。

また、拡張ノイゲバウア原色算出部１１１１は、シアンが下地で、マジェンダがシアンに重ね合わされる場合、シアンの分光反射率を下地の分光反射率Ｒ_０（λ）として、指令網点面積率で出現するマジェンダの濃度階調領域の膜厚を用い、（１）式により、重なり部分における濃度階調領域の分光反射率Ｒ_ＫＭ（λ）を算出する。例えば、拡張ノイゲバウア原色算出部１１１１は、領域Ｑ５における分光反射率Ｒ_ＫＭ（λ）を求める際、印刷媒体上のシアンのインキの膜厚１００％における濃度階調分光反射率Ｒ_ｉｍ（λ）を下地の分光反射率Ｒ_０（λ）として、マジェンダのインキの膜厚１００％を下地に重ねた際の重なり部分における印刷部分の分光反射率Ｒ_ＫＭ（λ）を算出する。そして、拡張ノイゲバウア原色算出部１１１１は、算出した分光反射率Ｒ_ＫＭ（λ）を、（５）式により、分光濃度階調光学濃度ＯＤ_ＫＭ（λ）に変換する。

同様に、拡張ノイゲバウア原色算出部１１１１は、領域Ｑ７における分光反射率Ｒ_ＫＭ（λ）を求める際、印刷媒体上のシアンのインキの膜厚５０％における濃度階調分光反射率Ｒ_ｉｍ（λ）を下地の分光反射率Ｒ_０（λ）とし、マジェンダのインキの膜厚１００％を重ねた際の分光反射率Ｒ_ＫＭ（λ）を算出する。そして、拡張ノイゲバウア原色算出部１１１１は、算出した分光反射率Ｒ_ＫＭ（λ）を、（５）式により、分光光学濃度ＯＤ_ＫＭ（λ）に変換する。同様に、拡張ノイゲバウア原色の全ての組み合わせについて分光光学濃度ＯＤ_ＫＭ（λ）を算出する。

そして、分光光学濃度算出部１１１３は、上述した領域Ｑ１から領域Ｑ９の拡張ノイゲバウア原色の各々の分光光学濃度に拡張ノイゲバウア原色出現率を乗算し、波長毎に加算する。
分光光学濃度算出部１１１３は、シアンのインキの網点が印刷された印刷媒体に対し、マジェンダのインキの網点を重ねて印刷した際における印刷部分の予測分光光学濃度ＯＤ_Ｄ（λ）を算出する。

分光反射率算出部１１１４は、光学濃度算出部１１１３が算出した予測分光光学濃度ＯＤ_Ｄ（λ）を、以下の（８）式を用い、第２予測分光反射率Ｒ_Ｄ２（λ）に変換する。

＜混合予測部１１４の動作＞
混合予測部１１４は、予測パラメータデータベース１２６から、第１予測分光反射率Ｒ_Ｄ１（λ）及び第２予測分光反射率Ｒ_Ｄ２（λ）に乗算する重み係数ｗを読み出す。例えば、混合予測部１１４は、ｗＲ_Ｄ１（λ）と（１−ｗ）Ｒ_Ｄ２（λ）とを加算して、統合予測分光率Ｒ_Ｄ（λ）を求める。
色予測テーブル作成部１０９は、この統合予測分光反射率Ｒ_Ｄ（λ）に基づき、観測光源を定めて、三刺激値ＸＹＺやＣＩＥＬＡＢ値等を算出し、再現色の予測を行う色予測テーブルを作成する。つまり、この色予測テーブル作成部１０９は、色予測処理の対象である入力データ（本実施例ではＣＭＹ値）の情報と、統合予測分光反射率Ｒ_Ｄ（λ）に基づき予測される再現色を表現する再現色情報とを対応付ける色予測テーブルを作成する。なお、本実施例においては、色予測テーブル作成部１０９がこの色予測テーブルを公知のＩＣＣ（International Color Consortium）プロファイルフォーマットで作成し、出力部１１０を介して出力する。

＜色予測テーブルの生成＞
図８は、本実施形態における色予測テーブル作成部１０９の色予測テーブルの生成処理の動作例を説明するフローチャートである。

ステップＳ４０１：
拡張ノイゲバウア原色出現率算出部１０８２は、一時記憶部１２５のノイゲバウア原色テーブルから、重なり合う領域の特色インキのノイゲバウア原色の分光反射率Ｒ_ＫＭ（λ）を読み込む。
同様に、拡張ノイゲバウア原色出現率算出部１１１２は、一時記憶部１２５のノイゲバウア原色テーブルから、重なり合う領域の特色インキのノイゲバウア原色の分光光学濃度ＯＤ_ＫＭ（λ）を読み込む。

ステップＳ４０２：
混合予測部１１４は、分光反射率予測部１０８が予測した第１予測分光反射率Ｒ_Ｄ１（λ）と、分光光学濃度予測部１１１が予測した第２予測分光反射率Ｒ_Ｄ２（λ）を統合する際、第１予測分光反射率Ｒ_Ｄ１（λ）及び第２予測分光反射率Ｒ_Ｄ２（λ）の各々に乗算する重み係数ｗを、予測パラメータデータベース１２６から読み出す。

ステップＳ４０３：
次に、拡張ノイゲバウア原色出現率算出部１０８２は、一時記憶部１２５における色予測テーブルにおける重ね合わせるインキの指令網点面積率のマトリクスにおいて、未計算の指令網点面積率を抽出する。
そして、拡張ノイゲバウア原色出現率算出部１０８２は、抽出した指令網点面積率の組み合わせを、計算対象として設定する。上述した処理は、拡張ノイゲバウア原色出現率算出部１１１２が行う構成としても良い。

ステップＳ４０４：
分光反射率予測部１０８は、上記指令網点面積率の組み合わせにより、第１予測分光反射率Ｒ_Ｄ１（λ）を算出する。
分光光学濃度予測部１１１は、上記指令網点面積率の組み合わせにより、第２予測分光反射率Ｒ_Ｄ２（λ）を算出する。

また、混合予測部１１４は、第１予測分光反射率Ｒ_Ｄ１（λ）及び第２予測分光反射率Ｒ_Ｄ２（λ）にステップＳ４０２で読み出した重み係数ｗと（１−ｗ）を積算し、ｗＲ_Ｄ１（λ）と（１−ｗ）Ｒ_Ｄ２（λ）とを加算して、統合予測分光反射率Ｒ_Ｄ（λ）を求める。さらに、算出した統合予測分光反射率Ｒ_Ｄ（λ）に対して、観察光源の分光分布と標準観測者を設定することにより、測色値を算出する。

ステップＳ４０５：
そして、混合予測部１１４は、算出した測色値を、重ね合わせるインキ（特色インキあるいは原色インキ）の指令網点面積率の組み合わせに対応させ、一時記憶部１２５の色予測テーブルに書き込んで記憶させる。

ステップＳ４０６：
拡張ノイゲバウア原色出現率算出部１０８２は、一時記憶部１２５の色予測テーブルにおける全ての指令網点面積率の組み合わせの統合予測分光反射率Ｒ_Ｄ（λ）の算出が終了したか否かの判定を行う。
そして、拡張ノイゲバウア原色出現率算出部１０８２は、全ての指令網点面積率の組み合わせの統合予測分光率Ｒ_Ｄ（λ）の算出が終了した場合、処理をステップＳ４０７へ進める。一方、拡張ノイゲバウア原色出現率算出部１０８２は、全ての指令網点面積率の組み合わせの統合予測分光反射率Ｒ_Ｄ（λ）の算出が終了していない場合、処理をステップＳ４０３へ進める。上述した処理は、拡張ノイゲバウア原色出現率算出部１１１２が行う構成としても良い。

ステップＳ４０７：
色予測テーブス作成部１０９は、一時記憶部１２５の色予測テーブルのデータを読み出し、公知のＩＣＣプロファイルフォーマットに変換し、変換後の色予測プロファイルのデータを出力部１１０を介して出力する。

＜特色分解テーブルの生成＞
特色分解テーブル作成部１３０は、特色分解テーブルを作成する。特色分解テーブルとは、重ね刷りする特色インキの組み合わせ毎に生成され、特色インキを組合わせて印刷される色空間（本実施形態においてはＬ*ａ*ｂ*の色空間）における各座標に対して、その座標が表す測色値を再現するために必要な特色インキ毎の指令網点面積率の組合わせが記憶されているテーブルである。この座標点は、色空間を所定の周期で分割して離散的に設定されている。
すなわち、特色分解テーブル作成部１３０は、ＸＹＺやＣＩＥＬＡＢ値等の入力値と、これを再現する特色インキの指令網点面積率の出力値とを対応付ける特色分解テーブルを作成し、一時記憶部１２５の特色分解テーブルに書き込んで記憶させる。また、本実施例においては、特色分解テーブル作成部１３０は、この一時記憶部１２５の特色分解テーブルを公知のＩＣＣプロファイルフォーマットに変換し、変換後の特色分解プロファイルのデータを出力部１１０を介して出力する。

図９は、本実施形態における特色分解テーブル作成部１３０の特色分解テーブルの生成処理の動作例を説明するフローチャートである。

ステップＳ５０１：
特色分解テーブル作成部１３０は、一時記憶部１２５のノイゲバウア原色テーブルから、重なり合う領域の特色インキのノイゲバウア原色の分光反射率Ｒ_ＫＭ（λ）を読み込む。

ステップＳ５０２：
特色分解テーブル作成部１３０は、予測パラメータデータベース１２６から、分光反射率予測部１０８が予測した第１予測分光反射率Ｒ_Ｄ１（λ）と、分光光学濃度予測部１１１が予測した第２予測分光反射率Ｒ_Ｄ２（λ）を統合する際、第１予測分光反射率Ｒ_Ｄ１（λ）及び第２予測分光反射率Ｒ_Ｄ２（λ）の各々に乗算する重み係数ｗを読み出す。

ステップＳ５０３：
特色分解テーブル作成部１３０は、組合わせる特色インキ各々の指令網点面積率の初期値が予めセットとして記憶されている網点面積率初期値セットデータベース１２７から、組合わせる特色インキ各々の指令網点面積率の初期設定値のセットを読み出す。
そして、特色分解テーブル作成部１３０は、読み出した初期設定値を以下の計算を行う指令網点面積率の初期値として設定する。

ステップＳ５０４：
特色分解テーブル作成部１３０は、指令網点面積率を修正（調整）する際における指令網点面積率の変化量の許容範囲を設定する。

ステップＳ５０５：
特色分解テーブル作成部１３０は、特色分解テーブルの色空間において、未計算である座標のＬ*ａ*ｂ*値を目標Ｌ*ａ*ｂ*値（目的のＬ*ａ*ｂ*値）として選択する。このとき、色分解テーブル作成部１３０は、修正回数が「０」の場合、読み出した初期設定値を指令網点面積率として用いる。

ステップＳ５０６：
特色分解テーブル作成部１３０は、指令網点面積率に対して、分光反射率予測部１０８と分光光学濃度予測部１１１との各々と同様の処理により、第１予測分光反射率Ｒ_Ｄ１（λ）及び第２予測分光反射率Ｒ_Ｄ２（λ）の各々を生成する。

そして、特色分解テーブル作成部１３０は、第１予測分光反射率Ｒ_Ｄ１（λ）及び第２予測分光反射率Ｒ_Ｄ２（λ）の各々に、読み出した重み係数ｗ、（１−ｗ）それぞれを乗算する。特色分解テーブル作成部１３０は、乗算結果のｗＲ_Ｄ１（λ）と（１−ｗ）Ｒ_Ｄ２（λ）とを加算して、統合予測分光率Ｒ_Ｄ（λ）を求める。さらに、算出した統合予測分光反射率Ｒ_Ｄ（λ）に対して、観察光源の分光分布と標準観測者を設定することにより、測色値を算出する。

ステップＳ５０７：
次に、特色分解テーブル作成部１３０は、目標Ｌ*ａ*ｂ*値と統合予測分光率Ｒ_Ｄ（λ）から算出したＬ*ａ*ｂ*値との色差を算出する。

ステップＳ５０８：
特色分解テーブル作成部１３０は、目標Ｌ*ａ*ｂ*値と算出されたＬ*ａ*ｂ*値との色差が予め設定された範囲内に含まれているか否かの判定を行う。
このとき、特色分解テーブル作成部１３０は、色差が予め設定された範囲内に含まれている場合、処理をステップＳ５０９へ進め、一方、色差が予め設定された範囲内に含まれていない場合、処理をステップＳ５１４へ進める。

ステップＳ５０９：
特色分解テーブル作成部１３０は、指令網点面積率の修正回数のカウンタをリセットし、修正回数を「０」とする。

ステップＳ５１０：
特色分解テーブル作成部１３０は、現在の統合予測分光率Ｒ_Ｄ（λ）を算出した際の指令網点面積率を、一時記憶部１２５の特色分解テーブルにおける目標Ｌ*ａ*ｂ*値に対応させて書き込んで記憶させる。

ステップＳ５１１：
特色分解テーブル作成部１３０は、特色分解テーブルの色空間において、所定周期の座標値におけるＬ*ａ*ｂ*値の全てを算出したか否かの判定を行う。
このとき、特色分解テーブル作成部１３０は、色空間におけるＬ*ａ*ｂ*値の全てを算出した場合、処理をステップＳ５１２へ進め、一方、色空間におけるＬ*ａ*ｂ*値の全てを算出していない場合、処理をステップＳ５０５へ進める。

ステップＳ５１２：
特色分解テーブル作成部１３０は、特色分解テーブルの作成ループ回数が規定内であるか否かの判定を行う。
このとき、特色分解テーブル作成部１３０は、特色分解テーブルの作成ループ回数が規定内である場合、処理をステップＳ５１６へ進め、一方、特色分解テーブルの作成ループ回数が規定内でない場合、処理をステップＳ５１３へ進める。

ステップＳ５１３：
特色分解テーブル作成部１３０は、一時記憶部１２５の特色分解テーブルのデータを読み出し、公知のＩＣＣプロファイルフォーマットに変換し、変換後の特色分解プロファイルのデータを出力部１１０を介して出力する。

ステップＳ５１４：
特色分解テーブル作成部１３０は、指令網点面積率の修正回数が予め設定された規定回数であるか否かの判定を行う。
特色分解テーブル作成部１３０は、指令網点面積率の修正回数が規定回数内である場合、処理をステップＳ５１５へ進め、指令網点面積率の修正回数が規定回数を超えた場合、処理をステップＳ５０９へ進める。

ステップＳ５１５：
特色分解テーブル作成部１３０は、特色インキ各々の指令網点面積率を修正許容範囲における変化量により修正する。
このとき、特色分解テーブル作成部１３０は、指令網点面積率の修正回数を計数するカウンタをインクリメント（「１」を加算）する。

ステップＳ５１６：
特色分解テーブル作成部１３０は、特色分解テーブルの色空間において、全てのＬ*ａ*ｂ*値に対応する指令網点面積率を特色インキ毎にそれぞれ平滑化する。平滑化の方法としては、隣接する座標値に対して設定されている指令網点面積率を、特色インキ毎に平均化した値を用いたり、あるいは、所定周期で分割されたＬ*ａ*ｂ*値の三次元空間に対して、一般的な平滑化フィルタを適用するなど、どのような方法を用いても良い。
そして、特色分解テーブル作成部１３０は、指令網点面積率の初期設定値のセットに代え、平滑化した指令網点面積率の数値のセットを新たな初期値とする。
このとき、特色分解テーブル作成部１３０は、テーブル作成ループ回数を計数するカウンタをインクリメントする。

ステップＳ５１７：
特色分解テーブル作成部１３０は、指令網点面積率の変化量の許容範囲を現在の数値に対し、所定の割合（例えば１／２）を乗じて、より小さな許容範囲に変更（再定義）する。さらに、色分解テーブルの色空間において設定した全ての座標を、未計算の状態に戻す。
そして、特色分解テーブル作成部１３０は、処理をステップＳ５０５に進める。

上述したように、本実施形態によれば、原色インキから特色インキの吸収特性及び散乱特性を推定し、特色インキの分光反射率を算出し、特色インキの重ね刷りされた際の重なり領域の分光反射率を求めることができる。
このため、本実施形態によれば、面積変調階調表現と、濃度変調階調表現との双方の色表現を有するグラビア印刷などの印刷において特色インキを重ね刷りした印刷物における再現色の色予測を高い精度で容易に行うことができる。

また、本実施形態によれば、拡張ノイゲバウア原色の分光反射率と、その出現率から算出した予測分光反射率Ｒ_Ｄ１（λ）と、拡張ノイゲバウア原色の分光光学濃度と、その出現率から算出した予測分光反射率Ｒ_Ｄ２（λ）とを、予め求めた重み係数ｗにより混合するため、より実測値の分光反射率に近い予測分光反射率として混合予測分光反射率Ｒ_Ｗ（λ）を求めることができる。

以下、本発明の第２の実施形態について、図面を参照して説明する。
図１０は、本発明の第２の実施形態による色予測システム２の構成例を示すブロック図である。図１０の構成において、図１における第１の実施形態と同様の構成には同一の符号を付してある。図１０において、色予測システム２は、入力部１０１、濃度階調出現率テーブルデータベース１０４、吸収係数・散乱係数データベース１０６、色予測テーブル作成部１０９、出力部１１０、特色分解部１２４、一時記憶部１２５、網点面積率初期値セットデータベース１２７、近似色データベース１２８、特色分解テーブル作成部１３０、ノイゲバウア原色分光反射率予測部１５０、特色インキ任意階調分光反射率算出部２００、網点掛け合わせ分光反射率計算部２０４、印刷色予測部２０５を備えている。

ノイゲバウア原色分光反射率算出部１５０は、吸収係数・散乱係数データベース１０６に記憶されている原色インキの散乱係数Ｓ（λ）及び吸収係数Ｋ（λ）から、特色インキの吸収・散乱特性を算出し、これを基にノイゲバウア原色の分光反射率を算出する。
特色インキ任意階調分光反射率算出部２００は、原色インキの階調特性テーブルデータベース２０１から、特色インキの階調特性を算出し、これを基に特色インキの任意の階調の分光反射率を算出する。
以下、ノイゲバウア原色分光反射率算出部１５０及び特色インキ任意階調分光反射率算出部２００における各部の動作を説明する。

ノイゲバウア原色分光反射率算出部１５０は、吸収係数・散乱係数算出部１０７、特色インキ配合決定部１５１、ノイゲバウア原色分光反射率予測部１５２を備えている。また、特色インキ任意階調分光反射率算出部２００は、特色インキ階調特性パラメータ決定部２０３、特色インキ階調分光反射率予測部２０３を備えている。
特色インキ配合決定部１５１は、図３のフローチャートに従い、見本色の測色値が与えられた特色に対して、その特色を構成する原色の組み合わせとその配合比を決定する。

すなわち、特色インキ配合決定部１５１は、第１の実施形態における特色インキ分光反射率算出部１２１及び特色インキ配合比決定部１２２の双方と同様の動作を行う。
吸収係数・散乱係数算出部１０７は、特色インキ配合決定部１５１が決定した特色インキにおける原色インキの配合比から、第１の実施形態と同様に、特色インキの散乱係数Ｓ_ｔ（λ）及び吸収係数Ｋ_ｔ（λ）を算出する。

ノイゲバウア原色分光反射率予測部１５２は、特色インキの散乱係数Ｓ_ｔ（λ）及び吸収係数Ｋ_ｔ（λ）、印刷媒体の下地分光反射率Ｒ_０（λ）、膜厚係数Ｘ_ｍを（１）式に代入することにより、特色ベタの分光反射率を算出する。また、ノイゲバウア原色分光反射率予測部１５２は、複数の使用する特色（原色を含んでも良い）の組み合わせと、その刷り順が指定された際、全ての使用する特色の組み合わせに対して、その重ね合わせの分光反射率（すなわちノイゲバウア原色の分光反射率）を（１）式により算出する。算出したノイゲバウア原色の分光反射率は、一時記憶部１２５のノイゲバウア原色テーブルに書き出して記憶する。

階調特性テーブルデータベース２０１には、原色インキの階調特性を表すパラメータ（後述する階調特性パラメータ）があらかじめ計算されて記憶されている。本実施形態では、単色の階調特性を（９）式で示すモデルで表すものとし、（９）式のａ_１（s）、ａ_２（s）を階調特性パラメータとする。原色インキの階調特性パラメータは、第１の実施形態と同様、指令網点面積率毎の測定分光反射率Ｒ_ｓ（λ）と、（９）式により算出される算出分光反射率の誤差が最小となるよう、事前に求めておくことができる。

特色インキ階調特性パラメータ決定部２０２は、第１の実施形態における特色濃度階調出現率算出部１２３と同様に、上記特色インキ配合決定部１５１が決定した特色の配合比に基づいて、階調特性テーブルデータベース２０１から原色の階調特性を読み出し、特色の階調特性を決定する。第１の実施形態と同様に、読み出したいずれかの原色の階調特性を特色の階調特性としても良いし、配合比に基づいてそれぞれの原色の階調特性を組み合わせて用いても良い。また、特色インキ階調分光反射率予測部２０３は、ノイゲバウア原色分光反射率予測部１５２が算出した特色ベタの分光反射率を読み出し、求めた特色の階調特性とともに（９）式に代入することにより、特色インキ毎に、各指令網点面積率に対応する分光反射率Ｒ’（ｓ，λ）を求め、一時記憶部１２５の特色インキ階調分光反射率テーブルに記録する。

図１１は、特色インキの階調毎の分光反射率を算出する処理を示すフローチャートである。
ステップＳ３０１：
特色インキ階調分光反射率予測部２０３は、階調毎の分光反射率を算出する特色インキのベタの分光反射率を、一時記憶部１２５のノイゲバウア原色テーブルから読み込む。

ステップＳ３０２：
特色インキ階調分光反射率予測部２０３は、一時記憶部１２５の特色インキの配合比テーブルから、上記設定した特色インキの配合比を読み出す。

ステップＳ３０３：
特色インキ階調特性パラメータ決定部２０２は、特色インキを構成する原色インキの階調特性パラメータａ_１（s）、ａ_２（s）を階調特性テーブルデータベース２０１から読み出す。

ステップＳ３０４：
特色インキ濃度階調特性パラメータ決定部２０２は、階調特性テーブルデータベース２０１から読み出した原色インキの階調特性パラメータａ_１（s）、ａ_２（s）を基にし、すでに述べた処理によって特色インキの階調特性パラメータａ_１（s）、ａ_２（s）を求める。

ステップＳ３０５：
特色インキ階調分光反射率予測部２０３は、原色インキを所定の比率で混ぜ合わせて生成した特色インキの指令網点面積率に対応する分光反射率Ｒ’（ｓ，λ）を、特色インキ階調特性パラメータ決定部２０２が求めた階調特性パラメータａ_１（s）、ａ_２（s）に基づき、上記（９）式により求める。

ステップＳ３０６：
そして、特色インキ階調分光反射率予測部２０３は、求めた分光反射率Ｒ’（ｓ，λ）を一時記憶部１２５の特色階調分光反射率テーブルに対して書き込んで記憶させる。

網点掛け合わせ分光反射率計算部２０４は、特色インキの指令網点面積率に対応する分光反射率Ｒ’（ｓ，λ）を、一時記憶部１２５の特色階調分光反射率テーブルから読み出す。
そして、網点掛け合わせ分光反射率計算部２０４は、分光反射率Ｒ’（ｓ，λ）、特色インキベタの分光反射率Ｒ_ＫＭ（λ）、印刷下地の紙の分光反射率Ｒ₀（λ）を、以下の（１０）式に対して代入して、分光実効網点面積率ａ_{ｅｆｆ，ｉ}（λ）を求める。

網点掛け合わせ分光反射率計算部２０４は、求めた実効網点面積率ａ_{ｅｆｆ，ｉ}（λ）から所定のノイゲバウア式により、ノイゲバウア原色ｌの分光実効網点面積率Ｆ_ａ，ｌ（λ）を求める。ここで、ｌはノイゲバウア原色を示している。
また、網点掛け合わせ分光反射率計算部２０４は、ノイゲバウア原色ｌの分光実効網点面積率Ｆ_ａ，ｌ（λ）及びノイゲバウア原色の分光反射率Ｒ_ｌ（λ）を、以下の（１１）式に代入して、再現色の分光反射率Ｒ_j（λ）を算出する。

印刷色予測部２０５は、この再現色の分光反射率Ｒ_j（λ）に対して、観察光源の分光分布と標準観測者を設定することにより、測色値を求める。
そして、印刷色予測部２０５は、ノイゲバウア原色分光反射率予測部１５２の算出した測色値を、一時記憶部１２５の色予測テーブルに対し、重ね合わせる特色インキの指令網点面積率の組み合わせに対応させて書き込み、記憶させる。

色予測テーブル作成部１０９は、第１の実施形態と同様に、一時記憶部１２５の予測テーブルにおける指令網点面積率の未計算の組み合わせの測色値を網点掛け合わせ分光反射率計算部２０４および印刷色予測部２０５に対して算出させる。そして、色予測テーブル作成部１０９は、所定の特色インキの組み合わせに対応した色予測テーブルにおける指令網点面積率の掛け合わせの全てに対して測色値を求める。
また、出力部１１０、特色分解部１２４及び特色分解テーブル作成部１３０の各々は、第１の実施形態と同様のため、説明を省略する。

第２の実施形態では、構成を簡便にするために、ノイゲバウア原色の分光反射率から求めた予測分光反射率を、印刷色予測に用いているが、第１の実施形態のようにノイゲバウア原色の分光光学濃度を算出し、分光光学濃度から求めた第２の予測分光反射率を求める構成を追加し、分光反射率から求めた第１の予測分光反射率と、分光光学濃度から求めた第２の予測分光反射率から、混合モデルにより統合予測分光反射率を求め、この統合予測分光反射率から印刷色の測色値を算出する構成としても良い。

上述したように、本実施形態によれば、オフセット印刷等のような面積変調階調表現の印刷において、原色インキから特色インキの吸収特性及び散乱特性を推定し、特色インキの分光反射率を算出し、特色インキの重ね刷りされた際の重なり領域の分光反射率を求めることができる。
また、本実施形態によれば、ノイゲバウア原色の分光反射率と特色インキの任意階調の分光反射率から、分光拡張ノイゲバウア混合モデルにより、網点の掛け合せにより再現されるあらゆる色の測色値や分光反射率を算出可能である。

すなわち、本実施形態によれば、原色インキの特性から特色インキの特性を求め、これを基に特色インキで再現されるあらゆる色の測色値や分光反射率を高精度に求めるという本発明に関わる基本的な機能はいずれの実施形態においても同様である。また、本実施形態は、上述したように、主にオフセット印刷など、網点の面積により階調を表現する面積変調階調表現の印刷方式に適用されると高い予測精度を得ることが可能である。
また、本実施形態に基づけば、伝統的なノイゲバウア方程式の発展により展開されるいずれの色予測方法においても容易に適用可能な特色インキの色予測方法を提供することが可能となる。

なお、本発明における図１の色予測システム１及び図１０の色予測システム２の各々の機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより色予測の制御を行ってもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。
また、「コンピュータシステム」は、ホームページ提供環境（あるいは表示環境）を備えたＷＷＷシステムも含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムが送信された場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリ（ＲＡＭ）のように、一定時間プログラムを保持しているものも含むものとする。

また、上記プログラムは、このプログラムを記憶装置等に格納したコンピュータシステムから、伝送媒体を介して、あるいは、伝送媒体中の伝送波により他のコンピュータシステムに伝送されてもよい。ここで、プログラムを伝送する「伝送媒体」は、インターネット等のネットワーク（通信網）や電話回線等の通信回線（通信線）のように情報を伝送する機能を有する媒体のことをいう。また、上記プログラムは、前述した機能の一部を実現するためのものであっても良い。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であっても良い。

１，２…色予測システム
１０１…入力部
１０２，２０２…濃度階調分光反率算出部
１０３，２０３…濃度階調出現率算出部
１０４，２０４…濃度階調出現率テーブルデータベース
１０６…吸収係数・散乱係数データベース
１０７…吸収係数・散乱係数算出部
１０８…分光反射率予測部
１０９…色予測テーブル作成部
１１０…出力部
１１１…分光光学濃度予測部
１１４…混合予測部
１２１…特色インキ分光反射率算出部
１２２…特色インキ配合比決定部
１２３…特色インキ濃度階調出現率算出部
１２４…特色分解部
１２５…一時記憶部
１２６…予測パラメータデータベース
１２７…網点面積率初期値セットデータベース
１２８…近似色データベース
１３０…特色分解テーブル作成部
１５０…ノイゲバウア原色分光反射率算出部
１５１…特色インキ配合決定部
１５２…ノイゲバウア原色分光反射率予測部
２００…特色インキ任意階調分光反射率算出部
２０２…特色インキ階調特性パラメータ決定部
２０３…特色インキ階調分光反射率予測部
１０８１，１１１１…拡張ノイゲバウア原色算出部
１０８２，１１１２…拡張ノイゲバウア原色出現率算出部
１０８３，１１１４…分光反射率算出部
１１１３…分光光学濃度算出部

Claims (8)

1. 原色インキ毎の吸収係数及び散乱係数が予め書き込まれて記憶されている吸収係数・散乱係数データベースと、
   色を再現したい色見本から測定した測色値を再現するのに好適な原色インキの組合せで構成される特色インキを選択し、当該特色インキに含まれる原色インキの吸収係数及び散乱係数を前記吸収係数・散乱係数データベースから読み出し、指定の配合比に応じて特色インキの吸収係数及び散乱係数を求め、求めた吸収係数及び散乱係数から再現したい特色インキの分光反射率を算出する特色インキ分光反射率算出部と、
   特色インキ分光反射率算出部が算出した分光反射率から再現色を求める特色インキ再現色算出部と、
   前記色見本から測定した測色値との差分を求め、適宜前記指定の配合比を修正して当該差分が所定の許容範囲以内となる配合比を算出する特色インキ配合比決定部と、
   前記特色インキを重ね刷りした際のノイゲバウア原色を、前記特色インキ配合比決定部が決定した特色の配合比に基づいて算出した特色インキ各々の吸収係数及び散乱係数から算出する分光データ算出部と
   を備えることを特徴とする色予測システム。
2. 原色インキ毎に指令網点面積率に応じた濃度階調領域の出現率を示す出現率関数が記憶された濃度階調出現率テーブルデータベースと、
   前記特色インキに含まれる前記原色インキの前記出現率関数から、前記原色インキの配合情報に基づいて前記特色インキの出現率関数を求める特色濃度階調出現率算出部と、
   特色インキの濃度階調領域の分光反射率を求める特色濃度階調分光反射率算出部と
   を備え、
   前記分光データ算出部が、前記特色インキの濃度階調領域及び当該濃度階調領域の各々の出現率と、分光反射率との対応関係を示す計算モデルにより、前記指令網点面積率により形成された網点による前記特色の再現色を示す分光反射率を算出する
   ことを特徴とする請求項１に記載の色予測システム。
3. 複数の特色インキの網点から形成される再現色を推定する場合、前記分光データ算出部により求めた前記分光反射率や前記出現率関数に基づいて各々の網点における重なり合う領域の組合わせとその出現率を求め、網点が他の網点の上部に印刷された分光反射率を、推定する
   ことを特徴とする請求項２に記載の色予測システム。
4. 前記分光データ算出部が前記計算モデルから求めた前記特色インキの分光反射率と、前記分光データ算出部が分光光学濃度を用いて求めた分光反射率とを混合する混合比を含む色予測パラメータが記憶されている予測パラメータデータベースをさらに有する
   ことを特徴とする請求項２または請求項３に記載の色予測システム。
5. 前記分光データ算出部に対し、重ね合わせて印刷する前記特色インキ毎の指令網点面積率の組合せ毎に分光反射率を算出させ、色空間における再現色のプロファイルを作成する色予測テーブル作成部をさらに有する
   ことを特徴とする請求項１から請求項４のいずれか一項に記載の色予測システム。
6. 色空間における連続的な座標値毎に、当該座標値の再現色を再現する特色インキ各々の指令網点面積率を求める特色分解テーブル作成部をさらに有する
   ことを特徴とする請求項１から請求項５のいずれか一項に記載の色予測システム。
7. 特色インキを印刷する印刷媒体、及び重ね刷りにおける前記特色インキの印刷順、及び特色インキのインキ膜厚を任意に設定する
   ことを特徴とする請求項１から請求項６のいずれか一項に記載の色予測システム。
8. 特色インキ分光反射率算出部が、色を再現したい色見本から測定した測色値を再現するのに好適な原色インキの組合せで構成される特色インキを選択し、当該特色インキに含まれる原色インキの吸収係数及び散乱係数を、原色インキ毎の吸収係数及び散乱係数が予め書き込まれて記憶されている吸収係数・散乱係数データベースから読み出し、指定の配合比に応じて特色インキの吸収係数及び散乱係数を求め、求めた吸収係数及び散乱係数から再現したい特色インキの分光反射率を算出する特色インキ分光反射率算出過程と、
   特色インキ再現色算出部が、特色インキ分光反射率算出部が算出した分光反射率から再現色を求める特色インキ再現色算出過程と、
   特色インキ配合比決定部が、前記色見本から測定した測色値との差分を求め、適宜前記指定の配合比を修正して当該差分が所定の許容範囲以内となる配合比を算出する特色インキ配合比決定過程と、
   分光データ算出部が、前記特色インキを重ね刷りした際のノイゲバウア原色を、前記特色インキ配合比決定部が決定した特色の配合比に基づいて算出した特色インキ各々の吸収係数及び散乱係数から算出する分光データ算出過程と
   を含むことを特徴とする色予測方法。

Images