(57)【要約】【課題】データ信号の信号劣化を高精度に補償するための信号処理装置、信号処理方法およびフィルタ設計方法を提供する。【解決手段】インパルス応答が畳み込まれたデータ信号に補償フィルタを畳み込む畳み込み部を備え、前記補償フィルタの周波数応答は、高域になるほど０に漸近する、信号処理装置。【選択図】図１
2012141726.tif 000002
【特許請求の範囲】
【請求項１】
インパルス応答が畳み込まれたデータ信号に補償フィルタを畳み込む畳み込み部を備え、
前記補償フィルタの周波数応答は、前記インパルス応答の推定結果の周波数応答Ｈ、ＳＮ比の周波数応答Ｎ／Ｓを用いる数式１に基づく演算により表現される、信号処理装置。
【数１】
2012141726.tif 000007
【請求項２】
前記データ信号は撮像光学系により撮像された画像信号であり、
前記インパルス応答は、前記撮像光学系の点像分布関数である、請求項１に記載の信号処理装置。
【請求項３】
前記インパルス応答の推定結果は、ガウシアンを仮定する推定処理により得られる、請求項２に記載の信号処理装置。
【請求項４】
インパルス応答が畳み込まれたデータ信号に補償フィルタを畳み込むステップを有し、
前記補償フィルタの周波数応答は、前記インパルス応答の推定結果の周波数応答Ｈ、ＳＮ比の周波数応答Ｎ／Ｓを用いる数式２に基づく演算により表現される、信号処理方法。
【数２】
2012141726.tif 000008
【請求項５】
データ信号に畳み込まれるインパルス応答を推定するステップと；
前記インパルス応答をフーリエ変換するステップと；
前記インパルス応答のフーリエ変換結果Ｈ、ＳＮ比の周波数応答Ｎ／Ｓを用いる数式３に基づく演算により周波数応答が表現される補償フィルタを設計するステップと；
を含む、フィルタ設計方法。
【数３】
2012141726.tif 000009
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、信号処理装置、信号処理方法およびフィルタ設計方法に関する。
【背景技術】
【０００２】
近日、ＣＣＤやＣＭＯＳなどの撮像素子を用い、撮像レンズにより撮像素子の受光面上に結像された被写体の光学像を撮像する撮像装置が普及している。このような撮像装置により得られる画像は、受光面上の受光画素数や撮像レンズの解像力に応じた解像度を有する。
【０００３】
このため、撮像装置により得られる画像の解像度を向上するためには、受光画素の数を増大させることや、撮像レンズの解像力を高めることが効果的である。例えば、受光面上の受光画素の密度を高くすると共に、撮像レンズを通して受光面上に投影される点像が１つの受光画素の範囲内に収まるように撮像レンズの解像力を高めることにより、撮像装置により得られる画像の解像度を向上することが可能である。
【０００４】
ここで、撮像素子を構成する受光画素の密度を高めることは近年の技術向上により比較的容易である。一方、撮像レンズの解像力を高めるためには、撮像レンズを構成する各レンズの形状誤差や組立誤差などを低減することが有効であるが、加工、組立および調整制度などの製作精度をさらに高めることの困難性は高い。
【０００５】
また、解像度を高める他のアプローチとして、撮像装置の撮像により得られた画像を点像分布関数（ＰＳＦ：Ｐｏｉｎｔ  Ｓｐｒｅａｄ  Ｆｕｎｃｔｉｏｎ）に基づいて信号処理することにより解像度を高める方法が知られている。撮像装置により得られた画像は、撮像レンズのＭＴＦ（Ｍｏｄｕｌａｔｉｏｎ  Ｔｒａｎｓｆｅｒ  Ｆｕｎｃｔｉｏｎ）特性により劣化しているが、点像分布関数に基づいて設計される補償フィルタによりこの信号劣化を補償することができる。
【０００６】
また、例えば特許文献１に記載されているように、補償フィルタとして、周波数応答Ｒｗが以下の数式１のように表現されるウィナーフィルタが広く知られている。このウィナーフィルタによれば、信号対雑音比の周波数特性Ｎ／Ｓが既知である場合に、劣化前の原信号に対して二乗誤差が最小になる出力が得られる。なお、数式１におけるＨは点像分布関数の周波数応答を示す。
【０００７】
【数１】
2012141726.tif 000003
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特開２００５−２５２７６０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００９】
上記のウィナーフィルタに関し、信号劣化した画像に対してウィナーフィルタを適用すると、空間周波数の高い信号成分を補償でなく減衰させてしまうという性質があった。これは、点像分布関数の周波数応答Ｈは高い周波数ほど小さくなるので、ウィナーフィルタの周波数応答Ｒｗも周波数が高いほど０に近づくためである。
【００１０】
原信号がどのようなインパルス応答（点像分布関数）で劣化したが既知である場合にはウィナーフィルタの上述した性質は問題にならない。なぜなら、ウィナーフィルタが減衰させてしまう高域成分は、ノイズに埋もれているので、減衰させた方が２乗誤差が小さくなるからである。
【００１１】
しかし、現実にはどのようなインパルス応答によって原信号が劣化したのか正確に把握できない場合がある。このため、ガウシアンなどの典型的なインパルス応答を仮定して設計したウィナーフィルタを用いて補償を行うことも考えられるが、現実の原信号の高域成分がガウシアンほど減衰していない場合、ガウシアンを仮定して設計したウィナーフィルタを用いると高域成分を必要以上に減衰させてしまうという問題があった。
【００１２】
そこで、本発明は、上記問題に鑑みてなされたものであり、本発明の目的とするところは、インパルス応答の畳み込みによるデータ信号の劣化を高精度に補償することが可能な、新規かつ改良された信号処理装置、信号処理方法およびフィルタ設計方法を提供することにある。
【課題を解決するための手段】
【００１３】
上記課題を解決するために、本発明のある観点によれば、インパルス応答が畳み込まれたデータ信号に補償フィルタを畳み込む畳み込み部を備え、前記補償フィルタの周波数応答は、前記インパルス応答の推定結果の周波数応答Ｈ、ＳＮ比の周波数応答Ｎ／Ｓを用いる数式２に基づく演算により表現される、信号処理装置が提供される。
【００１４】
【数２】
2012141726.tif 000004
【００１５】
前記データ信号は撮像光学系により撮像された画像信号であり、前記インパルス応答は、前記撮像光学系の点像分布関数であってもよい。
【００１６】
前記インパルス応答の推定結果は、ガウシアンを仮定する推定処理により得られてもよい。
【００１７】
また、上記課題を解決するために、本発明の別の観点によれば、インパルス応答が畳み込まれたデータ信号に補償フィルタを畳み込むステップを有し、前記補償フィルタの周波数応答は、前記インパルス応答の推定結果の周波数応答Ｈ、ＳＮ比の周波数応答Ｎ／Ｓを用いる数式２に基づく演算により表現される、信号処理方法が提供される。
【００１８】
また、上記課題を解決するために、本発明の別の観点によれば、データ信号に畳み込まれるインパルス応答を推定するステップと、前記インパルス応答をフーリエ変換するステップと、前記インパルス応答のフーリエ変換結果Ｈ、ＳＮ比の周波数応答Ｎ／Ｓを用いる数式２に基づく演算により周波数応答が表現される補償フィルタを設計するステップと、を含む、フィルタ設計方法が提供される。
【発明の効果】
【００１９】
以上説明したように本発明によれば、インパルス応答の畳み込みによるデータ信号の劣化を高精度に補償することができる。また、現実より大きいボケを仮定した場合にも必要以上に高域が減衰されることを防止できる。
【図面の簡単な説明】
【００２０】
【図１】本発明の実施形態による画像処理装置の構成を示した説明図である。
【図２】本発明の実施形態によるフィルタ設計装置の構成を示した説明図である。
【図３】ＰＳＦ、ウィナーフィルタ、および本発明の実施形態による補償フィルタのインパルス応答の関係を示した説明図である。
【図４】本発明の実施形態による補償フィルタの周波数特性を示した説明図である。
【図５】本発明の実施形態による補償フィルタによる出力の周波数特性を示した説明図である。
【図６】補償前の入力画像の具体例を示した説明図である。
【図７】入力画像にウィナーフィルタを適用した場合の出力画像を示した説明図である。
【図８】入力画像に本発明の実施形態による補償フィルタを適用した場合の出力画像を示した説明図である。
【図９】本発明の実施形態の動作を示したシーケンス図である。
【発明を実施するための形態】
【００２１】
以下に添付図面を参照しながら、本発明の実施の形態について詳細に説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。
【００２２】
また、本明細書及び図面において、実質的に同一の機能構成を有する複数の構成要素を、同一の符号の後に異なるアルファベットを付して区別する場合もある。ただし、実質的に同一の機能構成を有する複数の構成要素の各々を特に区別する必要がない場合、同一符号のみを付する。
【００２３】
また、以下に示す項目順序に従って当該「発明を実施するための形態」を説明する。
１．画像処理装置の構成
２．フィルタ設計装置の構成
３．本発明の実施形態による効果
４．本発明の実施形態の動作
５．むすび
【００２４】
＜１．画像処理装置の構成＞
【００２５】
図１は、本発明の実施形態による画像処理装置１の構成を示した説明図である。図１に示したように、本発明の実施形態による画像処理装置１は、撮像光学系１０と、アイリス２と、ＣＭＯＳ撮像素子１４と、Ｓ／Ｈ＿ＡＧＣ１６と、Ａ／Ｄ変換部１８と、ＡＥ検波ブロック２４と、マイクロコンピュータ２６と、アイリスドライバ２８と、フィルタ係数保持部３０と、を備える。なお、本発明の実施形態による画像処理装置１は、撮像装置、携帯機器、車載機器および医療機器などの多様な機器に適用可能である。
【００２６】
撮像光学系１０は、複数の撮像レンズからなり、アイリス１２を介して被写体の光学像をＣＭＯＳ撮像素子１４の受光面に結像させる。
【００２７】
ＣＭＯＳ（Ｃｏｍｐｌｅｍｅｎｔａｒｙ  Ｍｅｔａｌ  Ｏｘｉｄｅ  Ｓｅｍｉｃｏｎｄｕｃｔｏｒ）撮像素子１４は、受光面に結像された被写体の光学像を電気信号に変換する。なお、ＣＭＯＳ撮像素子１４は撮像素子の一例に過ぎず、画像処理装置１は、ＣＭＯＳ撮像素子１４に代えて例えばＣＣＤ（Ｃｈａｒｇｅ  Ｃｏｕｐｌｅｄ  Ｄｅｖｉｃｅ）撮像素子を備えてもよい。
【００２８】
Ｓ／Ｈ＿ＡＧＣ１６は、ＣＭＯＳ撮像素子１４からの出力に対してサンプルホールド（Ｓ／Ｈ）およびＡＧＣ（Ａｕｔｏｍａｔｉｃ  Ｇａｉｎ  Ｃｏｎｔｒｏｌ）を施し、Ａ／Ｄ変換部１８は、Ｓ／Ｈ＿ＡＧＣ１６からの出力をデジタル形式に変換する。Ａ／Ｄ変換部１８からのデジタル信号は、信号処理部２０およびＡＥ検波ブロック２４に供給される。
【００２９】
ＡＥ（Ａｕｔｏ  Ｅｘｐｏｓｕｒｅ）検波ブロック２４は、Ａ／Ｄ変換部１８から入力されるデジタル信号を検波する。マイクロコンピュータ２６は、ＡＥ検波ブロック２４から供給される検波値に基づき、アイリス１２の開度を制御するための制御信号をアイリスドライバ２８に供給する。アイリスドライバ２８は、マイクロコンピュータ２６から供給される制御信号に従ってアイリス１２を駆動する。
【００３０】
信号処理部２０は、Ａ／Ｄ変換部１８からのデジタル出力に対して各種の信号処理を施し、後段回路に処理後の信号を出力する。例えば、信号処理部２０にＡ／Ｄ変換部１８からデジタル出力として入力される入力画像は、撮像光学系１０のＭＴＦ特性によりボケている。このため、信号処理部２０は、フィルタ係数保持部３０に保持されているフィルタ係数からなる補償フィルタを入力画像に畳み込むことにより、撮像光学系１０のＭＴＦ特性による像の劣化を補償する畳み込み部としての機能を有する。
【００３１】
フィルタ係数保持部３０は、信号処理部２０における畳み込みに用いられる補償フィルタのフィルタ係数を保持する。このフィルタ係数は、後述のフィルタ設計装置４０により設計され、フィルタ係数保持部３０に格納される。
【００３２】
（本発明の実施形態に至る経緯）
ところで、補償フィルタとして、周波数応答Ｒｗが以下の数式３のように表現されるウィナーフィルタが広く知られている。このウィナーフィルタによれば、信号対雑音比の周波数特性Ｎ（ｕ，ｖ）／Ｓ（ｕ，ｖ）が既知である場合に、劣化前の原信号に対して二乗誤差が最小になる出力が得られる。なお、数式１におけるＨ（ｕ，ｖ）は点像分布関数の周波数応答（ＯＴＦ：Ｏｐｔｉｃａｌ  Ｔｒａｎｓｆｅｒ  Ｆｕｎｃｔｉｏｎ）を示す。
【００３３】
【数３】
2012141726.tif 000005
【００３４】
上記のウィナーフィルタに関し、信号劣化した画像に対してウィナーフィルタを適用すると、空間周波数の高い信号成分を補償でなく減衰させてしまうという性質があった。これは、点像分布関数の周波数応答Ｈは高い周波数ほど小さくなるので、ウィナーフィルタの周波数応答Ｒｗも周波数が高いほど０に近づくためである。
【００３５】
原信号がどのようなインパルス応答（点像分布関数）で劣化したが既知である場合にはウィナーフィルタの上述した性質は問題にならない。なぜなら、ウィナーフィルタが減衰させてしまう高域成分は、ノイズに埋もれているので、減衰させた方が２乗誤差が小さくなるからである。
【００３６】
しかし、現実にはどのようなインパルス応答によって原信号が劣化したのか正確に把握できない場合がある。このため、ガウシアンなどの典型的なインパルス応答を仮定して設計したウィナーフィルタを用いて補償を行うことも考えられるが、現実の原信号の高域成分がガウシアンほど減衰していない場合、ガウシアンを仮定して設計したウィナーフィルタを用いると高域成分を必要以上に減衰させてしまうという問題があった。
【００３７】
そこで、上記事情を一着眼点にして本発明の実施形態を創作するに至った。本発明の実施形態によれば、高域になるにつれてゲインが０ｄＢに漸近する補償フィルタを用いることにより、インパルス応答の畳み込みによるデータ信号の劣化を高精度に補償することが可能となる。以下、このような補償フィルタを設計するフィルタ設計装置について詳細に説明する。
【００３８】
＜２．フィルタ設計装置の構成＞
図２は、本発明の実施形態によるフィルタ設計装置４０の構成を示した説明図である。図２に示したように、フィルタ設計装置４０は、ＰＳＦ推定部４２と、フーリエ変換部４４と、周波数応答設計部４６と、逆フーリエ変換部４８と、を備える。
【００３９】
ＰＳＦ推定部４２は、撮像光学系１０のＭＴＦ特性によるボケの大きさを、ボケがガウシアンであると仮定して推定し、推定結果であるＰＳＦ（２次元のインパルス応答）を出力する。このようにガウシアンを仮定してＰＳＦを推定する場合、実際のボケよりも大きなボケを推定してしまうことも想定されるが、本発明の実施形態による補償フィルタによれば、ボケを大きく推定してしまった場合にも画像の高域成分を保つことが可能である。なお、上記ではガウシアンを仮定してＰＳＦを推定する方法を説明したが、ＰＳＦの推定方法はかかる例に限定されず、ＰＳＦ推定部４２は、他の方法によりＰＳＦを推定してもよい。
【００４０】
フーリエ変換部４４は、ＰＳＦ推定部４２により推定されたＰＳＦを、周波数応答であるＯＴＦに変換する。
【００４１】
周波数応答設計部４６は、フーリエ変換部４４により得られたＯＴＦに基づき、ＯＴＦを示すＨ、ＳＮ比の周波数応答Ｎ／Ｓを用いる数式４に基づく演算により、補償フィルタの周波数応答Ｒｐを設計する。なお、周波数応答設計部４６は、撮像光学系１０により試験的に得られた画像を見ながら調整して得られたＳＮ比の周波数応答Ｎ／Ｓを用いてもよいし、測定により得られたＳＮ比の周波数応答Ｎ／Ｓを用いてもよい。また、ＳＮ比の大きさは必ずしも周波数に依存しなくてもよく、例えば、ＳＮ比は定数であってもよい。
【００４２】
また、本発明の補償フィルタの周波数応答は、数式４で示される周波数応答Ｒｐに限定されず、例えば、数式４のいずれかの項に定数を乗じる演算、数式４の分子または分母に新たな項を追加する演算など、数式４に基づく演算により表現される周波数応答であってもよい。
【００４３】
【数４】
2012141726.tif 000006
【００４４】
逆フーリエ変換部４８は、周波数応答設計部４６により設計された補償フィルタの周波数応答Ｒｐを、逆フーリエ変換によりインパルス応答に変換する。こうして逆フーリエ変換部４８によりインパルス応答に対応するフィルタ係数が得られ、この補償フィルタのフィルタ係数が画像処理装置１のフィルタ係数保持部３０に格納され、信号処理部２０における畳み込み処理に用いられる。
【００４５】
ここで、数式４に示したように、本発明の実施形態による補償フィルタは、数式３に示したウィナーフィルタの周波数応答Ｒｗの分子にＳＮ比の周波数応答Ｎ／Ｓが加えられた周波数応答Ｒｐを有する。このため、ＯＴＦが０に近づくと、ウィナーフィルタの周波数応答Ｒｗも０に近づくのに対し、本発明の実施形態による補償フィルタの周波数応答Ｒｐは１に近づく。すなわち、ＯＴＦによる減衰が多い帯域は、ウィナーフィルタを用いるとさらに減衰されるのに対し、本発明の実施形態による補償フィルタを用いると素通しになる。したがって、本発明の実施形態によれば、推定されたＯＴＦによる信号劣化ほど実際の高域成分が失われていなかった場合、高域成分を減衰させることなく他の周波数成分を補償することが可能である。以下、このような本発明の実施形態による効果を具体的に説明する。
【００４６】
＜３．本発明の実施形態による効果＞
図３は、ＰＳＦ、ウィナーフィルタ、および本発明の実施形態による補償フィルタのインパルス応答の関係を示した説明図である。これらのインパルス応答を有するＰＳＦ、ウィナーフィルタ、および本発明の実施形態による補償フィルタは、図４に示す周波数応答を有する。
【００４７】
図４に示したように、ウィナーフィルタは、周波数Ｓ２までは信号を増幅し、周波数Ｓ２よりも高域側は信号を減衰する。一方、本発明の実施形態による補償フィルタは、周波数Ｓ１までゲインが増加し、周波数Ｓ１よりも高域側はゲインが減少して０ｄＢに漸近する。すなわち、減衰量が補償範囲を超えて大きくなると周波数応答が１（＝０ｄＢ）に漸近する。
【００４８】
なお、本発明の実施形態による補償フィルタによるゲインが増加傾向から減少傾向に転じる周波数Ｓ１は、数式４におけるＳＮ比の周波数応答Ｎ／Ｓに依存する。例えば、ＳＮ比の周波数応答Ｎ／Ｓが小さいほど周波数Ｓ１は高くなり、ＳＮ比の周波数応答Ｎ／Ｓが大きいほど周波数Ｓ１は低くなる。
【００４９】
図５は、本発明の実施形態による補償フィルタによる出力の周波数特性を示した説明図である。図５に示したように、ＭＴＦによる実際の減衰が図４に推定結果として示した減衰よりも弱く、高域成分が予想以上に残っている場合でも、ウィナーフィルタの出力では高域成分が大きく減衰されてしまう。一方、図５に示したように、本発明の実施形態による補償フィルタは、高域成分を素通しするので、他の周波数成分の減衰を補償しつつ、高域成分を維持することができる。以下、図６〜図８を参照し、本発明の実施形態による補償フィルタの効果を具体的な画像例を用いてより詳細に説明する。
【００５０】
図６は、補償前の入力画像の具体例を示した説明図である。図７は、図６に示した入力画像にウィナーフィルタを適用した場合の出力画像を示した説明図である。図８は、図６に示した入力画像に本発明の実施形態による補償フィルタを適用した場合の出力画像を示した説明図である。
【００５１】
図６に示したように、補償前の入力画像は、撮像光学系１０のＭＴＦ特性により全体的にボケている。この入力画像にウィナーフィルタを適用すると、図７に示したように、入力画像のボケが全体的にある程度解消される。しかし、上述したようにウィナーフィルタは高域成分を減衰させてしまうので、図７において点線で囲った領域のように、高域成分が強いズボンの縦縞模様がつぶれてしまう。
【００５２】
一方、図８に示したように、本発明の実施形態による補償フィルタによれば、入力画像のボケを解消しつつ、図８において点線で囲った領域のように、高域成分が強いズボンの縦縞模様を鮮明な状態で出力することができる。
【００５３】
＜４．本発明の実施形態の動作＞
以上、本発明の実施形態の構成および効果を説明した。続いて、図９を参照し、本発明の実施形態の動作を簡単にまとめる。
【００５４】
図９は、本発明の実施形態の動作を示したシーケンス図である。図９に示したように、まず、フィルタ設計装置４０のＰＳＦ推定部４２が画像処理装置１の撮像光学系１０のＰＳＦを推定し（Ｓ１１０）、ＰＳＦ推定部４２により推定されたＰＳＦをフーリエ変換部４４がフーリエ変換する（Ｓ１２０）。
【００５５】
続いて、周波数応答設計部４６が、フーリエ変換部４４によるフーリエ変換結果であるＯＴＦ（ＰＳＦの周波数応答）を用い、数式４に基づく演算により補償フィルタの周波数応答を設計する（Ｓ１３０）。そして、逆フーリエ変換部４８が、周波数応答設計部４６により設計された周波数応答を、逆フーリエ変換によりインパルス応答に変換する（Ｓ１４０）。その後、逆フーリエ変換部４８により得られた補償フィルタのフィルタ係数を画像処理装置１のフィルタ係数保持部３０に格納する（Ｓ１５０）。
【００５６】
次に、画像処置装置１の撮像光学系１０およびＣＭＯＳ１４などにより被写体が撮像されて信号処理部２０に画像が入力されると（Ｓ１６０）、信号処理部２０は、フィルタ係数保持部３０に保持されているフィルタ係数からなる補償フィルタを入力画像に畳み込むことにより、撮像光学系１０のＭＴＦ特性による像の劣化を補償する（Ｓ１７０）。
【００５７】
＜５．むすび＞
以上説明したように、本発明の実施形態による補償フィルタの周波数応答は、推定されたＯＴＦが０に近づくにつれて１に近づく。すなわち、推定されたＯＴＦによる減衰が多い帯域は、ウィナーフィルタを用いるとさらに減衰されるのに対し、本発明の実施形態による補償フィルタを用いると素通しになる。したがって、本発明の実施形態によれば、推定されたＯＴＦによる信号劣化ほど実際の高域成分が失われていなかった場合、高域成分を減衰させることなく他の周波数成分を補償することが可能である。また、本発明の実施形態によれば、現実より大きいボケを仮定した場合にも必要以上に高域が減衰されることを防止できる。
【００５８】
なお、添付図面を参照しながら本発明の好適な実施形態について詳細に説明したが、本発明はかかる例に限定されない。本発明の属する技術の分野における通常の知識を有する者であれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本発明の技術的範囲に属するものと了解される。
【００５９】
例えば、上記ではＰＳＦが畳み込まれた入力画像を補償する画像処理装置を説明したが、本発明はかかる例に限定されない。ＰＳＦはインパルス応答の一例に過ぎず、入力画像はデータ信号の一例に過ぎず、画像処理装置はインパルス応答が畳み込まれたデータ信号を補償する信号処理装置の一例に過ぎない。本発明は、インパルス応答が畳み込まれたデータ信号を補償する信号処理装置全般に適用可能である。
【００６０】
また、本明細書の画像処理装置１、およびフィルタ設計装置４０の処理における各ステップは、必ずしもシーケンス図として記載された順序に沿って時系列に処理する必要はない。例えば、画像処理装置１、およびフィルタ設計装置４０の処理における各ステップは、シーケンス図として記載した順序と異なる順序で処理されても、並列的に処理されてもよい。
【００６１】
また、画像処理装置１、およびフィルタ設計装置４０に内蔵されるＣＰＵ、ＲＯＭおよびＲＡＭなどのハードウェアを、上述した受信装置２０、無線通信装置１０および２０の各構成と同等の機能を発揮させるためのコンピュータプログラムも作成可能である。また、該コンピュータプログラムを記憶させた記憶媒体も提供される。
【符号の説明】
【００６２】
１      画像処理装置
１０    撮像光学系
１４    ＣＭＯＳ
２０    信号処理部
３０    フィルタ係数保持部
４０    フィルタ設計部
４２    ＰＳＦ推定部
４４    フーリエ変換部
４６    周波数応答設計部
４８    逆フーリエ変換部
【図１】
2012141726.tif 000010
【図２】
2012141726.tif 000011
【図３】
2012141726.tif 000012
【図４】
2012141726.tif 000013
【図５】
2012141726.tif 000014
【図９】
2012141726.tif 000015
【図６】
2012141726.tif 000016
【図７】
2012141726.tif 000017
【図８】
2012141726.tif 000018