JP2011223434A - Color conversion table creation method, color conversion table creation program, and printer - Google Patents
Color conversion table creation method, color conversion table creation program, and printer Download PDFInfo
- Publication number
- JP2011223434A JP2011223434A JP2010092027A JP2010092027A JP2011223434A JP 2011223434 A JP2011223434 A JP 2011223434A JP 2010092027 A JP2010092027 A JP 2010092027A JP 2010092027 A JP2010092027 A JP 2010092027A JP 2011223434 A JP2011223434 A JP 2011223434A
- Authority
- JP
- Japan
- Prior art keywords
- ink
- color
- ink amount
- conversion table
- value
- 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.)
- Withdrawn
Links
Images
Landscapes
- Facsimile Image Signal Circuits (AREA)
- Color Image Communication Systems (AREA)
Abstract
Description
この発明は、色変換テーブルを作成する技術に関し、特に、入力表色系の座標値を複数種類のインクで構成されるインク表色系のインク量に変換するための色変換テーブルを作成する技術に関する。 The present invention relates to a technique for creating a color conversion table, and in particular, a technique for creating a color conversion table for converting an input color system coordinate value into an ink color system ink amount composed of a plurality of types of ink. About.
近年、普及しているカラーインクジェットプリンタ等の印刷装置では、シアン（Ｃ）、マゼンタ（Ｍ）、イエロー（Ｙ）、ブラック（Ｋ）の各インクに加えて、色の再現範囲を広げるために特色インクと呼ばれるオレンジ（Ｏｒ）、グリーン（Ｇｒ）、ブルー（Ｂ）、レッド（Ｒ）、バイオレット（Ｖ）等のインクを搭載するものもある。これらは、シアン（Ｃ）、マゼンタ（Ｍ）、イエロー（Ｙ）が有彩１次色インクと呼ばれ、Ｏｒ、Ｇｒ、Ｂ、Ｒ、Ｖ等の特色インクが有彩２次色インクと呼ばれる。有彩２次色とは、２つの有彩１次色成分に分解できる色を意味する。 In recent years, printing apparatuses such as color ink jet printers, which are widely used, feature special colors to expand the color reproduction range in addition to cyan (C), magenta (M), yellow (Y), and black (K) inks. Some inks such as orange (Or), green (Gr), blue (B), red (R), and violet (V), which are called inks, are mounted. Of these, cyan (C), magenta (M), and yellow (Y) are called chromatic primary color inks, and special color inks such as Or, Gr, B, R, and V are called chromatic secondary color inks. . A chromatic secondary color means a color that can be separated into two chromatic primary color components.
印刷装置では、カラー画像の任意の色に応じて使用可能なインクのインク量を決定するが、このような色再現のために印刷時に用いる各インクのインク量を決定する分版処理を行う。カラー画像の色データと各色インク量との関係は、一般に予め色変換テーブル（色変換ルックアップテーブル（LUT:Look Up Table））として記憶されており、印刷時には、色変換テーブルにしたがって各画素位置における各色のインク量を決定する。この分版処理は、インク色分解処理とも呼ばれる。特許文献１には、一次色（ＣＭＹ）で構成された入力色を、１次色インクと２次色インクとを含んで構成されたインク量セットに分版する分版処理について、記載されている。
In the printing apparatus, the ink amount of ink that can be used is determined according to an arbitrary color of the color image, and for such color reproduction, a color separation process is performed to determine the ink amount of each ink used during printing. The relationship between the color data of the color image and the amount of each color ink is generally stored in advance as a color conversion table (LUT: Look Up Table), and each pixel position is determined according to the color conversion table during printing. The amount of ink for each color is determined. This separation processing is also called ink color separation processing.
上述したように、特許文献１では入力色をＣＭＹで表現される１次色表色系へと変換し、このＣＭＹで表現される色をＣＭＹＲＶで表現される再現色表色系に分版している。この分版を行う際に、特許文献１では、粒状性画質の向上を行うために、インクの置換条件を詳細に定義し、各定義、条件による処理の分岐を多用するなどの複雑な処理を実施している。
As described above, in
本発明は、入力表色系の座標値を特色インクを含むインク量セットへ分版する際に、より簡易な方法でインクの発生を制御し、効果的に粒状性画質を向上させた印刷結果を得られるように色対応関係を規定した色変換テーブルを作成することが可能な色変換テーブル作成方法、色変換テーブル作成プログラムおよび印刷装置の提供を目的とする。 The present invention controls the generation of ink by a simpler method when separating the coordinate values of the input color system into an ink amount set including special color ink, and effectively improves the graininess image quality. An object of the present invention is to provide a color conversion table creation method, a color conversion table creation program, and a printing apparatus that can create a color conversion table that defines a color correspondence relationship so as to obtain a color conversion table.
上記課題の少なくとも一部を解決するために、本発明にかかる色変換テーブル作成方法は、入力表色系の座標値を複数種類のインクで構成されるインク表色系のインク量セットに変換するための色変換テーブル作成方法において、上記入力表色系の座標値に対応する上記インク表色系のインク量を、所定の目的関数を用いた最適化によって決定する最適化工程と、上記複数種類のインクの中から選択した３つのインクを特定の階調値で印刷して測色した測色値の発色特性を示すベクトルを基底として構成された仮想色空間へ上記複数種類のインクのインク量をそれぞれ変換するための変換比率を示す置換比率ベクトルの中から選択された複数の置換比率ベクトルを線形結合して各単色インクを表したときの結合係数に基づいて作成された上記インク量セットの発生の制限条件を満たすインク量セットが上記最適化工程において最適化対象となるように制御するインク量セット制御工程と、上記最適化工程によって決定されたインク量に基づいて、上記入力表色系の座標値を上記インク表色系のインク量に変換するための色変換テーブルを作成する色変換テーブル作成工程と、を備える構成としてある。 In order to solve at least a part of the above problems, a color conversion table creation method according to the present invention converts an input color coordinate system coordinate value into an ink color set ink amount set composed of a plurality of types of ink. In the color conversion table creating method, an optimization step for determining the ink amount of the ink color system corresponding to the coordinate value of the input color system by optimization using a predetermined objective function, and the plurality of types Ink amounts of the plurality of types of inks to a virtual color space configured based on a vector indicating color development characteristics of colorimetric values obtained by printing three inks selected from the above-described inks with specific gradation values A plurality of replacement ratio vectors selected from replacement ratio vectors indicating conversion ratios for converting each of the above are created based on a coupling coefficient when each single color ink is represented by linearly combining the replacement ratio vectors. Based on the ink amount set control step for controlling the ink amount set that satisfies the restriction condition for the generation of the ink amount set to be the optimization target in the optimization step, and the ink amount determined by the optimization step, A color conversion table creating step for creating a color conversion table for converting the coordinate value of the input color system into the ink amount of the ink color system.
このように構成すると、最適化工程の中でインクの発生を制御することが可能であり、最適化工程の対象となるインク量セットを粒状性が悪化しないようなインク量セットに限定することが可能になり、上記入力表色系の座標値に対応づけられるインク量セットの粒状性が良好になるように制御することが可能となる。すなわち、簡易な方法で、効果的に粒状性画質を向上した印刷結果を得られるような色変換が可能な色変換テーブルを作成することができる。 With this configuration, it is possible to control the generation of ink during the optimization process, and it is possible to limit the ink quantity set that is the target of the optimization process to an ink quantity set that does not deteriorate the graininess. It becomes possible, and it becomes possible to control so that the granularity of the ink amount set associated with the coordinate value of the input color system is good. In other words, it is possible to create a color conversion table capable of color conversion so as to obtain a printing result with improved granularity image quality effectively by a simple method.
また、本発明の選択的な一態様として、上記インク量セット制御工程では、上記複数種類のインクのうち、濃色インクについて上記制限条件を適用する構成としてもよい。このように構成すると、粒状性を悪化させやすい濃色インクが他のインクで置換えられたインク量セットが最適化工程の対象となりやすくなり、粒状性を効果的に改善することができる。 Further, as a selective aspect of the present invention, in the ink amount set control step, the restriction condition may be applied to dark ink among the plurality of types of ink. With such a configuration, an ink amount set in which dark ink that easily deteriorates graininess is replaced with another ink is likely to be an object of the optimization process, and graininess can be effectively improved.
また、本発明の選択的な一態様として、上記インク量セット制御工程では、上記複数種類のインクのうち、淡色インクについては上記制限条件を適用しない構成としてもよい。このように構成すると、粒状性を悪化させにくい淡色インクが制限されずに最適化対象となるため、淡色インクを含むインク量セットが最適化工程の対象となりやすくなり、粒状性を効果的に改善することができる。 Further, as a selective aspect of the present invention, the ink amount set control step may be configured such that the restriction condition is not applied to the light color ink among the plurality of types of ink. With this configuration, light color inks that do not easily deteriorate graininess are subject to optimization without being restricted, so an ink amount set that includes light color inks is likely to be subject to optimization processes, effectively improving graininess. can do.
また、本発明の選択的な一態様として、上記インク量セット制御工程では、上記複数種類のインクのうち、特色インクについては上記制限条件を適用しない構成としてもよい。このように構成すると、特色インクについては制限されずに最適化対象となるため、複数種類のインクによって表現可能な色再現域をいたずらに抑制することが無くなる。よって、粒状性を良好に保つための制限条件の中で色再現も良好に保つことが可能となる。 Further, as a selective aspect of the present invention, in the ink amount setting control step, the restriction condition may not be applied to the special color ink among the plurality of types of ink. With this configuration, since the special color ink is not limited and is subject to optimization, the color reproduction range that can be expressed by a plurality of types of ink is not unnecessarily suppressed. Therefore, it is possible to maintain good color reproduction within the limiting conditions for maintaining good graininess.
また、本発明の選択的な一態様として、上記インク量セット制御工程では、インク量セットを構成する濃色インクを他のインクで置換可能なときは該インク量セットが発生しないように上記最適化工程を制御する構成としてもよい。このように構成すると、淡色インクで再現可能な上記入力表色系の座標値は優先的に淡色インクのインク量セットが対応付けられることになる。逆に、淡色インクで置換えられない量の濃色インクを含むインク量セットは、上記入力表色系の座標値に対応付けられるようになる。よって、粒状性画質を効果的に向上した印刷結果を得られるような色変換が可能な色変換テーブルを作成することができる。 Further, as a selective aspect of the present invention, in the ink amount set control step, when the dark color ink constituting the ink amount set can be replaced with another ink, the optimal ink amount setting is not generated. It is good also as a structure which controls a formation process. With this configuration, the coordinate value of the input color system that can be reproduced with the light color ink is preferentially associated with the ink amount set of the light color ink. Conversely, an ink amount set that includes an amount of dark color ink that cannot be replaced with light color ink is associated with the coordinate values of the input color system. Therefore, it is possible to create a color conversion table capable of color conversion so as to obtain a printing result that effectively improves the graininess image quality.
なお、本発明は、種々の形態で実現することが可能であり、例えば、色変換テーブル作成装置、そのための平滑化／最適化処理方法および装置、色変換テーブルを印刷装置に組み込む印刷装置の製造方法および製造システム、それらの方法または装置の機能を実現するためのコンピュータプログラム、そのコンピュータプログラムを記録した記録媒体等の形態で実現することができる。さらに、色変換テーブルが組み込まれた印刷装置においても本発明の有用性を実現することができる。 The present invention can be realized in various forms. For example, a color conversion table creation apparatus, a smoothing / optimization method and apparatus therefor, and a printing apparatus that incorporates the color conversion table into the printing apparatus. The present invention can be realized in the form of a method and manufacturing system, a computer program for realizing the functions of the method or apparatus, a recording medium on which the computer program is recorded, and the like. Furthermore, the usefulness of the present invention can be realized also in a printing apparatus in which a color conversion table is incorporated.
次に、本発明の実施の形態を以下の順序で説明する。
（１）装置構成と全体処理手順：
（２）置換比率マトリックスの作成手順：
（３）インク発生点の制限条件の作成手順：
（４）力学モデル：
（５）スムージング処理（平滑化／最適化処理）の処理手順：
（６）最適化処理の内容：
（７）変形例：
Next, embodiments of the present invention will be described in the following order.
(1) Device configuration and overall processing procedure:
(2) Procedure for creating a substitution ratio matrix:
(3) Procedure for creating ink generation point restriction conditions:
(4) Mechanical model:
(5) Smoothing processing (smoothing / optimization processing) processing procedure:
(6) Details of optimization processing:
(7) Modification:
（１）装置構成と全体処理手順：
図１は、本発明の一実施例におけるルックアップテーブル作成装置の構成を示すブロック図である。この装置は、ベースＬＵＴ作成モジュール１００と、色補正ＬＵＴ作成モジュール２００と、コンバーター３００と、ＬＵＴ格納部４００と、を備えている。「ＬＵＴ」は、色変換プロファイルとしてのルックアップテーブルの略語である。これらのモジュール１００，２００やコンバーター３００の機能は、メモリに格納されたコンピュータプログラムをコンピューターが実行することによってそれぞれ実現される。また、ＬＵＴ格納部４００は、ハードディスク装置などの記録媒体によって実現される。
(1) Device configuration and overall processing procedure:
FIG. 1 is a block diagram showing the configuration of a lookup table creation apparatus in an embodiment of the present invention. This apparatus includes a base
ベースＬＵＴ作成モジュール１００は、スムージング処理初期値設定モジュール１２０と、テーブル作成モジュール１４０とを有している。スムージング処理モジュール１３０は、色点移動モジュール１３２と、インク量最適化モジュール１３４と、画質評価指数算出モジュール１３６とを有している。コンバーター３００は、後述のコンバーター３００に基づいてインク量セットを仮想ＣＭＹデータに変換する。これらの各部の機能については後述する。
The base
ＬＵＴ格納部４００は、インバースモデル初期ＬＵＴ４１０や、ベース３Ｄ−ＬＵＴ５１０，ベース４Ｄ−ＬＵＴ５２０，色補正３Ｄ−ＬＵＴ６１０，色補正４Ｄ−ＬＵＴ６２０などを格納するためのものである。但し、インバースモデル初期ＬＵＴ４１０以外のＬＵＴは、ベースＬＵＴ作成モジュール１００や色補正ＬＵＴ作成モジュール２００によって作成されるものである。
ベース３Ｄ−ＬＵＴ５１０は、ＲＧＢ表色系を入力とし、インク量を出力とする色変換ルックアップテーブルである。一方、ベース４Ｄ−ＬＵＴ５２０は、ＣＭＹＫ表色系を入力とし、インク量を出力とする色変換ルックアップテーブルである。なお、「３Ｄ」や「４Ｄ」は、入力値の数を意味している。これらのベースＬＵＴ５１０，５２０は、例えば色補正ＬＵＴ６１０，６２０を作成する際に使用される。「ベースＬＵＴ」という名前は、色補正ＬＵＴを作成する基礎として用いられるからである。
色補正ＬＵＴ６１０，６２０は、標準的な機器依存表色系（例えばｓＲＧＢ表色系やＪＡＰＡＮ ＣＯＬＯＲ ２００１表色系）を、特定のプリンターのインク量に変換するためのルックアップテーブルである。インバースモデル初期ＬＵＴ４１０については後述する。
The
The base 3D-
The
図２は、実施例の全体処理手順を示すフローチャートである。図３（Ａ）〜（Ｃ）は、図２のステップＳ１００〜Ｓ３００によってベース３Ｄ−ＬＵＴを作成する場合の処理内容を示す説明図である。ステップＳ１００では、置換比率マトリックス３１０とインバースモデル初期ＬＵＴ４１０とが準備される。ここで、「インバースモデル」とは、後述する仮想空間上の値をインク量に変換する変換モデルを意味する。実施例では、機器非依存表色系としてＣＩＥ−Ｌａｂ表色系を使用する。なお、以下では、ＣＩＥ−Ｌａｂ表色系の色彩値を、単に「Ｌ*ａ*ｂ*値」または「Ｌａｂ値」とも呼ぶ。
FIG. 2 is a flowchart showing an overall processing procedure of the embodiment. 3A to 3C are explanatory diagrams showing processing contents when a base 3D-LUT is created by steps S100 to S300 of FIG. In step S100, a
図３（Ａ）に示すように、コンバーター３００は、置換比率マトリックス３１０を参照して、複数種類のインクのインク量を、仮想的な色空間である仮想ＣＭＹ空間の色点であるＶC，ＶM，ＶYに変換する。以下、仮想ＣＭＹ空間の色点ＶC，ＶM，ＶYを仮想ＣＭＹ値と呼ぶことにする。本実施例では、シアン（Ｃ）、マゼンタ（Ｍ）、イエロー（Ｙ）、ブラック（Ｋ）、淡シアン（Ｌｃ）、淡マゼンタ（Ｌｍ）、淡ブラック（Ｌｋ）、極淡ブラック（ＬＬｋ）、オレンジ（Ｏｒ）、グリーン（Ｇｒ）の１０種類のインクを利用可能なカラープリンタを想定しており、コンバーター３００もこの１０種類のインクのインク量を仮想ＣＭＹ空間の色点に変換する。但し、プリンターで使用する複数種類のインクとしては、任意のインクセットを利用することが可能である。なお、コンバーター３００の出力となる仮想的な色空間には、ＣＭＹ以外にもＬｃＬｍＹ等で構成した仮想的な色空間を採用してもよく、プリンターに搭載される複数のインクの中から任意に選択された３つのインク種により構成した仮想的な色空間を利用することができる。仮想的な色空間の構築については、後に詳述する。
As shown in FIG. 3A, the
インバースモデル初期ＬＵＴ４１０は、仮想ＣＭＹ値を入力とし、インク量を出力とするルックアップテーブルである。この初期ＬＵＴ４１０は、例えば、仮想ＣＭＹ空間を複数の小セルに区分し、各小セル毎に最適なインク量を選択して登録したものである。この選択は、例えば、そのインク量で印刷されるカラーパッチの画質を考慮して行われる。一般に、或る１つの仮想ＣＭＹ値を再現するインク量の組合せは多数存在する。そこで、初期ＬＵＴ４１０では、ほぼ同じ仮想ＣＭＹ値を再現する多数のインク量の組合せの中から、画質等の所望の観点から最適なインク量を選択したものが登録されている。この初期ＬＵＴ４１０の入力値である仮想ＣＭＹ値は各小セルの代表値である。一方、出力値であるインク量はそのセル内のいずれかの仮想ＣＭＹ値を再現するものである。従って、この初期ＬＵＴ４１０では、入力値である仮想ＣＭＹ値と出力値であるインク量とが厳密に対応したものとなっておらず、出力値のインク量をコンバーター３００で仮想ＣＭＹ値に変換すると、初期ＬＵＴ４１０の入力値とは多少異なる値が得られる。但し、初期ＬＵＴ４１０として、入力値と出力値とが完全に対応するものを利用してもよい。また、初期ＬＵＴ４１０を用いずにベースＬＵＴを作成することも可能である。なお、小セル毎に最適なインク量を選択して初期ＬＵＴ４１０を作成する方法としては、各仮想ＣＭＹ空間の小セルに対応するＬ*ａ*ｂ*空間の小セルについて、例えば前記特表２００７−５１１１７５号公報に記載された方法を採用して選択したインク量を、仮想ＣＭＹ空間の各小セルに対応づけることにより作成することが可能である。なお、インバースモデル初期ＬＵＴへ登録するインク量についても、後述の置換比率マトリックスから作成される発生点制御条件を適用しつつ作成されていてもよい。
The inverse model
図２のステップＳ２００では、ベースＬＵＴ作成のための初期入力値がユーザによって設定される。図３（Ｂ）は、ベース３Ｄ−ＬＵＴ５１０の構成とその初期入力値設定の例を示している。ベース３Ｄ−ＬＵＴ５１０の入力値としては、ＲＧＢの各値として予め定められたほぼ等間隔の値が設定される。１組のＲＧＢ値はＲＧＢ色空間内の点を表していると考えられるので、１組のＲＧＢ値を「入力格子点」とも呼ぶ。ステップＳ２００においては、複数の入力格子点のうちから予め選択されたいくつかの少数の入力格子点に対するインク量の初期値がユーザによって入力される。この初期入力値が設定される入力格子点としては、ＲＧＢ色空間における３次元色立体の頂点に相当する入力格子点を少なくとも選択することが好ましい。この３次元色立体の頂点では、ＲＧＢの各値がその定義範囲の最小値または最大値を取る。具体的には、ＲＧＢの各値を８ビットで表現した場合には、（Ｒ，Ｇ，Ｂ）＝（０，０，０）、（０，０，２５５）、（０，２５５，０）、（２５５，０，０）、（０，２５５，２５５）、（２５５，０，２５５）、（２５５，２５５，０）、（２５５，２５５，２５５）である８つの入力格子点に関してインク量の初期入力値が設定される。なお、（Ｒ，Ｇ，Ｂ）＝（２５５，２５５，２５５）の入力格子点に対するインク量は、すべて０に設定される。他の入力格子点に対するインク量の初期入力値は任意であり、例えば０に設定される。図３（Ｂ）の例では、（Ｒ，Ｇ，Ｂ）＝（０，０，３２）の入力格子点に対するインク量が０以外の値になっているが、これはこのＬＵＴ５１０が完成したときの値である。
In step S200 of FIG. 2, an initial input value for creating a base LUT is set by the user. FIG. 3B shows an example of the configuration of the base 3D-
図２のステップＳ３００では、スムージング処理モジュール１３０（図１）が、ステップＳ２００で設定された初期入力値に基づいてスムージング処理（平滑化／最適化処理）を実行する。図３（Ｃ）は、ステップＳ３００の処理内容を示している。図３（Ｃ）の左側には、スムージング処理前の状態における複数の色点の分布が２重丸と白丸とで示されている。これらの色点は、仮想ＣＭＹ空間における３次元色立体ＣＳを構成している。なお、この３次元色立体ＣＳはＬａｂ空間の３次元色立体と相関しており、同図には、仮想ＣＭＹ空間に対応するＬａｂ空間の各軸も併記してある。各色点の仮想ＣＭＹ座標値は、ベース３Ｄ−ＬＵＴ５１０の複数の入力格子点におけるインク量を、コンバーター３００（図３（Ａ））を用いて仮想ＣＭＹ値に変換した値である。上述したように、ステップＳ２００では一部の少数の入力格子点についてのみインク量の初期入力値が設定される。そこで、他の入力格子点に対するインク量の初期値は、初期入力値からスムージング処理初期値設定モジュール１２０（図１）によって設定される。この初期値設定方法については後述する。
In step S300 of FIG. 2, the smoothing processing module 130 (FIG. 1) executes smoothing processing (smoothing / optimization processing) based on the initial input value set in step S200. FIG. 3C shows the processing content of step S300. On the left side of FIG. 3C, the distribution of a plurality of color points in a state before the smoothing process is indicated by double circles and white circles. These color points constitute a three-dimensional color solid CS in the virtual CMY space. Note that this three-dimensional color solid CS correlates with the three-dimensional color solid in the Lab space, and each axis of the Lab space corresponding to the virtual CMY space is also shown in FIG. The virtual CMY coordinate value of each color point is a value obtained by converting the ink amount at a plurality of input grid points of the base 3D-
仮想ＣＭＹ空間の３次元色立体ＣＳは、以下の８つの頂点（図３（Ｃ）の２重丸の点）を有している。
・点ＰK：（Ｒ，Ｇ，Ｂ）＝（０，０，０）に対応する紙黒点。
・点ＰW：（Ｒ，Ｇ，Ｂ）＝（２５５，２５５，２５５）に対応する紙白点。
・点ＰC：（Ｒ，Ｇ，Ｂ）＝（０，２５５，２５５）に対応するシアン点。
・点ＰM：（Ｒ，Ｇ，Ｂ）＝（２５５，０，２５５）に対応するマゼンタ点。
・点ＰY：（Ｒ，Ｇ，Ｂ）＝（２５５，２５５，０）に対応するイエロー点。
・点ＰR：（Ｒ，Ｇ，Ｂ）＝（２５５，０，０）に対応するレッド点。
・点ＰG：（Ｒ，Ｇ，Ｂ）＝（０，２５５，０）に対応するグリーン点。
・点ＰB：（Ｒ，Ｇ，Ｂ）＝（０，０，２５５）に対応するブルー点。
The three-dimensional color solid CS in the virtual CMY space has the following eight vertices (double circle points in FIG. 3C).
Point P K : Paper black point corresponding to (R, G, B) = (0, 0, 0).
Point P W : Paper white point corresponding to (R, G, B) = (255, 255, 255).
Point P C : cyan point corresponding to (R, G, B) = (0, 255, 255).
Point P M : Magenta point corresponding to (R, G, B) = (255, 0, 255).
Point P Y : Yellow point corresponding to (R, G, B) = (255, 255, 0).
Point P R : Red point corresponding to (R, G, B) = (255, 0, 0).
Point P G : Green point corresponding to (R, G, B) = (0, 255, 0).
Point P B : Blue point corresponding to (R, G, B) = (0, 0, 255).
図３（Ｃ）の右側は、スムージング処理後の色点の分布を示している。スムージング処理は、仮想ＣＭＹ空間における複数の色点を移動させて、それらの色点の分布を等間隔に近い平滑なものにする処理である。スムージング処理では、さらに、移動後の各色点の仮想ＣＭＹ値を再現するために最適なインク量も決定される。この最適なインク量がベースＬＵＴ５１０の出力値として登録されると、ベースＬＵＴ５１０が完成する。
The right side of FIG. 3C shows the distribution of color points after the smoothing process. The smoothing process is a process of moving a plurality of color points in the virtual CMY space and smoothing the distribution of these color points close to equal intervals. In the smoothing process, an optimum ink amount is further determined in order to reproduce the virtual CMY value of each color point after movement. When this optimum ink amount is registered as the output value of the
図４（Ａ）〜（Ｃ）は、入力表色系の色点（すなわち入力格子点）と仮想ＣＭＹ空間の色点との対応関係を示している。仮想ＣＭＹ空間の３次元色立体ＣＳの頂点は、ベースＬＵＴ５１０の入力表色系の３次元色立体の頂点と一対一に対応している。また、各頂点を結ぶ辺（稜線）も、両方の色立体で互いに対応しているものと考えることができる。スムージング処理前の仮想ＣＭＹ空間の各色点は、ベースＬＵＴ５１０の入力格子点にそれぞれ対応付けられており、従って、スムージング処理後の仮想ＣＭＹ空間の各色点もベースＬＵＴ５１０の入力格子点にそれぞれ対応付けられる。なお、ベースＬＵＴ５１０の入力格子点はスムージング処理によって変化しない。スムージング処理後の仮想ＣＭＹ空間の３次元色立体ＣＳは、ベースＬＵＴ５１０の出力表色系を構成するインクセットで再現可能な色域（カラーガマット）の全体に対応している。従って、ベースＬＵＴ５１０の入力表色系は、このインクセットで再現可能な色域の全体を表す表色系としての意義を有している。
4A to 4C show the correspondence between the color points of the input color system (that is, the input grid points) and the color points of the virtual CMY space. The vertices of the three-dimensional color solid CS in the virtual CMY space have a one-to-one correspondence with the vertices of the three-dimensional color solid of the input color system of the
ベースＬＵＴ５１０を作成する際に、仮想ＣＭＹ空間においてスムージング処理を行う理由は以下の通りである。ベースＬＵＴ５１０では、なるべく大きな色域を再現できるように出力表色系のインク量を設定したいという要望がある。特定のインクセットで再現可能な色域は、インクデューティ制限（一定面積に吐出可能なインク量の制限）などの所定の制限条件を考慮して決定される。一方、上述した置換比率マトリックス３１０は、これらの制限条件が考慮されておらず、再現可能な色域とは無関係に作成されている。そこで、スムージング処理の際にインクデューティ制限等の制限条件を考慮して仮想ＣＭＹ空間内の色点の取り得る範囲を決定すれば、特定のインクセットで再現可能な色域を決定することが可能となる。なお、色点の移動を行うアルゴリズムとしては、例えば、後述する力学モデルを使用したものが利用される。
The reason for performing the smoothing process in the virtual CMY space when creating the
図２のステップＳ４００では、スムージング処理の結果を用いて、テーブル作成モジュール１４０がベースＬＵＴ５１０を作成する。すなわち、テーブル作成モジュール１４０は、各入力格子点に対応付けられた仮想ＣＭＹ空間の色点を再現するための最適なインク量をベースＬＵＴ５１０（図３（Ｃ））の出力値として登録する。なお、スムージング処理では、その計算負荷を軽減するために、ベースＬＵＴ５１０の入力格子点の一部のみに対応する色点のみを処理対象として選択することも可能である。例えば、ベースＬＵＴ５１０の入力格子点におけるＲＧＢ値の間隔が１６である場合に、スムージング処理の対象となる入力格子点におけるＲＧＢ値の間隔を３２に設定すれば、スムージング処理の負荷を半減することができる。この場合には、テーブル作成モジュール１４０は、スムージング処理結果を補間することによってベースＬＵＴ５１０のすべての入力格子点に対するインク量を決定して登録する。
In step S400 of FIG. 2, the
図５（Ａ）〜（Ｃ）は、図２のステップＳ１００〜Ｓ３００によってベース４Ｄ−ＬＵＴ５２０を作成する場合の処理内容を示す説明図である。図５（Ａ）は、図３（Ａ）と同じである。図５（Ｂ）に示すベース４Ｄ−ＬＵＴ５２０は、入力がＣＭＹＫ表色系である点が図３（Ｂ）に示したベース３Ｄ−ＬＵＴ５１０と異なっている。このベース４Ｄ−ＬＵＴ５２０の初期入力値としては、（Ｃ，Ｍ，Ｙ，Ｋ）＝（０，０，０，０），（０，０，２５５，０），（０，２５５，０，０），（０，２５５，２５５，０），（２５５，０，０，０），（２５５，０，２５５，０），（２５５，２５５，０，０），（２５５，２５５，２５５，０），（０，０，０，２５５），（０，０，２５５，２５５），（０，２５５，０，２５５），（０，２５５，２５５，２５５），（２５５，０，０，２５５），（２５５，０，２５５，２５５），（２５５，２５５，０，２５５），（２５５，２５５，２５５，２５５）である１６個の入力格子点に関してインク量の初期値が設定される。他の入力格子点に対するインク量の初期入力値は任意であり、例えば０に設定される。
FIGS. 5A to 5C are explanatory diagrams showing processing contents when the
図５（Ｃ）は、スムージング処理の様子を示している。なお、仮想ＣＭＹ空間においてベース４Ｄ−ＬＵＴ５２０に対応する色立体としては、図５（Ｃ）の右端に示すように、入力値のうちのＫ値のそれぞれの値に対して１つの３次元色立体ＣＳが存在する。この例では、Ｋ＝０の色立体とＫ＝３２の色立体とを含む複数の色立体ＣＳが図示されている。本明細書では、これらの個々の色立体ＣＳを「Ｋレイヤ」とも呼ぶ。この理由は、各色立体ＣＳが、ＣＭＹＫ値のうちのＫ値が一定でＣ，Ｍ，Ｙ値が可変である入力層に対応するものと考えることができるからである。複数の色立体ＣＳは、Ｋ値が大きいほど暗い色域を表現するものとなっている。これらの複数の色立体ＣＳは、入力表色系のＫ値が大きいほど濃ブラックインクＫのインク量が多くなるように濃ブラックインクＫのインク量を決定することによって実現できる。上述したように、再現可能な色域はインクデューティ制限値によって制限される。このインクデューティ制限値としては、個々のインクのインク量と、全インクの合計インク量と、の２つの制限値が課されるのが普通である。一方、暗い色を再現する方法としては、濃ブラックインクＫなどの無彩色インクを用いる方法と、コンポジットブラックを用いる方法とがある。しかし、コンポジットブラックは合計インク量が多くなるので、濃ブラックインクＫに比べてインクデューティ制限値に抵触する可能性が高く、暗い色を再現するのには不利である。従って、入力表色系のＫ値が大きく濃ブラックインクＫのインク量が多い色立体の方が、入力表色系のＫ値が小さく濃ブラックインクＫのインク量が少ない色立体に比べてより暗い色を再現することが可能となる。
FIG. 5C shows a state of the smoothing process. Note that as the color solid corresponding to the
図６（Ａ），（Ｂ）は、ベースＬＵＴを用いた色補正ＬＵＴの作成方法を示す説明図である。図６（Ａ）に示すように、ベース３Ｄ−ＬＵＴ５１０は、ＲＧＢ値をインク量Ｉjに変換する。このインク量Ｉjは、図３（Ｂ）に示した１０種類のインクのインク量を表している。このとき、インク量Ｉjの添え字ｊは１〜１０である。変換後のインク量Ｉjは、測色によってＬ*ａ*ｂ*値に変換される。すなわち、ベース３Ｄ−ＬＵＴ５１０もしくはベース４Ｄ−ＬＵＴ５２０の対応する印刷用紙に変換後のインク量Ｉjにてカラーパッチを印刷し、印刷されたカラーパッチを測色機等にて測色することにより、そのインク量に応じて印刷されるカラーパッチのＬａｂ表色系の色彩値を取得する。この色彩値の取得には、予め選択された光源（例えば標準の光Ｄ５０）がカラーパッチの観察条件として使用される。なお、本明細書において「カラーパッチ」という用語は、有彩色のパッチに限らず、無彩色のパッチも含む広い意味で使用される。
6A and 6B are explanatory diagrams illustrating a method for creating a color correction LUT using a base LUT. As shown in FIG. 6A, the base 3D-
一方、ｓＲＧＢ値は、既知の変換式に従ってＬ*ａ*ｂ*値に変換される。この変換後のＬ*ａ*ｂ*値は、その色域が、ベース３Ｄ−ＬＵＴによって変換されたインク量Ｉjで印刷されたカラーパッチの測色により得られるＬ*ａ*ｂ*値の色域と一致するようにガマットマッピングされる。一方、ベース３Ｄ−ＬＵＴ５１０と上述した測色により、ＲＧＢ値から変換したＬ*ａ*ｂ*値を、逆方向ルックアップテーブルとして、逆変換ＬＵＴ５１１を作成する。ガマットマッピングされたＬ*ａ*ｂ*値は、この逆変換ＬＵＴ５１１によってＲＧＢ値に変換される。このＲＧＢ値は、さらに、ベース３Ｄ−ＬＵＴ５１０によってインク量Ｉjに再度変換される。この最後のインク量Ｉjと最初のｓＲＧＢ値の対応関係をルックアップテーブルに登録することによって、色補正３Ｄ−ＬＵＴ６１０を作成することができる。この色補正３Ｄ−ＬＵＴ６１０は、ｓＲＧＢ表色系をインク表色系に変換する色変換テーブルである。
On the other hand, the sRGB values are converted into L * a * b * values according to a known conversion formula. L * a * b * value after the conversion, the color gamut, based 3D-LUT obtained by colorimetry of the printed color patches in the converted ink amount I j by L * a * b * values of Gamut mapping to match the color gamut. On the other hand, a
図６（Ｂ）は、色補正４Ｄ−ＬＵＴ６２０を作成する方法を示している。図６（Ａ）との違いは、ベース３Ｄ−ＬＵＴ５１０およびその逆変換ＬＵＴ５１１の代わりに、ベース４Ｄ−ＬＵＴ５２０およびその逆変換ＬＵＴ５２１を利用している点と、ｓＲＧＢ表色系をＬ*ａ*ｂ*値に変換する既知変換式の代わりにＪＡＰＡＮ ＣＯＬＯＲ表色系（図中では「ｊＣＭＹＫ」と記したもの）をＬ*ａ*ｂ*値に変換する既知変換式を使用している点である。良く知られているように、ＪＡＰＡＮ ＣＯＬＯＲは、ＣＭＹＫの４色で構成される表色系である。なお、図６（Ｂ）の方法では、逆変換ＬＵＴ５２１において、Ｌ*ａ*ｂ*値からＣＭＹＫ値に変換する際に、既知変換前の最初のｊＣＭＹＫ値のＫ値から、逆変換ＬＵＴ５２１のＫレイヤ（Ｋ値が一定を取る部分）が選択される。従って、色補正４Ｄ−ＬＵＴ６２０として、ベース４Ｄ−ＬＵＴ５２０のうちのＫレイヤにおける特性を反映したものを作成することが可能である。以上のようにテーブル作成モジュール１４０が実行するステップＳ４００が、本実施例において色変換テーブル作成工程を構成する。
FIG. 6B shows a method for creating the
なお、通常は、ベースＬＵＴ５１０，５２０がプリンタドライバに実装されており、色補正ＬＵＴの作成処理以外の処理にも活用されているが、ここでは他の活用例の説明は省略する。以下では、実施例の置換比率マトリックス３１０の作成手順および、置換比率マトリックス３１０に基づいたインクの発生点の制限条件の作成手順について説明した後に、スムージング処理（平滑化／最適化処理）に利用される力学モデル、スムージング処理の処理手順、および、最適化処理の内容について順次説明する。
Normally, the
（２）置換比率マトリックスの作成手順：
図７は、置換比率マトリックスの作成手順を示すフローチャートである。ステップＳ５００では、プリンターに搭載される各色のインクについて、特定のインク階調値を印刷し、その印刷結果を測色機で測色する。ここで得られる測色値は機器非依存色空間である色彩値空間を用いて取得するものとし、本実施例ではＬ*ａ*ｂ*値として取得されるものとする。この特定のインク階調値は、例えばプリンターが２５６階調で印刷可能であるときには、０〜２５５のうち任意の階調を採用する事ができる。ただし、特定の階調値としては、各色のインクについて同じ階調値を採用するものとする。
(2) Procedure for creating a substitution ratio matrix:
FIG. 7 is a flowchart showing a procedure for creating a replacement ratio matrix. In step S500, a specific ink gradation value is printed for each color ink mounted on the printer, and the print result is measured by a colorimeter. The colorimetric values obtained here are obtained using a color value space which is a device-independent color space. In this embodiment, the colorimetric values are obtained as L * a * b * values. As this specific ink gradation value, for example, when the printer can print with 256 gradations, any gradation of 0 to 255 can be adopted. However, the same gradation value is adopted for each color ink as the specific gradation value.
ステップＳ５１０では、各色の発色を代表するベクトルを作成する。発色特性ベクトルは、Ｌ*ａ*ｂ*空間において、紙白における色彩値と、ある特定のインク量階調値に対して得られた測色値と、の差分ベクトルとして表すことができる。むろん、ベクトルの基点は、紙白に限らず他のポイントであってもよい。このようにして作成された発色特性ベクトルは、各インクの発色特性を代表するベクトルであり、以下では発色特性ベクトルと呼ぶことにする。なお、一般にインク量階調値に対する色彩値は非線形になるが、本実施例では、特定のインク量階調値の発色特性により各インクの発色特性を代表させているため、この非線形性の影響を受けずに済む。よって、後述のスムージング処理において、仮想ＣＭＹ空間で格子点分布を平滑化すると、Ｌ*ａ*ｂ*空間などの色彩値空間においても階調性の優れた平滑な格子点分布を得ることができる。 In step S510, a vector representing the color development of each color is created. The color development characteristic vector can be represented as a difference vector between the color value of paper white and the colorimetric value obtained for a specific ink amount gradation value in the L * a * b * space. Of course, the base point of the vector is not limited to paper white but may be another point. The color development characteristic vector thus created is a vector representing the color development characteristics of each ink, and will be referred to as a color development characteristic vector hereinafter. In general, the color value with respect to the ink amount gradation value is non-linear. However, in this embodiment, the color development characteristic of each ink is represented by the color development characteristic of a specific ink amount gradation value. You do n’t have to. Therefore, in the smoothing process described later, if the grid point distribution is smoothed in the virtual CMY space, a smooth grid point distribution with excellent gradation can be obtained even in a color value space such as the L * a * b * space. .
ステップＳ５２０では、仮想色空間を代表する３種類のインクを選定する。ここで選定された３種類のインクは、後述のスムージング処理を行う際の作業空間の次元を構成することになる。一般的には、減法混色の３原色である濃シアン（Ｃ）、濃マゼンタ（Ｍ）、イエロー（Ｙ）の３色を選択し、これらの３色の特定のインク量階調値に対して得られた測色値である発色特性ベクトルを利用して仮想色空間の基底となる単位ベクトルを構築する。そこで、本実施例においてもこの３色を選択し、この仮想色空間を仮想ＣＭＹ空間と呼ぶことにする。また、以下の説明では、仮想ＣＭＹ空間の各成分である仮想濃シアン、仮想濃マゼンタ、仮想イエローをそれぞれＶC、ＶM、ＶYとして表すことにし、これらの各成分で指定される格子点を仮想ＣＭＹと呼ぶことにする。
むろん、３次元空間を表現するために必要な基底を構成することができる組合せであれば、プリンターに搭載されているインク種であるＣＭＹＫＬｃＬｍＬｋＬｌｋＯｒＧｒの中から任意に選択した３色の単位ベクトルを基底として構築した仮想色空間を利用することもできる。
In step S520, three types of ink representing the virtual color space are selected. The three types of inks selected here constitute the dimensions of the work space when performing the smoothing process described later. In general, three colors of dark cyan (C), dark magenta (M), and yellow (Y), which are the three subtractive primary colors, are selected, and the specific ink amount gradation value of these three colors is selected. A unit vector that is the basis of the virtual color space is constructed using the color development characteristic vector that is the obtained colorimetric value. Therefore, in this embodiment, these three colors are selected, and this virtual color space is called a virtual CMY space. In the following description, virtual dark cyan, virtual dark magenta, and virtual yellow, which are components of the virtual CMY space, are represented as V C , V M , and V Y , respectively. Will be referred to as virtual CMY.
Of course, as long as it is a combination that can form a base necessary for expressing a three-dimensional space, a unit vector of three colors arbitrarily selected from CMYKLcLmLkLlkOrGr, which are ink types mounted on the printer, is used as a base. The constructed virtual color space can also be used.
ステップＳ５３０では、各色のインク量を仮想ＣＭＹ空間に置換するための置換比率マトリックス３１０を作成する。まず、ＣＭＹの発色特性ベクトルを仮想ＣＭＹ空間の基底ベクトルに変換するための置換行列Ｍは、下記(１)式で表すことができる。
上記（１）式において、ｘｃはＣインクの発色特性ベクトルであり、ｘｍはＭインクの発色特性ベクトルであり、ｘｙはＹインクの発色特性ベクトルである。また、ベクトルの右上に付した「ｔ」は、行列・ベクトルの転置を表し、発色特性行ベクトルを転置して得られる列ベクトルであることを示す。
In step S530, a
In the above equation (1), xc is the color development characteristic vector of C ink, xm is the color development characteristic vector of M ink, and xy is the color development characteristic vector of Y ink. Further, “t” added to the upper right of the vector represents transposition of a matrix / vector, and indicates a column vector obtained by transposing the coloring characteristic row vector.
このようにして表された置換行列Ｍは、下記（２）式に示すように、仮想ＣＭＹ空間の各単位ベクトルＵＣ，ＵＭ，ＵＹを発色特性ベクトルｘｃ，ｘｍ，ｘｙに変換する行列と見做すことができる。
また、上記（２）式から、逆行列Ｍ-1は、仮想ＣＭＹ空間の基底となっているＣ、Ｍ、Ｙ各色の発色特性ベクトルを仮想ＣＭＹ空間の各基底の単位ベクトルへと規格化する行列である事を意味し、この逆行列Ｍ-1で各色のインクの発色特性ベクトルを変換して得られるベクトルは、各色のインク量を仮想ＣＭＹに置換えるための置換比率ベクトルとなる。このようにして各色について得られた置換比率ベクトルを並べると、下記（３）式に示す置換比率マトリックスＸとなる。この置換比率マトリックスＸが、上述したコンバーター３００の備える置換比率マトリックス３１０となる。
上記（３）式において、ＩC、ＩM、ＩY、ＩK、ＩLc、ＩLm、ＩLk、ＩLlk、ＩOr、ＩGrは、インク量空間の格子点における各色のインク量を示す。上記（３）式に示す置換比率マトリックスＸによれば、インク量データを仮想ＣＭＹに置換することができる。
Further, from the above equation (2), the inverse matrix M −1 normalizes the color development characteristic vector of each color of C, M, and Y, which is the basis of the virtual CMY space, to the unit vector of each basis of the virtual CMY space. This means a matrix, and a vector obtained by converting the color development characteristic vector of each color ink with the inverse matrix M −1 becomes a replacement ratio vector for replacing the ink amount of each color with virtual CMY. When the replacement ratio vectors obtained for the respective colors are arranged in this way, a replacement ratio matrix X shown in the following equation (3) is obtained. This replacement ratio matrix X becomes the
In the above equation (3), I C , I M , I Y , I K , I Lc , I Lm , I Lk , I Llk , I Or , and I Gr represent the ink amount of each color at the lattice point of the ink amount space. Show. According to the replacement ratio matrix X shown in the above equation (3), the ink amount data can be replaced with virtual CMY.
また、置換比率マトリックスＸおよびこの行列を構成する各置換比率ベクトルは、上記(３)式の下にかっこ書で示したように各ベクトル要素に負値が許容されている。本実施例にかかるプリンターにはＣＭＹの濃色インクやＬｃＬｍＬｋＬｌｋ等の単色インクのほかに、ＯｒＧｒといった特色インクを含んでいるが、これは図３，４，５に示しているＣＭＹの発色特性ベクトルを利用して構築された仮想ＣＭＹ空間では、特色インクによって表現される領域の一部が、仮想ＣＭＹの負値を含む形でしか表現できないからである。このように、仮想ＣＭＹに負値を許容することにより、特色インクが本来表す特性を表現できるようになっている。 Further, as for the replacement ratio matrix X and each replacement ratio vector constituting this matrix, negative values are allowed for each vector element as shown in parentheses below the above equation (3). The printer according to this embodiment includes a special color ink such as OrGr in addition to the CMY dark color ink and the single color ink such as LcLmLkLlk. This is the CMY color development characteristic vector shown in FIGS. This is because, in the virtual CMY space constructed using the above, a part of the area expressed by the spot color ink can be expressed only in a form including the negative value of the virtual CMY. In this way, by allowing negative values in the virtual CMY, it is possible to express the characteristics that the spot color ink originally represents.
（３）インク発生点の制御条件の作成手順：
以上のようにして作成された置換比率マトリックスＸを利用して、各インクの発生点の制御条件を作成する。図８は、インク発生点の制御条件の作成手順を示すフローチャートである。ステップＳ６００では、インク発生点の制御条件を作成する対象のインクを１つ選択する。ここでは、Ｃインクが選択された場合を例にとって以下の説明を行う。
(3) Procedure for creating ink generation point control conditions:
Using the replacement ratio matrix X created as described above, a control condition for the generation point of each ink is created. FIG. 8 is a flowchart showing a procedure for creating the ink generation point control condition. In step S600, one ink for which the ink generation point control condition is to be created is selected. Here, the following description will be given by taking the case where C ink is selected as an example.
ステップＳ６１０では、Ｃインクを除いた全てのインクの中から３色の組合せの１つを選択する。本実施例では１０色のインクを搭載するプリンターであることから、３色の組合せは９Ｃ３＝８４通りの組合せがある。本ステップ６１０では、これら８４通りの中から組合せを１つ選択し、以下のステップＳ６２０〜Ｓ６４０の処理の対象とする。なお、本実施例では、Ｌｃ，Ｍ，Ｌｍのインク組合せが選択された場合を例にとって説明を行うことにする。
In step S610, one of the three color combinations is selected from all the inks except for C ink. In this embodiment, since the printer is equipped with ink of 10 colors, there are 9 C 3 = 84 combinations of 3 colors. In this
ステップＳ６２０では、選択された３色について、下記（４）式で表される係数ベクトルαを算出する。
上記（４）式において、ｘはＣインクの置換比率ベクトルであり、Ｎは上記選択した３色の置換比率ベクトルを組み合わせた部分置換比率マトリックスであり、αは、上記ステップＳ６１０で選択された３色の置換比率ベクトルを基底としてＣインクの置換比率ベクトルであるベクトルｘを表すための係数ベクトルである。すなわち、この係数ベクトルαは、仮想ＣＭＹ空間の単位ベクトルの１つでもあるＣインクの置換比率ベクトルｘを、ステップＳ６１０で選択された３種のインクで置換しようとする時の各インクの組合せ比率（使用割合）を示すことになる。すなわち、使用割合が正であれば置換可能であることを示し、負であれば置換できないことを意味する。
In step S620, a coefficient vector α expressed by the following equation (4) is calculated for the selected three colors.
In the above equation (4), x is a C ink replacement ratio vector, N is a partial replacement ratio matrix combining the above-described three color replacement ratio vectors, and α is 3 selected in step S610. This is a coefficient vector for representing a vector x which is a C ink replacement ratio vector based on a color replacement ratio vector. That is, the coefficient vector α is a combination ratio of each ink when the C ink replacement ratio vector x, which is also one of the unit vectors in the virtual CMY space, is to be replaced with the three types of ink selected in step S610. (Use ratio). That is, if the use ratio is positive, it indicates that the replacement is possible, and if it is negative, it means that the replacement cannot be performed.
ステップＳ６３０では、係数ベクトルαに負値の要素を含むか否かを判断する。係数ベクトルαに負値の要素を含む場合は、発生制御条件の作成の元にできないためステップＳ６１０に戻って次の３色の組を選択し、係数ベクトルαに負値の要素を含まない場合はステップＳ６４０に進んで、発生制御条件を作成することになる。ただし、計算誤差や測色誤差などの要因を考慮して、これらの要因による誤差に相当する絶対値の小さい負値であれば許容して、ステップＳ６４０に進んでもよい。 In step S630, it is determined whether or not the coefficient vector α includes a negative value element. When the coefficient vector α includes a negative value element, the generation control condition cannot be created. Therefore, the process returns to step S610 to select the next set of three colors, and the coefficient vector α does not include a negative value element. In step S640, the generation control condition is created. However, in consideration of factors such as calculation error and colorimetric error, a negative value having a small absolute value corresponding to the error due to these factors may be allowed, and the process may proceed to step S640.
ステップＳ６４０では、すべての３色の組合せについて係数ベクトルαの算出を行ったか否か判断する。算出を完了していれば、ステップＳ６５０へ進み、未完了であれば他の全ての３色の組合せについての係数ベクトルの算出と負値の要素を含むか否かの判断を完了するまで、ステップＳ６１０〜Ｓ６４０の処理を繰り返し実行する。 In step S640, it is determined whether the coefficient vector α has been calculated for all three color combinations. If the calculation has been completed, the process proceeds to step S650. If the calculation has not been completed, the calculation of coefficient vectors for all other three color combinations and the determination of whether or not to include a negative element are completed. The processes of S610 to S640 are repeatedly executed.
ステップＳ６５０では、Ｃインクの発生制限条件を作成する。図９は、係数ベクトルαと、３色の組合せでＣインクを置換えた時のデューティー制限との関係を示す図である。同図に示す関係は、下記（５）式で表す事ができる。
上記（５）式において、ｉ1、ｉ2、ｉ3、は、上記３色のインクのみでデューティー制限ｄまで使えると仮定した時に使える最大のインク量を示しており、α1、α2、α3、は係数ベクトルαの各要素を示し、ｓｕｍ（α）は係数ベクトルαの要素の和を示している。以下、これら3色の最大のインク量ｉ1、ｉ2、ｉ3で表されるベクトルをｉで表すことにする。上記（５）式からも分かるように、ベクトルｉは、Ｃインクを上記３色の組合せで置換えた時に最大に置換えることのできるインク量を示している。また、上記（５）式から分かるように、ベクトルｉは係数ベクトルαの定数倍したものと成っている。
In step S650, a C ink generation restriction condition is created. FIG. 9 is a diagram showing the relationship between the coefficient vector α and the duty limit when C ink is replaced with a combination of three colors. The relationship shown in the figure can be expressed by the following equation (5).
In the above equation (5), i 1 , i 2 , i 3 indicate the maximum ink amount that can be used when it is assumed that only the above three colors of ink can be used up to the duty limit d, and α 1 , α 2 , α 3 represents each element of the coefficient vector α, and sum (α) represents the sum of the elements of the coefficient vector α. Hereinafter, the vector represented by the maximum ink amounts i 1 , i 2 , i 3 of these three colors is represented by i. As can be seen from the above equation (5), the vector i indicates the amount of ink that can be replaced to the maximum when the C ink is replaced with the combination of the three colors. Further, as can be seen from the above equation (5), the vector i is a constant multiplied by the coefficient vector α.
次に、上記３色の組合せの最大のインク量でＣインクを置換えた時の、仮想ＣＭＹ空間におけるベクトルを求める。このベクトルｖは、下記（６）式で表すことができる。
上記（６）式において、ベクトルｖは、上記３色の組合せの基底ベクトルから成る部分置換比率マトリックスＮと、上記最大のインク量のベクトルとの積で表されている。この式（６）から分かるように、ベクトルｖは、置換比率ベクトルｘと同じ向きを持っている。
Next, a vector in the virtual CMY space when the C ink is replaced with the maximum ink amount of the combination of the three colors is obtained. This vector v can be expressed by the following equation (6).
In the equation (6), the vector v is represented by the product of the partial replacement ratio matrix N composed of the basis vectors of the three colors and the vector of the maximum ink amount. As can be seen from this equation (6), the vector v has the same direction as the replacement ratio vector x.
そこで、置換比率ベクトルｘからベクトルｖへの倍率βを、下記（７）式に示すようにして求める。
このようにして求められた倍率βは、Ｃインクを３色の組合せで置換した時の最大の倍率を示すことになり、この倍率βはＣインクを、選択した３色の組合せで置換えることができる最大のインク量を表す。また、上述した８４通りの組合せでそれぞれ倍率βを求めたときに、Ｃインク単色のデューティー制限の中で最大になる倍率βが、Ｃインクを他のインクで置換えた時に最大に置換えることができるインク量を示す。以上説明した例では、Ｃインクについて最大に置換えることができるインク量の倍率βを求めたが、同様の手順で各インクについての倍率βを求めることができる。
Therefore, the magnification β from the replacement ratio vector x to the vector v is obtained as shown in the following equation (7).
The magnification β thus obtained indicates the maximum magnification when the C ink is replaced with a combination of three colors, and this magnification β can be used to replace the C ink with a selected combination of three colors. This represents the maximum amount of ink that can be printed. Further, when the magnification β is obtained for each of the 84 combinations described above, the magnification β that becomes the maximum within the duty limit of the C ink single color can be replaced to the maximum when the C ink is replaced with another ink. Indicates the amount of ink that can be produced. In the example described above, the magnification β of the ink amount that can be replaced to the maximum for C ink is obtained, but the magnification β for each ink can be obtained by the same procedure.
なお、倍率βを算出した時に、淡色インクや特色インクの場合に、単色デューティー制限内での最大の倍率βが０に成る場合があるが、これは置換える対象となるインクの組合せが存在しない事を意味する。すなわち、画一的で特定のインクセットに捉われないインク発生点の制御条件の作成を行いながらも、上述したように係数ベクトルαに負値を許容しないことから特色インクを特別扱いしたインク発生点の制御条件を作成する事ができる。すなわち、特色インクであるＯｒやＧｒに領域制限をかける必要がない。また、係数ベクトルαに負値を許容しないようにしたため、ガマットを確保すべき特色インクの使用が制限されることを回避できる。 Note that when the magnification β is calculated, in the case of light color ink or special color ink, the maximum magnification β within the single-color duty limit may be 0, but this does not include a combination of inks to be replaced. Means things. In other words, while creating uniform ink generation point control conditions that are not trapped by a specific ink set, as described above, negative values are not allowed for the coefficient vector α. Point control conditions can be created. That is, there is no need to limit the area to the special color inks Or and Gr. Further, since a negative value is not allowed for the coefficient vector α, it is possible to avoid limiting the use of special color ink that should ensure gamut.
このようにして作成されたインク量の発生制限条件を利用して、後述の最適化処理においては、下記（８）式を満たすか否かにより、選択されたインク種（ここではＣインク）を最適化処理の対象とするか否かを判断する。
上記（８）式に示すｕは後述するターゲット仮想ＣＭＹを表すベクトルであり、ｕ’はベクトルｕをＣインクの置換比率ベクトルｘの方向に射影したベクトルである。
In the optimization process described later using the ink amount generation restriction condition thus created, the selected ink type (here, C ink) is selected depending on whether or not the following expression (8) is satisfied. It is determined whether or not to be targeted for optimization processing.
In the above equation (8), u is a vector representing a target virtual CMY described later, and u ′ is a vector obtained by projecting the vector u in the direction of the C ink replacement ratio vector x.
図１０は、上記（８）式を説明する図である。同図に示すように、ベクトルｕをＣインクの置換比率ベクトルｘ方向へ射影した結果としてのベクトルをベクトルｕ’としたときに、ベクトルｘに対するベクトルｕ’の比が倍率β以下のときは、ターゲット仮想ＣＭＹを仮想ＣＭＹ空間の基底に対応するインク種（本実施例ではＣ，Ｍ，Ｙ）で表現したときに必要になるＣインクを他のインクで置換えて表現することができることを意味する。そこで、上記比が倍率β以下のときはＣインクを最適化処理の対象外とし、Ｃインクを発生させないようにする（Ｃインクの階調値を０とする）。一方、上記比が倍率βより大きければ、ターゲット仮想ＣＭＹを仮想ＣＭＹ空間の基底に対応するインク種（本実施例ではＣ，Ｍ，Ｙ）で表現したときにＣインクを使わなければ表現できないことを意味する。そこで、上記比が倍率βより大きいときは、Ｃインクを最適化処理の対象とし、Ｃインクが発生すること（Ｃインクの階調値が０以外となること）を許容する。 FIG. 10 is a diagram for explaining the equation (8). As shown in the figure, when the vector u ′ is a vector u ′ obtained by projecting the vector u in the C ink replacement ratio vector x direction, the ratio of the vector u ′ to the vector x is equal to or less than the magnification β. This means that the C ink that is required when the target virtual CMY is expressed by the ink type (C, M, Y in this embodiment) corresponding to the base of the virtual CMY space can be replaced with another ink. . Therefore, when the ratio is equal to or less than the magnification β, the C ink is excluded from the optimization process, and the C ink is not generated (the gradation value of the C ink is set to 0). On the other hand, if the ratio is larger than the magnification β, the target virtual CMY cannot be expressed unless C ink is used when the target virtual CMY is expressed by the ink type (C, M, Y in this embodiment) corresponding to the base of the virtual CMY space. Means. Therefore, when the ratio is larger than the magnification β, the C ink is targeted for optimization processing, and the generation of the C ink is allowed (the gradation value of the C ink is other than 0).
このように、最適化処理の対象とするインク種を上述したインク量の発生制限条件を利用して予め制限する事により、粒状性回避を効果的に行いつつ最適化処理の速度を向上する事ができる。むろん、限界ギリギリまで他のインクでＣインクの置換えを強要する必要がないのであれば、所定の定数ｒを用いて下記（９）式のように、βによる制限を緩和させてもよい。
（４）力学モデル：
図１１は、本実施例のスムージング処理（平滑化／最適化処理）に利用される力学モデルを示す説明図である。ここでは、仮想ＣＭＹ色空間内に複数の色点（白丸および２重丸）が配列されている様子を示している。ただし、ここでは説明の便宜上、色点の配置を２次元的に描いている。この力学モデルでは、着目色点ｇに対して次式の仮想的な力Ｆpgが係るものと仮定する。
ここで、Ｆgは着目色点ｇが隣接色点ｇｎ（ｎは１〜Ｎ）から受ける引力の合計値、Ｖgは着目色点ｇの速度ベクトル、−ｋvＶgは速度に応じた抵抗力、Ｘgは着目色点ｇの位置ベクトル、Ｘgnは隣接色点ｇｎの位置ベクトル、ｋp，ｋgは係数である。係数ｋp，ｋgは予め一定の値に設定される。なお、文中では、ベクトルを示す矢印は省略される。
(4) Mechanical model:
FIG. 11 is an explanatory diagram showing a dynamic model used for the smoothing process (smoothing / optimization process) of the present embodiment. Here, a state is shown in which a plurality of color points (white circles and double circles) are arranged in the virtual CMY color space. However, for convenience of explanation, the arrangement of the color points is drawn two-dimensionally here. In this dynamic model, it is assumed that virtual force Fp g of the following formula with respect to focus color point g is according.
Here, F g is a total value of the attractive force received by the target color point g from the adjacent color points gn (n is 1 to N), V g is a velocity vector of the target color point g, and −k v V g is a velocity. Resistance, X g is a position vector of the target color point g, X gn is a position vector of the adjacent color point gn, and k p and k g are coefficients. The coefficients k p and k g are set to constant values in advance. In the text, the arrow indicating the vector is omitted.
このモデルは、バネで互いに結ばれた質点の減衰振動モデルである。すなわち、着目色点ｇに係る仮想合力Ｆpgは、着目色点ｇと隣接色点ｇｎとの距離が大きいほど大きくなるバネ力Ｆgと、着目色点ｇの速度が大きいほど大きくなる抵抗力−ｋvＶgとの合計値である。この力学モデルでは、各色点について、位置ベクトルＸgと速度ベクトルＦgの初期値を設定した後に、微小時間経過の位置ベクトルＸgと速度ベクトルＦgを順次算出してゆく。なお、複数の色点の速度ベクトルＶgの初期値は、例えば０に設定される。このような力学モデルを利用すれば、各色点が徐々に移動して、平滑な色点分布を得ることが可能である。 This model is a damped vibration model of mass points connected to each other by springs. That is, the virtual force Fp g according to the target color point g is a focusing color point g and the spring force F g, which distance increases larger with adjacent color points gn, increases as the speed of the target color point g is large resistance -K v V g and the total value. In this dynamic model, for each color point, after setting the initial value of the position vector X g and the velocity vector F g, it slides into successively calculates the position vector X g and the velocity vector F g minute time lapse. Note that the initial value of the velocity vector V g for a plurality of color points is set to 0, for example. By using such a dynamic model, each color point moves gradually, and a smooth color point distribution can be obtained.
なお、各色点に係る力としては、バネ力Ｆgと抵抗力−ｋvＶg以外の力を用いても良い。例えば、本出願人により開示された特開２００６−１９７０８０号公報で説明されている他の種々の力をこの力学モデルで利用してもよい。また、力学モデルを適用して各色点を移動させる際に、特定の色点は、力学モデルによって移動しない拘束点として取り扱うことも可能である。 Note that forces other than the spring force F g and the resistance force −k v V g may be used as the force related to each color point. For example, various other forces described in Japanese Patent Laid-Open No. 2006-197080 disclosed by the present applicant may be used in this dynamic model. Further, when moving each color point by applying a dynamic model, a specific color point can be handled as a constraint point that does not move by the dynamic model.
（５）スムージング処理（平滑化／最適化処理）の処理手順：
図１２は、スムージング処理（図２のステップＳ３００）の典型的な処理手順を示すフローチャートである。ステップＴ１００では、スムージング処理初期値設定モジュール１２０（図１）が、スムージング処理の対象とする複数の色点を初期設定する。
(5) Smoothing processing (smoothing / optimization processing) processing procedure:
FIG. 12 is a flowchart showing a typical processing procedure of the smoothing process (step S300 in FIG. 2). In step T100, the smoothing process initial value setting module 120 (FIG. 1) initializes a plurality of color points to be smoothed.
図１３は、ステップＴ１００の詳細手順を示すフローチャートである。ステップＴ１０２では、インク量の初期入力値（図３（Ｂ），図５（Ｂ））から、スムージング処理の対象となる各色点の仮インク量が決定される。例えば、３Ｄ−ＬＵＴ用のスムージング処理では、次の（１１）式、（１２）式に従って、各入力格子点に対する仮インク量Ｉ(R,G,B)が決定される。
ここで、Ｉ(R,G,B)は、入力格子点のＲＧＢ値に対するインクセット全体のインク量（図３の例では１０種類のインクのインク量）を表している。ＲＧＢ値が０または２５５を取る入力格子点に対するインク量は、図２のステップＳ２００においてユーザによって予め入力された値である。前記（１１）式および（１２）式によれば、任意のＲＧＢ値における仮インク量Ｉ(R,G,B)を求めることが可能である。
FIG. 13 is a flowchart showing a detailed procedure of Step T100. In step T102, a temporary ink amount for each color point to be subjected to the smoothing process is determined from the initial input value of the ink amount (FIGS. 3B and 5B). For example, in the smoothing process for 3D-LUT, the temporary ink amount I (R, G, B) for each input grid point is determined according to the following equations (11) and (12).
Here, I (R, G, B) represents the ink amount of the entire ink set with respect to the RGB values of the input grid points (ink amounts of 10 types of ink in the example of FIG. 3). The ink amount for the input grid point where the RGB value is 0 or 255 is a value input in advance by the user in step S200 of FIG. According to the equations (11) and (12), the temporary ink amount I (R, G, B) at an arbitrary RGB value can be obtained.
４Ｄ−ＬＵＴ用のスムージング処理では、次の（１３）式、（１４）式に従って、各入力格子点に対する仮インク量Ｉ(C,M,Y,K)が決定される。
なお、（１３）式からも理解できるように、４Ｄ−ＬＵＴ用のインク量の初期入力値は１６個存在するので、初期入力値の設定が煩雑である。そこで、例えば、インク量の初期入力値を設定する入力格子点を、Ｋ＝０の８個の頂点、すなわち、（Ｃ，Ｍ，Ｙ，Ｋ）＝（０，０，０，０），（０，０，２５５，０），（０，２５５，０，０），（０，２５５，２５５，０），（２５５，０，０，０），（２５５，０，２５５，０），（２５５，２５５，０，０），（２５５，２５５，２５５，０）の８個の頂点と、Ｋ＝２５５の１個の頂点、例えば、（Ｃ，Ｍ，Ｙ，Ｋ）＝（０，０，０，２５５）の頂点のみとし、Ｋ＝２５５の色点のインク量を次の（１５）式または（１６）式で決定してもよい。
ここで、Ｉ(C,M,Y,0)は、Ｋ＝０の８個の頂点におけるインク量の初期入力値から、前記（１１）式と同様の式で算出されたインク量である。（１５）式の関数ｆ D1は値Ｉ(C,M,Y,0)と値Ｉ(0,0,0,255)の合計値がインクデューティ制限値をオーバーする場合に、値Ｉ(C,M,Y,0)を減じることによって、インク量Ｉ(C,M,Y,255)がインクデューティ制限値内に納まるようにする関数である。また、（１６）式の関数ｆ D2は、値Ｉ(C,M,Y,0)と値Ｉ(0,0,0,255)の合計値がインクデューティ制限値をオーバーする場合に、合計値（Ｉ(C,M,Y,0)＋Ｉ(0,0,0,255)）の全体を減じることによって、インク量Ｉ(C,M,Y,255)がインクデューティ制限値内に納まるようにする関数である。なお、インク量Ｉj（Ｉ j(R,G,B)，ΔＩj，Ｉjr，ｈjも含む）を下付文字ｊを付すことなく示す場合は、各インクのインク量Ｉjを各行要素として有する行列（ベクトル）を意味することとする。
As can be understood from the equation (13), since there are 16 initial input values of the ink amount for 4D-LUT, setting of the initial input values is complicated. Therefore, for example, the input grid points for setting the initial input value of the ink amount are eight vertices of K = 0, that is, (C, M, Y, K) = (0, 0, 0, 0), ( 0,0,255,0), (0,255,0,0), (0,255,255,0), (255,0,0,0), (255,0,255,0), ( 8 vertices of 255, 255, 0, 0), (255, 255, 255, 0) and 1 vertex of K = 255, for example, (C, M, Y, K) = (0, 0) , 0, 255) only, and the ink amount of the color point of K = 255 may be determined by the following equation (15) or (16).
Here, I (C, M, Y, 0) is an ink amount calculated from the initial input value of the ink amount at the eight vertices of K = 0 by the same expression as the expression (11). (15) of the function f D1 is the value I (C, M, Y, 0) if the total value of the values I (0, 0, 0, 255) is over the ink duty limit value I (C, M , Y, 0) is a function for keeping the ink amount I (C, M, Y, 255) within the ink duty limit value. Further, the function f D2 in the equation (16) is calculated when the total value of the value I (C, M, Y, 0) and the value I ( 0, 0, 0, 255) exceeds the ink duty limit value ( A function for keeping the ink amount I (C, M, Y, 255) within the ink duty limit value by subtracting the total of I (C, M, Y, 0) + I ( 0, 0, 0 , 255)) It is. In addition, when the ink amount I j (including I j (R, G, B) , ΔI j , I jr , h j ) is indicated without the subscript j, the ink amount I j of each ink is indicated in each row. It means a matrix (vector) having elements.
図１３のステップＴ１０４では、コンバーター３００を用いて、仮インク量に対応する仮想ＣＭＹを求める。この演算は、以下の（１７）式または（１８）式で表すことができる。
ここで、Ｖ C(R,G,B)、Ｖ M(R,G,B)、Ｖ Y(R,G,B)、Ｖ C(C,M,Y,K)、Ｖ M(C,M,Y,K)、Ｖ Y(C,M,Y,K)は変換後の仮想ＣＭＹ値を示しており、Ｘは上述した置換比率マトリックス３１０による変換を意味している。なお、これらの式からも理解できるように、この変換後の仮想ＣＭＹ値は、ベースＬＵＴの入力値であるＲＧＢ値またはＣＭＹＫ値に対応付けられている。
In Step T104 of FIG. 13, the
Here, V C (R, G, B) , V M (R, G, B) , V Y (R, G, B) , V C (C, M, Y, K) , V M (C, M, Y, K) and V Y (C, M, Y, K) indicate virtual CMY values after conversion, and X means conversion by the
図１３のステップＴ１０６では、ステップＴ１０４で得られた仮想ＣＭＹ値を、インバースモデル初期ＬＵＴ４１０（図３（Ａ））を用いてインク量に再度変換する。ここで、インバースモデル初期ＬＵＴ４１０を用いてインク量に再度変換する理由は、インク量の初期入力値や、ステップＴ１０２で決定された仮インク量が、仮想ＣＭＹ値を再現するインク量として必ずしも好ましいインク量では無いからである。一方、インバースモデル初期ＬＵＴ４１０では、画質等を考慮した好ましいインク量が登録されているので、これを用いて仮想ＣＭＹ値をインク量に再度変換すれば、その仮想ＣＭＹ値を実現するための好ましいインク量を初期値として得ることができる。但し、ステップＴ１０６を省略してもよい。
In step T106 in FIG. 13, the virtual CMY values obtained in step T104 are converted again into ink amounts using the inverse model initial LUT 410 (FIG. 3A). Here, the reason why the ink amount is converted again using the inverse model
上述のステップＴ１００の処理の結果、スムージング処理の対象となる色点について、以下の初期値が決定される。
（１）ベースＬＵＴの入力格子点の値：（Ｒ，Ｇ，Ｂ）または（Ｃ，Ｍ，Ｙ，Ｋ）
（２）各入力格子点に対応する仮想ＣＭＹ空間の色点の初期座標値：（Ｖ C(R,G,B)、Ｖ M(R,G,B)、Ｖ Y(R,G,B)）または（、Ｖ C(C,M,Y,K)、Ｖ M(C,M,Y,K)、Ｖ Y(C,M,Y,K)）
（３）各入力格子点に対応する初期インク量：Ｉ(R,G,B)またはＩ(C,M,Y,K)
以上の説明から理解できるように、初期値設定モジュール１２０は、代表的な入力格子点に関する入力初期値から他の入力格子点に関する初期値を設定する機能を有している。なお、初期値設定モジュール１２０は、スムージング処理モジュール１３０に含まれるものとしてもよい。
As a result of the process in step T100 described above, the following initial values are determined for the color points to be subjected to the smoothing process.
(1) Input grid point value of base LUT: (R, G, B) or (C, M, Y, K)
(2) Initial coordinate values of color points in the virtual CMY space corresponding to each input grid point: (V C (R, G, B) , V M (R, G, B) , V Y (R, G, B ) ) Or (, V C (C, M, Y, K) , V M (C, M, Y, K) , V Y (C, M, Y, K) )
(3) Initial ink amount corresponding to each input grid point: I (R, G, B) or I (C, M, Y, K)
As can be understood from the above description, the initial value setting module 120 has a function of setting initial values related to other input grid points from input initial values related to typical input grid points. The initial value setting module 120 may be included in the smoothing
図１２のステップＴ１２０では、色点移動モジュール１３２が、上述した力学モデルに従って仮想ＣＭＹ空間内の色点を移動させる。
In step T120 of FIG. 12, the color
図１４（Ａ）〜（Ｄ）は、図１２のステップＴ１２０〜Ｔ１５０の処理内容を示す説明図である。図１４（Ａ）に示すように、スムージング処理前には、色点の分布にはかなりの偏りがある。図１４（Ｂ）は、微小時間経過後の各色点の位置を示している。この移動後の各色点の仮想ＣＭＹ値を「ターゲット値（ＶＣＭＹt）」と呼ぶ。「ターゲット」という修飾語は、このターゲット値ＶＣＭＹtが、以下で説明するインク量の最適値の探索処理の際の目標値として使用されるからである。 14A to 14D are explanatory diagrams showing the processing contents of steps T120 to T150 in FIG. As shown in FIG. 14A, there is a considerable bias in the distribution of color points before the smoothing process. FIG. 14B shows the position of each color point after the minute time has elapsed. The virtual CMY value of each color point after the movement is referred to as “target value (VCMY t )”. The modifier "target" is because this target value VCMY t is used as a target value during the process of searching for the optimum value of the ink amount described below.
ステップＴ１３０では、インク量最適化モジュール１３４が、予め設定された目的関数Ｅを用いて、ターゲット値ＶＣＭＹtに対するインク量の最適値を探索する（図１４（Ｃ）参照）。この目的関数Ｅを用いた最適化では、力学モデルで微小量だけ移動した後の色点の座標値ＶＣＭＹtに近い仮想ＣＭＹ値を再現するインク量Ｉjの中で、複数のパラメータΔＶC，ΔＶM，ΔＶY，ΔＧＩ，ΔＣＩＩ，ΔＴＩの２乗誤差の和がなるべく小さいインク量が最適なインク量として決定される。また、最適なインク量の探索は、ステップＴ１００で設定された各入力格子点の初期インク値から開始される。従って、探索で得られるインク量は、この初期インク量を修正した値となる。後で詳述するように、（ＥＱ１）式で与えられる目的関数Ｅは、（ＥＱ２）式のようなインク量ベクトルＩに関する２次形式の関数として書き表すことができる。インク量の最適化は、このような２次形式の目的関数Ｅを用いて、２次計画法に従って実行される。このインク量最適化モジュール１３４が実行するステップＴ１３０の処理が、本実施例における最適化工程を構成する。なお、ステップＴ１３０の詳細手順や目的関数Ｅの内容については後述する。
In step T130, the ink
図１２のステップＴ１４０では、ステップＴ１３０で探索されたインク量Ｉjに対応する仮想ＣＭＹ値が、コンバーター３００で再計算される（図１４（Ｄ）参照）。ここで、仮想ＣＭＹ値を再計算する理由は、探索されたインク量Ｉjが目的関数Ｅを最小とするインク量なので、そのインク量Ｉjで再現される仮想ＣＭＹ値は、最適化処理のターゲット値ＶＣＭＹtから多少ずれているからである。こうして再計算された仮想ＣＭＹ値が、各色点の移動後の座標値として使用される。 In step T140 of FIG. 12, the virtual CMY value corresponding to the ink amount I j searched in step T130 is recalculated by the converter 300 (see FIG. 14D). Here, the reason why the virtual CMY value is recalculated is that the searched ink amount I j is the ink amount that minimizes the objective function E, so the virtual CMY value reproduced with the ink amount I j is the value of the optimization process. This is because it is slightly deviated from the target value VCMY t . The virtual CMY value recalculated in this way is used as the coordinate value after the movement of each color point.
ステップＴ１５０では、各色点の座標値の移動量の平均値（ΔＶＣＭＹ） aveが、予め設定された閾値ε以下であるか否かが判定される。移動量の平均値（ΔＶＣＭＹ） aveが閾値εよりも大きい場合には、ステップＴ１２０に戻りステップＴ１２０〜Ｔ１５０のスムージング処理が継続される。一方、移動量の平均値（ΔＶＣＭＹ） aveが閾値ε以下の場合には、色点の分布が十分に平滑になっているので、スムージング処理が終了する。なお、閾値εは、予め適切な値が実験的に決定される。 In Step T150, it is determined whether or not the average value (ΔVCMY) ave of the movement amount of the coordinate value of each color point is equal to or less than a preset threshold value ε. When the average value (ΔVCMY) ave of the movement amount is larger than the threshold value ε, the process returns to step T120 and the smoothing process of steps T120 to T150 is continued. On the other hand, when the average value (ΔVCMY) ave of the moving amount is equal to or less than the threshold value ε, the color point distribution is sufficiently smoothed, and thus the smoothing process ends. Note that an appropriate value for the threshold ε is experimentally determined in advance.
このように、本実施例の典型的なスムージング処理（平滑化／最適化処理）では、力学モデルによって各色点を微小時間間隔毎に移動させつつ、移動後の色点に対応する最適なインク量を最適化手法で探索する。そして、色点の移動量が十分に小さくなるまでそれらの処理が継続される。この結果、図３（Ｃ）に示したように、スムージング処理によって、平滑な色点分布を得ることが可能である。 As described above, in the typical smoothing process (smoothing / optimization process) of the present embodiment, the optimum ink amount corresponding to the color point after the movement is obtained while moving each color point at every minute time interval by the dynamic model. Is searched with an optimization method. Then, these processes are continued until the movement amount of the color point becomes sufficiently small. As a result, as shown in FIG. 3C, a smooth color point distribution can be obtained by the smoothing process.
（６）最適化処理の内容：
最適化処理の目的関数Ｅ（図１４（Ｃ）参照）は、インク量の関数である仮想ＣＭＹ値および画質評価指数に関するヤコビ行列Ｊを用いて表現することが可能である。ヤコビ行列Ｊは、例えば以下の（１９）式で表される。
(6) Details of optimization processing:
The objective function E (see FIG. 14C) of the optimization process can be expressed using a Jacobian matrix J relating to a virtual CMY value and an image quality evaluation index that are functions of the ink amount. The Jacobian matrix J is expressed by the following equation (19), for example.
（１９）式の右辺の第１行〜第３行は、仮想ＣＭＹ値を個々のインク量Ｉjで偏微分した値を示している。また、第４行以下は、１組のインク量Ｉj（ｊ＝１〜１０）で印刷されるカラーパッチの画質を表す画質評価指数（粒状性指数ＧＩ（Graininess Index）と、色非恒常性指数ＣＩＩ（Color Inconstancy Index）と、合計インク量ＴＩを個々のインク量Ｉjで偏微分した値を示している。なお、画質評価指数ＧＩ，ＣＩＩ，ＴＩは、その値が小さいほど、インク量Ｉjで再現されるカラーパッチの画質が良い傾向にあることを示す指数である。 The first row to the third row on the right side of the equation (19) indicate values obtained by partial differentiation of the virtual CMY values with individual ink amounts I j . In the fourth and subsequent lines, an image quality evaluation index (graininess index GI (Graininess Index) representing color image quality of a color patch printed with one set of ink amounts I j (j = 1 to 10) and color non-constancy is shown. An index CII (Color Inconstancy Index) and a value obtained by partial differentiation of the total ink amount TI with respect to each ink amount I j are shown, and the image quality evaluation indexes GI, CII, and TI indicate that the smaller the value, the more the ink amount. This is an index indicating that the image quality of the color patch reproduced with I j tends to be good.
仮想ＣＭＹ値は、コンバーター３００を用いて、以下の（２０）式でインク量Ｉjから変換される。
画質評価指数ＧＩ，ＣＩＩも、一般にインク量Ｉjの関数としてそれぞれ表現できる。
なお、（２２）式の色非恒常性指数ＣＩＩillの下付文字「ill」は、光源の種類を表している。上述した（１９）式では、光源の種類として、標準の光Ａと標準の光Ｆ１２とを用いている。なお、色非恒常性指数ＣＩＩの計算方法の例は後述するが、色非恒常性指数ＣＩＩとしては一種類または複数種類の任意の標準光源に関するものを利用することが可能である。 Note that the subscript “ill” of the color non-constancy index CII ill in the equation (22) represents the type of light source. In the above equation (19), standard light A and standard light F12 are used as the types of light sources. An example of a method for calculating the color non-constancy index CII will be described later. As the color non-constant index CII, one relating to one or a plurality of arbitrary standard light sources can be used.
粒状性指数ＧＩは、各種の粒状性予測モデルを用いて算出可能であり、例えば以下の（２４）式で算出することができる。
ここで、ａＬは明度補正係数、ＷＳ（ｕ）はカラーパッチの印刷に利用されるハーフトーンデータが示す画像のウイナースペクトラム、ＶＴＦ（ｕ）は視覚の空間周波数特性、ｕは空間周波数である。ハーフトーンデータは、カラーパッチのインク量Ｉjからハーフトーン処理（プリンター１０が実行するハーフトーン処理と同一のものとする）によって決定される。前記（２４）式は一次元で表現しているが、空間周波数の関数として二次元画像の空間周波数を算出することは容易である。粒状性指数ＧＩの計算方法としては、例えば、本出願人により開示された特表２００７−５１１１６１号公報に記載された方法や、Makoto Fujino， Ｉmage Quality Evaluation of Inkjet Prints， Japan Hardcopy '99， p.291-294に記載された方法を利用することができる。
The graininess index GI can be calculated using various graininess prediction models. For example, the graininess index GI can be calculated by the following equation (24).
Here, aL is a brightness correction coefficient, WS (u) is a winner spectrum of an image indicated by halftone data used for printing a color patch, VTF (u) is a visual spatial frequency characteristic, and u is a spatial frequency. The halftone data is determined from the ink amount I j of the color patch by halftone processing (same as the halftone processing executed by the printer 10). The equation (24) is expressed in one dimension, but it is easy to calculate the spatial frequency of the two-dimensional image as a function of the spatial frequency. As a method for calculating the graininess index GI, for example, the method described in Japanese Patent Application Publication No. 2007-511161 disclosed by the present applicant, Makoto Fujino, Image Quality Evaluation of Inkjet Prints, Japan Hardcopy '99, p. The method described in 291-294 can be used.
色非恒常性指数ＣＩＩは、例えば以下の（２５）式で与えられる。
ここで、ΔＬ*は２つの異なる観察条件下（異なる光源下）におけるカラーパッチの明度差、ΔＣ* abは彩度差、ΔＨ* abは色相差を示す。色非恒常性指数ＣＩＩの計算時には、２つの異なる観察条件下でのＬ*ａ*ｂ*値は、色順応変換（ＣＡＴ）を用いて標準観察条件（例えば標準の光Ｄ６５の観察下）に変換される。ＣＩＩについては、Billmeyer and Saltzman's Principles of Color Technology， 3rd edition， John Wiley & Sons， Ｉnc， 2000， p.129， pp. 213-215を参照。
The color non-constancy index CII is given by the following equation (25), for example.
Here, ΔL * is the lightness difference of the color patch under two different observation conditions (under different light sources), ΔC * ab is the chroma difference, and ΔH * ab is the hue difference. When calculating the color non-constant index CII, the L * a * b * values under two different viewing conditions are converted to standard viewing conditions (eg, under the observation of standard light D65) using chromatic adaptation transformation (CAT). Converted. For the CII, Billmeyer and Saltzman's Principles of Color Technology, 3rd edi t ion, John Wiley & Sons, Inc, 2000, p.129, see pp. 213-215.
ヤコビ行列Ｊの複数の成分（「要素」とも呼ぶ）のうち、例えばＶC値に関する成分は、（２６）式で与えられる。
ここで、Ｘ（Ｉｒ＋ｈｊ），Ｘ（Ｉｒ）はコンバーター３００によるインク量ＩからＶCへ変換た値、Ｉrはインク量Ｉの現在値（平滑化／最適化処理前のインク量）、ｈjはｊ番目のインク量Ｉjの微小変動量である。他の成分も同じ形式で表される。
Among a plurality of components (also referred to as “elements”) of the Jacobian matrix J, for example, a component related to the V C value is given by Expression (26).
Here, X (I r + h j ), X (I r ) are values converted from the ink amount I to V C by the
最適化の目的関数Ｅは、例えば以下の（２７）式で与えられる。
ここで、右辺の各項の最初に記載されているｗVC，ｗVM等は、各項の重みである。各項の重みｗVC、ｗVM…は、予め設定されている。
The optimization objective function E is given by the following equation (27), for example.
Here, w VC , w VM, etc. described at the beginning of each term on the right side are the weights of each term. The weights w VC , w VM ... For each term are set in advance.
（２７）式の右辺第１項ｗVC（ΔＶC−ΔＶCt）2は、仮想シアンＶCの変動量ΔＶC，ΔＶCtに関する２乗誤差である。これらの変動量ΔＶC，ΔＶCtは、次の式で与えられる。
前記（２８）式の右辺における偏微分値はヤコビ行列（（１９）式）で与えられる値であり、Ｉjは最適化処理の結果として得られるインク量、Ｉjrは現在のインク量である。第１の変動量ΔＶCは、最適化処理によるインク量の変動量ΔＩjを、ヤコビ行列の成分である偏微分値で線形変換した量である。一方、第２の変動量ΔＶCtは、ステップＴ１２０の平滑化処理で得られたターゲット値ＶCtと、現在インク量Ｉrで与えられる仮想シアンＶC（Ｉr）との差分である。なお、第２の変動量ΔＶCtは、平滑化処理の前後におけるＶC値の差分と考えることが可能である。 The partial differential value on the right side of the equation (28) is a value given by the Jacobian matrix (equation (19)), I j is the ink amount obtained as a result of the optimization process, and I jr is the current ink amount. . The first fluctuation amount ΔV C is an amount obtained by linearly converting the fluctuation amount ΔI j of the ink amount by the optimization process with a partial differential value that is a component of the Jacobian matrix. On the other hand, the second variation amount ΔV Ct is the difference between the target value V Ct obtained by the smoothing process in step T120 and the virtual cyan V C (I r ) given by the current ink amount I r . Note that the second variation amount ΔV Ct can be considered as a difference between the V C values before and after the smoothing process.
前記（２７）式の右辺の第２項以降の各項も、前記（２８）式および（２９）式と同様の式で与えられる。すなわち、目的関数Ｅは、最適化処理によるインク量の変動量ΔＩjをヤコビ行列の成分で線形変換して得られる第１の変動量ΔＶC，ΔＶM，ΔＶＹ，ΔＧＩ…と、パラメータＶC，ＶＭ，ＶＹ，ＧＩ…に関する平滑化処理の前後における第２の変動量ΔＶCt，ΔＶＭt，ΔＶYt，ΔＧＩt…と、の２乗誤差の和として与えられている。 The terms after the second term on the right side of the equation (27) are also given by the same equations as the equations (28) and (29). That is, the objective function E includes the first variation amounts ΔV C , ΔV M , ΔV Y , ΔGI... Obtained by linearly transforming the variation amount ΔI j of the ink amount by the optimization processing with the Jacobian matrix components, and the parameter V. This is given as a sum of square errors of second fluctuation amounts ΔV Ct , ΔV Mt , ΔV Yt , ΔGI t ... Before and after the smoothing process for C 1 , V M , V Y , GI.
ところで、第１の変動量ΔＶC，ΔＶM，ΔＶＹ，ΔＧＩ…は、行列を用いて以下の（３０）式および（３１）式の形式に書き表すことが可能である。
前記（２７）式は、行列を用いて（３２）式のように表記できる。
ここで、Ｔは行列の転置を表している。行列ＷMはそれぞれ対角要素に重みを配置した対角行列（（３３）式参照）であり、行列ΔＭは各パラメータに関する目標変動量ベクトル（（３４）式参照）である。
（３４）式の右辺は、各パラメータＶC，ＶM，ＶＹ，ＣＩＩ…（「要素」とも呼ぶ）に関するターゲット値と、現在のインク量Ｉrで与えられる各パラメータ値との差分である。各パラメータのターゲット値のうち、仮想ＣＭＹのターゲット値ＶＣＭＹtは平滑化処理（ステップＴ１２０）で決定される。画質評価指数のターゲット値と現在の画質評価指数から求められる目標変動量ΔＧＩt，ΔＣＩＩt，ΔＴＩtについては、いくつかの決定方法がある。第１の方法は、目標変動量ΔＧＩt，ΔＣＩＩt，ΔＴＩtとして所定の定数（例えばΔＧＩt＝−２，ΔＣＩＩt＝−１，ΔＴＩt＝１）を使用する方法である。なお、定数としてマイナスの値を使用する理由は、これらの画質評価指数は、より小さいほど高画質であることを示す指数だからである。また、粒状性指数ＧＩのターゲット値ＧＩtは、ゼロとすることも好ましい。第２の方法は、ターゲット値ＧＩt，ＣＩＩt，ＴＩtを仮想ＣＭＹ値のターゲット値ＶＣＭＹtの関数として定義しておく方法である。以上のように、各パラメータのターゲット値は最適化処理前に決められているので、目標変動量ベクトルΔＭの各成分はすべて定数である。 The right side of the equation (34) is the difference between the target value for each parameter V C , V M , V Y , CII (also referred to as “element”) and each parameter value given by the current ink amount I r. . Of Target values for each parameter, the target value VCMY t virtual CMY are determined by the smoothing processing (step T120). Image quality evaluation index of the target value and the current target variation DerutaGI t obtained from image quality evaluation index, DerutaCII t, for .DELTA.TI t, there are several determination methods. The first method, the target variation ΔGI t, ΔCII t, predetermined constant as .DELTA.TI t (e.g. ΔGI t = -2, ΔCII t = -1, ΔTI t = 1) is a method to use. The reason why negative values are used as constants is that these image quality evaluation indices are indices indicating that the smaller the image quality evaluation index, the higher the image quality. In addition, the target value GI t of the graininess index GI is preferably set to zero. The second method is a method in which the target values GI t , CII t , and TI t are defined as a function of the target value VCMY t of the virtual CMY value. As described above, since the target value of each parameter is determined before the optimization process, all components of the target variation vector ΔM are constants.
前記（３２）式の右辺の各項のうち、第３項（Ｉr TＪT+ΔＭT）ＷM（ＪＩr+ΔＭ）、は、最適化処理の結果として得られるインク量Ｉを含まないので定数である。一般に、最適化の目的関数Ｅとしては定数項は不要である。そこで、前記（３２）式から定数項を削除して全体に１／２を乗じると、次の（３５）式が得られる。 Of the terms on the right side of the equation (32), the third term (I r T J T + ΔM T ) W M (JI r + ΔM) includes the ink amount I obtained as a result of the optimization process. Since there is no constant. In general, no constant term is required as the optimization objective function E. Therefore, if the constant term is deleted from the equation (32) and the whole is multiplied by ½, the following equation (35) is obtained.
ここで、以下の（３６）式および（３７）式のように行列Ａおよびベクトルｇを定義すると、前記（３５）式は（３８）式のように書き表せる。
（３８）式で与えられる目的関数Ｅは、最適化で得られるインク量ベクトルＩに関する２次形式であることが理解できる。図１４（Ｃ）に示した（ＥＱ１）式と（ＥＱ２）式は、（２７）式と（３８）式とそれぞれ同じものである。 It can be understood that the objective function E given by the equation (38) is a quadratic form related to the ink amount vector I obtained by the optimization. The expressions (EQ1) and (EQ2) shown in FIG. 14C are the same as the expressions (27) and (38), respectively.
本実施例の最適化処理では、（３８）式のような２次形式の目的関数Ｅを用いるので、最適化手法として２次計画法を使用することが可能である。ここで、「２次計画法」とは、逐次２次計画法を含まない狭義の２次計画法を意味している。２次形式の目的関数を用いた２次計画法を利用すれば、準ニュートン法や逐次２次計画法などの他の非線形計画法を利用する場合に比べて、処理を大幅に高速化することが可能である。 In the optimization processing of the present embodiment, the quadratic form objective function E as shown in the equation (38) is used, so that the quadratic programming method can be used as the optimization method. Here, “secondary programming” means narrowly-ordered secondary programming that does not include sequential quadratic programming. Using quadratic programming with a quadratic objective function will greatly speed up processing compared to using other nonlinear programming methods such as quasi-Newton and sequential quadratic programming. Is possible.
ところで、本実施例における最適化処理によるインク量の探索は、以下の条件の下で実行される。
（最適化条件）目的関数Ｅを最小とする。
（制約条件１）インクデューティ制限を守る。
（制約条件２）インク発生点の制御条件を守る。
By the way, the search for the ink amount by the optimization process in the present embodiment is executed under the following conditions.
(Optimization condition) The objective function E is minimized.
(Restriction condition 1) The ink duty limit is observed.
(Restriction condition 2) The control condition of the ink generation point is observed.
インクデューティ制限値としては、例えば、個々のインクのインク量の最大許容値と、合計インク量の最大許容値とが使用される。例えば、各インクのインク量を８ビットで表現したとき、１０種類の個々のインクのインク量Ｉjの最大許容値を１８０に設定し、合計インク量ΣＩjの最大許容値を２４０に設定してもよい。 As the ink duty limit value, for example, the maximum allowable value of the ink amount of each ink and the maximum allowable value of the total ink amount are used. For example, when the ink amount of each ink is expressed by 8 bits, the maximum allowable value of the ink amount I j of 10 types of individual inks is set to 180, and the maximum allowable value of the total ink amount ΣI j is set to 240. May be.
インクデューティ制限は、次の（３９）式で表すことができる。
ここで、ベクトルｂは、デューティー制限の対象となるインク種類を識別するための係数であり、要素に０か１を持つベクトルである。例えば、１種類のインクに関するデューティー制限の場合には、ベクトルｂの１個の要素のみが１となる。一方、全インクの合計インク量に関するデューティー制限の場合には、ベクトルｂのすべての要素が１となる。（３９）式の右辺のlimIは、デューティー制限値である。
The ink duty limit can be expressed by the following equation (39).
Here, the vector b is a coefficient for identifying the ink type subject to duty limitation, and is a vector having 0 or 1 as an element. For example, in the case of a duty limit relating to one type of ink, only one element of the vector b is 1. On the other hand, in the case of the duty limit related to the total ink amount of all inks, all elements of the vector b are 1. Lim I on the right side of the equation (39) is a duty limit value.
各インク量Ｉjには、負でないという制約も存在する。この非負制限は、以下の（４０）式で表せる。
前記（３９）式と（４０）式とを合体すると、インクデューティ制限は、次の（４１）式で与えられる。
この（４１）式で表される制約は、線形不等号制約である。一般に、２次計画法は線形制約の下で実行することが可能である。すなわち、本実施例における最適化処理では、（４１）式の制約の下で、前記（３８）式で与えられる２次形式の目的関数Ｅを用いた２次計画法を実行することによって、最適なインク量を探索する。この結果、この線形制約を厳密に満足しつつ、インク量探索を高速に実行することが可能である。 The constraint expressed by the equation (41) is a linear inequality constraint. In general, quadratic programming can be performed under linear constraints. That is, in the optimization process in the present embodiment, the optimization is performed by executing quadratic programming using the quadratic objective function E given by the equation (38) under the constraint of the equation (41). Search for the correct ink amount. As a result, it is possible to execute the ink amount search at high speed while strictly satisfying this linear constraint.
図１５は、最適化処理（図１２のステップＴ１３０）の詳細手順を示すフローチャートである。ステップＴ１３２では、まず、前記（３４）式で与えられる目標変動量ΔＭを求める。この目標変動量ΔＭは、前述したように、ステップＴ１２０（平滑化処理）で得られたターゲット値ＶＣＭＹtと現在インク量Ｉrに基づいて決定される。 FIG. 15 is a flowchart showing a detailed procedure of the optimization process (step T130 in FIG. 12). In step T132, first, a target fluctuation amount ΔM given by the above equation (34) is obtained. As described above, this target fluctuation amount ΔM is determined based on the target value VCMY t obtained in step T120 (smoothing process) and the current ink amount I r .
ステップＴ１３４では、前記（１９）式で与えられるヤコビ行列Ｊを算出する。なお、ヤコビ行列Ｊの各成分は、前記（２６）式で例示されるように、インク量の現在値Ｉr（平滑化／最適化前の値）に関して算出される値である。 In step T134, the Jacobian matrix J given by the equation (19) is calculated. Each component of the Jacobian matrix J is a value calculated with respect to the current value I r (value before smoothing / optimization) of the ink amount, as exemplified by the equation (26).
ステップＴ１３６では、ヤコビ行列Ｊによる線形変換の結果ΔＶC，ΔＶM，ΔＶＹ，ΔＧＩ…と、目標変動量ΔＭ（ΔＶC，ΔＶM，ΔＶＹ，ΔＧＩt…）との差が最小になるように、インク量の最適化を実行する。この最適化は、前記（３８）式で与えられる２次形式の目的関数Ｅを用いた２次計画法を実行することによって実現される。 In step T136, the difference between the linear transformation results ΔV C , ΔV M , ΔV Y , ΔGI... And the target fluctuation amount ΔM (ΔV C , ΔV M , ΔV Y , ΔGI t ...) Is minimized. In this way, the ink amount is optimized. This optimization is realized by executing a quadratic programming method using the objective function E of the quadratic form given by the equation (38).
なお、図１２のフローチャートにおいて既に説明したように、ステップＴ１３０の最適化処理の後、収束が不十分と判断される場合には、平滑化処理（ステップＴ１２０）および最適化処理（ステップＴ１３０）が再度実行される。この際、平滑化／最適化処理の初期値としては、その前の平滑化／最適化処理で得られた値が利用される。なお、このような繰り返し処理は必須ではなく、少なくとも１回の平滑化／最適化処理を行えばよい。 As already described in the flowchart of FIG. 12, if it is determined that the convergence is insufficient after the optimization process in step T130, the smoothing process (step T120) and the optimization process (step T130) are performed. Will be executed again. At this time, the value obtained by the previous smoothing / optimization process is used as the initial value of the smoothing / optimization process. Note that such repeated processing is not essential, and at least one smoothing / optimization processing may be performed.
このように、本実施例では、２次形式の目的関数Ｅを用いた２次計画法を実行することによって最適なインク量を探索するので、インク量探索を高速に実行することが可能である。発明者らの実測では、従来の準ニュートン法を用いた場合に比べて約１／１０の時間で処理を完了できることが見いだされた。 As described above, in this embodiment, since the optimal ink amount is searched by executing the quadratic programming method using the objective function E of the quadratic form, the ink amount search can be executed at high speed. . The inventors have found that the processing can be completed in about one-tenth of the time when using the conventional quasi-Newton method.
（７）変形例：
なお、この発明は前記の実施例や実施形態に限られるものではなく、その要旨を逸脱しない範囲において種々の態様において実施することが可能であり、例えば次のような変形も可能である。
(7) Modification:
The present invention is not limited to the above-described examples and embodiments, and can be implemented in various modes without departing from the gist thereof. For example, the following modifications are possible.
（７−１）変形例１：
前記実施例では、平滑化処理として力学モデルを利用した処理を採用していたが、他の種類の平滑化処理を採用してもよい。例えば、隣接する色点同士の間隔を測定し、その平均値になるべく近づくように個々の間隔を調整する平滑化処理を採用することも可能である。
(7-1) Modification 1:
In the embodiment, the process using the dynamic model is adopted as the smoothing process, but other kinds of smoothing processes may be adopted. For example, it is possible to employ a smoothing process in which the distance between adjacent color points is measured and the individual distance is adjusted so as to be as close as possible to the average value.
（７−２）変形例２：
本明細書において「インク」とは、インクジェットプリンタやオフセット印刷等に用いられる液体状インクに限らず、レーザプリンタに用いられるトナーも含む広い意味で使用されている。このような「インク」の広い意味を有する他の用語としては、「色材」や「着色材」、「着色剤」を用いることも可能である。
(7-2) Modification 2:
In this specification, “ink” is not limited to liquid ink used for ink jet printers and offset printing, but is used in a broad sense including toner used for laser printers. As other terms having such a broad meaning of “ink”, “coloring material”, “coloring material”, and “coloring agent” can also be used.
（７−３）変形例３：
前記実施例では、色変換テーブルを作成する方法および装置に関して説明したが、本発明は、こうして得られた色変換テーブルを印刷装置に組み込む組み込み部を備える印刷装置製造システムにも適用可能である。色変換テーブルを作成する色変換テーブル作成装置は、この印刷装置製造システムに含まれるものとしてもよく、他のシステムや装置に含まれるものとしてもよい。なお、この製造システムの組み込み部は、例えば、プリンタドライバのインストーラ（インストールプログラム）として実現することができる。
(7-3) Modification 3:
In the above-described embodiments, the method and apparatus for creating a color conversion table has been described. However, the present invention can also be applied to a printing apparatus manufacturing system including a built-in unit that incorporates the color conversion table obtained in this way into the printing apparatus. The color conversion table creation device that creates the color conversion table may be included in the printing apparatus manufacturing system, or may be included in another system or apparatus. The built-in part of the manufacturing system can be realized as a printer driver installer (installation program), for example.
（７−４）変形例４：
上記実施例では、色変換テーブルを作成する方法および装置に関して説明したが、本発明は、こうして得られた色変換テーブルそのものや、色変換テーブルを記憶する記憶部を備えており、入力された印刷データを該色変換テーブルに基づいて変換して印刷を実行する印刷装置によっても実現することができる。
(7-4) Modification 4:
In the above embodiment, the method and apparatus for creating the color conversion table has been described. However, the present invention includes the color conversion table itself obtained in this way and a storage unit for storing the color conversion table, and the input print. It can also be realized by a printing apparatus that performs printing by converting data based on the color conversion table.
１０…プリンター、１３…ＣＰＵ、１４…ＲＡＭ、１５…ＲＯＭ、２０…メモリーカードスロット、２１…バス、２２…ＡＳＩＣ、１５ａ…プログラムデータ、１００…ベースＬＵＴ作成モジュール、１２０…スムージング処理初期値設定モジュール、１３０…スムージング処理モジュール、１３２…色点移動モジュール、１３４…インク量最適化モジュール、１３６…画質評価指数算出モジュール、１４０…テーブル作成モジュール、２００…色補正ＬＵＴ作成モジュール、３００…コンバーター、３１０…置換比率マトリックス、４００…ＬＵＴ格納部、４１０…インバースモデル初期ＬＵＴ、５１０…ベース３Ｄ−ＬＵＴ、５１１…逆変換ＬＵＴ、５２０…ベース４Ｄ−ＬＵＴ、５２１…逆変換ＬＵＴ、６１０…色補正３Ｄ−ＬＵＴ、６２０…色補正４Ｄ−ＬＵＴ
DESCRIPTION OF
Claims (7)
上記入力表色系の座標値に対応する上記インク表色系のインク量を、所定の目的関数を用いた最適化によって決定する最適化工程と、
上記複数種類のインクの中から選択した３つのインクを特定の階調値で印刷して測色した測色値の発色特性を示すベクトルを基底として構成された仮想色空間へ上記複数種類のインクのインク量をそれぞれ変換するための変換比率を示す置換比率ベクトルの中から選択された複数の置換比率ベクトルを線形結合して各単色インクを表したときの結合係数に基づいて作成された上記インク量セットの発生の制限条件を満たすインク量セットが上記最適化工程において最適化処理の対象となるように制御するインク量セット制御工程と、
上記最適化工程によって決定されたインク量に基づいて、上記入力表色系の座標値を上記インク表色系のインク量に変換するための色変換テーブルを作成する色変換テーブル作成工程と、
を備えることを特徴とする色変換テーブル作成方法。 A color conversion table creation method for converting an input color coordinate system coordinate value into an ink color set ink amount set composed of a plurality of types of ink,
An optimization step of determining the ink amount of the ink color system corresponding to the coordinate value of the input color system by optimization using a predetermined objective function;
The plurality of types of inks to a virtual color space configured based on a vector indicating color development characteristics of colorimetric values obtained by printing three inks selected from the plurality of types of inks with specific gradation values and measuring the colors. The ink created based on a combination coefficient when linearly combining a plurality of replacement ratio vectors selected from replacement ratio vectors indicating conversion ratios for converting the respective ink amounts to represent each single color ink An ink amount set control step for controlling an ink amount set that satisfies the restriction condition for the generation of an amount set to be an object of optimization processing in the optimization step;
A color conversion table creating step for creating a color conversion table for converting the coordinate value of the input color system into the ink amount of the ink color system based on the ink amount determined by the optimization step;
A color conversion table creation method characterized by comprising:
上記入力表色系の座標値に対応する上記インク表色系のインク量を、所定の目的関数を用いた最適化によって決定する最適化機能と、
上記複数種類のインクの中から選択した３つのインクを特定の階調値で印刷して測色した測色値の発色特性を示すベクトルを基底として構成された仮想色空間へ上記複数種類のインクのインク量をそれぞれ変換するための変換比率を示す置換比率ベクトルの中から選択された複数の置換比率ベクトルを線形結合して各単色インクを表したときの結合係数に基づいて作成された上記インク量セットの発生の制限条件を満たすインク量セットが上記最適化機能の実行する最適化処理の対象となるように制御するインク量セット制御機能と、
上記最適化機能によって決定されたインク量に基づいて、上記入力表色系の座標値を上記インク表色系のインク量に変換するための色変換テーブルを作成する色変換テーブル作成機能と、
を備えることを特徴とする色変換テーブル作成プログラム。 A color conversion table creation program for causing a computer to realize a function for converting an input color system coordinate value into an ink color set ink amount set composed of a plurality of types of ink,
An optimization function for determining the ink amount of the ink color system corresponding to the coordinate value of the input color system by optimization using a predetermined objective function;
The plurality of types of inks to a virtual color space configured based on a vector indicating color development characteristics of colorimetric values obtained by printing three inks selected from the plurality of types of inks with specific gradation values and measuring the colors. The ink created based on a combination coefficient when linearly combining a plurality of replacement ratio vectors selected from replacement ratio vectors indicating conversion ratios for converting the respective ink amounts to represent each single color ink An ink amount set control function for controlling an ink amount set that satisfies a restriction condition of the amount set to be an object of optimization processing executed by the optimization function;
A color conversion table creating function for creating a color conversion table for converting the coordinate value of the input color system into the ink amount of the ink color system based on the ink amount determined by the optimization function;
A color conversion table creating program comprising:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2010092027A JP2011223434A (en) | 2010-04-13 | 2010-04-13 | Color conversion table creation method, color conversion table creation program, and printer |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2010092027A JP2011223434A (en) | 2010-04-13 | 2010-04-13 | Color conversion table creation method, color conversion table creation program, and printer |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2011223434A true JP2011223434A (en) | 2011-11-04 |
Family
ID=45039781
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2010092027A Withdrawn JP2011223434A (en) | 2010-04-13 | 2010-04-13 | Color conversion table creation method, color conversion table creation program, and printer |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2011223434A (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8711432B2 (en) | 2011-06-14 | 2014-04-29 | Seiko Epson Corporation | Image processing device, printing apparatus, image processing method, and method of producing printing apparatus |
| JP2017216635A (en) * | 2016-06-01 | 2017-12-07 | 凸版印刷株式会社 | Reproduction color profile correction system, reproduction color profile correction method, and program |
| JP2018191132A (en) * | 2017-05-02 | 2018-11-29 | キヤノン株式会社 | Color conversion parameter creation device, color conversion parameter modification method, image processing device and program |
| US10375276B2 (en) | 2016-03-17 | 2019-08-06 | Seiko Epson Corporation | Color correction table creation method, color correction table creation apparatus, and non-transitory computer readable medium for storing program |
| US10540582B2 (en) | 2015-12-18 | 2020-01-21 | Seiko Epson Corporation | Method for creating color conversion table |
-
2010
- 2010-04-13 JP JP2010092027A patent/JP2011223434A/en not_active Withdrawn
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8711432B2 (en) | 2011-06-14 | 2014-04-29 | Seiko Epson Corporation | Image processing device, printing apparatus, image processing method, and method of producing printing apparatus |
| US10540582B2 (en) | 2015-12-18 | 2020-01-21 | Seiko Epson Corporation | Method for creating color conversion table |
| US10375276B2 (en) | 2016-03-17 | 2019-08-06 | Seiko Epson Corporation | Color correction table creation method, color correction table creation apparatus, and non-transitory computer readable medium for storing program |
| JP2017216635A (en) * | 2016-06-01 | 2017-12-07 | 凸版印刷株式会社 | Reproduction color profile correction system, reproduction color profile correction method, and program |
| WO2017208935A1 (en) * | 2016-06-01 | 2017-12-07 | 凸版印刷株式会社 | Reproduction color profile correction device, reproduction color profile correction method, reproduction color profile correction program, and computer readable medium including reproduction color profile correction program |
| US20190098178A1 (en) * | 2016-06-01 | 2019-03-28 | Toppan Printing Co., Ltd. | Reproduction color profile correction device, reproduction color profile correction method, reproduction color profile correction program, and computer-readable medium including reproduction color profile correction program |
| US10587779B2 (en) | 2016-06-01 | 2020-03-10 | Toppan Printing Co., Ltd. | Reproduction color profile correction device, reproduction color profile correction method, and computer-readable medium including reproduction color profile correction program |
| JP2018191132A (en) * | 2017-05-02 | 2018-11-29 | キヤノン株式会社 | Color conversion parameter creation device, color conversion parameter modification method, image processing device and program |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2011223392A (en) | Color conversion table creation method and color conversion table creation device | |
| JP5736884B2 (en) | Profile creation method, profile creation program, and printer | |
| US8335015B2 (en) | Systems and methods for building a color lookup table for a printer | |
| JP5760521B2 (en) | Profile creation method, profile creation program, and printer | |
| US8593692B2 (en) | Systems and methods for building a color lookup table for a printer | |
| JP4355964B2 (en) | Print control according to the combination of color materials | |
| JP2012049823A (en) | Gamut shape prediction method, gamut shape prediction program, and printing apparatus loaded with color conversion table created with the use of gamut predicted by the gamut shape prediction method | |
| JP2011223434A (en) | Color conversion table creation method, color conversion table creation program, and printer | |
| JP3912486B2 (en) | Color processing method, storage medium, color processing apparatus, color conversion apparatus, and image forming apparatus | |
| JP2015142250A (en) | Color conversion correspondence information creation device, method and program, and manufacturing system and manufacturing method for printer | |
| JP2010245966A (en) | Color conversion profile generation device, method, program, and printer | |
| JP2011223345A (en) | Apparatus, method, and program for creating color conversion profile, and printer | |
| JP4941152B2 (en) | Color conversion lookup table creation device and method, and program therefor | |
| JP4924458B2 (en) | Color conversion correspondence information creating apparatus and method, program therefor, and printing apparatus manufacturing system and method | |
| JP5903935B2 (en) | Profile creation method, profile creation program, and printer | |
| JP2012199832A (en) | Profile creation method, profile creation program, and printing device | |
| JP2011077869A (en) | Profile creating method, program, printer and display controller | |
| JP2010221438A (en) | Printing apparatus, printing method and method for forming lookup table | |
| JP5699776B2 (en) | Profile creation method, profile creation program, and printer | |
| JP2011139213A (en) | Print control device, print control program, and printer | |
| JP2012217075A (en) | Manufacturing method of printer, profile creating method and printing method | |
| JP2011223298A (en) | Method and program for color conversion table creation, and printer | |
| JP2012124730A (en) | Profile creation method, profile creation apparatus, and printer | |
| JP2012129912A (en) | Printer | |
| JP2012129905A (en) | Color conversion profile creation apparatus, color conversion profile creation method, color conversion profile creation program, and printer |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A300 | Withdrawal of application because of no request for examination |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20130702 |