JP6121664B2 - Magnetic resonance apparatus and program - Google Patents
Magnetic resonance apparatus and program Download PDFInfo
- Publication number
- JP6121664B2 JP6121664B2 JP2012166187A JP2012166187A JP6121664B2 JP 6121664 B2 JP6121664 B2 JP 6121664B2 JP 2012166187 A JP2012166187 A JP 2012166187A JP 2012166187 A JP2012166187 A JP 2012166187A JP 6121664 B2 JP6121664 B2 JP 6121664B2
- Authority
- JP
- Japan
- Prior art keywords
- signal
- navigator
- imaging data
- data
- magnetic resonance
- 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.)
- Active
Links
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Description
本発明は、所定の部位を撮影する磁気共鳴装置、およびその磁気共鳴装置に適用されるプログラムに関する。 The present invention relates to a magnetic resonance apparatus for imaging a predetermined part, and a program applied to the magnetic resonance apparatus.
磁気共鳴装置を用いて被検体の診断を行う方法として、被検体に造影剤を投与し、造影剤を投与した後に現れる各時相の画像データを取得する方法が知られている。 As a method for diagnosing a subject using a magnetic resonance apparatus, a method is known in which a contrast agent is administered to a subject and image data of each time phase appearing after the contrast agent is administered is acquired.
近年、肝臓のMRI検査ではガドキセト酸ナトリウム(Gd-EOB-DTPA:以下、EOBと呼ぶ)を用いた造影ダイナミック撮像検査が頻繁に行われるようになったが、従前より使用されていたガドリニウムキレート(Gd-DTPA)と比較して、投与量が半分である。したがって、EOBを用いた撮像では造影剤の濃度が一定になる時間が短く、リンギングやブラーなどのアーチファクトが現れやすいという問題がある。このため、リンギングやブラーなどのアーチファクトを低減することが望まれている。 Recently, contrast-enhanced dynamic imaging tests using sodium gadoxetate (Gd-EOB-DTPA: hereinafter referred to as EOB) have been frequently performed in liver MRI examinations, but gadolinium chelates ( Compared with Gd-DTPA), the dose is halved. Therefore, in the imaging using EOB, there is a problem that the time for which the concentration of the contrast agent becomes constant is short, and artifacts such as ringing and blurring are likely to appear. For this reason, it is desired to reduce artifacts such as ringing and blur.
本発明の第1の態様は、被検体の所定の部位を含む第1の領域からナビゲータ信号を収集するとともに、前記所定の部位を含む第2の領域から前記所定の部位の画像を再構成するときに使用されるイメージングデータを収集する磁気共鳴装置であって、
前記ナビゲータ信号と前記イメージングデータとを収集するためのスキャンを実行するスキャン手段と、
前記ナビゲータ信号に基づいて、前記イメージングデータの信号強度を補正する補正手段と、
信号強度が補正されたイメージングデータに基づいて画像再構成を行う再構成手段と、
を有する磁気共鳴装置である。
In the first aspect of the present invention, a navigator signal is collected from a first region including a predetermined part of a subject, and an image of the predetermined part is reconstructed from a second region including the predetermined part. A magnetic resonance apparatus for collecting imaging data used sometimes,
Scanning means for executing a scan for collecting the navigator signal and the imaging data;
Correction means for correcting the signal intensity of the imaging data based on the navigator signal;
Reconstruction means for performing image reconstruction based on imaging data whose signal intensity is corrected;
Is a magnetic resonance apparatus.
本発明の第2の態様は、被検体の所定の部位を含む第1の領域からナビゲータ信号を収集するとともに、前記所定の部位を含む第2の領域から前記所定の部位の画像を再構成するときに使用されるイメージングデータを収集する磁気共鳴装置であって、
前記ナビゲータ信号と前記イメージングデータとを収集するためのスキャンを実行するスキャン手段と、
前記ナビゲータ信号に基づいて、前記イメージングデータを収集するときに使用されるシーケンスのRFパルスのフリップ角を補正する補正手段と、
を有する、磁気共鳴装置である。
In the second aspect of the present invention, navigator signals are collected from a first region including a predetermined part of the subject, and an image of the predetermined part is reconstructed from the second region including the predetermined part. A magnetic resonance apparatus for collecting imaging data used sometimes,
Scanning means for executing a scan for collecting the navigator signal and the imaging data;
Correction means for correcting a flip angle of an RF pulse of a sequence used when acquiring the imaging data based on the navigator signal;
This is a magnetic resonance apparatus.
本発明の第3の態様は、被検体の所定の部位を含む第1の領域からナビゲータ信号を収集するとともに、前記所定の部位を含む第2の領域から前記所定の部位の画像を再構成するときに使用されるイメージングデータを収集するためのスキャンを実行する磁気共鳴装置のプログラムであって、
前記ナビゲータ信号に基づいて、前記イメージングデータの信号強度を補正する補正処理と、
信号強度が補正されたイメージングデータに基づいて画像再構成を行う再構成処理と、
を計算機に実行させるためのプログラムである。
In a third aspect of the present invention, navigator signals are collected from a first region including a predetermined part of the subject, and an image of the predetermined part is reconstructed from the second region including the predetermined part. A program of a magnetic resonance apparatus that performs a scan to collect imaging data that is sometimes used,
Correction processing for correcting the signal intensity of the imaging data based on the navigator signal;
Reconstruction processing to perform image reconstruction based on imaging data whose signal intensity is corrected;
Is a program for causing a computer to execute.
本発明の第4の態様は、被検体の所定の部位を含む第1の領域からナビゲータ信号を収集するとともに、前記所定の部位を含む第2の領域から前記所定の部位の画像を再構成するときに使用されるイメージングデータを収集するためのスキャンを実行する磁気共鳴装置のプログラムであって、
前記ナビゲータ信号に基づいて、前記イメージングデータを収集するときに使用されるシーケンスのRFパルスのフリップ角を補正する補正処理、
を計算機に実行させるためのプログラムである。
According to a fourth aspect of the present invention, navigator signals are collected from a first region including a predetermined part of the subject, and an image of the predetermined part is reconstructed from the second region including the predetermined part. A program of a magnetic resonance apparatus that performs a scan to collect imaging data that is sometimes used,
Correction processing for correcting a flip angle of an RF pulse of a sequence used when acquiring the imaging data based on the navigator signal;
Is a program for causing a computer to execute.
信号強度に基づいてイメージングデータ又はフリップ角を補正するので、リンギングやブラーなどのアーチファクトを低減することができる。 Since the imaging data or the flip angle is corrected based on the signal intensity, artifacts such as ringing and blur can be reduced.
以下、発明を実施するための形態について説明するが、本発明は、以下の形態に限定されることはない。 Hereinafter, although the form for inventing is demonstrated, this invention is not limited to the following forms.
(1)第1の形態
図1は、本発明の第1の形態の磁気共鳴装置の概略図である。
磁気共鳴装置(以下、「MR装置」と呼ぶ。MR:Magnetic Resonance)100は、マグネット2、テーブル3、受信コイル4などを有している。
(1) First Embodiment FIG. 1 is a schematic view of a magnetic resonance apparatus according to a first embodiment of the present invention.
A magnetic resonance apparatus (hereinafter referred to as “MR apparatus”, MR: Magnetic Resonance) 100 includes a
マグネット2は、被検体12が収容されるボア21と、超伝導コイル22と、勾配コイル23と、RFコイル24とを有している。超伝導コイル22は静磁場を印加し、勾配コイル23は勾配磁場を印加し、RFコイル24はRFパルスを送信する。尚、超伝導コイル22の代わりに、永久磁石を用いてもよい。
The
テーブル3は、クレードル3aを有している。クレードル3aは、ボア21内に移動できるように構成されている。クレードル3aによって、被検体12はボア21に搬送される。
The table 3 has a
受信コイル4は、被検体12の腹部に取り付けられている。受信コイル4は、被検体12からの磁気共鳴信号を受信する。
The
MR装置100は、更に、造影剤注入装置5、送信器6、勾配磁場電源7、受信器8、制御部9、操作部10、および表示部11などを有している。
The
造影剤注入装置5は、被検体に造影剤を注入する。
送信器6はRFコイル24に電流を供給し、勾配磁場電源7は勾配コイル23に電流を供給する。
受信器8は、受信コイル4から受け取った信号に対して、検波などの信号処理を実行する。
The contrast
The
The
制御部9は、表示部11に必要な情報を伝送したり、受信器8から受け取ったデータに基づいて画像を再構成するなど、MR装置100の各種の動作を実現するように、MR装置100の各部の動作を制御する。制御部9は、例えばコンピュータ(computer)によって構成される。制御部9は、補正手段91および再構成手段92などを有している。
The control unit 9 transmits necessary information to the
補正手段91は、イメージングデータの信号強度を補正する。補正手段91は、信号強度データ作成手段91a、検出手段91b、演算手段91cなどを有している。
The correcting unit 91 corrects the signal intensity of the imaging data. The correction unit 91 includes a signal intensity data creation unit 91a, a
信号強度データ作成手段91aは、ナビゲータ信号Anに基づいて、ナビゲータ領域Rnav内のSI方向の信号強度の変化を表す信号強度データDnを作成する(例えば、図5参照)。
検出手段91bは、信号強度データに基づいて肝臓の位置における信号強度を検出する。
演算手段91cは、後述する式を用いて、補正後のイメージングデータを求める。
Signal strength data generating means 91a, based on the navigator signal A n, to create a signal strength data D n representing the variation of the SI direction of the signal strength of the navigator region R nav (e.g., see FIG. 5).
The detection means 91b detects the signal intensity at the position of the liver based on the signal intensity data.
The
再構成手段92は、補正後のイメージングデータを用いて画像を再構成する。
The
制御部9は、補正手段91および再構成手段92を構成する一例であり、所定のプログラムを実行することにより、これらの手段として機能する。
The control unit 9 is an example constituting the correcting unit 91 and the reconfiguring
操作部10は、オペレータにより操作され、種々の情報を制御部9に入力する。表示部11は種々の情報を表示する。
MR装置100は、上記のように構成されている。
The operation unit 10 is operated by an operator and inputs various information to the control unit 9. The
The
図2は第1の形態で実行されるスキャンを示す図、図3は撮影領域を概略的に示す図である。
第1の形態では、息止めスキャンBHが実行される。
息止めスキャンBHは、被検体に造影剤を注入し、被検体が息止めをした状態で実行されるスキャンである。息止めスキャンBHでは、ナビゲータシーケンスNSn(n=1〜z)と、イメージングシーケンスISn(n=1〜z)とが交互に実行される。
FIG. 2 is a diagram showing a scan executed in the first form, and FIG. 3 is a diagram schematically showing an imaging region.
In the first mode, a breath holding scan BH is executed.
The breath-hold scan BH is a scan that is executed in a state where a contrast medium is injected into the subject and the subject holds the breath. In the breath holding scan BH, the navigator sequence NS n (n = 1 to z) and the imaging sequence IS n (n = 1 to z) are alternately executed.
ナビゲータシーケンスNSnは、肺と肝臓との境界部分を含むナビゲータ領域Rnav内のSI方向の信号強度の変化を表す信号強度データを取得するためのシーケンスである。イメージングシーケンスISnは肝臓を含む撮影領域Rの画像を再構成するときに使用されるイメージングデータを取得するためのシーケンスである。
以下に、息止めスキャンBHを実行するときのフローについて説明する。
The navigator sequence NS n is a sequence for acquiring signal strength data representing a change in signal strength in the SI direction in the navigator region R nav including the boundary between the lung and the liver. The imaging sequence IS n is a sequence for acquiring imaging data used when reconstructing an image of the imaging region R including the liver.
The flow when executing the breath holding scan BH will be described below.
図4は、第1の形態において被検体を撮影するときのフローを示す図である。
ステップST1では、被検体に造影剤を注入し、息止めスキャンBHを実行する。息止めスキャンBHでは、図5に示すように、ナビゲータシーケンスNSnとイメージングシーケンスISnとが交互に実行される。ナビゲータシーケンスNSnを実行することにより、ナビゲータ領域Rnavからナビゲータ信号An(n=1〜z)が収集され、イメージングシーケンスISnを実行することにより、撮影領域Rの画像を再構成するときに使用されるイメージングデータSn(n=1〜z)が収集される。信号強度データ作成手段91a(図1参照)は、ナビゲータ信号Anを一次元フーリエ変換(1D−FT)し、ナビゲータ領域Rnav内のSI方向の信号強度の変化を表す信号強度データD1〜Dzを作成する。信号強度データD1〜Dzを得た後、ステップST2に進む。
FIG. 4 is a diagram showing a flow when the subject is imaged in the first embodiment.
In step ST1, a contrast medium is injected into the subject, and a breath-hold scan BH is executed. In the breath holding scan BH, as shown in FIG. 5, the navigator sequence NS n and the imaging sequence IS n are executed alternately. When navigator sequence NS n is executed, navigator signals A n (n = 1 to z) are acquired from navigator region R nav, and when imaging sequence IS n is executed to reconstruct an image of imaging region R Imaging data S n (n = 1 to z) used in the above is collected. Signal strength data generating means 91a (see FIG. 1) is a one-dimensional Fourier transform of the navigator signals A n (1D-FT) and the
ステップST2では、検出手段91b(図1参照)が、息止めスキャンBHの間における肝臓のエッジのSI方向の位置P0を検出する。 At step ST2, the detection means 91b (see FIG. 1) detects the position P 0 of the SI direction of the edge of the liver during the breath-hold scan BH.
図6は、息止めスキャンBHの間における肝臓のエッジのSI方向の位置P0を検出するときの説明図である。 Figure 6 is an explanatory view when detecting the position P 0 of the SI direction of the edge of the liver during the breath-hold scan BH.
本形態では、肺側が低信号となり肝臓側が高信号となるようなナビゲータシーケンスを用いている。したがって、肺と肝臓との境界には信号強度の段差が現れるので、この段差を検出することにより、肝臓のエッジの位置P0を検出することができる。尚、本形態では、被検体は息を止めているので、スキャンの間、肝臓のエッジの位置は動かないと考えることができる。したがって、信号強度データD1〜Dzのうちのいずれかの信号強度データを用いて検出された肝臓のエッジの位置を、息止めスキャンの間における肝臓のエッジの位置P0とすることができる。ここでは、信号強度データD1を用いて検出された肝臓のエッジの位置を、息止めスキャンの間における肝臓のエッジの位置P0とする。 In this embodiment, a navigator sequence is used in which the lung side has a low signal and the liver side has a high signal. Accordingly, a step in signal intensity appears at the boundary between the lung and the liver. By detecting this step, the position P 0 of the edge of the liver can be detected. In this embodiment, since the subject holds his / her breath, it can be considered that the position of the edge of the liver does not move during the scan. Therefore, the position of the liver edge detected using any one of the signal intensity data D 1 to D z can be set as the liver edge position P 0 during the breath holding scan. . Here, the detected position of the liver edge with a signal intensity data D 1, and the position P 0 of the liver edge between the breath-hold scan.
肝臓のエッジの位置P0を検出した後、ステップST3に進む。
ステップST3では、検出手段91bが、信号強度データD1〜Dzごとに、肝臓のエッジの位置P0から肝臓側(I側)に距離Δdだけ離れた位置P0+Δdの信号強度N1〜Nzを検出する。図7に、位置P0+Δdの信号強度N1〜Nzを示す。信号強度N1〜Nzを検出することにより、肝臓における信号強度の時間変化を表すプロファイルFを得ることができる。図8に、そのプロファイルFを概略的に示す。肝臓に流入する造影剤の濃度が増加するに伴って、信号強度Nnの値は増加する。したがって、プロファイルFから、時間の経過とともに、肝臓に流入する造影剤の濃度が増加していることがわかる。尚、Δdの値が小さすぎると、肝臓と肺との境界に現れる段差の信号強度を検出してしまう可能性が高くなるので、肺側の領域の信号強度を検出してしまう恐れがある。一方、Δdの値が大きすぎると、腸の領域の信号強度を検出してしまう恐れがある。したがって、Δdは、肝臓の内側の信号強度を検出することができるように設定する必要がある。一例としては、Δd=5cmとすることにより、肝臓とは別の部位の信号強度を検出する危険を回避することができる。
プロファイルFを得た後、ステップST4に進む。
After detecting the position P 0 of the liver edge, the process proceeds to step ST3.
In step ST3, the detection means 91b detects the signal intensity N 1 to the signal intensity N 1 to the position P 0 + Δd separated from the liver edge position P 0 by the distance Δd from the liver edge position P 0 for each signal intensity data D 1 to D z. Nz is detected. FIG. 7 shows signal intensities N 1 to N z at the position P 0 + Δd. By detecting the signal intensities N 1 to N z , it is possible to obtain a profile F representing the time change of the signal intensity in the liver. FIG. 8 schematically shows the profile F thereof. As the concentration of the contrast agent flowing into the liver increases, the value of the signal intensity N n increases. Therefore, it can be seen from the profile F that the concentration of the contrast agent flowing into the liver increases with time. If the value of Δd is too small, there is a high possibility that the signal strength of the step appearing at the boundary between the liver and the lung will be detected, so there is a possibility that the signal strength of the lung side region will be detected. On the other hand, if the value of Δd is too large, the signal intensity of the intestinal region may be detected. Therefore, Δd needs to be set so that the signal intensity inside the liver can be detected. As an example, by setting Δd = 5 cm, it is possible to avoid the danger of detecting the signal intensity of a part other than the liver.
After obtaining the profile F, the process proceeds to step ST4.
ステップST4では、信号強度N1〜Nzに基づいて、イメージングシーケンスにより得られた撮影領域のイメージングデータS1〜Szを補正する。この補正は、以下の式を用いて行われる。
Sn,corr=Sn*Ni/Nn ・・・(1)
Sn:n番目のナビゲータ信号の直後に得られたイメージングデータ
Sn,corr:Snを補正した後のイメージングデータ
Ni:信号強度N1〜Nzの中から選択されたi番目の信号強度
Nn:n番目の信号強度
In step ST4, based on the
S n, corr = S n * N i / N n (1)
S n : Imaging data obtained immediately after the n-th navigator signal S n, corr : Imaging data after correcting S n N i : i-th signal selected from signal intensities N 1 to N z strength n n: n-th signal intensity
演算手段911c(図1参照)は、式(1)を用いて、補正後のイメージングデータSn,corrを計算する。式(1)では、イメージングデータSnをNi/Nn倍している。したがって、肝臓に流入する造影剤の濃度の違いにより生じるイメージングデータS1〜Szの信号強度の差を補正することができる。式(1)のNiは、信号強度N1〜Nzの中から選択されたi番目の信号強度であるが、Niは、例えば、i=1、すなわち、N1とすることができる。i=1の場合、式(1)は以下のように書き換えられる。
Sn,corr=Sn*N1/Nn ・・・(1)′
The
S n, corr = S n * N 1 / N n (1) ′
以下に、式(1)′を用いてイメージングデータS1〜Szを補正する方法について、図9を参照しながら説明する。 Hereinafter, a method of correcting the imaging data S 1 to S z using the expression (1) ′ will be described with reference to FIG.
図9は、式(1)′を用いてイメージングデータS1〜Szを補正するときの説明図である。 FIG. 9 is an explanatory diagram when the imaging data S 1 to S z are corrected using the equation (1) ′.
イメージングデータS1を補正する場合、式(1)′のnに、n=1を代入すればよい。
S1,corr=S1*N1/N1
=S1
したがって、補正後のイメージングデータS1,corrは、補正前のイメージングデータS1と同じである。
When correcting imaging data S 1, the n of the formula (1) 'may be substituted for n = 1.
S 1, corr = S 1 *
= S 1
Thus, imaging data S 1, corr corrected is the same as the imaging data S 1 before correction.
しかし、イメージングデータS2を補正する場合、式(1)′のnに、n=2が代入されるので、補正後のイメージングデータS2,corrは、以下の式で表される。
S2,corr=S2*N1/N2
However, when correcting the imaging data S 2, the n of the formula (1) ', since n = 2 is substituted, the imaging data S 2, corr corrected is expressed by the following equation.
S 2, corr = S 2 * N 1 / N 2
したがって、補正後のイメージングデータS2,corrは、補正前のイメージングデータS2のN1/N2倍になる。 Therefore, the corrected imaging data S 2 and corr are N 1 / N 2 times the uncorrected imaging data S 2 .
以下同様に、式(1)′のnの値を変更し、SnをN1/Nn倍することにより、補正後のイメージングデータSn,corrを得ることができる。イメージングデータを補正した後、ステップST5に進む。 Similarly, by changing the value of n in formula (1) ', the S n by multiplying N 1 / N n, imaging data S n after correction, can be obtained corr. After correcting the imaging data, the process proceeds to step ST5.
ステップST5では、補正後のイメージングデータS1,corr〜Sz,corrを用いて、画像再構成を行い、フローを終了する。 In step ST5, image reconstruction is performed using the corrected imaging data S1 , corr to Sz , corr , and the flow ends.
第1の形態では、ナビゲータシーケンスで得られた肝臓の位置における信号強度N1〜Nzに基づいて、イメージングシーケンスで取得されたイメージングデータS1〜Szを補正している。したがって、肝臓に流入する造影剤の濃度の違いにより生じるイメージングデータS1〜Szの信号強度の差を補正することができるので、リンギングやブラーなどのアーチファクトを低減することができる。 In the first embodiment, the imaging data S 1 to S z acquired in the imaging sequence are corrected based on the signal intensities N 1 to N z at the liver position obtained in the navigator sequence. Therefore, it is possible to correct the difference in signal intensity of the imaging data S 1 to S z caused by the difference in concentration of the contrast medium flowing into the liver, it is possible to reduce artifacts such as ringing and blur.
尚、上記の例では、信号強度N1〜Nzは、それぞれ信号強度データD1〜Dzから検出されている。しかし、信号強度N1〜Nzの中の一部の信号強度を、他の信号強度から推定してもよい。第1の形態では、信号強度N1〜Nzは、時間とともに増加する傾向があるので(図8参照)、信号強度N1とNzとが分かれば、信号強度が時間とともに増加するプロファイルFを推定することができる。したがって、信号強度N1とNzとに基づいてプロファイルFを推定する推定手段を設けておき、推定されたプロファイルFから、他の信号強度N2〜Nz−1を推定してもよい。この場合、ナビゲータシーケンスNS1およびNSzのみを実行するだけで信号強度N1〜Nzを求めることができるので、息止めスキャンBHの時間を更に短くすることができる。 In the above example, the signal intensities N 1 to N z are detected from the signal intensity data D 1 to D z , respectively. However, some signal strengths among the signal strengths N 1 to N z may be estimated from other signal strengths. In the first embodiment, since the signal strengths N 1 to N z tend to increase with time (see FIG. 8), if the signal strengths N 1 and N z are known, the profile F in which the signal strength increases with time. Can be estimated. Therefore, an estimation means for estimating the profile F based on the signal strengths N 1 and N z may be provided, and the other signal strengths N 2 to N z−1 may be estimated from the estimated profile F. In this case, since the signal intensities N 1 to N z can be obtained simply by executing only the navigator sequences NS 1 and NS z , the time of the breath holding scan BH can be further shortened.
第1の形態では、n番目の信号強度Nnを用いて、n番目のナビゲータ信号の直後に得られたn番目のイメージングデータSnを補正している。しかし、n番目の信号強度Nnを用いて、n番目のナビゲータ信号の直前に収集されたn−1番目のイメージングデータSn−1を補正してもよい。 In the first embodiment, by using the n-th signal intensity N n, it is corrected n-th imaging data S n obtained immediately after the n-th navigator signal. However, by using the n-th signal intensity N n, the n-th navigator imaging data S n-1 n-1 th collected immediately before the signal may be corrected.
尚、上記の説明では、ナビゲータシーケンスNSnとNSn+1との間に1つのイメージングシーケンスISnが実行される場合について述べられている。しかし、ナビゲータシーケンスNSnとNSn+1との間に複数のイメージングシーケンスを実行してもよい(図10参照)。 In the above description, the case where one imaging sequence IS n is executed between the navigator sequences NS n and NS n + 1 is described. However, a plurality of imaging sequences may be executed between the navigator sequences NS n and NS n + 1 (see FIG. 10).
図10は、ナビゲータシーケンスNSnとNSn+1との間に複数のイメージングシーケンスが実行される場合の一例を示す図である。 FIG. 10 is a diagram illustrating an example in which a plurality of imaging sequences are executed between the navigator sequences NS n and NS n + 1 .
図10では、ナビゲータシーケンスNSnとNSn+1との間に、複数のイメージングデータが収集される。このような場合は、ナビゲータシーケンスNSnとNSn+1との間に収集された複数のイメージングデータの各々をN1/Nn倍すればよい。例えば、ナビゲータシーケンスNS2とNS3との間では、イメージングデータS21〜S2qが収集されるので、これらのイメージングデータS21〜S2qの各々を、N1/N2倍すればよい。これにより、補正後のイメージングデータS21,corr〜S2q,corrを得ることができる。
In FIG. 10, a plurality of imaging data are collected between navigator sequences NS n and NS n + 1 . In such a case, each of the plurality of imaging data collected between the navigator sequences NS n and NS n + 1 may be multiplied by N 1 / N n . For example, since the imaging data S 21 to S 2q are collected between the navigator sequences NS 2 and NS 3 , each of these imaging data S 21 to S 2q may be multiplied by N 1 / N 2 . Thus,
尚、図10では、n番目の信号強度Nnを用いて、n番目のナビゲータ信号とn+1番目のナビゲータ信号との間に収集されたイメージングデータを補正している。しかし、n番目の信号強度Nnを用いて、n−1番目のナビゲータ信号とn番目のナビゲータ信号との間に収集されたイメージングデータを補正してもよい。 In FIG. 10, imaging data collected between the nth navigator signal and the (n + 1) th navigator signal is corrected using the nth signal intensity Nn. However, imaging data collected between the (n-1) th navigator signal and the nth navigator signal may be corrected using the nth signal intensity Nn.
また、第1の形態では、息止めスキャンBHを終了した後でイメージングデータを補正しているが、息止めスキャンBHを実行しながら、イメージングデータを補正してもよい。 In the first embodiment, the imaging data is corrected after the breath holding scan BH is completed. However, the imaging data may be corrected while the breath holding scan BH is executed.
(2)第2の形態
第1の形態では、ナビゲータ信号とイメージングデータとを別々のシーケンスで収集していたが、第2の形態では、ナビゲータ信号とイメージングデータとを同じシーケンスで収集する場合について説明する。
(2) Second form In the first form, the navigator signal and the imaging data are collected in separate sequences. In the second form, the navigator signal and the imaging data are collected in the same sequence. explain.
図11は、第2の形態のMR装置の概略図である。
第2の形態のMR装置200は、第1の形態のMR装置100の信号強度データ作成手段91aおよび検出手段91bの代わりに、信号強度算出手段91dを有している。尚、その他の構成は、第1の形態のMR装置100と同じである。
FIG. 11 is a schematic diagram of an MR apparatus according to the second embodiment.
The
図12は第2の形態で実行されるスキャンを示す図、図13は撮影領域を概略的に示す図である。
第2の形態では、息止めスキャンBHが実行される。第2の形態の息止めスキャンBHでは、シーケンスCn(n=1〜z)が実行される。シーケンスCnは、撮影領域Rの画像を再構成するときに使用されるイメージングデータSnを取得するためのシーケンスである。また、シーケンスCnには、勾配磁場GxおよびGyが印加される直前に待ち時間Twaitが設けられており、待ち時間Twaitのときにナビゲータ信号Anが取得される。ナビゲータ信号Anは、撮影領域Rのk空間の中心のデータを表している。待ち時間Twaitは、例えば20μmである。
FIG. 12 is a diagram showing a scan executed in the second form, and FIG. 13 is a diagram schematically showing an imaging region.
In the second form, a breath holding scan BH is executed. In the breath holding scan BH of the second form, the sequence C n (n = 1 to z) is executed. Sequence C n is a sequence for acquiring the imaging data S n to be used when reconstructing an image of the imaging region R. The sequence C n is provided with a waiting time T wait immediately before the gradient magnetic fields Gx and Gy are applied, and the navigator signal An is acquired at the waiting time T wait . Navigator signals A n denotes the data of the center of k-space imaging region R. The waiting time T wait is, for example, 20 μm.
したがって、シーケンスCnを実行するたびに、イメージングデータSnの他に、ナビゲータ信号Anを得ることができる。第2の形態では、シーケンスCnにより得られたナビゲータ信号Anを用いてイメージングデータSnを補正する。以下に、第2の形態において被検体を撮影するときのフローについて説明する。 Thus, each execution of the sequence C n, in addition to the imaging data S n, it is possible to obtain a navigator signal A n. In the second embodiment, to correct the imaging data S n by using the navigator signal A n obtained by the sequence C n. Hereinafter, a flow when imaging a subject in the second embodiment will be described.
図14は、第2の形態のフローを示す図である。
ステップST1では、被検体に造影剤を注入し、息止めスキャンBHを実行する。息止めスキャンBHでは、図15に示すように、シーケンスCnが実行される。シーケンスCnを実行することにより、ナビゲータ信号An(n=1〜z)とイメージングデータSn(n=1〜z)が収集される。息止めスキャンBHを実行した後、ステップST2に進む。
FIG. 14 is a diagram showing a flow of the second embodiment.
In step ST1, a contrast medium is injected into the subject, and a breath-hold scan BH is executed. In the breath holding scan BH, the sequence C n is executed as shown in FIG. By executing the sequence C n , the navigator signal A n (n = 1 to z) and the imaging data S n (n = 1 to z) are collected. After executing the breath holding scan BH, the process proceeds to step ST2.
ステップST2では、信号強度算出手段91d(図11参照)が、ナビゲータ信号A1〜Azごとに撮影領域Rの信号強度を算出する(図16参照)。
In step ST2, the signal
図16は、撮影領域Rの信号強度を算出するときの説明図である。
第2の形態では、信号強度算出手段91dは、撮影領域Rから得られたナビゲータ信号A1〜Azの積分値を計算する。これにより、撮影領域Rの信号強度N1〜Nzを算出することができる。信号強度N1〜Nzを算出した後、ステップST3に進む。
FIG. 16 is an explanatory diagram when the signal intensity of the imaging region R is calculated.
In the second mode, the signal
ステップST3では、信号強度N1〜Nzに基づいて、シーケンスC1〜Czにより得られた撮影領域RのイメージングデータS1〜Szを補正する。この補正は、以下の式を用いて行われる。
Sn,corr=Sn*Ni/Nn ・・・(2)
Sn:n番目のイメージングデータ
Sn,corr:Snを補正した後のイメージングデータ
Ni:信号強度N1〜Nzの中から選択されたi番目の信号強度
Nn:n番目の信号強度
At step ST3, the on the basis of the
S n, corr = S n * N i / N n (2)
S n : n-th imaging data S n, corr : imaging data after correcting S n N i : i-th signal intensity selected from signal intensities N 1 to N z N n : n-th signal Strength
演算手段91c(図11参照)は、式(2)を用いて、補正後のイメージングデータSn,corrを計算する。式(2)では、イメージングデータSnをNi/Nn倍している。したがって、肝臓に流入する造影剤の濃度の違いにより生じるイメージングデータS1〜Szの信号強度の差を補正することができる。式(2)のNiは、信号強度N1〜Nzの中から選択されたi番目の信号強度であるが、Niは、例えば、i=1、すなわち、N1とすることができる。i=1の場合、式(2)は以下のように書き換えられる。
Sn,corr=Sn*N1/Nn ・・・(2)′
The
S n, corr = S n * N 1 / N n (2) ′
以下に、式(2)′を用いてイメージングデータS1〜Szを補正する方法について、図17を参照しながら説明する。 Hereinafter, a method of correcting the imaging data S 1 to S z using the expression (2) ′ will be described with reference to FIG.
図17は、式(2)′を用いてイメージングデータS1〜Szを補正するときの説明図である。
イメージングデータS1を補正する場合、式(2)′のnに、n=1を代入すればよい。
S1,corr=S1*N1/N1
=S1
FIG. 17 is an explanatory diagram when the imaging data S 1 to S z are corrected using Expression (2) ′.
When correcting imaging data S 1, the n of the formula (2) ', may be substituted for n = 1.
S 1, corr = S 1 *
= S 1
したがって、補正後のイメージングデータS1,corrは、補正前のイメージングデータS1と同じである。 Thus, imaging data S 1, corr corrected is the same as the imaging data S 1 before correction.
しかし、イメージングデータS2を補正する場合、式(2)′のnに、n=2が代入されるので、補正後のイメージングデータS2,corrは、以下の式で表される。
S2,corr=S2*N1/N2
However, when correcting the imaging data S 2, the n of the formula (2) ', since n = 2 is substituted, the imaging data S 2, corr corrected is expressed by the following equation.
S 2, corr = S 2 * N 1 / N 2
したがって、補正後のイメージングデータS2,corrは、補正前のイメージングデータS2のN1/N2倍になる。 Therefore, the corrected imaging data S 2 and corr are N 1 / N 2 times the uncorrected imaging data S 2 .
以下同様に、式(2)′のnの値を変更し、SnをN1/Nn倍することにより、補正後のイメージングデータSn,corrを得ることができる。イメージングデータを補正した後、ステップST4に進む。 Similarly, by changing the value of n in formula (2) ', the S n by multiplying N 1 / N n, imaging data S n after correction, can be obtained corr. After correcting the imaging data, the process proceeds to step ST4.
ステップST4では、補正後のイメージングデータS1,corr〜Sz,corrを用いて、画像再構成を行い、フローを終了する。 In step ST4, image reconstruction is performed using the corrected imaging data S1 , corr to Sz , corr , and the flow ends.
第2の形態では、シーケンスCnで得られた肝臓の位置における信号強度N1〜Nzに基づいて、イメージングデータS1〜Szを補正している。したがって、肝臓に流入する造影剤の濃度の違いにより生じるイメージングデータS1〜Szの信号強度の差を補正することができるので、リンギングやブラーなどのアーチファクトを低減することができる。 In the second embodiment, the imaging data S 1 to S z are corrected based on the signal intensities N 1 to N z at the position of the liver obtained in the sequence C n . Therefore, it is possible to correct the difference in signal intensity of the imaging data S 1 to S z caused by the difference in concentration of the contrast medium flowing into the liver, it is possible to reduce artifacts such as ringing and blur.
また、第2の形態では、シーケンスCnを実行するときに励起される撮影領域Rから、ナビゲータ信号とイメージングデータとを取得することができる。したがって、ナビゲータ信号を取得するときに励起される領域と、イメージングデータを取得するときに励起される領域とが同じであるので、イメージングデータの補正の精度を向上させることができる。 In the second embodiment, the navigator signal and the imaging data can be acquired from the imaging region R excited when the sequence C n is executed. Therefore, since the region excited when acquiring the navigator signal is the same as the region excited when acquiring the imaging data, the accuracy of correction of the imaging data can be improved.
尚、第2の形態では、n番目のシーケンスCnにより収集された信号強度Nnを用いて、n番目のシーケンスCnにより収集されたイメージングデータSnを補正している。しかし、n番目のシーケンスCnにより収集された信号強度Nnを用いて、n−1番目のシーケンスCn−1により収集されたイメージングデータSn−1を補正してもよい。 In the second embodiment, by using the n-th sequence C n signal strength N n collected by, and correcting the imaging data S n collected by the n-th sequence C n. However, by using the n-th sequence C n signal strength N n collected by the imaging data S n-1 collected by the n-1 th sequence C n-1 may be corrected.
(3)第3の形態
第1および第2の形態では、被検体に息止めさせた状態でスキャンを実行する例について説明したが、第3の形態では、被検体に自由に呼吸をさせた状態でスキャンを実行する例について説明する。尚、MR装置のハードウェア構成は、第1の形態と同じである。
(3) Third Embodiment In the first and second embodiments, the example in which the scan is executed in a state where the subject holds his / her breath was explained, but in the third embodiment, the subject was allowed to breathe freely. An example in which scanning is executed in a state will be described. The hardware configuration of the MR apparatus is the same as that in the first embodiment.
図18は、第3の形態で実行されるスキャンの説明図である。
第3の形態では、プリスキャンAと、本スキャンBとが実行される。
FIG. 18 is an explanatory diagram of a scan executed in the third mode.
In the third mode, pre-scan A and main scan B are executed.
プリスキャンAは、呼吸信号を検出するためのスキャンである。プリスキャンAでは、ナビゲータシーケンスNS1〜NSxが実行される。ナビゲータシーケンスNS1〜NSxは、第1の形態のナビゲータシーケンスNSn(n=1〜z)と同様に、ナビゲータ領域Rnav内のSI方向の信号強度の変化を表す信号強度データを取得するためのシーケンスである。 Pre-scan A is a scan for detecting a respiratory signal. In the prescan A, navigator sequences NS 1 to NS x are executed. The navigator sequences NS 1 to NS x obtain signal strength data representing changes in signal strength in the SI direction within the navigator region R nav , as with the navigator sequence NS n (n = 1 to z) of the first embodiment. It is a sequence for.
本スキャンBは、肝臓を含む撮影領域Rの画像を再構成するときに使用されるイメージングデータを取得するためのスキャンである。本スキャンBでは、第1の形態の息止めスキャンBHと同様に、ナビゲータシーケンスNSn(n=1〜z)と、イメージングシーケンスISn(n=1〜z)とが交互に実行される。 The main scan B is a scan for acquiring imaging data used when reconstructing an image of the imaging region R including the liver. In the main scan B, the navigator sequence NS n (n = 1 to z) and the imaging sequence IS n (n = 1 to z) are alternately executed as in the breath holding scan BH of the first embodiment.
以下に、プリスキャンAおよび本スキャンBにより被検体を撮影するときのフローについて説明する。 Hereinafter, a flow when the subject is imaged by the pre-scan A and the main scan B will be described.
図19は、第3の形態において被検体を撮影するときのフローを示す図である。
ステップST1では、プリスキャンAが実行される(図20参照)。
FIG. 19 is a diagram illustrating a flow when the subject is imaged in the third embodiment.
In step ST1, prescan A is executed (see FIG. 20).
図20は、プリスキャンAの説明図である。
プリスキャンAでは、ナビゲータシーケンスNS1〜NSxが実行される。ナビゲータシーケンスNS1〜NSxを実行することにより、ナビゲータ領域Rnavからナビゲータ信号An(n=1〜x)が収集される。信号強度データ作成手段91a(図1参照)は、ナビゲータ信号Anを一次元フーリエ変換(1D−FT)し、ナビゲータ領域Rnav内のSI方向の信号強度の変化を表す信号強度データD1〜Dxを作成する。肺と肝臓との境界には、信号強度の段差が現れるので、この信号強度の段差を検出することにより、肝臓のSI方向のエッジの位置P1〜Pxを検出することができる。信号強度データD1〜Dxごとに肝臓のエッジ位置を検出することにより、肝臓のエッジ位置が時間とともにどのように変化するかを知ることができる。したがって、肝臓のエッジ位置の時間変化を表すエッジ位置プロファイルを得ることができる。図21に、エッジ位置プロファイルを概略的に示す。プリスキャンAを実行した後、ステップST2に進む。
FIG. 20 is an explanatory diagram of the pre-scan A.
In the prescan A, navigator sequences NS 1 to NS x are executed. By executing the navigator sequences NS 1 to NS x , navigator signals An (n = 1 to x) are collected from the navigator region R nav . Signal strength data generating means 91a (see FIG. 1) is a one-dimensional Fourier transform of the navigator signals A n (1D-FT) and the
ステップST2では、エッジ位置プロファイルに基づいて、被検体11の呼吸の吐き終わりの範囲AWを求める。図22に吐き終りの範囲AWの一例を示す。呼吸の吐き終りの範囲AWは、例えば、肝臓のエッジ位置の最小値から一定の範囲として求めることができる。呼吸の吐き終りの範囲AWを求めた後、ステップST3に進む。 In step ST2, based on the edge position profile, a breathing end range AW of the subject 11 is obtained. FIG. 22 shows an example of the range AW at the end of spitting. The range AW at the end of exhalation of breathing can be obtained as a fixed range from the minimum value of the edge position of the liver, for example. After obtaining the range AW at the end of breathing, the process proceeds to step ST3.
ステップST3では、本スキャンBが実行される。
図23は、本スキャンBの説明図である。
本スキャンBでは、ナビゲータシーケンスNSn(n=1〜z)と、イメージングシーケンスISn(n=1〜z)とが交互に実行される。
In step ST3, the main scan B is executed.
FIG. 23 is an explanatory diagram of the main scan B.
In the main scan B, the navigator sequence NS n (n = 1 to z) and the imaging sequence IS n (n = 1 to z) are executed alternately.
本スキャンBでは、ナビゲータシーケンスで検出された肝臓のエッジ位置が、呼吸の吐き終わりの範囲AWに含まれている場合、その直後のイメージングシーケンスで取得されたイメージングデータを、画像再構成のデータとして採用する。一方、ナビゲータシーケンスで検出された肝臓のエッジ位置が、呼吸の吐き終わりの範囲AWに含まれていない場合、その直後のイメージングシーケンスで取得されたイメージングデータは、画像再構成のデータとしては採用されない。例えば、図23では、ナビゲータシーケンスNS1により検出された肝臓のエッジ位置P1は、呼吸の吐き終わりの範囲AWから外れているので、その直後のイメージングシーケンスIS1で取得されたイメージングデータS1は、画像再構成のデータとしては採用されず、破棄される。 In the main scan B, when the edge position of the liver detected in the navigator sequence is included in the range AW of the end of breathing exhalation, the imaging data acquired in the immediately subsequent imaging sequence is used as image reconstruction data. adopt. On the other hand, when the edge position of the liver detected in the navigator sequence is not included in the end-of-breathing end range AW, the imaging data acquired in the immediately following imaging sequence is not adopted as image reconstruction data. . For example, in FIG. 23, since the liver edge position P 1 detected by the navigator sequence NS 1 is out of the range AW at the end of breathing exhalation, the imaging data S 1 acquired in the imaging sequence IS 1 immediately after that Are not adopted as data for image reconstruction and are discarded.
一方、ナビゲータシーケンスNS2により検出された肝臓のエッジ位置P2は、呼吸の吐き終わりの範囲AWに含まれているので、その直後のイメージングシーケンスIS2で取得されたイメージングデータS2は、画像再構成のデータとして採用される。同様に、ナビゲータシーケンスNS3により検出された肝臓のエッジ位置P3も、呼吸の吐き終わりの範囲AWに含まれているので、その直後のイメージングシーケンスIS3で取得されたイメージングデータS3も、画像再構成のデータとして採用される。
本スキャンBを実行した後、ステップST4に進む。
On the other hand, the edge position P 2 of the liver detected by the navigator sequence NS 2, so are included in the scope AW of end spitting respiratory, imaging data S 2 acquired by the imaging sequence IS 2 immediately thereafter, the image Adopted as reconstruction data. Similarly, since the liver edge position P 3 detected by the navigator sequence NS 3 is also included in the end-of-breathing end range AW, the imaging data S 3 acquired in the imaging sequence IS 3 immediately after that is also Used as data for image reconstruction.
After executing the main scan B, the process proceeds to step ST4.
ステップST4では、信号強度データD1〜Dzの中から、肝臓のエッジの位置が呼吸の吐き終りの範囲AWに含まれている信号強度データを特定する。信号強度データD1〜Dzにおける肝臓のエッジの位置P1〜Pzは既知であり、呼吸の吐き終りの範囲AWも既知であるので、肝臓のエッジの位置P1〜Pzが呼吸の吐き終りの範囲AWに含まれているか否かを判断することができる。したがって、信号強度データD1〜Dzの中から、肝臓のエッジの位置が呼吸の吐き終りの範囲AWに含まれている信号強度データを特定することができる。図23では、例えば、信号強度データD1は、肝臓のエッジの位置P1が呼吸の吐き終わりの範囲AWから外れているが、信号強度データD2、D3、Dzなどは、肝臓のエッジの位置が呼吸の吐き終わりの範囲AWに含まれている。信号強度データを特定した後、ステップST5に進む。
In step ST4, from the signal strength data D 1 to D z, the position of the liver edge is to identify the signal strength data contained in the range AW of spitting breath end.
ステップST5では、検出手段91b(図1参照)が、ステップST4で特定された信号強度データに対して、肝臓のエッジの位置から肝臓側(I側)に距離Δdだけ離れた位置における信号強度を検出する。図24に、検出された信号強度を示す。図24では、代表して、信号強度データD2、D3、Dzに対して検出された信号強度N2、N3、Nzが示されている。信号強度N2、N3、Nzは、それぞれ、肝臓のエッジの位置P2、P3、Pzから肝臓側(I側)に距離Δdだけ離れた位置における信号強度を表している。信号強度を検出した後、ステップST6に進む。 In step ST5, the detection means 91b (see FIG. 1) obtains the signal intensity at a position separated by a distance Δd from the position of the liver edge to the liver side (I side) with respect to the signal intensity data specified in step ST4. To detect. FIG. 24 shows the detected signal strength. FIG. 24 representatively shows signal intensities N 2 , N 3 , and N z detected for the signal intensity data D 2 , D 3 , and D z . The signal intensities N 2 , N 3 , and N z represent the signal intensities at positions that are separated from the liver edge positions P 2 , P 3 , and P z by a distance Δd from the liver side (I side), respectively. After detecting the signal strength, the process proceeds to step ST6.
ステップST6では、画像再構成用のデータとして採用されたイメージングデータを補正する。この補正は、以下の式を用いて行われる。
Sn,corr=Sn*N2/Nn ・・・(3)
Sn:画像再構成用のデータとして採用されたn番目のイメージングデータ
Sn,corr:Snを補正した後のイメージングデータ
N2:データD2の位置P2+Δdにおける信号強度
Nn:データDnの位置Pn+Δdにおける信号強度
In step ST6, the imaging data adopted as data for image reconstruction is corrected. This correction is performed using the following equation.
S n, corr = S n * N 2 / N n (3)
S n : n-th imaging data adopted as data for image reconstruction S n, corr : Imaging data after correcting S n N 2 : signal intensity at position P 2 + Δd of data D 2 N n : data Signal strength at position P n + Δd of D n
演算手段91c(図11参照)は、式(3)を用いて、補正後のイメージングデータSn,corrを計算する。以下に、式(3)を用いてイメージングデータを補正する方法について説明する。
The
図25は、式(3)を用いてイメージングデータを補正するときの説明図である。
イメージングデータS2を補正する場合、式(3)のnに、n=2を代入すればよい。
S2,corr=S2*N2/N2
=S2
FIG. 25 is an explanatory diagram when the imaging data is corrected using Expression (3).
When correcting imaging data S 2, the n of the formula (3) may be substituted for n = 2.
S 2, corr = S 2 *
= S 2
したがって、補正後のイメージングデータS2,corrは、補正前のイメージングデータS2と同じである。 Thus, imaging data S 2, corr corrected is the same as the imaging data S 2 before correction.
しかし、イメージングデータS3を補正する場合、式(3)のnに、n=3が代入されるので、補正後のイメージングデータS3,corrは、以下の式で表される。
S3,corr=S3*N2/N3
したがって、補正後のイメージングデータS3,corrは、補正前のイメージングデータS3のN2/N3倍になる。
However, when correcting the imaging data S 3, the n of the formula (3), since n = 3 is substituted, the imaging data S 3, corr corrected is expressed by the following equation.
S 3, corr = S 3 * N 2 / N 3
Thus, imaging data S 3, corr corrected will N 2 / N 3 times the uncorrected imaging data S 3.
以下同様に、式(3)のnの値を変更し、画像再構成用のデータとして採用されたイメージングデータSnをN2/Nn倍することにより、補正後のイメージングデータSn,corrを得ることができる。イメージングデータを補正した後、ステップST7に進む。 Similarly, by changing the value of n in formula (3), the imaging data S n, which was adopted as data for image reconstruction by multiplying N 2 / N n, imaging data S n corrected, corr Can be obtained. After correcting the imaging data, the process proceeds to step ST7.
ステップST7では、補正後のイメージングデータを用いて画像再構成を行い、フローを終了する。 In step ST7, image reconstruction is performed using the corrected imaging data, and the flow ends.
第3の形態では、画像再構成用のデータとして採用されたイメージングデータSnをN2/Nn倍している。したがって、肝臓に流入する造影剤の濃度の違いにより生じるイメージングデータSnの信号強度の差を補正することができるので、リンギングやブラーなどのアーチファクトを低減することができる。 In a third embodiment, and the imaging data S n, which was adopted as data for image reconstruction N 2 / N n multiplied. Therefore, it is possible to correct the difference in signal intensity of the imaging data S n caused by the difference in concentration of the contrast medium flowing into the liver, it is possible to reduce artifacts such as ringing and blur.
また、第3の形態では、肝臓のエッジの位置が呼吸の吐き終りの範囲AWに含まれているときに取得されたイメージングデータのみを用いて画像再構成が行われるので、自由呼吸下で被検体を撮影しても、ゴーストなどの体動アーチファクトを低減することができる。 In the third embodiment, image reconstruction is performed using only the imaging data acquired when the position of the edge of the liver is included in the end-of-breathing end range AW. Even when a specimen is photographed, body movement artifacts such as ghosts can be reduced.
尚、式(3)の信号強度N2の代わりに、別の信号強度Niを用いてイメージングデータを補正してもよい。 Instead of the signal intensity N 2 of formula (3), it may be corrected imaging data using another signal strength N i.
(4)第4の形態
第4の形態では、被検体に自由に呼吸をさせた状態でスキャンを実行する別の例について説明する。尚、MR装置のハードウェア構成は、第1の形態と同じである。
(4) Fourth Embodiment In the fourth embodiment, another example will be described in which a scan is executed in a state where the subject freely breathes. The hardware configuration of the MR apparatus is the same as that in the first embodiment.
図26は、第4の形態で実行されるスキャンの説明図である。
第4の形態では、プリスキャンAと、本スキャンBとが実行される。
FIG. 26 is an explanatory diagram of a scan executed in the fourth mode.
In the fourth mode, pre-scan A and main scan B are executed.
プリスキャンAでは、第3の形態と同様に、ナビゲータシーケンスNS1〜NSxが実行される。 In prescan A, navigator sequences NS 1 to NS x are executed as in the third embodiment.
本スキャンBでは、肝臓を含む撮影領域の画像を再構成するときに使用されるイメージングデータを収集するためのイメージングシーケンスIS1〜ISwが実行される。また、イメージングシーケンスIS1〜ISwの他に、ナビゲータシーケンスNSが実行される。ナビゲータシーケンスNSは、各イメージングシーケンスIS1〜ISwの直前に1回以上実行される。 In the main scan B, imaging sequences IS 1 to IS w for collecting imaging data used when reconstructing an image of an imaging region including the liver are executed. In addition to the imaging sequences IS 1 to IS w , a navigator sequence NS is executed. The navigator sequence NS is executed once or more immediately before each imaging sequence IS 1 to IS w .
以下に、プリスキャンAおよび本スキャンBにより被検体を撮影するときのフローについて説明する。 Hereinafter, a flow when the subject is imaged by the pre-scan A and the main scan B will be described.
図27は、第4の形態において被検体を撮影するときのフローを示す図である。
ステップST1では、プリスキャンAが実行される。プリスキャンAを実行することにより、第3の形態と同様に、肝臓のエッジ位置の時間変化を表すエッジ位置プロファイル(図21参照)を得ることができる。プリスキャンAを実行した後、ステップST2に進む。
FIG. 27 is a diagram showing a flow when the subject is imaged in the fourth embodiment.
In step ST1, prescan A is executed. By executing the pre-scan A, an edge position profile (see FIG. 21) representing a temporal change in the liver edge position can be obtained as in the third embodiment. After performing the pre-scan A, the process proceeds to step ST2.
ステップST2では、エッジ位置プロファイルに基づいて、トリガレベルTLが設定される。図28に、設定されたトリガレベルTLの一例を示す。トリガレベルTLは、本スキャンBを実行するときに使用されるものであり、具体的には、本スキャンBにおいてイメージングシーケンスIS1〜ISw(図26参照)を実行するときの肝臓のエッジの基準位置を表している。トリガレベルTLが本スキャンBを実行するときにどのように使用されるかについては、後述する。トリガレベルTLは、例えば、肝臓のエッジの位置の最大値と最小値の中間に設定することができる。トリガレベルTLを設定した後、ステップST3に進む。 In step ST2, the trigger level TL is set based on the edge position profile. FIG. 28 shows an example of the set trigger level TL. The trigger level TL is used when the main scan B is executed. Specifically, the trigger level TL is the level of the edge of the liver when the imaging sequences IS 1 to IS w (see FIG. 26) are executed in the main scan B. It represents the reference position. How the trigger level TL is used when executing the main scan B will be described later. The trigger level TL can be set, for example, between the maximum value and the minimum value of the position of the liver edge. After setting the trigger level TL, the process proceeds to step ST3.
ステップST3では、本スキャンBが実行される。
図29は、本スキャンBの説明図である。
本スキャンBでは、先ず、ナビゲータシーケンスNS11〜NS1Pが実行される。ナビゲータシーケンスNS11〜NS1Pにより収集されたナビゲータ信号A11〜A1Pに基づいてナビゲータ領域Rnavの信号強度データD11〜D1Pを求め、信号強度データD11〜D1Pごとに肝臓のエッジの位置P11〜P1Pを検出する。そして、肝臓のエッジの位置がトリガレベルTLの下側(I側)から上側(S側)に移動したときに、イメージングシーケンスIS1を実行する。イメージングシーケンスIS1を実行することにより、撮影領域RのイメージングデータS1が収集される。
In step ST3, the main scan B is executed.
FIG. 29 is an explanatory diagram of the main scan B.
In the main scan B, first, navigator sequences NS 11 to NS 1P are executed.
イメージングシーケンスIS1を実行した後、ナビゲータシーケンスNS21〜NS2qが実行される。そして、肝臓のエッジの位置がトリガレベルTLの下側(I側)から上側(S側)に移動したときに、次のイメージングシーケンスIS2を実行する。イメージングシーケンスIS2を実行することにより、撮影領域RのイメージングデータS2が収集される。
After executing the imaging sequence IS 1,
以下同様にナビゲータシーケンスとイメージングシーケンスを実行する。最後にナビゲータシーケンスNSz1〜NSzrを実行し、肝臓のエッジの位置がトリガレベルTLの下側(I側)から上側(S側)に移動したときに、イメージングシーケンスISwを実行する。イメージングシーケンスISwを実行することにより、撮影領域RのイメージングデータSwが収集され、本スキャンBが終了する。本スキャンBが終了した後、ステップST4に進む。 Similarly, the navigator sequence and the imaging sequence are executed. Finally running navigator sequence NS z1 ~NS zr, when the position of the liver edge is moved from the lower side of the trigger level TL (I side) to the upper (S side), it executes the imaging sequence IS w. By executing the imaging sequence IS w, is collected imaging data S w of the imaging region R, the scan B is completed. After the main scan B is completed, the process proceeds to step ST4.
ステップST4では、検出手段91b(図1参照)は、イメージングシーケンスIS1〜ISwの直前に取得された信号強度データD1p〜Dwrに対して、肝臓のエッジの位置から肝臓側(I側)に距離Δdだけ離れた位置における信号強度を検出する。図30に、検出された信号強度を示す。図30では、代表して、信号強度データD1p、D2q、Dwrに対して検出された信号強度N1、N2、Nwが示されている。信号強度N1、N2、Nwは、それぞれ、肝臓のエッジの位置P1p、P2q、Pwrから肝臓側(I側)に距離Δdだけ離れた位置における信号強度を表している。信号強度を検出した後、ステップST5に進む。 In step ST4, the detection means 91b (see FIG. 1) detects the signal intensity data D 1p to D wr acquired immediately before the imaging sequences IS 1 to IS w from the liver edge position to the liver side (I side). ) Is detected at a position separated by a distance Δd. FIG. 30 shows the detected signal strength. FIG. 30 representatively shows signal intensities N 1 , N 2 , and N w detected for the signal intensity data D 1p , D 2q , and D wr . The signal intensities N 1 , N 2 , and N w represent the signal intensities at positions separated from the liver edge positions P 1p , P 2q , and P wr by the distance Δd on the liver side (I side), respectively. After detecting the signal strength, the process proceeds to step ST5.
ステップST5では、イメージングデータS1〜Swを補正する。この補正は、以下の式を用いて行われる。
Sn,corr=Sn*N1/Nn ・・・(4)
Sn:n番目のイメージングデータ
Sn,corr:Snを補正した後のイメージングデータ
N1:1番目のイメージングデータS1の直前の信号強度データから検出された信号強度
Nn:n番目のイメージングデータSnの直前の信号強度データから検出された信号強度
In step ST5, to correct the
S n, corr = S n * N 1 / N n (4)
S n : n-th imaging data S n, corr : imaging data after correcting S n N 1 : signal intensity detected from signal intensity data immediately before the first imaging data S 1 N n : n-th signal intensity detected from the immediately preceding signal strength data of the imaging data S n
以下に、式(4)を用いてイメージングデータを補正する方法について説明する。
図31は、式(4)を用いてイメージングデータを補正するときの説明図である。
イメージングデータS1を補正する場合、式(4)のnに、n=1を代入すればよい。
S1,corr=S1*N1/N1
=S1
Hereinafter, a method for correcting imaging data using Expression (4) will be described.
FIG. 31 is an explanatory diagram when the imaging data is corrected using Expression (4).
When correcting imaging data S 1, the n of formula (4) may be substituted for n = 1.
S 1, corr = S 1 *
= S 1
したがって、補正後のイメージングデータS1,corrは、補正前のイメージングデータS1と同じである。 Thus, imaging data S 1, corr corrected is the same as the imaging data S 1 before correction.
しかし、イメージングデータS2を補正する場合、式(4)のnに、n=2が代入されるので、補正後のイメージングデータS2,corrは、以下の式で表される。
S2,corr=S2*N1/N2
However, when correcting the imaging data S 2, the n of formula (4), since n = 2 is substituted, the imaging data S 2, corr corrected is expressed by the following equation.
S 2, corr = S 2 * N 1 / N 2
したがって、補正後のイメージングデータS2,corrは、補正前のイメージングデータS2のN1/N2倍になる。 Therefore, the corrected imaging data S 2 and corr are N 1 / N 2 times the uncorrected imaging data S 2 .
以下同様に、式(4)のnの値を変更し、イメージングデータSnをN1/Nn倍することにより、補正後のイメージングデータSn,corrを得ることができる。イメージングデータを補正した後、ステップST6に進む。 Similarly, by changing the value of n in formula (4), the imaging data S n by multiplying N 1 / N n, imaging data S n after correction, can be obtained corr. After correcting the imaging data, the process proceeds to step ST6.
ステップST6では、補正後のイメージングデータを用いて画像再構成を行い、フローを終了する。 In step ST6, image reconstruction is performed using the corrected imaging data, and the flow ends.
第4の形態では、イメージングデータSnをN1/Nn倍しているので、肝臓に流入する造影剤の濃度の違いにより生じるイメージングデータSnの信号強度の差を補正することができるので、リンギングやブラーなどのアーチファクトを低減することができる。 In the fourth embodiment, since the imaging data S n are N 1 / N n times, it is possible to correct the difference in signal intensity of the imaging data S n caused by the difference in concentration of the contrast medium flowing into the liver Artifacts such as ringing and blur can be reduced.
また、第4の形態では、肝臓のエッジの位置がトリガレベルTLに到達したときにイメージングデータを取得するので、イメージングデータS1〜SWを取得するときの肝臓のエッジの位置のばらつきを小さくすることができる。したがって、ゴーストなどの体動アーチファクトを低減することができる。 In the fourth embodiment, since the position of the liver edge is to acquire imaging data when it reaches the trigger level TL, the variation in the position of the liver edge when acquiring the imaging data S 1 to S W smaller can do. Therefore, body motion artifacts such as ghosts can be reduced.
尚、式(4)の信号強度N1の代わりに、別の信号強度Niを用いてイメージングデータを補正してもよい。 Note that the imaging data may be corrected using another signal intensity N i instead of the signal intensity N 1 in Expression (4).
(5)第5の形態
第1〜第4の形態では、イメージングデータを補正することによりリンギングやブラーなどのアーチファクトを低減する例について説明したが、第5の形態では、RFパルスのフリップ角を補正することによりリンギングやブラーなどのアーチファクトを低減する例について説明する。
(5) Fifth Embodiment In the first to fourth embodiments, an example of reducing artifacts such as ringing and blur by correcting imaging data has been described. However, in the fifth embodiment, the flip angle of the RF pulse is set. An example of reducing artifacts such as ringing and blur by correcting will be described.
図32は、第5の形態のMR装置の概略図である。
第5の形態のMR装置300は、第1の形態のMR装置100の補正手段91(図1参照)の代わりに、RFパルスのフリップ角を補正する補正手段911を有している。尚、その他の構成は、第1の形態のMR装置100と同じである。制御部9は、補正手段911および再構成手段92を構成する一例であり、所定のプログラムを実行することにより、これらの手段として機能する。
FIG. 32 is a schematic view of an MR apparatus according to the fifth embodiment.
The
図33は、第5の形態で実行されるスキャンを示す図である。
第5の形態では、息止めスキャンBHが実行される。息止めスキャンBHでは、ナビゲータシーケンスNSn(n=1〜z)と、イメージングシーケンスISn(n=1〜z)とが交互に実行される。図33の下側に、イメージングシーケンスISnの一例が示されている。ここでは、イメージングシーケンスISnの一例として、SPGR(Spoiled Gradient Recalled)シーケンスが示されている。RFパルスXnのフリップ角αは、例えば90°である。
以下に、息止めスキャンBHを実行する方法について説明する。
FIG. 33 is a diagram illustrating a scan executed in the fifth mode.
In the fifth mode, a breath holding scan BH is executed. In the breath holding scan BH, the navigator sequence NS n (n = 1 to z) and the imaging sequence IS n (n = 1 to z) are alternately executed. An example of the imaging sequence IS n is shown on the lower side of FIG. Here, as an example of the imaging sequence IS n, SPGR (Spoiled Gradient Recalled ) sequence is shown. The flip angle α of the RF pulse Xn is 90 °, for example.
Hereinafter, a method for executing the breath holding scan BH will be described.
息止めスキャンBHでは、先ず、ナビゲータシーケンスNS1とイメージングシーケンスIS1とが実行される。図34に、ナビゲータシーケンスNS1により収集されたナビゲータ信号A1と、イメージングシーケンスIS1により収集されたイメージングデータS1とを示す。信号強度データ作成手段911a(図32参照)は、ナビゲータ信号A1を一次元フーリエ変換し、これにより、ナビゲータ領域Rnav内のSI方向の信号強度の変化を表す信号強度データD1が得られる。 In the breath holding scan BH, first, a navigator sequence NS 1 and an imaging sequence IS 1 are executed. FIG. 34 shows the navigator signal A 1 collected by the navigator sequence NS 1 and the imaging data S 1 collected by the imaging sequence IS 1 . Signal strength data generating means 911a (see FIG. 32), the navigator signal A 1 Fourier transform one dimension, thereby, the signal intensity data D 1 representing a change in the SI direction of the signal strength of the navigator region R nav is obtained .
検出手段911b(図32参照)は、信号強度データD1の信号強度の段差を検出し、肝臓のエッジのSI方向の位置P1を検出する。そして、肝臓のエッジの位置P1から肝臓側(I側)に距離Δdだけ離れた位置P1+Δdの信号強度N1を検出する。 Detection means 911b (see FIG. 32) detects the level difference of the signal intensity of the signal intensity data D 1, to detect the position P 1 of the SI direction liver edge. Then, the signal intensity N 1 at a position P 1 + Δd that is a distance Δd away from the liver edge position P 1 on the liver side (I side) is detected.
イメージングシーケンスIS1を実行した後、ナビゲータシーケンスNS2が実行される。図35に、ナビゲータシーケンスNS2により収集されたナビゲータ信号A2を示す。ナビゲータ信号A2は一次元フーリエ変換され、ナビゲータ領域Rnav内のSI方向の信号強度の変化を表す信号強度データD2が得られる。そして、信号強度データD2の位置P1+Δdにおける信号強度N2を検出する。 After executing the imaging sequence IS 1, navigator sequence NS 2 is executed. FIG. 35 shows the navigator signal A 2 collected by the navigator sequence NS 2 . The navigator signal A 2 is subjected to a one-dimensional Fourier transform to obtain signal intensity data D 2 representing a change in signal intensity in the SI direction within the navigator region R nav . Then, to detect the signal intensity N 2 at position P 1 + [Delta] d of the signal intensity data D 2.
次に、イメージングシーケンスIS2のRFパルスX2のフリップ角αを補正する。以下に、図36を参照しながら、フリップ角αの補正方法について説明する。フリップ角αは、以下の式を用いて補正される。
αn_corr=arcsin{(N1/Nn)*sinα} ・・・(5)
α:イメージングシーケンスISnの補正前のフリップ角
αn_corr:イメージングシーケンスISnの補正後のフリップ角
N1:信号強度データD1の位置P1+Δdにおける信号強度
Nn:信号強度データDnの位置P1+Δdにおける信号強度
Next, the flip angle α of the RF pulse X 2 of the imaging sequence IS 2 is corrected. Hereinafter, a method of correcting the flip angle α will be described with reference to FIG. The flip angle α is corrected using the following equation.
α n_corr = arcsin {(N 1 / N n ) * sin α} (5)
α: Flip angle before correction of imaging sequence IS n α n_corr : Flip angle after correction of imaging sequence IS n N 1 : Signal intensity at position P 1 + Δd of signal intensity data D 1 N n : Signal intensity data D n Signal strength at position P 1 + Δd
演算手段911c(図32参照)は、式(5)を用いて、補正後のフリップ角αn,corrを計算する。ここでは、n=2、つまり、イメージングシーケンスIS2のRFパルスX2のフリップ角αを補正するので、式(5)のnに、n=2を代入すればよい。
α2_corr=arcsin{(N1/N2)*sinα}
The computing means 911c (see FIG. 32) calculates the corrected flip angle α n, corr using equation (5). Here, since n = 2, that is, the flip angle α of the RF pulse X 2 of the imaging sequence IS 2 is corrected, n = 2 may be substituted for n in the equation (5).
α 2 — corr = arcsin {(N 1 / N 2 ) * sin α}
したがって、補正後のフリップ角α2_corrを求めることができる。補正後のフリップ角α2_corrを求めたら、イメージングシーケンスIS2のRFパルスX2のフリップ角αをα2_corrに補正し、撮影領域Rのスキャンを実行する。したがって、フリップ角α=α2_corrのときのイメージングデータS2を得ることができる。
Therefore, the corrected
以下同様に、ナビゲータシーケンスNSnを実行するたびに、式(5)を用いて補正後のフリップ角αn_corrを計算し、イメージングシーケンスISnを実行する。したがって、フリップ角α=αn_corrのときのイメージングデータSnを得ることができる。図37に、息止めスキャンBHにより得られたイメージングデータS1〜Szを示す。イメージングデータS1〜Szを収集した後、これらのイメージングデータS1〜Szを用いて画像再構成を行う。これにより、撮影領域Rの画像データを得ることができる。
Similarly, every time the navigator sequence NS n is executed, the corrected flip angle α n_corr is calculated using the equation (5), and the imaging sequence IS n is executed. Therefore, it is possible to obtain the imaging data S n when the flip angle α = α n_corr. FIG. 37 shows the imaging data S 1 to S z obtained by the breath-hold scan BH. After collecting
SPGRシーケンスでは、(N1/Nn)*sinαのarcsineを求めることにより、造影剤の濃度の違いが補正されたイメージングデータS1〜Szを収集することができるので、リンギングやブラーなどのアーチファクトを低減することができる。 In the SPGR sequence, by obtaining the arcsine of (N 1 / N n ) * sin α, it is possible to collect imaging data S 1 to S z in which the difference in contrast agent concentration is corrected. Artifacts can be reduced.
(6)第6の形態
第5の形態では、ナビゲータ信号とイメージングデータとを別々のシーケンスで収集する例について示されているが、第6の形態では、ナビゲータ信号とイメージングデータとを一つのシーケンスCnで収集する例について説明する。
(6) Sixth Embodiment In the fifth embodiment, an example is shown in which navigator signals and imaging data are collected in separate sequences, but in the sixth embodiment, navigator signals and imaging data are combined into one sequence. An example of collecting with C n will be described.
図38は、第6の形態のMR装置を示す概略図である。
第6の形態のMR装置400は、第5の形態のMR装置300(図32参照)の信号強度データ作成手段911aおよび検出手段911bの代わりに、信号強度算出手段911dを備えているが、その他の構成は、第5の形態と同じである。
FIG. 38 is a schematic diagram showing an MR apparatus according to the sixth embodiment.
The
図39は、第6の形態で実行されるスキャンを示す図である。
第6の形態では、息止めスキャンBHが実行される。息止めスキャンBHでは、第2の形態と同様に、シーケンスCn(n=1〜z)が実行される。
FIG. 39 is a diagram illustrating a scan executed in the sixth mode.
In the sixth embodiment, a breath holding scan BH is executed. In the breath holding scan BH, the sequence C n (n = 1 to z) is executed as in the second embodiment.
先ず、シーケンスC1が実行される。図40に、シーケンスC1により収集されたナビゲータ信号A1およびイメージングデータS1を示す。シーケンスC1を実行することによりナビゲータ信号A1およびイメージングデータS1が収集される。信号強度算出手段911d(図38参照)は、撮影領域R(図13参照)から得られたナビゲータ信号A1の積分値を計算する。これにより、撮影領域Rの信号強度N1を算出することができる。
First, the sequence C 1 is executed. FIG. 40 shows the navigator signal A 1 and the imaging data S 1 collected by the sequence C 1 . Navigator signals A 1 and imaging data S 1 is collected by performing a sequence C 1. Signal
シーケンスC1を実行した後、次のシーケンスC2が実行される。図41に、シーケンスC2により収集されたナビゲータ信号A2およびイメージングデータS2を示す。信号強度算出手段911dは、撮影領域Rから得られたナビゲータ信号A2の積分値を計算する。これにより、撮影領域Rの信号強度N2を算出することができる。
After executing the sequence C 1, the following sequence C 2 is performed. Figure 41 shows the navigator signal A 2 and the imaging data S 2 collected by the sequence C 2. Signal
シーケンスC2を実行した後、シーケンスC3が実行される。ただし、シーケンスC3を実行する前に、シーケンスC3のRFパルスX3のフリップ角αが補正される。以下に、シーケンスC3のRFパルスX3のフリップ角αの補正方法について述べる。 After executing the sequence C 2, a sequence C 3 is performed. However, before executing the sequence C 3, flip angle α of RF pulses X 3 of the sequence C 3 is corrected. Hereinafter, a method for correcting the flip angle α of the RF pulse X 3 of the sequence C 3 will be described.
図42は、シーケンスC3のRFパルスX3のフリップ角αの補正方法の説明図である。
シーケンスC3のRFパルスX3のフリップ角αは、以下の式を用いて補正される。
α3_corr=arcsin{(N1/N2)*sinα} ・・・(6)
α:シーケンスC3のRFパルスX3の補正前のフリップ角
α3_corr:シーケンスC3のRFパルスX3の補正後のフリップ角
N1:ナビゲータ信号A1により得られた撮影領域Rの信号強度
N2:ナビゲータ信号A2により得られた撮影領域Rの信号強度
Figure 42 is an explanatory view of a method of correcting a flip angle α of RF pulses X 3 of the sequence C 3.
The flip angle α of the RF pulse X 3 of the sequence C 3 is corrected using the following equation.
α 3 — corr = arcsin {(N 1 / N 2 ) * sin α} (6)
alpha: uncorrected flip angle alpha 3_Corr of RF pulses X 3 of the sequence C 3: Sequence C 3 of the RF pulses X 3 correction flip angle after N 1: signal intensity of the imaging region R obtained by the navigator signal A 1 N 2 : Signal intensity of the imaging region R obtained by the navigator signal A 2
演算手段911c(図38参照)は、式(6)を用いて、補正後のフリップ角α3_corrを計算する。式(6)の信号強度N1およびN2は既に算出されており、フリップ角αは事前に決まっている値である。したがって、式(6)にN1、N2、およびα3を代入することにより、補正後のフリップ角α3_corrを求めることができる。補正後のフリップ角α3_corrを求めたら、シーケンスC3のRFパルスX3のフリップ角αをα3_corrに補正し、シーケンスC3を実行する。したがって、ナビゲータ信号A3とイメージングデータS3とを得ることができる。信号強度算出手段911dは、撮影領域Rから得られたナビゲータ信号A3の積分値を計算し、撮影領域Rの信号強度N3を算出する。
The computing means 911c (see FIG. 38) calculates the corrected
シーケンスC3を実行した後、シーケンスC4が実行される。ただし、シーケンスC4を実行する前に、シーケンスC4のRFパルスX4のフリップ角αが補正される。以下に、シーケンスC4のRFパルスX4のフリップ角αの補正方法について述べる。 After executing the sequence C 3, sequence C 4 is performed. However, before executing the sequence C 4, flip angle α of RF pulses X 4 of the sequence C 4 is corrected. The following describes a correction method of the flip angle α of RF pulses X 4 of the sequence C 4.
シーケンスC3のフリップ角は、αからα3_corrに補正されているので、シーケンスC3の信号強度N3は、フリップ角α3_corrにより得られた値である。一方、シーケンスC1の信号強度N1はフリップ角αにより得られた値である。したがって、信号強度N1およびN3は異なるフリップ角により得られている。そこで、信号強度N3を、信号強度N1と同様に、フリップ角αにより得られた値に補正する。信号強度N3は、以下の式(7)を用いて補正することができる。
Nn_corr=Nn(sinα/sinαn_corr) ・・・(7)
α:シーケンスCnのRFパルスXnの補正前のフリップ角
αn_corr:シーケンスCnのRFパルスXnの補正後のフリップ角
Nn:ナビゲータ信号Anにより得られた撮影領域Rの信号強度
Nn_corr:Nnの補正後の値
Flip angle sequence C 3, since it is corrected to alpha 3_Corr from alpha, signal strength N 3 of the sequence C 3 is a value obtained by a flip angle α 3_corr. On the other hand, the signal intensity N 1 sequence C 1 is a value obtained by flip angle alpha. Therefore, the signal strengths N 1 and N 3 are obtained with different flip angles. Therefore, the signal intensity N 3 is corrected to a value obtained by the flip angle α, similarly to the signal intensity N 1 . Signal strength N 3 can be corrected using the following equation (7).
N n — corr = N n (sin α / sin α n — corr ) (7)
alpha: Sequence C n of RF pulses X n pre-correction flip angle alpha N_corr: sequence C n RF pulses X n correction flip angle after N n of: signal intensity of the imaging region R obtained by the navigator signal A n N n — corr : Value after correction of N n
ここでは、n=3、つまり、信号強度N3を補正するので、式(7)のnに、n=3を代入すればよい。
N3_corr=Nn(sinα/sinα3_corr)
Here, since n = 3, that is, the signal intensity N 3 is corrected, n = 3 may be substituted for n in the equation (7).
N 3 — corr = N n (sin α /
したがって、フリップ角α3_corrにより得られた信号強度N3を、フリップ角αにより得られた信号強度N3_corrに補正することができる。図43に、補正後の信号強度N3_corrが算出された様子を示す。補正後の信号強度N3_corrを求めた後、次のシーケンスC4のフリップ角を補正する。
Therefore, the signal intensity N 3 obtained by the flip angle α 3_corr can be corrected to the signal intensity N 3_corr obtained by the flip angle α. FIG. 43 shows a state in which the corrected
図44は、シーケンスC4のフリップ角を補正するときの説明図である。
シーケンスC4のRFパルスX4のフリップ角αは、以下の式を用いて補正される。
αn+1_corr=arcsin{(N1/Nn_corr)*sinα} ・・・(8)
Figure 44 is an explanatory view of correcting the flip angle of the sequence C 4.
The flip angle α of the RF pulse X 4 of the sequence C 4 is corrected using the following equation.
α n + 1_corr = arcsin {(N 1 / N n_corr ) * sin α} (8)
演算手段911cは、式(8)を用いて、補正後のフリップ角αn_corrを計算する。ここでは、シーケンスC4のRFパルスX4のフリップ角αを補正するので、式(8)のnに、n=3を代入すればよい。
α4_corr=arcsin{(N1/N3_corr)*sinα}
The calculating means 911c calculates the corrected flip angle α n_corr using the equation (8). Here, since the flip angle α of the RF pulse X 4 of the sequence C 4 is corrected, n = 3 may be substituted for n in the equation (8).
α 4 — corr = arcsin {(N 1 / N 3 —corr ) * sin α}
したがって、補正後のフリップ角α4_corrを求めることができる。補正後のフリップ角α4_corrを求めたら、シーケンスC4のRFパルスX4のフリップ角αをα4_corrに補正し、シーケンスC4を実行する。したがって、ナビゲータ信号A4とイメージングデータS4とを得ることができる。信号強度算出手段911dは、撮影領域Rから得られたナビゲータ信号A4の積分値を計算し、撮影領域Rの信号強度N4を算出する。
Therefore, the corrected
以下同様に、シーケンスCnを実行するたびに、式(7)を用いて信号強度Nnを補正してから、式(8)を用いて補正後のフリップ角αn+1_corrを求め、シーケンスを実行する。最後のシーケンスCzが終了したら、画像再構成を行う。これにより、撮影領域Rの画像データを得ることができる。 Similarly, each time the sequence C n is executed, the signal intensity N n is corrected using the equation (7), and then the corrected flip angle α n + 1_corr is obtained using the equation (8), and the sequence is executed. To do. When the last sequence C z is completed, image reconstruction is performed. Thereby, the image data of the imaging region R can be obtained.
第6の形態では、一つのシーケンスCnで、ナビゲータ信号AnとイメージングデータSnとの両方を収集している。この場合は、式(7)を用いて信号強度を補正した後で、フリップ角を補正することにより、息止めスキャンBHにより得られるイメージングデータの信号強度の差を小さくすることができるので、リンギングやブラーなどのアーチファクトを低減することができる。 In the sixth embodiment, a single sequence C n, are collecting both the navigator signal A n and imaging data S n. In this case, by correcting the flip angle after correcting the signal intensity using the equation (7), the difference in the signal intensity of the imaging data obtained by the breath-hold scan BH can be reduced. And artifacts such as blur can be reduced.
第5および第6の形態では息止めスキャンの場合について説明したが、第3および第4の形態と同様に、被検体に自由に呼吸をさせた状態でスキャンを実行し、フリップ角を補正してもよい。 In the fifth and sixth embodiments, the case of breath-hold scanning has been described. However, as in the third and fourth embodiments, the scan is executed with the subject freely breathing to correct the flip angle. May be.
尚、第1〜第6の形態では、造影剤を用いてスキャンを行う場合について説明されているが、本発明は、イメージングデータの信号強度が時間とともに変化する場合であれば、造影剤を用いないスキャンを実行する場合にも適用できる。 In the first to sixth embodiments, the case where scanning is performed using a contrast medium has been described. However, the present invention uses a contrast medium if the signal intensity of imaging data changes with time. It can also be applied when not scanning.
更に、第1〜第6の形態では、体動する部位の画像再構成を行う場合について説明されているが、本発明は、体動しない部位の画像再構成を行う場合にも適用することができる。 Furthermore, in the first to sixth embodiments, the case of performing image reconstruction of a part that moves is explained. However, the present invention can also be applied to the case of performing image reconstruction of a part that does not move. it can.
2 マグネット
3 テーブル
3a クレードル
4 受信コイル
5 造影剤注入装置
6 送信器
7 勾配磁場電源
8 受信器
9 制御部
10 操作部
11 表示部
12 被検体
21 ボア
22 超伝導コイル
23 勾配コイル
24 RFコイル
91a、911a 信号強度データ作成手段
91b、911b 検出手段
91c、911c 演算手段
92 再構成手段
100、200、300、400 MR装置
2
Claims (21)
前記ナビゲータ信号と前記イメージングデータとを収集するためのスキャンを実行するスキャン手段と、
前記ナビゲータ信号に基づいて、前記第1の領域の所定の方向における信号強度の変化を表す信号強度データを作成し、前記信号強度データの前記所定の方向における第1の位置の信号強度に基づいて前記イメージングデータの信号強度を補正する補正手段と、
信号強度が補正されたイメージングデータに基づいて画像再構成を行う再構成手段と、
を有する磁気共鳴装置。 The navigator signal is collected from the first area including the predetermined part of the subject, and the imaging data used when reconstructing the image of the predetermined part from the second area including the predetermined part is collected. A magnetic resonance apparatus,
Scanning means for executing a scan for collecting the navigator signal and the imaging data;
Based on the navigator signal, signal strength data representing a change in signal strength in a predetermined direction of the first region is created, and based on the signal strength of the first position in the predetermined direction of the signal strength data. Correction means for correcting the signal intensity of the imaging data;
Reconstruction means for performing image reconstruction based on imaging data whose signal intensity is corrected;
A magnetic resonance apparatus.
前記信号強度データを作成する信号強度データ作成手段と、
前記信号強度データの前記第1の位置の信号強度を検出する検出手段とを有し、
前記信号強度に基づいて前記イメージングデータの信号強度を補正する、請求項1に記載の磁気共鳴装置。 The correction means includes
And the signal intensity data creating means for creating the signal strength data,
Detecting means for detecting the signal intensity at the first position of the signal intensity data;
The magnetic resonance apparatus according to claim 1, wherein the signal intensity of the imaging data is corrected based on the signal intensity.
前記信号強度データに基づいて前記所定の部位のエッジの位置を検出し、前記所定の部位のエッジの位置から所定の距離離れた位置の信号強度を、前記第1の位置の信号強度として検出する、請求項2に記載の磁気共鳴装置。 The detection means includes
The position of the edge of the predetermined part is detected based on the signal intensity data, and the signal intensity at a position away from the edge position of the predetermined part by a predetermined distance is detected as the signal intensity of the first position. The magnetic resonance apparatus according to claim 2.
前記ナビゲータ信号と前記イメージングデータとを交互に収集し、
前記補正手段は、
i番目のナビゲータ信号により得られた信号強度データの前記第1の位置の信号強度と、n番目のナビゲータ信号により得られた信号強度データの前記第1の位置の信号強度とに基づいて、前記n番目のナビゲータ信号の直前又は直後に収集されたイメージングデータを補正する、請求項2又は3に記載の磁気共鳴装置。 The scanning means includes
Alternately collecting the navigator signal and the imaging data;
The correction means includes
Based on the signal strength at the first position of the signal strength data obtained by the i-th navigator signal and the signal strength at the first position of the signal strength data obtained by the n-th navigator signal, The magnetic resonance apparatus according to claim 2, wherein the imaging data collected immediately before or after the nth navigator signal is corrected.
前記ナビゲータ信号と複数のイメージングデータとを交互に収集し、
前記補正手段は、
i番目のナビゲータ信号により得られた信号強度データの前記第1の位置と、n番目のナビゲータ信号により得られた信号強度データの前記第1の位置の信号強度とに基づいて、前記n番目のナビゲータ信号とn+1番目のナビゲータ信号との間に収集されたイメージングデータ、又は前記n番目のナビゲータ信号とn−1番目のナビゲータ信号との間に収集されたイメージングデータを補正する、請求項2又は3に記載の磁気共鳴装置。 The scanning means includes
Alternately collecting the navigator signal and a plurality of imaging data;
The correction means includes
Based on the first position of the signal strength data obtained from the i-th navigator signal and the signal strength at the first position of the signal strength data obtained from the n-th navigator signal, the n th The imaging data collected between the navigator signal and the (n + 1) th navigator signal or the imaging data collected between the nth navigator signal and the (n-1) th navigator signal is corrected. 3. The magnetic resonance apparatus according to 3.
前記ナビゲータ信号に基づいて前記第1の領域の信号強度を算出する信号強度算出手段を有し、
前記信号強度に基づいて前記イメージングデータの信号強度を補正する、請求項1に記載の磁気共鳴装置。 The correction means includes
Signal intensity calculating means for calculating the signal intensity of the first region based on the navigator signal;
The magnetic resonance apparatus according to claim 1, wherein the signal intensity of the imaging data is corrected based on the signal intensity.
前記ナビゲータ信号と前記イメージングデータとを交互に収集し、
前記補正手段は、
i番目のナビゲータ信号により得られた信号強度データの前記第1の位置の信号強度と、n番目のナビゲータ信号により得られた信号強度データの前記第1の位置の信号強度とに基づいて、前記n番目のナビゲータ信号の直前又は直後に収集されたイメージングデータを補正する、請求項8に記載の磁気共鳴装置。 The scanning means includes
Alternately collecting the navigator signal and the imaging data;
The correction means includes
Based on the signal strength at the first position of the signal strength data obtained by the i-th navigator signal and the signal strength at the first position of the signal strength data obtained by the n-th navigator signal, The magnetic resonance apparatus according to claim 8, wherein the imaging data collected immediately before or after the n-th navigator signal is corrected.
前記ナビゲータ信号と複数のイメージングデータとを交互に実行し、
前記補正手段は、
i番目のナビゲータ信号により得られた信号強度データの前記第1の位置の信号強度と、n番目のナビゲータ信号により得られた信号強度データの前記第1の位置の信号強度とに基づいて、前記n番目のナビゲータ信号とn+1番目のナビゲータ信号との間に収集されたイメージングデータ、又は前記n番目のナビゲータ信号とn−1番目のナビゲータ信号との間に収集されたイメージングデータを補正する、請求項8に記載の磁気共鳴装置。 The scanning means includes
Alternately executing the navigator signal and a plurality of imaging data;
The correction means includes
Based on the signal strength at the first position of the signal strength data obtained by the i-th navigator signal and the signal strength at the first position of the signal strength data obtained by the n-th navigator signal, correcting the imaging data collected between the nth navigator signal and the n + 1th navigator signal, or the imaging data collected between the nth navigator signal and the (n-1) th navigator signal. Item 9. The magnetic resonance apparatus according to Item 8.
前記スキャン手段は、前記被検体が息止めをしている間に、前記ナビゲータ信号と前記イメージングデータとを収集するためのスキャンを実行する、請求項1〜12のうちのいずれか一項に記載の磁気共鳴装置。 The predetermined part is a part that moves,
The scanning unit performs scanning for collecting the navigator signal and the imaging data while the subject is holding a breath. Magnetic resonance device.
前記スキャン手段は、呼吸信号を取得するためのプリスキャンと、前記ナビゲータ信号と前記イメージングデータとを収集するための本スキャンとを実行し、
前記ナビゲータ信号は、前記体動する部位のエッジの位置と、前記体動する部位の信号強度とを求めるために使用される、請求項1〜3のうちのいずれか一項に記載の磁気共鳴装置。 The predetermined part is a part that moves,
The scanning means performs a pre-scan for acquiring a respiratory signal and a main scan for collecting the navigator signal and the imaging data,
The magnetic resonance according to any one of claims 1 to 3, wherein the navigator signal is used to obtain a position of an edge of the body moving part and a signal intensity of the body moving part. apparatus.
本スキャンでは、前記ナビゲータ信号と前記イメージングデータとを交互に収集し、前記ナビゲータ信号により求められた前記体動する部位のエッジの位置が、前記呼吸の範囲に含まれている場合、前記ナビゲータ信号の直前又は直後に収集されたイメージングデータを、画像再構成用のデータとして採用する、請求項14に記載の磁気共鳴装置。 Based on the respiration signal, set a respiration range used when determining whether to employ the imaging data as data for image reconstruction,
In the main scan, when the navigator signal and the imaging data are alternately collected, and the position of the edge of the body moving part determined by the navigator signal is included in the respiration range, the navigator signal The magnetic resonance apparatus according to claim 14, wherein imaging data collected immediately before or immediately after is adopted as data for image reconstruction.
前記本スキャンでは、前記体動する部位のエッジの位置が前記基準位置に到達したときに、前記イメージングデータを収集する、請求項14に記載の磁気共鳴装置。 Set the reference position of the edge of the body moving part based on the respiratory signal,
The magnetic resonance apparatus according to claim 14, wherein in the main scan, the imaging data is collected when a position of an edge of the body moving part reaches the reference position.
前記ナビゲータ信号に基づいて、前記第1の領域の所定の方向における信号強度の変化を表す信号強度データを作成し、前記信号強度データの前記所定の方向における第1の位置の信号強度に基づいて前記イメージングデータの信号強度を補正する補正処理と、
信号強度が補正されたイメージングデータに基づいて画像再構成を行う再構成処理と、
を計算機に実行させるためのプログラム。 The navigator signal is collected from the first area including the predetermined part of the subject, and the imaging data used when reconstructing the image of the predetermined part from the second area including the predetermined part is collected. A magnetic resonance apparatus program for executing a scan for performing
Based on the navigator signal, signal strength data representing a change in signal strength in a predetermined direction of the first region is created, and based on the signal strength of the first position in the predetermined direction of the signal strength data. Correction processing for correcting the signal intensity of the imaging data;
Reconstruction processing to perform image reconstruction based on imaging data whose signal intensity is corrected;
A program to make a computer execute.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2012166187A JP6121664B2 (en) | 2012-07-26 | 2012-07-26 | Magnetic resonance apparatus and program |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2012166187A JP6121664B2 (en) | 2012-07-26 | 2012-07-26 | Magnetic resonance apparatus and program |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2014023703A JP2014023703A (en) | 2014-02-06 |
| JP6121664B2 true JP6121664B2 (en) | 2017-04-26 |
Family
ID=50197948
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2012166187A Active JP6121664B2 (en) | 2012-07-26 | 2012-07-26 | Magnetic resonance apparatus and program |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6121664B2 (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6411843B2 (en) * | 2014-09-30 | 2018-10-24 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Magnetic resonance equipment |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH09262220A (en) * | 1996-03-29 | 1997-10-07 | Hitachi Medical Corp | Magnetic resonance imaging equipment |
| US6424153B1 (en) * | 1999-11-23 | 2002-07-23 | Koninklijke Philips Electronics, N.V. | On-the-fly removal of data inconsistency with k-space oversampling and demodulation in MRI acquisitions |
| US7467006B2 (en) * | 2003-10-07 | 2008-12-16 | Hitachi Medical Corporation | Magnetic resonance imaging system and contrast-enhanced angiography |
| JP4350679B2 (en) * | 2005-05-26 | 2009-10-21 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Magnetic resonance imaging apparatus and image photographing apparatus |
| JP2010075634A (en) * | 2008-09-29 | 2010-04-08 | Toshiba Corp | Magnetic resonance imaging system |
| JP2010075573A (en) * | 2008-09-29 | 2010-04-08 | Ge Medical Systems Global Technology Co Llc | Magnetic resonance imaging instrument |
-
2012
- 2012-07-26 JP JP2012166187A patent/JP6121664B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2014023703A (en) | 2014-02-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ludwig et al. | Pilot tone–based motion correction for prospective respiratory compensated cardiac cine MRI | |
| CN101991417B (en) | Magnetic resonance imaging apparatus and controlling method thereof | |
| US8126237B2 (en) | Magnetic resonance imaging and correcting device | |
| CN101455565B (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
| US7432710B2 (en) | Apparatus and method for reducing image artifact | |
| CN106233154A (en) | Use prepulsing and the nuclear magnetic resonance with motion correction of omniselector | |
| JP6408594B2 (en) | Time-efficient 4D magnetic resonance imaging | |
| CN102078196A (en) | Method to detect breathing movement of an examination subject corresponding to signal data by magnetic resonance | |
| WO2006064582A1 (en) | Magnetic resonance imaging system | |
| JP6549612B2 (en) | Improved multiphase dynamic contrast magnetic resonance imaging method | |
| JPH05154130A (en) | Body motion artifact eliminating method | |
| WO2011034004A1 (en) | Magnetic resonance imaging device and method for applying a gradient magnetic field | |
| JPWO2004080301A1 (en) | Magnetic resonance imaging system | |
| US20180368721A1 (en) | Medical imaging device and magnetic resonance imaging device, and control method therefor | |
| US20150346307A1 (en) | Method and magnetic resonance apparatus for acquiring magnetic resonance data with a prospective motion correction | |
| JP5023127B2 (en) | Magnetic resonance imaging system | |
| JP4133348B2 (en) | Inspection device using nuclear magnetic resonance | |
| CN108663646A (en) | 2D airmanships | |
| JP5304987B2 (en) | Magnetic resonance imaging system | |
| JP6121664B2 (en) | Magnetic resonance apparatus and program | |
| JP6356412B2 (en) | Magnetic resonance apparatus and program | |
| JP2641486B2 (en) | NMR imaging equipment | |
| JP2008086343A (en) | Magnetic resonance imaging apparatus | |
| JP6411843B2 (en) | Magnetic resonance equipment | |
| JP5627903B2 (en) | Magnetic resonance imaging system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20150525 |
|
| RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20150602 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20150701 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20160519 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20160607 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20160906 |
|
| 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: 20170228 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20170330 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6121664 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 |