JP4639526B2 - ディジタルサブトラクション装置 - Google Patents

Description

【０００１】
【発明の属する技術分野】
この発明は、医療分野、工業分野などに用いられる、被検体の撮像部位をＸ線透視撮像してその撮像部位のサブトラクション像を得るディジタルサブトラクション装置に係り、特に、被検体の撮像部位のサブトラクション像を好適に高速に得る技術に関する。
【０００２】
【従来技術】
従来のディジタルサブトラクション装置としては、例えば、医療分野で用いられる、被検体の所定の撮像部位のサブトラクション像を得るディジタルアンギオグラフィ装置がある。このディジタルアンギオグラフィ装置としては、例えば、特開平8-308823号公報に記載されたもののように、被検体の撮像部位を１回Ｘ線透視撮像するだけでその撮像部位のサブトラクション像を得れるものがある。以下に、このディジタルアンギオグラフィ装置により、被検体の所定の撮像部位のサブトラクション像を得る動作について、説明する。
【０００３】
まず、造影剤が投与された被検体の所定の撮像部位を、Ｘ線透視撮像装置（Ｘ線管と、イメージインテンシファイアとテレビカメラなどで構成される撮像系とを備えたもの）でもってＸ線透過像として撮像し、この撮像したＸ線透視像をディジタルデータに変換してＸ線透過像を取得する。このＸ線透過像は、造影剤が投与された血管像などの高周波数成分が残っている画像であり、これをライブ像として用いる。一方、周波数特性変換回路によって、Ｘ線透過像（ライブ像）を構成する周波数成分を空間／周波数変換処理で取り出し、所定のしきい値周波数以上の周波数成分を除去し、これに周波数／空間変換処理を施して、前記Ｘ線透過像（ライブ像）から前記しきい値周波数以上の周波数成分を除去した高周波数成分除去画像を得ている。この高周波数成分除去画像は、造影剤が投与された血管像などの高周波数成分が除去された画像であり、これをマスク像として用いる。前記の空間／周波数変換処理としては、ＦＦＴ（高速フーリエ変換）、カルーネン・レーベ変換、ＤＣＴ（離散コサイン変換）、アダマール変換などの各変換方式がある。また、前記の周波数／空間変換処理としては、前記空間／周波数変換処理の逆変換（逆ＦＦＴ、逆カルーネン・レーベ変換、逆ＤＣＴ、逆アダマール変換など）がある。
【０００４】
次に、遅延回路によって、後段の演算器へのＸ線透過像（ライブ像）の供給を前記周波数特性変換回路の処理時間分だけ遅延させることで、Ｘ線透過像（ライブ像）と高周波数成分除去画像（マスク像）とを後段の演算器に同期して供給している。そして、演算器は、Ｘ線透過像（ライブ像）から高周波数成分除去画像（マスク像）をサブトラクションしてサブトラクション像（サブトラクション画像）を求めている。なお、上述の所定のしきい値周波数は、サブトラクション像に残したい関心物（血管像など）を好適に除去し得る周波数値のことであり、このしきい値周波数を予め理論的あるいは実験的に求めておいて設定することで、サブトラクション像を得るようにしている。
【０００５】
このように、造影剤が投与された被検体の所定の撮像部位を１回だけＸ線透視撮像するだけで、その撮像部位のサブトラクション像を得ることができるので、造影剤の投与前後の２回にわたってＸ線透視撮像する場合に比べて、被検体へのＸ線曝射線量が軽減できるとともに、Ｘ線透過像（ライブ像）から高周波数成分除去画像（マスク像）を生成していることから、造影剤投与前後の２回の撮像における被検体の体動などに起因するマスク像とライブ像の画像ずれを完全に無くすことができる。
【０００６】
【発明が解決しようとする課題】
しかしながら、このような構成を有する従来例の場合には、次のような問題がある。
すなわち、上述した従来例では、撮影して得られたＸ線透過像（ライブ像）から所定のしきい値周波数以上の周波数成分を除去して高周波数成分除去画像（マスク像）を生成するという周波数特性変換処理に一定の時間がかかるため、Ｘ線透過像（ライブ像）を撮影した後、このＸ線透過像（ライブ像）から高周波数成分除去画像（マスク像）をサブトラクションすることで得られるサブトラクション像が生成されるまでには一定時間の遅延が生じ、被検体の動きに対して、サブトラクション像表示の追従性が悪いという問題がある。
【０００７】
また、撮影して得られたＸ線透過像（ライブ像）から所定のしきい値周波数以上の周波数成分を除去して高周波数成分除去画像（マスク像）を生成するという周波数特性変換処理を、撮影して得られたＸ線透過像（ライブ像）ごとに行わなければならないので、撮影レートをこれ以上に向上させることができないという問題もある。
【０００８】
この発明は、このような事情に鑑みてなされたものであって、被検体の撮像部位のサブトラクション像を好適に高速に得ることができるディジタルサブトラクション装置を提供することを目的とする。
【０００９】
【課題を解決するための手段】
この発明は、このような目的を達成するために、次のような構成をとる。
すなわち、請求項１に記載のディジタルサブトラクション装置は、被検体の所定の撮像部位のサブトラクション像を得るためのディジタルサブトラクション装置であって、（ａ）前記撮像部位にＸ線を照射し、その部位の複数枚のＸ線透過像を撮像するＸ線透視手段と、（ｂ）前記Ｘ線透過像をディジタルデータに変換するデータ変換手段と、（ｃ）前記ディジタルデータに変換されたＸ線透過像ごとに所定のしきい値周波数以上の周波数成分を除去した画像（以下、高周波数成分除去画像）を得る周波数特性変換手段と、（ｄ）Ｘ線透過像である各基本画像と各高周波数成分除去画像とのサブトラクションを行い、その撮像部位のサブトラクション像を求める演算手段と、（ｅ）前記高周波数成分除去画像と、この高周波数成分除去画像を生成するためのＸ線透視像の所定回後に別途撮影されたＸ線透過像である基本画像とを前記演算手段に同期して供給する供給制御手段とを備えたことを特徴とするものである。
【００１０】
また、請求項２に記載のディジタルサブトラクション装置は、請求項１に記載のディジタルサブトラクション装置において、前記供給制御手段は、前記高周波数成分除去画像と、別途撮影されたＸ線透過像としての、この高周波数成分除去画像が生成された時刻より以降に撮影されたＸ線透過像である基本画像とを前記演算手段に同期して供給することを特徴とするものである。
【００１１】
また、請求項３に記載のディジタルサブトラクション装置は、請求項１または請求項２に記載のディジタルサブトラクション装置において、前記供給制御手段は、前記高周波数成分除去画像と、別途撮影されたＸ線透過像としての、この高周波数成分除去画像が生成されるまでに撮影されたＸ線透過像である基本画像とを前記演算手段に同期して供給することを特徴とするものである。
【００１２】
【作用】
この発明の作用は次の通りである。
すなわち、請求項１に記載の発明によれば、被検体の所定の撮像部位に対して、Ｘ線透視手段で複数枚のＸ線透過像を撮像し、その撮像したＸ線透過像をデータ変換手段でディジタルデータに変換したＸ線透過像を得る。このＸ線透過像は、関心物の高周波数成分が残っている画像であり、これをライブ像として用いる。一方、周波数特性変換手段は、前記Ｘ線透過像ごとに所定のしきい値周波数以上の周波数成分を除去した高周波数成分除去画像を得る。この高周波数成分除去画像は、関心物の高周波数成分が除去された画像であり、これをマスク像として用いる。そして、供給制御手段は、高周波数成分除去画像（マスク像）と、この高周波数成分除去画像を生成するためのＸ線透視像の所定回後に別途撮影されたＸ線透過像である基本画像（ライブ像）とを演算手段に同期して供給する。演算手段は、基本画像（ライブ像）から高周波数成分除去画像（マスク像）をサブトラクションしてサブトラクション像を求める。
【００１３】
したがって、高周波数成分除去画像（マスク像）の生成元であるＸ線透過像（ライブ像）ではなく、その高周波数成分除去画像（マスク像）が生成されたＸ線透過像とは別途に撮影されたＸ線透過像である基本画像（ライブ像）から、高周波数成分除去画像（マスク像）をサブトラクションしているので、周波数特性変換処理の所要時間分の遅延を含まない基本画像（ライブ像）をサブトラクションの対象とすることができ、除去されるべき低周波数成分の時間追従性は悪くなるものの、関心物である高周波数成分は別途に撮影された基本画像（ライブ像）に含まれているため、最小限の遅延でサブトラクション像などが表示される。
【００１４】
また、請求項２に記載の発明によれば、請求項１に記載の発明と同様の作用でＸ線透過像（ライブ像）と高周波数成分除去画像（マスク像）とが得られる。そして、供給制御手段は、高周波数成分除去画像（マスク像）と、別途撮影されたＸ線透過像としての、この高周波数成分除去画像（マスク像）が生成された時刻より以降に撮影されたＸ線透過像である基本画像（ライブ像）とを演算手段に同期して供給する。演算手段は、基本画像（ライブ像）から高周波数成分除去画像（マスク像）をサブトラクションしてサブトラクション像を求める。したがって、高周波数成分除去画像（マスク像）の生成元であるＸ線透過像（ライブ像）ではなく、その高周波数成分除去画像（マスク像）が生成された時刻以降のＸ線透過像である基本画像（ライブ像）から、高周波数成分除去画像（マスク像）をサブトラクションしているので、周波数特性変換処理の所要時間分の遅延を含まない基本画像（ライブ像）をサブトラクションの対象とすることができ、除去されるべき低周波数成分の時間追従性は悪くなるものの、関心物である高周波数成分は高周波数成分除去画像（マスク像）が生成された時刻以降の基本画像（ライブ像）に含まれているため、最小限の遅延でサブトラクション像などが表示される。
【００１５】
また、請求項３に記載の発明によれば、請求項１に記載の発明と同様の作用でＸ線透過像（ライブ像）と高周波数成分除去画像（マスク像）とが得られる。そして、供給制御手段は、高周波数成分除去画像（マスク像）と、別途撮影されたＸ線透過像としての、この高周波数成分除去画像（マスク像）が生成されるまでに撮影されたＸ線透過像（ライブ像）とを演算手段に同期して供給する。したがって、複数枚のＸ線透過像（ライブ像）に対して個別に同じ高周波数成分除去画像（マスク像）でサブトラクションすることになり、撮影して得られたＸ線透過像（ライブ像）ごとに高周波数成分除去画像（マスク像）を生成する必要が無いので、その時間分周波数変換処理時間を低減することができ、撮影レートを向上させる、すなわち、表示間隔を短くすることができる。
【００１６】
【発明の実施の形態】
以下、図面を参照してこの発明のディジタルサブトラクション装置に係る一実施例としてのＸ線ディジタルアンギオグラフィ装置について、図面を参照しながら説明する。
【００１７】
＜第１実施例＞
図１は、この発明の第１実施例に係るＸ線ディジタルアンギオグラフィ装置の全体構成を示す正面図であり、図２は、Ｘ線透視装置を側面から見た図であり、図３は、画像処理部の構成を示すブロック図である。
【００１８】
この第１実施例のＸ線ディジタルアンギオグラフィ装置は、ベッド１、Ｘ線透視手段としてのＸ線透視装置２、画像処理部３、モニタ４、制御部５、操作盤６などを備えて構成されている。
【００１９】
ベッド１は、床面に設置されたベッド基台１１と天板１２を備えている。被検体Ｍは天板１２上に載置される。この天板１２はモータ１３の駆動で水平移動可能であり、天板１２上の被検体ＭとＸ線透視装置２との相対的な位置関係を被検体Ｍの体軸方向に変位することができる。モータ１３の駆動制御は、制御部５により行われる。
【００２０】
Ｘ線透視装置２は、Ｘ線管２１、撮像系２２を支持するＣ型アーム２３がベッド１の近傍に定置された装置基台２４の上部に支持されて構成されている。Ｃ型アーム２３は、モータ２５の駆動で図２の矢印方向に変位可能に装置基台２４に支持されており、Ｘ線管２１、撮像系２２を被検体Ｍの体軸回りに変位可能に構成し、Ｘ線透過像の撮像方向の調整が可能となっている。モータ２５の駆動制御は制御部５により行われる。
【００２１】
Ｘ線管２１と撮像系２２とはＣ型アーム２３の両端部に取り付けられており、天板１２上の被検体Ｍを挟み込んだ状態で対向配置されている。Ｘ線管２１から被検体Ｍの任意の撮像部位に向けて照射され、被検体Ｍを透過したＸ線は、撮像系２２で受像され、その部位のＸ線透過像が撮像される。Ｘ線管２１からのＸ線の照射は、Ｘ線高電圧発生装置２６から所定の電力（Ｘ線管電圧およびＸ線管電流）がＸ線管２１に供給されて行われる。Ｘ線高電圧発生装置２６からＸ線管２１への所定の電力の供給は制御部５に制御されて行われる。撮像系２２は、イメージインテンシファイアとテレビカメラ、ＦＰＤ（Ｆｌａｔ Ｐａｎｅｌ Ｄｅｔｅｃｔｏｒ）などで構成されている。撮像されたＸ線透過像は画像処理部３に与えられる。
【００２２】
画像処理部３は、図３に示すように、データ変換手段としてのＡ／Ｄ（アナログtoディジタルデータ）変換器３１と、高周波除去手段としての周波数特性変換回路３２を備えた供給制御手段としての供給制御回路３３と、演算手段としての演算器３４と、階調変換回路３５と、Ｄ／Ａ（ディジタルtoアナログ）変換器３６とで構成されている。
【００２３】
造影剤が投与された被検体Ｍの所定の撮像部位のＸ線透過像が撮像されると、撮像系２２からのその像の画像信号（アナログ信号）は、Ａ／Ｄ変換器３１でディジタルデータに変換されて、Ｘ線透過像が得られる。このＸ線透過像は、骨格などの低周波数成分や、造影剤が投与された血管像などの高周波数成分を含んだ画像であり、これをライブ像として用いる。このＸ線透過像（ライブ像）は、周波数特性変換回路３２に与えられる。周波数特性変換回路３２は、後述する処理によってＸ線透過像（ライブ像）から血管像などの高周波数成分を除去しマスク像を得る。供給制御回路３３は、高周波数成分除去画像（マスク像）と、別途撮影されたＸ線透過像としての、この高周波数成分除去画像（マスク像）が生成された時刻より以降に撮影されたＸ線透過像である基本画像（ライブ像）とを演算器３４に同期して供給する。この供給制御回路３３の機能がこの発明の１つの特徴である。
【００２４】
そして、演算器３４では、直接供給される基本画像（ライブ像）と、周波数特性変換回路３２を経て供給される高周波数成分除去画像（マスク像）とのサブトラクションを行いサブトラクション像を求めて階調変換回路３５に与える。階調変換回路３５では、サブトラクション像をモニタ５に表示したとき見やすい画像にするために、サブトラクション像を構成する各画素の濃度を調整（全画素を対象に、各画素の濃度に所定濃度を加算したり減算する）する。階調変換されたサブトラクション像はＤ／Ａ変換器３６に与えられ、そこでＤ／Ａ変換されてモニタ４に表示される。なお、画像処理部３を構成する各部の動作制御は、制御部５により行われる。
【００２５】
この第１実施例の周波数特性変換回路３２は、Ｘ線透過像（ライブ像）から血管像などの高周波数成分を除去して高周波数成分除去画像（マスク像）を得ることを、実空間データのままで高速処理することを目的とするものである。
【００２６】
周波数特性変換回路３２は、実空間上で直線形状で表されるテンプレートフィルタでもってＸ線透過像（ライブ像）の高周波数成分を実空間上で除去する実空間フィルタリング部３２ａを備えている。この実空間フィルタリング部３２ａのテンプレートフィルタ形状は、例えば、図７に示すように、互いに直交するＸ，Ｙ，Ｚ軸で表される、実空間上の３次元直交座標系において、Ｘ，Ｙ軸方向をテンプレートフィルタの矩形形状である底面大きさとし、Ｚ軸方向をこのテンプレートのゲインとした四角柱形状で表されるものとしている。この実空間上で四角柱形状のテンプレートフィルタ（以下、適宜に「ボックスフィルタ」と呼ぶ）による処理方法は、図８に示すように、入力画像（入力される基本画像）における画素（ｉ，ｊ）の近傍（Ｎ×Ｎ点の正方形分）の平均濃度を出力画像の画素（ｉ，ｊ）の値とし、この平均処理を全画素について行うことで、Ｘ線透過像（ライブ像）から所定のしきい値周波数以上の周波数成分（高周波数成分）を除去した高周波数成分除去画像（マスク像）を得るものであり、移動平均フィルタ法と適宜に呼ぶこととする。
【００２７】
なお、図８は、図７に示したボックスフィルタを、フィルタ処理しようとするＸ線透過像（ライブ像）の所定の複数個の画素上に位置させた状態をＺ軸方向から見下ろした図である。図８では、ボックスフィルタのフィルタサイズを、説明の便宜上、３×３の正方形として図示し、画素（ｉ，ｊ）とその周囲の８画素とからなる９画素（９点）の平均値を画素（ｉ，ｊ）の値としているが、このフィルタサイズは、次に説明するように所定のしきい値周波数以上の周波数成分（高周波数成分）を除去する大きさに設定される。すなわち、このボックスフィルタは、そのフィルタサイズに応じてしきい値周波数が変更されることになる。つまり、フィルタサイズを大きくすればしきい値周波数を下げることになりぼかしの程度が大きくなるし、フィルタサイズを小さくすればしきい値周波数を上げることになりぼかしの程度が小さくなるのである。ボックスフィルタのフィルタサイズは、図７に示すように、例えば２１×２１点（ＴＡＰ数×ＴＡＰ数）の正方形としている。このＴＡＰ数はカーネルの点数のことでもある。また、標準的な空間的形状の血管像を好適に除去し得るしきい値周波数は、実験的あるいは理論的に求めておくことができるので、標準的な空間的形状の血管像を好適に除去し得るボックスフィルタのフィルタサイズも同様に実験的あるいは理論的に求めておいても良いし、そのフィルタサイズは初期値として設定しておいても良い。
【００２８】
なお、上述の移動平均フィルタ法は、次に説明するようにすることで、高速に演算処理ができるように改良されたものである。例えば、フィルタサイズがＮ点×Ｎ点であるボックスフィルタを用いて移動平均を求めるには、ボックスフィルタが位置するＮ×Ｎ点の画素を正直に加算した合計値をボックスフィルタの点数分（Ｎ×Ｎ点）で平均することで、ボックスフィルタの中心の画素の値を算出し、この演算を全画素について個別に行なっていたのでは、演算量が膨大になり高速に演算処理することができないし、フィルタサイズが大きくなるにつれてその演算量は増加することになる。
【００２９】
そこで、次に説明するようにして、上述の移動平均フィルタ法を高速に演算処理可能に改良している。すなわち、図９に示すように、Ｎ×Ｎ点の移動平均を求めるには、各画素についてＮ×Ｎ点の加算を正直に計算することなく、前回の結果に対して、更新される部分についての加減算を行うことで、簡単に計算でき、演算時間が短縮できる。例えば、１次元状の「Ａ」〜「Ｆ」の画素に対して、１次元状の３点のボックスフィルタで移動平均する場合を例に挙げて説明する。「Ａ」〜「Ｃ」の画素にボックスフィルタがある場合は、「Ａ＋Ｂ＋Ｃ」の３点平均値が「Ｂ」画素の値となる。そして、ボックスフィルタを次の位置（「Ｂ」〜「Ｄ」の画素の位置）に移動させたときには、前回の値（「Ａ＋Ｂ＋Ｃ」）から更新される部分（「Ｄ」が加わり、「Ａ」が外される）についての加減算（−Ａ＋Ｄ）を行うだけでよく、前回と今回とで重複する加算を行なうという無駄が排除できるし、フィルタサイズが大きくなってもその演算量（演算の総量）は低減したまま一定である。この例では、説明の便宜上、ボックスフィルタを３点としているが、ボックスフィルタの点数が大きくなればなる程、全点について正直に順次加算していく演算量と、更新される部分の加減算のみを行なう演算量との差は大きくなり、効果的に改善されることがわかる。なお、図７，図８に示すようにＮ×Ｎ点のボックスフィルタの２次元移動平均の場合には、縦方向のＮ行累算値を、図９と同様の方法で求め、さらに横方向に図９と同様の方法を繰り返していけばよい。すなわち、Ｎ×Ｎ点のボックスフィルタの２次元移動平均の場合であっても、更新される部分の加減算は、１画素の加算および１画素の減算のみである。上述したように演算の総量を低減しているので、汎用のチップを用いて実空間フィルタリング部３２ａを構成することができ、撮像して得られたＸ線透過像（ライブ像）から極めて短時間で（リアルタイムに近く）高周波数成分除去画像を生成することができる。
【００３０】
図１，図２に戻って、制御部５は、操作盤６からの各種の指示などによって、各装置、各部の駆動制御や動作制御を行う。この制御部５は、例えば、後述する動作を実現するプログラムを遂行するＣＰＵ（中央処理装置）で構成されている。
【００３１】
操作盤６は、撮像部位や条件の設定、処理開始指示などを、操作者が行うためのものである。
【００３２】
上記構成を有する実施例装置の動作を以下に説明する。
まず、被検体Ｍのある１箇所の撮像部位（例えば胸部）のサブトラクション像を得る場合の動作を説明する。
【００３３】
この場合、まず、操作者により操作盤６から設定された撮像部位や条件（撮像方向など）に従って、制御部５は、モータ１３を駆動制御して被検体Ｍを載置した天板１２を水平移動させ、設定された撮像部位（胸部とする）を、Ｘ線管２１、撮像系２２の間の撮像位置に位置させ、モータ２５を駆動制御してＸ線管２１、撮像系２２を被検体Ｍ（の撮像部位）の体軸回りに変位させ、撮像方向を調節する。この状態を図４に示す。図４では、撮像部位ＳＢの下方からＸ線を照射してＸ線透過像を撮像するように撮像方向が調節されている。
【００３４】
次に、被検体Ｍに造影剤を投与する。なお、造影剤を投与してから上記位置合わせ動作などを行ってもよい。いずれにしても、以下の撮像動作の前に、被検体Ｍに造影剤を投与しておき、撮像部位ＳＢに造影剤が拡散した状態で操作者が操作盤６から処理開始を指示し、以下の撮像動作が実行される。
【００３５】
処理開始が指示されると、制御部５はＸ線高電圧発生装置２６を制御して、Ｘ線管２１に所定の電力を供給させてＸ線を照射させ、造影剤が拡散された撮像部位ＳＢのＸ線透過像を撮像させる。そして、制御部５は、画像処理部３の各部を制御して、Ｘ線透過像（ライブ像）を得るとともに、そのＸ線透過像（ライブ像）から高周波数成分除去画像（マスク像）を求め、この求めたマスク像と、このマスク像が生成された時刻より以降に撮影されたＸ線透過像である基本画像（ライブ像）とを演算器３４に同期して供給して、演算器３４でマスク像と基本画像（ライブ像）とのサブトラクションを行わせ、サブトラクション像をモニタ４に表示させる。
【００３６】
具体的には、供給制御回路３３は、図１０（ｂ）に示すように、撮像系２２で撮像されてＡ／Ｄ変換器３１でディジタルデータに変換されたＸ線透過像（ライブ像）「Ａ」から高周波数成分除去画像（マスク像）「Ａ´」を周波数特性変換回路３２により求め、この求めたマスク像「Ａ´」と、このマスク像「Ａ´」が生成された時刻ｔ２より以降に撮影されたＸ線透過像である基本画像（ライブ像）例えば「Ｂ」とを演算器３４に同期して供給する。演算器３４は、基本画像（ライブ像）「Ｂ」からマスク像「Ａ´」をサブトラクションしてサブトラクション像「Ｂ−Ａ´」を生成する。なお、次のサブトラクション像は、Ｘ線透過像（ライブ像）「Ｂ」から求めた高周波数成分除去画像「Ｂ´」と、このマスク像「Ｂ´」が生成された時刻ｔ６より以降に撮影されたＸ線透過像である基本画像（ライブ像）「Ｃ」とを演算器３４に同期して供給することで、マスク像「Ｂ´」と基本画像（ライブ像）「Ｃ」とによるサブトラクション像「Ｃ−Ｂ´」を生成している。あるＸ線透過像（ライブ像）から次のＸ線透過像（ライブ像）を撮影するまで時間（例えば、時刻ｔ２〜ｔ６までの時間など）は、種々の時間に設定されるが、この実施例では、例えば１／３０秒としている。以降も同様にしてサブトラクション像の生成が行われる。
【００３７】
したがって、高周波数成分除去画像（マスク像）の生成元であるＸ線透過像（ライブ像）ではなく、その高周波数成分除去画像（マスク像）が生成された時刻以降のＸ線透過像である基本画像（ライブ像）から、高周波数成分除去画像（マスク像）をサブトラクションしているので、周波数特性変換処理の所要時間分の遅延を含まない基本画像（ライブ像）をサブトラクションの対象とすることができ、除去されるべき低周波数成分の時間追従性は悪くなるものの、関心物である高周波数成分は高周波数成分除去画像が生成された時刻以降のＸ線透過像である基本画像（ライブ像）に含まれているため最小限の遅延でサブトラクション像などを表示することができ、被検体Ｍの動きに追従できるサブトラクション像を提供できる。
【００３８】
具体的には、従来例装置では、図１０（ａ）に示すように、Ｘ線透過像（ライブ像）ごとに高周波数成分除去画像を生成しているので、例えばＸ線透過像（ライブ像）「Ａ」の周波数特性変換処理の開始時刻ｔ１からサブトラクション像生成までの時刻ｔ５までの遅延時間（ｔ５−ｔ１の時間）が必要であったので、被検体Ｍの動きに追従したサブトラクション像を得ることができなかった。しかしながら、この第１実施例装置によれば、図１０（ｂ）に示すように、高周波数成分除去画像（マスク像）「Ａ´」と、この高周波数成分除去画像（マスク像）「Ａ´」が生成された時刻ｔ２より以降に撮影されたＸ線透過像である基本画像（ライブ像）例えば「Ｂ」とでサブトラクションするので、例えばＸ線透過像（ライブ像）「Ｂ」の撮影完了時刻ｔ３からサブトラクション像生成までの時刻ｔ５までの遅延時間（ｔ５−ｔ３の時間）で済み、被検体Ｍの動きに追従したサブトラクション像を得ることができる。
【００３９】
次に、被検体ＭとＸ線透視装置２との相対的な位置関係を被検体Ｍの体軸方向に変位させて、例えば、図５に示すように、胸部から腹部にかけての領域ＳＲ内の複数の撮像部位のサブトラクション像を得る場合の動作を説明する。なお、この実施例では、Ｘ線透視装置２を固定し、これに対して被検体Ｍを載置した天板１２を水平移動するように構成しているが、撮像状況をわかり易くするために、図５では、天板１２上の被検体Ｍを固定し、これに対してＸ線透視装置２（Ｘ線管２１、撮像系２２）が変位しているように描いている。
【００４０】
この場合、制御部５は、最初の撮像部位（図では、撮像領域ＳＲの左端部側の撮像部位）を撮像位置に位置させ、撮像方向を調節する。そして、以下の撮像動作の前に被検体Ｍに造影剤を投与しておく。
【００４１】
被検体Ｍの各撮像部位（領域ＳＲ）に造影剤が拡散し、処理開始が指示されると、上記１箇所の撮像領域ＳＢのサブトラクション像を求めた手順と同様の手順で、最初の撮像部位のサブトラクション像を求め、天板１２を図５の左方向に定速で移動させながら、以降の各撮像部位が撮像位置に位置するごとに、その撮像部位のサブトラクション像を順次求めていく。
【００４２】
したがって、被検体ＭとＸ線透視装置との相対的な位置関係を被検体Ｍの体軸方向に変位させて、複数の連続する撮像部位に対するサブトラクション像を得る場合であっても、前述と同様の効果を有する。
【００４３】
なお、天板１２を固定し、Ｘ線透視装置２を天板１２上の被検体Ｍの体軸方向に移動させることで、被検体ＭとＸ線透視装置との相対的な位置関係を被検体Ｍの体軸方向に変位させるように構成してもよい。
【００４４】
また、被検体Ｍのある部位（例えば、胸部）を撮像位置に位置させた状態で、図６に示すように、Ｘ線管２１、撮像系２２をその部位の回り（体軸回り）に回転変位させながら、各撮像方向からのサブトラクション像を求めることもあるが、このような場合であっても、上記各動作と同様に、被検体Ｍの動きに追従したサブトラクション像を得ることができる。
【００４５】
また、実空間フィルタリング部３２ａは、実空間上で直線形状で表されるテンプレートフィルタでもって、Ｘ線透過像（ライブ像）の高周波数成分を実空間上で除去しており、次に説明するような効果がある。
【００４６】
例えば、従来のガウス関数形状のテンプレートフィルタでは、図１７に示すそのガウス曲線形状からもわかるように、Ｘ線透過像（ライブ像）を単純な加減算のみで処理することできないし、フィルタサイズ（Ｎ×Ｎ点）が例えば５１（ＴＡＰ数）×５１（ＴＡＰ数）点のように大きくなることでフィルタリングの演算量が膨大になり、従来例の場合における演算量は次に示す式（１）で表されて、リアルタイムに近い短時間での処理は不可能であった。なお、式（１）中のＰは、１回の乗算と１回の加算との２回であり、式（１）中のＧ² は、Ｘ線透過像（ライブ像）の全画素数である。例えば、参考までに式（１）に具体的な数値を代入してみる。

【００４７】
これに対して、上述した第１実施例では、実空間フィルタリング部３２ａは、実空間上で直線形状で表される、四角柱形状のテンプレートフィルタでもって、Ｘ線透過像（ライブ像）をフィルタリングするので、Ｘ線透過像（ライブ像）を単純な加減算のみで処理することできるし、さらに、フィルタリング処理における重複する演算（加算）を繰り返し実行することがないように演算量を軽減することができる。この第１実施例の場合におけるフィルタリングの演算量は次に示す式（２）で表される。なお、式（２）中のＱは、１回の加算と１回の減算との２回であり、式（２）中のＧ² は、前述の式（１）と同様に、Ｘ線透過像（ライブ像）の全画素数である。例えば、参考までに式（２）に具体的な数値を代入してみる。

【００４８】
式（１），（２）を比してわかるように、第１実施例の場合におけるフィルタリングの演算量は、フィルタ移動によって更新される部分についての単純な加減算を実行するだけで良く、式（２）中にＴＡＰ数が含まれていないことから、フィルタサイズに依らずに一定量でしかも式（１）に比べて大幅に低減され、リアルタイムに近い短時間での処理が可能となることがわかる。
【００４９】
＜第２実施例＞
続いて、この発明のＸ線ディジタルサブトラクション装置の第２実施例に係るＸ線ディジタルアンギオグラフィ装置について説明する。図１１は、この発明の第２実施例に係る実空間フィルタリング部４０のテンプレートフィルタ形状を示す模式図である。図１２は、この第２実施例に係る実空間フィルタリング部４０のブロック図である。この第２実施例のＸ線ディジタルアンギオグラフィ装置は、第１実施例の実空間フィルタリング部３２ａに替えて、図１２に示す実空間フィルタリング部４０を採用し、第２実施例の実空間フィルタリング部４０のテンプレートフィルタ形状を、四角柱形状から後述するような八角錐形状に替えた点以外については、前述の第１実施例と同様であるので、特にテンプレートフィルタ形状が実空間上で後述するような八角錐形状である実空間フィルタリング部４０の構成および機能について詳細に説明するものとする。
【００５０】
まず、実空間フィルタリング部４０のテンプレートフィルタ形状について説明する。実空間フィルタリング部４０のテンプレートフィルタ形状は、図１１に示すように、互いに直交するＸ，Ｙ，Ｚ軸で表される、実空間上の３次元直交座標系において、Ｘ，Ｙ軸方向をテンプレートフィルタの矩形形状である底面大きさとし、Ｚ軸方向をこのテンプレートのゲインとし、Ｚ軸方向に先細りとなる立体形状であって、かつ、テンプレートフィルタのＺ方向の少なくとも中央から先細りの先端部にかけて八角錐形状で表されるものとしている。なお、このような実空間フィルタリング部４０のテンプレートフィルタ形状を、説明の便宜上、略八角錐形状とも呼ぶこととする。
【００５１】
実空間フィルタリング部４０は、図１２に示すように、入力されるＸ線透過像（ライブ像）の行方向に後述する所定のフィルタをかける行方向三角形フィルタ４１と、この行方向三角形フィルタ４１でフィルタ処理されたデータを記憶する記憶手段４２と、この記憶手段４２に記憶された行方向処理後のＸ線透過像（ライブ像）の列方向に後述する所定のフィルタをかける列方向三角形フィルタ４３と、行方向三角形フィルタ４１と記憶手段４２と列方向三角形フィルタ４３とを制御する制御手段４４とを備えている。
【００５２】
行方向三角形フィルタ４１は、図１３に示すように、更新される部分のうちで新たに加えられる、後述する三角形状の右側の下がり勾配部分の画素の値（図１５（ｂ）参照）を生成する第１の１次元フィルタブロック５１と、更新される部分のうちで削除される、後述する三角形状の左側の上がり勾配部分の画素の値（図１５（ｂ）参照）を生成する第２の１次元フィルタブロック５２と、これらの第１の１次元フィルタブロック５１と第２の１次元フィルタブロック５２とからの画像データを減算する減算器５３と、この減算器５３からの出力を累算して出力する累算回路５４と、入力される基本画像を第１の１次元フィルタブロック５１よりも所定量遅延させて第２の１次元フィルタブロック５２に入力するための遅延回路５５とを備えている。
【００５３】
第１の１次元フィルタブロック５１は、一方の入力ポートに入力される画像データ（基本画像の行方向についての一連の画像データ）から、他方の入力ポートから入力される画像データ（基本画像の行方向についての一連の画像データ）を減算して出力する減算器６１と、入力される基本画像を減算器６１の一方の入力ポートよりも所定量遅延させて減算器６１の他方の入力ポートに入力するための遅延回路６２と、この減算器６１からの出力を累算して出力する累算回路６３とを備えている。なお、第２の１次元フィルタブロック５２も、第１の１次元フィルタブロック５１と同様に構成されており、減算器６１と累算回路６２と遅延回路６３とを備えている。
【００５４】
なお、列方向三角形フィルタ４３は、前述の行方向三角形フィルタ４１と同様に構成されており、行方向三角形フィルタ４１で処理されて記憶手段４２に記憶されたＸ線透過像（ライブ像）の行方向のデータを、列方向にフィルタリング処理するものである。
【００５５】
ここで、テンプレートフィルタ形状が実空間上で略八角錐形状である実空間フィルタリング部４０について、そのフィルタリング機能について説明する。この実空間フィルタリング部４０は、実空間の１次元形状が三角形をしている図１４（ａ）に示す１次元フィルタを、図１１に示すように２次元に拡張したものと等価である。図１４（ｃ）に示すように、行列（縦横）分離型のフィルタで構成することで、高速化に有利となる。なお、フィルタ係数の一例を図１４（ｂ），（ｃ）に示す。例えば、図１４（ａ）に示すように実空間形状が三角形である１次元フィルタを５点のものとすると、図１４（ｂ）に示すように左から順に「１，２，３，２，１」の係数となる。図１４（ｂ）に示す行（横）方向の１次元フィルタを列（縦）方向に並べると、図１４（ｃ）の右上側に示すような２次元フィルタとなる。また、図１４（ｂ）に示す行（横）方向の１次元フィルタを縦方向にしたもの、つまり上から順に「１，２，３，２，１」の係数としたものを、行（横）方向に並べると、図１４（ｃ）の左上側に示すような２次元フィルタとなる。そして、これらの２次元フィルタをその同一点同士の係数をかけることで、図１４（ｃ）の下側に示すような、行列（縦横）分離型のフィルタ、つまり、実空間上で略八角錐形状（略ピラミッド形状）であるテンプレートフィルタと等価な効果が得られる。
【００５６】
なお、図１４に示してきた係数（「１，２，３，４，６，９」）は、図１５に説明するように、重み度を示すものである。すなわち、図１５（ａ）に示すように、図１４（ｂ）に示す１次元フィルタがＸ線透過像（ライブ像）の「Ａ」〜「Ｅ」の画素に位置する場合に、それらの画素の値を何倍するかを示している。具体的に、画素「Ａ」は図１４（ｂ）に示す１次元フィルタの係数「１」にあるので「Ａ」を１倍したもの、つまり「Ａ」のままとなり、画素「Ｂ」は図１４（ｂ）に示す１次元フィルタの係数「２」にあるので「Ｂ」を２倍したもの、つまり「２Ｂ」となり、画素「Ｃ」は図１４（ｂ）に示す１次元フィルタの係数「３」にあるので「Ｃ」を３倍したもの、つまり「３Ｃ」となり、画素「Ｄ」は図１４（ｂ）に示す１次元フィルタの係数「２」にあるので「Ｄ」を２倍したもの、つまり「２Ｄ」となり、画素「Ｅ」は図１４（ｂ）に示す１次元フィルタの係数「１」にあるので「Ｅ」を１倍したもの、つまり「Ｅ」のままとなる。
【００５７】
図１４（ｃ）の下側に示すような、行列（縦横）分離型のフィルタ、つまり、実空間上で略八角錐形状（略ピラミッド形状）であるテンプレートフィルタは、列（縦）方向と行（横）方向とに分離すると、図１５に示すように、前述の第１実施例での移動平均のアルゴリズムの延長で計算でき、フィルタリング処理における重複する演算（加算）を繰り返し実行することがないように演算量を軽減することができる。
【００５８】
すなわち、第１の１次元フィルタブロック５１は、図１５（ｂ）に示すように、更新される部分のうちで新たに加えられる、三角形状の下がり勾配部分の画素の値を生成する。図１５（ｂ）では、更新される部分のうちで新たに加えられる、三角形状の下がり勾配部分の画素の値は、「Ｄ＋Ｅ＋Ｆ」である。
【００５９】
また、第２の１次元フィルタブロック５２は、図１５（ｂ）に示すように、更新される部分のうちで削除される、三角形状の上がり勾配部分の画素の値を生成する。図１５（ｂ）では、更新される部分のうちで削除される、三角形状の上がり勾配部分の画素の値は、「Ａ＋Ｂ＋Ｃ」である。
【００６０】
図１４（ｂ）に示す１次元フィルタがＸ線透過像（ライブ像）の「Ａ」〜「Ｅ」の画素に位置する場合には、図１５（ａ）に示すように、重み付けされたもの「Ａ＋２Ｂ＋３Ｃ＋２Ｄ＋Ｅ」が得られる。次に、図１４（ｂ）に示す１次元フィルタが画素１つ分移動してＸ線透過像（ライブ像）の「Ｂ」〜「Ｆ」の画素に位置する場合には、図１５（ａ）に示したものから「Ａ＋Ｂ＋Ｃ」を引くとともに「Ｄ＋Ｅ＋Ｆ」を加えることで、図１５（ｂ）に示すように、重み付けされたもの「Ｂ＋２Ｃ＋３Ｄ＋２Ｅ＋Ｆ」が得られる。したがって、図１４（ｂ）に示す１次元フィルタが画素１つ分移動してＸ線透過像（ライブ像）の「Ｂ」〜「Ｆ」の画素に位置する場合に、最初から正直に演算して、つまり、前回（図１５（ａ）に示すもの）と今回（図１５（ｂ）に示すもの）とで重複する演算を行なうことで、図１５（ｂ）に示す重み付けされたもの「Ｂ＋２Ｃ＋３Ｄ＋２Ｅ＋Ｆ」を得る必要はない。次に、図１４（ｂ）に示す１次元フィルタがさらに画素１つ分移動してＸ線透過像（ライブ像）の「Ｃ」〜「Ｇ」の画素に位置する場合には、同様に、図１５（ｂ）に示したもの「Ｂ＋２Ｃ＋３Ｄ＋２Ｅ＋Ｆ」から「Ｂ＋Ｃ＋Ｄ」を引くとともに「Ｅ＋Ｆ＋Ｇ」を加えて重み付けされたもの「Ｃ＋２Ｄ＋３Ｅ＋２Ｆ＋Ｇ」が得られる。
【００６１】
このように、フィルタ移動によって更新される部分の加減算は、３画素分の加算および３画素の減算のみであるので、この更新される部分の加減算を行なうだけで良い。
【００６２】
次に、行方向三角形フィルタ４１によって、図１５（ｂ）に示す「Ｂ＋２Ｃ＋３Ｄ＋２Ｅ＋Ｆ」を得る動作を、具体的に説明する。なおこの図１５では、説明の便宜上、フィルタサイズを５点（ＴＡＰ数＝５点）としているので、第１の１次元フィルタブロック５１の遅延回路６２と、第２の１次元フィルタブロック５２の遅延回路６２と遅延回路５５との遅延量は、次に示す式（３）により、それぞれ「３」に設定されている。

【００６３】
これらの遅延回路のパラメータ（遅延量）などがデジタルフィルタで言うところのＴＡＰ数に相当する。ＴＡＰ数を変えることで、しきい値周波数の値を変更できる。なお、図１５では、説明の便宜上、フィルタサイズを５点としていたが、この第２実施例では、図１１に示すようにＴＡＰ数を５１×５１点とした場合には、第１の１次元フィルタブロック５１の遅延回路６２と、第２の１次元フィルタブロック５２の遅延回路６２と遅延回路５５との遅延量は、「２６」に設定される。制御手段４４は、第１の１次元フィルタブロック５１の累算回路６３と、第２の１次元フィルタブロック５２の累算回路６３と、累算回路５４との累算値の初期化や、第１の１次元フィルタブロック５１の遅延回路６２と、第２の１次元フィルタブロック５２の遅延回路６２と遅延回路５５との遅延量の設定を行なう。
【００６４】
例えば、Ｘ線透過像（ライブ像）の「Ａ」〜「Ｆ」が行方向三角形フィルタ４１に入力された時点での、この行方向三角形フィルタ４１の各構成での処理状況について見てみる。第１の１次元フィルタブロック５１の減算器６１には「Ｆ」と「Ｃ」とが入力されて「Ｆ−Ｃ」が第１の１次元フィルタブロック５１の累算回路６３に入力される。第１の１次元フィルタブロック５１の累算回路６３は、（「Ｃ＋Ｄ＋Ｅ」＋「Ｆ−Ｃ」）により「Ｄ＋Ｅ＋Ｆ」が出力される。また、遅延回路５５の遅延量が「３」であるので、第２の１次元フィルタブロック５２の減算器６１には未だ「Ｃ」のみが入力されるだけで「Ｃ」が第２の１次元フィルタブロック５２の累算回路６３に入力される。第２の１次元フィルタブロック５２の累算回路６３は、（「Ａ＋Ｂ」＋「Ｃ」）により「Ａ＋Ｂ＋Ｃ」が出力される。減算器５３では、第１の１次元フィルタブロック５１の累算回路６３から出力された「Ｄ＋Ｅ＋Ｆ」と、第２の１次元フィルタブロック５２の累算回路６３から出力された「Ａ＋Ｂ＋Ｃ」との減算が行なわれ、（「Ｄ＋Ｅ＋Ｆ」−「Ａ＋Ｂ＋Ｃ」）が出力される。累算回路５４は、前回の値「Ａ＋２Ｂ＋３Ｃ＋２Ｄ＋Ｅ」と、減算器５３からの（「Ｄ＋Ｅ＋Ｆ」−「Ａ＋Ｂ＋Ｃ」）とを累算し、「Ｂ＋２Ｃ＋３Ｄ＋２Ｅ＋Ｆ」を記憶手段４２に出力する。
【００６５】
以上、上述した第２実施例では、実空間フィルタリング部４０は、実空間上で直線形状で表され、略八角錐形状のテンプレートフィルタでもって、Ｘ線透過像（ライブ像）の高周波数成分を実空間上で除去しており、更新される部分の加減算のみを行なえば良いので、前述の第１実施例の場合と同等に、フィルタリングの演算量は、式（２）で表され、この式（２）中にＴＡＰ数が含まれていないことから、フィルタサイズに依らずに一定量でしかも、従来例の式（１）に比べて大幅に低減され、リアルタイムに近い短時間での処理が可能となる効果を有している。さらに、この第２実施例では、次に説明する点で、前述の第１実施例よりもさらに優れている。
【００６６】
すなわち、第２実施例のテンプレートフィルタ形状（図１１参照）は、前述の第１実施例のテンプレートフィルタ形状（図７参照）に比べて、理想的なテンプレートフィルタ形状（図１７参照）に近似しているので、図１６に示すように、前述の第１実施例の場合よりも、周波数特性に優れている。具体的には、前述の第１実施例のテンプレートフィルタでは０．２５ｌｐ／ｍｍ（ラインペア／ミリ）付近でリンギング（ＭＡＸが０．１５程度の山）が生じているが、このようなリンギングは第２実施例のテンプレートフィルタではかなり小さくなっていて僅かに存在しているだけであることから、前述の第１実施例に比べて、所定のしきい値周波数以上の周波数成分を良好に除去できることがわかる。
【００６７】
また、しきい値周波数の値の変更が、遅延回路の遅延量の調整で実現でき、演算量は増減しないため、しきい値周波数の値に依らず、低減した一定時間で処理できる。したがって、所定の処理時間を確保したまま、しきい値周波数の値の変更の自由度が高いディジタルサブトラクション装置が提供できる。
【００６８】
なお、この第２実施例では、実空間フィルタリング部４０のテンプレートフィルタ形状を、図１１に示すように、テンプレートフィルタのＺ方向の少なくとも中央から先細りの先端部にかけて八角錐形状で表されるものとした略八角錐形状としているが、底面から先細りの先端部までの全部を八角錐形状としても良い。
【００６９】
＜第３実施例＞
図１８は、第３実施例装置の画像処理部の構成を示すブロック図である。この第３実施例は、前述の第１，第２実施例装置の供給制御回路３３に替えて、後述する供給制御回路３３ａを設けたことを特徴とする。この供給制御回路３３ａは、高周波数成分除去画像（マスク像）と、別途撮影されたＸ線透過像としての、この高周波数成分除去画像（マスク像）が生成されるまでに撮影されたＸ線透過像である基本画像（ライブ像）とを演算器３４に同期して供給するものである。この第３実施例装置では、供給制御回路３３ａのみが前述の第１，第２実施例装置とは異なるので、この第３実施例ではこの供給制御回路３３ａを詳細に説明するものとする。
【００７０】
図１８に示すように、供給制御回路３３ａは、例えば、前述の第１，第２実施例と同様の周波数特性変換回路３２と、この周波数特性変換回路３２で生成された高周波数成分除去画像（マスク像）を所定時間記憶する記憶回路７１と、基本画像（ライブ像）を所定時間遅延させて出力する遅延回路７２とを備えている。
【００７１】
ここで、この第３実施例装置のサブトラクション処理動作を、図１９を用いて説明する。
【００７２】
図１９に示すように、Ｘ線透過像（ライブ像）「Ａ」が撮像して得られると（時刻ｔ１１）、このＸ線透過像（ライブ像）「Ａ」の周波数特性を変換して高周波数成分除去画像（マスク像）「Ａ´」を生成する周波数特性変換処理が行われる。この周波数特性変換処理の終了時刻がｔ１４であったとすると、この第３実施例では、時刻ｔ１１〜ｔ１４の間に、Ｘ線透過像（ライブ像）「Ａ」とは別のＸ線透過像（ライブ像）である基本画像（ライブ像）「Ａ＊」が撮影されている（時刻ｔ１３）。供給制御回路３３ａは、Ｘ線透過像（ライブ像）「Ａ」から高周波数成分除去画像（マスク像）「Ａ´」をサブトラクションしてサブトラクション像「Ａ−Ａ´」が生成され、基本画像（ライブ像）「Ａ＊」から前記と同じ高周波数成分除去画像（マスク像）「Ａ´」をサブトラクションしてサブトラクション像「Ａ＊−Ａ´」が生成されるように、記憶回路７１と遅延回路７２とを制御することで、演算器３４へのライブ像とマスク像との同期供給を制御している。なお、次のサブトラクション像は、Ｘ線透過像（ライブ像）「Ｂ」と、このＸ線透過像（ライブ像）「Ｂ」から求めた高周波数成分除去画像「Ｂ´」とをサブトラクションしてサブトラクション像「Ｂ−Ｂ´」を生成し、Ｘ線透過像（ライブ像）である基本画像「Ｂ＊」と、前記と同じ高周波数成分除去画像「Ｂ´」とをサブトラクションしてサブトラクション像「Ｂ＊−Ｂ´」を生成している。あるＸ線透過像（ライブ像）から次のＸ線透過像（ライブ像）を撮影するまで時間（例えば、時刻ｔ１０〜ｔ１２までの時間など）は、種々の時間に設定されるが、この実施例では、例えば１／６０秒としている。以降も同様にしてサブトラクション像の生成が行われる。
【００７３】
したがって、複数枚のＸ線透過像（ライブ像）に対して個別に同じ高周波数成分除去画像（マスク像）でサブトラクションすることができ、撮影して得られたＸ線透過像（ライブ像）ごとに高周波数成分除去画像（マスク像）を生成する必要が無いので、その時間分周波数変換処理時間を低減することができ、撮影レートを向上させる、すなわち、表示間隔を短くすることができ、被検体の動きに追従できるサブトラクション像を提供できる。
【００７４】
この発明は、上記の各実施例に限られるものではなく、下記のように変形実施することができる。
【００７５】
（１）上述の第１実施例の供給制御回路３３と第３実施例の供給制御回路３３ａとを組み合わせて構成しても良い。すなわち、この組み合わせた供給制御回路は、図２０に示すように、高周波数成分除去画像（マスク像）と、別途撮影されたＸ線透過像としての、この高周波数成分除去画像（マスク像）が生成されるまでに撮影されたＸ線透過像である基本画像（ライブ像）とを演算器３４に同期して供給するとともに、高周波数成分除去画像（マスク像）と、別途撮影されたＸ線透過像としての、この高周波数成分除去画像（マスク像）が生成された時刻より以降に撮影されたＸ線透過像である基本画像（ライブ像）とも演算器３４に同期して供給するものである。この場合は、図２０に示すように、遅延を低減し、かつ、撮影レートを向上させたディジタルサブトラクション装置を構成することができる。
【００７６】
（２）上述の各実施例のディジタルサブトラクション装置における実空間フィルタリング３２ａ，４０を汎用チップでもってソフトウエア的に実現することもできる。
【００７７】
（３）上述の各実施例のディジタルサブトラクション装置では、周波数特性変換回路３２に実空間フィルタリング部３２ａ，４０を用いているが、Ｘ線透過像（ライブ像）を構成する周波数成分を空間／周波数変換処理で取り出し、所定のしきい値周波数以上の周波数成分を除去し、これに周波数／空間変換処理を施して、前記Ｘ線透過像（ライブ像）から前記しきい値周波数以上の周波数成分を除去した高周波数成分除去画像（マスク像）を得る周波数特性変換回路３２を採用しても良い。前記の空間／周波数変換処理としては、ＦＦＴ（高速フーリエ変換）、カルーネン・レーベ変換、ＤＣＴ（離散コサイン変換）、アダマール変換などの各変換方式がある。また、前記の周波数／空間変換処理としては、前記空間／周波数変換処理の逆変換（逆ＦＦＴ、逆カルーネン・レーベ変換、逆ＤＣＴ、逆アダマール変換など）がある。
【００７８】
（４）上述の各実施例のディジタルサブトラクション装置は、上述の各実施例のように、被検体Ｍを人体などとして医療用に用いることもできるし、被検体ＭをＢＧＡ（Ball Grid Array）基板やプリント配線基板など各種の電子部品などとして非破壊検査用に用いることもできる。
【００７９】
【発明の効果】
以上の説明から明らかなように、請求項１に記載のディジタルサブトラクション装置によれば、高周波数成分除去画像の生成元であるＸ線透過像ではなく、別途撮影されたＸ線透過像である基本画像から、高周波数成分除去画像をサブトラクションしているので、周波数特性変換処理の所要時間分の遅延を含まない基本画像をサブトラクションの対象とすることができ、除去されるべき低周波数成分の時間追従性は悪くなるものの、関心物である高周波数成分は別途に撮影された基本画像に含まれているため最小限の遅延でサブトラクション像などを表示することができ、被検体の動きに追従できるサブトラクション像を提供できる。
【００８０】
また、請求項２に記載のディジタルサブトラクション装置によれば、高周波数成分除去画像の生成元であるＸ線透過像ではなく、その高周波数成分除去画像が生成された時刻以降のＸ線透過像である基本画像から、高周波数成分除去画像をサブトラクションしているので、周波数特性変換処理の所要時間分の遅延を含まない基本画像をサブトラクションの対象とすることができ、除去されるべき低周波数成分の時間追従性は悪くなるものの、関心物である高周波数成分は高周波数成分除去画像が生成された時刻以降の基本画像に含まれているため最小限の遅延でサブトラクション像などを表示することができ、被検体の動きに追従できるサブトラクション像を提供できる。
【００８１】
また、請求項３に記載のディジタルサブトラクション装置によれば、複数枚の基本画像に対して個別に同じ高周波数成分除去画像でサブトラクションすることができ、撮影して得られたＸ線透過像ごとに高周波数成分除去画像を生成する必要が無いので、その時間分周波数変換処理時間を低減することができ、撮影レートを向上させる、すなわち、表示間隔を短くすることができ、被検体の動きに追従できるサブトラクション像を提供できる。
【図面の簡単な説明】
【図１】この発明の第１実施例に係るディジタルアンギオグラフィ装置の全体構成を示す正面図である。
【図２】Ｘ線透視装置を側面から見た図である。
【図３】第１実施例装置に備えられた画像処理部の構成を示すブロック図である。
【図４】ある１箇所の撮像部位のサブトラクション像を得る場合の動作を説明するための図である。
【図５】被検体の体軸方向の複数の撮像部位のサブトラクション像を得る場合の動作を説明するための図である。
【図６】ある部位に対する複数の撮像方向からのサブトラクション像を得る場合の動作を説明するための図である。
【図７】第１実施例のテンプレートフィルタ形状を示す模式図である。
【図８】第１実施例の移動平均フィルタ法を説明するための模式図である。
【図９】第１実施例の移動平均フィルタの計算アルゴリズムを説明するための模式図である。
【図１０】（ａ）は従来例によるサブトラクション処理動作シーケンスを示す模式図であり、（ｂ）は第１実施例によるサブトラクション処理動作シーケンスを示す模式図である。
【図１１】第２実施例のテンプレートフィルタ形状を示す模式図である。
【図１２】この発明の第２実施例に係る実空間フィルタリング部のブロック図である。
【図１３】この発明の第２実施例に係る行方向三角形フィルタの構成を示すブロック図である。
【図１４】（ａ）〜（ｃ）は第２実施例のピラミッド形フィルタを説明するための模式図である。
【図１５】（ａ）、（ｂ）は第２実施例のピラミッド形フィルタの計算アルゴリズムを説明するための模式図である。
【図１６】各種フィルタの周波数特性を示す特性図である。
【図１７】従来のガウス関数形状で表される理想的なテンプレートフィルタ形状を示す模式図である。
【図１８】第３実施例装置の画像処理部の構成を示すブロック図である。
【図１９】第３実施例装置によるサブトラクション処理動作シーケンスを示す模式図である。
【図２０】第１，第３実施例の供給制御回路を組み合わせた場合のサブトラクション処理動作シーケンスを示す模式図である。
【符号の説明】
２ … Ｘ線透視装置
３ … 画像処理部
３１ … Ａ／Ｄ変換器
３２ … 周波数特性変換回路
３３ … 供給制御回路
３３ａ… 供給制御回路
３４ … 演算器
Ｍ … 被検体

Claims (3)

1. 被検体の所定の撮像部位のサブトラクション像を得るためのディジタルサブトラクション装置であって、（ａ）前記撮像部位にＸ線を照射し、その部位の複数枚のＸ線透過像を撮像するＸ線透視手段と、（ｂ）前記Ｘ線透過像をディジタルデータに変換するデータ変換手段と、（ｃ）前記ディジタルデータに変換されたＸ線透過像ごとに所定のしきい値周波数以上の周波数成分を除去した画像（以下、高周波数成分除去画像）を得る周波数特性変換手段と、（ｄ）Ｘ線透過像である各基本画像と各高周波数成分除去画像とのサブトラクションを行い、その撮像部位のサブトラクション像を求める演算手段と、（ｅ）前記高周波数成分除去画像と、この高周波数成分除去画像を生成するためのＸ線透視像の所定回後に別途撮影されたＸ線透過像である基本画像とを前記演算手段に同期して供給する供給制御手段とを備えたことを特徴とするディジタルサブトラクション装置。
2. 請求項１に記載のディジタルサブトラクション装置において、前記供給制御手段は、前記高周波数成分除去画像と、別途撮影されたＸ線透過像としての、この高周波数成分除去画像が生成された時刻より以降に撮影されたＸ線透過像である基本画像とを前記演算手段に同期して供給することを特徴とするディジタルサブトラクション装置。
3. 請求項１または請求項２に記載のディジタルサブトラクション装置において、前記供給制御手段は、前記高周波数成分除去画像と、別途撮影されたＸ線透過像としての、この高周波数成分除去画像が生成されるまでに撮影されたＸ線透過像である基本画像とを前記演算手段に同期して供給することを特徴とするディジタルサブトラクション装置。

Images