JPH06109544A - 合成着色料の色を電子ディスプレー上に再現する方法及び装置 - Google Patents

合成着色料の色を電子ディスプレー上に再現する方法及び装置

Info

Publication number
JPH06109544A
JPH06109544A JP5093685A JP9368593A JPH06109544A JP H06109544 A JPH06109544 A JP H06109544A JP 5093685 A JP5093685 A JP 5093685A JP 9368593 A JP9368593 A JP 9368593A JP H06109544 A JPH06109544 A JP H06109544A
Authority
JP
Japan
Prior art keywords
values
rgb
xyz
value
color
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.)
Ceased
Application number
JP5093685A
Other languages
English (en)
Inventor
Jr Louis W Adams
ルイス・ウイリアム・アダムス・ジュニア
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Milliken Research Corp
Original Assignee
Milliken Research Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Milliken Research Corp filed Critical Milliken Research Corp
Publication of JPH06109544A publication Critical patent/JPH06109544A/ja
Ceased legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/46—Measurement of colour; Colour measuring devices, e.g. colorimeters
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/46—Measurement of colour; Colour measuring devices, e.g. colorimeters
    • G01J3/462—Computing operations in or between colour spaces; Colour management systems
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/46—Measurement of colour; Colour measuring devices, e.g. colorimeters
    • G01J3/463—Colour matching
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/46—Measurement of colour; Colour measuring devices, e.g. colorimeters
    • G01J3/465—Measurement of colour; Colour measuring devices, e.g. colorimeters taking into account the colour perception of the eye; using tristimulus detection
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/46—Measurement of colour; Colour measuring devices, e.g. colorimeters
    • G01J2003/466—Coded colour; Recognition of predetermined colour; Determining proximity to predetermined colour

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Spectrometry And Color Measurement (AREA)
  • Controls And Circuits For Display Device (AREA)
  • Image Processing (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)
  • Facsimile Image Signal Circuits (AREA)
  • Processing Or Creating Images (AREA)
  • Image Analysis (AREA)
  • Color Image Communication Systems (AREA)
  • Coloring (AREA)

Abstract

(57)【要約】 【目的】本発明は、陰極線管、液晶ディスプレー、又は
その他のＲＧＢ値を用いるデバイスのごとき電子ディス
プレー上に合成着色料の色を再現する装置及び方法を提
供する。 【構成】基質上又は中の合成着色料の予測値を、着色料
無し、１着色料、及び１対の着色料で調製された試料の
ＸＹＺ測定値から作ることが出来る。その計算方法は光
吸収、光散乱及び光吸収合成の各係数を使用する。イメ
ージディジタイザーは試料からＸＹＺ値を得るために使
用することが出来る。更に、イメージディジタイザーの
ＲＧＢ値が、簡単な方法を用いた非直線モデルで、ＸＹ
Ｚ値に変換される。更に、イメージディジタイザーＲＧ
Ｂ値からＸＹＺ値を作り出す上記の方法が電子ディスプ
レーのためにＸＹＺ値からＲＧＢ値を作るのに用いられ
る。

Description

【発明の詳細な説明】
【０００１】
【産業上の利用分野】本発明は合成着色料の色を電子デ
ィスプレー上に再現する方法及び装置に関する。
【０００２】
【従来の技術】色に関する式を最も正確に計算する方法
はその使用方法が難しかったり、コンピューター的に実
用性が無かったりする。混合体の色を予想する最も簡単
で成功した数学理論はクーベルカ・ムンク・モデル（Ku
belka-Munk）である。適用範囲を広くするために、この
モデルは過度に簡単に割切り過ぎている。クーベルカ・
ムンク・モデルは、着色料を含む完全に平らな媒体に光
が正確に垂直に当たると仮定している。着色料は完全に
基質媒体に混合されていなければならず、出来上がった
着色基質は等方性でなければならない。媒体及び着色料
の屈折率は空気と同じと仮定されており、従って、内部
及び外部における鏡面反射及び屈折は無視される。
【０００３】これらの条件が満足され、且つ、基質が光
学的に厚いと仮定する。この場合は、クーベルカ・ムン
ク理論は次の簡単な関係式を予想している、即ち、Ｋ／
Ｓ＝（１−Ｒ2 ）／２Ｒ。Ｋ及びＳは着色媒体の物理的
性質である。Ｒは測定した色である。この式で表されて
いる関係は、可視スーペクトル・バンドの中の各波長の
光について成立する。Ｒは試料から反射される光の一部
を示す。Ｋ及びＳは、それぞれ、着色料混合体の光の吸
収率及び光の散乱率である。Ｋ／Ｓの方がＲよりも扱い
やすい。これは、混合体の物理的性質（Ｋ及びＳ）が、
かなりの近似をもって、混合体の各着色料の物理性質、
即ち着色料ｉの対応する係数Ｋi 及びＳi 、に比例する
からである。その比例定数が成分濃度Ｃi である。従っ
て、着色料の構成計算を容易にするために、Ｎ種類の着
色料のために、Ｋ＝Ｋ1 Ｃ1 ＋Ｋ2 Ｃ2 ＋…ＫN ＣN ；
Ｓ＝Ｓ1 Ｃ1 ＋Ｓ2 Ｃ2 ＋…ＳN ＣN 、と仮定する。前
と同様に、これらの式は各波長の光に適用される。Ｋ／
ＳとＲとを結ぶ前の式Ｋ／Ｓ＝（１−Ｒ2 ）／２Ｒを変
換し、Ｋ及びＳをＫi 及びＳi に結ぶ上記式を用いるこ
とにより、着色料濃度Ｃi と測定した色Ｒとの間の関係
が得られる。着色料の吸収率及び散乱率の値が、主とし
て、サンプルの色の測定値を含む最小二乗計算法から引
き出される。
【０００４】高い精度を必要とする場合は、この単純な
クーベルカ・ムンク理論は変形されねばならない。基質
表面の反射、内部屈折、及び着色料間の相互作用に対す
る修正が必要である。時には、スペクトル測定を紫外線
の領域にまで広げ、着色料の燐光性を処理しなければな
らない。目標（例えば織物）によっては、その組織が光
沢を持ち、測定によってこれを差し引いたり、数学的に
モデル化することによって相殺したりすることが容易で
ない。これは、結果として得られる着色料の式による予
想値が正確である以前に、Ｋi 及びＳi の演算値を慎重
に読取り、恐らくは更に変形する必要があることを意味
する。
【０００５】コンピュータを利用した設計（ＣＡＤ）に
おいては、色をフォーミュレーションする間、ピジュア
ル・フィードバック（画像によるフィードバック）が好
ましい。これを行う１つの方法は電子ディスプレー上で
演算結果をシミュレートすることである。クーベルカ・
ムンクの計算を行うことには、上述した修正と共に、膨
大な演算が含まれる。多くの波長におけるスペクトル・
データがコンピュータに記憶されねばならない。色の測
定は従来から分光測光器で行われている。これらのデバ
イスは比較的高価で、走査するのに特別な技術専門家を
必要とする。スペクトル・データを電子ディスプレーに
適した色座標に変換する一般的方法には複雑な非線形式
が入ってくる。ＣＡＤは、色の予測条件が、例えば織物
の染料の色座標ほど厳格でない場合の１適用例である。
本発明は、余り厳格でない適用に対して、従来技術に開
示されていない形で、これらの問題を解決する。
【０００６】
【課題を解決するための手段】本出願は、電子ディスプ
レー上に試料の合成彩色を再現する方法及びその装置を
開示している。この電子ディスプレーは陰極線管、液晶
ディスプレー、又は赤、緑、青（ＲＧＢ）の色座標を用
いたその他の型の電子デイスプレーのいずれでも良い。
本出願では、一般的に、着色料は合成されていて、基質
上に単に隣り合わせに並べられているものではないとす
る。この着色料が殆ど透明であり、余り不透明でない場
合は、着色料は層状であってもよい。カラー・イメージ
・デジタイザーは色々なＣＡＤの間に一般的に使用され
る。我々は、従来の色の測定装置よりも安価なカラー・
イメージ・デジタイザーが色の測定値を得るためにどの
ように使われるかを示す。本発明はシュミレーション作
業のためのみのものであり、微妙な着色料調合作業に使
うことは出来ない。
【０００７】我々は色の分析のためにＣＩＥの色のＸＹ
Ｚ三刺激座標を選択する。多くの異なった波長について
測定を行う代わりに、色の三刺激測定値Ｘ、Ｙ、及びＺ
は、それぞれ、赤、緑及び青のスペクトル・バンドの平
均値である。これは色を定量するために必要とする最低
のスペクトル情報である。なんとなれば、色を見ること
によって、網膜円錐細胞を介して、脳に赤、緑及び青の
スペクトルバンドの平均値が送られるからである。又、
これが、電子デイスプレーが三色光放射システムを使用
する理由でもある。ＣＲＴカラー・モニターは赤、緑及
び青の燐光を使用する。多くの波長の色を測定するデバ
イス（例えば分光測光器及び放射測定器）は赤、緑及び
青のスペクトルバンドの荷重平均からＸＹＺ値を演算す
る。測色計、輝度測色計及び卓上デジタイザーのごとき
デバイスは、ＸＹＺ値が三色光フィルターを使用して直
接測定されるので、比較的費用が掛からない。各フィル
ターはスペクトルバンドの平均値を直接出す。即ち、色
彩は、これらのデバイスの場合、実質的に、３又は４種
類の波長のみを測定する。
【０００８】可能な場合は常に、カラー・イメージ・デ
ジタイザーは測色計に好ましいものである。これは、安
価で操作に余り特殊技術を必要としないからである。カ
ラー・イメージ・デジタイザーは、例えば画像捕捉のご
ときその他の機能のために必要と考えられているようで
ある。
【０００９】カラー・イメージ・デジタイザーは全スペ
クトル測定デバイスよりも正確性に欠けるが、我々は高
い正確性が必要でない場合だけを考えている。例えば、
ＣＲＴカラー・モニターのための画像フィードバック
は、例えば、製造業における製品色彩の品質管理ほど正
確性を必要としない。人間の色彩感覚は色の正確性から
のシステマチックな逸脱にすぐ順応する。
【００１０】ＸＹＺ値の使用はクーベルカ・ムンク・モ
デルの仮定を破る。これは、この派生した使用が各波長
での放射散乱を扱うからである。波長の荷重平均はこの
モデルにおいては物理的意味を持っていない。クーベル
カ・ムンク・モデルにＸＹＺ値を使用すると、着色料混
合体の予測を不正確なものにする。満足の行く予測を得
るためには、モデルに更に別の項を加える必要がある。
上述したＫ及びＳの式にはＫi Ｃi 及びＳi Ｃi の形の
項が含まれている。我々は、Ｋの式にＫijＣiＣj の形
の項を加えれば十分なことを見出だした。我々は、計数
Ｋijを着色料ｉ及びｊの光吸収合成係数と呼ぶ。この項
はＳには加える必要はない。これらの新しい係数は、一
般的に、混合媒体の中の着色料の間の分子間相互作用に
は関係しない。このような作用がある場合は、これらの
項を加えたほうが相互作用によく適応する。本発明にお
いては、我々の変形したクーベルカ・ムンクの式に最小
二乗法を適用することによって、鏡面反射、滑らかでな
い表面（例えば織物）、色の三刺激値測定、及び、簡単
なクーベルカ・ムンクの理論には含まれていないその他
の因子のごとき因子が部分的に補正される。
【００１１】着色料が基質に十分に混ざっておらず、基
質の上に薄い層としてある場合も、我々の技術は色をう
まく予測することが出来る。着色料がほとんど透明であ
ったり、余り光沢がないような場合は、着色料ｉに着色
料ｊを使用するときの計算にＫijＣi Ｃj の項を使用
し、又、着色料ｊに着色料ｉを使用するときの計算にＫ
jiＣi Ｃj の項を使用することが出来る。Ｋij及びＫji
の値は、着色料の不透明度に相関する程度には異なって
いる。明らかに、暗い着色料に加えられた明るい着色料
は、一般的に、明るい着色料に加えられた暗い着色料よ
りも明るく見える。この出願のその他の部分に対して
は、出願を単純化するために、この区別は不必要と仮定
する。着色料が十分に混合されていれば、着色料が不透
明であると否とに係わらず、ＫijはＫjiに等しい。
【００１２】吸収及び散乱係数を演算するのに必要な第
１の手順は、試料の測定値を集めることである。我々
は、無彩色基質から色のＸ，Ｙ、Ｚ三刺激測定値を測定
する。（以下に説明する全ての試料は同じ型の基質を用
いて調製される。使用する基質が異なる場合は、各基質
はそれぞれ別に処理されねばならない。）従って、色の
Ｘ，Ｙ、Ｚ三刺激測定値は、１つの色の異なる濃度の試
料から作られる。これは合成する全ての着色料について
行われ、濃度は、実際的限界の範囲内になければならな
い。最後に、色のＸＹＺ三刺激測定値がいく種類かの濃
度の何組かの着色料を使用した試料から作られ、この場
合、その合成濃度の合計が、関係項の中で１００％に決
められる一定範囲にあるようにする。本出願では全ての
濃度がパーセンテージで表される。この相対スケールは
着色料の重量又は体積の如き完全に物理的測定値に基づ
くものでなければならない。
【００１３】合計濃度の限界は、一般的に、着色料使用
工程における物理的制限によっている。例えば、織物の
ポリマーの中に拡散することの出来る着色料の量には上
限がある。製造目的のために選ばれた実際的限界が真の
物理的限界より小さい場合は、１００％を僅かに越えた
外挿値で十分予測される。
【００１４】次に、この第１の手順で得られた測定値の
利用が始まる。第２の手順は、無彩色基質からの測定値
を用い、無彩色基質の光吸収係数Ｋo を演算することで
ある。第３の手順は、単一の着色料によって彩色された
基質からの測定値を利用して、着色料i の光吸収係数Ｋ
i 及び光散乱係数Ｓi を演算する。最後の手順は、２種
類の着色料の合成測定値を用いて、各１対の着色料ｉ及
びｊの光吸収合成係数Ｋijを演算する。Ｘ，Ｙ及びＺ
（赤、緑、青）のスペクトル・バンドのためにこれらの
全ての係数が演算される。我々は、このモデルを高次の
項にまで拡げる必要のないことを見出だした。我々の方
法においては、この基質の光散乱係数が他の係数の中に
組み込まれているので、着色料基質にはＳo 項はないか
らである。
【００１５】Ｋo 、Ｋi 、Ｓi 及びＫijは全ての試料測
定値のコード化された合計値を表している。コンピュー
タの機能はこれらの係数を記憶するためには余り必要で
なく、これらの係数を得るために用いられた測定値の記
憶のため必要なのである。これらの記憶された係数は色
の予測のためのコンパクトなデータベースを含んでい
る。最小二乗法の適用によって、後の計算での試料測定
値の変動がなくなる。これは、色勾配を計算するために
これらの係数を使用することにより、常に、可視的に均
一な色シリーズが作られることを意味する。このこと
は、このような方法が不適当な場合、多くの色の測定値
に基礎を置く内挿法にとって重要な利点である。
【００１６】これらの係数を任意の配合のＫ／Ｓ値の計
算に用いれば（この場合、これらの値はＸＹＺ値に変換
される）、これらのＸＹＺ値をＲＧＢ色座標を計算する
ために使用し、電子ディスプレー上に合成結果を示すこ
とが出来る。
【００１７】本発明の利点の１つは、実際に合成着色料
の試料を作ることなく、基質上の合成着色料の色を予測
することが出来ることである。
【００１８】本発明の別の利点は、データを色の標準Ｘ
ＹＺ測定値に変換するために、イメージ・デジタイザー
を使用することが出来ることである。
【００１９】本発明の更に別の利点は、２以上の着色料
の合成結果を予測するのに、２以上の着色料で試料を作
る必要がないことである。
【００２０】本発明の更に別の利点は、鏡面反射、平滑
でない表面（例えば織物）、殆ど透明な着色料の層（例
えば、コンピュータ・ハードコピー着色料）、及び、色
の三刺激測定値が適用出来ることで、これらの条件は従
来のクーベルカ・ムンク・モデルには不適当なものであ
る。
【００２１】更にその他の利点は、非直線関係を直線化
することなく、Ｘ，Ｙ，Ｚ値をＲ，Ｇ，Ｂ値に、及びそ
の逆に変換する新規な方法であることである。この利点
は、ＣＲＴカラー・デイスプレーのＲＧＢ値及びイメー
ジデジタイザーから得られたＲＧＢ値に適用される。
【００２２】これらの利点及びその他の利点については
以下の説明で明らかにする。
【００２３】
【実施例】添付したフローチャート及びグラフを用いて
説明する。図１は、コンピュータ技術を用いて、合成着
色料を電子ディスブレー上に再現する基本的構成要素を
模式図的に示す図である。図の中の箱書きはコンピュー
タ情報モジュール（算法、データベース）又は外部測定
デバイス（測色計、分光測光器）を表す。矢印はコンピ
ュータ・モジュール又は外部測定デバイス間の情報の流
れる方向を示す。
【００２４】コンピュータ・エイデッド・デザイン又は
製造システムのための１つの色の測定スタンダードはＣ
ＩＥの色のＸＹＺ三刺激座標である。これは、カラー・
インプット（イメージ・デジタイザ、オンライン測色
計、着色料調合のデータベース）、及びカラー・アウト
プット（カラーモニター、カラープリンター、着色料調
合データベース）のいずれにも使用される。この標準表
色系は長く使われている国際標準である。ＸＹＺ値は色
を含むコンピュータ処理の基礎として段々使用されるよ
うになった。特定の適用に利点を持つ専門の表色系を持
つ会社でも、一般的に、彼等の色座標をＸＹＺ色座標に
変換することが出来る。ＸＹＺ値がその他の標準表色系
にも好まれる１つの理由は、これらが、イメージ・デジ
タイザー及び電子ディスプレーのごときデバイスにＲＧ
Ｂ表色系に直接該当するからである。
【００２５】測色計で測定されるＸＹＺ値は無次元で、
普通パーセンテージで表される。習慣として％表示は省
略される。ＸＹＺ測定値は幾つかの標準光源、例えばＤ
６５（６５，０００°Ｋの黒体温度における昼光）、に
対してなされる。これらの色の相対値は、光を出さない
物体の色を定量するために使用される。パーセンテージ
は、赤、緑及び青のスペクトルバンドの目標から反射さ
れる白色光のフラクションを示す。多くの物体はＣＲＴ
の色モニターの燐光のごとき光を放射する。絶対ＸＹＺ
値はキャンデラス／平方メートル又はフィート・ラムバ
ーツのごとき単位で測定される。これらを白色光の幾つ
かの標準エミッターに対して測ることによって、ＸＹＺ
の絶対値はＸＹＺの相対値に変換される。
【００２６】普通、測定値は色度ｘ及びｙと、「輝度」
Ｙとで表される。ＸＹＺ値とｘｙＹ値との間の変換は、
ｘ＝Ｘ／（Ｘ＋Ｙ＋Ｚ）、及びｙ＝Ｙ／（Ｘ＋Ｙ＋Ｚ）
で与えられる。Ｙは両座標系で同じである。色度座標ｘ
及びｙは無次元で、色の中の赤と緑の相対比率を表す。
青の相対量はｚ＝１−ｘ−ｙである。等量の赤と緑を持
つ色では、ｘ＝ｙ＝ｚ＝１／３である。Ｙが、しばし
ば、輝度（強度、明るさ）として選ばれるのは、人間の
視覚が赤や青よりも緑に敏感だからである。一般の色の
色度が図３に示されている。同時に、コンピュータ・エ
イドのデザインに用いられる代表的なＣＲＴカラーモニ
ターの色度の範囲が示されている。
【００２７】色の測定値は伝統的に測色計２又は分光測
光器３、又はクロマ・メータ（図示無し）から得られて
いるが、本発明の好ましい方法では、図１に示すごと
く、イメージ・デジタイザ１の使用が含まれる。シャー
プ電子会社（Ｍａｈｗａｈ、Ｎｅｗ Ｊｅｒｓｅｙ）で
作られているようなイメージ・デジタイザは、試料の上
に赤、緑及び青のレセプタ又は赤、緑及び青の閃光を持
っている。然し、従来のイメージデジタイザはＸＹＺ
値、又はその他の標準表色系をリターンしない。このよ
うなデバイスは、一般的に、赤、緑及び青の強度を示す
１組のＲＧＢ値をリターンする。赤、緑及び青のイメー
ジデジタイザーの色の成分はＲ、Ｇ及びＢでそれぞれ示
されている。ＲＧＢ値は同じ０から２５５の範囲にあ
る。本出願は、スキャナーのＲＧＢ値から相対ＸＹＺ値
への変換の方法を開示する。
【００２８】図１のイメージデジタイザ１が試料を走査
し、試料の各測定範囲に対するイメージデジタイザーの
ＲＧＢ値を出す。このＲＧＢのデータの演算がイメージ
デジタイザ・モデル算法５によって行われる。ＲＧＢ値
がＸＹＺ値に変換される前に、ある測定と演算とを予め
行わねばならない。先行作業の第１段階において、同じ
色度ｘｙを持つ目標灰色シェードがイメージデジタイザ
ーで走査される。これらの測定値は赤、緑及び青のイメ
ージディジタイザーの光測定チャンネルの特徴であるモ
デル・パラメータを演算するために使用される。第２段
階において、全ての目標の色がミクシング・マトリック
スを演算するために使用される。このミクシング・マト
リックスが、全てのイメージデジタイザのＲＧＢ値をＸ
ＹＺ値に変換するために使用される。
【００２９】前の段落で説明したごとく使用するため
に、目標カラーをイメージデジタイザーのために選択し
なければならない。目標は同じ特性を持つ灰色の幾つか
のシェードを持っていなければならない。又、主なる色
度に跨がる幾つかの色を持っていなければならない。色
が重要な工業のために、非常に念入りな目標カラーが開
発されている。このような目標カラー基準がアメリカ規
格協会（ＡＮＳＩ）ＩＴ８分科会によって、グラフィッ
ク・アートのために定められている。透明体（ＡＮＳ
Ｉ，ＩＴ８．７−１）及び反射体（ＡＮＳＩ，ＩＴ８．
７−２）の目標が１９９２年中頃、主だった写真フィル
ム・メーカーによって要求されるであろう。我々の目的
のためには、Ｍａｃｂｅｔｈ（登録商標） Ｃｏｌｏｒ
Ｃｈｅｃｋｅｒ（登録商標）の色翻訳チャートを使えば
十分である。
【００３０】Ｍａｃｂｅｔｈ ＣｏｌｏｒＣｈｅｃｋｅ
ｒの色翻訳チャートは、過去１５年にわたる写真及びビ
デオの信頼性の高い色基準でる。Ｍａｃｂｅｔｈチャー
トは表１に示すごとく、２４色で、ＣＩＥ光源Ｃに基づ
く色度座標を持っている。色の中には６つの灰色シェー
ドと、３つの加法混色の原色（赤、緑、青）と、３つの
減法混色の原色（黄、マジェンタ、シアン）と、２つの
肌色と、及び、１０の種々雑多な色とが含まれている。
チャートの色の色度は、図３に示すごとく、殆どの色彩
空間に跨がっている。図３は、又、主なＣＲＴのカラー
・モニターに使用できる特性の色域を示している。
【００３１】ＣＩＥの色測定値は、Ｍａｃｂｅｔｈチャ
ートの各色のためにＭａｃｂｅｔｈによって提供されて
おり、表１に示されている。この代わりに、特別の測色
計を使用して、Ｍａｃｂｅｔｈの色を測定することが出
来る。演算のためＭａｃｂｅｔｈによって与えられる測
定値を使用するのが便利であるが、好ましい方法ではな
い。標準測色計を選択して、全ての測定値がコンピュー
タ・エイデッド・デザインの役に立つようにしなければ
ならない。例えば、Ｈｕｎｔｅｒ Ｌａｂ ｏｆ Ｒｅ
ｓｔｏｎ（Ｖｉｒｇｉｎｉａ）製のＨｕｎｔｅｒＬａｂ
（登録商標）ＬａｂＳｃａｎが全ての商品試料の測定に
使用されるとすると、Ｍａｃｂｅｔｈカラー・チャート
の色はＨｕｎｔｅｒＬａｂ ＬａｂＳｃａｎで測定した
方が好ましい。次に、後の試料測定値とイメージデジタ
イザーの測定値との一致は出来るだけ近いほうがよい。
これは、たとえ同じ会社で作られているにしても、違っ
たモデルの測色計は違ったＣＩＥ色測定値を出すことが
あるからである。例えば、測色計によって鏡面反射の扱
い方が違う。又、測色計によって試料照明の幾何学的配
列が違う。この目標の色のＲＧＢ及びＸＹＺ値を集める
段階が図２に参照符号２０で示されている。この明細書
では、ＲＧＢ及びＸＹＺ値が０から１の間になるように
正規化される。この手順が図２に符号３０で示されてい
る。
【００３２】赤、緑、青のイメージデジタイザー・チャ
ンネルのための三刺激値をそれぞれＸr 、Ｙg 、及びＺ
b で表し、ここではこれらを「Ｘr Ｙg Ｚb 値」と呼
ぶ。これらは同じ中間色度を持つ灰色シェードのみを規
定している。これらが灰色シェード光測定値から演算さ
れ、イメージデジタイザーの光センサーの機能を特徴付
ける次の計算に使用される。灰色シェードはＭａｃｂｅ
ｔｈ ＣｏｌｏｒＣｈｅｃｋｅｒによって与えられる。
これらの灰色シェードは全て色度ｘ＝０．３１０及びｙ
＝０．３１６を持ち、表１に示す色１９〜２４である。
【００３３】イメージデジタイザーのＲＧＢ値とＸＹＺ
値との非直線性を特徴付けるイメージデジタイザー・パ
ラメータをそれぞれｇr 、ｇg 、及びｇb で表す。これ
らを総称して「ガンマ」又は「ｇ」と呼ぶ。
【００３４】赤、緑、青のチャンネル・コントラストを
特徴付けるイメージデジタイザー・パラメータをそれぞ
れＧr 、Ｇg 、Ｇb で表す。これらを総称して「ゲイ
ン」又は「Ｇ」と呼ぶ。
【００３５】赤、緑、青のチャンネルの明るさを特徴付
けるイメージデジタイザー・パラメータをそれぞれＯr
、Ｏg 、Ｏb で表す。これらを総称して「オフセッ
ト」又は「Ｏ」と呼ぶ。
【００３６】これらの符号及びガンマ、ゲイン、オフセ
ットと言う名称は、最初、陰極線管ディスブレー・デバ
イスを数学的にモデリングするために使われた言葉を引
用したものである。
【００３７】イメージデジタイザー・モデルにおけるゲ
イン、オフセット及びガンマは測定した色（ＸＹＺ三刺
激値）と、赤、緑、青のイメージデジタイザー・チャン
ネルのイメージデジタイザーの色の成分（ＲＧＢ値）と
の間の定量的関係を示し、次の式で示される： Ｘr ＝（Ｇr ＋Ｏr ＊Ｒ）gr 式 １．０ Ｙg ＝（Ｇg ＋Ｏg ＊Ｇ）gg 式 １．１ Ｚb ＝（Ｇb ＋Ｏb ＊Ｂ）gb 式 １．２ 次に、Ｍａｃｂｅｔｈ ＣｏｌｏｒＣｈｅｃｋｅｒの色
の測定値からイメージデジタイザー・モデル・パラメー
タの演算方法を説明する。ＸＹＺ測定値が、同じ色度ｘ
ｙを持つ灰色シェードのために選択される。色１９〜２
４（表１）は同じ色度を持つ灰色シェードである。これ
らがＸr Ｙg Ｚb 値として同定される。各灰色シェード
はそれに対応するイメージデジタイザーのＲＧＢ値を持
つ。この手順が図２に符号４０で示されている。次の手
順が図２に符号５０で示されており、Ｘ、Ｙ、又はＺデ
ータのいずれを扱うか選択する。普通は、Ｘ成分から初
め、次にＹ成分、そしてＺ成分の順で行う。この順番は
重要ではない。表２はＭａｃｂｅｔｈの色に対するイメ
ージデジタイザーのＲＧＢ値の１例を示す。式１．０、
１．１、及び１．２が次の形になる（ここでは、ｘ及び
ｙは色度の変数ではない）： ｙ＝（Ｇ＋Ｏ＊ｘ）g 式 ２．０ 式２．０はｙ1/g にリニヤーである１つの式に書き替え
られる。即ち、ｙ1/gの関数としてのｘのグラフは直線
である。
【００３８】 ｙ1/g ＝Ｇ＋Ｏ＊ｘ 式 ２．１ ｇ、Ｇ、及びＯの特定値のために、最小二乗誤差Ｅを規
定する。下付け文字ｍはここの灰色シェードの測定値ｘ
m （Ｒm 又はＧm 又はＢm を代表する）及びｙ1/g （Ｘ
rm 1/gr又はＹgm 1/gg又はＺbm 1/gbを代表する）を示す。
ｍが１からＮの範囲にあるＮ個の灰色シェードの測定値
の合計がある。
【００３９】
【数１７】 演算目的のためには、式２．２の展開形が好ましい。式
２．３は最終算法の中の先行演算手順の中に現れる。 Ｅ2 ＝Ｓyy−２ＧＳy −２ＯＳxy＋Ｏ2 Ｓxx＋２ＧＯＳx ＋Ｇ2 Ｓ 式 ２．３ 次の合計式が“Ｓ”可変値を決定する。
【００４０】
【数１８】
【００４１】
【数１９】 注：式２．４の変数Ｓは灰色シェードの番号に等しい。
【００４２】式２．１はｇが非直線である。この非直線
性のため、我々の試みは、同時にｇ、Ｇ、及びＯに対し
て式２．３の中の最小二乗誤差を最少にすることではな
く、その代わりに、ｇのリーズナブルの値をピックし、
最小二乗誤差を最少にするＧ及びＯを見出だすことであ
る。ｇを選択する手順が図２に符号６０で示されてい
る。このｇのために、Ｇ及びＯが式３．０、３．１及び
３．２から演算される。これが図２に符号７０で示され
ている。これらの式が、与えられたｇの最少二乗フィッ
トに対する解である。
【００４３】 Ｄ＝ＳxxＳ−Ｓx 2 式 ３．０ Ｇ＝（ＳxxＳy −Ｓx Ｓxy）／Ｄ 式 ３．１ Ｏ＝（−Ｓx Ｓy ＋ＳＳxy）／Ｄ 式 ３．２ 我々が選択したｇとＧ及びＯの演算した最適値とを式
２．３に代入し、図２に符号８０で示すごとく、最少二
乗誤差Ｅ（又はＥ2 ）を演算する。次に、ｇの別の値を
ピックし、これにより、別の１組のＧ及びＯの最適値を
演算し、新しい最少二乗誤差を得ることが出来る。この
別のｇの選択決定が図２に符号９０で示されている。２
つのｇから選択する場合、より良い値はより小さい最少
二乗誤差Ｅ（又は、Ｅ2 ）を与えるものである。最少二
乗誤差Ｅはｇが０．２５違っても余り敏感に変化しな
い。算法は増分０．２５における１から５までのｇの値
を選択する。概略を上述したごとく、各ｇに対して、式
２．４〜２．９、及び３．０〜３．２からＧ及びＯを、
又、式２．３から最少二乗誤差Ｅ2 を演算する。ｇの最
良値はＥ2 を最少にするものである。このｇの最終選択
が図２に符号１００で示されている。
【００４４】図２のフローチャートの符号１００に関し
て前に説明したごとく、この計算が灰色シェードの赤、
緑、青の成分に対して別々に繰り返される。この決定手
順が図２に符号１１０で示されており、これが、各スペ
クトル・バンドに対して、図２の手順６０、７０、８
０、９０及び１００の繰り返しに繋がる。
【００４５】ガンマ、ゲイン及びオフセットを演算する
と、式１．０、１．１及び１．２が、オリジナルの灰色
ーシェードにマッチする色度を持つ灰色シェードの
Ｘr 、Ｙg 及びＺb を予測する。Ｍａｃｂｅｔｈチャー
トに対しては、式１．０、１．１及び１．２が、ｘ＝
０．３１０、ｙ＝０．３１６に等しい色度を持つ灰色シ
ェードを全て正確に予測する。式１．０、１．１及び
１．２は、それ自体では、任意の色を正確に予測はしな
い。
【００４６】次の手順は、任意の色の測定値を用いるこ
とで、これにより任意の色のイメージディジタイザーの
ＲＧＢ値をＸＹＺ値に（又はその逆）変換することが出
来る。１つのＧｒａｓｓｍａｎの法則が、任意の色の三
刺激値を計算する近似式を提供する。これは、イメージ
ディジタイザーのチャンネルが十分に独立しているから
この法則を適用できるのである（即ち、チャンネルは一
次の色を測定する）。従って、付加的混色が充当され
る。ミクシング・マトリックスＭによって、灰色シェー
ドの三刺激値のリニヤーな組み合わせが行われる。任意
の色のＸＹＺ測定値の次のごときコラム・マトリックス
と、
【００４７】
【数２０】 式１．０、１．１及び１．２から演算されたＸＹＺの予
測値の次のごときコラム・マトリックスとを規定する。
【００４８】
【数２１】 次に、ＸとＸrgb とが直線関係で結ばれる。
【００４９】 Ｘ＝ＭＸrgb 式 ４．２ ここで、Ｍは３ｘ３マトリックスである。Ｍが灰色シェ
ードを含む全ての目標の色から、次の如き最少二乗（疑
似逆数）の形で、演算される。コラムが式４．０の中に
示される型のＮ個のカラー・ターゲットの測定値である
３ｘＮマトリックスＱを作る。即ち、 Ｑ＝｜Ｘ1 Ｘ2 Ｘ3 …ＸN ｜ 式 ４．３ 又、同じＮ測定値に対して式１．０、１．１及び１．２
から得られる対応するイメージディジタイザー・チャン
ネルの三刺激値の予測値を含む３ｘＮマトリックスＱ
rgb を作る。次のごとし、 Ｑrgb ＝｜Ｘrgb1Ｘrgb2Ｘrgb3…ＸrgbN｜ 式 ４．４ ＸＹＺ測定値及び全ての色の目標予測値からのＱ及びＱ
rgb の演算が図２に符号１２０で示されている。Ｑ及び
Ｑrgb のコラムが式４．２からマトリックスＭによって
繋がれ、従って、 Ｑ＝ＭＱrgb 式 ４．５ この式は、イメージディジタイザー・モデル（式１．
０、１．１及び１．２）が正確に正しく、色の測定値に
誤差が無ければ、成立する。モデルとデータが正確でな
いと、疑似逆数を用いて式４．５を解くことにより、
「ベスト」（最少二乗）の解が得られる。次のごとし、 Ｍ＝Ｑ（Ｑrgb Ｑrgb T ）-1 式 ４．６ 肩の文字Ｔはマトリックスの入れ替えを示し、肩文字−
１はマトリックスの逆数を示す。Ｍの演算が図２に符号
１３０で示されている。この逆数が３ｘ３のマトリック
スについて演算される。ゲイン、オフセット、ガンマ及
び混合マトリックスＭが解ければ、式１．０、１．１、
１．２及び４．２を用いて、イメージディジタイザー
Ｒ，Ｇ，及びＢを任意の色に対してＸ，Ｙ及びＺに変換
する。この手順は図２に符号１４０で示されている。こ
の変換は図１に符号５で示すイメージディジタイザー・
モデルの算法を構成している。ある適用例においては、
この変換が図１に符号７で示す色のデータベースのＸＹ
Ｚ値を得るために使うことが出来る。
【００５０】Ｘ，Ｙ及びＺをＲ，Ｇ及びＢに変換する逆
数変形が次のごとくに演算される。
【００５１】 Ｘrgb ＝Ｍ-1Ｘ 式 ５．０ Ｒ＝（Ｘr 1/gr−Ｇr ）／Ｏr 式 ５．１ Ｇ＝（Ｙg 1/gg−Ｇg ）／Ｏg 式 ５．２ Ｂ＝（Ｚb 1/gb−Ｇb ）／Ｏb 式 ５．３ ガンマ、ゲイン及びオフセットを１つのデバイスのため
に演算すると、これらが今後の全ての色の変換のために
コンピュータのデータ・ストラクチャの中にセーブされ
る。図２に符号１４０で示す手順さえ行えば、次の色の
変換が出来る。これらのイメージディジタイザーのパラ
メータはイメージディジタイザーによって測定される目
標の型に対して独立している。布地、紙製品、写真等の
ＸＹＺ値は同じパラメータの値を使ってＲＧＢ値から演
算することが出来る。
【００５２】表２はイメージディジタイザーのＲＧＢ値
の１例を示し、表３及び表４は対応する演算されたイメ
ージディジタイザー・モデルのパラメータを示す。詳し
く言うと、表３は、Ｒ，Ｇ及びＢのイメージディジタイ
ザー・チャンネルのガンマ、ゲイン及びオフセットを示
し、表４は３ｘ３ミクシング・マトリックスＭを示す。
これらのデータはＳｈａｒｐ（登録商標） ＪＸ−４５
０イメージディジタイザーから得られた。図４はＭａｃ
ｂｅｔｈの灰色シェードに対する測定及び予測されたＸ
ＹＺ三刺激値のグラフで、イメージディジタイザーのＲ
ＧＢ値に対してプロットされている。予測値はガンマ、
ゲイン及びオフセット値の最少二乗値から演算された。
図５はイメージディジタイザーのＸＹＺ測定値と全ての
Ｍａｃｂｅｔｈの色に対するＸＹＺ予測値との関係を示
す。予測値はガンマ、ゲイン、オフセット、及びミクシ
ングマトリックスＭの最少二乗値から演算された。
【００５３】非常に重要な本発明の第２の態様は、光学
的に厚い基質にブレンドされた着色料の色の予測能力で
ある。着色料が殆ど透明なものと、余り不透明でなもの
とで、これらが層状に施されている場合でも、本発明を
適用することが出来る。後者の場合の説明のための非限
定の例として、絨毯基質の上に複数の染料をスプレーし
た場合がある。層状に重ねられた着色料は、時に、着色
料の吹き付け順による色を作り出す場合があるが、この
影響は本発明によって対応することが出来る。このこと
は後程説明する。
【００５４】Ｋｕｂｅｌｋａ−Ｍｕｎｋ理論において
は、記号Ｒは、特定の光の波長で、試料から反射される
光のフラクションを現す。本出願では、Ｒは一般的に、
三刺激値Ｘ，Ｙ又はＺの１つのスケール・バージョンを
示す。前に説明したごとく、三刺激値座標Ｘ，Ｙ又はＺ
は赤、緑、青のスペクトルバンドの平均値をそれぞれ表
わしている。従って、我々がＲのためにＸＹＺ値を使用
するのはＫｕｂｅｌｋａ−Ｍｕｎｋ理論とは異なってい
る。ＸＹＺ値の選ばれたスケーリングは１以下のＲ値を
与えなければならない。
【００５５】着色料混合体の光吸収及び光散乱係数の比
はＫ／Ｓで示される。Ｋ／Ｓは無次元で、赤、緑、青の
スペクトルバンドに該当する３つの成分を持っている。
Ｋ／Ｓは式６．０によってＲ（スケールされたＸＹＺ
値）に関係する。式６．０は光学的に厚い基質に対する
Ｋｕｂｅｌｋａ−Ｍｕｎｋ理論から引き出される。式
６．１は式６．０の数学的に逆の式である。式６．０は
Ｒが測定によって既知のときＫ／Ｓを演算するために使
われる。式６．１はＫ／Ｓが既知のときＲを予測するた
めに使われる。
【００５６】 Ｋ／Ｓ＝（１−Ｒ）2 ／２Ｒ 式 ６．０ Ｒ＝（１＋Ｋ／Ｓ）−［（１＋Ｋ／Ｓ）2 −１］1/2 式 ６．１ 次の式７．０は、個々の着色料の吸収及び散乱係数がど
のようにして混合体の中で組み合わされ、この明細書に
おけるＫ／Ｓを作り出しているかを示す。
【００５７】
【数２２】 Ｋ0 は着色料を持たない基質の光吸収係数である。ゼロ
の下付き文字は着色料が０濃度であることを示す。後に
示す式８．０はＲ0 とＫ0 との間に関連を持たせる式で
ある。無彩色の基質の反射率はＲ0 によって示される。
この関係は、着色料の濃度比がＣi 及びＣj がゼロに設
定されたとき式７．０から引き出される。この表記方法
は標準的測色計とは若干異なり、基質の光散乱係数Ｓ0
が式７．０には見当たらない。基質の光散乱係数は本出
願においては他の係数の中に因子として組み込まれてお
り、無次元項「１」が式７．０の分母に出て来る理由を
説明する。この文書の中の全ての光の吸収及び散乱係数
は、この正規化の故に無次元である。
【００５８】光の吸収及び散乱係数Ｋi 及びＳi はＫｕ
ｂｅｌｋａ−Ｍｕｎｋのモデルの中の類似の項に該当す
る。着色料ｉ及びｊの光吸収の合成係数はＫijで示され
る。我々はこの発明においては、これらの後者の係数を
Ｋｕｂｅｌｋａ−Ｍｕｎｋのモデルに加え、ＸＹＺ値が
特定の波長におけるスペクトル反射率の代わりに用いら
れると言う事実を相殺する。一般的に、これらの係数は
混合媒体の中の着色料の間の分子間反応に関係しない。
但し、これらの係数は相互反応があるとすればこれらを
補正することが出来る。着色料が層状に施されている場
合は、下付け文字の順序が重要である。例えば、着色料
ｉが最初に基質に施され、これに着色料ｊが続く場合
は、Ｋijを使用し、この逆の場合は、Ｋjiを使用する。
一般的に、ＫijとＫjiとは等しくない。
【００５９】式７．１は式７．０の配列変えから引き出
される。これは着色料濃度Ｃi に対してリニヤーで、リ
ニヤーの最小二乗適合度計算の基礎を構成する。Ｋ／Ｓ
が等式の両側に現れる点に注意。適合度処理の間、Ｋ／
Ｓが式６．０から演算される。全ての計算において、濃
度は分数で示され、パーセンテージではない。
【００６０】
【数２３】 次に、試料の測定値が光の吸収、散乱、及び吸収合成係
数の演算に使用される方法を点検する。第１の手順は、
無彩色基質のＸ，Ｙ，Ｚ値の測定で、図６に符号１７０
で示されている。着色料が基質に無いときは、全ての着
色料濃度はゼロであり、Ｋ0 を式８．０から演算する。
Ｒ0 は基質のＸＹＺ測定値を表す。
【００６１】 （Ｋ／Ｓ）0 ＝Ｋ0 ＝（１−Ｒ0 ）2 ／２Ｒ0 式 ８．０ 混合体の中の全ての着色料の濃度が非常に小さいとき
は、この制約が基質の予測した色が無彩色基質の色に近
付くことを保証する。この計算が図６に符号２００で示
されている。図７の手順２参照。計算したならば、着色
料濃度が正確にゼロのときの基質の色を演算するために
Ｋ0 が用いられる。
【００６２】この方法の第２の手順は、１種類の着色料
が色々な濃度で存在する試料から色の測定値を得ること
である。これが図６に符号１８０で示されている。色は
１種類なので、式７．０及び７．１の中のＫijＣi Ｃj
項は全てゼロである。Ｋ1 とＳ1 は着色料１の濃度Ｃ1
における光吸収及び散乱係数を示す。Ｋ／Ｓが式６．０
から計算され、Ｋ0 は既知である。式９．０から９．５
に示す最小二乗適合度を用い、Ｋ1 及びＳ1 を演算す
る。
【００６３】 ｙ＝Ｋ／Ｓ−Ｋ0 式 ９．０ ｘ1 ＝Ｃ1 式 ９．１ ｘ2 ＝（Ｋ／Ｓ）Ｃ1 式 ９．２ ｗ＝［Ｒ／（Ｒ0 −Ｒ）2 ］2 式 ９．３
【００６４】
【数２４】 変数ｙは従属変数を示し、ｘ1 及びｘ2 は独立変数を示
す。変数ｗは最小二乗ウエイトホーシング相対誤差で、
ＸＹＺに対して均一である。インデックスｑは異なる着
色料濃度に関し、Ｐは異なる着色料濃度の総数（試料の
総数）を示す。この計算が図６に符号２１０で示されて
いる。図７の手順３参照。計算が終わったならば、Ｋ1
及びＳ1 が着色料１の異なる色によって作り出される色
を予測するために使われる。
【００６５】この方法の第３の手順は、異なった濃度で
存在する２種類の着色料で作られた色測定値を得ること
である。この２つの着色料を下付け文字１及び２で表す
ことにする。２つの色濃度Ｃ1 及びＣ2 は、最小二乗適
合度に使用されるデータ・セットの中でゼロであっては
ならない。濃度計Ｃ1 ＋Ｃ2 は１００％でなければなら
ない。これらの濃度が出来るだけ広い範囲に跨がるのが
最善で、Ｃ1 ＋Ｃ2 ＝１００％を選ぶのに便利である
（但し必ずしも必要ではない）。Ｋ0 、Ｋ1 、Ｓ1 、Ｋ
2 及びＳ2 は着色料１及び２について既知とする。普
通、Ｋ／Ｓは式６．０を使用して計算される。Ｋ12を決
定する。何種類かの濃度で施された１対の着色料を持つ
基質のＸ，Ｙ，Ｚ値を測定する手順が図６に符号１９０
で示されている。
【００６６】 ｙ＝（１＋Ｓ1 Ｃ1 ＋Ｓ2 Ｃ2 ）（Ｋ／Ｓ）−（Ｋ0 ＋Ｋ1 Ｃ1 ＋Ｋ2 Ｃ2 ） 式 １０．０ ｘ＝Ｃ1 Ｃ2 式 １０．１ ｗ＝［Ｒ／（Ｒ0 −Ｒ）2 ］2 式 １０．２
【００６７】
【数２５】 インデックスｑは異なる着色料の濃度に関し、Ｐは合成
の組数の計（試料の総数）を示す。この計算が図６に符
号２２０で示されている。図７の手順４参照。計算した
後、Ｋ12（Ｋ0 、Ｋ1 、Ｓ1 、Ｋ2 及びＳ2 と共に）が
着色料１及び２の異なった配合により作られる色の予測
に使われる。
【００６８】同じ手順を用いて、別の着色料のパラメー
タを決定する。１回決定すれば、Ｋ0 、Ｋi 、Ｓi 、Ｋ
ijは、式７．０及び６．１により合成された任意の着色
料のＸＹＺ値の演算に使われる。これが図６に符号２３
０で示す最終手順である。
【００６９】１つの所定の基質の上の１つの着色料の色
は９つの数字：Ｋ0 、Ｋ1 、Ｓ1 、によって規定され
る。所定の基質の上の着色料は、１８の数字：Ｋ0 、Ｋ
1 、Ｓ1 、Ｋ2 、Ｓ2 、Ｋ12によって規定される色を持
つ。（本出願の初めの部分で説明したごとくＫ21が必要
な場合は除く。）所定の基質上の３つの着色料の色は３
０の数字Ｋ0 、Ｋ1 、Ｓ1 、Ｋ2 、Ｓ2 、Ｋ3 、Ｓ3 、
Ｋ12、Ｋ13、Ｋ23によって規定される。これらの係数が
全ての試料の測定値を集約し、普通、次の合成計算のた
めにコンピュータのデータペースに保存される。例えば
図１の７。
【００７０】この型の着色料合成分析の１つの適用例
は、コンピュータ制御されたジェット染色技術を使って
染料を絨毯の上に施す場合である。いかなる光学的に厚
い基質もこの方法を使用することが出来るので、上記例
に限定して考えるべきではない。測定は絨毯が最終製品
の形になった後に行われねばならない。例えば、絨毯の
染料が固定された後、裁断された後、主だった処理が行
われた後、等々。この条件には非染色の試料について行
われる測定を含み、これが絨毯基質のＫ0 を演算するた
めに使われる。測色計がイメージディジタイザーの代わ
りに使われる場合は、測定値を測色計だけに限定した方
が良い。測色計がＣＩＥのＬ* ａ* ｂ* 測定値を提供す
る場合は、これらの色座標は、適宜の測色計装置を使っ
て、ＸＹＺ三刺激値に変換しなければならない。この測
色計は可能な最大アパーチャ、好ましくは少なくとも直
径２”、を使ってキャリブレーションしなければならな
い。ガラスは、絨毯のパイルをガラスに押し付けると、
普通の使用状態では判らない光沢が出るので、測色計の
アパーチャには使われない。試料上の絨毯パイルは、そ
れが好ましい方向に寝るように測定前に手でならす。
【００７１】図７は「染料合成計算のための段階的回帰
手順」と銘打った４つのグラフによる実施例からなる。
着色料の光吸収、散乱、及び合成の係数を得るために必
要とする計算をグラフ的に示したものである。この例で
は、合成着色料が絨毯に施された２種類の染料から成り
立っているとする。染料１及び２が、それぞれ、相対濃
度Ｃ1 及びＣ2 で合成される。Ｃ1 、Ｃ2 及びＲは直交
座標軸で示される。
【００７２】図７の手順１参照。前にも述べたごとく、
比Ｋ／Ｓの使用はＲよりもコンピュータ的に単純であ
る。Ｋ／ＳとＲとの数学的関係が図７の手順１に示され
ている。全ての計算予測値が非染色又は染色試料の測定
から得られる。１つの試料は着色料が施されておらず
（Ｃ1 ＝Ｃ2 ＝０）、「ゼロ染料」と示した点である。
この点は最も明るい測定色を持ち、この測定値がＲ方向
の一番高いデータ点として示されている。
【００７３】幾つかの試料は異なる量の染料１を持たね
ばならない。これらは点「３３％染料１」、「６６％染
料１」及び「１００％染料１」で、これらは全てＣ1 及
びＲ軸で形成される面の中に横たわっている。これらの
点は「ゼロ染料」測定値より低い（暗い）。染料１の量
が増えるのにつれて、測定値はＣ1 軸に向かって近付く
（より暗くなる）。異なる量の染料２を持つ試料につい
ても同じことが言える。この試料においては、染料２の
Ｒの測定値が染料１よりも小さい。これは色成分Ｒに対
して染料１が染料２よりも明るい色であることを意味す
る。
【００７４】手順１は、又、染料１と２との合成割合の
違うものについての３種類の測定値を示している。１つ
の測定値だけが「５０％染料１＋５０％染料２」で示さ
れている。Ｃ1 及びＣ2 軸によって形成される面に１本
の破線が横たわっている。この破線上の点は 式Ｃ1 ＋
Ｃ2 ＝１００％ を満足する。この式を満足するように
合成試料を作ると便利である。３つの２染料合成色が図
面の中でこの関係を満足しており、従って、この三次元
空間の中で破線上に横たわっている。例えば、他の２つ
の点は「３３％染料１＋６７％染料２」と「６７％染料
１＋３３％染料２」である。従って、合成濃度の合計は
いずれの場合も１００％である。
【００７５】要約すると、図面の座標系をＲ軸に沿って
下に向かって見て行くと、１染料の測定値はＣ1 又はＣ
2 軸上に投影し、２染料の測定値は上記破線の上に投影
する。制約条件 Ｃ1 ＋Ｃ2 ≦１００％ を満足する任
意の合成はこれら３本の線で囲まれる三角形の中に落ち
るであろう。
【００７６】図７の手順２参照。一連の回帰の中での第
１段階は、着色料の無い試料を使ってＫ0 を演算する。
着色料の無い１００％濡れた溶液がしばしば無彩色基質
（例えば生機（きはた）織物）に施され、たとえ色を与
えたり、外観を変えても、混合絨毯試料の１部として処
理されることがある。この「クリヤーな」溶液は着色料
を基質上で合成するときにしばしば使われる。即ち、こ
れは、着色料の希釈及び二色合成試料作成の間に使われ
る。
【００７７】ゼロ着色料試料の測定値が基質の光吸収計
数Ｋ0 を演算するのに使われる。この時点で計算を止
め、Ｋ／Ｓ＝Ｋ0 を使って染料１及び２の合成体の色
を予測しようとすると、この予測値は、Ｃ1 ＝Ｃ2 ＝０
のときの測定値にマッチするだけである。予測値が、
図７の手順２の中に、「ゼロ染料」点と交差する三角形
の予測面によって示される。この回帰の第１段階におい
ては、「ゼロ染料」予測値のみが測定値と一致する。
【００７８】図７、手順３参照。この手順は着色した試
料を使用し、Ｋ1 、Ｓ1 、Ｋ2 及びＳ2 を演算する。各
着色料が少なくとも５種類の異なった濃度を持つように
したほうが良い、図７には３種類しか示されていない。
０％濃度はこの手順では使用されないが、１００％濃度
は含んでいる。考えられる組合わせは１５％、２０％、
３０％、５０％、及び１００％である。等しくない濃度
増分を選択すると、視覚的に均一な色勾配を作り出しや
すい。少量の着色料が最終の色に強い影響を与える。１
０種類の異なった濃度が使える場合、考えられる組合わ
せは１５％、２０％、３０％、３５％、４０％、５０
％、６０％、８０％、及び１００％である。
【００７９】１つの着色料の希釈測定値が光吸収及び散
乱係数の演算に使われる。これらの新しい項が、手順３
に示すごとく、Ｋ／Ｓの式に加えられる。新しい予測面
が実線で示されている。新しい予測値は単一の染料測定
値に近接して横たわっている。これらは、演算が最小二
乗適合度のため、完全には一致しない。適合度の良さは
測定精度及び我々の数学モデルの自由度によって限定さ
れる。然し予測値は「ゼロ染料」測定値とは一致する。
回帰の第２段階は最初の回帰段階を邪魔することはない
が、予測面は染料の合成測定値とは十分には一致しな
い。
【００８０】図７、手順４参照。この最終手順は１対の
着色料を使って、Ｋijを確立する。各着色料が少なくと
も５種類の異なった濃度を持つようにしたほうが良い、
図７には３種類しか示されていない。０％及び１００％
濃度は使ってはいけない。考えられる組合わせは１５％
／８５％、３３％／６７％、５０％／５０％、６７％／
３３％、及び８５％／１５％である。１０種類の異なっ
た濃度が使える場合、考えられる組合わせは１５％／８
５％、２０％／８０％、２５％／７５％、３０％／７０
％、４０％／６０％、６０％／４０％、７０％／３０
％、７５％／２５％、８０％／２０％、及び８５％／１
５％である。等しくない濃度増分を選択すると、視覚的
に均一な色勾配を作り出しやすい。少量の着色料が最終
の色に強い影響を与える。
【００８１】着色料の１対の合成測定値を使って、光の
吸収合成係数を演算する。手順４に示すごとく、これが
Ｋ／Ｓの式に１つの項を加える。この予測曲面は全ての
染料測定値を十分に予測する。再び、演算が最小二乗適
合度のため、最終予測値は完全には一致しない。この第
３の回帰段階はゼロ染料又は１対の染料合成予測値には
影響しない。
【００８２】要約すると、変形されたＫｕｂｅｌｋａ−
Ｍｕｎｋモデルが段階における測定値にフィットされ、
各段階がより多くの測定データを含み、更に、全体の最
小二乗誤差を小さくする。前の予測値の精度による影響
は無い。このデカップリングの代数的理由は、Ｋ／Ｓの
等式における分子Ｋ及び分母Ｓの式が着色料濃度の逐次
高くなるオーダーの積を含むと言うことである。若し、
両方の着色料の濃度がゼロの場合は、Ｃ1 、Ｃ2 及びＣ
1 ＊Ｃ2 を持つ全ての項が消える。残る項は手順２に示
すものだけである。１つの着色料がゼロとすると、Ｃ1
＊Ｃ2 を持つ項が消える。残る項は手順３に示すものだ
けである。両方の濃度がゼロでなく、従って、Ｃ1 ＊Ｃ
2 の項がある場合のみ、手順４に示す全ての項が適用さ
れる。
【００８３】合成着色料の予測値は図の手順４に示す曲
面の上に落ちる。予測値は全表面にわたって我々が特定
した仕様を満足する。任意の予測値は、図示した式（式
７．０と同じ）で作られたＫ／Ｓモデルに基づく内挿値
である。３種類の着色料が含まれる場合は、内挿範囲は
体積、等である。この数学的方法は高次の項（例えばＳ
12、Ｋ123 、Ｓ123 等）に拡大することが出来るが、こ
れは三刺激値座標予想には必要無いことが判った。
【００８４】第３の着色料の追加はＫ0 、Ｋ1 、Ｓ1 、
Ｋ2 又はＫ12の値に影響しない。但し、第３の着色料を
含む合成体の予測値が必要な場合は、係数Ｋ3 、Ｓ3 、
Ｋ13及びＫ23は計算しなければならない。再び、染料１
及び２の合成体の前の全ての予測値は、たとえ、Ｋ／Ｓ
の等式に別の項が加えられたとしても、影響されずに残
る。
【００８５】表５から７は、７種類の織物染料の光吸
収、光散乱及び合成光吸収の各係数のリストを示す。こ
れらの染料はコンピュータ制御の下での染料スプレーに
よって絨毯に施された。３つの暗い（「濃い」）染料
と、３つの中間の（「薄い」）染料と、１つの明るい染
料（黄色）とがある。４種類の染料（薄い赤、薄い緑、
薄い青及び黄色）の部分集合体の合成予測値が表８〜１
１に示されている。
【００８６】本発明の第３及び最後の態様は、計算され
た合成色をＲＧＢベースの電子ディスプレーに表示する
ことである。コンピュータ化された設計システムが、色
を陰極線管（ＣＲＴ）上に表示するのが普通である。Ｒ
ＧＢ値を使った任意の電子デイスプレー、例えば液晶デ
イスプレー、をこのために使うことが出来る。陰極線管
（ＣＲＴ）は光を三原色で放射する。これが人間の網膜
の赤、緑、青の感覚器官を刺激する。便宜上、ＲＧＢベ
ースの電子ディスプレーはエレクトロルミネッセンス・
ディスプレー又はプラズマ・ディスプレーとして使うこ
とができるものであるが、本出願ではＣＲＴディスプレ
ーを扱うことにする。電子ディスプレーによって放射さ
れる色はラジオメータ又はクロマメータ（例えば、Ｍｉ
ｎｏｌｔａ（登録商標） ＴＶカラー分析器II、Ｍｉｎ
ｏｌｔａ ＣＲＴカラー分析器ＣＡ−１００、Ｍｉｎｏ
ｌｔａ クロマメータＣＳ−１００）によって測定され
る。色の測定デバイスに使うことの出来る各種の色の座
標の中から、前に述べたごとく、本出願ではＸＹＺ三刺
激値を選択する。
【００８７】光を放射するデバイスから測定されるＸＹ
Ｚ値は例えばカンデラ／平方メートルのごとき単位を持
っている。これらの測定値は絶対であると言われるもの
で、ＸＹＺ相対値は以後Ｘ′Ｙ′Ｚ′で示す。これらの
値は無次元であり、普通パーセントで示される。同じ白
色媒体（同じ色度を持つ）の異なった絶対測定値は、普
通、電子可視ディスプレー又はコンピュータ・イメージ
ング・ハードコピーの適用に対しては、同じ相対値で示
される。ディスプレー・デバイスに対して可能な最大Ｙ
値を選択して、絶対ＸＹＺ測定値を相対Ｘ′Ｙ′Ｚ′値
に変換することが多い。例えば、あるＣＲＴに可能な最
大Ｙ値が８０cd/m2 で、１つの表示された色のメジャー
がＸ＝６０cd/m2 、Ｙ＝７０cd/m2 、Ｚ＝７５cd/m2 で
あったとすると、これに対応する相対Ｘ′Ｙ′Ｚ′値は
Ｘ′＝７５．００、Ｙ′＝８７．５０、及びＺ′＝９
３．７５である。これらはパーセントだが、慣例として
パーセント表示は省略される。
【００８８】明細書の初めに説明したごとく、色の測定
値を色度ｘ，ｙ及び「輝度」Ｙで表すことが一般的であ
る。ＸＹＺ値とｘｙＹ値と間の変換は明細書の初めに説
明したごとくにして行われる。Ｘ及びＹのディメンショ
ンを持った値が、これによって、ｘ，ｙのディメンショ
ンの無い値に変換される。コンピュータエイド設計のた
めの代表的ＣＲＴカラー・モニターの色度範囲が図３に
示されており、Ｍａｃｂｅｔｈ ＣｏｌｏｒＣｈｅｃｋ
ｅｒの色と比較される。
【００８９】従来の電子カラー・ディスプレーは色を表
示するのに直接ＸＹＺ値を使用しない。このようなデバ
イスは一般的に、赤、緑、青の燐光放射を示すＲＧＢ値
を使用する。赤、緑、青の色の成分が、それぞれ、Ｒ、
Ｇ、Ｂで表される。ＲＧＢ値の一般的範囲は０から２５
５である。次に、絶対ＸＹＺ三刺激値をＲＧＢ電子ディ
スプレー値に変換する方法を示す。図１はＣＲＴディス
プレーを符号４で示し、これが色をＲＧＢ値として授受
する。ＸＹＺ値とＲＧＢ値との間の変換が図１の中の符
号６、ＣＲＴモデル算法によって行われる。
【００９０】コンピュータ演算の第１段階において、
赤、緑、青のＣＲＴ燐光に対して別々の測定が異なった
輝度で行われる。これは幾つかのレベルのＲＧＢ値につ
いて行われねばならない。燐光測定値に対するＲＧＢ値
の代表的な１組は５０，１００，１５０，２００，２２
５及び２５５である。信頼性のある測定値を提供する最
小のＲＧＢ値を含むようにしたほうが良い。これは、Ｍ
ｉｎｏｌｔａ ＴＶ−カラー分析器IIに対してはＹ＝
０．６cd/m 2 で、Ｍｉｎｏｌｔａ色彩分析器ＣＡ−１０
０に対してはＹ＝０．３cd/m 2 である。このデータ収集
手順が図８に符号２６０で示されている。本明細書にお
いては、ＲＧＢ及びＸＹＺ値が正規化され、各値が０か
ら１の間になるようにしている。この手順が図８に符号
２７０で示されている。これらの測定値を使用して、
赤、緑、青の燐光に特徴的であるモデル・パラメータを
演算する。演算の第２段階において、混合マトリックス
を演算し、１つの燐光がオンのとき作られる色だけでな
く、任意の色が変換されるようにする。
【００９１】本明細書の初めに、イメージディジタイザ
ーのモデル算法（図１の５）の説明をした。電子カラー
・ディスプレーのためのＣＲＴのモデル算法（図１の
６）はこれに非常に似ている。事実、ＣＲＴのモデル算
法は本発明のディジタイザーのモデル算法を作るために
ディジタイザーの色測定に初めて適用された。ＣＲＴの
モデル算法において使用される用語（例えば、ガンマ、
ゲイン、オフセット）はＣＲＴの内部エレクトロニック
スと関連している。例えば、エレクトロン・ビーム・イ
ンテンシティー及び増幅器電圧である。本出願の説明に
当たって、モデルの数学的態様は、例えばイメージディ
ジタイザーやその他の電子カラー・ディスプレー・デバ
イスのごときＲＧＢベースのデバイスに適用することが
出来る。
【００９２】赤、緑、青のＣＲＴ燐光の三刺激値が
Ｘr 、Ｙg 、Ｚb でそれぞれ表され、全体的には「Ｘr
Ｙg Ｚb 値」と呼ばれる。Ｘr は赤の燐光のみが点灯し
ているときに測定されるＸ値であり（ＧとＢはゼロ）、
Ｙg は緑の燐光のみが点灯しているときに測定されるＹ
値であり（ＲとＢはゼロ）、Ｚb は緑の燐光のみが点灯
しているときに測定されるＺ値である（ＲとＧはゼ
ロ）。
【００９３】１つのＣＲＴに対して、ガンマは、電子ピ
ーム加速電圧とその結果の色の輝度との非直線関係を特
徴付けるパラメータである。赤、緑、青の燐光のガンマ
値が、それぞれ、ｇr 、ｇg 、及びｇb で示される。こ
れらを集合して「ガンマ」又は「ｇ」と呼ぶ。
【００９４】１つのＣＲＴに対して、ゲインは、結果と
して得られた色の知覚コントラスト・レベルを特徴付け
るパラメータである。赤、緑、青の燐光のゲイン値が、
それぞれ、Ｇr 、Ｇg 、及びＧb で示される。これらを
集合して「ゲイン」又は「Ｇ」と呼ぶ。
【００９５】１つのＣＲＴに対して、オフセットは、結
果として得られた色の知覚輝度を特徴付けるパラメータ
である。赤、緑、青の燐光のオフセット値が、それぞ
れ、Ｏr 、Ｏg 、及びＯb で示される。これらを集合し
て「オフセット」又は「Ｏ」と呼ぶ。
【００９６】ＣＲＴモデル算法におけるゲイン、オフセ
ット及びガンマのパラメータは、各燐光に対する、測定
した色（絶対ＸＹＺ三刺激値）とＣＲＴ色彩座標の色成
分（ＲＧＢ値）との間の定量的関係を規定する。これが
次の等式に示されている。これらの等式は、イメージデ
ィジタイザーのモデル算法について本明細書の初めに提
出したものと同じである。
【００９７】 Ｘr ＝（Ｇr ＋Ｏr ＊Ｒ）gr 式 １．０ Ｙg ＝（Ｇg ＋Ｏg ＊Ｇ）gg 式 １．１ Ｚb ＝（Ｇb ＋Ｏb ＊Ｂ）gb 式 １．２ 次に、燐光の色測定値からＣＲＴモデル・パラメータを
演算する方法の説明をする。Ｘr が赤の燐光のみが点灯
したときに測定され（Ｇ及びＢはゼロ）、Ｙgが緑の燐
光のみが点灯したときに測定され（Ｒ及びＢはゼロ）、
Ｚb が青の燐光のみが点灯したときに測定される（Ｒ及
びＧはゼロ）。この第１の手順は、図８に符号２８０で
示されているごとく、次の演算のために、赤、緑、青の
燐光のいずれか１つを選択することである。等式１．
０、１．１及び１．２が次の形になる（ｘ及びｙは色度
変数ではない）。即ち、 ｙ＝（Ｇ＋Ｏ＊ｘ）g 式 ２．０ 式２．０はｙ1/g にリニヤーである１つの式に書き替え
られる。即ち、ｙ1/gの関数としてのｘのグラフは直線
である。
【００９８】 ｙ1/g ＝Ｇ＋Ｏ＊ｘ 式 ２．１ ｇ、Ｇ、及びＯの特定値のために、最小二乗誤差Ｅを規
定する。下付け文字ｍはここの灰色シェードの測定値ｘ
m （Ｒm 又はＧm 又はＢm を代表する）及びｙ1/g （Ｘ
rm 1/gr又はＹgm 1/gg又はＺbm 1/gbを代表する）を示す。
Ｎ個の灰色シェード測定値の合計があり、従って、ｍは
１からＮの範囲にある。
【００９９】
【数２６】 演算目的のためには、式２．２の展開形が好ましい。式
２．３は最終算法の中の先行演算手順の中に現れる。 Ｅ2 ＝Ｓyy−２ＧＳy −２ＯＳxy＋Ｏ2 Ｓxx＋１ＧＯＳx ＋Ｇ2 Ｓ 式 ２．３ 次の合計式が“Ｓ”可変値を決定する。
【０１００】
【数２７】 式２．４の変数Ｓは灰色シェードの番号に等しい。式
２．１はｇが非直線である。この非直線性のため、我々
の試みは、同時にｇ、Ｇ、及びＯに対して式２．３の中
の最小二乗誤差を最少にすることではなく、その代わり
に、ｇのリーズナブルの値をピックし（図８の２９
０）、最小二乗誤差を最少にするＧ及びＯを見出だす
（図８の３００）。このｇのために、Ｇ及びＯが式３．
０、３．１及び３．２から演算される。これらの式が、
与えられたｇの最少二乗フィットに対する解である。
【０１０１】 Ｄ＝ＳxxＳ−Ｓx 2 式 ３．０ Ｇ＝（ＳxxＳy −Ｓx Ｓxy）／Ｄ 式 ３．１ Ｏ＝（−Ｓx Ｓy ＋ＳＳxy）／Ｄ 式 ３．２ 我々が選択したｇとＧ及びＯの演算した最適値とを式
２．３に代入し、図８に符号３１０で示すごとく、最少
二乗誤差Ｅ（又はＥ2 ）を演算する。次に、ｇの別の値
をピックし、これにより、別の１組のＧ及びＯの最適値
を演算し、新しい最少二乗誤差を得ることが出来る。こ
の手順が図８に符号３２０で示されている。２つのｇの
推計値に対して、より良い値はより小さい最少二乗誤差
Ｅ（又は、Ｅ2 ）を与えるものである。最少二乗誤差Ｅ
はｇが０．２５違っても余り敏感に変化しない。算法は
増分０．２５における１から５までのｇの値を選択す
る。概略を上述したごとく、各ｇに対して、式２．４〜
２．９、及び３．０〜３．２から、Ｇ及びＯを、又、式
２．３から最少二乗誤差Ｅ2 を演算する。ｇの最良の値
はＥ2 を最少にするものである。このｇの最終選択が図
８に符号３３０で示されている。
【０１０２】図８のフローチャートの符号２８０に関し
て前に説明したごとく、この計算が対応する測定された
燐光色のの赤、緑、青の成分に対して別々に繰り返され
る。この決定手順が図８に符号３４０で示されており、
これが、各ＣＲＴ燐光に対して、図８の手順２９０、３
００、３１０、３２０及び３３０を繰り返えす。
【０１０３】ガンマ、ゲイン及びオフセットのパラメー
タを演算したならば、式１．０が、最小誤差をもって、
Ｇ及びＢがゼロの時のＸr をＲから予測する。緑及び青
の燐光に対しても同様である（１つの燐光のみがオン
で、他の２つはオフ）。式１．０、１．１及び１．２
は、それ自体では、任意の色を正確に予測はしない。
【０１０４】次の手順は、すべての燐光測定値を使用し
て、任意の色のＲＧＢ値をＸＹＺ値に（又はその逆）変
換することである。１つのＧｒａｓｓｍａｎの法則が、
任意の色の三刺激値を計算する近似的方法を提供する。
これは、ＣＲＴモニターが付加的色彩混合によって色を
作り出すので、適用可能である。ミクシング・マトリッ
クスＭが燐光の三刺激値の直線組み合わせを作り出す。
任意の色の測定された三刺激値の下記のコラム・マトリ
ックスと、
【０１０５】
【数２８】 式１．０、１．１及び１．２から演算されたＸＹＺの予
測値の下記のコラム・マトリックスとを規定する。
【０１０６】
【数２９】 次に、ＸとＸrgb とが直線関係で結ばれる。
【０１０７】 Ｘ＝ＭＸrgb 式 ４．２ ここで、Ｍは３ｘ３マトリックスである。Ｍが、次のご
とく、最少二乗（疑似逆数）の形で、ゲイン、オフセッ
ト及びガンマの各パラメータを演算するために使われる
全ての測定値から演算される。そのコラムが式４．０の
中に規定される種類の３Ｎ燐光測定値Ｘm である３ｘ３
ＮマトリックスＱを規定する。
【０１０８】 Ｑ＝｜Ｘ1 Ｘ2 Ｘ3 …ＸN ｜ 式 ４．３ 又、同じ３Ｎ個の測定値に対して式１．０、１．１及び
１．２から得られる対応する予測値を含む３ｘ３Ｎマト
リックスＱrgb を規定する。
【０１０９】 Ｑrgb ＝｜Ｘrgb1Ｘrgb2Ｘrgb3…ＸrgbN｜ 式 ４．４ ＸＹＺ測定値及び全ての燐光色測定値からのＱ及びＱ
rgb の演算が図８に符号３５０で示されている。Ｑ及び
Ｑrgb のコラムが式４．２のマトリックスＭによって結
ばれ、従って、 Ｑ＝ＭＱrgb 式 ４．５ この式は、ＣＲＴモデル（式１．０、１．１及び１．
２）が正確に正しく、色の測定値に誤差が無ければ、成
立する。モデルとデータが正確でないと、疑似逆数を用
いて式４．５を解くことにより、「ベスト」（最少二
乗）の解が得られる。次のごとし、 Ｍ＝Ｑ（Ｑrgb Ｑrgb T ）-1 式 ４．６ 上付き文字Ｔはマトリックスの入れ替えを示し、上付き
文字−１はマトリックスの逆数を示す。この手順が図８
に符号３６０で示されている。この逆数計算が３ｘ３の
マトリックスについて行われる。ゲイン、オフセット、
ガンマ及び混合マトリックスＭが解ったならば、式１．
０、１．１、１．２及び４．２を用いて、ＣＲＴのＲ，
Ｇ，及びＢを任意の色に対してＸ，Ｙ及びＺに変換す
る。Ｘ，Ｙ及びＺをＲ，Ｇ及びＢに変換するこの逆数変
換が次のごとくに演算される。
【０１１０】 Ｘrgb ＝Ｍ-1Ｘ 式 ５．０ Ｒ＝（Ｘr 1/gr−Ｇr ）／Ｏr 式 ５．１ Ｇ＝（Ｙg 1/gg−Ｇg ）／Ｏg 式 ５．２ Ｂ＝（Ｚb 1/gb−Ｇb ）／Ｏb 式 ５．３ 図９は、コンピュータ・エイドの設計に用いられる代表
的ＣＲＴの赤、緑、青の燐光のＣＲＴモデル計算のグラ
フを示す。ＸＹＺ三刺激値がＲＧＢ値に対してプロット
されている。図１０は全ての燐光色に対するＣＲＴ測定
値と予想値との対比を示す。ＸＹＺ予想値がＸＹＺ測定
値に対してプロットされている。
【０１１１】本実施例においては、絶対色測定値（ＸＹ
Ｚ三刺激値）と電子カラー・ディスプレーの成分（ＲＧ
Ｂ値）との間の変換手段が設けられる。絶対色測定値
（ＸＹＺ三刺激値）と相対色測定値（Ｘ′Ｙ′Ｚ′三刺
激値）との間の変換も必要である。電子ディスプレー上
に色を再現する場合は、ここでは記載しないが、人間の
視覚を含め、微妙な問題があり、ここでは、この問題を
避けるため、ＸＹＺとＸ′Ｙ′Ｚ′の直線関係を適用す
る。マトリックス表示法で、 Ｘ＝Ｗ（Ｘ′−Ｂ） 式 １１．０ Ｘ′＝（Ｘ／Ｗ）＋Ｂ 式 １１．１ ここで、Ｗ及びＢはスカラー値である。スカラー値を使
って、赤、緑、青色の成分が同じ用にして測られる。Ｗ
及びＢが明暗両色のケベルをそれぞれ調整する。Ｘ及び
Ｘ′は、それぞれ、ＸＹＺ及びＸ′Ｙ′Ｚ′を表す３成
分のコラム・マトリックスである。
【０１１２】Ｂ＝０％のとき、色は、普通、代表的ＣＲ
Ｔカラーディスプレー上でウオッシュ・アウトされ、黒
は余り暗くないると考えられる。Ｂを増すと色のコント
ラストが向上する。我々はＢ＝２％が最良の結果を与え
ることを見出だした。必要ならば、Ｂを負数にし、Ｗを
これより大きくすることが出来る。Ｗは、ＲＧＢの計算
値がその限界、普通２５５の値、を決して越えないよう
にして計算される。適切な白色標準値（例えば、Ｍａｃ
ｂｅｔｈ ＣｏｉｏｒＣｈｅｃｋｅｒのカラー・セッ
ト）を選択し、これらの相対Ｘ′Ｙ′Ｚ′値と、ＲＧＢ
値が最大レベルにあるときのＣＲＴで得られたＸＹＺ値
とを比較することにより、式１１．０をＷの演算に使う
ことができる。（この目的のために明白な公式を引き出
すことが出来る。）これによって、フルレンジのＣＲＴ
輝度を表示している色に用いることが出来る。Ｗが絶対
色測定値ＸＹＺを相対色測定値Ｘ′Ｙ′Ｚ′に変換する
のに用いられ、従って、これは無次元ではない。輝度Ｗ
及びコントラストＢの選択又は演算が図８に符号３７０
で示されている。
【０１１３】次に、本発明は、色の測定装置から得られ
た相対色測定値（Ｘ′Ｙ′Ｚ′三刺激値）と、コンピュ
ータの色成分（ＲＧＢ値）との間の変換手段を備えてい
る。この手順が図８に符号３８０で示されている。これ
は、又、ＣＲＴモデル算法６を用いて上記作業を行い、
次に、ＲＧＢ値をＣＲＴディスプレー４上に表示する、
ＸＹＺ色彩データベース７からのＸＹＺ値の検索として
図１に示されている。
【０１１４】本発明の範囲は、以上に示し説明した特定
の実施例に限定されるものではなく、本発明の範囲は添
付した請求項及びそれに相当するものによって規定され
るものと考える。
【０１１５】
【表１】
【０１１６】
【表２】
【０１１７】
【表３】
【０１１８】
【表４】
【０１１９】
【表５】
【０１２０】
【表６】
【０１２１】
【表７】
【０１２２】
【表８】
【０１２３】
【表９】
【０１２４】
【表１０】
【０１２５】
【表１１】
【０１２６】
【表１２】
【０１２７】
【表１３】
【０１２８】
【表１４】
【図面の簡単な説明】
【図１】図１は電子ディスプレー上に基質の合成色を再
現するための基本的構成要素を示す模式図的ブロック図
である。
【図２】図２はイメージディジタイザーにより試料の色
を測定するのに用いられる手順のフローチャートであ
る。
【図３】図３は、６つの灰色シェード、３つの加法混色
の原色（赤、緑、青）、３つの減法混色の原色（黄、マ
ジェンタ、シアン）、２つの肌色、及びほとんどの色空
間に集合的に拡がる１０の雑色とからなるＭａｃｂｅｔ
ｈ ＣｏｌｏｒＣｈｅｃｋｅｒの演色評価チャートの色
度のグラフ図であり、このグラフはＣＲＴ色モニタに使
われる色度の色域を示す。
【図４】図４は、Ｍａｃｂｅｔｈの灰色シェードに対す
る、ＸＹＺの測定及び予測値を持つイメージディジタイ
ザーのＲＧＢ値と、イメージディジタイザーのモデルパ
ラメータのグラフ図である。
【図５】図５は、全てのＭａｃｂｅｔｈカラーの測定値
対予測値のグラフ図である。
【図６】図６は、基質上の着色料合成体を予測するのに
用いる手順を示すフローチャートである。
【図７】図７は、染料である着色料のＫ0 、Ｋi 、
Ｓi 、及びＫijの演算に必要な１連の４つの計算を示す
図である。
【図８】図８は、陰極線管又はその他のＲＧＢ電子イメ
ージ・ディスプレー上に色測定値を表示するのに用いら
れる手順を示すフローチャートである。
【図９】図９は、ＸＹＺの測定及び予測値とＣＲＴモデ
ルパラメータとを持つ陰極線管（ＣＲＴ）ＲＧＢ値のグ
ラフ図である。
【図１０】図１０は、全ての測定された燐光色に対す
る、陰極線管の測定値対予測値のグラフ図である。
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.5 識別記号 庁内整理番号 ＦＩ 技術表示箇所 Ｇ０６Ｆ 15/70 ３１０ 9071−5Ｌ

Claims (54)

    【特許請求の範囲】
  1. 【請求項１】 イメージ・ディジタイザーを用いて試料
    のＸＹＺ値を測定する方法で、この方法が： （ａ）色のＲＧＢ値を集めるためにイメージ・ディジタ
    イザーで目標カラーチャートを走査する； （ｂ）ＸＹＺ値を集めるために目標カラーチャートを測
    定する； （ｃ）上記目標の色の値及び上記目標の色の三刺戟値を
    正規化する； （ｄ）同じ色度を持つ目標から目標の灰色シェードを選
    択する； （ｅ）上記目標の灰色シェードのＸ三刺戟値のガンマ
    （ｇ）値をピッキングする； （ｆ）下記の式で、上記ピックされたガンマ（ｇ）値の
    最小二乗フィットを演算する； Ｄ＝ＳxxＳ−Ｓ x 2 Ｇ＝（ＳxxＳy −Ｓx Ｓxy）／Ｄ Ｏ＝（−Ｓx Ｓy ＋ＳＳxy）／Ｄ 【数１】 【数２】 （ｇ）次の式で最小二乗誤差を計算する； Ｅ2 ＝Ｓyy−２ＧＳy −２ＯＳxy＋Ｏ2 Ｓxx＋２ＧＯＳ
    x ＋Ｇ2 Ｓ （ｈ）ガンマ（ｇr ）、ゲイン（Ｇr ）及びオフセット
    （Ｏr ）の最適値を見出だすために、手順（ｅ）から
    （ｇ）までを繰り返し、Ｘ三刺戟値の誤差（Ｅ）を最小
    にする； （ｉ）ガンマ（ｇg ）、ゲイン（Ｇg ）及びオフセット
    （Ｏg ）の最適値を見出だすために、手順（ｅ）から
    （ｇ）までを繰り返し、Ｙ三刺戟値の誤差（Ｅ）を最小
    にする； （ｊ）ガンマ（ｇb ）、ゲイン（Ｇb ）及びオフセット
    （Ｏb ）の最適値を見出だすために、手順（ｅ）から
    （ｇ）までを繰り返し、Ｚ三刺戟値の誤差（Ｅ）を最小
    にする； （ｋ）次の式から、マトリックスＱ及びＱrgb を組み立
    てる； 【数３】 （ｌ）次の式から、ミクシング・マトリックスＭを組み
    立てる； Ｍ＝Ｑ（Ｑrgb Ｑrgb T ）-1 ； （ｍ）ＲＧＢ値を得るために、イメージ・ディジタイザ
    ーで試料を走査する；及び、 （ｎ）次の式で、ゲイン（Ｇ）、ガンマ（ｇ）、オフセ
    ット（Ｏ）及びミクシング・マトリックス（Ｍ）の上記
    最適値を使用する： Ｘr ＝（Ｇr ＋Ｏr ＊Ｒ）gr Ｙg ＝（Ｇg ＋Ｏg ＊Ｇ）gg Ｚb ＝（Ｇb ＋Ｏb ＊Ｂ）gb Ｘ＝ＭＸrgb 以上の手順を含み、これにより、イメージ・ディジタイ
    ザーのＲＧＢ値がＸＹＺ測定値に変換される、イメージ
    ・ディジタイザーを用いて試料のＸＹＺ値を測定する方
    法。
  2. 【請求項２】 上記試料が織物材料である、請求項１記
    載の方法。
  3. 【請求項３】 上記織物材料が敷物材料である、請求項
    ２記載の方法。
  4. 【請求項４】 上記試料が紙と着色料とである、請求項
    １記載の方法。
  5. 【請求項５】 ＸＹＺ値を集めるための目標カラーチャ
    ートを測定する上記手順が測色計の利用を含む、請求項
    １記載の方法。
  6. 【請求項６】 ＸＹＺ値を集めるための目標カラーチャ
    ートを測定する上記手順が分光測光器の利用を含む、請
    求項１記載の方法。
  7. 【請求項７】 ＸＹＺ値を集めるための目標カラーチャ
    ートを測定する上記手順が彩度計の利用を含む、請求項
    １記載の方法。
  8. 【請求項８】 イメージ・ディジタイザーを用いて試料
    のＸＹＺ値を測定する装置で、この装置が： （ａ）色のＲＧＢ値を集めるためにイメージ・ディジタ
    イザーで目標カラーチャートを走査する手段と； （ｂ）色のＸＹＺ値を集めるために目標カラーチャート
    を測定する手段と； （ｃ）上記目標カラーチャート及び上記目標の色の三刺
    戟値を正規化する手段と； （ｄ）同じ色度を持つ目標の灰色シェードを選択する手
    段と； （ｅ）上記目標の灰色シェードのＸ三刺戟値のためのガ
    ンマ（ｇ）の値をピッキングする手段と； （ｆ）上記ピックされたガンマ（ｇ）の値の最小二乗フ
    ィットを、下記の式で演算する手段と； Ｄ＝ＳxxＳ−Ｓ x 2 Ｇ＝（ＳxxＳy −Ｓx Ｓxy）／Ｄ Ｏ＝（−Ｓx Ｓy ＋ＳＳxy）／Ｄ 【数４】 （ｇ）次の式で最小二乗誤差を計算する手段と； Ｅ2 ＝Ｓyy−２ＧＳy −２ＯＳxy＋Ｏ2 Ｓxx＋２ＧＯＳ
    x ＋Ｇ2 Ｓ （ｈ）ガンマ（ｇr ）、ゲイン（Ｇr ）及びオフセット
    （Ｏr ）の最適値を見出だすために、手順（ｅ）から
    （ｇ）までを繰り返し、Ｘ三刺戟値の誤差（Ｅ）を最小
    にする手段と； （ｉ）ガンマ（ｇg ）、ゲイン（Ｇg ）及びオフセット
    （Ｏg ）の最適値を見出だすために、手順（ｅ）から
    （ｇ）までを繰り返し、Ｙ三刺戟値の誤差（Ｅ）を最小
    にする手段と； （ｊ）ガンマ（ｇb ）、ゲイン（Ｇb ）及びオフセット
    （Ｏb ）の最適値を見出だすために、手順（ｅ）から
    （ｇ）までを繰り返し、Ｚ三刺戟値の誤差（Ｅ）を最小
    にする手段と； （ｋ）次の式からマトリックスＱ及びＱrgb を組み立て
    る手段と； 【数５】 （ｌ）次の式からミクシング・マトリックスＭを組み立
    てる手段と； Ｍ＝Ｑ（Ｑrgb Ｑrgb T ）-1 ； （ｍ）ＲＧＢ値を得るために、イメージ・ディジタイザ
    ーで試料を走査する手段と；及び、 （ｎ）次の式で、ゲイン（Ｇ）、ガンマ（ｇ）、オフセ
    ット（Ｏ）及びミクシング・マトリックス（Ｍ）の上記
    最適値を使用する手段と； Ｘr ＝（Ｇr ＋Ｏr ＊Ｒ）gr Ｙg ＝（Ｇg ＋Ｏg ＊Ｇ）gg Ｚb ＝（Ｇb ＋Ｏb ＊Ｂ）gb Ｘ＝ＭＸrgb を含み、これにより、イメージ・ディジタイザーのＲＧ
    Ｂ値がＸＹＺ測定値に変換される、イメージ・ディジタ
    イザーを用いて試料のＸＹＺ値を測定する装置。
  9. 【請求項９】 上記試料が織物材料である、請求項８記
    載の装置。
  10. 【請求項１０】 上記織物材料が敷物材料である、請求
    項９記載の装置。
  11. 【請求項１１】 上記試料が紙と着色料とから成り立っ
    ている、請求項１記載の装置。
  12. 【請求項１２】 ＸＹＺ値を集めるための目標カラーチ
    ャートを測定する上記手段が測色計である、請求項８記
    載の装置。
  13. 【請求項１３】 ＸＹＺ値を集めるための目標カラーチ
    ャートを測定する上記手段が分光測光器である、請求項
    ８記載の装置。
  14. 【請求項１４】 ＸＹＺ値を集めるための目標カラーチ
    ャートを測定する上記手段が彩度計である、請求項８記
    載の装置。
  15. 【請求項１５】 着色料合成体を電子ディスプレー上に
    再現する方法で、この方法が： （ａ）非彩色状態の試料のＸＹＺ値を測定する； （ｂ）その上に単一の着色料を持つ上記試料のＸＹＺ値
    を、複数の着色料濃度で測定する； （ｃ）その上に一対の着色料を持つ上記試料のＸＹＺ値
    を、複数の着色料濃度で測定する； （ｄ）非彩色の状態の上記試料のＫo を次の式を用いて
    演算する， （Ｋ／Ｓ）o ＝Ｋo ＝（１−Ｒo ）2 ／２Ｒo ； （ｅ）次の式で、ｉで示されている各単一の着色料のＫ
    i 及びＳi を演算する； ｙ＝Ｋ／Ｓ−Ｋo ｘ1 ＝Ｃi ｘ2 ＝（Ｋ／Ｓ）Ｃi ｗ＝［Ｒ／（Ｒo −Ｒ）2 ］2 【数６】 （ｆ）次の式で、ｉ及びｊで示される各１対の着色料の
    Ｋijを演算する， ｙ＝（１＋Ｓi Ｃi ＋Ｓj Ｃj ）（Ｋ／Ｓ）−（Ｋo ＋
    Ｋi Ｃi ＋Ｋj Ｃj） ｘ＝Ｃi Ｃj ｗ＝［Ｒ／（Ｒo −Ｒ）2 ］2 【数７】 ここで、インデックスｑは異なる着色料の濃度及び反射
    率測定値に関し、文字Ｐはその合計を示す； （ｇ）算出したＫo 、Ｋi 、Ｓi 及びＫijの各値を使用
    して、次の式から色のＸＹＺ三刺戟値を演算する、この
    場合、文字Ｎは着色料の総数を示す、 【数８】 Ｒ＝（１＋Ｋ／Ｓ）−［（１＋Ｋ／Ｓ）2 −１］1/2 ；
    及び、 （ｈ）上記演算した色のＸＹＺ三刺戟値を電子デイスプ
    レー上に表示する；以上の手順を含む、着色料合成体を
    電子ディスプレー上に再現する方法。
  16. 【請求項１６】 上記試料が織物材料である、請求項１
    ５記載の方法。
  17. 【請求項１７】 上記織物材料が敷物材料である、請求
    項１６記載の方法。
  18. 【請求項１８】 上記試料が紙と着色料とから成り立っ
    ている、請求項１５記載の方法。
  19. 【請求項１９】 上記電子ディスプレーが陰極線管であ
    る、請求項１５記載の方法。
  20. 【請求項２０】 上記電子ディスプレーが液晶ディスプ
    レーである、請求項１５記載の方法。
  21. 【請求項２１】 上記電子ディスプレーがエレクトロル
    ミネッセント・ディスプレーである、請求項１５記載の
    方法。
  22. 【請求項２２】 上記電子ディスプレーがプラズマ・デ
    ィスプレーである、請求項１５記載の方法。
  23. 【請求項２３】 ＸＹＺ値を測定する上記手順がイメー
    ジ・ディジタイザの利用を含む、請求項１５記載の方
    法。
  24. 【請求項２４】 ＸＹＺ値を測定する上記手順が測色計
    の利用を含む、請求項１５記載の方法。
  25. 【請求項２５】 ＸＹＺ値を測定する上記手順がクロマ
    ・メーターの利用を含む、請求項１５記載の方法。
  26. 【請求項２６】 ＸＹＺ値を測定する上記手順が分光測
    光器の利用を含む、請求項１５記載の方法。
  27. 【請求項２７】 着色料合成体を電子ディスプレー上に
    再現する装置で、この装置が： （ａ）非彩色状態の試料の色のＸＹＺ値を測定する手段
    と； （ｂ）その上に単一の着色料を持つ上記試料のＸＹＺ値
    を、複数の着色料濃度で測定する手段と； （ｃ）その上に一対の着色料を持つ上記試料のＸＹＺ値
    を、複数の着色料濃度で測定する手段と； （ｄ）非彩色状態の上記試料のＫo を次の式を用いて演
    算する手段と、 （Ｋ／Ｓ）o ＝Ｋo ＝（１−Ｒo ）2 ／２Ｒo ； （ｅ）次の式で、ｉで示されている各単一の着色料のＫ
    i 及びＳi を演算する手段と； ｙ＝Ｋ／Ｓ−Ｋo ｘ1 ＝Ｃi ｘ2 ＝（Ｋ／Ｓ）Ｃi ｗ＝［Ｒ／（Ｒo −Ｒ）2 ］2 【数９】 （ｆ）次の式で、ｉ及びｊで示されている各１対の着色
    料のＫijを演算する手段と、 ｙ＝（１＋Ｓi Ｃi ＋Ｓj Ｃj ）（Ｋ／Ｓ）−（Ｋo ＋
    Ｋi Ｃi ＋Ｋj Ｃj） ｘ＝Ｃi Ｃj ｗ＝［Ｒ／（Ｒo −Ｒ）2 ］2 【数１０】 ここで、インデックスｑは異なる着色料の濃度及び反射
    率測定値に関し、文字Ｐはその合計を示す； （ｇ）算出したＫo 、Ｋi 、Ｓi 及びＫijの各値を使用
    して、次の式から色のＸＹＺ三刺戟値を演算する手段
    と、この場合、文字Ｎは着色料の総数を示す、 【数１１】 Ｒ＝（１＋Ｋ／Ｓ）−［（１＋Ｋ／Ｓ）2 −１］1/2 ；
    及び、 （ｈ）上記演算した色のＸＹＺ三刺戟値を電子デイスプ
    レー上に表示する手段と；を含む、着色料合成体を電子
    ディスプレー上に再現する装置。
  28. 【請求項２８】 上記試料が織物材料である、請求項２
    ７記載の装置。
  29. 【請求項２９】 上記織物材料が敷物材料である、請求
    項２８記載の装置。
  30. 【請求項３０】 上記試料が紙と着色料とである、請求
    項２７記載の装置。
  31. 【請求項３１】 上記電子ディスプレーが陰極線管であ
    る、請求項２７記載の装置。
  32. 【請求項３２】 上記電子ディスプレーが液晶ディスプ
    レーである、請求項２７記載の装置。
  33. 【請求項３３】 上記電子ディスプレーがエレクトロル
    ミネッセント・ディスプレーである、請求項２７記載の
    装置。
  34. 【請求項３４】 上記電子ディスプレーがプラズマ・デ
    ィスプレーである、請求項２７記載の装置。
  35. 【請求項３５】 ＸＹＺ値を測定する上記手段がイメー
    ジ・ディジタイザの利用を含む、請求項２７記載の装
    置。
  36. 【請求項３６】 ＸＹＺ値を測定する上記手段が測色計
    の利用を含む、請求項２７記載の装置。
  37. 【請求項３７】 ＸＹＺ値を測定する上記手段がクロマ
    ・メーターの利用を含む、請求項２７記載の装置。
  38. 【請求項３８】 ＸＹＺ値を測定する上記手段が分光測
    光器の利用を含む、請求項２７記載の装置。
  39. 【請求項３９】 試料のＸＹＺ値を、電子デイスプレー
    のためのＲＧＢ値に変換する方法で、この方法が： （ａ）電子デイスプレーのためのＲＧＢ値を選択する； （ｂ）上記ＲＧＢ値のためにＸＹＺ値を集める； （ｃ）上記ＲＧＢ値及び上記ＸＹＺ値を正規化する； （ｄ）Ｘ三刺戟値のためのガンマ（ｇ）の値をピッキン
    グする； （ｅ）上記ピックされたガンマ（ｇ）の値の最小二乗フ
    ィットを、下記の式で演算する； Ｄ＝ＳxxＳ−Ｓ x 2 Ｇ＝（ＳxxＳy −Ｓx Ｓxy）／Ｄ Ｏ＝（−Ｓx Ｓy ＋ＳＳxy）／Ｄ 【数１２】 （ｆ）次の式で最小二乗誤差を計算する； Ｅ2 ＝Ｓyy−２ＧＳy −２ＯＳxy＋Ｏ2 Ｓxx＋２ＧＯＳ
    x ＋Ｇ2 Ｓ （ｇ）ガンマ（ｇr ）、ゲイン（Ｇr ）及びオフセット
    （Ｏr ）の最適値を見出だすために、手順（ｅ）から
    （ｇ）までを繰り返し、Ｘ三刺戟値の誤差（Ｅ）を最小
    にする； （ｈ）ガンマ（ｇg ）、ゲイン（Ｇg ）及びオフセット
    （Ｏg ）の最適値を見出だすために、手順（ｅ）から
    （ｇ）までを繰り返し、Ｙ三刺戟値の誤差（Ｅ）を最小
    にする； （ｉ）ガンマ（ｇb ）、ゲイン（Ｇb ）及びオフセット
    （Ｏb ）の最適値を見出だすために、手順（ｅ）から
    （ｇ）までを繰り返し、Ｚ三刺戟値の誤差（Ｅ）を最小
    にする； （ｊ）次の式からマトリックスＱ及びＱrgb を組み立て
    る； 【数１３】 （ｋ）次の式からミクシング・マトリックスＭを組み立
    てる、 Ｍ＝Ｑ（Ｑrgb Ｑrgb T ）-1 ；及び、 （ｌ）次の式で、ゲイン（Ｇ）、ガンマ（ｇ）、オフセ
    ット（Ｏ）及びミクシング・マトリックス（Ｍ）の上記
    最適値を使用する、 Ｘr ＝（Ｇr ＋Ｏr ＊Ｒ）gr Ｙg ＝（Ｇg ＋Ｏg ＊Ｇ）gg Ｚb ＝（Ｇb ＋Ｏb ＊Ｂ）gb Ｘ＝ＭＸrgb ； 以上の手順を含み、これにより、イメージ・ディジタイ
    ザーのＲＧＢ値がＸＹＺ測定値に変換される、試料のＸ
    ＹＺ値を電子デイスプレーのためのＲＧＢ値に変換する
    方法。
  40. 【請求項４０】 上記試料が織物材料である、請求項３
    ９記載の方法。
  41. 【請求項４１】 上記織物材料が敷物材料である、請求
    項４０記載の方法。
  42. 【請求項４２】 上記試料が紙と着色料とである、請求
    項３９記載の方法。
  43. 【請求項４３】 上記電子ディスプレーが陰極線管であ
    る、請求項３９記載の方法。
  44. 【請求項４４】 上記電子ディスプレーが液晶ディスプ
    レーである、請求項３９記載の方法。
  45. 【請求項４５】 上記電子ディスプレーがエレクトロル
    ミネッセント・ディスプレーである、請求項３９記載の
    方法。
  46. 【請求項４６】 上記電子ディスプレーがプラズマ・デ
    ィスプレーである、請求項３９記載の方法。
  47. 【請求項４７】 電子ディスプレーのために、合成媒体
    の色のＸＹＺ三刺戟値をＲＧＢ値に変換する装置で、こ
    の装置が： （ａ）電子デイスプレーのためのＲＧＢ値を選択する手
    段と； （ｂ）上記ＲＧＢ値のためにＸＹＺ値を集める手段と； （ｃ）上記ＲＧＢ値及び上記ＸＹＺ値を正規化する手段
    と； （ｄ）Ｘ三刺戟値のためのガンマ（ｇ）の値をピッキン
    グする手段と； （ｅ）上記ピックされたガンマ（ｇ）の値の最小二乗フ
    ィットを、下記の式で演算する手段と； Ｄ＝ＳxxＳ−Ｓ x 2 Ｇ＝（ＳxxＳy −Ｓx Ｓxy）／Ｄ Ｏ＝（−Ｓx Ｓy ＋ＳＳxy）／Ｄ 【数１４】 【数１５】 （ｆ）次の式で最小二乗誤差を計算する手段と、 Ｅ2 ＝Ｓyy−２ＧＳy −２ＯＳxy＋Ｏ2 Ｓxx＋２ＧＯＳ
    x ＋Ｇ2 Ｓ； （ｇ）ガンマ（ｇr ）、ゲイン（Ｇr ）及びオフセット
    （Ｏr ）の最適値を見出だすために、手順（ｅ）から
    （ｇ）までを繰り返し、Ｘ三刺戟値の誤差（Ｅ）を最小
    にする手段と； （ｈ）ガンマ（ｇg ）、ゲイン（Ｇg ）及びオフセット
    （Ｏg ）の最適値を見出だすために、手順（ｅ）から
    （ｇ）までを繰り返し、Ｙ三刺戟値の誤差（Ｅ）を最小
    にする手段と； （ｉ）ガンマ（ｇb ）、ゲイン（Ｇb ）及びオフセット
    （Ｏb ）の最適値を見出だすために、手順（ｅ）から
    （ｇ）までを繰り返し、Ｚ三刺戟値の誤差（Ｅ）を最小
    にする手段と； （ｊ）次の式からマトリックスＱ及びＱrgb を組み立て
    る手段と、 【数１６】 （ｋ）次の式からミクシング・マトリックスＭを組み立
    てる手段と、 Ｍ＝Ｑ（Ｑrgb Ｑrgb T ）-1 ；及び、 （ｌ）次の式で、ゲイン（Ｇ）、ガンマ（ｇ）、オフセ
    ット（Ｏ）及びミクシング・マトリックス（Ｍ）の上記
    最適値を使用する手段と、 Ｘr ＝（Ｇr ＋Ｏr ＊Ｒ）gr Ｙg ＝（Ｇg ＋Ｏg ＊Ｇ）gg Ｚb ＝（Ｇb ＋Ｏb ＊Ｂ）gb Ｘ＝ＭＸrgb ； を含み、これにより、ＲＧＢ値がＸＹＺ測定値に変換さ
    れ、電子デイスプレーのための利用される、電子ディス
    プレーのために、合成媒体の色のＸＹＺ三刺戟値をＲＧ
    Ｂ値に変換する装置。
  48. 【請求項４８】 上記試料が織物材料である、請求項４
    ７記載の装置。
  49. 【請求項４９】 上記織物材料が敷物材料である、請求
    項４８記載の装置。
  50. 【請求項５０】 上記試料が紙と着色料とである、請求
    項４７記載の装置。
  51. 【請求項５１】 上記電子ディスプレーが陰極線管であ
    る、請求項４７記載の装置。
  52. 【請求項５２】 上記電子ディスプレーが液晶ディスプ
    レーである、請求項４７記載の装置
  53. 【請求項５３】 上記電子ディスプレーがエレクトロル
    ミネッセント・ディスプレーである、請求項４７記載の
    装置。
  54. 【請求項５４】 上記電子ディスプレーがプラズマ・デ
    ィスプレーである、請求項４７記載の装置。
JP5093685A 1992-03-27 1993-03-29 合成着色料の色を電子ディスプレー上に再現する方法及び装置 Ceased JPH06109544A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US85933992A 1992-03-27 1992-03-27
US859339 1992-03-27

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2002272681A Division JP2003202267A (ja) 1992-03-27 2002-09-19 合成着色料の色を電子ディスプレー上に再現する方法及び装置

Publications (1)

Publication Number Publication Date
JPH06109544A true JPH06109544A (ja) 1994-04-19

Family

ID=25330668

Family Applications (2)

Application Number Title Priority Date Filing Date
JP5093685A Ceased JPH06109544A (ja) 1992-03-27 1993-03-29 合成着色料の色を電子ディスプレー上に再現する方法及び装置
JP2002272681A Pending JP2003202267A (ja) 1992-03-27 2002-09-19 合成着色料の色を電子ディスプレー上に再現する方法及び装置

Family Applications After (1)

Application Number Title Priority Date Filing Date
JP2002272681A Pending JP2003202267A (ja) 1992-03-27 2002-09-19 合成着色料の色を電子ディスプレー上に再現する方法及び装置

Country Status (6)

Country Link
US (2) US5375193A (ja)
EP (2) EP0562745A3 (ja)
JP (2) JPH06109544A (ja)
AU (2) AU656057B2 (ja)
CA (1) CA2091155C (ja)
NZ (1) NZ247247A (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006510010A (ja) * 2002-12-12 2006-03-23 コーニンクレッカ　フィリップス　エレクトロニクス　エヌ　ヴィ 等色関数を概算する変換機構
KR100709687B1 (ko) * 2000-04-04 2007-04-19 제이 앤드 피 코우츠 리미티드 컬러 사양
CN109637441A (zh) * 2019-01-08 2019-04-16 武汉精立电子技术有限公司 一种基于卡尔曼滤波的模组Gamma校正方法

Families Citing this family (53)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3110746B2 (ja) * 1990-10-09 2000-11-20 キヤノン株式会社 画像処理方法
US5883632A (en) * 1992-05-04 1999-03-16 Hewlett-Packard Company Coordinating color produced by two devices--using a hue-controlled machine color space, or surface scaling
US5852444A (en) * 1992-12-07 1998-12-22 Intel Corporation Application of video to graphics weighting factor to video image YUV to RGB color code conversion
US5877754A (en) * 1993-06-16 1999-03-02 Intel Corporation Process, apparatus, and system for color conversion of image signals
US6061153A (en) * 1993-12-09 2000-05-09 Canon Kabushiki Kaisha Image processing apparatus and method
US5493518A (en) * 1994-04-14 1996-02-20 Cone Mills Corporation Method and apparatus for simulating colored material
US6002498A (en) * 1994-06-15 1999-12-14 Konica Corporation Image processing method and image forming method
EP0715754A1 (en) * 1994-07-01 1996-06-12 Apple Computer, Inc. Method and system for embedding a device profile into a document and extracting a device profile from a document in a color management system
US6147671A (en) * 1994-09-13 2000-11-14 Intel Corporation Temporally dissolved dithering
US5732205A (en) * 1994-12-30 1998-03-24 Intel Corporation Color conversion using 4.5 bit palette
US5900861A (en) * 1995-09-28 1999-05-04 Intel Corporation Table-driven color conversion using interleaved indices
US5673065A (en) * 1995-12-29 1997-09-30 Intel Corporation Color reduction and conversion using an ordinal lookup table
US6043909A (en) 1996-02-26 2000-03-28 Imagicolor Corporation System for distributing and controlling color reproduction at multiple sites
US7728845B2 (en) 1996-02-26 2010-06-01 Rah Color Technologies Llc Color calibration of color image rendering devices
US6459425B1 (en) * 1997-08-25 2002-10-01 Richard A. Holub System for automatic color calibration
US5798767A (en) * 1996-03-15 1998-08-25 Rendition, Inc. Method and apparatus for performing color space conversion using blend logic
FR2749077B1 (fr) 1996-05-23 1999-08-06 Oreal Procede et dispositif de mesure de la couleur
US5864345A (en) * 1996-05-28 1999-01-26 Intel Corporation Table-based color conversion to different RGB16 formats
US6278533B1 (en) * 1996-11-29 2001-08-21 Fuji Photo Film Co., Ltd. Method of processing image signal
US6509979B2 (en) * 1997-04-03 2003-01-21 Milliken & Company Printing method using inter-pixel blending on an absorbent substrate
US5933578A (en) * 1997-04-08 1999-08-03 Barco Graphics, N.V. Method and device for determining the color appearance of color overprints
US6232954B1 (en) * 1997-05-08 2001-05-15 Imation Corp. Arrangement for high-accuracy colorimetric characterization of display devices and method therefor
JP3525997B2 (ja) * 1997-12-01 2004-05-10 富士通株式会社 文字認識方法
US6611249B1 (en) * 1998-07-22 2003-08-26 Silicon Graphics, Inc. System and method for providing a wide aspect ratio flat panel display monitor independent white-balance adjustment and gamma correction capabilities
GB2346964A (en) * 1999-01-29 2000-08-23 Richard May Method and apparatus for colour sample identification and display
JP3853105B2 (ja) * 1999-05-24 2006-12-06 富士写真フイルム株式会社 カラーモニタのモノクロ画像表示方法およびこれに用いる画像表示装置
US6654150B1 (en) 1999-06-29 2003-11-25 Kodak Polychrome Graphics Colorimetric characterization of scanned media using spectral scanner and basis spectra models
US6650446B1 (en) * 1999-06-29 2003-11-18 Kodak Polychrome Graphics Colorimetric characterization of scanned media using spectral scanner and media models
EP1079606A3 (en) * 1999-08-11 2003-07-09 Canon Kabushiki Kaisha Color space adjustment for multiple different substrates
US7221373B2 (en) * 2000-04-04 2007-05-22 J & P Coats Limited Colour specification
US7102648B1 (en) 2000-04-11 2006-09-05 Rah Color Technologies Llc Methods and apparatus for calibrating a color display
US6772151B1 (en) * 2000-06-16 2004-08-03 Polyone Corporation Color matching system and method
US6907634B2 (en) * 2001-05-01 2005-06-21 Milliken & Company Patterning system using a limited number of process colors
US8812968B2 (en) * 2001-05-03 2014-08-19 Milliken & Company Systems and methods for displaying and manipulating images of floor covering elements
US6806856B2 (en) * 2001-08-09 2004-10-19 Microsoft Corporation Reflective displays with color filter cross-talk compensation
US6792329B2 (en) * 2001-08-22 2004-09-14 Milliken & Company Construction of colored images on absorbent substrates using a computer-aided design system
US6704610B2 (en) * 2001-08-22 2004-03-09 Milliken & Company Reproduction of colored images on absorbent substrates using color blending techniques
US7070846B2 (en) * 2002-05-03 2006-07-04 Milliken & Company Carpet constructions, systems, and methods
JP4013657B2 (ja) * 2002-06-03 2007-11-28 セイコーエプソン株式会社 画像表示装置、画像表示方法および画像表示プログラムが記録されたコンピュータ読み取り可能な記録媒体
US20060069794A1 (en) * 2003-01-03 2006-03-30 Thomson Licensing Inc. System for maintaining white uniformity in a displayed video image by predicting and compensating for display register changes
DE10335145A1 (de) * 2003-07-31 2005-02-17 Giesecke & Devrient Gmbh Verfahren für die Überprüfung und Auswahl von Druckfarben
JP4551068B2 (ja) * 2003-08-19 2010-09-22 関西ペイント株式会社 近似配合検索システム
WO2005104659A2 (en) * 2004-05-04 2005-11-10 Advanced Vision Technology (Avt) Ltd. Method and system for approximating the spectrum of a plurality of color samples
CA2492961C (en) * 2004-12-23 2010-03-23 Chris Frosztega Colour coatings blender apparatus, production of colour coatings gradients and application methods and uses therefor
KR100710302B1 (ko) * 2005-05-17 2007-04-23 엘지전자 주식회사 디스플레이 장치의 색보정 장치 및 방법
US20070059844A1 (en) * 2005-09-12 2007-03-15 Techmer Pm, Llc Non-metameric color matching system and method
US20080007574A1 (en) * 2006-07-05 2008-01-10 Fujifilm Corporation Image display method
JP5029187B2 (ja) * 2007-07-20 2012-09-19 株式会社ニコン 色座標変換装置、撮像装置、色座標変換プログラム、色座標変換方法
WO2010131499A1 (ja) * 2009-05-15 2010-11-18 シャープ株式会社 画像処理装置および画像処理方法
RU2012125614A (ru) * 2009-11-20 2013-12-27 Шарп Кабусики Кайся Устройство обработки изображения и способ обработки изображения
KR101803571B1 (ko) * 2011-06-17 2017-11-30 엘지디스플레이 주식회사 입체영상표시장치와 이의 구동방법
JP7235039B2 (ja) * 2018-03-16 2023-03-08 コニカミノルタ株式会社 光沢色の色調定量化装置、光沢色の色調測定装置および光沢色の色調定量化方法
WO2020147580A1 (zh) * 2019-01-18 2020-07-23 欧普照明股份有限公司 光源参数测量方法、装置、照明系统和终端设备

Family Cites Families (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4159892A (en) * 1977-08-30 1979-07-03 World Carpets, Inc. Method of print dyeing textile materials from a photograph and products made therefrom
US4188216A (en) * 1977-08-30 1980-02-12 World Carpets, Inc. Method of print dyeing textile materials
US4232366A (en) * 1978-10-25 1980-11-04 Digital Equipment Corporation Bus for a data processing system with overlapped sequences
NZ202895A (en) * 1982-12-22 1987-09-30 Univ Auckland Color discriminator
US4635078A (en) * 1983-04-28 1987-01-06 Canon Kabushiki Kaisha Intermediate gradient image producing method
DE3371693D1 (en) * 1983-12-14 1987-06-25 Hell Rudolf Dr Ing Gmbh Method of and circuit arrangement for the recognition of chrominances and colours
US4823108A (en) * 1984-05-02 1989-04-18 Quarterdeck Office Systems Display system and memory architecture and method for displaying images in windows on a video display
JPS6163893A (ja) * 1984-09-06 1986-04-02 株式会社日立製作所 デイスプレイ装置における擬似中間調画像の表示方法
US4654721A (en) * 1985-04-12 1987-03-31 International Business Machines Corporation System for reproducing multi-level digital images on a bi-level printer of fixed dot size
US4758897A (en) * 1985-04-30 1988-07-19 Konishiroku Photo Industry Co., Ltd. Method and apparatus for estimating halftone image from binary image
US4670780A (en) * 1985-05-28 1987-06-02 Tektronix, Inc. Method of matching hardcopy colors to video display colors in which unreachable video display colors are converted into reachable hardcopy colors in a mixture-single-white (MSW) color space
US4721943A (en) * 1985-10-11 1988-01-26 Tektronix, Inc. Digital-to-analog converter for video application
US4807182A (en) * 1986-03-12 1989-02-21 Advanced Software, Inc. Apparatus and method for comparing data groups
US4680645A (en) * 1986-08-25 1987-07-14 Hewlett-Packard Company Method for rendering gray scale images with variable dot sizes
JPH0691612B2 (ja) * 1987-03-11 1994-11-14 株式会社ヤマトヤ商会 画像の階調変換方法
US4833623A (en) * 1987-07-21 1989-05-23 Burlington Industries, Inc. Orifice plate inspection apparatus and method
US4875032A (en) * 1987-10-26 1989-10-17 Mcmanus Paul A Method and apparatus for processing colorimetric parameters of a color sample
US5150199A (en) * 1990-01-02 1992-09-22 Megatronics, Inc. Method for correlating color measuring scales

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100709687B1 (ko) * 2000-04-04 2007-04-19 제이 앤드 피 코우츠 리미티드 컬러 사양
JP2006510010A (ja) * 2002-12-12 2006-03-23 コーニンクレッカ　フィリップス　エレクトロニクス　エヌ　ヴィ 等色関数を概算する変換機構
CN109637441A (zh) * 2019-01-08 2019-04-16 武汉精立电子技术有限公司 一种基于卡尔曼滤波的模组Gamma校正方法
CN109637441B (zh) * 2019-01-08 2020-09-11 武汉精立电子技术有限公司 一种基于卡尔曼滤波的模组Gamma校正方法

Also Published As

Publication number Publication date
EP0762093A2 (en) 1997-03-12
EP0562745A2 (en) 1993-09-29
JP2003202267A (ja) 2003-07-18
AU656057B2 (en) 1995-01-19
EP0562745A3 (en) 1994-07-13
AU670012B2 (en) 1996-06-27
AU3396693A (en) 1993-10-07
US5375193A (en) 1994-12-20
CA2091155A1 (en) 1993-09-28
NZ247247A (en) 1995-07-26
AU8022494A (en) 1995-02-09
US5428720A (en) 1995-06-27
CA2091155C (en) 2004-01-06
EP0762093A3 (en) 1997-07-30

Similar Documents

Publication Publication Date Title
JPH06109544A (ja) 合成着色料の色を電子ディスプレー上に再現する方法及び装置
US5798943A (en) Apparatus and process for a digital swatchbook
JP4834130B2 (ja) Ｃｃｍによる調色サンプルの再現性の評価方法
Luo Applying colour science in colour design
US5150199A (en) Method for correlating color measuring scales
Setchell Jr Colour description and communication
JPH08297054A (ja) 色感測定装置
Luo A review of chromatic adaptation transforms
JP7436453B2 (ja) 塗色検索装置
JP2004258854A (ja) 変更色の生成及び表示方法及び装置
Rougeron et al. Color fidelity in computer graphics: A survey
US20060066097A1 (en) Method of color specification and, using the same, color column structure and color chart
Cauwerts et al. A color graphic informing on the impact of electric lighting and coated glazing in complex architectural scenes
Koirala et al. Color mixing and color separation of pigments with concentration prediction
Hirschler Electronic colour communication in the textile and apparel industry
US11900507B2 (en) Visualizing wood staining
CA2444218A1 (en) Method and apparatus for reproducing blended colorants on an electronic display
Mangine Variability in experimental color matching conditions: effects of observers, daylight simulators, and color inconstancy
Rossi Real material texture color management in CAD systems for Spatial Design
WYSZECKI 30 RECENT DEVELOPMENTS ON COLOR-DIFFERENCE EVALUATIONS
Pointer et al. Food colour appearance judged using images on a computer display
Golunov et al. Automated control and interpretation of colour reproduction
Kuo Quantification of metamerism and colour constancy
Wang et al. Color Graphics in the Service of Light-Source Visualization and Design
Shaw Evaluating the 1931 CIE color matching functions

Legal Events

Date Code Title Description
A045 Written measure of dismissal of application [lapsed due to lack of payment]

Free format text: JAPANESE INTERMEDIATE CODE: A045

Effective date: 20040427