JP2007050278A - Nuclear magnetic resonance photographing device - Google Patents
Nuclear magnetic resonance photographing device Download PDFInfo
- Publication number
- JP2007050278A JP2007050278A JP2006292684A JP2006292684A JP2007050278A JP 2007050278 A JP2007050278 A JP 2007050278A JP 2006292684 A JP2006292684 A JP 2006292684A JP 2006292684 A JP2006292684 A JP 2006292684A JP 2007050278 A JP2007050278 A JP 2007050278A
- Authority
- JP
- Japan
- Prior art keywords
- magnetic field
- gradient magnetic
- image
- echo signal
- phase encoding
- 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.)
- Granted
Links
- 238000005481 NMR spectroscopy Methods 0.000 title abstract description 5
- 230000035945 sensitivity Effects 0.000 claims abstract description 69
- 238000009826 distribution Methods 0.000 claims abstract description 65
- 230000005415 magnetization Effects 0.000 claims abstract description 28
- 238000000034 method Methods 0.000 claims abstract description 24
- 238000001514 detection method Methods 0.000 claims abstract description 10
- 238000005259 measurement Methods 0.000 claims description 36
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 claims description 8
- 230000001678 irradiating effect Effects 0.000 claims description 6
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 claims 12
- 238000003745 diagnosis Methods 0.000 abstract description 2
- 230000006866 deterioration Effects 0.000 abstract 1
- 238000003384 imaging method Methods 0.000 description 18
- 238000012545 processing Methods 0.000 description 12
- 230000015572 biosynthetic process Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 9
- 238000003786 synthesis reaction Methods 0.000 description 9
- 230000008859 change Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000004904 shortening Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000002194 synthesizing effect Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 238000012307 MRI technique Methods 0.000 description 1
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 1
- 230000001771 impaired effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052698 phosphorus Inorganic materials 0.000 description 1
- 239000011574 phosphorus Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000000264 spin echo pulse sequence Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Description
本発明は、連続的に被検体中の水素や燐等からの核磁気共鳴(以下、「NMR」という)信号を測定し、核の密度分布や緩和時間分布等を映像化する核磁気共鳴撮影(MRI)装置に関する。 The present invention continuously measures nuclear magnetic resonance (hereinafter referred to as "NMR") signals from hydrogen, phosphorus, etc. in a subject and visualizes nuclear density distribution, relaxation time distribution, etc. (MRI) relates to equipment.
MRI装置において被検体から発生するNMR信号を検出する受信コイルとして、「マルチプルRFコイル」もしくは「フェイズドアレイコイル」と呼ばれる高感度コイルが近年多用され始めている(特許文献1など)。マルチプルRFコイルは、相対的に高感度な小型RFコイルを複数個並べて、各RFコイルで受信した信号を合成することにより、小型RFコイルの高い感度を保ったまま視野を拡大し高感度化を図る受信専用RFコイルであり、静磁場方式や検出部位に応じて種々のものが提案されている。
In recent years, high-sensitivity coils called “multiple RF coils” or “phased array coils” have begun to be frequently used as receiving coils for detecting NMR signals generated from a subject in an MRI apparatus (
一方、近年、マルチプルコイルを用いて、位相エンコード方向のデータを間引くことにより撮影時間を短縮する手法が提案されている(例えば、非特許文献1、非特許文献2など)。このような技術は、空間エンコード法、もしくはパラレルMRIと呼ばれ、マルチプルRFコイルの感度分布が互いに空間的に異なることを用いて位相エンコードデータを間引いた際の折り返しを除去する。この除去には、高精度なRFコイルの感度分布を使った高精度な演算が必要である。上掲の非特許文献1に記載された方法では、この演算を計測空間(k空間)で行い、非特許文献2に記載された方法では、演算をフーリエ変換後の実空間で行う。
On the other hand, in recent years, methods have been proposed for shortening the imaging time by thinning out data in the phase encoding direction using multiple coils (for example, Non-Patent
一般にRFコイルの感度分布は、各RF受信信号から求めることができ、具体的には、事前に均一濃度のファントムを投影し画像の空間的シェーディングがRFコイルの感度分布であるとする方法、被検体を別途撮影した画像に低周波フィルタを作用させ計算する方法などが知られている。
上記RFコイルの感度分布を求める処理は、通常、撮影に先立ち行われる。従って、例えば腹部では呼吸により被検体の形状が時間的に変動する場合、感度分布を計算した時刻と本撮影を行うときで被検体形状が異なることがある。また、これに伴い被検体に密着したRFコイルの空間的配置が変化することも有りうる。また、手術中に撮影を行うときは、時々刻々とRFコイルの位置が変わることもある。 The processing for obtaining the sensitivity distribution of the RF coil is usually performed prior to imaging. Therefore, for example, when the shape of the subject fluctuates with time due to breathing in the abdomen, the shape of the subject may differ between the time when the sensitivity distribution is calculated and the time of actual imaging. As a result, the spatial arrangement of the RF coil that is in close contact with the subject may change. Also, when taking an image during surgery, the position of the RF coil may change from moment to moment.
このように撮影した画像をリアルタイムで表示することが要求されるような撮影では、事前に取得された感度データを用いる従来方法は、誤差を招き画質が劣化する。また、撮影に先立ち事前にRFコイルの感度分布を計測することは、トータルの撮影時間を延長することになり、本技術の特徴である短時間撮影の効果を低下させる。 In such shooting that requires real-time display of a shot image, the conventional method using sensitivity data acquired in advance introduces an error and deteriorates image quality. In addition, measuring the sensitivity distribution of the RF coil in advance prior to shooting extends the total shooting time and reduces the effect of short-time shooting, which is a feature of the present technology.
上記課題を解決するために、本発明は、特定の計測データ収集の制御(計測空間(k空間)が領域によって異なる密度となるような制御)を行うために好適なMRI装置を提供する。すなわち本発明のMRI装置は、検知する磁化を含む被検体に高周波パルスを照射し横磁化を発生する手段と、横磁化が付与された被検体に位相エンコード傾斜磁場パルスを印加する手段と、読み出し傾斜磁場パルスを印加する手段と、前記読み出し傾斜磁場の印加中に発生するエコー信号を検出する手段と、前記各手段を制御する制御手段と、エコー信号から画像を再構成する画像再構成手段を含み、
前記制御手段は、前記高周波パルスによる横磁化発生に続き、極性の反転する読み出し傾斜磁場を連続して印加するとともに読み出し傾斜磁場の印加に同期して位相エンコード傾斜磁場を印加し、その際、連続するエコー信号計測の所定の時間領域では前記位相エンコード傾斜磁場の強度を第1の値に設定し、前記所定の時間領域以外の時間領域では前記位相エンコード傾斜磁場強度を第2の値に設定することを特徴とする(請求項1)。
In order to solve the above-described problems, the present invention provides an MRI apparatus suitable for performing control of specific measurement data collection (control so that the measurement space (k space) has different densities depending on regions). That is, the MRI apparatus of the present invention includes means for irradiating a subject including magnetization to be detected with a high frequency pulse to generate transverse magnetization, means for applying a phase encoding gradient magnetic field pulse to a subject to which transverse magnetization is applied, and reading. A means for applying a gradient magnetic field pulse; a means for detecting an echo signal generated during application of the readout gradient magnetic field; a control means for controlling each of the means; and an image reconstruction means for reconstructing an image from the echo signal. Including
The control means continuously applies a read gradient magnetic field whose polarity is reversed following the generation of transverse magnetization by the high-frequency pulse and applies a phase encode gradient magnetic field in synchronization with the application of the read gradient magnetic field. The intensity of the phase encoding gradient magnetic field is set to a first value in a predetermined time region of echo signal measurement, and the phase encoding gradient magnetic field strength is set to a second value in a time region other than the predetermined time region. (Claim 1).
また本発明のMRI装置は、検知する磁化を含む被検体に高周波パルスを照射し横磁化を発生する手段と、横磁化が付与された被検体に位相エンコード傾斜磁場パルスを印加する手段と、読み出し傾斜磁場パルスを印加する手段と、前記読み出し傾斜磁場の印加中に発生するエコー信号を検出する手段と、前記各手段を制御する制御手段と、エコー信号から画像を再構成する画像再構成手段を含み、
前記エコー信号を検出する手段は、互いに空間的に異なる検出感度分布を有する複数の受信コイルを備え、
前記制御手段は、前記高周波パルスによる横磁化発生に続き、極性の反転する読み出し傾斜磁場を連続して印加するとともに読み出し傾斜磁場の印加に同期して位相エンコード傾斜磁場を印加し、エコー信号を連続して計測し、
前記画像再構成手段は、前記連続して計測したエコー信号の一部を用いて、各受信コイルの感度分布を計算し、計算した感度分布と前記連続計測したエコー信号とを用いて画像を再構成することを特徴とする(請求項2)。
The MRI apparatus of the present invention also includes means for irradiating a subject including magnetization to be detected with a radio frequency pulse to generate transverse magnetization, means for applying a phase encoding gradient magnetic field pulse to a subject to which transverse magnetization has been applied, and readout. A means for applying a gradient magnetic field pulse; a means for detecting an echo signal generated during application of the readout gradient magnetic field; a control means for controlling each of the means; and an image reconstruction means for reconstructing an image from the echo signal. Including
The means for detecting the echo signal includes a plurality of receiving coils having spatially different detection sensitivity distributions,
The control means continuously applies a read gradient magnetic field whose polarity is reversed following the generation of transverse magnetization by the high-frequency pulse, and applies a phase encode gradient magnetic field in synchronization with the application of the read gradient magnetic field, and continuously echo signals. And measure
The image reconstruction means calculates a sensitivity distribution of each receiving coil using a part of the continuously measured echo signals, and reconstructs an image using the calculated sensitivity distribution and the continuously measured echo signals. It is characterized by comprising (claim 2).
また本発明のMRI装置は、検知する磁化を含む被検体に高周波パルスを照射し横磁化を発生する手段と、横磁化が付与された被検体に位相エンコード傾斜磁場パルスを印加する手段と、読み出し傾斜磁場パルスを印加する手段と、前記読み出し傾斜磁場の印加中に発生するエコー信号を検出する手段と、前記各手段を制御する制御手段と、エコー信号から画像を再構成する画像再構成手段を含み、
前記エコー信号を検出する手段は、互いに空間的に異なる検出感度分布を有する複数の受信コイルを備え、
前記制御手段は、前記高周波パルスによる横磁化発生に続き、極性の反転する読み出し傾斜磁場を連続して印加するとともに読み出し傾斜磁場の印加に同期して位相エンコード傾斜磁場を印加し、その際、前記位相エンコード傾斜磁場の強度を異ならせて、エコー信号を連続して計測し、
前記画像再構成手段は、前記連続して計測したエコー信号の一部を用いて、各受信コイルの感度分布を計算し、計算した感度分布と前記連続計測したエコー信号とを用いて画像を再構成することを特徴とする(請求項3)。
The MRI apparatus of the present invention also includes means for irradiating a subject including magnetization to be detected with a radio frequency pulse to generate transverse magnetization, means for applying a phase encoding gradient magnetic field pulse to a subject to which transverse magnetization has been applied, and readout. A means for applying a gradient magnetic field pulse; a means for detecting an echo signal generated during application of the readout gradient magnetic field; a control means for controlling each of the means; and an image reconstruction means for reconstructing an image from the echo signal. Including
The means for detecting the echo signal includes a plurality of receiving coils having spatially different detection sensitivity distributions,
The control means continuously applies a read gradient magnetic field whose polarity is reversed following the generation of transverse magnetization by the high-frequency pulse and applies a phase encode gradient magnetic field in synchronization with the application of the read gradient magnetic field. Measure echo signals continuously with different phase encoding gradient magnetic field strengths,
The image reconstruction means calculates a sensitivity distribution of each receiving coil using a part of the continuously measured echo signals, and reconstructs an image using the calculated sensitivity distribution and the continuously measured echo signals. It is characterized by comprising (claim 3).
本発明の制御手段が行う撮影方法は、例えば、エコープレナーイメージング(EPI)であり、EPIにおいて連続して印加する位相エンコードの値を時間領域によって値を変えることにより、k空間における計測データの配列を一部の領域について粗にすることができる。 The imaging method performed by the control means of the present invention is, for example, echo planar imaging (EPI). By changing the value of phase encoding applied continuously in EPI according to the time domain, the arrangement of measurement data in k space Can be roughened for some areas.
本発明のMRI装置では、計算に必要なRFコイルの感度分布を本計測中の取得したデータから計算するため、感度分布のデータと計測データに時刻の差が無い。したがって、時々刻々と状況が変化するMRI撮影でも安定して、折り返し除去が可能になる。すなわちリアルタイム性が要求される撮影でも、誤差を招かず画質が劣化しない。 In the MRI apparatus of the present invention, since the sensitivity distribution of the RF coil necessary for the calculation is calculated from the acquired data during the main measurement, there is no time difference between the sensitivity distribution data and the measurement data. Therefore, even with MRI imaging in which the situation changes from moment to moment, the aliasing can be removed stably. That is, even when shooting is required to have real-time properties, no error is caused and the image quality does not deteriorate.
また、投影に先立ち事前にRFコイルの感度分布を計測しないので、トータルの撮影時間を延長することにならず、短時間撮影の効果を損ねない。 Also, since the sensitivity distribution of the RF coil is not measured in advance prior to projection, the total imaging time is not extended and the effect of short-time imaging is not impaired.
以下、実施例を用いて本発明を説明する。図1は、本発明が適用される典型的なMRI装置の構成を示す図で、このMRI装置は、被検体101の周囲に静磁場を発生する磁石102と、該空間に傾斜磁場を発生する傾斜磁場コイル103と、この領域に高周波磁場を発生するRFコイル104と、被検体101が発生するMR信号を検出するRFコイル105と、被検体101が横たわるベッド112とを備えている。
Hereinafter, the present invention will be described using examples. FIG. 1 is a diagram showing a configuration of a typical MRI apparatus to which the present invention is applied. This MRI apparatus generates a
傾斜磁場コイル103は、X、Y、Zの3方向の傾斜磁場コイルで構成され、傾斜磁場電源109からの信号に応じてそれぞれ傾斜磁場を発生する。傾斜磁場の印加の仕方により被検体101の任意の断面を選択することができ、またMR信号に位置情報を与えることができる。MR信号に位置情報を与える傾斜磁場は、位相エンコード傾斜磁場、読み出し傾斜磁場と呼ばれ、これによって計測データが配置される計測空間(k空間)が規定される。
The gradient
RFコイル104はRF送信部110の信号に応じて高周波磁場を発生する。この高周波磁場の周波数は撮影対象である核スピンの共鳴周波数に同調されている。通常MRIの撮影対象は、被検体の主たる構成物質、プロトンである。RFコイル105の信号は、信号検出部106で検出され、信号処理部107で信号処理され、また計算により画像信号に変換される。
The
RFコイル104は図2に示すように小型RFコイル202を複数個(図示する例では、4個)並べたマルチプルコイル201である。各小型RFコイル202は、例えば隣接して配置された一辺15cmの矩形サーフェスコイルとすることができ、矩形コイルは撮影の位相エンコード方向に一部がオーバーラップするように配置され、互に区分された検出感度領域を有する。
As shown in FIG. 2, the
各小型RFコイル202は、それぞれプリアンプ203に接続されており、各コイルからの出力はプリアンプ203でそれぞれ増幅される。各プリアンプ203はそれぞれAD変換・直交検波回路204に並列に接続され、信号検出部106を構成する。AD変換直交検波器204で検波された信号は、各コイル毎のフーリエ変換、フィルタリング、合成演算などを行うために信号処理部107に送られる。信号処理部107で行う処理は、予めプログラムとして組込まれている。
Each small RF coil 202 is connected to a
図1において、傾斜磁場電源109、RF送信部110、信号検出部106は、一般にパルスシーケンスと呼ばれるタイムチャートに従い、制御部111で制御される。本発明のMRI装置では、制御部111が、k空間の一部の領域の位相エンコードが粗であるように制御する。
In FIG. 1, a gradient magnetic
次に上記構成のMRI装置を用いた撮影方法および画像再構成(合成)について図3〜図9を参照して説明する。
本発明のMRI装置は、パルスシーケンスとして、1回後に連続してエコー信号を計測するパルスシーケンスを採用している。図3はそのようなパルスシーケンスの一実施例を示す図で、スピンエコー型のEPIシーケンスを示している。即ち、スライス選択傾斜磁場パルスGs1103と同時に被検体組織の核スピンを励起するRFパルス1101を印加し、TE/2時間後に最初のRFパルス1101によって発生した横磁化を反転するRFパルス1102をスライス選択傾斜磁場パルスGs1104と共に印加し、その後極性が反転する読み出し傾斜磁場パルスGr1105を連続して印加しながらエコー信号1106を計測する。この際、読み出し傾斜磁場パルスGr毎にエコー信号を位相エンコードする位相エンコード傾斜磁場Ge1107、1108を印加する。
Next, an imaging method and image reconstruction (synthesis) using the MRI apparatus having the above configuration will be described with reference to FIGS.
The MRI apparatus of the present invention employs a pulse sequence that continuously measures echo signals after one time as a pulse sequence. FIG. 3 is a diagram showing an example of such a pulse sequence, and shows a spin echo type EPI sequence. That is, simultaneously with the slice selection gradient magnetic field pulse Gs1103, the
通常のEPIシーケンスでは、例えばk空間のマトリクスサイズが128×128であるとすると、読み出し傾斜磁場パルスを128回反転し、128個のエコー信号を計測する。位相エンコード傾斜磁場パルスも最初のオフセットパルスの他、同じ大きさのパルスを128個印加する。これに対し本実施例のEPIシーケンスでは、例えば位相エンコード0のエコー信号から所定の位相エンコードまでのエコー信号は位相エンコード数が1ずつインクリメントするような大きさの位相エンコード傾斜磁場パルス1107を印加する。それ以降に計測するエコー信号については位相エンコード数が複数、例えば4ずつインクリメントするような大きさの位相エンコード傾斜磁場パルス1108を印加する。
In a normal EPI sequence, for example, assuming that the matrix size of the k space is 128 × 128, the readout gradient magnetic field pulse is inverted 128 times and 128 echo signals are measured. For the phase encode gradient magnetic field pulse, 128 pulses of the same magnitude are applied in addition to the first offset pulse. On the other hand, in the EPI sequence of this embodiment, for example, a phase encode gradient
図4は、上記シーケンスの1回または複数回の繰り返しによって測定された測定データのk空間データ配列(kトラジェクトリ)を示す図であり、1回の信号取得で、k空間横方向(kx方向)の全列または複数列のデータが取得される。 FIG. 4 is a diagram showing a k-space data array (k trajectory) of measurement data measured by one or a plurality of repetitions of the above-described sequence. In one signal acquisition, k-space lateral direction (kx direction) All or multiple columns of data are acquired.
図3に示す例では簡単にするために8個のエコー信号を計測する場合を示しており、5番目の信号計測まで通常の大きさの位相エンコード傾斜磁場パルス1107を印加し、6番目から8番目の信号計測まで通常の大きさの4倍の位相エンコード傾斜磁場パルス1108を印加している。1〜5番目までの信号が図4に示すk空間の領域1(402)に対応し、6〜8番目までの信号が領域2(403)に対応する。従って6〜8番目について、従来は12(3×4)個の信号を計測する時間が必要であったのに対し、このシーケンスでは9(12−3)個の信号計測時間分短縮を図ることができる。kトラジェクトリで言えば、上記シーケンスの実行により領域1について密に計測をし、領域2については粗に計測することができる。
For the sake of simplicity, the example shown in FIG. 3 shows a case where eight echo signals are measured, and a phase-encoding gradient
一例として、k空間の中心から上下各10%のデータまでは、1ステップずつ位相エンコードは増えていくが、10%を超えると2ステップずつもしくは4ステップずつ位相エンコードが増えていく。この結果、全体の撮影時間は間引かれた分だけ短縮する。 As an example, the phase encoding increases by 1 step from the center of the k space up to 10% of the upper and lower data, but when it exceeds 10%, the phase encoding increases by 2 steps or 4 steps. As a result, the entire shooting time is shortened by the thinned-out amount.
領域1(402)がk空間全体に対し占める割合(ここでは約20%)は、コイル感度および撮影時間の短縮の両者を考慮して適当に設定する。コイル感度が急峻の場合には割合を増やしたほうが望ましく、一方、撮影時間の短縮という観点からは少ないほうが良い。kトラジェクトリは、典型的には256×256マトリクスであり、領域1(402)では、この密度でデータが収集される。領域2(403)では、位相エンコード方向に2倍から4倍粗くサンプリングされることになる。 The ratio of the area 1 (402) to the entire k space (about 20% in this case) is appropriately set in consideration of both the coil sensitivity and the reduction of the photographing time. When the coil sensitivity is steep, it is desirable to increase the ratio, while from the viewpoint of shortening the photographing time, it is preferable to decrease the ratio. The k trajectory is typically a 256 × 256 matrix, and data is collected at this density in region 1 (402). In the region 2 (403), the sampling is performed roughly 2 to 4 times in the phase encoding direction.
このようにEPIパルスシーケンスを用い領域によって位相エンコードステップの密度を異ならせて計測されたデータは、各受信コイル202毎に信号処理部107に送られ、ここで図5に示すように受信コイル毎の感度分布の計算503と、各受信コイルからの信号の合成処理504が行なわれる。即ち、各コイルからの信号en(kx,ky)501を用いて各コイルの感度分布画像Wn(x,y)を求める。nはコイル番号でこの実施例ではn=1,2,3又は4である。また(kx,ky)はk空間の座標を、(x,y)は、実空間上の位置を表わす。感度分布画像と信号を使って、全体画像S(x,y)505を合成する。これら感度分布計算503と合成処理504についてさらに詳述する。
Thus, the data measured by using the EPI pulse sequence and varying the density of the phase encoding step depending on the region is sent to the
図6及び図7は、各受信コイル毎に信号en(kx,ky)501からRFコイルの感度分布を計算する手順を示す図である。まず領域1のデータのみが残るように位相エンコード方向に低周波通過フィルタ(LPF)601を作用させる。フィルタは、位相エンコード方向の折返しエリアシングを除去する目的では1次元フィルタでもよいが、生体の微細構造を排除するためには1次元よりも2次元フィルタが好ましい。この場合、例えば、ガウス型、ハミング型、ハニング型などが適している。さらに高精度なフィルタリング法として、画像空間でフライングウィンドウを用いた手法も可能である。
6 and 7 are diagrams showing a procedure for calculating the sensitivity distribution of the RF coil from the signal en (kx, ky) 501 for each receiving coil. First, a low frequency pass filter (LPF) 601 is applied in the phase encoding direction so that only the data of
このようなフィルタリングの結果、図7に示すように領域2(403)のデータはすべてゼロになる。また、領域1(402)と領域2(403)の境界部分は滑らかにゼロにつながる。尚、図7中、右側の曲線はフィルタプロファイルである。 As a result of such filtering, the data in area 2 (403) is all zero as shown in FIG. Further, the boundary between the region 1 (402) and the region 2 (403) is smoothly connected to zero. In FIG. 7, the right curve is the filter profile.
次に、k空間の全体を2次元フーリエ変換(FT)602する。この結果、画像の低周波成分のみを取り出すことができる。また、エリアシングアーチファクトが無い画像が得られる。これらの画像は、RFコイルの感度分布とみなせることが知られており、本実施例でもこれを感度分布Wn(x,y)とする。尚、マルチスライス撮影の場合には、コイルの感度分布は各コイルの各スライス毎に求めておく必要がある。また3次元撮影の場合には、各コイル毎に得られる信号は3次元計測空間(kx,ky,kz)の関数となるので、低周波通過フィルタに3次元フィルタを用い、フーリエ変換を3次元に拡張することによりコイルの3次元感度分布Wn(x,y,z)を求めることができる。 Next, a two-dimensional Fourier transform (FT) 602 is performed on the entire k space. As a result, only the low frequency component of the image can be extracted. In addition, an image free from aliasing artifacts can be obtained. It is known that these images can be regarded as the sensitivity distribution of the RF coil, and this is also referred to as a sensitivity distribution Wn (x, y) in this embodiment. In the case of multi-slice imaging, the sensitivity distribution of the coil needs to be obtained for each slice of each coil. In the case of three-dimensional imaging, since the signal obtained for each coil is a function of the three-dimensional measurement space (kx, ky, kz), a three-dimensional filter is used as a low-frequency pass filter, and Fourier transformation is performed in three dimensions. The three-dimensional sensitivity distribution Wn (x, y, z) of the coil can be obtained by extending to.
次にこの感度分布Wn(x,y)を用いた合成演算504について説明する。合成演算は、画像空間(実空間)で行なう方法と、計測空間で行なう方法とがあるが、いずれも採用できる。
Next, the
図8に画像空間で行なう方法を概念的に説明する図を示す。ここでは説明を簡単にするために2つの小型コイル801、802からなるマルチプルコイルの場合を説明する。図4に示すk空間の領域2のように位相エンコード方向に粗に(間引いて)データを収集した場合、2つのコイルからの信号は、図8(b)、(c)に示すように位相エンコード方向に折返しアーチファクト804が生じた画像803になる。このような画像に対し、各小型コイルの感度分布を掛合わせることにより折返しアーチファクト804を除去することができ、同図(a)に示すように折返しによる画像の重なりがない画像が得られる。
FIG. 8 conceptually illustrates a method performed in the image space. Here, in order to simplify the description, a case of a multiple coil composed of two
計測空間で信号を合成する方法では、小型RFコイルの感度分布を適当な重みで合成して得られる合成感度分布が所望の周波数となるように重みを決めることにより、計測空間上の不足しているデータを作成する。 In the method of synthesizing signals in the measurement space, there is a shortage in the measurement space by determining the weight so that the combined sensitivity distribution obtained by combining the sensitivity distribution of the small RF coil with the appropriate weight becomes the desired frequency. To create data.
例えば、小型RFコイルの感度分布Wn(x,y)を適当な重みCnで合成して、exp(i・mΔky・y)[mは整数]の形の合成感度分布Wcomp(x,y)が得られるとする。
図9は、この方法を概念的に示す図である。図示するように、まず領域1(402)の計測データからLPF、ゼロフィルング、FT処理によって感度分布を求める。次に、間引きによって計測データが粗に配列している領域2(403)について、上記計算により既に求められたデータ(拡大図の実線に対応)から不足しているデータ(拡大図の点線に対応)を作成する。このように作成することができる新たなデータの数(点線の本数)は、小型RFコイルの数に依存し、小型RFコイルが4個の場合には新たなデータを4本まで作成できる。図示する例では、間引き率が4で、3本のデータを作成する場合を示している。 FIG. 9 is a diagram conceptually showing this method. As shown in the figure, first, sensitivity distribution is obtained from the measurement data of region 1 (402) by LPF, zero filling, and FT processing. Next, for the area 2 (403) where the measurement data is roughly arranged by thinning, the data (corresponding to the solid line in the enlarged view) already obtained by the above calculation (corresponding to the dotted line in the enlarged view) ). The number of new data (number of dotted lines) that can be created in this way depends on the number of small RF coils, and up to four new data can be created when there are four small RF coils. In the example shown in the figure, the thinning rate is 4 and three data are created.
こうして新たなデータを補充した後のデータは、k空間の全領域を間引かずに測定したデータと同じであり、これを2次元フーリエ変換することにより折返しアーチファクトのない画像を得ることができる。 The data after supplementing with new data in this way is the same as the data measured without thinning out the entire area of k-space, and an image free from aliasing artifacts can be obtained by performing a two-dimensional Fourier transform.
このように本実施例では、位相エンコードの一部を間引いて時間を短縮した撮影を行なうとともにその撮影で得られた計測データのうち密である部分のデータを用いて折り返しのないRFコイル感度分布を得るようにしたので、感度分布を求めるための計測と本計測との間の時間差をなくすことができ、それに伴う演算誤差を無くすことができる。 As described above, in this embodiment, a part of the phase encode is thinned out to shorten the time, and the measurement data obtained by the photographing is used for the dense part of the data, so that the RF coil sensitivity distribution is not turned back. Therefore, the time difference between the measurement for obtaining the sensitivity distribution and the main measurement can be eliminated, and the calculation error associated therewith can be eliminated.
特に計測データのうち、主な画像情報が含まれている低位相エンコード成分のデータを密に取得していることにより、S/Nが劣化しにくい画像を得ることができ、臨床応用において診断に極めて有効な画像を提供することができる。また予めコイルの感度分布を求めるための計測を行なう必要がないので、本計測における時間短縮効果に加え、全体としての計測時間の短縮を図ることができる。したがって、特に連続撮影の場合に好適である。 In particular, by acquiring densely the low-phase encoded component data that contains the main image information in the measurement data, it is possible to obtain an image that does not deteriorate the S / N, and is useful for diagnosis in clinical applications. An extremely effective image can be provided. Further, since it is not necessary to perform measurement for obtaining the sensitivity distribution of the coil in advance, in addition to the time shortening effect in the main measurement, the entire measurement time can be shortened. Therefore, it is particularly suitable for continuous shooting.
図10に、撮影を時系列に繰り返しリアルタイムの連続画像を得るダイナミック撮影に、本発明を適用した実施例を示す。ダイナミック撮影では、マルチプルコイルの各小型RFコイルでそれぞれ取得した信号を合成し、画像表示/転送を連続して繰り返し、時系列に連続する画像(画像番号1,2,3・・・1000)を取得する。ここで1枚目の画像取得では、図5に示すステップ502、503、504と同様に、信号取得1001、感度分布計算1002、信号合成1003の各ステップを行ない、感度分布計算1002で求めた感度分布を用いて折返しアーチファクトを除去した画像を再構成する。この画像はMRI装置のディスプレイに表示され、或いは外部の表示装置や記憶装置等に転送される(1004)。1枚目で求めた感度分布計算結果は信号処理部のメモリの特定アドレスに格納される。
FIG. 10 shows an embodiment in which the present invention is applied to dynamic shooting for obtaining continuous images in real time by repeating shooting in time series. In dynamic shooting, the signals acquired by each small RF coil of the multiple coil are combined, and image display / transfer is repeated continuously to create time-sequential images (
次いで2枚目の画像では、信号取得後、このメモリ内に格納された感度分布データを用いて折返しアーチファクトを除去する合成処理1003を行なう。以後、受信コイルの配置が変らない限り、感度分布計算は行なわず、メモリ内の感度分布データを用いて合成処理する。受信コイルの配置が変らなければ感度分布も変らないので、上述のように最初に求めた計算結果を用いることができる。これによりダイナミック撮影の撮影間隔(画像取得と次の画像取得との間隔)が短くなっても画像化のための演算を短くでき、表示のリアルタイム性(時間分解能)を向上することができる。例えば撮影シーケンスとしてEPIのような高速シーケンスを採用した場合、連続高速撮影を実効あるものにするためには演算の高速化が重要であるが、図10に示す実施例を採用することにより、リアルタイムで被検体の状況をモニタすることができる。
Next, for the second image, after the signal is acquired, a
本発明は、上述したスピンエコー型のEPIのみならずグラディエントエコー型のEPIにも、また1回の励起で必要な計測データを収集するワンショットのEPIにもマルチショット(分割型)のEPIにも適用することができる。 The present invention can be applied not only to the above-mentioned spin echo type EPI but also to a gradient echo type EPI, to a one shot EPI that collects measurement data necessary for one excitation, or to a multi-shot (split type) EPI. Can also be applied.
またダイナミック撮影の場合には、図10に示したように、感度分布計算は第1枚目の画像の再構成時にのみ行ない、それ以降の画像についてはその感度分布を用いて信号取得から画像表示/転送までの処理を連続して行なう。この場合、既に述べたように本発明のEPIシーケンスは通常のEPIシーケンスに比べ信号取得時間が短縮されており、しかも感度分布計算を本計測とは別に或いは信号取得毎に行なう必要がないので極めて高速の連続撮影が可能となる。 In the case of dynamic shooting, as shown in FIG. 10, the sensitivity distribution calculation is performed only when the first image is reconstructed, and the subsequent images are displayed from the signal acquisition to the image display using the sensitivity distribution. / Continues processing up to transfer. In this case, as described above, the EPI sequence of the present invention has a shorter signal acquisition time than a normal EPI sequence, and it is not necessary to perform sensitivity distribution calculation separately from the main measurement or for each signal acquisition. High-speed continuous shooting is possible.
以上、本発明のMRI装置の各実施の形態を説明したが、本発明はこれら実施形態に限定されることなく種々の変更が可能である。例えば、以上の説明ではパラレルMRIに用いるパルスシーケンスとしてEPIシーケンスを例示したが、FSE(ファーストスピンエコーシーケンス)、Burstシーケンス、スパイラルシーケンスなどの公知のシーケンスにも適用できる。また、3次元撮影にも拡張できる。 Although the embodiments of the MRI apparatus of the present invention have been described above, the present invention is not limited to these embodiments, and various modifications can be made. For example, the EPI sequence is exemplified as the pulse sequence used for the parallel MRI in the above description, but the present invention can also be applied to known sequences such as FSE (First Spin Echo Sequence), Burst sequence, and spiral sequence. It can also be extended to 3D photography.
図11はスパイラルシーケンスの場合のkトラジェクトリ1200を示す図である。この場合、k空間中央の円形の領域1(1202)は信号が密に配列する領域であり、その周囲の領域2(1203)は信号が粗になっている。このようなスパイラルシーケンスの計測データから感度分布を求めるためには、図12に示すようなフィルタプロファイル1301、1302の2次元フィルタを用いる。このような2次元フィルタとしては、2次元ガウスフィルタ、2次元ハニングフィルタ、2次元ハミングフィルタ等がある。領域1(1202)の計測データから求めた各小型コイル毎の感度分布を用いて、折返しアーチファクトを排除した画像を合成する点は他のシーケンスの場合と同様である。
FIG. 11 is a diagram showing the
本発明のMRI装置によれば、マルチプルコイルを用いたパラレルMRIの実施にあたり、k空間の一部の領域を間引きして計測時間を短縮した撮影を行なうとともに計測データが密である領域のデータを用いて感度分布を求め、信号を合成するので、リアルタイム性が要求される撮影で画質が劣化しない。特にk空間の低周波成分を密にデータ取得することにより、S/Nが高く診断価値の高い画像を得ることができる。 According to the MRI apparatus of the present invention, when performing parallel MRI using multiple coils, imaging is performed by shortening the measurement time by thinning out a part of the k-space, and the data of the area where the measurement data is dense is obtained. Since the sensitivity distribution is obtained and the signals are synthesized, the image quality is not deteriorated in photographing that requires real-time characteristics. In particular, it is possible to obtain an image with a high S / N and a high diagnostic value by densely acquiring low-frequency components in the k space.
また、撮影に先立ってRFコイルの感度分布を計測する必要がないので、トータルの撮影時間を延長しない。したがって、パラレルMRI技術の特徴である短時間撮影の効果が発揮できる。 In addition, since it is not necessary to measure the sensitivity distribution of the RF coil prior to photographing, the total photographing time is not extended. Therefore, the effect of short-time imaging, which is a feature of the parallel MRI technique, can be exhibited.
101・・・・・・被検体
103・・・・・・傾斜磁場コイル
104・・・・・・送信RFコイル
105・・・・・・受信RFコイル
107・・・・・・信号処理部
111・・・・・・制御部
202・・・・・・小型RFコイル
101 ・ ・ ・ ・ ・ ・ Subject
103 ・ ・ ・ ・ ・ ・ Gradient magnetic field coil
104 ・ ・ ・ ・ ・ ・ Transmitting RF coil
105 ・ ・ ・ ・ ・ ・ Receiver RF coil
107 ・ ・ ・ ・ ・ ・ Signal processing section
111 ・ ・ ・ ・ ・ ・ Control part
202 ・ ・ ・ ・ ・ ・ Small RF coil
Claims (6)
前記制御手段は、前記高周波パルスによる横磁化発生に続き、読み出し傾斜磁場および位相エンコード傾斜磁場を印加し、エコー信号を連続して計測し、その際、連続するエコー信号計測の所定の時間領域では前記位相エンコード傾斜磁場の強度を第1の値に設定し、前記所定の時間領域以外の時間領域では前記位相エンコード傾斜磁場強度を第2の値に設定することを特徴とする核磁気共鳴撮影装置。 Means for irradiating a subject including a magnetization to be detected with a radio frequency pulse to generate transverse magnetization; means for applying a phase encoding gradient magnetic field pulse to a subject to which transverse magnetization is applied; and means for applying a read gradient magnetic field pulse In the nuclear magnetic resonance imaging apparatus, comprising: means for detecting an echo signal generated during the application of the readout gradient magnetic field; control means for controlling each of the means; and image reconstruction means for reconstructing an image from the echo signal.
The control means applies a readout gradient magnetic field and a phase encoding gradient magnetic field following the generation of transverse magnetization by the high-frequency pulse, continuously measures the echo signal, and at that time, in a predetermined time region of continuous echo signal measurement A nuclear magnetic resonance imaging apparatus, wherein the intensity of the phase encoding gradient magnetic field is set to a first value, and the phase encoding gradient magnetic field intensity is set to a second value in a time region other than the predetermined time region. .
前記エコー信号を検出する手段は、互いに空間的に異なる検出感度分布を有する複数の受信コイルを備え、
前記制御手段は、前記高周波パルスによる横磁化発生に続き、読み出し傾斜磁場および位相エンコード傾斜磁場を印加し、エコー信号を連続して計測し、
前記画像再構成手段は、前記連続して計測したエコー信号の一部を用いて、各受信コイルの感度分布を計算し、計算した感度分布と前記連続計測したエコー信号とを用いて画像を再構成することを特徴とする核磁気共鳴撮影装置。 Means for irradiating a subject including a magnetization to be detected with a radio frequency pulse to generate transverse magnetization; means for applying a phase encoding gradient magnetic field pulse to a subject to which transverse magnetization is applied; and means for applying a read gradient magnetic field pulse In the nuclear magnetic resonance imaging apparatus, comprising: means for detecting an echo signal generated during the application of the readout gradient magnetic field; control means for controlling each of the means; and image reconstruction means for reconstructing an image from the echo signal.
The means for detecting the echo signal includes a plurality of receiving coils having spatially different detection sensitivity distributions,
The control means, following the generation of transverse magnetization by the high-frequency pulse, applies a readout gradient magnetic field and a phase encoding gradient magnetic field, continuously measures the echo signal,
The image reconstruction means calculates a sensitivity distribution of each receiving coil using a part of the continuously measured echo signals, and reconstructs an image using the calculated sensitivity distribution and the continuously measured echo signals. A nuclear magnetic resonance imaging apparatus characterized by comprising.
前記エコー信号を検出する手段は、互いに空間的に異なる検出感度分布を有する複数の受信コイルを備え、
前記制御手段は、前記高周波パルスによる横磁化発生に続き、読み出し傾斜磁場および位相エンコード傾斜磁場を印加し、その際、前記位相エンコード傾斜磁場の強度を異ならせて、エコー信号を連続して計測し、
前記画像再構成手段は、前記連続して計測したエコー信号の一部を用いて、各受信コイルの感度分布を計算し、計算した感度分布と前記連続計測したエコー信号とを用いて画像を再構成することを特徴とする核磁気共鳴撮影装置。 Means for irradiating a subject including a magnetization to be detected with a radio frequency pulse to generate transverse magnetization; means for applying a phase encoding gradient magnetic field pulse to a subject to which transverse magnetization is applied; and means for applying a read gradient magnetic field pulse In a nuclear magnetic resonance imaging apparatus, comprising: means for detecting an echo signal generated during application of the readout gradient magnetic field; control means for controlling the means; and image reconstruction means for reconstructing an image from the echo signal.
The means for detecting the echo signal includes a plurality of receiving coils having spatially different detection sensitivity distributions,
The control means applies a readout gradient magnetic field and a phase encoding gradient magnetic field following the generation of transverse magnetization by the high frequency pulse, and continuously measures echo signals by varying the intensity of the phase encoding gradient magnetic field. ,
The image reconstruction means calculates a sensitivity distribution of each receiving coil using a part of the continuously measured echo signals, and reconstructs an image using the calculated sensitivity distribution and the continuously measured echo signals. A nuclear magnetic resonance imaging apparatus characterized by comprising.
前記制御手段は、連続するエコー信号計測の所定の時間領域では前記位相エンコード傾斜磁場の強度を第1の値に設定し、前記所定の時間領域以外の時間領域では前記位相エンコード傾斜磁場強度を第2の値に設定することを特徴とする核磁気共鳴撮影装置。 The nuclear magnetic resonance imaging apparatus according to claim 3,
The control means sets the intensity of the phase encoding gradient magnetic field to a first value in a predetermined time region of continuous echo signal measurement, and sets the phase encoding gradient magnetic field strength to a first value in a time region other than the predetermined time region. A nuclear magnetic resonance imaging apparatus, wherein the value is set to 2.
前記画像再構成手段は、1回または2回以上の横磁化発生に続いて連続計測したエコー信号を用いて1枚の画像を再構成し、その際、前記位相エンコード傾斜磁場の強度を第1又は第2の値に設定して計測したエコー信号を用いて、各受信コイルの感度分布を計算し、計算した感度分布と前記連続計測したエコー信号とを用いて画像を再構成することを特徴とした核磁気共鳴撮影装置。 The nuclear magnetic resonance imaging apparatus according to claim 4,
The image reconstruction means reconstructs one image using echo signals continuously measured after one or more transverse magnetization generations, and at this time, the intensity of the phase encoding gradient magnetic field is set to a first value. Alternatively, the sensitivity distribution of each receiving coil is calculated using the echo signal measured with the second value set, and the image is reconstructed using the calculated sensitivity distribution and the continuously measured echo signal. Nuclear magnetic resonance imaging equipment.
前記制御手段によるエコー信号の連続計測は、エコプラナー法またはファーストスピンエコー法のパルスシーケンスに基づき実行されることを特徴とする核磁気共鳴撮影装置。 The nuclear magnetic resonance imaging apparatus according to any one of claims 1 to 5,
The nuclear magnetic resonance imaging apparatus, wherein the echo signal is continuously measured by the control means based on a pulse sequence of an ecoplanar method or a first spin echo method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2006292684A JP3983792B2 (en) | 2006-10-27 | 2006-10-27 | Nuclear magnetic resonance imaging system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2006292684A JP3983792B2 (en) | 2006-10-27 | 2006-10-27 | Nuclear magnetic resonance imaging system |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP34932399A Division JP3952247B2 (en) | 1999-12-08 | 1999-12-08 | Nuclear magnetic resonance imaging system |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2007151364A Division JP4675936B2 (en) | 2007-06-07 | 2007-06-07 | Nuclear magnetic resonance imaging system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2007050278A true JP2007050278A (en) | 2007-03-01 |
| JP3983792B2 JP3983792B2 (en) | 2007-09-26 |
Family
ID=37915090
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2006292684A Expired - Fee Related JP3983792B2 (en) | 2006-10-27 | 2006-10-27 | Nuclear magnetic resonance imaging system |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3983792B2 (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7877129B2 (en) | 2001-09-13 | 2011-01-25 | Hitachi Medical Corporation | Magnetic resonance imaging apparatus and RF reception coil apparatus |
| CN113554566A (en) * | 2021-07-29 | 2021-10-26 | 杭州电子科技大学 | A system and method for removing moiré based on learnable frequency domain priors |
| CN116559740A (en) * | 2023-03-16 | 2023-08-08 | 中国科学院精密测量科学与技术创新研究院 | NMR method and system for accurately measuring peak field intensity distribution of pulse strong magnetic field |
-
2006
- 2006-10-27 JP JP2006292684A patent/JP3983792B2/en not_active Expired - Fee Related
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7877129B2 (en) | 2001-09-13 | 2011-01-25 | Hitachi Medical Corporation | Magnetic resonance imaging apparatus and RF reception coil apparatus |
| CN113554566A (en) * | 2021-07-29 | 2021-10-26 | 杭州电子科技大学 | A system and method for removing moiré based on learnable frequency domain priors |
| CN113554566B (en) * | 2021-07-29 | 2024-04-02 | 杭州电子科技大学 | Moire removing system and method based on learnable frequency domain priori |
| CN116559740A (en) * | 2023-03-16 | 2023-08-08 | 中国科学院精密测量科学与技术创新研究院 | NMR method and system for accurately measuring peak field intensity distribution of pulse strong magnetic field |
| CN116559740B (en) * | 2023-03-16 | 2024-01-12 | 中国科学院精密测量科学与技术创新研究院 | NMR method and system for accurately measuring peak field intensity distribution of pulse strong magnetic field |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3983792B2 (en) | 2007-09-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3952247B2 (en) | Nuclear magnetic resonance imaging system | |
| US6777934B2 (en) | Magnetic resonance imaging method and apparatus | |
| JP4152381B2 (en) | Magnetic resonance imaging system | |
| CN102597795B (en) | MR imaging using navigators | |
| JP2003079595A (en) | Magnetic resonance imaging device and rf receiving coil therefor | |
| US7319324B2 (en) | MRI method and apparatus using PPA image reconstruction | |
| WO2002065907A1 (en) | Magnetic resonance imaging apparatus and method | |
| JP4072879B2 (en) | Nuclear magnetic resonance imaging system | |
| JP4047553B2 (en) | Magnetic resonance imaging system | |
| JP5536358B2 (en) | Magnetic resonance imaging apparatus and sensitivity correction method | |
| JP5722212B2 (en) | Magnetic resonance imaging apparatus and method | |
| US9535148B2 (en) | Dynamic contrast enhanced magnetic resonance imaging with high spatial-temporal resolution | |
| JP4202855B2 (en) | Magnetic resonance imaging system | |
| JP4679158B2 (en) | Magnetic resonance imaging system | |
| JP4675936B2 (en) | Nuclear magnetic resonance imaging system | |
| JP4698231B2 (en) | Magnetic resonance diagnostic equipment | |
| JP3983792B2 (en) | Nuclear magnetic resonance imaging system | |
| JP4906952B2 (en) | Nuclear magnetic resonance imaging system | |
| JP3952310B2 (en) | Nuclear magnetic resonance imaging system | |
| JP4980662B2 (en) | Magnetic resonance imaging system | |
| JP2008055023A (en) | Magnetic resonance imaging apparatus | |
| JP4817381B2 (en) | Magnetic resonance imaging system | |
| JP4326910B2 (en) | Magnetic resonance imaging system | |
| Xu et al. | Integrating parallel imaging with generalized series for accelerated dynamic imaging |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20061204 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20061204 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070410 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070607 |
|
| 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: 20070703 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20070704 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100713 Year of fee payment: 3 |
|
| 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: 20110713 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120713 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120713 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130713 Year of fee payment: 6 |
|
| 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 |