JP3353826B2 - 磁場不均一測定装置、位相補正装置および磁気共鳴撮像装置 - Google Patents
磁場不均一測定装置、位相補正装置および磁気共鳴撮像装置Info
- Publication number
- JP3353826B2 JP3353826B2 JP17765899A JP17765899A JP3353826B2 JP 3353826 B2 JP3353826 B2 JP 3353826B2 JP 17765899 A JP17765899 A JP 17765899A JP 17765899 A JP17765899 A JP 17765899A JP 3353826 B2 JP3353826 B2 JP 3353826B2
- Authority
- JP
- Japan
- Prior art keywords
- phase
- image
- magnetic field
- pixel data
- unit
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/24—Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56563—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/24—Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
- G01R33/243—Spatial mapping of the polarizing magnetic field
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/4828—Resolving the MR signals of different chemical species, e.g. water-fat imaging
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- General Health & Medical Sciences (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Description
法および装置、位相補正方法および装置、並びに、磁気
共鳴撮像装置に関し、特に、磁気共鳴撮像した画像のピ
クセルデータの位相に基づいて撮像空間の磁場不均一を
測定する方法および装置、測定した磁場不均一に基づい
てピクセルデータの位相を補正する方法および装置、並
びに、位相を補正したピクセルデータに基づいて水と脂
肪を分離した画像を得る磁気共鳴撮像装置に関する。
する空間に静磁場を形成し、静磁場空間に勾配磁場と高
周波磁場を形成し、そのような磁場空間から発生する磁
気共鳴信号に基づいて画像を生成(再構成)するように
なっている。脂肪の磁気共鳴信号は、ケミカルシフト
(chemical shift)により水の磁気共鳴
信号とは周波数が異なるので、周波数の相違に基づく位
相差を利用して水と脂肪を別々に画像化することが行わ
れる。
の影響を受けるので、磁場不均一に影響されずに水と脂
肪を別々に画像化するために、いわゆるディクソン(D
ixon)法による撮像が行われる。この技法は、水の
磁気共鳴信号の位相(以下、単に水の位相という)と脂
肪の磁気共鳴信号の位相(以下、単に脂肪の位相とい
う)が同相および逆相となる画像をそれぞれ撮像し、両
画像の和によって水画像を求め、差によって脂肪画像を
求めるものである。
枚の画像を得るのに通常の2枚分の撮像を必要とするの
で、能率が悪いという問題があった。また、水に対する
脂肪の位相の変化速度は静磁場強度が低くなるほど遅く
なるので、例えばグラディエントエコー(gradie
nt echo)等を利用して水と脂肪が同相となる信
号を得るためには、長いTE(エコータイム:echo
time)を必要とし、このため信号の減衰が大きく
なり、例えば0.2T程度の低磁場においてはDixo
n法は実用に耐えないという問題があった。
されたもので、その目的は、水と脂肪を分離した撮像を
能率良く行うための磁場不均一測定方法および装置、磁
場不均一に基づく位相補正方法および装置、並びに、水
と脂肪を分離した撮像を能率良く行う磁気共鳴撮像装置
を実現することである。
する第1の観点での発明は、静磁場空間中の対象物体に
ついて磁気共鳴を利用して水のピクセルデータと脂肪の
ピクセルデータが2π/n(n≧2)の位相差を持つ断
層像を撮像し、前記断層像のピクセルデータの位相をn
倍し、前記n倍した位相が±πの範囲を超えることによ
る折り返しを補正し、前記折り返しを補正した位相を1
/nして前記静磁場空間の磁場不均一を表す位相を求め
ることを特徴とする磁場不均一測定方法である。
の発明は、静磁場空間中の対象物体について磁気共鳴を
利用して水のピクセルデータと脂肪のピクセルデータが
2π/n(n≧2)の位相差を持つ断層像を撮像する撮
像手段と、前記断層像のピクセルデータの位相をn倍す
る位相逓倍手段と、前記n倍した位相が±πの範囲を超
えることによる折り返しを補正する折り返し補正手段
と、前記折り返しを補正した位相を1/nして前記静磁
場空間の磁場不均一を表す位相を求める磁場不均一測定
手段とを具備することを特徴とする磁場不均一測定装置
である。
の発明は、静磁場空間中の対象物体について磁気共鳴を
利用して水のピクセルデータと脂肪のピクセルデータが
2π/n(n≧2)の位相差を持つ断層像を撮像し、前
記断層像のピクセルデータの位相をn倍し、前記n倍し
た位相が±πの範囲を超えることによる折り返しを補正
し、前記折り返しを補正した位相を1/nして前記静磁
場空間の磁場不均一を表す位相を求め、前記求めた位相
を用いて前記断層像のピクセルデータの位相を補正する
ことを特徴とする位相補正方法である。
の発明は、静磁場空間中の対象物体について磁気共鳴を
利用して水のピクセルデータと脂肪のピクセルデータが
2π/n(n≧2)の位相差を持つ断層像を撮像する撮
像手段と、前記断層像のピクセルデータの位相をn倍す
る位相逓倍手段と、前記n倍した位相が±πの範囲を超
えることによる折り返しを補正する折り返し補正手段
と、前記折り返しを補正した位相を1/nして前記静磁
場空間の磁場不均一を表す位相を求める磁場不均一測定
手段と、前記求めた位相を用いて前記断層像のピクセル
データにつき位相補正を行う位相補正手段とを具備する
ことを特徴とする位相補正装置である。
発明は、静磁場空間中の対象物体について磁気共鳴を利
用して水のピクセルデータと脂肪のピクセルデータが2
π/n(n≧2)の位相差を持つ断層像を撮像する撮像
手段と、前記断層像のピクセルデータの位相をn倍する
位相逓倍手段と、前記n倍した位相が±πの範囲を超え
ることによる折り返しを補正する折り返し補正手段と、
前記折り返しを補正した位相を1/nして前記静磁場空
間の磁場不均一を表す位相を求める磁場不均一測定手段
と、前記求めた位相を用いて前記断層像のピクセルデー
タにつき位相補正を行う位相補正手段と、前記位相補正
したピクセルデータの位相差を利用して水画像と脂肪画
像を別々に生成する画像生成手段とを具備することを特
徴とする磁気共鳴撮像装置である。
の発明は、前記位相補正したピクセルデータの位相差を
前記ピクセルデータの位相のヒストグラムに基づいて修
正する位相修正手段を有し、前記画像生成手段は前記修
正後の位相差を利用して水画像と脂肪画を別々に生成す
ることを特徴とする(5)に記載の磁気共鳴撮像装置で
ある。
発明は、静磁場空間中の対象物体について磁気共鳴を利
用して水のピクセルデータと脂肪のピクセルデータが2
π/n(n≧2)の位相差を持つ断層像を撮像し、前記
断層像のピクセルデータの位相をn倍し、前記n倍した
位相が±πの範囲を超えることによる折り返しを補正
し、前記折り返しを補正した位相を1/nして前記静磁
場空間の磁場不均一を表す位相を求め、前記求めた位相
を用いて前記断層像のピクセルデータにつき位相を補正
し、前記位相補正したピクセルデータの位相差を利用し
て水画像と脂肪画像を別々に生成することを特徴とする
磁気共鳴撮像方法である。
発明は、前記位相補正したピクセルデータの位相差を前
記ピクセルデータの位相のヒストグラムに基づいて修正
し、前記修正後の位相差を利用して水画像と脂肪画を別
々に生成することを特徴とする(7)に記載の磁気共鳴
撮像方法である。
相をn倍して水と脂肪を同相にし、±πの範囲を超える
分の折り返しを補正して1/nし、磁場不均一に対応す
る位相分布を得る。この位相分布でピクセルデータの位
相を補正して磁場不均一の影響を除く。磁場不均一の影
響を除いたピクセルデータから、位相差を利用して水画
像と脂肪画像を分離する。
施の形態を詳細に説明する。なお、本発明は実施の形態
に限定されるものではない。図1に磁気共鳴撮像装置の
ブロック(block)図を示す。本装置は本発明の実
施の形態の一例である。本装置の構成によって、本発明
の装置に関する実施の形態の一例が示される。本装置の
動作によって、本発明の方法に関する実施の形態の一例
が示される。
ステム(magnet system)100を有す
る。マグネットシステム100は主磁場コイル部10
2、勾配コイル部106およびRF(radio fr
equency)コイル部108を有する。これら各コ
イル部は概ね円筒状の外形を有し、互いに同軸的に配置
されている。マグネットシステム100の内部空間に、
撮像対象300がクレードル(cradle)500に
搭載されて図示しない搬送手段により搬入および搬出さ
れる。撮像対象300は、本発明における対象物体の実
施の形態の一例である。
ム100の内部空間に静磁場を形成する。静磁場の方向
は概ね撮像対象300の体軸の方向に平行である。すな
わちいわゆる水平磁場を形成する。主磁場コイル部10
2は例えば超伝導コイルを用いて構成される。なお、超
伝導コイルに限らず常伝導コイル等を用いて構成しても
良いのはもちろんである。
持たせるための勾配磁場を生じる。発生する勾配磁場
は、スライス(slice)勾配磁場、リードアウト
(read out)勾配磁場およびフェーズエンコー
ド(phase encode)勾配磁場の3種であ
り、これら3種類の勾配磁場に対応して勾配コイル部1
06は図示しない3系統の勾配コイルを有する。
象300の体内のスピン(spin)を励起するための
高周波磁場を形成する。以下、高周波磁場を形成するこ
とをRF励起信号の送信という。RFコイル部108
は、また、励起されたスピンが生じる電磁波すなわち磁
気共鳴信号を受信する。RFコイル部108は図示しな
い送信用のコイルおよび受信用のコイルを有する。送信
用のコイルおよび受信用のコイルは、同じコイルを兼用
するかあるいはそれぞれ専用のコイルを用いる。
が接続されている。勾配駆動部130は勾配コイル部1
06に駆動信号を与えて勾配磁場を発生させる。勾配駆
動部130は、勾配コイル部106における3系統の勾
配コイルに対応して、図示しない3系統の駆動回路を有
する。
が接続されている。RF駆動部140はRFコイル部1
08に駆動信号を与えてRF励起信号を送信し、撮像対
象300の体内のスピンを励起する。RFコイル部10
8には、また、データ収集部150が接続されている。
データ収集部150はRFコイル部108が受信した受
信信号を取り込み、それをディジタルデータ(digi
tal data)として収集する。
びデータ収集部150には制御部160が接続されてい
る。制御部160は、勾配駆動部130ないしデータ収
集部150をそれぞれ制御する。
部170に接続されている。データ処理部170は、デ
ータ収集部150から取り込んだデータを図示しないメ
モリ(memory)に記憶する。メモリ内にはデータ
空間が形成される。データ空間は2次元フーリエ(Fo
urier)空間を構成する。データ処理部170は、
これら2次元フーリエ空間のデータを2次元逆フーリエ
変換して撮像対象300の画像を再構成する。
されている。データ処理部170は制御部160の上位
にあってそれを統括する。データ処理部170には、表
示部180および操作部190が接続されている。表示
部180は、データ処理部170から出力される再構成
画像および各種の情報を表示する。操作部190は、操
作者によって操作され、各種の指令や情報等をデータ処
理部170に入力する。
示す。本装置は本発明の実施の形態の一例である。本装
置の構成によって、本発明の装置に関する実施の形態の
一例が示される。本装置の動作によって、本発明の方法
に関する実施の形態の一例が示される。
異なるマグネットシステム100’を有する。マグネッ
トシステム100’以外は図1に示した装置と同様な構
成になっており、同様な部分に同一の符号を付して説明
を省略する。
ネット部102’、勾配コイル部106’およびRFコ
イル部108’を有する。これら主磁場マグネット部1
02’および各コイル部は、いずれも空間を挟んで互い
に対向する1対のものからなる。また、いずれも概ね円
盤状の外形を有し中心軸を共有して配置されている。マ
グネットシステム100’の内部空間に、撮像対象30
0がクレードル500に搭載されて図示しない搬送手段
により搬入および搬出される。
システム100’の内部空間に静磁場を形成する。静磁
場の方向は概ね撮像対象300の体軸方向と直交する。
すなわちいわゆる垂直磁場を形成する。主磁場マグネッ
ト部102’は例えば永久磁石等を用いて構成される。
なお、永久磁石に限らず超伝導電磁石あるいは常伝導電
磁石等を用いて構成しても良いのはもちろんである。
を持たせるための勾配磁場を生じる。発生する勾配磁場
は、スライス勾配磁場、リードアウト勾配磁場およびフ
ェーズエンコード勾配磁場の3種であり、これら3種類
の勾配磁場に対応して勾配コイル部106’は図示しな
い3系統の勾配コイルを有する。
対象300の体内のスピンを励起するためのRF励起信
号を送信する。RFコイル部108’は、また、励起さ
れたスピンが生じる磁気共鳴信号を受信する。RFコイ
ル部108’は図示しない送信用のコイルおよび受信用
のコイルを有する。送信用のコイルおよび受信用のコイ
ルは、同じコイルを兼用するかあるいはそれぞれ専用の
コイルを用いる。
ケンス(pulse sequence)の一例を示
す。このパルスシーケンスは、スピンエコー(SE:s
pinecho)法のパルスシーケンスである。
起用の90°パルスおよび180°パルスのシーケンス
であり、(2)、(3)、(4)および(5)は、同じ
くそれぞれ、スライス勾配Gs、リードアウト勾配G
r、フェーズエンコード勾配GpおよびスピンエコーM
Rのシーケンスである。なお、90°パルスおよび18
0°パルスはそれぞれ中心信号で代表する。パルスシー
ケンスは時間軸tに沿って左から右に進行する。
ピンの90°励起が行われる。このときスライス勾配G
sが印加され所定のスライスについての選択励起が行わ
れる。90°励起から所定の時間後に、180°パルス
による180°励起すなわちスピン反転が行われる。こ
のときもスライス勾配Gsが印加され、同じスライスに
ついての選択的反転が行われる。
ードアウト勾配Grおよびフェーズエンコード勾配Gp
が印加される。リードアウト勾配Grによりスピンのデ
ィフェーズ(dephase)が行われる。フェーズエ
ンコード勾配Gpによりスピンのフェーズエンコードが
行われる。
ピンをリフェーズ(rephase)してスピンエコー
MRを発生させる。スピンエコーMRはデータ収集部1
50によりビューデータ(view data)として
収集される。このようなパスルシーケンスが周期TR
(repetition time)で128〜256
回繰り返される。繰り返しのたびにフェーズエンコード
勾配Gpを変更し、毎回異なるフェーズエンコードを行
う。これによって、128〜256ビューのビューデー
タが得られる。
対称的な波形を持つRF信号となる。中心エコーは90
°励起からTE(echo time)後に生じる。時
間TEを適切に選ぶことにより、水のエコーと脂肪のエ
コーの位相差を2π/n(n≧2)とすることができ
る。例えば、n=4としたとき、位相差をπ/2にする
TEは静磁場強度が0.2Tの場合で2τ+8.6ms
または2τ−8.6ms程度である。なお、τは90°
励起から180°励起までの時間間隔である。この程度
のTEで得られるスピンエコーは十分な信号強度を有す
る。
を図4に示す。このパルスシーケンスは、グラディエン
トエコー(GRE:gradient echo)法の
パルスシーケンスである。
励起用のα°パルスのシーケンスであり、(2)、
(3)、(4)および(5)は、同じくそれぞれ、スラ
イス勾配Gs、リードアウト勾配Gr、フェーズエンコ
ード勾配GpおよびスピンエコーMRのシーケンスであ
る。なお、α°パルスは中心信号で代表する。パルスシ
ーケンスは時間軸tに沿って左から右に進行する。
ンのα°励起が行われる。αは90以下である。このと
きスライス勾配Gsが印加され所定のスライスについて
の選択励起が行われる。
によりスピンのフェーズエンコードが行われる。次に、
リードアウト勾配Grにより先ずスピンをディフェーズ
し、次いでスピンをリフェーズして、グラディエントエ
コーMRを発生させる。グラディエントエコーMRはデ
ータ収集部150によりビューデータとして収集され
る。このようなパスルシーケンスが周期TRで128〜
256回繰り返される。繰り返しのたびにフェーズエン
コード勾配Gpを変更し、毎回異なるフェーズエンコー
ドを行う。これによって、128〜256ビューのビュ
ーデータが得られる。
に関して対称的な波形を持つRF信号となる。中心エコ
ーはα°励起からTE後に生じる。時間TEを適切に選
ぶことにより、水のエコーと脂肪のエコーの位相差を2
π/n(n≧2)とすることができる。例えば、n=4
としたとき、位相差をπ/2にするTEは静磁場強度が
0.2Tの場合で8.6ms程度である。この程度のT
Eで得られるグラディエントエコーは十分な信号強度を
有する。
て得られたビューデータが、データ処理部170のメモ
リに収集される。なお、パルスシーケンスはSE法また
はGRE法に限るものではなく、例えばファーストスピ
ンエコー(FSE:fastspin echo)法
等、他の適宜の技法のものであって良いのはいうまでも
ない。
次元逆フーリエ変換して撮像対象300の断層像を再構
成する。再構成した画像はメモリに記憶する。ここまで
の本装置の構成および機能は、本発明における撮像手段
の実施の形態の一例である。
ら、水分を画像化した像および脂肪分を画像化した像を
それぞれ生成する。以下、水分を画像化した像を水像、
脂肪分を画像化した像を脂肪像という。
ータ処理部170は静磁場の強度分布を求める。データ
処理部170は、本発明の磁場不均一測定装置の実施の
形態の一例である。データ処理部170の構成によっ
て、本発明の装置に関する実施の形態の一例が示され
る。データ処理部170の動作によって、本発明の方法
に関する実施の形態の一例が示される。
に基づいて、磁場不均一の影響を除去する位相補正を行
う。データ処理部170は、本発明の位相補正装置の実
施の形態の一例である。データ処理部170の構成によ
って、本発明の装置に関する実施の形態の一例が示され
る。データ処理部170の動作によって、本発明の方法
に関する実施の形態の一例が示される。
見たデータ処理部170のブロック図を図5に示す。同
図の各ブロックの機能は、例えばコンピュータプログラ
ム(computer program)等により実現
される。
(power)画像形成部702および位相画像形成部
704を有する。パワー画像形成部702および位相画
像形成部704には、前段の画像再構成部700から再
構成画像が入力される。再構成画像のピクセルデータ
(pixel data)は複素数で表されている。す
なわち、ピクセルデータは実数成分と虚数成分を有す
る。以下、実数成分および虚数成分をそれぞれリアルパ
ート(real part)およびイマジナリパート
(imaginary part)という。
の複素数データのパワーを求め、このパワーをピクセル
値とする画像すなわちパワー画像を形成する。位相画像
形成部704は、ピクセルごとに複素数データの位相す
なわちリアルパートとイマジナリパートのアークタンジ
ェント(arc tangent)を求め、この位相を
ピクセル値とする画像すなわち位相画像を形成する。
同図は、断層像が脂肪像とその周囲を囲む水像からなる
場合の、位相画像の1次元プロファイル(profil
e)である。ただし、水像と脂肪像の位相差をπ/2と
した場合である。すなわち、位相差2π/nにおけるn
を4とした場合である。以下、n=4の場合について説
明するが、他の値の場合もこれに準じる。
に位相画像という)は、静磁場が均一であるとすると、
水像の位相が0になることにより、同図の一点鎖線で示
すような図形になるべきであるが、例えば静磁場がリニ
ア(linear)に傾斜す不均一性を持つとすると、
実線で示すような位相画像となる。
る。位相n倍部706は、本発明における位相逓倍手段
の実施の形態の一例である。位相n倍部706は位相画
像について、そのピクセル値をn(=4)倍する。これ
により、図6の(b)に示すような位相画像が得られ
る。同図に示すように、4倍したことにより水と脂肪の
位相差が2πになり両者は同相となる。位相画像には、
±πの範囲を逸脱した分が±πの範囲内に折り返すいわ
ゆるラッピング(wrapping)が生じる。それに
加えて、水と脂肪の境界部分では位相の不連続ないし急
変が生じる。
8に入力される。複素画像形成部708にはパワー画像
形成部702からパワー画像も入力される。複素画像形
成部708は、位相画像とパワー画像に基づいて複素画
像を形成する。
ータのコサイン(cosine)として求められる。複
素画像のイマジナリパートは、パワー画像データのサイ
ン(sine)として求められる。なお、コサインおよ
びサインの演算に用いる角度は位相角度である。
ass filter)部710を通して位相画像形成
部712に入力する。位相画像形成部712は、ローパ
スフィルタリングされた複素画像から位相画像を形成す
る。ローパスフィルタリングにより、位相画像は、図7
の(a)に示すような位相の不連続ないし急変部分が、
例えば(b)に示すように連続化ないし急変緩和された
ものとなる。
(unwrapping)部714に入力される。アン
ラッピング部714は、本発明における折り返し補正手
段の実施の形態の一例である。アンラッピング部714
は、図8の(a)に示すようにラッピングしている位相
を(b)のようにアンラッピングする。
n倍部716に入力される。位相1/n倍部716は、
本発明における磁場不均一測定手段の実施の形態の一例
である。位相1/n倍部716は入力位相画像につい
て、そのピクセル値を1/n(=4)倍する。これによ
り、図8の(c)に示すような位相画像が得られる。こ
の位相画像は、撮像対象300が水だけからなる場合の
位相画像に相当する。したがって、この位相画像は静磁
場の強度分布すなわち静磁場不均一を表すものとなる。
718に入力される。位相補正部718は、本発明にお
ける位相補正手段の実施の形態の一例である。位相補正
部718には、位相画像形成部704からも位相画像が
入力される。この位相画像は、静磁場不均一による位相
誤差を含んだ水と脂肪の位相画像である。位相補正部7
18は、水と脂肪の位相画像における位相誤差を、静磁
場不均一を表す位相画像で補正する。
含んだ水と脂肪の位相画像(a)から静磁場不均一を表
す位相画像(b)を対応するピクセル同士で減算し、
(c)に示すように静磁場不均一による位相誤差を含ま
ない位相画像を得る。
0に入力される。複素画像形成部720は、また、パワ
ー画像形成部702からパワー画像が入力される。複素
画像形成部720は、位相画像とパワー画像に基づいて
複素画像を形成する。複素画像のリアルパートおよびイ
マジナリパートは、パワー画像データのコサインおよび
サインとしてそれぞれ求められる。なお、コサインおよ
びサインの演算に用いる角度は位相角度である。
る。位相修正部722は、本発明における位相修正手段
の実施の形態の一例である。位相修正部722は、ヒス
トグラム(histogram)生成部724から入力
されるヒストグラムに基づいて複素画像の位相を修正す
る。
位相のヒストグラムを生成する。これによって、例えば
図10に示すようなヒストグラムが得られる。同図に示
すように、ヒストグラムは2つのピーク(peak)を
持つ。位相0の位置にピークを持つものは水であり、位
相π/2付近にピークを持つものは脂肪である。
内に水も含むので、位相のヒストグラムのピークはπ/
2からずれる。その理由は、例えば図11に示すよう
に、脂肪信号のベクトル(vector)Fは、脂肪ベ
クトルfと水ベクトルwの和となるからである。位相修
正部722は、脂肪ヒストグラムのピークのずれδに応
じて脂肪像の位相を修正する。これによって、水像と脂
肪像の位相差をπ/2とする。
726に入力される。水・脂肪分離部726は、本発明
における画像生成手段の実施の形態の一例である。水・
脂肪分離部726は、複素画像のリアルパートを用いて
水像を生成し、イマジナリパートを用いて脂肪像を生成
する。生成した水像は水像メモリ728に記憶し、脂肪
像は脂肪像メモリ730に記憶する。
ら静磁場の不均一を表す情報を得て、静磁場不均一によ
る位相誤差を補正し、位相誤差補正後の画像から水像と
脂肪像を分離するので、水・脂肪分離像の撮像をきわめ
て能率良く行うことができる。
れば、水と脂肪を分離した撮像を能率良く行うための磁
場不均一測定方法および装置、磁場不均一に基づく位相
補正方法および装置、並びに、水と脂肪を分離した撮像
を能率良く行う磁気共鳴撮像装置を実現することができ
る。
である。
である。
シーケンスの一例を示す図である。
シーケンスの一例を示す図である。
理部のブロック図である。
である。
である。
である。
である。
ある。
Claims (6)
- 【請求項1】 静磁場空間中の対象物体について磁気共
鳴を利用して水のピクセルデータと脂肪のピクセルデー
タが2π/n(n≧2)の位相差を持つ断層像を撮像す
る撮像手段と、 前記断層像のピクセルデータの位相をn倍する位相逓倍
手段と、 前記n倍した位相が±πの範囲を超えている位相画像か
ら複素画像を形成する複素画像形成手段と、 前記形成された複素画像に対してローパスフィルタリン
グを施すローパスフィルタリング手段と、 前記ローパスフィルタリングが施された複素画像から位
相画像を形成する位相画像形成手段と、 前記形成された位相画像においてラッピングしている位
相をアンラッピングすることにより、位相が±πの範囲
を超えることによる折り返しを補正する位相アンラッピ
ング手段と、 前記折り返しを補正した位相を1/nして前記静磁場空
間の磁場不均一を表す位相を求める磁場不均一測定手段
とを具備することを特徴とする磁場不均一測定装置。 - 【請求項2】 静磁場空間中の対象物体について磁気共
鳴を利用して水のピクセルデータと脂肪のピクセルデー
タが2π/n(n≧2)の位相差を持つ断層像を撮像す
る撮像手段と、 前記断層像のピクセルデータの位相をn倍する位相逓倍
手段と、 前記n倍した位相が±πの範囲を超えている位相画像か
ら複素画像を形成する複素画像形成手段と、 前記形成された複素画像に対してローパスフィルタリン
グを施すローパスフィルタリング手段と、 前記ローパスフィルタリングが施された複素画像から位
相画像を形成する位相画像形成手段と、 前記形成された位相画像においてラッピングしている位
相をアンラッピングすることにより、位相が±πの範囲
を超えることによる折り返しを補正する位相アンラッピ
ング手段と、 前記折り返しを補正した位相を1/nして前記静磁場空
間の磁場不均一を表す位相を求める磁場不均一測定手段
と、 前記求めた位相を用いて前記断層像のピクセルデータに
つき位相補正を行う位相補正手段とを具備することを特
徴とする位相補正装置。 - 【請求項3】 静磁場空間中の対象物体について磁気共
鳴を利用して水のピクセルデータと脂肪のピクセルデー
タが2π/n(n≧2)の位相差を持つ断層像を撮像す
る撮像手段と、 前記断層像のピクセルデータの位相をn倍する位相逓倍
手段と、 前記n倍した位相が±πの範囲を超えている位相画像か
ら複素画像を形成する複素画像形成手段と、 前記形成された複素画像に対してローパスフィルタリン
グを施すローパスフィルタリング手段と、 前記ローパスフィルタリングが施された複素画像から位
相画像を形成する位相画像形成手段と、 前記形成された位相画像においてラッピングしている位
相をアンラッピングすることにより、位相が±πの範囲
を超えることによる折り返しを補正する位相アンラッピ
ング手段と、 前記折り返しを補正した位相を1/nして前記静磁場空
間の磁場不均一を表す位相を求める磁場不均一測定手段
と、 前記求めた位相を用いて前記断層像のピクセルデータに
つき位相補正を行う位相補正手段と、 前記位相補正したピクセルデータの位相差を利用して水
画像と脂肪画像を別々に生成する画像生成手段とを具備
することを特徴とする磁気共鳴撮像装置。 - 【請求項4】 n=4 であることを特徴とする請求項
3に記載の磁気共鳴撮像装置。 - 【請求項5】 前記位相補正したピクセルデータの位相
差を前記ピクセルデータの位相のヒストグラムに基づい
て修正する位相修正手段を有し、 前記画像形成手段は、前記修正後の位相差を利用して水
画像と脂肪画像を別々に生成することを特徴とする請求
項3又は請求項4に記載の磁気共鳴撮像装置。 - 【請求項6】 静磁場空間中の対象物体について磁気共
鳴を利用して水のピクセルデータと脂肪のピクセルデー
タが2π/n(n≧2)の位相差を持つ断層像を撮像す
る撮像手段と、 前記断層像のピクセルデータの位相をn倍する位相逓倍
手段と、 前記n倍した位相が±πの範囲を超えることによる折り
返しを補正する折り返し補正手段と、 前記折り返しを補正した位相を1/nして前記静磁場空
間の磁場不均一を表す位相を求める磁場不均一測定手段
と、 前記求めた位相を用いて前記断層像のピクセルデータに
つき位相補正を行う位相補正手段と、 前記位相補正したピクセルデータの位相差を利用して水
画像と脂肪画像を別々に生成する画像生成手段と、 前記位相補正したピクセルデータの位相差を前記ピクセ
ルデータの位相のヒストグラムに基づいて修正する位相
修正手段とを具備し、 前記画像形成手段は、前記修正後の位相差を利用して水
画像と脂肪画像を別々に生成することを特徴とする磁気
共鳴撮像装置。
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP17765899A JP3353826B2 (ja) | 1999-06-24 | 1999-06-24 | 磁場不均一測定装置、位相補正装置および磁気共鳴撮像装置 |
US09/575,595 US6515476B1 (en) | 1999-06-24 | 2000-05-22 | Magnetic field inhomogeneity measurement method and apparatus, phase correction method and apparatus, and magnetic resonance imaging apparatus |
EP00305092A EP1063532A3 (en) | 1999-06-24 | 2000-06-16 | Magnetic field inhomogeneity measurement method and apparatus, Phase correction method and apparatus |
KR10-2000-0034734A KR100405923B1 (ko) | 1999-06-24 | 2000-06-23 | 자장 불균일 측정 방법 및 장치, 위상 보정 방법 및 장치, 및 자기 공명 촬상 방법 및 장치 |
CNB001222651A CN1278647C (zh) | 1999-06-24 | 2000-06-24 | 磁场非均匀性测量方法和设备,相位校正方法和设备,以及磁共振成像设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP17765899A JP3353826B2 (ja) | 1999-06-24 | 1999-06-24 | 磁場不均一測定装置、位相補正装置および磁気共鳴撮像装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2001000414A JP2001000414A (ja) | 2001-01-09 |
JP3353826B2 true JP3353826B2 (ja) | 2002-12-03 |
Family
ID=16034849
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP17765899A Expired - Fee Related JP3353826B2 (ja) | 1999-06-24 | 1999-06-24 | 磁場不均一測定装置、位相補正装置および磁気共鳴撮像装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US6515476B1 (ja) |
EP (1) | EP1063532A3 (ja) |
JP (1) | JP3353826B2 (ja) |
KR (1) | KR100405923B1 (ja) |
CN (1) | CN1278647C (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7199583B2 (en) | 2004-09-10 | 2007-04-03 | Ge Medical Systems Global Technology Company, Llc | MR image production method and MRI apparatus |
Families Citing this family (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10122874B4 (de) * | 2001-05-11 | 2004-09-23 | Siemens Ag | Verfahren zur Extraktion von Spinkollektiven mit unterschiedlicher chemischer Verschiebung aus phasenkodierten Einzelbildern unter Berücksichtigung von Feldinhomogenitäten sowie Vorrichtung dazu |
JP3796426B2 (ja) * | 2001-10-04 | 2006-07-12 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 磁気共鳴撮像装置 |
US6703835B2 (en) * | 2002-04-11 | 2004-03-09 | Ge Medical Systems Global Technology Co. Llc | System and method for unwrapping phase difference images |
EP1623246A1 (en) * | 2003-04-14 | 2006-02-08 | Koninklijke Philips Electronics N.V. | Magnetic resonance imaging with histogram-based phase correction |
JP2005152114A (ja) * | 2003-11-21 | 2005-06-16 | Ge Medical Systems Global Technology Co Llc | Mri方法およびmri装置 |
JP3968352B2 (ja) * | 2004-02-03 | 2007-08-29 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Mri装置 |
JP3968353B2 (ja) * | 2004-02-18 | 2007-08-29 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Mri装置 |
JP4685496B2 (ja) * | 2005-04-11 | 2011-05-18 | 株式会社日立メディコ | 磁気共鳴イメージング装置 |
DE102006050551A1 (de) * | 2006-10-26 | 2008-06-26 | Robert Bosch Gmbh | Vorrichtung zum Bestimmen einer Phasenlage von zwei Magnetfeldern und Verfahren zum Überprüfen von Statorwicklungen |
US7557572B1 (en) * | 2006-12-19 | 2009-07-07 | Fonar Corporation | Method of detecting and compensating for faults in an MRI process |
CN101612042B (zh) * | 2008-06-26 | 2012-03-07 | 株式会社东芝 | 磁共振成像装置及磁共振成像方法 |
DE102009015006B4 (de) * | 2009-03-26 | 2011-09-22 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zum automatischen Unterscheiden von silikondominiertem, wasserdominiertem und fettdominiertem Gewebe |
RU2521751C2 (ru) * | 2009-03-30 | 2014-07-10 | Конинклейке Филипс Электроникс Н.В. | Двухточечный метод диксона с гибким выбором времени эха |
AU2010327289B2 (en) * | 2009-12-02 | 2015-05-28 | Nanalysis Corp. | Method and apparatus for producing homogeneous magnetic fields |
CN102621509B (zh) * | 2011-01-26 | 2014-11-26 | 西门子(深圳)磁共振有限公司 | 一种匀场装置的校正方法及装置 |
CN102727201B (zh) * | 2011-04-13 | 2014-02-12 | 深圳迈瑞生物医疗电子股份有限公司 | 磁共振系统及其水脂分离成像方法、装置 |
US9625549B2 (en) * | 2011-04-19 | 2017-04-18 | Diego Hernando | Method for error compensated chemical species signal separation with magnetic resonance imaging |
CN103371818B (zh) | 2012-04-27 | 2015-05-13 | 西门子(深圳)磁共振有限公司 | 水脂图像辨析方法及装置、一致化方法及装置 |
JP6013161B2 (ja) * | 2012-09-06 | 2016-10-25 | 株式会社日立製作所 | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 |
CN104918546B (zh) * | 2013-02-12 | 2018-01-02 | 株式会社日立制作所 | 磁共振成像装置及其处理方法 |
DE102013217651B4 (de) * | 2013-09-04 | 2016-03-24 | Siemens Aktiengesellschaft | Mehrpunkt Dixon-Technik |
JP5908878B2 (ja) | 2013-09-18 | 2016-04-26 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 磁気共鳴装置およびプログラム |
CN104155621B (zh) * | 2014-09-11 | 2017-01-18 | 中国科学院武汉物理与数学研究所 | 一种准确测量静磁场b0分布的方法 |
EP3314284B1 (en) * | 2015-06-26 | 2021-10-20 | Koninklijke Philips N.V. | Phase corrected dixon magnetic resonance imaging |
US9935646B2 (en) | 2015-07-17 | 2018-04-03 | Schweitzer Engineering Laboratories, Inc. | Systems and methods for identifying a failure in an analog to digital converter |
US9584148B1 (en) * | 2016-01-19 | 2017-02-28 | Schweitzer Engineering Laboratories, Inc. | Systems and methods for analog to digital converter failure identification |
DE102016204863B3 (de) * | 2016-03-23 | 2017-06-14 | Siemens Healthcare Gmbh | Verfahren zu einer Bestimmung einer Abweichung einer Homogenität eines Magnetfeldes eines Magnetresonanzgerätes |
CN108013876B (zh) * | 2016-10-31 | 2023-05-23 | 通用电气公司 | 磁共振图像均匀度校正方法及装置 |
EP3364207A1 (en) | 2017-02-17 | 2018-08-22 | Koninklijke Philips N.V. | Phase error correction for bipolar read out gradients |
US10033400B1 (en) | 2017-10-18 | 2018-07-24 | Schweitzer Engineering Laboratories, Inc. | Analog-to-digital converter verification using quantization noise properties |
US9985646B1 (en) | 2017-10-18 | 2018-05-29 | Schweitzer Engineering Laboratories, Inc. | Analog-to-digital converter verification using quantization noise properties |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4661775A (en) * | 1985-07-15 | 1987-04-28 | Technicare Corporation | Chemical shift imaging with field inhomogeneity correction |
US4720679A (en) * | 1985-12-31 | 1988-01-19 | Picker International, Inc. | Magnetic resonance imaging with phase encoded chemical shift correction |
IL80814A (en) * | 1986-11-30 | 1990-07-12 | Elscint Ltd | Spectral component separation in mr imaging |
US4797615A (en) * | 1987-09-30 | 1989-01-10 | Elscint Ltd. | Determining and correcting for phase jumps |
US5602480A (en) * | 1994-04-19 | 1997-02-11 | Hitachi Medical Corporation | Inspection method and apparatus using nuclear magnetic resonance |
US5498963A (en) * | 1994-08-05 | 1996-03-12 | General Electric Company | Correction of geometric distortion in MRI phase and phase difference images |
GB2304419A (en) * | 1995-08-18 | 1997-03-19 | Toshiba America Mri Inc | Method and apparatus for providing separate fat and water MRI images in a single acquisition scan |
US5617029A (en) * | 1995-12-14 | 1997-04-01 | General Electric Company | Method of magnet shimming |
US5652514A (en) * | 1996-03-25 | 1997-07-29 | Toshiba America Mri, Inc. | Correction for field variation in steady-state MRI by repeated acquisition of zero k-space line |
US5701074A (en) * | 1996-04-25 | 1997-12-23 | Eiscint Ltd. | Spectral component separation including unwrapping of the phase via a poisson equation utilizing a weighting map |
US6091243A (en) * | 1997-11-13 | 2000-07-18 | The University Of British Columbia | Water-fat imaging with direct phase encoding (DPE) |
US6177795B1 (en) * | 1998-05-19 | 2001-01-23 | Elscint Ltd. | Spectral component imaging using phased array coils |
US6263228B1 (en) * | 1998-08-27 | 2001-07-17 | Toshiba America, Mri, Inc. | Method and apparatus for providing separate water-dominant and fat-dominant images from single scan single point dixon MRI sequences |
US6147492A (en) * | 1998-10-28 | 2000-11-14 | Toshiba America Mri, Inc. | Quantitative MR imaging of water and fat using a quadruple-echo sequence |
US6259253B1 (en) * | 1999-06-10 | 2001-07-10 | Ohio State University Research Foundation | MRI transceiver |
-
1999
- 1999-06-24 JP JP17765899A patent/JP3353826B2/ja not_active Expired - Fee Related
-
2000
- 2000-05-22 US US09/575,595 patent/US6515476B1/en not_active Expired - Fee Related
- 2000-06-16 EP EP00305092A patent/EP1063532A3/en not_active Withdrawn
- 2000-06-23 KR KR10-2000-0034734A patent/KR100405923B1/ko not_active IP Right Cessation
- 2000-06-24 CN CNB001222651A patent/CN1278647C/zh not_active Expired - Fee Related
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7199583B2 (en) | 2004-09-10 | 2007-04-03 | Ge Medical Systems Global Technology Company, Llc | MR image production method and MRI apparatus |
Also Published As
Publication number | Publication date |
---|---|
EP1063532A3 (en) | 2009-03-18 |
EP1063532A2 (en) | 2000-12-27 |
JP2001000414A (ja) | 2001-01-09 |
CN1279052A (zh) | 2001-01-10 |
CN1278647C (zh) | 2006-10-11 |
US6515476B1 (en) | 2003-02-04 |
KR20010007498A (ko) | 2001-01-26 |
KR100405923B1 (ko) | 2003-11-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3353826B2 (ja) | 磁場不均一測定装置、位相補正装置および磁気共鳴撮像装置 | |
Benkert et al. | Free‐breathing volumetric fat/water separation by combining radial sampling, compressed sensing, and parallel imaging | |
JP3454760B2 (ja) | 位相分布測定方法および装置、位相補正方法および装置、並びに、磁気共鳴撮像装置 | |
JP4049649B2 (ja) | 磁気共鳴撮影装置 | |
US8148981B2 (en) | MRI apparatus and MRI method for SSFP with center frequency and 1st order gradient moments zeroed | |
EP0383631A2 (en) | Magnetic resonance imaging system | |
US6906515B2 (en) | Magnetic resonance imaging device and method | |
US20070285090A1 (en) | Phase cycling method and magnetic resonance imaging apparatus | |
JP3510542B2 (ja) | 位相分布測定方法および装置、位相補正方法および装置、並びに、磁気共鳴撮像装置 | |
JP4192139B2 (ja) | Ssfpパルスシーケンスにおけるフェーズサイクリング法および磁気共鳴撮影装置 | |
JP3884227B2 (ja) | 磁気共鳴撮影装置 | |
JP2001112735A (ja) | 磁気共鳴イメージング装置 | |
JP4901031B2 (ja) | 位相矛盾検出方法および装置、位相矛盾解消方法および装置、並びに、磁気共鳴撮影装置 | |
JP4558866B2 (ja) | 位相分布測定方法および装置、位相補正方法および装置、並びに、磁気共鳴撮像装置 | |
JP3907944B2 (ja) | 磁気共鳴イメージング方法及び装置 | |
JP2002052005A (ja) | 磁気共鳴イメージング方法 | |
US4896111A (en) | Method and system for improving resolution of images in magnetic resonance imaging | |
JPS61106140A (ja) | 高分解能核磁気共鳴イメ−ジング方式 | |
JP3557286B2 (ja) | Mr画像生成方法及びmri装置 | |
JP2022094087A (ja) | 磁気共鳴イメージング装置およびその制御方法 | |
US4996480A (en) | Image reconstructing device for use in magnetic resonance imaging apparatus | |
JP3532307B2 (ja) | 磁気共鳴イメージング装置 | |
JPS61194339A (ja) | 核磁気共鳴撮像装置 | |
JPH07213502A (ja) | 磁気共鳴イメージング方法 | |
JPH0377738B2 (ja) |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080927 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080927 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090927 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090927 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090927 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090927 Year of fee payment: 7 |
|
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 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100927 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100927 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110927 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110927 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120927 Year of fee payment: 10 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120927 Year of fee payment: 10 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130927 Year of fee payment: 11 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |