(57)【特許請求の範囲】
【請求項１】
３次元の第１の色空間における色域から４次元の第２の色空間における色域への変換を行うカラー画像処理方法において、前記第１の色空間の色座標を表す３つの独立変数のうち少なくとも明度に対して重みを付けるとともに該重みの大きさを色相角に応じて変化させ、該重みを付けた変数を含む３つの独立変数により前記第１の色空間における色域の外郭から前記第２の色空間における色域を前記第１の色空間に変換したときの色域の外郭への距離が最小になる関係を求め、該関係に基づいて前記第１の色空間から前記第２の色空間への色変換を行うことを特徴とするカラー画像処理方法。
【請求項２】
前記第１の色空間は、均等知覚色空間であることを特徴とする請求項１に記載のカラー画像処理方法。
【請求項３】
前記第１の色空間から前記第２の色空間への色変換は、色相角により圧縮または伸長または不変を含むことを特徴とする請求項１または請求項２に記載のカラー画像処理方法。
【請求項４】
前記第１の色空間における色域は、入力機器に依存した入力色信号の色空間を前記第１の色空間へ色変換した場合の色域であることを特徴とする請求項１ないし請求項３のいずれか１項に記載のカラー画像処理方法。
【請求項５】
前記第１の色空間における色域は、標準的な第３の色空間における入力色信号を前記第１の色空間へ色変換した場合の色域であることを特徴とする請求項１ないし請求項３のいずれか１項に記載のカラー画像処理方法。
【請求項６】
前記第２の色空間における色域は、出力装置において再現可能な色域であることを特徴とする請求項１ないし請求項５のいずれか１項に記載のカラー画像処理方法。
【請求項７】
前記関係は、前記第１の色空間における色域中及び色域の外郭上の複数の代表点について予め求めておき、前記第１の色空間の色信号に応じて予め求めておいた複数の前記関係から補間によって前記第２の色空間における色信号を決定することを特徴とする請求項１ないし請求項６のいずれか１項に記載のカラー画像処理方法。
【請求項８】
３次元の第１の色空間における色域から４次元の第２の色空間における色域への変換を行うための色変換係数を生成する色変換係数生成方法において、前記第２の色空間における色域を前記第１の色空間に変換したときの色域の外郭である変換外郭を求め、前記第１の色空間における色域の外郭の色から前記第２の色空間における変換外郭への重み付きの距離が最小になるように前記第１の色空間の色座標を表す３つの独立変数のうち複数の変数に付ける重みを算出しておき、前記第１の色空間における色信号を順次入力し、算出しておいた複数の前記重みから該色信号に対応する重みを算出し、算出した色信号に対応する重みに基づいて入力した前記色信号を前記第２の色空間における色信号に変換し、該変換結果と入力した前記色信号とを関係付けることによって順次色変換係数を生成することを特徴とする色変換係数生成方法。
【請求項９】
前記第１の色空間における色信号を順次入力したとき、入力した前記色信号が前記第２の色空間における色域内か否かを判定し、前記第２の色空間における色域外であるときに前記重みに基づいて前記変換外郭の色に変換することを特徴とする請求項８に記載の色変換係数生成方法。
【請求項１０】
前記重みは、前記第１の色空間における少なくとも明度に対するものであることを特徴とする請求項８または請求項９に記載の色変換係数生成方法。
【請求項１１】
３次元の第１の色空間における色域から４次元の第２の色空間における色域への変換を行うための色変換係数を生成する色変換係数生成装置において、前記第２の色空間における色域を前記第１の色空間に変換したときの色域の外郭である変換外郭を求める色域算出手段と、前記第１の色空間における色域の外郭の色から前記第２の色空間における変換外郭への重み付きの距離が最小になるように前記第１の色空間の色座標を表す３つの独立変数のうち複数の変数に付ける重みを算出する重み決定手段と、入力された第１の色空間における色信号に対応する重みを前記重み決定手段で算出しておいた複数の前記重みから算出する重み算出手段と、入力された前記色信号を前記重み算出手段で算出した複数の前記重みを用いて前記第２の色空間における色信号に変換し該変換結果と入力した前記色信号とを関係付けることによって順次色変換係数を生成する変換手段を有することを特徴とする色変換係数生成装置。
【請求項１２】
前記重み決定手段は、官能評価的手法を用いて前記重みを算出することを特徴とする請求項１１に記載の色変換係数生成装置。
【請求項１３】
さらに、入力された前記第１の色空間における色信号が前記第２の色空間における色域内か否かを判定する色域内外判定手段を有し、前記変換手段は、前記色域内外判定手段で色域外であると判定された色信号について前記重みを用いた変換を行うことを特徴とする請求項１１または請求項１２に記載の色変換係数生成装置。
【請求項１４】
３次元の第１の色空間における色域の色信号から４次元の第２の色空間における色域の色信号への変換を行うカラー画像処理装置において、前記第１の色空間における色域中及び色域の外郭上の複数の代表点について請求項１１ないし請求項１３のいずれか１項に記載の色変換係数生成装置によって求められた色変換係数を格納し前記第１の色空間における色域の色信号が入力されると１ないし複数の色変換係数を出力する変換テーブルと、該変換テーブルから出力される１ないし複数の色変換係数に基づいて入力された前記第１の色空間における色域の色信号に対応する第２の色空間における色域の色信号を算出する出力色信号算出手段を有することを特徴とするカラー画像処理装置。
【請求項１５】
３次元の第１の色空間における色域の色信号から４次元の第２の色空間における色域の色信号への変換を行うカラー画像処理装置において、前記第１の色空間の色座標を表す３つの独立変数のうち少なくとも明度に対して重みを付けるとともに該重みの大きさを色相角に応じて変化させ該重みを付けた変数を含む３つの独立変数により前記第１の色空間における色域の外郭から前記第２の色空間における色域を前記第１の色空間に変換したときの色域の外郭への距離が最小になる関係を代表色について求め前記第１の色空間における色域の色信号が入力されると１ないし複数の関係を出力する変換テーブルと、該変換テーブルから出力される１ないし複数の前記関係に基づいて入力された前記第１の色空間における色域の色信号に対応する第２の色空間における色域の色信号を算出する出力色信号算出手段を有することを特徴とするカラー画像処理装置。
【請求項１６】
３次元の第１の色空間における色域から４次元の第２の色空間における色域への変換を行う処理をコンピュータに実行させるプログラムを格納したコンピュータが読取可能な記憶媒体において、前記第１の色空間の色座標を表す３つの独立変数のうち少なくとも明度に対して重みを付けるとともに該重みの大きさを色相角に応じて変化させ、該重みを付けた変数を含む３つの独立変数により前記第１の色空間における色域の外郭から前記第２の色空間における色域を前記第１の色空間に変換したときの色域の外郭への距離が最小になる関係を求め、該関係に基づいて前記第１の色空間から前記第２の色空間への色変換を行う処理をコンピュータに実行させるプログラムを格納したことを特徴とするコンピュータが読取可能な記憶媒体。
【請求項１７】
３次元の第１の色空間における色域から４次元の第２の色空間における色域への変換を行うための色変換係数を生成する処理をコンピュータに実行させるプログラムを格納したコンピュータが読取可能な記憶媒体において、前記第２の色空間における色域を前記第１の色空間に変換したときの色域の外郭である変換外郭を求め、前記第１の色空間における色域の外郭の色から前記第２の色空間における変換外郭への重み付きの距離が最小になるように前記第１の色空間の色座標を表す３つの独立変数のうち少なくとも１つの変数に付ける重みを算出しておき、前記第１の色空間における色信号を順次入力し、算出しておいた複数の前記重みから該色信号に対応する重みを算出し、算出した色信号に対応する重みに基づいて入力した前記色信号を前記第２の色空間における色信号に変換し、該変換結果と入力した前記色信号とを関係付けることによって順次色変換係数を生成する処理をコンピュータに実行させるプログラムを格納したことを特徴とするコンピュータが読取可能な記憶媒体。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、３次元の第１の色空間における色域から４次元の第２の色空間における色域への変換を行うカラー画像処理方法、カラー画像処理装置、及び、変換の際の色変換係数を生成する色変換係数生成方法及び色変換係数生成装置と、それらの処理をコンピュータに実行させるプログラムを格納したコンピュータが読取可能な記憶媒体に関するものである。
【０００２】
【従来の技術】
ＤＴＰ（Ｄｅｓｋｔｏｐ  Ｐｕｂｌｉｓｈｉｎｇ）システム等の編集システムでは、例えば画像入力装置より得られた画像、あるいは、コンピュータグラフィックスなど生成された画像をＣＲＴ（Ｃａｔｈｏｄｅ  Ｒａｙ  Ｔｕｂｅ）などの表示装置上で画像の作成や編集を行ない、最終的なドキュメントをハードコピーの形で記録装置から出力することが行われている。用いられる画像入力装置や、表示装置または記録装置などの出力装置として、各種の装置が開発されている。それぞれの画像入力装置や出力装置は、その特性により、入力あるいは出力可能な色域（色再現範囲）が異なっている。そのため、例えば表示装置で編集していた色が記録装置では再現できない場合も発生する。このような場合でも記録装置では何らかの色で再現する必要があり、再現可能な色に変換して出力を行っている。
【０００３】
出力装置へ出力するために入力画像信号を出力画像信号に変換する場合、通常、入力画像信号の属する色空間から出力画像信号の属する色空間への色変換と、色域圧縮と呼ばれる入出力間の色再現範囲の相違を吸収する変換が行われる。典型的な例として、ディスプレイに表示するための標準的な色空間であるｓＲＧＢ色空間で表現された画像を、シアン（Ｃ）、マゼンタ（Ｍ）、イエロー（Ｙ）、ブラック（Ｋ）の４色の色材を用いて画像を形成するＣＭＹＫプリンタに出力する場合を考える。ＣＭＹＫプリンタでは、各色の色材に対応した４つの成分を有するＣＭＹＫ色空間の色信号を入力とする。そのため、ｓＲＧＢ色空間からＣＭＹＫ色空間への色空間の変換が必要になる。通常は、各種の装置に対応するため、ｓＲＧＢ色空間の色信号を、一旦、機器に独立な色信号、例えば１９７６ＣＩＥＬ* ａ* ｂ* 色空間（以降単にＬａｂまたはＬａｂ色空間と表記し、各成分をＬ、ａ、ｂと表記する）などの均等知覚色空間に変換した後、ＣＭＹＫ色空間に変換してＣＭＹＫプリンタに出力している。このとき、ｓＲＧＢでは再現可能な（概ねディスプレイでは再現可能な）色であっても、ＣＭＹＫプリンタでは再現できない場合がある。このような色に関しては、色空間の変換の際に色域圧縮処理を行って、少なくともＣＭＹＫプリンタでは再現できない色について、再現可能な色に変換している。
【０００４】
色域圧縮処理としては、従来より種々の方法が案出されている。例えば、特開平５−２２７４１８号公報に開示されている方法では、明度を圧縮し、次いで色相を変えずに彩度方向に圧縮している。これに類似したものも多いが、基本的に、明度（Ｌ）、彩度（Ｃ）、色相（Ｈ）からなるＬＣＨ色空間で、色相を固定するなどの拘束条件を設定して、幾何学的にマッピングするものである。
【０００５】
一方、特公平５−４８６５７号公報に記載されている装置では、色域外の色については、色差距離を最小にする色に変換することが開示されている。さらに、特開平１０−８４４８７号公報や特開平６−１８９１２１号公報に開示されている方法は、基本的には明度、彩度、色相に重みを乗じることにより、より好ましい色に変換するものである。
【０００６】
上述の色域圧縮は、いかに入力画像中に含まれる再現できない色を出力装置の色域内の色で再現し、かつ、再現された色あるいは再現された画像が官能的に好ましいものであるかが重要である。例えば、黄色等では色相がずれてしまっても純色で再現されたほうが好ましい。また、色差が最小でなくとも、彩度が高めに再現された方が好ましい。このような場合、上述の色差最小、あるいは、ＬＣＨ空間において幾何学的規則で決定する方法では、必ずしも好適な色域圧縮が行えるわけではない。
【０００７】
この点で上述の特開平６−１８９１２１号公報に開示されている方法では、幾何学的方法、あるいは、一義的な色差最小による色域圧縮ではなく、上述のような好ましい色域圧縮を意図したものである。しかしながらこの文献に記載されている方法は、Ｃ，Ｍ，Ｙの３色を用いるプリンタの色域のみを取り扱うものである。また、色域内外の判定及び明度、彩度、色相の重みを定義する色空間が、基本的に機器依存のＲＧＢ色空間であり、そのために入力色信号に依存した色域の外側については決められないものである。つまり、ディスプレイの色域からＣＭＹプリンタへの色域圧縮に限定した場合はよいが、ＩＣＣ  Ｐｒｏｆｉｌｅ  Ｆｏｒｍａｔ（色変換用の色変換係数の形式の共通フォーマット、すでに、デファクトスタンダードと呼ばれ広く流通している）のような一般に用いられているＬａｂレンジ（０≦Ｌ* ≦１００、−１２８≦ａ* ≦１２７、−１２８≦ｂ* ≦１２７、以下ＩＣＣＬａｂレンジと呼ぶ）に対して、全ての領域で画像出力装置の色再現範囲への色変換を行うには不向きである。
【０００８】
また、上述の特開平１０−８４４８７号公報に開示されている方法は、全色領域で単一のＬＣＨ重み係数を用いるものであり、ＩＣＣＬａｂレンジ全域での色域圧縮を行うことはできる。しかし単一のＬＣＨ重み係数を用いているため、色によっては最適な色域圧縮を行うことができない。また、やはりこの文献においてもＣＭＹプリンタを想定しており、Ｃ，Ｍ，Ｙ，Ｋの４色で再現できる色域を有効に使うことは記載されていない。
【０００９】
さらに別の色域圧縮手法として、例えば特開平１０−２５８５４８号公報に記載されているように、入力色信号が色再現範囲に対して過明、等明度過彩、過暗のいずれかであるかを認識し、さらに、再現色の純色度合い、特に黄色の純色度合いに応じて色域圧縮を行うことも考えられている。特に黄色はＣＭＹＫプリンタにおいて高い明度再現が可能であるが、混色による明度低下も大きい。そのため、黄色などの色材色について、良好に再現するものである。しかし、このような色再現については、上述のように再現された色あるいは再現された画像が官能的に好ましいものであることが重要であり、純色度合いに限らず、いずれの色についてもさらに良好に再現されることが望まれていた。
【００１０】
【発明が解決しようとする課題】
本発明は、上述した事情に鑑みてなされたもので、好ましい良好な色域の変換を行うことができるカラー画像処理方法及びカラー画像処理装置と、変換の際に用いる色変換係数を生成する色変換係数生成方法及び色変換係数生成装置、及び、それらの処理をコンピュータに実行させるプログラムを格納した記憶媒体を提供することを目的とするものである。
【００１１】
【課題を解決するための手段】
本発明は、３次元の第１の色空間における色域から４次元の第２の色空間における色域への変換を行う際に、第１の色空間の色座標を表す３つの独立変数のうち少なくとも明度に対して重みを付けるとともに、重みの大きさを色相角に応じて変化させる。そして、重みを付けた変数を含む３つの独立変数により、第１の色空間における色域の外郭から、第２の色空間における色域を前記第１の色空間に変換したときの色域の外郭への距離が最小になる関係を求める。この関係から色変換係数を求めることができる。また、この関係に基づいて、第１の色空間から第２の色空間への色変換を行うことができる。
【００１２】
また、上述の関係あるいは色変換係数を第１の色空間における色域中及び色域の外郭上の複数の代表点について予め求めておき、例えばテーブル化しておき、第１の色空間の色信号に応じて予め求めておいた複数の関係あるいは色変換係数から補間によって第２の色空間における色信号を決定することができる。
【００１３】
３次元の第１の色空間から４次元の第２の色空間への変換を行う場合、第１の色空間における色は、第２の色空間における色に一対一に対応しない。しかし、第１の色空間における色に官能的に対応する第２の色空間の色は存在する。特に色域圧縮を行う場合には、第１の色空間における色が第２の色空間における色域外の場合には、強制的に他の色に変換することになるが、その場合でも官能的に対応する第２の色空間の色に変換することが望ましい。
【００１４】
本発明では、変換先として４次元の第２の色空間を想定して変換を行うことによって、第２の色空間の色域全体を用いた変換を行うことができる。例えばＣＭＹＫプリンタに出力する場合、Ｋを除くＣＭＹ色空間への変換では再現可能な色域を十分に利用できない。４次元のＣＭＹＫ色空間における色域を対象とすることによってＣＭＹＫ色空間で再現可能な色域を十分利用して、良好な色変換を行うことができる。
さらに本発明では、上述のように、第１の色空間の色座標を表す３つの独立変数のうち、少なくとも明度に対して、色相角に応じた重みを付ける。例えば重みを官能評価的手法などによって、色相角ごとに決定しておけば、官能的に好ましい良好な色域の変換を行うことができる。特に色域圧縮においては、明度・彩度・色相のどの成分を重視するかを、色相角毎に最適化できる。例えば、棒グラフなどのビジネスグラフと呼ばれるような画像で重要である特性、すなわち、黄色の領域では色相は多少犠牲にしてもよいから彩度を余り落とさずに、また暗い領域では彩度はほどほどで明度を重視して再現するといった要求を、明度・彩度・色相の重み係数の最適値を選択することで満足させることができる。
【００１５】
また本発明は、３次元の第１の色空間における色域から４次元の第２の色空間における色域への変換を行うための色変換係数を生成する色変換係数生成方法及び色変換係数生成装置において、前記第２の色空間における色域を前記第１の色空間に変換したときの色域の外郭である変換外郭を求め、前記第１の色空間における色域の外郭の色から前記第２の色空間における変換外郭への重み付きの距離が最小になるように前記第１の色空間の色座標を表す３つの独立変数のうち複数の変数に付ける重みを算出しておき、前記第１の色空間における色信号を順次入力し、算出しておいた複数の前記重みから該色信号に対応する重みを算出し、算出した色信号に対応する重みに基づいて入力した前記色信号を前記第２の色空間における色信号に変換し、該変換結果と入力した前記色信号とを関係付けることによって順次色変換係数を生成することを特徴とするものである。
さらに本発明のカラー画像処理装置は、上述のような本発明の色変換係数生成装置を用い、第１の色空間における色域中及び色域の外郭上の複数の代表点について求め、その色変換係数を格納し前記第１の色空間における色域の色信号が入力されると１ないし複数の色変換係数を出力する変換テーブルと、該変換テーブルから出力される１ないし複数の色変換係数に基づいて入力された前記第１の色空間における色域の色信号に対応する第２の色空間における色域の色信号を算出する出力色信号算出手段を有し、３次元の第１の色空間における色域の色信号から４次元の第２の色空間における色域の色信号への変換を行うことを特徴とするものである。
【００１６】
【発明の実施の形態】
まず本発明の概要について説明してゆく。本発明の説明に先立って、色差最小による色域の変換の原理を説明する。色域の変換を行う色空間をＬａｂ色空間とし、第１の色域が第２の色域より広いとする。第１の色域に属する点Ｔ（ＬT ，ａT ，ｂT ）が与えられ、第２の色域の表面（外郭）上にあるｎ個の色のうちのｉ番目の色のＬａｂ値をＰｉ（ＬPi，ａPi，ｂPi）とする。色差ΔＥは、
ΔＥｉ＝｛（ＬT −ＬPi）2 ＋（ａT −ａPi）2 ＋（ｂT −ｂPi）2 ｝1/2 …式１
で表される。
【００１７】
色差最小とは、与えられた点Ｔに対して、色差ΔＥｉが最も小さくなるような点Ｐｉを求めることである。このとき、第２の色域の表面上にあるｎ個の点の数は多い方がより正確になることは言うまでもない。
【００１８】
このようにすれば、第１の色域から、第２の色域の表面へのマッピングは一意に定まる。ただし、色差最小によるマッピングでは、明度をできるだけ保つ、あるいは、彩度をできるだけ保つ、あるいは、色相をできるだけ保つなどの調整がまったくできない。すなわち、自由度がないことがわかる。
【００１９】
本発明では、まず、Ｌａｂ色空間と式１で定義される色差ΔＥより導出される、明度の差（ΔＬ* ）、メトリッククロマの差（ΔＣ* ）と呼ばれる彩度の差、メトリックヒューの差（ΔＨ* ）と呼ばれる色相の差に相当する量に変換する。その後、さらに各項に重み係数を掛けることにより、明度をできるだけ保つ、あるいは、彩度をできるだけ保つ、あるいは、色相をできるだけ保つなどの調整を可能としている。
【００２０】
２点Ｔ（ＬT ，ａT ，ｂT ）、Ｐｉ（ＬPi，ａPi，ｂPi）間の色差をΔＥｉ、明度の差（ＬT −ＬPi）をΔＬｉとすれば、メトリッククロマの差（ΔＣ* ｉ）、メトリックヒューの差（ΔＨ* ｉ）は次のように表される。
ΔＣ* ｉ＝（ａT2＋ｂT2）1/2 −（ａPi2 ＋ｂPi2 ）1/2 
ΔＨ* ｉ＝｛ΔＥｉ2 −ΔＬｉ2 −ΔＣ* ｉ2 ｝1/2 
【００２１】
このようにして求められたΔＬ、ΔＣ* 、ΔＨ* と、重み係数α、β、γを用いて、重み付き色差をΔｗＥｉで表すと式２となる。但し、α、β、γともに０以上とする。
ΔｗＥｉ＝｛（γ×ΔＬｉ）2 ＋（α×ΔＣ* ｉ）2 ＋（β×ΔＨ* ｉ）2 ｝1/2     …式２
【００２２】
式２においてα＝β＝γ＝１であれば式１と同じで一般に言われている色差である。α＜γかつβ＜γであれば、相対的に、ΔＬｉの寄与が大きくなり、重み付き色差を最小にすることは明度の誤差を小さくする方向に働く。すなわち、明度をできるだけ保つことになる。αだけを、少なくとも１よりは大きくし、大きくすればするほど彩度をできるだけ保つ効果が得られる。同様にβだけを大きくすれば色相をできるだけ保つような効果が得られる。その程度は、重み係数α、β、γの大きさを調節することによって制御できる。ただし、その強さは相対的なものであるため、以下の説明ではγ＝１として説明する。
【００２３】
重み係数α、βの大きさの最適値は、Ｌａｂ色空間の領域毎に決定することが好ましい。色領域毎に重要とする指標が異なるからである。例えば、黄色領域であれば、彩度が重視され、白や黒の付近では明度・彩度・色相が均等に重視されるなどである。官能的にどのようなα、βの組み合わせがよいかを決定するには、第１の色域が再現可能な出力装置でＴ（ＬT ，ａT ，ｂT ）に相当する色を再現し、（α，β）を振って、第２の色域に属する複数の色を再現し、最も好ましいものを選択するという官能評価（Ｐｒｅｆｅｒｅｎｃｅ  Ｔｅｓｔ）を行って決定できる。
【００２４】
このような重み係数α、βを、代表的な色領域、例えば、第１の色域のＷ（白）、Ｋ（黒）と、Ｒ（赤）、Ｇ（緑）、Ｂ（青）、Ｃ（シアン）、Ｍ（マゼンタ）、Ｙ（イエロー）のもっとも彩度が高い色、及び、Ｗ，Ｋ，Ｒ，Ｇ，Ｂ，Ｃ，Ｍ，Ｙの中間的な彩度の色などについて求める。
【００２５】
これらの色について最適な重み係数α、βが求められた後、任意の与えられた点Ｔ（ＬT ，ａT ，ｂT ）に対する重み係数α，βを決定する。もとより、全ての与えられた点Ｔに対して前記の官能評価手法によって求められればよいが、相当の数について官能評価手法を行うのは非効率である。従って、既に求められている代表点をもとに、与えられた点Ｔでの重み係数α、βを補間により求める方法が有効である。
【００２６】
いま、予め最適な重み係数が求められている代表点の座標をＡj （Ｌj ，ａj ，ｂj ）、その時の重み係数をｗj （αj ，βj ）で表し、任意の与えられた点Ｔ（ＬT ，ａT ，ｂT ）、求めるべき重み係数をｗT （αT ，βT ）とする。また、補間には距離の重み付け平均を用いる。点Ａと点Ｔの距離をｄＴj で表すと、以下の式により  ｗT （αT ，βT ）を決定できる。
αT ＝Σ（αj ×（１／（１＋ｄＴj ／ｄｃ））m ）／（Σ（１／（１＋ｄＴj ／ｄｃ））m ）  …式３
βT ＝Σ（βj ×（１／（１＋ｄＴj ／ｄｃ））m ）／（Σ（１／（１＋ｄＴj ／ｄｃ））m ）  …式４
ここで、ｄｃとｍは常数である。重み付け平均は、代表点が与えられた点Ｔに近いほどその寄与が高くなる平均化手法である。また、重み付け平均は、外挿にも強く、官能評価で決定できる（現実に存在する画像出力装置で再現できる色域内の色でしか、官能評価はできない）重み係数をもとに、さらにその外側に存在する色座標の重みを決定する（外挿する）のに好適である。
【００２７】
このようにして任意の点Ｔ（ＬT ，ａT ，ｂT ）に対する重み係数α，βを決定した後、上述の式２において求められる重み付き色差ΔｗＥｉが最も小さくなるような点Ｐｉ（ＬPi，ａPi，ｂPi）を求めればよい。これによって重み付きの色域圧縮処理を行うことができる。
【００２８】
上述の説明では、第１の色域及び第２の色域ともＬａｂ色空間として説明した。基本的には同じ色空間内で色域圧縮処理を行うが、例えばＣＭＹＫプリンタなどのように４次元の色空間の色信号を扱う機器に対して画像を出力する場合、第２の色域はＣＭＹＫ色空間における色域が反映されていなければならない。Ｌａｂ色空間とＣＭＹＫ色空間は一対一には対応しないため、ＣＭＹＫ色空間における色域をＬａｂ色空間で表現するためには、例えば探索的な手法等を用いて色域を決定する必要がある。
【００２９】
以下、このような４次元のＣＭＹＫ色空間から３次元の色空間における色域の表面（外郭）を求め、上述のようにして重みを求めてゆく過程について、詳細に説明する。
【００３０】
図１は、本発明のカラー画像処理方法及び色変換係数生成方法の実施の一形態を示す工程図である。この工程図では、予め行っておく事前工程と、点Ｔｉ（ＬTi，ａTi，ｂTi）が与えられると、色域圧縮した後の点Ｏｉ（ＬOi，ａOi，ｂOi）を逐次算出する逐次工程からなる。事前工程は、▲１▼モデル生成工程、▲２▼色域表面算出工程、▲４▼ＬＣＨ重み代表点決定工程を含んでいる。また逐次工程は、▲３▼色域内外判定工程、▲５▼ＬＣＨ重み係数算出工程、▲６▼変換工程を含んでいる。以下、各工程について順に説明してゆく。
【００３１】
▲１▼モデル生成工程
ここではまず、デバイスモデルとインバースモデルについて説明する。デバイスモデルとは、機器に依存した色空間から機器に独立な色空間へ変換するものである。ＣＲＴなどでは加法混色が成立するので、ＲＧＢ各色ごとのべき乗による非線型変換、または、各色毎の１次元のルックアップテーブル（ＬＵＴ）による非線型変換と、ＣＲＴのＲＧＢ各１色毎での最大発光時の三刺激値に基づくマトリクスにより、デバイスモデルが表されることは広く知られている。これをＣＲＴ加法混色モデルと呼ぶことにする。一方、減法混色系の画像出力装置においては、濃度変調型のプリンタではランバート・ベア則で、面積変調型のプリンタではノイゲバウアー式などで表すことができることが知られている。
【００３２】
また、複数の色票を生成し、測色に基づき、デバイスモデルを生成することができる。例えば、特公昭６３−３２３１３号公報に公開されているように高次多項式で近似する方法、特開平２−２４１２７１号公報に開示されているようにニューラルネットワークで近似する方法等を用いることができる。あるいは、特開平１０−２６２１５７号公報に開示されているように、複数の色票を生成し、測色に基づき、重み付け平均により予測する方法を用いて、デバイスモデルを生成することもできる。これらの高次多項式、ニューラルネットワーク、重み付け平均によるデバイスモデルは、一般に数十から数千の色票数が必要となるものの、デバイスの種類、画像形成（表示）方法に制限がなく、汎用的であるとともに、デバイス毎に適切な色票数を選択すれば、高精度にデバイスモデルを構築することができる。
【００３３】
以下の説明において、特別のことわりがなければ、ＣＲＴでは、ＣＲＴ加法混色モデルをデバイスモデルとして用い、プリンタでは重み付け平均によるデバイスモデルを用いる。
【００３４】
なお、デバイスの色域内に限れば、ＣＲＴでは算術的に逆変換、すなわち、機器に独立な色信号からＣＲＴのＲＧＢ色信号への予測が可能である。ＣＭＹＫプリンタでは、色域内に限っても、一般に逆変換は多義となり解を求めることはできない。しかし、例えば特開平１０−２６２１５７号公報、あるいは、特開平７−８７３４６号公報に開示されているように、機器に独立な色信号と適当なＫ（黒）の値を与えることで、逆変換を解くことができる。この逆変換を解くモデルをインバースモデルと呼ぶことにする。ＣＲＴのデバイスモデルを式５に、ＣＲＴのインバースモデルを式６に、プリンタのデバイスモデルを式７に、プリンタのインバースモデルを式８による形式で表すことにする。
（Ｘ，Ｙ，Ｚ）＝ＦCRT （Ｒ，Ｇ，Ｂ）          …式５
（Ｒ，Ｇ，Ｂ）＝ＦCRT-1 （Ｘ，Ｙ，Ｚ）        …式６
（Ｌ，ａ，ｂ）＝ＦPRINT （Ｃ，Ｍ，Ｙ，Ｋ）    …式７
（Ｃ，Ｍ，Ｙ）＝ＦPRINT-1 （Ｌ，ａ，ｂ，Ｋ）  …式８
【００３５】
なお、式５〜式８における機器に独立な色信号は、ＣＲＴについては三刺激値ＸＹＺ、プリンタについては均等知覚色空間Ｌａｂとして例示している。三刺激値ＸＹＺと均等知覚色空間Ｌａｂの関係はＪＩＳ  Ｚ８７３０−１９８０によって定義されており、光源を決めれば双方向の変換が可能である。もちろん、他の表色系、例えば均等知覚色空間は、Ｌａｂ色空間の他、Ｌ* ｕ* ｖ* ，ＬＣＨなどを用いてもよい。
【００３６】
また、３色のＣＭＹプリンタの場合は、式７、式８の右辺においてＫ＝０とするか、Ｋそのものをなくした場合であって、式７、式８に含まれるものである。この場合、機器に独立な色信号とＣＭＹ色空間の関係は、一般に一対一であり、インバースモデルも簡単に解くことができる。
【００３７】
▲２▼色域表面算出工程
上述したデバイスモデルを用いて、色域を表現する方法を述べる。ＣＲＴのＲＧＢや、３色のＣＭＹプリンタでは、次の規則に基づいて簡単に色域表面を求めることができる。すなわち、Ｃ，Ｍ，Ｙの各色信号値を網点面積率とするとき、色域表面上の点（色）は、Ｃ，Ｍ，Ｙのうち少なくとも１色は０％または１００％であるという規則に着目する。図２は、ＣＭＹ色空間において色域表面上の点の一例の説明図である。例えばＣ色に着目すると、Ｃ＝０％の場合およびＣ＝１００％の場合に図２に示すようなＣ、Ｍ、Ｙの組み合わせが色域表面上の点である。なお図２では説明のため、Ｍ，Ｙの振り幅は５０％としている。この振り幅は細かい方がよい。Ｍ色，Ｙ色に関しても同様にしてデータが作成できる。
【００３８】
このようにして、色域表面の点を表すＣ，Ｍ，Ｙの組が得られ、このＣ，Ｍ，ＹをＫ＝０として、式７に代入する。これによって、対応するＬ，ａ，ｂを求めることができる。ＣＲＴのＲＧＢ色信号についてもＲ，Ｇ，Ｂの少なくとも１色は、無発光または最高輝度での発光として、まったく同様に色域表面のＲ，Ｇ，Ｂの組を生成でき、式５によりその時のＸ，Ｙ，Ｚを求めることができる。
【００３９】
さて、４色のＣＭＹＫプリンタの場合は、３色プリンタの場合とは異なり、探索による手法を用いないとＬａｂ色空間上での色域表面を求めることはできない。図３は、４色のＣＭＹＫプリンタにおけるＬａｂ色空間上での色域表面を求める方法の一例の説明図である。図３に示すように、ＣＭＹＫプリンタのＬａｂ色空間における色域内の定点Ａを起点として、任意のＩＣＣＬａｂレンジのへりの点をＢとする。この線分ＡＢ上の点Ｃは色域内外の境界である。なお、Ｌａｂ色空間上の色域は３次元であるが、図３ではその一断面としてＬ* −ａ* 平面のみを示している。
【００４０】
さて、点Ｃを見つけるためには、線分ＡＢ上の点をＰ（Ｌi ，ａi ，ｂi ）として式８を用いてＣＭＹを解く（式９）。ただし、ＣＭＹの範囲は０〜１００％である。
（Ｃj ，Ｍj ，Ｙj ）＝ＦPRINT-1 （Ｌi ，ａi ，ｂi ，Ｋj ）  …式９
ここで、Ｋj は点Ｐが決まる毎に０〜１００％の範囲で探索し、その都度求まる（Ｃj ，Ｍj ，Ｙj ，Ｋj ）を式７に代入して（Ｌj ，ａj ，ｂj ）を求める。求めた（Ｌj ，ａj ，ｂj ）と、目標色である（Ｌi ，ａi ，ｂi ）との色差を評価することで、点Ｐが色域内にあるのか外にあるのかが分かる。色差が小さければ、色域内と判定できる。例えば、色差０．３以内であれば色域内と判定することができる。
【００４１】
さて、点Ｐ（Ｌi ，ａi ，ｂi ）で色域内に属する（Ｃj ，Ｍj ，Ｙj ，Ｋj ）が複数存在することがあるが、その時は、Ｋj が最も大きいものを採用すればよい。この最も大きいＫj を最大墨量と呼ぶことにする。また、効率的な探索を行うには、Ｋj を求めるとき、さらに、点Ｐを移動するときに２分探査を用いるとよい。実際に、Ｌａｂ色空間中にＣＭＹＫプリンタの色域を求めるには、点Ａを中心に３次元的に放射状の線分を多数作成し、その線分毎に色域の境界を求めていけばよい。
【００４２】
このようにして、ＣＭＹプリンタ、ＣＭＹＫプリンタとも、色域表面の表色系の色座標で表される点の集合が求められる。なお、ここでの説明は表色系にＬａｂを用いたが、Ｌａｂ以外の表色系であっても同様に色域表面を求めることができる。
【００４３】
▲３▼色域内外判定工程
上述の▲２▼色域表面算出工程でも多少色域の内外判定について触れたが、ここでは色域の内外判定方法について詳しく述べる。例えばこの実施の形態を用いて色変換用の３次元テーブルを生成する場合などでは、色域の内外判定を行いながら、与えられた目標色の色座標が想定される出力装置の色域内であればマッピング先を自分自身にするなど、マッピングを制御するのに色域の内外判定の結果を用いることがある。▲２▼色域表面算出工程とは独立に色域の内外判定ができれば、より汎用性は高まる。
【００４４】
色域の内外判定は、デバイスモデルとインバースモデルを若干変形したものを用いて以下の手順で行うことができる。与えられた目標色を（Ｌ，ａ，ｂ）、想定する出力装置をＣＭＹＫプリンタとして説明する。インバースモデルの変形の考え方を述べると、▲２▼色域表面算出工程で述べたように、一般に目標色（Ｌ，ａ，ｂ）に対して、目標色を満たす（Ｃ，Ｍ，Ｙ，Ｋ）は複数存在することがある。しかし最大墨量Ｋｍａｘはただ１つ存在する、ということに着目すると、式１０の形式で記述できる。
（Ｃ，Ｍ，Ｙ，Ｋｍａｘ）＝ＦPRINT-1 （Ｌ，ａ，ｂ）  …式１０
もとより、式１０の関数ＦPRINT-1 は、前記のごとく多分に探索的手法によるもので、それらを概念的に表現したものである。また、式１０は、色域外に属する目標色に対しても、上述の▲２▼色域表面算出工程による線分上の探索により、色域表面上の（Ｃ，Ｍ，Ｙ，Ｋｍａｘ）を返すような機能を追加している。なお、色域外に属する目標色に対する扱いの変形例として、色差最小、あるいは、ＬＣＨ重み付き色差最小の色域表面上の（Ｃ，Ｍ，Ｙ，Ｋｍａｘ）を返すものであってもよい。
【００４５】
このようにして求められた（Ｃ，Ｍ，Ｙ，Ｋｍａｘ）をデバイスモデルに代入し、式１１に示すように（Ｌ’，ａ’，ｂ’）を求める。
（Ｌ’，ａ’，ｂ’）＝ＦPRINT （Ｃ，Ｍ，Ｙ，Ｋｍａｘ）  …式１１
さらに、式１２に示すように、目標色（Ｌ，ａ，ｂ）とデバイスモデルによって求められた（Ｌ’，ａ’，ｂ’）の色差ΔＥをとる。
ΔＥ＝｛（Ｌ−Ｌ’）2 ＋（ａ−ａ’）2 ＋（ｂ−ｂ’）2 ｝1/2   …式１２
式１２で求められたΔＥがεより大きければ色域外と判定することができる。ここで、εは色域内外を判定する閾値であって、デバイスモデル、インバースモデルの精度に依存して決めるべきものである。例えばεを１とすることができる。
【００４６】
さらに簡易に行う方法として、３次元テーブルを用いる方法が知られている。例えば、特公昭５８−１６１８０号公報に開示されている３次元テーブルと補間によって色変換を行う方法を色域内外判定に応用することが可能である。図４は、３次元テーブルで用いる色空間の分割の一例の説明図である。例えば、ＩＣＣＬａｂレンジ内の色空間を、例えば図４に示すように分割し、その格子点に、色域内なら０、色域外なら２５５といった値を記憶させておく。与えられたＬａｂ（格子点のアドレスに相当する）に対して、近傍の格子点に記憶されている色域内外を示すデータを読み出して補間演算を行い、与えられたＬａｂに相当する色域内外を示すデータを算出する。算出された色域内外を示すデータと、設定されている閾値ｒとを比較することによって、色域の内外を判定することができる。例えば、閾値ｒを１２８とすると、１２８よりも算出された色域内外を示すデータが大きければ色域外と判定し、１２８以下であれば色域内と判定することができる。
【００４７】
この方法は、予め、３次元テーブルに記憶しておくべき色域内外を示すデータを求めておけば、高速に色域内外判定ができるという利点がある。反面、補間演算は、色域の形状を考慮したものではないため、その判定には、ある程度の誤差が含まれる。３次元テーブルを細かくとることによって、この誤差を小さくできる。基本的には、画像処理時に、画素毎に色域の判定を行う必要がある場合など、高速性が求められるときに有効な方法である。
【００４８】
▲４▼ＬＣＨ重み代表点決定工程
代表点について、官能評価手法により、ＬＣＨに対する重み係数を求める。ここでは、ＣＲＴとゼログラフィープリンタの間での好ましいＬＣＨに対する重み係数を求める。図５は、ＣＲＴに表示させるソフト色表の信号値の一例の説明図である。図５に示す値は、コンピュータに接続されたＣＲＴを発光させるためのＲＧＢ値である。ここで行った官能評価は、ＣＲＴを設置したブース（ＣＲＴブース）、後述する方法でゼログラフイープリンタにより作成したハード色票を例えばＤ５０の蛍光燈光源で照らすことのできるブース（プリントブース）を併置する。そして、両眼隔壁法により、ソフト色票に対して好ましいハード色票を選定するという方法を採った。なお、コンピュータに接続されたＣＲＴは、各種調整が可能なものであり、白色表示の色温度はＤ５０に設定し、白色の輝度は、プリントブースに完全拡散面をおいた場合を想定した輝度に合わせた。このような条件のもとでは、図５に示すＲＧＢの値によるソフト色票は、対象としたゼログラフィープリンタの色域外であった。
【００４９】
次に、ハード色票の作成方法について説明する。ハード色票を作成するために、その目標色の色座標を算出する。図５の各ＲＧＢ値を（Ｒi ，Ｇi ，Ｂi ）とすると、目標色（Ｌi ，ａi ，ｂi ）はＣＲＴのデバイスモデルを用いて算出することができる（式１３）。
（Ｌi ，ａi ，ｂi ）＝ＦCRT （Ｒi ，Ｇi ，Ｂi ）  …式１３
【００５０】
次に、▲２▼色域表面算出工程で述べた方法により、色域表面の（Ｌj ，ａj ，ｂj ）、（Ｃj ，Ｍj ，Ｙj ，Ｋj ）の対を多数求める。これら、（Ｌj ，ａj ，ｂj ）、（Ｃj ，Ｍj ，Ｙj ，Ｋj ）の対全てを集合ＰＣＧＳ（Ｐｏｉｎｔｓ  ｏｆ  Ｃｏｌｏｒ  Ｇａｍｕｔ  Ｓｕｒｆａｃｅの略）、その要素をＰＣＧＳｊと表記することとする。
【００５１】
図６は、ＬＣＨ重み係数の組み合わせの一例の説明図である。図６に示すＬＣＨ重み全ての場合について、集合ＰＣＧＳのなかから、式２に従ってＬＣＨ重み付き色差ΔｗＥｉが最も小さくなるＰＣＧＳｊを求める。色差ΔｗＥｉが最も小さくなるＰＣＧＳｊが見つかることは、プリンタで色票を作成するための（Ｃj ，Ｍj ，Ｙj ，Ｋj ）が求められたことに他ならない。従って、プリンタで印刷すればハード色票を作成することができる。
【００５２】
図７は、官能評価によって得られたＬＣＨ重み係数の一例の説明図である。上述のような官能評価の結果、好ましいＬＣＨに対する重み係数として、例えば図７に示すように得られる。なお、図７の上６色は、図５に示したＲＧＢの組み合わせの色に対応するものである。また“ｈａｌｆ”と記されている色は、Ｒ、Ｇ、Ｂ、Ｃ、Ｍに対して、中程度の白（グレー）を足したものであり、例えばＲｈａｌｆは（Ｒ，Ｇ，Ｂ）＝（２５５，１２８，１２８）である。これらは、中間的な明るい色について、目標色を追加したものである。色空間中でより緻密に好ましいマッピングを得たいのであれば、このように目標色を加えて官能評価によりＬＣＨに対する重み係数を求め、追加すればよい。このようにして、代表的な色におけるＬＣＨ重み係数を決定することができた。
【００５３】
なお、図７に示した１１色に加えて、（Ｌ，ａ，ｂ）＝（１００，０，０）の白（Ｗ）、（Ｌ，ａ，ｂ）＝（０，０，０）の黒（Ｋ）については、どちらもＬＣＨ重み係数（α，β）＝（１，１）を追加しておく。白や黒の付近のマッピングは色差最小で充分である。また、この追加するＬＣＨ重み係数の色座標は、白であれば紙の白、黒であればプリンタの再現しうる最高濃度の色座標であってもよい。例えば、（Ｌ，ａ，ｂ）＝（５０，０，０）といった色域内部の点について、適当な重み係数（α，β）を追加してもよい。以上の追加は、後述する補間演算によるＬＣＨ重み係数予測を安定にするためのものである。このようにして、ＬＣＨに対する重みを求めたものをＬＣＨ重み代表点と呼ぶことにする。
【００５４】
さらに、プリントする紙が色付いていたり、３色または４色でプリントした最高濃度色が多少色づいていたりした場合は、多少の補正が必要になる場合がある。図８は、ＬＣＨ重み係数の補正が必要な場合の一例の説明図である。例えば図８に示すように、最大明度（Ｌ* 最大）のとき、紙の色などによって多少色付き、Ｌ* 軸からはずれている。また同様に、最小明度（Ｌ* 最小、最大濃度による印刷時）のとき、色材の色によってＬ* 軸からはずれている。このような場合には、図８に示すように直線Ｐより上の部分および直線Ｑより下の部分を予め切り取っておくことが好ましい場合がある。例えば、切り取らなかった場合の色差最小でのマッピング例がベクトルＭ，Ｎである。
【００５５】
このように、Ｌ* 軸と白、黒の頂点が著しくずれてしまう場合は、無彩色が色づいてしまうという不自然なマッピングになってしまう。そのため、図８に示すように上下の部分を切り取っておくことにより、このような不自然なマッピングを防止することができる。なお、図８は、説明の都合上、２次元で行ったが、実際には３次元である。
【００５６】
図８に示すように色域の上または下を切り取るためには、▲２▼色域表面算出工程で述べた方法に、次の過程を追加すれば良い。まず第１に、Ｌ* 軸とプリンタの色域との交点を求める。この交点の求め方は、白方向、黒方向にＬ* 軸上を探索すれば求めることができる。白方向でのＬ* 軸との交点をＬｗ、黒方向でのＬ* 軸との交点をＬｄとする。次に、色域表面の色座標の集合ＰＣＧＳを作成する際に、Ｌｗより大きいＬに対しては全てＬｗとし、Ｌｄより小さいＬはＬｄとし、ａ，ｂに関してはそのままとしておけばよい。この操作を色域整形と呼ぶ。なお、ＬＣＨ重み代表点を官能評価で求めるとき、白、黒の近傍を入れた場合は、事前に色域整形を行っておくことが必要である。
【００５７】
上述のＬＣＨ重み係数の決定方法では、ＣＲＴの色域を第１の色域、ゼログラフィープリンタの色域を第２の色域とみなすことができる。このときのマッピングの対応関係について考えてみる。図９は、第１の色域と第２の色域の関係の一例の説明図である。図９では、第１の色域と第２の色域を、Ｌ* 軸と、ある色相におけるＣ* 軸の２次元断面で模式的に示している。上述の説明では、第２の色域が第１の色域に包含されるものとしていたが、実際には図９に示すように、第２の色域の一部が第１の色域に含まれない場合も存在する。図９において、第２の色域の中であって第１の色域との共通領域を領域Ａ、第２の色域の中であって第１の色域外の領域を領域Ｂとする。上述のＬＣＨ重み係数によるマッピング方法は、第１の色域から第２の色域へのＬＣＨ重み係数を用いた最短距離でマッピングする。そのため、図９に示すような場合には、領域Ａにマッピングされる部分については色域の圧縮方向に、領域Ｂにマッピングされる部分については色域の伸長方向にマッピングされることが分かる。
【００５８】
また、領域Ａのように第２の色域が第１の色域に包含される場合は、色域の圧縮しかあり得ない。しかし、領域Ｂのように第２の色域のうち第１の色域に含まれない部分については、伸長を行わない方が好ましいことがある。例えば忠実な色再現が望まれるような場合がこれに当たる。その場合には、第１の色域と第２の色域の共通領域、すなわち、領域Ａを第２の色域とすることで、色域の伸長を行わないマッピングが得られる。
【００５９】
Ｌａｂ色空間上での第１の色域と第２の色域の共通領域は、例えば、（Ｌ，ａ，ｂ）＝（５０，０，０）のような第１の色域と第２の色域の共通領域に属する収束点Ｐ０を決めて、収束点Ｐ０より、半直線を任意の方向に伸ばす。その半直線と第１の色域外郭との交点Ｐ１と、同じく、その半直線と第２の色域外郭との交点Ｐ２とを求める。求めた交点Ｐ１，Ｐ２のうち、収束点Ｐ０に近い方を色域外郭のデータＰとする。前記半直線をあらゆる方向に向かって引いて、その都度、Ｐを求めてゆけばよい。これによって、第１の色域と第２の色域の共通領域が求められたことになる。このように、色域の伸長が可能な場合は、実際に伸長を行うか否かの選択が可能である。
【００６０】
▲５▼ＬＣＨ重み係数算出工程
ここでは、色空間中の任意の色座標におけるＬＣＨ重み係数を決定する。上述の▲４▼ＬＣＨ重み代表点決定工程で説明したように、図７に示した１１色、あるいはそれに加えて白（Ｗ）、黒（Ｋ）、その他いくつかの色について追加されたＬＣＨ重み代表点が求められている。
【００６１】
Ｌａｂ色空間中の任意の色座標におけるＬＣＨ重み係数を決定するには、ＬＣＨ重み代表点をもとに、上述の式３、式４を用いて重み付け平均により求めることができる。任意に与えられたＬａｂは（ＬT ，ａT ，ｂT ）に相当し、式３、式４中のｄTjは（ＬT ，ａT ，ｂT ）とＬＣＨ重み代表点（Ｌj ，ａj ，ｂj ）との距離である。また常数ｄｃ、ｍとしては、例えばｄｃ＝４０、ｍ＝４とすることができる。
【００６２】
このように、ＬＣＨ重み代表点との距離に応じたＬＣＨ重み計数を用いることによって、代表色の個数は少なくても色空間中の任意の座標で重み係数が決定できる。また、色空間中の代表色の疎密があって色空間中の任意の座標で重み係数が決定でき、そして、それらは、どのような方向に対しても重み係数が滑らかに推移する。さらに、外挿適応性が高いために、実際に官能評価ができないような色座標でも重み係数を決定することができる。
【００６３】
▲６▼変換工程
上述の▲５▼ＬＣＨ重み係数算出工程において、色空間中の任意の色座標におけるＬＣＨ重み係数（α，β）を決定することができた。このＬＣＨ重み係数（α，β）を用い、例えば上述の式２に適用して、任意の色座標Ｔ（ＬT ，ａT ，ｂT ）から重み付き色差ΔｗＥｉが最小となるＰｉ（ＬPi，ａPi，ｂPi）を求めればよい。このようにして、Ｌａｂ色空間中の任意の色について、４色のＣＭＹＫプリンタにおける色域への色域圧縮処理を行うことができる。
【００６４】
なお変換工程における処理は、例えば▲３▼色域内外判定工程において色域内であると判定された場合には、その色座標をそのまま出力し、色域外であると判定された場合にのみ、色域圧縮処理による色変換を行うことができる。
【００６５】
上述の説明では、ＬＣＨ重み係数α及びβを変化させる例を示した。上述の特殊な例として、α＝βとして変化させることもできる。すなわち、上述の式２においてα＝β＝１としてγを変化させることと等価である。図１０は、設定する明度重み係数γと色相角との関係の一例の説明図である。図１０では、色相角に応じて明度重み係数γを設定する場合の一例を示している。上述のような▲４▼ＬＣＨ重み代表点決定工程における官能評価手法により、明度重み係数γを求めることにより、例えば図１０に示すように明度重み係数γを決定することができる。図１０に示す例では、ｓＲＧＢのＣＲＴとＣＭＹＫプリンタの間の色域変換を行う場合を想定している。このとき、緑（Ｇ：色相角１３２°）の明度重み係数γを１．５にし、シアン（Ｃ：色相角１９６°）へ向けて明度重み係数γを次第に大きくした関数として、明度重み係数γが得られた。もちろん、重み係数α、βを用いる場合には、この関数とは上下が逆の形状の関数として得られる。
【００６６】
図１１は、明度重み係数の変化によるマッピングの変化の一例の説明図である。図１１では、ＣＲＴの色域とプリンタの色域を、Ｌ* 軸と、ある色相におけるＣ* 軸の２次元断面で模式的に示している。明度重み係数γを大きくすると、相対的に明度を保存するように、上述の式２における重み付き色差ΔｗＥｉが計算される。そのため、図１１において短い破線で示したように明度変化の小さい色へのマッピングが行われる。また、明度重み係数γを小さくすると、逆に他の要素、ここでは彩度も考慮したマッピングが行われる。そのため、図１１において長い破線で示したようにマッピングされる。
【００６７】
図１０に示した明度重み係数γの設定では、上述のように緑の明度重み係数γを小さくしている。ＣＲＴでは緑（Ｇ）や黄（Ｙ）の明度が非常に高いので、明度を保存して色域圧縮を行うと、図１１の短い破線で示したように、明度は保存されるものの、彩度が非常に低下してしまう。そのため、明度重み係数γを小さくして、ある程度、彩度を確保している。
【００６８】
ここでは明度重み係数の大小によるマッピングの変化を示したが、上述のように重み係数α、βを独立して変化させる場合には、明度方向及び彩度方向とともに、色相方向にも変化することになる。その変化量は、重み係数αとβの値の違いが小さければ色相が保存される方向にマッピングされ、α＝βであれば図１１に示すように色相は変化しない。また、αとβの値の違いが大きくなれば、色相の変化量も大きくなる。
【００６９】
図１２は、設定する明度重み係数γと色相角との関係の別の例の説明図である。図１０に示した明度重み係数γの変化は一例であり、必ずしも図１０に示すように変化するわけではない。また、官能評価のみで決定するほか、人為的に一部を調整することも可能である。例えば図１２（Ａ）に示すように、ある特定の色付近において明度重み係数γの値を変化させることができる。また、図１２（Ｂ）に示すように、広い色範囲において明度重み係数γの値を変化させてもよい。さらには、例えば図１２（Ｃ）に示すように、さらに複雑な形状で明度重み係数γの値を変化させてもよい。
【００７０】
また、上述の例では、彩度による明度重み係数γの変化は考慮していないが、もちろん、彩度によって明度重み係数γを変化させてもよい。なお、重み係数α、βも独立に変化させる場合についても同様である。
【００７１】
上述のようにして、画像処理時に本発明を直接適用し、入力された色信号に対してＬＣＨ重み係数を求め、色域表面のＬＣＨ重み係数による重み付き距離最小点を探索することによって、色域圧縮後の色信号を得ることができる。しかし、このような処理を、入力された画像信号に対して逐次行うよりも、例えば図４を用いて上述したように予め３次元テーブルを作成しておき、この３次元テーブルを用いて実際の画像を処理すると、より効率的に画像処理を行うことができる。ここでは３次元テーブルの作成方法を説明する。なお、以下に示す３次元テーブルの作成方法では、単純に色域外の入力値（Ｌ，ａ，ｂ）はプリンタの色域表面にマッピングし、色域内に関しては、忠実に再現するものとする。
【００７２】
図４に示すように、ＩＣＣＬａｂレンジを分割し、格子点ごとに上述の色域圧縮処理を行ってゆく。ＩＣＣＬａｂレンジの分割は、例えば、各軸８分割することができ、この分割により得られる９代表点で、総計７２９格子点について、それぞれの色域圧縮後のＣＭＹＫ色空間による色信号値を取得し、色変換係数として格子点に対応付けて格納する。もちろん、分割数及び格子点数は任意である。
【００７３】
各格子点毎の処理を説明する。まず、格子点（ＬTi，ａTi，ｂTi）を与え、▲３▼色域内外判定工程において色域内外判定を行う。色域の内部であれば、格子点そのものを色域圧縮後のＬａｂ色空間における色信号とする。また、色域外であれば、格子点（ＬTi，ａTi，ｂTi）におけるＬＣＨ重み係数を▲５▼ＬＣＨ重み係数算出工程で求め、そのＬＣＨ重み係数による色域表面上の距離最小点（ＬPi，ａPi，ｂPi）を▲６▼変換工程で見つける。このようにして、各格子点におけるＬａｂ色空間での色域圧縮後の色信号が得られる。得られた色信号は、全てプリンタの色域内に圧縮されているので、上述の式１０により、（ＣPi，ＭPi，ＹPi，ＫｍａｘPi）を求めることができる。つまり、ＩＣＣＬａｂレンジ全域に対して（ＣPi，ＭPi，ＹPi，ＫｍａｘPi）が求められたことになる。
【００７４】
なお、３次元テーブルには、各格子点に対応して格納する色変換係数として、Ｌａｂ色空間におけるデータ、すなわち格子点（ＬTi，ａTi，ｂTi）そのものまたは距離最小点（ＬPi，ａPi，ｂPi）を格納してもよいし、あるいはＣＭＹＫ色空間に変換した（ＣPi，ＭPi，ＹPi，ＫｍａｘPi）を格納してもよい。
【００７５】
また、ＵＣＲ（Ｕｎｄｅｒ  Ｃｏｌｏｒ  Ｒｅｍｏｖａ１：下色除去）を任意に設定したければ、０≦ρ≦１となるようなＵＣＲ関数ρ＝ＵＣＲ（Ｌ，ａ，ｂ）を定義して、以下のように（ＣPi’，ＭPi’，ＹPi’，ＫPi）を決め直せばよい。
ρPi＝ＵＣＲ（ＬPi，ａPi，ｂPi）      …式１４
ＫPi＝ρPi×ＫｍａｘPi                …式１５
（ＣPi’，ＭPi’，ＹPi’）＝ＦPRINT-1 （ＬPi，ａPi，ｂPi，ＫPi）  …式１６
なお、関数ＵＣＲは、Ｌ* 軸（グレー軸）に近いほどＵＣＲ率を高く、明度が低くなるに従ってＵＣＲ率を高く、などのように予め決められているものである。
【００７６】
なお、厳密には、色域表面は４色色域で求めているので、３色色域と４色色域の間に属する領域では、式１６が求まらない場合がある。この場合は、ＫPiを固定して、４色色域上の（ＬPi，ａPi，ｂPi）から、色差最小などの方法で、（ＣPi’，ＭPi’，ＹPi’）がどの要素も０から１００％に収まるような点を探索することにより決定できる。簡易的には、そのような点をＫｍａｘとしておいてもよい。
【００７７】
上述のカラー画像処理方法及び色変換係数生成方法は、コンピュータプログラムによって実現することが可能である。その場合、そのプログラムおよびそのプログラムが用いるデータなどは、コンピュータが読み取り可能な記憶媒体に記憶することも可能である。記憶媒体とは、コンピュータのハードウェア資源に備えられている読取装置に対して、プログラムの記述内容に応じて、磁気、光、電気等のエネルギーの変化状態を引き起こして、それに対応する信号の形式で、読取装置にプログラムの記述内容を伝達できるものである。例えば、磁気ディスク、光ディスク、ＣＤ−ＲＯＭ、コンピュータに内蔵されるメモリ等である。
【００７８】
図１３は、本発明のカラー画像処理装置及び色変換係数生成装置の一例を示すブロック図である。図中、１はモデル生成部、２は色域算出部、３は色域内外判定部、４は重み決定部、５は重み算出部、６は変換部である。各部は図１に示した各工程に対応しており、モデル生成部１は▲１▼モデル生成工程を、色域算出部２は▲２▼色域表面算出工程を、色域内外判定部３は▲３▼色域内外判定工程を、重み決定部４は▲４▼ＬＣＨ重み代表点決定工程を、重み算出部５はＬＣＨ重み係数算出工程を、変換部６は▲６▼変換工程を、それぞれ実現する。
【００７９】
カラー画像処理装置として用いる場合には、任意の画像信号（ＬTi，ａTi，ｂTi）を色域内外判定部３に与え、変換部６から色域圧縮後の色信号（ＬOi，ａOi，ｂOi）を得ればよい。
【００８０】
また、色変換係数生成装置として用いる場合には、作成する３次元テーブルの格子点における色信号（ＬTi，ａTi，ｂTi）を順次色域内外判定部３に入力し、変換部６から得られた色域圧縮後の色信号（ＬOi，ａOi，ｂOi）あるいはこれをＣＭＹＫ色空間に変換した色値（ＣPi，ＭPi，ＹPi，ＫｍａｘPi）を各格子点における色変換係数として格納すればよい。なお、各格子点における色変換係数を順次求めるため、添え字ｉをインクリメントしてゆく操作も色変換係数生成装置で行って、すべての格子点データが順次得られるように構成してもよい。
【００８１】
なお、図１より分かるように、▲１▼モデル生成工程、▲２▼色域表面算出工程、▲４▼ＬＣＨ重み代表点決定工程は事前に行って準備しておくことができる。従って、プリンタなどの出力装置が決まっていれば、事前工程に属する各工程は一度だけ行っておけばよい。そのため、図１３において、これらの工程に対応するモデル生成部１，色域算出部２，重み決定部４については、予め実行しておけば装置内に存在しなくてもよい。例えば重み決定部４で決定される重み代表点のデータをテーブルなどで重み算出部５に供給するようにして構成することもできる。
【００８２】
図１４は、本発明のカラー画像処理装置の別の例を示すブロック図である。図中、１１は３次元変換テーブル、１２は出力色信号算出部、１３は記録部である。３次元変換テーブル１１は、上述のような色変換係数生成方法及び色変換係数生成装置を用いて作成した３次元テーブルである。例えば図４に示したようにＬａｂ色空間における各格子点について、その格子点に対応するアドレスに色変換係数が格納されている。画像信号が入力されると、その画像信号がアドレスとして３次元変換テーブル１１に入力される。これによって、画像信号に近い１ないし複数の格子点の色変換係数が出力色信号算出部１２に出力される。
【００８３】
出力色信号算出部１２は、３次元変換テーブル１１から出力される１ないし複数の色変換係数を受け取り、入力された画像信号に従って例えば補間演算を行い、入力された画像信号に対応する出力色信号を算出して出力する。なお、３次元変換テーブル１１に色変換係数として例えばＬａｂ色空間における値が格納されている場合には、ＣＭＹＫ色空間への変換処理も行うとよい。
【００８４】
このようにして出力色信号算出部１２で算出された出力色信号は、記録部１３に渡されてＣ，Ｍ，Ｙ，Ｋの４色、あるいはこの４色から導出されるさらに多数の色を用いて画像を形成することができる。
【００８５】
図１５は、本発明のカラー画像処理装置の実現例を含むシステムを示すブロック図である。図中、２１はコンピュータ、２２はディスプレイ、２３はプリンタＡ、２４はスキャナ、２５はネットワーク、２６はネットワークプリンタＢ、２７はネットワークプリンタＣ、２８は測色計Ａ、２９は測色計Ｂである。図１５に示すシステムでは、コンピュータ２１にディスプレイ２２、プリンタＡ２３、スキャナ２４が画像入出力装置として接続され、また、ネットワーク２５を介してネットワークプリンタＢ２６及びネットワークプリンタＣ２７と接続されている。さらにコンピュータ２１には、ディスプレイ２２の発する色を測色するための測色計Ａ２８と、反射物（記録物等）を測色するための測色計Ｂ２９が接続されている。このシステムは大きく２つの機能を有している。１つは色変換係数を生成するものであり、もう一つは入力された画像を出力装置用に実際に色変換する機能である。
【００８６】
まず、ディスプレイ２２に表示された画像をプリンタＡ２３に出力するための色変換係数を生成する場合について説明する。コンピュータ２１は、１７３色のソフト色票をディスプレイに順次表示し、その都度、測色計Ａ２８で、その時の測色値を取り込んでいく。ソフト色表は、例えばＲ，Ｇ，Ｂ単色で１６レベルずつの４８色と、Ｒ，Ｇ，Ｂ各５レベルの直積１２５色による１７３色について作成することができる。もちろん、色数はこれに限らない。
【００８７】
全てのソフト色表について測色計Ａ２８で測色値を取り込んだら、ブラックボックスモデルにて、ディスプレイ２２のデバイスモデル及びインバースモデルを作成する。ブラックボックスモデルを用いるのは、ディスプレイがＣＲＴでない場合も想定してのことである。
【００８８】
デバイスモデルを用いて、ディスプレイ２２に依存した色信号の色空間であるＲＧＢ色空間から、機器に独立したＬａｂ色空間への色変換係数を作成する。色変換係数は、Ｒ，Ｇ，Ｂの各色信号をそれぞれ例えば１次元ＬＵＴにより明度がリニアのＲ’Ｇ’Ｂ’色空間の色信号Ｒ’，Ｇ’，Ｂ’に変換する。そして、この色信号Ｒ’，Ｇ’，Ｂ’を３次元テーブルにてＬａｂ色空間の色信号Ｌ，ａ，ｂに変換する。このように、ＲＧＢ色空間からＬａｂ色空間への変換は、１次元ＬＵＴと３次元ＬＵＴのセットにより行うことができる。３次元ＬＵＴは代表点のみを記憶しておき、その間の値は、近傍の代表点より補間して求める。このとき行う補間処理は基本的に線形であるので、代表点と求めるＬａｂ色空間の色信号との間は、線形に近い方が補間誤差を小さくすることができる。そのため、ＲＧＢ色空間から、一旦、明度がリニアなＲ’Ｇ’Ｂ’色空間の色信号Ｒ’，Ｇ’，Ｂ’に変換している。なお、３次元テーブルを作る際の、３次元テーブルのサイズ（分割数）や、Ｌａｂレンジは、例えばオペレータが指示する等により、任意に設定することができる。
【００８９】
次に、プリンタＡ２３のデバイスモデルとインバースモデルを生成するための手順を述べる。まず、Ｃ，Ｍ，Ｙ，Ｋ各単色階調についてプリンタＡ２３から出力し、ハード色表を得る。このハード色表は、例えばＣ，Ｍ，Ｙ，Ｋの各色について３２色ずつ、計１２８色について作成することができる。もちろん、色数は任意である。これらのハード色表を作成するためのデータは、予めコンピュータ２１にファイルとして記憶させておくことができ、コンピュータ２１がプリンタＡ２３にデータを転送してハード色表を作成することができる。
【００９０】
作成されたハード色表は、測色計Ｂ２９により測色し、コンピュータ２１に取り込む。この測色したデータに基づいて、紙の白の測色値と各階調の測色値の色差（ΔＥｗｐ：ΔＥ  ｗｈｉｔｅ  ｐｏｉｎｔ）をとる。そして、プリンタＡ２３のＣ，Ｍ，Ｙ，Ｋ（面積変調型のプリンタなら網点面積率、濃度変調型のプリンタなら濃度値の信号）が色差ΔＥｗｐに対してリニアになるように、Ｃ，Ｍ，Ｙ，Ｋの各単色の１次元ＬＵＴを設計する。画像を処理するための１次元ＬＵＴの入出力の関係は、色差ΔＥｗｐに対してリニアなＣＭＹＫ色空間からプリンタのＣＭＹＫ色空間への１次元ＬＵＴである。これをΔＥｗｐリニアＬＵＴと呼ぶ。通常、印刷機やプリンタのＣＭＹＫによるＴＲＣ（Ｔｏｎｅ  Ｒｅｐｒｏｄｕｃｔｉｏｎ  Ｃｕｒｖｅ）において、階調の出始めでの飛び（Ｗａｓｈｏｕｔ）や高濃度部でのつぶれが発生する。この現象は、電子写真学会年次大会（通算７３回）“Ｊａｐａｎ  Ｈａｒｄｃｏｐｙ  ‘９４”論文集，Ｐ．１８２にあるように、インバースモデルを生成する際に都合が悪い。ΔＥｗｐリニアＬＵＴは、この現象を回避するためと、ＬａｂとΔＥｗｐに対してリニアなＣＭＹＫとの間の線形性を高めるのに好適である。
【００９１】
次に、所定のＣＭＹＫ色空間の色信号をΔＥｗｐリニアＬＵＴを通してプリンタＡ２３で印刷し、印刷結果を測色計Ｂ２９で測色して、その結果をコンピュータ２１に取り込む。与えるＣＭＹＫ色空間の色信号は、例えば予めコンピュータ２１にファイルとして記憶しておくことができる。色数としては、例えば９２８色とすることができる。
【００９２】
コンピュータ２１に取り込んだ測色結果について、プリンタＡ２３のデバイスモデルとインバースモデルが構築できる。色変換係数を作るオペレータは、ＵＣＲ関数を選択あるいはＵＣＲ関数を生成するためのパラメータ、テーブルのサイズ（分割数）に関わるパラメータ、Ｌａｂレンジを決めるためのパラメータ等である。これらのパラメータを入力することにより、ＬａｂからＣＭＹＫへの色変換用の３次元テーブルを作成することができる。
【００９３】
このようにして、プリンタ用の色変換係数はＬａｂ色空間からΔＥｗｐに対してリニアなＣＭＹＫ色空間への３次元テーブルと、ΔＥｗｐに対してリニアなＣＭＹＫ色空間からプリンタのＣＭＹＫ色空間への１次元ＬＵＴのセットが作成される。なお、プリンタの色域外の色について色域圧縮を行うためのＬＣＨ重み代表点は、例えば図１における事前工程で求めたものを予め記憶しておき、上述の手順でΔＥｗｐに対してリニアなＣＭＹＫ色空間への３次元テーブルを求めている。その際、図１中の事前工程に属する各工程は、プリンタのデバイスモデル、インバースモデルを構築する際に行ない、デバイスモデルなどと一緒に記憶される。そして、採取したデバイスモデル用の色票の測色値が更新されるまで、繰り返し呼び出し使うことができる。厳密には、使用するディスプレイとプリンタの色域は、それぞれの機種、設定条件等により、また同じ機種でも経時変動などによって、その形は異なる。そのため、使用するディスプレイとプリンタにおいて、そのときの状態で官能評価を行い、ＬＣＨ重み代表点を決定することが好ましい。しかし実用上は、上述のようにディスプレイ及びプリンタを特定してＬＣＨ重み代表点を決定したり、あるいは、代表的なディスプレイ（例えばＣＲＴ）と、代表的な印字方式、例えば、印刷、銀塩写真、ゼログラフイープリンタ、インクジェットなどでのＬＣＨ重み代表点を作成しておいて、プリンタの印字方式に応じて選択することも有効である。
【００９４】
このようにして、ディスプレイ２２に表示するためのＲＧＢ信号をＬａｂ信号に変換して、このＬａｂ信号をプリンタＡ２３のＣＭＹＫ信号に変換するための色変換係数が作成される。これらの色変換係数はコンピュータ２１に記録される。
【００９５】
次に、実際に画像処理する場合に付いて説明する。オペレータは、ディスプレイのＲＧＢ色空間における色信号で構成される画像ファイルと、その画像ファイルが特定のディスプレイ（例えばディスプレイ２２）に対するものであるとの指示を、プリンタＡ２３を駆動するためのプリンタドライバに渡す。なお、プリンタドライバは、コンピュータ２１内で実行される。プリンタドライバには、予め上述のようにして生成された色変換係数が登録されており、順次色変換を実施し、画像データをプリンタＡ２３に送出して、プリント出力が得られる。この過程で順次行われる色変換は、１次元ＬＵＴ、３次元ＬＵＴ、３次元ＬＵＴ、１次元ＬＵＴによって行われる。しかし、途中の２つの３次元ＬＵＴは合成が可能であり、実際には、１次元ＬＵＴ、３次元ＬＵＴ、１次元ＬＵＴで画像処理を行うことができる。
【００９６】
ここでは、コンピュータ２１内のプリンタドライバが色変換を行ったが、例えば、色変換係数と未処理の画像データをプリンタＡ２３に送出して、プリンタＡ２３側で処理を行うことも可能である。また、コンピュータ２１で実行されるアプリケーションソフトの中で処理して、ＣＭＹＫの色信号になった画像データをプリンタドライバに渡す構成も可能である。色変換係数と未処理の画像データをプリンタに送出する場合は、例えばＰｏｓｔＳｃｒｉｐｔ（商品名）のようなＰＤＬ（Ｐａｇｅ  Ｄｅｓｃｒｉｐｔｉｏｎ  Ｌａｎｇｕａｇｅ）により画像データに色変換係数を添付し、一つのファイルとして取り扱ってもよい。
【００９７】
上述の説明では、ディスプレイ２２のＲＧＢ色空間の色信号で構成された画像データをプリンタＡ２３で印刷する場合について示した。このほか、例えば図１５に示した構成では、同じくディスプレイ２２のＲＧＢ色空間の色信号で構成された画像データを、ネットワーク２５を介してネットワークプリンタＢ２６やネットワークプリンタＣ２７で印刷出力する場合についても、同様にしてそれぞれのプリンタのデバイスモデル、インバースモデルに基づいて色変換係数を求め、色変換を行って印刷出力させることができる。さらに、スキャナ２４で取り込んだ画像をディスプレイ２２やプリンタＡ２３，ネットワークプリンタＢ２６，ネットワークプリンタＣ２７に出力する場合についても同様である。
【００９８】
もちろん、システム構成は図１５に示した構成に限られるものではなく、種々の構成のシステムにおいて、本発明を適用することが可能である。
【００９９】
【発明の効果】
以上の説明から明らかなように、本発明によれば、入力画像あるいは入力画像の属する３次元の色空間から、４次元の色空間の色信号で画像を再現する出力装置において再現可能な色域への変換を行う際に、色空間の部分空間毎に最適な方向（明度方向、彩度方向、色相方向）へ変換を行うことができる。そのため、出力装置から好ましい出力画像を得ることができる。また、テーブル型の色変換係数を生成する場合などにおいては、表色系の色座標としては存在するが現実には存在しないような色まで無理なく、好適な色変換係数を得ることができるという効果がある。
【図面の簡単な説明】
【図１】  本発明のカラー画像処理方法及び色変換係数生成方法の実施の一形態を示す工程図である。
【図２】  ＣＭＹ色空間において色域表面上の点の一例の説明図である。
【図３】  ４色のＣＭＹＫプリンタにおけるＬａｂ色空間上での色域表面を求める方法の一例の説明図である。
【図４】  ３次元テーブルで用いる色空間の分割の一例の説明図である。
【図５】  ＣＲＴに表示させるソフト色表の信号値の一例の説明図である。
【図６】  ＬＣＨ重み係数の組み合わせの一例の説明図である。
【図７】  官能評価によって得られたＬＣＨ重み係数の一例の説明図である。
【図８】  ＬＣＨ重み係数の補正が必要な場合の一例の説明図である。
【図９】  第１の色域と第２の色域の関係の一例の説明図である。
【図１０】  設定する明度重み係数γと色相角との関係の一例の説明図である。
【図１１】  明度重み係数の変化によるマッピングの変化の一例の説明図である。
【図１２】  設定する明度重み係数γと色相角との関係の別の例の説明図である。
【図１３】  本発明のカラー画像処理装置及び色変換係数生成装置の一例を示すブロック図である。
【図１４】  本発明のカラー画像処理装置の別の例を示すブロック図である。
【図１５】  本発明のカラー画像処理装置の実現例を含むシステムを示すブロック図である。
【符号の説明】
１…モデル生成部、２…色域算出部、３…色域内外判定部、４…重み決定部、５…重み算出部、６…変換部、１１…３次元変換テーブル、１２…出力色信号算出部、１３…記録部、２１…コンピュータ、２２…ディスプレイ、２３…プリンタＡ、２４…スキャナ、２５…ネットワーク、２６…ネットワークプリンタＢ、２７…ネットワークプリンタＣ、２８…測色計Ａ、２９…測色計Ｂ。
【図１】
0003981790.tif 000002
【図２】
0003981790.tif 000003
【図３】
0003981790.tif 000004
【図４】
0003981790.tif 000005
【図５】
0003981790.tif 000006
【図６】
0003981790.tif 000007
【図７】
0003981790.tif 000008
【図８】
0003981790.tif 000009
【図９】
0003981790.tif 000010
【図１０】
0003981790.tif 000011
【図１１】
0003981790.tif 000012
【図１２】
0003981790.tif 000013
【図１３】
0003981790.tif 000014
【図１４】
0003981790.tif 000015
【図１５】
0003981790.tif 000016