JP5624896B2 - 画像処理装置、画像処理プログラム、及び、画像処理方法 - Google Patents

Description

本発明は、画像処理装置、画像処理プログラム、及び、画像処理方法に関する。

通常のイメージセンサで取得した画像の輝度範囲、即ちダイナミックレンジは、人の知覚できるダイナミックレンジと比べると狭い。このため、近年では異なる輝度範囲を撮影した複数の画像を合成することにより広いダイナミックレンジの画像(HDR（high dynamic range）画像)を生成する技術が使われている。その結果、より見たままに近いシーンが記録できるようになった。しかし、表示デバイスである各種のディスプレイの多くは、依然低いダイナミックレンジを有する。低いダイナミックレンジの表示デバイスにHDR画像を表示する際には、不自然にならないように階調を圧縮する必要がある。

一方、人の視覚は広い輝度範囲を認識できるが、実際には見ている領域の明暗によってゲインを変えて知覚しており、知覚のダイナミックレンジはそれほど広くない。このような視覚特性をモデル化して、階調圧縮に応用する方法が提案されているが、その方法の一つとして、良好な階調圧縮ができる「iCAM06」が知られている（非特許文献１、２参照）。

また、「iCAM06」のような視覚モデルに基づいた階調補正は、低いダイナミックレンジで記録したLDR（low dynamic range）画像に対する階調補正処理に応用されても、良好な階調補正の結果を生じる。なお、LDRに対する階調補正には、例えば、暗部や明部を適正な明るさにする階調変換がある。

なお、以下、階調補正は、画像の暗部や明部を適正な明るさにしたりする階調変換の他、画像のダイナミックレンジを狭める階調圧縮などの階調変換も含む。

Kuang, J., Johnson, G. M., Fairchild M. D.,「iCAM06: A refined image appearance model for HDR image rendering」, Journal of Visual Communication and Image Representation, Volume 18, Issue 5, October 2007, Pages 406-414 Mark D. Fairchild著、"Color Appearance Models"、（英国）、第２版、John Wiley & Sons社、２００５年１月１日

しかし、視覚モデルを応用した階調補正は、計算が複雑であり、多くの演算量を必要とする。このため、小さな回路規模の画像処理装置やソフトウェアが階調補正処理を行う場合に、静止画の階調補正処理に要する時間は多くなり、また、リアルタイムでの動画像の階調補正処理は困難である。

例えば、iCAM06において、国際標準のカラーアピアランスモデルであるCIE CAM02の錐体応答モデルを数式（１）及び（２）のように応用して、階調補正が行われている。

ここで、tcは、錐体応答の出力（階調補正値）で、tcを求める関数は錐体応答のモデルである。また、baseは、エッジ保存平滑化による平滑化画像（ベース成分画像）のリニア画素値であり、whiteは、平滑化した画像（平滑化画像）のリニア画素値で、基準の順応白色の輝度を示す。FLは、局所的な順応度合いを決定する中間計算値で、whiteから算出される。pは、コントラストを調整するユーザパラメータで標準値は0.75である。

実際には、エッジ保存平滑化の出力値及び錐体応答の出力値は対数であることが望ましい。このため、baseの対数値をbase’、錐体応答の出力の対数値をtc’とすると、錐体応答の演算式は、数式（３）となる。

このように錐体応答に基づく階調補正処理には、通常、多くの乗除算、冪演算、対数演算が含まれ、多くの演算量を必要とする。

本発明の目的は、錐体応答に基づく階調補正処理を高速に且つ／又は小さい演算回路で実施できるようにすることである。

本発明のある態様に係る画像処理装置は、錐体応答に基づく階調補正を行う画像処理装置であって、階調補正後の画像の画素値の対数tc'を、数式tc' ＝ a・base'＋b・white'＋cを満たすよう演算する階調補正手段を備える。ここで、base'は、階調補正前の画像に対しエッジ保存平滑化及び対数変換を行うことにより得られる第一の画像の画素値を表し、white'は、階調補正前の画像に対し平滑化及び対数変換を行なうことにより得られる第二の画像の画素値を表し、a、b及びcは、それぞれ、第一の定数、第二の定数、及び、第三の定数を表す。

本発明の別の態様に係る画像処理方法は、錐体応答に基づく階調補正を行う画像処理方法であって、階調補正後の画像の画素値の対数tc'を、数式tc' ＝ a・base'＋b・white'＋cを満たすよう演算する階調補正ステップを備える。

本発明のさらに別の態様に係る画像処理プログラムは、錐体応答に基づく階調補正を行う画像処理プログラムであって、階調補正後の画像の画素値の対数tc'を、数式tc' ＝ a・base'＋b・white'＋cを満たすよう演算する階調補正手順をコンピュータに実行させる。

本発明によると、錐体応答の演算が簡易に行える。従って、錐体応答に基づく階調補正処理が、高速に或いは/且つ小さい演算回路で実施できる。

第一実施形態に係る画像処理装置を示すブロック図である。 第一実施形態に係る補正処理部の詳細を示すブロック図である。 錐体応答の出力値tc’とAの関係を示すグラフである。 基準の順応白色の対数輝度値white’とlog(FL/white)との関係を示すグラフである。 第一実施形態に係る錐体応答演算部の機能ブロック図である。 第二実施形態に係る補正処理部の詳細を示すブロック図である。 （ａ）画像処理プログラムのメインルーチンを示すフローチャートである。（ｂ）補正処理のサブルーチンを示すフローチャートである。（ｃ）錐体応答演算処理のサブルーチンを示すフローチャートである。

［第一実施形態］
<<画像処理装置>>
図１は、第一実施形態に係る画像処理装置１（又は画像処理システム）のブロック図である。画像処理装置１は電子機器に搭載される。電子機器の例として、テレビジョン受像機、プリンタ、スチルカメラ、ビデオカメラなどが挙げられる。画像処理装置１は、入力した画像に対して、簡略化した錐体応答の演算に基づく階調補正を行う。

画像処理装置１は、入力されたカラー画像（例えばＲＧＢ画像）に対して、階調補正処理を行う。画像処理装置１は、輝度画像生成部１０、補正処理部１１、ゲイン生成部１２、ゲイン適用部１３を備える。例えば、各部１０−１３は論理回路から構成してよいし、ＣＰＵ（中央演算処理装置）と演算プログラムを格納するメモリなどから構成してもよい。

輝度画像生成部１０は、カラー画像を輝度画像に変換するもので、例えば、画素ごとにＲＧＢ−ＸＹＺ変換などにより輝度値（Ｙ）を得て、輝度値（Ｙ）を画素値として有する輝度画像Ｉを生成する。補正処理部１１は、輝度画像に対して階調補正処理を行う。ゲイン生成部１２は、カラー画像の各色の階調変換などに使用するゲインを生成する。例えば、ゲイン生成部１２は、補正処理部１１による処理後の画像と処理前の画像との間の各画素における輝度値の比を、カラー画像の各色のゲインに設定する。ゲイン適用部１３は、ゲイン生成部１２で生成されたゲインを、カラー画像の各色の階調変換などに適用し、ゲイン適用部１３は、入力したカラー画像の画素値にゲインを印加する。

なお、画像処理装置１が入力された単色（モノクロ）画像に対して、階調補正処理を行う場合、輝度生成部１０、ゲイン生成部１２、及び、ゲイン適用部１３を省略して、補正処理部１１の出力をそのまま画像処理装置１の出力としてよい。

<<補正処理部>>
図２は、補正処理部１１の詳細を示すブロック図である。補正処理部１１は、対数変換部２０、エッジ保存平滑化部２１、平滑化部２２、ディテール成分生成部２３、錐体応答演算部２４、合成処理部２５、リニア変換部２６を備える。例えば、各部２０−２６は論理回路から構成してよいし、ＣＰＵ（中央演算処理装置）と演算プログラムを格納するメモリなどから構成してもよい。

対数変換部２０は、輝度画像生成部１０から入力された輝度画像Ｉを対数変換して対数の輝度画像Ｉ’を生成する。エッジ保存平滑化部２１は、対数変換部２０から入力された対数輝度画像Ｉ’に対してエッジ保存平滑化を実施して、エッジが保存された平滑化画像（対数のベース成分画像Ｂ’）を生成する。例えば、エッジ保存平滑化には、公知のバイラテラルフィルタが使用される。

平滑化部２２は、対数変換部２０から入力された対数輝度画像Ｉ’に対して平滑化を行い、基準の順応白色の対数輝度画像として対数の平滑化画像Ｗ’を生成する。平滑化部２２の平滑化は、エッジを保存しないものでよく、エッジ保存平滑化部２１のエッジ保存平滑化とは異なる度合いで行われる。例えば、非特許文献１のようにガウシアンフィルタによる平滑化が行われる。

ディテール成分生成部２３は、対数変換部２０から入力された対数輝度画像Ｉ’とエッジ保存平滑化部２１から入力された対数ベース成分画像Ｂ’の差を計算して、細部の成分に関するディテール成分画像Ｄ’（対数）を得る。錐体応答演算部２４は、エッジ保存平滑化部２１から入力された対数ベース成分画像Ｂ’と平滑化部２２からの対数平滑化画像Ｗ’（基準の白色対数画像）を使用して、錐体応答モデルに基づいて対数画素値の階調の変換を行う。錐体応答演算部２４は、階調処理手段又は階調処理部とも呼ばれる。

合成処理部２５は、ディテール成分生成部２３からのディテール成分画像Ｄ’と錐体応答演算部２４からの対数階調補正画像Ｓ’との画素値の和を演算して、これら画像を合成する。リニア変換部２６は、合成処理部２５で生成された合成画像Ｔ’をリニア変換（対数変換の逆変換）して出力画像Ｒとして出力する。

<<錐体応答>>
次に、錐体応答演算部２４が実施する錐体応答の演算について説明する。本実施形態では、錐体応答の演算式である前述の数式（３）の計算が簡易に行えるような工夫が施される。

演算の簡易化を考えるために、まず数式（３）は、数式（４）と数式（５）のように分離される。各式のlog関数は常用対数であるものとする。

図３のグラフの実線は、式（５）のtc’とAの関数を示す。しかし、実際には、tc’とAの関数のうち直線部が使われる場合がほとんどである。従って、図３のグラフの点線のように関数を直線近似して、直線部のみを使用しても実用上の問題は生じない。

即ち、数式（５）の代わりに、点線を表す数式（６）の直線近似式、即ち一次関数を使用してもよい。

一方、数式（４）は、次のように展開できる。

ここで、log(FL/white)は、whiteのみから計算される値であり、数式（２）を用いて、数式（８）のように展開できる。

さらに、whiteの対数値であるwhite’とlog(FL/white)との関係は、図４のグラフの実線のようにほぼ直線（一次関数）となる。white’がごく小さい場合には、white’とlog(FL/white)との関係は、曲線となる。しかし、図４のグラフ中に点線のように関数を直線近似して、直線部のみを使用しても実用上の問題は生じない。

即ち、数式（８）の代わりに、点線に対応する数式（９）の直線近似式、即ち一次関数を使用してもよい。

従って、数式（７）に示したAは、式（９）を代入して、数式（１０）のように近似的に計算できる。

更にAの近似式（１０）を錐体応答値tc’の近似計算である式（６）に代入すると式（１１）が得られる。

ここで、a,b,cは、ユーザが決定してメモリに保存されるパラメータpによって計算される定数であるので、画素ごとの計算は不要である。

このように、錐体応答の近似を行うことで、視覚モデルに基づいた多くの演算を必要とする処理である数式（３）が、わずかな演算量で計算が行える数式（１１）に近似できる。そして、錐体応答演算部２４は、数式（１１）の演算を行う。

<<錐体応答演算部>>
図５は、錐体応答演算部（階調処理手段）２４の詳細の一例を示すブロック図である。錐体応答演算部２４は、第一乗算部３０、第二乗算部３１、加算部３２を備える。例えば、各部３０−３２は論理回路から構成してよいし、ＣＰＵ（中央演算処理装置）と、演算プログラムを格納するメモリなどから構成してもよい。

第一乗算部３０は、第一の入力画像としての対数ベース成分画像Ｂ’に第一の定数aを乗算する乗算処理を行う。即ち、第一乗算部３０は、対数ベース成分画像Ｂ’の各画素の画素値base'に第一の定数aを乗算して第一の積a・base'を算出する。第二乗算部３１は、第二の入力画像としての対数平滑化画像Ｗ’に第二の定数ｂを乗算する乗算処理を行う、即ち、第二乗算部３１は、対数平滑化画像Ｗ’の各画素の画素値white'に第一の定数ｂを乗算して第二の積ｂ・white'を算出する。

加算部３２は、第一乗算部３０の演算結果aＢと第二乗算部３１の演算結果ｂＷと第三の定数cとを足し合わせる加算処理を行い、対数階調補正画像Ｓ’＝aＢ’＋ｂＷ’＋ｃを出力する。即ち、加算部３２は、数式（１１）に基づいて、階調補正画像の画素値tc'（＝a・base'＋ｂ・white'＋ｃ）を求める。このように加算部３２は、base’を定数倍（ａ倍）したものとbase’を定数倍（ｂ倍）したものを足し合わせ、base’とwhite’の線形和を計算する加算手段となる。

なお、数式（１１）をtc'＝a・｛base'＋（ｂ／ａ）・white' ｝＋ｃと変形して計算するような場合にも、錐体応答演算部２４は、base’を定数倍(１倍）したものとbase’を定数倍（ｂ／ａ倍）したものを足し合わせ、線形和｛base'＋（ｂ／ａ）・white' ｝を計算する加算手段を有することになる。

数式（１１）において、定数ｃは、対数値である各画素値に対して一様に加算（ｃがマイナスであれば減算）をする値である。対数値を逆演算してリニアに変換した画像に関して、定数ｃの加算は、一様に画素値を変倍（スケーリング）することに相当する。即ち、定数ｃの加算は、階調の補正には寄与しない。このため、加算部３２は、定数ｃの加算を行わなくてもよい。この場合でも、階調補正の効果は得られる。

−作用効果−
階調補正する画像処理装置は、階調補正後の画像の画素値の対数tc'を、数式tc' ＝ a・base'＋b・white'＋c（数式（１１））を満たすよう演算する階調処理手段（錐体応答演算部２４）を備える。ここで、base'は、階調処理前の画像に対しエッジ保存平滑化及び対数変換を行うことにより得られる第一の画像(Ｂ’)の画素値を表す。white'は、階調処理前の画像に対し平滑化及び対数変換を行なうことにより得られる第二の画像(Ｗ’)の画素値を表す。

これにより、錐体応答の演算式の計算が簡易に行える。従って、視覚モデルに基づいた階調補正処理が、高速に或いは/且つ小さい演算回路で実施できる。典型的には、数式（１１）を用いる本実施形態の階調補正処理は、数式（３）を用いる階調補正処理より、処理時間は７０分の１程度になる。また、数式（１１）による階調補正処理を行う本実施形態の画像処理装置は、数式（３）による階調補正処理を行う画像処理装置より、回路規模は５分の１程度になる。

階調処理手段（錐体応答演算部２４）は、第一の画像の各画素の画素値を定数倍した値を取得する第一の取得手段（例えば、第一乗算部３０）と、第二の画像の各画素の画素値を定数倍した値を取得する第二の取得手段（例えば、第二乗算部３１）と、を備える。さらに、階調処理手段は、第一の画像と第二の画像とで対応する画素同士について第一の積と第二の積を足し合わせる加算手段（例えば、加算部３２）を備える。このため、好適に数式（１１）が計算できる。

なお、第三の定数ｃがゼロであれば、数式（１１）の演算が簡略化でき、錐体応答の演算式の計算がさらに高速に行える。

［第二実施形態］
第一実施形態では、輝度画像の階調補正の結果を得るために、ディテール成分画像Ｄ’を生成する。しかし、カラー画像を処理したい場合はゲインのみが取得できれば良いので、ディテール成分画像Ｄ’を生成する必要はない。

図６は、第二実施形態に係る補正処理部４１とゲイン生成部４２の詳細を示すブロック図である。補正処理部４１では、第一実施形態に係る補正処理部１１からディテール成分生成部２３、合成処理部２５、リニア変換部２６が省略されている。ゲイン生成部４２は、対数値ゲイン生成部４４とリニア変換部４５を備える。

ゲイン生成部４２の対数値ゲイン生成部４４は、錐体応答演算部２４による階調補正処理前後のゲインの対数値log(base/tc)を演算する。対数値log(base/tc)は、base'−tc'であるため、ゲインを除算でなく減算で計算できる。リニア変換部４５は、対数値log(base/tc)を逆対数変換してゲインbase/tcを求める。ゲイン適用部１３は、ゲイン生成部４２のリニア変換部４５で生成されたゲインを、カラー画像の各色の階調変換などに適用する。

第二実施形態において、階調補正処理前後のゲインを除算でなく減算で計算できる。このため、ゲインの演算量が減少する。

［第三実施形態］
第三実施形態は、錐体応答による階調変換の数式（１１）において定数ａ，ｂ，ｃがとりうる値の範囲に関する。

数式（１１）において、コントラスト調整の役割を持つユーザパラメータｐは、０から２程度の範囲で調整可能である。また、数式（１１）の定数ｂの計算に現れる値２／３は、順応度合い、即ち階調の圧縮の度合いを調整する値として機能する。このため、値２／３の代わりに順応度合いを調整するユーザパラメータｄを使用しても良い。ユーザパラメータｄは、極端な順応を行わせる値０から、順応させない値２程度の範囲で調整可能である。従って、ユーザパラメータｐとｄの積から決定される定数ｂ（＝−ｐ・ｄ）は、マイナス４から０程度までの値をとる。なお、ユーザパラメータｐとｄは、ユーザにより設定され、メモリに保存されている。

また、数式（１１）は、近似計算であり、定数値ｃの計算で用いている値１．１５８などは厳密に用いずに、１から１．５程度の幅があっても問題ない。即ち、数式（１１）での定数ａ，ｂ，ｃがとりうる値の範囲は、数式（１２）の通りに設定される。

第三実施形態によると、エッジ保存平滑化した画像が、階調補正を行いたい画像を常用対数で対数変換して生成した画像である。 第一の定数aは、０から２までの値であり、第二の定数ｂは、第一の定数に対してマイナス２倍から０倍した値であり、第三の定数ｃは、第一の定数をマイナス１倍したものに１から１．５までの値を加算した値である。従って、コントラストや階調の圧縮の度合いなどを適切に調整できる。

［第四実施形態］
第四実施形態において、第三の定数ｃは、ベース成分画像（リニア）の画素値baseと階調補正画像（リニア）の画素値tcの比が小さくなるように、階調補正画像の画素値tcをスケーリングするための調整値として用いられる。第三の定数ｃは、ユーザにより設定され、メモリに記憶しておくことができる。

カラー画像に対して階調補正する場合に、輝度成分の階調補正処理の前後のゲイン値が、カラー画像の各色の階調変換に適用される。この場合、ゲイン値が小さいほど、ビット数が少ない乗算を行え、処理時間の短縮、あるいはハードウェア規模の削減ができる。なお、ベース成分画像の画素値baseと階調補正画像の画素値tcの比が小さいほど、ゲイン値は小さい。

例えば、base及びwhiteが１から１００００の値をとる場合に、その対数値であるbase'及びwhite'は０から４となる。さらに、pを標準値である０．７５（ｃ＝０．４０８）として数式（１１）を適用すると、base'=４且つwhite'=４のとき、baseとtcの比（base/tc）が、最大値３９１．７４となる（数式（１３）)。

しかし、第三の定数ｃの値を増加して調整することにより、錐体応答の対数値tc’を対数ベース成分画像の画素値base'の範囲に近づけることができる。例えば、数式（１１）において、第三の定数ｃを（-ｐ+１．１５８）の代わりに１．５とすると、base'=４且つwhite'=４のとき、baseとtcの比（base/tc）は最大値３１．６３倍となり(数式（１４）)、よりビット数が少ないゲイン値を用いることができる。

第四実施形態によると、第三の定数ｃは、錐体応答の対数値tc’を第一の入力画像の画素値base'の範囲に近づける調整値である。従って、第三の定数ｃの調整により、ゲイン値を小さくして、ゲイン値を使用する演算の際にビット数が少ない高速の演算が行える。従って、階調補正処理の処理時間の短縮、あるいはハードウェア規模の削減ができる。

［その他の実施形態］
上述した各実施形態の説明では、画像処理装置が行う処理としてハードウェアによる処理を前提としていたが、このような構成に限定される必要はない。例えば、別途ソフトウェアにて処理する構成も可能である。この場合、画像処理装置は、コンピュータに相当し、ＣＰＵ、ＲＡＭ等の主記憶装置と、上記処理の全て或いは一部を実現させるためのプログラムが記憶されたコンピュータ読み取り可能な記憶媒体とを備えている。ここでは、このプログラムを画像処理プログラムと呼ぶ。そして、ＣＰＵが上記記憶媒体に記憶されている画像処理プログラムを読み出して、情報の加工・演算処理を実行することにより、上述の画像処理部と同様の処理を実現させる。

ここで、コンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、半導体メモリ等をいう。また、この画像処理プログラムを通信回線によってコンピュータに配信し、この配信を受けたコンピュータが当該画像処理プログラムを実行するようにしても良い。

図７（ａ）−（ｃ）は、画像処理装置（コンピュータ）が実行するプログラムを例示する。図７（ａ）は、画像処理のメインルーチンを示す。図７（ｂ）は、補正処理のサブルーチンを示す。図７（ｃ）は、エッジ保存平滑化処理のサブルーチンを示す。図７（ａ）のステップＳ１０―Ｓ１３で行われる処理は、それぞれ、図１の各部１０−１３が行う処理に対応する。図７（ｂ）のステップＳ２０―Ｓ２６で行われる処理は、それぞれ、図２の各部２０−２６が行う処理に対応する。図７（ｃ）のステップＳ３０―Ｓ３２で行われる処理は、それぞれ、図５の各部３０−３２が行う処理に対応する。

本発明は上記の各実施形態に限定されずに、その技術的な思想の範囲内において種々の変更がなしうることは明白である。

１０ 輝度画像生成部（輝度画像生成手段）
１１ 補正処理部（補正処理手段）
１２ ゲイン生成部（ゲイン生成手段）
１３ ゲイン適用部（ゲイン適用手段）
２０ 対数変換部（対数変換手段）
２１ エッジ保存平滑化部（エッジ保存平滑化手段）
２２ 平滑化部（平滑化手段）
２３ ディテール成分生成部（ディテール成分生成手段）
２４ 錐体応答演算部（錐体応答演算手段）
２５ 合成処理部（合成処理手段）
２６ リニア変換部（リニア変換手段）
３０ 第一乗算部（第一乗算手段）
３１ 第二乗算部（第二乗算手段）
３２ 加算部（加算手段）

Claims (8)

1. 錐体応答に基づく階調補正を行う画像処理装置であって、
   階調補正後の画像の画素値の対数tc'を、数式
   tc' ＝ a・base'＋b・white'＋c
   を満たすよう演算する階調補正手段を備え、
   ここで、base'は、階調補正前の画像に対しエッジ保存平滑化及び対数変換を行うことにより得られる第一の画像の画素値を表し、
   white'は、階調補正前の画像に対し平滑化及び対数変換を行なうことにより得られる第二の画像の画素値を表し、
   a、b及びcは、それぞれ、第一の定数、第二の定数、及び、第三の定数を表すことを特徴とする画像処理装置。
2. 前記階調補正手段は、
   前記第一の画像の各画素の画素値を定数倍した値を取得する第一の取得手段と、
   前記第二の画像の各画素の画素値を定数倍した値を取得する第二の取得手段と、
   前記第一の画像と前記第二の画像とで対応する画素同士について前記第一の積と前記第二の積を足し合わせる加算手段と、
   を備えることを特徴とする請求項１記載の画像処理装置。
3. 前記階調補正は、視覚モデルに基づいたものであることを特徴とする請求項１記載の画像処理装置。
4. 前記第三の定数がゼロであることを特徴とする請求項１記載の画像処理装置。
5. 前記エッジ保存平滑化した画像が、前記階調補正を行いたい画像を常用対数で対数変換して生成した画像であり、
   前記第一の定数は、０から２までの値であり、
   前記第二の定数は、前記第一の定数に対してマイナス２倍から０倍した値であり、
   前記第三の定数は、第一の定数をマイナス１倍したものに１から１．５までの値を加算した値であることを特徴とする請求項１に記載の画像処理装置。
6. 前記第三の定数は、前記対数tc'を前記第一の画像の画素値の範囲に近づける調整値であることを特徴とする請求項１に記載の画像処理装置。
7. 錐体応答に基づく階調補正を行う画像処理方法であって、
   階調補正後の画像の画素値の対数tc'を、数式
   tc' ＝ a・base'＋b・white'＋c
   を満たすよう演算する階調補正ステップを備え、
   ここで、base'は、階調補正前の画像に対しエッジ保存平滑化及び対数変換を行うことにより得られる第一の画像の画素値を表し、
   white'は、階調補正前の画像に対し平滑化及び対数変換を行なうことにより得られる第二の画像の画素値を表し、
   a、b及びcは、それぞれ、第一の定数、第二の定数、及び、第三の定数を表すことを特徴とする画像処理方法。
8. 錐体応答に基づく階調補正を行う画像処理プログラムであって、
   階調補正後の画像の画素値の対数tc'を、数式
   tc' ＝ a・base'＋b・white'＋c
   を満たすよう演算する階調補正手順をコンピュータに実行させ、
   ここで、base'は、階調補正前の画像に対しエッジ保存平滑化及び対数変換を行うことにより得られる第一の画像の画素値を表し、
   white'は、階調補正前の画像に対し平滑化及び対数変換を行なうことにより得られる第二の画像の画素値を表し、
   a、b及びcは、それぞれ、第一の定数、第二の定数、及び、第三の定数を表すことを特徴とする画像処理プログラム。