JP2019169000A - 画像処理装置、画像処理方法、およびプログラム - Google Patents
画像処理装置、画像処理方法、およびプログラム Download PDFInfo
- Publication number
- JP2019169000A JP2019169000A JP2018057092A JP2018057092A JP2019169000A JP 2019169000 A JP2019169000 A JP 2019169000A JP 2018057092 A JP2018057092 A JP 2018057092A JP 2018057092 A JP2018057092 A JP 2018057092A JP 2019169000 A JP2019169000 A JP 2019169000A
- Authority
- JP
- Japan
- Prior art keywords
- image data
- estimated
- processing apparatus
- phase
- image processing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000012545 processing Methods 0.000 title claims abstract description 83
- 238000003672 processing method Methods 0.000 title claims 3
- 238000012937 correction Methods 0.000 claims abstract description 75
- 238000001228 spectrum Methods 0.000 claims abstract description 32
- 238000006243 chemical reaction Methods 0.000 claims abstract description 4
- 238000006073 displacement reaction Methods 0.000 claims description 7
- 230000002194 synthesizing effect Effects 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims 2
- 230000002123 temporal effect Effects 0.000 abstract 1
- 238000000034 method Methods 0.000 description 43
- 238000003384 imaging method Methods 0.000 description 20
- 238000010586 diagram Methods 0.000 description 15
- 230000006870 function Effects 0.000 description 9
- 238000013519 translation Methods 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000015556 catabolic process Effects 0.000 description 2
- 238000005314 correlation function Methods 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
Images
Classifications
-
- G06T5/73—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration by the use of more than one image, e.g. averaging, subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/10—Image enhancement or restoration by non-spatial domain filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20224—Image subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30232—Surveillance
Abstract
【課題】ブロック処理を行うゆらぎ補正において、補正ずれを抑制すること。【解決手段】時間的に連続して撮像された複数の画像データを用いて大気のゆらぎによる影響を補正する画像処理装置を提供する。画像処理装置は、複数の画像データの処理対象の領域の画像に対応する複数の部分画像データを取得する取得手段と、取得された複数の部分画像データのそれぞれをフーリエスペクトルに変換する変換手段と、変換された複数のフーリエスペクトルの低周波成分の位相を用いて高周波成分を含む位相を推定する推定手段と、推定手段によって推定された結果と、取得手段によって取得された部分画像データに基づいて決定される基準データとに基づいて、推定手段によって推定された少なくとも一つの部分画像データに対応する結果を修正する修正手段と、修正手段によって修正された結果に基づいて画像データに対応する出力画像データを生成する生成手段とを備える。【選択図】図3
Description
本発明は、大気のゆらぎに起因する画像への影響を補正する技術に関する。
従来、監視カメラシステムなどにおいて、大気のゆらぎに起因する画像の劣化を補正する技術(以下、「ゆらぎ補正」と呼ぶ)が用いられている。ここで大気ゆらぎとは、大気の乱流や温度差により、屈折率が空間的および時間的に変化する自然現象を指す。大気ゆらぎは陽炎とも呼ばれる。大気ゆらぎは光の屈折や位相のばらつきを引き起こし、その結果として、特に望遠レンズで遠方を撮像した場合に、被写体の像に時間的に変化するぼけ及び歪みが発生する。このような画像の劣化は、被写体の視認性を大きく低下させるため、除去されることが望ましい。
複数フレームの画像から、ぼけ及び歪みをまとめて低減した画像を生成(再構成)する技術として、スペックルイメージングが知られている。以降の説明では、動画像を構成する時系列の画像データをフレームと呼ぶ。スペックルイメージングは、元来天文分野で開発された手法であり、短時間露光された画像からオリジナルの画像を再構成する技術である。スペックルイメージングにおいては、再構成される画像のフーリエスペクトルの振幅成分と位相成分とを独立に推定する処理が行われる。そして、推定された振幅成分と位相成分とを用いて逆フーリエ変換をすることで、ぼけ及び歪みを低減した画像が再構成される。特許文献1では、スペックルイメージングをゆらぎ補正に応用する技術が開示されている。
また、スペックルイメージングを行う際には、画像の領域ごとの部分画像に対してスペックルイメージングを行い、領域ごとに得られた部分画像を合成して最終出力画像を生成する処理(以下、「ブロック処理」と呼ぶ)が行われる。特許文献1においてもブロック処理を行う例が示されている。
しかしながら、特許文献1の技術では、スペックルイメージングによって生成される一部の部分画像に位置ずれなどの補正ずれが生じることがある。すると、それらを合成して得られる最終出力画像において特有のアーティファクトが生じてしまい、視認性が低下する。
本発明は、ブロック処理を行うゆらぎ補正において、補正ずれを抑制することを目的とする。
本発明の一態様に係る画像処理装置は、時間的に連続して撮像された複数の画像データを用いて大気のゆらぎによる影響を補正する画像処理装置であって、前記複数の画像データの処理対象の領域の画像に対応する複数の部分画像データを取得する取得手段と、前記取得された複数の部分画像データのそれぞれをフーリエスペクトルに変換する変換手段と、前記変換された複数のフーリエスペクトルの低周波成分の位相を用いて高周波成分を含む位相を推定する推定手段と、前記推定手段によって推定された結果と、前記取得手段によって取得された部分画像データに基づいて決定される基準データとに基づいて、前記推定手段によって推定された少なくとも一つの部分画像データに対応する結果を修正する修正手段と、前記修正手段によって修正された結果に基づいて前記画像データに対応する出力画像データを生成する生成手段とを備えることを特徴とする。
本発明によれば、ブロック処理を行うゆらぎ補正において、補正ずれを抑制することができる。
以下、本発明の実施形態について、図面を参照して説明する。なお、以下の実施形態は本発明を限定するものではなく、また、本実施形態で説明されている特徴の組み合わせの全てが本発明の解決手段に必須のものとは限らない。なお、同一の構成については、同じ符号を付して説明する。
<<実施形態1>>
本実施形態では、スペックルイメージングの過程で導出される位相データを用いて、ゆらぎ補正によって補正ずれが生じているかを判定する。本実施形態では、補正ずれの一例として画像位置の平行移動ずれ(位置ずれ)を扱う。そして、判定結果に基づいて、ゆらぎ補正によって生じる位置ずれを修正する例を述べる。
本実施形態では、スペックルイメージングの過程で導出される位相データを用いて、ゆらぎ補正によって補正ずれが生じているかを判定する。本実施形態では、補正ずれの一例として画像位置の平行移動ずれ(位置ずれ)を扱う。そして、判定結果に基づいて、ゆらぎ補正によって生じる位置ずれを修正する例を述べる。
<画像処理装置のハードウェア構成>
図1は、本実施形態の画像処理装置のハードウェア構成の一例を示す図である。本実施形態の画像処理装置100は、CPU101、RAM102、ROM103、二次記憶装置104、入力インターフェース105、および出力インターフェース106を含む。画像処理装置100の各構成要素は、システムバス107によって相互に接続されている。画像処理装置100は、入力インターフェース105を介して外部記憶装置108および操作部110に接続されている。画像処理装置100は、出力インターフェース106を介して外部記憶装置108および表示装置109に接続されている。
図1は、本実施形態の画像処理装置のハードウェア構成の一例を示す図である。本実施形態の画像処理装置100は、CPU101、RAM102、ROM103、二次記憶装置104、入力インターフェース105、および出力インターフェース106を含む。画像処理装置100の各構成要素は、システムバス107によって相互に接続されている。画像処理装置100は、入力インターフェース105を介して外部記憶装置108および操作部110に接続されている。画像処理装置100は、出力インターフェース106を介して外部記憶装置108および表示装置109に接続されている。
CPU101は、RAM102をワークメモリとして、ROM103に格納されたプログラムを実行し、システムバス107を介して画像処理装置100の各構成要素を統括的に制御する。これにより、後述する様々な処理が実行される。二次記憶装置104は、画像処理装置100で取り扱われる種々のデータを記憶する記憶装置であり、本実施形態ではHDDが用いられる。CPU101は、システムバス107を介して二次記憶装置104へのデータの書き込みおよび二次記憶装置104に記憶されたデータの読出しを行う。なお、二次記憶装置104にはHDDの他に、光ディスクドライブやフラッシュメモリなど、様々な記憶デバイスを用いることが可能である。
入力インターフェース105は、例えばUSBやIEEE1394等のシリアルバスインターフェースである。画像処理装置100は、入力インターフェース105を介して、外部装置からデータや命令等を入力する。本実施形態では、画像処理装置100は、入力インターフェース105を介して、外部記憶装置108(例えば、ハードディスク、メモリカード、CFカード、SDカード、USBメモリなどの記憶媒体)からデータを取得する。また本実施形態では、画像処理装置100は、操作部110に入力されたユーザの指示を、入力インターフェース105を介して取得する。操作部110は、マウスやキーボードなどの入力装置であり、ユーザの指示を入力する。
出力インターフェース106は、入力インターフェース105と同様にUSBやIEEE1394等のシリアルバスインターフェースである。なお、出力インターフェース106は、例えばDVIやHDMI(登録商標)等の映像出力端子であってもよい。画像処理装置100は、出力インターフェース106を介して、外部装置にデータ等を出力する。本実施形態では、画像処理装置100は、出力インターフェース106を介して表示装置109(液晶ディスプレイなどの各種画像表示デバイス)に、CPU101によって処理されたデータ(例えば、画像データ)を出力する。なお、画像処理装置100の構成要素は上記以外にも存在し得るが、ここでは説明を省略する。
<画像処理装置の機能ブロック図>
図2は、画像処理装置100の機能ブロックの例を示す図である。画像処理装置100は、画像データ取得部201、画像分割部202、ゆらぎ補正部203、および画像結合部204を有する。
図2は、画像処理装置100の機能ブロックの例を示す図である。画像処理装置100は、画像データ取得部201、画像分割部202、ゆらぎ補正部203、および画像結合部204を有する。
図3は、ゆらぎ補正部203の詳細を示すブロック図である。ゆらぎ補正部203は、フーリエ変換部301、振幅推定部302、位相推定部303、位置補償部304、および逆フーリエ変換部305を有する。位置補償部304は、基準位相決定部351、位置ずれ判定部352、および位置ずれ修正部353を有する。各部で行われる処理の説明については、後述するフローチャートの説明とともに行うこととする。
図2および図3に示す各部の機能は、CPU101がROM103内に格納された制御プログラムを読み込み実行することで実現される。なお、各部に相当する専用の処理回路を備えるように画像処理装置100を構成するようにしてもよい。
<一般的なゆらぎ補正の説明>
図4は、一般的なゆらぎ補正の一例を説明する図である。先に説明したように、複数の画像データ(フレーム画像)を用いて、大気のゆらぎによる、画像のぼけ及び歪み等をまとめて低減した画像を生成することが行われている。図4では、複数のフレーム401(入力画像データ)が画像処理装置100に入力された場合に、1枚のゆらぎ補正済みの出力画像データ441が得られる例を示している。ゆらぎ補正を行う手法として、「スペックルイメージング処理」が知られている。
図4は、一般的なゆらぎ補正の一例を説明する図である。先に説明したように、複数の画像データ(フレーム画像)を用いて、大気のゆらぎによる、画像のぼけ及び歪み等をまとめて低減した画像を生成することが行われている。図4では、複数のフレーム401(入力画像データ)が画像処理装置100に入力された場合に、1枚のゆらぎ補正済みの出力画像データ441が得られる例を示している。ゆらぎ補正を行う手法として、「スペックルイメージング処理」が知られている。
ゆらぎ補正処理においては、フレーム401の処理対象の領域に対応する部分画像をフレームごとに切り出して、複数の部分画像データ402が取得される。そして、複数の部分画像データ402をそれぞれフーリエ変換して、周波数空間の成分(位相成分411と振幅成分421)にそれぞれ変換する。そして、周波数空間において、再構成される画像の位相成分412と振幅成分422とをそれぞれ独立して推定する処理が行われる。振幅成分422は、例えば入力された各部分画像のパワースペクトルの平均を導出し、導出された平均に対し撮像装置の光学系や大気の影響による劣化モデルに基づいて決定された高周波成分を増幅するフィルタをかけることで推定される。位相成分412は、例えば入力された各部分画像のバイスペクトラムに基づき、既知の低周波成分の位相を用いて高周波成分の位相を逐次的に推定していくことで求められる。位相成分412の推定処理の詳細は後述する。
このように推定された位相成分412と振幅成分422とを逆フーリエ変換(逆変換)することで、ゆらぎ補正済みの部分画像データ431が得られる。このような処理を、フレームの各領域を処理対象の領域として処理することで、各領域に対応するゆらぎ補正済み部分画像データが得られることになる。そして、ゆらぎ補正済み部分画像データ431を合成することで、ゆらぎ補正がされた出力画像データ441が得られることになる。
このように、部分画像データからぼけ及び歪みを低減した後に、各部分画像データを合成して、出力画像データを生成する処理(ブロック処理)が行われる。このとき一部の部分画像において、スペックルイメージングによって生成される画像に位置ずれなどの補正ずれが生じることがある。すると、それらを合成して得られる出力画像データにおいて特有のアーティファクトが生じてしまい、視認性が低下してしまうことがある。
図5は、ゆらぎ補正と位置ずれを説明する図である。図5(a)は、理想的な画像510を示している。画像510は、周波数成分515−1、515−2、515−3・・・から構成されている。しかしながら、実際には、大気のゆらぎの影響によって、撮像された画像は、図5(b)の画像520となってしまっている。大気のゆらぎの影響を受けた画像520は、周波数毎に異なる方向で異なる量の位置ずれが生じている。このため、異なる方向で異なる量の位置ずれが生じている周波数成分525−1、525−2、525−3・・・によって構成される画像520の境界がボケて見える。図4に示すゆらぎ補正処理(スペックルイメージング処理)は、大気のゆらぎの影響を補正して理想的な画像510を得るための処理である。
しかしながら、図4に示すようなゆらぎ補正処理(スペックルイメージング処理)を行った場合、図5(c)の周波数成分535−1、535−2、535−3・・・のように、全ての周波数に同じ方向で同じ量の位置ずれが残存することがある。これらの周波数成分535−1、535−2、535−3・・・によって構成される画像530のボケは解消されているが、理想的な画像510に比べて位置ずれが生じてしまっている。
本実施形態では、このような位置ずれを修正する処理を行うものである。より詳細には、本実施形態では、図4に示すようなスペックルイメージングの過程で導出される位相成分411を用いて位置ずれ判定を行い、判定結果に基づいて位置ずれ修正を行う形態を説明する。
<本実施形態のゆらぎ補正の説明>
図6は、本実施形態で説明するゆらぎ補正における位置ずれを修正する処理を説明する図である。部分画像データ402をフーリエ変換する処理は、図4の例と同様である。本実施形態では、スペックルイメージングの過程で導出される位相成分411を用いて位置ずれ判定を行い、判定結果に基づいて位置ずれ修正を行う形態である。よって、図6において周波数空間の成分は、位相成分のみを図示している。振幅成分については、図4で説明した例と同様であり、図示を省略している。
図6は、本実施形態で説明するゆらぎ補正における位置ずれを修正する処理を説明する図である。部分画像データ402をフーリエ変換する処理は、図4の例と同様である。本実施形態では、スペックルイメージングの過程で導出される位相成分411を用いて位置ずれ判定を行い、判定結果に基づいて位置ずれ修正を行う形態である。よって、図6において周波数空間の成分は、位相成分のみを図示している。振幅成分については、図4で説明した例と同様であり、図示を省略している。
本実施形態では、ゆらぎ補正部203は、部分画像データの位相成分411を用いて、位置ずれの判定に用いられる基準データ(基準位相成分610)を決定する。ゆらぎ補正部203は、推定された位相成分412と、基準位相成分610とを比較して位置ずれ量を取得する。そして、位置ずれ量が所定の閾値より大きい場合、位置ずれが生じていると判定し、ゆらぎ補正部203は、位置ずれを修正した位相成分615を導出する。ゆらぎ補正部203は、この位置ずれを修正した位相成分615と、推定された振幅成分422(図6では不図示)とを逆フーリエ変換して、位置ずれを修正した、ゆらぎ補正済みの部分画像データ620を生成する。そして、各注目領域に対応する部分画像データ620を合成して、出力画像データ630が生成される。以下、本実施形態の処理を詳細に説明する。
<全体のフローチャート>
図7は、本実施形態における処理の流れを示すフローチャートである。図7に示す処理は、CPU101が、ROM103内に格納された制御プログラムを読み込んでRAM102に展開し、これをCPU101が実行することで実現される。あるいはまた、図7におけるステップの一部または全部の機能をASICや電子回路等のハードウェアで実現してもよい。なお、各処理の説明における記号「S」は、当該フローチャートにおけるステップであることを意味する(本明細書において同様である)。以下、図2と図7とを用いて説明する。
図7は、本実施形態における処理の流れを示すフローチャートである。図7に示す処理は、CPU101が、ROM103内に格納された制御プログラムを読み込んでRAM102に展開し、これをCPU101が実行することで実現される。あるいはまた、図7におけるステップの一部または全部の機能をASICや電子回路等のハードウェアで実現してもよい。なお、各処理の説明における記号「S」は、当該フローチャートにおけるステップであることを意味する(本明細書において同様である)。以下、図2と図7とを用いて説明する。
S701において画像データ取得部201は、ゆらぎ補正を行う処理対象である画像データを取得する。本実施形態では、動画像を構成する時間的に連続した複数のフレームを取得する。時間的に連続した複数の静止画像のデータを取得する形態でもよい。
S702において画像分割部202は、S701で取得した画像データを複数の小領域データ(部分画像データ)に分割する。即ち、各画像データを、複数の部分画像データに分割する。ここで複数の部分画像データは、領域の重複が許容される。部分画像データのサイズは、大気ゆらぎの影響が一定であるとみなせる範囲(アイソプラナティックパッチ)よりも小さいことが好ましい。S702の処理により、S701で取得された複数の画像データのそれぞれが、各小領域に対応する部分画像データに分割される。以降の処理では、小領域ごとに処理が行われる。
S703においてゆらぎ補正部203は、全ての小領域に対して処理をしたかを判定する。全ての小領域に対して処理をした場合、S706に進む。未処理の小領域がある場合、S704に進み、ゆらぎ補正部203は、未処理の小領域の中から処理対象とする小領域を決定する。
S705においてゆらぎ補正部203は、処理対象の小領域の部分画像データを取得し、ゆらぎ補正処理を行う。即ち、ゆらぎ補正部203は、処理対象の小領域に対応する、複数フレーム分の複数の部分画像データを用いてゆらぎ補正処理を行う。詳細は後述する。
S704において画像結合部204は、S705でゆらぎ補正された各小領域に対応する推定画像(部分画像)を合成して、出力画像を生成する。
<ゆらぎ補正処理のフローチャート>
次に、図3と図8とを用いてS705のゆらぎ補正処理の詳細を説明する。S801においてフーリエ変換部301は、各フレームの処理対象の小領域に対応する部分画像データのそれぞれに対してフーリエ変換を行ってフーリエスペクトル(周波数成分)を生成する。S801により、フレーム数分のフーリエスペクトルが得られる。
次に、図3と図8とを用いてS705のゆらぎ補正処理の詳細を説明する。S801においてフーリエ変換部301は、各フレームの処理対象の小領域に対応する部分画像データのそれぞれに対してフーリエ変換を行ってフーリエスペクトル(周波数成分)を生成する。S801により、フレーム数分のフーリエスペクトルが得られる。
S802において振幅推定部302は、S801で得られた各フレーム(各部分画像)のフーリエスペクトルを用いて、推定画像のフーリエスペクトルOの振幅を推定する。具体的には、例えば特許文献1に示されているように、式(1)に従いフーリエスペクトルOの振幅が導出される。
S803において位相推定部303は、S801で得られた各フレームのフーリエスペクトルを用いて、推定画像のフーリエスペクトルOの位相(推定位相φ)を推定する。特許文献1に示されている通り、この推定位相φは式(2)によって決定される。
式(2)によれば、uとvの位相を用いてu+vの位相を推定できる。例えば、u=(0,1)、v=(1,0)とすると、これらを用いてu+v=(1,1)の位相を推定することができる。一方、u=(0,1)の位相は−u=(0,−1)のマイナス1倍である。よって、uの位相が既知ということであれば−uの位相も既知となる。以上の点を考慮すると、(0,0)、(0,1)、(1,0)、(0,−1)、(−1,0)、の5個を初期位相として推定済みの位相とすることができる。
以降のステップでは、直前のステップまでで推定済みの位相から任意に2個を選択して推定を行い、推定済みの位相をさらに用いて以降のステップで用いる。例えば、以下の推定が可能である。
(0,1)、(0,1)から(0,2)の位相を推定
(0,1)、(1,0)から(1,1)の位相を推定
(0,1)、(0,−1)から(0,0)の位相を推定
(0,1)、(−1,0)から(−1,1)の位相を推定
(1,0)、(1,0)から(2,0)の位相を推定
(1,0)、(0,−1)から(1,−1)の位相を推定
(1,0)、(−1,0)から(0,0)の位相を推定
(0,−1)、(0,−1)から(0,−2)の位相を推定
(0,−1)、(−1,0)から(−1,−1)の位相を推定
(−1,0)、(−1,0)から(−2,0)の位相を推定
このような処理を繰り返すことで、全ての周波数座標における位相が推定される。なお、座標ベクトルの絶対値が周波数の高低を示している。即ち、上述の例では、絶対値が小さい(低周波)座標から絶対値の大きい(高周波)座標の位相を推定していることになる。
(0,1)、(1,0)から(1,1)の位相を推定
(0,1)、(0,−1)から(0,0)の位相を推定
(0,1)、(−1,0)から(−1,1)の位相を推定
(1,0)、(1,0)から(2,0)の位相を推定
(1,0)、(0,−1)から(1,−1)の位相を推定
(1,0)、(−1,0)から(0,0)の位相を推定
(0,−1)、(0,−1)から(0,−2)の位相を推定
(0,−1)、(−1,0)から(−1,−1)の位相を推定
(−1,0)、(−1,0)から(−2,0)の位相を推定
このような処理を繰り返すことで、全ての周波数座標における位相が推定される。なお、座標ベクトルの絶対値が周波数の高低を示している。即ち、上述の例では、絶対値が小さい(低周波)座標から絶対値の大きい(高周波)座標の位相を推定していることになる。
上述した初期位相として、例えば複数フレームを平均した部分画像データのフーリエスペクトルの位相の値を用いることができる。なお、初期位相は上記と異なる3点に対して定めても良いし、周波数平面上の4点以上に対して定めても良い。
S804において位置補償部304は、S803で得られた推定位相φに対して位置ずれを判定し、判定結果に基づいて、位置ずれを修正する処理を行う。位置ずれを判定し、判定結果に基づいて、位置ずれを修正する処理のことを位置補償処理ともいう。位置補償処理の詳細は後述する。
S805において逆フーリエ変換部305は、推定画像のフーリエスペクトルOに対して逆フーリエ変換を行って推定画像を生成する。即ち、逆フーリエ変換部305は、振幅推定され、かつ位相推定(必要に応じて位置ずれ修正)がされたフーリエスペクトルOに対して逆フーリエ変換を行う。以上の処理により、処理対象の小領域に対応する部分画像データのゆらぎ補正が行われる。なお、図8の例では、振幅推定の後に位相推定を行うフローチャートを示しているが、この限りではない。位相推定の後に振幅推定を行っても良いし、振幅推定と位相推定とを並行して行っても良い。
図9は、S804の位置補償処理の詳細を示すフローチャートである。S901において基準位相決定部351は、S801で得られた各フレームのフーリエスペクトルを用いて、基準位相を決定する。基準位相は、例えば複数フレームの平均を用いることができる。この平均は、ゆらぎ補正で用いるフレーム全てについてでも良く、一部でも良い。また、任意の1フレームの位相でも良い。あるいは、平均値でなく、中央値など別の統計量で基準位相を求めても良い。また、別途の位置推定処理を行って基準画像を生成しても良い。本実施形態では、u=(u1、u2)とした場合、基準位相として、式(4)で導出されるnフレーム分の平均のφaveを用いる。
推定位相φと基準位相φaveとの変化量が所定の範囲に含まれる場合、例えば、位相差分が予め定める閾値を超える場合に、位置ずれを修正する理由を説明する。ゆらぎ補正は、元来、局所的な位置ずれを補正するものである。つまり、位置ずれが小さい場合とは、本来的なゆらぎ補正をした結果のものと考えられる。従って、位相差分が所定の閾値を超えていない場合には、位置ずれを修正しない。また、ごく小さい位置ずれであれば、見た目がほとんど変わらないので、処理時間短縮のために、位置ずれを修正しない場合もある。また、1画素以下の位置ずれを修正しようとすると、補間処理によるボケの様な弊害が生じる場合もある。よって、本実施形態では、所定の閾値を超える場合に位置ずれを修正する処理を行っている。
ここまで簡単のために1次元のデータで説明したが、実際には2次元のデータに対して本処理は行われる。すなわち、位置ずれ判定部352は、位相差分を平面(を2πで割った余り)で近似し、近似した平面の傾きを変化量として算出し、傾きが所定の範囲に含まれるかを判定する。例えば、位置ずれ判定部352は、位相差分を平面(を2πで割った余り)で近似し、その平面を特徴づけるパラメータを閾値と比較して位置ずれが生じているか否かを判定する。具体例として、位相差分を適当なパラメータAとBを用いて式(6)で近似し、このときのAとBの絶対値の内で最小のものと閾値とを比較して判定を行うなどの方法がある。
以上説明したように、本実施形態では、スペックルイメージングにおいて推定された位相の位置ずれが生じている場合に、周波数空間において位相成分を修正することで、位置ずれを修正する処理を行う。このような処理によれば、大気のゆらぎを低減した画像を生成し、かつ、位置ずれが生じている場合には位置ずれを修正した画像を生成することができる。
<<実施形態2>>
実施形態1では、周波数空間において、推定された位相成分と基準位相との差分に基づいて位置ずれが生じているかを判定し、位置ずれが生じている場合に、周波数空間において位置ずれを修正する形態を説明した。
実施形態1では、周波数空間において、推定された位相成分と基準位相との差分に基づいて位置ずれが生じているかを判定し、位置ずれが生じている場合に、周波数空間において位置ずれを修正する形態を説明した。
本実施形態では、スペックルイメージングの結果導出される画像データ(即ち、逆フーリエ変換した画像データ)を用いて位置ずれ判定を行う。そして、判定結果に基づいて実空間において位置ずれ修正を行う形態を説明する。
図11は、本実施形態の位置ずれ修正の概要を説明する図である。複数の部分画像データ402からゆらぎ補正済み部分画像データ431が得られる処理については図4に示す例と同様である。本実施形態では、部分画像データ402から基準データとなる画像データ(基準画像データ1110という)が決定される。そして、ゆらぎ補正済み部分画像データ431と基準画像データ1110とを比較して位置ずれ量が取得される。そして、位置ずれ量が所定の閾値より大きい場合、位置ずれが修正される。このように位置ずれが修正されたゆらぎ補正済みの部分画像データ1120を用いて出力画像データ1130が生成される。
<ゆらぎ補正部のブロック図>
図12は、本実施形態のゆらぎ補正部203のブロック図である。本実施形態のゆらぎ補正部203は、実施形態1のゆらぎ補正部と同様の機能を有するが、各機能部の接続が異なる。具体的には、位置補償部1210の位置が異なる。また、位置補償部1210が有する機能も異なる。位置補償部1210は、基準画像決定部1211、位置ずれ判定部1212、位置ずれ修正部1213としての機能を有する。
図12は、本実施形態のゆらぎ補正部203のブロック図である。本実施形態のゆらぎ補正部203は、実施形態1のゆらぎ補正部と同様の機能を有するが、各機能部の接続が異なる。具体的には、位置補償部1210の位置が異なる。また、位置補償部1210が有する機能も異なる。位置補償部1210は、基準画像決定部1211、位置ずれ判定部1212、位置ずれ修正部1213としての機能を有する。
図13は、本実施形態における処理の流れを示すフローチャートである。S1301、S1302、S1303は、S801、S802、S803と同等の処理である。S1304において逆フーリエ変換部305は、推定画像のフーリエスペクトルに対して逆フーリエ変換を行って推定画像を生成する。S1305において位置補償部1210は、推定画像(部分画像)に対して位置補償を行う。
図14は、S1305の位置補償処理の詳細を説明する図である。S1401において基準画像決定部1211は、S702で得られた、処理対象の小領域に対応する部分画像データを用いて、基準画像を決定する。本実施形態では、基準画像として、複数フレーム分の部分画像データのフレーム平均を用いる。なお、基準画像は、例えばゆらぎ補正で用いるフレーム全てについてでも良く、一部でも良い。また、任意の1フレームの画像でも良い。あるいは、平均値でなく、中央値など別の統計量で基準画像を求めても良い。
S1402において位置ずれ判定部1212は、S1304で得られた推定画像とS1401で得られた基準画像とを用いて位置ずれ判定を行う。位置ずれ判定の一例を示す。まず、S1304で得られた推定画像とS1401で得られた基準画像との相互相関関数を最大にする平行移動量dを導出する。そして位置ずれ判定部1212は、dが予め定める所定の範囲に含まれる場合、例えば予め定める閾値より大きい場合、位置ずれが生じていると判定する。そうでない場合、位置ずれが生じていないと判定する。
S1403において位置ずれ修正部353は、S1402の判定の結果、平行移動量が閾値よりも大きい場合には、S1404に進み、位置ずれ修正処理を行う。この処理は、推定画像を平行移動量dだけ平行移動することで実現できる。平行移動量が閾値よりも大きくない場合には、S1404の処理はスキップする。
以上説明したように、本実施形態によれば、大気のゆらぎを低減した画像を生成し、かつ、位置ずれが生じている場合には位置ずれを修正した画像を生成することができる。
<<その他の実施形態>>
上述した実施形態においては、位置ずれの例として平行移動ずれの場合を例に挙げて説明したが、この限りではない。他の位置ずれを修正することもできる。例として、回転(または拡縮)ずれを補正してもよい。この場合、ステップS1402で推定画像と基準画像の相互相関関数を最大にする回転量r(または拡縮量e)を求め、ステップS1404で回転量r(または拡縮量e)だけ回転(または拡大もしくは縮小)させれば良い。
上述した実施形態においては、位置ずれの例として平行移動ずれの場合を例に挙げて説明したが、この限りではない。他の位置ずれを修正することもできる。例として、回転(または拡縮)ずれを補正してもよい。この場合、ステップS1402で推定画像と基準画像の相互相関関数を最大にする回転量r(または拡縮量e)を求め、ステップS1404で回転量r(または拡縮量e)だけ回転(または拡大もしくは縮小)させれば良い。
本発明は、上述の実施形態の1以上の機能を実現するプログラムを、ネットワーク又は記憶媒体を介してシステム又は装置に供給し、そのシステム又は装置のコンピュータにおける1つ以上のプロセッサーがプログラムを読出し実行する処理でも実現可能である。また、1以上の機能を実現する回路(例えば、ASIC)によっても実現可能である。
301 フーリエ変換部
302 振幅推定部
303 位相推定部
304 位置補償部
305 逆フーリエ変換部
302 振幅推定部
303 位相推定部
304 位置補償部
305 逆フーリエ変換部
Claims (14)
- 時間的に連続して撮像された複数の画像データを用いて大気のゆらぎによる影響を補正する画像処理装置であって、
前記複数の画像データの処理対象の領域の画像に対応する複数の部分画像データを取得する取得手段と、
前記取得された複数の部分画像データのそれぞれをフーリエスペクトルに変換する変換手段と、
前記変換された複数のフーリエスペクトルの低周波成分の位相を用いて高周波成分を含む位相を推定する推定手段と、
前記推定手段によって推定された結果と、前記取得手段によって取得された部分画像データに基づいて決定される基準データとに基づいて、前記推定手段によって推定された少なくとも一つの部分画像データに対応する結果を修正する修正手段と、
前記修正手段によって修正された結果に基づいて前記画像データに対応する出力画像データを生成する生成手段と
を備えることを特徴とする画像処理装置。 - 前記修正手段は、前記推定手段によって推定された結果と、前記取得手段によって取得された部分画像データに基づいて決定される基準データとの変化量が所定の範囲に含まれる場合、前記推定手段によって推定された結果を修正することを特徴とする請求項1に記載の画像処理装置。
- 前記修正手段は、前記推定手段によって推定された位相を修正し、
前記生成手段は、
前記位相が修正されたフーリエスペクトルを、前記処理対象の領域に対応する部分画像データに逆変換し、
各領域に対応する部分画像データを合成することで前記出力画像データを生成することを特徴とする請求項2に記載の画像処理装置。 - 前記修正手段は、
前記変換手段によって変換された複数のフーリエスペクトルの位相から前記基準データを決定し、
前記推定された位相と前記基準データとの差を平面で近似し、近似した平面の傾きを前記変化量として算出し、前記傾きが前記所定の範囲に含まれるかを判定することを特徴とする請求項3に記載の画像処理装置。 - 前記修正手段は、前記推定された位相から前記近似した平面を引くことで前記修正を行うことを特徴とする請求項4に記載の画像処理装置。
- 前記推定手段は、前記フーリエスペクトルの振幅をさらに推定し、
前記生成手段は、前記推定された振幅を含み、かつ前記位相が修正されたフーリエスペクトルを逆変換することを特徴とする請求項3から5のいずれか一項に記載の画像処理装置。 - 前記推定手段は、前記フーリエスペクトルの振幅をさらに推定し、
前記修正手段は、
前記推定された位相および振幅を含むフーリエスペクトルを前記処理対象の領域に対応する部分画像データに逆変換し、
前記逆変換した部分画像データを修正することを特徴とする請求項2に記載の画像処理装置。 - 前記修正手段は、
前記複数の部分画像データから前記基準データを決定し、
前記逆変換した部分画像データと前記基準データとから算出される変化量が前記所定の範囲に含まれるかを判定することを特徴とする請求項7に記載の画像処理装置。 - 前記生成手段は、前記修正された部分画像データを含む各領域に対応する部分画像データを合成することで前記出力画像データを生成することを特徴とする請求項7または8に記載の画像処理装置。
- 前記推定手段によって推定された結果と前記基準データとの変化量は、画像の平行移動ずれに相当することを特徴とする請求項2から9のいずれか一項に記載の画像処理装置。
- 前記推定手段によって推定された結果と前記基準データとの変化量は、画像の回転ずれに相当することを特徴とする請求項2から9のいずれか一項に記載の画像処理装置。
- 前記推定手段によって推定された結果と前記基準データとの変化量は、画像の拡大ずれまたは縮小ずれに相当することを特徴とする請求項2から9のいずれか一項に記載の画像処理装置。
- 時間的に連続して撮像された複数の画像データを用いて大気のゆらぎによる影響を補正する画像処理方法であって、
前記複数の画像データの処理対象の領域の画像に対応する複数の部分画像データを取得する取得ステップと、
前記取得された複数の部分画像データのそれぞれをフーリエスペクトルに変換する変換ステップと、
前記変換された複数のフーリエスペクトルの低周波成分の位相を用いて高周波成分を含む位相を推定する推定ステップと、
前記推定ステップにおいて推定された結果と、前記取得ステップにおいて取得された部分画像データに基づいて決定される基準データとに基づいて、前記推定ステップにおいて推定された少なくとも一つの部分画像データに対応する結果を修正する修正ステップと、
前記修正ステップにおいて修正された結果に基づいて前記画像データに対応する出力画像データを生成する生成ステップと
を備えることを特徴とする画像処理方法。 - コンピュータを、請求項1から12のいずれか一項に記載の画像処理装置の各手段として機能させるためのプログラム。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018057092A JP2019169000A (ja) | 2018-03-23 | 2018-03-23 | 画像処理装置、画像処理方法、およびプログラム |
US16/351,197 US10846839B2 (en) | 2018-03-23 | 2019-03-12 | Image processing apparatus, image processing method, and storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018057092A JP2019169000A (ja) | 2018-03-23 | 2018-03-23 | 画像処理装置、画像処理方法、およびプログラム |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2019169000A true JP2019169000A (ja) | 2019-10-03 |
Family
ID=67985462
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2018057092A Pending JP2019169000A (ja) | 2018-03-23 | 2018-03-23 | 画像処理装置、画像処理方法、およびプログラム |
Country Status (2)
Country | Link |
---|---|
US (1) | US10846839B2 (ja) |
JP (1) | JP2019169000A (ja) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7292123B2 (ja) * | 2019-06-20 | 2023-06-16 | キヤノン株式会社 | 撮像装置及びその制御方法、プログラム、記憶媒体 |
JP7263149B2 (ja) * | 2019-06-26 | 2023-04-24 | キヤノン株式会社 | 画像処理装置、画像処理方法、およびプログラム |
JP2022156693A (ja) * | 2021-03-31 | 2022-10-14 | キヤノン株式会社 | 補正装置、画像読取装置、画像形成装置、それらの制御方法及びプログラム |
CN115713475B (zh) * | 2023-01-10 | 2023-05-05 | 深圳市格灵精睿视觉有限公司 | 图像处理方法、装置、设备及存储介质 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7245742B2 (en) * | 2002-07-01 | 2007-07-17 | The Regents Of The University Of California | Video surveillance with speckle imaging |
JPWO2011132279A1 (ja) * | 2010-04-21 | 2013-07-18 | キヤノン株式会社 | 画像処理装置、方法、及び記録媒体 |
WO2013130629A1 (en) * | 2012-02-27 | 2013-09-06 | Greywolf Technical Services, Inc. | Methods and systems for modified wavelength diversity image compensation |
JP6258209B2 (ja) * | 2012-09-27 | 2018-01-10 | パナソニックIpマネジメント株式会社 | ステレオ画像処理装置およびステレオ画像処理方法 |
CN106165392B (zh) * | 2014-03-31 | 2019-01-22 | 富士胶片株式会社 | 图像处理装置、摄影装置以及图像处理方法 |
JP2017010095A (ja) * | 2015-06-17 | 2017-01-12 | キヤノン株式会社 | 画像処理装置、撮像装置、画像処理方法、画像処理プログラム、および、記憶媒体 |
-
2018
- 2018-03-23 JP JP2018057092A patent/JP2019169000A/ja active Pending
-
2019
- 2019-03-12 US US16/351,197 patent/US10846839B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
US10846839B2 (en) | 2020-11-24 |
US20190295234A1 (en) | 2019-09-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9692939B2 (en) | Device, system, and method of blind deblurring and blind super-resolution utilizing internal patch recurrence | |
JP5385969B2 (ja) | 画像の超解像のための方法及び装置 | |
JP2019169000A (ja) | 画像処理装置、画像処理方法、およびプログラム | |
KR100907120B1 (ko) | 열화 정보 복원 방법, 복원 장치 및 프로그램이 기록된 기록 매체 | |
JP2017010092A (ja) | 画像処理装置、撮像装置、画像処理方法、画像処理プログラム、および、記憶媒体 | |
KR101671391B1 (ko) | 레이어 블러 모델에 기반한 비디오 디블러링 방법, 이를 수행하기 위한 기록 매체 및 장치 | |
US11145032B2 (en) | Image processing apparatus, method and storage medium for reducing color noise and false color | |
JP2017010095A (ja) | 画像処理装置、撮像装置、画像処理方法、画像処理プログラム、および、記憶媒体 | |
Jeong et al. | Multi-frame example-based super-resolution using locally directional self-similarity | |
Khattab et al. | Regularization-based multi-frame super-resolution: a systematic review | |
WO2016193393A1 (en) | Real-time temporal filtering and super-resolution of depth image sequences | |
JP6819794B2 (ja) | レーダ画像処理装置、レーダ画像処理方法およびレーダ画像処理プログラム | |
JP2017010093A (ja) | 画像処理装置、撮像装置、画像処理方法、画像処理プログラム、および、記憶媒体 | |
JPWO2008102898A1 (ja) | 画質改善処理装置、画質改善処理方法及び画質改善処理プログラム | |
Maiseli et al. | A robust super-resolution method with improved high-frequency components estimation and aliasing correction capabilities | |
KR101341617B1 (ko) | 단일 영상의 오차모델을 기반으로 한 고해상도 영상 복원장치 및 방법 | |
JPWO2016098323A1 (ja) | 情報処理装置、情報処理方法、及び、プログラム | |
KR101362183B1 (ko) | 카메라 자세 기반의 깊이 영상 노이즈 제거 장치 및 방법 | |
Al Ismaeil et al. | Real-time non-rigid multi-frame depth video super-resolution | |
JP6541454B2 (ja) | 画像処理装置、撮像装置、画像処理方法、画像処理プログラム、および、記憶媒体 | |
KR100805802B1 (ko) | 모션 블러가 심한 영상의 카메라 자동보정 장치 및 그방법과, 그를 이용한 증강현실 시스템 | |
Mohan | Adaptive super-resolution image reconstruction with lorentzian error norm | |
JP2012169701A (ja) | 画像処理装置、画像処理方法及びプログラム | |
El Abbadi et al. | A Survey on Blind De-Blurring of Digital Image. | |
JP2021005324A (ja) | 画像処理装置、画像処理方法、およびプログラム |