(57)【要約】【課題】デジタル画像を低ダイナミック・レンジ（ＬＤＲ）画像から高ダイナミック・レンジ（ＨＤＲ）画像に変換する方法およびシステムを提供する。【解決手段】本発明は、画像を、所与のカメラを使って得た低ダイナミック・レンジ（ＬＤＲ）画像から高ダイナミック・レンジ（ＨＤＲ）画像に変換する方法を提供し、本方法は、  −  所与のカメラに対する露光量−ピクセル応答カーブを取得するステップ（２１）と、  −  ＬＤＲ画像をＨＳＢ色空間アレイに変換するステップ（２２）であって、前記ＨＳＢ色空間アレイは、色相アレイ、彩度アレイ、および明度アレイを含む、該変換するステップと、  −  露光量−ピクセル応答カーブの逆関数（ｆ−１）を使って、前記明度アレイ中の各ピクセルを逆マッピングすることによって放射輝度アレイを算定するステップ（２３、２４）と、を含む。【選択図】図１
2013507674.tif 000002
【特許請求の範囲】
【請求項１】
デジタル画像を、所与のカメラを使って得た低ダイナミック・レンジ（ＬＤＲ）画像から、高ダイナミック・レンジ（ＨＤＲ）画像に変換する方法であって、前記方法は、
−  前記所与のカメラに対する露光量−ピクセル応答カーブを取得するステップ（２１）と、
−  前記ＬＤＲ画像をＨＳＢ色空間アレイに変換するステップ（２２）であって、前記ＨＳＢ色空間アレイは、色相アレイ、彩度アレイ、および明度アレイを含む、前記変換するステップと、
−  前記露光量−ピクセル応答カーブの逆関数を使って、前記明度アレイ中の各ピクセルを逆マッピングすることによって放射輝度アレイを算定するステップ（２３、２４）と、
を含む方法。
【請求項２】
−  前記放射輝度アレイ中の各ピクセルに対し局部輝度平均を算定し（５０）、前記ピクセルに対する前記放射輝度アレイの値と前記局部輝度平均とから、各ピクセルに対する調整された放射輝度アレイを計算するステップ（５１）と、
−  前記放射輝度アレイを、前記ＨＳＢアレイの前記色相および彩度アレイと組み合わせるステップ（５２）と、
−  前記画像を、前記色相、彩度、および放射輝度アレイを使い画像色空間に変換し、前記ＨＤＲ画像を提供するステップ（５３）と、
をさらに含む、請求項１に記載の方法。
【請求項３】
前記露光量−ピクセル応答カーブはカメラのデータシートから取得される、先行する請求項のいずれかに記載の方法。
【請求項４】
前記露光量−ピクセル応答カーブは、同一の光景に対し異なる露光をされた一連の画像から取得される、請求項１または２に記載の方法。
【請求項５】
各ピクセル（ｕ，ｖ）に対し、前記調整された放射輝度Ｒａｄｉａｎｃｅ’［ｕ，ｖ］は次式により計算され、
Ｒａｄｉａｎｃｅ’［ｕ，ｖ］  ＝  Ｒａｄｉａｎｃｅ［ｕ，ｖ］＊Ｒａｄｉａｎｃｅ［ｕ，ｖ］／Ｌｏｃａｌ＿Ｌｕｍｉｎａｎｃｅ＿Ａｖｅｒａｇｅ［ｕ，ｖ］、
Ｒａｄｉａｎｃｅ［ｕ，ｖ］は、ピクセル（ｕ，ｖ）に対する前記放射輝度の値を指し、Ｌｏｃａｌ＿Ｌｕｍｉｎａｎｃｅ＿Ａｖｅｒａｇｅ［ｕ，ｖ］は、ピクセル（ｕ，ｖ）における前記局部輝度平均を指す、先行する請求項のいずれかに記載の方法。
【請求項６】
前記画像を画像色空間に変換する前記ステップは、次式により遂行され、
Ｃｏｎｖｅｒｔ（Ｈｕｅ［  ］，Ｓａｔｕｒａｔｉｏｎ［  ］，Ｒａｄｉａｎｃｅ［  ］／ｍａｘ（Ｒａｄｉａｎｃｅ［  ］））＊ｍａｘ（Ｒａｄｉａｎｃｅ）、
Ｈｕｅ［  ］は前記色相アレイを、Ｓａｔｕｒａｔｉｏｎ［  ］は前記彩度アレイを、Ｒａｄｉａｎｃｅ［  ］は前記放射輝度アレイを指す、先行する請求項のいずれかに記載の方法。
【請求項７】
前記局部輝度平均は次式で定義されるｉ畳み込みカーネルから算定され、
ｋｅｒｎｅｌｉ［  ］＝ＧａｕｓｓｉａｎＫｅｒｎｅｌ（ｒｉ）［  ］、
ｒｉはＧａｕｓｓｉａｎＫｅｒｎｅｌの半径であり、ｉは０〜８の間で変化する前記局部コントラストのスケール指標である、先行する請求項のいずれかに記載の方法。
【請求項８】
ｒｉは１／（２＊Ｓｑｒｔ（２））＊１．６ｉに設定される、請求項７に記載の方法。
【請求項９】
コンピュータ・プログラムが適切なコンピュータ・デバイス上で実行されたとき請求項１〜７のいずれか一つに記載の前記方法の前記ステップを実行するための命令を含む前記コンピュータ・プログラム。
【請求項１０】
中に符号化された、請求項９に記載のコンピュータ・プログラムのセットを有するコンピュータ可読媒体。
【請求項１１】
請求項１〜８のいずれか一つに記載の前記方法の前記ステップを実行するようになされた手段を含むシステム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、一般に画像処理に関し、さらに具体的には、デジタル画像を低ダイナミック・レンジ（ＬＤＲ）画像から高ダイナミック・レンジ（ＨＤＲ）画像に変換するための方法およびシステムに関する。
【背景技術】
【０００２】
デジタル写真術の分野における新興技術の一つに高ダイナミック・レンジ画像化（ＨＤＲＩ）がある。ＨＤＲＩは実世界の輝度のほとんどのキャプチャを提供し、適切なディスプレイを使用すれば、可能な限り現実に近い写真を再生することを可能にする。しかして、高ダイナミック・レンジ画像化は、実世界の光レベルに相応する値を有する光景の表現を提供する。実世界は、１２桁の範囲に及ぶ光度の広がりを生成し、これは、現在の多くのデジタル画像化の３桁よりもはるかに大きい。デジタル画像で各ピクセルが現在表現できる値の範囲は、通常、カラー・チャネル当たり２５６の値であり（最大６５５３６の値）、多くの光景を表現するのには適していない。ＨＤＲ画像を用いれば、光景を表現する光度の範囲で、何らかのディスプレイ・デバイスに合わせた範囲でなく人間の視力の限界に合った値の範囲によって、光景をキャプチャすることができる。現在のディスプレイ技術による表示に適した画像は、低ダイナミック・レンジ（ＬＤＲ）と呼ばれている。高ダイナミック・レンジ画像の表示品質は、在来の低ダイナミック・レンジ画像のそれよりもはるかに良好である。ＨＤＲ画像は、かかる画像のキャプチャ、保管、処理、および表示に関してＬＤＲ画像とは異なっており、写真術において急速に広い支持を得つつある。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】Ａ．Ｇ．Ｒｅｍｐｅｌ、Ｍ．Ｔｒｅｎｔａｃｏｓｔｅ、Ｈ．Ｓｅｅｔｚｅｎ、Ｈ．Ｄ．Ｙｏｕｎｇ、Ｗ．Ｈｅｉｄｒｉｃｈ、Ｌ．Ｗｈｉｔｅｈｅａｄ、およびＧ．Ｗａｒｄ著、論文題名「Ｌｄｒ２Ｈｄｒ：ｏｎ−ｔｈｅ−ｆｌｙ  ｒｅｖｅｒｓｅ  ｔｏｎｅ  ｍａｐｐｉｎｇ  ｏｆ  ｌｅｇａｃｙ  ｖｉｄｅｏ  ａｎｄ  ｐｈｏｔｏｇｒａｐｈｓ」、ＡＣＭ  ＳＩＧＧＲＡＰＨ  ２００７論文集（カリフォルニア州サンディエゴ、２００７年８月０５〜０９日）。
【非特許文献２】Ｆ．Ｂａｎｔｅｒｌｅ、Ｐ．Ｌｅｄｄａ、Ｋ．Ｄｅｂａｔｔｉｓｔａ、およびＡ．Ｃｈａｌｍｅｒｓ著、論文題名「Ｉｎｖｅｒｓｅ  ｔｏｎｅ  ｍａｐｐｉｎｇ」、第４回Ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｃｏｎｆｅｒｅｎｃｅ  ｏｎ  Ｃｏｍｐｕｔｅｒ  Ｇｒａｐｈｉｃｓ  ａｎｄ  ｉｎｔｅｒａｃｔｉｖｅ  Ｔｅｃｈｎｉｑｕｅｓ  ｉｎ  Ａｕｓｔｒａｌａｓｉａ  ａｎｄ  Ｓｏｕｔｈｅａｓｔ  Ａｓｉａ発表抄録（マーレシア、クアランプール、２００６年１１月２９日〜１２月０２日）、ＧＲＡＰＨＩＴＥ‘０６、ＡＣＭ、ＮＹ、ニューヨーク、３４９〜３５６頁。
【非特許文献３】Ｅｒｉｋ  Ｒｅｉｎｈａｒｄ、Ｍｉｃｈａｅｌ  Ｓｔａｒｋ、Ｐｅｔｅｒ  Ｓｈｉｒｌｅｙ、およびＪａｍｅｓ  Ｆｅｒｗｅｒｄａ著、論文題名「Ｐｈｏｔｏｇｒａｐｈｉｃ  ｔｏｎｅ  ｒｅｐｒｏｄｕｃｔｉｏｎ  ｆｏｒ  ｄｉｇｉｔａｌ  ｉｍａｇｅｓ」、ＡＣＭ  Ｔｒａｎｓ．Ｇｒａｐｈ．２１，３（２００２年７月）、２６７〜２７６頁。
【非特許文献４】Ｐ．Ｄｅｂｅｖｅｃ著、論文題名「Ａ  ｍｅｄｉａｎ  ｃｕｔ  ａｌｇｏｒｉｔｈｍ  ｆｏｒ  ｌｉｇｈｔ  ｐｒｏｂｅ  ｓａｍｐｌｉｎｇ」、ＡＣＭ  ＳＩＧＧＲＡＰＨ  ２００６  Ｃｏｕｒｓｅｓ（マサチューセッツ州ボストン、２００６年７月３０日〜８月０３日）、ＳＩＧＧＲＡＰＨ‘０６、ＡＣＭ、ＮＹ、ニューヨーク、６。
【非特許文献５】Ｔ．Ｏ．Ａｙｄｉｎ、Ｒ．Ｍａｎｔｉｕｋ、Ｋ．Ｍｙｓｚｋｏｗｓｋｉ、およびＨ．Ｓｅｉｄｅｌ著、論文題名「Ｄｙｎａｍｉｃ  ｒａｎｇｅ  ｉｎｄｅｐｅｎｄｅｎｔ  ｉｍａｇｅ  ｑｕａｌｉｔｙ  ａｓｓｅｓｓｍｅｎｔ」、ＡＣＭ  ＳＩＧＧＲＡＰＨ‘０８．（カリフォルニア州ロサンジェルス、２００８年８月１１〜１５日）論文集、ＡＣＭ．ＮＹ、ニューヨーク、１−１０。
【発明の概要】
【発明が解決しようとする課題】
【０００４】
デジタル写真術の分野でＨＤＲＩの使用が普及するにつれ、静止画像およびビデオの両方を表示することのできるＨＤＲＩディスプレイに対する必要性が増大している。これは、在来のディスプレイを超える表示品質の大幅なシフトを意味する。しかしながら、既存のメディアは高ダイナミック・レンジ（ＨＤＲ）のものでないので、ＨＤＲＩディスプレイの実用は、ＨＤＲＩセンサを使って新しく撮影されたＨＤＲ画像に限られている。既存の低ダイナミック・レンジ（ＬＤＲ）画像を相当するＨＤＲ画像に変換するための既存の対処技術は、「逆トーン・マッピング」として広く知られている。逆トーン・マッピングは、一般に２段階のフェーズを必要とする。第一フェーズでは、インプットＬＤＲ画像の輝度を、拡大ＨＤＲ輝度（ＨＤＲ放射輝度とも呼ばれる）に逆マップが行われる。このフェーズでは、画像量子化に起因して、高輝度の領域でディテールの喪失がもたらされ、ノイズが導入される。第二フェーズでは、ダイナミック・レンジをさらに増大する可能性も保ちながら、かかる領域を平滑化してこの欠陥を改善する。
【０００５】
この第一フェーズを実施するための一つの知られた対処技術に、“Ａ．Ｇ．Ｒｅｍｐｅｌ、Ｍ．Ｔｒｅｎｔａｃｏｓｔｅ、Ｈ．Ｓｅｅｔｚｅｎ、Ｈ．Ｄ．Ｙｏｕｎｇ、Ｗ．Ｈｅｉｄｒｉｃｈ、Ｌ．Ｗｈｉｔｅｈｅａｄ、およびＧ．Ｗａｒｄ著、論文題名「Ｌｄｒ２Ｈｄｒ：ｏｎ−ｔｈｅ−ｆｌｙ  ｒｅｖｅｒｓｅ  ｔｏｎｅ  ｍａｐｐｉｎｇ  ｏｆ  ｌｅｇａｃｙ  ｖｉｄｅｏ  ａｎｄ  ｐｈｏｔｏｇｒａｐｈｓ」、ＡＣＭ  ＳＩＧＧＲＡＰＨ  ２００７論文集（カリフォルニア州サンディエゴ、２００７年８月０５〜０９日）”で取られたアプローチがある。このアプローチは、リアルタイムのビデオ処理に適した高速の反転法に依存している。このアプローチによれば、逆ガンマ・マッピングが行われ、次いで、ダイナミック・レンジが５０００に拡大される。さらに、量子化の影響を低減するために平滑化フィルタリングが行われる。
【０００６】
逆トーン・マッピングの第一フェーズを実施する別の対処技術が“Ｆ．Ｂａｎｔｅｒｌｅ、Ｐ．Ｌｅｄｄａ、Ｋ．Ｄｅｂａｔｔｉｓｔａ、およびＡ．Ｃｈａｌｍｅｒｓ著、論文題名「Ｉｎｖｅｒｓｅ  ｔｏｎｅ  ｍａｐｐｉｎｇ」、第４回Ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｃｏｎｆｅｒｅｎｃｅ  ｏｎ  Ｃｏｍｐｕｔｅｒ  Ｇｒａｐｈｉｃｓ  ａｎｄ  ｉｎｔｅｒａｃｔｉｖｅ  Ｔｅｃｈｎｉｑｕｅｓ  ｉｎ  Ａｕｓｔｒａｌａｓｉａ  ａｎｄ  Ｓｏｕｔｈｅａｓｔ  Ａｓｉａ発表抄録（マーレシア、クアランプール、２００６年１１月２９日〜１２月０２日）、ＧＲＡＰＨＩＴＥ‘０６、ＡＣＭ、ＮＹ、ニューヨーク、３４９〜３５６頁”に記載されている。この対処技術は、以前に“Ｅ．Ｒｅｉｎｈａｒｄ、Ｍ．Ｓｔａｒｋ、Ｐ．Ｓｈｉｒｌｅｙ、およびＪ．Ｆｅｒｗｅｒｄａ著、論文題名「Ｐｈｏｔｏｇｒａｐｈｉｃ  ｔｏｎｅ  ｒｅｐｒｏｄｕｃｔｉｏｎ  ｆｏｒ  ｄｉｇｉｔａｌ  ｉｍａｇｅｓ」、ＡＣＭ  Ｔｒａｎｓ．Ｇｒａｐｈ．２１，３（２００２年７月）、２６７〜２７６頁”に記載されたグローバル・トーン・マッピング演算子に基づいた逆マッピング関数を使用する。次いで、二次元方程式を解いて逆数が得られ、これにより、相当により大きなダイナミック・レンジが生成され、特定のピクセルでは選択的にそのレンジが縮小される。しかしながら、これら既存の対処技術は、第一フェーズに対し、精度が十分でない逆トーン・マッピング関数を用いている。得られた放射輝度は、「包括的な」逆マッピング関数に起因して実世界の放射輝度とは正確に整合しない。これは実世界の輝度値を大まかに近似する。
【０００７】
逆トーン・マッピングの第二フェーズを実行するための２つの異なるアプローチがある。第一のアプローチは、Ｒａｍｐｅｌら著、論文題名「Ｌｄｒ２Ｈｄｒ：ｏｎ−ｔｈｅ−ｆｌｙ  ｒｅｖｅｒｓｅ  ｔｏｎｅ  ｍａｐｐｉｎｇ  ｏｆ  ｌｅｇａｃｙ  ｖｉｄｅｏ  ａｎｄ  ｐｈｏｔｏｇｒａｐｈｓ」、ＡＣＭ  ＳＩＧＧＲＡＰＨ  ２００７論文集（カリフォルニア州サンディエゴ、２００７年８月０５〜０９日）”に記載しているもので、ある高値を超えるピクセルを覆うガウス・マスクを生成する。加えて、このアプローチは、エッジにおける局部コントラストを改善するために「エッジ停止」関数を使用する。この合成明度関数を使って明暗度が大幅に拡大される。さらに複雑な技法は、“Ｆ．Ｂａｎｔｅｒｌｅら著、論文題名「Ｉｎｖｅｒｓｅ  ｔｏｎｅ  ｍａｐｐｉｎｇ」、第４回Ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｃｏｎｆｅｒｅｎｃｅ  ｏｎ  Ｃｏｍｐｕｔｅｒ  Ｇｒａｐｈｉｃｓ  ａｎｄ  ｉｎｔｅｒａｃｔｉｖｅ  Ｔｅｃｈｎｉｑｕｅｓ  ｉｎ  Ａｕｓｔｒａｌａｓｉａ  ａｎｄ  Ｓｏｕｔｈｅａｓｔ  Ａｓｉａ発表抄録（マーレシア、クアランプール、２００６年１１月２９日〜１２月０２日）、ＧＲＡＰＨＩＴＥ‘０６、ＡＣＭ、ＮＹ、ニューヨーク、３４９〜３５６頁”に記載されたものである。この第二のアプローチは、中央値切捨てアルゴリズム（Ｐ．Ｄｅｂｅｖｅｃ著、論文題名「Ａ  ｍｅｄｉａｎ  ｃｕｔ  ａｌｇｏｒｉｔｈｍ  ｆｏｒ  ｌｉｇｈｔ  ｐｒｏｂｅ  ｓａｍｐｌｉｎｇ」、ＡＣＭ  ＳＩＧＧＲＡＰＨ  ２００６  Ｃｏｕｒｓｅｓ（マサチューセッツ州ボストン、２００６年７月３０日〜８月０３日）、ＳＩＧＧＲＡＰＨ‘０６、ＡＣＭ、ＮＹ、ニューヨーク、６）を使い、インプット画像を、等しい光度の領域によって区分けすることを含んでいる。これら領域の中心部を使って、光密度が見積もられ「拡大」マップが構築される。次いで、このマップを使い、インプットＬＤＲと逆マップされたＬＤＲ画像との間で補間演算を誘導することによって最終的なＨＤＲ画像が生成される。逆トーン・マッピングの第二フェーズに対するこれらの対処技術は、高輝度値を有するピクセルを見出すことと、それを使ってそれらのピクセルのダイナミック・レンジを拡大することに依存している。しかしながら、かかる補間は、ホットスポット（ハイライト部）および近辺領域の輝度を拡大するためだけに行われ、暗い領域（陰影部）の輝度が低減されることはない。従って、これら技術は、局部的演算を使って一方だけのダイナミック・レンジの拡大を効果的に行い（陰影部は包括的に拡大される）、それが得られたＨＤＲ画像中の陰影領域の品質に影響を及ぼしている。
【課題を解決するための手段】
【０００８】
これらのおよび他の問題に対処するために、添付の独立請求項１によって、方法が、添付の請求項９、１０、および１１によって、それぞれコンピュータ・プログラム、コンピュータ可読媒体、およびシステムが提供される。好適な実施形態が、添付の従属請求項に定義されている。
【０００９】
概して言えば、本発明は、ＬＤＲ画像をＨＤＲ画像に変換するための改良された逆トーン・マッピングを提供する。
【００１０】
本発明のある態様によれば、従来の対処技術で提供されるような逆ガンマまたは標準の固定的な逆関数を実行する代わりに、カメラの応答カーブを使って画像の放射輝度マップを再構築する。
【００１１】
本発明の別の態様は、画像の輝度を選択的に、増大または低減するためそれぞれ覆い焼きおよび焼き込み操作に依存する。本発明は、ＨＤＲドメインにおいて覆い焼きおよび焼き込み様の操作を適用し、画像のダイナミック・レンジを拡大する。また、この態様は一般的に局部的コントラストを拡大し、ＬＤＲ画像では見えなかったより多くのディテールを視認可能にする。さらに、ＨＤＲドメインにおいて覆い焼きおよび焼き込み様の操作を適用することによって平滑化の実施が可能になり、これにより量子化の影響が低減される。
【００１２】
本発明の利点の一部には、以下が含まれる。
１−  単一のＬＤＲ画像からＨＤＲ値の再生を行うための単純化された技法。この技法は、従来技術のように一連の異なる露光の画像から画像を再構築するよりも容易である。
２−  輝度チャネルの２つのテール（高輝度および低輝度部分）の両方からのダイナミック・レンジの拡大。従来の対処技術は高輝度部分を拡大する。
３−  画像の中間／陰影部分に亘ってディテール可視性レベルの向上。
４−  カメラ応答カーブを使った、画像のより現実に近い放射輝度マップ。
５−  実証済みの写真技法を使った、仕上がり画像の向上された観賞品質。
６−  医療画像など重要な用途に対し、既存のＬＤＲ画像を見るために新規のＨＤＲディスプレイを用いることができる。
７−  逆トーン・マッピングを使い、高品質のＨＤＲ画像を生成し、順方向にトーン・マップしてＬＤＲに戻すことによる、可視コントラストの大幅な改善が実現された、ＬＤＲカメラ・アウトプットの品質の向上。
８−  さらなる画像処理作業のための、より容易なエッジの検出。
９−  医療画像化分野におけるエッジ検出作業などさらなる画像処理でのＨＤＲドメインにおける作業の改良など、デジタル画像に対する向上されたオペレーションの適用。
【００１３】
当業者は、図面および詳細説明を検討することにより本発明のさらなる利点を明らかにできよう。本明細書は、一切の追加的な利点も組み込むように意図されている。
【００１４】
以下に、例示として添付の図面を参照しながら本発明の実施形態を説明するものとし、図面中の同じ参照記号は同様なエレメントを表す。
【００１５】
これら図面は、本発明の典型的な実施形態だけを示すよう意図されており、従って、本発明の範囲を限定するものと見なすべきではない。
【図面の簡単な説明】
【００１６】
【図１】本発明の実施形態による、低ダイナミック・レンジ（ＬＤＲ）画像から高ダイナミック・レンジ（ＨＤＲ）画像を生成するためのシステム・ブロック図である。
【図２】本発明の実施形態による、低ダイナミック・レンジ（ＬＤＲ）画像からの高ダイナミック・レンジ（ＨＤＲ）画像の生成を表すハイレベル・フローチャートである。
【図３】例示的カメラ応答カーブを示す。
【図４】本発明の実施形態による、図３のカメラ応答カーブから得られた逆応答カーブを示す。
【図５】ＨＤＲドメインにおける覆い焼きおよび焼き込み操作を表すフローチャートを示す。
【図６】本発明の実施形態による、逆トーン・マッピングの第一フェーズに対する初期放射輝度アレイの例示的ヒストグラムを示す。
【図７】本発明の実施形態による、逆トーン・マッピングの第二フェーズに対する初期輝度アレイの例示的ヒストグラムを示す。
【図８】本発明の実施形態による、画像生成を表すフローチャートである。
【図９】画像を比較することによって得られた、青／緑ピクセルの和の例示的値を示す表を示す。
【図１０】図８の表中に包含された値を表現する図である。
【図１１】本発明の実施形態による、ＬＤＲ画像をＨＤＲ画像に変換するため用いられるあるコンピュータ・システムを示す。
【発明を実施するための形態】
【００１７】
本発明は、デジタル画像を低ダイナミック・レンジ（ＬＤＲ）画像から高ダイナミック・レンジ（ＨＤＲ）画像に変換するための、逆トーン・マッピングの方法およびシステムを提供する。本発明による逆トーン・マッピング対処技術は、デジタル計算および関連ロジックを実行するため使用可能な任意の種類のプロセッサに適応する。
【００１８】
図１は、インプットＬＤＲ画像からＨＤＲ画像を生成するためのシステム１００を示す。図１は、システム・コンポーネント間のデータ・フローを示し、各コンポーネントは、文字「Ｃ」で始まり後に番号が続く称号で標識されている。これらシステム・コンポーネントは、以降、Ｃ１、Ｃ２、…、Ｃ５の標記を使って呼ばれる。このシステムは、まず、ＬＤＲ画像を入力し、それを色空間コンバータＣ１に送り込むことから作動され、画像は、色相、彩度、および明度の色空間（ＨＳＢ：Ｈｕｅ、Ｓａｔｕｒａｔｉｏｎ、Ｂｒｉｇｈｔｎｅｓｓ）に変換され、３つの対応するアレイが生成される。明度アレイおよびカメラ・センサ応答カーブが逆マッパーＣ２に供給される。次いで該マッパーは供給された応答カーブを使って明度アレイを放射輝度アレイに逆マップする。局部コントラスト・エンジンＣ３は、放射輝度アレイに対する局部的コントラストを計算し、局部輝度平均（ＬＬＡ：Ｌｏｃａｌ  Ｌｕｍｉｎａｎｃｅ  Ａｖｅｒａｇｅ）アレイを生成する。覆い焼き／焼き込みエンジンＣ４は、ＬＬＡアレイを使って、放射輝度アレイに覆い焼き／焼き込みを行い、拡大放射輝度アレイを生成する。
【００１９】
生成されたアレイは、彩度および色相アレイとともに、色空間コンバータＣ５に送り込まれ、画像は元のインプットＬＤＲの色空間（例えばＲＧＢ）に変換され、アウトプットＨＤＲ画像が生成される。
【００２０】
ＬＤＲ画像に対し逆トーン・マッピングを実施するため、第一フェーズの過程で、インプットＬＤＲ画像の輝度は拡大ＨＤＲ輝度に逆マップされる。この初期ステップは、高輝度の領域でディテールの喪失およびノイズの導入をもたらし、これは、これらの領域を平滑化する第二フェーズの過程で改善される。
【００２１】
本発明の実施形態によれば、逆トーン・マッピングの第一フェーズは、キャプチャしたセンサ／デバイスの情報を使って実施される。フィルム・カメラと同様、デジタル・カメラを使ってキャプチャされたほとんどの画像に対し、そのセンサの識別は容易に入手可能である。この第一フェーズでは、カメラ応答カーブを使い、インプット画像の逆トーン・マッピングによって初期放射輝度マップが構築される。
【００２２】
図２は、本発明の実施形態による、この、逆トーン・マッピングの第一フェーズを記述したハイレベルのフローチャートを示す。
【００２３】
ステップ２０で、最初にＬＤＲ画像が入力され、次いでそのＬＤＲ画像がメモリ中に格納される。かかるＬＤＲ画像は、さまざまな色空間で表現されることができよう。標準的な色空間はＲＧＢ空間であり、この色空間では、各ピクセルの色は、３つの成分、すなわち赤（Ｒ：Ｒｅｄ）、緑（Ｇ：Ｇｒｅｅｎ）、および青（Ｂ：Ｂｌｕｅ）で表現される。別の色空間はＬｕｖであり、このＬは輝度成分、ｕおよびｖは各々クロミナンス成分である。逆トーン・マッピングは、輝度チャネル上で実行される。
【００２４】
ステップ２１において、所与のカメラに対する露光量−ピクセル応答カーブが取得される。露光量−ピクセル応答カーブは各カメラに内在的なもので、センサ・データシートから直接得ることができよう。あるいは、同一光景に対する一連の異なる露光量の画像を解析することによって、露光量−ピクセル応答カーブを算定することもできる。露光量応答カーブは、測定し、計算し、または見積もることができ、あるいは遠隔サイトから入手することすら可能である。このカーブは、カメラごとにほぼ一定であり、しかして、同じカメラによって生成されたさらなる画像に対し応答関数の再構築を繰り返す必要はない。以下の説明において、露光量−ピクセル応答カーブを関数「ｆ（ｘ）」で表すこととし、このｘは所与の露光量値を示す。このとき、ｙ＝ｆ（ｘ）は、所与の露光量ｘに対するピクセル輝度値を表すものとする。露光量−ピクセル応答カーブは、対応写真に対する実世界の放射輝度値を提供する。これらの値は、従来技術の単純な逆ガンマ手法を使って得られるものよりも正確である。
【００２５】
ステップ２２において、ＬＤＲインプット画像はＨＳＢ色空間アレイに変換され、このＨは色相を表し、Ｓは彩度を表し、Ｂは明度を表す。ＨＳＢ（ＨＳＬ表現としても知られる）は、知覚的色関係をＲＧＢよりも正確に表現する、ＲＧＢ色モデルにおけるポイント群の表現である以下の説明では、色相アレイを「Ｈｕｅ［  ］」で、彩度アレイを「Ｓａｔｕｒａｔｉｏｎ［  ］」で、明度アレイを「Ｂｒｉｇｈｔｎｅｓｓ［  ］」で表すことにする。
【００２６】
ステップ２３において、明度アレイ（Ｂｒｉｇｈｔｎｅｓｓ［  ］）の各ピクセルに対し、応答カーブを使ってその値が逆マップされ、得られた露光量がＲａｄｉａｎｃｅ（放射輝度）［ｉ］に格納される。露光量―ピクセル応答関数ｆ（ｘ）の逆関数を表すのに「ｆ−１（ｙ）」の表記を使えば、このとき、インプットピクセル値ｙに対し、ｆ−１（ｙ）は、露光量値ｘを戻してくる。上記から、チャネル反復子ｉに対する放射輝度は次式により定義される。
Ｒａｄｉａｎｃｅ［ｉ］＝ｆ−１（Ｂｒｉｇｈｔｎｅｓｓ［ｉ］）、
このＲａｄｉａｎｃｅ［  ］アレイが得られた放射輝度マップである。
【００２７】
斯く得られた放射輝度マップは、ステップ２４においてメモリに格納される。本発明は、カメラの露光量−ピクセル応答カーブを用いて画像の放射輝度マップを再構築することにより、逆ガンマまたは標準の固定的な逆関数に依存する従来の対処技術によって得られるものよりも現実に近い放射輝度マップを提供する。
【００２８】
図３は、例示的カメラの応答カーブを示す。「露光量」と題されたｘ−軸は実世界の露光量値を表し、「ピクセル明度値」と題されたｙ−軸はカメラによって記録された対応するピクセル放射輝度値を表す。「露光量」は、センサが受光した放射照度に露光時間を乗じた量として定義される。単位は平方メートル当たりワット秒である。ピクセル放射輝度の値は、０〜２５５の整数を取る一方、露光量値は実数である。
【００２９】
図４は、図３に表された例示的カメラ応答カーブに対する逆応答カーブを示す。「ピクセル明度値」と題されたｘ−軸はピクセル放射輝度値を指し、「露光量」と題されたｙ−軸は対応する実世界の露光量値を指す。
【００３０】
この応答カーブは、現実をできるだけ近似して模擬する。この段階で、ステップ２４で得られた放射輝度アレイ「Ｒａｄｉａｎｃｅ［  ］」を、色相アレイ「Ｈｕｅ［  ］」および彩度アレイ「Ｓａｔｕｒａｔｉｏｎ［  ］」と組み合わせ、これらを一つの画像に関連付け、次いでその画像を元の画像の色空間に変換してＨＤＲ画像を提供することも可能ではあろう。しかしながら、カメラ・カーブが、それ自体で十分な高ダイナミック・レンジを生成できないことがあり、量子化によるアーチファクトを取り込む可能性がある。本発明による逆トーン・マッピングの第二フェーズは、カメラ応答カーブの使用により導入される可能性のあるマッピング円滑性の不十分さおよびアーチファクトを補償する。
【００３１】
図５は、本発明の実施形態による、逆トーン・マッピングの第二フェーズを図解する。この第二フェーズの過程では、画像の輝度を選択的に増大または低減させるため、それぞれ覆い焼きおよび焼き込み様の操作が行われる。「覆い焼き」操作は輝度を増大させ、一方「焼き込み」操作はピクセルの輝度を低減させる。本発明は、ＨＤＲドメインにおいて覆い焼きおよび焼き込み様の操作を適用し、画像のダイナミック・レンジを拡大する。ハイライト部に対し、覆い焼き操作は、低幅な局部コントラストの領域のダイナミック・レンジをそれ以上に拡大することになろう。陰影部に対しては、焼き込み操作が左側のダイナミック・カーブをもっと拡大し、ピクセルの最小輝度を低下させることになる（これにより画像のダイナミック・レンジはさらに拡大される）。また、これは全般に局部コントラストを拡大し、これにより、ＬＤＲ画像では見えなかったより多くのディテールが視認可能になる。さらに、ＨＤＲドメインに覆い焼きおよび焼き込み様の操作を適用することによって、平滑化が遂行でき、これにより、量子化の影響が低減される。
【００３２】
ステップ５０において、Ｒａｄｉａｎｃｅ［  ］アレイ中の各ピクセル（ｕ，ｖ）に対し、局部輝度平均、Ｌｏｃａｌ＿Ｌｕｍｉｎａｎｃｅ＿Ａｖｅｒａｇｅ［ｕ，ｖ］が計算される。ピクセル（ｕ，ｖ）におけるＬｏｃａｌ＿Ｌｕｍｉｎａｎｃｅ＿Ａｖｅｒａｇｅは、“Ｅｒｉｋ  Ｒｅｉｎｈａｒｄ、Ｍｉｃｈａｅｌ  Ｓｔａｒｋ、Ｐｅｔｅｒ  Ｓｈｉｒｌｅｙ、およびＪａｍｅｓ  Ｆｅｒｗｅｒｄａ著、論文題名「Ｐｈｏｔｏｇｒａｐｈｉｃ  ｔｏｎｅ  ｒｅｐｒｏｄｕｃｔｉｏｎ  ｆｏｒ  ｄｉｇｉｔａｌ  ｉｍａｇｅｓ」、ＡＣＭ  Ｔｒａｎｓ．Ｇｒａｐｈ．２１，３（２００２年７月）、２６７〜２７６頁”において展開されているアプローチから、以下のように計算することができる。
−  畳み込みカーネルが、ｋｅｒｎｅｌｉ［  ］＝ＧａｕｓｓｉａｎＫｅｒｎｅｌ（ガウス・カーネル）（ｒｉ）［  ］として設定され、このｒｉはＧａｕｓｓｉａｎＫｅｒｎｅｌの半径であり、ｉは局部コントラストのスケール指標を指す。ｒｉの値は変化する。本発明のある特定の実施形態において、この値は１／（２＊Ｓｑｒｔ（２））＊１．６ｉに設定される。ｉの値は、０〜８の範囲で変化する。
−  値ｉに対するピクセル（ｕ，ｖ）における局部輝度平均は、Ｌｏｃａｌ＿Ｌｕｍｉｎａｎｃｅ＿Ａｖｅｒａｇｅｉ［ｕ，ｖ］＝ｋｅｒｎｅｌｉ  Ｒａｄｉａｎｃｅ［  ］として計算される。
−  次いで、パラメータｉの最小値ｍが、Ａｂｓ（Ｌｏｃａｌ＿Ｌｕｍｉｎａｎｃｅ＿Ａｖｅｒａｇｅｉ［ｕ，ｖ］−Ｌｏｃａｌ＿Ｌｕｍｉｎａｎｃｅ＿Ａｖｅｒａｇｅｉ＋１［ｕ，ｖ］）＜εにより計算され、このεは閾値を指し、ｉの値は０〜７の範囲で変化する。
−  Ｌｏｃａｌ＿Ｌｕｍｉｎａｎｃｅ＿Ａｖｅｒａｇｅ［ｕ，ｖ］は最終的にＬｏｃａｌＬｕｍｉｎａｎｃｅＡｖｅｒａｇｅｍ［ｕ，ｖ］に設定され、これによりピクセル（ｕ，ｖ）における局部輝度平均が与えられる。
【００３３】
ステップ５１において、各ピクセルの輝度は、局部輝度平均Ｌｏｃａｌ＿Ｌｕｍｉｎａｎｃｅ＿Ａｖｅｒａｇｅ［ｕ，ｖ］を使って調整される。新規の放射輝度値「Ｒａｄｉａｎｃｅ’［ｕ，ｖ］」は次式により定義される。
Ｒａｄｉａｎｃｅ’［ｕ，ｖ］  ＝  Ｒａｄｉａｎｃｅ［ｕ，ｖ］＊Ｒａｄｉａｎｃｅ［ｕ，ｖ］／Ｌｏｃａｌ＿Ｌｕｍｉｎａｎｃｅ＿Ａｖｅｒａｇｅ［ｕ，ｖ］
【００３４】
本発明によれば、この演算は、ＨＤＲドメインにおいて行われる。実際に、本出願人は、ピクセルの周囲（ｕ，ｖ）が原ピクセル（ｕ，ｖ）よりも明度が高い場合、ピクセル（ｕ，ｖ）の放射輝度は低減され、これにより局部コントラストが増大されるのを観察した。同様に、周囲のピクセルの明度が原ピクセルより低い場合、原ピクセル（ｕ，ｖ）の放射輝度が増大され、これにより局部コントラストが増大されることになるのも観察されている。この放射輝度の低減および増大の双方は、Ｒａｄｉａｎｃｅ［ｕ，ｖ］／Ｌｏｃａｌ＿Ｌｕｍｉｎａｎｃｅ＿Ａｖｅｒａｇｅ［ｕ，ｖ］の比率によって設定される。この比率は、放射輝度のスケーリング係数の機能を果たし、圧縮調整なしの任意のスケーリングを可能にする。
【００３５】
本発明のある実施形態による逆トーン・マッピングの第二フェーズは、従って、覆い焼きおよび焼き込みを介して画像の局部コントラストを高めながら、放射輝度のダイナミック・レンジを増大する。このようにして、本発明は覆い焼き／焼き込みの写真術コンセプトを活用して写真品質のＨＤＲ画像を生成する。
【００３６】
ステップ５２において、放射輝度アレイ「Ｒａｄｉａｎｃｅ［  ］」は、色相アレイ「Ｈｕｅ［  ］」および彩度アレイ「Ｓａｔｕｒａｔｉｏｎ［  ］」と組み合わされ、次いでこれらは一つの画像に関連付けられ、次いでステップ５３でこの画像は元画像の色空間に変換される。（例えばｃｉｍｇ．ｓｏｕｒｃｅｆｏｒｇｅ．ｏｒｇのｃｉｍｇライブラリ中のｇｅｔ＿ＲＧＢｔｏＨＳＶ（  ）などの）標準的ライブラリ・ルーチンを用い、この画像変換は次式を使って計算され、このＣｏｎｖｅｒｔ（  ）が画像を元画像色の空間に逆変換する。
Ｃｏｎｖｅｒｔ（Ｈｕｅ［  ］，Ｓａｔｕｒａｔｉｏｎ［  ］，Ｒａｄｉａｎｃｅ［  ］／ｍａｘ（Ｒａｄｉａｎｃｅ［  ］））＊ｍａｘ（Ｒａｄｉａｎｃｅ）。
【００３７】
通常ピクセル値は０から１の間で変化するように正規化されているので、この式の除算部「Ｒａｄｉａｎｃｅ［  ］／ｍａｘ（Ｒａｄｉａｎｃｅ［  ］）」が使われている。
【００３８】
図６は、本発明の実施形態による、初期放射輝度アレイの例示的ヒストグラムを示し、これは逆トーン・マッピングの第一フェーズだけを使って得られたものである。この放射輝度アレイはインプットＬＤＲ画像から生成される。このヒストグラムのｘ−軸はｌｏｇ２露光量値を表し、ｙ−軸は発生頻度を表す。図示のように、最大ｌｏｇ２露光量は８．３であり最小は−９．３８で、これは２１０，３８１（５桁）のダイナミック・レンジとなる。これは実世界の通常の光景における典型的な値である。
【００３９】
図７は、本発明の実施形態による、逆トーン・マッピングの第二フェーズ後のヒストグラムを示す。最大のｌｏｇ２露光量は、これでは８．７１であり最小は−１５．５２で、このダイナミック・レンジは１．９８×１０７である。従って、第二フェーズによってダイナミック・レンジは２桁増大したことになる。また、このヒストグラムがより平滑で、陰影部を強調しながら両側に広く広がっていることは留意に値する。
【００４０】
生成されたＨＤＲ画像の品質を評定するために、ＭＰＩ  ＨＤＲ（ＭＰＩはＭａｘ  Ｐｌａｎｃｋ  Ｉｎｓｔｉｔｕｔｅ  ｉｎｆｏｒｍａｔｉｋ（マックス・プランク境界情報科学）の頭字語）メトリックを使うことができる。このメトリックは、“Ｔ．Ｏ．Ａｙｄｉｎ、Ｒ．Ｍａｎｔｉｕｋ、Ｋ．Ｍｙｓｚｋｏｗｓｋｉ、およびＨ．Ｓｅｉｄｅｌ著、の論文題名「Ｄｙｎａｍｉｃ  ｒａｎｇｅ  ｉｎｄｅｐｅｎｄｅｎｔ  ｉｍａｇｅ  ｑｕａｌｉｔｙ  ａｓｓｅｓｓｍｅｎｔ」、ＡＣＭ  ＳＩＧＧＲＡＰＨ  ２００８年論文（カリフォルニア州ロサンジェルス、２００８年８月１１〜１５日）ＳＩＧＧＲＡＰＨ‘０８．ＡＣＭ，ＮＹ、ニューヨーク、１−１０”に定義されている。この画像品質メトリックは、両方の画像が任意のダイナミック・レンジを有する画像のペアに対して使われる。このメトリックによれば、青、緑、および赤のピクセルを使って要約画像が生成される。最も高いコントリビュータに応じて各ピクセルに対する色が決定される。青ピクセルは、向上されたコントラストを有するピクセルを表し（インプット画像では見えずアウトプット画像で見える）、緑ピクセルはコントラストの喪失を表し、赤ピクセルはコントラストの逆転を表す。各色の値は、その影響の確率が見えるように表される。本発明人は、メトリックの青および緑パラメータが、画像品質を評定する上で特に重要に見えることに気付いた。実際に、できるだけ緑を低減し青を増加させると、コントラストが向上でき、可視のディテールが喪失されないように見えた。結果として、各ピクセルに対する正規化された優位色の和も観測することによって、画像品質評定の改良が達成された。
【００４１】
図８は、いくつかの画像の比較に基づいて、生成されたＨＤＲ画像の品質を評定するために行われるステップを表したフローチャートを示す。後記の実験の理解を容易にするため、後記で、画像の識別に使われる特定の記号の定義を以下に示す。
１−  「Ｒｅａｌ−ＨＤＲ」は、露光を変えた画像の大きなシーケンスを使って得られた画像を指す。言い換えれば「現実の」ＨＤＲ画像である。
２−  「Ｉｎｐｕｔ−ＬＤＲ」はインプットＬＤＲ画像を指す。
３−  「Ｇｅｎ−ＨＤＲ」は、Ｉｎｐｕｔ−ＬＤＲ画像を逆処理して得られた、生成ＨＤＲ画像を指す。
４−  「Ｇｅｎ−ＬＤＲ」は、Ｒｅｉｎｈａｒｄらの写真トーン・マッピング演算子（Ｅｒｉｋ  Ｒｅｉｎｈａｒｄ、Ｍｉｃｈａｅｌ  Ｓｔａｒｋ、Ｐｅｔｅｒ  Ｓｈｉｒｌｅｙ、およびＪａｍｅｓ  Ｆｅｒｗｅｒｄａ著、論文題名「Ｐｈｏｔｏｇｒａｐｈｉｃ  ｔｏｎｅ  ｒｅｐｒｏｄｕｃｔｉｏｎ  ｆｏｒ  ｄｉｇｉｔａｌ  ｉｍａｇｅｓ」、ＳＩＧＧＲＡＰＨ’０２：Ｃｏｍｐｕｔｅｒ  Ｇｒａｐｈｉｃｓ  ａｎｄ  Ｉｎｔｅｒａｃｔｉｖｅ  Ｔｅｃｈｎｉｑｕｅｓ第２９回年度協議会抄録、２６７〜２７６頁、米国ＮＹ州ニューヨーク、２００２年。ＡＣＭ  Ｐｒｅｓｓ）を使って得られた、トーン・マップされたＧｅｎ−ＨＤＲを指す。
５−  「Ｇｅｎ−ＨＤＲ−Ｒａｄ」は、本発明による、逆トーン・マッピングの第一フェーズだけを使った、生成ＨＤＲ画像を指す。
【００４２】
以下の図８の説明は、上記で定義された画像１〜５を比較することによって得られた青／緑ピクセルの和の値を示す、図９および１０を参照しながらこれらの説明と併せ行うこととする。図９は、右側２列が各々の比較に対する青と緑との値を表している表であり、図１０は、各比較に対するピクセル当たりのコントラスト・ゲイン（正規化された青ピクセルの和）を表した図である。本発明の実施形態による、逆トーン・マッピング各フェーズの、生成ＨＤＲ画像の品質における影響を評定するために、最初にステップ８０において、Ｇｅｎ−ＨＤＲ−ＲａｄとＲｅａｌ−ＨＤＲとが比較される。この最初の比較は、本発明による、逆トーン・マッピングの第一フェーズ適用の影響を評定する。得られた結果のＡが図９および１０に示されている。この結果は、４０．０％の大幅なピクセル当たりコントラスト・ゲイン、および０．１％のごくわずかなコントラスト喪失があることを表している。
【００４３】
本発明による、逆トーン・マッピングの第二フェーズ適用の影響を評定するために、ステップ８１において、Ｇｅｎ−ＨＤＲとＲｅａｌ−ＨＤＲとの比較が行われる。得られた結果Ｂは、図９および１０に示され、７７．４％のずっと大きなピクセル当たりコントラスト・ゲインおよび０．０％に近い無視可能なコントラスト喪失を表している。これは、第二フェーズが、コントラスト・ゲインに３７％のさらなる増大を効果的に加えたことを示している。
【００４４】
ＨＤＲ画像の典型的な用途の一つはＬＤＲ画像の品質の向上なので、ステップ８２では、Ｇｅｎ−ＬＤＲとＩｎｐｕｔ−ＬＤＲとを比較することによってその影響を検証する。得られた結果Ｃは、図９および１０に表され、大幅なコントラスト・ゲインの向上が示されており、これでは４７．２％のゲインでコントラストの喪失はない。これによれば、ＨＤＲに変換する際に元画像のディテールの喪失はない。これに加えて、生成ＨＤＲは、元画像では見えなかったディテールも包含している。この後者の結果は、例えば医療画像分割などの画像改良アプリケーションの助力となり得るであろう。
【００４５】
本発明は、例えば、既存のＬＤＲビデオおよび画像ライブラリを、新規のＨＤＲディスプレイに用いるためのＨＤＲビデオおよび画像に変換するアプリケーションなど、多くの画像処理アプリケーションに応用することができよう。また、本発明は、ＬＤＲ画像の改善にも適用でき、ＬＤＲ画像はまずＨＤＲ画像に変換され、次いで（平滑化、エッジ検出などの）標準的な画像処理による改善が適用され、しかる後、（Ｒｅｉｎｈａｒｄのトーン・マッピング演算子などの）標準的トーン・マッピング技法を使って、ＨＤＲ画像からＬＤＲに戻し変換が行われる。
【００４６】
斯くのごとく、本発明は単一のＬＤＲ画像からＨＤＲ値を再生するための効率的な技法を提供する。
【００４７】
従来の対処技術は高輝度部分だけの拡大なのに対し、本発明では、ダイナミック・レンジは、輝度チャネルの両方のテール（高輝度および低輝度部分）から拡大される。これは、画像の中間部／陰影部に亘っての視認可能なディテールのレベルを増加させる。これは、画像のダイナミック・レンジを、従来技術によるものよりさらに増大する。
【００４８】
逆トーン・マッピングの第一フェーズで、カメラの露光量−ピクセル・カーブを使用することによって、画像のより現実に近い放射輝度が提供され、仕上がり画像の観賞品質が向上される。
【００４９】
本発明によって、医療画像化など多くのアプリケーションにおいて、最新のＨＤＲディスプレイを使って既存のＬＤＲ画像を見ることができる。
【００５０】
さらに、本発明の実施形態による逆トーン・マッピングを用い、より高品質のＨＤＲ画像を生成してそれを順方向トーン・マップしてＬＤＲに戻すことによって、ＬＤＲカメラのアウトプットの品質を向上することができ、これによりコントラストの可視性の大幅な改善が実現できる。その結果として、さらなる処理作業のためのエッジの検出を簡素化することができる。
【００５１】
本発明は多くのアプリケーションを有する。例えば、医療画像化におけるエッジ検出などのさらなる画像処理でのＨＤＲドメインにおける作業の改良など、デジタル画像に対する向上されたオペレーションとして本発明を用いることができる。
【００５２】
さらに一般的に言えば、本発明は、低品質の音声信号のダイナミック・レンジの増大など、任意のデジタル信号に適用できる。
【００５３】
図１１は、本発明の実施形態による、ビデオ画像をＬＤＲ画像からＨＤＲ画像に変換するため使われるコンピュータ・システム９０を示す。コンピュータ・システム９０は、プロセッサ９１と、プロセッサ９１に連結された入力デバイス９２と、プロセッサ９１に連結された出力デバイス９３と、各々がプロセッサ９１連結されたメモリ・デバイス９４および９５とを含む。プロセッサ９１は、中央処理装置（ＣＰＵ）などの処理装置である。入力デバイス９２は、とりわけ、キーボード、マウスなどとすることができる。出力デバイス９３は、とりわけ、プリンタ、プロッタ、ディスプレイ・デバイス（例えばコンピュータ・スクリーン）、磁気テープ、着脱可能ハード・ディスク、フレキシブル・ディスクなどとすることができる。このディスプレイ・デバイスには、図１の１００のディスプレイ領域を含めることができる。メモリ・デバイス９４および９５は、とりわけ、ハード・ディスク、フレキシブル・ディスク、磁気テープ、コンパクト・ディスク（ＣＤ：ｃｏｍｐａｃｔ  ｄｉｓｃ）またはデジタル・ビデオ・ディスク（ＤＶＤ：ｄｉｇｉｔａｌ  ｖｉｄｅｏ  ｄｉｓｃ）などの光記憶デバイス、ダイナミック・ランダム・アクセス・メモリ（ＤＲＡＭ：ｄｙｎａｍｉｃ  ｒａｎｄｏｍ  ａｃｃｅｓｓ  ｍｅｍｏｒｙ）、読み取り専用メモリ（ＲＯＭ：ｒｅａｄ−ｏｎｌｙ  ｍｅｍｏｒｙ）などとすることができる。メモリ・デバイス９５は、コンピュータ実行可能命令を含むコンピュータ・プログラムである、コンピュータ・コード９７を包含する。コンピュータ・コード９７は、ビデオ画像をＬＤＲ画像からＨＤＲ画像に変換するためのアルゴリズムを包含する。プロセッサ９１はコンピュータ・コード９７を実行する。メモリ・デバイス９４はインプット・データ９６を包含する。インプット・データ９６は、コンピュータ・コード９７に要求されるインプットを包含する。出力デバイス９３は、コンピュータ・コード９７からのアウトプットを表示する。メモリ・デバイス９４および９５（または図１１に示されていない一つ以上の追加のメモリ・デバイス）のいずれかまたは両方は、中に具現されたコンピュータ可読プログラムを有する、もしくはその中に格納された他のデータを有する、またはその両方を有するコンピュータ可用記憶媒体（またはプログラム記憶デバイス）として用いることができ、このコンピュータ可読プログラムはコンピュータ・コード９７を含む。一般に、コンピュータ・システム９０のコンピュータ・プログラム製品（あるいは、製造品）は、上記コンピュータ可用記憶媒体（または上記プログラム記憶デバイス）を含み得る。
【００５４】
図１１は、コンピュータ・システム９０をハードウエアおよびソフトウエアの特定の構成として示しているが、当業者なら周知のように、前述の目的のために、図１１の特定のコンピュータ・システム９０と組み合わせて、ハードウエアおよびソフトウエアの任意の構成を用いることができる。例えば、メモリ・デバイス９４と９５とを、別個のメモリ・デバイスでなく単一のメモリ・デバイスの部分とすることもできる。
【００５５】
本明細書では、例示説明を目的として本発明の特定の実施形態を説明してきたが、当業者には多くの変更および変形が明白になっていよう。
【図１】
2013507674.tif 000003
【図２】
2013507674.tif 000004
【図３】
2013507674.tif 000005
【図４】
2013507674.tif 000006
【図５】
2013507674.tif 000007
【図６】
2013507674.tif 000008
【図７】
2013507674.tif 000009
【図８】
2013507674.tif 000010
【図９】
2013507674.tif 000011
【図１０】
2013507674.tif 000012
2013507674.tif 000013
【図１１】
2013507674.tif 000014
【手続補正書】
【提出日】平成24年4月10日(2012.4.10)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
デジタル画像を、所与のカメラを使って得た低ダイナミック・レンジ（ＬＤＲ）画像から、高ダイナミック・レンジ（ＨＤＲ）画像に変換する方法であって、前記方法は、
−  前記所与のカメラに対する露光量−ピクセル応答カーブを取得するステップ（２１）と、
−  前記ＬＤＲ画像をＨＳＢ色空間アレイに変換するステップ（２２）であって、前記ＨＳＢ色空間アレイは、色相アレイ、彩度アレイ、および明度アレイを含む、前記変換するステップと、
−  前記露光量−ピクセル応答カーブの逆関数を使って、前記明度アレイ中の各ピクセルを逆マッピングすることによって放射輝度アレイを算定するステップ（２３、２４）と、
を含む方法。
【請求項２】
−  前記放射輝度アレイ中の各ピクセルに対し局部輝度平均を算定し（５０）、前記ピクセルに対する前記放射輝度アレイの値と前記局部輝度平均とから、各ピクセルに対する調整された放射輝度アレイを計算するステップ（５１）と、
−  前記放射輝度アレイを、前記ＨＳＢアレイの前記色相および彩度アレイと組み合わせるステップ（５２）と、
−  前記画像を、前記色相、彩度、および放射輝度アレイを使い画像色空間に変換し、前記ＨＤＲ画像を提供するステップ（５３）と、
をさらに含む、請求項１に記載の方法。
【請求項３】
前記露光量−ピクセル応答カーブはカメラのデータシートから取得される、請求項１または２に記載の方法。
【請求項４】
前記露光量−ピクセル応答カーブは、同一の光景に対し異なる露光をされた一連の画像から取得される、請求項１または２に記載の方法。
【請求項５】
各ピクセル（ｕ，ｖ）に対し、前記調整された放射輝度Ｒａｄｉａｎｃｅ’［ｕ，ｖ］は次式により計算され、
Ｒａｄｉａｎｃｅ’［ｕ，ｖ］  ＝  Ｒａｄｉａｎｃｅ［ｕ，ｖ］＊Ｒａｄｉａｎｃｅ［ｕ，ｖ］／Ｌｏｃａｌ＿Ｌｕｍｉｎａｎｃｅ＿Ａｖｅｒａｇｅ［ｕ，ｖ］、
Ｒａｄｉａｎｃｅ［ｕ，ｖ］は、ピクセル（ｕ，ｖ）に対する前記放射輝度の値を指し、Ｌｏｃａｌ＿Ｌｕｍｉｎａｎｃｅ＿Ａｖｅｒａｇｅ［ｕ，ｖ］は、ピクセル（ｕ，ｖ）における前記局部輝度平均を指す、請求項１〜４のいずれかに記載の方法。
【請求項６】
前記画像を画像色空間に変換する前記ステップは、次式により遂行され、
Ｃｏｎｖｅｒｔ（Ｈｕｅ［  ］，Ｓａｔｕｒａｔｉｏｎ［  ］，Ｒａｄｉａｎｃｅ［  ］／ｍａｘ（Ｒａｄｉａｎｃｅ［  ］））＊ｍａｘ（Ｒａｄｉａｎｃｅ）、
Ｈｕｅ［  ］は前記色相アレイを、Ｓａｔｕｒａｔｉｏｎ［  ］は前記彩度アレイを、Ｒａｄｉａｎｃｅ［  ］は前記放射輝度アレイを指す、請求項１〜５のいずれかに記載の方法。
【請求項７】
前記局部輝度平均は次式で定義されるｉ畳み込みカーネルから算定され、
ｋｅｒｎｅｌｉ［  ］＝ＧａｕｓｓｉａｎＫｅｒｎｅｌ（ｒｉ）［  ］、
ｒｉはＧａｕｓｓｉａｎＫｅｒｎｅｌの半径であり、ｉは０〜８の間で変化する前記局部コントラストのスケール指標である、請求項１〜６のいずれかに記載の方法。
【請求項８】
ｒｉは１／（２＊Ｓｑｒｔ（２））＊１．６ｉに設定される、請求項７に記載の方法。
【請求項９】
コンピュータ・プログラムがコンピュータ・デバイス上で実行されたとき請求項１〜７のいずれか一つに記載の前記方法の前記ステップを実行するための命令を含む、前記コンピュータ・プログラム。
【請求項１０】
請求項１〜８のいずれか一つに記載の前記方法の前記ステップを実行するようになされた手段を含む、システム。
【国際調査報告】
2013507674.tif 000015
2013507674.tif 000016
2013507674.tif 000017