JPH01160544A - 磁気共鳴イメージング装置 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 〔産業上の利用分野〕 本発明は、磁気共鳴現象を用いて被検体中に存在する特
定原子核のスピン密度、及び緩和時定数等の反映された
画像を得る磁気共鳴イメージング装置に関するものであ
る。

〔従来の技術〕

磁気共鳴は、ある原子核（水素原子核など）を一様な静
磁場中においたとき、原子核のスピンが静磁場の強度に
比例した周波数で磁場のまわシを歳差運動し、この状態
で物体に歳差運動の周波数と同一の周波数をもつ高周波
パルス（ＲＦバルスンを印加すると、磁気共鳴現象を生
ずる性質を利用したものである。このとき、特定の方向
に沿って強さが変化する傾斜磁場を前記静磁場に重畳さ
せると、物体の中の各原子核は、磁場強度の変化に応じ
て異なった周波数で共鳴及び緩和の運動を行うことにな
多、その結果、周波数の分析によ）物体の各原子核のマ
ツピング（画像）を得ることができる。測定の方法とし
ては一般に、前記高周波パルスのうち９０＠や１８０＠
パルスの一連のシーケンスの組み合わせによシ、物体の
原子核のスピンの共鳴から緩和へ至る自由誘尋減衰信号
を観測し、これをフーリエ変換することによシ、画佇化
する手法が用いられている（フーリエ変換法）。

磁気共０）３イメージング装置のブロック図を第１２図
に示し、前記フーリエ変換法で用いられる典Ｍ的なパル
スシーケンスを第１３図に示す。第１３図のシーケンス
において、同図（イ）に示す９０°パルスと同図（ロ）
に示す傾斜磁場ＧＺはＺ軸に垂直な断層面（スライス）
を選択的に励起するために用いる。また、同図（イ）に
示す１８０°パルスと同図（ハ）に示す断層面内のＸ方
向の傾斜磁場Ｇｘは、励起されたスピンからの信号をエ
コー（Ｅｃｈｏ）の形で集めるためのもので、このとき
のエコーの周波数成分は断層面内のＸ方向に対するスピ
ンの分布に依存する。エコー信号は、マグネットボアの
中に配置され九ＲＦ受信コイルで検出された後、９０°
位相の異なる参照波で直交２相検波され、２チヤンネル
（Ｒｅａｔ、　Ｉｍａｇｉｎａｒｙ　）の検波信号とな
る。さらに、上記２つの傾斜磁場と互いに直交する勾配
磁場（第１３図（ニ）参照）ＧＹの振幅を、各９０°、
１８０°パルスの印加毎に正から負（又は負から正）へ
順以変化させて、Ｘ方向の位相エンコーディングを行う
。なお、第１２図の説明は第１図（実施例）と重複する
ので、こ工では省略する。

このようにして収集された２チヤンネルの検波信号は、
一般に以下のような式で表現できる。

Ｓ（ｔ、Ｇ）＝に／ｆＭ（ｘ＊ｙ）ｅｘｐ（ｉｒ（Ｇ−
ｙ−ｔ。

＋Ｇｘ−ｘ−１）〕ｄｘｄｙ　　　・・・・・・（１）
ただし、Ｍ（ｘ＊ｙ）’断層面内のスピン密度分布Ｇｚ
：ｘ方向傾剰磁場勾配（定数） Ｇ：ｙ方向傾斜磁場勾配（変数） ｔ：ｘ方向傾斜磁場印加時間（変数） ｔｙ：ｙ方向傾斜磁場印加時間（定数）ｋ：比例定数、
ｒ：核磁気回転比。

ｉ：Ｖ：１ である。したがって、２チヤンネルの検波信号全体Ｓ（
ｔ、Ｇ）を２次元フーリエ変換すれば、スピン密度分布
Ｍ（ｘ、ｙ）が求められる。しかし、実しプには検波信
号の位相ずれ、Ｘ方向傾斜磁場の動特性の影響によるＸ
方向の位置に比例した位相ひずみ、または静磁場の空間
的不均一性による位相ひずみなどによシ、信号の位相ず
れが生じる。この位相ずれｆｉ”ｔ−Δθとすると、（
１）式は次式の如く表わされる。

Ｓ’（ｔ、Ｇ）”ｋｆＪ’Ｍ（ｘ＊ｙ）ｅｘｐ（ｉｒ（
Ｇ−ｙ−ｔｙ（２）式は、フーリエ変換の位相シフト定
理よシ、ｒｔＩ’（Ｘ＊　ｙ　）−Ｍ（ｘ、　ｙ　）　
ｅ’Δ’　　　　　−・−（３）となシ、位相ずれによ
る画像の誤差が生じることになる。

一般に、現在用いられている磁気共鳴イメージング装置
（以下、ＭＨＩと記す）では、再構成された画像Ｍ（ｘ
、ｙ）の画素値は、絶対値で表現する方式が用いられて
おシ、シたがって、ＩＭ’（ｘ、ｙ）ｌ−ＩＭ（ｘ、ｙ
）ｅ’Δ’ｌ−ＩＭ（ｘ、ｙ）１・・・・・・（４） となシ、位相ずれの影響を消すことができる。すなわち
、フーリエ変換によシ求められたスピン密度分布Ｍ（ｘ
、ｙ）のＲｅａｊ（実軸）成分とＩｍａｇｉｎａｒｙ（
ｍ軸）成分の２乗和の平方根をとればよいわけである。

〔発明が解決しようとする問題点〕

しかし、そのためには計算機で各画素毎に２乗和の平方
根計算をしなければならず、これがスキャンが終了して
からデイスプレィにｌｉｌ！ｉｉ低が出てくるまでの時
間の短縮のための障害になっていた。

また、スピンエコー法のように、正の画素値をとるスキ
ャン方式は絶対値画像で問題はないが、インバージョン
リカバリー（反転回復）法のように最初に１８０°パル
スでスピンを反転してから信号を収集する方式では、負
の値をもつ画素も生ずるため、前記のような絶対値処理
では画像の正確な表現ができないという問題点があった
。

このような状況に対応するため、例えばインバージョン
リカバリーの画像についてはｓ　Ｒｅａ７成分の画ｇ、
をそのまま用いる方法もおるが、本質的な解決策ではな
い。また、収集された検波信号のＲｅａｔパートとＩｍ
ａｇｌｎａｒｙパートを用いて、計算処理によル位相角
を求め、位相補正を行う手法がいくつか考えられている
が、この場合に問題になるのは、ＲｅａｔとＩｍａｇｉ
ｎａｒｙの位相差を検出するために、各信号の簡対値の
最大値を見つける必要があることである。つまり、ＭＨ
Ｉの信号はノイズ成分が多く、単純な最大値検出方式で
は、かえって誤葺を大きくする可能性があると云うわけ
である。

本発明はこのような事情を鑑みてなされたものでアシ、
磁侃共鳴イメージングにおいて、コイルで検出され検波
された磁気共鳴信号のノイズ成分およびひずみ成分を除
去し、安定したｉｉ！ｌｉＩ似化が可能なようにすると
〜もに、不要な位相差（Δθ）を除去して反転回復法等
、負の画像値をもつシーケンスに対しても正確な画像化
を可能にすること、さらには大きなノイズ等による補正
不能な信号波形に対しても、再度スキャンすることなぐ
残シの正常信号波形から、画像化することができるよう
にすることを目的とするものである。

〔問題点を解決するための手段〕

磁気共鳴イメージングにおいて、一連のパルスシーケン
スによ）励起されたスピンが放出する信号をコイルで検
知し、増幅して位相敏感検波した後生成される検波信号
に対し、突出ノイズ（スパイク状信号）を検知するスパ
イク検出、及び信号のクリップ（あるいは飽和）状態を
検知するクリップ検出、ならびに信号を移動子均等によ
シ干滑化するスムージング手段と、信号の直流レベルを
検出し、補正するＤＣ補正手段と、信号の不要な位相差
成分を算出し全体の位相量から該算出した位相差を除去
する位相補正手段とを設ける。

〔作用〕

スムージング手段、ＤＣ補正手段お上び位相補正手段金
膜けてノイズ、直流レベルずれおよびタイミングずれを
除去するととＫよシ、安定かつ正確な画像化を可能にす
る。

〔発明の実施例〕

第１図は本発明の一実施例を示すブロック図である。１
は一様な静磁場Ｈａ　′ｆ？、発生する静磁場コイルで
あり、２は各種傾叫磁場（Ｇｘ、Ｇｙ、Ｇ工）の鈍化系
である。３は送受信用の高周波コイル系、４は高周波送
信回路系、５は高周波受信回路系、６は直交２相検波（
ＱＤ　’　Ｑｕａｄｒａｔｕｒｅ　Ｄｅｔｅｃｔｉｏｎ
）を行う検波回路系、７は（アナログ／ディジタル）Ａ
／Ｄ変換系、８はスパイク検出及びクリップ検出及び平
滑化を行う前処理系、９はＤＣ捕正処理系、１０は位相
補正処理系、１１は呵憾再構成系、１２は上記処理全体
を管理・制御するコンビエータ系、１６は画像を表示す
るデイスプレィ系である。

次に、上記処理装置の動作を説明する。

いま、一様な静磁場ＨＯを被検体Ｐにかけた状態で、送
信回路系４から送受信コイル５ｔ−経由して高周波（Ｒ
Ｆ）パルスを印加すると、被検体の特定原子核と高周波
との相互作用による磁気共鳴信号が送受信コイル３に誘
起されるが、このときスライス選択のための傾斜磁場Ｇ
Ｚｙ位相位相コンコーディングめの傾斜磁場ＧＹｐエコ
ー信号読み出しのための傾斜磁場Ｇｘ　ｔ−１傾斜磁場
発生系２によ〕必要なタイミングで発生させる。誘起さ
れた磁気共鳴信号は、送信回路系４からの信号ω。

を用いて検波回路系６で検波された後、Ａ／Ｄ変換系７
でディジタル信号化される。乙のディジタル信号は各種
のノイズを含んだ不安定なものであるため、以下に述べ
るような各種の検出、補正を行ったのち、画像再構成系
１１によシ＠匂化されてデイスプレィ系１３に画像とし
て表示される。

次に、本発明の本質的な部分である、ノイズを含んだ磁
気共鳴信号のノイズ検出及び補正処理について説明する
。まず、始めに磁気共鳴信号の突出ノイズ（スパイク状
信号）の検出、及び信号のクリップ（あるいは飽和）状
態の検知、ならびに平滑化を行う前処理系８の動作につ
いてのべる。

一般に、磁気共鳴信号は第２図に示すように、理想的に
は一定の周波数成分を含み、なだらかな包絡線Ｅをもっ
ている。しかし、システムの各部分から混入するノイズ
やシステムそのものの不安定性により、現実の磁気共鳴
信号は、第３人図に示すような単発的に突出した電圧と
してあられれるスパイク状ノイズＳＰが時々発生し、゛
また、被検体の動きやシステムの動的変動によシ、第３
Ｂ図に示すような信号のクリッピング（又は飽和状ＩＬ
Ｌ　）　ＣＬが発生する。また、スパイクやクリッピン
グが無くても、一般に現実の磁気共鳴信号は、第３Ｃ図
に示すようになめらかな波形にならず、小さいノイズ成
分を含んだ波形になシ、以後の位相補正ＪＰＤＣ補正等
の際の障害になシ、また究極的には画像のノイズ又はア
ーチファクト（擬像）としてあられれる。そこで、前処
理系８でこれらの検出及び除去を行う。

前処理系８では、まずＡ／Ｄ変換系７において、画像の
再構成に必要なデータサンプリング間隔よシも短いサン
プリング間隔でＡ／Ｄ変換を行う。

たとえば、画像の作成に必要なサンプリング間隔がｔ。

ｃμｓ〕であるとすると、本方式においては、実際のＡ
／Ｄ変換のサンプリングを（ｔｏ／ｎ）μＳ間隔で行う
。この場合、必要表データサンプリング点数よ）１倍多
いデータがサンプリングされるわけで、このデータを前
処理系８でｎ点を１点に変換する処理を行う。このとき
、前記の如きスパイクやクリッピングの検出、及びデー
タの平滑化を行う。以下、このだめの手法を記す。

第４Ａ図はｎ　＝　４の場合、すなわち４倍多いサンプ
リングを行った後、４点から１点にθ退させるときの例
である。こ〜では、例えば点Ｐ１を基準にして次の点Ｐ
２）その次の点Ｐ３、その次の点Ｐ４の各点との差分の
呟Δ１を求める。すなわち、各点Ｐ１の信号の値をｆ（
Ｐｘ）とすると、ｊＩ　＝　ｆ（Ｐｉ＋、）　−ｆ（Ｐ
、　）　　　　　・・曲（５）である。つぎに各Δｉの
平均をとシ、これを改めてＰ−ｔの値ｆ（Ｐｌ）に加え
ることによシ、平滑化されたＰｌの値を得ることができ
る。っぎの対象点はピッチ全４つずらした点Ｐ　ｓ　ｅ
　Ｐ　６１　Ｐ　７　＃　Ｐ　ｓであり、これらについ
ても前記と同様の処理を行う。なお、本方式においては
、平滑化の対象点（Ｐｌなど）に対して補正点（Ｐ２　
ｅ　Ｐｓ　ｓ　Ｐ４など）が時間軸で前にくるような例
を示したが、これは任意である。すなわち、平滑化対象
点ｅＰｓとし、補正点をＰ　１＊　Ｐ　２１　Ｐ　ａな
どの如くしてもよい。

以上のようにして平滑化処理が実現できるが、このとき
各点間の差分筺ΔｉＱ監視することによシ、スパイク検
出及びクリッピング検出を行う。

すなわち、スパイク検出については、システムの特性を
考慮して、一定の閾値αを設けて差分呟ΔＩの１絶対稙
１Δ口を監視し、 １ΔｉＤα　　　　　　　　　　　　・・・・・・（６
）の場合にスパイクが発生したと見なし、このときのΔ
ｉを前記平滑化演算の対象から除くようにする。

また、クリッピング検出は、上記と同様にΔｉを監視す
ることにより実行できる。すなわち、すべての歪につい
て １Δｌ−Δ汁、１＝０　　　　　　　　・・・・・・（
７）のとき、信号はクリッピング状態であると判断する
。

第４Ｂ図にスパイク状ノイズ検出方法の例、及び第４Ｃ
図にデータクリッピング検出方法の例を示す。データク
リッピングの場合の信号の処理については、その部分の
信号を除去するケース、また再度スキャンして信号をと
シ直すケースなどいろいろ考えられる。

次に、上記した前処理系８により平滑化された信号に対
する、ＤＣ補正系９によるＤＣ補正の処理について説明
する。

第５Ａ図に示すように、直交２相検波回路系６によシ検
波され、Ａ／Ｄ変換系７を通ってディジタル化されたＲ
ｅａｔ及びＩｙｌａｇｉｎａｒ）’信号は送受信コイル
３ｊＰ受信回路５．検波回路６．Ａ／Ｄ変換系７の回路
特性や伝送特性などによシ、直流レベルがシフトしてい
ることが多い。したがって、測定された信号のうち、直
流レベルと見なされる部分を検出して、この部分のレベ
ル金もとに信号全体から差し引くことによシ、直流成分
が実効的にゼロである信号を得ることができる。

直流レベルの検出は次のように行う。例えば、収集した
信号のうち時間軸で最も右寄りの部分は、共鳴信号がほ
とんど派衰しきった状態であシ、直流レベルと見なして
もよいことになる。この様子を第５Ｂ図に示す。同図の
人の部分が直流レベルと見なされる部分でアシ、この部
分のサンプリング点数をｎとし、ＤＣレベルをＳｄｃと
すると、・・・・・・（８） となシ、ＤＣ？１ｉ正後の信号Ｓ（ｔ、Ｇ）は、Ｓ（ｔ
　、　Ｇ　）　−Ｓ　（ｔ　＊　Ｇ　）　　Ｓ　ｄｃＯ
＜ｔ＜ｔｍａｘ　　　　　　　　・・・・・・（９）と
して、Ｒｅａｔ　＊　Ｉｍａｇｉｎａｒｙの各信号のＤ
Ｃ補正をおとガう。

次に、以上のべてきた前処理及びＤＣ補正された信号Ｓ
（ｔ、Ｇ）に対して、位相補正系１０によ）位相補正を
行う処理について、以下に説明する。

まず、第一に位相エンコード量がゼロのとき（Ｇ＝０）
の検波信号を用いて、システムに固有の、不要な位相差
Δθを算出する。

前記（２）式においてＧ＝Ｏとすると、Ｓ’Ｃｔｅ　Ｏ
）＝　ｋＪｆ　Ｍ（ｘｅ　ｙ）　ｅｘｐ（１ｒ（Ｇｘ−
ｘ−ｔ＋Δθ））ｄｘｄｙ　　　　　　　・・・・・・
（１０）とな、６、ｔ＝ｏのとき、Δθは上式より、と
なる。したがって、実際の検波信号Ｓ’（ｔ、０）から
１−０の点を見つけることができれば、（１１）式よシ
ネ要な位相ＭΔθを求めることができる。

（１０）式よ多信号の絶対亀は、 Ｉｓ’（ｔ、０）１ ＝　ｌｋＪｆＭ（ｘ、　ｙ）ｅｘｐ　（ムｒＧｘ−ｘ−
ｔ）ｄｘｄｙ−ｅｘｐ（ｉｒΔθ）１−１に」了Ｍ（ｘ
ｅ　ｙ）ｅｘｐ　（ｌｒ　Ｇ）（−ｘ−１）ｄｘｄｙ　
１−Ｉｓ’（ｔ、ｏ）１　　　　　　　　　　曲−（１
２）となり、位相量に依存しないことがわ必・る。ここ
で、信号ｓ’（ｔ　＊　ｏ）は前述したように、一般に
は多くのノイズを含んでいるので、ｔ；０の点を確実に
求めることはできない。しかし、前処理系８及びＤＣ補
正系９により平滑化、補正された信号Ｓ（ｔ、Ｇ）はノ
イズがカットされ、直流のシフト分も補正されているの
で、　　ｌ５（ｔ、０）Ｉは１−０のときに最大ｆｉ＋
ｆｋとると見なしてよいことになる。したがって、Ｉｓ
　（ｔ、　ｏ月が最大［直をとる点がｔ＝Ｕで必り、し
たがって、ｌ　Ｓ　（０，０月の肱が求めら扛る。

次に、（１１）式より求められたΔθを用いて、Ｇ〜０
（位相エンコーディング時）のときの信号からΔｔ／ｌ
−除去する。信号′ｆ！：実都、虚部にわけ、位相補正
後の信号、すなわち（１）式でろられされる信号の実部
ｐ　ＩＬｔＪｔ−Ｒ６（’ｉ（、ｔ、Ｇ））　＊　Ｉｍ
（Ｓ（ｔ−Ｇ））とする。また、前処理・ＤＣ補正後で
位相補正前の信号ｔｓ’（ｔｓＧ）とし、その実部、虚
部をそれＳ　’（ｔ、　Ｇ）−Ｒｅ（Ｓ’（ｔ、　Ｇ）
）＋　ｉ　Ｉｒ、ｔ（Ｓ’（ｔ、　Ｇ））＝１６（Ｓ（
ｔｅθ）〕ω３Δθ−Ｉｍ（Ｓ（ｔ　、θ））ｓｉｎθ
十目〕□（Ｓ（ｔ、θ）　：ｌ　ｓｉｎΔθ＋Ｒｅ（Ｓ
（ｔ、θ））ｃｏｓΔθ） ・・・・・・（１６） Ｒｅ［８（ｔｌＧ）、１ −　ｋｆｆＭ（ｘｅｙ）ｃｏｓ［ｒ（Ｇ）（−ｘ−を十
Ｇ−３−ｔ、））Ｃｉｘａｙ工□（Ｓ（ｔ、Ｇ）） ＝　ｋｆｆＭ（ｘｅ７）ｓｉｎＬｒ（Ｇｘ−ｘ・ｔ＋Ｇ
−ｙ−ｔｙ）〕ｄｘｄｙ・・・・・・（１４） したがって、 ・・・・・・（１５） このようにして、測定された磁気共鳴信号のＲｅａｔｔ
ｌｍａｇｉｎａｒｙの各信号から不要なΔθ全除くこと
ができる。

さらに、スピンの密度分布Ｍ（ｘ、ｙ）は実数であるか
ら、（１）式より、 ５（−ｔ、−Ｇ） −１ｃｆｆ　Ｍ（ｘｅ　ｙ）　ｅｘｐ　（−１ｒ（Ｇ−
ｙ−ｔｙ＋Ｇｘ−ｘ−ｔ）：１ｄｘｄｙとなシ、信号Ｓ
　（ｔ　、　Ｇ）において、Ｇの負の部分の信号とＧの
正の部分の信号とは複窓共役の関係になっておシ、第６
図に示すように、複考平面において点対称である。した
がって、この性質を用いると前記スパイク及びクリップ
検出系において、異常なスパイク又はクリップが検出さ
れて、補正不能と見られる信号（Ｂａｄ　　Ｖｉｅｗと
呼ぶ。）については、そのＶｉｅｗを再度スキャンする
ことなく、複素共役な関係にある信号でｅきかえること
により信号全体を補償することができ、スキャン時間の
延長を防ぐことができる。

次に、本発明の主要な部分である前処理系８、ＤＣ補正
系９、位相補正系１０に関し、今まで述べてきたことを
呆視するだめの具体的なブロック図又はフローチャート
について説明する。

第７図に前処理系８のハードウェアブロック図を示す。

こ〜でも、前記の説明と同様に、必要なサンプリング点
数の４倍のデータから、平滑化処理ｔ−施して１点に縮
退させる場合で記述する。すなわち、ＰＨ−ＰＩ＋３１
での共鳴信号をデータバッファ８１に格納し、そのうち
差分のための基準データとなるものを基準データバッフ
ァ８２へ格納する（例えば、Ｐｌとする）。次に、差分
器８３へ均＊ｐｉ＋１・Ｐｉ＋２１　Ｐｉ＋５を入力し
て差分随Δ五〜Δｉ＋２　ｔ−求め、このデータをクリ
ッピング検出系８４及びスパイク検出系８５へ転送する
。クリッピング検出系８４では、差分器Δｉ〜Δｉ＋２
　間のさらに差分をとる差分器８４Ａによシ差分直ヲ監
視し、差分器の差分がゼロの場合はカフ／り８４Ｂｔ−
インクリメントし、ゼロでない場合はカウンタをリセッ
トする。カウンタの直が一定値以上になるということは
差分器の差分がゼロ、すなわち、信号のレベルが一定の
状態がつづくことであシ、シたがって、カウンタにプリ
セットされた匝を越えた場合はクリップ検出として、ア
ラーム全発生するようにする。一方、スパイク検出系８
５は、前記差分器Δｉの絶対呟そのものを一定の基準直
と比較し、基準呟よシ大きい場合はアラーム信号（ｉ−
発生し、そＯときの差分器Δｉは棄却し、以後の平滑化
処理には用いないようにする。なお、クリップ検出系８
４及びスパイク検出系８５においてアラームが発生しな
い場合は、差分器Δｌ〜Δｉ＋２に加算器８６へ転送し
、次の割算器（ビットシフトによる）８７によシ平均化
処理ｋｔａした後、その（μを基準データバッファ（こ
の場合Ｐｉ）の便に加算することによシ、平滑化された
後のＰｉ′を得る。

次に、ＤＣ補正系９の処理のブロック図を第８Ａ図に示
す。こ瓦では、磁気共鳴信号の変化分から直流部分を検
出するために、前処理系８で平滑化された信号Ｐｉ′を
用いる。

すなわち、信号Ｐｉ′と次の平滑化信号Ｐｉ＋１′から
差分器９１でその差分器Δｉ′を求め、この絶幻匝１Δ
／＋をコンパレータ９２で基準（直Ｃと比較し、Ｃよシ
小さい場合はＤＣ成分データバッファ９３へ　ｐ　Ｉ／
のデータを格納する。Ｃよシ大きい場合は、ＤＣ成分デ
ータバッファの内容をすべてクリアする。ＤＣ成分デー
タバッファは、例えば１６点分のデータがそろったとこ
ろで直流部分が検出されたものと見表し、次の加算器９
４及び割算器（ビットシフト）９５で上記１６点分の信
号の平均直を求め、Ｉ　Ｖｉｅｗ　ＯＲｅａｔ又はＩｍ
ａｇｉｎａｒｙ信号の直流レベルＳｄｃとする。直流レ
ベルＳｄｃが求まった後は、第８Ｂ図に示すようなＩＶ
ｉｅｗ分のＲｅａｔ又はＩｍａｇｉｎａｒｙのデータメ
モリから直流レベルＳｄｃを引くことにょシ、ＤＣ補正
を行う。

次に、位相補正に関し、第９図及び第１０図にそのフロ
ーチャートを示す。第９図は位相差Δθを求めるフロー
でる）、第１０５４はΔθを用いて位相補正を行うフロ
ーである。これらは、この手順に従ってハードウェアで
実現することも可能であるが、ζこではフローを示すだ
けにとどめておく。

次に、第１１図に本発明の最後のポイントであるＢａｄ
　Ｖｉｅｗデータの置きかえを示す。

前述したように、前処」系においてデータのクリッピン
グやスパイクが検出された場合、これをそのまま再栴成
演算に用いるとノイズやアーチファクト（擬像）の原因
になる丸め、前処理系にてアラームが発生した場合は、
とのｖｔ　ｅｗのデータをそれと複素共役なデータのＶ
ｉｅｗで置きかえる処理全行う。

〔発明の効果〕

本発明によれば、磁気共鳴イメージング装置において、
システムの各種の要因からくる共鳴信号のノイズを除去
し、さらに信号を平滑化し、ＤＣＣ正正位相１００正を
行９ようにしたので、生成される画像のノイズやアーチ
ファクトが低減し、また突発的に生ずるＢａｄ　Ｖｌｅ
ｗによる画像の１−チファクトについても再スキャンす
ることなく、既に収集し九Ｖｉ　ｅｗのデータで置きか
えることによ）、システムの安定性を向上させることが
できる。

【図面の簡単な説明】

第１図は本発明の実施例を示すブロック図、第２図は理
想的な磁気共鳴信号の例を示す波形図、第３Ａ図は信号
にｍ畳したスパイク状ノイズの例を示す波形図、第３Ｂ
図は信号のクリッピングの例を示す波形図、第３Ｃ図は
小さいノイズを含む磁気共鳴信号の例を示す波形図、第
４Ａ図は本発明における信号のスムージング方法を説明
するための説明図、第４Ｂ図は本発明におけるスパイク
状ノイズの検出方法を説明するための説明図、第４Ｃ図
は本発明における信号クリッピングの検出方法ｔ−説明
するための説明図、第５Ａ図は検波信号の直流分を説明
するための波形図、第５Ｂ図は本発明における直流レベ
ルの検出方法を説明するための説明図、第６図は磁気共
鳴信号の対称性を説明するための説明図、第７図は本発
明による前処理系を示すブロック図、Ｑ’ｇ　Ｂ　Ａ図
は本発明によるＤＣ補正系を示すブロック図、第８Ｂ図
は本発明によるＤＣ補正方法を説明するための説明図、
第９図は位相ずれの算出過程を示すフローチャート、第
１０図は位相補正方法を示すフローチャート、第１１図
はＢａｄ　Ｖｉｅｗの置き換え処理を説明するための説
明図、第１２図は磁気共鳴イメージング装置の従来例を
示すブロック図、第１３図はパルスシーケンスの例を示
すタイムチャートである。 符号説明 １・・・・・・￥Ｐ磁場コイル、２・・・・・・傾θ）
磁場発生系、３・・・・・・高周波送受信コイル、４・
・・・・・高周波送信回路系、５・・・・・・高周波受
信回路系、６・・・・・・検波回路系、７・・・・・・
Ａ／Ｄ変僕系、８・・・・・・前処理系、９・・・・・
・ＤＣ補正処理系、１０・・・・・・位相補正処理系、
１１・・・・・・画像再構成系、１２・・・・・・制御
コンピュータ系、１３・・・・・・デイスプレィ系、８
１・・・・・・バッファ、８２・・・・・・基準データ
バッファ、８３，９１・・・・・・差分器、８４・・・
・・・クリッピング検出系、８４Ａ・・・・・・差分値
差分回路、８４Ｂ・・・・・・カウンタ、８５・・・・
・・スパイク検出系、８６．９４・・・・・・加算器、
８７，９５・・・・・・”；！１３ＸＲ５，９２・・・
・・・コンパレータ、９３・・・・・・ＤＣ成分データ
バッファ。 代理人　弁理士　並　木　昭　夫 代理人　弁理士　松　崎　　　　済 簿　１　図 育２　■ ・冒３Ｎ［Ｆ］ ＮＴ３Ｂ図 薯３６図 π４　Ａ　ｒＡ ｆｘＡＢ図 ＊ａＣ：、図 膚５Ａ■ Ｒｅａｌ　　　　　　　　　　　　　　　　　　　　　
　　　　　丁ｍａｇｉｎａｒｙ硼５Ｂ図 Ｗ６　図 逼　７【４ 薯８Ｂ口 １９　口 笛１１図 第１２【′３

Claims (1)

1. 【特許請求の範囲】 １）被検体に対し、一様な静磁場と線形の勾配をもつ静
   磁場と一定周期の高周波磁場パルスとを印加して得られ
   る実部、虚部の２成分からなる磁気共鳴信号を受信コイ
   ルを介して検出し、２相検波した後アナログ−ディジタ
   ル変換することにより、少なくとも特定の原子核のスピ
   ン密度および緩和時定数の分布を映像化する磁気共鳴イ
   メージング装置において、 前記ディジタル化された磁気共鳴信号をディジタル的に
   平滑化する平滑化手段と、 該磁気共鳴信号の実部、虚部の各直流分を検出しその各
   々の補正を行う直流分補正手段と、磁気共鳴信号の実部
   と虚部の位相差を求めてその補正を行う位相補正手段と
   、 を設け、ノイズ、直流レベルずれおよびタイミングずれ
   を除去して画像化することを特徴とする磁気共鳴イメー
   ジング装置。 ２）特許請求の範囲第１項に記載の磁気共鳴イメージン
   グ装置において、前記平滑化の段階で補正不能な信号成
   分については、他の信号との置き換え処理にて補完を図
   ることを特徴とする磁気共鳴イメージング装置。