JP4339611B2

JP4339611B2 - コントラスト強調方法及びコントラスト強調装置

Info

Publication number: JP4339611B2
Application number: JP2003049952A
Authority: JP
Inventors: レイ、チューチュン; サルトル、ペイエルジョルジオ; イネスコ、ヤルシン; ワグナー、ペーター
Original assignee: ソニーインターナショナル（ヨーロッパ）ゲゼルシャフトミットベシュレンクテルハフツング
Priority date: 2002-02-26
Filing date: 2003-02-26
Publication date: 2009-10-07
Anticipated expiration: 2023-02-26
Also published as: JP2003281527A; US7102697B2; US20030161549A1; KR20030070845A; EP1345172A1

Description

【０００１】
【発明の属する技術分野】
本発明は、デジタル画像のコントラストを強調するコントラスト強調方法に関し、特にテレビジョン画像のコントラストを強調する方法に関する。更に、本発明は、デジタル画像のコントラストを向上させるコントラスト強調装置に関する。
【０００２】
【従来の技術】
ヒストグラム均等化法（histogram equalization）は、単純且つ効果的な技術であるため、デジタル画像のコントラストを向上させる処理に広く用いられている。ヒストグラム均等化法は、包括的手法（global approach）と局所的手法（local approach）といった２つのカテゴリに分類される。尚、ここでは、演算負荷を考慮して、包括的ヒストグラム均等化法が好ましいと判断し、特にテレビジョン画像のコントラストの強調に包括的ヒストグラム均等化法を用いる具体例のみを説明する。
【０００３】
ヒストグラム均等化法は、画像情報を最大限知覚（maximal Image information perception）するには、画像の輝度値の分布（ヒストグラム）が画像の輝度値の範囲全体に亘って均一でなくてはならないという原理に基づいている。ヒストグラム均等化変換曲線（histogram equalization transformation curve）によって画素の輝度値をマッピングすることにより、変換画像の画素の輝度値は、ダイナミックレンジ全体に亘って可能な限り均等に分布する。この「ヒストグラム均等化変換曲線」を以下「ＨＥ変換曲線」と呼ぶ。
【０００４】
ＨＥ変換曲線は、オリジナル画像のヒストグラムから得ることができる。適切なＨＥ変換曲線を生成するために、ヒストグラム値を累算することにより、累積ヒストグラム（cumulative histogram）、すなわち累積分布関数（Cumulative Distribution Function：以下、ＣＤＦともいう。）が算出される。そして、累積ヒストグラムを正規化し、輝度レベルをマッピングするために最適な変換関数を求める。
【０００５】
ここで、多くの用途において、ヒストグラム均等化に基づく画像の性質の変化の度合いが好ましくない場合がある。すなわち、ヒストグラム均等化を行うと、画像の明度（brightness）が大きく変化しすぎることがある。例えば、夜間に撮像された画像が昼間に撮影された画像のように見えてしまうこともある。また、人間の顔に対応する領域における自然な印象が損なわれることもある。
【０００６】
【発明が解決しようとする課題】
欧州特許ＥＰ３８３２６９Ａ２号「部分的画像適応化法（Regionally adaptive Imaging techniques）」には、画像雑音を強調しすぎることを回避するためのヒストグラムフィルタリング法が開示されている。
【０００７】
また、ワイ・ティー・キム（Y.-T.Kim.）に付与されている米国特許第５８５７０３３号「定量化された平均分離ヒストグラム均等化法を用いた画像強調方法及び回路（Method for image enhancement using quantified mean separate histogram equalization and a circuit thereof）」では、画像の平均値に基づいて画像を２つのサブ画像に分割することにより、明度が極端に変化しないようにしている。ここでは、各サブ画像は、個々にヒストグラム均等化される。
【０００８】
また、ワイ・ティー・キム（Y.-T.Kim.）に付与されている米国特許第５９２３３８３号「ヒストグラム均等化を用いた画像強調方法（Image enhancement method using histogram equalization）」では、明度の極端な変化を回避する更なる手法が開示されている。ここでは、ヒストグラムを算出する際に、各輝度値の出現数が所定数を超えないように制限される。しかしながら、この文献では、異なる画像に対して、この所定数をどのように決定するかは開示されていない。更に、欧州特許ＥＰ３８３２６９Ａ２号にも指摘されているように、この種の制限の結果、輝度値の範囲が小さくなり、画像のコントラストが劣化するという問題もある。
【０００９】
コーネリアス・エー（Cornelius A.）、エム・ヤスパース（M.Jaspers）に付与されている米国特許第５５３７０７１号「入力信号を振幅セグメント信号に分割し、各振幅セグメント信号の値に基づいて各振幅セグメント信号を非線形的に処理する方法及び回路（Method and circuit for dividing an input signal into amplitude segment signals and for non-linearity processing the amplitude segment signals on the basis of the value of each amplitude segment signal）」に開示される手法によれば、ヒストグラム均等化における問題を改善することができる。この手法では、ヒストグラムを複数（例えば３２個）のセグメントに分割し、これらセグメントのヒストグラムエントリからこれらの平均値を減算する。続いて、これらのセグメントに亘る積分を行う。そして、係数による重み付けを行い、積分結果を表す曲線と単位直線（unit straight line）とを混合することにより、画像コントラストが強調される。しかしながら、この手法では、実現されるヒストグラム均等化の度合いが限られている。特に、コントラストが弱い画像に対しては、強調量が、画像の性質により許容される望ましい強調量に比べて遙かに少ない。
【００１０】
欧州特許ＥＰ１０２２６７９Ａ２号には、「画像処理システム（Image processing system）」が開示されている。ここに開示される手法では、ヒストグラムは複数のクラスタに分割され、各クラスタに対してヒストグラム均等化又は伸張（stretching）が実行される。このヒストグラムの分割は、パターンマッチング法（pattern matching technique）に基づいて行われ、均一なガウス分布（Gaussian and uniform distribution）は、個々のクラスタに分割される。しかしながら、ヒストグラムパターンは、必ずしも均一なガウス分布を有するとは限らず、他の種類の分布パターンも存在する。
【００１１】
そこで、本発明は上述の課題に鑑みてなされたものであり、本発明の目的は、アーティファクトの発生を回避し、変更画像がより自然に見える、ヒストグラム均等化を用いたデジタル画像のコントラスト強調方法及びコントラスト強調装置を提供することである。
【００１２】
【課題を解決するための手段】
上述の課題を解決するために、本発明に係るコントラスト強調方法は、それぞれに輝度レベルが割り当てられた複数の画素を含むデジタル画像のコントラストを強調するコントラスト強調方法において、オリジナル画像の輝度レベルの分布を示すヒストグラムを判定するステップと、オリジナル画像のコントラストの特性測定値を判定するステップと、オリジナル画像のヒストグラムを均等化し、コントラストを向上させるステップとを有する。ヒストグラムの均等化の度合いは、特性測定値の関数として変更される。
【００１３】
従来の手法のように全体的なヒストグラムと均等化を行う代わりに、本発明では、部分的なヒストグラム均等化を実現することができる。ヒストグラム均等化の度合いは、ヒストグラム均等化を行わない処理から全体的なヒストグラム均等化を行う処理までの間で変更され、これにより画素の輝度値は、そのダイナミックレンジ全体に亘って可能な限り均等に分布される。本発明では、ヒストグラム均等化の度合いは、自動的に画像のコントラストに適応化される。オリジナル画像を分析することにより、画像コントラストの特性測定値が算出される。
【００１４】
画像のコントラストが高い場合、ヒストグラム均等化は不要であり、この場合、ヒストグラム均等化をオフにできる。一方、画像のコントラストが低い場合、ヒストグラム均等化をオンにし、全体的なヒストグラム均等化を行うことができる。画像コントラストが中間的である場合、適正量のヒストグラム均等化が実行される。
【００１５】
この方法の利点は、オリジナル画像の質が高い場合、そのオリジナル画像にはヒストグラム均等化が施されず、したがって質が高い画像を劣化させることがない点にある。更に、本発明では、コントラストの強調が適量に制限されるため、例えば人間の顔が不自然に見えてしまうような、ヒストグラム均等化に起因するアーチファクトの発生を回避することができる。すなわち、本発明により、ヒストグラム均等化は必要最小限に抑制される。したがって、ヒストグラム均等化後の画像の全体的な画質が向上する。
【００１６】
ヒストグラムは、好ましくは、変更輝度レベルを算出するために、ヒストグラム均等化変換曲線に基づいて、オリジナル画像の輝度レベルを変換することにより、均等化される。ヒストグラム均等化変換曲線は、オリジナル画像の輝度レベルを変更し、コントラストが強調された画像を生成するために使用される。ヒストグラム均等化関数は、輝度値の範囲内の初期の各輝度レベルを対応するターゲット輝度レベルにマッピングするマッピング曲線とみなすことができる。ターゲット輝度レベルは、画素毎に初期の輝度レベルを対応するターゲット輝度レベルにマッピングすることにより、画像のコントラストが高まるように選択される。
【００１７】
ヒストグラム均等化変換曲線は、好ましくは、ヒストグラムを積分及び換算することにより生成される。オリジナル画像の輝度レベルの分布を示すオリジナルヒストグラムの値を累積することにより、累積分布関数（Cumulative Distribution Function：以下、ＣＤＦという。）が得られる。累積分布関数のターゲット値の範囲は、輝度レベルの範囲に対応していないため、累積分布関数に対する換算を行う必要がある。累積分布関数を換算することにより、全体的なヒストグラム均等化に適切に用いることができるヒストグラム均等化変換曲線が得られる。
【００１８】
ここで、特性測定値は、オリジナル画像の輝度レベルの平均値と所定の定数との差を評価することにより算出され、この差が大きいほど画像のコントラストが低いと判定する。通常、低いコントラストは、照明が不適切であるという条件に起因し、コントラストが低い画像は、明るすぎるか暗すぎるかのいずれかである。このため、コントラストが低い画像の平均輝度レベルは、通常の画像の平均輝度レベルに比べて著しく小さい又は大きい。この平均輝度値と、例えば輝度値が０〜２５５の範囲を取る場合は１２８といった定数との差は、画像コントラストの特性測定値として適しており、また、この値は容易に算出することができる。
【００１９】
本発明の他の好ましい具体例においては、特性測定値は、輝度レベル０から最初の主要ヒストグラムピークまでの重み付けされた距離と、最後の主要ヒストグラムピークから最大の輝度レベルまでの重み付けされた距離との差を評価することにより算出され、この重み付けされた距離の差が大きい程画像のコントラストが低いと判定する。このような判断基準を設けることにより、画像全体が暗すぎるか、明るすぎるかを判定することができる。ヒストグラムのピークが輝度レベル範囲の中間にない場合、この画像はコントラストが低い画像であると知覚される。この場合、強いヒストグラム均等化が適切である。
【００２０】
また、好ましくは、特性測定値の算出において主要ヒストグラムピーク間の空き領域の総和を考慮し、この主要ヒストグラムピーク間の空き領域の総和が大きい程画像のコントラストが高いと判定する。画像内のどの画素にも割り当てられていない全ての輝度値は、この空き領域に対応する。ヒストグラムの空き領域は、画像の特徴を区別するのに有用である。すなわち、ヒストグラムのピーク間に大きな空き領域を有する画像は、コントラストが高い画像であり、したがって、比較的弱いヒストグラム均等化を行えばよい。このような追加的な基準を考慮に入れることにより、ヒストグラム均等化の度合いを少なく維持できる。
【００２１】
ここで、主要ヒストグラムピークは、好ましくは、ヒストグラム又はヒストグラム均等化変換曲線の傾きが所定の閾値をどこで超えているかを判断することにより判定される。この傾きは、隣接するヒストグラム値を減算することにより容易に求められる。この減算の結果が所定の閾値を超える場合、主要ヒストグラムピークが検出される。
【００２２】
更に、好ましくは、オリジナル画像に対してエッジ検出処理を施すことによりエッジ画像を生成し、特性測定値を算出する際、オリジナル画像のコントラストと、エッジ画像のコントラストの両方を考慮する。ここでは、オリジナル画像とエッジ画像の両方についてヒストグラムを生成する。エッジ画像については、エッジ周辺の画素のみを用いてヒストグラムを生成する。したがって、画像のエッジ程多くの情報を含まない、例えば背景等の均質な領域内の画素は、ヒストグラムに大きな影響を与えない。ここでは、オリジナル画像の特性測定値とエッジ画像の特性測定値といった２つの異なる特性測定値が得られる。この場合、高いコントラストを示している方の特性測定値に基づいてヒストグラム均等化の度合いを制御するとよい。これにより、画像内の均質な部分が大きい場合にヒストグラム均等化の度合いが強くなってしまうといった問題を回避できる。
【００２３】
好ましくは、ヒストグラムに対しローパスフィルタリング処理を施した後にヒストグラム均等化変換曲線を生成する。ローパスフィルタリングにより、ヒストグラムの高周波数成分が除去される。すなわち、ローパスフィルタのカットオフ周波数以下の周波数成分及び特にＤＣ成分のみがフィルタを通過する。これにより、ヒストグラムが平滑化され、ＤＣ成分及び低い周波数成分が強調される。この平滑化されたヒストグラムを積分することにより、ヒストグラム均等化変換曲線が生成される。これにより、ヒストグラムの高周波変動に起因するヒストグラム均等化変換曲線の急激な変化を回避することができる。更に、ローパスフィルタリングにより、ヒストグラムからの外乱雑音（noise disturbance）を除去でき、したがって、ヒストグラム均等化変換曲線がこれらの外乱雑音によって歪むことがなくなる。ヒストグラムの周波数スペクトルをカットオフ周波数以下の周波数成分に制限することにより、コントラストを過剰に強調してしまうことが回避される。
【００２４】
本発明の第１の具体例では、ヒストグラム均等化の度合いは、ローパスフィルタリング処理のカットオフ周波数を変化させることにより変更され、カットオフ周波数が高い程、ヒストグラム均等化の度合いが強くなる。
【００２５】
ここで、ローパスフィルタのカットオフ周波数を比較的高い周波数に設定したとする。この場合、ヒストグラム内の多くの周波数成分がフィルタを通過する。すなわち、高い周波数成分は大きくは減衰されず、ヒストグラムは余り変化しない。ヒストグラム均等化変換曲線は、高周波成分を含み、したがって比較的粗く見えるこのヒストグラムを積分し換算することにより生成される。生成されたヒストグラム均等化変換関数は、オリジナル画像の輝度レベルに大きな影響を与える。このように、ローパスフィルタのカットオフ周波数を高くすると、略全体的なヒストグラム均等化が行われる。
【００２６】
ここで、逆にローパスフィルタのカットオフ周波数を比較的低い周波数に設定したとする。この場合、主にヒストグラムのＤＣ成分がフィルタを通過する。ＤＣ成分を積分し換算すると、傾きが１の線形の変換曲線が得られる。このような傾きが１の線形の変換曲線をヒストグラム均等化変換曲線として用いると、オリジナル画像の輝度レベルは全く変化しない。すなわち、カットオフ周波数を非常に低くすれば、オリジナル画像は全く変化しない。したがって、ヒストグラムの高周波成分を強く減衰させることにより、ヒストグラム均等化の度合いを低減することができる。
【００２７】
特性測定値が高い画像コントラストを示している場合、カットオフ周波数はより低い周波数にシフトされ、ヒストグラムはヒストグラム均等化変換曲線を生成する前に大きく減衰され、特性測定値が低い画像コントラストを示している場合、カットオフ周波数はより高い周波数にシフトされ、ヒストグラムはヒストグラム均等化変換曲線を生成する前に小さく減衰される。本発明においては、高い特性測定値は、画像のコントラストが低いことを示している。したがって、オリジナル画像のコントラストの関数としてヒストグラム均等化の度合いを自動的に制御するためには、上述のカットオフ周波数と特性測定値との直接的な対応関係を確立すればよいことになる。
【００２８】
ここで、ローパスフィルタリング処理を実行するために、例えば、フーリエ変換、好ましくは高速フーリエ変換を実行してローパスフィルタリング処理を周波数領域で行ってもよい。オリジナルヒストグラムをフーリエ変換することにより、周波数領域におけるヒストグラムのフーリエスペクトルが得られる。このフーリエスペクトルが既知になると、カットオフ周波数以上の全ての周波数成分を減衰又は除去することによりローパスフィルタリングを実行できる。この手法では、周波数領域においてフィルタ特性特にカットオフ周波数を変更できる。このようなローパスフィルタリングに代えて、例えば順方向及び逆方向ＩＩＲフィルタを用いて、順方向及び逆方向フィルタリングを実行してもよい。
【００２９】
ローパスフィルタリング手段は、好ましくは、フーリエ変換されたローパスフィルタ特性が登録されているルックアップテーブルを備える。ルックアップテーブルには、特定の平滑化フィルタによるフーリエ変換の結果が登録される。これにより、例えば、ルックアップテーブルに登録されているフーリエ変換の結果を右にシフトするだけでカットオフ周波数を高めることができる。一方、ルックアップテーブルに登録されているフーリエ変換の結果を右にシフトすれば、カットオフ周波数を低めることができる。
【００３０】
本発明の第２の具体例においては、ヒストグラム均等化の度合いは、ヒストグラム均等化変換曲線を生成する際、第１のカットオフ周波数に対応する第１のローパスフィルタリングされたヒストグラムと、第２のカットオフ周波数に対応する第２のローパスフィルタリングされたヒストグラムとの間を補間することにより変更される。ここでは、固定された２つのカットオフ周波数を有する２つのローパスフィルタリングが行われ、これらの処理を並列に行ってもよい。この第２の具体例によれば、カットオフ周波数をシフトさせる必要がなくなる。
【００３１】
本発明の第３の具体例では、ヒストグラム均等化の度合いは、ヒストグラム均等化変換曲線と、傾きが１の線形の変換曲線との差をゲイン係数により換算することにより変更され、ゲイン係数が小さいほどヒストグラム均等化の度合いが弱くなる。この具体例では、ヒストグラム均等化変換曲線が直接変更される。ここでは、ヒストグラム又はヒストグラム均等化変換曲線を周波数領域に変換する必要がない。この手法では、ヒストグラム均等化変換曲線から傾きが１の線形の変換曲線を減算し、求められた差を換算し、換算された曲線を傾きが１の線形の変換曲線に加算するだけでよく、演算が単純である。この単純な手法により、値が０〜１の範囲をとるゲイン係数を用いて、ヒストグラム均等化の度合いを連続的に変化させることができる。
【００３２】
明るさが変化してしまう問題を解決する第１の手法では、ヒストグラム均等化の前後の画像の平均輝度レベルの比に基づいて、ヒストグラム均等化後の画像の輝度レベルを換算する。この平均輝度レベルの比は、画像全体に対して１回のみ算出すればよい。
【００３３】
明るさが変化してしまう問題を解決する第２の手法では、ヒストグラム均等化の前後の画像の平均輝度レベルの差に基づいて、ヒストグラム均等化後の画像の輝度レベルをシフトする。この第２の手法では、輝度レベルの適応化に加算及び減算のみを用いるため、演算負荷を著しく軽減できる。
【００３４】
明るさが変化してしまう問題を解決する第３の手法では、ヒストグラム均等化により生じたヒストグラムピーク又は選択された基準点の輝度レベルのシフトに基づいて、ヒストグラム均等化後の画像の輝度レベルをシフトする。ヒストグラムの主要ピークを追跡することにより、明るさの変化の度合いを容易に検出することができる。
【００３５】
均等化変換曲線は、好ましくは、輝度レベルの任意の範囲について、均等化変換曲線を傾きが１の線形の変換曲線に置換することにより均等化変換曲線を変更する。例えば大部分の領域が暗く、他の小さな領域が明るい画像に対しては、包括的なヒストグラム均等化によってもこの種の画像の暗い領域のコントラストを向上させることができるが、このとき、明るい領域の細部がある程度失われる可能性がある。したがって、この場合、暗い領域のコントラストをのみを向上させ、明るい領域の画素は、変更しない方がよい。これは、ヒストグラム均等化変換曲線のうち、明るい領域に対応する部分を傾きが１の線形の変換曲線に置換することにより実現される。これにより、明るい領域において細部が失われることが防止される。また、大部分が明るく残りの僅かな部分が暗い画像に対してもこれと同様の処理を施すことができる。
【００３６】
【発明の実施の形態】
以下、本発明に係るコントラスト強調方法及びコントラスト強調装置について、図面を参照して詳細に説明する。
【００３７】
図１Ａは、夜間に撮像された画像とこの画像に対応するヒストグラムを示している。主要なピーク（significant peaks）は、輝度値が１２８以下のグレーレベルの領域内に存在する。この画像に対し、従来の全体的なヒストグラム均等化（histogram equalization）を適用すると、図１Ｂに示すような変更画像が得られる。この変更画像に対応するヒストグラムからわかるように、この処理により、主要なピークは、輝度の全範囲に亘って分布するようになる。ここで、図１Ｂに示す画像の明るさの印象は、十分に満足できるものではない。この画像は夜間に撮像されたものであるが、変更画像は、昼間に撮像された画像のように見える。
【００３８】
図２Ａは、他の画像を対応するヒストグラムとともに示している。この画像に従来の全体的なヒストグラム均等化を施すことにより、図２Ｂに示すような画像が得られる。図２Ｂに示す画像では、コントラストを強調しすぎたときに生じる典型的なアーティファクト（artefact）が見られる。すなわち、この女性の顔は、コントラストが誇張されているために不自然に見える。このような場合、ヒストグラム均等化の度合い（amount）を減らすことにより、より自然な結果が得られる。
【００３９】
本発明では、ヒストグラム均等化の度合いを、自動的に画像のコントラストに適応させる。画像のコントラストが高い場合、ヒストグラム均等化は不要であり、したがってヒストグラム均等化をオフにすることができる。一方、コントラストが低い場合、ヒストグラム均等化をオンにし、全体的なヒストグラム均等化を行う。また、画像のコントラストが中間的である場合、ヒストグラム均等化の効果が適正な量となるように調整する。
【００４０】
この方法では、まず、オリジナル画像のコントラストを分析する。コントラストが低い画像では、ヒストグラムの主要なピークは、低い輝度値の領域又は高い輝度値の領域に集中する傾向がある。これは、照明が不適切な条件下で画像が撮影されたため、画像のコントラストが低下し、画像が明るすぎる又は暗すぎることを示している。ヒストグラムピークが中間的な輝度値の範囲に分布している場合、この画像のコントラストは低くは見えない。また、ヒストグラムピークが存在しない空き領域がある場合（すなわち、その領域に対応する輝度値を有する画素がない場合）、画像のコントラストは高く見える。実際、画像のコントラストを高めるために、ヒストグラム均等化ではヒストグラムピークが存在しない空き領域が生じるために、ヒストグラム均等化は、このような視覚効果をサポートしている。
【００４１】
画像コントラストの特性測定値（characteristic measure）を算出する第１の方法では、画像コントラストの特性測定値Δを以下の式に基づいて算出する。
Δ＝ｋ_０×ｍａｘ（λ_０，λ_２）−ｋ_２×ｍｉｎ（λ_０，λ_２）−ｋ_１×λ_１（１）
ここで、λ_０は輝度レベル０から最初の主要ヒストグラムピークまでの距離を表し、λ_２は輝度レベル２５５から最後の主要ヒストグラムピークまでの距離を表し、λ_１は主要ヒストグラムピーク間でヒストグラムピークが存在しない領域の和を表し、ｋ_０、ｋ_１、ｋ_２は、３つの所定の係数を表す。これらの係数ｋ_０、ｋ_１、ｋ_２は、例えば１に設定してもよく、この場合、以下の式が成り立つ。
Δ＝ｍａｘ（λ_０，λ_２）−ｍｉｎ（λ_０，λ_２）−λ_１（２）
実際には、特性測定値Δには上限及び下限を設ける必要がある。
【００４２】
例えば、殆どの画素が低い輝度値の範囲に集中している場合、λ_２は、λ_０より大きくなり、ｍａｘ（λ_０，λ_２）は、λ_２の値となり、ｍｉｎ（λ_０，λ_２）は、λ_０の値となる。ｍｉｎ（λ_０，λ_２）を減算しているのは、殆どの画素が輝度の中間的な範囲に分布していることを考慮したものである。この結果、式（１）又は式（２）に基づいて算出される特性測定値Δの値は、照明が適切な条件の下ではより小さくなる。また、λ_１を減算しているのは、ヒストグラムピークが存在しない領域を考慮したものである。λ_１≠０である場合、特性測定値Δはより小さくなり、したがって、ヒストグラム均等化の度合いも減らされる。このλ_１の他に、ヒストグラム均等化の度合いは、ｋ_１に大きく依存する。ｋ_１を大きくする程、ヒストグラム均等化の度合いは小さくなる。
【００４３】
主要ヒストグラムピークの位置は、ＨＥ変換曲線の傾斜を分析することにより求めることができる。整数フォーマットＨＥ変換曲線（integer format HE transformation curve）の傾斜が２より小さい場合、この曲線は、オリジナルの輝度値の変更に貢献しない。したがって、傾斜が２以上である主要ピークを検出する。また、ヒストグラムから同様の基準（criterion）を得ることもできる。
【００４４】
次に、画像コントラストの特性測定値Δを算出する第２の方法について説明する。コントラストが低い画像では、画像の輝度値の平均値は、通常の画像の輝度値の平均値より極端に小さく又は大きくなる。この画像の輝度値の平均値と定数、例えば８ビット量子化画像では１２８と間の差を画像コントラストの特性測定値として用いることができ、更にこの特性測定値を用いてヒストグラム均等化の度合いを制御することができる。
【００４５】
ヒストグラムは、通常、オリジナル画像から生成される。ここで、１９９２年１１月ＣＶＧＩＰ：グラフィカルモデルズアンドイメージプロセッシング（Graphical Models and Image Processing）Ｖｏｌ．５４，Ｎｏ．６、第４９７〜５０６頁、ジェイ・ジー・リウ（J.-G.Leu）著「エッジ画素の強度に基づいた画像コントラストの強調（Image contrast enhancement based on the intensities of edge pixels）」では、雑音増幅の抑制とコントラストの向上の観点からは、周知のエッジ検出処理によってオリジナル画像からエッジ画像を検出し、エッジ画像からヒストグラムを作成する方が効果的であることが指摘されている。この手法では、エッジ位置周辺の画素のみを用いてヒストグラムを作成している。したがって、画像のエッジ程多くの情報を含まない、例えば背景等の均質な領域内の画素は、ヒストグラムに大きな影響を与えない。このヒストグラムのエントリは、画像の構造に対応している。したがって、この手法により、ヒストグラム均等化による情報の欠落を防止することができる。更に、画像の中央領域のエッジのみを用いてヒストグラムを求めることにより、ＰＡＬｐｌｕｓテレビジョン信号（PALplus TV signal）のブラックバー等の幾つかのコンテンツがヒストグラムに影響を与えることを防ぐことができる。
【００４６】
画像のコントラストは、オリジナル画像のヒストグラム又はＣＤＦから測定することができる。画像のコントラストを高い信頼度で測定するために、ヒストグラムを、オリジナル画像とエッジ画像の両方から算出してもよい。すなわち、オリジナル画像のヒストグラムと、エッジ画像のヒストグラムの両方に対して、特性測定値Δを算出する。ここで、２つの特性測定値Δの値が異なる場合、小さい方の値を用いてヒストグラム均等化の度合いを制御する。このように、２つの特性測定値Δの値のうちの小さい方の値を用いることにより、画像のコントラストを強調しすぎることを防止することができる。
【００４７】
画像の性質が過度に変化することを回避するためには、ヒストグラム均等化の度合いを抑制する必要がある。これは、周波数領域のヒストグラムを変更することにより実現できる。
【００４８】
本発明に基づき、ヒストグラム均等化の度合いを適応化する方法を図３を用いて説明する。まず、オリジナル画像１に対し、好ましくは高速フーリエ変換であるフーリエ変換２を施し、オリジナル画像１の周波数スペクトル３を求める。周波数スペクトル３は、ＤＣ成分４と、低周波数成分５と高周波数成分６とから構成される。ヒストグラム均等化の度合いを下げるには、高周波数成分６をローパスフィルタリングして減衰させる。図３には、このとき用いるローパスフィルタの周波数特性７も示されている。このローパスフィルタリングにより、カットオフ周波数ｆ_ｇ以上の全ての周波数成分が減衰され、カットオフ周波数ｆ_ｇ以下の周波数成分は変化しない。
【００４９】
ヒストグラム均等化の度合いは、矢印８に示すように、カットオフ周波数ｆ_ｇをシフトさせることにより変更することができる。高周波数成分を制限することにより、画像のコントラストを過度に強調してしまうことを避けることができる。すなわち、周波数がより高い成分をより多く減衰させることにより、ヒストグラム均等化の度合いがより少なくなる。ここで、ＤＣ成分のみを通過させるフィルタリングを行った場合、傾きが１の直線（unity straight line）が得られる。すなわち、カットオフ周波数ｆ_ｇの値を大きくすると、ヒストグラムの減衰量が少なくなり、ヒストグラム均等化の度合いが強くなる。一方、カットオフ周波数ｆ_ｇの値を小さくすると、ヒストグラムの減衰量が多くなり、ヒストグラム均等化の度合いが弱くなる。
【００５０】
矢印９で示すローパスフィルタリングの結果、周波数スペクトル１０が得られる。すなわち、このローパスフィルタリングにより、ＤＣ成分４と低周波数成分５のみが残される。そして、逆フーリエ変換１１を実行することにより、フィルタリングされたヒストグラム１２が得られる。フィルタリングされたヒストグラム１２は、オリジナルのヒストグラム１に比べて平滑化されている。このローパスフィルタリングにより、ヒストグラム均等化は、雑音による妨害を受けづらくなる。
【００５１】
これに代えて、上述したローパスフィルタリング処理を順方向及び逆方向ＩＩＲフィルタにより実行してもよく、すなわち、順方向でのフィルタリングを行った後、フィルリングされたシーケンスを逆にして、再びフィルタに戻してもよい。
【００５２】
図４Ａは、ある輝度レベルの分布、すなわちヒストグラムに対する累積分布関数（Cumulative Distribution Function：以下、ＣＤＦともいう。）を示している。この累積分布関数を換算（scaling）することにより、全体的なヒストグラム均等化を実行するためのＨＥ変換関数を得ることができる。
【００５３】
図４Ｂは、部分的なヒストグラム均等化を実行するための３つの異なるＨＥ変換曲線１３、１４、１５を示しており、これらは、図４Ａに示す累積分布関数と同じ輝度レベル分布に対応している。曲線１５は、カットオフ周波数が比較的低いローパスフィルタリングにより得られる強く減衰されたヒストグラムを積分することにより求められている。曲線１５は、傾きが１の線形の変換曲線（unity straight transformation curve）に非常に近く、したがって、画像に曲線１５を適用することにより得られるヒストグラム均等化の度合いは少ない。このＨＥ変換曲線１５は、コントラストが高い画像に最適である。
【００５４】
一方、ＨＥ変換曲線１３は、減衰量が多くないヒストグラムを積分することにより得られている。このような少ない減衰量は、カットオフ周波数が高いローパスフィルタリングに対応している。画像にＨＥ変換曲線１３を適用することにより、輝度レベルの分布は大きく変化する。すなわち、輝度値２００は、より明るい輝度値である輝度値約２４０にシフトする。これにより、コントラストが大きく強調されるため、ＨＥ変換曲線１３は、コントラストが低い画像に最適である。すなわち、曲線１３を用いることにより、ヒストグラム均等化の度合いが強くなる。
【００５５】
本発明をハードウェアとして実現したコントラスト強調装置の構成を図５に示す。このコントラスト強調装置は、入力画像１６に対し、ヒストグラム１７を算出する。ヒストグラム１７は、高速フーリエ変換（fast fourier transform：以下、ＦＦＴという。）ユニット１９に供給され、ＦＦＴユニット１９は、ヒストグラム１７の周波数スペクトルを算出して、この周波数スペクトルを乗算器２０に供給する。乗算器２０は、この周波数スペクトルをルックアップテーブル２１に記憶されているローパスフィルタ特性に乗算する。このようにしてローパスフィルタリングされた周波数スペクトルは、逆高速フーリエ変換（inverse fast fourier transform：以下、ＩＦＦＴという。）ユニット２２に供給される。ＩＦＦＴユニット２２は、ヒストグラムを平滑化し、平滑化されたヒストグラムを積分器２３に供給する。積分器２３は、平滑化されたヒストグラムを積分し、これにより平滑化された累積分布関数２４を生成する。累積分布関数２４は、換算ユニット２５において換算され、これによりＨＥ変換曲線２６が生成される。ＨＥ変換曲線２６は、再マッピングユニット２７において、入力画像１６の再マッピングに使用される。
【００５６】
ここで、カットオフ周波数ｆ_ｇを特性測定値Δに比例させることにより、又はフィルタの次数を特性測定値Δに反比例させることにより、ヒストグラム均等化の度合いを自動的に制御することができる。すなわち、カットオフ周波数ｆ_ｇを適応的に調整することにより、コントラストが異なる各画像に対してそれぞれ最適なヒストグラム均等化の度合いを得ることができる。
【００５７】
本発明の第２の具体例では、２つの異なるカットオフ周波数ｆ_ｇ１、ｆ_ｇ２を用いて２つのローパスフィルタリングを行う。これにより、２つの平滑化されたヒストグラムが得られる。最も強い周波数減衰応答は、低いカットオフ周波数及び高い画像コントラストに対応し、最も弱い周波数減衰応答は、高いカットオフ周波数及び低い画像コントラストに対応する。測定された画像コントラストが中間的な値である場合、ヒストグラム均等化の度合いは、最も強い周波数減衰応答と最も弱い周波数減衰応答とを補間することにより決定される。
【００５８】
本発明の第３の具体例では、オリジナルのヒストグラムを平滑化して得られたＨＥ変換曲線から傾きが１の線形の変換曲線（unity straight transformation curve）を減算する。次に、この減算により求められた差に０〜１の範囲をとるゲイン係数を乗算し、更にこの積を傾きが１の線形の変換曲線（unity straight transformation curve）に加算する。ここで、ゲイン係数を変更すると、ヒストグラム均等化の度合いが変化し、これにより、ヒストグラム均等化の度合いを測定された画像コントラストに適応化することができる。
【００５９】
上述したように、ヒストグラム均等化における問題は、ヒストグラム均等化の前後で画像の明るさが過剰に変化してしまうことである。以下では、この問題を解決する方法について説明する。
【００６０】
以下、明るさが過剰に変化してしまう問題を解決するための第１の方法について説明する。まず、ヒストグラム均等化の前及び後の画像の平均値を算出する。これらの平均値は、画像の画素の輝度値の総和を求めることにより直接算出してもよく、ヒストグラムから間接的に算出してもよい。平均値をヒストグラムから間接的に算出する方法は、演算負荷の軽減の観点から優れている。ここで、ヒストグラム均等化の前後の画像の輝度値の平均値の比を用いて、ヒストグラム均等化により得られた画像の画素の輝度値を換算する。ここで、ヒストグラム均等化により得られた画像がオリジナル画像より明るい場合、ヒストグラム均等化により得られた画像の全ての画素を小さく換算（down-scaled）し、逆にヒストグラム均等化により得られた画像がオリジナル画像より暗い場合、ヒストグラム均等化により得られた画像の全ての画素を大きく換算（up-scaled）する。これにより、ヒストグラム均等化の前後の画像の明るさを略同じにすることができる。
【００６１】
次に、明るさが過剰に変化してしまう問題を解決するための第２の方法について説明する。ここでは、ヒストグラム均等化の前後の画像の輝度値の平均値の比ではなく、これらの平均値間の差を求める。これにより、時間がかかる除算及び乗算処理が少なくなるため、演算負荷が軽減される。
【００６２】
ヒストグラム均等化により得られた画像の平均値がオリジナル画像の平均値より大きい場合、ヒストグラム均等化により得られた画像の全ての画素の輝度値は、この差に相当する量だけ小さくされる。逆に、ヒストグラム均等化により得られた画像の平均値がオリジナル画像の平均値より小さい場合、ヒストグラム均等化により得られた画像の全ての画素の輝度値は、この差に相当する量だけ大きくされる。これにより、ヒストグラム均等化による明るさの変化の問題を加算及び減算のみで解決することができる。本願発明者らによるシミュレーションの結果、この種の加算及び減算によっては、情報の欠落は観察されなかった。これは、以下のように説明される。まず、ヒストグラム均等化の後の画像の明るさが増した場合について検討する。この場合、オリジナル画像の画素の輝度値、特に暗い画素の輝度値が増大したことになる。したがって、輝度値が増大した画素からオフセット値を減算することにより、画像の性質が維持される。一方、ヒストグラム均等化により画像の明るさが低下した場合、これはオリジナル画像の画素の輝度値、特に明るい画素の輝度値が減少したことを意味する。したがって、この場合、画素の輝度値にオフセット値を加算することが望まれる。このオフセット値は、ヒストグラム均等化の前後の画像間の輝度値の差の平均値として求められる。ここで、経験的に、比較的明るい画像が良好な画像に見えるため、この輝度値の差の平均値に正の小数点の係数（positive fractional factor）を乗算することにより、画像の印象を向上させることができる。これは、特に暗い画像に適用すると効果的である。
【００６３】
次に、明るさが過剰に変化してしまう問題を解決するための第３の方法について説明する。画像の明るさの変化の問題は、ヒストグラムのピーク又は選択された基準点（reference point）を調べることによっても改善される。明るさの過剰な変化は、通常、ヒストグラムのピークにより生じる。ヒストグラムが比較的均等に分布していれば、ヒストグラム均等化により、明るさが過剰に変化することはない。したがって、ヒストグラム均等化の前後におけるヒストグラムピーク又は選択された基準点の位置シフトを検出し、この検出されたシフトによってヒストグラム均等化により得られた画像の画素の輝度を修正することができる。
【００６４】
上述した３つの方法は、いずれも全ての画素の輝度値を変更するものであり、したがって画像のコントラストは変化しない。
【００６５】
上述したように、ヒストグラム均等化は、例えば細部が闇に「隠された」画像等、コントラストが低い画像に対して効果を発揮する。ここで、画像には、明るい領域と暗い領域とがある場合がある。例えば大部分の領域が暗く、他の小さな領域が明るい画像があるとする。この画像のヒストグラムでは、主要な累積点（major accumulation）、すなわちピーク間に値が存在しない空間が観察される。包括的なヒストグラム均等化（global histogram equalization）によってもこの種の画像の暗い領域のコントラストを向上させることができるが、この場合、明るい領域の細部がある程度失われる可能性がある。すなわち、ここでは、暗い領域のコントラストを向上させることが目的であり、明るい領域の画素は、変更しない方がよい。
【００６６】
図６は、ＨＥ変換曲線２８の置換を説明する図である。ここでは、輝度レベル１２８〜２５５の領域２９におけるＨＥ変換曲線２８を対応する傾きが１の直線３０に置換している。これにより、ＨＥ変換曲線を平滑化した後に、この平滑化されたＨＥ変換曲線を用いて画像の輝度レベルを変換することができる。
【００６７】
これと同様の方法を大部分が明るい領域であり残りの僅かな領域が暗い画像に適用し、暗い領域において細部が失われることを防ぐこともできる。
【００６８】
【発明の効果】
以上のように、本発明に係るコントラスト強調方法は、それぞれに輝度レベルが割り当てられた複数の画素を含むデジタル画像のコントラストを強調するコントラスト強調方法において、オリジナル画像の輝度レベルの分布を示すヒストグラムを判定するステップと、オリジナル画像のコントラストの特性測定値を判定するステップと、オリジナル画像のヒストグラムを均等化し、コントラストを向上させるステップとを有し、ヒストグラムの均等化の度合いは、特性測定値の関数として変更される。これにより、アーティファクトの発生を回避し、変更画像をより自然な画像にすることができる。
【００６９】
また、本発明に係るコントラスト強調装置は、それぞれに輝度レベルが割り当てられた複数の画素を含むデジタル画像のコントラストを強調するコントラスト強調装置において、オリジナル画像の輝度レベルの分布を示すヒストグラムを判定するヒストグラム判定手段と、オリジナル画像のコントラストの特性測定値を判定するコントラスト判定手段と、オリジナル画像のヒストグラムを均等化し、コントラストを向上させるヒストグラム均等化手段とを備え、ヒストグラム均等化手段は、ヒストグラムの均等化の度合いを特性測定値の関数として変更する。これにより、アーティファクトの発生を回避し、変更画像をより自然な画像にすることができる。
【図面の簡単な説明】
【図１】図１Ａは、夜間に撮像された画像と、ヒストグラム均等化前の対応するヒストグラムを示す図であり、図１Ｂは、図１Ａに示す画像に対し、従来の全体的なヒストグラム均等化を行った後の画像及びヒストグラムを示す図である。
【図２】図２Ａは、コントラストが高い人間の顔の画像と、ヒストグラム均等化前の対応するヒストグラムを示す図であり、図２Ｂは、図２Ａに示す画像に対し、従来の全体的なヒストグラム均等化を行った後の画像及びヒストグラムを示す図である。
【図３】ヒストグラム均等化前の度合いを変更するためにオリジナルヒストグラムに対して行われるローパスフィルタリングを説明する図である。
【図４】図４Ａは、オリジナルヒストグラムを積分することにより得られる累積分布関数の具体例を示す図であり、図４Ｂは、ヒストグラムを部分的に均等化する平滑化されたヒストグラム均等化変換曲線を示す図である。
【図５】ローパスフィルタを備える本発明に基づくヒストグラム均等化装置の構成例を示すブロック図である。
【図６】所定の輝度レベル範囲内のヒストグラム均等化変換曲線を傾きが１の線形の変換曲線に置換する手法を説明する図である。
【符号の説明】
１オリジナル画像、２フーリエ変換、３周波数スペクトル、４ＤＣ成分、５低周波数成分、６高周波数成分、７周波数特性、８カットオフ周波数シフト、９ローパスフィルタリング、１０周波数スペクトル、１１逆フーリエ変換、１２フィルタリングされたヒストグラム

Claims

それぞれに輝度レベルが割り当てられた複数の画素を含むデジタル画像のコントラストを強調するコントラスト強調方法において、
（ａ）オリジナル画像の輝度レベルの分布を示すヒストグラムを判定するステップと、
（ｂ）上記オリジナル画像のコントラストを示す特性測定値を上記ヒストグラムから判定するステップと、
（ｃ）上記オリジナル画像のヒストグラムを均等化し、コントラストを向上させるステップと
を有し、
上記ヒストグラムの均等化の度合いは、上記特性測定値の関数として変更され、
上記（ｂ）ステップは、オリジナル画像の輝度レベルの平均値と所定の定数との絶対差を評価することにより実行されて、該差が大きいほど画像のコントラストが低いと判定し、
上記（ｃ）ステップは、変更輝度レベルを算出するために、ヒストグラム均等化変換曲線に基づいて、オリジナル画像の輝度レベルを変換することにより実行され、
上記ヒストグラム均等化変換曲線は、ヒストグラムを積分及び換算することにより生成されることを特徴とするコントラスト強調方法。
それぞれに輝度レベルが割り当てられた複数の画素を含むデジタル画像のコントラストを強調するコントラスト強調方法において、
（ａ）オリジナル画像の輝度レベルの分布を示すヒストグラムを判定するステップと、
（ｂ）上記オリジナル画像のコントラストを示す特性測定値を上記ヒストグラムから判定するステップと、
（ｃ）上記オリジナル画像のヒストグラムを均等化し、コントラストを向上させるステップと
を有し、
上記ヒストグラムの均等化の度合いは、上記特性測定値の関数として変更され、
上記（ｂ）ステップは、輝度レベル０から最初の主要ヒストグラムピークまでの重み付けされた距離と、最後の主要ヒストグラムピークから最大の輝度レベルまでの重み付けされた距離との絶対差を評価することにより実行されて、該重み付けされた距離の差が大きい程画像のコントラストが低いと判定し、
上記（ｃ）ステップは、変更輝度レベルを算出するために、ヒストグラム均等化変換曲線に基づいて、オリジナル画像の輝度レベルを変換することにより実行され、
上記ヒストグラム均等化変換曲線は、ヒストグラムを積分及び換算することにより生成されることを特徴とするコントラスト強調方法。
上記（ｂ）ステップの実行において、主要ヒストグラムピーク間の空き領域の総和が評価されて、該主要ヒストグラムピーク間の空き領域の総和が大きい程画像のコントラストが高いと判定することを特徴とする請求項２記載のコントラスト強調方法。
上記主要ヒストグラムピークは、ヒストグラム又はヒストグラム均等化変換曲線の傾きが所定の閾値をどこで超えているかを判断することにより判定されることを特徴とする請求項２又は３記載のコントラスト強調方法。
上記（ｂ）ステップの実行において、エッジ画像がオリジナル画像に対してエッジ検出処理を施すことにより判定され、上記オリジナル画像の上記特性測定値と、上記エッジ画像の上記特性測定値の両方が算出され、上記算出された特性測定値のうち小さいほうの値が上記ステップ（ｃ）において使用されることを特徴とする請求項１乃至４いずれか１項記載のコントラスト強調方法。
上記ヒストグラムに対しローパスフィルタリング処理を施した後に上記ヒストグラム均等化変換曲線を生成することを特徴とする請求項１乃至５いずれか１項記載のコントラスト強調方法。
上記（ｃ）ステップでは、上記ヒストグラム均等化の度合いは更に、上記ローパスフィルタリング処理のカットオフ周波数を変化させることにより変更され、該カットオフ周波数が高い程、ヒストグラム均等化の度合いが強くなることを特徴とする請求項６記載のコントラスト強調方法。
上記特性測定値が高い画像コントラストを示している場合、カットオフ周波数はより低い周波数にシフトされ、上記ヒストグラムは上記ヒストグラム均等化変換曲線を生成する前に大きく減衰され、上記特性測定値が低い画像コントラストを示している場合、上記カットオフ周波数はより高い周波数にシフトされ、上記ヒストグラムは上記ヒストグラム均等化変換曲線を生成する前に小さく減衰されることを特徴とする請求項６又は７記載のコントラスト強調方法。
上記ローパスフィルタリング処理を実行するために、高速フーリエ変換を含むフーリエ変換を実行して上記ローパスフィルタリング処理を周波数領域で行い、又は順方向及び逆方向ＩＩＲフィルタリング処理を行うことを特徴とする請求項６乃至８いずれか１項記載のコントラスト強調方法。
上記（ｃ）ステップでは、上記ヒストグラム均等化の度合いは更に、上記ヒストグラム均等化変換曲線を生成する際、第１のカットオフ周波数に対応する第１のローパスフィルタリングされたヒストグラムと、第２のカットオフ周波数に対応する第２のローパスフィルタリングされたヒストグラムとの間を補間することにより変更されることを特徴とする請求項１乃至５いずれか１項記載のコントラスト強調方法。
上記（ｃ）ステップでは、上記ヒストグラム均等化の度合いは更に、ヒストグラム均等化変換曲線と、傾きが１の線形の変換曲線との差をゲイン係数により換算することにより変更され、該ゲイン係数が小さいほど上記ヒストグラム均等化の度合いが弱くなることを特徴とする請求項１乃至５いずれか１項記載のコントラスト強調方法。
ヒストグラム均等化の前後の画像の平均輝度レベルの比に基づいて、ヒストグラム均等化後の画像の輝度レベルを換算することを特徴とする請求項１乃至１１いずれか１項記載のコントラスト強調方法。
ヒストグラム均等化の前後の画像の平均輝度レベルの差に基づいて、ヒストグラム均等化後の画像の輝度レベルをシフトすることを特徴とする請求項１乃至１１いずれか１項記載のコントラスト強調方法。
ヒストグラム均等化により生じたヒストグラムピーク又は選択された基準点の輝度レベルのシフトに基づいて、ヒストグラム均等化後の画像の輝度レベルをシフトすることを特徴とする請求項１乃至１１いずれか１項記載のコントラスト強調方法。
輝度レベルの任意の範囲について、上記均等化変換曲線を傾きが１の線形の変換曲線に置換することにより該均等化変換曲線を変更することを特徴とする請求項１記載のコントラスト強調方法。
それぞれに輝度レベルが割り当てられた複数の画素を含むデジタル画像のコントラストを強調するコントラスト強調装置において、
オリジナル画像の輝度レベルの分布を示すヒストグラムを判定するヒストグラム判定手段と、
上記オリジナル画像のコントラストを示す上記ヒストグラムから特性測定値を判定するコントラスト判定手段と、
上記オリジナル画像のヒストグラムを均等化し、コントラストを向上させるヒストグラム均等化手段と
を備え、
上記ヒストグラム均等化手段は、上記ヒストグラムの均等化の度合いを上記特性測定値の関数として変更し、
上記ヒストグラム均等化手段は、変更輝度レベルを算出するために、ヒストグラム均等化変換曲線に基づいて、オリジナル画像の輝度レベルを変換することにより、ヒストグラムを均等化し、
上記コントラスト判定手段は、オリジナル画像に対してエッジ検出処理を施すことによりエッジ画像を生成する手段を備え、上記オリジナル画像の上記特性測定値と、上記エッジ画像の上記特性測定値の両方が算出され、上記両方の算出された特性測定値のうち値が小さい方が前記コントラスト判定手段によって前記特性測定値として考慮され、
上記ヒストグラム均等化手段は、上記ヒストグラム均等化曲線を生成する前に上記ヒストグラムをローパスフィルタリングするローパスフィルタリング手段を備えることを特徴とするコントラスト強調装置。
前記ヒストグラム均等手段は更に、上記特性測定値に応じて変化する上記ローパスフィルタリング処理のカットオフ周波数を変化させることにより上記ヒストグラム均等化の度合いを変更し、該カットオフ周波数が高い程、ヒストグラム均等化の度合いが強くなることを特徴とする請求項１６記載のコントラスト強調装置。
上記ヒストグラム均等化手段は、上記ヒストグラムを周波数領域に変換するフーリエ変換手段を備え、上記ローパスフィルタリング処理を周波数領域で行い、又は上記ヒストグラム均等化手段は、ＩＩＲフィルタを備え、順方向及び逆方向ＩＩＲフィルタリング処理により上記ローパスフィルタリング処理を行うことを特徴とする請求項１６又は１７いずれか１項記載のコントラスト強調装置。
上記ローパスフィルタリング手段は、フーリエ変換されたローパスフィルタ特性が登録されているルックアップテーブルを備えることを特徴とする請求項１６乃至１８いずれか１項記載のコントラスト強調装置。
前記ヒストグラム均等手段は更に、上記ヒストグラム均等化変換曲線を生成する際、第１のカットオフ周波数に対応する第１のローパスフィルタリングされたヒストグラムと、第２のカットオフ周波数に対応する第２のローパスフィルタリングされたヒストグラムとの間を補間することにより上記ヒストグラム均等化の度合いを変更することを特徴とする請求項１６記載のコントラスト強調装置。
前記ヒストグラム均等手段は更に、ヒストグラム均等化変換曲線と、傾きが１の線形の変換曲線との差をゲイン係数により換算することにより上記ヒストグラム均等化の度合いを変更し、該ゲイン係数が小さいほど上記ヒストグラム均等化の度合いが弱くなることを特徴とする請求項１６記載のコントラスト強調装置。
ヒストグラム均等化後の画像の平均輝度レベルがオリジナル画像の平均輝度レベルに等しくなるように、ヒストグラム均等化により生じた明度の変化を補償する明度補償手段を備える請求項１６乃至２１いずれか１項記載のコントラスト強調装置。
輝度レベルの任意の範囲について上記輝度レベルが変化しないように該均等化変換曲線を変更する均等化変換曲線変更手段を備える請求項１６乃至２２いずれか１項記載のコントラスト強調装置。
コンピュータ又はデジタルシグナルプロセッサにより実行されて請求項１乃至１５いずれか１項記載のコントラスト強調方法を実行するコンピュータプログラムを備えるコンピュータプログラム。