JP3835059B2 - カラー画像処理方法およびその装置 - Google Patents
カラー画像処理方法およびその装置 Download PDFInfo
- Publication number
- JP3835059B2 JP3835059B2 JP13821199A JP13821199A JP3835059B2 JP 3835059 B2 JP3835059 B2 JP 3835059B2 JP 13821199 A JP13821199 A JP 13821199A JP 13821199 A JP13821199 A JP 13821199A JP 3835059 B2 JP3835059 B2 JP 3835059B2
- Authority
- JP
- Japan
- Prior art keywords
- interpolation
- color
- gradation
- image
- output
- 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
Description
【０００１】
【発明の属する技術分野】
本発明は、カラー画像処理技術一般、特にカラーマネジメント技術に関する。
【０００２】
【従来の技術】
近年、カラーマネジメント技術の分野では、スキャナ、モニタ、プリンタなどのデバイスの間のカラーマッチングを図るためコンピュータのOSレベルでソフトウエア的な仕組みが組み込まれるようになっている。 このCMS（カラーマネジメントシステム）の機能は、あらかじめ各カラーデバイス毎に「プロファイル」とよばれるデバイス色変換特性情報を用意しておき、色変換エンジンがカラー画像を前記プロファイルに基づき色変換実行することにより行われる。プロファイルと色変換エンジンの構成として最近使われているのが多次元補間法である。
【０００３】
この方法は色変換における入出力関係をＭ次元の入力からN次元の出力空間への数学的変換と見なし、色変換を計算でなく色変換テーブルのルックアップで実行するものである。現在使用されている多くのカラー画像の場合、画素値が（RGB）あるいは（CIE-LAB）あるいは（CMYK）など３変数、４変数にて表現されるためＭ、Nは、３あるいは４である。しかし将来使用が想定されるカラーマルチスペクトル画像の場合にはＭ，Ｎは１０次元以上の可能性もある。ここでは説明の都合上、色変換を、３入力１出力の変換と考えることとする。テーブル法は、３次元の場合各色信号が８ビットあるとすると、そのままでは1600万色にもおよぶフルカラーの色の各々に対する膨大なサイズのテーブルを構成する必要がある。そこで例えば入力色が（RGB）空間内にある場合、まず、入力色空間を格子に区切り、格子点上での出力値を３次元テーブルf(i,j,k)に蓄積しこれを「プロファイル」とする。
【０００４】
実際の色変換に際しては、色変換エンジンが画像の各画素のカラー値（RGB）を前記３次元テーブルの３次元空間内で近傍格子点での出力値を探し、それらを重み付け加算する補間方法にて出力を計算する。この手法はテーブルサイズを大きくすることにより変換精度がいくらでも高められること、変換速度が高速でかつ色変換内容によらず一定であること、などの理由で広く使われるに至った。図１７に、この色変換エンジンの典型的な構成図を示す。この図では、色空間の３次元補間に四面体（三角錐）を用いているが、同図の下部に示すように立方体（８点）、三角柱、斜三角柱（６点）、ピラミッド（５点）など様々な３次元補間方法が従来から開発されている。
【０００５】
さて、３次元補間方法では補間誤差の発生は避けられない。特に色の値同志の絶対的な誤差よりも、画像の滑らかな階調を色変換した結果に生じる補間結果の暴れ現象（リップル）が人間の視覚上非常に影響が大きい。このリップルについては、例えば、本出願人による文献１：金森：「３Dテーブルを用いた色変換における補間誤差解析」（カラーフォーラム JAPAN’98）に解説されている。グレイの階調上ではおなじガンマ変換を立方体を用いる３次元補間する場合でも色空間全体を定義域とする大域的（グローバル）ガンマ変換では変換後の階調カーブが折れ線近似される補間特性を示すが色空間内のグレイ階調部分だけを局所的にガンマ変換する場合には補間結果で階調が反転する程の補間曲線の暴れ現象を呈し（リップル）しかもテーブルサイズを大きくしてもリップル振幅は変化しないという特異現象も見つかっている。
【０００６】
この例からもわかるように一般に線形の多次元補間方法では、多くの格子点情報を用いる必要があり、それが画像の２次元平面内の滑らかな階調上に隠れた高次元特性を顕現する結果リップルが出現し階調特性が悪化することが多かった。また、リップル現象は、現在３次元補間において既に発生しているが、より高次元の４次元色空間やマルチスペクトル画像を扱う際の１０次元色空間になれば一層深刻な問題になることは容易に想像できる。
【０００７】
一方、文献２：James M.Kasson：“Performing color space conversions with three-dimentional linear interpolation”、Journal of Electronic imaging ／Jul 1995 Vol.4(3) pp.226-249)におけるディスフェノイド補間や特開平８−８３３３０号公報での体心構造を用いた四面体補間方法では四面体群を同一形状としたこと、最も補間誤差の大きい立方体の体心に格子点を位置させること、格子点数を各軸２倍に増加すると全体で８倍になるところを２倍の増加で済むことなどの特徴を持っており、良好な補間特性を持つとされている。特に、四面体稜線が色空間内で多くの方向に発生するため、画像の入力方向に対して稜線方向が合致する度合いが増え、先に述べたように大きなリップルを発生することが少ない効果がある。
【０００８】
しかし、体心構造では、RGB立方体の外部に不必要な外挿点を持つ必要があること、空間分割が複雑で処理時間がかかる、などの欠点がある。また、稜線方向がいかに増えても、それをはずれる入力は常に作り得るため問題の完全な解決にはなっていない。また３次元においても複雑な空間分割手法を使用する必要がありマルチスペクトルを用いた高次元でのカラー処理においては計算量が多大なものとなりやすい課題もある。
【０００９】
また、従来、特開平６−３１１３５５号公報では、入力色によってトライ・リニア補間と四面体補間を切り替えたり、特開平７−９５４２３号公報ではトライ・リニア補間、三角柱補間、ピラミッド補間を切り替えたり、色空間内で入力する色の位置によって異なる補間方法を採用することによって全体での色変換誤差を小さくする方法もある。
【００１０】
しかし、これらの方法は、入力した色の空間内の位置により補間方法が決定されるため種々の非線形性色変換において汎用的に対処できるものではない。
【００１１】
【発明が解決しようとする課題】
以上のように従来の多次元色空間補間法を用いた色変換においては、様々な色変換、様々な画像入力に対応して常に階調性が良好な補間方法を実行することが要求されているが、あらゆる画像、あるいは今後予想される多次元カラー処理において充分対応可能な解決策がなかった。
【００１２】
本発明は、上記の課題を解決するもので、カラーマネジメント分野でカラー画像を多次元補間方法で色変換する場合において、色変換自身が非線形性をもつ場合には、カラー画像上の色空間内における階調（グラデーション）方向を画像の画面内の部分領域毎に検出し、該部分領域ごとに補間方法を逐次制御すること、より詳細には色空間内の補間において入力点を囲む全ての格子点を使うのではなく画像の階調方向に存在する少数の格子点を選択して補間することにより多次元補間における非線形性を下げ、リップルを生じない階調性を重視した色変換を実施することができるカラー画像処理方法を提供することを目的とする。
【００１３】
【課題を解決するための手段】
この課題を解決するため本発明は、カラー画像を走査し各画素値を入力する画素入力部と、前記画素を含む画像上の部分領域毎で、各色版毎に輝度勾配の大きさを判定し当該階調の方向性を色空間内で離散化した方向軸ごとの方向重みを計算する階調方向判定部と、前記階調方向判定部で計算した方向重みを画像面内で連続的に補間し階調補間係数を出力する階調方向係数補間部と、色変換テーブルを用いて前記画素値に対して色変換を実施する複数の相異なる色空間補間部と、前記複数の色空間補間部からの出力を前記階調方向係数補間部から出力される階調方向係数により混合あるいは選択する出力制御部とを備える。
【００１４】
これにより、いかなる色変換、いかなる画像入力についても画像全体でリップルを生じない階調性を重視した色変換を実施することができるものである。
【００１５】
本発明の第１の実施の形態は、カラー画像を走査し各画素値を入力し、前記画素値の色空間内の位置情報から色変換テーブルを参照し多次元補間方法を用いて色変換する方法において、前記画素値を含むカラー画像を複数の分割領域内で各色版毎に輝度勾配の大きさを判定し、当該階調の方向について色空間内で離散化した方向軸ごとの方向重み係数を計算することにより、画像の階調方向情報を求め、前記階調方向情報を用いて前記多次元補間方法の制御することを特徴とするカラー画像処理方法であり、複数の補間方法から階調方向に合致した補間方法を特に重視して出力するため不必要に多くの格子点を使うために発生していた非線形性を低下させ、出力階調にリップルなどのアーティファクトが低減されるという作用を有する。また、画像の階調の方向を色空間内の方向として表現できる。
【００１７】
好ましくは、カラー画像を分割領域内での大評点近傍の階調の離散化した方向軸への方向重みを前記分割領域間で連続的に補間することにより、複数の異なる多次元補間方法を画像の階調方向によって画像面内で逐次切り替えても画像の階調の連続性が崩れないという作用を有する。
【００１８】
本発明の第２の実施の形態は、カラー画像を走査し各画素値を入力する画素入力部と、前記画素を含む画像上の部分領域毎で、各色版毎に輝度勾配の大きさを判定し当該階調の方向性を色空間内で離散化した方向軸ごとの方向重みを計算する階調方向判定部と、前記階調方向判定部で計算した方向重みを画像面内で連続的に補間し階調補間係数を出力する階調方向係数補間部と、色変換テーブルを用いて前記画素値に対して色変換を実施する複数の相異なる色空間補間部と、前記複数の色空間補間部からの出力を前記階調方向係数補間部から出力される階調方向係数により混合あるいは選択する出力制御部とからなるものであり、複数の補間方法から階調方向に合致した補間方法を特に重視して出力するため出力階調においてリップルが削減されるという作用を有する。
【００１９】
本発明の第３の実施の形態は、カラー画像を走査し各画素値を入力する画素入力部と、前記画素を含む画像上の部分領域毎で、各色版毎に輝度勾配の大きさを判定し当該階調の方向性を色空間内で離散化した方向軸ごとの方向重みを計算する階調方向判定部と、前記階調方向判定部で計算した方向重みを画像面内で連続的に補間し、階調補間係数を出力する階調方向係数補間部と、色変換テーブルを用いて前記画素値に対して色変換を実施する複数の相異なる色空間補間部と、前記複数の色空間補間部からの補間出力を比較する出力比較部と、前記複数の色空間補間部からの補間出力を出力比較部での判定結果に基づき前記階調方向係数により混合あるいは選択する出力制御部とからなるものであり、色変換自体に非線形性がある場合に限り複数の補間方法から階調方向に合致した補間方法を特に重視して出力するため、処理速度を落とさずに出力階調にリップルが削減されるという作用を有する。
【００２０】
本発明の第４の実施の形態は、カラー画像を走査し各画素値を入力する画素入力部と、前記画素を含む画像上の部分領域毎で、各色版毎に輝度勾配の大きさを判定し当該階調の方向性を色空間内で離散化した方向軸ごとの方向重みを計算する調方向判定部と、前記方向判定部で計算した方向重みを画像面内で連続的に補間し階調補間係数を出力する階調方向係数補間部と、色変換テーブルを用いて前記画素値に対して色変換を実施する複数の相異なる色空間補間部と、前記色変換テーブルのうち非線形性が強い領域を記述するテーブル非線形性マップと、複数の色空間補間部からの出力を前記テーブル非線形性マップに基づき前記階調方向係数により混合あるいは選択する出力制御部とからなるものであり、色変換自体に非線形性がある場合に限り複数の補間方法から階調方向に合致した補間方法を特に重視して出力するため、処理速度を落とさずに出力階調にリップルが削減されるという作用を有する。
【００２２】
好ましくは、多次元テーブルに対してあらかじめ用意したテスト信号を入力して異なる多次元補間を実施し、それら出力のばらつきに基づいて色変換テーブル内の非線形強度の大きい部分を判定し、ファイル化するものであり、あらかじめ多次元色変換テーブル内で非線型性強度の強い部分を判定し対処することができるという作用を有する。
【００２３】
以下、本発明の実施の形態について図１から図１５を用いて説明する。
【００２４】
(実施の形態１)
図１に本発明の実施の形態１におけるカラー画像処理装置のブロック構成図を示し説明する。本実施例において、カラー入力画像１０１として想定しているものは、１画素が（R,G,B）の３変数信号で規定される画像であり、３種の画像としてR版、G版、B版と呼ぶものとする。
【００２５】
図１において、１０２は、カラー画像１０１を走査し各画素値を入力する画素入力部、１０４は画素を含む画像上の部分領域毎で階調の色空間内での方向性を判定し各領域の方向重みｕを求める階調方向判定部、１０５は階調方向判定部１０４での方向重みｕを画像面内で連続的に補間し階調方向係数Wi(m,n)を出力する階調方向係数補間部、１０７は多次元色変換テーブル１０３を用いて前記画素値に対して色変換を実施する複数の相異なる種類の色空間補間部、１０６は複数の色空間補間部１０７からの出力を階調方向係数補間部１０５から出力される階調方向係数Wi(m,n)によって混合あるいは選択し、出力カラー画像１０８を出力する出力制御部である。
【００２６】
上記のように構成されたカラー画像処理装置の動作を以下に説明する。
【００２７】
カラー入力画像１０１における注目画素の画面上のｘｙ座標を(m,n)、カラー値を(R,G,B)とする。画素入力部１０２は、注目画素のカラー値（R,G,B）を各信号を８ビット構成とした場合、上位信号３ビットの組（Rh,Gh,Bh）と下位信号５ビットの組（Rl,Gl,Bl）に分割し、それぞれ色変換テーブル１０３と色空間補間部１０７に供給する。
【００２８】
画素入力部１０２にて分割された上位信号（Rh,Gh,Bh）にて色変換テーブル１０３をルックアップし、下位信号（Rl,Gl,Bl）とで補間演算を実施する。色変換テーブル１０３には、上位信号（Rh,Gh,Bh）に対応して全RGB色空間を８×８×８＝２５６に分割した頂点である９×９×９個の格子点での色変換出力が格納されている。ここで用いた３ビットと５ビットというビット分割方法は色変換テーブル１０３のサイズ、すなわち格子点数（格子点間距離：ここでは２5＝３２）を決定しているが実際には任意である。色変換テーブル１０３は、色再現を実行する各種色補正を目的とする数式（モデル式）にて表現できる色変換を表現できるのみならず、カラー機器の測色によって数式モデルを持たない自由度の高いテーブルを作成することも可能である。
【００２９】
また、色変換の入出力の定義も様々であり、従来は積分濃度空間（Dr,DgDb）からインク（トナー）信号（C,M,Y,K）へのプリンタ用の色変換が一般的であったが、カラーマネジメントの世界では(R,G,B)から標準色空間（L,A,B）（X,Y,Z）への変換、あるいは逆変換などが対象とされる。本発明に関しては（R,G,B）から任意のある出力色ｆへの色変換を考えることとする。すなわち色変換テーブルには出力１色分についてはｆ(R,G,B)という３変数関数が格納されていることになる。
【００３０】
さて、画素入力部１０２からの上位信号（Rh,Gh,Bh）は、色変換テーブル１０３の８×８×８の分割された立方体の１つを指定する。つまり上位信号によって入力色を内包する最近傍点８個が見つかり、そこでの出力ｆが得られる。この８点を用いて補間する方法がもっとも一般的な３次元補間方法であるトライ・リニア補間である。
【００３１】
色空間補間部１０７のうち、（０）は、このトライ・リニア補間方法を実行する部分であり以下その機能について図２を用いて簡単に説明する。得られた８点で構成される立方体内で入力信号を表現するのは下位信号（Rl,Gl,Bl）であり、３次元補間処理では、これらを重みとして用いる。ここで図２に示すように入力色をPとし、Pを含む色空間内の立方体をABCD-EFGHとし、Ｐの位置を（Ｒh＋Ｒl，Ｇh＋Ｂl，Ｂh＋Ｂl）とすると、
【００３２】
【数１】
【００３３】
であるから、Ｐでの補間出力ｆIが以下の（数２）で表現できる。
【００３４】
【数２】
【００３５】
ただし、その際の重み係数は（表１）の通りである。
【００３６】
【表１】
【００３７】
トライ・リニア補間は、公知技術として従来から使用されている。しかし、８点を用いる補間法であるため一般には精度が高いといわれながら、非線型色変換では補間特性にリップルを生ずることが知られている。この原因は表１の各重み係数が３次式となっていることに関係がある。この課題の解決のため従来から利用されてきた方法が四面体補間である。四面体補間とは、色空間で入力色を含む４点を用いた線形補間の総称であるが、最も一般的なのは前記の立方体ABCD-EFGHを対角線を境界の陵線として６個の体積が等しい四面体分割を行い、いずれの四面体領域に含まれるかを判定した上で当該の４点を用いて補間するダイアゴナル・テトラヘドロン補間方法である。ダイアゴナル・テトラヘドロン補間では補間立体内では線形補間が行われるためリップルが発生しない上、用いる点数が３次元補間では最小数である４点であるため計算量も少ないという良好な性質を持っている。しかし３次元空間で等方的な８点補間と異なり、上記の良好な性質は方向性により特性が異なる。そこで本発明ではダイアゴナル・テトラヘドロン補間について方向性を考えて４種類を用いている。つまり、図１の色空間補間部１０７の（１）から色空間補間部（４）までは、全てダイアゴナル・テトラヘドロン補間を実行するが、図３、および（表２）の(D1)、（D2）,（D3）,(D4)、に示すようにダイアゴナルテトラヘドロンの分割方法である対角線方向が４種類に異なっている。
【００３８】
【表２】
【００３９】
ここでは（D１）タイプ則ち立方体ABCD-EFGHのA-G対角線に沿った分割を例として補間方法を説明することとするが、他の３タイプの場合にも説明は同様である。
【００４０】
図４は、ダイアゴナル・テトラヘドロン補間方法（D１）の６分割方法を説明する図である。（表３）に従い、下位信号（Rl,Gl,Bl）の大小関係は６通り存在し、これが６分割四面体のいずれかを決定する。Rl=Gl=Blの時にはArea０からArea５までのいずれに含ませても良いがここではArea０とした。
【００４１】
【表３】
【００４２】
分割された四面体が決定した後に、以下のように４点補間が実行される。
【００４３】
【数３】
【００４４】
この際の重み係数は以下の（表４）のように決定する。
【００４５】
【表４】
【００４６】
次に、階調方向判定部１０４、階調方向係数補間部１０５および出力制御部１０６について説明する
【００４７】
階調方向判定部１０４にて入力画像の階調の方向を判定し、方向重みｕを求める。階調方向係数補間部１０５は、階調方向判定部からの方向重みｕを画像画面内で連続的に補間し、階調補間係数Ｗi(m,n)を出力する。
【００４８】
出力制御部１０６において当該方向に近い方向の分割線を持つ補間方法に大きな階調方向係数を与え、相異なる方向の分割線を持つ補間方法には小さな階調方向係数を与えて出力結果を加算する。出力制御部１０６では、色空間補間部（０）から（４）までの出力をｆI(i)(i=0-4)と表現し、階調方向係数Wiを用いて混合する。
【００４９】
【数４】
【００５０】
この階調方向係数Wiは、各補間方法ごとに与えられるものであり、
【００５１】
【数５】
【００５２】
を満たすものとする。
【００５３】
次に、階調方向係数Wiの計算の仕方について、図８を用いて説明する。まず、カラー画像を、Red,Green,Blueの３つの色版からなる単色版の画像と考え、各々の版毎に局所的な画像の階調の大きさを判定する。はじめに図９に示すように、画像を部分領域に分割しておき、これから色変換を行う対象画素位置(m,n)を含む領域内で代表点（m0,n0）を決定し、その場所で階調方向を算出する。すなわち画像に対して輝度勾配ベクトルDを算出するフィルタリングを実施すればよく、ここでは偏微分を用いるが現実にはこの操作はデジタル的な差分操作となる。
【００５４】
【数６】
【００５５】
この階調方向ベクトルＤは、階調の急激さとは無関係な量にするため長さ＝１に正規化するのが望ましい。次に階調方向ベクトルＤの方向性を調べる。これは簡易的には、７個の方向ベクトルVjと内積をとることで実現できる。７個の方向ベクトルVとは図８で示すV1からV7の７方向であり、各ベクトルのRGB座標成分は以下の（表５）の通りである。
【００５６】
【数７】
【００５７】
【表５】
【００５８】
以上求めた各成分ごとの７個の重み係数ujについて、u1からu4までは色空間補間部（１）から（４）までの重みとしてそのまま用いられ、u5からu7までは加算され色空間補間部（０）の重みとして用いられる。
【００５９】
【数８】
【００６０】
そして、W0とW1、W2, W3, W4の階調方向係数Wiが領域Bnごとに決定される。
【００６１】
この値をそのまま出力し、分割領域ごとに用いても良いが、階調方向係数Wiが領域間で大きく変化した場合、分割領域ごとに色や階調の段差が発生する可能性があり、これを防ぐため階調方向係数Wiが画像内で連続性を保つように、階調方向係数補間部１０５で画面内のバイ・リニア補間によって階調方向係数Wi(m,n)を補間出力している。
【００６２】
【数９】
【００６３】
【数１０】
【００６４】
ただし、階調方向係数の重み係数ωP、ωQ、ωR、ωSは、以下の（表６）に示す通りである。ここでBX,BYは分割領域の幅、高さを示す。また上記の内挿補間は、画像内で各Wiごとに実施するものとする。
【００６５】
【表６】
【００６６】
図９、図１０に示すようにこの方法では、重み係数は画像の周辺に内挿不可能な領域を残すが、ここは外挿あるいは分割領域での代表点(m0,n0)の値をそのまま使用することも考えられる。また代表点(m0,n0)を各分割領域の左上にとり、画像の右端では仮想的な固定値を採用することにより外挿を不必要にすることも可能である。
【００６７】
次に、実際の画像処理の流れについて図１１を用いて説明する。この例では、カラー画像は入力出力とも（RGB）点順次画像を想定しており、色変換処理も点順次に実行されていくこととする。
【００６８】
処理開始後、Ｓ１１０１のステップにおいて、カラー画像は図９のように各領域に分割される。Ｓ１１０２のステップで領域内の代表点(m0,n0)における階調方向判定が行われ、各分割領域ごとの階調方向係数Wiを求める。ステップＳ１１０３において色変換に用いる３次元テーブルを読み込む。ここまでは画像処理の前処理として実行される。
【００６９】
次に、カラー画像の各画素を走査しながら、全画素の処理が終了するまでステップＳ１１０６以降ステップＳ１１１０までの処理を各画素ごとに繰り返す。最初にステップＳ１１０６において、画素の位置(m,n)および画素値（RGB）が入力される。ステップS１１０７において、当該画素位置(m,n)での階調方向係数Wi(m,n)が、前記の各分割領域ごとの階調方向係数Wiを補間して求められる。
【００７０】
次に、ステップＳ１１０８において、異なる色空間補間処理を並列に実行して補間出力fI(0)からfI(4)を求める。ステップＳ１１０９においてＳ１１０７で求めた階調方向係数Wi(m,n)を用いて補間出力が混合される。ステップ１１１０において出力画素(R’G’B’)が出力画像に挿入される。
【００７１】
処理は、そのまま繰り返しのためステップＳ１１０４へ戻り、全画素の処理が終了した場合にはステップＳ１１０５の判定により処理が終了する。
【００７２】
なお、以上の説明では、色空間補間部１０７をトライリニア補間とダイアゴナル・テトラヘドロン補間４種類ということで説明したが、方向性を持つ分割空間であれば他の補間方法での組み合わせでも同様に実施可能であり種類も４種類には限らないことはもちろんである。
【００７３】
以上まで、本発明の構成および動作を説明してきたが、ここで階調方向によって補間内容を変更するという本発明の原理の説明を行う。まず何故トライリニア補間とダイアゴナルテトラヘドロン補間（しかも４種）を併用しているのか、さらに画像の階調方向を考慮してそれら重み付けしているのかを説明する。
【００７４】
従来、トライ・リニア補間とダイアゴナル・テトラヘドロン補間のどちらが高精度な補間方法であるかについては用いる色空間、用いる色変換が異なった場合、明解な説明がなかった。しかし、前述の文献１に解説されているように用いる色変換自身が１次元変換の線形結合にて表現される場合には補間結果は全ての線形補間で同一となり、出力結果が各成分ごとに折れ線状態になることが知られている。つまり、この場合にはトライ・リニア補間とダイアゴナル・テトラヘドロン補間両者の結果の差はなくリップルも存在しない。
【００７５】
このことを２次元の補間を例にとって図形的に説明する。はじめに２次元では、トライリニア補間は４点を用いるバイリニア補間に対応し、ダイアゴナル・テトラヘドロン補間は、正方形を２分割した三角形の３点補間に相当することに注意する。図５の５００において正方形ＡＢＣＤはＲＧ平面に設定された１つの補間区間であり、出力１成分をｆとしてこの２次元平面上に立てた矢印にて表現する。色変換が線形であるとは、この正方形内ではｆ（Ａ）、ｆ（Ｂ）、ｆ（Ｃ）、ｆ（Ｄ）が同一平面π上に存在することに相当する。従ってこの４点からの線形補間であるバイリニア補間による補間値は５０１に示すように同一平面π上に拘束されており、正方形ＡＢＣＤ内をいかなる軌跡で入力された色についても出力は直線となる。
【００７６】
次に、図５の５０２，５０３は、A-C分割、B-D分割という異なる種類の三角形分割における補間結果をしめす図である。いずれの分割方法の場合でも平面内の入力点が２種の三角形のいずれに属するかを判定した上で各三角形内で補間される。三角形内での補間値はこの平面πを三角形ABC、または三角形ADCに分割して補間することに相当するが、元々平面を分割するのであるからA-C分割、B-D分割いずれの三角形分割でもうバイリニア補間と同一の補間結果となることは明らかである。このように色変換自身が非線形性を持たない場合には、補間方法の如何にかかわらず、結果が理想値からずれたり出力が曲線になる現象（リップル）が出現することもない。
【００７７】
次に、特開平６−３１１３５５号公報にてダイアゴナル・テトラヘドロン補間がトライリニア補間よりも良好な結果を得たと言われている色変換がMIN演算の場合、特にダイアゴナル・テトラヘドロン補間の中でも原点を通る分割線を採用する（D1）型の場合を取り上げる。同じように２次元で考えると、MIN演算とは、２変数の組(R,G)の小さい方を出力するものであるから図６の６００に示すようにR≧Gの領域に相当する三角形ABC内では、MIN(R,G)＝Gであり、R<Gの領域に相当する三角形ADC内ではMIN(R,G)=Rであるから直線R=Gを稜線として２枚の三角形平面が折れて接続されているような平面を構成する。これがMIN演算の場合における真の計算値である。
【００７８】
次に、これを補間で演算することを考える。まず６０１に示すようにABCDの４点補間法では補間面が曲面を形成してしまいこれがリップルを発生する原因になる。次に、図６の６０２に示すA-C分割での三角形分割では、稜線がA-C方向に形成されるため補間結果は正しい結果と完全に同一になる。しかし、６０３のB-D分割ではB-D方向で分割された１つの三角形では値が０に固定されたままとなり、２つの分割方法で補間出力が大きく異ってしまう。
【００７９】
前記特開平６−３１１３５５号公報では、この場合を取り上げてA-C分割の三角形補間がABCD４点より良好な結果を得るとしている。しかし正確に言えばＭＩＮ演算の真値の演算において、A-Cの分割線が全く違う特性を示す２種の平面が接続された稜線となっており、この稜線と三角形の分割境界とが一致しているため補間誤差が０になり良好な補間が行われるというのが正しい。結局、非線形性の強い色変換であってもMIN演算のように真の計算値が知られている場合には最適な補間方法が何らかの指標によって定まる点が重要である。逆にいえば非線形性が強い場合において真の計算値が不明な場合にはいかなる分割・補間をすれば良いかを知る指標が無くなるため即座に困難に直面することになる。
【００８０】
この困難を生む典型的な例が、前述の文献１で取り上げた局所ガンマ変換のような非線形変換である。この変換では、非線形性の強い色変換テーブルだけが与えられるので、色変換に当たってはこれを３次元補間する必要がある。すなわち先の例では、ｆ（A）、ｆ（B）、ｆ（C）、ｆ（D）という格子点上での値だけが非線形性が強い状態で与えられ、その間を接続する真の計算値は全く与えられない。非線形性が強いという意味は２次元では図７の７００に示すようにABCDの各点で対角線位置にあるA，Cで小さな値、B,Dでは大きな値をとるような場合である。これをバイリニア４点補間すると７０１に示すようにｆ軸Ｒ軸Ｇ軸で作られる空間において複雑な２次曲面を形成し、リップルの原因となる。
【００８１】
また、３点補間では７０２と７０３に示すようにA-C分割とB-D分割では、前者がA-C方向が深い谷を形成するのに対して後者はB-Dが稜線を形成して全く異なる補間結果となるが、いずれもリップルの代わりに激しい凹凸が発生する原因となる。つまり４点補間でも２種類の３点補間でも補間結果は大きく異なる上に、どれも良好な補間とはいえないという困難な状況となる。
【００８２】
本発明の基本的な考え方は、ここに入力画像の階調性の方向の情報を入れて最適な補間を決定するものである。たとえば図７の７００において鎖線(1)(2)が入力階調を示すものとすると、(1)のようなA-C方向の階調が入力された場合にはA-C分割の３点補間を採用し、(2)のようなB-D方向の階調が入力された場合にはB-D分割の３点補間を採用するものということができる。
【００８３】
次に、入力階調方向が求まった後、いかなる補間方法を使用するべきか、おなじく２次元の場合を例に図１２を用いて説明する。図１２では、入力階調として正方形ＡＢＣＤの中心点を通る４種の方向を設定し、入力点を黒丸で、その点での補間出力値を黒丸上の矢印で表示している。１２００と１２０１は入力階調方向がそれぞれＥＦ方向、ＧＨ方向の場合であり、このときは両方とも同一のバイリニア補間を使用すればよい。
【００８４】
図７の７０１では、複雑な２次曲面を呈しているかのように見えるバイリニア補間面であるが、１２００で示すように入力方向がＥＦ方向のときは補間結果がＡＢの中点ＥとＣＤの中点Ｆ上での値ｆ（Ｅ）とｆ（Ｆ）を結ぶ直線となりリップルが現れない性質がある。同様に１２０１の入力階調方向がＧＨの時もバイリニア補間の結果、補間値がｆ（Ｇ）とｆ（Ｈ）を結ぶ直線上に位置する。
【００８５】
次に１２０３では入力階調がＡＣ方向の場合なので、ＡＣ分割の三角形補間を用い、１２０４では入力階調方向がＢＤなのでＢＤ分割の三角形補間を用いることにより、いずれの場合も補間結果がそれぞれｆ（Ａ）とｆ（Ｃ）、あるいはｆ（Ｂ）とｆ（Ｄ）を結ぶ直線上に位置するためリップルは現れない。以上のように２次元の場合には、階調方向４種類（１）ＥＦ、（２）ＧＨ、（３）ＡＣ、（４）ＢＤについて（１）（２）はバイリニア補間、（３）と（４）は異なる三角形補間を使用すればよい。もちろん画像上の各画素で異なる補間方法を次々に変更していくと画像上にノイズを生む原因となる。このため補間方法の制御としては、最適な補間を１種選択するよりは重みとして用いるほが望ましい。つまり画像の階調方向から考えて最適な補間方法の結果には大きな重みを与え、そうでない補間結果には小さい重みを与えるように何種類からの補間方法に対して一種の重みづけを行う。
【００８６】
以上の議論を３次元に拡張すると、図８の方向ベクトルでＶ５とＶ６とＶ７のような立方体軸に平行の入力階調の場合にはいずれもトライリニア補間を用いるのがよく、その他の対角方向の入力Ｖ１、Ｖ２、Ｖ３、Ｖ４では、各々図３のＤ１，Ｄ２、Ｄ３、Ｄ４方法の分割のダイアゴナルテトラテドロン補間を用いるのが良いという結論が導かれる。このため、（数８）のように方向重みｕ5、ｕ6、ｕ7、ｕ8についてはそれらを加算したものをトライリニア補間用の階調方向係数Ｗ0として用い、他のＷ1、Ｗ2、Ｗ3、Ｗ4はｕ1、ｕ2、ｕ3、ｕ4をそのまま用いている。
【００８７】
以上のように、本発明の第１の実施形態では、いかなる非線形色変換においても、入力画像の階調を求め、数種類の補間方法を実施してその結果に重み付けすることで出力階調でのリップルを抑制するものである。
【００８８】
（実施の形態２）
図１３に本発明の実施の形態２の画像処理装置のブロック図を示し説明する。
【００８９】
実施の形態１では、異なる色空間補間部１０７からの出力は入力画像１０１の画素位置から計算される階調方向係数によって画素ごとに常に重み付け混合されていたが、この重みづけ混合の計算は画像の全画素について実行されるため計算処理負担が大きい。また実際には非線型の高い色変換やそれによる画像の階調の劣化は常に発生するわけではない。
【００９０】
そこで、本実施の形態では、５種の色空間補間部(0)〜(4)での出力を非線形性判定部１３０１にて比較しこれらが大きく異なる場合、すなわち非線形性が高い場合にのみ、重み付け混合を実施することによって、処理の負担を下げるものである。非線形判定部１３０１では、以下のような演算を行う。ここでTHRESHは固定の閾値、ｆI（i）、ｆI（ｊ）はｆI(0)から ｆI(4)のいずれかを示す。MAXは最大値演算を示す。
【００９１】
つまり、ある入力色について異なる分割方向や異なる種類の補間を用いたとき最大でどの程度の変動があるかを調査し、ある閾値以上であったら非線形性のが強いと判断し、結果resultをYESと判定し、そのほかはNOと判定し、この結果を一緒に出力制御部１０６に送ることにより出力混合を行うのかどうかを決定する。
【００９２】
【数１１】
【００９３】
この方法では、非線形性判定部１３０１がYES/NOの判定を行う境界が生じるため出力の切り替えによって階調飛びが生じるように思われるるが、閾値を人間の視覚が階調とびを認識できる限界値に選んでおけば補間方法による差が閾値を超える前までは出力混合を停止して、超えた時点で出力混合を開始してもそのギャップは前記閾値よりは小さいことが保証されるため認識できる大きな階調飛びは発生しない。
【００９４】
（実施の形態３）
図１４に本発明の実施の形態３における画像処理装置のブロック図を示し説明する。実施の形態３は、実施の形態２における非線形性判定部を無くし、その代わりにあらかじめ用意したテーブル非線形性マップ１４０１の情報を元に、入力画素の色信号値が色変換テーブル内での非線形性の大きい部分に入った場合にのみ、重み付け混合を実施することによって処理の負担を下げるものである。この場合、色変換テーブルを作成した段階でテーブル非線形性マップが一緒に作られ、この情報はICCプロファイルのようなカラープロファイルの１つの情報として蓄積しておくこともできる。
【００９５】
次に、図１５を用いてテーブル非線形性マップの作成方法を説明する。図１５はテーブル非線形性マップ作装置を示すブロック図である。色変換のための３次元テーブルに対応して入力格子点の色空間内の位置において色変換の非線形性、主として分割方向による補間結果の差の大きさを判定し当該色での補間出力の混合が必要かどうかを判定するものである。作成方法は、第２の実施例に似た構成で、対象となる色変換テーブルを用いてテスト用カラー画像を色変換する。異なる色空間補間部で出力された補間出力どうしを非線形性判別部１３０１で（数１１）と同一の判定を行い、（Rh,Gh,Bh）をアドレス入力とした時の非線形性の結果として、判定結果YES/NOをテーブル非線形性マップ１５０１に格納する。
【００９６】
（実施の形態４）
図１６に本発明の実施の形態４における画像処理装置のブロック図を示し説明する。実施の形態４は、実施の形態１における３次元色変換テーブルを異なる複数の色空間補間部１０７（０）〜（４）ごとに専用の色変換テーブル１６０１としたものである。本発明の実施の形態１から３では、１つの３次元色変換テーブルに対して画像の階調方向によって方向性をもつ色空間補間部１０７を使いわけることに利点を有するが、本実施の形態は３次元色変換テーブル自体をも使い分けることにより、画像の階調方向性の階調性を一層良好なものとすることができる。
【００９７】
尚、本発明は３次元のカラー画像を対象にして実施の形態を書いているが入力空間は一般に４次元以上のＭ次元でも同一の考え方でかまわない。出力空間がＮ次元の場合にはＭ入力１出力の個々の色変換部をＮ個並べるだけで実現できる。
【００９８】
従って本発明は、Ｍ入力Ｎ出力のマルチスペクトルカラー画像の色変換においても有効である。
【００９９】
なお、本発明は、階調判定の単位を画面上のブロック分割領域としたが、画像内の別のセグメンテーション方法によるものであってもよい。
【０１００】
【発明の効果】
以上のように本発明によれば、従来技術が色変換を各画素値のみを用いて行っていたため生ずる階調の乱れの問題に対して、カラー画像上の階調（グラデーション）方向を画面内の部分領域毎に検出し、該部分領域ごとに補間方法を逐次制御することにより、いかなる色変換、いかなる画像入力についても画像全体でリップルを生じない階調性を重視した色変換を実施することができるカラー画像処理装置を提供することができる。
多次元カラーマルチスペクトル画像の補間方法においては、入力次元が増大するため高次項が補間結果により深刻な影響を与え未知のリップルが発生しやくなることが考えられるが本発明の画像の階調方向を考慮するという考え方を利用すれば現実的に良好な補間演算が可能になる。
【図面の簡単な説明】
【図１】本発明の実施の形態１によるカラー画像処理装置のブロック構成図
【図２】トライ・リニア補間方法を示す図
【図３】４種の四面体分割方法を示す図
【図４】 D1分割方法でのダイアゴナル・テトラヘドロン補間方法を示す図
【図５】２次元の場合における非線形性が無い色変換補間を示す図
【図６】２次元の場合におけるMIN演算の色変換補間を示す図
【図７】２次元の場合における非線形性が強い色変換補間を示す図
【図８】カラー画像の階調方向判定を示す図
【図９】カラー画像ブロック分割領域と輝度勾配算出の代表点位置を示す図
【図１０】階調方向係数の補間を示す図
【図１１】本発明の実施の形態１における処理の流れを示す図
【図１２】２次元の場合における異なる階調方向で使用すべき補間方法を示す図
【図１３】本発明の実施の形態２におけるカラー画像処理装置のブロック構成図
【図１４】本発明の実施の形態３におけるカラー画像処理装置のブロック構成図
【図１５】本発明の実施の形態３におけるテーブル非線型性マップ作成のためのブロック構成図
【図１６】本発明の実施の形態４におけるカラー画像処理装置のブロック構成図
【図１７】従来の色変換エンジンのブロック構成図
【符号の説明】
１０１ 入力カラー画像
１０２ 画素入力部
１０３ ３次元色変換テーブル
１０４ 階調方向判定部
１０５ 階調方向係数補間部
１０６ 出力混合部
１０７ 色空間補間部
１０８ 出力カラー画像
１３０１ 非線形性判定部
１４０１ テーブル非線形性マップ
１６０１ ５種の色変換テーブル
【発明の属する技術分野】
本発明は、カラー画像処理技術一般、特にカラーマネジメント技術に関する。
【０００２】
【従来の技術】
近年、カラーマネジメント技術の分野では、スキャナ、モニタ、プリンタなどのデバイスの間のカラーマッチングを図るためコンピュータのOSレベルでソフトウエア的な仕組みが組み込まれるようになっている。 このCMS（カラーマネジメントシステム）の機能は、あらかじめ各カラーデバイス毎に「プロファイル」とよばれるデバイス色変換特性情報を用意しておき、色変換エンジンがカラー画像を前記プロファイルに基づき色変換実行することにより行われる。プロファイルと色変換エンジンの構成として最近使われているのが多次元補間法である。
【０００３】
この方法は色変換における入出力関係をＭ次元の入力からN次元の出力空間への数学的変換と見なし、色変換を計算でなく色変換テーブルのルックアップで実行するものである。現在使用されている多くのカラー画像の場合、画素値が（RGB）あるいは（CIE-LAB）あるいは（CMYK）など３変数、４変数にて表現されるためＭ、Nは、３あるいは４である。しかし将来使用が想定されるカラーマルチスペクトル画像の場合にはＭ，Ｎは１０次元以上の可能性もある。ここでは説明の都合上、色変換を、３入力１出力の変換と考えることとする。テーブル法は、３次元の場合各色信号が８ビットあるとすると、そのままでは1600万色にもおよぶフルカラーの色の各々に対する膨大なサイズのテーブルを構成する必要がある。そこで例えば入力色が（RGB）空間内にある場合、まず、入力色空間を格子に区切り、格子点上での出力値を３次元テーブルf(i,j,k)に蓄積しこれを「プロファイル」とする。
【０００４】
実際の色変換に際しては、色変換エンジンが画像の各画素のカラー値（RGB）を前記３次元テーブルの３次元空間内で近傍格子点での出力値を探し、それらを重み付け加算する補間方法にて出力を計算する。この手法はテーブルサイズを大きくすることにより変換精度がいくらでも高められること、変換速度が高速でかつ色変換内容によらず一定であること、などの理由で広く使われるに至った。図１７に、この色変換エンジンの典型的な構成図を示す。この図では、色空間の３次元補間に四面体（三角錐）を用いているが、同図の下部に示すように立方体（８点）、三角柱、斜三角柱（６点）、ピラミッド（５点）など様々な３次元補間方法が従来から開発されている。
【０００５】
さて、３次元補間方法では補間誤差の発生は避けられない。特に色の値同志の絶対的な誤差よりも、画像の滑らかな階調を色変換した結果に生じる補間結果の暴れ現象（リップル）が人間の視覚上非常に影響が大きい。このリップルについては、例えば、本出願人による文献１：金森：「３Dテーブルを用いた色変換における補間誤差解析」（カラーフォーラム JAPAN’98）に解説されている。グレイの階調上ではおなじガンマ変換を立方体を用いる３次元補間する場合でも色空間全体を定義域とする大域的（グローバル）ガンマ変換では変換後の階調カーブが折れ線近似される補間特性を示すが色空間内のグレイ階調部分だけを局所的にガンマ変換する場合には補間結果で階調が反転する程の補間曲線の暴れ現象を呈し（リップル）しかもテーブルサイズを大きくしてもリップル振幅は変化しないという特異現象も見つかっている。
【０００６】
この例からもわかるように一般に線形の多次元補間方法では、多くの格子点情報を用いる必要があり、それが画像の２次元平面内の滑らかな階調上に隠れた高次元特性を顕現する結果リップルが出現し階調特性が悪化することが多かった。また、リップル現象は、現在３次元補間において既に発生しているが、より高次元の４次元色空間やマルチスペクトル画像を扱う際の１０次元色空間になれば一層深刻な問題になることは容易に想像できる。
【０００７】
一方、文献２：James M.Kasson：“Performing color space conversions with three-dimentional linear interpolation”、Journal of Electronic imaging ／Jul 1995 Vol.4(3) pp.226-249)におけるディスフェノイド補間や特開平８−８３３３０号公報での体心構造を用いた四面体補間方法では四面体群を同一形状としたこと、最も補間誤差の大きい立方体の体心に格子点を位置させること、格子点数を各軸２倍に増加すると全体で８倍になるところを２倍の増加で済むことなどの特徴を持っており、良好な補間特性を持つとされている。特に、四面体稜線が色空間内で多くの方向に発生するため、画像の入力方向に対して稜線方向が合致する度合いが増え、先に述べたように大きなリップルを発生することが少ない効果がある。
【０００８】
しかし、体心構造では、RGB立方体の外部に不必要な外挿点を持つ必要があること、空間分割が複雑で処理時間がかかる、などの欠点がある。また、稜線方向がいかに増えても、それをはずれる入力は常に作り得るため問題の完全な解決にはなっていない。また３次元においても複雑な空間分割手法を使用する必要がありマルチスペクトルを用いた高次元でのカラー処理においては計算量が多大なものとなりやすい課題もある。
【０００９】
また、従来、特開平６−３１１３５５号公報では、入力色によってトライ・リニア補間と四面体補間を切り替えたり、特開平７−９５４２３号公報ではトライ・リニア補間、三角柱補間、ピラミッド補間を切り替えたり、色空間内で入力する色の位置によって異なる補間方法を採用することによって全体での色変換誤差を小さくする方法もある。
【００１０】
しかし、これらの方法は、入力した色の空間内の位置により補間方法が決定されるため種々の非線形性色変換において汎用的に対処できるものではない。
【００１１】
【発明が解決しようとする課題】
以上のように従来の多次元色空間補間法を用いた色変換においては、様々な色変換、様々な画像入力に対応して常に階調性が良好な補間方法を実行することが要求されているが、あらゆる画像、あるいは今後予想される多次元カラー処理において充分対応可能な解決策がなかった。
【００１２】
本発明は、上記の課題を解決するもので、カラーマネジメント分野でカラー画像を多次元補間方法で色変換する場合において、色変換自身が非線形性をもつ場合には、カラー画像上の色空間内における階調（グラデーション）方向を画像の画面内の部分領域毎に検出し、該部分領域ごとに補間方法を逐次制御すること、より詳細には色空間内の補間において入力点を囲む全ての格子点を使うのではなく画像の階調方向に存在する少数の格子点を選択して補間することにより多次元補間における非線形性を下げ、リップルを生じない階調性を重視した色変換を実施することができるカラー画像処理方法を提供することを目的とする。
【００１３】
【課題を解決するための手段】
この課題を解決するため本発明は、カラー画像を走査し各画素値を入力する画素入力部と、前記画素を含む画像上の部分領域毎で、各色版毎に輝度勾配の大きさを判定し当該階調の方向性を色空間内で離散化した方向軸ごとの方向重みを計算する階調方向判定部と、前記階調方向判定部で計算した方向重みを画像面内で連続的に補間し階調補間係数を出力する階調方向係数補間部と、色変換テーブルを用いて前記画素値に対して色変換を実施する複数の相異なる色空間補間部と、前記複数の色空間補間部からの出力を前記階調方向係数補間部から出力される階調方向係数により混合あるいは選択する出力制御部とを備える。
【００１４】
これにより、いかなる色変換、いかなる画像入力についても画像全体でリップルを生じない階調性を重視した色変換を実施することができるものである。
【００１５】
本発明の第１の実施の形態は、カラー画像を走査し各画素値を入力し、前記画素値の色空間内の位置情報から色変換テーブルを参照し多次元補間方法を用いて色変換する方法において、前記画素値を含むカラー画像を複数の分割領域内で各色版毎に輝度勾配の大きさを判定し、当該階調の方向について色空間内で離散化した方向軸ごとの方向重み係数を計算することにより、画像の階調方向情報を求め、前記階調方向情報を用いて前記多次元補間方法の制御することを特徴とするカラー画像処理方法であり、複数の補間方法から階調方向に合致した補間方法を特に重視して出力するため不必要に多くの格子点を使うために発生していた非線形性を低下させ、出力階調にリップルなどのアーティファクトが低減されるという作用を有する。また、画像の階調の方向を色空間内の方向として表現できる。
【００１７】
好ましくは、カラー画像を分割領域内での大評点近傍の階調の離散化した方向軸への方向重みを前記分割領域間で連続的に補間することにより、複数の異なる多次元補間方法を画像の階調方向によって画像面内で逐次切り替えても画像の階調の連続性が崩れないという作用を有する。
【００１８】
本発明の第２の実施の形態は、カラー画像を走査し各画素値を入力する画素入力部と、前記画素を含む画像上の部分領域毎で、各色版毎に輝度勾配の大きさを判定し当該階調の方向性を色空間内で離散化した方向軸ごとの方向重みを計算する階調方向判定部と、前記階調方向判定部で計算した方向重みを画像面内で連続的に補間し階調補間係数を出力する階調方向係数補間部と、色変換テーブルを用いて前記画素値に対して色変換を実施する複数の相異なる色空間補間部と、前記複数の色空間補間部からの出力を前記階調方向係数補間部から出力される階調方向係数により混合あるいは選択する出力制御部とからなるものであり、複数の補間方法から階調方向に合致した補間方法を特に重視して出力するため出力階調においてリップルが削減されるという作用を有する。
【００１９】
本発明の第３の実施の形態は、カラー画像を走査し各画素値を入力する画素入力部と、前記画素を含む画像上の部分領域毎で、各色版毎に輝度勾配の大きさを判定し当該階調の方向性を色空間内で離散化した方向軸ごとの方向重みを計算する階調方向判定部と、前記階調方向判定部で計算した方向重みを画像面内で連続的に補間し、階調補間係数を出力する階調方向係数補間部と、色変換テーブルを用いて前記画素値に対して色変換を実施する複数の相異なる色空間補間部と、前記複数の色空間補間部からの補間出力を比較する出力比較部と、前記複数の色空間補間部からの補間出力を出力比較部での判定結果に基づき前記階調方向係数により混合あるいは選択する出力制御部とからなるものであり、色変換自体に非線形性がある場合に限り複数の補間方法から階調方向に合致した補間方法を特に重視して出力するため、処理速度を落とさずに出力階調にリップルが削減されるという作用を有する。
【００２０】
本発明の第４の実施の形態は、カラー画像を走査し各画素値を入力する画素入力部と、前記画素を含む画像上の部分領域毎で、各色版毎に輝度勾配の大きさを判定し当該階調の方向性を色空間内で離散化した方向軸ごとの方向重みを計算する調方向判定部と、前記方向判定部で計算した方向重みを画像面内で連続的に補間し階調補間係数を出力する階調方向係数補間部と、色変換テーブルを用いて前記画素値に対して色変換を実施する複数の相異なる色空間補間部と、前記色変換テーブルのうち非線形性が強い領域を記述するテーブル非線形性マップと、複数の色空間補間部からの出力を前記テーブル非線形性マップに基づき前記階調方向係数により混合あるいは選択する出力制御部とからなるものであり、色変換自体に非線形性がある場合に限り複数の補間方法から階調方向に合致した補間方法を特に重視して出力するため、処理速度を落とさずに出力階調にリップルが削減されるという作用を有する。
【００２２】
好ましくは、多次元テーブルに対してあらかじめ用意したテスト信号を入力して異なる多次元補間を実施し、それら出力のばらつきに基づいて色変換テーブル内の非線形強度の大きい部分を判定し、ファイル化するものであり、あらかじめ多次元色変換テーブル内で非線型性強度の強い部分を判定し対処することができるという作用を有する。
【００２３】
以下、本発明の実施の形態について図１から図１５を用いて説明する。
【００２４】
(実施の形態１)
図１に本発明の実施の形態１におけるカラー画像処理装置のブロック構成図を示し説明する。本実施例において、カラー入力画像１０１として想定しているものは、１画素が（R,G,B）の３変数信号で規定される画像であり、３種の画像としてR版、G版、B版と呼ぶものとする。
【００２５】
図１において、１０２は、カラー画像１０１を走査し各画素値を入力する画素入力部、１０４は画素を含む画像上の部分領域毎で階調の色空間内での方向性を判定し各領域の方向重みｕを求める階調方向判定部、１０５は階調方向判定部１０４での方向重みｕを画像面内で連続的に補間し階調方向係数Wi(m,n)を出力する階調方向係数補間部、１０７は多次元色変換テーブル１０３を用いて前記画素値に対して色変換を実施する複数の相異なる種類の色空間補間部、１０６は複数の色空間補間部１０７からの出力を階調方向係数補間部１０５から出力される階調方向係数Wi(m,n)によって混合あるいは選択し、出力カラー画像１０８を出力する出力制御部である。
【００２６】
上記のように構成されたカラー画像処理装置の動作を以下に説明する。
【００２７】
カラー入力画像１０１における注目画素の画面上のｘｙ座標を(m,n)、カラー値を(R,G,B)とする。画素入力部１０２は、注目画素のカラー値（R,G,B）を各信号を８ビット構成とした場合、上位信号３ビットの組（Rh,Gh,Bh）と下位信号５ビットの組（Rl,Gl,Bl）に分割し、それぞれ色変換テーブル１０３と色空間補間部１０７に供給する。
【００２８】
画素入力部１０２にて分割された上位信号（Rh,Gh,Bh）にて色変換テーブル１０３をルックアップし、下位信号（Rl,Gl,Bl）とで補間演算を実施する。色変換テーブル１０３には、上位信号（Rh,Gh,Bh）に対応して全RGB色空間を８×８×８＝２５６に分割した頂点である９×９×９個の格子点での色変換出力が格納されている。ここで用いた３ビットと５ビットというビット分割方法は色変換テーブル１０３のサイズ、すなわち格子点数（格子点間距離：ここでは２5＝３２）を決定しているが実際には任意である。色変換テーブル１０３は、色再現を実行する各種色補正を目的とする数式（モデル式）にて表現できる色変換を表現できるのみならず、カラー機器の測色によって数式モデルを持たない自由度の高いテーブルを作成することも可能である。
【００２９】
また、色変換の入出力の定義も様々であり、従来は積分濃度空間（Dr,DgDb）からインク（トナー）信号（C,M,Y,K）へのプリンタ用の色変換が一般的であったが、カラーマネジメントの世界では(R,G,B)から標準色空間（L,A,B）（X,Y,Z）への変換、あるいは逆変換などが対象とされる。本発明に関しては（R,G,B）から任意のある出力色ｆへの色変換を考えることとする。すなわち色変換テーブルには出力１色分についてはｆ(R,G,B)という３変数関数が格納されていることになる。
【００３０】
さて、画素入力部１０２からの上位信号（Rh,Gh,Bh）は、色変換テーブル１０３の８×８×８の分割された立方体の１つを指定する。つまり上位信号によって入力色を内包する最近傍点８個が見つかり、そこでの出力ｆが得られる。この８点を用いて補間する方法がもっとも一般的な３次元補間方法であるトライ・リニア補間である。
【００３１】
色空間補間部１０７のうち、（０）は、このトライ・リニア補間方法を実行する部分であり以下その機能について図２を用いて簡単に説明する。得られた８点で構成される立方体内で入力信号を表現するのは下位信号（Rl,Gl,Bl）であり、３次元補間処理では、これらを重みとして用いる。ここで図２に示すように入力色をPとし、Pを含む色空間内の立方体をABCD-EFGHとし、Ｐの位置を（Ｒh＋Ｒl，Ｇh＋Ｂl，Ｂh＋Ｂl）とすると、
【００３２】
【数１】
であるから、Ｐでの補間出力ｆIが以下の（数２）で表現できる。
【００３４】
【数２】
ただし、その際の重み係数は（表１）の通りである。
【００３６】
【表１】
トライ・リニア補間は、公知技術として従来から使用されている。しかし、８点を用いる補間法であるため一般には精度が高いといわれながら、非線型色変換では補間特性にリップルを生ずることが知られている。この原因は表１の各重み係数が３次式となっていることに関係がある。この課題の解決のため従来から利用されてきた方法が四面体補間である。四面体補間とは、色空間で入力色を含む４点を用いた線形補間の総称であるが、最も一般的なのは前記の立方体ABCD-EFGHを対角線を境界の陵線として６個の体積が等しい四面体分割を行い、いずれの四面体領域に含まれるかを判定した上で当該の４点を用いて補間するダイアゴナル・テトラヘドロン補間方法である。ダイアゴナル・テトラヘドロン補間では補間立体内では線形補間が行われるためリップルが発生しない上、用いる点数が３次元補間では最小数である４点であるため計算量も少ないという良好な性質を持っている。しかし３次元空間で等方的な８点補間と異なり、上記の良好な性質は方向性により特性が異なる。そこで本発明ではダイアゴナル・テトラヘドロン補間について方向性を考えて４種類を用いている。つまり、図１の色空間補間部１０７の（１）から色空間補間部（４）までは、全てダイアゴナル・テトラヘドロン補間を実行するが、図３、および（表２）の(D1)、（D2）,（D3）,(D4)、に示すようにダイアゴナルテトラヘドロンの分割方法である対角線方向が４種類に異なっている。
【００３８】
【表２】
ここでは（D１）タイプ則ち立方体ABCD-EFGHのA-G対角線に沿った分割を例として補間方法を説明することとするが、他の３タイプの場合にも説明は同様である。
【００４０】
図４は、ダイアゴナル・テトラヘドロン補間方法（D１）の６分割方法を説明する図である。（表３）に従い、下位信号（Rl,Gl,Bl）の大小関係は６通り存在し、これが６分割四面体のいずれかを決定する。Rl=Gl=Blの時にはArea０からArea５までのいずれに含ませても良いがここではArea０とした。
【００４１】
【表３】
分割された四面体が決定した後に、以下のように４点補間が実行される。
【００４３】
【数３】
この際の重み係数は以下の（表４）のように決定する。
【００４５】
【表４】
次に、階調方向判定部１０４、階調方向係数補間部１０５および出力制御部１０６について説明する
【００４７】
階調方向判定部１０４にて入力画像の階調の方向を判定し、方向重みｕを求める。階調方向係数補間部１０５は、階調方向判定部からの方向重みｕを画像画面内で連続的に補間し、階調補間係数Ｗi(m,n)を出力する。
【００４８】
出力制御部１０６において当該方向に近い方向の分割線を持つ補間方法に大きな階調方向係数を与え、相異なる方向の分割線を持つ補間方法には小さな階調方向係数を与えて出力結果を加算する。出力制御部１０６では、色空間補間部（０）から（４）までの出力をｆI(i)(i=0-4)と表現し、階調方向係数Wiを用いて混合する。
【００４９】
【数４】
この階調方向係数Wiは、各補間方法ごとに与えられるものであり、
【００５１】
【数５】
を満たすものとする。
【００５３】
次に、階調方向係数Wiの計算の仕方について、図８を用いて説明する。まず、カラー画像を、Red,Green,Blueの３つの色版からなる単色版の画像と考え、各々の版毎に局所的な画像の階調の大きさを判定する。はじめに図９に示すように、画像を部分領域に分割しておき、これから色変換を行う対象画素位置(m,n)を含む領域内で代表点（m0,n0）を決定し、その場所で階調方向を算出する。すなわち画像に対して輝度勾配ベクトルDを算出するフィルタリングを実施すればよく、ここでは偏微分を用いるが現実にはこの操作はデジタル的な差分操作となる。
【００５４】
【数６】
この階調方向ベクトルＤは、階調の急激さとは無関係な量にするため長さ＝１に正規化するのが望ましい。次に階調方向ベクトルＤの方向性を調べる。これは簡易的には、７個の方向ベクトルVjと内積をとることで実現できる。７個の方向ベクトルVとは図８で示すV1からV7の７方向であり、各ベクトルのRGB座標成分は以下の（表５）の通りである。
【００５６】
【数７】
【表５】
以上求めた各成分ごとの７個の重み係数ujについて、u1からu4までは色空間補間部（１）から（４）までの重みとしてそのまま用いられ、u5からu7までは加算され色空間補間部（０）の重みとして用いられる。
【００５９】
【数８】
そして、W0とW1、W2, W3, W4の階調方向係数Wiが領域Bnごとに決定される。
【００６１】
この値をそのまま出力し、分割領域ごとに用いても良いが、階調方向係数Wiが領域間で大きく変化した場合、分割領域ごとに色や階調の段差が発生する可能性があり、これを防ぐため階調方向係数Wiが画像内で連続性を保つように、階調方向係数補間部１０５で画面内のバイ・リニア補間によって階調方向係数Wi(m,n)を補間出力している。
【００６２】
【数９】
【数１０】
ただし、階調方向係数の重み係数ωP、ωQ、ωR、ωSは、以下の（表６）に示す通りである。ここでBX,BYは分割領域の幅、高さを示す。また上記の内挿補間は、画像内で各Wiごとに実施するものとする。
【００６５】
【表６】
図９、図１０に示すようにこの方法では、重み係数は画像の周辺に内挿不可能な領域を残すが、ここは外挿あるいは分割領域での代表点(m0,n0)の値をそのまま使用することも考えられる。また代表点(m0,n0)を各分割領域の左上にとり、画像の右端では仮想的な固定値を採用することにより外挿を不必要にすることも可能である。
【００６７】
次に、実際の画像処理の流れについて図１１を用いて説明する。この例では、カラー画像は入力出力とも（RGB）点順次画像を想定しており、色変換処理も点順次に実行されていくこととする。
【００６８】
処理開始後、Ｓ１１０１のステップにおいて、カラー画像は図９のように各領域に分割される。Ｓ１１０２のステップで領域内の代表点(m0,n0)における階調方向判定が行われ、各分割領域ごとの階調方向係数Wiを求める。ステップＳ１１０３において色変換に用いる３次元テーブルを読み込む。ここまでは画像処理の前処理として実行される。
【００６９】
次に、カラー画像の各画素を走査しながら、全画素の処理が終了するまでステップＳ１１０６以降ステップＳ１１１０までの処理を各画素ごとに繰り返す。最初にステップＳ１１０６において、画素の位置(m,n)および画素値（RGB）が入力される。ステップS１１０７において、当該画素位置(m,n)での階調方向係数Wi(m,n)が、前記の各分割領域ごとの階調方向係数Wiを補間して求められる。
【００７０】
次に、ステップＳ１１０８において、異なる色空間補間処理を並列に実行して補間出力fI(0)からfI(4)を求める。ステップＳ１１０９においてＳ１１０７で求めた階調方向係数Wi(m,n)を用いて補間出力が混合される。ステップ１１１０において出力画素(R’G’B’)が出力画像に挿入される。
【００７１】
処理は、そのまま繰り返しのためステップＳ１１０４へ戻り、全画素の処理が終了した場合にはステップＳ１１０５の判定により処理が終了する。
【００７２】
なお、以上の説明では、色空間補間部１０７をトライリニア補間とダイアゴナル・テトラヘドロン補間４種類ということで説明したが、方向性を持つ分割空間であれば他の補間方法での組み合わせでも同様に実施可能であり種類も４種類には限らないことはもちろんである。
【００７３】
以上まで、本発明の構成および動作を説明してきたが、ここで階調方向によって補間内容を変更するという本発明の原理の説明を行う。まず何故トライリニア補間とダイアゴナルテトラヘドロン補間（しかも４種）を併用しているのか、さらに画像の階調方向を考慮してそれら重み付けしているのかを説明する。
【００７４】
従来、トライ・リニア補間とダイアゴナル・テトラヘドロン補間のどちらが高精度な補間方法であるかについては用いる色空間、用いる色変換が異なった場合、明解な説明がなかった。しかし、前述の文献１に解説されているように用いる色変換自身が１次元変換の線形結合にて表現される場合には補間結果は全ての線形補間で同一となり、出力結果が各成分ごとに折れ線状態になることが知られている。つまり、この場合にはトライ・リニア補間とダイアゴナル・テトラヘドロン補間両者の結果の差はなくリップルも存在しない。
【００７５】
このことを２次元の補間を例にとって図形的に説明する。はじめに２次元では、トライリニア補間は４点を用いるバイリニア補間に対応し、ダイアゴナル・テトラヘドロン補間は、正方形を２分割した三角形の３点補間に相当することに注意する。図５の５００において正方形ＡＢＣＤはＲＧ平面に設定された１つの補間区間であり、出力１成分をｆとしてこの２次元平面上に立てた矢印にて表現する。色変換が線形であるとは、この正方形内ではｆ（Ａ）、ｆ（Ｂ）、ｆ（Ｃ）、ｆ（Ｄ）が同一平面π上に存在することに相当する。従ってこの４点からの線形補間であるバイリニア補間による補間値は５０１に示すように同一平面π上に拘束されており、正方形ＡＢＣＤ内をいかなる軌跡で入力された色についても出力は直線となる。
【００７６】
次に、図５の５０２，５０３は、A-C分割、B-D分割という異なる種類の三角形分割における補間結果をしめす図である。いずれの分割方法の場合でも平面内の入力点が２種の三角形のいずれに属するかを判定した上で各三角形内で補間される。三角形内での補間値はこの平面πを三角形ABC、または三角形ADCに分割して補間することに相当するが、元々平面を分割するのであるからA-C分割、B-D分割いずれの三角形分割でもうバイリニア補間と同一の補間結果となることは明らかである。このように色変換自身が非線形性を持たない場合には、補間方法の如何にかかわらず、結果が理想値からずれたり出力が曲線になる現象（リップル）が出現することもない。
【００７７】
次に、特開平６−３１１３５５号公報にてダイアゴナル・テトラヘドロン補間がトライリニア補間よりも良好な結果を得たと言われている色変換がMIN演算の場合、特にダイアゴナル・テトラヘドロン補間の中でも原点を通る分割線を採用する（D1）型の場合を取り上げる。同じように２次元で考えると、MIN演算とは、２変数の組(R,G)の小さい方を出力するものであるから図６の６００に示すようにR≧Gの領域に相当する三角形ABC内では、MIN(R,G)＝Gであり、R<Gの領域に相当する三角形ADC内ではMIN(R,G)=Rであるから直線R=Gを稜線として２枚の三角形平面が折れて接続されているような平面を構成する。これがMIN演算の場合における真の計算値である。
【００７８】
次に、これを補間で演算することを考える。まず６０１に示すようにABCDの４点補間法では補間面が曲面を形成してしまいこれがリップルを発生する原因になる。次に、図６の６０２に示すA-C分割での三角形分割では、稜線がA-C方向に形成されるため補間結果は正しい結果と完全に同一になる。しかし、６０３のB-D分割ではB-D方向で分割された１つの三角形では値が０に固定されたままとなり、２つの分割方法で補間出力が大きく異ってしまう。
【００７９】
前記特開平６−３１１３５５号公報では、この場合を取り上げてA-C分割の三角形補間がABCD４点より良好な結果を得るとしている。しかし正確に言えばＭＩＮ演算の真値の演算において、A-Cの分割線が全く違う特性を示す２種の平面が接続された稜線となっており、この稜線と三角形の分割境界とが一致しているため補間誤差が０になり良好な補間が行われるというのが正しい。結局、非線形性の強い色変換であってもMIN演算のように真の計算値が知られている場合には最適な補間方法が何らかの指標によって定まる点が重要である。逆にいえば非線形性が強い場合において真の計算値が不明な場合にはいかなる分割・補間をすれば良いかを知る指標が無くなるため即座に困難に直面することになる。
【００８０】
この困難を生む典型的な例が、前述の文献１で取り上げた局所ガンマ変換のような非線形変換である。この変換では、非線形性の強い色変換テーブルだけが与えられるので、色変換に当たってはこれを３次元補間する必要がある。すなわち先の例では、ｆ（A）、ｆ（B）、ｆ（C）、ｆ（D）という格子点上での値だけが非線形性が強い状態で与えられ、その間を接続する真の計算値は全く与えられない。非線形性が強いという意味は２次元では図７の７００に示すようにABCDの各点で対角線位置にあるA，Cで小さな値、B,Dでは大きな値をとるような場合である。これをバイリニア４点補間すると７０１に示すようにｆ軸Ｒ軸Ｇ軸で作られる空間において複雑な２次曲面を形成し、リップルの原因となる。
【００８１】
また、３点補間では７０２と７０３に示すようにA-C分割とB-D分割では、前者がA-C方向が深い谷を形成するのに対して後者はB-Dが稜線を形成して全く異なる補間結果となるが、いずれもリップルの代わりに激しい凹凸が発生する原因となる。つまり４点補間でも２種類の３点補間でも補間結果は大きく異なる上に、どれも良好な補間とはいえないという困難な状況となる。
【００８２】
本発明の基本的な考え方は、ここに入力画像の階調性の方向の情報を入れて最適な補間を決定するものである。たとえば図７の７００において鎖線(1)(2)が入力階調を示すものとすると、(1)のようなA-C方向の階調が入力された場合にはA-C分割の３点補間を採用し、(2)のようなB-D方向の階調が入力された場合にはB-D分割の３点補間を採用するものということができる。
【００８３】
次に、入力階調方向が求まった後、いかなる補間方法を使用するべきか、おなじく２次元の場合を例に図１２を用いて説明する。図１２では、入力階調として正方形ＡＢＣＤの中心点を通る４種の方向を設定し、入力点を黒丸で、その点での補間出力値を黒丸上の矢印で表示している。１２００と１２０１は入力階調方向がそれぞれＥＦ方向、ＧＨ方向の場合であり、このときは両方とも同一のバイリニア補間を使用すればよい。
【００８４】
図７の７０１では、複雑な２次曲面を呈しているかのように見えるバイリニア補間面であるが、１２００で示すように入力方向がＥＦ方向のときは補間結果がＡＢの中点ＥとＣＤの中点Ｆ上での値ｆ（Ｅ）とｆ（Ｆ）を結ぶ直線となりリップルが現れない性質がある。同様に１２０１の入力階調方向がＧＨの時もバイリニア補間の結果、補間値がｆ（Ｇ）とｆ（Ｈ）を結ぶ直線上に位置する。
【００８５】
次に１２０３では入力階調がＡＣ方向の場合なので、ＡＣ分割の三角形補間を用い、１２０４では入力階調方向がＢＤなのでＢＤ分割の三角形補間を用いることにより、いずれの場合も補間結果がそれぞれｆ（Ａ）とｆ（Ｃ）、あるいはｆ（Ｂ）とｆ（Ｄ）を結ぶ直線上に位置するためリップルは現れない。以上のように２次元の場合には、階調方向４種類（１）ＥＦ、（２）ＧＨ、（３）ＡＣ、（４）ＢＤについて（１）（２）はバイリニア補間、（３）と（４）は異なる三角形補間を使用すればよい。もちろん画像上の各画素で異なる補間方法を次々に変更していくと画像上にノイズを生む原因となる。このため補間方法の制御としては、最適な補間を１種選択するよりは重みとして用いるほが望ましい。つまり画像の階調方向から考えて最適な補間方法の結果には大きな重みを与え、そうでない補間結果には小さい重みを与えるように何種類からの補間方法に対して一種の重みづけを行う。
【００８６】
以上の議論を３次元に拡張すると、図８の方向ベクトルでＶ５とＶ６とＶ７のような立方体軸に平行の入力階調の場合にはいずれもトライリニア補間を用いるのがよく、その他の対角方向の入力Ｖ１、Ｖ２、Ｖ３、Ｖ４では、各々図３のＤ１，Ｄ２、Ｄ３、Ｄ４方法の分割のダイアゴナルテトラテドロン補間を用いるのが良いという結論が導かれる。このため、（数８）のように方向重みｕ5、ｕ6、ｕ7、ｕ8についてはそれらを加算したものをトライリニア補間用の階調方向係数Ｗ0として用い、他のＷ1、Ｗ2、Ｗ3、Ｗ4はｕ1、ｕ2、ｕ3、ｕ4をそのまま用いている。
【００８７】
以上のように、本発明の第１の実施形態では、いかなる非線形色変換においても、入力画像の階調を求め、数種類の補間方法を実施してその結果に重み付けすることで出力階調でのリップルを抑制するものである。
【００８８】
（実施の形態２）
図１３に本発明の実施の形態２の画像処理装置のブロック図を示し説明する。
【００８９】
実施の形態１では、異なる色空間補間部１０７からの出力は入力画像１０１の画素位置から計算される階調方向係数によって画素ごとに常に重み付け混合されていたが、この重みづけ混合の計算は画像の全画素について実行されるため計算処理負担が大きい。また実際には非線型の高い色変換やそれによる画像の階調の劣化は常に発生するわけではない。
【００９０】
そこで、本実施の形態では、５種の色空間補間部(0)〜(4)での出力を非線形性判定部１３０１にて比較しこれらが大きく異なる場合、すなわち非線形性が高い場合にのみ、重み付け混合を実施することによって、処理の負担を下げるものである。非線形判定部１３０１では、以下のような演算を行う。ここでTHRESHは固定の閾値、ｆI（i）、ｆI（ｊ）はｆI(0)から ｆI(4)のいずれかを示す。MAXは最大値演算を示す。
【００９１】
つまり、ある入力色について異なる分割方向や異なる種類の補間を用いたとき最大でどの程度の変動があるかを調査し、ある閾値以上であったら非線形性のが強いと判断し、結果resultをYESと判定し、そのほかはNOと判定し、この結果を一緒に出力制御部１０６に送ることにより出力混合を行うのかどうかを決定する。
【００９２】
【数１１】
この方法では、非線形性判定部１３０１がYES/NOの判定を行う境界が生じるため出力の切り替えによって階調飛びが生じるように思われるるが、閾値を人間の視覚が階調とびを認識できる限界値に選んでおけば補間方法による差が閾値を超える前までは出力混合を停止して、超えた時点で出力混合を開始してもそのギャップは前記閾値よりは小さいことが保証されるため認識できる大きな階調飛びは発生しない。
【００９４】
（実施の形態３）
図１４に本発明の実施の形態３における画像処理装置のブロック図を示し説明する。実施の形態３は、実施の形態２における非線形性判定部を無くし、その代わりにあらかじめ用意したテーブル非線形性マップ１４０１の情報を元に、入力画素の色信号値が色変換テーブル内での非線形性の大きい部分に入った場合にのみ、重み付け混合を実施することによって処理の負担を下げるものである。この場合、色変換テーブルを作成した段階でテーブル非線形性マップが一緒に作られ、この情報はICCプロファイルのようなカラープロファイルの１つの情報として蓄積しておくこともできる。
【００９５】
次に、図１５を用いてテーブル非線形性マップの作成方法を説明する。図１５はテーブル非線形性マップ作装置を示すブロック図である。色変換のための３次元テーブルに対応して入力格子点の色空間内の位置において色変換の非線形性、主として分割方向による補間結果の差の大きさを判定し当該色での補間出力の混合が必要かどうかを判定するものである。作成方法は、第２の実施例に似た構成で、対象となる色変換テーブルを用いてテスト用カラー画像を色変換する。異なる色空間補間部で出力された補間出力どうしを非線形性判別部１３０１で（数１１）と同一の判定を行い、（Rh,Gh,Bh）をアドレス入力とした時の非線形性の結果として、判定結果YES/NOをテーブル非線形性マップ１５０１に格納する。
【００９６】
（実施の形態４）
図１６に本発明の実施の形態４における画像処理装置のブロック図を示し説明する。実施の形態４は、実施の形態１における３次元色変換テーブルを異なる複数の色空間補間部１０７（０）〜（４）ごとに専用の色変換テーブル１６０１としたものである。本発明の実施の形態１から３では、１つの３次元色変換テーブルに対して画像の階調方向によって方向性をもつ色空間補間部１０７を使いわけることに利点を有するが、本実施の形態は３次元色変換テーブル自体をも使い分けることにより、画像の階調方向性の階調性を一層良好なものとすることができる。
【００９７】
尚、本発明は３次元のカラー画像を対象にして実施の形態を書いているが入力空間は一般に４次元以上のＭ次元でも同一の考え方でかまわない。出力空間がＮ次元の場合にはＭ入力１出力の個々の色変換部をＮ個並べるだけで実現できる。
【００９８】
従って本発明は、Ｍ入力Ｎ出力のマルチスペクトルカラー画像の色変換においても有効である。
【００９９】
なお、本発明は、階調判定の単位を画面上のブロック分割領域としたが、画像内の別のセグメンテーション方法によるものであってもよい。
【０１００】
【発明の効果】
以上のように本発明によれば、従来技術が色変換を各画素値のみを用いて行っていたため生ずる階調の乱れの問題に対して、カラー画像上の階調（グラデーション）方向を画面内の部分領域毎に検出し、該部分領域ごとに補間方法を逐次制御することにより、いかなる色変換、いかなる画像入力についても画像全体でリップルを生じない階調性を重視した色変換を実施することができるカラー画像処理装置を提供することができる。
多次元カラーマルチスペクトル画像の補間方法においては、入力次元が増大するため高次項が補間結果により深刻な影響を与え未知のリップルが発生しやくなることが考えられるが本発明の画像の階調方向を考慮するという考え方を利用すれば現実的に良好な補間演算が可能になる。
【図面の簡単な説明】
【図１】本発明の実施の形態１によるカラー画像処理装置のブロック構成図
【図２】トライ・リニア補間方法を示す図
【図３】４種の四面体分割方法を示す図
【図４】 D1分割方法でのダイアゴナル・テトラヘドロン補間方法を示す図
【図５】２次元の場合における非線形性が無い色変換補間を示す図
【図６】２次元の場合におけるMIN演算の色変換補間を示す図
【図７】２次元の場合における非線形性が強い色変換補間を示す図
【図８】カラー画像の階調方向判定を示す図
【図９】カラー画像ブロック分割領域と輝度勾配算出の代表点位置を示す図
【図１０】階調方向係数の補間を示す図
【図１１】本発明の実施の形態１における処理の流れを示す図
【図１２】２次元の場合における異なる階調方向で使用すべき補間方法を示す図
【図１３】本発明の実施の形態２におけるカラー画像処理装置のブロック構成図
【図１４】本発明の実施の形態３におけるカラー画像処理装置のブロック構成図
【図１５】本発明の実施の形態３におけるテーブル非線型性マップ作成のためのブロック構成図
【図１６】本発明の実施の形態４におけるカラー画像処理装置のブロック構成図
【図１７】従来の色変換エンジンのブロック構成図
【符号の説明】
１０１ 入力カラー画像
１０２ 画素入力部
１０３ ３次元色変換テーブル
１０４ 階調方向判定部
１０５ 階調方向係数補間部
１０６ 出力混合部
１０７ 色空間補間部
１０８ 出力カラー画像
１３０１ 非線形性判定部
１４０１ テーブル非線形性マップ
１６０１ ５種の色変換テーブル
Claims (6)
- カラー画像を走査し各画素値を入力し、前記画素値の色空間内の位置情報から色変換テーブルを参照し多次元補間方法を用いて色変換する方法において、前記画素値を含むカラー画像を複数の分割領域内で各色版毎に輝度勾配の大きさを判定し、当該階調の方向について色空間内で離散化した方向軸ごとの方向重み係数を計算することにより、画像の階調方向情報を求め、前記階調方向情報を用いて前記多次元補間方法の制御することを特徴とするカラー画像処理方法。
- 階調方向情報の算出は、前記方向重み係数を画像内の領域間で連続的に補間することを特徴とする請求項１記載のカラー画像処理方法。
- カラー画像を走査し各画素値を入力する画素入力部と、前記画素を含む画像上の部分領域毎で、各色版毎に輝度勾配の大きさを判定し当該階調の方向性を色空間内で離散化した方向軸ごとの方向重みを計算する階調方向判定部と、前記階調方向判定部で計算した方向重みを画像面内で連続的に補間し階調補間係数を出力する階調方向係数補間部と、色変換テーブルを用いて前記画素値に対して色変換を実施する複数の相異なる色空間補間部と、前記複数の色空間補間部からの出力を前記階調方向係数補間部から出力される階調方向係数により混合あるいは選択する出力制御部とからなることを特徴とするカラー画像処理装置。
- カラー画像を走査し各画素値を入力する画素入力部と、前記画素を含む画像上の部分領域毎で、各色版毎に輝度勾配の大きさを判定し当該階調の方向性を色空間内で離散化した方向軸ごとの方向重みを計算する階調方向判定部と、前記階調方向判定部で計算した方向重みを画像面内で連続的に補間し、階調補間係数を出力する階調方向係数補間部と、色変換テーブルを用いて前記画素値に対して色変換を実施する複数の相異なる色空間補間部と、前記複数の色空間補間部からの補間出力を比較する出力比較部と、前記複数の色空間補間部からの補間出力を出力比較部での判定結果に基づき前記階調方向係数により混合あるいは選択する出力制御部とからなることを特徴とするカラー画像処理装置。
- カラー画像を走査し各画素値を入力する画素入力部と、前記画素を含む画像上の部分領域毎で、各色版毎に輝度勾配の大きさを判定し当該階調の方向性を色空間内で離散化した方向軸ごとの方向重みを計算する調方向判定部と、前記方向判定部で計算した方向重みを画像面内で連続的に補間し階調補間係数を出力する階調方向係数補間部と、色変換テーブルを用いて前記画素値に対して色変換を実施する複数の相異なる色空間補間部と、前記色変換テーブルのうち非線形性が強い領域を記述するテーブル非線形性マップと、複数の色空間補間部からの出力を前記テーブル非線形性マップに基づき前記階調方向係数により混合あるいは選択する出力制御部からなることを特徴とするカラー画像処理装置。
- 多次元テーブルに対してあらかじめ用意したテスト信号を入力して異なる多次元補間を実施し、それら出力のばらつきに基づいて色変換テーブル内の非線形強度の大きい部分を判定し、ファイル化する請求項５記載のテーブル非線型性マップの作成方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP13821199A JP3835059B2 (ja) | 1999-05-19 | 1999-05-19 | カラー画像処理方法およびその装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP13821199A JP3835059B2 (ja) | 1999-05-19 | 1999-05-19 | カラー画像処理方法およびその装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2000333028A JP2000333028A (ja) | 2000-11-30 |
| JP3835059B2 true JP3835059B2 (ja) | 2006-10-18 |
Family
ID=15216688
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP13821199A Expired - Fee Related JP3835059B2 (ja) | 1999-05-19 | 1999-05-19 | カラー画像処理方法およびその装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3835059B2 (ja) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006154452A (ja) | 2004-11-30 | 2006-06-15 | Toshiba Corp | 映像信号処理装置とグラデーション段差検出方法 |
| JP4481881B2 (ja) * | 2005-06-07 | 2010-06-16 | キヤノン株式会社 | 色分解テーブル作成方法および画像処理装置 |
| JP2007134750A (ja) * | 2005-10-12 | 2007-05-31 | Seiko Epson Corp | 色変換装置、色変換方法、色変換プログラム、画像処理装置、及び画像表示装置 |
-
1999
- 1999-05-19 JP JP13821199A patent/JP3835059B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2000333028A (ja) | 2000-11-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4209599B2 (ja) | ２次元非線形補間システム及び補間方法、並びに２次元混合補間システム及び混合補間方法 | |
| US6058208A (en) | Color correction device and color correction method | |
| US6058248A (en) | Computerized method for improving data resolution | |
| US7046862B2 (en) | Image processing apparatus and program | |
| JP3664364B2 (ja) | 色変換方法 | |
| JPH0730772A (ja) | 画像処理装置 | |
| EP1221812B1 (en) | Fast interpolation of large color lookup tables | |
| US6898332B2 (en) | Image processing device and image processing method | |
| JP3830304B2 (ja) | 補間方法及び装置 | |
| US5442717A (en) | Sharpness processing apparatus | |
| JP3583924B2 (ja) | 解像度変換方法及び変換装置 | |
| JP3835059B2 (ja) | カラー画像処理方法およびその装置 | |
| JP4164944B2 (ja) | 色域圧縮装置及び色域圧縮方法 | |
| JP4265363B2 (ja) | 画像処理装置 | |
| US7127118B2 (en) | System and method for compressing image files while preserving visually significant aspects | |
| JP3045134B2 (ja) | 画像データ補間装置、画像データ補間方法および画像データ補間プログラムを記録した媒体 | |
| JPH0799587A (ja) | 色変換方法および色変換装置 | |
| JPH0814843B2 (ja) | 色変換方法および色変換装置 | |
| JPH0467670B2 (ja) | ||
| JP3777672B2 (ja) | 画像処理方法および画像処理装置 | |
| US20040091173A1 (en) | Method, apparatus and system for the spatial interpolation of color images and video sequences in real time | |
| JP3743036B2 (ja) | 色変換方法および色変換装置 | |
| JP2001525623A (ja) | 第１の解像度のラスタのデジタルデータを第２の解像度のデジタルデータに変換する方法 | |
| JPH0946542A (ja) | 補間演算方法およびデータ変換装置 | |
| JP3403912B2 (ja) | 色変換装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| RD01 | Notification of change of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7421 Effective date: 20050629 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20051130 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060214 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060306 |
|
| 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: 20060704 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20060717 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090804 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100804 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110804 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110804 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120804 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130804 Year of fee payment: 7 |
|
| LAPS | Cancellation because of no payment of annual fees |