ローリングシャッタ方式で駆動する撮像素子を備える撮像装置であって、前記撮像素子を用いて画像を撮像する場合に、当該撮像素子の露光期間における前記撮像装置の振動を検出する検出手段と、前記検出手段により検出された前記振動の時間特性を取得する取得手段と、前記振動の時間特性を、前記振動の周波数について複数の帯域に分離し、前記振動に起因して前記画像に生じる被写体像の歪み量の周波数特性を、各帯域について算出する算出手段と、前記複数の帯域の前記歪み量の周波数特性のうち、当該歪み量の周波数特性が最大となる帯域を、前記歪み量の中心周波数が含まれる帯域として決定する決定手段と、少なくとも前記中心周波数が含まれる帯域の1つの周波数について、前記取得手段により生じた位相ズレがなくなるように、前記振動の時間特性を補償する補償手段と、前記補償手段により補償された前記振動から算出された補正量を用いて、前記撮像素子のラインごとに前記振動に起因して生じた被写体像の歪みを補正する補正手段と、を備えることを特徴とする撮像装置。前記算出手段は、前記複数の帯域に分離された前記振動の時間特性を各帯域で前記露光期間について積算することにより、前記複数の帯域の歪み量の周波数特性を算出することを特徴とする請求項1に記載の撮像装置。前記補償手段は、少なくとも前記中心周波数が含まれる帯域の中央の周波数について、前記IIRフィルタにより生じた位相ズレがなくなるように、前記振動の時間特性を補償することを特徴とする請求項1または2に記載の撮像装置。前記算出手段は、前記複数の帯域のそれぞれが、少なくとも他の帯域と重複する周波数帯を含むように構成される場合は、前記画像を撮像時の前記撮像装置のパラメータに基づいて決定された前記複数の帯域それぞれについての重み付け係数を、前記複数の帯域の歪み量の周波数特性に乗じることを特徴とする請求項1乃至3のいずれか1項に記載の撮像装置。前記撮像装置のパラメータは、撮像時の振動、焦点距離、主被写体情報、シャッタスピード、及びシーンモードの少なくともいずれかであり、前記算出手段は、前記撮像時のパラメータに基づいて前記画像を撮像時の前記撮像装置の撮影状況を推定し、当該推定された撮影状況ごとに予め定められた重み付け係数を用いて、前記複数の帯域それぞれについての重み付け係数を決定することを特徴とする請求項4に記載の撮像装置。前記算出手段は、前記推定された撮影状況ごとに予め定められた重み付け係数と当該撮影状況の推定の信頼度とに基づいて、前記複数の帯域それぞれについての重み付け係数を決定することを特徴とする請求項5に記載の撮像装置。前記算出手段は、前記撮像素子の連続するライン間で前のラインが露光開始してから次のラインが露光開始するまでの期間の前記振動に起因して前記画像に生じる被写体像の歪み量の周波数特性を、各帯域について算出することを特徴とする請求項1乃至6のいずれか1項に記載の撮像装置。前記算出手段は、前記振動の時間特性のサンプリング周波数をダウンサンプルするダウンサンプリング手段を含むことを特徴とする請求項1乃至7のいずれか1項に記載の撮像装置。ローリングシャッタ方式で駆動する撮像素子を備える撮像装置の制御方法であって、検出手段が、前記撮像素子を用いて画像を撮像する場合に、当該撮像素子の露光期間における前記撮像装置の振動を検出する検出工程と、取得手段が、前記検出工程において検出された前記振動の時間特性を取得する取得工程と、算出手段が、前記振動の時間特性を、前記振動の周波数について複数の帯域に分離し、前記振動に起因して前記画像に生じる被写体像の歪み量の周波数特性を、各帯域について算出する算出工程と、決定手段が、前記複数の帯域の前記歪み量の周波数特性のうち、当該歪み量の周波数特性が最大となる帯域を、前記歪み量の中心周波数が含まれる帯域として決定する決定工程と、補償手段が、少なくとも前記中心周波数が含まれる帯域の1つの周波数について、前記取得工程において生じた位相ズレがなくなるように、前記振動の時間特性を補償する補償工程と、補正手段が、前記補償工程において補償された前記振動から算出された補正量を用いて、前記撮像素子のラインごとに前記振動に起因して生じた被写体像の歪みを補正する補正工程と、を備えることを特徴とする撮像装置の制御方法。コンピュータを、請求項1乃至8のいずれか1項に記載の撮像装置の各手段として機能させるためのプログラム。しかしながら、角速度センサのサンプリング周波数が数kHz~数十kHzであるのに対し、手ぶれの周波数は0Hz~15Hz程度であり、両者の周波数はオーダーが著しく異なる。このような角速度センサの出力信号から手ぶれの時間特性を取得するために、特許文献1に示されるFIR有限インパルス応答フィルタを用いて信号処理する場合、多くのタップ数が必要となるため、回路規模を増大させてしまうことになる。回路規模の増大を回避するためには、FIRフィルタではなくIIR無限インパルス応答フィルタを用いることが考えられるが、IIRフィルタは直線位相ではないため、周波数帯ごとに異なる位相遅れあるいは位相進みを生じる。即ち、このようにIIRフィルタを用いて得られた手ぶれの時間特性にはIIRフィルタの位相特性による位相ズレが生じているため、手ぶれの周波数によってはローリングシャッタ問題による歪みを補正するための補正量が正しく算出されないことがある。つまり、IIRフィルタを用いて得られた手ぶれの時間特性から算出した補正量を用いて、ローリングシャッタ問題による歪みを補正した場合、画像に補整残り、あるいは過補正といった現象が生じることになる。本発明は、上述の問題点に鑑みてなされたものであり、ローリングシャッタ問題による歪みを補正する場合の補正残り、あるいは過補正を低減することを目的とする。前述の目的を達成するために、本発明の撮像装置の一態様は、以下の構成を備える。ローリングシャッタ方式で駆動する撮像素子を備える撮像装置であって、撮像素子を用いて画像を撮像する場合に、当該撮像素子の露光期間における撮像装置の振動を検出する検出手段と、検出手段により検出された振動の時間特性を取得する取得手段と、振動の時間特性を、振動の周波数について複数の帯域に分離し、振動に起因して画像に生じる被写体像の歪み量の周波数特性を、各帯域について算出する算出手段と、複数の帯域の歪み量の周波数特性のうち、当該歪み量の周波数特性が最大となる帯域を、歪み量の中心周波数が含まれる帯域として決定する決定手段と、少なくとも中心周波数が含まれる帯域の1つの周波数について、取得手段により生じた位相ズレがなくなるように、振動の時間特性を補償する補償手段と、補償手段により補償された振動から算出された補正量を用いて、撮像素子のラインごとに振動に起因して生じた被写体像の歪みを補正する補正手段と、を備えることを特徴とする。このような構成により本発明によれば、ローリングシャッタ問題による歪みを補正する場合の補正残り、あるいは過補正を低減することが可能となる。以下、本発明の好適な一実施形態について、図面を参照して詳細に説明する。なお、以下に説明する一実施形態は、撮像装置の一例としての、ローリングシャッタ問題による歪みのうち、手ぶれに起因して生じる歪みを補正可能なデジタルビデオカメラに、本発明を適用した例を説明する。しかし、本発明は、ローリングシャッタ問題による歪みのうち、手ぶれに起因して生じる歪みを補正することが可能な任意の機器に適用可能である。また、本明細書において、ローリングシャッタ問題により生じる歪みのうち、手ぶれに起因して生じる歪みをローリング歪みとして記載し、以下に説明するものとする。デジタルビデオカメラ100の機能構成図1は、本発明の実施形態に係るデジタルビデオカメラ100の機能構成を示すブロック図である。制御部101は、例えばCPUであり、デジタルビデオカメラ100が備える各ブロックの動作を制御する。具体的には制御部101は、例えばROM102に記憶されている撮影処理等の処理プログラムを読み出して、不図示の揮発性メモリに展開して実行することにより、デジタルビデオカメラ100が備える各ブロックの動作を制御する。ROM102は、書き換え可能な不揮発性メモリであり、上述した処理プログラムに加え、各ブロックの動作において必要なパラメータ等が記憶される。なお、本実施形態ではハードウェアとしてデジタルビデオカメラ100が備える各ブロックにおいて各処理が実現されるものとして説明するが、本発明の実施はこれに限らず、各ブロックの処理は当該各ブロックと同様の処理を行うプログラムで実現されてもよい。撮像部103は、例えばCMOSセンサ等のローリングシャッタ方式の撮像素子を備えるブロックである。撮像部103は、不図示の光学系により撮像素子に結像された光学像を光電変換して得られたアナログ画像信号に対してA/D変換処理を適用し、デジタル画像信号画像データを出力する。撮像部103より出力された画像データは、例えば画像メモリ104に一時的に格納される。画像処理部107は、撮像部103により撮像されて画像メモリ104に格納された画像データに対して、ゲイン調整、ホワイトバランス調整等の画像処理を適用するブロックである。また画像処理部107は、後述する歪み補正量算出部106により算出された歪み補正量を用いて、画像メモリ104に格納された画像データに生じている、ローリング歪みをラインごとに補正する。デジタルビデオカメラ100を用いて動画撮影を行う場合、画像処理部107において各種画像処理あるいはローリング歪み補正処理が実行された画像データは、所定の符号化方式に従って画像処理部107において符号化される。そして符号化されて得られた動画データは、記録媒体108に伝送されて記録される。ブレ検出部105は、例えば角速度センサであって、デジタルビデオカメラ100に加わったぶれ手ぶれにより生じる角速度を検出するセンサである。ブレ検出部105は、検出したデジタルビデオカメラ100の3軸の角速度を、例えば歪み補正量算出部106に出力する。歪み補正量算出部106は、手ぶれにより画像データ内の像に生じたローリング歪みを補正する補正量を算出するブロックであり、算出した補正量の情報を画像処理部107に出力する。歪み補正量算出部106の内部構成ここで、歪み補正量算出部106の内部構成、及び歪み補正量算出部106内で行われる処理について、図を用いて以下に詳細に説明する。図2は、本実施形態に係る歪み補正量算出部106の内部構成を示した図である。歪み補正量算出部106には、ブレ検出部105により検出された、1フレームの画像の撮像中に生じた角速度変化の情報に加え、当該画像の撮像時のデジタルビデオカメラ100の状態を示す補助情報が制御部101により入力される。補助情報とは、例えば不図示のジャイロセンサの最大加速度ジャイロデータ、焦点距離、主被写体情報、シャッタスピード、あるいはシーンモード等の情報であり、本実施形態では当該補助情報は補助情報解析部209に入力される。ブレ量算出部201は、入力された1フレームの画像の撮像中に生じた角速度変化の情報から、当該角速度変化による撮像素子上の画素を単位とした像のブレ量の時間特性を時間特性する。ブレ量算出部201は、図3に示すような構成となっており、まず入力された角速度変化の情報はIIRフィルタ301で手ぶれの周波数に合わせて信号処理される。さらにブレ量算出部201により出力された角速度変化情報は、積分器302において積分処理がなされることにより、手ぶれの角変位情報θとして出力される。なお、手ぶれの角変位情報θは、例えばローリング歪みが生じうる、撮像素子の撮像面に対するピッチ方向及びヨー方向の角度について出力される。ブレ量変換部303は、積分器302から入力された角変位情報θを、当該角変異により生じる、撮像素子上の画素を単位とした像のブレ量に変換する。具体的にはブレ量変換部303は、以下の式を用いて角変位情報θから、撮像素子上の像のブレ量に変換する。fLtanθ/pここで、fLは焦点距離、pは撮像素子上の画素ピッチである。このように算出されたブレ量の時間特性は、位相補償部207及びダウンサンプル部202に出力される。なお、ブレ量算出部201では、IIRフィルタ301による手ぶれの周波数に応じた位相ズレ、及び積分器302の積分演算によるπ/2の位相遅れが生じることになる。このうち後者の積分演算による位相遅れについては手ぶれの周波数によらず固定の位相ズレであるため、π/2位相進めるように処理することにより影響をなくすことができる。本実施形態では、手ぶれの周波数に応じた位相ズレを低減するために、1フレームの撮影中に生じた手ぶれの周波数分布から手ぶれの中心周波数を特定する。そして、後述のブロックの処理において、少なくとも当該中心周波数について位相ズレが積分演算によるπ/2のみとなるように位相補償パラメータを算出し、位相補償部207に出力する。ダウンサンプル部202は、入力された手ぶれの時間特性ブレ量のサンプリング周波数を低減する。ダウンサンプル部202は、例えば手ぶれの時間特性のサンプリング周波数を12kHzから60Hzに低減する場合は、200データにつき1つのデータを取得する、所謂間引きによりサンプリング周波数を変更すればよい。フィルタ処理部203は、ダウンサンプリングされた手ぶれの時間特性に対して、複数の帯域制限フィルタを適用することにより、周波数特性を簡易的に解析する。フィルタ処理部203は図4に示すような構成であり、入力されたダウンサンプリング後の手ぶれの時間特性に対して、まず3種類の帯域フィルタを適用し、周波数帯ごとに手ぶれの時間特性を分離する。各帯域制限フィルタを通過したダウンサンプリング後の手ぶれの時間特性は、第1応答算出部404、第2応答算出部405、及び第3応答算出部406において、1フレームの露光期間について積算される。即ち、応答算出部からの出力フィルタ応答は、手ぶれの時間特性についてのブレ量の周波数応答を帯域ごとに積分した結果に相当しており、手ぶれの時間特性の3つの帯域についての周波数特性とみなすことができる。図4では、第1応答算出部404から出力されたフィルタ応答をA、第2応答算出部405から出力されたフィルタ応答をB、第3応答算出部406から出力されたフィルタ応答をCとして示している。なお、帯域制限フィルタとしてIIRフィルタを用いると、手ぶれの時間特性の周波数ごとにさらに位相ズレを生じさせてしまうため、本実施形態ではFIRフィルタによる帯域制限を行うために、ダウンサンプル部202を設けている。即ち、手ぶれの時間特性のサンプリング周波数を制限帯域に対して十分低くすることにより、FIRフィルタを用いても回路規模の増大を抑えつつ帯域制限を行うことができる。なお、各帯域制限フィルタの帯域については、手ぶれの時間特性の周波数成分を考慮して、フィルタ処理部203において決定すればよい。本発明は、IIRフィルタを用いてローリングシャッタ問題による歪みを補正する場合の補正残り、あるいは過補正を低減することが課題であり、当該課題の解決のために振動の中心周波数を動的に決定するものである。振動の中心周波数は、手ぶれの周波数が分布する帯域を細分化し、重複する帯域を有さずに多数の帯域制限フィルタを設けることにより、より正確に特定することができる。しかしながら、多数の帯域制限フィルタを設けることは回路規模の増大につながるため、本実施形態では帯域制限フィルタの数を制限することにより、回路規模を抑えて動的に振動の中心周波数を決定する方法について以下に説明する。上述のような重複する帯域を有さない帯域制限フィルタの数を制限すると、2つの帯域制限フィルタ間の帯域の閾値付近に中心周波数がある場合、例えば動画撮影時など、フレーム間で中心周波数があると判別される帯域が異なることがある。このような場合、フレーム間で異なる帯域に合わせた位相補償が行われることになり、最終的に得られる動画において閲覧者に違和感を与える可能性がある。そこで本実施形態では、手ぶれの時間特性が0~12Hzの周波数成分を主に有するものとして、帯域制限フィルタであるFIRローパスフィルタ401、第1FIRバンドパスフィルタ402、及び第2FIRバンドパスフィルタ403を以下のように設定する。FIRローパスフィルタ401:6Hz以下第1FIRバンドパスフィルタ402:3~9Hz第2FIRバンドパスフィルタ403:6~12Hz即ち、フレーム間で異なる位相補償が行われるような状況を回避するために、各帯域制限フィルタは重複する帯域を有するように設定される。なお、帯域制限フィルタが重複する帯域を有する場合は、当該重複する帯域の振動については、帯域ごとのフィルタ応答の差に現れない。このため、本実施形態では後述する重み付け部204において、撮影時の撮影状況の推定結果を用いて帯域ごとに異なる重み付けすることにより、重複する帯域の振動を帯域ごとのフィルタ応答の差として現れるようにする。また、撮影状況のパラメータを考慮した重み付けを行うことで、フレーム間で異なる帯域に振動の中心周波数が分類されることを回避することができる。重み付け部204は、例えば乗算器であり、フィルタ処理部203から入力された各帯域のフィルタ応答に対して、推定した撮影状況に応じた重み付けを乗算し、得られた重み付きフィルタ応答を中心周波数決定部205に出力する。ここで、重み付け部204で各帯域のフィルタ応答に乗じる重み付け係数について説明する。各帯域のフィルタ応答に対して乗じる重み付け係数は、歪み補正量算出部106に入力された補助情報を用いて補助情報解析部209及び重み決定部210により決定される。補助情報解析部209は、入力された補助情報から、デジタルビデオカメラ100の撮影状況が、両手持ち撮影、片手持ち撮影、歩き撮り撮影、及び不明な撮影状況のいずれであるかを推定する。一般的に両手持ち撮影、片手持ち撮影、歩き撮り撮影の撮影状況では、撮影者による把持方法や移動状態により、デジタルビデオカメラ100に生じうる手ぶれの周波数分布が図5aのように異なる。具体的には、移動しながらの撮影である歩き撮り撮影では、手ぶれによるブレ量は相対的に大きくなり、これに対して撮影者が移動せずに両手でデジタルビデオカメラ100を構える両手持ち撮影では、手ぶれによるブレ量は相対的に小さくなると考えられる。ジャイロデータについては、例えば2~3Hz程度の1つの低周波を検波することで、当該周波数におけるブレ量を、閾値を用いて大きい、標準、及び小さいのいずれかに分類する。そして補助情報解析部209は、ブレ量が大きい場合は歩き撮り撮影、小さい場合は両手持ち撮影と推定する。焦点距離については、ズーム位置がテレ寄りであるほど、手ぶれが像の歪みとして現れやすいため、撮影者はデジタルビデオカメラ100をしっかりと固定して撮影する状況が推定される。即ち、補助情報解析部209は、焦点距離が長い場合は両手持ち撮影、焦点距離が短い場合は歩き撮り撮影と推定する。主被写体情報については、ズーム位置がワイド時の撮像画像における主被写体の大きさが小さい、あるいは主被写体までの距離が遠い場合、即ちズーム位置をテレ寄りにした撮影がなされると思われる。即ち、補助情報解析部209は、主被写体の大きさが小さい、あるいは主被写体までの距離が遠い場合は両手持ち撮影、主被写体の大きさが大きい、あるいは主被写体までの距離が近い場合は歩き撮り撮影と推定する。シャッタスピードについては、露光時間が長い場合は蓄積ブレの影響が大きくなり、低周波のブレに起因する撮像画像の歪みは目立ちにくくなると考えられる。即ち、補助情報解析部209は、シャッタスピードが遅い場合は低周波成分のブレ量が少ない両手持ち撮影と推定し、位相補償周波数を高周波側に寄せる。シーンモードについては、ポートレートモードに設定された場合は、撮影者はデジタルビデオカメラ100を両手でしっかりと固定して撮影すると考えられるため、補助情報解析部209は両手持ち撮影と推定する。またスポーツモードに設定された場合は、動いている被写体を追いやすいように、撮影者はデジタルビデオカメラ100を片手持ちで撮影すると考えられるため、補助情報解析部209は片手持ち撮影と推定する。さらに、風景モードに設定された場合は、撮影者が野山等を散策しながら撮影することが想定されるため、補助情報解析部209は歩き撮り撮影と推定する。また、撮影状況が限定できないようなシーンモードについては、補助情報解析部209は撮影状況を不明と推定する。即ち、補助情報解析部209は、入力された補助情報のそれぞれについて、予め定められた状態や閾値を用いて、それぞれのパラメータを図7のように分類し、補助情報のパラメータごとに撮影状況を推定する。このようにすることで、一般的なブレ量の周波数分布から中心周波数が存在すると思われる周波数帯に大きい重み付けがなされるように、重み決定部210で重み付け係数を決定することができる。重み決定部210は、補助情報解析部209で推定された補助情報のパラメータごとの撮影状況から、フィルタ処理部203の出力である各帯域のフィルタ応答に乗じる重み付け係数を決定する。撮像画像における歪み量の周波数分布を考えると、中心周波数が撮影状況ごとに例えば図5bのような傾向になる。このため、当該傾向に基づいて、実際の手ぶれの時間特性の中心周波数が存在する帯域に大きい重み付けがなされるようにすることで、中心周波数が帯域制限フィルタの重複する帯域にある場合であっても、適切な帯域を中心周波数のある帯域とすることができる。重み決定部210は、まず補助情報のそれぞれのパラメータについての各帯域の重み付け係数を、当該それぞれのパラメータについて推定された撮影状況に応じて選択する。例えば、推定された撮影状況が歩き撮り撮影である場合は図5bのように中心周波数は低周波領域にあると考えられるため、図6に示すように当該帯域について大きい重み付けとなり、周波数が大きい帯域ほど小さい重み付けとなるように係数が設定される。そして重み決定部210は、補助情報のそれぞれのパラメータに対して選択された各帯域の重み付け係数を組み合わせて、重み付け部204で各帯域のフィルタ応答に乗じる最終的な各帯域の重み付け係数を決定する。例えば、補助情報の2つのパラメータを用いて最終的な各帯域の重み付け係数を決定する場合を考える。1つのパラメータについて選択された各帯域の重み付け係数が図8aであり、もう一方のパラメータについて選択された各帯域の重み付け係数が図8bである場合、最終的な各帯域の重み付け係数は帯域ごとに平均をとり、図8cのようになる。重み付け部204は、重み決定部210が決定した最終的な各帯域の重み付け係数を、対応するフィルタ応答に乗じることにより、重み付きフィルタ応答を出力する。例えば重み付け部204に入力された各帯域のフィルタ応答が図9aのようである場合、重み決定部210により図9bのような重み付け係数が決定されると、各帯域の重み付きフィルタ応答は図9cのようになる。中心周波数決定部205は、入力された各帯域の重み付きフィルタ応答から、1フレームの撮像中に生じた手ぶれの中心周波数を決定する。手ぶれの時間特性のうち、単一周波数fの振動は、当該周波数における最大ブレ量である振幅Afを用いて、Afsinωtω=2πfと表すことができる。即ち、周波数ごとの振動の速度は、手ぶれの時間特性を微分してAfωcosωtと表すことができる。例えば、撮像素子の連続するライン間の露光開始の時刻の差がΔtである場合、前のラインの撮像が開始されてから次のラインの撮像が開始されるまでの期間Δtでの像の変位量は、図10aのようなΔtの中央の時刻t0の速度を平均速度とみなし、Afωcosωt0Δtと表すことができる。上述したようにフィルタ処理部203より得られた各帯域のフィルタ応答図10bは、手ぶれの時間特性についてのブレ量の周波数応答を帯域ごとに積分した結果に相当している。このため、当該式からブレ量の項Afを分離し、ωt0が一定として近似すると、残りの項はωcosωt0Δt⇒2πcosωt0Δt×f=c×fc=const.と変形でき、図10cのような周波数に対する一次関数で表すことができる。即ち、中心周波数決定部205では、入力された図10bのような各帯域の重み付きフィルタ応答に対し、図10cから各帯域の平均周波数に応じた変位量を乗じる。これにより、図10dのような、撮像画像における歪み量を考慮した各帯域のフィルタ応答を得ることができる。本実施形態では、撮像画像における歪み量を考慮したフィルタ応答が最大となる帯域の平均周波数fa周波数帯域の中央に相当する周波数3~9Hzの場合は6Hzを中心周波数として決定する。そして、当該中心周波数においてIIRフィルタ301による位相遅れがπ/2となるように後述する位相補償制御部206に位相補償パラメータを設定させる。このようにすることで、歪み量が最も大きい、即ち位相ズレが生じた場合に像の補整残りまたは過補正による歪みが目立ちやすい周波数中心周波数において歪みを軽減することができる。なお、本実施形態では、撮像画像における歪み量を考慮したフィルタ応答が極大となる帯域を中心周波数が含まれる帯域として決定し、当該帯域の中央に相当する周波数を中心周波数として決定するものとして説明するが、本発明の実施はこれに限らない。即ち、中心周波数として決定する周波数は、中心周波数が含まれる帯域のうちの1つの周波数、あるいは周波数帯であってよく、他の帯域のフィルタの傾向により動的に選択されてもよい。例えば、図10dの場合は、低周波領域と高周波領域のフィルタ応答の比に従って、中間周波領域のうちの1つの周波数を中心周波数として決定してもよい。位相補償制御部206は、位相補償部207での手ぶれの時間特性の位相補償処理において、中心周波数決定部205で決定した中心周波数近傍の位相遅れがπ/2となるような位相補償パラメータを算出して、位相補償部207に出力する。なお、位相補償パラメータは、中心周波数近傍の周波数帯において、位相補償後のブレ量の位相が目標位相である-π/2となるように、公知の方法で決定される。位相補償部207は、位相補償フィルタであり、例えば当該位相補償フィルタの位相特性は図12aのようになる。位相補償部207は、入力された手ぶれの時間特性を、中心周波数決定部205で決定された中心周波数でIIRフィルタ301により生じた位相ズレを補償し、当該周波数での位相遅れが積分器302により生じたπ/2となるように位相補償を行う。位相補償部207の内部構成は、例えば図11のようになる。加算器1101は、入力された手ぶれの時間特性ブレ量と、後述する増幅器1103の出力とを足し合わせ、位相補償後のブレ量を撮像画像歪み量算出部208に出力する。またシフトレジスタ1102は、入力された位相補償後のブレ量を遅延させ、増幅器1103に出力する。増幅器1103は、シフトレジスタ1102からの出力に対して、位相補償制御部206より入力された位相補償パラメータであるフィルタ係数Kを乗算して加算器1101に出力する。このような構成により、位相補償制御部206で適切に設定された位相補償パラメータを用いて、位相補償部207の出力であるブレ量の位相特性を、図12bのように中心周波数において-π/2とすることができる。撮像画像歪み量算出部208は、入力された位相補償後の手ぶれの時間特性を用いて、ブレ量算出部201における位相ズレが補償された、撮像素子のラインごとの撮像画像の歪み量を算出する。撮像素子のライン間の撮像画像の歪み量は、例えばNライン目とN1ライン目の間に生じる歪み量を例にとると、図13に示すように露光開始時刻におけるブレ量の差から算出することができる。このように算出された歪み量をなくすように決定された補正量を用いて、画像処理部107は、撮像素子のラインごとに、ローリングシャッタ問題による撮像画像の歪みを補正することができる。まとめ以上説明したように、本実施形態の撮像装置は、ローリングシャッタ問題による歪みを補正する場合の補正残り、あるいは過補正を低減することができる。具体的には撮像装置は、ローリングシャッタ方式で駆動する撮像素子を用いて画像を撮像する場合に、当該撮像素子の露光期間における撮像装置の振動を検出し、検出された振動の時間特性を取得する。振動の時間特性を、振動の周波数について複数の帯域に分離し、振動に起因して画像に生じる被写体像の歪み量の周波数特性を、各帯域について算出し、当該歪み量の周波数特性が極大となる帯域を、歪み量の中心周波数が含まれる帯域として決定する。そして、少なくとも中心周波数が含まれる帯域の1つの周波数について、振動の時間特性を取得する際に生じた位相ズレがなくなるように振動の時間特性を補償する。さらに補償後の振動の時間特性から算出された補正量を用いて、撮像素子のラインごとに振動に起因して生じた被写体像の歪みを補正する。このようにすることで、撮像画像における歪み量が最も大きく、補整残りあるいは過補正による歪みが目立ちやすい中心周波数近傍の周波数について、IIRフィルタを用いたことによる位相ズレを補償することができる。このため、ローリングシャッタ問題による歪みを補正した場合に、撮像画像全体において歪みを目立たなくすることができる。さらに、1フレームの画像を撮像中の手ぶれの時間特性に応じて動的に中心周波数を決定するため、適切にローリングシャッタ問題による歪みを補正することができる。変形例上述した実施形態では、重み決定部210は補助情報のパラメータそれぞれについて選択された各帯域の重み付け係数の平均をとって、各帯域のフィルタ応答に乗じる最終的な重み付け係数を決定するものとして説明した。本変形例では、補助情報のパラメータそれぞれについて推定した撮影状況について信頼度を設け、各撮影状況の信頼度に基づいて各帯域の重み付け係数を決定する方法について説明する。例えば、補助情報のパラメータそれぞれについて推定した撮影状況について、図14に示すような推定の信頼度が設定されている場合を考える。重み決定部210は、まず上述した実施形態と同様に、補助情報解析部209により補助情報のパラメータそれぞれについて推定された撮影状況に対して、予め定められた各帯域の重み付け係数を選択する。そして重み決定部210は、選択した各帯域の重み付け係数について、当該選択に用いた補助情報のパラメータの信頼度に応じて、調整を行う。例えば補助情報の1つのパラメータについて推定された撮影状況から、図15a各帯域の重み付け係数が選択された場合を考える。重み決定部210は、選択された各帯域の重み付け係数について、重み付け係数が高い順に選択された2つの帯域間の差分について信頼度を乗じて得られた値が差分となるように、係数が最も高い低周波領域の係数を基準として各帯域の重み付け係数を調整する。具体的には撮影状況の信頼度が0.8であるとすると、低周波領域と中間周波領域の係数の差分0.3、及び中間周波数領域と高周波領域の係数の差分0.2は、それぞれ0.24、0.16となる。即ち、図15bのように各帯域の係数は1.0、0.76、0.6となる。そして重み決定部210は、補助情報のパラメータそれぞれについて撮影状況の信頼度に基づいて調整された各帯域の重み付け係数を平均することにより、各帯域のフィルタ応答に乗ずる、最終的な各帯域の重み付け係数を決定する。このように、補助情報である撮影時のデジタルビデオカメラ100の状態に応じて、撮影状況を推定する際に、推定の信頼度を用いることで、より撮影時の状態に柔軟に対応した各帯域の重み付け係数を用いて、手ぶれの時間特性の中心周波数を決定できる。その他の実施形態また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェアプログラムを、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータまたはCPUやMPU等がプログラムを読み出して実行する処理である。