JPH099086A - 色変換方法 - Google Patents
色変換方法Info
- Publication number
- JPH099086A JPH099086A JP7156555A JP15655595A JPH099086A JP H099086 A JPH099086 A JP H099086A JP 7156555 A JP7156555 A JP 7156555A JP 15655595 A JP15655595 A JP 15655595A JP H099086 A JPH099086 A JP H099086A
- Authority
- JP
- Japan
- Prior art keywords
- color
- signal
- virtual
- color signal
- cmy
- 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.)
- Granted
Links
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/603—Colour correction or control controlled by characteristics of the picture signal generator or the picture reproducer
- H04N1/6033—Colour correction or control controlled by characteristics of the picture signal generator or the picture reproducer using test pattern analysis
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
- Color Image Communication Systems (AREA)
- Color, Gradation (AREA)
- Complex Calculations (AREA)
Abstract
(57)【要約】 （修正有）
【目的】色信号を第１表色系から第２表色系に高精度に
変換することのできる色変換方法を提供する。 【構成】段階的に色の異なる複数のカラーパッチを測色
することで、色信号ＣＭＹと刺激値信号ＸＹＺとの関係
を求めるＳ１１、関係に最小自乗法を適用して第１仮想
色信号ＣＭＹ１* と第１仮想刺激値信号ＸＹＺ１* との
関係を求めＳ１２、これらの関係を順変換テーブルとし
て設定しＳ１３、次いで、ニュートン法を用いて順変換
テーブルから第１逆変換テーブルを生成しＳ１４、さら
に、第１逆変換テーブルにおける第１仮想色信号ＣＭＹ
１* と第１仮想刺激値信号ＸＹＺ１* との関係を物理的
に整合性のある第２仮想色信号ＣＭＹ２* と第２仮想刺
激値信号ＸＹＺ２* との関係に置き換えることによりＳ
１５、第１逆変換テーブルから第２逆変換テーブルを生
成するＳ１６。第２逆変換テーブルをマッピングして得
られた逆変換テーブルを用いてＳ１７、刺激値信号ＸＹ
Ｚを色信号ＣＭＹに変換する。
変換することのできる色変換方法を提供する。 【構成】段階的に色の異なる複数のカラーパッチを測色
することで、色信号ＣＭＹと刺激値信号ＸＹＺとの関係
を求めるＳ１１、関係に最小自乗法を適用して第１仮想
色信号ＣＭＹ１* と第１仮想刺激値信号ＸＹＺ１* との
関係を求めＳ１２、これらの関係を順変換テーブルとし
て設定しＳ１３、次いで、ニュートン法を用いて順変換
テーブルから第１逆変換テーブルを生成しＳ１４、さら
に、第１逆変換テーブルにおける第１仮想色信号ＣＭＹ
１* と第１仮想刺激値信号ＸＹＺ１* との関係を物理的
に整合性のある第２仮想色信号ＣＭＹ２* と第２仮想刺
激値信号ＸＹＺ２* との関係に置き換えることによりＳ
１５、第１逆変換テーブルから第２逆変換テーブルを生
成するＳ１６。第２逆変換テーブルをマッピングして得
られた逆変換テーブルを用いてＳ１７、刺激値信号ＸＹ
Ｚを色信号ＣＭＹに変換する。
Description
【０００１】
【産業上の利用分野】本発明は、色信号を第１表色系か
ら第２表色系に変換する色変換方法に関する。
ら第２表色系に変換する色変換方法に関する。
【０００２】
【従来の技術】カラープリンタ等の出力装置を用いてカ
ラー画像の記録を行う場合、例えば、Ｃ、Ｍ、Ｙに係る
色信号を操作することで所望の色を有したカラー画像を
実現することが可能である。しかしながら、前記色信号
は、出力装置の出力特性に依存するため、特性の異なる
外部装置から供給された色信号は、前記出力特性を考慮
して色変換処理を行う必要がある。
ラー画像の記録を行う場合、例えば、Ｃ、Ｍ、Ｙに係る
色信号を操作することで所望の色を有したカラー画像を
実現することが可能である。しかしながら、前記色信号
は、出力装置の出力特性に依存するため、特性の異なる
外部装置から供給された色信号は、前記出力特性を考慮
して色変換処理を行う必要がある。
【０００３】そこで、当該出力装置を用いて色の異なる
複数の既知のカラーパッチを作成し、前記カラーパッチ
を測色することにより、例えば、前記カラーパッチの既
知の色信号ＣＭＹを出力装置に依存しない刺激値信号Ｘ
ＹＺに変換する変換関係（以下、この変換関係を「順変
換関係」という）を得、次いで、前記順変換関係から、
刺激値信号ＸＹＺを色信号ＣＭＹに変換する変換関係
（以下、この変換関係を「逆変換関係」という）を求
め、この逆変換関係を用いて前記色変換処理を行う方法
がある。
複数の既知のカラーパッチを作成し、前記カラーパッチ
を測色することにより、例えば、前記カラーパッチの既
知の色信号ＣＭＹを出力装置に依存しない刺激値信号Ｘ
ＹＺに変換する変換関係（以下、この変換関係を「順変
換関係」という）を得、次いで、前記順変換関係から、
刺激値信号ＸＹＺを色信号ＣＭＹに変換する変換関係
（以下、この変換関係を「逆変換関係」という）を求
め、この逆変換関係を用いて前記色変換処理を行う方法
がある。
【０００４】ここで、図９は、規則的な格子間隔で作成
したカラーパッチにおける既知の色信号ＣＭＹと刺激値
信号ＸＹＺとの順変換関係および逆変換関係を便宜的に
２次元表示したもので、ＣＭＹ表色系での各格子点の色
信号ＣＭＹが、ＸＹＺ表色系での各格子点の刺激値信号
ＸＹＺに対応している。この場合、ＸＹＺ表色系では、
非常に不規則な格子間隔となるため、例えば、Ａ点の刺
激値信号ＸＹＺに対する色信号ＣＭＹを求める際、Ａ
１、Ａ２、Ａ３、Ａ４点の色信号ＣＭＹを用いて補間計
算を行うべきところ、Ａ点に近いＡ３、Ａ４、Ａ５、Ａ
６点の色信号ＣＭＹを用いて補間計算を行ってしまうお
それがある。このように、逆変換を行う場合、任意の刺
激値信号ＸＹＺを囲む格子点を求めることは、極めて困
難である。
したカラーパッチにおける既知の色信号ＣＭＹと刺激値
信号ＸＹＺとの順変換関係および逆変換関係を便宜的に
２次元表示したもので、ＣＭＹ表色系での各格子点の色
信号ＣＭＹが、ＸＹＺ表色系での各格子点の刺激値信号
ＸＹＺに対応している。この場合、ＸＹＺ表色系では、
非常に不規則な格子間隔となるため、例えば、Ａ点の刺
激値信号ＸＹＺに対する色信号ＣＭＹを求める際、Ａ
１、Ａ２、Ａ３、Ａ４点の色信号ＣＭＹを用いて補間計
算を行うべきところ、Ａ点に近いＡ３、Ａ４、Ａ５、Ａ
６点の色信号ＣＭＹを用いて補間計算を行ってしまうお
それがある。このように、逆変換を行う場合、任意の刺
激値信号ＸＹＺを囲む格子点を求めることは、極めて困
難である。
【０００５】ところで、前記刺激値信号ＸＹＺから色信
号ＣＭＹを求める方法として、例えば、４点の刺激値信
号ＸＹＺを頂点とする四面体を設定し、この四面体によ
り刺激値信号ＸＹＺの空間を分割するとともに、色信号
ＣＭＹの空間も同様にして四面体で分割し、対応する四
面体中の任意の刺激値信号ＸＹＺに対する色信号ＣＭＹ
を線形演算によって求める方法がある（特開平３−１８
１７９号公報参照）。また、ニュートン法（Newton's m
ethod ）を用いて、任意の刺激値信号ＸＹＺに対応する
色信号ＣＭＹを繰り返し演算によって求める方法がある
（PHOTOGRARHICSCIENCE AND ENGINEERING・ Volume 16,
Number 2. March-April 1972 pp136-143 "Metameric
color matching in subtractive color photography"
参照）。
号ＣＭＹを求める方法として、例えば、４点の刺激値信
号ＸＹＺを頂点とする四面体を設定し、この四面体によ
り刺激値信号ＸＹＺの空間を分割するとともに、色信号
ＣＭＹの空間も同様にして四面体で分割し、対応する四
面体中の任意の刺激値信号ＸＹＺに対する色信号ＣＭＹ
を線形演算によって求める方法がある（特開平３−１８
１７９号公報参照）。また、ニュートン法（Newton's m
ethod ）を用いて、任意の刺激値信号ＸＹＺに対応する
色信号ＣＭＹを繰り返し演算によって求める方法がある
（PHOTOGRARHICSCIENCE AND ENGINEERING・ Volume 16,
Number 2. March-April 1972 pp136-143 "Metameric
color matching in subtractive color photography"
参照）。
【０００６】
【発明が解決しようとする課題】しかしながら、前者の
方法では、刺激値信号ＸＹＺの空間が歪んでいればいる
ほど、四面体を構成する４点の刺激値信号ＸＹＺの選択
が不適切となるおそれが増大するため、対応する正確な
色信号ＣＭＹを探索することが困難となってしまう。
方法では、刺激値信号ＸＹＺの空間が歪んでいればいる
ほど、四面体を構成する４点の刺激値信号ＸＹＺの選択
が不適切となるおそれが増大するため、対応する正確な
色信号ＣＭＹを探索することが困難となってしまう。
【０００７】また、後者の方法では、カラーパッチから
得られる順変換関係の範囲が制限されているため、繰り
返し演算が演算途中において不可能となるおそれがあ
る。なお、この対策として、演算途中で算出される演算
値が必ず前記範囲内となるようにダンピング係数を導入
することが考えられるが、この方法では、前記ダンピン
グ係数を試行錯誤して設定しなければならないだけでな
く、前記範囲の境界近傍で演算が収束しなくなるおそれ
もある。また、前記範囲外に仮想的な刺激値信号ＸＹＺ
および色信号ＣＭＹを設定することで演算の続行を可能
とすることが考えられるが、その設定の仕方によっては
正確な色信号ＣＭＹが得られなくなるおそれもある。
得られる順変換関係の範囲が制限されているため、繰り
返し演算が演算途中において不可能となるおそれがあ
る。なお、この対策として、演算途中で算出される演算
値が必ず前記範囲内となるようにダンピング係数を導入
することが考えられるが、この方法では、前記ダンピン
グ係数を試行錯誤して設定しなければならないだけでな
く、前記範囲の境界近傍で演算が収束しなくなるおそれ
もある。また、前記範囲外に仮想的な刺激値信号ＸＹＺ
および色信号ＣＭＹを設定することで演算の続行を可能
とすることが考えられるが、その設定の仕方によっては
正確な色信号ＣＭＹが得られなくなるおそれもある。
【０００８】本発明は、前記の不具合を解消するもので
あって、第２表色系から第１表色系への順変換関係を用
いて、色信号を第１表色系から第２表色系に高精度に変
換することのできる色変換方法を提供することを目的と
する。
あって、第２表色系から第１表色系への順変換関係を用
いて、色信号を第１表色系から第２表色系に高精度に変
換することのできる色変換方法を提供することを目的と
する。
【０００９】
【課題を解決するための手段】前記の課題を解決するた
めに、本発明は、色信号を第１表色系から第２表色系に
変換する色変換方法において、前記第２表色系の既知の
実色信号から得られる前記第１表色系の実色信号の関係
を第１順変換関係として求める第１ステップと、前記第
１順変換関係を単調関数で近似して、前記実色信号から
なる領域外に仮想色信号を設定する第２ステップと、前
記第２表色系における前記実色信号および前記仮想色信
号より構成される色信号から得られる前記第１表色系の
当該色信号の関係を第２順変換関係として求める第３ス
テップと、前記第２順変換関係から、繰り返し演算法を
用いて、前記第１表色系の色信号から得られる前記第２
表色系の色信号の関係を逆変換関係として求める第４ス
テップと、からなり、前記逆変換関係を用いて色信号を
第１表色系から第２表色系に変換することを特徴とす
る。
めに、本発明は、色信号を第１表色系から第２表色系に
変換する色変換方法において、前記第２表色系の既知の
実色信号から得られる前記第１表色系の実色信号の関係
を第１順変換関係として求める第１ステップと、前記第
１順変換関係を単調関数で近似して、前記実色信号から
なる領域外に仮想色信号を設定する第２ステップと、前
記第２表色系における前記実色信号および前記仮想色信
号より構成される色信号から得られる前記第１表色系の
当該色信号の関係を第２順変換関係として求める第３ス
テップと、前記第２順変換関係から、繰り返し演算法を
用いて、前記第１表色系の色信号から得られる前記第２
表色系の色信号の関係を逆変換関係として求める第４ス
テップと、からなり、前記逆変換関係を用いて色信号を
第１表色系から第２表色系に変換することを特徴とす
る。
【００１０】
【作用】本発明の色変換方法では、第２表色系の既知の
実色信号（例えば、色信号ＣＭＹ）に対応する第１表色
系の実色信号（例えば、刺激値信号ＸＹＺ）を求めた
後、これらの実色信号間における第１順変換関係を単調
関数で近似し、前記実色信号で構成される領域外に仮想
色信号を設定する。そして、前記実色信号および前記仮
想色信号からなる第２表色系と第１表色系との間の第２
順変換関係より、ニュートン法に代表される繰り返し演
算により前記第１表色系を前記第２表色系に変換する逆
変換関係を求め、この逆変換関係を用いて色変換処理を
行う。この場合、前記順変換関係は実色信号で構成され
る領域外の仮想色信号の関係を含んでいるため、前記繰
り返し演算での繰り返し計算中に解が色再現域外に出る
ことがあっても、当該演算を収束させて演算値を確実に
求めることができ、これによって、高精度な色変換処理
が可能となる。
実色信号（例えば、色信号ＣＭＹ）に対応する第１表色
系の実色信号（例えば、刺激値信号ＸＹＺ）を求めた
後、これらの実色信号間における第１順変換関係を単調
関数で近似し、前記実色信号で構成される領域外に仮想
色信号を設定する。そして、前記実色信号および前記仮
想色信号からなる第２表色系と第１表色系との間の第２
順変換関係より、ニュートン法に代表される繰り返し演
算により前記第１表色系を前記第２表色系に変換する逆
変換関係を求め、この逆変換関係を用いて色変換処理を
行う。この場合、前記順変換関係は実色信号で構成され
る領域外の仮想色信号の関係を含んでいるため、前記繰
り返し演算での繰り返し計算中に解が色再現域外に出る
ことがあっても、当該演算を収束させて演算値を確実に
求めることができ、これによって、高精度な色変換処理
が可能となる。
【００１１】
【実施例】図１は、本実施例の色変換方法によって作成
された逆変換テーブル（逆変換関係）が適用される画像
出力装置の構成ブロックを示す。この装置は、Ｃ（シア
ン）、Ｍ（マゼンタ）、Ｙ（イエロー）の組み合わせか
らなるカラー画像を記録する装置であり、カラー原稿を
読み取って色信号を生成し、あるいは、外部装置から供
給される色信号を受信する画像入力部１０と、前記画像
入力部１０からの色信号を当該画像出力装置の出力特性
に依存しない共通色空間の色信号である刺激値信号ＸＹ
Ｚに変換する共通色空間変換部１２と、前記刺激値信号
ＸＹＺに対して所望の画像処理を施す画像処理部１４
と、画像処理された前記刺激値信号ＸＹＺを、後述する
変換テーブル作成装置において求めた逆変換テーブルを
用いて当該画像出力装置の固有の色信号ＣＭＹに変換す
る色変換部１６と、前記色信号ＣＭＹに基づく画像を出
力する画像出力部１８とを備える。
された逆変換テーブル（逆変換関係）が適用される画像
出力装置の構成ブロックを示す。この装置は、Ｃ（シア
ン）、Ｍ（マゼンタ）、Ｙ（イエロー）の組み合わせか
らなるカラー画像を記録する装置であり、カラー原稿を
読み取って色信号を生成し、あるいは、外部装置から供
給される色信号を受信する画像入力部１０と、前記画像
入力部１０からの色信号を当該画像出力装置の出力特性
に依存しない共通色空間の色信号である刺激値信号ＸＹ
Ｚに変換する共通色空間変換部１２と、前記刺激値信号
ＸＹＺに対して所望の画像処理を施す画像処理部１４
と、画像処理された前記刺激値信号ＸＹＺを、後述する
変換テーブル作成装置において求めた逆変換テーブルを
用いて当該画像出力装置の固有の色信号ＣＭＹに変換す
る色変換部１６と、前記色信号ＣＭＹに基づく画像を出
力する画像出力部１８とを備える。
【００１２】なお、本実施例では、色信号ＣＭＹと刺激
値信号ＸＹＺとの間での変換処理について説明するが、
他の表色系、例えば、Ｌａｂ表色系、Ｌｕｖ表色系、Ｒ
ＧＢ表色系等の間での変換処理にも同様にして適用する
ことができることは勿論である。
値信号ＸＹＺとの間での変換処理について説明するが、
他の表色系、例えば、Ｌａｂ表色系、Ｌｕｖ表色系、Ｒ
ＧＢ表色系等の間での変換処理にも同様にして適用する
ことができることは勿論である。
【００１３】図２は、変換テーブル作成装置の構成ブロ
ックを示す。この装置は、測色計２０、刺激値信号算出
部２２、色信号供給部２４、第１仮想色信号生成部２
６、順変換テーブル生成部２８、第１逆変換テーブル生
成部３０、第２仮想色信号生成部３２、第２逆変換テー
ブル生成部３４およびマッピング部３６を備える。
ックを示す。この装置は、測色計２０、刺激値信号算出
部２２、色信号供給部２４、第１仮想色信号生成部２
６、順変換テーブル生成部２８、第１逆変換テーブル生
成部３０、第２仮想色信号生成部３２、第２逆変換テー
ブル生成部３４およびマッピング部３６を備える。
【００１４】測色計２０は、図１に示す画像出力装置を
用いて作成したカラーパッチの反射光の分布を測定す
る。刺激値信号算出部２２は、前記反射光の分布と与え
られた観察光源の分光分布から刺激値信号ＸＹＺを算出
する。色信号供給部２４は、前記カラーパッチを生成す
るための色信号ＣＭＹを第１仮想色信号生成部２６およ
び順変換テーブル生成部２８に供給する。第１仮想色信
号生成部２６は、前記色信号供給部２４から供給される
色信号ＣＭＹと、前記刺激値信号算出部２２からの刺激
値信号ＸＹＺとに基づき、最小自乗法を用いて、前記色
信号ＣＭＹおよび刺激値信号ＸＹＺにより構成される領
域外の十分広い領域において、第１仮想色信号である第
１仮想色信号ＣＭＹ１* およびそれに対応する第１仮想
刺激値信号ＸＹＺ１* を設定する。順変換テーブル生成
部２８は、色信号供給部２４からの色信号ＣＭＹおよび
刺激値信号算出部２２からの刺激値信号ＸＹＺ（以下、
これらを「実色信号」という）と、前記第１仮想色信号
生成部２６からの前記第１仮想色信号ＣＭＹ１* とに基
づき、前記第１仮想色信号ＣＭＹ１* を含む色信号ＣＭ
Ｙから前記第１仮想刺激値信号ＸＹＺ１* を含む刺激値
信号ＸＹＺを得る順変換テーブルを生成する。第１逆変
換テーブル生成部３０は、前記順変換テーブルより、繰
り返し演算法であるニュートン法を用いて、第１仮想刺
激値信号ＸＹＺ１* を含む刺激値信号ＸＹＺから第１仮
想色信号ＣＭＹ１* を含む色信号ＣＭＹを得る第１逆変
換テーブルを生成する。第２仮想色信号生成部３２は、
前記第１逆変換テーブルから第１仮想色信号ＣＭＹ１*
および第１仮想刺激値信号ＸＹＺ１* の関係を除去した
第１逆変換テーブルを用いて、前記色信号ＣＭＹおよび
刺激値信号ＸＹＺにより構成される領域外の十分広い領
域において、第２仮想色信号ＣＭＹ２* およびそれに対
応する第２仮想刺激値信号ＸＹＺ２* を設定する。第２
逆変換テーブル生成部３４は、前記第２仮想刺激値信号
ＸＹＺ２* を含む刺激値信号ＸＹＺから前記第２仮想色
信号ＣＭＹ２* を含む色信号ＣＭＹを得るための第２逆
変換テーブルを生成する。マッピング部３６は、前記第
２逆変換テーブルを図１に示す画像出力装置で表現可能
な固有の色空間の範囲の逆変換テーブルにマッピングす
る。
用いて作成したカラーパッチの反射光の分布を測定す
る。刺激値信号算出部２２は、前記反射光の分布と与え
られた観察光源の分光分布から刺激値信号ＸＹＺを算出
する。色信号供給部２４は、前記カラーパッチを生成す
るための色信号ＣＭＹを第１仮想色信号生成部２６およ
び順変換テーブル生成部２８に供給する。第１仮想色信
号生成部２６は、前記色信号供給部２４から供給される
色信号ＣＭＹと、前記刺激値信号算出部２２からの刺激
値信号ＸＹＺとに基づき、最小自乗法を用いて、前記色
信号ＣＭＹおよび刺激値信号ＸＹＺにより構成される領
域外の十分広い領域において、第１仮想色信号である第
１仮想色信号ＣＭＹ１* およびそれに対応する第１仮想
刺激値信号ＸＹＺ１* を設定する。順変換テーブル生成
部２８は、色信号供給部２４からの色信号ＣＭＹおよび
刺激値信号算出部２２からの刺激値信号ＸＹＺ（以下、
これらを「実色信号」という）と、前記第１仮想色信号
生成部２６からの前記第１仮想色信号ＣＭＹ１* とに基
づき、前記第１仮想色信号ＣＭＹ１* を含む色信号ＣＭ
Ｙから前記第１仮想刺激値信号ＸＹＺ１* を含む刺激値
信号ＸＹＺを得る順変換テーブルを生成する。第１逆変
換テーブル生成部３０は、前記順変換テーブルより、繰
り返し演算法であるニュートン法を用いて、第１仮想刺
激値信号ＸＹＺ１* を含む刺激値信号ＸＹＺから第１仮
想色信号ＣＭＹ１* を含む色信号ＣＭＹを得る第１逆変
換テーブルを生成する。第２仮想色信号生成部３２は、
前記第１逆変換テーブルから第１仮想色信号ＣＭＹ１*
および第１仮想刺激値信号ＸＹＺ１* の関係を除去した
第１逆変換テーブルを用いて、前記色信号ＣＭＹおよび
刺激値信号ＸＹＺにより構成される領域外の十分広い領
域において、第２仮想色信号ＣＭＹ２* およびそれに対
応する第２仮想刺激値信号ＸＹＺ２* を設定する。第２
逆変換テーブル生成部３４は、前記第２仮想刺激値信号
ＸＹＺ２* を含む刺激値信号ＸＹＺから前記第２仮想色
信号ＣＭＹ２* を含む色信号ＣＭＹを得るための第２逆
変換テーブルを生成する。マッピング部３６は、前記第
２逆変換テーブルを図１に示す画像出力装置で表現可能
な固有の色空間の範囲の逆変換テーブルにマッピングす
る。
【００１５】本実施例の色変換方法が適用される画像出
力装置および変換テーブル作成装置は、基本的には以上
のように構成される。次に、図３に示すフローチャート
に従って、逆変換テーブルを作成する方法について説明
する。
力装置および変換テーブル作成装置は、基本的には以上
のように構成される。次に、図３に示すフローチャート
に従って、逆変換テーブルを作成する方法について説明
する。
【００１６】先ず、画像出力装置を構成する画像出力部
１８に対して、規則的な間隔からなる複数の色信号ＣＭ
Ｙを与えることにより、Ｃ、Ｍ、Ｙの濃度および混合割
合が段階的に変化するカラーパッチを作成する（ステッ
プＳ１０）。次いで、図２に示す変換テーブル作成装置
を用いて前記各カラーパッチの測色を行う（ステップＳ
１１）。
１８に対して、規則的な間隔からなる複数の色信号ＣＭ
Ｙを与えることにより、Ｃ、Ｍ、Ｙの濃度および混合割
合が段階的に変化するカラーパッチを作成する（ステッ
プＳ１０）。次いで、図２に示す変換テーブル作成装置
を用いて前記各カラーパッチの測色を行う（ステップＳ
１１）。
【００１７】すなわち、各カラーパッチの分光分布が測
色計２０により読み取られ、刺激値信号算出部２２にお
いて、前記分光分布と与えられた観察光源の分光分布か
ら刺激値信号ＸＹＺが算出される。なお、前記刺激値信
号ＸＹＺは、観察光源の分光分布をＳ（λ）（λは光の
波長）、物体の分光反射率分布をρ（λ）、人間の目に
対する等色関数をｘ（λ）、ｙ（λ）、ｚ（λ）、可視
光の範囲を３８０ｎｍ〜７８０ｎｍとして、
色計２０により読み取られ、刺激値信号算出部２２にお
いて、前記分光分布と与えられた観察光源の分光分布か
ら刺激値信号ＸＹＺが算出される。なお、前記刺激値信
号ＸＹＺは、観察光源の分光分布をＳ（λ）（λは光の
波長）、物体の分光反射率分布をρ（λ）、人間の目に
対する等色関数をｘ（λ）、ｙ（λ）、ｚ（λ）、可視
光の範囲を３８０ｎｍ〜７８０ｎｍとして、
【００１８】
【数１】
【００１９】の関係から導き出すことができる。ここ
で、基準化係数Ｋは、
で、基準化係数Ｋは、
【００２０】
【数２】
【００２１】で与えられる。
【００２２】この場合、Ｃ、Ｍ、Ｙの各色のカラーパッ
チの最小濃度および最大濃度を０および２５５に規格化
して８段階に分割して出力すると、図４に示す５１２
（＝８×８×８）のカラーパッチを得ることができる。
そして、各カラーパッチの色信号ＣＭＹに対応する刺激
値信号ＸＹＺの関係である第１順変換テーブルを得るこ
とができる。
チの最小濃度および最大濃度を０および２５５に規格化
して８段階に分割して出力すると、図４に示す５１２
（＝８×８×８）のカラーパッチを得ることができる。
そして、各カラーパッチの色信号ＣＭＹに対応する刺激
値信号ＸＹＺの関係である第１順変換テーブルを得るこ
とができる。
【００２３】次に、本実施例では、前記カラーパッチに
よる色再現域（図４の斜線で示す部分）外の十分離れた
位置に第１仮想色信号ＣＭＹ１* を設定し、前記第１仮
想色信号ＣＭＹ１* に対応する第１仮想刺激値信号ＸＹ
Ｚ１* を求める（ステップＳ１２）。
よる色再現域（図４の斜線で示す部分）外の十分離れた
位置に第１仮想色信号ＣＭＹ１* を設定し、前記第１仮
想色信号ＣＭＹ１* に対応する第１仮想刺激値信号ＸＹ
Ｚ１* を求める（ステップＳ１２）。
【００２４】第１仮想色信号生成部２６は、刺激値信号
算出部２２からの刺激値信号ＸＹＺと、色信号供給部２
４からの前記刺激値信号ＸＹＺに対応する色信号ＣＭＹ
との関係が単調性を持つことを仮定し、最小自乗法を用
いて前記第１仮想刺激値信号ＸＹＺ１* を求める。この
場合、「単調性を持つ」という前記の仮定は、色信号Ｃ
ＭＹが増加するに従ってカラーパッチからの反射光量が
減少し、（１）〜（３）式から得られる刺激値信号ＸＹ
Ｚが単調に減少する、という事実に基づいて保証され
る。
算出部２２からの刺激値信号ＸＹＺと、色信号供給部２
４からの前記刺激値信号ＸＹＺに対応する色信号ＣＭＹ
との関係が単調性を持つことを仮定し、最小自乗法を用
いて前記第１仮想刺激値信号ＸＹＺ１* を求める。この
場合、「単調性を持つ」という前記の仮定は、色信号Ｃ
ＭＹが増加するに従ってカラーパッチからの反射光量が
減少し、（１）〜（３）式から得られる刺激値信号ＸＹ
Ｚが単調に減少する、という事実に基づいて保証され
る。
【００２５】そこで、カラーパッチから得られる色信号
ＣＭＹおよび刺激値信号ＸＹＺの全てのデータセット
（Ｃｉ，Ｍｉ，Ｙｉ，Ｘｉ，Ｙｉ，Ｚｉ）（ｉはカラー
パッチの番号）を用いて最小自乗法によりＸ−ＣＭＹ、
Ｙ−ＣＭＹ、Ｚ−ＣＭＹの４次元平面を求める。この場
合、前記４次元平面を、 Ｔ＝Ａ・Ｄ …（５） と定義する。なお、前記（５）式は、次の（６）式の関
係を表すものとする。
ＣＭＹおよび刺激値信号ＸＹＺの全てのデータセット
（Ｃｉ，Ｍｉ，Ｙｉ，Ｘｉ，Ｙｉ，Ｚｉ）（ｉはカラー
パッチの番号）を用いて最小自乗法によりＸ−ＣＭＹ、
Ｙ−ＣＭＹ、Ｚ−ＣＭＹの４次元平面を求める。この場
合、前記４次元平面を、 Ｔ＝Ａ・Ｄ …（５） と定義する。なお、前記（５）式は、次の（６）式の関
係を表すものとする。
【００２６】
【数３】
【００２７】そして、前記（５）式の関係を満たす係数
Ａを、次の（７）式に示すＥが最小となるように、最小
自乗法を用いて求める。なお、（７）式において、ｉは
カラーパッチの番号を表し、また、上付記号のＴはマト
リクスの行と列を入れ換えた転置を表す。
Ａを、次の（７）式に示すＥが最小となるように、最小
自乗法を用いて求める。なお、（７）式において、ｉは
カラーパッチの番号を表し、また、上付記号のＴはマト
リクスの行と列を入れ換えた転置を表す。
【００２８】
【数４】
【００２９】（７）式から係数Ａが求められると、各４
次元平面Ｘ−ＣＭＹ、Ｙ−ＣＭＹ、Ｚ−ＣＭＹが決定す
る。次に、このようにして求められた４次元平面Ｘ−Ｃ
ＭＹ、Ｙ−ＣＭＹ、Ｚ−ＣＭＹ上において、カラーパッ
チの色信号ＣＭＹから十分に遠い位置にある第１仮想色
信号ＣＭＹ１* およびそれに対応した第１仮想刺激値信
号ＸＹＺ１* を求める。
次元平面Ｘ−ＣＭＹ、Ｙ−ＣＭＹ、Ｚ−ＣＭＹが決定す
る。次に、このようにして求められた４次元平面Ｘ−Ｃ
ＭＹ、Ｙ−ＣＭＹ、Ｚ−ＣＭＹ上において、カラーパッ
チの色信号ＣＭＹから十分に遠い位置にある第１仮想色
信号ＣＭＹ１* およびそれに対応した第１仮想刺激値信
号ＸＹＺ１* を求める。
【００３０】図５は、色信号ＣＭＹおよび第１仮想色信
号ＣＭＹ１* と、刺激値信号ＸＹＺおよび前記のように
して生成された第１仮想刺激値信号ＸＹＺ１* との関係
を２次元的模式図として示したものである。すなわち、
カラーパッチから得られる色信号ＣＭＹと刺激値信号Ｘ
ＹＺとの関係は、前述した知見から、例えば、点ａ１〜
ａ４に示すように、単調減少の関係にある。そして、こ
れらの点ａ１〜ａ４を用いて最小自乗法により算出され
た４次元平面を点線で示すと、前記４次元平面上の第１
仮想色信号ＣＭＹ１* に対応する第１仮想刺激値信号Ｘ
ＹＺ１* を表す点ｂ１、ｂ２と、点ａ１〜ａ４とを結ぶ
面は、実線で示すように単調減少の関係となる。従っ
て、前記のようにして、最小自乗法を用いて第１仮想色
信号ＣＭＹ１* に対する第１仮想刺激値信号ＸＹＺ１*
を生成することにより、第１仮想色信号ＣＭＹ１* を含
む色信号ＣＭＹと、第１仮想刺激値信号ＸＹＺ１* を含
む刺激値信号ＸＹＺとの間で単調関係を維持することが
できる。
号ＣＭＹ１* と、刺激値信号ＸＹＺおよび前記のように
して生成された第１仮想刺激値信号ＸＹＺ１* との関係
を２次元的模式図として示したものである。すなわち、
カラーパッチから得られる色信号ＣＭＹと刺激値信号Ｘ
ＹＺとの関係は、前述した知見から、例えば、点ａ１〜
ａ４に示すように、単調減少の関係にある。そして、こ
れらの点ａ１〜ａ４を用いて最小自乗法により算出され
た４次元平面を点線で示すと、前記４次元平面上の第１
仮想色信号ＣＭＹ１* に対応する第１仮想刺激値信号Ｘ
ＹＺ１* を表す点ｂ１、ｂ２と、点ａ１〜ａ４とを結ぶ
面は、実線で示すように単調減少の関係となる。従っ
て、前記のようにして、最小自乗法を用いて第１仮想色
信号ＣＭＹ１* に対する第１仮想刺激値信号ＸＹＺ１*
を生成することにより、第１仮想色信号ＣＭＹ１* を含
む色信号ＣＭＹと、第１仮想刺激値信号ＸＹＺ１* を含
む刺激値信号ＸＹＺとの間で単調関係を維持することが
できる。
【００３１】以上のようにして第１仮想色信号ＣＭＹ１
* に対する第１仮想刺激値信号ＸＹＺ１* を求めた後、
順変換テーブル生成部２８において、第１仮想色信号Ｃ
ＭＹ１* を含む色信号ＣＭＹを、第１仮想刺激値信号Ｘ
ＹＺ１* を含む刺激値信号ＸＹＺに変換するための第２
順変換テーブルを求める（ステップＳ１３）。
* に対する第１仮想刺激値信号ＸＹＺ１* を求めた後、
順変換テーブル生成部２８において、第１仮想色信号Ｃ
ＭＹ１* を含む色信号ＣＭＹを、第１仮想刺激値信号Ｘ
ＹＺ１* を含む刺激値信号ＸＹＺに変換するための第２
順変換テーブルを求める（ステップＳ１３）。
【００３２】次に、第１逆変換テーブル生成部３０にお
いて、前記第２順変換テーブルを用いて、ニュートン法
により、第１仮想刺激値信号ＸＹＺ１* を含む刺激値信
号ＸＹＺを、第１仮想色信号ＣＭＹ１* を含む色信号Ｃ
ＭＹに変換するための第１逆変換テーブルを求める（ス
テップＳ１４）。
いて、前記第２順変換テーブルを用いて、ニュートン法
により、第１仮想刺激値信号ＸＹＺ１* を含む刺激値信
号ＸＹＺを、第１仮想色信号ＣＭＹ１* を含む色信号Ｃ
ＭＹに変換するための第１逆変換テーブルを求める（ス
テップＳ１４）。
【００３３】この場合、ＣＭＹ色空間は、図４に示すよ
うに、格子状に配列された第１仮想色信号ＣＭＹ１* を
含む色信号ＣＭＹによって構成されているため、前記第
２順変換テーブルを用いて、任意の色信号ＣＭＹに対す
る刺激値信号ＸＹＺを補間によって求めることができ
る。一方、前記ＣＭＹ色空間から得られるＸＹＺ色空間
は、図９に示すように歪んでいるため、任意の刺激値信
号ＸＹＺに対する色信号ＣＭＹを求めることは困難であ
る。そのため、前記ＸＹＺ色空間に、図９の点線で示す
ように規則的な間隔で配列した刺激値信号ＸＹＺを設定
し、前記刺激値信号ＸＹＺに対する色信号ＣＭＹをニュ
ートン法を用いて算出することにより第１逆変換テーブ
ルを求める。
うに、格子状に配列された第１仮想色信号ＣＭＹ１* を
含む色信号ＣＭＹによって構成されているため、前記第
２順変換テーブルを用いて、任意の色信号ＣＭＹに対す
る刺激値信号ＸＹＺを補間によって求めることができ
る。一方、前記ＣＭＹ色空間から得られるＸＹＺ色空間
は、図９に示すように歪んでいるため、任意の刺激値信
号ＸＹＺに対する色信号ＣＭＹを求めることは困難であ
る。そのため、前記ＸＹＺ色空間に、図９の点線で示す
ように規則的な間隔で配列した刺激値信号ＸＹＺを設定
し、前記刺激値信号ＸＹＺに対する色信号ＣＭＹをニュ
ートン法を用いて算出することにより第１逆変換テーブ
ルを求める。
【００３４】図６は、前記第１逆変換テーブルを求める
処理のサブルーチンである。そこで、ＸＹＺ色空間の格
子上における目標値を（Ｘ０，Ｙ０，Ｚ０）とし、繰り
返し演算での許容誤差をΔＥmin に設定する（ステップ
Ｓ２０）。次いで、ＣＭＹ色空間での既知の初期値（Ｃ
１，Ｍ１，Ｙ１）を設定し（ステップＳ２１）、前記第
２順変換テーブルを用いて、前記初期値（Ｃ１，Ｍ１，
Ｙ１）に対する刺激値（Ｘ１，Ｙ１，Ｚ１）を求める
（ステップＳ２２）。そして、目標値（Ｘ０，Ｙ０，Ｚ
０）と前記刺激値（Ｘ１，Ｙ１，Ｚ１）との誤差量ΔＥ
を求め（ステップＳ２３）、前記誤差量ΔＥと許容誤差
ΔＥmin とを比較する（ステップＳ２４）。なお、前記
許容誤差ΔＥmin は、例えば、Ｌａｂ表色系における色
差式に従って算出された色差として設定することができ
る。この場合、｜ΔＥ｜＜ΔＥminでなければ、修正値
（ΔＣ，ΔＭ，ΔＹ）を算出し（ステップＳ２５）、初
期値（Ｃ１，Ｍ１，Ｙ１）を前記修正値（ΔＣ，ΔＭ，
ΔＹ）だけ修正した後（ステップＳ２６）、ステップＳ
２２〜Ｓ２４の処理を繰り返す。
処理のサブルーチンである。そこで、ＸＹＺ色空間の格
子上における目標値を（Ｘ０，Ｙ０，Ｚ０）とし、繰り
返し演算での許容誤差をΔＥmin に設定する（ステップ
Ｓ２０）。次いで、ＣＭＹ色空間での既知の初期値（Ｃ
１，Ｍ１，Ｙ１）を設定し（ステップＳ２１）、前記第
２順変換テーブルを用いて、前記初期値（Ｃ１，Ｍ１，
Ｙ１）に対する刺激値（Ｘ１，Ｙ１，Ｚ１）を求める
（ステップＳ２２）。そして、目標値（Ｘ０，Ｙ０，Ｚ
０）と前記刺激値（Ｘ１，Ｙ１，Ｚ１）との誤差量ΔＥ
を求め（ステップＳ２３）、前記誤差量ΔＥと許容誤差
ΔＥmin とを比較する（ステップＳ２４）。なお、前記
許容誤差ΔＥmin は、例えば、Ｌａｂ表色系における色
差式に従って算出された色差として設定することができ
る。この場合、｜ΔＥ｜＜ΔＥminでなければ、修正値
（ΔＣ，ΔＭ，ΔＹ）を算出し（ステップＳ２５）、初
期値（Ｃ１，Ｍ１，Ｙ１）を前記修正値（ΔＣ，ΔＭ，
ΔＹ）だけ修正した後（ステップＳ２６）、ステップＳ
２２〜Ｓ２４の処理を繰り返す。
【００３５】ここで、前記修正値（ΔＣ，ΔＭ，ΔＹ）
は、次のようにして求められる。すなわち、任意の色信
号ＣＭＹが与えられたとき、この色信号ＣＭＹ（点ａで
示す）に対する刺激値信号ＸＹＺは、図７に示すよう
に、８つの格子点ｃ０〜ｃ７における色信号（ｃ０、ｍ
０、ｙ０）〜（ｃ７、ｍ７、ｙ７）に対応する刺激値信
号（ｘ０、ｙ０、ｚ０）〜（ｘ７、ｙ７、ｚ７）、格子
点ｃ０〜ｃ７で囲まれる直方体の体積Ｖ、前記直方体内
の任意の補間点ｃにより８分割された体積Ｖ０〜Ｖ７を
用いて、
は、次のようにして求められる。すなわち、任意の色信
号ＣＭＹが与えられたとき、この色信号ＣＭＹ（点ａで
示す）に対する刺激値信号ＸＹＺは、図７に示すよう
に、８つの格子点ｃ０〜ｃ７における色信号（ｃ０、ｍ
０、ｙ０）〜（ｃ７、ｍ７、ｙ７）に対応する刺激値信
号（ｘ０、ｙ０、ｚ０）〜（ｘ７、ｙ７、ｚ７）、格子
点ｃ０〜ｃ７で囲まれる直方体の体積Ｖ、前記直方体内
の任意の補間点ｃにより８分割された体積Ｖ０〜Ｖ７を
用いて、
【００３６】
【数５】
【００３７】として求めることができる。この場合、
（８）〜（１０）式の関係において、色信号ＣＭＹに対
する刺激値信号ＸＹＺが微小範囲内で線形であると仮定
すると、前記色信号ＣＭＹの微小変化量である修正値
（ΔＣ，ΔＭ，ΔＹ）と前記刺激値信号ＸＹＺの微小変
化量（ΔＸ，ΔＹ，ΔＺ）とは、ニュートン法から、
（８）〜（１０）式の関係において、色信号ＣＭＹに対
する刺激値信号ＸＹＺが微小範囲内で線形であると仮定
すると、前記色信号ＣＭＹの微小変化量である修正値
（ΔＣ，ΔＭ，ΔＹ）と前記刺激値信号ＸＹＺの微小変
化量（ΔＸ，ΔＹ，ΔＺ）とは、ニュートン法から、
【００３８】
【数６】
【００３９】の関係を満たすことになる。なお、Ｊはヤ
コビアン行列であり、また、Ｄｃ、Ｄｍ、ＤｙはＣ、
Ｍ、Ｙを表し、ＸＹＺのＹとＣＭＹのＹとの混同を避け
るために表記を変えたものである。（１１）式におい
て、ヤコビアン行列Ｊが求まれば、色信号ＣＭＹの修正
値（ΔＣ，ΔＭ，ΔＹ）に対する刺激値信号ＸＹＺの微
小変化量（ΔＸ，ΔＹ，ΔＺ）を予測することができ
る。この場合、前記ヤコビアン行列Ｊは、（８）〜（１
０）式を色信号ＣＭＹで偏微分することで求められる。
従って、色信号ＣＭＹの修正値（ΔＣ，ΔＭ，ΔＹ）
は、
コビアン行列であり、また、Ｄｃ、Ｄｍ、ＤｙはＣ、
Ｍ、Ｙを表し、ＸＹＺのＹとＣＭＹのＹとの混同を避け
るために表記を変えたものである。（１１）式におい
て、ヤコビアン行列Ｊが求まれば、色信号ＣＭＹの修正
値（ΔＣ，ΔＭ，ΔＹ）に対する刺激値信号ＸＹＺの微
小変化量（ΔＸ，ΔＹ，ΔＺ）を予測することができ
る。この場合、前記ヤコビアン行列Ｊは、（８）〜（１
０）式を色信号ＣＭＹで偏微分することで求められる。
従って、色信号ＣＭＹの修正値（ΔＣ，ΔＭ，ΔＹ）
は、
【００４０】
【数７】
【００４１】として求められる。
【００４２】以上のようにして得られるヤコビアン行列
Ｊを用いて繰り返し計算を行うことにより、任意の目標
値（Ｘ０，Ｙ０，Ｚ０）に対する色信号ＣＭＹを求める
ことができる。同様の処理をＸＹＺ色空間における格子
上の全目標値に対して行い、刺激値信号ＸＹＺを色信号
ＣＭＹに変換するための第１逆変換テーブルが生成され
る（ステップＳ１４）。
Ｊを用いて繰り返し計算を行うことにより、任意の目標
値（Ｘ０，Ｙ０，Ｚ０）に対する色信号ＣＭＹを求める
ことができる。同様の処理をＸＹＺ色空間における格子
上の全目標値に対して行い、刺激値信号ＸＹＺを色信号
ＣＭＹに変換するための第１逆変換テーブルが生成され
る（ステップＳ１４）。
【００４３】ここで、ニュートン法では、収束の条件と
して、求める解を有する方程式が単調関数であることが
必要である。この場合、第１仮想色信号ＣＭＹ１* を含
む色信号ＣＭＹと第１仮想刺激値信号ＸＹＺ１* を含む
刺激値信号ＸＹＺとの間には、単調関係が成立し、しか
も、ＣＭＹ色空間およびＸＹＺ色空間が仮想色信号の部
分にまで拡大されており、繰り返し演算中において算出
される全ての色信号ＣＭＹおよび刺激値信号ＸＹＺが存
在するため、発散することなく確実に繰り返し演算を実
行し、正確な値を得ることができる。
して、求める解を有する方程式が単調関数であることが
必要である。この場合、第１仮想色信号ＣＭＹ１* を含
む色信号ＣＭＹと第１仮想刺激値信号ＸＹＺ１* を含む
刺激値信号ＸＹＺとの間には、単調関係が成立し、しか
も、ＣＭＹ色空間およびＸＹＺ色空間が仮想色信号の部
分にまで拡大されており、繰り返し演算中において算出
される全ての色信号ＣＭＹおよび刺激値信号ＸＹＺが存
在するため、発散することなく確実に繰り返し演算を実
行し、正確な値を得ることができる。
【００４４】ところで、前記のようにして生成された第
１逆変換テーブルには、第１仮想色信号ＣＭＹ１* およ
び第１仮想刺激値信号ＸＹＺ１* から導かれた関係が含
まれている。これらの信号は、ニュートン法による処理
に際して確実に解が得られるようにするために、図５に
示すように、最小自乗法を用いて設定されたものであ
り、実際の物理特性とは異なっている。実際の物理特性
を加味した仮想色信号は、図５に一点鎖線で示すよう
に、点ａ１〜ａ４の関係に対して滑らかに連結するよう
な点ｂ１* 、ｂ２* として設定すべきである。しかしな
がら、図８に示すように、前記点ｂ１* 、ｂ２* を、例
えば、Ｘ−ＣＭ色空間で表した場合、単調関係を維持し
ていないことがあるため、ニュートン法による処理には
適用できない。
１逆変換テーブルには、第１仮想色信号ＣＭＹ１* およ
び第１仮想刺激値信号ＸＹＺ１* から導かれた関係が含
まれている。これらの信号は、ニュートン法による処理
に際して確実に解が得られるようにするために、図５に
示すように、最小自乗法を用いて設定されたものであ
り、実際の物理特性とは異なっている。実際の物理特性
を加味した仮想色信号は、図５に一点鎖線で示すよう
に、点ａ１〜ａ４の関係に対して滑らかに連結するよう
な点ｂ１* 、ｂ２* として設定すべきである。しかしな
がら、図８に示すように、前記点ｂ１* 、ｂ２* を、例
えば、Ｘ−ＣＭ色空間で表した場合、単調関係を維持し
ていないことがあるため、ニュートン法による処理には
適用できない。
【００４５】そこで、点ｂ１、ｂ２のデータを用いてニ
ュートン法を実行した後、第２仮想色信号生成部３２に
おいて、前記一点鎖線で示す関係を得るため、例えば、
点ｂ１* に近い点ａ１、ａ２のみを用いて、最小自乗法
により点ｂ１* における第２仮想色信号ＣＭＹ２* およ
び第２仮想刺激値信号ＸＹＺ２* を求め、また、点ｂ２
* に近い点ａ３、ａ４のみを用いて、最小自乗法により
点ｂ１* における第２仮想色信号ＣＭＹ２* および第２
仮想刺激値信号ＸＹＺ２* を求める（ステップＳ１
５）。そして、第２逆変換テーブル生成部３４におい
て、前記第２仮想色信号ＣＭＹ２* および第２仮想刺激
値信号ＸＹＺ２* を第１仮想色信号ＣＭＹ１*および第
１仮想刺激値信号ＸＹＺ１* と置き換え、これらの第２
仮想色信号ＣＭＹ２* および第２仮想刺激値信号ＸＹＺ
２* を含む色信号ＣＭＹおよび刺激値信号ＸＹＺの関係
を第２逆変換テーブルとして設定する（ステップＳ１
６）。
ュートン法を実行した後、第２仮想色信号生成部３２に
おいて、前記一点鎖線で示す関係を得るため、例えば、
点ｂ１* に近い点ａ１、ａ２のみを用いて、最小自乗法
により点ｂ１* における第２仮想色信号ＣＭＹ２* およ
び第２仮想刺激値信号ＸＹＺ２* を求め、また、点ｂ２
* に近い点ａ３、ａ４のみを用いて、最小自乗法により
点ｂ１* における第２仮想色信号ＣＭＹ２* および第２
仮想刺激値信号ＸＹＺ２* を求める（ステップＳ１
５）。そして、第２逆変換テーブル生成部３４におい
て、前記第２仮想色信号ＣＭＹ２* および第２仮想刺激
値信号ＸＹＺ２* を第１仮想色信号ＣＭＹ１*および第
１仮想刺激値信号ＸＹＺ１* と置き換え、これらの第２
仮想色信号ＣＭＹ２* および第２仮想刺激値信号ＸＹＺ
２* を含む色信号ＣＭＹおよび刺激値信号ＸＹＺの関係
を第２逆変換テーブルとして設定する（ステップＳ１
６）。
【００４６】最後に、前記のようにして設定された第２
逆変換テーブルを当該画像出力装置における色信号ＣＭ
Ｙの再現可能範囲にマッピング（Gamut Mapping ）する
ことで、刺激値信号ＸＹＺから色信号ＣＭＹに変換する
ことのできる所望の逆変換テーブルを得ることができる
（ステップＳ１７）。
逆変換テーブルを当該画像出力装置における色信号ＣＭ
Ｙの再現可能範囲にマッピング（Gamut Mapping ）する
ことで、刺激値信号ＸＹＺから色信号ＣＭＹに変換する
ことのできる所望の逆変換テーブルを得ることができる
（ステップＳ１７）。
【００４７】以上のようにして生成された逆変換テーブ
ルは、画像出力装置の色変換部１６に設定される。そこ
で、当該画像出力装置において、画像入力部１０で読み
込まれた画像信号は、共通色空間変換部１２において刺
激値信号ＸＹＺに変換された後、画像処理部１４で所望
の画像処理が施される。次いで、色変換部１６におい
て、体積補間を用いて、前記刺激値信号ＸＹＺが前記逆
変換テーブルに従って色信号ＣＭＹに変換されると同時
に、マッピング処理が施される。そして、前記色信号Ｃ
ＭＹは、画像出力部１８より可視画像として出力され
る。
ルは、画像出力装置の色変換部１６に設定される。そこ
で、当該画像出力装置において、画像入力部１０で読み
込まれた画像信号は、共通色空間変換部１２において刺
激値信号ＸＹＺに変換された後、画像処理部１４で所望
の画像処理が施される。次いで、色変換部１６におい
て、体積補間を用いて、前記刺激値信号ＸＹＺが前記逆
変換テーブルに従って色信号ＣＭＹに変換されると同時
に、マッピング処理が施される。そして、前記色信号Ｃ
ＭＹは、画像出力部１８より可視画像として出力され
る。
【００４８】なお、前記逆変換テーブルには、当該画像
出力装置において再現可能なＣＭＹ色空間よりも広範囲
の色信号ＣＭＹに対する刺激値信号ＸＹＺが第２仮想刺
激値信号ＸＹＺ２* として設定されている。従って、例
えば、前記刺激値信号ＸＹＺが当該画像出力装置におい
て再現可能な領域の境界近傍にある場合であっても、前
記第２仮想刺激値信号ＸＹＺ２* を用いて体積補間によ
り高精度に色信号ＣＭＹを求めることができる。
出力装置において再現可能なＣＭＹ色空間よりも広範囲
の色信号ＣＭＹに対する刺激値信号ＸＹＺが第２仮想刺
激値信号ＸＹＺ２* として設定されている。従って、例
えば、前記刺激値信号ＸＹＺが当該画像出力装置におい
て再現可能な領域の境界近傍にある場合であっても、前
記第２仮想刺激値信号ＸＹＺ２* を用いて体積補間によ
り高精度に色信号ＣＭＹを求めることができる。
【００４９】ここで、前記逆変換テーブルでは、色信号
ＣＭＹへの変換処理とマッピング処理とを同時に行うよ
うにしているため、所望の色信号ＣＭＹを高速に得るこ
とができる。これに対して、マッピング処理を行った
後、刺激値信号ＸＹＺから色信号ＣＭＹへの変換処理を
行うようにすると、色信号ＣＭＹに対応する刺激値信号
ＸＹＺの精度を維持することができるため、変換処理の
精度をより向上させることができる。また、刺激値信号
ＸＹＺから色信号ＣＭＹに変換を行った後、マッピング
処理を行うようにすると、色再現可能な領域の境界近傍
で変換処理の精度が上記実施例の場合よりも向上すると
ともに、マッピング処理が容易となる効果がある。
ＣＭＹへの変換処理とマッピング処理とを同時に行うよ
うにしているため、所望の色信号ＣＭＹを高速に得るこ
とができる。これに対して、マッピング処理を行った
後、刺激値信号ＸＹＺから色信号ＣＭＹへの変換処理を
行うようにすると、色信号ＣＭＹに対応する刺激値信号
ＸＹＺの精度を維持することができるため、変換処理の
精度をより向上させることができる。また、刺激値信号
ＸＹＺから色信号ＣＭＹに変換を行った後、マッピング
処理を行うようにすると、色再現可能な領域の境界近傍
で変換処理の精度が上記実施例の場合よりも向上すると
ともに、マッピング処理が容易となる効果がある。
【００５０】なお、上述した実施例では、単調関数によ
る近似法として最小自乗法を用いているが、他の近似法
として重回帰分析の手法や重み付き最小自乗法を用いる
ことができる。また、繰り返し演算法として、ニュート
ン法を用いる代わりに、最急降下法、準ニュートン法等
を用いることができる。
る近似法として最小自乗法を用いているが、他の近似法
として重回帰分析の手法や重み付き最小自乗法を用いる
ことができる。また、繰り返し演算法として、ニュート
ン法を用いる代わりに、最急降下法、準ニュートン法等
を用いることができる。
【００５１】
【発明の効果】以上のように、本発明では、第１表色系
と第２表色系との関係を単調関数で近似することによ
り、前記第２表色系の既知の実色信号で構成される領域
外に仮想色信号を設定し、前記仮想色信号を含む順変換
関係を求め、前記順変換関係から繰り返し演算法を用い
て、前記第１表色系を前記第２表色系に変換する逆変換
関係を求めている。この場合、前記逆変換関係は実色信
号で構成される領域外の仮想色信号の関係を含んでいる
ため、繰り返し演算において演算値を確実に求めること
ができ、これによって、高精度な色変換処理が可能とな
る。
と第２表色系との関係を単調関数で近似することによ
り、前記第２表色系の既知の実色信号で構成される領域
外に仮想色信号を設定し、前記仮想色信号を含む順変換
関係を求め、前記順変換関係から繰り返し演算法を用い
て、前記第１表色系を前記第２表色系に変換する逆変換
関係を求めている。この場合、前記逆変換関係は実色信
号で構成される領域外の仮想色信号の関係を含んでいる
ため、繰り返し演算において演算値を確実に求めること
ができ、これによって、高精度な色変換処理が可能とな
る。
【図１】画像出力装置の構成ブロック図である。
【図２】変換テーブル作成装置の構成ブロック図であ
る。
る。
【図３】変換テーブルの作成処理のフローチャートであ
る。
る。
【図４】仮想色信号を含む色信号ＣＭＹの説明図であ
る。
る。
【図５】最小自乗法による仮想色信号生成の説明図であ
る。
る。
【図６】ニュートン法による色信号算出処理のフローチ
ャートである。
ャートである。
【図７】体積補間の説明図である。
【図８】Ｘ−ＣＭ色空間で表した色信号および仮想色信
号の説明図である。
号の説明図である。
【図９】ＣＭＹ色空間とＸＹＺ色空間の変換関係の説明
図である。
図である。
１０…画像入力部 １２…共通色空
間変換部 １４…画像処理部 １６…色変換部 １８…画像出力部 ２０…測色計 ２２…刺激値信号算出部 ２４…色信号供
給部 ２６…第１仮想色信号生成部 ２８…順変換テ
ーブル生成部 ３０…第１逆変換テーブル生成部 ３２…第２仮想
色信号生成部 ３４…第２逆変換テーブル生成部 ３６…マッピン
グ部
間変換部 １４…画像処理部 １６…色変換部 １８…画像出力部 ２０…測色計 ２２…刺激値信号算出部 ２４…色信号供
給部 ２６…第１仮想色信号生成部 ２８…順変換テ
ーブル生成部 ３０…第１逆変換テーブル生成部 ３２…第２仮想
色信号生成部 ３４…第２逆変換テーブル生成部 ３６…マッピン
グ部
Claims (3)
- 【請求項１】色信号を第１表色系から第２表色系に変換
する色変換方法において、 前記第２表色系の既知の実色信号から得られる前記第１
表色系の実色信号の関係を第１順変換関係として求める
第１ステップと、 前記第１順変換関係を単調関数で近似して、前記実色信
号からなる領域外に仮想色信号を設定する第２ステップ
と、 前記第２表色系における前記実色信号および前記仮想色
信号より構成される色信号から得られる前記第１表色系
の当該色信号の関係を第２順変換関係として求める第３
ステップと、 前記第２順変換関係から、繰り返し演算法を用いて、前
記第１表色系の色信号から得られる前記第２表色系の色
信号の関係を逆変換関係として求める第４ステップと、 からなり、前記逆変換関係を用いて色信号を第１表色系
から第２表色系に変換することを特徴とする色変換方
法。 - 【請求項２】請求項１記載の方法において、 第２ステップにおける前記仮想色信号は、最小自乗法を
用いて設定されることを特徴とする色変換方法。 - 【請求項３】請求項１記載の方法において、 第４ステップにおける前記繰り返し演算法は、ニュート
ン法であることを特徴とする色変換方法。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP15655595A JP3714701B2 (ja) | 1995-06-22 | 1995-06-22 | 色変換方法 |
| US08/655,235 US5619427A (en) | 1995-06-22 | 1996-06-05 | Color conversion method and color conversion table generating apparatus |
| DE1996633262 DE69633262T2 (de) | 1995-06-22 | 1996-06-18 | Farbumwandlungsverfahren und Vorrichtung zum Herstellen von Farbumwandlungstabellen |
| EP96109789A EP0750419B1 (en) | 1995-06-22 | 1996-06-18 | Color conversion method and color conversion table generating apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP15655595A JP3714701B2 (ja) | 1995-06-22 | 1995-06-22 | 色変換方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH099086A true JPH099086A (ja) | 1997-01-10 |
| JP3714701B2 JP3714701B2 (ja) | 2005-11-09 |
Family
ID=15630362
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP15655595A Expired - Lifetime JP3714701B2 (ja) | 1995-06-22 | 1995-06-22 | 色変換方法 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US5619427A (ja) |
| EP (1) | EP0750419B1 (ja) |
| JP (1) | JP3714701B2 (ja) |
| DE (1) | DE69633262T2 (ja) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6229916B1 (en) | 1997-09-30 | 2001-05-08 | Fuji Photo Film Co., Ltd. | Color transformation look-up table |
| US6377366B1 (en) | 1996-11-15 | 2002-04-23 | Fuji Photo Film Co., Ltd. | Color reproduction data converting method |
| JP2009207098A (ja) * | 2008-02-29 | 2009-09-10 | Seiko Epson Corp | 印刷制御装置、印刷システムおよび印刷制御プログラム |
| WO2012108304A1 (en) * | 2011-02-09 | 2012-08-16 | Ricoh Company, Ltd. | Control device, image forming apparatus, and control method |
| US8428347B2 (en) | 2011-02-10 | 2013-04-23 | Fuji Xerox Co., Ltd. | Color processing apparatus, color processing method, and non-transitory computer readable medium |
Families Citing this family (31)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU695071B2 (en) * | 1995-02-14 | 1998-08-06 | Canon Information Systems Research Australia Pty Ltd | Colour conversion method |
| US5872898A (en) * | 1995-09-15 | 1999-02-16 | Agfa Gevaert N.V. | Method and apparatus for calculating color gamuts |
| US5791781A (en) * | 1995-11-01 | 1998-08-11 | Samsung Electronics Co., Ltd. | Method of determining color temperature for color display device |
| US5650942A (en) * | 1996-02-02 | 1997-07-22 | Light Source Computer Images, Inc. | Appearance-based technique for rendering colors on an output device |
| US6229915B1 (en) * | 1996-08-22 | 2001-05-08 | Fuji Photo Film Co., Ltd. | Method of converting data of a first colorimetric system to data of a second colorimetric system, including a mapping of out of boundary data |
| US6134029A (en) * | 1996-09-09 | 2000-10-17 | Light Source Acquisition Company | Scanner calibration technique |
| US6281984B1 (en) * | 1997-03-25 | 2001-08-28 | International Business Machines Corporation | Enhanced system, method and program for converting an externally defined four dimensional colorant (CMYK) into an equivalent four dimensional colorant defined in terms of the four inks (C′M′Y′K′) that are associated with a given printer |
| US6005968A (en) * | 1997-08-29 | 1999-12-21 | X-Rite, Incorporated | Scanner calibration and correction techniques using scaled lightness values |
| JPH11313219A (ja) * | 1998-01-20 | 1999-11-09 | Fujitsu Ltd | カラーデータ変換方法 |
| JP2000013619A (ja) * | 1998-06-19 | 2000-01-14 | Fuji Photo Film Co Ltd | 色変換方法 |
| EP1001611A1 (en) * | 1998-11-13 | 2000-05-17 | Seiko Epson Corporation | Improving color calibration accuracy by reiterative mapping |
| DE19930688A1 (de) * | 1999-07-02 | 2001-01-04 | Byk Gardner Gmbh | Vorrichtung und Verfahren zur Bestimmung der Qualität von Oberflächen |
| JP2001045316A (ja) * | 1999-08-04 | 2001-02-16 | Fuji Photo Film Co Ltd | 色修正関係抽出方法および色修正方法 |
| DE10056723A1 (de) | 1999-12-15 | 2001-06-21 | Heidelberger Druckmasch Ag | Verfahren zur Erzeugung eines an einen realen Prozeß angepaßten zweiten Bilddatensatzes aus einem ersten Bilddatensatz und Colour-Management-Verfahren |
| PL200158B1 (pl) * | 2000-12-06 | 2008-12-31 | Avt Advanced Vision Technology | Sposób sterowania kolorem na arkuszu drukarskim |
| US7136523B2 (en) * | 2001-10-04 | 2006-11-14 | Canon Kabushiki Kaisha | Color correction table forming method and apparatus, control program and storage medium |
| US7304766B2 (en) * | 2002-04-01 | 2007-12-04 | Hewlett-Packard Development Company, L.P. | Method and apparatus for data adjustment |
| JP2004128585A (ja) * | 2002-09-30 | 2004-04-22 | Canon Inc | 色補正テーブル作成方法及び制御プログラム及び記憶媒体及び装置 |
| US7375851B2 (en) | 2003-01-21 | 2008-05-20 | Xerox Corporation | State-space based modeling of pixel elements of a dynamically varying color marking device |
| US7218358B2 (en) * | 2004-06-15 | 2007-05-15 | Coretronic Corporation | Method and apparatus for calibrating color temperature of color display devices |
| US7626742B2 (en) * | 2005-02-16 | 2009-12-01 | Samsung Electronics Co., Ltd. | Color data conversion apparatus and method using gamut boundary descriptors |
| US7505173B2 (en) * | 2005-06-15 | 2009-03-17 | Xerox Corporation | System and method for spatial gray balance calibration using hybrid sensing systems |
| US8699103B2 (en) * | 2005-06-15 | 2014-04-15 | Xerox Corporation | System and method for dynamically generated uniform color objects |
| US7800779B2 (en) * | 2005-12-21 | 2010-09-21 | Xerox Corporation | System and method for image based control using inline sensors |
| JP4652980B2 (ja) * | 2006-01-06 | 2011-03-16 | キヤノン株式会社 | 色処理装置、色処理方法、及びプログラム |
| DE102006040609B4 (de) * | 2006-08-30 | 2009-07-02 | OCé PRINTING SYSTEMS GMBH | Verfahren und Anordnung zum Ermitteln einer Farbzuordnung von Quellfarben und Zielfarben |
| JP5230118B2 (ja) * | 2007-03-30 | 2013-07-10 | キヤノン株式会社 | 階調性評価装置および階調性評価方法 |
| US7945093B2 (en) | 2008-01-16 | 2011-05-17 | Eastman Kodak Company | Simplified color workflow |
| EP2091227A1 (en) * | 2008-02-15 | 2009-08-19 | Thomson Licensing | Method for adjusting the settings of a reproduction color device |
| JP5693192B2 (ja) * | 2010-12-13 | 2015-04-01 | キヤノン株式会社 | 色処理装置、色処理方法、プログラム及び記憶媒体 |
| EP2672718A1 (en) | 2012-06-07 | 2013-12-11 | Thomson Licensing | Color calibration of an image capture device in a way that is adaptive to the scene to be captured |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4929978A (en) * | 1987-10-23 | 1990-05-29 | Matsushita Electric Industrial Co., Ltd. | Color correction method for color copier utilizing correction table derived from printed color samples |
| JP2872690B2 (ja) * | 1989-06-15 | 1999-03-17 | 三菱電機株式会社 | 色変換方法 |
| EP0454275B1 (en) * | 1990-04-26 | 1996-07-03 | Bayer Corporation | Color separation method |
-
1995
- 1995-06-22 JP JP15655595A patent/JP3714701B2/ja not_active Expired - Lifetime
-
1996
- 1996-06-05 US US08/655,235 patent/US5619427A/en not_active Expired - Lifetime
- 1996-06-18 DE DE1996633262 patent/DE69633262T2/de not_active Expired - Lifetime
- 1996-06-18 EP EP96109789A patent/EP0750419B1/en not_active Expired - Lifetime
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6377366B1 (en) | 1996-11-15 | 2002-04-23 | Fuji Photo Film Co., Ltd. | Color reproduction data converting method |
| US6229916B1 (en) | 1997-09-30 | 2001-05-08 | Fuji Photo Film Co., Ltd. | Color transformation look-up table |
| JP2009207098A (ja) * | 2008-02-29 | 2009-09-10 | Seiko Epson Corp | 印刷制御装置、印刷システムおよび印刷制御プログラム |
| WO2012108304A1 (en) * | 2011-02-09 | 2012-08-16 | Ricoh Company, Ltd. | Control device, image forming apparatus, and control method |
| CN103370927A (zh) * | 2011-02-09 | 2013-10-23 | 株式会社理光 | 控制装置、图像形成设备和控制方法 |
| US8947742B2 (en) | 2011-02-09 | 2015-02-03 | Ricoh Company, Ltd. | Image forming apparatus for increasing combined-color reproduction accuracy for superimposed colors |
| CN103370927B (zh) * | 2011-02-09 | 2015-10-14 | 株式会社理光 | 控制装置、图像形成设备和控制方法 |
| US8428347B2 (en) | 2011-02-10 | 2013-04-23 | Fuji Xerox Co., Ltd. | Color processing apparatus, color processing method, and non-transitory computer readable medium |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3714701B2 (ja) | 2005-11-09 |
| EP0750419A2 (en) | 1996-12-27 |
| DE69633262T2 (de) | 2005-01-13 |
| US5619427A (en) | 1997-04-08 |
| EP0750419B1 (en) | 2004-09-01 |
| DE69633262D1 (de) | 2004-10-07 |
| EP0750419A3 (en) | 1997-07-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH099086A (ja) | 色変換方法 | |
| KR100902700B1 (ko) | 색 처리 방법 및 색 처리 장치 | |
| JPH1141475A (ja) | 色変換データ作成方法およびそれを用いた色変換装置 | |
| JP3720691B2 (ja) | 色処理方法および装置 | |
| JPH09508768A (ja) | カラー・スキャナ信号を測色値に変換する方法および装置 | |
| JP2002112055A (ja) | 色再現特性表示装置および色再現特性表示プログラム記憶媒体 | |
| JP6156401B2 (ja) | 色変換方法、プログラム、及び画像処理装置 | |
| WO1998052347A1 (en) | Arrangement for efficient characterization of printing devices and method therefor | |
| JPH07170421A (ja) | 色変換方法 | |
| JP7297547B2 (ja) | 画像処理装置、画像形成装置 | |
| JP2003125210A (ja) | 色処理方法、記憶媒体、色処理装置、色変換装置、および画像形成装置 | |
| JP4189188B2 (ja) | 色変換定義補正方法、色変換定義補正装置、および色変換定義補正プログラム | |
| JP4077938B2 (ja) | カラーパッチに関する測色値の推定方法、該方法を利用したデバイス・プロファイルの作成方法および画像処理装置 | |
| JP2000253270A (ja) | 色変換テーブル作成装置、作成方法、色変換テーブル作成プログラムを記録した記録媒体、および、色変換装置 | |
| JP2003259141A (ja) | 色変換定義作成方法、色変換定義作成装置、および色変換定義作成プログラム | |
| JP2003289446A (ja) | 色再現特性作成方法、色再現特性作成装置、および色再現特性作成プログラム | |
| JPH08116460A (ja) | カラーマスキングパラメータ決定装置 | |
| JP2008312119A (ja) | カラーパッチに関する測色値の推定方法、該方法を利用したデバイス・プロファイルの作成方法および画像処理装置 | |
| CN111800557A (zh) | 颜色转换信息生成方法、以及颜色转换信息生成装置 | |
| US20080144114A1 (en) | Method and system for dynamic printer profiling | |
| JP3910323B2 (ja) | プロファイル作成方法およびプロファイル作成装置 | |
| JP4459127B2 (ja) | 色処理方法およびその装置 | |
| JP2001245172A (ja) | 階調補正曲線評価結果表示方法、階調補正曲線評価結果表示装置、および階調補正曲線評価結果表示プログラム記憶媒体 | |
| JP2004242072A (ja) | 色変換定義修正装置および色変換定義修正プログラム | |
| JP2003018419A (ja) | 画像処理方法及び画像処理装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050201 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050408 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050606 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20050816 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050823 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080902 Year of fee payment: 3 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080902 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090902 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090902 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100902 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100902 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110902 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120902 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130902 Year of fee payment: 8 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| EXPY | Cancellation because of completion of term |