JP3671097B2 - 画像信号処理方法 - Google Patents
画像信号処理方法 Download PDFInfo
- Publication number
- JP3671097B2 JP3671097B2 JP31974996A JP31974996A JP3671097B2 JP 3671097 B2 JP3671097 B2 JP 3671097B2 JP 31974996 A JP31974996 A JP 31974996A JP 31974996 A JP31974996 A JP 31974996A JP 3671097 B2 JP3671097 B2 JP 3671097B2
- Authority
- JP
- Japan
- Prior art keywords
- signal
- colorimetric
- color
- rgb
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Editing Of Facsimile Originals (AREA)
- Facsimile Image Signal Circuits (AREA)
- Color Image Communication Systems (AREA)
Description
【発明の属する技術分野】
この発明は、グレー条件が規定されている信号であるEND信号に基づき行われるセットアップ処理をデバイス非依存の系へ適用することを可能とする画像信号処理方法に関する。
【0002】
【従来の技術】
例えば、線走査読取装置であるカラースキャナ(以下、単にスキャナともいう。)は、その入力部において、原稿としてのポジ画像が担持されたカラーリバーサル原稿(カラーリバーサルフイルムが露光され現像されて画像が担持された原稿であって、以下、単にリバーサル原稿ともいう。)に照明光を照射することにより、前記ポジ画像情報を含む透過光を得、この透過光を結像光学系と3色分解光学系を介してCCDリニアイメージセンサ等のリニアイメージセンサに導き、リニアイメージセンサの出力信号である光電変換信号をAD変換することで、前記カラーリバーサル原稿上の画像を光電的に読み取るように構成されている。
【0003】
この場合、例えば、RGB用の各リニアイメージセンサにより前記リバーサル原稿を主走査方向に線走査(ライン走査)して読み取るとともに、前記リバーサル原稿を前記主走査方向と略直交する副走査方向に相対的に搬送することでRGB各色についての2次元的な画像情報であるデジタルRGB画像信号(単にRGB信号ともいう。)を得ることができる。
【0004】
そして、このようにして得られたRGB画像信号は、c(シアン)m(マゼンタ)y(黄)のEND信号に変換され、このEND信号に基づき、ハイライト点とシャドー点がセットアップ処理により決定され、このセットアップポイント(ハイライト点とシャドー点)を基準に網%信号に変換される。
【0005】
網%信号に変換されたcmy信号は、印刷会社のノウハウになっている3色4色変換処理により、cmyk画像信号(単にcmyk信号ともいう。)に変換される。
【0006】
cmyk信号は、例えば、外部の網点化処理装置に供給され、この網点化処理装置において、cmyk各信号が、それぞれ、閾値テンプレート中の閾値と比較され、いわゆる網点化処理(2値化処理)が行われる。
【0007】
そして、この網点化処理後の2値画像信号に基づいて、フイルム作成装置において製版用フイルムが作成され、この製版用フイルムから刷版作成装置により刷版が作成される。作成された刷版(cmyk各版)がカラー印刷機に装着され、インキが付けられて、本紙上に対して多色刷りが行われる。これにより、本紙上に転移されたインキによる印刷画像が形成された印刷物が完成する。
【0008】
なお、最近は、網点化処理後の2値画像信号に基づいて、直接、刷版を作成するCTP(Computer To Plate )装置も実用化されている。
【0009】
【発明が解決しようとする課題】
ところで、RGB信号を網%信号であるcmyk信号に変換する場合には、上記セットアップ処理が必須であるが、従来のセットアップ処理は、いわゆるデバイス依存(デバイスディペンデント)な信号系での処理であるため、そのままでは、デバイス非依存(デバイスインディペンデント)な信号系で使用することはできない。
【0010】
この発明は、このような課題を考慮してなされたものであって、デバイス依存の信号系においてグレー条件が規定されている信号であるEND信号に基づき行われるセットアップ処理を、デバイス非依存の系へ適用することを可能とする画像信号処理方法に関する。
【0011】
【課題を解決するための手段】
第1のこの発明は、測色値信号をRGB3原色の表色系の信号であるRGB信号に変換した後、求めたRGB信号の各々の基礎刺激に対する比率を対数変換して測色的END信号に変換する表色系第1変換過程と、
前記測色的END信号から階調変換カーブを設定するセットアップ過程と、
この設定に基づき、前記測色的END信号を階調変換処理する階調変換過程と、
階調変換処理後の測色的END信号を逆対数変換してRGB信号にもどし、このRGB信号を測色値信号に変換する表色系第2変換過程と、
を有することを特徴とする。
【0012】
第1のこの発明によれば、測色値信号を一旦測色的END信号に変換し、この測色的END信号によりセットアップ処理を行い、階調変換カーブを設定する。次に、この階調変換カーブに基づき、測色的END信号を階調変換し、階調変換された測色的END信号を求める。最後に、この階調変換された測色的END信号を測色値に戻すようにしている。このため、デバイス非依存な信号である測色値信号に対して、結果としてセットアップおよび階調変換処理を行うことができる。
【0015】
第2のこの発明は、第1のこの発明において、前記表色系第1変換過程におけるRGB3原色の表色系は、リバーサルフイルムの色再現域を包含する原刺激を選択した領域とすることを特徴とする。
【0016】
第2のこの発明によれば、リバーサルフイルムについての既存のスキャナ(画像読取装置)のセットアップ処理を利用することができる。
【0017】
第3のこの発明は、非リバーサル原稿より得られたデバイス依存の画像信号を測色値信号に変換する過程と、
前記測色値信号をRGB3原色の表色系の信号であるRGB信号に変換した後、求めたRGB信号の各々の基礎刺激に対する比率を対数変換して測色的END信号に変換する測色的END信号変換過程と、
前記測色的END信号から階調変換カーブを設定するセットアップ過程と、
この設定に基づき、前記測色的END信号を階調変換処理する階調変換過程と、
階調変換処理後の測色的END信号を逆対数変換してRGB信号にもどし、このRGB信号を測色値信号に変換する過程とを有し、
前記測色的END信号変換過程におけるRGB3原色の表色系は、リバーサルフイルムの色再現域を包含する原刺激を選択した領域または前記非リバーサルフイルム原稿の色再現域を包含する原刺激を選択した領域とすることを特徴とする。
【0018】
第3のこの発明によれば、非リバーサル原稿より得られたデバイス依存の画像信号を一旦測色値信号に変換し、その測色値信号を測色的END信号に変換することでセットアップ処理を行い、階調変換カーブを設定する。次に、この階調変換カーブに基づき、測色的END信号を階調変換し、階調変換された測色的END信号を求める。最後に、この階調変換された測色的END信号を測色値に戻すようにしているので、非リバーサル原稿より得られたデバイス依存の画像信号に対して、リバーサル原稿を取り扱う系でのセットアップおよび階調変換処理を行うことができる。
第4のこの発明によれば、結果として、測色値信号をセットアップし、階調変換したRGB信号を得ることができる。
【0019】
【発明の実施の形態】
以下、この発明の一実施の形態について図面を参照して説明する。
【0020】
図1は、この発明の一実施の形態が適用された画像信号処理装置10の全体的な構成を示している。
【0021】
この画像信号処理装置10は、測色変換部16を有し、この測色変換部16には、撮像装置であるデジタルカメラ(不図示)で撮影したシーン(場面)の画像情報を有するRGB信号(デジタル画像信号)12が供給されるとともに、必要に応じて、画像読取装置であるスキャナ(不図示)を構成する画像入力部(線走査読取装置またはスキャナ入力部ともいう。)で読み取られたリバーサル原稿の画像情報を有するRGB信号(デジタル画像信号)14が供給される。なお、デジタルカメラはデジタルビデオカメラでもよく、その他、撮像機能を有したデジタル画像信号(デジタル画像データ)を出力する装置であれば、この発明を適用することができる。
【0022】
また、スキャナとしては、リニアイメージセンサを搭載するカラースキャナである線走査読取装置でもよく、また、エリアイメージセンサを搭載するカラースキャナである面走査読取装置でも、この発明を適用することができる。
【0023】
測色変換部16は、デジタルカメラのRGB信号12をXYZまたはL* a* b* 等の測色値信号(XYZ信号ともいう。)18に変換する測色変換マトリクス(以下、単にマトリクスともいう。)20と、スキャナ入力部のRGB信号14をXYZ信号またはL* a* b* 等の測色値信号22に変換する測色変換テーブル(測色変換ルックアップテーブル、ルックアップテーブルまたは単にテーブルともいう。)24とを有している。
【0024】
なお、通常、RGB信号12、14は、デバイスに依存する(デバイス依存またはデバイスディペンデントともいう。)信号(データ)といい、これに対して測色値信号18、22は、デバイスに依存しない(デバイス非依存またはデバイスインディペンデントともいう。)信号(データ)という。
【0025】
デジタルカメラの測色値信号18は、撮影光源と観察光源の違いを吸収する光源変更部25と、ハイライト部濃度とシャドー部濃度とをセットアップする測色的セットアップ部(単に、セットアップ部ともいう。)26と、原シーン(撮影されたもとの場面)の前記測色値信号18をリバーサル原稿上の色素濃度信号28に変換する色素濃度変換部30と、cmyk変換部32を構成する原シーン忠実再現テーブル42に必要に応じて供給される。
【0026】
一方、スキャナ入力部の測色値信号22は、前記測色的セットアップ部26と、標準条件再現テーブル部43を構成し前記測色値信号22を色素濃度信号34に変換する色素濃度変換部36に供給される。
【0027】
色素濃度変換部30、36から出力される色素濃度信号(cmy信号)28、34は、スイッチ45を通じて選択され、標準条件再現テーブル44に供給される。
【0028】
cmyk変換部32は、基本的には、原シーン忠実再現テーブル42と標準条件再現テーブル部43とを有している。原シーン忠実再現テーブル42は、供給されるデジタルカメラの測色値信号18を測色的に保存された網%信号であるcmyk信号46に変換するルックアップテーブルである。
【0029】
標準条件再現テーブル部43は、測色値信号22を色素濃度信号34に変換する色素濃度変換部36と、色素濃度変換部30または色素濃度変換部36のいずれかの出力信号を選択するスイッチ(マルチプレクサ、選択手段)45と、標準条件再現テーブル44とを有している。標準条件再現テーブル44は、スイッチ45により選択された色素濃度信号28、34のいずれかを網%信号であるcmyk信号48に変換する処理を行う。cmyk変換部32は、3色表色系の信号を4色表色系の信号に変換する、いわゆる3色4色変換機能を有する。
【0030】
cmyk変換部32から出力される網%信号であるcmyk信号46、48は、このcmyk信号46、48に基づく画像を出力する画像出力部35に供給される。画像出力部35は、例えば、図示しない2値化変換部と、レーザ露光走査部(イメージセッタ等)と、現像部と、刷版作成部と、印刷部とからなる公知の構成を採用することができる。2値化変換部では、cmyk信号46またはcmyk信号48をスクリーン線数、網角度等の出力条件に応じて選択されるCMYKの各閾値マトリクスと比較して、2値化処理する。レーザ露光走査部では、この2値信号(2値画像信号ともいう。)に基づいてオンオフするレーザビームによりフイルムを露光走査して潜像を形成する。現像部では、この潜像が形成されたフイルムを現像して画像を顕像化し、製版用フイルムを作成する。刷版作成部では、この製版用フイルムから刷版を作成する。印刷部では、前記刷版、この場合、4版分のCMYK用刷版を印刷機に装着し、刷版に付けられた4色のインキが本紙(印刷用紙)に転写されることで、画像が形成されたハードコピーとしての印刷物が作成される。
【0031】
なお、画像出力部35として、フイルムの現像処理が不要であり、網点やスクリーン線数やスクリーン角度を直接本紙に網点画像として印刷してシミュレーションすることのできる、ダイレクトデジタルカラー校正(DDCP)システムを用いることもできる。
【0032】
前記測色的セットアップ部26は、測色値信号18、22を測色的END(等価中性濃度)信号52に変換する測色的END順変換・逆変換マトリクス(単に、マトリクスともいう。)50と、測色的END信号52から間引き処理により作成した測色的END信号56を出力するラフデータ作成部54と、作成された測色的END信号56に基づいて、ハイライト濃度信号60とシャドー濃度信号62とを自動的に決定するオートセットアップ部58と、測色的END信号52に対して階調変換処理を施して測色的END信号66に変換するEND・END変換部64と、変換された測色的END信号66を測色値信号22に逆変換する測色的END順変換・逆変換マトリクス50とを有する。
【0033】
なお、画像信号処理装置10は、図示していないコンピュータ(CPU、ROM、RAM、外部記憶装置、モニタ、その他、入出力機器等を含む。)により制御される構成となっており、画像信号処理装置10を構成する各ブロックは、ハードウエアばかりでなく、ソフトウエアで構成される部分をも有する。コンピュータは、制御、判断、演算、比較手段等として機能する。
【0034】
次に、画像信号処理装置10を構成する各ブロックの詳細な構成と動作について説明する。
【0035】
測色変換部16を構成し、RGB信号12をXYZ信号18に変換するマトリクス20は、図2に示すフローチャートに基づいて作成される。なお、以下の説明において、マトリクス20は、RGB色空間からCIE−L* a* b* (光源:補助標準の光CIE−D50)色空間への変換を例として説明する。この場合、CIE−L* a* b* 色空間とXYZ色空間との相互の変換は、次に示す公知の(1)式により一意に行うことが可能である。したがって、以下の全ての説明において、XYZ色空間(またはCIE−L* a* b* 色空間)での処理はCIE−L* a* b* 色空間(またはXYZ色空間)での処理に置き換えることができる。また、その他、これらと等価な測色的色空間での処理に置き換えることができる。
【0036】
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)
まず、代表的な色を彩度、明度、色相に応じて変化させた複数のカラーパッチ72(図2参照)を有する一種の色票であるカラーチャート70を準備する(ステップS1)。この実施の形態では、カラーチャートとして、マクベスカラーチェッカー{登録商標:米国のコールモージェン社のマクベス部門(Macbeth A division Kollmorgen )製}を用いる。マクベスカラーチェッカーは、公知のように、CIE(1931)xyY値、色相、マンセル表記値、彩度が規定されたカラーチャートである。
【0037】
24色は、具体的には、
1.暗い肌(dark skin )
2.明るい肌(light skin)
3.青い空(blue sky)
4.葉(foliage )
5.青い花(blue flower )
6.青みの緑(bluish green)
7.オレンジ(orange)
8.紫みの青(purplish green)
9.明度彩度のほどよい赤(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.中性8(neutral 8:明るい灰色で、8はマンセル表記値の8)
21.中性6.5(neutral 6.5:ライトミディアム灰色)
22.中性5(neutral 5:中間の灰色)
23.中性3.5(neutral 3.5:暗い灰色)
24.黒(black )
である。
【0038】
なお、カラーチャートとしては、マクベスカラーチェッカーに限らず、例えば、JIS標準色標等、色空間を概ね均等に網羅したカラーチャートを用いることができる。
【0039】
次に、CIE−D50の撮影光源下でデジタルカメラを用いてカラーチャート70の24色の各色、すなわち24個の各パッチ72を撮影して、パッチ72のそれぞれのRGB信号12を得、得られた各RGB信号12を輝度値に変換した後、1/3乗する(ステップS2)。なお、輝度値への変換は、例えば、デバイス内部で施されているγ補正を解除することで求めることができる。また、1/3乗するのは、上記(1)式からも理解されるように、得られたRGB信号12をCIE−L* a* b* 表色系で処理するためである。
【0040】
次いで、ステップS2で得られた各パッチ72についてのRGB信号12、具体的には、各パッチ72について、輝度変換され1/3乗されたR値、G値、B値から2次項までの各値R、G、B、R2 、G2 、B2 、RG、GB、BR(9変数)を計算する(ステップS3)。
【0041】
次に、後に説明する重回帰分析を行った場合に、マルチコ(多重共線性)現象を起こさないように、ステップS3で得られた9変数の24色のデータに対して主成分分析を行い、9変数の主成分スコア(主成分得点)Vを求める。各色毎に次の(2)式で示す主成分スコアV(Vは、ベクトルと考える。)が求められる。
【0042】
V=(v1 ,v2 ,v3 ,v4 ,v5 ,v6 ,v7 ,v8 ,v9 )…(2)
なお、この(2)式において、主成分スコアVの成分v1 ,v2 ,v3 ,v4 ,v5 ,v6 ,v7 ,v8 ,v9 は、相互に相関が全くなく前記マルチコ現象を起こさない。
【0043】
次に、図示していない測色計により、カラーチャート70の各パッチ72、すなわち24色の各色についての測色値C(L* a* b* )(Cも、24色分得られるので、ベクトルと考える。)を求める(ステップS5)。この測色値Cを求める過程は、ステップS2〜S4の過程に対してどの時点で求めてもよい。
【0044】
次に、測色値Cを目的変数(従属変数)とし、主成分スコアVを説明変数(独立変数)として、重回帰分析により偏回帰係数A(Aもベクトルと考える。)を求める(ステップS6)。
【0045】
この重回帰分析を行う際に、目的変数群となっている測色値Cを構成する24色に1:1に対応した重み付けマトリクス(単に、マトリクスともいう。)P=[Pi](i=1,2,…24)を作用させる(ステップS7)。Piは、各色の重みであり、例えば、上述の24色中、20.中性8(neutral 8:明るい灰色で、8はマンセル表記値の8)、21.中性6.5(neutral 6.5:ライトミディアム灰色)、22.中性5(neutral 5:中間の灰色)、23.中性3.5(neutral 3.5:暗い灰色)の灰色についての重みを大きくすることで、灰色についての色の再現性を向上させることができる。灰色に限らず、目的に応じて所望の色、例えば、人であれば肌色、外の景色であれば空色について限定して重みを大きくすれば、肌色または空色をより忠実に再現することが可能となり、好適である。また、重みは、入力RGB信号12を統計的に処理し、一画面分の画像信号、言い換えれば、画素分割された画像データについて、例えば、カラーチャートの各パッチの色を中心とした領域に色空間を分割し、各領域に存在する画素の頻度値に応じた重み付け、例えば、最頻値の重みを最も大きくし、以下、順に重みを小さくするようにすれば自動的な重み付け処理を行うことができる。なお、色空間の分割は、RGB色空間であってもよく、重みなしのマトリクスで変換されたXYZ(L* a* b* )色空間であってもよい。
【0046】
この場合、重みを大きくした色を画像の主要色と呼ぶ。なお、マトリクスPの各要素の合計値が1になるように、各要素の値を各要素の合計値で割った値に規格化しておく。
【0047】
偏回帰係数Aを求めるステップ6の重回帰分析について詳しく説明する。
【0048】
目的変数である測色値(ベクトル)Cと、求めようとする偏回帰係数(ベクトル)Aと、主成分スコア(ベクトル)Vとの間で、24色の各測色値Cに対して次の(3)式に示す線形1次式が成立するものとする。
【0049】
C=AV …(3)
この(3)式と等価なマトリクスによる表現式を(4)式に、和の記号Σによる表現式を(5)式に各々示す。
【0050】
【数1】
【0051】
【数2】
【0052】
なお、(5)式のように和の記号Σを使用する式を、必要に応じて、L* =Σ(j=1→10)a1jvj と表すこととする。
【0053】
偏回帰係数Aは、L* a* b* のそれぞれ別個に最小二乗法を用いて求める。例えば、L* については、次の(6)式によるeL を最小とする偏回帰係数a1jを求めればよい。
【0054】
eL =Σ(i=1→24)Pi{Ci−Σ(j=1→10)a1jvj }2
…(6)
ここで、iは、カラーチャート70のパッチ番号、Piは、各色の重み、jは、変数の番号(1、2、…、10)である。
【0055】
(6)式をベクトルとマトリクスで表すと(7)式となる。ただし、(7)式において、測色値Cと偏回帰係数aとはベクトル、主成分スコア[V]と重み付けマトリクス[P]とはマトリクスである。tはマトリクスの転置を表す。
【0056】
eL =
(ベクトルC−ベクトルa[V])t [P](ベクトルC−ベクトルa[V])…(7)
以下、ベクトルCは単にCと記載し、ベクトルaは単にaと記載する。先ず、(7)式は、以下のように変形することができる。
【0057】
eL =(Ct −[V]t at )t [P](Ct −[V]t a)
ここで、aと[V]とを入れ換えれば、
eL =(C−a[V])[P](Ct −[V]t a)となる。
【0058】
一般的に、(ABC)t =Ct Bt At であるので、
eL =C[P]Ct +a[V][P][V]t at
−a[V][P]Ct −C[P][V]t at となる。
【0059】
ここで、[V][P][V]t =[N]、[V][P]Ct =[U]と置けば、
eL =C[P]Ct +a[N]at −2a[U] …(8)
となる。(8)式において、eL を最小にするためには、偏回帰係数aの各要素についての微分が0に等しくならなければならない。従って、次の(9)式が成立する。
【0060】
【数3】
【0061】
この(9)式から、以下に示すように偏回帰係数aを(10)式により求めることができる。
【0062】
([ ]-1は逆行列を示す。)
このようにしてL* についての偏回帰係数a(a1j)が10個求められ、以下、同様の方法で、残りのa* 、b* についても偏回帰係数a2j、a3jを各10個求めることができる。求められた偏回帰係数a1j、a2j、a3jの合計は、3×10個である{(4)式参照}。
【0063】
結局、L* a* b* についての偏回帰係数Aをまとめて表現すると、(11)式で表すことができる。
【0064】
A=CPVt (VPVt )-1 …(11)
次に、各目的変数L* a* b* について求めた偏回帰係数Aの例えば5%有意で検定を行い(ステップS8)、5%有意(信頼度95%)な説明変数Vを記憶手段に格納しておくとともに、偏回帰係数Aを図1に示したマトリクス20として格納する(ステップS9)。なお、1%有意(信頼度99%)で検定してもよい。
【0065】
ステップS8の検定では、(12)式〜(16)式に示すように、回帰平方和SR {求めた推定値L* i(Lの上に記号「^」を付けている。)と測色値の平均L* i(Lの上に記号「- 」を付けている。)の差の2乗和:(12)式参照}と残差平方和SE {測色値L* iと推定値L* iの差の2乗和:(13)式参照}とを求め、さらに回帰平方和SR と残差平方和SE の不偏分散VR 、VE {(14)式、(15)式参照}を各々求める。さらに、偏F値Fj{(16)式参照}を求める。なお、(16)式において、aの頭に「^」を付けたものは、(4)式の右辺のマトリクスの係数の推定値、Sjjは、説明変数vj の分散共分散行列の逆行列の対角項を意味している。
【0066】
【数4】
【0067】
この偏F値Fjが5%有意のF分布を参照して求めた値=F′n-p-1 (0.05)=F′24-9-1(0.05)=F′14(0.05)=4.60011より大きい場合には、5%有意水準で回帰が有意であり、予測に役立つ偏回帰係数と判断してマトリクス20(図1参照)に格納しておく。
【0068】
このようにして、デジタルカメラにより得られたRGB信号12(図1参照)をXYZ測色値信号18に変換するL* a* b* を求めるためのマトリクス20を作成することができる。このマトリクス20に係る式を(17)式に示す。
【0069】
L* =Σa1jvj (j:5%有意な変数)
a* =Σa2jvj (j:5%有意な変数)
b* =Σa3jvj (j:5%有意な変数) …(17)
この(17)式において、(j:5%有意な変数)の意味は、変数として前記検定が5%有意であったもののみを使用するという意味である。その点で上述の(5)式と異なる。なお、a1j、a2j、a3jは、偏回帰係数マトリクスAの要素であり、vj は、ステップS4で求めた説明変数の主成分スコアである。
【0070】
このマトリクス20を使用することにより、既成のカラーチャート70を使用して測色変換を行うことができる。この場合、ルックアップテーブルを使用しないのでメモリ容量が少なくて済み、換言すれば、メモリ容量が少なくても、精度よく測色変換を行うことができる。さらに、重み付けマトリクスPの使用により画像の所望色(主要色)に限定して精度よく色変換を行うこともできる。なお、説明変数の数が9箇の場合には、ステップS8の検定処理を行わず、全ての偏回帰係数Aを用いて(17)式に対応するマトリクスを作成してもよい。
【0071】
このようにして得られたCIE−L* a* b* 色空間の値が、上述の(1)式によりXYZ色空間の値へ変換され、マトリクス20の出力信号である測色値信号18とされる(図1参照)。
【0072】
次に、測色値信号18は、必要に応じて光源変更部25により測色的に変更され、新たな測色値信号18(同一の符号を用いる。)とされる。なお、観察光源をデジタルカメラの撮影時における撮影光源(CIE−D50と同一)とした場合には、この光源変更処理は不要である。
【0073】
測色値XYZを光源変更部25による光源変更処理後の新たな測色値XYZ(符号は、X′Y′Z′とする。)に変更する場合、その変換は、次の(18)式、(19)式、(20)式に基づいて行うことができる。
【0074】
X′=XX2 /X1 …(18)
Y′=YY2 /Y1 …(19)
Z′=ZZ2 /Z1 …(20)
ただし、(18)式〜(20)式において、XYZは各々撮影光源下における測色値、X′Y′Z′は各々観察光源下における測色値、X1 Y1 Z1 は各々撮影光源の白色点、X2 Y2 Z2 は各々観察光源の白色点である。
【0075】
以下、測色値信号18は、必要に応じて光源変更処理のなされた信号であるものとする。
【0076】
一方、スキャナ入力部から出力されたRGB信号14は測色変換テーブル24によりXYZ信号である測色値信号22に変換される。
【0077】
この測色変換テーブル24は、例えば、cmyの色パッチの各濃度を13段階規則正しく振った、合計13×13×13=2197個からなるカラーパッチを有するカラーリバーサル原稿を準備する。そして、このカラーリバーサル原稿を構成する前記各カラーパッチをスキャナの入力部で読み取るとともに、測色計で読み取る。そして、スキャナで読み取ったRGB値と測色計で読み取ったXYZ値との対応関係を求めてテーブルとしたものである。このテーブルに存在しない読取値間の値は補間処理により求める。
【0078】
そして、スキャナの入力部の出力信号であるRGB信号14を測色変換テーブル24により測色値信号22に変換すること、およびデジタルカメラの出力信号であるRGB信号12を測色変換マトリクス20により測色値信号18に変換することにより、次に説明する測色的セットアップ部26を共通に使用して、いわゆるオートセットアップ処理を行うことができる。すなわち、測色値信号18、22を、一旦、測色的ENDに変換することにより、オートセットアップ処理ソフトウエアを共通に使用することができるという利点がある。
【0079】
次に、測色的セットアップ部26について図3のフローチャートをも参照して説明する。
【0080】
従来から、階調変換は、直感的な理解がし易いために濃度空間で行われており、また、現在、市販されているカラースキャナにおける階調変換処理および色修正処理等の画像処理も濃度信号を基準として実施されている。そこで、まず、測色値信号18または測色値信号22を濃度信号の表色系であるcmyのEND(等価中性濃度)信号52に変換する(ステップS11)。この変換処理は、図3に示すように、XYZ表色系からRGB表色系への線形変換処理(ステップS11a)とRGB表色系からcmy表色系への非線形変換処理(ステップS11b)の2段階の処理になっている。なお、測色値信号18または測色値信号22から変換されたEND信号を測色的END信号52と呼ぶ。
【0081】
例えば、リバーサル原稿の測色値信号22を測色的END信号52に変換する場合、図4に示すCIEの色度図上、リバーサルフイルムの色再現領域71(ハッチングで示す領域)を描き、その色再現領域71を包含する領域73の3つの原刺激RGB(それぞれRxyz、Gxyz、Bxyzとする。)の各点の色度図上の座標、すなわち色度座標が(21)式〜(23)式に示す座標であるものとする。この場合、リバーサルフイルムの色再現域71を包含する領域73は、図4に示す色度図上、原刺激Rxyz、Gxyz、Bxyzを頂点とする三角形の領域である。
【0082】
Rxyz=Rxyz(xR ,yR ,zR ) …(21)
Gxyz=Gxyz(xG ,yG ,zG ) …(22)
Bxyz=Bxyz(xB ,yB ,zB ) …(23)
また、色度図上、XYZ表色系の基礎刺激(白色刺激)Wxyzの座標を(24)式に示す。
【0083】
Wxyz=Wxyz(xW ,yW ,zW ) …(24)
この場合、次の(25)式によりXYZ表色系の測色値信号22(右辺の右側のマトリクス)は、変換マトリクス(右辺の左側のマトリクス)を介してRGB表色系の色信号RGB(左辺のマトリクス)に変換することができる(ステップS11a)。
【0084】
【数5】
【0085】
余因数Aijは、下式で求められる。
【0086】
Aij=(−1)i+j Dij
ただし、Dijは、i行j列を除いた小行列式である。
【0087】
測色的END信号52のcmy値は、基礎刺激(白色刺激)WXYZ の座標を(25)式に代入して求めたR、G、Bの値Rw、Gw、Bwに対する(25)式で求めたR、G、Bの値の各比率を(26)式〜(28)式に示すように各々対数変換することにより求めることができる。
【0088】
c=−log(R/Rw) …(26)
m=−log(G/Gw) …(27)
y=−log(B/Bw) …(28)
ラフデータ作成部54は、セットアップ処理を短時間に行うため、例えば、原稿1枚分の測色的END信号(実際には、デジタルデータ)52を対象とするのではなく、原稿中、オペレータからCPUを通じて領域指定された画像の存在する部分のみのデータを選択すること、あるいは画像が原稿の全面に存在する場合にはデータを間引いて作成すること等、いわゆるラフスキャン処理を行う。
【0089】
次に、オートセットアップ部58において、ラフデータ作成部54により選択されたラフデータである測色的END信号56にもとづいてオートセットアップ処理を行う(ステップS12)。このオートセットアップ処理では、例えば、特開平2−105676号公報にも公知のように、測色的END信号56についてヒストグラムを作成した後、累積ヒストグラムを作成する。
【0090】
そして、図1中、END・END変換部64の動作を示す特性図(図1中に描いている図)に示すように、前記累積ヒストグラムの、例えば、0.1%点データ(HL濃度)D1に対応する測色的END信号56の値をDHに設定し、98%点データ(SD濃度)D2に対応する測色的END信号56の値をDSに設定する(ステップS12)。
【0091】
分かり易く仮想的な値で説明すると、例えば、D1の濃度がD1=1.0、D2の濃度がD2=2.0であった場合、D1の濃度に対応するDHをDH=0.1に設定し、D2の濃度に対応するDSをDS=3.0に設定する。なお、実際上は、0.1%点データ(HL濃度)D1が網%の0%に対応する濃度に変換され、98%点データ(SD濃度)D2が網%の100%に対応する濃度に変換される。
このようにして設定された直線74の式から、END・END変換部64により本スキャンデータ(ラフスキャン処理の対象となったデータ)の全ての測色的END信号52を測色的END信号66に変換することができる(ステップS13)。すなわち、階調特性を変換すること、言い換えれば階調補正処理を行うことができる。なお、変換式は、ハイライトポイント75(D1,DH)およびシャドーポイント76(D2,DS)を通る曲線の式とすることもできる。なお、HL濃度D1とSD濃度D2(またはセットアップポイント75、76)を上述のように予め定めた一定の条件に基づいて自動的に決定することをオートセットアップというが、HL濃度D1とSD濃度D2の値をマニュアルで決定することも可能であり、あるいは、オートセットアップ後、DHとDSとをマニュアルで修正することも可能である。マニュアルにより決定する動作をマニュアルセットアップという。
【0092】
次いで、測色的END逆変換マトリクス50により、END・END変換後の測色的END信号66のcmy値をXYZ値である測色値信号22に逆変換する(ステップS14)。
【0093】
この逆変換処理は、測色的END信号66のcmy値をRGB表色系のRGB値に変換する処理(ステップS14a)と、このRGB値をXYZ表色系のXYZ値をとる測色値信号22に変換する処理(ステップS14b)とからなる。測色的END信号66のcmy値をRGB表色系のRGB値に変換する処理は、前記の(26)式〜(28)式をRGBについて解いた次の(29)式〜(31)式により行われ、RGB値をXYZ表色系のXYZ値に変換する処理は、前記の(25)式をマトリクスXYZについて解いた次の(32)式により行われる。
【0094】
R=Rw10-c …(29)
G=Gw10-m …(30)
B=Bw10-y …(31)
【0095】
【数6】
【0096】
なお、上述のセットアップ処理は、リバーサル原稿の測色値信号22を例に説明したが、同様にデジタルカメラの測色値信号18に対しても適用することができる。
【0097】
すなわち、一般的に、印刷製版分野においては、オリジナル原稿としてリバーサル原稿を使用する。このリバーサル原稿上に原シーンが露光され、そのリバーサル原稿上で発色している色再現域71内を対象として印刷は仕上げられる。このことは、取り扱う画像信号の色再現域として、リバーサル原稿が発色する色再現域71で十分であることを示唆している。
【0098】
このように、測色的セットアップ部26では、測色的XYZをRGB3原色の表色系へ線形変換し、基礎刺激(光源)Rw、Gw、Bwに対する求めたR、G、Bの比率を対数変換することにより、END値を求め、これを用いてオートセットアップする。END・END変更されたEND値を逆対数変換処理することにより、変更されたRGB値を求め、逆変換マトリクスにより変更された測色値XYZを求めるようにしている。このため、既存のいわゆる濃度信号に基づいたオートセットアップ処理を行うことができるという効果が達成される。なお、セットアップ処理は、マニュアルセットアップ処理として行うことができることはいうまでもない。
【0099】
次に、色素濃度変換部30について図5、図6のフローチャートをも参照して説明する。
【0100】
この色素濃度変換部30では、デジタルカメラにより撮影された原シーンの測色値信号(XYZ)18をリバーサル原稿上の色素濃度信号(cmy)28に変換する。
【0101】
この変換過程は、図5に示す方法と図6に示す方法の2通りの方法が考えられる。
【0102】
すなわち、図5に示す方法は、原シーンの測色値信号(XYZ)18を、後述するルックアップテーブルにより原シーンの色素濃度信号(cmy)81に変換し(ステップS21:第1段階の処理)、この原シーンの色素濃度信号(cmy)81を後述するマトリクスによりリバーサル(単に、RVともいう。)原稿上の色素濃度信号(cmy)28に変換する方法である(ステップS22:第2段階の処理)。また、図6に示す方法は、原シーンの測色値信号(XYZ)18を、後述するマトリクスによりリバーサル原稿上の測色値信号(XYZ)82に変換し(ステップS31:第1段階の処理)、このリバーサル原稿上の測色値信号(XYZ)82を後述するルックアップテーブルによりリバーサル原稿上の色素濃度信号(cmy)28に変換する(ステップS32:第2段階の処理)方法である。
【0103】
まず、図5に示す方法の詳細について、図7のフローチャートをも参照して説明する。図5中、ステップS21の処理に係る測色値XYZを色素濃度cmyに変換するためには、色空間を概ね均等に網羅したカラーチャート、例えば、上述のマクベスカラーチェッカー70の24色の各測色値XYZを測色計により求める(ステップS21a:図7参照)。
【0104】
次に、求めた測色値XYZから色素濃度値cmyを算出し(ステップS21b)、測色値XYZに対する色素濃度値cmyのルックアップテーブルを作成する。このルックアップテーブルを色素濃度変換部30における第1段階の処理であるステップS21のルックアップテーブルとして色素濃度変換部30に設定記憶しておく。なお、ステップS21の処理では、このルックアップテーブルを用いるとともに、上記24色の測色値XYZ以外の色については、その24色の測色値XYZを用いて補間法により色素濃度値cmyを算出することができる。
【0105】
ステップS21bの算出処理、すなわち測色値(三刺激値)XYZから色素濃度cmyを求める処理を説明する。
【0106】
測色値XYZと色素濃度cmyとの間では、次の(33)式〜(36)式が成立する。
【0107】
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)
ただし、
h={cDc(λ)+mDm(λ)+yDy(λ)+base(λ)}
k=100/∫visP(λ)y(λ)dλ(λは光の波長)
∫vis:可視波長域(380nm〜780nm)での定積分
P(λ):観察光源の分光特性データ
T(λ):透過物体の色素の分光透過率データ
{(36)式はランベルト・ベールの法則が成立すると仮定}
x(λ)、y(λ)、z(λ):等色関数
Dc(λ)、Dm(λ)、Dy(λ):cmy色素の分光濃度データ
base(λ):フイルムベースの分光濃度データ
これら(33)式〜(36)式から色素濃度cmyを求める場合には、逆関数を求めればよいのであるが、直接求めることができない。そこで、公知のニュートン(ニュートン・ラフソン)法(例えば、太田登著「色彩工学」pp254−260 東京電機大学出版局 1993年12月20日 第1版1刷発行)のような逐次近似法を用いればよい。前記参考書を引用して、このニュートン・ラフソン法(N・R法と略記する。)について簡単に説明する。
【0108】
一般的な方程式y=f(x)をf(x)=0の根に近いx=x0でテーラ展開して1次の項のみを取ると、xの微小変化Δxに対して、f(x0+Δx)=f(x0)+f′(x0)・Δxが成立する。ただし、f′(x0)は、f(x)の微分係数f′(x)にx=x0を代入したものである。したがって、f(x)=0のより正しい値x1は、f(x0+Δx)=0として、x1=x0+Δx=x0−f(x0)/f′(x0)で求められる。これは、図8に示すように、関数y=f(x)上の点(x0,y0)で接線83を引き、その接線83とx軸との交点x1を求めることに相当する。
【0109】
これを(33)式〜(36)式に適用するためには、(36)式を(33)式〜(35)式に代入した後、ある関数fx、fy、fzを用いれば、(33)式〜(35)式は、(37)式〜(39)式のように表すことができる。
【0110】
X=fx(c,m,y) …(37)
Y=fy(c,m,y) …(38)
Z=fz(c,m,y) …(39)
これら(37)式〜(39)式において、初期値をc0,m0,y0とし、そのときの三刺激値をX0,Y0,Z0とする。今、c0,m0,y0に微小変化Δc,Δm,Δyを加えたとき、三刺激値がΔX,ΔY,ΔZだけ変化したとすれば、次の(40)式が得られる。
【0111】
ただし、例えば、∂fx/∂cは、関数fxのcに関する偏微分係数を表す。
【0112】
この(40)式を整理すれば、(41)式が得られる。同様に、ΔY、ΔZも(42)式、(43)式に示すように得られる。
【0113】
Δ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)
(41)式〜(43)式を(44)式に示すように行列表示する。
【0114】
(Q)=(J)(P) …(44)
ただし、(Q)は、要素が1行目から順にΔX,ΔY,ΔZである1行3列のマトリクス、(J)は偏微分係数が3行3列からなるヤコビアンマトリクス、(P)は、要素が1行目から順にΔc,Δm,Δyである1行3列のマトリクスである。
【0115】
(44)式の両辺に、ヤコビアンマトリクス(J)の逆マトリクス(J)-1を掛ければ、(45)式が得られる。
【0116】
(P)=(J)-1(Q) …(45)
したがって、初期値c0,m0,y0を(46)式のようにそれぞれc1,m1,y1に修正すれば、より正しい近似値を得ることができる。
【0117】
c1=c0+Δc
m1=m0+Δm
y1=y0+Δy …(46)
以上のようにして得られるヤコビアンマトリクス(J)を用いて繰り返し計算を行うことにより、任意の目標値である測色値XYZに対する色素濃度信号cmyを求めることができる。同様の処理をXYZ色空間における格子上の全目標値に対して行うことで、原シーンに対する測色値信号(XYZ)18を原シーンに対する色素濃度信号(cmy)81に変換するための逆変換テーブルを生成し、これをルックアップテーブルとして色素濃度変換部30に第1段階の処理用として保持しておく(ステップS21)。
【0118】
次に、この第1段階の処理に続いて行われるステップS22の第2段階処理用のマトリクスの作成処理について説明する。なお、このマトリクスの作成処理は、図2を参照して既に説明した重回帰分析処理を含む処理であり、簡潔に説明する。
【0119】
まず、ステップS21bで算出した色素濃度cmyについて、c、m、y、c2 、m2 、y2 、cm、my、ycの9変数24色のデータを算出する(ステップS22a)。
【0120】
次いで、この9変数24色のデータを主成分分析し、9主成分の主成分スコアを求める(ステップS22b)。
【0121】
一方、リバーサルフイルムに露光されたマクベスカラーチェッカー24色についての測色値を測色計により求める(ステップS22c)。
【0122】
次に、上述したN・R法によりリバーサルフイルム上の色素濃度cmy(RVcmyとも記載する。)を求める(ステップS22d)。
【0123】
次に、リバーサルフイルム上の色素濃度RVcmyの各々(RVc、RVm、RVy)を目的変数とし、ステップS22bで求めた主成分スコア(定数項を含む)を説明変数として、重回帰分析処理により3×10の偏回帰係数マトリクスを求める(ステップS22e)。
【0124】
なお、重回帰分析処理を行う際に、上述したように、目的変数群となっている24色に1:1に対応した重み付けマトリクスを作用させてもよい(ステップS22f)。
【0125】
次に、各目的変数RVc、RVm、RVyについて求めた偏回帰係数について、例えば5%有意で検定を行い(ステップS22g)、5%有意な説明変数を記憶手段に格納しておくとともに、偏回帰係数を第2段階処理用のマトリクスとして、色素濃度変換部30に格納しておく(ステップS22h)。なお、この場合においても、検定を行わず、全ての係数を用いてもよい。
【0126】
このようにして、原シーンの測色値信号(XYZ)18をリバーサルフイルム上の色素濃度信号cmyに変換する色素濃度変換部30が構築される。
【0127】
次に、図6に示す方法の詳細について、図9のフローチャートをも参照して説明する。なお、図9に示す処理も、図2を参照して既に説明した重回帰分析処理を含む処理であり、簡潔に説明する。
【0128】
まず、マクベスカラーチェッカー70の24色の各測色値XYZを測色計により求める(ステップS31a:図9参照)。
【0129】
次に、求めた測色値XYZについて、X、Y、Z、X2 、Y2 、Z2 、XY、YZ、ZXの9変数24色のデータを算出する(ステップS31b)。
【0130】
次いで、この9変数24色のデータを主成分分析し、9主成分の主成分スコアを求める(ステップS31c)。
【0131】
一方、リバーサルフイルムに露光されたマクベスカラーチェッカー24色についての測色値XYZを測色計により求める(ステップS31d)。なお、リバーサルフイルムに露光されたマクベスカラーチェッカーとは、カラーリバーサルフイルムが装着された光学式カメラにより、所定の光源下で前記マクベスカラーチェッカーを撮影してリバーサルフイルムを露光した後、この露光されたリバーサルフイルムを現像して得られるリバーサル原稿をいう。
【0132】
次に、リバーサル原稿上の測色値XYZの各々を目的変数(RVX、RVY、RVZ)とし、ステップS31cで求めた主成分スコア(定数項を含む)を説明変数として、重回帰分析処理により3×10の偏回帰係数マトリクスを求める(ステップS31e)。なお、重回帰分析処理を行う際に、上述したように、目的変数群となっている24色に1:1に対応した重み付けマトリクスを作用させてもよい(ステップS31f)。
【0133】
次に、各目的変数RVX、RVY、RVZについて求めた偏回帰係数の例えば5%有意で検定を行い(ステップS31g)、5%有意な説明変数を記憶手段に格納しておくとともに、偏回帰係数を第1段階の処理用マトリクスとして、色素濃度変換部30に格納しておく(ステップS31h、ステップS31)。この場合においても、検定を行わず、全ての係数を用いてもよい。
【0134】
次いで、ステップS32の第2段階の処理用テーブルは、図7を参照して説明したステップS21bあるいはステップS22dの方法と同一の方法で作成されるので、その説明を省略する。
【0135】
このようにして、原シーンの測色値信号(XYZ)18をリバーサル原稿上の色素濃度信号cmyに変換する図6の処理に基づく色素濃度変換部30が構築される。
【0136】
次に、cmyk変換部32を構成する原シーン忠実再現テーブル42の作成処理について、図10に示すフローチャートに基づいて説明する。
【0137】
原シーン忠実再現テーブル42は、測色値信号(XYZ)18を測色的に保存された各網%データであるcmyk信号(cmykデータ)46に変換するルックアップテーブルである。
【0138】
なお、この原シーン忠実再現テーブル42の作成処理は、この出願人の出願に係る特願平8−154584号明細書に記載された作成処理と同じ処理である。
【0139】
この場合、まず、画像出力部35に対して、規則的な間隔からなる複数のcmyk網%データを与えることにより、cmykの濃度および混合割合が段階的に変化するカラーパッチを有するcmykカラーチャートを作成する(ステップS41)。
【0140】
この場合、例えば、cmykの規則的な間隔とは、cmykの各色について網%で0、20、…、100%等、20%きざみで増加させた6段階の間隔とすれば、カラーパッチの総数は、46 =1296になる。
【0141】
次いで、画像出力部35により作成されたcmykカラーチャートの各パッチを測色計を用いて測色し(ステップS42)、この測色データから測色値(刺激値)XYZを求め、前記cmyk網%データから測色値XYZデータへの変換テーブル(順変換テーブルという。)を作成する(ステップS43)。
【0142】
なお、この順変換テーブルは、補間用のテーブルとしても使用するので、前記規則的な間隔は、補間の正確性を考慮した場合には細かいほどよいが、その分、測色作業が膨大となるので、その測色作業の繁雑さと前記正確性と以下に説明するコンピュータの処理時間との、いわゆるトレードオフで決定すればよい。
【0143】
ところで、前記原シーン忠実再現テーブル42において、入力される任意の測測色値信号(測色値XYZデータ、刺激値データXYZ、または単にXYZともいう。)18からこれに対応するcmyk信号(cmyk網%データ、色データcmyk、または単にcmykともいう。)46を求めることは、変数が3変数から4変数と増加するため、1つの測色値XYZデータ18に対してcmyk網%データ46の複数の解が存在する可能性がある。この問題を解決するためには、3変数同士の関係とする必要があり、ここでは、cmyk網%データのうち、色データkについて画像出力部35が採りうる最大値Kmax(k=Kmax)に固定する(ステップS44)。最大値Kmaxとは、例えば、cmyk網%データのkの値が100%である。
【0144】
そして、3変数である任意のXYZの値を3変数である対応するcmy(kは固定)の値に変換する(ステップS45)。
【0145】
この場合、任意のXYZの値である目標値X0、Y0、Z0に対するcmyk(ここでは、kはk=Kmaxに固定しているのでkは定数であり、その意味でcmyの3変数)の値である、c0、m0、y0、k0=Kmaxを求める場合には、k0=Kmaxであるときの上記順変換テーブルを利用して、回帰式の偏回帰係数を求める。
【0146】
このときの回帰式は、係数項となる3行4列のマトリクスをA、測色値XYZの1行3列のマトリクスをT、kを固定した、換言すれば定数と考えて、これを含めた残りのcmyの1行4列のマトリクスをDと表せば、次の(47)式のように表すことができる。
【0147】
T=AD …(47)
この(47)式は、次の(48)式の関係を表している。
【0148】
【数7】
【0149】
(47)式、(48)式において、マトリクスDを構成する「1」は、cmyの3次元平面式における定数項を与えるために設定した値である。
【0150】
この(48)式における係数Aは、k=Kmaxのときに得られている上述の順変換テーブルの各データセットを代入することで、上述した重回帰分析により求めることができる。
【0151】
次に、k=Kmaxのときに求められた回帰式を用いて、上述のN・R法により、目標値X0、Y0、Z0に対応するcmy(kはKmaxとして設定されている。)の値を求めることができる。
【0152】
次に、求められたcmyの値が画像出力部35の再現範囲内の値(色)であるかどうかの判断を行う(ステップS46)。すなわち、網%データ(色データ)cmyの再現可能な最小濃度をCmin、Mmin、Ymin、最大濃度をCmax、Mmax、Ymaxとして求められたデータcmyが次の(49)式の関係を満足するかどうかを判定する。
【0153】
Cmin≦c≦Cmax
Mmin≦m≦Mmax
Ymin≦y≦Ymax …(49)
求めたcmyの各値が、(49)式を満足している場合には、刺激値データXYZの目標値(X0,Y0,Z0)に対する色データcmykを各々c=csol、m=msol、y=ysol、k=ksol(この場合、ksol=Kmax)として逆変換テーブルを作成し、これを原シーン忠実再現テーブル42として設定する(ステップS47)。以下の説明において、必要に応じて(csol,msol,ysol,ksol)のデータセットをcmyksolという。
【0154】
ステップS45の過程で求めたcmyデータの値が(49)式を満足しない値であった場合には、k=Kmaxに固定してある色データkをk=k−Δk、この場合k=Kmax−Δkとした後(ステップS48)、色データkが所定の最小値k=Kminよりも小さくならない範囲で(ステップS49)、ステップS45の処理を繰り返す。なお、微小変化分Δkは、第1逆変換テーブルを構成する色データkの任意のデータ間隔であり、例えば、色データkが0%〜100%の範囲のデータとして設定される場合には、1%きざみ等に設定すればよく、色データkが0〜255の範囲のデータとして設定されるのであれば、値1きざみに設定すればよい。ステップS45の2回目の処理を行う場合には、色データk=Kmax−1=100−1=99に対する、(48)式の左辺の測色値であるXYZ値は、k=Kmax=100%とk=80%の値で測色され順変換テーブルとして格納されている各XYZデータの補間処理(内挿処理)により求めることができる。
【0155】
一方、ステップS49の処理において、Kmin>kと判定された場合には、目標値(X0,Y0,Z0)に対する色データcmykを画像出力部35の色再現域外のデータに指定し、ここでは、色データcmyksolの算出は行わない(ステップS50)。
【0156】
以上の処理を原シーン忠実再現テーブル42に供給される全ての刺激値データXYZ18を目標値(X0,Y0,Z0)として行うことにより、画像出力部35の色再現域内の色データである網%データcmyk46を得ることのできる刺激値データXYZ18に対して、最大となる色データkを与えた場合の色データcmyksolを求めることができる(ステップS51)。
【0157】
なお、ステップS44における色データkの固定は、最大値Kmaxではなく、その最小値Kmin(k=Kmin)に固定してもよく、その場合には、ステップS48の処理をk=k+Δkとし、ステップS49の処理をKmax<kとすればよい。また、色データkを任意の値とすることも容易に考えられる。例えば、ステップS48の処理は、任意の値に設定したkについて、k=k−Δkとk=k+Δkを交互に行えばよく、ステップS49の処理は、ステップS48の処理に対応させて行えばよい。
【0158】
さらに、ステップS49、S50において色再現域外データに指定された色データcmykは、この発明の要旨ではないので詳しくは説明しないが、前記特願平8−154584号に記載されているように、いわゆるガマットマッピング(Gamat Mapping )処理による色データCMYKの圧縮処理やクリッピング処理により、目標値(X0,Y0,Z0)に対する色データcmykの逆変換テーブルを作成することができる。
【0159】
このようにして、測色値信号18を構成する全ての目標値(X0,Y0,Z0)に対して求められたcmyk信号46の対応表を原シーン忠実再現テーブル42として格納しておくことで、任意の測色値信号18を画像出力部35の色再現範囲内のcmyk信号46に変換することができる。
【0160】
次に、cmyk変換部32を構成する標準条件再現テーブル44の作成について説明する。
【0161】
この標準条件再現テーブル44は、図1に示すように、デジタルカメラにより取り込まれたRGB信号12がマトリクス20により測色値信号(XYZ)18に変換され、この測色値信号XYZが色素濃度変換部30により色素濃度信号(cmy)28に変換され、この変換後の色素濃度信号(cmy)28を網%信号であるcmyk信号48に変換するためのテーブルである。
【0162】
また、標準条件再現テーブル44は、スキャナにより取り込まれたRGB信号14がテーブル24により測色値信号(XYZ)22に変換され、この変換後の測色値信号(XYZ)22が色素濃度変換部36により色素濃度信号(cmy)34に変換され、この変換後の色素濃度信号(cmy)34を網%信号であるcmyk信号に変換するためのテーブルである。この場合、色素濃度変換部36は、前述の色素濃度変換部30の作成過程で示した図7中、ステップS21b、S22dの方法と同様の過程で作成することができるので、その作成過程についての説明は省略する。
【0163】
標準条件再現テーブル44の作成処理を、図11に示すフローチャートを参照して説明する。
【0164】
まず、リバーサル原稿上に色素濃度cmyが格子状に配列された13×13×13個のカラーパッチを有するカラーチャートを作成する(ステップS61)。このカラーチャートは、c(シアン)m(マゼンタ)y(黄)が、それぞれ最小濃度から最大濃度まで13段階振られた各カラーパッチを有するものである。
【0165】
次に、このカラーチャートをスキャナのデフォルト分解条件で分解する、言い換えれば、透過原稿であるカラーチャートをスキャナで読み取り、デジタルデータ化する(ステップS62)。スキャナのデフォルト分解条件は、少なくとも、階調特性変換処理、色修正処理、UCR(下色除去)処理およびK版生成処理を含んでいる。ただし、分解するとき、カラーパッチのグレーの最小濃度を、cmyk網%が全て0%となるようにし、かつカラーパッチのグレーの最大濃度を、いわゆるベタとなるようなcmyk網%となるように設定する。
【0166】
そして、各カラーパッチの色素濃度とスキャナのデフォルト分解条件で分解したときのcmyk網%値を対応させて変換テーブル(ルックアップテーブル)を作成し(ステップS63)、この変換テーブルを標準条件再現テーブル44として設定する。なお、実際には、この標準条件再現テーブル44と補間処理により、最小濃度から最大濃度までの任意の色素濃度信号(cmy)28または色素濃度信号(cmy)34を所望のcmyk網%信号48に変換することができる。
【0167】
以上詳述したように、上述の実施の形態によれば、デジタルカメラで撮影した原シーンのについてのデバイス依存の画像信号であるRGB信号12および(または)スキャナで読み取ったリバーサルフイルム原稿についてのデバイス依存の画像信号であるRGB信号14を測色変換部16により、一旦、各々、デバイス非依存の画像信号である測色値信号(XYZ)18、22に変換する。そして、変換後の測色値信号(XYZ)18、22を測色的END順変換・逆変換マトリクス50により、測色的END信号52に順変換する。
【0168】
この順変換過程では、色度図上における3つの原刺激RGBの色度座標Rxyz、Gxyz、Bxyzを、リバーサルフイルムの色再現領域71を包含する領域73(図4参照)に設定している。
【0169】
そして、オートセットアップの計算を短時間で行うために、測色的END信号52のうち、間引き処理した信号のみをラフデータ作成部54により作成している。
【0170】
オートセットアップ部58では、このグレー条件が規定された測色的END信号52により既存のオートセットアップ処理を行う。例えば、測色的END信号52についてのヒストグラムを作成した後、累積ヒストグラムを作成し、この累積ヒストグラムの所定値、例えば、0.1%値をハイライト点濃度(ハイライトセットアップ点)D1に設定するとともに、98%値をシャドー点濃度(シャドーセットアップ点)D2に設定する。
【0171】
そして、このハイライトセットアップ点D1とシャドーセットアップ点D2とを、それぞれ、ハイライト網%、シャドー網%に対応する測色的END値、例えば、DHを0.1およびDSを2.0となるような階調変換カーブ74に基づき、END・END変換部64により、測色的END信号52を測色的END信号66に変換する。
【0172】
そして、オートセットアップによって設定された階調変換カーブ74に基づき、階調変換された測色値END信号66を測色的END順変換・逆変換マトリクス50により測色的信号18に逆変換して戻すようにしている。
【0173】
したがって、結果として、デジタルカメラから得られるRGB信号12に対して、スキャナに搭載されている既存の濃度信号に基づいたオートセットアップ処理を行うことができるという効果が達成される。なお、ハイライト点濃度(ハイライトセットアップ点)D1とシャドー点濃度(シャドーセットアップ点)D2の設定は、任意に所望の値に設定すること、言い換えれば、マニュアルで設定することも可能であり、いわゆる従来通りの濃度感覚のセットアップも可能である。
【0174】
なお、この発明は上述の実施の形態に限らず、この発明の要旨を逸脱することなく種々の構成を採り得ることはもちろんである。
【0175】
【発明の効果】
この発明によれば、測色値信号を一旦測色的END信号に変換しているので、従来通りの濃度の感覚でセットアップ処理を行うことができる。
【0176】
また、測色的END信号に変換する際、リバーサルフイルムの再現域内で変更するようにしているので、リバーサル原稿から得られる画像信号に対してセットアップ処理を行う、線走査読取装置、例えば、スキャナの既存のオート(またはマニュアル)セットアップ処理を利用することができるという効果が達成される。
【図面の簡単な説明】
【図1】この発明の一実施の形態が適用された画像信号処理装置の全体的な構成を示すブロック図である。
【図2】測色変換マトリクスの作成処理の説明に供されるフローチャートである。
【図3】測色的セットアップ部の説明に供されるフローチャートである。
【図4】CIE色度図上でのリバーサルフイルムの色再現域等の説明に供される色度図である。
【図5】デジタルカメラの測色信号をリバーサル原稿上の色素濃度信号に変換する変換過程の説明に供されるフローチャートである。
【図6】デジタルカメラの測色信号をリバーサル原稿上の色素濃度信号に変換する他の変換過程の説明に供されるフローチャートである。
【図7】色素濃度変換マトリクスの作成処理の説明に供されるフローチャートである。
【図8】ニュートン・ラフソン法の説明に供される線図である。
【図9】色素濃度変換マトリクスの作成処理の一部の説明に供されるフローチャートである。
【図10】原シーン忠実再現テーブルの作成処理の説明に供されるフローチャートである。
【図11】標準条件再現テーブルの作成処理の説明に供されるフローチャートである。
【符号の説明】
10…画像信号処理装置 12、14…RGB信号
16…測色変換部 18、22…測色値信号(XYZ)
20…測色変換マトリクス 24…測色変換テーブル
25…光源変更部 26…測色的セットアップ部
28、34…色素濃度信号(cmy) 30、36…色素濃度変換部
32…cmyk変換部 35…画像出力部
42…原シーン忠実再現テーブル 44…標準条件再現テーブル
46、48…cmyk信号 64…END・END変換部
Claims (4)
- 測色値信号をRGB3原色の表色系の信号であるRGB信号に変換した後、求めたRGB信号の各々の基礎刺激に対する比率を対数変換して測色的END信号に変換する表色系第1変換過程と、
前記測色的END信号から階調変換カーブを設定するセットアップ過程と、
この設定に基づき、前記測色的END信号を階調変換処理する階調変換過程と、
階調変換処理後の測色的END信号を逆対数変換してRGB信号にもどし、このRGB信号を測色値信号に変換する表色系第2変換過程と、
を有することを特徴とする画像信号処理方法。 - 請求項1記載の方法において、
前記表色系第1変換過程におけるRGB3原色の表色系は、リバーサルフイルムの色再現域を包含する原刺激を選択した領域とすることを特徴とする画像信号処理方法。 - 非リバーサル原稿より得られたデバイス依存の画像信号を測色値信号に変換する過程と、
前記測色値信号をRGB3原色の表色系の信号であるRGB信号に変換した後、求めたRGB信号の各々の基礎刺激に対する比率を対数変換して測色的END信号に変換する測色的END信号変換過程と、
前記測色的END信号から階調変換カーブを設定するセットアップ過程と、
この設定に基づき、前記測色的END信号を階調変換処理する階調変換過程と、
階調変換処理後の測色的END信号を逆対数変換してRGB信号にもどし、このRGB信号を測色値信号に変換する過程とを有し、
前記測色的END信号変換過程におけるRGB3原色の表色系は、リバーサルフイルムの色再現域を包含する原刺激を選択した領域または前記非リバーサルフイルム原稿の色再現域を包含する原刺激を選択した領域とすることを特徴とする画像信号処理方法。 - 測色値信号を変換して求められたRGB3原色の表色系の信号であるRGB信号の各々の基礎刺激に対する比率を対数変換して測色的END信号に変換する過程と、
前記測色的END信号から階調変換カーブを設定するセットアップ過程と、
この設定に基づき、前記測色的END信号を階調変換処理する階調変換過程と、
階調変換処理後の測色的END信号を逆対数変換してRGB信号にもどす過程と
を有することを特徴とする画像信号処理方法。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP31974996A JP3671097B2 (ja) | 1996-11-29 | 1996-11-29 | 画像信号処理方法 |
US08/980,367 US6278533B1 (en) | 1996-11-29 | 1997-11-28 | Method of processing image signal |
US09/852,633 US6668079B2 (en) | 1996-11-29 | 2001-05-11 | Method of processing image signal |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP31974996A JP3671097B2 (ja) | 1996-11-29 | 1996-11-29 | 画像信号処理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH10164352A JPH10164352A (ja) | 1998-06-19 |
JP3671097B2 true JP3671097B2 (ja) | 2005-07-13 |
Family
ID=18113749
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP31974996A Expired - Fee Related JP3671097B2 (ja) | 1996-11-29 | 1996-11-29 | 画像信号処理方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3671097B2 (ja) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002027272A (ja) | 2000-07-11 | 2002-01-25 | Fuji Photo Film Co Ltd | カラーチャートおよび色変換定義作成方法 |
JP4835900B2 (ja) * | 2004-11-30 | 2011-12-14 | Nkワークス株式会社 | デジタルカメラからの画像データのための画像処理方法と画像処理装置 |
JP5438944B2 (ja) * | 2007-10-23 | 2014-03-12 | 富士フイルム株式会社 | 濃度算出装置、濃度設定装置、濃度算出プログラム、および濃度設定プログラム |
JP5618211B2 (ja) | 2011-02-09 | 2014-11-05 | 株式会社リコー | 制御装置、画像形成装置、及び制御方法 |
JP5994460B2 (ja) * | 2012-07-31 | 2016-09-21 | 株式会社富士通ゼネラル | 二乗和ルート演算装置 |
-
1996
- 1996-11-29 JP JP31974996A patent/JP3671097B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JPH10164352A (ja) | 1998-06-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6668079B2 (en) | Method of processing image signal | |
US7206100B2 (en) | Image processing method and apparatus | |
Vrhel et al. | Color device calibration: A mathematical formulation | |
EP0675636B1 (en) | Extended density color printing | |
US5317425A (en) | Technique for use in conjunction with an imaging system for providing an appearance match between two images and for calibrating the system thereto | |
JPH06501144A (ja) | マルチトーン画像の表示方法 | |
US5781709A (en) | Method of and apparatus for generating proof | |
US20080158580A1 (en) | Method of converting image signal, method of determining original color of color image, and apparatus for converting image signal | |
JP3803441B2 (ja) | 色変換方法 | |
JP3304740B2 (ja) | 色シミュレーション装置 | |
JP3671097B2 (ja) | 画像信号処理方法 | |
US8064112B1 (en) | Color separation and reproduction method to control a printing process | |
GB2234410A (en) | Converting image densities between different types of reproduction process | |
JP3770671B2 (ja) | 色変換方法 | |
JPH10191077A (ja) | 画像信号処理方法 | |
JPH10215384A (ja) | 色信号処理方法 | |
US6909523B2 (en) | Information color system and printing method in the information color system | |
JP3360531B2 (ja) | 色再現域判定方法及び装置 | |
JP2001157073A (ja) | 画像信号変換方法、カラー画像の原色決定方法および画像信号変換装置 | |
JPH08163387A (ja) | カラー画像記録装置およびカラー画像記録システム | |
JPH11262033A (ja) | カラー画像出力装置における階調補正方法 | |
JPH08214157A (ja) | 印刷プルーフ作成方法および装置 | |
JP3874976B2 (ja) | 色信号変換装置 | |
JP2001277463A (ja) | カラー原稿の色再現性に優れた4色カラー網点印刷画像を製作するための各色版の網点面積比(c,m,y,bk)の設定法 | |
JP3305070B2 (ja) | カラー画像再現システム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20040831 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050111 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050308 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20050408 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050418 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080422 Year of fee payment: 3 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090422 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090422 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100422 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110422 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120422 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130422 Year of fee payment: 8 |
|
LAPS | Cancellation because of no payment of annual fees |