ノイズ低減装置、ノイズ低減方法及びプログラム
【課題】ノイズを含んでいる一枚の画像に対し、位置毎に異なるノイズ低減処理を簡単に行えるようにしたノイズ低減装置を提供する。
【解決手段】ノイズ低減装置は、画像から画像パッチ集合を生成する、画像パッチ集合生成部と、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部と、ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成する、画像合成部と、
を備える。
【解決手段】ノイズ低減装置は、画像から画像パッチ集合を生成する、画像パッチ集合生成部と、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部と、ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成する、画像合成部と、
を備える。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、デジタル画像処理技術に関し、特に、デジタル画像のノイズを低減し、デノイズされた画像を得るためのノイズ低減装置、ノイズ低減方法及びプログラムに関する。
【背景技術】
【０００２】
近年、画像取得装置の一つとして、デジタルカメラは急速に普及している。高品位な画質を得るために、デジタルカメラで撮影されたデジタル画像に対し、様々なデジタル画像処理が施される。数多くのデジタル画像処理（以下、単に、「画像処理」と言う。）の中で、ノイズ低減処理が画質に大きな影響を与えるため、大変重要な画像処理の一つであることが広く知られている。
【０００３】
従来では、例えば、非特許文献１、非特許文献２、及び特許文献１に開示されたように、デジタル画像（以下、単に、「画像」とも言う。）に含まれるノイズを低減するための画像ノイズ低減技術が数多く提案されている。
【０００４】
特許文献１には、ノイズ低減効果の低い部分の出現を抑え、画像全体として高いノイズ低減効果を得ることができる画像信号処理装置が開示されている。
【０００５】
特許文献１の画像信号処理装置は、画像信号に対して第１の平滑化処理を行うことで、ノイズ低減された第１の平滑化信号を生成する第１の平滑化信号生成部と、画像信号に対して第１の平滑化処理と異なる第２の平滑化処理を行うことで、ノイズ低減された第２の平滑化信号を生成する第２の平滑化信号生成部と、画像信号、第１の平滑化信号、及び第２の平滑化信号のいずれか１つの値に応じた混合割合に従って第１の平滑化信号と第２の平滑化信号とを混合処理する混合処理部とを備えることを主な特徴としている。
【０００６】
つまり、特許文献１では、画像に対してノイズ低減効果の異なる２つの平滑化処理を行い、それら２つの平滑化処理により生成された２つの平滑化信号を所定の混合割合で混合処理することにより、当該画像のノイズを低減するようにしている。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２０１１−８７０８７号公報
【非特許文献】
【０００８】
【非特許文献１】ジェー. ポーティラ（J. Portilla）、ブィ. ストレラ（V. Strela）、エム. ジェー. ウェインライト（M. J. Wainwright）、イー. ピー. シモンセリ（E. P. Simoncelli）共著,「イメージ デノイジング ユージング スケール ミックスチャー オフ ガウシアン イン ザ ウェーブレット ドメイン （Image Denoising using Scale Mixture of Gaussians in the Wavelet Domain）」,ＩＥＥＥ トランスアクションズ オン イメージ プロセシング（IEEE Transactions on Image Processing）,第12巻,第11号,p.1338-1351,2003年
【非特許文献２】エイ. ブアデス（A. Buades）、ビー. コル（B. Coll）、ジェー. エム. モレル（J.M. Morel）共著,「ア ノンローカル アルゴリズム フォー イメージ デノイジング （A non-local algorithm for image denoising）」,ＩＥＥＥ コンピュータ ソサイエティー カンファレンス オン コンピュータ ビジョン アンド パターン レコグニション （IEEE Computer Society Conference on Computer Vision and Pattern Recognition）,第2巻,p.60-65,2005年
【発明の概要】
【発明が解決しようとする課題】
【０００９】
ところで、自然画像は、その画素位置（又は画像領域）によりテクスチャやエッジが大きく異なるという性質を有している。
【００１０】
しかし、従来のノイズ低減技術（例えば、非特許文献１及び非特許文献２を参照）の殆どは、テクスチャ領域・エッジ領域を問わず、画像全体に対して均質的にノイズ低減処理を行うものであり、つまり、画像領域毎（又は画素位置毎）に異なるノイズ低減処理を行っていない。
【００１１】
また、特許文献１では、２つの平滑化処理により生成された２つの平滑化信号を混合処理する際の混合割合は、画像の箇所毎（画素位置毎）に異なる値となっているものの、基本的に、特許文献１において、画像のノイズを低減するために、２種類の平滑化処理しか行っておらず、画像の箇所毎（画素位置毎）に異なる処理とはなっていない。
【００１２】
本発明は、上述のような事情よりなされたものであり、本発明の目的は、ノイズを含んでいる一枚の画像に対し、位置毎に異なるノイズ低減処理を簡単に行えるようにした、ノイズ低減装置、ノイズ低減方法及びプログラムを提供することにある。
【課題を解決するための手段】
【００１３】
本発明は、ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置に関し、本発明の上記目的は、前記画像から画像パッチ集合を生成する、画像パッチ集合生成部と、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部と、前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成部とを備え、前記ノイズ低減処理部は、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出部と、前記画像パッチ集合、及び前記平均ベクトル算出部で算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出部と、固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出部で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出部と、前記固有値固有ベクトル算出部で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成部と、前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出部で算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成部で生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影部とを備えることによって効果的に達成される。
【００１４】
また、本発明の上記目的は、ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置であって、前記画像に対し、フィルタリング処理を行うことにより、フィルタリングされた画像を生成する、フィルタリング処理部と、前記画像から画像パッチ集合を生成するとともに、前記フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成する、画像パッチ集合生成部と、前記画像パッチ集合及び前記フィルタリングされた画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部と、前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成部とを備え、前記ノイズ低減処理部は、前記画像パッチ集合及び前記フィルタリングされた画像パッチ集合に基づき、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出部と、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記平均ベクトル算出部で算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出部と、固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出部で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出部と、前記固有値固有ベクトル算出部で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成部と、前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出部で算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成部で生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影部とを備えることによって効果的に達成される。
【００１５】
また、本発明の上記目的は、ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置であって、予め撮影された一枚以上の画像を格納する画像データベースと、前記画像から画像パッチ集合を生成するとともに、前記画像データベースに格納されている画像からデータベース画像パッチ集合を生成する、画像パッチ集合生成部と、前記画像パッチ集合及び前記データベース画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部と、前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成部とを備え、前記ノイズ低減処理部は、前記画像パッチ集合及び前記データベース画像パッチ集合に基づき、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出部と、前記画像パッチ集合、前記データベース画像パッチ集合、及び前記平均ベクトル算出部で算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出部と、固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出部で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出部と、前記固有値固有ベクトル算出部で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成部と、前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出部で算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成部で生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影部とを備えることによって効果的に達成される。
【００１６】
また、本発明の上記目的は、ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置であって、前記画像に対し、フィルタリング処理を行うことにより、フィルタリングされた画像を生成する、フィルタリング処理部と、予め撮影された一枚以上の画像を格納する画像データベースと、前記フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成するとともに、前記画像データベースに格納されている画像からデータベース画像パッチ集合を生成し、さらに、前記画像から画像パッチ集合を生成する、画像パッチ集合生成部と、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部と、前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成部とを備え、前記ノイズ低減処理部は、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づき、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出部と、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、前記データベース画像パッチ集合、及び前記平均ベクトル算出部で算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出部と、固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出部で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出部と、前記固有値固有ベクトル算出部で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成部と、前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出部で算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成部で生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影部とを備えることによって効果的に達成される。
【００１７】
また、本発明の上記目的は、ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置であって、フィルタリング処理部と、画像データベースと、画像パッチ集合生成部と、ノイズ低減処理部と、画像合成部とを備え、前記ノイズ低減装置では、前記フィルタリング処理部が、前記画像に対し、フィルタリング処理を行うことにより、フィルタリングされた画像を生成し、前記画像パッチ集合生成部が、前記画像から画像パッチ集合を生成するとともに、前記フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成し、前記ノイズ低減処理部が、前記画像パッチ集合及び前記フィルタリングされた画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成し、前記画像合成部が、前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成し、前記デノイズされた画像が、前記画像データベースの画像として、前記画像データベースに格納され、前記画像パッチ集合生成部が、前記画像データベースに格納されている画像に対して、データベース画像パッチ集合生成処理を行うことにより、データベース画像パッチ集合を生成し、前記ノイズ低減処理部が、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づき、前記全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成し、前記画像合成部が、前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成し、繰り返し処理終了条件を満たしたと判定された場合に、前記デノイズされた画像を前記ノイズ低減装置が生成したデノイズされた画像とし、前記繰り返し処理終了条件を満たさないと判定された場合に、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づいて行われるノイズ低減処理を繰り返し行うことにより、デノイズされた画像パッチを生成し、生成したデノイズされた画像パッチを前記画像データベースの画像とすることによって効果的に達成される。
【発明の効果】
【００１８】
本発明に係るノイズ低減装置によれば、ノイズを含んでいる一枚の画像に対し、画像パッチ毎に異なるパラメータによるノイズ低減処理を行った上で、デノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を得るようにしているので、ノイズを含む画像の位置毎に異なるノイズ低減処理を簡単に行うことが可能になり、当該画像の領域（テクスチャ領域・エッジ領域）別に最適なノイズ低減効果を得ることができる。
【００１９】
本発明の発明者らは、異なるレベルのノイズを付加した１２枚のデジタル画像に対し、非特許文献１に記載された従来のノイズ低減方法を用いてノイズ低減を行った結果、非特許文献２に記載された従来のノイズ低減方法を用いてノイズ低減を行った結果、及び、本発明に係るノイズ低減方法を用いてノイズ低減を行った結果を比較することにより、本発明のノイズ低減効果を検証した。
【００２０】
表１は、異なるレベルのノイズが付加されたデジタル画像に対して、非特許文献１によるノイズ低減方法、非特許文献２によるノイズ低減方法、及び本発明によるノイズ低減方法を用いて、ノイズ低減を行って得られた結果画像（即ち、デノイズされた画像）のＰＳＮＲ(Peak Signal-to-Noise Ratio)を示している。
【００２１】
ＰＳＮＲは画質指標の一つであり、ＰＳＮＲの値が高いほど、結果画像（デノイズされた画像）が原画像（ノイズが付加されていないデジタル画像）に近いことを意味する。
【００２２】
表１の結果から、本発明に係るノイズ低減方法によるノイズ低減効果が一番高いことは、はっきりと確認された。
【００２３】
【表１】
【図面の簡単な説明】
【００２４】
【図１】本発明に係るノイズ低減装置の第１実施形態（ノイズ低減装置１）を示すブロック構成図である。
【図２】本発明に係るノイズ低減装置の画像パッチ集合生成部で行われる画像パッチ集合生成処理を説明するための模式図である。
【図３】本発明に係るノイズ低減装置１のノイズ低減処理部２０の実施形態を示すブロック構成図である。
【図４】本発明に係るノイズ低減装置の第２実施形態（ノイズ低減装置２）を示すブロック構成図である。
【図５】本発明に係るノイズ低減装置２のノイズ低減処理部２１の実施形態を示すブロック構成図である。
【図６】本発明に係るノイズ低減装置の第３実施形態（ノイズ低減装置３）を示すブロック構成図である。
【図７】本発明に係るノイズ低減装置３のノイズ低減処理部２２の実施形態を示すブロック構成図である。
【図８】本発明に係るノイズ低減装置の第４実施形態（ノイズ低減装置４）を示すブロック構成図である。
【図９】本発明に係るノイズ低減装置４のノイズ低減処理部２３の実施形態を示すブロック構成図である。
【図１０】本発明に係るノイズ低減装置の第５実施形態（ノイズ低減装置５）を示すブロック構成図である。
【図１１】本発明に係るノイズ低減装置５の処理流れを示すフロー図である。
【発明を実施するための形態】
【００２５】
本発明は、ノイズを含んでいる一枚の画像に対し、位置毎に異なるノイズ低減処理を簡単に行えるようにした、ノイズ低減装置、ノイズ低減方法及びプログラムに関する。
【００２６】
以下、本発明を実施するための形態について、図面を参照しながら詳細に説明する。
【００２７】
図１は、本発明に係るノイズ低減装置の第１実施形態（以下、単に「ノイズ低減装置１」とも言う。）を示すブロック構成図である。以下、図１を用いて、ノイズ低減装置１を詳細に説明する。
【００２８】
図１に示すように、ノイズ低減装置１は、ノイズを含む一枚のデジタル静止画像（以下、単に、「入力画像」又は「画像」とも言う。）から、デノイズされた画像を得るための装置であり、入力画像から画像パッチ集合を生成する画像パッチ集合生成部１０と、画像パッチ集合生成部１０で生成された画像パッチ集合に属する全ての画像パッチに対して、ノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部２０と、ノイズ低減処理部２０で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を得る、画像合成部３０とから構成される。
【００２９】
本発明では、デジタル静止画像であれば、複数チャネル画像・単一チャネル画像を問わず、本発明に係るノイズ低減技術を適用し、その画像のデノイズされた画像を簡単に得ることができる。
【００３０】
また、デジタル静止画像に限らず、動画像のフレームに、本発明に係るノイズ低減技術を適用し、そのフレームのデノイズされた画像（デノイズされたフレーム）を簡単に得ることが可能であり、本発明を動画像のノイズ低減にも適用することができると言える。
【００３１】
以下、説明上の便宜から、本発明のノイズ低減装置に入力されるデジタル画像は、例えば、グレースケール画像のような単一チャネル画像（即ち、単一のカラーチャネルを持つ画像）であることを前提とする。
【００３２】
また、本発明に係るノイズ低減技術を複数チャネル画像（複数のカラーチャネルを持つカラー画像）に適用する場合に、各カラーチャネルに対して、後述のような、単一チャネル画像を入力画像とする本発明に係るノイズ低減装置で行われる処理を繰り返せばよい。
【００３３】
図１に示すように、ノイズ低減装置１では、まず、画像パッチ集合生成部１０が、入力画像から、当該入力画像の画像パッチ集合を生成する。
【００３４】
具体的に、画像パッチ集合生成部１０は、入力画像から画像パッチ集合を得るために、入力画像に対して、画像パッチ集合生成処理を行うことにより、画像パッチ集合を生成する。画像パッチ集合生成処理は、図２に示すように、入力画像に対して、画素位置をずらしながら、複数の画像パッチを生成する処理である。
【００３５】
本発明では、画像パッチ集合生成処理において、画素位置のずらし量を任意に設定することができるが、ずらし量の一具体例として、１画素であり、つまり、画素位置を１画素ずらしながら、複数の画像パッチを生成する。
【００３６】
本発明では、画像パッチ集合生成処理により得られた複数の画像パッチを一つの集合（即ち、画像パッチ集合）としている。画像パッチは、その画像パッチに属する全ての画素の画素値を一列に並べることにより、列ベクトルに変換することができる。つまり、一つの画像パッチは、一つの列ベクトルに対応する。
【００３７】
以下、画像パッチの一具体例として、画像内のｍ画素×ｎ画素の矩形領域を画像パッチとする（図２を参照する。）。説明上の便宜から、以下の説明では、画像パッチを矩形領域（ｍ画素×ｎ画素）としているが、本発明において、画像パッチは矩形領域に限定されることがなく、任意の形状（例えば、多角形）を有する領域を画像パッチとすることができる。
【００３８】
次に、ノイズ低減装置１では、図１に示すように、ノイズ低減処理部２０が、画像パッチ集合生成部１０で生成された画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【００３９】
そして、ノイズ低減装置１では、図１に示すように、画像合成部３０が、ノイズ低減処理部２０で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成する。
【００４０】
ここで、画像合成部３０にて行われる画像合成処理とは、ノイズ低減された全ての画像パッチ（即ち、ノイズ低減処理部２０で得られた全てのデノイズされた画像パッチ）を画像の元の画素位置に戻す処理であり、つまり、デノイズされた画像パッチの画素値を画像の元の画素位置に対応するデノイズされた画像の画素値とする処理である。
【００４１】
ところで、図２にも示されるように、画像パッチ集合生成部１０では、画像パッチ集合生成処理により入力画像から複数の画像パッチを生成する際に、画素位置を所定個数の画素（例えば、１画素）ずらしながら、画像パッチ同士が重なりあうように、画像パッチが生成される。
【００４２】
このように重なりあって生成された各画像パッチに対し、それぞれノイズ低減処理を行うと、画像パッチが重なっている画素において、画素値が一致しないことになる。つまり、画素位置によって、複数のデノイズされた画像パッチに対応する画素があれば、一つのデノイズされた画素パッチに対応する画素もある。
【００４３】
そこで、本発明は、このような画像パッチが重なっている画素（即ち、複数のデノイズされた画像パッチに対応する画素）に対し、これら複数のデノイズされた画像パッチの画素値を平均し、その平均値を画像合成処理後の画素値（即ち、デノイズされた画像の画素値）とするように、画像合成処理を行う。なお、本発明では、複数のデノイズされた画像パッチの画素値を平均する際に、単純平均のかわりに、重み付き平均を用いても良い。
【００４４】
このようにして、本発明のノイズ低減装置１によれば、ノイズを含む入力画像から、その画像のデノイズされた画像を得ることができる。
【００４５】
次に、本発明に係るノイズ低減装置１のノイズ低減処理部２０の実施形態について説明する。ノイズ低減処理部２０は、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【００４６】
図３は、本発明に係るノイズ低減装置１のノイズ低減処理部２０の実施形態を示すブロック構成図である。
【００４７】
図３に示されたように、ノイズ低減処理部２０は、平均ベクトル算出部２１０と、分散共分散行列算出部２２０と、固有値固有ベクトル算出部２３０と、基底行列生成部２４０と、画像パッチ射影部２５０とから構成され、画像パッチ集合に属する全ての画像パッチに対し、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【００４８】
既述したように、一つの画像パッチは、一つの列ベクトルに対応する。よって、画像パッチ集合を列ベクトルの集合とすることが可能になり、画像パッチ集合を意味する列ベクトル並べた行列
は、下記数１によって表すことができる。以下、この行列
を単に「画像パッチ集合」とも言う。
【００４９】
【数１】
ここで、
は入力画像におけるｋ番目の画像パッチに対応する列ベクトルを表す。以下、この列ベクトル
を単に「ｋ番目の画像パッチ」、又は
とも言う。また、Ｎは画像パッチ集合の要素数を表す。つまり、Ｎは入力画像から得られた複数の画像パッチの個数である。
【００５０】
図３に示すように、ノイズ低減処理部２０では、まず、平均ベクトル算出部２１０が、画像パッチ集合に属する全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する。
【００５１】
ここで、説明の便宜上、ｋ番目の画像パッチに対し、平均ベクトル算出部２１０にて行われる平均ベクトル算出処理だけを説明する。
【００５２】
平均ベクトル算出部２１０では、下記数２及び数３に基づき、平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出する。
【００５３】
【数２】
【００５４】
【数３】
ここで、
は、画像パッチ
と画像パッチ
の間の類似度に基づく重み関数である。画像パッチ
と画像パッチ
の間の類似度が大きければ、重み
が大きくなるように設定される。即ち、ノイズ低減装置１で使用される重み
は、入力画像の画像パッチ間の類似度に基づいて設定されている。重み関数
の具体例については、後述する。
【００５５】
次に、ノイズ低減処理部２０では、分散共分散行列算出部２２０が、画像パッチ集合、及び平均ベクトル算出部２１０で算出された平均ベクトルに基づき、画像パッチ集合に属する全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する。
【００５６】
分散共分散行列算出部２２０では、下記数４に基づき、分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する。
【００５７】
【数４】
ここで、重み関数
の具体実施例について説明する。
【００５８】
第１例として、下記数５で表す関数を重み関数
として用いることができる。
【００５９】
【数５】
ただし、
は類似度評価における標準偏差を表す。また、
はユークリッドノルムを表す。
【００６０】
また、第２例として、下記数６で表す二値出力の重み関数を重み関数
として用いることができる。
【００６１】
【数６】
ただし、τは所定の定数である。
【００６２】
更に、第３例として、下記数７で表す重み関数を重み関数
として用いることができる。
【００６３】
【数７】
次に、ノイズ低減処理部２０では、固有値固有ベクトル算出部２３０が、固有値固有ベクトル算出処理を行うことにより、分散共分散行列算出部２２０で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する。
【００６４】
次に、ノイズ低減処理部２０では、基底行列生成部２４０が、固有値固有ベクトル算出部２３０で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する。
【００６５】
ここで、基底行列生成部２４０の２つの実施形態について説明する。
【００６６】
基底行列生成部２４０の第１実施形態では、下記数８に基づき、基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列
を生成する。
【００６７】
【数８】
ただし、
は画像パッチ
に対応する分散共分散行列
に大きな固有値に対応する固有ベクトルである。また、Ｌは予め設定された所定のパラメータであり、基底の個数を表す。
【００６８】
また、基底行列生成部２４０の第２実施形態では、固有値の値を利用して、上記数８に基づき、基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成する。
【００６９】
ただし、基底行列生成部２４０の第２実施形態では、
を満足するとする。ここで、
は固有ベクトル
に対応する固有値である。また、ρは予め設定された所定のパラメータである。
【００７０】
なお、基底行列生成部２４０の第２実施形態では、パラメータρを予め設定せず、
に基づき、パラメータρを設定するようにしても良い。ここで、
は入力画像のノイズの分散を表す。また、ｒは予め設定された所定のパラメータである。一例として、例えば、ｒの値を５と設定する。つまり、パラメータρを入力画像のノイズの分散の定数倍と設定するようにしても良い。
【００７１】
次に、ノイズ低減処理部２０では、画像パッチ射影部２５０が、画像パッチ集合に属する全ての画像パッチに対し、各画像パッチ、平均ベクトル算出部２１０で算出された各画像パッチに対応する平均ベクトル、及び、基底行列生成部２４０で生成された各画像パッチに対応する基底行列に基づき、画像パッチの射影を算出する画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチ（即ち、射影された画像パッチ）を生成する。
【００７２】
画像パッチ射影部２５０では、画像パッチ
平均ベクトル
及び基底行列
を用いて、下記数９に基づき、画像パッチ射影処理を行うことにより、画像パッチ
のデノイズされた画像パッチ
を生成する。
【００７３】
【数９】
上述のようにして、本発明のノイズ低減装置１のノイズ低減処理部２０が、画像パッチ集合に属する全ての画像パッチのデノイズされた画像パッチを生成することができる。
【００７４】
上記から分かるように、本発明では、画像パッチに対して行われるノイズ低減処理は、画像パッチをアフィン空間に射影する、画像パッチ射影処理を意味し、アフィン空間に射影された画像パッチは、デノイズされた画像パッチとなる。
【００７５】
つまり、画像パッチ
のデノイズされた画像パッチ
は、画像パッチ
をアフィン空間に射影することで得られた、アフィン空間に射影された画像パッチである。
【００７６】
ただし、平均ベクトル算出部２１０で算出された平均ベクトル
は、画像パッチ
に対応するアフィン空間の原点となる。また、基底行列生成部２４０で生成された基底行列
は、画像パッチ
に対応するアフィン空間の基底行列となる。
【００７７】
要するに、ノイズ低減処理部２０にて行われる「画像パッチ毎に、異なるパラメータによるノイズ低減処理」とは、各画像パッチを異なるアフィン空間（異なる原点と基底行列を持つアフィン空間）に射影することである。
【００７８】
図４は、本発明に係るノイズ低減装置の第２実施形態（以下、単に「ノイズ低減装置２」とも言う。）を示すブロック構成図である。以下、図４を用いて、ノイズ低減装置２を詳細に説明する。
【００７９】
図４に示すように、ノイズ低減装置２は、ノイズを含む入力画像から、デノイズされた画像を得るための装置であり、入力画像に対し、フィルタリング処理を行うフィルタリング処理部４０と、入力画像から画像パッチ集合を生成するとともに、フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成する画像パッチ集合生成部１１と、画像パッチ集合及びフィルタリングされた画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成するノイズ低減処理部２１と、ノイズ低減処理部２１で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を得る画像合成部３０とから構成される。
【００８０】
図４に示すように、ノイズ低減装置２では、まず、フィルタリング処理部４０が、入力画像に対し、所定のフィルタを用いてフィルタリング処理を行うことにより、フィルタリングされた画像を生成する。
【００８１】
ここで、所定のフィルタとして、例えば、ローパスフィルタや、バイラテラルフィルタを用いることができる。
【００８２】
また、離散コサイン変換（ＤＣＴ）を用いて入力画像を周波数領域へ変換し、変換された周波数領域で閾値処理を施す処理を、フィルタリング処理部４０にて行われるフィルタリング処理とすることもできる。
【００８３】
次に、ノイズ低減装置２では、図４に示すように、画像パッチ集合生成部１１が、入力画像から当該入力画像の画像パッチ集合を生成するとともに、フィルタリング処理部４０で得られたフィルタリングされた画像からフィルタリングされた画像パッチ集合を生成する。
【００８４】
具体的に、画像パッチ集合生成部１１は、入力画像に対して、画像パッチ集合生成処理を行うことにより、画像パッチ集合を生成するとともに、フィルタリングされた画像に対して、フィルタリングされた画像パッチ集合生成処理を行うことにより、フィルタリングされた画像パッチ集合を生成する。
【００８５】
ここで、画像パッチ集合生成処理は、前述したノイズ低減装置１の画像パッチ集合生成部１０にて行われる画像パッチ集合生成処理と同じであるため、その詳細な説明は省略する。
【００８６】
また、フィルタリングされた画像パッチ集合生成処理は、基本的に、画像パッチ集合生成処理と同じであるため、その詳細な説明は省略する。
【００８７】
ただし、画像パッチ集合生成処理は、入力画像に対して、その画素位置をずらしながら、複数の画像パッチを生成する処理であるのに対し、フィルタリングされた画像パッチ集合生成処理は、フィルタリングされた画像に対して、その画素位置をずらしながら、複数のフィルタリングされた画像パッチを生成する処理である。
【００８８】
次に、ノイズ低減装置２では、図４に示すように、ノイズ低減処理部２１が、画像パッチ集合生成部１１で生成された画像パッチ集合及びフィルタリングされた画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【００８９】
そして、ノイズ低減装置２では、図４に示すように、画像合成部３０が、ノイズ低減処理部２１で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成する。
【００９０】
ここで、画像合成処理は、前述したノイズ低減装置１の画像合成部３０にて行われる画像合成処理と同じであるため、その詳細な説明は省略する。
【００９１】
このようにして、本発明のノイズ低減装置２によれば、ノイズを含む入力画像から、その画像のデノイズされた画像を得ることができる。
【００９２】
次に、本発明に係るノイズ低減装置２のノイズ低減処理部２１の実施形態について説明する。ノイズ低減処理部２１は、画像パッチ集合及びフィルタリングされた画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【００９３】
図５は、本発明に係るノイズ低減装置２のノイズ低減処理部２１の実施形態を示すブロック構成図である。
【００９４】
図５に示されたように、ノイズ低減処理部２１は、平均ベクトル算出部２１１と、分散共分散行列算出部２２１と、固有値固有ベクトル算出部２３０と、基底行列生成部２４０と、画像パッチ射影部２５０とから構成され、画像パッチ集合及びフィルタリングされた画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対し、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【００９５】
図５に示すように、ノイズ低減処理部２１では、まず、平均ベクトル算出部２１１が、画像パッチ集合及びフィルタリングされた画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する。
【００９６】
ここで、説明の便宜上、ｋ番目の画像パッチに対し、平均ベクトル算出部２１１にて行われる平均ベクトル算出処理だけを説明する。
【００９７】
平均ベクトル算出部２１１では、下記数１０及び数１１に基づき、平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出する。
【００９８】
【数１０】
【００９９】
【数１１】
ここで、
は、画像パッチ集合生成部１１で生成されたフィルタリングされた画像パッチ集合におけるｋ番目のフィルタリングされた画像パッチを表し、以下、単に
とも言う。
はフィルタリングされた画像パッチ
とフィルタリングされた画像パッチ
の間の類似度に基づく重み関数である。フィルタリングされた画像パッチ
とフィルタリングされた画像パッチ
の間の類似度が大きければ、重み
が大きくなるように設定される。即ち、重み関数
は、重み関数
と同じように設定すれば良いので、その具体例の説明は省略する。ただし、ノイズ低減装置２で使用される重み
は、フィルタリングされた画像パッチ間の類似度に基づいて設定されている。
【０１００】
次に、ノイズ低減処理部２１では、分散共分散行列算出部２２１が、画像パッチ集合、フィルタリングされた画像パッチ集合、及び平均ベクトル算出部２１１で算出された平均ベクトルに基づき、画像パッチ集合に属する全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する。
【０１０１】
分散共分散行列算出部２２１では、下記数１２に基づき、分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する。
【０１０２】
【数１２】
次に、ノイズ低減処理部２１では、固有値固有ベクトル算出部２３０が、固有値固有ベクトル算出処理を行うことにより、分散共分散行列算出部２２１で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する。
【０１０３】
次に、ノイズ低減処理部２１では、基底行列生成部２４０が、固有値固有ベクトル算出部２３０で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する。
【０１０４】
ここで、ノイズ低減処理部２１の基底行列生成部２４０にて行われる処理は、前述したノイズ低減処理部２０の基底行列生成部２４０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１０５】
次に、ノイズ低減処理部２１では、画像パッチ射影部２５０が、画像パッチ集合に属する全ての画像パッチに対し、各画像パッチ、平均ベクトル算出部２１１で算出された各画像パッチに対応する平均ベクトル、及び、基底行列生成部２４０で生成された各画像パッチに対応する基底行列に基づき、画像パッチの射影を算出する画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチ（即ち、射影された画像パッチ）を生成する。
【０１０６】
ここで、ノイズ低減処理部２１の画像パッチ射影部２５０にて行われる処理は、前述したノイズ低減処理部２０の画像パッチ射影部２５０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１０７】
上述のようにして、本発明のノイズ低減装置２のノイズ低減処理部２１が、画像パッチ集合及びフィルタリングされた画像パッチ集合に基づいて、画像パッチ集合に属する全ての画像パッチのデノイズされた画像パッチを生成することができる。
【０１０８】
図６は、本発明に係るノイズ低減装置の第３実施形態（以下、単に「ノイズ低減装置３」とも言う。）を示すブロック構成図である。以下、図６を用いて、ノイズ低減装置３を詳細に説明する。
【０１０９】
図６に示すように、ノイズ低減装置３は、ノイズを含む入力画像から、デノイズされた画像を得るための装置であり、予め撮影された一枚以上の画像を格納する画像データベース５０と、入力画像から画像パッチ集合を生成するとともに、画像データベース５０に格納されている画像からデータベース画像パッチ集合を生成する画像パッチ集合生成部１２と、画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、ノイズ低減処理を行うことにより、デノイズされた画像パッチを生成するノイズ低減処理部２２と、ノイズ低減処理部２２で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を得る画像合成部３０とから構成される。
【０１１０】
図６に示すように、ノイズ低減装置３では、まず、画像パッチ集合生成部１２が、入力画像から当該入力画像の画像パッチ集合を生成するとともに、画像データベース５０に格納されている画像からデータベース画像パッチ集合を生成する。
【０１１１】
具体的に、画像パッチ集合生成部１２は、入力画像に対して、画像パッチ集合生成処理を行うことにより、画像パッチ集合を生成するとともに、画像データベース５０に格納されている画像に対して、データベース画像パッチ集合生成処理を行うことにより、データベース画像パッチ集合を生成する。
【０１１２】
ここで、画像パッチ集合生成処理は、前述したノイズ低減装置１の画像パッチ集合生成部１０にて行われる画像パッチ集合生成処理と同じであるため、その詳細な説明は省略する。
【０１１３】
また、データベース画像パッチ集合生成処理は、基本的に、画像パッチ集合生成処理と同じであるため、その詳細な説明は省略する。
【０１１４】
ただし、画像パッチ集合生成処理は、入力画像に対して、その画素位置をずらしながら、複数の画像パッチを生成する処理であるのに対し、データベース画像パッチ集合生成処理は、画像データベース５０に格納されている画像に対して、その画素位置をずらしながら、複数のデータベース画像パッチを生成する処理である。
【０１１５】
次に、ノイズ低減装置３では、図６に示すように、ノイズ低減処理部２２が、画像パッチ集合生成部１２で生成された画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【０１１６】
そして、ノイズ低減装置３では、図６に示すように、画像合成部３０が、ノイズ低減処理部２２で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成する。
【０１１７】
ここで、画像合成処理は、前述したノイズ低減装置１の画像合成部３０にて行われる画像合成処理と同じであるため、その詳細な説明は省略する。
【０１１８】
このようにして、本発明のノイズ低減装置３によれば、ノイズを含む入力画像から、その画像のデノイズされた画像を得ることができる。
【０１１９】
次に、本発明に係るノイズ低減装置３のノイズ低減処理部２２の実施形態について説明する。ノイズ低減処理部２２は、画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【０１２０】
図７は、本発明に係るノイズ低減装置３のノイズ低減処理部２２の実施形態を示すブロック構成図である。
【０１２１】
図７に示されたように、ノイズ低減処理部２２は、平均ベクトル算出部２１２と、分散共分散行列算出部２２２と、固有値固有ベクトル算出部２３０と、基底行列生成部２４０と、画像パッチ射影部２５０とから構成され、画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対し、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【０１２２】
図７に示すように、ノイズ低減処理部２２では、まず、平均ベクトル算出部２１２が、画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する。
【０１２３】
ここで、説明の便宜上、ｋ番目の画像パッチに対し、平均ベクトル算出部２１２にて行われる平均ベクトル算出処理だけを説明する。
【０１２４】
平均ベクトル算出部２１２では、下記数１３及び数１４に基づき、平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出する。
【０１２５】
【数１３】
【０１２６】
【数１４】
ここで、
は、画像パッチ集合生成部１２で生成されたデータベース画像パッチ集合におけるｉ番目のデータベース画像パッチを表し、以下、単に
とも言う。Ｍはデータベース画像パッチ集合の要素数である。
はデータベース画像パッチ
と画像パッチ
の間の類似度に基づく重み関数である。データベース画像パッチ
と画像パッチ
の間の類似度が大きければ、重み
が大きくなるように設定される。即ち、重み関数
は、重み関数
と同じように設定すれば良いので、その具体例の説明は省略する。ただし、ノイズ低減装置３で使用される重み
は、データベース画像パッチと画像パッチの間の類似度に基づいて設定されている。
【０１２７】
次に、ノイズ低減処理部２２では、分散共分散行列算出部２２２が、画像パッチ集合、データベース画像パッチ集合、及び平均ベクトル算出部２１２で算出された平均ベクトルに基づき、画像パッチ集合に属する全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する。
【０１２８】
分散共分散行列算出部２２２では、下記数１５に基づき、分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する。
【０１２９】
【数１５】
次に、ノイズ低減処理部２２では、固有値固有ベクトル算出部２３０が、固有値固有ベクトル算出処理を行うことにより、分散共分散行列算出部２２２で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する。
【０１３０】
次に、ノイズ低減処理部２２では、基底行列生成部２４０が、固有値固有ベクトル算出部２３０で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する。
【０１３１】
ここで、ノイズ低減処理部２２の基底行列生成部２４０にて行われる処理は、前述したノイズ低減処理部２０の基底行列生成部２４０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１３２】
次に、ノイズ低減処理部２２では、画像パッチ射影部２５０が、画像パッチ集合に属する全ての画像パッチに対し、各画像パッチ、平均ベクトル算出部２１２で算出された各画像パッチに対応する平均ベクトル、及び、基底行列生成部２４０で生成された各画像パッチに対応する基底行列に基づき、画像パッチの射影を算出する画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する。
【０１３３】
ここで、ノイズ低減処理部２２の画像パッチ射影部２５０にて行われる処理は、前述したノイズ低減処理部２０の画像パッチ射影部２５０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１３４】
上述のようにして、本発明のノイズ低減装置３のノイズ低減処理部２２が、画像パッチ集合及びデータベース画像パッチ集合に基づいて、画像パッチ集合に属する全ての画像パッチのデノイズされた画像パッチを生成することができる。
【０１３５】
図８は、本発明に係るノイズ低減装置の第４実施形態（以下、単に「ノイズ低減装置４」とも言う。）を示すブロック構成図である。以下、図８を用いて、ノイズ低減装置４を詳細に説明する。
【０１３６】
図８に示すように、ノイズ低減装置４は、ノイズを含む入力画像から、デノイズされた画像を得るための装置であり、入力画像に対し、フィルタリング処理を行うフィルタリング処理部４０と、予め撮影された一枚以上の画像を格納する画像データベース５０と、フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成するとともに、画像データベース５０に格納されている画像からデータベース画像パッチ集合を生成し、さらに、入力画像から画像パッチ集合を生成する画像パッチ集合生成部１３と、画像パッチ集合、フィルタリングされた画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成するノイズ低減処理部２３と、ノイズ低減処理部２３で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を得る画像合成部３０とから構成される。
【０１３７】
図８に示すように、ノイズ低減装置４では、まず、フィルタリング処理部４０が、入力画像に対し、所定のフィルタを用いてフィルタリング処理を行うことにより、フィルタリングされた画像を生成する。
【０１３８】
ここで、ノイズ低減装置４のフィルタリング処理部４０にて行われる処理は、前述したノイズ低減装置２のフィルタリング処理部４０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１３９】
次に、ノイズ低減装置４では、図８に示すように、画像パッチ集合生成部１３が、フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成するとともに、画像データベース５０に格納されている画像からデータベース画像パッチ集合を生成し、さらに、入力画像から画像パッチ集合を生成する。
【０１４０】
具体的に、画像パッチ集合生成部１３は、フィルタリングされた画像に対して、フィルタリングされた画像パッチ集合生成処理を行うことにより、フィルタリングされた画像パッチ集合を生成するとともに、画像データベース５０に格納されている画像に対して、データベース画像パッチ集合生成処理を行うことにより、データベース画像パッチ集合を生成し、さらに、入力画像に対して、画像パッチ集合生成処理を行うことにより、画像パッチ集合を生成する。
【０１４１】
ここで、フィルタリングされた画像パッチ集合生成処理は、前述したノイズ低減装置２の画像パッチ集合生成部１１にて行われるフィルタリングされた画像パッチ集合生成処理と同じであるため、その詳細な説明は省略する。
【０１４２】
また、データベース画像パッチ集合生成処理は、前述したノイズ低減装置３の画像パッチ集合生成部１２にて行われるデータベース画像パッチ集合生成処理と同じであるため、その詳細な説明は省略する。
【０１４３】
また、画像パッチ集合生成処理は、前述したノイズ低減装置１の画像パッチ集合生成部１０にて行われる画像パッチ集合生成処理と同じであるため、その詳細な説明は省略する。
【０１４４】
次に、ノイズ低減装置４では、図８に示すように、ノイズ低減処理部２３が、画像パッチ集合生成部１３で生成された、画像パッチ集合、フィルタリングされた画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【０１４５】
そして、ノイズ低減装置４では、図８に示すように、画像合成部３０が、ノイズ低減処理部２３で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成する。
【０１４６】
ここで、画像合成処理は、前述したノイズ低減装置１の画像合成部３０にて行われる画像合成処理と同じであるため、その詳細な説明は省略する。
【０１４７】
このようにして、本発明のノイズ低減装置４によれば、ノイズを含む入力画像から、その画像のデノイズされた画像を得ることができる。
【０１４８】
次に、本発明に係るノイズ低減装置４のノイズ低減処理部２３の実施形態について説明する。ノイズ低減処理部２３は、画像パッチ集合、フィルタリングされた画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【０１４９】
図９は、本発明に係るノイズ低減装置４のノイズ低減処理部２３の実施形態を示すブロック構成図である。
【０１５０】
図９に示されたように、ノイズ低減処理部２３は、平均ベクトル算出部２１３と、分散共分散行列算出部２２３と、固有値固有ベクトル算出部２３０と、基底行列生成部２４０と、画像パッチ射影部２５１とから構成され、画像パッチ集合、フィルタリングされた画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対し、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【０１５１】
図９に示すように、ノイズ低減処理部２３では、まず、平均ベクトル算出部２１３が、画像パッチ集合、フィルタリングされた画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する。
【０１５２】
ここで、説明の便宜上、ｋ番目の画像パッチに対し、平均ベクトル算出部２１３にて行われる平均ベクトル算出処理だけを説明する。
【０１５３】
平均ベクトル算出部２１３では、下記数１６及び数１７に基づき、平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出する。
【０１５４】
【数１６】
【０１５５】
【数１７】
ここで、
はデータベース画像パッチ
とフィルタリングされた画像パッチ
の間の類似度に基づく重み関数である。データベース画像パッチ
とフィルタリングされた画像パッチ
の間の類似度が大きければ、重み
が大きくなるように設定される。即ち、重み関数
は、重み関数
と同じように設定すれば良いので、その具体例の説明は省略する。ただし、ノイズ低減装置４で使用される重み
は、データベース画像パッチとフィルタリングされた画像パッチの間の類似度に基づいて設定されている。
【０１５６】
次に、ノイズ低減処理部２３では、分散共分散行列算出部２２３が、画像パッチ集合、フィルタリングされた画像パッチ集合、データベース画像パッチ集合、及び平均ベクトル算出部２１３で算出された平均ベクトルに基づき、画像パッチ集合に属する全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する。
【０１５７】
分散共分散行列算出部２２３では、下記数１８に基づき、分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する。
【０１５８】
【数１８】
次に、ノイズ低減処理部２３では、固有値固有ベクトル算出部２３０が、固有値固有ベクトル算出処理を行うことにより、分散共分散行列算出部２２３で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する。
【０１５９】
次に、ノイズ低減処理部２３では、基底行列生成部２４０が、固有値固有ベクトル算出部２３０で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する。
【０１６０】
ここで、ノイズ低減処理部２３の基底行列生成部２４０にて行われる処理は、前述したノイズ低減処理部２０の基底行列生成部２４０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１６１】
次に、ノイズ低減処理部２３では、画像パッチ射影部２５０が、画像パッチ集合に属する全ての画像パッチに対し、各画像パッチ、平均ベクトル算出部２１３で算出された各画像パッチに対応する平均ベクトル、及び、基底行列生成部２４０で生成された各画像パッチに対応する基底行列に基づき、画像パッチの射影を算出する画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する。
【０１６２】
ここで、ノイズ低減処理部２３の画像パッチ射影部２５０にて行われる処理は、前述したノイズ低減処理部２０の画像パッチ射影部２５０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１６３】
上述のようにして、本発明のノイズ低減装置４のノイズ低減処理部２３が、画像パッチ集合、フィルタリングされた画像パッチ集合及びデータベース画像パッチ集合に基づいて、画像パッチ集合に属する全ての画像パッチのデノイズされた画像パッチを生成することができる。
【０１６４】
図１０は、本発明に係るノイズ低減装置の第５実施形態（以下、単に「ノイズ低減装置５」とも言う。）を示すブロック構成図である。また、図１１はノイズ低減装置５にて行われる処理の流れを示すフロー図である。以下、図１０及び図１１を用いて、ノイズ低減装置５を詳細に説明する。
【０１６５】
ノイズ低減装置５は、本発明に係るノイズ低減装置の第２実施形態と本発明に係るノイズ低減装置の第４実施形態を組み合わせた、繰り返し処理を有する実施形態である。
【０１６６】
図１０に示すように、ノイズ低減装置５は、ノイズを含む入力画像から、デノイズされた画像を得るための装置であり、フィルタリング処理部４０と、画像データベース５１と、画像パッチ集合生成部１４と、ノイズ低減処理部２４と、画像合成部３１とから構成される。
【０１６７】
ここで、ノイズを含む入力画像からデノイズされた画像を得るために、図１０と図１１を参照しながら、ノイズ低減装置５にて行われる処理の流れについて説明する。
【０１６８】
ノイズ低減装置５では、まず、繰り返し処理の回数ｋの値を０と設定する（図１１のステップＳ１００を参照）。
【０１６９】
次に、ノイズ低減装置５では、フィルタリング処理部４０が、入力画像に対し、所定のフィルタを用いてフィルタリング処理を行うことにより、フィルタリングされた画像を生成する（図１１のステップＳ１１０を参照）。
【０１７０】
ここで、ノイズ低減装置５のフィルタリング処理部４０にて行われる処理は、前述したノイズ低減装置２のフィルタリング処理部４０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１７１】
次に、ノイズ低減装置５では、画像パッチ集合生成部１４が、入力画像から当該入力画像の画像パッチ集合を生成するとともに、フィルタリング処理部４０で得られたフィルタリングされた画像からフィルタリングされた画像パッチ集合を生成する（図１１のステップＳ１２０を参照）。
【０１７２】
ここで、ノイズ低減装置５の画像パッチ集合生成部１４にて行われる処理は、前述したノイズ低減装置２の画像パッチ集合生成部１１にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１７３】
次に、ノイズ低減装置５では、ノイズ低減処理部２４が、画像パッチ集合生成部１４で生成された画像パッチ集合及びフィルタリングされた画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する（図１１のステップＳ１３０を参照）。
【０１７４】
ここで、ノイズ低減装置５のノイズ低減処理部２４にて行われる処理は、前述したノイズ低減装置２のノイズ低減処理部２１にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１７５】
次に、ノイズ低減装置５では、画像合成部３１が、ノイズ低減処理部２４で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成する（図１１のステップＳ１４０を参照）。
【０１７６】
ここで、ノイズ低減装置５の画像合成部３１にて行われる処理は、前述したノイズ低減装置２の画像合成部３０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１７７】
次に、ノイズ低減装置５では、デノイズされた画像を、画像データベース５１の画像として、画像データベース５１に格納する（図１１のステップＳ１５０を参照）。
【０１７８】
次に、ノイズ低減装置５では、画像パッチ集合生成部１４が、画像データベース５１に格納されている画像に対して、データベース画像パッチ集合生成処理を行うことにより、データベース画像パッチ集合を生成する（図１１のステップＳ１６０を参照）。
【０１７９】
ここで、ノイズ低減装置５の画像パッチ集合生成部１４にて行われるデータベース画像パッチ集合生成処理は、前述したノイズ低減装置４の画像パッチ集合生成部１３にて行われるデータベース画像パッチ集合生成処理と同じであるため、その詳細な説明は省略する。
【０１８０】
次に、ノイズ低減装置５では、ノイズ低減処理部２４が、画像パッチ集合生成部１４で生成された、画像パッチ集合、フィルタリングされた画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する（図１１のステップＳ１７０を参照）。
【０１８１】
ここで、ノイズ低減装置５のノイズ低減処理部２４にて行われる処理は、前述したノイズ低減装置４のノイズ低減処理部２３にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１８２】
次に、ノイズ低減装置５では、画像合成部３１が、ノイズ低減処理部２４で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成する（図１１のステップＳ１８０を参照）。
【０１８３】
ここで、ノイズ低減装置５の画像合成部３１にて行われる処理は、前述したノイズ低減装置４の画像合成部３０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１８４】
次に、ノイズ低減装置５では、繰り返し処理の回数ｋの値を１増加させるようにする（図１１のステップＳ１９０を参照）。
【０１８５】
次に、ノイズ低減装置５では、繰り返し処理終了条件（即ち、ｋ＞所定の値）を満たすか否かを判定し（図１１のステップＳ２００を参照）、繰り返し処理終了条件（ｋ＞所定の値）を満たしたと判定した場合に、デノイズされた画像をノイズ低減装置５が生成したデノイズされた画像として外部に出力する（図１１のステップＳ２１０を参照）。これで、ノイズ低減装置５での処理が終了する。
【０１８６】
一方、ステップＳ２００において、ノイズ低減装置５が、繰り返し処理終了条件（ｋ＞所定の値）を満たさないと判定した場合に、ステップＳ１５０に戻り、ノイズ低減装置５では、ステップＳ１５０での処理、ステップＳ１６０での処理、ステップＳ１７０での処理、ステップＳ１８０での処理、ステップＳ１９０での処理、及びステップＳ２００での処理を繰り返すようにする。
【０１８７】
このようにして、本発明のノイズ低減装置５によれば、ノイズを含む入力画像から、その画像のデノイズされた画像を得ることができる。
【０１８８】
なお、本発明に係るノイズ低減装置は、コンピュータシステムを利用し、ソフトウェア（コンピュータプログラム）により実装されることができ、そして、ＡＳＩＣ（Application Specific Integrated Circuit）、ＧＰＵ（Graphics Processing Unit）やＦＰＧＡ（Field Programmable Gate Array）などのハードウェアにより実装されることも勿論できる。
【符号の説明】
【０１８９】
ノイズ低減装置 １,２,３,４,５
画像パッチ集合生成部 １０,１１,１２,１３,１４
ノイズ低減処理部 ２０,２１,２２,２３,２４
画像合成部 ３０,３１
フィルタリング処理部 ４０
画像データベース ５０,５１
平均ベクトル算出部 ２１０,２１１,２１２,２１３
分散共分散行列算出部 ２２０,２２１,２２２,２２３
固有値固有ベクトル算出部 ２３０
基底行列生成部 ２４０
画像パッチ射影部 ２５０
【技術分野】
【０００１】
本発明は、デジタル画像処理技術に関し、特に、デジタル画像のノイズを低減し、デノイズされた画像を得るためのノイズ低減装置、ノイズ低減方法及びプログラムに関する。
【背景技術】
【０００２】
近年、画像取得装置の一つとして、デジタルカメラは急速に普及している。高品位な画質を得るために、デジタルカメラで撮影されたデジタル画像に対し、様々なデジタル画像処理が施される。数多くのデジタル画像処理（以下、単に、「画像処理」と言う。）の中で、ノイズ低減処理が画質に大きな影響を与えるため、大変重要な画像処理の一つであることが広く知られている。
【０００３】
従来では、例えば、非特許文献１、非特許文献２、及び特許文献１に開示されたように、デジタル画像（以下、単に、「画像」とも言う。）に含まれるノイズを低減するための画像ノイズ低減技術が数多く提案されている。
【０００４】
特許文献１には、ノイズ低減効果の低い部分の出現を抑え、画像全体として高いノイズ低減効果を得ることができる画像信号処理装置が開示されている。
【０００５】
特許文献１の画像信号処理装置は、画像信号に対して第１の平滑化処理を行うことで、ノイズ低減された第１の平滑化信号を生成する第１の平滑化信号生成部と、画像信号に対して第１の平滑化処理と異なる第２の平滑化処理を行うことで、ノイズ低減された第２の平滑化信号を生成する第２の平滑化信号生成部と、画像信号、第１の平滑化信号、及び第２の平滑化信号のいずれか１つの値に応じた混合割合に従って第１の平滑化信号と第２の平滑化信号とを混合処理する混合処理部とを備えることを主な特徴としている。
【０００６】
つまり、特許文献１では、画像に対してノイズ低減効果の異なる２つの平滑化処理を行い、それら２つの平滑化処理により生成された２つの平滑化信号を所定の混合割合で混合処理することにより、当該画像のノイズを低減するようにしている。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２０１１−８７０８７号公報
【非特許文献】
【０００８】
【非特許文献１】ジェー. ポーティラ（J. Portilla）、ブィ. ストレラ（V. Strela）、エム. ジェー. ウェインライト（M. J. Wainwright）、イー. ピー. シモンセリ（E. P. Simoncelli）共著,「イメージ デノイジング ユージング スケール ミックスチャー オフ ガウシアン イン ザ ウェーブレット ドメイン （Image Denoising using Scale Mixture of Gaussians in the Wavelet Domain）」,ＩＥＥＥ トランスアクションズ オン イメージ プロセシング（IEEE Transactions on Image Processing）,第12巻,第11号,p.1338-1351,2003年
【非特許文献２】エイ. ブアデス（A. Buades）、ビー. コル（B. Coll）、ジェー. エム. モレル（J.M. Morel）共著,「ア ノンローカル アルゴリズム フォー イメージ デノイジング （A non-local algorithm for image denoising）」,ＩＥＥＥ コンピュータ ソサイエティー カンファレンス オン コンピュータ ビジョン アンド パターン レコグニション （IEEE Computer Society Conference on Computer Vision and Pattern Recognition）,第2巻,p.60-65,2005年
【発明の概要】
【発明が解決しようとする課題】
【０００９】
ところで、自然画像は、その画素位置（又は画像領域）によりテクスチャやエッジが大きく異なるという性質を有している。
【００１０】
しかし、従来のノイズ低減技術（例えば、非特許文献１及び非特許文献２を参照）の殆どは、テクスチャ領域・エッジ領域を問わず、画像全体に対して均質的にノイズ低減処理を行うものであり、つまり、画像領域毎（又は画素位置毎）に異なるノイズ低減処理を行っていない。
【００１１】
また、特許文献１では、２つの平滑化処理により生成された２つの平滑化信号を混合処理する際の混合割合は、画像の箇所毎（画素位置毎）に異なる値となっているものの、基本的に、特許文献１において、画像のノイズを低減するために、２種類の平滑化処理しか行っておらず、画像の箇所毎（画素位置毎）に異なる処理とはなっていない。
【００１２】
本発明は、上述のような事情よりなされたものであり、本発明の目的は、ノイズを含んでいる一枚の画像に対し、位置毎に異なるノイズ低減処理を簡単に行えるようにした、ノイズ低減装置、ノイズ低減方法及びプログラムを提供することにある。
【課題を解決するための手段】
【００１３】
本発明は、ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置に関し、本発明の上記目的は、前記画像から画像パッチ集合を生成する、画像パッチ集合生成部と、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部と、前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成部とを備え、前記ノイズ低減処理部は、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出部と、前記画像パッチ集合、及び前記平均ベクトル算出部で算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出部と、固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出部で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出部と、前記固有値固有ベクトル算出部で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成部と、前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出部で算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成部で生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影部とを備えることによって効果的に達成される。
【００１４】
また、本発明の上記目的は、ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置であって、前記画像に対し、フィルタリング処理を行うことにより、フィルタリングされた画像を生成する、フィルタリング処理部と、前記画像から画像パッチ集合を生成するとともに、前記フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成する、画像パッチ集合生成部と、前記画像パッチ集合及び前記フィルタリングされた画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部と、前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成部とを備え、前記ノイズ低減処理部は、前記画像パッチ集合及び前記フィルタリングされた画像パッチ集合に基づき、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出部と、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記平均ベクトル算出部で算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出部と、固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出部で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出部と、前記固有値固有ベクトル算出部で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成部と、前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出部で算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成部で生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影部とを備えることによって効果的に達成される。
【００１５】
また、本発明の上記目的は、ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置であって、予め撮影された一枚以上の画像を格納する画像データベースと、前記画像から画像パッチ集合を生成するとともに、前記画像データベースに格納されている画像からデータベース画像パッチ集合を生成する、画像パッチ集合生成部と、前記画像パッチ集合及び前記データベース画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部と、前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成部とを備え、前記ノイズ低減処理部は、前記画像パッチ集合及び前記データベース画像パッチ集合に基づき、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出部と、前記画像パッチ集合、前記データベース画像パッチ集合、及び前記平均ベクトル算出部で算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出部と、固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出部で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出部と、前記固有値固有ベクトル算出部で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成部と、前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出部で算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成部で生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影部とを備えることによって効果的に達成される。
【００１６】
また、本発明の上記目的は、ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置であって、前記画像に対し、フィルタリング処理を行うことにより、フィルタリングされた画像を生成する、フィルタリング処理部と、予め撮影された一枚以上の画像を格納する画像データベースと、前記フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成するとともに、前記画像データベースに格納されている画像からデータベース画像パッチ集合を生成し、さらに、前記画像から画像パッチ集合を生成する、画像パッチ集合生成部と、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部と、前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成部とを備え、前記ノイズ低減処理部は、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づき、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出部と、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、前記データベース画像パッチ集合、及び前記平均ベクトル算出部で算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出部と、固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出部で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出部と、前記固有値固有ベクトル算出部で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成部と、前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出部で算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成部で生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影部とを備えることによって効果的に達成される。
【００１７】
また、本発明の上記目的は、ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置であって、フィルタリング処理部と、画像データベースと、画像パッチ集合生成部と、ノイズ低減処理部と、画像合成部とを備え、前記ノイズ低減装置では、前記フィルタリング処理部が、前記画像に対し、フィルタリング処理を行うことにより、フィルタリングされた画像を生成し、前記画像パッチ集合生成部が、前記画像から画像パッチ集合を生成するとともに、前記フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成し、前記ノイズ低減処理部が、前記画像パッチ集合及び前記フィルタリングされた画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成し、前記画像合成部が、前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成し、前記デノイズされた画像が、前記画像データベースの画像として、前記画像データベースに格納され、前記画像パッチ集合生成部が、前記画像データベースに格納されている画像に対して、データベース画像パッチ集合生成処理を行うことにより、データベース画像パッチ集合を生成し、前記ノイズ低減処理部が、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づき、前記全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成し、前記画像合成部が、前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成し、繰り返し処理終了条件を満たしたと判定された場合に、前記デノイズされた画像を前記ノイズ低減装置が生成したデノイズされた画像とし、前記繰り返し処理終了条件を満たさないと判定された場合に、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づいて行われるノイズ低減処理を繰り返し行うことにより、デノイズされた画像パッチを生成し、生成したデノイズされた画像パッチを前記画像データベースの画像とすることによって効果的に達成される。
【発明の効果】
【００１８】
本発明に係るノイズ低減装置によれば、ノイズを含んでいる一枚の画像に対し、画像パッチ毎に異なるパラメータによるノイズ低減処理を行った上で、デノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を得るようにしているので、ノイズを含む画像の位置毎に異なるノイズ低減処理を簡単に行うことが可能になり、当該画像の領域（テクスチャ領域・エッジ領域）別に最適なノイズ低減効果を得ることができる。
【００１９】
本発明の発明者らは、異なるレベルのノイズを付加した１２枚のデジタル画像に対し、非特許文献１に記載された従来のノイズ低減方法を用いてノイズ低減を行った結果、非特許文献２に記載された従来のノイズ低減方法を用いてノイズ低減を行った結果、及び、本発明に係るノイズ低減方法を用いてノイズ低減を行った結果を比較することにより、本発明のノイズ低減効果を検証した。
【００２０】
表１は、異なるレベルのノイズが付加されたデジタル画像に対して、非特許文献１によるノイズ低減方法、非特許文献２によるノイズ低減方法、及び本発明によるノイズ低減方法を用いて、ノイズ低減を行って得られた結果画像（即ち、デノイズされた画像）のＰＳＮＲ(Peak Signal-to-Noise Ratio)を示している。
【００２１】
ＰＳＮＲは画質指標の一つであり、ＰＳＮＲの値が高いほど、結果画像（デノイズされた画像）が原画像（ノイズが付加されていないデジタル画像）に近いことを意味する。
【００２２】
表１の結果から、本発明に係るノイズ低減方法によるノイズ低減効果が一番高いことは、はっきりと確認された。
【００２３】
【表１】
【図面の簡単な説明】
【００２４】
【図１】本発明に係るノイズ低減装置の第１実施形態（ノイズ低減装置１）を示すブロック構成図である。
【図２】本発明に係るノイズ低減装置の画像パッチ集合生成部で行われる画像パッチ集合生成処理を説明するための模式図である。
【図３】本発明に係るノイズ低減装置１のノイズ低減処理部２０の実施形態を示すブロック構成図である。
【図４】本発明に係るノイズ低減装置の第２実施形態（ノイズ低減装置２）を示すブロック構成図である。
【図５】本発明に係るノイズ低減装置２のノイズ低減処理部２１の実施形態を示すブロック構成図である。
【図６】本発明に係るノイズ低減装置の第３実施形態（ノイズ低減装置３）を示すブロック構成図である。
【図７】本発明に係るノイズ低減装置３のノイズ低減処理部２２の実施形態を示すブロック構成図である。
【図８】本発明に係るノイズ低減装置の第４実施形態（ノイズ低減装置４）を示すブロック構成図である。
【図９】本発明に係るノイズ低減装置４のノイズ低減処理部２３の実施形態を示すブロック構成図である。
【図１０】本発明に係るノイズ低減装置の第５実施形態（ノイズ低減装置５）を示すブロック構成図である。
【図１１】本発明に係るノイズ低減装置５の処理流れを示すフロー図である。
【発明を実施するための形態】
【００２５】
本発明は、ノイズを含んでいる一枚の画像に対し、位置毎に異なるノイズ低減処理を簡単に行えるようにした、ノイズ低減装置、ノイズ低減方法及びプログラムに関する。
【００２６】
以下、本発明を実施するための形態について、図面を参照しながら詳細に説明する。
【００２７】
図１は、本発明に係るノイズ低減装置の第１実施形態（以下、単に「ノイズ低減装置１」とも言う。）を示すブロック構成図である。以下、図１を用いて、ノイズ低減装置１を詳細に説明する。
【００２８】
図１に示すように、ノイズ低減装置１は、ノイズを含む一枚のデジタル静止画像（以下、単に、「入力画像」又は「画像」とも言う。）から、デノイズされた画像を得るための装置であり、入力画像から画像パッチ集合を生成する画像パッチ集合生成部１０と、画像パッチ集合生成部１０で生成された画像パッチ集合に属する全ての画像パッチに対して、ノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部２０と、ノイズ低減処理部２０で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を得る、画像合成部３０とから構成される。
【００２９】
本発明では、デジタル静止画像であれば、複数チャネル画像・単一チャネル画像を問わず、本発明に係るノイズ低減技術を適用し、その画像のデノイズされた画像を簡単に得ることができる。
【００３０】
また、デジタル静止画像に限らず、動画像のフレームに、本発明に係るノイズ低減技術を適用し、そのフレームのデノイズされた画像（デノイズされたフレーム）を簡単に得ることが可能であり、本発明を動画像のノイズ低減にも適用することができると言える。
【００３１】
以下、説明上の便宜から、本発明のノイズ低減装置に入力されるデジタル画像は、例えば、グレースケール画像のような単一チャネル画像（即ち、単一のカラーチャネルを持つ画像）であることを前提とする。
【００３２】
また、本発明に係るノイズ低減技術を複数チャネル画像（複数のカラーチャネルを持つカラー画像）に適用する場合に、各カラーチャネルに対して、後述のような、単一チャネル画像を入力画像とする本発明に係るノイズ低減装置で行われる処理を繰り返せばよい。
【００３３】
図１に示すように、ノイズ低減装置１では、まず、画像パッチ集合生成部１０が、入力画像から、当該入力画像の画像パッチ集合を生成する。
【００３４】
具体的に、画像パッチ集合生成部１０は、入力画像から画像パッチ集合を得るために、入力画像に対して、画像パッチ集合生成処理を行うことにより、画像パッチ集合を生成する。画像パッチ集合生成処理は、図２に示すように、入力画像に対して、画素位置をずらしながら、複数の画像パッチを生成する処理である。
【００３５】
本発明では、画像パッチ集合生成処理において、画素位置のずらし量を任意に設定することができるが、ずらし量の一具体例として、１画素であり、つまり、画素位置を１画素ずらしながら、複数の画像パッチを生成する。
【００３６】
本発明では、画像パッチ集合生成処理により得られた複数の画像パッチを一つの集合（即ち、画像パッチ集合）としている。画像パッチは、その画像パッチに属する全ての画素の画素値を一列に並べることにより、列ベクトルに変換することができる。つまり、一つの画像パッチは、一つの列ベクトルに対応する。
【００３７】
以下、画像パッチの一具体例として、画像内のｍ画素×ｎ画素の矩形領域を画像パッチとする（図２を参照する。）。説明上の便宜から、以下の説明では、画像パッチを矩形領域（ｍ画素×ｎ画素）としているが、本発明において、画像パッチは矩形領域に限定されることがなく、任意の形状（例えば、多角形）を有する領域を画像パッチとすることができる。
【００３８】
次に、ノイズ低減装置１では、図１に示すように、ノイズ低減処理部２０が、画像パッチ集合生成部１０で生成された画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【００３９】
そして、ノイズ低減装置１では、図１に示すように、画像合成部３０が、ノイズ低減処理部２０で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成する。
【００４０】
ここで、画像合成部３０にて行われる画像合成処理とは、ノイズ低減された全ての画像パッチ（即ち、ノイズ低減処理部２０で得られた全てのデノイズされた画像パッチ）を画像の元の画素位置に戻す処理であり、つまり、デノイズされた画像パッチの画素値を画像の元の画素位置に対応するデノイズされた画像の画素値とする処理である。
【００４１】
ところで、図２にも示されるように、画像パッチ集合生成部１０では、画像パッチ集合生成処理により入力画像から複数の画像パッチを生成する際に、画素位置を所定個数の画素（例えば、１画素）ずらしながら、画像パッチ同士が重なりあうように、画像パッチが生成される。
【００４２】
このように重なりあって生成された各画像パッチに対し、それぞれノイズ低減処理を行うと、画像パッチが重なっている画素において、画素値が一致しないことになる。つまり、画素位置によって、複数のデノイズされた画像パッチに対応する画素があれば、一つのデノイズされた画素パッチに対応する画素もある。
【００４３】
そこで、本発明は、このような画像パッチが重なっている画素（即ち、複数のデノイズされた画像パッチに対応する画素）に対し、これら複数のデノイズされた画像パッチの画素値を平均し、その平均値を画像合成処理後の画素値（即ち、デノイズされた画像の画素値）とするように、画像合成処理を行う。なお、本発明では、複数のデノイズされた画像パッチの画素値を平均する際に、単純平均のかわりに、重み付き平均を用いても良い。
【００４４】
このようにして、本発明のノイズ低減装置１によれば、ノイズを含む入力画像から、その画像のデノイズされた画像を得ることができる。
【００４５】
次に、本発明に係るノイズ低減装置１のノイズ低減処理部２０の実施形態について説明する。ノイズ低減処理部２０は、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【００４６】
図３は、本発明に係るノイズ低減装置１のノイズ低減処理部２０の実施形態を示すブロック構成図である。
【００４７】
図３に示されたように、ノイズ低減処理部２０は、平均ベクトル算出部２１０と、分散共分散行列算出部２２０と、固有値固有ベクトル算出部２３０と、基底行列生成部２４０と、画像パッチ射影部２５０とから構成され、画像パッチ集合に属する全ての画像パッチに対し、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【００４８】
既述したように、一つの画像パッチは、一つの列ベクトルに対応する。よって、画像パッチ集合を列ベクトルの集合とすることが可能になり、画像パッチ集合を意味する列ベクトル並べた行列
は、下記数１によって表すことができる。以下、この行列
を単に「画像パッチ集合」とも言う。
【００４９】
【数１】
ここで、
は入力画像におけるｋ番目の画像パッチに対応する列ベクトルを表す。以下、この列ベクトル
を単に「ｋ番目の画像パッチ」、又は
とも言う。また、Ｎは画像パッチ集合の要素数を表す。つまり、Ｎは入力画像から得られた複数の画像パッチの個数である。
【００５０】
図３に示すように、ノイズ低減処理部２０では、まず、平均ベクトル算出部２１０が、画像パッチ集合に属する全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する。
【００５１】
ここで、説明の便宜上、ｋ番目の画像パッチに対し、平均ベクトル算出部２１０にて行われる平均ベクトル算出処理だけを説明する。
【００５２】
平均ベクトル算出部２１０では、下記数２及び数３に基づき、平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出する。
【００５３】
【数２】
【００５４】
【数３】
ここで、
は、画像パッチ
と画像パッチ
の間の類似度に基づく重み関数である。画像パッチ
と画像パッチ
の間の類似度が大きければ、重み
が大きくなるように設定される。即ち、ノイズ低減装置１で使用される重み
は、入力画像の画像パッチ間の類似度に基づいて設定されている。重み関数
の具体例については、後述する。
【００５５】
次に、ノイズ低減処理部２０では、分散共分散行列算出部２２０が、画像パッチ集合、及び平均ベクトル算出部２１０で算出された平均ベクトルに基づき、画像パッチ集合に属する全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する。
【００５６】
分散共分散行列算出部２２０では、下記数４に基づき、分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する。
【００５７】
【数４】
ここで、重み関数
の具体実施例について説明する。
【００５８】
第１例として、下記数５で表す関数を重み関数
として用いることができる。
【００５９】
【数５】
ただし、
は類似度評価における標準偏差を表す。また、
はユークリッドノルムを表す。
【００６０】
また、第２例として、下記数６で表す二値出力の重み関数を重み関数
として用いることができる。
【００６１】
【数６】
ただし、τは所定の定数である。
【００６２】
更に、第３例として、下記数７で表す重み関数を重み関数
として用いることができる。
【００６３】
【数７】
次に、ノイズ低減処理部２０では、固有値固有ベクトル算出部２３０が、固有値固有ベクトル算出処理を行うことにより、分散共分散行列算出部２２０で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する。
【００６４】
次に、ノイズ低減処理部２０では、基底行列生成部２４０が、固有値固有ベクトル算出部２３０で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する。
【００６５】
ここで、基底行列生成部２４０の２つの実施形態について説明する。
【００６６】
基底行列生成部２４０の第１実施形態では、下記数８に基づき、基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列
を生成する。
【００６７】
【数８】
ただし、
は画像パッチ
に対応する分散共分散行列
に大きな固有値に対応する固有ベクトルである。また、Ｌは予め設定された所定のパラメータであり、基底の個数を表す。
【００６８】
また、基底行列生成部２４０の第２実施形態では、固有値の値を利用して、上記数８に基づき、基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成する。
【００６９】
ただし、基底行列生成部２４０の第２実施形態では、
を満足するとする。ここで、
は固有ベクトル
に対応する固有値である。また、ρは予め設定された所定のパラメータである。
【００７０】
なお、基底行列生成部２４０の第２実施形態では、パラメータρを予め設定せず、
に基づき、パラメータρを設定するようにしても良い。ここで、
は入力画像のノイズの分散を表す。また、ｒは予め設定された所定のパラメータである。一例として、例えば、ｒの値を５と設定する。つまり、パラメータρを入力画像のノイズの分散の定数倍と設定するようにしても良い。
【００７１】
次に、ノイズ低減処理部２０では、画像パッチ射影部２５０が、画像パッチ集合に属する全ての画像パッチに対し、各画像パッチ、平均ベクトル算出部２１０で算出された各画像パッチに対応する平均ベクトル、及び、基底行列生成部２４０で生成された各画像パッチに対応する基底行列に基づき、画像パッチの射影を算出する画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチ（即ち、射影された画像パッチ）を生成する。
【００７２】
画像パッチ射影部２５０では、画像パッチ
平均ベクトル
及び基底行列
を用いて、下記数９に基づき、画像パッチ射影処理を行うことにより、画像パッチ
のデノイズされた画像パッチ
を生成する。
【００７３】
【数９】
上述のようにして、本発明のノイズ低減装置１のノイズ低減処理部２０が、画像パッチ集合に属する全ての画像パッチのデノイズされた画像パッチを生成することができる。
【００７４】
上記から分かるように、本発明では、画像パッチに対して行われるノイズ低減処理は、画像パッチをアフィン空間に射影する、画像パッチ射影処理を意味し、アフィン空間に射影された画像パッチは、デノイズされた画像パッチとなる。
【００７５】
つまり、画像パッチ
のデノイズされた画像パッチ
は、画像パッチ
をアフィン空間に射影することで得られた、アフィン空間に射影された画像パッチである。
【００７６】
ただし、平均ベクトル算出部２１０で算出された平均ベクトル
は、画像パッチ
に対応するアフィン空間の原点となる。また、基底行列生成部２４０で生成された基底行列
は、画像パッチ
に対応するアフィン空間の基底行列となる。
【００７７】
要するに、ノイズ低減処理部２０にて行われる「画像パッチ毎に、異なるパラメータによるノイズ低減処理」とは、各画像パッチを異なるアフィン空間（異なる原点と基底行列を持つアフィン空間）に射影することである。
【００７８】
図４は、本発明に係るノイズ低減装置の第２実施形態（以下、単に「ノイズ低減装置２」とも言う。）を示すブロック構成図である。以下、図４を用いて、ノイズ低減装置２を詳細に説明する。
【００７９】
図４に示すように、ノイズ低減装置２は、ノイズを含む入力画像から、デノイズされた画像を得るための装置であり、入力画像に対し、フィルタリング処理を行うフィルタリング処理部４０と、入力画像から画像パッチ集合を生成するとともに、フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成する画像パッチ集合生成部１１と、画像パッチ集合及びフィルタリングされた画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成するノイズ低減処理部２１と、ノイズ低減処理部２１で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を得る画像合成部３０とから構成される。
【００８０】
図４に示すように、ノイズ低減装置２では、まず、フィルタリング処理部４０が、入力画像に対し、所定のフィルタを用いてフィルタリング処理を行うことにより、フィルタリングされた画像を生成する。
【００８１】
ここで、所定のフィルタとして、例えば、ローパスフィルタや、バイラテラルフィルタを用いることができる。
【００８２】
また、離散コサイン変換（ＤＣＴ）を用いて入力画像を周波数領域へ変換し、変換された周波数領域で閾値処理を施す処理を、フィルタリング処理部４０にて行われるフィルタリング処理とすることもできる。
【００８３】
次に、ノイズ低減装置２では、図４に示すように、画像パッチ集合生成部１１が、入力画像から当該入力画像の画像パッチ集合を生成するとともに、フィルタリング処理部４０で得られたフィルタリングされた画像からフィルタリングされた画像パッチ集合を生成する。
【００８４】
具体的に、画像パッチ集合生成部１１は、入力画像に対して、画像パッチ集合生成処理を行うことにより、画像パッチ集合を生成するとともに、フィルタリングされた画像に対して、フィルタリングされた画像パッチ集合生成処理を行うことにより、フィルタリングされた画像パッチ集合を生成する。
【００８５】
ここで、画像パッチ集合生成処理は、前述したノイズ低減装置１の画像パッチ集合生成部１０にて行われる画像パッチ集合生成処理と同じであるため、その詳細な説明は省略する。
【００８６】
また、フィルタリングされた画像パッチ集合生成処理は、基本的に、画像パッチ集合生成処理と同じであるため、その詳細な説明は省略する。
【００８７】
ただし、画像パッチ集合生成処理は、入力画像に対して、その画素位置をずらしながら、複数の画像パッチを生成する処理であるのに対し、フィルタリングされた画像パッチ集合生成処理は、フィルタリングされた画像に対して、その画素位置をずらしながら、複数のフィルタリングされた画像パッチを生成する処理である。
【００８８】
次に、ノイズ低減装置２では、図４に示すように、ノイズ低減処理部２１が、画像パッチ集合生成部１１で生成された画像パッチ集合及びフィルタリングされた画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【００８９】
そして、ノイズ低減装置２では、図４に示すように、画像合成部３０が、ノイズ低減処理部２１で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成する。
【００９０】
ここで、画像合成処理は、前述したノイズ低減装置１の画像合成部３０にて行われる画像合成処理と同じであるため、その詳細な説明は省略する。
【００９１】
このようにして、本発明のノイズ低減装置２によれば、ノイズを含む入力画像から、その画像のデノイズされた画像を得ることができる。
【００９２】
次に、本発明に係るノイズ低減装置２のノイズ低減処理部２１の実施形態について説明する。ノイズ低減処理部２１は、画像パッチ集合及びフィルタリングされた画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【００９３】
図５は、本発明に係るノイズ低減装置２のノイズ低減処理部２１の実施形態を示すブロック構成図である。
【００９４】
図５に示されたように、ノイズ低減処理部２１は、平均ベクトル算出部２１１と、分散共分散行列算出部２２１と、固有値固有ベクトル算出部２３０と、基底行列生成部２４０と、画像パッチ射影部２５０とから構成され、画像パッチ集合及びフィルタリングされた画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対し、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【００９５】
図５に示すように、ノイズ低減処理部２１では、まず、平均ベクトル算出部２１１が、画像パッチ集合及びフィルタリングされた画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する。
【００９６】
ここで、説明の便宜上、ｋ番目の画像パッチに対し、平均ベクトル算出部２１１にて行われる平均ベクトル算出処理だけを説明する。
【００９７】
平均ベクトル算出部２１１では、下記数１０及び数１１に基づき、平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出する。
【００９８】
【数１０】
【００９９】
【数１１】
ここで、
は、画像パッチ集合生成部１１で生成されたフィルタリングされた画像パッチ集合におけるｋ番目のフィルタリングされた画像パッチを表し、以下、単に
とも言う。
はフィルタリングされた画像パッチ
とフィルタリングされた画像パッチ
の間の類似度に基づく重み関数である。フィルタリングされた画像パッチ
とフィルタリングされた画像パッチ
の間の類似度が大きければ、重み
が大きくなるように設定される。即ち、重み関数
は、重み関数
と同じように設定すれば良いので、その具体例の説明は省略する。ただし、ノイズ低減装置２で使用される重み
は、フィルタリングされた画像パッチ間の類似度に基づいて設定されている。
【０１００】
次に、ノイズ低減処理部２１では、分散共分散行列算出部２２１が、画像パッチ集合、フィルタリングされた画像パッチ集合、及び平均ベクトル算出部２１１で算出された平均ベクトルに基づき、画像パッチ集合に属する全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する。
【０１０１】
分散共分散行列算出部２２１では、下記数１２に基づき、分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する。
【０１０２】
【数１２】
次に、ノイズ低減処理部２１では、固有値固有ベクトル算出部２３０が、固有値固有ベクトル算出処理を行うことにより、分散共分散行列算出部２２１で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する。
【０１０３】
次に、ノイズ低減処理部２１では、基底行列生成部２４０が、固有値固有ベクトル算出部２３０で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する。
【０１０４】
ここで、ノイズ低減処理部２１の基底行列生成部２４０にて行われる処理は、前述したノイズ低減処理部２０の基底行列生成部２４０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１０５】
次に、ノイズ低減処理部２１では、画像パッチ射影部２５０が、画像パッチ集合に属する全ての画像パッチに対し、各画像パッチ、平均ベクトル算出部２１１で算出された各画像パッチに対応する平均ベクトル、及び、基底行列生成部２４０で生成された各画像パッチに対応する基底行列に基づき、画像パッチの射影を算出する画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチ（即ち、射影された画像パッチ）を生成する。
【０１０６】
ここで、ノイズ低減処理部２１の画像パッチ射影部２５０にて行われる処理は、前述したノイズ低減処理部２０の画像パッチ射影部２５０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１０７】
上述のようにして、本発明のノイズ低減装置２のノイズ低減処理部２１が、画像パッチ集合及びフィルタリングされた画像パッチ集合に基づいて、画像パッチ集合に属する全ての画像パッチのデノイズされた画像パッチを生成することができる。
【０１０８】
図６は、本発明に係るノイズ低減装置の第３実施形態（以下、単に「ノイズ低減装置３」とも言う。）を示すブロック構成図である。以下、図６を用いて、ノイズ低減装置３を詳細に説明する。
【０１０９】
図６に示すように、ノイズ低減装置３は、ノイズを含む入力画像から、デノイズされた画像を得るための装置であり、予め撮影された一枚以上の画像を格納する画像データベース５０と、入力画像から画像パッチ集合を生成するとともに、画像データベース５０に格納されている画像からデータベース画像パッチ集合を生成する画像パッチ集合生成部１２と、画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、ノイズ低減処理を行うことにより、デノイズされた画像パッチを生成するノイズ低減処理部２２と、ノイズ低減処理部２２で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を得る画像合成部３０とから構成される。
【０１１０】
図６に示すように、ノイズ低減装置３では、まず、画像パッチ集合生成部１２が、入力画像から当該入力画像の画像パッチ集合を生成するとともに、画像データベース５０に格納されている画像からデータベース画像パッチ集合を生成する。
【０１１１】
具体的に、画像パッチ集合生成部１２は、入力画像に対して、画像パッチ集合生成処理を行うことにより、画像パッチ集合を生成するとともに、画像データベース５０に格納されている画像に対して、データベース画像パッチ集合生成処理を行うことにより、データベース画像パッチ集合を生成する。
【０１１２】
ここで、画像パッチ集合生成処理は、前述したノイズ低減装置１の画像パッチ集合生成部１０にて行われる画像パッチ集合生成処理と同じであるため、その詳細な説明は省略する。
【０１１３】
また、データベース画像パッチ集合生成処理は、基本的に、画像パッチ集合生成処理と同じであるため、その詳細な説明は省略する。
【０１１４】
ただし、画像パッチ集合生成処理は、入力画像に対して、その画素位置をずらしながら、複数の画像パッチを生成する処理であるのに対し、データベース画像パッチ集合生成処理は、画像データベース５０に格納されている画像に対して、その画素位置をずらしながら、複数のデータベース画像パッチを生成する処理である。
【０１１５】
次に、ノイズ低減装置３では、図６に示すように、ノイズ低減処理部２２が、画像パッチ集合生成部１２で生成された画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【０１１６】
そして、ノイズ低減装置３では、図６に示すように、画像合成部３０が、ノイズ低減処理部２２で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成する。
【０１１７】
ここで、画像合成処理は、前述したノイズ低減装置１の画像合成部３０にて行われる画像合成処理と同じであるため、その詳細な説明は省略する。
【０１１８】
このようにして、本発明のノイズ低減装置３によれば、ノイズを含む入力画像から、その画像のデノイズされた画像を得ることができる。
【０１１９】
次に、本発明に係るノイズ低減装置３のノイズ低減処理部２２の実施形態について説明する。ノイズ低減処理部２２は、画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【０１２０】
図７は、本発明に係るノイズ低減装置３のノイズ低減処理部２２の実施形態を示すブロック構成図である。
【０１２１】
図７に示されたように、ノイズ低減処理部２２は、平均ベクトル算出部２１２と、分散共分散行列算出部２２２と、固有値固有ベクトル算出部２３０と、基底行列生成部２４０と、画像パッチ射影部２５０とから構成され、画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対し、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【０１２２】
図７に示すように、ノイズ低減処理部２２では、まず、平均ベクトル算出部２１２が、画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する。
【０１２３】
ここで、説明の便宜上、ｋ番目の画像パッチに対し、平均ベクトル算出部２１２にて行われる平均ベクトル算出処理だけを説明する。
【０１２４】
平均ベクトル算出部２１２では、下記数１３及び数１４に基づき、平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出する。
【０１２５】
【数１３】
【０１２６】
【数１４】
ここで、
は、画像パッチ集合生成部１２で生成されたデータベース画像パッチ集合におけるｉ番目のデータベース画像パッチを表し、以下、単に
とも言う。Ｍはデータベース画像パッチ集合の要素数である。
はデータベース画像パッチ
と画像パッチ
の間の類似度に基づく重み関数である。データベース画像パッチ
と画像パッチ
の間の類似度が大きければ、重み
が大きくなるように設定される。即ち、重み関数
は、重み関数
と同じように設定すれば良いので、その具体例の説明は省略する。ただし、ノイズ低減装置３で使用される重み
は、データベース画像パッチと画像パッチの間の類似度に基づいて設定されている。
【０１２７】
次に、ノイズ低減処理部２２では、分散共分散行列算出部２２２が、画像パッチ集合、データベース画像パッチ集合、及び平均ベクトル算出部２１２で算出された平均ベクトルに基づき、画像パッチ集合に属する全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する。
【０１２８】
分散共分散行列算出部２２２では、下記数１５に基づき、分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する。
【０１２９】
【数１５】
次に、ノイズ低減処理部２２では、固有値固有ベクトル算出部２３０が、固有値固有ベクトル算出処理を行うことにより、分散共分散行列算出部２２２で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する。
【０１３０】
次に、ノイズ低減処理部２２では、基底行列生成部２４０が、固有値固有ベクトル算出部２３０で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する。
【０１３１】
ここで、ノイズ低減処理部２２の基底行列生成部２４０にて行われる処理は、前述したノイズ低減処理部２０の基底行列生成部２４０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１３２】
次に、ノイズ低減処理部２２では、画像パッチ射影部２５０が、画像パッチ集合に属する全ての画像パッチに対し、各画像パッチ、平均ベクトル算出部２１２で算出された各画像パッチに対応する平均ベクトル、及び、基底行列生成部２４０で生成された各画像パッチに対応する基底行列に基づき、画像パッチの射影を算出する画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する。
【０１３３】
ここで、ノイズ低減処理部２２の画像パッチ射影部２５０にて行われる処理は、前述したノイズ低減処理部２０の画像パッチ射影部２５０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１３４】
上述のようにして、本発明のノイズ低減装置３のノイズ低減処理部２２が、画像パッチ集合及びデータベース画像パッチ集合に基づいて、画像パッチ集合に属する全ての画像パッチのデノイズされた画像パッチを生成することができる。
【０１３５】
図８は、本発明に係るノイズ低減装置の第４実施形態（以下、単に「ノイズ低減装置４」とも言う。）を示すブロック構成図である。以下、図８を用いて、ノイズ低減装置４を詳細に説明する。
【０１３６】
図８に示すように、ノイズ低減装置４は、ノイズを含む入力画像から、デノイズされた画像を得るための装置であり、入力画像に対し、フィルタリング処理を行うフィルタリング処理部４０と、予め撮影された一枚以上の画像を格納する画像データベース５０と、フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成するとともに、画像データベース５０に格納されている画像からデータベース画像パッチ集合を生成し、さらに、入力画像から画像パッチ集合を生成する画像パッチ集合生成部１３と、画像パッチ集合、フィルタリングされた画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成するノイズ低減処理部２３と、ノイズ低減処理部２３で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を得る画像合成部３０とから構成される。
【０１３７】
図８に示すように、ノイズ低減装置４では、まず、フィルタリング処理部４０が、入力画像に対し、所定のフィルタを用いてフィルタリング処理を行うことにより、フィルタリングされた画像を生成する。
【０１３８】
ここで、ノイズ低減装置４のフィルタリング処理部４０にて行われる処理は、前述したノイズ低減装置２のフィルタリング処理部４０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１３９】
次に、ノイズ低減装置４では、図８に示すように、画像パッチ集合生成部１３が、フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成するとともに、画像データベース５０に格納されている画像からデータベース画像パッチ集合を生成し、さらに、入力画像から画像パッチ集合を生成する。
【０１４０】
具体的に、画像パッチ集合生成部１３は、フィルタリングされた画像に対して、フィルタリングされた画像パッチ集合生成処理を行うことにより、フィルタリングされた画像パッチ集合を生成するとともに、画像データベース５０に格納されている画像に対して、データベース画像パッチ集合生成処理を行うことにより、データベース画像パッチ集合を生成し、さらに、入力画像に対して、画像パッチ集合生成処理を行うことにより、画像パッチ集合を生成する。
【０１４１】
ここで、フィルタリングされた画像パッチ集合生成処理は、前述したノイズ低減装置２の画像パッチ集合生成部１１にて行われるフィルタリングされた画像パッチ集合生成処理と同じであるため、その詳細な説明は省略する。
【０１４２】
また、データベース画像パッチ集合生成処理は、前述したノイズ低減装置３の画像パッチ集合生成部１２にて行われるデータベース画像パッチ集合生成処理と同じであるため、その詳細な説明は省略する。
【０１４３】
また、画像パッチ集合生成処理は、前述したノイズ低減装置１の画像パッチ集合生成部１０にて行われる画像パッチ集合生成処理と同じであるため、その詳細な説明は省略する。
【０１４４】
次に、ノイズ低減装置４では、図８に示すように、ノイズ低減処理部２３が、画像パッチ集合生成部１３で生成された、画像パッチ集合、フィルタリングされた画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【０１４５】
そして、ノイズ低減装置４では、図８に示すように、画像合成部３０が、ノイズ低減処理部２３で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成する。
【０１４６】
ここで、画像合成処理は、前述したノイズ低減装置１の画像合成部３０にて行われる画像合成処理と同じであるため、その詳細な説明は省略する。
【０１４７】
このようにして、本発明のノイズ低減装置４によれば、ノイズを含む入力画像から、その画像のデノイズされた画像を得ることができる。
【０１４８】
次に、本発明に係るノイズ低減装置４のノイズ低減処理部２３の実施形態について説明する。ノイズ低減処理部２３は、画像パッチ集合、フィルタリングされた画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【０１４９】
図９は、本発明に係るノイズ低減装置４のノイズ低減処理部２３の実施形態を示すブロック構成図である。
【０１５０】
図９に示されたように、ノイズ低減処理部２３は、平均ベクトル算出部２１３と、分散共分散行列算出部２２３と、固有値固有ベクトル算出部２３０と、基底行列生成部２４０と、画像パッチ射影部２５１とから構成され、画像パッチ集合、フィルタリングされた画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対し、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する。
【０１５１】
図９に示すように、ノイズ低減処理部２３では、まず、平均ベクトル算出部２１３が、画像パッチ集合、フィルタリングされた画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する。
【０１５２】
ここで、説明の便宜上、ｋ番目の画像パッチに対し、平均ベクトル算出部２１３にて行われる平均ベクトル算出処理だけを説明する。
【０１５３】
平均ベクトル算出部２１３では、下記数１６及び数１７に基づき、平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出する。
【０１５４】
【数１６】
【０１５５】
【数１７】
ここで、
はデータベース画像パッチ
とフィルタリングされた画像パッチ
の間の類似度に基づく重み関数である。データベース画像パッチ
とフィルタリングされた画像パッチ
の間の類似度が大きければ、重み
が大きくなるように設定される。即ち、重み関数
は、重み関数
と同じように設定すれば良いので、その具体例の説明は省略する。ただし、ノイズ低減装置４で使用される重み
は、データベース画像パッチとフィルタリングされた画像パッチの間の類似度に基づいて設定されている。
【０１５６】
次に、ノイズ低減処理部２３では、分散共分散行列算出部２２３が、画像パッチ集合、フィルタリングされた画像パッチ集合、データベース画像パッチ集合、及び平均ベクトル算出部２１３で算出された平均ベクトルに基づき、画像パッチ集合に属する全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する。
【０１５７】
分散共分散行列算出部２２３では、下記数１８に基づき、分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する。
【０１５８】
【数１８】
次に、ノイズ低減処理部２３では、固有値固有ベクトル算出部２３０が、固有値固有ベクトル算出処理を行うことにより、分散共分散行列算出部２２３で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する。
【０１５９】
次に、ノイズ低減処理部２３では、基底行列生成部２４０が、固有値固有ベクトル算出部２３０で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する。
【０１６０】
ここで、ノイズ低減処理部２３の基底行列生成部２４０にて行われる処理は、前述したノイズ低減処理部２０の基底行列生成部２４０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１６１】
次に、ノイズ低減処理部２３では、画像パッチ射影部２５０が、画像パッチ集合に属する全ての画像パッチに対し、各画像パッチ、平均ベクトル算出部２１３で算出された各画像パッチに対応する平均ベクトル、及び、基底行列生成部２４０で生成された各画像パッチに対応する基底行列に基づき、画像パッチの射影を算出する画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する。
【０１６２】
ここで、ノイズ低減処理部２３の画像パッチ射影部２５０にて行われる処理は、前述したノイズ低減処理部２０の画像パッチ射影部２５０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１６３】
上述のようにして、本発明のノイズ低減装置４のノイズ低減処理部２３が、画像パッチ集合、フィルタリングされた画像パッチ集合及びデータベース画像パッチ集合に基づいて、画像パッチ集合に属する全ての画像パッチのデノイズされた画像パッチを生成することができる。
【０１６４】
図１０は、本発明に係るノイズ低減装置の第５実施形態（以下、単に「ノイズ低減装置５」とも言う。）を示すブロック構成図である。また、図１１はノイズ低減装置５にて行われる処理の流れを示すフロー図である。以下、図１０及び図１１を用いて、ノイズ低減装置５を詳細に説明する。
【０１６５】
ノイズ低減装置５は、本発明に係るノイズ低減装置の第２実施形態と本発明に係るノイズ低減装置の第４実施形態を組み合わせた、繰り返し処理を有する実施形態である。
【０１６６】
図１０に示すように、ノイズ低減装置５は、ノイズを含む入力画像から、デノイズされた画像を得るための装置であり、フィルタリング処理部４０と、画像データベース５１と、画像パッチ集合生成部１４と、ノイズ低減処理部２４と、画像合成部３１とから構成される。
【０１６７】
ここで、ノイズを含む入力画像からデノイズされた画像を得るために、図１０と図１１を参照しながら、ノイズ低減装置５にて行われる処理の流れについて説明する。
【０１６８】
ノイズ低減装置５では、まず、繰り返し処理の回数ｋの値を０と設定する（図１１のステップＳ１００を参照）。
【０１６９】
次に、ノイズ低減装置５では、フィルタリング処理部４０が、入力画像に対し、所定のフィルタを用いてフィルタリング処理を行うことにより、フィルタリングされた画像を生成する（図１１のステップＳ１１０を参照）。
【０１７０】
ここで、ノイズ低減装置５のフィルタリング処理部４０にて行われる処理は、前述したノイズ低減装置２のフィルタリング処理部４０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１７１】
次に、ノイズ低減装置５では、画像パッチ集合生成部１４が、入力画像から当該入力画像の画像パッチ集合を生成するとともに、フィルタリング処理部４０で得られたフィルタリングされた画像からフィルタリングされた画像パッチ集合を生成する（図１１のステップＳ１２０を参照）。
【０１７２】
ここで、ノイズ低減装置５の画像パッチ集合生成部１４にて行われる処理は、前述したノイズ低減装置２の画像パッチ集合生成部１１にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１７３】
次に、ノイズ低減装置５では、ノイズ低減処理部２４が、画像パッチ集合生成部１４で生成された画像パッチ集合及びフィルタリングされた画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する（図１１のステップＳ１３０を参照）。
【０１７４】
ここで、ノイズ低減装置５のノイズ低減処理部２４にて行われる処理は、前述したノイズ低減装置２のノイズ低減処理部２１にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１７５】
次に、ノイズ低減装置５では、画像合成部３１が、ノイズ低減処理部２４で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成する（図１１のステップＳ１４０を参照）。
【０１７６】
ここで、ノイズ低減装置５の画像合成部３１にて行われる処理は、前述したノイズ低減装置２の画像合成部３０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１７７】
次に、ノイズ低減装置５では、デノイズされた画像を、画像データベース５１の画像として、画像データベース５１に格納する（図１１のステップＳ１５０を参照）。
【０１７８】
次に、ノイズ低減装置５では、画像パッチ集合生成部１４が、画像データベース５１に格納されている画像に対して、データベース画像パッチ集合生成処理を行うことにより、データベース画像パッチ集合を生成する（図１１のステップＳ１６０を参照）。
【０１７９】
ここで、ノイズ低減装置５の画像パッチ集合生成部１４にて行われるデータベース画像パッチ集合生成処理は、前述したノイズ低減装置４の画像パッチ集合生成部１３にて行われるデータベース画像パッチ集合生成処理と同じであるため、その詳細な説明は省略する。
【０１８０】
次に、ノイズ低減装置５では、ノイズ低減処理部２４が、画像パッチ集合生成部１４で生成された、画像パッチ集合、フィルタリングされた画像パッチ集合及びデータベース画像パッチ集合に基づき、画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する（図１１のステップＳ１７０を参照）。
【０１８１】
ここで、ノイズ低減装置５のノイズ低減処理部２４にて行われる処理は、前述したノイズ低減装置４のノイズ低減処理部２３にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１８２】
次に、ノイズ低減装置５では、画像合成部３１が、ノイズ低減処理部２４で得られた全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成する（図１１のステップＳ１８０を参照）。
【０１８３】
ここで、ノイズ低減装置５の画像合成部３１にて行われる処理は、前述したノイズ低減装置４の画像合成部３０にて行われる処理と同じであるため、その詳細な説明は省略する。
【０１８４】
次に、ノイズ低減装置５では、繰り返し処理の回数ｋの値を１増加させるようにする（図１１のステップＳ１９０を参照）。
【０１８５】
次に、ノイズ低減装置５では、繰り返し処理終了条件（即ち、ｋ＞所定の値）を満たすか否かを判定し（図１１のステップＳ２００を参照）、繰り返し処理終了条件（ｋ＞所定の値）を満たしたと判定した場合に、デノイズされた画像をノイズ低減装置５が生成したデノイズされた画像として外部に出力する（図１１のステップＳ２１０を参照）。これで、ノイズ低減装置５での処理が終了する。
【０１８６】
一方、ステップＳ２００において、ノイズ低減装置５が、繰り返し処理終了条件（ｋ＞所定の値）を満たさないと判定した場合に、ステップＳ１５０に戻り、ノイズ低減装置５では、ステップＳ１５０での処理、ステップＳ１６０での処理、ステップＳ１７０での処理、ステップＳ１８０での処理、ステップＳ１９０での処理、及びステップＳ２００での処理を繰り返すようにする。
【０１８７】
このようにして、本発明のノイズ低減装置５によれば、ノイズを含む入力画像から、その画像のデノイズされた画像を得ることができる。
【０１８８】
なお、本発明に係るノイズ低減装置は、コンピュータシステムを利用し、ソフトウェア（コンピュータプログラム）により実装されることができ、そして、ＡＳＩＣ（Application Specific Integrated Circuit）、ＧＰＵ（Graphics Processing Unit）やＦＰＧＡ（Field Programmable Gate Array）などのハードウェアにより実装されることも勿論できる。
【符号の説明】
【０１８９】
ノイズ低減装置 １,２,３,４,５
画像パッチ集合生成部 １０,１１,１２,１３,１４
ノイズ低減処理部 ２０,２１,２２,２３,２４
画像合成部 ３０,３１
フィルタリング処理部 ４０
画像データベース ５０,５１
平均ベクトル算出部 ２１０,２１１,２１２,２１３
分散共分散行列算出部 ２２０,２２１,２２２,２２３
固有値固有ベクトル算出部 ２３０
基底行列生成部 ２４０
画像パッチ射影部 ２５０
【特許請求の範囲】
【請求項１】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置であって、
前記画像から画像パッチ集合を生成する、画像パッチ集合生成部と、
前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部と、
前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成部と、
を備え、
前記ノイズ低減処理部は、
前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出部と、
前記画像パッチ集合、及び前記平均ベクトル算出部で算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出部と、
固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出部で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出部と、
前記固有値固有ベクトル算出部で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成部と、
前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出部で算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成部で生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影部とを備えることを特徴とするノイズ低減装置。
【請求項２】
前記平均ベクトル算出部では、次の数式に基づき、前記平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出し、
ここで、
は、画像パッチ
と画像パッチ
の間の類似度に基づく重み関数であり、Ｎは前記画像パッチ集合の要素数を表し、
前記分散共分散行列算出部では、
に基づき、前記分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する請求項１に記載のノイズ低減装置。
【請求項３】
前記基底行列生成部では、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列
を生成し、
ただし、
は画像パッチ
に対応する分散共分散行列
に大きな固有値に対応する固有ベクトルであり、Ｌは基底の個数を表す所定のパラメータである請求項２に記載のノイズ低減装置。
【請求項４】
前記基底行列生成部では、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、ρは予め設定された所定のパラメータである請求項２に記載のノイズ低減装置。
【請求項５】
前記基底行列生成部では、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、
に基づいてρを設定し、
は前記画像のノイズの分散を表し、ｒは予め設定された所定のパラメータである請求項２に記載のノイズ低減装置。
【請求項６】
前記画像パッチ射影部では、画像パッチ
平均ベクトル
及び基底行列
を用いて、
に基づき、前記画像パッチ射影処理を行うことにより、画像パッチ
のデノイズされた画像パッチ
を生成する請求項３乃至請求項５のいずれかに記載のノイズ低減装置。
【請求項７】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置であって、
前記画像に対し、フィルタリング処理を行うことにより、フィルタリングされた画像を生成する、フィルタリング処理部と、
前記画像から画像パッチ集合を生成するとともに、前記フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成する、画像パッチ集合生成部と、
前記画像パッチ集合及び前記フィルタリングされた画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部と、
前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成部と、
を備え、
前記ノイズ低減処理部は、
前記画像パッチ集合及び前記フィルタリングされた画像パッチ集合に基づき、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出部と、
前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記平均ベクトル算出部で算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出部と、
固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出部で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出部と、
前記固有値固有ベクトル算出部で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成部と、
前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出部で算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成部で生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影部とを備えることを特徴とするノイズ低減装置。
【請求項８】
前記平均ベクトル算出部では、次の数式に基づき、前記平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出し、
ここで、
はフィルタリングされた画像パッチ
とフィルタリングされた画像パッチ
の間の類似度に基づく重み関数であり、Ｎは前記画像パッチ集合の要素数を表し、
前記分散共分散行列算出部では、
に基づき、前記分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する請求項７に記載のノイズ低減装置。
【請求項９】
前記基底行列生成部では、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列
を生成し、
ただし、
は画像パッチ
に対応する分散共分散行列
に大きな固有値に対応する固有ベクトルであり、Ｌは基底の個数を表す所定のパラメータである請求項８に記載のノイズ低減装置。
【請求項１０】
前記基底行列生成部では、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、ρは予め設定された所定のパラメータである請求項８に記載のノイズ低減装置。
【請求項１１】
前記基底行列生成部では、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、
に基づいてρを設定し、
は前記画像のノイズの分散を表し、ｒは予め設定された所定のパラメータである請求項８に記載のノイズ低減装置。
【請求項１２】
前記画像パッチ射影部では、画像パッチ
平均ベクトル
及び基底行列
を用いて、
に基づき、前記画像パッチ射影処理を行うことにより、画像パッチ
のデノイズされた画像パッチ
を生成する請求項９乃至請求項１１のいずれかに記載のノイズ低減装置。
【請求項１３】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置であって、
予め撮影された一枚以上の画像を格納する画像データベースと、
前記画像から画像パッチ集合を生成するとともに、前記画像データベースに格納されている画像からデータベース画像パッチ集合を生成する、画像パッチ集合生成部と、
前記画像パッチ集合及び前記データベース画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部と、
前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成部と、
を備え、
前記ノイズ低減処理部は、
前記画像パッチ集合及び前記データベース画像パッチ集合に基づき、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出部と、
前記画像パッチ集合、前記データベース画像パッチ集合、及び前記平均ベクトル算出部で算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出部と、
固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出部で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出部と、
前記固有値固有ベクトル算出部で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成部と、
前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出部で算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成部で生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影部とを備えることを特徴とするノイズ低減装置。
【請求項１４】
前記平均ベクトル算出部では、次の数式に基づき、前記平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出し、
ここで、
はデータベース画像パッチ
と画像パッチ
の間の類似度に基づく重み関数であり、Ｍは前記データベース画像パッチ集合の要素数を表し、
前記分散共分散行列算出部では、
に基づき、前記分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する請求項１３に記載のノイズ低減装置。
【請求項１５】
前記基底行列生成部では、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列
を生成し、
ただし、
は画像パッチ
に対応する分散共分散行列
に大きな固有値に対応する固有ベクトルであり、Ｌは基底の個数を表す所定のパラメータである請求項１４に記載のノイズ低減装置。
【請求項１６】
前記基底行列生成部では、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、ρは予め設定された所定のパラメータである請求項１４に記載のノイズ低減装置。
【請求項１７】
前記基底行列生成部では、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、
に基づいてρを設定し、
は前記画像のノイズの分散を表し、ｒは予め設定された所定のパラメータである請求項１４に記載のノイズ低減装置。
【請求項１８】
前記画像パッチ射影部では、画像パッチ
平均ベクトル
及び基底行列
を用いて、
に基づき、前記画像パッチ射影処理を行うことにより、画像パッチ
のデノイズされた画像パッチ
を生成する請求項１５乃至請求項１７のいずれかに記載のノイズ低減装置。
【請求項１９】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置であって、
前記画像に対し、フィルタリング処理を行うことにより、フィルタリングされた画像を生成する、フィルタリング処理部と、
予め撮影された一枚以上の画像を格納する画像データベースと、
前記フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成するとともに、前記画像データベースに格納されている画像からデータベース画像パッチ集合を生成し、さらに、前記画像から画像パッチ集合を生成する、画像パッチ集合生成部と、
前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部と、
前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成部と、
を備え、
前記ノイズ低減処理部は、
前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づき、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出部と、
前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、前記データベース画像パッチ集合、及び前記平均ベクトル算出部で算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出部と、
固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出部で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出部と、
前記固有値固有ベクトル算出部で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成部と、
前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出部で算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成部で生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影部とを備えることを特徴とするノイズ低減装置。
【請求項２０】
前記平均ベクトル算出部では、次の数式に基づき、前記平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出し、
ここで、
はデータベース画像パッチ
とフィルタリングされた画像パッチ
の間の類似度に基づく重み関数であり、Ｍは前記データベース画像パッチ集合の要素数を表し、
前記分散共分散行列算出部では、
に基づき、前記分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する請求項１９に記載のノイズ低減装置。
【請求項２１】
前記基底行列生成部では、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列
を生成し、
ただし、
は画像パッチ
に対応する分散共分散行列
に大きな固有値に対応する固有ベクトルであり、Ｌは基底の個数を表す所定のパラメータである請求項２０に記載のノイズ低減装置。
【請求項２２】
前記基底行列生成部では、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、ρは予め設定された所定のパラメータである請求項２０に記載のノイズ低減装置。
【請求項２３】
前記基底行列生成部では、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、
に基づいてρを設定し、
は前記画像のノイズの分散を表し、ｒは予め設定された所定のパラメータである請求項２０に記載のノイズ低減装置。
【請求項２４】
前記画像パッチ射影部では、画像パッチ
平均ベクトル
及び基底行列
を用いて、
に基づき、前記画像パッチ射影処理を行うことにより、画像パッチ
のデノイズされた画像パッチ
を生成する請求項２１乃至請求項２３のいずれかに記載のノイズ低減装置。
【請求項２５】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置であって、
フィルタリング処理部と、画像データベースと、画像パッチ集合生成部と、ノイズ低減処理部と、画像合成部とを備え、
前記ノイズ低減装置では、
前記フィルタリング処理部が、前記画像に対し、フィルタリング処理を行うことにより、フィルタリングされた画像を生成し、
前記画像パッチ集合生成部が、前記画像から画像パッチ集合を生成するとともに、前記フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成し、
前記ノイズ低減処理部が、前記画像パッチ集合及び前記フィルタリングされた画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成し、
前記画像合成部が、前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成し、
前記デノイズされた画像が、前記画像データベースの画像として、前記画像データベースに格納され、
前記画像パッチ集合生成部が、前記画像データベースに格納されている画像に対して、データベース画像パッチ集合生成処理を行うことにより、データベース画像パッチ集合を生成し、
前記ノイズ低減処理部が、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づき、前記全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成し、
前記画像合成部が、前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成し、
繰り返し処理終了条件を満たしたと判定された場合に、前記デノイズされた画像を前記ノイズ低減装置が生成したデノイズされた画像とし、
前記繰り返し処理終了条件を満たさないと判定された場合に、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づいて行われるノイズ低減処理を繰り返し行うことにより、デノイズされた画像パッチを生成し、生成したデノイズされた画像パッチを前記画像データベースの画像とすることを特徴とするノイズ低減装置。
【請求項２６】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減方法であって、
前記画像から画像パッチ集合を生成する、画像パッチ集合生成ステップと、
前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理ステップと、
前記ノイズ低減処理ステップで生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成ステップと、
を有し、
前記ノイズ低減処理ステップは、
前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出ステップと、
前記画像パッチ集合、及び前記平均ベクトル算出ステップで算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出ステップと、
固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出ステップで算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出ステップと、
前記固有値固有ベクトル算出ステップで算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成ステップと、
前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出ステップで算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成ステップで生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影ステップとを有することを特徴とするノイズ低減方法。
【請求項２７】
前記平均ベクトル算出ステップでは、次の数式に基づき、前記平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出し、
ここで、
は、画像パッチ
と画像パッチ
の間の類似度に基づく重み関数であり、Ｎは前記画像パッチ集合の要素数を表し、
前記分散共分散行列算出ステップでは、
に基づき、前記分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する請求項２６に記載のノイズ低減方法。
【請求項２８】
前記基底行列生成ステップでは、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列
を生成し、
ただし、
は画像パッチ
に対応する分散共分散行列
に大きな固有値に対応する固有ベクトルであり、Ｌは基底の個数を表す所定のパラメータである請求項２７に記載のノイズ低減方法。
【請求項２９】
前記基底行列生成ステップでは、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、ρは予め設定された所定のパラメータである請求項２７に記載のノイズ低減方法。
【請求項３０】
前記基底行列生成ステップでは、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、
に基づいてρを設定し、
は前記画像のノイズの分散を表し、ｒは予め設定された所定のパラメータである請求項２７に記載のノイズ低減方法。
【請求項３１】
前記画像パッチ射影ステップでは、画像パッチ
平均ベクトル
及び基底行列
を用いて、
に基づき、前記画像パッチ射影処理を行うことにより、画像パッチ
のデノイズされた画像パッチ
を生成する請求項２８乃至請求項３０のいずれかに記載のノイズ低減方法。
【請求項３２】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減方法であって、
前記画像に対し、フィルタリング処理を行うことにより、フィルタリングされた画像を生成する、フィルタリング処理ステップと、
前記画像から画像パッチ集合を生成するとともに、前記フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成する、画像パッチ集合生成ステップと、
前記画像パッチ集合及び前記フィルタリングされた画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理ステップと、
前記ノイズ低減処理ステップで生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成ステップと、
を有し、
前記ノイズ低減処理ステップは、
前記画像パッチ集合及び前記フィルタリングされた画像パッチ集合に基づき、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出ステップと、
前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記平均ベクトル算出ステップで算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出ステップと、
固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出ステップで算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出ステップと、
前記固有値固有ベクトル算出ステップで算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成ステップと、
前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出ステップで算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成ステップで生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影ステップとを有することを特徴とするノイズ低減方法。
【請求項３３】
前記平均ベクトル算出ステップでは、次の数式に基づき、前記平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出し、
ここで、
はフィルタリングされた画像パッチ
とフィルタリングされた画像パッチ
の間の類似度に基づく重み関数であり、Ｎは前記画像パッチ集合の要素数を表し、
前記分散共分散行列算出部では、
に基づき、前記分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する請求項３２に記載のノイズ低減方法。
【請求項３４】
前記基底行列生成ステップでは、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列
を生成し、
ただし、
は画像パッチ
に対応する分散共分散行列
に大きな固有値に対応する固有ベクトルであり、Ｌは基底の個数を表す所定のパラメータである請求項３３に記載のノイズ低減方法。
【請求項３５】
前記基底行列生成ステップでは、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、ρは予め設定された所定のパラメータである請求項３３に記載のノイズ低減方法。
【請求項３６】
前記基底行列生成ステップでは、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、
に基づいてρを設定し、
は前記画像のノイズの分散を表し、ｒは予め設定された所定のパラメータである請求項３３に記載のノイズ低減方法。
【請求項３７】
前記画像パッチ射影ステップでは、画像パッチ
平均ベクトル
及び基底行列
を用いて、
に基づき、前記画像パッチ射影処理を行うことにより、画像パッチ
のデノイズされた画像パッチ
を生成する請求項３４乃至請求項３６のいずれかに記載のノイズ低減方法。
【請求項３８】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減方法であって、
予め撮影された一枚以上の画像を画像データベースに格納するステップと、
前記画像から画像パッチ集合を生成するとともに、前記画像データベースに格納されている画像からデータベース画像パッチ集合を生成する、画像パッチ集合生成ステップと、
前記画像パッチ集合及び前記データベース画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理ステップと、
前記ノイズ低減処理ステップで生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成ステップと、
を有し、
前記ノイズ低減処理ステップは、
前記画像パッチ集合及び前記データベース画像パッチ集合に基づき、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出ステップと、
前記画像パッチ集合、前記データベース画像パッチ集合、及び前記平均ベクトル算出ステップで算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出ステップと、
固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出ステップで算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出ステップと、
前記固有値固有ベクトル算出ステップで算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成ステップと、
前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出ステップで算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成ステップで生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影ステップとを有することを特徴とするノイズ低減方法。
【請求項３９】
前記平均ベクトル算出ステップでは、次の数式に基づき、前記平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出し、
ここで、
はデータベース画像パッチ
と画像パッチ
の間の類似度に基づく重み関数であり、Ｍは前記データベース画像パッチ集合の要素数を表し、
前記分散共分散行列算出ステップでは、
に基づき、前記分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する請求項３８に記載のノイズ低減方法。
【請求項４０】
前記基底行列生成ステップでは、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列
を生成し、
ただし、
は画像パッチ
に対応する分散共分散行列
に大きな固有値に対応する固有ベクトルであり、Ｌは基底の個数を表す所定のパラメータである請求項３９に記載のノイズ低減方法。
【請求項４１】
前記基底行列生成ステップでは、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、ρは予め設定された所定のパラメータである請求項３９に記載のノイズ低減方法。
【請求項４２】
前記基底行列生成ステップでは、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、
に基づいてρを設定し、
は前記画像のノイズの分散を表し、ｒは予め設定された所定のパラメータである請求項３９に記載のノイズ低減方法。
【請求項４３】
前記画像パッチ射影ステップでは、画像パッチ
平均ベクトル
及び基底行列
を用いて、
に基づき、前記画像パッチ射影処理を行うことにより、画像パッチ
のデノイズされた画像パッチ
を生成する請求項４０乃至請求項４２のいずれかに記載のノイズ低減方法。
【請求項４４】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減方法であって、
前記画像に対し、フィルタリング処理を行うことにより、フィルタリングされた画像を生成する、フィルタリング処理ステップと、
予め撮影された一枚以上の画像を画像データベースに格納するステップと、
前記フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成するとともに、前記画像データベースに格納されている画像からデータベース画像パッチ集合を生成し、さらに、前記画像から画像パッチ集合を生成する、画像パッチ集合生成ステップと、
前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理ステップと、
前記ノイズ低減処理ステップで生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成ステップと、
を有し、
前記ノイズ低減処理ステップは、
前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づき、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出ステップと、
前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、前記データベース画像パッチ集合、及び前記平均ベクトル算出ステップで算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出ステップと、
固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出ステップで算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出ステップと、
前記固有値固有ベクトル算出ステップで算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成ステップと、
前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出ステップで算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成ステップで生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影ステップとを有することを特徴とするノイズ低減方法。
【請求項４５】
前記平均ベクトル算出ステップでは、次の数式に基づき、前記平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出し、
ここで、
はデータベース画像パッチ
とフィルタリングされた画像パッチ
の間の類似度に基づく重み関数であり、Ｍは前記データベース画像パッチ集合の要素数を表し、
前記分散共分散行列算出ステップでは、
に基づき、前記分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する請求項４４に記載のノイズ低減方法。
【請求項４６】
前記基底行列生成ステップでは、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列
を生成し、
ただし、
は画像パッチ
に対応する分散共分散行列
に大きな固有値に対応する固有ベクトルであり、Ｌは基底の個数を表す所定のパラメータである請求項４５に記載のノイズ低減方法。
【請求項４７】
前記基底行列生成ステップでは、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、ρは予め設定された所定のパラメータである請求項４５に記載のノイズ低減方法。
【請求項４８】
前記基底行列生成ステップでは、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、
に基づいてρを設定し、
は前記画像のノイズの分散を表し、ｒは予め設定された所定のパラメータである請求項４５に記載のノイズ低減方法。
【請求項４９】
前記画像パッチ射影ステップでは、画像パッチ
平均ベクトル
及び基底行列
を用いて、
に基づき、前記画像パッチ射影処理を行うことにより、画像パッチ
のデノイズされた画像パッチ
を生成する請求項４６乃至請求項４８のいずれかに記載のノイズ低減方法。
【請求項５０】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減方法であって、
フィルタリング処理ステップと、画像パッチ集合生成ステップと、ノイズ低減処理ステップと、画像合成ステップとを有し、
前記ノイズ低減方法では、
前記フィルタリング処理部ステップが、前記画像に対し、フィルタリング処理を行うことにより、フィルタリングされた画像を生成し、
前記画像パッチ集合生成ステップが、前記画像から画像パッチ集合を生成するとともに、前記フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成し、
前記ノイズ低減処理ステップが、前記画像パッチ集合及び前記フィルタリングされた画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成し、
前記画像合成ステップが、前記ノイズ低減処理ステップで生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成し、
前記デノイズされた画像が、画像データベースの画像として、前記画像データベースに格納され、
前記画像パッチ集合生成ステップが、前記画像データベースに格納されている画像に対して、データベース画像パッチ集合生成処理を行うことにより、データベース画像パッチ集合を生成し、
前記ノイズ低減処理ステップが、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づき、前記全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成し、
前記画像合成ステップが、前記ノイズ低減処理ステップで生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成し、
繰り返し処理終了条件を満たしたと判定された場合に、前記デノイズされた画像を前記ノイズ低減装置が生成したデノイズされた画像とし、
前記繰り返し処理終了条件を満たさないと判定された場合に、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づいて行われるノイズ低減処理を繰り返し行うことにより、デノイズされた画像パッチを生成し、生成したデノイズされた画像パッチを前記画像データベースの画像とすることを特徴とするノイズ低減方法。
【請求項５１】
請求項２６乃至請求項５０の何れかに記載のノイズ低減方法をコンピュータに実行させるためのプログラム。
【請求項１】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置であって、
前記画像から画像パッチ集合を生成する、画像パッチ集合生成部と、
前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部と、
前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成部と、
を備え、
前記ノイズ低減処理部は、
前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出部と、
前記画像パッチ集合、及び前記平均ベクトル算出部で算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出部と、
固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出部で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出部と、
前記固有値固有ベクトル算出部で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成部と、
前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出部で算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成部で生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影部とを備えることを特徴とするノイズ低減装置。
【請求項２】
前記平均ベクトル算出部では、次の数式に基づき、前記平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出し、
ここで、
は、画像パッチ
と画像パッチ
の間の類似度に基づく重み関数であり、Ｎは前記画像パッチ集合の要素数を表し、
前記分散共分散行列算出部では、
に基づき、前記分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する請求項１に記載のノイズ低減装置。
【請求項３】
前記基底行列生成部では、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列
を生成し、
ただし、
は画像パッチ
に対応する分散共分散行列
に大きな固有値に対応する固有ベクトルであり、Ｌは基底の個数を表す所定のパラメータである請求項２に記載のノイズ低減装置。
【請求項４】
前記基底行列生成部では、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、ρは予め設定された所定のパラメータである請求項２に記載のノイズ低減装置。
【請求項５】
前記基底行列生成部では、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、
に基づいてρを設定し、
は前記画像のノイズの分散を表し、ｒは予め設定された所定のパラメータである請求項２に記載のノイズ低減装置。
【請求項６】
前記画像パッチ射影部では、画像パッチ
平均ベクトル
及び基底行列
を用いて、
に基づき、前記画像パッチ射影処理を行うことにより、画像パッチ
のデノイズされた画像パッチ
を生成する請求項３乃至請求項５のいずれかに記載のノイズ低減装置。
【請求項７】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置であって、
前記画像に対し、フィルタリング処理を行うことにより、フィルタリングされた画像を生成する、フィルタリング処理部と、
前記画像から画像パッチ集合を生成するとともに、前記フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成する、画像パッチ集合生成部と、
前記画像パッチ集合及び前記フィルタリングされた画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部と、
前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成部と、
を備え、
前記ノイズ低減処理部は、
前記画像パッチ集合及び前記フィルタリングされた画像パッチ集合に基づき、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出部と、
前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記平均ベクトル算出部で算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出部と、
固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出部で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出部と、
前記固有値固有ベクトル算出部で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成部と、
前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出部で算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成部で生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影部とを備えることを特徴とするノイズ低減装置。
【請求項８】
前記平均ベクトル算出部では、次の数式に基づき、前記平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出し、
ここで、
はフィルタリングされた画像パッチ
とフィルタリングされた画像パッチ
の間の類似度に基づく重み関数であり、Ｎは前記画像パッチ集合の要素数を表し、
前記分散共分散行列算出部では、
に基づき、前記分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する請求項７に記載のノイズ低減装置。
【請求項９】
前記基底行列生成部では、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列
を生成し、
ただし、
は画像パッチ
に対応する分散共分散行列
に大きな固有値に対応する固有ベクトルであり、Ｌは基底の個数を表す所定のパラメータである請求項８に記載のノイズ低減装置。
【請求項１０】
前記基底行列生成部では、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、ρは予め設定された所定のパラメータである請求項８に記載のノイズ低減装置。
【請求項１１】
前記基底行列生成部では、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、
に基づいてρを設定し、
は前記画像のノイズの分散を表し、ｒは予め設定された所定のパラメータである請求項８に記載のノイズ低減装置。
【請求項１２】
前記画像パッチ射影部では、画像パッチ
平均ベクトル
及び基底行列
を用いて、
に基づき、前記画像パッチ射影処理を行うことにより、画像パッチ
のデノイズされた画像パッチ
を生成する請求項９乃至請求項１１のいずれかに記載のノイズ低減装置。
【請求項１３】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置であって、
予め撮影された一枚以上の画像を格納する画像データベースと、
前記画像から画像パッチ集合を生成するとともに、前記画像データベースに格納されている画像からデータベース画像パッチ集合を生成する、画像パッチ集合生成部と、
前記画像パッチ集合及び前記データベース画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部と、
前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成部と、
を備え、
前記ノイズ低減処理部は、
前記画像パッチ集合及び前記データベース画像パッチ集合に基づき、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出部と、
前記画像パッチ集合、前記データベース画像パッチ集合、及び前記平均ベクトル算出部で算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出部と、
固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出部で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出部と、
前記固有値固有ベクトル算出部で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成部と、
前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出部で算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成部で生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影部とを備えることを特徴とするノイズ低減装置。
【請求項１４】
前記平均ベクトル算出部では、次の数式に基づき、前記平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出し、
ここで、
はデータベース画像パッチ
と画像パッチ
の間の類似度に基づく重み関数であり、Ｍは前記データベース画像パッチ集合の要素数を表し、
前記分散共分散行列算出部では、
に基づき、前記分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する請求項１３に記載のノイズ低減装置。
【請求項１５】
前記基底行列生成部では、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列
を生成し、
ただし、
は画像パッチ
に対応する分散共分散行列
に大きな固有値に対応する固有ベクトルであり、Ｌは基底の個数を表す所定のパラメータである請求項１４に記載のノイズ低減装置。
【請求項１６】
前記基底行列生成部では、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、ρは予め設定された所定のパラメータである請求項１４に記載のノイズ低減装置。
【請求項１７】
前記基底行列生成部では、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、
に基づいてρを設定し、
は前記画像のノイズの分散を表し、ｒは予め設定された所定のパラメータである請求項１４に記載のノイズ低減装置。
【請求項１８】
前記画像パッチ射影部では、画像パッチ
平均ベクトル
及び基底行列
を用いて、
に基づき、前記画像パッチ射影処理を行うことにより、画像パッチ
のデノイズされた画像パッチ
を生成する請求項１５乃至請求項１７のいずれかに記載のノイズ低減装置。
【請求項１９】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置であって、
前記画像に対し、フィルタリング処理を行うことにより、フィルタリングされた画像を生成する、フィルタリング処理部と、
予め撮影された一枚以上の画像を格納する画像データベースと、
前記フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成するとともに、前記画像データベースに格納されている画像からデータベース画像パッチ集合を生成し、さらに、前記画像から画像パッチ集合を生成する、画像パッチ集合生成部と、
前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理部と、
前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成部と、
を備え、
前記ノイズ低減処理部は、
前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づき、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出部と、
前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、前記データベース画像パッチ集合、及び前記平均ベクトル算出部で算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出部と、
固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出部で算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出部と、
前記固有値固有ベクトル算出部で算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成部と、
前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出部で算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成部で生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影部とを備えることを特徴とするノイズ低減装置。
【請求項２０】
前記平均ベクトル算出部では、次の数式に基づき、前記平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出し、
ここで、
はデータベース画像パッチ
とフィルタリングされた画像パッチ
の間の類似度に基づく重み関数であり、Ｍは前記データベース画像パッチ集合の要素数を表し、
前記分散共分散行列算出部では、
に基づき、前記分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する請求項１９に記載のノイズ低減装置。
【請求項２１】
前記基底行列生成部では、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列
を生成し、
ただし、
は画像パッチ
に対応する分散共分散行列
に大きな固有値に対応する固有ベクトルであり、Ｌは基底の個数を表す所定のパラメータである請求項２０に記載のノイズ低減装置。
【請求項２２】
前記基底行列生成部では、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、ρは予め設定された所定のパラメータである請求項２０に記載のノイズ低減装置。
【請求項２３】
前記基底行列生成部では、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、
に基づいてρを設定し、
は前記画像のノイズの分散を表し、ｒは予め設定された所定のパラメータである請求項２０に記載のノイズ低減装置。
【請求項２４】
前記画像パッチ射影部では、画像パッチ
平均ベクトル
及び基底行列
を用いて、
に基づき、前記画像パッチ射影処理を行うことにより、画像パッチ
のデノイズされた画像パッチ
を生成する請求項２１乃至請求項２３のいずれかに記載のノイズ低減装置。
【請求項２５】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減装置であって、
フィルタリング処理部と、画像データベースと、画像パッチ集合生成部と、ノイズ低減処理部と、画像合成部とを備え、
前記ノイズ低減装置では、
前記フィルタリング処理部が、前記画像に対し、フィルタリング処理を行うことにより、フィルタリングされた画像を生成し、
前記画像パッチ集合生成部が、前記画像から画像パッチ集合を生成するとともに、前記フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成し、
前記ノイズ低減処理部が、前記画像パッチ集合及び前記フィルタリングされた画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成し、
前記画像合成部が、前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成し、
前記デノイズされた画像が、前記画像データベースの画像として、前記画像データベースに格納され、
前記画像パッチ集合生成部が、前記画像データベースに格納されている画像に対して、データベース画像パッチ集合生成処理を行うことにより、データベース画像パッチ集合を生成し、
前記ノイズ低減処理部が、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づき、前記全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成し、
前記画像合成部が、前記ノイズ低減処理部で生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成し、
繰り返し処理終了条件を満たしたと判定された場合に、前記デノイズされた画像を前記ノイズ低減装置が生成したデノイズされた画像とし、
前記繰り返し処理終了条件を満たさないと判定された場合に、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づいて行われるノイズ低減処理を繰り返し行うことにより、デノイズされた画像パッチを生成し、生成したデノイズされた画像パッチを前記画像データベースの画像とすることを特徴とするノイズ低減装置。
【請求項２６】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減方法であって、
前記画像から画像パッチ集合を生成する、画像パッチ集合生成ステップと、
前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理ステップと、
前記ノイズ低減処理ステップで生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成ステップと、
を有し、
前記ノイズ低減処理ステップは、
前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出ステップと、
前記画像パッチ集合、及び前記平均ベクトル算出ステップで算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出ステップと、
固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出ステップで算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出ステップと、
前記固有値固有ベクトル算出ステップで算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成ステップと、
前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出ステップで算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成ステップで生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影ステップとを有することを特徴とするノイズ低減方法。
【請求項２７】
前記平均ベクトル算出ステップでは、次の数式に基づき、前記平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出し、
ここで、
は、画像パッチ
と画像パッチ
の間の類似度に基づく重み関数であり、Ｎは前記画像パッチ集合の要素数を表し、
前記分散共分散行列算出ステップでは、
に基づき、前記分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する請求項２６に記載のノイズ低減方法。
【請求項２８】
前記基底行列生成ステップでは、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列
を生成し、
ただし、
は画像パッチ
に対応する分散共分散行列
に大きな固有値に対応する固有ベクトルであり、Ｌは基底の個数を表す所定のパラメータである請求項２７に記載のノイズ低減方法。
【請求項２９】
前記基底行列生成ステップでは、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、ρは予め設定された所定のパラメータである請求項２７に記載のノイズ低減方法。
【請求項３０】
前記基底行列生成ステップでは、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、
に基づいてρを設定し、
は前記画像のノイズの分散を表し、ｒは予め設定された所定のパラメータである請求項２７に記載のノイズ低減方法。
【請求項３１】
前記画像パッチ射影ステップでは、画像パッチ
平均ベクトル
及び基底行列
を用いて、
に基づき、前記画像パッチ射影処理を行うことにより、画像パッチ
のデノイズされた画像パッチ
を生成する請求項２８乃至請求項３０のいずれかに記載のノイズ低減方法。
【請求項３２】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減方法であって、
前記画像に対し、フィルタリング処理を行うことにより、フィルタリングされた画像を生成する、フィルタリング処理ステップと、
前記画像から画像パッチ集合を生成するとともに、前記フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成する、画像パッチ集合生成ステップと、
前記画像パッチ集合及び前記フィルタリングされた画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理ステップと、
前記ノイズ低減処理ステップで生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成ステップと、
を有し、
前記ノイズ低減処理ステップは、
前記画像パッチ集合及び前記フィルタリングされた画像パッチ集合に基づき、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出ステップと、
前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記平均ベクトル算出ステップで算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出ステップと、
固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出ステップで算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出ステップと、
前記固有値固有ベクトル算出ステップで算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成ステップと、
前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出ステップで算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成ステップで生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影ステップとを有することを特徴とするノイズ低減方法。
【請求項３３】
前記平均ベクトル算出ステップでは、次の数式に基づき、前記平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出し、
ここで、
はフィルタリングされた画像パッチ
とフィルタリングされた画像パッチ
の間の類似度に基づく重み関数であり、Ｎは前記画像パッチ集合の要素数を表し、
前記分散共分散行列算出部では、
に基づき、前記分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する請求項３２に記載のノイズ低減方法。
【請求項３４】
前記基底行列生成ステップでは、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列
を生成し、
ただし、
は画像パッチ
に対応する分散共分散行列
に大きな固有値に対応する固有ベクトルであり、Ｌは基底の個数を表す所定のパラメータである請求項３３に記載のノイズ低減方法。
【請求項３５】
前記基底行列生成ステップでは、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、ρは予め設定された所定のパラメータである請求項３３に記載のノイズ低減方法。
【請求項３６】
前記基底行列生成ステップでは、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、
に基づいてρを設定し、
は前記画像のノイズの分散を表し、ｒは予め設定された所定のパラメータである請求項３３に記載のノイズ低減方法。
【請求項３７】
前記画像パッチ射影ステップでは、画像パッチ
平均ベクトル
及び基底行列
を用いて、
に基づき、前記画像パッチ射影処理を行うことにより、画像パッチ
のデノイズされた画像パッチ
を生成する請求項３４乃至請求項３６のいずれかに記載のノイズ低減方法。
【請求項３８】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減方法であって、
予め撮影された一枚以上の画像を画像データベースに格納するステップと、
前記画像から画像パッチ集合を生成するとともに、前記画像データベースに格納されている画像からデータベース画像パッチ集合を生成する、画像パッチ集合生成ステップと、
前記画像パッチ集合及び前記データベース画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理ステップと、
前記ノイズ低減処理ステップで生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成ステップと、
を有し、
前記ノイズ低減処理ステップは、
前記画像パッチ集合及び前記データベース画像パッチ集合に基づき、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出ステップと、
前記画像パッチ集合、前記データベース画像パッチ集合、及び前記平均ベクトル算出ステップで算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出ステップと、
固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出ステップで算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出ステップと、
前記固有値固有ベクトル算出ステップで算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成ステップと、
前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出ステップで算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成ステップで生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影ステップとを有することを特徴とするノイズ低減方法。
【請求項３９】
前記平均ベクトル算出ステップでは、次の数式に基づき、前記平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出し、
ここで、
はデータベース画像パッチ
と画像パッチ
の間の類似度に基づく重み関数であり、Ｍは前記データベース画像パッチ集合の要素数を表し、
前記分散共分散行列算出ステップでは、
に基づき、前記分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する請求項３８に記載のノイズ低減方法。
【請求項４０】
前記基底行列生成ステップでは、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列
を生成し、
ただし、
は画像パッチ
に対応する分散共分散行列
に大きな固有値に対応する固有ベクトルであり、Ｌは基底の個数を表す所定のパラメータである請求項３９に記載のノイズ低減方法。
【請求項４１】
前記基底行列生成ステップでは、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、ρは予め設定された所定のパラメータである請求項３９に記載のノイズ低減方法。
【請求項４２】
前記基底行列生成ステップでは、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、
に基づいてρを設定し、
は前記画像のノイズの分散を表し、ｒは予め設定された所定のパラメータである請求項３９に記載のノイズ低減方法。
【請求項４３】
前記画像パッチ射影ステップでは、画像パッチ
平均ベクトル
及び基底行列
を用いて、
に基づき、前記画像パッチ射影処理を行うことにより、画像パッチ
のデノイズされた画像パッチ
を生成する請求項４０乃至請求項４２のいずれかに記載のノイズ低減方法。
【請求項４４】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減方法であって、
前記画像に対し、フィルタリング処理を行うことにより、フィルタリングされた画像を生成する、フィルタリング処理ステップと、
予め撮影された一枚以上の画像を画像データベースに格納するステップと、
前記フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成するとともに、前記画像データベースに格納されている画像からデータベース画像パッチ集合を生成し、さらに、前記画像から画像パッチ集合を生成する、画像パッチ集合生成ステップと、
前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成する、ノイズ低減処理ステップと、
前記ノイズ低減処理ステップで生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、前記デノイズされた画像を合成する、画像合成ステップと、
を有し、
前記ノイズ低減処理ステップは、
前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づき、前記全ての画像パッチに対し、平均ベクトル算出処理を行うことにより、各画像パッチに対応する平均ベクトルを算出する、平均ベクトル算出ステップと、
前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、前記データベース画像パッチ集合、及び前記平均ベクトル算出ステップで算出された平均ベクトルに基づき、前記全ての画像パッチに対し、分散共分散行列算出処理を行うことにより、各画像パッチに対応する分散共分散行列を算出する、分散共分散行列算出ステップと、
固有値固有ベクトル算出処理を行うことにより、前記分散共分散行列算出ステップで算出された、各画像パッチに対応する分散共分散行列の固有値固有ベクトルをそれぞれ算出する、固有値固有ベクトル算出ステップと、
前記固有値固有ベクトル算出ステップで算出された、各画像パッチに対応する固有値固有ベクトルに基づき、基底行列生成処理を行うことにより、各画像パッチに対応する基底行列をそれぞれ生成する、基底行列生成ステップと、
前記全ての画像パッチに対し、各画像パッチ、前記平均ベクトル算出ステップで算出された各画像パッチに対応する平均ベクトル、及び、前記基底行列生成ステップで生成された各画像パッチに対応する基底行列に基づき、画像パッチ射影処理を行うことにより、各画像パッチのデノイズされた画像パッチを生成する、画像パッチ射影ステップとを有することを特徴とするノイズ低減方法。
【請求項４５】
前記平均ベクトル算出ステップでは、次の数式に基づき、前記平均ベクトル算出処理を行うことにより、画像パッチ
に対応する平均ベクトル
を算出し、
ここで、
はデータベース画像パッチ
とフィルタリングされた画像パッチ
の間の類似度に基づく重み関数であり、Ｍは前記データベース画像パッチ集合の要素数を表し、
前記分散共分散行列算出ステップでは、
に基づき、前記分散共分散行列算出処理を行うことにより、画像パッチ
に対応する分散共分散行列
を算出する請求項４４に記載のノイズ低減方法。
【請求項４６】
前記基底行列生成ステップでは、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列
を生成し、
ただし、
は画像パッチ
に対応する分散共分散行列
に大きな固有値に対応する固有ベクトルであり、Ｌは基底の個数を表す所定のパラメータである請求項４５に記載のノイズ低減方法。
【請求項４７】
前記基底行列生成ステップでは、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、ρは予め設定された所定のパラメータである請求項４５に記載のノイズ低減方法。
【請求項４８】
前記基底行列生成ステップでは、固有値の値を利用して、
に基づき、前記基底行列生成処理を行うことにより、画像パッチ
に対応する基底行列を生成し、
ただし、
が成立し、
は固有ベクトル
に対応する固有値であり、
に基づいてρを設定し、
は前記画像のノイズの分散を表し、ｒは予め設定された所定のパラメータである請求項４５に記載のノイズ低減方法。
【請求項４９】
前記画像パッチ射影ステップでは、画像パッチ
平均ベクトル
及び基底行列
を用いて、
に基づき、前記画像パッチ射影処理を行うことにより、画像パッチ
のデノイズされた画像パッチ
を生成する請求項４６乃至請求項４８のいずれかに記載のノイズ低減方法。
【請求項５０】
ノイズを含む一枚の画像から、当該画像のデノイズされた画像を生成するノイズ低減方法であって、
フィルタリング処理ステップと、画像パッチ集合生成ステップと、ノイズ低減処理ステップと、画像合成ステップとを有し、
前記ノイズ低減方法では、
前記フィルタリング処理部ステップが、前記画像に対し、フィルタリング処理を行うことにより、フィルタリングされた画像を生成し、
前記画像パッチ集合生成ステップが、前記画像から画像パッチ集合を生成するとともに、前記フィルタリングされた画像からフィルタリングされた画像パッチ集合を生成し、
前記ノイズ低減処理ステップが、前記画像パッチ集合及び前記フィルタリングされた画像パッチ集合に基づき、前記画像パッチ集合に属する全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成し、
前記画像合成ステップが、前記ノイズ低減処理ステップで生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成し、
前記デノイズされた画像が、画像データベースの画像として、前記画像データベースに格納され、
前記画像パッチ集合生成ステップが、前記画像データベースに格納されている画像に対して、データベース画像パッチ集合生成処理を行うことにより、データベース画像パッチ集合を生成し、
前記ノイズ低減処理ステップが、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づき、前記全ての画像パッチに対して、画像パッチ毎に、異なるパラメータによるノイズ低減処理を行うことにより、デノイズされた画像パッチを生成し、
前記画像合成ステップが、前記ノイズ低減処理ステップで生成された全てのデノイズされた画像パッチに対し、画像合成処理を行うことにより、デノイズされた画像を合成し、
繰り返し処理終了条件を満たしたと判定された場合に、前記デノイズされた画像を前記ノイズ低減装置が生成したデノイズされた画像とし、
前記繰り返し処理終了条件を満たさないと判定された場合に、前記画像パッチ集合、前記フィルタリングされた画像パッチ集合、及び前記データベース画像パッチ集合に基づいて行われるノイズ低減処理を繰り返し行うことにより、デノイズされた画像パッチを生成し、生成したデノイズされた画像パッチを前記画像データベースの画像とすることを特徴とするノイズ低減方法。
【請求項５１】
請求項２６乃至請求項５０の何れかに記載のノイズ低減方法をコンピュータに実行させるためのプログラム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【公開番号】特開２０１３−２６６６９（Ｐ２０１３−２６６６９Ａ）
【公開日】平成２５年２月４日（２０１３．２．４）
【国際特許分類】
【出願番号】特願２０１１−１５６６５９（Ｐ２０１１−１５６６５９）
【出願日】平成２３年７月１５日（２０１１．７．１５）
【出願人】（３０４０２１４１７）国立大学法人東京工業大学 (1,821)
【Ｆターム（参考）】
【公開日】平成２５年２月４日（２０１３．２．４）
【国際特許分類】
【出願日】平成２３年７月１５日（２０１１．７．１５）
【出願人】（３０４０２１４１７）国立大学法人東京工業大学 (1,821)
【Ｆターム（参考）】
[ Back to top ]