JP6579908B2 - Magnetic resonance imaging apparatus and diffusion weighted image calculation method - Google Patents

Magnetic resonance imaging apparatus and diffusion weighted image calculation method Download PDF

Info

Publication number
JP6579908B2
JP6579908B2 JP2015206108A JP2015206108A JP6579908B2 JP 6579908 B2 JP6579908 B2 JP 6579908B2 JP 2015206108 A JP2015206108 A JP 2015206108A JP 2015206108 A JP2015206108 A JP 2015206108A JP 6579908 B2 JP6579908 B2 JP 6579908B2
Authority
JP
Japan
Prior art keywords
adc
diffusion coefficient
diffusion
imaging apparatus
resonance imaging
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
JP2015206108A
Other languages
Japanese (ja)
Other versions
JP2017077316A (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 Ltd
Original Assignee
Hitachi Ltd
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 Ltd filed Critical Hitachi Ltd
Priority to JP2015206108A priority Critical patent/JP6579908B2/en
Publication of JP2017077316A publication Critical patent/JP2017077316A/en
Application granted granted Critical
Publication of JP6579908B2 publication Critical patent/JP6579908B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)

Description

本発明は、磁気共鳴イメージング(以下、「MRI」という)装置に関し、特に、複数(2以上)のb値を用いてそれぞれ取得した複数の拡散強調画像を用いて、任意のb値の拡散強調画像を求める技術に関する。   The present invention relates to a magnetic resonance imaging (hereinafter referred to as `` MRI '') apparatus, and in particular, diffusion enhancement of an arbitrary b value using a plurality of diffusion weighted images respectively acquired using a plurality (two or more) of b values. The present invention relates to a technique for obtaining an image.

MRI装置は、被検体、特に人体の組織を構成する原子核スピンが発生するNMR信号を計測し、その頭部、腹部、四肢等の形態や機能を2次元的に或いは3次元的に画像化する装置である。撮像においては、NMR信号には、傾斜磁場によって異なる位相エンコードが付与されるとともに周波数エンコードされて、時系列データとして計測される。計測されたNMR信号は、2次元又は3次元フーリエ変換されることにより画像に再構成される。   The MRI device measures NMR signals generated by the spins of the subject, especially the tissues of the human body, and visualizes the form and function of the head, abdomen, limbs, etc. in two or three dimensions Device. In imaging, the NMR signal is given different phase encoding depending on the gradient magnetic field and is frequency-encoded and measured as time-series data. The measured NMR signal is reconstructed into an image by two-dimensional or three-dimensional Fourier transform.

MRI装置で得られた拡散強調画像を用いた任意のb値の拡散強調画像の計算方法として、特許文献1に開示されている方法がある。ここでb値とは、拡散強調シーケンスに付加されるMPG(Motion Probing Gradient)の強さを表すパラメータであり、b値が大きいほど拡散が強調された拡散強調画像となる。なお、b値=0とは、拡散強調画像撮像シーケンスにMPGが付加されない(つまり、MPGの振幅がゼロである)ことを意味する。   As a method for calculating a diffusion-weighted image having an arbitrary b value using a diffusion-weighted image obtained by an MRI apparatus, there is a method disclosed in Patent Document 1. Here, the b value is a parameter representing the strength of MPG (Motion Probing Gradient) added to the diffusion enhancement sequence, and the diffusion enhancement image with enhanced diffusion as the b value increases. Note that b value = 0 means that no MPG is added to the diffusion weighted image capturing sequence (that is, the amplitude of the MPG is zero).

特許文献1に記載の方法は、2以上のb値を用いてそれぞれ撮像領域を撮像することにより、2以上の元画像を得て、2以上の元画像から見かけの拡散係数(ADC)を算出し、見かけの拡散係数(ADC)に基づいて任意のb値の拡散強調画像を計算する。この方法によって、SNRの低下やアーチファクトの増加のない拡散強調画像を計算することができる。以下、誤解が生じない限り、見かけの拡散係数(ADC)を、単に拡散係数(ADC)と記載する。   The method described in Patent Document 1 obtains two or more original images by imaging each imaging region using b values of two or more, and calculates an apparent diffusion coefficient (ADC) from the two or more original images. Then, a diffusion-weighted image having an arbitrary b value is calculated based on the apparent diffusion coefficient (ADC). By this method, it is possible to calculate a diffusion weighted image without a decrease in SNR or an increase in artifacts. Hereinafter, an apparent diffusion coefficient (ADC) is simply referred to as a diffusion coefficient (ADC) unless misunderstanding occurs.

特開2013-255854号公報JP 2013-255854 A

しかしながら、特許文献1では、撮像するMRI装置の静磁場強度による計算値の違いに関して配慮がされていない。特許文献1の方法で計算された拡散強調画像では、撮像に使用するMRI装置の静磁場強度によって得られる元拡散強調画像の信号値が異なるため、計算された拡散係数(ADC)も静磁場強度に依存することになる。   However, in Patent Document 1, no consideration is given to the difference in the calculated values depending on the static magnetic field strength of the MRI apparatus to be imaged. In the diffusion-weighted image calculated by the method of Patent Document 1, since the signal value of the original diffusion-weighted image obtained by the static magnetic field strength of the MRI apparatus used for imaging differs, the calculated diffusion coefficient (ADC) is also the static magnetic field strength. Will depend on.

そこで本発明の目的は、静磁場強度に依存することなく、見かけの拡散係数(ADC)を算出し、任意のb値で拡散強調画像を計算することのできるMRI装置及び拡散強調画像計算方法を提供することである。   Therefore, an object of the present invention is to provide an MRI apparatus and a diffusion weighted image calculation method capable of calculating an apparent diffusion coefficient (ADC) and calculating a diffusion weighted image with an arbitrary b value without depending on the static magnetic field strength. Is to provide.

上記目的を達成するために、本発明のMRI装置及び拡散強調画像計算方法は、b値=0を含む複数のb値でそれぞれ取得した各拡散強調画像から計測拡散係数(ADCmrs)を算出する。そして、組織毎に該計測拡散係数(ADCmrs)を基準拡散係数(ADCrfr)を用いて補正して補正拡散係数(ADCcrt)を得、該補正拡散係数(ADCcrt)を用いて任意のb値の拡散強調画像を算出する。或いは、計測拡散係数(ADCmrs)の修正するための情報に基づいて計測拡散係数(ADCmrs)を補正して補正拡散係数(ADCcrt)を得、該補正拡散係数(ADCcrt)を用いて任意のb値の拡散強調画像を算出する。 In order to achieve the above object, the MRI apparatus and the diffusion weighted image calculation method of the present invention calculate a measurement diffusion coefficient (ADC mrs ) from each diffusion weighted image acquired with a plurality of b values including b value = 0. . Then, by correcting using the measured diffusion coefficient for each tissue (ADC mrs) the reference diffusion coefficient (ADC RFR) to give the correction diffusion coefficient (ADC crt), any by using the correction diffusion coefficient (ADC crt) b-value diffusion weighted image is calculated. Alternatively, to obtain a measurement diffusion coefficient information measurement diffusion coefficient based on the order to correct the (ADC mrs) (ADC mrs) correction diffusion coefficient corrected (ADC crt), using the correction diffusion coefficient (ADC crt) A diffusion weighted image having an arbitrary b value is calculated.

本発明のMRI装置及び拡散強調画像の計算方法によれば、静磁場強度に依存することなく、見かけの拡散係数(ADC)を算出し、任意のb値で拡散強調画像を計算することができる。   According to the MRI apparatus and the diffusion weighted image calculation method of the present invention, an apparent diffusion coefficient (ADC) can be calculated without depending on the static magnetic field strength, and a diffusion weighted image can be calculated with an arbitrary b value. .

本発明に係るMRI装置の全体構成の一例を示す図The figure which shows an example of the whole structure of the MRI apparatus which concerns on this invention 本発明の実施例1の処理フローを示すフローチャートThe flowchart which shows the processing flow of Example 1 of this invention. 本発明の実施例1に用いるGUIの図GUI diagram used in Example 1 of the present invention 本発明の実施例1の拡散係数を示す図The figure which shows the diffusion coefficient of Example 1 of this invention 本発明の実施例2の処理フローを示すフローチャートThe flowchart which shows the processing flow of Example 2 of this invention. 本発明の実施例2に用いるGUIの図GUI diagram used in Embodiment 2 of the present invention 本発明の実施例2の拡散係数を示す図The figure which shows the diffusion coefficient of Example 2 of this invention 本発明の実施例3の処理フローを示すフローチャートThe flowchart which shows the processing flow of Example 3 of this invention. 本発明の実施例3に用いるGUIの図GUI diagram used in Example 3 of the present invention 画素ごとの補正値を画像化して表示する例を示す図The figure which shows the example which images and displays the correction value for every pixel

以下、添付図面に従って本発明のMRI装置の好ましい実施形態について詳説する。なお、発明の実施形態を説明するための全図において、同一機能を有するものは同一符号を付け、その繰り返しの説明は省略する。   Hereinafter, preferred embodiments of the MRI apparatus of the present invention will be described in detail 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 of the invention, and the repetitive description thereof is omitted.

最初に、本発明に係るMRI装置の一例の全体概要を図1に基づいて説明する。図1は、本発明に係るMRI装置の一実施例の全体構成を示すブロック図である。このMRI装置は、NMR現象を利用して被検体の断層画像を得るもので、図1に示すように、MRI装置は静磁場発生系2と、傾斜磁場発生系3と、送信系5と、受信系6と、信号処理系7と、シーケンサ4と、演算処理部(CPU)8とを備えて構成される。   First, an overall outline of an example of an MRI apparatus according to the present invention will be described with reference to FIG. FIG. 1 is a block diagram showing the overall configuration of an embodiment of an MRI apparatus according to the present invention. This MRI apparatus uses a NMR phenomenon to obtain a tomographic image of a subject.As shown in FIG. 1, the MRI apparatus includes a static magnetic field generation system 2, a gradient magnetic field generation system 3, a transmission system 5, A reception system 6, a signal processing system 7, a sequencer 4, and an arithmetic processing unit (CPU) 8 are provided.

静磁場発生系2は、垂直磁場方式であれば被検体1の体軸と直交する方向に、水平磁場方式であれば被検体1の体軸方向に均一な静磁場を発生させるもので、被検体1の周りに永久磁石方式、常電導方式あるいは超電導方式の静磁場発生源が配置されている。   The static magnetic field generation system 2 generates a uniform static magnetic field in the direction perpendicular to the body axis of the subject 1 in the case of the vertical magnetic field method and in the direction of the body axis of the subject 1 in the case of the horizontal magnetic field method. A permanent magnet type, normal conduction type or superconducting type static magnetic field generation source is arranged around the specimen 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 that applies a gradient magnetic field 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 that drives each gradient magnetic field coil. Gradient magnetic fields Gx, Gy, and Gz are applied in the X, Y, and Z axis directions by driving the gradient magnetic field power supply 10 of each coil according to a command from the sequencer 4 described later. . At the time of imaging, a slice selection 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 encode gradient magnetic field pulse (Gp) and a frequency encode gradient magnetic field pulse (Gf) are applied in one direction, and position information in each direction is encoded in the echo signal.

シーケンサ4は、高周波磁場パルス(以下、「RFパルス」という)と傾斜磁場パルスをある所定のパルスシーケンスで繰り返し印加する制御手段で、演算処理部8の制御で動作し、被検体1の断層画像のデータ収集に必要な種々の命令を送信系5、傾斜磁場発生系3、および受信系6に送る。   The sequencer 4 is a control means that repeatedly applies a high-frequency magnetic field pulse (hereinafter referred to as an “RF pulse”) and a gradient magnetic field pulse in a predetermined pulse sequence, and operates under the control of the arithmetic processing unit 8 to obtain a tomographic image of the subject 1. Various commands necessary for data collection are sent 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とから成る。高周波発振器11から出力されたRFパルスをシーケンサ4からの指令によるタイミングで変調器12により振幅変調し、この振幅変調されたRFパルスを高周波増幅器13で増幅した後に被検体1に近接して配置された高周波コイル14aに供給することにより、RFパルスが被検体1に照射される。   The transmission system 5 irradiates the subject 1 with RF pulses 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 RF pulse output from the high-frequency oscillator 11 is amplitude-modulated by the modulator 12 at the timing according to the command from the sequencer 4, and the amplitude-modulated RF pulse is amplified by the high-frequency amplifier 13 and then placed close to the subject 1. By supplying to the high frequency coil 14a, the subject 1 is irradiated with the RF pulse.

受信系6は、被検体1の生体組織を構成する原子核スピンの核磁気共鳴により放出されるエコー信号(NMR信号)を検出するもので、受信側の高周波コイル(受信コイル)14bと信号増幅器15と直交位相検波器16と、A/D変換器17とから成る。送信側の高周波コイル14aから照射された電磁波によって誘起された被検体1の応答のNMR信号が被検体1に近接して配置された高周波コイル14bで検出され、信号増幅器15で増幅された後、シーケンサ4からの指令によるタイミングで直交位相検波器16により直交する二系統の信号に分割され、それぞれがA/D変換器17でディジタル量に変換されて、信号処理系7に送られる。   The receiving system 6 detects an echo signal (NMR signal) emitted by nuclear magnetic resonance of nuclear spins constituting the biological tissue of the subject 1, and receives a high-frequency coil (receiving coil) 14b on the receiving side and a signal amplifier 15 And a quadrature phase detector 16 and an A / D converter 17. After the NMR signal of the response of the subject 1 induced by the electromagnetic wave irradiated from the high frequency coil 14a on the transmission side is detected by the high frequency coil 14b arranged close to the subject 1 and amplified by the signal amplifier 15, The quadrature phase detector 16 divides the signal into two orthogonal signals at the timing according to the command from the sequencer 4, and each signal is converted into a digital quantity by the A / D converter 17 and sent to the signal processing system 7.

信号処理系7は、各種データ処理と処理結果の表示及び保存等を行うもので、光ディスク19、磁気ディスク18等の外部記憶装置と、CRT等からなるディスプレイ20とを有する。受信系6からのデータが演算処理部(CPU)8に入力されると、演算処理部8が信号処理、画像再構成等の処理を実行し、その結果である被検体1の断層画像をディスプレイ20に表示すると共に、外部記憶装置の磁気ディスク18等に記録する。
The signal processing system 7 performs various data processing and display and storage of processing results, and includes an external storage device such as an optical disk 19 and a magnetic disk 18, and a display 20 including a CRT or the like. When data from the reception system 6 is input to the arithmetic processing unit (CPU) 8, the arithmetic processing unit 8 executes processing such as signal processing and image reconstruction, and displays the tomographic image of the subject 1 as a result of the processing. 20 and recorded on the magnetic disk 18 of the external storage device.

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

なお、図1において、送信側の高周波コイル14aと傾斜磁場コイル9は、被検体1が挿入される静磁場発生系2の静磁場空間内に、垂直磁場方式であれば被検体1に対向して、水平磁場方式であれば被検体1を取り囲むようにして設置されている。一方、受信側の高周波コイル14bは、被検体1に対向して、或いは取り囲むように設置されている。なお、上記傾斜磁場発生系3、送信系5、受信系6、及びシーケンサ4を纏めて計測制御部ともいう。   In FIG. 1, the high-frequency coil 14a and the gradient magnetic field coil 9 on the transmission side face the subject 1 in the static magnetic field space of the static magnetic field generation system 2 into which the subject 1 is inserted, in the case of the vertical magnetic field method. If the horizontal magnetic field method is used, the subject 1 is installed so as to surround it. On the other hand, the high-frequency coil 14b on the reception side is disposed so as to face or surround the subject 1. The gradient magnetic field generation system 3, the transmission system 5, the reception system 6, and the sequencer 4 are collectively referred to as a measurement control unit.

現在MRI装置の撮像対象核種は、臨床で普及しているものとしては、被検体の主たる構成物質である水素原子核(プロトン)である。プロトン密度の空間分布や、励起状態の緩和時間の空間分布に関する情報を画像化することで、人体頭部、腹部、四肢等の形態または、機能を2次元もしくは3次元的に撮像する。   At present, the radionuclide to be imaged by the MRI apparatus is a hydrogen nucleus (proton) which is a main constituent material of the subject as being widely used clinically. By imaging information on the spatial distribution of proton density and the spatial distribution of relaxation time in the excited state, the form or function of the human head, abdomen, limbs, etc. is imaged two-dimensionally or three-dimensionally.

以下、本発明のMRI装置及び拡散強調画像計算方法の各実施例を説明する。
本発明は、b値=0を含む複数のb値でそれぞれ取得した拡散強調画像から計測拡散係数(ADCmrs)を算出する。
Embodiments of the MRI apparatus and the diffusion weighted image calculation method of the present invention will be described below.
In the present invention, a measured diffusion coefficient (ADC mrs ) is calculated from diffusion-weighted images acquired with a plurality of b values including b value = 0.

そして、組織毎に該計測拡散係数(ADCmrs)を基準拡散係数(ADCrfr)を用いて補正して補正拡散係数(ADCcrt)を得、該補正拡散係数(ADCcrt)を用いて任意のb値の拡散強調画像を算出する。 Then, by correcting using the measured diffusion coefficient for each tissue (ADC mrs) the reference diffusion coefficient (ADC RFR) to give the correction diffusion coefficient (ADC crt), any by using the correction diffusion coefficient (ADC crt) b-value diffusion weighted image is calculated.

或いは、計測拡散係数(ADCmrs)の修正するための情報に基づいて計測拡散係数(ADCmrs)を補正して補正拡散係数(ADCcrt)を得、該補正拡散係数(ADCcrt)を用いて任意のb値の拡散強調画像を算出する。 Alternatively, to obtain a measurement diffusion coefficient information measurement diffusion coefficient based on the order to correct the (ADC mrs) (ADC mrs) correction diffusion coefficient corrected (ADC crt), using the correction diffusion coefficient (ADC crt) A diffusion weighted image having an arbitrary b value is calculated.

本発明のMRI装置及び拡散強調画像計算方法に係る実施例1を説明する。実施例1では、b値=0を含む複数(2以上)のb値を用いてそれぞれ取得した拡散強調画像から、拡散係数(ADC)を算出し、該算出した拡散係数(ADC)をより正しい拡散係数(ADC)を用いて補正し、該補正した拡散係数(ADC)を用いて任意のb値の拡散強調画像を計算して表示する。以下、拡散強調画像から算出した拡散係数(ADC)を計測拡散係数(ADCmrs)、該補正後の計測拡散係数(ADCmrs)を補正拡散係数(ADCcrt)と記載する。 Example 1 according to the MRI apparatus and the diffusion weighted image calculation method of the present invention will be described. In Example 1, a diffusion coefficient (ADC) is calculated from diffusion weighted images acquired using a plurality of (two or more) b values including b value = 0, and the calculated diffusion coefficient (ADC) is more correct. Correction is performed using the diffusion coefficient (ADC), and a diffusion-weighted image having an arbitrary b value is calculated and displayed using the corrected diffusion coefficient (ADC). Hereinafter, the diffusion coefficient (ADC) calculated from the diffusion weighted image is referred to as a measured diffusion coefficient (ADC mrs ), and the corrected measured diffusion coefficient (ADC mrs ) is referred to as a corrected diffusion coefficient (ADC crt ).

図2に示す実施例1の処理フローを表すフローチャートに基づいて、実施例1を詳細に説明する。本処理フローは、予めプログラムとして磁気ディスク18に記憶されており、演算処理部8が磁気ディスク18からそのプログラムを読み込んで実行することにより実施される。   The first embodiment will be described in detail based on the flowchart showing the processing flow of the first embodiment shown in FIG. This processing flow is stored in advance on the magnetic disk 18 as a program, and is executed by the arithmetic processing unit 8 reading the program from the magnetic disk 18 and executing it.

(ステップS201)
操作者は、被検体を静磁場発生系2にセットアップする。そして、操作部25を介して、b値=0を含む複数のb値の各拡散強調画像と、T1強調画像及びT2強調画像の撮像を行うための拡散強調画像撮像シーケンスの撮像条件を設定する。
(Step S201)
The operator sets up the subject in the static magnetic field generation system 2. Then, the imaging condition of the diffusion weighted image capturing sequence for capturing each of the diffusion weighted images having a plurality of b values including b value = 0 and the T1 weighted image and the T2 weighted image is set via the operation unit 25. .

(ステップS202)
MRI装置は、ステップS201で設定された撮像条件で、b値=0を含む複数のb値の各拡散強調画像と、T1強調画像及びT2強調画像を被検体の同一断面で取得し、取得した各画像データを磁気ディスク18に記憶する。
(Step S202)
The MRI apparatus acquires and acquires each of the diffusion-weighted images having a plurality of b values including b value = 0, the T1-weighted images, and the T2-weighted images on the same cross section of the subject under the imaging conditions set in step S201. Each image data is stored in the magnetic disk 18.

(ステップS203)
操作者は、操作部25を介して、拡散強調画像の計算に用いる任意のb値を入力する。図3に、任意のb値の入力を受け付けるためのディスプレイ20に表示されるGUIの一例を示す。操作者は、図3に示すGUI上で入力エリア301にb値(bfactor)を直接入力するか、又は、スライダーバー302の操作によりb値(bfactor)の入力を行う。
(Step S203)
The operator inputs an arbitrary b value used for calculation of the diffusion weighted image via the operation unit 25. FIG. 3 shows an example of a GUI displayed on the display 20 for accepting input of an arbitrary b value. The operator directly inputs the b value (bfactor) into the input area 301 on the GUI shown in FIG. 3 or inputs the b value (bfactor) by operating the slider bar 302.

(ステップS204)
演算処理部8は、ステップS202で取得されたb値=0の拡散強調画像S0と、複数のb値の内の一つのb値の拡散強調画像Sとを用いて、計測拡散係数(ADCmrs)の算出を行う。計測拡散係数(ADCmrs)の算出は以下の式に従う。

Figure 0006579908
(Step S204)
The arithmetic processing unit 8 uses the diffusion-weighted image S 0 with b value = 0 acquired in step S202 and the diffusion-weighted image S with one b value among a plurality of b values, to measure a diffusion coefficient (ADC). mrs ). The measurement diffusion coefficient (ADC mrs ) is calculated according to the following equation.
Figure 0006579908

ここで、
ADCmrs:計測された見かけの拡散係数(計測拡散係数)
S:b値の拡散強調画像
S0:b値=0の拡散強調画像
b:任意のb値
他のb値の拡散強調画像(S)についても、上記計算式に基づいて計測拡散係数(ADCmrs)を算出する。図4に複数のb値で算出された計測拡散係数(ADCmrs)の一例401を示す。
here,
ADC mrs : Measured apparent diffusion coefficient (measured diffusion coefficient)
S: b value diffusion weighted image
S 0 : diffusion weighted image with b value = 0
b: Arbitrary b value The measurement diffusion coefficient (ADC mrs ) is also calculated based on the above formula for the diffusion weighted image (S) of other b values. FIG. 4 shows an example 401 of the measured diffusion coefficient (ADC mrs ) calculated with a plurality of b values.

複数のb値を用いてADCを求める場合は、上記式に各値の組(bi,Si)を代入して以下の式群を作り、それらを最小二乗法で解くことで、ADCを算出する。
1n(Sl)=1n(S0)−ADCXb1
1n(S2)=1n(S0)−ADCXb2
・・・
When calculating ADC using multiple b values, substituting each set of values (b i , S i ) into the above equation to create the following group of equations and solving them using the least squares method, calculate.
1n (S l ) = 1n (S 0 ) −ADCXb 1
1n (S 2 ) = 1n (S 0 ) −ADCXb 2
...

(ステップS205)
演算処理部8は、ステップS204で算出した計測拡散係数(ADCmrs)を、より正しい拡散係数(ADC)として事前に用意しておいた基準拡散係数(ADCrfr)を用いて補正し、補正拡散係数(ADCcrt)を得る。このように基準拡散係数(ADCrfr)に基づいて計測拡散係数(ADCmrs)を補正することにより、該計測拡散係数(ADCmrs)を得たMRI装置の静磁場強度に依存しない、より適切な補正拡散係数(ADCcrt)を得ることができる。
(Step S205)
The arithmetic processing unit 8 corrects the measured diffusion coefficient (ADC mrs ) calculated in step S204 using a reference diffusion coefficient (ADC rfr ) prepared in advance as a more correct diffusion coefficient (ADC), and corrects diffusion. Get the coefficient (ADC crt ). In this way, by correcting the measurement diffusion coefficient (ADC mrs ) based on the reference diffusion coefficient (ADC rfr ), the measurement diffusion coefficient (ADC mrs ) does not depend on the static magnetic field strength of the MRI apparatus from which the measurement diffusion coefficient (ADC mrs ) was obtained. A corrected diffusion coefficient (ADC crt ) can be obtained.

具体的には、演算処理部8は、ステップS202で取得されたT1強調画像、T2強調画像、及び拡散強調画像の画素毎に画素値の比を求める。そして、求めた比を磁気ディスク18に事前に保存しておいた組織毎の画素値の比と比較し、比の値が最も近い組織を画素毎に決定する。   Specifically, the arithmetic processing unit 8 obtains a pixel value ratio for each pixel of the T1-weighted image, the T2-weighted image, and the diffusion-weighted image acquired in step S202. Then, the obtained ratio is compared with the ratio of the pixel value for each tissue stored in advance on the magnetic disk 18, and the tissue having the closest ratio value is determined for each pixel.

そして、決定した組織毎にステップS204で算出した計測拡散係数(ADCmrs)の平均値を算出し、事前に用意しておいた組織毎の基準拡散係数(ADCrfr)との差を算出する。最後に、画素毎の計測拡散係数(ADCmrs)に対して、該画素の組織に対応する算出した拡散係数(ADC)の差を加算する。このようにして、演算処理部8は、組織毎に、基準拡散係数(ADCrfr)を用いて計測拡散係数(ADCmrs)を補正して補正拡散係数(ADCcrt)を得る。この画素毎の拡散係数(ADC)の差を、図10に示すように、拡散係数(ADC)補正値画像1001としてディスプレイ20に表示しても良い。 Then, the average value of the measured diffusion coefficient (ADC mrs ) calculated in step S204 is calculated for each determined tissue, and the difference from the reference diffusion coefficient (ADC rfr ) for each tissue prepared in advance is calculated. Finally, the difference of the calculated diffusion coefficient (ADC) corresponding to the organization of the pixel is added to the measured diffusion coefficient (ADC mrs ) for each pixel. In this way, the arithmetic processing unit 8 corrects the measured diffusion coefficient (ADC mrs ) using the reference diffusion coefficient (ADC rfr ) for each tissue to obtain a corrected diffusion coefficient (ADC crt ). The difference in diffusion coefficient (ADC) for each pixel may be displayed on the display 20 as a diffusion coefficient (ADC) correction value image 1001, as shown in FIG.

なお、上記のように磁気ディスク18に事前に保存しておく組織毎の画素値の比と基準拡散係数(ADCrfr)については後述する。 The ratio of pixel values for each tissue and the reference diffusion coefficient (ADC rfr ) stored in advance on the magnetic disk 18 as described above will be described later.

(ステップS206)
演算処理部8は、ステップS205で取得した補正拡散係数(ADCcrt)と、ステップS203で入力された任意のb値(bfactor)を用いて、任意のb値の拡散強調画像(Sbfactor)を以下の式を用いて計算する。

Figure 0006579908
(Step S206)
The arithmetic processing unit 8 uses the corrected diffusion coefficient (ADC crt ) acquired in step S205 and the arbitrary b value (bfactor) input in step S203 to generate a diffusion weighted image (S bfactor ) of an arbitrary b value. Calculate using the following formula:
Figure 0006579908

ここで、
ADCcrt:補正拡散係数
S0:b値=0の拡散強調画像
Sbfactor:任意のb値の拡散強調画像
bfactor:任意のb値
here,
ADC crt : Corrected diffusion coefficient
S 0 : diffusion weighted image with b value = 0
S bfactor : diffusion-weighted image with an arbitrary b value
bfactor: any b value

(ステップS207)
演算処理部8は、ステップ206で算出した任意のb値の拡散強調画像をディスプレイ20に表示する。
(Step S207)
The arithmetic processing unit 8 displays the diffusion weighted image of the arbitrary b value calculated in step 206 on the display 20.

以上迄が、実施例1の処理フローの説明である。   The above is the description of the processing flow of the first embodiment.

次に、上記ステップS205で必要となる、磁気ディスク18に事前に保存しておく、T1強調画像、T2強調画像、及び拡散強調画像の組織毎の画素値の比と基準拡散係数(ADCrfr)について説明する。実施例1では、高磁場(例えば、1.2テスラ、1.5テスラ、3テスラ、及び7テスラ等のいずれか)MRI装置で撮像した拡散強調画像から算出した計測拡散係数(ADCmrs)を基準拡散係数(ADCrfr)とする。図4に、事前に用意した基準拡散係数(ADCrfr)の一例402を示す。 Next, the ratio of pixel values for each tissue of the T1-weighted image, the T2-weighted image, and the diffusion-weighted image and the reference diffusion coefficient (ADC rfr ), which are required to be stored in the magnetic disk 18 in advance in step S205. Will be described. In Example 1, a measurement diffusion coefficient (ADC mrs ) calculated from a diffusion-weighted image captured by a high magnetic field (for example, any of 1.2 Tesla, 1.5 Tesla, 3 Tesla, 7 Tesla, etc.) MRI apparatus is used as a reference diffusion coefficient ( ADC rfr ). FIG. 4 shows an example 402 of a reference diffusion coefficient (ADC rfr ) prepared in advance.

具体的には次の通りに求める。高磁場MRI装置で撮像されたT1強調画像とT2強調画像を用いて、画素毎に組織を判別する。この画素毎の組織の判別は、操作者が視覚的に指定しても良いし、公知のアルゴリズムを用いて演算処理部8が行っても良い。そして、演算処理部8は、上記画素毎の組織を判別したT1強調画像とT2強調画像、及び、同じ高磁場MRI装置で撮像された拡散強調画像の画素値の平均値の比を組織毎に算出し、磁気ディスク18に保存しておく。また、演算処理部8は、上記高磁場MRI装置で撮像された拡散強調画像から算出した拡散係数(ADC)の組織毎の平均値を求めて、組織毎の基準拡散係数(ADCrfr)として磁気ディスク18に保存しておく。 Specifically, it is calculated as follows. A tissue is discriminated for each pixel using a T1-weighted image and a T2-weighted image captured by a high magnetic field MRI apparatus. The discrimination of the tissue for each pixel may be visually specified by the operator, or may be performed by the arithmetic processing unit 8 using a known algorithm. Then, the arithmetic processing unit 8 determines the ratio of the average value of the pixel values of the T1-weighted image and the T2-weighted image obtained by determining the tissue for each pixel and the diffusion-weighted image captured by the same high magnetic field MRI apparatus for each tissue. Calculated and stored in the magnetic disk 18. In addition, the arithmetic processing unit 8 obtains an average value for each tissue of the diffusion coefficient (ADC) calculated from the diffusion weighted image captured by the high magnetic field MRI apparatus, and uses the magnetic field as a reference diffusion coefficient (ADC rfr ) for each tissue. Save to disk 18.

以上説明したように、本発明の実施例1では、演算処理部8が、b値=0を含む複数(2以上)のb値の拡散強調画像から計測拡散係数(ADCmrs)を算出し、組織毎に該計測拡散係数(ADCmrs)を基準拡散係数(ADCrfr)を用いて補正して補正拡散係数(ADCcrt)を得、該補正拡散係数(ADCcrt)を用いて任意のb値の拡散強調画像を算出する。これにより、静磁場強度に依存することなく、より正確な任意のb値の拡散強調画像を算出することができる。特に、高出力の傾斜磁場電源を備えていない低磁場MRI装置で、高いb値の拡散強調画像の撮像が困難な場合でも、高磁場MRI装置と同等の拡散強調画像を得ることが可能になる。 As described above, in the first embodiment of the present invention, the arithmetic processing unit 8 calculates a measured diffusion coefficient (ADC mrs ) from a plurality of (two or more) b-value diffusion weighted images including b value = 0, the measured diffusion coefficient for each tissue (ADC mrs) corrected using the reference diffusion coefficient (ADC RFR) to give the correction diffusion coefficient (ADC crt) with any b value using the correction diffusion coefficient (ADC crt) Is calculated. As a result, a more accurate diffusion-weighted image having an arbitrary b value can be calculated without depending on the static magnetic field strength. In particular, it is possible to obtain a diffusion-weighted image equivalent to that of a high-field MRI device even when it is difficult to capture a diffusion-weighted image with a high b value with a low-field MRI device that does not have a high-output gradient magnetic field power supply .

本発明のMRI装置及び拡散強調画像計算方法に係る実施例2を説明する。実施例2では、2以上のb値を用いて取得した拡散強調画像から計測拡散係数(ADCmrs)を算出し、補正の基準とする基準拡散係数(ADCrfr)を選択し、選択した基準拡散係数(ADCrfr)を用いて計測拡散係数(ADCmrs)を補正し、該補正により得た補正拡散係数(ADCcrc)と任意のb値とを用いて、拡散強調画像を算出する。 A second embodiment according to the MRI apparatus and the diffusion weighted image calculation method of the present invention will be described. In Example 2, a measurement diffusion coefficient (ADC mrs ) is calculated from a diffusion weighted image acquired using a b value of 2 or more, a reference diffusion coefficient (ADC rfr ) as a correction reference is selected, and the selected reference diffusion The measurement diffusion coefficient (ADC mrs ) is corrected using the coefficient (ADC rfr ), and a diffusion weighted image is calculated using the corrected diffusion coefficient (ADC crc ) obtained by the correction and an arbitrary b value.

図5は、本発明の実施例2の処理フローを示すフローチャートである。本処理フローは、予めプログラムとして磁気ディスク18に記憶されており、演算処理部8が磁気ディスク18からそのプログラムを読み込んで実行することにより実施される。以下、図5の各ステップについて詳細に説明する。なお、実施例1と同じステップには同じステップ符号を付け、説明を省略する。   FIG. 5 is a flowchart showing the processing flow of the embodiment 2 of the present invention. This processing flow is stored in advance on the magnetic disk 18 as a program, and is executed by the arithmetic processing unit 8 reading the program from the magnetic disk 18 and executing it. Hereinafter, each step of FIG. 5 will be described in detail. The same steps as those in the first embodiment are denoted by the same step symbols, and the description thereof is omitted.

(ステップS201〜S204)
実施例1と同じである。
(Steps S201 to S204)
Same as Example 1.

(ステップS301)
演算処理部8は、事前に用意しておいた複数の基準拡散係数(ADCrfr)の中から目的とする基準拡散係数(ADCrfr)を選択して、該選択した基準拡散係数(ADCrfr)に基づいて、ステップS204で算出した計測拡散係数(ADCmrs)を補正し補正拡張係数(ADCcrc)を得る。
(Step S301)
Processing unit 8, a plurality of reference spreading factor which had been prepared beforehand by selecting a reference spreading factor of interest from the (ADC rfr) (ADC rfr) , the selected criteria diffusion coefficient (ADC RFR) Based on the above, the measured diffusion coefficient (ADC mrs ) calculated in step S204 is corrected to obtain a corrected expansion coefficient (ADC crc ).

図6に基準拡散係数(ADCrfr)の選択に用いるGUIの一例を示す。操作者は、cDWI画像上に表示された直交2軸の内の横軸601をトラックボール又はマウス23でドラッグすることで、任意のb値を入力でき、縦軸602をトラックボール又はマウス23でドラッグすることで、任意の基準選択係数(ADCrfr)を選択できる。 FIG. 6 shows an example of a GUI used for selecting the reference diffusion coefficient (ADC rfr ). The operator can input an arbitrary b value by dragging the horizontal axis 601 of the two orthogonal axes displayed on the cDWI image with the trackball or the mouse 23, and the vertical axis 602 with the trackball or the mouse 23. An arbitrary reference selection coefficient (ADC rfr ) can be selected by dragging .

基準選択係数(ADCrfr)としては、例えば、1.2テスラ、1.5テスラ、3テスラ、及び7テスラ等の静磁場強度のMRI装置で取得した拡散係数(ADC)を用いることができ、これらの基準拡散係数(ADCrfr)を事前に磁気ディスク18に保存しておく。図7に各拡散係数の一例を示す。点線701が1.5テスラMRI装置で取得した拡散係数(ADC)を示し、一点鎖線702が3テスラMRI装置で取得した拡散係数(ADC)を示す。 As the reference selection coefficient (ADC rfr ), for example, the diffusion coefficient (ADC) acquired with an MRI apparatus having a static magnetic field strength such as 1.2 Tesla, 1.5 Tesla, 3 Tesla, and 7 Tesla can be used. The coefficient (ADC rfr ) is stored in advance on the magnetic disk 18. FIG. 7 shows an example of each diffusion coefficient. A dotted line 701 indicates the diffusion coefficient (ADC) acquired by the 1.5 Tesla MRI apparatus, and a dashed line 702 indicates the diffusion coefficient (ADC) acquired by the 3 Tesla MRI apparatus.

そして、図6に示すようなGUIを介して、操作者により、複数の基準拡散係数(ADCrfr)の中から目的とする基準拡散係数(ADCrfr)の選択がされると、演算処理部8は、選択された基準拡散係数(ADCrfr)に基づいて計測拡散係数(ADCmrs)を補正し、補正拡散係数(ADCcrt)を得る。補正方法は実施例1と同じである。 Then, when the target reference diffusion coefficient (ADC rfr ) is selected from a plurality of reference diffusion coefficients (ADC rfr ) by the operator via the GUI as shown in FIG. 6, the arithmetic processing unit 8 Corrects the measured diffusion coefficient (ADC mrs ) based on the selected reference diffusion coefficient (ADC rfr ) to obtain a corrected diffusion coefficient (ADC crt ). The correction method is the same as in the first embodiment.

なお、実施例2においても、前述の実施例1のステップS205で説明したように、上記のように選択された基準拡散件数(ADCrfr)に基づいて定まる画素毎の拡散係数(ADC)の差を、図11に示すように、拡散係数(ADC)補正値画像1101としてディスプレイ20に表示しても良い。 In Example 2, as described in Step S205 of Example 1 above, the difference in diffusion coefficient (ADC) for each pixel determined based on the reference diffusion number (ADC rfr ) selected as described above. 11 may be displayed on the display 20 as a diffusion coefficient (ADC) correction value image 1101 as shown in FIG.

(ステップS302)
演算処理部8は、ステップS301で算出した補正拡散係数(ADCcrt)と、ステップS203で設定された任意のb値とを用いて、拡散強調画像を計算する。計算方法は実施例1と同じである。
(Step S302)
The arithmetic processing unit 8 calculates a diffusion weighted image using the corrected diffusion coefficient (ADC crt ) calculated in step S301 and the arbitrary b value set in step S203. The calculation method is the same as in Example 1.

(ステップS206)
実施例1と同じである。
(Step S206)
Same as Example 1.

以上までが、実施例2の処理フローの説明である。   The above is the description of the processing flow of the second embodiment.

以上説明したように本発明の実施例2では、演算処理部8が、b値=0を含む複数(2以上)のb値の拡散強調画像を用いて計測拡散係数(ADCmrs)を算出し、複数の基準拡散係数(ADCrfr)の内から選択された基準拡散係数(ADCrfr)を用いて該計測拡散係数(ADCmrs)を補正して補正拡散係数(ADCcrt)を得、該補正拡散係数(ADCcrt)を用いて任意のb値の拡散強調画像を算出する。これにより、前述の実施例1と同様に、静磁場強度に依存することなく、より正確な任意のb値の拡散強調画像を算出することができる。特に、補正に用いる基準拡散係数(ADCrfr)を選択することから、任意の静磁場強度の高磁場MRI装置と同等の拡散強調画像を得ることが可能になる。 As described above, in the second embodiment of the present invention, the arithmetic processing unit 8 calculates a measurement diffusion coefficient (ADC mrs ) using a plurality of (two or more) b-value diffusion weighted images including b value = 0. Correcting the measured diffusion coefficient (ADC mrs ) using a reference diffusion coefficient (ADC rfr ) selected from among a plurality of reference diffusion coefficients (ADC rfr ) to obtain a corrected diffusion coefficient (ADC crt ) A diffusion-weighted image having an arbitrary b value is calculated using the diffusion coefficient (ADC crt ). As a result, as in the first embodiment, a more accurate diffusion-weighted image having an arbitrary b value can be calculated without depending on the static magnetic field strength. In particular, since the reference diffusion coefficient (ADC rfr ) used for correction is selected, a diffusion-weighted image equivalent to a high magnetic field MRI apparatus having an arbitrary static magnetic field strength can be obtained.

本発明のMRI装置及び拡散強調画像計算方法に係る実施例3を説明する。実施例3では、2以上のb値の拡散強調画像から計測拡散係数(ADCmrs)を算出し、操作者の入力により該計測拡散係数(ADCmrs)を補正して補正拡散係数(ADCcrt)を得、該補正により得た補正拡散係数(ADCcrc)と任意のb値とを用いて拡散強調画像を算出する。 A third embodiment according to the MRI apparatus and the diffusion weighted image calculation method of the present invention will be described. In Example 3, 2 or more to calculate a measured diffusion coefficient from the diffusion-weighted images of the b value (ADC mrs), the correction diffusion coefficient by correcting the measured diffusion coefficient (ADC mrs) by the operator's input (ADC crt) A diffusion weighted image is calculated using the corrected diffusion coefficient (ADC crc ) obtained by the correction and an arbitrary b value.

図8は、本発明の実施例3の処理フローを示すフローチャートである。本処理フローは、予めプログラムとして磁気ディスク18に記憶されており、演算処理部8が磁気ディスク18からそのプログラムを読み込んで実行することにより実施される。以下、図8の各ステップについて詳細に説明する。   FIG. 8 is a flowchart showing a process flow of the third embodiment of the present invention. This processing flow is stored in advance on the magnetic disk 18 as a program, and is executed by the arithmetic processing unit 8 reading the program from the magnetic disk 18 and executing it. Hereinafter, each step of FIG. 8 will be described in detail.

(ステップS201〜S204)
実施例1と同じである。
(Steps S201 to S204)
Same as Example 1.

(ステップS401)
操作者は、S204で算出した計測拡散係数(ADCmrs)を補正するための情報を、操作部25を介して入力する。そして、演算処理部8は、該入力情報に基づいて計測拡散係数(ADCmrs)を補正して補正拡散係数(ADCcrt)を得る。そのための、GUIの一例を図9に示す。演算処理部8は、計測拡散係数(ADCmrs)のグラフ901をディスプレイ20に表示する。そして、操作者は、操作部25を介して、計測拡散係数(ADCmrs)のグラフ901の形状を修正し、演算処理部8は修正されたグラフ形状902に基づいて補正拡散係数(ADCcrt)を求める。
(Step S401)
The operator inputs information for correcting the measured diffusion coefficient (ADC mrs ) calculated in S204 via the operation unit 25. Then, the arithmetic processing unit 8 corrects the measured diffusion coefficient (ADC mrs ) based on the input information to obtain a corrected diffusion coefficient (ADC crt ). An example of a GUI for this purpose is shown in FIG. The arithmetic processing unit 8 displays a graph 901 of the measured diffusion coefficient (ADC mrs ) on the display 20. Then, the operator corrects the shape of the graph 901 of the measured diffusion coefficient (ADC mrs ) via the operation unit 25, and the arithmetic processing unit 8 corrects the corrected diffusion coefficient (ADC crt ) based on the corrected graph shape 902. Ask for.

(ステップS402)
演算処理部8は、ステップS401で算出した補正拡散係数(ADCcrt)と、ステップS203で設定された任意のb値を用いて、拡散強調画像を計算する。計算方法は実施例1と同じである。
(Step S402)
The arithmetic processing unit 8 calculates a diffusion weighted image using the corrected diffusion coefficient (ADC crt ) calculated in step S401 and the arbitrary b value set in step S203. The calculation method is the same as in Example 1.

(ステップS206)
実施例1と同じである。
(Step S206)
Same as Example 1.

以上までが、実施例3の処理フローの説明である。   The above is the description of the processing flow of the third embodiment.

以上説明したように本発明の実施例3では、計測拡散係数(ADCmrs)を修正するための情報を受け付けるGUIを表示し、演算処理部8が、GUIを介して入力された情報に基づいて計測拡散係数(ADCmrs)を補正して補正拡散係数(ADCcrt)を得、該補正拡散係数(ADCcrt)を用いて任意のb値の拡散強調画像を算出する。これにより、前述の実施例1,2と同様に、静磁場強度に依存することなく、拡散強調画像における拡散の程度を強調することができ、診断を容易にすることができる。 As described above, in the third embodiment of the present invention, the GUI for receiving information for correcting the measurement diffusion coefficient (ADC mrs ) is displayed, and the arithmetic processing unit 8 is based on the information input via the GUI. correcting the measured diffusion coefficient (ADC mrs) to give the correction diffusion coefficient (ADC crt) and to calculate the diffusion weighted image of any b value using the correction diffusion coefficient (ADC crt). As a result, as in the first and second embodiments, the degree of diffusion in the diffusion weighted image can be enhanced without depending on the static magnetic field strength, and diagnosis can be facilitated.

以上、本発明の実施例を述べたが、本発明はこれらに限定されるものではないことは言うまでも無い。   As mentioned above, although the Example of this invention was described, it cannot be overemphasized that this invention is not limited to these.

1 被検体、2 静磁場発生系、3 傾斜磁場発生系、4 シーケンサ、5 送信系、6 受信系、7 信号処理系、8 演算処理部(CPU)、9 傾斜磁場コイル、10 傾斜磁場電源、11 高周波発振器、12 変調器、13 高周波増幅器、14a 高周波コイル(送信コイル)、14b 高周波コイル(受信コイル)、15 信号増幅器、16 直交位相検波器、17 A/D変換器、18 磁気ディスク、19 光ディスク、20 ディスプレイ、23 トラックボール又はマウス、24 キーボード   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 arithmetic processing unit (CPU), 9 gradient magnetic field coil, 10 gradient magnetic field power supply, 11 High-frequency oscillator, 12 modulator, 13 high-frequency amplifier, 14a high-frequency coil (transmitting coil), 14b high-frequency coil (receiving coil), 15 signal amplifier, 16 quadrature detector, 17 A / D converter, 18 magnetic disk, 19 Optical disc, 20 display, 23 trackball or mouse, 24 keyboard

Claims (10)

被検体を収容する空間に均一な静磁場を発生させる静磁場発生部と、
拡散強調シーケンスに従い、前記被検体にb値=0を含む任意のb値のMPGを印加して、該被検体から発生するNMR信号を検出する計測制御部と、
前記検出されたNMR信号から前記b値毎の拡散強調画像を取得する演算処理部と、
を備えた磁気共鳴イメージング装置であって、
前記演算処理部は、前記拡散強調画像から計測拡散係数(ADCmrs)を算出し、組織毎に該計測拡散係数(ADCmrs)を基準拡散係数(ADCrfr)を用いて補正して補正拡散係数(ADCcrt)を得、該補正拡散係数(ADCcrt)を用いて任意のb値の拡散強調画像を算出することを特徴とする磁気共鳴イメージング装置。
A static magnetic field generator for generating a uniform static magnetic field in a space for accommodating a subject;
In accordance with a diffusion weighted sequence, a measurement control unit that detects an NMR signal generated from the subject by applying an MPG having an arbitrary b value including b value = 0 to the subject;
An arithmetic processing unit for acquiring a diffusion weighted image for each b value from the detected NMR signal;
A magnetic resonance imaging apparatus comprising:
The arithmetic processing unit, the calculated measurement diffusion coefficient (ADC mrs) from the diffusion weighted image, correction diffusion coefficient is corrected using the measured diffusion coefficient for each tissue (ADC mrs) the reference diffusion coefficient (ADC RFR) A magnetic resonance imaging apparatus characterized in that (ADC crt ) is obtained, and a diffusion weighted image having an arbitrary b value is calculated using the corrected diffusion coefficient (ADC crt ).
請求項1記載の磁気共鳴イメージング装置において、
前記基準拡散係数(ADCrfr)は、1.2テスラ、1.5テスラ、3テスラ、及び7テスラのいずれかの静磁場強度で得られた拡散強調画像から得られたものであることを特徴とする磁気共鳴イメージング装置。
The magnetic resonance imaging apparatus according to claim 1.
The reference diffusion coefficient (ADC rfr ) is obtained from a diffusion weighted image obtained with a static magnetic field strength of 1.2 Tesla, 1.5 Tesla, 3 Tesla, or 7 Tesla. Magnetic resonance imaging apparatus.
請求項1又は2記載の磁気共鳴イメージング装置において、
前記演算処理部は、組織毎の前記計測拡散係数(ADCmrs)の平均値と組織毎の前記基準拡散係数(ADCrfr)との差を求め、画素毎に該画素の組織に対応する前記差を該計測拡散係数(ADCmrs)に加算して、前記補正拡散係数(ADCcrt)を求めることを特徴とする磁気共鳴イメージング装置。
The magnetic resonance imaging apparatus according to claim 1 or 2,
The arithmetic processing unit obtains a difference between an average value of the measured diffusion coefficient (ADC mrs ) for each tissue and the reference diffusion coefficient (ADC rfr ) for each tissue, and the difference corresponding to the tissue of the pixel for each pixel. Is added to the measured diffusion coefficient (ADC mrs ) to obtain the corrected diffusion coefficient (ADC crt ).
請求項3記載の磁気共鳴イメージング装置において、
前記組織毎の計測拡散係数(ADCmrs)の平均値と、前記組織毎の基準拡散係数(ADCrfr)との差を画像として表示することを特徴とする磁気共鳴イメージング装置。
The magnetic resonance imaging apparatus according to claim 3.
A magnetic resonance imaging apparatus, wherein a difference between an average value of measured diffusion coefficients (ADC mrs ) for each tissue and a reference diffusion coefficient (ADC rfr ) for each tissue is displayed as an image.
請求項1乃至4のいずれか一項に記載の磁気共鳴イメージング装置において、
前記演算処理部は、異なる静磁場強度で得られた複数の前記基準拡散係数(ADCrfr)の内から選択された基準拡散係数(ADCrfr)を用いて前記補正拡散係数(ADCcrt)を求めることを特徴とする磁気共鳴イメージング装置。
The magnetic resonance imaging apparatus according to any one of claims 1 to 4,
The arithmetic processing unit obtains the corrected diffusion coefficient (ADC crt ) using a reference diffusion coefficient (ADC rfr ) selected from among the plurality of reference diffusion coefficients (ADC rfr ) obtained with different static magnetic field strengths. A magnetic resonance imaging apparatus.
請求項5に記載の磁気共鳴イメージング装置において、
前記複数の基準拡散係数(ADCrfr)の内から一つの基準拡散係数(ADCrfr)の選択を受け付けるためのGUIを表示することを特徴とする磁気共鳴イメージング装置。
The magnetic resonance imaging apparatus according to claim 5.
A magnetic resonance imaging apparatus for displaying a GUI for receiving selection of one reference diffusion coefficient (ADC rfr ) from among the plurality of reference diffusion coefficients (ADC rfr ).
請求項1乃至6のいずれか一項に記載の磁気共鳴イメージング装置において、
組織毎の前記基準拡散係数(ADCrfr)を記憶する記憶手段を備えていることを特徴とする磁気共鳴イメージング装置。
The magnetic resonance imaging apparatus according to any one of claims 1 to 6,
A magnetic resonance imaging apparatus comprising storage means for storing the reference diffusion coefficient (ADC rfr ) for each tissue.
請求項1乃至のいずれか一項に記載の磁気共鳴イメージング装置において、
前記任意のb値の入力を受け付けるためのGUIを表示することを特徴とする磁気共鳴イメージング装置。
The magnetic resonance imaging apparatus according to any one of claims 1 to 7 ,
A magnetic resonance imaging apparatus for displaying a GUI for receiving an input of the arbitrary b value.
被検体を収容する空間に均一な静磁場を発生させる静磁場発生部と、
拡散強調シーケンスに従い、前記被検体にb値=0を含む任意のb値のMPGを印加して、該被検体から発生するNMR信号を検出する計測制御部と、
前記検出されたNMR信号から前記b値毎の拡散強調画像を取得する演算処理部と、
を備えた磁気共鳴イメージング装置における拡散強調画像計算方法であって、
前記拡散強調画像から計測拡散係数(ADCmrs)を算出するステップと、
組織毎に前記計測拡散係数(ADCmrs)を基準拡散係数(ADCrfr)を用いて補正して補正拡散係数(ADCcrt)を得るステップと、
前記補正拡散係数(ADCcrt)を用いて任意のb値の拡散強調画像を算出するステップと
を有することを特徴とする拡散強調画像計算方法。
A static magnetic field generator for generating a uniform static magnetic field in a space for accommodating a subject;
In accordance with a diffusion weighted sequence, a measurement control unit that detects an NMR signal generated from the subject by applying an MPG having an arbitrary b value including b value = 0 to the subject;
An arithmetic processing unit for acquiring a diffusion weighted image for each b value from the detected NMR signal;
A diffusion-weighted image calculation method in a magnetic resonance imaging apparatus comprising:
Calculating a measured diffusion coefficient (ADC mrs ) from the diffusion weighted image;
Correcting the measured diffusion coefficient (ADC mrs ) for each tissue using a reference diffusion coefficient (ADC rfr ) to obtain a corrected diffusion coefficient (ADC crt );
And calculating a diffusion-weighted image having an arbitrary b value using the corrected diffusion coefficient (ADC crt ).
請求項9記載の拡散強調画像計算方法において、前記補正拡散係数(ADCcrt)を得るステップは、
前記拡散強調画像の組織毎に前記計測拡散係数(ADCmrs)の平均値を求めるステップと、
組織毎の前記計測拡散係数(ADCmrs)の平均値と組織毎の前記基準拡散係数(ADCrfr)との差を求めるステップと、
組織毎に前記差を前記計測拡散係数(ADCmrs)に加算して、前記補正拡散係数(ADCcrt)を求めるステップと
を有して成ることを特徴とする拡散強調画像計算方法。
10. The diffusion weighted image calculation method according to claim 9 , wherein the step of obtaining the corrected diffusion coefficient (ADC crt ) includes:
Obtaining an average value of the measured diffusion coefficient (ADC mrs ) for each tissue of the diffusion weighted image;
Determining a difference between an average value of the measured diffusion coefficient (ADC mrs ) for each tissue and the reference diffusion coefficient (ADC rfr ) for each tissue;
Adding the difference to the measured diffusion coefficient (ADC mrs ) for each tissue to obtain the corrected diffusion coefficient (ADC crt ).
JP2015206108A 2015-10-20 2015-10-20 Magnetic resonance imaging apparatus and diffusion weighted image calculation method Active JP6579908B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2015206108A JP6579908B2 (en) 2015-10-20 2015-10-20 Magnetic resonance imaging apparatus and diffusion weighted image calculation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2015206108A JP6579908B2 (en) 2015-10-20 2015-10-20 Magnetic resonance imaging apparatus and diffusion weighted image calculation method

Publications (2)

Publication Number Publication Date
JP2017077316A JP2017077316A (en) 2017-04-27
JP6579908B2 true JP6579908B2 (en) 2019-09-25

Family

ID=58665753

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2015206108A Active JP6579908B2 (en) 2015-10-20 2015-10-20 Magnetic resonance imaging apparatus and diffusion weighted image calculation method

Country Status (1)

Country Link
JP (1) JP6579908B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117233676B (en) * 2023-11-15 2024-03-26 之江实验室 Echo time-dependent magnetic resonance diffusion imaging signal generation method and device

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61276073A (en) * 1985-05-31 1986-12-06 Shimadzu Corp Formation of picture for nmr-ct device
JP5619448B2 (en) * 2009-08-20 2014-11-05 株式会社東芝 Magnetic resonance imaging system
WO2012096105A1 (en) * 2011-01-11 2012-07-19 株式会社 東芝 Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP2014195532A (en) * 2013-03-29 2014-10-16 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Estimation device, magnetic resonance device, program, and estimation method
DE102013224264B4 (en) * 2013-11-27 2018-12-27 Siemens Healthcare Gmbh Method for processing magnetic resonance diffusion image data

Also Published As

Publication number Publication date
JP2017077316A (en) 2017-04-27

Similar Documents

Publication Publication Date Title
KR101301490B1 (en) Magnetic resonance imaging device and diffusion weighted image acquiring method thereof
US9297876B2 (en) Magnetic resonance imaging apparatus and eddy current compensation method
JP6162142B2 (en) Magnetic resonance imaging apparatus and SAR prediction method
JP6464088B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
EP2762910A1 (en) Susceptibility-weighted magnetic resonance imaging
US10203387B2 (en) MR imaging with enhanced susceptibility contrast
EP2526440B1 (en) Susceptibility gradient mapping
JP5666470B2 (en) Nuclear magnetic resonance imaging apparatus and SAR estimation method
JP2016533215A (en) MR imaging with temperature mapping
JP6618988B2 (en) Magnetic resonance imaging apparatus and RF shimming parameter setting method
JP5618683B2 (en) Magnetic resonance imaging apparatus and brightness non-uniformity correction method
JP6579908B2 (en) Magnetic resonance imaging apparatus and diffusion weighted image calculation method
JP2013121437A (en) MAGNETIC RESONANCE IMAGING APPARATUS AND T1ρ IMAGING METHOD
JP6718764B2 (en) Magnetic resonance imaging apparatus and control method thereof
JP5718148B2 (en) Magnetic resonance imaging apparatus and dual slice measurement method
JP5421600B2 (en) Nuclear magnetic resonance imaging apparatus and method of operating nuclear magnetic resonance imaging apparatus
JP2006061235A (en) Magnetic resonance imaging device
JP6157976B2 (en) Magnetic resonance imaging apparatus and method
JP2012095891A (en) Magnetic resonance imaging apparatus
JP2013183967A (en) Magnetic resonance imaging apparatus and phase correction method
US11914017B2 (en) Magnetic resonance imaging apparatus and image processing method
JP6169909B2 (en) Magnetic resonance imaging apparatus and real number component image acquisition method
JP5837354B2 (en) Magnetic resonance imaging apparatus and magnetic resonance spectroscopy imaging method
JP2016131847A (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP5638324B2 (en) Magnetic resonance imaging apparatus and image correction method

Legal Events

Date Code Title Description
A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20160427

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20171031

RD04 Notification of resignation of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7424

Effective date: 20171107

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180905

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20190417

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190521

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190718

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190827

R150 Certificate of patent or registration of utility model

Ref document number: 6579908

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350