JP2007228153A - 色変換装置、色変換方法および色変換プログラム - Google Patents

Abstract

【課題】 画像データに適合した修正を高速に行いつつ、画像の見えを正確に再現する。
【解決手段】 第一色変換工程においてはディスプレイ６０に表示された画像の色を、その観察環境下においてどのように知覚されるかを推定することができる。そして、第二色変換工程においては印刷結果の観察環境において、第一色変換工程にて推定された知覚色と同様に知覚される色が印刷できるＣＭＹＫ値を特定することができる。また、ガマット修正をテーブルデータのＧＭＴ１３ｇを参照して行うことにより、処理効率を向上させることができる。さらに、画像データ１３ｂの属性Ａ〜Ｄに応じてＧＭＴ１３ｇ１〜１３ｇ４を選択して使用するため、適切なガマット修正を実現することができる。
【選択図】 図１３

Description

本発明は、色変換装置、色変換方法および色変換プログラムに関する。

モニターを見ながらレタッチした写真画像データや、モニターを見ながら作成したＣＧ画像データを、印刷用紙に印刷することが一般的に行われている。この場合、モニターに映し出された画像データを見たときに感じる色と、印刷用紙に印刷された画像データを見たときに感じる色とが一致していることが望ましい。人が感じる色は周囲の照明等によって変動するため、この変動を抑えるために周囲光に基づき色変換を行う画像処理装置が提案されている（例えば、特許文献１、参照。）。
かかる構成において、測光装置によって周囲光情報を取得し、その周囲光情報に基づいて色変換を行うことができるため、周囲光による色の見えの変動を抑制することが可能であった。
特開平９−２１４７８７号公報

上述した技術において周囲光の状況に応じた色変換が実現できる反面、周囲光の変動に応じて毎回色変換の演算をしなければならない。また、周辺光を考慮した色変換には数多くの行列変換が必要であるため、画像を出力するたびに膨大な演算処理を行わなければならないという問題があった。さらに、画像データの特性に適合した色修正が並行して行われるように、画像データの特性に適した色変換を行うことが望ましいが、この場合も画像データの特性ごとに膨大な演算処理を行わなければならないという問題があった。
本発明は、画像データに適合した修正を高速に行いつつ、画像の見えを正確に再現することができる色変換装置、色変換方法および色変換プログラムの提供を目的とする。

課題を解決するための手段および発明の効果

上記課題を解決するため請求項１にかかる発明では、色変換装置において元画像データが色変換画像データに変換される。上記元画像データは各画素の色が第一出力デバイスに入力可能な第一表色系で表現された画像データであり、上記色変換画像データは各画素の色が第二出力デバイスに入力可能な第二表色系で表現された画像データである。このような色変換を行うことにより、第一出力デバイスに入力されている元画像データと同等の色変換画像データを得ることができ、同色変換画像データを上記第二出力デバイスに入力し、出力させることができる。

環境パラメータ取得手段は、上記第一出力デバイスおよび上記第二出力デバイスの出力結果を観察するときの環境パラメータを取得する。第一色変換手段は、所定の変換プロファイルを利用して、上記元画像データを各画素の色が第三表色系で表現された画像データに変換する。この変換プロファイルは、上記環境パラメータを代入し、これらを変換結果に反映させることができるため、同環境パラメータを考慮した画像データを得ることができる。すなわち、上記環境パラメータによって表される観察環境のもとで、上記第一出力デバイスにて出力したときにおける上記元画像データの見えを上記第三表色系にて予測することができる。属性取得手段は、色変換の対象となる上記元画像データを取得するとともに、その属性を取得する。修正前後の対応関係を特定した複数の修正プロファイルが予め用意されており、上記修正手段はこのなかから上記元画像データの属性に適合するものを選択する。そして、選択された上記修正プロファイルに基づいて上記第一色変換手段にて変換された上記第三表色系の画像データを上記第三表色系にて修正する。

第二色変換手段は、所定の変換プロファイルを利用して、上記第一色変換手段にて得られた画像データであって、さらに上記修正手段によって修正された画像データを上記色変換画像データに変換する。すなわち、各画素の色が上記第三表色系で表現された画像データを、各画素の色が上記第二表色系で表現された上記色変換画像データに変換する。上記変換プロファイルは、上記環境パラメータが代入可能であるため、これらを変換結果に反映させることができる。従って、上記環境パラメータを考慮した色変換画像データを得ることができる。すなわち、上記環境パラメータによって表される観察環境のもとで上記色変換画像データを上記第二出力デバイスにて出力したときの見えが、上記第三表色系において上記第一色変換手段にて得られた画像データと基本的には一致するような色変換を行うことができる。従って、上記第一出力デバイスと上記第二出力デバイス間での見えを基本的に一致させることができる。ただし、上記修正手段による修正が加わっているため、基本的には見えを一致させつつも、意図的な色修正を実現させることができる。この修正は複数用意された上記修正プロファイルのうち変換対象の上記元画像データの属性に対応したものを使用して行われるため、上記元画像データに適した修正を行うことができる。

さらに、請求項２にかかる発明のように、上記修正プロファイルによって上記第一表色系による色空間の色域と上記第二表色系による色空間の色域との差に基づく修正を行うようにしてもよい。例えば、上記第一表色系による色空間の色域の方が上記第二表色系による色空間の色域よりも広い場合は、上記第二表色系による色空間の色に収まるように色を修正する上記修正プロファイルを適用すればよい。すなわち、いわゆるガマットマッピングを実現することができ、上記第一出力デバイスと上記第二出力デバイスで使用できるガマットを有効に利用することができる。

また、請求項３にかかる発明のように、上記修正プロファイルを、上記第三表色系における修正前後の色座標値が複数記述されたテーブルとして用意してもよい。このようにすることにより、演算処理の負担を軽減することができ、高速な色変換を実現することができる。

さらに、請求項４にかかる発明のように、上記属性としての上記第一表色系の種類に応じた色域の大きさに基づいて上記修正プロファイルを選択することにより、上記第一表色系による色空間の色域の大きさに応じたガマットマッピングを行うことができる。

また、請求項５にかかる発明のように、上記元画像データを画像解析することにより上記元画像データの属性を取得してもよい。すなわち、上記元画像データにて表される画像の特徴を画像解析によって抽出し、同抽出された特徴に応じた属性を特定することができる。

さらに、請求項６にかかる発明のように、上記元画像データに添付された同画像データの付帯情報に基づいて上記元画像データの属性を取得してもよい。例えば、上記画像データがＥＸＩＦ規格に準拠していれば、そのタグ情報に基づいて属性を特定することができる。ＥＸＩＦ規格のタグ情報にはデジタルカメラで撮影したことを示すデータが格納されるため、このデータに基づいて当該元画像データがＣＧ画でなく自然画であることが判別できる。

また、請求項７にかかる発明では、上記元画像データの属性が自然画であるとき修正前後で彩度の変動量が少ない上記修正プロファイルが選択される。これにより、上記修正手段による修正前後で上記元画像データの彩度を維持することができ、上記第二デバイスの出力においても自然画に含まれる多様な色彩を忠実に再現することができる。

さらに、請求項８にかかる発明では、上記元画像データの属性がＣＧ画であるとき修正前後で明度の変動量が少ない上記修正プロファイルが選択される。これにより、上記修正手段による修正前後で上記元画像データの明度を維持することができ、上記第二デバイスの出力においてもＣＧ画における明度を作成者の意図通りに再現することができる。

また、請求項９にかかる発明では、派生テーブル作成手段が備えられ、同派生テーブル作成手段によって上記第一表色系と上記第二表色系との対応関係がテーブルとして記述される。なお、上記第一表色系と上記第二表色系との対応関係は、上記環境パラメータと上記修正プロファイルを使用して上記第一色変換手段と上記修正手段と上記第二色変換手段にて変換および修正を順次行うことにより特定することができる。以上のようにして特定された色変換前後の色の対応関係を記述することにより、上記派生テーブルを作成することができる。さらに、第三色変換手段が備えられ、同第三色変換手段が上記環境パラメータ取得手段と上記属性取得手段が取得した上記環境パラメータと上記属性との組み合わせと、上記派生テーブルを作成したときの上記環境パラメータと上記修正プロファイルとの組み合わせとを比較する。そして、これらの組み合わせが一致したときには、当該派生テーブルを参照して上記元画像データを上記色変換画像データに変換することにより、上記第一色変換手段と上記第二色変換手段とを介することなく上記色変換データを得ることができ、処理負担を軽減させることができる。

また、請求項１０にかかる発明では、上記環境パラメータ取得手段がテストパターンを印刷媒体に出力する。そして、同テストパターンの視覚結果に基づいて、上記環境パラメータが取得される。これにより、観察環境を測定するための測定装置を用意しなくても済む。

むろん、以上の発明は、装置のみならず、請求項１１のような色変換方法によって実現することも可能であるし、請求項１２のように上記方法に従った処理を実行する色変換プログラムによって実現することも可能である。また、本発明にかかる装置、方法、プログラムは単独で実施される場合もあるし、ある機器に組み込まれた状態で他の装置、方法、プログラムとともに実施されることもあるなど、発明の思想としてはこれに限らず、各種の態様を含むものであり、適宜、変更可能である。

さらに、本発明のプログラムを記録した記録媒体として提供することも可能である。このプログラムの記録媒体は、磁気記録媒体であってもよいし光磁気記録媒体であってもよいし、今後開発されるいかなる記録媒体においても全く同様に考えることができる。また、一次複製品、二次複製品などの複製段階については全く問う余地無く同等である。さらに、一部がソフトウェアであって、一部がハードウェアで実現されている場合においても発明の思想において全く異なるものではなく、一部を記録媒体上に記憶しておいて必要に応じて適宜読み込まれるような形態のものとしてあってもよい。また、必ずしも全部の機能を単独のプログラムで実現するのではなく、複数のプログラムにて実現させるようなものであってもよい。この場合、各機能を複数のコンピュータに実現させるものであればよい。

ここでは、下記の順序に従って本発明の実施形態について説明する。
（１）コンピュータの構成：
（２）環境パラメータ取得処理の流れ：
（３）色変換処理の流れ：
（４）変換プロファイル（ＣＩＥＣＡＭ０２）について：
（５）まとめと変形例：

（１）コンピュータの構成：
図１は、本発明の色変換装置としてのコンピュータの概略構成を示している。同図において、コンピュータ１０には、内部バス１０ａによって接続されたＣＰＵ１１とＲＡＭ１２とＨＤＤ１３とＵＳＢインターフェイス（Ｉ／Ｆ）１４と入力機器インターフェイス（Ｉ／Ｆ）１５とビデオインターフェイス（Ｉ／Ｆ）１６とが備えられており、ＨＤＤ１３には各種プログラムデータ１３ａと複数の画像データ１３ｂとテストパターンデータ１３ｃと色変換ＬＵＴ１３ｄとガマットデータ１３ｅと環境パラメータ１３ｆとガマット修正テーブル１３ｇと分光分布データ１３ｈが記憶されている。ＣＰＵ１１は、このプログラムデータ１３ａを読み出して、同プログラムデータ１３ａに基づいた処理をＲＡＭ１２をワークエリアとして利用しながら実行する。ＵＳＢＩ／Ｆ１４にはプリンタ２０が接続されており、入力機器インターフェイス１５にはマウス４０およびキーボード５０が接続されている。さらに、ビデオＩ／Ｆ１６にはディスプレイ６０が接続されている。

図２は、コンピュータ１０にて実行されるプログラムのソフトウェア構成を示している。同図において、プリンタドライバＰが図示しないオペレーティングシステム（Ｏ／Ｓ）上にて実行されている。プリンタドライバＰは、画像データ取得部Ｐ１と環境パラメータ取得部Ｐ２と属性取得部Ｐ３と第一色変換部Ｐ４とガマット修正部Ｐ５と第二色変換部Ｐ６とハーフトーン処理部Ｐ７と印刷データ生成部Ｐ８とから構成されている。画像データ取得部Ｐ１は、印刷すべき画像データ１３ｂの指定を受け付けるとともに、指定された画像データ１３ｂをＨＤＤ１３から取得する。画像データ１３ｂは、他のアプリケーションで作成しされたり、デジタルスチルカメラ等の画像入力機器から入力されたりして予め用意されている。

画像データ１３ｂは、ドットマトリクス状に配列する複数の画素で構成されており、各画素の色がｓＲＧＢ表色系のデジタル階調によって表現されている。ｓＲＧＢ表色系の画像データ１３ｂは、ビデオＩ／Ｆ１６が入力可能な形式であり、同画像データ１３ｂを入力した画像データ１３ｂをディスプレイ６０にて出力することが可能となっている。例えば、画質調整アプリケーション等によって画像データ１３ｂの色調を調整する場合には、ユーザーはディスプレイ６０にて出力した画像データ１３ｂを見ながら色調を補正することとなる。色調を補正した画像データ１３ｂはＨＤＤ１３にて更新され、その後、プリンタ２０にて印刷するように指示されることとなる。

環境パラメータ取得部Ｐ２は、ＨＤＤ１３に記憶されたテストパターンデータ１３ｃを取得し、プリンタ２０において印刷用紙にテストパターンを印刷する。テストパターンデータ１３ｃは、プリンタ２０が印刷に使用するインク色の表色系で各画素の色が表現された画像データであり、同テストパターンデータ１３ｃをそのままハーフトーン処理等することによりプリンタ２０にてテストパターンを印刷することができる。なお、本実施形態においてプリンタ２０はインクジェット方式を採用しており、インク色としてＣ（シアン）Ｍ（マゼンタ）Ｙ（イエロー）Ｋ（ブラック）を使用している。従って、テストパターンデータ１３ｃでは各画素の色がＣＭＹＫの各階調によって表現されている。なお、本実施形態においてディスプレイ６０が本発明の第一出力デバイスに相当し、ｓＲＧＢ表色系が本発明の第一表色系に相当する。一方、プリンタ２０が本発明の第二出力デバイスに相当し、ＣＭＹＫ表色系が本発明の第二表色系に相当する。

環境パラメータ取得部Ｐ２はテストパターンを印刷すると、マウス４０やキーボード５０の入力を受け付ける。環境パラメータ取得部Ｐ２はユーザーによるテストパターンの視覚結果として照度と照明光の分光分布を取得する。環境パラメータ取得部Ｐ２は、取得した照度および照明光の分光分布を指標として、ディスプレイ６０での画像出力結果を観察する際の環境パラメータと、プリンタ２０によって印刷用紙上に出力した画像出力結果を観察する際の環境パラメータを算出する。算出した環境パラメータは、ＨＤＤに環境パラメータ１３ｆとして記憶される。

属性取得部Ｐ３は、画像データ１３ｂを取得するとともに、同画像データ１３ｂに対する空間周波数解析を行う。本実施形態においては属性取得部Ｐ３が画像データ１３ｂ１をフーリエ変換することにより、画像データ１３ｂの輝度に関する空間周波数のスペクトルを得ることができる。これにより、画像データ１３ｂにおいて輝度が変動する周期が高周波となっているか低周波となっているかを特定することができる。一般的に、自然画の画像データ１３ｂにおいては空間的に短い周期のノイズが多く含まれるため、高周波成分が多く含まれている画像データ１３ｂが自然画であると判定することができる。反対に、ＣＧ画のような画像においてはノイズが含まれないため、高周波成分が少ない画像データ１３ｂがＣＧ画であると判定することができる。ここでは、属性取得部Ｐ３が属性として自然画かＣＧ画のいずれかの画像タイプを特定する。

さらに、属性取得部Ｐ３は、画像データ１３ｂに添付されたカラープロファイルを取得する。例えば、本実施形態において画像データ１３ｂにＡｄｏｂｅＲＧＢ（以下、本明細書においてａＲＧＢと表記するものとする。ＡｄｏｂｅはＡｄｏｂｅ ｓｙｓｔｅｍｓ社の登録商標）のカラープロファイルが添付されている場合には、当該画像データ１３ｂが採用する表色系の種類がａＲＧＢ表色系であると特定する。一方、画像データ１３ｂがＲＧＢのビットマップデータであり、カラープロファイルが添付されていない場合には、当該画像データ１３ｂが採用する表色系の種類がｓＲＧＢ表色系であると特定する。すなわち、属性取得部Ｐ３は、自然画かＣＧ画のいずれかの画像タイプのほかに、ａＲＧＢかｓＲＧＢのいずれかの表色系の種類を画像データ１３ｂの属性として特定する。従って、各画像データ１３ｂが帰属する属性の組み合わせは下記の属性Ａ〜Ｄのいずれかとなる。
属性Ａ：自然画−ａＲＧＢ
属性Ｂ：自然画−ｓＲＧＢ
属性Ｃ：ＣＧ画−ａＲＧＢ
属性Ｄ：ＣＧ画−ｓＲＧＢ

第一色変換部Ｐ４は、環境パラメータ取得部Ｐ２が取得・記憶したディスプレイ６０に関する環境パラメータと画像データ取得部Ｐ１が取得した画像データ１３ｂをＨＤＤ１３から入力し、同画像データ１３ｂを各画素がＸＹＺ表色系で表現された画像データに変換し、さらに画像データ１３ｂを各画素が非機器依存色空間のＪＣｈ表色系で表現された画像データに変換する。

後者の変換においては、ＣＩＥＣＡＭ０２という変換手法を用いる。ＣＩＥＣＡＭ０２の詳細については後述するが、ＣＩＥＣＡＭ０２では環境パラメータを変数として使用することができ、環境パラメータに応じた変換結果を得ることができる。具体的には、与えられた画像データの色が、入力された環境パラメータのもとで、どのように知覚されるかをＪＣｈ表色系にて特定するための変換を行うことができる。これにより、画像データ１３ｂにおける注目物がディスプレイ６０にて実際にどのように見えているかを特定することができる。

ガマット修正部Ｐ５は、本発明の修正手段に相当し、属性取得部Ｐ３が取得した画像データ１３ｂの属性Ａ〜Ｄに基づいて参照するガマット修正テーブル（ＧＭＴ）１３ｇを選択する。なお、本実施形態において４種類の属性Ａ〜Ｄが考えられ、それぞれの属性Ａ〜Ｄに対応してＧＭＴ１３ｇ（１３ｇ１〜１３ｇ４）が予め用意されている。ガマット修正部Ｐ５は選択されたＧＭＴ１３ｇ１〜１３ｇ４を参照することにより、第一色変換部Ｐ４によって予めＪＣｈ表色系に色変換されている画像データ１３ｂをＪＣｈ表色系にて修正する。

第二色変換部Ｐ６は、環境パラメータ取得部Ｐ２が取得したプリンタ２０（印刷物）に関する環境パラメータを使用したＣＩＥＣＡＭ０２によって、ガマット修正部Ｐ５が修正を行った画像データを各画素の色がＸＹＺ表色系で表現される画像データに色変換する。ここでは、ＪＣｈ表色系の画像データをＸＹＺ表色系に色変換しているため、第一色変換部Ｐ４の逆変換を行っていることになる。なお、ＣＩＥＣＡＭ０２は可逆性を有する変換プロファイルであるため、逆変換を行うことが可能である。さらに、第二色変換部Ｐ６は、予め作成されている色変換ＬＵＴ１３ｄを参照することによりＸＹＺ表色系に色変換された画像データを各画素の色がＣＭＹＫ表色系で表現される画像データに色変換する。なお、第一色変換部Ｐ４におけるＸＹＺ値と区別するために、ここでのＸＹＺ値をＸ’Ｙ’Ｚ’値と表記するものとする。色変換ＬＵＴ１３ｄは、ＣＭＹＫ値とＸ’Ｙ’Ｚ’値との等色対応関係が予め行われた測色結果に基づいて定義されたテーブルである。

ハーフトーン処理部Ｐ７は、第二色変換部Ｐ６にて色変換されたＣＭＹＫ表色系の画像データを入力し、ディザ法や誤差拡散法等によって同画像データをハーフトーンデータに変換する。これにより、各画素において、ＣＭＹＫインクを吐出させるか吐出させないかを特定できる画像データを得ることができる。印刷データ生成部Ｐ８は、ハーフトーンデータを入力し、ラスター化等の処理を行い、プリンタ２０に出力可能な印刷データを生成する。印刷データ生成部Ｐ８が生成した画像データは、プリンタ２０に出力され、同プリンタ２０にて印刷媒体としての印刷用紙に印刷される。

（２）環境パラメータ取得処理の流れ：
図３は、環境パラメータ取得部Ｐ２が環境パラメータを取得する処理の流れを示している。ステップＳ１００にて環境パラメータを設定する旨の指示をマウス４０やキーボード５０を介して受け付ける。ステップＳ１１０においては、環境パラメータ取得部Ｐ２がテストパターンデータ１３ｃをＨＤＤ１３から取得し、同テストパターンデータ１３ｃをプリンタ２０にて印刷用紙上に出力する。テストパターンデータ１３ｃはプリンタ２０が採用するＣＭＹＫ表色系によって各画素の色が表現されているため、そのままハーフトーン処理部Ｐ７と印刷データ生成部Ｐ８にて変換を行い、プリンタ２０に出力することができる。

図４は、テストパターンの一例を示している。同図において、テストパターンは印刷媒体としての印刷用紙上に形成されており、照明光の波長分布を評価するための色度エリアＳ１と、照度を評価するための照度エリアＳ２とから構成されている。色度エリアＳ１と照度エリアＳ２はそれぞれ略矩形状に形成されており、印刷用紙上にて互いに独立して形成されている。プリンタ２０はＣＭＹＫインクを吐出可能なインクジェットプリンタであり、テストパターンは印刷用紙上にインクを被覆させることにより形成されている。

図５は、色度エリアＳ１を拡大して示している。同図において、色度エリアＳ１は、それぞれ矩形状に形成された基準カラーパッチＥ１〜Ｅ４とカラーパッチＧ１〜Ｇ４，Ｈ１〜Ｈ４，Ｉ１〜Ｉ４，Ｊ１〜Ｊ４の組から構成されている。基準カラーパッチＥ１〜Ｅ４とカラーパッチＧ１〜Ｇ４，Ｈ１〜Ｈ４，Ｉ１〜Ｉ４，Ｊ１〜Ｊ４では、それぞれ略無彩色が再現されており、１〜４行目になるほど明度（Ｌ^*＝８０，６０，４０，２０相当）が低くなるようなグラデーションパターンとなっている。基準カラーパッチＥ１〜Ｅ４は、印刷用紙上にＫインクのみを吐出させることにより形成されている。基準カラーパッチＥ１〜Ｅ４においては、Ｋインクによる被覆率を徐々に高くしていくことにより、次第に明度が低くなるグラデーションパターンが再現されている。

カラーパッチＧ１〜Ｇ４，Ｈ１〜Ｈ４，Ｉ１〜Ｉ４，Ｊ１〜Ｊ４は、Ｋインクによるドットの他に有彩色のＣＭＹインクによるドットも使用して形成された、いわゆる有彩色インクによるコンポジットグレーとされている。ただし、カラーパッチＧ１〜Ｇ４，Ｈ１〜Ｈ４，Ｉ１〜Ｉ４，Ｊ１〜Ｊ４は、無彩色に近い色ではあるが、無彩色となるＣＭＹインクのドット構成比率から意図的にずらしたＣＭＹインクのドット構成比率で形成されており、厳密な無彩色ではない。ただし、カラーパッチＧ１〜Ｇ４，Ｈ１〜Ｈ４，Ｉ１〜Ｉ４，Ｊ１〜Ｊ４は観察する照明光によっては、無彩色に知覚されるように各色インクドットの構成比率が設定されている。

例えば、各明度に対応するカラーパッチＧ１〜Ｇ４の組は太陽光（ＣＩＥ−Ｄ５０光）の下で無彩色に知覚されるように各色インクドットの構成比率が設定されており、各明度に対応するカラーパッチＨ１〜Ｈ４の組は白熱灯（ＣＩＥ−Ａ光）の下で無彩色に知覚されるように各色インクドットの構成比率が設定されている。また、各明度に対応するカラーパッチＩ１〜Ｉ４の組はＣＩＥ−Ｄ６０光の下で無彩色に知覚されるように各色インクドットの構成比率が設定されており、各明度に対応するカラーパッチＪ１〜Ｊ４の組はＦ２光の下で無彩色に知覚されるように各色インクドットの構成比率が設定されている。

実際に印刷物を観察する照明光の下で、色度エリアＳ１を観察したユーザーは、各明度において最も基準カラーパッチＥ１〜Ｅ４と似た色に知覚されるカラーパッチＧ１〜Ｇ４，Ｈ１〜Ｈ４，Ｉ１〜Ｉ４，Ｊ１〜Ｊ４をマウス４０やキーボード５０によって選択し、この選択を環境パラメータ取得部Ｐ２がステップＳ１２０にて受け付ける。これにより、テストパターンがどのような照明光の下で観察されているかを判断することができる。環境パラメータ取得部Ｐ２は各照明光の分光分布を、ＨＤＤ１３に記憶された分光分布データ１３ｈから取得する。

図６は、照明光（Ｄ５０光，Ａ光）の分光分布と等色関数をグラフに示している。同図に示すように、各波長λのエネルギー分布（分光分布）が照明光（Ｄ５０光，Ａ光）によって大きく異なっていることが分かる。Ｄ５０光には各波長のエネルギーが均等に含まれ、Ａ光には長波長のエネルギーが多く含まれている。このような分光分布は、照明光ごとに分光分布データ１３ｈに格納されている。分光分布データ１３ｈにおいては各照明光の分光分布の値がテーブルに記述されていてもよいし、対応関係が関数等によって定義されていてもよい。各照明光の分光分布はＣＩＥ規格等に定められているため、入手することができる。一方、等色関数ｘ（λ），ｙ（λ），ｚ（λ）は、波長毎に赤・緑・青の３刺激をどれくらい感じるかを示す関数であり、実験・経験的に得られている。また、照明光を観察物に反射させた場合に知覚される色の３刺激値ＸＹＺは下記式（１）によって表すことができる。

なお、上記式（１）においてＲ（λ）は観察物の分光反射率であり、Ｐ（λ）は照明光の分光分布を示している。ステップＳ１２０にて照明光が特定できると、ステップＳ１３０ではその分光分布Ｐ（λ）を取得し、上記式（３）に代入することにより、その照明光における白色点の色を３刺激値のＸＹＺ値として算出する。なお、本実施形態においては白色点の色は完全白色板における色を想定しており、同完全白色板は全波長領域において分光反射率Ｒ（λ）＝１となる。全可視光領域において分光分布Ｐ（λ）と分光反射率Ｒ（λ）とが与えられるため、上記式（１）にて白色点の色を算出することができる。

白色点の色は本発明における照明光の分光分布を指標とした環境パラメータに相当し、それぞれ３刺激値をＸ_W，Ｙ_W，Ｚ_Wと表記するものとする。なお、ディスプレイ６０上の白色点の色をＸ_WＹ_WＺ_Wと表すものとし、印刷用紙上の白色点の色をＸ_WＹ_WＺ_Wと表すものとする。また、照明光と印刷用紙上の白色点の色Ｘ_WＹ_WＺ_Wとの間には一義的な対応関係があるため、この対応関係を規定したテーブルを記憶しておき、環境パラメータ取得部Ｐ２が同テーブルを参照して照明光に対応する白色点の色Ｘ_WＹ_WＺ_Wを取得するようにしてもよい。さらに、完全白色板の分光反射率Ｒ（λ）を全波長領域において１としたが、印刷媒体に応じて白色点の分光反射率Ｒ（λ）を設定してもよい。ところで、上記式（１）では、分光分布Ｐ（λ）と分光反射率Ｒ（λ）とを相乗しており、ＸＹＺ値が観察物の分光反射率Ｒ（λ）に依存することが分かる。

図７は、Ｋインクドットのみを使用した所定明度のグレーと、ＣＭＹインクドットも併用した所定明度のコンポジットグレーの分光反射率Ｒ（λ）を示している。Ｋインクドットによるグレーにおいては分光反射率Ｒ（λ）が各波長領域においてほぼ均一である。一方、コンポジットグレーは、個々の分光反射率が不均一であるＣＭＹインクの合成であるため、全体としての分光反射率Ｒ（λ）が不均一となっている。このように、分光反射率Ｒ（λ）が不均一なコンポジットグレーにおいては、上記式（１）によるＸＹＺ値は照明光の差によって大きく変動することとなる。一方、各波長領域において分光反射率Ｒ（λ）がほぼ均一なＫインクドットのグレーにおいては照明光の差によるＸＹＺ値の変動は小さい。このことは、Ｋインクドットのグレーは照明光の変動に応じた知覚色の変動が少なく、常に色味を帯びることがなく無彩色に感じられることを意味する。

上述したとおり基準カラーパッチＥ１〜Ｅ４はＫインクドットのみを使用して再現されており、基準カラーパッチＥ１〜Ｅ４の分光反射率Ｒ（λ）が各波長領域においてほぼ均一であるということができる。従って、コンポジットグレーのカラーパッチＦ１〜Ｆ４，Ｇ１〜Ｇ４，Ｈ１〜Ｈ４，Ｉ１〜Ｉ４は照明光に応じて色味を帯びるのに対して、基準カラーパッチＥ１〜Ｅ４は常に無彩色に感じられることとなる。従って、基準カラーパッチＥ１〜Ｅ４と最も近い色に知覚されるカラーパッチＦ１〜Ｆ４，Ｇ１〜Ｇ４，Ｈ１〜Ｈ４，Ｉ１〜Ｉ４の組を選択することにより、その照明光において最も無彩色に近い色に知覚されるカラーパッチＦ１〜Ｆ４，Ｇ１〜Ｇ４，Ｈ１〜Ｈ４，Ｉ１〜Ｉ４の組を容易に選択することができる。

以上のように無彩色からの変動が少ない基準カラーパッチＥ１〜Ｅ４をカラーパッチＦ１〜Ｆ４，Ｇ１〜Ｇ４，Ｈ１〜Ｈ４，Ｉ１〜Ｉ４と各明度において並列させておくことにより、ユーザーは正確に無彩色に近いカラーパッチＦ１〜Ｆ４，Ｇ１〜Ｇ４，Ｈ１〜Ｈ４，Ｉ１〜Ｉ４の組を選択することができる。また、基準カラーパッチＥ１〜Ｅ４とカラーパッチＦ１〜Ｆ４，Ｇ１〜Ｇ４，Ｈ１〜Ｈ４，Ｉ１〜Ｉ４を明度のグラデーションとすることにより、いずれの明度においてもバランスよく無彩色に見えるカラーパッチＦ１〜Ｆ４，Ｇ１〜Ｇ４，Ｈ１〜Ｈ４，Ｉ１〜Ｉ４の組を選択することができる。これにより、環境パラメータ取得部Ｐ２は正確に現在の照明光の分光分布Ｐ（λ）を取得することができ、正確に印刷用紙上の白色点の色Ｘ_WＹ_WＺ_Wを算出することができる。

ただし、無彩色は人間が最もよく記憶している記憶色であり、人間は無彩色に付近における色差には敏感である。従って、ユーザーは自己の記憶に基づく絶対的な基準によって無彩色に近いカラーパッチＦ１〜Ｆ４，Ｇ１〜Ｇ４，Ｈ１〜Ｈ４，Ｉ１〜Ｉ４を選択することも可能である。そのため、必ずしも基準カラーパッチＥ１〜Ｅ４が形成されている必要はなく、基準カラーパッチＥ１〜Ｅ４を省略してもよい。従って、例えばＫインクを使用しないプリンタにおいても本発明を適用することは可能である。なお、以上において算出した印刷用紙上の白色点の色Ｘ_WＹ_WＺ_Wは、実際に印刷物を観察する照明光の下でテストパターンを観察した結果得られたものであり、ディスプレイ６０上の白色点の色Ｘ_WＹ_WＺ_Wは別途取得しておく必要がある。

ただし、ディスプレイ６０は自発光デバイスであるため、ディスプレイ６０に関する白色点の色Ｘ_WＹ_WＺ_Wは照明光に依存することなく、ディスプレイの仕様から得ることができる。すなわち、ディスプレイ６０に関しては、ディスプレイ６０の観察環境下における照明光による白色点の色Ｘ_WＹ_WＺ_Wの依存度が０であるため、ディスプレイ６０が発光する仕様上の白色の色度そのものをディスプレイの観察環境下における白色点の色Ｘ_WＹ_WＺ_Wとすることができる。従って、ステップＳ１２０にてディスプレイ６０の白色仕様の指定入力を受け付けることにより、環境パラメータ取得部Ｐ２がディスプレイの観察環境下における白色点の色Ｘ_WＹ_WＺ_Wを取得することができる。以上のようにして環境パラメータ取得部Ｐ２が各環境下における白色点の色度を示す環境パラメータＸ_WＹ_WＺ_Wを取得すると、次にパラメータ取得部Ｐ２は照度に関する環境パラメータを取得する。

図８は、照度エリアＳ２を拡大して示している。照度エリアＳ２はＡ〜Ｄ行×１〜５列のパッチＡ１，Ａ２・・，Ｂ１，Ｂ２・・，Ｃ１，Ｃ２・・，Ｃ１，Ｃ２・・，Ｄ１，Ｄ２・・から構成されている。Ａ行のパッチＡ１，Ａ２・・は、それぞれ略矩形状の黒帯と白帯を交互に配列させた模様となっている。Ａ行のパッチＡ１，Ａ２・・において黒帯と白帯が配列する空間周波数は、Ａ１＞Ａ２＞Ａ３＞Ａ４＞Ａ５となっている。Ｂ行のパッチＢ１，Ｂ２・・は、それぞれ略矩形状の黒帯と灰帯を交互に配列させた模様となっており、交互に配列する帯の空間周波数はＢ１＞Ｂ２＞Ｂ３＞Ｂ４＞Ｂ５となっている。

Ｃ行，Ｄ行についても同様の空間周波数とされており、Ｃ行，Ｄ行については順に灰帯が濃くなっている。すなわち、Ａ〜Ｄ行のいずれにおいても異色の帯が交互に配列する模様となっており、その明度コントラストがＡ＞Ｂ＞Ｃ＞Ｄとなっている。同図の下段においては、各パッチＡ１，Ａ２・・・の濃度の変動を示しており、各パッチＡ１，Ａ２・・・の濃度が矩形波状に変動し、空間周波数が徐々に変化させられていることが分かる。なお、同図においては図示の都合上、各帯間に明度コントラストを付けたテストパターンを例示したが、各帯間に色相や彩度等のコントラストを有するテストパターンを用意してもよい。

図９は、一定の照度における人間の視覚の空間周波数特性をグラフにして示している。同図において、縦軸に知覚レベルを示しており、横軸に空間周波数[cycle/mm]を示している。同図から空間周波数が大きいほど知覚レベルが低下することが分かる。図５に示した照度エリアＳ２では、各パッチの黒帯と白帯が配列する空間周波数は、１＞２＞３＞４＞５列目となっているため、１＞２＞３＞４＞５列目の順に知覚しづらいこととなる。また、隣接する異色の明度を異ならせた場合と、色相（黄−青，赤−緑）を異ならせた場合とでは、明度を異ならせた場合の方が知覚レベルの空間周波数依存性が高く、色相差は明度差ほど知覚できないことが分かる。図５に示した照度エリアＳ２では、隣接する帯間で明度コントラストを異ならせており、その明度差がＡ＞Ｂ＞Ｃ＞Ｄ行目となっているため、Ａ＞Ｂ＞Ｃ＞Ｄの順に境界が知覚しやすくなっている。このように、一定の照度においては、視覚の空間周波数特性に図９のような傾向が見られるが、照度が変動すると視覚の空間周波数も変動することが分かっている。

図１０は、実験によって得られた照度と分解限界周波数との関係を示している。また、分解限界周波数とは、与えられた条件下で交互に配列する帯の境界を認識できる最大の空間周波数を意味し、具体的には図８の照度エリアＳ２と同様にいくつかの配列周波数を有するパッチを用意しておき、各照度において帯間の色の差異が視覚できる限界のパッチの配列周波数を分解限界周波数として特定する。なお、同図において、縦軸が分解限界周波数を示し、横軸は照度を示している。同図において、Ａ〜Ｄ行のパッチと同様の明度コントラストとなる帯についての各照度における分解限界周波数の傾向をそれぞれ線で結んでいる。また、各パッチＡ１，Ａ２・・・の帯が配列する周波数と、実験で得られた分解限界周波数とが一致する点をプロットしている。いずれの明度コントラストにおいても、照度が低下すると、分解限界周波数が低下することが分かる。すなわち、いずれの明度コントラストにおいても、分解限界周波数を特定することにより、一義的に照度を特定することができる。

図８の照度エリアＳ２を観察し、帯間の色の差を視覚できる限界のパッチＡ１，Ａ２・・・を特定することにより、そのときの照度を特定することができる。例えば、パッチＡ２の帯間の色の差が視認でき、パッチＡ１の帯間の色の差が視認できない場合には、分解限界周波数がパッチＡ２の配列周波数となり、照度がパッチＡ２に対応するＩ_Aであると特定することができる。この照度において、Ｂ行のパッチＢ１〜Ｂ５を観察すると、パッチＢ４の帯間の色の差が視認でき、パッチＢ３の帯間の色の差が視認できず、パッチＢ４に対応する照度Ｉ_Bを特定することができる。さらに、Ｃ行のパッチＣ１〜Ｃ５を観察すると、パッチＣ５の帯間の色の差が視認でき、パッチＣ４の帯間の色の差が視認できず、パッチＣ５に対応する照度Ｉ_Cを特定することができる。

従って、ステップ１４０にて視覚できた限界のパッチＡ１，Ａ２・・・の指定を受け付けることにより、環境パラメータ取得部Ｐ２は、照度Ｉ_A，Ｉ_B，Ｉ_Cを特定することができる。環境パラメータ取得部Ｐ２は図１０の対応関係をテーブルとして記憶しており、同テーブルを参照することにより、照度Ｉ_A，Ｉ_B，Ｉ_Cを特定することができる。各パッチＡ１，Ａ２・・・の周辺に対応する照度Ｉ_A，Ｉ_B，Ｉ_Cを示す文字を印刷しておき、ユーザーから知覚できたパッチＡ１，Ａ２・・・に対応する照度Ｉ_A，Ｉ_B，Ｉ_Cの入力を受け付けるようにしてもよい。

また、図８においては、パッチＡ１，Ａ２・・・の帯間に明度コントラストを設けるテストパターンを例示したが、色相や彩度のコントラストを各パッチの帯間に設けるようにしてもよい。例えば、Ｌ^*ａ^*ｂ^*表色系においてａ^*ｂ^*成分を変動させてもよい。図９に示すように、色相コントラストは明度コントラストよりも知覚されにくい。従って、コントラストが知覚されやすい高照度環境において微妙な照度を特定するために色相コントラストのテストパターンを使用してもよい。さらに、図８の例では、明度が矩形波状に変動するものを例示したが、明度、色相、彩度をサイン波状に変動させたり、のこぎり波状に変動させたりしてもよい。この場合、色の変動が緩やかとなるため、コントラストが知覚されにくくなる。従って、高照度環境において微妙な照度を特定するために好適であるということができる。

ステップＳ１４０では、ディスプレイ６０を観察する周囲照明の環境下でテストパターンを視認した結果が入力される。これにより、ディスプレイ６０を観察する差異の照度Ｉ_A，Ｉ_B，Ｉ_Cを特定することができる。ステップＳ１５０においては、ステップＳ１４０にて特定した照度Ｉ_A，Ｉ_B，Ｉ_Cの平均値を照度Ｉとして算出する。例えば、下記式（２）のように、相加平均によって照度Ｉを算出する。
Ｉ＝（Ｉ_A，Ｉ_B，Ｉ_C）／３ ・・・（２）

Ａ〜Ｄ行のいずれかについて視覚できる限界のパッチを特定すれば照度を得ることができるが、各明度コントラストにおいてそれぞれ照度Ｉ_A，Ｉ_B，Ｉ_Cを特定し、平均を照度Ｉとして特定することにより、照度Ｉの精度を向上させることができる。以上のようにして照度Iが特定できると、ステップS１６０にて環境パラメータ取得部Ｐ２が照度Iを指標として、環境パラメータを算出する。以下、ステップＳ１６０にて算出される環境パラメータについて説明する。
まず、下記式（３）によってディスプレイ観察環境における白色点の輝度Ｌ_sw［cd/m²］を算出する。
Ｌ_sw＝I／π ・・・（３）

例えば、照度I＝５００［lux］であった場合には、白色点輝度Ｌ_sw≒１５９［cd/m²］であると算出される。次に、ディスプレイ６０が白色を表示させたときのデバイス輝度Ｌ_dwを取得する。ディスプレイ６０は自発光デバイスであるため、デバイス輝度Ｌ_dwはディスプレイ６０の仕様に依存する。従って、ディスプレイ６０のデバイス輝度Ｌ_dwはディスプレイ６０の仕様書等から得ることができる。ここでは、ディスプレイ６０の仕様上のデバイス輝度Ｌ_dwが８０［cd/m²］であったものとして説明する。次に、順応輝度Ｌ_Aを下記式（４）によって算出する。なお、順応輝度Ｌ_Aは本発明の色変換プロファイルに使用される環境パラメータの１つを構成する。
Ｌ_A＝Ｌ_dw／５ ・・・（４）

白色輝度が８０［cd/m²］であった場合には、順応輝度Ｌ_A＝１６［cd/m²］となる。さらに、下記式（５）によって周囲比Ｓｒを算出する。
Ｓｒ＝Ｌ_sw／Ｌ_dw ・・・（５）
上記式（５）において、周囲比Ｓｒはディスプレイ観察環境における完全白色板輝度Ｌ_swとディスプレイ６０の仕様上の白色輝度Ｌ_dwとの比として算出される。上記の例では、周囲比Ｓｒ≒２となる。以上のようにして周囲比Ｓｒを算出すると、周囲比Ｓｒの値に応じて環境パラメータｃ，Ｎ_c，Ｆを特定する。

図１１は、周囲比Ｓｒの値と、環境パラメータｃ，Ｎ_c，Ｆの値との対応関係を示すテーブルである。周囲比Ｓｒの値に応じて、観察環境が、暗黒の周囲環境と、薄暗い周囲環境と、平均的な周囲環境とに分類され、各周囲環境に応じてｃ，Ｎ_c，Ｆが決められている。周囲比Ｓｒを算出した環境パラメータ取得部Ｐ２は、このテーブルを参照して、環境パラメータｃ，Ｎ_c，Ｆを特定する。上記の例では、周囲比Ｓｒ≒２となっているため、環境パラメータｃ，Ｎ_c，Ｆは、それぞれ平均的な周囲環境の０．６９，１．０，１．０であると特定される。例えば、同一仕様のディスプレイ６０においてテストパターンの観察結果によって得られる照度I＝３８［lux］であった場合には、周囲比Ｓｒ≒０．１５となり、環境パラメータｃ，Ｎ_c，Ｆは、それぞれ薄暗い周囲環境の０．５９，０．９，０．９であると特定される。

以上の手順によって環境パラメータ取得部Ｐ２が環境パラメータとしてのＬ_A，ｃ，Ｎ_c，Ｆを取得すると、ステップＳ１７０にて印刷結果を観察する環境下でのテストパターンの観察結果の入力を受け付ける。すなわち、ステップＳ１４０ではディスプレイ６０の観察環境におけるテストパターンの観察結果を受け付けるのに対して、ステップＳ１７０ではプリンタ２０にて印刷した印刷物の観察環境におけるテストパターンの観察結果を受け付ける。例えば、プリンタ２０にて店頭ＰＯＰ用のポスターを印刷したい場合には、店頭の照明環境のもとでステップＳ１１０にて印刷したテストパターンの観察を行い、その観察結果をステップＳ１７０では受け付ける。具体的には、ステップＳ１４０と同様に、視覚できた限界のパッチＡ１，Ａ２・・・の指定を受け付ける。

ステップＳ１８０においては、ステップＳ１５０と同様に上記式（２）によって平均の照度Ｉを算出する。さらに、ステップＳ１９０においては、ステップＳ１６０と同様に上記式（３）〜（５）によって印刷物の観察環境に関しての環境パラメータＬ_A，ｃ，Ｎ_c，Ｆを取得する。ただし、印刷物はディスプレイ６０のような自発光デバイスではないため、デバイス輝度Ｌ_dwは印刷物の観察環境における白色点輝度Ｌ_swと一致すると考えることができる。印刷用紙は完全白色板であると考えることができ、印刷物にて白色を表現する場合には、印刷用紙には何ら色材が付着されないからである。

例えば、ステップＳ１７０にて取得した照度Ｉが１０００［lux］であった場合には、Ｌ_sw＝Ｌ_dw＝３１８，Ｌ_A＝６３．７，Ｓｒ＝１．０となり、図１１のテーブルからｃ，Ｎ_c，Ｆはそれぞれ０．６９，１．０，１．０であると特定される。印刷物においては常にＬ_sw＝Ｌ_dwが成り立つため、ｃ，Ｎ_c，Ｆも常に０．６９，１．０，１．０となる。このようにステップＳ１４０，Ｓ１７０を行うことにより、ディスプレイ６０の観察環境に関する環境パラメータＬ_A，ｃ，Ｎ_c，Ｆと、プリンタ２０（印刷用紙）の観察環境に関する環境パラメータＬ_A，ｃ，Ｎ_c，Ｆをそれぞれ環境パラメータ取得部Ｐ２が取得することができる。環境パラメータ取得部Ｐ２は、ステップＳ１９５にてディスプレイ６０の観察環境に関する環境パラメータ（Ｘ_WＹ_WＺ_W，Ｌ_A，ｃ，Ｎ_c，Ｆ）と、印刷物の観察環境に関する環境パラメータ（Ｘ_WＹ_WＺ_W，Ｌ_A，ｃ，Ｎ_c，Ｆ）をそれぞれＨＤＤ１３に環境パラメータ１３ｆとして更新記憶させる。なお、環境パラメータＹ_Wはそれぞれディスプレイ６０と印刷用紙の白色点の輝度をそれぞれ意味している。

（３）色変換処理の流れ：
図１２は色変換処理の流れをフローチャートにより示し、図１３は色変換処理の流れを模式的に示している。ステップＳ２００においては、画像データ取得部Ｐ１がＨＤＤ１３から画像データ１３ｂを取得する。この時点で画像データ１３ｂは、各画素の色が本発明の第一表色系のｓＲＧＢ表色系またはａＲＧＢ表色系で表現されており、元画像データに相当する。ディスプレイ６０は、ｓＲＧＢ表色系の画像データ１３ｂを入力することが可能であり、ｓＲＧＢ表色系の画像データ１３ｂに基づく表示を行っている。

ステップＳ２１０では属性取得部Ｐ３が画像データ１３ｂを取得し、同画像データ１３ｂに対してフーリエ変換を行う。画像データ１３ｂの各画素のＲＧＢに基づいて各画素の輝度を算出し、この輝度について空間に関するフーリエ変換を行う。なお、フーリエ変換における空間軸方向は画像データ１３ｂの縦方向または横方向のいずれかにしてもよいし、２次元フーリエ変換を行ってもよい。このようにフーリエ変換を行うことにより、各空間周波数のごとの輝度強度に応じたスペクトルを得ることができる。従って、画像データ１３ｂにおいてどの空間周波数成分が多く含まれるかを特定することができる。

ステップＳ２２０においては、画像データ１３ｂが自然画であるかＣＧ画であるかを判定する。各空間周波数のごとの輝度強度に応じたスペクトルが得られており、このスペクトルが高周波側に偏っていれば自然画であると判定する。例えば、ある平均的な空間周波数を設定しておき、この空間周波数よりも高周波側のスペクトル強度の積算値と、この空間周波数よりも低周波側のスペクトル強度の積算値とを比較し、前者の方が大きければ自然画であると判定するようにしてもよい。なお、自然画でない場合には、ＣＧ画の画像データ１３ｂであると判定する。すなわち、ＣＧ画においては、基本的に高周波ノイズが含まれることがないため、高周波のスペクトルが少ないことをもってＣＧ画であると特定することができる。

ステップＳ２２０において画像データ１３ｂの画像タイプが特定されると、ステップＳ２３０においては画像データ１３ｂの表色系を特定する。本実施形態において、画像データ１３ｂがａＲＧＢ表色系であるかｓＲＧＢ表色系であるかを、画像データ１３ｂに添付されたａＲＧＢカラープロファイルの有無によって特定する。ａＲＧＢカラープロファイルは、ＲＧＢで表現される色座標と非機器依存色空間における色座標とを対応づけるプロファイルデータであり、画像データ１３ｂにおけるＲＧＢの色座標値が絶対的にどの色を意味しているかを特定することができる。従って、画像データ１３ｂにａＲＧＢカラープロファイルを添付することにより、非機器依存色空間において画像データ１３ｂの色が存在し得る色域も特定することができる。

ａＲＧＢカラープロファイルが添付されている画像データ１３ｂが採用する表色系の種類はａＲＧＢ表色系となる。一方、カラープロファイルが添付されていない画像データ１３ｂが採用する表色系の種類はｓＲＧＢ表色系であると特定される。ｓＲＧＢ表色系は汎用的な表色系であり、特に指定がない場合には、ｓＲＧＢ表色系であるとして画像データ１３ｂが取り扱われる。以上のようにして画像データ１３ｂの属性である画像タイプと表色系とが特定されると、ステップＳ２４０にて、これらの組み合わせが上述した属性Ａ〜Ｄのいずれに属するかを特定する。

ステップＳ２５０においては、画像データ１３ｂから一つの画素を選択し、同画素のＲＧＢ階調を取得する。ステップＳ２６０においては、画像データ１３ｂに添付されたカラープロファイルを参照したり、公知の等色変換式等により、ｓＲＧＢやａＲＧＢ表色系で表された当該画素の色を、ＸＹＺ表色系に変換する。ステップＳ２７０においては、第一色変換部Ｐ４がＸＹＺ値を取得するとともに、環境パラメータ取得部Ｐ２を介してディスプレイ６０に関する環境パラメータＸ_WＹ_WＺ_W，Ｌ_A，ｃ，Ｎ_c，Ｆを取得する。上述したとおりディスプレイ６０に関する環境パラメータＸ_WＹ_WＺ_W，Ｌ_A，ｃ，Ｎ_c，Ｆは、予めＨＤＤ１３の環境パラメータ１３ｆに格納されているため、ＣＰＵ１１がこれを読み出して使用することができる。さらに、第一色変換部Ｐ４が、ディスプレイ６０に関する環境パラメータＸ_WＹ_WＺ_W，Ｌ_A，ｃ，Ｎ_c，ＦをＣＩＥＣＡＭ０２に代入しつつ、ＣＩＥＣＡＭ０２を用いてＸＹＺ値をＪＣｈ値に変換する。

ＣＩＥＣＡＭ０２は可逆変換プロファイルであり、順方向の変換においてＸＹＺ値をＪＣｈ値に変換することができる。その際に、環境パラメータＸ_WＹ_WＺ_W，Ｌ_A，ｃ，Ｎ_c，Ｆを変換式に代入することが可能であり、ディスプレイ６０の観察環境に応じた変換結果をＪＣｈ表色系にて得ることができる。ＣＩＥＣＡＭ０２モデルによれば、環境パラメータＸ_WＹ_WＺ_W，Ｌ_A，ｃ，Ｎ_c，Ｆで表されるディスプレイ６０の観察環境下において、ｓＲＧＢ値やａＲＧＢ値に基づいてディスプレイ６０が表示した色が、どのような色で人間に知覚されるかをＪＣｈ表色系にて特定することができる。なお、ＪＣｈ表色系は本発明の第三表色系に相当し、ステップＳ２６０〜Ｓ２７０が本発明の第一色変換工程に相当することとなる。

以上のようにしてＪＣｈ値を算出すると、ステップＳ２８０にてガマット修正が行われる。ステップＳ２８０では、ガマット修正部Ｐ５が予めステップＳ２４０にて特定しておいた画像データ１３ｂの属性Ａ〜Ｄを取得し、同属性Ａ〜Ｄに適合するＧＭＴ１３ｇ１〜１３ｇ４を選択する。そして、ステップＳ２８０においてガマット修正部Ｐ５は選択されたＧＭＴ１３ｇ１〜１３ｇ４を参照することにより、上述した第一色変換工程にて変換して得られたＪＣｈ表色系の画像データ１３ｂを修正する。

図１４はＧＭＴ１３ｇ１〜１３ｇ４を示している。同図において、各ＧＭＴ１３ｇ１〜１３ｇ４においては、それぞれ対応する属性Ａ〜Ｄが関連づけられており。各属性Ａ〜Ｄに対応するＧＭＴ１３ｇ１〜１３ｇ４を識別することが可能となっている。例えば、属性Ａの画像データ１３ｂについてはＧＭＴ１３ｇ１が選択される。このようにすることにより、当該画像データ１３ｂに適合した修正を行うことが可能となる。各ＧＭＴ１３ｇ１〜１３ｇ４においては、修正前のＪＣｈ値と修正後のＪＣｈ値（以下、Ｊ’Ｃ’ｈ’値と表記するものとする。）との対応関係が記述されている。本実施形態において、ＪＣｈ値はそれぞれ２５６階調であるが、一部の代表的なＪＣｈ値の組み合わせについて、対応するＪ’Ｃ’ｈ’値の組み合わせが記述されている。

図１５は各ＧＭＴ１３ｇ１〜１３ｇ４による修正態様を示している。同図においては、ＪＣｈ表色空間を所定の色相ｈ軸にてスライスした断面を示しており、縦軸に明度Ｊが示され、横軸に彩度Ｃが示されている。修正態様は修正前後のＪＣｈ−Ｊ’Ｃ’ｈ’値を結んだベクトルｖによって示されており、同ベクトルｖの長さと方向が修正量および修正方向に対応している。また、同図においては第一表色系であるａＲＧＢ色空間とｓＲＧＢ色空間の色域の最外郭、および、プリンタ２０が再現可能な色域の最外郭が波線によって示されている。各色空間の設計や各出力デバイスのハード上の制限上、各色域には相違が見られる。また、それぞれの色域は後述するＣＩＥＣＡＭ０２における環境パラメータＸ_WＹ_WＺ_W，Ｌ_A，ｃ，Ｎ_c，Ｆの変動に応じて変動する。図１５においては、一般的な観察環境における環境パラメータＸ_WＹ_WＺ_W，Ｌ_A，ｃ，Ｎ_c，Ｆの変動に起因した各色域の最外郭の変動幅をハッチングによって示している。

属性Ａに対応するＧＭＴ１３ｇ１において、ベクトルｖは明度Ｊ軸に平行になるような方向とされており、彩度Ｃが修正前後で変化しないようになっている。属性Ａの画像タイプは、多様な色彩が存在しうる自然画であり、極力、彩度Ｃを損なわないようにＧＭＴ１３ｇ１が作成されている。ＧＭＴ１３ｇ１において、プリンタ２０が再現可能な色域よりも修正後のＪ’Ｃ’ｈ’値が内側となるようにベクトルｖの長さが設定されている。ａＲＧＢ表色空間は比較的広い色域を有しているため、ベクトルｖを長めに設定しておく必要がある。また、上述したとおり色域の最外郭は、環境パラメータＸ_WＹ_WＺ_W，Ｌ_A，ｃ，Ｎ_c，Ｆの変動に起因して変動するが、その変動幅を見越してプリンタ２０が再現可能な色域よりも修正後のＪ’Ｃ’ｈ’値が内側となるようにベクトルｖの長さが設定されている。

かかるＧＭＴ１３ｇ１を使用して、属性Ａの画像データ１３ｂを修正することにより、自然画の色彩を維持しつつ、ａＲＧＢ表色空間の色域とプリンタ２０が再現可能な色域との不整合を緩和するガマットマッピングを行うことができる。これにより、基本的には等色性を維持しつつも、プリンタ２０が再現できるような画像データ１３ｂを得ることができる。なお、ＧＭＴ１３ｇ１〜１３ｇ４を作成するにあたり、予めＪＣｈ表色系においてａＲＧＢ表色系の色域とｓＲＧＢ表色系の色域とプリンタ２０が再現可能な色域とを予め特定しておく必要がある。これらの色域の最外郭の色を非機器依存色空間にて特定したデータを用意しておけば、後述するＣＩＥＣＡＭ０２変換によって環境パラメータＸ_WＹ_WＺ_W，Ｌ_A，ｃ，Ｎ_c，Ｆを考慮した色域を得ることができる。例えば、図２に示すように各色域の最外郭を非機器依存色空間にて特定したガマットデータ１３ｅを用意しておき、これを第一色変換部Ｐ３がＣＩＥＣＡＭ０２の順変換を行うことにより、観察環境における色域を特定することができる。ａＲＧＢ表色系の色域とｓＲＧＢ表色系の色域については、これらのＲＧＢ値と非機器依存空間における色彩値との対応を記述したカラープロファイルが用意されているため、これらのガマットデータ１３ｅを得ることができる。プリンタ２０に関しては、実際にカラーパッチを印刷して測色を行うことにより、ガマットデータ１３ｅを得ることができる。

一方、属性Ｂに対応するＧＭＴ１３ｇ２においても、ベクトルｖは明度Ｊ軸に平行になるような方向となっている。属性Ｂの画像タイプも自然画であり、極力、彩度Ｃを損なわないようにする必要があるからである。ＧＭＴ１３ｇ２においても、プリンタ２０が再現可能な色域よりも修正後のＪ’Ｃ’ｈ’値が内側となるようにベクトルｖの長さが設定されているが、ＧＭＴ１３ｇ１よりもベクトルｖの長さが短く設定されている。これは、ａＲＧＢ表色空間の色域よりもｓＲＧＢ表色空間の色域の方がプリンタ２０の再現可能な色域に近い大きさであり、修正量も少なくて済むからである。これによって、過剰に修正することが防止でき、プリンタ２０が再現可能な色域を無駄なく使用させることができる。

属性Ｃに対応するＧＭＴ１３ｇ３において、ベクトルｖは彩度Ｃ軸に平行になるような方向とされており、明度Ｊが修正前後で変化しないようになっている。属性Ｂの画像タイプは、ＣＧ画である。ＣＧ画においては明度のグラデーションが多用されるため、極力、明度Ｊを損なわないようにＧＭＴ１３ｇ３が作成されている。ＧＭＴ１３ｇ３においても、ａＲＧＢ表色空間の広範な色域とプリンタ２０が再現可能な色域との不整合を解消するため、ベクトルｖが長めに設定されている。さらに、属性Ｄに対応するＧＭＴ１３ｇ４においても、ＧＭＴ１３ｇ３と同様にベクトルｖは彩度Ｃ軸に平行になるような方向とされており、明度Ｊが修正前後で変化しないようになっている。ＧＭＴ１３ｇ４においても、ＧＭＴ１３ｇ２と同様にベクトルｖが短めに設定されている。

このように、属性Ａ〜ＤごとにＧＭＴ１３ｇ１〜１３ｇ４を用意しておくことにより、入力された元の画像データ１３ｂの画像タイプや表色系に適合した修正を行うことができる。なお、図１５に示したＧＭＴ１３ｇ１〜１３ｇ４はあくまでも一例であり、色相ｈが異なれば色域の大小関係も異なるものとなるため、色相ｈに応じてＧＭＴ１３ｇ１〜１３ｇ４による修正態様は異なるものとなる。むろん、第一表色系と第二表色系が異なるものとなれば、色域の大小関係も違ったものとなるため、第一表色系と第二表色系に応じてもＧＭＴ１３ｇ１〜１３ｇ４による修正態様は異なるものとなる。ステップＳ２７０にて変換されたＪＣｈ値が、ＧＭＴ１３ｇ１〜１３ｇ４に記述された代表的なＪＣｈ値であれば、対応するＪ’Ｃ’ｈ’値をＧＭＴ１３ｇ１〜１３ｇ４から読み取って、修正後の値として出力する。一方、ステップＳ２７０にて変換されたＪＣｈ値が、ＧＭＴ１３ｇ１〜１３ｇ４に直接記述されていないＪＣｈ値であれば、対応するＪ’Ｃ’ｈ’値をＧＭＴ１３ｇ１〜１３ｇ４に記述された代表的なＪＣｈ値とＪ’Ｃ’ｈ’値との対応関係に基づいて補間する。具体的には、ＪＣｈ表色空間において修正対象のＪＣｈ値の周囲の代表的なＪＣｈ値との位置関係に基づいて体積補間や面積補間等を行う。むろん、非線形補間も適用することができる。

ＧＭＴ１３ｇ１〜１３ｇ４を予め用意しておくことにより、画像データ１３ｂの色変換を行うたびに、同画像データ１３ｂの第一表色系の色域を計算し、ベクトルｖの長さを決定するような計算を行う必要がなく、高速に色変換を行うことができる。また、複数のＧＭＴ１３ｇ１〜１３ｇ４を用意しておくことにより、種々の第一表色系の画像データ１３ｂが入力された場合でも適合するガマット修正を行うことができる。さらに、複数のＧＭＴ１３ｇ１〜１３ｇ４を用意しておくことにより、画像データ１３ｂの画像タイプに適合したガマット修正を行うことができる。画像データ１３ｂの属性Ａ〜Ｄは予め取得されるため、各画素について色変換を行う際に、適切なＧＭＴ１３ｇ１〜１３ｇ４を参照することができる。

ガマット修正部Ｐ５がガマット修正後のＪ’Ｃ’ｈ’値を取得すると、ステップＳ２９０において、第二色変換部Ｐ６が環境パラメータ取得部Ｐ２を介してプリンタ２０（印刷物）に関する環境パラメータＸ_WＹ_WＺ_W，Ｌ_A，ｃ，Ｎ_c，Ｆを取得する。上述したとおりプリンタ２０に関する環境パラメータＸ_WＹ_WＺ_W，Ｌ_A，ｃ，Ｎ_c，Ｆは、予めＨＤＤ１３の環境パラメータ１３ｆに格納されているため、ＣＰＵ１１がこれを読み出して使用することができる。さらに、第二色変換部Ｐ６が、プリンタ２０（印刷物）に関する環境パラメータＸ_WＹ_WＺ_W，Ｌ_A，ｃ，Ｎ_c，ＦをＣＩＥＣＡＭ０２に代入しつつ、ＣＩＥＣＡＭ０２を用いてＪＣｈ値をＸ’Ｙ’Ｚ’値に変換する。ＣＩＥＣＡＭ０２は可逆変換プロファイルであるため、逆方向の変換においてＪＣｈ値をＸＹＺ値に変換することができる。

ステップＳ３００では、第二色変換部Ｐ６が予めＨＤＤ１３に記憶された色変換ＬＵＴ１３ｄを取得し、同取得した色変換ＬＵＴ１３ｄを参照してＸ’Ｙ’Ｚ’値に対応するＣＭＹＫ値を特定する。色変換ＬＵＴ１３ｄは、ＸＹＺ表色系とＣＭＹＫ表色系との等色関係を規定したテーブルであり、プリンタ２０にて印刷した際のＣＭＹＫ値が測色的にどのＸ’Ｙ’Ｚ’値に対応しているかを検証していくことにより予め作成されている。なお、ＣＭＹＫ表色系は本発明の第二表色系に相当し、ステップＳ２９０〜Ｓ３００が本発明の第二色変換工程に相当することとなる。

ステップＳ３１０においては画像データ１３ｂの全画素についてステップＳ２５０の選択が完了したかどうかが判断され、全て選択していない場合にはステップＳ２５０に戻り次の画素が選択される。すなわち、ステップＳ２５０〜Ｓ３１０を繰り返すことにより、順に画素をシフトさせていき、最終的には全画素について対応するＣＭＹＫ値を特定していくことができる。そして、ステップＳ３１０にて全画素についてＣＭＹＫ値の特定が完了したことが確認されると、各画素の色がＣＭＹＫ値で表現された色変換画像データを次の工程に出力する。なお、次の工程においては、色変換データがハーフトーン処理部Ｐ７と印刷データ生成部Ｐ８にて変換され、プリンタ２０に出力可能な印刷データが生成される。

このように、第一色変換工程においてはディスプレイ６０に表示された画像の色を、その観察環境下においてどのように知覚されるかを推定することができる。そして、第二色変換工程においては印刷結果の観察環境において、第一色変換工程にて推定された知覚色と同様に知覚される色が印刷できるＣＭＹＫ値を特定することができる。すなわち、本発明の色変換を行うことにより、ディスプレイ６０を観察しているときに知覚する色と、印刷結果を観察しているときに知覚する色とを同じにすることができ、ディスプレイ６０を見ながら作成した画像データをイメージどおりに出力することができる。

観察環境下においてどのように知覚されるかを推定するにあたってはＣＩＥＣＡＭ０２モデルが使用され、ＣＩＥＣＡＭ０２では環境パラメータＸ_WＹ_WＺ_W，Ｌ_A，ｃ，Ｎ_c，Ｆを変換式に代入することができる。このうち環境パラメータＸ_WＹ_WＺ_W，Ｌ_A，ｃ，Ｎ_c，Ｆは、照明光の分光分布Ｐ（λ）および照度Ｉから得ることができ、分光分布Ｐ（λ）および照度Ｉはテストパターンの視覚結果に基づいて特定することが可能となっている。従って、環境パラメータＸ_WＹ_WＺ_W，Ｌ_A，ｃ，Ｎ_c，Ｆを得るために分光光度計や照度計を用意する必要はなく、一般のユーザーにおいても容易に本発明を利用することができる。

なお、本実施形態の色変換処理においてはｓＲＧＢ値やａＲＧＢ値をＣＭＹＫ値に変換する第一色変換処理および第二色変換処理を画素ごとに一貫して行うようにしたが、ｓＲＧＢ表色系やａＲＧＢ表色系の元画像データの全画素について第一色変換部Ｐ４がｓＲＧＢ値をＪＣｈ値に変換することにより、一旦、ＪＣｈ表色系の画像データに変換し、さらにＪＣｈ表色系の画像データの全画素を第二色変換部Ｐ６がＣＭＹＫ値に変換することにより最終的にＣＭＹＫ表色系の色変換画像データを得るようにしてもよい。両者の差異は、第一色変換および第二色変換を画素単位で行うか画像データ単位で行うかという便宜的な事象に過ぎず、実質的には同じ処理を行っているということができる。さらに、一画素を画像データと捉えれば、本実施形態においても、ｓＲＧＢ表色系やａＲＧＢ表色系の元画像データを、一度、ＪＣｈ表色系の画像データに変換し、さらにＪＣｈ表色系の画像データをＣＭＹＫ表色系の色変換画像データに変換しているということができる。

（４）変換プロファイル（ＣＩＥＣＡＭ０２）について：
次に、ＣＩＥＣＡＭ０２について説明する。上述したとおりＣＩＥＣＡＭ０２は環境パラメータＬ_A，ｃ，Ｎ_c，Ｆを代入することにより、その観察環境にて知覚される色変換結果を得ることが可能なＸＹＺ−ＪＣｈ可逆変換プロファイルである。具体的には、下記に説明する手順によってＸＹＺ値をＪＣｈ値に変換する。

まず、下記式（６）の行列変換によって、ＸＹＺ表色系をＲＧＢ表色系に変換する。

次に、等エネルギー白色からデバイス白色への順応の度合いとなる順応ファクタＤを下記式（７）によって算出する。

ここで、環境パラメータ取得部Ｐ２が予め取得した環境パラメータＬ_A，Ｆが代入されることとなる。順応ファクタＤは、０〜１の値となり、１となる場合には完全順応となる。次に、上記式（６）によって算出したＲＧＢ値を下記式（８）の色順応式に代入する。
上記式（８）によって、順応ファクタＤにおける色順応を考慮したＲＧＢ値（Ｒ_cＧ_cＢ_c）を算出することができる。なお、上記式（８）におけるＹ_wは白色点のＹ値を意味し、Ｒ_wＧ_wＢ_wは白色点のＲＧＢ値をそれぞれ意味している。ここで、上述したように環境パラメータ取得部Ｐ２が取得した白色点の色Ｘ_wＹ_wＺ_wがＨＤＤ１３から読み出されて上記式（８）に代入される。なお、Ｒ_wＧ_wＢ_wは上記式（６）によってＸ_wＹ_wＺ_wから変換することができる。なお、ステップＳ２６０におけるＣＩＥＣＡＭ０２の順変換においてはディスプレイに関する環境パラメータＸ_WＹ_WＺ_Wが代入され、ステップＳ２９０におけるＣＩＥＣＡＭ０２の逆変換においてはプリンタ（印刷物）に関する環境パラメータＸ_WＹ_WＺ_Wが代入される。

次に、下記式（９）ＲＧＢ表色系で表されたＲ_cＧ_cＢ_c値をHunt-Pointer-Estevezの錐体刺激空間の座標値Ｒ’Ｇ’Ｂ’に変換する。

さらに、下記式（１０）〜（１４）に基づいて後の計算に使用するパラメータｋ，Ｆ_L，ｎ，Ｎ_bb，Ｎ_cb，ｚを算出する。
なお、本実施形態においては背景輝度Ｙ_b＝２０が代入される。
上記式（１１）においては、環境パラメータ取得部Ｐ２が予め取得した環境パラメータＬ_Aが代入されることとなる。

そして、上記式（９）によって算出したＲ’Ｇ’Ｂ’値を下記式（１５）に代入することにより、順応後の非線形圧縮をする変換を行う。これにより、錐体応答の入出力特性を考慮したＲ’Ｇ’Ｂ’値（Ｒ’_aＧ’_aＢ’_a）を算出することができる。
なお、上記式（１５）に代入するＲ’Ｇ’Ｂ’値が負の値である場合には、その絶対値を上記式（１５）に代入する。以上によりの環境パラメータＬ_A，ｃ，Ｎ_c，Ｆおよび背景輝度Ｙ_b1，Ｙ_b2による色順応を考慮したＲ’_aＧ’_aＢ’_a値を特定することができる。

次に、ここまでで得られたＲ’_aＧ’_aＢ’_a値からＪＣｈ表色系の座標値を特定する計算を行う。まず、Ｒ’_aＧ’_aＢ’_a値を下記式（１６），（１７）に代入することにより座標変換を行いａ，ｂを算出する。
さらに、ａ，ｂを下記式（１８）に代入することより色相角ｈを算出する。

色相角ｈが算出できると、次に下記式（１９）により離心率ｅ_tを算出する。
図１６は、代表的なユニーク色についての離心率ｅ_iと色相角ｈ_iと色相成分（hue quadrature）Ｈ_iとの対応関係を示している。この対応関係に基づいて、上記式（１５）によって得られたＲ’_aＧ’_aＢ’_a値に対応する色相成分Ｈを直線補間により算出する。具体的には、下記式（２０）によってＲ’_aＧ’_aＢ’_a値に対応する色相成分Ｈを算出する。
なお、上記式（２０）における色相角ｈ’は、原則的に上記式（１８）によって算出した色相角ｈと同じ値とし、赤の色相角ｈ_i（ｉ＝１）よりも色相角ｈが小さい場合に限りｈ’＝（ｈ＋３６０）とする。これにより、補間に使用する色相角ｈ’が必ずいずれかのユニーク色の色相角ｈ_iの間の値を取ることとなり、補間が可能となる。なお、補間に使用するユニーク色は、ｈ_i≦ｈ’＜ｈ_iとなるｉが選択される。

次に、上記式（１５）によって得られたＲ’_aＧ’_aＢ’_a値についての無彩色応答Ａを下記式（２１）によって算出する。
Ｒ’_aＧ’_aＢ’_a値の無彩色応答Ａが得られたら、次に明度Ｊを下記式（２２）によって算出する。
上記式（２２）においては、Ｒ’_aＧ’_aＢ’_a値の無彩色応答Ａと、刺激の無彩色応答Ａ_wとの比を算出し、この比を環境パラメータｃと上記式（１４）で得られたｚの積算値で累乗することにより明度Jを算出している。

次に、クロマＣを算出する。まず、クロマＣを算出するために必要な係数ｔを下記式（２３）によって算出する。
上記式（２３）では、環境パラメータ取得部Ｐ２が予め取得した環境パラメータＮ_cが代入されることとなる。係数ｔが得られると、下記式（２４）に係数ｔを代入して、クロマＣを算出する。
上記式（２４）によりクロマＣが算出できると、ＪＣｈ表色系の各座標値Ｊ，Ｃ，ｈがそれぞれ算出できたこととなる。これにより、元画像データ１３ｂにおいて各画素の色を表現するｓＲＧＢ値を、環境パラメータＸ_wＹ_wＺ_w，Ｌ_A，ｃ，Ｎ_c，Ｆおよび背景輝度Ｙ_bを代入可能な変換プロファイル（ＣＩＥＣＡＭ０２）を用いてＪＣｈ値に変換できたこととなる。

（５）まとめと変形例：
本発明では、第一色変換工程においてはディスプレイ６０に表示された画像の色を、その観察環境下においてどのように知覚されるかを推定することができる。そして、第二色変換工程においては印刷結果の観察環境において、第一色変換工程にて推定された知覚色と同様に知覚される色が印刷できるＣＭＹＫ値を特定することができる。また、ガマット修正をテーブルデータのＧＭＴ１３ｇを参照して行うことにより、処理効率を向上させることができる。さらに、画像データ１３ｂの属性Ａ〜Ｄに応じてＧＭＴ１３ｇ１〜ｇ４を選択して使用するため、適切なガマット修正を実現することができる。

上述した実施形態においては、属性取得部Ｐ３が画像データ１３ｂを空間周波数解析することにより画像データ１３ｂの属性を判定していたが、他の画像解析手法によって属性を判定してもよい。
例えば、デジタルカメラにて生成された写真画像においては、ＥＸＩＦ規格（ＪＥＩＴＡ ｖｅｒ．２．２）に準拠した画像データが生成される場合が多い。従って、画像データ１３ｂにＥＸＩＦ特有のヘッダ（付帯情報）が添付されている場合には、当該画像データ１３ｂの画像タイプが自然画であると判断することができる。一方、ＣＧを作成して保存する場合、ＪＥＰＧ形式やＢＭＰ形式等の汎用形式で保存することが多く、ＥＸＩＦ特有のヘッダが添付される可能性は低いということができる。従って、画像データ１３ｂのヘッダを解析することにより、画像タイプが自然画であるかＣＧ画であるかを容易かつ高速に特定することができる。

図１７は変形例にかかる色変換処理の流れを模式的に示し、図１８は本変形例にかかるコンピュータ１１０にて実行されるプログラムのソフトウェア構成を示している。両図において、本変形例においてはＧＭＴ１３ｇ１〜１３ｇ４から派生した派生ＧＭＴ１３ｉが色変換に使用される。前実施形態にて実行された各モジュールＰ１〜Ｐ８の他に、派生ＧＭＴ１３ｉを作成するための派生テーブル作成部Ｐ９と、派生ＧＭＴ１３ｉを使用して色変換を行うための第三色変換部Ｐ１０とがＣＰＵ上にて実行される。図１７に示したように派生ＧＭＴ１３ｉは、第一色変換工程におけるＸＹＺ値と第二色変換工程におけるＸ’Ｙ’Ｚ’値との対応関係を規定したテーブルであり、第三色変換部Ｐ１０は派生ＧＭＴ１３ｉを参照してＸＹＺ値をＸ’Ｙ’Ｚ’値に変換する。

図１９は、本変形例にかかる色変換処理の流れを示している。同図において、前実施形態と同様にステップＳ１２４０までの処理によって画像データ１３ｂの属性Ａ〜Ｄが取得される。次に、ディスプレイ６０およびプリンタ２０（印刷物）に関する環境パラメータＬ_A，ｃ，Ｎ_c，Ｆをそれぞれ取得する。本変形例においても予め環境パラメータ取得処理にてテストパターンを印刷することにより、環境パラメータＬ_A，ｃ，Ｎ_c，Ｆが取得される。ただし、本変形例においてテストパターンには照度エリアＳ２のみが形成されており、環境パラメータＬ_A，ｃ，Ｎ_c，Ｆのみが観察環境に応じた環境パラメータとして取得される。本変形例では白色点の色Ｘ_wＹ_wＺ_wが一定の値であると仮定され、ＣＩＥＣＡＭ０２の変換式には常に一定のＸ_wＹ_wＺ_wを代入するものとする。

ステップＳ１２５０においては、第三色変換部Ｐ１０が環境パラメータＬ_A，ｃ，Ｎ_c，Ｆと属性Ａ〜Ｄの組み合わせを取得する。次に、ステップＳ１２６０において、ＨＤＤ１３に記憶されている派生ＧＭＴ１３ｉを取得し、同派生ＧＭＴ１３ｉを作成したときの環境パラメータＬ_A，ｃ，Ｎ_c，Ｆと属性Ａ〜Ｄの組み合わせを取得する。そして、双方の組み合わせが一致する派生ＧＭＴ１３ｉを検索する。なお、派生ＧＭＴ１３ｉは基本的に色変換処理を行うごとに作成されるため、複数の派生ＧＭＴ１３ｉがＨＤＤ１３に蓄積されている。

図２０は、環境パラメータＬ_A，ｃ，Ｎ_c，Ｆと属性Ａ〜Ｄの組み合わせを表にして示している。同図上段において、順応輝度Ｌ_Aに関するテーブルが設けられており、順応輝度Ｌ_Aが５０［cd/m²］幅の区間に区切られている。そして、各区間に対して、中間の代表値２５，７５，１２５・・・４７５が示されている。まず、このテーブルを参照することにより、連続的な順応輝度Ｌ_Aは離散的な１０個の代表値２５，７５，１２５・・・４７５のいずれかに帰属させられる。このようにすることにより、環境パラメータＬ_A，ｃ，Ｎ_c，Ｆの組み合わせの個数を低減させることができる。なお、図７にて示すように他の環境パラメータｃ，Ｎ_c，Ｆの組み合わせは３通りしか存在しないため、ディスプレイ６０に関する環境パラメータＬ_A，ｃ，Ｎ_c，Ｆの組み合わせは１０×３＝３０通りになる。プリンタ２０（印刷物）に関する環境パラメータＬ_A，ｃ，Ｎ_c，Ｆの組み合わせも１０×１＝１０通りになる。環境パラメータ取得処理で述べたように、プリンタ２０（印刷物）に関しては環境パラメータＬ_A，ｃ，Ｎ_cは常に０．６９，１．０，１．０となるからである。属性Ａ〜Ｄは４通りであるため、ディスプレイ６０およびプリンタ２０の環境パラメータＬ_A，ｃ，Ｎ_c，Ｆと属性Ａ〜Ｄとの組み合わせは、３０×１０×４＝１２００通りとなる。

図２１は、派生ＧＭＴ１３ｉを示している。派生ＧＭＴ１３ｉは、色空間を全体的に網羅するグリッドに関してＸＹＺ値とＸ’Ｙ’Ｚ’値との対応関係を記述したボディと、派生ＧＭＴ１３ｉを作成したときのディスプレイ６０およびプリンタ２０の環境パラメータＬ_A，ｃ，Ｎ_c，Ｆと属性Ａ〜Ｄの組み合わせを記述したヘッダとから構成されている。このようにすることにより、各派生ＧＭＴ１３ｉがどのような条件で作成されたかを特定することができる。第三色変換部Ｐ１０は、ステップＳ１２６０においてＨＤＤ１３に記憶された複数の派生ＧＭＴ１３ｉのなかからステップＳ１２５０にて取得したディスプレイ６０およびプリンタ２０の環境パラメータＬ_A，ｃ，Ｎ_c，Ｆと属性Ａ〜Ｄの組み合わせに合致するものを検索する。

ステップＳ１２７０においては、ディスプレイ６０およびプリンタ２０の環境パラメータＬ_A，ｃ，Ｎ_c，Ｆと属性Ａ〜Ｄの組み合わせに合致する派生ＧＭＴ１３ｉが検索されたかどうかが判定され、合致するものが検索された場合には、ステップＳ１２９０にて検索された派生ＧＭＴ１３ｉについて使用可能である旨のフラグをオンにする。ステップＳ１２９０，Ｓ１３００においては、前実施形態と同様に画素を選択し、第一色変換部Ｐ３がその画素のＸＹＺ値を取得する。ステップＳ１３１０においては、第三色変換部Ｐ１０が、フラグがオンとなっている使用可能の派生ＧＭＴ１３ｉの有無を判断し、使用可能な派生ＧＭＴ１３ｉがある場合には、ステップＳ１３２０にて、その派生ＧＭＴ１３ｉを参照してＸＹＺ値をＸ’Ｙ’Ｚ’値に変換する。

一方、使用可能な派生ＧＭＴ１３ｉがない場合には、前実施形態と同様にステップＳ１３３０〜Ｓ１３５０にて第一色変換部Ｐ４によるＣＩＥＣＡＭ０２の順変換と、ガマット修正部Ｐ５によるＪ’Ｃ’ｈ’への修正と、第二色変換部Ｐ６によるＣＩＥＣＡＭ０２の逆変換がディスプレイ６０およびプリンタ２０の環境パラメータＬ_A，ｃ，Ｎ_c，Ｆと属性Ａ〜Ｄに応じて行われる。これにより、元のＸＹＺ値に基づいてＸ’Ｙ’Ｚ’値を得ることができる。ステップＳ１３６０においては、元のＸＹＺ値とステップＳ１３３０〜Ｓ１３５０による変換後のＸ’Ｙ’Ｚ’値との対応関係を記憶する。ステップＳ１３２０またはステップＳ１３３０〜Ｓ１３５０によって得られたＸ’Ｙ’Ｚ’値は、ステップＳ１３７０にて前実施形態と同様にＣＭＹＫ値に変換される。また、ステップＳ１３８０を実行することにより前実施形態と同様に全ての画素に関してＣＭＹＫ値が得られるまで各画素の変換が繰り返される。

以上の処理を行うことにより、ＣＭＹＫ表色系の画像データに色変換を行うことができる。また、すでに作成された派生ＧＭＴ１３ｉとディスプレイ６０およびプリンタ２０の環境パラメータＬ_A，ｃ，Ｎ_c，Ｆと属性Ａ〜Ｄの組み合わせが合致した場合には、ＣＩＥＣＡＭ０２の演算を行うことなく色変換を実現することができるため、処理の効率を向上させることができる。ステップＳ１３９０においては、オンになっていたフラグをオフにする。さらに、ステップＳ１４００においては、ステップＳ１３５０において記憶したＸＹＺ値とＸ’Ｙ’Ｚ’値との対応関係に基づいて派生テーブル作成部Ｐ９が派生ＧＭＴ１３ｉを作成する。多数の画素で構成される画像データ１３ｂの個々の画素についてＸＹＺ値とＸ’Ｙ’Ｚ’値との対応関係が得られているため、これらに基づいて派生ＧＭＴ１３ｉを生成することができる。具体的には、色空間の全体を網羅するグリッド上にてＸＹＺ値とＸ’Ｙ’Ｚ’値との対応関係が記述できるように派生ＧＭＴ１３ｉを作成する。その際、ＸＹＺ値とＸ’Ｙ’Ｚ’値との対応関係が得られていない領域については、近隣の対応関係に基づいて補間を行うことにより、グリッド上におけるＸＹＺ値とＸ’Ｙ’Ｚ’値との対応関係を推測する。

以上のようにして派生ＧＭＴ１３ｉのボディを作成することができる。派生ＧＭＴ１３ｉは、ステップＳ１３３０〜Ｓ１３５０を順次行うことによって得られたＸ’Ｙ’Ｚ’値とＸＹＺ値との対応関係に基づいて作成されるため、ディスプレイ６０およびプリンタ２０の環境パラメータＬ_A，ｃ，Ｎ_c，Ｆと属性Ａ〜Ｄを考慮したものとすることができる。さらに、このボディに対してステップＳ１３２０〜Ｓ１３４０を行った際のディスプレイ６０およびプリンタ２０の環境パラメータＬ_A，ｃ，Ｎ_c，Ｆと属性Ａ〜Ｄの組み合わせを特定するヘッダを添付する。これにより、図２１にて示した派生ＧＭＴ１３ｉを作成することができ、以降の色変換において第三色変換部Ｐ１０が使用することができる。また、新たな派生ＧＭＴ１３ｉについて、ディスプレイ６０およびプリンタ２０の環境パラメータＬ_A，ｃ，Ｎ_c，Ｆと属性Ａ〜Ｄの組み合わせの適合性を判断することができる。さらに、図２１に示した派生ＧＭＴ１３ｉに対して、さらにａＲＧＢ表色系のカラープロファイルや色変換ＬＵＴ１３ｄをマージすることも可能である。この場合、ａＲＧＢ表色系からＣＭＹＫ表色系へ直接色変換を行うことができ、第三色変換部Ｐ１０のみで色変換を実現することができる。

コンピュータのハードウェア構成図である。 プログラムのソフトウェア構成図である。 環境パラメータ取得処理のフローチャートである。 テストパターンの一例を示す図である。 色度エリアの一例を示す図である。 照明光の分光分布を示すグラフである。 グレーの分光反射率を示すグラフである。 照度エリアの一例を示す図である。 知覚レベルと空間周波数の関係を示すグラフである。 分解限界周波数の照度依存性を示すグラフである。 周囲比と環境パラメータとの対応関係を示す表である。 色変換処理の流れを示すフローチャートである。 色変換処理の流れを示す模式図である。 ＧＭＴを示す図である。 ＧＭＴを示すグラフである。 ユニーク色の色相角と色相成分と離心率を示す表である。 変形例にかかる色変換処理の流れを示す模式図である。 変形例にかかるプログラムのソフトウェア構成図である。 変形例にかかる色変換処理の流れを示すフローチャートである。 属性と環境パラメータの組み合わせを示す表である。 派生ＧＭＴを示す図である。

符号の説明

１０…コンピュータ，１０ａ…バス，１１…ＣＰＵ，１２…ＲＡＭ，１３…ＨＤＤ，１３ａ…プログラムデータ，１３ｂ…画像データ，１３ｃ…テストパターンデータ，１３ｄ…色変換ＬＵＴ，１３ｅ…ガマットデータ，１３ｆ…環境パラメータ，１３ｇ…ガマット修正テーブル，１３ｈ…分光分布データ，派生ＧＭＴ１３ｉ，１４…ＵＳＢＩ／Ｆ，１５…入力機器Ｉ／Ｆ，１６…ビデオＩ／Ｆ，２０…プリンタ，４０…マウス，５０…キーボード，６０…ディスプレイ，Ｐ…プリンタドライバ，Ｐ１…画像データ取得部，Ｐ２…環境パラメータ取得部，Ｐ３…属性取得部，Ｐ４…第一色変換部，Ｐ５…ガマット修正部，Ｐ６…第二色変換部，Ｐ７…ハーフトーン処理部，Ｐ８…印刷データ生成部，Ｐ９…派生テーブル作成部，Ｐ１０…第三色変換部，Ｓ１…色度エリア，Ｓ２…照度エリア

Claims (12)

1. 第一出力デバイスに入力可能な第一表色系で各画素の色が表現された元画像データを、第二出力デバイスに入力可能な第二表色系で各画素の色が表現された色変換画像データに変換する色変換装置において、
   上記第一出力デバイスおよび上記第二出力デバイスの出力結果を観察するときの環境パラメータを取得する環境パラメータ取得手段と、
   色変換の対象となる上記元画像データの属性を取得する属性取得手段と、
   上記第一出力デバイスの出力結果を観察するときの上記環境パラメータを代入した変換プロファイルを利用して、上記元画像データを各画素の色が第三表色系で表現された画像データに変換する第一色変換手段と、
   修正前後の対応関係を特定した複数の修正プロファイルから上記元画像データの属性に適合するものを選択し、同選択された修正プロファイルに基づいて上記第一色変換手段にて変換された画像データを上記第三表色系にて修正する修正手段と、
   上記第二出力デバイスの出力結果を観察するときの上記環境パラメータを代入した変換プロファイルを利用して、上記修正手段にて修正された画像データを上記色変換画像データに変換する第二色変換手段とを具備することを特徴とする色変換装置。
2. 上記修正プロファイルは、上記第一表色系による色空間の色域と上記第二表色系による色空間の色域との差に基づく修正を行うこと特徴とする請求項１に記載の色変換装置。
3. 上記修正プロファイルは、上記第三表色系における修正前後の色座標値が複数記述されたテーブルであることを特徴とする請求項１または請求項２のいずれかに記載の色変換装置。
4. 上記属性取得手段は、上記元画像データが採用する上記第一表色系の種類を上記属性として取得することを特徴とする請求項１から請求項３のいずれかに記載の色変換装置。
5. 上記属性取得手段は、上記元画像データを画像解析することにより上記元画像データの属性を取得することを特徴とする請求項１から請求項４のいずれかに記載の色変換装置。
6. 上記属性取得手段は、上記元画像データに添付された同画像データの付帯情報に基づいて上記元画像データの属性を取得することを特徴とする請求項１から請求項５のいずれかに記載の色変換装置。
7. 上記修正手段は、上記元画像データの属性が自然画であるとき修正前後で彩度の変動量が少ない上記修正プロファイルを選択することを特徴とする請求項１から請求項６のいずれかに記載の色変換装置。
8. 上記修正手段は、上記元画像データの属性がＣＧ画であるとき修正前後で明度の変動量が少ない上記修正プロファイルを選択することを特徴とする請求項１から請求項７のいずれかに記載の色変換装置。
9. 上記環境パラメータと上記修正プロファイルを使用して上記第一色変換手段と上記修正手段と上記第二色変換手段にて変換および修正を順次行うことにより特定される上記第一表色系と上記第二表色系との対応関係がテーブルとして記述された派生テーブルを作成する派生テーブル作成手段と、
   上記環境パラメータ取得手段と上記属性取得手段が取得した上記環境パラメータと上記属性との組み合わせが、上記派生テーブルを作成したときの上記環境パラメータと上記修正プロファイルとの組み合わせと合致する場合に、当該派生テーブルを参照して上記元画像データを上記色変換画像データに変換する第三色変換手段とを具備することを特徴とする請求項１から請求項８のいずれかに記載の色変換装置。
10. 上記環境パラメータ取得手段は、テストパターンを印刷媒体に出力し、その視覚結果に基づいて、上記環境パラメータを取得することを特徴とする請求項１から請求項９のいずれかに記載の色変換装置。
11. 第一出力デバイスに入力可能な第一表色系で各画素の色が表現された元画像データを、第二出力デバイスに入力可能な第二表色系で各画素の色が表現された色変換画像データに変換する色変換方法において、
    上記第一出力デバイスおよび上記第二出力デバイスの出力結果を観察するときの環境パラメータを取得し、
    色変換の対象となる上記元画像データの属性を取得し、
    上記第一出力デバイスの出力結果を観察するときの上記環境パラメータを代入した変換プロファイルを利用して、上記元画像データを各画素の色が第三表色系で表現された画像データに変換し、
    修正前後の対応関係を特定した複数の修正プロファイルから上記元画像データの属性に適合するものを選択し、同選択された修正プロファイルに基づいて、第三表色系に変換された画像データを上記第三表色系にて修正し、
    上記第二出力デバイスの出力結果を観察するときの上記環境パラメータを代入した変換プロファイルを利用して、修正された画像データを上記色変換画像データに変換することを特徴とする色変換方法。
12. 第一出力デバイスに入力可能な第一表色系で各画素の色が表現された元画像データを、第二出力デバイスに入力可能な第二表色系で各画素の色が表現された色変換画像データに変換する機能をコンピュータ上にて実現させる色変換プログラムにおいて、
    上記第一出力デバイスおよび上記第二出力デバイスの出力結果を観察するときの環境パラメータを取得する環境パラメータ取得機能と、
    色変換の対象となる上記元画像データの属性を取得する属性取得機能と、
    上記第一出力デバイスの出力結果を観察するときの上記環境パラメータを代入した変換プロファイルを利用して、上記元画像データを各画素の色が第三表色系で表現された画像データに変換する第一色変換機能と、
    修正前後の対応関係を特定した複数の修正プロファイルから上記元画像データの属性に適合するものを選択し、同選択された修正プロファイルに基づいて上記第一色変換機能にて変換された画像データを上記第三表色系にて修正する修正機能と、
    上記第二出力デバイスの出力結果を観察するときの上記環境パラメータを代入した変換プロファイルを利用して、上記修正機能にて修正された画像データを上記色変換画像データに変換する第二色変換機能とをコンピュータ上にて実現させることを特徴とする色変換プログラム。