"短時間露光画像と長時間露光画像の少なくとも2種類の異なる露光時間の撮影画像を入力して、短時間露光画像に含まれるフリッカ成分を除去するフリッカ補正処理を実行する信号処理部を有し、前記信号処理部は、前記短時間露光画像の行単位の信号量の積分値からなる短時間露光画像行プロファイルと、前記長時間露光画像の行単位の信号量の積分値からなる基準プロファイルを生成し、前記短時間露光画像行プロファイルに対する前記基準プロファイルによる除算処理を含むデータ処理によってフリッカ補正波形を生成し、前記フリッカ補正波形に対応する行単位の係数を、短時間露光画像の各行の画素値に乗算する処理によりフリッカ成分を除去した補正画像を生成する画像処理装置。前記長時間露光画像は、画像撮影環境の輝度変化周期の整数倍にほぼ一致する露光時間の設定の下で撮影された画像である請求項1に記載の画像処理装置。前記信号処理部は、前記フリッカ補正波形の生成処理に際して、前記短時間露光画像行プロファイルに対する前記基準プロファイルによる除算処理によって得られるフリッカ成分データ、または、該フリッカ成分データの位相反転データに対する正弦波近似処理を実行して前記フリッカ補正波形を生成する請求項1に記載の画像処理装置。前記信号処理部は、前記フリッカ補正波形の生成処理に際して、前記短時間露光画像行プロファイルに対する前記基準プロファイルによる除算処理によって得られるフリッカ成分データ、または、該フリッカ成分データの位相反転データに対するスプライン補間処理を実行して前記フリッカ補正波形を生成する請求項1に記載の画像処理装置。前記信号処理部は、前記フリッカ補正波形の生成処理に際して、前記短時間露光画像行プロファイルに対する前記基準プロファイルによる除算処理によって得られるフリッカ成分データ、または、該フリッカ成分データの位相反転データに対する線型補間処理を実行して前記フリッカ補正波形を生成する請求項1に記載の画像処理装置。前記信号処理部は、前記フリッカ補正波形の生成処理に際して、前記短時間露光画像行プロファイルに対する前記基準プロファイルによる除算処理によって得られるフリッカ成分データ、または、該フリッカ成分データの位相反転データに対する離散フーリエ変換を実行して、2次以降の高次成分を低次に振り分ける処理を行う請求項1に記載の画像処理装置。前記信号処理部は、フリッカ成分を除去した補正画像の生成処理に際して、補正対象画像の行プロファイルのピーク位置と、前記フリッカ補正波形のピーク位置を一致させる位相合わせ後に、前記フリッカ補正波形に対応する行単位の係数を、短時間露光画像の各行の画素値に乗算する処理を実行する請求項1に記載の画像処理装置。前記信号処理部は、各画像に動被写体を含むか否かを判定し、動被写体を含まないと判定した画像を適用して前記フリッカ補正波形の生成処理を行う請求項1に記載の画像処理装置。前記信号処理部は、補正対象とする短時間露光画像に動被写体を含むか否かを判定し、動被写体を含まないと判定した画像に対して、前記フリッカ補正波形を適用した画像補正を行う請求項1に記載の画像処理装置。前記信号処理部は、補正対象とする短時間露光画像に動被写体を含むか否かを判定し、動被写体を含むと判定した画像については、動被写体の影響による行プロファイルの位相変化を考慮して、前記フリッカ補正波形を修正し、修正したフリッカ補正波形を適用した画像補正を行う請求項1に記載の画像処理装置。前記信号処理部は、補正対象とする短時間露光画像に動被写体を含むか否かを判定し、動被写体を含むと判定した画像については、静止画像に基づいて算出したフレーム間位相ずれに基づいて、動被写体を含む画像と位相が一致するように前記フリッカ補正波形の切り出しを実行して、切り出したフリッカ補正波形を適用した画像補正を行う請求項1に記載の画像処理装置。画像処理装置において実行する画像処理方法であり、信号処理部が、短時間露光画像と長時間露光画像の少なくとも2種類の異なる露光時間の撮影画像を入力して、短時間露光画像に含まれるフリッカ成分を除去するフリッカ補正処理を実行する信号処理ステップを実行し、前記信号処理ステップにおいては、前記短時間露光画像の行単位の信号量の積分値からなる短時間露光画像行プロファイルと、前記長時間露光画像の行単位の信号量の積分値からなる基準プロファイルの生成処理と、前記短時間露光画像行プロファイルに対する前記基準プロファイルによる除算処理を含むデータ処理によるフリッカ補正波形の生成処理と、前記フリッカ補正波形に対応する行単位の係数を、短時間露光画像の各行の画素値に乗算する処理によりフリッカ成分を除去した補正画像の生成処理を実行する画像処理方法。画像処理装置において画像に含まれるフリッカ成分を除去するフリッカ補正処理を実行させるプログラムであり、信号処理部に、短時間露光画像と長時間露光画像の少なくとも2種類の異なる露光時間の撮影画像を入力し、前記信号処理部に、前記短時間露光画像の行単位の信号量の積分値からなる短時間露光画像行プロファイルと、前記長時間露光画像の行単位の信号量の積分値からなる基準プロファイルの生成処理と、前記短時間露光画像行プロファイルに対する前記基準プロファイルによる除算処理を含むデータ処理によるフリッカ補正波形の生成処理と、前記フリッカ補正波形に対応する行単位の係数を、短時間露光画像の各行の画素値に乗算する処理によりフリッカ成分を除去した補正画像の生成処理を実行させるプログラム。本開示は、例えば上記の問題に鑑みてなされたものであり、より高精度なフリッカ補正を行う画像処理装置、および画像処理方法、並びにプログラムを提供することを目的とする。また、本開示の一実施例の構成においては、例えば過去画像の積算処理等を行うことなく、処理負荷を軽減して高精度なフリッカ補正を行う画像処理装置、および画像処理方法、並びにプログラムを提供する。本開示の第1の側面は、短時間露光画像と長時間露光画像の少なくとも2種類の異なる露光時間の撮影画像を入力して、短時間露光画像に含まれるフリッカ成分を除去するフリッカ補正処理を実行する信号処理部を有し、前記信号処理部は、前記短時間露光画像の行単位の信号量の積分値からなる短時間露光画像行プロファイルと、前記長時間露光画像の行単位の信号量の積分値からなる基準プロファイルを生成し、前記短時間露光画像行プロファイルに対する前記基準プロファイルによる除算処理を含むデータ処理によってフリッカ補正波形を生成し、前記フリッカ補正波形に対応する行単位の係数を、短時間露光画像の各行の画素値に乗算する処理によりフリッカ成分を除去した補正画像を生成する画像処理装置にある。さらに、本開示の画像処理装置の一実施態様において、前記長時間露光画像は、画像撮影環境の輝度変化周期の整数倍にほぼ一致する露光時間の設定の下で撮影された画像である。さらに、本開示の画像処理装置の一実施態様において、前記信号処理部は、前記フリッカ補正波形の生成処理に際して、前記短時間露光画像行プロファイルに対する前記基準プロファイルによる除算処理によって得られるフリッカ成分データ、または、該フリッカ成分データの位相反転データに対する正弦波近似処理を実行して前記フリッカ補正波形を生成する。さらに、本開示の画像処理装置の一実施態様において、前記信号処理部は、前記フリッカ補正波形の生成処理に際して、前記短時間露光画像行プロファイルに対する前記基準プロファイルによる除算処理によって得られるフリッカ成分データ、または、該フリッカ成分データの位相反転データに対するスプライン補間処理を実行して前記フリッカ補正波形を生成する。さらに、本開示の画像処理装置の一実施態様において、前記信号処理部は、前記フリッカ補正波形の生成処理に際して、前記短時間露光画像行プロファイルに対する前記基準プロファイルによる除算処理によって得られるフリッカ成分データ、または、該フリッカ成分データの位相反転データに対する線型補間処理を実行して前記フリッカ補正波形を生成する。さらに、本開示の画像処理装置の一実施態様において、前記信号処理部は、前記フリッカ補正波形の生成処理に際して、前記短時間露光画像行プロファイルに対する前記基準プロファイルによる除算処理によって得られるフリッカ成分データ、または、該フリッカ成分データの位相反転データに対する離散フーリエ変換を実行して、2次以降の高次成分を低次に振り分ける処理を行う。さらに、本開示の画像処理装置の一実施態様において、前記信号処理部は、フリッカ成分を除去した補正画像の生成処理に際して、補正対象画像の行プロファイルのピーク位置と、前記フリッカ補正波形のピーク位置を一致させる位相合わせ後に、前記フリッカ補正波形に対応する行単位の係数を、短時間露光画像の各行の画素値に乗算する処理を実行する。さらに、本開示の画像処理装置の一実施態様において、前記信号処理部は、各画像に動被写体を含むか否かを判定し、動被写体を含まないと判定した画像を適用して前記フリッカ補正波形の生成処理を行う。さらに、本開示の画像処理装置の一実施態様において、前記信号処理部は、補正対象とする短時間露光画像に動被写体を含むか否かを判定し、動被写体を含まないと判定した画像に対して、前記フリッカ補正波形を適用した画像補正を行う。さらに、本開示の画像処理装置の一実施態様において、前記信号処理部は、補正対象とする短時間露光画像に動被写体を含むか否かを判定し、動被写体を含むと判定した画像については、動被写体の影響による行プロファイルの位相変化を考慮して、前記フリッカ補正波形を修正し、修正したフリッカ補正波形を適用した画像補正を行う。さらに、本開示の画像処理装置の一実施態様において、前記信号処理部は、補正対象とする短時間露光画像に動被写体を含むか否かを判定し、動被写体を含むと判定した画像については、静止画像に基づいて算出したフレーム間位相ずれに基づいて、動被写体を含む画像と位相が一致するように前記フリッカ補正波形の切り出しを実行して、切り出したフリッカ補正波形を適用した画像補正を行う。さらに、本開示の第2の側面は、画像処理装置において実行する画像処理方法であり、信号処理部が、短時間露光画像と長時間露光画像の少なくとも2種類の異なる露光時間の撮影画像を入力して、短時間露光画像に含まれるフリッカ成分を除去するフリッカ補正処理を実行する信号処理ステップを実行し、前記信号処理ステップにおいては、前記短時間露光画像の行単位の信号量の積分値からなる短時間露光画像行プロファイルと、前記長時間露光画像の行単位の信号量の積分値からなる基準プロファイルの生成処理と、前記短時間露光画像行プロファイルに対する前記基準プロファイルによる除算処理を含むデータ処理によるフリッカ補正波形の生成処理と、前記フリッカ補正波形に対応する行単位の係数を、短時間露光画像の各行の画素値に乗算する処理によりフリッカ成分を除去した補正画像の生成処理を実行する画像処理方法にある。さらに、本開示の第3の側面は、画像処理装置において画像に含まれるフリッカ成分を除去するフリッカ補正処理を実行させるプログラムであり、信号処理部に、短時間露光画像と長時間露光画像の少なくとも2種類の異なる露光時間の撮影画像を入力し、前記信号処理部に、前記短時間露光画像の行単位の信号量の積分値からなる短時間露光画像行プロファイルと、前記長時間露光画像の行単位の信号量の積分値からなる基準プロファイルの生成処理と、前記短時間露光画像行プロファイルに対する前記基準プロファイルによる除算処理を含むデータ処理によるフリッカ補正波形の生成処理と、前記フリッカ補正波形に対応する行単位の係数を、短時間露光画像の各行の画素値に乗算する処理によりフリッカ成分を除去した補正画像の生成処理を実行させるプログラムにある。なお、本開示のプログラムは、例えば、様々なプログラムコードを実行可能な情報処理装置やコンピュータシステムに対して例えば記憶媒体によって提供されるプログラムである。このようなプログラムを情報処理装置やコンピュータシステム上のプログラム実行部で実行することでプログラムに応じた処理が実現される。本開示のさらに他の目的、特徴や利点は、後述する本開示の実施例や添付する図面に基づくより詳細な説明によって明らかになるであろう。なお、本明細書においてシステムとは、複数の装置の論理的集合構成であり、各構成の装置が同一筐体内にあるものには限らない。本開示の一実施例の構成によれば、効果的なフリッカの除去を可能とする装置、方法が実現される。具体的には、輝度変化のある蛍光灯のような照明環境下における撮影画像に発生するフリッカを効果的に除去または削減する。短時間露光画像と長時間露光画像の少なくとも2種類の異なる露光時間の撮影画像を入力し、短時間露光画像の行単位の信号量の積分値からなる短時間露光画像行プロファイルと、長時間露光画像の行単位の信号量の積分値からなる基準プロファイルを生成し、短時間露光画像行プロファイルに対する基準プロファイルによる除算処理等によってフリッカ補正波形を生成し、フリッカ補正波形で規定される行単位の係数を、短時間露光画像の各行の画素値に乗算する処理によりフリッカ成分を除去した補正画像を生成する。これらの処理によって、効果的なフリッカの除去を可能とする装置、方法が実現される。以下、図面を参照しながら本発明の画像処理装置、および画像処理方法、並びにプログラムの詳細について説明する。説明は、以下の項目に従って行う。1.本開示の画像処理装置の全体構成と処理について2.フリッカ補正処理の詳細について2-1.行プロファイルRowProfile生成処理について2-2.フリッカ補正係数算出処理について2-3.動被写体検出処理について2-4.静止画と判定された画像に対するフリッカ補正処理について2-5.動被写体を含む画像に対するフリッカ補正処理について3.本開示の構成のまとめ[1.本開示の画像処理装置の全体構成と処理について]まず、本開示の画像処理装置の一実施例である撮像装置の構成例について、図2を参照して説明する。図2は、本開示の一実施例にかかる撮像装置100の全体構成を示すブロック図である。撮像装置100は、撮像部101、信号処理部102、コーデック符号化復号処理部103、記録再生部104、記録メディア120、表示部130、制御部150、メモリ151を有する。撮像部101は、レンズや撮像素子等によって構成される。撮像素子は、XYアドレス走査型の撮像素子であるCMOSイメージセンサである。XYアドレス走査型の撮像素子であるCMOSイメージセンサは、先に図1を参照して説明したように、時間経過に従って順次、上段行から下段行に向かって露光処理が行われる。CMOSイメージセンサによる受光信号は、光電変換され電気信号として信号処理部102に入力される。信号処理部102は、一般的なカメラ信号処理として実行されるホワイトバランス調整、ガンマ補正、あるいはデモザイク処理等の処理の他、本発明に従ったフリッカ補正処理を実行する。フリッカ補正処理を実行するフリッカ補正処理部の構成、および処理の詳細については、図3以下を参照して後段で説明する。信号処理部102における信号処理の結果として得られたデータはコーデック符号化復号処理部103に入力される。また、表示部130に出力され表示処理が実行される。コーデック103では例えばMPEG符号化処理等、予め規定された符号化処理を実行し、処理結果を記録再生部104に出力する。記録再生部104は、予め規定された記録フォーマットに従って、記録メディア120に対する撮影画像の記録処理を実行する。記録メディア120はたとえばフラッシュメモリ、磁気ディスク、光ディスクなどのデータ記録可能なメディアである。なお、記録メディア120に記録された画像の再生を行う場合は、記録再生部部104が記録メディア120からのデータ読み出しを実行し、読み出されたデータがコーデック103に提供され、復号処理が実行される。復号処理の結果として得られる画像が表示部130に表示される。なお、これらの各構成部の制御は、制御部150がメモリ151に予め記録されたプログラムに従って実行する。制御部150はプログラム実行機能を持つCPUを有する。メモリ151はRAM,ROM等によって構成される。図2に示す撮像装置100は、短時間露光画像、長時間露光画像、これら少なくとも2種類の異なる露光時間の画像を連続的に撮影する。さらに1つの露光時間の設定を蛍光灯の輝度周期、あるいは輝度周期の整数倍に設定する。なお、このような複数の異なる露光時間の撮影処理は、広ダイナミックレンジの画像撮影を行う処理として行われる手法である。広ダイナミックレンジ撮像においては、イメージセンサの露光時間を周期的に変えて撮像を行い、各画素の輝度レベル等に応じた画像合成処理によって広ダイナミックレンジ画像を得る。すなわち低輝度画素部分は、長時間露光画像の画素値を適用し、高輝度画素部分は短時間露光画像の画素値を用いるといった画像合成処理によって、広ダイナミックレンジの画像を生成するものである。なお、各画素値の設定の際には、全体的な整合性を保持させるための画素値算出処理が行われる。このような処理により、一定時間の露光時間の撮影画像では得られないダイナミックレンジの広い画像を得ることができる。このような広ダイナミックレンジ画像生成処理に際して、異なる露光時間の画像を合成する前の処理として、以下に説明するフリッカ補正を施すことで、フリッカを抑制した広ダイナミックレンジ画像を生成することができる。異なる露光時間の撮影シーケンスと、蛍光灯の波形との対応関係の例について図3を参照して説明する。図3は、電源周波数=60Hzの電源によって点灯した蛍光灯の下で、撮像装置が1秒間に120枚、すなわち120fpsのフレームレートで画像撮影を行った場合のシーケンスを示している。なお、電源周波数=60Hzの電源によって点灯した蛍光灯の輝度周期は120Hzであり、1秒間に120回の明暗パターンを繰り返す。撮像装置は、露光時間t1の短時間露光画像、露光時間t2の長時間露光画像、少なくともこれら2種類以上の異なる露光時間の画像を繰り返し撮影する。本開示の撮像装置は、少なくとも1つの画像は、画像撮影環境の輝度変化周期の整数倍にほぼ一致する露光時間の設定の下で撮影された画像とする。図3に示す設定は、1秒間に60枚の短時間露光画像と、60枚の長時間露光画像を交互に撮影する設定である。ここで、長時間露光画像の露光時間t2は、蛍光灯の輝度周期[1/120sec]にほぼ一致している。すなわち、t2≒1/120secこの設定である。長時間露光画像の露光時間t2は、1/120secに完全に一致させることは困難であるが、長時間露光画像の露光時間t2を、蛍光灯の輝度周期[1/120sec]にほぼ一致させることで、長時間露光画像は、フリッカの影響をほとんど含まない画像として撮影される。一方、短時間露光画像の露光時間t1は、蛍光灯の輝度周期、あるいはその倍数とは異なる設定であり、短時間露光画像は、フリッカの影響を受けた画像として撮影される。図2に示す撮像装置100の信号処理部102は、これらの複数の露光時間の設定の下で撮影された画像を適用したフリッカ補正を行う。以下、信号処理部102において実行するフリッカ補正の詳細について説明する。[2.フリッカ補正処理の詳細について]信号処理部102において実行するフリッカ補正の詳細について説明する。図4は、信号処理部102において実行するフリッカ補正処理のシーケンスを示す図である。なお、図4に示すシーケンスは、例えば、図4に示す各ステップS101~の処理を規定したプログラムに従って、信号処理部102内のプロセッサによって実行される。あるいは各処理を実行するための専用のハードウェアを利用して実行してもよい。信号処理部102は、撮像部101を構成するイメージセンサで取得されるRAW画像を入力画像201として入力し、RAW画像に対するフリッカ補正処理を行う。なお、RAW画像は撮像部101の出力データであり、デモザイク処理等の信号処理が施される前の生データである。デモザイク処理は、撮像部101を構成するイメージセンサの各画素対応のデータとして例えばRGBのすべての画素値を設定する処理であり、デモザイク処理を実行する前のRAWデータは、イメージセンサの各画素に対して例えばRGBいずれかの画素値のみが設定されたデータである。図4に示すように、信号処理部102は、デモザイク処理前のRAW画像を入力画像201として、入力画像201に対するフリッカ補正が施された出力画像202を生成する。ただし、この処理は一例であり、信号処理部102は、デモザイク処理が実行されたカラー画像を入力してカラー画像に対して、以下に説明すると同様の処理を行うことでフリッカ補正処理を実行する構成としてもよい。以下、図4に示すフリッカ補正処理シーケンスについて、各処理ステップに従って順次説明する。[2-1.行プロファイルRowProfile生成処理について]まず、ステップS101において、入力画像201に対する行プロファイルRowProfileを生成する。行プロファイルRowProfileとは、撮影画像の各行単位の信号量強度を示すデータである。RAW画像201の水平方向にカラーチャンネルRGBごとに信号値画素値の積分処理を行い、行単位の1次元の信号量Intensityを得る。前述したように、RAWデータは、イメージセンサの各画素に対して例えばRGBいずれかの画素値のみが設定されたデータである。例えば図5に示すようなデータを行プロファイルRowProfileとして生成する。行プロファイルRowProfileは、フリッカ補正としての画素値補正を行単位で実行するための基礎データとして利用する。行プロファイルの生成と行単位のフリッカ補正は、フリッカの撮影画像に対する影響が画像のスキャンライン水平方向でほぼ同一という仮定のもとに行われる処理であり、このような行単位の処理とすることで、フリッカ補正の計算量削減や処理の効率化が実現される。例えばベイヤBayer配列型のイメージセンサであれば、行単位でR,G,Bチャンネルそれぞれ独立に積分処理を行い、RGB各チャンネルの行単位の1次元の信号量RPとして、信号量:RP_Rky,RP_Gky,RP_Bkyを得る。なお、kは、画像のフレーム番号であり、yは画像の縦方向の座標値を表す。信号処理部102は、このように、入力RAW画像フレームkに対して、各行y単位で、信号量:RP_Rky,RP_Gky,RP_Bkyを算出する。なお、処理手法としては、1色信号成分、例えばRGBのカラーチャンネル単位で個別に処理を実行する手法と、2カラーチャンネル例えばRGBを区別せず、行単位で1つの信号量RPkyを算出してこれを用いて全てのカラーチャンネル対応画素に対して共通の処理を実行する手法、これら2つの手法がある。いずれの処理を行っても効果的なフリッカの削減効果を発揮することが可能である。なお、1のRGBのカラーチャンネル単位の処理は、色ムラの削減効果を高めることができる。行プロファイルRowProfile生成処理は、上記1,2のいずれの処理として実行してもよい。先に、図3を参照して説明したように、本開示の撮像装置は、露光時間t1の短時間露光画像、露光時間t2の長時間露光画像、少なくともこれら2種類以上の異なる露光時間の画像を繰り返し撮影する。図3に示す設定は、1秒間に60枚の短時間露光画像と、60枚の長時間露光画像を交互に撮影する設定である。ここで、長時間露光画像の露光時間t2は、蛍光灯の輝度周期にほぼ一致している。すなわち、t2≒1/120secこの設定である。従って、長時間露光画像は、フリッカの影響をほとんど含まない画像として撮影される。一方、短時間露光画像の露光時間t1は、蛍光灯の輝度周期、あるいはその倍数とは異なる設定であり、短時間露光画像は、フリッカの影響を受けた画像として撮影される。この結果、短時間露光画像の行プロファイルは、フリッカの影響を含む行プロファイルとなる。一方、長時間露光画像の行プロファイルは、フリッカの影響をほとんど含まない行プロファイルとなる。図6を参照して短時間露光画像の行プロファイルと長時間露光画像の行プロファイルの例について説明する。図6には、a短時間露光画像の輝度分布特徴データと、短時間露光画像の行プロファイルb長時間露光画像の輝度分布特徴データと、長短時間露光画像の行プロファイルこれらを示している。図6aに示す短時間露光画像の行プロファイルは、フリッカの影響を大きく受けた行プロファイルとなる。RGBの各信号量が行の中央部分で落ち込み、両端が高くなっているのは、フリッカの影響によるものと推定される。一方、図6bに示す長時間露光画像の行プロファイルは、フリッカの影響をほとんど含まない行プロファイルとなる。RGBの各信号量は、各行に渡ってほぼ均一であり、顕著な変化は見られない。フリッカの影響による信号量の変化はほとんど発生していないと推定される。図4に示すシーケンス図のステップS101では、短時間露光画像と長時間露光画像の各撮影画像について、行プロファイルを順次、生成する。例えば先に図3を参照して説明した撮影シーケンス、すなわち、露光時間t1の短時間露光画像、露光時間t2の長時間露光画像、これらの画像を交互に撮影するシーケンスでは、短時間露光画像の行プロファイルの生成、長時間露光画像の行プロファイルの生成、短時間露光画像の行プロファイルの生成、長時間露光画像の行プロファイルの生成、これらの処理を繰り返し実行する。[2-2.フリッカ補正係数算出処理について]次に、ステップS102において、フリッカ補正係数を算出する。なお、本実施例において、フリッカの影響が大きく発生するのは短時間露光画像であり、短時間露光画像各々に対応したフリッカ補正係数を算出する。フリッカ補正係数の算出処理について、図7、図8を参照して説明する。図7には、aフリッカ補正係数算出対象となる短時間露光画像の行プロファイルbフリッカ補正係数算出に適用する基準プロファイルとしての長時間露光画像の行プロファイルca/bの除算処理によって得られるa短時間露光画像に対応するフリッカ成分データこれらを示している。なお、図7cの縦軸に示すフリッカゲインは、フリッカゲイン=1/補正係数である。さらに、図8には、c図7cと同様のa/bの除算処理によって得られるa短時間露光画像に対応するフリッカ成分データdcに示すフリッカ成分の位相の反転処理によって生成するフリッカ補正係数≒フリッカ補正波形これらを示している。図7aに示すフリッカ補正係数算出対象の短時間露光画像の行プロファイルを、その短時間露光画像の直前または直後に撮影された長時間露光画像の行プロファイル基準プロファイルで除算各色、各対応行の信号量同士で除算する。すなわち、図7bに示す長時間露光画像の行プロファイル基準プロファイルで除算する。この除算処理によって、図7cに示すフリッカ成分データが得られる。なお、図7に示す、a/bの除算処理、この除算処理の実行前に、長時間露光画像と短時間露光画像の露光比を調整する処理を行う。いずれかのプロファイルの信号強度信号量に露光比に応じた係数を乗算する露光比調整を行う。この露光比調整後に、除算処理を行う。具体的な処理例について説明する。例えば、RGB中、R赤の行プロファイルのある行pの信号量を、長時間露光画像の信号量RLp、短時間露光画像の信号量RSp、さらに、長時間露光画像と短時間露光画像の露光比をn:1、とする。この場合、R赤の行pの補正係数RApは以下の式で算出する。RAp=RSp/RLp/n同様の除算処理を全ての行で算出して、1つの短時間露光画像に対応するR赤のフリッカ成分を求め、このフリッカ成分の位相反転処理によって、図8dに示すフリッカ補正係数≒フリッカ補正波形を算出することができる。この処理をRGBすべての色について行い、1つの短時間露光画像に対応するR赤、G緑、B青のフリッカ成分の算出処理によって、cフリッカ成分データを算出する。cフリッカ成分データは、a短時間露光画像に含まれるフリッカ成分のみに相当する輝度分布データとなる。さらに、図7に示すcに示すフリッカ成分データに対して、図8に示すように位相反転処理を行うことで、dフリッカ補正係数≒フリッカ補正波形が生成される。dフリッカ補正係数≒フリッカ補正波形は、a短時間露光画像に含まれるフリッカ成分の逆の強度分布を持つデータに相当する。この処理を少なくとも複数の短時間露光画像に対して実行し、複数のフリッカ補正係数≒フリッカ補正波形を算出する。図8dに示すフリッカ補正係数≒フリッカ補正波形は、フリッカのみの信号成分の逆の位相を持つ信号パターンであり、このdフリッカ補正係数≒フリッカ補正波形の行単位の補正係数を、a短時間露光画像の各行の信号画素値に乗算することで、フリッカ成分を取り除く処理が可能となる。ただし、この図8dに示すフリッカ補正係数≒フリッカ補正波形は、被写体による影響は図7の除算処理によりほぼ取り除かれているが、画像中の飽和領域、暗部領域による影響などが含まれ、このまま画像補正処理に適用してもよい補正結果は得られない。本開示の画像処理装置では、この図8dに示すフリッカ補正係数≒フリッカ補正波形に対してさらに処理を行い、最終的なフリッカ補正波形を算出する。この処理については後述する。[2-3.動被写体検出処理について]図4に示すフリッカ補正シーケンスについて、説明を続ける。ステップS102におけるフリッカ補正係数の算出処理が終了すると、次に、ステップS103に進む。ステップS103では、動被写体検出処理を実行する。これは撮影画像に動被写体が含まれるか否かの検出処理である。動被写体が含まれる場合と、含まれない場合とでは、その後の処理が異なるためである。ステップS103における動被写体検出処理は、ステップS101において求めた各画像フレームの行プロファイルの差分と予め規定した閾値との比較処理として実行する。なお、比較処理を行う行プロファイルは、フリッカの影響の少ない長時間露光時間画像とする。例えば近接する撮影時間に撮影された複数の長時間露光画像の行プロファイルを比較し、その差分データを算出する。具体的には、例えば、撮影時間txに撮影された長時間露光画像xの行プロファイルXproと、時間txの次の長時間露光撮影タイミングである撮影時間tyに撮影された長時間露光画像yの行プロファイルYproを比較して、これらの行プロファイルの差分を算出する。行プロファイルは、図5、図6を参照して説明したようにRGB各色ごとのデータとして算出する。従って、比較対象となる2つの行プロファイルXpro,Yproについて、各色、各行単位の信号量の差分を算出し、これらをすべて加算した積算値Zと、予め設定した閾値Thとを比較する。積算値Zが閾値Th以上、すなわち、Z≧Th上記式が成立した場合は、動被写体が含まれると判定する。上記式が成立しない場合は、動被写体が含まれないと判定する。[2-4.静止画と判定された画像に対するフリッカ補正処理について]次のステップS111~S116の処理は、ステップS103において動被写体が含まれないと判定された長時間露光画像の撮影タイミングの間に撮影された静止画であると判定された短時間露光画像に対する処理として実行する。まず、ステップS111において、ステップS102で算出した複数のフリッカ補正係数≒フリッカ補正波形を合成して、少なくとも一周期分以上のフリッカ補正波形を生成する。この合成処理に際しては、各画像対応のフリッカ補正係数≒フリッカ補正波形の位相を考慮して、正しく連結する処理を行う。図9に複数のフリッカ補正係数≒フリッカ補正波形の合成処理例を示す。図9には、a1フレームf1の短時間露光画像のフリッカ補正係数G緑、a2フレームf2の短時間露光画像のフリッカ補正係数G緑、これらを示している。これらは異なるタイコミングで撮影された短時間露光画像であり、フリッカの影響による行単位の信号量の波がずれている。すなわち位相差を持つデータとなっている。これらの位相差のある複数のフリッカ補正波形の部分データを、位相差を考慮して合成することで、少なくとも一周期分以上の補正波形を生成する。この合成処理によって、図9bに示すフリッカ補正波形が得られる。図9bに示すフリッカ補正波形は、フリッカによる周期的な輝度変化、すなわち蛍光灯の輝度変化パターンに対応する波形である。しかし、このフリッカ補正波形は、画像中の飽和領域、暗部領域による影響などが含まれ、純粋なフリッカに対応した補正波形とは一致しない。すなわちノイズを含むフリッカ補正波形である。図4のシーケンス図に戻り、フリッカ補正処理シーケンスについての説明を続ける。ステップS111において、図9を参照して説明したように1周期以上のフリッカ補正波形を生成すると、次に、ステップS112に進む。ステップS112では、ステップS111において生成した1周期以上のフリッカ補正波形から1周期分のフリッカ補正波形を切り出す処理を行う。図10に示すように1周期分のフリッカ補正波形を切り出す。次に、ステップS113~S114において、1周期分のフリッカ補正波形に対する離散フーリエ変換DFTを行う。さらに、DFT変換結果を用いたノイズ除去処理等を行い、正弦波近似処理を実行する。具体的には、例えば1周期分のフリッカ補正波形に対する離散フーリエ変換DFTにより、信号成分を1~N次成分のデータに区分する。次に、2~N-1次成分のデータを1次成分データとN次成分データに振り分ける。このように、2次以降の高次成分を例えば低次に振り分ける処理を行うことで、波形振幅の変動を防止してノイズと推定される成分を除去することができる。さらに、フリッカ波形の振幅調整としてのゲイン調整を実行して、ノイズを除去し正弦波近似されたフリッカ補正波形を生成する。正弦波近似を行ったフリッカ補正波形の1周期分を繰り返すことで複数周期のフリッカ補正波形が生成される。例えば図11に示すフリッカ補正波形を生成する。このフリッカ補正波形は蛍光灯の明暗周期に対応する波形である。すらわち、蛍光灯下での撮影画像の画素値に対するフリッカの影響を示すデータである。なお、このフリッカ補正波形生成処理に際して、例えば一部のデータが欠損している場合などには、例えば前後の値の平均値を利用した補間処理やスプライン補間処理、線形補間処理などを適用して欠損部分を埋める処理を行う。この正弦波近似によって得られるフリッカ補正波形は、先に図8を参照して説明した図8cのフリッカ成分データを位相反転した図8dに示すフリッカ補正係数を適用して生成する処理として説明した。しかし、図8cのフリッカ成分データの位相反転を行わずに、図8cのフリッカ成分データをそのまま利用して制限は近似を行っても、位相がずれるのみで波形としては同じ波形を得ることが可能であり、図8dの位相反転処理は行わずに図8cのフリッカ成分データを利用して図11に示す正弦波近似されたフリッカ補正波形を生成する構成としてもよい。次に、ステップS115において、ステップS114で生成したフリッカ補正波形を適用して、各短時間露光画像のフリッカ補正を実行する。フリッカ補正波形を適用したフリッカ補正処理の具体例について図12を参照して説明する。図12には、1Aフリッカ補正対象となる短時間露光画像2フリッカ補正波形=図11に示すフリッカ補正波形3フリッカ補正後の補正された短時間露光画像これらを示している。図121Aに示す短時間露光画像は、フリッカの影響により、行単位で明らかなフリッカの影響による信号値画素値の高低の波が存在する。画像の上端を行a、下端を行bとする。行a~bにはフリッカの影響と推定される明らかな信号値画素値の高低の波が存在する。このフリッカによる影響を除去するために、ステップS113~S114において生成したフリッカ補正波形を適用した信号処理を実行する。先に、図8を参照して説明したように、フリッカ補正波形は、フリッカ成分信号の逆位相の信号である。従って、フリッカの影響を受けた短時間露光画像の各行の画素値に対して、フリッカ補正係数からなるフリッカ補正波形に対応する係数を乗算すれば、フリッカの影響を排除した補正画像を生成することができる。なお、このフリッカ補正波形に対応する係数の乗算処理に際しては、フリッカ補正対象とする短時間露光画像におけるフリッカの位相と、フリッカ補正波形の位相を整合させることが必要である。この位相調整のためには、例えばステップS102において生成した短時間露光画像に対応するフリッカ補正係数データと、フリッカ補正波形のピークを合わせ、さらに、短時間露光画像の露光時間に応じてフリッカ補正波形の切り出し範囲を設定するといった処理が有効である。具体的には、図13に示すように、1A短時間露光画像に対応する1Bフリッカ補正係数データと、2フリッカ補正波形のピークを合わせ、さらに、短時間露光画像の露光時間に応じてフリッカ補正波形の切り出し範囲を設定する。このような処理によってフリッカ補正波形の切り出し範囲を決定し、切り出したフリッカ補正波形に対応する乗算係数を各行単位で設定して、1A短時間露光画像の各行の画素値に乗算して補正画素値を算出する。ステップS115では、このような処理によって、補正画像を生成する。このように、信号処理部102は、補正対象画像の行プロファイルのピーク位置と、フリッカ補正波形のピーク位置を一致させる位相合わせを行い、フリッカ補正波形に対応する行単位の係数を、短時間露光画像の各行の画素値に乗算する処理を実行して、フリッカ成分を除去した補正画像を生成する。なお、ステップS115におけるフリッカ補正処理は、同一環境の下で撮影された画像中、動被写体を含まないと判定された画像に対して、フリッカ補正波形を適用して実行することができる。なお、同じ蛍光灯の下で撮影された画像には同じ周期のフリッカパターンが生じ、1つのフリッカ補正波形を生成後は、そのフリッカ補正波形を適用して、同一環境下でその後に撮影された画像に対する補正を行うことが可能だからである。ただし、各短時間露光画像に対するフリッカ補正の適用範囲の抽出に際しては、図13を参照して説明したように、各処理対象画像の行プロファイルのピーク位置と、フリッカ補正波形のピーク一を合わせるといった位相調整を行う必要がある。ステップS114において生成したフリッカ補正波形は、ステップS116においてメモリに格納され、各撮影画像の補正に利用される。[2-5.動被写体を含む画像に対するフリッカ補正処理について]次に、図4に示すシーケンス図のステップS121~S122の処理について説明する。この処理は、ステップS103において撮影画像に動被写体が含まれると判定された場合の処理である。ステップS121~S122の処理は、ステップS103における動被写体検出処理に適用された複数の長時間画像の撮影期間に撮影された短時間露光画像に対する処理として実行される。例えば明るい被写体が画像の縦方向に動いた場合、あるいは暗い被写体が縦方向に動いた場合、撮影画像には、フリッカの影響以外に被写体の動きに応じた行単位の輝度変化が生じる。このような動被写体による輝度変化を含む場合、長時間露光画像から生成した行プロファイルを基準プロファイルとして、先に図7を参照して説明した、短時間露光画像の行プロファイル/長時間露光画像の行プロファイルこのような除算処理を実行しても、図7cフリッカ成分データや、図8dに示すフリッカ補正波形は、動被写体の移動に伴う誤差を含むデータとなってしまう。このような動被写体の移動の影響を含むフリッカ補正波形を適用して、先に図12、図13を参照して説明した補正処理、すなわち、補正対象とする動被写体を含む画像の行プロファイルのピーク位置と、静止画に基づいて生成したフリッカ補正波形のピーク位置を合わせて補正処理を行っても、動被写体の影響によって良好なフリッカ補正が行えない場合がある。ステップS121~S122の処理は、このような問題を解決するために、動被写体を含むと判断された短時間露光画像に対して行われる処理である。ステップS121では、フリッカ補正波形をメモリから読み出す。このフリッカ補正波形は、ステップS111~S114の処理において、静止画と判断された画像に基づいて生成されたフリッカ補正波形である。次に、ステップS122において、メモリから読み出したフリッカ補正波形を、動被写体撮影画像に対して適用可能なフリッカ補正波形に設定する調整処理を行う。この調整処理には、動被写体が撮影された所定期間の複数の短時間露光画像を利用する。動被写体が撮影された所定期間の複数の短時間露光画像の各々の行プロファテイルを合成して生成したフリッカ補正波形=動被写体画像フリッカ補正波形と、メモリから読み出したフリッカ補正波形を比較して、メモリから読み出したフリッカ補正波形の位相を、動被写体画像フリッカ補正波形の位相に合わせる位相制御を行う。ステップS115では、動被写体を含む短時間露光画像に対して、ステップS122で生成した位相制御されたフリッカ補正波形を適用した画像補正処理を実行する。この画像補正処理は、位相制御されたフリッカ補正波形を適用する以外は、先に図12、図13を参照して説明した処理と同様の処理である。このように、補正対象とする短時間露光画像に動被写体を含むか否かを判定し、動被写体を含むと判定した画像については、静止画像に基づいて算出したフレーム間位相ずれに基づいて、動被写体を含む画像と位相が一致するように、フリッカ補正波形の切り出しを実行して、切り出したフリッカ補正波形を適用した画像補正を行うこのようにして、フリッカ補正された出力画像202を生成して出力する。なお、前述したようにこの処理は色信号カラーチャンネル毎に独立した補正係数を算出して色信号単位で補正処理を行ってもよいし、色信号を区別することなく行単位に設定した共通の補正係数を全ての行構成画素に適用する処理を行う構成としてもよい。色信号単位で補正処理を行う場合は、色信号カラーチャンネル個別に前記行単位の信号量の積分値を算出し、色信号個別のフリッカ成分、フリッカ補正波形を算出する。色信号を区別することなく行単位に設定した共通の補正係数を利用する場合は、色信号カラーチャンネルを区別することなく行単位の信号量の積分値を算出し、行単位のフリッカ成分、フリッカ補正波形を算出する。上述したように、本開示の構成では、実際の照明の輝度変化を再現したフリッカ成分を解析し、この正確なフリッカ波形に基づくフリッカ補正係数による補正処理が実行されるので、フリッカの除去処理をより効果的に行うことが可能となる。また、フリッカ補正波形を生成するために、フリッカ補正対象となる短時間露光画像の近傍時間で撮影された長時間露光画像の行プロファイルを基準プロファイルとして利用する構成としたので、従来の過去画像の積算による基準プロファイルの生成処理が不要となり、効率的な補正処理が実現される。また、フリッカの影響を多く含む過去画像の積算による基準プロファイルを利用した場合、長周期のフリッカが除去できないという問題があったが、本開示の構成では、もともとフリッカの影響のほとんどない長時間露光画像を参照画像として利用したことで、長周期フリッカに対する除去効果も高めることが可能となる。[3.本開示の構成のまとめ]以上、特定の実施例を参照しながら、本開示の構成について詳解してきた。しかしながら、本開示の要旨を逸脱しない範囲で当業者が実施例の修正や代用を成し得ることは自明である。すなわち、例示という形態で本発明を開示してきたのであり、限定的に解釈されるべきではない。本発明の要旨を判断するためには、特許請求の範囲の欄を参酌すべきである。なお、本明細書において開示した技術は、以下のような構成をとることができる。1短時間露光画像と長時間露光画像の少なくとも2種類の異なる露光時間の撮影画像を入力して、短時間露光画像に含まれるフリッカ成分を除去するフリッカ補正処理を実行する信号処理部を有し、前記信号処理部は、前記短時間露光画像の行単位の信号量の積分値からなる短時間露光画像行プロファイルと、前記長時間露光画像の行単位の信号量の積分値からなる基準プロファイルを生成し、前記短時間露光画像行プロファイルに対する前記基準プロファイルによる除算処理を含むデータ処理によってフリッカ補正波形を生成し、前記フリッカ補正波形に対応する行単位の係数を、短時間露光画像の各行の画素値に乗算する処理によりフリッカ成分を除去した補正画像を生成する画像処理装置。2前記長時間露光画像は、画像撮影環境の輝度変化周期の整数倍にほぼ一致する露光時間の設定の下で撮影された画像である前記1に記載の画像処理装置。3前記信号処理部は、前記フリッカ補正波形の生成処理に際して、前記短時間露光画像行プロファイルに対する前記基準プロファイルによる除算処理によって得られるフリッカ成分データ、または、該フリッカ成分データの位相反転データに対する正弦波近似処理を実行して前記フリッカ補正波形を生成する前記1または2に記載の画像処理装置。4前記信号処理部は、前記フリッカ補正波形の生成処理に際して、前記短時間露光画像行プロファイルに対する前記基準プロファイルによる除算処理によって得られるフリッカ成分データ、または、該フリッカ成分データの位相反転データに対するスプライン補間処理を実行して前記フリッカ補正波形を生成する前記1~いずれかに記載の画像処理装置。5前記信号処理部は、前記フリッカ補正波形の生成処理に際して、前記短時間露光画像行プロファイルに対する前記基準プロファイルによる除算処理によって得られるフリッカ成分データ、または、該フリッカ成分データの位相反転データに対する線型補間処理を実行して前記フリッカ補正波形を生成する前記1~4いずれかに記載の画像処理装置。6前記信号処理部は、前記フリッカ補正波形の生成処理に際して、前記短時間露光画像行プロファイルに対する前記基準プロファイルによる除算処理によって得られるフリッカ成分データ、または、該フリッカ成分データの位相反転データに対する離散フーリエ変換を実行して、2次以降の高次成分を低次に振り分ける処理を行う前記1~5いずれかに記載の画像処理装置。7前記信号処理部は、フリッカ成分を除去した補正画像の生成処理に際して、補正対象画像の行プロファイルのピーク位置と、前記フリッカ補正波形のピーク位置を一致させる位相合わせ後に、前記フリッカ補正波形に対応する行単位の係数を、短時間露光画像の各行の画素値に乗算する処理を実行する前記1~6いずれかに記載の画像処理装置。8前記信号処理部は、各画像に動被写体を含むか否かを判定し、動被写体を含まないと判定した画像を適用して前記フリッカ補正波形の生成処理を行う前記1~7いずれかに記載の画像処理装置。9前記信号処理部は、補正対象とする短時間露光画像に動被写体を含むか否かを判定し、動被写体を含まないと判定した画像に対して、前記フリッカ補正波形を適用した画像補正を行う前記1~9いずれかに記載の画像処理装置。10前記信号処理部は、補正対象とする短時間露光画像に動被写体を含むか否かを判定し、動被写体を含むと判定した画像については、動被写体の影響による行プロファイルの位相変化を考慮して、前記フリッカ補正波形を修正し、修正したフリッカ補正波形を適用した画像補正を行う前記1~9いずれかに記載の画像処理装置。11前記信号処理部は、補正対象とする短時間露光画像に動被写体を含むか否かを判定し、動被写体を含むと判定した画像については、静止画像に基づいて算出したフレーム間位相ずれに基づいて、動被写体を含む画像と位相が一致するように前記フリッカ補正波形の切り出しを実行して、切り出したフリッカ補正波形を適用した画像補正を行う。さらに、上記した装置等において実行する処理の方法や、処理を実行させるプログラムも本開示の構成に含まれる。また、明細書中において説明した一連の処理はハードウェア、またはソフトウェア、あるいは両者の複合構成によって実行することが可能である。ソフトウェアによる処理を実行する場合は、処理シーケンスを記録したプログラムを、専用のハードウェアに組み込まれたコンピュータ内のメモリにインストールして実行させるか、あるいは、各種処理が実行可能な汎用コンピュータにプログラムをインストールして実行させることが可能である。例えば、プログラムは記録媒体に予め記録しておくことができる。記録媒体からコンピュータにインストールする他、LANLocalAreaNetwork、インターネットといったネットワークを介してプログラムを受信し、内蔵するハードディスク等の記録媒体にインストールすることができる。なお、明細書に記載された各種の処理は、記載に従って時系列に実行されるのみならず、処理を実行する装置の処理能力あるいは必要に応じて並列的にあるいは個別に実行されてもよい。また、本明細書においてシステムとは、複数の装置の論理的集合構成であり、各構成の装置が同一筐体内にあるものには限らない。"