JP2013111097A - Magnetic resonance spectroscopy apparatus - Google Patents

Magnetic resonance spectroscopy apparatus Download PDF

Info

Publication number
JP2013111097A
JP2013111097A JP2011257122A JP2011257122A JP2013111097A JP 2013111097 A JP2013111097 A JP 2013111097A JP 2011257122 A JP2011257122 A JP 2011257122A JP 2011257122 A JP2011257122 A JP 2011257122A JP 2013111097 A JP2013111097 A JP 2013111097A
Authority
JP
Japan
Prior art keywords
magnetic resonance
ratio
unit
nuclear magnetic
resonance signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2011257122A
Other languages
Japanese (ja)
Other versions
JP5908704B2 (en
Inventor
Hitoshi Ikeda
仁 池田
Haruyuki Fukuchi
治之 福地
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.)
GE Medical Systems Global Technology Co LLC
Original Assignee
GE Medical Systems Global Technology Co LLC
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 GE Medical Systems Global Technology Co LLC filed Critical GE Medical Systems Global Technology Co LLC
Priority to JP2011257122A priority Critical patent/JP5908704B2/en
Publication of JP2013111097A publication Critical patent/JP2013111097A/en
Application granted granted Critical
Publication of JP5908704B2 publication Critical patent/JP5908704B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

PROBLEM TO BE SOLVED: To provide a magnetic resonance spectroscopy apparats providing a spectrogram having a fixed or more level of an SN (Signal Noise) ratio.SOLUTION: An apparatus (10) includes: a reception part (150) for receiving a first nuclear magnetic resonance signal from a specific area; a storage part (172) for storing the first nuclear magnetic resonance signal received by the reception part; and an image processing part (174) for calculating a frequency spectrum, based on the first nuclear magnetic resonance signal. The apparatus also includes: an SN ratio determination part (178) for determining whether the SN ratio of a specific substance in the frequency spectrum is equal to or less than a threshold or not; and a control part (171) causing the reception part to receive a second nuclear magnetic resonance signal from the specific area when the SN ratio is equal to or less than the threshold. The image processing part calculated a frequency spectrum, based on the first nuclear magnetic resonance signal and the second nuclear magnetic resonance signal.

Description

本発明は、磁気共鳴スペクトロスコピー装置に関する。特に代謝物質毎の分布を測定するスペクトログラムのSN比の向上に関する。   The present invention relates to a magnetic resonance spectroscopy apparatus. In particular, the present invention relates to an improvement in the SN ratio of a spectrogram for measuring the distribution of each metabolite.

近年、磁気共鳴現象を利用することで、MRI(Magnetic Resonance Imaging)画像では得られない代謝物質の濃度分布をスペクトルにより視覚的に捉えることができる磁気共鳴スペクトロスコピー(MRS:Magnetic Resonance Spectroscopy)装置が開発されている。MRSとは、様々な分子の化学結合の違いによる磁気共鳴周波数の差異(ケミカルシフト)に基づいて、分子ごとの信号を分離する方法である。ケミカルシフトは、ある分子中の原子核の磁場が付近に存在する原子核によって作られる磁場の影響を受け、それにより同じ原子核であっても周囲の環境により共鳴周波数の違いが生じることで発生する。   2. Description of the Related Art In recent years, a magnetic resonance spectroscopy (MRS) apparatus that can visually grasp the concentration distribution of a metabolite that cannot be obtained with an MRI (Magnetic Resonance Imaging) image by using a spectrum by utilizing a magnetic resonance phenomenon has been developed. Has been developed. MRS is a method of separating signals for each molecule based on a difference (chemical shift) in magnetic resonance frequency due to a difference in chemical bonds of various molecules. A chemical shift occurs when the magnetic field of a nucleus in a molecule is affected by the magnetic field created by nearby nuclei, thereby causing a difference in resonance frequency depending on the surrounding environment even in the same nucleus.

MRSにおいて、病態の状況を把握するためには、代謝物質の相対量的情報が重要となる。正確な情報を得るためには撮影回数を増やせばよいが、臨床の現場においては、限られた検査時間に、最良の撮影を行わなくてはならない。現状のMRSでは、撮影が終了してから代謝物質の情報を確認するが、その時点でSN(Signal Noize)比が低い場合でも、代謝物質の量的情報を検討しなくてはならない。   In MRS, relative quantitative information of metabolites is important for grasping the state of a disease state. In order to obtain accurate information, the number of imaging may be increased. However, in clinical practice, the best imaging must be performed within a limited examination time. In the current MRS, information on metabolites is confirmed after imaging is completed, but quantitative information on metabolites must be examined even if the SN (Signal Noize) ratio is low at that time.

特許文献1の装置では、エッジフィルタをSN比が低いスペクトログラムに対して処理することで、高分解能のスペクトログラムを得ている。   In the apparatus of Patent Document 1, a high-resolution spectrogram is obtained by processing an edge filter on a spectrogram having a low S / N ratio.

特開2007−301118号公報JP 2007-301118 A

しかし、エッジフィルタで処理しても、結果としてスペクトログラムが不正確になることもあり、虚部信号によって代謝物質量が実際に存在しているものより多く見えてしまうこともある。そこで操作者又は医師に誤解を与えないスペクトログラム、つまり一定以上のSN比を有するスペクトログラムを得ることが必要である。   However, processing with an edge filter may result in an inaccurate spectrogram, and the imaginary part signal may make the amount of metabolite appear more than what is actually present. Therefore, it is necessary to obtain a spectrogram that does not give a misunderstanding to the operator or doctor, that is, a spectrogram having a certain S / N ratio or more.

そこで、本発明は、一定以上のSN比を有するスペクトログラムを得る磁気共鳴スペクトロスコピー装置を提供する。   Therefore, the present invention provides a magnetic resonance spectroscopy apparatus that obtains a spectrogram having a certain S / N ratio or more.

第1の観点の磁気共鳴スペクトロスコピー装置は、被検体の特定領域に高周波磁場を照射する送信部と、特定領域から第1の核磁気共鳴信号を受信する受信部と、受信部で受信した第1の核磁気共鳴信号を記憶する記憶部と、第1の核磁気共鳴信号に基づいて、周波数スペクトラムを算出する算出部とを備える。さらに、磁気共鳴スペクトロスコピー装置は、周波数スペクトラムの特定物質のSN比が閾値以下であるか否かを判定するSN比判定部と、SN比が閾値以下である場合に、送信部に再び高周波磁場を照射させ、受信部に特定領域から第2の核磁気共鳴信号を受信させる制御部と、を備える。そして画像処理部は、第1の核磁気共鳴信号及び第2の核磁気共鳴信号に基づいて、周波数スペクトラムを算出する。
つまり、SN比が低かった場合でも、さらに第1の核磁気共鳴信号に加えて、必要な受信回数だけ第2の核磁気共鳴信号を得て、周波数スペクトラムを算出する。
A magnetic resonance spectroscopy apparatus according to a first aspect includes a transmitter that irradiates a specific region of a subject with a high-frequency magnetic field, a receiver that receives a first nuclear magnetic resonance signal from the specific region, and a receiver that is received by the receiver. A storage unit that stores one nuclear magnetic resonance signal; and a calculation unit that calculates a frequency spectrum based on the first nuclear magnetic resonance signal. Further, the magnetic resonance spectroscopy apparatus includes an SN ratio determination unit that determines whether or not the SN ratio of a specific substance in the frequency spectrum is equal to or less than a threshold value, and when the SN ratio is equal to or less than the threshold value, the transmission unit again receives a high frequency magnetic field. And a control unit that causes the receiving unit to receive the second nuclear magnetic resonance signal from the specific region. The image processing unit calculates a frequency spectrum based on the first nuclear magnetic resonance signal and the second nuclear magnetic resonance signal.
That is, even when the SN ratio is low, in addition to the first nuclear magnetic resonance signal, the second nuclear magnetic resonance signal is obtained as many times as necessary and the frequency spectrum is calculated.

第2の観点の磁気共鳴スペクトロスコピー装置の制御部は、第1の核磁気共鳴信号の受信回数とSN比とに基づいて、閾値に達する不足の受信回数を算出し、受信部に第2の核磁気共鳴信号を不足の受信回数分だけ受信させる。
第3の観点の磁気共鳴スペクトロスコピー装置の制御部は、第2の核磁気共鳴信号の撮影時間に基づいて、受信回数を算出する。
The control unit of the magnetic resonance spectroscopy apparatus according to the second aspect calculates the insufficient number of receptions reaching the threshold based on the number of receptions of the first nuclear magnetic resonance signal and the SN ratio, The nuclear magnetic resonance signal is received for the number of times of insufficient reception.
The control unit of the magnetic resonance spectroscopy apparatus of the third aspect calculates the number of receptions based on the imaging time of the second nuclear magnetic resonance signal.

第4の観点の磁気共鳴スペクトロスコピー装置は、SN比が閾値を入力する入力部をさらに備える。
第5の観点の磁気共鳴スペクトロスコピー装置は、SN比を計算するSN比計算部を備え、SN比計算部は、周波数スペクトラムを信号領域とノイズ領域とに判別分析し、特定物質のピークとノイズ成分の平均値とに基づいてSN比を計算する。
The magnetic resonance spectroscopy apparatus according to the fourth aspect further includes an input unit for inputting a threshold value for the SN ratio.
A magnetic resonance spectroscopy apparatus according to a fifth aspect includes an S / N ratio calculation unit that calculates an S / N ratio, and the S / N ratio calculation unit discriminates and analyzes a frequency spectrum into a signal region and a noise region, and the peak and noise of a specific substance. The SN ratio is calculated based on the average value of the components.

第6の観点の磁気共鳴スペクトロスコピー装置の画像処理部は、複数のボクセル毎の周波数スペクトラムを算出する際には、複数のボクセルの少なくとも1つが基準ボクセルとして設定され、SN比判定部は、基準ボクセルにおいて、SN比が閾値以下であるか否かを判定する。
第7の観点の磁気共鳴スペクトロスコピー装置は、基準ボクセルが二以上の場合には、SN比判定部は、基準ボクセルにおいて、SN比を平均化する。
In the image processing unit of the magnetic resonance spectroscopy apparatus of the sixth aspect, when calculating the frequency spectrum for each of the plurality of voxels, at least one of the plurality of voxels is set as a reference voxel, and the SN ratio determination unit In the voxel, it is determined whether or not the SN ratio is equal to or less than a threshold value.
In the magnetic resonance spectroscopy apparatus according to the seventh aspect, when the number of reference voxels is two or more, the SN ratio determination unit averages the SN ratios in the reference voxels.

本発明の磁気共鳴スペクトロスコピー装置は、一定以上のSN比を有するスペクトログラムを得ることができる。   The magnetic resonance spectroscopy apparatus of the present invention can obtain a spectrogram having an S / N ratio above a certain level.

磁気共鳴イメージング装置10の概略構成図である。1 is a schematic configuration diagram of a magnetic resonance imaging apparatus 10. FIG. 第1実施形態のスペクトラムの算出フローチャートである。It is a calculation flowchart of the spectrum of a 1st embodiment. 第1の核磁気共鳴信号に基づいたスペクトログラムSP1の一例である。It is an example of spectrogram SP1 based on the 1st nuclear magnetic resonance signal. スペクトログラムSP1の信号領域SAとノイズ領域NAとを示す図である。It is a figure which shows signal area SA and noise area NA of spectrogram SP1. 第1及び第2の核磁気共鳴信号に基づいたスペクトログラムSP2の一例である。It is an example of spectrogram SP2 based on the 1st and 2nd nuclear magnetic resonance signal. 第2実施形態のスペクトログラムの算出フローチャートである。It is a calculation flowchart of the spectrogram of the second embodiment. 第3実施形態のマルチボクセルに基づいた、設定とスペクトログラムとを示した図である。It is the figure which showed the setting and spectrogram based on the multi-voxel of 3rd Embodiment. 第3実施形態のスペクトログラムの生成フローチャートである。It is a production | generation flowchart of the spectrogram of 3rd Embodiment.

第1から第3実施形態の磁気共鳴スペクトロスコピー装置(MRS:Magnetic Resonance Spectroscopy)について説明する。また、以下では、本発明に係る磁気共鳴スペクトロスコピー装置の機能を搭載した磁気共鳴イメージング装置を実施形態として説明し、磁気共鳴イメージング装置を「MRI(Magnetic Resonance Imaging)装置」と呼ぶ。   The magnetic resonance spectroscopy apparatus (MRS: Magnetic Resonance Spectroscopy) of the first to third embodiments will be described. In the following, a magnetic resonance imaging apparatus equipped with the function of a magnetic resonance spectroscopy apparatus according to the present invention will be described as an embodiment, and the magnetic resonance imaging apparatus will be referred to as an “MRI (Magnetic Resonance Imaging) apparatus”.

<<第1実施形態>>
<磁気共鳴イメージング装置の構成>
第1実施形態に係るMRI装置の構成について説明する。図1は、第1実施形態に係るMRI装置の構成を説明するための図である。なお、第1実施形態に係るMRI装置は、MRI画像を撮影するとともに、関心領域における磁気共鳴スペクトルを収集して、複数の代謝物質の濃度分布強度を示すスペクトル画像を生成する装置である。図1を参照して、第1実施形態の磁気共鳴イメージング装置10の構成及びその基本動作について説明する。
<< First Embodiment >>
<Configuration of magnetic resonance imaging apparatus>
A configuration of the MRI apparatus according to the first embodiment will be described. FIG. 1 is a diagram for explaining the configuration of the MRI apparatus according to the first embodiment. The MRI apparatus according to the first embodiment is an apparatus that captures an MRI image and collects magnetic resonance spectra in a region of interest to generate a spectrum image indicating the concentration distribution intensities of a plurality of metabolites. With reference to FIG. 1, the configuration and basic operation of the magnetic resonance imaging apparatus 10 of the first embodiment will be described.

磁気共鳴イメージング装置10は、マグネットシステム100、勾配コイル駆動部130、RFコイル駆動部140、データ収集部150、パルスシーケンス制御部160、演算部170、表示部180及び操作部190を有する。   The magnetic resonance imaging apparatus 10 includes a magnet system 100, a gradient coil drive unit 130, an RF coil drive unit 140, a data collection unit 150, a pulse sequence control unit 160, a calculation unit 170, a display unit 180, and an operation unit 190.

マグネットシステム100は、主磁場コイル部102、勾配コイル部106及びRFコイル部108を有している。これら各コイル部は概ね円筒状の形状を有し、概ね円柱状のボアに互いに同軸状に配置されている。ボア内には被検者SBが寝台110に載置されており、寝台110は、撮影部位に応じて、マグネットシステム100内のボア内を移動可能になっている。   The magnet system 100 includes a main magnetic field coil unit 102, a gradient coil unit 106, and an RF coil unit 108. Each of these coil portions has a substantially cylindrical shape, and is arranged coaxially with each other in a substantially columnar bore. The subject SB is placed on the bed 110 in the bore, and the bed 110 is movable in the bore in the magnet system 100 according to the imaging region.

主磁場コイル部102は、マグネットシステム100の内部空間に静磁場を形成する。静磁場の方向は、概ね被検者SBの体軸の方向に平行であり水平磁場を形成する。主磁場コイル部102は、通常、超伝導コイルを用いて構成されるが、超伝導コイルに限らず永久磁石等を用いて構成してもよい。   The main magnetic field coil unit 102 forms a static magnetic field in the internal space of the magnet system 100. The direction of the static magnetic field is generally parallel to the direction of the body axis of the subject SB and forms a horizontal magnetic field. The main magnetic field coil unit 102 is normally configured using a superconducting coil, but may be configured using a permanent magnet or the like without being limited to the superconducting coil.

勾配コイル部106は、互いに直交する3軸、すなわち、スライス軸、位相軸及び周波数軸の方向において、それぞれ主磁場コイル部102によって形成された静磁場強度に勾配を持たせるための3種の勾配磁場を発生する。このような勾配磁場の発生を可能にするために、勾配コイル部106は、図示しない3系統の勾配コイルを有する。勾配コイル部106には勾配コイル駆動部130が接続されており、勾配コイル駆動部130は勾配コイル部106に駆動信号を与えて勾配磁場を発生させる。勾配コイル駆動部130は、勾配コイル部106における3系統の勾配コイルに対応して、図示しない3系統の駆動回路を有する。   The gradient coil unit 106 has three types of gradients for imparting gradients to the static magnetic field strength formed by the main magnetic field coil unit 102 in the three axes orthogonal to each other, that is, in the direction of the slice axis, the phase axis, and the frequency axis. Generate a magnetic field. In order to make it possible to generate such a gradient magnetic field, the gradient coil unit 106 has three gradient coils (not shown). A gradient coil drive unit 130 is connected to the gradient coil unit 106, and the gradient coil drive unit 130 gives a drive signal to the gradient coil unit 106 to generate a gradient magnetic field. The gradient coil drive unit 130 has three drive circuits (not shown) corresponding to the three gradient coils in the gradient coil unit 106.

RFコイル部108は、静磁場空間に被検者SBの体内のスピンを励起するための高周波磁場を形成する。高周波磁場を形成することをRF励起信号の送信といい、RF励起信号をRFパルスという。RFコイル部108にはRFコイル駆動部140が接続されており、RFコイル駆動部140はRFコイル部108に駆動信号を与え、その駆動信号に基づいてRFコイル部108はRFパルスを送信する。励起されたスピンが生じる電磁波すなわち核磁気共鳴信号は、RFコイル部108によって受信される。RFコイル部108にはデータ収集部150が接続されている。データ収集部150は、RFコイル部108が受信した核磁気共鳴信号をデジタルデータとして収集する。   The RF coil unit 108 forms a high-frequency magnetic field for exciting spins in the body of the subject SB in the static magnetic field space. Formation of a high-frequency magnetic field is called transmission of an RF excitation signal, and the RF excitation signal is called an RF pulse. An RF coil drive unit 140 is connected to the RF coil unit 108. The RF coil drive unit 140 provides a drive signal to the RF coil unit 108, and the RF coil unit 108 transmits an RF pulse based on the drive signal. An electromagnetic wave generated by the excited spin, that is, a nuclear magnetic resonance signal is received by the RF coil unit 108. A data collection unit 150 is connected to the RF coil unit 108. The data collection unit 150 collects the nuclear magnetic resonance signals received by the RF coil unit 108 as digital data.

勾配コイル駆動部130、RFコイル駆動部140及びデータ収集部150にはパルスシーケンス制御部160が接続されている。   A pulse sequence control unit 160 is connected to the gradient coil drive unit 130, the RF coil drive unit 140, and the data collection unit 150.

パルスシーケンス制御部160は、操作者が入力した撮影条件、すなわち撮影プロトコルに従い、勾配コイル駆動部130及びRFコイル駆動部140を駆動させる。より具体的には、パルスシーケンス制御部160は、「MRI画像を撮影するためのシーケンス情報」および「スペクトログラムを撮影するためのシーケンス情報」を生成する。それらシーケンス情報は、勾配コイル駆動部130及びRFコイル駆動部140に送信される。   The pulse sequence control unit 160 drives the gradient coil driving unit 130 and the RF coil driving unit 140 in accordance with the imaging conditions input by the operator, that is, the imaging protocol. More specifically, the pulse sequence control unit 160 generates “sequence information for photographing an MRI image” and “sequence information for photographing a spectrogram”. The sequence information is transmitted to the gradient coil driving unit 130 and the RF coil driving unit 140.

表示部180は、グラフィックディスプレー等で構成されている。表示部180は演算部170に接続されている。表示部180は、GUIの操作画面、MRI画像再構成用の核磁気共鳴信号に基づいて画像再構成された磁気共鳴画像、及びスペクトログラム用の核磁気共鳴信号に基づいてスペクトログラムなどを表示することができる。   The display unit 180 is configured by a graphic display or the like. The display unit 180 is connected to the calculation unit 170. The display unit 180 may display a GUI operation screen, a magnetic resonance image reconstructed based on a nuclear magnetic resonance signal for MRI image reconstruction, a spectrogram based on a nuclear magnetic resonance signal for a spectrogram, and the like. it can.

操作部190は、ポインティングデバイスを有するキーボード等で構成される。操作部190は演算部170に接続されている。操作部190は、操作者によって表示部180を介して操作される。操作部190は、キーボード等の代わりに表示部180にタッチパネルを配置してもよい。   The operation unit 190 includes a keyboard having a pointing device. The operation unit 190 is connected to the calculation unit 170. The operation unit 190 is operated by the operator via the display unit 180. The operation unit 190 may arrange a touch panel on the display unit 180 instead of a keyboard or the like.

演算部170は、制御部171、記憶部172、画像処理部174、SN比計算部176及びSN比判定部178を有する。
演算部170は、各種データの処理及びプログラムを実行する。制御部171は、データ収集部150及びパルスシーケンス制御部160等を制御する。記憶部172は、各種撮影プロトコル、各種プログラム及び各種データを記憶する。また記憶部172は、データ収集部150が収集した核磁気共鳴信号を記憶する。
The calculation unit 170 includes a control unit 171, a storage unit 172, an image processing unit 174, an SN ratio calculation unit 176, and an SN ratio determination unit 178.
The arithmetic unit 170 executes various data processes and programs. The control unit 171 controls the data collection unit 150, the pulse sequence control unit 160, and the like. The storage unit 172 stores various shooting protocols, various programs, and various data. The storage unit 172 stores the nuclear magnetic resonance signals collected by the data collection unit 150.

画像処理部174は、MRI画像再構成用の核磁気共鳴信号に基づいて磁気共鳴画像を画像再構成する。また、画像処理部174はスペクトログラム用の核磁気共鳴信号をフーリエ変換して周波数空間のデータに変換する。画像処理部174はその周波数空間のデータを周波数に対してプロットしてスペクトログラムを生成する。   The image processing unit 174 reconstructs a magnetic resonance image based on the nuclear magnetic resonance signal for MRI image reconstruction. Also, the image processing unit 174 performs Fourier transform on the spectrogram nuclear magnetic resonance signal to convert it into frequency space data. The image processing unit 174 plots the data in the frequency space against the frequency to generate a spectrogram.

SN比計算部176は、このスペクトログラムをまず信号領域とノイズ領域とに二分化する。例えば判別分析法を適用して、スペクトログラムの信号値が、基準値より大きい領域と、信号値が基準値より小さい領域とに、スペクトログラムの周波数の領域を二分化する。信号値が基準値以上の領域を信号領域、信号値が基準値以下の領域をノイズ領域とする。そして、SN比計算部176は、信号領域のピークとノイズ領域の平均値とのSN比を計算する。   The SN ratio calculation unit 176 first divides this spectrogram into a signal region and a noise region. For example, by applying a discriminant analysis method, the spectrogram frequency region is divided into a region where the spectrogram signal value is larger than the reference value and a region where the signal value is smaller than the reference value. An area where the signal value is equal to or greater than the reference value is defined as a signal area, and an area where the signal value is equal to or less than the reference value is defined as a noise area. Then, the SN ratio calculation unit 176 calculates the SN ratio between the peak of the signal region and the average value of the noise region.

SN比判定部178は、SN比と閾値とを比較し、SN比が閾値を越えているか否かを判定する。   The SN ratio determination unit 178 compares the SN ratio with a threshold value, and determines whether or not the SN ratio exceeds the threshold value.

<スペクトログラムの生成工程>
図2は、第1実施形態のスペクトログラムの生成フローチャートである。
ステップS11では、操作者は、スペクトログラムに要求するSN比の閾値を入力する。大脳のスペクトログラムを生成する際には、代謝物質として、例えばクロアチン(Creatine)、コリン(Cholin)又はNAA(N-Aceryl-L―aspartate)が選ばれる。部位によっては、代謝物質として、乳酸又はクエン酸等が選ばれることもある。第1実施形態では、操作者は、操作部190を介して、例えば、代謝物質としてクロアチンのピークとノイズ領域の平均値とのSN比の閾値を“8”と入力する。
<Spectrogram generation process>
FIG. 2 is a spectrogram generation flowchart of the first embodiment.
In step S11, the operator inputs a threshold value of the SN ratio required for the spectrogram. When generating a cerebral spectrogram, for example, creatine, choline, or NAA (N-Aceryl-L-aspartate) is selected as a metabolite. Depending on the site, lactic acid or citric acid may be selected as a metabolite. In the first embodiment, the operator inputs, for example, “8” as the SN ratio threshold between the peak of croatin and the average value of the noise region as a metabolite via the operation unit 190.

なお、磁気共鳴イメージング装置10の記憶部172が、これまでの経験値等からデフォルトとして閾値を記憶している場合には、ステップS11を省略してもよい。また、第1実施形態のSN比は、代謝物質のピークの高さとノイズ領域の平均値とのピーク信号対雑音比(PSNR)で説明する。しかし、本発明はこれに限られず、代謝物質の平均値とノイズ領域の平均値とのSN比を適用してもよい。   Note that when the storage unit 172 of the magnetic resonance imaging apparatus 10 stores a threshold value as a default from the experience values so far, step S11 may be omitted. The SN ratio of the first embodiment will be described as a peak signal-to-noise ratio (PSNR) between the peak height of the metabolite and the average value of the noise region. However, the present invention is not limited to this, and an SN ratio between the average value of the metabolite and the average value of the noise region may be applied.

ステップS13では、まず表示部180に撮影条件設定用の初期画面が表示される。操作者は、表示された被検者SBの大脳内に撮影領域を設定し、ボクセルサイズを設定する。さらに操作者は、データ収集部150が収集するスペクトログラム用の核磁気共鳴信号の受信回数を設定する。第1実施形態では、操作者は、1箇所のボクセル(シングルボクセル)を設定し、操作部190を介して、例えば受信回数を“30”と入力する。   In step S13, an initial screen for setting shooting conditions is first displayed on the display unit 180. The operator sets an imaging region in the cerebrum of the displayed subject SB and sets a voxel size. Further, the operator sets the number of receptions of spectrogram nuclear magnetic resonance signals collected by the data collection unit 150. In the first embodiment, the operator sets one voxel (single voxel), and inputs, for example, “30” as the number of receptions via the operation unit 190.

ステップS15では、演算部170の指令に応じて、パルスシーケンス制御部160が、撮影プロトコルに従い、スペクトログラムを撮影するためのパルスシーケンスを送信する。   In step S15, the pulse sequence control unit 160 transmits a pulse sequence for imaging the spectrogram according to the imaging protocol in accordance with a command from the arithmetic unit 170.

ステップS17では、データ収集部150が、RFコイル部108が受信した第1の核磁気共鳴信号を収集する。収集された第1の核磁気共鳴信号は、記憶部172に記憶される。記憶部172に記憶された第1の核磁気共鳴信号は、少なくとも本フローチャートの終了ステップに至るまでは記憶されている。   In step S17, the data collection unit 150 collects the first nuclear magnetic resonance signal received by the RF coil unit. The collected first nuclear magnetic resonance signal is stored in the storage unit 172. The first nuclear magnetic resonance signal stored in the storage unit 172 is stored at least until the end step of this flowchart.

ステップS19では、画像処理部174が、記憶された第1の核磁気共鳴信号に基づいて、スペクトログラムを生成する。詳述すると画像処理部174は、第1の核磁気共鳴信号をフーリエ変換し、周波数空間のデータに変換する。画像処理部174は、その周波数空間のデータを周波数に対してプロットしてスペクトログラムSP1を生成する。例えば図3は、その一例である。図3の横軸は基準物質からの周波数のずれの割合(ケミカルシフト量)を示し、縦軸はその周波数の波の波高、すなわち水素原子の量を示している。   In step S19, the image processing unit 174 generates a spectrogram based on the stored first nuclear magnetic resonance signal. More specifically, the image processing unit 174 performs Fourier transform on the first nuclear magnetic resonance signal to convert it into frequency space data. The image processing unit 174 plots the data in the frequency space against the frequency to generate a spectrogram SP1. For example, FIG. 3 shows an example. The horizontal axis of FIG. 3 shows the ratio of frequency deviation from the reference substance (chemical shift amount), and the vertical axis shows the wave height of that frequency, that is, the amount of hydrogen atoms.

ステップS21では、SN比計算部176が、代謝物質としてクロアチンのピークとノイズ領域の平均値とのSN比を計算する。詳細は図4を参照しながら説明する。SN比計算部176は、図4に表示されたスペクトログラムSP1を、例えば判別分析法を適用して、まず信号領域とノイズ領域とに二分化する。スペクトログラムSP1の信号値が、基準値STより大きい領域と、信号値が基準値STより小さい領域とに、スペクトログラムの周波数の領域を二分化する。信号値が基準値以上の領域を信号領域SA、信号値が基準値以下の領域をノイズ領域NA(網掛け領域)とする。   In step S21, the SN ratio calculation unit 176 calculates the SN ratio between the peak of croatin as a metabolite and the average value of the noise region. Details will be described with reference to FIG. The SN ratio calculation unit 176 first divides the spectrogram SP1 displayed in FIG. 4 into a signal region and a noise region, for example, by applying a discriminant analysis method. The region of the spectrogram frequency is divided into a region where the signal value of the spectrogram SP1 is larger than the reference value ST and a region where the signal value is smaller than the reference value ST. An area where the signal value is equal to or greater than the reference value is defined as a signal area SA, and an area where the signal value is equal to or less than the reference value is defined as a noise area NA (shaded area).

さらに、SN比計算部176は、ノイズ領域NAの平均値を計算し、ノイズ領域NAの平均値との信号領域SAのピークとSN比を計算する。例えば第1実施形態では、代謝物質がケミカルシフト量3.02ppmのクロアチンであり、クロアチンのピークの波高PK1とノイズ領域NAの平均値NV1とのSN比を計算する。代謝物質がクロアチンと決まっていれば、すべてのノイズ領域NAの平均値ではなく、クロアチンの周辺のノイズ領域NAの平均値を使用してもよい。第1実施形態では、SN比が“6”であったと仮定する。   Furthermore, the SN ratio calculation unit 176 calculates the average value of the noise area NA, and calculates the peak of the signal area SA and the SN ratio with respect to the average value of the noise area NA. For example, in the first embodiment, the metabolite is croatin having a chemical shift amount of 3.02 ppm, and the SN ratio between the peak height PK1 of croatin and the average value NV1 of the noise region NA is calculated. If the metabolite is determined to be croatin, the average value of noise regions NA around croatin may be used instead of the average value of all noise regions NA. In the first embodiment, it is assumed that the SN ratio is “6”.

ステップS23では、SN比判定部178がステップS11で設定された閾値よりSN比が大きいかを判定する。SN比が大きいということは、操作者又は医師が充分にそのスペクトログラムSP1を使って代謝物質を判断することができることを意味する。そのため、スペクトログラムSP1の生成はここで終了する。SN比が閾値を越えない場合には、ステップS25に進む。第1実施形態では、ステップS11で設定された閾値が“8”であり、SN比が“6”であるため、ステップS25に進む。   In step S23, the SN ratio determination unit 178 determines whether the SN ratio is larger than the threshold set in step S11. A high SN ratio means that the operator or doctor can sufficiently determine the metabolite using the spectrogram SP1. Therefore, the generation of spectrogram SP1 ends here. If the SN ratio does not exceed the threshold value, the process proceeds to step S25. In the first embodiment, since the threshold set in step S11 is “8” and the SN ratio is “6”, the process proceeds to step S25.

ステップS25では、SN比計算部176が、さらに核磁気共鳴振動を受信する回数を計算する。SN比が“6”のスペクトログラムSP1は、受信回数が30回のパルスシーケンスで生成された。このため、SN比計算部176は、30(受信回数)×8(閾値)/6(SN比)−30(受信回数)=10(不足分の受信回数)と計算する。つまり、10回以上、スペクトログラム用の核磁気共鳴信号を受信すれば、SN比の高いスペクトログラムを得ることができることがわかる。   In step S25, the SN ratio calculation unit 176 further calculates the number of times of receiving nuclear magnetic resonance vibration. The spectrogram SP1 with an S / N ratio of “6” was generated with a pulse sequence of 30 receptions. Therefore, the SN ratio calculation unit 176 calculates 30 (number of receptions) × 8 (threshold) / 6 (SN ratio) −30 (number of receptions) = 10 (number of receptions for the shortage). That is, it is understood that a spectrogram having a high S / N ratio can be obtained by receiving a nuclear magnetic resonance signal for spectrogram ten times or more.

ステップS27では、パルスシーケンス制御部160が、10回以上の核磁気共鳴信号を受信する撮影プロトコルに従い、パルスシーケンスを送信する。
ステップS29では、データ収集部150が、RFコイル部108が受信した第2の核磁気共鳴信号を収集する。収集された第2の核磁気共鳴信号は、記憶部172に記憶される。
In step S27, the pulse sequence control unit 160 transmits a pulse sequence according to an imaging protocol for receiving 10 or more nuclear magnetic resonance signals.
In step S29, the data collection unit 150 collects the second nuclear magnetic resonance signal received by the RF coil unit 108. The collected second nuclear magnetic resonance signal is stored in the storage unit 172.

ステップS31では、画像処理部174が、記憶された第1の核磁気共鳴信号及び第2の核磁気共鳴信号に基づいて、スペクトログラムSP2を生成する。例えば図5は、第1の核磁気共鳴信号及び第2の核磁気共鳴信号に基づいて生成されたスペクトログラムSP2の一例である。図5では、スペクトログラムSP2がスペクトログラムSP1(点線)と併記して描いてある。スペクトログラムSP2は、SN比が大きいため、操作者又は医師が充分にそのスペクトログラムSP2を使って代謝物質を判断することができる。   In step S31, the image processing unit 174 generates a spectrogram SP2 based on the stored first nuclear magnetic resonance signal and second nuclear magnetic resonance signal. For example, FIG. 5 is an example of a spectrogram SP2 generated based on the first nuclear magnetic resonance signal and the second nuclear magnetic resonance signal. In FIG. 5, the spectrogram SP2 is drawn together with the spectrogram SP1 (dotted line). Since the spectrogram SP2 has a large S / N ratio, the operator or doctor can sufficiently determine the metabolite using the spectrogram SP2.

<<第2実施形態>>
<スペクトログラムの生成工程>
図6は、第2実施形態のスペクトログラムの生成フローチャートである。第2実施形態では、第1実施形態のフローチャートに撮影時間の条件が加えられている。第1実施形態の生成フローチャートとは、ステップS11、S13、S25及びS27が異なっている。このため、新たなステップS11R、S13R、S25R及びS27Rのみを説明する。磁気共鳴イメージング装置の構成は第1実施形態と同じである。
<< Second Embodiment >>
<Spectrogram generation process>
FIG. 6 is a flowchart of spectrogram generation according to the second embodiment. In the second embodiment, a shooting time condition is added to the flowchart of the first embodiment. Steps S11, S13, S25, and S27 are different from the generation flowchart of the first embodiment. For this reason, only new steps S11R, S13R, S25R and S27R will be described. The configuration of the magnetic resonance imaging apparatus is the same as that of the first embodiment.

ステップS11Rでは、操作者は、スペクトログラムに要求するSN比の閾値を入力する。第2実施形態でも、操作者は、代謝物質としてクロアチンのピークとノイズ領域の平均値とのSN比の閾値を“8”と入力する。さらに、操作者は、上限の撮影時間を入力する。例えば操作者は、被検者SBの体調などを考慮して、撮影時間の上限を5分と設定する。   In step S11R, the operator inputs a threshold value of the S / N ratio required for the spectrogram. Also in the second embodiment, the operator inputs “8” as the threshold value of the SN ratio between the peak of croatin and the average value of the noise region as a metabolite. Furthermore, the operator inputs an upper limit shooting time. For example, the operator sets the upper limit of the imaging time to 5 minutes in consideration of the physical condition of the subject SB.

ステップS13Rでは、操作者は、表示された被検者SBの大脳内に撮影領域を1箇所に設定し、ボクセルサイズ及びスペクトログラム用の核磁気共鳴信号の受信回数を設定する。これにより、演算部170は、パルスシーケンスの送信(撮影開始)から表示部180にスペクトログラムを表示する(撮影終了)までの時間を算出する。例えば、第2実施形態では、操作者は、操作部190を介して、例えば受信回数を“30”と入力したため、撮影時間が4分と算出する。   In step S13R, the operator sets one imaging region in the cerebrum of the displayed subject SB, and sets the voxel size and the number of reception times of the nuclear magnetic resonance signal for the spectrogram. Thereby, the arithmetic unit 170 calculates the time from the transmission of the pulse sequence (start of imaging) to the display of the spectrogram on the display unit 180 (end of imaging). For example, in the second embodiment, since the operator inputs, for example, “30” as the number of receptions via the operation unit 190, the shooting time is calculated as 4 minutes.

ステップS25Rでは、SN比計算部176が、さらに核磁気共鳴振動を受信する回数を計算する。SN比が“6”のスペクトログラムSP1は、受信回数が30回のパルスシーケンスで生成された。このため、SN比計算部176は、30(受信回数)×8(閾値)/6(SN比)−30(受信回数)=10(不足分の受信回数)と計算する。つまり、10回以上、スペクトログラム用の核磁気共鳴信号を受信すれば、SN比の高いスペクトログラムを得ることができることがわかる。   In step S25R, the SN ratio calculation unit 176 further calculates the number of times of receiving nuclear magnetic resonance vibration. The spectrogram SP1 with an S / N ratio of “6” was generated with a pulse sequence of 30 receptions. Therefore, the SN ratio calculation unit 176 calculates 30 (number of receptions) × 8 (threshold) / 6 (SN ratio) −30 (number of receptions) = 10 (number of receptions for the shortage). That is, it is understood that a spectrogram having a high S / N ratio can be obtained by receiving a nuclear magnetic resonance signal for spectrogram ten times or more.

一方、演算部170は、スペクトログラム用の核磁気共鳴信号を10回分受信する撮影時間を算出する。演算部170は、例えば、パルスシーケンスの送信(撮影開始)から表示部180にスペクトログラムを表示する(撮影終了)までに2分かかると算出する。すでに、ステップS15からステップS23までに撮影時間が4分経過しているため、スペクトログラム用の核磁気共鳴信号を10回分受信すると合計撮影時間が6分になってしまう。これでは、ステップS11Rで設定した上限の撮影時間が5分を越えてしまう。そこで、演算部170は、残りの撮影時間1分で得られるパルスシーケンスをパルスシーケンス制御部160に指示する。例えば、スペクトログラム用の核磁気共鳴信号を5回分受信できるパルスシーケンスになる。   On the other hand, the calculation unit 170 calculates an imaging time for receiving the spectrogram nuclear magnetic resonance signal for 10 times. For example, the calculation unit 170 calculates that it takes 2 minutes from the transmission of the pulse sequence (start of imaging) to the display of the spectrogram on the display unit 180 (end of imaging). Since the imaging time has already passed from step S15 to step S23, if the nuclear magnetic resonance signal for the spectrogram is received 10 times, the total imaging time becomes 6 minutes. In this case, the upper limit photographing time set in step S11R exceeds 5 minutes. Therefore, the calculation unit 170 instructs the pulse sequence control unit 160 to obtain a pulse sequence that can be obtained in the remaining imaging time of 1 minute. For example, the pulse sequence can receive five nuclear magnetic resonance signals for the spectrogram.

ステップS27Rでは、パルスシーケンス制御部160が、5回の核磁気共鳴信号を受信する撮影プロトコルに従い、パルスシーケンスを送信する。   In step S27R, the pulse sequence control unit 160 transmits a pulse sequence according to an imaging protocol for receiving five nuclear magnetic resonance signals.

第2実施形態では、被検者SBの体調などの許す範囲で、充分にSN比が改善されたスペクトログラムSP2を使って代謝物質を操作者又は医師が判断することができる。   In the second embodiment, the operator or doctor can determine the metabolite using the spectrogram SP2 in which the SN ratio is sufficiently improved within the range allowed by the physical condition of the subject SB.

<<第3実施形態>>
第1実施形態及び第2実施形態では、シングルボクセルを設定した際に、複数回、核磁気共鳴信号を収集する場合について説明した。第3実施形態では、マルチボクセルを設定する例を説明する。
<< Third Embodiment >>
In the first embodiment and the second embodiment, the case where a nuclear magnetic resonance signal is collected a plurality of times when a single voxel is set has been described. In the third embodiment, an example in which multi-voxels are set will be described.

図7(a)は、操作者は、撮影条件設定用の初期画面にマルチボクセルの撮影領域を設定した状態を示した図である。図7(b)は、各ボクセルサイズのスペクトログラムを示した図である。   FIG. 7A is a diagram illustrating a state in which the operator sets a multi-voxel shooting area on the initial screen for setting shooting conditions. FIG. 7B shows a spectrogram of each voxel size.

<スペクトログラムの生成工程>
図8は、第3実施形態のスペクトログラムの生成フローチャートである。第1実施形態の生成フローチャートとは、ステップS13、及びS21が異なっている。このため、新たなステップS13T、及びS21Tのみを説明する。磁気共鳴イメージング装置の構成は第1実施形態と同じである。
<Spectrogram generation process>
FIG. 8 is a spectrogram generation flowchart of the third embodiment. Steps S13 and S21 are different from the generation flowchart of the first embodiment. Therefore, only new steps S13T and S21T will be described. The configuration of the magnetic resonance imaging apparatus is the same as that of the first embodiment.

ステップS13Tでは、まず表示部180に図7(a)に示されるような、撮影条件設定用の初期画面が表示される。操作者は、表示された被検者SBの大脳内に撮影領域を設定し、ボクセルサイズ及びボクセル数を設定する。図7(a)では9つのボクセルVX(VX1〜VX9)が設定されている。また、第1及び第2実施形態と同様に、操作者は、データ収集部150が収集するスペクトログラム用の核磁気共鳴信号の受信回数を設定する。   In step S13T, an initial screen for setting shooting conditions as shown in FIG. 7A is displayed on the display unit 180 first. The operator sets an imaging region in the cerebrum of the displayed subject SB, and sets the voxel size and the number of voxels. In FIG. 7A, nine voxels VX (VX1 to VX9) are set. Similarly to the first and second embodiments, the operator sets the number of receptions of the nuclear magnetic resonance signal for the spectrogram collected by the data collection unit 150.

さらに、操作者は、ボクセルVX(VX1〜VX9)から基準ボクセルSVXを1つ設定する。例えばボクセルVX5が基準ボクセルSVXと設定される。代謝物質のピークの高さとノイズ領域の平均値とのSN比を計算する場合、9つのボクセルVX(VX1〜VX9)があると、撮影領域の違い等から各々のSN比が異なる。このため、マルチボクセルの撮影領域の場合には、基準ボクセルSVXを1つ設定する。   Further, the operator sets one reference voxel SVX from the voxels VX (VX1 to VX9). For example, the voxel VX5 is set as the reference voxel SVX. When calculating the SN ratio between the peak height of the metabolite and the average value of the noise region, if there are nine voxels VX (VX1 to VX9), the respective SN ratios are different due to differences in imaging regions. Therefore, one reference voxel SVX is set in the case of a multi-voxel imaging region.

ステップS21Tでは、SN比計算部176が、基準ボクセルSVXであるボクセルVX5における代謝物質のピークとノイズ領域の平均値とのSN比を計算する。   In step S21T, the SN ratio calculation unit 176 calculates the SN ratio between the peak of the metabolite in the voxel VX5 which is the reference voxel SVX and the average value of the noise region.

なお、第3実施形態の応用例として、基準ボクセルSVXを2つ以上設定するように応用してもよい。例えばステップS13Tで、基準ボクセルSVXをボクセルVX5とボクセルVX6との2つに設定したい場合を考える。この場合には、ステップS21Tで、SN比計算部176が、ボクセルVX5における代謝物質のピークとノイズ領域の平均値とのSN比を計算し、且つボクセルVX6における代謝物質のピークとノイズ領域の平均値とのSN比を計算する。そしてSN比計算部176は、ボクセルVX5のSN比とボクセルVX6のSN比とを平均化する。SN比判定部178は、その平均化されたSN比が閾値より大きいかを判定すればよい。   Note that, as an application example of the third embodiment, it may be applied to set two or more reference voxels SVX. For example, consider a case in which it is desired to set the reference voxel SVX to two voxels VX5 and VX6 in step S13T. In this case, in step S21T, the SN ratio calculation unit 176 calculates the SN ratio between the peak of the metabolite in the voxel VX5 and the average value of the noise region, and the average of the peak of the metabolite in the voxel VX6 and the noise region. The signal-to-noise ratio with the value is calculated. Then, the SN ratio calculation unit 176 averages the SN ratio of the voxel VX5 and the SN ratio of the voxel VX6. The SN ratio determination unit 178 may determine whether the averaged SN ratio is larger than the threshold value.

10 … 磁気共鳴イメージング装置
100 … マグネットシステム
102 … 主磁場コイル部
106 … 勾配コイル部
108 … RFコイル部
110 … 寝台
130 … 勾配コイル駆動部
140 … RFコイル駆動部
150 … データ収集部
160 … パルスシーケンス制御部
170 … 演算部
171 … 制御部
172 … 記憶部
174 … 画像処理部
176 … SN比計算部
178 … SN比判定部
180 … 表示部
190 … 操作部
NA … ノイズ領域
NV1 … 平均値
PK1 … 波高
SA … 信号領域
SB … 被検者
SP1 … スペクトログラム
SP2 … スペクトログラム
ST … 基準値
SVX … 基準ボクセル
VX … ボクセル
DESCRIPTION OF SYMBOLS 10 ... Magnetic resonance imaging apparatus 100 ... Magnet system 102 ... Main magnetic field coil part 106 ... Gradient coil part 108 ... RF coil part 110 ... Bed 130 ... Gradient coil drive part 140 ... RF coil drive part 150 ... Data acquisition part 160 ... Pulse sequence Control unit 170 ... Calculation unit 171 ... Control unit 172 ... Storage unit 174 ... Image processing unit 176 ... SN ratio calculation unit 178 ... SN ratio determination unit 180 ... Display unit 190 ... Operation unit NA ... Noise region NV1 ... Average value PK1 ... Wave height SA ... Signal region SB ... Subject SP1 ... Spectrogram SP2 ... Spectrogram ST ... Reference value SVX ... Reference voxel VX ... Voxel

Claims (8)

被検体の特定領域に高周波磁場を照射する送信部と、
前記特定領域から第1の核磁気共鳴信号を受信する受信部と、
前記受信部で受信した前記第1の核磁気共鳴信号を記憶する記憶部と、
前記第1の核磁気共鳴信号に基づいて、周波数スペクトラムを算出する算出部と、
前記周波数スペクトラムの特定物質のSN比が閾値以下であるか否かを判定するSN比判定部と、
前記SN比が閾値以下である場合に、前記送信部に再び高周波磁場を照射させ、前記受信部に前記特定領域から第2の核磁気共鳴信号を受信させる制御部と、を備え、
前記画像処理部は、前記第1の核磁気共鳴信号及び前記第2の核磁気共鳴信号に基づいて、周波数スペクトラムを算出する磁気共鳴スペクトロスコピー装置。
A transmitter that irradiates a specific region of the subject with a high-frequency magnetic field;
A receiver for receiving a first nuclear magnetic resonance signal from the specific region;
A storage unit for storing the first nuclear magnetic resonance signal received by the reception unit;
A calculation unit for calculating a frequency spectrum based on the first nuclear magnetic resonance signal;
An SN ratio determination unit that determines whether an SN ratio of a specific substance in the frequency spectrum is equal to or less than a threshold;
A control unit that, when the SN ratio is equal to or less than a threshold value, causes the transmission unit to irradiate a high-frequency magnetic field again and causes the reception unit to receive a second nuclear magnetic resonance signal from the specific region;
The image processing unit is a magnetic resonance spectroscopy apparatus that calculates a frequency spectrum based on the first nuclear magnetic resonance signal and the second nuclear magnetic resonance signal.
前記制御部は、前記第1の核磁気共鳴信号の受信回数と前記SN比とに基づいて、前記閾値に達する不足の受信回数を算出し、前記受信部に前記第2の核磁気共鳴信号を不足の受信回数分だけ受信させる請求項1に記載の磁気共鳴スペクトロスコピー装置。   The control unit calculates an insufficient number of receptions that reach the threshold based on the number of receptions of the first nuclear magnetic resonance signal and the SN ratio, and sends the second nuclear magnetic resonance signal to the reception unit. The magnetic resonance spectroscopy apparatus according to claim 1, wherein reception is performed for an insufficient number of receptions. 前記制御部は、前記第2の核磁気共鳴信号の撮影時間に基づいて、前記受信回数を算出する請求項2に記載の磁気共鳴スペクトロスコピー装置。   The magnetic resonance spectroscopy apparatus according to claim 2, wherein the control unit calculates the number of receptions based on an imaging time of the second nuclear magnetic resonance signal. 前記前記SN比が閾値を入力する入力部をさらに備える請求項1から請求項3のいずれか一項に記載の磁気共鳴スペクトロスコピー装置。   The magnetic resonance spectroscopy apparatus according to claim 1, further comprising an input unit that inputs a threshold value for the SN ratio. 前記SN比を計算するSN比計算部を備え、
前記SN比計算部は、前記周波数スペクトラムを信号領域とノイズ領域とに判別分析し、前記特定物質のピークと前記ノイズ成分の平均値とに基づいて前記SN比を計算する請求項1から請求項4のいずれか一項に記載の磁気共鳴スペクトロスコピー装置。
An SN ratio calculation unit for calculating the SN ratio;
The SN ratio calculation unit discriminates and analyzes the frequency spectrum into a signal region and a noise region, and calculates the SN ratio based on a peak of the specific substance and an average value of the noise component. 5. The magnetic resonance spectroscopy apparatus according to any one of 4 above.
前記画像処理部が、複数のボクセル毎の前記周波数スペクトラムを算出する際には、前記複数のボクセルの少なくとも1つが基準ボクセルとして設定され、
前記SN比判定部は、前記基準ボクセルにおいて、前記SN比が閾値以下であるか否かを判定する請求項1から請求項5のいずれか一項に記載の磁気共鳴スペクトロスコピー装置。
When the image processing unit calculates the frequency spectrum for each of a plurality of voxels, at least one of the plurality of voxels is set as a reference voxel,
The magnetic resonance spectroscopy apparatus according to any one of claims 1 to 5, wherein the SN ratio determination unit determines whether the SN ratio is equal to or less than a threshold value in the reference voxel.
前記基準ボクセルが二以上の場合には、前記SN比判定部は、前記基準ボクセルにおいて、前記SN比を平均化する請求項6に記載の磁気共鳴スペクトロスコピー装置。   The magnetic resonance spectroscopy apparatus according to claim 6, wherein when the number of reference voxels is two or more, the SN ratio determination unit averages the SN ratio in the reference voxels. 前記磁気共鳴スペクトロスコピー装置において、前記算出された周波数スペクトラムをグラフにして表示する表示部を備える請求項1から請求項8のいずれかに記載の磁気共鳴スペクトロスコピー装置。
The magnetic resonance spectroscopy apparatus according to claim 1, further comprising a display unit that displays the calculated frequency spectrum in a graph.
JP2011257122A 2011-11-25 2011-11-25 Magnetic resonance spectroscopy equipment Active JP5908704B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2011257122A JP5908704B2 (en) 2011-11-25 2011-11-25 Magnetic resonance spectroscopy equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2011257122A JP5908704B2 (en) 2011-11-25 2011-11-25 Magnetic resonance spectroscopy equipment

Publications (2)

Publication Number Publication Date
JP2013111097A true JP2013111097A (en) 2013-06-10
JP5908704B2 JP5908704B2 (en) 2016-04-26

Family

ID=48707294

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2011257122A Active JP5908704B2 (en) 2011-11-25 2011-11-25 Magnetic resonance spectroscopy equipment

Country Status (1)

Country Link
JP (1) JP5908704B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101774888B1 (en) 2016-10-31 2017-09-06 가천대학교 산학협력단 Pre-treatment apparatus for magnetic resonance spectroscopic data and method of pre-treatment for magnetic resonance spectroscopic data

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05126927A (en) * 1991-11-08 1993-05-25 Jeol Ltd Nmr device with automatic scan switching function
JPH06235760A (en) * 1993-02-10 1994-08-23 Hitachi Ltd Nuclear magnetic resonator and automatically setting method for measurement condition of the resonator
JPH08191816A (en) * 1995-01-19 1996-07-30 Ge Yokogawa Medical Syst Ltd Mr spctroscopy and mris apparatus
JP2006235760A (en) * 2005-02-22 2006-09-07 Nippon Telegr & Teleph Corp <Ntt> Information browsing system, method, and program
JP2006259326A (en) * 2005-03-17 2006-09-28 Shimadzu Corp Infrared microscope
JP2009240767A (en) * 2008-03-10 2009-10-22 Toshiba Corp Magnetic resonance imaging apparatus
JP2010110567A (en) * 2008-11-10 2010-05-20 Hitachi Medical Corp Magnetic resonance imaging apparatus

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH05126927A (en) * 1991-11-08 1993-05-25 Jeol Ltd Nmr device with automatic scan switching function
JPH06235760A (en) * 1993-02-10 1994-08-23 Hitachi Ltd Nuclear magnetic resonator and automatically setting method for measurement condition of the resonator
JPH08191816A (en) * 1995-01-19 1996-07-30 Ge Yokogawa Medical Syst Ltd Mr spctroscopy and mris apparatus
JP2006235760A (en) * 2005-02-22 2006-09-07 Nippon Telegr & Teleph Corp <Ntt> Information browsing system, method, and program
JP2006259326A (en) * 2005-03-17 2006-09-28 Shimadzu Corp Infrared microscope
JP2009240767A (en) * 2008-03-10 2009-10-22 Toshiba Corp Magnetic resonance imaging apparatus
JP2010110567A (en) * 2008-11-10 2010-05-20 Hitachi Medical Corp Magnetic resonance imaging apparatus

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101774888B1 (en) 2016-10-31 2017-09-06 가천대학교 산학협력단 Pre-treatment apparatus for magnetic resonance spectroscopic data and method of pre-treatment for magnetic resonance spectroscopic data
WO2018080293A1 (en) * 2016-10-31 2018-05-03 가천대학교 산학협력단 Magnetic resonance spectroscopic signal pre-processing apparatus and magnetic resonance spectroscopic signal pre-processing method

Also Published As

Publication number Publication date
JP5908704B2 (en) 2016-04-26

Similar Documents

Publication Publication Date Title
JP5559223B2 (en) Magnetic resonance imaging system
JP6462306B2 (en) Medical image processing apparatus and magnetic resonance imaging apparatus
EP2506026A1 (en) Method of dynamically compensating for magnetic field heterogeneity in magnetic resonance imaging
JP6084392B2 (en) Magnetic resonance imaging system
JP5479427B2 (en) Magnetic resonance imaging device
JP2009050615A (en) Magnetic resonance imaging apparatus and magnetic resonance image displaying method
JP2014057861A (en) Method and apparatus for controlling magnetic resonance system
EP2592432A1 (en) Chemical-shift selective multislice magnetic resonance imaging using the gradient reversal technique
US9720061B2 (en) Systems, methods and GUI for chemical exchange saturation transfer (CEST) analysis
JP2008173175A (en) Mri apparatus
JP2016198392A (en) Magnetic resonance imaging apparatus
JP2009034152A (en) Mri apparatus
JP5908704B2 (en) Magnetic resonance spectroscopy equipment
EP2511696A1 (en) MRI of chemical species having different resonance frequencies using an ultra-short echo time sequence
JP5636058B2 (en) Magnetic resonance imaging device
JP4862069B2 (en) Magnetic resonance diagnostic equipment
JP6445297B2 (en) Magnetic resonance imaging apparatus and medical data analysis apparatus
JP7242472B2 (en) MAGNETIC RESONANCE IMAGING APPARATUS AND MAGNETIC RESONANCE IMAGING METHOD
US20150351659A1 (en) Magnetic resonanceimaging apparatus and measurement method thereof
EP1995604B1 (en) High speed, high resolution, silent, real-time MRI method
EP2898337A2 (en) Metabolite detection system and operation thereof
JP2013240537A (en) Mri apparatus
JP2010099539A (en) Magnetic resonance device
JP2017121463A (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP2004148024A (en) Quantification method for n-acetyl asparate, glutamine, and gultamate and magnetic resonance imaging apparatus

Legal Events

Date Code Title Description
A625 Written request for application examination (by other person)

Free format text: JAPANESE INTERMEDIATE CODE: A625

Effective date: 20140612

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20150323

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150420

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20150609

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20150616

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20150721

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20150908

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20151208

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: 20160301

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160324

R150 Certificate of patent or registration of utility model

Ref document number: 5908704

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250