[go: up one dir, main page]

JP3276669B2 - 磁気共鳴イメージング装置 - Google Patents

磁気共鳴イメージング装置

Info

Publication number
JP3276669B2
JP3276669B2 JP13854092A JP13854092A JP3276669B2 JP 3276669 B2 JP3276669 B2 JP 3276669B2 JP 13854092 A JP13854092 A JP 13854092A JP 13854092 A JP13854092 A JP 13854092A JP 3276669 B2 JP3276669 B2 JP 3276669B2
Authority
JP
Japan
Prior art keywords
echo
magnetic field
space
data
gradient magnetic
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
Application number
JP13854092A
Other languages
English (en)
Other versions
JPH05329127A (ja
Inventor
孝治 梶山
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Healthcare Manufacturing Ltd
Original Assignee
Hitachi Medical Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hitachi Medical Corp filed Critical Hitachi Medical Corp
Priority to JP13854092A priority Critical patent/JP3276669B2/ja
Priority to US08/067,624 priority patent/US5357200A/en
Publication of JPH05329127A publication Critical patent/JPH05329127A/ja
Priority to US08/274,281 priority patent/US5543709A/en
Application granted granted Critical
Publication of JP3276669B2 publication Critical patent/JP3276669B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Description

【発明の詳細な説明】
【0001】
【産業上の利用分野】本発明は、磁気共鳴イメージング
装置に係り、特に、いわゆる高速スピンエコー法を用い
た磁気共鳴イメージング装置に関する。
【0002】
【従来の技術】磁気共鳴イメージング装置は、磁気共鳴
現象を利用して被検体中の所望の検査部位における原子
核スピン(以下、単にスピンと称す。)の密度分布、緩
和時間分布等を計測して、その計測データから、被検体
の断面を画像表示するものである。
【0003】この装置では、図9に示すように、0.0
2〜2テスラ程度の静磁場を発生させる静磁場発生装置
4の中に被検体7が置かれる。この時、被検体中のスピ
ンは静磁場の強さH0によって決まる周波数で静磁場の
方向を軸として歳差運動を行なうようになる。この周波
数はラーモア周波数と称されている。ここで、ラーモア
周波数ν0は、
【0004】
【数1】ν0=γ/2π・H0 (1) で表わせる。ここで、γは磁気回転比で原子核の種類毎
に固有の値を持つ。また、ラーモア歳差運動の角速度を
ω0とすると、
【0005】
【数2】ω0=2πν0 (2) の関係があるため、
【0006】
【数3】ω0=γH0 (3) で与えられる。
【0007】ここで、高周波照射コイル11によって計
測しようとする原子核のラーモア周波数ν0に等しい周
波数の高周波磁場(電磁波)を加えると、スピンが励起
され高いエネルギー状態に遷移する。この高周波磁場を
打ち切ると、スピンはもとの低いエネルギー状態に戻る
ようになる。そして、この際に放出される電磁波を高周
波受信コイル14で受信し、増幅器15で増幅および波
形整形した後、A/D変換器17(以下、ADCと称
す)でディジタル化して中央処理装置1(以下、CPU
と称す)に送る。CPU1では、このデータを基に再構
成演算し、この演算されたデータが被検体7の断層画像
としてディスプレイ18に表示される。上記の高周波磁
場は、CPU1により制御されるシーケンサ2が送り出
す信号を高周波送信コイル用増幅器10によって増幅し
たものを高周波送信コイル11に送ることで得られる。
【0008】MRI装置においては、以上の静磁場と高
周波磁場の他に、空間内の位置情報を得るための傾斜磁
場を作るために傾斜磁場コイル群21を備えている。こ
れらの傾斜磁場コイル13は、シーケンサ2からの信号
で動作する傾斜磁場コイル用電源12から電流が供給さ
れ、傾斜磁場を発生するようになっている。
【0009】ここで、MRI装置の撮影原理を説明す
る。図10(a)に示すようにZ方向の静磁場H0中に
置かれた原子核は、古典物理学的に見ると1個の棒磁石
のように振舞い、先に述べたラーモア周波数ν0でZ軸
の回りに歳差運動を行なっている。この周波数は前記
(2)式で与えられ静磁場の強度に比例するものとなっ
ている。(1)式及び(3)式におけるγは磁気回転比
と呼ばれ原子核に固有の値を持っている。一般には測定
対象の原子核は膨大な数にのぼり、それぞれが勝手な位
相で回転しているために、全体で見るとX−Y面内の成
分は打ち消しあい、Z方向成分のみの巨視的磁化が残る
ようになる。この状態でX方向にラーモア周波数ν0
等しい周波数の高周波磁場H1を印加すると(図10
(b))、巨視的磁化はY方向に倒れ始める。この倒れ
る角度はH1の振幅と印加時間との積にほぼ比例し、パ
ルス印加時点に対し90°倒れる時のH1を90°パル
ス、180°倒れるときのH1を180°パルスと称し
ている。
【0010】ここで、現在MRI装置による撮影で一般
的に用いられている方法に2次元フーリエイメージング
法がある。この方法のうち代表的なスピンエコー法の模
式的なパルスシーケンスを図11に示す。このパルスシ
ーケンスでは、まず、90°パルス28を印加した後、
エコー時間をTeとしたときTe/2の時間後に180
°パルス29を加える。90°パルス28を加えた後、
各スピンはそれぞれに固有の速度でX−Y面内で回転を
始めるため、時間の経過とともに各スピン間に位相差が
生じる。ここで180°パルス29が加わると、各スピ
ンはx’軸に対称に反転し、その後も同じ速度で回転を
続けるために時刻Teでスピンは集束しエコー信号を形
成する。
【0011】さらに図12にマルチエコー計測のパルス
シーケンスの模式図を示す。マルチエコー計測では、さ
らに180°パルスを印加し時刻2Te,3Te,・・で
もスピンを収束させエコー信号を計測するようになって
いる。
【0012】そして、このようにして信号は計測される
が、断層画像を構成するためには信号の空間的な分布を
求めねばならない。このために線形な傾斜磁場を用いて
いる。均一な静磁場に傾斜磁場を重畳することで空間的
な磁場勾配ができるようになる。先にも述べたようにス
ピンの回転周波数は磁場強度に比例しているから傾斜磁
場が加わった状態においては、各スピンの回転周波数は
空間的に異なる。従って、この周波数を調べることによ
って各スピンの位置を知ることができる。この目的のた
めに、位相エンコード方向傾斜磁場31、周波数エンコ
ード方向傾斜磁場33,34が用いられている。
【0013】以上に述べたパルスシーケンスを基本単位
として、位相エンコード方向傾斜磁場の強度を毎回変え
ながら一定の繰り返し時間(TR)毎に、所定回数、例
えば256回繰り返す。こうして得られた計測信号を2
次元逆フーリエ変換することで巨視的磁化の空間的分布
が求められる。なお、このようなMRI基本原理に関し
ては、「NMR医学」(基礎と臨床)(核磁気共鳴医学
研究会編,丸善株式会社,昭和59年1月20日発行)
に詳述されている。
【0014】次にいわゆる高速スピンエコー法について
記す。高速スピンエコー法は、上述のマルチエコー計測
の各エコーを生データのK空間上に振り分けることによ
り、撮像時間を短縮する手法である。図13は高速スピ
ンエコーのパルスシーケンスを示す図である。同図にお
いて、高速スピンエコー法は、マルチエコー計測のよう
に180度パルスを複数回かけることにより行う。そし
て、各エコー信号計測毎に位相エンコード方向に所定の
傾斜磁場を所定の時間印加することにより、各エコーの
信号を生データ空間上の位相方向に振り分けるようにな
っている。振り分け方は、使用するエコー数、及び画質
を決定するエコー番号により決められる。その決め方
は、K空間の直流部分(中央付近)を所望のエコー時間
のデータとし、その隣接する部分は、隣接するエコー番
号のデータとしている。また、K空間上のky方向の端
部のエコーと他の端部のエコー番号が同一または隣接す
るものとしている。
【0015】図14にたとえば第4エコーまで使用し
て、第4エコーのエコー時間の強調を得る時の、K空間
上のデータの並び方を示す。各エコー番号のデータはほ
ぼ同数となっている。まず、低域(中心付近)に所望の
第4エコーのデータ27を割り振る。その隣接する領域
は、第3エコーのデータ26とし、さらに高域(端部)
に向かうに従い、第2エコーのデータ25、第1エコー
のデータ24とする。これにより、ky方向の最高域
(端部)は、第1エコーとなり、また他の端部も第1エ
コーと同じエコー番号となる。所望のエコー番号が他の
エコー時間のエコーの時は、図14をky方向にバレル
シフトさせ、所望のエコー番号が中心に来るようにすれ
ばよい。
【0016】
【発明が解決しようとする課題】このように高速スピン
エコー法を用いた撮像においては、その計測時間をかな
り短くすることができるようになったが、より診断効率
を向上させるためには、さらに、計測時間を短くするこ
とが望まれている。
【0017】それ故、本発明は、このような事情に基づ
ていなされたものであり、その目的とするところのもの
は、さらに計測時間を短くすることのできる磁気共鳴イ
メージング装置を提供することにある。
【0018】
【課題を解決するための手段】本願において開示される
発明の内、代表的なものの概要を簡単に説明すれば、以
下の通りである。 (1)静磁場中の被検体への高周波パルスと傾斜磁場と
を逐次所定値で所定順に印加するパルスシーケンスを複
数回実行し、被検体の所定断面から放出されるエコー信
号からフーリエ変換により画像再構成を行う磁気共鳴イ
メージング装置において、1回のパルスシーケンスで複
数個のエコー信号を計測する計測手段と、前記計測手段
でのエコー信号計測を、位相エンコード方向傾斜磁場及
びエコー信号読みだし傾斜磁場の直前に印加する位相オ
フセット傾斜磁場の値の制御により、K空間を複数のエ
コー領域に分割し任意のエコー番号を前記K空間の中央
部を含むエコー領域に配置すると共に、前記K空間の一
方の端部から前記中央部のエコー領域を越える少なくと
も次のエコー領域までのエコー番号が前記任意のエコー
番号から連続し、前記K空間の一方の端部のエコー番号
と前記次のエコー領域のエコー番号とが連続するように
配置する制御手段と、前記K空間の他方の端部のエコー
番号が前記次のエコー領域のエコー番号から連続するよ
うに推定で埋める推定手段とを備えた。 (2)前述した(1)に記載の磁気共鳴イメージング装
置において、前記パルスシーケンスは、高周波パルスと
して、1回の90゜パルスと該90゜パルスに続く複数
個の180゜パルスとからなる。 (3)静磁場中の被検体への高周波パルスと傾斜磁場と
を逐次所定値で所定順に印加するパルスシーケンスを複
数回実行し、被検体の所定断面から放出されるエコー信
号からフーリエ変換により画像再構成を行う磁気共鳴イ
メージング装置において、単一のパルスシーケンスで複
数個のエコー信号を計測する計測手段と、位相エンコー
ド方向傾斜磁場またはエコー信号読みだし傾斜磁場の直
前に印加する位相オフセット傾斜磁場の値を制御するこ
とにより、前記計測手段でのエコー信号計測を行わせる
制御手段と、前記計測データからK空間上における未計
測データを推定により埋める推定手段とを備え、前記制
御手段は、前記K空間を位相エンコード方向に不均等に
2分割して大となる領域を計測範囲とし、該計測領域に
任意のエコー番号を配置すると共に、前記計測範囲のエ
コー番号が前記任意のエコー番号から連続し、かつ前記
計測範囲の一方の端部のエコー番号が他方の端部のエコ
ー番号に連続するように配置し、前記推定手段は前記K
空間の他方の端部のエコー番号が前記計測範囲の内で前
記K空間の端部となる側のエコー番号から連続するよう
に未計測データを推定する。
【0019】
【作用】前述した(1)及び(2)の手段によれば、本
発明の磁気共鳴イメージング装置は、1回のパルスシー
ケンスを実行することで複数のエコー信号が得られる。
例えば高速スピンエコー法を用いてK空間上の約半数の
データを取り入れる際に、制御手段が任意のエコー番号
をK空間の中央部を含む領域に配置すると共に、K空間
の一方の端部から中央部を越えた計測領域の端部に至る
エコー領域までのエコー番号が任意のエコー番号から連
続するように配置すると共に、K空間の一方の端部のエ
コー番号と中央部を越えた少なくとも次のエコー領域ま
でのエコー番号とが連続するように配置することによっ
て、K空間上ではエコー番号に従った配置となる複数の
エコー信号が得られる。その後、推定手段がいわゆるハ
ーフスキャン再構成法によって該K空間上に全てのデー
タを埋める際に、K空間の他方の端部のエコー番号が計
測範囲の内でK空間の端部となる側のエコー番号から連
続するように未計測データを推定し埋めるようにしたも
のである。この場合、ハーフスキャン再構成法を行う
はそれほどの時間がかからないことから、従来高速スピ
ンエコー法のみで行う場合と比べて、約半分の時間で済
むようになる。したがって、計測時間をさらに短くする
ことができるようになる。また、K空間上の中央部に位
置するエコー番号に対応した画質の画像を得ることがで
きる。
【0020】前述した(3)の手段によれば、本発明の
磁気共鳴イメージング装置は、単一のパルスシーケンス
を実行することで複数のエコー信号が得られる。K空間
を位相エンコード方向に不均等に2分割して大となる領
域を計測範囲とし、該計測範囲に任意のエコー番号を配
置すると共に、計測範囲のエコー番号が任意のエコー番
号から連続し、かつ計測範囲の一方の端部のエコー番号
が他方の端部のエコー番号に連続するように配置し、例
えば高速スピンエコー法を用いて、この計測範囲のデー
タを取り入れる。その後、いわゆるハーフスキャン再構
成法によって、K空間の他方の端部のエコー番号が計測
範囲の内でK空間の端部となる側のエコー番号から連続
するように未計測データを推定し、該K空間上に全ての
データを埋めるようにしたものである。この場合、ハー
フスキャン再構成法を行うにはそれほどの時間がかから
ないことから、従来高速スピンエコー法のみで行う場合
と比べて、約半分の時間で済むようになる。よって、計
測時間をさらに短くすることができるようになる。
【0021】
【0022】
【実施例】第9図は本発明が適用される磁気共鳴イメー
ジング装置の一実施例を示す全体構成のブロック説明図
である。
【0023】同図において、この磁気共鳴イメージング
装置は、大別すると、中央処理装置(CPU)1と、シ
ーケンサ2と、送信系3と、磁場発生系4と、受信系5
と、信号処理系6とを備えて構成されている。
【0024】中央処理装置(CPU)1は、予め定めら
れたプログラムに従って、シーケンサ2、送信系3、受
信系5、信号処理系6の各々を制御するものであり、シ
ーケンサ2は、中央処理装置1からの制御指令に基づい
て動作し、被検体7の断層画像のデータ収集に必要な種
々の命令を送信系3、磁場発生系4の傾斜磁場発生系2
1、受信系5に送るものである。
【0025】送信系3は、高周波発信器8と変調器9と
高周波コイルとしての照射コイル11を有し、シーケン
サ2の指令により高周波発信器8からの高周波パルスを
変調器9で振幅変調し、この振幅変調された高周波パル
スを高周波増幅器10によって増幅して照射コイル11
に供給することにより、所定のパルス状の電磁波を被検
体7に照射するものである。
【0026】磁場発生系4は、被検体7の回りに任意の
方向に均一な静磁場を発生させる静磁場発生磁石を有し
ている。この静磁場発生磁石の内部には、照射コイル1
1の他に、傾斜磁場を発生させる傾斜磁場コイル13
と、受信系5の受信コイル14が設置されている。傾斜
磁場発生系21は互いに直交するデカルト座標軸方向に
それぞれ独立に傾斜磁場を印加できる構成を有する傾斜
磁場コイル13と傾斜磁場コイルに電流を供給する傾斜
磁場電源12と、傾斜磁場電源12を制御するシーケン
サ2により構成されている。
【0027】受信系5は、高周波コイルとしての受信コ
イル14と該受信コイル14に接続された増幅器15と
直交位相検波器16とA/D変換器17とを有し、被検
体7からのNMR信号を受信コイル14が検出すると、
その信号を増幅器15、直交位相検波器16、A/D変
換器17を介しデジタル量に変換するとともに、シーケ
ンサ2からの指令によるタイミングで直交位相検波器1
6によってサンプリングされた二系列の収集データに変
換して中央処理装置1に送るようになっている。
【0028】信号処理系6は、磁気ディスク20、光デ
ィスク19等の外部記憶装置と、CRT等からなるディ
スプレイ18とを有し、受信系5からのデータが中央処
理装置1に入力されると、該中央処理装置1が信号処
理、画像再構成等の処理を実行し、その結果の被検体7
の所望の断面像をディスプレイ18に表示するととも
に、外部記憶装置の磁気ディスク20等に記録するよう
になっている。
【0029】そして、前記シーケンサ2には図2に示す
ようなパルスシーケンスが組み込まれ、このパルスシー
ケンスの繰返しで磁気共鳴シメージング装置が駆動さ
れ、各パルスシーケンス毎にたとえば4つのエコー信号
(第1エコー24、第2エコー25、第3エコー26、
第4エコー27)が検出される。
【0030】同図において、T2では被検体7内のスピ
ンを励起して信号を発生させるための90°パルス28
の印加と、所望の断層面の位置を特定するためのスライ
ス選択傾斜磁場(以下、Gsと称す)301とが同時に
印加される。T4では位相エンコード傾斜磁場(以下、
Gpと称す)31を印加するとともに、周波数エンコー
ド方向傾斜磁場(以下、Gfと称す)33を印加する。
Gpは後に得られる第1エコー信号24をK空間上の所
望の位置に配置するために印加し、Gf33は第1エコ
ー信号24において計測されるスピンの位相が信号収集
の中心で収束するように予め拡散させておくために印加
する。T6では180°パルス291を印加するととも
にGs302を印加する。T8ではGf341を印加しつ
つ、発生する第1エコー信号24を収集する。T12で
は第2エコー信号25を得るために第2の180°パル
ス292を印加するとともに、Gs303を印加する。T
14ではGf342を印加しつつ、発生する第2エコー
信号25を収集する。T16では第3エコー信号26を
得るために第3の180°パルス293を印加するとと
もに、Gs304を印加する。T18ではGf343を印
加しつつ、発生する第3エコー信号26を収集する。そ
して、第4エコー27の場合においても同様にして収集
する。
【0031】このような高速スピンエコー法の特に特徴
的な部分は、位相方向傾斜磁場(Gp)32にあり、位
相エンコード傾斜磁場31と周波数エンコード方向に印
加する読みだし傾斜磁場の直前に印加する位相オフセッ
トの傾斜磁場32とをそれぞれの磁場強度の大きさがK
空間上で隣合う異なったエコーの番号同士が連続するよ
うにするとともに、計測した生データのky方向(位相
エンコード方向)の中央部と一方の端とがエコー番号が
連続するように印加することにある。
【0032】このようなシーケンスの繰返しによって検
出される各エコー信号は計測生データとして、図1
(b)に示すK空間上に格納される。ここで、図1
(b)を用いて計測生データのK空間上の並び方を説明
する。
【0033】同図において、K空間のky方向(位相エ
ンコード方向)は離散値をとり、ky方向の中心はその
中心のエコーの計測までにスピンが感じた位相方向の傾
斜磁場の積分値が零となるようになっている。中心から
ky方向に1つずれた箇所は、スピンが撮像領域(F.
O.V.)の両端で1回転するような位相方向の傾斜磁
場を感じさせるようになっている。式にすると、
【0034】
【数4】
【0035】となる。ここで、γは磁気回転比、Gpn
はn番目の位相方向傾斜磁場強度(ただし、スピンが感
じる強度とするため、180°パルス1つに付き正負が
一回反転するものとする)、tnはn番目の位相方向傾
斜磁場の印加時間、nは目的のエコー番号のエコー信号
を計測するまでの位相方向に印加する傾斜磁場の個数、
(F.O.V.)は撮像領域の一辺の長さとする。同様
にK空間上でky方向に対し中心からm番目のデータ
は、
【0036】
【数5】
【0037】が成立するように、位相方向の傾斜磁場強
度、印加時間が決定されている。
【0038】そして、この実施例では、特にこのように
撮像した生データ(図1(b))をいわゆるハーフスキ
ャン法で再生する。
【0039】ここで、いわゆるハーフスキャン法の原理
について説明する。ハーフスキャン法とは、フーリエ変
換の計測生データの複素共役対称性と、磁気共鳴イメー
ジング装置の位相誤差は低域となる性質とを用いて、撮
像時間を短縮する方法である。磁気共鳴イメージング装
置の撮像法は、前記したように、計測生データ(以下、
生データと称す)を2次元フーリエ変換を行うことによ
り画像を得ることができる。逆にいうと、画像を逆フー
リエ変換したものが生データとなる。つまり、画像は実
部のみであるとすると虚部が零である。これを逆フーリ
エ変換したものは、実部は偶関数、虚部は奇関数とな
る。
【0040】
【数6】 Real[ρ(kx,ky)]= Real[ρ(−kx,−ky)] (6)
【0041】
【数7】 Imag[ρ(kx,ky)]=−Imag[ρ(−kx,−ky)] (7) ここで、ρ( )は生データ(複素数)、Real[ ]
は実部、Imag[ ]は虚部とする。
【0042】すなわち、理想的な生データは、K空間上
の原点を中心に点対称な位置のデータは、互いに複素共
役の関係が存在することになる。
【0043】しかし、実際の装置においては、様々な位
相誤差が存在することがあり、それが低域に存在するこ
とが分かっている。このため、低域の生データをローパ
スフィルタで切り出し、2次元フーリエ変換し、そのベ
クトル量の位相を求めることにより、画像空間上での位
相誤差が推定できる。この位相と前記の性質を用いるこ
とにより、約半分の計測データで画像が作成できる。
【0044】図3は上述したハーフスキャン法の概略を
示す説明図である。同図(1)は、計測データS(k
x,ky)を示している。(a)では、該計測データを
ky方向にローパスフィルタをかけ、これにより、該計
測データのほぼ中央付近22を切り出す(図3
(2))。
【0045】この際、ky方向の切り出しのデータは、
同一エコーとする。これは、各エコー間でも位相に差が
生じるためである。ここでは、ky方向には、全てのデ
ータを切り出しているが、ky方向もローパスフィルタ
をかけ、これにより中央付近を切り出すようにしてもよ
い。(b)では、前記のように切り出したデータ以外の
高域側データ空間に零を挿入して、2次元フーリエ変換
を行う。
【0046】生データは、複素数(ベクトル)である
が、ここでその位相θ(x,y)を求める(図3
(3))。求められたものは推定位相マップと称されて
いる。
【0047】次に、計測データS(kx,ky)の未計
測部分に零を詰めて、2次元フーリエ変換を行い画像を
作成し、それに(3)の推定位相マップθ(x,y)の
位相分だけ逆方向に回転させる。次にそれを2次元逆フ
ーリエ変換し、未計測部分23のみを切り出す(図3
(4))。
【0048】これが、未計測部分の推定データS(k
x,ky)23となる。(d)では、計測データ(図3
(1))と推定データ(図3(4))とをK空間上で合
わせ込み、その後2次元フーリエ変換を行い画像(図3
(5))を作成する。この場合、通常該合わせ込みの際
に、推定データと計測データとのつなぎが滑らかになる
ようにフィルタをかける。
【0049】なお、以上のハーフスキャン法に関して
は、「位相マップを用いた半数計測データからのMRI
画像再生アルゴリズム」(電子情報通信学会論文紙D
VolJ71−D No.1 pp182−187 1
988年1月)に詳述されている。
【0050】上述の実施例では、4つのエコー信号を使
用して各エコー信号のエコー時間の画質を作成するため
の説明を行ったが、他のエコー数の画質や何エコーまで
使用するかによりデータの並び順が変化することになる
ことはいうまでもない。
【0051】たとえば、図4(a)に2つのエコー信号
を使用した第1エコー信号のエコー時間の画質のデータ
並び順を、図4(b)に2つエコー使用の第2エコーの
エコー時間の画質のデータ並び順を示している。
【0052】また、図5(a)に3つのエコー信号を使
用した第1エコー信号のエコー時間の画質のデータ並び
順を、図5(b)に3つのエコー信号を使用した第2エ
コーのエコー時間の画質のデータ並び順を示している。
【0053】また、図6に3つのエコー信号使用の各エ
コーのエコー時間の画質のデータ並び順を示している。
【0054】また、図7(a)に4つのエコー信号を使
用した第1エコーのエコー時間の画質のデータ並び順
を、図7(b)に4つのエコー信号を使用した第2エコ
ーのエコー時間の画質のデータ並び順を示している。
【0055】さらに、図8(a)に4つのエコー信号を
使用した第3エコーのエコー時間の画質のデータ並び順
を、図8(b)に4つのエコー信号を使用した第4エコ
ーのエコー時間の画質のデータ並び順を示している。
【0056】上述した実施例で示す磁気共鳴イメージン
グ装置は、図2に示す高速スピンエコー法を用いてK空
間上の約半数のデータを取り入れ、その後、いわゆるハ
ーフスキャン再構成法によって該K空間上に全てのデー
タを埋めるようにしたものである。
【0057】この場合、ハーフスキャン再構成法を行な
うにはそれほどの時間がかからないことから、従来高速
スピンエコー法のみで行なう場合と比べて、約半分の時
間で済むようになる。
【0058】したがって、計測時間をさらに短くするこ
とができるようになる。
【0059】なお、上述した実施例においての並び順と
しては、使用するエコー数は2ないし4つのエコーのみ
を示したが、それ以上のエコー数に関しても同様に行う
ことができるのはもちろんである。
【0060】なお、本発明は図2に示した高速スピンエ
コー法、すなわちスピンの励起を90°−180°−1
80°−…系列で行なうもの以外にも適用できる。図1
5はCPMG法と称されるパルスシーケンスの模式図で
ある。このCPMG法はスライス方向傾斜磁場の印加と
ともに90°パルスを印加してスピンを選択励起し、ス
ピンが拡散した所定時間後、たとえばTe/2後にスラ
イス方向傾斜磁場とエコー信号発生のための高周波パル
スを印加して、スピンを90°励起面で180°位相回
転させ、その後位相エンコードをし、周波数方向に信号
を読みだし、その後図15に示すように、高周波パルス
と傾斜磁場の印加によりエコー信号を順次読みだすとい
うものである。このようなCPMG法にも本発明を適用
することができる。
【0061】
【発明の効果】以上説明したことから明らかなように、
本発明による磁気共鳴イメージング装置によれば、計測
時間をさらに短縮することができるようになる。
【図面の簡単な説明】
【図1】 本発明の一実施例を示す説明図で、第4エコ
ー信号までのデータを使用したK空間上のデータの並び
順を示した説明図である。
【図2】 本発明の一実施例を示す説明図で、高速スピ
ンエコー法のパルスシーケンスの模式的説明図である。
【図3】 本発明の一実施例を示す説明図で、ハーフス
キャン法を示した説明図である。
【図4】 本発明の他の実施例を示す説明図で、第2エ
コー信号までのデータを使用したK空間上のデータの並
び順を示した説明図である。
【図5】 本発明の他の実施例を示す説明図で、第3エ
コー信号までのデータを使用したK空間上のデータの並
び順を示した説明図である。
【図6】 本発明の他の実施例を示す説明図で、第3エ
コー信号までのデータを使用したK空間上のデータの並
び順を示した説明図である。
【図7】 本発明の他の実施例を示す説明図で、第4エ
コー信号までのデータを使用したK空間上のデータの並
び順を示した説明図である。
【図8】 本発明の他の実施例を示す説明図で、第4エ
コー信号までのデータを使用したK空間上のデータの並
び順を示した説明図である。
【図9】 本発明の一実施例を示す説明図で、磁気共鳴
イメージング装置の概略を示したブロック図である。
【図10】 巨視的磁化の模式図である。
【図11】 スピンエコー法のパルスシーケンスを示し
た説明図である。
【図12】 マルチエコー法のパルスシーケンスを示し
た説明図である。
【図13】 従来における高速スピンエコー法のパルス
シーケンスを示した説明図である。
【図14】 図13による高速スピンエコー法でK空間
にデータを格納する場合の説明図である。
【図15】 CPMG法のパルスシーケンスの模式的説
明図である。
【符号の説明】
28…90°パルス、29…180°パルス、30…ス
ライス方向傾斜磁場、31,32…位相方向傾斜磁場、
33,34…周波数方向傾斜磁場、24…第1エコー信
号、25…第2エコー信号、26…第3エコー信号、2
7…第4エコー信号。

Claims (3)

    (57)【特許請求の範囲】
  1. 【請求項1】 静磁場中の被検体への高周波パルスと傾
    斜磁場とを逐次所定値で所定順に印加するパルスシーケ
    ンスを複数回実行し、被検体の所定断面から放出される
    エコー信号からフーリエ変換により画像再構成を行う磁
    気共鳴イメージング装置において、 1回のパルスシーケンスで複数個のエコー信号を計測す
    る計測手段と、前記計測手段でのエコー信号計測を、位
    相エンコード方向傾斜磁場及びエコー信号読みだし傾斜
    磁場の直前に印加する位相オフセット傾斜磁場の値の制
    御により、K空間を複数のエコー領域に分割し任意のエ
    コー番号を前記K空間の中央部を含むエコー領域に配置
    すると共に、前記K空間の一方の端部から前記中央部の
    エコー領域を越える少なくとも次のエコー領域までのエ
    コー番号が前記任意のエコー番号から連続し、前記K空
    間の一方の端部のエコー番号と前記次のエコー領域のエ
    コー番号とが連続するように配置する制御手段と、前記
    K空間の他方の端部のエコー番号が前記次のエコー領域
    エコー番号から連続するように推定で埋める推定手段
    とを備えたことを特徴とする磁気共鳴イメージング装
    置。
  2. 【請求項2】 請求項1に記載の磁気共鳴イメージング
    装置において、 前記パルスシーケンスは、高周波パルスとして、1回の
    90゜パルスと該90゜パルスに続く複数個の180゜
    パルスとからなることを特徴とする磁気共鳴イメージン
    グ装置。
  3. 【請求項3】 静磁場中の被検体への高周波パルスと傾
    斜磁場とを逐次所定値で所定順に印加するパルスシーケ
    ンスを複数回実行し、被検体の所定断面から放出される
    エコー信号からフーリエ変換により画像再構成を行う磁
    気共鳴イメージング装置において、 単一のパルスシーケンスで複数個のエコー信号を計測す
    る計測手段と、位相エンコード方向傾斜磁場またはエコ
    ー信号読みだし傾斜磁場の直前に印加する位相オフセッ
    ト傾斜磁場の値を制御することにより、前記計測手段で
    のエコー信号計測を行わせる制御手段と、前記計測デー
    タからK空間上における未計測データを推定により埋め
    る推定手段とを備え、前記制御手段は、前記K空間を位
    相エンコード方向に不均等に2分割して大となる領域を
    計測範囲とし、該計測領域に任意のエコー番号を配置す
    ると共に、前記計測範囲のエコー番号が前記任意のエコ
    ー番号から連続し、かつ前記計測範囲の一方の端部のエ
    コー番号が他方の端部のエコー番号に連続するように配
    置し、前記推定手段は前記K空間の他方の端部のエコー
    番号が前記計測範囲の内で前記K空間の端部となる側の
    エコー番号から連続するように未計測データを推定する
    ことを特徴とする磁気共鳴イメージング装置。
JP13854092A 1992-05-29 1992-05-29 磁気共鳴イメージング装置 Expired - Fee Related JP3276669B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP13854092A JP3276669B2 (ja) 1992-05-29 1992-05-29 磁気共鳴イメージング装置
US08/067,624 US5357200A (en) 1992-05-29 1993-05-28 Magnetic resonance imaging method and apparatus
US08/274,281 US5543709A (en) 1992-05-29 1994-07-13 Magnetic resonance imaging method and apparatus for generating a plurality of single spin echoes

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP13854092A JP3276669B2 (ja) 1992-05-29 1992-05-29 磁気共鳴イメージング装置

Publications (2)

Publication Number Publication Date
JPH05329127A JPH05329127A (ja) 1993-12-14
JP3276669B2 true JP3276669B2 (ja) 2002-04-22

Family

ID=15224547

Family Applications (1)

Application Number Title Priority Date Filing Date
JP13854092A Expired - Fee Related JP3276669B2 (ja) 1992-05-29 1992-05-29 磁気共鳴イメージング装置

Country Status (2)

Country Link
US (2) US5357200A (ja)
JP (1) JP3276669B2 (ja)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3276669B2 (ja) * 1992-05-29 2002-04-22 株式会社日立メディコ 磁気共鳴イメージング装置
DE4415393B4 (de) * 1993-06-01 2010-10-14 Siemens Ag Verfahren zur Erzeugung von Bildern in einem Kernspintomographiegerät mit einer Spinecho-Pulssequenz
AUPP340698A0 (en) * 1998-05-07 1998-05-28 Canon Kabushiki Kaisha Method and apparatus for highly efficient representation and compression of images
WO2000037956A1 (en) * 1998-12-22 2000-06-29 Koninklijke Philips Electronics N.V. Enhanced t2*-contrast in half-fourier magnetic resonance imaging
JP4225620B2 (ja) * 1999-01-11 2009-02-18 株式会社東芝 Mri装置
US6486667B1 (en) * 2000-03-31 2002-11-26 Koninklijke Philips Electronics N.V. Combination of fluid-attenuated inversion-recovery complex images acquired using magnetic resonance imaging
WO2001076482A2 (en) * 2000-04-07 2001-10-18 Washington University Method and apparatus for compressing computed tomography raw projection data
JP3853585B2 (ja) * 2000-10-11 2006-12-06 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Mri装置
US7327866B2 (en) * 2001-04-09 2008-02-05 Bae Kyongtae T Method and apparatus for compressing computed tomography raw projection data
JP2003038456A (ja) * 2001-07-10 2003-02-12 Ge Medical Systems Global Technology Co Llc スピン励起方法、磁気共鳴撮影方法および磁気共鳴撮影装置
US7949169B2 (en) * 2002-12-04 2011-05-24 Bae Kyongtae T Method and apparatus for automated detection of target structures from medical images using a 3D morphological matching algorithm
JP5100374B2 (ja) * 2005-04-26 2012-12-19 株式会社日立メディコ 磁気共鳴イメージング装置および方法
JP5591500B2 (ja) * 2009-07-30 2014-09-17 株式会社日立メディコ 磁気共鳴イメージング装置、および、画像再構成プログラム
WO2014042026A1 (ja) 2012-09-11 2014-03-20 株式会社 東芝 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
JP6084430B2 (ja) * 2012-10-30 2017-02-22 株式会社日立製作所 再収束rfパルスのフリップ角制御法及び磁気共鳴イメージング装置

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3434161A1 (de) * 1984-09-18 1986-03-27 Bruker Medizintechnik Gmbh, 7512 Rheinstetten Verfahren zum messen der magnetischen kernresonanz
US4684891A (en) * 1985-07-31 1987-08-04 The Regents Of The University Of California Rapid magnetic resonance imaging using multiple phase encoded spin echoes in each of plural measurement cycles
US4833407A (en) * 1987-06-24 1989-05-23 Picker International, Inc. Scan time reduction using conjugate symmetry and recalled echo
US4952876A (en) * 1988-11-23 1990-08-28 General Electric Company Variable bandwidth multiecho NMR imaging
US4937526A (en) * 1988-11-23 1990-06-26 Mayo Foundation For Medical Education And Research Adaptive method for reducing motion and flow artifacts in NMR images
US5239266A (en) * 1990-08-03 1993-08-24 The Regents Of The University Of California MRI using variable imaging parameter(s) within a single image sequence
US5202632A (en) * 1990-08-03 1993-04-13 The Regents Of The University Of California Mri using variable imaging parameter(s) within a single image sequence
JP3153574B2 (ja) * 1991-08-23 2001-04-09 株式会社東芝 磁気共鳴映像装置
JPH05237067A (ja) * 1992-02-27 1993-09-17 Toshiba Corp 磁気共鳴イメージング装置
JP3276669B2 (ja) * 1992-05-29 2002-04-22 株式会社日立メディコ 磁気共鳴イメージング装置
US5311133A (en) * 1992-10-05 1994-05-10 Picker International, Inc. Concurrent generation of multiple images in fast (TR<T2) scanning sequences

Also Published As

Publication number Publication date
US5357200A (en) 1994-10-18
US5543709A (en) 1996-08-06
JPH05329127A (ja) 1993-12-14

Similar Documents

Publication Publication Date Title
US9983284B2 (en) MRI with dixon-type water/fat separation and prior knowledge about inhomogeneity of the main magnetic field
CN104685368B (zh) 用于耐金属mr成像的方法和设备
EP0515197A1 (en) Acquisition of multiple images in fast spin echo nmr scans
EP3295194B1 (en) Parallel mr imaging with rf coil sensitivity mapping
JPH0654828A (ja) Nmrシステム
JP3276669B2 (ja) 磁気共鳴イメージング装置
JPH0420618B2 (ja)
CN107076818B (zh) 零回波时间mr成像
JPH06500947A (ja) Nmrシステム及び別々に収集された1組のnmr信号から像を発生する方法
US4472683A (en) Imaging apparatus using nuclear magnetic resonance
US4684892A (en) Nuclear magnetic resonance apparatus
JP3516421B2 (ja) Mri装置
JPH11225995A (ja) 磁気共鳴イメージング装置
JP3526347B2 (ja) 磁気共鳴イメージング装置
JP3332951B2 (ja) 磁気共鳴イメージング装置
JP3167038B2 (ja) 磁気共鳴イメージング装置
EP4139699B1 (en) Mr imaging with t1 compensated b1 mapping
JP7645739B2 (ja) 磁気共鳴イメージング装置およびその制御方法
JP3163125B2 (ja) Mri装置
JP2000316830A (ja) 磁気共鳴イメージング方法及びそれを用いた磁気共鳴イメージング装置
JP3201649B2 (ja) 磁気共鳴イメージング装置
JP3322695B2 (ja) 磁気共鳴イメージング装置
JP2002143121A (ja) 磁気共鳴イメージング装置
JP3369586B2 (ja) 磁気共鳴イメージング装置
JP3317552B2 (ja) Mri装置

Legal Events

Date Code Title Description
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080208

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090208

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090208

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100208

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100208

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110208

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110208

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120208

Year of fee payment: 10

LAPS Cancellation because of no payment of annual fees