JP4649781B2 - 画像処理方法および装置 - Google Patents
画像処理方法および装置 Download PDFInfo
- Publication number
- JP4649781B2 JP4649781B2 JP2001187146A JP2001187146A JP4649781B2 JP 4649781 B2 JP4649781 B2 JP 4649781B2 JP 2001187146 A JP2001187146 A JP 2001187146A JP 2001187146 A JP2001187146 A JP 2001187146A JP 4649781 B2 JP4649781 B2 JP 4649781B2
- Authority
- JP
- Japan
- Prior art keywords
- value
- pixel
- threshold
- image
- level
- 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
- 238000003672 processing method Methods 0.000 title claims description 14
- 238000006243 chemical reaction Methods 0.000 claims description 73
- 238000000034 method Methods 0.000 claims description 37
- 230000008569 process Effects 0.000 claims description 21
- 238000001914 filtration Methods 0.000 claims description 13
- 230000009466 transformation Effects 0.000 claims description 8
- 238000005286 illumination Methods 0.000 abstract description 77
- 230000006835 compression Effects 0.000 abstract description 25
- 238000007906 compression Methods 0.000 abstract description 25
- 238000010586 diagram Methods 0.000 description 22
- 230000008859 change Effects 0.000 description 21
- 238000007796 conventional method Methods 0.000 description 10
- 230000001965 increasing effect Effects 0.000 description 8
- 230000005540 biological transmission Effects 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- 230000001186 cumulative effect Effects 0.000 description 4
- 238000009499 grossing Methods 0.000 description 3
- 239000000470 constituent Substances 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000004069 differentiation Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N5/00—Details of television systems
- H04N5/14—Picture signal circuitry for video frequency region
- H04N5/20—Circuitry for controlling amplitude response
- H04N5/205—Circuitry for controlling amplitude response for correcting amplitude versus frequency characteristic
- H04N5/208—Circuitry for controlling amplitude response for correcting amplitude versus frequency characteristic for compensating for attenuation of high frequency components, e.g. crispening, aperture distortion correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/12—Edge-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10016—Video; Image sequence
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20004—Adaptive image processing
- G06T2207/20012—Locally adaptive
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20172—Image enhancement details
- G06T2207/20192—Edge enhancement; Edge preservation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Image Processing (AREA)
- Facsimile Image Signal Circuits (AREA)
- Picture Signal Circuits (AREA)
Description
【発明の属する技術分野】
本発明は、テレビジョン、ビデオテープレコーダー、スチルカメラ、ビデオカメラおよびプリンタなど、各種の画像の入出力装置に好適に利用可能なものであり、特に、入力された画像を、相対的によりダイナミックレンジの狭い画像装置において再現するための画像処理方法および装置に関する。
【0002】
【従来の技術】
従来、例えば画像の階調特性の変換のために、入力画像の各画素に対し、そのレベルを図18の実線で示すような入出力関係を持つ関数(以下、「レベル変換関数」と記す。)で変換する方法(以下、「レベル変換」と記す。)がある。図18において、横軸は入力画像の画素レベル(入力レベル)lを、縦軸はレベル変換処理による出力画像の画素レベル(出力レベル)T(l)を表す。Lmaxは、入出力画像の各画素が取り得る最大レベルを表す。レベル変換後の画像のコントラストは、レベル変換関数の傾きが大きいほど増加することになる。図18の例では、入力レベルlbを境にして高レベル側と入力レベルlsを境にして低レベル側とにおけるレベル変換関数を示す直線の傾きが、中間レベル(入力レベルls〜lb)の傾きに比べて小さくなっている。従って、図18に示した関数を用いたレベル変換では、高レベルおよび低レベルにおけるコントラストを犠牲にすることで、中間レベルのコントラストを増加させていることになる。
【0003】
レベル変換関数は、図18に示したものに限らず、例えば図19の実線で示したものを用いることもできる。図19に示したレベル変換関数は、入力レベルlkを境にして、高レベル側における直線の傾きが、低レベル、中間レベルでの傾きに比べて小さくなっている。従って、図19に示した関数を用いたレベル変換では、高レベルでのコントラストを犠牲にすることで、低レベル、中間レベルでのコントラストを増加させることができる。また、以下の(1)式のガンマ関数や(2)式のLOG関数等、図18および図19に示した関数に比べて、より連続的なレベル変換関数が用いられることもある。なお、(1)式におけるgは関数の傾きを調整するパラメータである。
【0004】
【数1】
【0005】
【数2】
【0006】
また、さらに別の従来例としては、入力画像の画素レベルの頻度分布に応じて、レベル変換関数を適応的に変化させる方法があり、その代表例としてはヒストグラムイコライゼーションと呼ばれる方法が挙げられる。図20(A),(B)に、このヒストグラムイコライゼーションの原理を示す。図20(A)において、横軸は入力画像の画素レベル(入力レベル)lを、縦軸は度数(頻度または累積頻度)を表す。Fmaxは、累積頻度の最大値であり、頻度を算出するために用いる画素の総数である。この方法では、図20(A)に示したように、はじめに入力画像の画素レベルlに関する頻度分布H(l)が生成され、次に以下の(3)式を用いて累積頻度分布C(l)が生成される。
【0007】
【数3】
【0008】
この累積頻度分布C(l)の縦軸を、以下の(4)式を用いて出力画像が取り得るレベル範囲に正規化することにより、レベル変換関数T(l)が生成される(図20(B))。この関数T(l)を用いることにより、出現頻度の高いレベルによって構成される領域(面積が大きい領域)のコントラストを増加させることが可能となる。
【0009】
【数4】
【0010】
入力された画像を、よりダイナミックレンジの小さい、すなわち画素レベルを表現するビット数が少ない環境で利用する場合(ビット数の少ない伝送路で伝送する場合や、表示装置に表示する場合、あるいは記憶装置に保存する場合など)には、ダイナミックレンジの圧縮を行う必要がある。従来は、このような目的でのダイナミックレンジの圧縮処理にも、上述した方法と同様のレベル変換が用いられている。ただし、この場合には、レベル変換関数の出力画像の最大レベルが、入力画像のそれよりも小さい値となる。
【0011】
一方、文献「Z. Rahman, et, alt.:"A Multiscale retinex for color rendition and dynamic range compression in Applications of Digital image Processing", XIX Proc. SPIE 2847 (1996)」では、空間的に緩やかに変化する照明光の成分をローパスフィルタによって抽出し、これを圧縮することで全体的なダイナミックレンジを圧縮する方法が提案されている(以下、この方法を「Multiscale retinex法」と記す。)。照明成分の抽出には帯域の狭い線形ローパスフィルタが用いられている。この方法では、以下の(5)式に示すように、入力画素の値I(x,y)、およびローパスフィルタ出力LPF(I(x,y))の対数値を取り、前者から後者を差し引くことによりダイナミックレンジの圧縮が行われる。
【0012】
【数5】
【0013】
【発明が解決しようとする課題】
ところで、上述した従来のレベル変換方法では、不自然な画像が生成されるのを回避するため、単調増加性を有するレベル変換関数を用いている。このため、いずれかのレベル範囲のコントラスト(レベル変換関数の傾き)を増加させた場合、他のレベル範囲では逆にコントラストが低下してしまうといった問題がある。
【0014】
また、Multiscale retinex法では、単調増加性を犠牲にすることで、よりコントラストの高い画像を再現することを可能としているが、照明条件が急激に切り替わる場合、線形フィルタではその変化を抽出することができず、主観的に好ましくないノイズが発生してしまうという問題がある。
【0015】
例えば、図21に示すように、照明条件の異なる2つの領域が隣接している画像(図中、実線)に対して線形ローパスフィルタを施すと、図中、細い破線で示したように境界のぼけた信号がフィルタ出力として得られる。これを照明成分と見なした場合、照明境界の左側の領域(B領域)において、境界付近(BNB領域)は、境界より離れた部分(BFB領域)よりも照明レベルが低いことになる。上述の(5)式は入力信号を照明成分で割り算することと等価であり、照明成分が大きいほど大きく圧縮されることを意味するため、結果として再現画像(図中、太い破線)のBNB領域にはオーバーシュートが発生する。逆に、照明境界の右側の領域(D領域)において、境界付近(DNB領域)は、境界より離れた部分(DFB)に比べて照明レベルが高いと見なされてアンダーシュートが発生する。Multiscale retinex法ではこの問題を回避するためにスケールの異なる複数の線形ローパスフィルタを用い、それぞれのローパスフィルタによって得られる結果を線形荷重によって合成する方法を用いているが、各スケールに対する重みは固定されており、上記の問題を十分に抑制できていない。
【0016】
そこで、照明成分の抽出に線形ローパスフィルタではなく、例えばイプシロンフィルタなどの非線形フィルタを用いることが考えられる。イプシロンフィルタは、線形フィルタに比べエッジを保存する能力にすぐれ、異なる照明光が存在する画像においても照明成分をより有効に抽出することができる。しかしながら、ノイズ除去などを目的に一般的に用いられているしきい値固定のイプシロンフィルタでは、その出力にエッジ近傍において不連続な波形が形成されるため、これをダイナミックレンジの圧縮用に用いた場合、圧縮後の再現画像に原画像にはない不自然な画像パターンが発生する可能性がある。
【0017】
本発明はかかる問題点に鑑みてなされたもので、その目的は、イプシロンフィルタを用いて、複数の異なる照明が存在する場合にも、それらの境界を適切に抽出することを可能とし、不自然な画像パターンの発生を抑えて、主観的に好ましいダイナミックレンジの圧縮を実現できる画像処理方法および装置を提供することにある。
【0018】
【課題を解決するための手段】
本発明による画像処理方法は、入力画像上の各位置ごとにエッジ強度を算出するエッジ強度算出過程と、算出されたエッジ強度に基づいて、イプシロンフィルタのしきい値を制御するしきい値制御過程と、しきい値制御過程において制御されたしきい値を用いて、入力画像にイプシロンフィルタによるフィルタ処理を施すフィルタリング過程と、フィルタリング過程の出力値に応じて、画素値を変換するための係数を算出し、その算出した係数により、各画素ごとに画素値の変換を行う画素値変換過程とを含むものである。
【0019】
本発明による画像処理装置は、入力画像上の各位置ごとにエッジ強度を算出するエッジ強度算出手段と、設定されたしきい値を用いて、入力画像にフィルタ処理を施すイプシロンフィルタと、エッジ強度算出手段によって算出されたエッジ強度に基づいて、イプシロンフィルタで用いられるしきい値を制御するしきい値制御手段と、イプシロンフィルタからの出力値に応じて、画素値を変換するための係数を算出し、その算出した係数により、各画素ごとに画素値の変換を行う画素値変換手段とを備えたものである。
【0020】
本発明による画像処理方法および装置では、入力画像上の各位置ごとにエッジ強度が算出され、その算出されたエッジ強度に基づいて、イプシロンフィルタで用いられるしきい値が制御される。そして、イプシロンフィルタからの出力値に応じて、画素値を変換するための係数が算出され、その算出された係数により、各画素ごとに画素値の変換が行われる。これにより、複数の異なる照明が存在する場合にも、イプシロンフィルタによるフィルタ処理において、それらの境界の抽出を適切に行うことが可能とされる。
【0022】
また、本発明による画像処理方法および装置では、大きさの異なる2つのしきい値を算出して、しきい値の制御を行う。このとき、イプシロンフィルタでは、近傍画素の値が注目画素に比べて大きい場合と小さい場合とで異なるしきい値を用いてフィルタ処理を行う。これによって、フィルタ処理において、画素値の変化における照明レベルの影響が軽減され、より適切に照明成分の抽出がなされる。
【0023】
【発明の実施の形態】
以下、本発明の実施の形態について図面を参照して詳細に説明する。
【0024】
[第1の実施の形態]
まず、本実施の形態に係る画像処理装置において処理される入力画像信号につて説明する。本画像処理装置において処理される入力画像信号は、2次元ディジタル画像を図2に示すように水平方向、垂直方向の順に走査して得られた時系列な画素値の信号である。本実施の形態では、2次元画像上の任意の位置(x,y)に対応する画素値をI(x,y)と表し、これを入力画像信号として処理する。
【0025】
次に、本実施の形態に係る画像処理装置の構成について説明する。本画像処理装置は、図1に示したように、エッジ強度算出器10と、しきい値制御器11と、イプシロンフィルタ12と、除算器13と、レベル変換器14と、乗算器15とを備えている。
【0026】
エッジ強度算出器10は、入力画像の各位置において画素値I(x,y)のエッジ強度G(x,y)を算出する機能を有している。エッジ強度G(x,y)としては、例えば以下の(6)式で与えられるようなI(x,y)の1次微分値を用いることができる。
【0027】
【数6】
【0028】
あるいは、ノイズの影響を抑えるために平滑化効果を持つ以下の(7)式による値をエッジ強度G(x,y)として用いることも可能である。
【0029】
【数7】
【0030】
ここで、(6),(7)式において、dは微分を算出するための微小距離を示す定数である。エッジ強度算出器10において算出されたエッジ強度G(x,y)は、しきい値制御器11に送られる。
【0031】
しきい値制御器11は、エッジ強度算出器10によって算出されたエッジ強度G(x,y)に基づき、後述のイプシロンフィルタ12で用いるしきい値E(x,y)の大きさを画素ごとに決定する機能を有している。しきい値制御器11の機能により、しきい値E(x,y)は、例えば以下の(8)式を用いることにより、エッジ強度G(x,y)が大きいほど小さな値となるように制御される。
【0032】
【数8】
【0033】
(8)式において、Gmin,Gmax,Emin,Emaxは、エッジ強度G(x,y)をしきい値E(x,y)に変換するための定数であり、それぞれエッジ強度の最小値、最大値、しきい値E(x,y)の最小値、最大値を表している。しきい値制御器11において決定されたしきい値E(x,y)は、イプシロンフィルタ12に送られる。
【0034】
イプシロンフィルタ12は、図4に示したように、例えば、差分器20と、絶対値算出器21と、比較器22と、線形ローパスフィルタ(図ではLPFと記す。)23とを有して構成されている。このイプシロンフィルタ12は、2次元のフィルタであり、しきい値制御器11によって決定されたしきい値E(x,y)を用いて、入力画像に非線形なフィルタ処理を施す機能を有している。イプシロンフィルタ12の出力R(x,y)は、照明成分として除算器13およびレベル変換器14に送られる。
【0035】
除算器13は、入力画像からイプシロンフィルタ12で算出された照明成分を除去するために、以下の(9)式に示すように、入力画像の各画素値I(x,y)を照明成分R(x,y)で割り算するようになっている。割り算した結果得られる非照明成分S(x,y)は、乗算器15に送られる。
【0036】
【数9】
【0037】
レベル変換器14は、イプシロンフィルタ12で算出された照明成分R(x,y)を、以下の(10)式に示すように、レベル変換関数T(l)によってレベル変換することによって圧縮し、補正照明成分CR(x,y)を算出する機能を有している。
【0038】
【数10】
【0039】
レベル変換器14で用いるレベル変換関数T(l)としては、例えば図3に示すような関数を用いることができる。なお、図3において、Rmax,CRmaxは、それぞれ入力レベルおよび出力レベルの最大値を表している。
【0040】
乗算器15は、以下の(11)式に示すように、非照明成分S(x,y)に補正照明成分CR(x,y)を積算することで画像信号を復元するようになっている。復元された結果を示す画像信号O(x,y)は、図示しない伝送路、記憶装置、または表示装置などに出力される。
【0041】
【数11】
【0042】
なお、本実施の形態において、除算器13、レベル変換器14および乗算器15が、本発明における「画素値変換手段」の一具体例に対応する。
【0043】
次に、以上のように構成された画像処理装置の作用、動作を説明する。なお、以下の説明は、本実施の形態に係る画像処理方法の説明を兼ねている。
【0044】
本画像処理装置において、入力画像を示す信号は、エッジ強度算出器10、イプシロンフィルタ12および乗算器13に入力される。まず、エッジ強度算出器10では、入力画像の各位置ごとにエッジの大きさ、すなわちエッジ強度G(x,y)を算出する。このとき、エッジ強度算出器10は、例えば上述の(6)式または(7)式を用いることにより、注目画素の近傍領域内の画素値の1次微分値が大きいほど大きな値となるようにして、エッジ強度G(x,y)を算出する。エッジ強度算出器10は、算出したエッジ強度G(x,y)をしきい値制御器11に出力する。
【0045】
しきい値制御器11では、エッジ強度G(x,y)に基づいて、イプシロンフィルタ12のしきい値Eを制御する。より詳しくは、しきい値制御器11は、例えば上述の(8)式を用いて、しきい値E(x,y)の大きさを画素ごとに決定し、エッジ強度G(x,y)が大きいほどしきい値E(x,y)が小さくなるように制御する。しきい値制御器11は、決定されたしきい値E(x,y)をイプシロンフィルタ12に出力する。
【0046】
イプシロンフィルタ12では、しきい値制御器11によって決定されたしきい値E(x,y)を用いて、入力画像にフィルタ処理を施す。
【0047】
イプシロンフィルタ12でのフィルタ処理は、より詳しくは、例えば図4に示した構成により、以下のように行われる。イプシロンフィルタ12において、差分器20には、図4に示したように、現在の注目画素の値I(x,y)を示す信号とその近傍領域NB内の画素の値I(x+dx,y+dy)を示す信号とが入力される。差分器20では、注目画素の値I(x,y)とその近傍領域NB内の画素の値I(x+dx,y+dy)との差分を算出する。差分器20では、近傍領域NB内のすべての画素に対して順次この差分値を計算し、各近傍画素に対応づけて、その値D(dx,dy)を絶対値算出器21に出力する。
【0048】
絶対値算出器21では、差分器20より送られてくる各差分値D(dx,dy)の絶対値AD(dx,dy)を算出する。絶対値算出器21では、算出した絶対値AD(dx,dy)を比較器22に出力する。
【0049】
比較器22には、絶対値算出器21で算出された絶対値AD(dx,dy)が入力されると共に、注目画素の値I(x,y)を示す信号およびその近傍領域NB内の画素の値I(x+dx,y+dy)を示す信号、ならびに、しきい値制御器11で決定されたしきい値E(x,y)が入力される。比較器22では、以下の(12)式に示すように、絶対値AD(dx,dy)としきい値E(x,y)とを比較し、その結果に応じて、注目画素の値I(x,y)または近傍画素の値I(x+dx,y+dy)のいずれか一方を選択し、それを値J(dx,dy)として線形ローパスフィルタ23に出力する。
【0050】
【数12】
【0051】
線形ローパスフィルタ23では、近傍領域NB内のすべての画素に対応する値J(dx,dy)が比較器22によって算出された時点で、以下の(13)式による加重平均値R(x,y)を算出する。
【0052】
【数13】
【0053】
ここで、NBはフィルタリング処理における近傍領域を定義する相対的な座標の集合である。また、a(dx,dy)は各画素値に対する重み係数であり、この線形ローパスフィルタ23としては、例えば以下の(14)式に示すような平均値フィルタなどを用いることができる。
【0054】
【数14】
【0055】
(14)式において、Nは近傍領域内NBの画素の数を表している。なお、イプシロンフィルタ12の目的は、画像中の細かい構造を除去し、まとまった領域を抽出することであり、その近傍領域は大きい方が望ましい。
【0056】
以上のようにしてイプシロンフィルタ12で得られた値R(x,y)は、近似的に画像に含まれる照明成分を表しているものと考えられる。イプシロンフィルタ12は、値R(x,y)を照明成分として除算器13およびレベル変換器14に出力する。
【0057】
除算器13では、上述の(9)式に示したように、入力画像の各画素値I(x,y)を照明成分R(x,y)で除算することにより、入力画像からイプシロンフィルタ12で算出された照明成分を除去し、その結果得られる非照明成分S(x,y)を乗算器15に出力する。
【0058】
一方、レベル変換器14では、イプシロンフィルタ12で算出された照明成分R(x,y)を、例えば図3に示すようなレベル変換関数T(l)によってレベル変換することによって圧縮し、補正照明成分CR(x,y)を算出する。レベル変換器14は、算出した補正照明成分CR(x,y)を乗算器15に出力する。
【0059】
乗算器15では、除算器13からの出力である非照明成分S(x,y)にレベル変換器14からの出力である補正照明成分CR(x,y)を積算することで画像信号を復元する。ここで、以上の除算器13、レベル変換器14および乗算器15での全体の演算について考察すると、非照明成分S(x,y)に補正照明成分CR(x,y)を積算することは、後述する(16)式に示すように、イプシロンフィルタ12からの出力値R(x,y)に応じて画素値を変換するための係数F(R(x,y))を算出し、それを対応する入力画素値I(x,y)に積算することで、各画素ごとに画素値の変換を行ってダイナミックレンジの圧縮を行っていることに相当する。
【0060】
以上のようにして乗算器15から出力される画像信号O(x,y)は、入力画像よりも相対的にダイナミックレンジの狭い画像装置、すなわち画素レベルを表現するビット数が少ない環境(ビット数の少ない伝送路で伝送する場合や、表示装置に表示する場合、あるいは記憶装置に保存する場合など)において利用される。
【0061】
次に、図5および図6を参照して、本実施の形態におけるイプシロンフィルタ12を用いてダイナミックレンジの圧縮を行った場合の、従来法(単調増加性を有するレベル変換関数を用いて圧縮を行う場合)に対する有効性について説明する。
【0062】
図5(A)は、入力画像の画素値I(x,y)、および本実施の形態におけるイプシロンフィルタ12からの出力R(x,y)を1次元信号として表した図である。図5(B)は、図5(A)に示した入力画像に対して従来のレベル変換法によるダイナミックレンジの圧縮を行った結果(再現画像)を示し、図5(C)は本実施の形態におけるダイナミックレンジの圧縮を行った結果を示す。
【0063】
また、図6(A),(B)には、各領域の画素レベルとレベル変換関数T(l)、および係数算出関数F(l)との関係を示す。ここで、係数算出関数F(l)とは、レベル変換関数T(l)を用いて以下の(15)式のように定義される。
【0064】
【数15】
【0065】
この係数算出関数F(l)を用い、上述の(9)式および(10)式を考慮することで、出力画像O(x,y)を与える(11)式は、以下の(16)式のように書き換えることが可能となる。
【0066】
【数16】
【0067】
この(16)式は、照明成分R(x,y)の圧縮によるダイナミックレンジの圧縮が、画素ごとに算出される係数F(R(x,y))を、対応する入力画素値I(x,y)に積算することで実現できることを示している。このとき、係数算出関数F(l)は、イプシロンフィルタ12の出力値を各画素に施すゲイン係数に変換する機能を持つ。なお、図6(B)における係数算出関数F(l)の値の最小値Cminは、次の(17)式で与えられることになる。
【0068】
【数17】
【0069】
図5(B)からも分かるように、従来法は低レベル領域(図5(A)に示すレベルl1とレベルl3とで構成されている領域)でのコントラストを保存することができるが、高レベル領域(レベルl4とレベルl6とで構成されている領域)においてはコントラストの低下を招いている。これは、変曲点lk以上の高いレベルに対するレベル変換関数T(l)の傾きの影響を直接的に受けた結果である。従来法においてコントラストを向上させるためにはレベル変換関数T(l)の傾きを大きくする必要がある。
【0070】
これに対し、本実施の形態(図5(C))では、高レベル領域および低レベル領域のそれぞれには、係数算出関数F(l)によって与えられる単一の補正係数が施されるため、各領域内のコントラストはこの補正係数の大きさに依存することになる。本実施の形態では、低レベル領域に対してはその平均レベルl2によって決まる補正係数が一様に施されるが、その値は、レベルl1およびレベルl3に対するものと同じ1.0であり、従来法と同程度のコントラストが得られる。また、高レベル領域に対しては、その平均値l5で決まる一定の補正係数c5が施されるため、レベルl4の部分とレベルl6の部分の間のコントラストは、このゲインで確保されることになる。
【0071】
実際、図6(A)に示すような、傾きの異なる2本の直線からなるレベル変換関数T(l)を用い、その変曲点レベルlk以上に対応する直線が以下の(18)式で表される場合には、従来法における高レベル領域のコントラストはその直線の傾きaに依存し、「a(l6−l4)/l5」で与えられる。ただし、ここではコントラストを「(最大レベル最小レベル)/平均レベル」で定義している。
【0072】
【数18】
【0073】
一方、本実施の形態において高レベル領域に適用される補正係数c5は、以下の(19)式で与えられる。
【0074】
【数19】
【0075】
従って、本実施の形態では、この領域のコントラストは、c5(l6−l4)/l5={a(l6−l4)/l5}+{b/(l5*l5)}となるが、ダイナミックレンジの圧縮におけるレベル変換関数T(l)の高レベルにおける傾きは通常1.0より小さいため、切片bは常に正の値を取る。このことは、従来法に比べて、イプシロンフィルタ12を用いた本実施の形態による方法がより高いコントラストを実現できることを示している。
【0076】
このように、本実施の形態においては、イプシロンフィルタ12によって抽出された領域内のコントラストが、係数算出関数F(l)によって与えられる補正係数の値そのものによって決まり、レベル変換関数T(l)の傾きは領域間のコントラストに影響を及ぼすことになる。従って、本実施の形態によれば、領域間のコントラストを圧縮することにより、領域内のコントラストを保存することができ、主観的に好ましい出力画像を得ることができる。
【0077】
次に、線形ローパスフィルタを用いた従来法に対する本実施の形態の効果について説明する。ただしここでは、説明を簡単にするため画像を1次元信号として表す。従来法の問題点は既に[発明が解決しようとする課題]の項目において図21を用いて説明したとおりであるが、その問題点を解決するためには、照明条件の異なる領域の境界(照明境界)を保存したまま、同一照明下の領域内を平滑化する必要がある。ところで、経験上、照明強度の変化に起因する画素レベルの変化は、物体表面の反射率に起因する画素レベルの変化よりもはるかに大きく、その結果、照明境界では画素レベルの大きなエッジが発生することになる。
【0078】
図7(A)に示すように、このようなエッジ周辺では、イプシロンフィルタ12の注目画素の近傍領域NB内にしきい値E(x,y)を超える大きな差分絶対値AD(dx,dy)を与える画素値I(x+dx,y+dy)が存在することになる。この画素値I(x+dx,y+dy)は、比較器22(図4)によって、図7(A)の太い破線(I'(dx,dy))で示すように、現在の注目画素の値(近傍領域NB中央の値)I(x,y)に置き換えられるため、(13)式による平滑化には大きく寄与せず、結果的にエッジの形状が保存されることになる。これに対し、照明境界以外の部分では、図7(B)に示すように、画素レベルの変化はそれほど大きなものにはならず、注目画素の近傍領域NB全体に渡り、差分絶対値AD(dx,dy)がしきい値E(x,y)よりも小さくなる。この場合、(12)式のJ(dx,dy)はすべて入力画素値I(x+dx,y+dy)と等しく、イプシロンフィルタは単純な線形ローパスフィルタと等価になり、近傍領域NB全体に渡る平滑化が行われる。
【0079】
このようにイプシロンフィルタ12は、線形フィルタに比べエッジを保存する能力にすぐれ、異なる照明光が存在する画像においても照明成分をより有効に抽出することができる。しかしながら、ノイズ除去などを目的に一般的に用いられているしきい値固定の従来のイプシロンフィルタでは、その出力にエッジ近傍において不連続な波形が形成されるため、これをダイナミックレンジの圧縮用に用いた場合、圧縮後の再現画像に原画像にはない不自然な画像パターンが発生する可能性がある。
【0080】
この問題点を説明するために、図8(A),(B)にモデル化したエッジとその周辺におけるしきい値固定の従来のイプシロンフィルタの出力を示す。ここではエッジ中央部から高レベル側のみを考え、レベルが急激に変化するエッジ部分は傾きaの直線で、それ以外のエッジ周辺部分は傾き0の平坦な直線で近似している。図8(A),(B)では、入力信号を細い実線81で、イプシロンフィルタの出力を太い実線82で、対応する線形ローパスフィルタの出力を太い破線83で示す。また、用いるイプシロンフィルタでの注目画素の近傍領域の大きさをN、そのしきい値をEとする。図8(A)は、a,N,Eの間に以下の(20)式の関係が成立している場合のイプシロンフィルタの出力である。
【0081】
【数20】
【0082】
図8(A)において、横軸は空間的な位置座標を表し、右に行くほど座標値が大きくなるものとする。p0は、エッジの立下り点peからN/2だけ離れた位置を示している。この場合、イプシロンフィルタは線形ローパスフィルタと同等に振る舞う。しかしながら、照明光の変化による画像上のレベル変化は急峻に生じること(傾きaが大きい)、また上述したように効果的なダイナミックレンジの圧縮を行うためには大きなフィルタを用いる必要があること(Nが大きい)から、照明条件が変化する照明境界では通常(20)式の条件は成立しないと考えられる。
【0083】
一方、図8(B)には、(20)式が成立しない場合のイプシロンフィルタの出力を示す。図8(B)において、出力波形を右側から左側へと見ていくと、p0からp0−E/aまでは線形フィルタと同じ出力となり、その後peまでは一定の値を出力する。peからpe−E/aまでは2次曲線に従って下降するが、(20)式が成立しない条件のもとでは、この間に必ず入力信号を示す直線81と交差する。pe−E/aでは不連続的に入力信号と同じ値になり、その後入力信号をそのまま出力することになる。しきい値固定のイプシロンフィルタのこの振る舞いは、単純化したエッジモデルに対するものであるが、エッジ付近で複雑な波形を出力することは明らかである。特に、pe−E/aのあたりでは不連続に変化し、入力信号よりも大きな傾きを持った波形を出力する。
【0084】
出力画像O(x,y)が自然な再現画像となるためには、画像上の各位置における局所的な形状(空間的勾配の向き)が保存されていることが必要となる。すなわち、以下の(21)式に示すように、出力画像O(x,y)の微分値O'(x,y)の符号と入力画像I(x,y)の微分値I'(x,y)の符号とが一致しなければならない。
【0085】
【数21】
【0086】
ここで、sign(x)はxの符号を表す。(21)式の条件が満たされない場合、入力画像と出力画像との間でレベル勾配の逆転が生じ、入力画像には存在しない画像パターンが出力画像に現れることになる。例えば図9(A)に示す入力画像が、図9(B)に示すように出力される。図9(B)の出力画像では、レベル勾配の逆転により、入力画像には存在しないパターン90が出現している。しかし、図8(B)に示すエッジ付近におけるイプシロンフィルタの不連続な振る舞いは、この(21)式の条件の成立を困難にする可能性があり、レベル勾配の逆転が生じるおそれがある。
【0087】
(21)式が成立するか否かは、イプシロンフィルタの振る舞いとその出力に施されるレベル変換関数T(l)に依存する。このことを明らかにするために、はじめに(11)式に(9)式を代入して入力画像I(x,y)と出力画像O(x,y)との関係(22)式を得る。
【0088】
【数22】
【0089】
(22)式の右辺におけるCR(x,y)/R(x,y)は、先に述べた係数算出関数F(R(x,y))に相当する。(22)式の両辺を微分し、(21)式に代入することで、出力画像O(x,y)が満たすべき条件は以下の(23)式のように表されることになる。
【0090】
【数23】
【0091】
この条件に対するレベル変換関数T(l)の関与をさらに明確にするために、図10にイプシロンフィルタの出力に対するレベル変換関数T(l)の出力の関係を示す。横軸はイプシロンフィルタの出力(レベル変換関数T(l)への入力)、縦軸はレベル変換関数T(l)の出力を表している。いま、ある画素のイプシロンフィルタの出力をR、そのレベル変換後の値をCRとし、これらの比の値CR/Rをaとする。すなわち、aは原点と(R,CR)を通過する直線の傾きに相当する((24)式)。
【0092】
【数24】
【0093】
またこの同じ画素におけるイプシロンフィルタの出力の空間微分値をR'とすると、これは図10の横軸上におけるRを中心とする微小範囲の大きさに相当するものと考えられる(ただしRは微分値を算出する画像上の方向に応じて適当な符号を有する)。従って、Rに対応するレベル変換関数T(l)のlに関する微係数をbとすると、レベル変換後の空間微分値は以下の(25)式によって近似できる。
【0094】
【数25】
【0095】
(24)式および(25)式を(22)式に代入することで、イプシロンフィルタの出力に関する部分R'/Rと、レベル変換関数T(l)の性質に関する部分(1−b/a)とからなる以下の(26)式の条件を得ることができる。
【0096】
【数26】
【0097】
なお、照明成分圧縮のためのレベル変換関数T(l)は、通常単調増加性を有するため、aおよびbは正の値となり、以下の(27)式が成立する。
【0098】
【数27】
【0099】
なお、1−b/a≧0は、a≧bを意味するが、この場合には照明成分レベルRの近辺ではレベルが大きいほど大きな圧縮を施されることを意味している。なぜならば、レベルRの圧縮率はaであり、その近辺を一様に圧縮する場合にはレベルR近辺におけるレベル変換関数T(l)の傾きはaである必要があるが、実際の傾きはより小さいbであり、図11(A)に示すように高レベル側ほど大きく圧縮されることになる。逆に1−b/a、すなわちa<bの場合には、図11(B)に示すように低レベル側ほど大きく圧縮されることになる。なお、図11(A),(B)の縦軸および横軸は、図10と同様、それぞれイプシロンフィルタの出力(レベル変換関数T(l)への入力)およびレベル変換関数T(l)の出力を表す。
【0100】
エッジ付近以外の平坦部分では、イプシロンフィルタはローパスフィルタとして機能するため、一般的にその出力における空間的な変化の割合は、入力画像のそれよりも緩やかなものとなる。従って、これらの部分では以下の(28)式((28A),(28B))に示す2つの条件が成立していると仮定できる。
【0101】
【数28】
【0102】
I'≧0の場合、(28)式よりR'≧0となるため、1−b/aが負であれば常に(26)式が成立することになる。また、0≦1−b/a≦1である場合には、(28)式から得られるI'/I≧R'/Rより、以下の(29)式のようになる。
【0103】
【数29】
【0104】
同様に、I'<0の場合にも(26)式が成立することは容易に分かる。
【0105】
一方、図8(B)に示すとおり、エッジ中央部分においてはイプシロンフィルタは入力信号をそのまま出力することになるため、この部分では以下の(30)式が成立すると考えられるが、この場合にも(26)式が満足されることは明らかである。
【0106】
【数30】
【0107】
本実施の形態においてイプシロンフィルタ12を除去し、R(x,y)=I(x,y)となるようにすれば、その結果はレベル変換のみによる従来法と等価となるが、これは画像全体に渡って(30)式の条件が成立している場合に相当し、従来法ではレベル勾配の逆転が発生しないことになる。
【0108】
これらの条件が満たされない部分では、(26)式が成立するか否かはイプシロンフィルタ12の出力とレベル変換関数T(l)の特性に依存する。レベル変換関数T(l)の特性は照明成分をどのように圧縮するかによって決められるべきものであるため、ここではイプシロンフィルタ12の出力がどのように(26)式の成立に影響を及ぼすかを4つの場合(1)〜(4)に分けて説明する。
【0109】
(1)I'≧0かつR'が1−b/aと同じ符号の場合;
右辺、左辺がともに正となるため、R'/Rの絶対値が大きいほど(26)式の成立が困難となる。
【0110】
(2)I'≧0かつR'が1−b/aと異なる符号の場合;
右辺が正、左辺が負となるため、常に(26)式が成立する。
【0111】
(3)I'<0かつR'が1−b/aと同じ符号の場合;
右辺が負、左辺が正となるため、常に(26)式が成立する。
【0112】
(4)I'<0かつR'が1−b/aと異なる符号の場合;
右辺、左辺がともに負となるため、R'/Rが大きいほど(26)式の成立が困難となる。
【0113】
上記(2)および(3)の条件のもとでは、近傍領域内において入力画像のレベルが高いほど少ない圧縮が施されることになるため、勾配の逆転が発生する可能性がない。従ってR'/Rがどのような値になっても常に(26)式が成立することになる。これに対して(1),(4)の条件では、R'/Rの絶対値が大きいほど(26)が満たされない可能性が高くなるが、少なくともその比が入力画像と同じであれば(26)式が成立することは上記に述べたとおりである。本実施の形態では、イプシロンフィルタ12のしきい値Eを可変とし、エッジ部分においては入力画像にできる限り近い信号を出力することで、すなわち、エッジ強度G(x,y)が大きいほど、イプシロンフィルタ12のしきい値Eが小さくなるようにしきい値を制御することで、勾配の逆転を最小限に抑制し、図9(B)に示したような不自然なパターン90が生じることのないよう、自然な画像を再現することを可能としている。
【0114】
なお、図3に示したレベル変換関数T(l)はあくまでも一例であり、目的に応じて任意のものを用いることができることはいうまでもない。例えば、(1)式や(2)式に示した関数を用いても良い。
【0115】
以上説明したように、本実施の形態によれば、入力画像上の各位置ごとにエッジ強度G(x,y)を算出すると共に、このエッジ強度G(x,y)に基づいて、イプシロンフィルタ12のしきい値E(x,y)を制御して入力画像にフィルタ処理を施し、また、そのイプシロンフィルタ12からの出力値R(x,y)に応じて算出される係数F(R(x,y))を、(16)式に示したように入力画素値I(x,y)に積算することによって、各画素ごとに画素値の変換を行ってダイナミックレンジの圧縮を行うようにしたので、複数の異なる照明が存在する場合にも、それらの境界を適切に抽出することを可能とし、不自然な画像パターンの発生を抑えて、主観的に好ましいダイナミックレンジの圧縮を実現できる。すなわち、イプシロンフィルタ12を用いて、入力画像から照明成分を抽出し、その照明成分を圧縮することにより、局所的なコントラストを保存したまま全体的なダイナミックレンジを削減し、主観的に好ましい再現画像を得ることができる。このとき、エッジ強度G(x,y)に基づいて、イプシロンフィルタ12のしきい値Eを画素値I(x,y)の局所的な勾配に応じて適応的に変化させるようにしたので、線形ローパスフィルタを用いた場合やしきい値固定のイプシロンフィルタを用いた場合よりも、より正確に照明境界を抽出することができる。
【0116】
[第2の実施の形態]
次に、本発明の第2の実施の形態について説明する。なお、以下の説明では、上記第1の実施の形態における構成要素と実質的に同一の機能を有する部分には同一の符号を付し、適宜説明を省略する。
【0117】
図12は、本発明の第2の実施の形態に係る画像処理装置の構成を示している。本実施の形態に係る画像処理装置の全体的な機能は、上記第1の実施の形態と同様であるが、本実施の形態では、除算器13およびレベル変換器14(図1)の代わりに、これらの機能をあわせ持った係数算出器16が設置されている点が、上記第1の実施の形態とは異なる。すなわち、本実施の形態は、上記第1の実施の形態で示した(16)式に基づいてダイナミックレンジの圧縮を行うものである。
【0118】
本実施の形態における係数算出器16では、イプシロンフィルタ12の出力R(x,y)に対して、図13に示すような係数算出関数F(l)を施して係数C(x,y)を算出する。係数算出関数F(l)は、上記第1の実施の形態で述べたように、レベル変換関数T(l)を用いて(15)式によって得られる。係数算出器16によって算出された係数C(x,y)は、乗算器15に送られる。
【0119】
本実施の形態における乗算器15には、入力画像を示す信号が直接入力されると共に、算出器16によって算出された係数C(x,y)が入力される。乗算器15では、入力画像の各画素値I(x,y)に対して対応する係数C(x,y)を積算することで画像信号を復元し、上記第1の実施の形態と同様、その結果の画像信号O(x,y)を図示しない伝送路、記憶装置または表示装置などに出力する。
【0120】
本実施の形態においても、エッジ強度G(x,y)に基づいて、イプシロンフィルタ12のしきい値E(x,y)を制御して入力画像にフィルタ処理を施し、また、そのイプシロンフィルタ12からの出力値R(x,y)に応じて算出される係数F(R(x,y))を、入力画素値I(x,y)に積算することによって、各画素ごとに画素値の変換を行ってダイナミックレンジの圧縮を行うようにしているので、上記第1の実施の形態と同様の効果を得ることができる。
【0121】
[第3の実施の形態]
次に、本発明の第3の実施の形態について説明する。なお、以下の説明では、上記第1の実施の形態または上記第2の実施の形態における構成要素と実質的に同一の機能を有する部分には同一の符号を付し、適宜説明を省略する。
【0122】
図14は、本発明の第3の実施の形態に係る画像処理装置の構成を示している。本実施の形態に係る画像処理装置の構成は、上記第1の実施の形態(図1)とほぼ同じであるが、しきい値制御器11に、エッジ強度算出器10からの出力であるエッジ強度G(x,y)のほかに、さらに、入力画像の画素値I(x,y)が直接入力されている点が異なる。
【0123】
本実施の形態におけるしきい値制御器11Aでは、後段のイプシロンフィルタ12で用いるしきい値E(x,y)を、エッジ強度G(x,y)のみではなく、入力画像の画素レベルによっても制御する。より具体的には、しきい値制御器11Aでは、入力画像の画素値I(x,y)が大きいほどしきい値E(x,y)が大きくなると共に、エッジ強度G(x,y)が大きいほどしきい値E(x,y)が小さくなるようにしきい値E(x,y)を制御する。
【0124】
このようなしきい値制御は、以下のようにして実現できる。例えば、まず、あらかじめ設定された正の係数r(r≦1.0)を用いて、以下の(31)式によって、入力画像の画素値I(x,y)が大きいほどその値が大きくなるような、仮のしきい値Etmp(x,y)を設定する。
【0125】
【数31】
【0126】
その後、このしきい値Etmp(x,y)をエッジ強度G(x,y)により修正して実際に用いるしきい値E(x,y)を決定する。例えば、(8)式における正規化のための定数Emin,Emaxをそれぞれ0.0,1.0として、(8)式によってエッジ強度に応じた係数G(x,y)を算出するようにする。これを以下の(32)式のように、仮のしきい値Etmp(x,y)に積算して最終的なしきい値E(x,y)を求める。これにより、入力画像の画素値I(x,y)が大きいほどその値が大きくなると共に、エッジ強度G(x,y)が大きいほどその値が小さくなるよう、しきい値E(x,y)が制御される。
【0127】
【数32】
【0128】
イプシロンフィルタ12のしきい値E(x,y)は、入力画素の空間的変化が照明成分の変化によるものか、物体表面の反射率の変化によるものかを分離する役割を担うが、物体表面の反射率の変化が小さい場合でも照明レベル自体が大きいと画素値I(x,y)の変化は大きくなる。すなわち、照明レベルの大きな変化と強い照明下における小さな反射率の変化とをどのように区別するかが問題となる。本実施の形態では、画素レベルが大きいほどイプシロンフィルタ12のしきい値E(x,y)を大きく設定するようにしているが、これにより画素値I(x,y)の変化における照明レベルの影響を軽減でき、より適切に照明成分を抽出することが可能となる。
【0129】
以上説明したように、本実施の形態によれば、エッジ強度G(x,y)のみならず、入力画像の画素レベルを考慮して、イプシロンフィルタ12のしきい値E(x,y)の制御を行うようにしたので、イプシロンフィルタ12において、より適切に照明成分を抽出することが可能となる。
【0130】
[第4の実施の形態]
次に、本発明の第4の実施の形態について説明する。なお、以下の説明では、上記第1〜第3の実施の形態における構成要素と実質的に同一の機能を有する部分には同一の符号を付し、適宜説明を省略する。
【0131】
図15は、本発明の第4の実施の形態に係る画像処理装置の構成を示している。本実施の形態に係る画像処理装置の構成は、上記第3の実施の形態(図14)とほぼ同じであるが、特に、各画素ごとに2種類のしきい値Elo(x,y),Eup(x,y)を用いてイプシロンフィルタのしきい値制御を行うようにした点が上記第3の実施の形態とは異なる。
【0132】
本実施の形態におけるしきい値制御器11Bでは、各画素ごとに大きさの異なる2種類のしきい値Elo(x,y),Eup(x,y)を算出し、それらによってイプシロンフィルタ12Aのしきい値制御を行う。すなわち、しきい値制御器11Bでは、例えば、あらかじめ設定された大きさの異なる2つの係数rl,ru(0.0≦rl,ru≦1.0)を用いて、以下の(33)式((33A),(33B))のように、2種類の仮のしきい値Etmplo(x,y),Etmpup(x,y)を算出する。そして、それらを第1のしきい値Elo(x,y)および第2のEup(x,y)として、イプシロンフィルタ12Aに出力する。
【0133】
【数33】
【0134】
本実施の形態におけるイプシロンフィルタ12Aでは、しきい値制御器11Bで算出された2つのしきい値Elo(x,y),Eup(x,y)によりしきい値処理を行う。イプシロンフィルタ12Aでの処理は、より詳しくは、例えば、図16に示した構成により、以下のようにして行われる。差分器20によって算出された差分値D(x,y)は、本実施の形態においては絶対値算出器21のほかに符号判定器24にも送られる。
【0135】
符号判定器24では、差分値D(x,y)の符号を判定し、その結果をスイッチ27に送る。
【0136】
第1の比較器25では、しきい値制御器11Bから送られてくる第1のしきい値Elo(x,y)を用いて、第1の実施の形態と同様、(12)式により信号の選択を行う。すなわち、絶対値算出器21で算出された値AD(dx,dy)としきい値Elo(x,y)とを比較し、その結果に応じて、注目画素の値I(x,y)または近傍画素の値I(x+dx,y+dy)のいずれか一方を選択し、それを値J(dx,dy)として出力する。
【0137】
第2の比較器26では、しきい値制御器11Bから送られてくる第2のしきい値Eup(x,y)を用いて、第1の実施の形態と同様(12)式により信号の選択を行う。すなわち、絶対値算出器21で算出された値AD(dx,dy)としきい値Eup(x,y)とを比較し、その結果に応じて、注目画素の値I(x,y)または近傍画素の値I(x+dx,y+dy)のいずれか一方を選択し、それを値J(dx,dy)として出力する。
【0138】
スイッチ27では、符号判定器24の判定結果をもとに、第1の比較器25または第2の比較器26の出力のいずれか一方を選択し、線形ローパスフィルタ23に送る。スイッチ27では、例えば、符号判定結果が正であることを示している場合には、第1の比較器25の出力を選択する。逆に、判定結果が負である場合には、第2の比較器26の出力を選択する。
【0139】
本実施の形態では、近傍画素の値I(x+dx,y+dy)が現在の注目画素の値(近傍領域NBの中央の値)I(x,y)よりも大きい場合にはしきい値Eup(x,y)が用いられ、小さい場合にはElo(x,y)が用いられる。すなわち、イプシロンフィルタ12Aにおいて、高レベル側と低レベル側とで異なるしきい値を設定することが可能となっている。特に、第2のしきい値Eup(x,y)が第1のしきい値Elo(x,y)よりも大きくなるようにすることで、第3の実施の形態と同様、画素値の変化に含まれる照明成分の影響を軽減することができ、より適切に照明成分を抽出することが可能となる。
【0140】
[第5の実施の形態]
次に、本発明の第5の実施の形態について説明する。なお、以下の説明では、上記第1〜第4の実施の形態における構成要素と実質的に同一の機能を有する部分には同一の符号を付し、適宜説明を省略する。
【0141】
図17は、本発明の第5の実施の形態に係る画像処理装置の構成を示している。本実施の形態に係る画像処理装置の構成は、上記第1の実施の形態(図1)と類似しているが、特に、はじめに入力画像の画素レベルに例えば対数変換などの非線形な変換を施すようにした点が異なる。
【0142】
すなわち、本実施の形態においては、回路の入力段に対数変換器17を設置し、はじめに入力画像の各画素値I(x,y)に対して(1)式に示した対数変換を施すようになっている。またこれに伴い、イプシロンフィルタ12で得られた照明成分R(x,y)を差し引くための除算器13が減算器18に、圧縮された照明成分CR(x,y)を非照明成分S(x,y)に積算するための乗算器15が加算器19に置き換えられている。これらは積算、除算が対数変換後には加算、減算となるという良く知られた事実に基づくものである。
【0143】
対数変換器17によって対数変換された後の入力画像の画素値IL(x,y)は、エッジ強度算出器10、イプシロンフィルタ12、および減算器18に送られる。エッジ強度算出器10およびイプシロンフィルタ12では、その画素値IL(x,y)に基づいて、上記第1の実施の形態と同様の処理を行う。一方、減算器18では、各画素値IL(x,y)から、イプシロンフィルタ12で得られた照明成分R(x,y)を減算することにより、入力画像から照明成分を除去し、その結果得られる非照明成分S(x,y)を加算器19に出力する。加算器19では、非照明成分S(x,y)に補正照明成分CR(x,y)を加算することで画像信号を復元し、上記第1の実施の形態と同様、その結果の画像信号O(x,y)を図示しない伝送路、記憶装置または表示装置などに出力する。
【0144】
本実施の形態で用いた対数変換は、それ自体がダイナミックレンジの圧縮効果を持ち、入力レベルが高いほど画素値I(x,y)が大きく圧縮される。これにより、上記第3および第4の実施の形態と実質的に同様、画素値の空間的な変化における照明レベルの影響が軽減されることになり、より適切に照明成分を抽出することが可能となる。
【0145】
なお、本実施の形態では、非線形な変換として対数変換を行う例について説明したが、対数変換以外の他の非線形な変換を行うようにしても良い。
【0146】
【発明の効果】
以上説明したように、本発明の画像処理方法または画像処理装置によれば、入力画像上の各位置ごとに算出されたエッジ強度に基づいて、イプシロンフィルタで用いられるしきい値を制御すると共に、そのしきい値の制御されたイプシロンフィルタからの出力値に応じて、画素値を変換するための係数を算出し、その算出された係数により、各画素ごとに画素値の変換を行うようにしたので、複数の異なる照明が存在する場合にも、それらの境界を適切に抽出することが可能となり、不自然な画像パターンの発生を抑えて、主観的に好ましいダイナミックレンジの圧縮を実現できる。
【図面の簡単な説明】
【図1】本発明の第1の実施の形態に係る画像処理装置の構成を示すブロック図である。
【図2】画像の走査方向を示す説明図である。
【図3】レベル変換関数の例について示す説明図である。
【図4】図1に示した画像処理装置におけるイプシロンフィルタの構成を示すブロック図である。
【図5】従来法による再現画像と図1に示した画像処理装置による再現画像との違いについて示す説明図である。
【図6】レベル変換関数と係数算出関数との関係について示す説明図である。
【図7】イプシロンフィルタの効果について示す説明図である。
【図8】しきい値固定のイプシロンフィルタのエッジ周辺における振る舞いについて示す説明図である。
【図9】しきい値固定のイプシロンフィルタにおいて生ずるレベル勾配の逆転現象について示す説明図である。
【図10】イプシロンフィルタの出力に用いられるレベル変換関数とその微分値との関係を示す説明図である。
【図11】レベル変換曲線における入力レベルと圧縮率の大小関数について示す説明図である。
【図12】本発明の第2の実施の形態に係る画像処理装置の構成を示すブロック図である。
【図13】図12に示した画像処理装置における係数算出器において用いられる係数算出関数の例について示す説明図である。
【図14】本発明の第3の実施の形態に係る画像処理装置の構成を示すブロック図である。
【図15】本発明の第4の実施の形態に係る画像処理装置の構成を示すブロック図である。
【図16】図15に示した画像処理装置におけるイプシロンフィルタの構成を示す説明図である。
【図17】本発明の第5の実施の形態に係る画像処理装置の構成を示すブロック図である。
【図18】従来用いられているレベル変換関数の例について示す説明図である。
【図19】従来用いられている他のレベル変換関数の例について示す説明図である。
【図20】ヒストグラムイコライゼーションの原理について示す説明図である。
【図21】 Multiscale retinex法の問題点について説明するための図である。
【符号の説明】
10…エッジ強度算出器、11…しきい値制御器、12,12A…イプシロンフィルタ、13…除算器、14…レベル変換器、15…乗算器、16…係数算出器、17…対数変換器、18,20…差分器、19…加算器、21…絶対値算出器、22、25,26…比較器、23…線形ローパスフィルタ(LPF)、24…符号判定器、27…スイッチ。
Claims (10)
- 入力された画像を、相対的によりダイナミックレンジの小さい画像に変換するための画像処理方法であって、
入力画像上の各位置ごとにエッジ強度を算出するエッジ強度算出過程と、
算出された前記エッジ強度に基づいて、イプシロンフィルタのしきい値を制御するしきい値制御過程と、
前記しきい値制御過程において制御されたしきい値を用いて、前記入力画像に前記イプシロンフィルタによるフィルタ処理を施すフィルタリング過程と、
前記フィルタリング過程の出力値に応じて、画素値を変換するための係数を算出し、その算出した係数により、各画素ごとに画素値の変換を行う画素値変換過程と
を含み、
前記しきい値制御過程において、大きさの異なる2つのしきい値を算出し、
前記フィルタリング過程において、近傍画素の値が注目画素に比べて大きい場合と小さい場合とで異なるしきい値を用いる
画像処理方法。 - 前記エッジ強度算出過程では、前記エッジ強度として、注目画素の近傍領域内の画素値の1次微分値が大きいほど大きな値を算出する
請求項1記載の画像処理方法。 - 前記フィルタリング過程において、近傍画素の値が注目画素の値よりも大きい場合には、前記2つのしきい値のうち、値の大きい方のしきい値を用い、近傍画素の値が注目画素の値よりも小さい場合には、値の小さい方のしきい値を用いる
請求項1記載の画像処理方法。 - さらに、前記入力画像の画素レベルに、非線形な変換を施す過程を有し、
その非線形な変換が行われた後の入力画像に対して前記イプシロンフィルタによるフィルタ処理を施す
請求項1記載の画像処理方法。 - 前記非線形変換は対数変換である
請求項4記載の画像処理方法。 - 入力された画像を、相対的によりダイナミックレンジの小さい画像に変換するための画像処理装置であって、
入力画像上の各位置ごとにエッジ強度を算出するエッジ強度算出手段と、
設定されたしきい値を用いて、前記入力画像にフィルタ処理を施すイプシロンフィルタと、
前記エッジ強度算出手段によって算出された前記エッジ強度に基づいて、前記イプシロンフィルタで用いられるしきい値を制御するしきい値制御手段と、
前記イプシロンフィルタからの出力値に応じて、画素値を変換するための係数を算出し、その算出した係数により、各画素ごとに画素値の変換を行う画素値変換手段と
を備え、
前記しきい値制御手段は、大きさの異なる2つのしきい値を算出し、
前記イプシロンフィルタは、近傍画素の値が注目画素に比べて大きい場合と小さい場合とで異なるしきい値を用いるよう構成されている
画像処理装置。 - 前記エッジ強度算出手段は、前記エッジ強度として、注目画素の近傍領域内の画素値の1次微分値が大きいほど大きな値を算出するよう構成されている
請求項6記載の画像処理装置。 - 前記イプシロンフィルタは、近傍画素の値が注目画素の値よりも大きい場合には、前記2つのしきい値のうち、値の大きい方のしきい値を用い、近傍画素の値が注目画素の値よりも小さい場合には、値の小さい方のしきい値を用いるよう構成されている
請求項6記載の画像処理装置。 - さらに、前記入力画像の画素レベルに、非線形な変換を施す手段を有し、
前記イプシロンフィルタは、前記非線形な変換が行われた後の入力画像に対して前記フィルタ処理を施すよう構成されている
ことを特徴とする請求項6記載の画像処理装置。 - 前記非線形な変換を施す手段は、前記入力画像の画素レベルに対数変換を施すものである
ことを特徴とする請求項9記載の画像処理装置。
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001187146A JP4649781B2 (ja) | 2001-06-20 | 2001-06-20 | 画像処理方法および装置 |
US10/362,201 US7127122B2 (en) | 2001-06-20 | 2002-06-19 | Image processing method and image processing apparatus |
EP02743646A EP1404120A4 (en) | 2001-06-20 | 2002-06-19 | METHOD AND DEVICE FOR PROCESSING IMAGES |
PCT/JP2002/006106 WO2003001793A1 (fr) | 2001-06-20 | 2002-06-19 | Procédé et dispositif de traitement d'images |
KR1020037002412A KR100895207B1 (ko) | 2001-06-20 | 2002-06-19 | 화상 처리 방법 및 장치 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2001187146A JP4649781B2 (ja) | 2001-06-20 | 2001-06-20 | 画像処理方法および装置 |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2003008935A JP2003008935A (ja) | 2003-01-10 |
JP2003008935A5 JP2003008935A5 (ja) | 2008-04-24 |
JP4649781B2 true JP4649781B2 (ja) | 2011-03-16 |
Family
ID=19026482
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2001187146A Expired - Fee Related JP4649781B2 (ja) | 2001-06-20 | 2001-06-20 | 画像処理方法および装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US7127122B2 (ja) |
EP (1) | EP1404120A4 (ja) |
JP (1) | JP4649781B2 (ja) |
KR (1) | KR100895207B1 (ja) |
WO (1) | WO2003001793A1 (ja) |
Families Citing this family (60)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4045422B2 (ja) * | 2002-07-23 | 2008-02-13 | 株式会社日立メディコ | 画像処理装置 |
JP4329391B2 (ja) * | 2003-04-25 | 2009-09-09 | ソニー株式会社 | 信号処理装置および方法、記録媒体、並びにプログラム |
US7599578B2 (en) | 2004-09-17 | 2009-10-06 | Nikon Corporation | Apparatus, program, and method for image tone transformation, and electronic camera |
US7936941B2 (en) * | 2004-10-08 | 2011-05-03 | Panasonic Corporation | Apparatus for clearing an image and method thereof |
JP4127262B2 (ja) * | 2004-11-29 | 2008-07-30 | ソニー株式会社 | 画像処理装置及び画像処理方法 |
US7525579B2 (en) * | 2004-12-27 | 2009-04-28 | Konica Minolta Holdings, Inc. | Image sensing apparatus and image processing method for use therein |
US7961966B2 (en) * | 2005-01-04 | 2011-06-14 | Etron Technology, Inc. | Digitized image stabilization using energy analysis method |
JP4538358B2 (ja) * | 2005-03-30 | 2010-09-08 | パイオニア株式会社 | 画像処理装置 |
US8014034B2 (en) * | 2005-04-13 | 2011-09-06 | Acd Systems International Inc. | Image contrast enhancement |
JP2007011926A (ja) * | 2005-07-04 | 2007-01-18 | Fujifilm Holdings Corp | 画像処理方法および装置並びにプログラム |
JP4736598B2 (ja) * | 2005-07-27 | 2011-07-27 | ソニー株式会社 | ノイズ除去回路及びノイズ除去方法 |
US7791656B2 (en) | 2005-08-16 | 2010-09-07 | Konica Minolta Holdings, Inc. | Image sensing apparatus and image processing method |
JP4736939B2 (ja) * | 2005-08-16 | 2011-07-27 | コニカミノルタホールディングス株式会社 | 撮像装置及び画像処理方法 |
JP4135737B2 (ja) | 2005-08-30 | 2008-08-20 | ブラザー工業株式会社 | 画像処理プログラム |
US7590303B2 (en) * | 2005-09-29 | 2009-09-15 | Samsung Electronics Co., Ltd. | Image enhancement method using local illumination correction |
CN101288099B (zh) | 2005-10-12 | 2012-11-21 | 松下电器产业株式会社 | 视觉处理装置及方法 |
US7755670B2 (en) * | 2005-10-24 | 2010-07-13 | Nikon Corporation | Tone-conversion device for image, program, electronic camera, and tone-conversion method |
US7995058B2 (en) * | 2006-01-27 | 2011-08-09 | Tandent Vision Science, Inc. | Method and system for identifying illumination fields in an image |
US7843493B2 (en) * | 2006-01-31 | 2010-11-30 | Konica Minolta Holdings, Inc. | Image sensing apparatus and image processing method |
JP4977395B2 (ja) | 2006-04-14 | 2012-07-18 | 富士フイルム株式会社 | 画像処理装置及び方法 |
CN101427559B (zh) | 2006-04-19 | 2011-11-30 | 松下电器产业株式会社 | 视觉处理装置、视觉处理方法、显示装置和集成电路 |
US7903898B2 (en) | 2006-04-28 | 2011-03-08 | Panasonic Corporation | Visual processing apparatus, visual processing method, program, recording medium, display device, and integrated circuit |
WO2008020487A1 (fr) * | 2006-08-17 | 2008-02-21 | Fujitsu Limited | dispositif de traitement d'image, programme DE TRAITEMENT D'IMAGE, et procédé de traitement d'image |
JP4208909B2 (ja) | 2006-08-24 | 2009-01-14 | 株式会社東芝 | 画像処理装置と撮影装置 |
JP4859632B2 (ja) | 2006-11-15 | 2012-01-25 | 富士通セミコンダクター株式会社 | 画像処理装置及び画像処理方法 |
JP4288623B2 (ja) * | 2007-01-18 | 2009-07-01 | ソニー株式会社 | 撮像装置、ノイズ除去装置、ノイズ除去方法、ノイズ除去方法のプログラム及びノイズ除去方法のプログラムを記録した記録媒体 |
US8396290B1 (en) | 2007-03-09 | 2013-03-12 | Pixar | Saturation varying color space |
US8942476B1 (en) | 2007-03-09 | 2015-01-27 | Pixar | Saturation varying and lighting independent color color control for computer graphics |
JP4746575B2 (ja) | 2007-03-16 | 2011-08-10 | 株式会社東芝 | 画像処理装置及びその方法 |
JP4772721B2 (ja) | 2007-03-26 | 2011-09-14 | 株式会社東芝 | 画像処理装置及びその方法 |
JP4527750B2 (ja) * | 2007-05-30 | 2010-08-18 | 三菱電機株式会社 | 画像処理装置及び方法並びに画像表示装置 |
KR100901354B1 (ko) * | 2007-05-31 | 2009-06-05 | 주식회사 코아로직 | 영상의 에지 보정장치 및 그 방법 |
CN102325257A (zh) * | 2007-07-20 | 2012-01-18 | 富士胶片株式会社 | 图像处理设备和图像处理方法 |
US8330768B2 (en) * | 2007-07-27 | 2012-12-11 | Sharp Laboratories Of America, Inc. | Apparatus and method for rendering high dynamic range images for standard dynamic range display |
JP5012333B2 (ja) | 2007-08-30 | 2012-08-29 | コニカミノルタアドバンストレイヤー株式会社 | 画像処理装置および画像処理方法ならびに撮像装置 |
US8300974B2 (en) | 2007-11-12 | 2012-10-30 | Mitsubishi Electric Corporation | Image processing device, image display device, and image processing method therefor |
JP4502055B2 (ja) * | 2007-11-12 | 2010-07-14 | 三菱電機株式会社 | 画像処理装置、画像表示装置、および画像処理方法 |
KR20100099686A (ko) | 2007-12-04 | 2010-09-13 | 소니 주식회사 | 화상 처리 장치 및 방법, 프로그램, 및 기록 매체 |
US7995854B2 (en) * | 2008-03-28 | 2011-08-09 | Tandent Vision Science, Inc. | System and method for identifying complex tokens in an image |
US8139883B2 (en) * | 2008-07-29 | 2012-03-20 | Sony Corporation | System and method for image and video encoding artifacts reduction and quality improvement |
EP2166509A1 (en) * | 2008-09-22 | 2010-03-24 | Solomon Systech Limited | Method and apparatus of local contrast enhancement |
US8103120B2 (en) * | 2008-09-22 | 2012-01-24 | Solomon Systech Limited | Method and apparatus of local contrast enhancement |
KR100975221B1 (ko) * | 2008-11-05 | 2010-08-10 | 매그나칩 반도체 유한회사 | 샤프니스 보정장치 및 그 방법 |
GB2468304A (en) * | 2009-03-03 | 2010-09-08 | Sharp Kk | Video image processing method that smooths the jagged diagonal edges created by epsilon filtering |
JP5233757B2 (ja) * | 2009-03-09 | 2013-07-10 | 三菱電機株式会社 | 画像処理装置および画像表示装置 |
JP4585602B1 (ja) * | 2009-09-18 | 2010-11-24 | 株式会社東芝 | 画像処理装置、表示装置及び画像処理方法 |
JP2011090569A (ja) * | 2009-10-23 | 2011-05-06 | Sony Corp | 画像処理装置および画像処理方法 |
CN101930599B (zh) * | 2010-08-24 | 2013-02-13 | 黄伟萍 | 一种医学图像增强方法及*** |
DE102010039807A1 (de) * | 2010-08-26 | 2012-03-01 | Siemens Aktiengesellschaft | Vorrichtung und Verfahren zur Filterung in Röntgenbildern |
JP5470415B2 (ja) * | 2012-03-30 | 2014-04-16 | Eizo株式会社 | イプシロンフィルタの閾値決定方法およびローパスフィルタの係数決定方法 |
JP5267704B2 (ja) * | 2012-06-06 | 2013-08-21 | コニカミノルタアドバンストレイヤー株式会社 | 画像処理装置および画像処理方法ならびに撮像装置 |
JP2013017240A (ja) * | 2012-10-10 | 2013-01-24 | Konica Minolta Advanced Layers Inc | 画像処理装置、画像処理方法、及び撮像装置 |
JP6625165B2 (ja) | 2013-02-14 | 2019-12-25 | キヤノン株式会社 | 画像処理装置、画像処理方法、およびプログラム |
JP6362333B2 (ja) | 2013-02-14 | 2018-07-25 | キヤノン株式会社 | 画像処理装置、画像処理方法、およびプログラム |
JP6303816B2 (ja) * | 2014-05-28 | 2018-04-04 | 富士ゼロックス株式会社 | 画像処理装置及びプログラム |
JP2016086347A (ja) | 2014-10-28 | 2016-05-19 | 三星ディスプレイ株式會社Samsung Display Co.,Ltd. | 画像処理装置、画像処理方法、及びプログラム |
JP6548907B2 (ja) | 2015-02-24 | 2019-07-24 | 三星ディスプレイ株式會社Samsung Display Co.,Ltd. | 画像処理装置、画像処理方法、及びプログラム |
JP6548403B2 (ja) | 2015-02-24 | 2019-07-24 | 三星ディスプレイ株式會社Samsung Display Co.,Ltd. | 画像処理装置、画像処理方法、及びプログラム |
KR102495754B1 (ko) | 2016-10-25 | 2023-02-02 | 한화테크윈 주식회사 | 광역 영상의 동적 범위를 압축하는 영상 처리 장치 및 방법 |
JP6742231B2 (ja) * | 2016-12-09 | 2020-08-19 | キヤノン株式会社 | 画像処理装置及び方法、及び撮像装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62155679A (ja) * | 1985-12-27 | 1987-07-10 | Olympus Optical Co Ltd | 画像信号処理回路 |
JPH07212759A (ja) * | 1994-01-27 | 1995-08-11 | Sharp Corp | 画像信号復号器 |
JPH0865673A (ja) * | 1994-08-22 | 1996-03-08 | Kyocera Corp | 画像信号復号装置 |
JPH09247674A (ja) * | 1996-03-11 | 1997-09-19 | Matsushita Electric Ind Co Ltd | 輪郭および輪郭近傍歪み除去装置 |
JPH11266358A (ja) * | 1998-03-17 | 1999-09-28 | Fuji Photo Film Co Ltd | 画像処理方法 |
JP2001034748A (ja) * | 1999-07-26 | 2001-02-09 | Nippon Telegr & Teleph Corp <Ntt> | 画像補正方法及びその装置及びその方法を記録した記録媒体及びその装置を組み込んだ画像撮影装置及びその装置を組み込んだ画像表示装置 |
Family Cites Families (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4527265A (en) * | 1983-06-30 | 1985-07-02 | Staar S. A. | Position memory for laser reproduction unit |
BE902324A (fr) * | 1985-04-30 | 1985-08-16 | Staar Sa | Appareil assurant la reproduction et/ou l'enregistrement selectif et/ou successif de supports d'informations. |
JPS62214576A (ja) * | 1986-03-14 | 1987-09-21 | Pioneer Electronic Corp | ピツクアツプの位置制御方法 |
US4878129A (en) * | 1986-08-27 | 1989-10-31 | Pioneer Electronic Corporation | Method for reproduction of stored video signals with disk reproducing apparatus |
JP2609542B2 (ja) * | 1988-10-04 | 1997-05-14 | パイオニア株式会社 | ディスクプレーヤ |
CA2020059C (en) * | 1989-06-30 | 1995-05-23 | Shigemi Maeda | Modular ceramic igniter |
JPH04254959A (ja) * | 1991-02-05 | 1992-09-10 | Sony Corp | 記録媒体の再生装置 |
JP2970008B2 (ja) * | 1991-02-22 | 1999-11-02 | ソニー株式会社 | ディスク状記録媒体再生装置 |
JP3396894B2 (ja) * | 1992-05-14 | 2003-04-14 | ソニー株式会社 | 音響再生装置 |
JPH07326175A (ja) * | 1994-05-30 | 1995-12-12 | Sony Corp | 記録又は再生装置 |
JPH08147702A (ja) * | 1994-11-11 | 1996-06-07 | Mitsumi Electric Co Ltd | 光ディスク書き込み方法 |
JPH08241565A (ja) * | 1995-03-03 | 1996-09-17 | Toshiba Corp | ポータブルコンピュータ |
JPH10508138A (ja) * | 1995-08-22 | 1998-08-04 | フィリップス エレクトロニクス ネムローゼ フェンノートシャップ | 媒体のオーディオ信号及び/又はビデオ信号読取り用プレーヤ |
US6463601B1 (en) * | 1995-10-12 | 2002-10-15 | Randy E. Fetty | Self-contained waste and refuse disposal unit |
KR100279147B1 (ko) * | 1996-07-27 | 2001-01-15 | 구자홍 | 에지강조 오차확산법의 에지검출 방법 |
KR100283948B1 (ko) * | 1997-09-12 | 2001-05-02 | 윤종용 | 디스크플레이어제어장치및그방법 |
JP3873523B2 (ja) * | 1999-05-21 | 2007-01-24 | ソニー株式会社 | 再生装置 |
US6724943B2 (en) * | 2000-02-07 | 2004-04-20 | Sony Corporation | Device and method for image processing |
JP4556276B2 (ja) * | 2000-03-23 | 2010-10-06 | ソニー株式会社 | 画像処理回路及び画像処理方法 |
US6396777B1 (en) * | 2000-05-30 | 2002-05-28 | Delphi Technologies, Inc. | Compact disk player instant-on feature |
-
2001
- 2001-06-20 JP JP2001187146A patent/JP4649781B2/ja not_active Expired - Fee Related
-
2002
- 2002-06-19 EP EP02743646A patent/EP1404120A4/en not_active Withdrawn
- 2002-06-19 WO PCT/JP2002/006106 patent/WO2003001793A1/ja active Application Filing
- 2002-06-19 KR KR1020037002412A patent/KR100895207B1/ko not_active IP Right Cessation
- 2002-06-19 US US10/362,201 patent/US7127122B2/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62155679A (ja) * | 1985-12-27 | 1987-07-10 | Olympus Optical Co Ltd | 画像信号処理回路 |
JPH07212759A (ja) * | 1994-01-27 | 1995-08-11 | Sharp Corp | 画像信号復号器 |
JPH0865673A (ja) * | 1994-08-22 | 1996-03-08 | Kyocera Corp | 画像信号復号装置 |
JPH09247674A (ja) * | 1996-03-11 | 1997-09-19 | Matsushita Electric Ind Co Ltd | 輪郭および輪郭近傍歪み除去装置 |
JPH11266358A (ja) * | 1998-03-17 | 1999-09-28 | Fuji Photo Film Co Ltd | 画像処理方法 |
JP2001034748A (ja) * | 1999-07-26 | 2001-02-09 | Nippon Telegr & Teleph Corp <Ntt> | 画像補正方法及びその装置及びその方法を記録した記録媒体及びその装置を組み込んだ画像撮影装置及びその装置を組み込んだ画像表示装置 |
Also Published As
Publication number | Publication date |
---|---|
JP2003008935A (ja) | 2003-01-10 |
EP1404120A4 (en) | 2006-10-04 |
EP1404120A1 (en) | 2004-03-31 |
US20030156761A1 (en) | 2003-08-21 |
WO2003001793A1 (fr) | 2003-01-03 |
US7127122B2 (en) | 2006-10-24 |
KR100895207B1 (ko) | 2009-05-04 |
KR20030025294A (ko) | 2003-03-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4649781B2 (ja) | 画像処理方法および装置 | |
EP2833317B1 (en) | Image display device and/or method therefor | |
US7570390B2 (en) | Image processing device and method | |
US8654221B2 (en) | Image processing device and method, and program | |
JP4556276B2 (ja) | 画像処理回路及び画像処理方法 | |
US8248492B2 (en) | Edge preserving and tone correcting image processing apparatus and method | |
US8849056B2 (en) | Apparatus and method for image edge-preserving smoothing | |
US7421120B2 (en) | Apparatus correcting image by luminance histogram | |
JP4858610B2 (ja) | 画像処理方法 | |
EP1209621A2 (en) | A method for enhancing a digital image based upon pixel color | |
US6724942B1 (en) | Image processing method and system | |
US8831346B2 (en) | Image processing apparatus and method, and program | |
JP2004147333A (ja) | 選択的空間フィルタを使用するディジタル画像の階調及び空間特性の向上方法 | |
EP1111906A2 (en) | A method for enhancing the edge contrast of a digital image independently from the texture | |
KR101031153B1 (ko) | 신호 처리 장치 및 방법 | |
KR100768035B1 (ko) | 적응적 화질 개선 장치 및 그 방법 | |
JP5295854B2 (ja) | 画像処理装置及び画像処理プログラム | |
JP2007215188A (ja) | 画像変換方法、変換画像生成方法、および画像補正装置 | |
JP4277550B2 (ja) | 信号処理装置および方法、記録媒体、並びにプログラム | |
EP1622080B1 (en) | Signal processing device and method, recording medium, and program | |
US8665315B2 (en) | Signal processing device, reproducing device, signal processing method and program | |
JP2003009178A (ja) | 画像処理方法および装置 | |
JP2004289606A (ja) | 信号処理装置および方法、記録媒体、並びにプログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080310 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080310 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100928 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20101029 |
|
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: 20101116 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20101129 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131224 Year of fee payment: 3 |
|
LAPS | Cancellation because of no payment of annual fees |