JP2011177245A - Magnetic resonance imaging apparatus - Google Patents
Magnetic resonance imaging apparatus Download PDFInfo
- Publication number
- JP2011177245A JP2011177245A JP2010042679A JP2010042679A JP2011177245A JP 2011177245 A JP2011177245 A JP 2011177245A JP 2010042679 A JP2010042679 A JP 2010042679A JP 2010042679 A JP2010042679 A JP 2010042679A JP 2011177245 A JP2011177245 A JP 2011177245A
- Authority
- JP
- Japan
- Prior art keywords
- displacement amount
- sequence
- magnetic resonance
- imaging apparatus
- diaphragm
- 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
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Description
本発明は、体動を補正するための磁気共鳴イメージング装置に関する。 The present invention relates to a magnetic resonance imaging apparatus for correcting body movement.
患者の呼吸による体動アーチファクトを低減するために、患者に息止めを指示し、患者が息止めをしている間に撮影を行う場合がある。しかし、この方法では、息止めが困難な患者に対しては十分に体動アーチファクトを低減することができない。そこで、呼吸同期法によって撮影を行う方法が知られている（特許文献１参照）。 In order to reduce body movement artifacts due to patient breathing, the patient may be instructed to hold his / her breath, and imaging may be performed while the patient is holding his / her breath. However, this method cannot sufficiently reduce body movement artifacts for patients who are difficult to hold their breath. Therefore, a method of performing imaging by the respiratory synchronization method is known (see Patent Document 1).
特許文献１の方法では、体動アーチファクトを低減するために、呼吸が安定しているときにデータを収集している。したがって、呼吸が安定している間にデータを収集した後、再び呼吸が安定するまではデータを収集することができず、撮影時間が長くなるという問題があり、この問題を解決することが望まれている。 In the method of Patent Document 1, data is collected when breathing is stable in order to reduce body movement artifacts. Therefore, after collecting data while breathing is stable, data cannot be collected until breathing stabilizes again, and there is a problem that the imaging time becomes long, and it is hoped to solve this problem. It is rare.
位相エンコード方向に勾配磁場を印加し、被検体の体動に伴って変位する所定の部位を撮影するための本スキャンを実行する磁気共鳴イメージング装置であって、
前記所定の部位にタグが付されたタギング画像のデータを収集するためのトレーニングスキャンを実行するスキャン手段と、
前記タギング画像のデータに基づいて、前記所定の部位の変位量を算出する第１の変位量算出手段と、
前記第１の変位量算出手段により算出された前記所定の部位の変位量に基づいて、前記本スキャンのｎ回目の位相エンコードにおける前記所定の部位の変位量を算出する第２の変位量算出手段と、
前記第２の変位量算出手段が算出した前記所定の部位の変位量と、前記本スキャンにより得られるｋ空間のデータとに基づいて、前記所定の部位の画像データを算出する画像データ算出手段と、
を有する。
A magnetic resonance imaging apparatus that performs a main scan for applying a gradient magnetic field in a phase encoding direction and imaging a predetermined portion that is displaced in accordance with a body motion of a subject,
Scanning means for executing a training scan for collecting tagging image data with a tag attached to the predetermined site;
First displacement amount calculating means for calculating a displacement amount of the predetermined part based on the tagging image data;
Second displacement amount calculation means for calculating the displacement amount of the predetermined portion in the n-th phase encoding of the main scan based on the displacement amount of the predetermined portion calculated by the first displacement amount calculation means. When,
Image data calculating means for calculating image data of the predetermined part based on the displacement amount of the predetermined part calculated by the second displacement amount calculating means and k-space data obtained by the main scan; ,
Have
タギングデータのタグの歪を解析することにより、所定の部位がどれだけ変位しているかが分かるので、体動補正をすることができる。 By analyzing the tag distortion of the tagging data, it is possible to know how much the predetermined part has been displaced, so that the body movement can be corrected.
以下、発明の実施するための形態について説明するが、発明を実施するための形態は、以下の形態に限定されることはない。 Hereinafter, although the form for inventing is demonstrated, the form for inventing is not limited to the following forms.
（１）第１の実施形態
図１は、本発明の第１の実施形態の磁気共鳴イメージング装置を示す概略図である。
磁気共鳴イメージング（ＭＲＩ（Magnetic Resonance Imaging））装置１は、磁場発生装置２と、テーブル３と、クレードル４と、受信コイル５などを有している。
(1) First Embodiment FIG. 1 is a schematic view showing a magnetic resonance imaging apparatus according to a first embodiment of the present invention.
A magnetic resonance imaging (MRI) apparatus 1 includes a magnetic field generator 2, a table 3, a cradle 4, a receiving coil 5, and the like.
磁場発生装置２は、被検体１３が収容されるボア２１と、超伝導コイル２２と、勾配コイル２３と、送信コイル２４とを有している。超伝導コイル２２は静磁場Ｂ0を印加し、勾配コイル２３は、周波数エンコード方向、位相エンコード方向、およびスライス選択方向に勾配磁場を印加する。また、送信コイル２４はＲＦパルスを送信する。尚、超伝導コイル２２の代わりに、永久磁石を用いてもよい。 The magnetic field generator 2 includes a bore 21 in which the subject 13 is accommodated, a superconducting coil 22, a gradient coil 23, and a transmission coil 24. The superconducting coil 22 applies a static magnetic field B0, and the gradient coil 23 applies a gradient magnetic field in the frequency encoding direction, the phase encoding direction, and the slice selection direction. The transmission coil 24 transmits an RF pulse. In place of the superconducting coil 22, a permanent magnet may be used.
クレードル４は、テーブル３からボア２１に移動できるように構成されている。クレードル４によって、被検体１３はボア２１に搬送される。 The cradle 4 is configured to be movable from the table 3 to the bore 21. The subject 13 is transported to the bore 21 by the cradle 4.
受信コイル５は、被検体１３の胸部から腹部にかけて取り付けられている。受信コイル５は、撮影部位からの磁気共鳴信号を受信する。 The receiving coil 5 is attached from the chest to the abdomen of the subject 13. The receiving coil 5 receives a magnetic resonance signal from the imaging region.
ＭＲＩ装置１は、更に、シーケンサ６、送信器７、勾配磁場電源８、受信器９、中央処理装置１０、入力装置１１、および表示装置１２を有している。 The MRI apparatus 1 further includes a sequencer 6, a transmitter 7, a gradient magnetic field power supply 8, a receiver 9, a central processing unit 10, an input device 11, and a display device 12.
シーケンサ６は、中央処理装置１０の制御を受けて、後述するトレーニングスキャンおよび本スキャン（図２参照）を実行するための情報を送信器７および勾配磁場電源８に送る。具体的には、シーケンサ６は、中央処理装置１０の制御を受けて、ＲＦパルスの情報（中心周波数、バンド幅など）を送信器７に送り、勾配磁場の情報（勾配磁場の強度など）を勾配磁場電源８に送る。 Under the control of the central processing unit 10, the sequencer 6 sends information for executing a training scan and a main scan (see FIG. 2) described later to the transmitter 7 and the gradient magnetic field power supply 8. Specifically, under the control of the central processing unit 10, the sequencer 6 sends RF pulse information (center frequency, bandwidth, etc.) to the transmitter 7, and gradient magnetic field information (gradient magnetic field strength, etc.). Send to gradient magnetic field power supply 8.
送信器７は、シーケンサ６から送られた情報に基づいて、ＲＦコイル２４を駆動する駆動信号を出力する。 The transmitter 7 outputs a drive signal for driving the RF coil 24 based on the information sent from the sequencer 6.
勾配磁場電源８は、シーケンサ６から送られた情報に基づいて、勾配コイル２３を駆動する駆動信号を出力する。 The gradient magnetic field power supply 8 outputs a drive signal for driving the gradient coil 23 based on the information sent from the sequencer 6.
受信器９は、受信コイル５で受信された磁気共鳴信号を信号処理し、中央処理装置１０に伝送する。 The receiver 9 processes the magnetic resonance signal received by the receiving coil 5 and transmits it to the central processing unit 10.
中央処理装置１０は、シーケンサ６および表示装置１２に必要な情報を伝送したり、受信器９から受け取った信号に基づいて画像を再構成するなど、ＭＲＩ装置１の各種の動作を実現するように、ＭＲＩ装置１の各部の動作を総括する。中央処理装置１０は、例えばコンピュータ（computer）によって構成される。中央処理装置１０は、第１の変位量算出手段１０１、第２の変位量算出手段１０２、および画像データ算出手段１０３を有している。 The central processing unit 10 implements various operations of the MRI apparatus 1 such as transmitting necessary information to the sequencer 6 and the display unit 12 and reconstructing an image based on a signal received from the receiver 9. The operation of each part of the MRI apparatus 1 will be summarized. The central processing unit 10 is configured by, for example, a computer. The central processing unit 10 includes a first displacement amount calculation unit 101, a second displacement amount calculation unit 102, and an image data calculation unit 103.
第１の変位量算出手段１０１は、トレーニングスキャン（図２参照）により得られたタギング画像のデータに基づいて、所定の部位の変位量を算出する。 The first displacement amount calculation means 101 calculates a displacement amount of a predetermined part based on tagging image data obtained by the training scan (see FIG. 2).
第２の変位量算出手段１０２は、第１の変位量算出手段１０１により算出された所定の部位の変位量に基づいて、本スキャン（図２参照）のｎ回目の位相エンコードにおける所定の部位の変位量を算出する。 Based on the displacement amount of the predetermined part calculated by the first displacement amount calculation means 101, the second displacement amount calculation means 102 determines the predetermined part in the n-th phase encoding of the main scan (see FIG. 2). The displacement amount is calculated.
画像データ算出手段１０３は、第２の変位量算出手段１０２が算出した所定の部位の変位量と、本スキャンにより得られるｋ空間のデータとに基づいて、所定の部位の画像データを算出する。また、画像データ算出手段１０３は、タギング画像のデータも算出する。 The image data calculation unit 103 calculates image data of a predetermined part based on the displacement amount of the predetermined part calculated by the second displacement amount calculation unit 102 and k-space data obtained by the main scan. The image data calculation unit 103 also calculates tagging image data.
中央処理装置１０は、第１の変位量算出手段１０１、第２の変位量算出手段１０２、および画像データ算出手段１０３の一例であり、所定のプログラムを実行することにより、これらの手段として機能する。 The central processing unit 10 is an example of the first displacement amount calculation unit 101, the second displacement amount calculation unit 102, and the image data calculation unit 103, and functions as these units by executing a predetermined program. .
入力装置１１は、オペレータ１４の操作に応答して種々の命令を中央処理装置１０に入力する。表示装置１２は種々の情報を表示する。 The input device 11 inputs various commands to the central processing unit 10 in response to the operation of the operator 14. The display device 12 displays various information.
尚、磁場発生装置２と、シーケンサ６と、送信器７と、勾配磁場電源８とを合わせたものが、課題を解決するための手段に記載されたスキャン手段に相当する。
磁気共鳴イメージング装置１は、上記のように構成されている。
次に、被検体１３を撮影するときに実行されるスキャンについて説明する。
A combination of the magnetic field generator 2, the sequencer 6, the transmitter 7, and the gradient magnetic field power supply 8 corresponds to the scanning means described in the means for solving the problem.
The magnetic resonance imaging apparatus 1 is configured as described above.
Next, a scan executed when imaging the subject 13 will be described.
図２は、被検体１３を撮影するときに実行されるスキャンの説明図、図３は、撮影部位を概略的に示す図である。 FIG. 2 is an explanatory diagram of a scan executed when the subject 13 is imaged, and FIG. 3 is a diagram schematically showing an imaging region.
第１の実施形態では、トレーニングスキャンおよび本スキャンが実行される。
トレーニングスキャンは、変位量Δｘ（ｕ，ｘ，ｙ）およびΔｙ（ｕ，ｘ，ｙ）を決定するために実行されるスキャンである。変位量Δｘ（ｕ，ｘ，ｙ）、Δｙ（ｕ，ｘ，ｙ）については、後述する。
In the first embodiment, a training scan and a main scan are executed.
The training scan is a scan executed to determine the displacement amounts Δx (u, x, y) and Δy (u, x, y). The displacement amounts Δx (u, x, y) and Δy (u, x, y) will be described later.
トレーニングスキャンでは、タギングシーケンスＴＳａ（ａ＝１〜Ａ）と複数のデータ収集シーケンスＤＳとの組合せSetが、複数回実行される。タギングシーケンスＴＳａは、イメージング領域Ｒimにタグを付するためのシーケンスである。データ収集シーケンスＤＳは、タグが付されたイメージング領域Ｒimからｋ空間のデータを収集するためのシーケンスである。 In the training scan, the combination Set of the tagging sequence TS a (a = 1 to A) and the plurality of data collection sequences DS is executed a plurality of times. The tagging sequence TSa is a sequence for attaching a tag to the imaging region Rim. The data collection sequence DS is a sequence for collecting k-space data from the tagged imaging region Rim.
また、タギングシーケンスＴＳａが実行される前と、データ収集シーケンスＤＳが実行される前には、一回又は複数回のナビゲータシーケンスＮＴが実行される。ナビゲータシーケンスＮＴは、肝臓１３ａ、横隔膜１３ｂ、および肺１３ｃに跨る柱状のナビゲータ領域Ｒna（図３参照）をスライス選択し、ナビゲータ領域Ｒnaから、被検体の体軸方向に関する横隔膜１３ｂの位置ｕを検出するのに必要なデータを収集するためのシーケンスである。 Further, before and tagging sequence TS a is executed before the data acquisition sequence DS is performed one or more times of navigator sequence NT is executed. The navigator sequence NT selects a columnar navigator region Rna (see FIG. 3) straddling the liver 13a, the diaphragm 13b, and the lung 13c, and detects the position u of the diaphragm 13b in the body axis direction of the subject from the navigator region Rna. This is a sequence for collecting data necessary to perform the process.
図４は、トレーニングスキャンの説明図である。
図４（ａ）は、トレーニングスキャン時における横隔膜の変位を表すグラフである。
FIG. 4 is an explanatory diagram of the training scan.
FIG. 4A is a graph showing the displacement of the diaphragm during a training scan.
グラフの横軸は時間を表し、グラフの縦軸は横隔膜位置ｕを示している。グラフの黒丸は、ナビゲータシーケンスＮＴにより検出された横隔膜位置ｕを表している。 The horizontal axis of the graph represents time, and the vertical axis of the graph represents the diaphragm position u. A black circle in the graph represents the diaphragm position u detected by the navigator sequence NT.
トレーニングスキャンでは、イメージング領域Ｒimのタギング画像を作成するためのデータを収集する。図４（ｂ）には、横隔膜が位置ｕ’αに変位したときのイメージング領域Ｒimにおけるタギング画像ＴＧαの概略図が示されており、図４（ｃ）には、横隔膜が位置ｕ’βに変位したときのイメージング領域Ｒimにおけるタギング画像ＴＧβの概略図が示されている。以下に、タギング画像ＴＧαのデータと、タギング画像ＴＧβのデータとを収集する方法について説明する。 In the training scan, data for creating a tagging image in the imaging region Rim is collected. FIG. 4B shows a schematic diagram of the tagging image TG α in the imaging region Rim when the diaphragm is displaced to the position u ′ α , and FIG. 4C shows the position of the diaphragm at the position u ′ β. A schematic diagram of the tagging image TG β in the imaging region Rim when displaced to λ is shown. Hereinafter, a method for collecting the data of the tagging image TG α and the data of the tagging image TG β will be described.
トレーニングスキャンでは、ナビゲータシーケンスＮＴを実行し、横隔膜位置ｕを検出する。第１の実施形態では、検出された横隔膜位置の時間変化に基づいて、横隔膜位置ｕが、減少から増大に変化したか否かを判断する。図４（ａ）では、時刻ｔ１において、横隔膜位置ｕが、減少から増大に変化している。時刻ｔ１以降も、ナビゲータシーケンスＮＴによって横隔膜位置ｕを検出し、横隔膜が基準位置ｕrefに変位したら、タギングシーケンスＴＳ１を実行する。この場合、基準位置ｕrefに対して許容範囲ＡＷrefを予め決めておき、検出された横隔膜位置ｕが許容範囲ＡＷrefに含まれている場合は、横隔膜は基準位置ｕrefに変位したと判断して、タギングシーケンスＴＳ１を実行する。タギングシーケンスＴＳ１を実行することによりイメージング領域Ｒimにタグが付される。タギングシーケンスＴＳ１を実行した後、再びナビゲータシーケンスＮＴを実行し、横隔膜位置ｕを検出する。横隔膜が、位置ｕ’αに変位したら、データ収集シーケンスＤＳを実行する。この場合、位置ｕ’αに対して許容範囲ＡＷαを予め決めておき、検出された横隔膜位置ｕが許容範囲ＡＷαに含まれている場合は、横隔膜は位置ｕ’αに変位したと判断して、データ収集シーケンスＤＳを実行する。データ収集シーケンスＤＳを実行することにより、タギングシーケンスＴＳ１でタグが付されたイメージング領域Ｒimから、横隔膜がｕ’αに変位したときのｋ空間のデータＤα１を収集することができる。 In the training scan, the navigator sequence NT is executed to detect the diaphragm position u. In the first embodiment, it is determined whether or not the diaphragm position u has changed from a decrease to an increase based on the detected temporal change in the diaphragm position. In FIG. 4A, the diaphragm position u changes from decrease to increase at time t1. Time t1 and beyond, to detect the diaphragm position u by navigator sequence NT, diaphragmatic When displaced to the reference position u ref, executes the tagging sequence TS 1. In this case, an allowable range AW ref is determined in advance with respect to the reference position u ref , and when the detected diaphragm position u is included in the allowable range AW ref , it is determined that the diaphragm is displaced to the reference position u ref. and, to perform the tagging sequence TS 1. Tag is attached to the imaging region Rim by performing the tagging sequence TS 1. After executing the tagging sequence TS 1, perform the navigator sequence NT again detects the diaphragm position u. When the diaphragm is displaced to the position u ′ α , the data acquisition sequence DS is executed. In this case, it determines that the position u 'determined in advance an allowable range AW alpha relative to alpha, if the detected diaphragm position u is included in the allowable range AW alpha is the diaphragm position u' displaced alpha Then, the data collection sequence DS is executed. By executing the data acquisition sequence DS, it can be from the imaging region Rim which tagged with tagging sequence TS 1, collecting data D [alpha] 1 of the k space when the diaphragm is displaced in the u 'alpha.
ｋ空間のデータＤα１を収集した後、再びナビゲータシーケンスＮＴを実行し、横隔膜位置ｕを検出する。横隔膜が、位置ｕ’βに変位したら、データ収集シーケンスＤＳを実行する。この場合、位置ｕ’βに対して許容範囲ＡＷβを予め決めておき、検出された横隔膜位置ｕが許容範囲ＡＷβに含まれている場合は、横隔膜は位置ｕ’βに変位したと判断して、データ収集シーケンスＤＳを実行する。データ収集シーケンスＤＳを実行することにより、タギングシーケンスＴＳ１でタグが付されたイメージング領域Ｒimから、横隔膜がｕ’βに変位したときのｋ空間のデータＤβ１を収集することができる。 After collecting the k-space data Dα1 , the navigator sequence NT is executed again to detect the diaphragm position u. When the diaphragm is displaced to the position u ′ β , the data acquisition sequence DS is executed. In this case, an allowable range AW β is determined in advance for the position u ′ β , and if the detected diaphragm position u is included in the allowable range AW β , it is determined that the diaphragm is displaced to the position u ′ β. Then, the data collection sequence DS is executed. By executing the data acquisition sequence DS, from the imaging region Rim which tagged with tagging sequence TS 1, the diaphragm can be collected the data D .beta.1 of k space when displaced u 'beta.
ｋ空間のデータＤβ１を収集した後、再びナビゲータシーケンスＮＴを実行し、横隔膜位置ｕを検出する。そして、検出された横隔膜位置の時間変化に基づいて、横隔膜位置ｕが、減少から増大に変化したか否かを判断する。図４（ａ）では、時刻ｔ２において、横隔膜位置ｕが、減少から増大に変化している。時刻ｔ２以降も、ナビゲータシーケンスＮＴによって横隔膜位置ｕを検出し、横隔膜が基準位置ｕrefに変位したら、タギングシーケンスＴＳ２を実行する。タギングシーケンスＴＳ２を実行することによりイメージング領域Ｒimに再びタグが付される。タギングシーケンスＴＳ２を実行した後、ナビゲータシーケンスＮＴによって再び横隔膜位置ｕを検出する。横隔膜が、位置ｕ’αおよびｕ’βおよび変位したら、データ収集シーケンスＤＳを実行し、ｋ空間のデータＤα２およびＤβ２を収集する。 After collecting the k-space data Dβ1 , the navigator sequence NT is executed again to detect the diaphragm position u. Then, based on the detected temporal change in the diaphragm position, it is determined whether or not the diaphragm position u has changed from a decrease to an increase. In FIG. 4A, the diaphragm position u changes from decrease to increase at time t2. Time t2 beyond detects the diaphragm position u by navigator sequence NT, diaphragmatic When displaced to the reference position u ref, executes the tagging sequence TS 2. Again Tag imaging region Rim is assigned by executing the tagging sequence TS 2. After executing the tagging sequence TS 2, again detecting a diaphragm position u by navigator sequence NT. Once the diaphragm is in position u ′ α and u ′ β and displaced, a data acquisition sequence DS is performed to collect k-space data D α2 and D β2 .
以下、同様の手順で、タギング画像ＴＧαおよびＴＧβを作成するのに必要なｋ空間のデータが収集されるまでスキャンを実行する。したがって、タギング画像ＴＧαおよびＴＧβを作成することができる。 Thereafter, the scan is executed in the same procedure until k-space data necessary to create the tagging images TG α and TG β is acquired. Therefore, tagging images TG α and TG β can be created.
また、横隔膜位置ｕは時間とともに変位するので、イメージング領域Ｒimに付されたタグも、時間とともに歪んでいく。したがって、タギング画像ＴＧαのタグと、タギング画像ＴＧβのタグには、イメージング領域Ｒimにタグを付してからデータが収集されるまでの時間に応じた歪が生じる。尚、図４（ｂ）および（ｃ）では、説明の便宜上、タギング画像ＴＧαおよびＴＧβには、歪のないタグが示されているが、実際には、タグには、イメージング領域Ｒimにタグを付してからデータが収集されるまでの時間に応じた歪が生じている。 Since the diaphragm position u is displaced with time, the tag attached to the imaging region Rim is also distorted with time. Therefore, the tag of the tagging image TG α and the tag of the tagging image TG β are distorted according to the time from when the tag is added to the imaging region Rim until data is collected. In FIGS. 4B and 4C, the tagging images TG α and TG β show tags without distortion for convenience of explanation, but in actuality, the tags include the imaging region Rim. There is a distortion corresponding to the time from when the tag is attached until the data is collected.
次に、図４（ｂ）を参照しながら、イメージング領域Ｒim内の組織Ｔの変位について説明する。 Next, the displacement of the tissue T in the imaging region Rim will be described with reference to FIG.
横隔膜の基準位置ｕrefにおいて、イメージング領域Ｒim内の組織Ｔが、ピクセル位置Ｐ０＝（ｘ０，ｙ０）に位置しているとする。しかし、横隔膜位置ｕが、位置ｕ’αに変位すると、組織Ｔは、ピクセル位置Ｐ０＝（ｘ０，ｙ０）から、ピクセル位置Ｐα＝（ｘα，ｙα）に移動する。つまり、組織Ｔは、イメージング領域Ｒimにタグが付されてからデータが収集されるまでの間に、ｘ方向に変位量Δｘだけ変位し、ｙ方向に変位量Δｙだけ変位する（尚、ｚ方向の変位は無視する）。また、変位量ΔｘおよびΔｙは、横隔膜位置ｕ、ピクセル位置（ｘ，ｙ）に依存する値であるので、ΔｘおよびΔｙは、以下のように、ｕ、ｘ、およびｙの関数として表すことができる。
Δｘ（ｕ，ｘ，ｙ）、Δｙ（ｕ，ｘ，ｙ）
It is assumed that the tissue T in the imaging region Rim is located at the pixel position P0 = (x 0 , y 0 ) at the reference position u ref of the diaphragm. However, when the diaphragm position u is displaced to the position u ′ α , the tissue T moves from the pixel position P0 = (x 0 , y 0 ) to the pixel position Pα = (x α , y α ). That is, the tissue T is displaced by the displacement amount Δx in the x direction and is displaced by the displacement amount Δy in the y direction between the time when the imaging region Rim is tagged and the data is collected (the z direction). Ignore the displacement. Further, since the displacement amounts Δx and Δy are values depending on the diaphragm position u and the pixel position (x, y), Δx and Δy can be expressed as functions of u, x, and y as follows. it can.
Δx (u, x, y), Δy (u, x, y)
変位量Δｘ（ｕ，ｘ，ｙ）およびΔｙ（ｕ，ｘ，ｙ）の値は、タギング画像に付されたタグの歪を解析することによって算出することができる。例えば、横隔膜位置ｕ＝ｕ’αにおける変位量Δｘ（ｕ’α，ｘ，ｙ）およびΔｙ（ｕ’α，ｘ，ｙ）を算出するには、タギング画像ＴＧαに付されたタグの歪を解析すればよい。タギング画像ＴＧαに付されたタグの歪を解析することによって、変位量Δｘ（ｕ’α，ｘ，ｙ）およびΔｙ（ｕ’α，ｘ，ｙ）を算出することができる。 The values of the displacement amounts Δx (u, x, y) and Δy (u, x, y) can be calculated by analyzing the distortion of the tag attached to the tagging image. For example, in order to calculate the displacement amounts Δx (u ′ α , x, y) and Δy (u ′ α , x, y) at the diaphragm position u = u ′ α , the distortion of the tag attached to the tagging image TG α Can be analyzed. By analyzing the distortion of the tag attached to tagged image TG alpha, it can be calculated displacement amount Δx (u 'α, x, y) and Δy (u' α, x, y) a.
上記の説明では、横隔膜位置ｕ’αにおける変位量Δｘ（ｕ’α，ｘ，ｙ）およびΔｙ（ｕ’α，ｘ，ｙ）を算出する例について説明したが、横隔膜位置ｕ’βにおける変位量Δｘ（ｕ’β，ｘ，ｙ）およびΔｙ（ｕ’β，ｘ，ｙ）を算出する場合には、タギング画像ＴＧβに付されたタグの歪を解析すればよい。タギング画像ＴＧβに付されたタグの歪を解析することにより、横隔膜位置ｕ’βにおける変位量Δｘ（ｕ’β，ｘ，ｙ）およびΔｙ（ｕ’β，ｘ，ｙ）を算出することができる。 In the above description, the diaphragm position u 'displacement at α Δx (u' α, x , y) and Δy (u 'α, x, y) has been described an example of calculating a diaphragm position u' displacement at β In calculating the quantities Δx (u ′ β , x, y) and Δy (u ′ β , x, y), the distortion of the tag attached to the tagging image TG β may be analyzed. By calculating the distortion of the tag attached to the tagging image TG β , the displacement amounts Δx (u ′ β , x, y) and Δy (u ′ β , x, y) at the diaphragm position u ′ β are calculated. Can do.
図４では、横隔膜位置ｕ’αにおける変位量Δｘ（ｕ’α，ｘ，ｙ）およびΔｙ（ｕ’α，ｘ，ｙ）と、横隔膜位置ｕ’βにおける変位量Δｘ（ｕ’β，ｘ，ｙ）およびΔｙ（ｕ’β，ｘ，ｙ）とを算出する例について説明している。しかし、タギングシーケンスを実行してからデータを収集するまでの時間を更に細かく規定することにより、横隔膜位置ｕ’αおよびｕ’β以外の横隔膜位置における変位量も算出することができる。第１の実施形態では、横隔膜位置ｕ’ｋ（ｋ＝１〜ｚ）における変位量を算出している（図５参照）。 In Figure 4, the diaphragm position u 'displacement at α Δx (u' α, x , y) and Δy (u 'α, x, y) and, diaphragm position u' displacement at beta [Delta] x (u 'beta, x , Y) and Δy (u ′ β , x, y) are described. However, the displacement amount at the diaphragm positions other than the diaphragm positions u ′ α and u ′ β can also be calculated by further specifying the time from the execution of the tagging sequence to the collection of data. In the first embodiment, the amount of displacement at the diaphragm position u ′ k (k = 1 to z) is calculated (see FIG. 5).
図５は、横隔膜位置ｕ’１〜ｕ’ｚにおける変位量を算出するときの説明図である。
例えば、横隔膜位置ｕ＝ｕ’ｚにおける変位量Δｘ（ｕ’ｚ，ｘ，ｙ）およびΔｙ（ｕ’ｚ，ｘ，ｙ）を算出するには、横隔膜位置ｕ’ｚおけるタギング画像の歪を解析すればよい。横隔膜位置ｕ’ｚおけるタギング画像を作成するためのｋ空間のデータは、横隔膜が位置ｕ’ｚに変位したときにデータ収集シーケンスＤＳを実行することにより、収集することができる。図６に、横隔膜位置ｕ’１〜ｕ’ｚにおける変位量をまとめて示す。
FIG. 5 is an explanatory diagram when the displacement amount at the diaphragm positions u ′ 1 to u ′ z is calculated.
For example, the diaphragm position u = u 'displacement in z [Delta] x (u' z, x, y) and [Delta] y (u 'z, x, y) to calculate the the diaphragm position u' strain of z definitive tagging images Just analyze it. Diaphragm position u 'data k-space to create z definitive tagging images, the diaphragm is positioned u' by executing the data acquisition sequence DS when displaced z, can be collected. 6 collectively shows the displacement at the diaphragm position u '1 ~u' z.
トレーニングスキャンの後に、本スキャンが実行される（図２参照）。本スキャンでは、イメージングシーケンスＩＳｂ（ｂ＝１〜Ｂ）が実行される。 The main scan is executed after the training scan (see FIG. 2). In the main scan, the imaging sequence IS b (b = 1 to B) is executed.
イメージングシーケンスＩＳｂは、イメージング領域Ｒim（図３参照）のｋ空間の各ラインのデータを収集するためのシーケンスである。第１の実施形態では、イメージングシーケンスＩＳｂは、グラディエントエコー系のシーケンスである。ただし、イメージングシーケンスＩＳｂは、スピンエコー系など、別のシーケンスを用いてもよい。位相エンコード数Ｎは、例えば、Ｎ＝１２８である。 Imaging sequence IS b is a sequence for collecting data for each line of k-space imaging region Rim (see FIG. 3). In the first embodiment, the imaging sequence IS b is a gradient echo sequence. However, the imaging sequence IS b, such as spin echo system, it may be used another sequence. The number N of phase encodes is, for example, N = 128.
また、イメージングシーケンスＩＳｂが実行される前に、一回又は複数回のナビゲータシーケンスＮＶが実行される。ナビゲータシーケンスＮＶは、トレーニングスキャンにおいて実行されるナビゲータシーケンスＮＴと同様に、ナビゲータ領域Ｒna（図３参照）から、横隔膜の位置ｕを検出するのに必要なデータを収集するためのシーケンスである。 Further, before the imaging sequence IS b are executed one or more times of navigator sequence NV is executed. The navigator sequence NV is a sequence for collecting data necessary for detecting the position u of the diaphragm from the navigator region Rna (see FIG. 3), similarly to the navigator sequence NT executed in the training scan.
図７は、本スキャンの説明図である。
図７（ａ）は、横隔膜の変位を表すグラフである。
グラフの横軸は時間を表し、グラフの縦軸は横隔膜位置ｕを示している。グラフの黒丸は、ナビゲータシーケンスＮＶにより検出された横隔膜位置ｕを表している。
FIG. 7 is an explanatory diagram of the main scan.
FIG. 7A is a graph showing the displacement of the diaphragm.
The horizontal axis of the graph represents time, and the vertical axis of the graph represents the diaphragm position u. The black circle in the graph represents the diaphragm position u detected by the navigator sequence NV.
図７（ｂ）は、横隔膜が基準位置ｕrefに変位したときのイメージング領域Ｒimの画像の概略図であり、図７（ｃ）は、ｎ回目の位相エンコード（イメージングシーケンスＩＳｎ）を実行したときのイメージング領域Ｒimの画像の概略図が示されている。 FIG. 7B is a schematic diagram of an image in the imaging region Rim when the diaphragm is displaced to the reference position u ref , and FIG. 7C is a diagram illustrating the execution of n-th phase encoding (imaging sequence IS n ). A schematic diagram of the image of the imaging region Rim at the time is shown.
横隔膜が基準位置ｕrefに変位したときのイメージング領域Ｒimの位置（ｘ，ｙ）における画像データをｆ（ｘ，ｙ）とし（図７（ｂ）参照）、ｎ回目の位相エンコードにおけるイメージング領域Ｒimの位置（ｘ，ｙ）における画像データをｇｎ（ｘ，ｙ）とする（図７（ｃ）参照）。ｆ（ｘ，ｙ）は、Δｘ（ｕ，ｘ，ｙ）およびΔｙ（ｕ，ｘ，ｙ）を用いて、以下の式（１）で表される。
ｇｎ（ｘ，ｙ）
＝ｆ（ｘ＋Δｘ（ｕｎ，ｘ，ｙ），ｙ＋Δｙ（ｕｎ，ｘ，ｙ）） ・・・（１）
ただし、ｕｎ：ｎ回目の位相エンコードの際の横隔膜位置
The image data at the position (x, y) of the imaging region Rim when the diaphragm is displaced to the reference position u ref is f (x, y) (see FIG. 7B), and the imaging region Rim in the n-th phase encoding. The image data at the position (x, y) is defined as g n (x, y) (see FIG. 7C). f (x, y) is expressed by the following equation (1) using Δx (u, x, y) and Δy (u, x, y).
g n (x, y)
= F (x + Δx (u n , x, y), y + Δy (u n , x, y)) (1)
Where u n : diaphragm position at the time of n-th phase encoding
ｇｎ（ｘ，ｙ）をフーリエ変換したものが、ｋ空間のデータＳ（ｔ，ｎ）となるので、ｋ空間のデータＳ（ｔ，ｎ）は、以下の式で表される。
尚、式（２）では、画像データｆ（ｘ，ｙ）をベクトルｆで表し、体動補正を受けたフーリエ変換のカーネルを行列Ｋで表した。式（２）において、ｋ空間のデータＳ（ｔ，ｎ）をベクトルＳで表すと、以下の式が得られる。
ｆ＝Ｋ−１Ｓ ・・・（３）
In Expression (2), the image data f (x, y) is represented by a vector f, and the kernel of Fourier transform subjected to body motion correction is represented by a matrix K. In equation (2), when k-space data S (t, n) is represented by a vector S, the following equation is obtained.
f = K −1 S (3)
式（３）において、ベクトルＳによって表されるｋ空間のデータＳ（ｔ，ｎ）は、本スキャンのイメージングシーケンスＩＳｂを実行することにより得られるデータである。また、逆行列Ｋ−１は、変位量Δｘ（ｕｎ，ｘ，ｙ）およびΔｙ（ｕｎ，ｘ，ｙ）が分かれば、求めることができる。したがって、変位量Δｘ（ｕｎ，ｘ，ｙ）およびΔｙ（ｕｎ，ｘ，ｙ）を求めることができれば、式（３）により、ベクトルｆで表される画像データｆ（ｘ，ｙ）を算出することができる。横隔膜位置ｕｎにおける変位量Δｘ（ｕｎ，ｘ，ｙ）およびΔｙ（ｕｎ，ｘ，ｙ）を求める方法としては、例えば、以下の２つの方法（１）および（２）が考えられる。 In the formula (3), data S k space represented by a vector S (t, n) is data obtained by executing the imaging sequence IS b of the scan. Further, the inverse matrix K −1 can be obtained if the displacement amounts Δx (u n , x, y) and Δy (u n , x, y) are known. Therefore, if the displacement amounts Δx (u n , x, y) and Δy (u n , x, y) can be obtained, the image data f (x, y) represented by the vector f is obtained from the equation (3). Can be calculated. Displacement Δx at the diaphragm position u n (u n, x, y) and [Delta] y (u n, x, y) as a method for determining the, for example, two methods (1) and (2) are considered below.
（方法１）先ず、ｎ回目の位相エンコードの際の横隔膜位置ｕｎを求める。横隔膜位置ｕｎを求める方法としては、例えば、ナビゲータシーケンスＮＶｎにより検出された横隔膜位置ｕｖｎと、ナビゲータシーケンスＮＶｎ＋１により検出された横隔膜位置ｕｖｎ＋１との平均値として求める方法や、横隔膜位置ｕｖｎ又はｕｖｎ＋１を、横隔膜位置ｕｎとして採用する方法が考えられる。
横隔膜位置ｕｎを求めた後、トレーニングスキャンにより求められた横隔膜位置ｕ’ｋ（ｋ＝１〜ｚ）における変位量Δｘ（ｕ’ｋ，ｘ，ｙ）およびΔｙ（ｕ’ｋ，ｘ，ｙ）（図６参照）に基づいて、横隔膜位置ｕｎにおける変位量Δｘ（ｕｎ，ｘ，ｙ）およびΔｙ（ｕｎ，ｘ，ｙ）を算出する。算出方法としては、図６に示す変位量Δｘ（ｕ’ｋ，ｘ，ｙ）およびΔｙ（ｕ’ｋ，ｘ，ｙ）の中から、ｕ’ｋがｕｎに近い変位量Δｘ（ｕ’ｋ，ｘ，ｙ）およびΔｙ（ｕ’ｋ，ｘ，ｙ）を用いて、変位量Δｘ（ｕｎ，ｘ，ｙ）およびΔｙ（ｕｎ，ｘ，ｙ）を算出する方法が考えられる。図７（ａ）を参照すると、ｕｎは、横隔膜位置ｕ’αとｕ’α−１との間に位置しているので、横隔膜位置ｕ’αとｕ’α−１は、ｕｎに近い値を有している。したがって、横隔膜位置ｕ’αにおける変位量Δｘ（ｕ’α，ｘ，ｙ）およびΔｙ（ｕ’α，ｘ，ｙ）と、横隔膜位置ｕ’α−１における変位量Δｘ（ｕ’α−１，ｘ，ｙ）およびΔｙ（ｕ’α−１，ｘ，ｙ）とに基づいて、横隔膜位置ｕｎにおける変位量Δｘ（ｕｎ，ｘ，ｙ）およびΔｙ（ｕｎ，ｘ，ｙ）を算出することができる。この場合、例えば、横隔膜位置ｕ’αにおける変位量Δｘ（ｕ’α，ｘ，ｙ）およびΔｙ（ｕ’α，ｘ，ｙ）と、横隔膜位置ｕ’α−１における変位量Δｘ（ｕ’α−１，ｘ，ｙ）およびΔｙ（ｕ’α−１，ｘ，ｙ）との平均値を、横隔膜位置ｕｎにおける変位量Δｘ（ｕｎ，ｘ，ｙ）およびΔｙ（ｕｎ，ｘ，ｙ）として算出することができる。尚、ｕ’αとｕｎとの差、およびｕ’α−１とｕｎとの差に応じて、変位量Δｘ（ｕ’α，ｘ，ｙ），Δｙ（ｕ’α，ｘ，ｙ）、および変位量Δｘ（ｕ’α−１，ｘ，ｙ），Δｙ（ｕ’α−１，ｘ，ｙ）に重み付けしてもよい。重み付けられた変位量に基づいて、横隔膜位置ｕｎにおける変位量Δｘ（ｕｎ，ｘ，ｙ）およびΔｙ（ｕｎ，ｘ，ｙ）を算出することにより、変位量Δｘ（ｕｎ，ｘ，ｙ）およびΔｙ（ｕｎ，ｘ，ｙ）の値の精度を高めることができる。
(Method 1) First, the diaphragm position un in the n-th phase encoding is obtained. As a method of obtaining the diaphragm position u n, for example, a method of obtaining an average value of the diaphragm position uv n detected by the navigator sequence NV n, detected by the navigator sequence NV n + 1 and the diaphragm position uv n + 1, the diaphragm position uv n or uv n + 1, is considered a method of adopting a diaphragm position u n.
After obtaining the diaphragm position u n , the displacement amounts Δx (u ′ k , x, y) and Δy (u ′ k , x, y) at the diaphragm position u ′ k (k = 1 to z) obtained by the training scan. ) (based on FIG. 6), calculates the amount of displacement [Delta] x (u n in diaphragm position u n, x, y) and [Delta] y (u n, x, y) a. As calculation method, the amount of displacement Δx shown in FIG. 6 (u 'k, x, y) and Δy (u' k, x, y) from the, u 'k is close displacement Δx to u n (u' A method of calculating the displacement amounts Δx (u n , x, y) and Δy (u n , x, y) using k , x, y) and Δy (u ′ k , x, y) is conceivable. Referring to FIG. 7 (a), u n Because located between the diaphragm position u 'alpha and u' alpha-1, the diaphragm position u 'alpha and u' alpha-1 is the u n Has a close value. Accordingly, the diaphragm position u 'displacement at α Δx (u' α, x , y) and Δy (u 'α, x, y) and, diaphragm position u' displacement [Delta] x (u 'alpha-1 in alpha-1 , x, y) and Δy (u 'α-1, x, based on the y), the displacement amount [Delta] x (u n in diaphragm position u n, x, y) and [Delta] y (u n, x, y) and Can be calculated. In this case, for example, the diaphragm position u 'displacement at α Δx (u' α, x , y) and Δy (u 'α, x, y) and, diaphragm position u' displacement at alpha-1 [Delta] x (u ' α-1, x, y) and Δy (u 'α-1, x, the average value of the y), displacement at the diaphragm position u n Δx (u n, x , y) and Δy (u n, x , Y). Incidentally, according to the difference between the 'difference between the alpha and u n, and u' u alpha-1 and u n, the displacement amount Δx (u 'α, x, y), Δy (u' α, x, y ), And displacement amounts Δx (u ′ α−1 , x, y), Δy (u ′ α−1 , x, y) may be weighted. Based on the displacement amount weighted amount of displacement [Delta] x in the diaphragm position u n (u n, x, y) and Δy (u n, x, y ) by calculating the amount of displacement [Delta] x (u n, x, The accuracy of the values of y) and Δy (u n , x, y) can be increased.
（方法２） 先ず、横隔膜位置ｕｎを求める（横隔膜位置ｕｎを求める方法は、方法１と同じでよい）。横隔膜位置ｕｎを求めたら、横隔膜位置ｕ’１〜ｕ’ｚの中で、ｕｎに一番近い横隔膜位置（ここでは、横隔膜位置ｕ’αとする）における変位量Δｘ（ｕ’α，ｘ，ｙ），Δｙ（ｕ’α，ｘ，ｙ）（図６参照）を、横隔膜位置ｕｎにおける変位量Δｘ（ｕｎ，ｘ，ｙ）およびΔｙ（ｕｎ，ｘ，ｙ）として採用する。この場合、ｕ’αとｕｎとの差に応じて、変位量Δｘ（ｕ’α，ｘ，ｙ），Δｙ（ｕ’α，ｘ，ｙ）を重み付けしてもよい。 (Method 2) First, a diaphragm position u n (method of obtaining the diaphragm position u n may be the same as the method 1). After determining the diaphragm position u n, in the diaphragm position u '1 ~u' z, closest diaphragm positioned u n (here, the diaphragm position u 'and alpha) displacement [Delta] x (u in' alpha, x, y), adopted [Delta] y (u 'alpha, x, y) (see FIG. 6), the displacement amount [Delta] x (u n in diaphragm position u n, x, y) and [Delta] y (u n, x, a y) To do. In this case, 'according to the difference between the alpha and u n, the displacement amount Δx (u' u α, x , y), Δy (u 'α, x, y) may be weighted.
したがって、図４〜図７を参照しながら説明したように、トレーニングスキャンと、本スキャンとを実行することにより、画像データｆ（ｘ，ｙ）が得られることがわかる。 Therefore, as described with reference to FIGS. 4 to 7, it is understood that the image data f (x, y) can be obtained by executing the training scan and the main scan.
次に、画像データｆ（ｘ，ｙ）を算出するときのＭＲＩ装置１の処理フローについて説明する。 Next, a processing flow of the MRI apparatus 1 when calculating the image data f (x, y) will be described.
図８は、ＭＲＩ装置１の処理フローの一例を示す図である。
ステップＳ１では、画像データ算出手段１０３（図１参照）が、トレーニングスキャンにより得られたデータに基づいて、タギング画像を作成する。そして、第１の変位量算出手段１０１（図１参照）が、タギング画像に付されたタグの歪を解析することにより、変位量Δｘ（ｕ’ｋ，ｘ，ｙ）およびΔｙ（ｕ’ｋ，ｘ，ｙ）を算出する（図６参照）。変位量Δｘ（ｕ’ｋ，ｘ，ｙ）およびΔｙ（ｕ’ｋ，ｘ，ｙ）を算出した後、ステップＳ２に進む。
FIG. 8 is a diagram illustrating an example of a processing flow of the MRI apparatus 1.
In step S1, the image data calculation unit 103 (see FIG. 1) creates a tagging image based on the data obtained by the training scan. Then, the first displacement amount calculation means 101 (see FIG. 1) analyzes the distortion of the tag attached to the tagging image, so that the displacement amounts Δx (u ′ k , x, y) and Δy (u ′ k , X, y) is calculated (see FIG. 6). After calculating the displacement amount Δx (u ′ k , x, y) and Δy (u ′ k , x, y), the process proceeds to step S2.
ステップＳ２では、第２の変位量算出手段１０２（図１参照）が、本スキャンにおけるｎ回目の位相エンコーディングの際の変位量Δｘ（ｕｎ，ｘ，ｙ）およびΔｙ（ｕｎ，ｘ，ｙ）を算出する。変位量Δｘ（ｕｎ，ｘ，ｙ）およびΔｙ（ｕｎ，ｘ，ｙ）は、図７を参照しながら説明した方法で算出することができる。変位量Δｘ（ｕｎ，ｘ，ｙ）およびΔｙ（ｕｎ，ｘ，ｙ）を算出した後、ステップＳ３に進む。 In step S2, the second displacement amount calculation means 102 (see FIG. 1) determines the displacement amounts Δx (u n , x, y) and Δy (u n , x, y) during the n-th phase encoding in the main scan. ) Is calculated. The displacement amounts Δx (u n , x, y) and Δy (u n , x, y) can be calculated by the method described with reference to FIG. After calculating the displacement amount Δx (u n , x, y) and Δy (u n , x, y), the process proceeds to step S3.
ステップＳ３では、画像データ算出手段１０３が、ステップＳ２で算出した変位量Δｘ（ｕｎ，ｘ，ｙ）およびΔｙ（ｕｎ，ｘ，ｙ）と、本スキャンのイメージングシーケンスにより得られたｋ空間のデータＳ（ｔ，ｎ）とを、式（３）に代入する。これにより、画像データｆ（ｘ，ｙ）を得ることができる。 In step S3, the image data calculation unit 103 uses the displacement amounts Δx (u n , x, y) and Δy (u n , x, y) calculated in step S2 and the k-space obtained by the imaging sequence of the main scan. Data S (t, n) is substituted into equation (3). Thereby, the image data f (x, y) can be obtained.
第１の実施形態では、タギングデータに基づいて変位量Δｘ（ｕｎ，ｘ，ｙ）およびΔｙ（ｕｎ，ｘ，ｙ）を算出するので、イメージング領域の各位置（ｘ，ｙ）がどれだけ変位したかを求めることができる。したがって、被検体の呼吸による体動アーチファクトが低減されたＭＲ画像を得ることができる。 In the first embodiment, since the displacement amounts Δx (u n , x, y) and Δy (u n , x, y) are calculated based on the tagging data, which position (x, y) in the imaging region is determined. It can be determined whether it is displaced. Therefore, an MR image in which body motion artifacts due to breathing of the subject are reduced can be obtained.
また、本スキャンにおいては、体動が安定する期間を待たずにイメージングシーケンスを実行することができるので、スキャン時間を短縮することができる。更に、被検体は自由呼吸下で撮影をすることができるので、撮影中に、被検体に呼吸停止の指示をする必要もなく、被検体の負担を軽減することができる。 In the main scan, since the imaging sequence can be executed without waiting for a period in which the body movement is stabilized, the scan time can be shortened. Furthermore, since the subject can be imaged under free breathing, there is no need to instruct the subject to stop breathing during imaging, and the burden on the subject can be reduced.
更に、上記の説明では、コロナル画像の画像データについて説明されているが、本発明は、上記のサジタル画像、アキシャル画像、オブリーク画像など、任意の断面の画像データを算出する場合にも適用することができる。 Further, in the above description, the image data of the coronal image is described. However, the present invention is also applicable to the case of calculating image data of an arbitrary cross section such as the above sagittal image, axial image, oblique image. Can do.
尚、図４および図５では、横隔膜が基準位置ｕrefに変位したときにタギングシーケンスＴＳａおよびを実行している。しかし、横隔膜位置に対して複数の基準位置を設け、横隔膜が各基準位置に変位するたびにタギングシーケンスＴＳａを実行してもよい。以下では、横隔膜位置に対して２箇所の基準位置を設けた場合にタギングシーケンスＴＳａを実行する例について説明する。 In FIG. 4 and FIG. 5, running tagging sequence TS a and when the diaphragm is displaced to the reference position u ref. However, a plurality of reference positions with respect to diaphragm position, may perform the tagging sequence TS a each time the diaphragm is displaced to the reference position. In the following, an example to perform tagging sequence TS a described case in which the reference position of the two positions relative to the diaphragm position.
図９は、２箇所の基準位置ｕref１およびｕref２を設け、各基準位置ｕref１およびｕref２においてタギングシーケンスを実行する例を示す図である。 FIG. 9 is a diagram illustrating an example in which two reference positions u ref1 and u ref2 are provided and a tagging sequence is executed at each of the reference positions u ref1 and u ref2 .
図９において、グラフの横軸は時間を表し、グラフの縦軸は横隔膜位置ｕを示している。グラフの黒丸は、ナビゲータシーケンスＮＴにより検出された横隔膜位置ｕを表している。 In FIG. 9, the horizontal axis of the graph represents time, and the vertical axis of the graph represents the diaphragm position u. A black circle in the graph represents the diaphragm position u detected by the navigator sequence NT.
基準位置ｕref１においてタギングシーケンスを実行する場合、先ず、横隔膜位置ｕが、減少から増大に変化したか否かを判断する。図９では、時刻ｔ１において、横隔膜位置ｕが、減少から増大に変化している。時刻ｔ１以降も、ナビゲータシーケンスＮＴによって横隔膜位置ｕを検出し、横隔膜が基準位置ｕref１に変位したら、タギングシーケンスＴＳ１を実行する。タギングシーケンスＴＳ１を実行することによりイメージング領域Ｒimにタグが付される。タギングシーケンスＴＳ１を実行した後、再びナビゲータシーケンスＮＴを実行し、横隔膜位置ｕが、ｕ’α１、ｕ’α２、ｕ’α３、およびｕ’α４に変位したら、データ収集シーケンスＤＳを実行する。データ収集シーケンスＤＳを実行することにより、タギングシーケンスＴＳ１でタグが付されたイメージング領域Ｒimから、横隔膜がｕ’α１、ｕ’α２、ｕ’α３、およびｕ’α４に変位したときのｋ空間のデータＤα１１、Ｄα２１、Ｄα３１、およびＤα４１を収集することができる。 When the tagging sequence is executed at the reference position u ref1 , it is first determined whether or not the diaphragm position u has changed from decreasing to increasing. In FIG. 9, at time t1, the diaphragm position u changes from decreasing to increasing. Time t1 and beyond, to detect the diaphragm position u by navigator sequence NT, diaphragmatic When displaced to the reference position u ref1, executes the tagging sequence TS 1. Tag is attached to the imaging region Rim by performing the tagging sequence TS 1. After executing the tagging sequence TS 1 , the navigator sequence NT is executed again, and when the diaphragm position u is displaced to u ′ α1 , u ′ α2 , u ′ α3 , and u ′ α4 , the data acquisition sequence DS is executed. By executing the data acquisition sequence DS, the k-space when the diaphragm is displaced to u ′ α1 , u ′ α2 , u ′ α3 , and u ′ α4 from the imaging region Rim tagged with the tagging sequence TS 1 data D α11, D α21, it is possible to collect the D α31, and D α41.
一方、基準位置ｕref２においてタギングシーケンスを実行する場合、先ず、横隔膜位置ｕが、増大から減少に変化したか否かを判断する。図９では、時刻ｔ２において、横隔膜位置ｕが、減少から増大に変化している。時刻ｔ２以降も、横隔膜位置ｕを検出し、横隔膜が基準位置ｕref２に変位したら、タギングシーケンスＴＳ２を実行する。タギングシーケンスＴＳ２を実行することによりイメージング領域Ｒimにタグが付される。タギングシーケンスＴＳ２を実行した後、再びナビゲータシーケンスＮＴを実行し、横隔膜位置ｕが、ｕ’β１、ｕ’β２、およびｕ’β３に変位したら、データ収集シーケンスＤＳを実行する。データ収集シーケンスＤＳを実行することにより、タギングシーケンスＴＳ２でタグが付されたイメージング領域Ｒimから、横隔膜がｕ’β１、ｕ’β２、およびｕ’β３に変位したときのｋ空間のデータＤβ１１、Ｄβ２１、およびＤβ３１を収集することができる。 On the other hand, when the tagging sequence is executed at the reference position u ref2 , it is first determined whether or not the diaphragm position u has changed from increasing to decreasing. In FIG. 9, the diaphragm position u changes from decrease to increase at time t2. Time t2 beyond detects the diaphragm position u, diaphragmatic When displaced to the reference position u ref2, executes the tagging sequence TS 2. Tag is attached to the imaging region Rim by performing the tagging sequence TS 2. After executing the tagging sequence TS 2 , the navigator sequence NT is executed again. When the diaphragm position u is displaced to u ′ β1 , u ′ β2 and u ′ β3 , the data acquisition sequence DS is executed. By executing the data acquisition sequence DS, the k-space data D β11 when the diaphragm is displaced to u ′ β1 , u ′ β2 , and u ′ β3 from the imaging region Rim tagged with the tagging sequence TS 2. , D β21 , and D β31 can be collected.
以下、同様に、基準位置ｕref１およびｕref２が検出されたらタギングシーケンスを実行し、横隔膜位置ｕ’α１、ｕ’α２、ｕ’α３、およびｕ’α４におけるｋ空間のデータと、横隔膜位置ｕ’β１、ｕ’β２、およびｕ’β３におけるｋ空間のデータとを交互に収集する。このようにしてイメージング領域Ｒimからｋ空間のデータを収集してタギング画像を作成し、変位量Δｘ（ｕ，ｘ，ｙ）およびΔｙ（ｕ，ｘ，ｙ）を算出してもよい。ただし、基準位置ｕref１と基準位置ｕref２との間には位置のズレがあるので、基準位置ｕref２を基準にして算出された変位量は、基準位置ｕref１を基準にして算出された変位量に補正する必要がある。 Hereinafter, similarly, when the reference positions u ref1 and u ref2 are detected, a tagging sequence is executed, and k-space data at the diaphragm positions u ′ α1 , u ′ α2 , u ′ α3 , and u ′ α4, and the diaphragm position u The data of k-space in ' β1 , u'β2 , and u'β3 are collected alternately. In this way, k-space data may be collected from the imaging region Rim to create a tagging image, and the displacement amounts Δx (u, x, y) and Δy (u, x, y) may be calculated. However, since there is a positional deviation between the reference position u ref1 and the reference position u ref2 , the displacement calculated with reference to the reference position u ref2 is the displacement calculated with reference to the reference position u ref1. The amount needs to be corrected.
（２）第２の実施形態
第２の実施形態では、本スキャンにより収集されたｋ空間のデータを取り直すか否かを決定する方法について説明する。尚、ハードウェア構成は、第１の実施形態と同じである。
(2) Second Embodiment In the second embodiment, a method for determining whether to reacquire k-space data collected by the main scan will be described. The hardware configuration is the same as that of the first embodiment.
図１０は、第２の実施形態において、トレーニングスキャン時における横隔膜の変位を表すグラフである。グラフの黒丸は、トレーニングスキャン時のナビゲータシーケンスＮＴにより検出された横隔膜位置ｕを表している。 FIG. 10 is a graph showing the displacement of the diaphragm during a training scan in the second embodiment. The black circle in the graph represents the diaphragm position u detected by the navigator sequence NT during the training scan.
トレーニングスキャンにより検出された横隔膜位置に基づいて、被検体１３が息を吐き終わった時の横隔膜１３ｂの位置である最大呼気位置Ｘｅと、被検体１３が息を吸い終わった時の横隔膜１３ｂの位置である最大吸気位置Ｘｉとを決定する。 Based on the diaphragm position detected by the training scan, the maximum exhalation position Xe, which is the position of the diaphragm 13b when the subject 13 has exhaled, and the position of the diaphragm 13b when the subject 13 has finished breathing And the maximum intake position Xi is determined.
最大呼気位置Ｘｅを決定する方法としては、例えば、検出された横隔膜位置の中で、座標が最大（座標の正負の取り方によっては最小）となる位置を、最大呼気位置Ｘｅとして決定する方法が考えられる。また、最大吸気位置Ｘｉを決定する方法としては、例えば、検出された横隔膜位置の中で、座標が最小（座標の正負の取り方によっては最大）となる位置を、最大吸気位置Ｘｉとして決定する方法が考えられる。 As a method for determining the maximum exhalation position Xe, for example, there is a method for determining, as the maximum exhalation position Xe, a position where the coordinate is maximum (minimum depending on whether the coordinate is positive or negative) among the detected diaphragm positions. Conceivable. Further, as a method of determining the maximum intake position Xi, for example, a position where the coordinates are minimum (maximum depending on whether the coordinates are positive or negative) among the detected diaphragm positions is determined as the maximum intake position Xi. A method is conceivable.
最大呼気位置Ｘｅおよび最大吸気位置Ｘｉを決定した後、最大呼気位置Ｘｅおよび最大吸気位置Ｘｉに基づいて、横隔膜の位置の許容範囲ＡＷを設定する。ＡＷの幅は、例えば、最大呼気位置Ｘｅと最大吸気位置Ｘｉとの間の範囲に設定したり、最大呼気位置Ｘｅ±ΔＸｅ（ΔＸｅ＞０）と最大吸気位置Ｘｉ±ΔＸｉ（ΔＸｉ＞０）との間の範囲に設定することができる。
以上のようにして、横隔膜の位置の許容範囲ＡＷが設定される。
After determining the maximum exhalation position Xe and the maximum inhalation position Xi, an allowable range AW of the diaphragm position is set based on the maximum exhalation position Xe and the maximum inspiration position Xi. For example, the width of the AW is set to a range between the maximum expiration position Xe and the maximum inspiration position Xi, or the maximum expiration position Xe ± ΔXe (ΔXe> 0) and the maximum inspiration position Xi ± ΔXi (ΔXi> 0). Can be set to a range between.
As described above, the allowable range AW of the diaphragm position is set.
本スキャンでは、横隔膜位置の許容範囲ＡＷに基づいて、イメージングシーケンスＩＳｎで収集されたイメージングデータを取り直すか否かを決定する（図１１参照）。 In this scan, based on the allowable range AW diaphragm position determines whether take heart imaging data collected in the imaging sequence IS n (see FIG. 11).
図１１は、本スキャンにおける横隔膜の変位と、横隔膜位置の許容範囲ＡＷとの関係を示す図である。グラフの黒丸は、本スキャン時のナビゲータシーケンスＮＶにより検出された横隔膜位置ｕを表している。 FIG. 11 is a diagram illustrating the relationship between the displacement of the diaphragm and the allowable range AW of the diaphragm position in the main scan. The black circle in the graph represents the diaphragm position u detected by the navigator sequence NV during the main scan.
ｎ回目の位相エンコードにおける横隔膜位置ｕｎが、横隔膜位置の許容範囲ＡＷに含まれている場合、イメージングシーケンスＩＳｎで収集されたイメージングデータは、ｋ空間を埋めるデータとして採用される。一方、横隔膜位置ｕｎが、横隔膜位置の許容範囲ＡＷからはずれた場合、イメージングシーケンスＩＳｎで収集されたイメージングデータは、イメージングシーケンスＩＳｎ＋ｄ（ｄ：１以上の整数）おいて、データが取り直される。図１１では、ｎ回目の位相エンコードにおける横隔膜位置ｕｎは、横隔膜位置の許容範囲ＡＷに含まれているので、イメージングシーケンスＩＳｎで収集されたイメージングデータは、ｋ空間を埋めるデータとして採用される。 diaphragm position u n in the n-th phase encoding, if included in the allowable range AW diaphragm position, the imaging data collected in the imaging sequence IS n is employed as data to fill the k-space. On the other hand, the diaphragm position u n is, if out of the allowable range AW diaphragm position, the imaging data collected in the imaging sequence IS n is imaging sequence IS n + d (d: 1 or more integer) Oite, data Torinaosa It is. In Figure 11, n-th diaphragm position u n in the phase encoding, since it is included in the allowable range AW diaphragm position, the imaging data collected in the imaging sequence IS n is employed as data to fill the k-space .
一方、ｘ回目の位相エンコードにおけるイメージングシーケンスＩＳｘは、横隔膜位置の許容範囲ＡＷから外れている。したがって、イメージングシーケンスＩＳｘで収集されたイメージングデータは、ｋ空間を埋めるデータとして採用されず、イメージングシーケンスＩＳｘ＋ｄ（ｄ：１以上の整数）おいて、データが取り直される。したがって、第２の実施形態によれば、横隔膜位置の許容範囲ＡＷにおけるデータのみを、ｋ空間を埋めるデータとして採用することができ、体動アーチファクトが更に低減された画像を得ることができる。 On the other hand, the imaging sequence IS x in the x-th phase encoding is out of the allowable range AW of the diaphragm position. Thus, the imaging data collected in the imaging sequence IS x is not employed as data to fill the k-space, the imaging sequence IS x + d (d: 1 or more integer) Oite, data is Torinaosa. Therefore, according to the second embodiment, only data in the allowable range AW of the diaphragm position can be employed as data for filling the k space, and an image in which body motion artifacts are further reduced can be obtained.
尚、本発明の実施に際しては、上記した実施形態に限定されるものではなく、種々の変形形態を採用することができる。 In carrying out the present invention, the present invention is not limited to the above-described embodiment, and various modifications can be employed.
第１および第２の実施形態では、横隔膜位置は、ナビゲータエコー法で算出しているが、ナビゲータエコー法の代わりに、ベローズを用いて横隔膜位置を算出してもよい。 In the first and second embodiments, the diaphragm position is calculated by the navigator echo method, but the diaphragm position may be calculated by using a bellows instead of the navigator echo method.
第１および第２の実施形態では、ナビゲータエコーは、被検体のイメージング領域の内側から収集されている。しかし、体動アーチファクトが十分に低減された画像を得ることができるのであれば、ナビゲータエコーは、イメージング領域の外側に存在する部位から収集してもよい。 In the first and second embodiments, navigator echoes are collected from the inside of the imaging region of the subject. However, navigator echoes may be collected from a site existing outside the imaging region as long as an image with sufficiently reduced body motion artifacts can be obtained.
第１および第２の実施形態では、変位量Δｘ（ｕ，ｘ，ｙ）およびΔｙ（ｕ，ｘ，ｙ）を、横隔膜位置ｕの関数として求めている。しかし、イメージング領域Ｒimの位置（ｘ，ｙ）における変位量を求めることができるのであれば、横隔膜とは別の部位の位置の関数として、変位量Δｘ（ｕ，ｘ，ｙ）およびΔｙ（ｕ，ｘ，ｙ）を求めてもよい。 In the first and second embodiments, the displacement amounts Δx (u, x, y) and Δy (u, x, y) are obtained as a function of the diaphragm position u. However, if the displacement amount at the position (x, y) of the imaging region Rim can be obtained, the displacement amounts Δx (u, x, y) and Δy (u) as a function of the position of a part other than the diaphragm. , X, y) may be obtained.
また、第１および第２の実施形態では、式（３）に従って画像データｆ（ｘ，ｙ）を算出しているが、式（３）とは別の式に従って画像データｆ（ｘ，ｙ）を算出してもよい。 In the first and second embodiments, the image data f (x, y) is calculated according to the equation (3), but the image data f (x, y) according to an equation different from the equation (3). May be calculated.
更に、第１および第２の実施形態では、肝臓１３ａを撮影する例について説明されているが、被検体の体動に伴って変位する部位であれば、心臓などの他の部位であってもよい。 Furthermore, in the first and second embodiments, an example in which the liver 13a is imaged has been described. However, any other part such as the heart may be used as long as it is a part that is displaced with the body movement of the subject. Good.
１ ＭＲＩ装置
２ 磁場発生装置
３ テーブル
４ クレードル
５ 受信コイル
６ シーケンサ
７ 送信器
８ 勾配磁場電源
９ 受信器
１０ 中央処理装置
１１ 入力装置
１２ 表示装置
１３ 被検体
１４ オペレータ
２１ ボア
２２ 超伝導コイル
２３ 勾配コイル
２４ 送信コイル
DESCRIPTION OF SYMBOLS 1 MRI apparatus 2 Magnetic field generator 3 Table 4 Cradle 5 Receiving coil 6 Sequencer 7 Transmitter 8 Gradient magnetic field power supply 9 Receiver 10 Central processing unit 11 Input device 12 Display device 13 Subject 14 Operator 21 Bore 22 Superconducting coil 23 Gradient coil 24 Transmitting coil
Claims (12)
前記所定の部位にタグが付されたタギング画像のデータを収集するためのトレーニングスキャンを実行するスキャン手段と、
前記タギング画像のデータに基づいて、前記所定の部位の変位量を算出する第１の変位量算出手段と、
前記第１の変位量算出手段により算出された前記所定の部位の変位量に基づいて、前記本スキャンのｎ回目の位相エンコードにおける前記所定の部位の変位量を算出する第２の変位量算出手段と、
前記第２の変位量算出手段が算出した前記所定の部位の変位量と、前記本スキャンにより得られるｋ空間のデータとに基づいて、前記所定の部位の画像データを算出する画像データ算出手段と、
を有する、磁気共鳴イメージング装置。 A magnetic resonance imaging apparatus that performs a main scan for applying a gradient magnetic field in a phase encoding direction and imaging a predetermined portion that is displaced in accordance with a body motion of a subject,
Scanning means for executing a training scan for collecting tagging image data with a tag attached to the predetermined site;
First displacement amount calculating means for calculating a displacement amount of the predetermined part based on the tagging image data;
Second displacement amount calculation means for calculating the displacement amount of the predetermined portion in the n-th phase encoding of the main scan based on the displacement amount of the predetermined portion calculated by the first displacement amount calculation means. When,
Image data calculating means for calculating image data of the predetermined part based on the displacement amount of the predetermined part calculated by the second displacement amount calculating means and k-space data obtained by the main scan; ,
A magnetic resonance imaging apparatus.
前記タギング画像の各ピクセル位置ごとに前記変位量を算出する、請求項１に記載の磁気共鳴イメージング装置。 The first displacement amount calculating means includes:
The magnetic resonance imaging apparatus according to claim 1, wherein the displacement amount is calculated for each pixel position of the tagging image.
前記被検体の体動に伴って変位する第１の部位の位置を検出するためのシーケンスが実行される、請求項１又は２に記載の磁気共鳴イメージング装置。 In the training scan and the main scan,
The magnetic resonance imaging apparatus according to claim 1, wherein a sequence for detecting a position of a first part that is displaced in accordance with a body movement of the subject is executed.
前記所定の部位の変位量を、前記第１の部位の位置の関数として算出し、
前記第２の変位量算出手段は、
前記本スキャンのｎ回目の位相エンコードにおける前記所定の部位の変位量を、前記第１の部位の位置の関数として算出する、請求項３に記載の磁気共鳴イメージング装置。 The first displacement amount calculating means includes:
Calculating the amount of displacement of the predetermined part as a function of the position of the first part;
The second displacement amount calculating means includes:
The magnetic resonance imaging apparatus according to claim 3, wherein a displacement amount of the predetermined part in the n-th phase encoding of the main scan is calculated as a function of the position of the first part.
前記所定の部位にタグを付すためのタギングシーケンスと、
タグが付された前記所定の部位からデータを収集するデータ収集シーケンスと、
前記第１の部位の位置を検出するためのナビゲータシーケンスと、を有する、請求項３〜５のうちのいずれか一項に記載の磁気共鳴イメージング装置。 The training scan is
A tagging sequence for tagging the predetermined site;
A data collection sequence for collecting data from the predetermined part to which the tag is attached;
The magnetic resonance imaging apparatus according to claim 3, further comprising: a navigator sequence for detecting a position of the first part.
位相エンコード方向に勾配磁場を印加し、前記所定の部位からｋ空間のデータを収集するためのイメージングシーケンスと、
前記第１の部位の位置を検出するためのナビゲータシーケンスと、を有する、請求項３〜６のうちのいずれか一項に記載の磁気共鳴イメージング装置。 The main scan is
An imaging sequence for applying a gradient magnetic field in the phase encoding direction and collecting k-space data from the predetermined site;
The magnetic resonance imaging apparatus according to claim 3, further comprising: a navigator sequence for detecting a position of the first part.
前記本スキャンにより得られた前記所定の部位のｋ空間のデータと、体動補正を受けたフーリエ変換のカーネルの逆行列とに基づいて、前記所定の部位の画像データを算出する、請求項１〜１０のうちのいずれか一項に記載の磁気共鳴イメージング装置。 The image data calculation means includes
The image data of the predetermined part is calculated based on k-space data of the predetermined part obtained by the main scan and an inverse matrix of a Fourier transform kernel subjected to body motion correction. 10. The magnetic resonance imaging apparatus according to any one of 10.
前記第２の変位量算出手段が算出した前記所定の部位の変位量に基づいて、前記体動補正を受けたフーリエ変換のカーネルの逆行列を算出する、請求項１０に記載の磁気共鳴イメージング装置。 The image data calculation means includes
11. The magnetic resonance imaging apparatus according to claim 10, wherein an inverse matrix of a Fourier transform kernel subjected to the body motion correction is calculated based on a displacement amount of the predetermined part calculated by the second displacement amount calculation unit. .
前記トレーニングスキャンにより得られた磁気共鳴信号に基づいて、前記タギング画像のデータを算出する、請求項１〜１１のうちのいずれか一項に記載の磁器K表明イメージング装置。
The image data calculation means includes
The porcelain K manifest imaging apparatus according to any one of claims 1 to 11, wherein data of the tagging image is calculated based on a magnetic resonance signal obtained by the training scan.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2010042679A JP5627903B2 (en) | 2010-02-26 | 2010-02-26 | Magnetic resonance imaging system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2010042679A JP5627903B2 (en) | 2010-02-26 | 2010-02-26 | Magnetic resonance imaging system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2011177245A true JP2011177245A (en) | 2011-09-15 |
| JP5627903B2 JP5627903B2 (en) | 2014-11-19 |
Family
ID=44689480
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2010042679A Active JP5627903B2 (en) | 2010-02-26 | 2010-02-26 | Magnetic resonance imaging system |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP5627903B2 (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2013132333A (en) * | 2011-12-26 | 2013-07-08 | Ge Medical Systems Global Technology Co Llc | Magnetic resonance apparatus and program |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000296120A (en) * | 1999-04-13 | 2000-10-24 | Hitachi Medical Corp | Magnetic resonance image diagnostic system |
| WO2004080301A1 (en) * | 2003-03-14 | 2004-09-23 | Hitachi Medical Corporation | Magnetic resonance imaging device |
| JP2005040416A (en) * | 2003-07-24 | 2005-02-17 | Hitachi Medical Corp | Magnetic resonance imaging apparatus |
| JP2007190114A (en) * | 2006-01-18 | 2007-08-02 | Hitachi Medical Corp | Magnetic resonance imaging apparatus |
-
2010
- 2010-02-26 JP JP2010042679A patent/JP5627903B2/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000296120A (en) * | 1999-04-13 | 2000-10-24 | Hitachi Medical Corp | Magnetic resonance image diagnostic system |
| WO2004080301A1 (en) * | 2003-03-14 | 2004-09-23 | Hitachi Medical Corporation | Magnetic resonance imaging device |
| JP2005040416A (en) * | 2003-07-24 | 2005-02-17 | Hitachi Medical Corp | Magnetic resonance imaging apparatus |
| JP2007190114A (en) * | 2006-01-18 | 2007-08-02 | Hitachi Medical Corp | Magnetic resonance imaging apparatus |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2013132333A (en) * | 2011-12-26 | 2013-07-08 | Ge Medical Systems Global Technology Co Llc | Magnetic resonance apparatus and program |
Also Published As
| Publication number | Publication date |
|---|---|
| JP5627903B2 (en) | 2014-11-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4807825B2 (en) | Magnetic resonance imaging system | |
| JP5105848B2 (en) | Magnetic resonance imaging apparatus and imaging condition setting method in magnetic resonance imaging apparatus | |
| EP3633401A1 (en) | Prevention of compensating a wrongly detected motion in mri | |
| WO2006064582A1 (en) | Magnetic resonance imaging system | |
| CN103271740B (en) | Magnetic resonance imaging method employing and system | |
| JP6408594B2 (en) | Time-efficient 4D magnetic resonance imaging | |
| JP4953823B2 (en) | Magnetic resonance imaging method and apparatus | |
| JP5023127B2 (en) | Magnetic resonance imaging system | |
| WO2014027547A1 (en) | Magnetic resonance imaging device and magnetic resonance imaging method | |
| JPWO2004080301A1 (en) | Magnetic resonance imaging system | |
| JP2011147561A (en) | Magnetic resonance imaging apparatus | |
| JP5106697B2 (en) | Magnetic resonance imaging system | |
| JP2006320527A (en) | Magnetic resonance imaging system | |
| KR20170102656A (en) | Magnetic resonance imaging apparatus and scanning method for magnetic resonance image thereof | |
| JP2004209084A (en) | Inspection apparatus using nuclear magnetic resonance | |
| JP4558397B2 (en) | Magnetic resonance imaging system | |
| JP6109598B2 (en) | Magnetic resonance imaging system | |
| JP5627903B2 (en) | Magnetic resonance imaging system | |
| JP6356412B2 (en) | Magnetic resonance apparatus and program | |
| KR102232606B1 (en) | Device and method for dynamic tagged magnet resonance imaging | |
| JP5162126B2 (en) | Magnetic resonance imaging system | |
| JP6121664B2 (en) | Magnetic resonance apparatus and program | |
| JP2013154015A (en) | Magnetic resonance device and program | |
| JP5477948B2 (en) | Magnetic resonance imaging system | |
| JP6325267B2 (en) | Magnetic resonance imaging system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A625 | Written request for application examination (by other person) |
Free format text: JAPANESE INTERMEDIATE CODE: A625 Effective date: 20120604 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20131120 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140218 |
|
| 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: 20140901 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20141001 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 5627903 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |