JPH0647021A - Ｍｒイメージングにおける動き補正方法 - Google Patents

Ｍｒイメージングにおける動き補正方法

Info

Publication number
JPH0647021A
JPH0647021A JP4206434A JP20643492A JPH0647021A JP H0647021 A JPH0647021 A JP H0647021A JP 4206434 A JP4206434 A JP 4206434A JP 20643492 A JP20643492 A JP 20643492A JP H0647021 A JPH0647021 A JP H0647021A
Authority
JP
Japan
Prior art keywords
data
motion
phase
line
measurement data
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
Application number
JP4206434A
Other languages
English (en)
Inventor
Kunihiko Kido
邦彦 木戸
Koichi Sano
耕一 佐野
Junichi Taguchi
順一 田口
Hideaki Koizumi
英明 小泉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP4206434A priority Critical patent/JPH0647021A/ja
Publication of JPH0647021A publication Critical patent/JPH0647021A/ja
Pending legal-status Critical Current

Links

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

(57)【要約】 【目的】計測データの上半分もしくは下半分のどちらか
一方のみで、突発的に発生するスライス面内における平
行移動性の動きを検出し、動き補正する方法を提供す
る。 【構成】計測データを中央低域と高域に分けて、それぞ
れ異なる方法で動きを検出し動き補正を行った。中央低
域の動き補正処理を行った後、高域に対する動き補正
を、計測データを上半分と下半分とに分けたとき、動き
の起こらなかった側の半分のデータから、動きの発生し
た側の半分のデータを複素共役点対称性を利用して推定
し、この推定したデータと動きの発生した側の半分の計
測データとの位相差をとり、この位相差データから動き
による位相誤差を、最小二乗法で推定することによって
動きを検出することにより動き補正処理を行った。これ
により画像のＳ／Ｎを劣化させることなく、動き補正す
ることが出来る。 【効果】上記の構成をとることで、画像のＳ／Ｎを劣化
させることなく、平行移動性のモーションアーチファク
トを抑制することが出来る。

Description

【発明の詳細な説明】
【０００１】
【産業上の利用分野】本発明は、核磁気共鳴現象を利用
した体内断層撮影に係り、特に撮影中の検査対象物の動
き、特にスライス面内の平行移動性の動きによって画像
に生ずるアーチファクトを抑制する方法に関する。
【０００２】
【従来の技術】ＭＲＩは撮影時間が長いために、検査対
象物の動きによるアーチファクトが発生しやすい。この
ため動きを補正する方法が多く提案されてきている。平
行移動性の動きに関しては代表的なものに、下記のもの
がある。
【０００３】(１) Gerchberg−Saxton アルゴリズムを
利用するもの エス エム アール エム，ブック オブ アブストラ
クト(SMRM，Book ofAbstracts）(１９９０）ｐ５６１ エス エム アール エム，ブック オブ アブストラ
クト(SMRM，Book ofAbstracts）(１９９１）ｐ７４３ (２) プロジェクションデータを利用するもの ナビゲーションエコーを使うもの ラジオロジー（Radiology）１７３，ｐ２２５−２６３
（１９８９） ナビゲーションエコーを使わないもの ラジオロジー（Radiology）１７９，ｐ１３９（１９９
１） ここで、平行移動の計測信号への影響をまとめておく。
まず、以下で提案された方法はすべて、動きが各位相エ
ンコード信号の間に発生したものを対象にしている。こ
れは各位相エンコード信号の計測時間が数ミリｓｅｃな
のに対して、計測待ちの時間が数百ミリｓｅｃと長いの
で、動きはこの計測待ちの時間に発生するものと仮定し
てよいという、いわゆるView−to−View効果にもとづい
ている。この時オリジナルの信号Ｓ(ｋx，ｋy)と、スラ
イス面内の平行移動による影響を含んだ信号Ｓ′(ｋx，
ｋy)には、次のような関係がある。
【０００４】
【数１】 Ｓ′(ｋx，ｋy)＝exp〔φ(ｋx，ｋy)ｉ〕Ｓ(ｋx，ｋy) …(１)
【０００５】
【数２】 φ（ｋx，ｋy)＝−２π(ｐ(ｋy)ｋx／Ｎx＋ｑ(ｋy)ｋy／Ｎy) …(２) ここで、ｋx＝−Ｎx／２〜Ｎx／２−１，ｋx＝−Ｎy／
２〜Ｎy／２−１，Ｎx，Ｎy はサンプル点数、ｐ
(ｋy)，ｑ(ｋy)は、ラインｋy−１とｋyの間で起こった
読みだし方向，位相エンコード方向の位置ずれ量であ
り、ｉは虚数単位である。
【０００６】すなわち平行移動は計測信号に線形な位相
誤差として現れる。
【０００７】(１)は検査対象物が有界な領域にあること
を利用するもので、まず画像において検査対象物の境界
線を決定して内部と外部に分ける。そして次のような手
続きを行う。
【０００８】 画像において、境界線から外部にはみ
出した部分は偽像として切り捨てる。
【０００９】 で得られた画像データを、逆フーリ
エ変換して信号データに直す。
【００１０】 で得られた信号データとオリジナル
の計測データとの位相差をとり、線形な位相誤差を位相
差データに最小二乗法を適用することにより推定する。
【００１１】 で推定された位相誤差により、オリ
ジナルの計測データに位相補正処理を行う。
【００１２】 で得られた信号データをフーリエ変
換して画像にする。
【００１３】 信号データの位相誤差が十分補正され
るまで以上の操作を繰り返す。
【００１４】この方法は、１ピクセル以下の小さな動き
に対して有効である。
【００１５】(２)は、読みだし方向の各ラインのデータ
を一次元フーリエ変換することによって得られるハイブ
リッドスペースデータが、各ライン計測時の読みだし方
向にそっての検査対象物の位置情報を持っていることを
利用するものであり、次の２通りの方法がある。
【００１６】 ナビゲーションエコーを使うもの 各位相エンコード信号の計測の間に、ナビゲーションエ
コーと呼ばれる位相エンコードがゼロの信号（ゼロエン
コード信号）を計測し、ナビゲーションエコーから得ら
れるハイブリッドスペースデータの絶対値から、エッジ
を検出することによって読みだし方向の位置ずれ量を検
出するものである。
【００１７】 ナビゲーションエコーを使わないもの Ｓ／Ｎが高い計測空間中央部の−３２〜＋３２ラインに
おいては、ナビゲーションエコーを使わなくても、読み
だし方向の各ラインに対するハイブリッドスペースデー
タの絶対値から、エッジを検出することにより、各計測
間の読みだし方向の位置ずれ量を精度良く検出すること
ができる。しかし、高域部分のラインに行くに従って計
測信号のＳ／Ｎが低くなるので、エッジの検出誤差が大
きくなり、高域部分では正確な位置ずれ量を検出するこ
とが困難になってくる。そこで、画像を再生するのに必
要なデータを計測する際に、高域部分のハイブリッドス
ペースデータにも明確に現れるような、小さなマーカー
を検査対象物の外部に据えて計測することによって、こ
のマーカーから読みだし方向の位置ずれ量を検出するも
のである。
【００１８】(２)は、１ピクセル以上の大きな動きに対
して威力を発揮する。
【００１９】上記(１)，(２)の方法の問題点として次の
ようなものがある。
【００２０】(１) Gerchberg−Saxton アルゴリズム
を利用するもの １ピクセル以上の大きな平行移動に関しては、かな
りの回数アルゴリズムを繰り返さないと効果が現れな
い。従ってアルゴリズム中にフーリエ変換が二回含まれ
ているため、効果が現れるまでに多くの計算時間が必要
となる。
【００２１】 画像上の検査対象物の境界線は手作業
で決定するか、もしくはリージョン・グローイング等の
技術により半自動的に決定しなければならず手間がかか
る。
【００２２】 (２) プロジェクションデータを利用するもの (ｉ) ナビゲーションエコーを使うもの 通常の画像を再生するためのエコーに加えて、ナビ
ゲーションエコーを計測するので計測時間が増える。
【００２３】 読みだし方向の平行移動だけしか補正
できない。
【００２４】 (ii) ナビゲーションエコーを使わないもの マーカーをつけても高域部分で数十ピクセルの検出
誤差が発生するる場合が多い。
【００２５】 読みだし方向の平行移動だけしか補正
できない。
【００２６】上記のように、どの平行移動補正法も長所
短所を合わせ持っているのが普通である。
【００２７】
【発明が解決しようとする課題】本発明は、突発的に起
こるスライス面内の平行移動を、特別なシーケンスを使
わずに、計測データ自身から動きを検出して動き補正処
理を行うものであり、１ピクセル以上の動きに対しても
効果がある。患者の突発的な動きは、患者が落ち着くま
での計測前半や、疲れてきてじっとしていられなくなる
計測後半など、計測中の限られた時間内に起こるのが普
通であり、計測空間の上下どちらか半分のみで発生する
ものと考える。例えば計測空間の上半分を計測中に平行
移動が発生したとき、体動発生ラインで一番低域側のラ
インをｊとすると、ラインｊよりも下側のデータ（ｋy
＜ｊ）は動きの影響を含んでいないと考えられ、しかも
全計測データの内の半分を含んでいるから、このデータ
から複素共役点対称性を利用して、残りの動きの影響を
含んだデータ（ｋy ≧ｊ）を推定することができ、動き
補正画像が得られる。しかし、オリジナルの計測データ
から動きの影響を含んだ部分を使わないので画像のＳ／
Ｎが劣化する。そこで、Ｓ／Ｎの劣化を抑えて動き補正
画像を得るためには、動きの影響を含んだデータ（ｋy
≧ｊ）に、動き補正する必要がある。
【００２８】数(１)，(２)により、ｋy ≧ｊにおいて読
みだし方向および位相エンコード方向の位置ずれ量ｐ
(ｋy)，ｑ(ｋy)が推定できれば、位相補正処理を行うこ
とにより、平行移動を補正することができる。従って、
本発明の目的は、このようなｐ(ｋy)，ｑ(ｋy)（ｋy ≧
ｊ）を検出することである。
【００２９】
【課題を解決するための手段】上記の目的を達成するた
めに次のような手続きを行う。
【００３０】〔Ｉ〕 図３のように、低域中央３０２と
高域３０１，３０３に分ける。そして、低域中央３０２
の動き補正のステップと高域３０１の動き補正のステッ
プに分ける。
【００３１】〔１〕 低域中央の動き補正 低域中央３０２で発生した平行移動に関しては、公知の
技術を使って行う。例えば、次のような方法で低域中央
３０２の動き補正を行うことができる。まず計測の際、
低域中央３０２を計測するときだけ、位相エンコード方
向の動きを検出するために、ｋy 軸３０４方向のゼロエ
ンコード信号（ナビゲーションエコー）｛ＮＡＶ
ky(ｙ)｝を計測する。
【００３２】 低域３０２の読みだし方向の各ライン
ｋy の計測データを一次元フーリエ変換して、ハイブリ
ッドスペースデータ｛ｈky(ｘ)｝を得る。 各ラインｋy の｜ｈky(ｘ)｜からエッジＥkyを検出
する。 ゼロエンコードデータのエッジを基準として、各ラ
インｋy ごとに読みだし方向の位置ずれ量ΔＸ(ｋy)＝
Ｅky−Ｅ0を求める。 ｐ(ｋy)＝ΔＸ(ｋy)，ｑ(ｋy)＝０として、数
(１)，(２)にもとづいて位相補正処理を行う。 ナビゲーションエコー信号｛ＮＡＶky(ｙ)｝につい
ても同様に、上記〜の手続きを行い、位相エンコー
ド方向の位置ずれ量ΔＹ(ｋy）を求める。 ｐ(ｋy)＝０，ｑ(ｋy)＝ΔＹ(ｋy)として、数
(１)，(２)にもとづいて計測データに位相補正処理を行
う。
【００３３】〔２〕 高域の動き補正 〔１〕で補正をかけた低域中央３０２とオリジナル
の計測データの高域３０３から、ハーフスキャンイメー
ジング法により高域３０１を推定する。 で推定した高域部３０１とオリジナルの計測デー
タの高域部３０１のデータの位相差をとる。 で求めた位相差データのノイズを除去するため
に、平滑化処理を行う。 で得た高域３０１の平滑化された位相差データに
対して、読みだし方向の各ラインにおいて数（２）のｐ
(ｋy)，ｒ(ｋy)＝ｑ(ｋy)ｋyを最小二乗推定する。 で推定したｐ(ｋy),ｑ(ｋy)により、数(１),(２)
にもとづいて〔１〕で補正をかけたデータに位相補正処
理を行う。
【００３４】〔II〕 〔Ｉ〕の〔２〕において、位相
差データに最小二乗法を適用するが、位相は２πごとに
サイクリックに同一の値をとるので、読みだし方向の１
ピクセル以上の動きに対して、図５（ａ）のように、各
ｋy における数（２）は折り返しが起こり直線にはなら
ない。従って最小二乗法を適用するには折り返しを解消
して直線になおすアンラップ処理を行う必要がある。し
かし実際の計測データは、ノイズの影響でアンラップ処
理が困難である。そこで、読みだし方向各ラインにおい
て、折り返し点を含まない区間を各ラインの低周波中央
部〔−ｃ，ｃ〕から次のように決定し、この区間で最小
二乗法を適用する。
【００３５】 読みだし方向の各ラインの位相差デー
タに移動分散を施す。ここで移動分散とは、ラインの各
点ｋx においてその点と前後ｂ点による位相差データに
よって分散を求めることである。 各ラインの区間〔−ｃ，ｃ〕において、しきい値６
以上の分散値をとる点ａ1＜ａ2＜…＜ａkを求める。 −ｃ，ａ1，ａ2，…，ａk，ｃ の隣接する２点どう
しの長さを求めて、最も大きい区間を最小二乗法の適用
区間とする。
【００３６】
【作用】
〔１〕 低域中央の動き補正 低域中央の動き補正は公知の技術であり、前記のとおり
Ｓ／Ｎが高い計測空間中央部の−３２〜＋３２ラインに
おいては、ナビゲーションエコーを使わなくても、読み
だし方向の各ラインに対するハイブリッドスペースデー
タを調べることにより、各計測間の読みだし方向の位置
ずれ量を精度良く検出することができる。また位相エン
コード方向の動き補正に関しては、ナビゲーションエコ
ーを使うことによって行うが、低域中央３０２を計測す
るときだけナビゲーションエコーを計測するので、計測
時間の増加は少なくてすむ。
【００３７】〔２〕 高域の動き補正 平行移動は計測空間の上半分のみで発生しているのだか
ら、計測空間の下半分の計測データは動きの影響を含ん
でいないと考えられる。画像が実数値を取る理想的な場
合には、計測データは原点に対して点対称な位置どうし
複素共役の関係にあるので、動きの影響のない下半分の
計測データを使うことにより、動きの影響を含んでいな
い上半分のデータを推定することができる。しかし実際
は、血流や装置歪み等で、画像に位相歪みが生じ複素成
分が現れるので、計測データの対称性は成り立たたず、
精度の良い上半分のデータを推定するができない。公知
のハーフスキャンイメージング法は、計測データの低域
中央部から得られる位相マップにより画像の位相歪みを
近似して、これを補正することにより計測データの領域
３０１と３０２の対称性を回復させるものである。
【００３８】従って、〔１〕で動き補正をした低域中央
１０２のデータと動きの影響のない高域３０３のデータ
により、ハーフスキャンイメージング法を使うことによ
って、動きの影響を含んでいない精度の良い高域３０１
のデータが推定できる。ここで、オリジナルの計測デー
タの高域３０１とハーフスキャンイメージング法により
推定した高域３０１のデータの位相差をとり、読みだし
方向各ラインｋy において位相差データに最小二乗法を
適用することによって、数（２）の位相誤差を推定する
ことができ、位相補正処理を行うことによって、平行移
動の補正をすることができる。前記のとおり、得られた
位相補正画像は全計測データを使うので、ハーフスキャ
ンイメージング法によって得られる動き補正画像よりも
Ｓ／Ｎが優れている。
【００３９】前記の通り、最小二乗法を適用する際に、
位相の折り返しの影響を避けるために最小二乗法の適用
区間を制限することを考える。もし検査対象物の位相エ
ンコード方向の動きがなく、位相にオフセットが乗って
いない場合には、最小二乗法の適用区間を読みだし方向
の各ラインに対して一定に、しかもライン中央低域にと
ればよい。しかし、位相エンコード方向の動きがある場
合には、数（２）からわかるとおり、ラインが高域に行
くに従って、位相オフセットが大きくまたは小さくなっ
てくるので、図５（ｂ）のように折り返し点５０１が中
央によってくる。このために各ライン毎に折り返し点５
０１を検出して、折り返し点を含まない区間を求め、こ
の区間で最小二乗法の適用する必要がある。
【００４０】折り返し点を検出するためにまず、各ライ
ン毎の位相差データに前記定義した移動分散をとる。折
り返し点付近の位相値は、ノイズの影響による振動に加
えて折り返しによる−πからπへの振幅が加わるので、
折り返し点付近の分散値（２πに近い値）はその他の点
の分散値よりも大きい値をとる。従って、適当なしきい
値を定め、このしきい値以上の分散値をとる点が折り返
し点の候補である。各ラインの高域部分は、ノイズの影
響が強いので、最小二乗法の適用区間から除外する。従
ってライン中央の区間〔−ｃ，ｃ〕を設定して、この区
間内から最小二乗法の適用区間を決定する。この区間に
おいて、分散がしきい値以上になる点をａ1，ａ2，…，
ａl とすると、折り返し点付近にこのような点が局在す
るので、この付近の隣接する２点どうしの距離は短い。
よって、−ｃ＜ａ1＜ａ2＜…＜ａl＜ｃ の隣接する２点
どうしの距離を求めて一番大きい区間を選べば、この区
間に折り返し点は含まれていない。よって、この区間を
最小二乗法の適用区間とすることにより、最小二乗法適
用時の折り返し点の影響を避けることができる。
【００４１】
【実施例】
＜実施例１＞図２は、本発明を適用しうるＭＲＩシステ
ムの一例を示す。静磁場発生系201は、均一な静磁場を
発生し、送信系２０２は、スピンを励起するための高周
波パルス磁場を発生する。傾斜磁場発生系２０３は、
ｘ，ｙ，ｚ方向の各傾斜磁場を発生し、各方向の磁場を
変化させることができる。受信系２０４は、被検体200
から放射される電磁波を受信して、検波し、デジタル信
号に変換した後、処理装置２０５に供給する。処理装置
２０５は、受信系２０４からのデータに各種の演算を施
して、画像データを生成し、ＣＲＴ表示装置２０６にそ
の画像を表示する。パルスシーケンスファイル２０７
は、系２０１〜２０４の動作シーケンスを規制する制御
情報を保持しており、シーケンス制御部２０８は、パル
スシーケンスファイル２０７から制御情報を読みだし、
それに従って、計２０１〜２０４の動作を制御する。
【００４２】図１３は、図２のシステムにおけるパルス
シーケンスを示し、これは、パルスシーケンスファイル
２０７からの制御情報に従って、シーケンス制御部２０
８の制御のもとに発生される。ＲＦ１３０１は、送信系
２０２により発生させられる高周波パルス磁場のタイミ
ングを示し、Ｇｚ１３０２，Ｇｙ１３０３，Ｇｘ1304
は、それぞれ、ｚ，ｙ，ｘ方向の傾斜磁場が傾斜磁場発
生系２０３により発生されるタイミングを示し、Signal
１３０５は、受信系２０４が計測データ信号1307を計測
するタイミングを示す。ＲＦ１３０１の周波数とＧｚ１
３０２のレベルでｚ方向のスライス面を選択し、Ｇｙ１
３０３のレベルでｙ方向の位置分離を行い、Ｇｘ１３０
４のレベルでｘ方向の位置分離を行う。ｙ方向の位置分
離は一度にはできないので、Ｇｙ１３０３のレベルを変
えて、それぞれのレベルについて計測データ信号１３０
７を計測する。通常、Ｇｙ１３０３は２５６レベルにわ
たって変化し、したがって、一回のスキャンの間に２５
６回の計測が行われる。このようにスキャンの間にレベ
ルが変化される傾斜磁場(この例ではｙ方向傾斜磁場)
が、位相エンコードパルス１３０６である。図１３のパ
ルスシーケンスにおいて、位相エンコードレベルを逐次
増加させるか減少させることによって計測されたデータ
をＳ(ｋx，ｋy）とし、ｋx＝−Ｎx／２〜Ｎx／２−１，
ｋy＝−Ｎy／２〜Ｎy ／２−１とする。
【００４３】図１は、計測空間の上半分か下半分のどち
らか一方だけで、検査対象物がスライス面内で平行移動
を起こした時の動き補正処理の流れを示すものであり、
ステップ１０１で、計測空間を低域中央３０２と高域３
０１と３０２に分ける。ここで、低域中央３０２を−Ｎ
〜＋Ｎラインとする。通常、Ｎは８，６，３２のいずれ
かにとるのが普通である。次に、ステップ１０２で低域
中央３０２の動き補正した後、ステップ１０３で高域の
動き補正をするが、計測データ自身から、発生した動き
が計測データの上半分で起こったものなのか、下半分で
起こったものなのか決定できないので、発生した動きが
計測空間の上半分で起こったものとしての動き補正処理
と、下半分で起こったものとしての動き補正処理を並列
して行う。従って、ステップ１０４では、動き補正処理
をした２つのデータをフーリエ変換することによって、
２枚の再生画像が得られる。以下最初に、低域の動き補
正ステップ１０２を、動きの検出方法の異なる３つの例
について説明する。
【００４４】〔低域中央の動き補正処理例１〕低域中央
の動き補正処理例１を図６により説明する。この例で
は、最初に低域中央の計測データを読みだし方向に一次
元フーリエ変換して得られるハイブリッドスペースデー
タからエッジを検出し、読みだし方向の位置ずれ量を検
出して読みだし方向の動き補正を行う。この処理によっ
て、読みだし方向の動きによる計測データの位相成分の
傾きがフラットになり、各ライン毎の位相エンコード方
向の動きによる位相オフセットの変化を、読みだし方向
の動き補正処理をしたデータから、各ライン毎の位相の
平均を調べることにより検出できる。すなわち、ゼロエ
ンコードラインの位相の平均を基準として、各ライン毎
の位相の平均の変化量を検出することにより、位相エン
コード方向の位置ずれ量を検出できる。この方法が有効
であるためには、装置歪みによる計測データの位相成分
に一様に乗っている線形な位相歪みを、予め除去してお
く必要がある。
【００４５】〔ステップ６０１〕装置歪みによる線形な
位相歪みを除去したデータに対して、低域中央３０２に
属する各ラインｋy(ｋy＝−Ｎ〜Ｎ)を一次元フーリエ変
換して、ハイブリッドスペースデータ｛ｈky(ｘ)｝を得
る。
【００４６】〔ステップ６０２〕｛｜ｈky(ｘ)｜｝から
エッジを検出して、各ｋy(ｋy＝−Ｎ〜Ｎ)における位置
ずれ量ΔＸ(ｋy)を求める。
【００４７】〔ステップ６０３〕 ｐ(ｋy)＝ΔＸ(ｋy)，ｑ(ｋy)＝０ ｋy∈〔−Ｎ，Ｎ〕 として、数(１)，(２)にもとづいて区間〔−Ｎ，Ｎ〕で
計測データＳ(ｋx，ｋy)に位相補正処理を施したデータ
Ｓ′(ｋx，ｋy)を求める。
【００４８】〔ステップ６０４〕低域中央３０２に属す
る各ラインｋy について、ステップ６０３で求めたデー
タＳ′(ｋx，ｋy)の位相成分ａｒｇ(Ｓ′(ｋx，ｋy))の
平均ｗ(ｋy）を計算する。
【００４９】〔ステップ６０５〕低域中央３０２に属す
る各ラインｋy に対して、位相エンコード方向の位置ず
れ量を ΔＹ(ｋy)＝(ｗ(ｋy)−ｗ(０))／ｋy として計算する。
【００５０】〔ステップ６０５〕 ｐ(ｋy)＝０，ｑ(ｋy)＝ΔＹ(ｋy) ｋy∈〔−Ｎ，Ｎ〕 として、数式(１)，(２)にもとづいて区間〔−Ｎ，Ｎ〕
で位相補正処理をステップ６０３で求めたデータＳ′
(ｋx，ｋy)に施す。
【００５１】〔低域中央の動き補正処理例２〕低域中央
の動き補正処理例２を図７によって説明する。この例で
は、低域中央３０２における位相エンコード方向の動き
を検出するためにナビゲーションエコーを利用する。位
相エンコードレベルが−Ｎ〜Ｎときのみ図１４で示され
るパルスシーケンスを使い、画像再生に必要な低域中央
データに加えて、ｋy 軸３０４方向のゼロエンコード信
号（ナビゲーションエコー）を計測する。ここで、ナビ
ゲーションエコー信号を｛ＮＡＶk(ｋy)}(ｋ＝−Ｎ〜
Ｎ）とする。読みだし方向の動きに対しては、低域中央
の読みだし方向各ラインの計測データを一次元フーリエ
変換して得られるハイブリッドスペースデータからエッ
ジを検出することにより、位相エンコード方向の動きに
対しては、ナビゲーションエコーを一次元フーリエ変換
して得られるハイブリッドスペースデータからエッジを
検出することによって、それぞれ位置ずれ量を検出す
る。
【００５２】〔ステップ７０１〕低域中央３０２に属す
る各ラインｋy について計測データを一次元フーリエ変
換して、｛ｈky(ｘ)｝を得る。
【００５３】〔ステップ７０２〕｛｜ｈky(ｘ)｜｝から
エッジを検出して、各ｋy(ｋy＝−Ｎ〜Ｎ)における位置
ずれ量ΔＸ(ｋy)を求める。
【００５４】〔ステップ７０３〕各ｋ(ｋ＝−Ｎ〜Ｎ）
に対して、ＮＡＶk(ｋｙ)を一次元フーリエ変換してハ
イブリッドスペースデータ｛ＨＮＡＶk(ｙ)｝を得る。
【００５５】〔ステップ７０４〕｛｜ＨＮＡＶk(ｙ)
｜｝からエッジを検出して、各ｋy(ｋy＝−Ｎ〜Ｎ）に
おける位置ずれ量ΔＹ(ｋy)を求める。
【００５６】〔ステップ７０５〕 ｐ(ｋy)＝ΔＸ(ｋy)，ｑ(ｋy)＝ΔＹ(ｋy) ｋy∈〔−Ｎ，Ｎ〕 として、数(１)，(２)にもとづき区間〔−Ｎ，Ｎ〕で計
測データＳ(ｋx，ｋy)に位相補正処理を行う。
【００５７】〔低域中央の動き補正処理例３〕この例で
は、マルチスライスの場合を考える。図１５はマルチス
ライススピンエコーシーケンスの一例である。この例で
は、第２スライス撮影時に位相エンコード方向と読みだ
し方向を、他のスライス撮影時の位相エンコード方向と
読みだし方向と入れ替えてある。これによって、各スラ
イスの低域中心のエコー計測間に発生した平行移動に関
して、図３の計測空間におけるｋx 軸３０５方向の動き
を第１スライスのエコーから得られるハイブリッドスペ
ースデータから、またｋy 軸３０４方向の動きに対して
は、第２スライスのエコーから得られるハイブリッドス
ペースデータから位置ずれ量を検出できる。各スライス
から得られるスピンエコー信号を Ｓ１(ｋx，ｋy)，Ｓ２(ｋx，ｋy)，…，Ｓｎ(ｋx，ｋy) とする。以下、低域の動き補正処理例３を図８により説
明する。
【００５８】〔ステップ８０１〕第１スライスの各ｋy
＝−Ｎ〜Ｎエコーによる計測データＳ１(ｋx，ｋy)を一
次元フーリエ変換して、｛ｈ１ky(ｘ)｝を得る。
【００５９】〔ステップ８０２〕｛｜ｈ１ky(ｘ)｜｝か
らエッジを検出して、各ｋy(ｋy＝−Ｎ〜Ｎ)におけるｋ
x 軸３０５方向の位置ずれ量ΔＸ(ｋy)を求める。
【００６０】〔ステップ８０３〕第２スライスの各ｋx
＝−Ｎ〜Ｎエコーによる計測データＳ２(ｋx，ｋy)を一
次元フーリエ変換して、｛ｈ２kx(ｙ)｝を得る。
【００６１】〔ステップ８０４〕｛｜ｈ２kx(ｙ)｜｝か
らエッジを検出して、各ｋy(ｋy＝−Ｎ〜Ｎ)におけるｋ
y 軸３０４方向の位置ずれ量ΔＹ(ｋy)を求める。
【００６２】〔ステップ８０５〕 ｐ(ｋy)＝ΔＸ(ｋy)，ｑ(ｋy)＝ΔＹ(ｋy) ｋy∈〔Ｎ，Ｎ〕 として、数(１),(２)にもとづいて区間〔−Ｎ，Ｎ〕に
位相補正処理をＳ１(ｋx,ｋy)，Ｓ２(ｋx，ｋy)，…，
Ｓｎ(ｋx，ｋy)に施す。
【００６３】以下、低域の動き補正処理を行ったデータ
Ｓ１(ｋx,ｋy)，Ｓ２(ｋx,ｋy),…,Ｓｎ(ｋx，ｋy)それ
ぞれに、高域の動き補正処理を行う。
【００６４】次に、高域の動き補正ステップ１０３に関
して説明する。高域の動きの検出は、高域３０１と３０
３の計測データが、互いに複素共役点対称の関係にある
ことを利用する。まず、高域３０１のデータを高域３０
３の計測データから推定し、同様に高域３０３のデータ
を高域３０１の計測データから推定する。そして、この
推定した高域のデータと高域の計測データとの位相差を
とり、平行移動による位相誤差を位相差データに最小二
乗法を適用することによって求める。ここで、位相差デ
ータに最小二乗法を適用する時に、推定精度を上げるた
めに、位相差データにノイズの除去処理を行うが、以下
２つの高域の動き補正処理例では、それぞれ異なる方法
を適用する。
【００６５】〔高域の動き補正処理例１〕高域の動き補
正処理例１では、位相差データのノイズの除去処理を、
各ラインの位相差データに移動平均処理によって行う。
以下、高域の動き補正処理例１を図９により説明する。
【００６６】〔ステップ９０１〕低域の動き補正ステッ
プ１０２で動き補正した低域中央３０２のデータと高域
３０３のデータから、公知のハーフキャンイメージング
法により、高域３０１の推定データＩＵ(ｋx，ｋy)(ｋy
≧Ｎ＋１）を求める。同様にして低域中央３０２のデー
タと高域３０１のデータから高域３０３の推定データＩ
Ｄ(ｋx，ｋy)(ｋy≦−Ｎー１）を求める。
【００６７】〔ステップ９０２〕Ｓ(ｋx，ｋy)の高域３
０１と推定データＩＵ(ｋx，ｋy)との位相差をとる。
【００６８】 ＰＵ(ｋx，ｋy)＝arctan(Ｓ(ｋx，ｋy)＊conj（ＩＵ(ｋx，ｋy))） （ｋy≧Ｎ＋１） 同様にして、Ｓ(ｋx，ｋy)の高域３０３と推定データＩ
Ｄ(ｋx，ｋy)との位相差をとる。
【００６９】 ＰＤ(ｋx，ｋy)＝arctan(Ｓ(ｋx，ｋy)＊conj（(ＩＤ(ｋx，ｋy))） （ｋy≦−Ｎー１） ここで、conj( )は複素共役を表わす。
【００７０】〔ステップ９０３〕各ラインｋy(ｋy≧Ｎ
＋１，ｋy≦−Ｎ−１）における位相差データＰky(ｋx)
＝ＰＵ(ｋx，ｋy）またはＰＤ(ｋx，ｋy）に対して、Ｍ
ＰＵ(ｋx，ｋy）または ＭＰＤ(ｋx，ｋy)＝(Ｐky(ｋx−ｂ)＋Ｐky(ｋx−ｂ＋１)＋…＋Ｐky(ｋx)＋ …＋Ｐky(ｋx＋ｂ−１)＋Ｐky(ｋx＋ｂ))／(２ｂ＋１) （ｋx＝−Ｎx＋ｂ〜Ｎx−ｂ） として、位相差データに平滑化を行う。
【００７１】〔ステップ９０４〕以下で説明する最小二
乗法適用区間決定ステップに行き、高域３０１の各ｋ
y(ｋy≧Ｎ＋１）において、ＭＰＵ(ｋx，ｋy)の最小二
乗法適用区間〔ｓｕ(ｋy),ｅｕ(ｋy)〕を求める。同様
にして、高域３０３の各ｋy(ｋy≦−Ｎ−１）におい
て、ＭＰＤ(ｋx，ｋy)の最小二乗法適用区間〔ｓｄ
(ｋy)，ｅｄ(ｋy)〕を求める。
【００７２】〔ステップ９０５〕高域３０１の各ｋy(ｋ
y≧Ｎ＋１）において、区間〔ｓｕ(ｋy)，ｅｕ(ｋy)〕
でＭＰＵ(ｋx，ｋy）に最小二乗法を適用し、数（２）
のｐ(ｋy)，ｒ(ｋy)＝ｑ(ｋy)ｋy(ｋy≧Ｎ＋１）を推定
する。
【００７３】同様にして、高域３０３の各ｋy(ｋy≦−
Ｎ−１）において、区間〔ｓｄ(ｋy)，ｅｄ(ｋy)〕でＭ
ＰＤ(ｋx，ｋy)に最小二乗法を適用し、数（２）のｐ
(ｋy)，ｒ(ｋy)＝ｑ(ｋy)ｋy(ｋy≦−Ｎ−１）を推定す
る。
【００７４】〔ステップ９０６〕ステップ９０５で推定
したｐ(ｋy)，ｑ(ｋy)(ｋy≧Ｎ＋１）により、数(１)，
（２）にもとづきｋy≧Ｎ＋１ において、Ｓ(ｋx，ｋy)
に位相補正処理を行ったデータＳＵ(ｋx，ｋy)を得る。
【００７５】同様にして、ステップ９０５で推定したｐ
(ｋy)，ｑ(ｋy)（ｋy≦−Ｎ−１)により、数(１)，(２)
にもとづきｋy≦−Ｎ−１において、Ｓ(ｋx，ｋy）に位
相補正処理を行ったデータＳＤ(ｋx，ｋy)を得る。
【００７６】〔最小二乗法適用区間決定ステップ〕最小
二乗法適用区間決定ステップを図１０を使って説明す
る。
【００７７】〔ステップ１００１〕入力として、各ライ
ンｋy における位相差データＰＵ(ｋx，ｋy）またはＰ
Ｄ(ｋx，ｋy)を受け取る。
【００７８】〔ステップ１００２〕各ラインｋy におけ
る位相差データＰｋy(ｋx)＝ＰＵ(ｋx，ｋy）またはＰ
Ｄ(ｋx，ｋy）に対して、移動分散σ(ｋx)(ｋx＝−Ｎx
＋ｂ〜Ｎx−ｂ）をとる。
【００７９】 μ(ｋx)＝(Ｐｋy(ｋx−ｂ)＋Ｐｋy(ｋx−ｂ＋１)＋…＋Ｐｋy(ｋx) ＋…＋Ｐｋy(ｋx＋ｂ−１)＋Ｐｋy(ｋx＋ｂ))／(２ｂ＋１） σ（ｋx)＝((Ｐｋy(ｋx−ｂ)−μ(ｋx))2＋(Ｐｋy(ｋx−ｂ＋１)−μ(ｋx))
2 ＋…＋(Ｐｋy(ｋx)−μ(ｋx))2＋…＋(Ｐｋy(ｋx＋ｂ−１) −μ(ｋx))2＋(Ｐｋy(ｋx＋ｂ)−μ(ｋx))2)／(２ｂ＋１) 〔ステップ１００３〕移動分散値σ(ｋx）が〔−ｃ，
ｃ〕において、しきい値６をこえるｋx を求め、−ｃ＜
ａ1＜ａ2＜…＜ａl＜ｃとする。
【００８０】〔ステップ１００４〕length(０)＝ａ1＋
c，length(ｉ)＝ａi+1−ａi(ｉ＝１〜ｌ−１)，length
(１)＝c−ａlとする。このときlength(ｉ)（ｉ＝１〜
ｌ）が最大となるｉを求めimaxとする。
【００８１】〔ステップ１００５〕出力として、〔ａ
imax，ａimax+1〕を戻す。ここでimax＝０ ，imax＋１
＝l＋１のとき、それぞれａimax＝−ｃ,ａimax+1＝ｃと
する。
【００８２】〔高域の動き補正処理例２〕高域の動き補
正処理例２ではノイズの除去処理において、位相エンコ
ード方向に位相差データの平均をとる方法を採用する。
詳しく説明すると、あるラインにおいて、このラインと
高域側のｇ−１本のラインからなるブロックを考え、こ
のブロック内の位相差データから、位相エンコード方向
に位相差データの平均をとることによって、位相差デー
タを平滑化する。そして、この方法で平滑化処理をした
各ラインの位相差データに、最小二乗法を適用すること
によって、数（２）のｐ(ｋy)，ｑ(ｋy)を推定するが、
動き発生ラインと低域側のｇ本のラインにおけるｐ
(ｋy)またはｑ(ｋy)の値は、単調に増加または減少す
る。これは、上で述べたブロックが、動き発生ラインを
低域側から高域側に横切るに従って、ブロック内の動き
発生ライン前部と後部の平均に対する寄与が単調に変化
するので、動き発生ラインにおける、各ラインの位相差
データの傾きまたはオフセットの突発的な変化は、平滑
化した後の各ラインごとの位相差データの傾きまたはオ
フセットの単調な変化として現れる。そして、この単調
な変化は、ブロックの幅ｇ回連続して起こる。
【００８３】従って、動き発生ラインから低域側のｇ−
１本のラインにおける上記の方法によるｐ(ｋy)，ｑ(ｋ
y)の推定値は正しくないので、動き発生ラインから低域
側のｇ本目のラインにおけるｐ(ｋy)，ｑ(ｋy)の値を埋
め込む。また、２つの動き発生ラインどうしがｇライン
以上離れていない場合、その２つの発生ラインの間のｐ
(ｋy)，ｑ(ｋy)の値は、以上の方法では推定できないの
で、高域の動き補正処理例１の方法により、各ラインご
とｐ(ｋy)，ｑ(ｋy)の値を求めることにする。最後に、
低域中央と高域の境界付近の−Ｎ−ｇ〜−Ｎ−１ライン
とＮ＋１〜Ｎ＋ｇラインにおいて、動きが発生している
場合には、その発生ラインを特定できないので、上記の
処理ができない。そこで、低域の動き補正ステップ１０
２において、低域中央３０２に上下ｇラインずつ加えた
領域（−Ｎ−ｇ≦ｋy≦Ｎ＋ｇ）に動き補正を行うこと
によって、これを回避する。
【００８４】以下、高域の動き補正処理例２を図１１に
より説明する。
【００８５】〔ステップ１１０１〜１１０２〕高域の動
き補正処理例１のステップ９０１〜９０２と同じ。
【００８６】〔ステップ１１０３〕高域３０１のライン
ｋy((Ｎy／２−１)−ｇ＋１≧ｋy≧Ｎ＋１）において、 ＥＰＵ(ｋx，ｋy)＝ＰＤ(ｋx，ｋy)＋ＰＤ(ｋx，ｋy＋１)＋ …＋ＰＤ(ｋx，ｋy＋ｇ−１)／ｇ とすることにより、位相差データに平滑化を行う。
【００８７】同様にして、高域３０３のラインｋy(−Ｎ
y／２＋ｇ−１≦ｋy≦−Ｎ−１）において、 ＥＰＤ(ｋx，ｋy)＝ＰＵ(ｋx，ｋy)＋ＰＵ(ｋx，ｋy−１)＋ …＋ＰＵ(ｋx，ｋy−ｇ＋１)／ｇ とすることにより、位相差データに平滑化を行う。
【００８８】〔ステップ１１０４〕高域の動き補正処理
例１の最小二乗法適用区間決定ステップに行き、高域30
1のラインｋy((Ｎy／２−１)−ｇ＋１≧ｋy≧Ｎ＋１)に
おけるＥＰＵ(ｋx，ｋy)の最小二乗法適用区間〔ｓｕ
(ｋy)，ｅｕ(ｋy)〕を求める。
【００８９】同様にして、高域３０３のラインｋy(−Ｎ
y／２＋ｇ−１≦ｋy≦−Ｎ−１）におけるＥＰＤ(ｋx，
ｋy)の最小二乗法適用区間〔ｓｄ(ｋy)，ｅｄ(ｋy)〕を
求める。
【００９０】〔ステップ１１０５〕高域３０１のライン
ｋy((Ｎy／２−１)−ｇ＋１≧ｋy≧Ｎ＋１）において、
区間〔ｓｕ(ｋy)，ｅｕ(ｋy)〕でＥＰＵ(ｋx，ｋy)に最
小二乗法を適用し、数(２)のｐ(ｋy)，ｒ(ｋy)＝ｑ
(ｋy)ｋy((Ｎy／２−１)−ｇ＋１≧ｋy≧Ｎ＋１)を推定
する。
【００９１】同様にして、高域３０３の各ラインｋy(−
Ｎy／２＋ｇ−１≦ｋy≦−Ｎ−１)において、区間〔ｓ
ｄ(ｋy)，ｅｄ(ｋy)〕でＥＰＤ(ｋx，ｋy)に最小二乗法
を適用し、数(２）のｐ(ｋy)，ｒ(ｋy)＝ｑ(ｋy)ｋy(−
Ｎy／２＋ｇ−１≦ｋy≦−Ｎ−１）を推定する。
【００９２】〔ステップ１１０６〕高域３０１の部分領
域((Ｎy／２−１)−ｇ＋１≧ｋy≧Ｎ＋１）において、
ｐ(ｋy)，ｑ(ｋy)が位相エンコード正の方向にｇ回連続
して増加もしくは減少する開始ラインをＮ＋１≦ｊ1＜
ｊ2＜…＜ｊn＜Ｎyとする、同様に高域３０３の部分領
域（−Ｎy／２＋ｇ−１≦ｋy≦−Ｎ−１）においてｐ
(ｋy)，ｑ(ｋy)が位相エンコード負の方向にｇ回連続し
て増加もしくは減少する開始ラインを−Ｎy＜ｊ1′＜ｊ
2′＜…＜ｊn′≦−Ｎ−１とする。
【００９３】〔ステップ１１０７〕各区間〔ｊk，ｊk+1
−１〕(ｋ＝１〜ｎ)(ここで、ｊn+1＝Ｎyとする)におい
て、もしｊk+1−ｊk≧ｇならば、〔ｊk，ｊk+1−ｇ〕に
おいて、ステップ１１０５で推定したｐ(ｋy)，ｑ(ｋy)
により、〔ｊk+1−ｇ＋１，ｊk+1−１〕においては、ｐ
(ｊk+1−ｇ)，ｑ(ｊk+1−ｇ)により、Ｓ(ｋx，ｋy)に位
相補正処理を行ったデータＳＵ(ｋx，ｋy)を得る。
【００９４】もしｊk+1−ｊk＜ｇならば、〔ｊk，ｊk+1
−１〕においてステップ９０５を実行することによって
得られるｐ(ｋy)，ｑ(ｋy)により、Ｓ(ｋx，ｋy)に位相
補正処理を行ったデータＳＵ(ｋx，ｋy)を得る。
【００９５】同様にして、各区間〔ｊk′−１，
ｊk+1′〕(ｋ＝０〜ｎ−１）(ここで、ｊ0＝Ｎyとす
る）において、ｊk+1′−ｊk′≧ｇならば、〔ｊk′＋
ｇ，ｊk+1′〕においては、ステップ１１０５で推定し
たｐ(ｋy），ｑ(ｋy）により、〔ｊk′＋ｇ−１，
ｊk′〕においては、ｐ(ｊk+1′＋ｇ)，ｑ(ｊk+1′＋
ｇ）により、Ｓ(ｋx，ｋy）に位相補正処理を行ったデ
ータＳＤ(ｋx，ｋy)を得る。
【００９６】もしｊk+1′−ｊk′＜ｇならば、
〔ｊk′，ｊk+1′−１〕においてステップ９０５を実行
することにより得られるｐ(ｋy)，ｑ(ｋy)により、Ｓ
(ｋx，ｋy)に位相補正処理を行ったデータＳＤ(ｋx，ｋ
y)を得る。
【００９７】＜実施例２＞実施例２では、高速スピンエ
コー撮影法の場合を考える。図１６（ａ）は図２のシス
テムにおいて、５１２×２５６マトリックスの計測空間
を、図４のように５つのブロックに分けたときに、５エ
コーのマルチエコーシーケンスを利用して、各ブロック
にそれぞれ１つのエコーを割当てることによって各ブロ
ックを計測する高速スピンエコーシーケンスの１例であ
る。この場合、従来のスピンエコー撮影法の約１／５の
撮影時間になるので、撮影中に検査対象物が動くことは
少なくなり、検査対象物の突発的な動きが発生したとし
ても、計測の前半か後半のどちらか一方のみで起こると
いう状況は、いっそう現実的になる。
【００９８】実施例１において、高域で発生した動きを
検出するために、ハーフスキャンイメージング法を利用
したが、そのために計測データの低域中央部を動き補正
する必要があったが、この実施例では、ブロック数より
も１つ多いエコー数を持ったマルチエコーシーケンスを
利用し、低域中央部に相当するブロック４０３に２つの
マルチエコーを割当て計測することによって、計測の前
半，後半どちらで動きが発生しても、動きの影響のない
低域中央部のデータを得ることができるようにした高速
スピンエコー撮影法について説明する。その後、この高
速スピンエコー撮影法によって得られたデータからの動
き補正方法について説明する。
【００９９】図１６（ｂ）は、６エコーのマルチエコー
シーケンスを利用した高速スピンエコーシーケンスであ
り、第１エコーと第２エコーをそれぞれブロック４０１
とブロック４０２に割当て、第５エコーと第６エコーを
それぞれブロック４０４とブロック４０５に割当てる。
また、第３エコーと第４エコーをブロック４０３に割当
て、第３エコーでは、ブロック４０３を上から下に、第
４エコーでは、ブロック４０３を下から上に計測する。
このようにブロック４０３を計測しておけば、動きが計
測の前半か後半のどちらか一方のみで起こる場合は、第
３エコーにより計測されたブロック４０３の上半分と、
第４エコーにより計測されたブロック４０３の下半分、
または、第３エコーにより計測されたブロック４０３の
下半分と、第４エコーにより計測されたブロック４０３
の上半分の、どちらか一方の組合せから得られるブロッ
ク４０３のデータは、動きの影響を含んでいない。
【０１００】以下、この高速スピンエコー撮影法により
計測されたデータにおける動き検出方法について説明す
る。動きの検出は、ブロック４０２と４０４の計測デー
タの複素共役点対称性を利用して行う。１つのブロック
で動きが検出できれば、他のブロックにおいても、その
検出された動きから動き補正処理を行えば良い。第ｉエ
コー（ｉ＝１，２，…，６）によって計測された計測デ
ータをＳi(ｋx，ｋy)（ｉ＝１，２，…，６）とする
と、上記のような高速スピンエコー撮影法により計測さ
れたので、次のような対応関係がある。
【０１０１】 ブロック４０１の計測データ Ｓ1(ｋx,ｋy) ７７≦ｋy≦１２７ ブロック４０２の計測データ Ｓ2(ｋx,ｋy) ２６≦ｋy≦７６ ブロック４０３の計測データ Ｓ3(ｋx,ｋy)，Ｓ4(ｋx,ｋy) −２５≦ｋy≦２５ ブロック４０４の計測データ Ｓ5(ｋx,ｋy) −７６≦ｋy≦−２６ ブロック４０５の計測データ Ｓ6(ｋx,ｋy) −１２７≦ｋy≦−７８ また、ブロック４０３以外の各ブロックは、割当てられ
たエコーにより上から下に計測される。以上の準備のも
と、図１２により動き検出方法について説明する。
【０１０２】〔ステップ１２０１〕ブロック４０３の上
半分のデータＳ3(ｋx，ｋy）(０≦ｋy≦２５）と下半分
のデータＳ4(ｋx，ｋy)（−１≦ｋy≦−２５）による中
央低域の計測データと、ブロック４０４の上半分のデー
タＳ5(ｋx，ｋy)（−２６≦ｋy≦−５１）から、ハーフ
スキャンイメージング法により、ブロック４０２の下半
分の推定データＩＤ2(ｋx，ｋy)（２６≦ｋy≦５１）を
計算する。
【０１０３】同様にして、ブロック４０３の上半分のデ
ータＳ4(ｋx，ｋy)(０≦ｋy≦２５)と下半分のデータＳ
3(ｋx，ｋy)(−１≦ｋy≦−２５)による中央低域の計測
データと、ブロック４０４の下半分のデータＳ5(ｋx，
ｋy)(−７６≦ｋy≦−５１)から、ハーフスキャンイメ
ージング法により、ブロック４０２の上半分の推定デー
タＩＵ2(ｋx，ｋy）(５１≦ｋy≦７６）を計算する。
【０１０４】〔ステップ１２０２〕ブロック４０２の上
半分のデータＳ2(ｋx，ｋy)（５１≦ｋy≦７６）と推定
データＩＵ2(ｋx，ｋy)との位相差をとる。
【０１０５】 ＰＵ2(ｋx，ｋy)＝arctan(Ｓ2(ｋx，ｋy)＊conj(ＩＵ2(ｋx，ｋy))） （５１≦ｋy≦７６） 同様にして、ブロック４０２の下半分のデータＳ
2(ｋx，ｋy)（２６≦ｋy≦５１）と推定データＩＤ2(ｋ
x，ｋy)との位相差をとる。
【０１０６】 ＰＤ2(ｋx，ｋy)＝arctan(Ｓ2(ｋx，ｋy)＊conj((ＩＤ2(ｋx，ｋy))） （２６≦ｋy≦５１） ここで、conj（ ）は複素共役を表わす。
【０１０７】〔ステップ１２０３〕各ラインｋy(２６≦
ｋy≦７６）における位相差データＰＵ2(ｋx，ｋy）と
ＰＤ2(ｋx，ｋy）に対して実施例１の高域の動き補正処
理例１のステップ９０３の方法で平滑化を行い、ＭＰＵ
2(ｋx，ｋy)(５１≦ｋy≦７６）とＭＰＤ2(ｋx，ｋy)
（２６≦ｋy≦５１）を計算する。
【０１０８】〔ステップ１２０４〕実施例１の高域の動
き補正処理例１の最小二乗法適用区間決定ステップに行
き、各ｋy(５１≦ｋy≦７６）において、ＭＰＵ2(ｋx，
ｋy)の最小二乗法適用区間〔ｓｕ(ｋy)，ｅｕ(ｋy)〕を
求める。
【０１０９】同様にして、各ｋy(２６≦ｋy≦５１）に
おいて、ＭＰＤ2(ｋx，ｋy)の最小二乗法適用区間〔ｓ
ｄ(ｋy)，ｅｄ(ｋy)〕を求める。
【０１１０】〔ステップ１２０５〕各ｋy(５１≦ｋy≦
７６）において、区間〔ｓｕ(ｋy)０ｅｕ(ｋy)〕でＭＰ
Ｕ2（ｋx，ｋy)に最小二乗法を適用し、数(２)のｐ2(ｋ
y)，ｒ2(ｋy)＝ｑ2(ｋy)ｋy（５１≦ｋy≦７６）を推定
する。
【０１１１】同様にして、各ｋy(２６≦ｋy≦５１)にお
いて、区間〔ｓｄ(ｋy),ｅd(ｋy)〕でＭＰＤ2(ｋx，
ｋy)に最小二乗法を適用し、数（２）のｐ2(ｋy)，ｒ
2(ｋy)＝ｑ2(ｋy)ｋy（２６≦ｋy≦５１）を推定する。
【０１１２】〔ステップ１２０６〕ステップ１２０５に
より推定したｐ2(ｋy)，ｑ2(ｋy)(５１≦ｋy≦７６）に
もとづき、各ブロックに対する上半分の計測データに ブロック４０１の上半分 Ｓ1(ｋx,ｋy) １０２≦ｋy≦１２７ ｐ(ｋy)＝ｐ
2(ｋy−５１)，ｑ(ｋy)＝ｑ2(ｋy−５１) ブロック４０２の上半分 Ｓ2(ｋx,ｋy) ５１≦ｋy≦７６ ｐ(ｋy)＝ｐ
2(ｋy)，ｑ(ｋy)＝ｑ2(ｋy) ブロック４０３の上半分 Ｓ3(ｋx,ｋy) ０≦ｋy≦２５ ｐ(ｋy)＝ｐ
2(ｋy＋５１)，ｑ(ｋy)＝ｑ2(ｋy＋５１) ブロック４０４の上半分 Ｓ5(ｋx,ｋy) −５１≦ｋy≦−２６ ｐ(ｋy)＝ｐ
2(ｋy＋２×５１),ｑ(ｋy)＝ｑ2(ｋy＋２×５１) ブロック４０５の上半分 Ｓ6(ｋx,ｋy) −１０２≦ｋy≦−７７ ｐ(ｋy)＝ｐ
2(ｋy＋３×５１）,ｑ(ｋy)＝ｑ2(ｋy＋３×５１) とすることによって、位相補正処理を行った後、これら
のデータを集めることによって得られる計測空間全体の
データＳＵ(ｋx，ｋy)を得る。
【０１１３】同様にして、ステップ１２０５により推定
したｐ2(ｋy)，ｑ2(ｋy)（２６≦ｋy≦５１)にもとづ
き、各ブロックに対する下半分の計測データに ブロック４０１の下半分 Ｓ1(ｋx,ｋy) ７７≦ｋy≦１０２ ｐ(ｋy)＝ｐ
2(ｋy−５１)，ｑ(ｋy)＝ｑ2(ｋy−５１) ブロック４０２の下半分 Ｓ2(ｋx,ｋy) ２６≦ｋy≦５１ ｐ(ｋy)＝ｐ
2(ｋy)，ｑ(ｋy)＝ｑ2(ｋy) ブロック４０３の下半分 Ｓ3(ｋx,ｋy) −２５≦ｋy≦０ ｐ(ｋy)＝ｐ
2(ｋy＋５１），ｑ(ｋy)＝ｑ2(ｋy＋５１） ブロック４０４の下半分 Ｓ5(ｋx,ｋy) −７６≦ｋy≦−５１ ｐ(ｋy)＝ｐ
2(ｋy＋２×５１)，ｑ（ｋy)＝ｑ2(ｋy＋２×５１） ブロック４０５の下半分 Ｓ6(ｋx,ｋy) −１２７≦ｋy≦−１０２ ｐ(ｋy)＝ｐ
2(ｋy＋３×５１),ｑ(ｋy)＝ｑ2(ｋy＋３×５１) とすることによって、位相補正処理を行った後、これら
のデータを集めることによって得られる計測空間全体の
データＳＤ(ｋx，ｋy)を求める。
【０１１４】〔ステップ１２０７〕位相補正データＳＵ
(ｋx，ｋy)，ＳＤ(ｋx，ｋy)をフーリエ変換して画像再
生する。
【０１１５】
【発明の効果】本発明は、画像再生に必要な全データを
計測するまでの、前半もしくは後半だけで突発的に起こ
るスライス面内における平行移動性の動きを、特別なシ
ーケンスを用いなくても、計測データ自身から動きを検
出し、その検出した動きから計測データに動き補正する
ことにより、Ｓ／Ｎを劣化させることなく平行移動性の
モーションアーチファクトを抑制することができる。
【図面の簡単な説明】
【図１】本発明の実施例１における動き補正処理の処理
フロー図。
【図２】本発明を実施するシステム構成の一例を示す
図。
【図３】計測空間を低域中央と高域に分けた時の各領域
を示す図。
【図４】６エコーの高速スピンエコー撮影における各エ
コーの計測空間の割当てを示した図。
【図５】計測空間の各ラインにおける位相差データの折
り返しを説明する図。
【図６】本発明の実施例１における第１の低域中央の動
き補正処理例の処理フロー図。
【図７】本発明の実施例１における第２の低域中央の動
き補正処理例の処理フロー図。
【図８】本発明の実施例１における第３の低域中央の動
き補正処理例の処理フロー図。
【図９】本発明の実施例１における第１の高域の動き補
正処理例の処理フロー図。
【図１０】本発明の実施例１における、計測空間の各ラ
インの位相差データに最小二乗法を適用するときの適用
区間を決定する処理フロー図。
【図１１】本発明の実施例１における第２の高域の動き
補正処理例の処理フロー図。
【図１２】本発明の実施例２における動き補正処理例の
処理フロー図。
【図１３】本発明の実施例１において動き補正処理を行
う計測データを撮影するためのシーケンスの一例を示す
図。
【図１４】本発明の実施例１において動き補正処理を行
う計測データの低域中央を撮影するためのシーケンスの
一例を示す図。
【図１５】本発明の実施例１において動き補正処理を行
う計測データを撮影するためのシーケンスの一例を示す
図。
【図１６】高速スピンエコー撮影法における、従来シー
ケンスと本発明シーケンスを示す図。
【符号の説明】
２００…被検体、２０１…静磁場発生系、２０２…送信
系、２０３…傾斜磁場発生系、２０４…受信系、２０５
…処理装置、２０６…ＣＲＴ表示装置、２０７…パルス
シーケンスファイル、２０８…シーケンス制御部、３０
４…ｋy 軸、３０５…ｋx 軸、５０１…折り返し点、１
３０１…ＲＦパルス、１３０２…ｚ方向傾斜磁場、１３
０３…ｙ方向傾斜磁場、１３０４…ｘ方向傾斜磁場、１
３０５…Signal、１３０６…位相エンコードパルス、１
３０７…計測信号、１４０１…ナビゲーションエコー。
フロントページの続き (72)発明者 小泉 英明 茨城県勝田市市毛882番地 株式会社日立 製作所計測器事業部内

Claims (9)

    【特許請求の範囲】
  1. 【請求項１】静磁場，傾斜磁場，高周波磁場の発生装置
    と検査対象物からの核磁気共鳴信号を取り出す検出装置
    と、画像再生構成を含む各種演算を行う処理装置を有す
    る磁気共鳴イメージング装置において、異なる２種類以
    上の動き検出手段を用意しておき、計測データを２個以
    上のブロックに分け、各ブロックは、何れか１種類の動
    き検出手段によって動きを検出し、計測データの動き補
    正を行うことを特徴とする動き補正方法。
  2. 【請求項２】請求項１の動き補正方法において、計測デ
    ータを中央部の低域データとそれ以外の高域データに分
    け、高域データは計測データの位相情報にもとづいて動
    き検出することを特徴とする動き補正方法。
  3. 【請求項３】請求項２における低域中央部の動き補正処
    理は、各ラインの計測データを一次元フーリエ変換して
    得られるデータの位置情報にもとづき、ライン方向の動
    きを補正する処理を含むことを特徴とする動き補正方
    法。
  4. 【請求項４】請求項２における低域中央部の動き補正処
    理は、各ラインの位相成分の平均をとり、ライン毎の位
    相の平均値のずれ量により位相オフセットを検出して、
    ラインと垂直の方向（位相エンコード方向）の動きを補
    正する処理を含むことを特徴とする動き補正方法。
  5. 【請求項５】静磁場，傾斜磁場，高周波磁場の発生装置
    と検査対象物からの磁気共鳴信号を取り出す検出装置
    と、画像再生構成を含む各種演算を行う処理装置を有す
    る磁気共鳴イメージング装置において、計測データを上
    半分と下半分に分け、上半分の計測データおよび下半分
    の計測データの動きによる位相誤差を推定し、その位相
    誤差にもとづき位相補正処理を行うことを特徴とする動
    き補正方法。
  6. 【請求項６】請求項５の動き補正方法において位相誤差
    の推定は、片側半分のデータを、もう一方の側の半分の
    計測データから複素共役点対称性を利用して推定し、こ
    の推定したデータと推定した側の計測データとの位相差
    を計算し、動きによる位相誤差を位相差データに最小二
    乗法を適用することによって推定することを特徴とする
    動き補正方法。
  7. 【請求項７】請求項６において、位相差データに最小二
    乗法を適用する処理は、注目しているラインに、そのラ
    インに隣接する数ラインを加えたブロックを考え、この
    ブロック内でラインに垂直な方向（位相エンコード方
    向）に位相データの平均をとることによって位相差デー
    タに平滑化を行う処理を含むことを特徴とする動き補正
    方法。
  8. 【請求項８】請求項６における、位相差データに最小二
    乗法を適用する処理は、各ラインの中央低周波区間にお
    いて、位相差の値が−πとπをとる点を含まない区間を
    求めて、この区間で最小二乗法を適用する処理を含むこ
    とを特徴とする動き補正方法。
  9. 【請求項９】計測空間を複数のブロックに分け、各ブロ
    ック毎にマルチエコーシーケンスを用いた異なるエコー
    タイム（ＴＥ）で計測したスピンエコー信号を割り当て
    て、画像再生に必要な計測データを得る高速スピンエコ
    ー撮影法において、エコー数をブロック数よりも１つ多
    くして、位相エンコードがゼロのラインを含む計測空間
    中央のブロックを２つのマルチエコーで各々のエコーが
    該ブロックを互いに反対方向に計測することを特徴とす
    る高速スピンエコー撮影法。
JP4206434A 1992-08-03 1992-08-03 Ｍｒイメージングにおける動き補正方法 Pending JPH0647021A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP4206434A JPH0647021A (ja) 1992-08-03 1992-08-03 Ｍｒイメージングにおける動き補正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP4206434A JPH0647021A (ja) 1992-08-03 1992-08-03 Ｍｒイメージングにおける動き補正方法

Publications (1)

Publication Number Publication Date
JPH0647021A true JPH0647021A (ja) 1994-02-22

Family

ID=16523315

Family Applications (1)

Application Number Title Priority Date Filing Date
JP4206434A Pending JPH0647021A (ja) 1992-08-03 1992-08-03 Ｍｒイメージングにおける動き補正方法

Country Status (1)

Country Link
JP (1) JPH0647021A (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011072652A (ja) * 2009-09-30 2011-04-14 Ge Medical Systems Global Technology Co Llc 磁気共鳴イメージング装置
US8126237B2 (en) 2004-11-12 2012-02-28 Kabushiki Kaisha Toshiba Magnetic resonance imaging and correcting device

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8126237B2 (en) 2004-11-12 2012-02-28 Kabushiki Kaisha Toshiba Magnetic resonance imaging and correcting device
US8526703B2 (en) 2004-11-12 2013-09-03 Kabushiki Kaisha Toshiba Making different corrections for different areas of magnetic resonance image data
JP2011072652A (ja) * 2009-09-30 2011-04-14 Ge Medical Systems Global Technology Co Llc 磁気共鳴イメージング装置
US8564290B2 (en) 2009-09-30 2013-10-22 Ge Medical Systems Global Technology Company, Llc Magnetic resonance imaging apparatus and method

Similar Documents

Publication Publication Date Title
Block et al. Spiral imaging: a critical appraisal
Manke et al. Respiratory motion in coronary magnetic resonance angiography: a comparison of different motion models
US9629612B2 (en) Biomedical image reconstruction method and apparatus
US12523726B2 (en) Joint k-space and image-space reconstruction imaging method and device
Schäffter et al. Motion compensated projection reconstruction
EP1644751B1 (en) BACKGROUND SUPPRESSION METHOD FOR TEMPORARILY RELATED MAGNETIC RESONANCE ANGIOGRAPHY
US6518759B2 (en) Motion correction of magnetic resonance images
US20140210469A1 (en) Nonrigid motion correction in 3d using autofocusing with localized linear translations
Hedley et al. Motion artifact suppression: a review of post-processing techniques
Lai et al. A respiratory self‐gating technique with 3D‐translation compensation for free‐breathing whole‐heart coronary MRA
US10832452B2 (en) Method and apparatus for generating magnetic resonance image
US9366741B2 (en) Medical image imaging method, medical diagnostic apparatus using the same, and recording medium therefor
JPH10137215A (ja) 画像の分解能を向上させる方法、投影画像の一部を拡大させる方法及び画像の分解能を向上させる装置
CN113298902A (zh) 卷褶视野磁共振图像的重建方法、计算机设备及存储介质
CN100468459C (zh) 基于边缘检测的局域化抗混叠方法
Marchetto et al. Robust retrospective motion correction of head motion using navigator‐based and markerless motion tracking techniques
WO2002010790A1 (en) Magnetic resonance method for forming a fast dynamic imaging
JPS63183046A (ja) 磁気共鳴映像装置
US20090091321A1 (en) Self navigating cartesian trajectory
JP2641486B2 (ja) Ｎｍｒイメージング装置
Chiou et al. A simple simultaneous geometric and intensity correction method for echo-planar imaging by EPI-based phase modulation
US11333734B2 (en) Quality and signal-to-noise ratio of free-breathing quantitative measurement of magnetic resonance imaging parameters and related biomarkers
McLeish et al. Free‐breathing radial acquisitions of the heart
JP4975614B2 (ja) 磁気共鳴イメージング装置及び方法
JP2595006B2 (ja) Ｍｒイメージング方法