(57)【特許請求の範囲】
【請求項１】
画像データにより表わされる、光源により照らされたシーンの画像から偽画像を除去または減少する方法であって、
前記画像データは、複数行の画像データからなる画像データアレイからなり、
フリッカ関数を割り出すステップと、
行単位ベースで前記画像データを前記フリッカ関数により除算し、前記画像から偽画像を除去または減少する処理ステップ、を含み、
前記フリッカ関数は、前記光源のフリッカ周波数とフリッカ振幅とフリッカ位相の関数であり、
前記フリッカ関数の割り出しステップが、フレーム周期と画像フレームごとの交流線の繰り返し周波数との積からフリッカ周波数を割り出すことを含む
ことを特徴とする方法。
【請求項２】
請求項１記載の方法であって、
前記フリッカ関数の割り出しステップはさらに、前記画像データアレイの各行沿いに前記画像データアレイの画素を寄せ集めて列ベクトルを生成するステップを含み、
さらに前記補正列ベクトルのフリッカエネルギを最小化することで前記フリッカ位相と前記フリッカ振幅を割り出すステップを含む方法。
【請求項３】
請求項２記載の方法であって、
補正列ベクトルのフリッカエネルギの前記最小化ステップが、フリッカの振幅及び位相の推定ステップと、割り出されたフリッカ周波数と推定されたフリッカの振幅及び位相とを用いて計算したフリッカ関数により生成列ベクトルを除算することで補正列ベクトルを計算するステップと、補正列ベクトルのフリッカエネルギを計算するステップと、フリッカエネルギがさらに低減できるかどうか判定するステップと、フリッカエネルギがさらに低減できない場合は、推定されたフリッカの振幅及び位相を用いてフリッカ関数を割り出し、フリッカエネルギがさらに低減できる場合は、フリッカエネルギがこれ以上低減できなくなるまで、フリッカの振幅及び位相を推定し、補正列ベクトルを計算し、フリッカエネルギを計算するステップを繰り返すステップを含む方法。
【請求項４】
前記光源は、周期的に変動する光源からなる請求項１〜３のいずれか１つに記載の方法。
【請求項５】
前記周期的に変動する光源は、蛍光源からなる請求項４記載の方法。
【請求項６】
前記画像データアレイの行の行番号で表わされる前記フリッカ関数が、
ｆ（ｊ）＝１＋αsin2（ｊωＴf／Ｎ＋Φ）
であり、αがフリッカ振幅、ωがフリッカ周波数、Φがフリッカ位相、Ｔfがフレーム周期、Ｎが画像フレームごとの行数を表わす請求項５記載の方法。
【請求項７】
前記画像データは回転シャッタを用いてＣＭＯＳ画像センサから収集され、露光制御をもたらす請求項１〜６のいずれか１つに記載の方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、概ね画像処理分野に係り、より具体的には周期的に変動する光源によって生ずる偽画像を画像から除去する方法及び装置に関する。
【０００２】
【従来の技術】
ＣＭＯＳ（相補型金属酸化膜半導体）のアクティブ画素方式画像センサでは、露光制御に一般に回転シャッタが用いられる。回転シャッタは、画像センサ列のセンサ各行を順次リセットし、後刻各行の値を順次読み出すことで稼働する。露光期間はリセットと読出しの時間間隔に等しく、フレーム周期は配列内の各行の値を読み出すのに必要な時間となる。
【０００３】
図１は、本発明を明瞭に理解して貰えるよう、回転シャッタの動作について概略図示した線図である。１００行（図１には説明の便宜上９行だけを図示）を有する画像データアレイについて検討する。単一行の読み出しに必要な時間、すなわち行期間をＴとし、図１の矢印１０で示す露光期間をｎ×Ｔに設定したとする。行１の露光は時間ｔ＝０にて始まり、その行は時間ｔ＝ｎ×Ｔで読み出される。行２の露光は時間ｔ＝Ｔで開始され、行２は時間ｔ＝（ｎ＋１）×Ｔにて読み出される。配列の行１００は時間ｔ＝９９×Ｔで露光を開始し、その行は時間ｔ＝（ｎ＋９９）×Ｔにて読み出される。
【０００４】
回転シャッタの一つの特性は、図１に示す如く、各行ごとの露光が異なる時間期間に亙って起こることである。従って、シーンの照明が時間とともに変動する場合、撮像された画像の輝度は行ごとに変動する。シーンの照度レベルにおける変動期間が回転シャッタの集積期間より長くてフレーム周期よりも短い場合、上記変動は画像内に水平帯の形で偽画像(image artifacts)が現われる結果を招く。画像内の帯の数は、フレーム周期で割ったフリッカ周期に等しい。
【０００５】
光源の照度レベルにおけるこれらの時間に対する変動は、「フリッカ」として知られており、画像の水平帯は偽画像或いは「フリッカパターン」と呼ばれる。例えば、蛍光灯は非常に強いフリッカを有する。ランプの蛍光体は、電力線周波数の各ピークで励起され、ピーク間で減衰する。生成される光パルスは、交流電力線の２倍のフリッカ周波数を有する。白熱バルブはフリッカを示さないが、それはフィラメントの熱量が十分に大きく、フィラメントが交流ピーク間でかなり冷めるのを防止するからである。
【０００６】
ＣＭＯＳ画像センサの出力へのフリッカの影響は、図２の概略画像１４に例示してあり、同図では周期的な水平バー１６が非常に明白であり、画像の全般的な品質を相当に損なうものである。
【０００７】
図２に示したような画像内の水平バーの数は、フリッカ周期に対するフレーム周期の比に等しい。水平バーの振幅は、特定光源の特性とＣＭＯＳ画像センサの回転シャッタの集積期間に依存する。
【０００８】
幾つかの方法が、フレームごとのフリッカを除去するのに考案されてきた。例えば、米国特許第５，７９３，８８６号には異なる画像フレームのヒストグラムを比較し、これにより時間の関数としての平均照度の変動を取り除くようにした方法が開示されている。しかしながら、ここに記述された方法は、単一の画像フレームしか利用できず、しかもフリッカが生成画像内の一連の水平帯として明示される状況には適用し得ないのである。
【０００９】
フレームごとのフリッカを取り除く他の公知の方法は、光源の輝度をフレーム露光期間中に監視するものである。例えば、米国特許第６，２０８，４３３号に記載された方法では、光源の変動監視用に別個のセンサが用いられ、このセンサの出力が画像補正に用いられる。米国特許第５，９６０，１５３号では、撮像と光源変動の計測の両方にリニア画像センサが用いられる。光源変動を露光時間制御に用い、これにより最終画像にフリッカが現われないようにしてある。
【００１０】
さらに別の公知の方法では、露光期間をフリッカ周期の整数倍に制限することでフリッカを防止している。この種の状況下では、各行は光源変動の全周期の整数倍で光を集積する。その結果、集積された光レベルは各行ごとに同一となり、画像に帯が現われることはなくなる。本方法の欠点は、露光期間をフリッカ期間よりも短くできないことにある。光源が非常に明るく、露光期間が低減できない場合、生成画像は露出過多となる。この種の問題は、明るい蛍光源でシーンを照らすか、カメラを直接蛍光源へ向けたときに生ずることがある。
【００１１】
【発明が解決しようとする課題】
一般に、画像からフレームごとのフリッカを除去する公知の方法は、処理用に幾つかの連続画像フレームを必要としたり、フリッカを独立に監視する必要があったり、追加の入力を必要としたりする。従って、フリッカ除去用の公知の方法は完全に満足は出来ないのである。必要なものは、単一画像上に作用する周期的に変動する光源により生成されたフリッカを除去し、フリッカ周期とフレーム周期の値を除きどんな追加の入力も必要とせず、フリッカを独立して監視する必要のない方法と装置である。同じく必要とされるのは、非常に明るい蛍光照明下でも非常に低いフリッカをもって画像を撮像することに活用できる方法と装置である。
【００１２】
【課題を解決するための手段】
本発明は、蛍光源のような周期的に変動する光源により生成される画像から偽画像を除去する方法と装置を提供するものである。本発明による方法は、画像データにより表される画像から偽画像を除去するものである。本方法では、フリッカ関数を割り出し、このフリッカ関数を用いて画像データを処理して偽画像を除去する。
【００１３】
本発明にあっては、撮像シーンを照明する光源の光照射をモデル化したフリッカ関数を用いて画像を処理することにより画像から偽画像を除去或いは大幅に低減できることが判った。本発明は、回転シャッタを用いてＣＭＯＳ画像センサから収集した画像データにより表される画像から偽画像を除去するのに特に適したものである。
【００１４】
本発明の一実施形態によれば、画像を表わす画像データは、複数行の画像データからなる画像データアレイからなり、画像データの処理ステップは、行単位ベースで画像データをフリッカ関数により除算することからなる。
【００１５】
画像データのフリッカ関数の割り出しは、前以て用意したフリッカモデルに関連するパラメータ値の割り出しを含む。特に、フリッカ関数を割り出すため光源のフリッカ周波数やフリッカ振幅やフリッカ位相が割り出される。本発明の一実施形態では、フリッカ周波数とフリッカ振幅とフリッカ位相は、画像データアレイの各行沿いに画素を寄せ集めて列ベクトルを生成し、フレーム周期と画像フレームごとの交流線繰り返し周波数との積からフリッカ周波数を割り出し、補正列ベクトル内のフリッカエネルギを最小化してフリッカの振幅及び位相を割り出すことで割り出される。
【００１６】
より詳しくは、本発明のさらなる実施形態によれば、フリッカの振幅及び位相の初期推定を行ない、割り出されたフリッカ周波数と推定されたフリッカの振幅及び位相から計算したフリッカ関数により生成列ベクトルを除算して補正列ベクトルを計算し、補正列ベクトル内のフリッカエネルギを算出することにより、フリッカの振幅及び位相を割り出す。そこで、フリッカエネルギがさらに低減可能であるかどうか判断する。判断が否定された場合は、推定されたフリッカの振幅及び位相をフリッカ関数の割り出しに用いる。判断が肯定された場合は、フリッカの振幅及び位相の推定は改められ、さらなる補正列ベクトルが計算され、さらに補正された列ベクトルを用いてフリッカエネルギを計算する。この手順は、フリッカエネルギがこれ以上低減できなくなるまで繰り返される。
【００１７】
一旦フリッカエネルギが最小化されると、割り出されたフリッカ関数でもって画像データを除算することでフリッカが画像データから取り除かれる。
【００１８】
本発明方法は、蛍光源のような任意の周期的に変動する光源にもよってもたらされるフリッカを除去することに効果があるという長所がある。本発明方法には任意の露光期間を用いることができ、その結果、本方法は明るい蛍光源の下でさえフリッカを大幅に低減した画像の撮像に用いることができる。本発明は、単一画像上に作用する偽画像を画像から除去し、フリッカとフレーム周期の値以外のどんな追加の入力も必要とせず、フリッカを独立して監視する必要のない方法と装置を提供するものである。さらに、本発明は上記したそれらに加え或いはそれらに代わる他の特徴ならびに効果を備えた実施形態をもたらす。これらの特徴と効果の多くは、本発明の例示実施形態の下記の図面と詳細な説明とを併せ、以下に明らかとなろう。
【００１９】
【発明の実施の形態】
回転シャッタモードで動作するＣＭＯＳ画像センサでは、画像データアレイの露光期間は、図１に示したアレイ内で行ごとに変動する（ただし、特定の行では全ての画素は同一期間に亙って露光する）。撮像シーンを照らす光源の強度が時間とともに変動する場合は、画像データ内で行ごとに強度変動が生ずることがある。行期間が固定され、光源強度が時間とともに周期的に変動する場合は、強度変動は画像データから生成される画像内の周期的な水平バーの形をした偽画像として現われる。
【００２０】
集積期間がフリッカ周期の整数倍に等しい場合は、水平バーの振幅は零となろう。しかしながら、バー振幅を前以て割り出すことは一般に不可能であり、それは光源の特性が既知でないからである。加えて、画像についての水平バーパターンの垂直方向の並び方は露光期間当初の光フリッカの位相に依存する。しかしながら、振幅と同様、位相もまた前以て割り出すことはできない。
【００２１】
蛍光の点滅は、光の電源交播によって生じる。蛍光管は、水銀蒸気原子雲に電子を放射するフィラメントを含む。蒸気原子は励起され、紫外（ＵＶ）光を発する。紫外光が蛍光管の管球容器を叩くと、紫外光が管内面の蛍光体を励起する。そこで蛍光体が蛍光し、白色光を発する。
【００２２】
電子の放出は、電源の駆動電圧のピークで最大化される。１周期ごとに２つのピークが存在するため、フリッカ周波数は駆動電圧周波数の２倍となる。従って、６０Ｈｚ電気系により給電される蛍光管は１２０Ｈｚで点滅する。
【００２３】
関数（１＋αｓｉｎ２（ωｔ＋Φ））は、蛍光管の照射をモデル化するのに用いることができる。同式中、αはフリッカ振幅であり、ωはフリッカ周波数であり、Φはフリッカ位相である。このフリッカ関数はまた、画像アレイの行の行番号をもって次式の如く表わすことができる。
【００２４】
ｆ（ｊ）＝１＋αｓｉｎ２（ｊωＴｆ／Ｎ＋Φ）                    （１）
【００２５】
同式中、Ｔｆはフレーム周期であり、Ｎはフレームごとのライン数である。蛍光管の実際の発光は、図３の曲線２０に示すもっと複雑な時間変動を有する。同図は、ｘ軸が時間サンプルを表わし、ｙ軸が照度を表わすグラフである。蛍光管内の蛍光体とフィラメントの特性が原因で、この関数はゆっくりと減衰する。従って、所望とあらば、上記の式（１）により与えられるもっと複雑なフリッカ関数を生成し、曲線にもっと適合させることもできる。
【００２６】
本発明の一実施形態によれば、行単位ベースのフリッカ関数でもって画像を表わす画像データを除算することにより画像から偽画像が取り除かれる。画像が二次元データアレイＩ［ｉ，ｊ］により表される場合、偽画像を除去する処理は下記の如く表わされる。
【００２７】
Ｉ’［ｉ，ｊ］＝Ｉ［ｉ，ｊ］／ｆ（ｊ）                          （２）
【００２８】
本発明の一実施形態によれば、フリッカ関数は照明光源のフリッカ周波数とフリッカ振幅とフリッカ位相を算出することにより割り出される。フリッカ周波数とフリッカ振幅とフリッカ位相は、本発明の一実施形態によれば、デジタル画像の各行沿いに画素を寄せ集めて列ベクトルを生成し、フレーム周期と画像フレームごとの交流線繰り返し周波数との積からフリッカ周波数を算出し、補正列ベクトル内のフリッカエネルギを最小化してフリッカの振幅及び位相を割り出すことで算出される。
【００２９】
より詳しくは、列ベクトルを作成するため、画像データアレイの全ての列は併せ加算されて列平均を生成する。この処理が、二次元画像データアレイを一次元ベクトルに折りたたむ。この処理は、下式により表わされる。
【００３０】
【数１】
0004200516.tif 000002
（３）
【００３１】
この式において、Ｍは画像データの列数である。
【００３２】
フリッカの振幅及び位相を算出するため、補正列ベクトル内のフリッカエネルギを最小化するよう、初回の推定が先ずフリッカの振幅及び位相についてなされる。次に、割り出されたフリッカ周波数と推定されたフリッカの振幅及び位相から算出されたフリッカ関数により生成列ベクトルを除算することで、補正列ベクトルが計算される。そこで補正列ベクトル内でフリッカエネルギが計算され、フリッカエネルギがそれ以上低減できるかどうかの判断がなされる。この判断が否定された場合は、推定されたフリッカの振幅及び位相がフリッカ関数の割り出しに用いられる。判断が肯定された場合は、フリッカの振幅及び位相の改善推定がなされ、フリッカエネルギがそれ以上低減できなくなるまで、この手順が繰り返される。
【００３３】
補正列ベクトルのフリッカエネルギは、Ｃ［ｊ］／ｆ（ｊ）の正弦値及び余弦値を計算することにより算出される。エネルギは、画像データの平均値により正規化した２つの項の二乗の平方根である。
【００３４】
【数２】
0004200516.tif 000003
0004200516.tif 000004
0004200516.tif 000005
0004200516.tif 000006
【００３５】
式（４），（５）中、Ａ（ｊ）は、画像フレーム端での不連続部分がエネルギ値に寄与するのを防止するのに用いるアポディゼーション関数である。この種の関数例には、ハニング関数やハミング関数が含まれる。
【００３６】
フリッカエネルギは、下記で表わされる。
【００３７】
Ｅｆ＝√（ＸＳ＋ＸＣ）／＜Ｃ＞                                  （６）
【００３８】
式（６）では、平均画像データは列平均から割り出してある。
【００３９】
振幅と位相の関数としてのフリッカエネルギの変動により規定される面が図４にプロットしてあり、図ではｘ軸が振幅を表わしｙ軸が位相を表わす。図４に示したように、この関数は良好に規定されたユニークな最小値を備えた滑らかなものである。
【００４０】
一旦フリッカエネルギが最小化されると、フリッカ周波数の計算値と最小化から得られたフリッカの位相及び振幅の計算値を用いて割り出したフリッカ関数により画像データを除算することで、全画像から偽画像を取り除くことができる。
【００４１】
図５，６は、本発明の一実施形態になる方法の効果を示す画像である。図５の画像は、ＣＭＯＳ画像センサすなわち原画像の単色平面からの生データで構成してある。図６の画像は、本発明の一実施形態になるフリッカ補正、すなわちフリッカ補正画像を適用した結果である。図５，６から明らかなように、図５の画像に存在する水平線が図６の画像からは除去されている。
【００４２】
図７は、本発明の一実施形態によるシーンの画像から偽画像を取り除く方法を大略例示するフローチャートである。図７に示すように、本方法はステップ３０でのシーンを照らす光源の発光をモデル化するフリッカ関数の割り出しを含むものであり、その後にステップ４０でのフリッカ関数を用いてシーンの画像を表わす画像データの処理を含む。
【００４３】
図８Ａは、本発明の別の実施形態により図７のフリッカ関数を割り出すステップをより詳細に例示するフローチャートである。特に、前述した如く、撮像対象シーンを照明する光源のフリッカ関数を割り出すため、光源のフリッカ周波数とフリッカ振幅とフリッカ位相が割り出される。これは、ステップ５０に示したように先ず画像アレイの各行沿いに画像アレイの画素を寄せ集めて列ベクトルを生成することにより、図８Ａに開示した如く達成される。フリッカ周波数は、ステップ６０に示したように、フレーム周期と画像フレームごとの交流線繰り返し周波数との積から割り出される。フリッカ位相とフリッカ振幅は、ステップ７０で示すように補正列ベクトルのフリッカエネルギを最小化することにより割り出される。
【００４４】
図８Ａに示されるフリッカ位相とフリッカ振幅を割り出すステップを、図８Ｂでさらに詳細に説明する。具体的には、フリッカの振幅及び位相の初期推定はステップ７２に示した如くなされる。それ故に、補正列ベクトルは図８Ａのステップ５０で生成された列ベクトルをステップ７４の推定フリッカ関数で除算することで計算される。この推定フリッカ関数は、図８Ａのステップ６０で割り出したフリッカ周波数と図８Ｂのステップ７２で推定したフリッカ振幅とフリッカ位相を用いて計算される。補正列ベクトルのフリッカエネルギは次にステップ７６において計算され、ステップ７８においてフリッカエネルギがさらに低減できるかどうかの判断がなされる。この判断が否定された場合は、ステップ７２においてなされたフリッカ振幅とフリッカ位相の初期推定を用いてフリッカ関数が割り出される。判断が肯定された場合は、ステップ８０においてフリッカ振幅とフリッカ位相の改善推定がなされ、ステップ７４においてさらなる補正列ベクトルが計算され、ステップ７６において新たなフリッカエネルギが計算される。この手順は、フリッカエネルギがこれ以上低減できなくなるまで、すなわちステップ７８における判断がなされなくなるまで繰り返される。
【００４５】
図９は、本発明の一実施形態によりシーンの画像から偽画像を取り除く装置を図示するものである。この装置は、参照符合１００により概ね指し示したディジタルカメラなどの撮像装置を備える。撮像装置は、シーンから受光し、ＣＭＯＳ画像センサなどの画像検出ユニット１０４に光を合焦する図９のレンズにより表わされる光学システム１０２を含む。画像検出ユニットは、その上の光パターンを画像検出ユニット内に含まれるＡＤ変換器によりディジタル化される電圧サンプルアレイへ変換する。
【００４６】
フリッカ関数割り出し器１０６は、前述したように、光源のフリッカ周波数とフリッカ振幅とフリッカ位相を割り出すことでディジタル画像データから光源のフリッカ関数を割り出し、画像データプロセッサ１０８はフリッカ関数を用いて画像データを処理し、シーンの生成画像から偽画像を取り除く。
【００４７】
ここに説明した内容は本発明の例示実施形態を構成するものであるが、本発明はその範囲を逸脱することなく様々に変形可能であることは認識されたい。従って、本発明が特許請求の範囲により必要とされる範囲内にだけ限定さるべきことは理解されたい。
【００４８】
本発明の態様を以下に例示する。
【００４９】
１．画像データにより表わされる、光源により照らされたシーンの画像から偽画像を除去する方法であって、フリッカ関数を割り出すステップと、前記フリッカ関数を用いて前記画像データを処理し、前記画像から偽画像を除去するステップを含む方法。
【００５０】
２．上記１記載の方法であって、前記画像データは複数行の画像データからなる画像データアレイからなり、前記処理ステップが行単位ベースで前記画像データを前記フリッカ関数により除算するステップを含む方法。
【００５１】
３．上記２記載の方法であって、前記フリッカ関数は、前記光源のフリッカ周波数とフリッカ振幅とフリッカ位相の関数である方法。
【００５２】
４．上記３記載の方法であって、前記フリッカ関数の割り出しステップが、フレーム周期と画像フレームごとの交流線の繰り返し周波数との積からフリッカ周波数を割り出すことを含む方法。
【００５３】
５．上記４記載の方法であって、前記フリッカ関数の割り出しステップはさらに、前記画像データアレイの各行沿いに前記画像データアレイの画素を寄せ集めて列ベクトルを生成するステップを含み、さらに前記補正列ベクトルのフリッカエネルギを最小化することで前記フリッカ位相と前記フリッカ振幅を割り出すステップを含む方法。
【００５４】
６．上記５記載の方法であって、補正列ベクトルのフリッカエネルギの前記最小化ステップが、フリッカの振幅及び位相の推定ステップと、割り出されたフリッカ周波数と推定されたフリッカの振幅及び位相とを用いて計算したフリッカ関数により生成列ベクトルを除算することで補正列ベクトルを計算するステップと、補正列ベクトルのフリッカエネルギを計算するステップと、フリッカエネルギがさらに低減できるかどうか判定するステップと、フリッカエネルギがさらに低減できない場合は、推定されたフリッカの振幅及び位相を用いてフリッカ関数を割り出し、フリッカエネルギがさらに低減できる場合は、フリッカエネルギがこれ以上低減できなくなるまで、フリッカの振幅及び位相を推定し、補正列ベクトルを計算し、フリッカエネルギを計算するステップを繰り返すステップを含む方法。
【００５５】
７．前記光源は、周期的に変動する光源からなる上記２記載の方法。
【００５６】
８．前記周期的に変動する光源は、蛍光源からなる上記７記載の方法。
【００５７】
９．前記画像データアレイの行の行番号で表わされる前記フリッカ関数が、
ｆ（ｊ）＝１＋αｓｉｎ２（ｊωＴｆ／Ｎ＋Φ）
であり、αがフリッカ振幅、ωがフリッカ周波数、Φがフリッカ位相、Ｔｆがフレーム周期、Ｎが画像フレームごとの行数を表わす上記８記載の方法。
【００５８】
１０．前記画像データは回転シャッタを用いてＣＭＯＳ画像センサから収集され、露光制御をもたらす上記１記載の方法。
【００５９】
１１．複数行の画像データからなる画像データアレイにより表わされる、周期的に変動する光源により照明されたシーンの画像から偽画像を除去する方法であって、周期的に変動する光源の光照射をモデル化したフリッカ関数を割り出すステップと、前記フリッカ関数を用いて前記画像データを処理し、前記画像から偽画像を除去するステップを含む方法。
【００６０】
１２．上記１１記載の方法であって、前記フリッカ関数は、周期的に変動する光源のフリッカ振幅とフリッカ周波数とフリッカ位相の関数である方法。
【００６１】
１３．上記１２記載の方法であって、フリッカ関数を割り出す前記ステップが、画像データアレイの各行沿いに画素を寄せ集めて列ベクトルを生成し、フレーム周期と画像フレームごとの交流線の繰り返し周波数との積からフリッカ周波数を割り出すステップと、補正列ベクトルのフリッカエネルギを最小化することによりフリッカの振幅及び位相を割り出すステップとを含む方法。
【００６２】
１４．上記１３記載の方法であって、補正列ベクトルのフリッカエネルギを最小化する前記ステップが、フリッカの振幅及び位相の推定ステップと、割り出されたフリッカ周波数と推定されたフリッカの振幅及び位相とを用いて計算したフリッカ関数により生成列ベクトルを除算することで補正列ベクトルを計算するステップと、補正列ベクトルのフリッカエネルギを計算するステップと、フリッカエネルギがさらに低減できるかどうか判定するステップと、フリッカエネルギがさらに低減できない場合は、推定されたフリッカの振幅及び位相を用いてフリッカ関数を割り出し、フリッカエネルギがさらに低減できる場合は、フリッカエネルギがこれ以上低減できなくなるまで、フリッカの振幅及び位相を推定し、補正列ベクトルを計算し、フリッカエネルギを計算するステップを繰り返すステップを含む方法。
【００６３】
１５．上記１１記載の方法であって、前記周期的に変動する光源は蛍光源からなり、前記画像データアレイの行の行番号で表わされる前記フリッカ関数が、
ｆ（ｊ）＝１＋αｓｉｎ２（ｊωＴｆ／Ｎ＋Φ）
であり、αがフリッカ振幅、ωがフリッカ周波数、Φがフリッカ位相、Ｔｆがフレーム周期、Ｎが画像フレームごとの行数を表わす方法。
【００６４】
１６．上記１１記載の方法であって、前記画像データは回転シャッタを用いてＣＭＯＳ画像センサから収集され、画像制御をもたらす方法。
【００６５】
１７．光源によって照らされたシーンの画像から偽画像を除去する装置であって、前記シーンを表わす画像データを供給するユニットと、フリッカ関数を割り出すフリッカ関数割り出し器と、前記フリッカ関数を用いて前記画像データを処理し、前記画像から前記偽画像を除去する画像データプロセッサを備える装置。
【００６６】
１８．上記１７記載の装置であって、前記画像データが複数行の画像データからなる画像データアレイを含み、前記画像データプロセッサが行単位ベースで前記画像データを前記フリッカ関数により除算することで前記画像データを処理する装置。
【００６７】
１９．上記１８記載の装置であって、前記フリッカ関数は、前記光源のフリッカ周波数とフリッカ振幅とフリッカ位相の関数であり、前記フリッカ関数割り出し器が前記フリッカ周波数と前記フリッカ振幅と前記フリッカ位相を割り出す手段を含む装置。
【００６８】
２０．  前記光源は、蛍光源からなる上記１７記載の装置。
【００６９】
２１．ディジタル撮像方法であって、ディジタル画像を分析してフリッカ関数を割り出すステップと、前記フリッカ関数を用いて前記ディジタル画像を処理し、該ディジタル画像から偽画像を低減するステップとを含み、前記割り出しステップが前以て定めたフリッカモデルに関連するパラメータ値の割り出しを含む方法。
【００７０】
２２．前記処理ステップが、前記ディジタル画像を前記フリッカ関数により除算するステップを含む上記２１記載のディジタル撮像方法。
【００７１】
２３．前記割り出しステップが、フリッカ周波数とフリッカ振幅とフリッカ位相の算出を含む上記２１記載のディジタル撮像方法。
【００７２】
２４．上記２３記載のディジタル撮像方法であって、フリッカ周波数とフリッカ振幅とフリッカ位相の前記割り出しステップが、画像データアレイの各行沿いに画素を寄せ集めて列ベクトルを生成し、フレーム周期と画像フレームごとの交流線の繰り返し周波数との積からフリッカ周波数を割り出すステップと、補正列ベクトルのフリッカエネルギを最小化することによりフリッカの振幅及び位相を割り出すステップとを含む方法。
【００７３】
２５．上記２４記載のディジタル撮像方法であって、前記補正列ベクトルのフリッカエネルギを最小化するステップが、フリッカの振幅及び位相を推定するステップと、割り出されたフリッカ周波数と推定されたフリッカの振幅及び位相を用いて計算したフリッカ関数により生成列ベクトルを除算することにより補正列ベクトルを計算するステップと、補正列ベクトルのフリッカエネルギを算出するステップと、フリッカエネルギがさらに低減できるかどうか判定するステップと、フリッカエネルギがさらに低減できない場合は、推定されたフリッカの振幅及び位相を用いてフリッカ関数を割り出し、フリッカエネルギがさらに低減できる場合は、フリッカエネルギがこれ以上低減できなくなるまで、フリッカの振幅及び位相を推定し、補正列ベクトルを計算し、フリッカエネルギを計算するステップを繰り返すステップを含む方法。
【００７４】
２６．ディジタル画像処理システムであって、ディジタル画像のフリッカ関数を割り出す手段と、前記フリッカ関数を用いて前記ディジタル画像内の偽画像を低減する該ディジタル画像の処理手段とを備えるシステム。
【図面の簡単な説明】
【図１】本発明の明瞭な理解を助ける回転シャッタの動作を概略説明する図である。
【図２】回転シャッタを用いてＣＭＯＳ画像センサの出力上のフリッカの影響を例示する概略画像である。
【図３】蛍光管の発光特性を例示するグラフである。
【図４】本発明の一実施形態になるフリッカ補正を例示するグラフである。
【図５】本発明の一実施形態の方法と装置の効果を表わす図であり、回転シャッタを用いたＣＭＯＳ画像センサの単色平面からの生データの画像である。
【図６】本発明の一実施形態の方法と装置の効果を表わす図であり、本発明の一実施形態によりフリッカ補正を適用した結果得られる画像である。
【図７】本発明の一実施形態による偽画像を画像から除去する方法のステップを示すフローチャートである。
【図８Ａ】本発明の他の実施形態による図７のフリッカ関数の割り出しステップを詳細に示すフローチャートである。
【図８Ｂ】本発明の他の実施形態による図８Ａのフリッカの位相及び振幅の割り出しステップを詳細に示すフローチャートである。
【図９】本発明の他の実施形態による偽画像を画像から除去する装置を概略示す図である。
【符号の説明】
１０  露光期間
１４  概略画像
１６  水平バー
１００  撮像装置
１０２  光学システム
１０４  画像検出ユニット
１０６  フリッカ関数割り出し器
１０８  画像データプロセッサ
【図１】
0004200516.tif 000007
【図２】
0004200516.tif 000008
【図３】
0004200516.tif 000009
【図４】
0004200516.tif 000010
【図５】
0004200516.tif 000011
【図６】
0004200516.tif 000012
【図７】
0004200516.tif 000013
【図８Ａ】
0004200516.tif 000014
【図８Ｂ】
0004200516.tif 000015
【図９】
0004200516.tif 000016