JPS63135147A - Ｎｍｒ撮像方法 - Google Patents
Ｎｍｒ撮像方法Info
- Publication number
- JPS63135147A JPS63135147A JP61283687A JP28368786A JPS63135147A JP S63135147 A JPS63135147 A JP S63135147A JP 61283687 A JP61283687 A JP 61283687A JP 28368786 A JP28368786 A JP 28368786A JP S63135147 A JPS63135147 A JP S63135147A
- Authority
- JP
- Japan
- Prior art keywords
- magnetic field
- pulse
- image
- gradient
- chemical shift
- 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.)
- Pending
Links
- 238000003384 imaging method Methods 0.000 title claims description 17
- 239000000126 substance Substances 0.000 claims description 28
- 238000000034 method Methods 0.000 claims description 24
- 230000003068 static effect Effects 0.000 claims description 16
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 10
- 238000012937 correction Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 10
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 10
- 238000005481 NMR spectroscopy Methods 0.000 description 9
- 238000005259 measurement Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 239000002131 composite material Substances 0.000 description 2
- 230000005415 magnetization Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 241000600169 Maro Species 0.000 description 1
- 125000004429 atom Chemical group 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 239000012847 fine chemical Substances 0.000 description 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000002075 inversion recovery Methods 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000000425 proton nuclear magnetic resonance spectrum Methods 0.000 description 1
- 238000009738 saturating Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000009987 spinning Methods 0.000 description 1
- 125000004434 sulfur atom Chemical group 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
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/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/56527—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to chemical shift effects
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
（産業上の利用分野）
本発明はフーリエ変換法により断層像を得るＮＭＲ−Ｃ
Ｔのプロトンイメージングに関し、特に読み出し勾配を
低勾配にして装置のＳＮ比を改善するＮＭＲＩＩＩ像方
法に関する。
Ｔのプロトンイメージングに関し、特に読み出し勾配を
低勾配にして装置のＳＮ比を改善するＮＭＲＩＩＩ像方
法に関する。
（従来の技術）
核磁気共鳴（以下ＮＭＲという）現象を用いて特定原子
核に注目した被検体の断層像を得るＮＭＲ−ＣＴは従来
から知られている。このＮＭＲ−Ｃ１’の原理の概要を
簡単に説明する。
核に注目した被検体の断層像を得るＮＭＲ−ＣＴは従来
から知られている。このＮＭＲ−Ｃ１’の原理の概要を
簡単に説明する。
原子核は磁気を帯びた回転している独楽と見ることがで
きるが、それを例えばＺ軸方向の静＆ｆｉＪａＨｏの中
におくと、前記の原子核は次式で示す角速度ω０で歳差
運動をする。これをラーモアの歳差運動という。
きるが、それを例えばＺ軸方向の静＆ｆｉＪａＨｏの中
におくと、前記の原子核は次式で示す角速度ω０で歳差
運動をする。これをラーモアの歳差運動という。
ω０−γＨｏ 但し、γ：核磁気回転比今、静磁場
のあるｌ軸に垂直な軸、例えばＸ軸に高周波コイルを配
置し、ｘｙ面内で回転する前記の角周波数ω０の高周波
回転磁場を印加すると磁気共鳴が起り、静磁場１１ ｏ
のもとでゼーマン分裂をしていた原子核の集団は共鳴条
件を満足する高周波ｌｌ場によって単位間の遷移を生じ
、エネルギ一単位の高い方の単位に遷移する。ここで、
核磁気回転比γは原子核の種類によって異なるので共鳴
周波数によって当該原子核を特定することができる。更
にその共鳴の強さを測定すれば、その原子核の存在門を
知ることができる。共鳴後緩和時閣と呼ばれる時定数で
定まる時間の間に高い準位へ励起された原子核は低い準
位へ戻ってエネルギーの放射を行う。
のあるｌ軸に垂直な軸、例えばＸ軸に高周波コイルを配
置し、ｘｙ面内で回転する前記の角周波数ω０の高周波
回転磁場を印加すると磁気共鳴が起り、静磁場１１ ｏ
のもとでゼーマン分裂をしていた原子核の集団は共鳴条
件を満足する高周波ｌｌ場によって単位間の遷移を生じ
、エネルギ一単位の高い方の単位に遷移する。ここで、
核磁気回転比γは原子核の種類によって異なるので共鳴
周波数によって当該原子核を特定することができる。更
にその共鳴の強さを測定すれば、その原子核の存在門を
知ることができる。共鳴後緩和時閣と呼ばれる時定数で
定まる時間の間に高い準位へ励起された原子核は低い準
位へ戻ってエネルギーの放射を行う。
このＮＭＲの現象の観測方法の中スピンエコー（ＳＥ）
法について第７図を参照しながら説明する。
法について第７図を参照しながら説明する。
前述のように共鳴条件を満足する高周波パルス（）（１
）を静磁場（２軸）に垂直な（Ｘ軸）方向に印加すると
、第７図（イ）に示すように磁化ベクトルＭは回転座標
系でω′−γＨ１の角周波数でｚｙ面内で回転を始める
。今パルス幅をｔＤとするとＨＯからの回転角θは次式
で表わされる。
）を静磁場（２軸）に垂直な（Ｘ軸）方向に印加すると
、第７図（イ）に示すように磁化ベクトルＭは回転座標
系でω′−γＨ１の角周波数でｚｙ面内で回転を始める
。今パルス幅をｔＤとするとＨＯからの回転角θは次式
で表わされる。
θ＝γＨ＋ｊｏ ・・・（１）（
１）式においてθ＝９０°となるようなｔｏをもつパル
スを９０°パルスと呼ぶ。この９０°パルス直後では磁
化ベクトルＭは第７図（ロ）のようにｘｙ面をω０で回
転していることになり、例えばＸ軸においたコイルに誘
導起・電力を生じる。しかし、この信号は時間と共に減
衰していくので、この信号を自由誘導減衰信号（ＦｒＤ
信号）と呼ぶ。ＦＩＤ信号をフーリエ変換すれば周波数
領域での信号が得られる。次に第７図（ハ）に示すよう
に９０°パルスからτ時間後θ−１８０°になるような
パルス幅の第２のパルス（１８０°パルス）を加えると
ばらばらになっていた磁気モーメントがτ時間後−ｙ方
向で再び焦点を合せて信号が観測される。この信号をＳ
Ｅ倍信号呼んでいる。
１）式においてθ＝９０°となるようなｔｏをもつパル
スを９０°パルスと呼ぶ。この９０°パルス直後では磁
化ベクトルＭは第７図（ロ）のようにｘｙ面をω０で回
転していることになり、例えばＸ軸においたコイルに誘
導起・電力を生じる。しかし、この信号は時間と共に減
衰していくので、この信号を自由誘導減衰信号（ＦｒＤ
信号）と呼ぶ。ＦＩＤ信号をフーリエ変換すれば周波数
領域での信号が得られる。次に第７図（ハ）に示すよう
に９０°パルスからτ時間後θ−１８０°になるような
パルス幅の第２のパルス（１８０°パルス）を加えると
ばらばらになっていた磁気モーメントがτ時間後−ｙ方
向で再び焦点を合せて信号が観測される。この信号をＳ
Ｅ倍信号呼んでいる。
このＳＥ倍信号強度を測定して所望の像を得ることがで
きる。ＮＭＲの共鳴条件は ν８γＨＯ／２π ・・・（２）で
与えられる。ここで、νは共鳴周波数、Ｈａは静磁場の
強さである。従って共鳴周波数は磁場の強さに比例する
ことが分る。このため静磁場に線形の磁場勾配を重畳さ
せて、位置によって異なる強さの磁場を与え、共鳴周波
数を変化させて位置情報を得るＮＭＲイメージングの方
法がある。この内フーリエ変換法について説明する。こ
の手法に用いる高周波磁場及び勾配磁場印加のパルスシ
ーケンスを第８図に示す。（イ）図において、Ｘ。
きる。ＮＭＲの共鳴条件は ν８γＨＯ／２π ・・・（２）で
与えられる。ここで、νは共鳴周波数、Ｈａは静磁場の
強さである。従って共鳴周波数は磁場の強さに比例する
ことが分る。このため静磁場に線形の磁場勾配を重畳さ
せて、位置によって異なる強さの磁場を与え、共鳴周波
数を変化させて位置情報を得るＮＭＲイメージングの方
法がある。この内フーリエ変換法について説明する。こ
の手法に用いる高周波磁場及び勾配磁場印加のパルスシ
ーケンスを第８図に示す。（イ）図において、Ｘ。
ｙ、ｚ軸に夫々Ｇｘ 、ＧＶ、ＧＺの勾配磁場を与え、
高周波磁場をＸ軸に印加する状態を示している。（ロ）
図は夫々の磁場を印加するタイミングを示す図である。
高周波磁場をＸ軸に印加する状態を示している。（ロ）
図は夫々の磁場を印加するタイミングを示す図である。
図においてＲＦは高周波の回転磁場で９０″パルスと１
８０°パルスをＸ軸に印加する。ＧＸはＸ軸に印加する
固定の勾配磁場、Ｇｙはｙ軸に印加する時間によって振
幅を変化させる勾配磁場、Ｇｚはｚ軸に印加する固定の
勾配磁場である。信号は１８０”パルス後のＳＥ倍信号
示している。ＦＩＤ信号はＱｚを加えたためほとんど見
えない。期間は各軸に与える勾配磁場の信号の時期を示
すために設けである。期間１において９０°パルスと勾
配磁場Ｑｚ÷によって２−０を中心とするｚ軸に垂直な
断層撮影におけるスライス面内のスピンが選択的に励起
される。期間２のＱｘ＋はスピンの位相を乱れさせて１
８０゜パルスで反転させるためのもので、ディフェーズ
勾配と呼ばれる。Ｇｚ−はＱｚ◆によって乱れたスピン
の位相を元に戻すためのものである。期間２ではＧｙｎ
も印加する。これはｙ方向の位置に比例してスピンの位
相をずらしてやるためのもので、その強度は毎周期異な
るように制御される。期間３において１８０°パルスを
与えて再び磁気モーメントを揃え、その侵に現われるＳ
Ｅ倍信号観察する。期間４のＱｘ÷は乱れた位相を揃え
、５Ｅｆｉ号を生じさせるための勾配磁場で読み出し勾
配といい、読み出し勾配とディフェーズ勾配の面積が等
しくなったところにＳＥ倍信号現われる。
８０°パルスをＸ軸に印加する。ＧＸはＸ軸に印加する
固定の勾配磁場、Ｇｙはｙ軸に印加する時間によって振
幅を変化させる勾配磁場、Ｇｚはｚ軸に印加する固定の
勾配磁場である。信号は１８０”パルス後のＳＥ倍信号
示している。ＦＩＤ信号はＱｚを加えたためほとんど見
えない。期間は各軸に与える勾配磁場の信号の時期を示
すために設けである。期間１において９０°パルスと勾
配磁場Ｑｚ÷によって２−０を中心とするｚ軸に垂直な
断層撮影におけるスライス面内のスピンが選択的に励起
される。期間２のＱｘ＋はスピンの位相を乱れさせて１
８０゜パルスで反転させるためのもので、ディフェーズ
勾配と呼ばれる。Ｇｚ−はＱｚ◆によって乱れたスピン
の位相を元に戻すためのものである。期間２ではＧｙｎ
も印加する。これはｙ方向の位置に比例してスピンの位
相をずらしてやるためのもので、その強度は毎周期異な
るように制御される。期間３において１８０°パルスを
与えて再び磁気モーメントを揃え、その侵に現われるＳ
Ｅ倍信号観察する。期間４のＱｘ÷は乱れた位相を揃え
、５Ｅｆｉ号を生じさせるための勾配磁場で読み出し勾
配といい、読み出し勾配とディフェーズ勾配の面積が等
しくなったところにＳＥ倍信号現われる。
このフーリエ変換法による撮像の従来の方法におけるＲ
Ｆ軸とＧＸ軸とのパルスシーケンス及びＳＥ倍信号時間
関係を第６図に示す。図において、１はＲＦ軸に印加し
た高周波の回転磁場の９０゜パルスでＦＩＤ信号を誘起
させる。（スライス勾配を掛けた場合はとんど見えない
。）２はＸ軸に印加したスピンの位相を乱れさせるディ
フェーズ勾配で、ディフェーズ勾配後ＲＦ軸に１８０°
パルス３を印加してスピンを反転させ、読み出し勾配４
により位相を整えてＳＥ信号５を誘起させる。
Ｆ軸とＧＸ軸とのパルスシーケンス及びＳＥ倍信号時間
関係を第６図に示す。図において、１はＲＦ軸に印加し
た高周波の回転磁場の９０゜パルスでＦＩＤ信号を誘起
させる。（スライス勾配を掛けた場合はとんど見えない
。）２はＸ軸に印加したスピンの位相を乱れさせるディ
フェーズ勾配で、ディフェーズ勾配後ＲＦ軸に１８０°
パルス３を印加してスピンを反転させ、読み出し勾配４
により位相を整えてＳＥ信号５を誘起させる。
ディフェーズ勾配のパルス振幅をＧＯ，パルス幅をＴｏ
、ｎみ出し勾配のパルス振幅をＧＲＩパルス幅をＴ大
とずれば、次式を満足するＴＲの位置にＳＥ信＠５が誘
起される。
、ｎみ出し勾配のパルス振幅をＧＲＩパルス幅をＴ大
とずれば、次式を満足するＴＲの位置にＳＥ信＠５が誘
起される。
ＧＤ −ＴＤ −ＧＲ−ＴＩＩＩ ・・
・ （３）又、９０’パルスと１８０’パルスの間隔を
ＴＥ／２とし、この２倍の時間間１１１ＴＥにおける点
をＴＥ点６とする。
・ （３）又、９０’パルスと１８０’パルスの間隔を
ＴＥ／２とし、この２倍の時間間１１１ＴＥにおける点
をＴＥ点６とする。
このようなパルスシーケンスにより撮像するのであるが
、例えば最も重要な水素原子核のプロトンＮ Ｍ Ｒに
ついて考察する。
、例えば最も重要な水素原子核のプロトンＮ Ｍ Ｒに
ついて考察する。
一般に物質中ではＮＭＲの共鳴周波数νは外から加えた
静磁場Ｈｏによる共鳴条件式である（２）式からずれる
場合が多い。この原因の一つは電子の軌道運動によって
、Ｈｏと逆向きに生じる１１１（反磁場）によるもので
、電子による核の遮蔽効果と呼ばれている。従って共鳴
周波数は、この遮蔽効果による定数φを用いて、次式に
より与えられる。
静磁場Ｈｏによる共鳴条件式である（２）式からずれる
場合が多い。この原因の一つは電子の軌道運動によって
、Ｈｏと逆向きに生じる１１１（反磁場）によるもので
、電子による核の遮蔽効果と呼ばれている。従って共鳴
周波数は、この遮蔽効果による定数φを用いて、次式に
より与えられる。
シー（γＰＪ ／ ２π）Ｈｏ（１−φ）この様な共鳴
周波数のずれをケミカルシフトと呼び、低磁場側に共鳴
が観測される。そして、このケミカルシフトは化学結合
によって、電子密度の分布や電子の軌道運動が変化し、
φが異なってくる。そのため同じ原子でも化合物によっ
て共鳴条件が異なっている。
周波数のずれをケミカルシフトと呼び、低磁場側に共鳴
が観測される。そして、このケミカルシフトは化学結合
によって、電子密度の分布や電子の軌道運動が変化し、
φが異なってくる。そのため同じ原子でも化合物によっ
て共鳴条件が異なっている。
（Ｒ明が解決しようとする問題点）
ところでこのフーリエ変換法でｍ像する場合、ＴＥが長
いとき、読み出し勾配４の振幅ＧＲを小さく即ち低勾配
にするとノイズ帯域が狭くなるためＳＮ比が改善される
が、周波数分解能が細かくなった分だけ前記のケミカル
シフトによるケミカルシフトアーティファクト（水と油
がずれて映る）が現われたり、静磁場の不均一に起因す
る歪が出たりして成る程度までしか読み出し勾配４の振
幅ＧＲを小さく（低読出勾配に）することはできなかっ
た。
いとき、読み出し勾配４の振幅ＧＲを小さく即ち低勾配
にするとノイズ帯域が狭くなるためＳＮ比が改善される
が、周波数分解能が細かくなった分だけ前記のケミカル
シフトによるケミカルシフトアーティファクト（水と油
がずれて映る）が現われたり、静磁場の不均一に起因す
る歪が出たりして成る程度までしか読み出し勾配４の振
幅ＧＲを小さく（低読出勾配に）することはできなかっ
た。
本発明は、上記の点に鑑みてなされたもので、その目的
は、ＳＮ比改善のため読み出し勾配をできるだ１ノ低勾
配にしても、ケミカルシフトアーティファクトや静磁場
不均一による歪の無い画像を得ることのできるＮＭＲ撮
像方法を得ることにある。
は、ＳＮ比改善のため読み出し勾配をできるだ１ノ低勾
配にしても、ケミカルシフトアーティファクトや静磁場
不均一による歪の無い画像を得ることのできるＮＭＲ撮
像方法を得ることにある。
（問題点を解決するための手段）
前記の問題を解決する本発明は、フーリエ変換法により
断層像を得るＮＭＲ−ＣＴのプロトンイメージングにお
いて、読み出し勾配を可能な限りの低勾配としたケミカ
ルシフトイメージングを行い水と脂肪の分離像を得、そ
の結果に基づきケミカルシフトに対応した位置ずれ補正
を実施し、必要に応じて静磁場不均一による歪補正処理
を行うことを特徴と１６ものである。
断層像を得るＮＭＲ−ＣＴのプロトンイメージングにお
いて、読み出し勾配を可能な限りの低勾配としたケミカ
ルシフトイメージングを行い水と脂肪の分離像を得、そ
の結果に基づきケミカルシフトに対応した位置ずれ補正
を実施し、必要に応じて静磁場不均一による歪補正処理
を行うことを特徴と１６ものである。
（作用）
必要に応じて磁場不均一度を測定し、次いで読み出し勾
配を可能な限り低勾配にしてケミカルシフトイメージン
グを行い、その結果に基づいて周波数及び位相ずれを補
正して画像表示をする。
配を可能な限り低勾配にしてケミカルシフトイメージン
グを行い、その結果に基づいて周波数及び位相ずれを補
正して画像表示をする。
（実施例）
以下、図面を参照して本発明の実施例を詳細に説明する
。実施例のＮＭＲＩＩ像方法を説明する前に、磁場不均
一度とケミカルシフト量の測定法について説明する。
。実施例のＮＭＲＩＩ像方法を説明する前に、磁場不均
一度とケミカルシフト量の測定法について説明する。
・磁場不均−麿測定は次の通り行う。先ず、撮像エリア
を十分カバーする水ファントムを第４図に示すシーケン
スによって撮像する。図において、第６図と同じ部分に
は同じ符号を付しである。
を十分カバーする水ファントムを第４図に示すシーケン
スによって撮像する。図において、第６図と同じ部分に
は同じ符号を付しである。
（イ）図は通常のＳＥ信号を得るシーケンスで、ディフ
ェーズ勾配２の振幅Ｇｏと、読み出し勾配４の振幅ＧＲ
とは等しく取っである。（ロ）図は９０°パルスと１８
０”パルスとの間隔をτＡ。
ェーズ勾配２の振幅Ｇｏと、読み出し勾配４の振幅ＧＲ
とは等しく取っである。（ロ）図は９０°パルスと１８
０”パルスとの間隔をτＡ。
１８０°パルスとＴＥ点６との間隔をτＢとし、τ二τ
日−τ＾の値を適当に取ってスキャンした場合の図であ
る。この場合も振幅Ｇｏと振幅ＧＦＩは等しく取っであ
る。（イ）図によるイメージをＰ、（ロ）図によるイメ
ージをＱとすると、両者のイメージ間の位相差θは、Ｘ
、Ｖを画像上の各ピクセルの座標として、次式から求め
られる。Ｐ。
日−τ＾の値を適当に取ってスキャンした場合の図であ
る。この場合も振幅Ｇｏと振幅ＧＦＩは等しく取っであ
る。（イ）図によるイメージをＰ、（ロ）図によるイメ
ージをＱとすると、両者のイメージ間の位相差θは、Ｘ
、Ｖを画像上の各ピクセルの座標として、次式から求め
られる。Ｐ。
Ｑは複累データである。
ｅＪＱ（ゞ′ｙ）−Ｑ（ｘ、ｙ）／Ｐ（ｘ、ｙ）・・・
（４） この位相差から各点の磁場不均一度分布ΔＢＯ（Ｘ 、
Ｖ ）は次式により求められる。
（４） この位相差から各点の磁場不均一度分布ΔＢＯ（Ｘ 、
Ｖ ）は次式により求められる。
△５ｏ（Ｘ、Ｖ）−θ（ｘ、ｙ）／２πγτ・・・（５
） この磁気不均一度分布を予め測定しておけば、ビクセル
の位’１１．Ｑ度を元に戻すという手法で磁場不均一に
よる位置歪と濃度歪補正がイメージの領域においてでき
る。
） この磁気不均一度分布を予め測定しておけば、ビクセル
の位’１１．Ｑ度を元に戻すという手法で磁場不均一に
よる位置歪と濃度歪補正がイメージの領域においてでき
る。
次にＣ３ｌ（ケミカルシフトイメージ）の測定は次の通
り行う。第５図はそのパルスシーケンスの図で、符号は
第４図と同じに付しである。Ｃ８■測定ではＴＥを長く
取ってあり、そのため振幅ＧＲを振幅Ｇｏより小さいパ
ルスにしである。ＴＥは長く取った方が良い結果が出る
。この測定においてはスライス面は予め既述の方法によ
って磁場不均一度が分っている面とする。第５図におい
て、ケミカルシフトによって分離しようとする２本のＳ
Ｅ倍信号図示せず）のシフト同σ（Ｈｚ ）と、τ（−
τＢ−τ＾）との間には次式の関係がある。
り行う。第５図はそのパルスシーケンスの図で、符号は
第４図と同じに付しである。Ｃ８■測定ではＴＥを長く
取ってあり、そのため振幅ＧＲを振幅Ｇｏより小さいパ
ルスにしである。ＴＥは長く取った方が良い結果が出る
。この測定においてはスライス面は予め既述の方法によ
って磁場不均一度が分っている面とする。第５図におい
て、ケミカルシフトによって分離しようとする２本のＳ
Ｅ倍信号図示せず）のシフト同σ（Ｈｚ ）と、τ（−
τＢ−τ＾）との間には次式の関係がある。
σ−１／（２τ）
水と脂肪のケミカルシフト（３，５Ｄｌｌ■）を求める
場合、静磁場の強さが決まればτは一義的に決まる。τ
−〇で得たイメージをＩ＾、τ−τで得たイメージをＩ
ｓとし、水のイメージを■Ｒ９脂肪のイメージをＴＦと
すれば、水と脂肪のケミカルシフトはそれぞれ次式の計
算を各ビクセル毎に行って得ることができる。
場合、静磁場の強さが決まればτは一義的に決まる。τ
−〇で得たイメージをＩ＾、τ−τで得たイメージをＩ
ｓとし、水のイメージを■Ｒ９脂肪のイメージをＴＦと
すれば、水と脂肪のケミカルシフトはそれぞれ次式の計
算を各ビクセル毎に行って得ることができる。
ＩＲ＝（Ｉ＾＋Ｉｇ）／２
Ｉｐ＝（Ｉ＾−ｒｅ）／２
但し、■＾、１Ｂについては予め磁場不均一による位相
歪の補正をしておくものとする。
歪の補正をしておくものとする。
第１図は本発明の一実施例のＮＭＲ１１ｉ像方法にお行
方法ルスシーケンスを示す図である。図において、第４
図と同じ部分には同じ符号を付しである。図中、ＴＥは
特に長く取っである。下は読み出し勾配によるデータサ
ンプリング時間で、ＴＥに近い長さになるまで最大限に
取っておく。このシーケンスにおいて既述のＣ８Ｉの測
定において行った方法で、即ちＣ８Ｉに必要なτ−τＢ
−τ＾のアンバランスを生じさける第５図のスキャンと
同様な第１図のスキャンによってＣ８Ｉを得る。
方法ルスシーケンスを示す図である。図において、第４
図と同じ部分には同じ符号を付しである。図中、ＴＥは
特に長く取っである。下は読み出し勾配によるデータサ
ンプリング時間で、ＴＥに近い長さになるまで最大限に
取っておく。このシーケンスにおいて既述のＣ８Ｉの測
定において行った方法で、即ちＣ８Ｉに必要なτ−τＢ
−τ＾のアンバランスを生じさける第５図のスキャンと
同様な第１図のスキャンによってＣ８Ｉを得る。
この場合１尺が長くなるため、ＴＥ点６にエコーを生じ
る（３）式の条件を満足させる読み出し勾配４の振幅Ｇ
尺は必然的に小さくなり、ケミカルシフト即ち位置歪を
生じている。
る（３）式の条件を満足させる読み出し勾配４の振幅Ｇ
尺は必然的に小さくなり、ケミカルシフト即ち位置歪を
生じている。
上記のシーケンスで得られた水と脂肪のイメージはケミ
カルシフトによって３．５ｐｐｍの周波数シフトを受け
ているため、そのずれσだけ周波数をずらせてイメージ
を再合成する。第２図はその状態を示した図で、（イ）
は水のイメージＩＲ１（ロ）は脂肪のイメージＩＦ、（
ハ）はＩＲとＩｐの合成イメージを表わしている。図に
おいて、縦軸のθは位相の歪、横軸のｆは周波数のずれ
を示している。図において明らかなように１尺とＩｐの
間にはσの周波数のずれがあり、これを周波数をずらせ
て再合成し、（ハ）の合成イメージを得る。
カルシフトによって３．５ｐｐｍの周波数シフトを受け
ているため、そのずれσだけ周波数をずらせてイメージ
を再合成する。第２図はその状態を示した図で、（イ）
は水のイメージＩＲ１（ロ）は脂肪のイメージＩＦ、（
ハ）はＩＲとＩｐの合成イメージを表わしている。図に
おいて、縦軸のθは位相の歪、横軸のｆは周波数のずれ
を示している。図において明らかなように１尺とＩｐの
間にはσの周波数のずれがあり、これを周波数をずらせ
て再合成し、（ハ）の合成イメージを得る。
次に磁場不均一による像の歪を除去して正規の位置に正
規のＳ度の画像を得るように歪補正を行う。この状態を
第３図に示す。図において、符号は第２図と同様に付し
である。（イ）はケミカルシフトを補正したイメージ、
（ロ）は磁場不均一による像の歪を補正したイメージの
図である。図において、明らかなように（イ）のイメー
ジはケミカルシフトを補正してケミカルシフトアーティ
ファクトはないが、１１％不均一による像の歪のためイ
メージが変形し、又、濃度も変化を受けているので、既
述の磁場不均一度測定によって得られた（４）の位相差
から求められた（５）式の磁場不均一度分布によって補
正しくハ）の正しい位置、正しい１１１度のイメージを
得ることができる。
規のＳ度の画像を得るように歪補正を行う。この状態を
第３図に示す。図において、符号は第２図と同様に付し
である。（イ）はケミカルシフトを補正したイメージ、
（ロ）は磁場不均一による像の歪を補正したイメージの
図である。図において、明らかなように（イ）のイメー
ジはケミカルシフトを補正してケミカルシフトアーティ
ファクトはないが、１１％不均一による像の歪のためイ
メージが変形し、又、濃度も変化を受けているので、既
述の磁場不均一度測定によって得られた（４）の位相差
から求められた（５）式の磁場不均一度分布によって補
正しくハ）の正しい位置、正しい１１１度のイメージを
得ることができる。
以上のように本実施例の撮像法によれば、低勾配の読み
出し勾配によってもケミカルシフトアーティファクトが
なく、歪のないイメージが得られる。
出し勾配によってもケミカルシフトアーティファクトが
なく、歪のないイメージが得られる。
尚、本発明は上記実施例に限定されるものではない。
（１）シーケンスはＳＥ法に限ることはなく、２次元フ
ーリエ変換法を用いるものならば例えば反転回復法を用
いても良い。
ーリエ変換法を用いるものならば例えば反転回復法を用
いても良い。
（２）ＴＥが短ければケミカルシフト、位置歪自身が目
立たないから本手法は不要であるが、行っても特に害は
ない。
立たないから本手法は不要であるが、行っても特に害は
ない。
（３）磁場不均一度補正は本手法に限るものではなく、
他のアイディアで行っても本発明の範晴に入るものであ
る。
他のアイディアで行っても本発明の範晴に入るものであ
る。
（４）Ｃ８Ｉも本手法に限定しない。ｆｔ＋ｊｒＦの形
でイメージ実部に水、虚部に脂肪を出す方式でも良い。
でイメージ実部に水、虚部に脂肪を出す方式でも良い。
又、ファインケミカルシフト（Ｆｅ２）法や片方のスペ
クトルを飽和させる方式でも良い。
クトルを飽和させる方式でも良い。
（５）多成分のケミカルシフトを補正することも可能、
イのときは上記のＦｅ２法が適している。
イのときは上記のＦｅ２法が適している。
（６）静ｆｆｌ場不均一度が小さい場合には、当然のこ
とながら静磁場不均一度補正は省略しても良い。
とながら静磁場不均一度補正は省略しても良い。
（Ｒ明の効果）
以上詳細に説明したように本発明によれば、良好なＳＮ
比（ＳＮ比はデータサンプリング時間Ｔの平方根に比例
する。）を持ら、然もケミカルシフトアーティクトがな
く、又、歪もないイメージが得られる。
比（ＳＮ比はデータサンプリング時間Ｔの平方根に比例
する。）を持ら、然もケミカルシフトアーティクトがな
く、又、歪もないイメージが得られる。
第１図は本発明の一実施例のパルスシーケンスを示す図
、第２図はケミカルシフトを補正する状態を示す図、第
３図は静磁場不均一による像の歪を補正する状態を示す
図、第４図は静磁場不均一度測定のためのパルスシーケ
ンスを示す図、第５図はケミカルシフト・測定のための
パルスシーケンスを示す図、第６図は従来の撮像方法に
おけるパルスシーケンスを示す図、第７図はＮＭＲ−Ｃ
ＴのＳＥ法の原理の説明図、第８図はＮＭＲ−ＣＴの磁
場のパルスシーケンスを示す図である。 １・・・９０’パルス ２・・・ディフェーズ勾配
３・・・１８０°パルス ４・・・読み出し勾配５・
・・ＳＥ倍信号 ６・・・ＴＥ点特許出願人 横
河メディカルシステム株式会社繭１５く ２、ディフェーズ勾配 ４、読み出し勾配 ５ｉＳＥぽ８ ６ｉＴＥ慮 第２図 （イ） （ロ）
（ハ）角等３ 図 （イ） （ロ） 第４図 ２・ディフェーズ勾配 ４１読み出し勾配 ５；ＳＥ傷信 号５図 第６図 ２；ディフェーズ勾配 ４Ｉ亮み出し勾配 ５ｉＳＥ４Ａ号 ６ｊＴＥ点 第７図 （イ） （Ｏ）（
ハ） 鏑 ８図
、第２図はケミカルシフトを補正する状態を示す図、第
３図は静磁場不均一による像の歪を補正する状態を示す
図、第４図は静磁場不均一度測定のためのパルスシーケ
ンスを示す図、第５図はケミカルシフト・測定のための
パルスシーケンスを示す図、第６図は従来の撮像方法に
おけるパルスシーケンスを示す図、第７図はＮＭＲ−Ｃ
ＴのＳＥ法の原理の説明図、第８図はＮＭＲ−ＣＴの磁
場のパルスシーケンスを示す図である。 １・・・９０’パルス ２・・・ディフェーズ勾配
３・・・１８０°パルス ４・・・読み出し勾配５・
・・ＳＥ倍信号 ６・・・ＴＥ点特許出願人 横
河メディカルシステム株式会社繭１５く ２、ディフェーズ勾配 ４、読み出し勾配 ５ｉＳＥぽ８ ６ｉＴＥ慮 第２図 （イ） （ロ）
（ハ）角等３ 図 （イ） （ロ） 第４図 ２・ディフェーズ勾配 ４１読み出し勾配 ５；ＳＥ傷信 号５図 第６図 ２；ディフェーズ勾配 ４Ｉ亮み出し勾配 ５ｉＳＥ４Ａ号 ６ｊＴＥ点 第７図 （イ） （Ｏ）（
ハ） 鏑 ８図
Claims (1)
- フーリエ変換法により断層像を得るＮＭＲ−ＣＴのプロ
トンイメージングにおいて、読み出し勾配を可能な限り
の低勾配としたケミカルシフトイメージングを行い水と
脂肪の分離像を得、その結果に基づきケミカルシフトに
対応した位置ずれ補正を実施し、必要に応じて静磁場不
均一による歪補正処理を行うことを特徴とするＮＭＲ撮
像方法。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP61283687A JPS63135147A (ja) | 1986-11-28 | 1986-11-28 | Ｎｍｒ撮像方法 |
| EP19870907819 EP0339090A4 (en) | 1986-11-28 | 1987-11-27 | Nmr imaging method |
| PCT/JP1987/000918 WO1988003780A1 (en) | 1986-11-28 | 1987-11-27 | Nmr imaging method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP61283687A JPS63135147A (ja) | 1986-11-28 | 1986-11-28 | Ｎｍｒ撮像方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPS63135147A true JPS63135147A (ja) | 1988-06-07 |
Family
ID=17668773
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP61283687A Pending JPS63135147A (ja) | 1986-11-28 | 1986-11-28 | Ｎｍｒ撮像方法 |
Country Status (3)
| Country | Link |
|---|---|
| EP (1) | EP0339090A4 (ja) |
| JP (1) | JPS63135147A (ja) |
| WO (1) | WO1988003780A1 (ja) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0236846A (ja) * | 1988-07-27 | 1990-02-06 | Toshiba Corp | 磁気共鳴イメージング装置 |
| JP2009532163A (ja) * | 2006-04-06 | 2009-09-10 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 磁気共鳴装置及び方法 |
| JP2015531251A (ja) * | 2012-09-04 | 2015-11-02 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | ディクソン水脂肪分離を用いるプロペラ |
| JP2016514563A (ja) * | 2013-04-03 | 2016-05-23 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | 高ｓｎｒの同相画像及び低ｓｎｒの少なくとも部分的に異相の画像を使用するディクソン型の水／脂肪分離ｍｒｉ |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011063840A1 (en) | 2009-11-26 | 2011-06-03 | Brainlab Ag | Shift compensation in medical image processing |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS60179644A (ja) * | 1983-12-27 | 1985-09-13 | ゼネラル・エレクトリツク・カンパニイ | 空間解像度の高いｎｍｒ像を求める方法と装置 |
| JPS6176944A (ja) * | 1984-09-10 | 1986-04-19 | ゼネラル・エレクトリツク・カンパニイ | Ｎｍｒ像を得る方法 |
| JPS61258152A (ja) * | 1985-05-13 | 1986-11-15 | Hitachi Ltd | 核磁気共鳴を用いた検査装置 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB1601970A (en) * | 1978-05-31 | 1981-11-04 | Nat Res Dev | Methods of deriving image information from objects |
| JPH0785737B2 (ja) * | 1985-04-01 | 1995-09-20 | 株式会社日立製作所 | 核磁気共鳴を用いた検査装置 |
| US5455900A (en) * | 1992-10-20 | 1995-10-03 | Ricoh Company, Ltd. | Image processing apparatus |
| JPH06250650A (ja) * | 1993-02-23 | 1994-09-09 | Kawai Musical Instr Mfg Co Ltd | 電子楽器 |
| JPH06272346A (ja) * | 1993-03-24 | 1994-09-27 | Sumitomo Metal Ind Ltd | デッキプレート及び合成スラブ |
-
1986
- 1986-11-28 JP JP61283687A patent/JPS63135147A/ja active Pending
-
1987
- 1987-11-27 EP EP19870907819 patent/EP0339090A4/en not_active Withdrawn
- 1987-11-27 WO PCT/JP1987/000918 patent/WO1988003780A1/ja not_active Application Discontinuation
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS60179644A (ja) * | 1983-12-27 | 1985-09-13 | ゼネラル・エレクトリツク・カンパニイ | 空間解像度の高いｎｍｒ像を求める方法と装置 |
| JPS6176944A (ja) * | 1984-09-10 | 1986-04-19 | ゼネラル・エレクトリツク・カンパニイ | Ｎｍｒ像を得る方法 |
| JPS61258152A (ja) * | 1985-05-13 | 1986-11-15 | Hitachi Ltd | 核磁気共鳴を用いた検査装置 |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0236846A (ja) * | 1988-07-27 | 1990-02-06 | Toshiba Corp | 磁気共鳴イメージング装置 |
| JP2009532163A (ja) * | 2006-04-06 | 2009-09-10 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 磁気共鳴装置及び方法 |
| JP2015531251A (ja) * | 2012-09-04 | 2015-11-02 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | ディクソン水脂肪分離を用いるプロペラ |
| JP2016514563A (ja) * | 2013-04-03 | 2016-05-23 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | 高ｓｎｒの同相画像及び低ｓｎｒの少なくとも部分的に異相の画像を使用するディクソン型の水／脂肪分離ｍｒｉ |
Also Published As
| Publication number | Publication date |
|---|---|
| WO1988003780A1 (en) | 1988-06-02 |
| EP0339090A4 (en) | 1991-03-20 |
| EP0339090A1 (en) | 1989-11-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR100610990B1 (ko) | Ｍｒ 촬영 방법 및 ｍｒｉ 시스템 | |
| Graumann et al. | Multiple‐spin‐echo imaging with a 2D Fourier method | |
| JPH11267111A (ja) | Ｍｒｉシステムで拡散重み付き画像を成形する方法 | |
| Reeder et al. | Quantification and reduction of ghosting artifacts in interleaved echo‐planar imaging | |
| JP2010502301A (ja) | 交互読み出し勾配を用いて取得される複数ｍｒ画像の整列方法 | |
| US6906515B2 (en) | Magnetic resonance imaging device and method | |
| JPS635734A (ja) | 望ましくない応答信号を抑圧する方法 | |
| JPS63238856A (ja) | 核磁化分布決定方法及び装置 | |
| JPH11322A (ja) | 磁気共鳴映像法（ｍｒｉ）における横方向緩和速度への可逆性寄与分測定方法 | |
| JPH0315455B2 (ja) | ||
| JPH03264046A (ja) | 核磁気共鳴映像法及び装置 | |
| JPS63135147A (ja) | Ｎｍｒ撮像方法 | |
| JPS6134451A (ja) | 核磁気共鳴断層撮影のために試料を励起する方法 | |
| JPH01207044A (ja) | 核磁気共鳴画像診断装置の受信装置 | |
| JPH01170446A (ja) | 核磁気共鳴画像診断装置の領域制限方法 | |
| EP0223279B1 (en) | Method of and device for the phase correction of mr inversion recovery images | |
| JPH0576518A (ja) | 磁気共鳴イメージング装置 | |
| US6046590A (en) | MR imaging method with time dependence provided in multiple spatial encoding directions | |
| Altbach et al. | Chemical-shift imaging utilizing the positional shifts along the readout gradient direction | |
| JPH09154831A (ja) | Ｍｒイメージング方法及びｍｒｉ装置 | |
| JP2602443B2 (ja) | 核磁気共鳴画像撮影装置 | |
| Loecher et al. | k-Space | |
| JP2002301041A (ja) | 磁気共鳴イメージング装置 | |
| JP2824656B2 (ja) | Ｍｒｉ装置 | |
| Gichohi | Wavelet Based Denoising Methods For Magnetic Resonance Images |