Ｘ線診断装置及び画像処理装置
【課題】撮像時間が比較的長時間にわたる例えば血管形態観察目的と心筋灌流目的を兼ねた撮像において、被曝低減を実現すること。
【解決手段】Ｘ線診断装置は、パルスＸ線を発生するＸ線源２と、被検体を透過したＸ線を検出するＸ線検出器５と、Ｘ線検出器の出力に基づいて画像を発生する画像演算記憶部１０と、パルスＸ線を撮像期間内に繰り返し発生させるとともにパルスＸ線の発生周期を撮像期間内で変化させるためにＸ線源を制御するシステム制御部８とを具備する。
【解決手段】Ｘ線診断装置は、パルスＸ線を発生するＸ線源２と、被検体を透過したＸ線を検出するＸ線検出器５と、Ｘ線検出器の出力に基づいて画像を発生する画像演算記憶部１０と、パルスＸ線を撮像期間内に繰り返し発生させるとともにパルスＸ線の発生周期を撮像期間内で変化させるためにＸ線源を制御するシステム制御部８とを具備する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、パルスＸ線により被検体を撮像するＸ線診断装置及びＸ線の心臓冠状動脈造影診断に好適な画像処理装置に関する。
【背景技術】
【０００２】
血管内インターベンション治療直後に、血管造影撮像をおこない、造影剤が心筋へ流れ込むようすの時間変化を観測し、治療終了判断の支援をする手法がある。２０００年にこの手法で予後が予測できるという論文（Gibson）が発表されて以来、この手法で薬効や手技技法の評価をおこなう例が多くなってきている。
【０００３】
この指標は臨床的に有用であることがエビデンスとれているのにもかかわらず、まだルーチンの臨床で毎日使われているには至っていない。原因として、被曝が増大すること、定量性が低いこと、が挙げられている。
【０００４】
血管形態観察を検査目的とする撮像法では一般的に図６のように撮像している。すなわち、心拍運動にして３〜５心拍、撮像時間にして５秒程度だけ、Ｘ線をパルス状に照射して画像を得ている。
【０００５】
一方、心筋灌流の計測を検査目的とする場合には、長時間にわたる撮像が必要になる。少なくとも５秒、一般に３０秒程度、長くて６０秒程度の継続観察が好ましいと考えられる。そこで、血管形態観察を目的とした撮像にこの心筋灌流を目的とした撮像を加えると図７のような撮像法が必要になる。
【０００６】
ここで別の案として、血管形態観察を目的とした撮像と、心筋灌流を目的とした撮像を別々にして合計２回の撮像をおこなうという考え方もある。しかしその場合は、投与する造影剤が２倍になってしまい、最も好ましくない。
【０００７】
また、従来の血管観察用のＸ線診断装置、いわゆるＸ線アンジオ装置においては、冠状動脈造影は施行されるが、造影剤による心筋血流の定量的計測はできなかったので、別室や別の時間帯に核医学診断装置やＭＲＩ装置にて心筋血流を計測する必要があった。
【０００８】
冠状動脈に複数の狭窄が存在してそのいずれかの病変による虚血がある場合に、カテーテル室では心筋血流が計測できないので虚血の部位を簡単に特定できなかった。
【０００９】
冠状動脈の末梢血管の狭窄・塞栓・血栓の存在が疑われる場合に、その確認や血栓溶解剤による治療後の心筋血流回復確認が困難であった。
【００１０】
Ｘ線診断装置を用いて、心筋血流を計測し、画像表示することができなかった。
【非特許文献１】Relationship of TIMI Myocardial Perfusion Grade to Mortality After Administration of Thrombolytic Drugs, C. Michael Gibson, MS, MD; Christopher P. Cannon, MD; Sabina A. Murphy, MPH;Kathryn A. Ryan, BS; Rebecca Mesley, BS; Susan J. Marble, RN, MS; Carolyn H. McCabe, BS;Frans Van de Werf, MD, PhD; Eugene Braunwald, MD;, Circulation, 101, 125-130, 2000
【非特許文献２】Arnoud W.J. van 't Hof, MD; Aylee Liem, MD; Harry Suryapranata, MD; Jan C.A. Hoorntje, MD; Menko-Jan de Boer, MD; Felix Zijlstra, MD;, Angiographic Assessment of Myocardial Reperfusion in Patients Treated With Primary Angioplasty for Acute Myocardial Infarction, Myocardial Blush Grade, Circulation, 97, 2302-2306, 1998）
【発明の開示】
【発明が解決しようとする課題】
【００１１】
本発明の目的は、撮像時間が比較的長時間にわたる例えば血管形態観察目的と心筋灌流目的を兼ねた撮像において、被曝低減と造影剤量低減を実現することにある。
また本発明の目的は、心筋血流に関する有用な情報を提供することにある。
【課題を解決するための手段】
【００１２】
本発明によるＸ線診断装置は、パルスＸ線を発生するＸ線源と、被検体を透過したＸ線を検出するＸ線検出器と、前記Ｘ線検出器の出力に基づいて画像を発生する画像発生部と、前記パルスＸ線を撮像期間内に繰り返し発生させるとともに前記パルスＸ線の発生周期を前記撮像期間内で変化させるために前記Ｘ線源を制御する制御部とを具備する。
本発明に係る画像処理装置は、心臓冠状動脈造影撮影シーケンスのもとで発生された複数の画像のデータを記憶する記憶部と、前記複数の画像のデータに基づいて、心筋への血液供給領域に設定された基準領域に関する基準の時間濃度曲線と、心筋領域内に設定された複数の局所領域に関する複数の時間濃度曲線とを発生し、前記基準の時間濃度曲線に対する前記複数の時間濃度曲線の相関をそれぞれ表す複数の指標を演算する演算部とを具備する。
【発明の効果】
【００１３】
本発明によれば、撮像時間が比較的長時間にわたる例えば血管形態観察目的と心筋灌流目的を兼ねた撮像において、被曝低減と造影剤量低減を実現することができる。
また本発明によれば、心筋血流に関する有用な情報を提供することができる。
【発明を実施するための最良の形態】
【００１４】
（第１実施形態）
以下、図面を参照して、本発明の第１実施形態を説明する。
まず、以下の説明で用いる用語について次の通り定義する。
パルスレート；パルスＸ線を被検体に照射する頻度を規定する指標であり、単位時間（１秒又は１分）又は１心拍あたりのパルスＸ線の照射回数を、単位（回／秒、回／分）又は（回／心拍）で表したものであり、単位時間で撮像する画像枚数を表すフレームレートはパルスレートに実質的に等価である。なお、パルスＸ線の照射には、Ｘ線管球からパルスＸ線を発生してそのまま被検体に照射する形態だけでなく、連続的にＸ線を発生し、Ｘ線シャッター等でパルスＸ線を生成して被検体に照射する形態も含む。
灌流；心筋への血液の流れを表す。
Ｒ波；心電波形のピーク波を表す。
ＲＲ間隔；Ｒ波とＲ波の間の時間間隔を表す。
ｍＡ；Ｘ線管球の電極間に流れる管電流であり、パルスＸ線の高さを表す指標として用いられる。
パルス幅；パルスＸ線の継続時間を表し、単位はｍｓｅｃである。
ｍＡｓ；管電流とパルス幅を掛け算した管電流時間積であり、Ｘ線の強度を表す指標として用いられる。
プロトコル；各種条件をまとめて保存してあるパラメータのセット
ブラッシュ(Blush)；心筋に造影剤が流れ込む程度をGibsonらが半定量化した数値である。
コントロール画像；基準画像
図１は、本実施形態に係るＸ線撮影装置を示している。Ｘ線撮影装置は、ガントリ１００を有する。ガントリ１００は、Ｃアーム７を有する。Ｃアーム７は、支持機構６に回転自在に支持される。Ｃアーム７の一端にはＸ線発生部２が取り付けられる。Ｘ線発生部２は、Ｘ線管球２１とＸ線コリメータ２２を有する。高電圧発生部１は、Ｘ線管球２２の電極間に印加する高電圧（管電圧）を発生し、またＸ線管球２２のフィラメントに供給するフィラメント電流を発生する。高電圧制御部２０は、システム制御部８の制御に従って、高電圧発生部１で発生する管電圧及び／又はフィラメント電流を制御する。
【００１５】
Ｃアーム７の他端にはＸ線検出部５が取り付けられる。Ｘ線検出部５は、寝台４に載置される被検体３を挟んで、Ｘ線発生部２のＸ線管球２２に対峙する。Ｘ線検出部５は、典型的には、入射Ｘ線を直接的又は間接的に電荷に変換する複数の検出素子（画素）が２次元状に配列されてなる固体平面検出器である。Ｘ線検出部５は、システム制御部８の制御により典型的には電荷蓄積、電荷読出及びリセットからなる１サイクルの検出動作を一定周期で繰り返す。
【００１６】
画像演算・記憶部１０は、Ｘ線検出部５からの出力に基づいて画像のデータを発生する機能、画像のデータを記憶する機能、画像のデータを処理する機能を有している。システム制御部８は、インジェクタ１５から被検体３へ造影剤を注入を開始する時点でインジェクタ１５から出力される注入開始信号、被検体への造影剤注入を終了する時点でインジェクタ１５から出力される注入終了信号と、心電計１６で測定された被検体３の心電図（ＥＣＧ）とに基づいて、後述のように撮像動作を制御することを主な機能として有している。システム制御部８には操作部９が接続される。操作部９には、ハンドスイッチ１２、ディスプレイ及びタッチパネル等を有するユーザインタフェース１４が設けられる。
【００１７】
図２にはシステム制御部８による撮像動作を示している。撮像に先だって、被検体に対しカテーテルが挿入され、例えば心臓冠状動脈入り口部まで進められ、その位置で造影検査が開始される。
【００１８】
検査開始とともに、パルスＸ線の発生が開始される。パルスＸ線の発生開始からパルスＸ線の発生終了までの撮像期間は、システム制御部８による撮像制御、より具体的にはパルスＸ線の発生制御の相違により、順番に３つの区間（Ａ）、（Ｂ）、（Ｃ）に区分されることができる。
【００１９】
区間（Ａ）は造影剤がまだ投与されていない期間であり、その時間長は約１〜３心拍程度である。区間（Ａ）では、心電同期のもとでパルスＸ線は不定周期で繰り返し発生される。例えば、区間（Ａ）では、システム制御部８は心電図（ＥＣＧ）のＲ波を検出し、ＲＲ間隔から心臓拡張末期を同定し、その時点でのみパルスＸ線を発生させる。つまり、区間（Ａ）のパルスレートは１回／心拍であり、パルスＸ線は不定周期で繰り返し発生される。
【００２０】
Ｘ線検出部５の検出動作は、区間（Ａ）を含めて、区間（Ａ）、区間（Ｂ）及び区間（Ｃ）の全撮像期間にわたって、クロックに同期して、例えば３０サイクル／秒（３０ｆｐｓ）のフレームレートの逆数の一定周期で繰り返される。
【００２１】
なお、心電同期として且つ一定周期（例えば３０ｆｐｓの逆数）で収集方法を採用してもよい。具体的には、図８に示すように、例えば一定周期（例えば３０ｆｐｓの逆数）で収集している間にＲ波を検出したら、それまで一定周期で照射していたリズムをリセットし、再び一定周期（例えば３０ｆｐｓの逆数）で収集する。この手法によれば、Ｒ波直後のフレームはどの心拍でも心電同期がとれていることになり、且つその後は高速に例えば３０ｆｐｓで収集される。
【００２２】
また、一定の周期（例えば３０ｆｐｓの逆数）で収集する区間において、図９に示すように、心電同期とし且つ一定の周期（例えばＲＲ間隔の１／１０）で収集してもよい。具体的には、例えば一定周期（例えば３０ｆｐｓの逆数）で収集している間にＲ波を検出したとき、それまで一定周期で照射していたリズムをリセットし、ＲＲ間隔を計算してＲＲ間隔の１０分の１の時間を計算し、次の心拍周期では当該計算された周期（一定周期）で収集する。この手法によれば、Ｒ波直後のフレームはどの心拍でも心電同期がとれていることになり、且つその後は高速に収集される。
【００２３】
パルスＸ線の発生周期に対応する画像データが画像演算記憶部１０に記憶される。しかし、パルスＸ線が発生しない周期に対応する画像データは画像演算記憶部１０に記憶されない。画像のデータは、造影前の区間（Ａ）を後述の区間（Ｂ）、（Ｃ）と区別するコードとともに記憶される。
【００２４】
インジェクタ１６のトリガ操作により造影剤の注入が開始される。インジェクタ１６からシステム制御部８に造影剤注入開始信号が供給される。システム制御部８は、造影剤注入開始時点から予め決められた遅延時間Δｔ１を経過した時点で区間（Ａ）を終了し、区間（Ｂ）を開始する。遅延時間Δｔ１は、インジェクタ１６から造影剤注入信号が出力された時点から造影剤がカテーテル先端から出るまでに要する時間、例えば約１秒に設定される。
【００２５】
区間（Ｂ）では、パルスＸ線は一定周期で繰り返し発生される。例えば、区間（Ｂ）では、システム制御部８は、クロックに同期して、検出動作の３０サイクル／秒（３０ｆｐｓ）のフレームレートと同じパルスレート（３０回／秒）の逆数の一定周期で繰り返しパルスＸ線を発生させる。パルスＸ線の時間幅（継続時間）は、Ｘ線検出部５の検出動作における電荷蓄積期間長と等価に設定される。結果的に、パルスＸ線はＸ線検出部５の検出サイクルに同期して発生される。区間（Ｂ）では、画像データが３０ｆｐｓのフレームレートで発生され、全て画像演算記憶部１０に記憶される。画像のデータは、造影期間中の区間（Ｂ）を区間（Ａ）、（Ｃ）と区別するコードとともに記憶される。
【００２６】
インジェクタ１６は予め決められた量の造影剤の注入が完了した時点でシステム制御部８に造影剤注入終了信号を供給する。システム制御部８は、造影剤注入終了時点から予め決められた遅延時間Δｔ２を経過した時点で区間（Ｂ）を終了し、区間（Ｃ）を開始する。遅延時間Δｔ２は、遅延時間Δｔ１より長く、例えばインジェクタ１６から造影剤注入信号が出力された時点から造影剤がカテーテル先端から出るまでに要する時間の２倍の時間（約２秒）に設定される。区間（Ｂ）は、例えば約３〜５心拍程度又は５秒程度の時間である。
【００２７】
区間（Ｃ）は、区間（Ａ）と同様に、心電同期のもとでパルスＸ線は不定周期で繰り返し発生される。区間（Ｃ）でも、システム制御部８は心電図（ＥＣＧ）のＲ波を検出し、ＲＲ間隔から心臓拡張末期を同定し、その時点でのみパルスＸ線を発生させる。区間（Ｃ）のパルスレートは１回／心拍であり、パルスＸ線は不定周期で繰り返し発生される。画像のデータは、造影後の区間（Ｃ）を区間（Ａ）、（Ｂ）と区別するコードとともに記憶される。
【００２８】
なお、上述の説明では、区間（Ａ）、（Ｃ）は、心電同期による不定周期でパルスＸ線を発生しているが、図３に示すように、一定周期でパルスＸ線を繰り返し発生するものであってもよい。典型的には、区間（Ａ）、（Ｃ）のパルスＸ線は、区間（Ｂ）の周期（１／３０秒）より長い一定周期（１秒）で繰り返し発生される。この場合、区間（Ｂ）のパルスレートは、１回／秒となる。
【００２９】
上記パルスレート制御によれば、全期間にわたって高いパルスレートを継続するよりも被曝を低減することができる。しかも、関心部位に造影剤が流通する期間は、比較的高い時間分解能を確保することができる。それ以外の期間は、比較的低い時間分解能であっても、診断に要する心臓拡張末期の画像を収集することができる。
【００３０】
区間（Ｃ）の終了は、画像演算記憶部１０の画像処理結果を受けてシステム制御部８により判定される。例えば、区間（Ｃ）に収集される画像各々から、区間（Ａ）で収集した複数の画像から発生された基準画像（コントロール画像）を差分し、得られた差分画像の全体又は局所領域の濃度合計又は平均濃度（測定値）をリアルタイムで測定する。上記コントロール画像は、区間（Ａ）で収集された複数の画像の平均画像として画像演算記憶部１０により発生され、記憶される。平均処理によりランダムノイズが抑制され得る。
【００３１】
差分画像の測定値は、図４に例示するように、初期値Ｌ0から、撮影部位又は関心部位への造影剤流入に伴って時間経過とともに変化し（図４では低くなり）、造影剤流出に伴って初期値Ｌ0に戻っていく。そこで濃度の初期値Ｌ0に近似するしきい値を、例えば、０．９＜ａ＜１として、ａ・Ｌ0の値に予め設定しておき、測定値が当該しきい値に達した時刻ｔ４で、システム制御部８は高電圧制御部２０にＸ線照射中止信号を出力する。高電圧制御部２０は、システム制御部８からＸ線照射中止信号を受け、高電圧発生部１からＸ線管球２１への管電圧の印加を終了する。なお、区間（Ｃ）が例えば３０秒を経過した場合には、システム制御部８は、測定値がしきい値以下にならなくとも、Ｘ線照射中止信号を出力して、Ｘ線照射を強制的に終了させる。
【００３２】
区間（Ｃ）の終了は、図４に破線で示す濃度時間曲線の傾斜によって判定しても良い。画像演算記憶部１０で繰り返し計算した濃度時間曲線の傾斜が、水平に近似したしきい値に達した時刻ｔ３に区間（Ｃ）は終了する。
【００３３】
なお、画像演算記憶部１０では、様々な画像処理機能を有する。例えば、区間（Ｂ）で得られた画像から心筋灌流に必要な画像を抜き出す。これには区間（Ｂ）で収集した画像から心臓拡張末期の時刻に最も近い画像を各心拍ごとに特定する。区間（Ｂ）では３０ｆｐｓで画像収集されているため、最大でも時間ずれ誤差は３３ｍｓｅｃ以内となる。これは心筋灌流計測には十分に問題のないずれ量である。区間（Ｂ）で特定された複数の画像各々と、区間（Ｃ）で得られた複数の画像から選択された画像とを合成することにより、合成画像を生成する。複数の合成画像から、区間（Ａ）のコントロール画像を差分して、複数の差分画像を生成する。複数の差分画像を区間（Ｂ）内での時間経過に従って配列して表示することにより、造影剤の流出の様子を把握することができる。
【００３４】
なお、上記説明では、造影剤の注入開始／終了により、自動的に区間の切り替え、つまりパルスレート（フレームレート）を変化させるようにしているが、オペレータによる操作部９のパルスレート切り替え指示により、任意の時点でパルスレートを手動で切り替えるようにしてもよい。
【００３５】
上記説明では、３区間としたが、最初の区間（Ａ）をなくして２区間でもよい。この場合コントロール画像は区間（Ｂ）の最初のほうのフレームを用いる。
【００３６】
差分画像を得た後、定量的な評価をする場合がある。この場合、画像上に関心領域（ＲＯＩ）を設定し、その値を読む。
【００３７】
尚、厳密さを要する場合は、血管関心領域と心筋関心領域の間の時間差を考慮するのがよい。すなわち血流は太い血管を流れた後に小血管を得て心筋に達するため、時間差がうまれる。実際には３００ｍｓｅｃ〜１ｓｅｃ程度の時間差がある。病変を有する患者ではこの時間差がさらに長くなる。そこで、太い血管に血管関心領域を設定した場合は、結果、関心領域のデータは、時刻０から（Ｔ−τ）までを積算し、これに対する心筋関心領域は、時刻Ｔのデータを取ることが好ましい。こうした厳密なことをしたい場合には、先の画像収集において１心拍に２枚のデータを取ることが好ましく、すなわち１心拍に２回のＸ線パルス照射をする。
【００３８】
上記定量計測において、さらにノイズを低減したい場合は、近傍フレームで平均処理をするとよい。すなわち図５に示すように、心臓動きが少ない瞬間に、複数、例えば２つのパルスＸ線を一体として発生し、それに伴って収集される２画像を平均した画像から計測を行う。これにより１画像だけの場合よりもノイズが低減される。ただし心臓動きの要素が加算されてしまうので、心臓動きがほとんど無い時間で複数パルスを出すことが必須である。あるいは動き補正をするという手段もある。
【００３９】
この心筋灌流の画像および定量数値を得るのは、空間的に粗くても大丈夫である。そこで例えば生画像が１０２４×１０２４マトリクスの画像であっても、心筋灌流画像は２５６×２５６マトリクスでもよい。そこで心拍同期収集区間では画素加算して画像収集することにより、画像を粗くするかわりに、１画素あたりのＸ線吸収量が少なくててもよいようにする。例えば４画素加算する場合は、Ｘ線量を１／４にする。
【００４０】
心筋灌流を計算するためには、フレーム間でＸ線条件が変わってはならない。Ｘ線条件が変わってしまうと画像の濃度も変わってしまうので、何を測っているのか意味がなくなる。よって提案する収集方法のうち少なくとも１心拍に１枚の画像は、ｋＶおよびｍＡ、ｍｓｅｃを同一にして撮像する。
【００４１】
ここでさらなる工夫としては、ｋＶのみは固定とし、ｍＡｓは変動させてもよい。この場合はフレームごとのｍＡｓ値を記録しておき、心筋灌流画像を作成する場合は記録したｍＡｓ値に基づいて画像を補正してあげる。これが成り立つ理由は単純なＸ線物理式を参照すればよく、ｍＡｓの変動は単純に濃度の変動に比例する。
【００４２】
その他、以下の通り、様々な変形が可能である。
最終区間（Ｃ）の終了判定のために画像濃度値をしきい値と比較したが、所望の灌流計算値をリアルタイムに計算し、計算に十分なデータが得られたら終了判定をするものであってもよい。その他に、区間（Ｃ）において十分なデータが得られたら撮像終了する、近似曲線が描けるのに十分なデータが得られたら撮像終了するものとしてもよい。そのとき、区間（Ｃ）において十分なデータが得られたら撮像を終了してよい旨を通知し、オペレータの終了指示を待機するようにしてもい。当該通知は、メッセージ表示、ピピピ等の音を鳴らすことにより行われる。
【００４３】
画像演算記憶部１０で撮像期間中にリアルタイムで灌流計測値を計算し、表示部１１に表示する。
一定周期でパルスＸ線を発生する区間と、心電同期でパルスＸ線を発生する区間とをオペレータによる手動で切り替えるようにしてもよい。例えば、操作部９に設けられた撮像ボタンが２段階になっていて押し方で切り替わる。サブボタンを有し、メインボタンだけでは一定周期でパルスＸ線を発生するが、メインとサブを同時に押している間は心拍同期でパルスＸ線を発生する。サブボタンは、メインボタンの近傍にあってもよいし、離れていても良い（フットスイッチと手元スイッチ）。最終区間においてはオペレータが撮像ボタンを離したら終わりとする。
【００４４】
心拍同期においては、１心拍につき１照射（１パルス）とするとき、１照射は、心臓収縮末期、心臓拡張末期、心臓拡張中期のいずれかがオペレータにより選択される。
【００４５】
１心拍につき２照射（２パルス）のとき、例えば、第１照射はＲ波から所定後、第２照射はＲ波から異なる所定後とする。所定後とは、Ｒ波からの絶対時間で管理してもよいし、ＲＲ間隔を元にした相対時間で管理しても良い。例えば、第１照射は、心臓収縮末期、拡張末期、拡張中期、のいずれかとする。第２照射は、第１照射より所定時刻だけ前のタイミングとする。所定時刻とは計測点Ａから計測点Ｂまでに造影剤が流れ込むのに要する時間とする
１心拍につき３照射（３パルス）以上とし、３パルスは隣接していて、得られる複数画像は平均化のために利用される。複数心拍につき１照射（１パルス）とする。区間（Ｂ）の心拍同期撮像の時間が所定の時間を超えた場合には、自動的にＫ心拍につき１照射（１パルス）に移行する。
【００４６】
心拍同期の照射においては、Ｘ線エネルギーを同一にする。全てのパルスＸ線で同じｋＶとｍＡｓを用いることを基本とするが、ｋＶを同じとするが、ｍＡｓは異なるものとしてもよい。そのときｍＡｓ値を記録し、後段処理で補正する。
【００４７】
再生時には、見やすいように画像処理して表示する。パルスレートが変化する制御のもとで収集された複数の画像のデータから、一定間隔の画像のデータを生成することを考慮し、区間（Ａ、Ｃ）、（Ｂ）ともに、共通のクロックに対して、パルスＸ線の発生タイミングを同期させる。これによりデータ生成時は必ず同じ心位相のデータが生成できる。
【００４８】
区間（Ａ、Ｃ）においても、区間（Ｂ）と同じ周期で画像を収集してもよい。心電同期による画像列を生成する際には、最も近い時刻の画像を特定して生成する。一定周期の画像群と、心拍間隔の画像群、とを別々に再生する。再生する時は、実時間で再生表示する。
【００４９】
（第２実施形態）
以下図面を参照して第２実施形態を説明する。まず、以下説明で用いる用語について次の通り定義する。
心筋灌流；心筋濃染、心筋パーフュージョン（Perfusion）、心筋ブラッシュ(Blush)、これらは技術的には厳密には異なるが、血液が毛細血管から心筋へ流出入する現象を表現する。
微小灌流(micro circulation)；毛細血管における血流の流れ。
Ｘ線アンジオ装置；Ｘ線診断装置の中の特に血管撮像を主用途とするものである。
Ｘ線画像；被検体を透過したＸ線の強度分布を表す画像であり、Ｘ線画像ともいう。
Ｘ線動画像；２次元検出器によって時間tにわたって繰り返し撮影された一連のＸ線画像のデータセットである。
冠状動脈撮像；造影剤により強調した冠状動脈内腔に関するＸ線画像を撮像することをいう。
冠状動脈；冠状動脈(coronary)
ピクセル；収集されたＸ線画像における画素
フュージョン画像；Ｘ線画像とＸ線画像とをオーバーレイにより合成した画像をいう。
自然対数；Ln
心位相；心電図信号によって検出されたＲ波を目印とし、Ｒ波から次のＲ波までの時間を１００％で規格化し、現在のフレームの時刻を％で表現する。例えば心臓収縮末期の心位相は２５%付近である。
【００５０】
ＲＯＩ；関心領域(Region Of Interest)
カテーテル室；カテーテル検査室、血管造影検査室、血管インターベンション治療室とほぼ同義である。
【００５１】
図１０に本実施形態に係るＸ線診断装置、ここではＸ線アンジオ撮影装置の外観を示し、図１１に機能ブロック図を示す。Ｘ線撮影装置は、ガントリ１００を有する。ガントリ１００は、Ｃアーム７を有する。Ｃアーム７は、機構制御部６に回転自在に支持される。Ｃアーム７の一端にはＸ線発生部２が取り付けられる。Ｘ線発生部２は、Ｘ線管２０とＸ線絞り器２１を有する。高電圧発生部１は、Ｘ線管２０の電極間に印加する高電圧（管電圧）を発生し、またＸ線管２０のフィラメントに供給するフィラメント電流を発生する。高電圧制御部１７は、システム制御部８の制御に従って、高電圧発生部１で発生する管電圧及び／又はフィラメント電流を制御する。
【００５２】
Ｃアーム７の他端にはＸ線検出部５が取り付けられる。Ｘ線検出部５は、検出器１８と、検出器１８の出力を処理して画像データを生成する画像データ生成部１９とを有する。検出器１８は寝台４に載置される被検体３を挟んで、Ｘ線発生部２のＸ線管２０に対峙する。検出器１８は、典型的には、入射Ｘ線を直接的又は間接的に電荷に変換する複数の検出素子（画素）が２次元状に配列されてなる固体平面検出器である。Ｘ線検出部５は、システム制御部８の制御により電荷蓄積、電荷読出及びリセットからなる１サイクルの検出動作を一定周期で繰り返す。被検体３には心電図モニタ端子１５が装着される。心電図モニタ受信部１６は心電図モニタ端子１５からの信号を受信し、被検体３の心電図のデータを発生する。システム制御部８には操作部９が接続される。操作部９には、ハンドスイッチ１２、ディスプレイ及びタッチパネル等を有するユーザインタフェース１４が設けられる。
【００５３】
演算処理部２３は、画像収集保存部１０と画像処理部１１とを有する。画像収集保存部１０は、Ｘ線検出部５から出力される画像データを心位相データを関連付けて記憶する。画像処理部１１は、心臓冠状動脈造影撮影シーケンスのもとで発生された複数の画像のデータに基づいて、心筋への血液供給領域としての冠状動脈上に設定された基準領域に関する時間濃度曲線（基準時間濃度曲線）に対する、心筋領域内に設定された複数の局所領域に関する複数の時間濃度曲線の相関をそれぞれ表す複数の指標を演算するとともに、指標のマップを作成し、マップを冠状動脈像と合成して合成画像を発生する。なお、時間濃度曲線は、元画像の濃度を対象としても良いし、造影前のマスク画像と造影後のコントラスト画像との差分画像上の濃度（造影剤濃度の近似値）を対象としても良い。ここでは、差分画像上の濃度に関する時間濃度曲線として説明する。
【００５４】
本実施形態では、指標として３種類を提供する。
第１の指標Ｋ１は、造影剤注入開始から造影剤注入終了までの造影剤流入期間を対象として、冠状動脈の基準時間濃度曲線を入力関数、心筋局所領域の時間濃度曲線を出力関数とした時の、局所心筋に関する「血液の流入状態を表す指標Ｋ１」として計算される。第２の指標Ｋ２は、造影剤注入終了から撮影完了までの造影剤流出期間を対象として、冠状動脈の基準時間濃度曲線を入力関数、心筋局所領域の時間濃度曲線を出力関数とした時の、局所心筋に関する「血液の流出状態を表す指標Ｋ２」として計算される。第３の指標Ｋ３は、指標Ｋ１と指標Ｋ２とから計算される。尚、この時の時間濃度曲線とは、人体に投与された造影剤が注目部位（本実施形態では冠状動脈と心筋）に分布したり通過したりする際の造影剤量のみに比例した信号成分と定義する、。
【００５５】
図１２、図１３は一般的なＸ線アンジオ撮影で得えられる冠状動脈造影画像（以下ＣＡＧ画像という）を例示している。ＣＡＧ画像では、造影剤のＸ線の吸収線量が大きいので冠状動脈を通過中は冠状動脈の形状が他の組織と区別できる程にコントラストがついて識別できる（図１２参照）。図１３に示すように、カテーテル１０１が冠状動脈１０２まで挿入操作され、その位置でカテーテル１０１から一定時間継続的に造影剤が注入される。Ｘ線撮像は、少なくとも造影剤注入開始前から、造影剤注入終了後所定時間経過までの期間にわたって行われる。
【００５６】
造影剤が冠状動脈にカテーテル１０１を経由して注入されるときに図１０のＸ線診断装置にてＸ線画像を心電図とともに収集する。画像収集後、操作部９を介してオペレータによりＣＡＧ画像上に冠状動脈上の基準領域（心筋血液供給領域）１０３と、心筋上に複数の心筋局所領域１０４とが設定される（図３２Ａ、図３３参照）。
【００５７】
複数の心筋局所領域１０４は、心筋領域上に設定され、典型的には複数の画素を有する。心筋局所領域１０４の造影剤量に相当する濃度は典型的には画素平均値として計算される。しかし、心筋局所領域１０４は単一の画素を有するものであってもよい。心筋血液供給領域１０３は、典型的には、血管と略等価な幅又は少し小さい幅を有する矩形形状を有し、血管に沿って任意の向きに設定され、複数の画素が含まれる。心筋血液供給領域１０３の造影剤量に相当する濃度としては典型的には画素平均値として計算される。心筋血液供給領域１０３は、インジェクターと心筋注目領域の間の流路の任意の部分、より限定するとカテーテルの任意の部分、又はカテーテルの出口（冠状動脈入口に等価）から心筋注目領域の間のの任意の部分に設定される。
【００５８】
画像処理部１１は、心筋血液供給領域１０３に関する時間濃度曲線（Time Density curve；ＴＤＣ）を生成する。同様に、画像処理部１１は、複数の心筋局所給領域１０４に関する複数の時間濃度曲線を生成する。この生成過程の詳細を図３４に示す。
【００５９】
図３２のＣＡＧＸ線画像の座標(x,y)での時刻tでの画素値をCAG(x,y,t)とし、図３４の画素値をこの自然対数をとったｌｎ（CAG(x,y,t)）とする。これにより組織或いは造影剤の単位体積あたりの減弱係数をμ、図２１のＸ線の透過方向のその組織または造成剤の分布厚さをＬとすると、図１４の画素はｌｎ（CAG(x,y,t)）∝μＬの関係がある。２次元画像である図１３画像上で認識される冠状動脈や心筋部画像に分布する造影剤の総量はμＬに比例する。すなわち造影剤の無い時刻t=t0の画像Ｔ0の画素値をCAG(x,y,t0)とすると、
ｌｎ（CAG(x,y,t)）−ｌｎ（CAG(x,y,t0)）
は、図１３のＸ線画像上で認識される冠状動脈や心筋部画像に分布する造影剤の総量に比例して図１４に示す画像Ｔ１からＴ５で示す一連の画像内では相対的に冠状動脈や心筋部画像などに分布する造影剤総量を示している。
【００６０】
すなわち、ＣＡＧ画像は一方向からの透過画像であるので、画像中心付近、通常は画像中心に位置させる心筋の検査対象部分では心筋に対して略垂直にＸ線が当たるため厚み効果はあまり生じないが、画像中心から外れた周縁部分では心筋に対してＸ線がその接線に近い向きで当たるため、厚み効果が強く表れる。つまり、中心部分に対して周縁部分は造影剤濃度の積分距離が長くなる。従って、心臓の周縁部分で測定される造影剤濃度の信頼性は、中心部分で測定される造影剤濃度の信頼性よりも低い。しかし、上述したように、通常は画像中心付近に心筋の検査対象部分を位置合わせするためその部分に関しては高い信頼性で指標を得ることができる。なお、検査対象ではない画像中心から所定距離離れた信頼性の低い心臓の周縁部分に関する指標は表示対象から外してもよい。
【００６１】
なお、図１４の血液供給領域１０３は、イラスト上でわかりやすくするために円で描画したがその形状は図１３で示すようにカテーテル１０１を通過して冠状動脈に流入して動脈血と混合した時のＴＤＣを計測できるように冠状動脈を囲むように例えば矩形に近い形状に設定される。心筋上の局所領域１０４はカテーテル１０１にて造影剤が注入される冠状動脈の支配領域（その冠状動脈を通過した血液が供給されることにより活動している心筋の領域）であって狭窄やその先の毛細血管の微小循環の有無が疑われる心筋部位全体について、例えば関心領域１０４で示すＮ×Ｎピクセル（Ｎ＝１，３，５，．．．奇数）で例えばＮ＝５として画素をサンプリングし、例えば５×５サイズのスムージング処理後以下で述べる演算結果をその画素の中心位置を代表位置として記録する。
【００６２】
［観測］
上記のように、観測される値はｌｎ（CAG(x,y,t)）−ｌｎ（CAG(x,y,t0)）である。ここで心筋血液供給領域１０３に関する血管単位体積あたりの造影剤相対濃度の基準ＴＤＣをCa(t)とし、同様に心筋単位体積あたりの造影剤相対濃度をCmyo(t)とする。すると観測値とCa、Cmyoとの関係は次式で表される。なお、Caは造影剤注入前後の差分値を表しているので、初期値Ca(0)は０値になる。
【数１】
【００６３】
ここで(xa,ya)は図１３の冠状動脈の上流に設定した関心領域１０３内の画素であり、値の平均値、或いは関心領域内の総量を表す。どちらを用いても良いが以下の操作では選択した一方のみを使用する。混在使用はしない。ｔは時刻を表し、単位は秒である。ａは動脈を意味する。
【数２】
【００６４】
ここで(xmyo,ymyo)は図１３に示す心筋に設定した局所の関心領域１０４内の画素であり、値の平均値、或いは関心領域１０４内の総量を表す。複数の関心領域１０４は、図１５Ａ、図１５Ｂに示す操作部９を介してＣＡＧ画像上に設定されたパフュージョンの演算範囲１０５を図３２Ａ、図３３に示すように分割した相互に同一形状及びサイズの複数の局所領域である。なお、添え字のmyoはmyocardium（冠状動脈）を意味する。
【００６５】
［モデル］
次に、心筋への造影剤の流出入モデルとして後述の（理論１）を導入する。これにより次の関係を得る。
【数３】
【００６６】
［観測値のモデルへの適応］
画像Ｔ１からＴ５は一連の心電同期で撮像したＣＡＧ画像のうちの特定の心位相１０６に対応するＣＡＧ画像である。その収集時刻をｔ1,ｔ2,..ｔ5とする。かかる心電同期画像Ｔ１〜Ｔ５のセットについて、画像処理部１１では、［観測］と［モデル］で導かれた（１）（２）（３）式にて記述する処理を行う。
【数４】
【００６７】
画像処理部１１では、心筋血流に比例する実際的な心筋局所領域への血液の流入の様子を反映した指標Ｋ１'を算出する。まず、造影剤がカテーテルから注入開始された時点から注入完了の時点までの造影剤流入期間が、指標Ｋ１の計算対象期間として設定される。図１６に示すように、当該造影剤流入期間内に収集した複数のＸ線画像から、心筋への血液供給領域１０３に関する時間濃度曲線Ｃａ（ｔ）と、複数の心筋局所領域１０３に関する複数の時間濃度曲線Ｃmyo（ｔ）とが生成される。
【００６８】
次に、画像処理部１１では、図１７に示すように、血管供給領域１０３を通過する造影剤の時間積分（造影剤の流入量）∫Ｃa(t)を横軸とし、縦軸に、心筋局所領域１０４の造影剤存在量Ｃmyo（ｔ）として、時刻ごとの値をプロットする。つまり、心筋全体への血液供給量に対する心筋局所における血液の取り込み量の時間変化の離散分布が生成される。この離散分布に対して直線フィティング処理により得られた直線の傾きとして指標（第１指標）Ｋ１'が算出される。この指標Ｋ１'は、心筋全体への血液供給に対して心筋局所がどの程度追従して血液を受け入れているか否かを定量化したものとしての意義を有する。指標Ｋ１'が正常範囲から低値側に外れているとき、冠状動脈の血液流入に対して、心筋局所での血液流入が追従していない、つまり当該心筋局所領域は流入障害を起こしている可能性がある。指標Ｋ１'が正常範囲から高値側に外れているとき、当該心筋局所領域は流出障害を起こしている可能性がある。画像処理部１１では、すべての局所領域１０４について同様の処理により指標Ｋ１'を計算し、指標Ｋ１'のマップを作成する。この得られたマップをＣＡＧ画像に重ねあわせ（フュージョン）て合成画像を生成し、表示部２４に表示することができる。これにより心筋機能を心筋組織との位置関係とともに判断することができる。
【００６９】
変形例Ａ）
なお、上述では一連の心電図同期ＣＡＧ画像Ｔ１からＴ５は、連続した画像から抽出したが、患者への被曝を低減するために、撮像時に特定の心電図位相のみＸ線を発生させるように図１０の装置を制御して後述の（理論１）に必要な画像セットのみを収集することもできる。
変形例Ｂ）
さらに後述の（理論２）によれば造影剤が心筋部に到達してさらに排泄される流出期間における一連のＴＤＣカーブ（図１８）から図１９で示すグラフィプロット法で直線近似により得られたＫ１とＫ２により血流と造影剤の心筋からのバックディフュージョン（心筋から血管に造影剤が戻る量）を計測でき、図１５を用いて前述同様に表示することができる。同様に時間ｔｎは患者への被曝を低減するために、撮像時に特定の心電図位相のみＸ線を発生させるように図１０の装置を制御して理論１に必要な画像セットのみを収集することもできる。
【００７０】
変形例Ｃ）
さらにＫ1とＫ２の比（Ｋ１／Ｋ２）を計算して表示することにより、「C. Michael Gibson, MS, MD; Christopher P. Cannon, MD; Sabina A. Murphy, MPH; Kathryn A. Ryan, BS; Rebecca Mesley, BS; Susan J. Marble, RN, MS; Carolyn H. McCabe, BS; Frans Van de Werf, MD, PhD; Eugene Braunwald, MD; for the TIMI (Thrombolysis In Myocardial Infarction) Study Group, Relationship of TIMI Myocardial Perfusion Grade to Mortality After Administration of Thrombolytic Drugs, Circulation. 2000;101:125-130.」に示された分類法を定量化（自動化）することが可能になる。すなわち図２８の分類表に従って自動分類したり、図２９Ａ、図２９Ｂのような指標Ｋ１、Ｋ２の関係を表すグラフにプロットして表現することが可能になる。
【００７１】
変形例Ｄ）
本実施形態により、同一患者、同一方向の一連の心電図同期ＣＡＧ画像を二種類収集する。すなわち、図２０Ａ、図２０Ｂに示すように、安静時の心筋血流画像２０５（＝Ｐrest(x,y)）と例えばアデノシンのような心筋血流を増加させる効果を有する薬剤を投与後の薬剤負荷時の心筋血流画像２０６（＝Ｐstress(x,y)）を収集する。安静時と負荷時の比較画像として（心筋血流画像２０６）／（心筋血流画像２０５）＝Ｐstress(x,y)／Ｐrest(x,y)を定義すると、心筋血流画像２０５，２０６で同一方向の収集された同一位置(x,y)では撮像方向の心筋分厚さは同じなので割り算により近似的にキャンセルされ、単位心筋体積あたりの血流増加比（心筋Flow reserve）を得ることができる。
【数５】
【００７２】
また同様に、狭窄等の治療前の指標Ｋ'beforeと治療後の指標Ｋ'afterを計測し、割り算すると、単位心筋体積あたりの血流増加比（治療による回復比）を得ることができる。
【００７３】
変形例Ｅ）
また実施形態によれば、図２７に例示するように、後述の［理論１、２］を完全に演算しなくても、Ｃmyo（ｔ）＝ｌｎ（CAG(x,y,t)−ｌｎ（CAG(x,y,t0)）とおいたときに、Ｃmyo（ｔ）−Ｃmyo（ｔ０）の画像は心筋への血流を示しており、この差分動画像を表示することも同一検査時に限って可能である。ここで、ｔ０＝造影剤注入直前である。
【００７４】
（理論１）
時刻ｔにおける心筋部の造影剤相対濃度Cmyo(t)、冠状動脈血液中の造影剤相対濃度時間曲線をCa(t)、心筋流入血流をＫＫ１、心筋流出血流をＫＫ２と表す。図２３において、心筋部と冠状動脈血中の造影剤の質量バランスは式（１）で表現できる
【数６】
【００７５】
心筋部の血流あるいは血流と相関のあるパラメータを以下に述べる造影剤が心筋部へ流入する現象を利用するPatlak plot法で計算することができる。造影剤が心筋部へ流入し始める条件下（0≦t≦MTT）（MTT＝mean tangent time＝造影剤が動脈から心筋に入り静脈に出てくるまでの平均通過時間、心筋では通常５〜１０ｓｅｃ程度）では、流出量ｋ２Ｃｍｙｏ(t)は極めて少なく、Ｃｍｙｏ(t)≒０、Ｋ１Ｃa(t)＞＞Ｋ２Ｃｍｙｏ(t)の条件が成立し、このときEq.(1) は以下のように簡略化できる。
【数７】
【００７６】
Ｋ１は造影剤が血液から心筋間質部へ移行する移行常数（sec-1）或いは (ml/min/g)を表し、血流に比例する。
Rutland MD. A single injection technique for subtraction of blood background in 131I-hippuran renograms. Br J Radiol 1979;52:134-137.
Patlak CS, Blasberg RG, Fenstermacher JD. Graphical evaluation of blood-to-brain transfer constants from multiple-time uptake data. J Cereb Blood Flow Metab 1983;3:1-7.
（理論２）
時刻ｔにおける心筋部の造影剤相対濃度をCmyo(t)、冠状動脈血液中の造影剤相対濃度をCa(t)と表す。図１０において、心筋部と冠状動脈血中の造影剤の質量バランスは式（１）で表現できる
【数８】
【００７７】
心筋部の血流あるいは血流と相関のあるパラメータを以下に述べる造影剤が心筋から流出するクリヤランスを利用する方法で計算することができる。
【００７８】
Eq. (1)を積分し(2)を得る。
【数９】
【００７９】
Ｙ軸にY(t) Ｘ軸にX(t) をプロットすると直線の傾きがＫ２、Ｙ軸の切片はＫ１を表す。Ｋ１は造影剤が血液から心筋間質部へ移行する移行常数で（sec-1）或いは (ml/min/g)を表し血流に比例する。Ｋ２は造影剤が心筋間質部から血液へ移行する移行常数（sec-1）を表す（以下文献を図３１とともに参照）。
【００８０】
Yokoi T, Iida H, Itoh H, Kanno I. A new graphic plot analysis for cerebral blood flow and partition coefficient with iodine-123-iodoamphetamine and dynamic SPECT validation studies using oxygen-15-water and PET. J Nucl Med 1993;34:498-505.
（実際の処理フロー）
以下、実際の画像処理フローの一例を説明する。ここに示す例は一例であり、前述の理論に基づく他の処理方法であってもよい。
図３２を用いて、本実施形態の全体の具体的処理フローを説明する。Ｘ線画像は血管造影された被写体を写した動画像であることを特徴とする（Ｓ１１）。なお、例えば動画は３〜６０秒の継続時間を有する。１秒あたり１０〜３０枚の画像を有する動画像である。画像は同一方向、同一Ｘ線条件で撮像されたものである。撮像中には、患者の動き、寝台の動き、は無いものとする。
【００８１】
画像処理部１１では１心拍ごとに特定心位相に対応する１枚の画像を抽出し（Ｓ１４）、その画像を造影剤注入前の同心位相の画像から差分処理する（Ｓ１２）。具体的には、造影剤注入前の同心位相が複数枚ある場合は、平均した画像を差分処理する。
【００８２】
画像処理部１１では動き補償処理する（Ｓ１３）。画像パターン（血管、心壁、カテ）を基準とした、パターンマッチング処理により動き補償量を検出する。特徴的なパターンの無いフレームでは、他時刻同位相で得られた移動量で代用する。動き補償は、画像全体を動かす、もしくはＲＯＩ移動位置がわかるテーブルを保有する。動き補償では、患者の多少の体動き、呼吸運動、も補償される。
【００８３】
特定の心位相は、操作部９を介して操作者により予め指定される。例えば心臓拡張末期に指定される。パルス撮像された多フレームから、最も近い時刻のフレームを抜き出してくる。なお、血液供給領域１０３の濃度時間曲線Caと心筋局所領域１０４の濃度時間曲線Cmyoでは異なる心位相を用いてもよい。Ca(t)に対しては、Cmyo(t+T)を用いる。ここでＴは、血液供給領域１０３から心筋局所領域１０４まで造影剤が流れていくのに要する遅延値である。具体的には１〜１０フレーム（１／３０〜１／３ｓｅｃ）程度であることが多い。この遅延時間は、一般にTIMI frame count(TFC)もしくはcorrected TIMI frame count(CTFC)と呼ばれる。
【００８４】
画像処理部１１では、濃度時間積分Σca(t)に対するCmyo(t)の相関を表す傾斜としての指標Ｋ１を計算する（Ｓ１５）。図１７に例示したようにフィッティングして回帰直線を計算する。回帰は所定時間区間のみで計算させ、この所定区間は造影剤インジェクターと連動して決定する。すなわちインジェクターからの造影剤注入開始時刻をJ1、注入終了時刻をJ2とすると、所定区間の開始時刻t1は例えばt1=J1+1sec、所定区間の終了は例えばJ2+2secとする。また、所定区間は画像の値から判定するものであってもよい。すなわちＴＤＣが立ちあがる時刻をt1とし、ＴＤＣがフラットになる時刻を所定区間の終了と定める。回帰直線との相関値を計算する。不整脈を含む心拍は計算から除外することが好ましい。ひとつの動画像に対してひとつのＫ１を計算する。ただし高速化を目的として、フレームごとに随時Ｋ１を計算する手法でもよい。上述したように、画素ごとではなく、周囲N×M画素の平均値（Medianでもよい）をその画素の値とする。表示部２４は血流画像を表示する（Ｓ１７）。指標Ｋ１のマップを表示する。Ｋ１のマップは濃度値をカラー値に変換して表示してもよい。収集した元の動画像を表示し、その直後にＫ１のマップを重ねて表示しても良いし、元画像とＫ１のマップを左右に並べて表示してもよい。
【００８５】
上述した変形例Ａに対応して、本実施形態は、画像処理を、Ｘ線診断装置に一体化してもよいし（図１０、図２２）、もしくはＸ線診断装置から独立させた画像処理装置（図３０、図２３）として提供しても良い。血流画像（指標画像、マップ画像）に、元画像（血管画像、撮像された画像）を重ねて表示する（Ｓ１６）。指標をそのままカラーテーブルにあてはめて得られるカラー血流画像に、白黒の元画像（血管画像）を重ねて表示する。なお、元画像は動画であり、血管動画像に血流カラー静止画像を重ねるのもよい。実際的には、血管静止画像に血流カラー静止画像を重ねる。動画から静止画像を作成するには、動画のうち最も濃い１フレームを選択するか、各画素ごろに動画の画素値の最も小さい値を検出して画像とする。元画像とは、差分前の画像、もしくは差分後の画像とする。グラフを表示することを特長とする（Ｓ２１）。特定のＲＯＩ（x,y）について、図１６〜１９のグラフを表示する。図１７、図１８のグラフには、フィッティング関数（回帰直線）の傾きと切片と相関値を表示する。血液供給領域１０３はオペレーターが指定する。血液供給領域１０３を画像上に指定のためのグラフィカルユーザインタフェース１４を備える。１心拍の特定の心位相で１パルスＸ線を照射する（図２３）。
【００８６】
画像処理部１１は指標Ｋ１とＫ２を計算する（図２４）。指標Ｋ１に対するＫ２の比（Ｋ２）／Ｋ１を計算する（図２５Ａ）。画像処理部１１はＫ１とＫ２の比のマップを生成する。
【００８７】
また、画像処理部１１は、図２５Ｂに示すように、Ｋ１を横軸に、Ｋ２を縦軸に配置して、局所領域毎にＫ１とＫ２に対応する点をプロットする。分布上では、「正常」、「流入障害」、「流出障害」、「流出入障害」のような心筋機能を表す４種の区分に、各心筋局所領域を分類され得る。つまり、心筋局所領域を指標Ｋ１を第１閾値ＴＨ１に対する高低の区別と、指標Ｋ２の第２閾値ＴＨ２に対する高低の区別との組み合わせに分類する。これにより心筋機能を局所ごとに判定することができる。
【００８８】
さらに、画像処理部１１では、図３２Ｂに示すように４つの区分に４つのコード（１〜４）を割り当てて、コードを図３２Ａに示すように心筋局所領域に当てはめて元画像に重ねて表示するようにしてもよいし、図３３に示すように、４つの区分を異なる表示態様、例えば輝度又は色相で区別して元画像に重ねて表示するようにしてもよい。この表示により、心筋機能障害の状態をその区別とともに空間的に把握することができる。
また、区分は、図２８、図２９Ｂに示したように分類してもよい。
【００８９】
さらに画像処理部１１は、比較解析処理を装備する（図２６）。例えば同一手術中の時刻Aと時刻B（C,D,E,,,）の指標Ｋ１をそれぞれ計算し左右に並べて表示する。もしくは、時間的に切り替えて表示する。時刻Aと時刻B（C,D,E,,,）のＫ１の比を計算して、比を画像表示する。比のカラー画像に、元画像を重ねて表示する。例えば治療前後での比のカラー画像を、血管が描出されている元画像に重ねて表示することにより、心筋のどの領域は治療によって血流回復が得られたかが可視化されてわかるようになる。
【００９０】
以上説明した本実施形態の特徴としては以下の通りである。
Ｘ線ＣＴ装置から出されている造影剤を用いたPerfusion計算方法との差異としては、画像処理がＣＴとＸ線では大きく異なり、Ｘ線特有の処理が必要となる。Ｘ線は動画像であるため、特有の工夫、具体的には、１心拍に１枚の画像の抽出、背景差分、動き補償、厚み方向の処理、血管ROIの設定等の特有の処理を必要とした。
【００９１】
Ｘ線での時間濃度曲線の傾きを見る公知方法との差異としては、[理論１、２]で述べた考え方そのものが異なる。指標Ｋ１においては横軸を流入量としたグラフを描く考え方が従来にはない。これにより、後段の解析や、結果画像が異なってくる。処理方法は時間濃度曲線の傾きを得る処理とはまったく異なる。
【００９２】
文献の（或いは図２８の）MBG（TMP）手法は、画像数値じたいには意味を持っていなかった。本実施形態では画像数値は心筋への流れ込み流速であるという意味を持つ点で、MBGより明らかに有用性が高い。流入量でNormalizeしているので、注入する造影剤量に左右されず、臨床の場でも定量性が得られる。
【００９３】
目視ではないので定量性がある。予後予測できるので、この値を参考にしてインターベンション治療を終わって良いかどうかの判断ができる。処置（治療）前後で比較できるので、どれだけ流量が増えたか、すなわち心筋への血液供給がどれだけ改善したかが数値でわかる。処理自動化することにより、臨床でも簡単に使えるようになる。
【００９４】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【図面の簡単な説明】
【００９５】
【図１】本発明の第１実施形態に係るＸ線診断装置の構成を示す図。
【図２】図１のシステム制御部による撮像動作の一例を示す図。
【図３】図１のシステム制御部による他の撮像動作を示す図。
【図４】図１のシステム制御部による撮像終了判定処理を補足するための時間濃度曲線を示す図。
【図５】図１のシステム制御部によるさらに他の撮像動作を示す図。
【図６】従来の撮像動作を示す図。
【図７】従来の撮像動作を示す図。
【図８】図１のシステム制御部による他の撮像動作を示す図。
【図９】図１のシステム制御部による他の撮像動作を示す図。
【図１０】本発明の第２実施形態によるＸ線診断装置の構成を示す図。
【図１１】図１０の特定部分の詳細図。
【図１２】冠状動脈造影画像の一例を示す図。
【図１３】本実施形態において、冠状動脈造影画像上に設定されるＲＯＩを示す図。
【図１４】本実施形態において、画像収集シーケンスを示す図。
【図１５Ａ】本実施形態において、パフュージョンの演算範囲を示す図。
【図１５Ｂ】本実施形態において、ＣＡＧ画像上に設定操作されたパフュージョンの演算範囲を示す図。
【図１６】図１１の画像処理部により生成される造影剤流入期間における時間濃度曲線を示す図。
【図１７】図１１の画像処理部により計算される指標Ｋ１を示す図。
【図１８】図１１の画像処理部により生成される撮像期間における時間濃度曲線を示す図。
【図１９】図１１の画像処理部により計算される指標Ｋ２を示す図。
【図２０Ａ】本実施形態において、安静時と負荷時のパフュージョンの演算範囲を示す図。
【図２０Ｂ】図２の画像処理部により生成される安静画像と負荷画像との比較画像を示す図。
【図２１】本実施形態において、透過画像の厚み効果の補足図。
【図２２】本実施形態の処理手順を示す流れ図。
【図２３Ａ】本実施形態の他の処理手順を示す流れ図。
【図２３Ｂ】本実施形態の他の処理手順を示す流れ図。
【図２４】本実施形態の他の処理手順を示す流れ図。
【図２５Ａ】本実施形態の他の処理手順を示す流れ図。
【図２５Ｂ】本実施形態の２種の指標の分布を示す図。
【図２６】本実施形態の他の処理手順を示す流れ図。
【図２７】本実施形態の他の処理手順を示す流れ図。
【図２８】図１１の画像処理部による２種の指標を用いた局所領域の分類のための分類表を示す図。
【図２９Ａ】図１１の画像処理部により生成される２種の指標の分布を示す図。
【図２９Ｂ】図２８の分類表に対応する指標分布上の区分を示す図。
【図３０】本実施形態において、Ｘ線診断装置から独立させた画像処理装置を示す図。
【図３１】本実施形態において、指標Ｋ１、Ｋ２の説明図。
【図３２Ａ】図１１の画像処理部により２種の指標Ｋ１、Ｋ２から生成した心筋機能マップを示す図。
【図３２Ｂ】図３２Ａのマップコードを示す図。
【図３３】図３２Ａの他の心筋機能マップを示す図。
【符号の説明】
【００９６】
１…高電圧発生部、２…Ｘ線発生部、３…被検体、４…寝台、５…Ｘ線検出部、６…機構制御部、７…Ｃアーム、８…システム制御部、９…操作部、１０…画像演算記憶部、１１…表示部、１２…ハンドスイッチ、１３…スイッチホルダ、１４…ユーザインターフェース、１５…インジェクター、２０…高電圧制御部。
【技術分野】
【０００１】
本発明は、パルスＸ線により被検体を撮像するＸ線診断装置及びＸ線の心臓冠状動脈造影診断に好適な画像処理装置に関する。
【背景技術】
【０００２】
血管内インターベンション治療直後に、血管造影撮像をおこない、造影剤が心筋へ流れ込むようすの時間変化を観測し、治療終了判断の支援をする手法がある。２０００年にこの手法で予後が予測できるという論文（Gibson）が発表されて以来、この手法で薬効や手技技法の評価をおこなう例が多くなってきている。
【０００３】
この指標は臨床的に有用であることがエビデンスとれているのにもかかわらず、まだルーチンの臨床で毎日使われているには至っていない。原因として、被曝が増大すること、定量性が低いこと、が挙げられている。
【０００４】
血管形態観察を検査目的とする撮像法では一般的に図６のように撮像している。すなわち、心拍運動にして３〜５心拍、撮像時間にして５秒程度だけ、Ｘ線をパルス状に照射して画像を得ている。
【０００５】
一方、心筋灌流の計測を検査目的とする場合には、長時間にわたる撮像が必要になる。少なくとも５秒、一般に３０秒程度、長くて６０秒程度の継続観察が好ましいと考えられる。そこで、血管形態観察を目的とした撮像にこの心筋灌流を目的とした撮像を加えると図７のような撮像法が必要になる。
【０００６】
ここで別の案として、血管形態観察を目的とした撮像と、心筋灌流を目的とした撮像を別々にして合計２回の撮像をおこなうという考え方もある。しかしその場合は、投与する造影剤が２倍になってしまい、最も好ましくない。
【０００７】
また、従来の血管観察用のＸ線診断装置、いわゆるＸ線アンジオ装置においては、冠状動脈造影は施行されるが、造影剤による心筋血流の定量的計測はできなかったので、別室や別の時間帯に核医学診断装置やＭＲＩ装置にて心筋血流を計測する必要があった。
【０００８】
冠状動脈に複数の狭窄が存在してそのいずれかの病変による虚血がある場合に、カテーテル室では心筋血流が計測できないので虚血の部位を簡単に特定できなかった。
【０００９】
冠状動脈の末梢血管の狭窄・塞栓・血栓の存在が疑われる場合に、その確認や血栓溶解剤による治療後の心筋血流回復確認が困難であった。
【００１０】
Ｘ線診断装置を用いて、心筋血流を計測し、画像表示することができなかった。
【非特許文献１】Relationship of TIMI Myocardial Perfusion Grade to Mortality After Administration of Thrombolytic Drugs, C. Michael Gibson, MS, MD; Christopher P. Cannon, MD; Sabina A. Murphy, MPH;Kathryn A. Ryan, BS; Rebecca Mesley, BS; Susan J. Marble, RN, MS; Carolyn H. McCabe, BS;Frans Van de Werf, MD, PhD; Eugene Braunwald, MD;, Circulation, 101, 125-130, 2000
【非特許文献２】Arnoud W.J. van 't Hof, MD; Aylee Liem, MD; Harry Suryapranata, MD; Jan C.A. Hoorntje, MD; Menko-Jan de Boer, MD; Felix Zijlstra, MD;, Angiographic Assessment of Myocardial Reperfusion in Patients Treated With Primary Angioplasty for Acute Myocardial Infarction, Myocardial Blush Grade, Circulation, 97, 2302-2306, 1998）
【発明の開示】
【発明が解決しようとする課題】
【００１１】
本発明の目的は、撮像時間が比較的長時間にわたる例えば血管形態観察目的と心筋灌流目的を兼ねた撮像において、被曝低減と造影剤量低減を実現することにある。
また本発明の目的は、心筋血流に関する有用な情報を提供することにある。
【課題を解決するための手段】
【００１２】
本発明によるＸ線診断装置は、パルスＸ線を発生するＸ線源と、被検体を透過したＸ線を検出するＸ線検出器と、前記Ｘ線検出器の出力に基づいて画像を発生する画像発生部と、前記パルスＸ線を撮像期間内に繰り返し発生させるとともに前記パルスＸ線の発生周期を前記撮像期間内で変化させるために前記Ｘ線源を制御する制御部とを具備する。
本発明に係る画像処理装置は、心臓冠状動脈造影撮影シーケンスのもとで発生された複数の画像のデータを記憶する記憶部と、前記複数の画像のデータに基づいて、心筋への血液供給領域に設定された基準領域に関する基準の時間濃度曲線と、心筋領域内に設定された複数の局所領域に関する複数の時間濃度曲線とを発生し、前記基準の時間濃度曲線に対する前記複数の時間濃度曲線の相関をそれぞれ表す複数の指標を演算する演算部とを具備する。
【発明の効果】
【００１３】
本発明によれば、撮像時間が比較的長時間にわたる例えば血管形態観察目的と心筋灌流目的を兼ねた撮像において、被曝低減と造影剤量低減を実現することができる。
また本発明によれば、心筋血流に関する有用な情報を提供することができる。
【発明を実施するための最良の形態】
【００１４】
（第１実施形態）
以下、図面を参照して、本発明の第１実施形態を説明する。
まず、以下の説明で用いる用語について次の通り定義する。
パルスレート；パルスＸ線を被検体に照射する頻度を規定する指標であり、単位時間（１秒又は１分）又は１心拍あたりのパルスＸ線の照射回数を、単位（回／秒、回／分）又は（回／心拍）で表したものであり、単位時間で撮像する画像枚数を表すフレームレートはパルスレートに実質的に等価である。なお、パルスＸ線の照射には、Ｘ線管球からパルスＸ線を発生してそのまま被検体に照射する形態だけでなく、連続的にＸ線を発生し、Ｘ線シャッター等でパルスＸ線を生成して被検体に照射する形態も含む。
灌流；心筋への血液の流れを表す。
Ｒ波；心電波形のピーク波を表す。
ＲＲ間隔；Ｒ波とＲ波の間の時間間隔を表す。
ｍＡ；Ｘ線管球の電極間に流れる管電流であり、パルスＸ線の高さを表す指標として用いられる。
パルス幅；パルスＸ線の継続時間を表し、単位はｍｓｅｃである。
ｍＡｓ；管電流とパルス幅を掛け算した管電流時間積であり、Ｘ線の強度を表す指標として用いられる。
プロトコル；各種条件をまとめて保存してあるパラメータのセット
ブラッシュ(Blush)；心筋に造影剤が流れ込む程度をGibsonらが半定量化した数値である。
コントロール画像；基準画像
図１は、本実施形態に係るＸ線撮影装置を示している。Ｘ線撮影装置は、ガントリ１００を有する。ガントリ１００は、Ｃアーム７を有する。Ｃアーム７は、支持機構６に回転自在に支持される。Ｃアーム７の一端にはＸ線発生部２が取り付けられる。Ｘ線発生部２は、Ｘ線管球２１とＸ線コリメータ２２を有する。高電圧発生部１は、Ｘ線管球２２の電極間に印加する高電圧（管電圧）を発生し、またＸ線管球２２のフィラメントに供給するフィラメント電流を発生する。高電圧制御部２０は、システム制御部８の制御に従って、高電圧発生部１で発生する管電圧及び／又はフィラメント電流を制御する。
【００１５】
Ｃアーム７の他端にはＸ線検出部５が取り付けられる。Ｘ線検出部５は、寝台４に載置される被検体３を挟んで、Ｘ線発生部２のＸ線管球２２に対峙する。Ｘ線検出部５は、典型的には、入射Ｘ線を直接的又は間接的に電荷に変換する複数の検出素子（画素）が２次元状に配列されてなる固体平面検出器である。Ｘ線検出部５は、システム制御部８の制御により典型的には電荷蓄積、電荷読出及びリセットからなる１サイクルの検出動作を一定周期で繰り返す。
【００１６】
画像演算・記憶部１０は、Ｘ線検出部５からの出力に基づいて画像のデータを発生する機能、画像のデータを記憶する機能、画像のデータを処理する機能を有している。システム制御部８は、インジェクタ１５から被検体３へ造影剤を注入を開始する時点でインジェクタ１５から出力される注入開始信号、被検体への造影剤注入を終了する時点でインジェクタ１５から出力される注入終了信号と、心電計１６で測定された被検体３の心電図（ＥＣＧ）とに基づいて、後述のように撮像動作を制御することを主な機能として有している。システム制御部８には操作部９が接続される。操作部９には、ハンドスイッチ１２、ディスプレイ及びタッチパネル等を有するユーザインタフェース１４が設けられる。
【００１７】
図２にはシステム制御部８による撮像動作を示している。撮像に先だって、被検体に対しカテーテルが挿入され、例えば心臓冠状動脈入り口部まで進められ、その位置で造影検査が開始される。
【００１８】
検査開始とともに、パルスＸ線の発生が開始される。パルスＸ線の発生開始からパルスＸ線の発生終了までの撮像期間は、システム制御部８による撮像制御、より具体的にはパルスＸ線の発生制御の相違により、順番に３つの区間（Ａ）、（Ｂ）、（Ｃ）に区分されることができる。
【００１９】
区間（Ａ）は造影剤がまだ投与されていない期間であり、その時間長は約１〜３心拍程度である。区間（Ａ）では、心電同期のもとでパルスＸ線は不定周期で繰り返し発生される。例えば、区間（Ａ）では、システム制御部８は心電図（ＥＣＧ）のＲ波を検出し、ＲＲ間隔から心臓拡張末期を同定し、その時点でのみパルスＸ線を発生させる。つまり、区間（Ａ）のパルスレートは１回／心拍であり、パルスＸ線は不定周期で繰り返し発生される。
【００２０】
Ｘ線検出部５の検出動作は、区間（Ａ）を含めて、区間（Ａ）、区間（Ｂ）及び区間（Ｃ）の全撮像期間にわたって、クロックに同期して、例えば３０サイクル／秒（３０ｆｐｓ）のフレームレートの逆数の一定周期で繰り返される。
【００２１】
なお、心電同期として且つ一定周期（例えば３０ｆｐｓの逆数）で収集方法を採用してもよい。具体的には、図８に示すように、例えば一定周期（例えば３０ｆｐｓの逆数）で収集している間にＲ波を検出したら、それまで一定周期で照射していたリズムをリセットし、再び一定周期（例えば３０ｆｐｓの逆数）で収集する。この手法によれば、Ｒ波直後のフレームはどの心拍でも心電同期がとれていることになり、且つその後は高速に例えば３０ｆｐｓで収集される。
【００２２】
また、一定の周期（例えば３０ｆｐｓの逆数）で収集する区間において、図９に示すように、心電同期とし且つ一定の周期（例えばＲＲ間隔の１／１０）で収集してもよい。具体的には、例えば一定周期（例えば３０ｆｐｓの逆数）で収集している間にＲ波を検出したとき、それまで一定周期で照射していたリズムをリセットし、ＲＲ間隔を計算してＲＲ間隔の１０分の１の時間を計算し、次の心拍周期では当該計算された周期（一定周期）で収集する。この手法によれば、Ｒ波直後のフレームはどの心拍でも心電同期がとれていることになり、且つその後は高速に収集される。
【００２３】
パルスＸ線の発生周期に対応する画像データが画像演算記憶部１０に記憶される。しかし、パルスＸ線が発生しない周期に対応する画像データは画像演算記憶部１０に記憶されない。画像のデータは、造影前の区間（Ａ）を後述の区間（Ｂ）、（Ｃ）と区別するコードとともに記憶される。
【００２４】
インジェクタ１６のトリガ操作により造影剤の注入が開始される。インジェクタ１６からシステム制御部８に造影剤注入開始信号が供給される。システム制御部８は、造影剤注入開始時点から予め決められた遅延時間Δｔ１を経過した時点で区間（Ａ）を終了し、区間（Ｂ）を開始する。遅延時間Δｔ１は、インジェクタ１６から造影剤注入信号が出力された時点から造影剤がカテーテル先端から出るまでに要する時間、例えば約１秒に設定される。
【００２５】
区間（Ｂ）では、パルスＸ線は一定周期で繰り返し発生される。例えば、区間（Ｂ）では、システム制御部８は、クロックに同期して、検出動作の３０サイクル／秒（３０ｆｐｓ）のフレームレートと同じパルスレート（３０回／秒）の逆数の一定周期で繰り返しパルスＸ線を発生させる。パルスＸ線の時間幅（継続時間）は、Ｘ線検出部５の検出動作における電荷蓄積期間長と等価に設定される。結果的に、パルスＸ線はＸ線検出部５の検出サイクルに同期して発生される。区間（Ｂ）では、画像データが３０ｆｐｓのフレームレートで発生され、全て画像演算記憶部１０に記憶される。画像のデータは、造影期間中の区間（Ｂ）を区間（Ａ）、（Ｃ）と区別するコードとともに記憶される。
【００２６】
インジェクタ１６は予め決められた量の造影剤の注入が完了した時点でシステム制御部８に造影剤注入終了信号を供給する。システム制御部８は、造影剤注入終了時点から予め決められた遅延時間Δｔ２を経過した時点で区間（Ｂ）を終了し、区間（Ｃ）を開始する。遅延時間Δｔ２は、遅延時間Δｔ１より長く、例えばインジェクタ１６から造影剤注入信号が出力された時点から造影剤がカテーテル先端から出るまでに要する時間の２倍の時間（約２秒）に設定される。区間（Ｂ）は、例えば約３〜５心拍程度又は５秒程度の時間である。
【００２７】
区間（Ｃ）は、区間（Ａ）と同様に、心電同期のもとでパルスＸ線は不定周期で繰り返し発生される。区間（Ｃ）でも、システム制御部８は心電図（ＥＣＧ）のＲ波を検出し、ＲＲ間隔から心臓拡張末期を同定し、その時点でのみパルスＸ線を発生させる。区間（Ｃ）のパルスレートは１回／心拍であり、パルスＸ線は不定周期で繰り返し発生される。画像のデータは、造影後の区間（Ｃ）を区間（Ａ）、（Ｂ）と区別するコードとともに記憶される。
【００２８】
なお、上述の説明では、区間（Ａ）、（Ｃ）は、心電同期による不定周期でパルスＸ線を発生しているが、図３に示すように、一定周期でパルスＸ線を繰り返し発生するものであってもよい。典型的には、区間（Ａ）、（Ｃ）のパルスＸ線は、区間（Ｂ）の周期（１／３０秒）より長い一定周期（１秒）で繰り返し発生される。この場合、区間（Ｂ）のパルスレートは、１回／秒となる。
【００２９】
上記パルスレート制御によれば、全期間にわたって高いパルスレートを継続するよりも被曝を低減することができる。しかも、関心部位に造影剤が流通する期間は、比較的高い時間分解能を確保することができる。それ以外の期間は、比較的低い時間分解能であっても、診断に要する心臓拡張末期の画像を収集することができる。
【００３０】
区間（Ｃ）の終了は、画像演算記憶部１０の画像処理結果を受けてシステム制御部８により判定される。例えば、区間（Ｃ）に収集される画像各々から、区間（Ａ）で収集した複数の画像から発生された基準画像（コントロール画像）を差分し、得られた差分画像の全体又は局所領域の濃度合計又は平均濃度（測定値）をリアルタイムで測定する。上記コントロール画像は、区間（Ａ）で収集された複数の画像の平均画像として画像演算記憶部１０により発生され、記憶される。平均処理によりランダムノイズが抑制され得る。
【００３１】
差分画像の測定値は、図４に例示するように、初期値Ｌ0から、撮影部位又は関心部位への造影剤流入に伴って時間経過とともに変化し（図４では低くなり）、造影剤流出に伴って初期値Ｌ0に戻っていく。そこで濃度の初期値Ｌ0に近似するしきい値を、例えば、０．９＜ａ＜１として、ａ・Ｌ0の値に予め設定しておき、測定値が当該しきい値に達した時刻ｔ４で、システム制御部８は高電圧制御部２０にＸ線照射中止信号を出力する。高電圧制御部２０は、システム制御部８からＸ線照射中止信号を受け、高電圧発生部１からＸ線管球２１への管電圧の印加を終了する。なお、区間（Ｃ）が例えば３０秒を経過した場合には、システム制御部８は、測定値がしきい値以下にならなくとも、Ｘ線照射中止信号を出力して、Ｘ線照射を強制的に終了させる。
【００３２】
区間（Ｃ）の終了は、図４に破線で示す濃度時間曲線の傾斜によって判定しても良い。画像演算記憶部１０で繰り返し計算した濃度時間曲線の傾斜が、水平に近似したしきい値に達した時刻ｔ３に区間（Ｃ）は終了する。
【００３３】
なお、画像演算記憶部１０では、様々な画像処理機能を有する。例えば、区間（Ｂ）で得られた画像から心筋灌流に必要な画像を抜き出す。これには区間（Ｂ）で収集した画像から心臓拡張末期の時刻に最も近い画像を各心拍ごとに特定する。区間（Ｂ）では３０ｆｐｓで画像収集されているため、最大でも時間ずれ誤差は３３ｍｓｅｃ以内となる。これは心筋灌流計測には十分に問題のないずれ量である。区間（Ｂ）で特定された複数の画像各々と、区間（Ｃ）で得られた複数の画像から選択された画像とを合成することにより、合成画像を生成する。複数の合成画像から、区間（Ａ）のコントロール画像を差分して、複数の差分画像を生成する。複数の差分画像を区間（Ｂ）内での時間経過に従って配列して表示することにより、造影剤の流出の様子を把握することができる。
【００３４】
なお、上記説明では、造影剤の注入開始／終了により、自動的に区間の切り替え、つまりパルスレート（フレームレート）を変化させるようにしているが、オペレータによる操作部９のパルスレート切り替え指示により、任意の時点でパルスレートを手動で切り替えるようにしてもよい。
【００３５】
上記説明では、３区間としたが、最初の区間（Ａ）をなくして２区間でもよい。この場合コントロール画像は区間（Ｂ）の最初のほうのフレームを用いる。
【００３６】
差分画像を得た後、定量的な評価をする場合がある。この場合、画像上に関心領域（ＲＯＩ）を設定し、その値を読む。
【００３７】
尚、厳密さを要する場合は、血管関心領域と心筋関心領域の間の時間差を考慮するのがよい。すなわち血流は太い血管を流れた後に小血管を得て心筋に達するため、時間差がうまれる。実際には３００ｍｓｅｃ〜１ｓｅｃ程度の時間差がある。病変を有する患者ではこの時間差がさらに長くなる。そこで、太い血管に血管関心領域を設定した場合は、結果、関心領域のデータは、時刻０から（Ｔ−τ）までを積算し、これに対する心筋関心領域は、時刻Ｔのデータを取ることが好ましい。こうした厳密なことをしたい場合には、先の画像収集において１心拍に２枚のデータを取ることが好ましく、すなわち１心拍に２回のＸ線パルス照射をする。
【００３８】
上記定量計測において、さらにノイズを低減したい場合は、近傍フレームで平均処理をするとよい。すなわち図５に示すように、心臓動きが少ない瞬間に、複数、例えば２つのパルスＸ線を一体として発生し、それに伴って収集される２画像を平均した画像から計測を行う。これにより１画像だけの場合よりもノイズが低減される。ただし心臓動きの要素が加算されてしまうので、心臓動きがほとんど無い時間で複数パルスを出すことが必須である。あるいは動き補正をするという手段もある。
【００３９】
この心筋灌流の画像および定量数値を得るのは、空間的に粗くても大丈夫である。そこで例えば生画像が１０２４×１０２４マトリクスの画像であっても、心筋灌流画像は２５６×２５６マトリクスでもよい。そこで心拍同期収集区間では画素加算して画像収集することにより、画像を粗くするかわりに、１画素あたりのＸ線吸収量が少なくててもよいようにする。例えば４画素加算する場合は、Ｘ線量を１／４にする。
【００４０】
心筋灌流を計算するためには、フレーム間でＸ線条件が変わってはならない。Ｘ線条件が変わってしまうと画像の濃度も変わってしまうので、何を測っているのか意味がなくなる。よって提案する収集方法のうち少なくとも１心拍に１枚の画像は、ｋＶおよびｍＡ、ｍｓｅｃを同一にして撮像する。
【００４１】
ここでさらなる工夫としては、ｋＶのみは固定とし、ｍＡｓは変動させてもよい。この場合はフレームごとのｍＡｓ値を記録しておき、心筋灌流画像を作成する場合は記録したｍＡｓ値に基づいて画像を補正してあげる。これが成り立つ理由は単純なＸ線物理式を参照すればよく、ｍＡｓの変動は単純に濃度の変動に比例する。
【００４２】
その他、以下の通り、様々な変形が可能である。
最終区間（Ｃ）の終了判定のために画像濃度値をしきい値と比較したが、所望の灌流計算値をリアルタイムに計算し、計算に十分なデータが得られたら終了判定をするものであってもよい。その他に、区間（Ｃ）において十分なデータが得られたら撮像終了する、近似曲線が描けるのに十分なデータが得られたら撮像終了するものとしてもよい。そのとき、区間（Ｃ）において十分なデータが得られたら撮像を終了してよい旨を通知し、オペレータの終了指示を待機するようにしてもい。当該通知は、メッセージ表示、ピピピ等の音を鳴らすことにより行われる。
【００４３】
画像演算記憶部１０で撮像期間中にリアルタイムで灌流計測値を計算し、表示部１１に表示する。
一定周期でパルスＸ線を発生する区間と、心電同期でパルスＸ線を発生する区間とをオペレータによる手動で切り替えるようにしてもよい。例えば、操作部９に設けられた撮像ボタンが２段階になっていて押し方で切り替わる。サブボタンを有し、メインボタンだけでは一定周期でパルスＸ線を発生するが、メインとサブを同時に押している間は心拍同期でパルスＸ線を発生する。サブボタンは、メインボタンの近傍にあってもよいし、離れていても良い（フットスイッチと手元スイッチ）。最終区間においてはオペレータが撮像ボタンを離したら終わりとする。
【００４４】
心拍同期においては、１心拍につき１照射（１パルス）とするとき、１照射は、心臓収縮末期、心臓拡張末期、心臓拡張中期のいずれかがオペレータにより選択される。
【００４５】
１心拍につき２照射（２パルス）のとき、例えば、第１照射はＲ波から所定後、第２照射はＲ波から異なる所定後とする。所定後とは、Ｒ波からの絶対時間で管理してもよいし、ＲＲ間隔を元にした相対時間で管理しても良い。例えば、第１照射は、心臓収縮末期、拡張末期、拡張中期、のいずれかとする。第２照射は、第１照射より所定時刻だけ前のタイミングとする。所定時刻とは計測点Ａから計測点Ｂまでに造影剤が流れ込むのに要する時間とする
１心拍につき３照射（３パルス）以上とし、３パルスは隣接していて、得られる複数画像は平均化のために利用される。複数心拍につき１照射（１パルス）とする。区間（Ｂ）の心拍同期撮像の時間が所定の時間を超えた場合には、自動的にＫ心拍につき１照射（１パルス）に移行する。
【００４６】
心拍同期の照射においては、Ｘ線エネルギーを同一にする。全てのパルスＸ線で同じｋＶとｍＡｓを用いることを基本とするが、ｋＶを同じとするが、ｍＡｓは異なるものとしてもよい。そのときｍＡｓ値を記録し、後段処理で補正する。
【００４７】
再生時には、見やすいように画像処理して表示する。パルスレートが変化する制御のもとで収集された複数の画像のデータから、一定間隔の画像のデータを生成することを考慮し、区間（Ａ、Ｃ）、（Ｂ）ともに、共通のクロックに対して、パルスＸ線の発生タイミングを同期させる。これによりデータ生成時は必ず同じ心位相のデータが生成できる。
【００４８】
区間（Ａ、Ｃ）においても、区間（Ｂ）と同じ周期で画像を収集してもよい。心電同期による画像列を生成する際には、最も近い時刻の画像を特定して生成する。一定周期の画像群と、心拍間隔の画像群、とを別々に再生する。再生する時は、実時間で再生表示する。
【００４９】
（第２実施形態）
以下図面を参照して第２実施形態を説明する。まず、以下説明で用いる用語について次の通り定義する。
心筋灌流；心筋濃染、心筋パーフュージョン（Perfusion）、心筋ブラッシュ(Blush)、これらは技術的には厳密には異なるが、血液が毛細血管から心筋へ流出入する現象を表現する。
微小灌流(micro circulation)；毛細血管における血流の流れ。
Ｘ線アンジオ装置；Ｘ線診断装置の中の特に血管撮像を主用途とするものである。
Ｘ線画像；被検体を透過したＸ線の強度分布を表す画像であり、Ｘ線画像ともいう。
Ｘ線動画像；２次元検出器によって時間tにわたって繰り返し撮影された一連のＸ線画像のデータセットである。
冠状動脈撮像；造影剤により強調した冠状動脈内腔に関するＸ線画像を撮像することをいう。
冠状動脈；冠状動脈(coronary)
ピクセル；収集されたＸ線画像における画素
フュージョン画像；Ｘ線画像とＸ線画像とをオーバーレイにより合成した画像をいう。
自然対数；Ln
心位相；心電図信号によって検出されたＲ波を目印とし、Ｒ波から次のＲ波までの時間を１００％で規格化し、現在のフレームの時刻を％で表現する。例えば心臓収縮末期の心位相は２５%付近である。
【００５０】
ＲＯＩ；関心領域(Region Of Interest)
カテーテル室；カテーテル検査室、血管造影検査室、血管インターベンション治療室とほぼ同義である。
【００５１】
図１０に本実施形態に係るＸ線診断装置、ここではＸ線アンジオ撮影装置の外観を示し、図１１に機能ブロック図を示す。Ｘ線撮影装置は、ガントリ１００を有する。ガントリ１００は、Ｃアーム７を有する。Ｃアーム７は、機構制御部６に回転自在に支持される。Ｃアーム７の一端にはＸ線発生部２が取り付けられる。Ｘ線発生部２は、Ｘ線管２０とＸ線絞り器２１を有する。高電圧発生部１は、Ｘ線管２０の電極間に印加する高電圧（管電圧）を発生し、またＸ線管２０のフィラメントに供給するフィラメント電流を発生する。高電圧制御部１７は、システム制御部８の制御に従って、高電圧発生部１で発生する管電圧及び／又はフィラメント電流を制御する。
【００５２】
Ｃアーム７の他端にはＸ線検出部５が取り付けられる。Ｘ線検出部５は、検出器１８と、検出器１８の出力を処理して画像データを生成する画像データ生成部１９とを有する。検出器１８は寝台４に載置される被検体３を挟んで、Ｘ線発生部２のＸ線管２０に対峙する。検出器１８は、典型的には、入射Ｘ線を直接的又は間接的に電荷に変換する複数の検出素子（画素）が２次元状に配列されてなる固体平面検出器である。Ｘ線検出部５は、システム制御部８の制御により電荷蓄積、電荷読出及びリセットからなる１サイクルの検出動作を一定周期で繰り返す。被検体３には心電図モニタ端子１５が装着される。心電図モニタ受信部１６は心電図モニタ端子１５からの信号を受信し、被検体３の心電図のデータを発生する。システム制御部８には操作部９が接続される。操作部９には、ハンドスイッチ１２、ディスプレイ及びタッチパネル等を有するユーザインタフェース１４が設けられる。
【００５３】
演算処理部２３は、画像収集保存部１０と画像処理部１１とを有する。画像収集保存部１０は、Ｘ線検出部５から出力される画像データを心位相データを関連付けて記憶する。画像処理部１１は、心臓冠状動脈造影撮影シーケンスのもとで発生された複数の画像のデータに基づいて、心筋への血液供給領域としての冠状動脈上に設定された基準領域に関する時間濃度曲線（基準時間濃度曲線）に対する、心筋領域内に設定された複数の局所領域に関する複数の時間濃度曲線の相関をそれぞれ表す複数の指標を演算するとともに、指標のマップを作成し、マップを冠状動脈像と合成して合成画像を発生する。なお、時間濃度曲線は、元画像の濃度を対象としても良いし、造影前のマスク画像と造影後のコントラスト画像との差分画像上の濃度（造影剤濃度の近似値）を対象としても良い。ここでは、差分画像上の濃度に関する時間濃度曲線として説明する。
【００５４】
本実施形態では、指標として３種類を提供する。
第１の指標Ｋ１は、造影剤注入開始から造影剤注入終了までの造影剤流入期間を対象として、冠状動脈の基準時間濃度曲線を入力関数、心筋局所領域の時間濃度曲線を出力関数とした時の、局所心筋に関する「血液の流入状態を表す指標Ｋ１」として計算される。第２の指標Ｋ２は、造影剤注入終了から撮影完了までの造影剤流出期間を対象として、冠状動脈の基準時間濃度曲線を入力関数、心筋局所領域の時間濃度曲線を出力関数とした時の、局所心筋に関する「血液の流出状態を表す指標Ｋ２」として計算される。第３の指標Ｋ３は、指標Ｋ１と指標Ｋ２とから計算される。尚、この時の時間濃度曲線とは、人体に投与された造影剤が注目部位（本実施形態では冠状動脈と心筋）に分布したり通過したりする際の造影剤量のみに比例した信号成分と定義する、。
【００５５】
図１２、図１３は一般的なＸ線アンジオ撮影で得えられる冠状動脈造影画像（以下ＣＡＧ画像という）を例示している。ＣＡＧ画像では、造影剤のＸ線の吸収線量が大きいので冠状動脈を通過中は冠状動脈の形状が他の組織と区別できる程にコントラストがついて識別できる（図１２参照）。図１３に示すように、カテーテル１０１が冠状動脈１０２まで挿入操作され、その位置でカテーテル１０１から一定時間継続的に造影剤が注入される。Ｘ線撮像は、少なくとも造影剤注入開始前から、造影剤注入終了後所定時間経過までの期間にわたって行われる。
【００５６】
造影剤が冠状動脈にカテーテル１０１を経由して注入されるときに図１０のＸ線診断装置にてＸ線画像を心電図とともに収集する。画像収集後、操作部９を介してオペレータによりＣＡＧ画像上に冠状動脈上の基準領域（心筋血液供給領域）１０３と、心筋上に複数の心筋局所領域１０４とが設定される（図３２Ａ、図３３参照）。
【００５７】
複数の心筋局所領域１０４は、心筋領域上に設定され、典型的には複数の画素を有する。心筋局所領域１０４の造影剤量に相当する濃度は典型的には画素平均値として計算される。しかし、心筋局所領域１０４は単一の画素を有するものであってもよい。心筋血液供給領域１０３は、典型的には、血管と略等価な幅又は少し小さい幅を有する矩形形状を有し、血管に沿って任意の向きに設定され、複数の画素が含まれる。心筋血液供給領域１０３の造影剤量に相当する濃度としては典型的には画素平均値として計算される。心筋血液供給領域１０３は、インジェクターと心筋注目領域の間の流路の任意の部分、より限定するとカテーテルの任意の部分、又はカテーテルの出口（冠状動脈入口に等価）から心筋注目領域の間のの任意の部分に設定される。
【００５８】
画像処理部１１は、心筋血液供給領域１０３に関する時間濃度曲線（Time Density curve；ＴＤＣ）を生成する。同様に、画像処理部１１は、複数の心筋局所給領域１０４に関する複数の時間濃度曲線を生成する。この生成過程の詳細を図３４に示す。
【００５９】
図３２のＣＡＧＸ線画像の座標(x,y)での時刻tでの画素値をCAG(x,y,t)とし、図３４の画素値をこの自然対数をとったｌｎ（CAG(x,y,t)）とする。これにより組織或いは造影剤の単位体積あたりの減弱係数をμ、図２１のＸ線の透過方向のその組織または造成剤の分布厚さをＬとすると、図１４の画素はｌｎ（CAG(x,y,t)）∝μＬの関係がある。２次元画像である図１３画像上で認識される冠状動脈や心筋部画像に分布する造影剤の総量はμＬに比例する。すなわち造影剤の無い時刻t=t0の画像Ｔ0の画素値をCAG(x,y,t0)とすると、
ｌｎ（CAG(x,y,t)）−ｌｎ（CAG(x,y,t0)）
は、図１３のＸ線画像上で認識される冠状動脈や心筋部画像に分布する造影剤の総量に比例して図１４に示す画像Ｔ１からＴ５で示す一連の画像内では相対的に冠状動脈や心筋部画像などに分布する造影剤総量を示している。
【００６０】
すなわち、ＣＡＧ画像は一方向からの透過画像であるので、画像中心付近、通常は画像中心に位置させる心筋の検査対象部分では心筋に対して略垂直にＸ線が当たるため厚み効果はあまり生じないが、画像中心から外れた周縁部分では心筋に対してＸ線がその接線に近い向きで当たるため、厚み効果が強く表れる。つまり、中心部分に対して周縁部分は造影剤濃度の積分距離が長くなる。従って、心臓の周縁部分で測定される造影剤濃度の信頼性は、中心部分で測定される造影剤濃度の信頼性よりも低い。しかし、上述したように、通常は画像中心付近に心筋の検査対象部分を位置合わせするためその部分に関しては高い信頼性で指標を得ることができる。なお、検査対象ではない画像中心から所定距離離れた信頼性の低い心臓の周縁部分に関する指標は表示対象から外してもよい。
【００６１】
なお、図１４の血液供給領域１０３は、イラスト上でわかりやすくするために円で描画したがその形状は図１３で示すようにカテーテル１０１を通過して冠状動脈に流入して動脈血と混合した時のＴＤＣを計測できるように冠状動脈を囲むように例えば矩形に近い形状に設定される。心筋上の局所領域１０４はカテーテル１０１にて造影剤が注入される冠状動脈の支配領域（その冠状動脈を通過した血液が供給されることにより活動している心筋の領域）であって狭窄やその先の毛細血管の微小循環の有無が疑われる心筋部位全体について、例えば関心領域１０４で示すＮ×Ｎピクセル（Ｎ＝１，３，５，．．．奇数）で例えばＮ＝５として画素をサンプリングし、例えば５×５サイズのスムージング処理後以下で述べる演算結果をその画素の中心位置を代表位置として記録する。
【００６２】
［観測］
上記のように、観測される値はｌｎ（CAG(x,y,t)）−ｌｎ（CAG(x,y,t0)）である。ここで心筋血液供給領域１０３に関する血管単位体積あたりの造影剤相対濃度の基準ＴＤＣをCa(t)とし、同様に心筋単位体積あたりの造影剤相対濃度をCmyo(t)とする。すると観測値とCa、Cmyoとの関係は次式で表される。なお、Caは造影剤注入前後の差分値を表しているので、初期値Ca(0)は０値になる。
【数１】
【００６３】
ここで(xa,ya)は図１３の冠状動脈の上流に設定した関心領域１０３内の画素であり、値の平均値、或いは関心領域内の総量を表す。どちらを用いても良いが以下の操作では選択した一方のみを使用する。混在使用はしない。ｔは時刻を表し、単位は秒である。ａは動脈を意味する。
【数２】
【００６４】
ここで(xmyo,ymyo)は図１３に示す心筋に設定した局所の関心領域１０４内の画素であり、値の平均値、或いは関心領域１０４内の総量を表す。複数の関心領域１０４は、図１５Ａ、図１５Ｂに示す操作部９を介してＣＡＧ画像上に設定されたパフュージョンの演算範囲１０５を図３２Ａ、図３３に示すように分割した相互に同一形状及びサイズの複数の局所領域である。なお、添え字のmyoはmyocardium（冠状動脈）を意味する。
【００６５】
［モデル］
次に、心筋への造影剤の流出入モデルとして後述の（理論１）を導入する。これにより次の関係を得る。
【数３】
【００６６】
［観測値のモデルへの適応］
画像Ｔ１からＴ５は一連の心電同期で撮像したＣＡＧ画像のうちの特定の心位相１０６に対応するＣＡＧ画像である。その収集時刻をｔ1,ｔ2,..ｔ5とする。かかる心電同期画像Ｔ１〜Ｔ５のセットについて、画像処理部１１では、［観測］と［モデル］で導かれた（１）（２）（３）式にて記述する処理を行う。
【数４】
【００６７】
画像処理部１１では、心筋血流に比例する実際的な心筋局所領域への血液の流入の様子を反映した指標Ｋ１'を算出する。まず、造影剤がカテーテルから注入開始された時点から注入完了の時点までの造影剤流入期間が、指標Ｋ１の計算対象期間として設定される。図１６に示すように、当該造影剤流入期間内に収集した複数のＸ線画像から、心筋への血液供給領域１０３に関する時間濃度曲線Ｃａ（ｔ）と、複数の心筋局所領域１０３に関する複数の時間濃度曲線Ｃmyo（ｔ）とが生成される。
【００６８】
次に、画像処理部１１では、図１７に示すように、血管供給領域１０３を通過する造影剤の時間積分（造影剤の流入量）∫Ｃa(t)を横軸とし、縦軸に、心筋局所領域１０４の造影剤存在量Ｃmyo（ｔ）として、時刻ごとの値をプロットする。つまり、心筋全体への血液供給量に対する心筋局所における血液の取り込み量の時間変化の離散分布が生成される。この離散分布に対して直線フィティング処理により得られた直線の傾きとして指標（第１指標）Ｋ１'が算出される。この指標Ｋ１'は、心筋全体への血液供給に対して心筋局所がどの程度追従して血液を受け入れているか否かを定量化したものとしての意義を有する。指標Ｋ１'が正常範囲から低値側に外れているとき、冠状動脈の血液流入に対して、心筋局所での血液流入が追従していない、つまり当該心筋局所領域は流入障害を起こしている可能性がある。指標Ｋ１'が正常範囲から高値側に外れているとき、当該心筋局所領域は流出障害を起こしている可能性がある。画像処理部１１では、すべての局所領域１０４について同様の処理により指標Ｋ１'を計算し、指標Ｋ１'のマップを作成する。この得られたマップをＣＡＧ画像に重ねあわせ（フュージョン）て合成画像を生成し、表示部２４に表示することができる。これにより心筋機能を心筋組織との位置関係とともに判断することができる。
【００６９】
変形例Ａ）
なお、上述では一連の心電図同期ＣＡＧ画像Ｔ１からＴ５は、連続した画像から抽出したが、患者への被曝を低減するために、撮像時に特定の心電図位相のみＸ線を発生させるように図１０の装置を制御して後述の（理論１）に必要な画像セットのみを収集することもできる。
変形例Ｂ）
さらに後述の（理論２）によれば造影剤が心筋部に到達してさらに排泄される流出期間における一連のＴＤＣカーブ（図１８）から図１９で示すグラフィプロット法で直線近似により得られたＫ１とＫ２により血流と造影剤の心筋からのバックディフュージョン（心筋から血管に造影剤が戻る量）を計測でき、図１５を用いて前述同様に表示することができる。同様に時間ｔｎは患者への被曝を低減するために、撮像時に特定の心電図位相のみＸ線を発生させるように図１０の装置を制御して理論１に必要な画像セットのみを収集することもできる。
【００７０】
変形例Ｃ）
さらにＫ1とＫ２の比（Ｋ１／Ｋ２）を計算して表示することにより、「C. Michael Gibson, MS, MD; Christopher P. Cannon, MD; Sabina A. Murphy, MPH; Kathryn A. Ryan, BS; Rebecca Mesley, BS; Susan J. Marble, RN, MS; Carolyn H. McCabe, BS; Frans Van de Werf, MD, PhD; Eugene Braunwald, MD; for the TIMI (Thrombolysis In Myocardial Infarction) Study Group, Relationship of TIMI Myocardial Perfusion Grade to Mortality After Administration of Thrombolytic Drugs, Circulation. 2000;101:125-130.」に示された分類法を定量化（自動化）することが可能になる。すなわち図２８の分類表に従って自動分類したり、図２９Ａ、図２９Ｂのような指標Ｋ１、Ｋ２の関係を表すグラフにプロットして表現することが可能になる。
【００７１】
変形例Ｄ）
本実施形態により、同一患者、同一方向の一連の心電図同期ＣＡＧ画像を二種類収集する。すなわち、図２０Ａ、図２０Ｂに示すように、安静時の心筋血流画像２０５（＝Ｐrest(x,y)）と例えばアデノシンのような心筋血流を増加させる効果を有する薬剤を投与後の薬剤負荷時の心筋血流画像２０６（＝Ｐstress(x,y)）を収集する。安静時と負荷時の比較画像として（心筋血流画像２０６）／（心筋血流画像２０５）＝Ｐstress(x,y)／Ｐrest(x,y)を定義すると、心筋血流画像２０５，２０６で同一方向の収集された同一位置(x,y)では撮像方向の心筋分厚さは同じなので割り算により近似的にキャンセルされ、単位心筋体積あたりの血流増加比（心筋Flow reserve）を得ることができる。
【数５】
【００７２】
また同様に、狭窄等の治療前の指標Ｋ'beforeと治療後の指標Ｋ'afterを計測し、割り算すると、単位心筋体積あたりの血流増加比（治療による回復比）を得ることができる。
【００７３】
変形例Ｅ）
また実施形態によれば、図２７に例示するように、後述の［理論１、２］を完全に演算しなくても、Ｃmyo（ｔ）＝ｌｎ（CAG(x,y,t)−ｌｎ（CAG(x,y,t0)）とおいたときに、Ｃmyo（ｔ）−Ｃmyo（ｔ０）の画像は心筋への血流を示しており、この差分動画像を表示することも同一検査時に限って可能である。ここで、ｔ０＝造影剤注入直前である。
【００７４】
（理論１）
時刻ｔにおける心筋部の造影剤相対濃度Cmyo(t)、冠状動脈血液中の造影剤相対濃度時間曲線をCa(t)、心筋流入血流をＫＫ１、心筋流出血流をＫＫ２と表す。図２３において、心筋部と冠状動脈血中の造影剤の質量バランスは式（１）で表現できる
【数６】
【００７５】
心筋部の血流あるいは血流と相関のあるパラメータを以下に述べる造影剤が心筋部へ流入する現象を利用するPatlak plot法で計算することができる。造影剤が心筋部へ流入し始める条件下（0≦t≦MTT）（MTT＝mean tangent time＝造影剤が動脈から心筋に入り静脈に出てくるまでの平均通過時間、心筋では通常５〜１０ｓｅｃ程度）では、流出量ｋ２Ｃｍｙｏ(t)は極めて少なく、Ｃｍｙｏ(t)≒０、Ｋ１Ｃa(t)＞＞Ｋ２Ｃｍｙｏ(t)の条件が成立し、このときEq.(1) は以下のように簡略化できる。
【数７】
【００７６】
Ｋ１は造影剤が血液から心筋間質部へ移行する移行常数（sec-1）或いは (ml/min/g)を表し、血流に比例する。
Rutland MD. A single injection technique for subtraction of blood background in 131I-hippuran renograms. Br J Radiol 1979;52:134-137.
Patlak CS, Blasberg RG, Fenstermacher JD. Graphical evaluation of blood-to-brain transfer constants from multiple-time uptake data. J Cereb Blood Flow Metab 1983;3:1-7.
（理論２）
時刻ｔにおける心筋部の造影剤相対濃度をCmyo(t)、冠状動脈血液中の造影剤相対濃度をCa(t)と表す。図１０において、心筋部と冠状動脈血中の造影剤の質量バランスは式（１）で表現できる
【数８】
【００７７】
心筋部の血流あるいは血流と相関のあるパラメータを以下に述べる造影剤が心筋から流出するクリヤランスを利用する方法で計算することができる。
【００７８】
Eq. (1)を積分し(2)を得る。
【数９】
【００７９】
Ｙ軸にY(t) Ｘ軸にX(t) をプロットすると直線の傾きがＫ２、Ｙ軸の切片はＫ１を表す。Ｋ１は造影剤が血液から心筋間質部へ移行する移行常数で（sec-1）或いは (ml/min/g)を表し血流に比例する。Ｋ２は造影剤が心筋間質部から血液へ移行する移行常数（sec-1）を表す（以下文献を図３１とともに参照）。
【００８０】
Yokoi T, Iida H, Itoh H, Kanno I. A new graphic plot analysis for cerebral blood flow and partition coefficient with iodine-123-iodoamphetamine and dynamic SPECT validation studies using oxygen-15-water and PET. J Nucl Med 1993;34:498-505.
（実際の処理フロー）
以下、実際の画像処理フローの一例を説明する。ここに示す例は一例であり、前述の理論に基づく他の処理方法であってもよい。
図３２を用いて、本実施形態の全体の具体的処理フローを説明する。Ｘ線画像は血管造影された被写体を写した動画像であることを特徴とする（Ｓ１１）。なお、例えば動画は３〜６０秒の継続時間を有する。１秒あたり１０〜３０枚の画像を有する動画像である。画像は同一方向、同一Ｘ線条件で撮像されたものである。撮像中には、患者の動き、寝台の動き、は無いものとする。
【００８１】
画像処理部１１では１心拍ごとに特定心位相に対応する１枚の画像を抽出し（Ｓ１４）、その画像を造影剤注入前の同心位相の画像から差分処理する（Ｓ１２）。具体的には、造影剤注入前の同心位相が複数枚ある場合は、平均した画像を差分処理する。
【００８２】
画像処理部１１では動き補償処理する（Ｓ１３）。画像パターン（血管、心壁、カテ）を基準とした、パターンマッチング処理により動き補償量を検出する。特徴的なパターンの無いフレームでは、他時刻同位相で得られた移動量で代用する。動き補償は、画像全体を動かす、もしくはＲＯＩ移動位置がわかるテーブルを保有する。動き補償では、患者の多少の体動き、呼吸運動、も補償される。
【００８３】
特定の心位相は、操作部９を介して操作者により予め指定される。例えば心臓拡張末期に指定される。パルス撮像された多フレームから、最も近い時刻のフレームを抜き出してくる。なお、血液供給領域１０３の濃度時間曲線Caと心筋局所領域１０４の濃度時間曲線Cmyoでは異なる心位相を用いてもよい。Ca(t)に対しては、Cmyo(t+T)を用いる。ここでＴは、血液供給領域１０３から心筋局所領域１０４まで造影剤が流れていくのに要する遅延値である。具体的には１〜１０フレーム（１／３０〜１／３ｓｅｃ）程度であることが多い。この遅延時間は、一般にTIMI frame count(TFC)もしくはcorrected TIMI frame count(CTFC)と呼ばれる。
【００８４】
画像処理部１１では、濃度時間積分Σca(t)に対するCmyo(t)の相関を表す傾斜としての指標Ｋ１を計算する（Ｓ１５）。図１７に例示したようにフィッティングして回帰直線を計算する。回帰は所定時間区間のみで計算させ、この所定区間は造影剤インジェクターと連動して決定する。すなわちインジェクターからの造影剤注入開始時刻をJ1、注入終了時刻をJ2とすると、所定区間の開始時刻t1は例えばt1=J1+1sec、所定区間の終了は例えばJ2+2secとする。また、所定区間は画像の値から判定するものであってもよい。すなわちＴＤＣが立ちあがる時刻をt1とし、ＴＤＣがフラットになる時刻を所定区間の終了と定める。回帰直線との相関値を計算する。不整脈を含む心拍は計算から除外することが好ましい。ひとつの動画像に対してひとつのＫ１を計算する。ただし高速化を目的として、フレームごとに随時Ｋ１を計算する手法でもよい。上述したように、画素ごとではなく、周囲N×M画素の平均値（Medianでもよい）をその画素の値とする。表示部２４は血流画像を表示する（Ｓ１７）。指標Ｋ１のマップを表示する。Ｋ１のマップは濃度値をカラー値に変換して表示してもよい。収集した元の動画像を表示し、その直後にＫ１のマップを重ねて表示しても良いし、元画像とＫ１のマップを左右に並べて表示してもよい。
【００８５】
上述した変形例Ａに対応して、本実施形態は、画像処理を、Ｘ線診断装置に一体化してもよいし（図１０、図２２）、もしくはＸ線診断装置から独立させた画像処理装置（図３０、図２３）として提供しても良い。血流画像（指標画像、マップ画像）に、元画像（血管画像、撮像された画像）を重ねて表示する（Ｓ１６）。指標をそのままカラーテーブルにあてはめて得られるカラー血流画像に、白黒の元画像（血管画像）を重ねて表示する。なお、元画像は動画であり、血管動画像に血流カラー静止画像を重ねるのもよい。実際的には、血管静止画像に血流カラー静止画像を重ねる。動画から静止画像を作成するには、動画のうち最も濃い１フレームを選択するか、各画素ごろに動画の画素値の最も小さい値を検出して画像とする。元画像とは、差分前の画像、もしくは差分後の画像とする。グラフを表示することを特長とする（Ｓ２１）。特定のＲＯＩ（x,y）について、図１６〜１９のグラフを表示する。図１７、図１８のグラフには、フィッティング関数（回帰直線）の傾きと切片と相関値を表示する。血液供給領域１０３はオペレーターが指定する。血液供給領域１０３を画像上に指定のためのグラフィカルユーザインタフェース１４を備える。１心拍の特定の心位相で１パルスＸ線を照射する（図２３）。
【００８６】
画像処理部１１は指標Ｋ１とＫ２を計算する（図２４）。指標Ｋ１に対するＫ２の比（Ｋ２）／Ｋ１を計算する（図２５Ａ）。画像処理部１１はＫ１とＫ２の比のマップを生成する。
【００８７】
また、画像処理部１１は、図２５Ｂに示すように、Ｋ１を横軸に、Ｋ２を縦軸に配置して、局所領域毎にＫ１とＫ２に対応する点をプロットする。分布上では、「正常」、「流入障害」、「流出障害」、「流出入障害」のような心筋機能を表す４種の区分に、各心筋局所領域を分類され得る。つまり、心筋局所領域を指標Ｋ１を第１閾値ＴＨ１に対する高低の区別と、指標Ｋ２の第２閾値ＴＨ２に対する高低の区別との組み合わせに分類する。これにより心筋機能を局所ごとに判定することができる。
【００８８】
さらに、画像処理部１１では、図３２Ｂに示すように４つの区分に４つのコード（１〜４）を割り当てて、コードを図３２Ａに示すように心筋局所領域に当てはめて元画像に重ねて表示するようにしてもよいし、図３３に示すように、４つの区分を異なる表示態様、例えば輝度又は色相で区別して元画像に重ねて表示するようにしてもよい。この表示により、心筋機能障害の状態をその区別とともに空間的に把握することができる。
また、区分は、図２８、図２９Ｂに示したように分類してもよい。
【００８９】
さらに画像処理部１１は、比較解析処理を装備する（図２６）。例えば同一手術中の時刻Aと時刻B（C,D,E,,,）の指標Ｋ１をそれぞれ計算し左右に並べて表示する。もしくは、時間的に切り替えて表示する。時刻Aと時刻B（C,D,E,,,）のＫ１の比を計算して、比を画像表示する。比のカラー画像に、元画像を重ねて表示する。例えば治療前後での比のカラー画像を、血管が描出されている元画像に重ねて表示することにより、心筋のどの領域は治療によって血流回復が得られたかが可視化されてわかるようになる。
【００９０】
以上説明した本実施形態の特徴としては以下の通りである。
Ｘ線ＣＴ装置から出されている造影剤を用いたPerfusion計算方法との差異としては、画像処理がＣＴとＸ線では大きく異なり、Ｘ線特有の処理が必要となる。Ｘ線は動画像であるため、特有の工夫、具体的には、１心拍に１枚の画像の抽出、背景差分、動き補償、厚み方向の処理、血管ROIの設定等の特有の処理を必要とした。
【００９１】
Ｘ線での時間濃度曲線の傾きを見る公知方法との差異としては、[理論１、２]で述べた考え方そのものが異なる。指標Ｋ１においては横軸を流入量としたグラフを描く考え方が従来にはない。これにより、後段の解析や、結果画像が異なってくる。処理方法は時間濃度曲線の傾きを得る処理とはまったく異なる。
【００９２】
文献の（或いは図２８の）MBG（TMP）手法は、画像数値じたいには意味を持っていなかった。本実施形態では画像数値は心筋への流れ込み流速であるという意味を持つ点で、MBGより明らかに有用性が高い。流入量でNormalizeしているので、注入する造影剤量に左右されず、臨床の場でも定量性が得られる。
【００９３】
目視ではないので定量性がある。予後予測できるので、この値を参考にしてインターベンション治療を終わって良いかどうかの判断ができる。処置（治療）前後で比較できるので、どれだけ流量が増えたか、すなわち心筋への血液供給がどれだけ改善したかが数値でわかる。処理自動化することにより、臨床でも簡単に使えるようになる。
【００９４】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【図面の簡単な説明】
【００９５】
【図１】本発明の第１実施形態に係るＸ線診断装置の構成を示す図。
【図２】図１のシステム制御部による撮像動作の一例を示す図。
【図３】図１のシステム制御部による他の撮像動作を示す図。
【図４】図１のシステム制御部による撮像終了判定処理を補足するための時間濃度曲線を示す図。
【図５】図１のシステム制御部によるさらに他の撮像動作を示す図。
【図６】従来の撮像動作を示す図。
【図７】従来の撮像動作を示す図。
【図８】図１のシステム制御部による他の撮像動作を示す図。
【図９】図１のシステム制御部による他の撮像動作を示す図。
【図１０】本発明の第２実施形態によるＸ線診断装置の構成を示す図。
【図１１】図１０の特定部分の詳細図。
【図１２】冠状動脈造影画像の一例を示す図。
【図１３】本実施形態において、冠状動脈造影画像上に設定されるＲＯＩを示す図。
【図１４】本実施形態において、画像収集シーケンスを示す図。
【図１５Ａ】本実施形態において、パフュージョンの演算範囲を示す図。
【図１５Ｂ】本実施形態において、ＣＡＧ画像上に設定操作されたパフュージョンの演算範囲を示す図。
【図１６】図１１の画像処理部により生成される造影剤流入期間における時間濃度曲線を示す図。
【図１７】図１１の画像処理部により計算される指標Ｋ１を示す図。
【図１８】図１１の画像処理部により生成される撮像期間における時間濃度曲線を示す図。
【図１９】図１１の画像処理部により計算される指標Ｋ２を示す図。
【図２０Ａ】本実施形態において、安静時と負荷時のパフュージョンの演算範囲を示す図。
【図２０Ｂ】図２の画像処理部により生成される安静画像と負荷画像との比較画像を示す図。
【図２１】本実施形態において、透過画像の厚み効果の補足図。
【図２２】本実施形態の処理手順を示す流れ図。
【図２３Ａ】本実施形態の他の処理手順を示す流れ図。
【図２３Ｂ】本実施形態の他の処理手順を示す流れ図。
【図２４】本実施形態の他の処理手順を示す流れ図。
【図２５Ａ】本実施形態の他の処理手順を示す流れ図。
【図２５Ｂ】本実施形態の２種の指標の分布を示す図。
【図２６】本実施形態の他の処理手順を示す流れ図。
【図２７】本実施形態の他の処理手順を示す流れ図。
【図２８】図１１の画像処理部による２種の指標を用いた局所領域の分類のための分類表を示す図。
【図２９Ａ】図１１の画像処理部により生成される２種の指標の分布を示す図。
【図２９Ｂ】図２８の分類表に対応する指標分布上の区分を示す図。
【図３０】本実施形態において、Ｘ線診断装置から独立させた画像処理装置を示す図。
【図３１】本実施形態において、指標Ｋ１、Ｋ２の説明図。
【図３２Ａ】図１１の画像処理部により２種の指標Ｋ１、Ｋ２から生成した心筋機能マップを示す図。
【図３２Ｂ】図３２Ａのマップコードを示す図。
【図３３】図３２Ａの他の心筋機能マップを示す図。
【符号の説明】
【００９６】
１…高電圧発生部、２…Ｘ線発生部、３…被検体、４…寝台、５…Ｘ線検出部、６…機構制御部、７…Ｃアーム、８…システム制御部、９…操作部、１０…画像演算記憶部、１１…表示部、１２…ハンドスイッチ、１３…スイッチホルダ、１４…ユーザインターフェース、１５…インジェクター、２０…高電圧制御部。
【特許請求の範囲】
【請求項１】
パルスＸ線を発生するＸ線源と、
被検体を透過したＸ線を検出するＸ線検出器と、
前記Ｘ線検出器の出力に基づいて画像を発生する画像発生部と、
前記パルスＸ線を撮像期間内に繰り返し発生させるとともに前記パルスＸ線の発生周期を前記撮像期間内で変化させるために前記Ｘ線源を制御する制御部とを具備することを特徴とするＸ線診断装置。
【請求項２】
前記制御部は、前記撮像期間内の第１区間では前記パルスＸ線を一定周期で発生させ、前記撮像期間内の第２区間では前記パルスＸ線を心電同期により不定周期で発生させることを特徴とする請求項１記載のＸ線診断装置。
【請求項３】
前記制御部は、前記被検体に対する造影剤の注入開始又は注入終了の時点を基準として前記第１区間と前記第２区間とを切り替えることを特徴とする請求項２記載のＸ線診断装置。
【請求項４】
前記制御部は、前記造影剤の注入開始又は注入終了の時点から所定の遅れ時間を経過した時点に前記第１区間と前記第２区間とを切り替えることを特徴とする請求項３記載のＸ線診断装置。
【請求項５】
前記制御部は、前記画像の濃度に基づいて前記第１区間と前記第２区間とを切り替えることを特徴とする請求項２記載のＸ線診断装置。
【請求項６】
前記第１区間は１〜５秒に予め設定され、前記第２区間は１〜６０秒に予め設定されることを特徴とする請求項２記載のＸ線診断装置。
【請求項７】
前記第１区間は任意の時間長に予め設定され、前記第２区間は任意の心拍数に予め設定されることを特徴とする請求項２記載のＸ線診断装置。
【請求項８】
前記第１区間では前記パルスＸ線は１０〜３０回／秒の頻度で発生され、前記第２区間では前記パルスＸ線は１〜５回／心拍の頻度で発生されることを特徴とする請求項２記載のＸ線診断装置。
【請求項９】
前記制御部は、前記撮像期間内の第１区間では前記パルスＸ線を心電同期により不定周期で発生させ、前記撮像期間内の前記第１区間に続く第２区間では前記パルスＸ線を一定の周期で発生させ、前記撮像期間内の前記第２区間に続く第３区間では前記パルスＸ線を心電同期により不定周期で発生させることを特徴とする請求項１記載のＸ線診断装置。
【請求項１０】
前記制御部は、前記被検体に対する造影剤の注入開始の時点を基準として前記第１区間から前記第２区間に切り替え、前記造影剤の注入終了の時点を基準として前記第２区間から前記第３区間に切り替えることを特徴とする請求項９記載のＸ線診断装置。
【請求項１１】
前記制御部は、前記撮像期間内の第１区間では前記パルスＸ線を第１周期で発生させ、前記撮像期間内の第２区間では前記パルスＸ線を前記第１周期と異なる第２周期で発生させることを特徴とする請求項１記載のＸ線診断装置。
【請求項１２】
前記制御部は、前記被検体に対する造影剤の注入開始又は注入終了の時点を基準として前記第１区間と前記第２区間とを切り替えることを特徴とする請求項１１記載のＸ線診断装置。
【請求項１３】
前記制御部は、前記撮像期間内の第１区間では前記パルスＸ線を第１周期で発生させ、前記撮像期間内の前記第１区間に続く第２区間では前記パルスＸ線を前記第１周期より短い第２周期で発生させ、前記撮像期間内の第３区間では前記パルスＸ線を前記第２周期より長い第３周期で発生させることを特徴とする請求項１記載のＸ線診断装置。
【請求項１４】
前記制御部は、前記撮像期間内の第１区間では前記パルスＸ線を心電同期により不定周期で発生させ、前記撮像期間内の前記第１区間に続く第２区間では前記パルスＸ線を一定周期で発生させ、前記撮像期間内の第３区間では前記パルスＸ線を心電同期により不定周期で発生させることを特徴とする請求項１３記載のＸ線診断装置。
【請求項１５】
前記制御部は、前記パルスＸ線の発生周期を、前記撮像期間内の操作者により指示された時点で延長又は短縮することを特徴とする請求項１記載のＸ線診断装置。
【請求項１６】
前記制御部は、前記第１、第２区間にわたって一定周期でＸ線検出動作を繰り返すために前記Ｘ線検出器を制御することを特徴とする請求項１記載のＸ線診断装置。
【請求項１７】
前記画像のデータを前記第１、第２区間を区別して記憶する記憶部をさらに備えることを特徴とする請求項２又は９記載のＸ線診断装置。
【請求項１８】
前記画像を前記第１区間と前記第２区間とで別々に再生する再生部をさらに備えることを特徴とする請求項１７記載のＸ線診断装置。
【請求項１９】
冠状動脈造影撮影シーケンスのもとで発生された複数の画像のデータを記憶する記憶部と、
前記複数の画像のデータに基づいて、心筋への血液供給領域に設定された基準領域に関する基準の時間濃度曲線と、心筋領域内に設定された複数の局所領域に関する複数の時間濃度曲線とを発生し、前記基準の時間濃度曲線に対する前記複数の時間濃度曲線の相関をそれぞれ表す複数の指標を演算する演算部とを具備することを特徴とする画像処理装置。
【請求項２０】
前記演算部は、前記心筋への造影剤流入期間における前記基準の時間濃度曲線に対する前記時間濃度曲線の相関を前記指標として前記局所領域ごとに演算することを特徴とする請求項１９記載の画像処理装置。
【請求項２１】
前記演算部は、前記心筋からの造影剤流出期間における前記基準の時間濃度曲線に対する前記時間濃度曲線の相関を前記指標として前記局所領域ごとに演算することを特徴とする請求項１９記載の画像処理装置。
【請求項２２】
前記演算部は、前記心筋への造影剤流入期間における前記基準の時間濃度曲線に対する前記時間濃度曲線の相関を表す第１指標と、前記心筋からの造影剤流出期間における前記基準の時間濃度曲線に対する前記時間濃度曲線の相関を表す第２指標とを前記局所領域ごとに演算することを特徴とする請求項１９記載の画像処理装置。
【請求項２３】
前記演算部は、前記第１指標と前記第２指標との比を前記局所領域ごとに演算し、前記比のマップを生成することを特徴とする請求項２２記載の画像処理装置。
【請求項２４】
前記比のマップを前記画像に重ねて表示する表示部をさらに備えることを特徴とする請求項２２記載の画像処理装置。
【請求項２５】
前記基準の時間濃度曲線を前記時間濃度曲線とともにグラフとして表示する表示部をさらに備えることを特徴とする請求項１９記載の画像処理装置。
【請求項２６】
前記演算部は、前記局所領域を前記第１指標の第１閾値に対する高低の区別と、前記第２指標の第２閾値に対する高低の区別との組み合わせに分類することを特徴とする請求項２２記載の画像処理装置。
【請求項２７】
前記演算部は、前記分類した組み合わせに関するマップを生成することを特徴とする請求項２６記載の画像処理装置。
【請求項２８】
前記演算部は、治療前に収集した前記心筋への造影剤流入期間における前記基準の時間濃度曲線に対する、治療後に収集した前記心筋への造影剤流入期間における前記基準の時間濃度曲線の相関を前記指標として前記局所領域ごとに演算することを特徴とする請求項１９記載の画像処理装置。
【請求項２９】
前記演算部は、治療前に収集した前記画像のデータから求めた前記指標に対する治療後に収集した前記画像のデータから求めた前記指標の比較結果を演算することを特徴とする請求項１９記載の画像処理装置。
【請求項３０】
前記指標は、前記心筋への造影剤流入期間における前記画像のデータから求められることを特徴とする請求項２９記載の画像処理装置。
【請求項３１】
前記指標は、前記心筋への造影剤流出期間における前記画像のデータから求められることを特徴とする請求項２９記載の画像処理装置。
【請求項３２】
前記演算部は、薬剤投与前に収集した前記画像のデータから求めた前記指標に対する薬剤投与後に収集した前記画像のデータから求めた前記指標の比較結果を演算することを特徴とする請求項１９記載の画像処理装置。
【請求項３３】
前記演算部は、前記指標に関するマップを生成することを特徴とする請求項１９記載のＸ線診断装置。
【請求項３４】
前記指標に関するマップを前記画像にオーバーレイして表示する表示部をさらに備えることを特徴とする請求項３３記載の画像処理装置。
【請求項３５】
前記演算部は、前記第１指標に関するマップ又は前記第２指標に関するマップを生成することを特徴とする請求項３０記載の画像処理装置。
【請求項３６】
Ｘ線を発生するＸ線管と、
被検体を透過したＸ線を検出するＸ線検出器と、
前記Ｘ線検出器の出力に基づいて画像のデータを発生する画像発生部と、
冠状動脈造影撮影シーケンスのもとで発生された複数の画像のデータに基づいて、心筋への血液供給領域に設定された基準領域に関する基準の時間濃度曲線と、心筋領域内に設定された複数の局所領域に関する複数の時間濃度曲線とを発生し、前記基準の時間濃度曲線に対する前記複数の時間濃度曲線の相関をそれぞれ表す複数の指標を演算する演算部とを具備することを特徴とするＸ線診断装置。
【請求項１】
パルスＸ線を発生するＸ線源と、
被検体を透過したＸ線を検出するＸ線検出器と、
前記Ｘ線検出器の出力に基づいて画像を発生する画像発生部と、
前記パルスＸ線を撮像期間内に繰り返し発生させるとともに前記パルスＸ線の発生周期を前記撮像期間内で変化させるために前記Ｘ線源を制御する制御部とを具備することを特徴とするＸ線診断装置。
【請求項２】
前記制御部は、前記撮像期間内の第１区間では前記パルスＸ線を一定周期で発生させ、前記撮像期間内の第２区間では前記パルスＸ線を心電同期により不定周期で発生させることを特徴とする請求項１記載のＸ線診断装置。
【請求項３】
前記制御部は、前記被検体に対する造影剤の注入開始又は注入終了の時点を基準として前記第１区間と前記第２区間とを切り替えることを特徴とする請求項２記載のＸ線診断装置。
【請求項４】
前記制御部は、前記造影剤の注入開始又は注入終了の時点から所定の遅れ時間を経過した時点に前記第１区間と前記第２区間とを切り替えることを特徴とする請求項３記載のＸ線診断装置。
【請求項５】
前記制御部は、前記画像の濃度に基づいて前記第１区間と前記第２区間とを切り替えることを特徴とする請求項２記載のＸ線診断装置。
【請求項６】
前記第１区間は１〜５秒に予め設定され、前記第２区間は１〜６０秒に予め設定されることを特徴とする請求項２記載のＸ線診断装置。
【請求項７】
前記第１区間は任意の時間長に予め設定され、前記第２区間は任意の心拍数に予め設定されることを特徴とする請求項２記載のＸ線診断装置。
【請求項８】
前記第１区間では前記パルスＸ線は１０〜３０回／秒の頻度で発生され、前記第２区間では前記パルスＸ線は１〜５回／心拍の頻度で発生されることを特徴とする請求項２記載のＸ線診断装置。
【請求項９】
前記制御部は、前記撮像期間内の第１区間では前記パルスＸ線を心電同期により不定周期で発生させ、前記撮像期間内の前記第１区間に続く第２区間では前記パルスＸ線を一定の周期で発生させ、前記撮像期間内の前記第２区間に続く第３区間では前記パルスＸ線を心電同期により不定周期で発生させることを特徴とする請求項１記載のＸ線診断装置。
【請求項１０】
前記制御部は、前記被検体に対する造影剤の注入開始の時点を基準として前記第１区間から前記第２区間に切り替え、前記造影剤の注入終了の時点を基準として前記第２区間から前記第３区間に切り替えることを特徴とする請求項９記載のＸ線診断装置。
【請求項１１】
前記制御部は、前記撮像期間内の第１区間では前記パルスＸ線を第１周期で発生させ、前記撮像期間内の第２区間では前記パルスＸ線を前記第１周期と異なる第２周期で発生させることを特徴とする請求項１記載のＸ線診断装置。
【請求項１２】
前記制御部は、前記被検体に対する造影剤の注入開始又は注入終了の時点を基準として前記第１区間と前記第２区間とを切り替えることを特徴とする請求項１１記載のＸ線診断装置。
【請求項１３】
前記制御部は、前記撮像期間内の第１区間では前記パルスＸ線を第１周期で発生させ、前記撮像期間内の前記第１区間に続く第２区間では前記パルスＸ線を前記第１周期より短い第２周期で発生させ、前記撮像期間内の第３区間では前記パルスＸ線を前記第２周期より長い第３周期で発生させることを特徴とする請求項１記載のＸ線診断装置。
【請求項１４】
前記制御部は、前記撮像期間内の第１区間では前記パルスＸ線を心電同期により不定周期で発生させ、前記撮像期間内の前記第１区間に続く第２区間では前記パルスＸ線を一定周期で発生させ、前記撮像期間内の第３区間では前記パルスＸ線を心電同期により不定周期で発生させることを特徴とする請求項１３記載のＸ線診断装置。
【請求項１５】
前記制御部は、前記パルスＸ線の発生周期を、前記撮像期間内の操作者により指示された時点で延長又は短縮することを特徴とする請求項１記載のＸ線診断装置。
【請求項１６】
前記制御部は、前記第１、第２区間にわたって一定周期でＸ線検出動作を繰り返すために前記Ｘ線検出器を制御することを特徴とする請求項１記載のＸ線診断装置。
【請求項１７】
前記画像のデータを前記第１、第２区間を区別して記憶する記憶部をさらに備えることを特徴とする請求項２又は９記載のＸ線診断装置。
【請求項１８】
前記画像を前記第１区間と前記第２区間とで別々に再生する再生部をさらに備えることを特徴とする請求項１７記載のＸ線診断装置。
【請求項１９】
冠状動脈造影撮影シーケンスのもとで発生された複数の画像のデータを記憶する記憶部と、
前記複数の画像のデータに基づいて、心筋への血液供給領域に設定された基準領域に関する基準の時間濃度曲線と、心筋領域内に設定された複数の局所領域に関する複数の時間濃度曲線とを発生し、前記基準の時間濃度曲線に対する前記複数の時間濃度曲線の相関をそれぞれ表す複数の指標を演算する演算部とを具備することを特徴とする画像処理装置。
【請求項２０】
前記演算部は、前記心筋への造影剤流入期間における前記基準の時間濃度曲線に対する前記時間濃度曲線の相関を前記指標として前記局所領域ごとに演算することを特徴とする請求項１９記載の画像処理装置。
【請求項２１】
前記演算部は、前記心筋からの造影剤流出期間における前記基準の時間濃度曲線に対する前記時間濃度曲線の相関を前記指標として前記局所領域ごとに演算することを特徴とする請求項１９記載の画像処理装置。
【請求項２２】
前記演算部は、前記心筋への造影剤流入期間における前記基準の時間濃度曲線に対する前記時間濃度曲線の相関を表す第１指標と、前記心筋からの造影剤流出期間における前記基準の時間濃度曲線に対する前記時間濃度曲線の相関を表す第２指標とを前記局所領域ごとに演算することを特徴とする請求項１９記載の画像処理装置。
【請求項２３】
前記演算部は、前記第１指標と前記第２指標との比を前記局所領域ごとに演算し、前記比のマップを生成することを特徴とする請求項２２記載の画像処理装置。
【請求項２４】
前記比のマップを前記画像に重ねて表示する表示部をさらに備えることを特徴とする請求項２２記載の画像処理装置。
【請求項２５】
前記基準の時間濃度曲線を前記時間濃度曲線とともにグラフとして表示する表示部をさらに備えることを特徴とする請求項１９記載の画像処理装置。
【請求項２６】
前記演算部は、前記局所領域を前記第１指標の第１閾値に対する高低の区別と、前記第２指標の第２閾値に対する高低の区別との組み合わせに分類することを特徴とする請求項２２記載の画像処理装置。
【請求項２７】
前記演算部は、前記分類した組み合わせに関するマップを生成することを特徴とする請求項２６記載の画像処理装置。
【請求項２８】
前記演算部は、治療前に収集した前記心筋への造影剤流入期間における前記基準の時間濃度曲線に対する、治療後に収集した前記心筋への造影剤流入期間における前記基準の時間濃度曲線の相関を前記指標として前記局所領域ごとに演算することを特徴とする請求項１９記載の画像処理装置。
【請求項２９】
前記演算部は、治療前に収集した前記画像のデータから求めた前記指標に対する治療後に収集した前記画像のデータから求めた前記指標の比較結果を演算することを特徴とする請求項１９記載の画像処理装置。
【請求項３０】
前記指標は、前記心筋への造影剤流入期間における前記画像のデータから求められることを特徴とする請求項２９記載の画像処理装置。
【請求項３１】
前記指標は、前記心筋への造影剤流出期間における前記画像のデータから求められることを特徴とする請求項２９記載の画像処理装置。
【請求項３２】
前記演算部は、薬剤投与前に収集した前記画像のデータから求めた前記指標に対する薬剤投与後に収集した前記画像のデータから求めた前記指標の比較結果を演算することを特徴とする請求項１９記載の画像処理装置。
【請求項３３】
前記演算部は、前記指標に関するマップを生成することを特徴とする請求項１９記載のＸ線診断装置。
【請求項３４】
前記指標に関するマップを前記画像にオーバーレイして表示する表示部をさらに備えることを特徴とする請求項３３記載の画像処理装置。
【請求項３５】
前記演算部は、前記第１指標に関するマップ又は前記第２指標に関するマップを生成することを特徴とする請求項３０記載の画像処理装置。
【請求項３６】
Ｘ線を発生するＸ線管と、
被検体を透過したＸ線を検出するＸ線検出器と、
前記Ｘ線検出器の出力に基づいて画像のデータを発生する画像発生部と、
冠状動脈造影撮影シーケンスのもとで発生された複数の画像のデータに基づいて、心筋への血液供給領域に設定された基準領域に関する基準の時間濃度曲線と、心筋領域内に設定された複数の局所領域に関する複数の時間濃度曲線とを発生し、前記基準の時間濃度曲線に対する前記複数の時間濃度曲線の相関をそれぞれ表す複数の指標を演算する演算部とを具備することを特徴とするＸ線診断装置。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５Ａ】
【図１５Ｂ】
【図１６】
【図１７】
【図１８】
【図１９】
【図２０Ａ】
【図２０Ｂ】
【図２１】
【図２２】
【図２３Ａ】
【図２３Ｂ】
【図２４】
【図２５Ａ】
【図２５Ｂ】
【図２６】
【図２７】
【図２８】
【図２９Ａ】
【図２９Ｂ】
【図３０】
【図３１】
【図３２Ａ】
【図３２Ｂ】
【図３３】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５Ａ】
【図１５Ｂ】
【図１６】
【図１７】
【図１８】
【図１９】
【図２０Ａ】
【図２０Ｂ】
【図２１】
【図２２】
【図２３Ａ】
【図２３Ｂ】
【図２４】
【図２５Ａ】
【図２５Ｂ】
【図２６】
【図２７】
【図２８】
【図２９Ａ】
【図２９Ｂ】
【図３０】
【図３１】
【図３２Ａ】
【図３２Ｂ】
【図３３】
【公開番号】特開２０１３−１０６９９０（Ｐ２０１３−１０６９９０Ａ）
【公開日】平成２５年６月６日（２０１３．６．６）
【国際特許分類】
【出願番号】特願２０１３−４８４５２（Ｐ２０１３−４８４５２）
【出願日】平成２５年３月１１日（２０１３．３．１１）
【分割の表示】特願２００６−３５３４３０（Ｐ２００６−３５３４３０）の分割
【原出願日】平成１８年１２月２７日（２００６．１２．２７）
【出願人】（０００００３０７８）株式会社東芝 (54,554)
【出願人】（５９４１６４５４２）東芝メディカルシステムズ株式会社 (4,066)
【Ｆターム（参考）】
【公開日】平成２５年６月６日（２０１３．６．６）
【国際特許分類】
【出願日】平成２５年３月１１日（２０１３．３．１１）
【分割の表示】特願２００６−３５３４３０（Ｐ２００６−３５３４３０）の分割
【原出願日】平成１８年１２月２７日（２００６．１２．２７）
【出願人】（０００００３０７８）株式会社東芝 (54,554)
【出願人】（５９４１６４５４２）東芝メディカルシステムズ株式会社 (4,066)
【Ｆターム（参考）】
[ Back to top ]