JP2000009537A - カラー印刷物のシミュレーション方法およびそのシステム - Google Patents

Abstract

(57)【要約】 【課題】 ３次元空間に配置されたカラー印刷物を、出
力デバイスの色再現範囲に応じてなるべく忠実に再現す
ることのできる技術を提供する。 【解決手段】 カラー印刷物に所定の光を照射した際に
所定の観察点で観察される反射光の照度スペクトルは、
光源と印刷物と観察点との位置関係に依存する拡散反射
成分と鏡面反射成分との線形結合である所定の関数で表
される。この関数の拡散反射成分と鏡面反射成分とのそ
れぞれに所定の定数あるいは所定の関数を乗じておけ
ば、照度スペクトルの輝度Ｙの値を任意の値に調整する
ことができる。したがって、印刷物を再現する出力デバ
イスの色再現範囲が小さい場合でも、その出力デバイス
の色再現範囲内にほぼ収まるように調整することができ
る。これにより、出力デバイスの色再現範囲に応じてな
るべく忠実に３次元空間に配置されたカラー印刷物を再
現することが可能となる。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】この発明は、カラー印刷物を
ディスプレイデバイスやプリンタ等の出力デバイスで再
現する技術に関する。

【０００２】

【従来の技術】カラー印刷物をディスプレイデバイスや
プリンタなどの種々の出力デバイスによって再現する際
には、実物の印刷物になるべく近い色を再現することが
望ましい。カラー印刷物の色再現方法としては、本願出
願人によって開示された特開平７−２３４１５８号公報
に示されるような一般化された「拡張２色性反射モデ
ル」の式を用いることができる。この式に従って変換さ
れた色データを用いれば、出力デバイスにおいてカラー
印刷物をかなり忠実に再現することができる。

【０００３】

【発明が解決しようとする課題】しかし、上記の一般化
された拡張２色性反射モデルの式を用いる場合には、カ
ラー印刷物を忠実に再現できるか否かは、出力デバイス
の色再現の能力（色再現範囲）に依存するという問題が
あった。すなわち、出力デバイスによっては、変換され
た色データがその出力デバイスの色再現範囲を超えてし
まい、うまくカラー印刷物を再現できない場合があると
いう問題があった。

【０００４】この発明は、従来技術における上述の課題
を解決するためになされたものであり、３次元空間に配
置されたカラー印刷物を、出力デバイスの色再現範囲に
応じてなるべく忠実に再現することのできる技術を提供
することを目的とする。

【０００５】

【課題を解決するための手段およびその作用・効果】上
述の課題の少なくとも一部を解決するため、本発明の第
１の方法は、３次元空間に配置されたカラー印刷物に関
してレンダリング処理を行なうことによって、出力デバ
イスに前記カラー印刷物を再現するシミュレーション方
法であって、（ａ）前記カラー印刷物に所定の輝度スペ
クトルφ（λ）の光を照射した場合に所定の観察点で観
察される反射光の照度スペクトルＩ（λ，θ，ρ）を、
数式１１に従って求める工程と、（ｂ）前記反射光の照
度スペクトルＩ（λ，θ，ρ）から所定のデバイス・イ
ンディペンデント表色系における３刺激値を求める工程
と、（ｃ）前記所定のデバイス・インディペンデント表
色系における３刺激値を前記出力デバイスの所定の表色
系の色データに変換する工程と、を備え、前記工程
（ａ）は、前記数式１１の所定の定数Ｃｂ，Ｃｓの少な
くとも一方を１未満の値に設定することによって、前記
反射光の照度スペクトルＩ（λ，θ，ρ）の輝度の範囲
を調整する工程と、を備えることを特徴とする。

【数１１】 ここで、λは光の波長、Ｓｂ（λ），Ｓｓ（λ）はそれ
ぞれ所定の反射係数、Ｃｂ，Ｃｓは所定の定数、ρは前
記光の反射方向と観察方向のずれ角、ｆ（ρ）は前記ず
れ角ρに応じて単調に減少する所定の関数、θは反射
角、Ｉｅ（λ）は前記観察点で観察される環境光の所定
の照度スペクトルである。

【０００６】上記の第１のシミュレーション方法におい
ては、反射光の照度スペクトルＩ（λ，θ，ρ）が数式
１１によって表され、それぞれの反射係数Ｓｂ（λ），
Ｓｓ（λ）に定数Ｃｂ，Ｃｓが乗じられている。したが
って、定数Ｃｂ，Ｃｓの少なくとも一方を１未満の値に
設定することによって、照射スペクトルＩ（λ，θ，
ρ）の輝度を出力デバイスの色再現範囲に応じた所定の
レベルに調整することができるため、３次元空間に配置
されたカラー印刷物を、出力デバイスの色再現範囲に応
じてなるべく忠実に再現することが可能となる。

【０００７】本発明の第２の方法は、３次元空間に配置
されたカラー印刷物に関してレンダリング処理を行なう
ことによって、出力デバイスに前記カラー印刷物を再現
する方法であって、（ａ）前記カラー印刷物に所定の輝
度スペクトルφ（λ）の光を照射した場合に所定の観察
点で観察される反射光の照度スペクトルＩ（λ，θ，
ρ）を、数式１２に従って求める工程と、（ｃ）前記反
射光の照度スペクトルＩ（λ，θ，ρ）から所定のデバ
イス・インディペンデント表色系における３刺激値を求
める工程と、（ｄ）前記所定のデバイス・インディペン
デント表色系における３刺激値を前記出力デバイスの所
定の表色系の色データに変換する工程と、を備え、前記
工程（ａ）は、前記数式１２の所定の関数ｆｂ（ρ），
ｆｓ（ρ）のそれぞれを１以下の値の範囲をとるよう設
定することによって、前記反射光の照度スペクトルＩ
（λ，θ，ρ）の輝度の範囲を調整する工程と、を備え
ることを特徴とする。

【数１２】 ここで、λは光の波長、Ｓｂ（λ），Ｓｓ（λ）はそれ
ぞれ所定の反射係数、ρは前記光の反射方向と観察方向
のずれ角、ｆｂ（ρ）は前記ずれ角ρに応じて単調に増
加する所定の関数、ｆｓ（ρ）は前記ずれ角ρに応じて
単調に減少する所定の関数、θは反射角、Ｉｅ（λ）は
前記観察点で観察される環境光の所定の照度スペクトル
である。

【０００８】上記の第２のシミュレーション方法におい
ては、反射光の照度スペクトルＩ（λ，θ，ρ）が数式
１２によって表され、それぞれの反射係数Ｓｂ（λ），
Ｓｓ（λ）にずれ角ρに応じて単調に増加する関数ｆｂ
（ρ），単調に減少する関数ｆｓ（ρ）が乗じられてい
る。したがって、２つの関数ｆｂ（ρ），ｆｓ（ρ）の
それぞれを１以下の値となるように設定することによっ
て、照射スペクトルＩ（λ，θ，ρ）の輝度を出力デバ
イスの色再現範囲に応じた所定のレベルに調整すること
ができる。また、関数ｆｂ（ρ），ｆｓ（ρ）を乗じる
ことにより、ρに応じて反射係数Ｓｓ（λ）の照度スペ
クトルＩ（λ，θ，ρ）への寄与が大きいときには、反
射係数Ｓｂ（λ）の照度スペクトルＩ（λ，θ，ρ）へ
の寄与を小さくすることができる。一方、反射係数Ｓｓ
（λ）の照度スペクトルＩ（λ，θ，ρ）への寄与が小
さいときには、反射係数Ｓｓ（λ）の照度スペクトルＩ
（λ，θ，ρ）への寄与を大きくすることができる。し
たがって、主となる反射係数による色味を用いて３次元
空間に配置されたカラー印刷物をなるべく忠実に再現す
ることが可能となる。

【０００９】上記のカラー印刷物のシミュレーション方
法において、前記所定の関数ｆｂ（ρ），ｆｓ（ρ）の
それぞれは数式１３によって表されることが好ましい。

【数１３】 ここで、Ｃｂ，Ｃｓは所定の定数、ｍは１以下の所定の
正の定数、ｆ（ρ）は前記ずれ角ρに応じて単調に減少
する所定の関数である。

【００１０】こうすれば、関数ｆｂ（ρ），ｆｓ（ρ）
は、関数ｆ（ρ）で関連付けられるため、ずれ角ρに応
じて主となる反射係数を用いて３次元空間に配置された
カラー印刷物を再現することが容易となる。

【００１１】上記の第１および第２のシミュレーション
方法において、前記所定の関数ｆ（ρ）は、数式１４ま
たは数式１５であることが好ましい。

【数１４】 ここで、ｎは約３５０から約４００の間の定数である。

【数１５】 ここで、σは約７０から約９０の間の定数である。

【００１２】関数ｆ（ρ）として上記のようなずれ角ρ
に応じて単調に減少する関数を用いれば、反射係数の照
度スペクトルＩ（λ，θ，ρ）への寄与をうまく表すこ
とができる。

【００１３】本発明の第１の装置は、３次元空間に配置
されたカラー印刷物に関してレンダリング処理を行なう
ことによって、出力デバイスに前記カラー印刷物を再現
するシステムであって、前記カラー印刷物に所定の輝度
スペクトルφ（λ）の光を照射した場合に所定の観察点
で観察される反射光の照度スペクトルＩ（λ，θ，ρ）
を、数式１６に従って求める照度スペクトル算出部と、
前記反射光の照度スペクトルＩ（λ，θ，ρ）から所定
のデバイス・インディペンデント表色系における３刺激
値を求める３刺激値算出部と、前記所定のデバイス・イ
ンディペンデント表色系における３刺激値を前記出力デ
バイスの所定の表色系の色データに変換する色データ変
換部と、を備え、前記照度スペクトル算出部は、前記数
式１６の所定の定数Ｃｂ，Ｃｓの少なくとも一方を１未
満の値に設定することによって、前記反射光の照度スペ
クトルＩ（λ，θ，ρ）の輝度の範囲を調整することを
特徴とする。

【数１６】 ここで、λは光の波長、Ｓｂ（λ），Ｓｓ（λ）はそれ
ぞれ所定の反射係数、Ｃｂ，Ｃｓは所定の定数、ρは前
記光の反射方向と観察方向のずれ角、ｆ（ρ）は前記ず
れ角ρに応じて単調に減少する所定の関数、θは反射
角、Ｉｅ（λ）は前記観察点で観察される環境光の所定
の照度スペクトルである。

【００１４】上記の第１のシミュレーションシステムに
おいては、照度スペクトル算出部が反射光の照度スペク
トルＩ（λ，θ，ρ）を数式１６に従って求める。数式
１６の反射係数Ｓｂ（λ），Ｓｓ（λ）には、それぞれ
定数Ｃｂ，Ｃｓが乗じられている。したがって、照度ス
ペクトル算出部が、定数Ｃｂ，Ｃｓの少なくとも一方を
１未満の値に設定することによって、照射スペクトルＩ
（λ，θ，ρ）の輝度を出力デバイスの色再現範囲に応
じた所定のレベルに調整することができる。こうすれ
ば、色データ変換部によって変換された色データを用い
て、３次元空間に配置されたカラー印刷物を出力デバイ
スの色再現範囲に応じてなるべく忠実に再現することが
可能となる。

【００１５】本発明の第２の装置は、３次元空間に配置
されたカラー印刷物に関してレンダリング処理を行なう
ことによって、出力デバイスに前記カラー印刷物を再現
するシステムであって、前記カラー印刷物に所定の輝度
スペクトルφ（λ）の光を照射した場合に所定の観察点
で観察される反射光の照度スペクトルＩ（λ，θ，ρ）
を、数式１７に従って求める照度スペクトル算出部と、
前記反射光の照度スペクトルＩ（λ，θ，ρ）から所定
のデバイス・インディペンデント表色系における３刺激
値を求める３刺激値算出部と、前記所定のデバイス・イ
ンディペンデント表色系における３刺激値を前記出力デ
バイスの所定の表色系の色データに変換する色データ変
換部と、を備え、前記照度スペクトル算出部は、前記数
式１７の所定の関数ｆｂ（ρ），ｆｓ（ρ）のそれぞれ
を１以下の値の範囲をとるよう設定することによって、
前記反射光の照度スペクトルＩ（λ，θ，ρ）の輝度の
範囲を調整することを特徴とする。

【数１７】 ここで、λは光の波長、Ｓｂ（λ），Ｓｓ（λ）はそれ
ぞれ所定の反射係数、ρは前記光の反射方向と観察方向
のずれ角、ｆｂ（ρ）は前記ずれ角ρに応じて単調に増
加する所定の関数、ｆｓ（ρ）は前記ずれ角ρに応じて
単調に減少する所定の関数、θは反射角、Ｉｅ（λ）は
前記観察点で観察される環境光の所定の照度スペクトル
である。

【００１６】上記の第２のシミュレーションシステムに
おいては、照度スペクトル算出部が反射光の照度スペク
トルＩ（λ，θ，ρ）を数式１７に従って求める。数式
１７の反射係数Ｓｂ（λ），Ｓｓ（λ）には、それぞれ
ずれ角ρに応じて単調に増加する関数ｆｂ（ρ），単調
に減少する関数ｆｓ（ρ）が乗じられている。したがっ
て、照度スペクトル算出部が、２つの関数ｆｂ（ρ），
ｆｓ（ρ）のそれぞれを１以下の値となるように設定す
ることによって、照射スペクトルＩ（λ，θ，ρ）の輝
度を出力デバイスの色再現範囲に応じた所定のレベルに
調整することができる。また、関数ｆｂ（ρ），ｆｓ
（ρ）を乗じることにより、ρに応じて反射係数Ｓｓ
（λ）の照度スペクトルＩ（λ，θ，ρ）への寄与が大
きいときには、反射係数Ｓｂ（λ）の照度スペクトルＩ
（λ，θ，ρ）への寄与を小さくすることができる。一
方、反射係数Ｓｓ（λ）の照度スペクトルＩ（λ，θ，
ρ）への寄与が小さいときには、反射係数Ｓｓ（λ）の
照度スペクトルＩ（λ，θ，ρ）への寄与を大きくする
ことができる。したがって、主となる反射係数による色
味を用いて３次元空間に配置されたカラー印刷物をなる
べく忠実に再現することが可能となる。

【００１７】

【発明の他の態様】この発明は、以下のような態様も含
んでいる。第１の態様は、コンピュータに上記の発明の
各工程または各手段の機能を実行させるコンピュータプ
ログラムを記録した記録媒体である。記録媒体として
は、フレキシブルディスクやＣＤ−ＲＯＭなどのコンピ
ュータが読取り可能な携帯型の記憶媒体や、コンピュー
タシステムの内部記憶装置（ＲＡＭやＲＯＭなどのメモ
リ）および外部記憶装置、あるいは、これ以外のコンピ
ュータプログラムが記録された媒体であってコンピュー
タシステムが読取り可能な種々の媒体を利用できる。

【００１８】第２の態様は、コンピュータに上記の発明
の各工程または各手段の機能を実行させるコンピュータ
プログラムを通信経路を介して供給するプログラム供給
装置である。

【００１９】

【発明の実施の形態】Ａ．シミュレーション方法の概
要：次に、本発明の実施の形態を実施例に基づき説明す
る。図１は、印刷紙上に塗布されたインクによって入射
光が反射される様子を示す説明図である。図１に示すよ
うに、インクへの入射光Ｌｉは２つの異なった経路で反
射される。第１の反射光は、インクの表面と空気層との
境界において反射される鏡面反射光Ｌｓである。第２の
反射光は、インクの表面を通過した光がインク中の色素
粒子によって散乱された後に、外部に出射される拡散反
射光Ｌｂである。

【００２０】図２は、光源と印刷物と視覚系（観察者）
との関係を示す説明図である。図２に示すように、視覚
系の位置は、鏡面反射光Ｌｓの方向から角度ρ（以下、
「ずれ角」と呼ぶ）だけずれているのが普通である。印
刷物表面における光の反射点と視覚系とを結ぶ方向（観
察方向）は、入射光Ｌｉと反射光Ｌｓとで構成される平
面上に存在しない場合があるが、この場合にも、鏡面反
射光Ｌｓと観察方向とのなす角が、ずれ角ρとして定義
される。

【００２１】図２の観察者によって観察される反射光の
照度は、次の数式１８に示すように、鏡面反射光Ｌｓの
照度と拡散反射光Ｌｂの照度の線形結合で表わされると
考えられる。

【００２２】

【数１８】

【００２３】ここで、λは光の波長、Ｉ（λ，θ，ρ）
は観察される反射光の照度スペクトル（分光照度）、Ｉ
ｓ（λ，θ，ρ）は鏡面反射光Ｌｓの照度スペクトル、
Ｉｂ（λ，θ，ρ）は拡散反射光Ｌｂの照度スペクトル
である。すなわち、照度スペクトルＩ（λ，θ，ρ），
Ｉｓ（λ，θ，ρ）及びＩｂ（λ，θ，ρ）は反射角
θ、ずれ角ρ及び波長λに依存する。

【００２４】数式１８において、照度スペクトルＩｓ
（λ，θ，ρ），Ｉｂ（λ，θ，ρ）をそれぞれ角度成
分と波長成分に分離できると仮定すると、数式１８は次
の数式１９のように書き換えられる。

【００２５】

【数１９】

【００２６】ここで、φ（λ）は入射光の波長λに依存
する輝度スペクトル（分光輝度）である。また、Ｓｂ
（λ）は波長λの光に対する拡散反射係数、Ｓｓ（λ）
は波長λの光に対する鏡面反射係数である。

【００２７】ところで、鏡面反射光Ｌｓの照度は角度θ
に沿った反射方向（すなわち、ずれ角ρが０°）で観察
した場合に最も強く、この方向からずれるに従って急激
に低下する。このような現象に対応して、数式１９にお
ける鏡面反射光Ｌｓの角度成分Ａｓ（θ，ρ）をずれ角
ρに応じて単調に減少する関数ｆ（ρ）で表すことがで
きる。このような関数ｆ（ρ）としては、例えば、数式
２０や数式２１で表される関数が好ましいことが発明者
による実験によって明らかとなっている。

【００２８】

【数２０】

【００２９】

【数２１】

【００３０】ここで、数式２０のｎは定数であり、約３
５０〜４００の範囲の値で設定することが好ましい。ま
た、数式２１のσは定数であり、約７０〜約９０の範囲
の値で設定することが好ましい。

【００３１】一方、拡散反射光Ｌｂの角度成分Ａｂ
（θ，ρ）は cosθで表わすことができる。すなわち、
拡散反射光Ｌｂは、光の反射角θのみに依存して、ずれ
角ρ（観察方向）には依存しない。これらの仮定の下で
は、数式１９は次の数式２２のように書き換えられる。

【００３２】

【数２２】

【００３３】なお、本明細書において、観察される反射
光の照度スペクトルＩ（λ，θ，ρ）を構成する拡散反
射光Ｌｂの拡散反射成分Ｓｂ（λ）と角度成分cosθと
の積を「拡散反射成分」と呼び、鏡面反射光Ｌｓの鏡面
反射成分Ｓｓ（λ）と角度成分ｆ（ρ）との積を「鏡面
反射成分」と呼ぶ。

【００３４】実際に印刷物を観察する際には、標準的な
光源の他に自然光など、直接及び印刷物で反射してから
観察点で観察される環境光が存在する。そこで、数式２
２の右辺に観察点で観察される環境光の照度スペクトル
Ｉｅ（λ）を加算することによって、反射光の照度スペ
クトルＩ（λ，θ，ρ）を与える一般化された拡張２色
性反射モデルの式である数式２３が得られる。

【００３５】

【数２３】

【００３６】数式２３において、反射係数Ｓｓ（λ），
Ｓｂ（λ）、および関数ｆ（ρ）の定数ｎ（数式２０）
や定数σ（数式２１）の値は、輝度スペクトルφ（λ）
を有する光源を用いた反射光の照度スペクトルに関する
測定を行ない、その結果を解析することによって決定で
きる。なお、この測定については、本願出願人によって
開示された特開平７−２３４１５８号公報に詳述されて
いるので、本明細書での説明は省略する。また、環境光
の照度スペクトルＩｅ（λ）は、色再現の前提条件とし
て予め決定しておくことができる。このとき、数式２３
における未知数は、反射角θとずれ角ρだけになる。し
たがって、観察者と印刷物と光源との３次元的な位置関
係を設定すれば、反射角θとずれ角ρが決まるので、数
式２３を用いることによって反射光の照度スペクトルＩ
（λ，θ，ρ）を求めることができる。

【００３７】反射光の照度スペクトルＩ（λ，θ，ρ）
が得られれば、ＣＩＥ−ＸＹＺ表色系における３刺激値
Ｘ，Ｙ，Ｚは、その定義に従って数式２４で与えられ
る。

【００３８】

【数２４】

【００３９】ここで、ｘ（λ），ｙ（λ），ｚ（λ）は
等色関数である（なお、数式以外では、便宜上、ｘ，
ｙ，ｚの文字上のバーを省略している）。また、ｋは照
明する光の明るさに比例する量である。

【００４０】このように、数式２４で表される３刺激値
Ｘ，Ｙ，Ｚは、観察される反射光の照度スペクトルＩ
（λ，θ，ρ）に、それぞれ等色関数ｘ（λ），ｙ
（λ），ｚ（λ）を乗じて、可視光の波長領域（３８０
ｎｍ〜７８０ｎｍ）で積分して求められる。また、ｋは
入射光の輝度スペクトルφ（λ）に等色関数ｙ（λ）を
乗じて、可視光の波長領域で積分して求められる。した
がって、刺激値Ｙは輝度を表している。ところで、数式
２４に従えば刺激値Ｙの値は規格化され、０〜１の値を
取り得るが、本実施例においてはｋ＝１として刺激値Ｙ
を表す。したがって、刺激値Ｙの値は、入射光の輝度ス
ペクトルφ（λ）の大きさに依存する。

【００４１】数式２４によって求められた３刺激値Ｘ，
Ｙ，Ｚを出力デバイスの表色系（例えばＲＧＢ表色系）
の色データに変換すれば、その出力デバイスにおいて印
刷物の色を忠実に再現することが可能である。

【００４２】Ｂ．シミュレーションシステムの構成と処
理：図３は、この発明の実施例を適用して３次元空間に
配置されたカラー印刷物を再現するためのコンピュータ
システムを示すブロック図である。このコンピュータシ
ステムは、ＣＰＵ１０と、バスライン１２とを備えてい
る。バスライン１２には、ＲＯＭ１４と、ＲＡＭ１６と
が接続されており、また、入出力インタフェイス４０，
４１を介してキーボード３０と、マウス３２と、デジタ
イザ３４と、カラーＣＲＴ３６と、磁気ディスク３８
と、フルカラープリンタ４２とが接続されている。

【００４３】ＲＡＭ１６には、照度スペクトル算出部１
７と、３刺激値算出部１８と、色データ変換部１９との
機能を有するアプリケーションプログラムがロードさ
れ、ＣＰＵ１０がこれらのプログラムを実行することに
よって、以下に説明する処理が実現される。このコンピ
ュータシステムは、いわゆるレンダリング処理を行なう
ことによって３次元空間に配置されたカラー印刷物を再
現する。

【００４４】図４は、レイトレーシング方法によって印
刷物の色を再現する際の、印刷物５０と光源５２とスク
リーン（カラーＣＲＴ３６の画面）５４と観察点５６の
３次元的な配置を示す説明図である。レイトレーシング
法は、３次元コンピュータグラフィックスの分野におい
て周知のレンダリング手法であり、スクリーン５４上の
各画素ＰＸと観察点５６とを通る視線ＲＬを、観察点５
６から光源５２まで逆にたどることによってスクリーン
５４上の各画素ＰＸの色データ（ＲＧＢデータ）を求め
る方法である。なお、印刷物５０は予め複数の微小多角
形ＥＬに分割されている。個々の微小多角形ＥＬは平面
であるが、印刷物５０の全体は３次元的に湾曲している
ようにモデル化することが可能である。なお、入射角θ
とずれ角ρは、図４に示すように、視線ＲＬと印刷物５
０が交わる位置に存在する微小多角形ＥＬにおいて定義
される。

【００４５】図５は、３次元的に湾曲した印刷物を再現
する際のレンダリング処理を説明するための説明図であ
る。図５では、観察点５６と光源５２とが同じ位置Ａに
あり、また、印刷物６０は球体である。図５に示すモデ
ルにおいては、観察点５６で観察される拡散反射光Ｌｂ
のずれ角ρは、印刷物６０によって反射された鏡面反射
光Ｌｓと、印刷物６０に入射する入射光Ｌｉとのなす角
に等しい。すなわち、ずれ角ρは、ρ＝−２θ（θは反
射角）で表される。したがって、入射光Ｌｉの入射角θ
が０°の場合（垂直入射の場合）には、ずれ角ρも０°
となる。

【００４６】このような位置関係が与えられると、図３
のコンピュータシステムはレンダリングイメージを生成
することができる。照度スペクトル算出部１７は、ま
ず、インクの種類や色毎に異なる値を有する２つの反射
係数Ｓｓ（λ），Ｓｂ（λ）を、印刷物６０に塗布され
たインクに応じて選択する。また、与えられた印刷物５
０と光源５２とスクリーン（カラーＣＲＴ３６の画面）
５４と観察点５６の３次元的な配置から、球体表面の各
微小多角形における反射角θとずれ角ρとを求める。照
度スペクトル算出部１７は、反射角θとずれ角ρとの関
係に基づいて球体の印刷物６０の各微小多角形における
反射光の照度スペクトルＩ（λ，θ，ρ）を求める。な
お、種々のインクに応じた反射係数Ｓｓ（λ），Ｓｂ
（λ）と、光源の輝度スペクトルφ（λ）と、環境光の
照度スペクトルＩｅ（λ）とは、予め磁気ディスク３８
に登録されている。また、ずれ角ρに応じて単調に減少
する関数ｆ（ρ）も磁気ディスクに３８に登録されてい
る。

【００４７】３刺激値算出部１８は、照度スペクトル算
出部１７で求められた反射光の照度スペクトルＩ（λ，
θ，ρ）を用いて、前述のように数式２４を用いて計算
することによって、所定のデバイス・インディペンデン
ト表色系における３刺激値を求める。

【００４８】色データ変換部１９は、３刺激値を用い
て、カラーＣＲＴ３６に適した所定の表色系の色データ
（ＲＧＢデータ）に変換する。変換された色データは、
カラーＣＲＴ３６の各画素に割り当てられる。具体的に
は、フレームメモリの画素位置にＲＧＢデータを記憶す
る。なお、ＸＹＺ表色系からＲＧＢ表色系への変換は、
出力デバイスの色表現の特性に依存するので、その出力
デバイスに応じた変換式を用いることによって、その出
力デバイスに応じた色再現を行うことが可能である。こ
のように変換された色データに従って、カラーＣＲＴ３
６は印刷物６０のレンダリングイメージを表示する。

【００４９】上記のようにデバイス・インディペンデン
トなＣＩＥ−ＸＹＺ表色系の３刺激値Ｘ，Ｙ，Ｚを求め
れば、その３刺激値Ｘ，Ｙ，Ｚを出力デバイス（カラー
ＣＲＴ，プリンタ等）に応じた変換式で他の表色系（Ｒ
ＧＢ表色系，ＹＭＣＫ表色系等）に変換することができ
るという利点がある。

【００５０】図６は、図５に示す印刷物６０のレンダリ
ングイメージを模式的に示す説明図である。このレンダ
リングイメージは、カラーＣＲＴ３６に表示される。こ
のレンダリングイメージは、前述した数式２３から得ら
れる反射光の照度スペクトルＩ（λ，θ，ρ）に従って
再現されたものである。なお、このレンダリングイメー
ジには環境光の照度スペクトルＩｅ（λ）の寄与は含ま
れていない。図中に示す５つの円形あるいはリング状の
領域Ａ１〜Ａ５は、それぞれ輝度の異なる領域を示して
おり、中心から外側に向かうほど輝度が低下している。
図５では、各領域Ａ１〜Ａ５内の輝度がそれぞれの領域
内で一定であるかのように描かれているが、実際には各
領域の境目はなく、徐々にその輝度が変化している。ま
た、領域Ｂは、入射光Ｌｉを反射する光がない領域、す
なわち、輝度が０である背景領域を示している。領域Ａ
５と領域Ｂとの境界は、図５の点Ａ（光源５２）からの
入射光Ｌｉが球体の印刷物６０の表面と接する点（すな
わち、θ＝９０°，ρ＝−１８０°）に該当する。

【００５１】ところで、数式２３によって得られる反射
光の照度スペクトルＩ（θ，ρ，λ）をそのまま出力デ
バイスに適した表色系での色データに変換した場合に
は、鏡面反射光Ｌｓの輝度が大きく観察される領域（す
なわち、ずれ角ρが０°近傍の領域）において、出力デ
バイスによってはその色再現範囲外になってしまうこと
がある。例えば、図６のレンダリングイメージの中央部
の領域Ａ１は最も輝度が高く、カラーＣＲＴ３６の色再
現範囲外にある。

【００５２】図７は、数式２３に従って算出された反射
光の照度スペクトルＩ（λ，θ，ρ）から得られた刺激
値Ｘ，Ｙの分布を示すグラフである。図７は、図５に示
した球状の一様な色の印刷物に関する刺激値Ｘ，Ｙの分
布なので、ほぼ一本の曲線ＣＬ１（刺激値曲線と呼ぶ）
の上に乗っている。なお、環境光の照度スペクトルＩｅ
（λ）は０としている。また、ずれ角ρに応じて単調に
減少する関数ｆ（ρ）としては、数式２０に示す関数を
用いており、定数ｎを８０とした場合を示している。グ
ラフのＸ軸は照度スペクトルＩ（λ，θ，ρ）のＸＹＺ
表色系における３刺激値のうちの刺激値Ｘを示してお
り、Ｙ軸は刺激値Ｙを示している。これらの刺激値は、
数式２４に従って求められたものである。なお、ＸＹＺ
表色系における刺激値Ｙは輝度を表している。

【００５３】図７の刺激値曲線ＣＬ１の傾きの変化点Ｐ
は、数式２３の拡散反射成分（Ｓｂ（λ）・cosθ）の
寄与と鏡面反射成分（Ｓｓ（λ）・ｆ（ρ））の寄与の
割合が変化する点を示している。すなわち、輝度Ｙが、
変化点Ｐでの値Ｙ_P以上の領域については鏡面反射成分
の寄与が大きく、輝度Ｙ_P以下の領域については拡散反
射成分の寄与が大きい。したがって、輝度Ｙが値Ｙ_P以
上となる領域は、ずれ角ρが０°近傍の領域に対応し、
輝度Ｙが値Ｙ_P以下となる領域は、ずれ角ρが比較的大
きな値をとる領域に対応する。輝度Ｙが最大となる点は
ρ＝０°の点、すなわち図６の領域Ａ１の中心点に相当
し、輝度Ｙが０となる点はρ≦−１８０°の点、すなわ
ち図６の領域Ｂに含まれる点に相当する。

【００５４】また、図７には、出力デバイスとしてのカ
ラーＣＲＴ３６の色再現範囲の上限が破線で示されてい
る。このカラーＣＲＴ３６の色再現範囲は、輝度Ｙが０
〜８０程度までである。通常、輝度Ｙが出力デバイスの
色再現範囲を超える場合には、色再現範囲を超えた部分
に対応する出力デバイスの画素領域の色データは、その
色データが取り得る値の最大値（８ビットならば２５
５）となる。図７に示すように、反射光の照度スペクト
ルＩ（λ，θ，ρ）を与える式として数式２３を用いる
ときには、輝度ＹがカラーＣＲＴ３６の色再現範囲の上
限（輝度８０）を超える部分が生じる場合がある。この
場合には、色再現範囲を超える部分に対応するカラーＣ
ＲＴ３６の画素領域においては、領域内の色の違いを再
現することができないので、すべて白色に近い色として
表現される。図６のレンダリングイメージでは、領域Ａ
がカラーＣＲＴ３６の色再現範囲外となる領域に対応
し、本来の色味が表現されずに白く表示されてしまう。

【００５５】Ｃ．第１のシミュレーション方法：本実施
例の第１のシミュレーション方法においては、上記の出
力デバイスの色再現範囲に起因する問題を解決するため
に照射光の照度スペクトルＩ（λ，θ，ρ）を与える式
として、数式２３に代えて数式２５を用いる。

【００５６】

【数２５】

【００５７】ここで、Ｃｂ，Ｃｓは所定の定数である。
また、以下の説明では、数式２５の環境光の照度スペク
トルＩｅ（λ）が０の場合について説明する。

【００５８】数式２５では、数式２３の拡散反射成分
（Ｓｂ（λ）・cosθ）と鏡面反射成分（Ｓｓ（λ）・
ｆ（ρ））とにそれぞれ定数Ｃｂ，Ｃｓが乗じられてい
る。２つの定数Ｃｂ，Ｃｓの少なくとも一方を１未満の
値に設定すれば、拡散反射成分と鏡面反射成分との輝度
Ｙへの寄与を個別に調整することができるため、輝度Ｙ
をカラーＣＲＴ３６の色再現範囲に収まるように調整す
ることが可能となる。

【００５９】図８は、数式２５に従って、輝度Ｙを出力
デバイスの色再現範囲に応じて調整した場合の刺激値
Ｘ，Ｙの分布を示すグラフである。この刺激値曲線ＣＬ
２は、数式２５の２つの定数Ｃｂ，Ｃｓを共に０．７に
設定した場合の刺激値を示している。また、関数ｆ
（ρ）としては、数式２０に示す関数が用いられてお
り、定数ｎは８０に設定されている。

【００６０】図８では、輝度Ｙの最大値がほぼ８０程度
に調整されており、カラーＣＲＴ３６の色再現範囲内に
ほぼ抑制されている。この刺激値曲線ＣＬ２は、定数Ｃ
ｂ，Ｃｓを共に０．７に設定した場合であるため、図７
の刺激値曲線ＣＬ１を単純に０．７倍したものとなって
いるが、上記のように定数Ｃｂ，Ｃｓは個別に調整可能
である。例えば、鏡面反射成分の輝度Ｙへの寄与のみを
抑制する場合には、Ｃｂを１に設定し、Ｃｓを１未満の
値に設定すればよい。さらに、拡散反射成分の輝度Ｙへ
の寄与も抑制する場合には、Ｃｂも１未満の値に設定す
ればよい。このように、２つの定数Ｃｂ，Ｃｓは任意の
値に設定可能であるが、ＣｂとＣｓとの比Ｃｂ／Ｃｓは
約０．７ないし約１．４の範囲で設定することが好まし
い。この範囲で設定すれば、輝度Ｙを出力デバイスの色
再現範囲内にうまく抑制できるとともに、作成されたレ
ンダリングイメージにおいて印刷物の本来の色味を失わ
ずに表現できる可能性がある。

【００６１】なお、上記の実施例では、関数ｆ（ρ）と
して数式２０を用いているが、これに代えて数式２１を
用いても同様の効果を得ることができる。

【００６２】上記のように、照射光の照度スペクトルＩ
（λ，θ，ρ）として数式２５を用いれば、２つの定数
Ｃｂ，Ｃｓをうまく設定することによって、拡散反射成
分と鏡面反射成分との輝度Ｙへの寄与をそれぞれ調整す
ることができるので、出力デバイスの色再現範囲に応じ
たカラー印刷物の再現が可能となる。

【００６３】Ｄ．第２のシミュレーション方法：本実施
例の第２のシミュレーション方法においては、照射光の
照度スペクトルＩ（λ，θ，ρ）を与える式として、数
式２３に代えて数式２６を用いる。

【００６４】

【数２６】

【００６５】ここで、ｆｂ（ρ）はずれ角ρに応じて単
調に増加する関数であり、ｆｓ（ρ）はずれ角ρに応じ
て単調に減少する関数である。関数ｆｂ（ρ），ｆｓ
（ρ）は、それぞれ１以下の値となるよう調整される。
以下の説明では、数式２６の環境光の照度スペクトルＩ
ｅ（λ）が０の場合について説明する。

【００６６】数式２６では、数式２３の拡散反射成分
（Ｓｂ（λ）・cosθ）にずれ角ρに応じて単調に増加
する関数ｆｂ（ρ）が乗じられており、また、鏡面反射
係数Ｓｓ（λ）にずれ角ρに応じて単調に減少する関数
ｆｓ（ρ）が乗じられている。２つの関数ｆｂ（ρ），
ｆｓ（ρ）はいずれも０〜１の範囲の値をとる。数式２
３を用いた場合にはずれ角ρに応じて鏡面反射成分のみ
の輝度Ｙへの寄与が変化するのに対し、数式２６を用い
れば、ずれ角ρに応じて拡散反射成分の輝度Ｙへの寄与
をも変化させることが可能となる。

【００６７】２つの関数ｆｂ（ρ），ｆｓ（ρ）として
は、例えば数式２７を用いることができる。

【００６８】

【数２７】

【００６９】ここで、Ｃｂ，Ｃｓは１以下の所定の定
数、ｍは１以下の正の定数、ｆ（ρ）はずれ角ρに応じ
て単調に減少する所定の関数である。なお、定数ｍおよ
び定数Ｃｂ，Ｃｓは個別に調整可能である。

【００７０】図９は、数式２６，数式２７に従って、輝
度Ｙを出力デバイスの色再現範囲に応じて調整した場合
の刺激値Ｘ，Ｙの分布を示すグラフである。図９では、
定数Ｃｂ，Ｃｓは共に０．７に設定されており、また定
数ｍは０．８に、定数ｎは３８０に設定されている。

【００７１】なお、図９のグラフでは、関数ｆ（ρ）と
して数式２０が用いられているが、これに代えて数式２
１を用いても同様の効果を得ることができる。また、数
式２７では、２つの関数ｆｂ（ρ），ｆｓ（ρ）は、関
数ｆ（ρ）によって関連付けられているが、関数ｆｂ
（ρ），ｆｓ（ρ）ともに独立の関数を適用できる。図
９では、輝度Ｙの最大値が８０未満に調整されており、
カラーＣＲＴ３６の色再現範囲内に抑制されていること
が分かる。

【００７２】反射光の照度スペクトルＩ（λ，θ，ρ）
として数式２６に示すような関数を用いれば、前述のよ
うに、ずれ角ρに応じて２つの反射成分を変化させるこ
とができる。すなわち、ずれ角ρが０°近傍の鏡面反射
成分の輝度Ｙへの寄与が大きくなる領域においては、拡
散反射成分の輝度Ｙへの寄与を抑制することができる。
このとき数式２３を用いた場合と比べて、鏡面反射成分
の輝度Ｙへの寄与が、拡散反射成分の輝度Ｙへの寄与に
対して相対的に大きくなる。また、ずれ角ρが大きく鏡
面反射成分の輝度Ｙへの寄与が小さくなる領域において
は、拡散反射成分の輝度Ｙへの寄与を増大させることが
できる。このとき、数式２３を用いた場合と比べて、拡
散反射成分の輝度Ｙへの寄与が、鏡面反射成分の輝度Ｙ
への寄与に対して相対的に大きくなる。換言すれば、ず
れ角ρに応じて、鏡面反射成分の輝度Ｙへの寄与が大き
い場合には、鏡面反射成分を反射光の照度スペクトルＩ
（λ，θ，ρ）の主成分とすることができる。一方、鏡
面反射成分の輝度Ｙへの寄与が小さい場合には、拡散反
射成分を反射光の照度スペクトルＩ（λ，θ，ρ）の主
成分とすることができる。

【００７３】このように決定された反射光の照度スペク
トルを用いてシミュレーションを行えば、照射スペクト
ルＩ（λ，θ，ρ）への寄与が大きくなるべき反射成分
を主成分としたレンダリングイメージを作成することが
できる。このとき、鏡面反射成分の大きい領域（ずれ角
ρが０°近傍の領域）では、主として鏡面反射成分の色
味で色を再現でき、拡散反射成分の大きい領域（ずれ角
ρが大きい領域）では、主として拡散反射成分の色味で
色を再現することができる。

【００７４】上記のように、照射光の照度スペクトルＩ
（λ，θ，ρ）として数式２６を用いれば、２つの関数
ｆｂ（ρ），ｆｓ（ρ）をうまく設定することによっ
て、拡散反射成分と鏡面反射成分との輝度Ｙへの寄与を
ずれ角ρに応じてそれぞれ調整することができる。これ
により、輝度Ｙを出力デバイスの色再現範囲内になるべ
く収めるようにすることができるとともに、２つの反射
成分のうちの主となる反射成分の色味を用いて色を再現
することが可能となる。

【００７５】図１０は、図７〜図９の３つの刺激値曲線
ＣＬ１〜ＣＬ３をまとめて示したグラフである。図１０
から分かるように、本発明による数式２５，数式２６を
用いた場合の刺激値曲線ＣＬ２，ＣＬ３の輝度Ｙは、カ
ラーＣＲＴ３６の色再現範囲でうまく調整されている。
上記において説明したように、本発明においては、数式
２３の拡散反射成分と鏡面反射成分とに所定の定数、あ
るいは所定の関数を乗じることにより輝度Ｙの値を出力
デバイスの色再現範囲内にほぼ収まるように調整するこ
とが可能となっている。さらに、これらの所定の定数
や、所定の関数を適当に設定すれば、印刷物で反射され
る反射光の色味をうまく再現することができる。

【００７６】なお、この発明は上記の実施例や実施形態
に限られるものではなく、その要旨を逸脱しない範囲に
おいて種々の態様において実施することが可能であり、
例えば以下のような変形も可能である。

【００７７】（１）上記実施例では、ずれ角ρに応じて
単調に減少する関数ｆ（ρ）として数式２０や数式２１
で表される関数を用いているが、このような関数であれ
ば他の関数を用いてもよい。

【００７８】（２）上記実施例では、デバイス・インデ
ィペンデント表色系における３刺激値として、数式２４
に従いＣＩＥ−ＸＹＺ表色系における３刺激値ＸＹＺを
求めているが他の表色系を用いてもよい。他の表色系と
しては、ＣＩＥ−Ｌ*ａ*ｂ*表色系やＣＩＥ−Ｌ*ｕ*ｖ*
表色系などの任意のデバイス・インディペンデント表色
系を用いることができる。

【図面の簡単な説明】

【図１】印刷紙上に塗布されたインクによって入射光が
反射される様子を示す説明図。

【図２】光源と印刷物と視覚系（観察者）との関係を示
す説明図。

【図３】この発明の実施例を適用して３次元空間に配置
されたカラー印刷物を再現するためのコンピュータシス
テムを示すブロック図。

【図４】レイトレーシング方法によって印刷物の色を再
現する際の、印刷物５０と光源５２とスクリーン（カラ
ーＣＲＴ３６の画面）５４と観察点５６の３次元的な配
置を示す説明図。

【図５】３次元的に湾曲した印刷物を再現する際のレン
ダリング処理を説明するための説明図。

【図６】図５に示す印刷物６０のレンダリングイメージ
を模式的に示す説明図。

【図７】数式２３に従って算出された反射光の照度スペ
クトルＩ（λ，θ，ρ）から得られた刺激値Ｘ，Ｙの分
布を示すグラフ。

【図８】数式２５に従って、輝度Ｙを出力デバイスの色
再現範囲に応じて調整した場合の刺激値Ｘ，Ｙの分布を
示すグラフ。

【図９】数式２６，数式２７に従って、輝度Ｙを出力デ
バイスの色再現範囲に応じて調整した場合の刺激値Ｘ，
Ｙの分布を示すグラフ。

【図１０】図７〜図９の３つの刺激値曲線ＣＬ１〜ＣＬ
３をまとめて示したグラフ。

【符号の説明】

１０…ＣＰＵ １２…バスライン １４…ＲＯＭ １６…ＲＡＭ １７…照度スペクトル算出部 １８…３刺激値算出部 １９…色データ変換部 ３０…キーボード ３２…マウス ３４…デジタイザ ３６…カラーＣＲＴ ３８…磁気ディスク ４０，４１…入出力インタフェイス ４２…フルカラープリンタ ５０…印刷物 ５２…光源 ５４…スクリーン ５６…観察点 ６０…印刷物 Ｌｂ…拡散反射光 Ｌｉ…入射光 Ｌｓ…鏡面反射光 θ…入射角，反射角 ρ…ずれ角

フロントページの続き Ｆターム(参考） 2C262 AA29 AB11 BA01 BA16 2G020 AA08 DA05 DA14 DA15 DA23 DA33 DA34 DA43 5B050 AA09 BA13 EA09 EA30 FA02 FA03 FA05 5C077 LL12 MP08 PP31 PP32 PP37 PP43 PQ12 SS07 TT02 TT06 5C079 HA19 HB01 HB05 HB12 JA25 LA02 LB02 MA11 MA17 MA20 NA03 PA02 PA03

Claims (8)

【特許請求の範囲】
1. 【請求項１】 ３次元空間に配置されたカラー印刷物に
   関してレンダリング処理を行なうことによって、出力デ
   バイスに前記カラー印刷物を再現する方法であって、
   （ａ）前記カラー印刷物に所定の輝度スペクトルφ
   （λ）の光を照射した場合に所定の観察点で観察される
   反射光の照度スペクトルＩ（λ，θ，ρ）を、数式１に
   従って求める工程と、（ｂ）前記反射光の照度スペクト
   ルＩ（λ，θ，ρ）から所定のデバイス・インディペン
   デント表色系における３刺激値を求める工程と、（ｃ）
   前記所定のデバイス・インディペンデント表色系におけ
   る３刺激値を前記出力デバイスの所定の表色系の色デー
   タに変換する工程と、を備え、 前記工程（ａ）は、 前記数式１の所定の定数Ｃｂ，Ｃｓの少なくとも一方を
   １未満の値に設定することによって、前記反射光の照度
   スペクトルＩ（λ，θ，ρ）の輝度の範囲を調整する工
   程と、を備えることを特徴とするカラー印刷物のシミュ
   レーション方法。 【数１】 ここで、λは光の波長、Ｓｂ（λ），Ｓｓ（λ）はそれ
   ぞれ所定の反射係数、Ｃｂ，Ｃｓは所定の定数、ρは前
   記光の反射方向と観察方向のずれ角、ｆ（ρ）は前記ず
   れ角ρに応じて単調に減少する所定の関数、θは反射
   角、Ｉｅ（λ）は前記観察点で観察される環境光の所定
   の照度スペクトルである。
2. 【請求項２】 ３次元空間に配置されたカラー印刷物に
   関してレンダリング処理を行なうことによって、出力デ
   バイスに前記カラー印刷物を再現する方法であって、
   （ａ）前記カラー印刷物に所定の輝度スペクトルφ
   （λ）の光を照射した場合に所定の観察点で観察される
   反射光の照度スペクトルＩ（λ，θ，ρ）を、数式２に
   従って求める工程と、（ｃ）前記反射光の照度スペクト
   ルＩ（λ，θ，ρ）から所定のデバイス・インディペン
   デント表色系における３刺激値を求める工程と、（ｄ）
   前記所定のデバイス・インディペンデント表色系におけ
   る３刺激値を前記出力デバイスの所定の表色系の色デー
   タに変換する工程と、を備え、 前記工程（ａ）は、 前記数式２の所定の関数ｆｂ（ρ），ｆｓ（ρ）のそれ
   ぞれを１以下の値の範囲をとるよう設定することによっ
   て、前記反射光の照度スペクトルＩ（λ，θ，ρ）の輝
   度の範囲を調整する工程と、を備えることを特徴とする
   カラー印刷物のシミュレーション方法。 【数２】 ここで、λは光の波長、Ｓｂ（λ），Ｓｓ（λ）はそれ
   ぞれ所定の反射係数、ρは前記光の反射方向と観察方向
   のずれ角、ｆｂ（ρ）は前記ずれ角ρに応じて単調に増
   加する所定の関数、ｆｓ（ρ）は前記ずれ角ρに応じて
   単調に減少する所定の関数、θは反射角、Ｉｅ（λ）は
   前記観察点で観察される環境光の所定の照度スペクトル
   である。
3. 【請求項３】 請求項２記載のカラー印刷物のシミュレ
   ーション方法であって、 前記所定の関数ｆｂ（ρ），ｆｓ（ρ）のそれぞれは数
   式３によって表される、カラー印刷物のシミュレーショ
   ン方法。 【数３】 ここで、Ｃｂ，Ｃｓは所定の定数、ｍは１以下の所定の
   正の定数、ｆ（ρ）は前記ずれ角ρに応じて単調に減少
   する所定の関数である。
4. 【請求項４】 請求項１または３記載のカラー印刷物の
   シミュレーション方法であって、 前記所定の関数ｆ（ρ）は、数式４または数式５であ
   る、カラー印刷物のシミュレーション方法。 【数４】 ここで、ｎは約３５０から約４００の間の定数である。 【数５】 ここで、σは約７０から約９０の間の定数である。
5. 【請求項５】 ３次元空間に配置されたカラー印刷物に
   関してレンダリング処理を行なうことによって、出力デ
   バイスに前記カラー印刷物を再現するシステムであっ
   て、 前記カラー印刷物に所定の輝度スペクトルφ（λ）の光
   を照射した場合に所定の観察点で観察される反射光の照
   度スペクトルＩ（λ，θ，ρ）を、数式６に従って求め
   る照度スペクトル算出部と、 前記反射光の照度スペクトルＩ（λ，θ，ρ）から所定
   のデバイス・インディペンデント表色系における３刺激
   値を求める３刺激値算出部と、 前記所定のデバイス・インディペンデント表色系におけ
   る３刺激値を前記出力デバイスの所定の表色系の色デー
   タに変換する色データ変換部と、を備え、 前記照度スペクトル算出部は、前記数式６の所定の定数
   Ｃｂ，Ｃｓの少なくとも一方を１未満の値に設定するこ
   とによって、前記反射光の照度スペクトルＩ（λ，θ，
   ρ）の輝度の範囲を調整することを特徴とするカラー印
   刷物のシミュレーションシステム。 【数６】 ここで、λは光の波長、Ｓｂ（λ），Ｓｓ（λ）はそれ
   ぞれ所定の反射係数、Ｃｂ，Ｃｓは所定の定数、ρは前
   記光の反射方向と観察方向のずれ角、ｆ（ρ）は前記ず
   れ角ρに応じて単調に減少する所定の関数、θは反射
   角、Ｉｅ（λ）は前記観察点で観察される環境光の所定
   の照度スペクトルである。
6. 【請求項６】 ３次元空間に配置されたカラー印刷物に
   関してレンダリング処理を行なうことによって、出力デ
   バイスに前記カラー印刷物を再現するシステムであっ
   て、 前記カラー印刷物に所定の輝度スペクトルφ（λ）の光
   を照射した場合に所定の観察点で観察される反射光の照
   度スペクトルＩ（λ，θ，ρ）を、数式７に従って求め
   る照度スペクトル算出部と、 前記反射光の照度スペクトルＩ（λ，θ，ρ）から所定
   のデバイス・インディペンデント表色系における３刺激
   値を求める３刺激値算出部と、 前記所定のデバイス・インディペンデント表色系におけ
   る３刺激値を前記出力デバイスの所定の表色系の色デー
   タに変換する色データ変換部と、を備え、 前記照度スペクトル算出部は、前記数式７の所定の関数
   ｆｂ（ρ），ｆｓ（ρ）のそれぞれを１以下の値の範囲
   をとるよう設定することによって、前記反射光の照度ス
   ペクトルＩ（λ，θ，ρ）の輝度の範囲を調整すること
   を特徴とするカラー印刷物のシミュレーションシステ
   ム。 【数７】 ここで、λは光の波長、Ｓｂ（λ），Ｓｓ（λ）はそれ
   ぞれ所定の反射係数、ρは前記光の反射方向と観察方向
   のずれ角、ｆｂ（ρ）は前記ずれ角ρに応じて単調に増
   加する所定の関数、ｆｓ（ρ）は前記ずれ角ρに応じて
   単調に減少する所定の関数、θは反射角、Ｉｅ（λ）は
   前記観察点で観察される環境光の所定の照度スペクトル
   である。
7. 【請求項７】 請求項６記載のカラー印刷物のシミュレ
   ーションシステムであって、 前記所定の関数ｆｂ（ρ），ｆｓ（ρ）のそれぞれは数
   式８によって表される、カラー印刷物のシミュレーショ
   ンシステム。 【数８】 ここで、Ｃｂ，Ｃｓは所定の定数、ｍは１以下の所定の
   正の定数、ｆ（ρ）は前記ずれ角ρに応じて単調に減少
   する所定の関数である。
8. 【請求項８】 請求項５または７記載のカラー印刷物の
   シミュレーションシステムであって、 前記所定の関数ｆ（ρ）は、数式９または数式１０であ
   る、カラー印刷物のシミュレーションシステム。 【数９】 ここで、ｎは約３５０から約４００の間の定数である。 【数１０】 ここで、σは約７０から約９０の間の定数である。