JP3865821B2

JP3865821B2 - ディジタル画像のシーケンスの画像内のノイズの時間的フィルタリングの方法及び装置

Info

Publication number: JP3865821B2
Application number: JP17146796A
Authority: JP
Inventors: マクラン−エベイドシェリフ
Original assignee: Koninklijke Philips Electronics NV
Current assignee: Koninklijke Philips NV
Priority date: 1995-06-30
Filing date: 1996-07-01
Publication date: 2007-01-10
Anticipated expiration: 2016-07-01
Also published as: US5909515A; FR2736181A1; DE69627755D1; DE69627755T2; JPH09147103A; EP0751482B1; EP0751482A1

Description

【０００１】
【発明の属する技術分野】
本発明は、ディジタル画像のシーケンスの画像内のノイズの時間的フィルタリングの方法と、その方法を実行する装置とに関する。
本発明は、特に、Ｘ線ビームが低い明暗度を有し、ディテールを劣化することなく画像からノイズを除去するためにフィルタリングされるべきノイズ性の低コントラスト画像を生成するシステムを用いて、Ｘ線蛍光透視モードで形成された医用画像を処理する際に使用可能である。
【０００２】
時間的フィルタリングは、一次元の時間の関数として、所謂、シーケンスの画像内に所定の場所を有する画素の明暗度の値によって形成された時間的信号を平滑化することにより行なわれる。
【０００３】
【従来の技術】
時間的フィルタリング方法は、アールイーカルマン(R.E. Kalman) による論文：“線形フィルタリング及び予測問題の新しいアプローチ(A new approach to linear filtering and prediction problems)”、“アメリカ機械学会論文誌”、基礎工学ジャーナル、シリーズ８２Ｄ、３３−４５ページ、１９６０年によって周知である。
【０００４】
カルマンフィルタは、アプリオリに作成された仮説の関数、シーケンスの中の前の画像内の同じ場所の画素の明暗度の関数、及び、カルマン利得と呼ばれる係数の関数として、シーケンスの画像の瞬時的な画素のフィルタリングされた明暗度を生成する帰納的な式によって定義される。
上記式は二つの帰納的なアルゴリズムを生成する。二つの画像の間に僅かな動きが生じると即座に、上記動きによって、平滑化されるべき時間的な信号の曲線に現れ、明暗度不連続点と呼ばれる立ち上がり又は立ち下がりエッジが発生するという問題が存在する。
【０００５】
第１のアルゴリズムの場合、カルマン利得は完全に一定であるよう選択され；これは、動きによって誘起された明暗度の不連続点に影響を与える指数関数的な縞（ストリーキング）を発生させる。かくして、明暗度信号に段差を生じさせるノイズ性の元画像内の小さい対象物、例えば、素早く移動させ得るカテーテルは、段差の側面がフィルタリングによって変形されるので、フィルタリングされた画像から消える可能性がある。上記アルゴリズムは小さい対象物を消去する。
【０００６】
第２のアルゴリズムの場合、カルマン利得は、所定の時点で所定の場所を占める画素に対し観察されたノイズ性の明暗度と、前の時点でフィルタリングされた明暗度との差の関数である。その結果として、時間的信号は、不連続点の前で平滑化されるが、しかし、明暗度の不連続点の後でフィルタリングされることがないので、残留ノイズが明暗度の不連続点のエッジの後ろに存在する。
【０００７】
【発明が解決しようとする課題】
従って、周知の時間的フィルタリング方法は、動いている小さい対象物を表わす非常にノイズの多い画像の系列に効果的に適用できないという欠点がある。
かくして、周知の時間的フィルタリング法は、例えば、カテーテルのような非常に小さい径の器具が観察される領域内で挿入又は移動されている間に医学的な手術を実時間で追跡するため実行されたＸ線蛍光透視モードで得られた画像のシーケンスに、時間的フィルタリング法が適用されたとき生じる幾つかの主要な問題を解決しない。
【０００８】
第１の問題は、蛍光透視法のＸ線ビームの非常に低い明暗度に起因して、シーケンスの画像は非常に多くのノイズを含み、屡々、ノイズのピークを含むことである。
第２の問題は、シーケンスの各画像が短い時間的ギャップによって次の画像から分離されているので、カテーテルのような小さい器具の移動のような重要な事象が、画像と別の画像の間に生じ得ることである。小さい対象物の動きが出現する画像のフィルタリング対象物を歪め、又は、消すことは許されない。
【０００９】
本発明の目的は、
− 厳密に実時間的に実行され、即ち、観察された瞬時的な画像の出現の直後に、画像遅延を伴うことなく、かつ、シーケンス内の画像レートとは無関係にフィルタリングされた瞬時的な画像を出力し、
− 時間的明暗度信号の不連続性エッジの後で、不連続性エッジを弱めることなく残留ノイズを低減し、
− ノイズピークを低減することが可能であり、
− 動いている小さい対象物を消去又は歪めることなく、
ディジタル化された画像のシーケンスの連続する画像内のノイズを低減する時間的フィルタリング方法を提供することである。
【００１０】
【課題を解決するための手段】
上記目的は、ディジタル化されたサンプルと呼ばれるノイズ性の明暗度の値を有する画素の２次元マトリックスの形で画像のシーケンスの一部を形成する現在の画像と呼ばれる画像内のノイズの時間的フィルタリング方法であって：
上記現在の画像の所定の場所の画素に対応するノイズ性のサンプルを再生するため、重み付けされたノイズ性の因果的なサンプルと呼ばれる上記シーケンス内の現在のノイズ性サンプルと前のサンプルの線形結合によって、現在のフィルタリングされた因果的なサンプルと呼ばれるフィルタリングされたサンプルを評価する段階からなり、
上記重みは、上記重みが割り当てられた該ノイズ性の因果的なサンプルと、値１を有する重みが割り当てられた上記現在のノイズ性のサンプルの間で、０から１まで変わる明暗度の連続性の確率の係数である方法によって解決される。
【００１１】
上記フィルタリング方法を実行する装置は：
− ノイズ性の現在のサンプルと呼ばれ、現時点に得られる画素のマトリックスの形で画像に所定の場所を有する画素のノイズ性のディジタル化された明暗度を供給する画像処理装置と、
− 入力が現在のノイズ性のサンプルを受け、現在の画素よりも前の画素のノイズ性又はフィルタリングされた因果的サンプルと呼ばれ、前の画像のマトリックス内の同一場所にあるディジタル化された明暗度を供給するメモリと；因果的な利得倍率を供給するメモリと；因果的な積分式を評価し、上記ノイズ性の現在のサンプルを再生するためフィルタリングされた場合に、現在のフィルタリングされた因果的なサンプルの値を出力に発生する計算手段とからなる因果的なサブ組立体と呼ばれるサブ組立体とにより構成されている。
【００１２】
上記装置は、容易に作動させることができ、画像のディテールを劣化させることなく、ノイズの効果的、厳密かつ正確な実時間の時間的フィルタリングを提供する利点が得られる。
【００１３】
【発明の実施の形態】
以下、添付図面を参照して本発明を詳細に説明する。
（１）Ｘ線装置
図１には、Ｘ線源１と、患者のための移動台２と、マイクロプロセッサからなるディジタル画像処理装置５にデータを提供する撮像管４に接続された画像増倍器３とにより構成されたディジタルＸ線透視図システムが示されている。ディジタル画像処理装置５は、複数の出力を有し、その中の一つの出力６は、Ｘ線透視画像又は明暗度画像のシーケンスを表示するモニター７に接続されている。
【００１４】
ディジタル化されたＸ線透視画像は、８ビット又は１０ビットで符号化された５１２×５１２画素又は１０２４×１０２４画素によって構成可能である。各画素は、従って、２５６又は１０２４明暗度レベルの中の一つのレベルを有する。例えば、画像の暗い領域は低い明暗度レベルを有し、明るい領域は高い明暗度レベルを有する。
【００１５】
ディジタル化された画像は蛍光透視モードで得られる。本発明は、特に、血管造影画像のフィルタリングに使用することが可能である。
いずれにしても、本発明は、ディジタル画像のシーケンスを取得するため使用される方法、或いは、画像の中に表わされた対象物の性質を考慮していないが、ノイズを除去するためシーケンス内の画像をフィルタリングすることだけに関係している。
【００１６】
（２）時間的ノイズフィルタリング方法
本発明は、ノイズ性のディジタル画像のシーケンスの各画像内で連続的にノイズの時間的フィルタリングを行う方法を提案する。上記方法は、最後の観察されたノイズ性画像内でフィルタリング演算を行う。従って、上記方法は厳密に実時間で行われる。
【００１７】
上記方法はシーケンス内の画像に記録された対象物の動きを検出することが可能である。これは、適応的な反復フィルタリングステップに基づいている。
図２を参照すると、ノイズの時間的フィルタリング方法は、最初に、所定のレートで形成されたノイズ性画像のシーケンスの取得及びディジタル化を行う。
時間的に最も離れた時点から最も近い時点までの間で、上記シーケンスは：
− ｋ₀がシーケンス内の画像の数よりも１だけ小さい整数に一致し、ｔがフィルタリングされる画像の到達の時点を表わすとき、最初の時点ｊ＝ｔ−ｋ₀から時点ｊ＝ｔ−１までの間に形成され、Ｊ_j ^Cで示され、“因果的な画像”と呼ばれる多数の過去の明暗度画像と、
− 時点ｊ＝ｔに到達し、現在の画像と呼ばれるフィルタリングされる画像Ｊ_t ^Pとにより構成される。
【００１８】
ｊが画像の時点に対応するインデックスであるとき、ノイズ性のディジタル化された各画像Ｊ_jは、画素の２次元マトリックスＡ_j（ｘ，ｙ）からなり、各画素は、マトリックス内の座標ｘ，ｙと、明暗度レベルのスケール上の明暗度レベルＩ_j（ｘ，ｙ）とにより表わされる。
所謂ノイズ性のディジタル化された画像は、図３の（Ａ）及び（Ｂ）に示されているように、時間τの関数として、画素Ａ_j（ｘ，ｙ）の種々の明暗度レベルＩ_j（ｘ，ｙ）により形成され、同図において、明暗度Ｉが縦軸にプロットされ、時間τが横軸にプロットされている。本発明による時間的なフィルタリング方法は、図４の（Ａ）及び（Ｃ）に表わされているように、フィルタリングされた時間的信号Ｐ（τ）を得るため、時間的信号Ｉ（τ）に影響を与えるノイズをフィルタリング又は平滑化することを目的としている。以下、時間的信号Ｉ（τ）を構成する明暗度を“サンプル”と呼ぶ。過去又は因果的な明暗度はＩ_j ^Cと呼ばれ、現在の明暗度はＩ_t ^Pと呼ばれる。
【００１９】
本発明による方法は、画像シーケンスに発生した主な動きが、例えば、並進補正法のような従来技術によって周知の動き補償法によって既に補償されている場合に関係している。しかし、大量の計算時間を必要とする動き評価法は避けるべきである。本発明は厳密に実時間的な方法の提供を目的としているので、このようなタイプの方法は避ける必要がある。
【００２０】
たとえ、動き補償法が画像のシーケンスに適用されるとしても、動き補償法は完璧ではないので、全ての動きが現実的に本当に補償されることが保証された訳ではないことに注意すべきである。
この事実を考慮に入れるため、現在の方法は、画像シーケンス内の対象物の小さい残留又は補償不可能な局部的な動きを取り扱うことによって、時間的信号Ｉ（τ）のノイズ上のフィルタリング演算を実行する。
【００２１】
上記方法は、現時点ｔで得られたノイズ性サンプルＩ_t ^P上でフィルタリング演算を実行し、現時点ｔのノイズ性サンプルＩ_t ^Pと、時点ｔ−ｋ０からｔ−１の間に観察されたノイズ性因果サンプルＩ_t-k0 ^C，．．．，Ｉ_t-3 ^C，Ｉ_t-2 ^C，Ｉ_t-1 ^Cと呼ばれる前のサンプルとを利用することにより、フィルタリングされたサンプルＰ_t ^Cを生成する。
【００２２】
時間的フィルタリングは、好ましくは、２次元マトリックスの種々のｘ，ｙ座標の各画素Ａ_j（ｘ，ｙ）に対し別個に行なわれる。
図３の（Ａ）には、例えば、小さい局部的な動きに起因したＩ_t-3 ^CとＩ_t-4 ^Cの間の明暗度の不連続点Ｄからなるノイズ性時間的信号Ｉ（τ）が表わされている。上記小さい局部的な動きは、時点ｔ−４と時点ｔ−３の間で時間的信号Ｉ（τ）の立ち上がりエッジと、時点ｔ−３とｔの間の大きい信号“安定状態”と、時点ｔ−４より前の時点の“抑圧状態”として表わされる。
【００２３】
本発明による時間的フィルタリングは、Ｄのような明暗度の不連続点に適応する。かくして、図４の（Ａ）及び（Ｃ）は、本発明による方法の種々の変形例によって平滑化された時間的信号Ｐ（τ）を表わしている。
上記時間的信号は、ノイズに起因して、シーケンスの画像内で同一の座標（ｘ，ｙ）を有する画素の明暗度レベルが絶えず変化するため、時点ｔ−８とｔ−４の間に小さい鋸波を更に有し、これにより、降雪のような斑な局面が現れる。フィルタリングは、時間的信号Ｉ（τ）を所定の画素に関して一定の平均値に平滑化することにより、かかる局面を抑圧する。
【００２４】
更に、図３には、例えば、時点ｔ−４とｔ−２の間、即ち、時点ｔ−３に現れたノイズのピークＤ’を有するノイズ性の時間的信号Ｉ（τ）が表わされている。かかるノイズのピークＤ’は、実際の空間的な動きの現象に対応していない点で、図３の（Ａ）のＤのような明暗度の不連続点とは区別される。かかるノイズのピークは、本発明によって提案された方法を用いてフィルタリングすることができる。
【００２５】
現在のノイズ性サンプルＩ_t ^Pに対応しているフィルタリングされたサンプルＰ_t ^Cは、以下の式：
【００２６】
【数１】

【００２７】
のように表現される。フィルタリングされたサンプルＰ_t ^C、又は、時点ｔでフィルタリングされた時間的信号は、因果的な時点ｊ＝ｔ−ｋ₀からｊ＝ｔ−１で観察されたノイズ性サンプル及び時点ｔで観察されたノイズ性サンプルの線形結合である。
上記フィルタリングされた信号Ｐ_t ^Cの式（１）において、重みｂ_j ^Cが因果的なサンプルＩ_t-1 ^C、Ｉ_t-2 ^C等に割り当てられる。
【００２８】
上記（１）式の分母は正規化係数であり、その正規化係数の存在は、フィルタリングされた信号Ｐ_t ^Cの平均値がノイズ性信号Ｉ_t ^Pの平均値と一致することを保証するため、種々のサンプルに付けられた重みの合計は１に一致すべきであるという事実に基づいている。
時間的フィルタリング方法は、因果的なサンプルの現在のサンプルとの線形結合に基づいている。上記方法において、因果的なサンプルに付けられた重みｂ_j ^Cは特定の形式を有し、現在のサンプルに付けられた重みは１である。
【００２９】
所定の因果的なサンプルに付けられた重みは、所定の因果的なサンプルと現在のサンプルの間の明暗度の連続性の確率である。
一般的に言うと、因果的なサンプルＩ_t-k ^Cに割り当てられるべき因果的な重みｂ_t-k ^Cの式は、次式（２）：
【００３０】
【数２】

【００３１】
によって与えられ、式中、α_j ^Cは、Ｉ_t-k ^CからＩ_t ^Pまでの連続した因果的なサンプルの間の明暗度の連続性の確率である。上記関係式（２）によれば、因果的なサンプルＩ_t-k ^Cに関係する重みｂ_t-k ^Cは、因果的な時点ｊ＝ｔ−ｋ＋１から現在の時点ｔまでの明暗度の連続性の全確率α_j ^Cの積に一致する。
従って、式（１）に従って所定の因果的なサンプルに付けられた重みの決定は、現在のサンプルから大きく離れ過ぎることがなく、同一の対象物に関係することを示しているならば、上記所定のサンプルが確保され、線形結合に考慮されるという仮説に基づいている。
【００３２】
上記仮説に従って、例えば、因果的なサンプルＩ_t-1 ^Cに関係する因果的な重みｂ_t-1 ^Cは、ｂ_t-1 ^Cが、現時点のサンプルＩ_t ^Pと前の時点のサンプルＩ_t-1 ^Cの差の絶対値の関数であることを記述することにより定式化される。
好ましくは、因果的な重みｂ_t-1 ^Cは、以下の差：
Δ^C＝｜Ｉ_t ^P−Ｐ_t-1 ^C｜（６^C）
の関数であり、式中、Ｐ_t-1 ^Cは、前の時点ｔ−１で既にフィルタリングされたサンプルであるので、フィルタリングの結果としてＩ_t-1 ^Cよりノイズが少ないと考えられる。サンプルＩ_t ^PとＰ_t-1 ^Cの差が小さい場合、１に近い“大きい”値が対応する重みｂ_t-1 ^Cに割り当てられる。上記の差が大きい場合、ゼロに近い値が重みｂ_t-1 ^Cに割り当てられる。その場合、サンプルＩ_t-1 ^Cは、実際上、考慮されない。
【００３３】
次いで、上記例の場合、第２の因果的な重みｂ_t-2 ^Cは、時点ｔのサンプルと、時点ｔ−１のサンプルの間の差の関数だけではなく、時点ｔ−１のサンプルと時点ｔ−２のサンプルの間の関数である。かくして、現在のサンプルＩ_t ^Pに関し、ノイズに起因したサンプル以外の修正されていないサンプルだけが考慮されるならば、即ち、サンプルの間の差は小さくなければならないという条件が課されるならば、重みｂ_t-2 ^Cは１に近い大きい値を有する。これより、時間的信号に現れる明暗度差の関数の積として因果的な重みｂ_j ^Cの式が得られ、即ち、フィルタリングされるべき現在のサンプルに先行するサンプルの間で明暗度の連続性の積として上記重みの式が得られる。
【００３４】
従って、サンプルＩ_t-1 ^Cに関係する因果的な重みの式は：
ｂ_t-1 ^C＝α_t ^C （２ａ）
であり、式中、α_t ^Cは、現在のサンプルＩ_t ^PとフィルタリングされたサンプルＰ_t-1 ^Cの間の明暗度の連続性の確率である。時点ｔのサンプルと、前の時点ｔ−１のサンプルとの間の明暗度の連続性は、高い連続性の確率α_t ^Cとして表わされる。サンプルの間の明暗度の不連続点の確率は、０に近いα_t ^Cとして表わされる。
【００３５】
サンプルＩ_t-2 ^Cに関係する因果的な重みの式は：
ｂ_t-2 ^C＝α_t ^C×α_t-1 ^C （２ｂ）
であり、式中、α_t-1 ^Cは、時点ｔ−１のサンプルと、時点ｔ−２のサンプルの間の明暗度の連続性の確率である。
時点ｔと時点ｔ−２の間のサンプルの明暗度の連続性は、時点ｔとｔ−１の間の明暗度の連続性と、時点ｔ−１とｔ−２の間の明暗度の連続性とを意味し、時点ｔとｔ−１の間の明暗度の連続性の確率が高く、同様に時点ｔ−１とｔ−２の間の明暗度の連続性の確率が高い確率（２ｂ）の積として表わされる。
【００３６】
サンプルＩ_t-3 ^C等に関係する因果的な重みの式は：
ｂ_t-3 ^C＝α_t-2 ^C×α_t-1 ^C×α_t ^C （２ｃ）
等によって与えられる。
従って、図３の（Ａ）の例の場合、連続性の確率α_t-3 ^Cは、サンプルＩ_t-4 ^CとＩ_t-3 ^Cの間の明暗度の不連続点Ｄのため、非常にゼロに近い。確率α_t-3 ^Cを含む確率の積として得られる重みｂ_t-4 ^Cは、同様にゼロに近づく傾向がある。重みｂ_t-5 ^Cは、その式が上記ゼロに近い確率を含む積であるため、同様にゼロに近づく傾向がある。重みｂ_j ^Cを与える確率の積の中にゼロに近い確率が現れると直ぐに、上記確率の存在によって、上記積と、ゼロに近い上記確率が現れた時点より前の時点に対応する全ての積とが打ち消される。その結果として、上記前の時点に対応する因果的なサンプルの重みｂ_j ^Cは同様にゼロに近くなる。図３の（Ａ）の具体的な例の場合、因果的な重みは：
ｂ_t-1 ^C≒１；ｂ_t-2 ^C≒１；ｂ_t-3 ^C≒１；ｂ_t-4 ^C≒０；ｂ_t-5 ^C≒０
によって表わされる。
【００３７】
因果的な係数ｂ_j ^Cの具体的な式（２）の結果として、計算が非常に複雑な式（１）は、所望の時間的フィルタリングを実行する三つの簡単な関係式からなる一つのグループに変換される。
因果的な重みが特定の形式（２）を有する線形性の式（１）により、以下の帰納的な式：
Ｐ_t ^C＝Ｐ_t-1 ^C＋Ｋ_t ^C×（Ｉ_t ^P−Ｐ_t-1 ^C）（３^C）
が定式化される。
【００３８】
積分式と呼ばれる上記式（３^C）は、時点ｔ−１でフィルタリングされた因果的なサンプルＰ_t-1 ^Cを含み、この因果的なサンプルＰ_t-1 ^Cは、時点ｔ−１でフィルタリングされたサンプルＰ_t-1 ^Cと現時点ｔで観察されたノイズ性サンプルＩ_t ^Pの差によって形成された因果的な更新項と呼ばれる項によって修正されている。上記因果の更新項は、０と１の間を変わる因果の利得Ｋ_t ^Cと呼ばれる項で乗算される。上記式（３^C）は、所望のフィルタリングされたサンプルであるフィルタリングされた値Ｐ_t ^Cを生成する。
【００３９】
上記因果的な利得Ｋ_t ^Cは、カルマンの利得ではないが、以下の反復式（４^C）：
【００４０】
【数３】

【００４１】
によって与えられる全く別の利得倍率である。
式（１）及び（２）と、式（３^C）及び（４^C）の間の推移に対し、因果的な利得は、因果的な重みの逆数：
【００４２】
【数４】

【００４３】
であり、更に、上記の如く：
α_t ^C＝ｂ_t-1 ^C （２ａ）
のように表わされ、式中、α_t ^Cは、現在のサンプルＩ_t ^Pよりも前の因果的なサンプルＩ_t-1 ^Cの重みを表わし、その確率を以下連続性の因果係数と呼ぶ。
式（１）＋（２）と、式（３^C）の間の推移は、以下の因果式：
積分式Ｐ_t ^C＝Ｐ_t-1 ^C＋Ｋ_t ^C×（Ｉ_t ^P−Ｐ_t-1 ^C）（３^C）
と、因果的な利得倍率である
【００４４】
【数５】

【００４５】
と、因果的な連続性の係数（５^C）である
【００４６】
【数６】

【００４７】
からなるシステムに基づき、ここで、Ｓ_t ^Cは正規化係数を表わし、Ｆ^Cは以下で定義されている関数である。
時間的にｔ−ｋ₀で始まる式（３^C）の反復により、上記定式化されたように式（１）が厳密に得られる。
積分式の定式化（３^C）は、かくして、連続したサンプルの間の明暗度の連続性の関数であり、現在のサンプルに至るまで重みを計算された連続するサンプルの間の連続性の確率の積（２）によって定められた特定の重みと関係したサンプルの線形結合（１）に直接対応する。
【００４８】
・因果的な不連続性の係数α_t ^Cの定式化
因果的な不連続性の係数α_t ^Cは、差：
Δ^C＝｜Ｉ_t ^P−Ｐ_t-1 ^C｜（６^C）
の関数として定義される。
上記差Δ^Cが、図３の（Ａ）のＤのような不連続点に関係しているか、或いは、ノイズだけに関係しているかどうかを判定することにより利点が得られる。
【００４９】
ノイズの関与を判定するため、上記差Δ^Cは、上記不連続性Δ^Cの各サンプルに関するノイズの分散を考慮した係数Ｓ_t ^Cによって正規化される。σ_Bと呼ばれ、各サンプルに対する明暗度レベルで測定されたノイズの標準偏差は、あらゆる周知の方法を用いて評価し、或いは、アプリオリに評価することが可能である。
【００５０】
ノイズの分散σ_B ²は、画像内で想定された領域のノイズの平方の平均偏差であり、ノイズに起因した明暗度の偏差は、平均局部明暗度に関して測定されている点に注意する必要がある。
標準ノイズ偏差の判定は、図１に示したシステムにおいて、信号のディジタル化の前に、低明暗度レベル及び高明暗度レベルの圧縮の規則に従って画像の変換が導入されることによって複雑化するので、ノイズは、非常に暗い領域又は非常に明るい領域の場合、中間の明暗度の領域よりも見えにくくなるということを考慮する。標準ノイズ偏差σ_Bは、全ての画素の全ての時点における明暗度の関数である。
【００５１】
従って、標準ノイズ偏差の判定は、対応する領域の平均明暗度の関数として、画像の各画素上で局部的に実行する必要がある。かくして、画像自体又は信号に起因したノイズの分離が行なわれるように、ノイズの分散及び観察された信号の平均値が局部的に計算されるので、ノイズは空間的に表わされている。
図５の（Ａ）を参照して一例として提案された方法に従って、画像Ｊ_j内で、典型的に、各画素Ａ_j（ｘ，ｙ）の周囲の５×５又は３×３画素の区域からなる小さいウィンドウＷ_jが定義される。上記各小さいウィンドウＷ_j（ｘ，ｙ）において、ノイズの分散及び観察された平均明暗度が判定される。
【００５２】
このため、ｎが小さいウィンドウＷ_j（ｘ，ｙ）の画素数を表わす場合に、画素Ａ_j（ｘ，ｙ）の周辺に中心が置かれた小さいウィンドウの各位置に対し、以下の如く：
− インデックスｋ及びｌが、画像Ｊ_j内の矩形状の座標Ｘ及びＹに関してウィンドウ内にあると考えられる画素の座標を表わす場合に、小さいウィンドウＷ_j（ｘ，ｙ）の各画素の明暗度Ｉ_j（ｋ，ｌ）と、
− この位置で、関係式：
【００５３】
【数７】

【００５４】
によって与えられるウィンドウ内の平均明暗度ｍ_W（ｘ，ｙ）と、
− 上記ウィンドウ内で、総和の記号Σの後ろでウィンドウの各画素の明暗度と平均値ｍ_W（ｘ，ｙ）の差の平方が与えられ、（ウィンドウ内の画素数ｎ）分の（上記総和）である平均値が与えられる以下の関係式：
【００５５】
【数８】

【００５６】
によって与えられるノイズの分散σ_BW（ｘ，ｙ）の平方とが定められる。
図５の（Ｂ）を参照すると、次いで、２次元のヒストグラムが形成され、同図において、例えば、１０グレイレベル×１０グレイレベルの所定の精度で、明暗度の平均値ｍ_W（ｘ，ｙ）が水平軸上にプロットされ、一方、分散の平均値が垂直軸上にプロットされている。水平軸上にプロットされた平均明暗度ｍ_W（ｘ，ｙ）の各間隔には、分散の頻度、又は、所定の分散が上記平均明暗度の間隔に生じた回数の平均が記載されている。説明の簡単化のため図５の（Ｂ）には、一例として１０と２０の間の平均値ｍ_W（ｘ，ｙ）の値が以下のように示されている：
値１を有する分散の頻度Ｎ₀₁は０；
値４を有する分散の頻度Ｎ₁₁は１０；
値９を有する分散の頻度Ｎ₂₁は２０；
値１６を有する分散の頻度Ｎ₃₁は３５；
値２５を有する分散の頻度Ｎ₄₁は２５；以下同様である。
【００５７】
提案された方法は、分散σ² _BW（ｘ，ｙ）の関数として垂直軸上にプロットされた分散の出現頻度Ｎの曲線が、想定された平均明暗度ｍ_W（ｘ，ｙ）のあらゆる間隔に対し描かれている図５の（Ｃ）を参照して続けられる。分散の関数として、分散の可能な値を与えながら最大値を通過して延在する頻度Ｎの曲線は、平均明暗度の各間隔に対し定められる必要がある。
【００５８】
かくして、ノイズ分散が、各画素Ａ_j（ｘ，ｙ）に属する平均明暗度の値の関数として得られる。この方法によって、画像自体から抽出されるべきノイズに関する情報が得られる。
次いで、標準ノイズ偏差σ_Bが、ノイズ分散の平方根：
【００５９】
【数９】

【００６０】
の計算により求められる。
差Δ^Cに加わるノイズ寄与の分散は、安定状態の場合、次式：
【００６１】
【数１０】

【００６２】
のように定式化され、この値は、Δ^Cの正規化の近似として使用される。差Δ^Cの正規化係数は、以下の関係式：
【００６３】
【数１１】

【００６４】
によって与えられる。
式（５^C）において、連続性の係数α_t ^Cは、差Δ^Cがノイズによる影響を受けるということを考慮するため、関係式（６^C）によって与えられ、標準ノイズ偏差σ_Bの近似Ｓ_t ^Cによって正規化された差Δ^Cの関数である。このため、係数α_t ^Cの計算は、時間的信号の曲線Ｉ（τ）の各点に対し、平均値に関して標準ノイズ偏差σ_Bを判定することからなる。
【００６５】
本発明に従って提案され、関係式（４^C）及び（５^C）に従う因果的な積分関数（３^C）の計算に基づくフィルタリング方法の利点は、特に、表現α_t ^Cが提案されたように正規化される点にある。この正規化は、現在のサンプルＩ_t ^Pの分散だけではなく、フィルタリングされたサンプルＰ_t-1 ^Cのノイズの分散を考慮するので、利得倍率を与える式（５^C）によって、フィルタリングされた各サンプルに対する残留ノイズの評価が得られる。
【００６６】
帰納的な積分式（３^C）は、式（４^C）によって与えられた利得倍率Ｋ_t ^Cだけに依存しているので、容易に計算することができる。利得倍率Ｋ_t ^C自体は、過去に既に計算された利得倍率Ｋ_t-1 ^Cと、連続性係数α_t ^Cとだけに依存している。連続性係数α_t ^Cは、既に計算された利得倍率Ｋ_t-1 ^Cに依存しているので、連続性係数α_t ^Cを最初に計算する必要があり、その後、連続性係数α_t ^Cを用いて利得式Ｋ_t ^Cが簡単に計算され、フィルタリングされたサンプルＰ_t ^Cが関係式（３^C）を用いて容易に計算される。
【００６７】
・関数Ｆ^Cの定義
α_t ^Cを与える式（５^C）の計算のため、関数Ｆ^Cが以下のように定義される。
関数Ｆ^Cの引数は、以下のＺ^C：
【００６８】
【数１２】

【００６９】
によって表わされる。
関数Ｆ^C（Ｚ^C）は、例えば、図６の（Ａ）乃至（Ｄ）において、縦軸に関数の値をプロットし、横軸にＺ^Cをプロットすることにより、描かれている。
最初に、関数Ｆ^C（Ｚ^C）は、０≦Ｚ^C≦１のとき、一定値Ｆ^Cｍａｘを取り；次いで、関数Ｆ^C（Ｚ^C）は、１＜Ｚ^Cの場合に、値Ｆ^Cｍｉｎまで減少する。値Ｆ^Cｍａｘは、Ｆ^Cｍａｘ≦１を満たすような値である。値Ｆ^Cｍｉｎは、０≦Ｆ^Cｍｉｎを満たすような値である。このような関数Ｆ^C（Ｚ^C）の例を以下に説明する。
【００７０】
・例１：図６の（Ａ）に示された例
関数Ｆ^C（Ｚ^C）は一定であり、０≦Ｚ^C≦１のとき、値Ｆ^Cｍａｘ＝１に一致し；次いで、関数Ｆ^C（Ｚ^C）は、１≦Ｚ^Cの場合に、値Ｆ^Cｍｉｎ＝０まで減少するよう選択される。関数Ｆ^C（Ｚ^C）は、好ましくは、上記領域で、正規化倍率によって加えられるようなガウシアン形状を有するので：
【００７１】
【数１３】

【００７２】
のように表わされる。
・例２：図６の（Ｂ）に示された例
関数Ｆ^C（Ｚ^C）は一定であり、０≦Ｚ^C＜１のとき、１よりも小さい値Ｆ^Cｍａｘ、例えば、Ｆ^Cｍａｘ＝０．８５に一致し；次いで、関数Ｆ^C（Ｚ^C）は、１≦Ｚ^Cの場合に減少する。関数Ｆ^C（Ｚ^C）は、例えば、０．１０よりも減少することがない。関数Ｆ^C（Ｚ^C）は、この例の場合、シフトしたガウシアン曲線である。
【００７３】
・例３：図６の（Ｃ）に示された例
関数Ｆ^C（Ｚ^C）は一定であり、０≦Ｚ^C≦１のとき、１よりも小さいが１に近い値Ｆ^Cｍａｘ、従って、Ｆ^Cｍａｘ＝０．８５に一致する。次いで、関数Ｆ^C（Ｚ^C）は、１≦Ｚ^Cの場合に減少する。関数Ｆ^C（Ｚ^C）は、０に近い値Ｆ^Cｍｉｎ、例えば、Ｆ^Cｍｉｎ＝０．１０まで直線的に減少し；例えば、Ｚ^C＝２で、その値に達する。次いで、関数Ｆ^C（Ｚ^C）は、２＜Ｚ^Cの場合に、直線的に一定であり、Ｆ^Cｍｉｎ＝０．１０に一致する。
【００７４】
・例４：図６の（Ｄ）に示された例
関数Ｆ^C（Ｚ^C）は、例１の関数又は例２の関数を近似する折れ線関数である。
関数Ｆ^C（Ｚ^C）は、当業者に周知の定義を有する他の適当な形式から選択してもよい。
【００７５】
関数Ｆ^C（Ｚ^C）は、フィルタリング出力を制御する。Ｆ^Cｍａｘ＝１のとき、フィルタリング出力は最大になる。Ｆ^Cｍａｘが１よりも小さい場合、フィルタリング出力は僅かに制限される。Ｆ^Cｍｉｎが０よりも大きい場合、フィルタリング出力は決して０にならない。
・例１の関数Ｆ^Cによる本発明の実施例
図４の（Ａ）及び（Ｃ）は、共に、図３の（Ａ）に示された動きに起因して、少なくとも一つの明暗度の不連続点Ｄを含むノイズ性の時間的信号Ｉ（τ）に基づいて得られたフィルタリングされた時間的信号Ｐ（τ）を表わしている。
【００７６】
上記フィルタリング演算は、関数Ｆ^C（Ｚ^C）が図６の（Ａ）に示された例１の関数と類似した関数であるという状況で、フィルタリングシステム（３^C）、（４^C）及び（５^C）を用いて実行される。
図４の（Ａ）を参照すると、因果的なフィルタリングは、不連続点Ｄの左側で、ｔ−４からｔ−ｋ₀の方向に時間的信号の平滑化を行い；次いで、不連続点ＤのためサンプルＩ_t-4 ^C、Ｉ_t-5 ^C等によって形成された過去が記憶されていないので、因果的なフィルタリングは、時点ｔ−３でノイズ性サンプルＩ_t-3 ^Cのフィルタリングを実行できない。因果的なフィルタリングは、時間的信号の曲線のＩ_t-4 ^Cと、そのサンプルの左側で非常に効果的であり、その後、信号はｔ−３のノイズ性信号Ｉ_t-3 ^Cのままである。その理由は、因果的なフィルタリングが、時点ｔ−３で過去のデータによって修正されることなく観察されたノイズ性明暗度だけを生成するからである。次いで、因果的なフィルタリングは、時点ｔ−２において時間的信号の次の点から始めて、サンプルＩ_t-3 ^Cの右側で非常に良好に平滑化された信号をもう一度生成し、時点ｔ−３とｔ−２の間に明暗度の不連続点が存在しないので、時点ｔ−３のサンプルが考慮に入れられる。
【００７７】
従って、図４の（Ａ）に示されたような因果的なフィルタリングの後、ｔ−３とｔ−４の間の明暗度の差が大きいとき、時点ｔ−４の信号に割り当てられた係数ｂ_t-4 ^Cはゼロであるため、Ｉ_t-3 ^Cの最初のノイズ歯状突起が時点ｔ−３で残る。
一方、図３の（Ｂ）を参照すると、ノイズ性時間的信号Ｉ（τ）は、例えば、時点ｔ−３で出現するノイズピークＤ’からなる。因果的なフィルタリングは、図３の（Ｂ）のノイズ性時間的信号によって誘起された不連続点Ｄ’を、図３の（Ａ）のノイズ性時間的信号において動きによって誘起された不連続点Ｄと区別できない。従って、関数Ｆ^C（Ｚ^C）が図６の（Ａ）によって示された例１の関数に選択されたならば、反−因果的なフィルタリングは、図３の（Ｂ）に示されたようなノイズピークＤ’の適当なフィルタリングを作動可能にしない。
【００７８】
・例２の関数Ｆ^Cによる本発明の実施例
本発明によれば、上記例において、時間的信号の明暗度の不連続点Ｄの右側で発生した残留ノイズ性明暗度のピークをできる限り低減し、一方、上記時間的信号の平滑化をできる限り完璧に維持するため、フィルタリング出力は、Ｆ^Cｍｉｎ及びＦ^Cｍａｘの値を変えることにより変更可能である。
【００７９】
従って、本発明によれば、過去は、“ある程度まで”しか忘れられないので、因果的なフィルタリングは、Ｉ_t-3 ^Cで完全にはゼロになり得ない。
図３の（Ａ）の信号のフィルタリング結果を表わし、図６の（Ｂ）に示されたようにＦ^Cｍａｘ＝０．８５及びＦ^Cｍｉｎ＝０．１０の場合に関数Ｆ^C（Ｚ^C）を利用する図４の（Ｃ）を参照すると、上記例の場合、関数Ｆ^C（Ｚ^C）は、少なくとも０．１０に一致するため決してゼロにはならないので、Ｉ_t-3 ^Cの第１のノイズ歯状突起が、不連続点Ｄの後で低減されることが分かる。従って、時点ｔ−４のデータが常に考慮され、明暗度の不連続点Ｄの後、時点ｔ−３の信号のフィルタリングをある程度まで行なうことができる。
【００８０】
更に、上記の通り、図３の（Ｂ）に示されたようにノイズピークＤ’のフィルタリングのため値Ｆ^Cｍｉｎ及びＦ^Cｍａｘを変えることにより、フィルタリング出力を意のままに変更可能である。
かくして、図３の（Ｂ）の信号のフィルタリングの結果を表わし、図６の（Ｂ）に示されたようにＦ^Cｍａｘ＝０．８５及びＦ^Cｍｉｎ＝０．１０の場合に関数Ｆ^C（Ｚ^C）を利用する図４の（Ｂ）を想定することにより、時点ｔ−４のデータが常に考慮され、かつ、ノイズピークＤ’の信号のフィルタリングが時点ｔ−３で少なくともある程度まで行なえるため、関数Ｆ^C（Ｚ^C）が決してゼロではないので、ノイズピークＤ’が低減される。
【００８１】
瞬時的な画像のフィルタリングに利用可能な前の画像のデータを有することが必要とされるので、第１のフィルタリングされた画像は、実際的にシーケンスに得られた第２の画像であることに注意する必要がある。シーケンスの第１の画像のデータは、その実際の表示のためにフィルタリングされることなく、第２の画像として得られた画像のフィルタリングを実行するためフィルタリングされることはない。実際のフィルタリングは第２の画像から始まる。
【００８２】
（３）本発明の方法を実行する装置
図７には、上記フィルタリング方法を実行する装置が機能ブロックの形で表わされている。
時点ｔの現在の画像に対し、位置ｘ，ｙを有する画素のノイズ性明暗度Ｉ_t ^Pが装置の入力１００に供給される。
【００８３】
ノイズ性明暗度Ｉ_t ^Pは、二つの経路を辿る。
一方の経路において、ノイズ性明暗度Ｉ_t ^Pが、ノイズ性信号Ｉ_t ^Pの標準ノイズ偏差σ_Bを提供するルックアップテーブル１１内で表にされた関数σ_Bに与えられる。標準ノイズ偏差σ_Bは、テーブル１１の出力で得られる。利得Ｋはメモリ１７に記憶される。メモリ１７は、所定の位置ｘ，ｙを有する画素に対する前の時点の利得Ｋ_t-1 ^Cを供給する出力と、現在の時点ｔに対し計算された利得Ｋ_t ^Cを受ける入力とからなる。
【００８４】
上記装置は、ルックアップテーブル１１からの値σ_Bと、メモリ１７からの利得Ｋ_t-1 ^Cの値とが入力されると即座に、値：
【００８５】
【数１４】

【００８６】
を直接的に供給するため表にされた値を含むルックアップテーブル１３を更に有する。値Ｓ_t ^Cは、連続性係数α_t ^Cの計算を行なう式（５^C）の分母を構成する。分母の値Ｓ_t ^Cは、ルックアップテーブル１３の出力で得られ、ルックアップテーブル１４に入力される。
上記装置は、処理されている画素に対し最後にフィルタリングされた明暗度を記憶するメモリ１０を更に有する。
【００８７】
フィルタリングされた明暗度Ｐ_t ^Cが計算されるべき現時点ｔに、メモリ１０は、前の時点に対応するフィルタリングされた明暗度Ｐ_t-1 ^Cを既に含んでいるので、明暗度Ｐ_t-1 ^Cをその出力に発生する。
前の時点にフィルタリングされた明暗度Ｐ_t-1 ^Cは、時点ｔで観察されたノイズ性明暗度Ｉ_t ^Pと同時に加算器１２に供給される。加算器１２は、Ｉ_t ^Pと呼ばれる時点ｔのノイズ性明暗度と、Ｐ_t-1 ^Cと呼ばれる前の時点ｔ−１にフィルタリングされた別の明暗度とからなる二つの明暗度の間の差を出力する。
【００８８】
ルックアップテーブル１４は、加算器１２によって与えられた差の絶対値をとり、更に、分母Ｓ_t ^Cを受ける。ルックアップテーブル１４は、次に、関数Ｆ^C（Ｚ^C）の引数Ｚ^C：
【００８９】
【数１５】

【００９０】
を出力する。
ルックアップテーブル１５は、表にされた形で選択された関数Ｆ^C（Ｚ^C）を含んでいる。ルックアップテーブル１５の入力は、引数Ｚ^Cを受け、ルックアップテーブル１５は、連続性係数α_t ^C：
【００９１】
【数１６】

【００９２】
を出力する。
ルックアップテーブル１６は、メモリ１７から得られた前の時点の利得と、ルックアップテーブル１５から得られた連続性係数α_t ^Cとに基づいて、時点ｔの利得を供給する表にされた式を含む。上記利得：
【００９３】
【数１７】

【００９４】
は、ルックアップテーブル１６を介して与えられ、一方で、前の時点の利得Ｋ_t-1 ^Cの値を置き換えるためメモリ１７に供給され、他方で、乗算器１９に供給される。
次いで、積分関数（３^C）が計算される。
加算器１８は、最初に、時点ｔで観察され、第２のチャンネルを介して入力１００から得られたＩ_t ^Pと呼ばれるノイズ性信号と、前の時点ｔ−１にフィルタリングされ、メモリ１０の出力から発生したフィルタリングされた信号Ｐ_t-1 ^Cとを受け、差Ｉ_t ^P−Ｐ_t-1 ^Cを出力する。
【００９５】
この差は、ルックアップテーブル１６からの利得Ｋ_t ^Cと同時に加算器１８から乗算器１９に供給される。乗算器１９は積：
Ｋ_t ^C×（Ｉ_t ^P−Ｐ_t-1 ^C）
を形成し、その結果を出力する。
加算器２０は、一方で、乗算器１９からの積を受け、他方で、前の時点に既にフィルタリングされ、メモリ１０から発生されたＰ_t-1 ^Cと呼ばれる信号を受ける。加算器２０は、時点ｔにフィルタリングされた明暗度の値をフィルタの汎用出力２００に供給する。メモリ１０において、以下の積分式：
Ｐ_t ^C＝Ｐ_t-1 ^C＋Ｋ_t ^C×（Ｉ_t ^P−Ｐ_t-1 ^C）（３^C）
に従って、メモリ２１の入力１２０を介して入力されたフィルタリングされた値Ｐ_t ^Cが、時点ｔ−１にフィルタリングされた明暗度を置換する。
【００９６】
上記装置は、時間的フィルタリング方法を実行するため、フィルタリングされたサンプルＰ_t-1 ^C及び前の時点の利得Ｋ_t-1 ^Cの記憶だけを必要とし、更に、現時点のノイズ性明暗度Ｉ_t ^Pの値しか要求しない点に注意が必要である。入力１００で得られた明暗度を記憶することは必ずしも必要ではない。
現在の画像の種々の画素のフィルタリング後、この画像は、例えば、図１に示された表示装置のモニター７に表示される。かくして、シーケンスの種々の画像を夫々の処理の後に表示させることが可能である。
【図面の簡単な説明】
【図１】Ｘ線画像化装置の構成図である。
【図２】ディジタル画像のシーケンスを表わす図である。
【図３】（Ａ）は動きに起因する不連続点を有する座標ｘ，ｙにある画素Ａに対応するノイズ性の時間的な信号を表わし、（Ｂ）はノイズのピークを有する別のノイズ性の時間的な信号を表わす図である。
【図４】（Ａ）は図３の（Ａ）のノイズ性信号に対応するフィルタリングされた時間的な信号を表わし、（Ｂ）は図３の（Ｂ）のノイズ性信号に対応するフィルタリングされた時間的な信号を表わし、（Ｃ）は本発明の他の実施例において図３の（Ａ）のノイズ性の時間的な信号に対応するフィルタリングされた時間的な信号を表わす図である。
【図５】（Ａ）乃至（Ｃ）は平均ノイズｍ_Bに関するノイズの標準偏差σ_Bの決定法の説明図である。
【図６】（Ａ）乃至（Ｄ）は関数Ｆ^C（Ｚ^C）の種々の実現可能な例を表わす図である。
【図７】時間的なフィルタリング方法を実行する装置を機能的なブロックの形で表わした図である。
【符号の説明】
１Ｘ線源
２移動台
３画像増倍器
４撮像管
５ディジタル画像処理装置
６出力
７モニター
１０，１７メモリ
１１，１３，１４，１５，１６ルックアップテーブル
１２，１８，２０加算器
１９乗算器
１００入力
２００出力

Claims

ディジタル化されたノイズ性の明暗度の値を有し、サンプルＩ_j（ｘ，ｙ）と呼ばれる画素の２次元マトリックスの形で画像のシーケンスの一部を形成する現在の画像と呼ばれる画像（Ｊ_t ^P）内のノイズの時間的フィルタリング方法であって、
上記現在の画像の所定の場所（ｘ，ｙ）の画素に対応するノイズ性のサンプル（Ｉ_t ^P）を再生するため、重み（ｂ_j ^C）付けされたノイズ性の因果的なサンプル（Ｉ_j ^C）と呼ばれる上記シーケンス内の現在のノイズ性サンプル（Ｉ_t ^P）と前のサンプルの線形結合によって、現在のフィルタリングされた因果的なサンプルと呼ばれるフィルタリングされたサンプル（Ｐ_t ^C）を評価する段階からなり、
上記重みは、上記重みが割り当てられた該ノイズ性の因果的なサンプルと、値１を有する重みが割り当てられた上記現在のノイズ性のサンプルの間で、０から１まで変わる明暗度の連続性の確率の係数である、方法。
所定の因果的なノイズ性のサンプルに付けられた上記重みは、上記所定の因果的なノイズ性のサンプルから上記現在のノイズ性のサンプルまでの連続するノイズ性の因果的なサンプルの間の連続性の確率の係数の積であり、
上記因果的なノイズ性のサンプルと現在のサンプルの線形結合は、該サンプルに関係する重みの合計によって正規化されている、請求項１記載の方法。
現時点（ｔ）の上記フィルタリングされた因果的なサンプル（Ｐ_t ^C）に対し、上記線形結合は因果的な積分関係と呼ばれる帰納的な関係を評価することにより形成され、
前の時点（ｔ−１）の上記フィルタリングされた因果的なサンプル（Ｐ_t-1 ^C）と、因果的な利得倍率（Ｋ_t ^C）と呼ばれる倍率と上記現在のノイズ性のサンプル（Ｉ_t ^P）と前の時点の上記フィルタリングされた因果的なサンプル（Ｐ_t- ₁ ^C）との差の積との合計（Ｐ_t-1 ^C＋Ｋ_t ^C×（Ｉ_t ^P−Ｐ_t-1 ^C））を判定する段階を有する請求項２記載の方法。
現時点（ｔ）の該因果的な利得倍率（Ｋ_t ^C）は、上記前の時点（ｔ−１）の因果的な利得倍率（Ｋ_t-1 ^C）の該前の時点の因果的な利得倍率（Ｋ_t-1 ^C）と所謂因果的な連続性の係数（α_t ^C）との合計に対する比によって帰納的な方法で評価される請求項３記載の方法。
該因果的な利得倍率（Ｋ_t ^C）は、上記線形結合の因果的なサンプルに付けられた重みの合計の逆数であり、上記因果的な連続性の係数（α_t ^C）は、現在のサンプル（Ｉ_t ^P）よりも前の因果的なサンプル（Ｉ_t-1 ^C）の重みである、請求項４記載の方法。
上記因果的な積分関係の評価は：
− 上記現時点（ｔ）のノイズ性のサンプル（Ｉ_t ^P）と、上記前の時点のフィルタリングされた因果的なサンプル（Ｐ_t-1 ^C）との間の所謂因果的な差の絶対値（Δ^C）を判定する段階と、
− 上記因果的な連続性の係数（α_t ^C）を、該因果的な差（Δ^C）の減少関数（Ｆ^C）として判定する段階とからなる、請求項４又は５記載の方法。
上記因果的な連続性の係数（α_t ^C）の判定のため、上記減少関数は、引数（Ｚ^C）が上記因果的な差（Δ^C）によって与えられ、平均ノイズに関し上記各サンプルのノイズの差の分散の累積の平方根と一致する倍率（Ｓ_t ^C）によって正規化された所謂因果的な関数（Ｆ^C）である請求項６記載の方法。
上記因果的な減少関数（Ｆ^C）は、上記引数（Ｚ^C）が０と１の間にあるとき、１以下の一定の最大値を有し、上記引数が１より大きいとき、０以上の最小値まで減少する、請求項７記載の方法。
上記フィルタリングの出力は、ノイズのピークを平滑化し、かつ、空間的な動きに起因する不連続点に続くノイズの歯状突起を平滑化するべく選択された上記減少関数（Ｆ^C）の上記最大値及び最小値によって制御される、請求項８記載の方法。
− ノイズ性の現在のサンプルと呼ばれ、現時点（ｔ）に得られる画素のマトリックスの形で画像（Ｊ_t ^P）に所定の場所（ｘ，ｙ）を有する画素Ａ_t（ｘ，ｙ）のノイズ性のディジタル化された明暗度（Ｉ_t ^P）を供給する画像処理装置（５）と、
− 入力（１００）が現在のノイズ性のサンプル（Ｉ_t ^P）を受け、
現在の画素よりも前の画素のノイズ性又はフィルタリングされた因果的サンプル（Ｉ_t-1 ^C；Ｐ_t-1 ^C）と呼ばれ、前の画像のマトリックス内の同一場所（ｘ，ｙ）にあるディジタル化された明暗度を供給するメモリ（１０）と；
因果的な利得倍率（Ｋ_t-1 ^C）を供給するメモリ（１７）と；
因果的な積分式を評価し、上記ノイズ性の現在のサンプル（Ｉ_t ^P）を再生するためフィルタリングされた場合に、現在のフィルタリングされた因果的なサンプル（Ｐ_t ^C）の値を出力（２００）に発生する計算手段とからなる因果的なサブ組立体と呼ばれるサブ組立体とにより構成された請求項１乃至９のうちいずれか１項記載のフィルタリング方法を実行する装置。
上記因果的な積分式を評価し、上記フィルタリングされた因果的なサンプル（Ｐ_t ^C）を発生する計算手段は：
− 現時点ｔの上記ノイズ性のサンプル（Ｉ_t ^P）と、前の時点（ｔ−１）の上記記憶されたフィルタリングされた因果的なサンプル（Ｐ_t-1 ^C）の間の所謂因果的な差の絶対値（Δ^C）を判定する加算器（１２）と、ノイズの分散（σ_B）の倍率（Ｓ_t ^C）によって正規化された上記因果的な差（Δ^C）の減少関数（Ｆ^C）として、因果的な連続性の係数（α_t ^C）を判定するルックアップテーブル（１３，１４，１５）と、
− 記憶された前の時点（ｔ−１）の因果的な利得倍率（Ｋ_t-1 ^C）と、該因果的な利得倍率（Ｋ_t-1 ^C）と因果的な不連続性の係数（α_t ^C）の合計との比率に一致する現在（ｔ）の因果的な利得倍率と呼ばれる倍率（Ｋ_t ^C）を判定するルックアップテーブル（１６）と、
− 前の時点（ｔ−１）の上記フィルタリングされた因果的なサンプル（Ｐ_t-1 ^C）と、上記ノイズ性の現在のサンプル（Ｉ_t ^P）と前の時点の上記フィルタリングされた因果的なサンプル（Ｐ_t-1 ^C）の間の差の該因果的な利得倍率（Ｋ_t ^C）との積との合計（Ｐ_t-1 ^C＋Ｋ_t ^C×（Ｉ_t ^P−Ｐ_t-1 ^C））を判定する加算器（１８，２０）及び乗算器（１９）とからなる、請求項１０記載の装置。
− 画像のデータをディジタル化された明暗度の値を有する画素の２次元マトリックスの形で供給するシステムと、
− 上記画像を表示する表示システムと、
− 上記画像のデータ及び上記表示システムへのアクセスを有する画像処理システムとからなり、上記画像処理システムは請求項１０又は１１記載の装置を含む医用画像の検査装置。