(57)【要約】【課題】  従来技術として前述のパラメータ初期値に依存しないモンテカルロ法やＳＡ等の最適化手法は、計算時間が膨大になり実時間でのパラメータ算出が困難である。【解決手段】  従来技術にあるパラメータ値更新でのランダム性を排除した最適化手法を実現し、これによって効率的な色補正用パラメータを算出するものである。【選択図】  図１
2006339776.tif 000002
【特許請求の範囲】
【請求項１】
画像を出力する画像出力手段と、前記画像出力装置内での変換特性を算出する変換特性算出手段とからなる画像処理装置であって、
前記画像出力手段は、
第１の表色系で表現される色信号値を、同一または異なる第３の表色系の色信号値に変換する色信号変換手段と、
前記変換された第３の表色系の色信号値を、さらに第２の表色系で表現される色信号値へ変換する色信号変換手段と、
前記変換された第２の表色系の色信号値に基づいて画像を出力する画像出力手段
を有すると共に、
さらに変換特性算出装置で算出された変換特性に基づき、前記第１の色信号値に対する前記第２の色信号値の変換値を補正する補正手段
を有し、
前記変換特性算出手段は、変換特性の算出を
前記第３の表色系の色信号に色信号補正を施して画像出力した際の色を、第４の表色系で推定する補正色推定手段と、
前記第３の表色系の色信号に対しての色再現目標を第４の表色系で与える色再現目標設定手段、
を用い、
前記算出した変換特性と色再現目標とから新たに変換特性を算出する変換特性算出手段
を用いて前記推定値が目標値と一致するよう漸近的に変換特性を変更していくと共に
変更の途中段階に於いて、色再現目標の少なくとも一つを再現するよう、前記変換特性を修正する
ことを特徴とする画像処理装置。
【請求項２】
請求項１に記載の画像処理装置であって、前記画像出力手段における補正手段は
前記変換された第３の表色系の色信号値を、算出された変換特性に従って同じく第３の表色系の色信号値に変換する色信号変換手段により構成される
ことを特徴とする画像処理装置。
【請求項３】
請求項１に記載の画像処理装置であって、前記画像出力手段における補正手段は
第１の表色系の色信号値から同第３の表色系の色信号値への変換に、算出された変換特性に基き第３の表色系の色信号値から同一の第３の表色系への変換を重畳し、改めて第１の表色系の色信号値から同第３の表色系の色信号値への変換として設定する手段か、
算出された変換特性に基い第３の表色系の色信号値から同一の第３の表色系への変換に、第３の表色系の色信号値から同第２の表色系の色信号値への変換を重畳し、改めて第２の表色系の色信号値から同第４の表色系の色信号値への変換として設定する手段
のいずれかで構成される
ことを特徴とする画像処理装置。
【請求項４】
請求項１〜３に記載の画像処理装置であって、前記変換特性算出手段において前記推定値が目標値に近くなるよう変換特性を算出する際は、
前記目標色再現設定手段により設定された前記複数色の夫々につき、第４の表色系での推定色と目標色との差分を評価値として算出し、前記評価値の総和が減ずるよう変換特性を更新していくことで、変換特性を算出する
と共に、
前記変換特性に対し、前記目標色再現設定手段により設定された前記複数色の入力値により構成される前記第一の行列と、第４の表色系で推定する目標色再現と同一の色再現を実現する前記複数色の入力値および前記変換特性により変換された入力値により構成される前記第２の行列による行列演算から変換特性を算出する
と共に、
前記算出した変換特性と第４の表色系で推定する色再現目標とから新たに変換特性を算出する変換特性算出手段を用いて前記推定値が目標値と一致するよう変換特性を算出する
ことを特徴とする画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、色再現を管理する為の画像処理装置に関する。
【背景技術】
【０００２】
近年、パーソナルコンピュータ／ワークステーションの普及に伴い、これらを用いたＤＴＰ（デスクトップ・パブリッシング）やＣＡＤ（コンピュータ・エイディッド・デザイン）によるカラードキュメント作成が広く一般に用いられるようになってきた。このような中、種々のカラープリンタが混在する中でカラー印刷の色再現を一定に保つため、高精度かつ簡易な色再現管理技術が求められている。
【０００３】
まず高精度な色再現管理という意味では、印字デバイスの如何に依らない安定した一定の色再現、ならびに印字デバイスの経時変化あるいは環境変動に依らない安定した色再現が求められている。
【０００４】
他方で簡易な色再現管理技術という意味では、色再現管理に関する初心者でも使いこなせる色再現管理技術が求められている。安価なカラー印刷機器が普及してきた現在では、オフィス等でカラーカタログ印刷を行う例が増えつつある。そこで、従来色再現管理に携わってこなかったユーザでも色再現を簡易に管理できるような色再現管理技術が求められている。
【０００５】
このような高精度かつ簡易な色再現管理技術として、多項式近似を用いた画像処理装置が知られており、また出願人はこの画像処理装置に関連して色管理を行うべき色を指定し、その色の再現色を与えると共に重要度（より正確には重み）を指定するだけで、色補正パラメータを算出し、色再現管理を行う装置としての画像処理装置を出願している。
【０００６】
上記画像処理装置の色補正パラメータ算出は、線形計画法あるいは非線形計画法に代表される最適化手法によって行われている。一般に最適化手法は初期設定（パラメータ初期値、スケジューリング等）により収束先、収束までの挙動が異なることが知られている。例えば、ある初期値では数千回の探索を行っても目標収束値に到達し得なかったとしても、異なる初期値を与えた場合、数回の探索で目標収束値に達することができることがある。この様に、初期設定は探索効率及び探索による収束先に大きく影響する。しかしながら、これら初期設定は問題の特徴と目的から設計者の経験則により決定されることが一般的であり、当然ながら対象となる各最適化問題に対して異なる設定が必要である。これら初期設定のうちパラメータ初期値に依存しない最適化手法としてモンテカルロ法、ＳＡ（シミュレーテッドアニーリング）等が知られている。上記手法は任意のパラメータ初期値から最適値探索を行い、ある程度局所解に収束した時点でランダムに最適化パラメータを変更し、変更された最適化パラメータを初期値に置き換えて再度探索を継続する手法である。ただし、これら手法は最適化パラメータの変更をランダムに行うため、探索回数が増加することによる探索効率悪化は避けられないばかりか、必ずしも最適値に収束する保証もない。
【０００７】
又、従来例としては、例えば特許文献１をあげることが出来る。
【特許文献１】特開２００５−００５８０２号公報
【発明の開示】
【発明が解決しようとする課題】
【０００８】
従来技術として前述のパラメータ初期値に依存しないモンテカルロ法やＳＡ等の最適化手法は、計算時間が膨大になり実時間でのパラメータ算出が困難である。
【０００９】
そこで、本発明は、従来技術にあるパラメータ値更新でのランダム性を排除した最適化手法を実現し、これによって効率的な色補正用パラメータを算出するものである。
【課題を解決するための手段】
【００１０】
以上のような課題を解決する為、本発明は
画像を出力する画像出力手段と、前記画像出力装置内での変換特性を算出する変換特性算出手段とからなる画像処理装置であって、
前記画像出力手段は、
第１の表色系で表現される色信号値を、同一または異なる第３の表色系の色信号値に変換する色信号変換手段と、
前記変換された第３の表色系の色信号値を、さらに第２の表色系で表現される色信号値へ変換する色信号変換手段と、
前記変換された第２の表色系の色信号値に基づいて画像を出力する画像出力手段
を有すると共に、
さらに変換特性算出装置で算出された変換特性に基づき、前記第１の色信号値に対する前記第２の色信号値の変換値を補正する補正手段
を有し、
前記変換特性算出手段は、変換特性の算出を
前記第３の表色系の色信号に色信号補正を施して画像出力した際の色を、第４の表色系で推定する補正色推定手段と、
前記第３の表色系の色信号に対しての色再現目標を第４の表色系で与える色再現目標設定手段、
を用い、
前記算出した変換特性と色再現目標とから新たに変換特性を算出する変換特性算出手段
を用いて前記推定値が目標値と一致するよう漸近的に変換特性を変更していくと共に
変更の途中段階に於いて、色再現目標の少なくとも一つを再現するよう、前記変換特性を修正する
ことを特徴とする画像処理装置とを用い、前記推定値が目標値に近くなるよう変換特性を算出する構成をとる。
【発明の効果】
【００１１】
以上の説明から明らかなように本発明に依れば、色補正パラメータを最適化手法により求める際に、パラメータ値更新におけるランダム性を排除した最適化手法により、短時間かつ効率的に、色補正用パラメータの算出を可能にするものである。
【発明を実施するための最良の形態】
【００１２】
＜第１実施例＞
図１は本発明の第１の実施例としての画像処理装置を示すブロック図である。前記構成において、１０１はＣＰＵ、１０２はメインメモリ、１０３はＳＣＳＩインタフェース、１０４はネットワークインタフェース、１０５はＨＤＤ、１０６はグラフィックアクセラレータ、１０７はカラーモニタ、１０８はＵＳＢコントローラ、１０９はカラープリンタ、１１０はキーボード／マウスコントローラ、１１１はキーボード、１１２はマウス、１１３はローカルエリアネットワーク、１１４はＰＣＩバスである。
【００１３】
上記構成における、ディジタル画像のプリンタ出力動作について述べる。まず、ＨＤＤ１０５に格納されている画像アプリケーションが、ユーザの指示を受けたＯＳプログラムに基づき、ＣＰＵ１０１にて起動される。続いてユーザの指示による画像アプリケーション内の処理にしたがって、ＨＤＤ１０５に格納されているディジタル画像データが、ＣＰＵ１０１からの指令に基づきＳＣＳＩＩ／Ｆ１０３を介してＰＣＩバス１１４経由によりメインメモリ１０２に転送される。あるいは、ＬＡＮに接続されているサーバに格納されている画像データあるいはインターネット上のディジタル画像データが、ＣＰＵ１０１からの指令によりネットワークＩ／Ｆ１０４を介してＰＣＩバス１１４経由によりメインメモリ１０２に転送される。以下、前記メインメモリ１０２に保持されているディジタル画像データは、ＲＧＢ各色信号が符号無し８ｂｉｔで表現される画像データであるものとして実施例を説明する。前記メインメモリ１０２に保持されているディジタル画像データは、ＣＰＵ１０１からの指令によりＰＣＩバス１１４経由によってグラフィックアクセラレータ１０６に転送され、グラフィックアクセラレータ１０６はディジタル画像データをＤ／Ａ変換した後ディスプレイケーブルを通じてカラーモニタ１０７に送信する。これにより、カラーモニタ１０７上に画像が表示される。ここで、ユーザがメインメモリ１０２に保持されているディジタル画像をプリンタ１０９から出力するよう画像アプリケーションに指令すると、画像アプリケーションはディジタル画像データをプリンタドライバに転送する。ＣＰＵ１０１は、プリンタドライバの後述の処理フローチャートに基づき、ディジタル画像データをＣＭＹＫディジタル画像データに変換し、ＵＳＢコントローラ１０８を介して前記ＣＭＹＫ画像データをプリンタ１０９へ送信する。以上一連の動作の結果として、プリンタ１０９よりＣＭＹＫ画像が印字される。
【００１４】
以下では、上記構成におけるプリンタドライバ動作について、図２のフローチャートを用いて説明する。本実施例におけるプリンタドライバは、ユーザの好みに応じて色を変換する為の色補正ＬＵＴを用いた色変換、所望の色再現を得る為にＲＧＢ値を修正する為の色補正変換、ＲＧＢ色信号をＣＭＹＫ色信号に変換する為の色変換ＬＵＴを用いた変換、以上の３つの色変換を用いてＲＧＢ２４ｂｉｔ画像データをＣＭＹＫ３２ｂｉｔ画像データに変換する。
【００１５】
まず、ステップ２０１では、メインメモリ１０２に保持されているＲＧＢ２４ｂｉｔ画像データ（各色信号が符号無し８ｂｉｔの画像データ）へのポインタを取得すると共にＣＭＹＫ変換後の画像データを保持する為のメモリ領域を確保し、同時に他の初期化動作を行う。
【００１６】
次にステップ２０２では、色補正で用いるパラメータを、ＨＤＤ１０５より取得する。色補正パラメータと色補正変換の詳細に関しては、ステップ２０５の説明にて述べる。
【００１７】
ステップ２０３では、ラスタスキャンの順に従って、メインメモリ１０２内のＲＧＢ２４ｂｉｔ画像データから画素値であるＲＧＢ値を１つ取得する。
【００１８】
ステップ２０４では、色補正ＬＵＴを用いての四面体補間演算により、ＲＧＢ２４ｂｉｔ→ＲＧＢ２４ｂｉｔ変換を施す。色補正ＬＵＴはＲＧＢ色空間での格子点の色座標データと、前記格子点が再現するＲＧＢ色空間の座標値との対応を記した、図１７の様なデータ構造で構成される。データ構造の先頭には、Ｒ／Ｇ／Ｂ値のステップが記述され、この後に、各格子点に対応するＲＧＢ色座標が、Ｒ、Ｇ、Ｂの順でネストされて記述される。このデータ構造をＲＧＢ色空間で模式的に表すと、図３の様に表される。また、四面体補間は、隣接する８つの格子点で囲まれる６面体領域を、図４の模式図に示す様な６つに分割された四面体領域毎に、補間演算を行う。
【００１９】
ステップ２０５では、３次元ベクトルに対する２次非線形変換と各要素独立でのγ変換とを組み合わせ、図５の式の様に色補正変換を行う。ここで、図５の行列要素ｍij、γ変換係数γiが色補正パラメータである。演算は浮動小数点演算により実施されるが、図５の式による変換が終わった後は、各色信号が符号無し８ｂｉｔデータへと丸められる。また、この際２５５を超える値は２５５、０を下回る値は０へとクリップされる。
【００２０】
ステップ２０６では、色変換ＬＵＴに従った四面体補間により、ＲＧＢ２４ｂｉｔデータをＣＭＹＫ３２ｂｉｔ画像データ（ＣＭＹＫの各色信号が符号無し８ｂｉｔの画像データ）に変換する。ここで、以下の説明ではこの変換前のＲＧＢ色空間を、特にＰｒｉｖａｔｅ  ＲＧＢと呼ぶ。つまり、ステップ２０４でのＲＧＢ２４ｂｉｔ→ＲＧＢ２４ｂｉｔ変換以降ステップ２０５での色補正演算は、Ｐｒｉｖａｔｅ  ＲＧＢ上で演算である。
【００２１】
ステップ２０７では、このＣＭＹＫ３２ｂｉｔ値を画素位置に従ったしかるべきメモリアドレスに書き込む。
【００２２】
ステップ２０８では、必要な画素総てに対し変換を施したかどうかを判断し、施していればステップ２０９へ、施していなければステップ２０３へジャンプする。
【００２３】
ステップ２０９では、メインメモリに格納されたＣＭＹＫ３２ｂｉｔ画像データを、ＵＳＢコントローラ１０８を介してプリンタ１０９へ送信し、ステップ２１０で動作を終了する。
【００２４】
以下では図１の構成における、色補正パラメータ作成アプリケーションを用いたパラメータ作成動作について説明する。
【００２５】
最初、ＨＤＤ１０５に格納されている色補正パラメータ作成アプリケーションが、ユーザの指示を受けたＯＳプログラムに基づき、ＣＰＵ１０１にて起動される。色補正パラメータ作成アプリケーションが起動されると、図６に示すダイアログウィンドウが表示される。ダイアログ並びにアプリケーションの動作フローについては後で詳しく説明するが、ユーザは図６のウィンドウを用いて、Ｐｒｉｖａｔｅ  ＲＧＢ値（ステップ２０６の説明参照）と、そのＲＧＢ値が再現すべき色のＬ*ａ*ｂ*値とについて確認ならびに指定する。ここで、Ｐｒｉｖａｔｅ  ＲＧＢ値と、このＲＧＢ値で画像出力した際の色再現目標たるＬ*ａ*ｂ*値との対応セット（以下、目標色再現情報と呼ぶ）に問題ないとユーザが判断したユーザがパラメータ作成ボタンを押下することで、色補正パラメータが生成され、ＨＤＤ１０５に保存される。
【００２６】
以下では、図６のダイアログウィンドウについて説明する。６０１は目標色再現情報をリスト表示するウィンドウであり、ウィンドウ内の左側に対象となっている色のＰｒｉｖａｔｅ  ＲＧＢ値が、右側に画像出力した際の色再現目標たるＬ*ａ*ｂ*値が表示される。また、選択されている目標色再現情報が反転表示される。６０２はスライダバーであり、このバーを制御することで表示されている目標色再現情報をスクロールならびに選択する。６０３は再現色指定ボタンであり、このボタンを押下すると図７の設定ダイアログが表示され、この図７のダイアログから、ユーザは目標色再現情報を設定する。６０４はパラメータ作成ボタンであり、このボタンを押下すると、色補正パラメータ作成アプリケーションは色補正パラメータを作成し、ＨＤＤ１０５に保存する。６０５は色再現保証チェックボックスであり、当該チェックボックスにチェックすると色補正パラメータ作成時にチェックのある目標色再現情報を特に重要色としてパラメータが作成される。６０６は終了ボタンであり、このボタンを押下すると色補正パラメータ作成アプリケーションを終了する。
【００２７】
以下では、色補正パラメータ作成アプリケーション動作について、図８の状態図を用いて説明する。
【００２８】
ステート８０１では、目標色再現情報の初期設定値を読み込むといった初期化動作を行う。
【００２９】
次にステート８０２では、図６のウィンドウでのユーザ操作判断待ち状態となる。ここで、６０２のスライダバーが操作されるとステート８０３へ移行し、再現色指定ボタン６０３が押下されるとステート８０４へ移行し、パラメータ作成ボタン６０４が押下されるとステート８０５へ移行し、終了ボタン６０５が押下されるとステート８０６へ移行する。
【００３０】
ステート８０３では、スライダバーの制御量に応じて、目標色再現情報のウィンドウ６０１への表示をスクロールすると共に、選択されている目標色再現情報を変更する。
【００３１】
ステート８０４では、図６のダイアログを表示して、目標色再現情報の編集状態に入る。ここで、７０１のエディットボックスで、編集対象となっている色のＰｒｉｖａｔｅ  ＲＧＢ値を変更し、７０２のエディットボックスで色再現目標たるＬ*ａ*ｂ*値を編集する。また、７０３のエディットボックスでは、編集対象となっている色の重要度を整数で設定する。ただし、６０５のチェックボックスにチェックがあるなら当該エディットボックスは入力不能となる。この値は、後述の評価値算出処理のステップ１００８にて、評価値算出の際の重みとして用いる。７０４の修正ボタンを押下すると、現在選択されている目標色再現情報に上書きし、７０５の追加ボタンを押下すると、新たな目標色再現情報として追加する。７０６の削除ボタンを押下すると、現在選択されている目標色再現情報を削除する。また、これら７０４、７０５、７０６、ボタンを押すと、それぞれの動作を行った後に目標色再現情報に記憶し、図６のダイアログをクローズする。
【００３２】
ステート８０５では、次に説明する図９のフローチャートに従って色補正パラメータを作成し他の後にＨＤＤ１０５に保存する。
【００３３】
ステート８０６では、メモリ開放などの終了動作を行った後、色補正パラメータ作成アプリケーション動作を終了する。
【００３４】
以下では、ステート８０５における色補正パラメータ作成処理ついて、図９のフローチャートを用いて説明する。
【００３５】
ステップ９０１では、色補正パラメータ（図５に示した式の行列要素ｍij、γ変換係数γi）の初期設定やメモリ領域の確保といった初期化を行う。また、イタレーション回数を表す変数ｉに０を、補正回数を表す変数ｊに０を設定する。
【００３６】
ステップ９０２では、まず現在の色補正パラメータに対する評価値ｓｕｍを算出する。評価値ｓｕｍの算出に関しては、この後で図１０のフローチャートを用いて説明する。
【００３７】
ステップ９０３では、算出した評価値ｓｕｍが所定値より小さいかを判断し、所定値より小さければステップ９０８へ、所定値より大きければステップ９０４へ進む。
【００３８】
ステップ９０４では、イタレーション回数を表すｉが所定の回数以上かどうかを判定し、所定回数未満であればステップ９０５へ、所定回数以上であればステップ９０６へ進む。
【００３９】
ステップ９０５では、後述の色補正パラメータ算出アルゴリズムに基づいて色補正パラメータを更新する。
【００４０】
ステップ９０６では、後述の色補正パラメータ補正アルゴリズムに基づいて色補正パラメータを変更する。
【００４１】
ステップ９０７では、色補正パラメータの補正回数を表すｊが所定の回数以上かどうかを判定し、所定回数未満であればステップ９０２へ、所定回数以上であればステップ９０８へ進む。
【００４２】
ステップ９０８では、現在の色補正パラメータをメモリに記憶して、色補正パラメータ作成処理を終了する。
【００４３】
以下では、ステップ９０２で用いた評価値算出処理について、図１０のフローチャートを用いて説明する。このフローでは、ｉ番目の色再現目標色に対する評価値ｆiと、ｆiの総和による評価値ｓｕｍと算出する。
【００４４】
ステップ１００１では、まず評価すべき色補正パラメータを取得し、設定する。
【００４５】
ステップ１００２では、初期化動作としてループ回数を示す変数であるｉを１に設定する。また、ステップ１００６での色推定に必要な色推定用ＬＵＴを読み込む。当該ＬＵＴはステップ２０４の説明内で述べた、図３ならびに図１７のデータ構造と全く同じ形式である為、詳細に関しては割愛する。尚、このＬＵＴは、Ｐｒｉｖａｔｅ  ＲＧＢ値の色を実際に画像出力したものをＬ*ａ*ｂ*表色系で測定したものであり、あらかじめ準備してある。
【００４６】
ステップ１００３では、Ｐｒｉｖａｔｅ  ＲＧＢ値と、このＲＧＢ値で画像を出力した際の色再現目標たるＬ*ａ*ｂ*値との対応セットであるところのｉ番目の目標色再現情報を、メモリから読み込む。
【００４７】
ステップ１００４では、メモリから読み込んだ目標色再現情報より、Ｐｒｉｖａｔｅ  ＲＧＢ値  ＣｏｒｇとＬ*ａ*ｂ*値Ｄｒｅｐとの対応関係を一対取得する。また、当該色に対する重みｗも同時に取得する。
【００４８】
ステップ１００５では、ステップ２０５で用いたと同じ図５の式により、Ｐｒｉｖａｔｅ  ＲＧＢ値ＣｏｒｇをＣｃｍｐへ変換する。
【００４９】
ステップ１００７では、Ｐｒｉｖａｔｅ  ＲＧＢ値ＣｃｍｐよりＬ*ａ*ｂ*値Ｄｃａｌを、ステップ１００２で読み込んだ色推定用ＬＵＴを用いて推定する。推定演算には四面体補間を用い、図４の模式図に示す様な６つに分割された四面体領域毎に補間演算を行う。
【００５０】
ステップ１００８では、目標色再現値であるＬ*ａ*ｂ*値Ｄｒｅｐと、推定したＬ*ａ*ｂ*値Ｄｃａｌとの色差をＣＩＥΔＥ９４で算出し、重みｗを乗算して評価値ｆiとする。
【００５１】
ステップ１００９では、目標色再現情報内の総てのＰｒｉｖａｔｅ  ＲＧＢ値  ＣｏｒｇとＬ*a*b*値Ｄｒｅｐとの対応関係についてステップ１００４から１００８までの演算を施したかどうかを判断し、施したのであればステップ１０１０へ、施したのでなければループ回数を示す変数ｉに１を加算した後、１ステップ１００４へジャンプする。
【００５２】
ステップ１０１０では、評価値ｆiの総和を評価値ｓｕｍとして算出する。
【００５３】
ステップ１０１１では、評価値算出処理の終了動作を行う。
【００５４】
以下では、ステップ９０６の色補正パラメータ更新処理について図１１のフローチャートを用いて説明する。図１１のフローチャートは、評価値ｓｕｍを逐次に減少させる為、Ｄ．Ｌ．Ｓ．法（減衰最小自乗法）に基づく色補正パラメータ更新を行うものである。
【００５５】
ここで以下の説明では、行列要素ｍijとγ変換係数γiと一列に並べてベクトルｘとして扱い、色補正パラメータの各要素をあらためてベクトル要素ｘk（ｌ＝０〜ｌ）で表す。また、色補正パラメータｋｌ番目のパラメータｘkの更新すべき変化量をΔｘkとする。
【００５６】
ステップ１１０１では、更新すべき色補正パラメータ、これに基づく目標色再現情報総ての評価値、評価値空間と列ベクトルｘのｇｒａｄｉｅｎｔを計算する際に必要な微小変動幅、ダンピングファクタ（ρ）を取得し、設定する。ダンピングファクタは次に示す式１にて定義されるものとする。
ρ・Δx＝０                                               式１１−１
また微小変動幅の詳細に関しては、ステップ１１０６の説明にて述べる。
【００５７】
ステップ１１０２では、更新すべき色補正パラメータの各要素を次に示す式１１−２によってベクトル要素ｘiに変換する。
【００５８】
【数１】
2006339776.tif 000003
式１１−２
このときの列ベクトルｘをベクトルの初期値とする。また初期化動作としてループ回数を示す変数であるｓを１に設定する。
【００５９】
ステップ１１０３では、列ベクトルｘの第ｌ番目のパラメータｘｌをステップ１１０１にて設定した微小変動幅だけ加算し、変化させる。
【００６０】
ステップ１１０４では、ステップ１１０３にて変更した列ベクトルｘに基づき、図１０のフローチャートで説明した処理を用いて、評価値を算出する。
【００６１】
ステップ１１０５では、ステップ１１０３にて変更した列ベクトルｘをベクトルの初期値へ復元する。
【００６２】
ステップ１１０６では、ステップ１１０２にて設定した目標色再現情報総ての初期評価値とステップ１１０４にて算出した変動後評価値、ステップ１１０１にて設定した変動幅を用い、次に示す式１１−３から評価値空間と列ベクトルｘのｇｒａｄｉｅｎｔを計算する。
【００６３】
【数２】
2006339776.tif 000004
式１１−３
ｇｒａｄｉｅｎｔは目標色再現情報ｎ（ｎ＝０〜ｐ）に対して算出し、ａknとしてメモリに格納する。ここで微小変動幅は偏微分における変数の差分に相当する。
【００６４】
ステップ１１０７では、ステップ１１０３からステップ１１０６までの操作を列ベクトルｘの総てのパラメータに対して施したかを判断し、施したのであればステップ１２０８へ、施したのでなければ、ループ回数を示す変数ｉに１を加算した後、ステップ１１０３へジャンプする。
【００６５】
ステップ１１０８では、ステップ１１０１にて設定したダンピングファクタとステップおよび初期評価値、ステップ１１０６にて算出したｇｒａｄｉｅｎｔからＤ．Ｌ．Ｓ．法を適用し、次に示す式１１−４から列ベクトルｘのパラメータであるｘiの変動量Δｘiを算出する。
【００６６】
【数３】
2006339776.tif 000005
式１１−４
ステップ１１０９では、式１１−５、式１１−６から色補正パラメータを更新する。
ｘi＝ｘi＋Δｘi                                            式１１−５
ｍ＝｛ｘ1，ｘ2，…，ｘk｝，γ＝｛ｘk+1，ｘk+2，…，ｘl｝      式１１−６
ステップ１１１０では、色処理パラメータ更新処理の終了動作を行う。
【００６７】
以下では、ステップ９０７の色補正パラメータの変更について説明する。
【００６８】
当該処理は初期値に依らない最適化を行うための処理について図１２のフローチャートを用いて説明する。
【００６９】
ステップ１２０１では、総ての目標色再現情報、色推定用ＬＵＴ、を取得し設定する。
【００７０】
ステップ１２０２では、初期化動作としてループ回数を示す変数であるｉを１に設定する。
【００７１】
ステップ１２０３では、６０５のチェックボックスにチェックのある目標色再現情報を重要色情報として設定する。
【００７２】
ステップ１２０４では、ステップ１２０３にて設定された重要色とプライマリＲＧＢＣＭＹが、変換後にしかるべきＬ*ａ*ｂ*となるＰｒｉｖａｔｅＲＧＢ値をステップ１２０１にて設定された色推定用ＬＵＴから補間演算により算出する。
【００７３】
ステップ１２０５では、６０５のチェックボックスにチェックが入っていない目標色再現情報について、ステート８０５における色補正パラメータ作成処理にて作成された色補正パラメータを用いて変換を行う。
【００７４】
ステップ１２０６では、ステップ１２０５の操作を６０５のチェックボックスにチェックが入っていない総ての色再現目標色について施したかを判断し、施したのであれば１２０７へ、施したのでなければループ回数を示す変数ｉに１を加算した後、ステップ１２０５へジャンプする。
【００７５】
ステップ１２０７では、後述の色補正パラメータ変換処理により色補正パラメータを生成する。
【００７６】
ステップ１２０８では、色処理パラメータ更新処理の終了動作を行う。
【００７７】
以下では、ステップ１２０７の色補正パラメータの生成について図１３のフローチャートを用いて説明する。
【００７８】
ステップ１３０１では、ステップ１２０４にて取得された重要色とプライマリＲＧＢＣＭＹについて変換前ＰｒｉｖａｔｅＲＧＢ値および変換後ＰｒｉｖａｔｅＲＧＢ値を設定するとともに、ステップ１２０５にて変換された色再現目標色の変換後ＰｒｉｖａｔｅＲＧＢ値および変換前ＰｒｉｖａｔｅＲＧＢ値を設定する。
【００７９】
ステップ１３０２では、ステップ１３０１にて設定されたプライマリＲＧＢＣＭＹ、重要色および色再現目標色について、変換前ＰｒｉｖａｔｅＲＧＢ値を図１４に示す行列Ａの要素に設定された順にネストする。ここでは行列Ａは９×９行列とする。また変換前ＰｒｉｖａｔｅＲＧＢ値をＲ、Ｇ、Ｂにて表記する。
【００８０】
ステップ１３０３では、ステップ１３０１にて設定されたプライマリＲＧＢＣＭＹ、重要色および色再現目標色について、変換後ＰｒｉｖａｔｅＲＧＢ値を図１５に示す行列Ｂの要素に設定された順にネストする。ここでは行列Ｂは３×９行列とする。また変換後ＰｒｉｖａｔｅＲＧＢ値をＲ’、Ｇ’、Ｂ’にて表記する。
【００８１】
ステップ１３０４では、図１６の式に従い、色補正パラメータを算出する。ここでは色補正パラメータをＭ、行列ＡをＡ、行列ＢをＢにて表すものとする。また行列Ｍの要素［ｍij］は図５の色補正パラメータの各要素に対応するものである。
【００８２】
ステップ１３０５では、色補正パラメータ算出処理の終了動作を行う。
【００８３】
本実施例に依れば、簡便な操作で色再現目標を実現しうる色再現管理を行うことが可能である。
【００８４】
本実施例に依れば、ＳＡのランダム性を排除した最適化手法を用いて、パラメータを目標の方向へ変更させることが可能であり、これを用いて目標色再現を効率的に色補正パラメータを自動で作成することが可能である。
【００８５】
（他の実施の形態）
＜色再現目標設定＞
目標色再現情報は入力値と出力値の対応関係が明示されていれば良い。したがって、例えばステップ１００７にて用いた色推定用ＬＵＴを用いることも可能である。
【００８６】
＜初期値に依らない色補正パラメータ算出＞
前記実施例では初期値に依らない色補正パラメータ算出の際に用いる目標色情報は９色を用いたが９色以上の任意の色数を用いることが可能である。したがって、処理に用いる前期行列Ａはｍ×ｎ行列（ｍ＜８、ｎ＜８）、前期行列Ｂは３×ｍ行列を用いることが可能である。
【図面の簡単な説明】
【００８７】
【図１】本発明の第１の実施例としての画像処理装置のシステム構成を示すブロック図である。
【図２】本発明の第１の実施例におけるプリンタドライバでの画像出力動作を示すフローチャート図である。
【図３】ＲＧＢ色空間での格子点配置を表す模式図である。
【図４】四面体補間での領域分割を表す模式図である。
【図５】色補正の為の変換を表す式である。
【図６】色補正パラメータ作成アプリケーションのユーザインタフェースを示す図である。
【図７】再現色設定の為のユーザインタフェースを示す図である。
【図８】色補正パラメータ作成アプリケーションの動作を示す状態図である。
【図９】色補正パラメータ算出処理の動作を示すフローチャート図である。
【図１０】評価値算出の動作を示すフローチャート図である。
【図１１】ＤＬＳ法を用いた色補正パラメータ更新処理を示すフローチャート図である。
【図１２】作成した色補正パラメータに対する追加処理を示すフローチャート図である。
【図１３】最適化初期値に依らない色補正パラメータ算出処理を示すフローチャート図である。
【図１４】色補正パラメータ算出に用いる行列Ａの組成を示す模式図である。
【図１５】色補正パラメータ算出に用いる行列Ｂの組成を示す模式図である。
【図１６】色補正パラメータである行列Ｍを示す式である。
【図１７】ＬＵＴのデータ構造を表す図である。
【符号の説明】
【００８８】
１０１  ＣＰＵ
１０２  メインメモリ
１０３  ＳＣＳＩ  Ｉ／Ｆ
１０４  ネットワーク  Ｉ／Ｆ
１０５  ＨＤＤ
１０６  グラフィックアクセラレータ
１０７  カラーモニタ
１０８  ＵＳＢコントローラ
１０９  カラープリンタ
１１０  キーボード／マウスコントローラ
１１１  キーボード
１１２  マウス
１１３  ＬＡＮ
１１４  ＰＣＩバス
【図１】
2006339776.tif 000006
【図２】
2006339776.tif 000007
【図３】
2006339776.tif 000008
【図４】
2006339776.tif 000009
【図５】
2006339776.tif 000010
【図６】
2006339776.tif 000011
【図７】
2006339776.tif 000012
【図８】
2006339776.tif 000013
【図９】
2006339776.tif 000014
【図１０】
2006339776.tif 000015
【図１１】
2006339776.tif 000016
【図１２】
2006339776.tif 000017
【図１３】
2006339776.tif 000018
【図１４】
2006339776.tif 000019
【図１５】
2006339776.tif 000020
【図１６】
2006339776.tif 000021
【図１７】
2006339776.tif 000022