JP3644421B2 - Image conversion method and image conversion apparatus - Google Patents

Image conversion method and image conversion apparatus Download PDF

Info

Publication number
JP3644421B2
JP3644421B2 JP2001265895A JP2001265895A JP3644421B2 JP 3644421 B2 JP3644421 B2 JP 3644421B2 JP 2001265895 A JP2001265895 A JP 2001265895A JP 2001265895 A JP2001265895 A JP 2001265895A JP 3644421 B2 JP3644421 B2 JP 3644421B2
Authority
JP
Japan
Prior art keywords
luminance
image
range
low
coefficient
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
Application number
JP2001265895A
Other languages
Japanese (ja)
Other versions
JP2003076984A (en
Inventor
峰樹 曽我
啓一 山田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toyota Central R&D Labs Inc
Original Assignee
Toyota Central R&D Labs Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Toyota Central R&D Labs Inc filed Critical Toyota Central R&D Labs Inc
Priority to JP2001265895A priority Critical patent/JP3644421B2/en
Priority to US10/232,692 priority patent/US7065257B2/en
Publication of JP2003076984A publication Critical patent/JP2003076984A/en
Application granted granted Critical
Publication of JP3644421B2 publication Critical patent/JP3644421B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Image Processing (AREA)
  • Picture Signal Circuits (AREA)
  • Facsimile Image Signal Circuits (AREA)

Description

【0001】
【発明の属する技術分野】
本発明は各画素が濃淡を示す輝度値を有する画像の、明るさダイナミックレンジを圧縮する画像変換方法及び画像変換装置に関する。尚、本願明細書において2次元画像の輝度の「空間周波数」とは、各ライン(各行又は各列)毎の1次元の輝度の周波数をも包含する表現とする。
【0002】
【従来の技術】
陰の領域及び日向の領域の両方を有する明暗差の大きな画像を撮像し、明るさダイナミックレンジの小さい表示装置に表示させる場合、画像中の暗部が黒つぶれしたり、明部が白飛びすることとなり、画像中の物体等の把握が困難となる。明部や暗部などの特定の輝度範囲に多くの階調を割り当てるような階調変換により、特定の輝度範囲のコントラストを向上させることはできるが、その他の輝度範囲のコントラストが低下するという問題が生じる。そこでHomomorphic Filteringと呼ばれる画像のダイナミックレンジの圧縮方法が知られている。Homomorphic Filteringにおいては、輝度の低空間周波数成分を除くか、又は抑制することにより陰日向の差や照明ムラ等の照明強度の情報を抑え、陰の領域及び日向の領域いずれにおいても画像中の物体等の把握を容易とするものである。図18の(a)、(b)にその概念を示す。図18の(a)においては1画像の画素の輝度を取り込み、対数変換(log、91)したのち高域濾波器(HPF、92)で高空間周波数成分による画像データとし、逆対数変換(log-1、93)するものである。また、図18の(b)においては1画像の画素の輝度を取り込み、対数変換(log、91)したのち、低域濾波器(LPF、94)で低空間周波数成分による画像データを得て、符号を反転させて低域濾波されていない対数変換(log、91)の出力である画像データと加算器95により加算したのち逆対数変換(log-1、93)するものである。
【0003】
【発明が解決しようとする課題】
Homomorphic Filteringにおいては「低空間周波数成分」が一律に抑制されることとなり、例えば大きな白い物体(又はその表面)を撮像した画像に適用するとその画面内で大きくなるような白い物体が一律に灰色化してしまう。例えば道路標示を含む画像についてHomomorphic Filteringを適用すると、白線や横断歩道が灰色化し、道路標示の識別が極めて困難になる場合がある。これは、低空間周波数成分が抑制されることにより、高周波の白い物体よりも低周波の白い物体が相対的に暗くなるためである。低空間周波数成分として、画像の縦又は横の1方向の周波数を抑制するような簡便な方法を取る場合、特にこの問題が顕著となる。例えば画像中に縦方向の白線と横方向の白線が存在するとき、横方向に周波数分析すると、横方向の白線のみが灰色化してしまう。これは極めて違和感の大きな画像変換となる。同様の理由により、低輝度の低周波成分が明化することも有るが、白色物体の灰色化に比べ、その違和感は小さい。また、一律に低周波成分が抑制されることにより、低輝度部分においてノイズの目立つ画像となるという問題もある。これは、撮像画像のS/N比の小さい低輝度部分においても一律に高周波成分が相対的に強調されるためである。ダイナミックレンジの圧縮された画像の平均輝度を大きくして表示する場合、低輝度部分が明るく表示されるため、ノイズがより目立つ画像となる。
【0004】
そこで例えば医療用撮像装置等において、特許第1530832号、2509503号、2663189号公報記載の方法が提案されている。特許第1530832号に開示の方法は、画像の濃度が大きくなるほど高周波成分の強調度合いを大きくする画像の鮮鋭化方法であり、原画像をSorg、低周波画像をSus、強調係数をβとし、式S'=Sorg+β(Sorg+Sus)により、再生画像S'を得るものである。強調係数βを原画像又は低周波画像の濃度に応じて単調増加させることにより、ノイズの強調を防止している。この方法は鮮鋭化方法であるため、ノイズの強調を抑制することはできるが、原画像よりもノイズを抑制することはできない。また、ダイナミックレンジを圧縮することはできない。特許第2509503号に開示の方法は、画像信号の低空間周波数成分を用いてダイナミックレンジ圧縮用補正データを作成して、この補正データを用いて画像信号のダイナミックレンジを圧縮する方法である。ダイナミックレンジ圧縮用補正データは、画像信号に基づく信号値を閾値処理して、一方の場合は一定値とし、他方の場合は低空間周波数成分に基づいて求めるものである。つまり、画像中の明部又は暗部のいずれか一方は補正を行わないようにし、他方はダイナミックレンジを圧縮できるようにする方法である。この方法により、画像中の暗部のみでダイナミックレンジを圧縮し、明部ではダイナミックレンジを圧縮しないことにより、白い大きな物体の灰色化は防止できるが、低輝度部分でのノイズを抑制することはできない。また、明部と暗部の境界において周波数特性が急激に変化するため、儀輪郭が発生するという問題がある。特許第2663189号に開示の方法は、原画像の各画素においてその近傍の所定範囲内の画素での平均化により低周波画像を算出し、それに基づいて処理する方法である。具体的には、各画素において低周波画像の輝度を引数とする単調減少関数の出力を原画像の対応する画素の濃度値に加算することによりダイナミックレンジを圧縮する方法である。単調減少関数の特性の設定により、明部や暗部など所望の明るさ範囲のみを圧縮することができる方法である。単調減少関数の微係数を連続することにより儀輪郭の発生を防止することはできるが、原画像成分が常に一定で加算されているため、低輝度部分においてノイズを抑制することはできない。
【0005】
本発明は上記の課題を解決する為に成されたものであり、その目的は、暗部が黒つぶれすることなく、暗部のノイズが目立つことなく、画面内で大きな白い物体(又はその領域)が灰色化することのない、これらの特徴を全て併せ持つ画像変換方法又は画像変換装置を提供することである。更には、低空間周波数成分として画像の1方向(縦方向又は横方向)のような1次元の周波数のみを扱う場合においても、変換後の画像の違和感の少ない画像変換方法又は画像変換装置を提供することである。
【0006】
【課題を解決するための手段】
請求項1乃至請求項8に記載の手段は各々本願発明の画像変換方法であり、上記の課題を解決するために成されたものである。即ち、請求項1に記載の手段によれば、各画素が輝度値を有する元画像を処理して、明るさダイナミックレンジを圧縮する画像変換方法において、元画像の輝度の周波数の低い成分から成る低周波画像を抽出する低周波画像抽出工程と、元画像の画素の輝度Y1に対応する低周波画像の画素の輝度Y2から、2つの係数C1、C2を求める係数決定工程と、2つの係数C1、C2と、元画像の画素の輝度Y1と対応する低周波画像の画素の輝度Y2とから、対応する変換画像の画素の輝度Y3を式Y3=C1Y1+C2Y2により求める線形演算工程とを有し、低周波画像の画素の輝度Y2の関数である2つの係数C1、C2はいずれも連続関数であり、係数C1は正又は0であり、係数C2は、輝度Y2の最小値を含む連続した第1の範囲で正、0及び負の値をとる広義の単調減少関数であり、輝度Y2の最大値を含む連続した第2の範囲で負又は0の値をとる広義の単調増加関数であることを特徴とする。
【0007】
請求項2乃至請求項6に記載の手段は請求項1の従属請求項であって、各々2つの係数C1、C2の特徴を示すものである。即ち、請求項2に記載の手段によれば、2つの係数C1、C2の和が輝度Y2の全範囲のうちの少なくとも連続する第3の範囲において正の一定の値であることを特徴とする。請求項3に記載の手段によれば、2つの係数C1、C2の和が輝度Y2の全範囲において一定であることを特徴とする。また、請求項4に記載の手段によれば、輝度Y2の第3の範囲の内部である連続する第4の範囲において、係数C1が正の一定値であり、C2が負の一定値であることを特徴とする。
【0008】
請求項5に記載の手段によれば、輝度Y2の第1の範囲において輝度Y2の関数C2/C1が広義の単調減少関数であり、輝度Y2の第2の範囲において輝度Y2の関数C2/C1が広義の単調増加関数であることを特徴とする。また、請求項6に記載の手段によれば、輝度Y2の第4の範囲及びそれを越えた範囲の全てにおいて、係数C1が正の一定値であることを特徴とする。
【0009】
請求項7に記載の手段によれば、2つの係数C1とC2は、輝度Y2の最小値を含む連続した第5の範囲においてC1が0であり、輝度Y2の最大値を含む連続した第6の範囲においてC2が0であり、当該2つの範囲に挟まれた範囲の連続した少なくとも一部において係数C1が正の一定値、係数C2が負の一定値であることを特徴とする。
【0010】
請求項8に記載の手段は、請求項1乃至請求項7のいずれか1項に記載の画像変換方法において、低周波画像の抽出が、各画素ラインごとの1次元方向のみにて行われることを特徴とする。
【0011】
請求項9乃至請求項16に記載の手段は、各々本願発明の画像変換装置であり、それぞれ請求項1乃至請求項8に記載の画像変換方法に対応するものである。即ち、請求項9に記載の手段は、各画素が輝度値を有する元画像を処理して、明るさダイナミックレンジを圧縮する画像変換装置において、元画像の輝度の周波数の低い成分から成る低周波画像を抽出する低周波画像抽出手段と、低周波画像抽出手段から出力された低周波画像の画素の輝度Y2に対応して、各々係数C1、C2を求める第1、第2の係数決定手段と、当該第1、第2の係数決定手段により得られる係数C1、C2と、それに用いた低周波画像抽出手段から出力された低周波画像の画素の輝度Y2、及びそれに対応する元画像の画素の輝度Y1から、式Y3=C1Y1+C2Y2により変換画像の画素の輝度Y3を求める線形演算手段とを有し、第1、第2の係数決定手段における係数C1、C2について、低周波画像の画素の輝度Y2の関数である係数C1、C2はいずれも連続関数であり、係数C1は正又は0であり、係数C2は、輝度Y2の最小値を含む連続した第1の範囲で正、0及び負の値をとる広義の単調減少関数であり、輝度Y2の最大値を含む連続した第2の範囲で負又は0の値をとる広義の単調増加関数であることを特徴とする。
【0012】
請求項10に記載の手段は、係数C1、C2の和が前記輝度Y2の全範囲のうちの少なくとも連続する第3の範囲において正の一定の値であることを特徴とする。請求項11に記載の手段は、係数C1、C2の和が輝度Y2の全範囲において一定であることを特徴とする。請求項12に記載の手段は、輝度Y2の第3の範囲の内部である連続する第4の範囲において、係数C1が正の一定値であり、C2が負の一定値であることを特徴とする。
【0013】
請求項13に記載の手段は、輝度Y2の第1の範囲において輝度Y2の関数C2/C1が広義の単調減少関数であり、輝度Y2の第2の範囲において輝度Y2の関数C2/C1が広義の単調増加関数であることを特徴とする。請求項14に記載の手段は、輝度Y2の第4の範囲及びそれを越えた範囲の全てにおいて、係数C1が正の一定値であることを特徴とする。
【0014】
請求項15に記載の手段は、係数C1とC2は、輝度Y2の最小値を含む連続した第5の範囲においてC1が0であり、輝度Y2の最大値を含む連続した第6の範囲においてC2が0であり、当該2つの範囲に挟まれた範囲の連続した少なくとも一部において係数C1が正の一定値、係数C2が負の一定値であることを特徴とする。
【0015】
請求項16に記載の手段は、請求項9乃至請求項15のいずれか1項に記載の画像変換装置において、低周波画像抽出手段における低周波画像の抽出が、各画素ラインごとの1次元方向のみにて行われることを特徴とする。
【0016】
【作用及び発明の効果】
元画像の画素の輝度Y1に対応する低周波画像の画素の輝度Y2から、2つの係数C1、C2を求めることで、照明強度等の低周波成分の輝度に応じて、低周波成分の抑制量を変化させることができる。また、係数C1は連続関数で正又は0、係数C2は連続関数で当該輝度Y2の最小値を含む連続した第1の範囲で正、0及び負の値をとる広義の単調減少関数、輝度Y2の最大値を含む連続した第2の範囲で負又は0の値をとる広義の単調増加関数である。よって、輝度Y2の最小値を含む、係数C2が正の部分、即ち低輝度部分においては低周波成分が元画像よりも強調されるので、元画像よりもノイズ成分を抑制することが可能となる。係数C2が単調減少する負の部分においては、輝度Y2が小さいほど低周波成分の抑制量が小さいので、撮像画像のS/N比が小さくなるほど高周波成分が相対的に強調されなくなるのでノイズが目立つのを防止することができる。また、輝度Y2の最大値を含む、係数C2が負又は0で単調増加する部分においては、輝度Y2が大きくなるほど低周波成分の抑制量が小さくなる。よって高輝度の低周波成分の抑制量が小さくなるので、画面内で大きな面積となる白い物体の像が灰色化することを防止することが可能な画像変換方法又は画像変換装置とすることができる(請求項1、9)。
【0017】
また、2つの係数C1、C2の和が輝度Y2の全範囲のうちの少なくとも連続する第3の範囲において正の一定の値であるので、変換後の画像におけるゼロ空間周波数成分の伝達率が一定となるため、その範囲においてコントラストが低下するのを防止することが可能となる(請求項2、10)。更に、輝度Y2の全範囲において係数C1、C2の和が一定であれば、低周波成分のコントラスト低下防止が全輝度範囲において可能となる(請求項3、11)。また、輝度Y2の連続する第4の範囲において、係数C1が正の一定値であり係数C2が負の一定値であるので、当該範囲においては、低周波成分の抑制度合いが一定となるので、一定の圧縮率でダイナミックレンジを圧縮することができる。(請求項4、12)。
【0018】
また、輝度Y2の第1の範囲において輝度Y2の関数C2/C1が広義の単調減少関数であるので、第1の範囲(低輝度部分)においては輝度Y2が小さくなるほどより低周波成分が抑制度合いが小さくなる。従って、S/N比の小さい暗部においては、相対的に高周波成分が強調される度合いが小さくなるためノイズが目立つことはなく、且つ周波数特性が連続的に変化するため違和感の無い画像となる。また、輝度Y2の第2の範囲において輝度Y2の関数C2/C1が広義の単調増加関数であるので、第2の範囲(高輝度部分)においては輝度Y2が大きくなるほど低周波成分が抑制されなくなるので高輝度の大きな領域の暗化(白色領域の灰色化)を起こさない(請求項5、13)。尚、関数C2/C1が輝度Y2の第1の範囲において狭義の単調減少関数であり、輝度Y2の第2の範囲において狭義の単調増加関数であれば、これらの効果は尚一層顕著となる。更に、輝度Y2の第2の範囲及びそれを越えた範囲の全てにおいて、係数C1が正の一定値として、高輝度部分において輝度Y2が大きくなるほど低周波成分の抑制量を小さくすることにより、高輝度部分のコントラストを大きくすることができる(請求項6、14)。
【0019】
係数C1とC2は、輝度Y2の全範囲において和が正の一定値であり、輝度Y2の最小値を含む連続した第5の範囲においてC1が0、輝度Y2の最大値を含む連続した第6の範囲においてC2が0であり、当該2つの範囲に挟まれた範囲の連続した少なくとも一部において係数C1が正の一定値、係数C2が負の一定値であるようにすることによっても良い。即ち、輝度Y2の最小値を含む連続した第5の範囲において低周波成分のみを含んだ変換画像となるため、ノイズを大きく抑制することができ、低輝度部分のノイズが大きい場合に特に有効である。また、輝度Y2の最大値を含む連続した第6の範囲において低周波成分の抑制量は0となり、第6の範囲(高輝度部分)においては大きな領域の暗化(白色領域の灰色化)を起こさない。また、それらの間の領域では係数C1が正の一定値、係数C2が負の一定値であるので当該範囲においては一律に低周波成分が抑制されるので、一定の圧縮率でダイナミックレンジを圧縮することができる(請求項7、15)。
【0020】
上記各発明においては、低空間周波数として画像の1方向(縦方向又は横方向)のような1次元について低周波画像抽出を行う場合、低周波画像抽出を行わない方向のエッジ部分においては、低周波成分の抑制によるオーバーシュートやアンダーシュートの発生を抑えることができる。即ち、より自然な変換画像となる。更に低周波画像抽出を1次元で行うことで、簡便な画像変換方法、又は簡易で高速な画像変換装置とすることができる(請求項8、16)。
【0021】
【発明の実施の形態】
以下、本発明を具体的な実施例とともに説明する。尚、本発明は以下の実施例に限定されるものではない。
【0022】
〔第1実施例〕
図1は本発明の第1の実施例に係る画像変換装置10の構成を示すブロック図である。画像変換装置10は、対数変換器11、低周波画像抽出器12、遅延回路16、第1ルックアップテーブル13−1、第2ルックアップテーブル13−2、線形演算器14から成る。本実施例は請求項1、2、4乃至6、8乃至1012乃至14、16の具体的な実施例に当る。また、低周波画像抽出器12が低周波画像抽出手段に、第1ルックアップテーブル13−1及び第2ルックアップテーブル13−2がそれぞれ第1、第2の係数決定手段に、線形演算器14が線形演算手段に対応する。本実施例においては、図18のような逆対数変換93を用いないものを示した。また、本実施例においては、1ライン毎に周波数フィルタリングを行う、高速画像変換装置を示した。尚、本実施例においては図示しない制御装置により各構成要素が制御されるものである。
【0023】
撮像画像の各画素の輝度Y0が対数変換器11に出力され、各輝度Y0に対し、輝度Y1=(Y0maxlogY0)/logY0maxが計算され、低周波画像抽出器12と遅延回路16に出力される(図2)。ここでY0maxは、輝度Y0のとり得る値の最大値であり、この変換によりY0=Y0maxはY1=Y0maxに変換される。この、各画素が輝度Y1である画像を元画像と呼ぶ。各画素が輝度Y1である元画像は、次のようにして1ライン毎に変換される。
【0024】
まず、1ライン分の輝度Y1が、低周波画像抽出器12に読み込まれ、1ライン分の各画素の輝度をY2とする低周波画像が形成される(図3の(a))。ここにおいて、低周波画像の抽出は図3の(b)のような周波数特性を有するコンボルーションフィルタにより実現される。こうして各画素が輝度Y2である1ライン分の低周波画像が低周波画像抽出器12から出力される。
【0025】
次に、低周波画像抽出器12から各画素の輝度Y2が第1ルックアップテーブル13−1及び第2ルックアップテーブル13−2に順次入力される。第1ルックアップテーブル13−1及び第2ルックアップテーブル13−2においては、各々次の式による輝度Y2に対する係数C1、C2が記憶されている。ただし0<g<1とする。
【数1】

Figure 0003644421
【0026】
こうして、第1ルックアップテーブル13−1からは輝度Y2に対応する係数C1が、第2ルックアップテーブル13−2からは輝度Y2に対応する係数C2が線形演算器14に出力される。同時に、当該画素の輝度Y2が低周波画像抽出器12から線形演算器14に出力される。また、遅延回路16は、入力された輝度Y1を、低周波画像の出力に要する時間分遅延させ、低周波画像の画素に対応する元画像の画素の輝度Y1が線形演算器14に出力される。こうして、線形演算器14は、式Y3=C1Y1+C2Y2により輝度Y3を計算して出力する。各画素がこのような輝度Y3を有する変換画像は、各画素が輝度Y1である元画像に対し、低周波成分が抑制されたものである。
【0027】
各画素が輝度Y3の変換画像の特性について説明する。上記数1で示される、低周波画像の輝度Y2の関数である係数C1、C2及びそれらの商C2/C1のグラフを図4に示す。係数C1、C2は、Y2≦T2においてC1+C2が一定値g(ただし0<g<1)であるので、この範囲においてコントラストの低下を防止することができる。また、Y2≦T1においてC2が単調減少しているので、輝度Y2が小さくなるほど低周波成分の抑制量は小さくなる。従って、高周波成分であるノイズが強調されることを防止することができる。更にY2≦T0においてC2≧0なので、元画像よりも高周波成分が抑制され、ノイズを抑制することができる。特にT1≦Y2≦T2においてC2=g-1、C1=1で一定であるので、一律な低周波成分抑制となり、中輝度部分において一定の圧縮率でダイナミックレンジを圧縮することができる。また、Y2≧T2においては、Y2の増加とともにC1が一定値1で、C2が負でその絶対値が小さくなるので、Y2の増加とともに低周波成分の抑制量が小さくなる。よって、高輝度の低周波成分の暗化を抑えることができ、白い大きな物体(又はその領域、明部)が灰色化することがない。
【0028】
これを伝達特性として図5で更に説明する。C1が一定値1、C2が一定値g-1(ただし0<g<1)である範囲即ちY2≦T2においては、周波数伝達特性は図5の(a)のようになると言える。即ち、ゼロ空間周波数成分の伝達率がgであり、周波数が大きくなるほど1に近づくような、伝達率が徐々に大きくなる特性である。。
【0029】
次に、C1+C2が一定値gのまま、C2を正、0、負即ち、C1を小から大へと変化させると伝達特性は図5の(b)のようになる。即ち、C2が正の場合は相対的に高周波成分が抑制されることとなり、C2が負の場合は相対的に低周波成分が抑制されることとなる。C1が最大1であれば、図5の(a)と同様の特性となる。
【0030】
このように、本実施例の画像変換装置10は、低輝度部分においてノイズを抑制することができ、高輝度部分において画面内で大きな白い物体(又はその領域、明部)が灰色化することがない。中輝度部分では一定の圧縮率でダイナミックレンジが圧縮されるので、明暗差の大きな画像を黒つぶれや白飛びすることなく表示することができる。また、1つの参照値(Y2)によるルックアップテーブルを用いることでメモリを少量とし、また、高速化が可能である。また、1ライン毎の低周波成分抽出を行っているので、高速化及び装置全体の小型化が可能である。本実施例において、第2ルックアップテーブル13−2を、輝度Y2に対するC2Y2を記憶させて出力させるものとし、線形演算器14への輝度Y2の入力を省略した構成とすることも可能である。
【0031】
〔第2実施例〕
図6は、本発明の第2の実施例に係る画像変換装置20の構成を示すブロック図である。本実施例の画像変換装置20は、第1実施例の画像変換装置10に階調変換器15を加えたものである。階調変換器15は、線形演算器14の出力である変換画像の輝度Y3を階調変換してY4とするものである。階調変換器15はルックアップテーブルで構成することができる。
【0032】
階調変換器15の階調変換について図7乃至図9を用いて説明する。図7は、線形演算器14の出力である変換画像の輝度Y3において、低周波成分のダイナミックレンジが圧縮される様子を示しており、低周波画像の輝度Y2と変換画像の輝度Y3との関係を示したグラフである。今、元画像が低周波成分のみ有しているとして、元画像の輝度Y1と低周波画像の輝度Y2が一致すると考える。すると、変換画像の輝度Y3は(C1+C2)Y2とおくことができる。図7はこの仮定による低周波画像の輝度Y2と変換画像の輝度Y3との関係である。このグラフの傾きは、低周波画像の輝度Y2がT2未満のときg(0<g<1)であり、低周波画像の輝度Y2がT2を越えるときに1である。従って、低周波画像の輝度Y2がT2未満のときにダイナミックレンジはg倍に圧縮されていることがわかる。これは、低周波成分即ち、照明強度の成分のコントラストが低下したことを意味する。物体の反射率の成分は高周波成分であるのでそのコントラストは理想的には低下しないことになるが、実際には低周波成分も含まれるため、物体の反射率の成分のコントラストもわずかに低下する。そこで図8の実線に示すように、変換画像の輝度Y3がgT2以下のときにY4=Y3/gとなるよう階調変換する。このときのY4とY3の関係はF(Y2)=(C1+C2)Y2としたとき、Y4=F-1(Y3)である。つまり、図7に示したY2からY3への変換関数の逆関数である。この階調変換により、低輝度部分と高輝度部分で低周波成分のコントラストを等しくすることができる。つまり、低輝度部分において、わずかに低下した物体の反射率の成分のコントラストを回復することができる。また、図8の点線に示すように、輝度Y3がgT2以下のときにY4=Y3/gよりも傾きの小さい階調変換としても良い。この場合の階調変換の変換関数は、定数α(0<α<1)を用いてY4=αF-1(Y3)+(1-α)Y3と表すことができる。αを大きくするとダイナミックレンジの圧縮率は小さくなるが、低輝度部分での物体のコントラストは大きくすることができる。逆に、αを小さくすると、ダイナミックレンジの圧縮率は大きくなるが、低輝度部分での物体のコントラストは小さくなる。
【0033】
階調変換により、次のような効果もある。図9の(a)はステップ状に元画像の輝度Y1が変化している場合の、1ラインのプロフィールを、低周波画像の輝度Y2とともに示したものであり、横方向がラインの並び、縦方向が輝度(上が大)である。ここで輝度Y1及びY2が、わずかに閾値T2より低い部分に着目する。この部分の左端のステップ部分において、変換画像の輝度Y3は、図9の(b)のように、アンダーシュートUS3と、gT2を越えるオーバーシュートOS3を生じることとなる。このオーバーシュートOS3は、その左側のダイナミックレンジが圧縮されたステップ部分との差が比較的大きく視認されるため、擬似輪郭として把握されてしまう。そこで階調変換により、オーバーシュートOS3の左側のダイナミックレンジをひろげた画像の輝度Y4とすることで、ステップ部分との差が比較的小さくなるため、図9の(c)のように、オーバーシュートOS4は目立たないものとなる。尚、経験的に、アンダーシュートUS4は暗部となるため違和感は小さい。
【0034】
〔第3実施例〕
第1実施例の変形として、低周波画像の輝度Y2の関数である係数C1、C2として、次のようなものも採用できる。
【数2】
Figure 0003644421
【0035】
ここで、Y2maxはY2のとり得る値の最大値であり、γは1以上の定数である。数2は、Y3=C1Y1+C2Y2=C1(Y1+C2Y2/C1)と考え、C2Y2/C1のY2による微係数が連続となるようにしたものである。C1が一定であるときC2Y2/C1のY2による微係数の符号を反転させたものはダイナミックレンジ圧縮率であるため、輝度Y2に対してダイナミックレンジ圧縮率を連続的に変化させることができる。実際、係数C1、C2は連続であり、C2Y2/C1のY2=T2での微係数は左側も右側も(g-1)であり一致する。C2Y2/C1のY2≧T2での微係数は次のとおりである。
【数3】
Figure 0003644421
【0036】
数3から、Y2=T2でC2Y2/C1の右側微係数はg-1である。Y2=T2においてダイナミックレンジ圧縮率が滑らかに変化していることがわかる。γを大きくすると輝度Y2の変化に対してより急峻にダイナミックレンジ圧縮率が変化する。γ=1であるとき、Y2≧T2においてC2=g-1であり、輝度Y2が変化してもダイナミックレンジ圧縮率は変化しない。また、数3に示した微係数は輝度Y2の増加に対して単調に増加する、即ち、ダイナミックレンジ圧縮率は単調に減少し、Y2=Y2maxにおいてダイナミックレンジ圧縮率は0となる。従って、第1実施例と同様の効果を有する他、ダイナミックレンジ圧縮率が輝度に対して連続的に変化するため、より違和感の少ない変換画像を得ることができる。
【0037】
〔第4実施例〕
第1実施例の変形として、低周波画像の輝度Y2の関数である係数C1、C2として、次のようなものも採用できる。本実施例は請求項1乃至4、7乃至12、15、16の実施例に当たる。尚、h>1であり、第4の範囲はT1≦Y2≦T2、第5の範囲はY2≦T0、第6の範囲はY2≧T3である。
【数4】
Figure 0003644421
【0038】
数4のグラフを図11に示す。本実施例は低周波画像の輝度Y2の全範囲において係数C1、C2の和が一定であるので、低周波成分のコントラスト低下防止が全輝度において可能となる。また、その他も第1実施例と同様の効果を有する。本実施例の係数C1、C2は、輝度Y2の各範囲において直線的に変化するので、線形演算のみで算出することができる。第1実施例のルックアップテーブル13−1、13−2の替わりに、2つの演算器を設けて係数C1、C2を算出する場合に、簡単な演算で高速に算出することができる。
【0039】
〔第5実施例〕
第1実施例の変形として、低周波画像の輝度Y2の関数である係数C1、C2として、次のようなものも採用できる。本実施例は請求項1乃至16の実施例に当たる。尚、h>1である。
【数5】
Figure 0003644421
【0040】
数5のグラフを図12に示す。本実施例は、C2/C1が輝度Y2の各範囲において直線的に変化する例である。本実施例も、低輝度部分のノイズ抑制、高輝度部分の低周波成分の灰色化を防止するなど、第1実施例と第4実施例と同様の効果を併せ持つ。
【0041】
〔第6実施例〕
第1実施例、第3乃至第5実施例における、第1ルックアップテーブル13−1、第2ルックアップテーブル13−2、線形演算器14を1体のものとしても良い。更に第2実施例における階調変換器15をも含めて1体としても良い。構成を図13に示す。本実施例の画像変換装置30は、対数変換器11、低周波画像抽出器12、遅延回路16、並びに低空間周波数成分抑制演算器31若しくはルックアップテーブル31の構成である。低空間周波数成分抑制演算器31を構成要素とする場合は当該空間周波数成分抑制演算器31に係数C1、C2の演算機能並びに線形演算機能が組み込まれたものとする。或いは第2実施例と同様な効果を得るための構成とする場合は、階調変換機能をも組み込まれたものとする。また、ルックアップテーブル31を構成要素とする場合は、2つの参照値Y1、Y2から1つの値Y3或いはY4を出力するメモリ機能を組み込まれたものとする。
【0042】
〔第7、8、9実施例〕
第2実施例における階調変換器15を、対数変換する前の撮像画像の輝度Y0と変換画像の輝度Y3の1フレーム毎の平均値により変換度を調整する構成としても良い。これを第7実施例として画像変換装置40の構成を図14に示す。階調変換器43は、線形演算器14の出力Y3の他、撮像画像の輝度Y0の平均値を求める平均値算出器41の出力Y0aveと、変換画像の輝度Y3の平均値を求める平均値算出器42の出力Y3aveをも入力される。これにより階調変化後の輝度Y4の平均を所望の値にすることができる。ここから、平均値算出器41と平均値算出器42の一方を省略した構成でも良い。第8実施例として、第7実施例の画像変換装置40から平均値算出器42を省略した画像変換装置50の構成を図15に、第9実施例として、第7実施例の画像変換装置40から平均値算出器41を省略した画像変換装置60の構成を図16に示す。画像変換装置50の階調変換器51は、線形演算器14の出力Y3の他、撮像画像の輝度Y0の平均値を求める平均値算出器41の出力Y0aveにより階調変化後の輝度Y4を算出する。画像変換装置60の階調変換器61は、線形演算器14の出力Y3の他、変換画像の輝度Y3の平均値を求める平均値算出器42の出力Y3aveにより階調変化後の輝度Y4を算出する。これらいずれの実施例においても係数C1、C2の決定方法は数1、2、4、5のいずれでも良い。第7及び第8実施例では、対数変換する前の撮像画像の輝度Y0の平均値Y0aveを平均値算出器41により算出し、階調変換器43及び51の入力とする例を示したが、平均値算出器41を対数変換器11の後段に設け、対数変換後の元画像の輝度Y1の平均値を階調変換器43及び51の入力とすることも可能である。
【0043】
〔第10実施例〕
オートゲインコントローラ(AGC)付き撮像装置に接続する画像変換装置として、第1実施例の第1ルックアップテーブル、第2ルックアップテーブルに、撮像装置のAGCの情報を入力構成としても良い。第10実施例として画像変換装置70の構成を図17に示す。AGC81付き撮像装置80と組み合わせる画像変換装置70は、第1実施例の第1ルックアップテーブル13−1、第2ルックアップテーブル13−2に替えて第1ルックアップテーブル71−1、第2ルックアップテーブル71−2とし、各々AGC81からのゲイン値Gと低周波画像の輝度Y2の2つの参照値に対する係数C1、C2を記憶させておき、出力するものである。係数C1、C2の決定方法は数1、2、4、5のいずれでも良く、撮像装置80のAGCのゲイン値Gに応じて、次のように各式における閾値T0とT1を変化させる。閾値T0とT1は撮像装置のS/N比が小さくなる低輝度部分においてノイズを目立たなくするために、低周波成分の抑制量を小さくするY2の範囲を定めるための閾値である。この閾値を大きくし過ぎると高周波成分のコントラストが低下する輝度の範囲が広くなる。従って、閾値T0とT1は撮像画像のノイズレベルに応じた必要最小限の大きさであることが好ましい。一般に、撮像対象が暗い程AGCのゲイン値Gは大きくなり、撮像画像の輝度のノイズレベルは大きくなる。逆に、撮像対象が明るいほどAGCのゲイン値Gは小さくなり、撮像画像の輝度のノイズレベルは小さくなる。そこで閾値T0とT1をAGCのゲイン値Gに略比例させることにより、常に適切な輝度レベルでノイズの抑制をすることができ、ノイズの抑制が不足したり、暗部において必要以上にぼけた変換画像となることを防止することができる。
【図面の簡単な説明】
【図1】本発明の第1の実施例に係る画像変換装置の構成を示すブロック図。
【図2】対数変換器の作用を示す図。
【図3】低周波画像抽出器の作用を示す図。
【図4】本発明の第1の実施例に係る画像変換装置の線形変換における低周波画像の輝度Y2の関数である2つの係数の特性を示すグラフ図。
【図5】線形変換による周波数特性を示すグラフ図。
【図6】本発明の第2の実施例に係る画像変換装置の構成を示すブロック図。
【図7】第2実施例における階調変換器に入力される変換画像の特性を示す図。
【図8】第2実施例における階調変換器の作用を示す図。
【図9】第2実施例における階調変換器の作用を示す図。
【図10】本発明の第3の実施例に係る画像変換装置の線形変換における低周波画像の輝度Y2の関数である2つの係数の特性を示すグラフ図。
【図11】本発明の第4の実施例に係る画像変換装置の線形変換における低周波画像の輝度Y2の関数である2つの係数の特性を示すグラフ図。
【図12】本発明の第5の実施例に係る画像変換装置の線形変換における低周波画像の輝度Y2の関数である2つの係数の特性を示すグラフ図。
【図13】本発明の第6の実施例に係る画像変換装置の構成を示すブロック図。
【図14】本発明の第7の実施例に係る画像変換装置の構成を示すブロック図。
【図15】本発明の第8の実施例に係る画像変換装置の構成を示すブロック図。
【図16】本発明の第9の実施例に係る画像変換装置の構成を示すブロック図。
【図17】本発明の第10の実施例に係る画像変換装置の構成を示すブロック図。
【図18】 Homomorphic Filteringの2種類の構成を示すブロック図。
【符号の説明】
10、20、30、40、50、60、70 画像変換装置
11 対数変換器
12 低周波画像抽出器
13−1、71−1 第1ルックアップテーブル
13−2、71−2 第2ルックアップテーブル
14 線形演算器
15 階調変換器(ルックアップテーブル)
16 遅延回路
43、51、61 階調変換器
31 低空間周波数成分抑制演算器(ルックアップテーブル)
41、42 平均値算出器
80 AGC機能付撮像装置
81 AGC[0001]
BACKGROUND OF THE INVENTION
The present invention relates to an image conversion method and an image conversion apparatus for compressing a brightness dynamic range of an image in which each pixel has a luminance value indicating shading. In the present specification, the “spatial frequency” of the luminance of the two-dimensional image is also an expression including the frequency of the one-dimensional luminance for each line (each row or each column).
[0002]
[Prior art]
When taking an image with a large difference in brightness and darkness that has both a shadow area and a sunny area and displaying it on a display device with a low brightness dynamic range, the dark area in the image will be blacked out or the bright area will be whiteout. Thus, it is difficult to grasp an object or the like in the image. Contrast of a specific luminance range can be improved by gradation conversion that assigns many gradations to a specific luminance range such as a bright part and a dark part, but the contrast of other luminance ranges is reduced. Arise. Therefore, a method for compressing the dynamic range of an image called “Homomorphic Filtering” is known. In Homomorphic Filtering, by removing or suppressing low spatial frequency components of luminance, information on illumination intensity such as difference in shade and uneven illumination is suppressed, and objects in images in both shaded areas and sunny areas Etc. are easy to grasp. 18 (a) and 18 (b) show the concept. In FIG. 18 (a), the luminance of a pixel of one image is captured, logarithmic transformed (log, 91), then converted to image data by a high spatial frequency component by a high pass filter (HPF, 92), and inverse logarithmic transformed (log -1 93). Further, in FIG. 18B, after the luminance of one image pixel is captured and logarithmically converted (log, 91), image data based on low spatial frequency components is obtained with a low-pass filter (LPF, 94). The inverse logarithm transformation (log) is performed after adding the image data which is the output of the logarithmic transformation (log, 91) which has not been subjected to low-pass filtering by inverting the sign and the adder 95 -1 93).
[0003]
[Problems to be solved by the invention]
In Homomorphic Filtering, “low spatial frequency components” are uniformly suppressed. For example, when applied to an image of a large white object (or its surface), white objects that become large in the screen are grayed out uniformly. End up. For example, when Homomorphic Filtering is applied to an image including road markings, white lines and pedestrian crossings may become gray and identification of road markings may be extremely difficult. This is because the low-frequency white object becomes relatively darker than the high-frequency white object by suppressing the low spatial frequency component. This problem is particularly noticeable when a simple method is adopted in which the frequency in one direction of the vertical or horizontal direction is suppressed as the low spatial frequency component. For example, when a vertical white line and a horizontal white line exist in the image, if the frequency analysis is performed in the horizontal direction, only the horizontal white line is grayed out. This is an image conversion with a very uncomfortable feeling. For the same reason, a low-frequency low-frequency component may be brightened, but the uncomfortable feeling is small compared to the graying of a white object. There is also a problem that noise is conspicuous in a low luminance part by uniformly suppressing low frequency components. This is because the high-frequency component is uniformly relatively enhanced even in a low-luminance portion with a small S / N ratio of the captured image. When an image with a dynamic range compressed with a high average luminance is displayed, the low luminance portion is displayed brightly, so that the image becomes more conspicuous.
[0004]
Therefore, for example, methods described in Japanese Patent Nos. 1530832, 2509503, and 2663189 have been proposed for medical imaging devices and the like. The method disclosed in Japanese Patent No. 1530832 is an image sharpening method in which the enhancement degree of the high-frequency component increases as the image density increases. The original image is Sorg, the low-frequency image is Sus, and the enhancement coefficient is β. A reproduced image S ′ is obtained by S ′ = Sorg + β (Sorg + Sus). The enhancement of noise is prevented by monotonically increasing the enhancement coefficient β according to the density of the original image or the low-frequency image. Since this method is a sharpening method, noise enhancement can be suppressed, but noise cannot be suppressed more than the original image. Also, the dynamic range cannot be compressed. The method disclosed in Japanese Patent No. 2509503 is a method of creating correction data for dynamic range compression using low spatial frequency components of an image signal, and compressing the dynamic range of the image signal using this correction data. The dynamic range compression correction data is obtained by subjecting a signal value based on an image signal to a threshold value to be a constant value in one case, and based on a low spatial frequency component in the other case. That is, either the bright part or the dark part in the image is not corrected, and the other is a method in which the dynamic range can be compressed. By this method, by compressing the dynamic range only in the dark part of the image and not compressing the dynamic range in the bright part, graying of a large white object can be prevented, but noise in a low luminance part cannot be suppressed. . In addition, since the frequency characteristics change abruptly at the boundary between the bright part and the dark part, there is a problem that a ceremonial outline is generated. The method disclosed in Japanese Patent No. 2663189 is a method of calculating a low-frequency image by averaging the pixels in a predetermined range near each pixel of the original image, and processing based on the low-frequency image. Specifically, this is a method of compressing the dynamic range by adding the output of a monotone decreasing function with the luminance of the low-frequency image as an argument to each pixel to the density value of the corresponding pixel of the original image. This is a method in which only a desired brightness range such as a bright part and a dark part can be compressed by setting the characteristic of a monotonic decrease function. Although it is possible to prevent the appearance of a ceremonial outline by continuing the derivative of the monotone decreasing function, noise cannot be suppressed in a low-luminance portion because the original image components are always added at a constant value.
[0005]
The present invention has been made in order to solve the above-described problems. The purpose of the present invention is to make a large white object (or an area thereof) within the screen without causing dark areas to be blackened or causing dark noise to be noticeable. An object is to provide an image conversion method or an image conversion apparatus that has all of these characteristics without being grayed out. Furthermore, even when only a one-dimensional frequency such as one direction (vertical direction or horizontal direction) of an image is handled as a low spatial frequency component, an image conversion method or an image conversion apparatus with less sense of incongruity of the converted image is provided. It is to be.
[0006]
[Means for Solving the Problems]
Each of the means described in claims 1 to 8 is an image conversion method according to the present invention, and is made to solve the above-mentioned problems. That is, according to the first aspect of the present invention, in the image conversion method for processing the original image in which each pixel has a luminance value and compressing the brightness dynamic range, the image conversion method includes components having a low luminance frequency of the original image. Low-frequency image extraction process for extracting low-frequency images and the luminance Y of the pixels of the original image 1 Low-frequency image pixel brightness Y corresponding to 2 From the two coefficients C 1 , C 2 The coefficient determination process to find the two coefficients C 1 , C 2 And the luminance Y of the pixel of the original image 1 And luminance Y of the corresponding low frequency image pixel 2 And the luminance Y of the pixel of the corresponding converted image Three The formula Y Three = C 1 Y 1 + C 2 Y 2 And a linear calculation process obtained by 2 Two coefficients C that are functions of 1 , C 2 Are continuous functions, and the coefficient C 1 Is positive or zero and the coefficient C 2 Is the luminance Y 2 Is a monotonically decreasing function in a broad sense that takes positive, zero and negative values in a continuous first range including the minimum value of 2 It is a monotonically increasing function in a broad sense that takes a negative or zero value in a continuous second range including the maximum value of.
[0007]
The means according to claim 2 to claim 6 are dependent claims of claim 1, each comprising two coefficients C 1 , C 2 It shows the characteristics of. That is, according to the means of claim 2, two coefficients C 1 , C 2 Is the sum of brightness Y 2 It is characterized by being a positive constant value in at least a third continuous range of all the ranges. According to the means of claim 3, two coefficients C 1 , C 2 Is the sum of brightness Y 2 It is characterized in that it is constant over the entire range. Further, according to the means of claim 4, the luminance Y 2 In the fourth continuous range that is inside the third range of 1 Is a constant positive value and C 2 Is a negative constant value.
[0008]
According to the means of claim 5, the luminance Y 2 Luminance Y in the first range of 2 Function C 2 / C 1 Is a monotonically decreasing function in a broad sense, and luminance Y 2 Luminance Y in the second range of 2 Function C 2 / C 1 Is a monotonically increasing function in a broad sense. Further, according to the means of claim 6, the luminance Y 2 The coefficient C in all of the fourth range and beyond 1 Is a positive constant value.
[0009]
According to the measure of claim 7, two coefficients C 1 And C 2 Is the luminance Y 2 In a fifth continuous range including the minimum of 1 Is 0, luminance Y 2 C in a continuous sixth range containing the maximum value of 2 Is 0, and the coefficient C in at least a part of the continuous range between the two ranges 1 Is a positive constant value, coefficient C 2 Is a negative constant value.
[0010]
According to an eighth aspect of the present invention, in the image conversion method according to any one of the first to seventh aspects, the low-frequency image is extracted only in the one-dimensional direction for each pixel line. It is characterized by.
[0011]
The means described in claims 9 to 16 is an image conversion apparatus according to the present invention, and corresponds to the image conversion method described in claims 1 to 8, respectively. That is, according to the ninth aspect of the present invention, in the image conversion apparatus that processes the original image in which each pixel has a luminance value and compresses the brightness dynamic range, the low frequency composed of the low frequency component of the luminance of the original image. Low frequency image extraction means for extracting an image, and luminance Y of a pixel of the low frequency image output from the low frequency image extraction means 2 Corresponding to each coefficient C 1 , C 2 And the coefficient C obtained by the first and second coefficient determining means and the first and second coefficient determining means. 1 , C 2 And the luminance Y of the pixel of the low-frequency image output from the low-frequency image extraction means used therefor 2 , And the corresponding pixel brightness Y of the original image 1 From the formula Y Three = C 1 Y 1 + C 2 Y 2 By the brightness Y of the pixel of the converted image Three Linear calculation means for obtaining the coefficient C in the first and second coefficient determination means 1 , C 2 The luminance Y of the pixel of the low frequency image 2 Coefficient C which is a function of 1 , C 2 Are continuous functions, and the coefficient C 1 Is positive or zero and the coefficient C 2 Is the luminance Y 2 Is a monotonically decreasing function in a broad sense that takes positive, zero and negative values in a continuous first range including the minimum value of 2 It is a monotonically increasing function in a broad sense that takes a negative or zero value in a continuous second range including the maximum value of.
[0012]
The means according to claim 10 comprises a factor C 1 , C 2 Is the brightness Y 2 It is characterized by being a positive constant value in at least a third continuous range of all the ranges. The means according to claim 11 comprises a coefficient C 1 , C 2 Is the sum of brightness Y 2 It is characterized in that it is constant over the entire range. The means according to claim 12 is the luminance Y 2 In the fourth continuous range that is inside the third range of 1 Is a constant positive value and C 2 Is a negative constant value.
[0013]
The means according to claim 13 is the luminance Y 2 Luminance Y in the first range of 2 Function C 2 / C 1 Is a monotonically decreasing function in a broad sense, and luminance Y 2 Luminance Y in the second range of 2 Function C 2 / C 1 Is a monotonically increasing function in a broad sense. The means according to claim 14 is the luminance Y 2 The coefficient C in all of the fourth range and beyond 1 Is a positive constant value.
[0014]
The means according to claim 15 comprises a coefficient C 1 And C 2 Is the luminance Y 2 In a fifth continuous range including the minimum of 1 Is 0, luminance Y 2 C in a continuous sixth range containing the maximum value of 2 Is 0, and the coefficient C in at least a part of the continuous range between the two ranges 1 Is a positive constant value, coefficient C 2 Is a negative constant value.
[0015]
The means according to claim 16 is the image conversion apparatus according to any one of claims 9 to 15, wherein the low-frequency image extraction means performs the one-dimensional direction for each pixel line. It is characterized by being performed only by.
[0016]
[Operation and effect of the invention]
Original image pixel brightness Y 1 Low-frequency image pixel brightness Y corresponding to 2 From the two coefficients C 1 , C 2 Therefore, the suppression amount of the low frequency component can be changed according to the luminance of the low frequency component such as the illumination intensity. The coefficient C 1 Is a continuous function positive or zero, coefficient C 2 Is the continuous function Y 2 A monotonically decreasing function in a broad sense that takes positive, zero, and negative values in a continuous first range including the minimum value of 2 Is a monotonically increasing function in a broad sense that takes a negative or zero value in a second continuous range including the maximum value of. Therefore, luminance Y 2 Coefficient C, including the minimum of 2 Since the low frequency component is emphasized more than the original image in the positive portion, that is, the low luminance portion, the noise component can be suppressed more than the original image. Coefficient C 2 In the negative part where the monotonously decreases, the luminance Y 2 Since the suppression amount of the low frequency component is smaller as the S is smaller, the high frequency component is relatively less emphasized as the S / N ratio of the captured image is smaller, so that noise can be prevented from being noticeable. Luminance Y 2 Coefficient C including the maximum value of 2 When Y is negative or increases monotonically with 0, luminance Y 2 As the value increases, the amount of suppression of low frequency components decreases. Therefore, since the amount of suppression of the low-frequency component with high luminance is reduced, it is possible to provide an image conversion method or image conversion apparatus that can prevent a white object image having a large area on the screen from being grayed out. (Claims 1 and 9).
[0017]
Two coefficients C 1 , C 2 Is the sum of brightness Y 2 Since the transmission rate of the zero spatial frequency component in the converted image is constant in at least the third continuous range in the entire range, the contrast decreases in that range. This can be prevented (claims 2 and 10). In addition, brightness Y 2 Coefficient C over the entire range of 1 , C 2 If the sum of the two is constant, it is possible to prevent the low-frequency component from lowering the contrast in the entire luminance range (claims 3 and 11). Luminance Y 2 In the fourth continuous range of 1 Is a positive constant value and the coefficient C 2 Is a negative constant value, and in this range, the degree of suppression of low-frequency components is constant, so that the dynamic range can be compressed at a constant compression rate. (Claims 4 and 12).
[0018]
Luminance Y 2 Luminance Y in the first range of 2 Function C 2 / C 1 Is a monotonically decreasing function in a broad sense, so that the luminance Y in the first range (low luminance portion) 2 As the value of becomes smaller, the degree of suppression of lower frequency components becomes smaller. Therefore, in a dark part where the S / N ratio is small, the degree of high frequency component enhancement is relatively small, so that noise does not stand out, and the frequency characteristics continuously change, so that the image does not feel strange. Luminance Y 2 Luminance Y in the second range of 2 Function C 2 / C 1 Is a monotonically increasing function in a broad sense, so that the luminance Y in the second range (high luminance part) 2 Since the low-frequency component is not suppressed as the value of becomes larger, darkening of a high-brightness large region (greying of the white region) does not occur (claims 5 and 13). Function C 2 / C 1 Is Y 2 Is a monotonically decreasing function in a narrow sense in the first range of 2 In the second range, these effects become even more prominent if they are monotonically increasing functions in a narrow sense. In addition, brightness Y 2 The coefficient C in all of the second range and beyond 1 Is a positive constant value, and brightness Y 2 By reducing the amount of suppression of low-frequency components as the value of increases, the contrast of the high-luminance portion can be increased (claims 6 and 14).
[0019]
Coefficient C 1 And C 2 Is the luminance Y 2 The sum is positive and constant in the entire range of brightness Y 2 In a fifth continuous range including the minimum of 1 Is 0, brightness Y 2 C in a continuous sixth range containing the maximum value of 2 Is 0, and the coefficient C in at least a part of the continuous range between the two ranges 1 Is a positive constant value, coefficient C 2 May be a negative constant value. That is, luminance Y 2 Since the converted image includes only the low-frequency component in the continuous fifth range including the minimum value of the noise, the noise can be largely suppressed, and is particularly effective when the noise in the low luminance portion is large. Luminance Y 2 The suppression amount of the low-frequency component is 0 in the continuous sixth range including the maximum value of, and no darkening of the large region (greying of the white region) occurs in the sixth range (high luminance portion). Also, in the area between them, the coefficient C 1 Is a positive constant value, coefficient C 2 Since N is a negative constant value, low frequency components are uniformly suppressed in this range, so that the dynamic range can be compressed at a constant compression rate (claims 7 and 15).
[0020]
In each of the above inventions, when low-frequency image extraction is performed for one dimension such as one direction (vertical direction or horizontal direction) of the image as a low spatial frequency, the edge portion in the direction where low-frequency image extraction is not performed is low. Generation of overshoot and undershoot due to suppression of frequency components can be suppressed. That is, the converted image becomes more natural. Furthermore, by performing low-frequency image extraction in one dimension, a simple image conversion method or a simple and high-speed image conversion device can be obtained.
[0021]
DETAILED DESCRIPTION OF THE INVENTION
Hereinafter, the present invention will be described together with specific examples. The present invention is not limited to the following examples.
[0022]
[First embodiment]
FIG. 1 is a block diagram showing a configuration of an image conversion apparatus 10 according to the first embodiment of the present invention. The image conversion apparatus 10 includes a logarithmic converter 11, a low-frequency image extractor 12, a delay circuit 16, a first lookup table 13-1, a second lookup table 13-2, and a linear calculator 14. This embodiment corresponds to specific embodiments of claims 1, 2, 4 to 6, 8 to 1012 to 14, and 16. Further, the low-frequency image extractor 12 serves as the low-frequency image extractor, the first lookup table 13-1 and the second lookup table 13-2 serve as the first and second coefficient determination units, respectively, and the linear calculator 14. Corresponds to the linear operation means. In the present embodiment, the one not using the inverse logarithmic transformation 93 as shown in FIG. 18 is shown. In the present embodiment, a high-speed image conversion apparatus that performs frequency filtering for each line is shown. In the present embodiment, each component is controlled by a control device (not shown).
[0023]
Luminance Y of each pixel of the captured image 0 Is output to the logarithmic converter 11 and each luminance Y 0 Against luminance Y 1 = (Y 0max logY 0 ) / LogY 0max Is calculated and output to the low-frequency image extractor 12 and the delay circuit 16 (FIG. 2). Where Y 0max Is the luminance Y 0 Is the maximum value that can be 0 = Y 0max Y 1 = Y 0max Is converted to Each pixel has luminance Y 1 Is called an original image. Each pixel has luminance Y 1 The original image is converted for each line as follows.
[0024]
First, luminance Y for one line 1 Is read into the low-frequency image extractor 12, and the luminance of each pixel for one line is set to Y 2 Is formed ((a) in FIG. 3). Here, the extraction of the low-frequency image is realized by a convolution filter having frequency characteristics as shown in FIG. Thus, each pixel has luminance Y 2 A low frequency image for one line is output from the low frequency image extractor 12.
[0025]
Next, the luminance Y of each pixel is output from the low frequency image extractor 12. 2 Are sequentially input to the first lookup table 13-1 and the second lookup table 13-2. In the first look-up table 13-1 and the second look-up table 13-2, the luminance Y according to the following equation is used. 2 Coefficient C for 1 , C 2 Is remembered. However, 0 <g <1.
[Expression 1]
Figure 0003644421
[0026]
Thus, from the first lookup table 13-1, the luminance Y 2 Coefficient C corresponding to 1 However, from the second lookup table 13-2, the luminance Y 2 Coefficient C corresponding to 2 Is output to the linear calculator 14. At the same time, the luminance Y of the pixel 2 Is output from the low-frequency image extractor 12 to the linear calculator 14. In addition, the delay circuit 16 receives the input luminance Y 1 Is delayed by the time required to output the low-frequency image, and the luminance Y of the pixel of the original image corresponding to the pixel of the low-frequency image 1 Is output to the linear calculator 14. In this way, the linear computing unit 14 is converted to the equation Y Three = C 1 Y 1 + C 2 Y 2 Brightness Y Three Is calculated and output. Each pixel has such brightness Y Three Each pixel has a luminance Y 1 In the original image, the low frequency component is suppressed.
[0027]
Each pixel has luminance Y Three The characteristics of the converted image will be described. Luminance Y of the low-frequency image represented by the above equation 1 2 Coefficient C which is a function of 1 , C 2 And their quotient C 2 / C 1 This graph is shown in FIG. Coefficient C 1 , C 2 Y 2 ≦ T 2 At C 1 + C 2 Is a constant value g (where 0 <g <1), it is possible to prevent a decrease in contrast in this range. Y 2 ≦ T 1 At C 2 Is monotonously decreasing, so brightness Y 2 The smaller the is, the smaller the suppression amount of the low frequency component is. Therefore, it is possible to prevent noise that is a high-frequency component from being emphasized. Y 2 ≦ T 0 At C 2 Since ≧ 0, high frequency components are suppressed more than the original image, and noise can be suppressed. T in particular 1 ≦ Y 2 ≦ T 2 At C 2 = G-1, C 1 Since = 1 is constant, the low frequency component is uniformly suppressed, and the dynamic range can be compressed at a constant compression rate in the middle luminance portion. Y 2 ≧ T 2 In Y 2 C with increasing 1 Is constant value 1, C 2 Is negative and its absolute value is small, so Y 2 As the frequency increases, the amount of suppression of low-frequency components decreases. Therefore, darkening of high-luminance low-frequency components can be suppressed, and a large white object (or its area, bright part) does not become gray.
[0028]
This will be further described with reference to FIG. 5 as a transfer characteristic. C 1 Is constant value 1, C 2 Is a constant value g-1 (where 0 <g <1), ie Y 2 ≦ T 2 In FIG. 5, it can be said that the frequency transfer characteristic is as shown in FIG. That is, the transmission rate of the zero spatial frequency component is g, and the transmission rate gradually increases so as to approach 1 as the frequency increases. .
[0029]
Next, C 1 + C 2 Is constant value g, C 2 Is positive, 0, negative, ie C 1 5 is changed from small to large, the transfer characteristic is as shown in FIG. That is, C 2 When is positive, the high frequency component is relatively suppressed, and C 2 When is negative, the low-frequency component is relatively suppressed. C 1 If is a maximum of 1, the characteristics are the same as in FIG.
[0030]
As described above, the image conversion apparatus 10 according to the present embodiment can suppress noise in the low luminance portion, and a large white object (or a region or a bright portion thereof) may be grayed in the screen in the high luminance portion. Absent. Since the dynamic range is compressed at a constant compression rate in the medium luminance portion, an image with a large contrast can be displayed without being blacked out or overexposed. One reference value (Y 2 ) Can be used in a small amount of memory and can be speeded up. Further, since the low-frequency component extraction is performed for each line, it is possible to increase the speed and reduce the size of the entire apparatus. In the present embodiment, the second look-up table 13-2 is set as the luminance Y 2 Against C 2 Y 2 Is stored and output, and the luminance Y to the linear calculator 14 is 2 It is also possible to adopt a configuration in which the input of is omitted.
[0031]
[Second Embodiment]
FIG. 6 is a block diagram showing the configuration of the image conversion apparatus 20 according to the second embodiment of the present invention. The image conversion apparatus 20 of this embodiment is obtained by adding a gradation converter 15 to the image conversion apparatus 10 of the first embodiment. The gradation converter 15 is the luminance Y of the converted image that is the output of the linear calculator 14. Three Y Four It is what. The gradation converter 15 can be configured by a lookup table.
[0032]
The gradation conversion of the gradation converter 15 will be described with reference to FIGS. FIG. 7 shows the luminance Y of the converted image that is the output of the linear calculator 14. Three Shows how the dynamic range of the low-frequency component is compressed, and the brightness Y of the low-frequency image 2 And converted image brightness Y Three It is the graph which showed the relationship. Now, assuming that the original image has only low frequency components, the luminance Y of the original image 1 And low-frequency image brightness Y 2 Are considered to match. Then, the brightness Y of the converted image Three (C 1 + C 2 ) Y 2 It can be said. FIG. 7 shows the luminance Y of the low frequency image based on this assumption. 2 And converted image brightness Y Three Relationship. The slope of this graph is the luminance Y of the low-frequency image 2 T 2 G (0 <g <1) and the luminance Y of the low-frequency image 2 T 2 1 when crossing over. Therefore, the luminance Y of the low frequency image 2 T 2 It can be seen that the dynamic range is compressed to g times when the value is less than. This means that the contrast of the low frequency component, that is, the component of the illumination intensity has decreased. Since the reflectance component of the object is a high-frequency component, its contrast will not ideally decrease. However, since the low-frequency component is actually included, the contrast of the reflectance component of the object slightly decreases. . Therefore, as shown by the solid line in FIG. Three GT 2 Y when Four = Y Three Tone conversion is performed to obtain / g. Y at this time Four And Y Three The relationship is F (Y 2 ) = (C 1 + C 2 ) Y 2 Y Four = F -1 (Y Three ). That is, Y shown in FIG. 2 To Y Three This is the inverse function of the conversion function. By this gradation conversion, the contrast of the low frequency component can be made equal in the low luminance portion and the high luminance portion. That is, it is possible to recover the contrast of the reflectance component of the object slightly lowered in the low luminance portion. Further, as shown by the dotted line in FIG. Three GT 2 Y when Four = Y Three Gradation conversion with a slope smaller than / g may be used. In this case, the conversion function of the gradation conversion is Y using a constant α (0 <α <1). Four = αF -1 (Y Three ) + (1-α) Y Three It can be expressed as. Increasing α reduces the compression ratio of the dynamic range, but can increase the contrast of the object in the low luminance part. On the other hand, if α is reduced, the compression ratio of the dynamic range is increased, but the contrast of the object in the low luminance portion is decreased.
[0033]
The gradation conversion has the following effects. FIG. 9A shows the luminance Y of the original image in a stepped manner. 1 The profile of one line when the value of 2 The horizontal direction is line arrangement, and the vertical direction is luminance (upper is large). Where brightness Y 1 And Y 2 But slightly threshold T 2 Focus on the lower part. In the leftmost step part of this part, the brightness Y of the converted image Three Is undershoot US as shown in FIG. Three And gT 2 Overshoot OS beyond Three Will result. This overshoot OS Three Is recognized as a pseudo contour because a difference between the left dynamic range and the compressed step portion is relatively large. Therefore, overshoot OS by gradation conversion. Three Brightness Y of the image with a wide dynamic range on the left side of Four Since the difference from the step portion becomes relatively small, the overshoot OS as shown in FIG. Four Is inconspicuous. In addition, empirically, undershoot US Four Since it becomes a dark part, the discomfort is small.
[0034]
[Third embodiment]
As a modification of the first embodiment, the luminance Y of the low-frequency image 2 Coefficient C which is a function of 1 , C 2 The following can also be adopted.
[Expression 2]
Figure 0003644421
[0035]
Where Y 2max Y 2 Γ is a constant of 1 or more. Equation 2 is Y Three = C 1 Y 1 + C 2 Y 2 = C 1 (Y 1 + C 2 Y 2 / C 1 ), C 2 Y 2 / C 1 Y 2 The derivative by is made continuous. C 1 C when is constant 2 Y 2 / C 1 Y 2 The value obtained by inverting the sign of the derivative by is the dynamic range compression ratio. 2 In contrast, the dynamic range compression rate can be continuously changed. In fact, the coefficient C 1 , C 2 Is continuous and C 2 Y 2 / C 1 Y 2 = T 2 The derivative of at (g-1) is the same on both the left and right sides. C 2 Y 2 / C 1 Y 2 ≧ T 2 The derivative at is as follows.
[Equation 3]
Figure 0003644421
[0036]
From Equation 3, Y 2 = T 2 At C 2 Y 2 / C 1 The right derivative of is g-1. Y 2 = T 2 It can be seen that the dynamic range compression rate changes smoothly at. Increase γ to increase brightness Y 2 The dynamic range compression rate changes more steeply with respect to the change in the value. Y when γ = 1 2 ≧ T 2 At C 2 = g-1 and brightness Y 2 The dynamic range compression ratio does not change even if changes. In addition, the derivative shown in Equation 3 is the luminance Y 2 Increases monotonically with increasing, i.e., the dynamic range compression ratio decreases monotonically, Y 2 = Y 2max The dynamic range compression rate is zero. Therefore, in addition to having the same effect as in the first embodiment, the dynamic range compression rate continuously changes with respect to the luminance, so that a converted image with a less uncomfortable feeling can be obtained.
[0037]
[Fourth embodiment]
As a modification of the first embodiment, the luminance Y of the low-frequency image 2 Coefficient C which is a function of 1 , C 2 The following can also be adopted. The present embodiment corresponds to the embodiments of claims 1 to 4, 7 to 12, 15, and 16. Note that h> 1 and the fourth range is T 1 ≦ Y 2 ≦ T 2 , The fifth range is Y 2 ≦ T 0 , The sixth range is Y 2 ≧ T Three It is.
[Expression 4]
Figure 0003644421
[0038]
The graph of Formula 4 is shown in FIG. In this embodiment, the luminance Y of the low frequency image 2 Coefficient C over the entire range of 1 , C 2 Since the sum of the two is constant, it is possible to prevent low frequency component contrast reduction at all luminance levels. The other effects are the same as those of the first embodiment. Coefficient C of this example 1 , C 2 Is the luminance Y 2 Since it changes linearly in each range, it can be calculated only by linear calculation. Instead of the look-up tables 13-1 and 13-2 of the first embodiment, two arithmetic units are provided and the coefficient C 1 , C 2 Can be calculated at high speed with a simple calculation.
[0039]
[Fifth embodiment]
As a modification of the first embodiment, the luminance Y of the low-frequency image 2 Coefficient C which is a function of 1 , C 2 The following can also be adopted. This embodiment corresponds to the embodiments of claims 1 to 16. Note that h> 1.
[Equation 5]
Figure 0003644421
[0040]
The graph of Formula 5 is shown in FIG. In this example, C 2 / C 1 Is Y 2 It is an example which changes linearly in each range. This embodiment also has the same effects as the first embodiment and the fourth embodiment, such as suppressing noise in the low luminance portion and preventing graying out of the low frequency component in the high luminance portion.
[0041]
[Sixth embodiment]
The first lookup table 13-1, the second lookup table 13-2, and the linear computing unit 14 in the first embodiment and the third to fifth embodiments may be one. Furthermore, it is good also as one body including the gradation converter 15 in 2nd Example. The configuration is shown in FIG. The image conversion apparatus 30 of this embodiment has a configuration of a logarithmic converter 11, a low-frequency image extractor 12, a delay circuit 16, and a low spatial frequency component suppression calculator 31 or a lookup table 31. When the low spatial frequency component suppression calculator 31 is a constituent element, the spatial frequency component suppression calculator 31 has a coefficient C 1 , C 2 It is assumed that the calculation function and the linear calculation function are incorporated. Alternatively, in the case of a configuration for obtaining the same effect as in the second embodiment, it is assumed that a gradation conversion function is also incorporated. When the lookup table 31 is a component, two reference values Y 1 , Y 2 1 value from Three Or Y Four It is assumed that a memory function for outputting is incorporated.
[0042]
[Seventh, eighth and ninth embodiments]
The gradation Y of the captured image before logarithmic conversion is applied to the gradation converter 15 in the second embodiment. 0 And converted image brightness Y Three The degree of conversion may be adjusted by the average value for each frame. FIG. 14 shows the configuration of the image conversion apparatus 40 as a seventh embodiment. The gradation converter 43 outputs the output Y of the linear calculator 14. Three In addition, the brightness Y of the captured image 0 The output Y of the average value calculator 41 for calculating the average value of 0ave And brightness Y of the converted image Three The output Y of the average value calculator 42 for calculating the average value of 3ave Is also entered. As a result, brightness Y after gradation change Four Can be set to a desired value. From here, the structure which abbreviate | omitted one of the average value calculator 41 and the average value calculator 42 may be sufficient. As an eighth embodiment, FIG. 15 shows a configuration of an image conversion apparatus 50 in which the average value calculator 42 is omitted from the image conversion apparatus 40 of the seventh embodiment, and as an ninth embodiment, the image conversion apparatus 40 of the seventh embodiment. FIG. 16 shows a configuration of the image conversion apparatus 60 in which the average value calculator 41 is omitted. The gradation converter 51 of the image conversion device 50 outputs the output Y of the linear calculator 14. Three In addition, the brightness Y of the captured image 0 The output Y of the average value calculator 41 for calculating the average value of 0ave Luminance Y after gradation change by Four Is calculated. The gradation converter 61 of the image conversion device 60 outputs the output Y of the linear calculator 14. Three In addition, the brightness Y of the converted image Three The output Y of the average value calculator 42 for calculating the average value of 3ave Luminance Y after gradation change by Four Is calculated. In any of these examples, the coefficient C 1 , C 2 Any of the formulas 1, 2, 4, and 5 may be used. In the seventh and eighth embodiments, the luminance Y of the captured image before logarithmic conversion 0 Mean value Y 0ave The average value calculator 41 is used as an input to the gradation converters 43 and 51. However, the average value calculator 41 is provided in the subsequent stage of the logarithmic converter 11, and the luminance of the original image after logarithmic conversion is shown. Y 1 Can be used as an input to the gradation converters 43 and 51.
[0043]
[Tenth embodiment]
As an image conversion apparatus connected to an image pickup apparatus with an auto gain controller (AGC), AGC information of the image pickup apparatus may be input to the first look-up table and the second look-up table of the first embodiment. FIG. 17 shows the configuration of an image conversion apparatus 70 as a tenth embodiment. The image conversion device 70 combined with the imaging device 80 with the AGC 81 includes a first look-up table 71-1, a second look instead of the first look-up table 13-1, the second look-up table 13-2 of the first embodiment. Up table 71-2, gain value G from AGC 81 and luminance Y of low frequency image, respectively 2 Coefficient C for two reference values 1 , C 2 Is stored and output. Coefficient C 1 , C 2 Any of the formulas 1, 2, 4, and 5 may be used, and the threshold value T in each equation is expressed as follows according to the gain value G of the AGC of the imaging device 80. 0 And T 1 To change. Threshold T 0 And T 1 Reduces the amount of suppression of low-frequency components in order to make noise less noticeable in low-luminance parts where the S / N ratio of the imaging device is small 2 This is a threshold value for determining the range. If this threshold value is increased too much, the range of the luminance in which the contrast of the high frequency component is reduced is widened. Therefore, the threshold T 0 And T 1 Is preferably the minimum size according to the noise level of the captured image. Generally, the darker the imaging target, the larger the AGC gain value G, and the higher the noise level of the brightness of the captured image. Conversely, the brighter the imaging target, the smaller the AGC gain value G, and the lower the noise level of the brightness of the captured image. So threshold T 0 And T 1 Is approximately proportional to the gain value G of AGC, noise can be suppressed at an appropriate luminance level at all times, preventing noise from being insufficiently suppressed or resulting in a conversion image that is unnecessarily blurred in the dark portion. can do.
[Brief description of the drawings]
FIG. 1 is a block diagram showing a configuration of an image conversion apparatus according to a first embodiment of the present invention.
FIG. 2 is a diagram illustrating an operation of a logarithmic converter.
FIG. 3 is a diagram illustrating an operation of a low-frequency image extractor.
FIG. 4 is a luminance Y of a low-frequency image in linear conversion of the image conversion apparatus according to the first embodiment of the present invention. 2 The graph which shows the characteristic of the two coefficients which are the functions of.
FIG. 5 is a graph showing frequency characteristics by linear transformation.
FIG. 6 is a block diagram showing a configuration of an image conversion apparatus according to a second embodiment of the present invention.
FIG. 7 is a diagram illustrating characteristics of a converted image input to a gradation converter according to the second embodiment.
FIG. 8 is a diagram illustrating an operation of a gradation converter according to the second embodiment.
FIG. 9 is a diagram illustrating an operation of a gradation converter according to the second embodiment.
FIG. 10 is a luminance Y of a low-frequency image in linear conversion of the image conversion apparatus according to the third embodiment of the present invention. 2 The graph which shows the characteristic of the two coefficients which are the functions of.
FIG. 11 is a luminance Y of a low-frequency image in linear conversion of the image conversion apparatus according to the fourth embodiment of the present invention. 2 The graph which shows the characteristic of the two coefficients which are the functions of.
12 is a luminance Y of a low-frequency image in linear conversion of the image conversion apparatus according to the fifth embodiment of the present invention. FIG. 2 The graph which shows the characteristic of the two coefficients which are the functions of.
FIG. 13 is a block diagram showing a configuration of an image conversion apparatus according to a sixth embodiment of the present invention.
FIG. 14 is a block diagram showing a configuration of an image conversion apparatus according to a seventh embodiment of the present invention.
FIG. 15 is a block diagram showing a configuration of an image conversion apparatus according to an eighth embodiment of the present invention.
FIG. 16 is a block diagram showing a configuration of an image conversion apparatus according to a ninth embodiment of the present invention.
FIG. 17 is a block diagram showing a configuration of an image conversion apparatus according to a tenth embodiment of the present invention.
FIG. 18 is a block diagram showing two types of configurations of Homomorphic Filtering.
[Explanation of symbols]
10, 20, 30, 40, 50, 60, 70 Image converter
11 Logarithmic converter
12 Low frequency image extractor
13-1, 71-1 first lookup table
13-2, 71-2 Second lookup table
14 Linear operator
15 gradation converter (look-up table)
16 Delay circuit
43, 51, 61 gradation converter
31 Low spatial frequency component suppression calculator (look-up table)
41, 42 Average calculator
80 Imaging device with AGC function
81 AGC

Claims (16)

各画素が輝度値を有する元画像を処理して、明るさダイナミックレンジを圧縮する画像変換方法において、
元画像の輝度の周波数の低い成分から成る低周波画像を抽出する低周波画像抽出工程と、
元画像の画素の輝度Y1に対応する低周波画像の画素の輝度Y2から、2つの係数C1、C2を求める係数決定工程と、
前記2つの係数C1、C2と、前記元画像の画素の輝度Y1と対応する前記低周波画像の画素の輝度Y2とから、対応する変換画像の画素の輝度Y3を式Y3=C1Y1+C2Y2により求める線形演算工程とを有し、
前記低周波画像の画素の輝度Y2の関数である前記2つの係数C1、C2はいずれも連続関数であり、
前記係数C1は正又は0であり、
前記係数C2は、輝度Y2の最小値を含む連続した第1の範囲で正、0及び負の値をとる広義の単調減少関数であり、輝度Y2の最大値を含む第2の連続した範囲で負又は0の値をとる広義の単調増加関数であることを特徴とする画像変換方法。
In an image conversion method in which an original image in which each pixel has a luminance value is processed and the brightness dynamic range is compressed,
A low-frequency image extraction step for extracting a low-frequency image composed of low-frequency components of the original image; and
A coefficient determination step for obtaining two coefficients C 1 and C 2 from the luminance Y 2 of the pixel of the low-frequency image corresponding to the luminance Y 1 of the pixel of the original image;
The two coefficients C 1, C 2, wherein from the low frequency of the pixels of the image luminance Y 2 Metropolitan and the corresponding luminance Y 1 of the pixel of the original image, the luminance Y 3 pixels of the corresponding converted image wherein Y 3 = C 1 Y 1 + C 2 Y 2
The two coefficients C 1 and C 2 that are functions of the luminance Y 2 of the pixels of the low-frequency image are both continuous functions,
The coefficient C 1 is positive or 0,
The coefficient C 2, the positive in the first contiguous range including a minimum value of the luminance Y 2, a monotonic decreasing function of broad taking 0 and negative values, the second continuous including the maximum value of the luminance Y 2 An image conversion method characterized by being a monotonically increasing function in a broad sense that takes a negative or zero value within the range.
前記2つの係数C1、C2の和が前記輝度Y2の全範囲のうちの少なくとも連続する第3の範囲において正の一定の値であることを特徴とする請求項1に記載の画像変換方法。2. The image conversion according to claim 1, wherein the sum of the two coefficients C 1 and C 2 is a positive constant value in at least a third continuous range of the entire range of the luminance Y 2. Method. 前記2つの係数C1、C2の和が前記輝度Y2の全範囲において一定であることを特徴とする請求項2に記載の画像変換方法。3. The image conversion method according to claim 2, wherein the sum of the two coefficients C 1 and C 2 is constant over the entire range of the luminance Y 2 . 前記輝度Y2の前記第3の範囲の内部である連続する第4の範囲において、前記係数C1が正の一定値であり、前記係数C2が負の一定値であることを特徴とする請求項2又は請求項3に記載の画像変換方法。The coefficient C 1 is a positive constant value and the coefficient C 2 is a negative constant value in a continuous fourth range that is within the third range of the luminance Y 2. The image conversion method according to claim 2 or claim 3. 前記輝度Y2の第1の範囲において前記輝度Y2の関数C2/C1が広義の単調減少関数であり、
前記輝度Y2の第2の範囲において前記輝度Y2の関数C2/C1が広義の単調増加関数であることを特徴とする請求項1乃至請求項4のいずれか1項に記載の画像変換方法。
Wherein the function C 2 / C 1 of the luminance Y 2 in the first range of luminance Y 2 is a monotonically decreasing function of broad,
Image according to any one of claims 1 to 4, wherein the function C 2 / C 1 of the luminance Y 2 in the second range of the luminance Y 2 is a monotonically increasing function of the broad Conversion method.
前記輝度Y2の前記第4の範囲及びそれを越える全ての範囲において、前記係数C1が正の一定値であることを特徴とする請求項5に記載の画像変換方法。In all ranges exceeding said fourth range and its the luminance Y 2, image conversion method according to claim 5, wherein the coefficient C 1 is a positive constant value. 前記2つの係数C1とC2は、前記輝度Y2の最小値を含む連続した第5の範囲においてC1が0であり、前記輝度Y2の最大値を含む連続した第6の範囲においてC2が0であり、当該2つの範囲に挟まれた範囲の連続した少なくとも一部において前記係数C1が正の一定値、前記係数C2が負の一定値であることを特徴とする請求項3に記載の画像変換方法。The two coefficients C 1 and C 2 are, C 1 is 0 in consecutive fifth range including the minimum value of the luminance Y 2, in successive sixth range including the maximum value of the luminance Y 2 C 2 is 0, and the coefficient C 1 is a positive constant value and the coefficient C 2 is a negative constant value in at least a part of a continuous range between the two ranges. Item 4. The image conversion method according to Item 3. 前記低周波画像の抽出が、各画素ラインごとの1次元方向のみにて行われることを特徴とする請求項1乃至請求項7のいずれか1項に記載の画像変換方法。The image conversion method according to claim 1, wherein the extraction of the low-frequency image is performed only in a one-dimensional direction for each pixel line. 各画素が輝度値を有する元画像を処理して、明るさダイナミックレンジを圧縮する画像変換装置において、
元画像の輝度の周波数の低い成分から成る低周波画像を抽出する低周波画像抽出手段と、
低周波画像抽出手段から出力された低周波画像の画素の輝度Y2に対応して、各々係数C1、C2を求める第1、第2の係数決定手段と、
当該第1、第2の係数決定手段により得られる係数C1、C2と、それに用いた前記低周波画像抽出手段から出力された低周波画像の画素の輝度Y2、及びそれに対応する元画像の画素の輝度Y1から、式Y3=C1Y1+C2Y2により変換画像の画素の輝度Y3を求める線形演算手段とを有し、
前記第1、第2の係数決定手段における前記係数C1、C2について、
前記低周波画像の画素の輝度Y2の関数である前記係数C1、C2はいずれも連続関数であり、
前記係数C1は正又は0であり、
前記係数C2は、輝度Y2の最小値を含む連続した第1の範囲で正、0及び負の値をとる広義の単調減少関数であり、輝度Y2の最大値を含む連続した第2の範囲で負又は0の値をとる広義の単調増加関数であることを特徴とする画像変換装置。
In an image conversion apparatus that processes an original image in which each pixel has a luminance value and compresses the brightness dynamic range,
Low-frequency image extraction means for extracting a low-frequency image composed of low frequency components of the luminance of the original image;
First and second coefficient determination means for obtaining coefficients C 1 and C 2 corresponding to the luminance Y 2 of the pixels of the low-frequency image output from the low-frequency image extraction means;
The coefficients C 1 and C 2 obtained by the first and second coefficient determining means, the luminance Y 2 of the pixel of the low frequency image output from the low frequency image extracting means used therefor, and the original image corresponding thereto Linear calculation means for calculating the luminance Y 3 of the pixel of the converted image from the luminance Y 1 of the pixel by the formula Y 3 = C 1 Y 1 + C 2 Y 2 ,
Regarding the coefficients C 1 and C 2 in the first and second coefficient determining means,
The coefficients C 1 and C 2 that are functions of the luminance Y 2 of the pixels of the low-frequency image are both continuous functions,
The coefficient C 1 is positive or 0,
The coefficient C 2, the positive in the first contiguous range including a minimum value of the luminance Y 2, a monotonic decreasing function of broad taking 0 and negative values, the second consecutive including the maximum value of the luminance Y 2 An image conversion apparatus characterized by being a monotonically increasing function in a broad sense that takes a negative or zero value within a range of.
前記係数C1、C2の和が前記輝度Y2の全範囲のうちの少なくとも連続する第3の範囲において正の一定の値であることを特徴とする請求項9に記載の画像変換装置The image conversion apparatus according to claim 9, wherein the sum of the coefficients C 1 and C 2 is a positive constant value in at least a third continuous range of the entire range of the luminance Y 2 . 前記係数C1、C2の和が前記輝度Y2の全範囲において一定であることを特徴とする請求項10に記載の画像変換装置。The image conversion apparatus according to claim 10, wherein the sum of the coefficients C 1 and C 2 is constant over the entire range of the luminance Y 2 . 前記輝度Y2の前記第3の範囲の内部である連続する第4の範囲において、前記係数C1が正の一定値であり、前記係数C2が負の一定値であることを特徴とする請求項10又は請求項11に記載の画像変換装置。The coefficient C 1 is a positive constant value and the coefficient C 2 is a negative constant value in a continuous fourth range that is within the third range of the luminance Y 2. The image conversion apparatus according to claim 10 or 11. 前記輝度Y2の第1の範囲において前記輝度Y2の関数C2/C1が広義の単調減少関数であり、
前記輝度Y2の第2の範囲において前記輝度Y2の関数C2/C1が広義の単調増加関数であることを特徴とする請求項9乃至請求項12のいずれか1項に記載の画像変換装置。
Wherein the function C 2 / C 1 of the luminance Y 2 in the first range of luminance Y 2 is a monotonically decreasing function of broad,
Image according to any one of claims 9 to 12, wherein the function C 2 / C 1 of the luminance Y 2 in the second range of the luminance Y 2 is a monotonically increasing function of the broad Conversion device.
前記輝度Y2の前記第4の範囲及びそれを越えた範囲の全てにおいて、前記係数C1が正の一定値であることを特徴とする請求項13に記載の画像変換装置。In all ranges beyond the fourth range and its the luminance Y 2, an image conversion apparatus according to claim 13, wherein the coefficient C 1 is a positive constant value. 前記係数C1とC2は、前記輝度Y2の最小値を含む連続した第5の範囲においてC1が0であり、前記輝度Y2の最大値を含む連続した第6の範囲においてC2が0であり、当該2つの範囲に挟まれた範囲の連続した少なくとも一部において前記係数C1が正の一定値、前記係数C2が負の一定値であることを特徴とする請求項11に記載の画像変換装置。It said coefficients C 1 and C 2, the luminance C 1 in consecutive fifth range including the minimum value of Y 2 is 0, C 2 in consecutive sixth range including the maximum value of the luminance Y 2 The coefficient C 1 is a positive constant value and the coefficient C 2 is a negative constant value in at least a part of a continuous range between the two ranges. The image conversion apparatus described in 1. 前記低周波画像抽出手段における低周波画像の抽出が、各画素ラインごとの1次元方向のみにて行われることを特徴とする請求項9乃至請求項15のいずれか1項に記載の画像変換装置。The image conversion apparatus according to any one of claims 9 to 15, wherein the low-frequency image is extracted by the low-frequency image extraction unit only in a one-dimensional direction for each pixel line. .
JP2001265895A 2001-09-03 2001-09-03 Image conversion method and image conversion apparatus Expired - Fee Related JP3644421B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2001265895A JP3644421B2 (en) 2001-09-03 2001-09-03 Image conversion method and image conversion apparatus
US10/232,692 US7065257B2 (en) 2001-09-03 2002-09-03 Image processing method and apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2001265895A JP3644421B2 (en) 2001-09-03 2001-09-03 Image conversion method and image conversion apparatus

Publications (2)

Publication Number Publication Date
JP2003076984A JP2003076984A (en) 2003-03-14
JP3644421B2 true JP3644421B2 (en) 2005-04-27

Family

ID=19092293

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2001265895A Expired - Fee Related JP3644421B2 (en) 2001-09-03 2001-09-03 Image conversion method and image conversion apparatus

Country Status (1)

Country Link
JP (1) JP3644421B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5932392B2 (en) * 2012-02-28 2016-06-08 キヤノン株式会社 Image processing apparatus and image processing method
US9418403B2 (en) 2013-10-04 2016-08-16 Fuji Xerox Co., Ltd. Image processor and non-transitory computer readable medium for generating a reproduction image which is reproduced so that visibility of an original image is enhanced

Also Published As

Publication number Publication date
JP2003076984A (en) 2003-03-14

Similar Documents

Publication Publication Date Title
JP3465226B2 (en) Image density conversion processing method
EP2833317B1 (en) Image display device and/or method therefor
US7020332B2 (en) Method and apparatus for enhancing a digital image by applying an inverse histogram-based pixel mapping function to pixels of the digital image
EP3087730B1 (en) Method for inverse tone mapping of an image
KR100756318B1 (en) Image Processing Method and System Using Gain-Controllable Clipped Histogram Equalization
KR100782845B1 (en) A digital image enhancement method and system using non-log domain illumination correction
KR100925419B1 (en) Color Image Enhancement using Laplacian Pyramid and method thereof
US7953286B2 (en) Automatic contrast enhancement
JP5157753B2 (en) Image processing apparatus, image processing method, and image processing program
US20140153820A1 (en) Image processing apparatus and method of processing image
JP6097588B2 (en) Image processing apparatus and image processing method
JP3545979B2 (en) Contour correction device
JP4021261B2 (en) Image processing device
JP2008511048A (en) Image processing method and computer software for image processing
JP2001014460A (en) Method for emphasizing edge contrast of digital picture
JPH10210324A (en) Image-processing unit
JP6600216B2 (en) Image processing apparatus, image processing method, program, and storage medium
Albu et al. One scan shadow compensation and visual enhancement of color images
JP3298330B2 (en) Imaging device
JP2003060913A (en) Image processing method and image processor
JP3644421B2 (en) Image conversion method and image conversion apparatus
Asari et al. Nonlinear enhancement of extremely high contrast images for visibility improvement
JP5295854B2 (en) Image processing apparatus and image processing program
JP4664938B2 (en) Image processing apparatus, image processing method, and program
KR100701306B1 (en) The device of processing the shoot artifact in the image signal and the method thereof

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20040914

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20040924

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: 20050111

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20050124

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090210

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100210

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110210

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110210

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120210

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130210

Year of fee payment: 8

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313532

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130210

Year of fee payment: 8

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130210

Year of fee payment: 8

LAPS Cancellation because of no payment of annual fees