JP3901448B2 - 磁気共鳴イメージング装置 - Google Patents
磁気共鳴イメージング装置 Download PDFInfo
- Publication number
- JP3901448B2 JP3901448B2 JP2000381171A JP2000381171A JP3901448B2 JP 3901448 B2 JP3901448 B2 JP 3901448B2 JP 2000381171 A JP2000381171 A JP 2000381171A JP 2000381171 A JP2000381171 A JP 2000381171A JP 3901448 B2 JP3901448 B2 JP 3901448B2
- Authority
- JP
- Japan
- Prior art keywords
- measurement
- magnetic resonance
- resonance imaging
- imaging apparatus
- region
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Description
【発明の属する技術分野】
本発明は、核磁気共鳴現象を利用して被検体の断層画像を得る磁気共鳴イメージング装置(MRI)に関し、特に、磁気共鳴イメージング装置おけるダイナミック磁気共鳴アンジオグラフィ撮像に好適な計測方法に関するものである。
【0002】
【従来の技術】
MRI(磁気共鳴イメージング装置)は、NMR現象を利用して、被検体中の所望の検査部位における原子核スピンの密度分布、緩和時間分布等を計測して、その計測データから被検体の任意の断面を画像表示するものである。
【0003】
また、このような磁気共鳴イメージング装置を利用して血流画像を得る技術は、MRA(磁気共鳴アンジオグラフィ)と呼ばれる。このMRAは、造影剤を使用しない手法と造影剤を使用する手法に大別することができる。
【0004】
造影剤を使用しないMRAの手法としては、スライス面への流入効果を利用してて血流画像を抽出するタイムオブフライト(Time-Of-Flight; TOF)法、血流による位相拡散の有無を利用して血流画像を抽出する差分処理を施すフェイズセンシティブ法(Phase-Sensitive:PS法)、血流による位相拡散の極性を反転し血流画像を抽出する差分処理を施すフェイズコントラスト法(Phase-Contrast:PC法)などが知られている。
【0005】
なお、これらの手法については、"Magnetic Resonance Imaging; Stark DD et al. edited, The C.V.Mosby Company, pp108-137, 1988"、"Cerebral MR Angio-imaging(脳血管時期共鳴画像法)の研究-第1報-、福井啓二 他、CT研究、10(2) 1988年、133頁-142頁"、"Magnetic resonance angiography, Dumoulin CL, et al, Radiology 161:717-720, 1986"、"Simultaneous acquisition of Phase-contrast angiograms and stationary-tissue images with Hadamard encoding of flow-induced phase shifts, Dumoulin CL, et al. J. Magnetic Resonance Imaging 1:399-404, 1991"、"Encoding strategies for three-direcrion Phase-contrast MR imaging, Dumoulin CL, et al. J. Magnetic Resonance Imaging 1:405-413, 1991"等に詳しい。
【0006】
一方、造影剤を用いるMRAの手法としては、Gd-DTPA等のT1短縮型の造影剤とグラジエントエコー系のTRの短いパルスシーケンスを組み合わせる手法が一般的である。この手法では、数ms〜数10msの短時間間隔の高周波磁場による同一領域の励起の繰り返しによって、血流以外の部分のピンを飽和させ、造影剤を含むためにT1が短縮し、スピンの飽和が生じにくい血流部分からのみ高信号を得ることにより血流画像を抽出する。または、さらに、造影前後の画像差分をとるDSA; Digital Subtraction Angiographyによって、より血流部分のコントラストを上げた画像を得る。ここで、このような造影剤を用いるMRAの手法では、通常肘静脈から造影剤を注入する。注入された造影剤は、心臓、動脈系、毛細血管、静脈系と順に循環する。そこで、ダイナミックMRAと呼ばれる技術では、造影剤の循環の各段階において繰り返し計測を行い、各部の血流の時系列画像が得られるようにしている。なお、このような造影剤を用いるMRAの手法やダイナミックMRAについては、"3D Contrast MR Angiography 2nd edition. Prince MR, Grist TM and Debatin JF, Springer, pp3-39,1988"に詳しい。なお、ダイナミックMRAについては特に本書のpp16-19に詳しい。
【0007】
さて、このようなダイナミックMRAにおける1回の撮像時間は、2次元計測の場合は繰り返し時間TRと位相エンコードステップ数を乗じた時間となり、3次元計測の場合は繰り返し時間TRと位相エンコードステップ数にスライスエンコードステップ数を乗じた時間となる。したがって、空間分解能を向上するためには位相エンコードステップ数やスライスエンコードステップ数を大きくとることが望ましい一方、そのようにすると時間分解能は低下する。すなわち基本的に空間分解能と時間分解能はトレードオフの関係にある。
【0008】
そこで、このような問題を解決するために、米国特許5,713,358号、米国特許5,830,143号では、次のような計測方法を提案している。
すなわち、これらの技術では、一度の計測でk空間全てを計測するのではなく、k空間をky方向に複数の領域に分割し、分割した領域を単位として、k空間の中心領域(低周波数領域)がより高い頻度で計測されるように、各回で計測するk空間の部分を制御する。そして、今回計測しなかったk空間の部分については、他の回に計測した計測結果を流用する。または、他の回に計測した複数の計測結果から補間によって生成する。
【0009】
この技術によれば、各回の計測において位相エンコードステップの全てを実行しないために、各回の計測時間を短縮化し時間分解能を向上することができると共に、比較的、診断に重要なコントラストを決定するk空間の中心領域(低周波領域)については経時的変化を確実に捉えることができる。
【0010】
【発明が解決しようとする課題】
しかしながら、前記米国特許5,713,358号、米国特許5,830,143号の技術では、高周波数領域の時間分解能の犠牲が著しく大きかった。また、k空間の分割の仕方がky方向のみと一次元方向のみであるために、計測結果に方向依存性が生じ、これによって生成される画像の品質を劣化する場合があった。
【0011】
そこで、本発明は、MRIにおいて繰り返し撮像を行う際に、高周波数領域の時間分解能を大きく劣化させることなく、各回の計測時間を短縮化する計測方法を提供することを課題とする。また、このような計測方法において、さらに、計測結果の方向依存性を排除することを課題とする。
【0012】
【課題を解決するための手段】
上記課題を達成するために本発明のMRI装置は、画像の再構成に必要な3次元のk空間データを取得する計測制御手段と、前記k空間データから前記画像を再構成する信号処理手段とを備え、前記計測制御手段は、前記k空間をその低周波領域を含む必須計測領域と複数の周辺計測領域とに分割すると共に該分割を2つのエンコード方向に行い、前記必須計測領域と1以上の前記周辺計測領域のデータ取得を、該1以上の周辺計測領域の内の少なくとも1つを変えて繰り返し、前記信号処理手段は、前記繰り返し毎に、前記取得された必須計測領域と周辺計測領域のデータを含んで前記画像を再構成するものである。
このようにすれば、 k 空間を ky 方向のみに分割する前記従来の米国特許 5,713,358 号、米国特許 5,830,143 号に比べ、計測結果に生じる方向依存性を低減し、これによって生成される画像品質の低下を軽減することができる。
好ましくは、前記周辺計測領域は、前記k空間における前記必須計測領域以外の高周波領域の少なくとも一部とする。
【0013】
このような計測方法によれば、各回の計測ステップにおける測定ステップで測定されるデータに、必ず一つの周辺計測領域分の高周波領域の成分が含まれることになるため、前記従来の米国特許5,713,358号、米国特許5,830,143号に比べ、高周波数領域の時間分解能の劣化は小さく、したがって、生成画質の劣化も小さい。
【0015】
また、k空間の全領域を計測するのではなく、その一部を非計測とすることができる。例えば、前記必須計測領域および前記複数の周辺計測領域の和集合空間に外接する直方体状のk空間は、前記和集合空間と重複しない、非計測領域を有するものとすることができる。このように、直方体状のk空間のうちの、計測目的上重要でない空間周波数部分を非計測領域とすることにより、計測目的への適合度を大きく損なうことなく、時間分解能を向上することができる。この場合、前記非計測領域にはゼロを埋めて前記再構成を行っても良い。
【0016】
さらに、前記繰り返しの際に、前記分割を変更しても良い。例えば、前記繰り返しの間に、前記複数の周辺計測領域、または、前記必須計測領域および前記複数の周辺計測領域の設定を変更してもよい。このようにすれば、時間経過に伴い、空間分解能と時間分解能の重要性の程度が変化していくような場合でも、それに対応した計測を行うことができるようになる。例えば、造影剤が注入された被検体の撮像の場合に、前記造影剤が撮像領域に滞在する期間では前記周辺計測領域を狭くし、前記造影剤が撮像領域に滞在していない期間で前記周辺計測領域を広くして良い。
【0017】
さらに、、データ取得されなかった前記周辺計測領域のデータとして、同じ領域において最近の過去に取得されたデータを用いてもよい。
【0018】
【発明の実施の形態】
以下、本発明の実施形態について、ダイナミックMRAへの適用を例にとり説明する。
【0019】
図1に、本実施形態に係るMRI装置の構成を示す。
図示するように、本MRI装置は、NMR現象を利用して被検体の断層像を得るもので、静磁場発生磁石2と、傾斜磁場発生系3と、シーケンサ4と、送信系5と、受信系6と、信号処理系7と、中央処理装置(CPU)8とを備えて成る。
【0020】
静磁場発生磁石2は、被検体1の周りにその体軸方向または体軸と直交する方向に均一な静磁場を発生させるもので、被検体1の周りのある広がりをもった空間に永久磁石方式または常電導方式あるいは超電導方式の磁場発生手段が配置されている。
【0021】
傾斜磁場発生系3は、X,Y,Zの三軸方向に巻かれた傾斜磁場コイル9と、それぞれの傾斜磁場コイルを駆動する傾斜磁場電源10とから成り、後述のシーケンサ4からの命令に従ってそれぞれのコイルの傾斜磁場電源10を駆動することにより、X,Y,Zの三軸方向の傾斜磁場Gx,Gy,Gzを被検体1に印加するようになっている。この傾斜磁場の加え方により被検体1の特定のスライス又はスラブを選択的に励起することができ、また計測空間(k空間)における計測点(サンプリング点)の位置、計測順序を規定することができる。
【0022】
シーケンサ4は、CPU8の制御で動作し、被検体1の断層像のデータ収集に必要な種々の命令を、傾斜磁場発生系3、送信系5及び受信系6に送るようになっている。シーケンサ4が制御する傾斜磁場発生系3、送信系5及び受信系6の動作タイミングはパルスシーケンスと呼ばれ、ここではパルスシーケンスの一つとして三次元血流撮像のためのシーケンスが採用される。シーケンサ4の制御については後に詳述する。
【0023】
送信系5は、シーケンサ4から送り出される高周波パルスにより被検体1の生体組織を構成する原子の原子核に核磁気共鳴を起こさせるために高周波磁場を照射するもので、高周波発振器11と変調器12と高周波増幅器13と送信側の高周波コイル14aとから成る。高周波発振器11から出力された高周波パルスをシーケンサ4の命令にしたがって変調器12で振幅変調し、この振幅変調された高周波パルスを高周波増幅器13で増幅した後に被検体1に近接して配置された高周波コイル14aに供給することにより、電磁波が被検体1に照射されるようになっている。
【0024】
受信系6は、被検体1の生体組織の原子核の核磁気共鳴により放出されるエコー信号(NMR信号)を検出するもので、受信側の高周波コイル14bと増幅器15と直交位相検波器16とA/D変換器17とから成る。上記送信側の高周波コイル14aから照射された電磁波による被検体1の応答の電磁波(NMR信号)は被検体1に近接して配置された高周波コイル14bで検出される。検出されたエコー信号は、増幅器15及び直交位相検波器16を介してA/D変換器17に入力されディジタル量に変換され、さらにシーケンサ4からの命令によるタイミングで直交位相検波器16によりサンプリングされた二系列の収集データとされ、信号処理系7に送られる。
【0025】
信号処理系7は、CPU8と、磁気ディスク18及び磁気テープ19等の記録装置と、CRT等のディスプレイ20とから成り、CPU8で3次元フーリエ変換、補正係数計算像再構成等の処理を行い、任意断面の信号強度分布あるいは複数の信号に適当な演算を行って得られた分布を画像化してディスプレイ20に断層像として表示するようになっている。なお、図1において、送信側及び受信側の高周波コイル14a、14bと傾斜磁場コイル9は、被検体1の周りの空間に配置された静磁場発生磁石2の磁場空間内に設置されている。
【0026】
次に、このようなMRI装置における、ダイナミックMRA撮像動作について説明する。
まず、計測に用いるパルスシーケンスについて説明する。
【0027】
シーケンサ4は、所定のパルスシーケンス、ここでは3次元MRAシーケンスに従い傾斜磁場発生系3、送信系5及び受信系6の動作タイミングを制御する。このパルスシーケンスはCPU8に備えられたメモリに予めプログラムとして組み込まれており、他のパルスシーケンスと同様、使用者が撮影の目的に応じて適宜選択することにより実行することができる。即ち、CPU8の入力装置を介して造影剤を用いたMRAが選択されると、シーケンス4はCPU8によって制御され、3次元MRAシーケンスを実行する。
【0028】
このパルスシーケンスは、例えば図2に示すように、グラディエントエコー法を基本とするシーケンスで、3次元MRAシーケンスに一般的なものである。ただし、Gradient Moment Nulling 等のフローコンペンセーションなどが付加される場合もある。
【0029】
図示した例では、高周波磁場パルスRFを印加して目的血管を含む領域(スラブ)を励起した後、Z方向の位相エンコード傾斜磁場GzおよびY方向の位相エンコード傾斜磁場Gyを印加し、次いで読み出し/周波数エンコード傾斜磁場Gx印加するとともにその極性を反転させてエコー信号を計測している。3次元イメージングでは、このようなパルスシーケンスを、Z方向の位相エンコード傾斜磁場Gz、Y方向の位相エンコード傾斜磁場Gyの磁場強度の組み合わせを変えながら、所定の繰り返し時間TRで繰り返すことによって、三次元データを得る。
【0030】
ここで、Z、Y方向のエンコードステップ数はZ、Y方向の領域の分割数を決めるもので、例えばそれぞれ128、256などに設定されている。また、このZ、Y方向のエンコードステップ数によってk空間のky−kz平面の大きさが規定される。即ち図2のシーケンスにおいて、傾斜磁場強度Gzのある値と傾斜磁場強度Gyのある値の組み合わせのときに計測された信号は、Gz、Gyの値に対応するk空間の座標(ky,kz)に配置されることになる。逆に言えば、k空間のky−kz平面のある座標(ky1、kz1)をサンプルするということは、ky1に対応するGyとkz1に対応するGzの組み合わせにおいて、図2のパルスシーケンスを実行することであり、座標(ky1、kz1)をサンプルしないということは、ky1に対応するGyとkz1に対応するGzの組み合わせにおいて、図2のパルスシーケンスの実行を省略するということである。
【0031】
従来の基本的なダイナミックMRAでは、各回の計測におけるZ、Y方向のエンコードステップ数、ステップ幅を固定として、毎回の計測において、全てのGyとGzの組み合わせについて図2のパルスシーケンスを実行し、全ての座標にデータが充足された3次元のk空間を得る。
【0032】
これに対し、本実施形態では、ダイナミックMRA撮像において、各回の計測におけるZ、Y方向のエンコードステップ数、ステップ幅を固定とせずに、各回の計測を以下のように制御する。
【0033】
すなわち、本実施形態では、予め可能なZ、Y方向の位相エンコードステップ数に応じて大きさを定めたky, kzの座標軸を有する2次元のk空間上に、その中央領域である必須計測領域と、その周辺に複数の周辺計測領域を設定する。ここで、各周辺計測領域は他の周辺計測領域と部分的に同じ領域を含んでいてもよい。必須計測領域については、計測の各回において必ず計測し、各周辺計測領域については計測の各回において所望の順序、例えば巡回的に一つづつ計測する。また各回の計測において計測しなかったk空間の領域については、その回で計測しなかったk空間の部分を計測した回であって、その回に最も近い過去の回の計測結果を流用する。そして、前述のように各回の計測毎に、その回の計測と過去の回の計測結果の流用によって必須計測領域と全ての周辺計測領域のデータが充填されたk空間に対して3次元フーリエ変換を施し、3次元データを得、これを投影した2次元画像を得る。
【0034】
或いはこの再構成した3次元データと、造影剤投入前に計測しておいて3次元データとの差分をとり、これを最終的な3次元データとしてもよい。この差分は、複素差分であることが好ましいが、強度絶対値の差分であっても良い。また、各回の全領域充填後のk空間データと造影剤投入前に計測しておいたk空間データとのk空間データ同士で複素差分を取った後に再構成を行うようにしてもよい。
【0035】
さらに3次元データは、これを所定のまたは指定された配置関係において、仮想的な2次元スクリーンと共に仮想3次元座標系に配置し、3次元データを仮想的な2次元スクリーンに投影した2次元画像を生成し、これを表示する。ここで、投影の方向は、任意で良いが、一般的には、冠状断、矢状断、軸横断の方向に投影する。さらに、ある軸を中心として投影方向を回転させながら順次二次元画像を生成し、これらをつなげて動画像を生成するようにしてもよい。たとえば、3次元の血管データDbを、A軸の回りに45度ずつ回転させて得られる3方向に対して、3次元の血管データDbを投影して得られる3枚の二次元画像Ixは、図3に示すようになる。
【0036】
なお、この2次元画像の生成には、一般には、2次元スクリーンのある座標に投影される位置にあるボクセルの内の最大値のボクセルのみをその座標に投影する最大値投影法(Maximum Intensity Projection)とRay tracing (光線追跡法)の組み合わせの他、Surface rendering、Volume renderingなどの周知のレンダリング手法が用いられる。
【0037】
次に、このような本発明の実施形態において、周辺計測領域の設定および必須計測領域と周辺計測領域の計測の具体例について説明する。
【0038】
まず、第1の具体例について示す。
本例では、図4aに示すように、ky方向に、3次元のk空間を3分割し、図中上下の領域を周辺計測領域42、中央の領域を必須計測領域41とする。そして、図4bに示したように必須計測領域41については、計測の各回において必ず計測し、各周辺計測領域42については計測の各回において交互に一つづつ計測する。
【0039】
また各回の計測において計測しなかった3次元のk空間の領域については、その回で計測しなかった3次元のk空間の部分を計測した回であって、その回に最も近い過去の回の計測結果を流用する。画像再構成に際しては、前述のように各回の計測毎に、その回の計測と過去の回の計測結果の流用によって充填された3次元のk空間に対して3次元フーリエ変換を施し、3次元データを得、これを投影した2次元画像を得る。
【0040】
なお、図4b中、黒塗りの領域は、その回の計測で計測を行わない領域である。以下の各図についても同様である。
【0041】
本第1具体例を用いる場合の、計測の手順の一例を図5に示す。
すなわち、図5aに示す静脈への造影剤注入前から、計測を開始し、図5b、cに示すように、計測の各回において3次元のk空間の必須計測領域については必ず計測し、3次元のk空間の各周辺計測領域については計測の各回において交互に一つづつ計測するように、計測を繰り返す。各回の計測において計測しなかったk空間の領域については、その回で計測しなかったk空間の部分を計測した回であって、その回に最も近い過去の回の計測結果を流用し、3次元のk空間を充填する。充填されたk空間データに対して3次元フーリエ変換を施し、図5dに示すように3次元画像データ(時系列造影剤分布画像)を得、これを投影した図5eの2次元画像(投影血管画像)を得る。
【0042】
この結果、造影剤の血管A、血管B、血管Cへの流動に伴い、順次、各血管が高コントラストで表示されることになる。
【0043】
なお、各回の計測における3次元のk空間の測定の走査方向、すなわち、Z,Y方向の位相エンコーディングの順序は、図4に示す例では、図中k空間の内部に矢印で示したように、必須計測領域のうち、その回で計測されない周辺計測領域の側から始めて、その回で計測される周辺計測領域の側へと走査するようにしているが、これに限られるものではない。
【0044】
以上、説明した第1の具体例の計測によれば、各回の計測に、必ず一つの周辺計測領域分の高周波領域の成分が含まれることになるため、高周波領域を低頻度で計測する従来技術(例えば米国特許5,713,358号、米国特許5,830,143号)に比べ、高周波数領域の時間分解能の劣化は小さく、したがって、生成画質の劣化も小さい。
【0045】
次に、第2の具体例について説明する。
本例では、図6aに示すように、ky方向に3次元のk空間を3分割し図中上下の領域を周辺計測領域62とすると共に、kz方向に3次元のk空間を3分割し図中左右の領域を周辺計測領域63とし、これら4つの周辺計測領域に囲まれた中央の領域を必須計測領域61とする。図6bに示したように必須計測領域61については、計測の各回において必ず計測し、各周辺計測領域62、63については計測の各回において、所望の順序で一つづつ計測する。周辺計測領域を計測する順序は、例えば図6bに示したように巡回的に行ってもよいし、上下(或いは左右)の周辺計測領域を順次計測し、次いで左右(或いは上下)の周辺計測領域を順次計測してもよい。
【0046】
各回の計測において計測しなかった3次元のk空間の領域については、前記第1の具体例と同様に、その回で計測しなかった3次元のk空間の部分を計測した回であって、その回に最も近い過去の回の計測結果を流用し、各回の計測毎に、その回の計測と過去の回の計測結果の流用によって充填された3次元のk空間データに対して3次元フーリエ変換を施し、3次元画像データを得、これを投影した2次元画像を得る。
【0047】
なお、本例においても、各回の計測における3次元のk空間の測定の走査方向、すなわち、Z,Y方向の位相エンコーディングの順序は、図示するように、必須測定領域のうち、その回で計測されない周辺計測領域の側から始め、その回で計測される周辺計測領域の側へと走査するようにしてもよいが、これに限られるものではない。
【0048】
以上のように、本例によれば、必須計測領域、周辺計測領域の境界をkyとkzの2方向について設定するために、さらに、k空間をky方向のみに分割する従来技術(例えば米国特許5,713,358号、米国特許5,830,143号)に比べ、計測結果に生じる方向依存性を低減し、これによって生成される画像の品質を劣化を軽減することができる。
【0049】
本例では、図7aに示すように、ky方向に3次元のk空間を3分割し、kz方向に3次元のk空間を3分割して得られる9つの領域のうち、図中上側左右方向中央、下側左右方向中央、左側上下方向中央、右側上下方向中央の4領域を周辺計測領域72とすると共に、周辺計測領域に囲まれた中央の領域を必須計測領域71とする。また、残りの領域を非計測領域73として計測を行わない領域とする。この具体例でも、図7bに示したように必須計測領域については、計測の各回において必ず計測し、各周辺計測領域については計測の各回において所望の順序で一つづつ計測する。
【0050】
各回の計測において計測しなかった3次元のk空間の領域については、その回で計測しなかった3次元のk空間の部分を計測した回であって、その回に最も近い過去の回の計測結果を流用し、各回の計測毎に、その回の計測と過去の回の計測結果の流用によって非計測領域以外の領域が充填された3次元のk空間データを得る。非計測領域についてゼロフィリングした後、k空間データに対して3次元フーリエ変換を施し、3次元画像データを得、これを投影した2次元画像を得る。
【0051】
なお、本例においても、各回の計測における3次元のk空間の測定の走査方向、すなわち、Z,Y方向の位相エンコーディングの順序は、図示するように、必須計測領域のうち、その回で計測される周辺計測領域と反対の側から始め、その回で計測される周辺計測領域の側へと走査するようにしてもよいが、これに限定されない。
【0052】
以上のように、本例によれば、さらに、診断等に比較的重要でない3次元のk空間上の角の領域を非計測領域とすることにより、より短い時間で計測を繰り返せるようになるので、比較的小さな画質の低下において時間分解能を向上することができる。
【0053】
次に第4の具体例について説明する。
本例では、図8aに示すように、3次元k空間のうちky−kz平面に内接する円の内側の領域と、外側の領域に2分割し、外側の領域を非計測領域83とする。さらに円の内側の領域については、それに包含される四角形の内側の領域と外側の領域とに分割し、四角形の内側の領域を必須計測領域81とする。非計測領域83と必須計測領域81との間の領域、すなわち円の内側であって必須計測領域81の外側の領域については、これを4つに等分割し、得られる4領域を周辺計測領域82とする。なお、ky−kz平面が正方形の場合、円は真円となるが、ky軸とkz軸の長さが異なり長方形の場合には、円は楕円状となる。本明細書におい「円」はこのような「楕円」を含む意である。
【0054】
この場合にも、図8bに示したように必須計測領域81については、計測の各回において必ず計測し、各周辺計測領域82については計測の各回において所望の順序で一つづつ計測する。
【0055】
各回の計測において計測しなかった3次元のk空間の領域については、第1の具体例同様に、その回で計測しなかった3次元のk空間の部分を計測した回であって、その回に最も近い過去の回の計測結果を流用し、各回の計測毎に、その回の計測と過去の回の計測結果の流用によって非計測領域以外の領域が充填された3次元のk空間データを得る。ここでも非計測領域はゼロフィリングしてもよい。このk空間データに3次元フーリエ変換を施し、3次元画像データを得、これを投影した2次元画像を得る。
【0056】
なお、本例においても、各回の計測における3次元のk空間の測定の走査方向、すなわち、Z,Y方向の位相エンコーディングの順序は、図示するように、必須計測領域のうち、その回で計測される周辺計測領域と反対の側から始め、その回で計測される周辺計測領域の側へと走査するようにしてもよいが、これに限定されない。
【0057】
または、上述するようなラスター走査ではなく、図9に示すように、螺旋状に走査してもよい。この場合には、ky−kz平面の原点との距離が近いものより走査する。或いは、同じ螺旋状の走査であっても、図10に示すように、必須計測領域のうち、ky−kz平面の原点に達するまでは、原点との距離が遠いものより走査を始め、当該原点に達した後は、ky−kz平面の原点との距離が近いものより走査し、次いで周辺計測領域を走査するようにしてもよい。
【0058】
次に第5の具体例について説明する。
本例は、図11aに示すように、ky−kz平面に設定した円内の領域を必須計測領域111とし、さらにky−kz平面に少なくとも必須計測領域を包含するように円を設定し、二つの円で囲まれる領域を周方向に4つに等分割して得られる4領域を周辺計測領域112としたものである。周辺計測領域112の外側は非計測領域113である。
【0059】
図11bに示したように必須計測領域111については、計測の各回において必ず計測し、各周辺計測領域112については計測の各回において所望の順序で一つづつ計測する。
【0060】
第1の具体例同様に、各回の計測において計測しなかった3次元のk空間の領域については、その回で計測しなかった3次元のk空間の部分を計測した回であって、その回に最も近い過去の回の計測結果を流用し、各回の計測毎に、その回の計測と過去の回の計測結果の流用によって非計測領域以外の領域が充填された3次元のk空間データを得る。ここでも非計測領域113はゼロフィリングする。このk空間データに3次元フーリエ変換を施し、3次元画像データを得、これを投影した2次元画像を得る。
【0061】
なお、本例においても、各回の計測における3次元のk空間の測定の走査方向、すなわち、Z,Y方向の位相エンコーディングの順序は任意である。すなわち、図11に示すようにky−kz平面の原点との距離が近いものより螺旋状に走査するようにしてもよいし、図12に示すように、必須計測領域のうち、ky−kz平面の原点に達するまでは、当該原点との距離が遠いものより走査し、当該原点に達した後は、ky−kz平面の原点との距離が近いものより走査し、次いで周辺計測領域を走査するようにしてもよい。
【0062】
以下、第6の具体例を説明する。
本例では、図13に示すように、一連の計測手順の流れの中で、周辺計測領域132の大きさを変更する。
【0063】
すなわち、より時間分解能が必要とされる期間は個々の周辺計測領域を小さくとることにより、各回の計測で計測する領域の大きさを小さくして、より短い時間で計測を繰り返せるようにし、より空間分解能が必要とされる期間は個々の周辺計測領域を大きくとることにより、各回の計測で計測する領域の大きさを大きくして、より高い空間分解能で計測を行えるようにする。
【0064】
具体的には、たとえば、図13aに示すように造影剤の投入前後及び主要な大動脈へ造影剤が到達する時間帯には図13b、cに示すように高い時間分解能で計測が行えるように周辺計測領域132を小さく設定する。また図13aに示すように造影剤が動脈相を過ぎ静脈相に到達した時点では、図13b、cに示すように高い空間分解能で計測が行えるように周辺計測領域を大きく設定することにより、血流動態の変化に合わせた計測を可能とする。
【0065】
なお、図13に示した例では、ky−kz平面に設定した円内の領域を必須計測領域131とし、この必須計測領域131を包含するように設定した円と必須計測領域131との間の領域を周方向にnに等分割して得られるn個の領域を周辺計測領域132とすることとし、nの数を測定開始からの経過時間に応じて変化させることにより個々の周辺計測領域の大きさを変化させている。周辺計測領域の大きさを変化させる方法は、これに限らず、例えば図4に示す周辺計測領域の設定と、図6や図7に示す周辺計測領域の設定とを適宜組み合わせてもよい。
【0066】
以上、本発明の実施形態を説明したが、上述の実施形態における周辺計測領域の数は例示であって、その用途、目的に併せて任意に設定できる。
【0067】
また各回の計測において計測しなかった3次元のk空間の領域については、その回で計測しなかった3次元のk空間の部分を計測した回であって、その回に最も近い過去の回の計測結果を流用したが、他の手法で、たとえば、その回で計測しなかった3次元のk空間の部分を計測した回であって、その回に最も近い過去と次回の2つの計測結果の補間によって、その回で計測しなかった3次元のk空間の部分のデータを生成するようにしてもよい。
【0068】
また、以上では、Z方向とY方向について位相エンコーディングを行うことにより3次元計測を行う場合について説明したが、Z方向についてはスライス選択傾斜磁場とRFパルスの組み合わせによって、その位置をエンコードする3次元計測を行う場合についても、以上の必須計測領域および周辺計測領域を用いた計測は同様に適用することができる。
【0069】
また、以上では3次元計測を行う場合について説明したが、2次元計測(1スライス毎の計測)を行う場合についてもkx−ky方向について必須計測領域および周辺計測領域の境界を設定することにより同様に適用することができる、また、以上では、造影剤と用いたダイナミックMRAへの適用を例にとり説明したが、以上の必須計測領域および周辺計測領域を用いた計測は、MRIにおいて任意の繰り返し計測を行う場合について同様に適用することができる。
【0070】
【発明の効果】
以上のように、本発明によれば、MRIにおいて繰り返し撮像を行う際に、高周波数領域の時間分解能を大きく劣化させることなく、各回の計測時間を短縮化する計測方法を提供することができる。また、このような計測方法において、さらに、計測結果の方向依存性を排除することができる。
【図面の簡単な説明】
【図1】本発明の実施形態に係る磁気共鳴イメージング装置の構成を示すブロック図である。
【図2】本発明の実施形態に係る計測に用いることのできるパルスシーケンスの例を示すタイミングチャートである。
【図3】本発明の実施形態に係る3次元データを投影した二次元画像の生成法の例を示す図である。
【図4】本発明の実施形態に係る計測法の第1の具体例を示す図である。
【図5】本発明の実施形態に係る計測法の第1の具体例による計測手順を示す図である。
【図6】本発明の実施形態に係る計測法の第2の具体例を示す図である。
【図7】本発明の実施形態に係る計測法の第3の具体例を示す図である。
【図8】本発明の実施形態に係る計測法の第4の具体例を示す図である。
【図9】本発明の実施形態に係る計測法の第4の具体例における3次元k空間の走査順序の他の例を示す図である。
【図10】本発明の実施形態に係る計測法の第4の具体例における3次元k空間の走査順序の他の例を示す図である。
【図11】本発明の実施形態に係る計測法の第5の具体例を示す図である。
【図12】本発明の実施形態に係る計測法の第5の具体例における3次元k空間の走査順序の他の例を示す図である。
【図13】本発明の実施形態に係る計測法の第6の具体例における計測手順を示す図である。
【符号の説明】
1…被検体
2…磁場発生装置
3…傾斜磁場発生系
4…シーケンサ
5…送信系
6…受信系
7…信号処理系
8…CPU
9…傾斜磁場コイル
14a…送信側の高周波コイル
14b…受信側の高周波コイル
Claims (8)
- 画像の再構成に必要な3次元のk空間データを取得する計測制御手段と、前記k空間データから前記画像を再構成する信号処理手段とを備えた磁気共鳴イメージング装置において、
前記計測制御手段は、前記k空間をその低周波領域を含む必須計測領域と複数の周辺計測領域とに分割すると共に該分割を2つのエンコード方向に行い、前記必須計測領域と1以上の前記周辺計測領域のデータ取得を、該1以上の周辺計測領域の内の少なくとも1つを変えて繰り返し、
前記信号処理手段は、前記繰り返し毎に、前記取得された必須計測領域と周辺計測領域のデータを含んで前記画像を再構成することを特徴とする磁気共鳴イメージング装置。 - 請求項1記載の磁気共鳴イメージング装置において、
前記周辺計測領域は、前記k空間における前記必須計測領域以外の高周波領域の少なくとも一部であることを特徴とする磁気共鳴イメージング装置。 - 請求項1又は2記載の磁気共鳴イメージング装置において、
前記計測制御手段は、データ取得を行う前記周辺計測領域を前記必須計測領域の周りに巡回させて変えることを特徴とする磁気共鳴イメージング装置。 - 請求項1乃至3のいずれかに記載の磁気共鳴イメージング装置において、
前記計測制御手段は、更に非計測領域を少なくとも一つ含むように前記分割を行い、
前記信号処理手段は、前記非計測領域にはゼロを埋めて前記再構成を行うことを特徴とする磁気共鳴イメージング装置。 - 請求項1乃至4のいずれかに記載の磁気共鳴イメージング装置において、
前記計測制御手段は、前記繰り返しの際に、前記分割を変更することを特徴とする磁気共鳴イメージング装置。 - 請求項5記載の磁気共鳴イメージング装置において、
前記計測制御手段は、前記周辺計測領域を広狭させて、前記画像の空間分解能又は時間分解能を変更することを特徴とする磁気共鳴イメージング装置。 - 請求項6記載の磁気共鳴イメージング装置において、
前記計測制御手段は、造影剤が注入された被検体の撮像の場合に、前記造影剤が撮像領域に滞在する期間では前記周辺計測領域を狭くし、前記造影剤が撮像領域に滞在していない期間で前記周辺計測領域を広くすることを特徴とする磁気共鳴イメージング装置。 - 請求項1乃至7のいずれかに記載の磁気共鳴イメージング装置において、
前記信号処理手段は、データ取得されなかった前記周辺計測領域のデータとして、同じ領域において最近の過去に取得されたデータを用いることを特徴とする磁気共鳴イメージング装置。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000381171A JP3901448B2 (ja) | 2000-12-15 | 2000-12-15 | 磁気共鳴イメージング装置 |
PCT/JP2001/010538 WO2002045584A1 (fr) | 2000-12-04 | 2001-12-03 | Procede de mesure dans un dispositif d'imagerie par resonance magnetique et dispositif d'imagerie par resonance magnetique |
US10/433,736 US7167740B2 (en) | 2000-12-04 | 2001-12-03 | Measuring method in magnetic resonance imaging device and magnetic resonance imaging device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000381171A JP3901448B2 (ja) | 2000-12-15 | 2000-12-15 | 磁気共鳴イメージング装置 |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2002177240A JP2002177240A (ja) | 2002-06-25 |
JP2002177240A5 JP2002177240A5 (ja) | 2005-10-20 |
JP3901448B2 true JP3901448B2 (ja) | 2007-04-04 |
Family
ID=18849227
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2000381171A Expired - Fee Related JP3901448B2 (ja) | 2000-12-04 | 2000-12-15 | 磁気共鳴イメージング装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3901448B2 (ja) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4699739B2 (ja) * | 2004-10-28 | 2011-06-15 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Mri装置 |
US10098563B2 (en) | 2006-11-22 | 2018-10-16 | Toshiba Medical Systems Corporation | Magnetic resonance imaging apparatus |
US9538936B2 (en) | 2006-11-22 | 2017-01-10 | Toshiba Medical Systems Corporation | MRI apparatus acquires first and second MR data and generates therefrom third image data having higher contrast between blood and background tissues |
WO2012085810A2 (en) * | 2010-12-22 | 2012-06-28 | Koninklijke Philips Electronics N.V. | Rapid parallel reconstruction for arbitrary k-space trajectories |
JP5526168B2 (ja) * | 2012-02-20 | 2014-06-18 | 株式会社東芝 | 磁気共鳴イメージング装置 |
-
2000
- 2000-12-15 JP JP2000381171A patent/JP3901448B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2002177240A (ja) | 2002-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4646015B2 (ja) | 磁気共鳴イメージング(mri)システムの作動方法 | |
US7999543B2 (en) | MR method for spatially-resolved determination of relaxation parameters | |
JP3847512B2 (ja) | 磁気共鳴イメージング装置 | |
US7450982B2 (en) | Magnetic resonance imaging system and method | |
JP3386864B2 (ja) | 核磁気共鳴撮影方法及び装置 | |
JP4416221B2 (ja) | 磁気共鳴画像診断装置 | |
JP2014502910A (ja) | インターリーブスピンロッキングイメージング | |
US7167740B2 (en) | Measuring method in magnetic resonance imaging device and magnetic resonance imaging device | |
JPH06169896A (ja) | 磁気共鳴イメージング装置 | |
JP3847554B2 (ja) | 磁気共鳴イメージング装置 | |
JP3516421B2 (ja) | Mri装置 | |
JP3901448B2 (ja) | 磁気共鳴イメージング装置 | |
JP3514547B2 (ja) | 磁気共鳴イメージング装置 | |
JP3452400B2 (ja) | 磁気共鳴イメージング装置 | |
JP3847519B2 (ja) | 磁気共鳴イメージング装置 | |
JP3501182B2 (ja) | 流速画像を算出できる磁気共鳴イメージング装置 | |
JP2019122554A (ja) | 撮影装置 | |
JPH0576518A (ja) | 磁気共鳴イメージング装置 | |
JP2001252262A (ja) | Mri差分画像処理方法及びmri装置 | |
JP2021137501A (ja) | 磁気共鳴イメージング方法及び装置 | |
JP4297541B2 (ja) | 磁気共鳴イメージング装置 | |
JP3903276B2 (ja) | 磁気共鳴撮像装置 | |
JP2000316830A (ja) | 磁気共鳴イメージング方法及びそれを用いた磁気共鳴イメージング装置 | |
JP3440134B2 (ja) | 磁気共鳴イメージング装置 | |
JP7383509B2 (ja) | 医用画像データ処理装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050610 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050610 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050613 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20060919 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20061115 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20061222 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20061226 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100112 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110112 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110112 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120112 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130112 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140112 Year of fee payment: 7 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
LAPS | Cancellation because of no payment of annual fees |