JPH08173396A - Ｍｒイメージング方法およびｍｒｉ装置 - Google Patents
Ｍｒイメージング方法およびｍｒｉ装置Info
- Publication number
- JPH08173396A JPH08173396A JP6317872A JP31787294A JPH08173396A JP H08173396 A JPH08173396 A JP H08173396A JP 6317872 A JP6317872 A JP 6317872A JP 31787294 A JP31787294 A JP 31787294A JP H08173396 A JPH08173396 A JP H08173396A
- Authority
- JP
- Japan
- Prior art keywords
- data
- phase
- view
- imaging
- axis
- 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.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 114
- 238000002592 echocardiography Methods 0.000 claims abstract description 30
- 238000005070 sampling Methods 0.000 claims abstract description 17
- 238000002595 magnetic resonance imaging Methods 0.000 description 24
- 238000005314 correlation function Methods 0.000 description 19
- 238000010586 diagram Methods 0.000 description 19
- 230000033001 locomotion Effects 0.000 description 17
- 238000000034 method Methods 0.000 description 16
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 15
- 238000005481 NMR spectroscopy Methods 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
(57)【要約】 （修正有）
【目的】 移動する被検体を撮像したときのアーチファ
クトを抑制する。 【構成】 ｙ軸に位相エンコード勾配ｐｅを印加し、ｘ
軸に周波数エンコード勾配ｒｉを印加しながらイメージ
エコーＩＥをサンプリングしｐビューのイメージング用
データを収集する。次にリワインド勾配ｐｒを印加して
位相エンコード量を“０”に戻し、エンコード勾配ｒｉ
と極性を反転したサンプリング数の読み出し勾配ｒｎ１
を印加しながらナビゲーションエコーＮＥ１をサンプリ
ングしｐビューのイメージング用データに対応する周波
数軸方向のナビゲーション用データを収集する。これを
勾配ｐｅ及びｐｒを変えながら繰り返し、Ｍビュー分と
ナビゲーション用の各データを収集する。各ビューのナ
ビゲーション用データの相互関係から位相補正量を算出
して、各ビューのイメージング用データを位相補正しそ
のデータから画像再構成する。
クトを抑制する。 【構成】 ｙ軸に位相エンコード勾配ｐｅを印加し、ｘ
軸に周波数エンコード勾配ｒｉを印加しながらイメージ
エコーＩＥをサンプリングしｐビューのイメージング用
データを収集する。次にリワインド勾配ｐｒを印加して
位相エンコード量を“０”に戻し、エンコード勾配ｒｉ
と極性を反転したサンプリング数の読み出し勾配ｒｎ１
を印加しながらナビゲーションエコーＮＥ１をサンプリ
ングしｐビューのイメージング用データに対応する周波
数軸方向のナビゲーション用データを収集する。これを
勾配ｐｅ及びｐｒを変えながら繰り返し、Ｍビュー分と
ナビゲーション用の各データを収集する。各ビューのナ
ビゲーション用データの相互関係から位相補正量を算出
して、各ビューのイメージング用データを位相補正しそ
のデータから画像再構成する。
Description
【０００１】
【産業上の利用分野】この発明は、ＭＲ（Magnetic Res
onance）イメージング方法およびＭＲＩ（Magnetic Res
onance Imaging）装置に関する。さらに詳しくは、移動
する被検体を撮像し且つアーチファクト（Artifact）を
抑制してイメージングすることが出来るＭＲイメージン
グ方法およびＭＲＩ装置に関する。
onance）イメージング方法およびＭＲＩ（Magnetic Res
onance Imaging）装置に関する。さらに詳しくは、移動
する被検体を撮像し且つアーチファクト（Artifact）を
抑制してイメージングすることが出来るＭＲイメージン
グ方法およびＭＲＩ装置に関する。
【０００２】
【従来の技術】図１４に、従来のＭＲＩ装置の一例の構
成概念図を示す。このＭＲＩ装置５００において、５０
は撮像装置、５０ａはボア、６０は被検体Ｇを移動させ
るクレードルである。５１は撮像時にクレードル６０を
停止させる移動停止制御部である。このＭＲＩ装置５０
０では、被検体Ｇの移動を停止して撮像し、移動方向の
長さＴの第１の画像を得る。次に、距離Ｔだけ移動して
から再び停止して撮像し、移動方向の長さＴの第２の画
像を得る。これをｎ回繰り返し、得られた第１から第ｎ
までのｎ枚の画像を継ぎ合せて、長さｎＴの画像を得
る。上記のＭＲＩ装置５００は、例えば特開昭６３−２
７２３３５号公報に開示されている。
成概念図を示す。このＭＲＩ装置５００において、５０
は撮像装置、５０ａはボア、６０は被検体Ｇを移動させ
るクレードルである。５１は撮像時にクレードル６０を
停止させる移動停止制御部である。このＭＲＩ装置５０
０では、被検体Ｇの移動を停止して撮像し、移動方向の
長さＴの第１の画像を得る。次に、距離Ｔだけ移動して
から再び停止して撮像し、移動方向の長さＴの第２の画
像を得る。これをｎ回繰り返し、得られた第１から第ｎ
までのｎ枚の画像を継ぎ合せて、長さｎＴの画像を得
る。上記のＭＲＩ装置５００は、例えば特開昭６３−２
７２３３５号公報に開示されている。
【０００３】
【発明が解決しようとする課題】上記従来のＭＲＩ装置
５００では、撮像装置５０により撮像をする時はクレー
ドル６０は停止し、クレードル６０により被検体Ｇを移
動させる時は撮像装置５０は停止している。従って、撮
像効率がよくないという問題点がある。そこで、この発
明の目的は、移動する被検体を撮像し且つアーチファク
トを抑制してイメージングすることが出来るようにした
ＭＲイメージング方法およびＭＲＩ装置を提供すること
にある。
５００では、撮像装置５０により撮像をする時はクレー
ドル６０は停止し、クレードル６０により被検体Ｇを移
動させる時は撮像装置５０は停止している。従って、撮
像効率がよくないという問題点がある。そこで、この発
明の目的は、移動する被検体を撮像し且つアーチファク
トを抑制してイメージングすることが出来るようにした
ＭＲイメージング方法およびＭＲＩ装置を提供すること
にある。
【０００４】
【課題を解決するための手段】第１の観点では、この発
明は、ｋ−空間の周波数軸データ番号をｆ（−Ｎ／２≦
ｆ≦Ｎ／２−１）とし且つ位相軸データ番号をｐ（−Ｍ
／２≦ｐ≦Ｍ／２−１）とするとき、ＲＦパルスを印加
し且つ位相エンコードを加え且つ周波数エンコードを加
え、発生するエコーからイメージング用データＤ’
（ｆ，ｐ）を収集することを繰り返して、Ｍビュー分の
イメージング用データＤ’（ｆ，−Ｍ／２）〜Ｄ’
（ｆ，Ｍ／２−１）を収集するＭＲイメージング方法に
おいて、１つのＲＦパルスに対して少なくとも２つのエ
コーを発生させ、１つのエコーに対しては位相エンコー
ド量を“０”とし、周波数軸にサンプリング数Ｎ’の読
み出し勾配を印加しながら、データ番号ｆ’（−Ｎ’／
２≦ｆ’≦Ｎ’／２−１）のｐビューにおける周波数軸
方向についてのナビゲーション用データＳＲｐ（ｆ’）
を収集し、残りのエコーに対してはｋ−空間のｐビュー
に応じた位相エンコードを加え且つ周波数エンコードを
加え、これによりｋ−空間のｐビューのイメージング用
データＤ’（ｆ，ｐ）を収集することを繰り返し、Ｍビ
ュー分のイメージング用データＤ’（ｆ，−Ｍ／２）〜
Ｄ’（ｆ，Ｍ／２−１）を収集し、各ビューのナビゲー
ション用データＳＲ-M/2（ｆ’）〜ＳＲM/2-1（ｆ'）の
相互関係から位相補正量Δφｘ（ｐ）を算出し、その位
相補正量Δφｘ（ｐ）で各ビューのイメージング用デー
タＤ’（ｆ，ｐ）を位相補正することを特徴とするＭＲ
イメージング方法を提供する。
明は、ｋ−空間の周波数軸データ番号をｆ（−Ｎ／２≦
ｆ≦Ｎ／２−１）とし且つ位相軸データ番号をｐ（−Ｍ
／２≦ｐ≦Ｍ／２−１）とするとき、ＲＦパルスを印加
し且つ位相エンコードを加え且つ周波数エンコードを加
え、発生するエコーからイメージング用データＤ’
（ｆ，ｐ）を収集することを繰り返して、Ｍビュー分の
イメージング用データＤ’（ｆ，−Ｍ／２）〜Ｄ’
（ｆ，Ｍ／２−１）を収集するＭＲイメージング方法に
おいて、１つのＲＦパルスに対して少なくとも２つのエ
コーを発生させ、１つのエコーに対しては位相エンコー
ド量を“０”とし、周波数軸にサンプリング数Ｎ’の読
み出し勾配を印加しながら、データ番号ｆ’（−Ｎ’／
２≦ｆ’≦Ｎ’／２−１）のｐビューにおける周波数軸
方向についてのナビゲーション用データＳＲｐ（ｆ’）
を収集し、残りのエコーに対してはｋ−空間のｐビュー
に応じた位相エンコードを加え且つ周波数エンコードを
加え、これによりｋ−空間のｐビューのイメージング用
データＤ’（ｆ，ｐ）を収集することを繰り返し、Ｍビ
ュー分のイメージング用データＤ’（ｆ，−Ｍ／２）〜
Ｄ’（ｆ，Ｍ／２−１）を収集し、各ビューのナビゲー
ション用データＳＲ-M/2（ｆ’）〜ＳＲM/2-1（ｆ'）の
相互関係から位相補正量Δφｘ（ｐ）を算出し、その位
相補正量Δφｘ（ｐ）で各ビューのイメージング用デー
タＤ’（ｆ，ｐ）を位相補正することを特徴とするＭＲ
イメージング方法を提供する。
【０００５】第２の観点では、この発明は、ｋ−空間の
周波数軸データ番号をｆ（−Ｎ／２≦ｆ≦Ｎ／２−１）
とし且つ位相軸データ番号をｐ（−Ｍ／２≦ｐ≦Ｍ／２
−１）とするとき、ＲＦパルスを印加し且つ位相エンコ
ードを加え且つ周波数エンコードを加え、発生するエコ
ーからイメージング用データＤ’（ｆ，ｐ）を収集する
ことを繰り返して、Ｍビュー分のイメージング用データ
Ｄ’（ｆ，−Ｍ／２）〜Ｄ’（ｆ，Ｍ／２−１）を収集
するＭＲイメージング方法において、１つのＲＦパルス
に対して少なくとも２つのエコーを発生させ、１つのエ
コーに対しては周波数エンコード量を“０”とし、位相
軸にサンプリング数Ｍ’の読み出し勾配を印加しなが
ら、データ番号ｐ’（−Ｍ’／２≦ｐ’≦Ｍ’／２−
１）のｐビューにおける位相軸方向についてのナビゲー
ション用データＳＷｐ（ｐ’）を収集し、残りのエコー
に対してはｋ−空間のｐビューに応じた位相エンコード
を加え且つ周波数エンコードを加え、これによりｋ−空
間のｐビューのイメージング用データＤ’（ｆ，ｐ）を
収集することを繰り返し、Ｍビュー分のイメージング用
データＤ’（ｆ，−Ｍ／２）〜Ｄ’（ｆ，Ｍ／２−１）
を収集し、各ビューのナビゲーション用データＳＷ-M/2
（ｐ’）〜ＳＷM/2-1（ｐ'）の相互関係から位相補正量
Δφｙ（ｐ）を算出し、その位相補正量Δφｙ（ｐ）で
各ビューのイメージング用データＤ’（ｆ，ｐ）を位相
補正することを特徴とするＭＲイメージング方法を提供
する。
周波数軸データ番号をｆ（−Ｎ／２≦ｆ≦Ｎ／２−１）
とし且つ位相軸データ番号をｐ（−Ｍ／２≦ｐ≦Ｍ／２
−１）とするとき、ＲＦパルスを印加し且つ位相エンコ
ードを加え且つ周波数エンコードを加え、発生するエコ
ーからイメージング用データＤ’（ｆ，ｐ）を収集する
ことを繰り返して、Ｍビュー分のイメージング用データ
Ｄ’（ｆ，−Ｍ／２）〜Ｄ’（ｆ，Ｍ／２−１）を収集
するＭＲイメージング方法において、１つのＲＦパルス
に対して少なくとも２つのエコーを発生させ、１つのエ
コーに対しては周波数エンコード量を“０”とし、位相
軸にサンプリング数Ｍ’の読み出し勾配を印加しなが
ら、データ番号ｐ’（−Ｍ’／２≦ｐ’≦Ｍ’／２−
１）のｐビューにおける位相軸方向についてのナビゲー
ション用データＳＷｐ（ｐ’）を収集し、残りのエコー
に対してはｋ−空間のｐビューに応じた位相エンコード
を加え且つ周波数エンコードを加え、これによりｋ−空
間のｐビューのイメージング用データＤ’（ｆ，ｐ）を
収集することを繰り返し、Ｍビュー分のイメージング用
データＤ’（ｆ，−Ｍ／２）〜Ｄ’（ｆ，Ｍ／２−１）
を収集し、各ビューのナビゲーション用データＳＷ-M/2
（ｐ’）〜ＳＷM/2-1（ｐ'）の相互関係から位相補正量
Δφｙ（ｐ）を算出し、その位相補正量Δφｙ（ｐ）で
各ビューのイメージング用データＤ’（ｆ，ｐ）を位相
補正することを特徴とするＭＲイメージング方法を提供
する。
【０００６】第３の観点では、この発明は、ｋ−空間の
周波数軸データ番号をｆ（−Ｎ／２≦ｆ≦Ｎ／２−１）
とし且つ位相軸データ番号をｐ（−Ｍ／２≦ｐ≦Ｍ／２
−１）とするとき、ＲＦパルスを印加し且つ位相エンコ
ードを加え且つ周波数エンコードを加え、発生するエコ
ーからイメージング用データＤ’（ｆ，ｐ）を収集する
ことを繰り返して、Ｍビュー分のイメージング用データ
Ｄ’（ｆ，−Ｍ／２）〜Ｄ’（ｆ，Ｍ／２−１）を収集
するＭＲＩ装置において、１つのＲＦパルスに対して少
なくとも２つのエコーを発生させ、１つのエコーに対し
ては位相エンコード量を“０”とし、周波数軸にサンプ
リング数Ｎ’の読み出し勾配を印加しながら、データ番
号ｆ’（−Ｎ’／２≦ｆ’≦Ｎ’／２−１）のｐビュー
における周波数軸方向についてのナビゲーション用デー
タＳＲｐ（ｆ’）を収集し、残りのエコーに対してはｋ
−空間のｐビューに応じた位相エンコードを加え且つ周
波数エンコードを加え、これによりｋ−空間のｐビュー
のイメージング用データＤ’（ｆ，ｐ）を収集すること
を繰り返し、Ｍビュー分のイメージング用データＤ’
（ｆ，−Ｍ／２）〜Ｄ’（ｆ，Ｍ／２−１）を収集する
データ収集手段と、各ビューのナビゲーション用データ
ＳＲ-M/2（ｆ’）〜ＳＲM/2-1（ｆ'）の相互関係から位
相補正量Δφｘ（ｐ）を算出し、その位相補正量Δφｘ
（ｐ）で各ビューのイメージング用データＤ’（ｆ，
ｐ）を位相補正する位相補正手段とを具備したことを特
徴とするＭＲＩ装置を提供する。
周波数軸データ番号をｆ（−Ｎ／２≦ｆ≦Ｎ／２−１）
とし且つ位相軸データ番号をｐ（−Ｍ／２≦ｐ≦Ｍ／２
−１）とするとき、ＲＦパルスを印加し且つ位相エンコ
ードを加え且つ周波数エンコードを加え、発生するエコ
ーからイメージング用データＤ’（ｆ，ｐ）を収集する
ことを繰り返して、Ｍビュー分のイメージング用データ
Ｄ’（ｆ，−Ｍ／２）〜Ｄ’（ｆ，Ｍ／２−１）を収集
するＭＲＩ装置において、１つのＲＦパルスに対して少
なくとも２つのエコーを発生させ、１つのエコーに対し
ては位相エンコード量を“０”とし、周波数軸にサンプ
リング数Ｎ’の読み出し勾配を印加しながら、データ番
号ｆ’（−Ｎ’／２≦ｆ’≦Ｎ’／２−１）のｐビュー
における周波数軸方向についてのナビゲーション用デー
タＳＲｐ（ｆ’）を収集し、残りのエコーに対してはｋ
−空間のｐビューに応じた位相エンコードを加え且つ周
波数エンコードを加え、これによりｋ−空間のｐビュー
のイメージング用データＤ’（ｆ，ｐ）を収集すること
を繰り返し、Ｍビュー分のイメージング用データＤ’
（ｆ，−Ｍ／２）〜Ｄ’（ｆ，Ｍ／２−１）を収集する
データ収集手段と、各ビューのナビゲーション用データ
ＳＲ-M/2（ｆ’）〜ＳＲM/2-1（ｆ'）の相互関係から位
相補正量Δφｘ（ｐ）を算出し、その位相補正量Δφｘ
（ｐ）で各ビューのイメージング用データＤ’（ｆ，
ｐ）を位相補正する位相補正手段とを具備したことを特
徴とするＭＲＩ装置を提供する。
【０００７】第４の観点では、この発明は、ｋ−空間の
周波数軸データ番号をｆ（−Ｎ／２≦ｆ≦Ｎ／２−１）
とし且つ位相軸データ番号をｐ（−Ｍ／２≦ｐ≦Ｍ／２
−１）とするとき、ＲＦパルスを印加し且つ位相エンコ
ードを加え且つ周波数エンコードを加え、発生するエコ
ーからイメージング用データＤ’（ｆ，ｐ）を収集する
ことを繰り返して、Ｍビュー分のイメージング用データ
Ｄ’（ｆ，−Ｍ／２）〜Ｄ’（ｆ，Ｍ／２−１）を収集
するＭＲＩ装置において、１つのＲＦパルスに対して少
なくとも２つのエコーを発生させ、１つのエコーに対し
ては周波数エンコード量を“０”とし、位相軸にサンプ
リング数Ｍ’の読み出し勾配を印加しながら、データ番
号ｐ’（−Ｍ’／２≦ｐ’≦Ｍ’／２−１）のｐビュー
における位相軸方向についてのナビゲーション用データ
ＳＷｐ（ｐ’）を収集し、残りのエコーに対してはｋ−
空間のｐビューに応じた位相エンコードを加え且つ周波
数エンコードを加え、これによりｋ−空間のｐビューの
イメージング用データＤ’（ｆ，ｐ）を収集することを
繰り返し、Ｍビュー分のイメージング用データＤ’
（ｆ，−Ｍ／２）〜Ｄ’（ｆ，Ｍ／２−１）を収集する
データ収集手段と、各ビューのナビゲーション用データ
ＳＷ-M/2（ｐ’）〜ＳＷM/2-1（ｐ'）の相互関係から位
相補正量Δφｙ（ｐ）を算出し、その位相補正量Δφｙ
（ｐ）で各ビューのイメージング用データＤ’（ｆ，
ｐ）を位相補正する位相補正手段とを具備したことを特
徴とするＭＲＩ装置を提供する。
周波数軸データ番号をｆ（−Ｎ／２≦ｆ≦Ｎ／２−１）
とし且つ位相軸データ番号をｐ（−Ｍ／２≦ｐ≦Ｍ／２
−１）とするとき、ＲＦパルスを印加し且つ位相エンコ
ードを加え且つ周波数エンコードを加え、発生するエコ
ーからイメージング用データＤ’（ｆ，ｐ）を収集する
ことを繰り返して、Ｍビュー分のイメージング用データ
Ｄ’（ｆ，−Ｍ／２）〜Ｄ’（ｆ，Ｍ／２−１）を収集
するＭＲＩ装置において、１つのＲＦパルスに対して少
なくとも２つのエコーを発生させ、１つのエコーに対し
ては周波数エンコード量を“０”とし、位相軸にサンプ
リング数Ｍ’の読み出し勾配を印加しながら、データ番
号ｐ’（−Ｍ’／２≦ｐ’≦Ｍ’／２−１）のｐビュー
における位相軸方向についてのナビゲーション用データ
ＳＷｐ（ｐ’）を収集し、残りのエコーに対してはｋ−
空間のｐビューに応じた位相エンコードを加え且つ周波
数エンコードを加え、これによりｋ−空間のｐビューの
イメージング用データＤ’（ｆ，ｐ）を収集することを
繰り返し、Ｍビュー分のイメージング用データＤ’
（ｆ，−Ｍ／２）〜Ｄ’（ｆ，Ｍ／２−１）を収集する
データ収集手段と、各ビューのナビゲーション用データ
ＳＷ-M/2（ｐ’）〜ＳＷM/2-1（ｐ'）の相互関係から位
相補正量Δφｙ（ｐ）を算出し、その位相補正量Δφｙ
（ｐ）で各ビューのイメージング用データＤ’（ｆ，
ｐ）を位相補正する位相補正手段とを具備したことを特
徴とするＭＲＩ装置を提供する。
【０００８】第５の観点では、この発明は、ｋ−空間の
第１のエンコード軸のエンコード量を“０”とし、ｋ−
空間の第２のエンコード軸に読み出し勾配を印加しなが
らナビゲーション用データＳ１を収集し，その後，同様
にしてナビゲーション用データＳ２を収集するデータ収
集手段と、前記ナビゲーション用データＳ１の１次元フ
ーリエ変換の絶対値から求まる投影プロファイルと前記
ナビゲーション用データＳ２の１次元フーリエ変換の絶
対値から求まる投影プロファイルの相互相関から前記第
２のエンコード軸に対応する実空間の方向についての信
号源の移動量Δを算出する移動量算出手段とを具備した
ことを特徴とするＭＲイメージング方法を提供する。
第１のエンコード軸のエンコード量を“０”とし、ｋ−
空間の第２のエンコード軸に読み出し勾配を印加しなが
らナビゲーション用データＳ１を収集し，その後，同様
にしてナビゲーション用データＳ２を収集するデータ収
集手段と、前記ナビゲーション用データＳ１の１次元フ
ーリエ変換の絶対値から求まる投影プロファイルと前記
ナビゲーション用データＳ２の１次元フーリエ変換の絶
対値から求まる投影プロファイルの相互相関から前記第
２のエンコード軸に対応する実空間の方向についての信
号源の移動量Δを算出する移動量算出手段とを具備した
ことを特徴とするＭＲイメージング方法を提供する。
【０００９】
【作用】図１２に示すように、−Ｆｗ／２≦ｘ≦Ｆｗ／
２，−Ｐｗ／２≦ｙ≦Ｐｗ／２の実空間（ｘ，ｙ）に信
号源ｇ（ｘ，ｙ）があるとき、収集されるｋ−空間のデ
ータＤ（ｆ，ｐ）は、次式のようになる。但し、ｘ軸方
向を周波数軸方向とし、ｙ軸方向を位相軸方向とする
（なお、ｙ軸，ｘ軸は、位相軸または周波数軸のいずれ
かに対応して規定される実空間の軸である。一方、後述
する勾配磁場コイルのＺ軸，Ｙ軸，Ｘ軸は、マグネット
アセンブリ１の構造として規定される実空間の軸であ
る）。
２，−Ｐｗ／２≦ｙ≦Ｐｗ／２の実空間（ｘ，ｙ）に信
号源ｇ（ｘ，ｙ）があるとき、収集されるｋ−空間のデ
ータＤ（ｆ，ｐ）は、次式のようになる。但し、ｘ軸方
向を周波数軸方向とし、ｙ軸方向を位相軸方向とする
（なお、ｙ軸，ｘ軸は、位相軸または周波数軸のいずれ
かに対応して規定される実空間の軸である。一方、後述
する勾配磁場コイルのＺ軸，Ｙ軸，Ｘ軸は、マグネット
アセンブリ１の構造として規定される実空間の軸であ
る）。
【００１０】
【数１】
【００１１】次に、図１３に示すように、上記信号源ｇ
（ｘ，ｙ）がｘ軸方向にΔｘだけ移動した場合を想定す
ると、収集されるｋ−空間のデータＤ’（ｆ，ｐ）は、
次式のようになる。
（ｘ，ｙ）がｘ軸方向にΔｘだけ移動した場合を想定す
ると、収集されるｋ−空間のデータＤ’（ｆ，ｐ）は、
次式のようになる。
【００１２】
【数２】
【００１３】ここで、Ｘ＝ｘ−Δｘとすれば、（数２）
式は、次式のようになる。
式は、次式のようになる。
【００１４】
【数３】
【００１５】（数３）より、データＤ’（ｆ，ｐ）に対
して、 Ｄ（ｆ，ｐ）＝exp｛ｊ２πｆΔｘ／Ｆｗ｝×Ｄ’（ｆ，ｐ） のように位相補正を行い、２次元フーリエ変換を行え
ば、図１３に示すように、信号源ｇ（ｘ，ｙ）が移動し
ていない場合に相当するイメージングが得られることが
判る。ここで、移動量Δｘがビュー毎に異なるなら、ビ
ュー毎の移動量Δｘｐを用いて、 Ｄ（ｆ，ｐ）＝exp｛ｊ２πｆΔｘｐ／Ｆｗ｝×Ｄ’（ｆ，ｐ） …（４） により全てのビューのＤ’（ｆ，ｐ）を位相補正し、得
られたＤ（ｆ，ｐ）に対して２次元フーリエ変換を行え
ば、信号源ｇ（ｘ，ｙ）が移動していない場合に相当す
るイメージが得られる。すなわち、移動する被検体を撮
像し、且つ、アーチファクトを抑制して、イメージング
することが出来る。
して、 Ｄ（ｆ，ｐ）＝exp｛ｊ２πｆΔｘ／Ｆｗ｝×Ｄ’（ｆ，ｐ） のように位相補正を行い、２次元フーリエ変換を行え
ば、図１３に示すように、信号源ｇ（ｘ，ｙ）が移動し
ていない場合に相当するイメージングが得られることが
判る。ここで、移動量Δｘがビュー毎に異なるなら、ビ
ュー毎の移動量Δｘｐを用いて、 Ｄ（ｆ，ｐ）＝exp｛ｊ２πｆΔｘｐ／Ｆｗ｝×Ｄ’（ｆ，ｐ） …（４） により全てのビューのＤ’（ｆ，ｐ）を位相補正し、得
られたＤ（ｆ，ｐ）に対して２次元フーリエ変換を行え
ば、信号源ｇ（ｘ，ｙ）が移動していない場合に相当す
るイメージが得られる。すなわち、移動する被検体を撮
像し、且つ、アーチファクトを抑制して、イメージング
することが出来る。
【００１６】上記信号源ｇ（ｘ，ｙ）がｙ軸方向にΔｙ
ｐだけ移動した場合も上記と同様であり、 Ｄ（ｆ，ｐ）＝exp｛ｊ２πｐΔｙｐ／Ｐｗ｝×Ｄ’（ｆ，ｐ） …（５） により全てのビューのＤ’（ｆ，ｐ）を位相補正し、得
られたＤ（ｆ，ｐ）に対して２次元フーリエ変換を行え
ば、信号源ｇ（ｘ，ｙ）が移動していない場合に相当す
るイメージが得られる。
ｐだけ移動した場合も上記と同様であり、 Ｄ（ｆ，ｐ）＝exp｛ｊ２πｐΔｙｐ／Ｐｗ｝×Ｄ’（ｆ，ｐ） …（５） により全てのビューのＤ’（ｆ，ｐ）を位相補正し、得
られたＤ（ｆ，ｐ）に対して２次元フーリエ変換を行え
ば、信号源ｇ（ｘ，ｙ）が移動していない場合に相当す
るイメージが得られる。
【００１７】さらに、上記信号源ｇ（ｘ，ｙ）がｘ軸方
向にΔｘｐだけ移動し、ｙ軸方向にΔｙｐだけ移動した
場合も上記と同様であり、 Ｄ（ｆ，ｐ）＝exp｛ｊ２π（ｆΔｘｐ／Ｆｗ＋ｐΔｙｐ／Ｐｗ）｝ ×Ｄ’（ｆ，ｐ） …（６） により全てのビューのＤ’（ｆ，ｐ）を位相補正し、得
られたＤ（ｆ，ｐ）に対し２次元フーリエ変換を行え
ば、信号源ｇ（ｘ，ｙ）が移動していない場合に相当す
るイメージが得られる。
向にΔｘｐだけ移動し、ｙ軸方向にΔｙｐだけ移動した
場合も上記と同様であり、 Ｄ（ｆ，ｐ）＝exp｛ｊ２π（ｆΔｘｐ／Ｆｗ＋ｐΔｙｐ／Ｐｗ）｝ ×Ｄ’（ｆ，ｐ） …（６） により全てのビューのＤ’（ｆ，ｐ）を位相補正し、得
られたＤ（ｆ，ｐ）に対し２次元フーリエ変換を行え
ば、信号源ｇ（ｘ，ｙ）が移動していない場合に相当す
るイメージが得られる。
【００１８】さて、上記第１の観点によるＭＲイメージ
ング方法および上記第３の観点によるＭＲＩ装置では、
１つのＲＦパルスに対して少なくとも２つのエコーを発
生させる。その内の１つのエコーに対しては、位相エン
コード量を“０”とし、周波数軸にサンプリング数Ｎ’
の読み出し勾配を印加しながら、データ番号ｆ’（−
Ｎ’／２≦ｆ’≦Ｎ’／２−１）のｐビューにおける周
波数軸方向についてのナビゲーション用データＳＲｐ
（ｆ’）を収集する。また、残りのエコーに対しては、
ｋ−空間のｐビューに応じた位相エンコードを加え且つ
周波数エンコードを加え、これによりｋ−空間のｐビュ
ーのイメージング用データＤ’（ｆ，ｐ）を順に収集す
る。そして、各ビューのナビゲーション用データＳＲｐ
（ｆ’）の相互関係から位相補正量Δφｘ（ｐ）を算出
し、その位相補正量Δφｘ（ｐ）で各ビューのイメージ
ング用データＤ’（ｆ，ｐ）を位相補正する。上記ナビ
ゲーション用データＳＲｐ（ｆ’）は、常に位相エンコ
ード量“０”で収集しているので、その相互関係から各
ビューでの信号源のｘ軸方向の移動量Δｘｐを求めるこ
とが出来る。ここで、 ２πｆΔｘｐ／Ｆｗ＝Δφｘ（ｐ） とすれば、前記（４）式となり、位相補正して得られた
Ｄ（ｆ，ｐ）に対して２次元フーリエ変換を行えば、信
号源ｇ（ｘ，ｙ）が移動していない場合に相当するイメ
ージが得られる。すなわち、ｘ軸方向に移動する被検体
を撮像し、且つ、アーチファクトを抑制して、イメージ
ングすることが出来る。
ング方法および上記第３の観点によるＭＲＩ装置では、
１つのＲＦパルスに対して少なくとも２つのエコーを発
生させる。その内の１つのエコーに対しては、位相エン
コード量を“０”とし、周波数軸にサンプリング数Ｎ’
の読み出し勾配を印加しながら、データ番号ｆ’（−
Ｎ’／２≦ｆ’≦Ｎ’／２−１）のｐビューにおける周
波数軸方向についてのナビゲーション用データＳＲｐ
（ｆ’）を収集する。また、残りのエコーに対しては、
ｋ−空間のｐビューに応じた位相エンコードを加え且つ
周波数エンコードを加え、これによりｋ−空間のｐビュ
ーのイメージング用データＤ’（ｆ，ｐ）を順に収集す
る。そして、各ビューのナビゲーション用データＳＲｐ
（ｆ’）の相互関係から位相補正量Δφｘ（ｐ）を算出
し、その位相補正量Δφｘ（ｐ）で各ビューのイメージ
ング用データＤ’（ｆ，ｐ）を位相補正する。上記ナビ
ゲーション用データＳＲｐ（ｆ’）は、常に位相エンコ
ード量“０”で収集しているので、その相互関係から各
ビューでの信号源のｘ軸方向の移動量Δｘｐを求めるこ
とが出来る。ここで、 ２πｆΔｘｐ／Ｆｗ＝Δφｘ（ｐ） とすれば、前記（４）式となり、位相補正して得られた
Ｄ（ｆ，ｐ）に対して２次元フーリエ変換を行えば、信
号源ｇ（ｘ，ｙ）が移動していない場合に相当するイメ
ージが得られる。すなわち、ｘ軸方向に移動する被検体
を撮像し、且つ、アーチファクトを抑制して、イメージ
ングすることが出来る。
【００１９】上記第２の観点によるＭＲイメージング方
法および上記第４の観点によるＭＲＩ装置では、１つの
ＲＦパルスに対して少なくとも２つのエコーを発生させ
る。その内の１つのエコーに対しては、周波数エンコー
ド量を“０”とし、位相軸にサンプリング数Ｍ’の読み
出し勾配を印加しながら、データ番号ｐ’（−Ｍ’／２
≦ｐ’≦Ｍ’／２−１）のｐビューにおける位相軸方向
についてのナビゲーション用データＳＷｐ（ｐ’）を収
集する。また、残りのエコーに対しては、ｋ−空間のｐ
ビューに応じた位相エンコードを加え且つ周波数エンコ
ードを加え、これによりｋ−空間のｐビューのイメージ
ング用データＤ’（ｆ，ｐ）を順に収集する。そして、
各ビューのナビゲーション用データＳＷｐ（ｐ’）の相
互関係から位相補正量Δφｙ（ｐ）を算出し、その位相
補正量Δφｙ（ｐ）で各ビューのイメージング用データ
Ｄ’（ｆ，ｐ）を位相補正する。上記ナビゲーション用
データＳＷｐ（ｐ’）は、常に周波数エンコード量
“０”で収集しているので、その相互関係から各ビュー
での信号源のｙ軸方向の移動量Δｙｐを求めることが出
来る。ここで、 ２πｐΔｙｐ／Ｐｗ＝Δφｙ（ｐ） とすれば、前記（５）式となり、位相補正して得られた
Ｄ（ｆ，ｐ）に対して２次元フーリエ変換を行えば、信
号源ｇ（ｘ，ｙ）が移動していない場合に相当するイメ
ージが得られる。すなわち、ｙ軸方向に移動する被検体
を撮像し、且つ、アーチファクトを抑制して、イメージ
ングすることが出来る。
法および上記第４の観点によるＭＲＩ装置では、１つの
ＲＦパルスに対して少なくとも２つのエコーを発生させ
る。その内の１つのエコーに対しては、周波数エンコー
ド量を“０”とし、位相軸にサンプリング数Ｍ’の読み
出し勾配を印加しながら、データ番号ｐ’（−Ｍ’／２
≦ｐ’≦Ｍ’／２−１）のｐビューにおける位相軸方向
についてのナビゲーション用データＳＷｐ（ｐ’）を収
集する。また、残りのエコーに対しては、ｋ−空間のｐ
ビューに応じた位相エンコードを加え且つ周波数エンコ
ードを加え、これによりｋ−空間のｐビューのイメージ
ング用データＤ’（ｆ，ｐ）を順に収集する。そして、
各ビューのナビゲーション用データＳＷｐ（ｐ’）の相
互関係から位相補正量Δφｙ（ｐ）を算出し、その位相
補正量Δφｙ（ｐ）で各ビューのイメージング用データ
Ｄ’（ｆ，ｐ）を位相補正する。上記ナビゲーション用
データＳＷｐ（ｐ’）は、常に周波数エンコード量
“０”で収集しているので、その相互関係から各ビュー
での信号源のｙ軸方向の移動量Δｙｐを求めることが出
来る。ここで、 ２πｐΔｙｐ／Ｐｗ＝Δφｙ（ｐ） とすれば、前記（５）式となり、位相補正して得られた
Ｄ（ｆ，ｐ）に対して２次元フーリエ変換を行えば、信
号源ｇ（ｘ，ｙ）が移動していない場合に相当するイメ
ージが得られる。すなわち、ｙ軸方向に移動する被検体
を撮像し、且つ、アーチファクトを抑制して、イメージ
ングすることが出来る。
【００２０】なお、上記構成の組み合わせの構成とし
て、１つのＲＦパルスに対して少なくとも３つのエコー
を発生させる。その内の１つのエコーに対しては、位相
エンコード量を“０”とし、周波数軸にサンプリング数
Ｎ’の読み出し勾配を印加しながら、データ番号ｆ’
（−Ｎ’／２≦ｆ’≦Ｎ’／２−１）のｐビューにおけ
る周波数軸方向についてのナビゲーション用データＳＲ
ｐ（ｆ’）を収集する。また、他の１つのエコーに対し
ては、周波数エンコード量を“０”とし、位相軸にサン
プリング数Ｍ’の読み出し勾配を印加しながら、データ
番号ｐ’（−Ｍ’／２≦ｐ’≦Ｍ’／２−１）のｐビュ
ーにおける位相軸方向についてのナビゲーション用デー
タＳＷｐ（ｐ’）を収集する。さらに、残りのエコーに
対しては、ｋ−空間のｐビューに応じた位相エンコード
を加え且つ周波数エンコードを加え、これによりｋ−空
間のｐビューのイメージング用データＤ’（ｆ，ｐ）を
順に収集する。そして、各ビューのナビゲーション用デ
ータＳＲｐ（ｆ’）の相互関係から位相補正量Δφｘ
（ｐ）を算出すると共に各ビューのナビゲーション用デ
ータＳＷｐ（ｐ’）の相互関係から位相補正量Δφｙ
（ｐ）を算出し、その位相補正量Δφｘ（ｐ），Δφｙ
（ｐ）で各ビューのイメージング用データＤ’（ｆ，
ｐ）を位相補正する。上記ナビゲーション用データＳＲ
ｐ（ｆ’）は、常に位相エンコード量“０”で収集して
いるので、その相互関係から各ビューでの信号源のｘ軸
方向の移動量Δｘｐを求めることが出来る。また、上記
ナビゲーション用データＳＷｐ（ｐ’）は、常に周波数
エンコード量“０”で収集しているので、その相互関係
により各ビューでの信号源のｙ軸方向の移動量Δｙｐを
求めることが出来る。ここで、 ２πｆΔｘｐ／Ｆｗ＝Δφｘ（ｐ） ２πｐΔｙｐ／Ｐｗ＝Δφｙ（ｐ） とすれば、前記（６）式となり、位相補正して得られた
Ｄ（ｆ，ｐ）に対して２次元フーリエ変換を行えば、信
号源ｇ（ｘ，ｙ）が移動していない場合に相当するイメ
ージが得られる。すなわち、ｘ軸方向およびｙ軸方向に
移動する被検体を撮像し、且つ、アーチファクトを抑制
して、イメージングすることが出来る。
て、１つのＲＦパルスに対して少なくとも３つのエコー
を発生させる。その内の１つのエコーに対しては、位相
エンコード量を“０”とし、周波数軸にサンプリング数
Ｎ’の読み出し勾配を印加しながら、データ番号ｆ’
（−Ｎ’／２≦ｆ’≦Ｎ’／２−１）のｐビューにおけ
る周波数軸方向についてのナビゲーション用データＳＲ
ｐ（ｆ’）を収集する。また、他の１つのエコーに対し
ては、周波数エンコード量を“０”とし、位相軸にサン
プリング数Ｍ’の読み出し勾配を印加しながら、データ
番号ｐ’（−Ｍ’／２≦ｐ’≦Ｍ’／２−１）のｐビュ
ーにおける位相軸方向についてのナビゲーション用デー
タＳＷｐ（ｐ’）を収集する。さらに、残りのエコーに
対しては、ｋ−空間のｐビューに応じた位相エンコード
を加え且つ周波数エンコードを加え、これによりｋ−空
間のｐビューのイメージング用データＤ’（ｆ，ｐ）を
順に収集する。そして、各ビューのナビゲーション用デ
ータＳＲｐ（ｆ’）の相互関係から位相補正量Δφｘ
（ｐ）を算出すると共に各ビューのナビゲーション用デ
ータＳＷｐ（ｐ’）の相互関係から位相補正量Δφｙ
（ｐ）を算出し、その位相補正量Δφｘ（ｐ），Δφｙ
（ｐ）で各ビューのイメージング用データＤ’（ｆ，
ｐ）を位相補正する。上記ナビゲーション用データＳＲ
ｐ（ｆ’）は、常に位相エンコード量“０”で収集して
いるので、その相互関係から各ビューでの信号源のｘ軸
方向の移動量Δｘｐを求めることが出来る。また、上記
ナビゲーション用データＳＷｐ（ｐ’）は、常に周波数
エンコード量“０”で収集しているので、その相互関係
により各ビューでの信号源のｙ軸方向の移動量Δｙｐを
求めることが出来る。ここで、 ２πｆΔｘｐ／Ｆｗ＝Δφｘ（ｐ） ２πｐΔｙｐ／Ｐｗ＝Δφｙ（ｐ） とすれば、前記（６）式となり、位相補正して得られた
Ｄ（ｆ，ｐ）に対して２次元フーリエ変換を行えば、信
号源ｇ（ｘ，ｙ）が移動していない場合に相当するイメ
ージが得られる。すなわち、ｘ軸方向およびｙ軸方向に
移動する被検体を撮像し、且つ、アーチファクトを抑制
して、イメージングすることが出来る。
【００２１】上記第５の観点によるＭＲＩ装置では、ｋ
−空間の第１のエンコード軸のエンコード量を“０”と
し、ｋ−空間の第２のエンコード軸に読み出し勾配を印
加しながらナビゲーション用データＳ１を収集する。そ
の後、同様にしてナビゲーション用データＳ２を収集す
る。そして、前記ナビゲーション用データＳ１の１次元
フーリエ変換の絶対値から求まる投影プロファイルと前
記ナビゲーション用データＳ２の１次元フーリエ変換の
絶対値から求まる投影プロファイルの相互相関から前記
第２のエンコード軸に対応する実空間の方向についての
信号源の移動量Δを算出する。上記投影プロファイル
は、実空間における信号源の位置を反映しているから、
それらの相互相関関数の最大値を与える変数値は位置の
ずれを表している。従って、これより信号源の移動量Δ
が判る。
−空間の第１のエンコード軸のエンコード量を“０”と
し、ｋ−空間の第２のエンコード軸に読み出し勾配を印
加しながらナビゲーション用データＳ１を収集する。そ
の後、同様にしてナビゲーション用データＳ２を収集す
る。そして、前記ナビゲーション用データＳ１の１次元
フーリエ変換の絶対値から求まる投影プロファイルと前
記ナビゲーション用データＳ２の１次元フーリエ変換の
絶対値から求まる投影プロファイルの相互相関から前記
第２のエンコード軸に対応する実空間の方向についての
信号源の移動量Δを算出する。上記投影プロファイル
は、実空間における信号源の位置を反映しているから、
それらの相互相関関数の最大値を与える変数値は位置の
ずれを表している。従って、これより信号源の移動量Δ
が判る。
【００２２】
【実施例】以下、図に示す実施例に基づいてこの発明を
さらに詳細に説明する。なお、これによりこの発明が限
定されるものではない。
さらに詳細に説明する。なお、これによりこの発明が限
定されるものではない。
【００２３】−第１実施例− 第１実施例は、ｋ−空間の周波数軸方向に対応する実空
間の軸をｘ軸とし、ｋ−空間の位相軸方向に対応する実
空間の軸をｙ軸とし、ｘ軸およびｙ軸に直交する実空間
の軸をｚ軸としたときに、被検体の移動方向がｘ軸方向
である場合に対応した実施例である。
間の軸をｘ軸とし、ｋ−空間の位相軸方向に対応する実
空間の軸をｙ軸とし、ｘ軸およびｙ軸に直交する実空間
の軸をｚ軸としたときに、被検体の移動方向がｘ軸方向
である場合に対応した実施例である。
【００２４】図１は、この発明のＭＲＩ装置の第１実施
例のブロック図である。このＭＲＩ装置１００におい
て、マグネットアセンブリ１は、被検体（図示省略）を
内部に挿入するためのボア（孔）を有し、このボアを取
りまくようにして、被検体に一定の静磁場を印加する静
磁場コイルと、勾配磁場を発生するための勾配磁場コイ
ル（勾配磁場コイルは、直交するＺ軸，Ｙ軸，Ｘ軸のコ
イルを備えている）と、被検体内の原子核のスピンを励
起するためのＲＦパルスを与える送信コイルと、被検体
からのＮＭＲ（Nuclear Magnetic Resonance）信号を検
出する受信コイル等が配置されている。静磁場コイル，
勾配磁場コイル，送信コイルおよび受信コイルは、それ
ぞれ主磁場電源２，勾配磁場駆動回路３，ＲＦ電力増幅
器４および前置増幅器５に接続されている。
例のブロック図である。このＭＲＩ装置１００におい
て、マグネットアセンブリ１は、被検体（図示省略）を
内部に挿入するためのボア（孔）を有し、このボアを取
りまくようにして、被検体に一定の静磁場を印加する静
磁場コイルと、勾配磁場を発生するための勾配磁場コイ
ル（勾配磁場コイルは、直交するＺ軸，Ｙ軸，Ｘ軸のコ
イルを備えている）と、被検体内の原子核のスピンを励
起するためのＲＦパルスを与える送信コイルと、被検体
からのＮＭＲ（Nuclear Magnetic Resonance）信号を検
出する受信コイル等が配置されている。静磁場コイル，
勾配磁場コイル，送信コイルおよび受信コイルは、それ
ぞれ主磁場電源２，勾配磁場駆動回路３，ＲＦ電力増幅
器４および前置増幅器５に接続されている。
【００２５】シーケンス記憶回路８は、計算機７からの
指令に従い、記憶されているパルスシーケンスに基づい
て勾配磁場駆動回路３を操作し、前記マグネットアセン
ブリ１の勾配磁場コイルから勾配磁場を発生させると共
に、ゲート変調回路９を操作し、ＲＦ発振回路１０の搬
送波出力信号を所定タイミング・所定包絡線形状のパル
ス状信号に変調し、それをＲＦパルスとしてＲＦ電力増
幅器４に加え、ＲＦ電力増幅器４でパワー増幅した後、
前記マグネットアセンブリ１の送信コイルに印加し、被
検体の目的領域を励起する。
指令に従い、記憶されているパルスシーケンスに基づい
て勾配磁場駆動回路３を操作し、前記マグネットアセン
ブリ１の勾配磁場コイルから勾配磁場を発生させると共
に、ゲート変調回路９を操作し、ＲＦ発振回路１０の搬
送波出力信号を所定タイミング・所定包絡線形状のパル
ス状信号に変調し、それをＲＦパルスとしてＲＦ電力増
幅器４に加え、ＲＦ電力増幅器４でパワー増幅した後、
前記マグネットアセンブリ１の送信コイルに印加し、被
検体の目的領域を励起する。
【００２６】前置増幅器５は、マグネットアセンブリ１
の受信コイルで検出された被検体からのＮＭＲ信号を増
幅し、位相検波器１２に入力する。位相検波器１２は、
ＲＦ発振回路１０の搬送波出力信号を参照信号とし、前
置増幅器５からのＮＭＲ信号を位相検波して、Ａ／Ｄ変
換器１１に与える。Ａ／Ｄ変換器１１は、位相検波後の
アナログ信号をディジタル信号に変換して、計算機７に
入力する。計算機７は、Ａ／Ｄ変換器１１からのデジタ
ル信号に対する画像再構成演算を行い、被検体の目的領
域のイメージを生成する。このイメージは、表示装置６
にて表示される。また、計算機７は、操作卓１３から入
力された情報を受け取るなどの全体的な制御を受け持
つ。
の受信コイルで検出された被検体からのＮＭＲ信号を増
幅し、位相検波器１２に入力する。位相検波器１２は、
ＲＦ発振回路１０の搬送波出力信号を参照信号とし、前
置増幅器５からのＮＭＲ信号を位相検波して、Ａ／Ｄ変
換器１１に与える。Ａ／Ｄ変換器１１は、位相検波後の
アナログ信号をディジタル信号に変換して、計算機７に
入力する。計算機７は、Ａ／Ｄ変換器１１からのデジタ
ル信号に対する画像再構成演算を行い、被検体の目的領
域のイメージを生成する。このイメージは、表示装置６
にて表示される。また、計算機７は、操作卓１３から入
力された情報を受け取るなどの全体的な制御を受け持
つ。
【００２７】移動装置１４は、被検体を乗せて、前記マ
グネットアセンブリ１のボア内に被検体を搬入する。
グネットアセンブリ１のボア内に被検体を搬入する。
【００２８】前記シーケンス記憶回路８は、この発明の
ＭＲイメージング方法を実施するためのパルスシーケン
スを記憶している。また、前記計算機７は、この発明の
ＭＲイメージング方法を実施するためのパルスシーケン
スに基づいてイメージング用データおよびナビゲーショ
ン用データを収集すると共に、ナビゲーション用データ
に基づいてイメージング用データを位相補正し、位相補
正後のイメージング用データに対する画像再構成演算を
行い、被検体の目的領域のイメージを生成する。
ＭＲイメージング方法を実施するためのパルスシーケン
スを記憶している。また、前記計算機７は、この発明の
ＭＲイメージング方法を実施するためのパルスシーケン
スに基づいてイメージング用データおよびナビゲーショ
ン用データを収集すると共に、ナビゲーション用データ
に基づいてイメージング用データを位相補正し、位相補
正後のイメージング用データに対する画像再構成演算を
行い、被検体の目的領域のイメージを生成する。
【００２９】図２は、被検体をｘ軸方向に移動しながら
撮像する場合のパルスシーケンスである。このパルスシ
ーケンスＡでは、フリップ角α゜のＲＦパルスＲαを印
加すると共にｚ軸にスライス勾配ｓｓを印加する。次
に、ｙ軸に位相エンコード勾配ｐｅを印加する。次に、
ｘ軸に周波数エンコード勾配ｒｉを印加しながら、前記
ＲＦパルスＲαから時間ＴＥ１後に発生するイメージエ
コーＩＥをサンプリングして、ｐビューのイメージング
用データＤ’（ｆ，ｐ）を収集する。次に、リワインド
（Rewind）勾配ｐｒを印加して位相エンコード量を
“０”に戻す。次に、前記周波数エンコード勾配ｒｉと
極性を反転したサンプリング数Ｎ’の読み出し勾配ｒｎ
１を印加しながら、前記ＲＦパルスＲαから時間ＴＥ２
後に発生するナビゲーションエコーＮＥ１をサンプリン
グして、ｐビューのイメージング用データＤ’（ｆ，
ｐ）に対応するナビゲーション用データＳＲｐ（ｆ’）
を収集する。これを位相エンコード勾配ｐｅおよびリワ
インド勾配ｐｒを変えながら、繰り返し時間ＴＲでＭ回
繰り返し、Ｍビュー分のイメージング用データＤ’
（ｆ，−Ｍ／２）〜Ｄ’（ｆ，Ｍ／２−１）およびそれ
らに対応するＭビュー分のナビゲーション用データＳＲ
-M/2（ｆ’）〜ＳＲM/2-1（ｆ’）を収集する。なお、
図２で、ｄｐはディフェーズ（Dephase）勾配，ｃｒは
クラッシュ（Crush）勾配である。
撮像する場合のパルスシーケンスである。このパルスシ
ーケンスＡでは、フリップ角α゜のＲＦパルスＲαを印
加すると共にｚ軸にスライス勾配ｓｓを印加する。次
に、ｙ軸に位相エンコード勾配ｐｅを印加する。次に、
ｘ軸に周波数エンコード勾配ｒｉを印加しながら、前記
ＲＦパルスＲαから時間ＴＥ１後に発生するイメージエ
コーＩＥをサンプリングして、ｐビューのイメージング
用データＤ’（ｆ，ｐ）を収集する。次に、リワインド
（Rewind）勾配ｐｒを印加して位相エンコード量を
“０”に戻す。次に、前記周波数エンコード勾配ｒｉと
極性を反転したサンプリング数Ｎ’の読み出し勾配ｒｎ
１を印加しながら、前記ＲＦパルスＲαから時間ＴＥ２
後に発生するナビゲーションエコーＮＥ１をサンプリン
グして、ｐビューのイメージング用データＤ’（ｆ，
ｐ）に対応するナビゲーション用データＳＲｐ（ｆ’）
を収集する。これを位相エンコード勾配ｐｅおよびリワ
インド勾配ｐｒを変えながら、繰り返し時間ＴＲでＭ回
繰り返し、Ｍビュー分のイメージング用データＤ’
（ｆ，−Ｍ／２）〜Ｄ’（ｆ，Ｍ／２−１）およびそれ
らに対応するＭビュー分のナビゲーション用データＳＲ
-M/2（ｆ’）〜ＳＲM/2-1（ｆ’）を収集する。なお、
図２で、ｄｐはディフェーズ（Dephase）勾配，ｃｒは
クラッシュ（Crush）勾配である。
【００３０】図３の（ａ）は、ｐ＝−Ｍ／２のときのｋ
−空間でのイメージング用データＤ’（ｆ，−Ｍ／２）
と、ナビゲーション用データＳＲ-M/2（ｆ’）の関係を
示す概念図である。図３の（ｂ）は、ｐ＝−Ｍ／２＋１
のときのｋ−空間でのイメージング用データＤ’（ｆ，
−Ｍ／２＋１）と、ナビゲーション用データＳＲ-M/2+1
（ｆ’）の関係を示す概念図である。
−空間でのイメージング用データＤ’（ｆ，−Ｍ／２）
と、ナビゲーション用データＳＲ-M/2（ｆ’）の関係を
示す概念図である。図３の（ｂ）は、ｐ＝−Ｍ／２＋１
のときのｋ−空間でのイメージング用データＤ’（ｆ，
−Ｍ／２＋１）と、ナビゲーション用データＳＲ-M/2+1
（ｆ’）の関係を示す概念図である。
【００３１】図４は、上記ナビゲーション用データＳＲ
ｐ（ｆ’）に基づいてイメージング用データＤ’（ｆ，
ｐ）を位相補正する処理のフローチャートである。ステ
ップＶ１では、ビュー番号カウンタｐ＝−Ｍ／２に初期
化する。ステップＶ２では、０ビューに対応するナビゲ
ーション用データＳＲｏ（ｆ’）と、ｐビューでのナビ
ゲーション用データＳＲｐ（ｆ’）の共役ＳＲｐ
*（ｆ'）の乗算を行う。すなわち、 ＺＲｐ（ｆ’）＝ＳＲｐ*（ｆ’）・ＳＲｏ（ｆ’） を求める。
ｐ（ｆ’）に基づいてイメージング用データＤ’（ｆ，
ｐ）を位相補正する処理のフローチャートである。ステ
ップＶ１では、ビュー番号カウンタｐ＝−Ｍ／２に初期
化する。ステップＶ２では、０ビューに対応するナビゲ
ーション用データＳＲｏ（ｆ’）と、ｐビューでのナビ
ゲーション用データＳＲｐ（ｆ’）の共役ＳＲｐ
*（ｆ'）の乗算を行う。すなわち、 ＺＲｐ（ｆ’）＝ＳＲｐ*（ｆ’）・ＳＲｏ（ｆ’） を求める。
【００３２】ステップＶ３では、積ＺＲｐ（ｆ’）を１
次元フーリエ変換し、相互相関関数ｚｒｐ（ｘ）を求め
る。図５に、相互相関関数を示す。相互相関関数ｚｒｐ
（ｘ）は、ナビゲーション用データＳＲｏ（ｆ’）の１
次元フーリエ変換の絶対値から求まる投影プロファイル
ｓｒｏ（ｘ）と、ナビゲーション用データＳＲｐ
（ｆ’）の１次元フーリエ変換の絶対値から求まる投影
プロファイルｓｒｐ（ｘ）の相互相関関数になっている
（相関定理）。
次元フーリエ変換し、相互相関関数ｚｒｐ（ｘ）を求め
る。図５に、相互相関関数を示す。相互相関関数ｚｒｐ
（ｘ）は、ナビゲーション用データＳＲｏ（ｆ’）の１
次元フーリエ変換の絶対値から求まる投影プロファイル
ｓｒｏ（ｘ）と、ナビゲーション用データＳＲｐ
（ｆ’）の１次元フーリエ変換の絶対値から求まる投影
プロファイルｓｒｐ（ｘ）の相互相関関数になっている
（相関定理）。
【００３３】図４に戻り、ステップＶ４では、最小二乗
法を用いた補間近似により、相互相関関数ｚｒｐ（ｘ）
の最大値を与える変数値Δｘｐを求める。ステップＶ５
では、２πｆΔｘｐ／Ｆｗ＝Δφｘ（ｐ）により位相補
正量Δφｘ（ｐ）を算出し、イメージング用データＤ’
（ｆ，ｐ）に対して、 exp｛ｊΔφｘ（ｐ）｝ を乗算し、これをｐビューの位相補正後のイメージング
用データＤ（ｆ，ｐ）とする。ステップＶ６，Ｖ７で
は、ｐ＝−Ｍ／２＋１からｐ＝Ｍ／２−１まで上記ステ
ップＶ２〜Ｖ５を反復実行し、ｋ空間を構成する全ビュ
ーのイメージング用データＤ（ｆ，ｐ）を取得する。
法を用いた補間近似により、相互相関関数ｚｒｐ（ｘ）
の最大値を与える変数値Δｘｐを求める。ステップＶ５
では、２πｆΔｘｐ／Ｆｗ＝Δφｘ（ｐ）により位相補
正量Δφｘ（ｐ）を算出し、イメージング用データＤ’
（ｆ，ｐ）に対して、 exp｛ｊΔφｘ（ｐ）｝ を乗算し、これをｐビューの位相補正後のイメージング
用データＤ（ｆ，ｐ）とする。ステップＶ６，Ｖ７で
は、ｐ＝−Ｍ／２＋１からｐ＝Ｍ／２−１まで上記ステ
ップＶ２〜Ｖ５を反復実行し、ｋ空間を構成する全ビュ
ーのイメージング用データＤ（ｆ，ｐ）を取得する。
【００３４】上記位相補正したイメージング用データＤ
（ｆ，ｐ）から画像再構成し、イメージを生成すれば、
被検体をｘ軸方向に連続的に移動しながら撮像しても、
モーションアーチファクトを低減した高画質のイメージ
が得られる。
（ｆ，ｐ）から画像再構成し、イメージを生成すれば、
被検体をｘ軸方向に連続的に移動しながら撮像しても、
モーションアーチファクトを低減した高画質のイメージ
が得られる。
【００３５】−第２実施例− 第２実施例は、ｋ−空間の周波数軸方向に対応する実空
間の軸をｘ軸とし、ｋ−空間の位相軸方向に対応する実
空間の軸をｙ軸とし、ｘ軸およびｙ軸に直交する実空間
の軸をｚ軸としたときに、被検体の移動方向がｙ軸方向
である場合に対応した実施例である。ＭＲＩ装置のブロ
ック図は、図１と同様である。
間の軸をｘ軸とし、ｋ−空間の位相軸方向に対応する実
空間の軸をｙ軸とし、ｘ軸およびｙ軸に直交する実空間
の軸をｚ軸としたときに、被検体の移動方向がｙ軸方向
である場合に対応した実施例である。ＭＲＩ装置のブロ
ック図は、図１と同様である。
【００３６】図６は、被検体をｙ軸方向に移動しながら
撮像する場合のパルスシーケンスである。このパルスシ
ーケンスＢでは、フリップ角α゜のＲＦパルスＲαを印
加すると共にｚ軸にスライス勾配ｓｓを印加する。次
に、ｙ軸に位相エンコード勾配ｐｅを印加する。次に、
ｘ軸に周波数エンコード勾配ｒｉを印加しながら、前記
ＲＦパルスＲαから時間ＴＥ１後に発生するイメージエ
コーＩＥをサンプリングして、ｐビューのイメージング
用データＤ’（ｆ，ｐ）を収集する。次に、リワインド
勾配ｐｒを印加して位相エンコード量を“０”に戻す。
また、前記周波数エンコード勾配ｒｉと極性を反転した
リフェーズ（Ｒephase）勾配ｒｐを印加して周波数エン
コード量を“０”に戻す。次に、位相軸（ｙ軸）にディ
フェーズ勾配ｄｐおよびサンプリング数Ｍ’の読み出し
勾配ｒｎ２を印加しながら、前記ＲＦパルスＲαから時
間ＴＥ３後に発生するナビゲーションエコーＮＥ２をサ
ンプリングして、ｐビューのイメージング用データＤ’
（ｆ，ｐ）に対応するナビゲーション用データＳＷｐ
（ｐ’）を収集する。これを位相エンコード勾配ｐｅお
よびリワインド勾配ｐｒを変えながら、繰り返し時間Ｔ
ＲでＭ回繰り返し、Ｍビュー分のイメージング用データ
Ｄ’（ｆ，−Ｍ／２）〜Ｄ’（ｆ，Ｍ／２−１）および
それらに対応するＭビュー分のナビゲーション用データ
ＳＷ-M/2（ｆ’）〜ＳＷM/2-1（ｆ’）を収集する。図
７の（ａ）は、ｐ＝−Ｍ／２のときのｋ−空間でのイメ
ージング用データＤ’（ｆ，−Ｍ／２）と、ナビゲーシ
ョン用データＳＷ-M/2（ｆ’）の関係を示す概念図であ
る。図７の（ｂ）は、ｐ＝−Ｍ／２＋１のときのｋ−空
間でのイメージング用データＤ’（ｆ，−Ｍ／２＋１）
と、ナビゲーション用データＳＷ-M/2+1（ｆ’）の関係
を示す概念図である。
撮像する場合のパルスシーケンスである。このパルスシ
ーケンスＢでは、フリップ角α゜のＲＦパルスＲαを印
加すると共にｚ軸にスライス勾配ｓｓを印加する。次
に、ｙ軸に位相エンコード勾配ｐｅを印加する。次に、
ｘ軸に周波数エンコード勾配ｒｉを印加しながら、前記
ＲＦパルスＲαから時間ＴＥ１後に発生するイメージエ
コーＩＥをサンプリングして、ｐビューのイメージング
用データＤ’（ｆ，ｐ）を収集する。次に、リワインド
勾配ｐｒを印加して位相エンコード量を“０”に戻す。
また、前記周波数エンコード勾配ｒｉと極性を反転した
リフェーズ（Ｒephase）勾配ｒｐを印加して周波数エン
コード量を“０”に戻す。次に、位相軸（ｙ軸）にディ
フェーズ勾配ｄｐおよびサンプリング数Ｍ’の読み出し
勾配ｒｎ２を印加しながら、前記ＲＦパルスＲαから時
間ＴＥ３後に発生するナビゲーションエコーＮＥ２をサ
ンプリングして、ｐビューのイメージング用データＤ’
（ｆ，ｐ）に対応するナビゲーション用データＳＷｐ
（ｐ’）を収集する。これを位相エンコード勾配ｐｅお
よびリワインド勾配ｐｒを変えながら、繰り返し時間Ｔ
ＲでＭ回繰り返し、Ｍビュー分のイメージング用データ
Ｄ’（ｆ，−Ｍ／２）〜Ｄ’（ｆ，Ｍ／２−１）および
それらに対応するＭビュー分のナビゲーション用データ
ＳＷ-M/2（ｆ’）〜ＳＷM/2-1（ｆ’）を収集する。図
７の（ａ）は、ｐ＝−Ｍ／２のときのｋ−空間でのイメ
ージング用データＤ’（ｆ，−Ｍ／２）と、ナビゲーシ
ョン用データＳＷ-M/2（ｆ’）の関係を示す概念図であ
る。図７の（ｂ）は、ｐ＝−Ｍ／２＋１のときのｋ−空
間でのイメージング用データＤ’（ｆ，−Ｍ／２＋１）
と、ナビゲーション用データＳＷ-M/2+1（ｆ’）の関係
を示す概念図である。
【００３７】図８は、上記ナビゲーション用データＳＷ
ｐ（ｐ’）に基づいてイメージング用データＤ’（ｆ，
ｐ）を位相補正する処理のフローチャートである。ステ
ップＵ１では、ビュー番号カウンタｐ＝−Ｍ／２に初期
化する。ステップＵ２では、０ビューに対応するナビゲ
ーション用データＳＷｏ（ｐ’）と、ｐビューでのナビ
ゲーション用データＳＷｐ（ｐ’）の共役ＳＷｐ
*（ｐ'）の乗算を行う。すなわち、 ＺＷｐ（ｐ’）＝ＳＷｐ*（ｐ’）・ＳＷｏ（ｐ’） を求める。
ｐ（ｐ’）に基づいてイメージング用データＤ’（ｆ，
ｐ）を位相補正する処理のフローチャートである。ステ
ップＵ１では、ビュー番号カウンタｐ＝−Ｍ／２に初期
化する。ステップＵ２では、０ビューに対応するナビゲ
ーション用データＳＷｏ（ｐ’）と、ｐビューでのナビ
ゲーション用データＳＷｐ（ｐ’）の共役ＳＷｐ
*（ｐ'）の乗算を行う。すなわち、 ＺＷｐ（ｐ’）＝ＳＷｐ*（ｐ’）・ＳＷｏ（ｐ’） を求める。
【００３８】ステップＵ３では、積ＺＷｐ（ｐ’）を１
次元フーリエ変換し、相互相関関数ｚｗｐ（ｙ）を求め
る。図９に、相互相関関数ｚｗｐ（ｙ）を示す。相互相
関関数ｚｗｐ（ｙ）は、ナビゲーション用データＳＷｏ
（ｐ’）の１次元フーリエ変換の絶対値から求まる投影
プロファイルｓｗｏ（ｙ）と、ナビゲーション用データ
ＳＷｐ（ｐ’）の１次元フーリエ変換の絶対値から求ま
る投影プロファイルｓｒｐ（ｙ）の相互相関関数になっ
ている（相関定理）。
次元フーリエ変換し、相互相関関数ｚｗｐ（ｙ）を求め
る。図９に、相互相関関数ｚｗｐ（ｙ）を示す。相互相
関関数ｚｗｐ（ｙ）は、ナビゲーション用データＳＷｏ
（ｐ’）の１次元フーリエ変換の絶対値から求まる投影
プロファイルｓｗｏ（ｙ）と、ナビゲーション用データ
ＳＷｐ（ｐ’）の１次元フーリエ変換の絶対値から求ま
る投影プロファイルｓｒｐ（ｙ）の相互相関関数になっ
ている（相関定理）。
【００３９】図８に戻り、ステップＵ４では、最小二乗
法を用いた補間近似により、相互相関関数ｚｗｐ（ｙ）
の最大値を与える変数値Δｙｐを求める。ステップＵ５
では、２πｐΔｙｐ／Ｐｗ＝Δφｙ（ｐ）により位相補
正量Δφｙ（ｐ）を算出し、イメージング用データＤ’
（ｆ，ｐ）に対して、 exp｛ｊΔφｙ（ｐ）｝ を乗算し、これをｐビューの位相補正後のイメージング
用データＤ（ｆ，ｐ）とする。ステップＵ６，Ｕ７で
は、ｐ＝−Ｍ／２＋１からｐ＝Ｍ／２−１まで上記ステ
ップＵ２〜Ｕ５を反復実行し、ｋ空間を構成する全ビュ
ーのイメージング用データＤ（ｆ，ｐ）を取得する。
法を用いた補間近似により、相互相関関数ｚｗｐ（ｙ）
の最大値を与える変数値Δｙｐを求める。ステップＵ５
では、２πｐΔｙｐ／Ｐｗ＝Δφｙ（ｐ）により位相補
正量Δφｙ（ｐ）を算出し、イメージング用データＤ’
（ｆ，ｐ）に対して、 exp｛ｊΔφｙ（ｐ）｝ を乗算し、これをｐビューの位相補正後のイメージング
用データＤ（ｆ，ｐ）とする。ステップＵ６，Ｕ７で
は、ｐ＝−Ｍ／２＋１からｐ＝Ｍ／２−１まで上記ステ
ップＵ２〜Ｕ５を反復実行し、ｋ空間を構成する全ビュ
ーのイメージング用データＤ（ｆ，ｐ）を取得する。
【００４０】上記位相補正したイメージング用データＤ
（ｆ，ｐ）から画像再構成し、イメージを生成すれば、
被検体をｙ軸方向に連続的に移動しながら撮像しても、
モーションアーチファクトを低減した高画質のイメージ
が得られる。
（ｆ，ｐ）から画像再構成し、イメージを生成すれば、
被検体をｙ軸方向に連続的に移動しながら撮像しても、
モーションアーチファクトを低減した高画質のイメージ
が得られる。
【００４１】−第３実施例− 第３実施例は、ｋ−空間の周波数軸方向に対応する実空
間の軸をｘ軸とし、ｋ−空間の位相軸方向に対応する実
空間の軸をｙ軸とし、ｘ軸およびｙ軸に直交する実空間
の軸をｚ軸としたときに、被検体の移動方向がｘ軸およ
びｙ軸方向である場合に対応した実施例である。ＭＲＩ
装置のブロック図は、図１と同様である。
間の軸をｘ軸とし、ｋ−空間の位相軸方向に対応する実
空間の軸をｙ軸とし、ｘ軸およびｙ軸に直交する実空間
の軸をｚ軸としたときに、被検体の移動方向がｘ軸およ
びｙ軸方向である場合に対応した実施例である。ＭＲＩ
装置のブロック図は、図１と同様である。
【００４２】図１０は、被検体をｘ軸方向およびｙ軸方
向に移動しながら撮像する場合のパルスシーケンスであ
る。このパルスシーケンスＣでは、フリップ角α゜のＲ
ＦパルスＲαを印加すると共にｚ軸にスライス勾配ｓｓ
を印加する。次に、ｙ軸に位相エンコード勾配ｐｅを印
加する。次に、ｘ軸に周波数エンコード勾配ｒｉを印加
しながら、前記ＲＦパルスＲαから時間ＴＥ１後に発生
するイメージエコーＩＥをサンプリングして、ｐビュー
のイメージング用データＤ’（ｆ，ｐ）を収集する。次
に、リワインド勾配ｐｒを印加して位相エンコード量を
“０”に戻す。次に、前記周波数エンコード勾配ｒｉと
極性を反転したサンプリング数Ｎ’の読み出し勾配ｒｎ
１を印加しながら、前記ＲＦパルスＲαから時間ＴＥ２
後に発生するナビゲーションエコーＮＥ１をサンプリン
グして、ｐビューのイメージング用データＤ’（ｆ，
ｐ）に対応するナビゲーション用データＳＲｐ（ｆ’）
を収集する。次に、前記周波数エンコード勾配ｒｎ１と
極性を反転したリフェーズ勾配ｒｐを印加して周波数エ
ンコード量を“０”に戻す。また、位相軸（ｙ軸）にデ
ィフェーズ勾配ｄｐおよびサンプリング数Ｍ’の読み出
し勾配ｒｎ２を印加しながら、前記ＲＦパルスＲαから
時間ＴＥ３後に発生するナビゲーションエコーＮＥ２を
サンプリングして、ｐビューのイメージング用データ
Ｄ’（ｆ，ｐ）に対応するナビゲーション用データＳＷ
ｐ（ｐ’）を収集する。これを位相エンコード勾配ｐｅ
およびリワインド勾配ｐｒを変えながら、繰り返し時間
ＴＲでＭ回繰り返し、Ｍビュー分のイメージング用デー
タＤ’（ｆ，−Ｍ／２）〜Ｄ’（ｆ，Ｍ／２−１）およ
びそれらに対応するＭビュー分のナビゲーション用デー
タＳＲ-M/2（ｆ’）〜ＳＲM/2-1（ｆ’）およびナビゲ
ーション用データＳＷ-M/2（ｆ’）〜ＳＷ
M/2-1（ｆ’）を収集する。図１１は、上記ナビゲーシ
ョン用データＳＲｐ（ｆ’）およびＳＷｐ（ｐ’）に基
づいてイメージング用データＤ’（ｆ，ｐ）を位相補正
する処理のフローチャートである。ステップＳ１では、
ビュー番号カウンタｐ＝−Ｍ／２に初期化する。ステッ
プＳ２では、０ビューに対応するナビゲーション用デー
タＳＲｏ（ｆ’）と、ｐビューでのナビゲーション用デ
ータＳＲｐ（ｆ’）の共役ＳＲｐ*（ｆ'）の乗算を行
う。すなわち、 ＺＲｐ（ｆ’）＝ＳＲｐ*（ｆ’）・ＳＲｏ（ｆ’） を求める。ステップＳ３では、積ＺＲｐ（ｆ’）を１次
元フーリエ変換し、相互相関関数ｚｒｐ（ｘ）を求め
る。ステップＳ４では、最小二乗法を用いた補間近似に
より、相互相関関数ｚｒｐ（ｘ）の最大値を与える変数
値Δｘｐを求める。
向に移動しながら撮像する場合のパルスシーケンスであ
る。このパルスシーケンスＣでは、フリップ角α゜のＲ
ＦパルスＲαを印加すると共にｚ軸にスライス勾配ｓｓ
を印加する。次に、ｙ軸に位相エンコード勾配ｐｅを印
加する。次に、ｘ軸に周波数エンコード勾配ｒｉを印加
しながら、前記ＲＦパルスＲαから時間ＴＥ１後に発生
するイメージエコーＩＥをサンプリングして、ｐビュー
のイメージング用データＤ’（ｆ，ｐ）を収集する。次
に、リワインド勾配ｐｒを印加して位相エンコード量を
“０”に戻す。次に、前記周波数エンコード勾配ｒｉと
極性を反転したサンプリング数Ｎ’の読み出し勾配ｒｎ
１を印加しながら、前記ＲＦパルスＲαから時間ＴＥ２
後に発生するナビゲーションエコーＮＥ１をサンプリン
グして、ｐビューのイメージング用データＤ’（ｆ，
ｐ）に対応するナビゲーション用データＳＲｐ（ｆ’）
を収集する。次に、前記周波数エンコード勾配ｒｎ１と
極性を反転したリフェーズ勾配ｒｐを印加して周波数エ
ンコード量を“０”に戻す。また、位相軸（ｙ軸）にデ
ィフェーズ勾配ｄｐおよびサンプリング数Ｍ’の読み出
し勾配ｒｎ２を印加しながら、前記ＲＦパルスＲαから
時間ＴＥ３後に発生するナビゲーションエコーＮＥ２を
サンプリングして、ｐビューのイメージング用データ
Ｄ’（ｆ，ｐ）に対応するナビゲーション用データＳＷ
ｐ（ｐ’）を収集する。これを位相エンコード勾配ｐｅ
およびリワインド勾配ｐｒを変えながら、繰り返し時間
ＴＲでＭ回繰り返し、Ｍビュー分のイメージング用デー
タＤ’（ｆ，−Ｍ／２）〜Ｄ’（ｆ，Ｍ／２−１）およ
びそれらに対応するＭビュー分のナビゲーション用デー
タＳＲ-M/2（ｆ’）〜ＳＲM/2-1（ｆ’）およびナビゲ
ーション用データＳＷ-M/2（ｆ’）〜ＳＷ
M/2-1（ｆ’）を収集する。図１１は、上記ナビゲーシ
ョン用データＳＲｐ（ｆ’）およびＳＷｐ（ｐ’）に基
づいてイメージング用データＤ’（ｆ，ｐ）を位相補正
する処理のフローチャートである。ステップＳ１では、
ビュー番号カウンタｐ＝−Ｍ／２に初期化する。ステッ
プＳ２では、０ビューに対応するナビゲーション用デー
タＳＲｏ（ｆ’）と、ｐビューでのナビゲーション用デ
ータＳＲｐ（ｆ’）の共役ＳＲｐ*（ｆ'）の乗算を行
う。すなわち、 ＺＲｐ（ｆ’）＝ＳＲｐ*（ｆ’）・ＳＲｏ（ｆ’） を求める。ステップＳ３では、積ＺＲｐ（ｆ’）を１次
元フーリエ変換し、相互相関関数ｚｒｐ（ｘ）を求め
る。ステップＳ４では、最小二乗法を用いた補間近似に
より、相互相関関数ｚｒｐ（ｘ）の最大値を与える変数
値Δｘｐを求める。
【００４３】ステップＳ５では、０ビューに対応するナ
ビゲーション用データＳＷｏ（ｐ’）と、ｐビューでの
ナビゲーション用データＳＷｐ（ｐ’）の共役ＳＷｐ*
（ｐ'）の乗算を行う。すなわち、 ＺＷｐ（ｐ’）＝ＳＷｐ*（ｐ’）・ＳＷｏ（ｐ’） を求める。ステップＳ６では、積ＺＷｐ（ｐ’）を１次
元フーリエ変換し、相互相関関数ｚｗｐ（ｙ）を求め
る。ステップＳ７では、最小二乗法を用いた補間近似に
より、相互相関関数ｚｗｐ（ｙ）の最大値を与える変数
値Δｙｐを求める。ステップＳ８では、２πｆΔｘｐ／
Ｆｗ＝Δφｘ（ｐ）により位相補正量Δφｘ（ｐ）を算
出し、２πｐΔｙｐ／Ｐｗ＝Δφｙ（ｐ）により位相補
正量Δφｙ（ｐ）を算出し、イメージング用データＤ’
（ｆ，ｐ）に対して、 exp｛ｊ（Δφｘ（ｐ）＋Δφｙ（ｐ））｝ を乗算し、これをｐビューの位相補正後のイメージング
用データＤ（ｆ，ｐ）とする。ステップＳ９，Ｓ１０で
は、ｐ＝−Ｍ／２＋１からｐ＝Ｍ／２−１まで上記ステ
ップＳ２〜Ｓ８を反復実行し、ｋ空間を構成する全ビュ
ーのイメージング用データＤ（ｆ，ｐ）を取得する。
ビゲーション用データＳＷｏ（ｐ’）と、ｐビューでの
ナビゲーション用データＳＷｐ（ｐ’）の共役ＳＷｐ*
（ｐ'）の乗算を行う。すなわち、 ＺＷｐ（ｐ’）＝ＳＷｐ*（ｐ’）・ＳＷｏ（ｐ’） を求める。ステップＳ６では、積ＺＷｐ（ｐ’）を１次
元フーリエ変換し、相互相関関数ｚｗｐ（ｙ）を求め
る。ステップＳ７では、最小二乗法を用いた補間近似に
より、相互相関関数ｚｗｐ（ｙ）の最大値を与える変数
値Δｙｐを求める。ステップＳ８では、２πｆΔｘｐ／
Ｆｗ＝Δφｘ（ｐ）により位相補正量Δφｘ（ｐ）を算
出し、２πｐΔｙｐ／Ｐｗ＝Δφｙ（ｐ）により位相補
正量Δφｙ（ｐ）を算出し、イメージング用データＤ’
（ｆ，ｐ）に対して、 exp｛ｊ（Δφｘ（ｐ）＋Δφｙ（ｐ））｝ を乗算し、これをｐビューの位相補正後のイメージング
用データＤ（ｆ，ｐ）とする。ステップＳ９，Ｓ１０で
は、ｐ＝−Ｍ／２＋１からｐ＝Ｍ／２−１まで上記ステ
ップＳ２〜Ｓ８を反復実行し、ｋ空間を構成する全ビュ
ーのイメージング用データＤ（ｆ，ｐ）を取得する。
【００４４】上記位相補正したイメージング用データＤ
（ｆ，ｐ）から画像再構成し、イメージを生成すれば、
被検体をｘ軸方向およびｙ軸方向に移動しながら撮像し
ても、モーションアーチファクトを低減した高画質のイ
メージが得られる。
（ｆ，ｐ）から画像再構成し、イメージを生成すれば、
被検体をｘ軸方向およびｙ軸方向に移動しながら撮像し
ても、モーションアーチファクトを低減した高画質のイ
メージが得られる。
【００４５】−他の実施例− 第１実施例〜第３実施例では、グラディエントエコー法
のパルスシーケンスで説明したが、これに限定されず、
例えばスピンエコー法のパルスシーケンスにリワインド
勾配とナビゲーションエコーＮＥを読み出すための読み
出し勾配を加えたパルスシーケンスを用いてもよい。
のパルスシーケンスで説明したが、これに限定されず、
例えばスピンエコー法のパルスシーケンスにリワインド
勾配とナビゲーションエコーＮＥを読み出すための読み
出し勾配を加えたパルスシーケンスを用いてもよい。
【００４６】また、第１実施例〜第３実施例では、最小
二乗法により相互相関関数の最大値を与える変数値を求
めたが、Lagrange補間法，Chebyshev補間法，Hanning補
間法，Ｃubic補間法などを用いて相互相関関数の最大値
を与える変数値を求めてもよい。
二乗法により相互相関関数の最大値を与える変数値を求
めたが、Lagrange補間法，Chebyshev補間法，Hanning補
間法，Ｃubic補間法などを用いて相互相関関数の最大値
を与える変数値を求めてもよい。
【００４７】また、第１実施例〜第３実施例では、相関
定理を利用して位相補正量を算出したが、ナビゲーショ
ン用データの１次元フーリエ変換の絶対値から求まる投
影プロファイルの重心位置の変化や立ち上り／立ち下り
の位置の変化から位相補正量を算出してもよい。
定理を利用して位相補正量を算出したが、ナビゲーショ
ン用データの１次元フーリエ変換の絶対値から求まる投
影プロファイルの重心位置の変化や立ち上り／立ち下り
の位置の変化から位相補正量を算出してもよい。
【００４８】
【発明の効果】この発明のＭＲイメージング方法および
ＭＲＩ装置によれば、移動する被検体を撮像し且つアー
チファクトを抑制してイメージングすることが出来るよ
うになる。従って、被検体を連続的に移動させつつ撮像
し、高品質のイメージを得られるようになり、撮像効率
を向上することが出来る。
ＭＲＩ装置によれば、移動する被検体を撮像し且つアー
チファクトを抑制してイメージングすることが出来るよ
うになる。従って、被検体を連続的に移動させつつ撮像
し、高品質のイメージを得られるようになり、撮像効率
を向上することが出来る。
【図１】この発明のＭＲＩ装置の一実施例の構成図であ
る。
る。
【図２】被検体をｘ軸方向に移動しながら撮像するパル
スシーケンス図である。
スシーケンス図である。
【図３】ｋ−空間でのイメージング用データＤ’（ｆ，
ｐ）とナビゲーション用データＳＲｐ（ｆ’）の関係を
示す概念図である。
ｐ）とナビゲーション用データＳＲｐ（ｆ’）の関係を
示す概念図である。
【図４】被検体をｘ軸方向に移動しながら撮像する場合
の位相補正処理のフローチャートである。
の位相補正処理のフローチャートである。
【図５】被検体をｘ軸方向に移動しながら撮像する場合
の相互相関関数の概念図である。
の相互相関関数の概念図である。
【図６】被検体をｙ軸方向に移動しながら撮像するパル
スシーケンス図である。
スシーケンス図である。
【図７】ｋ−空間でのイメージング用データＤ’（ｆ，
ｐ）とナビゲーション用データＳＷｐ（ｐ’）の関係を
示す概念図である。
ｐ）とナビゲーション用データＳＷｐ（ｐ’）の関係を
示す概念図である。
【図８】被検体をｙ軸方向に移動しながら撮像する場合
の位相補正処理のフローチャートである。
の位相補正処理のフローチャートである。
【図９】被検体をｙ軸方向に移動しながら撮像する場合
の相互相関関数の概念図である。
の相互相関関数の概念図である。
【図１０】被検体をｘ軸方向およびｙ軸方向に移動しな
がら撮像するパルスシーケンス図である。
がら撮像するパルスシーケンス図である。
【図１１】被検体をｘ軸方向およびｙ軸方向に移動しな
がら撮像する場合の位相補正処理のフローチャートであ
る。
がら撮像する場合の位相補正処理のフローチャートであ
る。
【図１２】移動しない信号源とｋ−空間の概念図であ
る。
る。
【図１３】ｘ軸方向に移動した信号源とｋ−空間の概念
図である。
図である。
【図１４】従来のＭＲＩ装置の一例の構成図である。
１００ ＭＲＩ装置 １ マグネットアセンブリ ３ 勾配磁場駆動回路 ７ 計算機 ８ シーケンス記憶回路 １４ 移動装置 ＩＥ イメージエコー ＮＥ ナビゲーションエコー
Claims (5)
- 【請求項１】 ｋ−空間の周波数軸データ番号をｆ（−
Ｎ／２≦ｆ≦Ｎ／２−１）とし且つ位相軸データ番号を
ｐ（−Ｍ／２≦ｐ≦Ｍ／２−１）とするとき、ＲＦパル
スを印加し且つ位相エンコードを加え且つ周波数エンコ
ードを加え、発生するエコーからイメージング用データ
Ｄ’（ｆ，ｐ）を収集することを繰り返して、Ｍビュー
分のイメージング用データＤ’（ｆ，−Ｍ／２）〜Ｄ’
（ｆ，Ｍ／２−１）を収集するＭＲイメージング方法に
おいて、 １つのＲＦパルスに対して少なくとも２つのエコーを発
生させ、１つのエコーに対しては位相エンコード量を
“０”とし、周波数軸にサンプリング数Ｎ’の読み出し
勾配を印加しながら、データ番号ｆ’（−Ｎ'／２≦ｆ'
≦Ｎ’／２−１）のｐビューにおける周波数軸方向につ
いてのナビゲーション用データＳＲｐ（ｆ’）を収集
し、残りのエコーに対してはｋ−空間のｐビューに応じ
た位相エンコードを加え且つ周波数エンコードを加え、
これによりｋ−空間のｐビューのイメージング用データ
Ｄ’（ｆ，ｐ）を収集することを繰り返し、Ｍビュー分
のイメージング用データＤ’（ｆ，−Ｍ／２）〜Ｄ’
（ｆ，Ｍ／２−１）を収集し、 各ビューのナビゲーション用データＳＲ-M/2（ｆ’）〜
ＳＲM/2-1（ｆ'）の相互関係から位相補正量Δφｘ
（ｐ）を算出し、その位相補正量Δφｘ（ｐ）で各ビュ
ーのイメージング用データＤ’（ｆ，ｐ）を位相補正す
ることを特徴とするＭＲイメージング方法。 - 【請求項２】 ｋ−空間の周波数軸データ番号をｆ（−
Ｎ／２≦ｆ≦Ｎ／２−１）とし且つ位相軸データ番号を
ｐ（−Ｍ／２≦ｐ≦Ｍ／２−１）とするとき、ＲＦパル
スを印加し且つ位相エンコードを加え且つ周波数エンコ
ードを加え、発生するエコーからイメージング用データ
Ｄ’（ｆ，ｐ）を収集することを繰り返して、Ｍビュー
分のイメージング用データＤ’（ｆ，−Ｍ／２）〜Ｄ’
（ｆ，Ｍ／２−１）を収集するＭＲイメージング方法に
おいて、 １つのＲＦパルスに対して少なくとも２つのエコーを発
生させ、１つのエコーに対しては周波数エンコード量を
“０”とし、位相軸にサンプリング数Ｍ’の読み出し勾
配を印加しながら、データ番号ｐ’（−Ｍ’／２≦ｐ’
≦Ｍ’／２−１）のｐビューにおける位相軸方向につい
てのナビゲーション用データＳＷｐ（ｐ’）を収集し、
残りのエコーに対してはｋ−空間のｐビューに応じた位
相エンコードを加え且つ周波数エンコードを加え、これ
によりｋ−空間のｐビューのイメージング用データＤ’
（ｆ，ｐ）を収集することを繰り返し、Ｍビュー分のイ
メージング用データＤ’（ｆ，−Ｍ／２）〜Ｄ’（ｆ，
Ｍ／２−１）を収集し、 各ビューのナビゲーション用データＳＷ-M/2（ｐ’）〜
ＳＷM/2-1（ｐ'）の相互関係から位相補正量Δφｙ
（ｐ）を算出し、その位相補正量Δφｙ（ｐ）で各ビュ
ーのイメージング用データＤ’（ｆ，ｐ）を位相補正す
ることを特徴とするＭＲイメージング方法。 - 【請求項３】 ｋ−空間の周波数軸データ番号をｆ（−
Ｎ／２≦ｆ≦Ｎ／２−１）とし且つ位相軸データ番号を
ｐ（−Ｍ／２≦ｐ≦Ｍ／２−１）とするとき、ＲＦパル
スを印加し且つ位相エンコードを加え且つ周波数エンコ
ードを加え、発生するエコーからイメージング用データ
Ｄ’（ｆ，ｐ）を収集することを繰り返して、Ｍビュー
分のイメージング用データＤ’（ｆ，−Ｍ／２）〜Ｄ’
（ｆ，Ｍ／２−１）を収集するＭＲＩ装置において、 １つのＲＦパルスに対して少なくとも２つのエコーを発
生させ、１つのエコーに対しては位相エンコード量を
“０”とし、周波数軸にサンプリング数Ｎ’の読み出し
勾配を印加しながら、データ番号ｆ’（−Ｎ’／２≦
ｆ’≦Ｎ’／２−１）のｐビューにおける周波数軸方向
についてのナビゲーション用データＳＲｐ（ｆ’）を収
集し、残りのエコーに対してはｋ−空間のｐビューに応
じた位相エンコードを加え且つ周波数エンコードを加
え、これによりｋ−空間のｐビューのイメージング用デ
ータＤ’（ｆ，ｐ）を収集することを繰り返し、Ｍビュ
ー分のイメージング用データＤ’（ｆ，−Ｍ／２）〜
Ｄ’（ｆ，Ｍ／２−１）を収集するデータ収集手段と、 各ビューのナビゲーション用データＳＲ-M/2（ｆ’）〜
ＳＲM/2-1（ｆ'）の相互関係から位相補正量Δφｘ
（ｐ）を算出し、その位相補正量Δφｘ（ｐ）で各ビュ
ーのイメージング用データＤ’（ｆ，ｐ）を位相補正す
る位相補正手段とを具備したことを特徴とするＭＲＩ装
置。 - 【請求項４】 ｋ−空間の周波数軸データ番号をｆ（−
Ｎ／２≦ｆ≦Ｎ／２−１）とし且つ位相軸データ番号を
ｐ（−Ｍ／２≦ｐ≦Ｍ／２−１）とするとき、ＲＦパル
スを印加し且つ位相エンコードを加え且つ周波数エンコ
ードを加え、発生するエコーからイメージング用データ
Ｄ’（ｆ，ｐ）を収集することを繰り返して、Ｍビュー
分のイメージング用データＤ’（ｆ，−Ｍ／２）〜Ｄ’
（ｆ，Ｍ／２−１）を収集するＭＲＩ装置において、 １つのＲＦパルスに対して少なくとも２つのエコーを発
生させ、１つのエコーに対しては周波数エンコード量を
“０”とし、位相軸にサンプリング数Ｍ’の読み出し勾
配を印加しながら、データ番号ｐ’（−Ｍ’／２≦ｐ’
≦Ｍ’／２−１）のｐビューにおける位相軸方向につい
てのナビゲーション用データＳＷｐ（ｐ’）を収集し、
残りのエコーに対してはｋ−空間のｐビューに応じた位
相エンコードを加え且つ周波数エンコードを加え、これ
によりｋ−空間のｐビューのイメージング用データＤ’
（ｆ，ｐ）を収集することを繰り返し、Ｍビュー分のイ
メージング用データＤ’（ｆ，−Ｍ／２）〜Ｄ’（ｆ，
Ｍ／２−１）を収集するデータ収集手段と、 各ビューのナビゲーション用データＳＷ-M/2（ｐ’）〜
ＳＷM/2-1（ｐ'）の相互関係から位相補正量Δφｙ
（ｐ）を算出し、その位相補正量Δφｙ（ｐ）で各ビュ
ーのイメージング用データＤ’（ｆ，ｐ）を位相補正す
る位相補正手段とを具備したことを特徴とするＭＲＩ装
置。 - 【請求項５】 ｋ−空間の第１のエンコード軸のエンコ
ード量を“０”とし、ｋ−空間の第２のエンコード軸に
読み出し勾配を印加しながらナビゲーション用データＳ
１を収集し，その後，同様にしてナビゲーション用デー
タＳ２を収集するデータ収集手段と、前記ナビゲーショ
ン用データＳ１の１次元フーリエ変換の絶対値から求ま
る投影プロファイルと前記ナビゲーション用データＳ２
の１次元フーリエ変換の絶対値から求まる投影プロファ
イルの相互相関から前記第２のエンコード軸に対応する
実空間の方向についての信号源の移動量Δを算出する移
動量算出手段とを具備したことを特徴とするＭＲＩ装
置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP31787294A JP3384897B2 (ja) | 1994-12-21 | 1994-12-21 | Ｍｒｉ装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP31787294A JP3384897B2 (ja) | 1994-12-21 | 1994-12-21 | Ｍｒｉ装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08173396A true JPH08173396A (ja) | 1996-07-09 |
| JP3384897B2 JP3384897B2 (ja) | 2003-03-10 |
Family
ID=18093005
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP31787294A Expired - Lifetime JP3384897B2 (ja) | 1994-12-21 | 1994-12-21 | Ｍｒｉ装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3384897B2 (ja) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004000590A (ja) * | 2002-05-01 | 2004-01-08 | Ge Medical Systems Global Technology Co Llc | 高次の勾配磁場を用いて動いている被検体のｍｒデータを空間エンコーディングする方法 |
| US6707300B2 (en) | 2002-05-17 | 2004-03-16 | Ge Medical Systems Global Technology Co., Llc | Gradient non-linearity compensation in moving table MRI |
| US6794869B2 (en) | 2001-03-30 | 2004-09-21 | General Electric Company | Moving table MRI with frequency-encoding in the z-direction |
| US6801034B2 (en) | 2001-03-30 | 2004-10-05 | General Electric Company | Method and apparatus of acquiring large FOV images without slab-boundary artifacts |
| US6897655B2 (en) | 2001-03-30 | 2005-05-24 | General Electric Company | Moving table MRI with frequency-encoding in the z-direction |
| US6963768B2 (en) | 2002-05-16 | 2005-11-08 | General Electric Company | Whole body MRI scanning with moving table and interactive control |
| JP2006223869A (ja) * | 2005-02-18 | 2006-08-31 | Mayo Foundation For Medical Education & Research | 連続的テーブル移動により取得された磁気共鳴画像におけるアーチファクトの除去方法 |
| WO2007094174A1 (ja) * | 2006-02-13 | 2007-08-23 | Hitachi Medical Corporation | 磁気共鳴撮影装置および方法 |
-
1994
- 1994-12-21 JP JP31787294A patent/JP3384897B2/ja not_active Expired - Lifetime
Cited By (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6794869B2 (en) | 2001-03-30 | 2004-09-21 | General Electric Company | Moving table MRI with frequency-encoding in the z-direction |
| US6801034B2 (en) | 2001-03-30 | 2004-10-05 | General Electric Company | Method and apparatus of acquiring large FOV images without slab-boundary artifacts |
| US6891374B2 (en) | 2001-03-30 | 2005-05-10 | General Electric Company | Moving table MRI with frequency-encoding in the z-direction |
| US6897655B2 (en) | 2001-03-30 | 2005-05-24 | General Electric Company | Moving table MRI with frequency-encoding in the z-direction |
| JP2004000590A (ja) * | 2002-05-01 | 2004-01-08 | Ge Medical Systems Global Technology Co Llc | 高次の勾配磁場を用いて動いている被検体のｍｒデータを空間エンコーディングする方法 |
| US7558613B2 (en) | 2002-05-01 | 2009-07-07 | General Electric Company | Spatial encoding MR data of a moving subject using a higher-order gradient field |
| US6963768B2 (en) | 2002-05-16 | 2005-11-08 | General Electric Company | Whole body MRI scanning with moving table and interactive control |
| US7738944B2 (en) | 2002-05-16 | 2010-06-15 | General Electric Company | Whole body MRI scanning with moving table and interactive control |
| US6967479B2 (en) | 2002-05-17 | 2005-11-22 | Ge Medical Systems Global Technology Co., Llc | Gradient non-linearity compensation in moving table MRI |
| US6707300B2 (en) | 2002-05-17 | 2004-03-16 | Ge Medical Systems Global Technology Co., Llc | Gradient non-linearity compensation in moving table MRI |
| JP2006223869A (ja) * | 2005-02-18 | 2006-08-31 | Mayo Foundation For Medical Education & Research | 連続的テーブル移動により取得された磁気共鳴画像におけるアーチファクトの除去方法 |
| WO2007094174A1 (ja) * | 2006-02-13 | 2007-08-23 | Hitachi Medical Corporation | 磁気共鳴撮影装置および方法 |
| JP4991689B2 (ja) * | 2006-02-13 | 2012-08-01 | 株式会社日立メディコ | 磁気共鳴撮影装置および方法 |
| US8594766B2 (en) | 2006-02-13 | 2013-11-26 | Hitachi Medical Corporation | Magnetic resonance imaging apparatus and method |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3384897B2 (ja) | 2003-03-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3869337B2 (ja) | 磁気共鳴撮影装置 | |
| US7015696B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
| US4780675A (en) | Conjugate symmetry magnetic resonance imaging | |
| JP3117670B2 (ja) | マルチスライスｍｒイメージング方法およびｍｒｉ装置 | |
| JP4266577B2 (ja) | 磁気共鳴イメージング・システム | |
| US6222365B1 (en) | Magnetic resonance imaging apparatus and method | |
| CN103260510B (zh) | 磁共振成像装置以及对比度增强图像取得方法 | |
| JPH07323021A (ja) | Ｍｒイメージング装置 | |
| JP2755125B2 (ja) | Ｍｒイメージング装置 | |
| US6472873B2 (en) | Method and apparatus for capturing magnetic resonance signals, recording medium and system for magnetic resonance imaging using the same | |
| US20070285090A1 (en) | Phase cycling method and magnetic resonance imaging apparatus | |
| JP3384897B2 (ja) | Ｍｒｉ装置 | |
| JP3440049B2 (ja) | Ｍｒｉ装置 | |
| EP1191347A2 (en) | Tracking of the spin resonance frequency in magnetic resonance | |
| US4912413A (en) | Conjugate symmetry magnetic resonance imaging | |
| KR100559345B1 (ko) | 자기 공명 신호 획득 방법 및 장치, 기록 매체 및 자기공명 촬영 장치 | |
| JP3884227B2 (ja) | 磁気共鳴撮影装置 | |
| JP2008093022A (ja) | 磁気共鳴イメージング装置，磁気共鳴イメージング方法，スキャン装置，プログラム，および、記憶媒体 | |
| JP3538249B2 (ja) | Ｍｒｉｓ装置 | |
| JP4040745B2 (ja) | Ｍｒ装置 | |
| JPH06237910A (ja) | Ｍｒｉ装置 | |
| JP3322943B2 (ja) | Ｍｒｉ装置 | |
| JP7609622B2 (ja) | 磁気共鳴イメージング装置およびその制御方法 | |
| JP3720752B2 (ja) | ０次位相検出方法およびｍｒｉ装置 | |
| JPH11512957A (ja) | 磁気共鳴による画像化のための方法及び装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081227 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081227 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091227 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091227 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091227 Year of fee payment: 7 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091227 Year of fee payment: 7 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091227 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101227 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101227 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111227 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111227 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121227 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121227 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131227 Year of fee payment: 11 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| EXPY | Cancellation because of completion of term |