JP5230967B2 - Nuclear magnetic resonance imaging system - Google Patents

Nuclear magnetic resonance imaging system Download PDF

Info

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
Application number
JP2007147865A
Other languages
Japanese (ja)
Other versions
JP2008295925A (en
JP2008295925A5 (en
Inventor
康弘 鎌田
博幸 板垣
将宏 瀧澤
哲彦 高橋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Healthcare Manufacturing Ltd
Original Assignee
Hitachi Medical Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hitachi Medical Corp filed Critical Hitachi Medical Corp
Priority to JP2007147865A priority Critical patent/JP5230967B2/en
Publication of JP2008295925A publication Critical patent/JP2008295925A/en
Publication of JP2008295925A5 publication Critical patent/JP2008295925A5/ja
Application granted granted Critical
Publication of JP5230967B2 publication Critical patent/JP5230967B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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)の関係がある。

Figure 0005230967
ここでBは静磁場強度[T]、αはフリップ角[rad]、Dはデューティーサイクル(単位時間当たりのRFパルス照射時間)である。(式1)より、SARは静磁場強度の2乗に比例するため、MRI装置の高磁場化によりSARの影響はより顕著なものとなる。また、デューティーサイクルに比例するため、特に180度のリフォーカスRFパルスを高速に照射する高速スピンエコーシーケンスや、短い繰り返し時間(TR)で励起RFパルスを照射するSSFPシーケンスなどでSARが大きくなる。 In imaging using an MRI apparatus, if the specific absorption rate (SAR), which is the rate at which the RF pulse energy used for spin excitation is absorbed by the human body, increases, the human body is affected by heat generation and the like. . There is a relationship of the following (formula 1) between the SAR and the imaging parameter.
Figure 0005230967
Here, B 0 is the static magnetic field strength [T], α is the flip angle [rad], and D is the duty cycle (RF pulse irradiation time per unit time). From (Equation 1), since the SAR is proportional to the square of the static magnetic field strength, the influence of the SAR becomes more remarkable by increasing the magnetic field of the MRI apparatus. In addition, since it is proportional to the duty cycle, the SAR is increased particularly in a high-speed spin echo sequence that irradiates a 180-degree refocus RF pulse at a high speed or an SSFP sequence that irradiates an excitation RF pulse with a short repetition time (TR).

一般に、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 Document 1 changes the flip angle (flop angle) of the refocus pulse constituting the fast spin echo sequence, and Patent Document 1 changes the flip angle of the excitation pulse of the SSFP sequence according to the phase encoding amount. A technique for reducing the SAR of the entire image is disclosed. Here, since the echo signal arranged in the low spatial frequency region of the measurement space (so-called k space) strongly contributes to the SNR of the reconstructed image, a signal with a high SNR and a high flip angle is converted into a low spatial frequency. Arranged in the area prevents the SNR of the reconstructed image from decreasing.

“Hyperechoes”、 Magnetic Resonance in Medicine 46: 6-12 (2001)“Hyperechoes”, Magnetic Resonance in Medicine 46: 6-12 (2001) 特表2005−524453号公報JP 2005-524453 A

前述したように、高磁場の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 field generation system 2, a gradient magnetic field generation system 3, a transmission system 5, and a reception system 6. A signal processing system 7, a sequencer 4, and a central processing unit (CPU) 8.

静磁場発生系2は、垂直磁場方式であれば、被検体1の周りの空間にその体軸と直交する方向に、水平磁場方式であれば、体軸方向に均一な静磁場を発生させるもので、被検体1の周りに永久磁石方式、常電導方式あるいは超電導方式の静磁場発生源が配置されている。   The static magnetic field generation system 2 generates a uniform static magnetic field in the direction perpendicular to the body axis in the space around the subject 1 if the vertical magnetic field method is used, and in the body axis direction if the horizontal magnetic field method is used. Thus, a permanent magnet type, normal conduction type or superconducting type static magnetic field generation source is arranged around the subject 1.

傾斜磁場発生系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 field generation system 3 includes a gradient magnetic field coil 9 wound in the three-axis directions of X, Y, and Z, which is a coordinate system (stationary coordinate system) of the MRI apparatus, and a gradient magnetic field power source 10 that drives each gradient magnetic field coil. The gradient magnetic fields Gx, Gy, and Gz are applied in the three axial directions of X, Y, and Z by driving the gradient magnetic field power supply 10 of each coil in accordance with a command from the sequencer 4 described later. At the time of imaging, a slice direction gradient magnetic field pulse (Gs) is applied in a direction orthogonal to the slice plane (imaging cross section) to set a slice plane for the subject 1, and the remaining two orthogonal to the slice plane and orthogonal to each other. A phase encoding direction gradient magnetic field pulse (Gp) and a frequency encoding direction gradient magnetic field pulse (Gf) are applied in one direction, and position information in each direction is encoded in the echo signal.

シーケンサ4は、高周波磁場パルス(以下、RFパルスという)と傾斜磁場パルスとを所定のパルスシーケンスに従って繰り返し印加するよう各部を制御する。CPU8の制御に従って動作し、被検体1の断層画像のデータ収集に必要な種々の命令を送信系5、傾斜磁場発生系3、および受信系6に送る。   The sequencer 4 controls each unit so as to repeatedly apply a high-frequency magnetic field pulse (hereinafter referred to as an RF pulse) and a gradient magnetic field pulse according to a predetermined pulse sequence. It operates according to the control of the CPU 8, and sends various commands necessary for collecting tomographic image data of the subject 1 to the transmission system 5, the gradient magnetic field generation system 3, and the reception system 6.

送信系5は、被検体1の生体組織を構成する原子の原子核スピンに核磁気共鳴を起こさせるために、被検体1にRFパルスを照射するもので、高周波発振器11と変調器12と高周波増幅器13と送信側の高周波コイル(送信コイル)14aとを備える。送信系5は、高周波発振器11から出力された高周波パルスをシーケンサ4からの指令によるタイミングで変調器12によって振幅変調し、振幅変調された高周波パルスを高周波増幅器13で増幅した後に被検体1に近接して配置された送信コイル14aに供給することにより、RFパルスを被検体1に照射する。   The transmission system 5 irradiates the subject 1 with an RF pulse in order to cause nuclear magnetic resonance to occur in the nuclear spins of the atoms constituting the living tissue of the subject 1, and includes a high-frequency oscillator 11, a modulator 12, and a high-frequency amplifier. 13 and a high frequency coil (transmission coil) 14a on the transmission side. The transmission system 5 amplitude-modulates the high-frequency pulse output from the high-frequency oscillator 11 by a modulator 12 at a timing according to a command from the sequencer 4, amplifies the amplitude-modulated high-frequency pulse by the high-frequency amplifier 13, and then approaches the subject 1. Then, the object 1 is irradiated with an RF pulse by being supplied to the transmitting coil 14a.

受信系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 subject 1. The receiving system 6 receives a high-frequency coil (receiving coil) 14 b on the receiving side and a signal amplifier 15. And a quadrature phase detector 16 and an A / D converter 17. The echo signal of the response of the subject 1 induced by the electromagnetic wave irradiated from the transmission coil 14a is detected by the reception coil 14b arranged close to the subject 1, amplified by the signal amplifier 15, and then from the sequencer 4. Are divided into two orthogonal signals by the quadrature detector 16, converted into digital quantities by the A / D converter 17, and sent to the signal processing system 7.

信号処理系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 optical disk 19 and a magnetic disk 18, a memory such as a ROM 21 and a RAM 22, and a display 20 including a CRT. With. The CPU 8 displays a tomographic image of the subject 1 as a result of the signal processing, image reconstruction, and the like performed on the data received from the receiving system 6 on the display 20 and records it on the magnetic disk 18 of the external storage device. .

操作部25は、MRI装置100の各種制御情報や信号処理系7で行う処理の制御情報の入力インタフェースで、トラックボール又はマウス23、及び、キーボード24を備える。操作部25はディスプレイ20に近接して配置され、操作者がディスプレイ20を見ながら操作部25を通してインタラクティブにMRI装置100の各種処理を制御する。   The operation unit 25 is an input interface for various control information of the MRI apparatus 100 and control information of processing performed by the signal processing system 7, and includes a trackball or mouse 23 and a keyboard 24. The operation unit 25 is disposed in the vicinity of the display 20, and an operator controls various processes of the MRI apparatus 100 interactively through the operation unit 25 while looking at the display 20.

なお、送信コイル14aと傾斜磁場コイル9は、被検体1が挿入される静磁場発生系2の静磁場空間内に、垂直磁場方式であれば被検体1に対向して、水平磁場方式であれば被検体1を取り囲むようにして設置される。また、受信コイル14bは、被検体1に対向して、或いは取り囲むように設置されている。現在、MRI装置の撮像対象核種であって臨床で普及しているものは、被検体の主たる構成物質である水素原子核(プロトン)である。プロトン密度の空間分布や、励起状態の緩和時間の空間分布に関する情報を画像化することで、人体頭部、腹部、四肢等の形態または、機能を2次元もしくは3次元的に撮像する。   Note that the transmission coil 14a and the gradient magnetic field coil 9 may be opposed to the subject 1 in the static magnetic field space of the static magnetic field generation system 2 in which the subject 1 is inserted, in the horizontal magnetic field method. For example, it is installed so as to surround the subject 1. The receiving coil 14b is disposed so as to face or surround the subject 1. At present, the nuclide to be imaged by the MRI apparatus, which is widely used clinically, is a hydrogen nucleus (proton) which is a main constituent material of the subject. Information on the spatial distribution of the proton density and the spatial distribution of the relaxation time of the excited state is imaged, thereby imaging the form or function of the human head, abdomen, limbs, etc. two-dimensionally or three-dimensionally.

次に、一般的な非直交系計測におけるデータの再構成手順を説明する。ここでは、非直交系計測の中のラディアルサンプリング計測と呼ばれる、計測空間の略一点を回転中心として回転角を変えながら放射状にサンプリングを行い、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). Data 311, 312, 313 acquired in units of segments are arranged in one measurement space (step 301). Then, interpolation processing for converting non-orthogonal sampling data called gridding into orthogonal sampling data (hereinafter referred to as gridding processing) is performed on the data arranged in the measurement space, and the data is non-orthogonal coordinate system. Is converted into an orthogonal coordinate system (orthogonal grid) and rearranged in the k space (step 302), and Fourier transform is performed on the rearranged data (step 303) to form an image (step 304).

グリッディング処理では、直交座標系の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 segment 211 is FA1, the flip angle in the segment 212 is FA2, and the flip angle in the segment 213 is FA3, FA1> FA2> FA3 is changed. In this embodiment, since the flip angle FA is changed in each segment, the echo signals (measurement space data) 221, 222, and 223 acquired in the segments 211, 212, and 213 have different SNRs. Sampling data obtained from echo signals 221, 222, and 223 having different SNRs are arranged in measurement spaces as shown in FIG. Note that the change in the flip angle is not limited to that described above. However, in order to keep the magnetization motion stable, it is desirable to gradually increase or decrease it.

次に、本実施形態の可変フリップ角非直交系セグメント計測におけるデータの再構成手順を説明する。なお、以下の処理は、信号処理系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 ROM 21 or the RAM 22 in accordance with a program stored in the ROM 21 or the RAM 22 in advance, or in accordance with an instruction from the user input via the operation unit 25. Or it performs using the data input from the user via the operation part. FIG. 4 is a diagram for explaining a data reconstruction procedure in variable flip angle non-orthogonal segment measurement. The procedure of data reconstruction is basically the same as the procedure in the general non-orthogonal system measurement described above.

本実施形態では、セグメント毎にフリップ角を変化させているため、各セグメントで得られるエコー信号の信号強度は異なる。異なる信号強度のエコー信号から得られたサンプリングデータ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. Sampling data 1111, 1112, 1113 obtained from echo signals having different signal intensities are arranged in one measurement space (step 1101). Then, gridding is performed and rearranged in the k space (step 1102), and the rearranged data is subjected to Fourier transform (step 1103) and imaged (step 1104).

上述のように本実施形態ではセグメント毎にフリップ角を変化させているため、全セグメントの信号を配置した計測空間には、フリップ角の異なる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)を用いる。

Figure 0005230967
ここで、rmaxは、グリッディング最大半径である。 As the function w (r) of the distance r, a function whose value decreases as the distance from the orthogonal grid increases is adopted. For example, the following (Formula 2) is used.
Figure 0005230967
Here, r max is a gridding maximum radius.

図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 graph 501 of the function w (r) represented by (Expression 2). When such a function W (r) is used, the measurement points closer to the orthogonal grid at the time of gridding are more important. Note that as the distance R from the k-space origin to the orthogonal grid increases, the density of measurement points decreases, and the number of measurement points in a circle with a radius r decreases. Therefore, the weighting factor W (r) may be changed according to the distance R from the origin.

また、フリップ角FAの関数w(FA)は、高いフリップ角のRFパルスによるエコー信号に対しては大きく、低いフリップ角のそれに対しては小さく作用するものを採用する。これは、フリップ角の高いRFパルスによるエコー信号からの寄与を高め、画像SNRを向上させるためである。例えば、関数w(FA)として、例えば、以下の(式3)を用いる。

Figure 0005230967
ここで、FAmaxは、フリップ角FAの最大値である。 A function w (FA) of the flip angle FA that is large for an echo signal by an RF pulse with a high flip angle and a function that acts small for that of a low flip angle is adopted. This is to increase the contribution from the echo signal due to the RF pulse having a high flip angle and improve the image SNR. For example, the following (Formula 3) is used as the function w (FA), for example.
Figure 0005230967
Here, FA max is the maximum value of the flip angle FA.

図5(b)は、(式3)で示されるフリップ角FAの関数w(FA)のグラフ502である。ここで、wminはフリップ角FAが0のときの重み係数を示す。wminはFAmax時のw(FA)の値(ここでは1)より小さい値をとり、0としても良い。本図に示すw(FA)は、SNRの高い高フリップ角のRFパルスによるエコー信号に対して大きく、低フリップ角のエコー信号に関しては小さい。従って、フリップ角が大きい(高い)サンプリングデータほど重視される。 FIG. 5B is a graph 502 of the function w (FA) of the flip angle FA expressed by (Equation 3). Here, w min indicates a weighting coefficient when the flip angle FA is zero. w min takes a value smaller than the value of w (FA) at FA max (here, 1), and may be 0. The w (FA) shown in the figure is large with respect to an echo signal with a high flip angle RF pulse having a high SNR, and is small with respect to an echo signal with a low flip angle. Accordingly, sampling data with a larger (higher) flip angle is more important.

次に、本実施形態の重み係数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 orthogonal grid 408. To do.

以上説明したように、本実施形態では、非直交系計測により得られたサンプリングデータをグリッディングする際に、各サンプリングデータに、距離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, sampling data 1211, 1212, and 1213 obtained from echo signals having different signal intensities are arranged in one measurement space (step 1201). In the measurement space where the signals of all segments are arranged, data obtained from echo signals by RF pulses having different flip angles are mixed. After performing amplitude correction (step 1202) by the above-described method, gridding processing is performed (step 1203), Fourier transform is performed on the rearranged data (step 1204), and imaging is performed (step 1205).

なお、振幅補正と同様の効果を得られるように、重み係数を変更してもよい。この場合
の重み係数W’は以下のとおりである。各エコー信号の振幅Aiは、磁化のxy平面成分
(横磁化:Mxy=M×sin(FA))の関数になることから、振幅補正係数Ciは
、以下の(式4)で表すことができる。

Figure 0005230967
すなわち、振幅補正係数Ciは、フリップ角の関数C(FA)と表すことができる。従っ
て、重み係数のフリップ角FAの関数w(FA)にフリップ角の関数で表される振幅補正
係数を乗算することにより、振幅補正と同様の効果が得られる。ただし、(式4)は、フ
リップ角FAが0に近づくと発散するため、(式5)で表される発散防止の係数cを導入
する。
Figure 0005230967
上記発散係数を用い、フリップ角の関数として表される振幅補正係数C(FA)は、以下
の(式6)で表される。
Figure 0005230967
従って、重み係数W’は以下の(式7)で表される。
Figure 0005230967
Note that the weighting coefficient may be changed so as to obtain the same effect as the amplitude correction. The weighting factor W ′ in this case is as follows. Since the amplitude Ai of each echo signal is a function of the xy plane component of magnetization (transverse magnetization: M xy = M z × sin (FA)), the amplitude correction coefficient Ci is expressed by the following (formula 4). Can do.
Figure 0005230967
That is, the amplitude correction coefficient Ci can be expressed as a flip angle function C (FA). Therefore, by multiplying the function w (FA) of the flip angle FA of the weighting coefficient by the amplitude correction coefficient expressed by the function of the flip angle, the same effect as the amplitude correction can be obtained. However, since (Equation 4) diverges when the flip angle FA approaches 0, the divergence prevention coefficient c expressed by (Equation 5) is introduced.
Figure 0005230967
Using the divergence coefficient, the amplitude correction coefficient C (FA) expressed as a function of the flip angle is expressed by the following (formula 6).
Figure 0005230967
Therefore, the weight coefficient W ′ is expressed by the following (Equation 7).
Figure 0005230967

また、距離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)で表される。

Figure 0005230967
Similarly to the function w (r) of the distance r, the function w (FA) of the flip angle FA may be changed by the distance R404 from the k-space origin. FIGS. 5C and 5D are diagrams showing an example of w (FA) 503 and 504 using the distance R from the origin and the flip angle FA as variables. As indicated by w (FA) 503 (or 504), the function part w (FA) of the flip angle FA in the weighting factor W is a function of the flip angle FA and the distance R from the origin (w (FA) = (W (FA), w (R)), where R max is the gridding maximum radius, and w (FA) 503 is a function of R −2 in the R direction, which is a non-orthogonal system. This is because the measurement point density of the segment imaging decreases in proportion to R- 2 in the R direction, and w (R) is expressed by, for example, (Equation 8) below.
Figure 0005230967

重み関数に用いるフリップ角の関数に、原点からの距離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)のように定義する。

Figure 0005230967
ここで角周波数ωはTR(繰り返し時間)や位相エンコード数などの撮影パラメータ、または、心拍や呼吸周期などの生体情報から決定してもよい。図8は、本実施形態の計測を説明するための図である。図8(a)は、(式9)で定義されるフリップ角FAの変化関数のグラフ701である。また、図8(b)は、各フリップ角で取得するエコー信号と時間tとの関係を示す。 The continuous change of the flip angle FA of this embodiment is defined as a flip angle change function FA (t). It is desirable that the change function FA (t) has a continuity and regularity, for example, a regular shape such as a sine wave shape. This is because if the flip angle changes rapidly, it affects the steady state of magnetization and causes artifacts in the image. For example, assuming a sinusoidal shape in which the range of the flip angle FA is FA min to FA max , the change function FA (t) is defined as in the following (Equation 9).
Figure 0005230967
Here, the angular frequency ω may be determined from imaging parameters such as TR (repetition time) and the number of phase encodings, or biological information such as heartbeat and respiratory cycle. FIG. 8 is a diagram for explaining the measurement of the present embodiment. FIG. 8A is a graph 701 of the change function of the flip angle FA defined by (Equation 9). FIG. 8B shows the relationship between the echo signal acquired at each flip angle and time t.

本実施形態においても、第一の実施形態と同様に、フリップ角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)などが考えられる。

Figure 0005230967
Specifically, the time (609) when the RF intensity of the CHESS pre-pulse 601 becomes maximum is t = 0, and the time (610) when the main measurement is started is t = t start, and this period is a dummy pulse period (transition period) ). As a graph 608 shown in FIG. 9D, as a time factor of the weighting factor W, a signal in a dummy pulse period in which the magnetization motion is not stable is reduced, and a signal from a steady-state magnetization is reduced. A function W (t) that increases is used. As W (t), for example, the following (Formula 10) can be considered.
Figure 0005230967

さらに、本実施形態においても、第一の実施形態同様、各計測点について、直交系グリッドからの距離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 change function 802 of the flip angle FA is changed so that the echo signal by the RF pulse having the largest flip angle is arranged at the center of the k space. Here, the change function is such that the flip angle of the central echo signal (here, the third echo signal) among the echo signals acquired sequentially is the largest. The shape of the change function of the flip angle FA is determined such that when the data is arranged in the k space, the flip angle FA increases as the echo signal passes through the k space origin.

図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, reference numeral 805 denotes a number (echo number) assigned to each echo signal in accordance with the scanning order of FIG. Here, the scanning order is changed so that the echo signal (here, the fifth echo signal) obtained by the largest flip angle FA passes through the origin of the k space. The order in which the echo signals are arranged is determined such that when the data is arranged in the k space, the flip angle FA increases as the echo signal passes near the origin of the k space.

なお、変更を加える対象は、撮影の特性によって決定または選択する。連続計測の場合はフリップ角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 echo signal 901 is the measurement start point, the arrangement position of the (desired) echo signal 902 having the highest SNR deviates from the k-space origin. It becomes the position. At this time, the k-space scanning for all the 8 blades is represented as 903. FIG. 11B shows a case where the measurement start position is shifted. When the measurement start position (904) is shifted in the direction orthogonal to the scanning direction in the blade, the desired echo signal 905 is placed at the center of the k-space. At this time, the k-space scanning for all 8 blades is represented by 906.

本実施形態では、以上のように、所望のエコー信号を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.

第一の実施形態のMRI装置のブロック図である。It is a block diagram of the MRI apparatus of 1st embodiment. 非直交系計測のデータ再構成の手順を説明するための図である。It is a figure for demonstrating the procedure of the data reconstruction of non-orthogonal system measurement. 第一の実施形態の計測を説明するための図である。It is a figure for demonstrating the measurement of 1st embodiment. 第一の実施形態のデータ再構成の手順を説明するための図である。It is a figure for demonstrating the procedure of the data reconstruction of 1st embodiment. 第一の実施形態の重み係数のグラフである。It is a graph of the weighting coefficient of 1st embodiment. 第一の実施形態のグリッディング処理の模式図である。It is a schematic diagram of the gridding process of 1st embodiment. 第一の実施形態の振幅補正を含む場合のデータ再構成の手順を説明するための図である。It is a figure for demonstrating the procedure of the data reconstruction in case the amplitude correction of 1st embodiment is included. 第二の実施形態の計測を説明するための図である。It is a figure for demonstrating the measurement of 2nd embodiment. 第三の実施形態の計測を説明するための図である。It is a figure for demonstrating the measurement of 3rd embodiment. 第四の実施形態のエコー信号の配置を変換する処理を説明するための図である。It is a figure for demonstrating the process which converts the arrangement | positioning of the echo signal of 4th embodiment. 第四の実施形態のエコー信号の配置を変換する処理の別の例を説明するための図である。It is a figure for demonstrating another example of the process which converts the arrangement | positioning of the echo signal of 4th embodiment. パラレルイメージングに適用した実施形態のデータ再構成の手順を説明するための図である。It is a figure for demonstrating the procedure of the data reconstruction of embodiment applied to parallel imaging.

符号の説明Explanation of symbols

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 SYMBOLS 1 ... Subject, 2 ... Static magnetic field generation system, 3 ... Gradient magnetic field generation system, 4 ... Sequencer, 5 ... Transmission system, 6 ... Reception system, 7 ... Signal processing system, 8 ... Central processing unit (CPU), 9 ... Gradient magnetic field coil, 10 Gradient magnetic field power supply, 11 High frequency oscillator, 12 Modulator, 13 High frequency amplifier, 14 a High frequency coil (transmitting coil), 14 b High frequency coil (receiving coil), 15 Signal amplifier, 16 Quadrature phase detector, 17 ... A / D converter, 18 ... magnetic disk, 19 ... optical disk, 20 ... display, 21 ... ROM, 22 ... RAM, 23 ... trackball or mouse, 24 ... keyboard, 51 ... gantry, 52 ... Table, 53 ... Case, 54 ... Processing device

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.
請求項1記載の磁気共鳴イメージング装置であって、
前記グリッディング手段は、当該非直交系計測空間データの収集に用いた前記高周波磁場パルスのフリップ角に応じた重み係数を用いて前記重み付け加算を行うこと
を特徴とする磁気共鳴イメージング装置。
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 .
請求項1又は2に記載の磁気共鳴イメージング装置であって、
前記撮影手段は、前記エコー信号の計測を複数のセグメントに分けて実行し、
前記フリップ角設定手段は、前記フリップ角を、セグメント毎に変化させること
を特徴とする磁気共鳴イメージング装置。
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 .
請求項1又は2に記載の磁気共鳴イメージング装置であって、
前記フリップ角設定手段は、前記フリップ角を連続的に変化させること
を特徴とする磁気共鳴イメージング装置。
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 .
請求項1乃至4のいずれか一項に記載の磁気共鳴イメージング装置であって、
前記グリッディング手段は、さらに、前記エコー信号の計測時刻に応じた重み係数を用いて前記重み付け加算を行うこと
を特徴とする磁気共鳴イメージング装置。
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 .
請求項乃至5のいずれか一項に記載の磁気共鳴イメージング装置であって、
前記撮影手段は、
全計測を計測空間における回転角の異なる複数のブレードの計測に分割し、
前記複数のブレードの各々で位相エンコード量を変えた複数のエコー信号を取得する計測を行い、
前記フリップ角設定手段は、少なくとも一つの前記ブレード内でフリップ角を変更すること
を特徴とする磁気共鳴イメージング装置。
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 .
請求項1乃至9のいずれか一項記に載の磁気共鳴イメージング装置であって、
前記グリッディング手段は、前記重み付け加算に用いる重み係数を作用させる前に、前記エコー信号の振幅補正を行うこと
を特徴とする磁気共鳴イメージング装置。
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 .
請求項乃至10のいずれか一項に記載の磁気共鳴イメージング装置であって、
前記グリッディング手段は、前記フリップ角に前記計測空間の原点からの距離を加味し、前記重み付け加算に用いる重み係数を設定すること
を特徴とする磁気共鳴イメージング装置。
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 .
請求項乃至11のいずれか一項に記載の磁気共鳴イメージング装置であって、
前記グリッディング手段は、前記重み付け加算に用いる重み係数に、前記非直交系計測空間データと、当該非直交系計測空間データが再配置される直交系計測空間データとの間の距離をさらに加味すること
を特徴とする磁気共鳴イメージング装置。
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.
JP2007147865A 2007-06-04 2007-06-04 Nuclear magnetic resonance imaging system Active JP5230967B2 (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

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