JP3803441B2

JP3803441B2 - 色変換方法

Info

Publication number: JP3803441B2
Application number: JP31973696A
Authority: JP
Inventors: 文人竹本
Original assignee: Fuji Photo Film Co Ltd
Current assignee: Fujifilm Holdings Corp
Priority date: 1996-11-29
Filing date: 1996-11-29
Publication date: 2006-08-02
Anticipated expiration: 2016-11-29
Also published as: JPH10164381A

Description

【０００１】
【発明の属する技術分野】
この発明は、例えば、カラーデジタルカメラ（以下、単にデジタルカメラともいう。）から出力されるデジタル信号であるＲＧＢ（ＲＧＢは、それぞれ、赤、緑、青の意）信号を測色値信号に変換する色変換方法に関する。
【０００２】
【従来の技術】
カラーデジタルカメラの出力信号であるＲＧＢ信号のようなデバイス依存（デバイスディペンデント）の画像信号を、３刺激値信号ＸＹＺのようなデバイス非依存（デバイスインディペンデント）の画像信号に色変換する従来技術として、例えば、特開平２−２９１７７７号公報または特開平２−２９１７７８号公報に開示された色変換に係る技術（第１の技術という。）がある。
【０００３】
この第１の技術は、ＲＧＢ信号を説明変数とし、このＲＧＢ信号のそれぞれを１／３乗した後、１／３乗後のＲＧＢ信号に対して３×ｎ（ｎは３の倍数）のマトリクスを作用させて目的変数としてのＣＩＥ−Ｌ^*ａ^*ｂ^*信号を得る技術である。
【０００４】
また、デバイス依存の印刷のｃｍｙ網％信号からデバイス非依存の３刺激値信号ＸＹＺを予測したり、あるいは、ｃｍｙの３原色信号に基づいて再現色の濃度を予測する従来の技術として、例えば、特開平４−３３７９６５号公報または特開平４−３３７９６６号公報に開示された印刷色再現色予測に係る技術（第２の技術という。）がある。
【０００５】
この第２の技術は、基準の重回帰モデルに対して、予め設定した説明変数群の中から変数を一つずつ順に選択して重回帰モデルを再構築し、ある許容範囲内に入ったか否かを判定する手順を繰り返しながら、重回帰モデルの精度を上げていく技術である。
【０００６】
【発明が解決しようとする課題】
しかしながら、これらの第１および第２の技術で対象としている説明変数は、ｃｍｙ信号、あるいはＲＧＢ信号とその２乗項、クロス項等である。このような場合、すなわち、例えば、Ｒ信号と、このＲ信号の２乗信号等の関係のような場合には、説明変数相互の相関が高く、統計学上、不安定なモデルとなることが知られている。例えば、目的変数ａ ^*と、Ｒ信号項と、Ｒ信号の２乗項の説明変数とは、正の相関関係を有するが、Ｒ信号項の偏回帰係数は正の値となり、Ｒ信号の２乗項の偏回帰係数は負の値となる場合がある（これを統計学用語で多重共線性という。）。
【０００７】
したがって、このような要素を有するマトリクスを用いて色信号を変換する場合、マトリクス作成に使用した色以外については色変換の精度を保障することができないという問題がある。ただし、第２の技術では、説明変数を選択することで、この問題をある程度解決している。
【０００８】
なお、色を正確に変換する技術として、例えば、画像電子学会誌第１８巻第５号（１９８９年）に発表された、予め、３次元（ＸＹＺ）の色空間を例えば５１２個の立方体に領域分割し、合計７２９点の格子点での色修正値を出力デバイスの特性に最適化した方法で計算して、いわゆるルックアップテーブル（ＬＵＴ）として持ち、各格子点間の入力値についてはこのＬＵＴを３次元的に補間して求める補間方法を使用する技術（第３の技術またはＬＵＴ＋補間方式による技術という。）が知られているが、この第３の技術では、メモリを多量に使用するという欠点があり、また、ＬＵＴを作成するために、規則性のある多数のカラーパッチを作成し、かつ測定する必要があることから時間がかかるという欠点もある。
【０００９】
この発明はこのような課題を考慮してなされたものであり、マトリクス方式の精度不足、並びにＬＵＴ＋補間方式におけるパッチの作成、その測定およびメモリ上の問題を一掃することを可能とする色変換方法を提供することを目的とする。
【００１０】
【課題を解決するための手段】
この発明は、第１の色空間の任意の色信号を第２の色空間の色信号に変換する色変換方法において、所望の代表的な色を彩度、明度、色相に応じて変化させた複数のカラーパッチを有するカラーチャートを前記第１の色空間の色信号を出力するデバイスで読み取る過程と、読み取った第１の色空間の色信号を輝度値又は濃度値に変換し、変換した前記輝度値又は前記濃度値の少なくとも一次項の輝度値又は濃度値を第１の説明変数群として求める過程と、前記第１の説明変数群を主成分分析して求めた前記輝度値又は前記濃度値の主成分スコアを相互の説明変数が直交する第２の説明変数群とする過程と、測色計により求めた前記カラーチャートの各カラーパッチの色についてのＸＹＺまたはＬ ^* ａ ^* ｂ ^* の各測色値を目的変数群とし、この目的変数群と前記第２の説明変数群を重回帰分析して、偏回帰係数のマトリクスを求める過程と、前記偏回帰係数の検定を行い、前記第２の説明変数群の中、所定程度有意な説明変数を選択し、この選択した説明変数を使用可能な説明変数とし、この使用可能な説明変数に対応する前記偏回帰係数のマトリクスを所定程度有意な偏回帰係数のマトリクスとする過程とを有し、前記第１の色空間の任意の色信号に対して前記使用可能な説明変数に対応する前記所定程度有意な偏回帰係数のマトリクスを作用させて、前記第２の色空間の色信号に変換することを特徴とする。
【００１１】
この発明によれば、第１の説明変数群を主成分分析し、相互の説明変数が直交する第２の説明変数群に変換し、この第２の説明変数群と、カラーチャートの測色値である目的変数群とから偏回帰係数をデータ量の少ないマトリクスとして求めるようにしている。このため、色変換の精度が保障され、かつマトリクス方式であるので、メモリ容量を少なくすることができる。
【００１２】
また、この発明によれば、求めた偏回帰係数の検定を行って、所定程度（例えば、１％、５％等）有意な第２の説明変数群を選択し、この選択した第２の説明変数群と前記偏回帰係数とでマトリクスを作成するようにしている。このため、一層、精度よく色変換でき、かつマトリクス方式であるので、メモリ容量を少なくすることができる。
【００１３】
この場合、前記偏回帰係数のマトリクスを求める過程では、前記カラーチャート中、所望の特定色に対して重み付けを行うことで、その所望の特定色の色変換精度を向上させることができる。
【００１４】
また、重み付け用の所望の特定色を灰色としておくことで、一般的に重要な色である灰色を精度よく色変換することができる。
【００１５】
さらに所望の特定色を、第１の色空間の色信号を統計処理して得られる、例えば、頻度の高い色とすることで、原シーンの色を精度よく色変換することができる。また、統計処理用の色として、肌色、青色、緑色を抽出することで、顔の肌色、空の青色、葉の緑色を特に正確に色変換することができる。
【００１６】
さらに、統計処理を、前記カラーチャートの各カラーパッチの色を中心とした領域に分割して行うことにより、実際の画像に存在している色を中心として色変換することができるので、色変換精度を高くできる。
【００１７】
なお、第１の色空間をＲＧＢ空間とし、第２の色空間をＸＹＺまたはＬ ^* ａ ^* ｂ ^* の測色値空間とすることで、種々のデバイスのＲＧＢ信号について、共通の色空間で画像処理を行うことができる。
【００１８】
さらにまた、重み付けするための所望の特定色は、ＲＧＢ信号を、一旦ＸＹＺまたはＬ ^* ａ ^* ｂ ^* の測色値信号に変換した後、その測色値空間の中で統計処理を行なって決定することもできる。
【００１９】
なお、前記カラーチャートとして、市販の既製品であるマクベスカラーチェッカー（登録商標）を用いた場合には、カラーパッチの数は２４色分と少なくて済み、かつ市販であるので、容易にカラーチャートを入手することができる。
【００２０】
【発明の実施の形態】
以下、この発明の一実施の形態について図面を参照して説明する。
【００２１】
図１は、この発明の一実施の形態が適用された画像信号処理装置１０の全体的な構成を示している。
【００２２】
この画像信号処理装置１０は、測色変換部１６を有し、この測色変換部１６には、撮像装置であるデジタルカメラ（不図示）で撮影したシーン（場面）の画像情報を有するＲＧＢ信号（デジタル画像信号）１２が供給されるとともに、必要に応じて、画像読取装置であるスキャナ（不図示）を構成する画像入力部（線走査読取装置またはスキャナ入力部ともいう。）で読み取られたリバーサル原稿の画像情報を有するＲＧＢ信号（デジタル画像信号）１４が供給される。なお、デジタルカメラはデジタルビデオカメラでもよく、その他、撮像機能を有したデジタル画像信号（デジタル画像データ）を出力する装置であれば、この発明を適用することができる。
【００２３】
また、スキャナとしては、リニアイメージセンサを搭載するカラースキャナである線走査読取装置でもよく、また、エリアイメージセンサを搭載するカラースキャナである面走査読取装置でも、この発明を適用することができる。
【００２４】
測色変換部１６は、デジタルカメラのＲＧＢ信号１２をＸＹＺまたはＬ^*ａ^*ｂ^*等の測色値信号（ＸＹＺ信号ともいう。）１８に変換する測色変換マトリクス（以下、単にマトリクスともいう。）２０と、スキャナ入力部のＲＧＢ信号１４をＸＹＺ信号またはＬ^*ａ^*ｂ^*等の測色値信号２２に変換する測色変換テーブル（測色変換ルックアップテーブル、ルックアップテーブルまたは単にテーブルともいう。）２４とを有している。
【００２５】
なお、通常、ＲＧＢ信号１２、１４は、デバイスに依存する（デバイス依存またはデバイスディペンデントともいう。）信号（データ）といい、これに対して測色値信号１８、２２は、デバイスに依存しない（デバイス非依存またはデバイスインディペンデントともいう。）信号（データ）という。
【００２６】
デジタルカメラの測色値信号１８は、撮影光源と観察光源の違いを吸収する光源変更部２５と、ハイライト部濃度とシャドー部濃度とをセットアップする測色的セットアップ部（単に、セットアップ部ともいう。）２６と、原シーン（撮影されたもとの場面）の前記測色値信号１８をリバーサル原稿上の色素濃度信号２８に変換する色素濃度変換部３０と、ｃｍｙｋ変換部３２を構成する原シーン忠実再現テーブル４２に必要に応じて供給される。
【００２７】
一方、スキャナ入力部の測色値信号２２は、前記測色的セットアップ部２６と、標準条件再現テーブル部４３を構成し前記測色値信号２２を色素濃度信号３４に変換する色素濃度変換部３６に供給される。
【００２８】
色素濃度変換部３０、３６から出力される色素濃度信号（ｃｍｙ信号）２８、３４は、スイッチ４５を通じて選択され、標準条件再現テーブル４４に供給される。
【００２９】
ｃｍｙｋ変換部３２は、基本的には、原シーン忠実再現テーブル４２と標準条件再現テーブル部４３とを有している。原シーン忠実再現テーブル４２は、供給されるデジタルカメラの測色値信号１８を測色的に保存された網％信号であるｃｍｙｋ信号４６に変換するルックアップテーブルである。
【００３０】
標準条件再現テーブル部４３は、測色値信号２２を色素濃度信号３４に変換する色素濃度変換部３６と、色素濃度変換部３０または色素濃度変換部３６のいずれかの出力信号を選択するスイッチ（マルチプレクサ、選択手段）４５と、標準条件再現テーブル４４とを有している。標準条件再現テーブル４４は、スイッチ４５により選択された色素濃度信号２８、３４のいずれかを網％信号であるｃｍｙｋ信号４８に変換する処理を行う。ｃｍｙｋ変換部３２は、３色表色系の信号を４色表色系の信号に変換する、いわゆる３色４色変換機能を有する。
【００３１】
ｃｍｙｋ変換部３２から出力される網％信号であるｃｍｙｋ信号４６、４８は、このｃｍｙｋ信号４６、４８に基づく画像を出力する画像出力部３５に供給される。画像出力部３５は、例えば、図示しない２値化変換部と、レーザ露光走査部（イメージセッタ等）と、現像部と、刷版作成部と、印刷部とからなる公知の構成を採用することができる。２値化変換部では、ｃｍｙｋ信号４６またはｃｍｙｋ信号４８をスクリーン線数、網角度等の出力条件に応じて選択されるＣＭＹＫの各閾値マトリクスと比較して、２値化処理する。レーザ露光走査部では、この２値信号（２値画像信号ともいう。）に基づいてオンオフするレーザビームによりフイルムを露光走査して潜像を形成する。現像部では、この潜像が形成されたフイルムを現像して画像を顕像化し、製版用フイルムを作成する。刷版作成部では、この製版用フイルムから刷版を作成する。印刷部では、前記刷版、この場合、４版分のＣＭＹＫ用刷版を印刷機に装着し、刷版に付けられた４色のインキが本紙（印刷用紙）に転写されることで、画像が形成されたハードコピーとしての印刷物が作成される。
【００３２】
なお、画像出力部３５として、フイルムの現像処理が不要であり、網点やスクリーン線数やスクリーン角度を直接本紙に網点画像として印刷してシミュレーションすることのできる、ダイレクトデジタルカラー校正（ＤＤＣＰ）システムを用いることもできる。
【００３３】
前記測色的セットアップ部２６は、測色値信号１８、２２を測色的ＥＮＤ（等価中性濃度）信号５２に変換する測色的ＥＮＤ順変換・逆変換マトリクス（単に、マトリクスともいう。）５０と、測色的ＥＮＤ信号５２から間引き処理により作成した測色的ＥＮＤ信号５６を出力するラフデータ作成部５４と、作成された測色的ＥＮＤ信号５６に基づいて、ハイライト濃度信号６０とシャドー濃度信号６２とを自動的に決定するオートセットアップ部５８と、測色的ＥＮＤ信号５２に対して階調変換処理を施して測色的ＥＮＤ信号６６に変換するＥＮＤ・ＥＮＤ変換部６４と、変換された測色的ＥＮＤ信号６６を測色値信号２２に逆変換する測色的ＥＮＤ順変換・逆変換マトリクス５０とを有する。
【００３４】
なお、画像信号処理装置１０は、図示していないコンピュータ（ＣＰＵ、ＲＯＭ、ＲＡＭ、外部記憶装置、モニタ、その他、入出力機器等を含む。）により制御される構成となっており、画像信号処理装置１０を構成する各ブロックは、ハードウエアばかりでなく、ソフトウエアで構成される部分をも有する。コンピュータは、制御、判断、演算、比較手段等として機能する。
【００３５】
次に、画像信号処理装置１０を構成する各ブロックの詳細な構成と動作について説明する。
【００３６】
測色変換部１６を構成し、ＲＧＢ信号１２をＸＹＺ信号１８に変換するマトリクス２０は、図２に示すフローチャートに基づいて作成される。なお、以下の説明において、マトリクス２０は、ＲＧＢ色空間からＣＩＥ−Ｌ^*ａ^*ｂ^*（光源：補助標準の光ＣＩＥ−Ｄ５０）色空間への変換を例として説明する。この場合、ＣＩＥ−Ｌ^*ａ^*ｂ^*色空間とＸＹＺ色空間との相互の変換は、次に示す公知の（１）式により一意に行うことが可能である。したがって、以下の全ての説明において、ＸＹＺ色空間（またはＣＩＥ−Ｌ^*ａ^*ｂ^*色空間）での処理はＣＩＥ−Ｌ^*ａ^*ｂ^*色空間（またはＸＹＺ色空間）での処理に置き換えることができる。また、その他、これらと等価な測色的色空間での処理に置き換えることができる。
【００３７】
Ｌ^*＝１１６（Ｙ／Ｙｎ）^1/3−１６
ａ^*＝５００｛（Ｘ／Ｘｎ）^1/3−（Ｙ／Ｙｎ）^1/3｝
ｂ^*＝２００｛（Ｙ／Ｙｎ）^1/3−（Ｚ／Ｚｎ）^1/3｝ …（１）
まず、代表的な色を彩度、明度、色相に応じて変化させた複数のカラーパッチ７２（図２参照）を有する一種の色票であるカラーチャート７０を準備する（ステップＳ１）。この実施の形態では、カラーチャートとして、マクベスカラーチェッカー｛登録商標：米国のコールモージェン社のマクベス部門（Macbeth A division Kollmorgen ）製｝を用いる。マクベスカラーチェッカーは、公知のように、ＣＩＥ（１９３１）ｘｙＹ値、色相、マンセル表記値、彩度が規定されたカラーチャートである。
【００３８】
２４色は、具体的には、
１．暗い肌（dark skin ）
２．明るい肌（light skin）
３．青い空（blue sky）
４．葉（foliage ）
５．青い花（blue flower ）
６．青みの緑（bluish green）
７．オレンジ（orange）
８．紫みの青（purplish green）
９．明度彩度のほどよい赤（moderate red）
１０．紫（purple）
１１．黄緑（yellow green）
１２．黄だいだい（orange yellow ）
１３．青（Blue）
１４．緑（Green ）
１５．赤（Red ）
１６．黄（Yellow）
１７．マゼンタ（Magenta ）
１８．シアン（Cyan）
１９．白（White ）
２０．中性８（neutral ８：明るい灰色で、８はマンセル表記値の８）
２１．中性６．５（neutral ６．５：ライトミディアム灰色）
２２．中性５（neutral ５：中間の灰色）
２３．中性３．５（neutral ３．５：暗い灰色）
２４．黒（black ）
である。
【００３９】
なお、カラーチャートとしては、マクベスカラーチェッカーに限らず、例えば、ＪＩＳ標準色標等、色空間を概ね均等に網羅したカラーチャートを用いることができる。
【００４０】
次に、ＣＩＥ−Ｄ５０の撮影光源下でデジタルカメラを用いてカラーチャート７０の２４色の各色、すなわち２４個の各パッチ７２を撮影して、パッチ７２のそれぞれのＲＧＢ信号１２を得、得られた各ＲＧＢ信号１２を輝度値に変換した後、１／３乗する（ステップＳ２）。なお、輝度値への変換は、例えば、デバイス内部で施されているγ補正を解除することで求めることができる。また、１／３乗するのは、上記（１）式からも理解されるように、得られたＲＧＢ信号１２をＣＩＥ−Ｌ^*ａ^*ｂ^*表色系で処理するためである。
【００４１】
次いで、ステップＳ２で得られた各パッチ７２についてのＲＧＢ信号１２、具体的には、各パッチ７２について、輝度変換され１／３乗されたＲ値、Ｇ値、Ｂ値から２次項までの各値Ｒ、Ｇ、Ｂ、Ｒ²、Ｇ²、Ｂ²、ＲＧ、ＧＢ、ＢＲ（９変数）を計算する（ステップＳ３）。
【００４２】
次に、後に説明する重回帰分析を行った場合に、マルチコ（多重共線性）現象を起こさないように、ステップＳ３で得られた９変数の２４色のデータに対して主成分分析を行い、９変数の主成分スコア（主成分得点）Ｖを求める。各色毎に次の（２）式で示す主成分スコアＶ（Ｖは、ベクトルと考える。）が求められる。
【００４３】
Ｖ＝（ｖ₁，ｖ₂，ｖ₃，ｖ₄，ｖ₅，ｖ₆，ｖ₇，ｖ₈，ｖ₉）…（２）
なお、この（２）式において、主成分スコアＶの成分ｖ₁，ｖ₂，ｖ₃，ｖ₄，ｖ₅，ｖ₆，ｖ₇，ｖ₈，ｖ₉は、相互に相関が全くなく前記マルチコ現象を起こさない。
【００４４】
次に、図示していない測色計により、カラーチャート７０の各パッチ７２、すなわち２４色の各色についての測色値Ｃ（Ｌ^*ａ^*ｂ^*）（Ｃも、２４色分得られるので、ベクトルと考える。）を求める（ステップＳ５）。この測色値Ｃを求める過程は、ステップＳ２〜Ｓ４の過程に対してどの時点で求めてもよい。
【００４５】
次に、測色値Ｃを目的変数（従属変数）とし、主成分スコアＶを説明変数（独立変数）として、重回帰分析により偏回帰係数Ａ（Ａもベクトルと考える。）を求める（ステップＳ６）。
【００４６】
この重回帰分析を行う際に、目的変数群となっている測色値Ｃを構成する２４色に１：１に対応した重み付けマトリクス（単に、マトリクスともいう。）Ｐ＝［Ｐｉ］（ｉ＝１，２，…２４）を作用させる（ステップＳ７）。Ｐｉは、各色の重みであり、例えば、上述の２４色中、２０．中性８（neutral ８：明るい灰色で、８はマンセル表記値の８）、２１．中性６．５（neutral ６．５：ライトミディアム灰色）、２２．中性５（neutral ５：中間の灰色）、２３．中性３．５（neutral ３．５：暗い灰色）の灰色についての重みを大きくすることで、灰色についての色の再現性を向上させることができる。灰色に限らず、目的に応じて所望の色、例えば、人であれば肌色、外の景色であれば空色について限定して重みを大きくすれば、肌色または空色をより忠実に再現することが可能となり、好適である。また、重みは、入力ＲＧＢ信号１２を統計的に処理し、一画面分の画像信号、言い換えれば、画素分割された画像データについて、例えば、カラーチャートの各パッチの色を中心とした領域に色空間を分割し、各領域に存在する画素の頻度値に応じた重み付け、例えば、最頻値の重みを最も大きくし、以下、順に重みを小さくするようにすれば自動的な重み付け処理を行うことができる。なお、色空間の分割は、ＲＧＢ色空間であってもよく、重みなしのマトリクスで変換されたＸＹＺ（Ｌ^*ａ^*ｂ^*）色空間であってもよい。
【００４７】
この場合、重みを大きくした色を画像の主要色と呼ぶ。なお、マトリクスＰの各要素の合計値が１になるように、各要素の値を各要素の合計値で割った値に規格化しておく。
【００４８】
偏回帰係数Ａを求めるステップＳ６の重回帰分析について詳しく説明する。
【００４９】
目的変数である測色値（ベクトル）Ｃと、求めようとする偏回帰係数（ベクトル）Ａと、主成分スコア（ベクトル）Ｖとの間で、２４色の各測色値Ｃに対して次の（３）式に示す線形１次式が成立するものとする。
【００５０】
Ｃ＝ＡＶ …（３）
この（３）式と等価なマトリクスによる表現式を（４）式に、和の記号Σによる表現式を（５）式に各々示す。
【００５１】
【数１】

【００５２】
【数２】

【００５３】
なお、（５）式のように和の記号Σを使用する式を、必要に応じて、Ｌ^*＝Σ（ｊ＝１→１０）ａ_1jｖ_jと表すこととする。
【００５４】
偏回帰係数Ａは、Ｌ^*ａ^*ｂ^*のそれぞれ別個に最小二乗法を用いて求める。例えば、Ｌ^*については、次の（６）式によるｅ_Lを最小とする偏回帰係数ａ_1jを求めればよい。
【００５５】
ｅ_L＝Σ（ｉ＝１→２４）Ｐｉ｛Ｃｉ−Σ（ｊ＝１→１０）ａ_1jｖ_j｝²…（６）
ここで、ｉは、カラーチャート７０のパッチ番号、Ｐｉは、各色の重み、ｊは、変数の番号（１、２、…、１０）である。
【００５６】
（６）式をベクトルとマトリクスで表すと（７）式となる。ただし、（７）式において、測色値Ｃと偏回帰係数ａとはベクトル、主成分スコア［Ｖ］と重み付けマトリクス［Ｐ］とはマトリクスである。ｔはマトリクスの転置を表す。
【００５７】

以下、ベクトルＣは単にＣと記載し、ベクトルａは単にａと記載する。先ず、（７）式は、以下のように変形することができる。
【００５８】
ｅ_L＝（Ｃ^t−［Ｖ］^tａ^t）^t［Ｐ］（Ｃ^t−［Ｖ］^tａ）
ここで、ａと［Ｖ］とを入れ換えれば、
ｅ_L＝（Ｃ−ａ［Ｖ］）［Ｐ］（Ｃ^t−［Ｖ］^tａ）となる。
【００５９】
一般的に、（ＡＢＣ）^t＝Ｃ^tＢ^tＡ^tであるので、
ｅ_L＝Ｃ［Ｐ］Ｃ^t＋ａ［Ｖ］［Ｐ］［Ｖ］^tａ^t−ａ［Ｖ］［Ｐ］Ｃ^t−Ｃ［Ｐ］［Ｖ］^tａ^tとなる。
【００６０】
ここで、［Ｖ］［Ｐ］［Ｖ］^t＝［Ｎ］、［Ｖ］［Ｐ］Ｃ^t＝［Ｕ］と置けば、
ｅ_L＝Ｃ［Ｐ］Ｃ^t＋ａ［Ｎ］ａ^t−２ａ［Ｕ］ …（８）
となる。（８）式において、ｅ_Lを最小にするためには、偏回帰係数ａの各要素についての微分が０に等しくならなければならない。従って、次の（９）式が成立する。
【００６１】
【数３】

【００６２】
この（９）式から、以下に示すように偏回帰係数ａを（１０）式により求めることができる。
【００６３】

（［］^-1は逆行列を示す。）
このようにしてＬ^*についての偏回帰係数ａ（ａ_1j）が１０個求められ、以下、同様の方法で、残りのａ^*、ｂ^*についても偏回帰係数ａ_2j、ａ_3jを各１０個求めることができる。求められた偏回帰係数ａ_1j、ａ_2j、ａ_3jの合計は、３×１０個である｛（４）式参照｝。
【００６４】
結局、Ｌ^*ａ^*ｂ^*についての偏回帰係数Ａをまとめて表現すると、（１１）式で表すことができる。
【００６５】
Ａ＝ＣＰＶ^t（ＶＰＶ^t）^-1 …（１１）
次に、各目的変数Ｌ^*ａ^*ｂ^*について求めた偏回帰係数Ａの例えば５％有意で検定を行い（ステップＳ８）、５％有意（信頼度９５％）な説明変数Ｖを記憶手段に格納しておくとともに、偏回帰係数Ａを図１に示したマトリクス２０として格納する（ステップＳ９）。なお、１％有意（信頼度９９％）で検定してもよい。
【００６６】
ステップＳ８の検定では、（１２）式〜（１６）式に示すように、回帰平方和Ｓ_R｛求めた推定値Ｌ^*ｉ（Ｌの上に記号「＾」を付けている。）と測色値の平均Ｌ^*ｉ（Ｌの上に記号「^-」を付けている。）の差の２乗和：（１２）式参照｝と残差平方和Ｓ_E｛測色値Ｌ^*ｉと推定値Ｌ^*ｉの差の２乗和：（１３）式参照｝とを求め、さらに回帰平方和Ｓ_Rと残差平方和Ｓ_Eの不偏分散Ｖ_R、Ｖ_E｛（１４）式、（１５）式参照｝を各々求める。さらに、偏Ｆ値Ｆｊ｛（１６）式参照｝を求める。なお、（１６）式において、ａの頭に「＾」を付けたものは、（４）式の右辺のマトリクスの係数の推定値、Ｓ^jjは、説明変数ｖ_jの分散共分散行列の逆行列の対角項を意味している。
【００６７】
【数４】

【００６８】
この偏Ｆ値Ｆｊが５％有意のＦ分布を参照して求めた値＝Ｆ′_n-p-1（０．０５）＝Ｆ′_24-9-1（０．０５）＝Ｆ′₁₄（０．０５）＝４．６００１１より大きい場合には、５％有意水準で回帰が有意であり、予測に役立つ偏回帰係数と判断してマトリクス２０（図１参照）に格納しておく。
【００６９】
このようにして、デジタルカメラにより得られたＲＧＢ信号１２（図１参照）をＸＹＺ測色値信号１８に変換するＬ^*ａ^*ｂ^*を求めるためのマトリクス２０を作成することができる。このマトリクス２０に係る式を（１７）式に示す。
【００７０】
Ｌ^*＝Σａ_1jｖ_j （ｊ：５％有意な変数）
ａ^*＝Σａ_2jｖ_j （ｊ：５％有意な変数）
ｂ^*＝Σａ_3jｖ_j （ｊ：５％有意な変数） …（１７）
この（１７）式において、（ｊ：５％有意な変数）の意味は、変数として前記検定が５％有意であったもののみを使用するという意味である。その点で上述の（５）式と異なる。なお、ａ_1j、ａ_2j、ａ_3jは、偏回帰係数マトリクスＡの要素であり、ｖ_jは、ステップＳ４で求めた説明変数の主成分スコアである。
【００７１】
このマトリクス２０を使用することにより、既成のカラーチャート７０を使用して測色変換を行うことができる。この場合、ルックアップテーブルを使用しないのでメモリ容量が少なくて済み、換言すれば、メモリ容量が少なくても、精度よく測色変換を行うことができる。さらに、重み付けマトリクスＰの使用により画像の所望色（主要色）に限定して精度よく色変換を行うこともできる。なお、説明変数の数が９箇の場合には、ステップＳ８の検定処理を行わず、全ての偏回帰係数Ａを用いて（１７）式に対応するマトリクスを作成してもよい。
【００７２】
このようにして得られたＣＩＥ−Ｌ^*ａ^*ｂ^*色空間の値が、上述の（１）式によりＸＹＺ色空間の値へ変換され、マトリクス２０の出力信号である測色値信号１８とされる（図１参照）。
【００７３】
次に、測色値信号１８は、必要に応じて光源変更部２５により測色的に変更され、新たな測色値信号１８（同一の符号を用いる。）とされる。なお、観察光源をデジタルカメラの撮影時における撮影光源（ＣＩＥ−Ｄ５０と同一）とした場合には、この光源変更処理は不要である。
【００７４】
測色値ＸＹＺを光源変更部２５による光源変更処理後の新たな測色値ＸＹＺ（符号は、Ｘ′Ｙ′Ｚ′とする。）に変更する場合、その変換は、次の（１８）式、（１９）式、（２０）式に基づいて行うことができる。
【００７５】
Ｘ′＝ＸＸ₂／Ｘ₁ …（１８）
Ｙ′＝ＹＹ₂／Ｙ₁ …（１９）
Ｚ′＝ＺＺ₂／Ｚ₁ …（２０）
ただし、（１８）式〜（２０）式において、ＸＹＺは各々撮影光源下における測色値、Ｘ′Ｙ′Ｚ′は各々観察光源下における測色値、Ｘ₁Ｙ₁Ｚ₁は各々撮影光源の白色点、Ｘ₂Ｙ₂Ｚ₂は各々観察光源の白色点である。
【００７６】
以下、測色値信号１８は、必要に応じて光源変更処理のなされた信号であるものとする。
【００７７】
一方、スキャナ入力部から出力されたＲＧＢ信号１４は測色変換テーブル２４によりＸＹＺ信号である測色値信号２２に変換される。
【００７８】
この測色変換テーブル２４は、例えば、ｃｍｙの色パッチの各濃度を１３段階規則正しく振った、合計１３×１３×１３＝２１９７個からなるカラーパッチを有するカラーリバーサル原稿を準備する。そして、このカラーリバーサル原稿を構成する前記各カラーパッチをスキャナの入力部で読み取るとともに、測色計で読み取る。そして、スキャナで読み取ったＲＧＢ値と測色計で読み取ったＸＹＺ値との対応関係を求めてテーブルとしたものである。このテーブルに存在しない読取値間の値は補間処理により求める。
【００７９】
そして、スキャナの入力部の出力信号であるＲＧＢ信号１４を測色変換テーブル２４により測色値信号２２に変換すること、およびデジタルカメラの出力信号であるＲＧＢ信号１２を測色変換マトリクス２０により測色値信号１８に変換することにより、次に説明する測色的セットアップ部２６を共通に使用して、いわゆるオートセットアップ処理を行うことができる。すなわち、測色値信号１８、２２を、一旦、測色的ＥＮＤに変換することにより、オートセットアップ処理ソフトウエアを共通に使用することができるという利点がある。
【００８０】
次に、測色的セットアップ部２６について図３のフローチャートをも参照して説明する。
【００８１】
従来から、階調変換は、直感的な理解がし易いために濃度空間で行われており、また、現在、市販されているカラースキャナにおける階調変換処理および色修正処理等の画像処理も濃度信号を基準として実施されている。そこで、まず、測色値信号１８または測色値信号２２を濃度信号の表色系であるｃｍｙのＥＮＤ（等価中性濃度）信号５２に変換する（ステップＳ１１）。この変換処理は、図３に示すように、ＸＹＺ表色系からＲＧＢ表色系への線形変換処理（ステップＳ１１ａ）とＲＧＢ表色系からｃｍｙ表色系への非線形変換処理（ステップＳ１１ｂ）の２段階の処理になっている。なお、測色値信号１８または測色値信号２２から変換されたＥＮＤ信号を測色的ＥＮＤ信号５２と呼ぶ。
【００８２】
例えば、リバーサル原稿の測色値信号２２を測色的ＥＮＤ信号５２に変換する場合、図４に示すＣＩＥの色度図上、リバーサルフイルムの色再現領域７１（ハッチングで示す領域）を描き、その色再現領域７１を包含する領域７３の３つの原刺激ＲＧＢ（それぞれＲｘｙｚ、Ｇｘｙｚ、Ｂｘｙｚとする。）の各点の色度図上の座標、すなわち色度座標が（２１）式〜（２３）式に示す座標であるものとする。この場合、リバーサルフイルムの色再現域７１を包含する領域７３は、図４に示す色度図上、原刺激Ｒｘｙｚ、Ｇｘｙｚ、Ｂｘｙｚを頂点とする三角形の領域である。
【００８３】
Ｒｘｙｚ＝Ｒｘｙｚ（ｘ_R，ｙ_R，ｚ_R） …（２１）
Ｇｘｙｚ＝Ｇｘｙｚ（ｘ_G，ｙ_G，ｚ_G） …（２２）
Ｂｘｙｚ＝Ｂｘｙｚ（ｘ_B，ｙ_B，ｚ_B） …（２３）
また、色度図上、ＸＹＺ表色系の基礎刺激（白色刺激）Ｗｘｙｚの座標を（２４）式に示す。
【００８４】
Ｗｘｙｚ＝Ｗｘｙｚ（ｘ_W，ｙ_W，ｚ_W） …（２４）
この場合、次の（２５）式によりＸＹＺ表色系の測色値信号２２（右辺の右側のマトリクス）は、変換マトリクス（右辺の左側のマトリクス）を介してＲＧＢ表色系の色信号ＲＧＢ（左辺のマトリクス）に変換することができる（ステップＳ１１ａ）。
【００８５】
【数５】

【００８６】
余因数Ａ_ijは、下式で求められる。
【００８７】
Ａ_ij＝（−１）^i+jＤ_ij
ただし、Ｄ_ijは、ｉ行ｊ列を除いた小行列式である。
測色的ＥＮＤ信号５２のｃｍｙ値は、基礎刺激（白色刺激）Ｗ_XYZの座標を（２５）式に代入して求めたＲ、Ｇ、Ｂの値Ｒｗ、Ｇｗ、Ｂｗに対する（２５）式で求めたＲ、Ｇ、Ｂの値の各比率を（２６）式〜（２８）式に示すように各々対数変換することにより求めることができる。
【００８８】
ｃ＝−ｌｏｇ（Ｒ／Ｒｗ） …（２６）
ｍ＝−ｌｏｇ（Ｇ／Ｇｗ） …（２７）
ｙ＝−ｌｏｇ（Ｂ／Ｂｗ） …（２８）
ラフデータ作成部５４は、セットアップ処理を短時間に行うため、例えば、原稿１枚分の測色的ＥＮＤ信号（実際には、デジタルデータ）５２を対象とするのではなく、原稿中、オペレータからＣＰＵを通じて領域指定された画像の存在する部分のみのデータを選択すること、あるいは画像が原稿の全面に存在する場合にはデータを間引いて作成すること等、いわゆるラフスキャン処理を行う。
【００８９】
次に、オートセットアップ部５８において、ラフデータ作成部５４により選択されたラフデータである測色的ＥＮＤ信号５６にもとづいてオートセットアップ処理を行う（ステップＳ１２）。このオートセットアップ処理では、例えば、特開平２−１０５６７６号公報にも公知のように、測色的ＥＮＤ信号５６についてヒストグラムを作成した後、累積ヒストグラムを作成する。
【００９０】
そして、図１中、ＥＮＤ・ＥＮＤ変換部６４の動作を示す特性図（図１中に描いている図）に示すように、前記累積ヒストグラムの、例えば、０．１％点データ（ＨＬ濃度）Ｄ１に対応する測色的ＥＮＤ信号５６の値をＤＨに設定し、９８％点データ（ＳＤ濃度）Ｄ２に対応する測色的ＥＮＤ信号５６の値をＤＳに設定する（ステップＳ１２）。
【００９１】
分かり易く仮想的な値で説明すると、例えば、Ｄ１の濃度がＤ１＝１．０、Ｄ２の濃度がＤ２＝２．０であった場合、Ｄ１の濃度に対応するＤＨをＤＨ＝０．１に設定し、Ｄ２の濃度に対応するＤＳをＤＳ＝３．０に設定する。なお、実際上は、０．１％点データ（ＨＬ濃度）Ｄ１が網％の０％に対応する濃度に変換され、９８％点データ（ＳＤ濃度）Ｄ２が網％の１００％に対応する濃度に変換される。
【００９２】
このようにして設定された直線７４の式から、ＥＮＤ・ＥＮＤ変換部６４により本スキャンデータ（ラフスキャン処理の対象となったデータ）の全ての測色的ＥＮＤ信号５２を測色的ＥＮＤ信号６６に変換することができる（ステップＳ１３）。すなわち、階調特性を変換すること、言い換えれば階調補正処理を行うことができる。なお、変換式は、ハイライトポイント７５（Ｄ１，ＤＨ）およびシャドーポイント７６（Ｄ２，ＤＳ）を通る曲線の式とすることもできる。なお、ＨＬ濃度Ｄ１とＳＤ濃度Ｄ２（またはセットアップポイント７５、７６）を上述のように予め定めた一定の条件に基づいて自動的に決定することをオートセットアップというが、ＨＬ濃度Ｄ１とＳＤ濃度Ｄ２の値をマニュアルで決定することも可能であり、あるいは、オートセットアップ後、ＤＨとＤＳとをマニュアルで修正することも可能である。マニュアルにより決定する動作をマニュアルセットアップという。
【００９３】
次いで、測色的ＥＮＤ逆変換マトリクス５０により、ＥＮＤ・ＥＮＤ変換後の測色的ＥＮＤ信号６６のｃｍｙ値をＸＹＺ値である測色値信号２２に逆変換する（ステップＳ１４）。
【００９４】
この逆変換処理は、測色的ＥＮＤ信号６６のｃｍｙ値をＲＧＢ表色系のＲＧＢ値に変換する処理（ステップＳ１４ａ）と、このＲＧＢ値をＸＹＺ表色系のＸＹＺ値をとる測色値信号２２に変換する処理（ステップＳ１４ｂ）とからなる。測色的ＥＮＤ信号６６のｃｍｙ値をＲＧＢ表色系のＲＧＢ値に変換する処理は、前記の（２６）式〜（２８）式をＲＧＢについて解いた次の（２９）式〜（３１）式により行われ、ＲＧＢ値をＸＹＺ表色系のＸＹＺ値に変換する処理は、前記の（２５）式をマトリクスＸＹＺについて解いた次の（３２）式により行われる。
【００９５】
Ｒ＝Ｒｗ１０^-c …（２９）
Ｇ＝Ｇｗ１０^-m …（３０）
Ｂ＝Ｂｗ１０^-y …（３１）
【００９６】
【数６】

【００９７】
なお、上述のセットアップ処理は、リバーサル原稿の測色値信号２２を例に説明したが、同様にデジタルカメラの測色値信号１８に対しても適用することができる。
【００９８】
すなわち、一般的に、印刷製版分野においては、オリジナル原稿としてリバーサル原稿を使用する。このリバーサル原稿上に原シーンが露光され、そのリバーサル原稿上で発色している色再現域７１内を対象として印刷は仕上げられる。このことは、取り扱う画像信号の色再現域として、リバーサル原稿が発色する色再現域７１で十分であることを示唆している。
【００９９】
このように、測色的セットアップ部２６では、測色的ＸＹＺをＲＧＢ３原色の表色系へ線形変換し、基礎刺激（光源）Ｒｗ、Ｇｗ、Ｂｗに対する求めたＲ、Ｇ、Ｂの比率を対数変換することにより、ＥＮＤ値を求め、これを用いてオートセットアップする。ＥＮＤ・ＥＮＤ変更されたＥＮＤ値を逆対数変換処理することにより、変更されたＲＧＢ値を求め、逆変換マトリクスにより変更された測色値ＸＹＺを求めるようにしている。このため、既存のいわゆる濃度信号に基づいたオートセットアップ処理を行うことができるという効果が達成される。なお、セットアップ処理は、マニュアルセットアップ処理として行うことができることはいうまでもない。
【０１００】
次に、色素濃度変換部３０について図５、図６のフローチャートをも参照して説明する。
【０１０１】
この色素濃度変換部３０では、デジタルカメラにより撮影された原シーンの測色値信号（ＸＹＺ）１８をリバーサル原稿上の色素濃度信号（ｃｍｙ）２８に変換する。
【０１０２】
この変換過程は、図５に示す方法と図６に示す方法の２通りの方法が考えられる。
【０１０３】
すなわち、図５に示す方法は、原シーンの測色値信号（ＸＹＺ）１８を、後述するルックアップテーブルにより原シーンの色素濃度信号（ｃｍｙ）８１に変換し（ステップＳ２１：第１段階の処理）、この原シーンの色素濃度信号（ｃｍｙ）８１を後述するマトリクスによりリバーサル（単に、ＲＶともいう。）原稿上の色素濃度信号（ｃｍｙ）２８に変換する方法である（ステップＳ２２：第２段階の処理）。また、図６に示す方法は、原シーンの測色値信号（ＸＹＺ）１８を、後述するマトリクスによりリバーサル原稿上の測色値信号（ＸＹＺ）８２に変換し（ステップＳ３１：第１段階の処理）、このリバーサル原稿上の測色値信号（ＸＹＺ）８２を後述するルックアップテーブルによりリバーサル原稿上の色素濃度信号（ｃｍｙ）２８に変換する（ステップＳ３２：第２段階の処理）方法である。
【０１０４】
まず、図５に示す方法の詳細について、図７のフローチャートをも参照して説明する。図５中、ステップＳ２１の処理に係る測色値ＸＹＺを色素濃度ｃｍｙに変換するためには、色空間を概ね均等に網羅したカラーチャート、例えば、上述のマクベスカラーチェッカー７０の２４色の各測色値ＸＹＺを測色計により求める（ステップＳ２１ａ：図７参照）。
【０１０５】
次に、求めた測色値ＸＹＺから色素濃度値ｃｍｙを算出し（ステップＳ２１ｂ）、測色値ＸＹＺに対する色素濃度値ｃｍｙのルックアップテーブルを作成する。このルックアップテーブルを色素濃度変換部３０における第１段階の処理であるステップＳ２１のルックアップテーブルとして色素濃度変換部３０に設定記憶しておく。なお、ステップＳ２１の処理では、このルックアップテーブルを用いるとともに、上記２４色の測色値ＸＹＺ以外の色については、その２４色の測色値ＸＹＺを用いて補間法により色素濃度値ｃｍｙを算出することができる。
【０１０６】
ステップＳ２１ｂの算出処理、すなわち測色値（三刺激値）ＸＹＺから色素濃度ｃｍｙを求める処理を説明する。
【０１０７】
測色値ＸＹＺと色素濃度ｃｍｙとの間では、次の（３３）式〜（３６）式が成立する。
【０１０８】
Ｘ＝ｋ∫ｖｉｓＰ（λ）Ｔ（λ）ｘ（λ）ｄλ …（３３）
Ｙ＝ｋ∫ｖｉｓＰ（λ）Ｔ（λ）ｙ（λ）ｄλ …（３４）
Ｚ＝ｋ∫ｖｉｓＰ（λ）Ｔ（λ）ｚ（λ）ｄλ …（３５）
Ｔ（λ）＝１０^-h …（３６）
ただし、
ｈ＝｛ｃＤｃ（λ）＋ｍＤｍ（λ）＋ｙＤｙ（λ）＋ｂａｓｅ（λ）｝
ｋ＝１００／∫ｖｉｓＰ（λ）ｙ（λ）ｄλ（λは光の波長）
∫ｖｉｓ：可視波長域（３８０ｎｍ〜７８０ｎｍ）での定積分
Ｐ（λ）：観察光源の分光特性データ
Ｔ（λ）：透過物体の色素の分光透過率データ
｛（３６）式はランベルト・ベールの法則が成立すると仮定｝
ｘ（λ）、ｙ（λ）、ｚ（λ）：等色関数
Ｄｃ（λ）、Ｄｍ（λ）、Ｄｙ（λ）：ｃｍｙ色素の分光濃度データ
ｂａｓｅ（λ）：フイルムベースの分光濃度データ
これら（３３）式〜（３６）式から色素濃度ｃｍｙを求める場合には、逆関数を求めればよいのであるが、直接求めることができない。そこで、公知のニュートン（ニュートン・ラフソン）法（例えば、太田登著「色彩工学」ｐｐ２５４−２６０東京電機大学出版局１９９３年１２月２０日第１版１刷発行）のような逐次近似法を用いればよい。前記参考書を引用して、このニュートン・ラフソン法（Ｎ・Ｒ法と略記する。）について簡単に説明する。
【０１０９】
一般的な方程式ｙ＝ｆ（ｘ）をｆ（ｘ）＝０の根に近いｘ＝ｘ０でテーラ展開して１次の項のみを取ると、ｘの微小変化Δｘに対して、ｆ（ｘ０＋Δｘ）＝ｆ（ｘ０）＋ｆ′（ｘ０）・Δｘが成立する。ただし、ｆ′（ｘ０）は、ｆ（ｘ）の微分係数ｆ′（ｘ）にｘ＝ｘ０を代入したものである。したがって、ｆ（ｘ）＝０のより正しい値ｘ１は、ｆ（ｘ０＋Δｘ）＝０として、ｘ１＝ｘ０＋Δｘ＝ｘ０−ｆ（ｘ０）／ｆ′（ｘ０）で求められる。これは、図８に示すように、関数ｙ＝ｆ（ｘ）上の点（ｘ０，ｙ０）で接線８３を引き、その接線８３とｘ軸との交点ｘ１を求めることに相当する。
【０１１０】
これを（３３）式〜（３６）式に適用するためには、（３６）式を（３３）式〜（３５）式に代入した後、ある関数ｆｘ、ｆｙ、ｆｚを用いれば、（３３）式〜（３５）式は、（３７）式〜（３９）式のように表すことができる。
【０１１１】
Ｘ＝ｆｘ（ｃ，ｍ，ｙ） …（３７）
Ｙ＝ｆｙ（ｃ，ｍ，ｙ） …（３８）
Ｚ＝ｆｚ（ｃ，ｍ，ｙ） …（３９）
これら（３７）式〜（３９）式において、初期値をｃ０，ｍ０，ｙ０とし、そのときの三刺激値をＸ０，Ｙ０，Ｚ０とする。今、ｃ０，ｍ０，ｙ０に微小変化Δｃ，Δｍ，Δｙを加えたとき、三刺激値がΔＸ，ΔＹ，ΔＺだけ変化したとすれば、次の（４０）式が得られる。
【０１１２】

ただし、例えば、∂ｆｘ／∂ｃは、関数ｆｘのｃに関する偏微分係数を表す。
【０１１３】
この（４０）式を整理すれば、（４１）式が得られる。同様に、ΔＹ、ΔＺも（４２）式、（４３）式に示すように得られる。
【０１１４】
ΔＸ＝Δｃ・∂Ｘ／∂ｃ＋Δｍ・∂Ｘ／∂ｍ＋Δｙ・∂Ｘ／∂ｙ …（４１）
ΔＹ＝Δｃ・∂Ｙ／∂ｃ＋Δｍ・∂Ｙ／∂ｍ＋Δｙ・∂Ｙ／∂ｙ …（４２）
ΔＺ＝Δｃ・∂Ｚ／∂ｃ＋Δｍ・∂Ｚ／∂ｍ＋Δｙ・∂Ｚ／∂ｙ …（４３）
（４１）式〜（４３）式を（４４）式に示すように行列表示する。
【０１１５】
（Ｑ）＝（Ｊ）（Ｐ） …（４４）
ただし、（Ｑ）は、要素が１行目から順にΔＸ，ΔＹ，ΔＺである１行３列のマトリクス、（Ｊ）は偏微分係数が３行３列からなるヤコビアンマトリクス、（Ｐ）は、要素が１行目から順にΔｃ，Δｍ，Δｙである１行３列のマトリクスである。
【０１１６】
（４４）式の両辺に、ヤコビアンマトリクス（Ｊ）の逆マトリクス（Ｊ）^-1を掛ければ、（４５）式が得られる。
【０１１７】
（Ｐ）＝（Ｊ）^-1（Ｑ） …（４５）
したがって、初期値ｃ０，ｍ０，ｙ０を（４６）式のようにそれぞれｃ１，ｍ１，ｙ１に修正すれば、より正しい近似値を得ることができる。
【０１１８】
ｃ１＝ｃ０＋Δｃ
ｍ１＝ｍ０＋Δｍ
ｙ１＝ｙ０＋Δｙ …（４６）
以上のようにして得られるヤコビアンマトリクス（Ｊ）を用いて繰り返し計算を行うことにより、任意の目標値である測色値ＸＹＺに対する色素濃度信号ｃｍｙを求めることができる。同様の処理をＸＹＺ色空間における格子上の全目標値に対して行うことで、原シーンに対する測色値信号（ＸＹＺ）１８を原シーンに対する色素濃度信号（ｃｍｙ）８１に変換するための逆変換テーブルを生成し、これをルックアップテーブルとして色素濃度変換部３０に第１段階の処理用として保持しておく（ステップＳ２１）。
【０１１９】
次に、この第１段階の処理に続いて行われるステップＳ２２の第２段階処理用のマトリクスの作成処理について説明する。なお、このマトリクスの作成処理は、図２を参照して既に説明した重回帰分析処理を含む処理であり、簡潔に説明する。
【０１２０】
まず、ステップＳ２１ｂで算出した色素濃度ｃｍｙについて、ｃ、ｍ、ｙ、ｃ²、ｍ²、ｙ²、ｃｍ、ｍｙ、ｙｃの９変数２４色のデータを算出する（ステップＳ２２ａ）。
【０１２１】
次いで、この９変数２４色のデータを主成分分析し、９主成分の主成分スコアを求める（ステップＳ２２ｂ）。
【０１２２】
一方、リバーサルフイルムに露光されたマクベスカラーチェッカー２４色についての測色値を測色計により求める（ステップＳ２２ｃ）。
【０１２３】
次に、上述したＮ・Ｒ法によりリバーサルフイルム上の色素濃度ｃｍｙ（ＲＶｃｍｙとも記載する。）を求める（ステップＳ２２ｄ）。
【０１２４】
次に、リバーサルフイルム上の色素濃度ＲＶｃｍｙの各々（ＲＶｃ、ＲＶｍ、ＲＶｙ）を目的変数とし、ステップＳ２２ｂで求めた主成分スコア（定数項を含む）を説明変数として、重回帰分析処理により３×１０の偏回帰係数マトリクスを求める（ステップＳ２２ｅ）。
【０１２５】
なお、重回帰分析処理を行う際に、上述したように、目的変数群となっている２４色に１：１に対応した重み付けマトリクスを作用させてもよい（ステップＳ２２ｆ）。
【０１２６】
次に、各目的変数ＲＶｃ、ＲＶｍ、ＲＶｙについて求めた偏回帰係数について、例えば５％有意で検定を行い（ステップＳ２２ｇ）、５％有意な説明変数を記憶手段に格納しておくとともに、偏回帰係数を第２段階処理用のマトリクスとして、色素濃度変換部３０に格納しておく（ステップＳ２２ｈ）。なお、この場合においても、検定を行わず、全ての係数を用いてもよい。
【０１２７】
このようにして、原シーンの測色値信号（ＸＹＺ）１８をリバーサルフイルム上の色素濃度信号ｃｍｙに変換する色素濃度変換部３０が構築される。
【０１２８】
次に、図６に示す方法の詳細について、図９のフローチャートをも参照して説明する。なお、図９に示す処理も、図２を参照して既に説明した重回帰分析処理を含む処理であり、簡潔に説明する。
【０１２９】
まず、マクベスカラーチェッカー７０の２４色の各測色値ＸＹＺを測色計により求める（ステップＳ３１ａ：図９参照）。
【０１３０】
次に、求めた測色値ＸＹＺについて、Ｘ、Ｙ、Ｚ、Ｘ²、Ｙ²、Ｚ²、ＸＹ、ＹＺ、ＺＸの９変数２４色のデータを算出する（ステップＳ３１ｂ）。
【０１３１】
次いで、この９変数２４色のデータを主成分分析し、９主成分の主成分スコアを求める（ステップＳ３１ｃ）。
【０１３２】
一方、リバーサルフイルムに露光されたマクベスカラーチェッカー２４色についての測色値ＸＹＺを測色計により求める（ステップＳ３１ｄ）。なお、リバーサルフイルムに露光されたマクベスカラーチェッカーとは、カラーリバーサルフイルムが装着された光学式カメラにより、所定の光源下で前記マクベスカラーチェッカーを撮影してリバーサルフイルムを露光した後、この露光されたリバーサルフイルムを現像して得られるリバーサル原稿をいう。
【０１３３】
次に、リバーサル原稿上の測色値ＸＹＺの各々を目的変数（ＲＶＸ、ＲＶＹ、ＲＶＺ）とし、ステップＳ３１ｃで求めた主成分スコア（定数項を含む）を説明変数として、重回帰分析処理により３×１０の偏回帰係数マトリクスを求める（ステップＳ３１ｅ）。なお、重回帰分析処理を行う際に、上述したように、目的変数群となっている２４色に１：１に対応した重み付けマトリクスを作用させてもよい（ステップＳ３１ｆ）。
【０１３４】
次に、各目的変数ＲＶＸ、ＲＶＹ、ＲＶＺについて求めた偏回帰係数の例えば５％有意で検定を行い（ステップＳ３１ｇ）、５％有意な説明変数を記憶手段に格納しておくとともに、偏回帰係数を第１段階の処理用マトリクスとして、色素濃度変換部３０に格納しておく（ステップＳ３１ｈ、ステップＳ３１）。この場合においても、検定を行わず、全ての係数を用いてもよい。
【０１３５】
次いで、ステップＳ３２の第２段階の処理用テーブルは、図７を参照して説明したステップＳ２１ｂあるいはステップＳ２２ｄの方法と同一の方法で作成されるので、その説明を省略する。
【０１３６】
このようにして、原シーンの測色値信号（ＸＹＺ）１８をリバーサル原稿上の色素濃度信号ｃｍｙに変換する図６の処理に基づく色素濃度変換部３０が構築される。
【０１３７】
次に、ｃｍｙｋ変換部３２を構成する原シーン忠実再現テーブル４２の作成処理について、図１０に示すフローチャートに基づいて説明する。
【０１３８】
原シーン忠実再現テーブル４２は、測色値信号（ＸＹＺ）１８を測色的に保存された各網％データであるｃｍｙｋ信号（ｃｍｙｋデータ）４６に変換するルックアップテーブルである。
【０１３９】
なお、この原シーン忠実再現テーブル４２の作成処理は、この出願人の出願に係る特願平８−１５４５８４号明細書に記載された作成処理と同じ処理である。
【０１４０】
この場合、まず、画像出力部３５に対して、規則的な間隔からなる複数のｃｍｙｋ網％データを与えることにより、ｃｍｙｋの濃度および混合割合が段階的に変化するカラーパッチを有するｃｍｙｋカラーチャートを作成する（ステップＳ４１）。
【０１４１】
この場合、例えば、ｃｍｙｋの規則的な間隔とは、ｃｍｙｋの各色について網％で０、２０、…、１００％等、２０％きざみで増加させた６段階の間隔とすれば、カラーパッチの総数は、４⁶＝１２９６になる。
【０１４２】
次いで、画像出力部３５により作成されたｃｍｙｋカラーチャートの各パッチを測色計を用いて測色し（ステップＳ４２）、この測色データから測色値（刺激値）ＸＹＺを求め、前記ｃｍｙｋ網％データから測色値ＸＹＺデータへの変換テーブル（順変換テーブルという。）を作成する（ステップＳ４３）。
【０１４３】
なお、この順変換テーブルは、補間用のテーブルとしても使用するので、前記規則的な間隔は、補間の正確性を考慮した場合には細かいほどよいが、その分、測色作業が膨大となるので、その測色作業の繁雑さと前記正確性と以下に説明するコンピュータの処理時間との、いわゆるトレードオフで決定すればよい。
【０１４４】
ところで、前記原シーン忠実再現テーブル４２において、入力される任意の測色値信号（測色値ＸＹＺデータ、刺激値データＸＹＺ、または単にＸＹＺともいう。）１８からこれに対応するｃｍｙｋ信号（ｃｍｙｋ網％データ、色データｃｍｙｋ、または単にｃｍｙｋともいう。）４６を求めることは、変数が３変数から４変数と増加するため、１つの測色値ＸＹＺデータ１８に対してｃｍｙｋ網％データ４６の複数の解が存在する可能性がある。この問題を解決するためには、３変数同士の関係とする必要があり、ここでは、ｃｍｙｋ網％データのうち、色データｋについて画像出力部３５が採りうる最大値Ｋｍａｘ（ｋ＝Ｋｍａｘ）に固定する（ステップＳ４４）。最大値Ｋｍａｘとは、例えば、ｃｍｙｋ網％データのｋの値が１００％である。
【０１４５】
そして、３変数である任意のＸＹＺの値を３変数である対応するｃｍｙ（ｋは固定）の値に変換する（ステップＳ４５）。
【０１４６】
この場合、任意のＸＹＺの値である目標値Ｘ０、Ｙ０、Ｚ０に対するｃｍｙｋ（ここでは、ｋはｋ＝Ｋｍａｘに固定しているのでｋは定数であり、その意味でｃｍｙの３変数）の値である、ｃ０、ｍ０、ｙ０、ｋ０＝Ｋｍａｘを求める場合には、ｋ０＝Ｋｍａｘであるときの上記順変換テーブルを利用して、回帰式の偏回帰係数を求める。
【０１４７】
このときの回帰式は、係数項となる３行４列のマトリクスをＡ、測色値ＸＹＺの１行３列のマトリクスをＴ、ｋを固定した、換言すれば定数と考えて、これを含めた残りのｃｍｙの１行４列のマトリクスをＤと表せば、次の（４７）式のように表すことができる。
【０１４８】
Ｔ＝ＡＤ …（４７）
この（４７）式は、次の（４８）式の関係を表している。
【０１４９】
【数７】

【０１５０】
（４７）式、（４８）式において、マトリクスＤを構成する「１」は、ｃｍｙの３次元平面式における定数項を与えるために設定した値である。
【０１５１】
この（４８）式における係数Ａは、ｋ＝Ｋｍａｘのときに得られている上述の順変換テーブルの各データセットを代入することで、上述した重回帰分析により求めることができる。
【０１５２】
次に、ｋ＝Ｋｍａｘのときに求められた回帰式を用いて、上述のＮ・Ｒ法により、目標値Ｘ０、Ｙ０、Ｚ０に対応するｃｍｙ（ｋはＫｍａｘとして設定されている。）の値を求めることができる。
【０１５３】
次に、求められたｃｍｙの値が画像出力部３５の再現範囲内の値（色）であるかどうかの判断を行う（ステップＳ４６）。すなわち、網％データ（色データ）ｃｍｙの再現可能な最小濃度をＣｍｉｎ、Ｍｍｉｎ、Ｙｍｉｎ、最大濃度をＣｍａｘ、Ｍｍａｘ、Ｙｍａｘとして求められたデータｃｍｙが次の（４９）式の関係を満足するかどうかを判定する。
【０１５４】
Ｃｍｉｎ≦ｃ≦Ｃｍａｘ
Ｍｍｉｎ≦ｍ≦Ｍｍａｘ
Ｙｍｉｎ≦ｙ≦Ｙｍａｘ …（４９）
求めたｃｍｙの各値が、（４９）式を満足している場合には、刺激値データＸＹＺの目標値（Ｘ０，Ｙ０，Ｚ０）に対する色データｃｍｙｋを各々ｃ＝ｃｓｏｌ、ｍ＝ｍｓｏｌ、ｙ＝ｙｓｏｌ、ｋ＝ｋｓｏｌ（この場合、ｋｓｏｌ＝Ｋｍａｘ）として逆変換テーブルを作成し、これを原シーン忠実再現テーブル４２として設定する（ステップＳ４７）。以下の説明において、必要に応じて（ｃｓｏｌ，ｍｓｏｌ，ｙｓｏｌ，ｋｓｏｌ）のデータセットをｃｍｙｋｓｏｌという。
【０１５５】
ステップＳ４５の過程で求めたｃｍｙデータの値が（４９）式を満足しない値であった場合には、ｋ＝Ｋｍａｘに固定してある色データｋをｋ＝ｋ−Δｋ、この場合ｋ＝Ｋｍａｘ−Δｋとした後（ステップＳ４８）、色データｋが所定の最小値ｋ＝Ｋｍｉｎよりも小さくならない範囲で（ステップＳ４９）、ステップＳ４５の処理を繰り返す。なお、微小変化分Δｋは、第１逆変換テーブルを構成する色データｋの任意のデータ間隔であり、例えば、色データｋが０％〜１００％の範囲のデータとして設定される場合には、１％きざみ等に設定すればよく、色データｋが０〜２５５の範囲のデータとして設定されるのであれば、値１きざみに設定すればよい。ステップＳ４５の２回目の処理を行う場合には、色データｋ＝Ｋｍａｘ−１＝１００−１＝９９に対する、（４８）式の左辺の測色値であるＸＹＺ値は、ｋ＝Ｋｍａｘ＝１００％とｋ＝８０％の値で測色され順変換テーブルとして格納されている各ＸＹＺデータの補間処理（内挿処理）により求めることができる。
【０１５６】
一方、ステップＳ４９の処理において、Ｋｍｉｎ＞ｋと判定された場合には、目標値（Ｘ０，Ｙ０，Ｚ０）に対する色データｃｍｙｋを画像出力部３５の色再現域外のデータに指定し、ここでは、色データｃｍｙｋｓｏｌの算出は行わない（ステップＳ５０）。
【０１５７】
以上の処理を原シーン忠実再現テーブル４２に供給される全ての刺激値データＸＹＺ１８を目標値（Ｘ０，Ｙ０，Ｚ０）として行うことにより、画像出力部３５の色再現域内の色データである網％データｃｍｙｋ４６を得ることのできる刺激値データＸＹＺ１８に対して、最大となる色データｋを与えた場合の色データｃｍｙｋｓｏｌを求めることができる（ステップＳ５１）。
【０１５８】
なお、ステップＳ４４における色データｋの固定は、最大値Ｋｍａｘではなく、その最小値Ｋｍｉｎ（ｋ＝Ｋｍｉｎ）に固定してもよく、その場合には、ステップＳ４８の処理をｋ＝ｋ＋Δｋとし、ステップＳ４９の処理をＫｍａｘ＜ｋとすればよい。また、色データｋを任意の値とすることも容易に考えられる。例えば、ステップＳ４８の処理は、任意の値に設定したｋについて、ｋ＝ｋ−Δｋとｋ＝ｋ＋Δｋを交互に行えばよく、ステップＳ４９の処理は、ステップＳ４８の処理に対応させて行えばよい。
【０１５９】
さらに、ステップＳ４９、Ｓ５０において色再現域外データに指定された色データｃｍｙｋは、この発明の要旨ではないので詳しくは説明しないが、前記特願平８−１５４５８４号に記載されているように、いわゆるガマットマッピング（Gamat Mapping ）処理による色データＣＭＹＫの圧縮処理やクリッピング処理により、目標値（Ｘ０，Ｙ０，Ｚ０）に対する色データｃｍｙｋの逆変換テーブルを作成することができる。
【０１６０】
このようにして、測色値信号１８を構成する全ての目標値（Ｘ０，Ｙ０，Ｚ０）に対して求められたｃｍｙｋ信号４６の対応表を原シーン忠実再現テーブル４２として格納しておくことで、任意の測色値信号１８を画像出力部３５の色再現範囲内のｃｍｙｋ信号４６に変換することができる。
【０１６１】
次に、ｃｍｙｋ変換部３２を構成する標準条件再現テーブル４４の作成について説明する。
【０１６２】
この標準条件再現テーブル４４は、図１に示すように、デジタルカメラにより取り込まれたＲＧＢ信号１２がマトリクス２０により測色値信号（ＸＹＺ）１８に変換され、この測色値信号ＸＹＺが色素濃度変換部３０により色素濃度信号（ｃｍｙ）２８に変換され、この変換後の色素濃度信号（ｃｍｙ）２８を網％信号であるｃｍｙｋ信号４８に変換するためのテーブルである。
【０１６３】
また、標準条件再現テーブル４４は、スキャナにより取り込まれたＲＧＢ信号１４がテーブル２４により測色値信号（ＸＹＺ）２２に変換され、この変換後の測色値信号（ＸＹＺ）２２が色素濃度変換部３６により色素濃度信号（ｃｍｙ）３４に変換され、この変換後の色素濃度信号（ｃｍｙ）３４を網％信号であるｃｍｙｋ信号に変換するためのテーブルである。この場合、色素濃度変換部３６は、前述の色素濃度変換部３０の作成過程で示した図７中、ステップＳ２１ｂ、Ｓ２２ｄの方法と同様の過程で作成することができるので、その作成過程についての説明は省略する。
【０１６４】
標準条件再現テーブル４４の作成処理を、図１１に示すフローチャートを参照して説明する。
【０１６５】
まず、リバーサル原稿上に色素濃度ｃｍｙが格子状に配列された１３×１３×１３個のカラーパッチを有するカラーチャートを作成する（ステップＳ６１）。このカラーチャートは、ｃ（シアン）ｍ（マゼンタ）ｙ（黄）が、それぞれ最小濃度から最大濃度まで１３段階振られた各カラーパッチを有するものである。
【０１６６】
次に、このカラーチャートをスキャナのデフォルト分解条件で分解する、言い換えれば、透過原稿であるカラーチャートをスキャナで読み取り、デジタルデータ化する（ステップＳ６２）。スキャナのデフォルト分解条件は、少なくとも、階調特性変換処理、色修正処理、ＵＣＲ（下色除去）処理およびＫ版生成処理を含んでいる。ただし、分解するとき、カラーパッチのグレーの最小濃度を、ｃｍｙｋ網％が全て０％となるようにし、かつカラーパッチのグレーの最大濃度を、いわゆるベタとなるようなｃｍｙｋ網％となるように設定する。
【０１６７】
そして、各カラーパッチの色素濃度とスキャナのデフォルト分解条件で分解したときのｃｍｙｋ網％値を対応させて変換テーブル（ルックアップテーブル）を作成し（ステップＳ６３）、この変換テーブルを標準条件再現テーブル４４として設定する。なお、実際には、この標準条件再現テーブル４４と補間処理により、最小濃度から最大濃度までの任意の色素濃度信号（ｃｍｙ）２８または色素濃度信号（ｃｍｙ）３４を所望のｃｍｙｋ網％信号４８に変換することができる。
【０１６８】
以上詳述したように、上述の実施の形態によれば、第１の色空間の任意の色信号、例えば、デジタルカメラで撮影したＲＧＢ信号１２を、第２の色空間の色信号、例えば、測色値信号（ＸＹＺまたはＬ^*ａ^*ｂ^*）１８に変換する際に、まず、所望の代表的な色を彩度、明度、色相に応じて変化させた複数のカラーパッチを有するカラーチャート、例えば、市販のマクベスカラーチェッカーを前記第１の色空間の色信号を出力するデバイスであるデジタルカメラで撮像する。そして、撮像し読み取ったカラーチャートのＲＧＢ信号の、少なくとも一次項、二次項、クロス項を算出し、これらを第１の説明変数群とする。そして、前記第１の説明変数群を主成分分析し、相互の説明変数が直交する第２の説明変数群に変換する。次いで、前記カラーチャートの各カラーパッチの各測色値（前記第２の色空間の色信号）を目的変数群とし、この目的変数群と第２の説明変数群を重回帰分析して、偏回帰係数のマトリクスを求める。最後に、前記マトリクスを構成する偏回帰係数の検定を行い、前記第２の説明変数群の中、所定程度有意な説明変数を選択し、この選択した説明変数を使用可能な説明変数とする。
【０１６９】
このようにして得られた、使用可能な説明変数に係る偏回帰係数のマトリクスを前記ＲＧＢ信号１２に作用させて、前記測色値信号（ＸＹＺ）１８に変換するようにしている。
【０１７０】
したがって、この実施の形態によれば、色信号を変換する方法として、データ量の少ないマトリクスによる方法であるので、このマトリクス等を記憶するためのメモリの記憶容量を少なくすることができる。また、主成分分析した第２の説明変数群を用いるようにしているので、従来の技術の項で説明したマトリクス方式に比べて色変換精度が優れている。さらに、目的変数群を構成する測色値Ｃのうち、原稿のシーン種（人物の肌色、空の青色、植物の緑色）等に対応して所望の色に対する重みを大きくすることで、その所望の色に限定して色精度よく変換することができる。さらにまた、その所望の色は、デジタルカメラから供給される１画面についてのＲＧＢ信号を統計処理して得られる、最頻値または平均値等の色とすることで、この１画面の原稿全体を最も精度よく色変換できる可能性が高くなる。さらにまた、カラーチャートとして、色空間を概ね均等に網羅する２４色のカラーパッチを有するマクベスカラーチェッカーを用いることで、ＲＧＢ信号を精度よく測色値に色変換するための変換マトリクスを簡易に作成することができる。
【０１７１】
なお、この発明は上述の実施の形態に限らず、この発明の要旨を逸脱することなく種々の構成を採り得ることはもちろんである。
【０１７２】
【発明の効果】
以上説明したように、この発明によれば、改良されたマトリクス方式を採用しているので、少ないメモリ容量で精度よく色変換処理を行うことができるという効果が達成される。
【０１７３】
また、画像の所望の色に限定して精度よく色変換することも可能である。
【０１７４】
さらに、必要に応じて、既成のカラーチャートを使用することができる。例えば、マクベスカラーチェッカー（登録商標）をカラーチャートとして用いた場合には、２４色分のパッチで色空間全体に対して色変換精度のよいマトリクスを作成することができるという利点が得られる。
【図面の簡単な説明】
【図１】この発明の一実施の形態が適用された画像信号処理装置の全体的な構成を示すブロック図である。
【図２】測色変換マトリクスの作成処理の説明に供されるフローチャートである。
【図３】測色的セットアップ部の説明に供されるフローチャートである。
【図４】ＣＩＥ色度図上でのリバーサルフイルムの色再現域等の説明に供される色度図である。
【図５】デジタルカメラの測色信号をリバーサル原稿上の色素濃度信号に変換する変換過程の説明に供されるフローチャートである。
【図６】デジタルカメラの測色信号をリバーサル原稿上の色素濃度信号に変換する他の変換過程の説明に供されるフローチャートである。
【図７】色素濃度変換マトリクスの作成処理の説明に供されるフローチャートである。
【図８】ニュートン・ラフソン法の説明に供される線図である。
【図９】色素濃度変換マトリクスの作成処理の一部の説明に供されるフローチャートである。
【図１０】原シーン忠実再現テーブルの作成処理の説明に供されるフローチャートである。
【図１１】標準条件再現テーブルの作成処理の説明に供されるフローチャートである。
【符号の説明】
１０…画像信号処理装置１２、１４…ＲＧＢ信号
１６…測色変換部１８、２２…測色値信号（ＸＹＺ）
２０…測色変換マトリクス２４…測色変換テーブル
２５…光源変更部２６…測色的セットアップ部
２８、３４…色素濃度信号（ｃｍｙ）３０、３６…色素濃度変換部
３２…ｃｍｙｋ変換部３５…画像出力部
４２…原シーン忠実再現テーブル４４…標準条件再現テーブル
４６、４８…ｃｍｙｋ信号６４…ＥＮＤ・ＥＮＤ変換部

Claims

ＲＧＢ空間の任意の色信号をＸＹＺまたはＬ ^* ａ ^* ｂ ^* の測色値空間の色信号に変換する色変換方法において、
所望の代表的な色を彩度、明度、色相に応じて変化させた複数のカラーパッチを有するカラーチャートを前記ＲＧＢ空間の色信号を出力するデバイスで読み取る過程と、
読み取ったＲＧＢ空間の色信号を輝度値又は濃度値に変換し、変換した前記輝度値又は前記濃度値の少なくとも一次項の輝度値又は濃度値を第１の説明変数群として求める過程と、
前記第１の説明変数群を主成分分析して求めた前記輝度値又は前記濃度値の主成分スコアを相互の説明変数が直交する第２の説明変数群とする過程と、
測色計により求めた前記カラーチャートの各カラーパッチの色についてのＸＹＺまたはＬ ^* ａ ^* ｂ ^* の各測色値を目的変数群とし、この目的変数群と前記第２の説明変数群を重回帰分析して、偏回帰係数のマトリクスを求める過程と、
前記偏回帰係数の検定を行い、前記第２の説明変数群の中、所定程度有意な説明変数を選択し、この選択した説明変数を使用可能な説明変数とし、この使用可能な説明変数に対応する前記偏回帰係数のマトリクスを所定程度有意な偏回帰係数のマトリクスとする過程とを有し、
前記カラーチャートの各カラーパッチの各測色値を中心とした領域に色空間を分割し、かつ前記ＲＧＢ空間のＲＧＢ信号を、一旦、重み付けのない前記偏回帰係数のマトリクスにより前記ＸＹＺまたは前記Ｌ ^* ａ ^* ｂ ^* の測色値信号に変換し、前記分割した各領域に存在する前記測色値信号の画素の頻度値を求め、この頻度値に応じて、前記カラーチャート中、所望の特定色に対して重み付けを行い、前記所望の特定色の色変換精度を向上させ、
前記ＲＧＢ空間の任意の色信号に対して前記使用可能な説明変数に対応する前記所定程度有意な偏回帰係数のマトリクスを作用させて、前記測色値空間の色信号に変換することを特徴とする色変換方法。
請求項１記載の方法において、前記偏回帰係数のマトリクスを求める過程では、前記カラーチャート中、所望の特定色に対して重み付けを行い、前記所望の特定色の色変換精度を向上させることを特徴する色変換方法。
請求項２記載の方法において、所望の特定色は、灰色であることを特徴とする色変換方法。
請求項２記載の方法において、所望の特定色は、前記ＲＧＢ空間の色信号を統計処理して得られる色とすることを特徴とする色変換方法。
請求項４記載の方法において、前記統計処理は、前記ＲＧＢ空間の色信号中、肌色、青色、緑色を抽出して行うことを特徴とする色変換方法。
請求項４記載の方法において、前記統計処理は、前記カラーチャートの各カラーパッチの色を中心とした領域に色空間を分割し、分割した各領域に存在する、前記ＲＧＢ空間の色信号の画素の頻度値を求める処理であり、この頻度値に応じて前記所望の特定色の重み付けを行うことを特徴とする色変換方法。
請求項１〜７のいずれか１項に記載の方法において、前記カラーチャートはマクベスカラーチェッカーであることを特徴とする色変換方法。