(57)【特許請求の範囲】
【請求項１】
ディジタル化されたノイズ性の明暗度の値を有し、サンプルＩj （ｘ，ｙ）と呼ばれる画素の２次元マトリックスの形で画像のシーケンスの一部を形成する現在の画像と呼ばれる画像（Ｊt P ）内のノイズの時間的フィルタリング方法であって、
上記現在の画像の所定の場所（ｘ，ｙ）の画素に対応するノイズ性のサンプル（Ｉt P ）を再生するため、重み（ｂj C ）付けされたノイズ性の因果的なサンプル（Ｉj C ）と呼ばれる上記シーケンス内の現在のノイズ性サンプル（Ｉt P ）と前のサンプルの線形結合によって、現在のフィルタリングされた因果的なサンプルと呼ばれるフィルタリングされたサンプル（Ｐt C ）を評価する段階からなり、
上記重みは、上記重みが割り当てられた該ノイズ性の因果的なサンプルと、値１を有する重みが割り当てられた上記現在のノイズ性のサンプルの間で、０から１まで変わる明暗度の連続性の確率の係数である、方法。
【請求項２】
所定の因果的なノイズ性のサンプルに付けられた上記重みは、上記所定の因果的なノイズ性のサンプルから上記現在のノイズ性のサンプルまでの連続するノイズ性の因果的なサンプルの間の連続性の確率の係数の積であり、
上記因果的なノイズ性のサンプルと現在のサンプルの線形結合は、該サンプルに関係する重みの合計によって正規化されている、請求項１記載の方法。
【請求項３】
現時点（ｔ）の上記フィルタリングされた因果的なサンプル（Ｐt C ）に対し、上記線形結合は因果的な積分関係と呼ばれる帰納的な関係を評価することにより形成され、
前の時点（ｔ−１）の上記フィルタリングされた因果的なサンプル（Ｐt-1 C ）と、因果的な利得倍率（Ｋt C ）と呼ばれる倍率と上記現在のノイズ性のサンプル（Ｉt P ）と前の時点の上記フィルタリングされた因果的なサンプル（Ｐt-1 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 ））を判定する加算器（１８，２０）及び乗算器（１９）とからなる、請求項１０記載の装置。
【請求項１２】
−  画像のデータをディジタル化された明暗度の値を有する画素の２次元マトリックスの形で供給するシステムと、
−  上記画像を表示する表示システムと、
−  上記画像のデータ及び上記表示システムへのアクセスを有する画像処理システムとからなり、上記画像処理システムは請求項１０又は１１記載の装置を含む医用画像の検査装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、ディジタル画像のシーケンスの画像内のノイズの時間的フィルタリングの方法と、その方法を実行する装置とに関する。
本発明は、特に、Ｘ線ビームが低い明暗度を有し、ディテールを劣化することなく画像からノイズを除去するためにフィルタリングされるべきノイズ性の低コントラスト画像を生成するシステムを用いて、Ｘ線蛍光透視モードで形成された医用画像を処理する際に使用可能である。
【０００２】
時間的フィルタリングは、一次元の時間の関数として、所謂、シーケンスの画像内に所定の場所を有する画素の明暗度の値によって形成された時間的信号を平滑化することにより行なわれる。
【０００３】
【従来の技術】
時間的フィルタリング方法は、アール  イー  カルマン(R.E. Kalman) による論文：“線形フィルタリング及び予測問題の新しいアプローチ(A new approach to linear filtering and prediction problems)”、“アメリカ機械学会論文誌”、基礎工学ジャーナル、シリーズ８２Ｄ、３３−４５ページ、１９６０年によって周知である。
【０００４】
カルマンフィルタは、アプリオリに作成された仮説の関数、シーケンスの中の前の画像内の同じ場所の画素の明暗度の関数、及び、カルマン利得と呼ばれる係数の関数として、シーケンスの画像の瞬時的な画素のフィルタリングされた明暗度を生成する帰納的な式によって定義される。
上記式は二つの帰納的なアルゴリズムを生成する。二つの画像の間に僅かな動きが生じると即座に、上記動きによって、平滑化されるべき時間的な信号の曲線に現れ、明暗度不連続点と呼ばれる立ち上がり又は立ち下がりエッジが発生するという問題が存在する。
【０００５】
第１のアルゴリズムの場合、カルマン利得は完全に一定であるよう選択され；これは、動きによって誘起された明暗度の不連続点に影響を与える指数関数的な縞（ストリーキング）を発生させる。かくして、明暗度信号に段差を生じさせるノイズ性の元画像内の小さい対象物、例えば、素早く移動させ得るカテーテルは、段差の側面がフィルタリングによって変形されるので、フィルタリングされた画像から消える可能性がある。上記アルゴリズムは小さい対象物を消去する。
【０００６】
第２のアルゴリズムの場合、カルマン利得は、所定の時点で所定の場所を占める画素に対し観察されたノイズ性の明暗度と、前の時点でフィルタリングされた明暗度との差の関数である。その結果として、時間的信号は、不連続点の前で平滑化されるが、しかし、明暗度の不連続点の後でフィルタリングされることがないので、残留ノイズが明暗度の不連続点のエッジの後ろに存在する。
【０００７】
【発明が解決しようとする課題】
従って、周知の時間的フィルタリング方法は、動いている小さい対象物を表わす非常にノイズの多い画像の系列に効果的に適用できないという欠点がある。
かくして、周知の時間的フィルタリング法は、例えば、カテーテルのような非常に小さい径の器具が観察される領域内で挿入又は移動されている間に医学的な手術を実時間で追跡するため実行されたＸ線蛍光透視モードで得られた画像のシーケンスに、時間的フィルタリング法が適用されたとき生じる幾つかの主要な問題を解決しない。
【０００８】
第１の問題は、蛍光透視法のＸ線ビームの非常に低い明暗度に起因して、シーケンスの画像は非常に多くのノイズを含み、屡々、ノイズのピークを含むことである。
第２の問題は、シーケンスの各画像が短い時間的ギャップによって次の画像から分離されているので、カテーテルのような小さい器具の移動のような重要な事象が、画像と別の画像の間に生じ得ることである。小さい対象物の動きが出現する画像のフィルタリング対象物を歪め、又は、消すことは許されない。
【０００９】
本発明の目的は、
−  厳密に実時間的に実行され、即ち、観察された瞬時的な画像の出現の直後に、画像遅延を伴うことなく、かつ、シーケンス内の画像レートとは無関係にフィルタリングされた瞬時的な画像を出力し、
−  時間的明暗度信号の不連続性エッジの後で、不連続性エッジを弱めることなく残留ノイズを低減し、
−  ノイズピークを低減することが可能であり、
−  動いている小さい対象物を消去又は歪めることなく、
ディジタル化された画像のシーケンスの連続する画像内のノイズを低減する時間的フィルタリング方法を提供することである。
【００１０】
【課題を解決するための手段】
上記目的は、ディジタル化されたサンプルと呼ばれるノイズ性の明暗度の値を有する画素の２次元マトリックスの形で画像のシーケンスの一部を形成する現在の画像と呼ばれる画像内のノイズの時間的フィルタリング方法であって：
上記現在の画像の所定の場所の画素に対応するノイズ性のサンプルを再生するため、重み付けされたノイズ性の因果的なサンプルと呼ばれる上記シーケンス内の現在のノイズ性サンプルと前のサンプルの線形結合によって、現在のフィルタリングされた因果的なサンプルと呼ばれるフィルタリングされたサンプルを評価する段階からなり、
上記重みは、上記重みが割り当てられた該ノイズ性の因果的なサンプルと、値１を有する重みが割り当てられた上記現在のノイズ性のサンプルの間で、０から１まで変わる明暗度の連続性の確率の係数である方法によって解決される。
【００１１】
上記フィルタリング方法を実行する装置は：
−  ノイズ性の現在のサンプルと呼ばれ、現時点に得られる画素のマトリックスの形で画像に所定の場所を有する画素のノイズ性のディジタル化された明暗度を供給する画像処理装置と、
−  入力が現在のノイズ性のサンプルを受け、現在の画素よりも前の画素のノイズ性又はフィルタリングされた因果的サンプルと呼ばれ、前の画像のマトリックス内の同一場所にあるディジタル化された明暗度を供給するメモリと；因果的な利得倍率を供給するメモリと；因果的な積分式を評価し、上記ノイズ性の現在のサンプルを再生するためフィルタリングされた場合に、現在のフィルタリングされた因果的なサンプルの値を出力に発生する計算手段とからなる因果的なサブ組立体と呼ばれるサブ組立体とにより構成されている。
【００１２】
上記装置は、容易に作動させることができ、画像のディテールを劣化させることなく、ノイズの効果的、厳密かつ正確な実時間の時間的フィルタリングを提供する利点が得られる。
【００１３】
【発明の実施の形態】
以下、添付図面を参照して本発明を詳細に説明する。
（１）Ｘ線装置
図１には、Ｘ線源１と、患者のための移動台２と、マイクロプロセッサからなるディジタル画像処理装置５にデータを提供する撮像管４に接続された画像増倍器３とにより構成されたディジタルＸ線透視図システムが示されている。ディジタル画像処理装置５は、複数の出力を有し、その中の一つの出力６は、Ｘ線透視画像又は明暗度画像のシーケンスを表示するモニター７に接続されている。
【００１４】
ディジタル化されたＸ線透視画像は、８ビット又は１０ビットで符号化された５１２×５１２画素又は１０２４×１０２４画素によって構成可能である。各画素は、従って、２５６又は１０２４明暗度レベルの中の一つのレベルを有する。例えば、画像の暗い領域は低い明暗度レベルを有し、明るい領域は高い明暗度レベルを有する。
【００１５】
ディジタル化された画像は蛍光透視モードで得られる。本発明は、特に、血管造影画像のフィルタリングに使用することが可能である。
いずれにしても、本発明は、ディジタル画像のシーケンスを取得するため使用される方法、或いは、画像の中に表わされた対象物の性質を考慮していないが、ノイズを除去するためシーケンス内の画像をフィルタリングすることだけに関係している。
【００１６】
（２）時間的ノイズフィルタリング方法
本発明は、ノイズ性のディジタル画像のシーケンスの各画像内で連続的にノイズの時間的フィルタリングを行う方法を提案する。上記方法は、最後の観察されたノイズ性画像内でフィルタリング演算を行う。従って、上記方法は厳密に実時間で行われる。
【００１７】
上記方法はシーケンス内の画像に記録された対象物の動きを検出することが可能である。これは、適応的な反復フィルタリングステップに基づいている。
図２を参照すると、ノイズの時間的フィルタリング方法は、最初に、所定のレートで形成されたノイズ性画像のシーケンスの取得及びディジタル化を行う。
時間的に最も離れた時点から最も近い時点までの間で、上記シーケンスは：
−  ｋ0 がシーケンス内の画像の数よりも１だけ小さい整数に一致し、ｔがフィルタリングされる画像の到達の時点を表わすとき、最初の時点ｊ＝ｔ−ｋ0 から時点ｊ＝ｔ−１までの間に形成され、Ｊj C で示され、“因果的な画像”と呼ばれる多数の過去の明暗度画像と、
−  時点ｊ＝ｔに到達し、現在の画像と呼ばれるフィルタリングされる画像Ｊt P とにより構成される。
【００１８】
ｊが画像の時点に対応するインデックスであるとき、ノイズ性のディジタル化された各画像Ｊj は、画素の２次元マトリックスＡj （ｘ，ｙ）からなり、各画素は、マトリックス内の座標ｘ，ｙと、明暗度レベルのスケール上の明暗度レベルＩj （ｘ，ｙ）とにより表わされる。
所謂ノイズ性のディジタル化された画像は、図３の（Ａ）及び（Ｂ）に示されているように、時間τの関数として、画素Ａj （ｘ，ｙ）の種々の明暗度レベルＩj （ｘ，ｙ）により形成され、同図において、明暗度Ｉが縦軸にプロットされ、時間τが横軸にプロットされている。本発明による時間的なフィルタリング方法は、図４の（Ａ）及び（Ｃ）に表わされているように、フィルタリングされた時間的信号Ｐ（τ）を得るため、時間的信号Ｉ（τ）に影響を与えるノイズをフィルタリング又は平滑化することを目的としている。以下、時間的信号Ｉ（τ）を構成する明暗度を“サンプル”と呼ぶ。過去又は因果的な明暗度はＩj C と呼ばれ、現在の明暗度はＩt P と呼ばれる。
【００１９】
本発明による方法は、画像シーケンスに発生した主な動きが、例えば、並進補正法のような従来技術によって周知の動き補償法によって既に補償されている場合に関係している。しかし、大量の計算時間を必要とする動き評価法は避けるべきである。本発明は厳密に実時間的な方法の提供を目的としているので、このようなタイプの方法は避ける必要がある。
【００２０】
たとえ、動き補償法が画像のシーケンスに適用されるとしても、動き補償法は完璧ではないので、全ての動きが現実的に本当に補償されることが保証された訳ではないことに注意すべきである。
この事実を考慮に入れるため、現在の方法は、画像シーケンス内の対象物の小さい残留又は補償不可能な局部的な動きを取り扱うことによって、時間的信号Ｉ（τ）のノイズ上のフィルタリング演算を実行する。
【００２１】
上記方法は、現時点ｔで得られたノイズ性サンプルＩt P 上でフィルタリング演算を実行し、現時点ｔのノイズ性サンプルＩt P と、時点ｔ−ｋ０からｔ−１の間に観察されたノイズ性因果サンプルＩt-k0C ，．．．，Ｉt-3 C ，Ｉt-2 C ，Ｉt-1 C と呼ばれる前のサンプルとを利用することにより、フィルタリングされたサンプルＰt C を生成する。
【００２２】
時間的フィルタリングは、好ましくは、２次元マトリックスの種々のｘ，ｙ座標の各画素Ａj （ｘ，ｙ）に対し別個に行なわれる。
図３の（Ａ）には、例えば、小さい局部的な動きに起因したＩt-3 C とＩt-4 C の間の明暗度の不連続点Ｄからなるノイズ性時間的信号Ｉ（τ）が表わされている。上記小さい局部的な動きは、時点ｔ−４と時点ｔ−３の間で時間的信号Ｉ（τ）の立ち上がりエッジと、時点ｔ−３とｔの間の大きい信号“安定状態”と、時点ｔ−４より前の時点の“抑圧状態”として表わされる。
【００２３】
本発明による時間的フィルタリングは、Ｄのような明暗度の不連続点に適応する。かくして、図４の（Ａ）及び（Ｃ）は、本発明による方法の種々の変形例によって平滑化された時間的信号Ｐ（τ）を表わしている。
上記時間的信号は、ノイズに起因して、シーケンスの画像内で同一の座標（ｘ，ｙ）を有する画素の明暗度レベルが絶えず変化するため、時点ｔ−８とｔ−４の間に小さい鋸波を更に有し、これにより、降雪のような斑な局面が現れる。フィルタリングは、時間的信号Ｉ（τ）を所定の画素に関して一定の平均値に平滑化することにより、かかる局面を抑圧する。
【００２４】
更に、図３には、例えば、時点ｔ−４とｔ−２の間、即ち、時点ｔ−３に現れたノイズのピークＤ’を有するノイズ性の時間的信号Ｉ（τ）が表わされている。かかるノイズのピークＤ’は、実際の空間的な動きの現象に対応していない点で、図３の（Ａ）のＤのような明暗度の不連続点とは区別される。かかるノイズのピークは、本発明によって提案された方法を用いてフィルタリングすることができる。
【００２５】
現在のノイズ性サンプルＩt P に対応しているフィルタリングされたサンプルＰt C は、以下の式：
【００２６】
【数１】
0003865821.tif 000002
【００２７】
のように表現される。フィルタリングされたサンプルＰt C 、又は、時点ｔでフィルタリングされた時間的信号は、因果的な時点ｊ＝ｔ−ｋ0 からｊ＝ｔ−１で観察されたノイズ性サンプル及び時点ｔで観察されたノイズ性サンプルの線形結合である。
上記フィルタリングされた信号Ｐt C の式（１）において、重みｂj C が因果的なサンプルＩt-1 C 、Ｉt-2 C 等に割り当てられる。
【００２８】
上記（１）式の分母は正規化係数であり、その正規化係数の存在は、フィルタリングされた信号Ｐt C の平均値がノイズ性信号Ｉt P の平均値と一致することを保証するため、種々のサンプルに付けられた重みの合計は１に一致すべきであるという事実に基づいている。
時間的フィルタリング方法は、因果的なサンプルの現在のサンプルとの線形結合に基づいている。上記方法において、因果的なサンプルに付けられた重みｂj C は特定の形式を有し、現在のサンプルに付けられた重みは１である。
【００２９】
所定の因果的なサンプルに付けられた重みは、所定の因果的なサンプルと現在のサンプルの間の明暗度の連続性の確率である。
一般的に言うと、因果的なサンプルＩt-k C に割り当てられるべき因果的な重みｂt-k C の式は、次式（２）：
【００３０】
【数２】
0003865821.tif 000003
【００３１】
によって与えられ、式中、α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 ）：
【００４０】
【数３】
0003865821.tif 000004
【００４１】
によって与えられる全く別の利得倍率である。
式（１）及び（２）と、式（３C ）及び（４C ）の間の推移に対し、因果的な利得は、因果的な重みの逆数：
【００４２】
【数４】
0003865821.tif 000005
【００４３】
であり、更に、上記の如く：
α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 ）
と、因果的な利得倍率である
【００４４】
【数５】
0003865821.tif 000006
【００４５】
と、因果的な連続性の係数（５C ）である
【００４６】
【数６】
0003865821.tif 000007
【００４７】
からなるシステムに基づき、ここで、Ｓt C は正規化係数を表わし、ＦC は以下で定義されている関数である。
時間的にｔ−ｋ0 で始まる式（３C ）の反復により、上記定式化されたように式（１）が厳密に得られる。
積分式の定式化（３C ）は、かくして、連続したサンプルの間の明暗度の連続性の関数であり、現在のサンプルに至るまで重みを計算された連続するサンプルの間の連続性の確率の積（２）によって定められた特定の重みと関係したサンプルの線形結合（１）に直接対応する。
【００４８】
・因果的な不連続性の係数αt C の定式化
因果的な不連続性の係数αt C は、差：
ΔC ＝｜Ｉt P −Ｐt-1 C ｜                （６C ）
の関数として定義される。
上記差ΔC が、図３の（Ａ）のＤのような不連続点に関係しているか、或いは、ノイズだけに関係しているかどうかを判定することにより利点が得られる。
【００４９】
ノイズの関与を判定するため、上記差ΔC は、上記不連続性ΔC の各サンプルに関するノイズの分散を考慮した係数Ｓt C によって正規化される。σB と呼ばれ、各サンプルに対する明暗度レベルで測定されたノイズの標準偏差は、あらゆる周知の方法を用いて評価し、或いは、アプリオリに評価することが可能である。
【００５０】
ノイズの分散σB 2 は、画像内で想定された領域のノイズの平方の平均偏差であり、ノイズに起因した明暗度の偏差は、平均局部明暗度に関して測定されている点に注意する必要がある。
標準ノイズ偏差の判定は、図１に示したシステムにおいて、信号のディジタル化の前に、低明暗度レベル及び高明暗度レベルの圧縮の規則に従って画像の変換が導入されることによって複雑化するので、ノイズは、非常に暗い領域又は非常に明るい領域の場合、中間の明暗度の領域よりも見えにくくなるということを考慮する。標準ノイズ偏差σB は、全ての画素の全ての時点における明暗度の関数である。
【００５１】
従って、標準ノイズ偏差の判定は、対応する領域の平均明暗度の関数として、画像の各画素上で局部的に実行する必要がある。かくして、画像自体又は信号に起因したノイズの分離が行なわれるように、ノイズの分散及び観察された信号の平均値が局部的に計算されるので、ノイズは空間的に表わされている。
図５の（Ａ）を参照して一例として提案された方法に従って、画像Ｊj 内で、典型的に、各画素Ａj （ｘ，ｙ）の周囲の５×５又は３×３画素の区域からなる小さいウィンドウＷj が定義される。上記各小さいウィンドウＷj （ｘ，ｙ）において、ノイズの分散及び観察された平均明暗度が判定される。
【００５２】
このため、ｎが小さいウィンドウＷj （ｘ，ｙ）の画素数を表わす場合に、画素Ａj （ｘ，ｙ）の周辺に中心が置かれた小さいウィンドウの各位置に対し、以下の如く：
−  インデックスｋ及びｌが、画像Ｊj 内の矩形状の座標Ｘ及びＹに関してウィンドウ内にあると考えられる画素の座標を表わす場合に、小さいウィンドウＷj （ｘ，ｙ）の各画素の明暗度Ｉj （ｋ，ｌ）と、
−  この位置で、関係式：
【００５３】
【数７】
0003865821.tif 000008
【００５４】
によって与えられるウィンドウ内の平均明暗度ｍW （ｘ，ｙ）と、
−  上記ウィンドウ内で、総和の記号Σの後ろでウィンドウの各画素の明暗度と平均値ｍW （ｘ，ｙ）の差の平方が与えられ、（ウィンドウ内の画素数ｎ）分の（上記総和）である平均値が与えられる以下の関係式：
【００５５】
【数８】
0003865821.tif 000009
【００５６】
によって与えられるノイズの分散σBW（ｘ，ｙ）の平方とが定められる。
図５の（Ｂ）を参照すると、次いで、２次元のヒストグラムが形成され、同図において、例えば、１０グレイレベル×１０グレイレベルの所定の精度で、明暗度の平均値ｍW （ｘ，ｙ）が水平軸上にプロットされ、一方、分散の平均値が垂直軸上にプロットされている。水平軸上にプロットされた平均明暗度ｍW （ｘ，ｙ）の各間隔には、分散の頻度、又は、所定の分散が上記平均明暗度の間隔に生じた回数の平均が記載されている。説明の簡単化のため図５の（Ｂ）には、一例として１０と２０の間の平均値ｍW （ｘ，ｙ）の値が以下のように示されている：
値１を有する分散の頻度Ｎ01は０；
値４を有する分散の頻度Ｎ11は１０；
値９を有する分散の頻度Ｎ21は２０；
値１６を有する分散の頻度Ｎ31は３５；
値２５を有する分散の頻度Ｎ41は２５；以下同様である。
【００５７】
提案された方法は、分散σ2 BW（ｘ，ｙ）の関数として垂直軸上にプロットされた分散の出現頻度Ｎの曲線が、想定された平均明暗度ｍW （ｘ，ｙ）のあらゆる間隔に対し描かれている図５の（Ｃ）を参照して続けられる。分散の関数として、分散の可能な値を与えながら最大値を通過して延在する頻度Ｎの曲線は、平均明暗度の各間隔に対し定められる必要がある。
【００５８】
かくして、ノイズ分散が、各画素Ａj （ｘ，ｙ）に属する平均明暗度の値の関数として得られる。この方法によって、画像自体から抽出されるべきノイズに関する情報が得られる。
次いで、標準ノイズ偏差σB が、ノイズ分散の平方根：
【００５９】
【数９】
0003865821.tif 000010
【００６０】
の計算により求められる。
差ΔC に加わるノイズ寄与の分散は、安定状態の場合、次式：
【００６１】
【数１０】
0003865821.tif 000011
【００６２】
のように定式化され、この値は、ΔC の正規化の近似として使用される。差ΔC の正規化係数は、以下の関係式：
【００６３】
【数１１】
0003865821.tif 000012
【００６４】
によって与えられる。
式（５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 ：
【００６８】
【数１２】
0003865821.tif 000013
【００６９】
によって表わされる。
関数Ｆ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 ）は、好ましくは、上記領域で、正規化倍率によって加えられるようなガウシアン形状を有するので：
【００７１】
【数１３】
0003865821.tif 000014
【００７２】
のように表わされる。
・例２：図６の（Ｂ）に示された例
関数Ｆ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 ）を用いて実行される。
図４の（Ａ）を参照すると、因果的なフィルタリングは、不連続点Ｄの左側で、ｔ−４からｔ−ｋ0 の方向に時間的信号の平滑化を行い；次いで、不連続点ＤのためサンプルＩ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 の値とが入力されると即座に、値：
【００８５】
【数１４】
0003865821.tif 000015
【００８６】
を直接的に供給するため表にされた値を含むルックアップテーブル１３を更に有する。値Ｓ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 ：
【００８９】
【数１５】
0003865821.tif 000016
【００９０】
を出力する。
ルックアップテーブル１５は、表にされた形で選択された関数ＦC （ＺC ）を含んでいる。ルックアップテーブル１５の入力は、引数ＺC を受け、ルックアップテーブル１５は、連続性係数αt C ：
【００９１】
【数１６】
0003865821.tif 000017
【００９２】
を出力する。
ルックアップテーブル１６は、メモリ１７から得られた前の時点の利得と、ルックアップテーブル１５から得られた連続性係数αt C とに基づいて、時点ｔの利得を供給する表にされた式を含む。上記利得：
【００９３】
【数１７】
0003865821.tif 000018
【００９４】
は、ルックアップテーブル１６を介して与えられ、一方で、前の時点の利得Ｋ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 ）の種々の実現可能な例を表わす図である。
【図７】時間的なフィルタリング方法を実行する装置を機能的なブロックの形で表わした図である。
【符号の説明】
１    Ｘ線源
２    移動台
３    画像増倍器
４    撮像管
５    ディジタル画像処理装置
６    出力
７    モニター
１０，１７    メモリ
１１，１３，１４，１５，１６    ルックアップテーブル
１２，１８，２０    加算器
１９    乗算器
１００    入力
２００    出力
【図１】
0003865821.tif 000019
【図２】
0003865821.tif 000020
【図３】
0003865821.tif 000021
【図４】
0003865821.tif 000022
【図５】
0003865821.tif 000023
【図６】
0003865821.tif 000024
【図７】
0003865821.tif 000025