JPS6321048A - 核磁気共鳴画像の歪み補正方法 - Google Patents
核磁気共鳴画像の歪み補正方法Info
- Publication number
- JPS6321048A JPS6321048A JP61167560A JP16756086A JPS6321048A JP S6321048 A JPS6321048 A JP S6321048A JP 61167560 A JP61167560 A JP 61167560A JP 16756086 A JP16756086 A JP 16756086A JP S6321048 A JPS6321048 A JP S6321048A
- Authority
- JP
- Japan
- Prior art keywords
- distortion
- magnetic resonance
- correction
- nuclear magnetic
- image
- 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
- 238000000034 method Methods 0.000 title claims description 33
- 238000005481 NMR spectroscopy Methods 0.000 title claims description 17
- 238000012937 correction Methods 0.000 claims description 38
- 238000009828 non-uniform distribution Methods 0.000 claims description 17
- 230000003068 static effect Effects 0.000 claims description 16
- 239000011159 matrix material Substances 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 8
- 238000003384 imaging method Methods 0.000 claims description 6
- 230000004069 differentiation Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 description 5
- 238000009434 installation Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 3
- 230000010355 oscillation Effects 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 239000000696 magnetic material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000012546 transfer Methods 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/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
【発明の詳細な説明】
(産業上の利用分野)
本発明は、フーリエ法を用いた核磁気断層撮像装置(以
下、NMR−CTという)の静磁場の不均一性による画
像歪みを補正する方法に関する。
下、NMR−CTという)の静磁場の不均一性による画
像歪みを補正する方法に関する。
(従来の技術)
NMR−CTは、−様な静磁場H8を作る静磁′場コイ
ル及び静磁場H8と同一方向磁場でx、y。
ル及び静磁場H8と同一方向磁場でx、y。
2の各方向に夫々直線勾配をもつ磁場を作る勾配磁場コ
イルから成る磁石部、該磁石部で形成される磁場内に設
置する被検体に高周波パルス(高周波電磁波)を印加し
、被検体からのNMR信号を検出する送・受信部、該送
・受信部及び前記磁石部の動作を制御したり、検出デー
タの!l!IL!T!をして画像表示する制御・画像処
理部等を有している。
イルから成る磁石部、該磁石部で形成される磁場内に設
置する被検体に高周波パルス(高周波電磁波)を印加し
、被検体からのNMR信号を検出する送・受信部、該送
・受信部及び前記磁石部の動作を制御したり、検出デー
タの!l!IL!T!をして画像表示する制御・画像処
理部等を有している。
以上の構成において、静磁場コイルにより所定の空間(
被検体の設置箇所)に所定強度の静磁場が形成される。
被検体の設置箇所)に所定強度の静磁場が形成される。
一方、制御・画像処理部の制御の下で、送・受信部は、
スピン・エコー法(SpinEcho法:SE法)等に
よる励起パルスシーケンス(90°パルス及び180°
パルスを所定のタイミングで発生するシーケンス)及び
フーリエ変換法による勾配磁場シーケンスを出力し、そ
のときのNMR信号を検出する。又、1IIl制御・画
像処理部は、上記シーケンスに従って収集された多数の
NMR信号、即ち、生データ(raw data)を用
いて画像再構成処理をして断層画像の表示を行う。
スピン・エコー法(SpinEcho法:SE法)等に
よる励起パルスシーケンス(90°パルス及び180°
パルスを所定のタイミングで発生するシーケンス)及び
フーリエ変換法による勾配磁場シーケンスを出力し、そ
のときのNMR信号を検出する。又、1IIl制御・画
像処理部は、上記シーケンスに従って収集された多数の
NMR信号、即ち、生データ(raw data)を用
いて画像再構成処理をして断層画像の表示を行う。
ところで、静1itu場が不均一であると、その大きさ
に比例した位置歪み、濃度歪み等の画像歪みが生じるこ
とがよく知られている。一方、均一性の高い静磁場を形
成し得る磁石部は、一般的に、高価になる傾向にあるの
で、コスト面から理想的な磁石部を得ることが難しい。
に比例した位置歪み、濃度歪み等の画像歪みが生じるこ
とがよく知られている。一方、均一性の高い静磁場を形
成し得る磁石部は、一般的に、高価になる傾向にあるの
で、コスト面から理想的な磁石部を得ることが難しい。
しかも、磁石部の性能を確保するために、磁石部は、周
囲に磁性体が少ない場所を選んで設置される必要がある
(設置場所の制約がある)。又、画像のS/Nを上げる
には、できるだけ低い読みだし勾配を用いる方がよいこ
とが分っているが、画像の歪みの方は、読みだし勾配を
低くするほど大きくなるため、本来エコー時間や空間分
解能等で決まっている限界まで読みだし勾配を低くでき
ないのが現状である。
囲に磁性体が少ない場所を選んで設置される必要がある
(設置場所の制約がある)。又、画像のS/Nを上げる
には、できるだけ低い読みだし勾配を用いる方がよいこ
とが分っているが、画像の歪みの方は、読みだし勾配を
低くするほど大きくなるため、本来エコー時間や空間分
解能等で決まっている限界まで読みだし勾配を低くでき
ないのが現状である。
これらのことから、画像の歪み補正を行うことにより磁
石部のコスト低減、磁石部の設置場所の制限緩和、低勾
配スキャンによるS/Nの改善等様々な効果が期待でき
る。
石部のコスト低減、磁石部の設置場所の制限緩和、低勾
配スキャンによるS/Nの改善等様々な効果が期待でき
る。
従来から、核磁気共鳴画像の歪みを補正する方法として
、スキャンをする面と同一の面で、予め不均一を測って
おいて、それにより2次元面内での補正を行う方法があ
る( 0278−0062/85/1200−0193
$01.OO19851EEE 193〜199頁)。
、スキャンをする面と同一の面で、予め不均一を測って
おいて、それにより2次元面内での補正を行う方法があ
る( 0278−0062/85/1200−0193
$01.OO19851EEE 193〜199頁)。
又、3次元での補正を目的としたもので、3次元の不均
一分布を数式でもっていて、それを用いて任意面又は立
体の不均一を計算し、補正を行う方法もある。
一分布を数式でもっていて、それを用いて任意面又は立
体の不均一を計算し、補正を行う方法もある。
(発明が解決しようとする問題点)
しかし、上記の前者の方法にあっては、電気的に任意の
断面のイメージングが可能なNMR−CTにおいて実用
上重要な任意面への拡張が考虞されていないため、実用
的ではないという問題がある。又、後者の方法にあって
は、元来、シムコイルで補正できない高次項を数式とし
てもつことになるが、そのような数式は計1fflも多
くなるため、スピードが要求されるオンラインの再構成
には向かないという問題がある。更に、立体での不均一
分布を直接もち、かつ、イメージングのマトリックスと
同程度(1283〜2563)のサイズをもつ場合、デ
ータ世が膨大なものとなるので実用上の問題もある。
断面のイメージングが可能なNMR−CTにおいて実用
上重要な任意面への拡張が考虞されていないため、実用
的ではないという問題がある。又、後者の方法にあって
は、元来、シムコイルで補正できない高次項を数式とし
てもつことになるが、そのような数式は計1fflも多
くなるため、スピードが要求されるオンラインの再構成
には向かないという問題がある。更に、立体での不均一
分布を直接もち、かつ、イメージングのマトリックスと
同程度(1283〜2563)のサイズをもつ場合、デ
ータ世が膨大なものとなるので実用上の問題もある。
本発明は、かかる点に鑑みてなされたものであり、その
目的は、実用性の高い核磁気共鳴画像の歪み補正方法を
提供することにある。
目的は、実用性の高い核磁気共鳴画像の歪み補正方法を
提供することにある。
(問題点を解決するための手段)
上記目的を達成する本発明の核磁気共鳴画像の歪み補正
方法は、核磁気共鳴画像の歪み補正をするとぎ、通常の
イメージングに用いられるマトリックスよりも少ないマ
トリックスで立体の不均一分布データを用意し、実際の
補正用データは、前記不均一分布データから所定の補間
を用いて作成するようになっている。
方法は、核磁気共鳴画像の歪み補正をするとぎ、通常の
イメージングに用いられるマトリックスよりも少ないマ
トリックスで立体の不均一分布データを用意し、実際の
補正用データは、前記不均一分布データから所定の補間
を用いて作成するようになっている。
(実施例)
以下、本発明について詳細に説明する。
第1図は、本発明の補正方法が実施されるNMR−CT
の概念構成回である。マグネットアセンブリ1は、内部
に被検体を挿入するための空間部分(孔)を有し、この
空間部分を取巻くようにして、被検体に一定の静磁場を
印加する静磁場コイルと、勾配磁場を発生するための勾
配li場ココイル勾配磁場コイルは、x、y、zの各軸
のコイルを備えている)と、被検体内の原子核のスピン
を励起するための高周波パルスを与えるRF送信コイル
と、被検体からのNMR信号を検出する受信コイル等が
配置されている。静磁場コイル、勾配磁場コイル、RF
送信コイル及び受信コイルは、静磁揚重′IQ(主磁場
電源)2、勾配磁場駆動回路3、RF電力増幅器4及び
前置増幅器5夫々に接続されている。シーケンス記憶回
路6は、泪Ω磯7からの指令に従って任意のビューで、
ゲート変調回路8を操作(所定のタイミングによってR
F発振回路9の高周波出力信号を変調)し、SR法に基
づく高周波パルス信号をRFF力増幅34からRF送1
8コイルに印加する。又、シーケンス記憶回路6は、フ
ーリエ変換法に基づくシーケンス信号によって勾配磁場
駆動回路3、ゲート変調回路8及びA/D変換器10を
操作する。位相検波器11は、RF発発註回路9出力を
参照信号とし、前置増幅器5の出力信号(受信コイルで
検出されたNMR信号)を位相検波してA/D変換器1
0に与える。A/D変換器10は、位相検波器11を介
して得られるNMR信号をアナログ・ディジタル変換し
て計算機7に入力する。計算機7は、操作コンソール1
2との間で情報の授受や秤々のスキャンシーケンスを実
現するために、シーケンス記憶回路6の動作の切替えや
メモリの書替えを′したり、A/D変換器10データを
用いてスキャンシーケンスに応じた所定の処理をして、
画像化した後、ディスプレイ13に画像表示をする。
の概念構成回である。マグネットアセンブリ1は、内部
に被検体を挿入するための空間部分(孔)を有し、この
空間部分を取巻くようにして、被検体に一定の静磁場を
印加する静磁場コイルと、勾配磁場を発生するための勾
配li場ココイル勾配磁場コイルは、x、y、zの各軸
のコイルを備えている)と、被検体内の原子核のスピン
を励起するための高周波パルスを与えるRF送信コイル
と、被検体からのNMR信号を検出する受信コイル等が
配置されている。静磁場コイル、勾配磁場コイル、RF
送信コイル及び受信コイルは、静磁揚重′IQ(主磁場
電源)2、勾配磁場駆動回路3、RF電力増幅器4及び
前置増幅器5夫々に接続されている。シーケンス記憶回
路6は、泪Ω磯7からの指令に従って任意のビューで、
ゲート変調回路8を操作(所定のタイミングによってR
F発振回路9の高周波出力信号を変調)し、SR法に基
づく高周波パルス信号をRFF力増幅34からRF送1
8コイルに印加する。又、シーケンス記憶回路6は、フ
ーリエ変換法に基づくシーケンス信号によって勾配磁場
駆動回路3、ゲート変調回路8及びA/D変換器10を
操作する。位相検波器11は、RF発発註回路9出力を
参照信号とし、前置増幅器5の出力信号(受信コイルで
検出されたNMR信号)を位相検波してA/D変換器1
0に与える。A/D変換器10は、位相検波器11を介
して得られるNMR信号をアナログ・ディジタル変換し
て計算機7に入力する。計算機7は、操作コンソール1
2との間で情報の授受や秤々のスキャンシーケンスを実
現するために、シーケンス記憶回路6の動作の切替えや
メモリの書替えを′したり、A/D変換器10データを
用いてスキャンシーケンスに応じた所定の処理をして、
画像化した後、ディスプレイ13に画像表示をする。
次に、上記構成の動作について説明する。
計算様7における補正演算は、補正に必要な静磁場の不
均一分布データを得るためのキャリブレイションの部分
と、そのデータに基づいて実際の補正を行う補正部分と
に大別できる。キVリプレイジョンは、オフラインの作
業であり、頻度も低いのでこれにかかる時間や処理の複
雑さは余り問題ではない。これに対し、補正部分は、オ
ンラインの処理であり、できる限り簡単で、演算量の少
ないことが要求される。
均一分布データを得るためのキャリブレイションの部分
と、そのデータに基づいて実際の補正を行う補正部分と
に大別できる。キVリプレイジョンは、オフラインの作
業であり、頻度も低いのでこれにかかる時間や処理の複
雑さは余り問題ではない。これに対し、補正部分は、オ
ンラインの処理であり、できる限り簡単で、演算量の少
ないことが要求される。
第2図に不均一測定シーケンスの例を示す。ここに示し
た方法は、90°〜180°の時[1εと180゜〜エ
コー中心の時間とを等しくでとしてスキャンしたデータ
と、その間隔をτ−ε/2、τ+ε/2としてスキャン
したデータとから得られる2つの複素画像間の位相差φ
(x、 y、 z)から、(1)式に基づいて不均一分
布D(x、 y、 z)を求めるものである。
た方法は、90°〜180°の時[1εと180゜〜エ
コー中心の時間とを等しくでとしてスキャンしたデータ
と、その間隔をτ−ε/2、τ+ε/2としてスキャン
したデータとから得られる2つの複素画像間の位相差φ
(x、 y、 z)から、(1)式に基づいて不均一分
布D(x、 y、 z)を求めるものである。
φ(x、y、z)
D(x、y、z)= −(1)
2πγε
但し、γ・・・磁気回転比
この方法は、リードアウト方向へ不均一分布に応じた位
置歪みが生ずるので、それが十分小さくなるよう強い勾
配を用いるか、後で何かの位置補正を行う必要がある。
置歪みが生ずるので、それが十分小さくなるよう強い勾
配を用いるか、後で何かの位置補正を行う必要がある。
この方法の他に3方向ともワーブを行い、読みだし勾配
をかけずに信号を読みだして(従ってデータは4次元と
なる)、各ボクセルのだすNMR信号の周波数から不均
一分布を求める方法もある。これは、時間がかかるが位
置歪みはない。これらは、いずれも不均一分布を測りた
い領域全体をカバーするファントムを用いて測定を行う
。尚、本発明は、不均一測定方法を限定する必要はない
ので、この他の測定方法によって不均一分布を求めるよ
うにしてもよい。
をかけずに信号を読みだして(従ってデータは4次元と
なる)、各ボクセルのだすNMR信号の周波数から不均
一分布を求める方法もある。これは、時間がかかるが位
置歪みはない。これらは、いずれも不均一分布を測りた
い領域全体をカバーするファントムを用いて測定を行う
。尚、本発明は、不均一測定方法を限定する必要はない
ので、この他の測定方法によって不均一分布を求めるよ
うにしてもよい。
こうして特定された3次元の不均一分布を保存するとぎ
に、そのサンプリング点を余り細かくすると、データm
が膨大なものとなって実用的でなくなるので、これをイ
メージングに通常使用するマトリックスサイズの数分の
1の大きさで保存する。不均一分布データのシンブル間
隔は、5〜10酬程度あれば、特に磁石部設置を意識し
ないで作られた建屋において、通常のシミングで達成で
きる程度の不均一は十分に再現できる。そこで、例えば
、サンプル間隔7#、643マトリツクスとすれば、約
450.の立方体領域をカバーすることができる。64
39262 Kポイントであり、この程度のデータタ隋
は現在のNMR−CTシステムの情報処理系においては
十分に扱うことができる聞であると言える。
に、そのサンプリング点を余り細かくすると、データm
が膨大なものとなって実用的でなくなるので、これをイ
メージングに通常使用するマトリックスサイズの数分の
1の大きさで保存する。不均一分布データのシンブル間
隔は、5〜10酬程度あれば、特に磁石部設置を意識し
ないで作られた建屋において、通常のシミングで達成で
きる程度の不均一は十分に再現できる。そこで、例えば
、サンプル間隔7#、643マトリツクスとすれば、約
450.の立方体領域をカバーすることができる。64
39262 Kポイントであり、この程度のデータタ隋
は現在のNMR−CTシステムの情報処理系においては
十分に扱うことができる聞であると言える。
第3図は、通常の2次元スキセンの場合の歪み補正のフ
ローチャートである。まず、立体の不均一分布データか
ら、補正したいイメージのスライス面と同一面内の不均
一分布を所定の補間演算により切出す。次に、読みだし
勾配の方向とその強さの情報に基づいて不均一分布から
位置歪みテーブル(各々のビクセルが不均一によりどの
ような位置ヘシフトするかを格納したテーブル)を作る
。
ローチャートである。まず、立体の不均一分布データか
ら、補正したいイメージのスライス面と同一面内の不均
一分布を所定の補間演算により切出す。次に、読みだし
勾配の方向とその強さの情報に基づいて不均一分布から
位置歪みテーブル(各々のビクセルが不均一によりどの
ような位置ヘシフトするかを格納したテーブル)を作る
。
不均一による位置歪みは、スピン・ワーブ法においては
、読みだし方向へのみ生じ、その大きさS(u、v)は
、(2)式となる。
、読みだし方向へのみ生じ、その大きさS(u、v)は
、(2)式となる。
S(u、v)= D(u、v)/ gpr
(2)但し、D(u、v)・・・不均一の大きざQ
、r・・・読みだし勾配の大きさ いま、仮にU軸が読みだし軸とした場合、点(u、v)
は、(u+s(u、 v)、 v)という位置シフトす
ることになる。次に、位置歪みの情報から濃度歪みの補
正テーブルを作る。一般に、(u、v)→(u−、v−
)という位置歪みが生ずると、それに応じて(3)式に
示すJ(u、v) (Jacobian)に反比例した
濃度歪みが生ずる。
(2)但し、D(u、v)・・・不均一の大きざQ
、r・・・読みだし勾配の大きさ いま、仮にU軸が読みだし軸とした場合、点(u、v)
は、(u+s(u、 v)、 v)という位置シフトす
ることになる。次に、位置歪みの情報から濃度歪みの補
正テーブルを作る。一般に、(u、v)→(u−、v−
)という位置歪みが生ずると、それに応じて(3)式に
示すJ(u、v) (Jacobian)に反比例した
濃度歪みが生ずる。
au−Elv−9u′ 9v−
J(u、v)= −−−□
au eV av 9U・・・(3)
そこで、その位置歪みを補正するには、J(u、v)を
求めて乗ずればよい。そのとき、位置歪みは、数式で与
えられていないので数値微分を用いる。
求めて乗ずればよい。そのとき、位置歪みは、数式で与
えられていないので数値微分を用いる。
以上により得られた位置歪みと濃度歪みの補正データは
、一般にマトリックスサイズ及びピクセルサイズ共に被
補正イメージのそれとは異なるので、所定のりサンプル
(補間)をして、サイズを揃えた後、補正に使用される
。
、一般にマトリックスサイズ及びピクセルサイズ共に被
補正イメージのそれとは異なるので、所定のりサンプル
(補間)をして、サイズを揃えた後、補正に使用される
。
第4図は、位置歪みの補正を示したものである。
補正後の(正しい座標系での)座標(U、V)に位置す
るビクセルは、補正前は(歪んだ座標系では)座標(u
′、v″)に存在することが位置歪みテーブルから分る
ので、補正前のイメージで(U、V)の値をとってきて
補正済みのイメージブレーンの(U、V)に格納すれば
よい。一般に、(u′、v−)は、ビクセル位置とは正
確に一致しないので、周囲のビクセルから所定の補間に
より求められる。このようにして、位置歪みを補正した
後、各ビクセルにa度歪み補正係数J(u、v)をかけ
れば、位置歪み及び濃度歪みが補正された画像が得られ
る。位置補正のとぎには、補正前後のイメージのFOV
が等しいと、とってくるべきデータの位置がイメージの
外へ出てしまうことがある。このため、補正前のイメー
ジのFOVは補正後のそれに比べて多少大きく設定して
おくことが望ましい。又、以下に説明するロータリー式
の補正を単独、又は、前記の「OVの操作と組合せて用
いることができる。第5図は、その方法を示したもので
ある。ここではイメージの範囲外の部分には、同一のイ
メージが繰返し並んでるとして、(U“、V“)からと
ってきたデータを(11,V)に植えている。、これは
、フーリエ法画像が本来ロータリー的性質をもつことを
利用したものである。
るビクセルは、補正前は(歪んだ座標系では)座標(u
′、v″)に存在することが位置歪みテーブルから分る
ので、補正前のイメージで(U、V)の値をとってきて
補正済みのイメージブレーンの(U、V)に格納すれば
よい。一般に、(u′、v−)は、ビクセル位置とは正
確に一致しないので、周囲のビクセルから所定の補間に
より求められる。このようにして、位置歪みを補正した
後、各ビクセルにa度歪み補正係数J(u、v)をかけ
れば、位置歪み及び濃度歪みが補正された画像が得られ
る。位置補正のとぎには、補正前後のイメージのFOV
が等しいと、とってくるべきデータの位置がイメージの
外へ出てしまうことがある。このため、補正前のイメー
ジのFOVは補正後のそれに比べて多少大きく設定して
おくことが望ましい。又、以下に説明するロータリー式
の補正を単独、又は、前記の「OVの操作と組合せて用
いることができる。第5図は、その方法を示したもので
ある。ここではイメージの範囲外の部分には、同一のイ
メージが繰返し並んでるとして、(U“、V“)からと
ってきたデータを(11,V)に植えている。、これは
、フーリエ法画像が本来ロータリー的性質をもつことを
利用したものである。
尚、本発明は、上記実施例にのみ限定するものではなく
、3次元法による立体のイメージングの歪み補正にも適
用できる。その場合、第3図において、任意面内の不均
一データの切出しを省き、その後の処理を全て3次元で
行えばよい。但し、そのときは、位置歪み、濃度補正係
数の拡大後のテーブルは、非常に大きくなる場合がある
ので、テーブルとはしないで、補正計算中に値を求めな
がら補正を行う等の工夫が必要となる。
、3次元法による立体のイメージングの歪み補正にも適
用できる。その場合、第3図において、任意面内の不均
一データの切出しを省き、その後の処理を全て3次元で
行えばよい。但し、そのときは、位置歪み、濃度補正係
数の拡大後のテーブルは、非常に大きくなる場合がある
ので、テーブルとはしないで、補正計算中に値を求めな
がら補正を行う等の工夫が必要となる。
又、本発明は、勾配磁場の空間非直線性により生ずる画
像歪みの補正と、そのがなりの部分を共用することがで
きる。勾配による歪みは、不均一による歪みと加算の関
係になるため、第3図の位置歪みテーブルに勾配によっ
て生ずる歪みを加えれば、それ以後の処理は全て共用す
ることができる。
像歪みの補正と、そのがなりの部分を共用することがで
きる。勾配による歪みは、不均一による歪みと加算の関
係になるため、第3図の位置歪みテーブルに勾配によっ
て生ずる歪みを加えれば、それ以後の処理は全て共用す
ることができる。
(発明の効!J!)
以上、説明の通り、本発明の核磁気共鳴画像の歪み補正
方法によれば、実現可能なメモリ容量及び演算量で、任
意の切断面、又は、立体のイメージに対して、静磁場の
不均一性から生じる位置歪み及び濃度歪みを効果的に除
去することが可能となるみその結果、均一度特性の劣る
磁石部を用いても均一度の良い磁石部を用いた場合と同
等の画像が得られ、同−設置条件ならば、磁石部のコス
トを下げることが可能となる。又、同一の磁石部ならば
、より設置条件の悪い場所にら設置することが可能とな
る。更に、読みだし勾配を可能な限り(エコー時間■[
や分解能でで決まる限界まで)下げることにより、S/
Nを上げることができる。
方法によれば、実現可能なメモリ容量及び演算量で、任
意の切断面、又は、立体のイメージに対して、静磁場の
不均一性から生じる位置歪み及び濃度歪みを効果的に除
去することが可能となるみその結果、均一度特性の劣る
磁石部を用いても均一度の良い磁石部を用いた場合と同
等の画像が得られ、同−設置条件ならば、磁石部のコス
トを下げることが可能となる。又、同一の磁石部ならば
、より設置条件の悪い場所にら設置することが可能とな
る。更に、読みだし勾配を可能な限り(エコー時間■[
や分解能でで決まる限界まで)下げることにより、S/
Nを上げることができる。
従って、本発明によって実用性の高い補正方法が実現さ
れる。
れる。
第1図は、NMR−CTの構成図、第2図は、静磁場の
不均一を測定するシーケンスの説明図、第3図は、通常
の2次元スキャンの場合の歪み補正のフローチャート、
第4図は、位置歪みの補正の説明図、第5図は、ロータ
リー式の位置補正の説明図である。 1・・・マグネットアセンブリ、2・・・静磁場電源、
3・・・勾配磁場駆動回路、4・・・RF電力増幅器、
5・・・前置増幅器、6・・・シーケンス記憶回路、7
・・・計01fi、8=・・ゲート変調回路、9・・・
RF発振回路、10・・・A/D変換器、11・・・位
相検波回路、12・・・操作コンソール、13・・・表
示装置。
不均一を測定するシーケンスの説明図、第3図は、通常
の2次元スキャンの場合の歪み補正のフローチャート、
第4図は、位置歪みの補正の説明図、第5図は、ロータ
リー式の位置補正の説明図である。 1・・・マグネットアセンブリ、2・・・静磁場電源、
3・・・勾配磁場駆動回路、4・・・RF電力増幅器、
5・・・前置増幅器、6・・・シーケンス記憶回路、7
・・・計01fi、8=・・ゲート変調回路、9・・・
RF発振回路、10・・・A/D変換器、11・・・位
相検波回路、12・・・操作コンソール、13・・・表
示装置。
Claims (2)
- (1)静磁場内に設置される被検体に予め定めたシーケ
ンスに従って勾配磁場及び高周波電磁波を印加して収集
されるスキャンデータを用いた再構成処理によつて、所
定のマトリックスで表示される核磁気共鳴画像の歪みを
補正する方法において、 前記核磁気共鳴画像の歪み補正をするとき、通常のイメ
ージングに用いられるマトリックスよりも少ないマトリ
ックスで立体の不均一分布データを用意し、実際の補正
用データは、前記不均一分布データから所定の補間を用
いて作成することを特徴とする核磁気共鳴画像の歪み補
正方法。 - (2)濃度歪み補正係数は、前記の少ないマトリックス
サイズのままのデータから数値微分で求めた後、所定の
大きさに拡大して用いることを特徴とする特許請求範囲
第1項の核磁気共鳴画像の歪み補正方法。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP61167560A JPS6321048A (ja) | 1986-07-16 | 1986-07-16 | 核磁気共鳴画像の歪み補正方法 |
EP19870904731 EP0314790A4 (en) | 1986-07-16 | 1987-07-16 | Method of correcting image distortion for nmr imaging apparatus |
US07/320,266 US4982159A (en) | 1986-07-16 | 1987-07-16 | Image distortion correcting method applied to an NMR imaging apparatus |
PCT/JP1987/000518 WO1988000453A1 (en) | 1986-07-16 | 1987-07-16 | Method of correcting image distortion for nmr imaging apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP61167560A JPS6321048A (ja) | 1986-07-16 | 1986-07-16 | 核磁気共鳴画像の歪み補正方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
JPS6321048A true JPS6321048A (ja) | 1988-01-28 |
Family
ID=15851988
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP61167560A Pending JPS6321048A (ja) | 1986-07-16 | 1986-07-16 | 核磁気共鳴画像の歪み補正方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US4982159A (ja) |
EP (1) | EP0314790A4 (ja) |
JP (1) | JPS6321048A (ja) |
WO (1) | WO1988000453A1 (ja) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5484045A (en) * | 1994-03-16 | 1996-01-16 | Usui Kokusai Sangyo Kaisha Ltd. | Fluid clutch |
US6125981A (en) * | 1998-06-17 | 2000-10-03 | Usui Kokusai Sangyo Kaisha Limited | Temperature sensitive fluid type fan coupling apparatus |
US6814033B2 (en) | 2001-07-03 | 2004-11-09 | Usui Kokusai Sangyo Kaisha Limited | Temperature sensitive fluid type fan coupling apparatus |
JP2007159718A (ja) * | 2005-12-12 | 2007-06-28 | Toshiba Corp | 磁気共鳴イメージング装置および画像補正評価方法 |
JP2009160409A (ja) * | 2007-12-29 | 2009-07-23 | Siemens Ag | 磁気共鳴イメージングにおける歪み補正方法及び装置 |
WO2009107492A1 (ja) * | 2008-02-29 | 2009-09-03 | 株式会社東芝 | 磁気共鳴イメージング装置及び磁気共鳴イメージング方法 |
JP2009273501A (ja) * | 2008-05-12 | 2009-11-26 | Hitachi Medical Corp | 磁気共鳴イメージング装置、画像再構成方法およびプログラム |
US8991581B2 (en) | 2010-04-28 | 2015-03-31 | Usui Kokusai Sangyo Kaisha Limitd | Highly reactive fluid fan coupling device |
EP2947345A1 (en) | 2014-05-19 | 2015-11-25 | Usui Kokusai Sangyo Kaisha Limited | Fluid fan clutch device |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2741885B2 (ja) * | 1989-02-17 | 1998-04-22 | 株式会社日立製作所 | 磁気共鳴を用いた検査装置におけるデータ処理方法 |
US5356613A (en) * | 1991-03-18 | 1994-10-18 | Mallinckrodt Medical, Inc. | X-ray contrast agents, compositions and methods |
JP3701616B2 (ja) * | 2002-03-06 | 2005-10-05 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 磁気共鳴撮影装置 |
US7911208B2 (en) * | 2007-10-15 | 2011-03-22 | Siemens Aktiengesellschaft | Methods for rectification of B0 inhomogeneity effects in magnetic resonance images |
US7535227B1 (en) * | 2007-10-26 | 2009-05-19 | General Electric Company | Method and apparatus for correcting distortion in MR images caused by metallic implants |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS59148854A (ja) * | 1983-02-14 | 1984-08-25 | Hitachi Ltd | 核磁気共鳴を用いた検査装置 |
JPS60209154A (ja) * | 1983-12-23 | 1985-10-21 | ゼネラル・エレクトリツク・カンパニイ | Nmr像の歪みを補正する方法 |
JPS61151450A (ja) * | 1984-12-25 | 1986-07-10 | Yokogawa Electric Corp | Nmr断層撮像装置の画像感度補正方法 |
JPS6319138A (ja) * | 1986-07-14 | 1988-01-26 | 株式会社日立製作所 | 核磁気共鳴を用いた検査装置 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4284950A (en) * | 1978-08-05 | 1981-08-18 | E M I Limited | Imaging systems |
JPS59190643A (ja) * | 1983-04-14 | 1984-10-29 | Hitachi Ltd | 核磁気共鳴を用いた検査装置 |
JPS60161552A (ja) * | 1984-02-01 | 1985-08-23 | Hitachi Ltd | Νmr検査装置の静磁場強度分布測定方法 |
US4585992A (en) * | 1984-02-03 | 1986-04-29 | Philips Medical Systems, Inc. | NMR imaging methods |
JPS61280551A (ja) * | 1985-05-16 | 1986-12-11 | Sanyo Electric Co Ltd | 磁気共鳴イメ−ジング装置 |
US4724388A (en) * | 1985-06-07 | 1988-02-09 | Hitachi, Ltd. | Magnetic resonance imaging method |
JPH0636790B2 (ja) * | 1985-06-07 | 1994-05-18 | 株式会社日立製作所 | 核磁気共鳴イメ−ジング静磁場分布測定方式 |
-
1986
- 1986-07-16 JP JP61167560A patent/JPS6321048A/ja active Pending
-
1987
- 1987-07-16 US US07/320,266 patent/US4982159A/en not_active Expired - Fee Related
- 1987-07-16 EP EP19870904731 patent/EP0314790A4/en not_active Withdrawn
- 1987-07-16 WO PCT/JP1987/000518 patent/WO1988000453A1/ja not_active Application Discontinuation
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS59148854A (ja) * | 1983-02-14 | 1984-08-25 | Hitachi Ltd | 核磁気共鳴を用いた検査装置 |
JPS60209154A (ja) * | 1983-12-23 | 1985-10-21 | ゼネラル・エレクトリツク・カンパニイ | Nmr像の歪みを補正する方法 |
JPS61151450A (ja) * | 1984-12-25 | 1986-07-10 | Yokogawa Electric Corp | Nmr断層撮像装置の画像感度補正方法 |
JPS6319138A (ja) * | 1986-07-14 | 1988-01-26 | 株式会社日立製作所 | 核磁気共鳴を用いた検査装置 |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5484045A (en) * | 1994-03-16 | 1996-01-16 | Usui Kokusai Sangyo Kaisha Ltd. | Fluid clutch |
US6125981A (en) * | 1998-06-17 | 2000-10-03 | Usui Kokusai Sangyo Kaisha Limited | Temperature sensitive fluid type fan coupling apparatus |
US6814033B2 (en) | 2001-07-03 | 2004-11-09 | Usui Kokusai Sangyo Kaisha Limited | Temperature sensitive fluid type fan coupling apparatus |
JP2007159718A (ja) * | 2005-12-12 | 2007-06-28 | Toshiba Corp | 磁気共鳴イメージング装置および画像補正評価方法 |
JP2009160409A (ja) * | 2007-12-29 | 2009-07-23 | Siemens Ag | 磁気共鳴イメージングにおける歪み補正方法及び装置 |
JP2009226199A (ja) * | 2008-02-29 | 2009-10-08 | Toshiba Corp | 磁気共鳴イメージング装置 |
WO2009107492A1 (ja) * | 2008-02-29 | 2009-09-03 | 株式会社東芝 | 磁気共鳴イメージング装置及び磁気共鳴イメージング方法 |
US8368397B2 (en) | 2008-02-29 | 2013-02-05 | Kabushiki Kaisha Toshiba | Magnetic resonance imaging apparatus and magnetic resonance imaging method |
JP2009273501A (ja) * | 2008-05-12 | 2009-11-26 | Hitachi Medical Corp | 磁気共鳴イメージング装置、画像再構成方法およびプログラム |
US8991581B2 (en) | 2010-04-28 | 2015-03-31 | Usui Kokusai Sangyo Kaisha Limitd | Highly reactive fluid fan coupling device |
EP2947345A1 (en) | 2014-05-19 | 2015-11-25 | Usui Kokusai Sangyo Kaisha Limited | Fluid fan clutch device |
KR20150133126A (ko) | 2014-05-19 | 2015-11-27 | 우수이 고쿠사이 산교 가부시키가이샤 | 고반응형 유체식 팬·클러치 장치 |
US9664239B2 (en) | 2014-05-19 | 2017-05-30 | Usui Kokusai Sangyo Kaisha Limited | Highly reactive fluid fan clutch device |
Also Published As
Publication number | Publication date |
---|---|
EP0314790A1 (en) | 1989-05-10 |
EP0314790A4 (en) | 1990-12-12 |
US4982159A (en) | 1991-01-01 |
WO1988000453A1 (en) | 1988-01-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6160397A (en) | Fast spin echo prescan for magnetic resonance imaging systems | |
US5498963A (en) | Correction of geometric distortion in MRI phase and phase difference images | |
JP4030805B2 (ja) | コイル感度マップ作成方法およびmri装置 | |
JPH0213434A (ja) | 磁気共鳴イメージング方法およびシステム | |
JPH0616765B2 (ja) | Nmr装置 | |
JPS6321048A (ja) | 核磁気共鳴画像の歪み補正方法 | |
JPH0448451B2 (ja) | ||
US5617029A (en) | Method of magnet shimming | |
US4713614A (en) | Method of correcting the phase and shading in a nuclear magnetic resonance tomographic device | |
CN111190131B (zh) | 用于磁共振图像重建的方法和磁共振系统 | |
JPH0747023B2 (ja) | 核磁気共鳴を用いた検査装置 | |
JP3440049B2 (ja) | Mri装置 | |
JPS62254744A (ja) | 物体の1部分の核磁化分布を測定する方法及び装置 | |
US6636756B2 (en) | K-space synthesis for MR imaging in the presence of gradient field nonlinearity | |
US5184075A (en) | Method and apparatus for compensating for nonuniformity of static magnetic field in MRI system | |
JPS62189056A (ja) | 磁界の均質性を改善する方法 | |
JP3538249B2 (ja) | Mris装置 | |
JPH03292934A (ja) | 核磁気共鳴を用いた検査方法 | |
JP4912802B2 (ja) | 磁気共鳴イメージング装置,送信感度分布計測装置および送信感度分布計測方法 | |
JPH0556140B2 (ja) | ||
JP3112929B2 (ja) | 磁気共鳴装置 | |
JPH0576518A (ja) | 磁気共鳴イメージング装置 | |
JP3112930B2 (ja) | 磁気共鳴装置 | |
JP4901627B2 (ja) | 磁気共鳴撮影装置 | |
EP4343356A1 (en) | Mr imaging with water/fat/b0 mapping |