JP5665442B2 - Color processing apparatus and color processing method - Google Patents
Color processing apparatus and color processing method Download PDFInfo
- Publication number
- JP5665442B2 JP5665442B2 JP2010201340A JP2010201340A JP5665442B2 JP 5665442 B2 JP5665442 B2 JP 5665442B2 JP 2010201340 A JP2010201340 A JP 2010201340A JP 2010201340 A JP2010201340 A JP 2010201340A JP 5665442 B2 JP5665442 B2 JP 5665442B2
- Authority
- JP
- Japan
- Prior art keywords
- color
- value
- grid point
- values
- color material
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
- Color Image Communication Systems (AREA)
Description
本発明は、入力画像の色信号値を画像形成用の色材信号値へと変換する色分解処理の際に用いられる色分解プロファイルを作成する色処理装置および色処理方法に関する。 The present invention relates to a color processing apparatus and a color processing method for creating a color separation profile used in color separation processing for converting color signal values of an input image into color material signal values for image formation.
従来より、入力画像の色信号値を画像形成用の色材信号値へと変換する色分解処理を行う際には、色分解プロファイルとしてのルックアップテーブル(色分解LUT)が参照される。この色分解LUTは、代表的な色信号（代表色）に対応する色分解後の色材信号をユーザが設定しておき、この代表色についての色分解結果を補間することで、色空間全体についてのLUTとして作成されていた（例えば、特許文献１参照）。 Conventionally, when performing color separation processing for converting color signal values of an input image into color material signal values for image formation, a lookup table (color separation LUT) as a color separation profile is referred to. In this color separation LUT, the user sets a color material signal after color separation corresponding to a representative color signal (representative color), and interpolates the color separation result for this representative color, so that the entire color space (For example, refer to Patent Document 1).
ところが、上記色分解LUT作成方法においては、色材量を滑らかにするような補間が行われているものの、再現される色（例えばL*a*b*）に対する平滑性は考慮されていなかった。そのため、生成した色分解LUTを用いた画像形成を行った場合に、擬似輪郭が発生してしまうという課題があった。 However, in the above color separation LUT creation method, although interpolation is performed to smooth the amount of color material, smoothness for reproduced colors (for example, L * a * b *) has not been considered. . Therefore, there is a problem that a pseudo contour is generated when an image is formed using the generated color separation LUT.
この課題を解決するため、色材量から色を予測し、RGB空間に対し、各格子点の再現色が均等色空間上でより均等になるように、色分解LUTを平滑化する技術が提案されている（例えば、特許文献２参照）。 In order to solve this problem, a technology has been proposed that predicts the color from the amount of color material and smoothes the color separation LUT so that the reproduction color of each grid point in the RGB space is more uniform in the uniform color space. (For example, refer to Patent Document 2).
しかしながら、上記特許文献２に記載された色分解LUTの平滑化方法によれば、各格子点の配置は均等色空間上で均等になるものの、各格子点の間の点においては必ずしも均等な配置にならない、という課題がある。そのため、格子点の間に位置する点の色材量を四面体補間等により求めて出力する出力装置においては、色分解プロファイルの各格子点から四面体補間などにより予測した出力色と、実際に出力される色とに乖離が生じてしまう。このような乖離は、デバイス色とデバイス非依存の色（例えばL*a*b*）との対応関係を示したカラープロファイル（色分解LUTを含む）を用いて出力色を制御するカラーマネジメントシステムにおいて、色再現精度を低下させる要因となる。
However, according to the smoothing method of the color separation LUT described in
そこで本発明は、色分解プロファイルを用いたカラーマネジメントシステムにおいて出力色をより高精度に制御可能とするために、該色分解プロファイルを格子点間の色の補間精度が高くなるように作成することを目的とする。 Therefore, the present invention creates the color separation profile so that the color interpolation accuracy between the grid points becomes high in order to control the output color with higher accuracy in the color management system using the color separation profile. With the goal.
上記目的を達成するために、本発明の色処理装置は以下の構成を備える。 In order to achieve the above object, the color processing apparatus of the present invention comprises the following arrangement.
本発明にかかる色処理は、入力画像の色空間における格子点の色信号値と、画像形成装置の色材信号値との関係を示す色分解プロファイルを作成する際に、前記色空間における着目格子点の色信号値、および、該着目格子点に対して予め設定された複数方向に隣接する複数の隣接格子点の色信号値に対応する複数の色材信号値を取得し、前記取得した複数の色材信号値から、色予測モデルを用いて予測される均等色空間上の複数の色値を用いた補間により得られる、前記着目格子点と前記隣接格子点との間の複数の補間対象点それぞれに対する第一の色値を取得し、前記取得した複数の色材信号値を用いた補間から得られる色材信号値から、前記複数の補間対象点それぞれに対する第二の色値を取得し、前記複数の補間対象点に対する、前記第一の色値と前記第二の色値との色差から前記複数の補間対象点における補間精度を評価する評価値を算出し、前記評価値に基づき前記着目格子点の色材信号値を設定することを特徴とする。 The color processing according to the present invention is performed when a color separation profile indicating a relationship between a color signal value of a grid point in the color space of an input image and a color material signal value of the image forming apparatus is created. A plurality of color material signal values corresponding to the color signal values of the points and the color signal values of a plurality of adjacent grid points adjacent in a plurality of directions set in advance with respect to the target grid point; A plurality of interpolation targets between the target grid point and the adjacent grid point obtained by interpolation using a plurality of color values in a uniform color space predicted using a color prediction model from the color material signal value of Obtaining a first color value for each of the points, and obtaining a second color value for each of the plurality of interpolation target points from a color material signal value obtained by interpolation using the plurality of obtained color material signal values. , For the plurality of interpolation target points, Calculating an evaluation value for evaluating the interpolation accuracy of the plurality of interpolation target point from the color difference between one color value and the second color values, to set the color material signal values of the grid point of interest on the basis of the evaluation value It is characterized by that.
本発明によれば、格子点間の色の補間精度が高い色分解プロファイルを作成することができるため、該色分解プロファイルを用いたカラーマネジメントシステムにおいて出力色をより高精度に制御可能となる。 According to the present invention, it is possible to create a color separation profile with high accuracy for interpolating colors between grid points, so that it is possible to control output colors with higher accuracy in a color management system using the color separation profile.
以下、本発明の実施形態について、図面を参照して説明する。なお、以下の実施の形態は特許請求の範囲に関る本発明を限定するものではなく、また、本実施の形態で説明されている特徴の組み合わせの全てが本発明の解決手段に必須のものとは限らない。 Hereinafter, embodiments of the present invention will be described with reference to the drawings. The following embodiments do not limit the present invention related to the scope of claims, and all combinations of features described in the present embodiments are essential to the solution means of the present invention. Not necessarily.
＜第１実施形態＞
●システム構成
本実施形態においては、カラーマネジメントシステムにおいて用いられる、色分解プロファイルを、その格子点間の色の補間精度が高くなるように作成する方法を示す。ここで色分解プロファイルとは、入力画像の色空間における格子点の色信号値と、画像形成装置（プリンタ）に対する色材信号値との関係を示すものである。以下、本実施形態では、色処理システムにおいて該プロファイル作成方法を実現するプログラムをアプリケーションとして機能させた場合を例として説明する。
<First Embodiment>
System Configuration In this embodiment, a method for creating a color separation profile used in a color management system so as to increase the accuracy of color interpolation between grid points will be described. Here, the color separation profile indicates the relationship between the color signal value of the grid point in the color space of the input image and the color material signal value for the image forming apparatus (printer). Hereinafter, in the present embodiment, a case where a program for realizing the profile creation method in the color processing system functions as an application will be described as an example.
図１は、本実施形態における色処理システムのブロック構成を示す図である。同図において、100はパーソナルコンピュータ等の電子計算機（以下、PC）であり、ハードディスク109、CPU110、RAM111、ROM112等のハードウェアによって構成される。また、105はプリンタであり、画像信号をプリンタ105で印刷出力する際には、PC100においてアプリケーション101からOS102に印刷出力要求を行う。OS102はアプリケーション101の出力要求を受け、プリンタ105に対応するプリンタドライバ103に、描画命令群を発行する。プリンタドライバ103は、OS102から入力した印刷要求と描画命令群を処理し、プリンタ105で印刷可能な印刷データを作成してプリンタ105に転送する。プリンタ105がラスタープリンタである場合には、プリンタドライバ103は、OS102からの描画命令に対して順次画像補正処理を行い、例えばRGB24ビットページメモリにラスタライズする。そして、全ての描画命令をラスタライズした後に、RGB24ビットページメモリの内容をプリンタ105が印刷可能なデータ形式（例えば色材色に対応したC,M,Y,K,Lc,Lmデータ）に変換し、これをプリンタ105に転送する。また、108は測色器であり、107は測色器108の制御を行うための測色器ドライバである。測色器108および測色器ドライバ107は、後述する色分解プロファイル作成処理において、プリンタ105で印字されたパッチ画像のL*a*b*値を測定し、その結果をPC100に入力する。106はモニタであり、OS102からの描画命令にしたがって画像信号やユーザインタフェース等の表示を行う。
FIG. 1 is a diagram showing a block configuration of a color processing system in the present embodiment. In the figure,
上述したように、本実施形態の色処理システムにおいては、アプリケーション101のひとつとして色分解プロファイル作成用のアプリケーションが動作することにより、プリンタドライバ103において参照される色分解プロファイルが作成される。図２は、色分解プロファイル作成アプリケーション205とプリンタドライバ103における概略機能構成を示す図である。201はカラーマッチング処理部であり、画像信号についてモニタ106上で再現される色とプリンタ105で印刷した色とが合うように、該画像信号のRGB値を変換する。この変換は、モニタ106に入力されるRGB値と表示される色との関係、およびプリンタ105に入力されるRGB値と印字される色との関係、をそれぞれ示したプロファイルを参照することで行われる。具体的にはプロファイルとしてICC（International Color Consortium）プロファイルを用いれば良い。さらに、周知のCMM（Color Management Module）が適用可能であるが、本発明は特にその変換方法を限定するものではないため、詳細な説明については割愛する。
As described above, in the color processing system of the present embodiment, a color separation profile that is referred to by the
202は色分解処理部であり、カラーマッチング処理部201から出力されたR'G'B'の多値データを入力値とし、プリンタ105の各色材の色情報を出力値とした変換（色分解処理）を行う。この出力値としては例えばC（シアン）、M（マゼンタ）、Y（イエロー）、K（ブラック）、Lc（淡シアン）、Lm（淡マゼンタ）等、プリンタ105が有する色材種に対応した次元となる。色分解処理部202における色分解処理は、プリンタ105で使用される印刷用紙やモードに応じた色分解プロファイルを色分解プロファイル保持部204から読み込み、参照することで行われる。
ここで、本実施形態で用いられる色分解プロファイルについて、図３を用いて詳細に説明する。色分解プロファイルは、図３(a)に示すような、RGBで表現される3次元色空間上の立方体を格子状に分割した各格子点に対応する、C,M,Y,K,Lc,Lmのそれぞれ多値の色情報（以下、色材信号値と称する）を格納したルックアップテーブルである。より具体的には、図３(b)に示すように、RGBの3次元色空間での各格子点の位置（座標）情報と、その位置における各色材信号値が対応づけられた形式の情報である。図３(b)において、R,G,B列は、各格子点の位置情報を示す。また、C,M,Y,K,Lc,Lm列は、同一行の格子点位置情報(R,G,B)に対応した、各色材の色材信号値を示す。本実施形態では、入力のRGBデータとして0〜255の値をとる8ビットデータを想定している。また、図３(b)に示す色分解プロファイルは、RGBのそれぞれについて、0〜255の値を32間隔で分割して得られる9点でサンプリングした、全9×9×9＝729点の格子点に対応する色材信号値を格納している。なお、格子点のサンプリング間隔やサンプリング点数、入力ビット数は上記の例に限るものでなく、任意の値とすることが可能である。 Here, the color separation profile used in the present embodiment will be described in detail with reference to FIG. The color separation profile is represented by C, M, Y, K, Lc, corresponding to each lattice point obtained by dividing a cube in a three-dimensional color space expressed in RGB as shown in FIG. This is a look-up table storing multi-value color information (hereinafter referred to as color material signal values) of Lm. More specifically, as shown in FIG. 3B, the position (coordinate) information of each grid point in the RGB three-dimensional color space and information in a format in which each color material signal value at that position is associated. It is. In FIG. 3B, columns R, G, and B indicate position information of each lattice point. The C, M, Y, K, Lc, and Lm columns indicate the color material signal values of the respective color materials corresponding to the lattice point position information (R, G, B) in the same row. In the present embodiment, 8-bit data taking values of 0 to 255 is assumed as input RGB data. In addition, the color separation profile shown in FIG. 3B is a grid of 9 × 9 × 9 = 729 points obtained by sampling 9 points obtained by dividing the values of 0 to 255 at 32 intervals for each of RGB. The color material signal value corresponding to the point is stored. Note that the sampling interval of the grid points, the number of sampling points, and the number of input bits are not limited to the above example, and can be set to arbitrary values.
色分解処理部202では、入力されたR'G'B'データが、色分解プロファイルの格子点上にない場合、その近傍の格子点データを用いた補間を行うことで、該入力データに対応するC,M,Y,K,Lc,Lmの色材信号値を得る。この補間方法としては一般的に四面体補間が多く用いられるが、立方体補間等、その他の既知のアルゴリズムを適用することも可能である。
In the color
203はハーフトーン処理部であり、色分解処理部202から出力されたC,M,Y,K,Lc,Lm多値データを入力し、プリンタ105で再現可能なC,M,Y,K,Lc,Lmの各2ビットデータに変換する。具体的には、例えばベイヤー型の16×16のマトリクスをC,M,Y,K,Lc,Lmの画像それぞれにあてがい、マトリクス要素と対応する画像上の画素値を比較することで2値化を行う。
すなわち、画素値がマトリクス要素よりも大きい場合には当該画素値を1とし、逆にマトリクス要素以下である場合には当該画素値を0とする。なお、このハーフトーニング手法としては誤差拡散法等の他の手法を用いても良い。このようにして得られた、プリンタで印刷可能なC,M,Y,K,Lc,Lmの各2ビットデータは、プリンタ105に送られて印刷出力される。
That is, when the pixel value is larger than the matrix element, the pixel value is set to 1. Conversely, when the pixel value is equal to or less than the matrix element, the pixel value is set to 0. Note that other methods such as an error diffusion method may be used as the halftoning method. The 2-bit data of C, M, Y, K, Lc, and Lm that can be printed by the printer is sent to the
色分解プロファイル作成アプリケーション205は、本実施形態の色分解プロファイル作成処理を実行するアプリケーションであり、色分解プロファイル作成部206とプリンタ特性入力部207から構成される。プリンタ特性入力部207は、色分解プロファイル作成時に必要となるパッチ画像の測色データ等を、色分解プロファイル作成部206に入力する。色分解プロファイル作成部206では、入力されたパッチ画像の測色データ等に基づいて色分解プロファイルを作成し、該作成された色分解プロファイルは、プリンタドライバ103内の色分解プロファイル保持部204に保存される。
The color separation
●色分解プロファイル作成処理
以下、色分解プロファイル作成部206における色分解プロファイル作成処理について、図４のフローチャートを用いて詳細に説明する。まずS401において、初期の色分解プロファイル（以下、初期プロファイル）を作成する。この初期プロファイルはすなわち、ユーザが代表的な格子点の色情報（色材信号値）を調整することによって作成されるが、その作成方法の詳細については後述する。S401で設定された初期プロファイルは、各格子点の配置が、均等色空間（例えばCIＥL*a*b*空間）上でコントロールされていない。つまりは、RGB空間で均等に配置された各格子点が、均等色空間上で非均等に配置されてしまうことになる。このような色分解プロファイルを使用して印字を行った場合、擬似輪郭の発生等の弊害が生じる可能性があり、ひいては、上述したカラーマッチング処理部201の色変換精度が低下する要因となる。そこで本実施形態では後述するS403で格子点最適化処理を行うことによって、各格子点が均等色空間上でより好適に位置するよう、色情報を決定することを特徴とする。
Color Separation Profile Creation Processing Hereinafter, the color separation profile creation processing in the color separation
図４に戻り、次にS402においてプリンタ特性入力部207より、後段のS403における予測処理に必要となるカラーパッチ画像の測色値入力処理を行う。S403では格子点ごとに、取得した色材信号値に応じて均等色空間での色を予測し、均等色空間上で各格子点が好適な位置となるように、各格子点に対応する色材信号値を更新する格子点最適化処理を行う。このときの色予測の方法としては、周知のセル分割ユールニールセン分光ノイゲバウアモデル（Cellular Yule-Nielsen Spectral Neugebauer Model）を適用することができる。このモデルでは、まずC,M,Y,K,Lc,Lmの各多値データを任意の間隔でサンプリングして、C,M,Y,K,Lc,Lm空間上で格子点を形成する（セル分割）。次に、全格子点についての色材信号値の組合せからなる色票を作成してプリンタ105で印字し、該印字された色票を測色器108で測色することにより、全格子点の分光反射率を取得する。格子点外に位置する色材信号値の組合せについての分光反射率は、近傍の格子点の分光反射率からの線形補間によって求めることができる（分光ノイゲバウアモデル）。このとき、ドットゲインによる非線形性をユールニールセンモデルを用いて考慮することで、予測精度を向上させることが可能である。S403にてセル分割ユールニールセン分光ノイゲバウアモデルを用いた色予測を行う場合、S402においては上述した全格子点の色材信号値の組合せからなる色票の作成、印字、測色を行い、該測色結果をアプリケーションに入力する。なお、S403での色予測方法としては、上述したセル分割ユールニールセン分光ノイゲバウアモデルに限らず、例えば分光反射率に替えて三刺激値XYZを測色し、ノイゲバウアモデルを適用して任意の色材信号値の組合せに対するXYZ値を予測してもよい。なお、S403における格子点最適化処理の詳細については後述する。
Returning to FIG. 4, next, in step S402, the printer
最後にS404において、S403で格子点位置が最適化された色分解プロファイルを色分解プロファイル保持部204にアップロードし、処理を終了する。
Finally, in S404, the color separation profile whose grid point position has been optimized in S403 is uploaded to the color separation
●初期プロファイル作成処理(S401）
以下、上記S401における初期プロファイルの作成処理について、詳細に説明する。ここでは、初期プロファイルをユーザ指示に基づいて作成する例について説明する。
● Initial profile creation process (S401)
Hereinafter, the initial profile creation process in S401 will be described in detail. Here, an example of creating an initial profile based on a user instruction will be described.
まず、入力画像の色空間上での、色分解プロファイルの格子点における各頂点であるW,C,M,Y,R,G,B,Bkの色材信号値を決定する。図５に、各頂点W,C,M,Y,R,G,B,BkのRGB空間上での位置を示す。ここで、各頂点W,C,M,Y,R,G,B,BkとRGB空間上での格子点の座標(R,G,B)は、以下の対応関係があるものとする。 First, the color material signal values of W, C, M, Y, R, G, B, and Bk, which are the vertices of the grid points of the color separation profile, on the color space of the input image are determined. FIG. 5 shows the positions of the vertices W, C, M, Y, R, G, B, and Bk on the RGB space. Here, the vertexes W, C, M, Y, R, G, B, Bk and the coordinates (R, G, B) of the grid points in the RGB space have the following correspondence.
W＝(255,255,255)、ホワイト、即ちプリント用紙の色を示す
C＝(0,255,255)、シアン原色を示す
M＝(255,0,255)、マゼンタ原色を示す
Y＝(255,255,0)、イエロー原色を示す
R＝(255,0,0)、レッド原色を示す
G＝(0,255,0)、グリーン原色を示す
B＝(0,0,255)、ブルー原色を示す
Bk＝(0,0,0)、ブラック、即ちプリンタの最暗点を示す。
W = (255,255,255), indicating the color of white, that is, print paper
C = (0,255,255), indicating the cyan primary color
M = (255,0,255), indicating the magenta primary color
Y = (255,255,0), indicating the yellow primary color
R = (255,0,0), indicating the red primary color
G = (0,255,0), indicating the green primary color
B = (0,0,255), indicating the blue primary color
Bk = (0,0,0), black, that is, the darkest point of the printer.
W,C,M,Y,R,G,B,Bkの各頂点に対する色材信号値は、例えば次のように決定することができる。まず、Wは紙白であるので、色材信号値は全て0に設定する。次に、C,M,Y,Bkについては、混色による色の濁りを防ぐため、対応する1種類の色材のみの色材信号値を、印字する用紙やモードに応じて印字可能な最大の色材量とし、その他の色材の色材信号値は0とする。そして2次色であるR、G、Bについては、例えばRの場合、まず色材Mおよび色材Yの2種の色材を用い、合計が最大色材量となる複数の組合せのカラーパッチデータを生成し、プリンタ105にて印刷する。そしてその印刷結果をユーザが目視で確認し、Rに対応する色として好ましい色となる色材Mおよび色材Yの色材量の組合せを選択する。または、印刷したカラーパッチを測色器108を用いて測色し、M、Yの各1次色における色相の中間の色相となる色材の組合せを決定するようにしても良い。なお、G、BについてもRの場合と同様に、Gは色材Yと色材Cの組み合わせ、Bは色材Cと色材Mの組合せから、それぞれ決定できる。なお、各頂点の色材信号値の決定方法は上記方法に限るものではない。
The color material signal value for each vertex of W, C, M, Y, R, G, B, and Bk can be determined as follows, for example. First, since W is paper white, all color material signal values are set to zero. Next, for C, M, Y, and Bk, in order to prevent color turbidity due to mixed colors, the color material signal value of only one type of corresponding color material is the maximum that can be printed according to the printing paper and mode. The color material amount is assumed, and the color material signal values of other color materials are set to 0. For R, G, and B, which are secondary colors, for example, in the case of R, first, two color materials of color material M and color material Y are used, and a plurality of combinations of color patches whose total color material amount is the maximum. Data is generated and printed by the
次に、図５に示した各頂点を結ぶ19本の直線上に位置する格子点における色材信号値を決定する。具体的には、頂点W−C,W−M,W−Y,W−R,W−G,W−B,C−Bk,M−Bk,Y−Bk,R−Bk,G−Bk,B−Bk,R−Y,Y−G,G−C,C−B,B−M,M−R,W−Bkのそれぞれを結ぶ全19本の直線である。図６(a)は、頂点Wを含む直線上の格子点における色分解例を示す。同図において横軸は図５に示すW−Cの直線上にある各格子点を示し、縦軸は各格子点に対応した色材信号値および総色材信号値を示している。図中、Wでは全色材の色材信号値が0である。またCにおいては、色材Cが最大色材量となり、その他の色材は0となる。本実施形態では低濃度領域の粒状性を向上させるために淡シアン(Lc)を用いるため、まずWの次の格子点から色材Lcの使用を開始し、格子点ごとにLc量が増していくが、途中からより発色の良い色材Cに切り替えるよう色分解が決定されている。この決定は、例えば図６(a)に示す図をUIとしてモニタ106に表示し、各格子点に対する各色材のインク値をユーザが手動で調整することによって行う。また、該調整結果に基づいてカラーパッチデータを生成し、プリンタ105にて印刷し、該印刷結果における粒状性や色味をユーザが目視で判断できるよう構成してもよい。同様に、頂点Bkを含む直線上の格子点における色分解例を図６(b)に、頂点WおよびBkを含まない直線上の格子点における色分解例を図６(c)に、それぞれ示す。なお、図６(a)〜(c)に示した色分解例およびその決定方法はあくまでも例であって、本発明はこの例に限定されるものではない。
Next, the color material signal values at the lattice points located on the 19 straight lines connecting the vertices shown in FIG. 5 are determined. Specifically, vertices W-C, W-M, W-Y, W-R, W-G, W-B, C-Bk, M-Bk, Y-Bk, R-Bk, G-Bk, There are a total of 19 straight lines connecting B-Bk, RY, Y-G, GC, CB, B-M, MR, and W-Bk. FIG. 6A shows an example of color separation at a lattice point on a straight line including the vertex W. FIG. In the figure, the horizontal axis indicates each grid point on the WC line shown in FIG. 5, and the vertical axis indicates the color material signal value and the total color material signal value corresponding to each grid point. In the figure, for W, the color material signal values of all color materials are zero. In C, the color material C is the maximum color material amount, and the other color materials are 0. In this embodiment, since light cyan (Lc) is used to improve the granularity of the low density region, the use of the color material Lc is first started from the grid point next to W, and the amount of Lc increases for each grid point. However, the color separation is determined so as to switch to the color material C with better color development from the middle. This determination is performed by, for example, displaying the diagram shown in FIG. 6A on the
最後に、上記19本の直線上に位置する格子点の色材信号値に基づき、上記以外の全ての格子点における色材信号値を決定する。この方法としては、上記特許文献１に記載された方法等を適用することができる。該方法によれば、色空間をWとBkを含む6つの四面体（W−R−Y−Bk,W−Y−G−Bk,W−G−C−Bk,W−G−B−Bk,W−B−M−Bk、W−M−R−Bk）に分割し、各四面体で補間処理を行うものである。各四面体内部の色材信号値は、四面体を形成する4つの三角形において、上記処理で既知となる各辺の色分解結果から、色材ごとの色材信号値の等高線を引くことで決定される。
Finally, based on the color material signal values of the grid points located on the 19 straight lines, the color material signal values at all other grid points are determined. As this method, the method described in
以上、初期プロファイルをユーザの調整により作成する例を示したが、この作成方法は上記方法に限らず、例えば類似する用紙やモードに対応した既存の色分解プロファイルを、初期プロファイルとして適用しても良い。 The example in which the initial profile is created by user adjustment has been described above. However, this creation method is not limited to the above method, and for example, an existing color separation profile corresponding to a similar sheet or mode may be applied as the initial profile. good.
●格子点最適化処理(S403)
以下、上記S403における格子点最適化処理について、図７のフローチャートを用いて詳細に説明する。まずS701で変数r'g'b'を初期化する。具体的には、対象とする格子点（着目格子点）の位置情報であるRGBと同じ値、すなわち、初期プロファイルである図３(b)のRGB列の値をそのまま、変数r'g'b'の初期値とする。なお、本実施形態で作成される色分解プロファイルの各格子点における色信号値は、初期プロファイルの設定された値、すなわち変数r'g'b'の初期値と同じとなる。次にS702にて、変数r'g'b'に対する補間精度評価値Ｅhを算出するが、その詳細については後述する。そしてS703で、S702における処理によって算出された補間精度評価値Ｅhが十分小さいか否かに応じて処理を分岐する。すなわち、補間精度評価値Ｅhが所定の閾値（第１の閾値）以上である場合にはS704に進み、変数r'g'b'を所定のアルゴリズム（準ニュートン法や減衰最小二乗法など）にしたがって更新して設定した後、処理をS702に戻す。この変数r'g'b'の設定はすなわち、着目格子点のRGB値を中心として、該着目格子点に隣接する隣接格子点のRGB値を超えない範囲で行う。
Lattice point optimization processing (S403)
Hereinafter, the lattice point optimization process in S403 will be described in detail with reference to the flowchart of FIG. First, in step S701, a variable r'g'b 'is initialized. More specifically, the same value as RGB, which is the position information of the target grid point (target grid point), that is, the value of the RGB column in FIG. The initial value of '. Note that the color signal value at each grid point of the color separation profile created in the present embodiment is the same as the value set in the initial profile, that is, the initial value of the variable r′g′b ′. Next, in S702, an interpolation accuracy evaluation value Eh for the variable r′g′b ′ is calculated, and details thereof will be described later. In S703, the process branches depending on whether the interpolation accuracy evaluation value Eh calculated by the process in S702 is sufficiently small. That is, if the interpolation accuracy evaluation value Eh is equal to or greater than a predetermined threshold (first threshold), the process proceeds to S704, and the variable r'g'b 'is set to a predetermined algorithm (such as a quasi-Newton method or an attenuation least square method). Therefore, after updating and setting, the process returns to S702. That is, the variable r′g′b ′ is set in a range that does not exceed the RGB value of the adjacent grid point adjacent to the target grid point, with the RGB value of the target grid point as the center.
一方、補間精度評価値Ｅhが第１の閾値未満である場合にはS705に進み、変数r'g'b'に対する各色材の色材信号値を初期プロファイルを参照することで求める。そして、得られた色材信号値を、作成される色分解プロファイルにおける当該格子点（着目格子点）の更新前の色信号値に対する色材信号値として設定する。ここで、色分解プロファイルにおける当該格子点の更新前の色信号値とは、S701で設定された変数r'g'b'初期値である。なお、r'g'b'値に対応する各色材信号値は、後述するようにS702の補間精度評価値算出時にC1,M1,Y1,K1,Lc1,Lm1として取得されているため、S705ではそれらの値をそのまま用いれば良い。 On the other hand, if the interpolation accuracy evaluation value Eh is less than the first threshold value, the process proceeds to S705, and the color material signal value of each color material for the variable r'g'b 'is obtained by referring to the initial profile. Then, the obtained color material signal value is set as a color material signal value for the color signal value before the update of the grid point (target grid point) in the color separation profile to be created. Here, the color signal value before the update of the grid point in the color separation profile is the initial value of the variable r′g′b ′ set in S701. Since each color material signal value corresponding to the r'g'b 'value is acquired as C1, M1, Y1, K1, Lc1, Lm1 when calculating the interpolation accuracy evaluation value in S702 as described later, in S705 Those values may be used as they are.
なおS703において、S704およびS702の処理を所定回数繰り返しても補間精度評価値Ｅhが所定の閾値未満にならない場合には、S705に分岐してループ処理を抜ける。このときS705では、当該着目格子点についてS702で算出した補間精度評価値Ｅhのうち、最小であったものに対応するr'g'b'値を決定する。そして、該r'g'b'値について初期プロファイルから得られる各色材の色材信号値を、色分解プロファイルにおける当該格子点の色材信号値に設定する。 In S703, if the interpolation accuracy evaluation value Eh does not become less than the predetermined threshold even if the processes of S704 and S702 are repeated a predetermined number of times, the process branches to S705 and the loop process is exited. At this time, in S705, the r′g′b ′ value corresponding to the smallest one of the interpolation accuracy evaluation values Eh calculated in S702 for the target lattice point is determined. Then, the color material signal value of each color material obtained from the initial profile for the r′g′b ′ value is set to the color material signal value of the grid point in the color separation profile.
上記S701〜S705による最適化処理を、作成する色空間プロファイルの全格子点について行うため、S706で全ての格子点について処理を行ったか否かを判定する。全格子点についての処理が終了していれば格子点最適化処理を終了するが、未処理の格子点がある場合にはS701に戻り、該格子点についての処理を行う。 In order to perform the optimization process in S701 to S705 for all grid points of the color space profile to be created, it is determined whether or not the process has been performed for all grid points in S706. If the processing for all grid points has been completed, the grid point optimization processing is terminated. If there is an unprocessed grid point, the process returns to S701 to perform processing for the grid point.
なお、S706で全格子点についてS701〜S705による最適化処理が終了したと判定された後にも、作成した色分解プロファイルを新たな初期プロファイルとして、その全格子点についてS701〜S705の最適化処理を再度行うことも有効である。このように格子点最適化処理を所定回数繰り返すことによって、より補間精度の高い色分解プロファイルを作成することができる。 Even after it is determined in S706 that the optimization process of S701 to S705 has been completed for all grid points, the created color separation profile is used as a new initial profile, and the optimization process of S701 to S705 is performed for all grid points. It is also effective to do it again. In this way, by repeating the grid point optimization process a predetermined number of times, a color separation profile with higher interpolation accuracy can be created.
本実施形態ではこのように、r'g'b'値を介して各格子点の均等色空間上での位置を最適化することにより、S401で初期プロファイルを作成する際にユーザが調整した色分解の特性を維持したまま、各格子点の色材信号値を最適化することができる。例えば、S401でユーザが濃淡の色材（ＬcとC等）の使用割合を印字結果の粒状性等を考慮して調整することができるが、この割合についてはS403にて変更されない。したがって本実施形態では、粒状性等の画質要素を損なうことなく、格子点の均等色空間上での配置のみが最適化される。 In this embodiment, the color adjusted by the user when creating the initial profile in S401 is obtained by optimizing the position of each grid point in the uniform color space through the r'g'b 'value. The color material signal value at each grid point can be optimized while maintaining the decomposition characteristics. For example, in S401, the user can adjust the usage ratio of dark and light color materials (Lc and C, etc.) in consideration of the granularity of the print result, but this ratio is not changed in S403. Therefore, in the present embodiment, only the arrangement of the grid points on the uniform color space is optimized without impairing image quality elements such as graininess.
●補間精度評価値算出処理(S702)
以下、上記S702における補間精度評価値Ｅhの算出処理について、図８のフローチャートを用いて詳細に説明する。まずS801において補間精度評価値Ｅhを0で初期化する。そしてS802で、初期プロファイルに基づき、変数r'g'b'の格子点位置に対応するC,M,Y,K,Lc,Lmの色材信号値を求め、それぞれをC1,M1,Y1,K1,Lc1,Lm1とする。この色材信号値取得処理は、例えば図３(b)に示した初期プロファイルに対して既知の四面体補間などの補間アルゴリズムを適用することで、任意のr'g'b'に対する各色材信号値を求めることができる。次にS803では、S802で求めた色材信号値をプリンタ105にて印字した際の均等色空間上の色L*a*b*を予測し、該予測値を格子点色情報L1,a1,b1とする。具体的には、例えばS402で取得した測色情報に基づき、セル分割ユールニールセン分光ノイゲバウアモデルなどを用いて分光反射率を予測する。そして、該予測した分光反射率に対し、例えばD50等の所定の光源の分光特性、およびCIEXYZ表色系の等色関数を掛け合わせることで三刺激値XYZを取得し、これをL*a*b*値に変換する。なお、本実施形態では均等色空間としてCIEL*a*b*色空間を用い、L*a*b*値を予測する例を示したが、本発明はこの例に限らず、例えばCIE CAM97sのJab空間を用いても良い。
Interpolation accuracy evaluation value calculation processing (S702)
Hereinafter, the calculation processing of the interpolation accuracy evaluation value Eh in S702 will be described in detail with reference to the flowchart of FIG. First, in S801, the interpolation accuracy evaluation value Eh is initialized to zero. In S802, based on the initial profile, C, M, Y, K, Lc, and Lm color material signal values corresponding to the lattice point positions of the variable r'g'b 'are obtained. Let K1, Lc1, Lm1. In this color material signal value acquisition processing, for example, by applying an interpolation algorithm such as known tetrahedral interpolation to the initial profile shown in FIG. The value can be determined. Next, in S803, the color L * a * b * in the uniform color space when the color material signal value obtained in S802 is printed by the
次にS804のループ制御により、色分解プロファイルの色空間上において予め定められた全ての方向について、S805〜S811の一連の処理を行う。すなわちS805〜S811は、ある一方向に隣接する格子点を求め、その格子点情報に基づいて実行される。 Next, a series of processing from S805 to S811 is performed for all predetermined directions in the color space of the color separation profile by loop control of S804. That is, S805 to S811 are executed based on the lattice point information obtained by obtaining lattice points adjacent in a certain direction.
ここで図９を用いて、上記「予め定められた全ての方向」について説明する。図９(a)において、901は最適化の対象となっている着目格子点である。902、903は、着目格子点901に対し、R軸が増加・減少する方向にそれぞれ隣接する格子点である。同様に、904、905はG軸方向、908、909はB軸方向にそれぞれ隣接する格子点である。また、906、907はR軸方向とG軸方向に同時に増減させた斜め方向にそれぞれ隣接する格子点である。910、911は同様にR軸方向とB軸方向に同時に増減させた方向にそれぞれ隣接する格子点、912、913はG軸方向とB軸方向に同時に増減させた方向にそれぞれ隣接する格子点である。また、914、915は、R軸、G軸、B軸全ての方向で同時に増減させた方向にそれぞれ隣接する格子点である。以上説明した、着目格子点901を原点とし、該原点と902〜915の14点を結ぶ14方向を、本実施形態において「予め定められた全ての方向」とし、以下、単に「全方向」と称する。図９(b)はこの全方向を説明するための模式図であり、図９(a)からR、B平面のみを抽出したものである。このR、B平面においては、図中矢印で示した902と901、903と901、908と901、909と901、910と901、911と901のそれぞれを結ぶ6つのベクトル（方向）において、S805〜S811の処理を実行する。
Here, the “all predetermined directions” will be described with reference to FIG. In FIG. 9A,
再び図８に戻り、S805〜S811の処理を説明する。上述した全方向中の1方向について、まずS805で、当該方向に隣接する隣接格子点の色材信号値を取得する。詳細には上記S802と同様に、初期プロファイルから該当する隣接格子点のC,M,Y,K,Lc,Lmの各色材信号値を求め、それぞれを隣接格子点色材信号値C2,M2,Y2,K2,Lc2,Lm2とする。 Returning to FIG. 8 again, the processing of S805 to S811 will be described. For one direction out of all the directions described above, first, in S805, the color material signal value of the adjacent lattice point adjacent to the direction is acquired. Specifically, as in S802 above, the respective color material signal values of C, M, Y, K, Lc, and Lm of the corresponding adjacent grid points are obtained from the initial profile, and the respective color material signal values C2, M2, Y2, K2, Lc2, and Lm2.
次にS806,S807にて、第１の色値取得処理を行う。まずS806にて、隣接格子点色材信号値C2,M2,Y2,K2,Lc2,Lm2をプリンタ105にて印字した際の均等色空間上の色L*a*b*を上記S803と同様の色予測モデルを用いて予測し、該予測値を隣接格子点色情報L2,a2,b2とする。そしてS807で、着目格子点と隣接格子点との中間に位置する色情報であるL*a*b*値を、各格子点のL*a*b*値の補間演算により求める。具体的には、各格子点のL*a*b*値の平均値、すなわち上記L1,a1,b1とL2,a2,b2の平均をとり、これを格子点中間色情報L3,a3,b3とする。この格子点中間色情報L3,a3,b3はすなわち、着目格子点と隣接格子点を均等色空間上で線形補間することによって算出された第１の色値である。
Next, in S806 and S807, a first color value acquisition process is performed. First, in S806, the color L * a * b * in the uniform color space when the adjacent grid point color material signal values C2, M2, Y2, K2, Lc2, and Lm2 are printed by the
次にS808,S809にて、第２の色値取得処理を行う。すなわち、上記S806,S807とは異なる方法によって着目格子点と隣接格子点との中間に位置する色情報L*a*b*値を求める。上記S806,S807では、格子点中間色情報を各格子点のL*a*b*値の補間演算により求めたが、ここでは、各格子点の中間位置に対応する色材信号値を各格子点の色材信号値の補間演算により求め、該色材信号値に対するL*a*b*値を予測する。具体的には、まず、S808にて各格子点の色材信号値C,M,Y,K,Lc,Lmの平均値を求めて格子点中間色材信号値C3,M3,Y3,K3,Lc3,Lm3とする。次にS809で該色材信号値に対するL*a*b*値を上記S803と同様の色予測モデルを用いて予測し、該予測結果を格子点中間色情報L4,a4,b4とする。この格子点中間色情報L4,a4,b4はすなわち、着目格子点と隣接格子点で線形補間した色材信号値から予測される、均等色空間上の第２の色値である。 Next, in S808 and S809, a second color value acquisition process is performed. That is, the color information L * a * b * value located in the middle between the target lattice point and the adjacent lattice point is obtained by a method different from the above S806 and S807. In the above S806 and S807, the grid point intermediate color information is obtained by interpolation calculation of the L * a * b * values of each grid point, but here, the color material signal value corresponding to the intermediate position of each grid point is calculated for each grid point. The color material signal value is calculated by interpolation, and the L * a * b * value for the color material signal value is predicted. Specifically, first, in S808, the average value of the color material signal values C, M, Y, K, Lc, and Lm at each lattice point is obtained, and the lattice point intermediate color material signal values C3, M3, Y3, K3, and Lc3 are obtained. , Lm3. Next, in S809, the L * a * b * value for the color material signal value is predicted using the same color prediction model as that in S803, and the prediction result is set as lattice point intermediate color information L4, a4, and b4. That is, the grid point intermediate color information L4, a4, b4 is the second color value in the uniform color space predicted from the color material signal value linearly interpolated between the target grid point and the adjacent grid point.
そして最後にS810において、それぞれが異なる方法によって算出された格子点中間色情報L3,a3,b3とL4,a4,b4との色差ΔＥを、CIE1976等にしたがって算出する。そしてS811で、該算出された色差ΔＥを補間精度評価値Ｅhに加算する。ここで色差ΔＥとは、着目格子点と隣接格子点とを結ぶ線分の中間点において、線形補間により求めた色材信号値をプリンタ105に出力した場合に予測される第２の色値と、各格子点のL*a*b*値から線形補間により算出された第１の色値との差分である。
Finally, in S810, the color difference ΔE between the lattice point intermediate color information L3, a3, b3 and L4, a4, b4 calculated by different methods is calculated according to CIE1976 or the like. In step S811, the calculated color difference ΔE is added to the interpolation accuracy evaluation value Eh. Here, the color difference ΔE is the second color value predicted when the color material signal value obtained by linear interpolation is output to the
以上説明したS805〜S811の一連の処理が全方向について行われることにより、当該変数r'g'b'に対する補間精度評価値Ｅhが、全方向での色差ΔＥの総計として算出される。 By performing the series of processes of S805 to S811 described above for all directions, the interpolation accuracy evaluation value Eh for the variable r′g′b ′ is calculated as the sum of the color differences ΔE in all directions.
本実施形態では以上のように、S807にて求めたL*a*b*値（L3,a3,b3）と、S809にて求めたL*a*b*値（L4,a4,b4）との色差が小さくなるように、色分解プロファイルの各格子点を最適化する。これにより、着目格子点と各方向の隣接格子点とを結ぶ線分の中間点において、線形補間により求めた色材信号値をプリンタ105に出力した場合の色と、各格子点のL*a*b*値から線形補間により予測したL*a*b*値との色差が小さくなる。つまりは、各格子点のL*a*b*値のみから、プリンタ105にて出力される色が高精度に予測可能となる。
In the present embodiment, as described above, the L * a * b * value (L3, a3, b3) obtained in S807 and the L * a * b * value (L4, a4, b4) obtained in S809 Each grid point of the color separation profile is optimized so that the color difference of the color becomes smaller. As a result, the color when the color material signal value obtained by linear interpolation is output to the
図９を用いて説明した全方向についてS805〜S811の処理を行って作成された色分解プロファイルを用いることで、周知の四面体補間等の補間アルゴリズムによる補間をより高精度に行うことが可能となる。以下、本実施形態において格子点を最適化することで四面体補間精度が高まる例について、図１０を用いて説明する。図１０は、四面体補間に使用される格子点と線分の関係を説明する図であり、図１０(a)に示す立方体内の任意の点は、図１０(b)〜(g)に示す6つの四面体のいずれかに含まれる。四面体補間時には、これら6つの四面体から、まず対象とする点を含む1つの四面体を選択し、該選択された四面体における4頂点の格子点情報から、任意の点に対応する色情報L*a*b*値、あるいは色材信号値を決定する。したがって本実施形態のように、図１０(b)〜(g)の全ての頂点間を結ぶ直線上にて各線分の補間精度が高くなるように最適化することで、結果的に色分解プロファイルの四面体補間精度が高まることが分かる。本実施形態において図９を用いて説明した全方向は、全格子点で最適化がなされた場合に、色分解プロファイル中の全ての四面体における各辺を含んでいる。なお全方向とは、図９で説明した方向に限らず、例えば各軸の方向、すなわち±R,±G,±Bの6方向のみで構成することも可能であるし、着目格子点とその26近傍の格子点とを結ぶ全ての方向で構成することも可能である。 By using the color separation profile created by performing the processing of S805 to S811 for all directions described with reference to FIG. 9, it is possible to perform interpolation by a known interpolation algorithm such as tetrahedral interpolation with higher accuracy. Become. Hereinafter, an example in which the tetrahedral interpolation accuracy is increased by optimizing the lattice points in the present embodiment will be described with reference to FIG. FIG. 10 is a diagram for explaining the relationship between grid points and line segments used for tetrahedral interpolation. Arbitrary points in the cube shown in FIG. 10 (a) are shown in FIGS. 10 (b) to 10 (g). Included in any of the six tetrahedrons shown. At the time of tetrahedron interpolation, one tetrahedron including the target point is first selected from these six tetrahedrons, and color information corresponding to an arbitrary point is obtained from the lattice point information of the four vertices in the selected tetrahedron. The L * a * b * value or color material signal value is determined. Therefore, as in the present embodiment, optimization is performed so that the interpolation accuracy of each line segment is high on the straight line connecting all the vertices in FIGS. 10B to 10G, resulting in a color separation profile. It can be seen that the tetrahedral interpolation accuracy increases. The all directions described with reference to FIG. 9 in the present embodiment include each side in all tetrahedrons in the color separation profile when optimization is performed at all grid points. The all directions are not limited to the directions described with reference to FIG. 9, but may be configured with only the directions of the respective axes, that is, ± R, ± G, and ± B. It is also possible to configure in all directions connecting 26 neighboring grid points.
また本実施形態では、S807,S809にて格子点間の中間点の色情報L*a*b*を算出する例を示したが、格子点間の補間対象点としては中間点に限らず、他の分割点を用いても良い。例えば格子点間を3分割し、各分割位置においてS807,S809と同様にL*a*b*値を求め、それらの平均色差ΔＥを最小化するように構成しても良い。 In the present embodiment, the color information L * a * b * of the intermediate point between the lattice points is calculated in S807 and S809, but the interpolation target point between the lattice points is not limited to the intermediate point, Other division points may be used. For example, the lattice points may be divided into three, and L * a * b * values may be obtained in the same manner as in S807 and S809, and the average color difference ΔE may be minimized.
以上説明したように本実施形態によれば、補間精度評価値Ｅhを最小化するように色分解プロファイルを最適化することにより、擬似輪郭等の弊害が少ない色分解プロファイルを生成することが可能となる。したがって、より高精度に出力色を制御可能なカラーマネジメントシステムが実現される。 As described above, according to the present embodiment, by optimizing the color separation profile so as to minimize the interpolation accuracy evaluation value Eh, it is possible to generate a color separation profile with less adverse effects such as pseudo contours. Become. Therefore, a color management system capable of controlling the output color with higher accuracy is realized.
＜第２実施形態＞
以下、本発明に係る第２実施形態について説明する。上述した第１実施形態では、各格子点の補間精度評価値Ｅhが最小となるように色分解プロファイルを最適化する例を示した。第２実施形態においてはさらに階調性評価値Ｅcを設け、該Ｅcと補間精度評価値Ｅhを組み合わせた評価値を用いて最適化を行う構成例を示す。なお、第２実施形態において上述した第１実施形態と同様の構成および処理については同一符号を参照するとし、説明を省略する。
Second Embodiment
Hereinafter, a second embodiment according to the present invention will be described. In the first embodiment described above, the example in which the color separation profile is optimized so that the interpolation accuracy evaluation value Eh of each grid point is minimized has been described. In the second embodiment, a configuration example in which a gradation evaluation value Ec is further provided and optimization is performed using an evaluation value obtained by combining the Ec and the interpolation accuracy evaluation value Eh is shown. In the second embodiment, the same reference numerals are used for the same configurations and processes as those in the first embodiment described above, and the description thereof is omitted.
●格子点最適化処理(S403)
第２実施形態における色分解プロファイル作成処理も、上述した第１実施形態と同様に図４のフローチャートに従うが、S403の格子点最適化処理が第１実施形態と異なる。以下、第２実施形態における格子点最適化処理について、図１１のフローチャートを用いて詳細に説明する。なお、図１１においてS701,S702,S704〜S706の各工程は第１実施形態と同様であるため説明を省略し、以下、第１実施形態と異なる工程についてのみ説明する。
Lattice point optimization processing (S403)
The color separation profile creation processing in the second embodiment follows the flowchart of FIG. 4 as in the first embodiment described above, but the lattice point optimization processing in S403 is different from that in the first embodiment. Hereinafter, the lattice point optimization process in the second embodiment will be described in detail with reference to the flowchart of FIG. In FIG. 11, steps S701, S702, and S704 to S706 are the same as those in the first embodiment, and thus description thereof will be omitted. Hereinafter, only steps different from those in the first embodiment will be described.
S702で補間精度評価値Ｅhを算出した後、S1101において階調性評価値Ｅcを算出する。階調性評価値Ｅcは、RGB空間上で直線に並ぶ格子点が、L*a*b*空間においても直線上に並ぶように格子点を最適化するために用いられる評価値である。ここで、階調性評価値算出方法について図１２を用いて説明する。図１２(a)に示すように、まず着目格子点をＯとする。そして、RGB色空間上で着目格子点Ｏを通る直線ｄについて、該直線上で点Ｏを挟んで隣り合う格子点Ｐ，Ｍの均等色空間(L*a*b*)上での位置をそれぞれ点Ｐd，Ｍdとする。点Ｐd，Ｍdはすなわち、RGB色空間上における直線ｄの、均等色空間上での軌跡上に存在する。そして、点Ｏ，Ｐd，Ｍdのそれぞれに対するL*a*b*値を、第１実施形態におけるS802，S803と同様の方法で予測し、格子点色情報(Lo,ao,bo)，(Lp,ap,bp)，(Lm,am,bm)を得る。 After calculating the interpolation accuracy evaluation value Eh in S702, the gradation evaluation value Ec is calculated in S1101. The gradation evaluation value Ec is an evaluation value used for optimizing the lattice points so that the lattice points arranged in a straight line in the RGB space are arranged in a straight line in the L * a * b * space. Here, the gradation evaluation value calculation method will be described with reference to FIG. First, as shown in FIG. Then, for the straight line d passing through the target grid point O in the RGB color space, the positions of the grid points P and M adjacent to each other with the point O on the straight line in the uniform color space (L * a * b *) Let them be points Pd and Md, respectively. That is, the points Pd and Md exist on the locus of the straight line d on the RGB color space on the uniform color space. Then, the L * a * b * values for the points O, Pd, and Md are predicted by the same method as S802 and S803 in the first embodiment, and the lattice point color information (Lo, ao, bo), (Lp , ap, bp), (Lm, am, bm).
そして、点Ｏ，Ｐd，Ｍdの3点を結ぶ線分についての直線らしさ、すなわち該3点の並びが直線からずれる度合いを示す評価値Ｇdを下式により求める。 Then, an evaluation value Gd indicating the straightness of the line segment connecting the three points O, Pd, and Md, that is, the degree that the arrangement of the three points deviates from the straight line is obtained by the following equation.
そして下式に示すように、点Ｏを通る所定の複数方向についての評価値Ｇdの総和をとることで、着目格子点Ｏの階調性評価値Ｅcを算出する。 Then, as shown in the following equation, the gradation evaluation value Ec of the target lattice point O is calculated by taking the sum of the evaluation values Gd for a plurality of predetermined directions passing through the point O.
Ｅc＝ΣＧd
ここで所定の複数方向としては、例えば図９において以下の1〜7に示す7方向をとることができる。
Ec = ΣGd
Here, as the predetermined plural directions, for example, seven directions shown in the following 1 to 7 in FIG. 9 can be taken.
1．R軸方向：902−901−903を通る直線
2．G軸方向：904−901−905を通る直線
3．B軸方向：908−901−909を通る直線
4．R軸＋G軸方向：906−901−907を通る直線
5．R軸＋B軸方向：910−901−911を通る直線
6．G軸＋B軸方向：912−901−913を通る直線
7．R軸＋G軸＋B軸方向（明度方向）：914−901−915を通る直線
なお、階調性評価値Ｅcを算出するための複数の方向としては上記7方向に限らず、例えば上記1〜3の3方向だけであっても良いし、より多くの方向を定義しても良い。また、階調性評価値Ｅcの算出も上記式に限らず、3点の並びが直線からずれる度合いを示す評価値が得られれば良く、例えば周知の最小二乗法に基づく直線に対する誤差を評価値として適用することも可能である。
1． R axis direction: straight line passing through 902-901-903
2． G axis direction: straight line passing through 904-901-905
3． B axis direction: Straight line passing through 908-901-909
Four. R axis + G axis direction: straight line passing through 906-901-907
Five. R axis + B axis direction: straight line passing through 910-901-911
6． G axis + B axis direction: straight line passing through 912-901-913
7． R axis + G axis + B axis direction (lightness direction): straight line passing through 914-901-915 Note that the plurality of directions for calculating the gradation evaluation value Ec are not limited to the above seven directions, for example, the above 1-3 The three directions may be defined, or more directions may be defined. Further, the calculation of the gradation evaluation value Ec is not limited to the above formula, and an evaluation value indicating the degree to which the arrangement of the three points deviates from the straight line may be obtained. For example, an error with respect to a straight line based on a known least square method is evaluated. It is also possible to apply as.
第２実施形態においては、S1101で算出された階調性評価値Ｅcを最小化するように、色分解プロファイルの各格子点を最適化する。すなわちS1102において、補間精度評価値Ｅhと階調性評価値Ｅcの重み付け平均値Ｅallを算出し、最終的にS1103で該重み付け平均値Ｅallが最小もしくは十分小さくなるように、各格子点を最適化する。詳細には上述した第１実施形態と同様に、重み付け平均値Ｅallを所定の閾値（第２の閾値）と比較し、Ｅallが第２の閾値未満になるまで、着目格子点の色信号値を更新していく。 In the second embodiment, each grid point of the color separation profile is optimized so as to minimize the gradation evaluation value Ec calculated in S1101. That is, in S1102, the weighted average value Eall of the interpolation accuracy evaluation value Eh and the gradation evaluation value Ec is calculated, and finally, in S1103, each grid point is optimized so that the weighted average value Eall is minimized or sufficiently small. To do. Specifically, as in the first embodiment described above, the weighted average value Eall is compared with a predetermined threshold value (second threshold value), and the color signal value of the target lattice point is calculated until Eall is less than the second threshold value. Update.
第２実施形態ではこのように、2つの評価値ＥhとＥcの重み付け平均値Ｅallを最小化することで、補間精度と階調性のバランスのとれた色分解プロファイルが生成可能となる。なお、Ｅall算出時におけるＥh，Ｅcそれぞれに対する重みwh，wcとしては、予め設定された値であっても良いし、格子点の位置に応じて可変であっても良い。例えば、プリンタ105における色再現範囲の最外郭付近の格子点においてはＥhにかかる重みwhを大きく設定し、補間精度を重視する。一方、色再現範囲の中心付近に位置する格子点においてはＥcにかかる重みwcを大きく設定し、より階調性を重視する等、予め格子点ごとに異なる重みを設定しておいても良い。
In the second embodiment, as described above, by minimizing the weighted average value Eall of the two evaluation values Eh and Ec, it is possible to generate a color separation profile that balances interpolation accuracy and gradation. The weights wh and wc for Eh and Ec at the time of Eall calculation may be preset values or may be variable depending on the position of the lattice points. For example, at the grid points near the outermost contour of the color reproduction range in the
以上説明したように第２実施形態によれば、補間精度評価値Ｅhおよび階調性評価値Ｅcの重み付け平均値が最小となるように色分解プロファイルを最適化することにより、擬似輪郭等の弊害が少なく、かつ階調性の良好な色分解プロファイルが生成可能となる。したがって、より高精度に出力色を制御可能なカラーマネジメントシステムが実現される。 As described above, according to the second embodiment, the color separation profile is optimized so that the weighted average value of the interpolation accuracy evaluation value Eh and the gradation evaluation value Ec is minimized, thereby causing adverse effects such as pseudo contours. And a color separation profile with good gradation can be generated. Therefore, a color management system capable of controlling the output color with higher accuracy is realized.
＜その他の実施形態＞
上述した第１および第２実施形態においては、C,M,Y,K,Lc,Lmからなる６色の色材を有するプリンタを用いて説明したが、本発明の色処理装置において想定される色材数および種類はこの例に限定されない。
<Other embodiments>
In the first and second embodiments described above, the description has been given using the printer having six color materials composed of C, M, Y, K, Lc, and Lm. However, it is assumed in the color processing apparatus of the present invention. The number and type of color materials are not limited to this example.
また、上述した第１および第２実施形態においては色分解プロファイル作成アプリケーションに本発明を適用した例を説明したが、本発明はこれに限るものではない。例えばプリンタ本体に内蔵されたプロファイル作成装置として本発明を機能させることも可能である。 In the first and second embodiments described above, the example in which the present invention is applied to the color separation profile creation application has been described. However, the present invention is not limited to this. For example, it is possible to make the present invention function as a profile creation device built in the printer body.
また、上述した第１および第２実施形態においては、均等色空間上での色差が十分に小さくなるように色分解プロファイルの各格子点を最適化する例を示したが、例えばC,M,Y,K,Lc,Lmの各色材信号値を直接最適化することも可能である。 In the first and second embodiments described above, an example in which each lattice point of the color separation profile is optimized so that the color difference in the uniform color space is sufficiently small has been described. For example, C, M, It is also possible to directly optimize the color material signal values of Y, K, Lc, and Lm.
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（またはＣＰＵやＭＰＵ等）がプログラムを読み出して実行する処理である。 The present invention can also be realized by executing the following processing. That is, software (program) that realizes the functions of the above-described embodiments is supplied to a system or apparatus via a network or various storage media, and a computer (or CPU, MPU, or the like) of the system or apparatus reads the program. It is a process to be executed.
Claims (8)
前記色空間における着目格子点の色信号値、および、該着目格子点に対して予め設定された複数方向に隣接する複数の隣接格子点の色信号値に対応する複数の色材信号値を取得する信号値取得手段と、
前記取得した複数の色材信号値から、色予測モデルを用いて予測される均等色空間上の複数の色値を用いた補間により得られる、前記着目格子点と前記隣接格子点との間の複数の補間対象点それぞれに対する第一の色値を取得する第一の取得手段と、
前記取得した複数の色材信号値を用いた補間から得られる色材信号値から、前記複数の補間対象点それぞれに対する第二の色値を取得する第二の取得手段と、
前記複数の補間対象点に対する、前記第一の色値と前記第二の色値との色差から前記複数の補間対象点における補間精度を評価する評価値を算出し、前記評価値に基づき前記着目格子点の色材信号値を設定する設定手段とを有することを特徴とする色処理装置。 A color processing device that creates a color separation profile indicating a relationship between a color signal value of a grid point in a color space of an input image and a color material signal value of an image forming device,
A color signal value of a target grid point in the color space and a plurality of color material signal values corresponding to color signal values of a plurality of adjacent grid points adjacent to the target grid point in a plurality of directions set in advance are acquired. Signal value acquisition means for
From the acquired plurality of color material signal values, obtained by interpolation using a plurality of color values on a uniform color space predicted using a color prediction model, between the grid point of interest and the adjacent grid point First acquisition means for acquiring a first color value for each of a plurality of interpolation target points;
A second acquisition means for acquiring a second color value for each of the plurality of interpolation target points from a color material signal value obtained by interpolation using the acquired plurality of color material signal values;
An evaluation value for evaluating the interpolation accuracy at the plurality of interpolation target points is calculated from a color difference between the first color value and the second color value for the plurality of interpolation target points , and the focus is calculated based on the evaluation value A color processing apparatus comprising: setting means for setting color material signal values of lattice points.
前記着目格子点の色信号値を設定し、前記設定した色信号値について、前記信号値取得手段、並びに、前記第一および第二の取得手段を用いて前記評価値の算出を繰り返す手段と、Setting the color signal value of the grid point of interest, and for the set color signal value, repeating the calculation of the evaluation value using the signal value acquisition means, and the first and second acquisition means;
前記繰り返しにより所定の閾値未満の評価値が算出されると、当該算出において設定された色信号値に対応する色材信号値を、当該算出において前記信号値取得手段が取得した複数の色材信号値の補間によって算出し、前記算出した色材信号値を前記繰り返し前の前記着目格子点の色信号値に対応する色材信号値として設定する手段とを有することを特徴とする請求項1に記載された色処理装置。When an evaluation value less than a predetermined threshold is calculated by the repetition, a plurality of color material signals acquired by the signal value acquisition unit in the calculation are obtained as color material signal values corresponding to the color signal values set in the calculation. 2. The apparatus according to claim 1, further comprising means for calculating by color value interpolation and setting the calculated color material signal value as a color material signal value corresponding to the color signal value of the target grid point before the repetition. The described color processing apparatus.
前記色予測モデルは、前記画像形成装置が前記色材信号値に応じて画像形成を行った際に得られる前記均等色空間上の色値を、前記測色値を用いて予測するためのモデルであることを特徴とする請求項2に記載された色処理装置。 And an input unit for inputting a colorimetric value of a color chart formed by the image forming apparatus with respect to a lattice point in a color space represented by the color material signal value,
Model for the color prediction model, the color values on the uniform color space obtained when the image forming apparatus performs an image forming in accordance with the color material signal values predicted using the colorimetric values 3. The color processing apparatus according to claim 2, wherein
前記作成手段は、前記格子点ごとの色材信号値をユーザ指示に応じて設定することを特徴とする請求項1から請求項4の何れか一項に記載された色処理装置。 In addition, for each grid point in the color space of the input image, there is a creation means for creating an initial profile in which initial values of color material signal values for the color signal values are set,
5. The color processing apparatus according to claim 1, wherein the creating unit sets a color material signal value for each grid point in accordance with a user instruction.
前記着目格子点、および、該着目格子点を通る直線上で該着目格子点を挟んで隣接する複数の隣接格子点それぞれの色信号値に対応する複数の色材信号値を取得し、前記複数の色材信号値に対応する前記均等色空間上の複数の色値を予測し、前記複数の色値の並びが前記均等色空間上の直線からずれる度合いを、該着目格子点に対する階調性の評価値として算出する評価手段と、
前記着目格子点に対する前記色差と前記評価値の重み付け平均値を算出する平均値算出手段と、
前記着目格子点の色信号値を設定し、前記設定した色信号値について、前記信号値取得手段と前記第一および第二の取得手段による前記色差の算出、前記評価手段による評価値の算出、並びに、前記平均値算出手段による重み付け平均値の算出を繰り返す手段と、
前記繰り返しにより所定の閾値未満の重み付け平均値が算出されると、当該算出において設定された色信号値に対応する色材信号値を、当該算出において前記信号値取得手段が取得した複数の色材信号値の補間によって算出し、前記算出した色材信号値を前記繰り返し前の前記着目格子点の色信号値に対応する色材信号値として設定する手段とを有することを特徴とする請求項1に記載された色処理装置。 The setting means includes
The grid point of interest, and to obtain a plurality of color material signal values corresponding to a plurality of adjacent grid points each color signal values adjacent to each other across the remarked lattice points on a straight line passing through the remarked grid points, said plurality A plurality of color values in the uniform color space corresponding to the color material signal value of the color material are predicted, and the degree of gradation with respect to the grid point of interest is determined based on the degree of deviation of the arrangement of the plurality of color values from the straight line in the uniform color space . An evaluation means for calculating the evaluation value of
An average value calculating means for calculating a weighted average value of the color difference and the evaluation value for the grid point of interest;
Set the color signal value of the grid point of interest, and for the set color signal value, the calculation of the color difference by the signal value acquisition means and the first and second acquisition means, the calculation of the evaluation value by the evaluation means , and, means for to repeat the calculation of the weighted average value by the average value calculating means,
When the weighted average value less than the predetermined threshold is calculated by the repetition, the color material signal values corresponding to the color signal values set in the calculation are obtained by the signal value acquisition unit in the calculation. and characterized in that it has means for calculating I by the interpolated signal value, sets the colorant signal value the calculated as the colorant signal value that corresponds to the color signal values of the grid point of interest of the repeated before The color processing apparatus according to claim 1 .
前記色空間における着目格子点の色信号値、および、該着目格子点に対して予め設定された複数方向に隣接する複数の隣接格子点の色信号値に対応する複数の色材信号値を取得し、
前記取得した複数の色材信号値から、色予測モデルを用いて予測される均等色空間上の複数の色値を用いた補間により得られる、前記着目格子点と前記隣接格子点との間の複数の補間対象点それぞれに対する第一の色値を取得し、
前記取得した複数の色材信号値を用いた補間から得られる色材信号値から、前記複数の補間対象点それぞれに対する第二の色値を取得し、
前記複数の補間対象点に対する、前記第一の色値と前記第二の色値との色差から前記複数の補間対象点における補間精度を評価する評価値を算出し、前記評価値に基づき前記着目格子点の色材信号値を設定することを特徴とする色処理方法。 A color processing method for creating a color separation profile indicating a relationship between a color signal value of a grid point in a color space of an input image and a color material signal value of an image forming apparatus,
A color signal value of a target grid point in the color space and a plurality of color material signal values corresponding to color signal values of a plurality of adjacent grid points adjacent to the target grid point in a plurality of directions set in advance are acquired. And
From the acquired plurality of color material signal values, obtained by interpolation using a plurality of color values on a uniform color space predicted using a color prediction model, between the grid point of interest and the adjacent grid point Obtain the first color value for each of multiple interpolation target points,
From the color material signal value obtained from the interpolation using the obtained plurality of color material signal values, to obtain a second color value for each of the plurality of interpolation target points,
An evaluation value for evaluating the interpolation accuracy at the plurality of interpolation target points is calculated from a color difference between the first color value and the second color value for the plurality of interpolation target points , and the focus is calculated based on the evaluation value A color processing method characterized by setting color material signal values of lattice points.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2010201340A JP5665442B2 (en) | 2010-09-08 | 2010-09-08 | Color processing apparatus and color processing method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2010201340A JP5665442B2 (en) | 2010-09-08 | 2010-09-08 | Color processing apparatus and color processing method |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JP2012060397A JP2012060397A (en) | 2012-03-22 |
| JP2012060397A5 JP2012060397A5 (en) | 2013-10-24 |
| JP5665442B2 true JP5665442B2 (en) | 2015-02-04 |
Family
ID=46056959
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2010201340A Expired - Fee Related JP5665442B2 (en) | 2010-09-08 | 2010-09-08 | Color processing apparatus and color processing method |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP5665442B2 (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6349747B2 (en) * | 2014-01-30 | 2018-07-04 | セイコーエプソン株式会社 | Color conversion lookup table generation method, apparatus, and program |
| JP2020102748A (en) | 2018-12-21 | 2020-07-02 | セイコーエプソン株式会社 | Color conversion table creation method, color conversion table creation device, and computer program |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH08287226A (en) * | 1995-04-18 | 1996-11-01 | Ricoh Co Ltd | Hierarchical optimizing method for color correction table |
| JPH099080A (en) * | 1995-06-16 | 1997-01-10 | Canon Inc | Method for preparing transformation table and color transforming device |
| JP3733200B2 (en) * | 1997-05-12 | 2006-01-11 | キヤノン株式会社 | Image processing apparatus and image processing method |
| JP2003348365A (en) * | 2002-05-30 | 2003-12-05 | Matsushita Electric Ind Co Ltd | Apparatus and method for color conversion |
| JP4771538B2 (en) * | 2006-07-20 | 2011-09-14 | キヤノン株式会社 | Color conversion table generation method, color conversion table, and color conversion table generation apparatus |
| JP2009182637A (en) * | 2008-01-30 | 2009-08-13 | Canon Inc | Color processor and its method |
-
2010
- 2010-09-08 JP JP2010201340A patent/JP5665442B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2012060397A (en) | 2012-03-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5381380B2 (en) | Image processing apparatus, image forming system, and program | |
| JP4595734B2 (en) | Profile creation method, profile creation device, profile creation program, print control method, print control device, and print control program | |
| JP4845127B2 (en) | Image processing device | |
| US8335015B2 (en) | Systems and methods for building a color lookup table for a printer | |
| US8593692B2 (en) | Systems and methods for building a color lookup table for a printer | |
| JPH05191641A (en) | Picture processing method and its device | |
| US20090059326A1 (en) | Color space conversion apparatus and method for controlling the same | |
| US8531729B2 (en) | Color processing apparatus and method thereof | |
| JP2008022140A (en) | Device and method for preparing color-conversion coefficient and color-conversion coefficient preparation program | |
| US9396419B2 (en) | Data-processing apparatus generating color conversion data | |
| JP5777497B2 (en) | Color processing apparatus and method | |
| JP2008148275A (en) | Color processing apparatus and program | |
| JP5324405B2 (en) | Color processing apparatus and method, and image forming apparatus | |
| JP5440195B2 (en) | Color profile creation method, image processing apparatus for creating color profile, and control program for image processing apparatus | |
| JP5665442B2 (en) | Color processing apparatus and color processing method | |
| JP5777322B2 (en) | Color processing apparatus and color processing method | |
| US8564830B2 (en) | Sensitivity matrix determination via chain rule of differentiation | |
| JP4159912B2 (en) | Color conversion method with hue correction using multiple lookup tables and interpolation | |
| JP2012129912A (en) | Printer | |
| US20200329174A1 (en) | High precision gamut mapping | |
| JP4139590B2 (en) | CMYK color printer black generation method using multiple look-up tables and interpolation method | |
| JP5306110B2 (en) | Image processing apparatus, image processing method, program, and storage medium | |
| JP2009296139A (en) | Image processor, image processing method and computer program | |
| JP5235817B2 (en) | Image processing apparatus, image processing method, and computer program | |
| JP2024042649A (en) | Image processing device, image processing method and program |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20130909 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20130909 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20140728 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140801 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140929 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20141110 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20141209 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 5665442 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
| LAPS | Cancellation because of no payment of annual fees |