JPWO2017221376A1 - 画像処理装置、画像処理方法および画像処理プログラム - Google Patents

Abstract

少ない計算量で精度よく欠陥画素を検出することを目的として、本発明に係る画像処理装置（１）は、取得された画像の注目画素と周辺画素とを含むブロック内の各画素の画素値に対して、順序統計フィルタ処理を施すフィルタ処理部（２）と、フィルタ処理部（２）により順序統計フィルタ処理が施される前後の注目画素の画素値の差分値を算出する差分値算出部（３）と、フィルタ処理部（２）により順序統計フィルタ処理が施された周辺画素の画素値の統計的バラツキを算出するバラツキ算出部（４）と、差分値算出部（３）により算出された差分値をバラツキ算出部（４）により算出された統計的バラツキにより除算して得られた値が所定の閾値以上である場合に、注目画素を欠陥画素であると判定する判定部（５）とを備える。

Description

本発明は、画像処理装置、画像処理方法および画像処理プログラムに関するものである。

デジタルカメラ等の撮影装置を用いて取得された画像において、注目画素とその周囲の８画素の信号レベルの分散値の平方根、すなわち、標準偏差を算出し、注目画素の信号レベルがその周囲の信号レベルと比較してどれだけ突出しているのかを示す突出度を算出し、突出度が所定の上限値より大きい場合に、注目画素が白欠陥、所定の下限値より小さい場合に、注目画素が黒欠陥であると判定する画像処理方法が知られている（例えば、特許文献１参照。）。

特許第４４５３３３２号公報

しかしながら、特許文献１の画像処理方法は、突出度を算出するために周辺画素の画素値の平均値および標準偏差を算出しているため、周辺画素に複数の欠陥画素が存在する場合には、平均値および標準偏差が欠陥画素の影響を受けて正しい値を求めることができない。その結果、欠陥検出の精度が低下するという不都合がある。また、ハードウェアへの実装を考えると、標準偏差の計算は計算量が多いため、処理に時間がかかったり、消費電力が増大したりするという不都合がある。

本発明は上述した事情に鑑みてなされたものであって、少ない計算量で精度よく欠陥画素を検出することができる画像処理装置、画像処理方法および画像処理プログラムを提供することを目的としている。

本発明の一態様は、取得された画像の注目画素と周辺画素とを含むブロック内の各画素の画素値に対して、順序統計フィルタ処理を施すフィルタ処理部と、該フィルタ処理部により順序統計フィルタ処理が施される前後の前記注目画素の前記画素値の差分値を算出する差分値算出部と、前記フィルタ処理部により順序統計フィルタ処理が施された前記周辺画素の前記画素値の統計的バラツキを算出するバラツキ算出部と、前記差分値算出部により算出された差分値を前記バラツキ算出部により算出された統計的バラツキにより除算して得られた値が所定の閾値以上である場合に、前記注目画素を欠陥画素であると判定する判定部とを備える画像処理装置である。

本態様によれば、取得された画像に対しフィルタ処理部において順序統計フィルタ処理が施され、順序統計フィルタ処理前後の注目画素の画素値の差分値が差分値算出部により算出される。一方、順序統計フィルタ処理が施された周辺画素の画素値の統計的バラツキがバラツキ算出部により算出され、差分値を統計的バラツキで除算して得られた値が所定の閾値以上であるか否かが判定部により判定される。判定の結果閾値以上である場合には、注目画素が欠陥画素であると判定される。

この場合において、本態様によれば、注目画素が欠陥画素であるか否かを平均値を用いることなく判定できる。したがって、同じブロック内に複数の欠陥画素が存在しても、欠陥画素の影響を受けることがなく、また、ノイズが存在していてもノイズの影響を低減することができ、欠陥検出の精度を向上することができる。特に、標準偏差を算出しないので少ない計算量で精度よく欠陥画素を検出することができる。

また、単純に順序統計フィルタ処理の前後の画素値の差分値を閾値によって評価するのではなく、順序統計フィルタ処理後の周辺画素から算出した統計的バラツキとの比率により評価することにより、注目画素がエッジに隣接もしくはテクスチャ上に存在する場合の誤検出を抑制することができる。

上記態様においては、前記バラツキ算出部が、統計的バラツキとして四分位範囲を算出してもよい。
これにより、簡易に統計的バラツキを算出することができる。

また、上記態様においては、前記判定部により前記欠陥画素であると判定された画素の座標値を登録するアドレス登録部と、該アドレス登録部に登録された座標値に対応する画素の前記欠陥画素の評価値を記憶する評価値記憶部と、前記差分値算出部により算出された差分値と、前記評価値記憶部に記憶されている過去の前記評価値とに基づいて前記欠陥画素の前記評価値を算出する評価値算出部とを備え、前記アドレス登録部は、前記評価値算出部により算出された前記評価値が、所定の閾値以下である場合に登録を解除してもよい。

このようにすることで、誤検出によって欠陥画素であると判定された場合であっても、欠陥画素ではない場合には、処理の進行に従って欠陥画素であるとの判定を解除することができる。したがって、画像の画素数と同数アドレス数を検出用メモリに用意する必要がなくなる。

また、上記態様においては、前記評価値算出部は、前記差分値と過去の前記評価値とを重み付け合成して前記評価値を算出してもよい。
このようにすることで、正常な画素が誤検出により欠陥画素であると判定された場合であっても、処理の進行に従って欠陥画素であるとの判定を解除することができる。

また、上記態様においては、前記判定部により前記欠陥画素であると判定された画素の前記座標値の前記アドレス登録部への登録を間引く間引き部を備えていてもよい。
このようにすることで、欠陥画素であるとの誤検出数を低減でき、欠陥画素のアドレス登録用のメモリ量をさらに低減することができる。

また、上記態様においては、前記評価値算出部は、前記バラツキ算出部が算出した順序統計フィルタ処理が施された前記周辺画素の前記画素値の統計的バラツキが所定の閾値より大きい場合に、前記評価値の算出を行わないこととしてもよい。
このようにすることで、順序統計フィルタ処理が施された周辺画素の画素値の統計的バラツキが所定の閾値より大きい場合には、画素値の分布が平坦ではないと判断でき、評価値の算出を行わないことにより、欠陥画素の誤検出を防止することができる。

本発明の他の態様は、取得された画像の注目画素と周辺画素とを含むブロック内の各画素の画素値に対して、順序統計フィルタ処理を施すフィルタ処理ステップと、該フィルタ処理ステップにより順序統計フィルタ処理が施される前後の前記注目画素の前記画素値の差分値を算出する差分値算出ステップと、前記フィルタ処理ステップにより順序統計フィルタ処理が施された前記周辺画素の前記画素値の統計的バラツキを算出するバラツキ算出ステップと、前記差分値算出ステップにより算出された差分値を前記バラツキ算出ステップにより算出された統計的バラツキにより除算して得られた値が所定の閾値以上である場合に、前記注目画素を欠陥画素であると判定する判定ステップとを含む画像処理方法である。

本発明の他の態様は、取得された画像の注目画素と周辺画素とを含むブロック内の各画素の画素値に対して、順序統計フィルタ処理を施すフィルタ処理ステップと、該フィルタ処理ステップにより順序統計フィルタ処理が施される前後の前記注目画素の前記画素値の差分値を算出する差分値算出ステップと、前記フィルタ処理ステップにより順序統計フィルタ処理が施された前記周辺画素の前記画素値の統計的バラツキを算出するバラツキ算出ステップと、前記差分値算出ステップにより算出された差分値を前記バラツキ算出ステップにより算出された統計的バラツキにより除算して得られた値が所定の閾値以上である場合に、前記注目画素を欠陥画素であると判定する判定ステップとコンピュータに実行させる画像処理プログラムである。

本発明によれば、少ない計算量で精度よく欠陥画素を検出することができるという効果を奏する。

本発明の第１実施形態に係る画像処理装置を備える撮影装置を示すブロック図である。 図１の画像処理装置を用いて実施される順序統計フィルタ処理の対象画素を示す一例である。 図２の対象画素とは異なる画素を順序統計フィルタ処理の対象画素とした一例を示す図である。 図３の画素の画素値を順序統計フィルタ処理後、それらの大きさによってソーティングした一例を示す図である。 図２の画素の画素値を順序統計フィルタ処理後、それらの大きさによってソーティングした一例を示す図である。 図１の画像処理装置を用いた画像処理方法を示すフローチャートである。 本発明の第２実施形態に係る画像処理装置を備える撮影装置を示すブロック図である。 本発明の第３実施形態に係る画像処理装置を備える撮影装置を示すブロック図である。 図８の画像処理装置を用いた順序統計フィルタ処理前後のノイズ量を示す図である。

以下、本発明の第１実施形態に係る画像処理装置１、画像処理方法および画像処理プログラムについて、図面を参照して以下に説明する。
図１に、本実施形態に係る画像処理装置１を備えたカメラ（撮影装置）１００を示す。

カメラ１００は、被写体からの光を撮影する撮像部１１０と、本実施形態に係る画像処理装置１と、表示部（図示略）とを備えている。
撮像部１１０は、被写体からの光を集光する撮影レンズ１３０と、該撮影レンズ１３０により集光された光を撮影する撮像素子１４０と、該撮像素子１４０により取得された画像信号に基づいて画像を生成する画像生成部１５０を備えている。

本実施形態に係る画像処理装置１は、画像生成部１５０から入力されてきた画像に対し、順序統計フィルタ処理を施すフィルタ処理部２と、該フィルタ処理部２による順序統計フィルタ処理前後の注目画素の画素値の差分値を算出する差分値算出部３と、順序統計フィルタ処理後の画素値に基づいて統計的バラツキを算出するバラツキ算出部４と、該バラツキ算出部４により算出された統計的バラツキと、差分値算出部３により算出された差分値とに基づいて欠陥画素か否かを判定する判定部（欠陥判定部）５とを備えている。

フィルタ処理部２は、メディアンフィルタによって順序統計フィルタ処理を行うものである。フィルタ処理部２は、例えば、図２に示されるように、ハッチングで示される注目画素ｐｔとその周辺の実線で囲まれた８個の周辺画素ｐ１からｐ８とからなる３×３＝９画素のブロックにおいて、注目画素ｐｔおよび周辺画素ｐ１からｐ８のそれぞれに対して順序統計フィルタ処理を実施するようになっている。

さらに具体的には、フィルタ処理部２は、注目画素ｐｔが設定されると、図３に示されるように、注目画素ｐｔおよび各周辺画素ｐ１からｐ８を中心とするそれぞれ３×３＝９画素の画素値をそれらの大きさによってソーティングし、中央値（メディアン値）を各画素の画素値として出力するようになっている。例えば、ｐ１の画素に対してフィルタ処理を行う場合は、図３の太線で囲った領域のようにｐ１を中心とする３×３＝９画素の範囲の画素の画素値を用いて処理をおこなう。図４に示すように、各画素の画素値を大きさによりソーティングし、中央値を出力する。この出力された値がｐ１の順序統計フィルタ処理後の値Ｐ１となる。フィルタ処理部２においては、画像上の全ての画素について、上記処理が行われるようになっている。

順序統計フィルタ処理は、撮像素子１４０がいわゆるベイヤ配列である場合には、同色の色フィルタが配置された画素に対して実施されるようになっている。
また、ここでは、順序統計フィルタ処理において、注目画素も使用しているが使用しなくてもよい。この場合には、画素数が８個（偶数）となるので中央値の算出は、４番目と５番目の画素値の平均により行えばよい。

バラツキ算出部４は、統計的バラツキとして四分位範囲を算出するようになっている。
四分位範囲は、分布の代表値を中央値としたときの分布のバラツキのことをいう。
バラツキ算出部４においては、図５に示されるように、注目画素ｐｔおよびその周辺の８画素ｐ１からｐ８の順序統計フィルタ処理後の３×３＝９画素の画素値Ｐｔ，Ｐ１からＰ８が、それらの大きさによってソーティングされる。

バラツキ算出部４は、図５のソーティング結果の中央値Ｐ６を挟んで両側にそれぞれ４個ずつ存在する画素値群においてさらに中央値を求めるようになっている。中央値Ｐ６の両側の各画素値群のデータ数はそれぞれ４個なので、２番目と３番目の画素値Ｐ１，Ｐ４の平均値および７番目と８番目の画素値Ｐｔ，Ｐ７の平均値を中央値とするようになっている。

バラツキ算出部４は、２番目と３番目の画素値Ｐ１，Ｐ４から算出した中央値を第１四分位数Ｑ１、７番目と８番目の画素値Ｐｔ，Ｐ７から算出した中央値を第３四分位数Ｑ３とし、下式により四分位範囲ＩＱＲを算出するようになっている。
ＩＱＲ＝Ｑ３−Ｑ１

判定部５は、差分値算出部３により算出された差分値を、バラツキ算出部４により算出された統計的バラツキである四分位範囲によって除算して得られた値を所定の閾値と比較し、閾値より大きい場合には、欠陥画素であると判定するようになっている。
四分位範囲はバラツキであるため、バラツキが小さいと注目画素の周辺は平坦であり、バラツキが大きいと注目画素の周辺は平坦ではないことが分かる。四分位範囲を分母とすることにより、バラツキが小さく周辺が平坦な場所でのみ欠陥画素を判定することができる。また、注目画素が突出している場合には、差分値は大きくなるため、差分値を分子とすることにより、注目画素の突出度を容易に判定することができる。

このように構成された本実施形態に係る画像処理装置１を用いた画像処理方法について以下に説明する。
本実施形態に係る画像処理方法は、図６に示されるように、取得された画像の注目画素と周辺画素とを含むブロック内の各画素の画素値に対して、順序統計フィルタ処理を施すフィルタ処理ステップＳ１と、順序統計フィルタ処理が施される前後の注目画素の画素値の差分値を算出する差分値算出ステップＳ２と、順序統計フィルタ処理が施された周辺画素の画素値の統計的バラツキを算出するバラツキ算出ステップＳ３と、差分値を統計的バラツキにより除算して得られた値が所定の閾値以上である場合に、注目画素を欠陥画素であると判定する判定ステップＳ４とを含んでいる。

本実施形態に係る画像処理装置１および画像処理方法によれば、メディアンフィルタ処理のような順序統計フィルタ処理を施すことにより、ノイズにより突出している画素値や、複数存在している他の欠陥画素の画素値の影響を低減して、欠陥画素の検出精度を向上することができるという利点がある。

また、単純に順序統計フィルタ処理の前後の画素値の差分値を閾値によって評価するのではなく、順序統計フィルタ処理後の周辺画素から算出した統計的バラツキによって差分値を除算することにより得られた値を閾値によって評価するので、注目画素がエッジに隣接もしくはテクスチャ上に存在する場合の誤検出を抑制することができるという利点がある。

さらに、本実施形態に係る画像処理装置１および画像処理方法によれば、バラツキ算出部４が、統計的バラツキとして四分位範囲を用いているので、標準偏差と比較して計算量を大幅に低減することができ、消費電力の節約を図ることができるという利点もある。

次に、本発明の第２実施形態に係る画像処理装置１０について、図面を参照して以下に説明する。
本実施形態の説明において、上述した第１実施形態に係る画像処理装置１と構成を共通とする箇所には同一符号を付して説明を省略する。

本実施形態に係る画像処理装置１０は、撮像部１１０により取得された動画が入力される点で第１の実施形態に係る画像処理装置１と相違している。
本実施形態に係る画像処理装置１０は、図７に示されるように、判定部５により欠陥画素であると判定された場合には、その画素の座標値を登録するアドレス登録部１１と、前フレームで算出された評価値を記憶する評価値記憶部１２と、記憶されている前フレームの評価値と差分値とに基づいて評価値を算出する評価値算出部１３とを備え、アドレス登録部１１は、評価値算出部１３により算出された評価値に応じて登録を解除するようになっている。

すなわち、アドレス登録部１１により欠陥画素の座標値が登録されると、それ以降の登録画素による画像の検出は行われないようになる。
評価値算出部１３は、欠陥画素の評価値を、下式により、フレーム毎に算出するようになっている。
ｖｃ＝ｖｂ−ｖｂ＊ｎ＋ｖｎ＊ｎ

ここで、ｖｃは更新された評価値、ｖｂは前フレームの評価値、ｎは係数、ｖｎは現フレームにおける差分値である。上記式によれば、現フレームの差分値と前フレームの評価値とが重み付け合成される。
なお、初回の評価値の算出においては前フレームの評価値が存在しないので、初期値として、任意の値もしくは差分値が設定されればよい。

アドレス登録部１１は、評価値算出部１３により算出された評価値と任意の閾値とを比較し、閾値以下の場合には登録された座標の画素は欠陥画素ではないと判定して座標値の登録を解除するようになっている。この場合には評価値はゼロとなり、座標値の登録を解除しない場合には、評価値がそのまま出力されて、評価値記憶部１２に記憶されるようになっている。

このように構成された本実施形態に係る画像処理装置１０によれば、誤検出により欠陥画素であると判定されてアドレスに登録された場合であっても、複数フレームの処理の進行に従って、誤検出であることが明らかになり、アドレス登録部１１における登録が解除される。したがって、健全な画素がいつまでもアドレスに登録されたままの無駄な状態を解消することができる。また、事前に用意するアドレス登録用のメモリ量を削減することができるという利点もある。

次に、本発明の第３実施形態に係る画像処理装置２０について、図面を参照して以下に説明する。
本実施形態の説明において上述した第２の実施形態に係る画像処理装置１０と構成を共通とする箇所には同一符号を付して説明を省略する。

本実施形態に係る画像処理装置２０は、図８に示されるように、間引き部２１と、平坦判定部２２とを備える点で第２実施形態に係る画像処理装置１０と相違している。
間引き部２１は、判定部５による欠陥画素の判定がなされた場合に、数回に１回の割合でアドレス登録部１１に欠陥画素の座標値を登録するように間引く機能を有している。

具体的には、例えば、欠陥画素であると判定される都度にカウント数をインクリメントしていき、カウント数が偶数の場合にはアドレス登録部１１への登録を行わないようになっている。あるいは、カウント数が下式によって算出される数列Ｎに該当する場合に欠陥画素の画素値を登録することにしてもよい。下式は４枚置きであるが、２枚置き、３枚置き、５枚置き以上であってもよい。
Ｎ＝４＊ｘ＋１ ｘ＝０，１，２，３，…

平坦判定部２２は、順序統計フィルタから出力された順序統計フィルタ処理された周辺画素の画素値の統計的バラツキと、閾値とを比較して閾値を上回っている場合には平坦ではないと判定するようになっている。
また、評価値算出部１３は、平坦判定部２２から平坦ではないと判定された画素については、評価値を算出しない、すなわち、評価値の更新を行わないようになっている。

この場合において、閾値は、判定を行う画素の画素値毎に算出されるように、画素値毎の平坦な領域における光ショットノイズのノイズモデルを使用する。使用するノイズモデルとしては、順序統計フィルタ処理された画素値に応じたノイズ量のモデルを使用することが好ましい。ノイズ量は、画素値毎の平坦な領域における標準偏差や分散を使用することが一般的だが、この場合は画素値毎の平坦な領域における四分位範囲を使用する。

順序統計フィルタ処理の前後においては、図９に示されるように、各画素値におけるノイズ量に差が生ずるため、このノイズモデルから算出されるモデル値に下式のように任意の係数を乗算することで最終的な閾値として使用すればよい。
ＴＨ＝ＭＶ＊Ｋ
ここで、ＴＨは閾値、ＭＶはモデル値、Ｋは任意の係数である。

このように、本実施形態に係る画像処理装置２０によれば、判定部５により欠陥画素の判定がなされた全ての場合にアドレス登録部１１への登録が行われるのではなく、間引き部２１によって、登録される欠陥画素が間引かれるので、欠陥画素の誤検出数を低減することができるという利点がある。これにより、事前に用意するアドレス登録用のメモリ量を低減することができる。

また、本実施形態に係る画像処理装置２０によれば、平坦判定部２２により平坦ではないと判定された画素については欠陥画素の評価値が算出されないので、急激に変化するシーンでの欠陥画素の評価値の誤った更新を抑制することができ、これにより正検出された画素の検出状態の維持と誤検出された画素の検出状態の解除を安定的に行うことが出来るという利点がある。

なお、上記各実施形態においては、順序統計フィルタ処理としてメディアンフィルタによる処理を例示したが、これに代えて、ソート後の画素値群から指定の順位の画素値を選択するランクオーダーフィルタ処理、ソート後の画素値群の両端を指定した数だけ除外し残った値で平均値を算出するα−トリムド平均値フィルタ処理、ソート後の中央値を選択するだけでなく、中央値から一定の範囲に含まれる画素値も選択し、それらの平均値を算出する、ダブルウィンドウモディファイドトリムド平均値フィルタ処理、あるいは、欠陥画素の大きさに応じて最小値や最大値を選択する処理を採用してもよい。

また、統計的バラツキとして四分位範囲を例示したが、これに代えて、分散や標準偏差、平均絶対値差あるいは範囲を用いてもよい。範囲とは、入力画素の画素値の分布幅、すなわち、最大画素値と最小画素値との差である。

また、第２実施形態および第３実施形態においては、撮像部１１０により取得された動画が入力される場合について例示したが、これに代えて、時系列的に取得された複数枚の静止画が入力される場合に適用してもよい。

また、本実施形態においては、画像処理方法をハードウェアによって実現する構成について説明したが、コンピュータにより実行可能な画像処理プログラムによって実現することとしてもよい。この場合、コンピュータは、ＣＰＵ、ＲＡＭ等の主記憶装置、およびコンピュータ読み取り可能な記録媒体を備え、該記録媒体に上述した処理を実現させるための画像処理プログラムが記録される。そして、ＣＰＵが記録媒体に記録されている画像処理プログラムを読み出して、画像処理装置１，１０，２０と同一の処理を実現することができる。

また、本実施形態においては、画像処理装置１がカメラ１００に適用される場合について説明したが、これに代えて、顕微鏡、内視鏡等の任意の撮影装置に適用されることにしてもよい。

１，１０，２０ 画像処理装置
２ フィルタ処理部
３ 差分値算出部
４ バラツキ算出部
５ 判定部
１１ アドレス登録部
１２ 評価値記憶部
１３ 評価値算出部
２１ 間引き部
Ｓ１ フィルタ処理ステップ
Ｓ２ 差分値算出ステップ
Ｓ３ バラツキ算出ステップ
Ｓ４ 判定ステップ

Claims (8)

1. 取得された画像の注目画素と周辺画素とを含むブロック内の各画素の画素値に対して、順序統計フィルタ処理を施すフィルタ処理部と、
   該フィルタ処理部により順序統計フィルタ処理が施される前後の前記注目画素の前記画素値の差分値を算出する差分値算出部と、
   前記フィルタ処理部により順序統計フィルタ処理が施された前記周辺画素の前記画素値の統計的バラツキを算出するバラツキ算出部と、
   前記差分値算出部により算出された差分値を前記バラツキ算出部により算出された統計的バラツキにより除算して得られた値が所定の閾値以上である場合に、前記注目画素を欠陥画素であると判定する判定部とを備える画像処理装置。
2. 前記バラツキ算出部が、統計的バラツキとして四分位範囲を算出する請求項１に記載の画像処理装置。
3. 前記判定部により前記欠陥画素であると判定された画素の座標値を登録するアドレス登録部と、
   該アドレス登録部に登録された座標値に対応する画素の前記欠陥画素の評価値を記憶する評価値記憶部と、
   前記差分値算出部により算出された差分値と、前記評価値記憶部に記憶されている過去の前記評価値とに基づいて前記欠陥画素の前記評価値を算出する評価値算出部とを備え、
   前記アドレス登録部は、前記評価値算出部により算出された前記評価値が、所定の閾値以下である場合に登録を解除する請求項１または請求項２に記載の画像処理装置。
4. 前記評価値算出部は、前記差分値と過去の前記評価値とを重み付け合成して前記評価値を算出する請求項３に記載の画像処理装置。
5. 前記判定部により前記欠陥画素であると判定された画素の前記座標値の前記アドレス登録部への登録を間引く間引き部を備える請求項３または請求項４に記載の画像処理装置。
6. 前記評価値算出部は、前記バラツキ算出部が算出した順序統計フィルタ処理が施された前記周辺画素の前記画素値の統計的バラツキが所定の閾値より大きい場合に、前記評価値の算出を行わない請求項３から請求項５のいずれかに記載の画像処理装置。
7. 取得された画像の注目画素と周辺画素とを含むブロック内の各画素の画素値に対して、順序統計フィルタ処理を施すフィルタ処理ステップと、
   該フィルタ処理ステップにより順序統計フィルタ処理が施される前後の前記注目画素の前記画素値の差分値を算出する差分値算出ステップと、
   前記フィルタ処理ステップにより順序統計フィルタ処理が施された前記周辺画素の前記画素値の統計的バラツキを算出するバラツキ算出ステップと、
   前記差分値算出ステップにより算出された差分値を前記バラツキ算出ステップにより算出された統計的バラツキにより除算して得られた値が所定の閾値以上である場合に、前記注目画素を欠陥画素であると判定する判定ステップとを含む画像処理方法。
8. 取得された画像の注目画素と周辺画素とを含むブロック内の各画素の画素値に対して、順序統計フィルタ処理を施すフィルタ処理ステップと、
   該フィルタ処理ステップにより順序統計フィルタ処理が施される前後の前記注目画素の前記画素値の差分値を算出する差分値算出ステップと、
   前記フィルタ処理ステップにより順序統計フィルタ処理が施された前記周辺画素の前記画素値の統計的バラツキを算出するバラツキ算出ステップと、
   前記差分値算出ステップにより算出された差分値を前記バラツキ算出ステップにより算出された統計的バラツキにより除算して得られた値が所定の閾値以上である場合に、前記注目画素を欠陥画素であると判定する判定ステップとコンピュータに実行させる画像処理プログラム。

Images