JP4154116B2 - 色変換テーブル作成方法および色変換テーブル作成装置並びに色変換テーブル作成プログラムを記録したコンピュータ読取可能な記録媒体 - Google Patents

Description

【０００１】
【発明の属する技術分野】
本発明は、イメージスキャナ等の色入力機器から出力された色信号（例えばＲＧＢ信号）と、その色入力機器の色空間とは異なる色空間の色信号（例えばL*a*b*信号）とを対応させる色変換テーブルを作成する方法および装置、並びに、色変換テーブル作成プログラムを記録したコンピュータ読取可能な記録媒体に関する。
【０００２】
【従来の技術】
一般に、色空間とは、色を表現する座標系のことである。色空間には様々なものがあり、カラープリンタ（色出力機器）やイメージスキャナ（色入力機器）などは、それぞれ、ＣＭＹ（シアン，マゼンタ，イエロー）やＲＧＢ（赤，緑，青）といったデバイスに依存した色空間（デバイス依存色空間）を有している。色信号値の最小値および最大値がそれぞれ０および２５５であるカラープリンタの出力色を例にとれば、ＣＭＹの各色信号値が０から２５５の範囲で出力色を表現する。
【０００３】
このようなデバイスに依存した色空間のほか、デバイスに依存しない色空間（デバイス独立色空間）であるL*a*b*やＸＹＺなども広く用いられている。デバイスに依存しない色空間は、絶対的な色を表現することができるので、異なるデバイスの色空間の間で色変換を行なう際、デバイス依存色空間の間における中間的な色空間として用いられている。ここで、色変換とは、ある色空間で表現されている色信号値を、異なる色空間で表現される色信号値に変換することをいう。
【０００４】
例えば、カラー原稿の色の見えと、そのカラー原稿をカラースキャナで読み込んで得られたＲＧＢ画像をカラープリンタによって出力した印刷物の色の見えとを一致させるには、以下の手順(1),(2)で色変換を行なう。
(1)イメージスキャナにより得られたＲＧＢ画像を例えばL*a*b*形式で表現された画像に色変換する。
(2)(1)で得られたL*a*b*画像をカラープリンタ用のＣＭＹデータに変換し、そのＣＭＹデータに基づいてカラープリンタが出力を行なう。
【０００５】
上述した各手順(1),(2)における色変換に際して、色変換テーブル（プロファイル）が使用される。色変換テーブルには、異なる色空間どうしの対応関係が登録されており、色変換に際しては、これらの対応関係が用いられる。色変換テーブルに登録されている色を変換する場合には、色変換テーブルを参照し、その色について登録されている対応関係を読み出して色変換を行なうが、色変換テーブルに登録されていない色を変換する場合には、色変換テーブルに登録されているデータに基づいて補間処理を行ない、所望の色の変換を行なう。
【０００６】
イメージスキャナ等の色入力機器の入力特性を示す色変換テーブル（プロファイル）は、例えば表１を参照しながら後述するごとく、ＲＧＢ表色系における座標を示す色信号値（ＲＧＢ値）とL*a*b*色空間における色度座標値（L*a*b*値）との対を、約５０００組、集めて構成されている。
このような色変換テーブルを用いて色変換を行なう場合、色変換テーブル作成の手間を削減することや、精度良く色変換を行なえる色変換テーブルを作成することなどが重要なポイントとなる。
【０００７】
なお、本発明は、イメージスキャナ等の色入力機器の入力特性を示す色変換テーブルの作成を対象としているので、以下、イメージスキャナの色変換テーブルの作成について説明する。
イメージスキャナの色変換テーブル作成する場合、まず、対象となるイメージスキャナを用いて、ＩＳＯによって規定された標準色票を読み込み、この標準色票上の各色領域からの三原色成分についての反射光強度を示す色信号を得る。
【０００８】
ついで、測色器などを用いて、上述した標準色票上の各色領域における絶対色を示す色度座標値を計測し、これらの色度座標と該当する色領域に対応する色信号値とを対応付けることにより、色変換テーブルの基礎となる対応関係を得る。
その後、この基礎となる対応関係に基づいて補間演算処理を行なうことにより、ＲＧＢ空間において均一に分布する格子点に相当する色信号値と、この色信号値に対応すべき絶対色をL*a*b*色空間における色度座標値との対応関係を求めている。
【０００９】
つまり、標準色票上に表現された例えば２８８色について上述したような計測作業を行ない、その計測作業により得られた２８８組の対応関係に基づいて、補間演算処理を行ない、５０００組の対応関係を作成している。
上述した補間演算処理に関する技術としては、例えば特開平７−３０７８７２号公報に開示された技術がある。同公報では、非格子形状の色空間で補間を行なう方式として、距離の４乗に反比例した重みで補間を行なうシェパード方式が従来技術として挙げられ、この方式で生じるスカラッピング（scalloping：扇形切り欠き）という問題を解決している。
【００１０】
スカラッピングについて単純化した例が特開平７−３０７８７２号公報の図１に示されている。この図１に示すように、測定によりベクトルＡ，Ｂ，Ｃ，Ｄ，Ｅが得られている時に、ベクトルＡとＢとの間に存在する新しいベクトルＦは、シェパード方式で重み付けを行なった全てのベクトルＡ，Ｂ，Ｃ，Ｄ，Ｅを平均化することで得られる。その結果、補間されたベクトルＦは、このベクトルＦから少し離れているベクトルＣ，Ｄ，Ｅが寄与しているために、ベクトルＡ，Ｂに基づいて予想された点線で示すベクトルとは異なるものとなってしまう。このような現象をスカラッピングという。このようなスカラッピングは、格子間隔の狭い色変換テーブルを用いることでかなり回避することができる。
【００１１】
そこで、上記公報に開示された技術では、色変換テーブルを構成する対応関係を線形補間によって増やし、増やした対応関係を登録して、構成数の多い（間隔の短い）色変換テーブルを作成している。このようにして、実際の測定数よりも多い数の対応関係を登録された色変換テーブルを提供することができるので、変換テーブル作成の手間が削減され、且つ、スカラッピングの影響が回避されて精度良く色変換を行なうことが可能である。
【００１２】
一方、測定結果として得られた対応関係に基づいて、イメージスキャナの分光特性を推定し、得られた分光特性を利用して、色変換テーブルを構成する対応関係を得る手法も提案されている。イメージスキャナなどの色入力装置による入力特性を推定する技術としては、例えば特開平６−４６２５２号公報に開示された技術がある。同公報では、最小二乗法を用いて、測定値に基づいて導いた連立方程式を解くことにより、分光特性を推定する手法などが提案されている。
【００１３】
【発明が解決しようとする課題】
しかしながら、基本となる対応関係を補間して色変換テーブルを作成する場合（特開平７−３０７８７２号公報）においても、測定結果から推定した分光特性に基づいて色変換テーブルを作成する場合（特開平６−４６２５２号公報）においても、多数の色サンプルの色を測定する作業や、対象となる色入力機器によって多数の色サンプルを読み取る作業が必要であることには変わりはない。
【００１４】
特に、特開平７−３０７８７２号公報に開示された技術では、前述のごとく線形補間により対応関係を増やせば、スカラッピングの問題は低減される。しかし、元の対応関係に対し線形補間を行なっているため、元の対応関係が少ないと、精度良く対応関係を増やすことができない。つまり、精度良く色変換を行なえる色変換テーブルを作成するには、基礎となる対応関係をある程度の数（数百色分）だけ用意しておく必要があるが、その基礎となる対応関係を用意するには原稿を色入力機器で読み取って出力された色信号値と、その色入力機器の色空間とは異なる色空間での色信号値とを数百種類分測定する必要があるため、かなり手間がかかってしまう。
【００１５】
そこで、色入力機器の推定分光特性を用いることで、色票を色入力機器で読み込んで得られた画像色信号値（ＲＧＢ値）と分光反射率の主成分ベクトルに対する３つの重みとを１対１で対応させることを可能にして、画像色信号値から、その色票の分光反射率や共通色信号値（L*a*b*値，ＸＹＺ値等）を算出できるようにした技術（国際出願番号ＰＣＴ／ＪＰ９９／０６４４７の出願において開示された技術）も提案されている。また、この技術を用いれば、カラーパッチ数（色票上の色領域の数）が数十色と少なくても分光特性を推定することができるので、結果として、少ないカラーパッチ数で色変換テーブルを作成することができ、色変換テーブル作成の際に掛かる手間を大幅に削減することができる。
【００１６】
しかし、色票上のカラーパッチ（色領域）の持つ分光反射率は様々な形を有しているため、色空間全体に分布する多数の対応関係を求める際に、分光反射率の３つの主成分ベクトルに対する重みとして単純に１種類のものだけを用いると、つまり、全色空間について同一の色変換式を用い分光反射率を３つの主成分ベクトルの線形和として求めると、一部の領域に色変換精度が極度に悪くなる場合があり分光反射率を正確に表現しきれず、その結果として色変換テーブルの精度が悪化するという課題があった。
【００１７】
本発明は、このような課題に鑑み創案されたもので、簡易な手法で色信号の対応関係を精度良く得られるようにして、手間をかけることなく、精度良く色変換を行なえる色変換テーブルを作成することのできる、色変換テーブル作成方法および色変換テーブル作成装置並びに色変換テーブル作成プログラムを記録したコンピュータ読取可能な記録媒体を提供することを目的とする。
【００１８】
【課題を解決するための手段】
上記目的を達成すべく、本発明の色変換テーブル作成方法（請求項１）は、色入力機器から出力された、該色入力機器の色空間（以下、変換元色空間）における色信号と、該変換元色空間とは異なる色空間（以下、変換先色空間という）における色信号とを対応させる色変換テーブルを作成する方法であって、該変換先色空間の全体を分割して得られる複数の領域を、予め定義しておく分割ステップと、色票上の複数種の色領域を該色入力機器で読み込み、各色領域に対応する、該変換元色空間での色信号を該色入力機器から出力する読込ステップと、該複数種の色領域を測色器で測色し、各色領域に対応する分光反射率を該測色器から出力する測色ステップと、該分光反射率に対応する、該変換先色空間での色信号が該色変換先色空間において該複数の領域のうちのどの領域に属するかに応じて、該分光反射率を分類する分類ステップと、該読込ステップで該色入力機器から出力された該色信号と該測色ステップで該測色器から出力された該分光反射率とに基づいて、該色入力機器の分光特性を推定する分光特性推定ステップと、該分割ステップで定義された該複数の領域にそれぞれ対応した複数の色変換式を用いて、該色変換テーブルを作成する色変換テーブル作成ステップとを有し、該色変換テーブル作成ステップが、該分類ステップで分類された該分光反射率と該分光特性推定ステップで推定された前記分光特性とに基づいて、該変換先色空間での前記領域毎に該色変換式を作成しながら、各領域に応じた該色変換式を用いて、該色変換式毎に、該変換元色空間における色信号と該変換先色空間における色信号との対応関係を複数作成する対応関係作成ステップと、該対応関係作成ステップにおいて該対応関係を複数作成する際に各色変換式により該変換元色空間の一の色信号を変換した色変換結果が、各色変換式に対応した該領域に属している場合に、該色変換結果は正しいものと判断する判断ステップと、該判断ステップで正しいと判断された該色変換結果に基づいて該色変換テーブルを作成する作成処理ステップとを有するとともに、該対応関係作成ステップにおいて、該分類ステップで分類された分光反射率群毎に当該分光反射率群に対して特異値分解を施すことにより算出された複数の主成分ベクトルと、該分光特性推定ステップで推定された前記分光特性とに基づいて、該変換元色空間における色信号値を代入することにより該複数の主成分ベクトルのそれぞれについての複数の重み係数を算出する第１式を、該色変換式として該変換先色空間での前記領域毎に作成し、該変換元色空間における色信号値毎に該第１式によって算出された該複数の重み係数をそれぞれ該複数の主成分ベクトルに乗算しその線形和を求めることによって当該色信号値に対応した分光反射率を前記領域毎に算出する第２式を、該色変換式として作成し、該第２式によって前記領域毎に算出された該分光反射率を該変換先色空間における色信号に変換する第３式を、該色変換式として作成し、これらの第１式，第２式および第３式を用いて該変換元色空間の一の色信号を該変換先色空間の色信号に変換することにより、当該一の色信号について、該分類ステップで分類された複数の分光反射率群のそれぞれに対応する複数の該対応関係を作成することを特徴としている。
【００２１】
一方、本発明の色変換テーブル作成装置（請求項２）は、色入力機器から出力された、該色入力機器の色空間（以下、変換元色空間）における色信号と、該変換元色空間とは異なる色空間（以下、変換先色空間という）における色信号とを対応させる色変換テーブルを作成するものであって、色票上の複数種の色領域を該色入力機器により読み込んで得られた、各色領域に対応する、該変換元色空間での色信号を入力する入力部と、該複数種の色領域を測色し、各色領域に対応する分光反射率を得る測色器と、該分光反射率に対応する、該変換先色空間での色信号が該色変換先色空間において該複数の領域のうちのどの領域に属するかに応じて、該分光反射率を分類する分類部と、該入力部から入力された該色信号と該測色器により得られた該分光反射率とに基づいて、該色入力機器の分光特性を推定する分光特性推定部と、該変換先色空間を分割して得られる複数の領域にそれぞれ対応した複数の色変換式を用いて、該色変換テーブルを作成する色変換テーブル作成部をそなえ、該色変換テーブル作成部が、該分類部により分類された該分光反射率と該分光特性推定部により推定された前記分光特性とに基づいて、該変換先色空間での前記領域毎に該色変換式を作成しながら、各領域に応じた該色変換式を用いて、該色変換式毎に、該変換元色空間における色信号と該変換先色空間における色信号との対応関係を複数作成する対応関係作成部と、該対応関係作成部において該対応関係を複数作成する際に各色変換式により該変換元色空間の一の色信号を変換した色変換結果が、各色変換式に対応した該領域に属している場合に、該色変換結果は正しいものと判断する判断部と、該判断部により正しいと判断された該色変換結果に基づいて該色変換テーブルを作成する作成処理部とをそなえるとともに、該対応関係作成部が、該分類部により分類された分光反射率群毎に当該分光反射率群に対して特異値分解を施すことにより算出された複数の主成分ベクトルと、該分光特性推定部により推定された前記分光特性とに基づいて、該変換元色空間における色信号値を代入することにより該複数の主成分ベクトルのそれぞれについての複数の重み係数を算出する第１式を、該色変換式として該変換先色空間での前記領域毎に作成し、該変換元色空間における色信号値毎に該第１式によって算出された該複数の重み係数をそれぞれ該複数の主成分ベクトルに乗算しその線形和を求めることによって当該色信号値に対応した分光反射率を前記領域毎に算出する第２式を、該色変換式として作成し、該第２式によって前記領域毎に算出された該分光反射率を該変換先色空間における色信号に変換する第３式を、該色変換式として作成し、これらの第１式，第２式および第３式を用いて該変換元色空間の一の色信号を該変換先色空間の色信号に変換することにより、当該一の色信号について、該分類部により分類された複数の分光反射率群のそれぞれに対応する複数の該対応関係を作成することを特徴としている。
【００２２】
また、本発明のコンピュータ読取可能な記録媒体（請求項３）は、色入力機器から出力された、該色入力機器の色空間（以下、変換元色空間）における色信号と、該変換元色空間とは異なる色空間（以下、変換先色空間という）における色信号とを対応させる色変換テーブルを、色票上の複数種の色領域を該色入力機器により読み込んで得られた、各色領域に対応する、該変換元色空間での色信号であって入力部から入力される色信号と、測色器により該複数種の色領域を測色して得られる各色領域に対応する分光反射率とに基づいて作成する機能をコンピュータにより実現させるための色変換テーブル作成プログラムを記録したものであって、該色変換テーブル作成プログラムが、該分光反射率に対応する、該変換先色空間での色信号が該色変換先色空間において該複数の領域のうちのどの領域に属するかに応じて、該分光反射率を分類する分類部、該入力部から入力された該色信号と該測色器により得られた該分光反射率とに基づいて、該色入力機器の分光特性を推定する分光特性推定部、および、該変換先色空間を分割して得られる複数の領域にそれぞれ対応した複数の色変換式を用いて、該色変換テーブルを作成する色変換テーブル作成部として、該コンピュータを機能させ、該色変換テーブル作成プログラムが該コンピュータを該色変換テーブル作成部として機能させる際に、該色変換テーブル作成プログラムが、該分類部により分類された該分光反射率と該分光特性推定部により推定された前記分光特性とに基づいて、該変換先色空間での前記領域毎に該色変換式を作成しながら、各領域に応じた該色変換式を用いて、該色変換式毎に、該変換元色空間における色信号と該変換先色空間における色信号との対応関係を複数作成する対応関係作成部、該対応関係作成部において該対応関係を複数作成する際に各色変換式により該変換元色空間の一の色信号を変換した色変換結果が、各色変換式に対応した該領域に属している場合に、該色変換結果は正しいものと判断する判断部、および、該判断部により正しいと判断された該色変換結果に基づいて該色変換テーブルを作成する作成処理部として、該コンピュータをさらに機能させるとともに、該色変換テーブル作成プログラムが該コンピュータを該対応関係作成部として機能させる際に、色変換テーブル作成プログラムが、該コンピュータに、該分類部により分類された分光反射率群毎に当該分光反射率群に対して特異値分解を施すことにより算出された複数の主成分ベクトルと、該分光特性推定部により推定された前記分光特性とに基づいて、該変換元色空間における色信号値を代入することにより該複数の主成分ベクトルのそれぞれについての複数の重み係数を算出する第１式を、該色変換式として該変換先色空間での前記領域毎に作成させ、該変換元色空間における色信号値毎に該第１式によって算出された該複数の重み係数をそれぞれ該複数の主成分ベクトルに乗算しその線形和を求めることによって当該色信号値に対応した分光反射率を前記領域毎に算出する第２式を、該色変換式として作成させ、該第２式によって前記領域毎に算出された該分光反射率を該変換先色空間における色信号に変換する第３式を、該色変換式として作成させ、これらの第１式，第２式および第３式を用いて該変換元色空間の一の色信号を該変換先色空間の色信号に変換することにより、当該一の色信号について、該分類部により分類された複数の分光反射率群のそれぞれに対応する複数の該対応関係を作成させることを特徴としている。
【００２３】
上述の構成により、色入力機器の色変換テーブルの作成に際して、変換先色空間の分割領域に応じた色変換式（つまり分割領域毎に異なる色変換式）を用いることにより、色空間全体に亘る対応関係を精度良く算出することが可能になる。
【００２４】
このとき、分割領域毎に算出した、分光反射率に対する主成分ベクトルは、その各分割領域に存在する分光反射率に対する寄与率が高い。そこで、色票上の複数種の色領域を測色して得られた分光反射率を分割領域毎に分類し、分類された分光反射率と色入力機器についての推定分光特性とに基づき、各分割領域に応じた色変換式を用いて対応関係を算出することで、色空間全体から算出した、分光反射率に対する主成分ベクトルを用いるよりも精度良く対応関係を算出することが可能になる。
【００２５】
また、一の色信号について複数の色変換式により複数組の対応関係を作成し、その複数組の対応関係に基づいて色変換テーブルが作成されるので、その複数組の対応関係を考慮しながら、より精度の高い対応関係を算出することが可能になる。
【００２６】
【発明の実施の形態】
以下、図面を参照して本発明の実施の形態を説明する。
〔１〕本実施形態の構成の説明
図１は本発明の一実施形態としての色変換テーブル作成装置の機能構成を示すブロック図、図２は本発明の一実施形態としての色変換テーブル作成方法を実現するためのコンピュータシステム（色変換テーブル作成装置）のハードウェア構成を示すブロック図である。また、図３はＲＧＢ空間とL*a*b*空間との対応関係を示す図、図４（Ａ）〜図４（Ｃ）は、いずれも、本実施形態における色相角に応じた領域分割例を説明するための図（a*b*平面）である。
【００２７】
図１および図２に示す本実施形態のコンピュータシステム（パーソナルコンピュータ：ＰＣ）１は、色入力機器としてのイメージスキャナ６０の色空間（変換元色空間：本実施形態ではＲＧＢ空間）とデバイスに依存しない色空間（変換先色空間：本実施形態ではL*a*b*空間）との対応関係を示す色変換テーブルを作成する装置として機能するものである。この色変換テーブルにより、イメージスキャナ６０から出力されたＲＧＢ信号と、L*a*b*空間におけるとを対応させることができる。
【００２８】
通常、色入力機器用の色変換テーブルは、５０００（〜３００００）程度の色の対応関係を格納している。本実施形態では、例えば下記表１に示すごとく、ＲＧＢ各色について１７段階の値を有するＲＧＢ空間（１７×１７×１７個の格子点を有するＲＧＢ空間）における格子状の色変換テーブルを作成する場合について説明する。ここで、L*a*b*値の範囲は、例えば、−１００≦L*≦１００，−１００≦a*≦１００，−１００≦b*≦１００とする。
【００２９】
【表１】

【００３０】
ＲＧＢ空間における格子状の色変換テーブルとは、ＲＧＢ空間で格子状に分布する色（格子点の色）の信号値とその色のL*a*b*空間での色信号値との対応関係、つまり、色原稿をイメージスキャナで読み込んで出力されたＲＧＢ値と、その色原稿のL*a*b*測色値との対応関係を示すものである。この色変換テーブルには、ＲＧＢ値として、例えば表１のように規則的に増加するような「格子状のＲＧＢ値データセット」が格納されている。
【００３１】
ＲＧＢ値はＲＧＢ空間内において格子状に分布するが、ＲＧＢ空間とL*a*b*空間との間に線形な関係はないため、ＲＧＢ値をL*a*b*値に変換して得られるL*a*b*空間内での点は格子状に分布せず、その分布は歪んだものとなる。このようなＲＧＢ空間とL*a*b*空間との対応関係を二次元上で概念的に示した図が図３である。
【００３２】
なお、本実施形態では、変換先色空間が、人が見た際の色の差とその色空間での色の距離とがほぼ均等であるような、均等色空間（ＣＩＥＬＡＢ，ＣＩＥＬＵＶ，ＣＩＥＣＡＭなど）で、特にL*a*b*空間である場合について説明している。
L*a*b*空間は、ＣＩＥ（国際照明委員会）により定められた色空間であり、空間内の距離が、人間が感じる色の見えの差を表している。即ち、L*a*b*空間内のどのような色に対しても、そのL*a*b*空間内の距離を、色変換のための重み係数として用いることができる。これに対し、ＸＹＺ表色系（ｘｙ色度図）やマンセル表色系では、距離と色の見えの差とが対応していないため、色空間内の距離を色変換のための重み係数として用いることができない。ただし、L*a*b*空間以外にもL*u*v*空間やＣＩＥＣＡＭ空間など、L*a*b*空間と同様に距離と色の見えの差が等しい座標系があるため、これらを用いてもよい。
【００３３】
さて、図２に示すように、コンピュータシステム１においては、ＣＰＵ１０とこのＣＰＵ１０に接続されるバスライン５０とがそなえられている。バスライン５０には、メモリ（ＲＯＭ，ＲＡＭ）２０が接続されるとともに、入出力インタフェース５１を介してディスプレイ（例えばＣＲＴ，ＬＣＤ，ＰＤＰ等）３０，キーボード４０およびマウス４１が接続されている。
【００３４】
ここで、ディスプレイ３０は、ＣＰＵ１０によってその表示状態を制御され、色変換テーブルの作成にかかる各種情報を表示するものである。また、キーボード４０およびマウス４１は、ディスプレイ３０における画面を参照したオペレータ（ユーザ）によって操作され、ＣＰＵ１０（コンピュータシステム１）に対して各種指示や各種情報を入力するものである。
【００３５】
メモリ２０には、図１に示す分光特性推定部１１，分光反射率分類部１２，色変換テーブル作成部１４およびＲＧＢ空間格子点信号値出力部１５からなる色変換テーブル作成装置を実現するためのアプリケーションプログラム（色信号値抽出プログラム）２１が格納されている。
【００３６】
そして、ＣＰＵ１０がバスライン５０を介してメモリ２０から上記アプリケーションプログラム２１を読み出して実行することにより、分光特性推定部１１，分光反射率分類部１２，色変換テーブル作成部１４およびＲＧＢ空間格子点信号値出力部１５としての機能（その詳細については後述）、つまりは、本発明の色変換テーブル作成装置としての機能が実現され、本発明の色変換テーブル作成方法が実行されるようになっている。
【００３７】
上述したアプリケーションプログラム２１は、例えばフレキシブルディスク，ＣＤ−ＲＯＭ等の、コンピュータ読取可能な記録媒体に記録された形態で提供される。そして、コンピュータシステム１はその記録媒体からプログラム２１を読み取って内部記憶装置または外部記憶装置に転送し格納して用いる。また、そのプログラム２１を、例えば磁気ディスク，光ディスク，光磁気ディスク等の記憶装置（記録媒体）に記録しておき、その記憶装置から通信経路を介してコンピュータシステム１に提供するようにしてもよい。
【００３８】
分光特性推定部１１，分光反射率分類部１２，色変換テーブル作成部１４およびＲＧＢ空間格子点信号値出力部１５としての機能を実現する際には、内部記憶装置（本実施形態ではメモリ２０）に格納されたプログラム２１がコンピュータのマイクロプロセッサ（本実施形態ではＣＰＵ１０）によって実行される。このとき、記録媒体に記録されたプログラム２１をコンピュータシステム１が直接読み取って実行するようにしてもよい。
【００３９】
なお、本実施形態において、コンピュータとは、ハードウェアとオペレーションシステムとを含む概念であり、オペレーションシステムの制御の下で動作するハードウェアを意味している。また、オペレーションシステムが不要でアプリケーションプログラム単独でハードウェアを動作させるような場合には、そのハードウェア自体がコンピュータに相当する。ハードウェアは、少なくとも、ＣＰＵ等のマイクロプロセッサと、記録媒体に記録されたコンピュータプログラムを読み取るための手段とをそなえている。
【００４０】
アプリケーションプログラム２１は、このようなコンピュータ（コンピュータシステム１）に分光特性推定部１１，分光反射率分類部１２，色変換テーブル作成部１４およびＲＧＢ空間格子点信号値出力部１５としての機能を実現させるためのプログラムコードを含んでいる。また、その機能の一部は、アプリケーションプログラム２１ではなくオペレーションシステムによって実現されてもよい。
【００４１】
さらに、本実施形態における記録媒体としては、上述したフレキシブルディスク，ＣＤ−ＲＯＭ，ＤＶＤ，磁気ディスク，光ディスク，光磁気ディスクのほか、ＩＣカード，ＲＯＭカートリッジ，磁気テープ，パンチカード，コンピュータの内部記憶装置（ＲＡＭやＲＯＭなどのメモリ），外部記憶装置等や、バーコードなどの符号が印刷された印刷物等の、コンピュータ読取可能な種々の媒体を利用することができる。
【００４２】
一方、図１および図２に示すコンピュータシステム１においては、イメージスキャナ（色入力機器）６０および測色器７０が、入出力インタフェース５２を介してバス５０に接続されている。
本実施形態のコンピュータシステム１は、イメージスキャナ６０による色票１００（図１参照）の読取結果と、測色器７０による色票１００の測色結果とに基づいて、イメージスキャナ６０用の色変換テーブルを作成する。
【００４３】
このとき、入出力インタフェース５２は、色票１００上の数十種の色領域をイメージスキャナ６０により読み込んで得られた、各色領域に対応するＲＧＢ信号や、各色領域についての測色器７０による測色結果（後述する分光反射率およびL*a*b*値）を、コンピュータシステム１に入力する入力部として機能することになる。この入出力インタフェース５２から入力されたＲＧＢ信号や測色結果は、メモリ２０に一時的に格納されるようになっている。
【００４４】
測色器７０は、前述した通り、色入力機器用色変換テーブル作成時に色票の各色領域を測色するために用いられ、３８０〜７３０nmの分光反射率を１０nm分解能で出力するものである。このとき出力される分光反射率は、最大値を１．０として正規化されている。
【００４５】
また、本実施形態の測色器７０は、分光反射率を測定・出力する機能のほか、その分光反射率をL*a*b*値に変換して出力する機能も有しているものとする。L*a*b*値は、後述するごとくＣＰＵ１０（分光反射率分類部１２）により分光反射率の分類を行なう際に必要になるものである。
なお、本実施形態では、分光反射率に対応するL*a*b*値を、測色器７０で得ているが、ＣＰＵ１０（分光反射率分類部１２）において測色器７０からの分光反射率を変換して算出するように構成してもよい。
【００４６】
色票１００は、少なくとも数十種程度のカラーパッチ（色領域）を有するものである。なお、色票１００としては、ＡＮＳＩ(American National Standards Institute：アメリカ規格協会)規格の入力カラーターゲット(ANSI/IT8.7/2)などを用いてもよい。
【００４７】
この場合、その入力カラーターゲットは、２６４色の有彩色パッチからなるカラーパッチ領域と、２４段階の異なる階調レベルの無彩色パッチからなるグレイパッチ領域（グレイスケール）とを有しているので、本実施形態では、これらの色領域の中から適当な数十種の色領域を指定し、指定された色領域について、イメージスキャナ６０による読込や測色器７０による測色を行なうことになる。
【００４８】
さらに、メモリ２０には、前述したアプリケーションプログラム（色変換テーブル作成プログラム）２１や、ＣＰＵ１０により作成された色変換テーブル２２が格納される。
また、図１を参照しながら後述するごとく、メモリ２０は、色変換テーブル２２を作成する際のワーキングメモリとしても機能するもので、このメモリ２０には、イメージスキャナ６０からのＲＧＢ信号，測色器７０からの分光反射率（L*a*b*値），分光特性，分類情報，対応関係１〜３等が格納される。
【００４９】
さて、次に、ＣＰＵ１０により実現される各種機能（分光特性推定部１１，分光反射率分類部１２，色変換テーブル作成部１４およびＲＧＢ空間格子点信号値出力部１５としての機能）について、詳細に説明する。
分光特性推定部１１は、同一の色票１００（同一のカラーパッチ）についてイメージスキャナ６０および測色器７０によりそれぞれ得られたＲＧＢ信号および分光反射率をメモリ２０から読み出し、これらのＲＧＢ信号および分光反射率に基づいて、イメージスキャナ６０の分光特性を推定するものである。この分光特性の推定は(1)〜(10)式を参照しながら後述する手法によって行なわれる。
【００５０】
分光反射率分類部１２は、測色器７０によって得られた各分光反射率に対応するL*a*b*値（L*a*b*空間での色信号）がL*a*b*空間におけるどの領域に属するかを判断し、測色器７０によって得られた、数十種のカラーパッチについての分光反射率を、属する領域に応じてグループに分類するものである。その分類結果は、分類情報としてメモリ２０に出力されて格納されるようになっている。
【００５１】
このとき、本実施形態では、L*a*b*空間の全体を分割して得られる３つの領域１，領域２および領域３が、予め定義され設定されている（分割ステップ）。これらの３つの領域（以下、分割領域という場合がある）１〜３は、例えば図４（Ａ）〜図４（Ｃ）にそれぞれ示すごとく、L*a*b*空間を、色相角に応じた境界（実際にはＬ＊軸に平行な境界面）Ｂ１，Ｂ２およびＢ３により区分けして得られたものである。
【００５２】
図４（Ａ）に示すように、領域１は、一次色Ｃ（シアン）の色相角２１０度を中心とした１２０度の範囲の領域であり、図４（Ｂ）に示すように、領域２は、一次色Ｍ（マゼンタ）の色相角３３０度を中心とした１２０度の範囲の領域であり、図４（Ｃ）に示すように、領域３は、一次色Ｙ（イエロー）の色相角９０度を中心とした１２０度の範囲の領域である。
【００５３】
また、領域１〜３には、これらの領域１〜３を相互に重複させるべく、それぞれ図４（Ａ）〜図４（Ｃ）に示すような幅ｄ（一定）の領域が付加・形成されている。このような領域を形成することにより、図９〜図１２に示すような重複領域、即ち、領域１〜３が重複する領域と、領域１と領域２とが重複する幅２ｄの領域と、領域２と領域３とが重複する幅２ｄの領域と、領域３と領域１とが重複する幅２ｄの領域とが形成されるようになっている。
【００５４】
色変換テーブル作成部１４は、前述した３つの領域１〜３にそれぞれ対応した色変換式を用いて、色変換テーブル２２を作成するもので、対応関係作成部１４１，判断部１４２および作成処理部１４３を有して構成されている。また、ＲＧＢ空間格子点信号値出力部１５は、上記表１を参照しながら説明したような、ＲＧＢ空間で格子状に分布する色の信号値（格子点信号値）を生成して色変換テーブル作成部１４（対応関係作成１４１）へ出力するものである。
【００５５】
ここで、対応関係作成部１４１は、分光反射率分類部１２により分類された各分光反射率群や、分光特性推定部１２により推定されたイメージスキャナ６０の分光特性に基づいて、L*a*b*空間での領域毎に色変換式を作成し、３種類の色変換式のそれぞれを用いて、ＲＧＢ空間格子点信号値出力部１５からのＲＧＢ値（変換元色空間における色信号）と、L*a*b*値（変換先色空間における色信号）との対応関係を作成するものである。
【００５６】
なお、色変換式は、対応関係作成部１４１において、対応関係を作成するのに伴い、(11)〜(16)式を参照しながら後述する手法で作成されるものであり、ＲＧＢ値をL*a*b*値に変換する色変換式は、実際には、(14)〜(16)式と、ＸＹＺ値をL*a*b*値に変換する式（周知の一般的な式であるため表記せず）とを組み合わせたものである。
【００５７】
このとき、対応関係作成部１４１は、ＲＧＢ空間内の全ての格子点のそれぞれについて、３種類の色変換式を用いて３つのL*a*b*値を算出することにより３種類の対応関係を作成し、その３種類の対応関係をそれぞれ対応関係１〜３としてメモリ２０に出力して格納するようになっている。
【００５８】
判断部１４２は、対応関係作成部１４１が対応関係１〜３を作成する際に各色変換式によりＲＧＢ信号（一格子点信号値）を変換した色変換結果（L*a*b*値）が、各色変換式に対応した領域に属しているか否かを判断し、その領域に属している場合、その色変換結果は正しいものと判断するものである。
【００５９】
より具体的に説明すると、判断部１４２は、領域１に対応した色変換式により作成された対応関係１〔L*a*b*1（L1,a1,b1）〕が、図４（Ａ）に示す領域１に属している否かを判断し、領域１に属している場合、その対応関係１は正しいものと判断し、同様に、領域２に対応した色変換式により作成された対応関係２〔L*a*b*2（L2,a2,b2）〕が、図４（Ｂ）に示す領域２に属している否かを判断し、領域２に属している場合、その対応関係２は正しいものと判断し、領域３に対応した色変換式により作成された対応関係３〔L*a*b*3（L3,a3,b3）〕が、図４（Ｃ）に示す領域３に属している否かを判断し、領域３に属している場合、その対応関係３は正しいものと判断する。
【００６０】
作成処理部１４３は、対応関係作成部１４１で３種類の色変換式によりそれぞれ作成された３組の対応関係１〜３に基づいて、一組の対応関係を得て、色変換テーブル２２をメモリ２０上に作成するものである。このとき、作成処理部１４３は、判断部１４２で正しいと判断された色変換結果（L*a*b*値）に基づいて色変換テーブル２２を作成するようになっている。
【００６１】
より具体的に説明すると、ＲＧＢ空間の一格子点（一の色信号）について判断部１４２で正しいと判断された色変換結果が複数存在する場合（本実施形態では２または３の色変換結果が正しいと判断された場合）、作成処理部１４３は、図９，図１０および(17)式,(18)式を参照しながら後述するごとく、正しいと判断された複数の色変換結果（L*a*b*値）と各色変換結果の属する領域の境界との距離に基づいて、その一格子点についての色変換結果（L*a*b*値）を算出して出力するようになっている。
【００６２】
また、ＲＧＢ空間の一格子点について判断部１４２で正しいと判断された色変換結果が一つだけである場合、作成処理部１４３は、図１１および(19)式を参照しながら後述するごとく、正しいと判断された一の色変換結果（L*a*b*値）を、そのまま、その一格子点についての色変換結果（L*a*b*値）として出力するようになっている。
【００６３】
ＲＧＢ空間の一格子点について判断部１４２で正しいと判断された色変換結果が存在しない場合、作成処理部１４３は、図１２および(20)式を参照しながら後述するごとく、対応関係作成部１４１でその一格子点について得られた３つの色変換結果と各色変換結果の属する領域の境界との距離の逆数に基づいて、その一格子点についての色変換結果を算出して出力するようになっている。
【００６４】
なお、ＲＧＢ空間の一格子点について判断部１４２で正しいと判断された色変換結果が複数存在する場合、作成処理部１４３を、最も信頼性の高い色変換結果を選択して出力するように機能させてもよい。
この場合、作成処理部１４３は、正しいと判断された複数の色変換結果と各色変換結果の属する領域の境界との距離に基づいて、正しいと判断された複数の色変換結果のうちの一つを、その一格子点についての色変換結果として選択して出力するように構成される。このとき、作成処理部１４３は、例えば前記距離が最も大きい色変換結果（境界から最も離れている色変換結果）を選択する。
【００６５】
同様に、ＲＧＢ空間の一格子点について判断部１４２で正しいと判断された色変換結果が存在しない場合、作成処理部１４３を、最も信頼性の高い色変換結果を選択して出力するように機能させてもよい。
この場合、作成処理部１４３が、対応関係作成部１４１でその一格子点について得られた３つの色変換結果と各色変換結果の属する領域の境界との距離の逆数に基づいて、３つの色変換結果のうちの一つを、その一格子点についての色変換結果として選択して出力するように構成される。このとき、作成処理部１４３は、例えば前記距離の逆数が最も大きい色変換結果（境界に最も近い色変換結果）を選択する。
【００６６】
また、図１に示す例では、対応関係作成部１４１により作成された対応関係１〜３は、一旦、メモリ２０に格納されているが、その対応関係１〜３は、対応関係作成部１４１から、直接、判断部１４２および作成処理部１４３に入力するように構成してもよい。
【００６７】
〔２〕本実施形態の動作の説明
〔２−１〕色変換テーブル作成手順の説明
次に、本発明の一実施形態としての色変換テーブル作成方法の流れ（上述のごとく構成された色変換テーブル作成装置１の動作）について、図５に示すフローチャート（ステップＳ１１〜Ｓ２５）に従って説明する。
【００６８】
本実施形態では、色変換テーブル２２の作成を開始するに先立って、L*a*b*空間の全体を分割して得られる３つの領域１，領域２および領域３が、予め定義され設定されている（分割ステップ；図示せず）。
そして、本実施形態の色変換テーブル作成装置１によってイメージスキャナ６０の色変換テーブル２２を作成する際には、図５に示すように、まず、色票１００上の数十種のカラーパッチがイメージスキャナ６０で読み込まれ（ステップＳ１１；読込ステップ）、各カラーパッチに対応するＲＧＢ信号が、イメージスキャナ６０から色変換テーブル作成装置１に出力され、入出力インタフェース５２を通じてメモリ２０に格納される（ステップＳ１２；読込ステップ）。
【００６９】
また、イメージスキャナ６０で読み込んだものと同一の数十種のカラーパッチが測色器７０で測色され（ステップＳ１３；測色ステップ）、各カラーパッチに対応する分光反射率が、測色器７０から出力され、入出力インタフェース５２を通じてメモリ２０に格納される（ステップＳ１４；測色ステップ）。このとき、本実施形態の測色器７０では、測定された分光反射率が、各カラーパッチに対応するL*a*b*値に変換され、そのL*a*b*値が、入出力インタフェース５２を通じてメモリ２０に格納される（ステップＳ１５）。
【００７０】
そして、分光特性推定部１１において、メモリ２０に格納された各カラーパッチについてのＲＧＢ信号および分光反射率に基づき、(1)〜(10)式を参照しながら後述する手法に従って、イメージスキャナ６０の分光特性が推定され（ステップＳ１６；分光特性推定ステップ）、その分光特性がメモリ２０に出力されて格納される（ステップＳ１７）。
【００７１】
また、分光反射率分類部１２において、各カラーパッチについてのL*a*b*値が図４（Ａ）〜図４（Ｃ）にそれぞれ示す領域１〜３のいずれに属するかが判断され、各カラーパッチについての分光反射率が、その領域に応じてグループ１〜Ｋ（本実施形態ではＫ＝３）に分類され（ステップＳ１８；分類ステップ）、その分類結果（グループ１，２，３）は分類情報としてメモリ２０に格納される（ステップＳ１９）。
【００７２】
この後、色変換テーブル作成部１４において、L*a*b*空間での領域１〜３にそれぞれ対応した３種類の色変換式を、分光反射率分類部１１により分類された各分光反射率群（各グループ１〜３）や、分光特性推定部１２により推定されたイメージスキャナ６０の分光特性に基づいて作成しながら、色変換テーブル２２が作成される（ステップＳ２１〜Ｓ２５；色変換テーブル作成ステップ）。
【００７３】
このとき、ＲＧＢ空間格子点信号値出力部１５において、ＲＧＢ空間で格子状に分布する色のＲＧＢ値（格子点信号値）が生成され、そのＲＧＢ値が色変換テーブル作成部１４（対応関係作成部１４１）へ出力される（ステップＳ２０）。
そして、対応関係作成部１４１において、領域１〜３のそれぞれに対応した３種類の色変換式を用い、その色変換式毎に、ステップＳ２１からのＲＧＢ値（格子点信号値）と、L*a*b*値との対応関係が作成される（ステップＳ２１；対応関係作成ステップ）。
【００７４】
このとき、ＲＧＢ値をL*a*b*値に変換する色変換式としては、後述する(14)〜(16)式と、ＸＹＺ値をL*a*b*値に変換する式（周知の一般的な式）とが用いられる。また、ステップＳ２１では、ＲＧＢ空間内の全ての格子点のそれぞれについて、３種類の色変換式を用いて３つのL*a*b*値が算出され、３種類の対応関係１〜３が作成されて出力される（ステップＳ２２）。
【００７５】
ついで、判断部１４２において、前述したように、ステップＳ２１で作成された各対応関係１〜３に含まれるL*a*b*値（色変換結果）が、各色変換式に対応した領域に属しているか否かを判断し、その領域に属している場合、その色変換結果は正しいものと判断して、その判断結果を作成処理部１４３に出力する（ステップＳ２３；判断ステップ）。
【００７６】
そして、作成処理部１４３において、ステップＳ２１で色変換式毎に作成された対応関係１〜３に基づいて、一組の対応関係を決定し、色変換テーブル２２が作成される（ステップＳ２４，Ｓ２５）。その際、基本的には、ＲＧＢ空間の各格子点についての対応関係１〜３に含まれるL*a*b*値のうち、ステップＳ２４で正しいと判断されたL*a*b*値に基づいて、その格子点についての最終的なL*a*b*値が決定され（ステップＳ２４；作成処理ステップ）、このように決定されたL*a*b*値が、その格子点と対応関係を成すものとしてメモリ２０に出力され、メモリ２０上に色変換テーブル２２が作成される（ステップＳ２５）。
【００７７】
２または３の色変換結果が正しいと判断された場合、前述した通り、作成処理部１４３においては、正しいと判断された２または３の色変換結果（L*a*b*値）と各色変換結果の属する領域の境界との距離に基づいて、その格子点についての色変換結果（L*a*b*値）が算出もしくは選択されて決定される。
正しいと判断された色変換結果が一つだけである場合は、前述した通り、作成処理部１４３においては、正しいと判断された一の色変換結果（L*a*b*値）が、そのまま、その格子点についての色変換結果（L*a*b*値）として出力される。
【００７８】
これに対し、正しいと判断された色変換結果が存在しない場合には、前述した通り、作成処理部１４３においては、対応関係作成部１４１でその格子点について得られた３つの色変換結果（L*a*b*値）と各色変換結果の属する領域の境界との距離の逆数に基づいて、その格子点についての色変換結果が算出もしくは選択されて決定される。
【００７９】
〔２−２〕分光特性推定手順の説明
次に、図５のステップＳ１６における処理、つまり、本実施形態の分光特性推定手順（分光特性推定部１１の動作）について、図６に示すフローチャート（ステップＳ３１〜Ｓ３５）に従って詳細に説明する。
分光特性推定部１１（ステップＳ１６）では、数十種の色信号および数十種の分光反射率を用いて、以下のようにして分光特性を推定する。
【００８０】
〔２−２−１〕連立一次方程式の作成（ステップＳ３１）
ここでは、分光特性を、可視光波長範囲380〜730nmにおける１０nm毎の各波長（380,390,…,730nm）に対応して得られる離散データとして扱う。
一般に、分光特性にカラーパッチの分光反射率を乗算して得られる値が、ＲＧＢ値になる。
【００８１】
そこで、イメージスキャナ６０の分光特性を推定すべく、ステップＳ３１では、イメージスキャナ６０によりカラーパッチを読み込んで得られたＲＧＢ値と、測色器７０により得られたカラーパッチの分光反射率とを用いて、Ｒ，Ｇ，Ｂのそれぞれについて、分光特性を未知数とする３６以上の方程式からなる連立一次方程式を下式(1)のごとく作成する。
【００８２】
【数１】

【００８３】
上式(1)において、Ｒ_i，Ｇ_i，Ｂ_i(i=1,...,N)は、それぞれＮ個のカラーパッチをイメージスキャナ６０により読み取って得られた値（ＲＧＢ値）であり、REF_(i,j)(i=1,...,N；j=380,390,...,730)は、Ｎ個のカラーパッチについて測色器７０によりそれぞれ得られた、Ｎ組の分光反射率データであり、（REF_(i,380) REF_(i,390) … REF_(i,730)）は、ｉ番目のカラーパッチについて得られた分光反射率（３６個の離散データ）である。また、LS_Rj(j=380,390,...,730)は、Ｒに対する分光特性（３６個の離散データ）、LS_Gj(j=380,390,...,730)は、Ｇに対する分光特性（３６個の離散データ）、LS_Bj(j=380,390,...,730)は、Ｂに対する分光特性（３６個の離散データ）である。
【００８４】
カラーパッチの数Ｎが３６未満である場合は、３６個の連立一次方程式を作成することができない。この場合、ステップＳ３１では、下式(2)のような、３６個の一次方程式からなる連立一次方程式を作成する。
【００８５】
【数２】

【００８６】
〔２−２−２〕連立一次方程式の特異値分解（ステップＳ３２,Ｓ３３）
ここで、上式(1)もしくは(2)中の以下の部分〔つまり下式(3)〕を、分光反射率データセットと呼ぶ。
【００８７】
【数３】

【００８８】
ステップＳ３２では、この分光反射率データセットに対して特異値分解を施した後、ステップＳ３３で、特異値分解の結果に基づいて、分光反射率データに含まれるノイズを除去する。
ノイズ除去は以下の手順で行なわれる。まず、分光反射率データセットに対して特異値分解を行ない、上記(3)式に示す分光反射率データセットは、下式(4)に示す３つの行列に分解される。
【００８９】
【数４】

【００９０】
上式(4)において、REF_(i,j)(i=1,...,N；j=380,390,...,730)は、前述したように、測色器７０により得られたＮ組のカラーパッチの分光反射率データであり、Ｗ_m(m=1,...,36)は、主成分ベクトルに乗算される重みである。また、Ｐ_(m,j)(m=1,...,36;j=380,390,...,730)は、分光反射率に対する３６組の主成分ベクトルで、各主成分ベクトルは３６個の離散データとなっている。さらに、Ｕ_im(i=1,...,N；m=1,...,36)は係数である。
【００９１】
特異値分解法によって算出された重みＷ₁〜Ｗ₃₆は、それぞれ、主成分ベクトルＰ_(1,j)〜Ｐ_(36,j)(j=380,390,...,730)に対する重み（重要度）を示すものであって、各重みＷ₁〜Ｗ₃₆の大きさに比例して、それに対応する主成分ベクトルの重要度は高くなる。
【００９２】
そこで、分光反射率データセットの有効桁数から考えると、最も大きな重みに対して有効桁数以上離れている重みは、意味の無い主成分ベクトルであるとみなすることができる。ステップＳ３３では、そのような主成分ベクトルを除去することにより、分光反射率データに含まれるノイズを除去する。
つまり、ステップＳ３３では、まず、下式(5)によって閾値（通常、条件数と呼ぶ）を算出し、重みＷ₁〜Ｗ₃₆のうち、その閾値以下のものを、下式(6)のように“０”に置き換え、ノイズ除去を行なっている。
【００９３】
閾値（条件数）
＝〔重みＷ₁〜Ｗ₃₆の中の最大値〕/〔データセット（測色器）の有効桁数〕(5)
【００９４】
なお、下式(6)は、上記閾値を超えた重みがＭ個あった場合の式であり、これらＭ個の重みをＷ₁〜Ｗ_Mとし、Ｍ＋１番目〜３６番目の主成分ベクトルに対する重みＷ_M+1〜Ｗ₃₆が“０”に置き換えられている。当然、Ｍ＜３６である。
【００９５】
【数５】

【００９６】
〔２−２−３〕分光特性の算出（ステップＳ３４，Ｓ３５）
ここでは、簡単のため、上式(6)に示す行列を、下式(7)に示すような行列に置き換える。
【００９７】
【数６】

【００９８】
ここで、Ｕは、前記係数Ｕ_im(i=1,...,N；m=1,...,36)を要素とするＮ×３６の行列であり、Ｖは、主成分ベクトルＰ_(m,j)(m=1,...,36;j=380,390,...,730)を要素とする３６×３６の行列である。上式(7)中における行列Ｕ，Ｖには、下式(8-1),(8-2)で示すような性質がある。なお、下式(8-1),(8-2)において、Ｉは単位行列、Ｕ^Tは行列Ｕの転置行列、Ｖ^Tは行列Ｖの転置行列である。
【００９９】
【数７】

【０１００】
従って、(8-1)式および(8-2)式を用いることにより、上式(1)もしくは(2)は下式(9)のように展開されるので、その展開式(9)から得られた下式(10)を用いて、イメージスキャナ６０の分光特性が算出され（ステップＳ３４）、算出された分光特性は、メモリ２０に出力され格納される（ステップＳ３５）。
【０１０１】
【数８】

【０１０２】
【数９】

【０１０３】
〔２−３〕分光反射率分類手順の説明
次に、図５のステップＳ１８における処理、つまり、本実施形態の分光反射率分類手順（分光反射率分類部１２の動作）について、図４（Ａ）〜図４（Ｃ）を参照しながら説明する。
【０１０４】
本実施形態では、前述した通り、色票１００の一次色Ｃ，Ｍ，Ｙの色相角をそれぞれ２１０度，３３０度，９０度とし、予め、図４（Ａ）〜図４（Ｃ）に示すように、ＣＩＥＬＡＢ空間（L*a*b*空間）が、それぞれ、前記色相角２１０度，３３０度，９０度を中心とした３つの領域１〜３に分割されている。また、前述した通り、重複領域を形成すべく、各領域１〜３には幅ｄの領域が付加・形成されているが、この幅ｄは、例えば５．０に設定される。
【０１０５】
そして、分光反射率分類部１２（ステップＳ１８）では、各カラーパッチについてのL*a*b*値が図４（Ａ）〜図４（Ｃ）にそれぞれ示す領域１〜３のいずれに属するかが判断され、各カラーパッチについての分光反射率が、その領域に応じて３つのグループ１〜３に分類される。
つまり、そのL*a*b*値が領域１に属すると判断された分光反射率は、グループ１に分類され、そのL*a*b*値が領域２に属すると判断された分光反射率は、グループ２に分類され、そのL*a*b*値が領域３に属すると判断された分光反射率は、グループ３に分類される。
【０１０６】
なお、本実施形態では、領域分割時の判断に用いられるL*a*b*値を、測色器７０から得ているが、分光反射率分類部１２において、測色器７０からの分光反射率を変換することにより、各カラーパッチの分光反射率に対応するL*a*b*値を得てもよい。
【０１０７】
〔２−４〕色変換式作成手順および対応関係作成手順の説明
次に、図５のステップＳ２０〜Ｓ２２における処理、つまり、本実施形態の色変換式作成手順および対応関係作成手順（対応関係作成部１４１の動作）について、図７に示すフローチャート（ステップＳ４１〜Ｓ４７）に従って詳細に説明する。
【０１０８】
対応関係作成部１４１では、格子状のＲＧＢ値に対応するL*a*b*値が算出される。
つまり、格子状のＲＧＢ値と、イメージスキャナ６０の分光特性と、グループ１の分光反射率から算出された主成分ベクトルとに基づいて、色空間全体に亘る対応関係１が算出され、同様に、格子状のＲＧＢ値と、イメージスキャナ６０の分光特性と、グループ２の分光反射率から算出された主成分ベクトルとに基づいて、色空間全体に亘る対応関係２が算出されるとともに、格子状のＲＧＢ値と、イメージスキャナ６０の分光特性と、グループ３の分光反射率から算出された主成分ベクトルとに基づいて、色空間全体に亘る対応関係３が算出される。これにより、１つの格子状のＲＧＢ値に対して、３種類の色変換式を用い、色変換結果としての３つのL*a*b*値（L*a*b*1，L*a*b*2，L*a*b*3）が算出される。
【０１０９】
ここで、格子状のＲＧＢ値は、ＲＧＢ空間格子点信号値出力部１５から出力される。また、イメージスキャナ６０の分光特性は、分光特性推定部１１において上式(10)により推定・算出された、LS_Rj(j=380,390,...,730)， LS_Gj(j=380,390,...,730)およびLS_Bj(j=380,390,...,730)である。さらに、各グループ１〜３の主成分ベクトルは、各グループ１〜３の分光反射率に対して特異値分解を施すことにより、上述したような(4)〜(6)式を用いて、Ｐ_(m,j)(m=1,...,Ｍ;j=380,390,...,730)として算出される。
【０１１０】
以下に、図７に示すフローチャート（ステップＳ４１〜Ｓ４７）を参照しながら、ある一つのグループ（分光反射率群）に対する対応関係の算出手順について説明する。
【０１１１】
〔２−４−１〕主成分ベクトルに乗算する重みの算出（ステップＳ４１，Ｓ４２）
分光反射率は、分光反射率群（本実施形態ではグループ１，２，３）に対して特異値分解法を施すことににより算出された、Ｍ組の主成分ベクトルの線形和として、下式(11)のごとく再現される。なお、この(11)式は、上式(4),(6)におけるｉ行目の要素を抜き出したものである。
【０１１２】
【数１０】

【０１１３】
ここで、REF_(i,j)(j=380,390,...,730)は、各グループ１〜３に分類された、ｉ番目の分光反射率である。なお、ｉは１,...,Ｎであり、Ｎは、各グループ１〜３に分類された分光反射率REF_(i,j)の数となる。また、Ｐ_(m,j)(m=1,...,Ｍ;j=380,390,...,730)はＭ組の主成分ベクトルである。さらに、Ｕ_im・Ｗ_m(m=1,...,Ｍ)は、ｉ番目の分光反射率に対応する重み係数で、それぞれ、ｍ番目の主成分ベクトル（Ｐ_(m,380) Ｐ_(m,390) … Ｐ_(m,730)）に乗算されるものである。
【０１１４】
上式(11)に示すごとく、主成分ベクトルを用いることにより、３６次元データである分光反射率を、数次元程度（３〜７次元；Ｍ＝３〜７）の重み係数（Ｕ・Ｗ）で表現することが可能になる。
【０１１５】
あるカラーパッチをイメージスキャナ６０で読み込み、このイメージスキャナ６０から出力されたＲＧＢ値(R G B)は、推定されたイメージスキャナ６０の分光特性LS_Rj(j=380,390,...,730)，LS_Gj(j=380,390,...,730)，LS_Bj(j=380,390,...,730)と、そのカラーパッチの色原稿を測色して得られた分光反射率REF_j(j=380,390,...,730)とに基づいて、下式(12)のごとく算出することができる。
【０１１６】
【数１１】

【０１１７】
ここで、分光反射率REF_j(j=380,390,...,730)は、上式(11)に示した通り、Ｍ個の主成分ベクトルの線形結合として表されるので、上式(12)は、下式(13)に書き換えられる。
【０１１８】
【数１２】

【０１１９】
上式(13)における行列部分Ａの逆行列を算出することができれば、下式(14)のごとく、ＲＧＢ値から重み係数Ｕ_m・Ｗ_m(m=1,...,Ｍ)を算出することができる。
【０１２０】
【数１３】

【０１２１】
以下、本実施形態では、Ｍ＝３の場合（即ち主成分ベクトルの数が３の場合）について説明する。この場合、上記行列部分Ａが３×３行列となり、ＲＧＢ値(R G B)を上式(14)に代入することにより、そのＲＧＢ値に対応した重み係数Ｕ₁・Ｗ₁，Ｕ₂・Ｗ₂，Ｕ₃・Ｗ₃が一意に算出されることになる。
本実施形態では、グループ１〜３のそれぞれについての３組の主成分ベクトルＰ_(m,j)(m=1,2,3;j=380,390,...,730)が、各グループ１〜３の分光反射率に対して特異値分解を施すことにより、上述したような(4)〜(6)式を用いて算出される（ステップＳ４１）。
【０１２２】
そして、算出された主成分ベクトルＰ_(m,j)(m=1,2,3;j=380,390,...,730)と、分光特性LS_Rj(j=380,390,...,730)，LS_Gj(j=380,390,...,730)，LS_Bj(j=380,390,...,730)とにより、３種類の(14)式が作成される。この(14)式に、ＲＧＢ空間格子点信号値出力部１５からのＲＧＢ値を代入することにより、３組の主成分ベクトルに乗算される重み係数Ｕ₁・Ｗ₁，Ｕ₂・Ｗ₂，Ｕ₃・Ｗ₃が算出される（ステップＳ４２）。
【０１２３】
〔２−４−２〕分光反射率の算出
ついで、上式(14)を用いてＲＧＢ値毎に得られた重み係数Ｕ₁・Ｗ₁，Ｕ₂・Ｗ₂，Ｕ₃・Ｗ₃を、それぞれ３組の主成分ベクトルＰ_(1,j)(j=380,390,...,730)，Ｐ_(2,j)(j=380,390,...,730)，Ｐ_(3,j)(j=380,390,...,730)に乗算し、その線形和を求めることによって、下式(15)に示すごとく、ＲＧＢ値に対応した分光反射率REF_j(j=380,390,...,730)が算出される（ステップＳ４３）。
【０１２４】
【数１４】

【０１２５】
このとき、上式(15)によって算出された分光反射率REF_j(j=380,390,...,730)の値が１．０未満である場合（つまり、万一、負値となった場合）は、その値を０．０に書き換える一方、万一、１．０を超えた場合は、１．０に書き換えるように、分光反射率の補正を行なってから（ステップＳ４４）、補正後の分光反射率を用いて、下式(16)によりＸＹＺ値を算出する（ステップＳ４５）。
【０１２６】
【数１５】

【０１２７】
このようにしてＲＧＢ値からＸＹＺ値を算出した後、さらに、周知の変換式を用いて、上式(16)により算出されたＸＹＺ値をL*a*b*値に変換する（ステップＳ４６）。これにより、格子状に並ぶ、ＲＧＢ値とL*a*b*値との対応関係が作成される。
上述のような対応関係作成は、分光反射率群であるグループ１〜３のそれぞれについて行なわれ、その結果として、例えば下記表２に示すごとく、１つのＲＧＢ信号値（格子点信号値）について３種類の対応関係１〜３が作成されてメモリ２０に出力・格納される（ステップＳ４７）。
【０１２８】
なお、本実施形態では、前述した通り、ＲＧＢ値をL*a*b*値に変換する色変換式は、上式(14)〜(16)と、ＸＹＺ値をL*a*b*値に変換する式とを組み合わせたものであり、色変換テーブル作成部１４（対応関係作成部１４１）では、これらの色変換式が、分類された分光反射率群毎に即ち分割領域毎に作成されながら、ＲＧＢ値とL*a*b*値との対応関係１〜３が作成される。
【０１２９】
〔２−５〕判断手順および作成処理手順の説明
次に、図５のステップＳ２３およびＳ２４における処理、つまり、本実施形態の判断手順および作成処理手順（判断部１４２および作成処理部１４３の動作）について、図９〜図１２を参照しながら、図８に示すフローチャート（ステップＳ５１〜Ｓ５８）に従って詳細に説明する。なお、図９〜図１２は、作成処理部１４３での内挿処理の具体例を説明するための図である。
【０１３０】
作成処理部１４３では、例えば、対応関係作成部１４１において１つのＲＧＢ値について作成された３つのL*a*b*値に対し補間演算処理（内挿処理）を行なうことにより、１つのL*a*b*値を算出する。
その算出処理（内挿処理）に先立ち、本実施形態では、まず、判断部１４２において、対応関係作成部１４１でのＲＧＢ信号の色変換結果（L*a*b*値）が、各色変換式に対応した領域に属しているか否かが判断される。そして、対応領域に属している場合、その色変換結果は正しいものと判断される。
【０１３１】
対応関係作成部１４１により作成された３つの対応関係１〜３の例と、判断部１４２の判断結果の例とを下表２に示す。なお、この表２において、判断部１４２が正しいと判断した色変換結果（L*a*b*値）の判断結果の欄には「○」が記入され、判断部１４２が正しくないと判断した色変換結果（L*a*b*値）の判断結果の欄には「×」が記入されている。
【０１３２】
【表２】

【０１３３】
ここで、例えば図９〜図１２に示すように、対応関係１の色変換結果をL*a*b*1（L1,a1,b1）とし、このL*a*b*1（L1,a1,b1）と領域１の境界面Ｂ１との距離をＤ１とする。同様に、対応関係２の色変換結果をL*a*b*2（L2,a2,b2）とし、このL*a*b*2（L2,a2,b2）と領域２の境界面Ｂ２との距離をＤ２とする。また、対応関係３の色変換結果をL*a*b*3（L3,a3,b3）とし、このL*a*b*3（L3,a3,b3）と領域３の境界面Ｂ３との距離をＤ３とする。
【０１３４】
このとき、判断部１４２においては、まず、L*a*b*1（L1,a1,b1）と領域１の境界面Ｂ１との距離Ｄ１が算出されるほか、そのL*a*b*1（L1,a1,b1）が領域１の内部にあると判断された場合、つまり、L*a*b*1（L1,a1,b1）が正しいと判断された場合には、カウンタ値としての「iNumber#Of#True」に１を加算する（ステップＳ５１）。
【０１３５】
同様に、L*a*b*2（L2,a2,b2）と領域２の境界面Ｂ２との距離Ｄ２を算出し、L*a*b*2（L2,a2,b2）が領域２の内部にあると判断した場合、つまり、L*a*b*2（L2,a2,b2）が正しいと判断した場合には、カウンタ値としての「iNumber#Of#True」に１を加算し（ステップＳ５２）、さらに、L*a*b*3（L3,a3,b3）と領域３の境界面Ｂ３との距離Ｄ３を算出し、L*a*b*3（L3,a3,b3）が領域３の内部にあると判断した場合、つまり、L*a*b*3（L3,a3,b3）が正しいと判断した場合には、カウンタ値としての「iNumber#Of#True」に１を加算する（ステップＳ５３）。
【０１３６】
そして、作成処理部１４３は、判断部１４２による判断結果（「iNumber#Of#True」の値、即ち表２の「○」の数）に応じて、以下のような内挿式(17)〜(20)を用いて、一つの色変換結果L*a*b*4（L4,a4,b4）を算出・決定する（ステップＳ５４〜Ｓ５８）。
【０１３７】
例えば図９に示すごとく３つの色変換結果L*a*b*1(L1,a1,b1)，L*a*b*2(L2,a2,b2)およびL*a*b*3(L3,a3,b3)がそれぞれ対応領域１〜３に属しており、「iNumber#Of#True」＝３となった場合〔例えば表２のＲＧＢ値(0,0,0)の色変換結果参照〕、下式(17)に示す内挿式を用い、距離Ｄ１〜Ｄ３に基づいて、一つのL*a*b*4（L4,a4,b4）が算出・決定される（ステップＳ５４）。
L4 = (L1*D1+L2*D2+L3*D3)/(D1+D2+D3)
a4 = (a1*D1+a2*D2+a3*D3)/(D1+D2+D3) (17)
b4 = (b1*D1+b2*D2+b3*D3)/(D1+D2+D3)
【０１３８】
例えば図１０に示すごとく２つの色変換結果L*a*b*1(L1,a1,b1)およびL*a*b*2(L2,a2,b2)がそれぞれ対応領域１および２に属しており、「iNumber#Of#True」＝２となった場合〔例えば表２のＲＧＢ値(0,0,15)の色変換結果参照〕、下式(18)に示す内挿式を用い、距離Ｄ１およびＤ２に基づいて、一つのL*a*b*4（L4,a4,b4）が算出・決定される（ステップＳ５５）。
L4 = (L1*D1+L2*D2)/(D1+D2)
a4 = (a1*D1+a2*D2)/(D1+D2) (18)
b4 = (b1*D1+b2*D2)/(D1+D2)
【０１３９】
例えば図１１に示すごとく一つの色変換結果L*a*b*1(L1,a1,b1)だけが対応領域１に属しており、「iNumber#Of#True」＝１となった場合〔例えば表２のＲＧＢ値(0,0,31)の色変換結果参照〕、下式(19)に示す通り、L*a*b*1(L1,a1,b1)がそのままL*a*b*4（L4,a4,b4）として決定される（ステップＳ５６）。
L4 =L1
a4 =a1 (19)
b4 =b1
【０１４０】
例えば図１２に示すごとく３つの色変換結果L*a*b*1(L1,a1,b1)，L*a*b*2(L2,a2,b2)およびL*a*b*3(L3,a3,b3)がいずれも対応領域１〜３に属しておらず、「iNumber#Of#True」＝０となった場合〔例えば表２のＲＧＢ値(255,255,255)の色変換結果参照〕、下式(20)に示す内挿式を用い、距離Ｄ１〜Ｄ３の逆数に基づいて、一つのL*a*b*4（L4,a4,b4）が算出・決定される（ステップＳ５７）。
L4=(L1*(1/D1)+L2*(1/D2)+L3*(1/D3))/((1/D1)+(1/D2)+(1/D3))
a4=(a1*(1/D1)+a2*(1/D2)+a3*(1/D3))/((1/D1)+(1/D2)+(1/D3)) (20)
b4=(b1*(1/D1)+b2*(1/D2)+b3*(1/D3))/((1/D1)+(1/D2)+(1/D3))
【０１４１】
上述したステップＳ５４〜Ｓ５７のいずれかの処理によって得られた、内挿処理後のL*a*b*値〔L*a*b*4（L4,a4,b4）〕が、ある一つのＲＧＢ値に対応する最終的な色変換結果としてメモリ２０に出力され、そのＲＧＢ値についての対応関係が決定される。
以上のような処理を、格子状のＲＧＢ値（格子点信号値）の全てについて実行することにより、メモリ２０上に、イメージスキャナ６０の色変換テーブル２２が作成される。
【０１４２】
なお、上述した作成処理部１４３では、正しいと判断された色変換結果が複数存在する場合、(17)式や(18)式に示すごとく距離Ｄ１〜Ｄ３に基づいて一つのL*a*b*4（L4,a4,b4）を算出・決定しているが、その距離が最も大きい色変換結果（境界から最も離れている色変換結果）を、最も信頼性の高いものであると判断し、L*a*b*4（L4,a4,b4）として選択・決定してもよい。
このような選択処理によれば、図９に示す例では、Ｄ１＞Ｄ３＞Ｄ２であるので、L*a*b*1(L1,a1,b1)が最終的な色変換結果L*a*b*4（L4,a4,b4）として選択・決定される。また、図１０に示す例では、Ｄ１＞Ｄ２であるので、L*a*b*1(L1,a1,b1)が最終的な色変換結果L*a*b*4（L4,a4,b4）として選択・決定される。
【０１４３】
同様に、上述した作成処理部１４３では、正しいと判断された色変換結果が存在しない場合、(20)式に示すごとく距離Ｄ１〜Ｄ３の逆数に基づいて一つのL*a*b*4（L4,a4,b4）を算出・決定しているが、その距離の逆数が最も大きい色変換結果（つまりは、距離が最も小さく境界に最も近い色変換結果）を、最も信頼性の高いものであると判断し、L*a*b*4（L4,a4,b4）として選択・決定してもよい。
このような選択処理によれば、図１２に示す例では、Ｄ３＜Ｄ１＜Ｄ２であるので、L*a*b*3(L3,a3,b3)が最終的な色変換結果L*a*b*4（L4,a4,b4）として選択・決定される。
【０１４４】
〔３〕本実施形態の作用効果の説明
このように、本実施形態によれば、イメージスキャナ６０の色変換テーブル２２の作成に際して、L*a*b*空間の分割領域に応じた色変換式（つまり分割領域毎に異なる色変換式）を用いることにより、色空間全体に亘る対応関係を精度良く算出することが可能になる。
【０１４５】
本実施形態の対応関係作成部１４１では、予め設定された各分割領域に属する、カラーパッチの分光反射率群に基づいて、分割領域毎に主成分ベクトルが算出される。色空間上の距離が近い２つの色についての分光反射率は類似しているので、前述のように同一の領域に属する分光反射率群に基づいて算出された主成分ベクトルは、その各分割領域に存在する分光反射率に対する寄与率が高い。このような主成分ベクトルを用いれば分光反射率を表現しやすくなり、結果として極めて精度良く対応関係を算出することができる。つまり、色空間全体から算出した主成分ベクトルを用いるよりも、各分割領域から算出した主成分ベクトルを用いた方が、精度良く対応関係を算出することができる。
【０１４６】
従って、手間をかけることなく、少数のカラーパッチ（色票データ：画像色信号，分光反射率）によって、精度良く色変換を行なえる色変換テーブル２２を作成することができる。また、このように作成された色変換テーブル２２を用いることにより、色空間全体に亘る色変換精度を向上させることができる。
【０１４７】
また、色変換テーブル作成部１４では、一つのＲＧＢ信号（格子点信号値）について３種類の色変換式により３組の対応関係１〜３を作成し、これらの対応関係１〜３に基づいて色変換テーブル２２が作成されるので、その対応関係１〜３を考慮しながらより精度の高い対応関係を算出することが可能になる。従って、より精度良く色変換を行なえる色変換テーブル２２を作成することができるとともに、このように作成された色変換テーブル２２を用いることにより、色空間全体に亘る色変換精度をさらに向上させることができる。
【０１４８】
さらに、本実施形態によれば、以下のような作用効果が得られる。
本実施形態では、変換先色空間として、人が見た際の色の差とその色空間での色の距離とがほぼ均等であるような均等色空間であるL*a*b*空間（ＣＩＥＬＡＢ）を用いることで、色変換を行なう際に、色空間内のどのような色でも、色空間内の距離を重み係数として用いることができる。
【０１４９】
L*a*b*空間を分割して得られる３つの分割領域１〜３が、相互に重複する領域を有することにより、色変換結果の階調性が保証される。つまり、重複領域に存在する色変換結果L*a*b*4（L4,a4,b4）として、分割領域毎の色変換式により算出された複数の色変換結果に対して補間演算を行なった結果を用いることにより、分割領域１〜３相互間の色の階調性が保証されることになる。
【０１５０】
色相角に応じてL*a*b*空間を分割して得られた領域１〜３を用いることで、より精度良く対応関係を算出することが可能になる。つまり、色票の分光反射率を、色相角領域毎に分類された分光反射率群に基づいて主成分ベクトルを算出した場合、その主成分ベクトルの、各色相角領域１〜３に属する分光反射率群に対する寄与率は、当然、高くなる。従って、そのような主成分ベクトルは、各色相角領域１〜３に属する分光反射率群をより正確に表現することができ、結果として精度良く対応関係を算出できる。このような効果は、様々な色相角をもつ色をもつ色票を用いた場合に顕著に得られる。
【０１５１】
特に、本実施形態のごとく、一次色（シアン，マゼンタ，イエロー）の色相角を各領域の中心としてL*a*b*空間を分割して３つの領域１〜３を得ることで、より精度良く対応関係を算出することが可能になる。一次色中心付近に存在する色の分光反射率の形状は、一次色の分光反射率の形状に類似している。そのため、一次色中心付近に存在する分光反射率群に基づいて主成分ベクトルを算出した場合、その主成分ベクトルの、一次色中心付近における分光反射率群に対する寄与率は、当然、高くなり、精度良く対応関係を算出できる。
【０１５２】
各色変換式により一つのＲＧＢ信号（格子点信号値）を変換して得られた色変換結果（L*a*b*値）が各色変換式に対応した領域に属している場合に、その色変換結果は正しいものと判断し、正しいと判断された色変換結果に基づいて色変換テーブル２２を作成することで、より精度良く対応関係を算出することが可能になる。色空間の各分割領域から得られた色変換式による色変換精度は、当然ながら、自らの領域外の色信号に対して色変換を行なう際には、低くなる。そのため、精度の低い色変換結果と、精度の高い色変換結果とを判別することで、より高い精度の色変換結果を取得することができる。
【０１５３】
一つのＲＧＢ信号（格子点信号値）について正しいと判断された色変換結果（L*a*b*値）が複数存在する場合、正しいと判断された複数（２または３）の色変換結果と各色変換結果の属する領域の境界（境界面）との距離Ｄ１〜Ｄ３に基づいて、そのＲＧＢ信号についての色変換結果L*a*b*4（L4,a4,b4）を算出することで、色変換結果の階調性が保証される。各色変換式を用いて算出された複数の色変換結果が大きく異なり、それぞれの領域境界面に対してずれる場合があるが、このような場合、上述のような色変換結果の算出処理（内挿処理）を行なうことにより、階調性が保証されることになる。
【０１５４】
また、一つのＲＧＢ信号（格子点信号値）について正しいと判断された色変換結果（L*a*b*値）が複数存在する場合、前記距離Ｄ１〜Ｄ３に基づいて、正しいと判断された複数の色変換結果のうちの一つを、そのＲＧＢ信号についての色変換結果L*a*b*4（L4,a4,b4）として選択することで、対応関係作成にかかる処理時間を短縮することができる。
【０１５５】
分割領域の数が多い場合などは、それに比例して、作成処理部１４３での色変換結果の算出処理にかかる時間が多くなる。そこで、処理時間を短縮するために、正しい色変換結果のうち、境界面からの距離が最も離れているものが、複数の正しい色変換結果の中でも最も信頼性の高いものであると判断して、その色変換結果を最終的な色変換結果L*a*b*4（L4,a4,b4）とする。
【０１５６】
一方、一つのＲＧＢ信号（格子点信号値）について正しいと判断された色変換結果（L*a*b*値）が存在しない場合、そのＲＧＢ信号について得られた色変換結果と各色変換結果の属する領域の境界（境界面）との距離Ｄ１〜Ｄ３の逆数に基づいて、そのＲＧＢ信号についての色変換結果L*a*b*4（L4,a4,b4）を算出することで、色変換結果が全て正しくないとしても、より信頼性のある色変換結果L*a*b*4（L4,a4,b4）を算出することが可能になる。つまり、全ての色変換結果が正しくない場合は、正しくないものの中でも各境界面Ｂ１〜Ｂ３との距離が近いものの方が信頼性が高いものと判断している。
【０１５７】
また、一つのＲＧＢ信号（格子点信号値）について正しいと判断された色変換結果（L*a*b*値）が存在しない場合、前記距離Ｄ１〜Ｄ３の逆数に基づいて、色変換結果のうちの一つを、そのＲＧＢ信号についての色変換結果L*a*b*4（L4,a4,b4）として選択することで、色変換結果が全て正しくないとしても、より信頼性の高い色変換結果L*a*b*4(L4,a4,b4)を短時間で得ることが可能になる。この場合、境界面に最も近い色変換結果が最も信頼性が高いと判断し、前述した通り、距離の逆数が最も大きい色変換結果を選択する。
【０１５８】
〔４〕本実施形態の変形例の説明
次に、本実施形態の変形例について図１３〜図１７を参照しながら説明する。
〔４−１〕第１変形例の説明
上述した実施形態では、L*a*b*空間の全体を色相角に応じて分割して得られた３つの領域１〜３が、予め定義されているが、第１変形例では、L*a*b*空間の全体を彩度に応じて２つの領域１，２に分割した場合について、図１３〜図１６を参照しながら説明する。なお、図１３は、本実施形態の第１変形例における、彩度に応じた領域分割例を説明するための図、図１４〜図１６は本実施形態の第１実施形態における内挿処理の具体例を説明するための図である。
【０１５９】
この第１変形例では、L*a*b*空間の全体を分割して得られる２つの領域１および領域２が、予め定義され設定されている（分割ステップ）。これら２つの領域１，２は、例えば図１３に示すごとく、L*a*b*空間を、彩度C*〔=(a*²＋b*²)^1/2〕に応じた境界（実際にはL*軸に平行な境界面）Ｂ１１およびＢ１２により区分けして得られたものである。
【０１６０】
図１３に示すように、領域１は、L*a*b*空間において彩度C*が所定値C*1（例えば10.0）以上の領域であり、領域２は、L*a*b*空間において彩度C*が所定値C*2(例えば20.0)未満の領域である。従って、図１３に示す領域１および領域２には、相互に重複する幅C*2−C*1（ここでは10.0）の領域が付加・形成されている。
【０１６１】
このように、L*a*b*空間の全体を彩度C*に応じて分割して得られた領域１および領域２を用いる場合も、基本的には、図１〜図１２を参照しながら上述した色変換テーブル作成装置１により、上述した手順と同様の手順（図５のステップＳ１１〜Ｓ２５）に従って、イメージスキャナ６０の色変換テーブル２２が作成される。
【０１６２】
以下に、その色変換テーブル２２の作成手順を簡単に説明する。
〔４−１−１〕分光特性推定手順の説明
分光特性推定部１１（ステップＳ１６）では、項目〔２−２〕において前述した手順と同様にして、数十種の色信号および数十種の分光反射率を用いて分光特性が推定される〔上式(10)参照〕。
【０１６３】
〔４−１−２〕分光反射率分類手順の説明
図１３に示すような領域分割を行なった第１変形例における分光反射率分類部１２（ステップＳ１８）では、各カラーパッチについてのL*a*b*値が図１３に示す領域１および領域２のいずれに属するかが判断され、各カラーパッチについての分光反射率が、その領域に応じて２つのグループ１，２に分類される。
つまり、そのL*a*b*値が領域１に属する（C*≧C*1）と判断された分光反射率は、グループ１に分類され、そのL*a*b*値が領域２に属する（C*＜C*2）と判断された分光反射率は、グループ２に分類される。
【０１６４】
〔４−１−３〕色変換式作成手順および対応関係作成手順の説明
対応関係作成部１４１（ステップＳ２１）では、項目〔２−４〕において前述した手順とほぼ同様にして、格子状のＲＧＢ値に対応するL*a*b*値を算出する。
ただし、第１変形例では、格子状のＲＧＢ値と、イメージスキャナ６０の分光特性と、グループ１の分光反射率から算出された主成分ベクトルとに基づいて、色空間全体に亘る対応関係１が算出され、同様に、格子状のＲＧＢ値と、イメージスキャナ６０の分光特性と、グループ２の分光反射率から算出された主成分ベクトルとに基づいて、色空間全体に亘る対応関係２が算出される。
【０１６５】
これにより、１つの格子状のＲＧＢ値に対して、２種類の色変換式を用い、色変換結果としての２つのL*a*b*値（L*a*b*1，L*a*b*2）が算出される。
１つのグループ（分光反射率群）に対する対応関係の算出手順は、図７に示すフローチャート（ステップＳ４１〜Ｓ４７）を参照しながら前述した通りであるので、ここでは、その説明を省略する。
【０１６６】
〔４−１−４〕判断手順および作成処理手順の説明
作成処理部１４３（ステップＳ２４）では、例えば、対応関係作成部１４１において１つのＲＧＢ値について作成された２つのL*a*b*値から、より信頼性の高いものを選択して、１つのL*a*b*値を決定する。
選択処理に先立って、第１変形例においても、まず、判断部１４２（ステップＳ２３）において、対応関係作成部１４１でのＲＧＢ信号の色変換結果（L*a*b*値）が、各色変換式に対応した領域に属しているか否かが判断される。そして、対応領域に属している場合、その色変換結果は正しいものと判断される。
【０１６７】
ここで、例えば図１４〜図１６に示すように、対応関係１の色変換結果をL*a*b*1（L1,a1,b1）とし、このL*a*b*1（L1,a1,b1）と領域１の境界面Ｂ１１との距離をＤ１とし、同様に、対応関係２の色変換結果をL*a*b*2（L2,a2,b2）とし、このL*a*b*2（L2,a2,b2）と領域２の境界面Ｂ１２との距離をＤ２とする。
そして、作成処理部１４３は、判断部１４２による判断結果と距離Ｄ１，Ｄ２とに基づいて、以下のようにして、一つの色変換結果を選択・決定する。
【０１６８】
２つの色変換結果L*a*b*1(L1,a1,b1)およびL*a*b*2(L2,a2,b2)がそれぞれ対応領域１，２に属しており、２つとも正しいと判断された場合、境界Ｂ１１，Ｂ１２との距離Ｄ１，Ｄ２の大きい方が、より信頼性が高いものと判断する。つまり、以下の２つのケースに従って処理を行なう。
▲１▼Ｄ１＞Ｄ２の場合
色変換結果L*a*b*1(L1,a1,b1)を選択する。
▲２▼Ｄ１＜Ｄ２の場合
色変換結果L*a*b*2(L2,a2,b2)を選択する。
【０１６９】
従って、図１４に示す例では、２つの色変換結果L*a*b*1(L1,a1,b1)およびL*a*b*2(L2,a2,b2)がそれぞれ対応領域１，２に属しているが、Ｄ２＞Ｄ１であるので、上記ケース▲２▼に該当し、色変換結果L*a*b*2(L2,a2,b2)が選択される。
なお、Ｄ１＝Ｄ２の場合には、例えば、上式(18)のような内挿式を適用し、２つの色変換結果L*a*b*1(L1,a1,b1)およびL*a*b*2(L2,a2,b2)の平均値を算出し、その平均値を最終的な色変換結果として選択する。
【０１７０】
２つの色変換結果L*a*b*1(L1,a1,b1)およびL*a*b*2(L2,a2,b2)のうちの一方だけが正しいと判断された場合は、その正しいと判断された一つの色変換結果が、そのまま最終的な色変換結果として選択される。
従って、図１５に示す例では、色変換結果L*a*b*1(L1,a1,b1)だけがそれぞれ対応領域１に属しているので、色変換結果L*a*b*1(L1,a1,b1)が選択される。
【０１７１】
２つの色変換結果L*a*b*1(L1,a1,b1)およびL*a*b*2(L2,a2,b2)がいずれも対応領域１，２に属しておらず、２つとも正しくないと判断された場合、境界Ｂ１１，Ｂ１２との距離Ｄ１，Ｄ２の小さい方が、より信頼性が高いものと判断する。つまり、以下の２つのケースに従って処理を行なう。
▲１▼Ｄ１＜Ｄ２の場合（つまり１／Ｄ１＞１／Ｄ２の場合）
色変換結果L*a*b*1(L1,a1,b1)を選択する。
▲２▼Ｄ１＞Ｄ２の場合（つまり１／Ｄ１＜１／Ｄ２の場合）
色変換結果L*a*b*2(L2,a2,b2)を選択する。
【０１７２】
従って、図１６に示す例では、２つの色変換結果L*a*b*1(L1,a1,b1)およびL*a*b*2(L2,a2,b2)がいずれも対応領域１，２に属していないが、Ｄ１＜Ｄ２であるので、上記ケース▲１▼に該当し、色変換結果L*a*b*1(L1,a1,b1)が選択される。
なお、Ｄ１＝Ｄ２（１／Ｄ１＝１／Ｄ２）の場合には、例えば、２つの色変換結果L*a*b*1(L1,a1,b1)およびL*a*b*2(L2,a2,b2)の平均値を算出し、その平均値を最終的な色変換結果として選択する。
【０１７３】
上述のようにして選択されたL*a*b*値が、ある一つのＲＧＢ値に対応する最終的な色変換結果としてメモリ２０に出力され、そのＲＧＢ値についての対応関係が決定される。
以上のような処理を、格子状のＲＧＢ値（格子点信号値）の全てについて実行することにより、メモリ２０上に、イメージスキャナ６０の色変換テーブル２２が作成される。
【０１７４】
〔４−１−５〕第１変形例の作用効果の説明
このように、第１変形例によれば、上述した本実施形態と同様の作用効果が得られるほか、この第１変形例では、彩度に応じてL*a*b*空間の全体を分割して得られた領域１，２を用いることで、精度良く対応関係を算出することが可能になる。つまり、色票の分光反射率を、彩度領域毎に分類された分光反射率群に基づいて主成分ベクトルを算出した場合、その主成分ベクトルの、各彩度領域１，２に属する分光反射率群に対する寄与率は、当然、高くなる。従って、そのような主成分ベクトルは、各彩度領域１，２に属する分光反射率群をより正確に表現することができ、結果として精度良く対応関係を算出できる。このような効果は、様々な彩度をもつ色をもつ色票を用いた場合に顕著に得られる。
【０１７５】
〔４−２〕第２変形例の説明
上述した実施形態やその第１変形例では、L*a*b*空間の全体を色相角や彩度に応じて分割して得られた領域１〜３が、予め定義されているが、第２変形例では、L*a*b*空間の全体を明度に応じて２つの領域１，２に分割した場合について、図１７を参照しながら説明する。なお、図１７は、本実施形態の第２変形例における、明度に応じた領域分割例を説明するための図である。
【０１７６】
この第２変形例でも、第１変形例と同様、L*a*b*空間の全体を分割して得られる２つの領域１および領域２が予め定義され設定されている（分割ステップ）。これら２つの領域１，２は、例えば図１７に示すごとく、L*a*b*空間を、明度L*に応じた境界（実際にはa*b*面に平行な境界面）Ｂ２１およびＢ２２により区分けして得られたものである。
【０１７７】
図１７に示すように、領域１は、L*a*b*空間において明度L*が所定値L*1以上の領域であり、領域２は、L*a*b*空間において明度L*が所定値L*2未満の領域である。従って、図１７に示す領域１および領域２には、相互に重複する幅L*2−L*1の領域が付加・形成されている。
【０１７８】
このように、L*a*b*空間の全体を明度L*に応じて分割して得られた領域１および領域２を用いる場合も、基本的には、図１〜図１２を参照しながら上述した色変換テーブル作成装置１により、上述した手順と同様の手順（図５のステップＳ１１〜Ｓ２５）に従って、イメージスキャナ６０の色変換テーブル２２が作成される。このとき、分割領域の数が２つという点では上述した第１変形例と同じであり、この第２変形例では、本実施形態およびその第１変形例で説明した手順とほぼ同様にして色変換テーブル２２が作成されることになる。
【０１７９】
このように、第２変形例によれば、上述した本実施形態と同様の作用効果が得られるほか、この第２変形例では、明度に応じてL*a*b*空間の全体を分割して得られた領域１，２を用いることで、精度良く対応関係を算出することが可能になる。つまり、色票の分光反射率を、明度領域毎に分類された分光反射率群に基づいて主成分ベクトルを算出した場合、その主成分ベクトルの、各明度領域１，２に属する分光反射率群に対する寄与率は、当然、高くなる。従って、そのような主成分ベクトルは、各明度領域１，２に属する分光反射率群をより正確に表現することができ、結果として精度良く対応関係を算出できる。このような効果は、様々な明度をもつ色をもつ色票を用いた場合に顕著に得られる。
【０１８０】
〔４−３〕第３変形例の説明
ところで、上述した実施形態では、色変換テーブル２２は、ＲＧＢ信号とL*a*b*値とを対応させるものとして作成したが、第３変形例では、ＲＧＢ信号を、上式(15)によって算出された分光反射率（色変換結果としてのL*a*b*値に対応する分光反射率）に対応させる色変換テーブルを作成するようになっている。
【０１８１】
このように、ＲＧＢ信号を分光反射率に対応させる色変換テーブルを作成する場合も、基本的には、図１〜図１２を参照しながら上述した色変換テーブル作成装置１により、上述した手順と同様の手順（図５のステップＳ１１〜Ｓ２５）に従って、イメージスキャナ６０の色変換テーブルが作成される。
【０１８２】
ただし、対応関係作成部１４１（ステップＳ２１）で用いられる色変換式は上式(14)，(15)となる。また、判断部１４２（ステップＳ２３）では、分光反射率の正当性を判断するために、ＲＧＢ値に対応するL*a*b*値が必要となるので、この第３変形例においても、上述した実施形態と同様にして、ＲＧＢ値に対応するL*a*b*値が、色変換結果として、対応関係作成部１４１もしくは判断部１４２において算出される。
【０１８３】
また、作成処理部１４３（ステップＳ２４）で補間演算処理（内挿処理）を行なう場合には、上式(17)〜(20)式に代えて、下式(21)〜(24)が用いられる。
ここでは、ある一つのＲＧＢ値に対する３つの対応関係１〜３として、それぞれ、３つの分光反射率REF1(REF1₃₈₀,REF1₃₉₀,…,REF1₇₃₀)，REF2(REF2₃₈₀,REF2₃₉₀,…,REF2₇₃₀)およびREF3(REF3₃₈₀,REF3₃₉₀,…,REF3₇₃₀)が得られているものとする。また、これらの分光反射率に対応するL*a*b*値（色変換結果）を、それぞれ、L*a*b*1(L1,a1,b1)，L*a*b*2(L2,a2,b2)およびL*a*b*3(L3,a3,b3)とする。
【０１８４】
さらに、図９〜図１２に示した例と同様、L*a*b*1（L1,a1,b1）と領域１の境界面Ｂ１との距離をＤ１、L*a*b*2（L2,a2,b2）と領域２の境界面Ｂ２との距離をＤ２、L*a*b*3（L3,a3,b3）と領域３の境界面Ｂ３との距離をＤ３とする。
【０１８５】
例えば図９に示すごとくL*a*b*1(L1,a1,b1)，L*a*b*2(L2,a2,b2)およびL*a*b*3(L3,a3,b3)がそれぞれ対応領域１〜３に属しており、L*a*b*1(L1,a1,b1)，L*a*b*2(L2,a2,b2)およびL*a*b*3(L3,a3,b3)がいずれも正しいと判断された場合、下式(21)に示す内挿式を用い、距離Ｄ１〜Ｄ３に基づいて、一つの分光反射率REF4(REF4₃₈₀,REF4₃₉₀,…,REF4₇₃₀)が算出・決定される。

【０１８６】
例えば図１０に示すごとくL*a*b*1(L1,a1,b1)およびL*a*b*2(L2,a2,b2)がそれぞれ対応領域１および２に属しており、L*a*b*1(L1,a1,b1)およびL*a*b*2(L2,a2,b2)が正しいと判断された場合、下式(22)に示す内挿式を用い、距離Ｄ１およびＤ２に基づいて、一つの分光反射率REF4(REF4₃₈₀,REF4₃₉₀,…,REF4₇₃₀)が算出・決定される。

【０１８７】
例えば図１１に示すごとくL*a*b*1(L1,a1,b1)だけが対応領域１に属しており、L*a*b*1(L1,a1,b1)だけが正しいと判断された場合、下式(23)に示す通り、REF1(REF1₃₈₀,REF1₃₉₀,…,REF1₇₃₀)がそのままREF4(REF4₃₈₀,REF4₃₉₀,…,REF4₇₃₀)が決定される。

【０１８８】
例えば図１２に示すごとくL*a*b*1(L1,a1,b1)，L*a*b*2(L2,a2,b2)およびL*a*b*3(L3,a3,b3)がいずれも対応領域１〜３に属しておらず、L*a*b*1(L1,a1,b1)，L*a*b*2(L2,a2,b2)およびL*a*b*3(L3,a3,b3)がいずれも正しくないと判断された場合、下式(24)に示す内挿式を用い、距離Ｄ１〜Ｄ３の逆数に基づいて、一つの分光反射率REF4(REF4₃₈₀,REF4₃₉₀,…,REF4₇₃₀)が算出・決定される。

【０１８９】
上述のような内挿処理によって得られた分光反射率REF4(REF4₃₈₀,REF4₃₉₀,…,REF4₇₃₀)が、ある一つのＲＧＢ値に対応する最終的な分光反射率としてメモリ２０に出力され、そのＲＧＢ値についての対応関係が決定される。
以上のような処理を、格子状のＲＧＢ値（格子点信号値）の全てについて実行することにより、メモリ２０上に、ＲＧＢ値と分光反射率とを対応させる色変換テーブルが作成される。
【０１９０】
なお、ここでは、正しいと判断された色変換結果（L*a*b*値）が複数存在する場合、(21)式や(22)式に示すごとく距離Ｄ１〜Ｄ３に基づいて一つの分光反射率REF4(REF4₃₈₀,REF4₃₉₀,…,REF4₇₃₀)を算出・決定しているが、その距離が最も大きい色変換結果（境界から最も離れている色変換結果）を、最も信頼性の高いものであると判断し、分光反射率REF4(REF4₃₈₀,REF4₃₉₀,…,REF4₇₃₀)として選択・決定してもよい。
【０１９１】
このような選択処理によれば、図９に示す例では、Ｄ１＞Ｄ３＞Ｄ２であるので、REF1(REF1₃₈₀,REF1₃₉₀,…,REF1₇₃₀)が最終的な分光反射率REF4(REF4₃₈₀,REF4₃₉₀,…,REF4₇₃₀)として選択・決定される。また、図１０に示す例では、Ｄ１＞Ｄ２であるので、REF1(REF1₃₈₀,REF1₃₉₀,…,REF1₇₃₀)が最終的な分光反射率REF4(REF4₃₈₀,REF4₃₉₀,…,REF4₇₃₀)として選択・決定される。
【０１９２】
同様に、正しいと判断された色変換結果（L*a*b*値）が存在しない場合、(20)式に示すごとく距離Ｄ１〜Ｄ３の逆数に基づいて一つの分光反射率REF4(REF4₃₈₀,REF4₃₉₀,…,REF4₇₃₀)を算出・決定しているが、その距離の逆数が最も大きい色変換結果（つまりは、距離が最も小さく境界に最も近い色変換結果）を、最も信頼性の高いものであると判断し、分光反射率REF4(REF4₃₈₀,REF4₃₉₀,…,REF4₇₃₀)として選択・決定してもよい。
【０１９３】
このような選択処理によれば、図１２に示す例では、Ｄ３＜Ｄ１＜Ｄ２であるので、REF3(REF3₃₈₀,REF3₃₉₀,…,REF3₇₃₀)が最終的な分光反射率REF4(REF4₃₈₀,REF4₃₉₀,…,REF4₇₃₀)として選択・決定される。
【０１９４】
上述のように、第３変形例によれば、上述した本実施形態と同様の作用効果が得られるほか、この第３変形例では、ＲＧＢ信号を色変換結果としての分光反射率REF(REF₃₈₀,REF₃₉₀,…,REF₇₃₀)に対応させる色変換テーブルが作成されるので、色変換結果を分光反射率にすることができ、イメージスキャナ６０を簡易測色器の代用とすることが可能になる。
【０１９５】
一つのＲＧＢ信号（格子点信号値）について正しいと判断された色変換結果（L*a*b*値）が複数存在する場合、正しいと判断された、２または３の色変換結果と各色変換結果の属する領域の境界（境界面）との距離Ｄ１〜Ｄ３に基づいて、そのＲＧＢ信号についての分光反射率REF4(REF4₃₈₀,REF4₃₉₀,…,REF4₇₃₀)を算出することで、分光反射率の階調性が保証される。
【０１９６】
また、一つのＲＧＢ信号（格子点信号値）について正しいと判断された色変換結果（L*a*b*値）が複数存在する場合、正しいと判断された複数の色変換結果と各色変換結果の属する領域の境界（境界面）との距離Ｄ１〜Ｄ３に基づいて、正しいと判断された複数の色変換結果に応じた分光反射率のうちの一つを、そのＲＧＢ信号についての分光反射率REF4(REF4₃₈₀,REF4₃₉₀,…,REF4₇₃₀)として選択することで、対応関係作成にかかる処理時間を短縮することができる。
【０１９７】
一方、一つのＲＧＢ信号（格子点信号値）について正しいと判断された色変換結果（L*a*b*値）が存在しない場合、そのＲＧＢ信号について得られた複数の色変換結果と各色変換結果の属する領域の境界（境界面）との距離Ｄ１〜Ｄ３の逆数に基づいて、そのＲＧＢ信号についての分光反射率REF4(REF4₃₈₀,REF4₃₉₀,…,REF4₇₃₀)を算出することで、色変換結果が全て正しくないとしても、より信頼性のある分光反射率REF4(REF4₃₈₀,REF4₃₉₀,…,REF4₇₃₀)を算出することが可能になる。
【０１９８】
また、一つのＲＧＢ信号（格子点信号値）について正しいと判断された色変換結果（L*a*b*値）が存在しない場合、その一の色信号について得られた複数の色変換結果と各色変換結果の属する領域の境界（境界面）との距離Ｄ１〜Ｄ３の逆数に基づいて、複数の色変換結果に応じた分光反射率のうちの一つを、一の色信号についての分光反射率REF4(REF4₃₈₀,REF4₃₉₀,…,REF4₇₃₀)として選択することで、色変換結果が全て正しくないとしても、より信頼性の高い分光反射率REF4(REF4₃₈₀,REF4₃₉₀,…,REF4₇₃₀)を短時間で得ることが可能になる。
【０１９９】
〔４−４〕第４変形例の説明
なお、上述した本実施形態や第１〜第３変形例では、作成処理部１４３において色変換結果や分光反射率の内挿処理もしくは選択処理を行なう際に、色変換結果（L*a*b*値）と各領域の境界との距離を用いているが、その距離に代えて、その距離を変換して得られるΔＥ９４の値を用いてもよい。
【０２００】
この場合、人間の目から見て、より正確な色変換を行なうことのできる色変換テーブルが作成されることになる。ΔＥ９４は、L*a*b*空間上の距離（ΔＥａｂ）よりも人間の見た目に近い色差をもつため、このΔＥ９４を用いた色変換結果は、人間の目から見て、より正確なものになる。
【０２０１】
〔４−５〕その他
なお、本発明は上述した実施形態に限定されるものではなく、本発明の趣旨を逸脱しない範囲で種々変形して実施することができる。
例えば、上述した実施形態では、変換元色空間がＲＧＢ空間であり変換先色空間L*a*b*空間である場合について説明したが、本発明は、これに限定されるものではなく、各種色空間における色信号どうしを対応させる色変換テーブルを作成する際に用いることにより、上述した実施形態と同様の作用効果を得ることができる。
【０２０２】
上述した実施形態では、変換先色空間を３分割あるいは２分割した場合について説明したが、本発明はこれに限定されるものではなく、４以上の領域に分割してもよい。また、上述した実施形態では、変換先色空間を、色相角のみに応じて分割した場合、彩度のみに応じて分割した場合、明度のみに応じて分割した場合について説明したが、色相角，彩度および明度のうちの少なくとも２つのパラメータを組み合わせ、これらのパラメータに応じて変換先色空間を複数の領域に分割してもよい。
【０２０３】
さらに、上述した実施形態では、主成分ベクトルや分光特性や分光反射率が、可視光波長範囲380〜730nmにおける１０nm毎の各波長（380,390,…,730nm）に対応して得られる３６個の離散データである場合について説明したが、本発明は、主成分ベクトルや分光特性や分光反射率を、上述のような離散データに限定するものではない。
【０２０４】
また、上述した実施形態では、判断部１４２においては、対応関係作成部１４１で算出されたL*a*b*値を用いて判断を行なっているが、判断部１４２が、L*a*b*値を算出するように構成してもよい。
【０２０５】
〔５〕付記
〔付記１〕 色入力機器から出力された、該色入力機器の色空間（以下、変換元色空間）における色信号と、該変換元色空間とは異なる色空間（以下、変換先色空間という）における色信号とを対応させる色変換テーブルを作成する方法であって、
該変換先色空間の全体を分割して得られる複数の領域を、予め定義しておく分割ステップと、
該分割ステップで定義された該複数の領域にそれぞれ対応した複数の色変換式を用いて、該色変換テーブルを作成する色変換テーブル作成ステップとを有することを特徴とする、色変換テーブル作成方法。
【０２０６】
〔付記２〕 色票上の複数種の色領域を該色入力機器で読み込み、各色領域に対応する、該変換元色空間での色信号を該色入力機器から出力する読込ステップと、
該複数種の色領域を測色器で測色し、各色領域に対応する分光反射率を該測色器から出力する測色ステップと、
該分光反射率に対応する、該変換先色空間での色信号が該色変換先色空間において該複数の領域のうちのどの領域に属するかに応じて、該分光反射率を分類する分類ステップと、
該読込ステップで該色入力機器から出力された該色信号と該測色ステップで該測色器から出力された該分光反射率とに基づいて、該色入力機器の分光特性を推定する分光特性推定ステップとをさらに有し、
該色変換テーブル作成ステップにおいて、該分類ステップで分類された該分光反射率と該分光特性推定ステップで推定された前記分光特性とに基づいて、該変換先色空間での前記領域毎に該色変換式を作成することを特徴とする、付記１記載の色変換テーブル作成方法。
【０２０７】
〔付記３〕 該色変換テーブル作成ステップが、
各領域に応じた該色変換式を用いて、該色変換式毎に、該変換元色空間における色信号と該変換先色空間における色信号との対応関係を作成する対応関係作成ステップと、
該対応関係作成ステップで該色変換式毎に作成された該対応関係に基づいて、該色変換テーブルを作成する作成処理ステップとを有していることを特徴とする、付記１または付記２に記載の色変換テーブル作成方法。
【０２０８】
〔付記４〕 該変換先色空間が均等色空間であることを特徴とする、付記１記載の色変換テーブル作成方法。
〔付記５〕 該複数の領域が、相互に重複する領域を有していることを特徴とする、付記１記載の色変換テーブル作成方法。
〔付記６〕 該分割ステップにおいて、色相角に応じて該変換先色空間の全体を分割することにより、該複数の領域が得られることを特徴とする、付記１記載の色変換テーブル作成方法。
【０２０９】
〔付記７〕 該分割ステップにおいて、彩度に応じて該変換先色空間の全体を分割することにより、該複数の領域が得られることを特徴とする、付記１記載の色変換テーブル作成方法。
〔付記８〕 該分割ステップにおいて、明度に応じて該変換先色空間の全体を分割することにより、該複数の領域が得られることを特徴とする、付記１記載の色変換テーブル作成方法。
【０２１０】
〔付記９〕 該色変換テーブル作成ステップが、
該対応関係作成ステップにおいて該対応関係を作成する際に各色変換式により該変換元色空間の一の色信号を変換した色変換結果が、各色変換式に対応した該領域に属している場合に、該色変換結果は正しいものと判断する判断ステップをさらに有し、
該作成処理ステップにおいて、該判断ステップで正しいと判断された該色変換結果に基づいて該色変換テーブルを作成することを特徴とする、付記３記載の色変換テーブル作成方法。
【０２１１】
〔付記１０〕 該作成処理ステップにおいて、該一の色信号について該判断ステップで正しいと判断された色変換結果が複数存在する場合、正しいと判断された該複数の色変換結果と各色変換結果の属する該領域の境界との距離に関する値に基づいて、該一の色信号についての色変換結果を算出することを特徴とする、付記９記載の色変換テーブル作成方法。
【０２１２】
〔付記１１〕 該作成処理ステップにおいて、該一の色信号について該判断ステップで正しいと判断された色変換結果が複数存在する場合、正しいと判断された該複数の色変換結果と各色変換結果の属する該領域の境界との距離に関する値に基づいて、正しいと判断された該複数の色変換結果のうちの一つを、該一の色信号についての色変換結果として選択することを特徴とする、付記９記載の色変換テーブル作成方法。
【０２１３】
〔付記１２〕 該作成処理ステップにおいて、該一の色信号について該判断ステップで正しいと判断された色変換結果が存在しない場合、該対応関係作成ステップで該一の色信号について得られた複数の色変換結果と各色変換結果の属する該領域の境界との距離に関する値の逆数に基づいて、該一の色信号についての色変換結果を算出することを特徴とする、付記９記載の色変換テーブル作成方法。
【０２１４】
〔付記１３〕 該作成処理ステップにおいて、該一の色信号について該判断ステップで正しいと判断された色変換結果が存在しない場合、該対応関係作成ステップで該一の色信号について得られた複数の色変換結果と各色変換結果の属する該領域の境界との距離に関する値の逆数に基づいて、該複数の色変換結果のうちの一つを、該一の色信号についての色変換結果として選択することを特徴とする、付記９記載の色変換テーブル作成方法。
【０２１５】
〔付記１４〕 該色変換テーブルが、該変換元色空間の色信号を、該色変換先空間の色信号としての、該色変換結果に応じた分光反射率に対応させるものであることを特徴とする、付記９記載の色変換テーブル作成方法。
〔付記１５〕 該作成処理ステップにおいて、該一の色信号について該判断ステップで正しいと判断された色変換結果が複数存在する場合、正しいと判断された該複数の色変換結果と各色変換結果の属する該領域の境界との距離に関する値に基づいて、該一の色信号についての分光反射率を算出することを特徴とする、付記１４記載の色変換テーブル作成方法。
【０２１６】
〔付記１６〕 色入力機器から出力された、該色入力機器の色空間（以下、変換元色空間）における色信号と、該変換元色空間とは異なる色空間（以下、変換先色空間という）における色信号とを対応させる色変換テーブルを作成する装置であって、
該変換先色空間を分割して得られる複数の領域にそれぞれ対応した複数の色変換式を用いて、該色変換テーブルを作成する色変換テーブル作成部をそなえたことを特徴とする、色変換テーブル作成装置。
【０２１７】
〔付記１７〕 色票上の複数種の色領域を該色入力機器により読み込んで得られた、各色領域に対応する、該変換元色空間での色信号を入力する入力部と、該複数種の色領域を測色し、各色領域に対応する分光反射率を得る測色器と、該分光反射率に対応する、該変換先色空間での色信号が該色変換先色空間において該複数の領域のうちのどの領域に属するかに応じて、該分光反射率を分類する分類部と、
該入力部から入力された該色信号と該測色器により得られた該分光反射率とに基づいて、該色入力機器の分光特性を推定する分光特性推定部とをさらにそなえ、
該色変換テーブル作成部が、該分類部により分類された該分光反射率と該分光特性推定部により推定された前記分光特性とに基づいて、該変換先色空間での前記領域毎に該色変換式を作成することを特徴とする、付記１６記載の色変換テーブル作成装置。
【０２１８】
〔付記１８〕 該色変換テーブル作成部が、
各領域に応じた該色変換式を用いて、該色変換式毎に、該変換元色空間における色信号と該変換先色空間における色信号との対応関係を作成する対応関係作成部と、
該対応関係作成部で該複数の色変換式によりそれぞれ作成された複数組の該対応関係に基づいて、一組の対応関係を得て、該色変換テーブルを作成する作成処理部とをそなえていることを特徴とする、付記１６または付記１７に記載の色変換テーブル作成装置。
【０２１９】
〔付記１９〕 色入力機器から出力された、該色入力機器の色空間（以下、変換元色空間）における色信号と、該変換元色空間とは異なる色空間（以下、変換先色空間という）における色信号とを対応させる色変換テーブルを作成する機能をコンピュータにより実現させるための色変換テーブル作成プログラムを記録したコンピュータ読取可能な記録媒体であって、
該色変換テーブル作成プログラムが、
該変換先色空間を分割して得られる複数の領域にそれぞれ対応した複数の色変換式を用いて、該色変換テーブルを作成する色変換テーブル作成部として、該コンピュータを機能させることを特徴とする、色変換テーブル作成プログラムを記録したコンピュータ読取可能な記録媒体。
【０２２０】
〔付記２０〕 該色変換テーブル作成プログラムが、さらに、
色票上の複数種の色領域を測色器により測色して該色領域毎に得られた分光反射率に対応する、該変換先色空間での色信号が、該色変換先色空間において該複数の領域のうちのどの領域に属するかに応じて、該分光反射率を分類する分類部、および、
該複数種の色領域を該色入力機器により読み込んで該色領域毎に得られた該変換元色空間での色信号と、該測色器により得られた該分光反射率とに基づいて、該色入力機器の分光特性を推定する分光特性推定部として、該コンピュータを機能させるとともに、
該コンピュータを該色変換テーブル作成部として機能させる際に、該分類部により分類された該分光反射率と該分光特性推定部により推定された前記分光特性とに基づいて、該変換先色空間での前記領域毎に該色変換式を作成させることを特徴とする、付記１９記載の色変換テーブル作成プログラムを記録したコンピュータ読取可能な記録媒体。
【０２２１】
〔付記２１〕 該色変換テーブル作成プログラムが、
該コンピュータを該色変換テーブル作成部として機能させる際に、
各領域に応じた該色変換式を用いて、該色変換式毎に、該変換元色空間における色信号と該変換先色空間における色信号との対応関係を作成する対応関係作成部、および、
該対応関係作成部で該複数の色変換式によりそれぞれ作成された複数組の該対応関係に基づいて、一組の対応関係を得て、該色変換テーブルを作成する作成処理部として、該コンピュータを機能させることを特徴とする、付記１９または付記２０に記載の色変換テーブル作成プログラムを記録したコンピュータ読取可能な記録媒体。
【０２２２】
【発明の効果】
以上詳述したように、本発明の色変換テーブル作成方法（請求項１）および色変換テーブル作成装置（請求項２）並びに色変換テーブル作成プログラムを記録したコンピュータ読取可能な記録媒体（請求項３）によれば、色入力機器の色変換テーブルの作成に際して、変換先色空間の分割領域に応じた色変換式（つまり分割領域毎に異なる色変換式）を用いることにより、色空間全体に亘る対応関係を精度良く算出することが可能になる。
【０２２３】
従って、手間をかけることなく、少数の色票データによって、精度良く色変換を行なえる色変換テーブルを作成することができる。また、このように作成された色変換テーブルを用いることにより、色空間全体に亘る色変換精度を向上させることができる。
【０２２４】
このとき、各色変換式を、色入力機器の分光特性と色入力機器で読み取った色票の分光反射率（分割領域毎に分類されたもの）とに基づいて作成する。即ち、各分割領域に含まれる色領域（カラーパッチ）の分光反射率から主成分ベクトルを算出する。色空間上で距離の近い２つの色の分光反射率は類似しているので、主成分ベクトルにより分光反射率を表現しやすくなり、結果としてさらに精度良く対応関係を算出することができる（請求項１〜３）。
【０２２５】
また、一の色信号について複数の色変換式により複数組の対応関係を作成し、その複数組の対応関係に基づいて色変換テーブルが作成されるので、その複数組の対応関係を考慮しながら、より精度の高い対応関係を算出することが可能であり、より精度良く色変換を行なえる色変換テーブルを作成することができる。そして、このように作成された色変換テーブルを用いることにより、色空間全体に亘る色変換精度をさらに向上させることができる（請求項１〜３）。
【０２２６】
さらに、本発明によれば、以下のような作用効果も得られる。
〔１〕変換先色空間として、人が見た際の色の差とその色空間での色の距離とがほぼ均等であるような均等色空間（ＣＩＥＬＡＢ，ＣＩＥＬＵＶ，ＣＩＥＣＡＭなど）を用いることで、色変換を行なう際に、色空間内のどのような色でも、色空間内の距離を重み係数として用いることができる。
【０２２７】
〔２〕変換先色空間を分割して得られる複数の領域が、相互に重複する領域を有することにより、色変換結果の階調性が保証される。つまり、重複領域に存在する色変換結果として、分割領域毎の色変換式により算出された複数の色変換結果に対して補間演算を行なった結果を用いることにより、各領域間の色の階調性が保証されることになる。
【０２２８】
〔３〕色相角に応じて変換先色空間の全体を分割して得られた複数の領域を用いることで、より精度良く対応関係を算出することが可能になる。つまり、色票の分光反射率を、色相角領域毎に分類された分光反射率群に基づいて主成分ベクトルを算出した場合、その主成分ベクトルの、各色相角領域に属する分光反射率群に対する寄与率は、当然、高くなる。従って、そのような主成分ベクトルは、各色相角領域に属する分光反射率群をより正確に表現することができ、結果として精度良く対応関係を算出できる。このような効果は、様々な色相角をもつ色をもつ色票を用いた場合に顕著に得られる。
【０２２９】
〔４〕一次色（例えばシアン，マゼンタ，イエローなど）の色相角を各領域の中心とすることで、より精度良く対応関係を算出することが可能になる。一次色中心付近に存在する色の分光反射率の形状は、一次色の分光反射率の形状に類似している。そのため、一次色中心付近に存在する分光反射率群に基づいて主成分ベクトルを算出した場合、その主成分ベクトルの、一次色中心付近における分光反射率群に対する寄与率は、当然、高くなり、精度良く対応関係を算出できる。
【０２３０】
〔５〕彩度に応じて変換先色空間の全体を分割して得られた複数の領域を用いることで、より精度良く対応関係を算出することが可能になる。つまり、色票の分光反射率を、彩度領域毎に分類された分光反射率群に基づいて主成分ベクトルを算出した場合、その主成分ベクトルの、各彩度領域に属する分光反射率群に対する寄与率は、当然、高くなる。従って、そのような主成分ベクトルは、各彩度領域に属する分光反射率群をより正確に表現することができ、結果として精度良く対応関係を算出できる。このような効果は、様々な彩度をもつ色をもつ色票を用いた場合に顕著に得られる。
【０２３１】
〔６〕明度に応じて変換先色空間の全体を分割して得られた複数の領域を用いることで、より精度良く対応関係を算出することが可能になる。つまり、色票の分光反射率を、明度領域毎に分類された分光反射率群に基づいて主成分ベクトルを算出した場合、その主成分ベクトルの、各明度領域に属する分光反射率群に対する寄与率は、当然、高くなる。従って、そのような主成分ベクトルは、各明度領域に属する分光反射率群をより正確に表現することができ、結果として精度良く対応関係を算出できる。このような効果は、様々な明度をもつ色をもつ色票を用いた場合に顕著に得られる。
【０２３２】
〔７〕各色変換式により変換元色空間の一の色信号を変換した色変換結果が各色変換式に対応した領域に属している場合に、その色変換結果は正しいものと判断し、正しいと判断された色変換結果に基づいて色変換テーブルを作成することで、より精度良く対応関係を算出することが可能になる。色空間の各分割領域から得られた色変換式による色変換精度は、当然ながら、自らの領域外の色信号に対して色変換を行なう際には、低くなる。そのため、精度の低い色変換結果と、精度の高い色変換結果とを判別することで、より高い精度の色変換結果を取得することができる。
【０２３３】
〔８〕一の色信号について正しいと判断された色変換結果が複数存在する場合、正しいと判断された複数の色変換結果と各色変換結果の属する領域の境界（境界面）との距離に関する値に基づいて、その一の色信号についての色変換結果を算出することで、色変換結果の階調性が保証される。各色変換式を用いて算出された複数の色変換結果が大きく異なり、それぞれの領域境界面に対してずれる場合があるが、このような場合、上述のような色変換結果の算出処理（内挿処理）を行なうことにより、階調性が保証されることになる。
【０２３４】
〔９〕一の色信号について正しいと判断された色変換結果が複数存在する場合、正しいと判断された複数の色変換結果と各色変換結果の属する領域の境界（境界面）との距離に関する値に基づいて、正しいと判断された複数の色変換結果のうちの一つを、その一の色信号についての色変換結果として選択することで、対応関係作成にかかる処理時間を短縮することができる。分割領域の数が多い場合などは、それに比例して、色変換結果の算出処理（内挿処理）にかかる時間が多くなる。そこで、処理時間を短縮するために、正しい色変換結果のうち、境界面からの距離が最も離れているものが、複数の正しい色変換結果の中でも最も信頼性の高いものであると判断して、その色変換結果を最終的な色変換結果とする。
【０２３５】
〔１０〕一の色信号について正しいと判断された色変換結果が存在しない場合、その一の色信号について得られた複数の色変換結果と各色変換結果の属する領域の境界（境界面）との距離に関する値の逆数に基づいて、その一の色信号についての色変換結果を算出することで、色変換結果が全て正しくないとしても、より信頼性のある色変換結果を算出することが可能になる。つまり、全ての色変換結果が正しくない場合は、正しくないものの中でも各境界面との距離が近いものの方が信頼性が高いものと判断している。
【０２３６】
〔１１〕一の色信号について正しいと判断された色変換結果が存在しない場合、その一の色信号について得られた複数の色変換結果と各色変換結果の属する領域の境界（境界面）との距離に関する値の逆数に基づいて、複数の色変換結果のうちの一つを、その一の色信号についての色変換結果として選択することで、色変換結果が全て正しくないとしても、より信頼性の高い色変換結果を短時間で得ることが可能になる。
【０２３７】
〔１２〕色変換テーブルを、変換元色空間の色信号を、色変換先空間の色信号としての、色変換結果に応じた分光反射率に対応させるものとすることにより、色変換結果を分光反射率にすることができ、色入力機器を簡易測色器の代用とすることが可能になる。
〔１３〕一の色信号について正しいと判断された色変換結果が複数存在する場合、正しいと判断された複数の色変換結果と各色変換結果の属する領域の境界（境界面）との距離に関する値に基づいて、その一の色信号についての分光反射率を算出することで、分光反射率の階調性が保証される。
【０２３８】
〔１４〕一の色信号について正しいと判断された色変換結果が複数存在する場合、正しいと判断された複数の色変換結果と各色変換結果の属する領域の境界（境界面）との距離に関する値に基づいて、正しいと判断された複数の色変換結果に応じた分光反射率のうちの一つを、その一の色信号についての分光反射率として選択することで、対応関係作成にかかる処理時間を短縮することができる。
【０２３９】
〔１５〕一の色信号について正しいと判断された色変換結果が存在しない場合、その一の色信号について得られた複数の色変換結果と各色変換結果の属する領域の境界（境界面）との距離に関する値の逆数に基づいて、その一の色信号についての分光反射率を算出することで、色変換結果が全て正しくないとしても、より信頼性のある分光反射率を算出することが可能になる。
【０２４０】
〔１６〕一の色信号について正しいと判断された色変換結果が存在しない場合、その一の色信号について得られた複数の色変換結果と各色変換結果の属する領域の境界（境界面）との距離に関する値の逆数に基づいて、複数の色変換結果に応じた分光反射率のうちの一つを、一の色信号についての分光反射率として選択することで、色変換結果が全て正しくないとしても、より信頼性の高い分光反射率を短時間で得ることが可能になる。
【０２４１】
〔１７〕複数の色変換結果と各色変換結果の属する領域の境界（境界面）との距離に関する値として、その距離を変換して得られるΔＥ９４の値を用いることにより、人間の目から見て、より正確な色変換を行なうことのできる色変換テーブルが作成される。ΔＥ９４は、L*a*b*空間上の距離（ΔＥａｂ）よりも人間の見た目に近い色差をもつため、このΔＥ９４を用いた色変換結果は、人間の目から見て、より正確なものになる。
【図面の簡単な説明】
【図１】本発明の一実施形態としての色変換テーブル作成装置の機能構成を示すブロック図である。
【図２】本発明の一実施形態としての色変換テーブル作成方法を実現するためのコンピュータシステム（色変換テーブル作成装置）のハードウェア構成を示すブロック図である。
【図３】ＲＧＢ空間とL*a*b*空間との対応関係を示す図である。
【図４】（Ａ）〜（Ｃ）は、いずれも、本実施形態における、色相角に応じた領域分割例を説明するための図である。
【図５】本発明の一実施形態としての色変換テーブル作成方法の流れを説明するためのフローチャートである。
【図６】本実施形態における分光特性推定手順を説明するためのフローチャートである。
【図７】本実施形態における色変換式作成手順および対応関係作成手順を説明するためのフローチャートである。
【図８】本実施形態における判断手順および内挿処理手順を説明するためのフローチャートである。
【図９】本実施形態における内挿処理の具体例を説明するための図である。
【図１０】本実施形態における内挿処理の具体例を説明するための図である。
【図１１】本実施形態における内挿処理の具体例を説明するための図である。
【図１２】本実施形態における内挿処理の具体例を説明するための図である。
【図１３】本実施形態の第１変形例における、彩度に応じた領域分割例を説明するための図である。
【図１４】本実施形態の第１変形例における内挿処理の具体例を説明するための図である。
【図１５】本実施形態の第１変形例における内挿処理の具体例を説明するための図である。
【図１６】本実施形態の第１変形例における内挿処理の具体例を説明するための図である。
【図１７】本実施形態の第２変形例における、明度に応じた領域分割例を説明するための図である。
【符号の説明】
１ コンピュータシステム（パーソナルコンピュータ，色変換テーブル作成装置）
１０ ＣＰＵ
１１ 分光特性推定部
１２ 分光反射率分類部
１４ 色変換テーブル作成部
１４１ 対応関係作成部
１４２ 判断部
１４３ 作成処理部
１５ ＲＧＢ空間格子点信号値出力部
２０ メモリ
２１ アプリケーションプログラム（色変換テーブル作成プログラム）
２２ 色変換テーブル
３０ ディスプレイ
４０ キーボード
４１ マウス
５０ バス
５１ 入出力インタフェース
５２ 入出力インタフェース（入力部）
６０ イメージスキャナ（色入力機器）
７０ 測色器
１００ 色票

Claims (3)

1. 色入力機器から出力された、該色入力機器の色空間（以下、変換元色空間）における色信号と、該変換元色空間とは異なる色空間（以下、変換先色空間という）における色信号とを対応させる色変換テーブルを作成する方法であって、
   該変換先色空間の全体を分割して得られる複数の領域を、予め定義しておく分割ステップと、
   色票上の複数種の色領域を該色入力機器で読み込み、各色領域に対応する、該変換元色空間での色信号を該色入力機器から出力する読込ステップと、
   該複数種の色領域を測色器で測色し、各色領域に対応する分光反射率を該測色器から出力する測色ステップと、
   該分光反射率に対応する、該変換先色空間での色信号が該色変換先色空間において該複数の領域のうちのどの領域に属するかに応じて、該分光反射率を分類する分類ステップと、
   該読込ステップで該色入力機器から出力された該色信号と該測色ステップで該測色器から出力された該分光反射率とに基づいて、該色入力機器の分光特性を推定する分光特性推定ステップと、
   該分割ステップで定義された該複数の領域にそれぞれ対応した複数の色変換式を用いて、該色変換テーブルを作成する色変換テーブル作成ステップとを有し、
   該色変換テーブル作成ステップが、
   該分類ステップで分類された該分光反射率と該分光特性推定ステップで推定された前記分光特性とに基づいて、該変換先色空間での前記領域毎に該色変換式を作成しながら、各領域に応じた該色変換式を用いて、該色変換式毎に、該変換元色空間における色信号と該変換先色空間における色信号との対応関係を複数作成する対応関係作成ステップと、
   該対応関係作成ステップにおいて該対応関係を複数作成する際に各色変換式により該変換元色空間の一の色信号を変換した色変換結果が、各色変換式に対応した該領域に属している場合に、該色変換結果は正しいものと判断する判断ステップと、
   該判断ステップで正しいと判断された該色変換結果に基づいて該色変換テーブルを作成する作成処理ステップとを有するとともに、
   該対応関係作成ステップにおいて、
   該分類ステップで分類された分光反射率群毎に当該分光反射率群に対して特異値分解を施すことにより算出された複数の主成分ベクトルと、該分光特性推定ステップで推定された前記分光特性とに基づいて、該変換元色空間における色信号値を代入することにより該複数の主成分ベクトルのそれぞれについての複数の重み係数を算出する第１式を、該色変換式として該変換先色空間での前記領域毎に作成し、
   該変換元色空間における色信号値毎に該第１式によって算出された該複数の重み係数をそれぞれ該複数の主成分ベクトルに乗算しその線形和を求めることによって当該色信号値に対応した分光反射率を前記領域毎に算出する第２式を、該色変換式として作成し、
   該第２式によって前記領域毎に算出された該分光反射率を該変換先色空間における色信号に変換する第３式を、該色変換式として作成し、
   これらの第１式，第２式および第３式を用いて該変換元色空間の一の色信号を該変換先色空間の色信号に変換することにより、当該一の色信号について、該分類ステップで分類された複数の分光反射率群のそれぞれに対応する複数の該対応関係を作成することを特徴とする、色変換テーブル作成方法。
2. 色入力機器から出力された、該色入力機器の色空間（以下、変換元色空間）における色信号と、該変換元色空間とは異なる色空間（以下、変換先色空間という）における色信号とを対応させる色変換テーブルを作成する装置であって、
   色票上の複数種の色領域を該色入力機器により読み込んで得られた、各色領域に対応する、該変換元色空間での色信号を入力する入力部と、
   該複数種の色領域を測色し、各色領域に対応する分光反射率を得る測色器と、
   該分光反射率に対応する、該変換先色空間での色信号が該色変換先色空間において該複 数の領域のうちのどの領域に属するかに応じて、該分光反射率を分類する分類部と、
   該入力部から入力された該色信号と該測色器により得られた該分光反射率とに基づいて、該色入力機器の分光特性を推定する分光特性推定部と、
   該変換先色空間を分割して得られる複数の領域にそれぞれ対応した複数の色変換式を用いて、該色変換テーブルを作成する色変換テーブル作成部をそなえ、
   該色変換テーブル作成部が、
   該分類部により分類された該分光反射率と該分光特性推定部により推定された前記分光特性とに基づいて、該変換先色空間での前記領域毎に該色変換式を作成しながら、各領域に応じた該色変換式を用いて、該色変換式毎に、該変換元色空間における色信号と該変換先色空間における色信号との対応関係を複数作成する対応関係作成部と、
   該対応関係作成部において該対応関係を複数作成する際に各色変換式により該変換元色空間の一の色信号を変換した色変換結果が、各色変換式に対応した該領域に属している場合に、該色変換結果は正しいものと判断する判断部と、
   該判断部により正しいと判断された該色変換結果に基づいて該色変換テーブルを作成する作成処理部とをそなえるとともに、
   該対応関係作成部が、
   該分類部により分類された分光反射率群毎に当該分光反射率群に対して特異値分解を施すことにより算出された複数の主成分ベクトルと、該分光特性推定部により推定された前記分光特性とに基づいて、該変換元色空間における色信号値を代入することにより該複数の主成分ベクトルのそれぞれについての複数の重み係数を算出する第１式を、該色変換式として該変換先色空間での前記領域毎に作成し、
   該変換元色空間における色信号値毎に該第１式によって算出された該複数の重み係数をそれぞれ該複数の主成分ベクトルに乗算しその線形和を求めることによって当該色信号値に対応した分光反射率を前記領域毎に算出する第２式を、該色変換式として作成し、
   該第２式によって前記領域毎に算出された該分光反射率を該変換先色空間における色信号に変換する第３式を、該色変換式として作成し、
   これらの第１式，第２式および第３式を用いて該変換元色空間の一の色信号を該変換先色空間の色信号に変換することにより、当該一の色信号について、該分類部により分類された複数の分光反射率群のそれぞれに対応する複数の該対応関係を作成することを特徴とする、色変換テーブル作成装置。
3. 色入力機器から出力された、該色入力機器の色空間（以下、変換元色空間）における色信号と、該変換元色空間とは異なる色空間（以下、変換先色空間という）における色信号とを対応させる色変換テーブルを、色票上の複数種の色領域を該色入力機器により読み込んで得られた、各色領域に対応する、該変換元色空間での色信号であって入力部から入力される色信号と、測色器により該複数種の色領域を測色して得られる各色領域に対応する分光反射率とに基づいて作成する機能をコンピュータにより実現させるための色変換テーブル作成プログラムを記録したコンピュータ読取可能な記録媒体であって、
   該色変換テーブル作成プログラムが、
   該分光反射率に対応する、該変換先色空間での色信号が該色変換先色空間において該複数の領域のうちのどの領域に属するかに応じて、該分光反射率を分類する分類部、
   該入力部から入力された該色信号と該測色器により得られた該分光反射率とに基づいて、該色入力機器の分光特性を推定する分光特性推定部、および、
   該変換先色空間を分割して得られる複数の領域にそれぞれ対応した複数の色変換式を用いて、該色変換テーブルを作成する色変換テーブル作成部として、該コンピュータを機能させ、
   該色変換テーブル作成プログラムが該コンピュータを該色変換テーブル作成部として機能させる際に、
   該色変換テーブル作成プログラムが、
   該分類部により分類された該分光反射率と該分光特性推定部により推定された前記分光 特性とに基づいて、該変換先色空間での前記領域毎に該色変換式を作成しながら、各領域に応じた該色変換式を用いて、該色変換式毎に、該変換元色空間における色信号と該変換先色空間における色信号との対応関係を複数作成する対応関係作成部、
   該対応関係作成部において該対応関係を複数作成する際に各色変換式により該変換元色空間の一の色信号を変換した色変換結果が、各色変換式に対応した該領域に属している場合に、該色変換結果は正しいものと判断する判断部、および、
   該判断部により正しいと判断された該色変換結果に基づいて該色変換テーブルを作成する作成処理部として、該コンピュータをさらに機能させるとともに、
   該色変換テーブル作成プログラムが該コンピュータを該対応関係作成部として機能させる際に、
   色変換テーブル作成プログラムが、該コンピュータに、
   該分類部により分類された分光反射率群毎に当該分光反射率群に対して特異値分解を施すことにより算出された複数の主成分ベクトルと、該分光特性推定部により推定された前記分光特性とに基づいて、該変換元色空間における色信号値を代入することにより該複数の主成分ベクトルのそれぞれについての複数の重み係数を算出する第１式を、該色変換式として該変換先色空間での前記領域毎に作成させ、
   該変換元色空間における色信号値毎に該第１式によって算出された該複数の重み係数をそれぞれ該複数の主成分ベクトルに乗算しその線形和を求めることによって当該色信号値に対応した分光反射率を前記領域毎に算出する第２式を、該色変換式として作成させ、
   該第２式によって前記領域毎に算出された該分光反射率を該変換先色空間における色信号に変換する第３式を、該色変換式として作成させ、
   これらの第１式，第２式および第３式を用いて該変換元色空間の一の色信号を該変換先色空間の色信号に変換することにより、当該一の色信号について、該分類部により分類された複数の分光反射率群のそれぞれに対応する複数の該対応関係を作成させることを特徴とする、色変換テーブル作成プログラムを記録したコンピュータ読取可能な記録媒体。

Images