(57)【特許請求の範囲】
【請求項１】
被検体を載置する静磁場に傾斜磁場を重畳すると共にＲＦ波を照射することによりデータ収集を行うデータ収集部と、
前記データ収集をフーリエ空間を構成する複数のボクセルグループを対象として複数回繰り返す撮影動作を撮影開始時のボクセルグループを変えながらＫ回（Ｋは２以上の整数）繰り返すように前記データ収集部を制御する制御部と、
前記Ｋ回の撮影動作により得られた前記複数のボクセルグループから補間処理により前記撮影動作の開始又は造影開始からの経過時間が略同一を示すフーリエ空間を生成し、前記生成したフーリエ空間を３次元離散フーリエ変換する画像生成部とを具備し、
前記制御部は、前記撮影動作において撮影が行われるサンプリング点を、前記撮影動作ごとに前記フーリエ空間で略一様にサンプリングするために前回までの撮影動作において撮影が完了したサンプリング点の間隔を埋めてようにいくように配置することを特徴とする磁気共鳴イメージング装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、回転撮影の可能なＸ線撮影装置及び磁気共鳴イメージング装置に関する。
【背景技術】
【０００２】
血管撮影（アンギオグラフィ:angiographyは、動脈に挿入したカテーテルを通して血管内に造影剤を注入し、Ｘ線撮影を行う撮影法である。この血管撮影に対応するＸ線撮影装置は、脳塞栓溶解術、動脈瘤閉塞術、狭窄血管拡張術などの血管内手術に必須の装置であって、精密な診断、血管内でのカテーテルの動作の観察、治療過程の観察、結果観察などに用いられる。また、診断だけの目的に用いられることもある。なお、以下に述べるデジタル・サブトラクション・アンギオグラフィ（ＤＳＡ；digital subtraction angiography）、ダイナミックＤＳＡ、回転ＤＳＡ、３Ｄ−ＤＳＡも含めて、「アンギオグラフィ（血管撮影）」と総称することが多い。
【０００３】
デジタル・サブトラクション・アンギオグラフィ；  アンギオグラフィにより造影剤を注入する前に第一の撮影を行ってデジタル画像（マスク画像）を得る。さらに、造影剤を注入して第二の撮影を行う。得られたデジタル画像からマスク画像を減算することによって、造影剤によって生じた像（すなわち血管の像）以外の像を消去し、血管像の観察を容易にする。
【０００４】
ダイナミック・デジタル・サブトラクション・アンギオグラフィ（ダイナミック-ＤＳＡ）；  ＤＳＡの第二の撮影において、造影剤を短時間かつ高濃度で注入（「ボーラス注入」と言う）し、その後、撮影を繰り返すことによって、造影剤が血流に乗って流れて行く様（動態）の時系列的なデジタル画像を得る。そして、これら一連のデジタル画像からマスク画像を減算することによって、血管の像以外の像を消去した一連のデジタル画像を生成する。
【０００５】
回転デジタル・サブトラクション・アンギオグラフィ（回転ＤＳＡ）；  デジタル・サブトラクション・アンギオグラフィ（ＤＳＡ）において、被写体の周囲を巡る軌道上の所定の角度範囲に渡る多数の異なる方向から第一および第二の撮影を行って、対応する方向から撮影した後者のデジタル画像から前者のデジタル画像を減算することによって、多数の異なる方向から撮影した血管像を得る。
【０００６】
三次元デジタル・サブトラクション・アンギオグラフィ（３Ｄ−ＤＳＡ: three dimensional ＤＳＡ）；  回転デジタル・サブトラクション・アンギオグラフィにおいて、円軌道上の１８０度を越える所定の角度範囲に渡る多数の異なる方向から第一および第二の撮影を行って回転ＤＳＡを得て、これらの画像を画素ごとに対数変換し、コーンビーム再構成法（Ｘ線ＣＴで用いられる画像再構成法やそのバリエーション）を適用して、血管の三次元データ（volume data）を再構成する。
【０００７】
図１には、回転ＤＳＡ、３Ｄ−ＤＳＡの可能な血管Ｘ線撮影装置（アンギオグラフィ装置）の撮影部を概念的に示している。Ｘ線焦点と撮像デバイスは図示しないアームに取り付けられており、これらが一体となって、回転中心軸Ｏを中心として回転できる。図１では、Ｘ線焦点が角度θの位置にあり、撮影範囲の中に置かれた被写体の影を撮像デバイス上に投影している。Ｘ線の広がり角度は２αである。Ｘ線焦点の角度をθ＝θ0から少なくともθ0＋π＋２αまで動かしながら繰り返し撮影を行うことによって、三次元画像再構成に必要な投影データを収集できる。
【０００８】
図２には回転ＤＳＡ撮影の手順を示している。  
［１］撮影準備  
（１）操作盤から撮影条件を入力し、撮影準備を指令すると、（２）制御装置はアーム駆動回路に指令を発して、撮影条件で定められた所定のアーム角度を設定するよう求める。アーム駆動回路はアーム回転モータを動作させてアームを回転させ、指令された回転角度になるようにする。（３）制御装置は高電圧発生装置に指令を発して、撮影条件に定められたＸ線条件（Ｘ線管装置の管電流、管電圧などの情報）を設定させる。（４）制御装置はオートインジェクタに指令を発して、撮影条件に定められた造影剤の流量、造影剤注入量を設定させる。
【０００９】
［２］撮影  
（１）操作盤から撮影開始を指示すると、（２）制御装置はオートインジェクタに指令を発して、造影剤の注入を開始させる。オートインジェクタは予め設定された注入量（ml）と注入時間（sec）に従って造影剤をカテーテルに注入する。（３）制御装置は造影剤注入開始後、撮影条件に定められた所定の遅延時間（d）だけ待ってから、アーム駆動回路に指令を発して、アームの回転を開始させる。（４）制御装置はアーム角度エンコーダを参照して、アームが撮影条件に定められた所定の角度に到達したときに高電圧発生装置に指令を送り、Ｘ線曝射を開始させる。（５）制御装置はアーム角度エンコーダを参照して、アームが撮影条件に定められた多数の所定の角度に到達するたびに、データ収集装置に指令を発して、Ｘ線像のデジタルデータを収集させる。収集したデータは記憶装置（図示しない）に蓄えられるとともに、データ処理装置（図示しない）によって、撮影条件に定められた所定の処理を受ける。
【００１０】
［３］撮影終了  
（１）制御装置はアーム角度エンコーダを参照して、アームが撮影条件に定められた所定の角度に到達したときに、高電圧発生装置に指令を送り、Ｘ線曝射を停止させる。さらに、アーム駆動回路に指令を送り、アームの回転を停止させる。
【００１１】
コーンビームＣＴ；  Ｘ線ＣＴであって、二次元の検出器を持ち、一度の撮影でボリュームデータを得られる装置である。  
ＣＴＡ（ＣＴ angiography）；  造影剤注入前のＣＴ画像をマスク画像とし、造影剤を動脈もしくは静脈から注入して撮影したＣＴ画像からマスク画像を減算して、血管のＣＴ画像を得る手法である。ヘリカルＣＴやコーンビームＣＴを用いてＣＴＡを行うと、血管の三次元データ（volume data）が得られる。コーンビームＣＴでは１回の撮影が０．５秒程度でできるので、繰り返し撮影を行って、時系列的なＣＴＡの三次元データを得ることが可能である。しかしながら、空間分解能（解像度）に限界があって、直径１ｍｍ以下の細い血管を正確に描出することは困難である。
【００１２】
従来の問題点は次の通りである。  
（１）時系列的な撮影による血流動態観察が必要な場合が多い。動脈から注入された造影剤は、末梢の毛細血管系を通過し、静脈に至る。アンギオグラフィにおいては、造影剤注入から撮影までの遅延時間を変える事によって、主に撮影される血管が動脈であるか、末梢血管系であるか、静脈であるかを選択する。ＤＳＡの第二の撮影においても、この事情は全く同じである。
【００１３】
しかし、血管内手術の対象となるような、血管系の一部が狭窄や閉塞を起こしている症例では、本来あるべき血流路が詰まっていたり、血流が正常の経路以外の経路（側副血行路）を経由して流れていたりすることがしばしばある。その病態の正確な診断には、血管系のどの部分が正常なタイミングで造影され（すなわち造影剤が到達し）、どの部分は正常よりも遅延するか、また、どの部分が先に造影され、どの部分が後で造影されるか、という動態を観察することが本質的に重要である。
【００１４】
ダイナミックＤＳＡを使えば、造影剤注入後の数秒間ないし十数秒間に渡って反復して撮影を行うことによって、造影剤が動脈、末梢血管系、静脈を通過する動態を時系列的なデジタル画像として観察できるので、この必要性を満たす事ができる。
【００１５】
一方、回転ＤＳＡや３Ｄ−ＤＳＡにおいては、被写体の周囲を巡る円軌道上の所定の角度範囲に渡る多数の異なる方向から撮影を行うために、一定の時間が必要で、典型的には１０秒程度掛かるため、血流動態観察に適さない。
【００１６】
（２）三次元再構成が必要な場合が多い。脳のアンギオグラフィでは、三次元的に屈曲した多数の血管が観察される。この画像の上で、小さい動脈瘤と血管の屈曲部を見分けること、また、血管像が重なり合った部分にまぎれた狭窄部位を発見することは容易でない。
【００１７】
回転ＤＳＡは被写体を囲む円軌道上の多様な方向から血管の像を観察できるので、この困難を若干緩和する。さらに３Ｄ−ＤＳＡを用いれば、volume dataが得られるので、これを任意の方向から観察したり、特定の部分（例えば観察したい部分以外の、重なり合っている血管）を除去した像を作って観察したりでき、小さい動脈瘤や、血狭窄部位を検索するのが容易である。また、直達手術を行う際の進入方向から見た像を作ることによって、手術前に外科医が術野の構造を理解するのに役立つ。
【００１８】
ただし、毛細血管系は非常に微細でしかも多数の血管から成るため、３Ｄ−ＤＳＡやＣＴを以てしても個別の血管を識別することは不可能である。しかし、３Ｄ−ＤＳＡやコーンビームＣＴを用いれば、局所の組織における、造影剤の平均濃度を定量的に測定することが可能である。すなわち、血中造影剤濃度が至るところ一様であったとすると、造影剤の平均濃度はその組織に占める毛細血管系の内腔体積の比率（血管床体積比率）、言い換えれば、当該局所組織の含む血液量（blood volume）に比例する。この血液量は梗塞の診断指標として重要である。
【００１９】
（３）高速回転する撮影装置は、カバーがないと危険であるが、カバーを持つ撮影装置は血管内手術や、開頭手術中の撮影に適さない。コーンビームＣＴにおいては、被写体を非常に多数の方向から短時間の内に撮影するために、撮影部が大きな角速度で回転する（１回転を０．５秒で行う装置もある）。このため、回転する機構部分をすべてカバーで覆って、安全性を確保することが必須である。
【００２０】
回転ＤＳＡは通常カバーを持たない。そして、被写体の周囲を巡る円軌道上から、所定の角度範囲に渡る多数の異なる方向から撮影を行うために、典型的には１０秒程度掛ける。これよりも短い時間で撮影を行うと、大きな角速度で装置を回転させることになり、危険性が高まって臨床に適さない。一方、通常のアンギオグラフィ装置は、いつでも撮影が開始できるように患者を位置決めした状態（撮影を行っていないとき）において、患者の周囲に比較的広い空間（open space）があって、カテーテル操作等を容易に行えることが特長である。また、装置が軽量で可動であるため、開頭手術等の外科手術の途中で撮影を行う際に、（患者を動かさずに）装置を移動して撮影することができる。
【００２１】
一方、撮影時間を短縮し、かつ安全性を確保するために、回転する機構部分をすべてカバーで覆い、短時間での撮影を可能にした回転ＤＳＡ装置の製品もかつて存在した。このような回転ＤＳＡ装置においては、カバーの一部にトンネル状もしくはくぼみ状の穴があり、患者の撮影対象部位をこの穴の中に入れて撮影を行う必要がある。このため、カテーテル操作等を容易に行える、手術中に撮影できる、という特長が損なわれてしまう。さらに、このような装置、あるいはＸ線ＣＴ装置は非常に高価であって、しかも広い設置面積を必要とする。アンギオグラフィを使って行う必要がある血管内手術は時間が掛かる場合が多く、一日に多数の患者を捌くことは難しいため、これらの装置をアンギオグラフィやＤＳＡの代用として用いかつ合理的な経済性を持って運用するのは困難である。
【００２２】
（４）高速回転する撮影装置は空間分解能に限界がある。脳などのアンギオグラフィにおいて行われる診断では、特に太い動脈に生じた動脈瘤、狭窄、閉塞を観察する事が重要であるが、開頭手術の準備として撮影を行う場合、あるいは血管内手術の際に撮影を行う場合には、しばしば、脳底部の太い動脈から生えている穿通枝動脈を観察することが重要になる。穿通枝動脈は脳深部にあって、重要な機能を持つ部分を司る部分を栄養しており、直径１ｍｍ以下のものが多く、また側副血行路（穿通枝動脈の血流が失われた際に代替の働きをする血管）がほとんどない。すなわち、もし閉塞したり、破裂したりすると重篤かつ不可逆的な症状を呈しやすく、致命的である場合も多い。
【００２３】
このため、開頭手術においては、誤って穿通枝動脈を切断しないように、また、誤って動脈瘤と一緒に穿通枝動脈をクリップで挟むことによって閉塞してしまうことのないように十分注意しなくてはならない。一方、血管内手術においては、太い動脈の内壁からはがれ落ちたアテローム（コレステロール塊）や凝固した血の塊が穿通枝動脈に流入して塞栓や狭窄を起こすことのないように十分注意しなくてはならず、もしそのような状況が生じたら、即座に塞栓物を溶解するなどの処置を行わなくてはならない。
【００２４】
しかし一般に、高速回転する撮影装置では、穿通枝動脈を確実かつ鮮明に描出することは難しい。これらの装置は空間分解能に限界があるという弱点を持からである。すなわち、高速回転で短時間の内に多数の方向から撮影を行うためには、高出力のＸ線管装置を搭載する必要がある。その場合、Ｘ線管装置の陰極と陽極の間に、短時間のうちに大きな電流が流れることになる。ところが、一般にＸ線管装置においては、電流の持つエネルギーのほとんど全てが、陽極において熱エネルギーに変わる。そして、陽極の表面が溶融・蒸散することを防ぐためには熱が発生する場所を集中させないことが必要不可欠であるから、回転陽極を用い、さらにＸ線焦点（陽極上でＸ線が発生する領域）の寸法を大きくせざるを得ない。このために、撮影されたＸ線像にぼけが生じる。この現象は、点光源による影が鮮明であるのに対して、広がりのある光源による影は不鮮明になるという幾何光学の原理的問題であり、回避できない。
【００２５】
さらにＸ線ＣＴ装置の場合には、撮像デバイスが多数のＸ線検出チャネルによって構成されている。このため、Ｘ線ＣＴ装置の空間分解能は、各Ｘ線検出チャネルのＸ線受光面の開口サイズによっても制約されている。すなわち、非常に小さい開口サイズを持つＸ線検出チャネルを用いれば空間分解能を向上させることができるが、そうすると、一定の撮影領域をカバーするには極めて多数のＸ線検出チャネルが必要になり、また、これらのＸ線検出チャネルから出力されるデータの量も膨大となって、装置が大規模なシステムになり、技術的に難しく、また非常に高価な装置になってしまう。
【００２６】
以上から、臨床での利用において要求される、（１）動態観察ができ、（２）三次元再構成ができ、（３）血管内手術、術中撮影に適し、（４）高い空間解像度を持ち、（５）経済的合理性を持つ、という条件を全て満たす手段は存在しない。
【発明の開示】
【発明が解決しようとする課題】
【００２７】
本発明の目的は、動態観察ができ、三次元再構成ができ、血管内手術や術中撮影に適し、高い空間解像度を持つ、経済的合理性を持つ、という条件を全て満たすダイナミック３Ｄ−ＤＳＡ対応のＸ線撮影装置及び磁気共鳴イメージング装置を提供することにある。
【課題を解決するための手段】
【００２８】
本発明のある局面によるＸ線撮影装置は、被検体を載置する静磁場に傾斜磁場を重畳すると共にＲＦ波を照射することによりデータ収集を行うデータ収集部と、前記データ収集をフーリエ空間を構成する複数のボクセルグループを対象として複数回繰り返す撮影動作を撮影開始時のボクセルグループを変えながらＫ回（Ｋは２以上の整数）繰り返すように前記データ収集部を制御する制御部と、前記Ｋ回の撮影動作により得られた前記複数のボクセルグループから補間処理により前記撮影動作の開始又は造影開始からの経過時間が略同一を示すフーリエ空間を生成し、前記生成したフーリエ空間を３次元離散フーリエ変換する画像生成部とを具備する。
【発明の効果】
【００２９】
本発明によれば、動態観察ができ、三次元再構成ができ、血管内手術や術中撮影に適し、高い空間解像度を持つ、経済的合理性を持つ、という条件を全て満たすダイナミック３Ｄ−ＤＳＡ対応のＸ線撮影装置及び磁気共鳴イメージング装置を提供することができる。
【発明を実施するための最良の形態】
【００３０】
以下、本発明の実施形態について図面を参照しながら説明する。本実施形態に係るＸ線撮影装置をダイナミック３Ｄ−ＤＳＡとして適用したとき、次の作用効果を奏することができる。  
（１）時系列的な一連の三次元再構成データ（volume data）が得られるので、血流の動態を観察することができる。  
（２）三次元再構成データが得られるので、これを任意の方向から観察したり、特定の部分、例えば観察したい部分以外の、重なり合っている血管を除去した像を作って観察したりできる。  
（３）撮影時に高速回転する必要がなく、従って、カバーを持たない撮影装置として構成する事ができる。このため、機構的には、従来のアンギオグラフィとほとんど同等の装置として構成でき、ゆえに、従来の３Ｄ−ＤＳＡ装置と同等程度の価格、設置面積の装置が実現できる。そして、いつでも撮影が開始できるように患者を位置決めした状態（撮影を行っていないとき）において、患者の周囲に比較的広い空間（open space）があって、カテーテル操作等を容易に行えるという、アンギオグラフィの特長をそのまま継承することができる。
【００３１】
（４）高い空間分解能が実現できる。１回のデータ収集（以下これを「ＲＵＮ」と呼ぶ）に掛かる時間は従来の３Ｄ−ＤＳＡ装置とほぼ同等であり、特に短時間で多量のデータを集中的に収集する必要がなく、従って、撮影時に高速回転する必要もない。このため、Ｘ線管装置の出力を抑えて、Ｘ線焦点のサイズが小さい装置を構成することができる。以上から、従来のアンギオグラフィ装置と同等の空間分解能が実現できる。  
（５）経済的合理性を持つ装置が構成できる。撮像デバイスとしてイメージインテンシファイアや平面撮像デバイスなど、主要な高額デバイスは、従来のアンギオグラフィ装置と同じものを使用できる。また、従来のアンギオグラフィの特長をそのまま継承しているために、アンギオグラフィ装置として利用でき、対象となる症例が幅広い。
【００３２】
また、本実施形態では多様な観察手段を提供することができ、それにより次の作用効果を奏することができる。本実施形態により撮影されたダイナミック３Ｄ−ＤＳＡは、ダイナミック３Ｄ−ＣＴＡで撮影されたデータと同等の多様な様式で観察できる。
【００３３】
つまり、注目する部位を囲む体積領域（Volume of interest: VOI）以外の血管や組織の像を消去して観察することによって、所望の部位における血流の動態を、多様な方向から観察することができる。  
また、任意の平面あるいは曲面を選択して、三次元データ選択した面上にあるボクセルだけを集めて平面画像を構成し、この平面画像を各時刻に置ける三次元データについて生成して、これらを系列的に並べることによって、時間軸と空間的平面から構成される三次元データを構成することができる。この三次元データ上では、選択した面上にある各所における造影剤の濃度変化が一覧できる。
【００３４】
また、選択した血管上の選択した２箇所に造影剤が到達した時刻を測定することによって、当該血管の血流速を測定することができる。  
さらに、選択した体積領域（Volume of interest: VOI）内の造影剤含有量の経時変化を測定して、これを動脈内の造影剤濃度変化の測定と比較することによって、当該VOIにおける、血流の平均通過時間と血液量が測定でき、これらから血流量が計算できる。血流量は虚血の診断指標として非常に重要である。
【００３５】
図３、図４に示すように、Ｘ線撮影装置は、Ｘ線管装置１２と撮像デバイス１４とを有する。Ｘ線管装置１２は、高電圧発生装置１３から管電圧の印加を受けてＸ線をファン角αで発生する。撮像デバイス１４は、典型的には、マトリクス状に配列された半導体検出素子を有するフラットパネルディテクタ（ＦＰＤ：平面型撮像デバイス）で構成される。データ収集装置２３は、被検体１５０を透過したＸ線のデータを撮像デバイス１４を介して収集する。
【００３６】
Ｘ線管装置１２は、撮像デバイス１４とともに、寝台１８の天板１７上の被検体１５０を挟んで互いに対向する向きでＣアーム５１に搭載される。Ｃアーム５１は、回転軸Ｚ１を中心に旋回自在に床面に支持される床旋回アーム５４から、スタンド５３、アームホルダ５２を介して支持される。スタンド５３は略鉛直な回転軸Ｚ２まわりに回転自在に設けられ、アームホルダ５２は略水平な回転軸Ｚ３まわりに回転自在に設けられる。Ｃアーム５１は、アームホルダ５２に、回転軸Ｚ４まわりにスライド回転自在に設けられる。これら回転軸Ｚ１，Ｚ２，Ｚ３，Ｚ４各々周りの回転は、アーム駆動回路１６に駆動されるアーム回転モータ１８により個別に駆動される。アーム角度エンコーダ１９は、回転軸Ｚ１，Ｚ２，Ｚ３，Ｚ４周りの回転角度を個別に検出するために設けられている。
【００３７】
制御装置２０は、システム全体の制御とともに、撮影制御を担う。撮影制御に際しては、制御装置２０は、操作盤２１を介して管電圧、管電流、撮影角度範囲、サンプリングピッチ（撮影間隔）Δθ、回転撮影の回数、造影開始から撮影開始までの遅延時間等の撮影条件を決定する。操作盤２１を介して撮影トリガが入力される。制御装置２０は、撮影条件に従って撮影が行われるように高電圧発生装置１３、アーム回転モータ１６、オートインジェクタ２２を制御する。
【００３８】
上記装置により、回転アンギオグラフィを繰り返し撮影すること（造影は毎回の回転撮影ごとに行う）によって収集したデータに基づいて、任意の時刻における３Ｄアンギオ像（３Ｄ×時刻＝４Ｄ）を生成する。この機能を実現するための制御方法が従来とは違っている。
【００３９】
本実施形態の構成では、図５に示すように、  
（１）Ｋ回（Ｋは２以上の整数）造影剤注入を行い、その都度１回ずつ回転撮影を行う。  
（２）毎回の回転撮影（ＲＵＮ）は、造影開始後所定の遅延時間Ｔｓが経過した時点から、造影開始後所定の時間Ｔｅ（Ｔｅ＞Ｔｓ）が経過した時点まで行う。毎回の回転撮影において、Ｘ線像の撮影はアーム５１の回転角度が所定のサンプリングピッチΔθの整数倍になった毎時に行う。  
（３）回転撮影を行う際のアームの回転方向は全て同じである。  
（４）回転撮影開始時のアーム５１の角度（回転撮影開始角度）は毎回異なる。第ｋ回目の回転撮影の撮影開始時のアームの角度をｓ（ｋ）とするとき、Ｋ回の撮影によってs（Ｋ）（Ｋ＝１, ２, …, Ｋ）が、所定の撮影角度範囲の中に概ね一様に分布するようになるよう、s（Ｋ）を決める。撮影角度範囲は、いわゆるハーフ再構成の可能な角度範囲（π＋２α以上、αはＸ線のファン角）に設定される。
【００４０】
実際に（４）を実現するための最も単純な方法は、ｋ回目の撮影（ＲＵＮ ｋと呼ぶ）における、撮影開始時のアームの角度をs（ｋ）とし、βをβ＞π＋２αである定数とするとき、  
s（ｋ＋１）＝s（ｋ）＋β/（ｋ-１）  
とすることである。
【００４１】
図５は、このようにしてs（Ｋ）を決定したものである。より実用的には、Ｋ回の回転撮影を予定していても、何らかの事情（たとえば患者の容態が悪化する、患者が動くなど）によって途中で撮影を止めなくてはならないことを考慮して、途中で撮影をやめた場合にでも（品質は多少悪くなっても）任意の時刻における３Ｄアンギオ像が得られるようにs（ｋ）を設定することが望ましい。
【００４２】
画像生成装置２４では、Ｋ回の回転撮影（ＲＵＮ）で実際に収集されたＸ線像群Ｑ（θ,ｋ）, ｋ＝１,２,…,Ｋ（黒丸で示す）から任意の時刻tにおけるＸ線像群（白抜き四角）を補間により生成し、各撮影角度ａで撮影されたであろうＸ線像群Ｐ（θ,t）が生成される。画像Ｐ（ａ,t）は、同じ撮影角度aにおける最寄りのデータから推算によって生成する。推算された画像の撮影角度の範囲がπ＋２α以上あれば、これらの推算された画像から、時刻tにおける３Ｄ−ＣＴＡ画像を再構成できる。
【００４３】
図６に示すように、補間の最も簡単な方法は、一次補間である。一次補間の場合、実際に撮影できたＸ線像Ｑの中から、アームの回転角度が同じａであって、しかも造影開始後の経過時間がtに最も近いＸ線像Ｑ（a,Ｋ＋１）とＱ（a,Ｋ）から、Ｐ（a,t）を推算する。なお、ここで言うＸ線像とは、Ｘ線撮影で得られた平面画像を対数変換したものを指す。
【００４４】
図７には、一次補間によってＰ（a,t）を推算する場合に、計算に使用される「アームの回転角度が同じaであって、しかも造影開始後の経過時間がtに最も近いＸ線像」を丸印で囲んで示している。
【００４５】
一次補間の方法としては、Ｘ線像はデジタル画像であり、画素（pixel）が二次元に配列されたもので構成される。各画素は濃度値という一つの数値を持つ。そこで、図８に示すように、Ｘ線像Ｑ（a,Ｋ）上の、ヨコ方向の座標がｘ、タテ方向の座標がｙの位置にあるひとつの画素の値をＱ（a,Ｋ）[x,y]と表すことにする。Ｑ（a,Ｋ＋１）[x,y], Ｐ（a,t）[x,y]も同様の意味である。Ｘ線像Q（a,Ｋ）が撮影された時の、造影開始後の経過時間をu、Ｘ線像Q（a,Ｋ＋１）が撮影された時の、造影開始後の経過時間をvとする。
【００４６】
このとき、P（a,t）[x,y]を以下のようにして推算する。  
P（a,t）[x,y]＝（（t-u） Q（a,Ｋ＋１）[x,y] ＋ （v-t） Q（a,Ｋ）[x,y]）/（v-u）
こうして、所望のＸ線像（すなわちアームの回転角度がaであって造影開始後の経過時間がtであるＸ線像）の推定P（a,t）が得られる。
【００４７】
このようにして、造影開始後の経過時間tの任意の値について、アーム回転角度θの範囲がπ＋２α以上に渡って、画像Ｐ（θ,t）が推算できる。そこで、従来の方法（回転アンギオグラフィにおける、コーンビーム投影からの３次元再構成法）を用いて、Ｐ（θ,t）（θはπ＋２α以上の範囲に渡る）から３次元再構成を行って、造影開始後の経過時間がtである時点における３次元像を再構成できる。
【００４８】
従って、繰り返し撮影を行って得たデータひと組から、造影開始後の経過時間tの異なる一連の３次元像が得られる。
【００４９】
（２）毎回の回転撮影（ＲＵＮ）は、造影開始後所定の遅延時間Ｔｓが経過した時点から、第二の所定の遅延時間Ｔｅが経過した時点まで行う。毎回の回転撮影において、Ｘ線像の撮影はアームの回転角度が所定の定数Δθの整数倍になった時に行う。
【００５０】
ここで、Δθは回転角度のサンプリングピッチである。臨床使用に耐える画質の３Ｄ像を得るには、Δθは比較的小さく例えば３度が必要とされる。それよりも大きいと、再構成された３Ｄ像に含まれる偽像（artifact）が強くなって、診断を妨げる。
【００５１】
また、もしアーム５１が常に一定の角速度で回転するのであれば、（a）「毎回の回転撮影において、Ｘ線像の撮影はアームの回転角度が所定の定数Δθの整数倍になった時に行う」という規則と（b）「毎回の回転撮影において、Ｘ線像の撮影は所定の時間Δt周期で行う」という規則とは同じことである。しかし、実際の装置では角速度は誤差を持つ。すなわち、回転部分の荷重の偏り、モータ１８の出力のふらつき、回転の加速・減速に伴う角速度変動などの要因によって、必ずしも正確に一定の角速度を維持できるとは限らず、さらに、Ｋ回の繰り返し撮影に於いてどの回転撮影（ＲＵＮ）でも同じ角速度になるようにできるとは限らない。
【００５２】
一方、本実施形態で開示する補間法の方式が要請するのは、所定の時間間隔で撮影された一連のＸ線像が得られていることではなく、所定のアーム回転角度において撮影された一連のＸ線像が得られていることである。
【００５３】
従って、上記のような角速度の誤差があっても適切な撮影ができるように構成するためには、（a）の規則に従う必要があり、規則（b）では旨く行かない。
【００５４】
より実用的には、Ｋ回の回転撮影が完遂できなかった場合にも撮影が無駄にならないように構成するのが望ましい。Ｋ回の回転撮影を予定していても、何らかの事情（たとえば患者の容態が悪化する、患者が動くなど）によって途中で撮影をやめなくてはならないことがあり得る。途中で撮影をやめた場合にでも（品質は多少悪くなっても）任意の時刻における３Ｄアンギオ像が得られるようにs（Ｋ）を設定することが望ましい。
【００５５】
このためには、Ｋ回のうち最初の任意のＬ回（Ｌ≦Ｋ）の撮影における撮影開始時のアーム５１の角度ｓ（１）〜ｓ（Ｌ）が、角度範囲βの中に概ね一様に、かつ、角度範囲β全般に渡るように分布するよう、撮影開始時の角度ｓ（１）〜ｓ（Ｌ）を決定すれば良い。
【００５６】
具体的には、撮影開始時のアーム５１の角度を決定する方法は、例えば以下のように構成しても良い。  
１回目の撮影はs（１）＝0とする。
【００５７】
２回目の撮影はs（２）＝βとする。
【００５８】
３回目以降Ｋ回目の撮影のs（Ｋ）は、s（１）, s（２）, …, s（Ｋ-１）を小さい順に並べたときに、隣り合うものの間隔が最も大きい箇所に入るようにし、しかも、その間隔をa:（１-a）に内分する値とする。（ここにaは0.３５〜0.４８程度の定数である。）
例えばa＝0.４５とする。 s（１）＝0, s（２）＝β, s（３）＝0.４５β, s（４）＝0.67５β であるとき、小さい順に並べると、  
s（１） ＜ s（３） ＜ s（４） ＜ s（２）
であり、隣り合うものの間隔が最大なのはs（１）とs（３）の間（0.４５β）である。そこで、s（５）はs（１）とs（３）の間隔を0.４５:0.５５に内分する値、すなわち、   s（５） ＝ s（１）＋0.４５（s（３）-s（１））＝0.２0２５β  
と計算する。
【００５９】
図９は、上記手続きを繰り返すことによって、Ｋ回目の撮影までに、角度範囲βがどのようにサンプリングされていくか、その進行の様子を示している。図９では３０回撮影する場合までを示している。例えば、当初３０回撮影する予定でありながら、実際には１５回しか撮影できなかったとしても、角度範囲βをほぼ一様にサンプリングしたデータが得られる。（横軸が１５の箇所に縦に並んだ点列が概ね一様に分布している。）つまり、各回の撮影は角度範囲βをほぼ一様にサンプリングするとともに、直前までの撮影によるサンプリング点を補完するようにサンプリングする。要するに、各回転撮影動作において撮影が行われるサンプリング点は、それまでの回転撮影動作において撮影が完了したサンプリング点の間隙を埋めていくように配置される。
【００６０】
撮影対象組織における造影剤の濃度の経時的変化は、心拍の影響を受ける。この効果は造影剤を動脈から注入する場合に著しい。特に造影開始直後は、血中造影剤濃度の変化が速いので、心拍の影響が大きくなる。本実施形態では、Ｋ回繰り返して回転撮影を行うにあたって、造影剤注入開始からの経過時間につれて撮影対象組織における造影剤の濃度が変化する現象が、毎回ほぼ同じように起こることを利用している。従って、より精度の高い３Ｄ画像を得るためには、図１０に示すように、毎回の撮影を心拍と同期して行うことが望ましい。
【００６１】
既存の心電計の多くは、心電図のＲ波のピークを自動検出する機能を持つ。図１０では、患者に心電計を接続してＲ波がピークになった時刻を常時モニタする。そして、操作盤１５から撮影開始ボタンを押すなどして、Ｋ回目の撮影開始を指示したあと、最初にＲ波のピークが出現した時刻を基準にして、それから所定の遅延時間ｖが経過した時点で造影剤の注入を開始する。そして、造影開始から所定の遅延時間ｄが経過した時点で撮影を開始する。
【００６２】
この構成によれば、もし患者の（毎分当りの）心拍数が安定しているならば、繰り返して回転撮影を行ううちの毎回の撮影において、心拍の影響がどれもほぼ同程度になるようにすることができる。
【００６３】
本装置が何回転も続けて回転できる構造を有するとき、図１１に示すように、複数の周回を重ねながら撮影を繰り返すことができる。この場合、撮影は、アームの回転角がΔθの整数倍に至った各時点で行う。この構成では、少なくとも４π（２回転）の回転ができるアンギオグラフィ装置を使用し、Ｋ回の繰り返し撮影を行う。  
（１）Ｋ回造影剤注入を行い、その都度１回ずつ回転撮影を行う。  
（２）毎回の回転撮影（ＲＵＮ）は、造影開始後所定の時間Ｔｓが経過した時点から、第二の所定の時間Ｔｅが経過した時点まで行う。毎回の回転撮影において、Ｘ線像の撮影はアームの回転角度が所定の定数Δθの整数倍になった時に行う。  
（３）回転撮影を行う際のアームの回転方向は全て同じである。
【００６４】
（４）回転撮影開始時のアームの角度は毎回異なる。第Ｋ回目の回転撮影の撮影開始時のアームの角度をs（Ｋ）とするとき、Ｋ回の撮影によってs（Ｋ）（Ｋ＝１, ２, …, Ｋ）が角度の範囲２πの中に概ね一様に分布するようになるよう、s（Ｋ）を決める。
【００６５】
実際に（４）を実現するための最も単純な方法は、Ｋ回目の撮影（ＲＵＮ Ｋと呼ぶ）における、撮影開始時のアームの角度をs（Ｋ）とするとき、  
ｓ（Ｋ＋１）＝ｓ（Ｋ）＋２π／Ｋ  
とすることである。しかしより実用的には、Ｋ回の回転撮影を予定していても、何らかの事情（たとえば患者の容態が悪化する、患者が動くなど）によって途中で撮影を止めなくてはならないことを考慮して、途中で撮影をやめた場合にでも（品質は多少悪くなっても）任意の時刻における３Ｄアンギオ像が得られるようにs（Ｋ）を設定することが望ましい。この方法については、図５でβ＝πとしたものと同じである。
【００６６】
収集したＸ線像（対数変換したもの）から、所望のＰ（a,t）を推算する方法も補間法として上述した通りである。
【００６７】
アーム５１の１回転分に渡る回転角度ａについて、造影開始後ｔ秒におけるＸ線像が推算できるため、コーンビームからの再構成法を用いて、造影開始後任意の時間ｔ秒（Ｔｓ＜ｔ＜Ｔｅ）が経過した時点における３Ｄ像を生成できる。以上をまとめると；  
（１）Ｋ回造影剤注入を行い、その都度１回ずつ回転撮影を行う。  
（２）毎回の回転撮影（ＲＵＮ）は、造影開始後所定の時間Tsが経過した時点から、第二の所定の時間Teが経過した時点まで行う。毎回の回転撮影において、Ｘ線像の撮影はアーム５１の回転角度が所定の定数Δθの整数倍になった時に行う。  
（３）回転撮影を行う際のアーム５１の回転方向は必ずしも同じでない。順逆の両方向で交互に回転撮影を行う。  
（４）回転撮影開始時のアーム５１の角度は毎回異なる。第Ｋ回目の回転撮影の撮影開始時のアーム５１の角度をs（Ｋ）とするとき、Ｋ回の撮影によってｓ（Ｋ）（Ｋ＝１, ２, …, Ｋ）が角度の範囲２πの中に概ね一様に分布するようになるよう、ｓ（Ｋ）を決める。
【００６８】
図５の例と異なるのは動作（３）である。例えば、少なくとも４πの回転ができるアンギオグラフィ装置を使用し、奇数回目のＲＵＮと偶数回目のＲＵＮとで回転方向が逆になるように構成することができる。
【００６９】
Ｋ回目の回転撮影ＲＵＮ  Ｋのあと、Ｋ＋１回目のＲＵＮに先立って、Ｋ＋１回目のＲＵＮの開始位置近辺まで帰還させるためにアームを逆回転させる必要があり、図５の例ではこの帰還の回転はかなり大きい。本構成では、逆回転の最中にもＸ線像の収集を行えるので、検査時間を短縮でき、さらにアームの動きによる事故の危険性を低減する効果が期待できる。
【００７０】
しかしこの構成は、図１２からわかるように、推算しようとするＰ（a,t）において、例えば（t,a）のある組み合わせでは、推算に用いるＸ線像が近傍にないために、推算の精度が落ちるという難点も持っている。
【００７１】
一方、図１３に示すように、Ｋ回繰り返して回転撮影を行う代わりに、アーム５１の回転角Ｋ通りについて、回転をせずに繰り返し撮影を行うことによっても、造影開始後からｔ秒経過した時点での３Ｄ像を得られる。これは原理的には正しい。しかし、実用的な画質の３Ｄ像を得るには、Δθはせいぜい３度以内でなくてはならない。例えばΔθを３度とし、βを２１０度とすると、Ｋ＝７０である。すなわち、７０回に渡って造影剤注入と撮影を繰り返す必要がある。この図１２の例と対比すると、図５の例ではＫ回の繰り返し撮影（ＲＵＮ）のそれぞれにおいて、アーム５１の回転角度を変化させながら撮影を行うという構成の優位性が明らかである。
【００７２】
本実施形態では、Ｋ回の回転撮影を行う。従って、もし毎回の回転撮影（ＲＵＮ）を、従来の回転アンギオグラフィ検査と同等のＸ線量で行ったとすると、被験者の被曝するＸ線量はＫ倍になってしまう。そこで、従来の回転アンギオグラフィ検査と同等の被曝線量で検査を行おうとすると、毎回の回転撮影（ＲＵＮ）で照射するＸ線量を１／Ｋにする必要がある。
【００７３】
一般に、Ｘ線像（対数変換したもの）の各画素の値に含まれるノイズ成分は、主に「量子ノイズ」と呼ばれる、Ｘ線光子の発生量のばらつきに起因するノイズである。そして、照射するＸ線量をｘ倍にすると、ノイズ成分の振幅は１／ｘ１／２倍になる。従って、Ｘ線量を１／Ｋにするとノイズの振幅はＫ１／２倍に増えることになる。ノイズの振幅があまり大きいと、３Ｄ像のしたときに偽像（artifact）が多くなり、診断の妨げになる。
【００７４】
ところが、本実施形態において、図５の例でβ＝２πの場合や図１１、図１２の場合（もともとβ＝２πである）には、それぞれのアーム角度ａについて、同じアーム角度で撮影した多数のＸ線像が得られる。そしてこれらは、アーム角度をａに固定したまま反復撮影を行った場合に得られる一連のＸ線像と同等である。このため、動的画像（ダイナミック image）一般に対して適用可能なノイズ低減処理法、例えばコヒーレントフィルタを用いて、解像度を維持したままノイズを著しく低減できる。従って、本実施形態とノイズ低減処理法を組み合わせることによって、Ｋ回の回転撮影による被曝は、通常の回転アンギオグラフィ検査における被曝のＫ倍よりも小さくでき、実際上はほぼ同等程度の被曝に抑えることが可能である。
【００７５】
次に、本実施形態の磁気共鳴イメージング装置（ＭＲＩ）への適用例を説明する。ＭＲＩで３Ｄ像を得るためには、フーリエ変換法を用いるのが一般的である。すなわち、図１４に示すように、フーリエ空間を多数のボクセル（voxel）に分け、それぞれのボクセルに測定データ（複素数値）を割り当てる。ただし、一度の測定ではフーリエ空間の全ボクセルに値を割り当てることはできないので、フーリエ空間のボクセルをＮ個のグループ（グループ番号0, １, ２, … , N-１）に分け、それぞれのグループごとに、値を割り当てるための測定を行う。従って、測定をＮ回繰り返すことになる。この繰り返しの周期は繰り返し時間ＴＲと呼ばれている。例えばグループ数がＮ＝１２８で、繰り返し時間がＲ＝0.１秒の場合、撮影に要する時間（撮影時間）はＮ×ＴＲ＝１２.８秒である。
【００７６】
フーリエ空間のすべてのボクセルについて値が割り当てられたら、フーリエ空間を３次元離散フーリエ変換によって変換すると、３Ｄ像（これもボクセルから構成される）が得られる。
【００７７】
測定データは傾斜磁場の一連の操作とＲＦ波の照射を組み合わせて行い、被検体から発生するエコー波を測定して得る。このとき、値を割り当てようとするグループに対応して、傾斜磁場の一連の操作の仕方が決まる。
【００７８】
図１５に本実施形態をＭＲＩに適用する例を示す。この例では、アーム５１の回転角度の代わりにグループ番号を縦軸として表す。  
（１）Ｋ回造影剤注入を行い、その都度１回ずつ連続撮影を行う。  
（２）毎回の連続撮影（ＲＵＮ）は、造影開始後所定の時間Ｔｓが経過した時点から、第二の所定の時間Ｔｅが経過した時点まで行う。毎回の連続撮影は、周期ＴＲで繰り返される多数回のデータ収集から成る。一度のデータ収集で、一つのボクセルグループに対応する測定値が決定される。  
（３）連続撮影ごとに、Ｎ個あるボクセルグループに対応する測定をどの順で実施するかは異なる。
【００７９】
ＭＲＩ装置ではどのボクセルグループに対応するデータを収集するかは撮影シーケンス（傾斜磁場とＲＦ波の操作の組み合わせ）によって決まる。従ってアンギオグラフィ装置を使った例とは違って、撮影に際して機械的可動部を用いる必要がないので、角速度の誤差に関する議論に対応する問題はもともと存在しない。造影開始後Ｔｓ秒経過した時点以降、一定の周期（ＴＲ秒）でデータ収集を行うことは容易である。
【００８０】
図１５では、１回目の連続撮影（ＲＵＮ１）ではボクセルグループＳ（１）に対応する測定から始めて、ＴＲ秒後にはボクセルグループＳ（１）＋１、さらにＴＲ秒後にボクセルグループＳ（１）＋２と測定を繰り返し、ボクセルグループＮ−１に対応する測定した次には、ボクセルグループ0に対応する測定する。２回目の連続撮影では、ボクセルグループＳ（２）に対応する測定から始め、以下同様である。
【００８１】
要するに、被検体を載置する静磁場に傾斜磁場を重畳すると共にＲＦ波を照射することによりデータ収集を行うデータ収集部を、データ収集をフーリエ空間を構成する複数のボクセルグループを対象として複数回繰り返す撮影動作を撮影開始時のボクセルグループを変えながらＫ回（Ｋは２以上の整数）繰り返すように制御部で制御する。
【００８２】
このようにして、各ボクセルグループについて、多数の測定データが得られる。次に、造影開始後ｔ秒（ｔはＴｓ＜ｔ＜Ｔｅであれば任意）における３Ｄ像を生成する方法を述べる。各ボクセルグループの造影開始後ｔ秒における対応する測定値を、既に述べたような補間法を用いて推算する。こうして得たフーリエ空間を３次元離散フーリエ変換すると、造影開始後ｔ秒における３Ｄ像が得られる。
【００８３】
心拍に同期した撮影開始タイミングの決定を行うことが望ましいのも、アンギオグラフィ装置を用いた例と同様である。
【００８４】
なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【図面の簡単な説明】
【００８５】
【図１】回転ＤＳＡ、３Ｄ-ＤＳＡの可能なアンギオグラフィ装置の撮影部の概念図。
【図２】回転ＤＳＡ、３Ｄ-ＤＳＡの撮影動作の説明図。
【図３】本実施形態において、回転ＤＳＡ、３Ｄ-ＤＳＡができるアンギオグラフィ装置の構成を示す図。
【図４】図３のアンギオグラフィ装置のガントリの外観を示す図。
【図５】図３の制御装置の制御による複数回の回転撮影で得られるデータの配置を示す図。
【図６】図５の配置の一部を拡大図。
【図７】図３の画像生成装置による一次補間処理に関する説明補足図。
【図８】図３の画像生成装置による一次補間処理に関する説明補足図。
【図９】図３の制御装置の制御による複数回の回転撮影で得られるデータの配置の他の例を示す図。
【図１０】図３の制御装置の制御による心拍同期撮影開始タイミングの決定に関する説明補足図。
【図１１】図３の制御装置の制御による連続的な回転撮影で得られるデータの配置例を示す図。
【図１２】図３の制御装置の制御による順逆交互の回転撮影で得られるデータの配置例を示す図。
【図１３】図３の制御装置の制御による回転を伴わない撮影で得られるデータの配置例を示す図。
【図１４】本実施形態において、ＭＲＩへの適用例の説明図。
【図１５】図１４のＭＲＩにおけるデータ配置例を示す図。
【符号の説明】
【００８６】
１２…Ｘ線管装置、１３…高電圧発生装置、１４…撮像デバイス、１６…アーム駆動回路、１７…天板、１８…寝台、１９…アーム角度エンコーダ、２０…制御装置、２２…オートインジェクタ、２３…データ収集装置、２４…画像生成装置、１５０…被検体、５１…Ｃアーム、５２…アームホルダ、５３…スタンド、５４…床旋回アーム。
【図１】
0005684329.tif 000002
【図２】
0005684329.tif 000003
【図３】
0005684329.tif 000004
【図４】
0005684329.tif 000005
【図５】
0005684329.tif 000006
【図６】
0005684329.tif 000007
【図７】
0005684329.tif 000008
【図８】
0005684329.tif 000009
【図９】
0005684329.tif 000010
【図１０】
0005684329.tif 000011
【図１１】
0005684329.tif 000012
【図１２】
0005684329.tif 000013
【図１３】
0005684329.tif 000014
【図１４】
0005684329.tif 000015
【図１５】
0005684329.tif 000016