(57)【特許請求の範囲】
【請求項１】
ルックアップテーブルに基づいて、色彩信号とデバイス信号との間で色空間における座標点を変換する画像処理装置において、
色彩信号の色空間におけるグリッド点の座標のデータであるグリッドデータと所定の色のカラーパッチを測色して得られる測定データを含むデータとを外部から取り込む入力手段と、
前記グリッド点の座標位置をＳ字間数を用いて非線形なグリッド位置に変換するグリッド変換手段と、
前記グリッド点の座標に対応する前記非線形変換された座標を生成し、新たなルックアップテーブルを形成するルックアップテーブル生成手段と、
前記ルックアップテーブル生成手段によって生成されたルックアップテーブルを含む所定のデータを出力する出力手段と、
所定の色順応モデルを用いて、所定の光源下における擬似的な前記測定データを生成する色順応モデル処理手段とを有し、
前記ルックアップテーブル生成手段は、
前記色順応モデル処理手段で生成された擬似的な前記測定データ、及び、前記入力手段で取り込んだ前記測定データにそれぞれ所定の重み係数を掛け、該重み付けされた各測定データを和算して新たな測定データを算出し、その新たな測定データを用いて前記ルックアップテーブルを形成する色補正処理装置。
【請求項２】
前記画像処理装置において、前記グリッドデータの明度に関する信号のダイナミックレンジの補正と、コントラストの補正のためのγ補正を行うレンジγ補正手段をさらに有する請求項１記載の色補正処理装置。
【請求項３】
前記画像処理装置は、前記レンジγ補正手段での処理により生ずる彩度の変化を補正する彩度補正手段をさらに有する請求項２記載の色補正処理装置。
【請求項４】
前記ルックアップテーブル生成手段は、
前記非線形変換されたグリッド点の座標が、前記測定データから得られる色域内に含まれるか否かを判定する色域内グリッド判定手段と、
前記色域内に含まれると判定されたグリッド点については、前記色域内のグリッド点の座標と前記測定データを用いて所定の補間処理を行い、前記色域内に含まれるグリッド点の座標に対応するデバイス信号の座標を生成する４面体補間手段と、
前記色域内に含まれないと判定されたグリッド点については、所定の色域圧縮処理を行い、前記色域内に含まれないグリッド点の座標に対応するデバイス信号の座標を生成する色域圧縮手段とをさらに有する請求項１ないし３のいずれかに記載の色補正処理装置。
【請求項５】
ルックアップテーブルに基づいて、色彩信号とデバイス信号との間で色空間における座標点を変換する画像処理方法において、
色彩信号の色空間におけるグリッド点の座標のデータであるグリッドデータと所定の色のカラーパッチを測色して得られる測定データを含むデータとを外部から取り込む入力ステップと、
前記グリッド点の座標位置をＳ字間数を用いて非線形なグリッド位置に変換するグリッド変換ステップと、
前記グリッド点の座標に対応する前記非線形変換された座標を生成し、新たなルックアップテーブルを形成するルックアップテーブル生成ステップと、
前記ルックアップテーブル生成ステップで生成されたルックアップテーブルを含む所定のデータを出力する出力ステップと、
所定の色順応モデルを用いて、所定の光源下における擬似的な前記測定データを生成する色順応モデル処理ステップとを有し、
前記ルックアップテーブル生成ステップは、
前記色順応モデル処理ステップで生成された擬似的な前記測定データ及び前記入力ステップで取り込んだ前記測定データにそれぞれ所定の重み係数を掛け、該重み付けされた各測定データを和算して新たな測定データを算出し、その新たな測定データを用いて前記ルックアップテーブルを形成する色補正処理方法。
【請求項６】
前記画像処理方法において、前記グリッドデータの明度に関する信号のダイナミックレンジの補正と、コントラストの補正のためのγ補正を行うレンジγ補正ステップをさらに有する請求項５記載の色補正処理方法。
【請求項７】
前記画像処理方法は、前記レンジγ補正ステップでの処理により生ずる彩度の変化を補正する彩度補正ステップをさらに有する請求項６記載の色補正処理方法。
【請求項８】
前記ルックアップテーブル生成ステップは、
前記非線形変換されたグリッド点の座標が、前記測定データから得られる色域内に含まれるか否かを判定する色域内グリッド判定ステップと、
前記色域内に含まれると判定されたグリッド点については、前記色域内のグリッド点の座標と前記測定データを用いて所定の補間処理を行い、前記色域内に含まれるグリッド点の座標に対応するデバイス信号の座標を生成する４面体補間ステップと、
前記色域内に含まれないと判定されたグリッド点については、所定の色域圧縮処理を行い、前記色域内に含まれないグリッド点の座標に対応するデバイス信号の座標を生成する色域圧縮ステップとをさらに有する請求項５ないし７のいずれかに記載の色補正処理方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、カラー画像処理方法および装置に関する。特に、グレーバランスを向上させるために行う色補正処理方法および装置に関する。
【０００２】
【従来の技術】
近年、ＰＣやインターネット等の普及により、一般的なホームユースにおいてもカラー画像を扱う機会が増えてきている。例えば、デジタル・スチル・カメラで撮影した画像を手持ちのプリンタで出力するといったことが誰でも簡単にできるような時代になってきた。それに伴い、異種デバイス間でも画像を同じ色で再現しようという、いわゆるデバイス・インディペンデント・カラー（Device Independent Color : ＤＩＣ）の概念が発生し、その実現が求められてきた。これを実現するためのシステムのことを、一般に、カラー・マネージメント・システム（Color Management System:ＣＭＳ）と呼んでおり、その代表的なものとしてMacOSのColorsyncやWindows98のICM等がある。
【０００３】
ＣＭＳでは、入出力デバイスの色信号における物理的な測色値を一致させることによってＤＩＣを実現している（図１参照）。すなわち、ＣＭＳでは、入力装置で取り込まれた画像の色信号は、入力機器のデバイス・プロファイル（単に、プロファイルと呼ぶこともある）を用いてデバイスに依存しない色信号に変換され、さらに、デバイスに依存しない色信号は、出力機器のデバイス・プロファイルを用いて出力装置の画像の色信号に変換される（図２参照）。このように、入力系の色信号から出力系の色信号に変換するとき、デバイス・プロファイルと呼ばれる変換式または変換テーブルを用いてデバイスに依存しない色空間における色信号に、一度、変換することによってＤＩＣを実現している。ここで、デバイス・プロファイルとは、デバイス信号値（RGB、CMYK等）と色彩計等により測定したその信号値での色彩値（XYZ、L*a*b*等）との関係から算出したパラメータ群のファイルと考えて良い。
【０００４】
プリンタ等におけるように、デバイス信号値とそれによって実現される画像の色彩値との関係が非線形の関係にある場合は、デバイス・プロファイルとしてルックアップテーブル（ＬＵＴ）を使用するのが一般的である。ＬＵＴは、以下に示す２ステップにより作成される。
（１）入力色空間に均等格子を作成する。
（２）均等格子における出力値を算出する。
【０００５】
以下、CMY出力のプリンタで用いるＬＵＴの一般的な作成方法（従来手法）を説明する。まず、CMY空間で均等に配置されたＮ3（Ｎ×Ｎ×Ｎ）個のカラーパッチ等を分光測色器等により色彩値（CIE/XYZ、CIE/L*a*b*等）の測定を行う。カラーパッチはどのような配置にしても構わないが、デバイスの色空間を充分に満たすような配置とする必要がある。ＬＵＴは、この測定データの逆写像を算出することにより作成することができる（図３参照）。なお、以下の説明では、ＬＵＴのグリッド点をL*a*b*空間のグリッド点とし、それをL*a*b*グリッド点という。色彩値の測定データをL*a*b*色空間でのデータとしたとき（以下、この測定データをL*a*b*測定データという）、ＬＵＴは、L*a*b*空間（例えば、0≦L*≦100、-128≦a*≦128、-128≦b*≦128）を分割してＭ3（Ｍ×Ｍ×Ｍ）個のグリッドを形成し、各グリッドに対応するCMY出力値を算出し、それを表にして生成される。
【０００６】
図４は、測定データをもとに、各L*a*b*グリッド点に対応するCMY出力値を算出するまでの処理の流れを示す図である。Ｎ3個の測定データは、CMY空間上では均等に配置されているが、このデータ点を色彩空間L*a*b*上でプロットすると不規則な配置となってしまう（図５参照）。ＬＵＴ形成には、L*a*b*空間をＭ3（Ｍ×Ｍ×Ｍ）に分割し、そのグリッド点毎のCMY出力値を求めるわけであるが、図５からも分かるように全てのグリッド点が色域内に入っているわけではない。L*a*b*グリッド点が色域内かどうかの判定は、以下のようにして行う。
【０００７】
Ｎ3個の測定データは、(Ｎ-1)3 個の６面体で構成されている。CMY空間では全く歪みのない６面体、すなわち、立方体（図６（ａ）参照）は、L*a*b*空間では歪んだ６面体（図６（ｂ）参照）となっている。１つの６面体は５個の４面体に分割することが可能である（図７（ａ）、（ｂ）参照）。まず、L*a*b*グリッド点が、L*a*b*測定データの測定点がなす６面体を分割して形成される４面体（以下、測定データ４面体という）のうちのどの４面体に入っているかを探す。L*a*b*グリッド点が測定データ４面体に入っているかどうかは、次のように判別している。
【０００８】
図８に示すように、入力されたグリッド点を点Ｐ:（L*p, a*p, b*p）とし、ある測定データ４面体の頂点の座標（L*i, a*i, b*i）（i = 0, 1, 2, 3）をとする（図８（ａ）参照）。このとき、点Ｐは、以下の式（１）から算出されるα、β、χが以下の式（２）の条件を満たせば、その測定データ４面体に含まれていることになる。
【数１】
0004250863.tif 000002
【０００９】
いずれかの測定データ４面体において上式が成り立てば、その入力グリッド点に対応する点は色域内に含まれていることになる。入力されたグリッド点を含む測定データ４面体に対応するデバイス信号空間（ここでは、CMY）における４面体の各頂点のCMY値を（c*i, m*i, y*i）（i = 0, 1, 2, 3）とすると、入力されたグリッド点に対応するデバイス信号空間上のCMY値（c*p, m*p, y*p）は、以下の式（３）のように線形補間によって算出することが可能である（図８（ｂ）参照）。
【数２】
0004250863.tif 000003
【００１０】
分割されたどの４面体にも含まれなければ、その入力グリッド点に対応する点はデバイス色域外にあり、色域圧縮が施される。色域圧縮については、様々な手法が提案されており、例えば、次のような手法がある。
（１）色相を保存して明度-彩度の２次元平面上でする圧縮
（図９（ａ）、特開平０９−０９８２９８参照）
（２）人間の視覚に合った色差を定義し、その色差が最小となるような圧縮
（図９（ｂ）、特開平１０−０８４４８７参照）
このように測定データの逆写像を行い、色域内のグリッド点については４面体補間、色域外のグリッド点については色域圧縮を行うことによって、ＬＵＴを作成することができる。
【００１１】
【発明が解決しようとする課題】
しかしながら、このようなＬＵＴ作成方法では、a*、b*などの色度方向についても均等分割してグリッド位置を算出しているためにデバイス色域内に入るＬＵＴのグリッド点が少なく、その結果、色域内でデバイス信号と色彩信号を対応付ける情報が少ないため、色再現性、特に、グレーバランスが良好でないという問題がある。
【００１２】
また、デバイス相互で黒レベルが異なっているため、場合によっては、黒潰れや黒浮きなどが発生するという問題がある。
【００１３】
また、黒レベルの違いを補正する処理を行ったとしても、処理によっては彩度の変化も生じさせてしまい、同時に彩度方向の補正も必要になってしまうことがあるという問題がある。
【００１４】
また、ある光源下での測定データを逆写像することでＬＵＴを作成しているため、その光源下ではカラーマッチングするが、他の光源下では色が合わなくなってしまうという問題がある。
【００１５】
本発明は、上述の課題に鑑みてなされたもので、その目的とするところは、グレーバランスに優れ、または、複数の光源でも色再現が良好なカラー画像処理方法および装置を提供することにある。
【００１６】
【課題を解決するための手段】
上記の目的を達成するため、本発明は、ルックアップテーブルに基づいて、色彩信号とデバイス信号との間で色空間における座標点を変換する画像処理装置において、色彩信号の色空間におけるグリッド点の座標のデータであるグリッドデータと所定の色のカラーパッチを測色して得られる測定データを含むデータとを外部から取り込む入力手段と、前記グリッド点の座標位置をＳ字間数を用いて非線形なグリッド位置に変換するグリッド変換手段と、前記グリッド点の座標に対応する前記非線形変換された座標を生成し、新たなルックアップテーブルを形成するルックアップテーブル生成手段と、前記ルックアップテーブル生成手段によって生成されたルックアップテーブルを含む所定のデータを出力する出力手段と、所定の色順応モデルを用いて、所定の光源下における擬似的な前記測定データを生成する色順応モデル処理手段とを有し、前記ルックアップテーブル生成手段は、前記色順応モデル処理手段で生成された擬似的な前記測定データ、及び、前記入力手段で取り込んだ前記測定データにそれぞれ所定の重み係数を掛け、該重み付けされた各測定データを和算して新たな測定データを算出し、その新たな測定データを用いて前記ルックアップテーブルを形成するように構成される。以下、第１の発明という。
【００１７】
また、他の発明は、前記第１の発明において、前記画像処理装置が前記グリッドデータの明度に関する信号のダイナミックレンジの補正と、コントラストの補正のためのγ補正を行うレンジγ補正手段をさらに有するように構成される。以下、第２の発明という。
【００１８】
また、他の発明は、前記第２の発明において、前記画像処理装置は、前記レンジγ補正手段での処理により生ずる彩度の変化を補正する彩度補正手段をさらに有するように構成される。以下、第３の発明という。
【００２０】
また、他の発明は、前記第１ないし３のいずれかの発明において、前記ルックアップテーブル生成手段は、前記非線形変換されたグリッド点の座標が、前記測定データから得られる色域内に含まれるか否かを判定する色域内グリッド判定手段と、前記色域内に含まれると判定されたグリッド点については、前記色域内のグリッド点の座標と前記測定データを用いて所定の補間処理を行い、前記色域内に含まれるグリッド点の座標に対応するデバイス信号の座標を生成する４面体補間手段と、前記色域内に含まれないと判定されたグリッド点については、所定の色域圧縮処理を行い、前記色域内に含まれないグリッド点の座標に対応するデバイス信号の座標を生成する色域圧縮手段とをさらに有するように構成される。
【００２１】
【発明の実施の形態】
実施の形態１．
図１０は、本発明に係る一実施の形態の色補正処理装置１００の全体構成を示す図である。図１０において、色補正処理装置１００は、操作部１１０、制御部１２０、入出力部１３０、色空間変換部１４０、グリッド変換部１５０、色域内グリッド判定部１６０、４面体補間部１７０、色域圧縮部１８０、および記憶部１９０から構成される。制御部１２０は、他の全ての構成部と接続される。
【００２２】
以下に、上記各構成部の動作について説明する。以下の説明では、所定の色空間でのグリッドデータ（上記の均等格子のデータ等）と、所定のカラーパッチを測色して得られた測定データとが、外部から色補正処理装置１００に入力されることを前提とする。本実施の形態１においては、グリッドデータの色空間をCIE/L*a*b*とし、測定データの色空間をCIE/XYZとするが、これらに限られることは言うまでもない。なお、以下では、入力されたグリッドデータを入力グリッドデータと呼び、CIE/XYZ色空間の測定データをXYZ測定データと呼ぶ。
【００２３】
操作部１１０には、色補正処理を行うのに必要な設定等を行うためのメニュー画面が表示され、キー等を用いてその設定がされる。また、操作部１１０には、処理の開始・停止等の指示、データの入力、その他の入力がなされる。さらに、各構成部での処理の状態や処理結果等が表示される。操作部１１０の出力は、制御部１２０に、または、制御部１２０を介して各構成部に伝達される。
【００２４】
制御部１２０は、操作部１１０から出力される設定、指示、データ等に応じて各処理部に所定の指示やデータ等を出力し、各構成部から制御部１２０に出力される信号に応じて他の構成部に指示やデータ等を出力する等の処理を行う。例えば、後述するグリッド変換部１５０からの要求により、後述する記憶部１９０から所定のデータを読み出してグリッド変換部１５０に出力したり、グリッド変換部１５０が出力した所定のデータを４面体補間部１６０に出力したりする等の処理を行う。
【００２５】
入出力部１３０は、制御部１２０からの指示に従って、ＬＵＴのグリッドデータ（入力グリッドデータ）およびカラーパッチの測定データを含む所定のデータを取り込むと共に、本装置によって生成されたＬＵＴ及びその他の中間データを含む所定のデータを出力する。なお、この入出力部１３０に関する制御部１２０の制御は、操作部１１０に入力された指示に基づいてされる。また、データの入出力方法は、フレキシブル・ディスクやＣＤ等の所定の記録媒体を用いて行う方法でも、他の情報処理装置との通信によって行うものでも良い。
【００２６】
色空間変換部１４０は、制御部１２０の制御に従って、入出力部１３０から入力されたカラーパッチのXYZ測定データを公知の変換技術によりL*a*b*データに変換する。この変換には、一般に線形行列を用いるが、その他の方法で変換することを排除するものではない。以下では、この変換によって得られたL*a*b*データをL*a*b*測定データと呼ぶ。なお、L*a*b*測定データに対して公知のLagrange補間を施すのでも良く、この補間によって得られたデータを擬似L*a*b*測定データと呼ぶことにする。
【００２７】
グリッド変換部１５０の説明の前に、この構成部での処理の原理に関して、簡単な説明を行う。
一般に、ＬＵＴは、a*、b*等の色度方向においても均等分割してグリッド位置を算出しているため、グリッドに対応する点が色域内に入る割合が少なく、グレーバランスも良くない。ＬＵＴを作成する際に、グレー軸近傍の精度は、非常に重要である。これを実現するためには、ＬＵＴのグリッドの密度をグレー軸は細かくし、グレー軸から遠ざかる程荒くしていくことで精度の向上を図ることができる。例えば、ＬＵＴのグリッド位置を算出する際に、均等なグリッド位置に対してＳ字関数を用いた変換を施し、非線形なグリッド位置とすることなどは非常に効果的である。ここで、CIE 1994 Color Difference:ΔE94の彩度補正係数による関数とＳ字関数の代表的なシグモイド関数とを例にとって説明する。
【００２８】
上記のΔE94は、次式で表される。
【数３】
0004250863.tif 000004
このときの彩度方向の変化分をΔC*94とすると、彩度差は以下のように表すことができる。
【数４】
0004250863.tif 000005
ここで、C*76 >>ΔC*76 として両辺を積分すると、以下の式（６）が得られる。
【数５】
0004250863.tif 000006
以下では、この式を彩度補正関数と呼ぶこととする。
【００２９】
彩度補正関数を使用することによって、これまでのＬＵＴよりも視覚的にリニアなグリッド位置を持つことも可能となった。なお、上記の式（６）の逆関数を算出すれば逆Ｓ字関数として利用することができる。さらに、上記の式（６）において、定数０．００４５をパラメータkに置き換えて、以下の式（７）のように一般化することにより、様々なＳ字、逆Ｓ字関数を作成することができる。このパラメータkの値を大きくすることで、この関数で表される曲線の曲率が大きくなり、値を小さくすることで曲率が小さくなる（図１１参照）。
【数６】
0004250863.tif 000007
【００３０】
一方、シグモイド関数は、以下の式（８）のように表される。
【数７】
0004250863.tif 000008
シグモイド関数も式中のパラメータαを変更することによって、様々なＳ字関数を作成することができる（図１２参照）。このパラメータαを変更してΔE94における彩度補正関数にグレー軸で接するような値を求めたところ、α≒５であった。図１３は、α＝５としたときのシグモイド関数と彩度補正関数とをプロットした図である。彩度補正関数は、全体的に滑らかに変化しているのに対し、シグモイド関数は、グラフの端の部分で急激に変化している。図１３において、入力信号が８刻みで増加するとき、前者は１６．９増加し、後者は４４．９も増加していることことによるものである。
【００３１】
色補正処理装置１００の構成部の説明に戻って、グリッド変換部１５０は、制御部１２０の制御に従って、入出力部１３０からの入力グリッドデータをＳ字変換する。その変換は、例えば、上記の式（７）、（８）、すなわち、彩度補正関数やシグモイド関数に基づいて変換するのでも良い。その際用いるパラメータkやαは、操作部１１０から入力されるのでも、後述の記憶部１９０にあらかじめ記憶されているものを用いるのでも良い。以下では、この変換で得られたデータを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*グリッドデータに対して所定の補間を行い、色域内座標点に対応する点を出力装置の座標点（CMY、RGB等）を決定する。この補間には、上記の式（３）を用いるのでも良い。色域内L*a*b*グリッドデータと、この補間によって得られた出力装置の座標点データとの対応関係は、ＬＵＴの一部を構成する。これらのデータは、制御部１２０に出力される。
【００３４】
色域圧縮部１８０は、制御部１２０の制御に従って、色域内グリッド判定部１６０からの色域外L*a*b*グリッドデータに対して、上記、公知の色域圧縮処理を行う。既に述べたように色域圧縮処理には種々の公知の処理方法があり、どの処理方法を選ぶかは本発明の本質はなく、いずれの方法を適用するのでも良い。色域外L*a*b*グリッドデータと、この色域圧縮処理によって得られた出力装置の座標点データとの対応関係は、ＬＵＴの他の一部を構成する。これらのデータは、制御部１２０に出力される。
【００３５】
記憶部１９０は、制御部１２０からの制御に従って所定の情報を格納する。記憶される情報としては、制御部１２０に入力される入力情報のうち、例えば、操作部１１０におけるメニュー画面での設定内容、入出力部１３０から入力される上記のデータ、色空間変換部１４０からのL*a*b*測定データ、グリッド変換部１５０からのL*a*b*グリッドデータ、色域内グリッド判定部１６０から出力される所定のデータ、４面体補間部１７０及び色域圧縮部１８０から出力されるＬＵＴを構成するデータ等が含まれる。なお、これら以外のデータを、記憶されるデータから排除するものでないことは言うまでもない。
【００３６】
以下に、実施の形態１に係る色補正処理方法における処理の流れを説明する。図１４は、本発明の実施の形態１に係る色補正処理方法における処理の流れを示すフローチャートである。ステップＳ１４１０で、制御部１２０は、操作部１１０を介して入力された設定に従って、入力すべきデータの種類やＬＵＴ生成に用いられるパラメータkまたはαを設定する等の初期設定を行う。ステップＳ１４２０で、制御部１２０はステップＳ１４１０で設定された種類のデータを取り込む。
【００３７】
ステップＳ１４２０で、入出力部１３０は、ステップＳ１４１０で設定された、上記の入力グリッドデータおよび所定のカラーパッチの測定データ等を取り込む。本実施の形態１においては、入力グリッドデータをL*a*b*空間のデータ、カラーパッチの測定データをXYZ空間のデータとするが、それぞれのデータは、異なる空間のデータであることを排除するものではない。ステップＳ１４３０で、色空間変換部１４０は、カラーパッチのXYZ測定データからL*a*b*測定データに変換する。なお、このステップにおいて、L*a*b*測定データに対して、公知のLagrange補間等の所定の補間を行うのでも良い。
【００３８】
ステップＳ１４４０で、グリッド変換部１５０は、入力されたグリッドデータに所定の非線形の変換を施してL*a*b*グリッドデータを生成する。この変換は、例えば、上記に示すＳ字変換であっても良い。この場合は、例えば、上記の式（６）を用いた変換とすることでも良い。ステップＳ１４５０で、色域内グリッド判定部１６０は、上記のL*a*b*測定データ（または擬似L*a*b*測定データ）とL*a*b*グリッドデータを用いて、L*a*b*グリッドデータで示される各座標点がL*a*b*測定データで指定される色域内に入っているか否かを座標点毎に判断する。L*a*b*グリッドデータで示される各座標点が色域内と判断された場合は、処理はステップＳ１４６０に進み、色域外と判断された場合は、処理はステップＳ１４７０に移る。
【００３９】
ステップＳ１４６０で、４面体補間部１７０は、上記の色域内L*a*b*グリッドデータに対して所定の補間を行い、色域内座標点に対応する点を出力装置の座標点（CMY、RGB等）を決定する。この補間には、上記の式（３）を用いるのでも良い。色域内L*a*b*グリッドデータと、この補間によって得られた出力装置の座標点データとの対応関係は、ＬＵＴの一部を構成する。ステップＳ１４７０で、色域圧縮部１８０は、上記の色域外L*a*b*グリッドデータに対して、公知の色域圧縮処理を行う。どの色域圧縮処理方法を選ぶかは、本発明の本質はなく、いずれの方法を適用するのでも良い。色域外L*a*b*グリッドデータと、この色域圧縮処理によって得られた出力装置の座標点データとの対応関係は、ＬＵＴの他の一部を構成する。ステップＳ１４８０で、制御部１２０は、ステップＳ１４６０とステップＳ１４７０からの処理結果を統合してＬＵＴを形成する。そして、処理は終了する。
【００４０】
以上説明したように、本実施の形態１によれば、入力グリッドデータを、上記の彩度補正関数やシグモイド関数等のＳ字関数を用いて非線形変換することにしたので、従来のＬＵＴよりもグレー近傍におけるグリッドの密度を増やすことができ、グレーバランスの向上が可能な色補正処理方法および装置が実現できる。
また、色域内のグリッドの割合を増やすことができるためにＬＵＴ自体の精度の向上も可能である。
【００４１】
実施の形態２．
上記、実施の形態１における色補正処理方法および装置は、上記の彩度補正関数やシグモイド関数等のＳ字関数を使用して色域内のＬＵＴデータを増加させるという方法によって、従来のＬＵＴよりもグレー近傍におけるグリッドの密度を増やし、グレーバランスを向上させるものであった。しかし、このようにしても以下に示す黒潰れの問題は解決されずに残る。一般に、プリンタの最大濃度である黒（通常は、インクを全て打ったときの色）の明度L*、色度a*、b*を測定しても、（L*, a*, b*）＝（0, 0, 0）となることはない。通常のプリンタでは、明度L*は１０前後であり、黒の再現性の良いプリンタであってもせいぜい２〜３位である。また、a*、b*のずれはプリンタによって様々であり、これがグレーバランスの再現性の善し悪しに依存していることも少なくない。このことから、モニタなどの他のデバイスと繋いで測色的に忠実とされている比色分析的（Colorimetric）なＣＭＳを実行しても、最大濃度の黒以下の入力信号は、全て潰れてしまうことになる。これがいわゆる黒潰れを発生させてしまう原因の１つなのである。この黒潰れを改善するためには、明るさ方向におけるダイナミックレンジを最大限に再現できるような補正が有効である。本実施の形態２では、明るさ方向のダイナミックレンジを再現するための補正であるレンジ補正、およびγ補正を行って黒色補正する色補正処理方法および装置について記載する。
【００４２】
図１５は、実施の形態２に係る色補正処理装置１００の全体構成を示す図であり、その構成は、実施の形態１において示したものに、レンジγ補正部２００を追加したものである。以下の説明では、実施の形態１で用いた図１０に示す構成部と同様の構成部については同一の参照番号を付し、それらの説明を省略する。
【００４３】
実施の形態２において、操作部１１０は、実施の形態１に示された設定に加え、所定の表示をし、色信号のダイナミックレンジを補正するか否かの入力を待つ。補正するとされた場合には、ダイナミックレンジ補正に要する情報（以下、レンジγ補正情報と呼ぶ）の入力も待つ。または、レンジγ補正情報は、記憶部１９０に格納しておいたものを用いることも可能であり、他の方法を排除するものではない。入力があった場合にその情報を制御部１２０に出力し、制御部１２０では、この情報に基づいて各部の制御を行うと共に、レンジγ補正情報をレンジγ補正部２００に出力する。なお、レンジγ補正情報を記憶部１９０に格納するのでも良い。ダイナミックレンジを補正するとされた場合、レンジγ補正部２００は、以下に述べる処理を含むレンジγ補正を行う。
【００４４】
ここでは、以下に示す３種類の補正の手法を説明する（図１６参照）。
（１）プリンタの黒点ＢのY(1/3)チャンネルがL*=0となる補正式をＸ、Ｚチャンネルにも適用する方法
（２）プリンタの黒点ＢのX(1/3)、Y(1/3)、Z(1/3)チャンネルがすべて、L* = 0となるように補正する方法
（３）  プリンタの黒点ＢのL*が0となるように補正する方法
（１）、（２）、（３）のそれぞれの手法を適用することによって、黒点Ｂは、図中の矢印の方向B1、B2、B3にマップされることになる。手法（１）では、プリンタの白と黒を結んだ軸を延長させたときのL*=0 の点にマップしている。この手法が一般的に優れていると思われるが、プリンタの黒点の明度が高く色度が無彩軸から大きく離れている場合には、注意が必要である。手法（２）では、プリンタの黒点をL*=a*=b*=0 となるようにマップしている。一般的に、黒領域における色域の差を一番小さくすることが可能であるが、無彩軸を大きくずらしてしまうことが起きる可能性が高いことに注意が必要である。手法（３）については、明度L*だけの補正を行っているため色度点についての保持は可能である。
【００４５】
これらの黒色補正は、ダイナミックレンジの補正であるわけだが、レンジをリニアに補正してしまうとコントラストの点からオリジナルとずれてしまうことがある。このコントラストを補正するためには、式（９）および後述の式（１０）で表されるようなべき関数によるγ補正を行うのが効果的である。ここで、Drangeinは、入力側における明度のダイナミックレンジであり、Drangeoutは、出力側のレンジである。
【数８】
0004250863.tif 000009
ここで、Drangeは、それぞれの手法における色信号のダイナミックレンジを表しており、手法（１）では、Y(1/3)チャンネルにおけるダイナミックレンジを、手法（３）では、L*でのレンジを表している。また、手法（２）では、それぞれ、X(1/3)、Y(1/3)、Z(1/3)チャンネルにおけるダイナミックレンジを表しており、それぞれのチャンネルにおけるγ値も当然異なるようになる。なお、下付の文字「in」は入力デバイスを示し、下付の文字「out」は出力デバイスを示す。
【００４６】
図１７は、本実施の形態２に係る色補正処理方法における処理の流れを示すフローチャートであり、その処理の流れは、実施の形態１において示したもの（図１４参照）に、上記のレンジγ補正処理を追加したものである。以下の説明では、実施の形態１で用いた図１４に示す処理と同様の処理については同一の参照番号を付し、それらの説明を省略する。
【００４７】
図１７において、ステップＳ１４１０で、制御部１２０が初期設定するものに、上記の入力側のダイナミックレンジと出力側のダイナミックレンジとに関するデータ等のレンジγ補正処理に必要なデータの設定が含まれる。ステップＳ１４２０で、入出力部１３０が、ステップＳ１４１０で設定された、上記の入力グリッドデータおよび所定のカラーパッチの測定データ等を取り込んだ後、ステップＳ１７２５で、レンジγ補正部２００は、上記のレンジγ処理を行う。ステップＳ１４３０以降のステップで、実施の形態１同様の処理を行い、ＬＵＴを形成する。
【００４８】
図１８は、本実施の形態２に係るレンジγ補正処理（ステップＳ１７２５）の詳細な流れを示すフローチャートである。図１８において、ステップＳ１８１０で、制御部１２０は、明度のダイナミックレンジ補正を行うか否かを判断し、補正を行うと判断した場合は、処理はステップ１８２０に進み、行わないと判断した場合は、処理は終了する。ステップＳ１８２０で、制御部１２０は、操作部１１０からの指示等に基づき、どのような補正手法とするかを決定する。それらの補正手法には、上記の（１）から（３）の手法が含まれるのでも良く、ここでは、それらを含む場合について説明する。ここの判断結果に応じて、処理はステップＳ１８３０−１からステップＳ１８３０−ｎのいずれかのステップに進む。
【００４９】
制御部１２０は、ステップＳ１８３０−１で、上記（１）の補正手法を設定し、ステップＳ１８３０−２で、上記（２）の補正手法を設定し、ステップＳ１８３０−３で、上記（３）の補正手法を設定し、ステップＳ１８３０−４以降のステップで他の補正手法を設定する。補正手法が設定されたら、ステップＳ１８４０で、レンジ補正部２１０は、上記のステップＳ１８３０−１からステップＳ１８３０−ｎのステップで設定された補正手法でダイナミックレンジ補正を行う。ステップＳ１８５０で、レンジ補正部２１０は、上記の補正手法毎にコントラストを補正するためのγ補正を行う。以上の処理を持ってレンジγ補正処理のサブルーチンは終了する。
【００５０】
以上説明したように、本実施の形態２によれば、ダイナミックレンジを再現するようなレンジ補正処理を行うこととしたので、いわゆる黒潰れを防止できる。
また、レンジγ補正処理においてコントラストを補正するためのγ補正を行うこととしたので、ダイナミックレンジの補正によって生じるおそれがあるコントラストのずれを防止できる。
【００５１】
実施の形態３．
上記、実施の形態２における色補正処理方法および装置は、ダイナミックレンジを再現するような補正とコントラストを補正するためのγ補正を行って黒潰れを防止するものであった。しかし、上記の手法（１）、手法（２）では、上記の補正処理に伴って図１９に示すように彩度が変化してしまう。本実施の形態３では、彩度変化を抑制するために後述するべき関数を用いて彩度補正する色補正処理方法および装置について記載する。
【００５２】
図２０は、実施の形態３に係る色補正処理装置１００の全体構成を示す図であり、その構成は、実施の形態２において示したものに、彩度補正部２１０を追加したものである。以下の説明では、実施の形態２で用いた図１５に示す構成部と同様の構成部については同一の参照番号を付し、それらの説明を省略する。
【００５３】
本実施の形態３において、操作部１１０は、実施の形態２に示された設定に加え、所定の表示をし、彩度補正を行うか否かの入力を待つ。補正するとされた場合には、彩度補正に要する情報（以下、彩度補正情報と呼ぶ）の入力も待つ。または、彩度補正情報は、記憶部１９０に格納しておいたものを用いることも可能であり、他の方法を排除するものではない。入力があった場合にその情報を制御部１２０に出力し、制御部１２０では、この情報に基づいて各部の制御を行うと共に、彩度補正情報を彩度補正部２１０に出力する。なお、彩度補正情報を記憶部１９０に格納するのでも良い。彩度補正するとされた場合、彩度補正部２１０は、以下に述べる処理を含む彩度補正を行う。
【００５４】
上記の方法（１）および方法（２）によるレンジγ補正では、その補正に伴って図１９に示すように彩度が増加してしまう。そのため、これを抑制するような補正処理が必要となる。そこで、べき関数を使用した場合に、どれくらい変化するかを調べてみた。Macbeth ColorChecker２４色の無彩色を除く１８色に対して、式（１０）におけるγを０．２、０．４、０．６、０．８、１．０、１．２、１．４と変化させ、そのときの彩度方向の変化を式（１１）で近似し、この式のパラメータkを計算した。
【数９】
0004250863.tif 000010
【００５５】
結果は、表１に示すようになり、これを近似したところ式（１４）のように表せた（図２１参照）。
【表１】
0004250863.tif 000011
【数１０】
0004250863.tif 000012
【００５６】
以上のことから、彩度補正処理は、レンジγ補正で使用したべき関数のべき係数（γ）を使って式（１３）のように簡単に決定づけることも可能である。
【数１１】
0004250863.tif 000013
この彩度補正の効果は、図２２のようになり、充分に補正できていることが確認できる。
【００５７】
図２３は、本実施の形態３に係る色補正処理方法における処理の流れを示すフローチャートであり、その処理の流れは、実施の形態２において示したもの（図１７参照）に、上記の彩度補正処理（ステップＳ２３２９）を追加したものである。以下の説明では、実施の形態２で用いた図１７に示す処理と同様の処理については同一の参照番号を付し、それらの説明を省略する。
【００５８】
図２３において、ステップＳ１４１０で、制御部１２０が初期設定するものに、上記のべき関数のべき等の彩度正処理に必要なデータの設定が含まれる。ステップＳ１４２０で、入出力部１３０が、ステップＳ１４１０で設定された上記の入力グリッドデータおよび所定のカラーパッチの測定データ等を取り込み、ステップＳ１７２５で、レンジγ補正部２００が上記のレンジγ補正処理を行った後、ステップＳ２３２９で、彩度補正部２１０は、上記の彩度補正処理を行う。ステップＳ１４３０以降のステップで、実施の形態２同様の処理を行い、ＬＵＴを形成する。
【００５９】
以上説明したように、本実施の形態３によれば、デバイスのレンジ補正と明度方向のレンジγ補正とに加え、そこで用いたγ値と関連付けた彩度補正を行うこととしたため、黒レベルの異なるデバイスを想定しても黒潰れや黒浮きの無い、より再現性の良い画像出力を得ることができるようになった。
【００６０】
実施の形態４．
上記、実施の形態３における色補正処理方法および装置は、単一測定光源下で測色したカラーパッチの測定データを用いた。しかし、そのＬＵＴを用いたのでは、他の光源下ではカラーマッチしないことがあるという問題がある。本実施の形態４では、複数の光源下で適切にカラーマッチするＬＵＴを形成するための色補正処理方法および装置について記載する。
【００６１】
図２４は、本実施の形態４に係る色補正処理装置１００の全体構成を示す図であり、その構成は、実施の形態３において示したものに、色順応モデル処理部２２０を追加したものである。以下の説明では、実施の形態３で用いた図２０に示す構成部と同様の構成部については同一の参照番号を付し、それらの説明を省略する。
【００６２】
本実施の形態４において、操作部１１０は、実施の形態３に示された設定に加え、所定の表示をし、後述の色順応モデルを用いてＬＵＴを生成する処理（以下、色順応モデル処理という）を行うか否かの入力を待つ。色順応モデル処理を行うとされた場合には、モデルの種類等の色順応モデル処理に要する所定の情報（以下、色順応モデル情報と呼ぶ）の入力も待つ。または、色順応モデル情報は、記憶部１９０に格納しておいたものを用いることも可能であり、他の方法を排除するものではない。入力があった場合にその情報を制御部１２０に出力し、制御部１２０では、この情報に基づいて各部の制御を行うと共に、色順応モデル情報を色順応モデル処理部２２０に出力する。なお、色順応モデル情報を記憶部１９０に格納するのでも良い。色順応モデル処理を行うとされた場合、色順応モデル処理部２２０は、以下に述べる処理を含む色順応モデル処理を行う。
【００６３】
上記のように、ＬＵＴの１つの問題として、ＬＵＴを作成したときの測定光源下ではカラーマッチングはできるが、他の光源下では全く色が合わないといった演色性の問題がある。この問題を回避するためには、複数の光源下での測定データから、ある中間的なデータを作りだし、その中間的なデータからＬＵＴを作成すればよい。もしくは、複数の光源下におけるＬＵＴを作成し、それらのＬＵＴから中間的な色再現を提供できるＬＵＴを作成しても良い。例えば、太陽光（D65）と蛍光灯（F10）、双方の下での観察を考慮する場合には、D65光源での測定値XYZD65とF10光源での測定値XYZF10の組み合わせとして、新しく３刺激値を次のように定義する。
【数１２】
0004250863.tif 000014
このように算出した新しい３刺激値XYZからＬＵＴを作成すれば、D65光源下においてもF10光源下においても全く色は合わないといった問題はなくなる。
【００６４】
この例は、２種類の光源を想定したが、複数の光源を想定しても構わない。そのような場合も、測定値の係数の和が１になるように、それぞれの光源下における係数を決定すればよい。係数を決定する際には、想定する観察光源の中で最も多い可能性のある光源の係数を大きく設定すれば効果を大きくできる。また、他の光源の測定データが無い場合には、次に示すような色順応モデルを使用して、その光源下における３刺激値を算出しても良い。
【００６５】
測定データ（XYZデータ）を他の光源についてのデータに変換するのに様々な色順応モデルが提案されている。代表的なものに、Bradford変換、Hunt-Pointer-Estevas変換などがある。色順応モデルのフローは、図２５に示すように表すことができる。まず、特定の光源下の測定データ (XYZ)in は、３×３マトリックスによって錐体レベルへの刺激値 (LMS)adp に変換される。さらに、出力側の視環境から錐体レベルの刺激値 (LMS)out に変換され、そして、３×３マトリックスにより、 (XYZ)outに変換されることになる。ここで、XYZ３刺激値とLMS３刺激値を変換するための３×３マトリックスの係数が、Bradford変換とHunt-Pointer-Estevas変換で異なっており、それぞれ以下のようなマトリックス係数のなっている。
【数１３】
0004250863.tif 000015
【数１４】
0004250863.tif 000016
【００６６】
上記の３×３変換マトリックスをＭTとして表すと、色順応モデルによる変換後の (XYZ)outは、以下のように定義できる。
【数１５】
0004250863.tif 000017
ここで、Lmax_i、Mmax_i、Smax_iは、測定した光源のXYZをLMSに変換した値であり、Lmax_o、Mmax_o、Smax_oは、求めたい光源の光源のXYZをLMSに変換した値である。
【００６７】
図２６は、本実施の形態４に係る色補正処理方法における処理の流れを示すフローチャートであり、その処理の流れは、実施の形態３において示したもの（図２３参照）に、上記の色順応モデル処理を追加したものである。以下の説明では、実施の形態３で用いた図２３に示す処理と同様の処理については同一の参照番号を付し、それらの説明を省略する。
【００６８】
図２６において、ステップＳ１４１０で、制御部１２０が初期設定するものに、上記のα、β等の重み係数、色順応モデルの種類等の色順応モデル処理に必要なデータの設定が含まれる。ステップＳ１４２０で、入出力部１３０が、ステップＳ１４１０で設定された上記の入力グリッドデータおよび所定のカラーパッチの測定データ等を取り込んだ後、ステップＳ２６２３で、色順応モデル処理部２２０は、指定された色順応モデルを用いて上記の色順応モデル処理を行う。ステップＳ１７２５以降のステップで、実施の形態３同様の処理を行い、ＬＵＴを形成する。
【００６９】
以上説明したように、本実施の形態４によれば、色順応モデルを使用することによって、どのような光源についてもXYZ値を算出することが可能となり、観察光源が異なる場合では全く色が合わないという演色性の問題を回避することが可能となった。
【００７０】
【発明の効果】
以上説明したように、本発明によれば、入力グリッドデータを、上記の彩度補正関数やシグモイド関数等のＳ字関数を用いて非線形変換することにしたので、従来のＬＵＴよりもグレー近傍におけるグリッドの密度を増やすことができ、グレーバランスの向上が可能な色補正処理方法および装置が実現できる。
また、色域内のグリッドの割合を増やすことができるためにＬＵＴ自体の精度の向上も可能である。
【００７１】
また、本発明によれば、本ダイナミックレンジを再現するようなレンジ補正処理を行うこととしたので、いわゆる黒潰れを防止できる。
また、レンジγ補正処理においてコントラストを補正するためのγ補正を行うこととしたので、ダイナミックレンジの補正によって生じるおそれがあるコントラストのずれを防止できる。
【００７２】
また、本発明によれば、デバイスのレンジ補正と明度方向のレンジγ補正とに加え、そこで用いたγ値と関連付けた彩度補正を行うこととしたため、黒レベルの異なるデバイスを想定しても黒潰れや黒浮きの無い、より再現性の良い画像出力を得ることができるようになった。
【００７３】
また、本発明によれば、色順応モデルを使用することによって、どのような光源についてもXYZ値を算出することが可能となり、観察光源が異なる場合では全く色が合わないという演色性の問題を回避することが可能となった。
【図面の簡単な説明】
【図１】  ＤＩＣの概念を説明するための図である。
【図２】  ＤＩＣを実現するための処理の流れを模式的に示す図である。
【図３】  ＬＵＴ生成のためにする、測定データの逆写像を模式的に示す図である。
【図４】  測定データをもとに、L*a*b*空間における各グリッドに対応するCMY出力値を算出するまでの従来の処理の流れを示す図である。
【図５】  CMY空間上では均等に配置されているグリッドデータがL*a*b*空間上では不規則になることを模式的に示す図である。
【図６】  CMY空間では全く歪みのない６面体がL*a*b*空間では歪んだ６面体となることを模式的に示す図である。
【図７】  １つの６面体が５個の４面体に分割できることを説明するための図である。
【図８】  L*a*b*空間での点Ｐ:（L*p,a*p,b*p）が、ある４面体内にあるか否かを判断する処理を説明するための図である。
【図９】  色相を保存して明度-彩度の２次元平面上でする圧縮を説明するための図である。
【図１０】  本発明の実施の形態１に係る色補正処理装置の全体構成を示す図である。
【図１１】  本発明の実施の形態１に係るグリッド変換処理において用いるＳ字関数の一例を示す図である。
【図１２】  本発明の実施の形態１に係るグリッド変換処理において用いるＳ字関数の他の一例を示す図である。
【図１３】  α＝５としたときのシグモイド関数と彩度補正関数とを示す図である。
【図１４】  本発明の実施の形態１に係る色補正処理方法における処理の流れを示すフローチャートである。
【図１５】  実施の形態２に係る色補正処理装置１００の全体構成を示す図である。
【図１６】  実施の形態２に係るレンジγ補正処理の方法を模式的に示す図である。
【図１７】  実施の形態２に係る色補正処理方法における処理の流れを示すフローチャートである。
【図１８】  実施の形態２に係るレンジγ補正処理の詳細な流れを示すフローチャートである。
【図１９】  彩度補正を行う前のデータと測色データと彩度変化の一例を示す図である。
【図２０】  実施の形態３に係る色補正処理装置１００の全体構成を示す図である。
【図２１】  γとパラメータkとのフィッティング曲線を示す図である。
【図２２】  彩度補正後のデータと測色データとの一致性の一例を示す図である。
【図２３】  実施の形態３に係る色補正処理方法における処理の流れを示すフローチャートである。
【図２４】  実施の形態４に係る色補正処理装置１００の全体構成を示す図である。
【図２５】  実施の形態４に係る、色順応モデルを用いた色順応モデル処理の概念的な流れを示す図である。
【図２６】  実施の形態４に係る色補正処理方法における処理の流れを示すフローチャートである。
【符号の説明】
１００…色補正処理装置、１１０…操作部、１２０…制御部、１３０…入出力部、１４０…色空間変換部、１５０…グリッド変換部、１６０…色域内グリッド判定部、１７０…４面体補間部、１８０…色域圧縮部、１９０…記憶部、２００…レンジγ補正部、２１０…彩度補正部、２２０…色順応モデル処理部
【図１】
0004250863.tif 000018
【図２】
0004250863.tif 000019
【図３】
0004250863.tif 000020
【図４】
0004250863.tif 000021
【図５】
0004250863.tif 000022
【図６】
0004250863.tif 000023
【図７】
0004250863.tif 000024
【図８】
0004250863.tif 000025
【図９】
0004250863.tif 000026
【図１０】
0004250863.tif 000027
【図１１】
0004250863.tif 000028
【図１２】
0004250863.tif 000029
【図１３】
0004250863.tif 000030
【図１４】
0004250863.tif 000031
【図１５】
0004250863.tif 000032
【図１６】
0004250863.tif 000033
【図１７】
0004250863.tif 000034
【図１８】
0004250863.tif 000035
【図１９】
0004250863.tif 000036
【図２０】
0004250863.tif 000037
【図２１】
0004250863.tif 000038
【図２２】
0004250863.tif 000039
【図２３】
0004250863.tif 000040
【図２４】
0004250863.tif 000041
【図２５】
0004250863.tif 000042
【図２６】
0004250863.tif 000043