知財求人 - 知財ポータルサイト「IP Force」
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2024167946
(43)【公開日】2024-12-05
(54)【発明の名称】画像処理装置、画像処理方法及びプログラム
(51)【国際特許分類】
H04N 1/60 20060101AFI20241128BHJP
H04N 1/48 20060101ALI20241128BHJP
G06T 1/00 20060101ALI20241128BHJP
H04N 9/64 20230101ALI20241128BHJP
【ＦＩ】
H04N1/60
H04N1/48
G06T1/00 510
H04N9/64 A
【審査請求】未請求
【請求項の数】15
【出願形態】ＯＬ
(21)【出願番号】P 2023084285
(22)【出願日】2023-05-23
(71)【出願人】
【識別番号】000001007
【氏名又は名称】キヤノン株式会社
(74)【代理人】
【識別番号】100090273
【弁理士】
【氏名又は名称】國分 孝悦
(72)【発明者】
【氏名】小倉 暢
【テーマコード（参考）】
5B057
5C066
【Ｆターム（参考）】
5B057CA01
5B057CA08
5B057CA12
5B057CA16
5B057CB01
5B057CB08
5B057CB16
5B057CE17
5B057CE18
5B057DA17
5B057DB06
5B057DB09
5B057DC25
5C066AA01
5C066BA20
5C066CA08
(57)【要約】
【課題】画像を色変換する際に生じる誤差を低減させる。
【解決手段】デジタルカメラ１で複数のパッチ（色）をもつ色票３を撮像した画像（ソース色空間の画像）と、特定のカメラで色票３を撮像した画像（ターゲット色空間の画像）のそれぞれから計測された、パッチ（色）の計測値のデータを取得する。そして、色空間上の入力点を取得する。そして、ソース色空間の画像から計測された計測値のデータから、複数のパッチ（色）に対応する複数の基準点候補を取得し、複数の基準点候補の中から、入力点との色空間上の距離が近い複数の基準点候補を、基準点として選択する。そして、入力点の位置に基づいて、色変換の係数を取得し、取得した色変換の係数と、複数の基準点に対応する、ターゲット色空間の画像から計測された計測値のデータに基づいて、入力点における色変換結果を算出する。
【選択図】図４
【解決手段】デジタルカメラ１で複数のパッチ（色）をもつ色票３を撮像した画像（ソース色空間の画像）と、特定のカメラで色票３を撮像した画像（ターゲット色空間の画像）のそれぞれから計測された、パッチ（色）の計測値のデータを取得する。そして、色空間上の入力点を取得する。そして、ソース色空間の画像から計測された計測値のデータから、複数のパッチ（色）に対応する複数の基準点候補を取得し、複数の基準点候補の中から、入力点との色空間上の距離が近い複数の基準点候補を、基準点として選択する。そして、入力点の位置に基づいて、色変換の係数を取得し、取得した色変換の係数と、複数の基準点に対応する、ターゲット色空間の画像から計測された計測値のデータに基づいて、入力点における色変換結果を算出する。
【選択図】図４
【特許請求の範囲】
【請求項1】
複数の色票を撮像した画像において計測された、前記複数の色票の色を表す第１のデータと、前記第１のデータが計測された計測点における目標の色を表す第２のデータを取得するデータ取得手段と、
色空間上の入力点を取得する入力点取得手段と、
前記第１のデータから、前記複数の色票に対応する複数の基準点候補を取得し、複数の前記基準点候補の中から、前記入力点との色空間上の距離が近い複数の前記基準点候補を、基準点として選択する基準点選択手段と、
前記入力点の位置に基づいて、色変換の係数を取得する係数取得手段と、
前記色変換の係数と、複数の前記基準点に対応する前記第２のデータに基づいて、前記入力点における色変換結果を算出する算出手段と、
を有することを特徴とする画像処理装置。
色空間上の入力点を取得する入力点取得手段と、
前記第１のデータから、前記複数の色票に対応する複数の基準点候補を取得し、複数の前記基準点候補の中から、前記入力点との色空間上の距離が近い複数の前記基準点候補を、基準点として選択する基準点選択手段と、
前記入力点の位置に基づいて、色変換の係数を取得する係数取得手段と、
前記色変換の係数と、複数の前記基準点に対応する前記第２のデータに基づいて、前記入力点における色変換結果を算出する算出手段と、
を有することを特徴とする画像処理装置。
【請求項2】
前記係数取得手段は、複数の前記基準点と、前記入力点との相対的な位置に基づいて、前記色変換の係数を算出することを特徴とする請求項１に記載の画像処理装置。
【請求項3】
前記算出手段は、前記計測点における、前記第１のデータと前記第２のデータとの変分値を算出し、前記色変換の係数と、複数の前記基準点に対応する前記変分値に基づいて、前記入力点における色変換結果を算出することを特徴とする請求項１に記載の画像処理装置。
【請求項4】
前記画像は、撮像装置において撮像された画像であって、
前記撮像装置が３次元ルックアップテーブルで保持する色変換関数を生成する生成手段を更に有し、
前記生成手段は、前記３次元ルックアップテーブルの格子点に対応する前記入力点と、当該入力点における前記色変換結果とに基づいて、前記色変換関数を生成することを特徴とする請求項１に記載の画像処理装置。
前記撮像装置が３次元ルックアップテーブルで保持する色変換関数を生成する生成手段を更に有し、
前記生成手段は、前記３次元ルックアップテーブルの格子点に対応する前記入力点と、当該入力点における前記色変換結果とに基づいて、前記色変換関数を生成することを特徴とする請求項１に記載の画像処理装置。
【請求項5】
前記データ取得手段は、同一の色票から輝度の異なる複数の計測値を取得し、
前記複数の計測値の間の値を前記複数の色票ごとに補間する補間手段、を更に有し、
前記基準点選択手段は、前記補間手段により補間された前記第１のデータから、前記複数の色票に対応する複数の基準点候補を取得することを特徴とする請求項１に記載の画像処理装置。
前記複数の計測値の間の値を前記複数の色票ごとに補間する補間手段、を更に有し、
前記基準点選択手段は、前記補間手段により補間された前記第１のデータから、前記複数の色票に対応する複数の基準点候補を取得することを特徴とする請求項１に記載の画像処理装置。
【請求項6】
前記補間手段は、前記複数の計測値を輝度で表現する関数を前記複数の色票ごとに生成することを特徴とする請求項５に記載の画像処理装置。
【請求項7】
前記補間手段は、前記複数の計測値を多項式近似することにより前記関数を生成することを特徴とする請求項６に記載の画像処理装置。
【請求項8】
前記画像は、露光時間を多段階にして前記複数の色票を撮像した画像データ群、開口率を多段階にして前記複数の色票を撮像した画像データ群、及び照明強度を多段階にして前記複数の色票を撮像した画像データ群のうちの少なくとも何れか１つを含むことを特徴とする請求項５に記載の画像処理装置。
【請求項9】
前記基準点選択手段は、前記複数の色票ごとの前記関数に対して、前記入力点の輝度値を代入して得られる点を、複数の前記基準点候補として取得することを特徴とする請求項６に記載の画像処理装置。
【請求項10】
前記基準点選択手段は、複数の前記基準点候補の中から、前記入力点との色相の近い複数の前記基準点候補を、基準点として選択することを特徴とする請求項１に記載の画像処理装置。
【請求項11】
前記データ取得手段は、同一の色票から輝度の異なる複数の計測値を取得し、
前記複数の計測値の間の値を前記複数の色票ごとに補間する補間手段と、
前記補間手段により補間された前記第１のデータを表す形状データを、前記複数の色票ごとに表示するよう制御する表示制御手段と、
を更に有することを特徴する請求項１に記載の画像処理装置。
前記複数の計測値の間の値を前記複数の色票ごとに補間する補間手段と、
前記補間手段により補間された前記第１のデータを表す形状データを、前記複数の色票ごとに表示するよう制御する表示制御手段と、
を更に有することを特徴する請求項１に記載の画像処理装置。
【請求項12】
前記基準点選択手段は、前記複数の色票に対応する前記第１のデータのうち、所定の条件を満たす一部の色票に対応する前記第１のデータを、複数の基準点候補を取得する対象から除外することを特徴する請求項１に記載の画像処理装置。
【請求項13】
前記基準点選択手段は、前記複数の色票に対応する前記第１のデータのうち、ユーザより操作指示のあった色票に対応する前記第１のデータを、複数の基準点候補を取得する対象から除外することを特徴とする請求項１２に記載の画像処理装置。
【請求項14】
複数の色票を撮像した画像において計測された、前記複数の色票の色を表す第１のデータと、前記第１のデータが計測された計測点における目標の色を表す第２のデータを取得するデータ取得工程と、
色空間上の入力点を取得する入力点取得工程と、
前記第１のデータから、前記複数の色票に対応する複数の基準点候補を取得し、複数の前記基準点候補の中から、前記入力点との色空間上の距離が近い複数の前記基準点候補を、基準点として選択する基準点選択工程と、
前記入力点の位置に基づいて、色変換の係数を取得する係数取得工程と、
前記色変換の係数と、複数の前記基準点に対応する前記第２のデータに基づいて、前記入力点における色変換結果を算出する算出工程と、
を含むことを特徴とする画像処理方法。
色空間上の入力点を取得する入力点取得工程と、
前記第１のデータから、前記複数の色票に対応する複数の基準点候補を取得し、複数の前記基準点候補の中から、前記入力点との色空間上の距離が近い複数の前記基準点候補を、基準点として選択する基準点選択工程と、
前記入力点の位置に基づいて、色変換の係数を取得する係数取得工程と、
前記色変換の係数と、複数の前記基準点に対応する前記第２のデータに基づいて、前記入力点における色変換結果を算出する算出工程と、
を含むことを特徴とする画像処理方法。
【請求項15】
画像処理装置のコンピュータを、
複数の色票を撮像した画像において計測された、前記複数の色票の色を表す第１のデータと、前記第１のデータが計測された計測点における目標の色を表す第２のデータを取得するデータ取得手段と、
色空間上の入力点を取得する入力点取得手段と、
前記第１のデータから、前記複数の色票に対応する複数の基準点候補を取得し、複数の前記基準点候補の中から、前記入力点との色空間上の距離が近い複数の前記基準点候補を、基準点として選択する基準点選択手段と、
前記入力点の位置に基づいて、色変換の係数を取得する係数取得手段と、
前記色変換の係数と、複数の前記基準点に対応する前記第２のデータに基づいて、前記入力点における色変換結果を算出する算出手段、
として機能させるためのプログラム。
複数の色票を撮像した画像において計測された、前記複数の色票の色を表す第１のデータと、前記第１のデータが計測された計測点における目標の色を表す第２のデータを取得するデータ取得手段と、
色空間上の入力点を取得する入力点取得手段と、
前記第１のデータから、前記複数の色票に対応する複数の基準点候補を取得し、複数の前記基準点候補の中から、前記入力点との色空間上の距離が近い複数の前記基準点候補を、基準点として選択する基準点選択手段と、
前記入力点の位置に基づいて、色変換の係数を取得する係数取得手段と、
前記色変換の係数と、複数の前記基準点に対応する前記第２のデータに基づいて、前記入力点における色変換結果を算出する算出手段、
として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、画像処理装置、画像処理方法及びプログラムに関する。
【背景技術】
【0002】
デジタルカメラは、予め用意された複数の色処理パラメータを内部に保持し、例えばユーザが選択したモードに応じた色処理パラメータを使用して、撮像画像に対する色変換などの画像処理を行う機能を有する。特許文献１には、カラーチャートをデジタルカメラで撮像したときの撮像データを、カラーチャートの本来の色、または、ユーザが好ましいと思う色に近づけるような、色変換関数を求める技術が開示されている。この色変換関数は、３次元ルックアップテーブル（３ＤＬＵＴ）という形式で保持されることもある。
【先行技術文献】
【特許文献】
【0003】
【特許文献1】特開２００４－１５３６８４号公報
【発明の概要】
【発明が解決しようとする課題】
【0004】
しかしながら、特許文献１では、上記の色変換関数は、入力データのＲＧＢ値を説明変数とした３次の項までの和で表現されている。そのため、僅かな入力色の変化に対して出力色が大きく変化するような性質（カオス性）のある条件下においては、関数の自由度が不足し、色変換関数による変換誤差が大きくなる。このような変換誤差の影響は、デジタルカメラの色再現範囲が広くなるほど顕著になり、例えば最近のＨＤＲ（ハイダイナミックレンジ）再現が可能なデジタルカメラでは無視できないものになっている。
【0005】
そこで、本発明は、画像を色変換する際に生じる誤差を低減させることを目的とする。
【課題を解決するための手段】
【0006】
本発明の画像処理装置は、複数の色票を撮像した画像において計測された、前記複数の色票の色を表す第１のデータと、前記第１のデータが計測された計測点における目標の色を表す第２のデータを取得するデータ取得手段と、色空間上の入力点を取得する入力点取得手段と、前記第１のデータから、前記複数の色票に対応する複数の基準点候補を取得し、複数の前記基準点候補の中から、前記入力点との色空間上の距離が近い複数の前記基準点候補を、基準点として選択する基準点選択手段と、前記入力点の位置に基づいて、色変換の係数を取得する係数取得手段と、前記色変換の係数と、複数の前記基準点に対応する前記第２のデータに基づいて、前記入力点における色変換結果を算出する算出手段と、を有することを特徴とする。
【発明の効果】
【0007】
本発明によれば、画像を色変換する際に生じる誤差を低減させることができる。
【図面の簡単な説明】
【0008】
【図1A】画像処理装置を含むシステムの構成例を示す図である。
【図1B】画像処理装置のハードウェア構成例を示す図である。
【図2】色票の一例を示す図である。
【図3】画像処理部の機能構成例を示す図である。
【図4】実施形態１に係る画像処理装置による処理を示すフローチャートである。
【図5】色票の計測値をプロットした例を示す図である。
【図6】計測値間の補間処理を説明する図である。
【図7】基準点の選択処理を説明する図である。
【図8】変換係数の算出処理を説明する図である。
【図9】学習データの出力例を示す図である。
【図10】実施形態２に係る画像処理装置による処理を示すフローチャートである。
【図11】ＧＵＩの表示例を示す図である。
【発明を実施するための形態】
【0009】
以下、本発明の実施形態について、図面を参照して説明する。尚、以下の実施形態は本発明を必ずしも限定するものではない。
【0010】
［実施形態１］
図１Ａは、本実施形態に係る画像処理装置２が適用されるシステムの構成の例を示すブロック図である。図１に示すシステムは、撮像装置の一例であるデジタルカメラ１と、色変換結果を算出する画像処理装置２を含む。このシステムでは、デジタルカメラ１で色票３を撮像した画像（以下、ソース色空間の画像と呼ぶことがある）を、理想の発色をもつ特定のカメラで色票３を撮像した画像（以下、ターゲット色空間の画像と呼ぶことがある）へ色変換するための３ＤＬＵＴを生成する。３ＤＬＵＴは、３次元ルックアップテーブルのことである。図２は、色票３の一例を示す図である。色票３は、例えば肌色、空の色、草の色、原色など色彩値の異なる複数のパッチが印刷等されたカラーチャートである。色票３として、ｘ－ｒｉｔｅ社のＣｏｌｏｒ Ｃｈｅｃｋｅｒ ＳＧ（１４０色）やＣｏｌｏｒ Ｃｈｅｃｋｅｒ（２４色）など、色設計に使用可能なカラーチャートが使用される。
図１Ａは、本実施形態に係る画像処理装置２が適用されるシステムの構成の例を示すブロック図である。図１に示すシステムは、撮像装置の一例であるデジタルカメラ１と、色変換結果を算出する画像処理装置２を含む。このシステムでは、デジタルカメラ１で色票３を撮像した画像（以下、ソース色空間の画像と呼ぶことがある）を、理想の発色をもつ特定のカメラで色票３を撮像した画像（以下、ターゲット色空間の画像と呼ぶことがある）へ色変換するための３ＤＬＵＴを生成する。３ＤＬＵＴは、３次元ルックアップテーブルのことである。図２は、色票３の一例を示す図である。色票３は、例えば肌色、空の色、草の色、原色など色彩値の異なる複数のパッチが印刷等されたカラーチャートである。色票３として、ｘ－ｒｉｔｅ社のＣｏｌｏｒ Ｃｈｅｃｋｅｒ ＳＧ（１４０色）やＣｏｌｏｒ Ｃｈｅｃｋｅｒ（２４色）など、色設計に使用可能なカラーチャートが使用される。
【0011】
まずデジタルカメラ１について説明する。
デジタルカメラ１は、撮像部１０１、データ入出力部１０２、画像処理部１０３、保持部１０４、および表示部１０５などを有して構成されている。
撮像部１０１は、撮像レンズとＣＣＤなどの画像センサ等から構成されている。データ入出力部１０２は、画像データおよび色変換関数などのデータを、外部機器との間で送受信可能なインタフェイスである。保持部１０４は、撮像部１０１による撮像画像、画像処理部１０３による画像処理後の画像データ、色変換関数などのデータを保持する。画像処理部１０３は、撮像部１０１による撮像画像に対して色処理等の画像処理を施す。表示部１０５は、液晶ディスプレイなどの表示装置を含み、撮像中や撮像後、色処理後の画像などを表示する。画像処理部１０３は、保持部１０４で保持する色変換関数を用いて画像処理を行う。この色変換関数は、３ＤＬＵＴの形態で保持されているものとする。また、データ入出力部１０２には、ケーブルなどを介して画像処理装置２が接続され、またメモリカードなどの記録媒体を接続することができる。
デジタルカメラ１は、撮像部１０１、データ入出力部１０２、画像処理部１０３、保持部１０４、および表示部１０５などを有して構成されている。
撮像部１０１は、撮像レンズとＣＣＤなどの画像センサ等から構成されている。データ入出力部１０２は、画像データおよび色変換関数などのデータを、外部機器との間で送受信可能なインタフェイスである。保持部１０４は、撮像部１０１による撮像画像、画像処理部１０３による画像処理後の画像データ、色変換関数などのデータを保持する。画像処理部１０３は、撮像部１０１による撮像画像に対して色処理等の画像処理を施す。表示部１０５は、液晶ディスプレイなどの表示装置を含み、撮像中や撮像後、色処理後の画像などを表示する。画像処理部１０３は、保持部１０４で保持する色変換関数を用いて画像処理を行う。この色変換関数は、３ＤＬＵＴの形態で保持されているものとする。また、データ入出力部１０２には、ケーブルなどを介して画像処理装置２が接続され、またメモリカードなどの記録媒体を接続することができる。
【0012】
被写体等の撮像が行われる場合、デジタルカメラ１では以下のような処理が行われる。
ユーザによって電源スイッチ（不図示）がオンにされると、撮像部１０１は、継続的に画像を撮像する動作を開始する。画像処理部１０３は、保持部１０４に保持されている３ＤＬＵＴを使用して、撮像部１０１から継続的に入力される画像データに画像処理を施し、さらにその画像処理後の画像データから表示用画像を生成する。表示部１０５は、画像処理部１０３から継続的に送られてくる表示用画像を表示する。ユーザは、表示部１０５に表示された表示用画像を観察して画面の構図を決め、シャッタボタン（不図示）を押す。シャッタボタンが押されると、画像処理部１０３は、シャッタボタンが押されたタイミングで撮像部１０１が取得した画像データに対して、保持部１０４の３ＤＬＵＴを用いた画像処理を施す。３ＤＬＵＴを用いた画像処理は、色、明るさを調整する処理である。画像処理が施された後の画像データは、データ入出力部１０２を介して、外部機器（画像処理装置２や記録媒体など）に出力される。なお、撮像部１０１が取得した画像データを、そのままデータ入出力部１０２へ出力するような色変換関数が保持部１０４に設定されている場合には、撮像部１０１が取得した画像データそのものがデジタルカメラ１から出力される。
ユーザによって電源スイッチ（不図示）がオンにされると、撮像部１０１は、継続的に画像を撮像する動作を開始する。画像処理部１０３は、保持部１０４に保持されている３ＤＬＵＴを使用して、撮像部１０１から継続的に入力される画像データに画像処理を施し、さらにその画像処理後の画像データから表示用画像を生成する。表示部１０５は、画像処理部１０３から継続的に送られてくる表示用画像を表示する。ユーザは、表示部１０５に表示された表示用画像を観察して画面の構図を決め、シャッタボタン（不図示）を押す。シャッタボタンが押されると、画像処理部１０３は、シャッタボタンが押されたタイミングで撮像部１０１が取得した画像データに対して、保持部１０４の３ＤＬＵＴを用いた画像処理を施す。３ＤＬＵＴを用いた画像処理は、色、明るさを調整する処理である。画像処理が施された後の画像データは、データ入出力部１０２を介して、外部機器（画像処理装置２や記録媒体など）に出力される。なお、撮像部１０１が取得した画像データを、そのままデータ入出力部１０２へ出力するような色変換関数が保持部１０４に設定されている場合には、撮像部１０１が取得した画像データそのものがデジタルカメラ１から出力される。
【0013】
次に画像処理装置２の機能構成について説明する。
画像処理装置２は、ソース色空間の画像をターゲット色空間の画像へ変換するための色変換関数を生成する。画像処理装置２は、データ入出力部２０１、画像処理部２０２、平均値算出部２０３、表示制御部２０４、ユーザ操作取得部２０５、色処理パラメータ保持部２０６、データ保持部２０７、および３ＤＬＵＴ生成部２０８などの機能を有している。
画像処理装置２は、ソース色空間の画像をターゲット色空間の画像へ変換するための色変換関数を生成する。画像処理装置２は、データ入出力部２０１、画像処理部２０２、平均値算出部２０３、表示制御部２０４、ユーザ操作取得部２０５、色処理パラメータ保持部２０６、データ保持部２０７、および３ＤＬＵＴ生成部２０８などの機能を有している。
【0014】
データ入出力部２０１は、画像データ、色変換関数などの各種のデータを外部機器との間でやり取りする。データ入出力部２０１は、汎用Ｉ／Ｆ２１５（図１Ｂ）や通信Ｉ／Ｆ２１７（図１Ｂ）を介してデジタルカメラ１などとの間で各種のデータの送受信を行ってよく、不図示のメモリカードなどの記録媒体との間で各種のデータの入出力を行ってもよい。ここでデータ入出力部２０１から入力される画像データ（入力画像データ）は、デジタルカメラ１や特定のカメラで色票３を撮像した画像とする。ここで、デジタルカメラ１で色票３を撮像した画像（ソース色空間の画像）は、撮像部１０１が取得した画像データそのものある。また、特定のカメラで色票３を撮像した画像（ターゲット色空間の画像）は、ソース色空間の画像の目標とする色で構成された画像データの例である。なお、ターゲット色空間の画像は、カメラで撮像した画像でなくても、ＣＧにより得られたデータなど、ユーザが理想とする色彩値が直接的に与えられたデータであっても構わない。
【0015】
画像処理部２０２は、データ入出力部２０１からの入力画像データを用いて、色空間上の任意の入力点に対する、色変換結果を算出する。画像処理部２０２の機能の詳細については、図３で後述する。
平均値算出部２０３は、データ入出力部２０１からの入力画像データから得られる色彩値の平均値を、計測値として算出する。ここでは色票３を撮像した画像内のパッチごとに、Ｌ*ａ*ｂ*値の平均値を、計測値として算出する。以下、ソース色空間の画像とターゲット色空間の画像のそれぞれから得られた計測値を学習データと呼ぶ。本実施形態では、ソース色空間の画像とターゲット色空間の画像は、多段階の露光時間で色票３を撮像した画像データ群である。ここでパッチ数を１８、露光時間の段階数を３２とすると、ソース色空間の画像とターゲット色空間の画像のそれぞれから、１８×３２の計測点における計測値が得られる。
平均値算出部２０３は、データ入出力部２０１からの入力画像データから得られる色彩値の平均値を、計測値として算出する。ここでは色票３を撮像した画像内のパッチごとに、Ｌ*ａ*ｂ*値の平均値を、計測値として算出する。以下、ソース色空間の画像とターゲット色空間の画像のそれぞれから得られた計測値を学習データと呼ぶ。本実施形態では、ソース色空間の画像とターゲット色空間の画像は、多段階の露光時間で色票３を撮像した画像データ群である。ここでパッチ数を１８、露光時間の段階数を３２とすると、ソース色空間の画像とターゲット色空間の画像のそれぞれから、１８×３２の計測点における計測値が得られる。
【0016】
データ保持部２０７は、入力画像データ、画像処理後の画像データ、学習データなどのデータを保持する。
色処理パラメータ保持部２０６は、デジタルカメラ１における色変換関数を保持する。
３Ｄルックアップテーブル（ＬＵＴ）生成部２０８は、画像処理部２０２により算出された画像処理結果を用いて３ＤＬＵＴを生成する。
ユーザ操作取得部２０５は、マウス、キーボード、タッチパネルなどの入力装置２２０（図１Ｂ）からユーザの入力操作を取得する。
表示制御部２０４は、３ＤＬＵＴの後述するグリッド数（格子点数）等を入力する際のＧＵＩや、学習データに対する操作指示を受け付けるためのＧＵＩを、表示装置２１４（図１Ｂ）表示する。
色処理パラメータ保持部２０６は、デジタルカメラ１における色変換関数を保持する。
３Ｄルックアップテーブル（ＬＵＴ）生成部２０８は、画像処理部２０２により算出された画像処理結果を用いて３ＤＬＵＴを生成する。
ユーザ操作取得部２０５は、マウス、キーボード、タッチパネルなどの入力装置２２０（図１Ｂ）からユーザの入力操作を取得する。
表示制御部２０４は、３ＤＬＵＴの後述するグリッド数（格子点数）等を入力する際のＧＵＩや、学習データに対する操作指示を受け付けるためのＧＵＩを、表示装置２１４（図１Ｂ）表示する。
【0017】
続いて、画像処理部２０２の内部構成について説明する。
図３は、画像処理部２０２の内部構成の例を示すブロック図である。画像処理部２０２は、計測値取得部３０１、計測値補間部３０２、入力点取得部３０３、基準点選択部３０４、変換係数算出部３０５、変分値取得部３０６、および色変換結果算出部３０７などの機能を有している。
計測値取得部３０１は、平均値算出部２０３から計測値を取得し、ソース色空間の画像とターゲット色空間の画像について対応する各計測点における計測値間のＬ*ａ*ｂ*の成分ごとの変分値を算出する。
計測値補間部３０２は、各計測点における変分値の間を補間する。具体的には、計測値補間部３０２は、色票３のパッチ（色）ごとに、変分値のＬ*ａ*ｂ*の各成分を輝度ｌの関数で表現する式を求める。
図３は、画像処理部２０２の内部構成の例を示すブロック図である。画像処理部２０２は、計測値取得部３０１、計測値補間部３０２、入力点取得部３０３、基準点選択部３０４、変換係数算出部３０５、変分値取得部３０６、および色変換結果算出部３０７などの機能を有している。
計測値取得部３０１は、平均値算出部２０３から計測値を取得し、ソース色空間の画像とターゲット色空間の画像について対応する各計測点における計測値間のＬ*ａ*ｂ*の成分ごとの変分値を算出する。
計測値補間部３０２は、各計測点における変分値の間を補間する。具体的には、計測値補間部３０２は、色票３のパッチ（色）ごとに、変分値のＬ*ａ*ｂ*の各成分を輝度ｌの関数で表現する式を求める。
【0018】
入力点取得部３０３は、ソース色空間における入力点Ｐを取得する。
基準点選択部３０４は、ソース色空間の画像から得られた計測値に基づいて、入力点Ｐと色相の近い２つ以上の基準点を選択する。
変換係数算出部３０５は、入力点Ｐの位置を、基準点選択部３０４で選択した基準点の位置で表現するような係数（以下、変換係数という）を算出する。
変分値取得部３０６は、変換係数算出部３０５で算出した変換係数と、基準点選択部３０４で選択した基準点に対応する変分値を用いて、入力点Ｐにおける変分の推定値を算出する。
色変換結果算出部３０７は、入力点Ｐと変分値取得部３０６で算出した変分の推定値に基づいて、入力点Ｐにおける色変換結果を算出する。
基準点選択部３０４は、ソース色空間の画像から得られた計測値に基づいて、入力点Ｐと色相の近い２つ以上の基準点を選択する。
変換係数算出部３０５は、入力点Ｐの位置を、基準点選択部３０４で選択した基準点の位置で表現するような係数（以下、変換係数という）を算出する。
変分値取得部３０６は、変換係数算出部３０５で算出した変換係数と、基準点選択部３０４で選択した基準点に対応する変分値を用いて、入力点Ｐにおける変分の推定値を算出する。
色変換結果算出部３０７は、入力点Ｐと変分値取得部３０６で算出した変分の推定値に基づいて、入力点Ｐにおける色変換結果を算出する。
【0019】
図１Ｂは、画像処理装置２のハードウェア構成を示すブロック図である。画像処理装置２は、ＣＰＵ２１１、ＲＯＭ２１２、ＲＡＭ２１３、表示装置２１４、汎用Ｉ／Ｆ２１５、ＳＡＴＡ（シリアルＡＴＡ）Ｉ／Ｆ２１６、および通信Ｉ／Ｆ２１７などを有して構成されている。汎用Ｉ／Ｆ（インタフェイス）２１５には、シリアルバス２１９を介して、マウス、キーボード、タッチパネルなどの入力装置２２０やデジタルカメラ１が接続されている。ＳＡＴＡＩ／Ｆ（インタフェイス）２１６には、シリアルバス２２１を介して、ＨＤＤ（ハードディスクドライブ）２２２やメモリカードなどの記録媒体の読み書きを行う汎用ドライブ２２３が接続されている。
【0020】
ＣＰＵ２１１は、ＲＡＭ２１３をワークメモリとして、ＲＯＭ２１２やＨＤＤ２２２などに格納されたＯＳ（オペレーティングシステム）や各種のプログラムを実行する。また、ＣＰＵ２１１は、システムバス２１８を介して各構成を制御する。なお、ＲＯＭ２１２やＨＤＤ２２２などに格納されたプログラムがＲＡＭ２１３に展開され、ＣＰＵ２１１によって実行されることで、後述するフローチャートによる処理が実現する。また、ＲＯＭ２１２やＨＤＤ２２２などに格納されたプログラムがＲＡＭ２１３に展開され、ＣＰＵ２１１によって実行されることで、図１や図３に示す画像処理装置２の各種の機能が実現する。通信Ｉ／Ｆ２１７は、外部装置との間で情報の送受信を行う。ＣＰＵ２１１は、プログラムによって提供される各種のＧＵＩを表示装置２１４に表示し、入力装置２２０を介してユーザの入力操作を取得し、その入力操作に応じて、ＧＵＩの表示を制御する。表示装置２１４と入力装置２２０は、表示機能と入力機能を兼ね備えたタッチパネルディスプレイにより構成されても構わない。ＣＰＵ２１１は、ＨＤＤ２２２や汎用ドライブ２２３にマウントされた記録媒体を、入力画像データや色変換関数などの各種のデータの格納場所として使用する。
【0021】
図４は、本実施形態により実行される画像処理装置２の処理の流れを示すフローチャートである。以下の説明では、各工程（ステップ）について先頭にＳを付けて表記することで、工程（ステップ）の表記を省略する。
Ｓ４０１において、データ入出力部２０１は、デジタルカメラ１で色票３を撮像した画像データ群と、特定のカメラで色票３を撮像した画像データ群を入力する。これらの画像データ群は、ソース色空間の画像とターゲット色空間の画像のペア画像であり、かつ、均一照明下で撮像媒体の露光量を調整しながら複数回撮像されたデータセットであるとする。ここでは、様々な輝度レベルの画像データが得られるよう、多段階（ここでは、３２段階）の露光時間で色票３が撮像されている。このとき、撮像結果が色空間の低輝度側から高輝度側までを被覆するように露光時間が調整されている。
Ｓ４０１において、データ入出力部２０１は、デジタルカメラ１で色票３を撮像した画像データ群と、特定のカメラで色票３を撮像した画像データ群を入力する。これらの画像データ群は、ソース色空間の画像とターゲット色空間の画像のペア画像であり、かつ、均一照明下で撮像媒体の露光量を調整しながら複数回撮像されたデータセットであるとする。ここでは、様々な輝度レベルの画像データが得られるよう、多段階（ここでは、３２段階）の露光時間で色票３が撮像されている。このとき、撮像結果が色空間の低輝度側から高輝度側までを被覆するように露光時間が調整されている。
【0022】
なお、露光時間の段階数は３２段階に限られない。また、撮像媒体の露光量を調整する方法は、露光時間の調整に限定されず、絞り（開口率）の調整、照明強度の調整、別色票の使用などの方法で行っても構わない。また、低輝度から高輝度へグラデーションするような色のパッチが印刷等された色票を使用し、複数の輝度レベルの画像データの撮像を１枚の色票で実現しても構わない。また、輝度レベルの異なる複数の色票の撮像を行っても構わない。また、自発光する色票である自発光チャートを利用し、多段階の発光強度で撮像を行っても構わない。本ステップにおける画像データ群の入力は、ケーブルなどを介してデジタルカメラ１から直接読み込むことで行ってもよく、メモリカードなどの記録媒体を介して読み込むことで行ってもよい。読み込んだ画像データ群は、データ保持部２０７に格納される。
【0023】
次にＳ４０２において、平均値算出部２０３は、データ保持部２０７から読み出された画像データ群の各画像における、画像内のパッチごとのＬ*ａ*ｂ*値の平均値を、計測値として算出する。なお、各画像のＲＧＢ値をＬ*ａ*ｂ*値に変換する際には、後述する色空間変換式を利用することができる。ここで、パッチ数を１８、露光時間の段階数を３２とすると、ソース色空間の画像とターゲット色空間の画像のそれぞれから、１８×３２の計測点における計測値が得られることになる。図５は、ソース色空間の画像から得られた計測値を、Ｌ*ａ*ｂ*の色空間上にプロットした例を示す。以下、ソース色空間の画像から得られた計測値をソース色空間の学習データと呼び、ターゲット色空間の画像から得られた計測値をターゲット色空間の学習データと呼ぶ。これらの学習データは、データ保持部２０７に格納される。ソース色空間の画像とターゲット色空間の画像のそれぞれから計測値を取得する処理は、データ取得処理の一例である。
【0024】
次にＳ４０３において、計測値取得部３０１は、データ保持部２０７から学習データを読み出す。そして、読み出したソース色空間の学習データとターゲット色空間の学習データについて対応する各計測点における、計測値間のＬ*ａ*ｂ*の成分ごとの変分Δｌ、Δａ、Δｂを算出する。変分Δｌ、Δａ、Δｂは、下記の式（１）～（３）で求められる。
Δｌ＝ｌ´―ｌ・・・（１）
Δａ＝ａ´―ａ・・・（２）
Δｂ＝ｂ´―ｂ・・・（３）
但し、（ｌ、ａ、ｂ）は、ソース色空間における学習データの計測値（Ｌ*ａ*ｂ*値）を表す。また、（ｌ´、ａ´、ｂ´）は、ターゲット色空間における学習データの計測値（Ｌ*ａ*ｂ*値）を表す。（ｌ´、ａ´、ｂ´）は、（ｌ、ａ、ｂ）に対応する計測点を表している。
Δｌ＝ｌ´―ｌ・・・（１）
Δａ＝ａ´―ａ・・・（２）
Δｂ＝ｂ´―ｂ・・・（３）
但し、（ｌ、ａ、ｂ）は、ソース色空間における学習データの計測値（Ｌ*ａ*ｂ*値）を表す。また、（ｌ´、ａ´、ｂ´）は、ターゲット色空間における学習データの計測値（Ｌ*ａ*ｂ*値）を表す。（ｌ´、ａ´、ｂ´）は、（ｌ、ａ、ｂ）に対応する計測点を表している。
【0025】
次にＳ４０４において、計測値補間部３０２は、色票３のパッチ（色）ごとに、Ｓ４０３で得られた各計測点における変分Δｌ、Δａ、Δｂについて、輝度ｌの関数で表現する式を求める。例えば計測値補間部３０２は、成分ａに関して、色票３のｉ番目のパッチ（色）の目的変数Δａiを説明変数ｌで説明する関数Δａi（ｌ）を求める。また、計測値補間部３０２は、色票３のパッチ（色）ごとに、各計測点におけるソース色空間の学習データの成分ｌ、ａ、ｂを、輝度ｌの関数で表現する式も同様に求めておく。以上のようにして計測値補間部３０２は、関数ｌi（ｌ）、ａi（ｌ）、ｂi（ｌ）、Δｌi（ｌ）、Δａi（ｌ）、Δｂi（ｌ）を生成する。なお、ここでｉはパッチの番号を表し、パッチ数が１８である場合、１≦ｉ≦１８である。
【0026】
ここで関数Δａi（ｌ）を求める方法について具体的に説明する。まず、計測点は離散的にしか存在しない。例えば、露出時間を３２段階で撮像した場合、色票３のパッチ（色）ごとの計測点数は３２点である。そこで、離散的に存在する計測点における変分値を、順番に滑らかに繋ぐ曲線を求める。図６は、色票３のｉ番目のパッチ（色）についての関数Δａi（ｌ）を表すグラフの例を示す。図６のグラフの横軸は輝度ｌを表し、縦軸は変分Δａを表す。丸印は計測点を表し、曲線は補間結果を表す。図６の例では、計測点を誤差なく通過するような7次多項式を最小２乗法により求めることで計測点間の補間を行っている。なお、多項式の次数は７以外でも構わない。また、多項式近似を行う代わりに、計測点の間を折れ線補間、スプライン補間、ベジエ補間など他の方法により補間を行っても構わない。また、必ずしも補間を行わなくてもよく、例えば計測点間の値は、最近傍の計測点の値と同じ値であると定義しても構わない。
【0027】
次にＳ４０５において、入力点取得部３０３は、色変換対象の入力点Ｐを取得する。ここで入力点Ｐは、ｌ、ａ、ｂの３要素で表される。本実施形態において入力点Ｐは、３ＤＬＵＴの格子点に対応している。３ＤＬＵＴのグリッド数が１７と仮定した場合、１７×１７×１７＝４９１３個の入力点Ｐ1，Ｐ2，・・・、Ｐ4913が存在することになる。つまり、以下のＳ４０６～Ｓ４０９に示す色変換処理が、４９１３個の入力点のそれぞれに対して繰り返し実行されることになる。
【0028】
なお、入力点Ｐ1，Ｐ2，・・・、Ｐ4913の求め方は次の通りである。まず、デジタルカメラ１が出力する画像のＲＧＢ値が、例えばＲＧＢの各色１０ビット（０～１０２３）であり、デジタルカメラ１の３ＤＬＵＴのグリッド数が１７であると仮定した場合、当該３ＤＬＵＴの各格子点は、以下のように表される。
(Ｒ,Ｇ,Ｂ)＝(0, 0, 0) (0, 0, 64) (0, 0, 128)
… (0, 0, 1023) (0, 64, 0) (0, 64, 64)
… (0, 64, 1023) (0, 128, 0) (0, 128, 64)
… (1023, 1023, 1023)
続いて入力点取得部３０３は、上記の３ＤＬＵＴの各格子点のＲＧＢ値、Ｌ*ａ*ｂ*値に変換する。これにより、入力点Ｐ1，Ｐ2，・・・，Ｐ4913が得られる。ＲＧＢ値をＬ*ａ*ｂ*値に変換する際には、後述する色空間変換式を利用することができる。
(Ｒ,Ｇ,Ｂ)＝(0, 0, 0) (0, 0, 64) (0, 0, 128)
… (0, 0, 1023) (0, 64, 0) (0, 64, 64)
… (0, 64, 1023) (0, 128, 0) (0, 128, 64)
… (1023, 1023, 1023)
続いて入力点取得部３０３は、上記の３ＤＬＵＴの各格子点のＲＧＢ値、Ｌ*ａ*ｂ*値に変換する。これにより、入力点Ｐ1，Ｐ2，・・・，Ｐ4913が得られる。ＲＧＢ値をＬ*ａ*ｂ*値に変換する際には、後述する色空間変換式を利用することができる。
【0029】
次にＳ４０６において、基準点選択部３０４は、Ｓ４０４にてソース色空間における学習データから得られた関数ｌi（ｌ）、ａi（ｌ）、ｂi（ｌ）を用いて、Ｓ４０５で取得した入力点Ｐと色相の近い基準点を２点選択する。以下、図７を用いて、基準点を選択する方法について説明する。
まず、基準点選択部３０４は、入力点Ｐの輝度値ｌ0を、Ｓ４０４で求めた関数ａi（ｌ）、ｂi（ｌ）に代入することで、下記の式（４）に示すように、色票３の全てのパッチ（色）に対応する、色彩値Ｐi（ｌ0）を求める。
Ｐi（ｌ0）＝（ｌ0，ａi（ｌ0），ｂi（ｌ0）） ・・・（４）
図７は、輝度値ｌ0であるときのａ*ｂ*色平面に、色彩値Ｐi（ｌ0）をプロットした例を示す。図７の例では、ｉ＝１，２，・・・，１８についての色彩値Ｐｉ（ｌ0）が表れている。これらのＰi（ｌ0）を、基準点候補と呼ぶ。点７０２は、入力点Ｐを表す。
まず、基準点選択部３０４は、入力点Ｐの輝度値ｌ0を、Ｓ４０４で求めた関数ａi（ｌ）、ｂi（ｌ）に代入することで、下記の式（４）に示すように、色票３の全てのパッチ（色）に対応する、色彩値Ｐi（ｌ0）を求める。
Ｐi（ｌ0）＝（ｌ0，ａi（ｌ0），ｂi（ｌ0）） ・・・（４）
図７は、輝度値ｌ0であるときのａ*ｂ*色平面に、色彩値Ｐi（ｌ0）をプロットした例を示す。図７の例では、ｉ＝１，２，・・・，１８についての色彩値Ｐｉ（ｌ0）が表れている。これらのＰi（ｌ0）を、基準点候補と呼ぶ。点７０２は、入力点Ｐを表す。
【0030】
次に、基準点選択部３０４は、入力点Ｐと色相が近い基準点候補を２点選択する。具体的には、輝度値ｌ0であるときのａ*ｂ*色平面において、点７０２から反時計回りに見て最も色相が近い基準点候補である基準点候補７０３と、時計回りに見て最も色相が近い基準点候補７０４を選択し、これらを基準点とする。ここでは、上記の式（４）で表される基準点候補のうち、ζ番目の色票を表すＰζ（ｌ0）＝（ｌ0，ａζ（ｌ0），ｂζ（ｌ0））と、η番目の色票を表すＰη（ｌ0）＝（ｌ0，ａη（ｌ0），ｂη（ｌ0））が、基準点として選択されたとする。このように、複数の色票に対応する複数の基準点候補の中から、複数の基準点（ここでは、Ｐζ（ｌ0）、Ｐη（ｌ0））を選択する処理は、基準点選択処理の一例である。
【0031】
なお、入力点Ｐと色空間上の距離が近い複数の基準点を選択する方法であれば、基準点を選択する方法は上記の方法に限定されない。例えば、時計回りや反時計回りといった方向を無視して入力点Ｐと色相の近い基準点候補を２点選択してもよい。また、色相差でなく、彩度差や色差が小さい点を選択してもよい。また、入力点Ｐとの差分の計算を、別の色空間で行っても構わない。また、入力点Ｐと基準点候補の相対位置に基づいて基準点の選択を行うのではなく、入力点Ｐの位置する領域に応じて基準点を選択するなど、入力点Ｐの絶対位置に基づいて基準点の選択を行っても構わない。
【0032】
次にＳ４０７において、変換係数算出部３０５は、入力点Ｐの位置を、Ｓ４０６において選択した基準点の位置で表現する変換係数を算出する。以下、変換係数を算出する方法について具体的に説明する。本実施形態では、変換係数算出部３０５は、入力点Ｐの位置を２つの基準点と原点との相対位置で表す線形結合係数を、変換係数として算出する。
図８は、入力点Ｐと基準点Ｐζ（ｌ0）、Ｐη（ｌ0）と原点Ｏの相対的な位置関係を示す図である。図８では、（ｌ，ａ，ｂ）＝（ｌ0，０，０）の位置を原点Ｏとして表す。この時、原点Ｏから入力点Ｐに向かうベクトルＯＰは、原点Ｏから基準点Ｐζ（ｌ0）に向かうベクトルＯＰζと、原点Ｏから基準点Ｐη（ｌ0）に向かうベクトルＯＰηの線形結合で表される。このことを、下記の式（５）で表す。
図８は、入力点Ｐと基準点Ｐζ（ｌ0）、Ｐη（ｌ0）と原点Ｏの相対的な位置関係を示す図である。図８では、（ｌ，ａ，ｂ）＝（ｌ0，０，０）の位置を原点Ｏとして表す。この時、原点Ｏから入力点Ｐに向かうベクトルＯＰは、原点Ｏから基準点Ｐζ（ｌ0）に向かうベクトルＯＰζと、原点Ｏから基準点Ｐη（ｌ0）に向かうベクトルＯＰηの線形結合で表される。このことを、下記の式（５）で表す。
【0033】
【数1】
ａ，ｂ：入力点Ｐのａ成分、ｂ成分の値
ａζ，ｂζ：基準点Ｐζ（ｌ0）のａ成分、ｂ成分の値
ａη，ｂη：基準点Ｐη（ｌ0）のａ成分、ｂ成分の値
α：ベクトルＯＰζの線形結合係数
β：ベクトルＯＰηの線形結合係数
上記の式（５）の線形結合係数α、βは、下記の式（６）で導出される。
【数2】
【0034】
次にＳ４０８において、変分値取得部３０６は、入力点Ｐにおける変分ΔＰの推定値を算出する。まず、変分値取得部３０６は、入力点Ｐの輝度値ｌ0を、Ｓ４０４で求めた関数Δｌi（ｌ）、Δａi（ｌ）、Δｂi（ｌ）に代入することで、下記の式（７）に示すように、色票３の全てのパッチ（色）における、変分ΔＰiを求める。
ΔＰi（ｌ0）＝（Δｌi（ｌ0）、Δａi（ｌ0）、Δｂi（ｌ0）） ・・・（７）
そして、変分値取得部３０６は、Ｓ４０７で算出した線形結合係数α、βを用いて、以下の式（８）から、入力点Ｐにおける変分ΔＰの推定値を算出する。
ΔＰ＝αΔＰζ（ｌ0）＋βΔＰη（ｌ0） ・・・（８）
上記の式（８）のΔＰζ（ｌ0）は、上記の式（７）でｉ＝ζのときのΔＰζ（ｌ0）であり、基準点Ｐζ（ｌ0）に対応する変分値である。また、上記の式（８）のΔＰη（ｌ0）は、上記の式（７）でｉ＝ηのときのΔＰη（ｌ0）であり、基準点Ｐη（ｌ0）に対応する変分値である。
ΔＰi（ｌ0）＝（Δｌi（ｌ0）、Δａi（ｌ0）、Δｂi（ｌ0）） ・・・（７）
そして、変分値取得部３０６は、Ｓ４０７で算出した線形結合係数α、βを用いて、以下の式（８）から、入力点Ｐにおける変分ΔＰの推定値を算出する。
ΔＰ＝αΔＰζ（ｌ0）＋βΔＰη（ｌ0） ・・・（８）
上記の式（８）のΔＰζ（ｌ0）は、上記の式（７）でｉ＝ζのときのΔＰζ（ｌ0）であり、基準点Ｐζ（ｌ0）に対応する変分値である。また、上記の式（８）のΔＰη（ｌ0）は、上記の式（７）でｉ＝ηのときのΔＰη（ｌ0）であり、基準点Ｐη（ｌ0）に対応する変分値である。
【0035】
次にＳ４０９において、色変換結果算出部３０７は、Ｓ４０８で算出された変分ΔＰの推定値を用いて、入力点Ｐに対する色変換結果Ｐ´を算出する。算出式は、下記の式（９）で表される。
Ｐ´＝Ｐ＋ΔＰ ・・・（９）
続いて、画像処理部２０２は、Ｓ４０５の処理に戻り、４９１３個の入力点Ｐのうち、次の入力点Ｐを色変換対象として取得して、Ｓ４０６～Ｓ４０９に示す色変換処理を実行する。その後すべての入力点Ｐに対してＳ４０６～Ｓ４０９に示す色変換処理が完了すると、処理はＳ４１０に遷移する。
Ｐ´＝Ｐ＋ΔＰ ・・・（９）
続いて、画像処理部２０２は、Ｓ４０５の処理に戻り、４９１３個の入力点Ｐのうち、次の入力点Ｐを色変換対象として取得して、Ｓ４０６～Ｓ４０９に示す色変換処理を実行する。その後すべての入力点Ｐに対してＳ４０６～Ｓ４０９に示す色変換処理が完了すると、処理はＳ４１０に遷移する。
【0036】
Ｓ４１０において、３ＤＬＵＴ生成部２０８は、Ｓ４０５において取得された入力点Ｐと、その入力点Ｐに対してＳ４０９において算出された色変換結果Ｐ´の組を生成する。
Ｔ＝｛（Ｐ1，Ｐ1´），（Ｐ2，Ｐ2´），…，（Ｐ4913，Ｐ4913´）｝ ・・・（１０）
上記の式（１０）のＴは、３ＤＬＵＴの各格子点における色変換結果を保持する。上記の式（１０）のＴが、３ＤＬＵＴに相当する。
Ｔ＝｛（Ｐ1，Ｐ1´），（Ｐ2，Ｐ2´），…，（Ｐ4913，Ｐ4913´）｝ ・・・（１０）
上記の式（１０）のＴは、３ＤＬＵＴの各格子点における色変換結果を保持する。上記の式（１０）のＴが、３ＤＬＵＴに相当する。
【0037】
次にＳ４１１において、データ入出力部２０１は、Ｓ４１０で生成された３ＤＬＵＴを、汎用Ｉ／Ｆ２１５を介してデジタルカメラ１にアップロードする。以上のようにして、図４に示すフローチャートの一連の処理が終了する。
【0038】
以下、図４に示すフローチャートで生成された３ＤＬＵＴを用いて、デジタルカメラ１で実際に撮像して得られた画像を色変換する方法を示す。ここでは、色変換する処理の主体がデジタルカメラ１の画像処理部１０３であるとして説明するが、画像処理装置２の画像処理部２０２であっても構わない。
【0039】
まず、任意の入力画像データの画素値Ｑ（ＲＧＢ値）が入力された際、画像処理部１０３は、この画素値Ｑを含む四面体を構成する４点ＰAＰBＰCＰDを、上記の式（１０）のＴが表す格子点の中から選択して、それらを対応格子点４点として確定する。
次に画像処理部１０３は、補間比率を算出する。まず画像処理部１０３は、四面体を構成する対応格子点４点のうち１点を基準とし、４つのベクトル（ＰAＰB）、(ＰAＰC)、（ＰAＰD）、（ＰAＱ）を作成する。そして画像処理部１０３は、作成した４つのベクトル（ＰAＰB）、(ＰAＰC)、（ＰAＰD）、（ＰAＱ）と、任意の係数ｗB、ｗC、ｗDを用いた下記の式（１１）を満たすｗB、ｗC、ｗDを求める。このｗB、ｗC、ｗDは、補間比率である。
（ＰAＱ）＝ｗB（ＰAＰB）＋ｗC(ＰAＰC)＋ｗD（ＰAＰD） ・・・（１１）
続いて画像処理部１０３は、上記の式（１１）で求めた補間比率ｗB、ｗC、ｗDを用いて、色変換結果Ｑ´を下記の式（１２）で算出する。
Ｑ´＝ＰA´＋ｗB（ＰA´ＰB´）＋ｗC(ＰA´ＰC´)＋ｗD（ＰA´ＰD´） ・・・（１２）
以上のようにして、デジタルカメラ１で撮像して得られた画像を、３ＤＬＵＴを用いて色変換することができる。このように、３ＤＬＵＴとして色変換関数を保持することで、計算時間とメモリ量の削減を図ることができる。
次に画像処理部１０３は、補間比率を算出する。まず画像処理部１０３は、四面体を構成する対応格子点４点のうち１点を基準とし、４つのベクトル（ＰAＰB）、(ＰAＰC)、（ＰAＰD）、（ＰAＱ）を作成する。そして画像処理部１０３は、作成した４つのベクトル（ＰAＰB）、(ＰAＰC)、（ＰAＰD）、（ＰAＱ）と、任意の係数ｗB、ｗC、ｗDを用いた下記の式（１１）を満たすｗB、ｗC、ｗDを求める。このｗB、ｗC、ｗDは、補間比率である。
（ＰAＱ）＝ｗB（ＰAＰB）＋ｗC(ＰAＰC)＋ｗD（ＰAＰD） ・・・（１１）
続いて画像処理部１０３は、上記の式（１１）で求めた補間比率ｗB、ｗC、ｗDを用いて、色変換結果Ｑ´を下記の式（１２）で算出する。
Ｑ´＝ＰA´＋ｗB（ＰA´ＰB´）＋ｗC(ＰA´ＰC´)＋ｗD（ＰA´ＰD´） ・・・（１２）
以上のようにして、デジタルカメラ１で撮像して得られた画像を、３ＤＬＵＴを用いて色変換することができる。このように、３ＤＬＵＴとして色変換関数を保持することで、計算時間とメモリ量の削減を図ることができる。
【0040】
続いて、Ｓ４０２やＳ４０５などで利用される色空間変換式（ｓＲＧＢ→ＣＩＥＸＹＺ変換式、およびＣＩＥＸＹＺ→ＣＩＥＬＡＢ変換式）について説明する。
まず、ｓＲＧＢ→ＣＩＥＸＹＺ変換式について説明する。
画像処理部２０２は、入力データのＲＧＢ値を正規化する。ここでの正規化は、入力画像データのＲＧＢ値が１０ビットの場合には、各値を２10－１＝１０２３で除算する処理であるとする。そして、画像処理部２０２は、以下の式(１３)～式(１７)を用いてＸＹＺ値を算出する。
ＲsRGB≦0.04045 のとき Ｒliner＝ＲsRGB/12.92
ＲsRGB＞0.04045 のとき Ｒliner＝[(ＲsRGB+0.055)/1.055)]2.4 …（１３）
ＧsRGB≦0.04045 のとき Ｇliner＝ＧsRGB/12.92
ＧsRGB＞0.04045 のとき Ｇliner＝[(ＧsRGB+0.055)/1.055)]2.4 …（１４）
ＢsRGB≦0.04045 のとき Ｂliner＝ＢsRGB/12.92
ＢsRGB＞0.04045 のとき Ｂliner＝[(ＢsRGB+0.055)/1.055)]2.4 …（１５）
まず、ｓＲＧＢ→ＣＩＥＸＹＺ変換式について説明する。
画像処理部２０２は、入力データのＲＧＢ値を正規化する。ここでの正規化は、入力画像データのＲＧＢ値が１０ビットの場合には、各値を２10－１＝１０２３で除算する処理であるとする。そして、画像処理部２０２は、以下の式(１３)～式(１７)を用いてＸＹＺ値を算出する。
ＲsRGB≦0.04045 のとき Ｒliner＝ＲsRGB/12.92
ＲsRGB＞0.04045 のとき Ｒliner＝[(ＲsRGB+0.055)/1.055)]2.4 …（１３）
ＧsRGB≦0.04045 のとき Ｇliner＝ＧsRGB/12.92
ＧsRGB＞0.04045 のとき Ｇliner＝[(ＧsRGB+0.055)/1.055)]2.4 …（１４）
ＢsRGB≦0.04045 のとき Ｂliner＝ＢsRGB/12.92
ＢsRGB＞0.04045 のとき Ｂliner＝[(ＢsRGB+0.055)/1.055)]2.4 …（１５）
【0041】
【数3】
【0042】
Ｘ＝Ｘ・Ｘｎ／95.045
Ｙ＝Ｙ・Ｙｎ／100.000
Ｚ＝Ｚ・Ｚｎ／108.892 …(１７)
ここで、Ｘｎ、Ｙｎ、Ｚｎは撮像時の照明の色温度から算出した値である。
Ｙ＝Ｙ・Ｙｎ／100.000
Ｚ＝Ｚ・Ｚｎ／108.892 …(１７)
ここで、Ｘｎ、Ｙｎ、Ｚｎは撮像時の照明の色温度から算出した値である。
【0043】
次に、ＣＩＥＸＹＺ→ＣＩＥＬＡＢ変換式について説明する。
ＣＩＥＸＹＺ→ＣＩＥＬＡＢ変換式は、以下の式(１８)～式(２２)により表される。
Ｘ/Ｘｎ＞0.00856 のとき ＸRate＝(Ｘ/Ｘｎ)1/3
Ｘ/Ｘｎ≦0.00856 のとき ＸRate＝7.787×(Ｘ/Ｘｎ)+16.0/116.0 …（１８）
Ｙ/Ｙｎ＞0.00856 のとき ＹRate＝(Ｙ/Ｙｎ)1/3
Ｙ/Ｙｎ≦0.00856 のとき ＹRate＝7.787×(Ｙ/Ｙｎ)+16.0/116.0 …（１９）
Ｚ/Ｚｎ＞0.00856 のとき ＺRate＝(Ｚ/Ｚｎ)1/3
Ｚ/Ｚｎ≦0.00856 のとき ＺRate＝7.787×(Ｚ/Ｚｎ)+16.0/116.0 …（２０）
Ｙ/Ｙｎ＞0.00856 のとき Ｌ*＝116.0×(Ｙ/Ｙｎ)1/3-16.0
Ｙ/Ｙｎ≦0.00856 のとき Ｌ*＝903.29×(Ｙ/Ｙｎ) …（２１）
ａ*＝500×(ＸRate-ＹRate)
ｂ*＝200×(ＹRate-ＺRate) …（２２）
ここで、Ｘｎ、Ｙｎ、Ｚｎは撮像時の照明の色温度から算出した値である。
ＣＩＥＸＹＺ→ＣＩＥＬＡＢ変換式は、以下の式(１８)～式(２２)により表される。
Ｘ/Ｘｎ＞0.00856 のとき ＸRate＝(Ｘ/Ｘｎ)1/3
Ｘ/Ｘｎ≦0.00856 のとき ＸRate＝7.787×(Ｘ/Ｘｎ)+16.0/116.0 …（１８）
Ｙ/Ｙｎ＞0.00856 のとき ＹRate＝(Ｙ/Ｙｎ)1/3
Ｙ/Ｙｎ≦0.00856 のとき ＹRate＝7.787×(Ｙ/Ｙｎ)+16.0/116.0 …（１９）
Ｚ/Ｚｎ＞0.00856 のとき ＺRate＝(Ｚ/Ｚｎ)1/3
Ｚ/Ｚｎ≦0.00856 のとき ＺRate＝7.787×(Ｚ/Ｚｎ)+16.0/116.0 …（２０）
Ｙ/Ｙｎ＞0.00856 のとき Ｌ*＝116.0×(Ｙ/Ｙｎ)1/3-16.0
Ｙ/Ｙｎ≦0.00856 のとき Ｌ*＝903.29×(Ｙ/Ｙｎ) …（２１）
ａ*＝500×(ＸRate-ＹRate)
ｂ*＝200×(ＹRate-ＺRate) …（２２）
ここで、Ｘｎ、Ｙｎ、Ｚｎは撮像時の照明の色温度から算出した値である。
【0044】
なお、本実施形態では、デジタルカメラ１と画像処理装置２は、ケーブルを用いてデータのやり取りを行うが、勿論、メモリカードなどの記録媒体、無線や赤外線などの通信路を用いてデータをやり取りしてもよい。
【0045】
以上説明したように、実施形態１によれば、色票を撮像した撮像画像データが、理想の発色をもつカメラで色票を撮像したデータに色変換されるよう最適化された３ＤＬＵＴを生成することができる。この３ＤＬＵＴを用いて色変換処理を行うことで、撮像画像データを色変換する際に生じる変換誤差を低減することが可能になる。これにより、本実施形態によれば、僅かな入力色の変化に対して出力色が大きく変化するような条件下においても、変換誤差を低減することができる。
【0046】
［変形例１］
上述の実施形態では、デジタルカメラ１からの入力データのＲＧＢ値をｓＲＧＢと仮定しＬ*ａ*ｂ*色空間で各種の計算を行っているが、他の色空間で計算を行っても構わない。例えば、入力データのＲＧＢ値をＢＴ．２０２０ＲＧＢと仮定しＩＣｔＣｐ色空間で計算を行っても構わない。その場合、画像処理部２０２は、以下の色空間変換式を利用する。
上述の実施形態では、デジタルカメラ１からの入力データのＲＧＢ値をｓＲＧＢと仮定しＬ*ａ*ｂ*色空間で各種の計算を行っているが、他の色空間で計算を行っても構わない。例えば、入力データのＲＧＢ値をＢＴ．２０２０ＲＧＢと仮定しＩＣｔＣｐ色空間で計算を行っても構わない。その場合、画像処理部２０２は、以下の色空間変換式を利用する。
【0047】
ＢＴ．２０２０ＲＧＢ→ＣＩＥＸＹＺ変換式は、以下の式(２３)～式(２５)により表される。
Ｒliner＝((Ｒ1/m-c1)/(c2-c3*Ｒ1/m))1/n
Ｇliner＝((Ｇ1/m-c1)/(c2-c3*Ｇ1/m))1/n
Ｂliner＝((Ｂ1/m-c1)/(c2-c3*Ｂ1/m))1/n …（２３）
m＝78.8437500;
n＝0.15930176;
c1＝0.83593750;
c2＝18.8515625;
c3＝18.6875000;
Ｒliner＝((Ｒ1/m-c1)/(c2-c3*Ｒ1/m))1/n
Ｇliner＝((Ｇ1/m-c1)/(c2-c3*Ｇ1/m))1/n
Ｂliner＝((Ｂ1/m-c1)/(c2-c3*Ｂ1/m))1/n …（２３）
m＝78.8437500;
n＝0.15930176;
c1＝0.83593750;
c2＝18.8515625;
c3＝18.6875000;
【0048】
【数4】
【0049】
Ｘ＝Ｘ・Ｘｎ/95.045
Ｙ＝Ｙ・Ｙｎ/100.000
Ｚ＝Ｚ・Ｚｎ/108.892 …（２５）
ここで、Ｘｎ、Ｙｎ、Ｚｎは撮像時の照明の色温度から算出した値である。
Ｙ＝Ｙ・Ｙｎ/100.000
Ｚ＝Ｚ・Ｚｎ/108.892 …（２５）
ここで、Ｘｎ、Ｙｎ、Ｚｎは撮像時の照明の色温度から算出した値である。
【0050】
ＣＩＥＸＹＺ→ＩＣｔＣｐ変換式は、以下の式(２６)～式(２８)により表される。
【0051】
【数5】
【0052】
ＬＰＱ＝((c1+c2*Ｌn)/(1+c3*Ｌn))m
ＭＰＱ＝((c1+c2*Ｍn)/(1+c3*Ｍn))m
ＳＰＱ＝((c1+c2*Ｓn)/(1+c3*Ｓn))m …（２７）
m＝78.8437500;
n＝0.15930176;
c1＝0.83593750;
c2＝18.8515625;
c3＝18.6875000;
ＭＰＱ＝((c1+c2*Ｍn)/(1+c3*Ｍn))m
ＳＰＱ＝((c1+c2*Ｓn)/(1+c3*Ｓn))m …（２７）
m＝78.8437500;
n＝0.15930176;
c1＝0.83593750;
c2＝18.8515625;
c3＝18.6875000;
【0053】
【数6】
【0054】
［変形例２］
上述の実施形態では、実際に撮像した画像の画素値Ｑが入力された場合に、３ＤＬＵＴを用いて色変換結果Ｑ´を算出する例について説明した。本変形例では、Ｓ４０５～Ｓ４０９までの処理を実行して、色変換結果Ｑ´を直接算出する例を示す。本変形例は、Ｓ４０５～Ｓ４０９についてのみ説明し、他のステップについては実施形態１と同様であるとして説明を割愛する。
上述の実施形態では、実際に撮像した画像の画素値Ｑが入力された場合に、３ＤＬＵＴを用いて色変換結果Ｑ´を算出する例について説明した。本変形例では、Ｓ４０５～Ｓ４０９までの処理を実行して、色変換結果Ｑ´を直接算出する例を示す。本変形例は、Ｓ４０５～Ｓ４０９についてのみ説明し、他のステップについては実施形態１と同様であるとして説明を割愛する。
【0055】
まずＳ４０５において、入力点取得部３０３は、ソース色空間における色変換対象の画素の一つを点Ｑとして取得する。続いてＳ４０６～Ｓ４０８と同様の処理により、画像処理部２０２は、入力点Ｑにおける変分ΔＱの推定値を算出する。次のＳ４０９において、色変換結果算出部３０７は、色変換結果Ｑ´＝Ｑ＋ΔＱを算出する。その後、画像処理部２０２は、色変換対象の画素が無くなるまで、Ｓ４０５～Ｓ４０９を繰り返し実行する。なお、その後のＳ４１０～Ｓ４１１の処理は省略される。
このように本変形例では、３ＤＬＵＴを用いることなく、色変換結果を算出することができる。
このように本変形例では、３ＤＬＵＴを用いることなく、色変換結果を算出することができる。
【0056】
［変形例３］
上述の実施形態では、Ｓ４０４で求めた変分の関数Δｌi（ｌ）、Δａi（ｌ）、Δｂi（ｌ）を用いて色変換結果Ｐ´を算出する例について説明した。本変形例では、変分の関数Δｌi（ｌ）、Δａi（ｌ）、Δｂi（ｌ）を求めることなく、色変換結果Ｐ´を算出する例を示す。本変形例は、Ｓ４０４～Ｓ４０９についてのみ説明し、他のステップについては実施形態１と同様であるとして説明を割愛する。
上述の実施形態では、Ｓ４０４で求めた変分の関数Δｌi（ｌ）、Δａi（ｌ）、Δｂi（ｌ）を用いて色変換結果Ｐ´を算出する例について説明した。本変形例では、変分の関数Δｌi（ｌ）、Δａi（ｌ）、Δｂi（ｌ）を求めることなく、色変換結果Ｐ´を算出する例を示す。本変形例は、Ｓ４０４～Ｓ４０９についてのみ説明し、他のステップについては実施形態１と同様であるとして説明を割愛する。
【0057】
まずＳ４０４において、計測値補間部３０２は、色票３のパッチ（色）ごとに、各計測点におけるターゲット色空間の学習データの成分ｌ´、ａ´、ｂ´から、輝度ｌを説明変数とする関数ｌ´i（ｌ）、ａ´i（ｌ）、ｂ´i（ｌ）を求めておく。そして、計測値補間部３０２は、これら３要素を持つ関数Ｐ´i（ｌ）を生成する。続いてＳ４０５～Ｓ４０７が実行される。本変形例では、次のＳ４０８～Ｓ４０９に代えて、画像処理部２０２は、下記の式（２９）を用いて、色変換結果Ｐ´を算出する。なお、下記の式（２９）のα、βは、Ｓ４０７において算出された線形結合係数である。また、下記の式（２９）のＰζ´（ｌ0）、Ｐη´（ｌ0）は、Ｓ４０６で選択された基準点である。
Ｐ´＝αＰζ´（ｌ0）＋βＰη´（ｌ0） …（２９）
このように本変形例では、変分の関数Δｌi（ｌ）、Δａi（ｌ）、Δｂi（ｌ）を求めることなく、色変換結果を算出することができる。
Ｐ´＝αＰζ´（ｌ0）＋βＰη´（ｌ0） …（２９）
このように本変形例では、変分の関数Δｌi（ｌ）、Δａi（ｌ）、Δｂi（ｌ）を求めることなく、色変換結果を算出することができる。
【0058】
［実施形態２］
実施形態１では、計測値取得部３０１において取得された、すべての学習データを用いて色変換処理を行う例について説明した。一方で、不適切な一部の学習データを色変換処理の対象から除去することで、色変換の精度が向上する場合がある。図９は、計測値取得部３０１が取得した学習データの例を示す。図９に示すグラフは、学習データの色ごとの補間結果をａ*ｂ*色平面に射影したものであり、系列１本がそれぞれ色票のｉ番目の色の色彩値Ｐi（ｌ）の形状データを表す。ここで、系列９０１は、デジタルカメラ１での撮像時の誤差やソース色空間自体の歪みが原因で、他の系列と交叉してしまっている系列である。この系列を学習データとして使用すると、色変換結果画像に階調反転や疑似輪郭が発生するなどの画質弊害が起こりうる。従って、学習データから系列９０１を除去し、それ以外の系列のみを使用することが望ましい。以下では、実施形態１と異なる部分について主として説明する。また、実施形態１と同一の構成については、同じ符号を付して説明する。
実施形態１では、計測値取得部３０１において取得された、すべての学習データを用いて色変換処理を行う例について説明した。一方で、不適切な一部の学習データを色変換処理の対象から除去することで、色変換の精度が向上する場合がある。図９は、計測値取得部３０１が取得した学習データの例を示す。図９に示すグラフは、学習データの色ごとの補間結果をａ*ｂ*色平面に射影したものであり、系列１本がそれぞれ色票のｉ番目の色の色彩値Ｐi（ｌ）の形状データを表す。ここで、系列９０１は、デジタルカメラ１での撮像時の誤差やソース色空間自体の歪みが原因で、他の系列と交叉してしまっている系列である。この系列を学習データとして使用すると、色変換結果画像に階調反転や疑似輪郭が発生するなどの画質弊害が起こりうる。従って、学習データから系列９０１を除去し、それ以外の系列のみを使用することが望ましい。以下では、実施形態１と異なる部分について主として説明する。また、実施形態１と同一の構成については、同じ符号を付して説明する。
【0059】
図１０は、本実施形態により実行される画像処理装置２の処理の流れを示すフローチャートである。本フローチャートの処理は、図４のフローチャートとは、Ｓ４０４とＳ４０５との間に、下記に説明するＳ１００１～Ｓ１００３の処理が実行される点で相違する。
まずＳ１０００において、Ｓ４０１～Ｓ４０４と同様の処理が実行される。
まずＳ１０００において、Ｓ４０１～Ｓ４０４と同様の処理が実行される。
【0060】
次にＳ１００１において、計測値補間部３０２は、Ｓ４０４で生成した関数形式のソース色空間の学習データを、データ保持部２０７に格納する。続いて、表示制御部２０４は、データ保持部２０７からソース色空間の学習データ（関数形式）を読み出し、表示装置２１４に出力することで、ユーザに提示する。図１１は、表示装置２１４に出力されたＧＵＩの例を示す。このＧＵＩは、ウインドウ１１０１内に、学習データ形状１１０２を表示する。この学習データ形状１１０２は、ソース色空間の学習データ（関数形式）をａ*ｂ*色平面に射影したものであり、系列１本がそれぞれ色票のｉ番目のパッチ（色）の色彩値Ｐi（ｌ）の形状データを表す。
【0061】
次にＳ１００２において、ユーザ操作取得部２０５は、学習データの中からユーザが削除したい系列の操作指示を受け付ける。系列が選択されている間、表示制御部２０４は、系列１１０３のように、選択された系列を強調表示する。なお、ユーザ操作取得部２０５は、選択のやり直しや複数選択などの機能を備えていてもよい。本ステップの操作は、マウス、キーボード、タッチパネルなどの入力装置２２０を介して行われる。ユーザ操作取得部２０５は、削除ボタン１１０４が押下された場合、処理をＳ１００３に進める。なお、削除ボタン１１０４が押下されなかった場合には、Ｓ１００３をスキップして、処理はＳ１００４に進む。
Ｓ１００３において、データ保持部２０７は、Ｓ１００３で削除指示された系列に対応する学習データを削除する。なお、データ保持部２０７は、Ｓ１００３で削除指示された系列に対応する学習データを破棄するのではなく、別のメモリに別途保存しておいてもよい。その後処理はＳ１００２に戻る。
Ｓ１００３において、データ保持部２０７は、Ｓ１００３で削除指示された系列に対応する学習データを削除する。なお、データ保持部２０７は、Ｓ１００３で削除指示された系列に対応する学習データを破棄するのではなく、別のメモリに別途保存しておいてもよい。その後処理はＳ１００２に戻る。
【0062】
上記の削除処理は、任意の回数繰り返されてもよい。ユーザは削除すべき系列が無いならば完了ボタン１１０５を押下する。ユーザ操作取得部２０５は、完了ボタン１１０５が押下された場合、処理をＳ１００４に進める。
次にＳ１００４において、Ｓ４０５～Ｓ４１１と同様の処理が実行される。本実施形態では、Ｓ１００３で削除指示された系列に対応する基準点候補が、候補対象から除外される。以上のようにして、図１０に示すフローチャートの一連の処理が終了する。
次にＳ１００４において、Ｓ４０５～Ｓ４１１と同様の処理が実行される。本実施形態では、Ｓ１００３で削除指示された系列に対応する基準点候補が、候補対象から除外される。以上のようにして、図１０に示すフローチャートの一連の処理が終了する。
【0063】
なお、より利便性を高めるため、ユーザ操作取得部２０５は、系列を新たに追加する機能を備えていてもよい。追加ボタン１１０７が押下されると、ユーザ操作取得部２０５は、過去に撮像、編集した学習データを、入力欄１１０６で指定された格納場所から読み出して、現在編集中の学習データに追加する。また、再読み込みボタン１１０８が押下されると、ユーザ操作取得部２０５は、現在編集中の学習データの編集を終了して、入力欄１１０６で指定された格納場所から学習データが新たに読み出す。表示制御部２０４は、新たに読み出した学習データに基づき学習データ形状１１０２を表示し直す。
【0064】
なお、データ保持部２０７は、ソース色空間の学習データ（関数形式）の形状を色票３のパッチ（色）ごとに解析し、解析の結果、色変換精度低下の要因となり得る所定の条件を満たすと判定した場合に、そのパッチ（色）に対応する学習データを削除してもよい。色変換精度低下の要因となり得る所定の条件の例としては、他の系列と交差するなど撮影誤差が大きいと判断できる条件がある。具体的には、データ保持部２０７は、パッチ（色）に対応する系列が、他の系列と所定値以上交差するか否かを判定し、判定の結果、所定値以上交差すると判定した場合に、当該系列を削除してもよい。
【0065】
以上説明したように、実施形態２によれば、撮像誤差の大きい学習データを予め除去することが可能になり、色変換の誤差を低減することができる。また、学習データの形状を表示することで、誤差の大きいものを除去する作業を支援することができる。
【0066】
［その他の実施形態］
なお、各種のデバイスと接続され、各種のデバイスを動作させる装置あるいはシステム内のコンピュータに、上記の画像処理装置２の機能を実現するためのソフトウェアのプログラムコードを供給しておいても構わない。この場合、そのシステムあるいは装置のコンピュータ（ＣＰＵあるいはＭＰＵ）が、供給されたプログラムに従い、各種のデバイスを動作させることによって、上記の実施形態を実施する場合も本発明の範疇に含まれる。例えば、上記の画像処理をコンピュータにインストールされたドライバソフト上で行う場合が含まれる。その場合、前記のソフトウェアのプログラムコード自体が前述した実施形態の機能を実現することになり、そのプログラムコード自体、および、そのプログラムコードをコンピュータに供給するための、当該プログラムコードが格納された記憶媒体は本発明を構成する。
なお、各種のデバイスと接続され、各種のデバイスを動作させる装置あるいはシステム内のコンピュータに、上記の画像処理装置２の機能を実現するためのソフトウェアのプログラムコードを供給しておいても構わない。この場合、そのシステムあるいは装置のコンピュータ（ＣＰＵあるいはＭＰＵ）が、供給されたプログラムに従い、各種のデバイスを動作させることによって、上記の実施形態を実施する場合も本発明の範疇に含まれる。例えば、上記の画像処理をコンピュータにインストールされたドライバソフト上で行う場合が含まれる。その場合、前記のソフトウェアのプログラムコード自体が前述した実施形態の機能を実現することになり、そのプログラムコード自体、および、そのプログラムコードをコンピュータに供給するための、当該プログラムコードが格納された記憶媒体は本発明を構成する。
【0067】
上述の各実施形態の開示は、以下の構成、方法及びプログラムを含む。
（構成１）
複数の色票を撮像した画像において計測された、前記複数の色票の色を表す第１のデータと、前記第１のデータが計測された計測点における目標の色を表す第２のデータを取得するデータ取得手段と、
色空間上の入力点を取得する入力点取得手段と、
前記第１のデータから、前記複数の色票に対応する複数の基準点候補を取得し、複数の前記基準点候補の中から、前記入力点との色空間上の距離が近い複数の前記基準点候補を、基準点として選択する基準点選択手段と、
前記入力点の位置に基づいて、色変換の係数を取得する係数取得手段と、
前記色変換の係数と、複数の前記基準点に対応する前記第２のデータに基づいて、前記入力点における色変換結果を算出する算出手段と、
を有することを特徴とする画像処理装置。
（構成２）
前記係数取得手段は、複数の前記基準点と、前記入力点との相対的な位置に基づいて、前記色変換の係数を算出することを特徴とする構成１に記載の画像処理装置。
（構成３）
前記算出手段は、前記計測点における、前記第１のデータと前記第２のデータとの変分値を算出し、前記色変換の係数と、複数の前記基準点に対応する前記変分値に基づいて、前記入力点における色変換結果を算出することを特徴とする構成１又は２に記載の画像処理装置。
（構成４）
前記画像は、撮像装置において撮像された画像であって、
前記撮像装置が３次元ルックアップテーブルで保持する色変換関数を生成する生成手段を更に有し、
前記生成手段は、前記３次元ルックアップテーブルの格子点に対応する前記入力点と、当該入力点における前記色変換結果とに基づいて、前記色変換関数を生成することを特徴とする構成１～３の何れか１つに記載の画像処理装置。
（構成５）
前記データ取得手段は、同一の色票から輝度の異なる複数の計測値を取得し、
前記複数の計測値の間の値を前記複数の色票ごとに補間する補間手段、を更に有し、
前記基準点選択手段は、前記補間手段により補間された前記第１のデータから、前記複数の色票に対応する複数の基準点候補を取得することを特徴とする構成１～４の何れか１つに記載の画像処理装置。
（構成６）
前記補間手段は、前記複数の計測値を輝度で表現する関数を前記複数の色票ごとに生成することを特徴とする構成５に記載の画像処理装置。
（構成７）
前記補間手段は、前記複数の計測値を多項式近似することにより前記関数を生成することを特徴とする構成６に記載の画像処理装置。
（構成８）
前記画像は、露光時間を多段階にして前記複数の色票を撮像した画像データ群、開口率を多段階にして前記複数の色票を撮像した画像データ群、及び照明強度を多段階にして前記複数の色票を撮像した画像データ群のうちの少なくとも何れか１つを含むことを特徴とする構成５～７の何れか１つに記載の画像処理装置。
（構成９）
前記基準点選択手段は、前記複数の色票ごとの前記関数に対して、前記入力点の輝度値を代入して得られる点を、複数の前記基準点候補として取得することを特徴とする構成６～７の何れか１つに記載の画像処理装置。
（構成１０）
前記基準点選択手段は、複数の前記基準点候補の中から、前記入力点との色相の近い複数の前記基準点候補を、基準点として選択することを特徴とする構成１～９の何れか１つに記載の画像処理装置。
（構成１１）
前記データ取得手段は、同一の色票から輝度の異なる複数の計測値を取得し、
前記複数の計測値の間の値を前記複数の色票ごとに補間する補間手段と、
前記補間手段により補間された前記第１のデータを表す形状データを、前記複数の色票ごとに表示するよう制御する表示制御手段と、
を更に有することを特徴する構成１～１０の何れか１つに記載の画像処理装置。
（構成１２）
前記基準点選択手段は、前記複数の色票に対応する前記第１のデータのうち、所定の条件を満たす一部の色票に対応する前記第１のデータを、複数の基準点候補を取得する対象から除外することを特徴する構成１～１１の何れか１つに記載の画像処理装置。
（構成１３）
前記基準点選択手段は、前記複数の色票に対応する前記第１のデータのうち、ユーザより操作指示のあった色票に対応する前記第１のデータを、複数の基準点候補を取得する対象から除外することを特徴とする構成１２に記載の画像処理装置。
（方法）
複数の色票を撮像した画像において計測された、前記複数の色票の色を表す第１のデータと、前記第１のデータが計測された計測点における目標の色を表す第２のデータを取得するデータ取得工程と、
色空間上の入力点を取得する入力点取得工程と、
前記第１のデータから、前記複数の色票に対応する複数の基準点候補を取得し、複数の前記基準点候補の中から、前記入力点との色空間上の距離が近い複数の前記基準点候補を、基準点として選択する基準点選択工程と、
前記入力点の位置に基づいて、色変換の係数を取得する係数取得工程と、
前記色変換の係数と、複数の前記基準点に対応する前記第２のデータに基づいて、前記入力点における色変換結果を算出する算出工程と、
を含むことを特徴とする画像処理方法。
（プログラム）
画像処理装置のコンピュータを、
複数の色票を撮像した画像において計測された、前記複数の色票の色を表す第１のデータと、前記第１のデータが計測された計測点における目標の色を表す第２のデータを取得するデータ取得手段と、
色空間上の入力点を取得する入力点取得手段と、
前記第１のデータから、前記複数の色票に対応する複数の基準点候補を取得し、複数の前記基準点候補の中から、前記入力点との色空間上の距離が近い複数の前記基準点候補を、基準点として選択する基準点選択手段と、
前記入力点の位置に基づいて、色変換の係数を取得する係数取得手段と、
前記色変換の係数と、複数の前記基準点に対応する前記第２のデータに基づいて、前記入力点における色変換結果を算出する算出手段、
として機能させるためのプログラム。
（構成１）
複数の色票を撮像した画像において計測された、前記複数の色票の色を表す第１のデータと、前記第１のデータが計測された計測点における目標の色を表す第２のデータを取得するデータ取得手段と、
色空間上の入力点を取得する入力点取得手段と、
前記第１のデータから、前記複数の色票に対応する複数の基準点候補を取得し、複数の前記基準点候補の中から、前記入力点との色空間上の距離が近い複数の前記基準点候補を、基準点として選択する基準点選択手段と、
前記入力点の位置に基づいて、色変換の係数を取得する係数取得手段と、
前記色変換の係数と、複数の前記基準点に対応する前記第２のデータに基づいて、前記入力点における色変換結果を算出する算出手段と、
を有することを特徴とする画像処理装置。
（構成２）
前記係数取得手段は、複数の前記基準点と、前記入力点との相対的な位置に基づいて、前記色変換の係数を算出することを特徴とする構成１に記載の画像処理装置。
（構成３）
前記算出手段は、前記計測点における、前記第１のデータと前記第２のデータとの変分値を算出し、前記色変換の係数と、複数の前記基準点に対応する前記変分値に基づいて、前記入力点における色変換結果を算出することを特徴とする構成１又は２に記載の画像処理装置。
（構成４）
前記画像は、撮像装置において撮像された画像であって、
前記撮像装置が３次元ルックアップテーブルで保持する色変換関数を生成する生成手段を更に有し、
前記生成手段は、前記３次元ルックアップテーブルの格子点に対応する前記入力点と、当該入力点における前記色変換結果とに基づいて、前記色変換関数を生成することを特徴とする構成１～３の何れか１つに記載の画像処理装置。
（構成５）
前記データ取得手段は、同一の色票から輝度の異なる複数の計測値を取得し、
前記複数の計測値の間の値を前記複数の色票ごとに補間する補間手段、を更に有し、
前記基準点選択手段は、前記補間手段により補間された前記第１のデータから、前記複数の色票に対応する複数の基準点候補を取得することを特徴とする構成１～４の何れか１つに記載の画像処理装置。
（構成６）
前記補間手段は、前記複数の計測値を輝度で表現する関数を前記複数の色票ごとに生成することを特徴とする構成５に記載の画像処理装置。
（構成７）
前記補間手段は、前記複数の計測値を多項式近似することにより前記関数を生成することを特徴とする構成６に記載の画像処理装置。
（構成８）
前記画像は、露光時間を多段階にして前記複数の色票を撮像した画像データ群、開口率を多段階にして前記複数の色票を撮像した画像データ群、及び照明強度を多段階にして前記複数の色票を撮像した画像データ群のうちの少なくとも何れか１つを含むことを特徴とする構成５～７の何れか１つに記載の画像処理装置。
（構成９）
前記基準点選択手段は、前記複数の色票ごとの前記関数に対して、前記入力点の輝度値を代入して得られる点を、複数の前記基準点候補として取得することを特徴とする構成６～７の何れか１つに記載の画像処理装置。
（構成１０）
前記基準点選択手段は、複数の前記基準点候補の中から、前記入力点との色相の近い複数の前記基準点候補を、基準点として選択することを特徴とする構成１～９の何れか１つに記載の画像処理装置。
（構成１１）
前記データ取得手段は、同一の色票から輝度の異なる複数の計測値を取得し、
前記複数の計測値の間の値を前記複数の色票ごとに補間する補間手段と、
前記補間手段により補間された前記第１のデータを表す形状データを、前記複数の色票ごとに表示するよう制御する表示制御手段と、
を更に有することを特徴する構成１～１０の何れか１つに記載の画像処理装置。
（構成１２）
前記基準点選択手段は、前記複数の色票に対応する前記第１のデータのうち、所定の条件を満たす一部の色票に対応する前記第１のデータを、複数の基準点候補を取得する対象から除外することを特徴する構成１～１１の何れか１つに記載の画像処理装置。
（構成１３）
前記基準点選択手段は、前記複数の色票に対応する前記第１のデータのうち、ユーザより操作指示のあった色票に対応する前記第１のデータを、複数の基準点候補を取得する対象から除外することを特徴とする構成１２に記載の画像処理装置。
（方法）
複数の色票を撮像した画像において計測された、前記複数の色票の色を表す第１のデータと、前記第１のデータが計測された計測点における目標の色を表す第２のデータを取得するデータ取得工程と、
色空間上の入力点を取得する入力点取得工程と、
前記第１のデータから、前記複数の色票に対応する複数の基準点候補を取得し、複数の前記基準点候補の中から、前記入力点との色空間上の距離が近い複数の前記基準点候補を、基準点として選択する基準点選択工程と、
前記入力点の位置に基づいて、色変換の係数を取得する係数取得工程と、
前記色変換の係数と、複数の前記基準点に対応する前記第２のデータに基づいて、前記入力点における色変換結果を算出する算出工程と、
を含むことを特徴とする画像処理方法。
（プログラム）
画像処理装置のコンピュータを、
複数の色票を撮像した画像において計測された、前記複数の色票の色を表す第１のデータと、前記第１のデータが計測された計測点における目標の色を表す第２のデータを取得するデータ取得手段と、
色空間上の入力点を取得する入力点取得手段と、
前記第１のデータから、前記複数の色票に対応する複数の基準点候補を取得し、複数の前記基準点候補の中から、前記入力点との色空間上の距離が近い複数の前記基準点候補を、基準点として選択する基準点選択手段と、
前記入力点の位置に基づいて、色変換の係数を取得する係数取得手段と、
前記色変換の係数と、複数の前記基準点に対応する前記第２のデータに基づいて、前記入力点における色変換結果を算出する算出手段、
として機能させるためのプログラム。
【符号の説明】
【0068】
１：デジタルカメラ、２：画像処理装置、３：色票
【図1A】
【図1B】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】