JP4068771B2 - カラープリンタ出力の色再現方法および装置並びにカラープリンタ出力の色再現プログラムを記録したコンピュータ読取可能な記録媒体 - Google Patents
カラープリンタ出力の色再現方法および装置並びにカラープリンタ出力の色再現プログラムを記録したコンピュータ読取可能な記録媒体 Download PDFInfo
- Publication number
- JP4068771B2 JP4068771B2 JP30897099A JP30897099A JP4068771B2 JP 4068771 B2 JP4068771 B2 JP 4068771B2 JP 30897099 A JP30897099 A JP 30897099A JP 30897099 A JP30897099 A JP 30897099A JP 4068771 B2 JP4068771 B2 JP 4068771B2
- Authority
- JP
- Japan
- Prior art keywords
- color
- spectral reflectance
- calculated
- area ratio
- dot area
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N1/00—Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
- H04N1/46—Colour picture communication systems
- H04N1/56—Processing of colour picture signals
- H04N1/60—Colour correction or control
- H04N1/6011—Colour correction or control with simulation on a subsidiary picture reproducer
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Spectrometry And Color Measurement (AREA)
- Facsimile Image Signal Circuits (AREA)
- Color Image Communication Systems (AREA)
- Image Processing (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、例えば電子写真方式のカラープリンタ出力の色を、ＣＲＴ(Cathode-Ray Tube)，ＬＣＤ(Liquid Crystal Display)，ＰＤＰ(Plasma Display Panel)等のモニタで再現する方法および装置、並びに、カラープリンタ出力の色再現プログラムを記録したコンピュータ読取可能な記録媒体に関する。
【０００２】
【従来の技術】
一般に、カラープリンタとモニタとでは色を表現する方法が全く異なるので、カラープリンタの出力（印刷結果）と同じ色をモニタで表現するためには、カラープリンタ用の色信号（ＹＭＣＫ等の色信号）を、モニタ用の色信号（ＲＧＢ等の色信号）に変換する必要がある。このような変換は、カラープリンタのために事前に設定された所定の色票を測色し、ＣＩＥＸＹＺ（１９３１）やＣＩＥＬＡＢなどの標準色空間値を得てから、その標準色空間値とモニタ用の色信号との対応関係を利用することにより実施されるのが通例となっている。
【０００３】
カラープリンタ出力の色をモニタで再現する従来の手法としては、以下のようなものがある。
〔ａ〕ノイゲバウアの式を用いる手法
ノイゲバウアの式は、イエロ（Ｙ），マゼンタ（Ｍ），シアン（Ｃ）およびブラック（Ｋ）の４色のインクの網点面積率（色信号により指定される値）からＲＧＢの各色成分の値を決定する場合に用いられるものである。このようなノイゲバウアの式を用いる手法は、１６種類の色票の測色値を、任意のＹ，Ｍ，Ｃ，Ｋの色成分の値の積を重みとして、平均化するものである。この手法は、印刷の分野で一般的に用いられており、簡単ではあるが理想状態を念頭においているため、特に、カラープリンタについては、その実際の印刷結果とモニタで再現される色とが合わないことが知られている。また、この手法を変形したものも多く存在しているが、上述と同様な問題を含んでいる。
【０００４】
〔ｂ〕特開平１０−７３２８７号公報，特開平１０−７３２８８号公報および特開平１０−７３２８９号公報に開示された手法
この手法は、印刷物（光透過性を有するインクで印刷されたもの）に適用されており、この手法では、３５２色の色票の分光反射率が使用されている。また、４次元色立体（４次元色空間）から３次元色立体（３次元色空間）、３次元色立体から２次元色平面というように、順次、次元を落としていき、再びその逆をたどりながら、求める分光反射率を算出している。なお、色空間または色平面においてブラックグラデーションを表す辺（ブラックのみの辺）についての分光反射率の算出方法は、それ以外の辺についての分光反射率の算出方法とは異なっている。そして、最終的に、印刷物用の色信号がＣＲＴ用の色信号に変換されＣＲＴで表示されている。
【０００５】
【発明が解決しようとする課題】
しかしながら、上述した従来の手法〔ａ〕や〔ｂ〕には、それぞれ、以下のような課題がある。
〔ａ〕ノイゲバウアの式を用いる手法
(a-1) ディザまたは網点がランダムである場合には実際の印刷結果とモニタで再現される色とが一致する傾向にあるが、前述した通り、一般的には一致しないことが判明している。
【０００６】
(a-2) この式を変形したものもいくつか提案されているが、(1-1) の場合と同様、実際の印刷結果とモニタで再現される色とは一致しない。
(a-3) 図２を参照しながら後述するごとく、インキ（色材）が載っている部分の縁辺部において、紙およびベタ色票では再現できない光学現象が生じているが、ノイゲバウアの式では、この光学現象を考慮していない。このため、前述のように、実際の印刷結果とモニタで再現される色とが一致しないものと考えられる。
【０００７】
〔ｂ〕特開平１０−７３２８７号公報，特開平１０−７３２８８号公報および特開平１０−７３２８９号公報に開示された手法
(b-1) これらの従来手法は、いずれも、印刷物（光透過性を有するインキで印刷されたもの）のみを考慮したものである。電子写真方式等を採用したカラープリンタでは、一般に粉体トナーが用いられている。粉体トナーは光を透過させないため、色材として粉体トナーを用いた場合、ある色材の上に異なる色の色材を重ねると、下の色材は見えなくなる。従って、光透過性を有するインキのみを想定した従来手法をカラープリンタに適用しても、当然、モニタでの色再現はうまくいかない。
【０００８】
(b-2) 特に電子写真方式の場合には、ブラックの網点面積率が１００％である際、つまりブラックをベタで印刷すると、他の色材の値にかかわらず、分光反射率がほぼ一定になる。このため、光透過性を有するインキのみを想定した従来手法では、ブラックの網点面積率が１００％である時の分光反射率を推定することができない。
【０００９】
(b-3) 一般に、カラープリンタにおいては、色の変動が大きく、同一のページ内で同一の色信号で印刷された像であっても、その像には印刷箇所によって同じ色にならない場合がある。このような色の不一致をページ内誤差と呼ぶ。このページ内誤差を無視することができないために、高精細な推定を実施することが困難になる場合がある。そのような場合に対応し、分光反射率の推定精度を高めるべく、不一致となる色を測定し、その色について所定の分光反射率が強制的に算出されるようにすること（決め打ち）が望まれている。つまり、特定の色信号（網点面積率）の組み合わせに対して特定の分光反射率を指定できるようにすることが望まれているが、これらの従来手法には、上述の要望を実現するための技術は一切開示されていない。
【００１０】
本発明は、このような課題に鑑み創案されたもので、カラープリンタ出力特有の光学現象を考慮することにより、分光反射率を精度良く推定し、実際のカラープリンタ出力の色とモニタで再現される色とを一致させることができるようにした、カラープリンタ出力の色再現方法および装置並びにカラープリンタ出力の色再現プログラムを記録したコンピュータ読取可能な記録媒体を提供することを目的とする。
【００１１】
【課題を解決するための手段】
上記目的を達成するために、本発明のカラープリンタ出力の色再現方法（請求項１）は、ブラックを含むｎ色（ｎは２以上の自然数）の色材の組み合わせによるカラー印刷を、光透過性を有していない粉体トナーによって行ない且つブラックを最後に刷るカラープリンタの出力を、モニタで再現する方法であって、
ｎ次元色立体を構成するｎ・２n-1 本の辺上において事前に指定された複数の色票（各辺において網点面積率の異なる３以上の基準色票を含む）について測定された分光反射率に基づいて、任意の網点面積率に対応する分光反射率を算出するための係数を算出する係数算出ステップと、
該係数算出ステップで算出された係数もしくは前記複数の色票について測定された前記分光反射率に基づいて、前記任意の網点面積率に対応する分光反射率を算出する分光反射率算出ステップと、
該分光反射率算出ステップで算出された分光反射率を、該モニタ用の色信号に変換する変換ステップとを含み、
該係数算出ステップにおいて、前記ｎ・２n-1 本の辺のうち、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く辺毎に、前記複数の色票について測定された前記分光反射率に基づいて、前記３以上の基準色票の分光反射率にそれぞれ関連した、前記カラープリンタの出力の色特性を表す３以上の係数を算出するとともに、ブラックグラデーションを表す前記辺上の色票について測定された前記分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、その平均値に基づいて、ドットゲインを表す２つの係数ξ，ηを算出し、
該分光反射率算出ステップにおいて、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く前記辺上における、任意の網点面積率に対応する分光反射率を、該係数算出ステップで算出された前記３以上の係数に基づいて算出する一方、ブラックの網点面積率が１００％である前記辺上における、任意の網点面積率での分光反射率を、前記複数の色票について測定された前記分光反射率に基づいて算出し、ブラックグラデーションを表す前記辺上における任意の網点面積率に対応する分光反射率を、該係数算出ステップで算出された前記２つの係数ξ，ηに基づいて算出することを特徴としている。
【００１２】
なお、該モニタで再現すべき色情報が、該カラープリンタ用のものであるか否かを、前記複数の色票について測定された前記分光反射率に基づいて判別する判別ステップを含み、該判別ステップにおいて、ブラックの網点面積率が１００％である各辺上の色票について測定された前記分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、前記ブラックの網点面積率が１００％である各辺上の色票について算出された平均値の標準偏差を算出し、その標準偏差と所定値とを比較し、その標準偏差が前記所定値以下であれば該モニタで再現すべき色情報が該カラープリンタ用のものであると判別され、該判別ステップにおいて該モニタで再現すべき色情報が該カラープリンタ用のものであると判別された場合に、該係数算出ステップ，該分光反射率算出ステップおよび該変換ステップを行なってもよい（請求項２）。
【００１４】
また、該係数算出ステップにおいて、前記複数の色票について測定された前記分光反射率に基づいて、前記３以上の係数を最小二乗法により算出してもよい（請求項３）。
【００１６】
また、該分光反射率算出ステップにおいて、前記複数の色票について測定された前記分光反射率に基づいて算出・決定される一定値を、ブラックの網点面積率が１００％である前記辺上における任意の網点面積率での分光反射率として出力してもよい（請求項４）。
【００１８】
一方、カラープリンタ出力の色再現装置（請求項５）は、ブラックを含むｎ色（ｎは２以上の自然数）の色材の組み合わせによるカラー印刷を、光透過性を有していない粉体トナーによって行ない且つブラックを最後に刷るカラープリンタの出力を、モニタで再現する装置であって、
ｎ次元色立体を構成するｎ・２n-1 本の辺上において事前に指定された複数の色票（各辺において網点面積率の異なる３以上の基準色票を含む）について測定された分光反射率に基づいて、任意の網点面積率に対応する分光反射率を算出するための係数を算出する係数算出部と、
該係数算出部で算出された係数または測定された前記分光反射率に基づいて、前記任意の網点面積率に対応する分光反射率を算出する分光反射率算出部と、
該分光反射率算出部で算出された分光反射率を、該モニタ用の色信号に変換する変換部とをそなえて構成され、
該係数算出部が、前記ｎ・２n-1本の辺のうち、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く辺毎に、前記複数の色票について測定された前記分光反射率に基づいて、前記３以上の基準色票の分光反射率にそれぞれ関連した、前記カラープリンタの出力の色特性を表す３以上の係数を算出するとともに、ブラックグラデーションを表す前記辺上の色票について測定された前記分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、その平均値に基づいて、ドットゲインを表す２つの係数ξ，ηを算出し、
該分光反射率算出部が、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く前記辺上における、任意の網点面積率に対応する分光反射率を、該係数算出部で算出された前記３以上の係数に基づいて算出する一方、ブラックの網点面積率が１００％である前記辺上における、任意の網点面積率での分光反射率を、前記複数の色票について測定された前記分光反射率に基づいて算出し、ブラックグラデーションを表す前記辺上における任意の網点面積率に対応する分光反射率を、該係数算出部で算出された前記２つの係数ξ，ηに基づいて算出することを特徴としている。
【００１９】
このとき、該モニタで再現すべき色情報が、該カラープリンタ用のものであるか否かを、前記複数の色票について測定された前記分光反射率に基づいて判別する判別部であって、ブラックの網点面積率が１００％である各辺上の色票について測定された前記分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、前記ブラックの網点面積率が１００％である各辺上の色票について算出された平均値の標準偏差を算出し、その標準偏差と所定値とを比較し、その標準偏差が前記所定値以下であれば該モニタで再現すべき色情報が該カラープリンタ用のものであると判別する判別部をさらにそなえて構成し、該判別部が、該モニタで再現すべき色情報が該カラープリンタ用のものであると判別された場合に、該係数算出部，該分光反射率算出部および該変換部の動作を実行させるように構成してもよい（請求項６）。
【００２０】
また、本発明のコンピュータ読取可能な記録媒体（請求項７）は、ブラックを含むｎ色（ｎは２以上の自然数）の色材の組み合わせによるカラー印刷を、光透過性を有していない粉体トナーによって行ない且つブラックを最後に刷るカラープリンタの出力を、モニタで再現する機能をコンピュータに実行させるための、カラープリンタ出力の色再現プログラムを記録したものであって、
該色再現プログラムが、
ｎ次元色立体を構成するｎ・２n-1 本の辺上において事前に指定された複数の色票（各辺において網点面積率の異なる３以上の基準色票を含む）について測定された分光反射率に基づいて、任意の網点面積率に対応する分光反射率を算出するための係数を算出する係数算出部、
該係数算出部で算出された係数または前記複数の色票について測定された前記分光反射率に基づいて、前記任意の網点面積率に対応する分光反射率を算出する分光反射率算出部、および、
該分光反射率算出部で算出された分光反射率を、該モニタ用の色信号に変換する変換部として、該コンピュータを機能させ、
該色再現プログラムが、該コンピュータを該係数算出部として機能させる際に、前記ｎ・２n-1本の辺のうち、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く辺毎に、前記複数の色票について測定された前記分光反射率に基づいて、前記３以上の基準色票の分光反射率にそれぞれ関連した、前記カラープリンタの出力の色特性を表す３以上の係数を算出させるとともに、ブラックグラデーションを表す前記辺上の色票について測定された前記分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、その平均値に基づいて、ドットゲインを表す２つの係数ξ，ηを算出させ、
該色再現プログラムが、該コンピュータを該分光反射率算出部として機能させる際に、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く前記辺上における、任意の網点面積率に対応する分光反射率を、該係数算出部で算出された前記３以上の係数に基づいて算出させる一方、ブラックの網点面積率が１００％である前記辺上における、任意の網点面積率での分光反射率を、前記複数の色票について測定された前記分光反射率に基づいて算出させ、ブラックグラデーションを表す前記辺上における任意の網点面積率に対応する分光反射率を、該係数算出部で算出された前記２つの係数ξ，ηに基づいて算出させることを特徴としている。
【００２１】
このとき、該色再現プログラムが、該モニタで再現すべき色情報が、該カラープリンタ用のものであるか否かを、前記複数の色票について測定された前記分光反射率に基づいて判別する判別部であって、ブラックの網点面積率が１００％である各辺上の色票について測定された前記分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、前記ブラックの網点面積率が１００％である各辺上の色票について算出された平均値の標準偏差を算出し、その標準偏差と所定値とを比較し、その標準偏差が前記所定値以下であれば該モニタで再現すべき色情報が該カラープリンタ用のものであると判別する判別部として該コンピュータを機能させ、該モニタで再現すべき色情報が該カラープリンタ用のものであると判別された場合に、該コンピュータを該係数算出部，該分光反射率算出部および該変換部として機能させてもよい（請求項８）。
【００２２】
【発明の実施の形態】
以下、図面を参照して本発明の実施の形態を説明する。
図１は本発明の一実施形態としてのカラープリンタ出力の色再現方法を実現 するためのコンピュータシステム（色再現装置）の構成を示すブロック図であり、この図１に示すように、本実施形態のコンピュータシステムは、ＣＰＵ１０とこのＣＰＵ１０に接続されるバスライン１２とをそなえている。
【００２３】
バスライン１２には、ＲＯＭ１４とＲＡＭ１６とが接続されており、入出力インタフェース（Ｉ／Ｏ）４０を介してキーボード３０，マウス３２，デジタイザ３４，モニタ（例えばＣＲＴ，ＬＣＤ，ＰＤＰ等）３６および磁気ディスク３８が接続されるとともに、入出力インタフェース（Ｉ／Ｏ）４１を介してプリンタ５０が接続されている。
【００２４】
そして、ＲＡＭ１６には、判別部４２，係数算出部４４，分光反射率算出部４６および変換部４８を実現するためのアプリケーションプログラムが格納されており、ＣＰＵ１０が、バスライン１２を介してＲＡＭ１６から上記アプリケーションプログラムを実行することにより、判別部４２，係数算出部４４，分光反射率算出部４６および変換部４８としての機能（その詳細については後述）が実現され、本実施形態のカラープリンタ出力の色再現方法が実行されるようになっている。
【００２５】
これらの判別部４２，係数算出部４４，分光反射率算出部４６および変換部４８としての機能を実現するためのプログラムは、例えばフレキシブルディスク，ＣＤ−ＲＯＭ等の、コンピュータ読取可能な記録媒体に記録された形態で提供される。そして、コンピュータはその記録媒体からプログラムを読み取って内部記憶装置または外部記憶装置に転送し格納して用いる。また、そのプログラムを、例えば磁気ディスク，光ディスク，光磁気ディスク等の記憶装置（記録媒体）に記録しておき、その記憶装置から通信経路を介してコンピュータに提供するようにしてもよい。
【００２６】
判別部４２，係数算出部４４，分光反射率算出部４６および変換部４８としての機能を実現する際には、内部記憶装置（本実施形態ではＲＡＭ１６）に格納されたプログラムがコンピュータのマイクロプロセッサ（本実施形態ではＣＰＵ１０）によって実行される。このとき、記録媒体に記録されたプログラムをコンピュータが直接読み取って実行するようにしてもよい。
【００２７】
なお、本実施形態において、コンピュータとは、ハードウェアとオペレーションシステムとを含む概念であり、オペレーションシステムの制御の下で動作するハードウェアを意味している。また、オペレーションシステムが不要でアプリケーションプログラム単独でハードウェアを動作させるような場合には、そのハードウェア自体がコンピュータに相当する。ハードウェアは、少なくとも、ＣＰＵ等のマイクロプロセッサと、記録媒体に記録されたコンピュータプログラムを読み取るための手段とをそなえている。
【００２８】
上記アプリケーションプログラムは、このようなコンピュータに、判別部４２，係数算出部４４，分光反射率算出部４６および変換部４８としての機能を実現させるプログラムコードを含んでいる。また、その機能の一部は、アプリケーションプログラムではなくオペレーションシステムによって実現されてもよい。
さらに、本実施形態における記録媒体としては、上述したフレキシブルディスク，ＣＤ−ＲＯＭ，磁気ディスク，光ディスク，光磁気ディスクのほか、ＩＣカード，ＲＯＭカートリッジ，磁気テープ，パンチカード，コンピュータの内部記憶装置（ＲＡＭやＲＯＭなどのメモリ），外部記憶装置等や、バーコードなどの符号が印刷された印刷物等の、コンピュータ読取可能な種々の媒体を利用することができる。
【００２９】
さて、図１に示す本実施形態のコンピュータシステムは、ブラック（Ｋ）を含むｎ色（ｎは２以上の自然数；本実施形態ではＹ，Ｍ，Ｃ，Ｋの４色）の色材の組み合わせによるカラー印刷を行ない且つブラック（Ｋ）を最後に刷るカラープリンタの出力を、モニタ３６で再現するものであり、後述するような判別部４２，係数算出部４４，分光反射率算出部４６および変換部４８としての機能を有している。
【００３０】
判別部４２は、モニタ３６で再現すべき色情報が、カラープリンタ用のものであるか否かを、事前に指定された色票の分光反射率に基づいて判別するもので、モニタ３６で再現すべき色情報がカラープリンタ用のものであると判別した場合に、係数算出部４４，分光反射率算出部４６および変換部４８の動作を実行させるようになっている。
【００３１】
この判別部４２は、後で詳述するごとく、ブラックの網点面積率が１００％である辺上の色票についての分光反射率の標準偏差を算出し、その標準偏差と所定値（例えば０．１）とを比較し、その標準偏差が所定値以下であればモニタ３６で再現すべき色情報がカラープリンタ用のものであると判別している。
係数算出部４４は、ｎ次元色立体を構成するｎ・２n-1 本（本実施形態ではｎ＝４とするので、３２本）の辺上において事前に指定された複数の色票（各辺において網点面積率が０％，５０％，１００％である時の３つの基準色票を含む）について測定された分光反射率に基づいて、任意の網点面積率に対応する分光反射率を算出するための係数αi,βi,γi,ξi,ηi （ｉは網点面積率）を算出するものである。
【００３２】
なお、各辺において事前に指定される色票の数をＮとする。本実施形態では、各辺における色票は、網点面積率ｉが０％，１０％，２０％，３０％，４０％，５０％，６０％，７０％，８０％，９０％，１００％である時の１１、つまりＮ＝１１としている。従って、事前に指定される色票の総数は１１×３２＝３５２である。
【００３３】
ここで、係数算出部４４は、４次元色立体を構成する３２本の辺のうち、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く辺毎に、その辺上の色票について測定された分光反射率に基づいて、３つの基準色票の分光反射率Ｓ0%（＝Ｓpaper ) ，Ｓ50% , Ｓ100%にそれぞれ関連した、カラープリンタの出力の色特性を表す３種類の係数αi,βi,γi を最小二乗法により算出する。
【００３４】
また、係数算出部４４は、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く辺毎に、その辺上の色票について測定された分光反射率に基づいて、網点面積率が０％，１００％である時の２つの基準色票Ｓ0%，Ｓ100%の分光反射率にそれぞれ関連した、ドットゲインを表す２種類の係数ξi,ηi も最小二乗法により算出する。
【００３５】
さらに、係数算出部４４は、ブラックグラデーションを表す辺上の色票について測定された分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、その平均値に基づいて、ドットゲインを表す２種類の係数ξi,ηi を算出する。
一方、分光反射率算出部４６は、図３〜図６を参照しながら後述するごとく、係数算出部４４で算出された係数αi,βi,γi,ξi,ηi 、または、事前に指定された色票の分光反射率に基づいて、任意の網点面積率に対応する分光反射率を算出するものである。
【００３６】
ここで、分光反射率算出部４６は、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く辺上における、任意の網点面積率に対応する分光反射率を、係数算出部４４で算出された係数αi,βi,γi に基づいて算出する一方、ブラックの網点面積率が１００％である辺上における、任意の網点面積率での分光反射率を、事前に指定された色票の分光反射率に基づいて算出する。
【００３７】
また、分光反射率算出部４６は、任意の位置に対応する分光反射率を算出する際に、係数算出部４４で算出された係数ξi,ηi を用いるほか、ブラックグラデーションを表す辺上における任意の網点面積率に対応する分光反射率を、係数算出部４４で算出された係数ξi,ηi に基づいて算出する。
さらに、分光反射率算出部４６は、事前に指定された色票の分光反射率に基づいて算出・決定される一定値（例えば０．０３５）を、ブラックの網点面積率が１００％である辺上における任意の網点面積率での分光反射率として出力することもできる。
【００３８】
またさらに、分光反射率算出部４６は、図７および図８を参照しながら後述するごとく、ｎ次元立体から、カラープリンタで最後に刷る第４色材（ブラック）の網点面積率を一定化した３次元立体を算出し、この３次元立体から、カラープリンタで３番目に刷る第３色材（本実施形態ではシアン）の網点面積率を一定化した２次元色平面を算出し、最終的に、この２次元色平面から、係数ξi,ηi を用いて、任意の網点面積率に設定された４色の色材を組み合わせた場合の分光反射率を算出する。
【００３９】
さらにまた、本実施形態の分光反射率算出部４６は、後述するごとく、特定の網点面積率に設定された４色の色材を組み合わせた場合の分光反射率に、特定の値を指定すべく、その特定の網点面積率を与える色信号と前記任意の網点面積率を与える色信号との距離の逆数に基づいて重み係数を算出し、その重み係数を特定の値に乗算した値を、算出された分光反射率に加算するとともに、前記距離が０の場合には、分光反射率として特定の値を出力する機能もそなえている。
【００４０】
そして、変換部４８は、分光反射率算出部４６で算出された分光反射率を、モニタ３６用の色信号（例えばＲＧＢ信号）に変換し、その色信号に基づく色表示をモニタ３６に行なわせるものである。
次に、図２〜図１３を参照しながら、本実施形態の色再現方法について詳細に説明する。
【００４１】
本実施形態においては、任意数ｎの色材の組み合わせから選定されたｎ・２(n-1) 本の辺から得られる分光反射率を使用する。これらの辺を表現するために、以下のような規定▲１▼〜▲３▼を設ける。
▲１▼各辺において、色信号の採る値（即ち、網点面積率）は、(a) ０％、(b) １００％、および、(c) ０％〜１００％の間における複数（本実施形態の１１）の数値の組、のうちのいずれか一つに限り、且つ、(c) の複数の数値の組を採る色信号は１色のみとする。即ち、これらの辺は、モニタ３６上で色を再現されるべきカラープリンタの色域の縁辺を表すものとする。
【００４２】
▲２▼本実施形態では、ｎ＝４とし、例えば、ブラックを含む任意の４つの色材をシアン，マゼンタ，イエロおよびブラックとする。ただし、これ以外の数や色材であってもよいものとする。そして、各辺は、色符号Ｃ，Ｍ，Ｙ，Ｋを用いて表記するものとし、シアン１００％を“Ｃ”、マゼンタ１００％を“Ｍ”、イエロ１００％を“Ｙ”、ブラック１００％を“Ｋ”として表記し、０％ならば、その色符号の表記を行なわない。色符号の表記順序は刷り順通りとするが、０％〜１００％の間における複数の数値の組をとる色信号の色符号については、一番最後に表記する。なお、本実施形態では、刷り順は“ＹＭＣＫ”とするが、これ以外であってもよい。ただし、ブラックだけは最後に刷られるものとする。
【００４３】
▲３▼上記▲１▼，▲２▼の規定に従った場合、例えば下記表１のような数値の組み合わせの辺は、色符号により“ＹＫＣ”として表記されることになる。つまり、イエロおよびブラックの網点面積率が１００％、マゼンタの網点面積率が０％、シアンの網点面積率が、０％〜１００％の間で１０％刻みの値をとる。
【００４４】
【表１】
【００４５】
なお、前述した通り、各辺でのデータの個数をＮで表し、本実施形態では、例えばＮ＝１１とする。
上述のような表記に従い、以下に、本実施形態の色再現方法について詳細に説明する。
〔１〕まず、本実施形態の判別部４２は、カラープリンタ出力の色再現処理を開始するに当たって、一度だけ、下記 (1-1)〜(1-5) のような判定工程を実施する。
【００４６】
(1-1) ブラックの網点面積率が１００％である辺（即ち、ＫＣ，ＫＭ，ＫＹ，ＣＫＭ，ＭＫＣ，ＣＫＹ，ＹＫＣ，ＹＫＭ，ＭＣＫＹ，ＹＣＫＭ，ＹＭＫＣの１２辺）のそれぞれにおいて、事前に指定された色票毎に、可視波長域に含まれる波長についての分光反射率を平均した値を求める。例えば、可視波長域を 380nm〜730nm とした場合、10nm毎の分光反射率（つまり、波長380nm,390nm,400nm,...,730nm の分光反射率；全部で３６の分光反射率）の総和を３６で除算した値を、分光反射率の平均値として求める。このようにして、分光反射率の平均値を、各辺上における１１個の色票のそれぞれについて算出する。なお、可視波長域は 380nm〜730nm に限定されるものではない。
【００４７】
(1-2) 各辺の１１個の色票について算出された分光反射率の平均値の標準偏差を算出する。
(1-3) (1-2) で算出された標準偏差と、予め与えられた所定値（例えば０．１）とを比較する。
(1-4) (1-3) で比較した結果、標準偏差が所定値よりも大きい場合、モニタ３６で再現すべき色情報には、以下に説明する本実施形態の色再現方法は適用されない。
【００４８】
(1-5) (1-3) で比較した結果、標準偏差が所定値以下であれば、今回、モニタ３６で再現すべき色情報がカラープリンタ用のものであると判別され、その色情報には、以下に説明する本実施形態の色再現方法が適用される。つまり、前述したように、係数算出部４４，分光反射率算出部４６および変換部４８が動作する。
【００４９】
ここで、ブラックを最後に刷るカラープリンタにおいて、ブラックの網点面積率が１００％である場合には、例えば図１２に示すように、その他の色シアン，マゼンタ，イエロの網点面積率がどのような値を採ろうとも、波長によらず、分光反射率はほぼ一定になる。
カラープリンタで用いられるトナーは、前述した通り、固体であり、元々、光透過性を有していない。従って、カラープリンタによる印刷結果では、最後に刷られたトナーの色、つまり、一番上に載るトナーの色が見えることになり、その下に何色があるかは色の見え方に関係なくなる。カラープリンタはこのような特徴を有しているので、ブラックの網点面積率が１００％である場合、他の色に関係なく、全ての色票の分光反射率の値は、ほぼ等しくなり、当然、前述した標準偏差も小さくなる。本実施形態の判別部４２では、このような特徴を利用して、モニタ３６で再現すべき色情報がカラープリンタ用のものであるか否かを判別しているのである。
【００５０】
なお、図１２は、ブラック，シアンおよびマゼンタの網点面積率が１００％でイエロの網点面積率が０％〜１００％となるとき（４次色）の分光反射率の測定結果を示すグラフである。この図１２を見て明らかなように、ブラックの網点面積率が１００％である場合、非常に高い精度で、分光反射率を一定値に置き換えることが可能であり、その一定値として、図１２に示す例では“０．０３５”を採用することができる。
【００５１】
〔２〕上述のようにして、カラープリンタ出力の判定を行なった後、本実施形態の係数算出部４４は、一度だけ、下記 (2-1),(2-2)のような係数算出工程を実施する。
(2-1) ブラックの網点面積率が１００％である辺（前記１２本の辺）とブラックグラデーションを表す辺（“Ｋ”として表記される辺）とを除く辺は、下記の通り表記される１９本の辺である。
【００５２】
Ｃ，Ｍ，Ｙ，ＣＭ，ＭＹ，ＣＹ，ＭＣ，ＹＭ，ＹＣ，ＣＫ，ＭＫ，ＹＫ，ＭＣＹ，ＹＣＭ，ＹＭＣ，ＭＣＫ，ＹＣＫ，ＹＭＫ，ＹＭＣＫ
これらの辺の各々に対して、Ｎ個（本実施形態では１１個）の色票の分光反射率が各色信号に対応して存在する。そして、本実施形態の係数算出部４４では、上記１９辺におけるＮ個の分光反射率に対して、下式（１）および（２）を使用して最小２乗法により、各辺のＮ個の色票に対応する係数α，β，γおよび係数ξ，ηが算出される。
【００５３】
ここで、下式（１）および（２）の説明をする前に、係数α，β，γおよび係数ξ，ηのもつ物理的な意味について説明する。
(2-1-1) 係数α，β，γの物理的意味について
図２に示すように、カラープリンタで印刷されるべき色材が紙に載っている場合、目に入る光は、矢印Ａ１，Ａ２，Ａ３で示す３種類の光のいずれかに分類される。
【００５４】
矢印Ａ１は、入射時および反射時の両方とも色材を通過する光（物理的ドットゲイン）を示し、係数γは、このような光の量を表している。また、このような光の影響は、色材の網点面積率が１００％のときに最も顕著になる。
矢印Ａ２は、入射時および反射時の両方とも色材を通過しない光を示し、係数αは、このような光の量を表している。また、このような光の影響は、色材の網点面積率が０％のときに最も顕著になる。
【００５５】
矢印Ａ３は、入射時および反射時のうちの一方のみに色材を通過する光（光学的ドットゲイン）を示し、係数βは、このような光の量を表している。また、このような光の影響は、色材の網点面積率が５０％のときに最も顕著になる。
実際の印刷結果とモニタで再現される色とが、ノイゲバウアの式を用いた場合には一致せず、本発明の方法を用いた場合に一致する理由がここにある。即ち、ノイゲバウアの式では、矢印Ａ１，Ａ２で示すような光は考慮しているが、矢印Ａ３で示すような光の効果（色材が載っている部分の縁辺部において生じる、紙およびベタ色票では再現できない光学現象）を無視している。これに対し、本発明の方法においては、矢印Ａ３で示すような光の効果を考慮し、後述する仮想１次色についての分光反射率を正確に再現できるようにするために、係数α，β，γを使用している。
【００５６】
(2-1-2) 係数ξ，ηの物理的意味について
係数ξ，ηは、前述した係数α，γに対応するものであり、係数ηは、入射時および反射時の両方とも色材を通過する光の量を表し、係数ξは、入射時および反射時の両方とも色材を通過しない光の量を表している。
ここでは、色材の載っていない面積を係数ξで、載っている面積をηで表している。これにより、従来法〔例えばMurray-Davies(マーレイ・デービス）の式〕よりも正確に、色材の載っている面積を把握することが可能になる。なお、これは前記の１９辺の全てに適用可能であり、１次色しか適用できなかった従来法とは大きく異なる点である。即ち、本実施形態では、１次から４次色に至る全ての基本的な色の組み合わせに対して、従来法よりも正確に色材の載っている面積を算出することにより、推定精度を向上させているのである。
【００５７】
【数１】
【００５８】
上式（１），（２）において、
λ：波長(nm)〔入力情報〕
Ｓi （λ）：ｉは、各辺上での網点面積率０％〜１００％におけるＮ個の数値（本実施形態では０，１０，２０，３０，４０，５０，６０，７０，８０，１００の１１個）で、個々の色票の値（網点面積率）を表し、Ｓi （λ）はその色票の分光分布〔出力情報〕
Ｓpaper （λ）〔＝Ｓ0%（λ）〕：紙白の分光反射率〔入力情報〕
Ｓ50% （λ）：網点面積率５０％に対応する色材の載った色票の分光反射率〔入力情報〕
Ｓ100%（λ）：網点面積率１００％に対応する色材の載った色票の分光反射率〔入力情報〕
αi ，βi ，γi ：網点面積率ｉ％の色票に対応する係数α，β，γ〔出力情報〕
ξi ，ηi ：網点面積率ｉ％の色票に対応する係数ξ，η〔出力情報〕
である。なお、〔入力情報〕と記したデータは、上式（１），（２）を用いて最小二乗法を行なう際に、上式（１），（２）に入力されるデータであることを示し、〔出力情報〕と記したデータは、最小二乗法により算出・出力されるデータであることを示している。
【００５９】
(2-2) ブラックグラデーションを表す辺（Ｋ）に属するＮ個（本実施形態では１１個）の色票に対しては、図１１を見れば明らかなように、データの線形性のために最小２乗法を適用することができない。
そこで、本実施形態では、ドットゲインを表す係数ξ，ηを以下のように算出し、分光反射率の再現およびブラックの載っている面積を正確に把握するために、係数ξ，ηを使用する。なお、判別部４２により本発明の手法が適用可能であると判断された場合には、ブラックの透過率が非常に小さいため、図２に矢印Ａ３で示すようなの光は無視することができる。
【００６０】
さて、係数ξ，ηを算出する際には、まず、ブラックグラデーションを表す辺（Ｋ）上のＮ個の色票の分光反射率を、図１１に示すようにグラフ化してから、グラフの平坦部を決定して、その平坦部での分光反射率の平均値を算出する。平均化する範囲（平坦部）は、可視波長域すべてに亘ってもよいし、その可視波長域の一部に亘ってもよい。
【００６１】
そして、Ｎ個の色票の中で、ブラックの網点面積率がｉ％の色票についての平均値をＳi とすると（０≦ｉ≦１００）、下記の連立方程式（３）を解くことにより、その色票に対応する係数ξi ，ηi が算出される。
【００６２】
【数２】
【００６３】
ここで、Ｓpaper は紙白の分光反射率の波長域における平均値〔入力情報〕、Ｓ100%は網点面積率１００％に対応する色材の載った色票の分光反射率の波長域における平均値〔入力情報〕である。
〔３〕上述のようにして係数算出工程を実施し、係数αi ，βi ，γi ，ξi およびηi を算出した後、本実施形態の分光反射率算出部４６は、下記 (3-1)〜(3-3) の工程を実施し、分光反射率の算出・推定を行なう。
【００６４】
以下では、３２辺の各々において、事前に指定された色票以外の色信号（網点面積率）に対応する分光反射率を算出・推定する方法について、項目(3-1) で説明する。なお、以降の説明では、３２辺のそれぞれのことを「仮想１次色」と呼称する。
そして、前記仮想１次色から「仮想２次色」を算出・推定する方法について、項目(3-2) で説明する。「仮想２次色」とは、シアン，マゼンタ，イエロおよびブラックの色信号のうちの２つの網点面積率が任意で、残りの２つの網点面積率が０％または１００％で固定されている色のことを意味する。
【００６５】
その後、項目 (3-1)および(3-2) で説明された方法を用いて、シアン，マゼンタ，イエロおよびブラックの色信号で表される４次元色立体から、順次次元を落としていき、最後に項目(3-2) の方法を実施することにより、任意のシアン，マゼンタ，イエロおよびブラックの色信号に対応する分光反射率を算出．推定できることを項目(3-3) で説明する。
【００６６】
(3-1) 仮想１次色の算出・推定方法
仮想１次色の算出・推定方法に際しては、辺毎に、上述した工程〔２〕で算出された以下のデータを使用する。
▲１▼ブラックの値（網点面積率）が１００％に該当する辺およびブラックグラデーションを表す辺以外の辺
（ａ）Ｎ個の色票におけるαi,βi,γi の値
（ｂ）Ｎ個の色票におけるξi,ηi の値
（ｃ）色信号値（網点面積率）が０，５０，１００％における色票の分光反射率Ｓpaper （λ），Ｓ50% （λ），Ｓ100%（λ）
（ｄ）Ｎ個の色票におけるｉの値（ｉは０％〜１００％の間における複数の数値を表す）
▲２▼ブラックグラデーションを表す辺
（ａ）Ｎ個の色票におけるξi,ηi の値
（ｂ）色信号値が０，１００％における色票の分光反射率Ｓpaper （λ），Ｓ100%（λ）
ｃ）Ｎ個の色票におけるｉの値（ｉは０％〜１００％の間における複数の数値を表す）
▲３▼ブラックの値（網点面積率）が１００％に該当する辺
（ａ）Ｎ個の色票の分光反射率Ｓi （λ）
（ｂ）Ｎ個の色票におけるｉの値（ｉは０％〜１００％の間における複数の数値を表す）
これらのデータを使用して、仮想１次色の分光反射率を算出・推定する方法を以下に記す。
【００６７】
(3-1-1) 以下においては、事前にイエロ，マゼンタ，シアンおよびブラックの色信号（即ちＹ，Ｍ，Ｃ，Ｋの網点面積率）が、％単位で与えられているものとする。ただし、ここでは仮想１次色を扱うので、先に定義した通り、４色のうちのの１色だけを、０％〜１００％の範囲で任意に指定することができ、残りの３色は０％か１００％のいずれか一方であるとする。以後、この任意に指定可能な数値（網点面積率）をｊ（％）で表す。
【００６８】
(3-1-2) この任意の網点面積率ｊ％の色信号に対応する分光反射率Ｓj （λ）を推定するために、式（１）や（２）を先程とは逆に使用する。即ち、
【００６９】
【数３】
【００７０】
ここで、上式（４）は、ブラックの値が１００％に該当する辺およびブラックグラデーションを表す辺以外の辺に適用される。また、上式（５）は、ブラックグラデーションを表す辺に適用される。
各係数αj,βj,γj は、図３に示すように、網点面積率ｊに依存する。分光反射率Ｓpaper （λ）に対する係数αj の値は、網点面積率ｊが０％のときには１であり、５０％および１００％のときには０である。同様に、分光反射率Ｓ50% （λ）に対する係数βj の値は、網点面積率ｊが５０％のときには１であり、０％および１００％のときには０であり、分光反射率Ｓ100%（λ）に対する係数γj の値は、網点面積率ｊが１００％のときには１であり、０％および５０％のときには０である。
【００７１】
また、係数αj は、工程〔２〕で求めたＮ個の係数αi （本実施形態は、ｉ＝0,10,20,30,40,50,60,70,80,90,100の１１個）によって代表されている。同様に、係数βj は、工程〔２〕で求めたＮ個の係数βi によって代表され、係数γj は、工程〔２〕で求めたＮ個の係数γi によって代表されている。
そして、任意の網点面積率ｊにおける係数αj,βj,γj の値は、それぞれ、上述したＮ個の係数αi,βi,γi の間を補間することによって、例えば下式 (6-1)〜(6-3) のごとく決定される。
【００７２】
各係数ξj,ηj も、図４に示すように、網点面積率ｊに依存する。分光反射率Ｓpaper （λ）に対する係数ξj の値は、網点面積率ｊが０％のときには１であり、１００％のときには０である。同様に、分光反射率Ｓ100%（λ）に対する係数ηj の値は、網点面積率ｊが１００％のときには１であり、０％のときには０である。
【００７３】
また、係数ξj は、工程〔２〕で求めたＮ個の係数ξi （本実施形態は、ｉ＝0,10,20,30,40,50,60,70,80,90,100の１１個）によって代表されている。同様に、係数ηj は、工程〔２〕で求めたＮ個の係数ηi によって代表されている。
そして、任意の網点面積率ｊにおける係数ξj,ηj の値は、それぞれ、上述したＮ個の係数ξi,ηi の間を補間することによって、例えば下式(7-1),(7-2) のごとく決定される。
【００７４】
つまり、Ｎ個の色票に対応する、先に述べたｉの値のうちで、ｊに近接しｊをはさむ２つの値ｋとｌとを算出し、下式 (6-1)〜(6-3) および(7-1),(7-2) に代入することにより、係数αj,βj,γj およびξj,ηj を算出する。ここで、ｋ≦ｊ≦ｌとする。そして、下式 (6-1)〜(6-3) および(7-1),(7-2) によってそれぞれ算出された係数αj,βj,γj およびξj,ηj を、それぞれ上式（３）および（４）に代入することで、任意の網点面積率ｊに対応する分光反射率Ｓj （λ）が算出される。
【００７５】
【数４】
【００７６】
上述した式（４),（５), (6-1)〜(6-3) および(7-1),(7-2) を用いて算出・推定された分光反射率を、実測された分光反射率とともに、図９〜図１１に示す。これらの図９〜図１１において、太線が、実測された分光反射率を示し、細線が、本実施形態の方法により算出・推定された分光反射率を示す。
なお、図９は、マゼンタ１次色の分光反射率を、係数α，β，γ〔式（４）＆ (6-1)〜(6-3) 〕を用いて算出・推定した結果を示すグラフであり、図１０は、マゼンタ１次色の分光反射率を、係数ξ，η〔式（５）＆(7-1),(7-2) 〕を用いて算出・推定した結果を示すグラフであり、図１１は、ブラック１次色の分光反射率を、係数ξ，η〔式（５）＆(7-1),(7-2) 〕を用いて推定した結果を示すグラフである。
【００７７】
図９に示すように、係数α，β，γ〔式（４）＆ (6-1)〜(6-3) 〕を用いることにより、極めて高い精度で仮想１次色（マゼンタ１次色）の分光反射率を算出・推定することができる。
また、図１０に示すように、係数ξ，η〔式（５）＆(7-1),(7-2) 〕を用いても、ある程度高い精度で仮想１次色（マゼンタ１次色）の分光反射率を算出・推定することができる。
【００７８】
さらに、図１１に示すように、係数ξ，η〔式（５）＆(7-1),(7-2) 〕を用いることにより、かなり高い精度でブラックグラデーションを表す辺（ブラック１次色）上の任意の位置での分光反射率を算出・推定することができる。
実際に本実施形態の方法を用いて仮想１次色の分光反射率を算出・推定する際には、係数ξ，η〔式（５）＆(7-1),(7-2) 〕ではなく、係数α，β，γ〔式（４）＆ (6-1)〜(6-3) 〕を用いる。ただし、式(7-1),(7-2) により算出される係数ξj ，ηj は、仮想２次色の分光反射率を算出・推定する際に、各辺上の位置を表す数値（ドットゲイン）として使用される。
【００７９】
(3-1-3) ブラックの値（網点面積率）が１００％である辺上における、任意の位置の分光反射率Ｓj （λ）を算出・推定する際には、下式（８）を用いる。即ち、前述した２つの値ｋ，ｌと、既知の分光反射率Ｓk （λ）およびＳl （λ）とを下式（８）に代入し、分光反射率Ｓk （λ）およびＳl （λ）を直接的に比例配分することにより、分光反射率Ｓj （λ）を算出する。
【００８０】
【数５】
【００８１】
(3-2) 仮想２次色の算出・推定方法
仮想２次色の算出・推定方法に際し、本実施形態の分光反射率算出部４６は、項目(3-1) において前述した仮想１次色の算出・推定方法を用い、下記 (3-2-1)〜(3-2-7) の工程を実施する。
(3-2-1) ここでは、色材Ａおよび色材Ｂに関する仮想２次色を算出する方法について説明する。このとき、色材Ａの網点面積率がｘ％で色材Ｂの網点面積率がｙ％である仮想２次色を推定すべき色であるとし、２次元色平面上でのその色の位置を点（ｘ，ｙ）で、またその位置における分光反射率をＳx,y （λ）で表す。なお、点（ｘ，ｙ）はその色の位置の呼称とし、実際の位置を示すものではない。
【００８２】
(3-2-2) 点（ｘ，ｙ）を含む平面を構成する４つの辺を決定する。ここで、各々の辺と、図５に示すような平面とを決定すべく、以下のデータ（ａ）〜（ｄ）を使用する。下記データ（ａ）〜（ｄ）は、すべて、項目(3-1) において前述した方法で算出可能である。なお、図５は、本実施形態における、仮想２次色の算出・推定に必要なデータの位置関係を説明するための図である。
【００８３】
（ａ）（ｉ，０）で表され０≦ｉ≦１００である“色材Ａ辺”上における点（ｘ，０）に対応する分光反射率ＳX,0 （λ）と係数ξx,o ，ηx,0
（ｂ）（ｉ，１００）で表され０≦ｉ≦１００である“色材Ａ，Ｂ１００％辺”上における点（ｘ，１００）に対応する分光反射率ＳX,100 （λ）と係数ξx,100 ，ηx,100 （この“色材Ａ，Ｂ１００％辺”は、図５に示すように、“色材Ａ辺”と対向する）
（ｃ）（０，ｊ）で表され０≦ｊ≦１００である“色材Ｂ辺”上における点（０，ｙ）に対応する分光反射率Ｓo,y （λ）と係数ξo,y ，η0,y
（ｄ）（１００，ｊ）で表され０≦ｊ≦１００である“色材Ａ１００％，Ｂ辺”上における点（１００，ｙ）に対応する分光反射率Ｓ100,y （λ）と係数ξ100,y ，η100,y （この“色材Ａ１００％，Ｂ辺”は、図５に示すように、“色材Ｂ辺”と対向する）
ただし、上述した４辺のなかに、ブラックの値（網点面積率）が１００％である辺が含まれる場合には、その辺と対向する、ブラックの値（網点面積率）が０％である辺の係数ξ，ηをそのまま採用する。
【００８４】
(3-2-3) 前記４辺は、図６に示すように配置される。その配置方法について、以下に説明する。なお、図６は、本実施形態における、仮想２次色の分光反射率の算出・推定方法を説明するための図である。
（ａ）（０，０）を原点（Ｐ）とする。
（ｂ）“色材Ａ辺”および“色材Ｂ辺”を、それぞれ上記原点から水平，垂直に配置する。この際、“色材Ａ辺”の長さ，点（ｘ，０）の位置，“色材Ｂ辺”の長さ，点（０，ｙ）の位置を、以下のように設定する。
【００８５】
“色材Ａ辺”の長さ…ξx,0 ＋ηx,0
点（ｘ，０）の位置…“色材Ａ辺”をηx,0 ：ξx,0 の比で内分する位置
“色材Ｂ辺”の長さ…ξ0,y ＋η0,y
点（０，ｙ）の位置…“色材Ａ辺”をη0,y ：ξ0,y の比で内分する位置
（ｃ）“色材Ａ，Ｂ１００％辺”および“色材Ａ１００％，Ｂ辺”については、以下のように配置する。なお、これらの辺の長さ、および、各辺上における点の算出方法は、前記（ｂ）と同様であるので、その説明は省略する。
【００８６】
図６に示すように、点（１００，０）から長さ（ξ100,y ＋η100,y ）を半径とする円と、点（０，１００）から長さ（ξx,100 ＋ηx,100 ）を半径とする円とを描く。このとき、２つの円は、２点で交わることになる。これらの２つの交点のうち、“色材Ａ辺”と“色材Ｂ辺”とで形成される象限と同一に存在し、且つ、図６に示す平面上の座標値（１００，１００）の位置からの距離が近いものを選択する。そして、“色材Ａ，Ｂ１００％辺”および“色材Ａ１００％，Ｂ辺”は、選択された交点で交わるように配置される。
【００８７】
(3-2-4) 求めるべき分光反射率Ｓx,y （λ）に対応する点（ｘ，ｙ）の位置を、以下のようにして決定する。つまり、点（ｘ，０）と点（ｘ，１００）とを直線で結ぶとともに、点（０，ｙ）と点（１００，ｙ）とを直線で結び、これら２つの直線の交点を点（ｘ，ｙ）の位置とする。
(3-2-5) 各点間の距離を、以下のように定義する。ここで、距離とは図６に示す２次元色平面におけるユークリッド距離のことを意味する。
【００８８】
点（ｘ，ｙ）と点（ｘ，０）との間の距離…ｌ1
点（ｘ，ｙ）と点（０，ｙ）との間の距離…ｌ2
点（ｘ，ｙ）と点（ｘ，１００）との間の距離…ｌ3
点（ｘ，ｙ）と点（１００，ｙ）との間の距離…ｌ4
(3-2-6) 重み係数ｗi （ｉ＝１，２，３，４）を、下式 (9-1)〜(9-4) もしくは(10-1)〜(10-4)のように定義する。なお、前者の式 (9-1)〜(9-4) と後者の式(10-1)〜(10-4)とを混合して使用してはならない。即ち、ｗ1 で前者の式(9-1) を使用する場合には、ｗ2,ｗ3,ｗ4 でも前者の式 (9-2)〜(9-4) を使用する。当然、ｗ1 で後者の式(10-1)を使用する場合には、ｗ2,ｗ3,ｗ4 でも後者の式(10-2)〜(10-4)を使用する。
【００８９】
【数６】
【００９０】
(3-2-7) 分光反射率Ｓx,y （λ）は、上述のごとく算出された分光反射率Ｓx,0 （λ），Ｓ0,y （λ），Ｓx,100 （λ），Ｓ100,y （λ）および重み係数ｗ1,ｗ2,ｗ3,ｗ4 を下式（１１）に代入することにより算出される。
【００９１】
【数７】
【００９２】
(3-3) ４次までの色の分光反射率の算出・推定方法
４次色の算出・推定方法に際し、本実施形態の分光反射率算出部４６は、項目(3-1) において前述した仮想１次色の算出・推定方法と、項目(3-2) において前述した仮想２次色の算出・推定方法とを用い、下記 (3-3-1)〜(3-3-6) の工程を実施する。その工程について、図６および図７を参照しながら説明する。
【００９３】
なお、図７は、本実施形態における、次元数リダクション処理を説明するための図、図８は、本実施形態における、ＹＭＣ３次色の分光反射率の算出・推定方法を説明するための図である。
また、ここでは、事前にイエロ，マゼンタ，シアンおよびブラックの色信号により各色の網点面積率（Ｙ，Ｍ，Ｃ，Ｋ）が％単位で与えられているものとする。ただし、０≦Ｙ，Ｍ，Ｃ，Ｋ≦１００とする。
【００９４】
ところで、光透過性を有する色材を用いた印刷では、ブラックの網点面積率が１００％である時でも３次元色立体を構成されていた。しかし、ブラックの網点面積率が１００％である場合、カラープリンタではブラックが最後に重ねられるため、イエロ，マゼンタ，シアンがどのように印刷されていても、印刷結果は全て同じ色になってしまう。従って、カラープリンタでは、ブラック１００％の３次元色立体は、図７に示すように、ほぼ一点に縮退することになる。このようなカラープリンタの特異性を利用することにより、Ｙ，Ｍ，Ｃ，Ｋからなる４次元色立体の次元数のリダクションを以下のような簡素な手法で行なうことができる。なお、図７において、“Composite black ”とは、３つの色材Ｙ，Ｍ，Ｃを等量混合することによって描かれるブラックのことを言い、１つの色材Ｋで描かれるブラックと区別している。
【００９５】
(3-3-1) まず、図７に示す１２個の四辺形から、項目(3-2) において前述した仮想２次色の算出・推定方法を使用して、最後に刷るブラックの網点面積率が一定値Ｋ％である１２辺を算出することにより、図７に示すように、４次元色立体を、１２辺から成る、ブラックＫ％一定の３次元色立体にリデュース(reduce)する。
【００９６】
（ａ）Ｋ，Ｃ，ＣＫ，ＫＣの辺で構成される四辺形から、ブラックＫ％一定で、シアンの値を可変とする辺１を算出する。
（ｂ）Ｋ，Ｍ，ＭＫ，ＫＭの辺で構成される四辺形から、ブラックＫ％一定で、マゼンタの値を可変とする辺２を算出する。
（ｃ）Ｋ，Ｙ，ＹＫ，ＫＹの辺で構成される四辺形から、ブラックＫ％一定で、イエロの値を可変とする辺３を算出する。
【００９７】
（ｄ）ＣＫ，ＣＭ，ＭＣＫ，ＣＫＭの辺で構成される四辺形から、ブラックＫ％一定で且つシアン１００％一定で、マゼンタの値を可変とする辺４を算出する。
（ｅ）ＭＫ，ＭＣ，ＭＣＫ，ＭＫＣの辺で構成される四辺形から、ブラックＫ％一定で且つマゼンタ１００％一定で、シアンの値を可変とする辺５を算出する。
【００９８】
（ｆ）ＣＫ，ＣＹ，ＹＣＫ，ＣＫＹの辺で構成される四辺形から、ブラックＫ％一定で且つシアン１００％一定で、イエロの値を可変とする辺６を算出する。
（ｇ）ＭＣＫ，ＭＣＹ，ＹＭＣＫ，ＭＣＫＹの辺で構成される四辺形から、ブラックＫ％一定，マゼンタ１００％一定かつシアン１００％一定で、イエロの値を可変とする辺７を算出する。
【００９９】
（ｈ）ＭＫ，ＭＹ，ＹＭＫ，ＭＫＹの辺で構成される四辺形から、ブラックＫ％一定で且つマゼンタ１００％一定で、イエロの値を可変とする辺８を算出する。
（ｉ）ＹＫ，ＹＣ，ＹＣＫ，ＹＫＣの辺で構成される四辺形から、ブラックＫ％一定で且つイエロ１００％一定で、シアンの値を可変とする辺９を算出する。
【０１００】
（ｊ）ＹＣＫ，ＹＣＭ，ＹＭＣＫ，ＹＣＫＭの辺で構成される四辺形から、ブラックＫ％一定，シアン１００％一定およびイエロ１００％一定で、マゼンタの値を可変とする辺１０を算出する。
（ｋ）ＹＭＫ，ＹＭＣ，ＹＭＣＫ，ＹＭＫＣの辺で構成される四辺形から、ブラックＫ％一定，マゼンタ１００％一定およびイエロ１００％一定で、シアンの値を可変とする辺１１を算出する。
【０１０１】
（ｌ）ＹＫ，ＹＭ，ＹＭＫ，ＹＫＭの辺で構成される四辺形から、ブラックＫ％一定で且つイエロ１００％一定で、マゼンタの値を可変とする辺１２を算出する。
(3-3-2) 次に、以下に示す４個の四辺形から、項目(3-2) において前述した仮想２次色の算出・推定方法を使用して、３番目に刷るシアンの網点面積率が一定値Ｃ％である４辺を算出することにより、図７および図８に示すように、３次元色立体を、４辺から成る、ブラックＫ％一定で且つシアンＣ％一定の２次元色平面（仮想２次平面）にリダクションする。つまり、ブラックＫ％一定で且つシアンＣ％一定に対応する辺を各々算出し、任意のマゼンタＭ％および任意のイエロＹ％の組み合わせに対する色空間を表す１つの四辺形（２次元色平面）を得る。
【０１０２】
（ａ）上述のごとく算出された４つの辺１，６，９，３で構成される四辺形から、イエログラデーションを表す辺▲１▼を算出する。
（ｂ）上述のごとく算出された４つの辺９，１０，１１，１２で構成される四辺形から、イエロ１００％一定で且つマゼンタグラデーションを表す辺▲２▼を算出する。
【０１０３】
（ｃ）上述のごとく算出された４つの辺５，８，１１，７で構成される四辺形から、マゼンタ１００％一定で且つイエログラデーションを表す辺▲３▼を算出する。
（ｄ）上述のごとく算出された４つの辺１，４，５，２で構成される四辺形から、マゼンタグラデーションを表す辺▲４▼を算出する。
【０１０４】
なお、図８では、ブラック０％一定のＹＭＣ３次色についての３次元色立体から、シアンＣ％一定でスライスされた仮想２次平面（２次元色平面）が示されている。
(3-3-3) 以上のようにして算出された４つの辺▲１▼〜▲４▼で構成される四辺形から、項目(3-2) において前述した仮想２次色の算出・推定方法を使用して、所望の色信号（Ｋ，Ｃ，Ｍ，Ｙ；Target Color）に対応する分光反射率が、算出・推定される。
【０１０５】
(3-3-4) １次色の分光反射率は、項目(3-1) において前述した仮想１次色の算出・推定方法を使用して、辺Ｋ，Ｃ，Ｍ，Ｙから算出・推定される。
(3-3-5) ２次色の分光反射率は、項目(3-2) において前述した仮想２次色の算出・推定方法を使用して、算出・推定される。
(3-3-6) ３次色の分光反射率は、例えば図８に示すようにブラックＫ％の辺のみから３次元色立体を構成してから、前述した項目(3-3-2) 以降の工程を実施することにより、算出推定される。
【０１０６】
なお、一度、前述した式（４),（５), (6-1)〜(6-3) および(7-1),(7-2) を用いて分光反射率を算出しその分光反射率を式（１１）に代入した後は、式（１）〜（３）を用いて算出された係数α，β，γおよび係数ξ，ηは、精度上、使用できなくなる。それ故、これ以後は、新たに算出された辺については、式（８）を利用して、分光反射率の算出を行なう。
【０１０７】
〔４〕ところで、本実施形態の分光反射率算出部４６は、必要に応じて、決め打ちと呼ばれる工程を下記のように実施する。この工程は、特定の網点面積率（Ｙ，Ｍ，Ｃ，Ｋ）に設定された４色の色材を組み合わせた場合の分光反射率を、特定の値に合わせる必要がある場合（合わせるべきデータが存在する場合）にのみ実施される。それゆえ、必要ないと判断した場合（合わせるべきデータが存在しない場合）には、その工程は実施されない。
【０１０８】
なお、この工程は、項目(3-2) で前述した仮想２次色の算出・推定方法において、重み係数ｗi （ｉ＝１，２，３，４）として式(10-1)〜(10-4)を使用する場合に適用される。以下、その工程について説明する。
(4-1) 合わせるべき色信号の値およびその分光反射率を、以下に定義する。
▲１▼合わせるべき色信号…（Ｋm,Ｃm,Ｍm,Ｙm ），ｍ＝１，２,..,Ｍ
▲２▼合わせるべき分光反射率…Ｓm （λ），ｍ＝１，２,..,Ｍ
(4-2) 本工程は、以下のとおりに実施する。
【０１０９】
(4-2-1) 合わせるべき色信号が１次色の場合には、項目(3-1) で前述した仮想１次色の算出・推定方法において、推定に使用されたＮ個の色票の分光反射率に、さらに、合わせるべき分光反射率Ｓm （λ）を追加し、項目(3-1) で説明した通りの方法を実施すればよい。データの数が増えるだけで、基本的に行なわれる処理は、項目(3-1) で説明した方法と何ら変わらない。
【０１１０】
(4-2-2) 合わせるべき色信号が２次色，３次色，４次色の場合には、以下の処理（ａ）〜（ｄ）を行なう。
（ａ）以下に説明するデータ合わせのための処理は、前記式（１１）により、最終的に分光反射率を算出する際にのみ実施される。
（ｂ）合わせるべき色信号（Ｋm,Ｃm,Ｍm,Ｙm ）と事前に与えられた色信号（Ｋ，Ｃ，Ｍ，Ｙ）との距離を、下式（１２）のように定義する。
【０１１１】
【数８】
【０１１２】
（ｃ）このように定義された距離ｌm+4 に対応する重み係数ｗk を、下式（１３）で与える。
【０１１３】
【数９】
【０１１４】
（ｄ）上述のごとく算出された重み係数ｗk と、合わせるべき分光反射率Ｓm （λ）とから、求める分光反射率Ｓx,y （λ）は、下式（１４）により算出されることになる。なお、式（１２）において、ｌm+4 ＝０となる場合（つまり、合わせるべき色信号Ｋm,Ｃm,Ｍm,Ｙm と事前に与えられた色信号Ｋ，Ｃ，Ｍ，Ｙとが一致している場合）、下式（１５）のように、求める分光反射率Ｓx,y （λ）を合わせるべき分光反射率Ｓm （λ）とする。
【０１１５】
【数１０】
【０１１６】
このようにして、決め打ちと呼ばれる工程を実行することにより、特定の色信号（網点面積率）の組み合わせに対して特定の分光反射率を指定することができる。従って、ページ内誤差により不一致となる色を予め測定しておけば、その色については所定の分光反射率を強制的に算出するようにすることができ、ページ内誤差を無視できるので、高精細な推定を実施することが可能になる。
【０１１７】
〔５〕上述のようにして分光反射率算出部４６により分光反射率を算出した後、本実施形態の変換部４８は、分光反射率算出部４６で算出された分光反射率をモニタ３６用の色信号に変換し、その色信号に基づく色表示をモニタ３６に行なわせる。変換部４８において、算出された分光反射率をＣＩＥＸＹＺ（１９３１）およびｓＲＧＢに変換する方法については、以下の国際規格に準拠する。
【０１１８】
（ａ）分光反射率からＣＩＥＸＹＺ（１９３１）への変換…ＣＩＥ１５．２ − Colorimetry
（ｂ）ＣＩＥＸＹＺ（１９３１）からｓＲＧＢへの変換…ＩＥＣ６１９６６ Ｐａｒｔ２．０（ｓＲＧＢ）;1999
（ｃ）ＣＩＥＸＹＺ（１９３１）からＣＲＴ出力信号への変換…ＩＥＣ６１９６６ Ｐａｒｔ５（ＣＲＴ）;20XX
（ｄ）ＣＩＥＸＹＺ（１９３１）からＬＣＤ出力信号への変換…ＩＥＣ６１９６６ Ｐａｒｔ６（ＬＣＤ）;20XX
（ｅ）ＣＩＥＸＹＺ（１９３１）からＰＤＰ出力信号への変換…ＩＥＣ６１９６６ Ｐａｒｔ７（ＰＤＰ）;20XX
以上、項目〔１〕〜〔５〕により説明した本実施形態のカラープリンタ出力の色再現方法の流れを、図１３に示すフローチャート（ステップＳ１１〜Ｓ２６）に従って説明する。
【０１１９】
まず、事前に与えられた複数（本実施形態では３５２個）の色票について測定された分光反射率に基づいて、判別部４２が、項目〔１〕で説明した工程に従いモニタ３６で再現すべき色情報がカラープリンタ用のものであるか否かを判別する（ステップＳ１１）。その色情報がカラープリンタ用のものでない場合（ステップＳ１２からＮＯルート）、処理を終了する一方、その色情報がカラープリンタ用のものである場合（ステップＳ１２からＹＥＳルート）、係数算出部４４により、項目〔２〕で説明した工程に従って、係数αi ，βi ，γi ，ξi およびηi を算出する（ステップＳ１３）。
【０１２０】
そして、分光反射率算出部４６は、色の次数に応じて、項目〔３〕で説明した工程に従った分光反射率の算出・推定を行なう（ステップＳ１４〜Ｓ２４）。その際、分光反射率の推定を行なうべき色の次数を判断する（ステップＳ１４）。
１次色の場合（ステップＳ１４から“１次色”ルート）、分光反射率算出部４６は、項目(3-1) で説明した方法に従って、分光反射率を算出・推定する（ステップＳ１５）。
【０１２１】
また、２次色の場合（ステップＳ１４から“２次色”ルート）、分光反射率算出部４６は、分光反射率を特定の値に合わせる必要があるか否かを判断し（ステップＳ１６）、合わせる必要がない場合（ステップＳ１６からＮＯルート）、項目(3-2) で説明した方法に従って分光反射率を算出・推定する（ステップＳ１７）。合わせる必要がある場合（ステップＳ１６からＹＥＳルート）、分光反射率算出部４６は、項目〔４〕で説明した方法に従って分光反射率を算出・推定する（ステップＳ１８）。
【０１２２】
そして、３次色または４次色の場合（ステップＳ１４から“３，４次色”ルート）、分光反射率算出部４６は、４次色については、まず項目(3-3-1) で説明した手順で４次元色立体をブラック一定の３次元色立体にリデュースしてから（ステップＳ１９）、さらに、項目(3-3-2) で説明した手順で、その３次元色立体をブラック一定で且つシアン一定の２次元色平面にリデュース（スライス）する（ステップＳ２０）。３次色については、ステップＳ２１のリダクション処理のみ行なう。
【０１２３】
分光反射率算出部４６は、２次元色平面まで次元数のリダクションを行なったところで、ステップＳ１６〜Ｓ１８と同様、分光反射率を特定の値に合わせる必要があるか否かを判断し（ステップＳ２１）、合わせる必要がない場合（ステップＳ２１からＮＯルート）、項目(3-2) で説明した方法に従って分光反射率を算出・推定する（ステップＳ２２）。合わせる必要がある場合（ステップＳ２１からＹＥＳルート）、分光反射率算出部４６は、項目〔４〕で説明した方法に従って分光反射率を算出・推定する（ステップＳ２３）。
【０１２４】
この後、変換部４８は、例えば項目〔５〕の（ａ）に記載の方法により、分光反射率をＣＩＥＸＹＺ（１９３１）に一旦変換してから（ステップＳ２４）、例えば項目〔５〕の（ｂ）に記載の方法により、ＣＩＥＸＹＺ（１９３１）をｓＲＧＢに変換し（ステップＳ２５）、モニタ３６上で所望のカラープリンタ出力の色表示を行なわせている（ステップＳ２６）。
【０１２５】
ここで、分光反射率の実測値と本実施形態の手法による推定値との精度比較結果を、下記の表２および表３に示す。これらの表２，表３では、シアン（Ｃ），マゼンタ（Ｍ），イエロ（Ｙ），ブラック（Ｋ）のうちで、２色の組み合わせすべて（ＣＫ，ＭＣ，ＭＫ，ＹＣ，ＹＫ，ＹＭ）の分光反射率について、本実施形態の手法によって推定した結果〔ΔＥ94の平均値および最大値と、その最大値の位置（％単位）〕が示されている。
【０１２６】
そして、各辺は１１段階の色信号（色票）をもっているため、一つの組み合わせ、例えば、ＣＫの色票は１１×１１＝１２１個存在し、ＣＫ２次色の任意の分光反射率を、推定すべき分光反射率に対応する色位置を取り囲む４辺上における４０色票の分光反射率を用いて推定した。
なお、表２では、重み係数ｗi を式 (9-1)〜(9-4) で算出した時の推定精度結果が示され、表３では、重み係数ｗi を式(10-1)〜(10-4)で算出した時の推定精度結果が示されている。
【０１２７】
これらの表２および表３を見て明らかな通り、式 (9-1)〜(9-4) を用いた場合も式(10-1)〜(10-4)を用いた場合も、ほぼ同程度の推定精度であり、また、ΔＥ94の平均値および最大値からも極めて良好な精度であることが分かる。
【０１２８】
【表２】
【０１２９】
【表３】
【０１３０】
このように、本発明の一実施形態の色再現方法によれば、以下のような作用効果を得ることができる。
（Ｉ）ノイゲバウアの式では考慮されていなかった、図２の矢印Ａ３で示すような光の効果を、分光反射率Ｓ50% （λ）に関連した係数βとして考慮することにより、図９にも示した通り仮想１次色についての分光反射率を正確に算出・推定できるので、実際のカラープリンタ出力の色とモニタ３６で再現される色とを一致させることができる。
【０１３１】
（II）ブラックの網点面積率が１００％である際であっても、分光反射率を正確に算出・推定できるとともに、ブラックグラデーション（ブラック１次色）の分光反射率も図１１に示した通り高い精度で算出・推定できるほか、従来法では推定不可能であったドットゲインを算出することもできる。
（III)項目〔４〕で説明した決め打ちと呼ばれる工程を実行することにより、特定の色信号の組み合わせに対して特定の分光反射率を指定できるので、ページ内誤差により不一致となる色を予め測定しておけば、その色については所定の分光反射率を強制的に算出するようにすることができる。従って、ページ内誤差を無視でき、分光反射率の推定精度を高めることが可能になる。
【０１３２】
（IV）判別部４２の機能を用いてカラープリンタ出力と印刷物とを区別することにより、カラープリンタ出力に限定した高精度な色再現が可能になる。
（Ｖ）従来一般的に使用されている９２８色票(ISO12642, JIS X 9202)や７２８色票に比べ、半分以下の色票（本実施形態では３５２色）で、極めて高い精度で色値（分光反射率）の推定を行なうことができる（図９〜図１１や表２，表３参照）。このように推定に必要な色票の数を低減することができるので、色票の数が少ない程、色面積が小さくなり、ページ内誤差も小さくなる。従って、ページ内誤差に起因する推定精度の低下を解消することができる。また、本発明により、事前に指定される複数の色票が、ｎ次元色立体を構成する全ての辺（ｎ・２n-1 本）上に存在していなければ、分光反射率の推定精度が落ちることを理論的に証明することが可能である。
【０１３３】
（VI）本発明の手法を用いることにより、カラープリンタにおいて、試行錯誤に頼ることなくモニタ３６上で確認することで高精度のＵＣＲ(Under Color Removal: 下色除去）量を決定することが可能になるほか、色域圧縮，拡張および色の見えなどの確認を、モニタ３６上で、より容易に行なうことが可能になる。
なお、本発明は上述した実施形態に限定されるものではなく、本発明の趣旨を逸脱しない範囲で種々変形して実施することができる。
【０１３４】
例えば、上述した実施形態では、４色の色材（Ｙ，Ｍ，Ｃ，Ｋ）の組み合わせによるカラープリンタを対象として説明したが、本発明は、この色数に限定されるものではない。
【０１３５】
【発明の効果】
以上詳述したように、本発明のカラープリンタ出力の色再現方法（請求項１〜９）および装置（請求項１０，１１）、並びに、カラープリンタ出力の色再現プログラムを記録したコンピュータ読取可能な記録媒体（請求項１２，１３）によれば、カラープリンタ出力特有の光学現象を考慮することにより、分光反射率を精度良く推定できるので、実際のカラープリンタ出力の色とモニタで再現される色とを一致させることができる。
【図面の簡単な説明】
【図１】本発明の一実施形態としてのカラープリンタ出力の色再現方法を実現するためのコンピュータシステム（色再現装置）の構成を示すブロック図である。
【図２】本実施形態における、分光反射率を算出するための係数の物理的な意味を説明するための図である。
【図３】本実施形態における、係数α，β，γを用いた分光反射率算出手法を説明するための図である。
【図４】本実施形態における、係数ξ，ηを用いた分光反射率算出手法を説明するための図である。
【図５】本実施形態における、仮想２次色の算出・推定に必要なデータの位置関係を説明するための図である。
【図６】本実施形態における、仮想２次色の分光反射率の算出・推定方法を説明するための図である。
【図７】本実施形態における、次元数リダクション処理を説明するための図である。
【図８】本実施形態における、ＹＭＣ３次色の分光反射率の算出・推定方法を説明するための図である。
【図９】マゼンタ１次色の分光反射率を、係数α，β，γを用いて算出・推定した結果を示すグラフである。
【図１０】マゼンタ１次色の分光反射率を、係数ξ，ηを用いて算出・推定した結果を示すグラフである。
【図１１】ブラック１次色の分光反射率を、係数ξ，ηを用いて算出・推定した結果を示すグラフである。
【図１２】ブラック，シアンおよびマゼンタの網点面積率が１００％でイエロの網点面積率が０％〜１００％となるとき（４次色）の分光反射率の測定結果を示すグラフである。
【図１３】本実施形態のカラープリンタ出力の色再現方法の流れを説明するためのフローチャートである。
【符号の説明】
１０ ＣＰＵ
１２ バスライン
１４ ＲＯＭ
１６ ＲＡＭ（記録媒体）
３０ キーボード
３２ マウス
３４ デジタイザ
３６ モニタ
３８ 磁気ディスク
４０，４１ 入出力インタフェース
４２ 判別部
４４ 係数算出部
４６ 分光反射率算出部
４８ 変換部
５０ プリンタ
【発明の属する技術分野】
本発明は、例えば電子写真方式のカラープリンタ出力の色を、ＣＲＴ(Cathode-Ray Tube)，ＬＣＤ(Liquid Crystal Display)，ＰＤＰ(Plasma Display Panel)等のモニタで再現する方法および装置、並びに、カラープリンタ出力の色再現プログラムを記録したコンピュータ読取可能な記録媒体に関する。
【０００２】
【従来の技術】
一般に、カラープリンタとモニタとでは色を表現する方法が全く異なるので、カラープリンタの出力（印刷結果）と同じ色をモニタで表現するためには、カラープリンタ用の色信号（ＹＭＣＫ等の色信号）を、モニタ用の色信号（ＲＧＢ等の色信号）に変換する必要がある。このような変換は、カラープリンタのために事前に設定された所定の色票を測色し、ＣＩＥＸＹＺ（１９３１）やＣＩＥＬＡＢなどの標準色空間値を得てから、その標準色空間値とモニタ用の色信号との対応関係を利用することにより実施されるのが通例となっている。
【０００３】
カラープリンタ出力の色をモニタで再現する従来の手法としては、以下のようなものがある。
〔ａ〕ノイゲバウアの式を用いる手法
ノイゲバウアの式は、イエロ（Ｙ），マゼンタ（Ｍ），シアン（Ｃ）およびブラック（Ｋ）の４色のインクの網点面積率（色信号により指定される値）からＲＧＢの各色成分の値を決定する場合に用いられるものである。このようなノイゲバウアの式を用いる手法は、１６種類の色票の測色値を、任意のＹ，Ｍ，Ｃ，Ｋの色成分の値の積を重みとして、平均化するものである。この手法は、印刷の分野で一般的に用いられており、簡単ではあるが理想状態を念頭においているため、特に、カラープリンタについては、その実際の印刷結果とモニタで再現される色とが合わないことが知られている。また、この手法を変形したものも多く存在しているが、上述と同様な問題を含んでいる。
【０００４】
〔ｂ〕特開平１０−７３２８７号公報，特開平１０−７３２８８号公報および特開平１０−７３２８９号公報に開示された手法
この手法は、印刷物（光透過性を有するインクで印刷されたもの）に適用されており、この手法では、３５２色の色票の分光反射率が使用されている。また、４次元色立体（４次元色空間）から３次元色立体（３次元色空間）、３次元色立体から２次元色平面というように、順次、次元を落としていき、再びその逆をたどりながら、求める分光反射率を算出している。なお、色空間または色平面においてブラックグラデーションを表す辺（ブラックのみの辺）についての分光反射率の算出方法は、それ以外の辺についての分光反射率の算出方法とは異なっている。そして、最終的に、印刷物用の色信号がＣＲＴ用の色信号に変換されＣＲＴで表示されている。
【０００５】
【発明が解決しようとする課題】
しかしながら、上述した従来の手法〔ａ〕や〔ｂ〕には、それぞれ、以下のような課題がある。
〔ａ〕ノイゲバウアの式を用いる手法
(a-1) ディザまたは網点がランダムである場合には実際の印刷結果とモニタで再現される色とが一致する傾向にあるが、前述した通り、一般的には一致しないことが判明している。
【０００６】
(a-2) この式を変形したものもいくつか提案されているが、(1-1) の場合と同様、実際の印刷結果とモニタで再現される色とは一致しない。
(a-3) 図２を参照しながら後述するごとく、インキ（色材）が載っている部分の縁辺部において、紙およびベタ色票では再現できない光学現象が生じているが、ノイゲバウアの式では、この光学現象を考慮していない。このため、前述のように、実際の印刷結果とモニタで再現される色とが一致しないものと考えられる。
【０００７】
〔ｂ〕特開平１０−７３２８７号公報，特開平１０−７３２８８号公報および特開平１０−７３２８９号公報に開示された手法
(b-1) これらの従来手法は、いずれも、印刷物（光透過性を有するインキで印刷されたもの）のみを考慮したものである。電子写真方式等を採用したカラープリンタでは、一般に粉体トナーが用いられている。粉体トナーは光を透過させないため、色材として粉体トナーを用いた場合、ある色材の上に異なる色の色材を重ねると、下の色材は見えなくなる。従って、光透過性を有するインキのみを想定した従来手法をカラープリンタに適用しても、当然、モニタでの色再現はうまくいかない。
【０００８】
(b-2) 特に電子写真方式の場合には、ブラックの網点面積率が１００％である際、つまりブラックをベタで印刷すると、他の色材の値にかかわらず、分光反射率がほぼ一定になる。このため、光透過性を有するインキのみを想定した従来手法では、ブラックの網点面積率が１００％である時の分光反射率を推定することができない。
【０００９】
(b-3) 一般に、カラープリンタにおいては、色の変動が大きく、同一のページ内で同一の色信号で印刷された像であっても、その像には印刷箇所によって同じ色にならない場合がある。このような色の不一致をページ内誤差と呼ぶ。このページ内誤差を無視することができないために、高精細な推定を実施することが困難になる場合がある。そのような場合に対応し、分光反射率の推定精度を高めるべく、不一致となる色を測定し、その色について所定の分光反射率が強制的に算出されるようにすること（決め打ち）が望まれている。つまり、特定の色信号（網点面積率）の組み合わせに対して特定の分光反射率を指定できるようにすることが望まれているが、これらの従来手法には、上述の要望を実現するための技術は一切開示されていない。
【００１０】
本発明は、このような課題に鑑み創案されたもので、カラープリンタ出力特有の光学現象を考慮することにより、分光反射率を精度良く推定し、実際のカラープリンタ出力の色とモニタで再現される色とを一致させることができるようにした、カラープリンタ出力の色再現方法および装置並びにカラープリンタ出力の色再現プログラムを記録したコンピュータ読取可能な記録媒体を提供することを目的とする。
【００１１】
【課題を解決するための手段】
上記目的を達成するために、本発明のカラープリンタ出力の色再現方法（請求項１）は、ブラックを含むｎ色（ｎは２以上の自然数）の色材の組み合わせによるカラー印刷を、光透過性を有していない粉体トナーによって行ない且つブラックを最後に刷るカラープリンタの出力を、モニタで再現する方法であって、
ｎ次元色立体を構成するｎ・２n-1 本の辺上において事前に指定された複数の色票（各辺において網点面積率の異なる３以上の基準色票を含む）について測定された分光反射率に基づいて、任意の網点面積率に対応する分光反射率を算出するための係数を算出する係数算出ステップと、
該係数算出ステップで算出された係数もしくは前記複数の色票について測定された前記分光反射率に基づいて、前記任意の網点面積率に対応する分光反射率を算出する分光反射率算出ステップと、
該分光反射率算出ステップで算出された分光反射率を、該モニタ用の色信号に変換する変換ステップとを含み、
該係数算出ステップにおいて、前記ｎ・２n-1 本の辺のうち、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く辺毎に、前記複数の色票について測定された前記分光反射率に基づいて、前記３以上の基準色票の分光反射率にそれぞれ関連した、前記カラープリンタの出力の色特性を表す３以上の係数を算出するとともに、ブラックグラデーションを表す前記辺上の色票について測定された前記分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、その平均値に基づいて、ドットゲインを表す２つの係数ξ，ηを算出し、
該分光反射率算出ステップにおいて、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く前記辺上における、任意の網点面積率に対応する分光反射率を、該係数算出ステップで算出された前記３以上の係数に基づいて算出する一方、ブラックの網点面積率が１００％である前記辺上における、任意の網点面積率での分光反射率を、前記複数の色票について測定された前記分光反射率に基づいて算出し、ブラックグラデーションを表す前記辺上における任意の網点面積率に対応する分光反射率を、該係数算出ステップで算出された前記２つの係数ξ，ηに基づいて算出することを特徴としている。
【００１２】
なお、該モニタで再現すべき色情報が、該カラープリンタ用のものであるか否かを、前記複数の色票について測定された前記分光反射率に基づいて判別する判別ステップを含み、該判別ステップにおいて、ブラックの網点面積率が１００％である各辺上の色票について測定された前記分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、前記ブラックの網点面積率が１００％である各辺上の色票について算出された平均値の標準偏差を算出し、その標準偏差と所定値とを比較し、その標準偏差が前記所定値以下であれば該モニタで再現すべき色情報が該カラープリンタ用のものであると判別され、該判別ステップにおいて該モニタで再現すべき色情報が該カラープリンタ用のものであると判別された場合に、該係数算出ステップ，該分光反射率算出ステップおよび該変換ステップを行なってもよい（請求項２）。
【００１４】
また、該係数算出ステップにおいて、前記複数の色票について測定された前記分光反射率に基づいて、前記３以上の係数を最小二乗法により算出してもよい（請求項３）。
【００１６】
また、該分光反射率算出ステップにおいて、前記複数の色票について測定された前記分光反射率に基づいて算出・決定される一定値を、ブラックの網点面積率が１００％である前記辺上における任意の網点面積率での分光反射率として出力してもよい（請求項４）。
【００１８】
一方、カラープリンタ出力の色再現装置（請求項５）は、ブラックを含むｎ色（ｎは２以上の自然数）の色材の組み合わせによるカラー印刷を、光透過性を有していない粉体トナーによって行ない且つブラックを最後に刷るカラープリンタの出力を、モニタで再現する装置であって、
ｎ次元色立体を構成するｎ・２n-1 本の辺上において事前に指定された複数の色票（各辺において網点面積率の異なる３以上の基準色票を含む）について測定された分光反射率に基づいて、任意の網点面積率に対応する分光反射率を算出するための係数を算出する係数算出部と、
該係数算出部で算出された係数または測定された前記分光反射率に基づいて、前記任意の網点面積率に対応する分光反射率を算出する分光反射率算出部と、
該分光反射率算出部で算出された分光反射率を、該モニタ用の色信号に変換する変換部とをそなえて構成され、
該係数算出部が、前記ｎ・２n-1本の辺のうち、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く辺毎に、前記複数の色票について測定された前記分光反射率に基づいて、前記３以上の基準色票の分光反射率にそれぞれ関連した、前記カラープリンタの出力の色特性を表す３以上の係数を算出するとともに、ブラックグラデーションを表す前記辺上の色票について測定された前記分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、その平均値に基づいて、ドットゲインを表す２つの係数ξ，ηを算出し、
該分光反射率算出部が、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く前記辺上における、任意の網点面積率に対応する分光反射率を、該係数算出部で算出された前記３以上の係数に基づいて算出する一方、ブラックの網点面積率が１００％である前記辺上における、任意の網点面積率での分光反射率を、前記複数の色票について測定された前記分光反射率に基づいて算出し、ブラックグラデーションを表す前記辺上における任意の網点面積率に対応する分光反射率を、該係数算出部で算出された前記２つの係数ξ，ηに基づいて算出することを特徴としている。
【００１９】
このとき、該モニタで再現すべき色情報が、該カラープリンタ用のものであるか否かを、前記複数の色票について測定された前記分光反射率に基づいて判別する判別部であって、ブラックの網点面積率が１００％である各辺上の色票について測定された前記分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、前記ブラックの網点面積率が１００％である各辺上の色票について算出された平均値の標準偏差を算出し、その標準偏差と所定値とを比較し、その標準偏差が前記所定値以下であれば該モニタで再現すべき色情報が該カラープリンタ用のものであると判別する判別部をさらにそなえて構成し、該判別部が、該モニタで再現すべき色情報が該カラープリンタ用のものであると判別された場合に、該係数算出部，該分光反射率算出部および該変換部の動作を実行させるように構成してもよい（請求項６）。
【００２０】
また、本発明のコンピュータ読取可能な記録媒体（請求項７）は、ブラックを含むｎ色（ｎは２以上の自然数）の色材の組み合わせによるカラー印刷を、光透過性を有していない粉体トナーによって行ない且つブラックを最後に刷るカラープリンタの出力を、モニタで再現する機能をコンピュータに実行させるための、カラープリンタ出力の色再現プログラムを記録したものであって、
該色再現プログラムが、
ｎ次元色立体を構成するｎ・２n-1 本の辺上において事前に指定された複数の色票（各辺において網点面積率の異なる３以上の基準色票を含む）について測定された分光反射率に基づいて、任意の網点面積率に対応する分光反射率を算出するための係数を算出する係数算出部、
該係数算出部で算出された係数または前記複数の色票について測定された前記分光反射率に基づいて、前記任意の網点面積率に対応する分光反射率を算出する分光反射率算出部、および、
該分光反射率算出部で算出された分光反射率を、該モニタ用の色信号に変換する変換部として、該コンピュータを機能させ、
該色再現プログラムが、該コンピュータを該係数算出部として機能させる際に、前記ｎ・２n-1本の辺のうち、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く辺毎に、前記複数の色票について測定された前記分光反射率に基づいて、前記３以上の基準色票の分光反射率にそれぞれ関連した、前記カラープリンタの出力の色特性を表す３以上の係数を算出させるとともに、ブラックグラデーションを表す前記辺上の色票について測定された前記分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、その平均値に基づいて、ドットゲインを表す２つの係数ξ，ηを算出させ、
該色再現プログラムが、該コンピュータを該分光反射率算出部として機能させる際に、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く前記辺上における、任意の網点面積率に対応する分光反射率を、該係数算出部で算出された前記３以上の係数に基づいて算出させる一方、ブラックの網点面積率が１００％である前記辺上における、任意の網点面積率での分光反射率を、前記複数の色票について測定された前記分光反射率に基づいて算出させ、ブラックグラデーションを表す前記辺上における任意の網点面積率に対応する分光反射率を、該係数算出部で算出された前記２つの係数ξ，ηに基づいて算出させることを特徴としている。
【００２１】
このとき、該色再現プログラムが、該モニタで再現すべき色情報が、該カラープリンタ用のものであるか否かを、前記複数の色票について測定された前記分光反射率に基づいて判別する判別部であって、ブラックの網点面積率が１００％である各辺上の色票について測定された前記分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、前記ブラックの網点面積率が１００％である各辺上の色票について算出された平均値の標準偏差を算出し、その標準偏差と所定値とを比較し、その標準偏差が前記所定値以下であれば該モニタで再現すべき色情報が該カラープリンタ用のものであると判別する判別部として該コンピュータを機能させ、該モニタで再現すべき色情報が該カラープリンタ用のものであると判別された場合に、該コンピュータを該係数算出部，該分光反射率算出部および該変換部として機能させてもよい（請求項８）。
【００２２】
【発明の実施の形態】
以下、図面を参照して本発明の実施の形態を説明する。
図１は本発明の一実施形態としてのカラープリンタ出力の色再現方法を実現 するためのコンピュータシステム（色再現装置）の構成を示すブロック図であり、この図１に示すように、本実施形態のコンピュータシステムは、ＣＰＵ１０とこのＣＰＵ１０に接続されるバスライン１２とをそなえている。
【００２３】
バスライン１２には、ＲＯＭ１４とＲＡＭ１６とが接続されており、入出力インタフェース（Ｉ／Ｏ）４０を介してキーボード３０，マウス３２，デジタイザ３４，モニタ（例えばＣＲＴ，ＬＣＤ，ＰＤＰ等）３６および磁気ディスク３８が接続されるとともに、入出力インタフェース（Ｉ／Ｏ）４１を介してプリンタ５０が接続されている。
【００２４】
そして、ＲＡＭ１６には、判別部４２，係数算出部４４，分光反射率算出部４６および変換部４８を実現するためのアプリケーションプログラムが格納されており、ＣＰＵ１０が、バスライン１２を介してＲＡＭ１６から上記アプリケーションプログラムを実行することにより、判別部４２，係数算出部４４，分光反射率算出部４６および変換部４８としての機能（その詳細については後述）が実現され、本実施形態のカラープリンタ出力の色再現方法が実行されるようになっている。
【００２５】
これらの判別部４２，係数算出部４４，分光反射率算出部４６および変換部４８としての機能を実現するためのプログラムは、例えばフレキシブルディスク，ＣＤ−ＲＯＭ等の、コンピュータ読取可能な記録媒体に記録された形態で提供される。そして、コンピュータはその記録媒体からプログラムを読み取って内部記憶装置または外部記憶装置に転送し格納して用いる。また、そのプログラムを、例えば磁気ディスク，光ディスク，光磁気ディスク等の記憶装置（記録媒体）に記録しておき、その記憶装置から通信経路を介してコンピュータに提供するようにしてもよい。
【００２６】
判別部４２，係数算出部４４，分光反射率算出部４６および変換部４８としての機能を実現する際には、内部記憶装置（本実施形態ではＲＡＭ１６）に格納されたプログラムがコンピュータのマイクロプロセッサ（本実施形態ではＣＰＵ１０）によって実行される。このとき、記録媒体に記録されたプログラムをコンピュータが直接読み取って実行するようにしてもよい。
【００２７】
なお、本実施形態において、コンピュータとは、ハードウェアとオペレーションシステムとを含む概念であり、オペレーションシステムの制御の下で動作するハードウェアを意味している。また、オペレーションシステムが不要でアプリケーションプログラム単独でハードウェアを動作させるような場合には、そのハードウェア自体がコンピュータに相当する。ハードウェアは、少なくとも、ＣＰＵ等のマイクロプロセッサと、記録媒体に記録されたコンピュータプログラムを読み取るための手段とをそなえている。
【００２８】
上記アプリケーションプログラムは、このようなコンピュータに、判別部４２，係数算出部４４，分光反射率算出部４６および変換部４８としての機能を実現させるプログラムコードを含んでいる。また、その機能の一部は、アプリケーションプログラムではなくオペレーションシステムによって実現されてもよい。
さらに、本実施形態における記録媒体としては、上述したフレキシブルディスク，ＣＤ−ＲＯＭ，磁気ディスク，光ディスク，光磁気ディスクのほか、ＩＣカード，ＲＯＭカートリッジ，磁気テープ，パンチカード，コンピュータの内部記憶装置（ＲＡＭやＲＯＭなどのメモリ），外部記憶装置等や、バーコードなどの符号が印刷された印刷物等の、コンピュータ読取可能な種々の媒体を利用することができる。
【００２９】
さて、図１に示す本実施形態のコンピュータシステムは、ブラック（Ｋ）を含むｎ色（ｎは２以上の自然数；本実施形態ではＹ，Ｍ，Ｃ，Ｋの４色）の色材の組み合わせによるカラー印刷を行ない且つブラック（Ｋ）を最後に刷るカラープリンタの出力を、モニタ３６で再現するものであり、後述するような判別部４２，係数算出部４４，分光反射率算出部４６および変換部４８としての機能を有している。
【００３０】
判別部４２は、モニタ３６で再現すべき色情報が、カラープリンタ用のものであるか否かを、事前に指定された色票の分光反射率に基づいて判別するもので、モニタ３６で再現すべき色情報がカラープリンタ用のものであると判別した場合に、係数算出部４４，分光反射率算出部４６および変換部４８の動作を実行させるようになっている。
【００３１】
この判別部４２は、後で詳述するごとく、ブラックの網点面積率が１００％である辺上の色票についての分光反射率の標準偏差を算出し、その標準偏差と所定値（例えば０．１）とを比較し、その標準偏差が所定値以下であればモニタ３６で再現すべき色情報がカラープリンタ用のものであると判別している。
係数算出部４４は、ｎ次元色立体を構成するｎ・２n-1 本（本実施形態ではｎ＝４とするので、３２本）の辺上において事前に指定された複数の色票（各辺において網点面積率が０％，５０％，１００％である時の３つの基準色票を含む）について測定された分光反射率に基づいて、任意の網点面積率に対応する分光反射率を算出するための係数αi,βi,γi,ξi,ηi （ｉは網点面積率）を算出するものである。
【００３２】
なお、各辺において事前に指定される色票の数をＮとする。本実施形態では、各辺における色票は、網点面積率ｉが０％，１０％，２０％，３０％，４０％，５０％，６０％，７０％，８０％，９０％，１００％である時の１１、つまりＮ＝１１としている。従って、事前に指定される色票の総数は１１×３２＝３５２である。
【００３３】
ここで、係数算出部４４は、４次元色立体を構成する３２本の辺のうち、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く辺毎に、その辺上の色票について測定された分光反射率に基づいて、３つの基準色票の分光反射率Ｓ0%（＝Ｓpaper ) ，Ｓ50% , Ｓ100%にそれぞれ関連した、カラープリンタの出力の色特性を表す３種類の係数αi,βi,γi を最小二乗法により算出する。
【００３４】
また、係数算出部４４は、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く辺毎に、その辺上の色票について測定された分光反射率に基づいて、網点面積率が０％，１００％である時の２つの基準色票Ｓ0%，Ｓ100%の分光反射率にそれぞれ関連した、ドットゲインを表す２種類の係数ξi,ηi も最小二乗法により算出する。
【００３５】
さらに、係数算出部４４は、ブラックグラデーションを表す辺上の色票について測定された分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、その平均値に基づいて、ドットゲインを表す２種類の係数ξi,ηi を算出する。
一方、分光反射率算出部４６は、図３〜図６を参照しながら後述するごとく、係数算出部４４で算出された係数αi,βi,γi,ξi,ηi 、または、事前に指定された色票の分光反射率に基づいて、任意の網点面積率に対応する分光反射率を算出するものである。
【００３６】
ここで、分光反射率算出部４６は、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く辺上における、任意の網点面積率に対応する分光反射率を、係数算出部４４で算出された係数αi,βi,γi に基づいて算出する一方、ブラックの網点面積率が１００％である辺上における、任意の網点面積率での分光反射率を、事前に指定された色票の分光反射率に基づいて算出する。
【００３７】
また、分光反射率算出部４６は、任意の位置に対応する分光反射率を算出する際に、係数算出部４４で算出された係数ξi,ηi を用いるほか、ブラックグラデーションを表す辺上における任意の網点面積率に対応する分光反射率を、係数算出部４４で算出された係数ξi,ηi に基づいて算出する。
さらに、分光反射率算出部４６は、事前に指定された色票の分光反射率に基づいて算出・決定される一定値（例えば０．０３５）を、ブラックの網点面積率が１００％である辺上における任意の網点面積率での分光反射率として出力することもできる。
【００３８】
またさらに、分光反射率算出部４６は、図７および図８を参照しながら後述するごとく、ｎ次元立体から、カラープリンタで最後に刷る第４色材（ブラック）の網点面積率を一定化した３次元立体を算出し、この３次元立体から、カラープリンタで３番目に刷る第３色材（本実施形態ではシアン）の網点面積率を一定化した２次元色平面を算出し、最終的に、この２次元色平面から、係数ξi,ηi を用いて、任意の網点面積率に設定された４色の色材を組み合わせた場合の分光反射率を算出する。
【００３９】
さらにまた、本実施形態の分光反射率算出部４６は、後述するごとく、特定の網点面積率に設定された４色の色材を組み合わせた場合の分光反射率に、特定の値を指定すべく、その特定の網点面積率を与える色信号と前記任意の網点面積率を与える色信号との距離の逆数に基づいて重み係数を算出し、その重み係数を特定の値に乗算した値を、算出された分光反射率に加算するとともに、前記距離が０の場合には、分光反射率として特定の値を出力する機能もそなえている。
【００４０】
そして、変換部４８は、分光反射率算出部４６で算出された分光反射率を、モニタ３６用の色信号（例えばＲＧＢ信号）に変換し、その色信号に基づく色表示をモニタ３６に行なわせるものである。
次に、図２〜図１３を参照しながら、本実施形態の色再現方法について詳細に説明する。
【００４１】
本実施形態においては、任意数ｎの色材の組み合わせから選定されたｎ・２(n-1) 本の辺から得られる分光反射率を使用する。これらの辺を表現するために、以下のような規定▲１▼〜▲３▼を設ける。
▲１▼各辺において、色信号の採る値（即ち、網点面積率）は、(a) ０％、(b) １００％、および、(c) ０％〜１００％の間における複数（本実施形態の１１）の数値の組、のうちのいずれか一つに限り、且つ、(c) の複数の数値の組を採る色信号は１色のみとする。即ち、これらの辺は、モニタ３６上で色を再現されるべきカラープリンタの色域の縁辺を表すものとする。
【００４２】
▲２▼本実施形態では、ｎ＝４とし、例えば、ブラックを含む任意の４つの色材をシアン，マゼンタ，イエロおよびブラックとする。ただし、これ以外の数や色材であってもよいものとする。そして、各辺は、色符号Ｃ，Ｍ，Ｙ，Ｋを用いて表記するものとし、シアン１００％を“Ｃ”、マゼンタ１００％を“Ｍ”、イエロ１００％を“Ｙ”、ブラック１００％を“Ｋ”として表記し、０％ならば、その色符号の表記を行なわない。色符号の表記順序は刷り順通りとするが、０％〜１００％の間における複数の数値の組をとる色信号の色符号については、一番最後に表記する。なお、本実施形態では、刷り順は“ＹＭＣＫ”とするが、これ以外であってもよい。ただし、ブラックだけは最後に刷られるものとする。
【００４３】
▲３▼上記▲１▼，▲２▼の規定に従った場合、例えば下記表１のような数値の組み合わせの辺は、色符号により“ＹＫＣ”として表記されることになる。つまり、イエロおよびブラックの網点面積率が１００％、マゼンタの網点面積率が０％、シアンの網点面積率が、０％〜１００％の間で１０％刻みの値をとる。
【００４４】
【表１】
なお、前述した通り、各辺でのデータの個数をＮで表し、本実施形態では、例えばＮ＝１１とする。
上述のような表記に従い、以下に、本実施形態の色再現方法について詳細に説明する。
〔１〕まず、本実施形態の判別部４２は、カラープリンタ出力の色再現処理を開始するに当たって、一度だけ、下記 (1-1)〜(1-5) のような判定工程を実施する。
【００４６】
(1-1) ブラックの網点面積率が１００％である辺（即ち、ＫＣ，ＫＭ，ＫＹ，ＣＫＭ，ＭＫＣ，ＣＫＹ，ＹＫＣ，ＹＫＭ，ＭＣＫＹ，ＹＣＫＭ，ＹＭＫＣの１２辺）のそれぞれにおいて、事前に指定された色票毎に、可視波長域に含まれる波長についての分光反射率を平均した値を求める。例えば、可視波長域を 380nm〜730nm とした場合、10nm毎の分光反射率（つまり、波長380nm,390nm,400nm,...,730nm の分光反射率；全部で３６の分光反射率）の総和を３６で除算した値を、分光反射率の平均値として求める。このようにして、分光反射率の平均値を、各辺上における１１個の色票のそれぞれについて算出する。なお、可視波長域は 380nm〜730nm に限定されるものではない。
【００４７】
(1-2) 各辺の１１個の色票について算出された分光反射率の平均値の標準偏差を算出する。
(1-3) (1-2) で算出された標準偏差と、予め与えられた所定値（例えば０．１）とを比較する。
(1-4) (1-3) で比較した結果、標準偏差が所定値よりも大きい場合、モニタ３６で再現すべき色情報には、以下に説明する本実施形態の色再現方法は適用されない。
【００４８】
(1-5) (1-3) で比較した結果、標準偏差が所定値以下であれば、今回、モニタ３６で再現すべき色情報がカラープリンタ用のものであると判別され、その色情報には、以下に説明する本実施形態の色再現方法が適用される。つまり、前述したように、係数算出部４４，分光反射率算出部４６および変換部４８が動作する。
【００４９】
ここで、ブラックを最後に刷るカラープリンタにおいて、ブラックの網点面積率が１００％である場合には、例えば図１２に示すように、その他の色シアン，マゼンタ，イエロの網点面積率がどのような値を採ろうとも、波長によらず、分光反射率はほぼ一定になる。
カラープリンタで用いられるトナーは、前述した通り、固体であり、元々、光透過性を有していない。従って、カラープリンタによる印刷結果では、最後に刷られたトナーの色、つまり、一番上に載るトナーの色が見えることになり、その下に何色があるかは色の見え方に関係なくなる。カラープリンタはこのような特徴を有しているので、ブラックの網点面積率が１００％である場合、他の色に関係なく、全ての色票の分光反射率の値は、ほぼ等しくなり、当然、前述した標準偏差も小さくなる。本実施形態の判別部４２では、このような特徴を利用して、モニタ３６で再現すべき色情報がカラープリンタ用のものであるか否かを判別しているのである。
【００５０】
なお、図１２は、ブラック，シアンおよびマゼンタの網点面積率が１００％でイエロの網点面積率が０％〜１００％となるとき（４次色）の分光反射率の測定結果を示すグラフである。この図１２を見て明らかなように、ブラックの網点面積率が１００％である場合、非常に高い精度で、分光反射率を一定値に置き換えることが可能であり、その一定値として、図１２に示す例では“０．０３５”を採用することができる。
【００５１】
〔２〕上述のようにして、カラープリンタ出力の判定を行なった後、本実施形態の係数算出部４４は、一度だけ、下記 (2-1),(2-2)のような係数算出工程を実施する。
(2-1) ブラックの網点面積率が１００％である辺（前記１２本の辺）とブラックグラデーションを表す辺（“Ｋ”として表記される辺）とを除く辺は、下記の通り表記される１９本の辺である。
【００５２】
Ｃ，Ｍ，Ｙ，ＣＭ，ＭＹ，ＣＹ，ＭＣ，ＹＭ，ＹＣ，ＣＫ，ＭＫ，ＹＫ，ＭＣＹ，ＹＣＭ，ＹＭＣ，ＭＣＫ，ＹＣＫ，ＹＭＫ，ＹＭＣＫ
これらの辺の各々に対して、Ｎ個（本実施形態では１１個）の色票の分光反射率が各色信号に対応して存在する。そして、本実施形態の係数算出部４４では、上記１９辺におけるＮ個の分光反射率に対して、下式（１）および（２）を使用して最小２乗法により、各辺のＮ個の色票に対応する係数α，β，γおよび係数ξ，ηが算出される。
【００５３】
ここで、下式（１）および（２）の説明をする前に、係数α，β，γおよび係数ξ，ηのもつ物理的な意味について説明する。
(2-1-1) 係数α，β，γの物理的意味について
図２に示すように、カラープリンタで印刷されるべき色材が紙に載っている場合、目に入る光は、矢印Ａ１，Ａ２，Ａ３で示す３種類の光のいずれかに分類される。
【００５４】
矢印Ａ１は、入射時および反射時の両方とも色材を通過する光（物理的ドットゲイン）を示し、係数γは、このような光の量を表している。また、このような光の影響は、色材の網点面積率が１００％のときに最も顕著になる。
矢印Ａ２は、入射時および反射時の両方とも色材を通過しない光を示し、係数αは、このような光の量を表している。また、このような光の影響は、色材の網点面積率が０％のときに最も顕著になる。
【００５５】
矢印Ａ３は、入射時および反射時のうちの一方のみに色材を通過する光（光学的ドットゲイン）を示し、係数βは、このような光の量を表している。また、このような光の影響は、色材の網点面積率が５０％のときに最も顕著になる。
実際の印刷結果とモニタで再現される色とが、ノイゲバウアの式を用いた場合には一致せず、本発明の方法を用いた場合に一致する理由がここにある。即ち、ノイゲバウアの式では、矢印Ａ１，Ａ２で示すような光は考慮しているが、矢印Ａ３で示すような光の効果（色材が載っている部分の縁辺部において生じる、紙およびベタ色票では再現できない光学現象）を無視している。これに対し、本発明の方法においては、矢印Ａ３で示すような光の効果を考慮し、後述する仮想１次色についての分光反射率を正確に再現できるようにするために、係数α，β，γを使用している。
【００５６】
(2-1-2) 係数ξ，ηの物理的意味について
係数ξ，ηは、前述した係数α，γに対応するものであり、係数ηは、入射時および反射時の両方とも色材を通過する光の量を表し、係数ξは、入射時および反射時の両方とも色材を通過しない光の量を表している。
ここでは、色材の載っていない面積を係数ξで、載っている面積をηで表している。これにより、従来法〔例えばMurray-Davies(マーレイ・デービス）の式〕よりも正確に、色材の載っている面積を把握することが可能になる。なお、これは前記の１９辺の全てに適用可能であり、１次色しか適用できなかった従来法とは大きく異なる点である。即ち、本実施形態では、１次から４次色に至る全ての基本的な色の組み合わせに対して、従来法よりも正確に色材の載っている面積を算出することにより、推定精度を向上させているのである。
【００５７】
【数１】
上式（１），（２）において、
λ：波長(nm)〔入力情報〕
Ｓi （λ）：ｉは、各辺上での網点面積率０％〜１００％におけるＮ個の数値（本実施形態では０，１０，２０，３０，４０，５０，６０，７０，８０，１００の１１個）で、個々の色票の値（網点面積率）を表し、Ｓi （λ）はその色票の分光分布〔出力情報〕
Ｓpaper （λ）〔＝Ｓ0%（λ）〕：紙白の分光反射率〔入力情報〕
Ｓ50% （λ）：網点面積率５０％に対応する色材の載った色票の分光反射率〔入力情報〕
Ｓ100%（λ）：網点面積率１００％に対応する色材の載った色票の分光反射率〔入力情報〕
αi ，βi ，γi ：網点面積率ｉ％の色票に対応する係数α，β，γ〔出力情報〕
ξi ，ηi ：網点面積率ｉ％の色票に対応する係数ξ，η〔出力情報〕
である。なお、〔入力情報〕と記したデータは、上式（１），（２）を用いて最小二乗法を行なう際に、上式（１），（２）に入力されるデータであることを示し、〔出力情報〕と記したデータは、最小二乗法により算出・出力されるデータであることを示している。
【００５９】
(2-2) ブラックグラデーションを表す辺（Ｋ）に属するＮ個（本実施形態では１１個）の色票に対しては、図１１を見れば明らかなように、データの線形性のために最小２乗法を適用することができない。
そこで、本実施形態では、ドットゲインを表す係数ξ，ηを以下のように算出し、分光反射率の再現およびブラックの載っている面積を正確に把握するために、係数ξ，ηを使用する。なお、判別部４２により本発明の手法が適用可能であると判断された場合には、ブラックの透過率が非常に小さいため、図２に矢印Ａ３で示すようなの光は無視することができる。
【００６０】
さて、係数ξ，ηを算出する際には、まず、ブラックグラデーションを表す辺（Ｋ）上のＮ個の色票の分光反射率を、図１１に示すようにグラフ化してから、グラフの平坦部を決定して、その平坦部での分光反射率の平均値を算出する。平均化する範囲（平坦部）は、可視波長域すべてに亘ってもよいし、その可視波長域の一部に亘ってもよい。
【００６１】
そして、Ｎ個の色票の中で、ブラックの網点面積率がｉ％の色票についての平均値をＳi とすると（０≦ｉ≦１００）、下記の連立方程式（３）を解くことにより、その色票に対応する係数ξi ，ηi が算出される。
【００６２】
【数２】
ここで、Ｓpaper は紙白の分光反射率の波長域における平均値〔入力情報〕、Ｓ100%は網点面積率１００％に対応する色材の載った色票の分光反射率の波長域における平均値〔入力情報〕である。
〔３〕上述のようにして係数算出工程を実施し、係数αi ，βi ，γi ，ξi およびηi を算出した後、本実施形態の分光反射率算出部４６は、下記 (3-1)〜(3-3) の工程を実施し、分光反射率の算出・推定を行なう。
【００６４】
以下では、３２辺の各々において、事前に指定された色票以外の色信号（網点面積率）に対応する分光反射率を算出・推定する方法について、項目(3-1) で説明する。なお、以降の説明では、３２辺のそれぞれのことを「仮想１次色」と呼称する。
そして、前記仮想１次色から「仮想２次色」を算出・推定する方法について、項目(3-2) で説明する。「仮想２次色」とは、シアン，マゼンタ，イエロおよびブラックの色信号のうちの２つの網点面積率が任意で、残りの２つの網点面積率が０％または１００％で固定されている色のことを意味する。
【００６５】
その後、項目 (3-1)および(3-2) で説明された方法を用いて、シアン，マゼンタ，イエロおよびブラックの色信号で表される４次元色立体から、順次次元を落としていき、最後に項目(3-2) の方法を実施することにより、任意のシアン，マゼンタ，イエロおよびブラックの色信号に対応する分光反射率を算出．推定できることを項目(3-3) で説明する。
【００６６】
(3-1) 仮想１次色の算出・推定方法
仮想１次色の算出・推定方法に際しては、辺毎に、上述した工程〔２〕で算出された以下のデータを使用する。
▲１▼ブラックの値（網点面積率）が１００％に該当する辺およびブラックグラデーションを表す辺以外の辺
（ａ）Ｎ個の色票におけるαi,βi,γi の値
（ｂ）Ｎ個の色票におけるξi,ηi の値
（ｃ）色信号値（網点面積率）が０，５０，１００％における色票の分光反射率Ｓpaper （λ），Ｓ50% （λ），Ｓ100%（λ）
（ｄ）Ｎ個の色票におけるｉの値（ｉは０％〜１００％の間における複数の数値を表す）
▲２▼ブラックグラデーションを表す辺
（ａ）Ｎ個の色票におけるξi,ηi の値
（ｂ）色信号値が０，１００％における色票の分光反射率Ｓpaper （λ），Ｓ100%（λ）
ｃ）Ｎ個の色票におけるｉの値（ｉは０％〜１００％の間における複数の数値を表す）
▲３▼ブラックの値（網点面積率）が１００％に該当する辺
（ａ）Ｎ個の色票の分光反射率Ｓi （λ）
（ｂ）Ｎ個の色票におけるｉの値（ｉは０％〜１００％の間における複数の数値を表す）
これらのデータを使用して、仮想１次色の分光反射率を算出・推定する方法を以下に記す。
【００６７】
(3-1-1) 以下においては、事前にイエロ，マゼンタ，シアンおよびブラックの色信号（即ちＹ，Ｍ，Ｃ，Ｋの網点面積率）が、％単位で与えられているものとする。ただし、ここでは仮想１次色を扱うので、先に定義した通り、４色のうちのの１色だけを、０％〜１００％の範囲で任意に指定することができ、残りの３色は０％か１００％のいずれか一方であるとする。以後、この任意に指定可能な数値（網点面積率）をｊ（％）で表す。
【００６８】
(3-1-2) この任意の網点面積率ｊ％の色信号に対応する分光反射率Ｓj （λ）を推定するために、式（１）や（２）を先程とは逆に使用する。即ち、
【００６９】
【数３】
ここで、上式（４）は、ブラックの値が１００％に該当する辺およびブラックグラデーションを表す辺以外の辺に適用される。また、上式（５）は、ブラックグラデーションを表す辺に適用される。
各係数αj,βj,γj は、図３に示すように、網点面積率ｊに依存する。分光反射率Ｓpaper （λ）に対する係数αj の値は、網点面積率ｊが０％のときには１であり、５０％および１００％のときには０である。同様に、分光反射率Ｓ50% （λ）に対する係数βj の値は、網点面積率ｊが５０％のときには１であり、０％および１００％のときには０であり、分光反射率Ｓ100%（λ）に対する係数γj の値は、網点面積率ｊが１００％のときには１であり、０％および５０％のときには０である。
【００７１】
また、係数αj は、工程〔２〕で求めたＮ個の係数αi （本実施形態は、ｉ＝0,10,20,30,40,50,60,70,80,90,100の１１個）によって代表されている。同様に、係数βj は、工程〔２〕で求めたＮ個の係数βi によって代表され、係数γj は、工程〔２〕で求めたＮ個の係数γi によって代表されている。
そして、任意の網点面積率ｊにおける係数αj,βj,γj の値は、それぞれ、上述したＮ個の係数αi,βi,γi の間を補間することによって、例えば下式 (6-1)〜(6-3) のごとく決定される。
【００７２】
各係数ξj,ηj も、図４に示すように、網点面積率ｊに依存する。分光反射率Ｓpaper （λ）に対する係数ξj の値は、網点面積率ｊが０％のときには１であり、１００％のときには０である。同様に、分光反射率Ｓ100%（λ）に対する係数ηj の値は、網点面積率ｊが１００％のときには１であり、０％のときには０である。
【００７３】
また、係数ξj は、工程〔２〕で求めたＮ個の係数ξi （本実施形態は、ｉ＝0,10,20,30,40,50,60,70,80,90,100の１１個）によって代表されている。同様に、係数ηj は、工程〔２〕で求めたＮ個の係数ηi によって代表されている。
そして、任意の網点面積率ｊにおける係数ξj,ηj の値は、それぞれ、上述したＮ個の係数ξi,ηi の間を補間することによって、例えば下式(7-1),(7-2) のごとく決定される。
【００７４】
つまり、Ｎ個の色票に対応する、先に述べたｉの値のうちで、ｊに近接しｊをはさむ２つの値ｋとｌとを算出し、下式 (6-1)〜(6-3) および(7-1),(7-2) に代入することにより、係数αj,βj,γj およびξj,ηj を算出する。ここで、ｋ≦ｊ≦ｌとする。そして、下式 (6-1)〜(6-3) および(7-1),(7-2) によってそれぞれ算出された係数αj,βj,γj およびξj,ηj を、それぞれ上式（３）および（４）に代入することで、任意の網点面積率ｊに対応する分光反射率Ｓj （λ）が算出される。
【００７５】
【数４】
上述した式（４),（５), (6-1)〜(6-3) および(7-1),(7-2) を用いて算出・推定された分光反射率を、実測された分光反射率とともに、図９〜図１１に示す。これらの図９〜図１１において、太線が、実測された分光反射率を示し、細線が、本実施形態の方法により算出・推定された分光反射率を示す。
なお、図９は、マゼンタ１次色の分光反射率を、係数α，β，γ〔式（４）＆ (6-1)〜(6-3) 〕を用いて算出・推定した結果を示すグラフであり、図１０は、マゼンタ１次色の分光反射率を、係数ξ，η〔式（５）＆(7-1),(7-2) 〕を用いて算出・推定した結果を示すグラフであり、図１１は、ブラック１次色の分光反射率を、係数ξ，η〔式（５）＆(7-1),(7-2) 〕を用いて推定した結果を示すグラフである。
【００７７】
図９に示すように、係数α，β，γ〔式（４）＆ (6-1)〜(6-3) 〕を用いることにより、極めて高い精度で仮想１次色（マゼンタ１次色）の分光反射率を算出・推定することができる。
また、図１０に示すように、係数ξ，η〔式（５）＆(7-1),(7-2) 〕を用いても、ある程度高い精度で仮想１次色（マゼンタ１次色）の分光反射率を算出・推定することができる。
【００７８】
さらに、図１１に示すように、係数ξ，η〔式（５）＆(7-1),(7-2) 〕を用いることにより、かなり高い精度でブラックグラデーションを表す辺（ブラック１次色）上の任意の位置での分光反射率を算出・推定することができる。
実際に本実施形態の方法を用いて仮想１次色の分光反射率を算出・推定する際には、係数ξ，η〔式（５）＆(7-1),(7-2) 〕ではなく、係数α，β，γ〔式（４）＆ (6-1)〜(6-3) 〕を用いる。ただし、式(7-1),(7-2) により算出される係数ξj ，ηj は、仮想２次色の分光反射率を算出・推定する際に、各辺上の位置を表す数値（ドットゲイン）として使用される。
【００７９】
(3-1-3) ブラックの値（網点面積率）が１００％である辺上における、任意の位置の分光反射率Ｓj （λ）を算出・推定する際には、下式（８）を用いる。即ち、前述した２つの値ｋ，ｌと、既知の分光反射率Ｓk （λ）およびＳl （λ）とを下式（８）に代入し、分光反射率Ｓk （λ）およびＳl （λ）を直接的に比例配分することにより、分光反射率Ｓj （λ）を算出する。
【００８０】
【数５】
(3-2) 仮想２次色の算出・推定方法
仮想２次色の算出・推定方法に際し、本実施形態の分光反射率算出部４６は、項目(3-1) において前述した仮想１次色の算出・推定方法を用い、下記 (3-2-1)〜(3-2-7) の工程を実施する。
(3-2-1) ここでは、色材Ａおよび色材Ｂに関する仮想２次色を算出する方法について説明する。このとき、色材Ａの網点面積率がｘ％で色材Ｂの網点面積率がｙ％である仮想２次色を推定すべき色であるとし、２次元色平面上でのその色の位置を点（ｘ，ｙ）で、またその位置における分光反射率をＳx,y （λ）で表す。なお、点（ｘ，ｙ）はその色の位置の呼称とし、実際の位置を示すものではない。
【００８２】
(3-2-2) 点（ｘ，ｙ）を含む平面を構成する４つの辺を決定する。ここで、各々の辺と、図５に示すような平面とを決定すべく、以下のデータ（ａ）〜（ｄ）を使用する。下記データ（ａ）〜（ｄ）は、すべて、項目(3-1) において前述した方法で算出可能である。なお、図５は、本実施形態における、仮想２次色の算出・推定に必要なデータの位置関係を説明するための図である。
【００８３】
（ａ）（ｉ，０）で表され０≦ｉ≦１００である“色材Ａ辺”上における点（ｘ，０）に対応する分光反射率ＳX,0 （λ）と係数ξx,o ，ηx,0
（ｂ）（ｉ，１００）で表され０≦ｉ≦１００である“色材Ａ，Ｂ１００％辺”上における点（ｘ，１００）に対応する分光反射率ＳX,100 （λ）と係数ξx,100 ，ηx,100 （この“色材Ａ，Ｂ１００％辺”は、図５に示すように、“色材Ａ辺”と対向する）
（ｃ）（０，ｊ）で表され０≦ｊ≦１００である“色材Ｂ辺”上における点（０，ｙ）に対応する分光反射率Ｓo,y （λ）と係数ξo,y ，η0,y
（ｄ）（１００，ｊ）で表され０≦ｊ≦１００である“色材Ａ１００％，Ｂ辺”上における点（１００，ｙ）に対応する分光反射率Ｓ100,y （λ）と係数ξ100,y ，η100,y （この“色材Ａ１００％，Ｂ辺”は、図５に示すように、“色材Ｂ辺”と対向する）
ただし、上述した４辺のなかに、ブラックの値（網点面積率）が１００％である辺が含まれる場合には、その辺と対向する、ブラックの値（網点面積率）が０％である辺の係数ξ，ηをそのまま採用する。
【００８４】
(3-2-3) 前記４辺は、図６に示すように配置される。その配置方法について、以下に説明する。なお、図６は、本実施形態における、仮想２次色の分光反射率の算出・推定方法を説明するための図である。
（ａ）（０，０）を原点（Ｐ）とする。
（ｂ）“色材Ａ辺”および“色材Ｂ辺”を、それぞれ上記原点から水平，垂直に配置する。この際、“色材Ａ辺”の長さ，点（ｘ，０）の位置，“色材Ｂ辺”の長さ，点（０，ｙ）の位置を、以下のように設定する。
【００８５】
“色材Ａ辺”の長さ…ξx,0 ＋ηx,0
点（ｘ，０）の位置…“色材Ａ辺”をηx,0 ：ξx,0 の比で内分する位置
“色材Ｂ辺”の長さ…ξ0,y ＋η0,y
点（０，ｙ）の位置…“色材Ａ辺”をη0,y ：ξ0,y の比で内分する位置
（ｃ）“色材Ａ，Ｂ１００％辺”および“色材Ａ１００％，Ｂ辺”については、以下のように配置する。なお、これらの辺の長さ、および、各辺上における点の算出方法は、前記（ｂ）と同様であるので、その説明は省略する。
【００８６】
図６に示すように、点（１００，０）から長さ（ξ100,y ＋η100,y ）を半径とする円と、点（０，１００）から長さ（ξx,100 ＋ηx,100 ）を半径とする円とを描く。このとき、２つの円は、２点で交わることになる。これらの２つの交点のうち、“色材Ａ辺”と“色材Ｂ辺”とで形成される象限と同一に存在し、且つ、図６に示す平面上の座標値（１００，１００）の位置からの距離が近いものを選択する。そして、“色材Ａ，Ｂ１００％辺”および“色材Ａ１００％，Ｂ辺”は、選択された交点で交わるように配置される。
【００８７】
(3-2-4) 求めるべき分光反射率Ｓx,y （λ）に対応する点（ｘ，ｙ）の位置を、以下のようにして決定する。つまり、点（ｘ，０）と点（ｘ，１００）とを直線で結ぶとともに、点（０，ｙ）と点（１００，ｙ）とを直線で結び、これら２つの直線の交点を点（ｘ，ｙ）の位置とする。
(3-2-5) 各点間の距離を、以下のように定義する。ここで、距離とは図６に示す２次元色平面におけるユークリッド距離のことを意味する。
【００８８】
点（ｘ，ｙ）と点（ｘ，０）との間の距離…ｌ1
点（ｘ，ｙ）と点（０，ｙ）との間の距離…ｌ2
点（ｘ，ｙ）と点（ｘ，１００）との間の距離…ｌ3
点（ｘ，ｙ）と点（１００，ｙ）との間の距離…ｌ4
(3-2-6) 重み係数ｗi （ｉ＝１，２，３，４）を、下式 (9-1)〜(9-4) もしくは(10-1)〜(10-4)のように定義する。なお、前者の式 (9-1)〜(9-4) と後者の式(10-1)〜(10-4)とを混合して使用してはならない。即ち、ｗ1 で前者の式(9-1) を使用する場合には、ｗ2,ｗ3,ｗ4 でも前者の式 (9-2)〜(9-4) を使用する。当然、ｗ1 で後者の式(10-1)を使用する場合には、ｗ2,ｗ3,ｗ4 でも後者の式(10-2)〜(10-4)を使用する。
【００８９】
【数６】
(3-2-7) 分光反射率Ｓx,y （λ）は、上述のごとく算出された分光反射率Ｓx,0 （λ），Ｓ0,y （λ），Ｓx,100 （λ），Ｓ100,y （λ）および重み係数ｗ1,ｗ2,ｗ3,ｗ4 を下式（１１）に代入することにより算出される。
【００９１】
【数７】
(3-3) ４次までの色の分光反射率の算出・推定方法
４次色の算出・推定方法に際し、本実施形態の分光反射率算出部４６は、項目(3-1) において前述した仮想１次色の算出・推定方法と、項目(3-2) において前述した仮想２次色の算出・推定方法とを用い、下記 (3-3-1)〜(3-3-6) の工程を実施する。その工程について、図６および図７を参照しながら説明する。
【００９３】
なお、図７は、本実施形態における、次元数リダクション処理を説明するための図、図８は、本実施形態における、ＹＭＣ３次色の分光反射率の算出・推定方法を説明するための図である。
また、ここでは、事前にイエロ，マゼンタ，シアンおよびブラックの色信号により各色の網点面積率（Ｙ，Ｍ，Ｃ，Ｋ）が％単位で与えられているものとする。ただし、０≦Ｙ，Ｍ，Ｃ，Ｋ≦１００とする。
【００９４】
ところで、光透過性を有する色材を用いた印刷では、ブラックの網点面積率が１００％である時でも３次元色立体を構成されていた。しかし、ブラックの網点面積率が１００％である場合、カラープリンタではブラックが最後に重ねられるため、イエロ，マゼンタ，シアンがどのように印刷されていても、印刷結果は全て同じ色になってしまう。従って、カラープリンタでは、ブラック１００％の３次元色立体は、図７に示すように、ほぼ一点に縮退することになる。このようなカラープリンタの特異性を利用することにより、Ｙ，Ｍ，Ｃ，Ｋからなる４次元色立体の次元数のリダクションを以下のような簡素な手法で行なうことができる。なお、図７において、“Composite black ”とは、３つの色材Ｙ，Ｍ，Ｃを等量混合することによって描かれるブラックのことを言い、１つの色材Ｋで描かれるブラックと区別している。
【００９５】
(3-3-1) まず、図７に示す１２個の四辺形から、項目(3-2) において前述した仮想２次色の算出・推定方法を使用して、最後に刷るブラックの網点面積率が一定値Ｋ％である１２辺を算出することにより、図７に示すように、４次元色立体を、１２辺から成る、ブラックＫ％一定の３次元色立体にリデュース(reduce)する。
【００９６】
（ａ）Ｋ，Ｃ，ＣＫ，ＫＣの辺で構成される四辺形から、ブラックＫ％一定で、シアンの値を可変とする辺１を算出する。
（ｂ）Ｋ，Ｍ，ＭＫ，ＫＭの辺で構成される四辺形から、ブラックＫ％一定で、マゼンタの値を可変とする辺２を算出する。
（ｃ）Ｋ，Ｙ，ＹＫ，ＫＹの辺で構成される四辺形から、ブラックＫ％一定で、イエロの値を可変とする辺３を算出する。
【００９７】
（ｄ）ＣＫ，ＣＭ，ＭＣＫ，ＣＫＭの辺で構成される四辺形から、ブラックＫ％一定で且つシアン１００％一定で、マゼンタの値を可変とする辺４を算出する。
（ｅ）ＭＫ，ＭＣ，ＭＣＫ，ＭＫＣの辺で構成される四辺形から、ブラックＫ％一定で且つマゼンタ１００％一定で、シアンの値を可変とする辺５を算出する。
【００９８】
（ｆ）ＣＫ，ＣＹ，ＹＣＫ，ＣＫＹの辺で構成される四辺形から、ブラックＫ％一定で且つシアン１００％一定で、イエロの値を可変とする辺６を算出する。
（ｇ）ＭＣＫ，ＭＣＹ，ＹＭＣＫ，ＭＣＫＹの辺で構成される四辺形から、ブラックＫ％一定，マゼンタ１００％一定かつシアン１００％一定で、イエロの値を可変とする辺７を算出する。
【００９９】
（ｈ）ＭＫ，ＭＹ，ＹＭＫ，ＭＫＹの辺で構成される四辺形から、ブラックＫ％一定で且つマゼンタ１００％一定で、イエロの値を可変とする辺８を算出する。
（ｉ）ＹＫ，ＹＣ，ＹＣＫ，ＹＫＣの辺で構成される四辺形から、ブラックＫ％一定で且つイエロ１００％一定で、シアンの値を可変とする辺９を算出する。
【０１００】
（ｊ）ＹＣＫ，ＹＣＭ，ＹＭＣＫ，ＹＣＫＭの辺で構成される四辺形から、ブラックＫ％一定，シアン１００％一定およびイエロ１００％一定で、マゼンタの値を可変とする辺１０を算出する。
（ｋ）ＹＭＫ，ＹＭＣ，ＹＭＣＫ，ＹＭＫＣの辺で構成される四辺形から、ブラックＫ％一定，マゼンタ１００％一定およびイエロ１００％一定で、シアンの値を可変とする辺１１を算出する。
【０１０１】
（ｌ）ＹＫ，ＹＭ，ＹＭＫ，ＹＫＭの辺で構成される四辺形から、ブラックＫ％一定で且つイエロ１００％一定で、マゼンタの値を可変とする辺１２を算出する。
(3-3-2) 次に、以下に示す４個の四辺形から、項目(3-2) において前述した仮想２次色の算出・推定方法を使用して、３番目に刷るシアンの網点面積率が一定値Ｃ％である４辺を算出することにより、図７および図８に示すように、３次元色立体を、４辺から成る、ブラックＫ％一定で且つシアンＣ％一定の２次元色平面（仮想２次平面）にリダクションする。つまり、ブラックＫ％一定で且つシアンＣ％一定に対応する辺を各々算出し、任意のマゼンタＭ％および任意のイエロＹ％の組み合わせに対する色空間を表す１つの四辺形（２次元色平面）を得る。
【０１０２】
（ａ）上述のごとく算出された４つの辺１，６，９，３で構成される四辺形から、イエログラデーションを表す辺▲１▼を算出する。
（ｂ）上述のごとく算出された４つの辺９，１０，１１，１２で構成される四辺形から、イエロ１００％一定で且つマゼンタグラデーションを表す辺▲２▼を算出する。
【０１０３】
（ｃ）上述のごとく算出された４つの辺５，８，１１，７で構成される四辺形から、マゼンタ１００％一定で且つイエログラデーションを表す辺▲３▼を算出する。
（ｄ）上述のごとく算出された４つの辺１，４，５，２で構成される四辺形から、マゼンタグラデーションを表す辺▲４▼を算出する。
【０１０４】
なお、図８では、ブラック０％一定のＹＭＣ３次色についての３次元色立体から、シアンＣ％一定でスライスされた仮想２次平面（２次元色平面）が示されている。
(3-3-3) 以上のようにして算出された４つの辺▲１▼〜▲４▼で構成される四辺形から、項目(3-2) において前述した仮想２次色の算出・推定方法を使用して、所望の色信号（Ｋ，Ｃ，Ｍ，Ｙ；Target Color）に対応する分光反射率が、算出・推定される。
【０１０５】
(3-3-4) １次色の分光反射率は、項目(3-1) において前述した仮想１次色の算出・推定方法を使用して、辺Ｋ，Ｃ，Ｍ，Ｙから算出・推定される。
(3-3-5) ２次色の分光反射率は、項目(3-2) において前述した仮想２次色の算出・推定方法を使用して、算出・推定される。
(3-3-6) ３次色の分光反射率は、例えば図８に示すようにブラックＫ％の辺のみから３次元色立体を構成してから、前述した項目(3-3-2) 以降の工程を実施することにより、算出推定される。
【０１０６】
なお、一度、前述した式（４),（５), (6-1)〜(6-3) および(7-1),(7-2) を用いて分光反射率を算出しその分光反射率を式（１１）に代入した後は、式（１）〜（３）を用いて算出された係数α，β，γおよび係数ξ，ηは、精度上、使用できなくなる。それ故、これ以後は、新たに算出された辺については、式（８）を利用して、分光反射率の算出を行なう。
【０１０７】
〔４〕ところで、本実施形態の分光反射率算出部４６は、必要に応じて、決め打ちと呼ばれる工程を下記のように実施する。この工程は、特定の網点面積率（Ｙ，Ｍ，Ｃ，Ｋ）に設定された４色の色材を組み合わせた場合の分光反射率を、特定の値に合わせる必要がある場合（合わせるべきデータが存在する場合）にのみ実施される。それゆえ、必要ないと判断した場合（合わせるべきデータが存在しない場合）には、その工程は実施されない。
【０１０８】
なお、この工程は、項目(3-2) で前述した仮想２次色の算出・推定方法において、重み係数ｗi （ｉ＝１，２，３，４）として式(10-1)〜(10-4)を使用する場合に適用される。以下、その工程について説明する。
(4-1) 合わせるべき色信号の値およびその分光反射率を、以下に定義する。
▲１▼合わせるべき色信号…（Ｋm,Ｃm,Ｍm,Ｙm ），ｍ＝１，２,..,Ｍ
▲２▼合わせるべき分光反射率…Ｓm （λ），ｍ＝１，２,..,Ｍ
(4-2) 本工程は、以下のとおりに実施する。
【０１０９】
(4-2-1) 合わせるべき色信号が１次色の場合には、項目(3-1) で前述した仮想１次色の算出・推定方法において、推定に使用されたＮ個の色票の分光反射率に、さらに、合わせるべき分光反射率Ｓm （λ）を追加し、項目(3-1) で説明した通りの方法を実施すればよい。データの数が増えるだけで、基本的に行なわれる処理は、項目(3-1) で説明した方法と何ら変わらない。
【０１１０】
(4-2-2) 合わせるべき色信号が２次色，３次色，４次色の場合には、以下の処理（ａ）〜（ｄ）を行なう。
（ａ）以下に説明するデータ合わせのための処理は、前記式（１１）により、最終的に分光反射率を算出する際にのみ実施される。
（ｂ）合わせるべき色信号（Ｋm,Ｃm,Ｍm,Ｙm ）と事前に与えられた色信号（Ｋ，Ｃ，Ｍ，Ｙ）との距離を、下式（１２）のように定義する。
【０１１１】
【数８】
（ｃ）このように定義された距離ｌm+4 に対応する重み係数ｗk を、下式（１３）で与える。
【０１１３】
【数９】
（ｄ）上述のごとく算出された重み係数ｗk と、合わせるべき分光反射率Ｓm （λ）とから、求める分光反射率Ｓx,y （λ）は、下式（１４）により算出されることになる。なお、式（１２）において、ｌm+4 ＝０となる場合（つまり、合わせるべき色信号Ｋm,Ｃm,Ｍm,Ｙm と事前に与えられた色信号Ｋ，Ｃ，Ｍ，Ｙとが一致している場合）、下式（１５）のように、求める分光反射率Ｓx,y （λ）を合わせるべき分光反射率Ｓm （λ）とする。
【０１１５】
【数１０】
このようにして、決め打ちと呼ばれる工程を実行することにより、特定の色信号（網点面積率）の組み合わせに対して特定の分光反射率を指定することができる。従って、ページ内誤差により不一致となる色を予め測定しておけば、その色については所定の分光反射率を強制的に算出するようにすることができ、ページ内誤差を無視できるので、高精細な推定を実施することが可能になる。
【０１１７】
〔５〕上述のようにして分光反射率算出部４６により分光反射率を算出した後、本実施形態の変換部４８は、分光反射率算出部４６で算出された分光反射率をモニタ３６用の色信号に変換し、その色信号に基づく色表示をモニタ３６に行なわせる。変換部４８において、算出された分光反射率をＣＩＥＸＹＺ（１９３１）およびｓＲＧＢに変換する方法については、以下の国際規格に準拠する。
【０１１８】
（ａ）分光反射率からＣＩＥＸＹＺ（１９３１）への変換…ＣＩＥ１５．２ − Colorimetry
（ｂ）ＣＩＥＸＹＺ（１９３１）からｓＲＧＢへの変換…ＩＥＣ６１９６６ Ｐａｒｔ２．０（ｓＲＧＢ）;1999
（ｃ）ＣＩＥＸＹＺ（１９３１）からＣＲＴ出力信号への変換…ＩＥＣ６１９６６ Ｐａｒｔ５（ＣＲＴ）;20XX
（ｄ）ＣＩＥＸＹＺ（１９３１）からＬＣＤ出力信号への変換…ＩＥＣ６１９６６ Ｐａｒｔ６（ＬＣＤ）;20XX
（ｅ）ＣＩＥＸＹＺ（１９３１）からＰＤＰ出力信号への変換…ＩＥＣ６１９６６ Ｐａｒｔ７（ＰＤＰ）;20XX
以上、項目〔１〕〜〔５〕により説明した本実施形態のカラープリンタ出力の色再現方法の流れを、図１３に示すフローチャート（ステップＳ１１〜Ｓ２６）に従って説明する。
【０１１９】
まず、事前に与えられた複数（本実施形態では３５２個）の色票について測定された分光反射率に基づいて、判別部４２が、項目〔１〕で説明した工程に従いモニタ３６で再現すべき色情報がカラープリンタ用のものであるか否かを判別する（ステップＳ１１）。その色情報がカラープリンタ用のものでない場合（ステップＳ１２からＮＯルート）、処理を終了する一方、その色情報がカラープリンタ用のものである場合（ステップＳ１２からＹＥＳルート）、係数算出部４４により、項目〔２〕で説明した工程に従って、係数αi ，βi ，γi ，ξi およびηi を算出する（ステップＳ１３）。
【０１２０】
そして、分光反射率算出部４６は、色の次数に応じて、項目〔３〕で説明した工程に従った分光反射率の算出・推定を行なう（ステップＳ１４〜Ｓ２４）。その際、分光反射率の推定を行なうべき色の次数を判断する（ステップＳ１４）。
１次色の場合（ステップＳ１４から“１次色”ルート）、分光反射率算出部４６は、項目(3-1) で説明した方法に従って、分光反射率を算出・推定する（ステップＳ１５）。
【０１２１】
また、２次色の場合（ステップＳ１４から“２次色”ルート）、分光反射率算出部４６は、分光反射率を特定の値に合わせる必要があるか否かを判断し（ステップＳ１６）、合わせる必要がない場合（ステップＳ１６からＮＯルート）、項目(3-2) で説明した方法に従って分光反射率を算出・推定する（ステップＳ１７）。合わせる必要がある場合（ステップＳ１６からＹＥＳルート）、分光反射率算出部４６は、項目〔４〕で説明した方法に従って分光反射率を算出・推定する（ステップＳ１８）。
【０１２２】
そして、３次色または４次色の場合（ステップＳ１４から“３，４次色”ルート）、分光反射率算出部４６は、４次色については、まず項目(3-3-1) で説明した手順で４次元色立体をブラック一定の３次元色立体にリデュースしてから（ステップＳ１９）、さらに、項目(3-3-2) で説明した手順で、その３次元色立体をブラック一定で且つシアン一定の２次元色平面にリデュース（スライス）する（ステップＳ２０）。３次色については、ステップＳ２１のリダクション処理のみ行なう。
【０１２３】
分光反射率算出部４６は、２次元色平面まで次元数のリダクションを行なったところで、ステップＳ１６〜Ｓ１８と同様、分光反射率を特定の値に合わせる必要があるか否かを判断し（ステップＳ２１）、合わせる必要がない場合（ステップＳ２１からＮＯルート）、項目(3-2) で説明した方法に従って分光反射率を算出・推定する（ステップＳ２２）。合わせる必要がある場合（ステップＳ２１からＹＥＳルート）、分光反射率算出部４６は、項目〔４〕で説明した方法に従って分光反射率を算出・推定する（ステップＳ２３）。
【０１２４】
この後、変換部４８は、例えば項目〔５〕の（ａ）に記載の方法により、分光反射率をＣＩＥＸＹＺ（１９３１）に一旦変換してから（ステップＳ２４）、例えば項目〔５〕の（ｂ）に記載の方法により、ＣＩＥＸＹＺ（１９３１）をｓＲＧＢに変換し（ステップＳ２５）、モニタ３６上で所望のカラープリンタ出力の色表示を行なわせている（ステップＳ２６）。
【０１２５】
ここで、分光反射率の実測値と本実施形態の手法による推定値との精度比較結果を、下記の表２および表３に示す。これらの表２，表３では、シアン（Ｃ），マゼンタ（Ｍ），イエロ（Ｙ），ブラック（Ｋ）のうちで、２色の組み合わせすべて（ＣＫ，ＭＣ，ＭＫ，ＹＣ，ＹＫ，ＹＭ）の分光反射率について、本実施形態の手法によって推定した結果〔ΔＥ94の平均値および最大値と、その最大値の位置（％単位）〕が示されている。
【０１２６】
そして、各辺は１１段階の色信号（色票）をもっているため、一つの組み合わせ、例えば、ＣＫの色票は１１×１１＝１２１個存在し、ＣＫ２次色の任意の分光反射率を、推定すべき分光反射率に対応する色位置を取り囲む４辺上における４０色票の分光反射率を用いて推定した。
なお、表２では、重み係数ｗi を式 (9-1)〜(9-4) で算出した時の推定精度結果が示され、表３では、重み係数ｗi を式(10-1)〜(10-4)で算出した時の推定精度結果が示されている。
【０１２７】
これらの表２および表３を見て明らかな通り、式 (9-1)〜(9-4) を用いた場合も式(10-1)〜(10-4)を用いた場合も、ほぼ同程度の推定精度であり、また、ΔＥ94の平均値および最大値からも極めて良好な精度であることが分かる。
【０１２８】
【表２】
【表３】
このように、本発明の一実施形態の色再現方法によれば、以下のような作用効果を得ることができる。
（Ｉ）ノイゲバウアの式では考慮されていなかった、図２の矢印Ａ３で示すような光の効果を、分光反射率Ｓ50% （λ）に関連した係数βとして考慮することにより、図９にも示した通り仮想１次色についての分光反射率を正確に算出・推定できるので、実際のカラープリンタ出力の色とモニタ３６で再現される色とを一致させることができる。
【０１３１】
（II）ブラックの網点面積率が１００％である際であっても、分光反射率を正確に算出・推定できるとともに、ブラックグラデーション（ブラック１次色）の分光反射率も図１１に示した通り高い精度で算出・推定できるほか、従来法では推定不可能であったドットゲインを算出することもできる。
（III)項目〔４〕で説明した決め打ちと呼ばれる工程を実行することにより、特定の色信号の組み合わせに対して特定の分光反射率を指定できるので、ページ内誤差により不一致となる色を予め測定しておけば、その色については所定の分光反射率を強制的に算出するようにすることができる。従って、ページ内誤差を無視でき、分光反射率の推定精度を高めることが可能になる。
【０１３２】
（IV）判別部４２の機能を用いてカラープリンタ出力と印刷物とを区別することにより、カラープリンタ出力に限定した高精度な色再現が可能になる。
（Ｖ）従来一般的に使用されている９２８色票(ISO12642, JIS X 9202)や７２８色票に比べ、半分以下の色票（本実施形態では３５２色）で、極めて高い精度で色値（分光反射率）の推定を行なうことができる（図９〜図１１や表２，表３参照）。このように推定に必要な色票の数を低減することができるので、色票の数が少ない程、色面積が小さくなり、ページ内誤差も小さくなる。従って、ページ内誤差に起因する推定精度の低下を解消することができる。また、本発明により、事前に指定される複数の色票が、ｎ次元色立体を構成する全ての辺（ｎ・２n-1 本）上に存在していなければ、分光反射率の推定精度が落ちることを理論的に証明することが可能である。
【０１３３】
（VI）本発明の手法を用いることにより、カラープリンタにおいて、試行錯誤に頼ることなくモニタ３６上で確認することで高精度のＵＣＲ(Under Color Removal: 下色除去）量を決定することが可能になるほか、色域圧縮，拡張および色の見えなどの確認を、モニタ３６上で、より容易に行なうことが可能になる。
なお、本発明は上述した実施形態に限定されるものではなく、本発明の趣旨を逸脱しない範囲で種々変形して実施することができる。
【０１３４】
例えば、上述した実施形態では、４色の色材（Ｙ，Ｍ，Ｃ，Ｋ）の組み合わせによるカラープリンタを対象として説明したが、本発明は、この色数に限定されるものではない。
【０１３５】
【発明の効果】
以上詳述したように、本発明のカラープリンタ出力の色再現方法（請求項１〜９）および装置（請求項１０，１１）、並びに、カラープリンタ出力の色再現プログラムを記録したコンピュータ読取可能な記録媒体（請求項１２，１３）によれば、カラープリンタ出力特有の光学現象を考慮することにより、分光反射率を精度良く推定できるので、実際のカラープリンタ出力の色とモニタで再現される色とを一致させることができる。
【図面の簡単な説明】
【図１】本発明の一実施形態としてのカラープリンタ出力の色再現方法を実現するためのコンピュータシステム（色再現装置）の構成を示すブロック図である。
【図２】本実施形態における、分光反射率を算出するための係数の物理的な意味を説明するための図である。
【図３】本実施形態における、係数α，β，γを用いた分光反射率算出手法を説明するための図である。
【図４】本実施形態における、係数ξ，ηを用いた分光反射率算出手法を説明するための図である。
【図５】本実施形態における、仮想２次色の算出・推定に必要なデータの位置関係を説明するための図である。
【図６】本実施形態における、仮想２次色の分光反射率の算出・推定方法を説明するための図である。
【図７】本実施形態における、次元数リダクション処理を説明するための図である。
【図８】本実施形態における、ＹＭＣ３次色の分光反射率の算出・推定方法を説明するための図である。
【図９】マゼンタ１次色の分光反射率を、係数α，β，γを用いて算出・推定した結果を示すグラフである。
【図１０】マゼンタ１次色の分光反射率を、係数ξ，ηを用いて算出・推定した結果を示すグラフである。
【図１１】ブラック１次色の分光反射率を、係数ξ，ηを用いて算出・推定した結果を示すグラフである。
【図１２】ブラック，シアンおよびマゼンタの網点面積率が１００％でイエロの網点面積率が０％〜１００％となるとき（４次色）の分光反射率の測定結果を示すグラフである。
【図１３】本実施形態のカラープリンタ出力の色再現方法の流れを説明するためのフローチャートである。
【符号の説明】
１０ ＣＰＵ
１２ バスライン
１４ ＲＯＭ
１６ ＲＡＭ（記録媒体）
３０ キーボード
３２ マウス
３４ デジタイザ
３６ モニタ
３８ 磁気ディスク
４０，４１ 入出力インタフェース
４２ 判別部
４４ 係数算出部
４６ 分光反射率算出部
４８ 変換部
５０ プリンタ
Claims (8)
- ブラックを含むｎ色（ｎは２以上の自然数）の色材の組み合わせによるカラー印刷を、光透過性を有していない粉体トナーによって行ない且つブラックを最後に刷るカラープリンタの出力を、モニタで再現する方法であって、
ｎ次元色立体を構成するｎ・２n-1本の辺上において事前に指定された複数の色票（各辺において網点面積率の異なる３以上の基準色票を含む）について測定された分光反射率に基づいて、任意の網点面積率に対応する分光反射率を算出するための係数を算出する係数算出ステップと、
該係数算出ステップで算出された係数もしくは前記複数の色票について測定された前記分光反射率に基づいて、前記任意の網点面積率に対応する分光反射率を算出する分光反射率算出ステップと、
該分光反射率算出ステップで算出された分光反射率を、該モニタ用の色信号に変換する変換ステップとを含み、
該係数算出ステップにおいて、前記ｎ・２n-1本の辺のうち、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く辺毎に、前記複数の色票について測定された前記分光反射率に基づいて、前記３以上の基準色票の分光反射率にそれぞれ関連した、前記カラープリンタの出力の色特性を表す３以上の係数を算出するとともに、ブラックグラデーションを表す前記辺上の色票について測定された前記分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、その平均値に基づいて、ドットゲインを表す２つの係数ξ，ηを算出し、
該分光反射率算出ステップにおいて、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く前記辺上における、任意の網点面積率に対応する分光反射率を、該係数算出ステップで算出された前記３以上の係数に基づいて算出する一方、ブラックの網点面積率が１００％である前記辺上における、任意の網点面積率での分光反射率を、前記複数の色票について測定された前記分光反射率に基づいて算出し、ブラックグラデーションを表す前記辺上における任意の網点面積率に対応する分光反射率を、該係数算出ステップで算出された前記２つの係数ξ，ηに基づいて算出することを特徴とする、カラープリンタ出力の色再現方法。 - 該モニタで再現すべき色情報が、該カラープリンタ用のものであるか否かを、前記複数の色票について測定された前記分光反射率に基づいて判別する判別ステップを含み、
該判別ステップにおいて、ブラックの網点面積率が１００％である各辺上の色票について測定された前記分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、前記ブラックの網点面積率が１００％である各辺上の色票について算出された平均値の標準偏差を算出し、その標準偏差と所定値とを比較し、その標準偏差が前記所定値以下であれば該モニタで再現すべき色情報が該カラープリンタ用のものであると判別され、
該判別ステップにおいて該モニタで再現すべき色情報が該カラープリンタ用のものであると判別された場合に、該係数算出ステップ，該分光反射率算出ステップおよび該変換ステップを行なうことを特徴とする、請求項１記載のカラープリンタ出力の色再現方法。 - 該係数算出ステップにおいて、前記複数の色票について測定された前記分光反射率に基づいて、前記３以上の係数を最小二乗法により算出することを特徴とする、請求項１または請求項２に記載のカラープリンタ出力の色再現方法。
- 該分光反射率算出ステップにおいて、前記複数の色票について測定された前記分光反射率に基づいて算出・決定される一定値を、ブラックの網点面積率が１００％である前記辺上における任意の網点面積率での分光反射率として出力することを特徴とする、請求項１〜請求項３のいずれか１項に記載のカラープリンタ出力の色再現方法。
- ブラックを含むｎ色（ｎは２以上の自然数）の色材の組み合わせによるカラー印刷を、光透過性を有していない粉体トナーによって行ない且つブラックを最後に刷るカラープリンタの出力を、モニタで再現する装置であって、
ｎ次元色立体を構成するｎ・２n-1本の辺上において事前に指定された複数の色票（各辺において網点面積率の異なる３以上の基準色票を含む）について測定された分光反射率に基づいて、任意の網点面積率に対応する分光反射率を算出するための係数を算出する係数算出部と、
該係数算出部で算出された係数または前記複数の色票について測定された前記分光反射率に基づいて、前記任意の網点面積率に対応する分光反射率を算出する分光反射率算出部と、
該分光反射率算出部で算出された分光反射率を、該モニタ用の色信号に変換する変換部とをそなえて構成され、
該係数算出部が、前記ｎ・２n-1本の辺のうち、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く辺毎に、前記複数の色票について測定された前記分光反射率に基づいて、前記３以上の基準色票の分光反射率にそれぞれ関連した、前記カラープリンタの出力の色特性を表す３以上の係数を算出するとともに、ブラックグラデーションを表す前記辺上の色票について測定された前記分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、その平均値に基づいて、ドットゲインを表す２つの係数ξ，ηを算出し、
該分光反射率算出部が、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く前記辺上における、任意の網点面積率に対応する分光反射率を、該係数算出部で算出された前記３以上の係数に基づいて算出する一方、ブラックの網点面積率が１００％である前記辺上における、任意の網点面積率での分光反射率を、前記複数の色票について測定された前記分光反射率に基づいて算出し、ブラックグラデーションを表す前記辺上における任意の網点面積率に対応する分光反射率を、該係数算出部で算出された前記２つの係数ξ，ηに基づいて算出することを特徴とする、カラープリンタ出力の色再現装置。 - 該モニタで再現すべき色情報が、該カラープリンタ用のものであるか否かを、前記複数の色票について測定された前記分光反射率に基づいて判別する判別部であって、ブラックの網点面積率が１００％である各辺上の色票について測定された前記分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、前記ブラックの網点面積率が１００％である各辺上の色票について算出された平均値の標準偏差を算出し、その標準偏差と所定値とを比較し、その標準偏差が前記所定値以下であれば該モニタで再現すべき色情報が該カラープリンタ用のものであると判別する判別部をさらにそなえて構成され、
該判別部が、該モニタで再現すべき色情報が該カラープリンタ用のものであると判別した場合に、該係数算出部，該分光反射率算出部および該変換部の動作を実行させることを特徴とする、請求項５記載のカラープリンタ出力の色再現装置。 - ブラックを含むｎ色（ｎは２以上の自然数）の色材の組み合わせによるカラー印刷を、光透過性を有していない粉体トナーによって行ない且つブラックを最後に刷るカラープリンタの出力を、モニタで再現する機能をコンピュータに実行させるための、カラープリンタ出力の色再現プログラムを記録したコンピュータ読取可能な記録媒体であって、
該色再現プログラムが、
ｎ次元色立体を構成するｎ・２n-1本の辺上において事前に指定された複数の色票（各辺において網点面積率の異なる３以上の基準色票を含む）について測定された分光反射率に基づいて、任意の網点面積率に対応する分光反射率を算出するための係数を算出する係数算出部、
該係数算出部で算出された係数または前記複数の色票について測定された前記分光反射率に基づいて、前記任意の網点面積率に対応する分光反射率を算出する分光反射率算出部、および、
該分光反射率算出部で算出された分光反射率を、該モニタ用の色信号に変換する変換部として、該コンピュータを機能させ、
該色再現プログラムが、該コンピュータを該係数算出部として機能させる際に、前記ｎ・２n-1本の辺のうち、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く辺毎に、前記複数の色票について測定された前記分光反射率に基づいて、前記３以上の基準色票の分光反射率にそれぞれ関連した、前記カラープリンタの出力の色特性を表す３以上の係数を算出させるとともに、ブラックグラデーションを表す前記辺上の色票について測定された前記分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、その平均値に基づいて、ドットゲインを表す２つの係数ξ，ηを算出させ、
該色再現プログラムが、該コンピュータを該分光反射率算出部として機能させる際に、ブラックの網点面積率が１００％である辺とブラックグラデーションを表す辺とを除く前記辺上における、任意の網点面積率に対応する分光反射率を、該係数算出部で算出された前記３以上の係数に基づいて算出させる一方、ブラックの網点面積率が１００％である前記辺上における、任意の網点面積率での分光反射率を、前記複数の色票について測定された前記分光反射率に基づいて算出させ、ブラックグラデーションを表す前記辺上における任意の網点面積率に対応する分光反射率を、該係数算出部で算出された前記２つの係数ξ，ηに基づいて算出させることを特徴とする、カラープリンタ出力の色再現プログラムを記録したコンピュータ読取可能な記録媒体。 - 該色再現プログラムが、
該モニタで再現すべき色情報が、該カラープリンタ用のものであるか否かを、前記複数の色票について測定された前記分光反射率に基づいて判別する判別部であって、ブラックの網点面積率が１００％である各辺上の色票について測定された前記分光反射率を、可視波長域すべてに亘りもしくはその可視波長域の一部に亘り平均し、前記ブラックの網点面積率が１００％である各辺上の色票について算出された平均値の標準偏差を算出し、その標準偏差と所定値とを比較し、その標準偏差が前記所定値以下であれば該モニタで再現すべき色情報が該カラープリンタ用のものであると判別する判別部として該コンピュータを機能させ、該モニタで再現すべき色情報が該カラープリンタ用のものであると判別された場合に、該コンピュータを該係数算出部，該分光反射率算出部および該変換部として機能させることを特徴とする、請求項７記載のカラープリンタ出力の色再現プログラムを記録したコンピュータ読取可能な記録媒体。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP30897099A JP4068771B2 (ja) | 1999-10-29 | 1999-10-29 | カラープリンタ出力の色再現方法および装置並びにカラープリンタ出力の色再現プログラムを記録したコンピュータ読取可能な記録媒体 |
| US09/627,455 US6738168B1 (en) | 1999-10-29 | 2000-07-27 | Method and apparatus for reproducing color printer output, and computer-readable recording medium storing color-printer-output reproducing program |
| EP00306620A EP1096787A3 (en) | 1999-10-29 | 2000-08-03 | Method and apparatus for reproducing color printer output, and color-printer-output reproducing program |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP30897099A JP4068771B2 (ja) | 1999-10-29 | 1999-10-29 | カラープリンタ出力の色再現方法および装置並びにカラープリンタ出力の色再現プログラムを記録したコンピュータ読取可能な記録媒体 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2001128020A JP2001128020A (ja) | 2001-05-11 |
| JP4068771B2 true JP4068771B2 (ja) | 2008-03-26 |
Family
ID=17987421
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP30897099A Expired - Fee Related JP4068771B2 (ja) | 1999-10-29 | 1999-10-29 | カラープリンタ出力の色再現方法および装置並びにカラープリンタ出力の色再現プログラムを記録したコンピュータ読取可能な記録媒体 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US6738168B1 (ja) |
| EP (1) | EP1096787A3 (ja) |
| JP (1) | JP4068771B2 (ja) |
Families Citing this family (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3998875B2 (ja) * | 1999-10-21 | 2007-10-31 | 富士フイルム株式会社 | プルーフ作成方法および装置 |
| US20020033959A1 (en) * | 2000-08-02 | 2002-03-21 | Jitsuhiko Ando | Image printing system and image printing apparatus |
| US20030204410A1 (en) * | 2002-04-26 | 2003-10-30 | Clariant International, Ltd. | Method and apparatus for approving color samples |
| JP2003333355A (ja) | 2002-05-10 | 2003-11-21 | Canon Inc | 色評価装置及び方法 |
| US7433102B2 (en) | 2002-05-10 | 2008-10-07 | Canon Kabushiki Kaisha | Reproduction color prediction apparatus and method |
| WO2003095212A1 (en) * | 2002-05-10 | 2003-11-20 | Canon Kabushiki Kaisha | Reproduction color prediction apparatus and method |
| US7564602B2 (en) | 2003-07-28 | 2009-07-21 | Hoya Corporation | Color chart for adjustment of a color printer and a color monitor utilized with an endoscope |
| US20050024657A1 (en) * | 2003-07-28 | 2005-02-03 | Pentax Corporation | Color chart |
| JP4725190B2 (ja) * | 2005-05-25 | 2011-07-13 | 凸版印刷株式会社 | 色情報処理方法および色情報処理装置ならびにそのプログラム |
| JP4586668B2 (ja) * | 2005-08-01 | 2010-11-24 | 凸版印刷株式会社 | 色情報処理方法および色情報処理装置ならびにそのプログラム |
| JP5147390B2 (ja) * | 2007-12-28 | 2013-02-20 | キヤノン株式会社 | 色処理装置およびその方法 |
| EP2381230B1 (de) * | 2010-04-23 | 2018-06-13 | GMG GmbH & Co. KG | Verfahren zur Erstellung von Überdruckvorhersagen |
| JP5713594B2 (ja) * | 2010-07-15 | 2015-05-07 | キヤノン株式会社 | 装置、記録装置、方法およびプログラム |
| JP6952091B2 (ja) * | 2019-10-28 | 2021-10-20 | 株式会社トーエネック | オイル劣化診断装置およびオイル劣化診断方法 |
Family Cites Families (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5018085A (en) * | 1988-06-16 | 1991-05-21 | Hallmark Cards, Inc. | Color printing system usable for reproduction of computer-generated images |
| EP0462712B1 (en) * | 1990-05-31 | 1998-03-04 | Mitsubishi Denki Kabushiki Kaisha | Signal processor for video printer |
| US5081529A (en) * | 1990-12-18 | 1992-01-14 | Eastman Kodak Company | Color and tone scale calibration system for a printer using electronically-generated input images |
| US5172224A (en) * | 1990-12-18 | 1992-12-15 | Eastman Kodak Company | Printer calibration method using electronically-decoupled color and tone scale adjustments |
| US5313291A (en) * | 1991-11-25 | 1994-05-17 | Xerox Corporation | Method for matching color prints to color images on a monitor screen |
| US5309257A (en) * | 1991-12-31 | 1994-05-03 | Eastman Kodak Company | Method and apparatus for providing color matching between color output devices |
| US5459678A (en) * | 1992-07-02 | 1995-10-17 | Feasey; Michael F. | Method and calibration apparatus for calibrating computer monitors used in the printing and textile industries |
| JP3183769B2 (ja) | 1994-02-23 | 2001-07-09 | 大日本スクリーン製造株式会社 | カラー印刷物の再現方法 |
| JPH1073487A (ja) | 1996-06-28 | 1998-03-17 | Dainippon Screen Mfg Co Ltd | カラー印刷物のシミュレーション方法およびシステム、並びに、その処理を実行するためのプログラムを記録した記録媒体 |
| JPH1073488A (ja) | 1996-06-28 | 1998-03-17 | Dainippon Screen Mfg Co Ltd | カラー印刷物のシミュレーション方法およびシステム、並びに、その処理を実行するためのプログラムを記録した記録媒体 |
| JPH1073489A (ja) | 1996-07-02 | 1998-03-17 | Dainippon Screen Mfg Co Ltd | カラー印刷物のシミュレーション方法およびシステム、並びに、その処理を実行するためのプログラムを記録した記録媒体 |
| JPH11132850A (ja) | 1997-08-29 | 1999-05-21 | Dainippon Screen Mfg Co Ltd | 色表現データの変換およびこれを用いたカラー印刷物のシミュレーション |
-
1999
- 1999-10-29 JP JP30897099A patent/JP4068771B2/ja not_active Expired - Fee Related
-
2000
- 2000-07-27 US US09/627,455 patent/US6738168B1/en not_active Expired - Fee Related
- 2000-08-03 EP EP00306620A patent/EP1096787A3/en not_active Withdrawn
Also Published As
| Publication number | Publication date |
|---|---|
| EP1096787A2 (en) | 2001-05-02 |
| JP2001128020A (ja) | 2001-05-11 |
| EP1096787A3 (en) | 2004-10-27 |
| US6738168B1 (en) | 2004-05-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4068771B2 (ja) | カラープリンタ出力の色再現方法および装置並びにカラープリンタ出力の色再現プログラムを記録したコンピュータ読取可能な記録媒体 | |
| US6657746B1 (en) | Look-up table constructing method and color conversion apparatus | |
| US6705703B2 (en) | Determination of control points for construction of first color space-to-second color space look-up table | |
| CN100376105C (zh) | 调整图像基色成分的设备和方法 | |
| US7573607B2 (en) | Method of selecting inks for use in imaging with an imaging apparatus | |
| US20160360072A1 (en) | Color estimation system, separation data generation system, color estimation method and separation data generation method | |
| EP1526711A1 (en) | method of converting color data | |
| CN101453548A (zh) | 颜色处理方法和设备 | |
| US20050253866A1 (en) | Method and apparatus for creating profile | |
| US5425134A (en) | Print color material amount determining method | |
| CN101426072B (zh) | 图像处理装置和图像处理方法 | |
| EP2421240A1 (en) | Image processing apparatus and image processing method | |
| US20070253039A1 (en) | Image processing apparatus | |
| US11455130B2 (en) | Method of color mapping of oligochromatic patterns | |
| JP2015142250A (ja) | 色変換対応情報作成装置、方法、プログラム、印刷装置の製造システム及び製造方法 | |
| JP2003125210A (ja) | 色処理方法、記憶媒体、色処理装置、色変換装置、および画像形成装置 | |
| JP4002006B2 (ja) | 色校正装置、色校正方法、および色校正プログラムを記録した記録媒体 | |
| JP6349708B2 (ja) | ルックアップテーブル生成方法、及び、ルックアップテーブル生成装置 | |
| CN107465847A (zh) | 打印机及颜色转换控制程序以及颜色转换控制方法 | |
| US7023585B1 (en) | Multi-dimensional edge interpolation | |
| US7059701B2 (en) | Method for calibrating production printing cartridges for use in an imaging system | |
| US20030019381A1 (en) | Method of generating color conversion table of high accuracy, color conversion method, and color conversion program | |
| US6160635A (en) | Color conversion method using ordered color signal interpolation | |
| JP4518280B2 (ja) | 色域形成装置、色変換装置、色域形成プログラム、色変換プログラム | |
| JP4586668B2 (ja) | 色情報処理方法および色情報処理装置ならびにそのプログラム |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050830 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050920 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20051121 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20051220 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060213 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20060213 |
|
| A911 | Transfer of reconsideration by examiner before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A911 Effective date: 20060323 |
|
| A912 | Removal of reconsideration by examiner before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A912 Effective date: 20060428 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080111 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110118 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110118 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120118 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130118 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130118 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140118 Year of fee payment: 6 |
|
| LAPS | Cancellation because of no payment of annual fees |