(57)【要約】【課題】  画像の劣化過程を表す関数がガウス型関数である観測画像（劣化画像）を復元してぼけのない原画像にできるだけ近い復元画像を得ることができる画像復元方法を提供する。【解決手段】  劣化画像を復元してぼけのない原画像にできるだけ近い復元画像を得る画像復元方法であって、劣化画像をＹＵＶの色空間に変換する準備処理工程（ＳＰ１〜ＳＰ７）と、この準備処理工程によって得たＹＵＶの色空間のうちのＹ画像（Ｙプレーン）を複数の分解画像に分解する分解処理工程（ＳＰ８〜ＳＰ１０）と、この分解処理工程によって得た分解画像を再構成してＹ画像（Ｙプレーン）とすると共に前記準備処理工程で変換したＵＶ画像（ＵＶプレーン）と合成して復元画像とする再構成処理工程（ＳＰ１１〜ＳＰ１４）を備えた。【選択図】        図２
2018028868.tif 000002
【特許請求の範囲】
【請求項１】
劣化画像を復元してぼけのない原画像にできるだけ近い復元画像を得る画像復元方法であって、劣化画像をＹＵＶの色空間に変換する準備処理工程と、この準備処理工程によって得たＹＵＶの色空間のうちのＹ画像（Ｙプレーン）を複数の分解画像に分解する分解処理工程と、この分解処理工程によって得た分解画像を再構成してＹ画像（Ｙプレーン）とすると共に前記準備処理工程で変換したＵＶ画像（ＵＶプレーン）と合成して復元画像とする再構成処理工程を備え、
前記分解処理工程は、以下の工程１〜３からなり、
工程１．Ｙ画像ｆ0に拡散値σ02のガウシアンカーネルＧσ02を適用した画像Ｇｆ0*σ02）をＹ画像ｆ0から引いて１回分解後の分解画像ｆ1を生成する（ｆ1＝ｆ−ｆ0*Ｇσ02）。
工程２．１回分解後の分解画像ｆ1に拡散値σ12のＬＯＧカーネルＬＯＧσ12を適用した画像ｆ1*ＬＯＧσ12を１回分解後の分解画像ｆ1から引いて２回分解後の分解画像ｆ2を生成する（ｆ2＝ｆ1−ｆ1*ＬＯＧσ12）。
工程３．工程２と同様な作業を繰り返して分解画像ｆJまで生成し、分解画像をＪ枚得る。
前記再構成処理工程は、以下の工程４〜７からなることを特徴とする画像復元方法。
工程４．分解画像ｆ1〜ｆJにそれぞれ前記準備処理工程で作成した復元ＬＯＧカーネルＬＯＧσ12-t〜ＬＯＧσJ2-tを適用し、全ての画像の和を求めて総和画像を作成する。
工程５．Ｙ画像に前記準備処理工程で作成した復元ガウシアンカーネルＧσ02-tを適用し、工程４で作成した総和画像と加算する。
工程６．Ｙ画像を入力画像ｆ(x,y)のビット深度に戻す。
工程７．工程６で得た画像をＹプレーンとして入力画像ｆ(x,y)のＵＶプレーンと合成して復元画像Ｉ'(x,y)として出力する。
【請求項２】
請求項１に記載の画像復元方法において、前記分解処理工程及び前記再構成処理工程は、複数のラインバッファにより同時並行処理されることを特徴とする画像復元方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像の劣化過程を表す関数がガウス型関数である観測画像（劣化画像）を復元して、ぼけのない原画像にできるだけ近い復元画像を得る画像復元方法に関する。
【背景技術】
【０００２】
カメラのレンズを通して得た画像が劣化する主な原因としては、カメラと被写体の相対運動によるもの、焦点ずれあるいはピンぼけと呼ばれるもの、厚い空気層など媒質のランダムな揺らぎによるもの、レンズの収差によるものなどが挙げられる。それぞれに対して劣化過程を表す関数が異なる。ここでは、レンズの収差により劣化した画像を復元する場合を取り上げる。レンズの収差による画像の劣化過程を表す関数は、ガウス型関数である。
【０００３】
特許文献１には、既知のＰＳＦ（点拡がり関数）及び１枚の劣化画像に基づき、真の画像を復元することにより、リンギングやノイズを抑えた復元画像を生成するようにした画像復元装置が記載されている。
【０００４】
特許文献１に記載の画像復元装置は、劣化画像に基づき、ベース画像生成処理を行うことにより、ベース画像を生成するベース画像生成処理部と、ベース画像生成処理部からの生成済みベース画像とＰＳＦ（点拡がり関数）との畳み込み積分を行うことにより、ぼかし画像を生成するぼかし画像生成処理部と、生成されたぼかし画像と、前記劣化画像との差分を計算することにより、劣化残差画像を生成する劣化残差画像生成処理部と、生成された劣化残差画像に対して、生成済みベース画像に基づき、拘束の強さを変化させるように画像復元処理を行うことにより、復元残差画像を生成する画像復元処理部と、生成された復元残差画像に生成済みベース画像を加えることにより、復元画像を生成する復元画像生成処理部とを備えている。
【０００５】
また、非特許文献１には、画像処理/信号処理に関して、ウェーブレットの工学的応用について記載されている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００９-２７１７２５号公報
【非特許文献】
【０００７】
【非特許文献１】「ウェーブレットによる信号処理と画像処理」                          中野 宏毅・山本 鎮男・吉田 靖夫 著                          １９９９年８月発行  共立出版
【非特許文献２】Open ＣＶ 2.2 documentation
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかし、特許文献１に開示された技術では、リンギングやノイズを抑えた復元画像を生成することができても、レンズの収差により劣化した画像を復元してぼけのない原画像にできるだけ近い復元画像を得ることは難しい。
【０００９】
本発明は、従来の技術が有するこのような問題点に鑑みてなされたものであり、画像の劣化過程を表す関数がガウス型関数である観測画像（劣化画像）を復元して、ぼけのない原画像にできるだけ近い復元画像が得られる画像復元方法を提供することを目的とする。
【課題を解決するための手段】
【００１０】
上記課題を解決すべく本発明は、劣化画像を復元してぼけのない原画像にできるだけ近い復元画像を得る画像復元方法であって、劣化画像をＹＵＶの色空間に変換する準備処理工程と、この準備処理工程によって得たＹＵＶの色空間のうちのＹ画像（Ｙプレーン）を複数の分解画像に分解する分解処理工程と、この分解処理工程によって得た分解画像を再構成してＹ画像（Ｙプレーン）とすると共に前記準備処理工程で変換したＵＶ画像（ＵＶプレーン）と合成して復元画像とする再構成処理工程を備えた。
【００１１】
前記画像復元方法において、前記準備処理工程は、例えば以下の１〜７である。
１．劣化画像（入力画像）ｆ(x,y)をＹＵＶの色空間に変換する。
２．ＹＵＶの色空間のうちのＹ画像ｆを３２ビット化する。
３．入力パラメータ（復元深度ｔ、比率ｋ）からσ0〜Jの拡散値配列を作成する。
４．拡散値σ02のガウシアンカーネルＧσ02を作成する。
５．拡散値σ12〜σJ2のＬＯＧカーネルＬＯＧσ12〜ＬＯＧσJ2を作成する。
６．上記３の配列内容を単純に、σを−ｔとした復元深度配列を作成する。
７．上記４〜５と同様の作業を上記６の腹元深度配列で行い、復元カーネルを作成する。
また、前記分解処理工程は、例えば以下の８〜１０である。
８．Ｙ画像ｆ0に上記４の拡散値σ02のガウシアンカーネルＧσ02を適用した画像Ｇσ02*ｆ0をＹ画像ｆ0から引いて１回分解後の分解画像ｆ1を生成する（ｆ1＝ｆ0−ｆ0*Ｇσ02）。
９．１回分解後の分解画像ｆ1に上記５の拡散値σ12のＬＯＧカーネルＬＯＧσ12を適用した画像ｆ1*ＬＯＧσ12を１回分解後の分解画像ｆ1から引いて２回分解後の分解画像ｆ2を生成する（ｆ2＝ｆ1−ｆ1*ＬＯＧσ12）。
１０．上記９と同様な作業を繰り返して分解画像ｆJまで生成し、分解画像をＪ枚得る。
更に、前記再構成処理工程は、例えば以下の１１〜１４である。
１１．分解画像ｆ1〜ｆJにそれぞれ上記７で作成した復元ＬＯＧカーネルＬＯＧσ12-t〜ＬＯＧσJ2-tを適用し、全ての画像の和を求めて総和画像を作成する。
１２．Ｙ画像に上記７で作成した復元ガウシアンカーネルＧσ02-tを適用し、上記１１で作成した総和画像と加算する。
１３．Ｙ画像を入力画像ｆ(x,y)のビット深度に戻す。
１４．上記１３で得た画像をＹプレーンとして入力画像ｆ(x,y)のＵＶプレーンと合成して復元画像Ｉ'(x,y)として出力する。
【発明の効果】
【００１２】
本発明に係る画像復元方法によれば、レンズの収差による劣化画像や厚い空気層など媒質のランダムな揺らぎを通した劣化画像など、画像の劣化過程を表す関数がガウス型関数である観測画像（劣化画像）を復元してぼけのない原画像にできるだけ近い復元画像を得ることができる。
【図面の簡単な説明】
【００１３】
【図１】本発明に係る画像復元方法を実施する画像復元装置の構成図
【図２】本発明に係る画像復元方法の処理手順を示すフローチャート
【図３】分解処理・復元処理におけるラインバッファ操作方法の説明図
【図４】（ａ）はレンズ収差により劣化した画像、（ｂ）は復元後の画像
【図５ａ】はテストチャートの原画
【図５ｂ】は原画のぼかし画像
【図５ｃ】はぼかし画像の復元画像
【発明を実施するための形態】
【００１４】
以下に本発明の実施の形態を添付図面に基づいて説明する。本発明に係る画像復元方法を実施する画像復元装置１は、図１に示すように、劣化画像ｆ(x,y)をＹＵＶの色空間に変換する準備処理手段２と、この準備処理手段２によって得たＹＵＶの色空間のうちのＹ画像（Ｙプレーン）ｆ0を複数の分解画像ｆ1〜ｆJに分解する分解処理手段３と、この分解処理手段３によって得た分解画像ｆ1〜ｆJを再構成してＹ画像（Ｙプレーン）ｆとすると共に準備処理手段２で変換したＵＶ画像（ＵＶプレーン）ｆuvと合成して復元画像Ｉ'(x,y)とする再構成処理手段４を備える。
【００１５】
レンズと撮像素子（撮像センサ）による画像は、一般にレンズの収差によるガウス型のボケと高周波によるノイズによって劣化している。本来、１ピクセルに集まるべき光は、ガウス分布でボケており、周辺のピクセルに影響している。全く収差のないレンズによる画像は、ボケを生じないが、実際の画像はガウス分布でボケを生じた光の集合となる。
【００１６】
画像を数式で表すために、先ず二次元ガウス関数をＧt(x,y)とすると、
Ｇt(x,y)＝1/2πt*exp(-(x2+y2)/2t)        （式１）
ｔはガウス関数の広がりを示す拡散値である。
【００１７】
収差のない（ボケのない）画像（原画像）をＩ(x,y)とすると、ボケた画像（劣化画像）の劣化過程ｆ(x,y)は次式のように表される。
ｆ(x,y)＝ＫtＩ(x,y)＋ω(x,y)＝Ｇt(x,y)*Ｉ(x,y)＋ω(x,y)      （式２）
＊はたたみ込みを表す演算子であり、次式で定義される。
【００１８】
Ｇt(x,y)*Ｉ(x,y)＝∬Ｇt(x',y')*Ｉ(x・x',y・y')dx'dy'          （式３）
これを有限画素で計算するために、次式に変換する。
Ｇt(x,y)*Ｉ(x,y)＝Σy'Σx'(Ｇt(x',y')*Ｉ(x・x',y・y'))        （式４）
なお、ω(x,y)は、観測時の高周波雑音と考えられるので省略する。
【００１９】
ここで、Ｋtは関数を二次元ガウス関数Ｇt(x,y)でたたみ込みをする演算子である。ｔはガウス関数の広がりを表しているので、ｔを時刻と考えるとＫtは例えば一点にあるインクが次第に拡散していく様子を表すと考えることができる。拡散時間は当然ｔ≧０であるが、これをｔ＜０とするとガウス関数拡散の逆、つまり復元の過程と捉えることができる。
【００２０】
このような操作をする演算子をＫ-tと表すと、原画像Ｉ(x,y)の近似解Ｉ'(x,y)を次式のように表すことができる。
Ｉ'(x,y)＝Ｋ-tｆ(x,y)＝Ｉ(x,y)＋Ｋ-tω(x,y)      （ｔ＞０）  （式５）
また、ＫtとＫ-tの関係は次式となる。
ＫtＫ-tｆ(x,y)＝ｆ(x,y)            （式６）
【００２１】
このままでは利用することができないので、以下の関係式を用いる。
ひとつは、ガウス関数のたたみ込みに対しての関係式
Ｋ-tＧσ2*ｆ＝Ｇσ2-t*ｆ            （σ2＞ｔ）        （式７）
もうひとつは、ＬＯＧ関数のたたみ込みに対しての関係式
Ｋ-tＬＯＧσ2*ｆ＝ＬＯＧσ2-t*ｆ    （σ2＞ｔ）        （式８）
【００２２】
ＬＯＧ関数とは、画像処理でよく利用されるラプラシアンオブガウシアンフィルターであり、二次元のガウス関数にラプラシアン演算子Ｌを作用させることで得られる。
ラプラシアン演算子Ｌ＝∂2/∂x2＋∂2/∂y2              （式９）
【００２３】
ラプラシアン演算子Ｌをガウス関数に作用させると次式となる。
ＬＯＧσ2(x,y)＝１/(２πσ4)*(２−(x2＋y2)/σ2)*exp(−(x2＋y2)/２σ2)  （式１０）
この関数（式１０）をフィルタと捉えると、中間周波数を通す帯域通過フィルタとなる。
【００２４】
式７と式８を利用して劣化画像ｆ(x,y)を分解すると、
ｆ＝Ｇσ02*ｆ1＋ΣJj=1ＬＯＧσj2*ｆj＋ｆJ+1              （式１１）
本案では、劣化画像ｆ(x,y)を有限画素で計算するために次式に変換する。
ｆ(x,y)＝Σy'Σx'(Ｇσ02-t(x',y')*ｆ1(x・x',y・y'))＋
ΣJj=1(Σy'Σx'(ＬＯＧσj2-t(x',y')*ｆj(x・x',y・y')))＋ｆJ+1    （式１２）
ここで、σ2はガウス関数およびＬＯＧ関数のスケールパラメータである。
【００２５】
ｆjはｊ回分解後の残差画像であり、次式で得られる。
１段階目
ｆ1＝ｆ−Ｇσ02*ｆ        （式１３）
ｊ段階目
ｆj＝ｆj-1−ＬＯＧσj-12*ｆj-1      （１≦ｊ≦Ｊ）        （式１４）
【００２６】
式７、式８、式１１から復元を式で表すと
Ｋ-tｆ＝Ｇσ2-t*ｆ＋ΣJj=1ＬＯＧσ2-t*ｆj＋Ｋ-tｆJ+1          （式１５）
となるが、式を安定させるには
σJ2＞ｔ，σJ+12≦ｔ                （条件１）
でなければならない。
【００２７】
この時、式１５のＫ-tｆJ+1は不安定であるが、この項は劣化画像の最も高い周波数成分であるため、ほとんどが雑音である。よって、この項は無視することとする。よって、求める復元画像Ｉ'(x,y)は次式となる。
Ｉ'(x,y)＝Ｇσ2-t*ｆ＋ΣJj=1ＬＯＧσ2-t*ｆj          （式１６）
【００２８】
Ｊは劣化画像の分解回数、ｔは復元深度である。
σjはｊ回目の拡散(分解周波数)であり、σj+1＝ｋσj（ｋ＜１）  （式１７）と等比数列になる。
【００２９】
本実施例では、式１６を有限画素で計算するため、復元画像Ｉ'(x,y)を次式に変換する。
Ｉ'(x,y)＝Σy'Σx'(Ｇσ02-t(x',y')*ｆ(x・x',y・y'))＋
ΣJj=1(Σy'Σx'(ＬＯＧσj2-t(x',y')*ｆj(x・x',y・y')))  （式１８）
【００３０】
式１２および式１８は、離散的プログラムで実装可能な式となっているが、劣化画像ｆ(x,y)に対して６つのパラメータを決定する必要がある。即ち、分解回数Ｊ、復元深度ｔ、ガウシアン/ＬＯＧカーネルサイズｎ、拡散の初期値σ0，σ1と比率ｋである。
【００３１】
カーネルは有限画素のたたみ込み計算を行う(x',y')の行列であり、中心のあるｎ×ｎ（ｎは奇数）の正方行列とする。これらは上記条件１から特定の関係にある。
σ1≧1.44σ0
σj+1≧0.52σj
分解回数が小さいほど処理が終了するため、処理速度を優先すると、
σ1＝1.44σ0
ｋ＝0.52
となる。但し、比率ｋが粗い（小さい）と分解精度が悪くなり、十分な復元が得られない。
【００３２】
カーネルサイズ（正方行列の大きさ）を自動的に計算するには以下のように行う。
実施例のカーネルはガウシアン／ＬＯＧの２次元行列である。ガウシアン／ＬＯＧはいずれも行列の中心を原点として、中心から遠いほど値はゼロに近づく性質を持つ。
そこで、中心からの距離をｒとし、半径ｒの円周で1/256未満になるσは簡易的にσ＝０．３・ｒ＋０．８と近似できる。
σからカーネルサイズを求めるには、上記の式から、
ｒ＝（σ−０．８）/０．３＋１と近似できる。
上記の式の右端の（＋１）は製法行列とする場合に各外辺で丸め誤差がでないよう、１要素分だけ大きい正方行列にするために必要となる。
これから、カーネルサイズは（ｒ×２）＋１と算出できる。
【００３３】
比率ｋに対して忠実度をｄＢで表現すると、次のようになる。
ｋ＝0.7：４５ｄＢ、ｋ＝0.8：４９ｄＢ、ｋ＝0.9：５５ｄＢ、ｋ＝0.95：５９ｄＢとなる。忠実度が５０ｄＢを確保するには、ｋ＝0.8以上が必要である。
【００３４】
必要条件は上記の通りであるが、実用になるパラメータを迅速に決定するためには、パラメータを自動算出する必要がある。
本案ではパラメータを復元量（復元深度）ｔと精度（比率）ｋの２つを決定すると他のパラメータは算術的に決定されるアルゴリズムを用いて自動算出を行う。
１. 拡散の初期値σ0＝√(ｔ+0.5)    経験的に求めた数値
２. 拡散の初期値σ1＝1.44σ0        非特許文献１より
３. ガウシアン/ＬＯＧカーネルサイズｎ＝2((σ1-0.8)/0.3+1.0)+1
非特許文献２よりσ＝0.3(ｎ/2-1)+0.8から逆算して概算奇数
４. 分解回数Ｊ＝σj≦0.7となるｊ    ＬＯＧの中心が(0±0.5,0±0.5)以内に収まるσ＝0.7、精度ｋが1.0に近い場合、分解回数Ｊが非常に大きくなる可能性がある。
【００３５】
この場合、式１７であるσj+1＝ｋσj（ｋ＜１）を次式のように等差数列に置き換えて分解回数を削減する。
σj＝σ1・(σ1・ｋσ1)・(ｊ−1)              （式１９）
【００３６】
復元処理のアルゴリズムは、主に準備処理、分解処理と再構成処理の３つに分かれる。
先ず、準備処理では、以下の処理が必要となる。復元画像Ｉ'(x,y)を表す式１８を処理するには、ｆjが必要である。ｆjは式１２によって算出される。式１２と式１８では、拡散値σ02のガウシアンカーネルと拡散値σ12〜σJ2のＪ個のＬＯＧカーネルが必要である。これらのカーネルは、入力パラメータである復元量（復元深度）ｔと精度（比率）ｋから自動計算される拡散の初期値σ0,σ1より、式１及び式１０で算出する。
【００３７】
σ0〜σJは再構成処理のため、式１７により配列等に記憶しておく。但し、分解回数Ｊが一定以上大きい場合（実装のメモリ容量などで規定）は、式１９を利用する。
【００３８】
次に、分解処理では、以下の処理が行われる。カーネルが用意できたら式１２の前半に従って、劣化画像ｆ(x,y)と、劣化画像ｆ(x,y)に拡散値σ02のガウシアンカーネルＧσ02を適用したＧσ02*ｆ0との差を求め、これをｆ1とする。次いで、式１２の後半に従って、ｆ1と、ｆ1に拡散値σ12のＬＯＧカーネルＬＯＧσ12を適用したＬＯＧσ12*ｆ1との差を求め、これをｆ2とする。以下ｆJまでこの処理を繰り返し、劣化画像ｆ(x,y)の分解画像ｆ1〜ｆJをＪ枚得る。
【００３９】
次に、再構成処理では、以下の処理が行われる。式１８に従って、分解画像ｆ1〜ｆJに対して拡散値σ12-t〜σJ2-tの復元ＬＯＧカーネルを適用し、全ての画像の和を求める。次いで、劣化画像ｆ(x,y)に拡散値σ02-tの復元ガウシアンカーネルを適用し、その画像に上記の復元ＬＯＧカーネルを適用した全ての画像を加算して、復元画像Ｉ'(x,y)とする。
【００４０】
ここでは、グレースケールを前提としたアルゴリズムである。従って、入力画像（劣化画像）は先ずＹＵＶの色空間に変換され、Ｙ画像（Ｙプレーン）のみに対して上記復元処理が適用される。
【００４１】
次に、本発明に係る画像復元方法について、図２に示すフローチャートに従って説明する。先ず、劣化画像（入力画像）ｆ(x,y)をＹＵＶの色空間に変換する（ステップＳＰ１）。次いで、ＹＵＶの色空間のうちのＹ画像（Ｙプレーン）ｆ0を３２ビット化する（ステップＳＰ２）。
【００４２】
ステップＳＰ３では、入力パラメータ（復元深度ｔ、比率ｋ）からσ0〜Jの拡散値配列を作成する（式１７または式１９）。ステップＳＰ４では、拡散値σ02のガウシアンカーネルＧσ02を作成する（式１）。ステップＳＰ５では、拡散値σ12〜σJ2のＬＯＧカーネルＬＯＧσ12〜ＬＯＧσJ2を作成する（式１０）。
【００４３】
ステップＳＰ６では、拡散値配列の配列内容を単純に、σを−ｔとした復元深度配列を作成する。ステップＳＰ７では、ステップＳＰ４,ＳＰ５と同様の作業をステップＳＰ６の腹元深度配列で行い、復元カーネルを作成する。ここまで、ステップＳＰ１〜ステップＳＰ７が準備処理工程である。
【００４４】
次いで、ステップＳＰ８では、Ｙ画像ｆ0にステップＳＰ４の拡散値σ02のガウシアンカーネルＧσ02を適用した画像Ｇσ02*ｆ0をＹ画像ｆ0から引いて１回分解後の分解画像ｆ1を生成する（ｆ1＝ｆ0−Ｇσ02*ｆ0）。ステップＳＰ９では、（ｊ−１）回分解後の分解画像ｆj-1にステップＳＰ５の拡散値σj-12のＬＯＧカーネルＬＯＧσj-12を適用した画像ＬＯＧσj-12*ｆj-1を（ｊ−１）回分解後の分解画像ｆj-1から引いてｊ回分解後の分解画像ｆjを生成する（ｆj＝ｆj-1−ＬＯＧσ12*ｆj-1）。
【００４５】
ステップＳＰ１０では、ｊ（１≦ｊ≦Ｊ）がＪと等しいか否かが判断され、ｊ＜Ｊの場合にはステップＳＰ９に戻って分解画像ｆjを生成し、ｊ＝Ｊの場合にはステップＳＰ１１へ進む。ここまで、ステップＳＰ８〜ステップＳＰ１０が分解処理工程であり、分解画像をＪ枚得る。
【００４６】
次いで、ステップＳＰ１１では、分解画像ｆ1〜ｆJにそれぞれステップＳＰ７で作成した復元ＬＯＧカーネルＬＯＧσ12-t〜ＬＯＧσJ2-tを適用し、全ての画像の和を求めて総和画像を作成する。ステップＳＰ１２では、Ｙ画像にステップＳＰ７で作成した復元ガウシアンカーネルＧσ02-tを適用し、ステップＳＰ１１で作成した総和画像と加算する。
【００４７】
ステップＳＰ１３では、Ｙ画像を劣化画像（入力画像）ｆ(x,y)のビット深度に戻す。ステップＳＰ１４では、ステップＳＰ１３で得た画像をＹプレーンとして劣化画像（入力画像）ｆ(x,y)のＵＶプレーンと合成し、復元画像Ｉ'(x,y)として出力する。ステップＳＰ１１〜ステップＳＰ１４が再構成処理工程である。
【００４８】
図３は分解処理・復元処理におけるラインバッファ操作方法を示し、分解処理工程（ステップＳＰ８〜ステップＳＰ１０）及び再構成処理工程（ステップＳＰ１１）において、Ｊ＝３の場合である。
【００４９】
Ｙ画像ｆ0から分解画像ｆ1、分解画像ｆ2、分解画像ｆ3を生成する作業及び分解画像ｆ1〜ｆ3に復元ＬＯＧカーネルＬＯＧσ12-t〜ＬＯＧσ32-tを適用し、全ての画像の和を求めて総和画像を作成する作業が、処理時刻（１〜８）毎に各ラインバッファ（ＬＢＬＶ１〜ＬＢＬＶ４）において並行処理される。
ここで、ラインバッファは入力画像の水平画素数×カーネル縦サイズの大きさを持つ一時記憶メモリである。
【００５０】
具体的には、最初の時刻（図３の最上段）では、ラインバッファ（ＬＶ１）からＹ画像ｆ0を読み込み、ラインバッファ（ＬＶ２）に分解画像ｆ1を書き込む。
次の時刻（図３の上から２段目）では、ラインバッファ（ＬＶ１）にリードモディファイ（ｆ0−＝ｆ0*Ｇ−）を書き込み、ラインバッファ（ＬＶ３）に分解画像ｆ2を書き込む。
次の時刻（図３の上から３段目）では、ラインバッファ（ＬＶ２）にリードモディファイ（ｆ1−＝ｆ1−*Ｌ１−）を書き込み、ラインバッファ（ＬＶ４）に分解画像ｆ3を書き込む。
次の時刻（図３の上から４段目）では、ラインバッファ（ＬＶ３）にリードモディファイ（ｆ2−＝ｆ2−*Ｌ２−）を書き込み、ラインバッファ（ＬＶ４）にリードモディファイ（ｆ3−＝ｆ3−*Ｌ３−）を書き込む。
次の時刻（図３の上から５段目）では、ラインバッファ（ＬＶ４）から（ｆ3−）を読み取り、ラインバッファ（ＬＶ３）に（ｆ2−）＋（ｆ3−）を書き込む。
次の時刻（図３の上から６段目）では、ラインバッファ（ＬＶ３）から（ｆ2−）＋（ｆ3−）を読み取り、ラインバッファ（ＬＶ２）に（ｆ1−）＋（ｆ2−）＋（ｆ3−）を書き込む。
次の時刻（図３の上から７段目）では、ラインバッファ（ＬＶ２）から（ｆ1−）＋（ｆ2−）＋（ｆ3−）を読み取り、ラインバッファ（ＬＶ１）に（ｆ0−）＋（ｆ1−）＋（ｆ2−）＋（ｆ3−）を書き込む。
このように処理時刻毎に各ラインバッファにおいて並行処理されることにより、分解処理及び復元処理の迅速化が可能になる。
【００５１】
図４はレンズ収差により劣化した画像の（ａ）は復元前の画像、（ｂ）は復元後の画像である。
また、図５は一般的に用いられるテストチャートを用いたものであり、（ａ）は原画、（ｂ）は原画のぼかし画像、（ｃ）はぼかし画像を本発明方法で復元した画像である。
これらの画像から、本発明方法によりレンズの収差によるガウス型のボケによる劣化画像が復元されていることが分かる。
【産業上の利用可能性】
【００５２】
本発明によれば、画像の劣化過程を表す関数がガウス型関数である観測画像（劣化画像）を復元してぼけのない原画像にできるだけ近い復元画像を得ることができる画像復元方法を提供することができる。
【符号の説明】
【００５３】
１…画像復元装置、２…準備処理手段、３…分解処理手段、４…再構成処理手段。
【図１】
2018028868.tif 000003
【図２】
2018028868.tif 000004
【図３】
2018028868.tif 000005
【図４】
2018028868.tif 000006
【図５ａ】
2018028868.tif 000007
【図５ｂ】
2018028868.tif 000008
【図５ｃ】
2018028868.tif 000009