(57)【要約】【課題】印刷物の分光反射率を精度良く予測する方法及び装置を提供する。【解決手段】面積変調で色を再現する印刷物の任意の掛け合わせに対する分光反射率を予測し、再現色を予測する方法であって、用紙、各色ベタ、ベタの掛け合わせについて分光反射率を測定する第１の工程と、一次色を任意の網点面積率で変化させた場合の分光反射率を測定する第２の工程と、ユールニールセン修正ノイゲバウア方程式の係数ｎを決定する第３の工程と、実効網点面積率を波長の関数として求める第４の工程と、各一次色の設定網点面積率と実効網点面積率の関係を求める第５の工程と、任意の設定網点面積率をそれぞれ、各一次色の実効網点面積率に変換する第６の工程と、ユールニールセン修正ノイゲバウア方程式に、前記第３の工程で取得した係数ｎと前記第６の工程で取得した実効網点面積率を代入する第７の工程によって、予測分光反射率を計算する印刷物分光反射率予測方法。【選択図】図１
2014143522.tif 000002
【特許請求の範囲】
【請求項１】
面積変調で色を再現する印刷物の任意の掛け合わせに対する分光反射率を予測し、再現色を予測する方法であって、
用紙、各色ベタ、ベタの掛け合わせについて分光反射率（この場合の分光反射率を分光反射率その１と称す）を測定する第１の工程と、
一次色を任意の網点面積率で変化させた場合の分光反射率（この場合の分光反射率を分光反射率その２と称す）を測定する第２の工程と、
第１の工程で測定した分光反射率その１と、第２の工程で取得した分光反射率その２から、最適化によりユールニールセン修正ノイゲバウア方程式の係数ｎを決定する第３の工程と、
前記第３の工程で求めた係数ｎを用いて、第２の工程で測定した分光反射率その２に対して、実効網点面積率を波長の関数として求める第４の工程と、
更に前記第４の工程を全ての分光反射率その２に対して繰り返すことにより、各一次色の設定網点面積率と実効網点面積率の関係を求める第５の工程と、
任意の設定網点面積率の掛け合わせに対する分光反射率を、前記第４の工程で求めた関係から、任意の設定網点面積率をそれぞれ、各一次色の実効網点面積率に変換する第６の工程と、
ユールニールセン修正ノイゲバウア方程式に、前記第３の工程で取得した係数ｎと前記第６の工程で取得した実効網点面積率を代入する第７の工程によって、
予測分光反射率を計算することを特徴とする印刷物分光反射率予測方法。
【請求項２】
前記ユールニールセン修正ノイゲバウア方程式として式（５）、式（６）を用い、前記単色の実効網点面積率を波長の関数として式（７）を用い、更に前記予測分光反射率を計算する際には式（８）、式（９）を用いることを特徴とする請求項１に記載の印刷物分光反射率予測方法。
【数１６】
2014143522.tif 000018
【数１７】
2014143522.tif 000019
【数１８】
2014143522.tif 000020
【数１９】
2014143522.tif 000021
【数２０】
2014143522.tif 000022
【請求項３】
面積変調で色を再現する印刷物の任意の掛け合わせに対する分光反射率を予測し、再現色を予測する印刷物分光反射率予測装置であって、
第１分光反射率受付処理部と、第２分光反射率受付処理部と、最適化計算処理部と、面積率変換関数保持部と、設定網点面積率受付処理部と、網点面積率変換処理部と、分光反射率算出処理部と、を備え、
第１分光反射率受付処理部は、用紙、各色ベタ、ベタの掛け合わせについて分光反射率（この場合の分光反射率を分光反射率その１と称す）を取得し記憶する手段で、
第２分光反射率受付処理部は、一次色の任意の面積率で変化させた分光反射率（この場合の分光反射率を分光反射率その２と称す）を取得し記憶する手段で、
最適化計算処理部は、第１分光反射率受付処理部が取得し記憶した分光反射率その１と、第２分光反射率受付処理部が取得し記憶した分光反射率その２と、を読み出し、最適化によりユールニールセン修正ノイゲバウア方程式の係数ｎを決定する手段で、
面積率変換関数保持部は、最適化計算処理部が計算したｎを読み出して、第２分光反射率受付処理部が取得した分光反射率その２に対して、実効網点面積率を波長の関数として求め、この処理を第２分光反射率受付処理部が取得した全ての分光反射率その２に対して繰り返すことにより、各一次色の設定網点面積率と実効網点面積率の関係を算出して保持する手段で、
設定網点面積率受付処理部は、分光反射率を予測する対象として任意の設定網点面積率
の掛け合わせを入力として取得する手段で、
網点面積率変換処理部は、設定網点面積率受付処理部が受け取った設定網点面積率を、面積率変換関数保持部が保持する関係を用いて各一次色の実効網点面積率に変換する手段で、
分光反射率算出処理部によって、ユールニールセン修正ノイゲバウア方程式に、最適化計算処理部が算出したｎ値と、網点面積率変換処理部が変換した実効網点面積率を代入し、印刷物の予測分光反射率を計算することを特徴とする印刷物分光反射率予測装置。
【請求項４】
前記ユールニールセン修正ノイゲバウア方程式として式（１０）、式（１１）を用い、前記単色の実効網点面積率を波長の関数として式（１２）を用い、更に前記予測分光反射率を計算する際には式（１３）、式（１４）を用いることを特徴とする請求項３に記載の印刷物分光反射率予測装置。
【数２１】
2014143522.tif 000023
【数２２】
2014143522.tif 000024
【数２３】
2014143522.tif 000025
【数２４】
2014143522.tif 000026
【数２５】
2014143522.tif 000027
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、面積変調により色を再現する印刷物において再現される色の分光反射率を予測する方法及び予測する装置に関するものである。
【背景技術】
【０００２】
従来、印刷物の分光反射率を予測する方法として、設定網点面積率（例えばＣ（シアン）、Ｍ（マゼンタ）、Ｙ（イエロー）、Ｋ（ブラック））と、実際にその掛け合わせを出力したものを測定した分光反射率の関係から、ルックアップテーブル（ＬＵＴ）を作成する方法が用いられている。この方法においては、例えば各一次色を適当な間隔で区切った設定網点面積率の全ての掛け合わせを出力し、その結果を測定することで図４に示すようなＬＵＴを作成する。図４において各格子点は適当な間隔に区切った設定網点面積率に相当し、格子点上には、実際に出力して測定した分光反射率が格納されている。格子点が存在しない掛け合わせ、すなわち実際に測定を行っていない掛け合わせに対しては、周辺の格子点が持つ値から、補間等の手法を用いて、任意の掛け合わせに対する分光反射率を求めることができる。
【０００３】
もう一つの方法として分光ノイゲバウア方程式を利用する方法が知られている。この方法では、各原色のベタおよびベタ同士の全ての掛け合わせと、用紙の分光反射率をあらかじめ測定し、再現される分光反射率を式（１）で予測する。式（１）中のＦｐはノイゲバウア原色の網点面積率であって式（２）で示され、式（２）中のａｍはノイゲバウア原色ｐの掛け合わせに原色Ｉｍが含まれるときであって、一方（１−ａｍ）は含まれないときである。尚、ノイゲバウア原色とは、用紙および、単色のベタとその全ての掛け合わせで再現される色であって、ＣＭＹＫ（シアン、マゼンタ、イエロー、ブラック）の印刷の場合には、１次色Ｃ、  Ｍ、  Ｙ、  Ｋ、２次色ＣＭ、  ＣＹ、  ＣＫ、  ＭＹ、  ＭＫ、
ＹＫ、３次色ＣＭＹ、  ＣＭＫ、  ＣＹＫ、  ＭＹＫ、  ４次色ＣＭＹＫと用紙の１６色がノイゲバウア原色となる。
【０００４】
【数１】
2014143522.tif 000003
【０００５】
【数２】
2014143522.tif 000004
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開平５−１１５００号公報
【特許文献２】特開平７−２０８７８８号公報
【非特許文献】
【０００７】
【非特許文献１】色再現工学の基礎  太田登著  コロナ社  Ｐ．１８２
【発明の概要】
【発明が解決しようとする課題】
【０００８】
上述したＬＵＴを用いる方法によれば、測定するサンプルの数を増加させることにより、補間による誤差が減少し、比較的高い精度での予測が可能となる反面、測定数の増加による作業時間、手間の増大が問題となる。
【０００９】
また、原色の数が増加すると、測定が必要な格子点の数が急激に増大する。例えば各原色を５段階の面積率に区切った場合を想定すると、４原色の印刷では５４＝６２５点の格
子点であるのに対し、６原色の印刷では５６＝１５６２５点となるため、多色を扱う印刷には適用しにくく、仮に多数の格子点に対応する掛け合わせを全て印刷するとしても、その印刷には大きな面積が必要となり、面内ばらつきの大きい印刷機においては、その面内ばらつきの誤差がＬＵＴに含まれてしまうという問題があった。
【００１０】
また、上述した分光ノイゲバウア方程式を用いる方法によれば、測定が必要な掛け合わせは各色ベタとベタ同士の掛け合わせ、と用紙のみとなるため非常に少ない測定で予測が可能となるものの、実際の印刷ではこの方程式が成立することはほとんどなく、予測精度の面で問題があった。ノイゲバウア方程式の現実との乖離を補正する方法として、ユールニールセン修正ノイゲバウア方程式が提案され、広く用いられている。ユールニールセン修正ノイゲバウア方程式は式（３）、式（４）で表される。
【００１１】
【数３】
2014143522.tif 000005
【００１２】
【数４】
2014143522.tif 000006
【００１３】
この式（３）では係数ｎを導入し、実測値と予測値の誤差を最小化するように、ｎを経験的に求める方法をとっている。この手法を用いることにより、通常のノイゲバウア方程式と比較して、より高い精度で分光反射率を予測することができることが知られているが、式（３）の形式から見て明らかなように、一つのパラメータｎを用いて全ての原色の誤差を吸収しているため、原色の特性が異なるような印刷では予測に問題があるとされている。先に挙げた多色印刷の例では、ＡＭスクリーンとＦＭスクリーンのハイブリッドの網点を用いた印刷等が行われることがあり、原色の特性に大きな差があることが想定される。そのため、従来のユールニールセン修正ノイゲバウア方程式では精度上問題が発生することがあった。
【００１４】
そこで本発明は、上述の問題に鑑みてなされたもので、印刷物の分光反射率を精度良く予測する方法及び装置を提供することを目的とする。
【課題を解決するための手段】
【００１５】
そこで本発明に係る請求項１に記載の発明は、面積変調で色を再現する印刷物の任意の掛け合わせに対する分光反射率を予測し、再現色を予測する方法であって、
用紙、各色ベタ、ベタの掛け合わせについて分光反射率（この場合の分光反射率を分光反射率その１と称す）を測定する第１の工程と、
一次色を任意の網点面積率で変化させた場合の分光反射率（この場合の分光反射率を分光反射率その２と称す）を測定する第２の工程と、
第１の工程で測定した分光反射率その１と、第２の工程で取得した分光反射率その２から、最適化によりユールニールセン修正ノイゲバウア方程式の係数ｎを決定する第３の工程と、
前記第３の工程で求めた係数ｎを用いて、第２の工程で測定した分光反射率その２に対して、実効網点面積率を波長の関数として求める第４の工程と、
更に前記第４の工程を全ての分光反射率その２に対して繰り返すことにより、各一次色の設定網点面積率と実効網点面積率の関係を求める第５の工程と、
任意の設定網点面積率の掛け合わせに対する分光反射率を、前記第４の工程で求めた関係から、任意の設定網点面積率をそれぞれ、各一次色の実効網点面積率に変換する第６の工程と、
ユールニールセン修正ノイゲバウア方程式に、前記第３の工程で取得した係数ｎと前記第６の工程で取得した実効網点面積率を代入する第７の工程によって、
予測分光反射率を計算することを特徴とする印刷物分光反射率予測方法である。
【００１６】
本発明に係る請求項２に記載の発明は、前記ユールニールセン修正ノイゲバウア方程式として式（５）、式（６）を用い、前記単色の実効網点面積率を波長の関数として表す式として式（７）を用い、更に前記予測分光反射率を計算する際には式（８）、式（９）を用いることを特徴とする請求項１に記載の印刷物分光反射率予測方法である。
【００１７】
【数５】
2014143522.tif 000007
【００１８】
【数６】
2014143522.tif 000008
【００１９】
【数７】
2014143522.tif 000009
【００２０】
【数８】
2014143522.tif 000010
【００２１】
【数９】
2014143522.tif 000011
【００２２】
本発明に係る請求項３に記載の発明は、面積変調で色を再現する印刷物の任意の掛け合わせに対する分光反射率を予測し、再現色を予測する印刷物分光反射率予測装置であって、
第１分光反射率受付処理部と、第２分光反射率受付処理部と、最適化計算処理部と、面積率変換関数保持部と、設定網点面積率受付処理部と、網点面積率変換処理部と、分光反射率算出処理部と、を備え、
第１分光反射率受付処理部は、用紙、各色ベタ、ベタの掛け合わせについて分光反射率（この場合の分光反射率を分光反射率その１と称す）を取得し記憶する手段で、
第２分光反射率受付処理部は、一次色を任意の面積率で変化させた分光反射率（この場合の分光反射率を分光反射率その２と称す）を取得し記憶する手段で、
最適化計算処理部は、第１分光反射率受付処理部が取得し記憶した分光反射率その１と、第２分光反射率受付処理部が取得し記憶した分光反射率その２と、を読み出し、最適化によりユールニールセン修正ノイゲバウア方程式の係数ｎを決定する手段で、
面積率変換関数保持部は、最適化計算処理部が計算したｎを読み出して、第２分光反射率受付処理部が取得した分光反射率その２に対して、実効網点面積率を波長の関数として求め、この処理を第２分光反射率受付処理部が取得した全ての分光反射率その２に対して繰り返すことにより、各一次色の設定網点面積率と実効網点面積率の関係を算出して保持する手段で、
設定網点面積率受付処理部は、分光反射率を予測する対象として任意の設定網点面積率の掛け合わせを入力として取得する手段で、
網点面積率変換処理部は、設定網点面積率受付処理部が受け取った設定網点面積率を、面積率変換関数保持部が保持する関係を用いて各一次色の実効網点面積率に変換する手段で、
分光反射率算出処理部によって、ユールニールセン修正ノイゲバウア方程式に、最適化計算処理部が算出したｎ値と、網点面積率変換処理部が変換した実効網点面積率を代入し、印刷物の予測分光反射率を計算することを特徴とする印刷物分光反射率予測装置である。
【００２３】
本発明に係る請求項４に記載の発明は、前記ユールニールセン修正ノイゲバウア方程式として式（１０）、式（１１）を用い、前記単色の実効網点面積率を波長の関数として表す式として式（１２）を用い、更に前記予測分光反射率を計算する際には式（１３）、式（１４）を用いることを特徴とする請求項３に記載の印刷物分光反射率予測装置である。
【００２４】
【数１０】
2014143522.tif 000012
【００２５】
【数１１】
2014143522.tif 000013
【００２６】
【数１２】
2014143522.tif 000014
【００２７】
【数１３】
2014143522.tif 000015
【００２８】
【数１４】
2014143522.tif 000016
【発明の効果】
【００２９】
本発明によれば、少数の一次色ステップとベタの掛け合わせおよび用紙の測定値から単色の実効網点面積率を波長の関数として表すことによって、任意の掛け合わせの設定網点面積率に対する分光反射率を精度良く予測することができる。また、必要とする測定が少ないことから、少ない面積で測定のためのパッチを出力することができ、印刷の面内ばらつき等の影響を受けないという特徴もある。さらに、各一次色で発生する誤差を個別に解消するように各一次色の実効網点面積率を波長の関数として求めているため、一次色の特性が異なるような印刷でも高い精度で分光反射率を予測できるという特徴がある。
【図面の簡単な説明】
【００３０】
【図１】本発明の一実施形態による印刷物分光反射率予測方法の処理の流れを表すフローチャートである。
【図２】本発明の一実施形態による印刷物分光反射率予測装置の構成を示すブロック図である。
【図３】本発明の一実施形態による印刷物分光反射率予測装置の処理の流れを表すフローチャートである。
【図４】多次元ＬＵＴの例を示す図である。
【発明を実施するための形態】
【００３１】
以下、本発明の一実施形態による印刷物分光反射率予測方法を図面を参照して説明する。
【００３２】
図１は本発明の一実施形態による印刷物分光反射率予測方法の処理の流れを表すフローチャートである。このフローチャートに沿って各工程（ステップ）での処理の詳細を説明する。なお、この実施形態例ではＣＭＹＫの４原色の面積変調（網点の面積変化）により色を再現する印刷を対象とする。
【００３３】
ステップ１０１（第１の工程）：（用紙、単色ベタ、ベタの掛け合わせの分光反射率その１の測定）
用紙、単色ベタ（Ｃ、Ｍ、Ｙ、Ｋ）、ベタの掛け合わせ（ＣＭ、ＭＹ、ＣＹ、ＣＫ、ＭＫ、ＹＫ、ＣＭＹ、ＣＭＫ、ＭＹＫ、ＣＹＫ、ＣＭＹＫ）の計１６色について分光反射率その１（Ｒｐ（λ））を測定する。単色ベタ及びベタの掛け合わせは既に印刷されているものとする。
【００３４】
ステップ１０２（第２の工程）：（単色の網点面積率を任意の面積率で変化させたパッチの分光反射率その２の測定）
ＣＭＹＫの各単色の網点面積率を任意の面積率（設定網点面積率）で変化させたパッチ（このパッチは既に印刷されているものとする）の分光反射率その２を測定する。本実施例では各単色の面積率変化は１０％刻みとする。
【００３５】
ステップ１０３（第３の工程）（ｎ値の決定）
ステップ１０２で測定した各単色のパッチの実測値と、ユールニールセン修正ノイゲバウア方程式（前記式（５）、前記式（６））に設定網点面積率を代入して得られる予測分光
反射率との誤差が最小となるように、ユールニールセン修正ノイゲバウア方程式の係数ｎを決定する。係数ｎの決定に際しては、全てのパッチの実測値に対して、実測値と予測値のＲＭＳ（Ｒｏｏｔ  Ｍｅａｎ  Ｓｑｕａｒｅ）誤差を計算し、これらの和を評価値として、非線形最適化手法を用いて係数ｎを決定しても良い。また、他の方法としては、ｎを適当な大きさ（例えば０．１）で適当な範囲（例えば０．１〜１０）まで変化させ、上記評価値を随時計算し、評価値が最小となるｎを決定する方法をとっても良い。
【００３６】
ステップ１０４（第４の工程）（実効網点面積率を波長の関数として求める）
前記式（７）を用いて各単色の設定網点面積率に対する実効網点面積率を波長の関数として求める。このときｎはステップ１０３で求めた係数である。
【００３７】
ステップ１０５（第５の工程）（設定網点面積率と実効網点面積率の関係を求める）
ステップ１０２で測定した全ての網点面積率および、ベタと用紙の測定値について前記式（７）の式を計算することで、ステップ１０２で測定した設定網点面積率（この例では１０％刻みのステップ）に対応する実効網点面積率を求めることができる。
【００３８】
ステップ１０１〜１０５は、本分光反射率予測方法において、一度行えば良く、これに基づいて任意の網点面積率の掛け合わせに対する分光反射率を予測することができる。
【００３９】
以下に分光反射率予測の手順をフローチャートを参照しながら説明する。
【００４０】
ステップ１０６（設定網点面積率を実効網点面積率に変換する）
任意の設定網点面積率の掛け合わせ（Ｃ、  Ｍ、  Ｙ、  Ｋ）について、各原色の設定網点面積率を、実効網点面積率に変換する。具体的にはステップ１０５から単色で測定した設定網点面積率に対しては、それに対応する実効網点面積率が計算されている。与えられた任意の掛け合わせの設定網点面積率が、ステップ１０２で測定した設定網点面積率の中に存在する場合には、単にその設定網点面積率に対応する実効網点面積率を取得すれば良い。与えられた任意の掛け合わせの設定網点面積率が、ステップ１０２で測定した設定網点面積率の中に存在しない場合には、その設定網点面積率の前後でステップ１０２で測定した設定網点面積率と実効網点面積率の関係から、補間等の手法により実効網点面積率を求めることができる。（例えば、この例では１０％刻みの設定網点面積率については設定網点面積率と実効網点面積率の関係が求められているため、２５％の設定網点面積率に対する実効網点面積率を求めるためには、２０％と３０％の設定網点面積率に対する実効網点面積率の関係を用いる。）この補間には単純な線形補間を用いても良いし、実効網点面積率の各波長成分について、設定網点面積率との関係をスプライン等の曲線で近似したり、多項式で近似したりして、補間しても良い。
【００４１】
ステップ１０７（分光反射率の予測）
ステップ１０３で決定したｎ値と、ステップ１０６で求めた実効網点面積率をユールニールセン修正ノイゲバウア方程式（前記式（８）、前記式（９））に代入することにより、分光反射率を予測（計算）する。
【００４２】
次に本発明の一実施形態による印刷物分光反射率予測装置を、図面を参照して説明する。図２は、本実施形態による印刷物分光反射率予測装置の構成の一例を示すブロック図である。この図において、符号１は印刷物分光反射率予測装置であり、第１分光反射率受付処理部１１、第２分光反射率受付処理部１２、最適化計算処理部１３、面積率変換関数保持部１４、設定網点面積率受付処理部１５、網点面積率変換処理部１６、分光反射率算出処理部１７を備えている。
【００４３】
第１分光反射率受付処理部１１は、用紙、各色ベタ、ベタの掛け合わせについて分光反
射率その１を取得し記憶する手段で、
第２分光反射率受付処理部１２は、一次色を任意の面積率で変化させた分光反射率その２を取得し記憶する手段で、
最適化計算処理部１３は、第１分光反射率受付処理部１１が取得し記憶した分光反射率その１と、第２分光反射率受付処理部１２が取得し記憶した分光反射率その２と、を読み出し、最適化によりユールニールセン修正ノイゲバウア方程式の係数ｎを決定する手段で、
面積率変換関数保持部１４は、最適化計算処理部１３が計算したｎを読み出して、第２分光反射率受付処理部１２が取得した分光反射率その２に対して、実効網点面積率を波長の関数として求め、この処理を第２分光反射率受付処理部１２が取得した全ての分光反射率その２に対して繰り返すことにより、各一次色の設定網点面積率と実効網点面積率の関係を算出して保持する手段で、
設定網点面積率受付処理部１５は、分光反射率を予測する対象として任意の設定網点面積率の掛け合わせを入力として取得する手段で、
網点面積率変換処理部１６は、設定網点面積率受付処理部１５が受け取った設定網点面積率を、面積率変換関数保持部１４が保持する関係を用いて各一次色の実効網点面積率に変換する手段で、
分光反射率算出処理部１７は、ユールニールセン修正ノイゲバウア方程式に、最適化計算処理部１３が算出したｎ値と、網点面積率変換処理部１６が変換した実効網点面積率を代入し、印刷物の予測分光反射率を計算する。
【００４４】
次に本発明の印刷物分光反射率予測装置をＣＭＹＫの面積変調で色を再現する印刷を対象として説明する。図３は印刷物分光反射率予測装置における分光反射率予測処理の処理フローを示す図である。図３を用いて印刷物分光反射率予測装置における分光反射率予測処理について説明する。
【００４５】
（ステップＳ１）
先ず、ユーザーによって測定された分光反射率その１を印刷物分光反射率予測装置１の第１分光反射率受付処理部１１が取得し一時記憶する。
【００４６】
（ステップＳ２）
次にユーザーは分光反射率その２を測定し、印刷物分光反射率予測装置１に入力する。この実施例では各単色の設定網点面積率１０％刻みのステップを測定することとする。印刷物分光反射率予測装置１では、第２分光反射率受付処理部１２が分光反射率その２を取得し一時記憶する。
【００４７】
（ステップＳ３）
次に最適化計算処理部１３は、第１分光反射率受付処理部１１と、第二分光反射率受付処理部１２が取得し記憶した分光反射率その１と分光反射率その２を読み取る。そしてユールニールセン修正ノイゲバウア方程式（前記式（１０）、前記式（１１））により予測される分光反射率と、各単色ステップの実測の分光反射率の誤差が最も小さくなるように、最適化手法を用いてユールニールセン修正ノイゲバウア方程式の係数ｎを決定する。
【００４８】
（ステップＳ４）
次に面積率変換関数保持部１４は、第一分光反射率受付処理部１１と、第二分光反射率受付処理部１２が取得し記憶した分光反射率その１と分光反射率その２を読み取る。さらに最適化計算処理部１３が決定した係数ｎを読み取る。これらを元に、分光反射率を測定した一次色の設定網点面積率のそれぞれに対して、対応する実効網点面積率を波長の関数として前記式（１２）を用いて計算する。例えばシアン（Ｃ）の設定網点面積率を例にとって説明すると、本実施形態では１０％刻みのステップが測定されているため、例えば、シ
アンの設定網点面積率１０％に対応する実効網点面積率は式（１５）で計算される。
【００４９】
【数１５】
2014143522.tif 000017
同様にしてシアンの設定網点面積率２０％、３０％・・・、９０％に対応する実効網点面積率を計算することができる。他の原色についても同様の式から、実効網点面積率を計算することができる。なお、前記式（１５）から明らかなように、設定網点面積率が０％の時の実効網点面積率は全ての波長において０となり、設定網点面積率が１００％の時の実効網点面積率は全ての波長において１となる。このようにして得られた各原色の０％、１０％、・・・、１００％に対する実効網点面積率をルックアップテーブルとして記憶する。
【００５０】
（ステップＳ５）
次にユーザーは、分光反射率を予測する対象となる設定網点面積率の掛け合わせ（Ｃ、Ｍ、Ｙ、Ｋ）を印刷物分光反射率予測装置１に入力する。印刷物分光反射率予測装置１では設定網点面積率受付処理部１５が入力された設定網点面積率を取得して一時記憶する。
【００５１】
（ステップＳ６）
次に網点面積率変換処理部１６は、面積率変換関数保持部１４が記憶する設定網点面積率と実効網点面積率の関係を参照し、設定網点面積率受付処理部１５が取得して保持する設定網点面積率を、各原色毎に実効網点面積率に変換する。
【００５２】
（ステップＳ７）
最後に、分光反射率算出処理部１７は、修正ユールニールセンノイゲバウア方程式を計算するために必要な情報として、第１分光反射率受付処理部１１からノイゲバウア原色を、最適化計算処理部１３から係数ｎを、網点面積率変換処理部１６から各原色の実効網点面積率を読み取り、前記式（１３）、前記式（１４）に基づいて予測分光反射率を計算する。
【００５３】
このようにして得られた予測分光反射率は、モニタ等の表示装置（図示しない）に結果を表示しても良いし、ファイルとして外部記憶装置（図示しない）に書き出しても良い。なお、同じ印刷に対して複数の異なる掛け合わせに対する分光反射率を予測する際は、上記ステップＳ４までの操作は一度行えば良く、ステップＳ５以降を異なる掛け合わせに対して繰り返すことにより、複数の予測分光反射率を取得することができる。
【００５４】
以上のように、本発明の印刷物分光反射率予測方法及び印刷物分光反射率予測装置によれば、任意の掛け合わせの設定網点面積率に対する分光反射率を精度良く予測することができる。また、分光反射率予測方法、および分光反射率予測装置の一実施例として、ＣＭＹＫの４色の面積変調で色を再現する印刷を対象として説明したが、本発明による印刷物分光反射率予測方法、および印刷物分光反射率予測装置は、任意の原色数を用いる多色印刷にも適用することが出来る。
【符号の説明】
【００５５】
１・・・印刷物分光反射率予測装置
１１・・・第１分光反射率受付処理部
１２・・・第２分光反射率受付処理部
１３・・・最適化計算処理部
１４・・・面積率変換関数保持部
１５・・・設定網点面積率受付処理部
１６・・・網点面積率変換処理部
１７・・・分光反射率算出処理部
【図１】
2014143522.tif 000028
【図２】
2014143522.tif 000029
【図３】
2014143522.tif 000030
【図４】
2014143522.tif 000031