JPH10155767A - 動きアーティファクトを除去する磁気共鳴方法及び装置 - Google Patents

動きアーティファクトを除去する磁気共鳴方法及び装置

Info

Publication number
JPH10155767A
JPH10155767A JP9312090A JP31209097A JPH10155767A JP H10155767 A JPH10155767 A JP H10155767A JP 9312090 A JP9312090 A JP 9312090A JP 31209097 A JP31209097 A JP 31209097A JP H10155767 A JPH10155767 A JP H10155767A
Authority
JP
Japan
Prior art keywords
magnetic resonance
image
resolution
data set
raster
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
JP9312090A
Other languages
English (en)
Inventor
Volker Dr Rasche
ラッシェ フォルカー
Roland Proksa
プロスカ ロラント
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.)
Koninklijke Philips NV
Original Assignee
Philips Electronics NV
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 Philips Electronics NV filed Critical Philips Electronics NV
Publication of JPH10155767A publication Critical patent/JPH10155767A/ja
Pending legal-status Critical Current

Links

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/56509—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling
    • 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/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/4824—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory
    • 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/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56308—Characterization of motion or flow; Dynamic imaging

Landscapes

  • Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General 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)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

(57)【要約】 【課題】 本発明は検査対象に関する複数のＭＲ（磁気
共鳴）データセットが時間的に連続して取得されるＭＲ
方法において、高解像度ＭＲ画像再生用のＭＲデータセ
ットを取得するため要する時間間隔中の検査対象の動き
に起因して高解像度ＭＲ画像に生じる妨害性の動きアー
ティファクトを除去するＭＲ方法の提供を目的とする。 【解決手段】 本発明によれば、低解像度ＭＲ画像に対
し複数のＭＲデータセットが短縮された測定期間内に取
得され、取得されたＭＲデータセットから高解像度ＭＲ
画像を再生する際に考慮されるべき画像変換パラメータ
が低解像度画像の比較から決定される。

Description

【発明の詳細な説明】
【０００１】
【発明の属する技術分野】本発明は、検査対象に関する
複数の磁気共鳴（ＭＲ）データセットが時間的に連続し
て取得される磁気共鳴方法において動きアーティファク
トを除去する磁気共鳴（ＭＲ）方法に関する。また、本
発明は、磁気共鳴信号から磁気共鳴画像を形成する手段
と、上記手段を制御するプログラマブル制御ユニット
と、投影再生法を用いて磁気共鳴画像を再生する再生ユ
ニットとからなり上記方法を実施する磁気共鳴装置に関
する。
【０００２】
【従来の技術】上記の種類の磁気共鳴方法及び磁気共鳴
装置は、米国特許明細書第5,502,385号により公知であ
る。上記明細書に記載された方法によれば、連続的に変
化させられる傾斜磁界を備えた磁気共鳴シーケンスが検
査対象に作用する。磁気共鳴画像はこのようにして発生
された磁気共鳴信号から再生され、磁気共鳴画像の再生
のため考慮された最先の磁気共鳴信号は、次の磁気共鳴
画像の再生のため新たに発生された夫々の磁気共鳴信号
により置き換えられる。検査ゾーンで生ずる動的な過程
は、再生のため要求される傾斜磁界の向きの範囲が新た
に発生された夫々の磁気共鳴信号の間で均一に分布され
るように、あるシーケンスから別のシーケンスまで磁気
共鳴信号の発生中に印加される傾斜磁界の向きを変化さ
せることにより僅かなアーティファクトを伴って再現さ
れる。
【０００３】一般的に言うと、高解像度磁気共鳴画像を
形成するため運動中の対象の磁気共鳴データセットを取
得するには、対象の動きに対し無視できる程には短くな
い時間の間隔が必要である。データ取得中の検査対象の
動きのため、アーティファクトが従来の方法を用いて再
生された高解像度磁気共鳴画像に発生する傾向がある。
しかし、所定の磁気共鳴画像の解像度が指定されるなら
ば、測定時間は一般的に非常に短い程度まで短縮され得
る。
【０００４】
【発明が解決しようとする課題】本発明の目的は、検査
ゾーン内の動的な過程の高解像度磁気共鳴画像に生じる
アーティファクトが低減された上記の種類の磁気共鳴方
法及び磁気共鳴装置を提供することである。
【０００５】
【課題を解決するための手段】上記の目的は、各低解像
度磁気共鳴画像が磁気共鳴データセットから再生され、
画像変換パラメータが低解像度磁気共鳴画像の比較によ
り決定され、高解像度磁気共鳴画像が上記画像変換パラ
メータに依存して磁気共鳴データセットから再生される
本発明の磁気共鳴方法により達成される。
【０００６】本発明によれば、高解像度磁気共鳴画像を
再生可能な単一の磁気共鳴データセットを取得するため
必要とされる時間の間隔は複数の（短期）時間の間隔に
細分され、各時間の間隔において夫々の磁気共鳴データ
セットが取得される。低解像度磁気共鳴画像は各磁気共
鳴データセットから再生可能である。この点に関し、用
語「低解像度」は、再生のため必要とされる磁気共鳴デ
ータセットが、高解像度磁気共鳴画像の再生のための磁
気共鳴データセットを取得するために使用される磁気共
鳴シーケンスの個数よりも少ない個数の磁気共鳴シーケ
ンスを用いて取得されることを意味する。即ち、所謂ｋ
空間（空間周波数空間）において、ｋ空間は高解像度磁
気共鳴画像の再生よりも少ない経路に沿って走査される
ので、ｋ空間の角度方向の走査点の密度は小さい。
【０００７】画像変換パラメータは、個別の画素、画像
セグメント、又は、完全な画像対象（例えば、骨）の位
置が連続的な磁気共鳴画像で変化する様子を記述する。
各画素、画像セグメント又は画像対象には、好ましく
は、二つの時間的に連続した低解像度磁気共鳴画像の比
較により取得され、動きベクトルのような相互に比較さ
れた第１の磁気共鳴画像から第２の磁気共鳴画像までの
画素、画像セグメント又は画像対象の位置の変化を表わ
す夫々の画像変換パラメータが割り当てられる。二つの
低解像度磁気共鳴画像から取得された画像変換パラメー
タは、例えば、像変換マトリックス（所謂、動きベクト
ル場）の形式に結合される。
【０００８】高解像度磁気共鳴画像は、上記画像変換パ
ラメータと、低解像度磁気共鳴画像が再生された磁気共
鳴データセットとから再生可能であり、この高解像度磁
気共鳴画像は、磁気共鳴データセットの取得中に検査ゾ
ーン内の対象が移動、回転又は変形するが、新しい対象
が検査ゾーン内に入らないならば、全く動きアーティフ
ァクトを示さず、若しくは、僅かしか動きアーティファ
クト示さない。
【０００９】本発明による方法の一つの変形は、基準画
像としての役割を果たす低解像度磁気共鳴画像内の画素
の位置が決定され、他の低解像度磁気共鳴画像内の画素
の位置が判定され、基準画像内の画素の位置と、他の低
解像度磁気共鳴画像内の画素の位置との比較から画像変
換パラメータが決定されることを特徴とする。時間的に
最初に取得された磁気共鳴データセットから再生された
低解像度磁気共鳴画像は、好ましくは基準画像として使
用され、画像変換パラメータは後続の磁気共鳴画像内の
画素の移動を記述する。しかし、別の磁気共鳴画像を基
準画像として選択してもよい。
【００１０】本発明に基づく他の変形例は、基準画像内
において画素の位置は、高解像度磁気共鳴画像が再生さ
れたラスタのラスタ位置に対応した場所内で決定される
ことを特徴とする。高解像度磁気共鳴画像の再生のため
画像変換パラメータが考慮されるので、高解像度磁気共
鳴画像の再生のため使用されるラスタが画像変換パラメ
ータの判定のため使用される点が有利である。このた
め、（低解像度）基準画像は、（好ましくは、ラスタ像
並びに他の低解像度磁気共鳴画像がこのラスタ上で再生
される）ラスタに従って細分され、基準画像内のこのラ
スタのラスタ点の画素の位置は他の低解像度磁気共鳴画
像内で判定され、追跡される。このようにしなければ、
高解像度像が再生されるようになる前に画像変換パラメ
ータを高解像度像の再生用のラスタに適合させる必要が
あるので、高解像度磁気共鳴画像の再生用の計算作業は
付加的に増加する。
【００１１】本発明の特に魅力的な変形例において、高
解像度磁気共鳴画像の画像値は各細分画像の加算によっ
て決定され、但し、第１の細分画像の値は基本ラスタの
ラスタ点での再生により第１の磁気共鳴データセットか
ら決定され、別の細分画像の値は関連した画像変換パラ
メータを用いて基本ラスタから得られたラスタのラスタ
点での再生により他の各磁気共鳴データセットから決定
される。高解像度磁気共鳴画像の画像値は、複数の細分
画像の値によって構成され、細分画像の個数は磁気共鳴
データセットの個数と一致する。最初に、第１の細分画
像は基本ラスタのラスタ点で再生される。次に、このラ
スタ点は関連した画像変換パラメータを用いてシフトさ
れ、このラスタ点の新しい位置が決定される。その場所
で次の細分画像の値が関連した次の磁気共鳴データセッ
トから再生される。この動作は全ての磁気共鳴データセ
ットに対し繰り返される。続いて、決定された全ての細
分画像の値が加算され、高解像度磁気共鳴画像の画像値
が得られる。高解像度磁気共鳴画像の画像値の位置は、
第１の細分画像が再生された基本ラスタのラスタ点の位
置に対応する。最後に、磁気共鳴データセットの取得中
に検査ゾーン内で生じた動きが考慮され、後で高解像度
像の再生中に磁気共鳴データセットから除去される。
【００１２】本発明の他の変形例は、高解像度磁気共鳴
画像の再生のため、磁気共鳴基本画像がラスタに基づい
て磁気共鳴データセットから形成され、各磁気共鳴画像
は、関連した画像変換パラメータを用いて基本ラスタか
ら得られた別のラスタに基づいて残りの磁気共鳴データ
セットから形成され、磁気共鳴基本画像及び磁気共鳴画
像は高解像度磁気共鳴画像を形成するように重ね合わさ
れることを特徴とする。磁気共鳴データセットの取得中
に発生し、基本ラスタ内で考慮された動きは、画像変換
パラメータを用いて、変換されたラスタと関連した磁気
共鳴データセットから磁気共鳴画像が再生される別のラ
スタに変換される。一般的に言うと、変換されたラスタ
は、一般的に例えば、デカルトラスタのような規則性を
もたないので、上記磁気共鳴画像を再生し得ない。再生
の線形性のため、高解像度磁気共鳴画像は磁気共鳴基本
画像及び磁気共鳴画像の重ね合わせにより得られる。
【００１３】本発明の好ましい一変形例において、ｋ空
間は磁気共鳴データセットを取得するため、ｋ空間の零
点を通り放射状に延在する経路に沿って走査される。こ
のため、磁気共鳴シーケンスは異なる位相符号化傾斜を
有する。このようなｋ空間の放射状走査は、磁気共鳴画
像が投影再生法を用いて容易に形成され得る利点があ
る。
【００１４】上記の好ましい一変形例に基づく本発明の
更なる変形例は、磁気共鳴データセットを取得するため
辿られる経路は、別の磁気共鳴データセットを取得する
ため辿られる経路と所定の角度で交差することを特徴と
する。好ましくは、ｋ空間は、高解像度磁気共鳴画像の
再生のため必要とされる磁気共鳴データセットの取得中
に均一に走査され、二つの隣接した経路間の全ての角度
は上記目的のため一致する。
【００１５】また、上記本発明の目的は、磁気共鳴デー
タセットの取得と並行して、検査対象の所定の位置に配
置されたマイクロコイルの位置が決定され、画像変換パ
ラメータが時間的に連続したマイクロコイルの位置デー
タの比較から決定され、高解像度磁気共鳴画像が画像変
換パラメータに依存して磁気共鳴データセットから再生
される上記の種類の磁気共鳴方法により達成される。従
って、本発明の方法によれば、低解像度磁気共鳴画像内
で画素又は画像対象の位置を決定及び追跡する代わり
に、検査対象よりも前方又は検査対象内に配置されたマ
イクロコイルの位置が決定並びに追跡され、各コイルは
磁気共鳴データセットの取得中に磁気共鳴信号を取得す
る。次に、画像変換パラメータが磁気共鳴データセット
から得られる。高解像度磁気共鳴画像は、次に上記の方
法で再生される。
【００１６】本発明の更なる変形例において、一つ以上
の低解像度磁気共鳴画像及び／又は高解像度磁気共鳴画
像が表示される。観察者には時間的及び空間的に異なる
解像度の磁気共鳴画像が提示される。低解像度磁気共鳴
画像は、蛍光透視法応用（例えば、介挿的な応用の場
合）のため望まれるように、空間解像度を犠牲にして高
い時間的な解像度を有し、これに対し、高解像度磁気共
鳴画像は時間的な解像度を犠牲にして高い空間解像度を
実現する。表示される磁気共鳴画像の空間解像度及び／
又は時間的解像度と、数とを観察者に選択させてもよ
い。
【００１７】上記更なる変形例に基づく別のバージョン
は、高解像度磁気共鳴画像が画像変換パラメータを考慮
することなく磁気共鳴データセットから再生されること
を特徴とする。かくして、画像変換パラメータに依存し
て再生される高解像度磁気共鳴画像の代わりに、従来の
方法を用いて磁気共鳴データセットから再生された磁気
共鳴画像が観察者に提供され得る。この磁気共鳴画像は
動きアーティファクトを含む場合があるが、画像変換パ
ラメータの計算を省略することが可能であり、高解像度
磁気共鳴画像の再生に必要とされる計算作業は少なくな
る。また、低解像度磁気共鳴画像を高解像度磁気共鳴画
像のラスタ上で再生することが可能であり、この目的の
ため、低解像度磁気共鳴画像は、例えば、数学的手法を
用いて適切にズームされる。この本発明の別のバージョ
ンにおいて、磁気共鳴データセットはｋ空間の放射状走
査により取得される必要がなく、走査は異なる方法で行
ってもよい。
【００１８】磁気共鳴装置に関する本発明の目的は、検
査対象の複数の磁気共鳴データが時間的に連続して取得
されるように制御ユニットがプログラムされ、再生ユニ
ットが磁気共鳴データセットから各低解像度磁気共鳴画
像を再生し、低解像度磁気共鳴画像の比較から画像変換
パラメータを取得し、画像変換パラメータに依存して磁
気共鳴データセットから高解像度磁気共鳴画像を再生す
ることにより実現される。
【００１９】以下、添付図面を参照して本発明を詳細に
説明する。
【００２０】
【発明の実施の形態】図１は、本発明による方法を実施
する磁気共鳴装置のブロック図である。同図には、例え
ば、４個のコイルからなる均一、静止磁界を発生させる
システム１が示される。システム１の内部で検査される
患者１０が搬送台上部４に収容される。三つの傾斜コイ
ルシステム３、５、７は、所謂傾斜磁界Ｇx 、Ｇy 、Ｇ
z を発生させるため設けられる。上記の各傾斜コイルシ
ステムは、均一磁界の方向に広がり、かつ、上記の三つ
の空間方向の中のいずれかの向きに線形に変化する夫々
の磁界を発生させる。
【００２１】更に、ＲＦパルスを発生する発振器６、若
しくは、検査ゾーン内に発生された磁気共鳴信号を受信
する受信器９のいずれかに、スイッチ８を介して接続さ
れ得るＲＦコイル１１が設けられている。しかし、磁気
共鳴シーケンスの受信用に別個のＲＦ受信器コイル（図
示しない）を使用してもよい。上記の磁気共鳴装置の構
成要素は、特に磁気共鳴シーケンスの時間的な連続性を
制御する制御ユニット１２により制御される。
【００２２】高速フーリエ変換（ＦＦＴ）演算ユニット
１３は測定された磁気共鳴信号から投影Ｐi （α）を計
算し、アーティファクトのない高解像度磁気共鳴画像
は、本発明の方法を用いて再生ユニット１４内で上記投
影から再生される。この画像は次に出力ユニット１５上
に表示され得る。ユニット１３、１４、１５は制御ユニ
ット１２により制御される。
【００２３】以下、図２のフローチャートを参照して本
発明による方法を詳述する。特に、再生ユニット１４に
より実行される段階が、ステップ１４１乃至１４４とし
て示されている。受信器９でディジタル化され、ディジ
タル走査値として受信器出力に現れる各磁気共鳴信号
は、ステップ１３において、逆フーリエ変換を受ける。
これにより得られた投影Ｐi （α）は、次に、ステップ
１４１においてフィルタリングされた背景映写を用いて
磁気共鳴画像Ｉi を再生するため使用される。通常、ｋ
空間の操作中に所定の時間間隔で決定された全ての投影
Ｐi （α）は高解像度磁気共鳴画像を再生するため使用
されるのに対し、本発明によれば、最初に、予め取得さ
れた多数の投影Ｐi （α）だけを使用して複数の低解像
度磁気共鳴画像Ｉi が再生される。各投影Ｐi （α）は
厳密に１回しか使用されない。例えば、１８０個の投影
Ｐi （α）から高解像度磁気共鳴画像を再生するのでは
なく、本発明によれば、最初に、三つの低解像度磁気共
鳴画像Ｉi （ｉ＝１，．．．，３）が６０個の投影Ｐi
（α）毎に再生される。
【００２４】ｋ空間は、低解像度磁気共鳴画像Ｉi に対
する投影Ｐi （α）を得るために均一に走査される。こ
れに関し、本発明によるｋ空間の走査法を示す図３を参
照して説明する。かくして得られた各磁気共鳴シーケン
ス又は各磁気共鳴信号（そこから投影がＦＦＴを用いて
計算される）は、ｋ空間の零点を通過する直線に対応
し、直線の向きは傾斜磁界の傾斜方向と一致する。最初
の時間間隔にｋ空間は実線の経路Ｓ11、Ｓ12、Ｓ13、Ｓ
14方向に走査される。次に、走査は破線の経路Ｓ 21、Ｓ
22、Ｓ23、Ｓ24方向に行われ、その後、一点鎖線の経路
Ｓ31、Ｓ32、Ｓ33、Ｓ34方向に行われる。ある走査のグ
ループの傾斜方向は、他のグループの傾斜方向から常に
一定の角度（α21）の偏りがあることが分かる。即ち、
Ｎ／ｎ個の走査グループのＮ個の走査経路の向きは、ｋ
空間全体に規則的に分散している（Ｎは走査経路の数を
表わし、ｎが１走査グループ当たりの走査経路の数を表
わし、本実施例の場合にＮ＝１２及びｎ＝４であり、走
査グループの数Ｎ／ｎ＝１２／４＝３である）。
【００２５】第１の走査グループの経路Ｓ11乃至Ｓ14に
沿って取得され、一体的に第１の磁気共鳴データセット
を構成する投影Ｐ1 に基づいて、つぎに、第１の低解像
度磁気共鳴画像Ｉ1 が再生され、第２の低解像度共鳴像
Ｉ2 が第２の走査グループ（経路Ｓ21乃至Ｓ24）の投影
Ｐ2 （第２の磁気共鳴データセット）から再生され、第
３の低解像度共鳴像Ｉ3 が第３の走査グループ（経路Ｓ
31乃至Ｓ34）の投影Ｐ 3 （第３の磁気共鳴データセッ
ト）から再生される。上記の再生は、ステップ１４１に
おいてフィルタリングされた背景映写を用いて行われ、
即ち、ｋ空間内の円Ｍの外側の走査点の密度は非常に小
さいのでこの円Ｍの内側に在る投影のデータだけが再生
のため使用される。
【００２６】走査の時間的な連続性は、傾斜方向αが縦
軸、時間ｔが横軸にプロットされた図４に示されてい
る。同図の各点は磁気共鳴信号の取得に対応する。０乃
至αma x の角度範囲は、高解像度磁気共鳴画像の再生の
ため必要とされる傾斜方向の範囲である。走査グループ
の経路は、夫々の時間間隔Ｔg 内で順次に追跡されるこ
とが分かる。第１の時間間隔Ｔg の間に、例えば、第１
の走査グループの経路は投影Ｐ1 を取得すべく追跡され
る。、一般的には、Ｔm ＝Ｎ／ｎ Ｔg 、本実施例の場
合には全体として時間Ｔm ＝３Ｔg が必要とされる。
【００２７】図２のステップ１４２において、かくして
取得された低解像度磁気共鳴像Ｉiは比較され、画像変
換パラメータＴijがその比較から得られる。これに関し
て図５を参照して説明する。図５には、例えば、図３に
従ってｋ空間の走査の投影から取得された時間的に連続
した３個の低解像度磁気共鳴画像Ｉ1 、Ｉ2 、Ｉ3 が示
されている。二つの画像対象Ａ及びＢが識別され、画像
対象Ｂは磁気共鳴データの取得中に移動している。画像
変換パラメータＴ12は、第２の磁気共鳴画像Ｉ2 と、基
準画像Ｉrefとして作用する第１の磁気共鳴画像Ｉ1 と
の比較により取得され、画像変換パラメータＴ23は、第
３の磁気共鳴画像Ｉ3 と第２の磁気共鳴画像Ｉ2 との比
較により取得される。このため、基準画像Ｉref 内で画
素が、選択された位置、例えば、基準画像Ｉref の再生
のため使用されたラスタＲ1 のラスタ点で判定され、次
の磁気共鳴画像Ｉ2 、Ｉ3 内の上記画素の位置が決定さ
れる。この点に関して、例えば、画素ｂi を参照して説
明する。
【００２８】画像対象Ｂの画素ｂi は基準画像Ｉref 内
のラスタＲ1 のラスタ点ｒ31（Ｉ1）に在る。画像対象
Ｂの動きのため、画像Ｉ2 及びＩ3 内の画素ｂ1 の位置
は図５に示される如く変化する。三つの低解像度磁気共
鳴画像Ｉ1 、Ｉ2 、Ｉ3 内の画素ｂi の位置を比較する
ことにより、画像変換パラメータがラスタ点ｒ31に在る
画素ｂi に対し決定される。選択された画素の位置は、
例えば、基準画像内の画素と関係した強度（画像値）を
決定し、別の磁気共鳴画像内の画像値に基づいて画素の
位置を決定することにより、磁気共鳴画像内で決定する
ことが可能である。この目的のための数学的な手法は、
例えば、ＩＥＥＥ、ビデオ技術のための回路及びシステ
ム論文誌、第３巻、第５号（１９９３年）、３６８ペー
ジに提案されている。
【００２９】この演算は基準画像Ｉref 内のラスタＲ1
のラスタ点に在る全ての画素に対し実行され、これによ
り、全ての画像変換パラメータＴ12及びＴ23が得られ
る。次に、図２のステップ１４３において、高解像度磁
気共鳴画像Ｉの画像値が画像変換パラメータＴij及び投
影Ｐi （α）から計算される。高解像度磁気共鳴画像Ｉ
の画像値の計算は図６を参照して詳述する。図６には高
解像度磁気共鳴画像Ｉの空間が座標ｘ及びｙと共に示さ
れている。
【００３０】最初に、画素ｂ1 の第１の低解像度磁気共
鳴画像Ｉ1 における位置ｒ31（Ｉ1）で、第１の細分画
像の値が第１の投影Ｐ1 （α11）から再生される。これ
は、図６において、投影Ｐ1 （α11）から位置Ｒ31（Ｉ
1 ）までの点線によって示される。同じ位置で別の細分
画像の値が、図３に示された第１の走査グループの別の
経路Ｓ12乃至Ｓ14方向に取得された別の投影Ｐ
1 （α1V；ｖ＝２，３，４）から再生され、別の細分画
像の値は第１の細分画像の値に加算される。同図には、
投影Ｐ1 （α12）からの再生だけが点線により示されて
いる。
【００３１】次に、画素ｂ1 の第２の低解像度磁気共鳴
画像Ｉ2 における位置ｒ31（Ｉ2 ）で、第２の細分画像
の値が第２の走査グループの経路Ｓ21乃至Ｓ24方向に取
得された投影Ｐ2 （α2V；ｖ＝１，２，３，４）から再
生され、加算される。同様のことが、画素ｂ1 の第３の
低解像度磁気共鳴画像Ｉ3 における位置ｒ31（Ｉ3 ）
で、第３の走査グループの経路Ｓ31乃至Ｓ34方向に取得
された投影Ｐ3 （α3V；ｖ＝１，２，３，４）に基づい
て行われる。これらは、図６において、投影Ｐ
2（α21）から位置ｒ31（Ｉ2 ）までの破線、並びに、
投影Ｐ3 （α31）から位置ｒ31（Ｉ3 ）までの一点鎖線
によって夫々示されている。
【００３２】次に、決定された全ての細分画像の値が図
２のステップ１４４において加算され、位置ｒ
31（Ｉ1 ）に在る高解像度磁気共鳴画像の画像値が得ら
れる。上記の細分画像の値の再生は、低解像度磁気共鳴
画像の形成の場合のように円Ｍの内部に在る投影のデー
タだけではなく、ｋ空間内の経路に沿って判定された投
影の全てのデータを利用する。
【００３３】上記の方法は、ラスタＲ1 の全てのラスタ
点ｒkl（但し、ｋ＝１，２，．．．，Ｋ；ｌ＝１，
２，．．．，Ｌ；Ｋはｙ方向のラスタ点の数；Ｌはｘ方
向のラスタ点の数）に対し繰り返され、即ち、高解像度
磁気共鳴画像Ｉの全ての画像値は、別々のラスタ点にお
ける再生によって別個の夫々の投影から取得された細分
画像の値から得られる。換言すれば、第１の細分画像の
値は第１の走査グループの投影Ｐ1 から基本ラスタＲ1
に基づいて再生され、画像変換パラメータを用いた変換
により基本ラスタから得られたラスタに基づく別の細分
画像の値は、別の走査グループの投影から再生され、上
記細分画像の値は、次に基本ラスタＲ1 上で加算され、
再生される。
【００３４】図７には種々のラスタＲ1 乃至Ｒ3 が示さ
れる。高解像度磁気共鳴画像Ｉが再構成、再生されるべ
き基本ラスタＲ1 は、第１の低解像度磁気共鳴画像Ｉ1
のラスタＲ1 に対応する。ラスタＲ2 は、ラスタＲ2 の
新しいラスタ点の座標が関連した座標変換パラメータを
用いてラスタＲ1 の各ラスタ点の座標から計算されるこ
とにより、画像変換パラメータＴ12を用いてラスタＲ1
から得られる。同様に、ラスタＲ3 は画像変換パラメー
タＴ23を用いてラスタＲ2 から得られる。例えば、ラス
タＲ1 のラスタ点ｒ31（Ｉ1 ）は、磁気共鳴画像Ｉ2 内
の画素ｂ1 の位置であるラスタＲ2 のラスタ点ｒ31（Ｉ
2 ）にシフトされ、磁気共鳴画像Ｉ3 内の画素ｂ1 の位
置であるラスタＲ3 のラスタ点ｒ31（Ｉ3 ）にシフトさ
れる。上記の如く、高解像度磁気共鳴画像の画像値を形
成するため加算されるべき細分画像の値は、上記のラス
タＲ1 、Ｒ2 、Ｒ3 のラスタ点上に再生される。
【００３５】画像変換パラメータの計算のため、別の低
解像度磁気共鳴画像内で位置が追跡された基準画像内の
画素はラスタ上に存在する必要はない。この目的のた
め、任意の（しかし、既知の）位置で任意の数の画素を
選択してもよい。画像変換パラメータは、画素のサイズ
よりも大きい精度、若しくは、２個の画素の間の距離よ
りも小さい精度で決定し得る。
【００３６】画像変換パラメータの計算のため、数学的
手法を用いて低分解能磁気共鳴画像を急激に拡大（ズー
ム）し、高解像度磁気共鳴画像上に再生してもよい。画
像変換パラメータは、上記の方法を用いて高解像度磁気
共鳴画像のラスタ点で決定され得るので、より多くの画
像変換パラメータが取得され、最終的により良い精度が
得られる。
【００３７】また、第１の時間的に低解像度の磁気共鳴
画像の代わりに、別の磁気共鳴画像を基準画像として選
択し、その基準画像から画像変換パラメータを決定し、
最終的に、高解像度磁気共鳴画像を再生することが可能
である。高解像度磁気共鳴画像は、動きの後半の段階
中、即ち、基準画像に関係した磁気共鳴データセットが
取得された時点の画像対象の位置を表わす。
【００３８】投影再生方法の代わりに、低解像度磁気共
鳴画像及び／又は高解像度磁気共鳴画像を形成するため
離散フーリエ変換を使用してもよい。しかし、不規則な
ラスタ上で画像値を計算する必要があるので、著しく長
い時間を要する。本発明による方法は、新しい画像対象
が画像化動作の途中に検査ゾーン内に入ることはない
が、画像対象の相互の位置が変化しやすい状況で使用す
ることが可能である。この状況は、例えば、動的な関節
の研究の間に生ずる。この場合、解剖学的構造の部位は
動きの間に変位及び／又は変形するが、解剖学的構造自
体は変化しない。
【図面の簡単な説明】
【図１】本発明を実施するため適当な磁気共鳴装置のブ
ロック図である。
【図２】本発明による方法を説明するフローチャートで
ある。
【図３】本発明によるｋ空間の走査を説明する図であ
る。
【図４】図３の走査の時間的な連続性を示す図である。
【図５】複数の連続的な低解像度磁気共鳴画像を示す図
である。
【図６】本発明による高解像度磁気共鳴画像の再生を説
明する図である。
【図７】再生のため使用されるラスタの説明図である。
【符号の説明】
１ 静止磁界発生システム ３，５，７ 傾斜コイルシステム ４ 搬送台上部 ６ 発振器 ８ スイッチ ９ 受信器 １０ 患者 １１ ＲＦコイル １２ 制御ユニット １３ ＦＦＴ演算ユニット １４ 再生ユニット １５ 出力ユニット

Claims (11)

    【特許請求の範囲】
  1. 【請求項１】 検査対象に関係する複数の磁気共鳴デー
    タセットが時間的に連続して取得される動きアーティフ
    ァクトを除去する磁気共鳴方法において、 各低解像度磁気共鳴画像（Ｉi ）が上記磁気共鳴データ
    セットから再生され、 画像変換パラメータ（Ｔij）が上記低解像度磁気共鳴画
    像（Ｉi ）の比較から決定され、 高解像度磁気共鳴画像（Ｉ）が上記画像変換パラメータ
    （Ｔij）に依存して上記磁気共鳴データセットから再生
    される磁気共鳴方法。
  2. 【請求項２】 基準画像（Ｉref ）としての役割を果た
    す低解像度磁気共鳴画像（Ｉ1 ）内の画素の位置が決定
    され、 別の低解像度磁気共鳴画像（Ｉi ）内の上記画素の位置
    が決定され、 上記画像変換パラメータ（Ｔij）が上記基準画像（Ｉ
    ref ）内の画素の位置と、上記別の低解像度磁気共鳴画
    像（Ｉi ）内の上記画素の位置との比較から決定される
    ことを特徴とする請求項１記載の磁気共鳴方法。
  3. 【請求項３】 上記基準画像（Ｉref ）内において、画
    素の位置は上記高解像度磁気共鳴画像（Ｉ）が再生され
    るラスタ（Ｒl ）のラスタ点（Ｒkl）に対応した位置の
    中で決定されることを特徴とする請求項２記載の磁気共
    鳴方法。
  4. 【請求項４】 上記高解像度磁気共鳴画像（Ｉ）の画像
    値は各細分画像の値の加算により決定され、 ここで、第１の細分画像の値は基本ラスタ（Ｒ1 ）のラ
    スタ点での再生により第１の磁気共鳴データセットから
    決定され、別の細分画像の値は関連した画像変換パラメ
    ータ（Ｔij）を用いて上記基本ラスタ（Ｒ1 ）から得ら
    れたラスタ（Ｒ 2 ，Ｒ3 ）のラスタ点での再生により別
    の各磁気共鳴データセット毎に決定されることを特徴と
    する請求項１記載の磁気共鳴方法。
  5. 【請求項５】 上記高解像度磁気共鳴画像（Ｉ）を再生
    するため、磁気共鳴基本画像が基本ラスタ（Ｒ1 ）に関
    する磁気共鳴データセットから形成され、夫々の磁気共
    鳴画像は関連した画像変換パラメータ（Ｔij）を用いて
    上記基本ラスタ（Ｒ1 ）から得られたラスタ（Ｒ2 ，Ｒ
    3 ）上の残りの磁気共鳴データセットから形成され、 上記磁気共鳴基本画像及び上記磁気共鳴画像は高解像度
    磁気共鳴画像（Ｉ）を形成するため重ね合わされること
    を特徴とする請求項１記載の磁気共鳴方法。
  6. 【請求項６】 磁気共鳴データセットを取得するため、
    ｋ空間の零点を通過して放射状に延在する経路（Ｓij）
    方向にｋ空間が走査されることを特徴とする請求項１記
    載の磁気共鳴方法。
  7. 【請求項７】 磁気共鳴データセットを取得するため追
    跡された上記経路（Ｓij）は、別の磁気共鳴データセッ
    トを取得するため追跡される経路（Ｓij）とある角度
    （αij）で交差することを特徴とする請求項６記載の磁
    気共鳴方法。
  8. 【請求項８】 検査対象に関係する複数の磁気共鳴デー
    タセットが時間的に連続して取得される動きアーティフ
    ァクトを除去する磁気共鳴方法において、 上記磁気共鳴データセットの取得と並行して検査対象の
    所定の位置に配置されたマイクロコイルの位置が決定さ
    れ、 画像変換パラメータ（Ｔij）が上記マイクロコイルの時
    間的に連続した位置データの比較から判定され、 高解像度磁気共鳴データ（Ｉ）が上記画像変換パラメー
    タ（Ｔij）に依存して上記磁気共鳴データセットから再
    生されることを特徴とする磁気共鳴方法。
  9. 【請求項９】 一つ以上の低解像度磁気共鳴画像
    （Ｉi ）及び／又は高解像度磁気共鳴画像（Ｉ）が表示
    されることを特徴とする請求項１又は８記載の磁気共鳴
    方法。
  10. 【請求項１０】 上記高解像度磁気共鳴画像は、上記画
    像変換パラメータ（Ｔij）を考慮することなく上記磁気
    共鳴データセットから再生されることを特徴とする請求
    項９記載の磁気共鳴方法。
  11. 【請求項１１】 磁気共鳴信号から磁気共鳴画像を発生
    する手段と、 上記手段を制御するプログラマブル制御ユニット（１
    ２）と、 投影再生法を用いて磁気共鳴画像を再生する再生ユニッ
    ト（１４）とからなる磁気共鳴装置において、 上記制御ユニット（１２）は、検査対象の複数の磁気共
    鳴データセットが時間的に連続して得られるようにプロ
    グラムされ、 上記再生ユニット（１４）は、磁気共鳴データセットか
    ら各低解像度磁気共鳴画像（Ｉi ）を再生し、上記低解
    像度磁気共鳴画像（Ｉi ）の比較から画像変換パラメー
    タ（Ｔij）を取得し、上記画像変換パラメータ（Ｔij）
    に依存して上記磁気共鳴データセットから高解像度磁気
    共鳴画像（Ｉ）を再生することを特徴とする請求項１記
    載の磁気共鳴方法を実施する磁気共鳴装置。
JP9312090A 1996-11-16 1997-11-13 動きアーティファクトを除去する磁気共鳴方法及び装置 Pending JPH10155767A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE19647537:6 1996-11-16
DE19647537A DE19647537A1 (de) 1996-11-16 1996-11-16 MR-Verfahren zur Reduzierung von Bewegungsartefakten und Anordnung zur Durchführung des Verfahrens

Publications (1)

Publication Number Publication Date
JPH10155767A true JPH10155767A (ja) 1998-06-16

Family

ID=7811928

Family Applications (1)

Application Number Title Priority Date Filing Date
JP9312090A Pending JPH10155767A (ja) 1996-11-16 1997-11-13 動きアーティファクトを除去する磁気共鳴方法及び装置

Country Status (4)

Country Link
US (1) US5933006A (ja)
EP (1) EP0843178B1 (ja)
JP (1) JPH10155767A (ja)
DE (2) DE19647537A1 (ja)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001190522A (ja) * 1999-12-10 2001-07-17 Siemens Ag 磁気共鳴トモグラフィ装置を作動させる方法
JP2008529643A (ja) * 2005-02-11 2008-08-07 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 短いエコー時間での磁気共鳴イメージング
JP2009500116A (ja) * 2005-07-08 2009-01-08 ウイスコンシン アラムナイ リサーチ フオンデーシヨン アンダーサンプリングされたｍｒｉの限定的な逆投影再構成法
JP2009508660A (ja) * 2005-09-22 2009-03-05 ウイスコンシン アラムナイ リサーチ フオンデーシヨン 機能的磁気共鳴イメージング用の、画像の取得及び再構成の方法
JP2009508654A (ja) * 2005-09-22 2009-03-05 ウイスコンシン アラムナイ リサーチ フオンデーシヨン 高度に限定された磁気共鳴スペクトロスコピー画像再構成方法
JP2009508634A (ja) * 2005-09-22 2009-03-05 ウイスコンシン アラムナイ リサーチ フオンデーシヨン 拍動している心臓の画像の再構成法
JP2009508656A (ja) * 2005-09-22 2009-03-05 ウイスコンシン アラムナイ リサーチ フオンデーシヨン 高度に限定された画像再構成法を使用する拡散テンソル・イメージング
JP2009508635A (ja) * 2005-09-22 2009-03-05 ウイスコンシン アラムナイ リサーチ フオンデーシヨン 運動コード化ｍｒ画像の高度に限定された再構成

Families Citing this family (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6097977A (en) * 1997-12-31 2000-08-01 General Electric Company Method for controlling data acquisition and image reconstruction during continuous MR imaging
US6404196B1 (en) 1998-07-31 2002-06-11 Case Western Reserve University Method for correction of MRI motion artifacts and main field fluctuation
DE10036207B4 (de) 2000-07-25 2006-11-30 Siemens Ag Verfahren zum Durchführen einer Perfusionsmessung mittels Magnetresonanzbildgebung
DE10038669C1 (de) * 2000-08-08 2002-01-31 Siemens Ag Verfahren zum Betrieb eines Magnetresonanzgeräts mit einem Bewegungsmodell
GB0109892D0 (en) * 2001-04-20 2001-06-13 Secr Defence Method and apparatus for reducing the effects of motion in an image
US7283859B2 (en) * 2001-04-20 2007-10-16 Brigham And Womens' Hospital, Inc. Artifact suppression in dynamic magnetic resonance imaging
US6617850B2 (en) 2001-12-03 2003-09-09 Mayo Foundation For Medical Education And Research Motion correction of magnetic resonance images using phase difference of two orthogonal acquisitions
GB0129465D0 (en) * 2001-12-08 2002-01-30 Qinetiq Ltd Method for compensating for effects of object motion in an image
US6904118B2 (en) * 2002-07-23 2005-06-07 General Electric Company Method and apparatus for generating a density map using dual-energy CT
DE10337932B4 (de) * 2003-08-18 2009-02-05 Siemens Ag Gerät und Verfahren zur Minimierung von Streifenartefakten bei radialer oder spiralförmiger k-Raum-Abtastung in der Magnetresonanzbildgebung
CN101300600B (zh) * 2005-07-08 2016-01-20 威斯康星校友研究基金会 用于ct成像的反投影重构方法
ATE492006T1 (de) * 2005-07-08 2011-01-15 Wisconsin Alumni Res Found Rückprojektions-rekonstruktionsverfahren für ct- bildgebung
US7408347B2 (en) * 2005-09-22 2008-08-05 Wisconsin Alumni Research Foundation Highly constrained magnetic resonance spectroscopy image reconstruction method
EP1959397B1 (en) * 2007-02-19 2019-08-07 Wisconsin Alumni Research Foundation Iterative HYPR medical image reconstruction
ATE542196T1 (de) * 2007-02-19 2012-02-15 Wisconsin Alumni Res Found Verfahren zur lokalisierten und stark eingeschränkten bildrekonstruktion
DE102007009209B3 (de) * 2007-02-26 2008-08-28 Siemens Ag Verfahren zur Erzeugung von zeitlich hochaufgelösten MR-Aufnahmen durch Multiplikation unterabgetasteter Bilder mit einem Mittelwertbild
EP2233941A1 (en) * 2009-03-25 2010-09-29 Koninklijke Philips Electronics N.V. MR imaging using parallel signal acquisition
EP2824472A1 (en) 2013-07-10 2015-01-14 Siemens Aktiengesellschaft MRI with self-navigated motion compensation
US8620404B2 (en) 2011-07-26 2013-12-31 Wisconsin Alumni Research Foundation System and method of high-frame rate, time-resolved, three-dimensional magnetic resonance angiograpy
US10732244B2 (en) 2012-03-26 2020-08-04 Sirona Dental Systems Gmbh Systems, methods, apparatuses, and computer-readable storage media for performing diagnostic examinations using MRI
US10126457B2 (en) 2016-03-04 2018-11-13 Baker Hughes, A Ge Company, Llc Motion detection and correction of magnetic resonance data
EP3719526A1 (en) * 2019-04-03 2020-10-07 Siemens Healthcare GmbH Magnetic resonance imaging and motion detection

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IL86231A (en) * 1988-04-29 1991-07-18 Elscint Ltd Correction for eddy current caused phase degradation
US4937526A (en) * 1988-11-23 1990-06-26 Mayo Foundation For Medical Education And Research Adaptive method for reducing motion and flow artifacts in NMR images
JPH05154130A (ja) * 1991-12-06 1993-06-22 Hitachi Ltd 体動アーチファクト除去方法
JPH05245128A (ja) * 1992-03-05 1993-09-24 Hitachi Ltd 核磁気共鳴装置における体動補正方法
DE4319539A1 (de) * 1993-06-12 1994-12-15 Philips Patentverwaltung Verfahren zur Erzeugung einer MR-Bildfolge und Anordnung zur Durchführung des Verfahrens
US5377680A (en) * 1993-08-04 1995-01-03 General Electric Company MRI cardiac image produced by temporal data sharing
US5742163A (en) * 1996-04-26 1998-04-21 Picker International, Inc. Magnetic resonance scan calibration and reconstruction technique for multi-shot, multi-echo imaging
US5729140A (en) * 1996-10-04 1998-03-17 Mayo Foundation For Medical Education And Research Suppression of artifacts in NMR images by correlation of two NEX acquistions

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001190522A (ja) * 1999-12-10 2001-07-17 Siemens Ag 磁気共鳴トモグラフィ装置を作動させる方法
JP2008529643A (ja) * 2005-02-11 2008-08-07 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 短いエコー時間での磁気共鳴イメージング
JP2009500116A (ja) * 2005-07-08 2009-01-08 ウイスコンシン アラムナイ リサーチ フオンデーシヨン アンダーサンプリングされたｍｒｉの限定的な逆投影再構成法
JP2009508660A (ja) * 2005-09-22 2009-03-05 ウイスコンシン アラムナイ リサーチ フオンデーシヨン 機能的磁気共鳴イメージング用の、画像の取得及び再構成の方法
JP2009508654A (ja) * 2005-09-22 2009-03-05 ウイスコンシン アラムナイ リサーチ フオンデーシヨン 高度に限定された磁気共鳴スペクトロスコピー画像再構成方法
JP2009508634A (ja) * 2005-09-22 2009-03-05 ウイスコンシン アラムナイ リサーチ フオンデーシヨン 拍動している心臓の画像の再構成法
JP2009508656A (ja) * 2005-09-22 2009-03-05 ウイスコンシン アラムナイ リサーチ フオンデーシヨン 高度に限定された画像再構成法を使用する拡散テンソル・イメージング
JP2009508655A (ja) * 2005-09-22 2009-03-05 ウイスコンシン アラムナイ リサーチ フオンデーシヨン 心臓ゲート磁気共鳴イメージングのための画像再構成法
JP2009508635A (ja) * 2005-09-22 2009-03-05 ウイスコンシン アラムナイ リサーチ フオンデーシヨン 運動コード化ｍｒ画像の高度に限定された再構成

Also Published As

Publication number Publication date
DE59712195D1 (de) 2005-03-17
US5933006A (en) 1999-08-03
EP0843178A1 (de) 1998-05-20
DE19647537A1 (de) 1998-05-20
EP0843178B1 (de) 2005-02-09

Similar Documents

Publication Publication Date Title
JPH10155767A (ja) 動きアーティファクトを除去する磁気共鳴方法及び装置
US6118273A (en) Magnetic resonance imaging method and device therefor
JP4612000B2 (ja) 磁気共鳴イメージング装置
EP1377214B1 (en) System for acquiring mri data from a large field of view using continuous table motion
US9684979B2 (en) MRI 3D cine imaging based on intersecting source and anchor slice data
JP5960163B2 (ja) 校正スキャン、コイル感度マップ及びナビゲータを使用する、剛体動き補償のためのパラレルｍｒｉ方法
US5502385A (en) Method of forming a series of MR images, and device for carrying out the method
US20110105884A1 (en) Mri involving dynamic profile sharing such as keyhole and motion correction
JPS63214247A (ja) 像を発生する方法
US5923789A (en) Band limited interpolation and projection of spatial 3-D images
JPH06217960A (ja) 映画磁気共鳴撮像方法及び装置
JP2005144173A (ja) Ｍｒｔ画像形成方法
US6404196B1 (en) Method for correction of MRI motion artifacts and main field fluctuation
JP3699304B2 (ja) 磁気共鳴撮像装置
CA2469031A1 (en) Method for compensating for motion artifacts in magnetic resonance imaging by phase line prediction
JP2713160B2 (ja) Ｍｒイメージング装置
EP0217196B1 (en) Apparatus and method for initially capturing most significant low spatial frequency nmr imaging data
JP3051374B2 (ja) 磁気共鳴撮像装置
EP1444529B1 (en) Magnetic resonance method for forming a fast dynamic image
US5498961A (en) MR method for two-dimensional or three-dimensional imaging of an examination zone, and apparatus for carrying out the method
US7626388B2 (en) Magnetic resonance imager
JP2002177245A (ja) Ｍｒ画像を再構成する方法
US20070001674A1 (en) Object motion correction during MR imaging
JP4512587B2 (ja) アンダーサンプル磁気共鳴画像化
US4966149A (en) Reducing the effects of coherence in magnetic resonance imaging

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20041112

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20060831

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20060905

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20061204

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20061207

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20070302

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20070724

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20071218