(57)【特許請求の範囲】
【請求項１】
線状陰影及びノイズが含まれている画像に対して少なくとも２種類の空間的な画像処理を施す手段と、
前記２種類の空間的な画像処理の両方の結果に基づいて、前記線状陰影の存在確率を画素ごとに求める手段と、
前記画像又はそれから画像処理により導かれる画像に対して時間的に連続している複数フレームの画像を加算平均する時間的な画像処理を施す手段と、
前記時間的な画像処理の特性を前記線状陰影の存在確率に従って画素ごとに変える手段とを具備したことを特徴とする画像処理装置。
【請求項２】
線状陰影及びノイズが含まれている画像に対して少なくとも２種類の空間的な画像処理を施す空間的画像処理手段と、
前記２種類の空間的な画像処理の両方の結果に基づいて、前記線状陰影の存在確率を画素ごとに求める手段とを具備し、
前記空間的画像処理手段は、第１の画像処理系と第２の画像処理系とからなり、前記第１の画像処理系は、前記線状陰影を強調するフィルタを有し、前記第２の画像処理系は、前記ノイズを低減するフィルタを有することを特徴とする画像処理装置。
【請求項３】
線状陰影及びノイズが含まれている画像に対して少なくとも２種類の空間的な画像処理を施す空間的画像処理手段と、
前記２種類の空間的な画像処理の両方の結果に基づいて、前記線状陰影の存在確率を画素ごとに求める手段とを具備し、
前記空間的画像処理手段は、第１の画像処理系と第２の画像処理系とからなり、
前記第１の画像処理系は、前記線状陰影を強調するフィルタと、前記線状陰影が強調された画像から空間周波数が比較的高い成分を抽出するハイパスフィルタと、前記線状陰影が強調された画像から前記ハイパスフィルタで処理された画像を差分する第１の差分回路と、前記第１の差分回路で差分された画像を極性に従って２値化する回路とを有し、
前記第２の画像処理系は、前記ノイズを低減するフィルタと、前記ノイズが低減された画像から空間周波数が比較的低い成分を抽出するローパスフィルタと、前記ノイズが低減された画像から前記ローパスフィルタで処理された画像を差分する前記第２の差分回路と、前記第２の差分回路で差分された画像を極性に従って２値化する回路とを有することを特徴とする画像処理装置。
【請求項４】
前記存在確率の値に基づいて、前記第１の画像処理系と、前記第２の画像処理系とのいずれか一方を選択して出力する手段をさらに備えたことを特徴とする請求項２又は３記載の画像処理装置。
【請求項５】
線状陰影及びノイズが含まれている画像に対して少なくとも２種類の空間的な画像処理を施す手段と、
前記２種類の空間的な画像処理の両方の結果に基づいて、前記線状陰影の存在確率を画素ごとに求める手段と、
前記画像又はそれから画像処理により導かれる画像に対して時間的に連続している複数フレームの画像を加算平均する時間的な画像処理を施す手段と、
前記時間的な画像処理の特性を前記線状陰影の存在確率に従って画素ごとに変える手段とを具備し、
前記時間的な画像処理の特性を変える手段は、前記存在確率が比較的高い画素に対しては、前記時間的な画像処理の平滑化特性を比較的低くし、前記存在確率が比較的高い画素に対しては、前記時間的な画像処理の平滑化特性を比較的高くすることを特徴とする画像処理装置。
【請求項６】
線状陰影及びノイズが含まれている画像に対して少なくとも２種類の空間的な画像処理を施すステップと、
前記２種類の空間的な画像処理の両方の結果に基づいて、前記線状陰影の存在確率を画素ごとに求めるステップと、
前記画像又はそれから画像処理により導かれる画像に対して時間的な画像処理を施すステップと、前記時間的な画像処理の特性を前記確率分布に従って画素又は近隣グループごとに変えるステップとを具備したことを特徴とする画像処理方法。
【請求項７】
入力画像の中から線状陰影領域と他の領域とを選別する第１選別手段と、
前記第１選別手段とは異なる種類の選別方法で前記入力画像の中から線状領域を選別する第２選別手段と、
前記第１選別手段と前記第２選別手段との出力に基づいて、線状陰影の存在確率を求める存在確率決定手段と、
前記入力画像にノイズ低減処理を行うものであり、前記存在確率決定手段の出力に基づいて３種類以上のフィルタ係数を切り換えるノイズ低減処理手段とを具備することを特徴とする画像処理装置。
【請求項８】
前記ノイズ低減処理手段は、前記入力画像中で前記存在確率が低くなるに従ってノイズ低減効果が高くなるようにフィルタ係数を変えるものであることを特徴とする請求項７記載の画像処理装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、血管やガイドワイヤ等の線状陰影を含む例えばディジタルＸ線画像を処理する画像処理装置及び方法に関する。
【０００２】
【従来の技術】
Ｘ線を利用した医療技術として、透視下におけるカテーテル治療が盛んに行われているが、Ｘ線透視は撮影に比べて線量を低減して患者被曝を低減していることから、画像上に重畳するノイズが大きく、カテーテル、その他、ガイドワイヤや血管等の線状陰影が、背景ノイズに邪魔されて見えにくくなるという問題があった。
【０００３】
また、線量を増加して透視を行うと、背景ノイズは相対的に小さくなり、線状陰影の視認性は高くなるが、患者や医療技術者の被曝量が大きくなるという欠点がある。そのため、画像処理によってノイズを低減させたりコントラストを向上させる技術が現在に至っても研究され続けている。
【０００４】
最も一般に用いられるノイズ低減のための画像処理技術は、時間的に連続している複数フレームの画像を加算平均する手法である。この技術はすでに公知であり、リカーシブフィルタと呼ばれる回帰フィルタも、その応用技術として広く一般に使用されている。
【０００５】
しかし、このような時間的フィルタは、動きの少ない対象に対してはノイズだけを効果的に低減できるが、心臓血管に挿入されたガイドワイヤのような比較的動きのある対象に対しては、そのコントラストを薄めたり、残像を生じさせたりする欠点がある。そのため、被写体の動きを検出し、動きのあるところではフィルタ特性を変更するような技術が開発されている（特許第２５０８０７８号公報、特公平６−６９４４７号公報、特開平３−１９８８３６号公報、特開平６−４７３０５号公報、特開平７−７９９５６号公報、特開平８−２５５２３８号公報）。
【０００６】
しかしながら、これらの技術では、被写体の動きの検出を現時点の画像と過去の画像の差分に基づいて検出しているため、ノイズによる画素値の変動も被写体の動きとして検出されてしまうという問題があった。特に細径のカテーテルやコントラストの低いガイドワイヤを使用するために、動き検出精度を高くすると、ノイズ低減が効果的に行うことができなくなるという問題があった。
【０００７】
一方、ノイズ低減やコントラストを向上させる目的で、空間フィルタも一般に使用されている。この空間フィルタは、近傍局所の画素の加重平均をとってノイズを低減しスムージングしたり、空間的に微分することによってコントラストを強調する等の様々な使い方のできるものである。しかし、スムージングによってノイズだけでなく、肝心のガイドワイヤやカテーテルのような線状陰影までもがぼけてしまい、そのコントラストが損なわれてしまったり、あるいはコントラスト強調の場合には、ノイズも強調されてしまうという欠点がある。
【０００８】
このため、空間的に血管影をパターン認識する装置（特開平４−１２２３５５）、あるいは被写体のエッジ領域を検出してエッジ領域にはコントラスト強調を、そうでない領域にはスムージングを施す装置（特開昭６０−２４５０８４号公報）などが開発されている。
【０００９】
しかし、前者のパターン認識の方法は、１画素単位に多くの平均回路や比較回路を必要とするため、透視処理などのリアルタイム演算を行うには非常に多くの回路を必要とし、そのうえ、ノイズを低減させる効果はない。
【００１０】
また、後者のエッジ強調の方法は、単に１ビットの情報によってエッジ領域と背景領域を区別するため、誤認識があった場合には、極端にスムージングをかけられた部分に隣接して誤認識されたノイズがコントラスト強調されて出力されることなどがあり、誤認識が敏感に出力画像のアーチファクトとなって画質を損ないやすい。
【００１１】
【発明が解決しようとする課題】
上述したように、従来の技術では、動きのある被写体の画像について、残像低減とノイズ低減の効果を同時に高くすること、また、コントラスト強調とノイズ低減の効果を同時に高くすることが困難であった。
【００１２】
本発明は、残像低減とノイズ低減の効果を同時に高くすること、また、コントラスト強調とノイズ低減の効果を同時に高くすることを実現し、Ｘ線医療におけるガイドワイヤ、カテーテル、血管像等の陰影を良好に表示することのできる画像処理装置及び方法を提供することを目的とする。
【００１３】
【課題を解決するための手段】
本発明は、線状陰影及びノイズが含まれている画像に対して少なくとも２種類の空間的な画像処理を施し、この２種類の空間的な画像処理の結果に基づいて線状陰影の存在確率を画素ごとに求め、この存在確率に従って時間的な画像処理の特性を画素ごとに変えることを特徴としたものである。
【００１４】
【発明の実施の形態】
以下、図面を参照しながら、本発明による画像処理装置を好ましい実施形態により説明する。ここでは、画像処理の対象として、血管、ガイドワイヤ、カテーテル等の線状陰影と背景ノイズとを含むＸ線透視画像を例に説明するものとする。
【００１５】
図１に、本実施形態に係る画像処理装置を組み込んだＸ線透視装置の構成を示している。Ｘ線透視装置本体１は、Ｘ線源２からばく射されたＸ線３を寝台５上に載置された被検体４の例えば胸部４ａに向けて照射し、そして、胸部４ａを透過するうちに、各部のＸ線吸収率の違いによって陰影が付けられた透過Ｘ線３ａをイメージインテンシファイア６に導入し、そこで変換され増幅された光学像をＴＶカメラ７で撮像するように組まれている。このＴＶカメラ７の出力信号は、図示しないプリアンプ、アナログディジタル変換器等を通って、ディジタルＸ線画像信号８に変換される。
【００１６】
そして、このディジタルＸ線画像信号８は、画像処理装置９で適当に処理されて、この処理により、背景ノイズは低減され、線状陰影は相対的に強調された後、ＣＲＴディスプレイ１０に濃淡表示される。
【００１７】
ところで、血管、その他、体内に挿入されるカテーテルやガイドワイヤ等は、その外形に応じて細い線状の陰影でディジタルＸ線画像上に表現される。従来技術でも述べたように、このような細い線状陰影は、その性質、例えば空間周波数及び時間周波数がノイズと同様に高く、このためノイズ低減処理の影響を受けやすく、ノイズと一緒に、消えてしまったり、ぼけたり、そのコントラストが損なわれてしまうというやっかいな問題がある。このような問題に対して、従来では、Ｘ線の照射線量を高くして、信号雑音比（Ｓ／Ｎ）を向上させることで対処していた。
【００１８】
これに対して、本実施形態では、Ｘ線の照射線量を高くすることなく、さらには低下させても、線状陰影のコントラストを損なわず、背景ノイズだけを効果的に低減しようとするものである。このように、背景ノイズが低減され、しかも線状陰影が十分なコントラストで観察できると、この透視下で行われるカテーテル治療やバイオプシーなどの医療行為の精度が格段に向上するものと期待されている。
【００１９】
本実施形態は、基本的に、同じディジタルＸ線画像に対して、線状陰影を強調する効果の認められる少なくとも２種類の画像処理を施し、これら２種類の画像処理による２種類の結果から総合的に判断して、線状陰影の存在確率を画素ごとに求め、この存在確率に従って時間フィルタの特性を画素ごとに変えていこうとするものであり、つまり、線状陰影の存在確率が比較的高い画素に対しては、時間フィルタのスムージング効果を低く抑えて線状陰影のコントラストの低下を防ぎ、一方、線状陰影の存在確率が比較的低く、背景ノイズである確率が比較的高い画素に対しては、時間フィルタのスムージング効果を高くしてノイズを効果的に低減しようとするものである。
【００２０】
図２に、この画像処理装置９の処理手順を示している。なお、この処理は、コンピュータで読み取り可能な記憶媒体に書き込まれたプログラムによりソフト的に実現してもよいし、また、図２の各ブロックに対応する回路を処理手順に従って組み合わせてハード的に実現するようにしてもよい。なお、この画像処理装置９は、リアルタイムで入力画像を処理して表示できるように構成する方が好ましい。
【００２１】
画像処理装置９に対して、入力画像として、ディジタルＸ線画像Ｉ0 が、ＴＶカメラ７から所定の周期で連続的に供給されてくる。そして、まず、このディジタルＸ線画像Ｉ0 から、線状陰影の存在確率を画素ごとに求める（全ての画素の存在確率のまとまりを、「存在確率分布」と称する）。図２では、この存在確率分布を、“Ｆ３”として示している。本実施形態では、この存在確率分布Ｆ３を求めるために、Ａ系列とＢ系列との２種類の画像処理を採用している。まず、Ａ系列の処理手順について説明する。
【００２２】
Ａ系列は、処理順に、Ａ１，Ａ２，Ａ３，Ａ４，Ａ５の５つのステップからなる。まず、ステップＡ１では、ディジタルＸ線画像Ｉ0 に対して、最小値フィルタリングという処理を施して、最小値画像Ｉmin を得る。最小値フィルタリングは、周知の通り、図３に示すように、例えば２×２マトリクスのフィルタマスクを重ならないように移動し、その各位置で、マスク内の４つの近傍画素の中の最小の画素値を選出し、その最小値を最小値画像Ｉmin の中の対応する画素の画素値とするいう処理である。
【００２３】
なお、この処理により、マトリクスサイズを１／４に縮小する、例えばディジタルＸ線画像Ｉ0 のマトリクスサイズが１０２４×１０２４の大きさである場合には、５１２×５１２のマトリクスサイズに縮小されるが、これ以降の処理を５１２×５１２のマトリクスサイズで行うことにより回路の規模を縮小することができる。
このような最小値フィルタリングによると、線状陰影のコントラストが、多少強調されることになる。
【００２４】
次にステップＡ２では、Ａ１で得られた最小値画像Ｉmin に対して、ハイパスタイプ（高域通過型）の空間フィルタリングを施し、空間周波数が比較的低い成分が濾過された最小値画像Ｉhighを得る。
【００２５】
このステップＡ２でハイパスフィルタされた最小値画像Ｉhighと、先のステップＡ１で得られた最小値画像Ｉmin に対して、ステップＡ３で、次の式（１）ような線形演算、つまり加重サブトラクション（加重を掛けて差分を取る）を実行することにより、差分画像Ｉsub1を得る。
【００２６】
Ｉsub1＝Ｗ1 ×Ｉmin −Ｗ2 ×Ｉhigh＋Ｃ    …（１）
図４（ａ）に、ハイパスフィルタされた最小値画像Ｉhighと、最小値画像Ｉmin それぞれの空間的なプロフィール（画素値の位置関数）を示し、同図（ｂ）に、上記式（１）の差分処理により得られた差分画像Ｉsub1の空間的なプロフィールを示している。上記式（１）の係数Ｗ1 、Ｗ2 、Ｃを適当に設定することにより、図４（ｂ）に示すように、差分画像Ｉsub1上で、線状陰影である確率が比較的高いと推定される部分だけを、それ以外の部分に対して、極性を反転させることができる。
【００２７】
このように線状陰影である確率が比較的高いと推定される部分と、それ以外の部分とが極性で識別された差分画像Ｉsub1を、その極性に従って２値化し、線状陰影である確率が比較的高いと推定される部分が“１”で、またそれ以外の部分が“０”で表現された２値画像Ｉbi1 を得る（ステップＡ４）。
【００２８】
この２値画像Ｉbi1 には、殆どの場合で背景ノイズが含まれているので、これを取り除く処理が必要になる。ここで、図５（ａ）に示しているように、ノイズは孤立的であり、一方、線状陰影は連結性が高い。このような性質の違いを利用して、ここでは、２値画像Ｉbi1 に対してパターンマッチングを実行する（ステップＡ５）。パターンマッチングについては周知の通り、図５（ｂ）〜（ｅ）に示すような複数種類の線状陰影パターンを予め用意しておき、このパターンに一致する部分はそのままで、一致しない部分はその画素値を“０”に反転するような処理である。このようなパターンマッチングにより、図５（ｆ）に示すような、ノイズの大部分が除去され、線状陰影が存在する確率が比較的高いと推定される部分だけを抽出することができる。上述したＡ系列で得られた２値画像を、“Ｆ１”とする。
【００２９】
次に、Ｂ系列について説明する。Ｂ系列は、Ａ系列と同様に、処理順に、Ｂ１，Ｂ２，Ｂ３，Ｂ４，Ｂ５の５つのステップからなる。まず、ステップＢ１では、ディジタルＸ線画像Ｉ0 に対して、空間的な平滑化フィルタリングという処理を施して、平均値画像Ｉave を得る。平滑化フィルタリングは、周知の通り、図６に示すように、例えば２×２マトリクスのフィルタマスクを重ならないように移動し、その各位置で、マスク内の４つの近傍画素の画素値の平均値を計算し、その平均値を平均値画像Ｉave の中の対応する画素の画素値とするいう処理である。
【００３０】
なお、この処理により、マトリクスサイズを１／４に縮小する、例えばディジタルＸ線画像Ｉ0 のマトリクスサイズが１０２４×１０２４の大きさである場合には、５１２×５１２のマトリクスサイズに縮小されるが、これ以降の処理を５１２×５１２のマトリクスサイズで行うことにより回路の規模を縮小することができる。
このような平滑化フィルタリングによると、背景ノイズが低減されるので、その分、線状陰影が認識しやすくなる。
【００３１】
次にステップＢ２では、Ｂ１で得られた平均値画像Ｉave に対して、ローパスタイプ（低域通過型）の空間フィルタリングを施し、空間周波数が比較的低い成分が濾過された平均値画像Ｉlow を得る。
【００３２】
このステップＢ２でローパスフィルタされた平均値画像Ｉlow と、先のステップＢ１で得られた平均値画像Ｉave に対して、ステップＢ３で、次の式（２）ような線形演算、つまり加重サブトラクション（加重を掛けて差分を取る）を実行することにより、差分画像Ｉsub2を得る。
【００３３】
Ｉsub2＝Ｗ3 ×Ｉave −Ｗ4 ×Ｉlow ＋Ｄ    …（２）
図７（ａ）に、ローパスフィルタされた平均値画像Ｉlow と、平均値画像Ｉave それぞれの空間的なプロフィールを示し、同図（ｂ）に、上記式（２）の差分処理により得られた差分画像Ｉsub2の空間的なプロフィールを示している。上記式（２）の係数Ｗ3 、Ｗ4 、Ｄを適当に設定することにより、図７（ｂ）に示すように、差分画像Ｉsub2上で、線状陰影である確率が比較的高いと推定される部分だけを、それ以外の部分に対して、極性を反転させることができる。なお、図７（ｂ）の例では、４つの領域が線状陰影である可能性が比較的高いと推定されているが、この内、左から２つ目の領域が線状陰影を正しく識別した部分で、残りの３つの領域は、誤識別した部分である。
【００３４】
このように線状陰影である確率が比較的高いと推定される部分と、それ以外の部分とが極性で識別された差分画像Ｉsub2を、その極性に従って２値化し、線状陰影である確率が比較的高いと推定される部分が“１”で、またそれ以外の部分が“０”で表現された２値画像Ｉbi2 を得る（ステップＢ４）。
【００３５】
この２値画像Ｉbi2 にも、殆どの場合で背景ノイズが残存しているので、これを取り除く処理が必要になり、Ａ５と同様に２値画像Ｉbi2 に対してパターンマッチングを実行し、（ステップＡ５）。ノイズの大部分が除去され、線状陰影が存在する確率が比較的高いと推定される部分だけを抽出することができる。上述したＢ系列で得られた２値画像を、“Ｆ２”とする。
【００３６】
このように２系統の画像処理で得られた２種類の２値画像Ｆ１，Ｆ２から、ステップＡＢで、４値のフラグ画像Ｆ３を作成する。このフラグ画像Ｆ３の画素値は、同じ位置の画素に関する一方の２値画像Ｆ１の画素値と、他方の２値画像Ｆ２の画素値との組み合わせに従って、次の表１に示すように、割り当てられる。
【００３７】
【表１】
0004095140.tif 000002
【００３８】
２種類の２値画像Ｆ１，Ｆ２では両方共に、線状陰影である確率が比較的高いと推定される部分が“１”で、またそれ以外の部分が“０”で表現されているので、２種類の２値画像Ｆ１，Ｆ２の両方で画素値が“１”である場合、その画素は、線状陰影である確率が最も高いと考えられ、ここでは、フラグ画像Ｆ３の該当画素に“３”の画素値を与えるようにしている。
【００３９】
逆に、２種類の２値画像Ｆ１，Ｆ２の両方で画素値が“０”である場合、その画素は、線状陰影である確率が最も低いと考えられ、ここでは、フラグ画像Ｆ３の該当画素に“０”の画素値を与えるようにしている。
【００４０】
また、最小値フィルタリングの方が、平滑化フィルタリングよりも、線状陰影の抽出能が一般的に高いので、最小値フィルタリングを通した２値画像Ｆ１では“１”が与えられているが、平滑化フィルタリングを通した２値画像Ｆ２では“０”である場合、線状陰影である確率が２番目に高いとして、ここでは、フラグ画像Ｆ３の該当画素に“２”の画素値を与え、一方、最小値フィルタリングを通した２値画像Ｆ１では“０”が与えられているが、平滑化フィルタリングを通した２値画像Ｆ２では“１”である場合、線状陰影である確率が３番目に高いとして、ここでは、フラグ画像Ｆ３の該当画素に“１”の画素値を与えるようにしている。
【００４１】
次に、リカーシブフィルタ（再帰型フィルタ）と呼ばれる時間フィルタに、ディジタルＸ線画像Ｉ0 又はそれから画像処理により導かれる画像を通すのであるが、本実施形態では、このリカーシブフィルタを、背景ノイズ用と線状陰影用との２系統設けている（Ｃ２，Ｄ２）。リカーシブフィルタは、周知の通り、前回のフィルタ結果を、今回のフィルタリングに反映させるというタイプの時間フィルタであり、以下の式（３）で定義される。なお、“ａ”は、前回のフィルタ結果を、今回のフィルタリングにどの程度反映させるかというフィルタ特性を決定するリカーシブ係数であり、ここでは“ａ”が低いほど、前回のフィルタ結果を、今回のフィルタリングに強く反映させることになる。
【００４２】
Ｙi(x,y) ＝ａ×Ｘi(x,y) ＋（１−ａ）×Ｙi-1(x,y)
（ｉ＝１，２，…，∞）
ここで、Ｙiはｉ番目にリカーシブフィルタリングされた画像を表し、Ｘiは、ｉ番目にＴＶカメラ７から画像処理装置９に供給された現フレームのディジタルＸ線画像Ｉ0 又は、それから画像処理により導かれる画像、例えばＩ0 ，Ｉmin ，Ｉave の３画像を使って、次式（４）で与えられるフレーム間加重平均画像を表している。なお、Ｉmin 、Ｉave は、Ｉ0 と等しいマトリクスサイズに画素数を変換してから（４）式に用いる。また、係数α、β、γは、固定の係数としてもよいし、フラグ画像Ｆ３の値に応じて画素毎に変えるような構成としてもよい。
【００４３】
Ｘi(x,y) ＝α×Ｉ0 (x,y) ＋β×Ｉmin (x,y) ＋γ×Ｉave (x,y)   …（４）
（α＋β＋γ＝１）
この２系統のリカーシブル係数ａを、フラグ画像Ｆ３の画素値に応じて画素ごとに変えていくのであるが、その変え方を背景用リカーシブフィルタと線状陰影用リカーシブフィルタとで相違させている。これは、両フィルタの目的の違いによるものであり、背景用リカーシブフィルタは、背景ノイズを低減することを目的としているため、リカーシブル係数ａを比較的低い範囲の中から割り当てて前回のフィルタ結果をより強く反映させ（Ｃ２）、逆に、線状陰影用リカーシブフィルタは、線状陰影を強調することを目的とし、前回のフィルタ結果があまり強く反映されてしまうと、カテーテルなどの挿入時の動きによる“ぼけ”が生じてしまいかねないので、リカーシブル係数ａを比較的高く設定し、前回のフィルタ結果をあまり強く反映させないようにしている（Ｄ２）。
【００４４】
具体的には、上述した表１に示したように、背景用リカーシブフィルタのリカーシブル係数ａを、フラグが３のときには“０．３”、フラグが２のときには“０．７”、フラグが１のときには“０．６”、フラグが０のときには“０．５”を割り当てる。背景用リカーシブフィルタは、基本的には線状陰影の存在確率が高い部分ほどリカーシブ係数を高くして残像の影響が少なくなるように構成されているが、後述のステップＥ１，Ｅ２で線状陰影用リカーシブフィルタの出力画像Ｉ1 の画素の値が出力画像形成回路の出力画像Ｉ3 として出力される部分（本実施形態では、フラグが３の部分）については極めて小さいリカーシブ係数を割り当てて、ぼけた画像となるように構成されている。
【００４５】
また、線状陰影用リカーシブフィルタのリカーシブル係数ａを、フラグが３のときには“０．９”、それ以外のときには全て“０．８”を割り当てる。線状陰影用リカーシブフィルタは、後述のステップＥ１，Ｅ２で線状陰影用リカーシブフィルタの出力画像Ｉ1 の画素の値が出力画像形成回路の出力画像Ｉ3 として出力される部分（本実施形態では、フラグが３の部分）については極めて大きいリカーシブ係数を割り当てて、残像の影響が小さくなるように構成されている。
【００４６】
このようなリカーシブフィルタリング処理により、背景ノイズが非常に低減された画像Ｉ1 と、線状陰影がその動きによるぼけが抑えられて十分強調されている画像Ｉ2 との２種類の画像が得られる。なお、本実施形態では、線状陰影を表す画像Ｉ2 として、リカーシブフィルタの出力を用いているが、リカーシブフィルタの出力画像の代わりに、原画像Ｉ0 、最小値画像Ｉmin 、平均値画像Ｉave の何れかを用いるような構成としてもよい。
【００４７】
次にステップＥ１，Ｅ２で、このような２種類の画像Ｉ1 ，Ｉ2 のいずれかを、フラグ画像Ｆ３の線状陰影の存在確率を表しているフラグに従って、画素ごとに選択する。具体的には、表１に示すように、フラグが０，１，２の部分には、背景用リカーシブフィルタの出力画像Ｉ1 の画素値を画像Ｉ3 の画素値として選択し、フラグが３の部分には、線状陰影用リカーシブフィルタの出力画像Ｉ2 の画素値を画像Ｉ3 の画素値として選択するように構成されている。このとき、出力される画像Ｉ3 は、背景ノイズが十分低減されると共に、線状陰影は十分なコントラストで明瞭な画像となっており、また画像Ｉ0 中でカテーテルやガイドワイヤ等が移動した場合にもそれらの残像を殆ど含まない画像である。
【００４８】
なお、空間フィルタ係数とフィルタサイズ、パターンマッチドフィルタとフィルタサイズ、リカーシブ係数等は、認識の対象とする線状陰影の空間的太さと濃度、背景ノイズの周波数特性、Ｉ．Ｉ．視野サイズ、及び透視Ｘ線条件を考慮して適切に定めることにより、最適な処理を行うことができる。
【００４９】
また、一般に、ハードウェアの回路規模は同時に多数の平均化回路や比較回路が動作することや、画像上の位置指定信号である画素アドレスの発生が煩雑である場合などに大きくなり実用的なリアルタイム演算の可能性が小さくなる。従来技術に見るようなパターン認識は、リアルタイム認識に適さない構成であるが、本実施形態では線状陰影のパターン認識は付随情報のビット数と同じ個数の線形空間フィルタとその後段に位置する比較器、及び２値化画像に対して施されるバイナリーフィルタで構成されるため順次走査による処理が可能であり回路規模が小さくてすむ。
【００５０】
さらに、本実施形態では、２系統の画像処理を用いて線状陰影の存在確率を決定したが、２系統以上の画像処理の選別結果に基づいて線状陰影の存在確率を決定する構成としてもよい。
【００５１】
本発明は、上述した実施形態に限定されることなく、種々変形して実施可能である。例えば、上述の説明では、Ｘ線透視画像を対象としているが、ＤＳＡなどの撮影画像やワンショットのＸ線撮影ではリカーシブフィルタを構成することができないケースも多い。なぜならば、連続的に画像が入力されず、間欠的、あるいは１枚だけの画像入力となるからである。この場合には、リカーシブフィルタの代わりに空間的フィルタを使用することができる。画像Ｉ1 は比較的強いスムージングが施され、Ｉ2 はコントラスト強調された画像、もしくは、最小値画像を用いると線状陰影のコントラストは増強されるかまたは損なわれることがなく、背景は平滑化される。この場合でも背景用スムージングはＦ１＝０，Ｆ２＝０の付随情報がある画素はスムージング効果を大きくし、Ｆ１＝１，Ｆ２＝０、もしくはＦ１＝０，Ｆ２＝１の付随情報の画素はさほど大きくない設定にして背景画像に存在する可能性のある線状陰影のコントラスト劣化を押さえるようにする必要がある。このような設計に基づいて画像を形成することによって血管や胃壁などのコントラストが保たれたまま背景のざらつきが軽減された画像を観察者に提供することができる。
【００５２】
【発明の効果】
本発明は、線状陰影及びノイズが含まれている画像に対して少なくとも２種類の空間的な画像処理を施し、この２種類の空間的な画像処理の結果に基づいて線状陰影の存在確率を表す確率分布を求めるようにしたので、線状陰影の存在の有無を比較的高精度で判定することができる。この確率分布に従って、時間的な画像処理の特性を画素又は近隣グループごとに変えることにより、ノイズを効果的に低減しながら、それと同時に、Ｘ線医療におけるガイドワイヤ、カテーテル、血管影などの線状陰影のコントラストを効果的に強調できる。
【図面の簡単な説明】
【図１】本発明の実施形態に係る画像処理装置が組み込まれたＸ線透視装置の構成図。
【図２】本発明の実施形態に係る画像処理装置の処理手順を示すフローチャート。
【図３】図２の最小値フィルタリングの補足説明図。
【図４】図２の最小値画像Ｉmin 、ハイパス空間フィルタを受けた最小値画像Ｉhigh、これら両者のサブトラクション画像Ｉsub1それぞれのプロフィールを示す図。
【図５】図２のパターンマッチングの補足説明図。
【図６】  図２の平滑化フィルタリングの補足図。
【図７】図２の平均値画像Ｉave 、ローパス空間フィルタを受けた平均値画像Ｉlow 、これら両者のサブトラクション画像Ｉsub2それぞれのプロフィールを示す図。
【符号の説明】
１…Ｘ線透視装置本体、
２…Ｘ線源、
３…Ｘ線、
３ａ…透過Ｘ線、
４…被検体、
４ａ…胸部、
５…寝台、
６…イメージインテンシファイア、
７…ＴＶカメラ、
８…ディジタルＸ線画像信号、
９…画像処理装置、
１０…ＣＲＴディスプレイ。
【図１】
0004095140.tif 000003
【図２】
0004095140.tif 000004
【図３】
0004095140.tif 000005
【図４】
0004095140.tif 000006
【図５】
0004095140.tif 000007
【図６】
0004095140.tif 000008
【図７】
0004095140.tif 000009