(57)【特許請求の範囲】
【請求項１】
被検体に対してＸ線を照射するＸ線発生部と前記被検体を透過したＸ線を検出するＸ線検出部とを含む撮影部と、
前記被検体の呼吸周期に合わせて前記被検体を撮影し、造影剤なしで撮影したマスク画像群と、造影剤を投入して撮影したコントラスト画像群を収集する画像収集部と、
前記マスク画像群及びコントラスト画像群のそれぞれ前後するフレームの画像の差分画像を生成する差分画像生成部と、
前記被検体の呼吸によって変動する前記差分画像の移動方向を検出し、前記差分画像が第１の方向に移動し反転するタイミングを起点にして呼吸周期を決定し、前記マスク画像群及びコントラスト画像群の各フレームの呼吸１周期中の位置を示す呼吸位相値を算出し、前記呼吸位相値が近似する前記マスク画像と前記コントラスト画像を選択して減算処理しサブトラクション画像を生成する演算処理部と、
前記演算処理部で生成したサブトラクション画像を表示する表示部と、を具備したことを特徴とするＸ線診断装置。
【請求項２】
前記演算処理部は、前記差分画像生成部で生成した差分画像に線状のＲＯＩを設定し、前記ＲＯＩの線上にある画素データからプロファイルカーブを取得し、前記プロファイルカーブの変化特性をもとに前記差分画像の移動方向を検出することを特徴とする請求項１記載のＸ線診断装置。
【請求項３】
前記演算処理部は、前記線状のＲＯＩを撮影部位に応じて予め設定することを特徴とする請求項２記載のＸ線診断装置。
【請求項４】
前記演算処理部は、前記差分画像が前記第１の方向に移動して反転するタイミングを起点に、前記マスク画像群及びコントラスト画像群の各フレームにフレーム番号を付し、前記フレーム番号をインデックスとしてサブトラクションに使用するマスク画像とコントラスト画像を抽出する特徴とする請求項１記載のＸ線診断装置。
【請求項５】
被検体に対してＸ線を照射するＸ線発生部と前記被検体を透過したＸ線を検出するＸ線検出部とを含む撮影部を備え、
前記被検体の呼吸周期に合わせて前記被検体を撮影し、造影剤なしで撮影したマスク画像群と、造影剤を投入して撮影したコントラスト画像群を収集し、
前記マスク画像群及びコントラスト画像群のそれぞれ前後するフレームの画像の差分画像を生成し、
前記被検体の呼吸によって変動する前記差分画像の移動方向を検出して、前記差分画像が第１の方向に移動し反転するタイミングを起点にして呼吸周期を決定し、前記マスク画像群及びコントラスト画像群の各フレームの呼吸１周期中の位置を示す呼吸位相値を算出し、
前記呼吸位相値が近似する前記マスク画像と前記コントラスト画像とを減算処理してサブトラクション画像を生成し、
前記サブトラクション画像を表示部に表示することを特徴とするＸ線診断方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、自然呼吸下において撮影したマスク画像データと、造影剤を注入した状態で撮影したコントラスト画像データとを減算処理して、サブトラクション画像を得るＸ線診断装置及びＸ線診断方法に関する。
【背景技術】
【０００２】
従来、医用画像診断装置としてＸ線診断装置がある。例えば消化器や循環器の検査では、血管の微細な状態を観察するために、血管内に造影剤を注入する前に撮影したマスク画像データと、造影剤を注入した状態で撮影したコントラスト画像データとを減算（サブトラクト）処理して、サブトラクション画像を得る撮影方法がある。この撮影方法は、一般にＤＳＡ（Digital Subtraction Angiography）撮影と呼ばれ、サブトラクション画像を得ることで造影された血管のみを表示することができる。
【０００３】
ところで、ＤＳＡ撮影では、被検者が動くことによって、サブトラクションで背景をうまく消せずにアーチファクト（偽像）が発生し、造影された血管像を明瞭に観察できないことがある。アーチファクトを少なくするため、血管造影の際には被検者に対して体を動かさないように指示をして撮影を行っている。また腹部の血管造影では、呼吸による強いアーチファクトが発生し、診断能を著しく低下させることから、造影撮影の間、呼吸をしないように指示をすることもある。 
しかしながら腹部の血管造影は、部位によっては数十秒要する場合があり、また高齢者では呼吸を長く止められない等の問題があるため、通常通りに呼吸をした状態で撮影を行う方法もある（これは、自然呼吸下ＤＳＡ撮影と呼ばれる）。この場合、造影前に呼吸１周期分のマスク画像を収集し、それに引き続き必要な造影画像（コントラスト画像群）を収集し、ポストプロセスにおいてコントラスト画像と背景が一致する画像をマスク画像群の中からマスクとして選択（リマスク処理）して、サブトラクションしている。これにより、アーチファクトを軽減させることができるが、マスク画像群の中からコントラスト画像と背景が一致するマスク画像を選択するのは人手によって選んでいるため、診断・治療に適する血管像の作成に多くの手間と時間がかかる。
【０００４】
特許文献１には、ポストプロセスでの手間を軽減するため、体にセンサ等を装着して呼吸の状態をモニタし、センサからの情報に基づいて特定した呼吸位相に応じてＸ線照射のタイミングを制御し、造影画像とマスク画像とを同じ呼吸位相で収集するようにしたＸ線診断装置が開示されている。しかしながら、特許文献１の例では、自然呼吸下のＤＳＡ撮影を行う場合に呼吸センサのような特別な器具を必要とするため、更なる改善が求められている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平７−８４７９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
従来、サブトラクション画像のアーチファクトを軽減させるために、マスク画像群の中からコントラスト画像と背景が一致する画像を選択しているが、この選択作業は人手によって行っているため、診断、治療に適する血管像の作成に多くの時間と手間がかかるという問題があった。また特許文献１のように自然呼吸下でＤＳＡ撮影を行う例もあるが、呼吸センサのような特別な器具を必要とするため、更なる改善が求められている。
【０００７】
本発明は上記事情に鑑みて成されたもので、自然呼吸下におけるＤＳＡ撮影において、アーチファクトの少ないサブトラクション画像を表示することができるＸ線診断装置及びＸ線診断方法を提供することを目的とする。
【課題を解決するための手段】
【０００８】
請求項1記載の本発明のＸ線診断装置は、被検体に対してＸ線を照射するＸ線発生部と前記被検体を透過したＸ線を検出するＸ線検出部とを含む撮影部と、前記被検体の呼吸周期に合わせて前記被検体を撮影し、造影剤なしで撮影したマスク画像群と、造影剤を投入して撮影したコントラスト画像群を収集する画像収集部と、前記マスク画像群及びコントラスト画像群のそれぞれ前後するフレームの画像の差分画像を生成する差分画像生成部と、前記被検体の呼吸によって変動する前記差分画像の移動方向を検出し、前記差分画像が第１の方向に移動し反転するタイミングを起点にして呼吸周期を決定し、前記マスク画像群及びコントラスト画像群の各フレームの呼吸１周期中の位置を示す呼吸位相値を算出し、前記呼吸位相値が近似する前記マスク画像と前記コントラスト画像を選択して減算処理しサブトラクション画像を生成する演算処理部と、前記演算処理部で生成したサブトラクション画像を表示する表示部と、を具備したことを特徴とする。
【０００９】
また請求項５記載の本発明のＸ線診断方法は、被検体に対してＸ線を照射するＸ線発生部と前記被検体を透過したＸ線を検出するＸ線検出部とを含む撮影部を備え、前記被検体の呼吸周期に合わせて前記被検体を撮影し、造影剤なしで撮影したマスク画像群と、造影剤を投入して撮影したコントラスト画像群を収集し、
前記マスク画像群及びコントラスト画像群のそれぞれ前後するフレームの画像の差分画像を生成し、前記被検体の呼吸によって変動する前記差分画像の移動方向を検出して、前記差分画像が第１の方向に移動し反転するタイミングを起点にして呼吸周期を決定し、前記マスク画像群及びコントラスト画像群の各フレームの呼吸１周期中の位置を示す呼吸位相値を算出し、前記呼吸位相値が近似する前記マスク画像と前記コントラスト画像とを減算処理してサブトラクション画像を生成し、前記サブトラクション画像を表示部に表示することを特徴とする。
【発明の効果】
【００１０】
本発明のＸ線診断装置では、自然呼吸下におけるＤＳＡ撮影で、アーチファクトの少ないサブトラクション画像を自動的に表示することができる。また煩雑なリマスク処理を繰り返す必要がないため、診断効率を向上することができる。
【図面の簡単な説明】
【００１１】
【図１】本発明の一実施形態に係るＸ線診断装置の構成を示すブロック図。
【図２】一実施形態に係るＸ線診断装置の外観を示す斜視図。
【図３】一実施形態における差分画像の取得方法を示す動作説明図。
【図４】差分画像をもとに抽出したプロファイルカーブの一例を示す説明図。
【図５】呼吸位相と呼吸位相値（％）を示す説明図。
【図６】サブトラクション画像の生成処理を示す説明図。
【図７】呼吸位相の他の算出例を示す説明図。
【図８】呼吸位相画像演算部の動作を説明するフローチャート。
【発明を実施するための形態】
【００１２】
以下、この発明の一実施の形態について図面を参照して詳細に説明する。
【実施例１】
【００１３】
図１は、本発明のＸ線診断装置（例：アンギオ装置）の全体構成を示すブロック図である。図１において、Ｘ線診断装置１００は、被検体Ｐに対してＸ線を発生するためのＸ線発生部１０と、被検体Ｐを透過したＸ線を２次元的に検出するとともに検出結果に基づいてＸ線投影データを生成するＸ線検出部２０を備えている。
【００１４】
Ｘ線発生部１０は、Ｘ線管１１とＸ線絞り器１２を含むＸ線照射部１３と、高電圧発生器１４を備えている。Ｘ線管１１は、Ｘ線を発生する真空管であり、陰極（フィラメント）から放出された電子を高電圧により加速してタングステン陽極に衝突させＸ線を発生する。Ｘ線絞り器１２は、撮影対象部位のみにＸ線を照射し、被検体Ｐに対する被曝線量の低減と、画質向上を目的として設けられたものである。これにより、後述するＸ線検出器２１は、Ｘ線絞り器１２によって形成された開口を通過し、被検体Ｐの関心領域を透過したＸ線によるＸ線像を検出することができる。
【００１５】
Ｘ線検出部２０は、Ｘ線検出器（ＦＰＤ）２１と、データ変換部２２を備えている。Ｘ線検出器２１は、被検体Ｐの関心領域を透過したＸ線を電荷に変換して蓄積するものであり、Ｘ線を検出する微小な検出素子を列方向及びライン方向に２次元的に配列して構成される。データ変換部２２はＡ／Ｄ変換器と、パラレル・シリアル変換器を含み、Ｘ線検出器２１から読み出された電荷をもとにＸ線投影データを生成する。
【００１６】
Ｘ線検出部２０に入射したＸ線は、データ変換部２２によってＸ線のパルス照射に同期してデジタルデータ化され、１つのパルスに対して１フレームのデジタルデータ（Ｘ線投影データ）が生成される。デジタルデータは後述する画像処理部３４で処理される。
【００１７】
Ｘ線発生部１０とＸ線検出部２０は、Ｃアーム４１（支持装置）に支持されており、Ｃアーム４１は、例えば被検体Ｐの体軸方向に移動可能である。さらにＸ線診断装置１００は、コンピュータシステム３０を有している。コンピュータシステム３０は、システム制御部３１、操作部３２、ユーザインターフェース３３、画像処理部３４、表示部３９を備えている。画像処理部３４は、記憶部３５、呼吸位相画像演算部３６、画像データ処理部３７、Ｄ／Ａ変換器３８を含む。
【００１８】
システム制御部３１は、ＣＰＵ３１１等を有し、Ｘ線診断装置１００の各ユニットを統括的に制御し、画像データの処理や表示を行なう。操作部３２は医師等のユーザが操作するもので、Ｘ線照射スイッチ３２１を含む。操作部３２からの操作信号は、ユーザインターフェース３３を介してシステム制御部３１に供給されＣＰＵ３１１によって処理され、高電圧制御部４０、Ｘ線絞り器１２、被検体Ｐを載置する天板４５の位置の制御等を行う。またシステム制御部３１は、Ｘ線検出器２１およびデータ変換部２２の制御も行う。
【００１９】
高電圧発生器１４は、Ｘ線管１１の陰極から発生した熱電子を加速するために、陽極と陰極の間に印加する高電圧を発生させるものであり、高電圧制御部４０は、操作部３２から入力されたＸ線照射条件に基いて、高電圧発生器１４の管電流／管電圧、照射時間、照射繰り返し周期等を制御する。
【００２０】
図２は、Ｘ線診断装置１００（アンギオ装置）の外観を示す斜視図である。図２において、Ｘ線発生部１０とＸ線検出部２０がＣアーム４１に対向支持されており、Ｃアーム４１は、アーム支持部４２を介してスタンド４３に支持されている。尚、Ｘ線発生部１０とＸ線検出部２０を含む部分は撮影部を構成する。
【００２１】
また、Ｃアーム４１に対して寝台４４が配置されており、寝台４４の天板４５には、被検体（図示せず）が載置され、天板４５の位置及び高さは制御可能である。これによりＣアーム４１は、被検体の体軸方向に沿って相対的に移動可能であって、かつ被検体の周囲に沿って矢印ａ方向に回転可能である。また撮影部（Ｘ線発生部１０とＸ線検出部２０）を矢印ｂ方向にスライドすることができる。
【００２２】
撮影部で生成したＸ線投影データは、画像処理部３４（詳細は後述）によって処理され、画像データを表示部３９に表示される。表示部３９は、例えば天井部に取り付けられており、複数のモニタで構成されている。
【００２３】
また、寝台４４には、操作ハンドル４６が取り付けられており、操作ハンドル４６を操作することで、天板４５の高さ位置の制御、Ｃアーム４１の移動・回転制御、Ｘ線の照射範囲の調整、表示部３９の表示制御等の各種の制御を行うことができる。尚、Ｃアーム４１は、床面に置かれたスタンド４３に支持された例を示しているが、Ｃアーム４１を天井部に移動可能に設置したタイプのＸ線診断装置１００であっても良い。
【００２４】
次に本発明の特徴部である画像処理部３４について説明する。本発明では、ＤＳＡ（ Digital Subtraction Angiography ）撮影の手法に特徴がある。ＤＳＡ撮影とは、デジタル画像処理によって、血管造影画像から造影前の背景を減算（サブトラクト）して、造影された血管のみを表示する撮影方法を意味する。
【００２５】
先ずＸ線撮影の基本動作について説明する。操作部３２のＸ線照射スイッチ３２１により透視または撮影の開始操作を行うと、Ｘ線発生部１０からＸ線がパルス状に断続的に照射され始める。Ｘ線管１１で発生したＸ線は、Ｘ線絞り器１２によって照射野を制限することができ、照射されたＸ線は被検体Ｐを透過し、Ｘ線検出部２０に入射する。Ｘ線検出部２０に入射したＸ線は、Ｘ線のパルス照射に同期してデジタルデータ化され、１つのパルスに対して１フレームのデジタルデータが生成される。ここで、ＣＰＵ１３１は、Ｘ線の発生の制御と同期して、Ｘ線検出器２１及びデータ変換部２２の制御も併せて行う。
【００２６】
次に画像処理部３４の動作を図３〜図６を参照して説明する。尚、図３〜図６においては、各種の符号（Ａ〜Ｇ）を付しているが、同一項目については同一符号を付して説明する。
【００２７】
ＤＳＡ撮影では、血管内に造影剤を注入する前に撮影したマスク画像データと、造影剤を注入した状態で撮影したコントラスト画像データとを減算（サブトラクト）処理して、サブトラクション画像を得る。１呼吸周期を含むマスク画像を収集できるシーケンスで撮影を行うと、マスク画像群に引き続いてコントラスト画像群（造影画像群）が生成され、記憶部３５に記録される。マスク画像群が何フレームあり、コントラスト画像群が何フレームあるかは、付帯情報記憶メモリ（図示せず）に記憶される。Ｘ線検出部２０及び記憶部３５は、マスク画像群とコントラスト画像群を収集する画像収集部を構成する。
【００２８】
呼吸位相画像演算部３６は、差分画像生成部３６１とサブトラクション演算処理部３６２を含む。差分画像生成部３６１は、記憶部３５のデジタルデータから、マスク画像群及びコントラスト画像群の各フレームの呼吸位相を算出する。サブトラクション演算処理部３６２は、コントラスト画像群の或るコントラスト画像が有している呼吸位相に最も近似する呼吸位相を持つマスク画像をマスク画像群から選択して減算処理し、サブトラクション画像を生成する。
【００２９】
サブトラクションされたデータは、空間フィルタなどの画像処理が行われた後、画像データ処理部３７で表示階調処理が施され、Ｄ／Ａ変換器３８でビデオ信号に変換され、表示部３９にリアルタイムに表示される。
【００３０】
マスク画像及びコントラスト画像の呼吸位相を特定する手段としては、マスク画像群とコントラスト画像群のそれぞれに対して、各画像の１フレーム前の画像と１フレーム後の画像の差分画像を取得し、差分画像のプロファイルカーブから１フレームの間での臓器や骨などの移動方向を検出し、移動方向の変化によって各画像フレームの呼吸位相を特定し、コントラスト画像に対応するマスク画像を自動的に決定してサブトラクションを行う。
【００３１】
図３は、差分画像の取得方法を示す動作説明図である。被検者が呼吸をしながら、例えば１秒間に数フレーム程度の画像収集レートでの撮影を行うと、横隔膜が上下に大きく動き、それに伴って臓器の位置および形が、主に上下に移動・変形する。図３（Ａ）は、被検者の呼吸状態を示し、横軸のタイミングｔ０〜ｔ１が呼吸の１周期を示し、タイミングｔ１〜ｔ２は次の１周期を示す。
【００３２】
ＤＳＡ撮影においては、図３（Ｂ）で示すように、造影剤を注入する前の背景だけの画像（マスク画像）を呼吸の１周期分撮影してマスク画像群を収集し、そのあとで造影剤を注入し、引き続き呼吸をしながらコントラスト画像（造影画像）を撮影し、コントラスト画像群を収集する。
【００３３】
図３（Ｃ）はマスク画像群における各フレームの画像を示しており、先頭フレームを除くフレームに対して、１つ前のフレームと１フレーム後の画像をそれぞれ差分画像生成部３６１で減算処理して差分画像を生成する。図３（Ｄ）は差分画像を示している。差分画像は、２つのフレーム間の違いが画像化されたものであり、動きのあった部分として主に臓器の辺縁部が抽出される。
【００３４】
図４は、差分画像をもとに抽出したプロファイルカーブを示している。プロファイルカーブは、画像上に線を引いた時、その線上にある画素データに対する画素値をプロットしたものである。即ち、図４（Ｄ）は、図３（Ｄ）の差分画像のうち、隣接する差分画像Ｄ１，Ｄ２を拡大して示しており、差分画像Ｄ１，Ｄ２上に上下方向の線状ＲＯＩ（Ｌ１）を設定したとき、その線上にある画素データに対する画素値をプロットすると、図４（Ｅ）のＥ１，Ｅ２で示すプロファイルカーブが得られる。
【００３５】
プロファイルカーブＥ１，Ｅ２は臓器の辺縁部に相当する部分が移動方向に応じて正または負の値を持つカーブとなる。プロファイルカーブＥ１，Ｅ２は、矢印で示すように横隔膜の上下方向の動きに合わせて移動する。またコントラスト画像群についてもマスク画像群と同様に、差分画像を取得してプロファイルカーブを取得する（図示は省略）。
【００３６】
マスク画像群及びコントラスト画像群の各フレームに対する差分画像ＤからプロファイルカーブＥを取得すると、プロファイルカーブの変化から呼吸による臓器の移動を判別することができ、呼吸１周期内での呼吸位相を特定することができる。また呼吸位相は、最も臓器が上方向に存在するタイミングを基点（０％）にして呼吸１周期の中のどの位置にあるかを呼吸位相値（％）で示すことができる。
【００３７】
図５は、呼吸位相と呼吸位相値（％）を示す説明図である。図５（Ａ〜Ｄ）は、図３（Ａ〜Ｄ）と同様の図であるが、図５（Ｆ）では呼吸位相を矢印で示し、呼吸によって臓器が上又は下方向に動いていること示している。矢印の無い期間は移動方向が反転する期間であり、臓器が最も上または最も下にある状態である。
【００３８】
図５（Ｇ）は呼吸位相値（％）を示している。臓器が最も上方向に存在するタイミングを０％（１００％）とすると、呼吸位相Ｆが上向きの矢印で示される期間は１００に向かって徐々に数値が増加する。移動方向が反転し呼吸位相Ｆが下向きの矢印で示される期間は、０％を起点にしてゼロから徐々に数値が増加する。また臓器が最も下方向に存在するタイミング（本例では４８％）を過ぎると再び１００に向かって増加する。したがって、収集した画像Ｃのフレームに呼吸位相値（％）を付すことで、臓器がどの位置にあるときの画像であるが分かる。
【００３９】
図６は、サブトラクション画像の生成処理を示す説明図である。図６（Ｂ）は収集したマスク画像とコントラスト画像群を示し、Ｇは呼吸位相値（％）をそれぞれ示している。図６で示すように、サブトラクション画像を得る際は、マスク画像群とコントラスト画像群の中から同じ呼吸位相値を有するフレームの画像を選択し、サブトラクション演算処理部３６２で減算処理することで、サブトラクション画像を得ることができる。 呼吸位相値が最も近似するマスク画像およびコントラスト画像を選択して減算することでアーチファクトの少ないサブトラクション画像を自動的に表示することができる。
【００４０】
従来では、サブトラクション表示で、体動などにより背景がうまく消えない場合に、別の画像フレームをマスク画像として再選択（リマスク）する必要があったが、本発明の実施形態によれば、煩雑なリマスク処理を繰り返す必要がないため、診断効率を向上することができる。また呼吸位相を検出するため特別なセンサ等の器具を使用することなく、自然呼吸下におけるＤＳＡ撮影で、アーチファクトの少ないサブトラクション画像を自動的に表示することができる。
【００４１】
次に本発明の変形例について説明する。以上の実施形態では、隣り合うフレームの差分画像のプロファイルカーブをもとに、フレームの進みによって変化する方向を判断し、呼吸位相値（％）を算出している。 変形例は、差分画像のプロファイルカーブの正負の符号の分布によって、１枚の差分画像だけで臓器の移動方向を検出するようにしたものである。図７を参照して具体的に説明する。
【００４２】
図７（Ｃ）はマスク画像群における各フレームの画像を概略的に示しており、１つ前のフレームと１フレーム後の画像をそれぞれ減算処理して差分画像（Ｄ）を生成する。図７では、例えば臓器が上方向に移動しているときの差分画像Ｄ１と、下方向に移動しているときの差分画像Ｄ２を代表例として拡大して示している。
【００４３】
差分画像Ｄ１、Ｄ２に上下方向の線状ＲＯＩ（Ｌ１）を設定した時、その線上にある画素データに対する画素値をプロットすると、図７のＥ１，Ｅ２で示すプロファイルカーブが得られる。プロファイルカーブＥ１とＥ２はフレームの進みによって正・負の特性が逆になるため、プロファイルカーブＥ１，Ｅ２の特性から呼吸位相（Ｆ）を判断し、呼吸位相値（％）を算出する。
【００４４】
呼吸位相値（％）は図５（Ｇ）と同様に算出し、臓器が最も上方向に存在するタイミングを０％（１００％）とし、呼吸位相Ｆが上向きの期間は１００に向かって徐々に数値が増加し、呼吸位相Ｆが下向きの期間は０％を起点にしてゼロから徐々に数値が増加する。このようにプロファイルカーブＥ１，Ｅ２の正・負の特性から呼吸位相Ｆを判断することで呼吸位相値（％）を簡便に算出することができる。
【００４５】
尚、プロファイルカーブを検出する線状ＲＯＩの設定方法としては、ポストプロセスで適切な場所に操作者が指定する方法や、線状ＲＯＩをプリセットしておく方法を採用するとよい。またプリセットされた複数の線状ＲＯＩを撮影のパラメータとし、撮影部位に応じて線状ＲＯＩを自動設定する方法や、プリセットされた複数の線状ＲＯＩを使用したときの処理結果を比較して適正なものを選択する方法を採用してもよい。さらに画像の周期的な動き部分を検出してダイナミックに線状ＲＯＩを自動決定する方法を採用してもよい。
【００４６】
また以上の説明では、呼吸による横隔膜の移動は頭尾方向であり、移動方向の検出は主に上下方向として考えているため、上下方向の線状ＲＯＩを例示したが、平面あるいは２方向の線によるＲＯＩにより左右の動きを検出するようにしても良い。 
また本実施形態では、ＤＳＡ撮影後のポストプロセスにおいて呼吸位相を特定する例を述べたが、マスク画像群の収集後にマスク画像群の各マスク画像の呼吸位相を特定してコントラスト画像の収集に移行し、コントラスト画像に対しては差分画像に基いて呼吸位相の特定をリアルタイムに行うようにしてもよい。これにより、ＤＳＡ撮影中のリアルタイム・リマスキングへの応用も可能となる。
【００４７】
図８は、呼吸位相画像演算部３６の動作を説明するフローチャートである。ステップＳ１では、マスク画像群及びコントラスト画像群のそれぞれについて、隣り合うフレームの差分画像Ｄを生成して保管する。ステップＳ２では例えば線状ＲＯＩを設定し、差分画像のＲＯＩの位置のプロファイルカーブＥを作成し、プロファイルカーブを数値情報として保管する。
【００４８】
ステップＳ３では、保管されたプロファイルカーブの分布と、隣り合うプロファイルカーブの分布と位置の変化とから移動方向を特定し、呼吸位相Ｆを検出する。 ステップＳ４では呼吸位相Ｆをもとに移動方向が反転するフレームを検出し、そのフレームを基準に呼吸周期を決定し、各フレームの呼吸位相値％Ｇを算出する。
【００４９】
ステップＳ５では、マスク画像群とコントラスト画像群のうち、呼吸位相値が最も近似するマスク画像及びコントラスト画像を選択し、サブトラクション時に使用すべき両画像のフレームを特定し、フレーム番号Ｈをフレームテーブルに保管する。ステップＳ６では、フレーム番号Ｈをインデックスとしてサブトラクションに使用するマスク画像とコントラスト画像を取り出し、減算（サブトラクト）処理して、サブトラクション画像を得て表示部３９に表示する。
【００５０】
以上述べたように本発明のＸ線診断装置では、自然呼吸下におけるＤＳＡ撮影で、アーチファクトの少ないサブトラクション画像を自動的に表示することができる。また煩雑なリマスク処理を繰り返す必要がないため、診断効率を向上することができる。また呼吸位相を検出するため特別なセンサ等の器具を使用する必要もない。
【００５１】
尚、以上の説明では、Ｃアームを用いたＸ線診断装置（アンギオ装置）を例に説明したがＸ線ＣＴ装置等の他のＸ線診断装置にも適用することができる。また、特許請求の範囲を逸脱しない範囲で種々の変形が可能である。
【符号の説明】
【００５２】
１００…Ｘ線診断装置
１０…Ｘ線発生部
１１…Ｘ線管
２０…Ｘ線検出部
２１…Ｘ線検出器
２２…データ変換部
３０…コンピュータシステム
３１…システム制御部
３２…操作部
３３…ユーザインターフェース
３４…画像処理部
３５…記憶部
３６…呼吸位相画像演算部
３６１…差分画像生成部
３６２…サブトラクション演算処理部
３７…画像データ処理部
３８…Ｄ／Ａ変換部
３９…表示部
４０…高電圧制御部
４１…Ｃアーム
４５…天板
【図１】
0005523782.tif 000002
【図２】
0005523782.tif 000003
【図３】
0005523782.tif 000004
【図４】
0005523782.tif 000005
【図５】
0005523782.tif 000006
【図６】
0005523782.tif 000007
【図７】
0005523782.tif 000008
【図８】
0005523782.tif 000009