JP3929047B2 - 磁気共鳴イメージング装置 - Google Patents
磁気共鳴イメージング装置 Download PDFInfo
- Publication number
- JP3929047B2 JP3929047B2 JP2003119403A JP2003119403A JP3929047B2 JP 3929047 B2 JP3929047 B2 JP 3929047B2 JP 2003119403 A JP2003119403 A JP 2003119403A JP 2003119403 A JP2003119403 A JP 2003119403A JP 3929047 B2 JP3929047 B2 JP 3929047B2
- Authority
- JP
- Japan
- Prior art keywords
- magnetic resonance
- resonance imaging
- imaging apparatus
- echo signal
- measurement
- 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
- 238000002595 magnetic resonance imaging Methods 0.000 title claims description 90
- 238000005259 measurement Methods 0.000 claims description 91
- 238000012545 processing Methods 0.000 claims description 58
- 230000035945 sensitivity Effects 0.000 claims description 49
- 238000005481 NMR spectroscopy Methods 0.000 claims description 10
- 230000002194 synthesizing effect Effects 0.000 claims 2
- 238000000034 method Methods 0.000 description 107
- 238000005070 sampling Methods 0.000 description 42
- 238000003384 imaging method Methods 0.000 description 36
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 28
- 230000008569 process Effects 0.000 description 15
- 238000010586 diagram Methods 0.000 description 13
- 230000033001 locomotion Effects 0.000 description 13
- 238000004364 calculation method Methods 0.000 description 9
- 238000002592 echocardiography Methods 0.000 description 9
- 230000005540 biological transmission Effects 0.000 description 7
- 230000003068 static effect Effects 0.000 description 6
- 238000007796 conventional method Methods 0.000 description 4
- 229940079593 drug Drugs 0.000 description 4
- 239000003814 drug Substances 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 230000000747 cardiac effect Effects 0.000 description 2
- 230000003111 delayed effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 238000004904 shortening Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005281 excited state Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
- G01R33/4824—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image 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
- G01R33/5611—Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General 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)
- Magnetic Resonance Imaging Apparatus (AREA)
Description
【発明の属する技術分野】
本発明は、核磁気共鳴(NMR)現象を利用して被検体の検査部位の断層画像を得る磁気共鳴イメージング装置(以下「MRI装置」という)に関し、特にプロペラサンプリング法と複数の高周波受信コイルを用いたパラレルイメージング法とを組み合わせて、撮影時間と画像再構成時間を短縮する技術に関する。
【0002】
【従来の技術】
MRI装置では、撮影中に被検体が動いた場合、その影響は画像全体に及び、位相エンコード方向に画像が流れた様なアーチファクト (以下、「体動アーチファクト」と呼ぶ) が生じる。これは、一般的に計測空間(いわゆる「k−空間」といわれる空間)上の各格子点のエコー信号をサンプリングするに際し、周波数エンコード方向に平行なサンプリングを位相エンコード方向に繰り返す(以下、「直交系(Cartesian)サンプリング法」という)為である。
【0003】
この直交系(Cartesian)サンプリング法に対して、非直交系(Non-cartesian)サンプリングと呼ばれる手法があり、ラディアルサンプリング法(非特許文献1)、プロペラサンプリング法(非特許文献2)がその一例である。これらは高速に撮影を行うことができるため、体動アーチファクトを低減することができ、リアルタイムイメージングや高速に拍動する心臓撮影に適用されつつある。
【0004】
ラディアルサンプリング法は、計測空間の略一点を回転中心として回転角を変えながら放射状にサンプリングを行い、一枚の画像再構成に必要なエコー信号を得る技術である。ラディアルサンプリング法を用いて撮影を行った場合、放射状にサンプリングを行うことから、体動アーチファクトが画像の周辺に散らばる(つまり、注目すべき視野の外側へ出る)ため、他の一般的サンプリング法の撮影と比較して体動アーチファクトが目立たなくなり、体動に対してロバストといわれている。
【0005】
プロペラサンプリング法は、ラディアルサンプリング法に位相エンコードの概念を適用した手法である。このプロペラサンプリング法では、全計測を複数のブロックに分割し、各ブロックを異なる計測空間の回転角で計測する。各ブロック内では、1つの回転角に対して、通常計測の場合と同様に位相エンコードパルスを付加して複数のエコー信号を取得し、これを回転角毎に繰り返して一枚の画像再構成に必要な全てのエコー信号を取得する。このようにプロペラサンプリング法では、計測空間の低空間周波数領域を重複して計測するため、体動による信号の変動を抽出でき体動補正が可能となるとの報告(非特許文献2)もある。
【0006】
また、MRI装置における撮影時間短縮の手法として、パラレルイメージング法が公知(非特許文献3)である。パラレルイメージング法では、局所的で不均一な感度分布を有する複数の高周波受信コイルを組み合わせて成るマルチプル高周波受信コイルを用いて、再構成画像上に生じた折り返しアーチファクトを除去する手法であり、一枚の画像再構成に必要なエコー数を低減できることから、パルスシーケンスの繰り返し回数を減らせ、撮影時間を短縮することが可能である。
【0007】
さらなる高速化を実現するために、非直交系サンプリング法の一種であるスパイラルサンプリング法とパラレルイメージング法を融合する例も公知(非特許文献4)である。スパイラルサンプリング法では、計測空間の略一点を中心として螺旋状にサンプリングを行い、一枚の画像再構成に必要なエコー信号を得る。
【0008】
【非特許文献1】
G.H.Glover et.al.,Projection Reconstruction Techniques for Redution of Motion Effects in MRI,Magnetic Resonance in Medicine 28:275-289 (1992)
【0009】
【非特許文献2】
Magnetic Resonance in Medicine 42:963-969,1999, James G.Pipe,Motion Correction With PROPELLER MRI:Application to Head Motion and Free-Breathing Cardiac Imaging
【0010】
【非特許文献3】
Magnetic Resonance in Medicine 30:142-145,1993,J.B.Ra,C.Y.Rim,Fast Imaging Using Subencoding Data Sets from Multiple Detectors
【0011】
【非特許文献4】
Magnetic Resonance in Medicine 46:638-651,2001,Klaas P.Pruessmann, Markus Weiger,Peter Bo¨rnert,and Peter Boesiger,Advances in Sensitivity Encoding With Arbitrary k-Space Trajectories
【0012】
【発明が解決しようとする課題】
以上のような非直交サンプリング法にパラレルイメージング法を適用する場合、一般化した折り返しの展開アルゴリズムによる折り返しアーチファクトの除去が必要であり、そのため演算時間が非常に長くなる。
【0013】
さらに、非直交系サンプリング法で取得されたエコー信号データには、グリッディングと呼ばれる非直交系サンプリングデータを直交系サンプリングデータに変換する補間処理を実施した後、2次元/3次元高速フーリエ変換して画像を再構成する必要がある。このため、グリッディング処理にも多くの時間が必要である。
【0014】
以上のことから、非直交サンプリング法に一般化したパラレルイメージングを適用すると、計測は早くなるが演算が遅くなり、結果としてリアルタイム撮影が実現できない。あるいは、演算速度を上げようとすると、非常に高性能なコンピュータを導入する必要があるため、コストパフォーマンスが低下する。
【0015】
また、プロペラサンプリング法では、ラディアルサンプリング法と比較して、1つの回転角すなわち各ブロックで異なる位相エンコード量を印加して複数のエコー信号を取得するため、1枚の画像再構成に必要な全てのエコー信号を取得するための時間が延長する。この時間延長により、画像に被検体の動きによる体動アーチファクトが混入したり、リアルタイム性が低下する。
【0016】
そこで、本発明は、プロペラサンプリング法に関して、1つのブロックで取得するエコー数を低減しても折り返しアーチファクトを生ずることなく撮影時間の短縮を図り、かつ、演算量を低減して高速演算を可能とすることにより、全体として高速イメージングを可能とすることを目的とする。
【0017】
【課題を解決するための手段】
前記課題を解決するために、本発明は次の様に構成される。即ち、被検体から発生する核磁気共鳴信号を受信する受信手段と、k空間の原点を含み低空間周波数領域から高空間周波数領域に亘って所定の幅を有して成る単位領域に対応するエコー信号を計測する単位計測を、前記k空間の原点を中心にして前記単位領域の回転角度を変えて繰り返す計測制御手段と、前記単位領域の各々に対応するエコー信号を用いて前記被検体の画像を再構成する信号処理手段と、を有する磁気共鳴イメージング装置において、
前記受信手段は、複数の受信コイルを組み合わせて成るマルチプル受信コイルを有し、
前記計測制御手段は、前記受信コイル毎にエコー信号を計測すると共に、少なくとも一つの前記単位計測においてエコー信号の計測を間引き、
前記信号処理手段は、前記エコー信号の計測を間引かれた単位領域を含む複数の単位領域の各々に対応するエコー信号と前記複数の受信コイルの感度分布データとを用いて、折り返しアーチフアクトの無い被検体画像を再構成することを特徴とする。
【0027】
以上の手段を備えることにより、プロペラサンプリング法にパラレルイメージング法を適用して、1つのブロックで取得するエコー数を低減し、撮影時間の短縮を図り、かつ、高速演算が可能になる。
【0028】
【発明の実施の形態】
図4は本発明による磁気共鳴イメージング装置の全体構成を示すブロック図である。この磁気共鳴イメージング装置は、核磁気共鳴(NMR)現象を利用して被検体の断層画像を得るもので、図4に示すように、静磁場発生系2と、傾斜磁場発生系3と、送信系5と、受信系6と、信号処理系7と、シーケンサ4と、中央処理装置(CPU)8とを備えて構成される。
【0029】
静磁場発生系2は、被検体1の周りの空間にその体軸方向または体軸と直交する方向に均一な静磁場を発生させるもので、被検体1の周りに永久磁石方式または常電導方式あるいは超電導方式の磁場発生手段が配置されている。
【0030】
傾斜磁場発生系3は、X,Y,Zの3軸方向に巻かれた傾斜磁場コイル9と、それぞれの傾斜磁場コイルを駆動する傾斜磁場電源10とから成り、後述のシ−ケンサ4からの命令に従ってそれぞれのコイルの傾斜磁場電源10を駆動することにより、X,Y,Zの3軸方向の傾斜磁場Gs,Gp,Grを被検体1に印加する。より具体的には、X,Y,Zのいずれかの1方向にスライス傾斜磁場パルス(Gs)を印加して被検体9に対するスライス面を設定し、残り2つの方向に位相エンコード傾斜磁場パルス(Gp)と周波数エンコード傾斜磁場パルス(Gr)を印加して、エコー信号にそれぞれの方向の位置情報をエンコードする。
【0031】
シーケンサ4は、高周波磁場パルス(以下、「RFパルス」という)と傾斜磁場パルスをある所定のパルスシーケンスで繰り返し印加する制御手段で、CPU8の制御で動作し、被検体1の断層画像のデータ収集に必要な種々の命令を送信系5、傾斜磁場発生系3、および受信系6に送る。
【0032】
送信系5は、被検体1の生体組織を構成する原子の原子核スピンに核磁気共鳴を起こさせるためにRFパルスを照射するもので、高周波発振器11と変調器12と高周波増幅器13と送信側の高周波コイル14aとから成る。高周波発振器11から出力された高周波パルスをシーケンサ4からの指令によるタイミングで変調器12により振幅変調し、この振幅変調された高周波パルスを高周波増幅器13で増幅した後に被検体1に近接して配置された高周波コイル14aに供給することにより、電磁波(RFパルス)が被検体1に照射される。
【0033】
受信系6は、被検体1の生体組織を構成する原子核スピンの核磁気共鳴により放出されるエコー信号(NMR信号)を検出するもので、受信側の高周波コイル14bと増幅器15と直交位相検波器16と、A/D変換器17とから成る。送信側の高周波コイル14aから照射された電磁波によって誘起される被検体1の応答の電磁波(NMR信号)が被検体1に近接して配置された高周波コイル14bで検出され、増幅器15で増幅された後、シーケンサ4からの指令によるタイミングで直交位相検波器16により直交する二系統の信号に分割され、それぞれがA/D変換器17でディジタル量に変換されて、信号処理系7に送られる。
【0034】
信号処理系7は、光ディスク19、磁気ディスク18等の外部記憶装置と、CRT等からなるディスプレイ20とを有し、受信系6からのデータがCPU8に入力されると、CPU8が信号処理、画像再構成等の処理を実行し、その結果である被検体1の断層画像をディスプレイ20に表示すると共に、外部記憶装置の磁気ディスク18等に記録する。
なお、図4において、送信側及び受信側の高周波コイル14a,14bと傾斜磁場コイル9は、被検体1の周りの空間に配置された静磁場発生系2の静磁場空間内に設置されている。
【0035】
上記の様なMRI装置の撮影対象スピン種は、臨床で普及しているものとしては、被検体の主たる構成物質であるプロトンである。プロトン密度の空間分布や、励起状態の緩和現象の空間分布を画像化することで、人体頭部、腹部、四肢等の形態または、機能を2次元もしくは3次元的に撮影する。
【0036】
次に、撮影方法を説明する。図5は直交系サンプリング法のグラディエントエコーパルスシーケンスである。図5のRF,Gs,Gp,Gr,A/D,echoはそれぞれ、RFパルス,スライス傾斜磁場,位相エンコード傾斜磁場,周波数エンコード傾斜磁場,AD変換,エコー信号を表し、501は高周波パルス,502はスライス傾斜磁場パルス、503は位相エンコード傾斜磁場パルス、504は周波数エンコード傾斜磁場パルス、505はサンプリングウインド、506はエコー信号、507は繰り返し時間(501の間隔)である。
【0037】
MRIでは、繰り返し507毎に位相エンコード傾斜磁場パルス503の量(=傾斜磁場パルス波形と時間軸との囲む面積)を変えて異なる位相エンコード量を与え、それぞれの位相エンコードで得られるエコー信号506を検出する。この操作を位相エンコードの数だけ繰り返し、画像取得時間508で1枚の画像再構成に必要なエコー信号を取得する。位相エンコードの数は通常1枚の画像あたり64,128,256,512等の値が選ばれる。各エコー信号は通常128,256,512,1024個のサンプリングデータからなる時系列信号として得られる。これらのデータを2次元フーリエ変換して1枚のMR画像を作成する。
【0038】
図5のパルスシーケンスを用いてサンプリングしたエコーを計測空間に配置した結果を図8(a)に示す。図8(a)では、507を12回繰り返し、それぞれ506-1〜506-cのエコー信号を取得した場合であり、エコー信号の沿え字の順序は、エコー信号を取得した時間(繰り返し507の回数)の順序と対応する。すなわち、沿え字の番号の小さいものは、エコー信号を取得した時間が早く、番号の大きなものは、エコー信号を取得した時間が遅い(以下、本文中の記述は全て同様である)。図8(a)では、エコー信号群506-1〜506-cは、計測空間801のKy方向に上から順番に配置されるようにパルスシーケンスを実行して取得した場合である。
【0039】
次に、ラディアルサンプリング法(以下、「ラディアルMRI法」という)を用いたグラディエントパルスシーケンスを図6に示す。図5で示した直交系サンプリング法との違いは、位相エンコード傾斜磁場の軸Gpにも周波数エンコード傾斜磁場と同様の傾斜磁場601を印加することと、パルスシーケンスの各繰り返し507でGp,Gr軸に印加する傾斜磁場601,602の振幅を変えることである。パルスシーケンスの各繰り返し507で異なる量の601,602の振幅を用いることで、計測空間の略一点を中心として放射状に回転したデータを取得できる。
【0040】
通常、ラディアルMRI法で取得したデータは、計測空間の規則正しい格子座標にはのらない。そこで、サンプリングしたデータから補間処理を用いて、規則正しい格子座標上のデータを作成するグリッディング処理が必要である。図6のパルスシーケンスを用いてサンプリングしたエコー信号を計測空間に配置した結果を図8(b)に示す。図8(b)は、507を12回繰り返し、それぞれ603-1〜603-cのエコー信号を取得した場合である。この場合、計測空間を均等に埋めるために必要な回転Δθ(802)は、
Δθ=π/12 (1)
であり、このようなサンプリングを行うため、Gp,Gr軸の傾斜磁場出力GP,GRは、直交系サンプリング法で用いる周波数エンコード傾斜磁場の出力をG、繰り返し番号をn(1≦n≦12)とした時、
GP=G×sin(Δθ×(n-1)) (2)
GR=G×cos(Δθ×(n-1)) (3)
である。ラディアルMRI法のメリットは、直交系サンプリング法の場合と比較して、少ないエコー数でも分解能が良いことと、ゴーストアーチファクトが発生しにくいことである。また、撮影時の回転角を調整することで、被検体の体動による画質への影響を少なくできる。
【0041】
ここで、計測空間801へのデータの配置について述べる。図8(c)は、非直交系サンプリング法で取得したデータを、グリッディング処理して計測空間に配置した場合の模式図である。計測空間801は、図8(c)の黒丸で示すように規則正しい格子点の座標を有する。しかし、非直交系サンプリングで取得したデータは803-1〜803-3の様に計測空間に対して異なる軌跡(座標)を通るので、白丸で示すサンプリングされたデータは計測空間の格子点座標(図8(c)中の黒丸)と一致しない。グリッディング処理では、これらサンプリングされたデータ(図8(c)中の白丸)を用いて、補間処理により規則正しい格子点座標(図8(c)中の黒丸)にデータを再配置することをいう。なお、グリッディング処理は、例えばSinc関数やKaiser-Bessel関数の補間用関数を用いて行う(非特許文献5)。
【非特許文献5】
J.I Jackson et.al.,Selection of a Convolution Function for Fourier Inversion Using Gridding,IEEE Trans.Med.Imaging,vol.10,pp.473-478,1991
【0042】
さらに、ラディアルMRI法に位相エンコードを付加した、プロペラサンプリング法(以下、「プロペラMRI法」という)を用いたグラディエントエコーパルスシーケンスを図7に示す。図7では、計測を4つのブロック701〜704に分割し、それぞれのブロック内で5個のエコー信号を取得した場合である。図6の場合との違いは、ブロック701〜704毎に異なる計測空間上の回転角となるように、傾斜磁場出力を設定し、それぞれのブロック701〜704内では、位相エンコードパルスを付加していることである。図7では、回転角度に応じて出力が変化する傾斜磁場ブロック707と、707の出力前に更に位相エンコード量を付加して出力するブロック706とで構成される。
【0043】
このようにして取得したエコー信号を、計測空間に配置した結果を、図8(d)に示す。この場合では、各ブロック701〜704毎に、角度804で回転し、それぞれのブロック内では位相エンコード量が異なる5個のエコー信号を取得している。
【0044】
次に、取得するエコー信号と、画像の折り返しを説明する。図2(a)は、図5のパルスシーケンスを実行して取得したエコー信号を計測空間201に配置したものであり、201をフーリエ変換した画像が202である。これに対し、付加する位相エンコードパルス503を2倍のステップで実行し、エコー信号を半分だけ取得したのが図2(b)である。この場合、計測空間204は1ラインおきに配置される(204内の実線が計測したエコー信号、点線が計測しなかったエコー信号をそれぞれ表す)。この場合、取得したエコー信号が1ラインおきであることから、計測したエコー信号のラインのみを有効とし、計測空間を半分の大きさに減少したのが205であり、この減少した計測空間205をフーリエ変換した画像が206となる。
【0045】
この場合、撮影時の空間分解能が同じまま、再構成マトリクスが半分になったことから、画像206には折り返しのアーチファクトが発生する。この折り返しアーチファクトは、図2(a)の全てのエコー信号を取得した場合の画像202のうち、点線で区切られた領域の上側部分203-1と、下側部分203-2を重ね合わせたものと同様であると考えることができる。このように、MRI装置では、取得するエコー数を間引いて計測を行った場合、画像にはその間引きステップに応じた折り返しアーチファクトが発生することが知られている。以下、折り返しアーチファクトの発生した画像を折り返し画像と呼ぶ。
【0046】
この折り返し画像を、高周波受信コイルユニット(以下、「RF受信コイル」という)の感度分布データを用いて展開する手法がパラレルイメージング法(以下、「パラレルMRI法」という)である。図2(c)は2個のRF受信コイルを用いた場合のパラレルMRI法の信号処理である。パラレルMRI法では、各RF受信コイルで折り返し画像206-1,206-2を取得し、各RF受信コイルの感度分布データ207-1,207-2を用いて、行列演算208を行い、折り返しアーチファクトの無い画像209を得る。ここで、感度分布データ207と折り返し画像206の沿え字はRF受信コイルの番号を表す。パラレルMRI法を安定に行うためには、使用する各RF受信コイルの感度分布は、折り返しの発生した方向(位相エンコード方向)に異なっている必要がある(図2ではy方向)。
【0047】
一般的に、各RF受信コイルの感度分布データ207-1,207-2は、画像取得のためのパルスシーケンスに先立ち、専用のパルスシーケンスを実行して作成する。或いは、画像取得のためのパルスシーケンスを利用し、計測空間の低空間周波数領域を計測する時点の位相エンコードステップを密にし、この低空間周波数領域のエコー信号を用いて作成しても良い(プロペラMRI法に適用した場合は後述する)。一般的にコイルの感度分布は空間的に穏やかに変動するので、低空間周波数領域のエコー信号データのみで十分であり、逆にその方が高空間周波数領域のデータまで考慮することによるノイズの混入を防止でき、感度分布データのS/N比を向上できるので好都合である。ただし、この場合は、専用パルスシーケンスを用いた場合と比較して、画像取得のためのパルスシーケンスで取得するエコー数が増加するため、撮影時間は延長する。
【0048】
以上の、プロペラMRI法とパラレルMRI法の原理を踏まえて、本発明の第1の実施例を図1,図3を用いて説明する。
本実施例では、4つのブロックで計測を行うプロペラMRI法の場合を説明する。図3(a)は、4つのブロック701〜704内で、それぞれ10個のエコー信号を取得する、通常のプロペラMRI法の場合である。これに対し図3(b)は、それぞれのブロック内で付加する位相エンコード傾斜磁場強度を図3(a)の2倍に設定し、エコー信号数を半分(5個)だけ取得した場合である。図の実線は計測したエコー信号、点線は計測しないエコー信号の位置に対応する。
【0049】
このように計測を行った場合、画像には折り返しアーチファクトが発生するが、プロペラMRI法ではブロック毎に計測空間の角度が異なるため、折り返しアーチファクトは図2(b)のように単純な形状とはならない。そこで、本発明ではそれぞれのブロック701〜704を別々に取り出し計測空間を図3(c)〜(f)の様に分割して考える。通常、プロペラMRI法で取得したデータは、画像再構成のため最終的な座標系である任意の計測空間(図3の場合では301である。以下、この様な座標系を「画像用の座標系」という)にグリッディング処理される。
【0050】
本発明では、グリッディング処理を行わず、それぞれ異なる座標系301〜304のままフーリエ変換を行う(図3の301と303,302と304は90°ずれている)。これにより、それぞれの座標系で作成された折り返し画像305〜308ができる。これら画像は、最終的な画像の座標系(図3の場合は305)に対して計測空間301〜304の関係と同様に回転した形状となる。このように、ブロック毎のデータを別々の座標系を用いて再構成することで、グリッディング処理を行うことなく各ブロックの折り返し画像305〜308を得る。
【0051】
次に、分割した折り返し画像305〜308の再構成処理を図1に示す。図1では、図2(c)と同様に、2個のRF受信コイルを用いて撮影を行った場合の処理であり、RF受信コイルの感度分布データ101〜104及び折り返し画像305〜308の沿え字は、RF受信コイルの番号を表す。このとき、1番目のブロック701の回転角に応じて取得した第1の座標系での各RF受信コイルの感度分布データ101-1,101-2を用いて、701で取得した折り返し画像305-1,305-2を行列演算処理208-1し、折り返しの無い画像105を得る。次いで、逆フーリエ変換109-1によってパラレルMRI法の再構成を行った、第1の座標系での計測空間データ110を得る。
【0052】
同様に、他のブロック702〜704で取得した折り返し画像306〜308についても、それぞれのブロックの座標系に対応したRF受信コイルの感度分布データ102〜104を用いて行列演算処理208を行って折り返しの無い画像106〜108を得た後、逆フーリエ変換109-2〜109-4を行い、それぞれのブロックの座標系での計測空間データ111〜113を作成する。
【0053】
これらの計測空間データ110〜113は、当然であるが座標系301〜304にそれぞれ対応しており、110〜113は互いに異なる座標系である。
【0054】
このようにして作成したそれぞれの座標系での計測空間データ110〜113は、グリッディング処理114-1〜114-4を行い、最終的な画像用の座標系(図1では、第1の座標系301と同じであるが、必ずしも同じである必要は無い)に配置され、信号結合処理115によって画像用の座標系116上で合成(加算または加算平均等の処理)される。最終画像118は、この画像用の座標系116をフーリエ変換117して得る。
【0055】
本実施例の説明では、ブロック701〜704の回転角に応じた各計測空間上で、グリッディング処理せずに、パラレルMRI法による画像再構成を行っている。これは、グリッディング処理後の計測空間に対して一般的化アルゴリズムを使って折り返しの除去を行う従来の方法と比較して、グリッディング処理前の計測空間でそれぞれ折り返しの除去を行うので一般化アルゴリズムが不要となり、かつ、パラレルMRI法の再構成に補間が必要無いこととあわせて、計算量が大幅に低減される。
【0056】
本発明を処理の観点から説明すると、図10の様になる。図10(a)は従来技術のプロペラMRI法とパラレルMRI法の融合の例である。この場合、エコー信号計測121〜124を行って各ブロック701〜704で取得したそれぞれの座標系のエコー信号データを、画像用の座標系にグリッディング処理125した後にフーリエ変換126して折り返し画像を得て、一般化パラレルMRI法の折り返し展開処理127を行う。
【0057】
これに対し、図10(b)の本発明では、それぞれのブロック701〜704で取得したエコー信号データを、それぞれの座標系でフーリエ変換128-1〜128-4して折り返し画像を作成した後、それぞれの座標系でパラレルMRI法の折り返し除去208-1〜208-4を行って折り返しの無い画像を得て、それらを逆フーリエ変換109-1〜109-4してそれぞれの座標系での計測空間データとする。最終画像は、それぞれの座標系で作成した計測空間データを、画像用の座標系にグリッディング処理114-1〜114-4した後、信号結合処理115によって画像用の座標系での最終的な計測空間データを作成して、それをフーリエ変換126する。
【0058】
さらに、本発明の方法では、ブロック701〜704毎にパラレルMRI法の再構成処理が独立であるため、ブロック内のエコー信号データを取得した時点で、そのブロックの画像再構成処理を開始できるので、リアルタイム再構成が可能である。図11は、プロペラMRI法を用いた連続撮影に本発明を適用した場合である。図11の例では、4つのブロック701〜704から成るプロペラMRI法を、連続的に実行する場合であり、横方向が時間経過を表す。図11(a)は従来法に基づいた手順であり、図11(b)が本発明に基づいた手順である。
【0059】
図11(a)の従来法では、パルスシーケンスは時間間隔130-1にて、1枚目の画像用のデータを取得し、次いで時間間隔130-2で2枚目の画像用のデータを取得、以降、連続的に複数枚の画像を取得する。このとき、計測は第1のブロックから始まり、第4のブロックまで取得した後に、第1のブロック(但し、2枚目の画像用)に戻りこれを繰り返す。この従来法では、4つのブロック全てのデータ701〜704を取得した後にグリッディング処理125を行い、一般化パラレル展開131(内部の処理は、126と127を行う)しているため、時間間隔130毎に画像が作成される。
【0060】
これに対し、本発明では、図11(b)のように各ブロック701〜704毎に処理を行うため、パラレル展開処理135(内部の処理は128,208,109を行う)及びグリッディング処理114は、各ブロック701〜704の計測終了と同時にそれぞれ開始することができる。プロペラMRI法では、ブロック番号が同じデータは、計測空間の回転角度が同じ、もしくはほぼ同じとなるように計測するので、グリッディング処理114で作成した計測空間を取捨選択して画像作成137(内部の処理は、115と126)することで、画像更新間隔130の中間の時相の画像を作成可能である。
【0061】
例えば、ブロック701-1,702-1,703-1,704-1のデータを選択して画像作成137-1を行い画像140-1を得、次いで第1ブロックのデータのみ更新した、702-1,703-1,704-1,701-2のデータを用いて画像作成138-1を行い、画像140-1に対してブロック701-2の計測の分だけ時相が遅れた画像141-1を得ることができる。このようにして、画像の更新時間間隔を短くできる。
【0062】
次に、本発明の第2の実施例をブロック数が4のプロペラMRI法を例に説明する。本実施例では、エコー信号は図9(a)のように取得する(実線が計測するエコー信号、点線が計測しないエコー信号の位置に対応する)。図3(b)との違いは、各ブロック内で取得するエコー信号の数及び計測空間上の間隔である。すなわち、各ブロック701〜704内の計測空間を低空間周波数領域904〜907と、その外側の高空間周波数領域に分け、高空間周波数領域でエコー信号を取得する際の位相エンコードステップを低空間周波数領域で位相エンコードステップより大きくする。図9(a)は2倍に設定して、高空間周波数領域でのエコー信号の取得間隔を低空間周波数領域のエコー信号の取得間隔より2倍疎にして計測を行う例である。
【0063】
図9(a)のエコー信号群に対し、説明のために2つの群にエコーを分割する。
▲1▼低空間周波数領域904〜907内で取得したエコー信号群のみを選択したのが図9(b)
▲2▼ブロック内のエコー信号群のうち、1ラインおきのステップで取得したエコー信号群を選択したのが図9(c)
である。図9(b)は、低空間周波数領域のデータであるので、各ブロック701〜704でパラレルMRI法の再構成に必要なコイルの感度分布データ101〜104を算出することが可能である。また、図9(c)は、図3(c)と等価なデータであるので、各ブロック701〜704毎の折り返し画像305〜308を作成できる。このようにして作成したブロック毎の感度分布データ101〜104と折り返し画像305〜308を用いて図1と同様の処理を行い、最終画像を得ることができる。
【0064】
本発明のメリットは、パラレルMRI法で必要なRF受信コイルの感度分布データを、画像取得のための撮影(本計測)と同時に取得できることであり、専用パルスシーケンスで感度分布データを取得する場合と比較して、撮影中に被検体の位置が変化しても、感度分布データと本計測データ間で誤差が少なくて済み、安定して折り返しの展開が行えることである。心臓領域やその他被検体の動きの大きな部位を撮影する場合や、リアルタイムに撮影断面を変更しながら撮影を行う手法(インタラクティブスキャン法)等へ応用できる。
【0065】
プロペラMRI法では、ブロック毎に計測空間の回転角が異なるが、前述したとおり計測空間の回転と画像空間(画像データを表現する空間、計測空間をフーリエ変換によって写像した空間)の回転は対応する。そこで、あるブロックで取得した感度分布データを補間処理により画像空間上で回転することで座標系を変換し、他のブロックの座標系での再構成に使用することも可能である。図1を例にして説明すると、例えば、ブロック701に対する第1の座標系の感度分布データ101のみを計測しておき、残りのブロック702〜704については、感度分布データ取得のための計測は行わず、第1の感度分布データ101を用いて補間による回転処理を行い、それぞれのブロックの感度分布データ102〜104を作成する。このように処理することで、感度分布データ取得のための計測に必要なブロック数が低減し、撮影全体の時間を短縮することができる。
【0066】
あるいは、全てのブロックで感度分布データ用のエコー信号群を取得した場合、取得した感度分布データ用のエコー信号群をグリッディング処理によって一旦画像用の座標系に変換し、信号結合した後にフーリエ変換して感度分布データを作成し、回転処理により各ブロックに応じた座標系に再度変換し、パラレルMRI法の折り返しの展開に用いることもできる。この場合、感度分布データ用のエコー信号群が加算されるので、感度分布データのS/N比が向上する。パラレルMRI法では、折り返し展開に用いた感度分布データのS/N比は、結果の画像のS/N比に大きく影響することが知られており、感度分布のS/N比の向上は有用と考えられる。
【0067】
さらに、プロペラMRI法では、ブロック毎に計測空間の低空間周波数領域を取得しており、また、各ブロックの回転角度が異なることから、図9(c)のように間引いて計測を行っても、計測空間の低空間周波数領域は密に取得できる。このため、各ブロックで取得したエコー信号群をグリッディング処理によって一旦画像用の座標系に変換し、その低空間周波数領域を用いて感度分布データを作成することができる。この場合、感度分布データ取得用の専用パルスシーケンスを用いずに、画像取得のためのパルスシーケンスのみでパラレルMRI法の再構成を行うことができる。
【0068】
以上までが、本発明の各請求項を具体的に説明する実施例である。
しかし、本発明は、以上の実施例で開示された内容にとどまらず、本発明の趣旨を踏まえた上で各種形態を取り得る。本実施例ではグラディエントエコーパルスシーケンスについて記載したが、プロペラMRI法はパルスシーケンスの種類には依存せず、SEパルスシーケンス,FSEパルスシーケンス,EPIパルスシーケンスなどに適用できる。
【0069】
また、プロペラMRI法の実施例として、二次元面内のGr軸,Gp軸で計測空間の回転を行う場合について説明したが、Gr軸,Gp軸は撮影空間のX,Y,Zの任意の軸と対応させることが可能であり、オブリーク撮影やオフセンター撮影も実行できる。更に、三次元球内での回転を行うこともできる。
【0070】
また、ブロック数が4のプロペラMRI法の場合を説明したが、実際の撮影では、ブロック数とブロック内のエコー数は任意に設定することができ、この場合も同様の処理を行う。
【0071】
また、パラレルイメージングとして、RF受信コイルの数が2、位相エンコード間隔が2倍の場合を説明したが、RF受信コイルの数及び位相エンコード間隔の倍率は任意に設定できる(但し、RF受信コイルの数≧位相エンコード間隔の倍率)。また撮影する断面に併せてRF受信コイルで取得した信号を選択したり、合成して使用することもできる。
【0072】
また、本実施例では画像空間で折り返しの除去を行うパラレルMRI法について説明したが、計測空間で間引いた分のエコー信号を作成するパラレルMRI法にも適用可能である。
【0073】
【発明の効果】
本発明は、プロペラMRI法にパラレルMRI法を適用することにより、1つのブロックで取得するエコー数を低減しても折り返しアーチファクトを生ずることなく撮影時間を短縮化でき、かつ、パラレルMRI法を適用する際に必要な折り返しアーチファクト除去のための一般化アルゴリズムの適用や、プロペラMRI法に必要なグリッディング処理を行う必要が無くなるため、演算量を低減して高速演算が可能になり、全体として高速イメージングを行うことができる。
【図面の簡単な説明】
【図1】本発明の信号処理を説明する図。
【図2】パラレルMRI法を説明する図。
【図3】本発明に使用するデータを説明する図。
【図4】本発明が適用されるMRI装置の全体構成を示す図。
【図5】一般的なグラディエントエコーのパルスシーケンスを説明する図。
【図6】ラディアルMRI法のグラディエントエコーのパルスシーケンスを説明する図。
【図7】本発明が適用されるプロペラMRI法のパルスシーケンスを説明する図。
【図8】サンプリング方法と計測空間を説明する図。
【図9】本発明に使用するデータを説明する図。
【図10】本発明の信号処理を説明する図。
【図11】本発明の連続撮影を行う際の信号処理を説明する図。
【符号の説明】
1 被検体、2 静磁場発生系、3 傾斜磁場発生系、4 シーケンサ、5 送信系、6 受信系、7 信号処理系、8 中央演算装置(CPU)、9 傾斜磁場コイル、10 傾斜磁場電源、11 高周波発信器、12 変調器、13 高周波増幅器、14a 高周波照射コイル、14b 高周波受信コイル、15 増幅器、16 直交位相検波器、17 A/D変換器、18 磁気ディスク、19 光ディスク、20 ディスプレイ、501 高周波パルス、502 スライス選択傾斜磁場、503 位相エンコード傾斜磁場パルス、504 周波数エンコード傾斜磁場パルス、505 データサンプルウインド、506 エコー信号、507繰り返し時間間隔、508 画像取得時間、401 被検体、402 磁石、403 傾斜磁場コイル、404 RFコイル、405 RFプローブ、406 信号検出部、407 信号処理部、408 表示部、409 傾斜磁場電源、410 RF送信部、411 制御部、412 ベッド
Claims (16)
- 被検体から発生する核磁気共鳴信号を受信する受信手段と、
k空間の原点を含み低空間周波数領域から高空間周波数領域に亘って所定の幅を有して
成る単位領域に対応するエコー信号を計測する単位計測を、前記k空間の原点を中心にし
て前記単位領域の回転角度を変えて繰り返す計測制御手段と、
前記単位領域の各々に対応するエコー信号を用いて前記被検体の画像を再構成する信号処理手段と、を有する磁気共鳴イメージング装置において、
前記受信手段は、複数の受信コイルを組み合わせて成るマルチプル受信コイルを有し、
前記計測制御手段は、前記受信コイル毎にエコー信号を計測すると共に、少なくとも一
つの前記単位計測においてエコー信号の計測を間引き、
前記信号処理手段は、前記エコー信号の計測を間引かれた単位領域を含む複数の単位領域の各々に対応するエコー信号と前記複数の受信コイルの感度分布データとを用いて、折り返しアーチフアクトの無い被検体画像を再構成することを特徴とする磁気共鳴イメージング装置。 - 請求項1に記載の磁気共鳴イメージング装置において、
前記信号処理手段は、前記単位領域毎にそのエコー信号を用いて該単位領域の画像をそれぞれ再構成する単位領域画像生成と、複数の単位領域画像から前記被検体画像を再構成する被検体画像生成と、を行うことを特徴とする磁気共鳴イメージング装置。 - 請求項2に記載の磁気共鳴イメージング装置において、
前記信号処理手段は、前記エコー信号の計測が間引かれた単位領域において、前記受信コイル毎のエコー信号と感度分布データとを用いて、折り返しアーチフアクトを発生させない演算を施して該単位領域の画像を再構成することを特徴とする磁気共鳴イメージング装置。 - 請求項1乃至3のいずれか一項に記載の磁気共鳴イメージング装置において、前記計測制御手段は、
前記単位領域に対応するエコー信号の計測を、複数の平行軌跡に沿って行い、
前記間引き計測を、前記平行軌跡のうちの少なくとも一つの軌跡に対応するエコー信号の計測を間引いて行うことを特徴とする磁気共鳴イメージング装置。 - 請求項1乃至4のいずれか一項に記載の磁気共鳴イメージング装置において、前記計測
制御手段は、
前記間引き計測を、前記単位領域の高空間周波数領域側に対応するエコー信号の計測を間引いて行うことを特徴とする磁気共鳴イメージング装置。 - 請求項1乃至5のいずれか一項に記載の磁気共鳴イメージング装置において、
前記信号処理手段は、前記エコー信号の計測が間引かれた単位領域では、該単位領域に対応するエコー信号の少なくとも一部を用いて前記受信コイル毎の感度分布データを作成することを特徴とする磁気共鳴イメージング装置。 - 請求項6に記載の磁気共鳴イメージング装置において、
前記信号処理手段は、前記エコー信号の計測が間引かれた単位領域では、該単位領域の低空間周波数領域における受信コイル毎のエコー信号を用いて、前記受信コイル毎の感度分布データをそれぞれ作成することを特徴とする磁気共鳴イメージング装置。 - 請求項4に記載の磁気共鳴イメージング装置において、前記エコー信号の計測が間引かれた単位領域では、
前記計測制御手段は、前記複数の平行軌跡の低空間周波数領域における間隔を高空間周波数領域のおける間隔よりも密にして、前記受信コイル毎に前記エコー信号を計測し、
前記信号処理手段は、前記密に計測された低空間周波数領域における受信コイル毎のエコー信号を用いて、前記受信コイル毎の感度分布データをそれぞれ作成することを特徴とする磁気共鳴イメージング装置。 - 請求項8に記載の磁気共鳴イメージング装置において、
前記信号処理手段は、前記エコー信号の計測が間引かれた単位領域で前記密に計測された低空間周波数領域における受信コイル毎のエコー信号を用いて、他の単位領域における前記受信コイル毎の感度分布データをそれぞれ作成することを特徴とする磁気共鳴イメージング装置。 - 請求項1に記載の磁気共鳴イメージング装置において、前記エコー信号の計測が間引かれた単位領域では、
前記信号処理手段は、事前計測された受信コイル毎の感度分布データ用エコー信号を用いて、前記受信コイル毎の感度分布データをそれぞれ作成することを特徴とする磁気共鳴イメージング装置。 - 請求項1に記載の磁気共鳴イメージング装置において、
前記信号処理手段は、前記単位領域計測が終了する度に、その単位領域画像を再構成すると共に、該単位領域画像と他の少なくとも一つの単位領域画像とから前記被検体画像を再構成することを特徴とする磁気共鳴イメージング装置。 - 請求項1に記載の磁気共鳴イメージング装置において、
前記計測制御手段は、少なくとも1つの前記単位領域において受信コイル毎の感度分布データ用エコー信号を計測し、
前記信号処理手段は、少なくとも1つの前記単位領域における前記受信コイル毎の感度分布データ用エコー信号を合成して求めた受信コイル毎の感度分布データに基づいて、各単位領域における前記受信コイル毎の感度分布データを求めることを特徴とする磁気共鳴イメージング装置。 - 請求項1に記載の磁気共鳴イメージング装置において、
前記信号処理手段は、前記受信コイル毎に、複数の前記間引き計測された単位領域における低空間周波数領域のデータを合成して求めた感度分布データ用エコー信号に基づいて
、各単位領域における感度分布データを作成することを特徴とする磁気共鳴イメージング装置。 - 請求項1に記載の磁気共鳴イメージング装置において、
前記信号処理手段は、一つの単位領域で求めた前記受信コイル毎の感度分布データを用いて、他の単位領域における前記受信コイル毎の感度分布データを作成することを特徴とする磁気共鳴イメージング装置。 - 請求項1に記載の磁気共鳴イメージング装置において、
前記信号処理手段は、前記単位領域計測が終了する度に、該単位領域のエコー信号と他の単位領域のエコー信号とを用いて前記被検体画像を再構成することを特徴とする磁気共鳴イメージング装置。 - 請求項2に記載の磁気共鳴イメージング装置において、
前記信号処理手段は、前記単位領域画像の各々を逆フーリエ変換して求めたk空間データをグリッディングにより一つのk空間データに合成し、該合成された一つのk空間デー タを用いて前記被検体画像を再構成することを特徴とする磁気共鳴イメージング装置。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003119403A JP3929047B2 (ja) | 2003-04-24 | 2003-04-24 | 磁気共鳴イメージング装置 |
US10/553,900 US7372269B2 (en) | 2003-04-24 | 2004-04-23 | Magnetic resonance imaging method and apparatus |
EP04729259A EP1618843A4 (en) | 2003-04-24 | 2004-04-23 | MAGNETIC RESONANCE TOMOGRAPHY PROCESS AND DEVICE |
PCT/JP2004/005928 WO2004093682A1 (ja) | 2003-04-24 | 2004-04-23 | 磁気共鳴イメージング方法及び装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003119403A JP3929047B2 (ja) | 2003-04-24 | 2003-04-24 | 磁気共鳴イメージング装置 |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2004344183A JP2004344183A (ja) | 2004-12-09 |
JP2004344183A5 JP2004344183A5 (ja) | 2006-06-01 |
JP3929047B2 true JP3929047B2 (ja) | 2007-06-13 |
Family
ID=33308100
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003119403A Expired - Fee Related JP3929047B2 (ja) | 2003-04-24 | 2003-04-24 | 磁気共鳴イメージング装置 |
Country Status (4)
Country | Link |
---|---|
US (1) | US7372269B2 (ja) |
EP (1) | EP1618843A4 (ja) |
JP (1) | JP3929047B2 (ja) |
WO (1) | WO2004093682A1 (ja) |
Families Citing this family (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4443079B2 (ja) | 2001-09-13 | 2010-03-31 | 株式会社日立メディコ | 磁気共鳴イメージング装置及び磁気共鳴イメージング装置用rf受信コイル |
US7450982B2 (en) * | 2003-03-14 | 2008-11-11 | Hitachi Medical Corporation | Magnetic resonance imaging system and method |
CN100392424C (zh) * | 2004-11-15 | 2008-06-04 | 华东师范大学 | 一种用于图形化脉冲序列编译器中实现回波数据重组的方法 |
US7023207B1 (en) * | 2005-02-16 | 2006-04-04 | General Electric Company | Method and system of MR imaging with reduced radial ripple artifacts |
DE102005046732B4 (de) * | 2005-04-18 | 2010-04-15 | Siemens Ag | Verbessertes Rekonstruktionsverfahren bei der Propellerbildgebung in der Magnetresonanztomographie |
JP2006304955A (ja) * | 2005-04-27 | 2006-11-09 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴イメージング装置 |
US7265547B2 (en) * | 2005-09-16 | 2007-09-04 | General Electric Company | Method and apparatus for acquiring MR data with a segmented multi-shot radial fan beam encoding order |
KR100726739B1 (ko) * | 2006-01-16 | 2007-06-11 | 한국과학기술원 | 방사형 케이-공간 경로법을 이용한 병렬 자기 공명 영상획득 방법, 그 장치 및 이를 실행하기 위한 컴퓨터로 읽을수 있는 기록매체. |
DE102006033862B3 (de) * | 2006-07-21 | 2007-12-06 | Siemens Ag | Verfahren zur dynamischen Magnet-Resonanz-Bildgebung sowie Magnet-Resonanz-Gerät |
US7482806B2 (en) * | 2006-12-05 | 2009-01-27 | Siemens Aktiengesellschaft | Multi-coil magnetic resonance data acquisition and image reconstruction method and apparatus using blade-like k-space sampling |
BRPI0720535A2 (pt) * | 2006-12-19 | 2014-01-07 | Koninkl Philips Electronics Nv | Arranjo de formação de imagem de diagnóstico, e, arranjo terapêutico |
US7535222B2 (en) * | 2007-01-02 | 2009-05-19 | The Board Of Trustees Of The Leland Stanford Junior University | MRI data acquisition using propeller k-space data acquisition |
US7663364B2 (en) * | 2007-02-06 | 2010-02-16 | Kabushiki Kaisha Toshiba | Magnetic resonance imaging apparatus, image processing apparatus, computer program product, and data storing method |
US8207734B2 (en) | 2007-03-09 | 2012-06-26 | Hitachi Medical Corporation | Parallel imaging in non-cartesian magnetic resonance imaging (MRI) and MRI apparatus realizing the same |
US8314614B2 (en) * | 2007-05-03 | 2012-11-20 | Koninklijke Philips Electronics N.V. | PROPELLER MRI with phase correction |
CN101470180B (zh) * | 2007-12-29 | 2016-01-20 | 西门子(中国)有限公司 | 磁共振成像中失真校准的方法和装置 |
WO2009093517A1 (ja) * | 2008-01-23 | 2009-07-30 | Hitachi Medical Corporation | 磁気共鳴イメージング装置及びマルチコントラスト画像取得方法 |
US8022700B2 (en) * | 2008-11-07 | 2011-09-20 | General Electric Company | Method and apparatus for view ordering of magnetic resonance imaging data for dynamic studies |
WO2013047275A1 (ja) * | 2011-09-29 | 2013-04-04 | 株式会社 日立メディコ | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 |
US10534056B2 (en) * | 2012-06-27 | 2020-01-14 | Siemens Healthcare Gmbh | System for simultaneous dual-slab acquisition of MR images with asymmetric and time-reversed asymmetric, concatenated pulses |
CN103076583B (zh) * | 2012-12-28 | 2015-04-15 | 深圳先进技术研究院 | 磁共振快速成像方法和系统 |
CN103646410B (zh) * | 2013-11-27 | 2016-06-08 | 中国科学院深圳先进技术研究院 | 磁共振快速参数成像方法和系统 |
RU2016141324A (ru) * | 2014-03-24 | 2018-05-03 | Конинклейке Филипс Н.В. | Магнитно-резонансная визуализация по технологии propeller |
US20170307716A1 (en) * | 2014-10-10 | 2017-10-26 | Koninklijke Philips N.V. | Propeller mr imaging with artefact suppression |
US9983283B2 (en) | 2015-03-16 | 2018-05-29 | Toshiba Medical Systems Corporation | Accelerated MRI using radial strips and undersampling of k-space |
CN106842084B (zh) * | 2016-12-30 | 2019-11-12 | 上海联影医疗科技有限公司 | 一种磁共振成像方法及装置 |
DE102017203936A1 (de) * | 2017-03-09 | 2018-09-13 | Siemens Healthcare Gmbh | Leise Echo Planar Bildgebung |
US10809341B1 (en) * | 2019-04-19 | 2020-10-20 | Canon Medical Systems Corporation | Readout-segmented echo planar imaging with k-space averaging |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5243284A (en) * | 1991-07-24 | 1993-09-07 | Board Of Trustees Of The Leland Stanford Junior University | Method of magnetic resonance reconstruction imaging from projections using partial data collected in k-space |
US5833609A (en) | 1996-11-26 | 1998-11-10 | Picker International, Inc. | Rotating diffusion MR imaging reduced motion artifacts |
JP3952247B2 (ja) * | 1999-12-08 | 2007-08-01 | 株式会社日立メディコ | 核磁気共鳴撮影装置 |
US6411089B1 (en) * | 2000-11-22 | 2002-06-25 | Philips Medical Systems (Cleveland), Inc. | Two-dimensional phase-conjugate symmetry reconstruction for 3d spin-warp, echo-planar and echo-volume magnetic resonance imaging |
US6556009B2 (en) * | 2000-12-11 | 2003-04-29 | The United States Of America As Represented By The Department Of Health And Human Services | Accelerated magnetic resonance imaging using a parallel spatial filter |
JP4047553B2 (ja) * | 2001-04-20 | 2008-02-13 | 株式会社日立メディコ | 磁気共鳴イメージング装置 |
DE10119660B4 (de) * | 2001-04-20 | 2006-01-05 | Siemens Ag | Verfahren zur schnellen Gewinnung eines Magnetresonanzbildes |
JP4109841B2 (ja) | 2001-06-19 | 2008-07-02 | 株式会社東芝 | 半導体集積回路装置および半導体機器システム |
JP2003315731A (ja) | 2002-04-25 | 2003-11-06 | Fuji Photo Film Co Ltd | プロジェクタ表示装置 |
US6882148B2 (en) * | 2003-07-09 | 2005-04-19 | Catholic Healthcare West | Split-blade data collection for propeller MRI |
US7102348B2 (en) * | 2004-08-05 | 2006-09-05 | Siemens Aktiengesellschaft | MRI method and apparatus for faster data acquisition or better motion artifact reduction |
US7023207B1 (en) * | 2005-02-16 | 2006-04-04 | General Electric Company | Method and system of MR imaging with reduced radial ripple artifacts |
-
2003
- 2003-04-24 JP JP2003119403A patent/JP3929047B2/ja not_active Expired - Fee Related
-
2004
- 2004-04-23 EP EP04729259A patent/EP1618843A4/en not_active Withdrawn
- 2004-04-23 WO PCT/JP2004/005928 patent/WO2004093682A1/ja active Application Filing
- 2004-04-23 US US10/553,900 patent/US7372269B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
JP2004344183A (ja) | 2004-12-09 |
WO2004093682A1 (ja) | 2004-11-04 |
US7372269B2 (en) | 2008-05-13 |
US20060232273A1 (en) | 2006-10-19 |
EP1618843A4 (en) | 2011-05-25 |
EP1618843A1 (en) | 2006-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3929047B2 (ja) | 磁気共鳴イメージング装置 | |
JP6513398B2 (ja) | 事前情報に制約される正則化を使用するmr画像再構成 | |
US8154294B2 (en) | Magnetic resonance imaging apparatus and method for correcting error due to gradient magnetic field | |
US8344729B2 (en) | Magnetic resonance imaging apparatus | |
EP1444530B1 (en) | Three-dimensional phase contrast magnetic resonance imaging using interleaved projection-reconstruction data | |
CN103608693B (zh) | 具有用于脂肪抑制的化学位移编码的对比度增强磁共振血管造影 | |
JP5221570B2 (ja) | 磁気共鳴イメージング装置及びマルチコントラスト画像取得方法 | |
US9170313B2 (en) | Coronary magnetic resonance angiography with signal separation for water and fat | |
WO2010116772A1 (ja) | 磁気共鳴イメージング装置 | |
EP2350698B1 (en) | System and method for moving table mri | |
JPWO2003092497A1 (ja) | 磁気共鳴イメージング装置 | |
US6611143B2 (en) | MRI apparatus generating odd and/or even echo images with sensitivity distribution of coils | |
JP2017529960A (ja) | アーチファクト抑制を有するプロペラmrイメージング | |
JP2021512671A (ja) | 固有運動補正を有するスタック・オブ・スター収集を用いたmr撮像 | |
US11543482B2 (en) | Magnetic resonance imaging using motion-compensated image reconstruction | |
WO2005023108A1 (ja) | 磁気共鳴イメージング装置 | |
US11009576B2 (en) | Method for magnetic resonance imaging using slice quadratic phase for spatiotemporal encoding | |
CN111344588A (zh) | 使用运动相关径向或螺旋k空间采样进行MR成像 | |
JP4817381B2 (ja) | 磁気共鳴イメージング装置 | |
US20070268019A1 (en) | Imaging Procedure and Magnetic-Resonance Imaging System for the Acquistion of the Longitudinal Spin-Lattice Relaxation Time | |
JP2022045474A (ja) | 磁気共鳴撮像装置、及び、画像生成方法 | |
CN116724246A (zh) | Mr乳房摄影 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060329 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060329 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20061121 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20061226 |
|
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: 20070129 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070130 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070305 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 3929047 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110316 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110316 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130316 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140316 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 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |