JP5230967B2 - Nuclear magnetic resonance imaging system - Google Patents
Nuclear magnetic resonance imaging system Download PDFInfo
- Publication number
- JP5230967B2 JP5230967B2 JP2007147865A JP2007147865A JP5230967B2 JP 5230967 B2 JP5230967 B2 JP 5230967B2 JP 2007147865 A JP2007147865 A JP 2007147865A JP 2007147865 A JP2007147865 A JP 2007147865A JP 5230967 B2 JP5230967 B2 JP 5230967B2
- Authority
- JP
- Japan
- Prior art keywords
- magnetic resonance
- flip angle
- resonance imaging
- imaging apparatus
- orthogonal
- 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
本発明は、核磁気共鳴イメージング装置(以下、MRI装置)において、アーチファクトおよび比吸収率を低減する技術に関する。特に、非直交系の走査における同技術に関する。 The present invention relates to a technique for reducing artifacts and specific absorptance in a nuclear magnetic resonance imaging apparatus (hereinafter referred to as an MRI apparatus). In particular, it relates to the same technique in non-orthogonal scanning.
MRI装置は、被検体、特に人体の組織を構成する原子核スピンが発生するNMR信号(エコー信号)を計測し、その頭部、腹部、四肢等の形態や機能を2次元的に或いは3次元的に画像化する装置である。撮影においては、エコー信号には、傾斜磁場によって異なる位相エンコード、周波数エンコードが付与される。計測されたエコー信号は、2次元又は3次元フーリエ変換されることにより画像に再構成される。 An MRI apparatus measures an NMR signal (echo signal) generated by a nuclear spin that constitutes a subject, particularly a human tissue, and forms the shape and function of the head, abdomen, limbs, etc. in a two-dimensional or three-dimensional manner. It is a device that images. In imaging, the echo signal is given different phase encoding and frequency encoding depending on the gradient magnetic field. The measured echo signal is reconstructed into an image by two-dimensional or three-dimensional Fourier transform.
MRI装置を用いた撮影では、スピンの励起に用いるRFパルスのエネルギーが人体に吸収される割合である比吸収率(Specific Absorption Rate;SAR)が高くなると、人体に対して発熱などの影響を与える。SARと撮影パラメータとの間には以下の(式1)の関係がある。
一般に、SARを低減するためには、TRを長くする(Dを小さくする)かフリップ角を小さくすればよい。しかし、TRを長くすることにより撮影時間の延長やコントラストの変化、フリップ角を低くすることにより信号ノイズ比(Signal to Noise Ratio; SNR)の低下、などの影響がある。 In general, in order to reduce the SAR, it is only necessary to increase TR (decrease D) or decrease the flip angle. However, there are effects such as an increase in shooting time and a change in contrast by increasing TR, and a decrease in signal-to-noise ratio (SNR) by decreasing a flip angle.
そこで、このような問題に対し、2次元または3次元k空間の位相エンコード量に応じてフリップ角を変化させ、RFパワーを低減することでSARを低減する技術が提案されている(非特許文献1、特許文献1)。非特許文献1には、高速スピンエコーシーケンスを構成するリフォーカスパルスのフリップ角(フロップ角)を、特許文献1には、SSFPシーケンスの励起パルスのフリップ角を、それぞれ位相エンコード量に応じて変更し、撮影全体のSARを低減する技術が開示されている。ここでは、計測空間(いわゆるk空間と呼ばれる空間)の低空間周波領域に配置されるエコー信号は、再構成画像のSNRに強く寄与するため、SNRの高い、高フリップ角の信号を低空間周波領域に配置し、再構成画像のSNRの低下を防いでいる。
Therefore, a technique for reducing the SAR by changing the flip angle and reducing the RF power according to the phase encoding amount in the two-dimensional or three-dimensional k-space has been proposed (Non-Patent Document). 1, Patent Document 1). Non-Patent
前述したように、高磁場のMRI装置を用いる撮影ではk空間位置によってフリップ角を可変にし、RFパワーを低減することでSARを低減させる。上記文献に開示の方法では、SNRの低下を防ぐため、SNRへの寄与が小さいk空間の高空間周波数領域に配置されるエコー信号を低いフリップ角のRFパルスで計測するよう、位相エンコード量の絶対値(3次元撮影の場合は、位相エンコード量とスライスエンコード量の合成量)に応じてフリップ角を制御する。これは、k空間をkx軸またはky軸に沿ってサンプリングする直交座標系の走査(直交系計測)では有効である。 As described above, in imaging using a high magnetic field MRI apparatus, the SAR is reduced by making the flip angle variable according to the k-space position and reducing the RF power. In the method disclosed in the above document, in order to prevent a decrease in SNR, the phase encoding amount is set so that an echo signal arranged in a high spatial frequency region of k-space that has a small contribution to SNR is measured with a low flip angle RF pulse. The flip angle is controlled in accordance with the absolute value (in the case of three-dimensional imaging, the combined amount of the phase encoding amount and the slice encoding amount). This is effective in orthogonal coordinate system scanning (orthogonal measurement) in which the k space is sampled along the kx axis or the ky axis.
直交系計測に対し、原点を中心に放射線状やラセン状にk空間を走査する非直交系の走査(非直交系計測)がある。非直交系計測では、1つのエコー信号でもky軸(位相エンコード方向)について異なる値をとるため、上述の方法でSNRの低下を防ぐことができない。 For orthogonal system measurement, there is non-orthogonal system scanning (non-orthogonal system measurement) in which the k-space is scanned radially or helically around the origin. In non-orthogonal measurement, even one echo signal takes different values for the ky axis (phase encoding direction), and thus the SNR cannot be prevented by the above-described method.
本発明は、上記事情に鑑みてなされたもので、非直交系計測において、SARを低減させつつ、良好なSNRを持たせる技術を提供することを目的とする。 The present invention has been made in view of the above circumstances, and an object of the present invention is to provide a technique for providing good SNR while reducing SAR in non-orthogonal measurement.
本発明は、SAR低減のため、非直交系計測内でフリップ角を変化させる。取得したエコー信号の各サンプリングデータの直交系座標(直交系グリッド)への再配置(グリッディング)において、k空間の低空間周波数領域の格子点にSNRの良好なサンプリングデータが高い寄与率で再配置されるようフリップ角に応じた重み係数を用いる。 The present invention changes the flip angle within non-orthogonal measurement to reduce SAR. In the rearrangement (gridding) of each sampling data of the acquired echo signal to the orthogonal system coordinates (orthogonal system grid), sampling data with good SNR is re-established at a high contribution rate to the lattice points in the low spatial frequency region of k-space. A weighting coefficient corresponding to the flip angle is used so as to be arranged.
具体的には、高周波磁場パルスおよび傾斜磁場パルスを制御し、非直交系計測空間データとなるエコー信号を計測する撮影手段と、前記非直交系計測空間データを直交系計測空間データに再配置するグリッディング手段と、前記直交系計測空間データを再構成して画像化する画像形成手段とを備える磁気共鳴イメージング装置であって、前記撮影手段は、前記非直交系計測空間データの収集に用いる高周波磁場パルスに複数種のフリップ角を設定するフリップ角設定手段を備え、前記グリッディング手段は、フリップ角の異なる非直交系計測空間データを用いて前記再配置を行うことを特徴とする磁気共鳴イメージング装置を提供する。
Specifically, imaging means for controlling the high-frequency magnetic field pulse and the gradient magnetic field pulse to measure an echo signal as non-orthogonal measurement space data, and rearranging the non-orthogonal measurement space data in the orthogonal measurement space data A magnetic resonance imaging apparatus comprising a gridding means and an image forming means for reconstructing and imaging the orthogonal measurement space data, wherein the imaging means is a high frequency used for collecting the non-orthogonal measurement space data Magnetic resonance imaging, comprising: flip angle setting means for setting a plurality of types of flip angles in a magnetic field pulse, wherein the gridding means performs the rearrangement using non-orthogonal measurement space data having different flip angles. Providing equipment.
本発明によれば、アーチファクトの少ない画像を再構成可能な非直交系計測において、SARを低減させつつ、良好なSNRを持たせることができる。 According to the present invention, in non-orthogonal measurement capable of reconstructing an image with few artifacts, it is possible to provide a favorable SNR while reducing the SAR.
<<第一の実施形態>>
以下、添付図面に従って本発明の実施形態について説明する。なお、本施形態を説明するための全図において、同一機能を有するものは同一符号を付け、その繰り返しの説明は省略する。最初に、第一の実施形態のMRI装置の全体概要を図1に基づいて説明する。図1は、本実施形態のMRI装置100の全体構成を示すブロック図である。MRI装置100は、NMR現象を利用して被検体の断層画像を得るもので、図1に示すように、静磁場発生系2と、傾斜磁場発生系3と、送信系5と、受信系6と、信号処理系7と、シーケンサ4と、中央処理装置(CPU)8とを備える。
<< First Embodiment >>
Hereinafter, embodiments of the present invention will be described with reference to the accompanying drawings. Note that components having the same function are denoted by the same reference symbols throughout the drawings for describing the embodiments, and the repetitive description thereof is omitted. First, an overall outline of the MRI apparatus of the first embodiment will be described with reference to FIG. FIG. 1 is a block diagram showing the overall configuration of the MRI apparatus 100 of the present embodiment. The MRI apparatus 100 obtains a tomographic image of a subject using an NMR phenomenon. As shown in FIG. 1, the MRI apparatus 100 includes a static magnetic
静磁場発生系2は、垂直磁場方式であれば、被検体1の周りの空間にその体軸と直交する方向に、水平磁場方式であれば、体軸方向に均一な静磁場を発生させるもので、被検体1の周りに永久磁石方式、常電導方式あるいは超電導方式の静磁場発生源が配置されている。
The static magnetic
傾斜磁場発生系3は、MRI装置の座標系(静止座標系)であるX、Y、Zの3軸方向に巻かれた傾斜磁場コイル9と、それぞれの傾斜磁場コイルを駆動する傾斜磁場電源10とを備え、後述のシ−ケンサ4からの命令に従ってそれぞれのコイルの傾斜磁場電源10を駆動することにより、X、Y、Zの3軸方向に傾斜磁場Gx、Gy、Gzを印加する。撮影時には、スライス面(撮影断面)に直交する方向にスライス方向傾斜磁場パルス(Gs)を印加して被検体1に対するスライス面を設定し、そのスライス面に直交して且つ互いに直交する残りの2つの方向に位相エンコード方向傾斜磁場パルス(Gp)と周波数エンコード方向傾斜磁場パルス(Gf)とを印加して、エコー信号にそれぞれの方向の位置情報をエンコードする。
The gradient magnetic
シーケンサ4は、高周波磁場パルス(以下、RFパルスという)と傾斜磁場パルスとを所定のパルスシーケンスに従って繰り返し印加するよう各部を制御する。CPU8の制御に従って動作し、被検体1の断層画像のデータ収集に必要な種々の命令を送信系5、傾斜磁場発生系3、および受信系6に送る。
The
送信系5は、被検体1の生体組織を構成する原子の原子核スピンに核磁気共鳴を起こさせるために、被検体1にRFパルスを照射するもので、高周波発振器11と変調器12と高周波増幅器13と送信側の高周波コイル(送信コイル)14aとを備える。送信系5は、高周波発振器11から出力された高周波パルスをシーケンサ4からの指令によるタイミングで変調器12によって振幅変調し、振幅変調された高周波パルスを高周波増幅器13で増幅した後に被検体1に近接して配置された送信コイル14aに供給することにより、RFパルスを被検体1に照射する。
The
受信系6は、被検体1の生体組織を構成する原子核スピンの核磁気共鳴により放出されるNMR信号(エコー信号)を検出するもので、受信側の高周波コイル(受信コイル)14bと信号増幅器15と直交位相検波器16と、A/D変換器17とを備える。送信コイル14aから照射された電磁波によって誘起された被検体1の応答のエコー信号が被検体1に近接して配置された受信コイル14bで検出され、信号増幅器15で増幅された後、シーケンサ4からの指令によるタイミングで直交位相検波器16により直交する二系統の信号に分割され、それぞれがA/D変換器17でディジタル量に変換されて、信号処理系7に送られる。
The receiving system 6 detects an NMR signal (echo signal) emitted by nuclear magnetic resonance of nuclear spins constituting the biological tissue of the
信号処理系7は、各種データ処理と処理結果の表示及び保存等とを行うもので、光ディスク19、磁気ディスク18等の外部記憶装置と、ROM21、RAM22等のメモリと、CRT等からなるディスプレイ20とを備える。受信系6から受け取ったデータに対してCPU8が信号処理、画像再構成等を実行した結果である被検体1の断層画像をディスプレイ20に表示すると共に、外部記憶装置の磁気ディスク18等に記録する。
The signal processing system 7 performs various data processing and display and storage of processing results. The signal processing system 7 includes an external storage device such as an
操作部25は、MRI装置100の各種制御情報や信号処理系7で行う処理の制御情報の入力インタフェースで、トラックボール又はマウス23、及び、キーボード24を備える。操作部25はディスプレイ20に近接して配置され、操作者がディスプレイ20を見ながら操作部25を通してインタラクティブにMRI装置100の各種処理を制御する。
The
なお、送信コイル14aと傾斜磁場コイル9は、被検体1が挿入される静磁場発生系2の静磁場空間内に、垂直磁場方式であれば被検体1に対向して、水平磁場方式であれば被検体1を取り囲むようにして設置される。また、受信コイル14bは、被検体1に対向して、或いは取り囲むように設置されている。現在、MRI装置の撮像対象核種であって臨床で普及しているものは、被検体の主たる構成物質である水素原子核(プロトン)である。プロトン密度の空間分布や、励起状態の緩和時間の空間分布に関する情報を画像化することで、人体頭部、腹部、四肢等の形態または、機能を2次元もしくは3次元的に撮像する。
Note that the
次に、一般的な非直交系計測におけるデータの再構成手順を説明する。ここでは、非直交系計測の中のラディアルサンプリング計測と呼ばれる、計測空間の略一点を回転中心として回転角を変えながら放射状にサンプリングを行い、1枚の画像再構成に必要なデータを得る計測を例にあげて説明する。また、ここでは、1枚の画像を構成するために必要な位相エンコードデータの取得を複数のセグメントに分割して行うものとする。以下、このような計測を非直交系セグメント計測と呼ぶ。 Next, a data reconstruction procedure in general non-orthogonal measurement will be described. Here, the radial sampling measurement in non-orthogonal system measurement is used to obtain data necessary for image reconstruction by sampling radially while changing the rotation angle around one point in the measurement space. An example will be described. Further, here, it is assumed that the acquisition of the phase encoding data necessary for constituting one image is performed by dividing it into a plurality of segments. Hereinafter, such measurement is referred to as non-orthogonal segment measurement.
図2は、非直交系セグメント計測におけるデータの再構成の手順を説明するための図である。ここでは、説明を簡単にするために、セグメントの数を3とする。本図に示すように、非直交系セグメント計測においては、2軸方向、例えば、位相エンコード方向と周波数エンコード方向との傾斜磁場を組み合わせるとともに、両者の大きさの比率を変化させて計測を繰り返す。各セグメントにおいて取得したエコー信号は、k空間原点を中心としてGxとGyとの比率で決まる傾斜を持つ走査線上にサンプリングされる(311、312、313)。セグメント単位で取得したデータ311、312、313をひとつの計測空間内に配置する(ステップ301)。そして、計測空間に配置されたデータに対してグリッディングと呼ばれる非直交系サンプリングデータを直交系サンプリングデータに変換する補間処理(以後、グリッディング処理と呼ぶ)を実施し、データを非直交座標系から直交座標系(直交系グリッド)に変換してk空間に再配置し(ステップ302)、再配置後のデータに対してフーリエ変換を行い(ステップ303)、画像化する(ステップ304)。
FIG. 2 is a diagram for explaining a data reconstruction procedure in non-orthogonal segment measurement. Here, in order to simplify the description, the number of segments is assumed to be 3. As shown in this figure, in non-orthogonal segment measurement, measurement is repeated while combining gradient magnetic fields in two axial directions, for example, a phase encoding direction and a frequency encoding direction, and changing the ratio of the magnitudes of both. The echo signal acquired in each segment is sampled on a scanning line having an inclination determined by the ratio of Gx and Gy with the k-space origin as the center (311, 312, 313).
グリッディング処理では、直交座標系の1つの格子点について、その点から一定距離にある非直交座標のデータを用いて当該格子点のデータを作成する。このとき、各サンプリングデータについて、当該データの直交系グリッドからの距離rの関数を重み係数Wとして作用させる。 In the gridding process, the data of the grid point is created for one grid point in the Cartesian coordinate system using the data of the non-orthogonal coordinates at a certain distance from the grid point. At this time, for each sampling data, a function of the distance r from the orthogonal grid of the data is applied as a weighting factor W.
次に、本実施形態の、非直交系セグメント計測においてフリップ角を変化させて行う計測(可変フリップ角非直交系セグメント計測と呼ぶ。)について説明する。図3は、可変フリップ角非直交系セグメント計測において、セグメント単位でフリップ角FAを変化させる計測を説明するための図である。本図において、(a)は、フリップ角FAと時間tとの関係を、(b)は、各フリップ角で取得するエコー信号と時間tとの関係を、(c)は、取得したエコー信号を計測空間に配置したイメージを、それぞれ示す。以下、本実施形態では、説明を簡単にするために、セグメントを3つとする。 Next, measurement (referred to as variable flip angle non-orthogonal segment measurement) performed by changing the flip angle in non-orthogonal segment measurement according to the present embodiment will be described. FIG. 3 is a diagram for explaining measurement for changing the flip angle FA in segment units in variable flip angle non-orthogonal segment measurement. In this figure, (a) shows the relationship between the flip angle FA and time t, (b) shows the relationship between the echo signal acquired at each flip angle and time t, and (c) shows the acquired echo signal. Each image is shown in the measurement space. Hereinafter, in this embodiment, in order to simplify the description, there are three segments.
本実施形態では、撮影全体のSARを低減するため、セグメント毎にフリップ角を変化させて撮像を行う。ここでは、セグメント211におけるフリップ角をFA1、セグメント212におけるフリップ角をFA2、セグメント213におけるフリップ角をFA3とすると、FA1>FA2>FA3となるよう変化させる。本実施形態では、各セグメントにおいてフリップ角FAを変化させているため、それぞれのセグメント211、212、213で取得したエコー信号(計測空間データ)221、222、223は、異なるSNRを持つ。異なるSNRを持つエコー信号221、222、223から得られたサンプリングデータは、図3(c)231、232、233に示すようにそれぞれ計測空間に配置される。なお、フリップ角の変化は上述のものに限られない。ただし、磁化の運動を安定に保つために、徐々に高く、もしくは低くしていくのが望ましい。
In the present embodiment, in order to reduce the SAR of the entire shooting, imaging is performed by changing the flip angle for each segment. Here, assuming that the flip angle in the
次に、本実施形態の可変フリップ角非直交系セグメント計測におけるデータの再構成手順を説明する。なお、以下の処理は、信号処理系7のCPU8が予めROM21またはRAM22に記憶されたプログラム、または、操作部25を介して入力されたユーザからの指示に従って、ROM21またはRAM22に記憶されたデータ、または、操作部を介してユーザから入力されたデータを用いて実行する。図4は、可変フリップ角非直交系セグメント計測におけるデータの再構成の手順を説明するための図である。データ再構成の手順は、上述の一般的な非直交系計測における手順と基本的に同様である。
Next, a data reconstruction procedure in variable flip angle non-orthogonal segment measurement of the present embodiment will be described. In the following processing, the CPU 8 of the signal processing system 7 stores the data stored in the
本実施形態では、セグメント毎にフリップ角を変化させているため、各セグメントで得られるエコー信号の信号強度は異なる。異なる信号強度のエコー信号から得られたサンプリングデータ1111、1112、1113を、ひとつの計測空間内に配置する(ステップ1101)。そして、グリッディング処理を行ってk空間に再配置し(ステップ1102)、再配置後のデータに対してフーリエ変換を行い(ステップ1103)、画像化する(ステップ1104)。
In this embodiment, since the flip angle is changed for each segment, the signal strength of the echo signal obtained in each segment is different.
上述のように本実施形態ではセグメント毎にフリップ角を変化させているため、全セグメントの信号を配置した計測空間には、フリップ角の異なるRFパルスによるエコー信号から得られたデータが混在する。従って、SNRに差があるエコー信号が計測空間の同じ空間周波数領域に配置される。この計測空間を従来の方法で再構成すると、再構成画像のSNRに大きく寄与する低空間周波数領域にも低いフリップ角のRFパルスによるエコー信号も配置され、画像のSNRが低下するとともに、信号強度差により画像にアーチファクトが生じる。 As described above, since the flip angle is changed for each segment in the present embodiment, data obtained from echo signals by RF pulses having different flip angles are mixed in a measurement space in which signals of all segments are arranged. Therefore, echo signals having a difference in SNR are arranged in the same spatial frequency region of the measurement space. When this measurement space is reconstructed by the conventional method, an echo signal by an RF pulse with a low flip angle is also arranged in a low spatial frequency region that greatly contributes to the SNR of the reconstructed image, and the SNR of the image is lowered and the signal intensity is reduced. The difference causes artifacts in the image.
そこで、本実施形態では、グリッディングする際に、再構成画像にアーチファクトの発生を低減させるように、距離rに加え、信号のフリップ角FAに応じた重みをデータに付けて再配置する。各データに作用させる重みを重み関数Wと呼ぶ。本実施形態では、重み関数Wは、各サンプリングデータ(計測点)の直交系グリッドからの距離rの関数w(r)と、フリップ角FAの関数w(FA)からなる。 Therefore, in this embodiment, when gridding is performed, weights corresponding to the signal flip angle FA are added to the data and rearranged so as to reduce the occurrence of artifacts in the reconstructed image. The weight applied to each data is called a weight function W. In the present embodiment, the weighting function W includes a function w (r) of the distance r from the orthogonal grid of each sampling data (measurement point) and a function w (FA) of the flip angle FA.
距離rの関数w(r)は、直交系グリッドからの距離が大きくなるにつれて値が小さくなるものを採用する。例えば、以下の(式2)を用いる。
図5に本実施形態の重み係数Wに用いる各関数を示す。図5(a)は、(式2)で示される関数w(r)のグラフ501である。このような関数W(r)を用いると、グリッディング時に直交系グリッドに近い計測点ほど重視される。なお、k空間原点から直交系グリッドまでの距離Rが大きくなると計測点の密度が低くなり、半径rの円内の計測点数が減少する。従って、原点からの距離Rに応じて重み係数W(r)を変化させてもよい。
FIG. 5 shows each function used for the weighting coefficient W of this embodiment. FIG. 5A is a
また、フリップ角FAの関数w(FA)は、高いフリップ角のRFパルスによるエコー信号に対しては大きく、低いフリップ角のそれに対しては小さく作用するものを採用する。これは、フリップ角の高いRFパルスによるエコー信号からの寄与を高め、画像SNRを向上させるためである。例えば、関数w(FA)として、例えば、以下の(式3)を用いる。
図5(b)は、(式3)で示されるフリップ角FAの関数w(FA)のグラフ502である。ここで、wminはフリップ角FAが0のときの重み係数を示す。wminはFAmax時のw(FA)の値(ここでは1)より小さい値をとり、0としても良い。本図に示すw(FA)は、SNRの高い高フリップ角のRFパルスによるエコー信号に対して大きく、低フリップ角のエコー信号に関しては小さい。従って、フリップ角が大きい(高い)サンプリングデータほど重視される。
FIG. 5B is a
次に、本実施形態の重み係数Wを用いたグリッディング処理について説明する。図6は、本実施形態のグリッディング処理の模式図である。図6(a)は、非直交系セグメント計測における計測空間のデータ配置の一部を表した図である。黒丸(401)は非直交系セグメント計測で取得したサンプリングデータ(計測点群)、白丸(402)は再配置する直交系グリッドである。また、R404は、原点からの距離である。本図においても、3つの異なるフリップ角(FA1>FA2>FA3)を用いて、エコー信号を取得したものとする。 Next, the gridding process using the weighting coefficient W of this embodiment will be described. FIG. 6 is a schematic diagram of the gridding process of the present embodiment. FIG. 6A is a diagram showing a part of the data arrangement in the measurement space in the non-orthogonal segment measurement. Black circles (401) are sampling data (measurement point group) acquired by non-orthogonal segment measurement, and white circles (402) are orthogonal grids to be rearranged. R404 is the distance from the origin. Also in this figure, it is assumed that the echo signal is acquired using three different flip angles (FA1> FA2> FA3).
本実施形態では、グリッディング処理時に用いる重み係数W=WFAを、直交系グリッドからの距離r(403)の関数w(r)とフリップ角FAの関数w(FA)との関数W(w(r)、w(FA))=WFA(r、FA)(405)として定義する。図6(b)に示すように、グリッディング処理前の各計測点(406)に、距離rおよびフリップ角FAにより変化する重み係数WFA(407)を作用させ、直交系グリッド408に再配置する。
In the present embodiment, the weighting factor W = W FA used in the gridding process is a function W (w) of the function w (r) of the distance r (403) from the orthogonal grid and the function w (FA) of the flip angle FA. (R), w (FA)) = W FA (r, FA) (405). As shown in FIG. 6B, the weighting factor W FA (407) that changes depending on the distance r and the flip angle FA is applied to each measurement point (406) before the gridding process, and rearranged in the
以上説明したように、本実施形態では、非直交系計測により得られたサンプリングデータをグリッディングする際に、各サンプリングデータに、距離rの要素にフリップ角FAの要素を加味した重み係数Wを作用させる。フリップ角FAの要素として、フリップ角FAが大きくなるにつれて値が大きくなる関数を採用する。従って、グリッディング時に、高いフリップ角FAのRFパルスにより得られたエコー信号、つまり、SNRの良好なエコー信号の寄与が高まる。これにより、体動などの影響を受けにくい非直交系計測において、フリップ角を変化させてSARを低減させながら、画像SNRを向上させることができる。 As described above, in the present embodiment, when the sampling data obtained by the non-orthogonal system measurement is gridded, the weighting factor W in which the factor of the flip angle FA is added to the element of the distance r is added to each sampling data. Make it work. As an element of the flip angle FA, a function whose value increases as the flip angle FA increases is adopted. Therefore, at the time of gridding, the contribution of an echo signal obtained by an RF pulse having a high flip angle FA, that is, an echo signal having a good SNR is increased. As a result, in non-orthogonal measurement that is not easily affected by body movement or the like, the image SNR can be improved while the SAR is reduced by changing the flip angle.
なお、上述のように、本実施形態では、セグメント毎にフリップ角を変化させているため、各セグメントで得られるエコー信号の信号強度は異なる。このフリップ角の異なるエコー信号間の信号強度差を補正するため、グリッディング処理の前に振幅補正を行うよう構成してもよい。すなわち、異なるフリップ角の振幅同士を規格化する。ここでは、異なる信号強度i番目のエコー信号に補正係数Ci=Ast/Ai(ここで、Aiはi番目エコー信号の最大振幅、Astは基準とするエコー信号の最大振幅とする。)を乗算することで、全エコー信号を基準とする振幅に揃える。 As described above, in this embodiment, since the flip angle is changed for each segment, the signal intensity of the echo signal obtained in each segment is different. In order to correct a signal intensity difference between echo signals having different flip angles, an amplitude correction may be performed before the gridding process. That is, the amplitudes of different flip angles are normalized. Here, the i-th echo signal having a different signal intensity is multiplied by a correction coefficient Ci = Ast / Ai (where Ai is the maximum amplitude of the i-th echo signal and Ast is the maximum amplitude of the reference echo signal). As a result, the amplitudes are set so that all echo signals are used as a reference.
振幅補正を含む場合のデータ再構成の手順を図7に示す。ここでは、異なる信号強度のエコー信号から得られたサンプリングデータ1211、1212、1213を、ひとつの計測空間内に配置する(ステップ1201)。全セグメントの信号を配置した計測空間には、フリップ角の異なるRFパルスによるエコー信号から得られたデータが混在する。上述の手法で振幅補正(ステップ1202)を行った後、グリッディング処理を行い(ステップ1203)、再配置後のデータに対してフーリエ変換を行い(ステップ1204)、画像化する(ステップ1205)。
FIG. 7 shows the data reconstruction procedure when amplitude correction is included. Here,
なお、振幅補正と同様の効果を得られるように、重み係数を変更してもよい。この場合
の重み係数W’は以下のとおりである。各エコー信号の振幅Aiは、磁化のxy平面成分
(横磁化:Mxy=Mz×sin(FA))の関数になることから、振幅補正係数Ciは
、以下の(式4)で表すことができる。
て、重み係数のフリップ角FAの関数w(FA)にフリップ角の関数で表される振幅補正
係数を乗算することにより、振幅補正と同様の効果が得られる。ただし、(式4)は、フ
リップ角FAが0に近づくと発散するため、(式5)で表される発散防止の係数cを導入
する。
の(式6)で表される。
また、距離rの関数w(r)と同様に、k空間原点からの距離R404によってフリップ角FAの関数w(FA)を変化させてもよい。図5(c)、(d)は、原点からの距離Rと、フリップ角FAとを変数とするw(FA)503、504の一例を示す図である。w(FA)503(または504)に示すように重み係数Wの中のフリップ角FAの関数部分w(FA)がフリップ角FAと原点からの距離Rとの関数になる(w(FA)=(w(FA)、w(R))。ここでRmaxはグリッディング最大半径である。w(FA)503は、R方向にはR−2の関数になっている。これは非直交系セグメント撮像の計測点密度がR方向にR−2に比例して減少していくためである。w(R)は、例えば、以下の(式8)で表される。
重み関数に用いるフリップ角の関数に、原点からの距離Rをさらに加味することにより、k空間の低空間周波数領域では、各エコー信号はフリップ角に応じて重み付けされるとともに、高空間周波数領域に向かうに従って重み係数がフリップ角に関して一定になり、すべての信号を同等にグリッディングに用いることができる。 By further adding the distance R from the origin to the function of the flip angle used for the weighting function, each echo signal is weighted according to the flip angle in the low spatial frequency region of k space, and in the high spatial frequency region. The weighting factor becomes constant with respect to the flip angle as it goes, so that all signals can be equally used for gridding.
<<第二の実施形態>>
次に、本発明の第二の実施形態を説明する。第一の実施形態では、非直交系計測の中でも、1枚の画像を構成するために必要なデータの取得を複数のセグメントに分割して行い、セグメント内ではフリップ角FAを一定としている。しかし、本実施形態では、フリップ角FAは連続的に変化する。以下、非直交系の時間分解血管造影法を例にあげて説明する。また、本実施形態において、MRI装置の構成は基本的に第一の実施形態と同様である。以下、第一の実施形態と異なる構成に主眼をおいて説明する。
<< Second Embodiment >>
Next, a second embodiment of the present invention will be described. In the first embodiment, even in non-orthogonal measurement, data necessary for constructing one image is divided into a plurality of segments, and the flip angle FA is constant within the segments. However, in the present embodiment, the flip angle FA changes continuously. Hereinafter, a non-orthogonal time-resolved angiographic method will be described as an example. In this embodiment, the configuration of the MRI apparatus is basically the same as that of the first embodiment. Hereinafter, a description will be given focusing on the configuration different from the first embodiment.
時間分解血管造影法は、ガドリニウム(Gd)などの造影剤の静脈注入と並行して高速撮影を連続的に行い、血管造影の様子を撮影する手法であり、高速撮影が可能なSSFPシーケンスが用いられる。本実施形態では、この時間分解血管造影法において、SARを低減するために、連続撮影間のフリップ角を連続的に変化させる。 The time-resolved angiography is a technique for continuously performing high-speed imaging in parallel with intravenous injection of a contrast medium such as gadolinium (Gd) and imaging the state of angiography, and uses an SSFP sequence capable of high-speed imaging. It is done. In the present embodiment, in this time-resolved angiography, the flip angle between consecutive imaging is continuously changed in order to reduce SAR.
本実施形態のフリップ角FAの連続変化を、フリップ角の変化関数FA(t)として定義する。変化関数FA(t)は、連続性、規則性を有する、例えば、正弦波形状のように規則的な形状とするのが望ましい。これは、フリップ角が急激に変化すると磁化の定常状態に影響し、画像にアーチファクトを招くためである。例えば、フリップ角FAの変域をFAmin〜FAmaxとする正弦波形状とすると、変化関数FA(t)は、以下の(式9)のように定義する。
本実施形態においても、第一の実施形態と同様に、フリップ角FAを変化させて取得したエコー信号のサンプリングデータのグリッディング処理時に、重み係数Wを作用させる。重み係数Wは、距離rの関数w(r)にフリップ角FAの要素を加味する。フリップ角FAの関数w(FA)は、たとえば、第一の実施形態の(式3)を用いる。本実施形態では、(式3)のFAに(式9)のFA(t)を代入する。 Also in the present embodiment, as in the first embodiment, the weighting factor W is applied during the gridding process of the sampling data of the echo signal acquired by changing the flip angle FA. The weighting factor W takes into account the factor of the flip angle FA to the function w (r) of the distance r. The function w (FA) of the flip angle FA uses, for example, (Equation 3) of the first embodiment. In the present embodiment, FA (t) in (Expression 9) is substituted for FA in (Expression 3).
このようにして、本実施形態においても、グリッディング時に、高いフリップ角FAのRFパルスにより得られたエコー信号の寄与が高まり、SNRの良好なエコー信号の寄与が高まる。これにより、体動などの影響を受けにくい非直交系計測において、フリップ角を変化させてSARを低減させながら、画像SNRを向上させることができる。 Thus, also in the present embodiment, the contribution of the echo signal obtained by the RF pulse having the high flip angle FA is increased during the gridding, and the contribution of the echo signal having a good SNR is increased. As a result, in non-orthogonal measurement that is not easily affected by body movement or the like, the image SNR can be improved while the SAR is reduced by changing the flip angle.
<<第三の実施形態>>
次に、本発明の第三の実施形態を説明する。第一および第二の実施形態では、各計測点を直交座標にグリッディングする際、重み係数として、直交系グリッドからの距離rの関数に、フリップ角FAの関数を加味している。本実施形態では、さらに、各計測点のデータを取得する時刻も加味する。以下、非直交系SSFPシーケンスのセグメント計測を例にあげて説明する。なお、本実施形態において、MRI装置の構成は基本的に第一の実施形態と同様である。以下、第一の実施形態と異なる構成に主眼をおいて説明する。
<< Third Embodiment >>
Next, a third embodiment of the present invention will be described. In the first and second embodiments, the function of the flip angle FA is added to the function of the distance r from the orthogonal grid as a weighting factor when each measurement point is gridded to the orthogonal coordinates. In the present embodiment, the time for acquiring the data of each measurement point is also taken into account. Hereinafter, segment measurement of a non-orthogonal SSFP sequence will be described as an example. In the present embodiment, the configuration of the MRI apparatus is basically the same as that of the first embodiment. Hereinafter, a description will be given focusing on the configuration different from the first embodiment.
図9は、本実施形態の計測を説明するための図であり、脂肪抑制パルス(CHEmical Shift Selective;CHESS)を併用した可変フリップ角SSFPシーケンスの模式図である。CHESSなどのプリパルス(601)を用いる計測の場合、その効果を保持するために計測を複数のセグメント(602〜604)に分割する。この他にもセグメント計測は、撮影を心電図などの生体信号に同期させる場合などに用いられる。図9(a)は、セグメント毎の取得したエコー信号と時刻tとの関係を示す。図9(b)は、1つのセグメントにおけるRFパルスの照射シーケンスを示す。図9(c)は、図9(b)に示す照射シーケンスにおけるフリップ角FAの変化を示す。図9(d)は、図9(b)に示す照射シーケンスにおける重み係数の時間の要素の変化を示す。 FIG. 9 is a diagram for explaining the measurement of this embodiment, and is a schematic diagram of a variable flip angle SSFP sequence using a fat suppression pulse (CHESS) together. In the case of measurement using a pre-pulse (601) such as CHESS, the measurement is divided into a plurality of segments (602 to 604) in order to maintain the effect. In addition to this, segment measurement is used when imaging is synchronized with a biological signal such as an electrocardiogram. FIG. 9A shows the relationship between the acquired echo signal for each segment and time t. FIG. 9B shows an irradiation sequence of RF pulses in one segment. FIG. 9C shows a change in the flip angle FA in the irradiation sequence shown in FIG. FIG. 9D shows a change in the time factor of the weighting coefficient in the irradiation sequence shown in FIG.
SSFPシーケンスでは磁化を定常化するため、計測開始前に複数のダミーパルス(605)を印加し、その後本計測(606)を開始する。ここでは計測セグメント内のフリップ角FA(607)は線形変化とする。なお、フリップ角FAの変化は、磁化を安定にできるものであれば線形変化に限られない。 In order to stabilize the magnetization in the SSFP sequence, a plurality of dummy pulses (605) are applied before the measurement is started, and then the main measurement (606) is started. Here, the flip angle FA (607) in the measurement segment is assumed to change linearly. Note that the change in the flip angle FA is not limited to a linear change as long as the magnetization can be stabilized.
セグメント計測には、セグメント毎にダミーパルスを照射する移行期間がある。一般に、移行期間中は信号を取得しないため、非セグメント計測に比べ撮影時間は延長し、SARも増大する。そこで、本実施形態では、撮影時間を短縮して撮影全体のSARを低減するために、ダミーパルス(605)によるエコー信号も読み出し、再構成に使用する。 Segment measurement includes a transition period in which a dummy pulse is irradiated for each segment. In general, since a signal is not acquired during the transition period, the imaging time is extended and the SAR is increased as compared with non-segment measurement. Therefore, in this embodiment, in order to shorten the imaging time and reduce the SAR of the entire imaging, the echo signal by the dummy pulse (605) is also read and used for reconstruction.
しかし、ダミーパルス照射期間は磁化が安定していないため、収集したエコー信号(データ)をそのまま計測空間に配置すると画質低下につながる。そこでこの間のデータを使用はするが、重視しないようにするため、グリッディングの重み係数Wに時間の要素を加える(608)。 However, since the magnetization is not stable during the dummy pulse irradiation period, if the collected echo signals (data) are arranged as they are in the measurement space, the image quality is degraded. Therefore, although data during this period is used, a time element is added to the weighting coefficient W of the gridding so as not to place importance (608).
具体的には、CHESSのプリパルス601のRF強度が最大となる時刻(609)をt=0、本計測が開始される時刻(610)をt=tstartとし、この間をダミーパルス期間(移行期間)とする。図9(d)に示すグラフ608のように、重み係数Wの時間の要素として、磁化の運動が安定でないダミーパルス期間の信号に対しては小さくし、定常状態の磁化からの信号に対しては大きくなるような関数W(t)を採用する。W(t)としては、例えば、以下の(式10)などが考えられる。
さらに、本実施形態においても、第一の実施形態同様、各計測点について、直交系グリッドからの距離rの関数w(r)およびフリップ角FAの関数w(FA)を作用させる。 Furthermore, also in this embodiment, as in the first embodiment, the function w (r) of the distance r from the orthogonal grid and the function w (FA) of the flip angle FA are applied to each measurement point.
以上説明したように、本実施形態によれば、磁化が安定していないダミーパルスによるエコー信号も画質を低下させる(アーチファクトを招く)ことなく画像の再構成に用いることができる。従って、上記各実施形態で得られる効果に加え、信号取得効率が向上し、全体の撮影時間を短縮することができるるため、さらにSARを低減することができる。 As described above, according to the present embodiment, an echo signal generated by a dummy pulse whose magnetization is not stable can also be used for image reconstruction without degrading the image quality (inducing artifacts). Therefore, in addition to the effects obtained in each of the above embodiments, the signal acquisition efficiency is improved and the entire imaging time can be shortened, so that the SAR can be further reduced.
<<第四の実施形態>>
次に、本発明の第四の実施形態を説明する。本実施形態は、非直交系撮像の中でも、ハイブリッドラディアル走査(プロペラスキャン)に本発明を適用するものである。なお、本実施形態において、MRI装置の構成は基本的に第一の実施形態と同様である。以下、第一の実施形態と異なる構成に主眼をおいて説明する。
<< Fourth Embodiment >>
Next, a fourth embodiment of the present invention will be described. In the present embodiment, the present invention is applied to hybrid radial scanning (propeller scanning) among non-orthogonal imaging. In the present embodiment, the configuration of the MRI apparatus is basically the same as that of the first embodiment. Hereinafter, a description will be given focusing on the configuration different from the first embodiment.
プロペラスキャンは、全計測を複数のブロック(ブレード)に分割し、各ブレードが計測空間上で異なる回転角となるよう計測する。各ブレード内では、1つの回転角に対して、通常計測の場合と同様に位相エンコードパルスを付加して複数のエコー信号を取得し、これを回転角毎に、すなわち、ブレード毎に繰り返して一枚の画像再構成に必要な全てのエコー信号を取得する。 In the propeller scan, all measurements are divided into a plurality of blocks (blades), and each blade is measured to have a different rotation angle in the measurement space. Within each blade, a plurality of echo signals are obtained by adding a phase encode pulse to one rotation angle as in the case of normal measurement, and this is repeated for each rotation angle, that is, for each blade. All echo signals necessary for image reconstruction of a sheet are acquired.
本実施形態では、1つのブレード内で取得されるエコー信号毎にフリップ角を変化させてSARの低減を図る。プロペラスキャンにおいても上述のようにフリップ角(および信号取得時刻)に応じた重み係数Wを設定することは、前述の各実施形態と同様である。しかし、プロペラスキャンの場合、ブレード内のエコー信号を取得する順番(走査順)の関係でフリップ角FAの変化関数はいくつかの態様を取り得る。 In this embodiment, the SAR is reduced by changing the flip angle for each echo signal acquired in one blade. In the propeller scan, setting the weighting factor W according to the flip angle (and the signal acquisition time) as described above is the same as in the above-described embodiments. However, in the case of a propeller scan, the change function of the flip angle FA can take several forms depending on the order (scanning order) in which echo signals in the blade are acquired.
図10は、本実施形態のエコー信号の配置とフリップ角FAの変化関数との関係を説明するための図である。ここでは、ブレード内のエコー信号数が5つの場合を例にあげて説明する。図10(a)は、プロペラスキャンの1つのブレードにおけるフリップ角FAの変化関数の一例である。本図に示すように、フリップ角の変化関数FA(t)が801のような形状を有し、ブレード内の5つのエコー信号をシーケンシャルに取得する場合を考える。図10(b)は、取得したエコー信号をk空間に配置したものである。ここで、最も高いSNRを持つエコー信号(803)が5番目のエコー信号であるとすると、本図に示すように、これはk空間においてブレードの端に配置され、k空間中心(k空間原点)には配置されない。 FIG. 10 is a diagram for explaining the relationship between the arrangement of echo signals and the change function of the flip angle FA in the present embodiment. Here, a case where the number of echo signals in the blade is five will be described as an example. FIG. 10A is an example of a change function of the flip angle FA in one blade of the propeller scan. As shown in this figure, a case is considered where the flip angle change function FA (t) has a shape such as 801 and five echo signals in the blade are sequentially acquired. FIG. 10B shows the acquired echo signals arranged in the k space. Here, assuming that the echo signal (803) having the highest SNR is the fifth echo signal, as shown in this figure, this is arranged at the end of the blade in the k space, and the center of the k space (k space origin) ) Is not placed.
従って、所望の(最もSNRの高い)エコー信号をk空間原点に配置するためには、ブレード内走査順とフリップ角FAの変化関数とのうち、いずれか一方を変更する。以下、所望の信号をk空間原点に配置する処理について説明する。 Therefore, in order to place a desired (highest SNR) echo signal at the k-space origin, either one of the scanning order within the blade and the change function of the flip angle FA is changed. Hereinafter, processing for arranging a desired signal at the k-space origin will be described.
図10(c)はフリップ角FAの変化関数を変更する例である。本図に示すように、k空間中心に最もフリップ角の大きなRFパルスによるエコー信号が配置されるよう、フリップ角FAの変化関数802の形状を変更する。ここでは、シーケンシャルに取得するエコー信号の中の中央のエコー信号(ここでは、3番目のエコー信号)のフリップ角が最も大きくなるような変化関数とする。フリップ角FAの変化関数の形状は、データがk空間に配置された際、k空間原点を通るエコー信号ほどフリップ角FAが大きくなるように定める。
FIG. 10C shows an example in which the change function of the flip angle FA is changed. As shown in this figure, the shape of the
図10(d)はブレード内のエコー配置順を変更する例である。この場合、フリップ角FAの変化関数は図10(a)のものを用いる。本図に示すように、k空間中心に最もフリップ角FAの大きなRFパルスの照射により得られるエコー信号が配置されるよう変更する。ここで、805は図10(a)の走査順に従って各エコー信号に付与された番号(エコー番号)である。ここでは、最も大きいフリップ角FAにより得られるエコー信号(ここでは、5番のエコー信号)がk空間の原点を通るよう走査順を変更する。なお、エコー信号の配置順は、データがk空間に配置された際、k空間原点近くを通るエコー信号ほどフリップ角FAが大きくなるよう定める。
FIG. 10D shows an example of changing the order of echo arrangement within the blade. In this case, the change function of the flip angle FA uses the one shown in FIG. As shown in the figure, the echo signal obtained by the irradiation of the RF pulse having the largest flip angle FA is arranged at the center of the k space. Here,
なお、変更を加える対象は、撮影の特性によって決定または選択する。連続計測の場合はフリップ角FAの変化関数を、セグメント計測の場合はブレード内走査順を変更するのが望ましい。 The object to be changed is determined or selected according to the characteristics of photographing. It is desirable to change the change function of the flip angle FA in the case of continuous measurement, and to change the scanning order within the blade in the case of segment measurement.
また、フリップ角FAの変化関数またはブレード内走査順の変更以外に計測開始位置のシフトによっても所望のエコー信号をk空間原点に配置できる。図11は計測開始位置のシフトによって所望のエコー信号がk空間原点を通るように配置する手順を説明するための図である。ここでは、全計測を8bladeに分割するものとする。また、フリップ角FAの変化関数は図10(a)に示すような時間軸に沿って漸増する関数であり、ブレード内をシーケンシャルに計測するものとする。 In addition to the change function of the flip angle FA or the change in the scanning order within the blade, a desired echo signal can be arranged at the k-space origin by shifting the measurement start position. FIG. 11 is a diagram for explaining a procedure for arranging a desired echo signal so as to pass through the k-space origin by shifting the measurement start position. Here, it is assumed that all measurements are divided into 8 blades. The change function of the flip angle FA is a function that gradually increases along the time axis as shown in FIG. 10A, and the inside of the blade is measured sequentially.
図11(a)は、計測開始位置をシフトさせない場合を示し、ここでは、エコー信号901を計測開始点とすると、最もSNRの高い(所望の)エコー信号902の配置位置がk空間原点から外れた位置となる。このとき、全8ブレードについてk空間走査を表すと903のようになる。図11(b)は、計測開始位置をシフトさせる場合を示す。計測開始位置(904)をブレード内の走査方向に直交する方向にシフトさせると、所望のエコー信号905は、k空間中心に配置される。このとき、全8ブレードについてk空間走査を表すと906のようになる。
FIG. 11A shows a case where the measurement start position is not shifted. Here, assuming that the
本実施形態では、以上のように、所望のエコー信号をk空間の原点に配置するよう計測条件を設定して計測を行った後、サンプリングデータと直交系グリッドとの距離r、フリップ角FA、時間tに応じた重み係数Wを作用させてグリッディング処理を行い、フーリエ変換し、画像を再構成する。 In the present embodiment, as described above, after performing measurement by setting measurement conditions so that a desired echo signal is placed at the origin of the k space, the distance r between the sampling data and the orthogonal grid, the flip angle FA, A gridding process is performed by applying a weighting factor W corresponding to the time t, Fourier transform is performed, and an image is reconstructed.
以上説明したように、本実施形態では、プロペラスキャンの場合であっても、グリッディング時に、高いフリップ角FAのRFパルスにより得られたエコー信号の寄与が高まり、SNRの良好なエコー信号の寄与が高まる。これにより、体動などの影響を受けにくい非直交系計測において、フリップ角を変化させてSARを低減させながら、画像SNRを向上させることができる。 As described above, in this embodiment, even in the case of a propeller scan, the contribution of an echo signal obtained by an RF pulse with a high flip angle FA increases during gridding, and the contribution of an echo signal having a good SNR. Will increase. As a result, in non-orthogonal measurement that is not easily affected by body movement or the like, the image SNR can be improved while the SAR is reduced by changing the flip angle.
上記各実施形態によれば、フリップ角を変化させる非直交系計測において、良好なSNRを得ることができる。従って、SARを低減させつつ良好なSNRを有し、かつ、アーチファクトの少ない画像を再構成できる。 According to each of the above embodiments, good SNR can be obtained in non-orthogonal measurement in which the flip angle is changed. Therefore, it is possible to reconstruct an image having a good SNR while reducing the SAR and having few artifacts.
なお、本発明は、パラレルイメージングにも適用することができる。パラレルイメージングではマルチプルコイルを用いる。本計測でエコー信号を間引いて取得し、マルチプルコイルの感度差を用いて折り返しを展開する。そのためパラレルイメージングでは、マルチプルコイルを構成する各コイルの感度分布を計測する必要があり、低いフリップ角を用いたシーケンスで各コイルの感度分布を本計測に先立って計測する。 The present invention can also be applied to parallel imaging. In parallel imaging, multiple coils are used. In this measurement, echo signals are thinned out and acquired, and aliasing is developed using the sensitivity difference of multiple coils. Therefore, in parallel imaging, it is necessary to measure the sensitivity distribution of each coil constituting the multiple coil, and the sensitivity distribution of each coil is measured prior to the main measurement in a sequence using a low flip angle.
図12は、本発明をパラレルイメージングに適用する場合の処理を説明するための図である。フリップ角を変更して取得した計測データ(1001〜1003)を合成して本計測データ(1004)とする。また、本計測データの中で低いフリップ角で取得したセグメントデータ(1001)を各コイルの感度データとして用い、パラレル再構成(1005)を行い、展開画像(1006)を得る。このとき、感度データにアーチファクトが生じないように、k空間にローパスフィルタを適用し、低域データのみを用いる。この手法によれば、パラレルイメージングにおいても、本計測の画像コントラストとは無関係に、常に低いフリップ角を用いたコイル感度に近いデータを取得できる。 FIG. 12 is a diagram for explaining processing when the present invention is applied to parallel imaging. The measurement data (1001 to 1003) acquired by changing the flip angle is combined into the main measurement data (1004). Further, segment data (1001) acquired at a low flip angle in the main measurement data is used as sensitivity data of each coil, and parallel reconstruction (1005) is performed to obtain a developed image (1006). At this time, a low-pass filter is applied to the k-space so that no artifact occurs in the sensitivity data, and only the low-frequency data is used. According to this method, even in parallel imaging, data close to the coil sensitivity using a low flip angle can be obtained regardless of the image contrast of the main measurement.
また、上記第一の実施形態で説明した振幅補正については、第二〜第四のいずれの実施形態においても適用できる。また、フリップ角FAの関数をk空間原点からの距離Rに応じて変化させる構成についても同様である。 The amplitude correction described in the first embodiment can be applied to any of the second to fourth embodiments. The same applies to the configuration in which the function of the flip angle FA is changed according to the distance R from the k-space origin.
1…被検体、2…静磁場発生系、3…傾斜磁場発生系、4…シーケンサ、5…送信系、6…受信系、7…信号処理系、8…中央処理装置(CPU)、9…傾斜磁場コイル、10…傾斜磁場電源、11…高周波発振器、12…変調器、13…高周波増幅器、14a…高周波コイル(送信コイル)、14b…高周波コイル(受信コイル)、15…信号増幅器、16…直交位相検波器、17…A/D変換器、18…磁気ディスク、19…光ディスク、20…ディスプレイ、21…ROM、22…RAM、23…トラックボール又はマウス、24…キーボード、51…ガントリ、52…テーブル、53…筐体、54…処理装置
DESCRIPTION OF
Claims (12)
前記非直交系計測空間データを直交系計測空間データに再配置するグリッディング手段と、
前記直交系計測空間データを再構成して画像化する画像形成手段と、
を備える磁気共鳴イメージング装置であって、
前記撮影手段は、前記非直交系計測空間データの収集に用いる高周波磁場パルスに複数種のフリップ角を設定するフリップ角設定手段を備え、
前記グリッディング手段は、フリップ角の異なる非直交系計測空間データの重み付け加算を用いて前記再配置を行うこ
を特徴とする磁気共鳴イメージング装置。 An imaging means for controlling the high-frequency magnetic field pulse and the gradient magnetic field pulse and measuring an echo signal that becomes non-orthogonal measurement space data
Gridding means for rearranging the non-orthogonal measurement space data into orthogonal measurement space data;
Image forming means for reconstructing and imaging the orthogonal measurement space data ;
A magnetic resonance imaging apparatus comprising:
The imaging means includes a flip angle setting means for setting a plurality of types of flip angles in the high-frequency magnetic field pulse used for collecting the non-orthogonal measurement space data,
The magnetic resonance imaging apparatus characterized in that the gridding means performs the rearrangement using weighted addition of non-orthogonal measurement space data having different flip angles.
前記グリッディング手段は、当該非直交系計測空間データの収集に用いた前記高周波磁場パルスのフリップ角に応じた重み係数を用いて前記重み付け加算を行うこと
を特徴とする磁気共鳴イメージング装置。 The magnetic resonance imaging apparatus according to claim 1,
The magnetic resonance imaging apparatus characterized in that the gridding means performs the weighted addition using a weighting coefficient corresponding to a flip angle of the high-frequency magnetic field pulse used for collecting the non-orthogonal measurement space data .
前記撮影手段は、前記エコー信号の計測を複数のセグメントに分けて実行し、
前記フリップ角設定手段は、前記フリップ角を、セグメント毎に変化させること
を特徴とする磁気共鳴イメージング装置。 The magnetic resonance imaging apparatus according to claim 1 or 2,
The imaging means executes the measurement of the echo signal divided into a plurality of segments,
The magnetic resonance imaging apparatus characterized in that the flip angle setting means changes the flip angle for each segment .
前記フリップ角設定手段は、前記フリップ角を連続的に変化させること
を特徴とする磁気共鳴イメージング装置。 The magnetic resonance imaging apparatus according to claim 1 or 2 ,
The magnetic resonance imaging apparatus, wherein the flip angle setting means continuously changes the flip angle .
前記グリッディング手段は、さらに、前記エコー信号の計測時刻に応じた重み係数を用いて前記重み付け加算を行うこと
を特徴とする磁気共鳴イメージング装置。 A magnetic resonance imaging apparatus according to any one of claims 1 to 4 ,
The magnetic resonance imaging apparatus , wherein the gridding means further performs the weighted addition using a weighting coefficient corresponding to a measurement time of the echo signal .
前記撮影手段は、
全計測を計測空間における回転角の異なる複数のブレードの計測に分割し、
前記複数のブレードの各々で位相エンコード量を変えた複数のエコー信号を取得する計測を行い、
前記フリップ角設定手段は、少なくとも一つの前記ブレード内でフリップ角を変更すること
を特徴とする磁気共鳴イメージング装置。 A magnetic resonance imaging apparatus according to any one of claims 1 to 5,
The photographing means includes
Dividing all measurements into measurements of multiple blades with different rotation angles in the measurement space,
Perform measurement to obtain a plurality of echo signals with different phase encoding amounts in each of the plurality of blades,
The magnetic resonance imaging apparatus characterized in that the flip angle setting means changes the flip angle in at least one of the blades .
前記フリップ角設定手段は、前記フリップ角を前記ブレードの中央で所望の値となるように単調に変化させること
を特徴とする磁気共鳴イメージング装置。 The magnetic resonance imaging apparatus according to claim 6 ,
The magnetic resonance imaging apparatus characterized in that the flip angle setting means monotonously changes the flip angle so as to become a desired value at the center of the blade .
前記撮影手段は、前記ブレード内で変更されるフリップ角が最大となるときに計測されるエコー信号が計測空間の原点を通る軌跡のデータとなるように、前記ブレード内で計測するエコー信号の計測順序を制御すること
を特徴とする磁気共鳴イメージング装置。 The magnetic resonance imaging apparatus according to claim 6 ,
The imaging means measures the echo signal measured in the blade so that the echo signal measured when the flip angle changed in the blade becomes the maximum becomes the trajectory data passing through the origin of the measurement space. A magnetic resonance imaging apparatus characterized by controlling an order .
前記撮影手段は、前記ブレード内で変更されるフリップ角が最大となるときに計測されるエコー信号が、計測空間の原点を通る軌跡のデータとなるように、前記計測空間における前記ブレードの位置を制御すること
を特徴とする磁気共鳴イメージング装置。 The magnetic resonance imaging apparatus according to claim 6 ,
The imaging means sets the position of the blade in the measurement space so that an echo signal measured when the flip angle changed in the blade becomes maximum becomes data of a trajectory passing through the origin of the measurement space. A magnetic resonance imaging apparatus characterized by controlling .
前記グリッディング手段は、前記重み付け加算に用いる重み係数を作用させる前に、前記エコー信号の振幅補正を行うこと
を特徴とする磁気共鳴イメージング装置。 A magnetic resonance imaging apparatus according to any one of claims 1 to 9 ,
The magnetic resonance imaging apparatus characterized in that the gridding means corrects the amplitude of the echo signal before applying a weighting coefficient used for the weighted addition .
前記グリッディング手段は、前記フリップ角に前記計測空間の原点からの距離を加味し、前記重み付け加算に用いる重み係数を設定すること
を特徴とする磁気共鳴イメージング装置。 The magnetic resonance imaging apparatus according to any one of claims 1 to 10 ,
The magnetic resonance imaging apparatus characterized in that the gridding means sets a weighting coefficient used for the weighted addition by adding a distance from an origin of the measurement space to the flip angle .
前記グリッディング手段は、前記重み付け加算に用いる重み係数に、前記非直交系計測空間データと、当該非直交系計測空間データが再配置される直交系計測空間データとの間の距離をさらに加味すること
を特徴とする磁気共鳴イメージング装置。 A magnetic resonance imaging apparatus according to any one of claims 1 to 11,
The gridding means further adds a distance between the non-orthogonal measurement space data and the orthogonal measurement space data to which the non-orthogonal measurement space data is rearranged, to the weighting coefficient used for the weighted addition. magnetic resonance imaging apparatus characterized by.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007147865A JP5230967B2 (en) | 2007-06-04 | 2007-06-04 | Nuclear magnetic resonance imaging system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007147865A JP5230967B2 (en) | 2007-06-04 | 2007-06-04 | Nuclear magnetic resonance imaging system |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2008295925A JP2008295925A (en) | 2008-12-11 |
JP2008295925A5 JP2008295925A5 (en) | 2010-07-15 |
JP5230967B2 true JP5230967B2 (en) | 2013-07-10 |
Family
ID=40170009
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2007147865A Active JP5230967B2 (en) | 2007-06-04 | 2007-06-04 | Nuclear magnetic resonance imaging system |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5230967B2 (en) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5735916B2 (en) * | 2009-08-12 | 2015-06-17 | 株式会社日立メディコ | Magnetic resonance imaging apparatus and synchronous measurement method |
JP5613065B2 (en) * | 2010-01-21 | 2014-10-22 | 株式会社東芝 | Magnetic resonance imaging system |
WO2011149087A1 (en) | 2010-05-27 | 2011-12-01 | 株式会社東芝 | Magnetic resonance imaging device |
WO2012005137A1 (en) * | 2010-07-07 | 2012-01-12 | 株式会社 日立メディコ | Magnetic resonance imaging device and rf pulse control method |
WO2012137606A1 (en) * | 2011-04-06 | 2012-10-11 | 株式会社 日立メディコ | Magnetic resonance imaging device and radial sampling method |
JP2012254189A (en) * | 2011-06-09 | 2012-12-27 | Hitachi Medical Corp | Magnetic resonance imaging apparatus |
JP6411063B2 (en) | 2014-05-07 | 2018-10-24 | キヤノンメディカルシステムズ株式会社 | Magnetic resonance imaging apparatus and SAR calculation method |
US11893662B2 (en) * | 2020-08-04 | 2024-02-06 | Northwestern University | Density compensation function in filtered backprojection |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3405813B2 (en) * | 1994-05-31 | 2003-05-12 | 株式会社東芝 | Magnetic resonance imaging method and magnetic resonance imaging apparatus |
DE10219528A1 (en) * | 2002-05-02 | 2003-11-13 | Philips Intellectual Property | Fast nuclear magnetic resonance imaging with gradient echoes |
JP4386354B2 (en) * | 2004-07-12 | 2009-12-16 | 株式会社日立メディコ | Magnetic resonance imaging device |
JP5179182B2 (en) * | 2005-07-27 | 2013-04-10 | 株式会社日立メディコ | Magnetic resonance imaging system |
JP2007275481A (en) * | 2006-04-12 | 2007-10-25 | Hitachi Medical Corp | Magnetic resonance imaging device |
US7868618B2 (en) * | 2006-06-16 | 2011-01-11 | Hitachi Medical Corporation | Magnetic resonance imaging apparatus |
-
2007
- 2007-06-04 JP JP2007147865A patent/JP5230967B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
JP2008295925A (en) | 2008-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5230967B2 (en) | Nuclear magnetic resonance imaging system | |
JP6513398B2 (en) | MR image reconstruction using prior information constrained regularization | |
EP1444530B1 (en) | Three-dimensional phase contrast magnetic resonance imaging using interleaved projection-reconstruction data | |
US8466679B2 (en) | Magnetic resonance imaging apparatus and method configured for susceptibility-emphasized imaging with improved signal-to-noise ratio | |
JP5726203B2 (en) | Magnetic resonance imaging apparatus, irradiation magnetic field measurement method | |
WO2011007691A1 (en) | Magnetic resonance imaging device and magnetic resonance imaging method | |
JP5735916B2 (en) | Magnetic resonance imaging apparatus and synchronous measurement method | |
US9513356B2 (en) | Magnetic resonance imaging apparatus and reconstructed image acquisition method | |
JP5337406B2 (en) | Magnetic resonance imaging system | |
WO2010064572A1 (en) | Magnetic resonance imaging device and synchronous imaging method | |
US9594142B2 (en) | Magnetic resonance imaging apparatus and echo signal measurement method | |
JP2008055023A (en) | Magnetic resonance imaging apparatus | |
JP5808659B2 (en) | Magnetic resonance imaging apparatus and T1ρ imaging method | |
JP3847554B2 (en) | Magnetic resonance imaging system | |
JP5164443B2 (en) | Magnetic resonance imaging system | |
JP4018107B2 (en) | Magnetic resonance imaging system | |
JP5637694B2 (en) | Magnetic resonance imaging apparatus and non-orthogonal coordinate system scanning method | |
JP6084430B2 (en) | Flip angle control method of refocusing RF pulse and magnetic resonance imaging apparatus | |
JP5421600B2 (en) | Nuclear magnetic resonance imaging apparatus and method of operating nuclear magnetic resonance imaging apparatus | |
JP2011036455A (en) | Magnetic resonance imaging apparatus | |
JP2012254189A (en) | Magnetic resonance imaging apparatus | |
US11914017B2 (en) | Magnetic resonance imaging apparatus and image processing method | |
JP2014087442A5 (en) | ||
JP5758230B2 (en) | Magnetic resonance imaging apparatus and inversion RF pulse phase control method | |
US20220015661A1 (en) | Magnetic resonance imaging apparatus and method of controlling the same |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20100422 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100527 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20100527 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20120227 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20120229 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20120402 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20121029 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20121225 |
|
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: 20130311 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20130321 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20160329 Year of fee payment: 3 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5230967 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |