JP3770671B2 - Color conversion method - Google Patents

Description

【００５４】
【数２】

【００５５】
なお、（５）式のように和の記号Σを使用する式を、必要に応じて、Ｌ^* ＝Σ（ｊ＝１→１０）ａ_1jｖ_j と表すこととする。
【００５６】
偏回帰係数Ａは、Ｌ^* ａ^* ｂ^* のそれぞれ別個に最小二乗法を用いて求める。例えば、Ｌ^* については、次の（６）式によるｅ_L を最小とする偏回帰係数ａ_1jを求めればよい。
【００５７】
ｅ_L ＝Σ（ｉ＝１→２４）Ｐｉ｛Ｃｉ−Σ（ｊ＝１→１０）ａ_1jｖ_j ｝² …（６）
ここで、ｉは、カラーチャート７０のパッチ番号、Ｐｉは、各色の重み、ｊは、変数の番号（１、２、…、１０）である。
【００５８】
（６）式をベクトルとマトリクスで表すと（７）式となる。ただし、（７）式において、測色値Ｃと偏回帰係数ａとはベクトル、主成分スコア［Ｖ］と重み付けマトリクス［Ｐ］とはマトリクスである。ｔはマトリクスの転置を表す。
【００５９】
ｅ_L ＝
（ベクトルＣ−ベクトルａ［Ｖ］）^t ［Ｐ］（ベクトルＣ−ベクトルａ［Ｖ］）…（７）
以下、ベクトルＣは単にＣと記載し、ベクトルａは単にａと記載する。先ず、（７）式は、以下のように変形することができる。
【００６０】
ｅ_L ＝（Ｃ^t −［Ｖ］^t ａ^t ）^t ［Ｐ］（Ｃ^t −［Ｖ］^t ａ）
ここで、ａと［Ｖ］とを入れ換えれば、
ｅ_L ＝（Ｃ−ａ［Ｖ］）［Ｐ］（Ｃ^t −［Ｖ］^t ａ）となる。
【００６１】
一般的に、（ＡＢＣ）^t ＝Ｃ^t Ｂ^t Ａ^t であるので、
ｅ_L ＝Ｃ［Ｐ］Ｃ^t ＋ａ［Ｖ］［Ｐ］［Ｖ］^t ａ^t −ａ［Ｖ］［Ｐ］Ｃ^t −Ｃ［Ｐ］［Ｖ］^t ａ^t となる。
【００６２】
ここで、［Ｖ］［Ｐ］［Ｖ］^t ＝［Ｎ］、［Ｖ］［Ｐ］Ｃ^t ＝［Ｕ］と置けば、
ｅ_L ＝Ｃ［Ｐ］Ｃ^t ＋ａ［Ｎ］ａ^t −２ａ［Ｕ］ …（８）
となる。（８）式において、ｅ_L を最小にするためには、偏回帰係数ａの各要素についての微分が０に等しくならなければならない。従って、次の（９）式が成立する。
【００６３】
【数３】

【００６４】
この（９）式から、以下に示すように偏回帰係数ａを（１０）式により求めることができる。
【００６５】

（［ ］^-1は逆行列を示す。）
このようにしてＬ^* についての偏回帰係数ａ（ａ_1j）が１０個求められ、以下、同様の方法で、残りのａ^* 、ｂ^* についても偏回帰係数ａ_2j、ａ_3jを各１０個求めることができる。求められた偏回帰係数ａ_1j、ａ_2j、ａ_3jの合計は、３×１０個である｛（４）式参照｝。
【００６６】
結局、Ｌ^* ａ^* ｂ^* についての偏回帰係数Ａをまとめて表現すると、（１１）式で表すことができる。
【００６７】
Ａ＝ＣＰＶ^t （ＶＰＶ^t ）^-1 …（１１）
次に、各目的変数Ｌ^* ａ^* ｂ^* について求めた偏回帰係数Ａの例えば５％有意で検定を行い（ステップＳ８）、５％有意（信頼度９５％）な説明変数Ｖを記憶手段に格納しておくとともに、偏回帰係数Ａを図１に示したマトリクス２０として格納する（ステップＳ９）。なお、１％有意（信頼度９９％）で検定してもよい。
【００６８】
ステップＳ８の検定では、（１２）式〜（１６）式に示すように、回帰平方和Ｓ_R ｛求めた推定値Ｌ^* ｉ（Ｌの上に記号「＾」を付けている。）と測色値の平均Ｌ^* ｉ（Ｌの上に記号「^- 」を付けている。）の差の２乗和：（１２）式参照｝と残差平方和Ｓ_E ｛測色値Ｌ^* ｉと推定値Ｌ^* ｉの差の２乗和：（１３）式参照｝とを求め、さらに回帰平方和Ｓ_R と残差平方和Ｓ_E の不偏分散Ｖ_R 、Ｖ_E ｛（１４）式、（１５）式参照｝を各々求める。さらに、偏Ｆ値Ｆｊ｛（１６）式参照｝を求める。なお、（１６）式において、ａの頭に「＾」を付けたものは、（４）式の右辺のマトリクスの係数の推定値、Ｓ^jjは、説明変数ｖ_j の分散共分散行列の逆行列の対角項を意味している。
【００６９】
【数４】

【００７０】
この偏Ｆ値Ｆｊが５％有意のＦ分布を参照して求めた値＝Ｆ′_n-p-1 （０．０５）＝Ｆ′_24-9-1（０．０５）＝Ｆ′₁₄（０．０５）＝４．６００１１より大きい場合には、５％有意水準で回帰が有意であり、予測に役立つ偏回帰係数と判断してマトリクス２０（図１参照）に格納しておく。
【００７１】
このようにして、デジタルカメラにより得られたＲＧＢ信号１２（図１参照）をＸＹＺ測色値信号１８に変換するＬ^* ａ^* ｂ^* を求めるためのマトリクス２０を作成することができる。このマトリクス２０に係る式を（１７）式に示す。
【００７２】
Ｌ^* ＝Σａ_1jｖ_j （ｊ：５％有意な変数）
ａ^* ＝Σａ_2jｖ_j （ｊ：５％有意な変数）
ｂ^* ＝Σａ_3jｖ_j （ｊ：５％有意な変数） …（１７）
この（１７）式において、（ｊ：５％有意な変数）の意味は、変数として前記検定が５％有意であったもののみを使用するという意味である。その点で上述の（５）式と異なる。なお、ａ_1j、ａ_2j、ａ_3jは、偏回帰係数マトリクスＡの要素であり、ｖ_j は、ステップＳ４で求めた説明変数の主成分スコアである。
【００７３】
このマトリクス２０を使用することにより、既成のカラーチャート７０を使用して測色変換を行うことができる。この場合、ルックアップテーブルを使用しないのでメモリ容量が少なくて済み、換言すれば、メモリ容量が少なくても、精度よく測色変換を行うことができる。さらに、重み付けマトリクスＰの使用により画像の所望色（主要色）に限定して精度よく色変換を行うこともできる。なお、説明変数の数が９箇の場合には、ステップＳ８の検定処理を行わず、全ての偏回帰係数Ａを用いて（１７）式に対応するマトリクスを作成してもよい。
【００７４】
このようにして得られたＣＩＥ−Ｌ^* ａ^* ｂ^* 色空間の値が、上述の（１）式によりＸＹＺ色空間の値へ変換され、マトリクス２０の出力信号である測色値信号１８とされる（図１参照）。
【００７５】
次に、測色値信号１８は、必要に応じて光源変更部２５により測色的に変更され、新たな測色値信号１８（同一の符号を用いる。）とされる。なお、観察光源をデジタルカメラの撮影時における撮影光源（ＣＩＥ−Ｄ５０と同一）とした場合には、この光源変更処理は不要である。
【００７６】
測色値ＸＹＺを光源変更部２５による光源変更処理後の新たな測色値ＸＹＺ（符号は、Ｘ′Ｙ′Ｚ′とする。）に変更する場合、その変換は、次の（１８）式、（１９）式、（２０）式に基づいて行うことができる。
【００７７】
Ｘ′＝ＸＸ₂ ／Ｘ₁ …（１８）
Ｙ′＝ＹＹ₂ ／Ｙ₁ …（１９）
Ｚ′＝ＺＺ₂ ／Ｚ₁ …（２０）
ただし、（１８）式〜（２０）式において、ＸＹＺは各々撮影光源下における測色値、Ｘ′Ｙ′Ｚ′は各々観察光源下における測色値、Ｘ₁ Ｙ₁ Ｚ₁ は各々撮影光源の白色点、Ｘ₂ Ｙ₂ Ｚ₂ は各々観察光源の白色点である。
【００７８】
以下、測色値信号１８は、必要に応じて光源変更処理のなされた信号であるものとする。
【００７９】
一方、スキャナ入力部から出力されたＲＧＢ信号１４は測色変換テーブル２４によりＸＹＺ信号である測色値信号２２に変換される。
【００８０】
この測色変換テーブル２４は、例えば、ｃｍｙの色パッチの各濃度を１３段階規則正しく振った、合計１３×１３×１３＝２１９７個からなるカラーパッチを有するカラーリバーサル原稿を準備する。そして、このカラーリバーサル原稿を構成する前記各カラーパッチをスキャナの入力部で読み取るとともに、測色計で読み取る。そして、スキャナで読み取ったＲＧＢ値と測色計で読み取ったＸＹＺ値との対応関係を求めてテーブルとしたものである。このテーブルに存在しない読取値間の値は補間処理により求める。
【００８１】
そして、スキャナの入力部の出力信号であるＲＧＢ信号１４を測色変換テーブル２４により測色値信号２２に変換すること、およびデジタルカメラの出力信号であるＲＧＢ信号１２を測色変換マトリクス２０により測色値信号１８に変換することにより、次に説明する測色的セットアップ部２６を共通に使用して、いわゆるオートセットアップ処理を行うことができる。すなわち、測色値信号１８、２２を、一旦、測色的ＥＮＤに変換することにより、オートセットアップ処理ソフトウエアを共通に使用することができるという利点がある。
【００８２】
次に、測色的セットアップ部２６について図３のフローチャートをも参照して説明する。
【００８３】
従来から、階調変換は、直感的な理解がし易いために濃度空間で行われており、また、現在、市販されているカラースキャナにおける階調変換処理および色修正処理等の画像処理も濃度信号を基準として実施されている。そこで、まず、測色値信号１８または測色値信号２２を濃度信号の表色系であるｃｍｙのＥＮＤ（等価中性濃度）信号５２に変換する（ステップＳ１１）。この変換処理は、図３に示すように、ＸＹＺ表色系からＲＧＢ表色系への線形変換処理（ステップＳ１１ａ）とＲＧＢ表色系からｃｍｙ表色系への非線形変換処理（ステップＳ１１ｂ）の２段階の処理になっている。なお、測色値信号１８または測色値信号２２から変換されたＥＮＤ信号を測色的ＥＮＤ信号５２と呼ぶ。
【００８４】
例えば、リバーサル原稿の測色値信号２２を測色的ＥＮＤ信号５２に変換する場合、図４に示すＣＩＥの色度図上、リバーサルフイルムの色再現領域７１（ハッチングで示す領域）を描き、その色再現領域７１を包含する領域７３の３つの原刺激ＲＧＢ（それぞれＲｘｙｚ、Ｇｘｙｚ、Ｂｘｙｚとする。）の各点の色度図上の座標、すなわち色度座標が（２１）式〜（２３）式に示す座標であるものとする。この場合、リバーサルフイルムの色再現域７１を包含する領域７３は、図４に示す色度図上、原刺激Ｒｘｙｚ、Ｇｘｙｚ、Ｂｘｙｚを頂点とする三角形の領域である。
【００８５】
Ｒｘｙｚ＝Ｒｘｙｚ（ｘ_R ，ｙ_R ，ｚ_R ） …（２１）
Ｇｘｙｚ＝Ｇｘｙｚ（ｘ_G ，ｙ_G ，ｚ_G ） …（２２）
Ｂｘｙｚ＝Ｂｘｙｚ（ｘ_B ，ｙ_B ，ｚ_B ） …（２３）
また、色度図上、ＸＹＺ表色系の基礎刺激（白色刺激）Ｗｘｙｚの座標を（２４）式に示す。
【００８６】
Ｗｘｙｚ＝Ｗｘｙｚ（ｘ_W ，ｙ_W ，ｚ_W ） …（２４）
この場合、次の（２５）式によりＸＹＺ表色系の測色値信号２２（右辺の右側のマトリクス）は、変換マトリクス（右辺の左側のマトリクス）を介してＲＧＢ表色系の色信号ＲＧＢ（左辺のマトリクス）に変換することができる（ステップＳ１１ａ）。
【００８７】
【数５】

【００８８】
余因数Ａ_ijは、下式で求められる。
【００８９】
Ａ_ij＝（−１）^i+j Ｄ_ij
ただし、Ｄ_ijは、ｉ行ｊ列を除いた小行列式である。
【００９０】
測色的ＥＮＤ信号５２のｃｍｙ値は、基礎刺激（白色刺激）Ｗ_XYZ の座標を（２５）式に代入して求めたＲ、Ｇ、Ｂの値Ｒｗ、Ｇｗ、Ｂｗに対する（２５）式で求めたＲ、Ｇ、Ｂの値の各比率を（２６）式〜（２８）式に示すように各々対数変換することにより求めることができる。
【００９１】
ｃ＝−ｌｏｇ（Ｒ／Ｒｗ） …（２６）
ｍ＝−ｌｏｇ（Ｇ／Ｇｗ） …（２７）
ｙ＝−ｌｏｇ（Ｂ／Ｂｗ） …（２８）
ラフデータ作成部５４は、セットアップ処理を短時間に行うため、例えば、原稿１枚分の測色的ＥＮＤ信号（実際には、デジタルデータ）５２を対象とするのではなく、原稿中、オペレータからＣＰＵを通じて領域指定された画像の存在する部分のみのデータを選択すること、あるいは画像が原稿の全面に存在する場合にはデータを間引いて作成すること等、いわゆるラフスキャン処理を行う。
【００９２】
次に、オートセットアップ部５８において、ラフデータ作成部５４により選択されたラフデータである測色的ＥＮＤ信号５６にもとづいてオートセットアップ処理を行う（ステップＳ１２）。このオートセットアップ処理では、例えば、特開平２−１０５６７６号公報にも公知のように、測色的ＥＮＤ信号５６についてヒストグラムを作成した後、累積ヒストグラムを作成する。
【００９３】
そして、図１中、ＥＮＤ・ＥＮＤ変換部６４の動作を示す特性図（図１中に描いている図）に示すように、前記累積ヒストグラムの、例えば、０．１％点データ（ＨＬ濃度）Ｄ１に対応する測色的ＥＮＤ信号５６の値をＤＨに設定し、９８％点データ（ＳＤ濃度）Ｄ２に対応する測色的ＥＮＤ信号５６の値をＤＳに設定する（ステップＳ１２）。
分かり易く仮想的な値で説明すると、例えば、Ｄ１の濃度がＤ１＝１．０、Ｄ２の濃度がＤ２＝２．０であった場合、Ｄ１の濃度に対応するＤＨをＤＨ＝０．１に設定し、Ｄ２の濃度に対応するＤＳをＤＳ＝３．０に設定する。なお、実際上は、０．１％点データ（ＨＬ濃度）Ｄ１が網％の０％に対応する濃度に変換され、９８％点データ（ＳＤ濃度）Ｄ２が網％の１００％に対応する濃度に変換される。
【００９４】
このようにして設定された直線７４の式から、ＥＮＤ・ＥＮＤ変換部６４により本スキャンデータ（ラフスキャン処理の対象となったデータ）の全ての測色的ＥＮＤ信号５２を測色的ＥＮＤ信号６６に変換することができる（ステップＳ１３）。すなわち、階調特性を変換すること、言い換えれば階調補正処理を行うことができる。なお、変換式は、ハイライトポイント７５（Ｄ１，ＤＨ）およびシャドーポイント７６（Ｄ２，ＤＳ）を通る曲線の式とすることもできる。なお、ＨＬ濃度Ｄ１とＳＤ濃度Ｄ２（またはセットアップポイント７５、７６）を上述のように予め定めた一定の条件に基づいて自動的に決定することをオートセットアップというが、ＨＬ濃度Ｄ１とＳＤ濃度Ｄ２の値をマニュアルで決定することも可能であり、あるいは、オートセットアップ後、ＤＨとＤＳとをマニュアルで修正することも可能である。マニュアルにより決定する動作をマニュアルセットアップという。
【００９５】
次いで、測色的ＥＮＤ逆変換マトリクス５０により、ＥＮＤ・ＥＮＤ変換後の測色的ＥＮＤ信号６６のｃｍｙ値をＸＹＺ値である測色値信号２２に逆変換する（ステップＳ１４）。
【００９６】
この逆変換処理は、測色的ＥＮＤ信号６６のｃｍｙ値をＲＧＢ表色系のＲＧＢ値に変換する処理（ステップＳ１４ａ）と、このＲＧＢ値をＸＹＺ表色系のＸＹＺ値をとる測色値信号２２に変換する処理（ステップＳ１４ｂ）とからなる。測色的ＥＮＤ信号６６のｃｍｙ値をＲＧＢ表色系のＲＧＢ値に変換する処理は、前記の（２６）式〜（２８）式をＲＧＢについて解いた次の（２９）式〜（３１）式により行われ、ＲＧＢ値をＸＹＺ表色系のＸＹＺ値に変換する処理は、前記の（２５）式をマトリクスＸＹＺについて解いた次の（３２）式により行われる。
【００９７】
Ｒ＝Ｒｗ１０^-c …（２９）
Ｇ＝Ｇｗ１０^-m …（３０）
Ｂ＝Ｂｗ１０^-y …（３１）
【００９８】
【数６】

【００９９】
なお、上述のセットアップ処理は、リバーサル原稿の測色値信号２２を例に説明したが、同様にデジタルカメラの測色値信号１８に対しても適用することができる。
【０１００】
すなわち、一般的に、印刷製版分野においては、オリジナル原稿としてリバーサル原稿を使用する。このリバーサル原稿上に原シーンが露光され、そのリバーサル原稿上で発色している色再現域７１内を対象として印刷は仕上げられる。このことは、取り扱う画像信号の色再現域として、リバーサル原稿が発色する色再現域７１で十分であることを示唆している。
【０１０１】
このように、測色的セットアップ部２６では、測色的ＸＹＺをＲＧＢ３原色の表色系へ線形変換し、基礎刺激（光源）Ｒｗ、Ｇｗ、Ｂｗに対する求めたＲ、Ｇ、Ｂの比率を対数変換することにより、ＥＮＤ値を求め、これを用いてオートセットアップする。ＥＮＤ・ＥＮＤ変更されたＥＮＤ値を逆対数変換処理することにより、変更されたＲＧＢ値を求め、逆変換マトリクスにより変更された測色値ＸＹＺを求めるようにしている。このため、既存のいわゆる濃度信号に基づいたオートセットアップ処理を行うことができるという効果が達成される。なお、セットアップ処理は、マニュアルセットアップ処理として行うことができることはいうまでもない。
【０１０２】
次に、色素濃度変換部３０について図５、図６のフローチャートをも参照して説明する。
【０１０３】
この色素濃度変換部３０では、デジタルカメラにより撮影された原シーンの測色値信号（ＸＹＺ）１８をリバーサル原稿上の色素濃度信号（ｃｍｙ）２８に変換する。
【０１０４】
この変換過程は、図５に示す方法と図６に示す方法の２通りの方法が考えられる。
【０１０５】
すなわち、図５に示す方法は、原シーンの測色値信号（ＸＹＺ）１８を、後述するルックアップテーブルにより原シーンの色素濃度信号（ｃｍｙ）８１に変換し（ステップＳ２１：第１段階の処理）、この原シーンの色素濃度信号（ｃｍｙ）８１を後述するマトリクスによりリバーサル（単に、ＲＶともいう。）原稿上の色素濃度信号（ｃｍｙ）２８に変換する方法である（ステップＳ２２：第２段階の処理）。また、図６に示す方法は、原シーンの測色値信号（ＸＹＺ）１８を、後述するマトリクスによりリバーサル原稿上の測色値信号（ＸＹＺ）８２に変換し（ステップＳ３１：第１段階の処理）、このリバーサル原稿上の測色値信号（ＸＹＺ）８２を後述するルックアップテーブルによりリバーサル原稿上の色素濃度信号（ｃｍｙ）２８に変換する（ステップＳ３２：第２段階の処理）方法である。
【０１０６】
まず、図５に示す方法の詳細について、図７のフローチャートをも参照して説明する。図５中、ステップＳ２１の処理に係る測色値ＸＹＺを色素濃度ｃｍｙに変換するためには、色空間を概ね均等に網羅したカラーチャート、例えば、上述のマクベスカラーチェッカー７０の２４色の各測色値ＸＹＺを測色計により求める（ステップＳ２１ａ：図７参照）。
【０１０７】
次に、求めた測色値ＸＹＺから色素濃度値ｃｍｙを算出し（ステップＳ２１ｂ）、測色値ＸＹＺに対する色素濃度値ｃｍｙのルックアップテーブルを作成する。このルックアップテーブルを色素濃度変換部３０における第１段階の処理であるステップＳ２１のルックアップテーブルとして色素濃度変換部３０に設定記憶しておく。なお、ステップＳ２１の処理では、このルックアップテーブルを用いるとともに、上記２４色の測色値ＸＹＺ以外の色については、その２４色の測色値ＸＹＺを用いて補間法により色素濃度値ｃｍｙを算出することができる。
【０１０８】
ステップＳ２１ｂの算出処理、すなわち測色値（三刺激値）ＸＹＺから色素濃度ｃｍｙを求める処理を説明する。
【０１０９】
測色値ＸＹＺと色素濃度ｃｍｙとの間では、次の（３３）式〜（３６）式が成立する。
【０１１０】
Ｘ＝ｋ∫ｖｉｓＰ（λ）Ｔ（λ）ｘ（λ）ｄλ …（３３）
Ｙ＝ｋ∫ｖｉｓＰ（λ）Ｔ（λ）ｙ（λ）ｄλ …（３４）
Ｚ＝ｋ∫ｖｉｓＰ（λ）Ｔ（λ）ｚ（λ）ｄλ …（３５）
Ｔ（λ）＝１０^-h …（３６）
ただし、
ｈ＝｛ｃＤｃ（λ）＋ｍＤｍ（λ）＋ｙＤｙ（λ）＋ｂａｓｅ（λ）｝
ｋ＝１００／∫ｖｉｓＰ（λ）ｙ（λ）ｄλ（λは光の波長）
∫ｖｉｓ：可視波長域（３８０ｎｍ〜７８０ｎｍ）での定積分
Ｐ（λ）：観察光源の分光特性データ
Ｔ（λ）：透過物体の色素の分光透過率データ
｛（３６）式はランベルト・ベールの法則が成立すると仮定｝
ｘ（λ）、ｙ（λ）、ｚ（λ）：等色関数
Ｄｃ（λ）、Ｄｍ（λ）、Ｄｙ（λ）：ｃｍｙ色素の分光濃度データ
ｂａｓｅ（λ）：フイルムベースの分光濃度データ
これら（３３）式〜（３６）式から色素濃度ｃｍｙを求める場合には、逆関数を求めればよいのであるが、直接求めることができない。そこで、公知のニュートン（ニュートン・ラフソン）法（例えば、太田登著「色彩工学」ｐｐ２５４−２６０ 東京電機大学出版局 １９９３年１２月２０日 第１版１刷発行）のような逐次近似法を用いればよい。前記参考書を引用して、このニュートン・ラフソン法（Ｎ・Ｒ法と略記する。）について簡単に説明する。
【０１１１】
一般的な方程式ｙ＝ｆ（ｘ）をｆ（ｘ）＝０の根に近いｘ＝ｘ０でテーラ展開して１次の項のみを取ると、ｘの微小変化Δｘに対して、ｆ（ｘ０＋Δｘ）＝ｆ（ｘ０）＋ｆ′（ｘ０）・Δｘが成立する。ただし、ｆ′（ｘ０）は、ｆ（ｘ）の微分係数ｆ′（ｘ）にｘ＝ｘ０を代入したものである。したがって、ｆ（ｘ）＝０のより正しい値ｘ１は、ｆ（ｘ０＋Δｘ）＝０として、ｘ１＝ｘ０＋Δｘ＝ｘ０−ｆ（ｘ０）／ｆ′（ｘ０）で求められる。これは、図８に示すように、関数ｙ＝ｆ（ｘ）上の点（ｘ０，ｙ０）で接線８３を引き、その接線８３とｘ軸との交点ｘ１を求めることに相当する。
【０１１２】
これを（３３）式〜（３６）式に適用するためには、（３６）式を（３３）式〜（３５）式に代入した後、ある関数ｆｘ、ｆｙ、ｆｚを用いれば、（３３）式〜（３５）式は、（３７）式〜（３９）式のように表すことができる。
【０１１３】
Ｘ＝ｆｘ（ｃ，ｍ，ｙ） …（３７）
Ｙ＝ｆｙ（ｃ，ｍ，ｙ） …（３８）
Ｚ＝ｆｚ（ｃ，ｍ，ｙ） …（３９）
これら（３７）式〜（３９）式において、初期値をｃ０，ｍ０，ｙ０とし、そのときの三刺激値をＸ０，Ｙ０，Ｚ０とする。今、ｃ０，ｍ０，ｙ０に微小変化Δｃ，Δｍ，Δｙを加えたとき、三刺激値がΔＸ，ΔＹ，ΔＺだけ変化したとすれば、次の（４０）式が得られる。
【０１１４】

ただし、例えば、∂ｆｘ／∂ｃは、関数ｆｘのｃに関する偏微分係数を表す。
【０１１５】
この（４０）式を整理すれば、（４１）式が得られる。同様に、ΔＹ、ΔＺも（４２）式、（４３）式に示すように得られる。
【０１１６】
ΔＸ＝Δｃ・∂Ｘ／∂ｃ＋Δｍ・∂Ｘ／∂ｍ＋Δｙ・∂Ｘ／∂ｙ …（４１）
ΔＹ＝Δｃ・∂Ｙ／∂ｃ＋Δｍ・∂Ｙ／∂ｍ＋Δｙ・∂Ｙ／∂ｙ …（４２）
ΔＺ＝Δｃ・∂Ｚ／∂ｃ＋Δｍ・∂Ｚ／∂ｍ＋Δｙ・∂Ｚ／∂ｙ …（４３）
（４１）式〜（４３）式を（４４）式に示すように行列表示する。
【０１１７】
（Ｑ）＝（Ｊ）（Ｐ） …（４４）
ただし、（Ｑ）は、要素が１行目から順にΔＸ，ΔＹ，ΔＺである１行３列のマトリクス、（Ｊ）は偏微分係数が３行３列からなるヤコビアンマトリクス、（Ｐ）は、要素が１行目から順にΔｃ，Δｍ，Δｙである１行３列のマトリクスである。
【０１１８】
（４４）式の両辺に、ヤコビアンマトリクス（Ｊ）の逆マトリクス（Ｊ）^-1を掛ければ、（４５）式が得られる。
【０１１９】
（Ｐ）＝（Ｊ）^-1（Ｑ） …（４５）
したがって、初期値ｃ０，ｍ０，ｙ０を（４６）式のようにそれぞれｃ１，ｍ１，ｙ１に修正すれば、より正しい近似値を得ることができる。
【０１２０】
ｃ１＝ｃ０＋Δｃ
ｍ１＝ｍ０＋Δｍ
ｙ１＝ｙ０＋Δｙ …（４６）
以上のようにして得られるヤコビアンマトリクス（Ｊ）を用いて繰り返し計算を行うことにより、任意の目標値である測色値ＸＹＺに対する色素濃度信号ｃｍｙを求めることができる。同様の処理をＸＹＺ色空間における格子上の全目標値に対して行うことで、原シーンに対する測色値信号（ＸＹＺ）１８を原シーンに対する色素濃度信号（ｃｍｙ）８１に変換するための逆変換テーブルを生成し、これをルックアップテーブルとして色素濃度変換部３０に第１段階の処理用として保持しておく（ステップＳ２１）。
【０１２１】
次に、この第１段階の処理に続いて行われるステップＳ２２の第２段階処理用のマトリクスの作成処理について説明する。なお、このマトリクスの作成処理は、図２を参照して既に説明した重回帰分析処理を含む処理であり、簡潔に説明する。
【０１２２】
まず、ステップＳ２１ｂで算出した色素濃度ｃｍｙについて、ｃ、ｍ、ｙ、ｃ² 、ｍ² 、ｙ² 、ｃｍ、ｍｙ、ｙｃの９変数２４色のデータを算出する（ステップＳ２２ａ）。
【０１２３】
次いで、この９変数２４色のデータを主成分分析し、９主成分の主成分スコアを求める（ステップＳ２２ｂ）。
【０１２４】
一方、リバーサルフイルムに露光されたマクベスカラーチェッカー２４色についての測色値を測色計により求める（ステップＳ２２ｃ）。
【０１２５】
次に、上述したＮ・Ｒ法によりリバーサルフイルム上の色素濃度ｃｍｙ（ＲＶｃｍｙとも記載する。）を求める（ステップＳ２２ｄ）。
【０１２６】
次に、リバーサルフイルム上の色素濃度ＲＶｃｍｙの各々（ＲＶｃ、ＲＶｍ、ＲＶｙ）を目的変数とし、ステップＳ２２ｂで求めた主成分スコア（定数項を含む）を説明変数として、重回帰分析処理により３×１０の偏回帰係数マトリクスを求める（ステップＳ２２ｅ）。
【０１２７】
なお、重回帰分析処理を行う際に、上述したように、目的変数群となっている２４色に１：１に対応した重み付けマトリクスを作用させてもよい（ステップＳ２２ｆ）。
【０１２８】
次に、各目的変数ＲＶｃ、ＲＶｍ、ＲＶｙについて求めた偏回帰係数について、例えば５％有意で検定を行い（ステップＳ２２ｇ）、５％有意な説明変数を記憶手段に格納しておくとともに、偏回帰係数を第２段階処理用のマトリクスとして、色素濃度変換部３０に格納しておく（ステップＳ２２ｈ）。なお、この場合においても、検定を行わず、全ての係数を用いてもよい。
【０１２９】
このようにして、原シーンの測色値信号（ＸＹＺ）１８をリバーサルフイルム上の色素濃度信号ｃｍｙに変換する色素濃度変換部３０が構築される。
【０１３０】
次に、図６に示す方法の詳細について、図９のフローチャートをも参照して説明する。なお、図９に示す処理も、図２を参照して既に説明した重回帰分析処理を含む処理であり、簡潔に説明する。
【０１３１】
まず、マクベスカラーチェッカー７０の２４色の各測色値ＸＹＺを測色計により求める（ステップＳ３１ａ：図９参照）。
【０１３２】
次に、求めた測色値ＸＹＺについて、Ｘ、Ｙ、Ｚ、Ｘ² 、Ｙ² 、Ｚ² 、ＸＹ、ＹＺ、ＺＸの９変数２４色のデータを算出する（ステップＳ３１ｂ）。
【０１３３】
次いで、この９変数２４色のデータを主成分分析し、９主成分の主成分スコアを求める（ステップＳ３１ｃ）。
【０１３４】
一方、リバーサルフイルムに露光されたマクベスカラーチェッカー２４色についての測色値ＸＹＺを測色計により求める（ステップＳ３１ｄ）。なお、リバーサルフイルムに露光されたマクベスカラーチェッカーとは、カラーリバーサルフイルムが装着された光学式カメラにより、所定の光源下で前記マクベスカラーチェッカーを撮影してリバーサルフイルムを露光した後、この露光されたリバーサルフイルムを現像して得られるリバーサル原稿をいう。
【０１３５】
次に、リバーサル原稿上の測色値ＸＹＺの各々を目的変数（ＲＶＸ、ＲＶＹ、ＲＶＺ）とし、ステップＳ３１ｃで求めた主成分スコア（定数項を含む）を説明変数として、重回帰分析処理により３×１０の偏回帰係数マトリクスを求める（ステップＳ３１ｅ）。なお、重回帰分析処理を行う際に、上述したように、目的変数群となっている２４色に１：１に対応した重み付けマトリクスを作用させてもよい（ステップＳ３１ｆ）。
【０１３６】
次に、各目的変数ＲＶＸ、ＲＶＹ、ＲＶＺについて求めた偏回帰係数の例えば５％有意で検定を行い（ステップＳ３１ｇ）、５％有意な説明変数を記憶手段に格納しておくとともに、偏回帰係数を第１段階の処理用マトリクスとして、色素濃度変換部３０に格納しておく（ステップＳ３１ｈ、ステップＳ３１）。この場合においても、検定を行わず、全ての係数を用いてもよい。
【０１３７】
次いで、ステップＳ３２の第２段階の処理用テーブルは、図７を参照して説明したステップＳ２１ｂあるいはステップＳ２２ｄの方法と同一の方法で作成されるので、その説明を省略する。
【０１３８】
このようにして、原シーンの測色値信号（ＸＹＺ）１８をリバーサル原稿上の色素濃度信号ｃｍｙに変換する図６の処理に基づく色素濃度変換部３０が構築される。
【０１３９】
次に、ｃｍｙｋ変換部３２を構成する原シーン忠実再現テーブル４２の作成処理について、図１０に示すフローチャートに基づいて説明する。
【０１４０】
原シーン忠実再現テーブル４２は、測色値信号（ＸＹＺ）１８を測色的に保存された各網％データであるｃｍｙｋ信号（ｃｍｙｋデータ）４６に変換するルックアップテーブルである。
【０１４１】
なお、この原シーン忠実再現テーブル４２の作成処理は、この出願人の出願に係る特願平８−１５４５８４号明細書に記載された作成処理と同じ処理である。
【０１４２】
この場合、まず、画像出力部３５に対して、規則的な間隔からなる複数のｃｍｙｋ網％データを与えることにより、ｃｍｙｋの濃度および混合割合が段階的に変化するカラーパッチを有するｃｍｙｋカラーチャートを作成する（ステップＳ４１）。
【０１４３】
この場合、例えば、ｃｍｙｋの規則的な間隔とは、ｃｍｙｋの各色について網％で０、２０、…、１００％等、２０％きざみで増加させた６段階の間隔とすれば、カラーパッチの総数は、４⁶ ＝１２９６になる。
【０１４４】
次いで、画像出力部３５により作成されたｃｍｙｋカラーチャートの各パッチを測色計を用いて測色し（ステップＳ４２）、この測色データから測色値（刺激値）ＸＹＺを求め、前記ｃｍｙｋ網％データから測色値ＸＹＺデータへの変換テーブル（順変換テーブルという。）を作成する（ステップＳ４３）。
【０１４５】
なお、この順変換テーブルは、補間用のテーブルとしても使用するので、前記規則的な間隔は、補間の正確性を考慮した場合には細かいほどよいが、その分、測色作業が膨大となるので、その測色作業の繁雑さと前記正確性と以下に説明するコンピュータの処理時間との、いわゆるトレードオフで決定すればよい。
【０１４６】
ところで、前記原シーン忠実再現テーブル４２において、入力される任意の測測色値信号（測色値ＸＹＺデータ、刺激値データＸＹＺ、または単にＸＹＺともいう。）１８からこれに対応するｃｍｙｋ信号（ｃｍｙｋ網％データ、色データｃｍｙｋ、または単にｃｍｙｋともいう。）４６を求めることは、変数が３変数から４変数と増加するため、１つの測色値ＸＹＺデータ１８に対してｃｍｙｋ網％データ４６の複数の解が存在する可能性がある。この問題を解決するためには、３変数同士の関係とする必要があり、ここでは、ｃｍｙｋ網％データのうち、色データｋについて画像出力部３５が採りうる最大値Ｋｍａｘ（ｋ＝Ｋｍａｘ）に固定する（ステップＳ４４）。最大値Ｋｍａｘとは、例えば、ｃｍｙｋ網％データのｋの値が１００％である。
【０１４７】
そして、３変数である任意のＸＹＺの値を３変数である対応するｃｍｙ（ｋは固定）の値に変換する（ステップＳ４５）。
【０１４８】
この場合、任意のＸＹＺの値である目標値Ｘ０、Ｙ０、Ｚ０に対するｃｍｙｋ（ここでは、ｋはｋ＝Ｋｍａｘに固定しているのでｋは定数であり、その意味でｃｍｙの３変数）の値である、ｃ０、ｍ０、ｙ０、ｋ０＝Ｋｍａｘを求める場合には、ｋ０＝Ｋｍａｘであるときの上記順変換テーブルを利用して、回帰式の偏回帰係数を求める。
【０１４９】
このときの回帰式は、係数項となる３行４列のマトリクスをＡ、測色値ＸＹＺの１行３列のマトリクスをＴ、ｋを固定した、換言すれば定数と考えて、これを含めた残りのｃｍｙの１行４列のマトリクスをＤと表せば、次の（４７）式のように表すことができる。
【０１５０】
Ｔ＝ＡＤ …（４７）
この（４７）式は、次の（４８）式の関係を表している。
【０１５１】
【数７】

【０１５２】
（４７）式、（４８）式において、マトリクスＤを構成する「１」は、ｃｍｙの３次元平面式における定数項を与えるために設定した値である。
【０１５３】
この（４８）式における係数Ａは、ｋ＝Ｋｍａｘのときに得られている上述の順変換テーブルの各データセットを代入することで、上述した重回帰分析により求めることができる。
【０１５４】
次に、ｋ＝Ｋｍａｘのときに求められた回帰式を用いて、上述のＮ・Ｒ法により、目標値Ｘ０、Ｙ０、Ｚ０に対応するｃｍｙ（ｋはＫｍａｘとして設定されている。）の値を求めることができる。
【０１５５】
次に、求められたｃｍｙの値が画像出力部３５の再現範囲内の値（色）であるかどうかの判断を行う（ステップＳ４６）。すなわち、網％データ（色データ）ｃｍｙの再現可能な最小濃度をＣｍｉｎ、Ｍｍｉｎ、Ｙｍｉｎ、最大濃度をＣｍａｘ、Ｍｍａｘ、Ｙｍａｘとして求められたデータｃｍｙが次の（４９）式の関係を満足するかどうかを判定する。
【０１５６】
Ｃｍｉｎ≦ｃ≦Ｃｍａｘ
Ｍｍｉｎ≦ｍ≦Ｍｍａｘ
Ｙｍｉｎ≦ｙ≦Ｙｍａｘ …（４９）
求めたｃｍｙの各値が、（４９）式を満足している場合には、刺激値データＸＹＺの目標値（Ｘ０，Ｙ０，Ｚ０）に対する色データｃｍｙｋを各々ｃ＝ｃｓｏｌ、ｍ＝ｍｓｏｌ、ｙ＝ｙｓｏｌ、ｋ＝ｋｓｏｌ（この場合、ｋｓｏｌ＝Ｋｍａｘ）として逆変換テーブルを作成し、これを原シーン忠実再現テーブル４２として設定する（ステップＳ４７）。以下の説明において、必要に応じて（ｃｓｏｌ，ｍｓｏｌ，ｙｓｏｌ，ｋｓｏｌ）のデータセットをｃｍｙｋｓｏｌという。
【０１５７】
ステップＳ４５の過程で求めたｃｍｙデータの値が（４９）式を満足しない値であった場合には、ｋ＝Ｋｍａｘに固定してある色データｋをｋ＝ｋ−Δｋ、この場合ｋ＝Ｋｍａｘ−Δｋとした後（ステップＳ４８）、色データｋが所定の最小値ｋ＝Ｋｍｉｎよりも小さくならない範囲で（ステップＳ４９）、ステップＳ４５の処理を繰り返す。なお、微小変化分Δｋは、第１逆変換テーブルを構成する色データｋの任意のデータ間隔であり、例えば、色データｋが０％〜１００％の範囲のデータとして設定される場合には、１％きざみ等に設定すればよく、色データｋが０〜２５５の範囲のデータとして設定されるのであれば、値１きざみに設定すればよい。ステップＳ４５の２回目の処理を行う場合には、色データｋ＝Ｋｍａｘ−１＝１００−１＝９９に対する、（４８）式の左辺の測色値であるＸＹＺ値は、ｋ＝Ｋｍａｘ＝１００％とｋ＝８０％の値で測色され順変換テーブルとして格納されている各ＸＹＺデータの補間処理（内挿処理）により求めることができる。
【０１５８】
一方、ステップＳ４９の処理において、Ｋｍｉｎ＞ｋと判定された場合には、目標値（Ｘ０，Ｙ０，Ｚ０）に対する色データｃｍｙｋを画像出力部３５の色再現域外のデータに指定し、ここでは、色データｃｍｙｋｓｏｌの算出は行わない（ステップＳ５０）。
【０１５９】
以上の処理を原シーン忠実再現テーブル４２に供給される全ての刺激値データＸＹＺ１８を目標値（Ｘ０，Ｙ０，Ｚ０）として行うことにより、画像出力部３５の色再現域内の色データである網％データｃｍｙｋ４６を得ることのできる刺激値データＸＹＺ１８に対して、最大となる色データｋを与えた場合の色データｃｍｙｋｓｏｌを求めることができる（ステップＳ５１）。
【０１６０】
なお、ステップＳ４４における色データｋの固定は、最大値Ｋｍａｘではなく、その最小値Ｋｍｉｎ（ｋ＝Ｋｍｉｎ）に固定してもよく、その場合には、ステップＳ４８の処理をｋ＝ｋ＋Δｋとし、ステップＳ４９の処理をＫｍａｘ＜ｋとすればよい。また、色データｋを任意の値とすることも容易に考えられる。例えば、ステップＳ４８の処理は、任意の値に設定したｋについて、ｋ＝ｋ−Δｋとｋ＝ｋ＋Δｋを交互に行えばよく、ステップＳ４９の処理は、ステップＳ４８の処理に対応させて行えばよい。
【０１６１】
さらに、ステップＳ４９、Ｓ５０において色再現域外データに指定された色データｃｍｙｋは、この発明の要旨ではないので詳しくは説明しないが、前記特願平８−１５４５８４号に記載されているように、いわゆるガマットマッピング（Gamat Mapping ）処理による色データＣＭＹＫの圧縮処理やクリッピング処理により、目標値（Ｘ０，Ｙ０，Ｚ０）に対する色データｃｍｙｋの逆変換テーブルを作成することができる。
【０１６２】
このようにして、測色値信号１８を構成する全ての目標値（Ｘ０，Ｙ０，Ｚ０）に対して求められたｃｍｙｋ信号４６の対応表を原シーン忠実再現テーブル４２として格納しておくことで、任意の測色値信号１８を画像出力部３５の色再現範囲内のｃｍｙｋ信号４６に変換することができる。
【０１６３】
次に、ｃｍｙｋ変換部３２を構成する標準条件再現テーブル４４の作成について説明する。
【０１６４】
この標準条件再現テーブル４４は、図１に示すように、デジタルカメラにより取り込まれたＲＧＢ信号１２がマトリクス２０により測色値信号（ＸＹＺ）１８に変換され、この測色値信号ＸＹＺが色素濃度変換部３０により色素濃度信号（ｃｍｙ）２８に変換され、この変換後の色素濃度信号（ｃｍｙ）２８を網％信号であるｃｍｙｋ信号４８に変換するためのテーブルである。
【０１６５】
また、標準条件再現テーブル４４は、スキャナにより取り込まれたＲＧＢ信号１４がテーブル２４により測色値信号（ＸＹＺ）２２に変換され、この変換後の測色値信号（ＸＹＺ）２２が色素濃度変換部３６により色素濃度信号（ｃｍｙ）３４に変換され、この変換後の色素濃度信号（ｃｍｙ）３４を網％信号であるｃｍｙｋ信号に変換するためのテーブルである。この場合、色素濃度変換部３６は、前述の色素濃度変換部３０の作成過程で示した図７中、ステップＳ２１ｂ、Ｓ２２ｄの方法と同様の過程で作成することができるので、その作成過程についての説明は省略する。
【０１６６】
標準条件再現テーブル４４の作成処理を、図１１に示すフローチャートを参照して説明する。
【０１６７】
まず、リバーサル原稿上に色素濃度ｃｍｙが格子状に配列された１３×１３×１３個のカラーパッチを有するカラーチャートを作成する（ステップＳ６１）。このカラーチャートは、ｃ（シアン）ｍ（マゼンタ）ｙ（黄）が、それぞれ最小濃度から最大濃度まで１３段階振られた各カラーパッチを有するものである。
【０１６８】
次に、このカラーチャートをスキャナのデフォルト分解条件で分解する、言い換えれば、透過原稿であるカラーチャートをスキャナで読み取り、デジタルデータ化する（ステップＳ６２）。スキャナのデフォルト分解条件は、少なくとも、階調特性変換処理、色修正処理、ＵＣＲ（下色除去）処理およびＫ版生成処理を含んでいる。ただし、分解するとき、カラーパッチのグレーの最小濃度を、ｃｍｙｋ網％が全て０％となるようにし、かつカラーパッチのグレーの最大濃度を、いわゆるベタとなるようなｃｍｙｋ網％となるように設定する。
【０１６９】
そして、各カラーパッチの色素濃度とスキャナのデフォルト分解条件で分解したときのｃｍｙｋ網％値を対応させて変換テーブル（ルックアップテーブル）を作成し（ステップＳ６３）、この変換テーブルを標準条件再現テーブル４４として設定する。なお、実際には、この標準条件再現テーブル４４と補間処理により、最小濃度から最大濃度までの任意の色素濃度信号（ｃｍｙ）２８または色素濃度信号（ｃｍｙ）３４を所望のｃｍｙｋ網％信号４８に変換することができる。
【０１７０】
以上詳述したように、上述の実施の形態によれば、デジタルカメラにより撮影したシーン（原シーン）のＲＧＢ信号１２を測色変換部１６を構成するマトリクス２０により測色値信号（ＸＹＺ）１８に変換し、この測色値信号（ＸＹＺ）１８に対して必要に応じて測色的セットアップ部２６においてセットアップ処理を行い、測色値信号１８を色素濃度変換部３０によりリバーサル原稿上の色素濃度信号（ｃｍｙ）２８に変換している。
【０１７１】
また、スキャナ入力部でリバーサル原稿から得られたＲＧＢ信号１４を、測色変換部１６を構成するテーブル２４により測色値信号（ＸＹＺ）２２に変換し、この測色値信号（ＸＹＺ）２２に対して必要に応じて測色的セットアップ部２６においてセットアップ処理を行い、測色値信号２２を色素濃度変換部３６により色素濃度信号（ｃｍｙ）３４に変換している。
【０１７２】
この場合、ｃｍｙｋ変換部３２を構成するスイッチ４５の一方の入力ポートに色素濃度変換部３６を通じて供給される色素濃度信号（ｃｍｙ）３４は、リバーサル原稿上の色素濃度信号であり、スイッチ４５の他方の入力ポートに色素濃度変換部３０を通じて供給される色素濃度信号（ｃｍｙ）２８も上述したようにリバーサル原稿上の色素濃度信号である。
【０１７３】
標準条件再現テーブル４４には、このリバーサル原稿上の色素濃度信号を分解するための、階調補正処理条件、色修正条件、下色除去条件、Ｋ版生成処理条件等の標準のデフォルト分解条件が記憶されており、このデフォルト分解条件によりリバーサル原稿上の色素濃度信号２８、３４を各々分解して網％信号であるｃｍｙｋ信号に変換する。
【０１７４】
したがって、この実施の形態によれば、デジタルカメラの出力（ＲＧＢ信号１２）またはスキャナの入力部からの出力（ＲＧＢ信号１４）に応じて、画像信号処理装置１０を構成するスイッチ４５を対応する入力ポート側に切り換えることで、いずれの出力も、結果として、共通の標準のデフォルト分解条件を使用して網％信号であるｃｍｙｋ信号に簡易に変換することができるという効果が達成される。
【０１７５】
なお、この発明は上述の実施の形態に限らず、この発明の要旨を逸脱することなく種々の構成を採り得ることはもちろんである。
【０１７６】
【発明の効果】
以上説明したように、この発明によれば、デジタルカメラの出力ＲＧＢ信号を測色値信号に変換した後、この測色値信号を色素濃度信号に変換する。また、画像読取装置で読み取った出力ＲＧＢ信号を測色値信号に変換した後、この測色値信号を色素濃度信号に変換する。色素濃度信号を階調変換処理、色修正処理、下色除去処理およびＫ版生成処理を含む標準の画像処理条件で処理して網％信号であるｃｍｙｋ信号に変換する。
【０１７７】
このため、入力されるＲＧＢ信号がデジタルカメラのものであっても、画像読取装置のものであっても、経験的に最適化された特定の画像読取装置、例えば、スキャナに標準で搭載されている（記憶されている）、いわゆるデフォルトの分解条件をそのまま利用して所望の画像を出力することができるという効果が達成される。また、ＲＧＢ信号をデバイス非依存の測色値信号に変換しているので、デバイスに依存しない画像を出力することができる。
【図面の簡単な説明】
【図１】この発明の一実施の形態が適用された画像信号処理装置の全体的な構成を示すブロック図である。
【図２】測色変換マトリクスの作成処理の説明に供されるフローチャートである。
【図３】測色的セットアップ部の説明に供されるフローチャートである。
【図４】ＣＩＥ色度図上でのリバーサルフイルムの色再現域等の説明に供される色度図である。
【図５】デジタルカメラの測色信号をリバーサル原稿上の色素濃度信号に変換する変換過程の説明に供されるフローチャートである。
【図６】デジタルカメラの測色信号をリバーサル原稿上の色素濃度信号に変換する他の変換過程の説明に供されるフローチャートである。
【図７】色素濃度変換マトリクスの作成処理の説明に供されるフローチャートである。
【図８】ニュートン・ラフソン法の説明に供される線図である。
【図９】色素濃度変換マトリクスの作成処理の一部の説明に供されるフローチャートである。
【図１０】原シーン忠実再現テーブルの作成処理の説明に供されるフローチャートである。
【図１１】標準条件再現テーブルの作成処理の説明に供されるフローチャートである。
【符号の説明】
１０…画像信号処理装置 １２、１４…ＲＧＢ信号
１６…測色変換部 １８、２２…測色値信号（ＸＹＺ）
２０…測色変換マトリクス ２４…測色変換テーブル
２５…光源変更部 ２６…測色的セットアップ部
２８、３４…色素濃度信号（ｃｍｙ） ３０、３６…色素濃度変換部
３２…ｃｍｙｋ変換部 ３５…画像出力部
４２…原シーン忠実再現テーブル ４４…標準条件再現テーブル
４６、４８…ｃｍｙｋ信号 ６４…ＥＮＤ・ＥＮＤ変換部[0001]
BACKGROUND OF THE INVENTION
The present invention relates to RGB which are digital signals respectively output from an input unit of a color scanner (hereinafter also referred to simply as a scanner) and / or a color digital camera (hereinafter also simply referred to as a digital camera) as an image reading apparatus. (RGB stands for red, green, and blue, respectively). Signals for image output devices such as film creation devices, printing plate creation devices, and color printers (cmyk are meanings for cyan, magenta, yellow, and black, respectively. ) A color conversion method for converting into a signal.
[0002]
[Prior art]
The color scanner has a color reversal document carrying a positive image as a document in its input section (a document carrying a color reversal film exposed and developed, and hereinafter also simply referred to as a reversal document). The transmitted light including the positive image information is obtained by irradiating the image with the illumination light, and this transmitted light is passed through the imaging optical system and the three-color separation optical system to a linear image sensor such as a CCD linear image sensor. The photoelectric conversion signal that is an output signal of the linear image sensor is AD-converted to photoelectrically read the image on the color reversal document.
[0003]
In this case, for example, the linear image sensor for RGB reads the reversal document by performing line scanning in the main scanning direction (line scanning) and relatively reads the reversal document in the sub-scanning direction substantially orthogonal to the main scanning direction. By carrying it, a digital RGB image signal (also simply referred to as an RGB signal), which is two-dimensional image information for each RGB color, can be obtained.
[0004]
The RGB image signal obtained in this way is converted into a cmy signal as a dye density signal, and then a cmyk image signal (also simply referred to as a cmyk signal), each of which is a halftone signal, by a cmyk conversion unit inside the scanner. .).
[0005]
The cmyk signal is supplied to, for example, an external halftone dot processing apparatus. In this halftone dot processing apparatus, each of the cmyk signals is compared with the threshold value in the threshold value template, so-called halftone processing (binarization). Process).
[0006]
Then, based on the binary image signal after the halftoning process, a plate making film is created in the film making device, and a plate is made from the plate making film by the plate making device. The produced printing plate (each plate of cmyk) is mounted on a color printing machine, inked, and multicolor printing is performed on the paper. As a result, a printed matter on which a printed image is formed with the ink transferred onto the paper is completed.
[0007]
Recently, a CTP (Computer To Plate) apparatus that directly creates a printing plate based on a binary image signal after halftone processing has been put into practical use.
[0008]
[Problems to be solved by the invention]
By the way, in the above-described cmyk conversion unit, conversion corresponding to the color printing machine is performed, and default separation conditions for this conversion are stored as a conversion table. The separation conditions basically include a gradation conversion curve (also referred to as gradation curve or tone curve, which is a characteristic of the halftone dot area ratio with respect to the original density), and in addition, highlights of this gradation conversion curve. Consists of conditions such as density value and halftone dot area ratio, shadow density value and halftone dot area ratio, color correction (color correction), UCR (undercolor removal), K plate (black plate) generation, sharpness, and the like. This decomposition condition is stored in each scanner as a default decomposition condition optimized empirically for each scanner.
[0009]
Therefore, when the document is a reversal document, the operator can convert the RGB signal from the RGB density signal based on the default decomposition condition that is empirically optimized for each scanner and stored in a storage unit such as a memory. Can be easily converted into a cmyk signal which is a half-tone signal.
[0010]
However, in an imaging apparatus such as a digital camera, which has only a captured RGB image signal, does not have an image processing function attached to the scanner, and does not have a default decomposition condition as a reference decomposition condition. The disassembly conditions cannot be determined with the same feeling as the scanner that the operator is familiar with on a daily basis.
[0011]
In such a case, the obtained RGB image signal is once taken into a computer, and an image processing tool (software) such as Adobe Photoshop (registered trademark) is used on the computer to perform default decomposition for each imaging device. A method of setting conditions is conceivable. Further, a method of converting to a cmyk image signal by trial and error while looking at an image on a monitor that does not necessarily match the appearance of a printed matter that is a hard copy can be considered.
[0012]
However, when there are various image input devices (including an image capturing device and an image reading device) such as a digital camera and a scanner, it is very difficult to set default disassembly conditions for each device, and time There is a problem that it takes.
[0013]
The present invention has been made in consideration of such a problem, and an RGB image signal photographed by a digital camera, an RGB image signal read by a scanner, or the like is empirically optimized for a specific scanner default. It is an object of the present invention to provide a color conversion method for converting to a cmyk signal using the decomposition conditions.
[0014]
[Means for Solving the Problems]
  The first aspect of the present invention is a process for converting a device-dependent image signal obtained for an original scene into a colorimetric value signal.In a color conversion method comprising,
  Converting the colorimetric signal into a dye density signal on a reversal document;
  The dye density signal includes gradation conversion processing, color correction processing, under color removal processing, and K plate generation processing.The default separation condition for a specific color scannerA process of converting to a cmyk signal by processing under standard image processing conditions;
  It is characterized by having.
[0015]
According to the first aspect of the present invention, the image signal of the original scene can be converted into a cmyk signal processed under standard image processing conditions as a result.
[0016]
  The second aspect of the present invention is a process of converting a device-dependent image signal obtained for an original scene into a first colorimetric value signal.In a color conversion method includingA process of converting the first colorimetric value signal into a first dye density signal on a reversal original, and a process of converting an image signal of a reversal original read by an image reading device into a second colorimetric value signal. A process of converting the second colorimetric value signal into a second dye density signal, and the first or second dye density signal.The floorIncludes tone conversion processing, color correction processing, undercolor removal processing, and K version generation processingThe default separation condition for a specific color scannerAnd a process of converting to a cmyk signal by processing under standard image processing conditions.
[0017]
According to the second aspect of the present invention, since the image signal of the reversal document read by the image reading device is once changed to a so-called device-independent colorimetric value signal, for example, the model of the image reading device is different. In other words, a cmyk signal processed under desired image processing conditions can be obtained without depending on the device.
[0018]
Further, the third aspect of the present invention is characterized in that, in the first or second aspect of the present invention, a setup process for performing setup processing on the colorimetric value signal is provided.
[0019]
According to the third aspect of the present invention, a desired finish depending on the image can be obtained.
[0020]
  Furthermore, the fourth invention is the dye concentration signal according to any one of the first to third inventions.The floorIncludes tone conversion processing, color correction processing, undercolor removal processing, and K version generation processingThe default separation condition for a specific color scannerIn the process of converting to a cmyk signal after processing under standard image processing conditions,
  Colorimetric values or signals that can be converted to colorimetric valuesThe color patches on the reversal manuscript arranged regularlyThe default separation condition for a specific color scannerProcessed under standard image processing conditions to create a cmyk image, and for the patch reproduced on the cmyk image,Colorimetric values or convertible with colorimetric valuesPick out the correspondence with the signal as a data set, make a table of the relationship between the two,The default separation condition for the specific color scannerA standard image processing condition table is created, and the dye density signal is converted into a cmyk signal using the image processing condition table and an interpolation method.
[0021]
  According to the fourth aspect of the present invention, the dye density signal cmy can be easily converted into a cmyk signal that is a halftone signal.
  For example, the image reading apparatus is a color scanner, and the standard image processing conditions are default disassembly conditions of the color scanner.
[0022]
DETAILED DESCRIPTION OF THE INVENTION
Hereinafter, an embodiment of the present invention will be described with reference to the drawings.
[0023]
FIG. 1 shows the overall configuration of an image signal processing apparatus 10 to which an embodiment of the present invention is applied.
[0024]
The image signal processing apparatus 10 includes a colorimetric conversion unit 16, and the colorimetric conversion unit 16 includes an RGB signal having image information of a scene (scene) captured by a digital camera (not shown) that is an imaging device. (Digital image signal) 12 is supplied, and is read by an image input unit (also referred to as a line scanning reading device or a scanner input unit) constituting a scanner (not shown) which is an image reading device as necessary. An RGB signal (digital image signal) 14 having image information of a reversal document is supplied. Note that the digital camera may be a digital video camera, and the present invention can be applied to any other device that outputs a digital image signal (digital image data) having an imaging function.
[0025]
Further, the scanner may be a line scanning reading device that is a color scanner equipped with a linear image sensor, or may be applied to a surface scanning reading device that is a color scanner equipped with an area image sensor.
[0026]
The colorimetric conversion unit 16 converts the RGB signal 12 of the digital camera into XYZ or L^*a^*b^*A colorimetric conversion matrix (hereinafter also simply referred to as a matrix) 20 for converting a colorimetric value signal (also referred to as an XYZ signal) 18 such as a colorimetric value signal, and an RGB signal 14 of the scanner input unit into an XYZ signal or L^*a^*b^*And a colorimetric conversion table 24 (also referred to as a colorimetric conversion lookup table, a lookup table, or simply a table) for conversion into a colorimetric value signal 22 such as.
[0027]
In general, the RGB signals 12 and 14 are referred to as device-dependent (also referred to as device-dependent or device-dependent) signals (data), whereas the colorimetric value signals 18 and 22 are dependent on the device. No (device-independent or device-independent) signal (data).
[0028]
The colorimetric value signal 18 of the digital camera includes a light source changing unit 25 that absorbs the difference between the photographing light source and the observation light source, and a colorimetric setup unit that sets up the highlight part density and the shadow part density (also simply referred to as a setup part). .) 26, faithfulness of the original scene constituting the color density conversion unit 30 that converts the colorimetric value signal 18 of the original scene (original photographed scene) into the pigment density signal 28 on the reversal document, and the cmyk conversion unit 32. It is supplied to the reproduction table 42 as necessary.
[0029]
On the other hand, the colorimetric value signal 22 of the scanner input unit constitutes the colorimetric setup unit 26 and a standard condition reproduction table unit 43, and converts the colorimetric value signal 22 into a dye density signal 34. To be supplied.
[0030]
The dye density signals (cmy signals) 28 and 34 output from the dye density conversion units 30 and 36 are selected through the switch 45 and supplied to the standard condition reproduction table 44.
[0031]
The cmyk conversion unit 32 basically includes an original scene faithful reproduction table 42 and a standard condition reproduction table unit 43. The original scene faithful reproduction table 42 is a look-up table that converts the colorimetric value signal 18 of the supplied digital camera into a cmyk signal 46 that is a halftone signal stored colorimetrically.
[0032]
The standard condition reproduction table unit 43 includes a dye density conversion unit 36 that converts the colorimetric value signal 22 into a dye density signal 34, and a switch that selects an output signal of either the dye density conversion unit 30 or the dye density conversion unit 36 ( Multiplexer and selection means) 45 and a standard condition reproduction table 44. The standard condition reproduction table 44 performs processing for converting any one of the dye density signals 28 and 34 selected by the switch 45 into a cmyk signal 48 which is a halftone signal. The cmyk conversion unit 32 has a so-called three-color four-color conversion function for converting a three-color color system signal into a four-color color signal.
[0033]
The cmyk signals 46 and 48 that are halftone signals output from the cmyk conversion unit 32 are supplied to the image output unit 35 that outputs an image based on the cmyk signals 46 and 48. The image output unit 35 employs a known configuration including, for example, a binarization conversion unit (not shown), a laser exposure scanning unit (image setter or the like), a developing unit, a plate making unit, and a printing unit. Can do. In the binarization conversion unit, the cmyk signal 46 or the cmyk signal 48 is compared with each CMYK threshold matrix selected according to the output conditions such as the number of screen lines and the screen angle, and binarization processing is performed. The laser exposure scanning unit exposes and scans the film with a laser beam that is turned on and off based on the binary signal (also referred to as a binary image signal) to form a latent image. In the developing unit, the film on which the latent image is formed is developed to visualize the image, thereby creating a plate-making film. The plate making section creates a plate from this plate-making film. In the printing section, the printing plate, in this case, four printing plates for CMYK, is mounted on a printing machine, and the four colors of ink attached to the printing plate are transferred to the main paper (printing paper), whereby the image is printed. A printed matter is produced as a hard copy on which is formed.
[0034]
The image output unit 35 does not require film development processing, and can directly simulate halftone dots, screen line numbers, and screen angles by printing them as halftone images on a sheet of paper and performing direct digital color proofing (DDCP). A system can also be used.
[0035]
The colorimetric setup unit 26 converts a colorimetric END (equivalent neutral density) signal 52 into a colorimetric END forward / inverse conversion matrix (also simply referred to as a matrix). 50, a rough data generation unit 54 that outputs a colorimetric END signal 56 generated by thinning out from the colorimetric END signal 52, and a highlight density signal 60 based on the generated colorimetric END signal 56 An auto setup unit 58 that automatically determines the shadow density signal 62, an END / END conversion unit 64 that performs gradation conversion processing on the colorimetric END signal 52 and converts it into a colorimetric END signal 66; A colorimetric END forward conversion / inverse conversion matrix 50 which reversely converts the converted colorimetric END signal 66 into the colorimetric value signal 22;
[0036]
Note that the image signal processing apparatus 10 is configured to be controlled by a computer (including a CPU, ROM, RAM, external storage device, monitor, and other input / output devices) that is not shown in the figure. Each block constituting the apparatus 10 includes not only hardware but also a part constituted by software. The computer functions as control, determination, calculation, comparison means, and the like.
[0037]
Next, the detailed configuration and operation of each block constituting the image signal processing apparatus 10 will be described.
[0038]
A matrix 20 that constitutes the colorimetric conversion unit 16 and converts the RGB signal 12 into the XYZ signal 18 is created based on the flowchart shown in FIG. In the following description, the matrix 20 is determined from the RGB color space as CIE-L.^*a^*b^*(Light source: auxiliary standard light CIE-D50) The conversion to a color space will be described as an example. In this case, CIE-L^*a^*b^*Mutual conversion between the color space and the XYZ color space can be performed uniquely by the following known formula (1). Therefore, in all the descriptions below, the XYZ color space (or CIE-L^*a^*b^*Processing in color space) is CIE-L^*a^*b^*It can be replaced with processing in a color space (or XYZ color space). In addition, it can be replaced with processing in a colorimetric color space equivalent to these.
[0039]
L^*= 116 (Y / Yn)^1/3-16
a^*= 500 {(X / Xn)^1/3-(Y / Yn)^1/3}
b^*= 200 {(Y / Yn)^1/3-(Z / Zn)^1/3} (1)
First, a color chart 70 which is a kind of color chart having a plurality of color patches 72 (see FIG. 2) in which representative colors are changed according to saturation, lightness, and hue is prepared (step S1). In this embodiment, a Macbeth color checker (registered trademark: manufactured by Macbeth A division Kollmorgen, USA) is used as a color chart. The Macbeth color checker is a color chart in which CIE (1931) xyY value, hue, Munsell notation value, and saturation are defined as is well known.
[0040]
Specifically, 24 colors
1. Dark skin
2. Light skin
3. Blue sky
4). Foliage
5). Blue flower
6). Bluish green
7). Orange
8). Purple green
9. Red with moderate brightness (moderate red)
10. Purple
11. Yellow green
12 Orange yellow
13. Blue
14 Green
15. Red
16. Yellow
17. Magenta
18. Cyan
19. White
20. Neutral 8 (neutral 8: light gray, 8 is 8 in Munsell notation)
21. Neutral 6.5 (neutral 6.5: light medium gray)
22. Neutral 5 (neutral 5: intermediate gray)
23. Neutral 3.5 (neutral 3.5: dark gray)
24. Black
It is.
[0041]
The color chart is not limited to the Macbeth color checker, and for example, a color chart that covers the color space almost uniformly, such as a JIS standard color target, can be used.
[0042]
Next, each color of the 24 colors of the color chart 70, that is, each of the 24 patches 72 is photographed using a digital camera under a CIE-D50 photographing light source, and the respective RGB signals 12 of the patches 72 are obtained. Each RGB signal 12 is converted into a luminance value and then raised to the 1/3 power (step S2). Note that the conversion to the luminance value can be obtained by, for example, canceling the γ correction performed in the device. Moreover, the 1/3 power is obtained by converting the obtained RGB signal 12 into CIE-L as understood from the above equation (1).^*a^*b^*This is for processing in the color system.
[0043]
Next, the RGB signal 12 for each patch 72 obtained in step S2, more specifically, for each patch 72, each of the R, G, and B values from the R, G, and B values that have undergone luminance conversion and are raised to the third power. Values R, G, B, R², G², B², RG, GB, BR (9 variables) are calculated (step S3).
[0044]
Next, in order to prevent a multi-co (multicollinearity) phenomenon when a multiple regression analysis described later is performed, a principal component analysis is performed on the 24-color data of 9 variables obtained in step S3. A nine-variable principal component score (principal component score) V is obtained. A principal component score V (V is considered as a vector) represented by the following equation (2) is obtained for each color.
[0045]
V = (v₁, V₂, V_Three, V_Four, V_Five, V₆, V₇, V₈, V₉) ... (2)
In this equation (2), the component v of the principal component score V₁, V₂, V_Three, V_Four, V_Five, V₆, V₇, V₈, V₉Are not correlated with each other and do not cause the multi-co phenomenon.
[0046]
Next, using a colorimeter (not shown), each patch 72 of the color chart 70, that is, the colorimetric values C (L^*a^*b^*) (Since C is also obtained for 24 colors, it is considered a vector.) (Step S5). The process of obtaining the colorimetric value C may be obtained at any time with respect to the processes of steps S2 to S4.
[0047]
Next, using the colorimetric value C as an objective variable (dependent variable) and the principal component score V as an explanatory variable (independent variable), a partial regression coefficient A (A is also considered as a vector) is obtained by multiple regression analysis (step S6). ).
[0048]
When this multiple regression analysis is performed, a weighting matrix (also simply referred to as a matrix) P = [Pi] (i = 1, 2,... 24) is operated (step S7). Pi is a weight of each color. For example, 20. Neutral 8 (neutral 8: light gray, 8 is 8 in Munsell notation), 21. Neutral 6.5 (neutral 6.5: light medium gray), 22. Neutral 5 (neutral 5: intermediate gray), 23. By increasing the weight for neutral 3.5 (neutral 3.5: dark gray), the color reproducibility for gray can be improved. It is possible to reproduce the skin color or sky blue more faithfully by limiting the weight to a desired color according to the purpose, for example, skin color for humans and sky blue for outside scenery. This is preferable. In addition, the weight is obtained by statistically processing the input RGB signal 12 and color the image signal for one screen, in other words, the pixel-divided image data in an area centered on the color of each patch of the color chart, for example. Automatic weighting processing is performed by dividing the space and weighting according to the frequency value of the pixels existing in each region, for example, by maximizing the weight of the mode value and decreasing the weight in the following order. Can do. Note that the color space may be divided into RGB color spaces, and XYZ (L^*a^*b^*) Color space may be used.
[0049]
In this case, the color with the increased weight is called the main color of the image. Note that the value of each element is normalized to a value obtained by dividing the value of each element by the total value of each element so that the total value of each element of the matrix P becomes 1.
[0050]
The multiple regression analysis of step 6 for obtaining the partial regression coefficient A will be described in detail.
[0051]
Between the colorimetric value (vector) C, which is the objective variable, the partial regression coefficient (vector) A to be obtained, and the principal component score (vector) V, for each colorimetric value C of 24 colors, It is assumed that the linear linear expression shown in Equation (3) is established.
[0052]
C = AV (3)
An expression using a matrix equivalent to the expression (3) is shown in an expression (4), and an expression using a sum symbol Σ is shown in an expression (5).
[0053]
[Expression 1]

[0054]
[Expression 2]

[0055]
It should be noted that an expression using the sum symbol Σ, such as the expression (5), can be changed to L^*= Σ (j = 1 → 10) a_1jv_jIt shall be expressed as
[0056]
The partial regression coefficient A is L^*a^*b^*Are separately determined using the least square method. For example, L^*For e according to the following equation (6)_LThe partial regression coefficient a that minimizes_1jYou can ask for.
[0057]
e_L= Σ (i = 1 → 24) Pi {Ci−Σ (j = 1 → 10) a_1jv_j}²(6)
Here, i is a patch number of the color chart 70, Pi is a weight of each color, and j is a variable number (1, 2,..., 10).
[0058]
When Expression (6) is expressed by a vector and a matrix, Expression (7) is obtained. However, in the equation (7), the colorimetric value C and the partial regression coefficient a are vectors, and the principal component score [V] and the weighting matrix [P] are matrices. t represents the transposition of the matrix.
[0059]
e_L=
(Vector C-vector a [V])^t[P] (vector C-vector a [V]) (7)
Hereinafter, the vector C is simply written as C, and the vector a is simply written as a. First, equation (7) can be modified as follows.
[0060]
e_L= (C^t-[V]^ta^t)^t[P] (C^t-[V]^ta)
Here, if a and [V] are interchanged,
e_L= (C−a [V]) [P] (C^t-[V]^ta).
[0061]
Generally, (ABC)^t= C^tB^tA^tSo
e_L= C [P] C^t+ A [V] [P] [V]^ta^t-A [V] [P] C^t-C [P] [V]^ta^tIt becomes.
[0062]
Here, [V] [P] [V]^t= [N], [V] [P] C^t= [U]
e_L= C [P] C^t+ A [N] a^t-2a [U] (8)
It becomes. In equation (8), e_LTo minimize the derivative for each element of the partial regression coefficient a must be equal to zero. Therefore, the following equation (9) is established.
[0063]
[Equation 3]

[0064]
From this equation (9), the partial regression coefficient a can be obtained from equation (10) as shown below.
[0065]

([]^-1Indicates an inverse matrix. )
In this way L^*Partial regression coefficient a (a_1j) Are obtained, and in the same manner, the remaining a^*, B^*Also for partial regression coefficient a_2j, A_3j10 each can be obtained. Obtained partial regression coefficient a_1j, A_2j, A_3jIs 3 × 10 {see equation (4)}.
[0066]
After all, L^*a^*b^*When the partial regression coefficients A are collectively expressed, they can be expressed by equation (11).
[0067]
A = CPV^t(VPV^t)^-1              ... (11)
Next, each objective variable L^*a^*b^*The partial regression coefficient A obtained for the test is tested with 5% significance, for example (step S8), and the explanatory variable V having 5% significance (95% reliability) is stored in the storage means. 1 is stored as the matrix 20 shown in FIG. 1 (step S9). In addition, you may test | inspect by 1% significance (99% of reliability).
[0068]
In the test of step S8, as shown in the equations (12) to (16), the regression sum of squares S_R{The estimated value L^*i (the symbol “^” is added above L) and the average L of the colorimetric values^*i (the symbol “^-" ) Sum of squares of difference: see equation (12)} and residual sum of squares S_E{Colorimetric value L^*i and estimated value L^*sum of squares of differences of i: (see equation (13)), and the regression sum S_RAnd residual sum of squares S_EUnbiased variance V_R, V_E{Refer to (14) and (15)} respectively. Further, the partial F value Fj {refer to equation (16)} is obtained. In the equation (16), “^” added to the head of a is the estimated value of the coefficient of the matrix on the right side of the equation (4), S^jjIs the explanatory variable v_jMeans the diagonal term of the inverse of the variance-covariance matrix.
[0069]
[Expression 4]

[0070]
The value obtained by referring to the F distribution in which the partial F value Fj is 5% significant = F ′_np-1(0.05) = F ′_24-9-1(0.05) = F ′₁₄If it is greater than (0.05) = 4.6001, the regression is significant at the 5% significance level, and is determined to be a partial regression coefficient useful for prediction, and stored in the matrix 20 (see FIG. 1).
[0071]
In this way, the RGB signal 12 (see FIG. 1) obtained by the digital camera is converted into an XYZ colorimetric value signal 18.^*a^*b^*Can be created. An expression related to this matrix 20 is shown in Expression (17).
[0072]
L^*= Σa_1jv_j  (J: 5% significant variable)
a^*= Σa_2jv_j  (J: 5% significant variable)
b^*= Σa_3jv_j  (J: 5% significant variable) (17)
In this equation (17), the meaning of (j: 5% significant variable) means that only the variable that is 5% significant by the test is used as a variable. This is different from the above-described equation (5). A_1j, A_2j, A_3jAre elements of the partial regression coefficient matrix A and v_jIs the principal component score of the explanatory variable obtained in step S4.
[0073]
By using this matrix 20, colorimetric conversion can be performed using an existing color chart 70. In this case, since the lookup table is not used, the memory capacity is small. In other words, the colorimetric conversion can be performed with high accuracy even when the memory capacity is small. Further, by using the weighting matrix P, color conversion can be performed with high accuracy by limiting to the desired color (main color) of the image. If the number of explanatory variables is nine, the matrix corresponding to the equation (17) may be created using all the partial regression coefficients A without performing the test process in step S8.
[0074]
The CIE-L obtained in this way^*a^*b^*The value of the color space is converted into the value of the XYZ color space by the above-described equation (1), and is set as a colorimetric value signal 18 that is an output signal of the matrix 20 (see FIG. 1).
[0075]
Next, the colorimetric value signal 18 is colorimetrically changed by the light source changing unit 25 as necessary, and becomes a new colorimetric value signal 18 (the same code is used). If the observation light source is a photographing light source (same as CIE-D50) at the time of photographing with a digital camera, this light source changing process is unnecessary.
[0076]
When the colorimetric value XYZ is changed to a new colorimetric value XYZ (sign is X′Y′Z ′) after the light source changing process by the light source changing unit 25, the conversion is performed by the following equation (18). , (19) and (20).
[0077]
X '= XX₂/ X₁                                ... (18)
Y '= YY₂/ Y₁                                ... (19)
Z '= ZZ₂/ Z₁                                ... (20)
In the equations (18) to (20), XYZ is a colorimetric value under the photographic light source, X′Y′Z ′ is a colorimetric value under the observation light source, and X₁Y₁Z₁Is the white point of the light source, X₂Y₂Z₂Is the white point of the observation light source.
[0078]
Hereinafter, it is assumed that the colorimetric value signal 18 is a signal subjected to light source change processing as necessary.
[0079]
On the other hand, the RGB signal 14 output from the scanner input unit is converted into a colorimetric value signal 22 which is an XYZ signal by a colorimetric conversion table 24.
[0080]
For example, the colorimetric conversion table 24 prepares a color reversal document having a total of 13 × 13 × 13 = 2197 color patches in which the density of each of the cmy color patches is regularly distributed in 13 steps. The color patches constituting the color reversal document are read by the scanner input unit and read by the colorimeter. Then, a correspondence relationship between the RGB values read by the scanner and the XYZ values read by the colorimeter is obtained as a table. Values between reading values not existing in this table are obtained by interpolation processing.
[0081]
Then, the RGB signal 14 that is the output signal of the scanner input unit is converted into the colorimetric value signal 22 by the colorimetric conversion table 24, and the RGB signal 12 that is the output signal of the digital camera is measured by the colorimetric conversion matrix 20. By converting the color value signal 18, a so-called auto setup process can be performed by commonly using a colorimetric setup unit 26 described below. That is, there is an advantage that the auto setup processing software can be used in common by once converting the colorimetric value signals 18 and 22 into the colorimetric END.
[0082]
Next, the colorimetric setup unit 26 will be described with reference to the flowchart of FIG.
[0083]
Conventionally, tone conversion has been performed in a density space because it is easy to understand intuitively, and image processing such as tone conversion processing and color correction processing in currently available color scanners is also performed in density. It is implemented with reference to the signal. Therefore, first, the colorimetric value signal 18 or the colorimetric value signal 22 is converted into a cmy END (equivalent neutral density) signal 52 which is a colorimetric system of the density signal (step S11). As shown in FIG. 3, this conversion process includes a linear conversion process from the XYZ color system to the RGB color system (step S11a) and a non-linear conversion process from the RGB color system to the cmy color system (step S11b). It is a two-stage process. The END signal converted from the colorimetric value signal 18 or the colorimetric value signal 22 is referred to as a colorimetric END signal 52.
[0084]
For example, when converting the colorimetric value signal 22 of the reversal document into the colorimetric END signal 52, the color reproduction area 71 (area shown by hatching) of the reversal film is drawn on the CIE chromaticity diagram shown in FIG. The coordinates on the chromaticity diagram of each point of the three primary stimuli RGB (respectively Rxyz, Gxyz, and Bxyz) of the area 73 including the color reproduction area 71, that is, the chromaticity coordinates are represented by the equations (21) to (23). The coordinates shown in the equation are assumed. In this case, the region 73 including the color reproduction region 71 of the reversal film is a triangular region having vertices Rxyz, Gxyz, and Bxyz as vertices on the chromaticity diagram shown in FIG.
[0085]
Rxyz = Rxyz (x_R, Y_R, Z_R... (21)
Gxyz = Gxyz (x_G, Y_G, Z_G... (22)
Bxyz = Bxyz (x_B, Y_B, Z_B... (23)
In addition, on the chromaticity diagram, the coordinates of the basic stimulus (white stimulus) Wxyz in the XYZ color system are shown in equation (24).
[0086]
Wxyz = Wxyz (x_W, Y_W, Z_W(24)
In this case, the color measurement value signal 22 of the XYZ color system (the right side matrix on the right side) is converted into the RGB color system color signal RGB (through the conversion matrix (the left side matrix on the right side)) by the following equation (25). Left side matrix) (step S11a).
[0087]
[Equation 5]

[0088]
Cofactor A_ijIs obtained by the following equation.
[0089]
A_ij= (-1)^{i + j}D_ij
However, D_ijIs a small determinant excluding i rows and j columns.
[0090]
The cmy value of the colorimetric END signal 52 is the basic stimulus (white stimulus) W_XYZR, G, B values Rw, Gw, Bw obtained by substituting the coordinates of (25) into the expression (25) The ratios of the R, G, B values obtained by the expression (25) to the expressions (26) to ( It can be obtained by logarithmic transformation as shown in equation (28).
[0091]
c = -log (R / Rw) (26)
m = −log (G / Gw) (27)
y = -log (B / Bw) (28)
In order to perform the setup process in a short time, the rough data creation unit 54 does not target, for example, the colorimetric END signal (actually digital data) 52 for one document, but from the operator in the document. A so-called rough scan process is performed, such as selecting data for only a portion of an image whose area is designated through the CPU, or creating data by thinning out the data when the image exists on the entire surface of the document.
[0092]
Next, the auto setup unit 58 performs auto setup processing based on the colorimetric END signal 56 that is the rough data selected by the rough data creation unit 54 (step S12). In this auto setup process, for example, as well known in Japanese Patent Laid-Open No. 2-105676, a histogram is created for the colorimetric END signal 56 and then a cumulative histogram is created.
[0093]
In FIG. 1, for example, 0.1% point data (HL concentration) of the cumulative histogram, as shown in a characteristic diagram (drawing in FIG. 1) showing the operation of the END / END conversion unit 64. The value of the colorimetric END signal 56 corresponding to D1 is set to DH, and the value of the colorimetric END signal 56 corresponding to 98% point data (SD density) D2 is set to DS (step S12).
For example, when the D1 density is D1 = 1.0 and the D2 density is D2 = 2.0, the DH corresponding to the D1 density is set to DH = 0.1. The DS corresponding to the density of D2 is set to DS = 3.0. In practice, the 0.1% point data (HL density) D1 is converted to a density corresponding to 0% of the halftone percentage, and the 98% point data (SD density) D2 is a density corresponding to 100% of the halftone percentage. Is converted to
[0094]
From the equation of the straight line 74 set in this way, the END / END conversion unit 64 converts all the calorimetric END signals 52 of the main scan data (data subjected to rough scan processing) into the calorimetric END signal 66. (Step S13). That is, the tone characteristics can be converted, in other words, tone correction processing can be performed. The conversion formula may be a formula of a curve passing through the highlight point 75 (D1, DH) and the shadow point 76 (D2, DS). The automatic determination of the HL concentration D1 and the SD concentration D2 (or the setup points 75 and 76) based on a predetermined condition as described above is referred to as auto-setup, but the HL concentration D1 and the SD concentration D2. Can be determined manually, or DH and DS can be manually corrected after auto-setup. The operation determined by the manual is called manual setup.
[0095]
Next, the CMY value of the colorimetric END signal 66 after the END / END conversion is inversely converted into the colorimetric value signal 22 which is an XYZ value by the colorimetric END inverse conversion matrix 50 (step S14).
[0096]
The inverse conversion process includes a process of converting the cmy value of the colorimetric END signal 66 into an RGB value of the RGB color system (step S14a), and a colorimetric value signal that takes the RGB value of the XYZ color system. 22 (step S14b). The process of converting the cmy value of the colorimetric END signal 66 into the RGB value of the RGB color system includes the following equations (29) to (31) obtained by solving the equations (26) to (28) for RGB. The process of converting RGB values to XYZ values of the XYZ color system is performed by the following equation (32) obtained by solving the above equation (25) for the matrix XYZ.
[0097]
R = Rw10^-c                                  ... (29)
G = Gw10^-m                                  ... (30)
B = Bw10^-y                                  ... (31)
[0098]
[Formula 6]

[0099]
The above-described setup process has been described by taking the colorimetric value signal 22 of the reversal document as an example, but can also be applied to the colorimetric value signal 18 of the digital camera.
[0100]
That is, generally, in the printing plate making field, a reversal manuscript is used as an original manuscript. The original scene is exposed on the reversal original, and the printing is finished in the color reproduction area 71 that is colored on the reversal original. This suggests that the color reproduction range 71 where the reversal document is colored is sufficient as the color reproduction range of the image signal to be handled.
[0101]
As described above, the colorimetric setup unit 26 linearly converts the colorimetric XYZ into the color system of the three primary colors of RGB, and calculates the ratios of R, G, and B obtained for the basic stimuli (light sources) Rw, Gw, and Bw as logarithms. By performing the conversion, the END value is obtained, and auto setup is performed using the END value. The END / END changed END value is subjected to inverse logarithmic conversion processing to obtain the changed RGB value, and the colorimetric value XYZ changed by the inverse transformation matrix is obtained. For this reason, the effect that the auto setup process based on the existing so-called density signal can be performed is achieved. Needless to say, the setup process can be performed as a manual setup process.
[0102]
Next, the dye density conversion unit 30 will be described with reference to the flowcharts of FIGS.
[0103]
The dye density conversion unit 30 converts a colorimetric value signal (XYZ) 18 of an original scene photographed by a digital camera into a dye density signal (cmy) 28 on a reversal document.
[0104]
For this conversion process, there are two possible methods: the method shown in FIG. 5 and the method shown in FIG.
[0105]
That is, in the method shown in FIG. 5, the colorimetric value signal (XYZ) 18 of the original scene is converted into the dye density signal (cmy) 81 of the original scene using a look-up table described later (step S21: first stage processing) In this method, the dye density signal (cmy) 81 of the original scene is converted into a reversal (simply also referred to as RV) 28 dye density signal (cmy) 28 on the original by a matrix described later (step S22: second stage). Processing). In the method shown in FIG. 6, the colorimetric value signal (XYZ) 18 of the original scene is converted into the colorimetric value signal (XYZ) 82 on the reversal document by a matrix described later (step S31: first stage processing) ), A colorimetric value signal (XYZ) 82 on the reversal document is converted into a dye density signal (cmy) 28 on the reversal document by a look-up table described later (step S32: second stage process).
[0106]
First, details of the method shown in FIG. 5 will be described with reference to the flowchart of FIG. In FIG. 5, in order to convert the colorimetric values XYZ relating to the processing in step S21 into the dye density cmy, a color chart covering the color space substantially uniformly, for example, each measurement of 24 colors of the Macbeth color checker 70 described above. The color value XYZ is obtained by a colorimeter (step S21a: see FIG. 7).
[0107]
Next, the dye density value cmy is calculated from the obtained colorimetric value XYZ (step S21b), and a lookup table of the dye density value cmy with respect to the colorimetric value XYZ is created. This look-up table is set and stored in the dye density conversion unit 30 as a look-up table in step S21 which is the first stage processing in the dye density conversion unit 30. In the process of step S21, this lookup table is used, and for the colors other than the 24 colorimetric values XYZ, the dye density value cmy is calculated by interpolation using the 24 colorimetric values XYZ. can do.
[0108]
The calculation process in step S21b, that is, a process for obtaining the dye density cmy from the colorimetric values (tristimulus values) XYZ will be described.
[0109]
The following equations (33) to (36) are established between the colorimetric value XYZ and the dye density cmy.
[0110]
X = k∫visP (λ) T (λ) x (λ) dλ (33)
Y = k∫visP (λ) T (λ) y (λ) dλ (34)
Z = k∫visP (λ) T (λ) z (λ) dλ (35)
T (λ) = 10^-h                                      ... (36)
However,
h = {cDc (λ) + mDm (λ) + yDy (λ) + base (λ)}
k = 100 / ∫visP (λ) y (λ) dλ (λ is the wavelength of light)
∫vis: Definite integral in the visible wavelength region (380 nm to 780 nm)
P (λ): spectral characteristic data of observation light source
T (λ): Spectral transmittance data of pigment of transmission object
{(36) is assumed that Lambert-Beer's Law holds}
x (λ), y (λ), z (λ): color matching function
Dc (λ), Dm (λ), Dy (λ): spectral density data of the cmy dye
base (λ): film-based spectral density data
When obtaining the dye density cmy from these equations (33) to (36), an inverse function may be obtained, but it cannot be obtained directly. Therefore, a successive approximation method such as the well-known Newton (Newton Raphson) method (for example, Ota Noboru “Color Engineering” pp254-260, Tokyo Denki University Press, December 20, 1993, 1st edition, 1st edition) is used. That's fine. The Newton-Raphson method (abbreviated as N · R method) will be briefly described with reference to the above reference book.
[0111]
When a general equation y = f (x) is Taylor-expanded with x = x0 close to the root of f (x) = 0 and only a first-order term is taken, f (x0 + Δx) with respect to a small change Δx of x ) = F (x0) + f ′ (x0) · Δx. However, f ′ (x0) is obtained by substituting x = x0 for the differential coefficient f ′ (x) of f (x). Therefore, a more correct value x1 of f (x) = 0 is obtained as f (x0 + Δx) = 0, and x1 = x0 + Δx = x0−f (x0) / f ′ (x0). As shown in FIG. 8, this corresponds to drawing a tangent line 83 at a point (x0, y0) on the function y = f (x) and obtaining an intersection point x1 between the tangent line 83 and the x axis.
[0112]
In order to apply this to the equations (33) to (36), after substituting the equation (36) into the equations (33) to (35), if a certain function fx, fy, fz is used, (33 Equations (35) to (35) can be expressed as equations (37) to (39).
[0113]
X = fx (c, m, y) (37)
Y = fy (c, m, y) (38)
Z = fz (c, m, y) (39)
In these equations (37) to (39), the initial values are c0, m0, y0, and the tristimulus values at that time are X0, Y0, Z0. Now, if the tristimulus values are changed by ΔX, ΔY, ΔZ when the minute changes Δc, Δm, Δy are added to c0, m0, y0, the following equation (40) is obtained.
[0114]

However, for example, ∂fx / ∂c represents a partial differential coefficient with respect to c of the function fx.
[0115]
If this formula (40) is arranged, formula (41) is obtained. Similarly, ΔY and ΔZ are also obtained as shown in equations (42) and (43).
[0116]
ΔX = Δc · ∂X / ∂c + Δm · ∂X / ∂m + Δy · ∂X / ∂y (41)
ΔY = Δc · ∂Y / ∂c + Δm · ∂Y / ∂m + Δy · ∂Y / ∂y (42)
ΔZ = Δc · ∂Z / ∂c + Δm · ∂Z / ∂m + Δy · ∂Z / ∂y (43)
Expressions (41) to (43) are displayed in a matrix as shown in Expression (44).
[0117]
(Q) = (J) (P) (44)
However, (Q) is a matrix of 1 row 3 columns whose elements are ΔX, ΔY, ΔZ in order from the first row, (J) is a Jacobian matrix whose partial differential coefficients are 3 rows 3 columns, and (P) is This is a matrix of 1 row and 3 columns in which elements are Δc, Δm, and Δy in order from the first row.
[0118]
On both sides of equation (44), the inverse matrix (J) of the Jacobian matrix (J)^-1Is multiplied by (45).
[0119]
(P) = (J)^-1(Q) ... (45)
Therefore, if the initial values c0, m0, and y0 are corrected to c1, m1, and y1, respectively, as shown in equation (46), a more correct approximate value can be obtained.
[0120]
c1 = c0 + Δc
m1 = m0 + Δm
y1 = y0 + Δy (46)
By repeatedly calculating using the Jacobian matrix (J) obtained as described above, the dye density signal cmy for the colorimetric value XYZ that is an arbitrary target value can be obtained. By performing the same processing for all target values on the grid in the XYZ color space, inverse conversion for converting the colorimetric value signal (XYZ) 18 for the original scene into the dye density signal (cmy) 81 for the original scene. A table is generated and stored as a look-up table in the dye density conversion unit 30 for the first stage processing (step S21).
[0121]
Next, the process of creating a matrix for the second stage process in step S22 performed following the process of the first stage will be described. The matrix creation process is a process including the multiple regression analysis process already described with reference to FIG. 2, and will be briefly described.
[0122]
First, c, m, y, c for the dye density cmy calculated in step S21b.², M², Y², Cm, my, yc, 9 variables and 24 colors are calculated (step S22a).
[0123]
Next, the principal component analysis is performed on the data of the nine variables and 24 colors, and the principal component score of the nine principal components is obtained (step S22b).
[0124]
On the other hand, colorimetric values for the 24 colors of Macbeth color checker exposed on the reversal film are obtained by a colorimeter (step S22c).
[0125]
Next, the dye density cmy (also referred to as RVcmy) on the reversal film is obtained by the N · R method described above (step S22d).
[0126]
Next, each of the pigment concentrations RVcmy on the reversal film (RVc, RVm, RVy) is set as an objective variable, and the principal component score (including a constant term) obtained in step S22b is used as an explanatory variable to perform 3 × by multiple regression analysis processing. Ten partial regression coefficient matrices are obtained (step S22e).
[0127]
When performing the multiple regression analysis process, as described above, a weighting matrix corresponding to 1: 1 may be applied to the 24 colors that are the target variable group (step S22f).
[0128]
Next, the partial regression coefficients obtained for the respective objective variables RVc, RVm, and RVy are tested with, for example, 5% significance (step S22g), and 5% significant explanatory variables are stored in the storage means, and partial regression is performed. The coefficients are stored in the dye density conversion unit 30 as a matrix for the second stage processing (step S22h). In this case, all the coefficients may be used without performing the test.
[0129]
In this way, the dye density conversion unit 30 that converts the colorimetric value signal (XYZ) 18 of the original scene into the dye density signal cmy on the reversal film is constructed.
[0130]
Next, details of the method shown in FIG. 6 will be described with reference to the flowchart of FIG. The process shown in FIG. 9 is also a process including the multiple regression analysis process already described with reference to FIG. 2 and will be briefly described.
[0131]
First, the 24 color measurement values XYZ of the Macbeth color checker 70 are obtained by a colorimeter (step S31a: see FIG. 9).
[0132]
Next, for the obtained colorimetric values XYZ, X, Y, Z, X², Y², Z², XY, YZ, ZX, 9 variables, 24 colors of data are calculated (step S31b).
[0133]
Next, the principal component analysis is performed on the data of the nine variables and 24 colors, and the principal component score of the nine principal components is obtained (step S31c).
[0134]
On the other hand, colorimetric values XYZ for the 24 colors of Macbeth color checker exposed on the reversal film are obtained by a colorimeter (step S31d). Note that the Macbeth color checker exposed to the reversal film is an image obtained by exposing the reversal film by photographing the Macbeth color checker under a predetermined light source using an optical camera equipped with the color reversal film. Reversal manuscript obtained by developing a reversal film.
[0135]
Next, each of the colorimetric values XYZ on the reversal document is set as an objective variable (RVX, RVY, RVZ), and the principal component score (including a constant term) obtained in step S31c is used as an explanatory variable to obtain 3 by multiple regression analysis processing. A partial regression coefficient matrix of x10 is obtained (step S31e). When performing the multiple regression analysis process, as described above, a weighting matrix corresponding to 1: 1 may be applied to the 24 colors that are the target variable group (step S31f).
[0136]
Next, the partial regression coefficients obtained for each objective variable RVX, RVY, and RVZ are tested with 5% significance, for example (step S31g), and 5% significant explanatory variables are stored in the storage means, and the partial regression coefficients are stored. Is stored in the dye density conversion unit 30 as a first-stage processing matrix (step S31h, step S31). Even in this case, all the coefficients may be used without performing the test.
[0137]
Next, the processing table for the second stage of step S32 is created by the same method as the method of step S21b or step S22d described with reference to FIG.
[0138]
In this way, the dye density conversion unit 30 is constructed based on the process of FIG. 6 for converting the colorimetric value signal (XYZ) 18 of the original scene into the dye density signal cmy on the reversal document.
[0139]
Next, the process of creating the original scene faithful reproduction table 42 that constitutes the cmyk conversion unit 32 will be described based on the flowchart shown in FIG.
[0140]
The original scene faithful reproduction table 42 is a look-up table for converting the colorimetric value signal (XYZ) 18 into a cmyk signal (cmyk data) 46 which is each halftone data stored colorimetrically.
[0141]
The creation process of the original scene faithful reproduction table 42 is the same process as the creation process described in the specification of Japanese Patent Application No. 8-154484 related to the applicant's application.
[0142]
In this case, first, a cmyk color chart having color patches in which the density and mixing ratio of cmyk changes stepwise by giving a plurality of cmyk mesh% data having regular intervals to the image output unit 35. Create (step S41).
[0143]
In this case, for example, the regular interval of cmyk is the total number of color patches if the interval of 6 colors is increased in increments of 20%, such as 0, 20,. Is 4⁶= 1296.
[0144]
Next, each patch of the cmyk color chart created by the image output unit 35 is measured using a colorimeter (step S42), and a colorimetric value (stimulus value) XYZ is obtained from the colorimetric data, and the cmyk network A conversion table (referred to as a forward conversion table) from% data to colorimetric value XYZ data is created (step S43).
[0145]
Since this forward conversion table is also used as a table for interpolation, the regular interval should be finer when the accuracy of interpolation is taken into account, but the color measurement work becomes enormous accordingly. Therefore, what is necessary is just to determine by the so-called trade-off between the complexity of the color measurement work, the accuracy, and the processing time of the computer described below.
[0146]
By the way, in the original scene faithful reproduction table 42, an input arbitrary colorimetric value signal (colorimetric value XYZ data, stimulus value data XYZ, or simply XYZ) 18 to a corresponding cmyk signal (cmyk). (They are also referred to as halftone data, color data cmyk, or simply cmyk.) 46 is obtained by increasing the number of variables from three to four. There may be multiple solutions. In order to solve this problem, it is necessary to establish a relationship between the three variables. Here, among the cmyk halftone data, the maximum value Kmax (k = Kmax) that the image output unit 35 can take for the color data k is used. It is fixed (step S44). The maximum value Kmax is, for example, the value of k in the cmyk dot% data is 100%.
[0147]
Then, the arbitrary XYZ values that are the three variables are converted into the corresponding cmy values (k is fixed) that are the three variables (step S45).
[0148]
In this case, the value of cmyk with respect to target values X0, Y0, Z0 which are arbitrary XYZ values (here, since k is fixed at k = Kmax, k is a constant, and in that sense, three variables of cmy). When c0, m0, y0, and k0 = Kmax are obtained, the partial regression coefficient of the regression equation is obtained using the above-described forward conversion table when k0 = Kmax.
[0149]
The regression equation at this time is such that a matrix of 3 rows and 4 columns that is a coefficient term is A, a matrix of 1 row and 3 columns of colorimetric values XYZ is fixed to T, k, in other words, a constant, If the remaining 1 cm 4 matrix of cmy is expressed as D, it can be expressed as the following equation (47).
[0150]
T = AD (47)
This equation (47) represents the relationship of the following equation (48).
[0151]
[Expression 7]

[0152]
In the equations (47) and (48), “1” constituting the matrix D is a value set to give a constant term in the three-dimensional plane of cmy.
[0153]
The coefficient A in the equation (48) can be obtained by the multiple regression analysis described above by substituting each data set of the above-described forward conversion table obtained when k = Kmax.
[0154]
Next, using the regression equation obtained when k = Kmax, the value of cmy (k is set as Kmax) corresponding to the target values X0, Y0, Z0 by the N · R method described above. Can be requested.
[0155]
Next, it is determined whether or not the obtained cmy value is a value (color) within the reproduction range of the image output unit 35 (step S46). That is, does the data cmy obtained with the reproducible minimum density of halftone data (color data) Cmy as Cmin, Mmin, Ymin and the maximum density as Cmax, Mmax, Ymax satisfy the relationship of the following equation (49)? Determine if.
[0156]
Cmin ≦ c ≦ Cmax
Mmin ≦ m ≦ Mmax
Ymin ≦ y ≦ Ymax (49)
When the obtained values of cmy satisfy the equation (49), the color data cmyk for the target values (X0, Y0, Z0) of the stimulus value data XYZ are respectively c = csol, m = msol, y = Ysol, k = ksol (in this case, ksol = Kmax), an inverse conversion table is created and set as the original scene faithful reproduction table 42 (step S47). In the following description, a data set of (csol, msol, ysol, ksol) is referred to as “cmyksol” as necessary.
[0157]
If the value of the cmy data obtained in the process of step S45 is a value that does not satisfy the equation (49), the color data k fixed at k = Kmax is k = k−Δk, in this case k = Kmax. After setting −Δk (step S48), the process of step S45 is repeated within a range where the color data k does not become smaller than the predetermined minimum value k = Kmin (step S49). Note that the minute change Δk is an arbitrary data interval of the color data k constituting the first inverse conversion table. For example, when the color data k is set as data in the range of 0% to 100%, What is necessary is just to set to 1% increment etc., and if the color data k is set as data of the range of 0-255, it should just set to 1 value increment. When the second processing of step S45 is performed, the XYZ value that is the colorimetric value of the left side of the equation (48) for the color data k = Kmax-1 = 100-1 = 99 is k = Kmax = 100%. And k = 80%, and can be obtained by interpolation processing (interpolation processing) of each XYZ data stored as a forward conversion table.
[0158]
On the other hand, if it is determined in step S49 that Kmin> k, the color data cmyk for the target value (X0, Y0, Z0) is designated as data outside the color gamut of the image output unit 35. Here, The calculation of the color data cmyksol is not performed (step S50).
[0159]
By performing all the stimulus value data XYZ18 supplied to the original scene faithful reproduction table 42 as the target value (X0, Y0, Z0) by the above processing, halftone% which is color data in the color gamut of the image output unit 35 is obtained. For the stimulus value data XYZ18 from which the data cmyk46 can be obtained, the color data cmyksol when the maximum color data k is given can be obtained (step S51).
[0160]
The color data k in step S44 may be fixed to the minimum value Kmin (k = Kmin) instead of the maximum value Kmax. In this case, the process of step S48 is set to k = k + Δk, The process of S49 may be set to Kmax <k. Also, it is easily conceivable to set the color data k to an arbitrary value. For example, in the process of step S48, k = k−Δk and k = k + Δk may be alternately performed for k set to an arbitrary value, and the process of step S49 may be performed corresponding to the process of step S48. .
[0161]
Further, the color data cmyk designated as the out-of-gamut data in steps S49 and S50 is not the gist of the present invention and will not be described in detail. However, as described in Japanese Patent Application No. 8-154584, so-called An inverse conversion table of the color data cmyk for the target values (X0, Y0, Z0) can be created by compression processing and clipping processing of the color data CMYK by gamut mapping processing.
[0162]
In this manner, the correspondence table of the cmyk signals 46 obtained for all target values (X0, Y0, Z0) constituting the colorimetric value signal 18 is stored as the original scene faithful reproduction table 42. Any colorimetric value signal 18 can be converted into a cmyk signal 46 within the color reproduction range of the image output unit 35.
[0163]
Next, creation of the standard condition reproduction table 44 constituting the cmyk conversion unit 32 will be described.
[0164]
In the standard condition reproduction table 44, as shown in FIG. 1, the RGB signal 12 captured by the digital camera is converted into the colorimetric value signal (XYZ) 18 by the matrix 20, and the colorimetric value signal XYZ is converted to the dye density conversion. It is a table for converting the dye density signal (cmy) 28 after being converted into a dye density signal (cmy) 28 by the section 30 and converting the dye density signal (cmy) 28 after the conversion into a cmyk signal 48 which is a halftone signal.
[0165]
In the standard condition reproduction table 44, the RGB signal 14 captured by the scanner is converted into the colorimetric value signal (XYZ) 22 by the table 24, and the colorimetric value signal (XYZ) 22 after this conversion is converted into the dye density conversion unit. 36 is a table for converting the dye density signal (cmy) 34 after being converted into a dye density signal (cmy) 34 by 36 and converting the dye density signal (cmy) 34 after this conversion into a cmyk signal which is a halftone signal. In this case, the dye density conversion unit 36 can be created in the same process as the method of steps S21b and S22d in FIG. 7 shown in the creation process of the dye density conversion unit 30 described above. Description is omitted.
[0166]
The process for creating the standard condition reproduction table 44 will be described with reference to the flowchart shown in FIG.
[0167]
First, a color chart having 13 × 13 × 13 color patches in which the dye density cmy is arranged in a grid pattern on a reversal document is created (step S61). This color chart has each color patch in which c (cyan) m (magenta) y (yellow) is swung in 13 steps from the minimum density to the maximum density.
[0168]
Next, the color chart is decomposed under the default decomposition conditions of the scanner, in other words, the color chart which is a transparent original is read by the scanner and converted into digital data (step S62). The default separation conditions of the scanner include at least gradation characteristic conversion processing, color correction processing, UCR (under color removal) processing, and K plate generation processing. However, when disassembling, the minimum density of the gray color patch is such that the cmyk mesh% is all 0%, and the maximum gray density of the color patch is the so-called solid cmyk mesh%. Set.
[0169]
Then, a conversion table (lookup table) is created by associating the dye density of each color patch with the cmyk halftone value obtained when the scanner is decomposed under the default separation conditions of the scanner (step S63). 44 is set. In practice, an arbitrary dye density signal (cmy) 28 or dye density signal (cmy) 34 from the minimum density to the maximum density is converted into a desired cmyk net% signal 48 by this standard condition reproduction table 44 and interpolation processing. Can be converted.
[0170]
As described above in detail, according to the above-described embodiment, the RGB signal 12 of the scene (original scene) photographed by the digital camera is converted into the colorimetric value signal (XYZ) 18 by the matrix 20 constituting the colorimetric conversion unit 16. The colorimetric value signal (XYZ) 18 is subjected to setup processing in the colorimetric setup unit 26 as necessary. The signal (cmy) 28 is converted.
[0171]
Further, the RGB signal 14 obtained from the reversal original by the scanner input unit is converted into a colorimetric value signal (XYZ) 22 by the table 24 constituting the colorimetric conversion unit 16, and this colorimetric value signal (XYZ) 22 is converted into this colorimetric value signal (XYZ) 22. On the other hand, the colorimetric setup unit 26 performs setup processing as necessary, and the colorimetric value signal 22 is converted into a dye density signal (cmy) 34 by the dye density converter 36.
[0172]
In this case, the dye density signal (cmy) 34 supplied to one input port of the switch 45 constituting the cmyk converter 32 through the dye density converter 36 is a dye density signal on the reversal document. The dye density signal (cmy) 28 supplied to the input port through the dye density converter 30 is also the dye density signal on the reversal document as described above.
[0173]
The standard condition reproduction table 44 includes standard default separation conditions such as tone correction processing conditions, color correction conditions, under color removal conditions, and K plate generation processing conditions for decomposing the dye density signal on the reversal document. The dye density signals 28 and 34 on the reversal original are each decomposed and converted into a cmyk signal which is a halftone signal under the default decomposition conditions.
[0174]
Therefore, according to the present embodiment, the switch 45 constituting the image signal processing apparatus 10 corresponds to the input corresponding to the output (RGB signal 12) of the digital camera or the output (RGB signal 14) from the input unit of the scanner. By switching to the port side, the effect is achieved that any output can be easily converted to a cmyk signal, which is a half-tone signal, using a common standard default decomposition condition as a result.
[0175]
Note that the present invention is not limited to the above-described embodiment, and various configurations can be adopted without departing from the gist of the present invention.
[0176]
【The invention's effect】
As described above, according to the present invention, after the output RGB signal of the digital camera is converted into a colorimetric value signal, the colorimetric value signal is converted into a dye density signal. Further, after the output RGB signal read by the image reading apparatus is converted into a colorimetric value signal, the colorimetric value signal is converted into a dye density signal. The dye density signal is processed under standard image processing conditions including gradation conversion processing, color correction processing, undercolor removal processing, and K plane generation processing, and converted to a cmyk signal that is a halftone signal.
[0177]
For this reason, regardless of whether the input RGB signal is for a digital camera or an image reading device, it is standardly installed in a specific image reading device that has been optimized empirically, for example, a scanner. The effect that the desired image can be output using the so-called default decomposition condition as it is (stored) is achieved. Further, since the RGB signal is converted into a device-independent colorimetric value signal, an image independent of the device can be output.
[Brief description of the drawings]
FIG. 1 is a block diagram showing an overall configuration of an image signal processing apparatus to which an embodiment of the present invention is applied.
FIG. 2 is a flowchart for explaining a colorimetric conversion matrix creation process;
FIG. 3 is a flowchart for explaining a colorimetric setup unit.
FIG. 4 is a chromaticity diagram used for explaining a color reproduction range of a reversal film on a CIE chromaticity diagram.
FIG. 5 is a flowchart for explaining a conversion process for converting a colorimetric signal of a digital camera into a dye density signal on a reversal document.
FIG. 6 is a flowchart for explaining another conversion process for converting a colorimetric signal of a digital camera into a dye density signal on a reversal document.
FIG. 7 is a flowchart for explaining a dye density conversion matrix creation process;
FIG. 8 is a diagram for explaining the Newton-Raphson method.
FIG. 9 is a flowchart for explaining a part of a dye density conversion matrix creation process;
FIG. 10 is a flowchart for explaining an original scene faithful reproduction table creation process;
FIG. 11 is a flowchart for explaining a process for creating a standard condition reproduction table.
[Explanation of symbols]
DESCRIPTION OF SYMBOLS 10 ... Image signal processing apparatus 12, 14 ... RGB signal
16 ... Colorimetric conversion unit 18, 22 ... Colorimetric value signal (XYZ)
20 ... Colorimetric conversion matrix 24 ... Colorimetric conversion table
25 ... Light source changing unit 26 ... Colorimetric setup unit
28, 34 ... Dye concentration signal (cmy) 30, 36 ... Dye concentration conversion section
32 ... cmyk conversion unit 35 ... image output unit
42 ... Original scene faithful reproduction table 44 ... Standard condition reproduction table
46, 48 ... cmyk signal 64 ... END / END converter

Claims (4)

In a color conversion method comprising a process of converting a device-dependent image signal obtained for an original scene into a colorimetric value signal,
Converting the colorimetric signal into a dye density signal on a reversal document;
The dye density signal is processed under standard image processing conditions that are default separation conditions of a specific color scanner including gradation conversion processing, color correction processing, undercolor removal processing, and K plate generation processing, and converted into a cmyk signal. Process,
A color conversion method characterized by comprising: In a color conversion method including a step of converting a device-dependent image signal obtained for an original scene into a first colorimetric value signal,
Converting the first colorimetric value signal into a first dye density signal on a reversal document;
A process of converting an image signal of a reversal document read by the image reading device into a second colorimetric value signal;
Converting the second colorimetric value signal into a second dye density signal;
The first or second gradation conversion processing dye density signals, color correction processing, and treated by the standard image processing condition of the default cracking conditions for a particular color scanner including undercolor removal process and the K-plate generating process , The process of converting to a cmyk signal;
A color conversion method characterized by comprising: The color conversion method according to claim 1, further comprising a setup processing step of performing a setup processing immediately before the step of converting the colorimetric value signal into the dye density signal. The method according to any one of claims 1 to 3, wherein
The dye density signals gradation conversion processing, color correction processing, and treated by the standard image processing condition of the default cracking conditions for a particular color scanner including undercolor removal process and the K-plate generating process, into a cmyk signal In the process
A color patch on a reversal document in which colorimetric values or signals that can be converted into colorimetric values are regularly arranged is processed under standard image processing conditions that are default separation conditions of the specific color scanner, and a cmyk image is processed. For the patch created and reproduced on the cmyk image, the correspondence between the colorimetric value or the colorimetric value and a convertible signal is picked up as a data set, and the relationship between the two is tabulated and the specified A color conversion method characterized in that a standard image processing condition table which is a default separation condition of the color scanner is prepared, and the dye density signal is converted into a cmyk signal using the image processing condition table and an interpolation method.

Images