JP5283092B2 - Magnetic resonance imaging apparatus and image processing apparatus - Google Patents
Magnetic resonance imaging apparatus and image processing apparatus Download PDFInfo
- Publication number
- JP5283092B2 JP5283092B2 JP2011133906A JP2011133906A JP5283092B2 JP 5283092 B2 JP5283092 B2 JP 5283092B2 JP 2011133906 A JP2011133906 A JP 2011133906A JP 2011133906 A JP2011133906 A JP 2011133906A JP 5283092 B2 JP5283092 B2 JP 5283092B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- flow velocity
- unit
- imaging
- magnetic resonance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Description
この発明は、被検体内を撮像した画像を処理する磁気共鳴イメージング装置および画像処理装置に関する。 The present invention relates to a magnetic resonance imaging apparatus and an image processing apparatus that process an image captured in a subject.
従来、磁気共鳴イメージング装置を用いてCSFや血液等の体液の動態を観測する方法として、フローエンコード傾斜磁場を利用した「フェイズシフト(phase shift)法」または「フェイズコントラスト法」と呼ばれる方法がある(例えば、非特許文献1参照)。 Conventionally, as a method for observing the dynamics of body fluids such as CSF and blood using a magnetic resonance imaging apparatus, there is a method called “phase shift method” or “phase contrast method” using a flow encode gradient magnetic field. (For example, refer nonpatent literature 1).
図14は、従来のフェイズシフト法のパルスシーケンスを示す図である。図14に示すように、通常、フェイズシフト法では、グラジエントエコー法が用いられ、励起パルス(図14に示すP1)とエコー信号(図14に示すP2)との間でフローエンコード傾斜磁場(Gfe)(図14に示すP3)を加えることによって、流速に比例した位相変化が画像に与えられる。 FIG. 14 is a diagram showing a pulse sequence of a conventional phase shift method. As shown in FIG. 14, the gradient shift method is normally used in the phase shift method, and a flow encode gradient magnetic field (Gfe) is generated between an excitation pulse (P1 shown in FIG. 14) and an echo signal (P2 shown in FIG. 14). ) (P3 shown in FIG. 14), a phase change proportional to the flow velocity is given to the image.
このフェイズシフト法は、基本的には、対象の流速変化が心電波形(図14に示すECG)に強く相関していることを前提としている。すなわち、図14に示すように、心電同期波形のR波から一定の待ち時間(図14に示すTdelay)が経過した後に高周波励起が行われて、エコー信号が収集される。ここで、通常、1つの画像を再構成するためには、128個から256個程度のエコー信号が必要である。そのため、一般的には、フェイズシフト法では、1回のR波ごとに1つのエコー信号を収集する手順が、位相エンコード傾斜磁場パルスを段階的に変化させながら128回から256回程度繰り返される。 This phase shift method is based on the premise that the flow velocity change of the object is strongly correlated with the electrocardiogram waveform (ECG shown in FIG. 14). That is, as shown in FIG. 14, high-frequency excitation is performed after a certain waiting time (Tdelay shown in FIG. 14) has elapsed from the R wave of the electrocardiogram synchronization waveform, and echo signals are collected. Here, normally, in order to reconstruct one image, about 128 to 256 echo signals are required. Therefore, in general, in the phase shift method, the procedure of collecting one echo signal for each R wave is repeated about 128 to 256 times while gradually changing the phase encoding gradient magnetic field pulse.
こうして収集されたエコー信号に離散フーリエ変換などの再構成処理を施すことによって得られた画像S(ra)は、流速に比例した位相変化を有する。この画像S(ra)に対して、装置の不完全性や静磁場不均一性による補正処理を施した後に、例えばV(ra)=k(venc)・arg{S(ra)}で表される処理(ここで、k(venc)は、フローエンコード傾斜磁場のパルスの形状によって変化する比例係数)を施すことによって、流速の分布を表す流速画像V(ra)が得られる。 An image S (r a ) obtained by performing reconstruction processing such as discrete Fourier transform on the collected echo signals has a phase change proportional to the flow velocity. For example, V (r a ) = k (venc) · arg {S (r a )} after the image S (r a ) is subjected to correction processing due to incompleteness of the apparatus or static magnetic field inhomogeneity. (Where k (venc) is a proportional coefficient that varies depending on the pulse shape of the flow encode gradient magnetic field), thereby obtaining a flow velocity image V (r a ) representing the flow velocity distribution.
しかしながら、従来技術では、体液の流動による信号強度の変化を容易に把握することができない場合があった。 However, in the prior art, there is a case where a change in signal intensity due to the flow of body fluid cannot be easily grasped.
この発明は、上述した従来技術による課題を解決するためになされたものであり、体液の流動による信号強度の変化を容易に把握することができる磁気共鳴イメージング装置および画像処理装置を提供することを目的とする。 The present invention has been made to solve the above-described problems caused by the prior art, and provides a magnetic resonance imaging apparatus and an image processing apparatus capable of easily grasping a change in signal intensity due to a flow of body fluid. Objective.
上述した課題を解決し、目的を達成するため、本発明は、トリガ信号に同期して被検体の撮像領域に励起パルスを印加して該撮像領域に含まれるスピンに関するエコー信号を収集し、収集したエコー信号に基づいて磁気共鳴画像を生成する磁気共鳴イメージング装置であって、前記被検体内を流れる脳脊髄液を得るための撮像を複数回繰り返すことによって脳脊髄液に関する複数の再構成画像を作成する画像作成手段と、前記画像作成手段によって作成された複数の再構成画像を用いて、各再構成画像における同一位置ごとに前記脳脊髄液の流れに関する情報を算出し、当該算出された情報に基づく計算画像を作成する計算画像作成手段と、を備えたことを特徴とする。 In order to solve the above-described problems and achieve the object, the present invention collects echo signals related to spins included in the imaging region by applying an excitation pulse to the imaging region of the subject in synchronization with the trigger signal. A magnetic resonance imaging apparatus for generating a magnetic resonance image based on an echo signal obtained, wherein a plurality of reconstructed images relating to cerebrospinal fluid are obtained by repeating imaging for obtaining cerebrospinal fluid flowing in the subject a plurality of times. Using the image creation means to create and a plurality of reconstructed images created by the image creation means, information on the flow of the cerebrospinal fluid is calculated for each identical position in each reconstructed image, and the calculated information And a calculation image creating means for creating a calculation image based on the above.
また、本発明は、被検体内を流れる脳脊髄液を得るための撮像を複数回繰り返すことによって作成された脳脊髄液に関する複数の再構成画像を取得する画像取得手段と、前記画像取得手段によって取得された複数の再構成画像を用いて、各再構成画像における同一位置ごとに前記脳脊髄液の流れに関する情報を算出し、当該算出された情報に基づく計算画像を作成する計算画像作成手段と、を備えたことを特徴とする。 The present invention also provides an image acquisition means for acquiring a plurality of reconstructed images relating to cerebrospinal fluid created by repeating imaging for obtaining cerebrospinal fluid flowing in the subject a plurality of times, and the image acquisition means Calculation image creating means for calculating information on the flow of the cerebrospinal fluid for each same position in each reconstructed image using the plurality of acquired reconstructed images and creating a calculated image based on the calculated information; , Provided.
本発明によれば、体液の流動による信号強度の変化を容易に把握することができるという効果を奏する。 According to the present invention, there is an effect that it is possible to easily grasp a change in signal intensity due to the flow of body fluid.
以下に添付図面を参照して、この発明に係る磁気共鳴イメージング装置および画像処理装置の好適な実施例を詳細に説明する。なお、以下に示す実施例では、本発明を磁気共鳴イメージング装置(以下、「MRI(Magnetic Resonance Imaging)装置」と呼ぶ)に適用した場合について説明する。 Exemplary embodiments of a magnetic resonance imaging apparatus and an image processing apparatus according to the present invention will be described below in detail with reference to the accompanying drawings. In the following embodiments, the case where the present invention is applied to a magnetic resonance imaging apparatus (hereinafter referred to as “MRI (Magnetic Resonance Imaging) apparatus”) will be described.
まず、本実施例1に係るMRI装置の全体構成について説明する。図1は、本実施例1に係るMRI装置の全体構成を示す図である。図1に示すように、本実施例1に係るMRI装置100は、静磁場磁石1、傾斜磁場コイル2、傾斜磁場電源3、寝台4、寝台制御部5、送信RFコイル6、送信部7、受信RFコイル8、受信部9、シーケンス制御部10、架台部15および計算機システム20を備える。
First, the overall configuration of the MRI apparatus according to the first embodiment will be described. FIG. 1 is a diagram illustrating an overall configuration of the MRI apparatus according to the first embodiment. As shown in FIG. 1, the
静磁場磁石1は、中空の円筒形状に形成された磁石であり、内部の空間に一様な静磁場を発生する。この静磁場磁石1としては、例えば永久磁石、超伝導磁石等が使用される。
The static
傾斜磁場コイル2は、中空の円筒形状に形成されたコイルであり、静磁場磁石1の内側に配置される。この傾斜磁場コイル2は、互いに直交するX,Y,Zの各軸に対応する3つのコイルが組み合わされて形成されており、これら3つのコイルは、後述する傾斜磁場電源3から個別に電流供給を受けて、X,Y,Zの各軸に沿って磁場強度が変化する傾斜磁場を発生させる。なお、Z軸方向は、静磁場と同方向とする。傾斜磁場電源3は、傾斜磁場コイル2に電流を供給する装置である。
The gradient
ここで、傾斜磁場コイル2によって発生するX,Y,Z各軸の傾斜磁場は、例えば、スライス選択用傾斜磁場Gs、位相エンコード用傾斜磁場Geおよびリードアウト用傾斜磁場Grにそれぞれ対応している。スライス選択用傾斜磁場Gsは、任意に撮像断面を決めるために利用される。位相エンコード用傾斜磁場Geは、空間的位置に応じてエコー信号(磁気共鳴信号)の位相を変化させるために利用される。リードアウト用傾斜磁場Grは、空間的位置に応じてエコー信号の周波数を変化させるために利用される。
Here, the gradient magnetic fields of the X, Y, and Z axes generated by the gradient
寝台4は、被検体Pが載置される天板4aを備えた装置であり、後述する寝台制御部5による制御のもと、天板4aを、被検体Pが載置された状態で傾斜磁場コイル2の空洞(撮像口)内へ挿入する。通常、この寝台4は、長手方向が静磁場磁石1の中心軸と平行になるように設置される。寝台制御部5は、寝台4を制御する装置であり、計算機システム20による制御のもと、寝台4を駆動して天板4aを長手方向および上下方向へ移動する。
The
送信RFコイル6は、傾斜磁場コイル2の内側に配置されたコイルであり、送信部7から高周波パルスの供給を受けて、高周波磁場を発生する。送信部7は、ラーモア周波数に対応する高周波パルスを送信RFコイル6に送信する装置である。
The
受信RFコイル8は、傾斜磁場コイル2の内側に配置されたコイルであり、上記の高周波磁場の影響によって被検体Pから発せられるエコー信号を受信する。この受信RFコイル8は、エコー信号を受信すると、そのエコー信号を受信部9へ出力する。
The reception RF coil 8 is a coil disposed inside the gradient
受信部9は、受信RFコイル8から出力されるエコー信号に基づいてk空間データを生成する装置である。具体的には、この受信部9は、受信RFコイル8から出力されるエコー信号をデジタル変換することによってk空間データを生成する。このk空間データには、前述したスライス選択用傾斜磁場Gs、位相エンコード用傾斜磁場Geおよびリードアウト用傾斜磁場Grによって、PE方向、RO方向、SE方向の空間周波数の情報が対応付けられている。そして、k空間データを生成すると、受信部9は、そのk空間データをシーケンス制御部10へ送信する。
The receiving
シーケンス制御部10は、計算機システム20から送信されるシーケンス情報に基づいて、傾斜磁場電源3、送信部7および受信部9を駆動することによって、被検体Pのスキャンを行う装置である。ここで、シーケンス情報とは、傾斜磁場電源3が傾斜磁場コイル2に供給する電源の強さや電源を供給するタイミング、送信部7が送信RFコイル6に送信するRF信号の強さやRF信号を送信するタイミング、受信部9がエコー信号を検出するタイミングなど、スキャンを行うための手順を定義した情報である。
The
なお、シーケンス制御部10は、傾斜磁場電源3、送信部7および受信部9を駆動して被検体Pをスキャンした結果、受信部9からk空間データが送信されると、そのk空間データを計算機システム20へ転送する。
The
架台部15は、前述した静磁場磁石1、傾斜磁場コイル2、送信RFコイル6および受信RFコイル8を有し、静磁場中に置かれた被検体に高周波磁場を照射し、当該被検体から発せられるNMR信号を検出する。
The
計算機システム20は、MRI装置100の全体制御や、データ収集、画像再構成などを行う装置であり、インタフェース部21、画像再構成部22、記憶部23、入力部24、表示部25および制御部26を有している。
The
インタフェース部21は、シーケンス制御部10との間で授受される各種信号の入出力を制御する。例えば、このインタフェース部21は、シーケンス制御部10に対してシーケンス情報を送信し、シーケンス制御部10からk空間データを受信する。k空間データを受信すると、インタフェース部21は、各k空間データを被検体Pごとに記憶部23に格納する。
The
画像再構成部22は、記憶部23によって記憶されたk空間データに対して、後処理、すなわちフーリエ変換処理等の再構成処理を施すことによって、被検体P内における所望核スピンのスペクトラムデータあるいは画像データを生成する。なお、かかる画像再構成部22の構成については、後に詳細に説明する。
The
記憶部23は、インタフェース部21によって受信されたk空間データと、画像再構成部22によって生成された画像データなどを、被検体Pごとに記憶する。なお、かかる記憶部23の構成については、後に詳細に説明する。
The
入力部24は、操作者からの各種指示や情報入力を受け付ける装置である。この入力部24としては、マウスやトラックボールなどのポインティングデバイス、モード切替スイッチ等の選択デバイス、あるいはキーボード等の入力デバイスを適宜に利用可能である。
The
表示部25は、制御部26による制御のもと、スペクトラムデータあるいは画像データ等の各種の情報を表示する装置である。この表示部25としては、液晶表示器などの表示デバイスを利用可能である。
The
制御部26は、図示していないCPUやメモリ等を有し、MRI装置100の全体制御を行う。具体的には、この制御部26は、入力部24を介して操作者から入力される撮像条件に基づいてシーケンス情報を生成し、生成したシーケンス情報をシーケンス制御部10に送信することによってスキャンを制御したり、スキャンの結果としてシーケンス制御部10から送られるk空間データに基づいて行われる画像の再構成を制御したりする。なお、かかる制御部26の構成については、後に詳細に説明する。
The
次に、図1に示した画像再構成部22、記憶部23および制御部26の構成について説明する。図2は、図1に示した画像再構成部22、記憶部23および制御部26の構成を示す機能ブロック図である。
Next, the configuration of the
図2に示すように、記憶部23は、特に、k空間データ記憶部23aおよび画像データ記憶部23bを有する。k空間データ記憶部23aは、インタフェース部21によって受信されたk空間データを記憶する。画像データ記憶部23bは、画像再構成部22によって生成された画像データを記憶する。
As shown in FIG. 2, the
画像再構成部22は、特に、フーリエ変換処理部22a、絶対値画像作成部22b、流速画像作成部22c、流速分散画像作成部22dおよび重畳画像処理部22eを有する。
The
フーリエ変換処理部22aは、k空間データ記憶部23aに記憶されたk空間データに離散2次元フーリエ変換などの再構成処理を施すことによって画像を再構成する。
The Fourier
ここで、本実施例1では、被検体内を流れる体液の流速成分が得られる撮像法による撮像を複数回繰り返すことによって得られた複数の画像が用いられる。具体的には、本実施例1では、1回の励起パルスで1画像を再構成し、かつ、再構成した画像に対して流速に比例した位相変化を与えることが可能なEPI法による撮像を被検体の心電波形に同期して複数回繰り返すことによって得られた複数の画像が用いられる。 Here, in the first embodiment, a plurality of images obtained by repeating imaging by an imaging method that obtains the flow velocity component of the body fluid flowing in the subject are used a plurality of times. More specifically, in the first embodiment, one image is reconstructed by one excitation pulse, and imaging by the EPI method that can give a phase change proportional to the flow velocity to the reconstructed image. A plurality of images obtained by repeating a plurality of times in synchronization with the electrocardiographic waveform of the subject are used.
図3は、本実施例1で用いられるEPI法のパルスシーケンスを示す図である。図3に示すように、本実施例1で用いられるEPI法では、心電波形(図3に示すECG)のR波から一定の待ち時間(図3に示すTdelay)が経過した後に、90度励起パルス(図3に示すP11)を印加することによって撮像断面内が励起される。この90度励起パルスとエコー信号収集部分の間(図3に示すP13)でフローエンコード傾斜磁場(図3に示すP12)を加えることによって、得られる画像の位相が撮像対象の流速に比例して変化する。 FIG. 3 is a diagram showing a pulse sequence of the EPI method used in the first embodiment. As shown in FIG. 3, in the EPI method used in the first embodiment, 90 degrees after a certain waiting time (Tdelay shown in FIG. 3) has elapsed from the R wave of the electrocardiogram waveform (ECG shown in FIG. 3). By applying an excitation pulse (P11 shown in FIG. 3), the imaging section is excited. By applying a flow encode gradient magnetic field (P12 shown in FIG. 3) between the 90-degree excitation pulse and the echo signal collection part (P13 shown in FIG. 3), the phase of the obtained image is proportional to the flow velocity of the imaging target. Change.
なお、ここでは、SE(Spin Echo)タイプのEPIを一例として示したが、例えば、FE(Field Echo)タイプのEPIや、その他、位相画像が得られる撮像法であればいずれの撮像法が用いられてもよい。望ましくは、なるべく短い時間で画像が得られる撮像方法を用いるのがよい。例えば、「Proc. of Annual Meeting, Society of Magnetic Resonance in Medicine, No.2138, 1998.」に記載されている撮像法もその一例である。撮像方法や装置の状態によっては、静磁場不均一性や渦磁場による位相変化を低減するため、位相補正や他のデータ補正処理を施すのが望ましい。 Here, SE (Spin Echo) type EPI is shown as an example. For example, any imaging method can be used as long as it is an FE (Field Echo) type EPI or any other imaging method capable of obtaining a phase image. May be. It is desirable to use an imaging method that can obtain an image in as short a time as possible. For example, the imaging method described in “Proc. Of Annual Meeting, Society of Magnetic Resonance in Medicine, No. 2138, 1998” is an example. Depending on the imaging method and the state of the apparatus, it is desirable to perform phase correction or other data correction processing in order to reduce phase variations due to static magnetic field inhomogeneity or eddy magnetic fields.
図4は、図3に示したパルスシーケンスによる撮像の流れを示す図である。図4に示すように、図3に示したパルスシーケンスによれば、心電波形(図4に示すECG)のR波から一定の待ち時間(図4に示すTdelay)が経過した後に、パルスシーケンスの1ショット分のスキャンが行われる。 FIG. 4 is a diagram showing a flow of imaging by the pulse sequence shown in FIG. As shown in FIG. 4, according to the pulse sequence shown in FIG. 3, after a certain waiting time (Tdelay shown in FIG. 4) has elapsed from the R wave of the electrocardiogram waveform (ECG shown in FIG. 4), the pulse sequence Are scanned for one shot.
本実施例1では、フーリエ変換処理部22aは、上記で説明したように1ショットで得られたエコー信号のk空間データそれぞれに離散2次元フーリエ変換などの再構成処理を施すことによって、図4に示すように複数の画像S(ra,i)を再構成する。ここで、raは、画像上の位置を示すベクトルであり、iは、収集の順番を示す添え字である。各画像S(ra,i)では、流速に比例して画素の位相が変化している。
In the first embodiment, the Fourier
絶対値画像作成部22bは、フーリエ変換処理部22aによって再構成された画像から絶対値画像を作成する。具体的には、絶対値画像作成部22bは、フーリエ変換処理部22aによって再構成された複数の画像S(ra,i)それぞれについて、以下に示す式(1)で与えられる絶対値画像I(ra,i)を作成する。
The absolute value
I(ra,i)=abs{S(ra,i)} ・・・(1) I (r a , i) = abs {S (r a , i)} (1)
さらに、絶対値画像作成部22bは、作成した複数の絶対値画像I(ra,i)をもとに、以下に示す式(2)で与えられる平均絶対値画像Iavg(ra)を作成する。
Further, the absolute value
流速画像作成部22cは、フーリエ変換処理部22aによって再構成された画像から流速成分の分布を含む計算画像を作成する。本実施例では、流速画像作成部22cは、計算画像として、流速成分の分布を表す流速画像を作成する。具体的には、流速画像作成部22cは、フーリエ変換処理部22aによって再構成された複数の画像S(ra,i)それぞれについて、以下に示す式(3)で与えられる流速画像V(ra,i)を作成する。なお、式(3)において、k(venc)は、フローエンコード傾斜磁場のパルスの形状によって変化する比例係数である。
The flow velocity image creation unit 22c creates a calculation image including the distribution of flow velocity components from the image reconstructed by the Fourier
V(ra,i)=k(venc)・arg{S(ra,i)} ・・・(3) V (r a , i) = k (venc) · arg {S (r a , i)} (3)
さらに、流速画像作成部22cは、作成した複数の流速画像V(ra,i)をもとに、以下に示す式(4)で与えられる平均流速画像Vavg(ra)を作成する。 Furthermore, the flow velocity image creation unit 22c creates an average flow velocity image V avg (r a ) given by the following equation (4) based on the created plurality of flow velocity images V (r a , i).
流速分散画像作成部22dは、流速画像作成部22cによって作成された複数の流速画像を用いて、各流速画像の同一位置ごとに体液の流速変化を示す統計値を算出し、算出した統計値の分布を表す統計値画像を作成する。本実施例では、流速分散画像作成部22dは、統計値として、位相変化を利用して時系列に沿った流速成分の分散を算出し、統計値画像として、流速成分の分散の分布を表す流速分散画像を作成する。
The flow velocity dispersion
具体的には、流速分散画像作成部22dは、流速画像作成部22cによって作成された複数の流速画像V(ra,i)と平均流速画像Vavg(ra)とから、以下に示す式(5)で与えられる流速分散画像Vari,V(ra)を作成する。
Specifically, the flow velocity dispersion
重畳画像処理部22eは、流速分散画像作成部22dによって作成された流速分散画像における流速成分の分散の分布を表す重畳画像を作成する。具体的には、重畳画像処理部22eは、流速分散画像Vari,V(ra)における流速成分の分散の分布を、絶対値画像作成部22bによって作成された平均絶対値画像Iavg(ra)に重畳させる。
The superimposed
このとき、重畳画像処理部22eは、流速成分の分散値の値に応じて、流速成分の分散の分布の表示形態を変化させる。例えば、重畳画像処理部22eは、平均絶対値画像Iavg(ra)をグレースケールの濃淡画像とし、流速分散画像Vari,V(ra)における流速成分の分散の分布を以下に示す式(6)〜(8)で与えられるRGB(Red-Green-Blue)値で表す。なお、式(6)〜(8)において、P,redは赤(Red)の輝度を、P,greenは緑(Green)の輝度を、P,blueは青(Blue)の輝度を示している。また、kIおよびkVは、それぞれ所定の比例係数である。
At this time, the superimposed
P,red =kI・Iavg(ra)+kV・Vari,V(ra) ・・・(6)
P,green=kI・Iavg(ra) ・・・(7)
P,blue =kI・Iavg(ra) ・・・(8)
P, red = kI · I avg (r a ) + kV · V ari , V (r a ) (6)
P, green = kI · I avg (r a ) (7)
P, blue = kI · I avg (r a ) (8)
これにより、重畳画像において流速成分の分散の分布が赤色で表されるとともに、分散値の値に応じて画素ごとに輝度が変化するので、操作者は体液の流速変化を容易に把握することができるようになる。なお、ここでは、重畳画像処理部22eが、流速成分の分散の分布を赤色で表す場合について説明したが、例えば、緑色または青色で表すようにしてもよい。
Accordingly, the distribution of the dispersion of the flow velocity component is displayed in red in the superimposed image, and the luminance changes for each pixel according to the value of the dispersion value, so that the operator can easily grasp the flow velocity change of the body fluid. become able to. Here, although the case where the superimposed
図2にもどって、制御部26は、特に、画像表示制御部26aを有する。画像表示制御部26aは、重畳画像処理部22eによって作成された重畳画像を表示部25に表示する。
Returning to FIG. 2, the
次に、本実施例1に係るMRI装置100による画像処理の処理手順について説明する。図5は、本実施例1に係るMRI装置100による画像処理の処理手順を示すフローチャートであり、図6は、本実施例1に係るMRI装置100による画像作成の流れを示す図である。
Next, a processing procedure of image processing performed by the
図5に示すように、このMRI装置100では、制御部26が、入力部24を介して操作者からEPI法による撮像の撮像条件を受け付ける(ステップS101)。その後、操作者によって撮像開始が指示されると(ステップS102,Yes)、制御部26が、EPI法によるパルスシーケンスのシーケンス情報を生成し、生成したシーケンス情報をシーケンス制御部10に送信することによって、エコー信号に基づくk空間データを収集する(ステップS103)。
As shown in FIG. 5, in the
続いて、フーリエ変換処理部22aが、収集されたk空間データに離散2次元フーリエ変換などの再構成処理を施すことによって、図6の(1)に示すように、複数の画像S(ra,i)を再構成する(ステップS104)。
Subsequently, the Fourier
これに続いて、絶対値画像作成部22bが、フーリエ変換処理部22aによって再構成された複数の画像S(ra,i)それぞれについて、図6の(2)に示すように、絶対値画像I(ra,i)を作成する(ステップS105)。さらに、絶対値画像作成部22bは、作成した複数の絶対値画像I(ra,i)をもとに、図6の(3)に示すように、平均絶対値画像Iavg(ra)を作成する(ステップS106)。
Subsequently, as shown in (2) of FIG. 6, the absolute value
この一方で、流速画像作成部22cが、フーリエ変換処理部22aによって再構成された複数の画像S(ra,i)それぞれについて、図6の(4)に示すように、流速画像V(ra,i)を作成する(ステップS107)。さらに、流速画像作成部22cは、作成した複数の流速画像V(ra,i)をもとに、図6の(5)に示すように、平均流速画像Vavg(ra)を作成する(ステップS108)。
On the other hand, as shown in (4) of FIG. 6, the flow velocity image creation unit 22c applies a flow velocity image V (r) to each of the plurality of images S (r a , i) reconstructed by the Fourier
これに続いて、流速分散画像作成部22dが、流速画像作成部22cによって作成された複数の流速画像V(ra,i)と平均流速画像Vavg(ra)とから、図6の(6)に示すように、流速分散画像Vari,V(ra)を作成する(ステップS109)。
Subsequently, the flow velocity dispersion
そして、重畳画像処理部22eが、流速分散画像Vari,V(ra)の分散成分を色づけし(ステップS110)、図6の(7)に示すように、色づけした流速分散画像Vari,V(ra)を平均絶対値画像Iavg(ra)に重畳させ、画像表示制御部26aが、重畳画像処理部22eによって作成された重畳画像を表示部25に表示する(ステップS111)。
Then, the superimposed
上述してきたように、本実施例1では、流速画像作成部22cが、被検体内を流れる体液の流速成分が得られるEPIを複数回繰り返すことによって得られた複数の画像それぞれについて流速成分の分布を表す流速画像を作成する。また、流速分散画像作成部22dが、作成された複数の流速画像を用いて、各流速画像の同一位置ごとに時系列に沿った流速成分の分散を算出し、算出した流速成分の分散の分布を表す流速分散画像を作成する。そして、重畳画像処理部22eが、流速分散画像における流速成分の分散の分布を平均絶対値画像に重畳させ、画像表示制御部26aが、重畳画像を表示部25に表示する。したがって、本実施例1によれば、CSFの流速変化のように流速変化に周期性がない体液であっても、その流速変化の分布を忠実に画像化することができる。
As described above, in the first embodiment, the flow velocity image creation unit 22c distributes the flow velocity component for each of a plurality of images obtained by repeating EPI that obtains the flow velocity component of the body fluid flowing in the subject a plurality of times. Is created. Also, the flow velocity dispersion
なお、本実施例1では、重畳画像処理部22eが、平均絶対値画像Iavg(ra)をグレースケールの濃淡画像として用いる場合について説明したが、例えば、平均絶対値画像Iavg(ra)の替わりに他の撮像で得られた画像を用いてもよい。その場合、重畳画像処理部22eは、例えば、T1強調画像やT2強調画像などを用いる。
In the first embodiment, the case where the superimposed
また、本実施例1では、重畳画像処理部22eが、流速成分の分散の分布を1つの色で表す場合について説明したが、さらに、体液が流れる方向に応じて、統計値の分布の表示態様を変化させるようにしてもよい。その場合、重畳画像処理部22eは、例えば、超音波診断装置におけるパルスドップラーモードのように、血流の流れ方向を赤青の色彩情報で表し、血流の流速を明度で表し、血流の流速成分の分散を緑色で表すことによって、二次元血流画像を作成する。
In the first embodiment, the superimposed
また、本実施例1では、心電波形に同期させて撮像を繰り返すことによって得られた複数の画像を用いる場合について説明したが(図3を参照)、本発明はこれに限られるわけではない。図7は、本実施例1で用いられるパルスシーケンスの変形例を示す図である。例えば、図7に示すように、心電波形に同期させずに、一定の周期で撮像を複数回繰り返すことによって得られた複数の画像を用いてもよい。CSFの流速変化は心電同期との相関が低いため、心電波形に同期させずにエコー信号を収集しても支障は少ない。また、一定の周期で撮像を繰り返す場合には、同期撮像で必要な待ち時間(例えば、図3に示したTdelay)が不要となるので、より短時間で撮像を行うことができる。 In the first embodiment, the case where a plurality of images obtained by repeating imaging in synchronization with an electrocardiogram waveform is used (see FIG. 3), but the present invention is not limited to this. . FIG. 7 is a diagram illustrating a modification of the pulse sequence used in the first embodiment. For example, as shown in FIG. 7, a plurality of images obtained by repeating imaging a plurality of times at a constant cycle without being synchronized with an electrocardiographic waveform may be used. Since the change in the flow rate of the CSF has a low correlation with the electrocardiogram synchronization, there is little trouble even if the echo signal is collected without synchronizing with the electrocardiogram waveform. In addition, when imaging is repeated at a constant cycle, a waiting time necessary for synchronous imaging (for example, Tdelay shown in FIG. 3) is unnecessary, and thus imaging can be performed in a shorter time.
ところで、上記実施例1では、フローエンコード傾斜磁場を用いたEPI撮像を繰り返すことによって得られた複数の画像を用いて流速成分の分散の分布を画像化する場合について説明したが、本発明はこれに限られるわけではない。例えば、反転励起パルスを印加することによって撮像断面内にある一定の領域をラベリングし、当該領域から流れ出る体液の信号を検出して画像を再構成する撮像方法を用いてもよい(例えば、特開2001−252263号公報を参照)。 In the first embodiment, the case where the distribution of flow velocity component dispersion is imaged using a plurality of images obtained by repeating EPI imaging using a flow encode gradient magnetic field has been described. It is not limited to. For example, an imaging method may be used in which a certain region in the imaging cross section is labeled by applying an inversion excitation pulse, and a signal of body fluid flowing out from the region is detected to reconstruct an image (for example, JP 2001-252263).
そこで、以下では、かかる撮像方法を用いる場合を実施例2として説明する。なお、本実施例2では、この撮像方法を「Time−SLIP法」と呼ぶ。また、本実施例2に係るMRI装置の全体構成は図1に示したものと同様であるので、ここでは説明を省略し、本実施例2における画像再構成部、記憶部および制御部の詳細について説明する。 Therefore, hereinafter, a case where such an imaging method is used will be described as a second embodiment. In the second embodiment, this imaging method is referred to as “Time-SLIP method”. Further, since the overall configuration of the MRI apparatus according to the second embodiment is the same as that shown in FIG. 1, description thereof is omitted here, and details of the image reconstruction unit, storage unit, and control unit according to the second embodiment are omitted. Will be described.
図8は、本実施例2における画像再構成部、記憶部および制御部の構成を示す機能ブロック図である。なお、ここでは説明の便宜上、図2に示した各部と同様の役割を果たす機能部については、同一符号を付すこととしてその詳細な説明を省略する。 FIG. 8 is a functional block diagram illustrating configurations of an image reconstruction unit, a storage unit, and a control unit according to the second embodiment. Here, for convenience of explanation, functional units that play the same roles as the respective units shown in FIG. 2 are given the same reference numerals, and detailed descriptions thereof are omitted.
図8に示すように、本実施例2における画像再構成部122は、特に、フーリエ変換処理部122a、絶対値画像作成部122b、信号強度分散画像作成部122d、重畳画像処理部122eを有する。
As shown in FIG. 8, the
フーリエ変換処理部122aは、k空間データ記憶部23aに記憶されたk空間データに離散2次元フーリエ変換などの再構成処理を施すことによって画像を再構成する。
The Fourier transform processing unit 122a reconstructs an image by performing reconstruction processing such as discrete two-dimensional Fourier transform on k-space data stored in the k-space
ここで、本実施例2では、被検体内を流れる体液の流速成分が得られる撮像法として、Time−SLIP法が用いられる。図9は、本実施例2で用いられるTime−SLIP法のパルスシーケンスを示す図であり、図10は、本実施例2におけるTime−SLIP法によって行われるラベリングの位置の一例を示す図である。 Here, in the second embodiment, the Time-SLIP method is used as an imaging method for obtaining the flow velocity component of the body fluid flowing in the subject. FIG. 9 is a diagram illustrating a pulse sequence of the Time-SLIP method used in the second embodiment, and FIG. 10 is a diagram illustrating an example of a labeling position performed by the Time-SLIP method in the second embodiment. .
図9に示すように、本実施例2で用いられるTime−SLIP法では、心電波形(図9に示すECG)のR波から一定の待ち時間(図9に示すTDtag1)が経過した後に、例えば、図10に示すように、撮像断面内にある一定の領域を反転励起させるIRパルス(図9に示すP21)が印加される。このとき、励起位置に合わせて励起周波数(図9に示すΔf)、選択励起の傾斜磁化(図9に示すGss、Gro、Gpe)を変化させる。 As shown in FIG. 9, in the Time-SLIP method used in the second embodiment, after a certain waiting time (TD tag1 shown in FIG. 9) has elapsed from the R wave of the electrocardiogram waveform (ECG shown in FIG. 9). For example, as shown in FIG. 10, an IR pulse (P21 shown in FIG. 9) is applied to invert and excite a certain region in the imaging section. At this time, the excitation frequency (Δf shown in FIG. 9) and the gradient magnetization of selective excitation (Gss, Gro, Gpe shown in FIG. 9) are changed in accordance with the excitation position.
その後、一定時間(図9に示すTI1)が経過した後に、1画像を再構成するのに必要なエコー信号を収集する(図9に示すP22)。なお、図9は、高速スピンエコー法を用いてエコー信号を収集する場合を示している。本実施例2で用いられるTime−SLIP法では、この収集を複数回繰り返すことによって複数画像分のエコー信号を収集する。 After that, after a predetermined time (TI 1 shown in FIG. 9) has elapsed, echo signals necessary to reconstruct one image are collected (P22 shown in FIG. 9). FIG. 9 shows a case where echo signals are collected using the fast spin echo method. In the Time-SLIP method used in the second embodiment, echo signals for a plurality of images are collected by repeating this collection a plurality of times.
本実施例2では、フーリエ変換処理部122aは、上記で説明したように収集されたエコー信号のk空間データそれぞれに離散2次元フーリエ変換などの再構成処理を施すことによって、複数の画像S(ra,i)を再構成する。 In the second embodiment, the Fourier transform processing unit 122a performs a reconstruction process such as a discrete two-dimensional Fourier transform on each of the k-space data of the echo signals collected as described above, so that a plurality of images S ( Reconstruct r a , i).
絶対値画像作成部122bは、フーリエ変換処理部122aによって再構成された画像から信号強度の分布を表す絶対値画像を作成する。具体的には、絶対値画像作成部122bは、フーリエ変換処理部122aによって再構成された複数の画像S(ra,i)それぞれについて、以下に示す式(9)で与えられる絶対値画像I(ra,i)を作成する。
The absolute value
I(ra,i)=abs{S(ra,i)} ・・・(9) I (r a , i) = abs {S (r a , i)} (9)
図11は、本実施例2における絶対値画像作成部122bによって作成される絶対値画像I(ra,i)の一例を示す図である。図11に示すように、ラベリングされた領域から流れ出たCSFが画像化されるが、各絶対値画像のもとになるエコー信号が収集された際のCSFの流速がそれぞれ異なるため、各絶対値画像において、ラベリングされた領域から流れ出たCSFによって信号強度が変化している領域は異なる。
FIG. 11 is a diagram illustrating an example of the absolute value image I (r a , i) created by the absolute value
さらに、絶対値画像作成部122bは、作成した複数の絶対値画像I(ra,i)をもとに、以下に示す式(10)で与えられる平均絶対値画像Iavg(ra)を作成する。
Furthermore, the absolute value
図8にもどって、信号強度分散画像作成部122dは、絶対値画像作成部122bによって作成された複数の絶対値画像を用いて、各絶対値画像の同一位置ごとに体液の流速変化を示す統計値を算出し、算出した統計値の分布を表す統計値画像を作成する。本実施例では、信号強度分散画像作成部122dは、統計値として、体液の信号値に基づいて時系列に沿った信号強度の分散を算出し、統計値画像として、信号強度の分散の分布を表す信号強度分散画像を作成する。
Returning to FIG. 8, the signal intensity dispersion
具体的には、信号強度分散画像作成部122dは、絶対値画像作成部122bによって作成された複数の絶対値画像I(ra,i)と平均絶対値画像Iavg(ra)とから、以下に示す式(11)で与えられる信号強度分散画像Vari,I(ra)を作成する。
Specifically, the signal intensity dispersion
重畳画像処理部122eは、信号強度分散画像作成部122dによって作成された信号強度分散画像における信号強度の分散の分布を表す重畳画像を作成する。具体的には、重畳画像処理部122eは、信号強度分散画像Vari,I(ra)における信号強度の分散の分布を、絶対値画像作成部122bによって作成された平均絶対値画像Iavg(ra)に重畳させる。
The superimposed
このとき、重畳画像処理部122eは、信号強度の分散値の値に応じて、信号強度の分散の分布の表示形態を変化させる。例えば、重畳画像処理部122eは、平均絶対値画像Iavg(ra)をグレースケールの濃淡画像とし、信号強度分散画像Vari,I(ra)における信号強度の分散の分布を以下に示す式(12)〜(14)で与えられるRGB(Red-Green-Blue)値で表す。なお、式(12)〜(14)において、P,redは赤(Red)の輝度を、P,greenは緑(Green)の輝度を、P,blueは青(Blue)の輝度を示している。また、kIおよびkVは、それぞれ所定の比例係数である。
At this time, the superimposed
P,red =kI・Iavg(ra)+kV・Vari,I(ra) ・・・(12)
P,green=kI・Iavg(ra) ・・・(13)
P,blue =kI・Iavg(ra) ・・・(14)
P, red = kI · I avg (r a ) + kV · V ari , I (r a ) (12)
P, green = kI · I avg (r a ) (13)
P, blue = kI · I avg (r a ) (14)
これにより、重畳画像において信号強度の分散の分布が赤色で表されるとともに、分散値の値に応じて画素ごとに輝度が変化するので、操作者は体液の流動による信号強度の変化を容易に把握することができるようになる。なお、ここでは、重畳画像処理部122eが、信号強度の分散の分布を赤色で表す場合について説明したが、例えば、緑色または青色で表すようにしてもよい。
As a result, the distribution of the signal intensity dispersion is displayed in red in the superimposed image, and the luminance changes for each pixel according to the value of the dispersion value, so that the operator can easily change the signal intensity due to the flow of the body fluid. It becomes possible to grasp. Although the case where the superimposed
次に、本実施例2に係るMRI装置100による画像処理の処理手順について説明する。図12は、本実施例2に係るMRI装置100による画像処理の処理手順を示すフローチャートであり、図13は、本実施例2に係るMRI装置100による画像作成の流れを示す図である。
Next, a processing procedure of image processing by the
図12に示すように、このMRI装置100では、制御部26が、入力部24を介して操作者からTime−SLIP法による撮像の撮像条件を受け付ける(ステップS201)。その後、操作者によって撮像開始が指示されると(ステップS202,Yes)、制御部26が、Time−SLIP法によるパルスシーケンスのシーケンス情報を生成し、生成したシーケンス情報をシーケンス制御部10に送信することによって、エコー信号に基づくk空間データを収集する(ステップS203)。
As shown in FIG. 12, in this
続いて、フーリエ変換処理部122aが、収集されたk空間データに離散2次元フーリエ変換などの再構成処理を施すことによって、図13の(1)に示すように、複数の画像S(ra,i)を再構成する(ステップS204)。 Subsequently, the Fourier transform processing unit 122a performs a reconstruction process such as a discrete two-dimensional Fourier transform on the collected k-space data, so that a plurality of images S (r a are shown in FIG. , I) is reconfigured (step S204).
これに続いて、絶対値画像作成部122bが、フーリエ変換処理部122aによって再構成された複数の画像S(ra,i)それぞれについて、図13の(2)に示すように、絶対値画像I(ra,i)を作成する(ステップS205)。さらに、絶対値画像作成部122bは、作成した複数の絶対値画像I(ra,i)をもとに、図13の(3)に示すように、平均絶対値画像Iavg(ra)を作成する(ステップS206)。
Subsequently, as shown in (2) of FIG. 13, the absolute value
その後、信号強度分散画像作成部122dが、絶対値画像作成部122bによって作成された複数の絶対値画像I(ra,i)と平均絶対値画像Iavg(ra)とから、図13の(4)に示すように、信号強度分散画像Vari,I(ra)を作成する(ステップS207)。
After that, the signal intensity dispersion
そして、重畳画像処理部122eが、信号強度分散画像Vari,I(ra)の分散成分を色づけし(ステップS208)、図13の(5)に示すように、色づけした信号強度分散画像Vari,I(ra)を平均絶対値画像Iavg(ra)に重畳させ、画像表示制御部26aが、重畳画像処理部122eによって作成された重畳画像を表示部25に表示する(ステップS209)。
Then, the superimposed
上述してきたように、本実施例2では、絶対値画像作成部122bが、ラベリングされた領域から流れ出る体液の信号を検出して画像を再構成するTime―SLIP法によって得られた複数の画像それぞれについて信号強度の分布を表す絶対値画像を作成する。また、信号強度分散画像作成部122dが、作成された複数の絶対値画像を用いて、各絶対値画像の同一位置ごとに時系列に沿った信号強度の分散を算出し、算出した信号強度の分散の分布を表す信号強度分散画像を作成する。そして、重畳画像処理部122eが、信号強度分散画像における信号強度の分散の分布を平均絶対値画像に重畳させ、画像表示制御部26aが、重畳画像を表示部25に表示する。したがって、本実施例2によれば、実施例1と同様、CSFの流速変化のように流速変化に周期性がない体液であっても、その流速変化の分布を忠実に画像化することができる。
As described above, in the second embodiment, each of the plurality of images obtained by the Time-SLIP method in which the absolute value
なお、上記実施例1および2では、統計値として、流速成分あるいは信号強度の分散を算出する場合について説明したが、本発明はこれに限られるわけではない。例えば、統計値として、複数の計算画像上の同一位置における流速成分あるいは信号強度の最大値を算出し、つまり最大値画像Vmax(ra)を作成し、算出した最大値に応じて流速成分を色づけするようにしてもよい。 In the first and second embodiments, the case where the flow velocity component or the variance of the signal intensity is calculated as the statistical value has been described. However, the present invention is not limited to this. For example, the flow velocity component or the maximum value of the signal intensity at the same position on a plurality of calculation images is calculated as the statistical value, that is, the maximum value image V max (r a ) is created, and the flow velocity component is generated according to the calculated maximum value. You may make it color.
または、統計値として、各計算画像における同一位置ごとに大きい順に流速成分あるいは信号強度を並べた場合に、全体に対して上位から所定の割合(例えば、上位10%など)に入る流速成分あるいは信号強度の平均値を算出し、算出した平均値に応じて流速成分を色づけするようにしてもよい。この場合、具体的には、計算画像から同一ピクセルの画素値のヒストグラムを求め、全体に対して上位から所定の割合に入る画素値の平均を算出する。 Alternatively, as a statistical value, when flow velocity components or signal strengths are arranged in descending order at the same position in each calculation image, flow velocity components or signals that fall within a predetermined ratio (for example, upper 10%) from the top with respect to the whole. An average value of intensity may be calculated, and the flow velocity component may be colored according to the calculated average value. In this case, specifically, a histogram of pixel values of the same pixel is obtained from the calculated image, and an average of pixel values that fall within a predetermined ratio with respect to the whole is calculated.
なお、上記実施例1および2では、MRI装置によって撮像された画像を用いる場合について説明したが、本発明はこれに限定されるものではなく、例えばX線CT装置やX線診断装置など、他の診断装置によって撮像された画像を用いる場合にも同様に適用することができる。すなわち、被検体内を流れる体液の流速成分が得られる撮像を複数回繰り返すことによって得られた複数の画像を用いる場合には、本発明を同様に適用することが可能である。 In the first and second embodiments, the case where an image picked up by an MRI apparatus is used has been described. However, the present invention is not limited to this, for example, an X-ray CT apparatus, an X-ray diagnostic apparatus, and the like. The same can be applied to the case of using an image captured by the diagnostic apparatus. That is, the present invention can be similarly applied when using a plurality of images obtained by repeating a plurality of times of imaging that obtains the flow velocity component of the body fluid flowing in the subject.
また、上記実施例1および2に係る画像処理方法は、十分な繰り返し回数の撮像が前提となるが、流速の周期的な変化を前提としないため、CSFのように心電同期とはあまり相関のない動きの画像上での分布を簡便に、かつ安定して得ることができる。また、1ショットの撮像を使用しているため、通常のフェイズシフト法のように平均加算効果によって流速変化が過小評価されるのを防ぐことができる。また、同期撮像を行わない場合には、不規則な動きに関する俯瞰的な情報を比較的簡便に得ることができる。 In addition, the image processing methods according to the first and second embodiments are premised on imaging with a sufficient number of repetitions, but are not premised on periodic changes in the flow velocity, and thus are not very correlated with electrocardiographic synchronization like CSF. It is possible to easily and stably obtain a distribution on an image having no motion. In addition, since one-shot imaging is used, it is possible to prevent an underestimation of a change in flow velocity due to an average addition effect as in a normal phase shift method. In addition, when synchronous imaging is not performed, it is possible to relatively easily obtain bird's-eye information regarding irregular movement.
1 静磁場磁石
2 傾斜磁場コイル
3 傾斜磁場電源
4 寝台
4a 天板
5 寝台制御部
6 送信RFコイル
7 送信部
8 受信RFコイル
9 受信部
10 シーケンス制御部
15 架台部
20 計算機システム
21 インタフェース部
22,122 画像再構成部
22a,122a フーリエ変換処理部
22b,122b 絶対値画像作成部
22c 流速画像作成部
22d 流速分散画像作成部
122d 信号強度分散画像作成部
22e,122e 重畳画像処理部
23 記憶部
23a k空間データ記憶部
23b 画像データ記憶部
24 入力部
25 表示部
26 制御部
26a 画像表示制御部
100 MRI装置(磁気共鳴イメージング装置)
DESCRIPTION OF
Claims (6)
前記被検体内を流れる脳脊髄液を得るための撮像を複数回繰り返すことによって脳脊髄液に関する複数の再構成画像を作成する画像作成手段と、
前記画像作成手段によって作成された複数の再構成画像を用いて、各再構成画像における同一位置ごとに前記脳脊髄液の流速成分又は信号強度を大きい順に並べた場合に全体に対して上位から所定の割合に入る流速成分又は信号強度の統計値を算出し、当該算出された統計値に基づく計算画像を作成する計算画像作成手段と、
を備えたことを特徴とする磁気共鳴イメージング装置。 A magnetic resonance imaging apparatus that applies an excitation pulse to an imaging region of a subject in synchronization with a trigger signal, collects echo signals related to spins included in the imaging region, and generates a magnetic resonance image based on the collected echo signals There,
Image creating means for creating a plurality of reconstructed images related to cerebrospinal fluid by repeating imaging for obtaining cerebrospinal fluid flowing in the subject a plurality of times;
When a plurality of reconstructed images created by the image creating means are used and the flow velocity components or signal strengths of the cerebrospinal fluid are arranged in descending order at the same position in each reconstructed image , the whole is predetermined from the top. A calculation image creation means for calculating a statistical value of a flow velocity component or a signal intensity that falls within a ratio of the above, and creating a calculation image based on the calculated statistical value;
A magnetic resonance imaging apparatus comprising:
前記画像取得手段によって取得された複数の再構成画像を用いて、各再構成画像における同一位置ごとに前記脳脊髄液の流速成分又は信号強度を大きい順に並べた場合に全体に対して上位から所定の割合に入る流速成分又は信号強度の統計値を算出し、当該算出された統計値に基づく計算画像を作成する計算画像作成手段と、
を備えたことを特徴とする画像処理装置。 Image acquisition means for acquiring a plurality of reconstructed images related to cerebrospinal fluid created by repeating imaging for obtaining cerebrospinal fluid flowing in the subject a plurality of times;
When a plurality of reconstructed images acquired by the image acquiring means are used and the flow velocity components or signal intensities of the cerebrospinal fluid are arranged in descending order at the same position in each reconstructed image , the whole is predetermined from the top. A calculation image creation means for calculating a statistical value of a flow velocity component or a signal intensity that falls within a ratio of the above, and creating a calculation image based on the calculated statistical value;
An image processing apparatus comprising:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011133906A JP5283092B2 (en) | 2011-06-16 | 2011-06-16 | Magnetic resonance imaging apparatus and image processing apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011133906A JP5283092B2 (en) | 2011-06-16 | 2011-06-16 | Magnetic resonance imaging apparatus and image processing apparatus |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2008258963A Division JP4771490B2 (en) | 2008-10-03 | 2008-10-03 | Magnetic resonance imaging system |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2011173020A JP2011173020A (en) | 2011-09-08 |
JP5283092B2 true JP5283092B2 (en) | 2013-09-04 |
Family
ID=44686413
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2011133906A Active JP5283092B2 (en) | 2011-06-16 | 2011-06-16 | Magnetic resonance imaging apparatus and image processing apparatus |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5283092B2 (en) |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3018196B2 (en) * | 1990-07-04 | 2000-03-13 | 名古屋大学長 | Image display device |
JP3366398B2 (en) * | 1993-10-25 | 2003-01-14 | 株式会社日立メディコ | MRI equipment |
JPH10234695A (en) * | 1997-02-21 | 1998-09-08 | Hitachi Medical Corp | Magnetic resonant imaging system |
EP1182970B1 (en) * | 1999-12-15 | 2006-05-10 | Koninklijke Philips Electronics N.V. | Diagnostic imaging system with ultrasound probe |
JP4594482B2 (en) * | 2000-03-14 | 2010-12-08 | 株式会社東芝 | Magnetic resonance imaging method and magnetic resonance imaging apparatus |
JP4622715B2 (en) * | 2005-07-15 | 2011-02-02 | 株式会社日立製作所 | Image data analysis method and system |
JP5002222B2 (en) * | 2006-09-13 | 2012-08-15 | 学校法人東海大学 | Magnetic resonance imaging system |
-
2011
- 2011-06-16 JP JP2011133906A patent/JP5283092B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
JP2011173020A (en) | 2011-09-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4771490B2 (en) | Magnetic resonance imaging system | |
JP5105848B2 (en) | Magnetic resonance imaging apparatus and imaging condition setting method in magnetic resonance imaging apparatus | |
US9305376B2 (en) | Magnetic resonance imaging apparatus and method of acquiring functional image | |
US9301704B2 (en) | Magnetic resonance imaging system for non-contrast MRA and magnetic resonance signal acquisition method employed by the same | |
JP5433436B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
US10168405B2 (en) | Method and apparatus for quantifying properties of an object through magnetic resonance imaging (MRI) | |
JP2008067857A (en) | Magnetic resonance imaging apparatus | |
JP5372015B2 (en) | Magnetic resonance imaging apparatus and synchronous imaging method | |
KR20160011012A (en) | Magnetic resonance imaging device and control method thereof | |
JP6266325B2 (en) | Image analysis apparatus, image analysis program, and magnetic resonance imaging apparatus | |
JP5283092B2 (en) | Magnetic resonance imaging apparatus and image processing apparatus | |
JP4862069B2 (en) | Magnetic resonance diagnostic equipment | |
JP5319835B2 (en) | Magnetic resonance imaging system | |
JP2008148918A (en) | Mri system and its control method | |
JP2015002834A (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
JP2011143282A (en) | Mri apparatus | |
JP5689595B2 (en) | Magnetic resonance imaging system | |
JP5984274B2 (en) | Magnetic resonance imaging apparatus and display processing system | |
JP5514293B2 (en) | Magnetic resonance imaging system | |
JP5133711B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance image generation method | |
JP2015054098A (en) | Magnetic resonance imaging apparatus | |
JP2013063191A (en) | Magnetic resonance imaging apparatus and roi analyzing method | |
JP2011183212A (en) | Magnetic resonance imaging apparatus | |
JP2011212453A (en) | Magnetic resonance imaging apparatus |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20110616 |
|
RD01 | Notification of change of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7426 Effective date: 20110701 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20110701 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20121030 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20121227 |
|
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: 20130430 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20130521 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5283092 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |