画像診断装置、画像処理装置、及び画像処理プログラムを記録した記録媒体
【課題】 造影剤注入前後において体動が発生した場合であっても、アーチファクトがなく臨床的価値の高いＤＳＡ画像を生成する画像診断装置等を提供すること。
【解決手段】 被検体への造影剤の注入の前後にわたって連続的に収集された撮影対象部位に関する複数の画像を、マスク像とコントラスト像とに区分して記憶する記憶装置１５と、各コントラスト像とのサブトラクション処理に利用する最適マスク画像を、複数のマスク像から選択するマスク選択部１９と、最適画像を部分的に変形させるワーピング処理を実行することで、コントラスト画像に位置合わせされたワーピング画像を生成するワーピング処理部２２と、各コントラスト像からワーピング画像をサブトラクションすることでワーピングＤＳＡ画像を生成するサブトラクション処理部２１と、を具備する。
【解決手段】 被検体への造影剤の注入の前後にわたって連続的に収集された撮影対象部位に関する複数の画像を、マスク像とコントラスト像とに区分して記憶する記憶装置１５と、各コントラスト像とのサブトラクション処理に利用する最適マスク画像を、複数のマスク像から選択するマスク選択部１９と、最適画像を部分的に変形させるワーピング処理を実行することで、コントラスト画像に位置合わせされたワーピング画像を生成するワーピング処理部２２と、各コントラスト像からワーピング画像をサブトラクションすることでワーピングＤＳＡ画像を生成するサブトラクション処理部２１と、を具備する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、例えば ディジタルアンギオ装置においてＤＳＡ画像を作成、保存、表示するＸ線診断装置等に代表される画像診断装置、画像処理装置、及び画像処理プログラムを記録した記録媒体に関する。
【背景技術】
【０００２】
造影剤を注入する前とか、観察領域に造影剤がない状態で収集した画像（マスク像）と、造影剤の注入後で観察領域に造影剤が存在する状態で収集した血管コントラストが増強された画像（コントラスト画像）との間でサブトラクション処理を行い、そのサブトラクション像を用いて行う診断は、例えば血管の狭窄位置や状態を観察するのに非常に有用であり、今日では多くの医療現場において使用されている。
【０００３】
このサブトラクション処理では、サブトラクション処理に供する画像間で位置や濃度の変化がない構造物は消去され、位置や濃度の変化のあった構造物は残存する。従って、サブトラクション処理によれば、造影前後の期間で被検体の体動がなければ、造影前後で濃度変化が生じる血管のみが抽出されることになるので、背景構造物に影響を受けることなしに血管狭窄等の観察を行うことができる。
【０００４】
ところが、このサブトラクションに供される各画像を収集する間に被検体が複雑に動くことが多く、その場合は各画像間で位置ずれが生じ偽像として映出されてしまう（なお、この様な偽像を「アーチファクト」と呼ぶ。）。従来は、例えば特許文献１に示すように、このアーチファクトを抑制するために、サブトラクションに供する前にいずれか一方の画像を求められた位置ずれ量に基づいて平行移動させて位置合せを行うことを行っていた。
【特許文献１】特開昭６４−２０８３０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、これらタイミングの間で被検体の心拍や呼吸或いは咳などによって体動が生じた場合、被検体の構造物に位置の変化が生じてしまう。特に、高齢者や重症者、意識がない被検者では息止めができないため、こういった被検者で特に腹部の検査においてはこの位置変化が甚だしい。この様に位置変化が生じた前後の画像を用いてサブトラクション処理が行われると、造影された血管が映出されるのはもちろんであるが、位置の変化が生じて完全に差分されず残存した他の構造物もアーチファクトとして映出されてしまう。実際の臨床の場面では、時々咳が止まらなくなったり、造影剤注入に熱さを感じて条件反射的に動いてしまう被検体もある。この様な場合には、横隔膜の移動によるアーチファクトが生じてしまい診断に支障をきたすことになる。
【０００６】
また、このような事情の他に心臓や大血管の拍動、体動、腸管の蠕動等の要因により、サブトラクションに供する２つの画像における被写体形状は完全同一ではないため、特許文献１に示すように平行移動による位置合せを行ってもＤＳＡ画像上には多くのアーチファクトが発生することがあり、医師による診断に悪影響を与える可能性がある。
【０００７】
本発明は、上記事情を鑑みてなされたもので、マスク像やコントラスト像の撮影中において、又はマスク像撮影からコントラスト像撮影までの間に体動が発生した場合であっても、アーチファクトがなく臨床的価値の高いＤＳＡ画像を生成することができる画像診断装置、画像処理装置、及び画像処理プログラムを記録したコンピュータ読み取り可能な記録媒体を提供することを目的としている。
【課題を解決するための手段】
【０００８】
本発明は、上記目的を達成するため、次のような手段を講じている。
【０００９】
請求項１に記載の発明は、被検体の撮影対象部位に関する複数の画像を、当該被検体への造影剤の注入の前後にわたって連続的に収集する画像収集手段と、前記複数の画像のうちの一画像である第１の画像を少なくとも部分的に変形させる変形処理を実行することで、前記複数の画像のうちの前記第１の画像とは異なる第２の画像に表示された各種部位との位置合わせが施された第３の画像を生成する変形処理手段と、前記各第２の画像と前記第３の画像の間でサブトラクションすることで、第１のサブトラクション画像を生成する画像生成手段と、を具備することを特徴とする画像診断装置である。
【００１０】
請求項１１に記載の発明は、被検体への造影剤の注入の前後にわたって連続的に収集された、当該被検体の撮影対象部位に関する複数の画像を記憶する記憶手段と、前記複数の画像のうちの一画像である第１の画像を少なくとも部分的に変形させる変形処理を実行することで、前記複数の画像のうちの前記第１の画像とは異なる第２の画像に表示された各種部位との位置合わせが施された第３の画像を生成する変形処理手段と、前記各第２の画像と前記第３の画像の間でサブトラクションすることで、第１のサブトラクション画像を生成する画像生成手段と、を具備することを特徴とする画像処理装置である。
【００１１】
請求項１６に記載の発明は、コンピュータに、被検体への造影剤の注入の前後にわたって連続的に収集された、当該被検体の撮影対象部位に関する複数の画像のうちの一画像である第１の画像を少なくとも部分的に変形させる変形処理を実行することで、前記複数の画像のうちの前記第１の画像とは異なる第２の画像に表示された各種部位との位置合わせが施された第３の画像を生成する第１の機能と、前記各第２の画像と前記第３の画像の間でサブトラクションすることで、第１のサブトラクション画像を生成する第２の機能と、を実現させるための画像処理プログラムを記録したコンピュータ読み取り可能な記録媒体である。
【発明の効果】
【００１２】
以上本発明によれば、マスク像やコントラスト像の撮影中において、又はマスク像撮影からコントラスト像撮影までの間に体動が発生した場合であっても、アーチファクトがなく臨床的価値の高いＤＳＡ画像を生成することができる画像診断装置、画像処理装置、及び画像処理プログラムを記録した記録媒体を実現することができる。
【発明を実施するための最良の形態】
【００１３】
以下、本発明の実施形態を図面に従って説明する。以下の説明において、略同一の機能及び構成を有する構成要素については、同一符号を付し、重複説明は必要な場合にのみ行う。
【００１４】
なお、本発明は造影画像を収集可能な画像診断装置を含んでいる。この種の画像診断装置としては、Ｘ線診断装置、超音波診断装置、磁気共鳴イメージング装置、さらにＰＥＴ，ＳＰＥＣＴが含まれ、そのいずれでも本発明に適用可能である。以下では、Ｘ線診断装置を例に説明する。
【００１５】
図１は本実施形態に係るＸ線診断装置の構成を示している。本Ｘ線診断装置は、Ｘ線管１とＸ線検出器３とを有する画像収集部２、インジェクタ５、画像処理装置７を具備している。
【００１６】
画像収集部２は、寝台に載置された被検体Ｐを挟んで対向配置されたＸ線管１とＸ線検出器３とを有する。
【００１７】
Ｘ線管１は、Ｘ線を発生する真空管であり、高電圧発生装置からの高電圧により電子を加速させ、ターゲットに衝突させることでＸ線を発生する。発生したＸ線は、被検体に対して曝射される。
【００１８】
Ｘ線検出器３は、シンチレータとフォトダイオードアレイとを有するフラットパネルデテクタ（ＦＰＤ）が採用される。しかし、Ｘ線検出器３は、イメージインテンシファイアとＴＶカメラとの組み合わせが採用されても良い。なお、ＦＰＤは、被検体を透過したＸ線を光電膜に当てることで電子正孔を生成し、これを半導体スイッチにおいて蓄積し、電気信号として読み出すことでＸ線信号を検出するものである。また、Ｉ．Ｉ．は、被検者を透過したＸ線を入力蛍光面で光学情報に変換し、さらに当該光学情報に基づいて生成された光電子を出力蛍光面に衝突させて高輝度光学画像を生成する。このＸ線検出器３から出力されるアナログの映像信号は、アナログディジタル変換器１３から画像処理装置７に供給される。
【００１９】
インジェクタ５は、被検体Ｐに造影剤を自動注入するための装置である。造影剤が被検体Ｐに注入された時、トリガ信号又はその注入時刻を表す信号が、インジェクタインタフェース１１を介して画像処理ユニット７に供給される。
【００２０】
画像処理装置７は、インジェクタインタフェース１１、アナログディジタル変換器１３とともに、ＣＰＵ９、画像データや各種制御データ等を記憶する記憶装置１５、コンソール（操作卓）１７、マスク選択部１９、サブトラクション処理部２１、ワーピング処理部２２、ディジタルアナログ変換器２３、ディスプレイ２５から構成される。
【００２１】
ＣＰＵ９は、画像データの収集に関する制御、及び収集した画像データの画像処理、画像再生処理等に関する制御を行う中央処理装置である。
【００２２】
インジェクタインタフェース１１は、インジェクタ５とＣＰＵ９との間のデータ通信を可能にするための装置である。
【００２３】
アナログディジタル変換器１３は、Ｘ線検出器３から入力したアナログ信号を、ディジタル信号に変換し、フレーム毎の画像データとして記憶装置１５、マスク選択部１６等に送り出す。
【００２４】
記憶装置１５は、取得された画像データ、画像取得及び画像処理のための各種制御データ等を記憶する。また、記憶装置１５は、ＣＰＵ９による制御のもと、造影剤注入の前後にわたって収集された複数の画像を、後述する遅延時間に基づいて、マスク像とコントラスト像とに区分して記憶する。さらに、記憶装置１５は、後述するレビュー機能を実現するために、マスク像、コントラスト像、ワーピングＤＳＡ処理によって生成されたワーピング画像等を記憶する。
【００２５】
コンソール（操作卓）１７は、キーボードや各種スイッチ、マウス等を備えた入力装置であり、撮影指示や画像選択等を入力する際に使用される。
【００２６】
マスク選択部１９は、造影剤注入後に収集した各コントラスト像について、後述する遅延時間内に収集した複数のマスク像からＤＳＡ処理に最適な１枚を最適マスク像として選択する。その選択手法は種々のものを採用することができる。本実施形態では、後述するように、撮影プログラムに応じていずれの選択手法を採用するかを決定する。
【００２７】
サブトラクション処理部２１は、各コントラスト像と、コントラスト像毎に選択されたマスク像とを用いてサブトラクション処理を実行し、ＤＳＡ画像を生成する。また、サブトラクション処理部２１は、ワーピング処理部２２において最適マスク像に基づいて生成されたワーピング画像（後述）と、各コントラスト像とを用いてサブトラクション処理を実行し、ワーピングＤＳＡ画像を生成する。
【００２８】
ワーピング処理部２２は、サブトラクション処理に用いられるコントラスト像と最適マスク像との間で形状が異なる部位を、最適マスク画像において画素単位で部分的に変形（ワーピング）させることで位置合わせを行うワーピング処理を実行する。このワーピング処理、及び当該ワーピング処理を前提とするワーピングＤＳＡ処理については、後で詳しく説明する。
【００２９】
ディジタルアナログ変換器２３は、サブトラクション処理部２１等から入力した画像データのディジタル信号列を、アナログ信号列に変換する。
【００３０】
ディスプレイ２５は、ディジタルアナログ変換器２３から受け取った信号により、Ｘ線診断画像等を表示するＣＲＴ、プラズマディスプレイ、液晶ディスプレイ等の表示装置である。特に、本ディスプレイは医用機器に使用されるものであるため、走査線１０００本以上の高精細ディスプレイであることが好ましい。
【００３１】
（ワーピングＤＳＡ処理）
本Ｘ線診断装置は、被検体の体動に起因して発生するコントラスト像に対するマスク像の位置ずれを、ワーピングを利用して修正した後にＤＳＡ処理を実行する処理（ワーピングＤＳＡ処理）を有している。このワーピングＤＳＡ処理により、マスク像やコントラスト像の撮影中において、又はマスク像撮影からコントラスト像撮影までの間に体動が発生した場合であっても、アーチファクトがなく臨床的価値の高いＤＳＡ画像を生成することができる。以下、この処理について詳しく説明する。
【００３２】
なお、被検体の体動は、周期的なものと非周期的なものとのに分類することができる。周期的な体動とは、心臓や大血管の拍動、呼吸による横隔膜や肺の部位の動き等である。また、非周期的な体動とは、腸管の蠕動、撮影中の咳、造影剤流入を体感した場合における条件反射等である。本ワーピングＤＳＡ処理は、いずれの体動をも対象とすることができる。
【００３３】
図２は、本ワーピングＤＳＡ処理の流れを示したフローチャートである。同図に示すように、本ワーピングＤＳＡ処理の各ステップでは、人体構造の解析、画像濃度の補正、ローカルマッチング、画像のワーピング、差分画像の生成の各処理が実行される。
【００３４】
まず、人体構造の解析処理（ステップＳａ）では、撮影したマスク像に基づいて人体構造物を解釈し、人体解析データを作成する。ここで、人体解析データとは、例えば胸部の血管撮影の場合であれば、横隔膜、脊椎、心臓など、マスク像上に写っているランドマークとなる人体構造物の位置、形状に関する情報である。
【００３５】
図３は、腹部アンギオ画像においてランドマークとする人体構造物（横隔膜、脊椎、心臓、肋骨）の一例を示した図である。また、図４は、胸部アンギオ画像においてランドマークとする人体構造物（胸郭、横隔膜、脊椎）の一例を示した図である。
【００３６】
この人体構造の解析処理においては、複数枚のマスク像を用いて平均化等された情報に基づいて、人体解析データを作成することも可能である。この場合、同一の撮影シーケンスにおいて撮影されたマスク像に限定されない。例えば、同一患者の同一検査のアンギオ画像であれば、それほど大きな変動はしないことが予想されるから、過去に撮影された同一患者の同一検査のマスク像や人体解析データも利用することができる。特に、骨、肺、腸管ガス、体外の空気等の高コントラストを形成する部位等、また、差分に用いる２画像の位相差を類推する目的で血管内に挿入された、コイル、ステント、ＩＶＣフィルター、カテーテルのマーカーや過去の手術で体内に挿入されたクリップ、義歯などは、人体解析データとして利用可能である。また、後述するワーピング処理によって体動補正をより効率よく行うために、前もって被検者の体表に付着させ金属等のマーカーを人体解析データとして利用することも可能である。
【００３７】
なお、一般に、撮影対象部位（例えば、頭部、頚部、胸部、腹部、四肢）によって、ランドマークとなる人体構造物は異なる。従って、人体解析データを作成する対象部位を予め指定すれば、人体構造物の解釈がより容易に行うことができる。
【００３８】
また、同一の撮影対象部位であっても、正面、側面、斜位といった撮影方向によって、ランドマークとなる人体構造物は異なる。本解析処理によって得られた人体解析データは、同一の撮影対象部位であれば、いかなる撮影方向にも適応させることが可能である。
【００３９】
この様に生成される人体解析データは、後段のステップにおいて、人体構造の知識（解剖学的情報）に沿った処理ができるように利用される。
【００４０】
次に、画像濃度の補正（ステップＳｂ）では、コントラスト像、マスク像（今の場合、最適マスク像）に対して、画像濃度の補正処理を適用して、補正コントラスト像、補正マスク像を作成する。この補正処理により、腸管ガスの影響除去、ノイズ除去、アンギオ装置の撮影時のフィルター除去など、造影剤注入以外の原因に起因する、ＤＳＡ画像画質低下の要因を低減させることができる。
【００４１】
次に、ローカルマッチング（ステップＳｃ）では、ステップ２で生成された補正コントラスト像と補正マスク像について、両画像が局所的に一致するような関心領域を多数配置し、これらに基づいて移動ベクトルを導出する。具体的なローカルマッチングの手法は、次の様である。
【００４２】
図５（ａ）、図５（ｂ）は、ローカルマッチングの手法の一例を説明するための図である。図５（ａ）に示すように、まず、ある一定サイズのテンプレートＲＯＩ５０を補正コントラスト像上に配置し、当該テンプレートＲＯＩ５０に対応する（複数個の場合には、それぞれに対応する）探索領域ＲＯＩ５４を補正マスク像上に配置する。そして、探索領域ＲＯＩ５４とテンプレートＲＯＩ５０との間で、テンプレートマッチングと呼ばれる技法を用いて、局所的な移動量（移動ベクトル）を導出する。なお、一般的には、テンプレートＲＯＩが複数設定されるため、移動ベクトルは、個々のテンプレートＲＯＩに対応して導出される。
【００４３】
なお、移動ベクトルを導出するためのテンプレートＲＯＩ５０及び探索領域ＲＯＩ５４は、ステップＳａにおいて取得した人体解析データに基づいて決定された各画像の所定エリア内に、例えば一定の間隔をおいて配置される。
【００４４】
図６は、補正コントラスト像又は補正マスク像としての腹部アンギオ画像上において、複数のテンプレートＲＯＩ５０又は複数の探索領域ＲＯＩ５４を配置するエリアの一例を示した図である。同図に示す腹部アンギオ画像の例では、テンプレートＲＯＩ５０又は探索領域ＲＯＩ５３を配置しないのは、脊椎のような不動部分や、横隔膜より上の部分などである。
【００４５】
次に、画像のワーピング（ステップＳｄ）では、ステップ２で導出された補正マスク像について、ステップＳｃで導出された移動ベクトルを基に、例えば画素単位の移動ベクトルを導出し、これに基づいて補正マスク像上に表示された所定部位を変形（ワープ）させるワーピング処理を実行する（以下、このワーピング処理によって得られる画像を「ワーピング画像」と呼ぶ。）。
【００４６】
この様なワーピング処理を行うのは、次の理由による。すなわち、補正マスク像と補正コントラスト像に写っている人体構造物は、ほぼ同じ形状で写っているものと考えてよい。しかしながら、呼吸、体動、心臓や大血管の拍動、腸管の蠕動により変形が発生する部位については、ＤＳＡ対象となる補正マスク像と補正コントラスト像との間で形状が異なる場合がある。この様な変形が発生する部位については、補正マスク像を画素単位で部分的に変形させ（画像ワーピング）、補正コントラスト像に合せこむことによって、その変形の効果を相殺する。これにより、ＤＳＡ画像におけるアーチファクトの原因となる補正マスク像と補正コントラスト像との位置ずれをなくすことができる。
【００４７】
なお、ローカルマッチング、及びワーピング処理のアルゴリズムについては、例えば特開平７−３７０７４号公報に記載されている経時的差分画像作成と同じアルゴリズムを適用することができる。
【００４８】
次に、差分画像の生成（ステップＳｅ）では、ステップＳｂにおいて生成された補正コントラスト像から、ステップＳｄにおいて生成されたワーピング画像を差し引くことでＤＳＡ画像を生成し、表示部に出力する。
【００４９】
（撮影シーケンス）
次に、上記ワーピングＤＳＡ処理を含む一連の撮影シーケンスについて説明する。以下説明を簡単にするため、被検体に対して撮影系を移動させない静的（static）なＤＳＡ撮影を対象とする。しかしながら、本発明は、ステッピングＤＳＡ撮影、回転ＤＳＡ撮影といった動的な撮影手法においても適用可能である。
【００５０】
図７は、当該撮影シーケンスにおいて実行される処理の流れを示したフローチャートである。同図に示すように、まず、患者ＩＤ、撮影プログラム等の入力・選択等が実行される（ステップＳ１）。ここで、撮影プログラムとは、撮影領域（部位）別に固有の遅延時間（マスク像の撮影開始からコントラスト像の撮影開始までの期間）、撮影条件、後処理等を実行するための制御プログラムである。本実施形態においては、「腹部専用プログラム」、「心臓専用プログラム」、「その他の領域（例えば、頭部、四肢等）用プログラム」等の複数種の撮影プログラムが選択可能であるとする。
【００５１】
次に、ステップＳ１において選択された撮影プログラムに合わせて、管電流／管電圧等、遅延時間（マスク像の撮影開始からコントラスト像の撮影開始までの期間）の設定を実行する（ステップＳ２）。なお、遅延時間は、確実に一呼吸分のマスク像を収集するようにするため、一呼吸期間よりも少し長めに設定するのが一般的である。
【００５２】
図８は、ステップＳ２において利用される入力画面の一例を示した図である。同図において、例えばステップＳ１において腹部専用プログラムが選択された場合には、１呼吸期間に多少のマージンを考慮して、最初の数秒間は造影剤を注入しないように遅延時間を設定する。また、ワーピングＤＳＡ処理についてより好適なマスク像を選択可能とするために、遅延時間内の画像収集フレームレートは、遅延時間経過後の画像収集フレームレートより高くなるように、例えば遅延時間内の画像収集フレームレートを１５／３０ｆｐｓ、遅延時間経過後の画像収集フレームレートを７．５／３ｆｐｓとして設定する。
【００５３】
また、例えばステップＳ１において心臓専用プログラムが選択された場合には、呼吸動より速い心拍動を考慮して、最初の数秒間は造影剤を注入しないように遅延時間を設定する。また、同様に、遅延時間内の画像収集フレームレートが遅延時間経過後の画像収集フレームレートより高くなるように、例えば遅延時間内の画像収集フレームレートを３０／６０ｆｐｓ、遅延時間経過後の画像収集フレームレートを１５／３０ｆｐｓとして設定する。
【００５４】
なお、ステップＳ１においてその他の領域用プログラムが選択された場合には、本実施形態においては、遅延時間、フレームレート等の設定には特に条件を設けない。しかしながら、これに限定する趣旨ではなく、その他の領域用プログラムが選択された場合であっても、必要に応じて領域毎の固有の遅延時間、フレームレート等を設定可能な構成にしてもよい。また、当該ステップＳ２における遅延時間、フレームレート等の設定は、ステップＳ１において選択された撮影プログラムに応じて、自動的に設定される構成であってもよい。
【００５５】
次に、専用ＤＳＡ処理プログラム実行の実行／非実行を選択する（ステップＳ３）。ここで、専用ＤＳＡ処理プログラムとは、複数枚のマスク像の中から自動的にＤＳＡ処理に最適なマスク像（最適マスク像）を自動的に選択し、これを用いてワーピングＤＳＡ処理を実行するものである。
【００５６】
次に、複数枚のマスク像撮影及び複数枚のコントラスト像撮影が実行される（ステップＳ４）。このとき、コントラスト像の最初の撮影は、マスク像の最初の撮影から遅延時間経過後に開始される。
【００５７】
次に、ステップＳ３において専用ＤＳＡ処理プログラムが選択された場合には、当該プログラムに従って最適マスク像の自動選択処理が実行される（ステップＳ５ａ）。この自動選択処理としては、例えば特開２００４−１１２４６９号公報に開示されている手法を採用することができる。この手法では、各マスク像と各コントラスト像との「位置ずれの程度の類似性を示す一致度」を計算することで、位置ずれの程度が最も近いマスク像（最適マスク像）をコントラスト像毎に選択し、これを用いてＤＳＡ処理を実行するものである。これにより、例えば体動があった領域について、コントラスト像に最も近い形状を有するマスク像を選択することができる。
【００５８】
なお、最適マスク像は必ずしもコントラスト像毎に選択する必要はなく、複数のコントラスト像に共通する最適マスク像（共通最適マスク像）を選択する構成であってもよい。本実施形態では、説明を簡単にするため、共通最適マスク像を選択するものであるとする。
【００５９】
また、自動選択処理は、上記手法に限定されない。例えば、心臓専用プログラムを実行する場合には、心電波形の位相が略一致するフレーム同士を選択しサブトラクションする構成であることが好ましい。また、その他の自動選択処理として、心臓や大血管の拍動では指先でモニターする脈波を、呼吸では被検者に装着した呼吸同期装置を、体軸の移動では被検者の皮膚に塗布したマーカーをそれぞれ用いることも可能である。
【００６０】
次に、ステップＳ５ａにおいて自動選択された最適マスク像と各コントラスト像とを用いて、ワーピングＤＳＡ処理を実行する（ステップＳ６ａ）。ワーピングＤＳＡ処理の内容は、図２のステップＳａ〜Ｓｅに示した通りである、当該ワーピングＤＳＡ処理によって得られたＤＳＡ画像は、ディスプレイ２５に表示される。
【００６１】
なお、当該ワーピングＤＳＡ処理中においては、ＤＳＡ画像又はコントラスト像と（ワーピング処理前の）最適マスク像とＤＳＡ画像を、中間処理画像として表示することが好ましい。ワーピングＤＳＡ画像が生成されると、ディスプレイ２５上は、自動的に中間処理画像からワーピングＤＳＡ画像に切り換えて表示する。この様な構成によれば、操作者はワーピングＤＳＡ処理中であっても診断を行うことができ、また、中間処理画像とワーピングＤＳＡ画像との双方を臨床情報として有効に利用することができる。
【００６２】
一方、ステップＳ３において専用ＤＳＡ処理プログラムが選択されなかった場合には、ステップ５ｂ〜６ｅまでの流れに従って処理される。すなわち、ＤＳＡ処理用として初期設定されたマスク像（例えば、最初に撮影されたマスク像）が最適マスク像として選択される（ステップＳ５ｂ）。
【００６３】
次に、例えば図９に示すような操作画面から「Ｗａｒｐｉｎｇ」のボタンがクリックされると、初期設定による最適マスク像と現在表示中のコントラスト像（カレント像）とを用いて、ワーピングＤＳＡ処理が実行される（ステップＳ６ｂ）。当該ワーピングＤＳＡ処理によって得られたＤＳＡ画像は、ディスプレイ２５に表示される。
【００６４】
次に、操作者は、表示されたＤＳＡ画像を観察し、リマスク処理を実行する必要があるか否かを判定し（ステップＳ６ｃ）、実行する場合には、最適マスク像の変更が実行される（ステップ６ｅ）。ここで、「リマスク処理」とは、最適マスク像の変更を行うための処理であり、例えば次の様なものがある。
【００６５】
図１０は、リマスク処理の一例を説明するための図である。同図に示すように、初期設定により「収集画像１」が全コントラスト像に対して共通する最適マスク像として設定されているとする。この状態で、図９に示した入力画面において「Ｒｅｍａｓｋ」のボタンをクリックすると、他の任意の収集画像を最適マスク像として再設定することができる。実際の撮影においては、造影剤注入の際に患者が「熱い」と感じる場合があり、このときに体動が発生する場合が多い。従って、造影剤注入が終了し且つ造影剤が撮影対象に流入する手前である画像（例えば、「収集画像８」）を新たな最適マスク像として設定することで、好適なＤＳＡ処理を実現することができる。また、ステップＳ５ａではあらかじめマスク画像を一致度の高い順に並べ替え、リマスクの候補の選択に用いてもよい。なお、図１０に示した例は、例えば静的なＤＳＡ撮影、ステッピングＤＳＡ撮影等において、特に有益である。
【００６６】
また、図１１は、回転ＤＳＡ撮影におけるリマスク処理の例を説明するための図である。回転ＤＳＡ撮影では、初期設定として、同じ角度で収集したマスク像ｉ（ｉ＝１〜ｎ）とコントラスト像ｉとの間でサブトラクション処理が実行される。この状態で、図９に示した入力画面において「Ｒｅｍａｓｋ」のボタンをクリックすることで、所定のコントラスト像に対するマスク像だけを変更することができる。なお、図１１では、「コントラスト像１」に対する最適マスク像を「マスク像１」から「マスク像３」に変更する例を示してある。
【００６７】
また、図１２は、回転ＤＳＡ撮影におけるリマスク処理の他の例を説明するための図である。同図に示すように、リマスク処理によって新たに選択したマスク像を基準に、全てのマスク像をコントラスト像に対して一定フレーム数ずらしてサブトラクション処理を実行することも可能である。
【００６８】
ステップＳ６ｅにおいて最適マスク像が変更されると、操作画面からの「Ｗａｒｐｉｎｇ」ボタンの操作に応答して、リマスク処理後の最適マスク像とカレント像とを用いて、ワーピングＤＳＡ処理が実行され、得られたＤＳＡ画像がディスプレイ２５に表示される。操作者は、表示されたＤＳＡ画像を観察し、再度リマスク処理を実行する必要があるか否かを判定し、必要ならば最終的な最適マスク像を決定するまでリマスク処理を繰り返す。最終的な最適マスク像は、例えば図９に示す「ｃｏｎｆｉｒｍ」ボタンをクリックすることで決定され、以降のコントラスト像に対するワーピングＤＳＡ処理は、当該最適マスク像を用いて実行される。一方、「ｃｏｎｆｉｒｍ」ボタンをクリックしなければ、当該最適マスク像は、カレント像に対してのみ適用されることになる。
【００６９】
次に、以降のコントラスト像に対して、決定された最適マスク像を用いてワーピングＤＳＡ処理を実行する（ステップＳ６ｄ）。このとき、例えば図９に示す「ａｐｐｌｙ」ボタンをクリックすれば、マルチフレーム画像に対して一度にワーピングＤＳＡ処理を実行することができる。
【００７０】
なお、ステップＳ３において専用ＤＳＡ処理プログラムが選択された場合であっても、必要に応じてワーピングを適用しないＤＳＡ処理を実行することも可能である。すなわち、専用ＤＳＡ処理プログラムが選択された場合には、図９に示す「Ｗａｒｐｉｎｇ」ボタンは初期設定でＯＮ状態となっている。また、専用ＤＳＡ処理プログラム選択後に撮影された画像は、初期設定でＣｈｅｃｋ−ＢｏｘがＯＮ状態となっており、また、「Ｗａｒｐｉｎｇ」ボタンがＯＮ状態となっている時にＣｈｅｃｋ−Ｂｏｘが選択可能となる。このＣｈｅｃｋ−ＢｏｘをＯＦＦ状態とした画像については、ワーピング処理前の補正マスク像が最適マスク像として選択することができる。
【００７１】
以上述べた処理によって得られるＤＳＡ画像等は、自動的に記憶装置１５に記憶される（ステップＳ７）。
【００７２】
（レビュー機能）
次に、本Ｘ線診断装置が有するレビュー機能について説明する。本機能は、ワーピングＤＳＡ処理によって得られた画像（例えば、図７のステップＳ６ａ、Ｓ６ｂ、Ｓ６ｄにおいて得られた画像）を事後的に再生するものである。
【００７３】
まず、図７のステップＳ７において各コントラスト像と各マスク像、及び各コントラスト像に対応するワーピング画像を自動的に記憶装置１５に保存する。この様に各コントラスト像に対応するワーピング画像を保存することで、ワーピング画像生成を省略したワーピングＤＳＡ処理を事後的に実行でき、ワーピングＤＳＡ画像を迅速に再生することができる。なお、メモリ資源をできるだけ有効に活用する観点から、保存するワーピング画像については非可逆圧縮を施す構成であってもよい。
【００７４】
また、ワーピング画像そのものを保存する構成ではなく、各コントラスト像に対応する最適マスク像のフレーム番号を保存し、ワーピング画像生成を含むワーピングＤＳＡ処理を事後的に実行することで、ワーピングＤＳＡ画像を再生することも可能である。この場合、最適マスク像からワーピング画像を生成する過程を保存しておくことで、ワーピング画像そのものを保存するよりもデータ量を小さくし、且つ最適マスク像のフレーム番号を保存するよりも事後再生において迅速化を図る構成としてもよい。また、この様にワーピング画像生成過程を保存する場合には、そのログの解析を元に、作成されたＤＳＡ画像の画質の推測をし、当該ＤＳＡ画像と供に表示することも可能である。
【００７５】
以上述べた構成によれば、以下の効果を得ることができる。
【００７６】
本Ｘ線診断装置においては、最適マスク像からワーピング画像を生成し、これとコントラスト像とを用いてＤＳＡ処理を実行する。従って、最適マスク像とコントラスト像との単純差分によりＤＳＡ処理を実行する従来技術に比して、より位置精度の高いサブトラクション処理を実行することができる。その結果、マスク像やコントラスト像の撮影中において体動が発生した場合であっても、従来よりもアーチファクトがなく臨床的価値の高いＤＳＡ画像を提供することができる。特に、ワーピング処理は、部分的な位置ずれを補正するのに好適であり、マスク画像とコントラスト像とが完全に同位相になるような撮影は困難である場合には、特に実益がある。
【００７７】
また、本Ｘ線診断装置においては、撮影プログラムに応じて遅延時間内の画像収集フレームレートを設定することができる。従って、コントラスト像より多くのマスク像を収集することができ、ＤＳＡ処理においてより位置精度の高いマスク像を選択することができる。
【００７８】
また、本Ｘ線診断装置においては、複数のマスク像の中から最適マスク像を、例えば位置ずれの程度の類似性を示す一致度に基づいて、自動的に選択することができる。従って、コントラスト像と位置対応性の高い最適マスク像を迅速且つ簡便に選択することができる。
【００７９】
また、本Ｘ線診断装置においては、必要に応じて最適マスク像を他の任意のマスク像に変更することができる。従って、自動選択された最適マスク像、又は初期設定されている最適マスク像によるＤＳＡ画像以外にも、他のマスク像を最適マスク像とすることで、種々の所望のＤＳＡ画像を取得することができ、診断の自由度を広げることができる。
【００８０】
また、本Ｘ線診断装置においては、レビュー機能により、ワーピングＤＳＡ処理によって得られた画像を事後的に再生することができる。これにより、撮影以後の診断や、画像再生による撮影・診断の事後的な検討が可能となり、医療行為の質の向上に寄与することができる。たとえば、血管奇形、動静脈瘻、肝細胞癌などでは動脈塞栓術の前後で、動脈の閉塞や狭窄病変では人工血管によるバイパス術、血管形成術、ステント留置術前後で、それぞれの造影剤注入後の画像同士の差分画像を作成し、治療の効果を客観的に表示することも可能となる。
【００８１】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。具体的には、次のような変形例を挙げることができる。
【００８２】
（１）ワーピングにより十分補正できなかった腸管の蠕動などのアーチファクトは、例えばさらにダイナミックレンジ圧縮の技法を用い、低吸収値部分のコントラストを減弱させ、アーチファクトを軽減させる構成としてもよい。
【００８３】
（２）胸部の血管撮影、特に冠動脈撮影では、心拍動によるアーチファクトが問題となる。係る撮影においては、例えば初期設定として、心電図同期によるＲ波からのトリガー遅延を最適マスク像の選択に使用する構成であってもよい。
【００８４】
（３）血管撮影装置のＸ線管球と検出器を移動させながら撮影する３次元撮影の場合、もしくは被検者を載せたテーブルを移動させながら撮影を行う場合では、被検者とＸ線管球および検出器の空間的位置関係を、最適マスク像の選択に利用することもできる。
【００８５】
（４）一般の血管撮影装置に併設されている、ＣＴ画像にも本ワーピングＤＳＡ処理は適用可能である。この場合、ＩＶＲＣＴ（Interventional Radiology CT）で撮影された造影前と造影後のボリュームデータを２次元に投影したＲａｙＳｕｍ（もしくはＲａｙ−Ｓｕｍ）表示画像やＭＩＰ（Maximum Intensity Projection）画像等を用いて本技術を応用することができる。なお、これらＲａｙＳｕｍ表示画像等についての詳細は、特開２００２−２１９１２３公報に記載されている。
【００８６】
（５）上記実施形態においては、Ｘ線診断装置において、ワーピングＤＳＡ処理等を実施する例を示した。これに対し、Ｘ線診断装置と高速なインターフェースで連結され、例えばＤＩＣＯＭ（Digital Imaging and Communications in Medicine）による通信を行う専用のワークステーションによって、ワーピングＤＳＡ処理等を実施する構成であってもよい。ＰＡＣＳ（Picture archiving and communications system）のネットワーク上にワークステーションを配置して、血管撮影のデータを取得し任意に本手法を適応することも可能である。
【００８７】
なお、この場合、ＤＩＣＯＭのプライベートタグ中に遅延時間、画像収集レート等の情報を記録することで、ワークステーション上においても、ワーピングＤＳＡ処理を含む一連の撮影シーケンス（すなわち、図７に示したシーケンス）を擬似的に再現することができる。
【００８８】
（６）本実施形態に係るワーピングＤＳＡ処理等は、当該処理を実行するプログラムをワークステーション等のコンピュータにインストールし、これらをメモリ上で展開することによっても実現することができる。このとき、コンピュータに当該手法を実行させることのできるプログラムは、磁気ディスク（フロッピー（登録商標）ディスク、ハードディスクなど）、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤなど）、半導体メモリなどの記録媒体に格納して頒布することも可能である。
【００８９】
また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【産業上の利用可能性】
【００９０】
以上本発明によれば、マスク像やコントラスト像の撮影中において体動が発生した場合であっても、アーチファクトがなく臨床的価値の高いＤＳＡ画像を生成することができるＸ線診断装置、画像処理装置、及び画像処理プログラムを格納する記憶媒体を実現できる。
【図面の簡単な説明】
【００９１】
【図１】図１は本実施形態に係るＸ線診断装置の構成を示している。
【図２】図２は、本ワーピングＤＳＡ処理の流れを示したフローチャートである。
【図３】図３は、腹部アンギオ画像においてランドマークとする横隔膜、脊椎、心臓、肋骨の例を示した図である。
【図４】図４は、腹部アンギオ画像においてランドマークとする横隔膜、脊椎、心臓、肋骨を示した図である。
【図５】図５は、ローカルマッチングの手法の一例を説明するための図である。
【図６】図６は、腹部アンギオ画像上において、複数のテンプレートＲＯＩ５０又は複数の探索領域ＲＯＩ５３を配置するエリアの一例を示した図である。
【図７】図７は、ワーピングＤＳＡ処理を含む一連の撮影シーケンスにおいて実行される処理の流れを示したフローチャートである。
【図８】図８は、遅延時間等の入力画面の一例を示した図である。
【図９】図９は、ワーピング処理、リマスク処理を指示するための入力画面の一例を示した図である。
【図１０】図１０は、リマスク処理の一例を説明するための図である。
【図１１】図１１は、回転ＤＳＡ撮影においてリマスク処理の例を説明するための図である。
【図１２】図１２は、回転ＤＳＡ撮影におけるリマスク処理の他の例を説明するための図である。
【符号の説明】
【００９２】
１…Ｘ線管、２…画像収集部、３…Ｘ線検出器、５…インジェクタ、７…画像処理装置、９…ＣＰＵ、１１…インジェクタインタフェース、１３…アナログディジタル変換器、１５…記憶装置、１７…コンソール、１９…マスク選択部、２１…サブトラクション処理部、２２…ワーピング処理部、２３…ディジタルアナログ変換器、２５…ディスプレイ
【技術分野】
【０００１】
本発明は、例えば ディジタルアンギオ装置においてＤＳＡ画像を作成、保存、表示するＸ線診断装置等に代表される画像診断装置、画像処理装置、及び画像処理プログラムを記録した記録媒体に関する。
【背景技術】
【０００２】
造影剤を注入する前とか、観察領域に造影剤がない状態で収集した画像（マスク像）と、造影剤の注入後で観察領域に造影剤が存在する状態で収集した血管コントラストが増強された画像（コントラスト画像）との間でサブトラクション処理を行い、そのサブトラクション像を用いて行う診断は、例えば血管の狭窄位置や状態を観察するのに非常に有用であり、今日では多くの医療現場において使用されている。
【０００３】
このサブトラクション処理では、サブトラクション処理に供する画像間で位置や濃度の変化がない構造物は消去され、位置や濃度の変化のあった構造物は残存する。従って、サブトラクション処理によれば、造影前後の期間で被検体の体動がなければ、造影前後で濃度変化が生じる血管のみが抽出されることになるので、背景構造物に影響を受けることなしに血管狭窄等の観察を行うことができる。
【０００４】
ところが、このサブトラクションに供される各画像を収集する間に被検体が複雑に動くことが多く、その場合は各画像間で位置ずれが生じ偽像として映出されてしまう（なお、この様な偽像を「アーチファクト」と呼ぶ。）。従来は、例えば特許文献１に示すように、このアーチファクトを抑制するために、サブトラクションに供する前にいずれか一方の画像を求められた位置ずれ量に基づいて平行移動させて位置合せを行うことを行っていた。
【特許文献１】特開昭６４−２０８３０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、これらタイミングの間で被検体の心拍や呼吸或いは咳などによって体動が生じた場合、被検体の構造物に位置の変化が生じてしまう。特に、高齢者や重症者、意識がない被検者では息止めができないため、こういった被検者で特に腹部の検査においてはこの位置変化が甚だしい。この様に位置変化が生じた前後の画像を用いてサブトラクション処理が行われると、造影された血管が映出されるのはもちろんであるが、位置の変化が生じて完全に差分されず残存した他の構造物もアーチファクトとして映出されてしまう。実際の臨床の場面では、時々咳が止まらなくなったり、造影剤注入に熱さを感じて条件反射的に動いてしまう被検体もある。この様な場合には、横隔膜の移動によるアーチファクトが生じてしまい診断に支障をきたすことになる。
【０００６】
また、このような事情の他に心臓や大血管の拍動、体動、腸管の蠕動等の要因により、サブトラクションに供する２つの画像における被写体形状は完全同一ではないため、特許文献１に示すように平行移動による位置合せを行ってもＤＳＡ画像上には多くのアーチファクトが発生することがあり、医師による診断に悪影響を与える可能性がある。
【０００７】
本発明は、上記事情を鑑みてなされたもので、マスク像やコントラスト像の撮影中において、又はマスク像撮影からコントラスト像撮影までの間に体動が発生した場合であっても、アーチファクトがなく臨床的価値の高いＤＳＡ画像を生成することができる画像診断装置、画像処理装置、及び画像処理プログラムを記録したコンピュータ読み取り可能な記録媒体を提供することを目的としている。
【課題を解決するための手段】
【０００８】
本発明は、上記目的を達成するため、次のような手段を講じている。
【０００９】
請求項１に記載の発明は、被検体の撮影対象部位に関する複数の画像を、当該被検体への造影剤の注入の前後にわたって連続的に収集する画像収集手段と、前記複数の画像のうちの一画像である第１の画像を少なくとも部分的に変形させる変形処理を実行することで、前記複数の画像のうちの前記第１の画像とは異なる第２の画像に表示された各種部位との位置合わせが施された第３の画像を生成する変形処理手段と、前記各第２の画像と前記第３の画像の間でサブトラクションすることで、第１のサブトラクション画像を生成する画像生成手段と、を具備することを特徴とする画像診断装置である。
【００１０】
請求項１１に記載の発明は、被検体への造影剤の注入の前後にわたって連続的に収集された、当該被検体の撮影対象部位に関する複数の画像を記憶する記憶手段と、前記複数の画像のうちの一画像である第１の画像を少なくとも部分的に変形させる変形処理を実行することで、前記複数の画像のうちの前記第１の画像とは異なる第２の画像に表示された各種部位との位置合わせが施された第３の画像を生成する変形処理手段と、前記各第２の画像と前記第３の画像の間でサブトラクションすることで、第１のサブトラクション画像を生成する画像生成手段と、を具備することを特徴とする画像処理装置である。
【００１１】
請求項１６に記載の発明は、コンピュータに、被検体への造影剤の注入の前後にわたって連続的に収集された、当該被検体の撮影対象部位に関する複数の画像のうちの一画像である第１の画像を少なくとも部分的に変形させる変形処理を実行することで、前記複数の画像のうちの前記第１の画像とは異なる第２の画像に表示された各種部位との位置合わせが施された第３の画像を生成する第１の機能と、前記各第２の画像と前記第３の画像の間でサブトラクションすることで、第１のサブトラクション画像を生成する第２の機能と、を実現させるための画像処理プログラムを記録したコンピュータ読み取り可能な記録媒体である。
【発明の効果】
【００１２】
以上本発明によれば、マスク像やコントラスト像の撮影中において、又はマスク像撮影からコントラスト像撮影までの間に体動が発生した場合であっても、アーチファクトがなく臨床的価値の高いＤＳＡ画像を生成することができる画像診断装置、画像処理装置、及び画像処理プログラムを記録した記録媒体を実現することができる。
【発明を実施するための最良の形態】
【００１３】
以下、本発明の実施形態を図面に従って説明する。以下の説明において、略同一の機能及び構成を有する構成要素については、同一符号を付し、重複説明は必要な場合にのみ行う。
【００１４】
なお、本発明は造影画像を収集可能な画像診断装置を含んでいる。この種の画像診断装置としては、Ｘ線診断装置、超音波診断装置、磁気共鳴イメージング装置、さらにＰＥＴ，ＳＰＥＣＴが含まれ、そのいずれでも本発明に適用可能である。以下では、Ｘ線診断装置を例に説明する。
【００１５】
図１は本実施形態に係るＸ線診断装置の構成を示している。本Ｘ線診断装置は、Ｘ線管１とＸ線検出器３とを有する画像収集部２、インジェクタ５、画像処理装置７を具備している。
【００１６】
画像収集部２は、寝台に載置された被検体Ｐを挟んで対向配置されたＸ線管１とＸ線検出器３とを有する。
【００１７】
Ｘ線管１は、Ｘ線を発生する真空管であり、高電圧発生装置からの高電圧により電子を加速させ、ターゲットに衝突させることでＸ線を発生する。発生したＸ線は、被検体に対して曝射される。
【００１８】
Ｘ線検出器３は、シンチレータとフォトダイオードアレイとを有するフラットパネルデテクタ（ＦＰＤ）が採用される。しかし、Ｘ線検出器３は、イメージインテンシファイアとＴＶカメラとの組み合わせが採用されても良い。なお、ＦＰＤは、被検体を透過したＸ線を光電膜に当てることで電子正孔を生成し、これを半導体スイッチにおいて蓄積し、電気信号として読み出すことでＸ線信号を検出するものである。また、Ｉ．Ｉ．は、被検者を透過したＸ線を入力蛍光面で光学情報に変換し、さらに当該光学情報に基づいて生成された光電子を出力蛍光面に衝突させて高輝度光学画像を生成する。このＸ線検出器３から出力されるアナログの映像信号は、アナログディジタル変換器１３から画像処理装置７に供給される。
【００１９】
インジェクタ５は、被検体Ｐに造影剤を自動注入するための装置である。造影剤が被検体Ｐに注入された時、トリガ信号又はその注入時刻を表す信号が、インジェクタインタフェース１１を介して画像処理ユニット７に供給される。
【００２０】
画像処理装置７は、インジェクタインタフェース１１、アナログディジタル変換器１３とともに、ＣＰＵ９、画像データや各種制御データ等を記憶する記憶装置１５、コンソール（操作卓）１７、マスク選択部１９、サブトラクション処理部２１、ワーピング処理部２２、ディジタルアナログ変換器２３、ディスプレイ２５から構成される。
【００２１】
ＣＰＵ９は、画像データの収集に関する制御、及び収集した画像データの画像処理、画像再生処理等に関する制御を行う中央処理装置である。
【００２２】
インジェクタインタフェース１１は、インジェクタ５とＣＰＵ９との間のデータ通信を可能にするための装置である。
【００２３】
アナログディジタル変換器１３は、Ｘ線検出器３から入力したアナログ信号を、ディジタル信号に変換し、フレーム毎の画像データとして記憶装置１５、マスク選択部１６等に送り出す。
【００２４】
記憶装置１５は、取得された画像データ、画像取得及び画像処理のための各種制御データ等を記憶する。また、記憶装置１５は、ＣＰＵ９による制御のもと、造影剤注入の前後にわたって収集された複数の画像を、後述する遅延時間に基づいて、マスク像とコントラスト像とに区分して記憶する。さらに、記憶装置１５は、後述するレビュー機能を実現するために、マスク像、コントラスト像、ワーピングＤＳＡ処理によって生成されたワーピング画像等を記憶する。
【００２５】
コンソール（操作卓）１７は、キーボードや各種スイッチ、マウス等を備えた入力装置であり、撮影指示や画像選択等を入力する際に使用される。
【００２６】
マスク選択部１９は、造影剤注入後に収集した各コントラスト像について、後述する遅延時間内に収集した複数のマスク像からＤＳＡ処理に最適な１枚を最適マスク像として選択する。その選択手法は種々のものを採用することができる。本実施形態では、後述するように、撮影プログラムに応じていずれの選択手法を採用するかを決定する。
【００２７】
サブトラクション処理部２１は、各コントラスト像と、コントラスト像毎に選択されたマスク像とを用いてサブトラクション処理を実行し、ＤＳＡ画像を生成する。また、サブトラクション処理部２１は、ワーピング処理部２２において最適マスク像に基づいて生成されたワーピング画像（後述）と、各コントラスト像とを用いてサブトラクション処理を実行し、ワーピングＤＳＡ画像を生成する。
【００２８】
ワーピング処理部２２は、サブトラクション処理に用いられるコントラスト像と最適マスク像との間で形状が異なる部位を、最適マスク画像において画素単位で部分的に変形（ワーピング）させることで位置合わせを行うワーピング処理を実行する。このワーピング処理、及び当該ワーピング処理を前提とするワーピングＤＳＡ処理については、後で詳しく説明する。
【００２９】
ディジタルアナログ変換器２３は、サブトラクション処理部２１等から入力した画像データのディジタル信号列を、アナログ信号列に変換する。
【００３０】
ディスプレイ２５は、ディジタルアナログ変換器２３から受け取った信号により、Ｘ線診断画像等を表示するＣＲＴ、プラズマディスプレイ、液晶ディスプレイ等の表示装置である。特に、本ディスプレイは医用機器に使用されるものであるため、走査線１０００本以上の高精細ディスプレイであることが好ましい。
【００３１】
（ワーピングＤＳＡ処理）
本Ｘ線診断装置は、被検体の体動に起因して発生するコントラスト像に対するマスク像の位置ずれを、ワーピングを利用して修正した後にＤＳＡ処理を実行する処理（ワーピングＤＳＡ処理）を有している。このワーピングＤＳＡ処理により、マスク像やコントラスト像の撮影中において、又はマスク像撮影からコントラスト像撮影までの間に体動が発生した場合であっても、アーチファクトがなく臨床的価値の高いＤＳＡ画像を生成することができる。以下、この処理について詳しく説明する。
【００３２】
なお、被検体の体動は、周期的なものと非周期的なものとのに分類することができる。周期的な体動とは、心臓や大血管の拍動、呼吸による横隔膜や肺の部位の動き等である。また、非周期的な体動とは、腸管の蠕動、撮影中の咳、造影剤流入を体感した場合における条件反射等である。本ワーピングＤＳＡ処理は、いずれの体動をも対象とすることができる。
【００３３】
図２は、本ワーピングＤＳＡ処理の流れを示したフローチャートである。同図に示すように、本ワーピングＤＳＡ処理の各ステップでは、人体構造の解析、画像濃度の補正、ローカルマッチング、画像のワーピング、差分画像の生成の各処理が実行される。
【００３４】
まず、人体構造の解析処理（ステップＳａ）では、撮影したマスク像に基づいて人体構造物を解釈し、人体解析データを作成する。ここで、人体解析データとは、例えば胸部の血管撮影の場合であれば、横隔膜、脊椎、心臓など、マスク像上に写っているランドマークとなる人体構造物の位置、形状に関する情報である。
【００３５】
図３は、腹部アンギオ画像においてランドマークとする人体構造物（横隔膜、脊椎、心臓、肋骨）の一例を示した図である。また、図４は、胸部アンギオ画像においてランドマークとする人体構造物（胸郭、横隔膜、脊椎）の一例を示した図である。
【００３６】
この人体構造の解析処理においては、複数枚のマスク像を用いて平均化等された情報に基づいて、人体解析データを作成することも可能である。この場合、同一の撮影シーケンスにおいて撮影されたマスク像に限定されない。例えば、同一患者の同一検査のアンギオ画像であれば、それほど大きな変動はしないことが予想されるから、過去に撮影された同一患者の同一検査のマスク像や人体解析データも利用することができる。特に、骨、肺、腸管ガス、体外の空気等の高コントラストを形成する部位等、また、差分に用いる２画像の位相差を類推する目的で血管内に挿入された、コイル、ステント、ＩＶＣフィルター、カテーテルのマーカーや過去の手術で体内に挿入されたクリップ、義歯などは、人体解析データとして利用可能である。また、後述するワーピング処理によって体動補正をより効率よく行うために、前もって被検者の体表に付着させ金属等のマーカーを人体解析データとして利用することも可能である。
【００３７】
なお、一般に、撮影対象部位（例えば、頭部、頚部、胸部、腹部、四肢）によって、ランドマークとなる人体構造物は異なる。従って、人体解析データを作成する対象部位を予め指定すれば、人体構造物の解釈がより容易に行うことができる。
【００３８】
また、同一の撮影対象部位であっても、正面、側面、斜位といった撮影方向によって、ランドマークとなる人体構造物は異なる。本解析処理によって得られた人体解析データは、同一の撮影対象部位であれば、いかなる撮影方向にも適応させることが可能である。
【００３９】
この様に生成される人体解析データは、後段のステップにおいて、人体構造の知識（解剖学的情報）に沿った処理ができるように利用される。
【００４０】
次に、画像濃度の補正（ステップＳｂ）では、コントラスト像、マスク像（今の場合、最適マスク像）に対して、画像濃度の補正処理を適用して、補正コントラスト像、補正マスク像を作成する。この補正処理により、腸管ガスの影響除去、ノイズ除去、アンギオ装置の撮影時のフィルター除去など、造影剤注入以外の原因に起因する、ＤＳＡ画像画質低下の要因を低減させることができる。
【００４１】
次に、ローカルマッチング（ステップＳｃ）では、ステップ２で生成された補正コントラスト像と補正マスク像について、両画像が局所的に一致するような関心領域を多数配置し、これらに基づいて移動ベクトルを導出する。具体的なローカルマッチングの手法は、次の様である。
【００４２】
図５（ａ）、図５（ｂ）は、ローカルマッチングの手法の一例を説明するための図である。図５（ａ）に示すように、まず、ある一定サイズのテンプレートＲＯＩ５０を補正コントラスト像上に配置し、当該テンプレートＲＯＩ５０に対応する（複数個の場合には、それぞれに対応する）探索領域ＲＯＩ５４を補正マスク像上に配置する。そして、探索領域ＲＯＩ５４とテンプレートＲＯＩ５０との間で、テンプレートマッチングと呼ばれる技法を用いて、局所的な移動量（移動ベクトル）を導出する。なお、一般的には、テンプレートＲＯＩが複数設定されるため、移動ベクトルは、個々のテンプレートＲＯＩに対応して導出される。
【００４３】
なお、移動ベクトルを導出するためのテンプレートＲＯＩ５０及び探索領域ＲＯＩ５４は、ステップＳａにおいて取得した人体解析データに基づいて決定された各画像の所定エリア内に、例えば一定の間隔をおいて配置される。
【００４４】
図６は、補正コントラスト像又は補正マスク像としての腹部アンギオ画像上において、複数のテンプレートＲＯＩ５０又は複数の探索領域ＲＯＩ５４を配置するエリアの一例を示した図である。同図に示す腹部アンギオ画像の例では、テンプレートＲＯＩ５０又は探索領域ＲＯＩ５３を配置しないのは、脊椎のような不動部分や、横隔膜より上の部分などである。
【００４５】
次に、画像のワーピング（ステップＳｄ）では、ステップ２で導出された補正マスク像について、ステップＳｃで導出された移動ベクトルを基に、例えば画素単位の移動ベクトルを導出し、これに基づいて補正マスク像上に表示された所定部位を変形（ワープ）させるワーピング処理を実行する（以下、このワーピング処理によって得られる画像を「ワーピング画像」と呼ぶ。）。
【００４６】
この様なワーピング処理を行うのは、次の理由による。すなわち、補正マスク像と補正コントラスト像に写っている人体構造物は、ほぼ同じ形状で写っているものと考えてよい。しかしながら、呼吸、体動、心臓や大血管の拍動、腸管の蠕動により変形が発生する部位については、ＤＳＡ対象となる補正マスク像と補正コントラスト像との間で形状が異なる場合がある。この様な変形が発生する部位については、補正マスク像を画素単位で部分的に変形させ（画像ワーピング）、補正コントラスト像に合せこむことによって、その変形の効果を相殺する。これにより、ＤＳＡ画像におけるアーチファクトの原因となる補正マスク像と補正コントラスト像との位置ずれをなくすことができる。
【００４７】
なお、ローカルマッチング、及びワーピング処理のアルゴリズムについては、例えば特開平７−３７０７４号公報に記載されている経時的差分画像作成と同じアルゴリズムを適用することができる。
【００４８】
次に、差分画像の生成（ステップＳｅ）では、ステップＳｂにおいて生成された補正コントラスト像から、ステップＳｄにおいて生成されたワーピング画像を差し引くことでＤＳＡ画像を生成し、表示部に出力する。
【００４９】
（撮影シーケンス）
次に、上記ワーピングＤＳＡ処理を含む一連の撮影シーケンスについて説明する。以下説明を簡単にするため、被検体に対して撮影系を移動させない静的（static）なＤＳＡ撮影を対象とする。しかしながら、本発明は、ステッピングＤＳＡ撮影、回転ＤＳＡ撮影といった動的な撮影手法においても適用可能である。
【００５０】
図７は、当該撮影シーケンスにおいて実行される処理の流れを示したフローチャートである。同図に示すように、まず、患者ＩＤ、撮影プログラム等の入力・選択等が実行される（ステップＳ１）。ここで、撮影プログラムとは、撮影領域（部位）別に固有の遅延時間（マスク像の撮影開始からコントラスト像の撮影開始までの期間）、撮影条件、後処理等を実行するための制御プログラムである。本実施形態においては、「腹部専用プログラム」、「心臓専用プログラム」、「その他の領域（例えば、頭部、四肢等）用プログラム」等の複数種の撮影プログラムが選択可能であるとする。
【００５１】
次に、ステップＳ１において選択された撮影プログラムに合わせて、管電流／管電圧等、遅延時間（マスク像の撮影開始からコントラスト像の撮影開始までの期間）の設定を実行する（ステップＳ２）。なお、遅延時間は、確実に一呼吸分のマスク像を収集するようにするため、一呼吸期間よりも少し長めに設定するのが一般的である。
【００５２】
図８は、ステップＳ２において利用される入力画面の一例を示した図である。同図において、例えばステップＳ１において腹部専用プログラムが選択された場合には、１呼吸期間に多少のマージンを考慮して、最初の数秒間は造影剤を注入しないように遅延時間を設定する。また、ワーピングＤＳＡ処理についてより好適なマスク像を選択可能とするために、遅延時間内の画像収集フレームレートは、遅延時間経過後の画像収集フレームレートより高くなるように、例えば遅延時間内の画像収集フレームレートを１５／３０ｆｐｓ、遅延時間経過後の画像収集フレームレートを７．５／３ｆｐｓとして設定する。
【００５３】
また、例えばステップＳ１において心臓専用プログラムが選択された場合には、呼吸動より速い心拍動を考慮して、最初の数秒間は造影剤を注入しないように遅延時間を設定する。また、同様に、遅延時間内の画像収集フレームレートが遅延時間経過後の画像収集フレームレートより高くなるように、例えば遅延時間内の画像収集フレームレートを３０／６０ｆｐｓ、遅延時間経過後の画像収集フレームレートを１５／３０ｆｐｓとして設定する。
【００５４】
なお、ステップＳ１においてその他の領域用プログラムが選択された場合には、本実施形態においては、遅延時間、フレームレート等の設定には特に条件を設けない。しかしながら、これに限定する趣旨ではなく、その他の領域用プログラムが選択された場合であっても、必要に応じて領域毎の固有の遅延時間、フレームレート等を設定可能な構成にしてもよい。また、当該ステップＳ２における遅延時間、フレームレート等の設定は、ステップＳ１において選択された撮影プログラムに応じて、自動的に設定される構成であってもよい。
【００５５】
次に、専用ＤＳＡ処理プログラム実行の実行／非実行を選択する（ステップＳ３）。ここで、専用ＤＳＡ処理プログラムとは、複数枚のマスク像の中から自動的にＤＳＡ処理に最適なマスク像（最適マスク像）を自動的に選択し、これを用いてワーピングＤＳＡ処理を実行するものである。
【００５６】
次に、複数枚のマスク像撮影及び複数枚のコントラスト像撮影が実行される（ステップＳ４）。このとき、コントラスト像の最初の撮影は、マスク像の最初の撮影から遅延時間経過後に開始される。
【００５７】
次に、ステップＳ３において専用ＤＳＡ処理プログラムが選択された場合には、当該プログラムに従って最適マスク像の自動選択処理が実行される（ステップＳ５ａ）。この自動選択処理としては、例えば特開２００４−１１２４６９号公報に開示されている手法を採用することができる。この手法では、各マスク像と各コントラスト像との「位置ずれの程度の類似性を示す一致度」を計算することで、位置ずれの程度が最も近いマスク像（最適マスク像）をコントラスト像毎に選択し、これを用いてＤＳＡ処理を実行するものである。これにより、例えば体動があった領域について、コントラスト像に最も近い形状を有するマスク像を選択することができる。
【００５８】
なお、最適マスク像は必ずしもコントラスト像毎に選択する必要はなく、複数のコントラスト像に共通する最適マスク像（共通最適マスク像）を選択する構成であってもよい。本実施形態では、説明を簡単にするため、共通最適マスク像を選択するものであるとする。
【００５９】
また、自動選択処理は、上記手法に限定されない。例えば、心臓専用プログラムを実行する場合には、心電波形の位相が略一致するフレーム同士を選択しサブトラクションする構成であることが好ましい。また、その他の自動選択処理として、心臓や大血管の拍動では指先でモニターする脈波を、呼吸では被検者に装着した呼吸同期装置を、体軸の移動では被検者の皮膚に塗布したマーカーをそれぞれ用いることも可能である。
【００６０】
次に、ステップＳ５ａにおいて自動選択された最適マスク像と各コントラスト像とを用いて、ワーピングＤＳＡ処理を実行する（ステップＳ６ａ）。ワーピングＤＳＡ処理の内容は、図２のステップＳａ〜Ｓｅに示した通りである、当該ワーピングＤＳＡ処理によって得られたＤＳＡ画像は、ディスプレイ２５に表示される。
【００６１】
なお、当該ワーピングＤＳＡ処理中においては、ＤＳＡ画像又はコントラスト像と（ワーピング処理前の）最適マスク像とＤＳＡ画像を、中間処理画像として表示することが好ましい。ワーピングＤＳＡ画像が生成されると、ディスプレイ２５上は、自動的に中間処理画像からワーピングＤＳＡ画像に切り換えて表示する。この様な構成によれば、操作者はワーピングＤＳＡ処理中であっても診断を行うことができ、また、中間処理画像とワーピングＤＳＡ画像との双方を臨床情報として有効に利用することができる。
【００６２】
一方、ステップＳ３において専用ＤＳＡ処理プログラムが選択されなかった場合には、ステップ５ｂ〜６ｅまでの流れに従って処理される。すなわち、ＤＳＡ処理用として初期設定されたマスク像（例えば、最初に撮影されたマスク像）が最適マスク像として選択される（ステップＳ５ｂ）。
【００６３】
次に、例えば図９に示すような操作画面から「Ｗａｒｐｉｎｇ」のボタンがクリックされると、初期設定による最適マスク像と現在表示中のコントラスト像（カレント像）とを用いて、ワーピングＤＳＡ処理が実行される（ステップＳ６ｂ）。当該ワーピングＤＳＡ処理によって得られたＤＳＡ画像は、ディスプレイ２５に表示される。
【００６４】
次に、操作者は、表示されたＤＳＡ画像を観察し、リマスク処理を実行する必要があるか否かを判定し（ステップＳ６ｃ）、実行する場合には、最適マスク像の変更が実行される（ステップ６ｅ）。ここで、「リマスク処理」とは、最適マスク像の変更を行うための処理であり、例えば次の様なものがある。
【００６５】
図１０は、リマスク処理の一例を説明するための図である。同図に示すように、初期設定により「収集画像１」が全コントラスト像に対して共通する最適マスク像として設定されているとする。この状態で、図９に示した入力画面において「Ｒｅｍａｓｋ」のボタンをクリックすると、他の任意の収集画像を最適マスク像として再設定することができる。実際の撮影においては、造影剤注入の際に患者が「熱い」と感じる場合があり、このときに体動が発生する場合が多い。従って、造影剤注入が終了し且つ造影剤が撮影対象に流入する手前である画像（例えば、「収集画像８」）を新たな最適マスク像として設定することで、好適なＤＳＡ処理を実現することができる。また、ステップＳ５ａではあらかじめマスク画像を一致度の高い順に並べ替え、リマスクの候補の選択に用いてもよい。なお、図１０に示した例は、例えば静的なＤＳＡ撮影、ステッピングＤＳＡ撮影等において、特に有益である。
【００６６】
また、図１１は、回転ＤＳＡ撮影におけるリマスク処理の例を説明するための図である。回転ＤＳＡ撮影では、初期設定として、同じ角度で収集したマスク像ｉ（ｉ＝１〜ｎ）とコントラスト像ｉとの間でサブトラクション処理が実行される。この状態で、図９に示した入力画面において「Ｒｅｍａｓｋ」のボタンをクリックすることで、所定のコントラスト像に対するマスク像だけを変更することができる。なお、図１１では、「コントラスト像１」に対する最適マスク像を「マスク像１」から「マスク像３」に変更する例を示してある。
【００６７】
また、図１２は、回転ＤＳＡ撮影におけるリマスク処理の他の例を説明するための図である。同図に示すように、リマスク処理によって新たに選択したマスク像を基準に、全てのマスク像をコントラスト像に対して一定フレーム数ずらしてサブトラクション処理を実行することも可能である。
【００６８】
ステップＳ６ｅにおいて最適マスク像が変更されると、操作画面からの「Ｗａｒｐｉｎｇ」ボタンの操作に応答して、リマスク処理後の最適マスク像とカレント像とを用いて、ワーピングＤＳＡ処理が実行され、得られたＤＳＡ画像がディスプレイ２５に表示される。操作者は、表示されたＤＳＡ画像を観察し、再度リマスク処理を実行する必要があるか否かを判定し、必要ならば最終的な最適マスク像を決定するまでリマスク処理を繰り返す。最終的な最適マスク像は、例えば図９に示す「ｃｏｎｆｉｒｍ」ボタンをクリックすることで決定され、以降のコントラスト像に対するワーピングＤＳＡ処理は、当該最適マスク像を用いて実行される。一方、「ｃｏｎｆｉｒｍ」ボタンをクリックしなければ、当該最適マスク像は、カレント像に対してのみ適用されることになる。
【００６９】
次に、以降のコントラスト像に対して、決定された最適マスク像を用いてワーピングＤＳＡ処理を実行する（ステップＳ６ｄ）。このとき、例えば図９に示す「ａｐｐｌｙ」ボタンをクリックすれば、マルチフレーム画像に対して一度にワーピングＤＳＡ処理を実行することができる。
【００７０】
なお、ステップＳ３において専用ＤＳＡ処理プログラムが選択された場合であっても、必要に応じてワーピングを適用しないＤＳＡ処理を実行することも可能である。すなわち、専用ＤＳＡ処理プログラムが選択された場合には、図９に示す「Ｗａｒｐｉｎｇ」ボタンは初期設定でＯＮ状態となっている。また、専用ＤＳＡ処理プログラム選択後に撮影された画像は、初期設定でＣｈｅｃｋ−ＢｏｘがＯＮ状態となっており、また、「Ｗａｒｐｉｎｇ」ボタンがＯＮ状態となっている時にＣｈｅｃｋ−Ｂｏｘが選択可能となる。このＣｈｅｃｋ−ＢｏｘをＯＦＦ状態とした画像については、ワーピング処理前の補正マスク像が最適マスク像として選択することができる。
【００７１】
以上述べた処理によって得られるＤＳＡ画像等は、自動的に記憶装置１５に記憶される（ステップＳ７）。
【００７２】
（レビュー機能）
次に、本Ｘ線診断装置が有するレビュー機能について説明する。本機能は、ワーピングＤＳＡ処理によって得られた画像（例えば、図７のステップＳ６ａ、Ｓ６ｂ、Ｓ６ｄにおいて得られた画像）を事後的に再生するものである。
【００７３】
まず、図７のステップＳ７において各コントラスト像と各マスク像、及び各コントラスト像に対応するワーピング画像を自動的に記憶装置１５に保存する。この様に各コントラスト像に対応するワーピング画像を保存することで、ワーピング画像生成を省略したワーピングＤＳＡ処理を事後的に実行でき、ワーピングＤＳＡ画像を迅速に再生することができる。なお、メモリ資源をできるだけ有効に活用する観点から、保存するワーピング画像については非可逆圧縮を施す構成であってもよい。
【００７４】
また、ワーピング画像そのものを保存する構成ではなく、各コントラスト像に対応する最適マスク像のフレーム番号を保存し、ワーピング画像生成を含むワーピングＤＳＡ処理を事後的に実行することで、ワーピングＤＳＡ画像を再生することも可能である。この場合、最適マスク像からワーピング画像を生成する過程を保存しておくことで、ワーピング画像そのものを保存するよりもデータ量を小さくし、且つ最適マスク像のフレーム番号を保存するよりも事後再生において迅速化を図る構成としてもよい。また、この様にワーピング画像生成過程を保存する場合には、そのログの解析を元に、作成されたＤＳＡ画像の画質の推測をし、当該ＤＳＡ画像と供に表示することも可能である。
【００７５】
以上述べた構成によれば、以下の効果を得ることができる。
【００７６】
本Ｘ線診断装置においては、最適マスク像からワーピング画像を生成し、これとコントラスト像とを用いてＤＳＡ処理を実行する。従って、最適マスク像とコントラスト像との単純差分によりＤＳＡ処理を実行する従来技術に比して、より位置精度の高いサブトラクション処理を実行することができる。その結果、マスク像やコントラスト像の撮影中において体動が発生した場合であっても、従来よりもアーチファクトがなく臨床的価値の高いＤＳＡ画像を提供することができる。特に、ワーピング処理は、部分的な位置ずれを補正するのに好適であり、マスク画像とコントラスト像とが完全に同位相になるような撮影は困難である場合には、特に実益がある。
【００７７】
また、本Ｘ線診断装置においては、撮影プログラムに応じて遅延時間内の画像収集フレームレートを設定することができる。従って、コントラスト像より多くのマスク像を収集することができ、ＤＳＡ処理においてより位置精度の高いマスク像を選択することができる。
【００７８】
また、本Ｘ線診断装置においては、複数のマスク像の中から最適マスク像を、例えば位置ずれの程度の類似性を示す一致度に基づいて、自動的に選択することができる。従って、コントラスト像と位置対応性の高い最適マスク像を迅速且つ簡便に選択することができる。
【００７９】
また、本Ｘ線診断装置においては、必要に応じて最適マスク像を他の任意のマスク像に変更することができる。従って、自動選択された最適マスク像、又は初期設定されている最適マスク像によるＤＳＡ画像以外にも、他のマスク像を最適マスク像とすることで、種々の所望のＤＳＡ画像を取得することができ、診断の自由度を広げることができる。
【００８０】
また、本Ｘ線診断装置においては、レビュー機能により、ワーピングＤＳＡ処理によって得られた画像を事後的に再生することができる。これにより、撮影以後の診断や、画像再生による撮影・診断の事後的な検討が可能となり、医療行為の質の向上に寄与することができる。たとえば、血管奇形、動静脈瘻、肝細胞癌などでは動脈塞栓術の前後で、動脈の閉塞や狭窄病変では人工血管によるバイパス術、血管形成術、ステント留置術前後で、それぞれの造影剤注入後の画像同士の差分画像を作成し、治療の効果を客観的に表示することも可能となる。
【００８１】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。具体的には、次のような変形例を挙げることができる。
【００８２】
（１）ワーピングにより十分補正できなかった腸管の蠕動などのアーチファクトは、例えばさらにダイナミックレンジ圧縮の技法を用い、低吸収値部分のコントラストを減弱させ、アーチファクトを軽減させる構成としてもよい。
【００８３】
（２）胸部の血管撮影、特に冠動脈撮影では、心拍動によるアーチファクトが問題となる。係る撮影においては、例えば初期設定として、心電図同期によるＲ波からのトリガー遅延を最適マスク像の選択に使用する構成であってもよい。
【００８４】
（３）血管撮影装置のＸ線管球と検出器を移動させながら撮影する３次元撮影の場合、もしくは被検者を載せたテーブルを移動させながら撮影を行う場合では、被検者とＸ線管球および検出器の空間的位置関係を、最適マスク像の選択に利用することもできる。
【００８５】
（４）一般の血管撮影装置に併設されている、ＣＴ画像にも本ワーピングＤＳＡ処理は適用可能である。この場合、ＩＶＲＣＴ（Interventional Radiology CT）で撮影された造影前と造影後のボリュームデータを２次元に投影したＲａｙＳｕｍ（もしくはＲａｙ−Ｓｕｍ）表示画像やＭＩＰ（Maximum Intensity Projection）画像等を用いて本技術を応用することができる。なお、これらＲａｙＳｕｍ表示画像等についての詳細は、特開２００２−２１９１２３公報に記載されている。
【００８６】
（５）上記実施形態においては、Ｘ線診断装置において、ワーピングＤＳＡ処理等を実施する例を示した。これに対し、Ｘ線診断装置と高速なインターフェースで連結され、例えばＤＩＣＯＭ（Digital Imaging and Communications in Medicine）による通信を行う専用のワークステーションによって、ワーピングＤＳＡ処理等を実施する構成であってもよい。ＰＡＣＳ（Picture archiving and communications system）のネットワーク上にワークステーションを配置して、血管撮影のデータを取得し任意に本手法を適応することも可能である。
【００８７】
なお、この場合、ＤＩＣＯＭのプライベートタグ中に遅延時間、画像収集レート等の情報を記録することで、ワークステーション上においても、ワーピングＤＳＡ処理を含む一連の撮影シーケンス（すなわち、図７に示したシーケンス）を擬似的に再現することができる。
【００８８】
（６）本実施形態に係るワーピングＤＳＡ処理等は、当該処理を実行するプログラムをワークステーション等のコンピュータにインストールし、これらをメモリ上で展開することによっても実現することができる。このとき、コンピュータに当該手法を実行させることのできるプログラムは、磁気ディスク（フロッピー（登録商標）ディスク、ハードディスクなど）、光ディスク（ＣＤ−ＲＯＭ、ＤＶＤなど）、半導体メモリなどの記録媒体に格納して頒布することも可能である。
【００８９】
また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【産業上の利用可能性】
【００９０】
以上本発明によれば、マスク像やコントラスト像の撮影中において体動が発生した場合であっても、アーチファクトがなく臨床的価値の高いＤＳＡ画像を生成することができるＸ線診断装置、画像処理装置、及び画像処理プログラムを格納する記憶媒体を実現できる。
【図面の簡単な説明】
【００９１】
【図１】図１は本実施形態に係るＸ線診断装置の構成を示している。
【図２】図２は、本ワーピングＤＳＡ処理の流れを示したフローチャートである。
【図３】図３は、腹部アンギオ画像においてランドマークとする横隔膜、脊椎、心臓、肋骨の例を示した図である。
【図４】図４は、腹部アンギオ画像においてランドマークとする横隔膜、脊椎、心臓、肋骨を示した図である。
【図５】図５は、ローカルマッチングの手法の一例を説明するための図である。
【図６】図６は、腹部アンギオ画像上において、複数のテンプレートＲＯＩ５０又は複数の探索領域ＲＯＩ５３を配置するエリアの一例を示した図である。
【図７】図７は、ワーピングＤＳＡ処理を含む一連の撮影シーケンスにおいて実行される処理の流れを示したフローチャートである。
【図８】図８は、遅延時間等の入力画面の一例を示した図である。
【図９】図９は、ワーピング処理、リマスク処理を指示するための入力画面の一例を示した図である。
【図１０】図１０は、リマスク処理の一例を説明するための図である。
【図１１】図１１は、回転ＤＳＡ撮影においてリマスク処理の例を説明するための図である。
【図１２】図１２は、回転ＤＳＡ撮影におけるリマスク処理の他の例を説明するための図である。
【符号の説明】
【００９２】
１…Ｘ線管、２…画像収集部、３…Ｘ線検出器、５…インジェクタ、７…画像処理装置、９…ＣＰＵ、１１…インジェクタインタフェース、１３…アナログディジタル変換器、１５…記憶装置、１７…コンソール、１９…マスク選択部、２１…サブトラクション処理部、２２…ワーピング処理部、２３…ディジタルアナログ変換器、２５…ディスプレイ
【特許請求の範囲】
【請求項１】
被検体の撮影対象部位に関する複数の画像を、当該被検体への造影剤の注入の前後にわたって連続的に収集する画像収集手段と、
前記複数の画像のうちの一画像である第１の画像を少なくとも部分的に変形させる変形処理を実行することで、前記複数の画像のうちの前記第１の画像とは異なる第２の画像に表示された各種部位との位置合わせが施された第３の画像を生成する変形処理手段と、
前記各第２の画像と前記第３の画像の間でサブトラクションすることで、第１のサブトラクション画像を生成する画像生成手段と、
を具備することを特徴とする画像診断装置。
【請求項２】
前記第１の画像及び前記第２の画像が複数ある場合には、当該複数の第１の画像から前記各第２の画像との前記サブトラクション処理に利用する最適画像を選択する選択手段をさらに具備することを特徴とする請求項１記載の画像診断装置。
【請求項３】
前記変形処理手段は、前記第１の画像又は前記第２の画像上の標識となる被検体構造物の位置、形状その他の幾何学的情報に基づいて前記変形処理を実行し、前記第３の画像を生成することを特徴とする請求項１又は２記載の画像診断装置。
【請求項４】
前記変形処理手段は、
前記最適画像上に少なくとも一つの第１の関心領域を設定し、
前記各第２の画像上に前記第１の関心領域と対応させた少なくとも一つの第２の領域を設定し、
前記少なくとも一つの第１の関心領域と、前記少なくとも一つの第２の関心領域との間の位置ずれを検出することにより、前記各第２の画像に対する前記最適画像の位置ずれ量を取得し、当該前記位置ずれ量に基づいて、前記変形処理を実行すること、
を特徴とする請求項１乃至３のうちいずれか一項記載の画像診断装置。
【請求項５】
前記被検体への造影剤注入開始からの所定時間に基づいて、前記複数の画像を前記複数の第１画像と前記複数の第２画像とに区分する画像区分手段をさらに具備することを特徴とする請求項１乃至４のうちいずれか一項記載の画像診断装置。
【請求項６】
前記撮影対象部位に応じて、前記所定時間、及び当該所定時間内の画像収集レートの少なくとも一方を所望の値に設定するための設定手段をさらに具備し、
前記画像区分手段は、前記設定手段によって前記所定時間が設定された場合には、
設定後の前記所定時間に基づいて、前記複数の第１画像と前記複数の第２画像との区分を実行し、
前記画像収集手段は、前記設定手段によって前記画像収集レートが設定された場合には、設定後の前記画像収集レートに基づいて、前記画像収集を実行すること、
を特徴とする請求項５記載の画像診断装置。
【請求項７】
前記選択手段は、前記複数の第１の画像から前記最適画像を選択するための複数の選択手法を有し、指定される撮影プログラムに応じて、前記複数の選択手法から選択手法を決定することを特徴とする請求項１乃至６のうちいずれか一項記載の画像診断装置。
【請求項８】
前記画像生成手段は、前記各第２の画像から前記最適画像をサブトラクションすることで、第２のサブトラクション画像を生成し、
前記変形処理手段により前記第３の画像が生成されるまでの間、前記画像収集手段によって収集された画像又は前記第２のサブトラクション画像を表示する表示手段をさらに具備すること、
を特徴とする請求項１乃至７のうちいずか一項記載の画像診断装置。
【請求項９】
前記複数の第２の画像と、前記第２の画像毎の前記第３の画像とを記憶する記憶手段をさらに具備し、
前記画像生成手段は、前記記憶手段に記憶された前記各第２の画像と前記第３の画像とを用いて、前記第１のサブトラクション画像を生成すること、
を特徴とする請求項１乃至８のうちいずれか一項記載の画像診断装置。
【請求項１０】
前記複数の第１の画像と、前記複数の第２の画像と、前記第２の画像毎の前記最適画像としての前記第１の画像のフレーム番号と、を記憶する記憶手段をさらに具備し、
前記変形処理手段は、前記記憶手段に記憶された前記各第２の画像と、前記複数の第２の画像と、前記第２の画像毎の前記最適画像としての前記第１の画像のフレーム番号とを用いて前記変形処理を実行することで、前記第３の画像を事後的に生成し、
前記画像生成手段は、前記記憶手段に記憶された前記各第２の画像と事後的に生成された前記第３の画像とを用いて、前記第１のサブトラクション画像を生成すること、
を特徴とする請求項１乃至９のうちいずれか一項記載の画像診断装置。
【請求項１１】
被検体への造影剤の注入の前後にわたって連続的に収集された、当該被検体の撮影対象部位に関する複数の画像を記憶する記憶手段と、
前記複数の画像のうちの一画像である第１の画像を少なくとも部分的に変形させる変形処理を実行することで、前記複数の画像のうちの前記第１の画像とは異なる第２の画像に表示された各種部位との位置合わせが施された第３の画像を生成する変形処理手段と、
前記各第２の画像と前記第３の画像の間でサブトラクションすることで、第１のサブトラクション画像を生成する画像生成手段と、
を具備することを特徴とする画像処理装置。
【請求項１２】
前記第１の画像及び前記第２の画像が複数ある場合には、当該複数の第１の画像から前記各第２の画像との前記サブトラクション処理に利用する最適画像を選択する選択手段をさらに具備することを特徴とする請求項１１記載の画像処理装置。
【請求項１３】
前記変形処理手段は、前記第１の画像又は前記第２の画像上の標識となる被検体構造物の位置、形状その他の幾何学的情報に基づいて前記変形処理を実行し、前記第３の画像を生成することを特徴とする請求項１１記載の画像処理装置。
【請求項１４】
前記変形処理手段は、
前記最適画像上に少なくとも一つの第１の関心領域を設定し、
前記各第２の画像上に前記第１の関心領域と対応させた少なくとも一つの第２の領域を設定し、
前記少なくとも一つの第１の関心領域と、前記少なくとも一つの第２の関心領域との間の位置ずれを検出することにより、前記各第２の画像に対する前記最適画像の位置ずれ量を取得し、当該前記位置ずれ量に基づいて、前記変形処理を実行すること、
を特徴とする請求項１１又は１３記載の画像処理装置。
【請求項１５】
前記選択手段は、前記複数の第１の画像から、初期設定又は所定の選択手法により前記最適画像として選択された所定の画像を、他の前記第１の画像に変更するための変更手段を有し、
前記変形処理手段は、変更された前記最適画像を部分的に変形させる変形処理を実行することで、前記各第２の画像に表示された各種部位との位置合わせが施された第３の画像を生成すること、
を特徴とする請求項１１乃至１４のうちいずれか一項記載の画像処理装置。
【請求項１６】
コンピュータに、
被検体への造影剤の注入の前後にわたって連続的に収集された、当該被検体の撮影対象部位に関する複数の画像のうちの一画像である第１の画像を少なくとも部分的に変形させる変形処理を実行することで、前記複数の画像のうちの前記第１の画像とは異なる第２の画像に表示された各種部位との位置合わせが施された第３の画像を生成する第１の機能と、
前記各第２の画像と前記第３の画像の間でサブトラクションすることで、第１のサブトラクション画像を生成する第２の機能と、
を実現させるための画像処理プログラムを記録したコンピュータ読み取り可能な記録媒体。
【請求項１７】
前記第１の画像及び前記第２の画像が複数ある場合には、当該複数の第１の画像から前記各第２の画像との前記サブトラクション処理に利用する最適画像を選択する選択手段をさらに具備することを特徴とする請求項１６記載の画像処理プログラムを記録したコンピュータ読み取り可能な記録媒体。
【請求項１８】
前記変形処理は、前記第１の画像又は前記第２の画像上の標識となる被検体構造物の位置、形状その他の幾何学的情報に基づいて前記変形処理を実行し、前記第３の画像を生成することを特徴とする請求項１６又は１７記載の画像処理プログラムを記録したコンピュータ読み取り可能な記録媒体。
【請求項１９】
前記変形処理は、
前記最適画像上に少なくとも一つの第１の関心領域を設定し、
前記各第２の画像上に前記第１の関心領域と対応させた少なくとも一つの第２の領域を設定し、
前記少なくとも一つの第１の関心領域と、前記少なくとも一つの第２の関心領域との間の位置ずれを検出することにより、前記各第２の画像に対する前記最適画像の位置ずれ量を取得し、当該前記位置ずれ量に基づいて、前記変形処理を実行すること、
を特徴とする請求項１６乃至１８のうちいずれか一項記載の画像処理プログラムを記録したコンピュータ読み取り可能な記録媒体。
【請求項２０】
前記選択ステップにおいては、前記複数の第１の画像から前記最適画像を選択するための複数の選択手法を有し、指定される撮影プログラムに応じて、前記複数の選択手法から選択手法を決定することを特徴とする請求項１６乃至１９のうちいずれか一項記載の画像処理プログラムを記録したコンピュータ読み取り可能な記録媒体。
【請求項２１】
前記選択ステップにおいては、前記複数の第１の画像から、初期設定又は所定の選択手法により前記最適画像として選択された所定の画像を、他の前記第１の画像に変更するための変更処理を実行し、
前記変形処理ステップにおいては、変更された前記最適画像を部分的に変形させる変形処理を実行することで、前記各第２の画像に表示された各種部位との位置合わせが施された第３の画像を生成すること、
を特徴とする請求項１６乃至２０のうちいずれか一項記載の画像処理プログラムを記録したコンピュータ読み取り可能な記録媒体。
【請求項１】
被検体の撮影対象部位に関する複数の画像を、当該被検体への造影剤の注入の前後にわたって連続的に収集する画像収集手段と、
前記複数の画像のうちの一画像である第１の画像を少なくとも部分的に変形させる変形処理を実行することで、前記複数の画像のうちの前記第１の画像とは異なる第２の画像に表示された各種部位との位置合わせが施された第３の画像を生成する変形処理手段と、
前記各第２の画像と前記第３の画像の間でサブトラクションすることで、第１のサブトラクション画像を生成する画像生成手段と、
を具備することを特徴とする画像診断装置。
【請求項２】
前記第１の画像及び前記第２の画像が複数ある場合には、当該複数の第１の画像から前記各第２の画像との前記サブトラクション処理に利用する最適画像を選択する選択手段をさらに具備することを特徴とする請求項１記載の画像診断装置。
【請求項３】
前記変形処理手段は、前記第１の画像又は前記第２の画像上の標識となる被検体構造物の位置、形状その他の幾何学的情報に基づいて前記変形処理を実行し、前記第３の画像を生成することを特徴とする請求項１又は２記載の画像診断装置。
【請求項４】
前記変形処理手段は、
前記最適画像上に少なくとも一つの第１の関心領域を設定し、
前記各第２の画像上に前記第１の関心領域と対応させた少なくとも一つの第２の領域を設定し、
前記少なくとも一つの第１の関心領域と、前記少なくとも一つの第２の関心領域との間の位置ずれを検出することにより、前記各第２の画像に対する前記最適画像の位置ずれ量を取得し、当該前記位置ずれ量に基づいて、前記変形処理を実行すること、
を特徴とする請求項１乃至３のうちいずれか一項記載の画像診断装置。
【請求項５】
前記被検体への造影剤注入開始からの所定時間に基づいて、前記複数の画像を前記複数の第１画像と前記複数の第２画像とに区分する画像区分手段をさらに具備することを特徴とする請求項１乃至４のうちいずれか一項記載の画像診断装置。
【請求項６】
前記撮影対象部位に応じて、前記所定時間、及び当該所定時間内の画像収集レートの少なくとも一方を所望の値に設定するための設定手段をさらに具備し、
前記画像区分手段は、前記設定手段によって前記所定時間が設定された場合には、
設定後の前記所定時間に基づいて、前記複数の第１画像と前記複数の第２画像との区分を実行し、
前記画像収集手段は、前記設定手段によって前記画像収集レートが設定された場合には、設定後の前記画像収集レートに基づいて、前記画像収集を実行すること、
を特徴とする請求項５記載の画像診断装置。
【請求項７】
前記選択手段は、前記複数の第１の画像から前記最適画像を選択するための複数の選択手法を有し、指定される撮影プログラムに応じて、前記複数の選択手法から選択手法を決定することを特徴とする請求項１乃至６のうちいずれか一項記載の画像診断装置。
【請求項８】
前記画像生成手段は、前記各第２の画像から前記最適画像をサブトラクションすることで、第２のサブトラクション画像を生成し、
前記変形処理手段により前記第３の画像が生成されるまでの間、前記画像収集手段によって収集された画像又は前記第２のサブトラクション画像を表示する表示手段をさらに具備すること、
を特徴とする請求項１乃至７のうちいずか一項記載の画像診断装置。
【請求項９】
前記複数の第２の画像と、前記第２の画像毎の前記第３の画像とを記憶する記憶手段をさらに具備し、
前記画像生成手段は、前記記憶手段に記憶された前記各第２の画像と前記第３の画像とを用いて、前記第１のサブトラクション画像を生成すること、
を特徴とする請求項１乃至８のうちいずれか一項記載の画像診断装置。
【請求項１０】
前記複数の第１の画像と、前記複数の第２の画像と、前記第２の画像毎の前記最適画像としての前記第１の画像のフレーム番号と、を記憶する記憶手段をさらに具備し、
前記変形処理手段は、前記記憶手段に記憶された前記各第２の画像と、前記複数の第２の画像と、前記第２の画像毎の前記最適画像としての前記第１の画像のフレーム番号とを用いて前記変形処理を実行することで、前記第３の画像を事後的に生成し、
前記画像生成手段は、前記記憶手段に記憶された前記各第２の画像と事後的に生成された前記第３の画像とを用いて、前記第１のサブトラクション画像を生成すること、
を特徴とする請求項１乃至９のうちいずれか一項記載の画像診断装置。
【請求項１１】
被検体への造影剤の注入の前後にわたって連続的に収集された、当該被検体の撮影対象部位に関する複数の画像を記憶する記憶手段と、
前記複数の画像のうちの一画像である第１の画像を少なくとも部分的に変形させる変形処理を実行することで、前記複数の画像のうちの前記第１の画像とは異なる第２の画像に表示された各種部位との位置合わせが施された第３の画像を生成する変形処理手段と、
前記各第２の画像と前記第３の画像の間でサブトラクションすることで、第１のサブトラクション画像を生成する画像生成手段と、
を具備することを特徴とする画像処理装置。
【請求項１２】
前記第１の画像及び前記第２の画像が複数ある場合には、当該複数の第１の画像から前記各第２の画像との前記サブトラクション処理に利用する最適画像を選択する選択手段をさらに具備することを特徴とする請求項１１記載の画像処理装置。
【請求項１３】
前記変形処理手段は、前記第１の画像又は前記第２の画像上の標識となる被検体構造物の位置、形状その他の幾何学的情報に基づいて前記変形処理を実行し、前記第３の画像を生成することを特徴とする請求項１１記載の画像処理装置。
【請求項１４】
前記変形処理手段は、
前記最適画像上に少なくとも一つの第１の関心領域を設定し、
前記各第２の画像上に前記第１の関心領域と対応させた少なくとも一つの第２の領域を設定し、
前記少なくとも一つの第１の関心領域と、前記少なくとも一つの第２の関心領域との間の位置ずれを検出することにより、前記各第２の画像に対する前記最適画像の位置ずれ量を取得し、当該前記位置ずれ量に基づいて、前記変形処理を実行すること、
を特徴とする請求項１１又は１３記載の画像処理装置。
【請求項１５】
前記選択手段は、前記複数の第１の画像から、初期設定又は所定の選択手法により前記最適画像として選択された所定の画像を、他の前記第１の画像に変更するための変更手段を有し、
前記変形処理手段は、変更された前記最適画像を部分的に変形させる変形処理を実行することで、前記各第２の画像に表示された各種部位との位置合わせが施された第３の画像を生成すること、
を特徴とする請求項１１乃至１４のうちいずれか一項記載の画像処理装置。
【請求項１６】
コンピュータに、
被検体への造影剤の注入の前後にわたって連続的に収集された、当該被検体の撮影対象部位に関する複数の画像のうちの一画像である第１の画像を少なくとも部分的に変形させる変形処理を実行することで、前記複数の画像のうちの前記第１の画像とは異なる第２の画像に表示された各種部位との位置合わせが施された第３の画像を生成する第１の機能と、
前記各第２の画像と前記第３の画像の間でサブトラクションすることで、第１のサブトラクション画像を生成する第２の機能と、
を実現させるための画像処理プログラムを記録したコンピュータ読み取り可能な記録媒体。
【請求項１７】
前記第１の画像及び前記第２の画像が複数ある場合には、当該複数の第１の画像から前記各第２の画像との前記サブトラクション処理に利用する最適画像を選択する選択手段をさらに具備することを特徴とする請求項１６記載の画像処理プログラムを記録したコンピュータ読み取り可能な記録媒体。
【請求項１８】
前記変形処理は、前記第１の画像又は前記第２の画像上の標識となる被検体構造物の位置、形状その他の幾何学的情報に基づいて前記変形処理を実行し、前記第３の画像を生成することを特徴とする請求項１６又は１７記載の画像処理プログラムを記録したコンピュータ読み取り可能な記録媒体。
【請求項１９】
前記変形処理は、
前記最適画像上に少なくとも一つの第１の関心領域を設定し、
前記各第２の画像上に前記第１の関心領域と対応させた少なくとも一つの第２の領域を設定し、
前記少なくとも一つの第１の関心領域と、前記少なくとも一つの第２の関心領域との間の位置ずれを検出することにより、前記各第２の画像に対する前記最適画像の位置ずれ量を取得し、当該前記位置ずれ量に基づいて、前記変形処理を実行すること、
を特徴とする請求項１６乃至１８のうちいずれか一項記載の画像処理プログラムを記録したコンピュータ読み取り可能な記録媒体。
【請求項２０】
前記選択ステップにおいては、前記複数の第１の画像から前記最適画像を選択するための複数の選択手法を有し、指定される撮影プログラムに応じて、前記複数の選択手法から選択手法を決定することを特徴とする請求項１６乃至１９のうちいずれか一項記載の画像処理プログラムを記録したコンピュータ読み取り可能な記録媒体。
【請求項２１】
前記選択ステップにおいては、前記複数の第１の画像から、初期設定又は所定の選択手法により前記最適画像として選択された所定の画像を、他の前記第１の画像に変更するための変更処理を実行し、
前記変形処理ステップにおいては、変更された前記最適画像を部分的に変形させる変形処理を実行することで、前記各第２の画像に表示された各種部位との位置合わせが施された第３の画像を生成すること、
を特徴とする請求項１６乃至２０のうちいずれか一項記載の画像処理プログラムを記録したコンピュータ読み取り可能な記録媒体。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【公開番号】特開２００６−８７６３１（Ｐ２００６−８７６３１Ａ）
【公開日】平成１８年４月６日（２００６．４．６）
【国際特許分類】
【出願番号】特願２００４−２７５９９４（Ｐ２００４−２７５９９４）
【出願日】平成１６年９月２２日（２００４．９．２２）
【出願人】（８００００００３５）株式会社産学連携機構九州 (34)
【出願人】（５９１１０２０９５）三菱スペース・ソフトウエア株式会社 (148)
【出願人】（５９４１６４５４２）東芝メディカルシステムズ株式会社 (4,066)
【Ｆターム（参考）】
【公開日】平成１８年４月６日（２００６．４．６）
【国際特許分類】
【出願日】平成１６年９月２２日（２００４．９．２２）
【出願人】（８００００００３５）株式会社産学連携機構九州 (34)
【出願人】（５９１１０２０９５）三菱スペース・ソフトウエア株式会社 (148)
【出願人】（５９４１６４５４２）東芝メディカルシステムズ株式会社 (4,066)
【Ｆターム（参考）】
[ Back to top ]