JPH021235A - 磁気共鳴イメージング装置の画像再構成方式 - Google Patents

Abstract

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

Description

【発明の詳細な説明】 〔産業上の利用分野〕 本発明は核磁気共鳴現象を利用した生体内断層像撮影Ｍ
置に係り、特に傾斜磁場強度を時間的に変動させながら
共鳴信号を計測する高速イメージング法において、化学
シフト成分を分離して撮像するのに好適な画像再構成方
式に関する。

〔従来の技術〕

従来、磁気共鳴イメージング（ＭａｇｎｅｔｉｃＲｅｓ
ｏｎａｎｃｅ　Ｉｍａｇｊ、ｎｇ、以下ＭＲＩと酪す）
において化学シフト成分を分離して撮像するには、例え
ばマグネティック・レゾナンス・イン・メデイスン。

第４巻第５号、第４５２頁から第４６０頁（１，９８７
年）（阿ａｇｎｅｔｉｃ　Ｒｏｓｏｎａｎｃｅ　ｉｎ　
Ｍｅｄｉｃｉｎｅ、ｖｏ　Ｑ　４　。

Ｎα５ｐｐ４５２−４６０，１９８７）における１′大
きな磁場不均一の下での化学シフトイメージング″（”
Ｃｈｅｍｉｃａｌ−５ｈｉｆｔ　Ｉｍａｇｉｎｇ　ｗｉ
ｔｈ　ＬａｒｇｅＭａｇｎｅｔｉｃ　Ｆｉｅｌｄ　Ｉｎ
ｈｏｍｏｇｅｎｅｉｔｙ”　）と題する文献に論じられ
ている如く、一定の傾斜磁場の下で共鳴信号を計測する
２次元フーリエ変換法を用い、かつ複数の異なるタイミ
ングで計測した共鳴信号を２次元フーリエ変換して再構
成された断層像から各々の化学シフ１−成分像を求める
第１の方式が知られている。この方式は、異なる化学シ
フト成分は、わずかに異なる共鳴３周波数を持つことが
ら。

計測タイミングをずらすことにより各化学シフト成分像
に異なる１ケ相を与え、その差を利用して各成分像を分
離するものである。

一方、全く新しいイメージング方式として、例えばプロ
シーデインゲス・オブ・フィツス・ソサエテイオブマグ
ネティック・レゾナンス・イン・メデイスン、第１５６
頁から第１５７頁（１９８６年）（Ｐｒｏｃｅｅｄｉｎ
ｇｓ　　ｏｆ　　５　ｔｈ　　５ｏｅｉｅｔｙ　　ｏｆ
　　ＭａｇｎｅｔｉｃＲｅ５ｏｎａｎｃｅ　ｉｎ　Ｍａ
ｄｉｃｉｎｅ、ｐｐ、１５　Ｇ　−１５７，１９８６）
における１′新しい高速スキャン法″（“Ａ　Ｎｏｖｅ
ｌＦａｓｔ　Ｓｃａｎｎｉｎｇ　Ｓｙｓｔｅｍ”　）と
題する文献に論じられている如く、傾斜磁場強度を時間
的に変動させながら共鳴信号を計測し、撮像に要する時
間を大幅に短縮する方式が知られている。

また、このようにして計測された共鳴信号から画像を再
構成する一般的な方法として、相関法と呼ばれる方法が
、ソサエティ・オブ・マグネティック・レゾナンス・イ
ン・メデイスン第６回大会予稿集７８１頁（１９８７）
　　（Ｐｒｏｃｅｅｄｉｎｇｓ　ｏｆｔｈｅ　６ｔｈ　
５ｏｃｉｅｔｙ　ｏｆ　Ｍａｇｎｅｔｉｃ　Ｒｅ５ｏｎ
ａｎｃｅ　１ｎＭｅｄｊｃｉｎｅ、　ｐ　７８１　（１
９８７））における文献で論じられている。

さて、上記のような方法により、原理的に１回の共鳴信
号計測により、複数の化学シフト成分を同時に画像化す
る試みが、ソサエティ・オブ・マグネティック・レゾナ
ンス・イン・メデイスン第６回大会予稿集２３０頁０９
８７）（Ｐｒｏｃｅｅｄｉｎｇｓｏｆ　ｔｈｅ　６ｔｈ
　５ｏｃｉｅｔｙ　ｏｆ　Ｍａｇｎｅｔｉｃ　Ｒｅ５ｏ
ｎａｎｃｅ　ｉｎＭｅｄｉｃｉｎａ、　ｐ　２３０　（
１９８７））において論じられている。この第２の方式
は、２次元画像の撮影において、各次元に対応する２つ
の軸方向の傾斜磁場の両方とも時間的に強度を変化させ
た場合、化学シフト効果による共鳴周波数の違いが、画
像のぼけとなる事実を利用する。すなわち、水素原子核
に注目した場合、水に含まれる水素原子核の共鳴周波数
に合わせて画像を再構成すると、脂肪に含まれる水素原
子核による信号は、共鳴周波数の違いにより余分の位相
回転を受けるため、ぼけた画像として水成分の画像に重
置する。脂肪成分の共鳴周波数に合わせて画像を再構成
した場合も同様に、水成分がぼけた画像になる。いわば
、水成分、脂肪成分各々に焦点を合わせた画像が再構成
でき、その画像には、もう一方の成分の焦点はずれの画
像が混入する。この焦点は−ずれの画像を無視すれば、
水成分と脂肪成分の画像が同時に計測できることになる
。

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

上記第１の方式においては、化学シフト成分像を時間的
に変動する傾斜磁場を用いて高速に撮像することに関し
て考慮されていなかった。

化学シフト成分像をρＩ（ｘ、ｙ）（ｉは化学シフト成
分を表わす添字、Ｘ＋　ｙは断層面内の２次元座標）と
し、各成分の共鳴周波数の基準周波数からのずれをΔω
１とすると、共鳴信号Ｓ　（ｔ）は次の様に書ける。

５（ｔ）＝ＣΣｆｄｘｄｙρｔ（ｘ、　ｙ）・ｅｘｐ［
ｊｙ（ｋｘ（ｔ）ｘ＋ｋｙ（ｔ）ｙ）十ｊΔω＋１］　
　　　　　　・・・（１）ここで ｋ　ｘ（ｔ、）＝ｆ　Ｇｘ（ｔ’　　）ｄ　ｔ、’ｋｙ
（ｔ）＝ｆ　Ｇｙ（ｔ’　　）ｄ　ｔ、’Ｇｘ（ｔ）、
Ｇｙ（ｔ）：　Ｘ、ｙ方向の傾斜磁場強度γ：核磁気回
転比 Ｃ：比例定数 ５、Ｊ′Ｔ′Ｔ である。ただし、静磁場不均一性やＴｚ（よこ緩和時間
）効果は、以下の説明に無関係なため、省略した。ここ
で、２次元フーリエ変換法のように、Ｇｘ、Ｇｙが一定
１例えば０ｘ（ｔ）＝Ｇｘ。

Ｇｙ（ｔ）＝Ｏならば ５（ｔ）＝ＣΣｆ　ｄｘｄｙρ＋（ｘ　＋　　ｙ）　・
・　（２） となり、各成分は、Δω、／γＧｘ　だけすれた位置に
イメージングされる。しかしながら、Ｇｘ。

Ｇｙ共に時間的に変動する場合は、各化学シフト成分を
同時に再構成することができず、従来の第１の方式によ
る化学シフトイメージングの手法が適用できない。

本発明の第１の目的は、上記の問題点を解決し。

時間的に変動する傾斜磁場を用いた高速イメージング手
法で化学シフト成分像を得ることを可能とする画像再構
成方式を提供することにある。

他方、上記第２の方式は焦点が合った画像の点拡がり関
数（Ｐｏｉｎｔ　５ｐｒｅａｄ　Ｆｕｎｃｔｉｏｎ、以
下Ｉ）　Ｓ　Ｆと略す）と焦点はずれの画像のＰＳＦが
等しい平均パワーを持つという事実について考慮がなさ
れていなかった。すなわち、焦点はずれの画像を無視す
ることができず、化学シフト成分の同時画像化は不可能
であった。

式（１）の共鳴信号から２次元画像を再構成する方法に
はいろいろ知られているが、前述の相関法は、その中で
も最も一般的な方法であり、以下の説明でも相関法を＠
提として考える。

第に化学シフト成分に着目して再構成した画像をρｋＯ
’＋ｙ）とするとρ、は一般に複素数の値をとり ｐｍ（Ｘｓ幻＝７ｄｔＳ（ｔ）ｅｘｐ　［＋ｊ　ｙ（ｋ
ｘ（ｔ）ｘ＋ｋｙ（ｔ）ｙ）コ・ｅｘｐ［ｊΔωｋｔ］
・ω（ｔ）・・・（３）である。ここでω（１）はＰＳ
Ｆをコントロールするための重み係数である。

式（１）を式（３）に代入すると、 ｐＨ（ｘ、　ｙ）＝ＣΣｆ　ｄｘ’　ｄｙ’　ｐ＋（ｘ
’　＋ｙ′）・ｅｘｐ［ｊ　ｙ（ｋｘ（ｔ）（ｘ’　−
ｘ）＋ｋｙ（ｔ）（ｙ’　　ｙ））］・ｅｘｐ［＋ｊ　
　（Δω１−Δω、）　ｔコ　・　ω（１）＝（）、ｋ
ｋ＊ρｋ）（Ｘ、ｙ）＋Σ　（ｈｋｔ”ρ、）（ｘ、ｙ
）１≠に となる。ここで記号率は２次元畳み込み積分を表わし くｆ”ｇ）（ｘ、い＝ｆｄｘ’　　ｄｙ’　　ｆ（ｘ’
　　　ｙ’　　）ｇ（ｘ　　ｘ’　　＋ｙ−ｙ’　　）
ｅｆ ・・・（５） である。

式（４）でｈ＋＋＋はＰＳＦを表わし ｈｋｔ　（ｘ、　　ｙ）＝Ｃ／　ｄｔ　　ｅｘｐ［ｊ　
ｙ　（ｋｘ（ｔ）ｘ＋ｋｙ（ｔ）ｙ）コ・ｅｘｐ［＋ｊ
（Δω、−Δωｈ）ｔ］・ω（１）・・・（６） である。ｈｋｋは焦点の合ったＰＳＦ、ｈｋｔ　（ｋ≠
ｉ）は焦点はずれのｐ　Ｓ　Ｆを表わす。１１ｋｌは傾
斜磁場変動パターンおよび重み係数ω（１）から求めら
れ既知と考えてよい。

式（６）から ｆｄ　ｘｄ　ｙｌｈｋｋ（ｘ、ｙ）１２＝／ｄｘｄｙｌ
　ｈｋｔ（ｘ＋　　ｙ）１２ｆｏｒａＱ　　Ｑ　ｉ ・・（７） であることが分る。したがって式（４）におけるρ、で
、第に成分以外の成分画像の影響を無視することはでき
ない。水素原子核の場合のように、たとえ化学シフト成
分が２つしかないと見なせる場合でも同様である。

本発明の第２の目的は、式（４）における焦点はずれの
成分を推定・除去し、１回の計測から複数の化学シフ１
〜成分画像を得ることを可能とする画像再構成方式を提
供することにある。

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

上記第１の目的は、複数の異なるタイミングで計１１１
１１１　Ｌだ共鳴信号から、１つの化学シフト成分のみ
による共鳴信号を求め、その結果に対して、注目した化
学シフト成分固有の画像再構成アルゴリズムを適用する
ことにより達成される。

データ計ｄ１リタイミングを八Ｔだけずらした時の共鳴
信号を５ＡＴ（シ）は、 ＳΔｒ（ｔ）＝ＣΣｆｄＸｄｙρ、（ｘ、ｙ）ｅｘｐ［
−ｊ　ｙ（ｋ　ｘ（ｔ）ｘ＋ｋ　ｙ（ｔ）ｙ）＋ｊΔω
＋１］　　・ ｅｘｐ［ｊΔω、・ΔＴ］　　　　　　　　・・（８）
となる。（８）式の最後の因子は、各化学シフト成分の
寄与に、ｅｘｐ［ｊΔωオ・ΔＴ］なる位相を与えてい
ることになる。したがって、複数の八Ｔに対して共鳴信
号を計測し、それらの適当な線形和をとることにより、
第ｉ化学シフト成分のみによる共鳴信号Ｓ＋（ｔ） Ｓｔ（ｔ）＝Ｃｆｄｘｄｙρｔ（ｘ＋ｙ）ｅｘｐ［ｊ７
（ｋｘ（ｔ）ｘ十ｋｙ（ｔ）ｙ）ｔ＋ｊΔωＩｔ］　　
　　　　　・・（９）を求めることができる。

第１成分に対する画像再構成は、（１〕）式において、
共鳴周波数のずれによる位相回りを打消した共鳴信号Ｓ
＋（ｔ） Ｓｔ（ｔ）＝Ｓ＋（ｔ）・ｅｘｐ［ＪΔ（Ｊ　ｔ　ｔ］
　　　−・（１，０）を求めＳｔ（ｔ）を用いて画像再
構成を行う。

上記第２の目的は、各化学シフト成分に注目し画像を求
めることができる。

て再構成した画像ρ鰍を用い、ρｋに含まれる他の化学
シフ１−成分の影響を推定し、除去することにより達成
される。具体的には： （１）初期推定値を ρ　−−ρ　ｈ 画像ρにの近似画像と考える。式（１１）により、式（
４）におけるρ８の、第に成分以外の成分の影響が除去
でき、さらによい近似が得られる。これを（ｎ＋１）番
目の推定値ρ、　　　を １≠ｋ により求める。ただしＡｎは ここでρには第に化学シフト成分のスピン密度に比例す
るから、実数で、かつ非頁の値をとる。

式（１１）の関数Ａｎはθ。＝０のとき第ｎ番目の推う
条件を満足する成分を取り出すことに対応する。

なる関数とする。Ｒｅ　（Ｃ）は複索数Ｃの実数部を表
わしＱｎ（≧○）はｎによって定まるしきい値パラメー
タである。

実際には適当な終了条件を満たすまで式（１１）による
推定を繰り返すことにより、化学シフト成分より真の画
像ρｋに近いということができ、それだけ式（１１）に
よる推定の精度が向上する。

またｎ→大にしたがってθ。→小とすることにより、収
束速度を速めることができる。これはｎ：小のときは式
（４）の右辺第２項の寄与が太きいため、しきい値θ。

を大きくとれば１式（１１）による推定への寄与を小さ
くおさえておくことができるためである。

〔実施例〕

以下、本発明の第１の実施例を第１図〜第６図を用いて
説明する。

第２図は本発明を実施するＭＲＩ装置のブロック構成図
である。検査対象物から検磁気共鳴信号を検出するため
に、あらかじめ定められた手順に従って装置各部を制御
するシーケンス制御部２０１と、共鳴を起こさせるため
に発生する高周波パルスの送信器２０２と、傾斜磁場を
駆動する傾斜磁場１駆動部２０４およびそれを制御する
磁場制御部２０３と、検査対象物から発生する該磁気共
鳴信号を受信検波する受イ３器２０５と、画像再構成お
よび各種演算を行なう処理袋ｈ’ｌ　２０６と、画像表
示用ＣＲＴデイスプレィ２０７と、検出信号データ、再
構成画像データなどを記憶する外部記憶装置＆　２０８
とから成る。

以上の構成における本発明の第１の実施例のシーケンス
を第３図に示す。まずＺ方向の傾斜磁場パルス３０１の
印加と同時に９０°高周波磁場パルス３０２を発生させ
、撮影したいスライス内のスピン（撮影対象核種の原子
核）を共鳴させる。

次いでスピンエコーを発生させるための１８０゜高周波
磁場パルス３０３を、Ｚ方向傾斜磁場パルス３０４と同
時に印加する。その後１時間的に強度の変化するＸ方向
傾斜磁場３０５、Ｘ方向傾斜磁場３０６を印加しながら
、時間３０７の間、共１！θ信号を計測する。９０’高
周波磁場パルス３０２と１８０°高周波磁場パルス３０
３の時間間隔をＴＥ／２とすると１８０°高周波パルス
３０３印加後にＴＥ！／２後にスピンエコー信号が生成
される。Ｘ＋　Ｘ方向傾斜磁場３０５，３０６の印加お
よび共鳴信号の計ａｌｌｌは、スピンエコー生成時点か
らΔＴずれたタイミングで開始する。

第４図は、第３図におけるＸ＋”ｊ方向傾斜磁場３０５
．３０６によるｋｘ（ｔ）＋　ｋｙ（ｔ）の変動パター
ンの一例を示す。

とすれば、（Ａ、Ｂは定数） であり、（ｋ　ｘ　（ｔ）　ｐ　ｋ　ｙ　（ｔ）　）は
、第４図のような軌跡となる。黒点４０．１はサンプル
点で表わす。

第５図は、１８ｏ°高周波磁場パルスを使用せずに、傾
斜磁場反転によるエコー信号を計測する場合のシーケン
スである。第３図同様、まず２方向の傾斜＆１場パルス
５０１と同時に９０″高周波磁場パルス５０２を印加し
、撮影したいスライス内のスピンのみを共鳴させる。９
０６高周波磁場パルス印加からΔτ後にｙ方向傾斜磁場
５０３、Ｘ方向傾斜磁場５０４を印加しながら、時間５
０５の間共鳴信号を計測する。

第３図、第５図いずれのシーケンスを用いても以下の本
発明の実施手順は同様である。

第１図は、化学シフト成分が２成分のみの場合に本発明
の第１の実施例における画像再構成の処理手順を示すフ
ローチャートである。

ステップ１０１：第３図または第５図のシーケンスにお
いて、ΔＴ＝ΔＴ　１゜ ΔＴ２（ΔＴ１≠ΔＴ２）として共鳴信号ＳΔｒｔ（ｔ
）。

ＳΔｒｚ（ｔ）、　Ｏ≦ｔ≦Ｔを計測する。

ステップ１０２：次式に従って成分１．成分２各々の寄
与のみによる共鳴信号５ｚ（ｔ）、５ｚ（ｔ）を求める
。

ただし上式により５１（ｔ）、５ｚ（ｔ）が求まるため
には、（Δω２−Δω１）（八Ｔ２−ΔＴｌ）≠２ｎπ
（ｎ：整数）でなければならず、あらかじめΔＴ　ｓ　
。

ΔＴ２の値は、上の条件を満足する様に例えば、（Δω
２−Δω１）（ΔＴ２−ΔＴ１）＝πとなる様に選んで
おくものとする。

ステップ１０３：次式により、共鳴周波数のずれによる
位相回りを補正した共鳴信号を得る。

ステップ１０４　：　５（ｔ）より、成分１の画像ρｘ
（ｘ＋ｙ）を次式により再構成する。

ここで、Ｗ（ｔ）は点像応答関数を；ｉ整するためによ
り画像再構成を行う。

以上の手順により、化学シフト成分１，２の画像を分離
して再構成することができる。

第６図は、化学シフト成分が一般にｎ成分ある場合の処
理フローを示す。

ステップ６０１：第３図または第５図のシーケンスにお
いて、ΔＴ＝ΔＴ　ｚ　、ΔＴｚ、・・・、ΔＴｎとし
てｎ個の共Ｉｊ１信号ＳΔｒｄｔ）、ｘ＝＝１．”’　
　ｎｙＯ≦ｔ≦Ｔを計３１９する。

ステップ６０２：次式に従ってＳ、（ｔ）を算出する。

とすればよい。またＲｅ’［］は実数部をとる質算を表
す。

ステップ１０５：成分２に関しても、ステップ１０４と
同様 ここで行列Ａは で与えられる。ΔＴ　１　、・・・、ΔＴ、の値は逆行
列Ａ−１が存在する様に選ぶ。

ステップ６０３：次式により共鳴周波数のずれによる位
置回りを補正した共鳴信号を得る。

ステップ６０４：各成分ｉについて、共鳴信号Ｓ、（ｔ
）より画像ρ＋（ｘ、ｙ）を次式により再構成する。

以上の手順により、ｎ個の化学シフト成分像を分離して
再構成することができる。

以上の実施例を通じ、本発明の効果は傾斜磁場の時間変
動パターンや画像再構成のアルゴリズムには依存しない
ことは明らかであり、また３次元撮像の場合も、全く容
易に本発明の効果が得られることは言うまでもない。

また、画像再構成の手順は全て線形であるから計ｉ！＋
’ｌ信号上の演算、例えばＳΔ１１（１）からＳ、（ｔ
）を求める演算などは、画像再構成演算と順序を入れ替
えて、再構成画像の演算として実行しても、全く等価な
演算となり、したがって効果も変わる所がない。

つぎに、本発明の第２の実施例を第７図により説明する
。

第７図は、本発明の第２の実施例により水：Ｍ原子核を
対象としたＭＲＩにおいて、水／脂肪成分画像を分離す
る場合の処理手順を示すフローチャートである。

ステップ７０１：第３図のシーケンスにより共鳴信号Ｓ
（し）を計測する。

ステップ７０２：次式により水成分／脂肪成分ｅｘｐ［
ｊΔω１ｔ］ω（ｔ） ここで重み係数ω（１）は ｉ＝１．２ １／２ （、＋　（ｔ　）＝（ｋ　ｘ　（ｔ）２＋　ｋ　ｙ　Ｄ
　）”）とする。

ステップ７ｏ３：変数ｎの値をφとおく。

ステップ７０６：変数ｎの値に１を加える。

ステップ７０７：あらかじめ定められた定数Ｎ（例えば
Ｎ＝１０）と変数ｎの値を比較、ｎ　＜　Ｎならステッ
プ７０４へｎ≧Ｎならステップ７０８へ。

（Ｎ）　　　　（Ｎ） ステップ７０８：ρ１　、ρ２　を結果として出力し、
処理終了。

上記ステップ７０４において、しきい値パラメータＯ１
は、例えば次の様に定めればよい。

ここで、ｋｚｔは式（６）で与えられる焦点はずれＰＳ
Ｆ、Δ。は式（１２）で定義される関数である。

により求める。

ここでｍａｘ（）は最大値を表わす。

以上の実施例によれば、傾斜磁場強度の時間変動パター
ンや、画像再構成のアルゴリズムには依存せず本発明の
効果が得られることは明らかであり、また３次元撮像の
場合も、全く容易に本発明の効果が得られることは言う
までもない。

また、第７図のステップ７０４および７０５の処理は、
画像上の畳き込み積分で実行する代わりに計測信号上で
行なってもよい。すなわち、ステップ７０４の代わりに
、検出信号Ｓ（し）から。

＋ｋｙ（ｔ）ｙ）］　昏ｅｘｐ［＋ｊΔｗｅｔ］第１図
、第６第１木、第６第１の実施例における処理手順を示
すフローチャート、第２図はＭＨＩ装置のブロック構成
図、第３図および第５図は信号計測のためのパルスシー
ケンスを表わす図、第４図は空間周波数領域における信
号計測点の分布を示す図、第７図は本発明の第２の実施
例における処理手順を示すフローチャートである。

ツブ７０５も同様である。

またステップ７０４，７０５では、左辺の実数部分だけ
が次の繰り返しで必要となるため、右辺の計算でも実数
部だけを求めるようにして必要な計算量を削減すること
も可能である。

〔発明の効果〕

以上述べたように、本発明によれば、時間的に変動する
傾斜磁場を用いた高速撮像法において化学シフト成分像
を分離して再構成できるため、高速に化学シフト成分像
を撮像することが可能になる。

【図面の簡単な説明】

猜 図 遁 図 茅 図 第５図 第 図

Claims (1)

1. 【特許請求の範囲】 １、静磁場および傾斜磁場および高周波磁場の発生手段
   と、該発生手段を定められた手順に従つて制御する手段
   と、検査対象物からの磁気共鳴信号を検出する手段と、
   検出された信号に対して画像再構成を含む各種演算を行
   なう手段とを有する磁気共鳴イメージング装置において
   、複数の異なるタイミングで計測した共鳴信号の組から
   、１つの化学シフト成分のみの寄与による共鳴信号を推
   定することを特徴とする磁気共鳴イメージング装置の画
   像再構成方式。 ２、上記１つの化学シフト成分のみの寄与による共鳴信
   号の推定は、異なるタイミングで計測した共鳴信号の組
   の線形和によつて行う第１項の磁気共鳴イメージング装
   置の画像再構成方式。 ３、上記１つの化学シフト成分のみの寄与による共鳴信
   号を推定した結果に対し、該化学シフト成分固有の画像
   再構成アルゴリズムを適用する第１項または第２項の磁
   気共鳴イメージング装置の画像再構成方式。 ４、上記異なるタイミングで計測した共鳴信号の組に対
   し、１つの化学シフト成分固有の画像再構成アルゴリズ
   ムを適用する第１項の磁気共鳴イメージング装置の画像
   再構成方式。 ５、上記異なるタイミングで測定した共鳴信号の組のう
   ち、１つの化学シフト成分固有の画像再構成アルゴリズ
   ムにより再構成された画像の組を用いて、該化学シフト
   成分像を推定する第４項の磁気共鳴イメージング装置の
   画像再構成方式。 ６、上記化学シフト成分像の推定は、再構成画像の組の
   線形和によつて行なう第５項の磁気共鳴イメージング装
   置の画像再構成方式。 ７、静磁場および傾斜磁場および高周波磁場の発生手段
   と、該発生手段を定められた手順に従って制御する手段
   と、検査対象物からの磁気共鳴信号を検出する手段と、
   検出された信号に対し画像再構成を含む各種演算を行な
   う手段とを有する磁気共鳴イメージング装置において、 （１）時間的に強度の変化する傾斜磁場の存在下におけ
   る磁気共鳴信号を検出し、 （２）該検出信号に対し、複数の化学シフト成分の各成
   分に固有な画像再構成処理を適用することにより、各化
   学シフト成分に対応する複数の再構成画像を作成し、 （３）得られた複数の再構成画像から、各化学シフト成
   分画像を推定する処理からなることを特徴とする磁気共
   鳴イメージング装置の画像再構成方式。 ８、上記推定する処理は、 （１）特定の１つの化学シフト成分に着目し、（２）着
   目した成分以外の成分画像が、着目した化学シフト成分
   に対応する再構成画像上に及ぼす影響を推定し、除去す
   る処理からなることを特徴とする第７項の磁気共鳴イメ
   ージング装置の画像再構成方式。 ９、上記影響を推定する処理は、その影響を推定すべき
   化学シフト成分に対応する再構成画像と、所定の点拡が
   り関数との畳み込み積分を計算することにより行なうこ
   とを特徴とする第８項の磁気共鳴イメージング装置の画
   像再構成方式。 １０、上記影響を推定する処理は、その影響を推定すべ
   き化学シフト成分に対応する再構成画像の実数部分のみ
   を取り出した画像と、上記所定の点拡がり関数との畳み
   込み積分を計算することを特徴とする第８項の磁気共鳴
   イメージング装置の画像再構成方式。 １１、上記影響を推定する処理は、上記実数部分のみを
   取り出した画像の各画素毎に所定のしきい値との比較を
   行ない、該しきい値よりも小さな値をもつすべての画素
   の値を所定の値に設定した画像と、所定の点拡がり関数
   との畳み込み積分を計算することを特徴とする第１０項
   の磁気共鳴イメージング装置の画像再構成方式。 １２、上記影響を推定し、除去する処理は、複数回繰り
   返されることを特徴とする第８項の磁気共鳴イメージン
   グ装置の画像再構成処理方式。 １３、上記しきい値は、各繰り返しにおいて異なる値と
   することを特徴とする第１１項の磁気共鳴イメージング
   装置の画像再構成方式。 １４、上記影響を推定し、除去する処理は、検出信号に
   対する演算により行なうことを特徴とする第８項の磁気
   共鳴イメージング装置の画像再構成方式。