JPS61128953A - Ｎｍｒ像の人為効果を減少する方法 - Google Patents
Ｎｍｒ像の人為効果を減少する方法Info
- Publication number
- JPS61128953A JPS61128953A JP60252996A JP25299685A JPS61128953A JP S61128953 A JPS61128953 A JP S61128953A JP 60252996 A JP60252996 A JP 60252996A JP 25299685 A JP25299685 A JP 25299685A JP S61128953 A JPS61128953 A JP S61128953A
- Authority
- JP
- Japan
- Prior art keywords
- pulse
- image
- phase
- nmr
- spin
- 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
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/567—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution gated by physiological signals, i.e. synchronization of acquired MR data with periodical motion of an object of interest, e.g. monitoring or triggering system for cardiac or respiratory gating
- G01R33/5676—Gating or triggering based on an MR signal, e.g. involving one or more navigator echoes for motion monitoring and correction
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
- A61B5/0816—Measuring devices for examining respiratory frequency
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7242—Details of waveform analysis using integration
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Signal Processing (AREA)
- General Health & Medical Sciences (AREA)
- Physiology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Pulmonology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Cardiology (AREA)
- Power Engineering (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
１１匹里点五里謁
この出願は出願人の係属中の米国特許出願通し番号筒６
７３，６８８号（１９８４年１１月２１日出願）と関連
を有する。
７３，６８８号（１９８４年１１月２１日出願）と関連
を有する。
■ユム亘旦
この発明は核磁気共鳴（ＮＭＲ＞方法及び装置に関する
。更に具体的に云えば、この発明はＮＭＲ走査器を用い
て呼吸速度を監視する方法と装置に関する。
。更に具体的に云えば、この発明はＮＭＲ走査器を用い
て呼吸速度を監視する方法と装置に関する。
ＮＭＲ方式を用いて胴体上部並びに腹部の診断上役立つ
像を取（ｑすることは、動きによって誘起される人為効
果（アーティファクト）の為に面倒になる場合が多い。
像を取（ｑすることは、動きによって誘起される人為効
果（アーティファクト）の為に面倒になる場合が多い。
動きの主な原因は、患者の呼吸である。こういう人為効
果の源は、普通「スピン捩れ形」と呼ばれる様なフーリ
エ変換作像方式を用いる時の、位相符号化次元に動きに
よって誘起された位相及び振幅誤差であることが判って
いる。こういう人為効果をなるべく小さく抑える為に数
多くの方式が提案されている。提案された方式の中には
、呼吸を止めること、呼気／吸気の間にデータ収集する
ことを含む種々のゲート方式、又は制御された浅い呼吸
がある。然し、こういう全ての方式は患者の協力を必要
とするか、或いは全呼吸サイクルを基準として走査時間
が長くなる。
果の源は、普通「スピン捩れ形」と呼ばれる様なフーリ
エ変換作像方式を用いる時の、位相符号化次元に動きに
よって誘起された位相及び振幅誤差であることが判って
いる。こういう人為効果をなるべく小さく抑える為に数
多くの方式が提案されている。提案された方式の中には
、呼吸を止めること、呼気／吸気の間にデータ収集する
ことを含む種々のゲート方式、又は制御された浅い呼吸
がある。然し、こういう全ての方式は患者の協力を必要
とするか、或いは全呼吸サイクルを基準として走査時間
が長くなる。
従来の提案された方式に伴う具体的な欠点は、例えば呼
吸ゲート方式は、ゲート信号を発生する為に余分の装置
を必要とすると共に、一般的に作像時間゛が長引く。呼
吸を止めること並びに浅い呼吸は、ボランチアに対して
はよく働くが、ＮＭＲ作象を受ける可能性が最も高い被
検体である重病の患者又は能力を失った患者にとっては
、困難であるか或いは不可能である。
吸ゲート方式は、ゲート信号を発生する為に余分の装置
を必要とすると共に、一般的に作像時間゛が長引く。呼
吸を止めること並びに浅い呼吸は、ボランチアに対して
はよく働くが、ＮＭＲ作象を受ける可能性が最も高い被
検体である重病の患者又は能力を失った患者にとっては
、困難であるか或いは不可能である。
、 従って、この発明の主な目的は、上に述べた従来
の方式の欠点をなくしながら、呼吸ザイクルとの非ゲー
ト形同期を保つことにより、動きによって誘起された人
為効果を減少し又はなくす方法を提供することである。
の方式の欠点をなくしながら、呼吸ザイクルとの非ゲー
ト形同期を保つことにより、動きによって誘起された人
為効果を減少し又はなくす方法を提供することである。
及」五里１
フーリエ変換作像パルス順序を使う時の、位相符号化次
元に動きによって誘起された誤差によるＮＭＲ像の人為
効果を減少する方法を提供する。
元に動きによって誘起された誤差によるＮＭＲ像の人為
効果を減少する方法を提供する。
この方法では、走査する被検体の呼吸サイクルを決定し
、この情報を使って、パルス順序の繰返し時間ＴＲを呼
吸サイクルの１７４の大体奇数倍に選ぶ。
、この情報を使って、パルス順序の繰返し時間ＴＲを呼
吸サイクルの１７４の大体奇数倍に選ぶ。
呼吸周期は任意の適当な方法によって決定することが出
来るが、好ましい実施例では、ＮＭＲ装置によって取得
した走査データを用いて、時間の関数として走査データ
の積分した大きさのスペクトルを監視することによって
決定される。
来るが、好ましい実施例では、ＮＭＲ装置によって取得
した走査データを用いて、時間の関数として走査データ
の積分した大きさのスペクトルを監視することによって
決定される。
この発明の新規と考えられる特徴は特許請求の範囲に具
体的に記載しであるが、この発明自体の構成、作用並び
にその他の目的及び利点は、以下図面について説明する
所から最もよく理解されよう。
体的に記載しであるが、この発明自体の構成、作用並び
にその他の目的及び利点は、以下図面について説明する
所から最もよく理解されよう。
及皿ユ用旦公Ｅ鴬
第１図はＮＭＲ作像装置のブロック図であり、これにつ
いてこの発明を説明する。装置を全体的に１００で示し
であるが、これはパルス制御モジュール１１２を含んで
おり、これがホスト・コンピュータ１１４の制御の下に
勾配磁界電源１１６に対して適正なタイミングのパルス
信号を発生する。この勾配磁界電源は、全体をブロック
１１８で示した勾配コイル集成体の一部分を形成する勾
配コイルを付勢する為に使われる。集成体にあるコイル
は、付勢された時、デカルト座標系の夫々Ｘ軸、ｙ軸及
びｚ軸方向の磁界勾配Ｇ２ ｅ Ｇｙ＊Ｇ２を発生する
。作像に勾配Ｇ−、Ｇｙ 、Ｇ２を使うことは、後で第
２図及び第３図について説明する。
いてこの発明を説明する。装置を全体的に１００で示し
であるが、これはパルス制御モジュール１１２を含んで
おり、これがホスト・コンピュータ１１４の制御の下に
勾配磁界電源１１６に対して適正なタイミングのパルス
信号を発生する。この勾配磁界電源は、全体をブロック
１１８で示した勾配コイル集成体の一部分を形成する勾
配コイルを付勢する為に使われる。集成体にあるコイル
は、付勢された時、デカルト座標系の夫々Ｘ軸、ｙ軸及
びｚ軸方向の磁界勾配Ｇ２ ｅ Ｇｙ＊Ｇ２を発生する
。作像に勾配Ｇ−、Ｇｙ 、Ｇ２を使うことは、後で第
２図及び第３図について説明する。
第１図について説明を続けると、パルス制御モジュール
が無線周波（ＲＦ）合成器１２０に対して作動パルスを
供給する。無線周波合成器１２０は、ＲＦ送受信装置の
一部分であり、この装置の一部分が破線のブロック１２
２の中に入っている。
が無線周波（ＲＦ）合成器１２０に対して作動パルスを
供給する。無線周波合成器１２０は、ＲＦ送受信装置の
一部分であり、この装置の一部分が破線のブロック１２
２の中に入っている。
パルス制御モジュールは送受信装置にある変調器１２４
にも変調信号を供給する。この変調器がＲＦ周波数合成
器の出力を変調する。変調されたＲＦ倍信号ＲＦ電力増
幅器１２８及び発信／受信スイッチ１３０を介してＲＦ
コイル集成体１２６に印加される。ＲＦ倍信号使って、
検査を受けるサンプル物体内の核スピンを励起する。こ
のサンプル物体はＲＦコイルの磁界の中に配置されてい
る。
にも変調信号を供給する。この変調器がＲＦ周波数合成
器の出力を変調する。変調されたＲＦ倍信号ＲＦ電力増
幅器１２８及び発信／受信スイッチ１３０を介してＲＦ
コイル集成体１２６に印加される。ＲＦ倍信号使って、
検査を受けるサンプル物体内の核スピンを励起する。こ
のサンプル物体はＲＦコイルの磁界の中に配置されてい
る。
励起された核スピンからのＮＭＲ信号を、核スピンの励
起に使ったのと同じ又は別の無線周波コイルで受信する
。受信信号を発信／受信スイッチを介してＲＦ前置増幅
器１５２に印加し、その後直角位相検波器１３４に印加
する。この検波器が同相及び直角位相出力Ｉ及びＱを夫
々発生する。
起に使ったのと同じ又は別の無線周波コイルで受信する
。受信信号を発信／受信スイッチを介してＲＦ前置増幅
器１５２に印加し、その後直角位相検波器１３４に印加
する。この検波器が同相及び直角位相出力Ｉ及びＱを夫
々発生する。
検出された信号をＡ／Ｄ変換器１３６でディジタル化し
、コンピュータ１１４に印加して、周知の形で処理し、
例えばサンプル物体のＮＭＲ像を再生する。
、コンピュータ１１４に印加して、周知の形で処理し、
例えばサンプル物体のＮＭＲ像を再生する。
この発明で用いるＮＭＲパルス順序は、ＮＭＲサンプル
物体１５０を示す第２図を最初に見れば、最もよく理解
されよう。この場合、サンプル物体は、患者テーブル１
５２の上に乗る部分的に示した仰臥の患者である。典型
的には、患者は座標系の２軸の正の方向を向く分極磁界
Ｂｏの中に配置されていて、それを縦方向に整合してい
る。座標系の原点は、破線１５６．１５８で定められた
厚さΔ２を持つ作像スライス１５４の中心にとるのが典
型的である。この場合９、スライ・スは、患者の胴体の
上側２／３の領域に選ばれているが、その厚さ並びに位
置が分極磁界Ｂ、ｏの強さ、使われる勾配磁界の大きざ
、並びにＲＦ励振パルスを変調する為に使われる周波数
成分によって決定されるが、これは当業者によく知られ
ている。
物体１５０を示す第２図を最初に見れば、最もよく理解
されよう。この場合、サンプル物体は、患者テーブル１
５２の上に乗る部分的に示した仰臥の患者である。典型
的には、患者は座標系の２軸の正の方向を向く分極磁界
Ｂｏの中に配置されていて、それを縦方向に整合してい
る。座標系の原点は、破線１５６．１５８で定められた
厚さΔ２を持つ作像スライス１５４の中心にとるのが典
型的である。この場合９、スライ・スは、患者の胴体の
上側２／３の領域に選ばれているが、その厚さ並びに位
置が分極磁界Ｂ、ｏの強さ、使われる勾配磁界の大きざ
、並びにＲＦ励振パルスを変調する為に使われる周波数
成分によって決定されるが、これは当業者によく知られ
ている。
第３図は、現在では２次元フーリエ変換（２ＤＦＴ）の
名前で知られている様な形式の普通の作像パルス順序と
呼ばれるもの１一部分である。このパルス順序は、「ス
ピン捩れ形」とも呼ばれる場合が多い。このパルス順序
は、例えば第２図に示した作像容積１５４の像を再生す
る為の作像データを周知の方法で求めるのに役立つ。こ
のパルス順序は、位相交番ＲＦ励振パルスを用いる点で
、従来のパルス順序を改良したものである。米国特許第
４．４４３，７６０号に記載されている様に、位相交番
ＲＦ励振パルスは位相交番のＮＭＲ信号を発生する。こ
ういう信号を減算すると、位相が反転されている信号成
分は強め合うが、その位相が反転していない望ましくな
いベースライン誤差成分は相殺される。
名前で知られている様な形式の普通の作像パルス順序と
呼ばれるもの１一部分である。このパルス順序は、「ス
ピン捩れ形」とも呼ばれる場合が多い。このパルス順序
は、例えば第２図に示した作像容積１５４の像を再生す
る為の作像データを周知の方法で求めるのに役立つ。こ
のパルス順序は、位相交番ＲＦ励振パルスを用いる点で
、従来のパルス順序を改良したものである。米国特許第
４．４４３，７６０号に記載されている様に、位相交番
ＲＦ励振パルスは位相交番のＮＭＲ信号を発生する。こ
ういう信号を減算すると、位相が反転されている信号成
分は強め合うが、その位相が反転していない望ましくな
いベースライン誤差成分は相殺される。
次に第３図についてスピン捩れ形作像パルス順序を更に
詳しく説明する。第３図には４つの位相符号化の図１．
１’ 、２．２’が示されている。
詳しく説明する。第３図には４つの位相符号化の図１．
１’ 、２．２’が示されている。
各対の図１−１′及び２−２′は位相符号化勾配Ｇｙ（
後で説明する）の同じ振幅を用いる。実際には、勾配Ｇ
ｙは例えば１２８個、２５６個又は５１２個の異なる振
幅を持つ様に選ぶことが出来る。第３図では、期間１（
横軸に沿って示す）に、正の磁界勾配パルスＧ２の存在
の下に印加される選択性９０°ＲＦ励撮パルスが示され
ている。。第１図のパルス制御モジュール１１２が周波
数合成器及び変調器に必要な制御信号を供給して、この
結果得られる励振パルスが、第２図に示した作像スライ
ス１４５の中にだけある核スピンを励起する為の正しい
周波数及び適切な変調がなされる様にする。典型的には
、励振パルスは５ｉｎｘ ／ｘ関数によって振幅変調し
て、励振された作像スライスが略矩形の輪郭を持つ様に
することが出来る。
後で説明する）の同じ振幅を用いる。実際には、勾配Ｇ
ｙは例えば１２８個、２５６個又は５１２個の異なる振
幅を持つ様に選ぶことが出来る。第３図では、期間１（
横軸に沿って示す）に、正の磁界勾配パルスＧ２の存在
の下に印加される選択性９０°ＲＦ励撮パルスが示され
ている。。第１図のパルス制御モジュール１１２が周波
数合成器及び変調器に必要な制御信号を供給して、この
結果得られる励振パルスが、第２図に示した作像スライ
ス１４５の中にだけある核スピンを励起する為の正しい
周波数及び適切な変調がなされる様にする。典型的には
、励振パルスは５ｉｎｘ ／ｘ関数によって振幅変調し
て、励振された作像スライスが略矩形の輪郭を持つ様に
することが出来る。
励振パルスの周波数は周知のラーマ方程式に従って、印
加された磁界の強さ並びに作像するＮＭＲ種目に関係す
る。パルス制御モジュールが勾配電源に作動信号をも印
加して、今の場合は勾配パルスＧ２を発生する。
加された磁界の強さ並びに作像するＮＭＲ種目に関係す
る。パルス制御モジュールが勾配電源に作動信号をも印
加して、今の場合は勾配パルスＧ２を発生する。
第３図について説明を続けると、勾配パルスＧ、、Ｇン
及びＧ、が期間２に同時に印加される。
及びＧ、が期間２に同時に印加される。
期間２の勾配Ｇ２は位相戻しパルスであって、典型的に
は、期間２にわたる勾配Ｇ２の波形の時間積分が、期間
１にわたる勾配Ｇ２の波形の時間積分の１／２に負の符
号を付したものに大体等しくなる様に選ぶのが典型的で
ある。このパルスＧ２の機能は、期間１に励起された核
スピンの位相戻しをすることである。前にも触れたが、
勾配パルスＧｙは位相符号化パルスであって、各々の図
１゜２・・・等で異なる振幅を持っていて、空間位相情
報をｙ軸方向に符号化する。
は、期間２にわたる勾配Ｇ２の波形の時間積分が、期間
１にわたる勾配Ｇ２の波形の時間積分の１／２に負の符
号を付したものに大体等しくなる様に選ぶのが典型的で
ある。このパルスＧ２の機能は、期間１に励起された核
スピンの位相戻しをすることである。前にも触れたが、
勾配パルスＧｙは位相符号化パルスであって、各々の図
１゜２・・・等で異なる振幅を持っていて、空間位相情
報をｙ軸方向に符号化する。
期間２の勾配パルスＧ、は、予定量だけ核スピンを位相
外しして、スピンエコー信号Ｓ１ （ｊ＞の発生時点を
期間４に遅らせるのに必要な位相外しパルスである。ス
ピンエコーが、典型的には非選択性１８０°ＲＦパルス
を期間３に印加することによって発生される。公知の様
に、１８０°ＲＦパルスは時間反転パルスであって、ス
ピンの位相外れの方向を反転する。これと線形の磁界勾
配パルスＧ８との組合せにより、スピンエコー信号が期
間４に発生される。空間情報をＸ軸方向に符号化する為
の線形の勾配パルスＧオの存在の下に、このスピンエコ
ー信号を標本化する。
外しして、スピンエコー信号Ｓ１ （ｊ＞の発生時点を
期間４に遅らせるのに必要な位相外しパルスである。ス
ピンエコーが、典型的には非選択性１８０°ＲＦパルス
を期間３に印加することによって発生される。公知の様
に、１８０°ＲＦパルスは時間反転パルスであって、ス
ピンの位相外れの方向を反転する。これと線形の磁界勾
配パルスＧ８との組合せにより、スピンエコー信号が期
間４に発生される。空間情報をＸ軸方向に符号化する為
の線形の勾配パルスＧオの存在の下に、このスピンエコ
ー信号を標本化する。
図１′は、期間５に印加される９０°ＲＦパルスが、負
の符号で表わす様に、期間１の９０”パルスに対して１
８０°位相がずれている他は、図１と略同−である。図
２−２’ 、３−３’・・・等のの間にも同じ関係があ
る。位相交番の結果として、期間８のスピンエコー信号
Ｓ２’（ｉ＞もスピンエコー信号Ｓ＋ （ｉ）に対し
て位相交番になる。然し、スピンエコー信号に先立つ、
不完全な１８０’ＲＦパルスによる擬似ＦＩＤ信号は位
相交番にならない。この為、信号３２（’ｊ）を信号Ｓ
＋ （ｔ）から減算すれば、スピンエコー信号は強め
合うが、望ましくない擬似ＦＩＤ信号が、直流オフセッ
トの様な他の誤差成分と共に相殺される。然し、９０’
ＲＦパルスの位相を反転して各々の図を繰返す結果、合
計の走査時間が２倍になる。
の符号で表わす様に、期間１の９０”パルスに対して１
８０°位相がずれている他は、図１と略同−である。図
２−２’ 、３−３’・・・等のの間にも同じ関係があ
る。位相交番の結果として、期間８のスピンエコー信号
Ｓ２’（ｉ＞もスピンエコー信号Ｓ＋ （ｉ）に対し
て位相交番になる。然し、スピンエコー信号に先立つ、
不完全な１８０’ＲＦパルスによる擬似ＦＩＤ信号は位
相交番にならない。この為、信号３２（’ｊ）を信号Ｓ
＋ （ｔ）から減算すれば、スピンエコー信号は強め
合うが、望ましくない擬似ＦＩＤ信号が、直流オフセッ
トの様な他の誤差成分と共に相殺される。然し、９０’
ＲＦパルスの位相を反転して各々の図を繰返す結果、合
計の走査時間が２倍になる。
図２の期間６で、勾配Ｇｙの振幅を次の値に増数し、期
間１６の図２′で繰返す。勾配Ｇｙを順次その振幅の範
囲全体にわたって進めるまで、この過程を繰返す。走査
過程が完了した時、収集されたデータに２次元フーリエ
変換過程を適用することにより、像情報が得られる。
間１６の図２′で繰返す。勾配Ｇｙを順次その振幅の範
囲全体にわたって進めるまで、この過程を繰返す。走査
過程が完了した時、収集されたデータに２次元フーリエ
変換過程を適用することにより、像情報が得られる。
第３図について注意すべき別の特徴は、９０゜ＲＦパル
スがＴＲの長さを持つ時間間隔で発生することである。
スがＴＲの長さを持つ時間間隔で発生することである。
期間ＴＲを繰返し時間と呼び、これはパルス順序（１つ
の図）の初めとその後に続く略同−のパルス順序の初め
との間の期間として定義される。繰返し時間を適当に調
節することにより、スピンエコー信号にＴ１依存性を導
入することが出来、これからＴ１加重像を再生すること
が出来る。
の図）の初めとその後に続く略同−のパルス順序の初め
との間の期間として定義される。繰返し時間を適当に調
節することにより、スピンエコー信号にＴ１依存性を導
入することが出来、これからＴ１加重像を再生すること
が出来る。
後の説明から明らかになるが、繰返し時間を慎重に選ぶ
ことにより、呼吸運動による人為効果を大幅に減少し又
はなくすことが出来る。こういう人為効果は、スピン捩
れ形作像順序を使う時、動きによって誘起された、位相
符号化方向（第３図に示したパルス順序の場合はｙ軸）
に於ける位相及び振幅の誤差が原因で起ることが判った
。この様な誤差による人為効果は、画素情報が機内の真
の位置か１ら単に変位したことによる像のぼやけ（分解
能の低下）からは識別すべきである。むしろこの発明の
方法によって減少し又は除去される人為効果は、位相符
号化方向に変位した、実際の（所望の）像の現実のゴー
スト像（人為効果）である。ゴースト像の変位は、物体
のフーリエ変換を表わす走査データが歪んでいることに
よるものである。この場合、位相誤差が実際にに一空間
に於ける画素位置のずれを表わし、それが変位の原因で
ある。
ことにより、呼吸運動による人為効果を大幅に減少し又
はなくすことが出来る。こういう人為効果は、スピン捩
れ形作像順序を使う時、動きによって誘起された、位相
符号化方向（第３図に示したパルス順序の場合はｙ軸）
に於ける位相及び振幅の誤差が原因で起ることが判った
。この様な誤差による人為効果は、画素情報が機内の真
の位置か１ら単に変位したことによる像のぼやけ（分解
能の低下）からは識別すべきである。むしろこの発明の
方法によって減少し又は除去される人為効果は、位相符
号化方向に変位した、実際の（所望の）像の現実のゴー
スト像（人為効果）である。ゴースト像の変位は、物体
のフーリエ変換を表わす走査データが歪んでいることに
よるものである。この場合、位相誤差が実際にに一空間
に於ける画素位置のずれを表わし、それが変位の原因で
ある。
この発明では、像に於ける被検体の実際の位置からの人
為効果の空間的な変位が、繰返し時間ＴＲと呼吸サイク
ルＢＣの比（ＴＲ／ＢＣ）に関係することが判った。こ
の関係が第４図に示されている。第４図で、像の中心線
１８０からの人為効果の変位を縦軸に沿って示してあり
、比ＴＲ／ＢＣを横軸に示しである。これが実際に線形
関係であって、位相符号化（又は仮想時間）領域に於、
けるナイキストの折畳みによる反復的な３角形の関係と
なって現われる。この関係の直線性は、像の中心線から
のゴースト変位を示す線分１８２−１８５から明らかで
ある。実線部分は、中心線１８００両側にある破線１８
６．１８８で区切られた装置の視野内のゴーストの変位
を表わす。理想化したＮＭＲ装置では、変位線１８２−
１８５の破線の延長で示す様に、変位線が視野の線を越
えて続く。然し、実際的な装置では、標本化速度の制約
の為、夫々実線の線分コ８６−１８９で示す様に、変位
線１８２−１８５が視野の中に折返される。この様なみ
かけの周期的な挙動の為、第３図について説明した様な
パルス順序でＴＲ＝ｎＢＣ／４（こ１でｎ＝１．３．５
−・・、即ちＴＲ＝０゜２５ＢＣ，０，７５ＢＣ、１．
２５８Ｃ等）にすると、その結果得られる像では、人為
効果の像が視野の縁に変位し、所望の像と干渉しない。
為効果の空間的な変位が、繰返し時間ＴＲと呼吸サイク
ルＢＣの比（ＴＲ／ＢＣ）に関係することが判った。こ
の関係が第４図に示されている。第４図で、像の中心線
１８０からの人為効果の変位を縦軸に沿って示してあり
、比ＴＲ／ＢＣを横軸に示しである。これが実際に線形
関係であって、位相符号化（又は仮想時間）領域に於、
けるナイキストの折畳みによる反復的な３角形の関係と
なって現われる。この関係の直線性は、像の中心線から
のゴースト変位を示す線分１８２−１８５から明らかで
ある。実線部分は、中心線１８００両側にある破線１８
６．１８８で区切られた装置の視野内のゴーストの変位
を表わす。理想化したＮＭＲ装置では、変位線１８２−
１８５の破線の延長で示す様に、変位線が視野の線を越
えて続く。然し、実際的な装置では、標本化速度の制約
の為、夫々実線の線分コ８６−１８９で示す様に、変位
線１８２−１８５が視野の中に折返される。この様なみ
かけの周期的な挙動の為、第３図について説明した様な
パルス順序でＴＲ＝ｎＢＣ／４（こ１でｎ＝１．３．５
−・・、即ちＴＲ＝０゜２５ＢＣ，０，７５ＢＣ、１．
２５８Ｃ等）にすると、その結果得られる像では、人為
効果の像が視野の縁に変位し、所望の像と干渉しない。
この方式は人為効果を減らす効果があることが実証され
、判っている。１例では、この発明に従って、任意の繰
返し時間ＴＲ＝５５０ミリ秒及びこの発明に従って選ん
だＴＲ＝８３３ミリ秒＝ＢＣ／８を使って、像を求めた
。ＴＲ＝５５０ミリ秒で得られた像は、所望の像と重な
る、動きによって誘起された人為効果があった。ＴＲ＝
８３３ミリ秒にした２番目の場合、人為効果を強める様
に（像の強度パラメータを制御する為に使われる）窓及
びレベルの設定をして、ｉ像した時でも、像の人為効果
は大幅に少なくなった。
、判っている。１例では、この発明に従って、任意の繰
返し時間ＴＲ＝５５０ミリ秒及びこの発明に従って選ん
だＴＲ＝８３３ミリ秒＝ＢＣ／８を使って、像を求めた
。ＴＲ＝５５０ミリ秒で得られた像は、所望の像と重な
る、動きによって誘起された人為効果があった。ＴＲ＝
８３３ミリ秒にした２番目の場合、人為効果を強める様
に（像の強度パラメータを制御する為に使われる）窓及
びレベルの設定をして、ｉ像した時でも、像の人為効果
は大幅に少なくなった。
例えばＴＲ＝０．５８Ｃ及び１ＢＣ（所望の像にゴース
ト像が重畳する場合に対応する。第４図の点１９０及び
１９２）に選ぶことにより、人為効果のない像が発生さ
れる様に思われるが、実際にはうまく行かない。その理
由は、ゴースト像と所望の像の間の整合がないこと、並
びに各々の像をぼかす人為効果の相加効果の為である。
ト像が重畳する場合に対応する。第４図の点１９０及び
１９２）に選ぶことにより、人為効果のない像が発生さ
れる様に思われるが、実際にはうまく行かない。その理
由は、ゴースト像と所望の像の間の整合がないこと、並
びに各々の像をぼかす人為効果の相加効果の為である。
毎分呼吸数が１０乃至２５の患者の呼吸速度では、この
発明で必要なＴＲ時間は０．６乃至１゜２秒の範囲内で
あり、ＴＲ＝１８Ｃ／４に相当する。この範囲内の繰返
し時間はデータ収集時間を著しく引伸さず、像のコント
ラストを大きく悪化させない位に短い。速く呼吸する患
者の様な成る場合、ＴＲ＝３８Ｇ／４を使うのが有利で
あることがある。これはＴＲ＝１゜５乃至２．０秒に対
応する。５８Ｃ／４及び７８Ｇ／４に対応する繰返し時
間は患者の検査には長過ぎるＴＲの値であるが、典型的
には自然の呼吸速度が高い実験室の動物に役立つことが
ある。部分的な飽和のスピンエコ一方式にこの発明の方
式を使うことによってコントラストが低下する場合、代
案は、飽和回復及び反転回復パルス順序の様な周知の他
の方式を用いて、Ｔ１加重像を求めることである。
発明で必要なＴＲ時間は０．６乃至１゜２秒の範囲内で
あり、ＴＲ＝１８Ｃ／４に相当する。この範囲内の繰返
し時間はデータ収集時間を著しく引伸さず、像のコント
ラストを大きく悪化させない位に短い。速く呼吸する患
者の様な成る場合、ＴＲ＝３８Ｇ／４を使うのが有利で
あることがある。これはＴＲ＝１゜５乃至２．０秒に対
応する。５８Ｃ／４及び７８Ｇ／４に対応する繰返し時
間は患者の検査には長過ぎるＴＲの値であるが、典型的
には自然の呼吸速度が高い実験室の動物に役立つことが
ある。部分的な飽和のスピンエコ一方式にこの発明の方
式を使うことによってコントラストが低下する場合、代
案は、飽和回復及び反転回復パルス順序の様な周知の他
の方式を用いて、Ｔ１加重像を求めることである。
呼吸の振幅を２倍に増加する効果（５ｃｉ及び１０ｃｖ
ｎのファントムの変位を用いる）は、ゴーストの強度を
強めると共に、高次のゴーストが出現し、実際の像から
第１のゴースト像までの距離に変化はないことが判った
。この為、ゴーストの変位の点で重要なパラメータはＴ
Ｒ／ＢＣ比（第４図）である。呼吸の振幅はゴーストの
強度にだけ影響する。
ｎのファントムの変位を用いる）は、ゴーストの強度を
強めると共に、高次のゴーストが出現し、実際の像から
第１のゴースト像までの距離に変化はないことが判った
。この為、ゴーストの変位の点で重要なパラメータはＴ
Ｒ／ＢＣ比（第４図）である。呼吸の振幅はゴーストの
強度にだけ影響する。
第４図の変位線（１８２−１８５＞は１次ゴースト像を
表わす。実際には第４図について１次ゴーストを説明し
たのと同様に、中心線から同じく変位するそれ程強度の
強くない高次のゴーストがあるが、第４図の破線１９４
で示す様に、周期が一層短くなる。高次のゴーストの人
為効果では、最大変位点が、前に述べたＴＲ−ｎＢＧ／
４の１次ゴーストと一致する。この様にして、この発明
の方法は高次の人為効果の影響を少なくし又は最小限に
抑えるのにも有効である。
表わす。実際には第４図について１次ゴーストを説明し
たのと同様に、中心線から同じく変位するそれ程強度の
強くない高次のゴーストがあるが、第４図の破線１９４
で示す様に、周期が一層短くなる。高次のゴーストの人
為効果では、最大変位点が、前に述べたＴＲ−ｎＢＧ／
４の１次ゴーストと一致する。この様にして、この発明
の方法は高次の人為効果の影響を少なくし又は最小限に
抑えるのにも有効である。
動きによって誘起された人為効果をなくし又は最小限に
抑える為の最適の繰返し時間を計算する。
抑える為の最適の繰返し時間を計算する。
為の呼吸速度の決定は、普通の外部呼吸モニタの様な任
意の適当な手段によって行なうことが出来る。然し、こ
の発明の好ましい実施例では、ＮＭＲ分光計自体をこの
目的の為に利用する。簡単に云うと、この為、呼吸の動
きがＮＭＲ信号の振幅変調を誘起する様な形で、ＮＭＲ
データを収集する。この後、変調された信号を解析する
ことによって呼吸速度を決定することが出来る。これを
どの様に実施するかを次に最初は第２図及び第３図につ
いて詳しく説明する。
意の適当な手段によって行なうことが出来る。然し、こ
の発明の好ましい実施例では、ＮＭＲ分光計自体をこの
目的の為に利用する。簡単に云うと、この為、呼吸の動
きがＮＭＲ信号の振幅変調を誘起する様な形で、ＮＭＲ
データを収集する。この後、変調された信号を解析する
ことによって呼吸速度を決定することが出来る。これを
どの様に実施するかを次に最初は第２図及び第３図につ
いて詳しく説明する。
時間の関数として積分したスピンエコーの大きさのスペ
クトルを監視することにより、呼吸速度を決定すること
が出来る。この目的の為、例えばＴＲ＝２００ミリ秒に
して、位相符号化勾配Ｇｙを使わずに、第３図のパルス
順序を用いて腹部（例えば第２図のスライス１５４）の
軸横断走査を行なうことにより、−組のデータを求める
。位相交番の図（１’ 、 ２’・・・等）は６４図の
合計には算入しないが、前に述べた様にベースライン誤
差を除く為に使われる。９０’ＲＦ励振パルスは、胴体
の上部、例えば第２図に参照数字１６０で示した横隔膜
の近くにある作像スライス１５４（第２図）を励振する
様に変調する。患者が呼吸すると、普通は作像スライス
の平面の外側にある横隔１１！（又はその他の組織の質
量）が、横隔膜を示す破線１６２及び両方向の矢印１６
３で表わす様に、スライスの容積に出入りする。スライ
ス１５４／））ら出るスピンエコー信号の強度が、呼吸
によって誘起された陽子スピン密度の変化並びに監視す
る作像スライス内の材料の移動により、呼吸サイクルと
共に変化する。
クトルを監視することにより、呼吸速度を決定すること
が出来る。この目的の為、例えばＴＲ＝２００ミリ秒に
して、位相符号化勾配Ｇｙを使わずに、第３図のパルス
順序を用いて腹部（例えば第２図のスライス１５４）の
軸横断走査を行なうことにより、−組のデータを求める
。位相交番の図（１’ 、 ２’・・・等）は６４図の
合計には算入しないが、前に述べた様にベースライン誤
差を除く為に使われる。９０’ＲＦ励振パルスは、胴体
の上部、例えば第２図に参照数字１６０で示した横隔膜
の近くにある作像スライス１５４（第２図）を励振する
様に変調する。患者が呼吸すると、普通は作像スライス
の平面の外側にある横隔１１！（又はその他の組織の質
量）が、横隔膜を示す破線１６２及び両方向の矢印１６
３で表わす様に、スライスの容積に出入りする。スライ
ス１５４／））ら出るスピンエコー信号の強度が、呼吸
によって誘起された陽子スピン密度の変化並びに監視す
る作像スライス内の材料の移動により、呼吸サイクルと
共に変化する。
第５図は、前に述べた６４図から成る順序を用い、スピ
ンエコー信号を標本化する標本化期間にわたって積分し
て得られた一組のデータの複素スピンエコー指数（大き
ざ〉のグラフである。長さ８．１９２ミリ秒のこういう
期間の１例は期間４（第３図）であり、この時１２８．
２５６又は５１２回に等しい回数だけ、スピンエコー信
号を標本化する。第５図で、スピンエコーの大きさを縦
軸に示し、横軸は図の番号（今の場合は１乃至６４）又
は時間であって、１目盛が２．５６秒である゛。パルス
順序の繰返し時間ＴＲは２００ミリ秒に選んだ。ベース
ライン誤差成分を除く為の位相交番パルス（図１’ 、
２’・・・等）を使う為、この結果、位相符号化勾配
Ｇｙの振幅の増分の間の周期は４００ミリ秒になる。
ンエコー信号を標本化する標本化期間にわたって積分し
て得られた一組のデータの複素スピンエコー指数（大き
ざ〉のグラフである。長さ８．１９２ミリ秒のこういう
期間の１例は期間４（第３図）であり、この時１２８．
２５６又は５１２回に等しい回数だけ、スピンエコー信
号を標本化する。第５図で、スピンエコーの大きさを縦
軸に示し、横軸は図の番号（今の場合は１乃至６４）又
は時間であって、１目盛が２．５６秒である゛。パルス
順序の繰返し時間ＴＲは２００ミリ秒に選んだ。ベース
ライン誤差成分を除く為の位相交番パルス（図１’ 、
２’・・・等）を使う為、この結果、位相符号化勾配
Ｇｙの振幅の増分の間の周期は４００ミリ秒になる。
第５図から明らかな様に、スピンエコーの大きさは循環
的な性質を持ち、これが呼吸サイクルに従うことが判っ
た。大きざのピーク間の平均時間を決定することにより
、１目盛２．５６秒であることに注意すれば、呼吸周期
を評価することが出来る。こうして呼吸サイクルは約３
．５秒であることが判る。
的な性質を持ち、これが呼吸サイクルに従うことが判っ
た。大きざのピーク間の平均時間を決定することにより
、１目盛２．５６秒であることに注意すれば、呼吸周期
を評価することが出来る。こうして呼吸サイクルは約３
．５秒であることが判る。
呼吸周期を測定する更に精密な方法は、（普通の様に周
波数を横軸にせずに）呼吸周期の関数として、スピンエ
コー指数のフーリエ変換を図（又は時間）に対してグラ
フに描くことである。こういうグラフが第６図に示され
ており、これは第５図のデータに対応する。スペクトル
のピーク２００が約３．７５秒の所で発生するが、これ
を呼吸周期とする。一般的に、呼吸周期を決定するには
、最低周波数で最大の大きざを持つスペクトル成分を突
止めることが望ましい。ピークをみつける計算機プログ
ラムを用いてスペクトル・データを解析し、最低周波数
で最大の振幅を持つピークをみつけ、こうして呼吸サイ
クルを決定することが有利である。
波数を横軸にせずに）呼吸周期の関数として、スピンエ
コー指数のフーリエ変換を図（又は時間）に対してグラ
フに描くことである。こういうグラフが第６図に示され
ており、これは第５図のデータに対応する。スペクトル
のピーク２００が約３．７５秒の所で発生するが、これ
を呼吸周期とする。一般的に、呼吸周期を決定するには
、最低周波数で最大の大きざを持つスペクトル成分を突
止めることが望ましい。ピークをみつける計算機プログ
ラムを用いてスペクトル・データを解析し、最低周波数
で最大の振幅を持つピークをみつけ、こうして呼吸サイ
クルを決定することが有利である。
２０１乃至２０７に示す様な一層小ざいピークは患者の
複雑な動きのスペクトルである。振幅が大きいこの様な
ピークが沢山あることは、患者のかなりの動きがあるこ
とを示すことが判った。経験によると、この様な状態で
は、有用な像を求めるのが困難であるのが普通である。
複雑な動きのスペクトルである。振幅が大きいこの様な
ピークが沢山あることは、患者のかなりの動きがあるこ
とを示すことが判った。経験によると、この様な状態で
は、有用な像を求めるのが困難であるのが普通である。
この様な動きが存在すること並びに強いことを監視し、
満足な像が得られるかどうかの判断に使う。この様な表
示を、実際に作像走査をぜずに、呼吸サイクルの予備検
査で有利に求めることが出来る。これによって走査器の
時間を更に効率よく使うことが出来る。
満足な像が得られるかどうかの判断に使う。この様な表
示を、実際に作像走査をぜずに、呼吸サイクルの予備検
査で有利に求めることが出来る。これによって走査器の
時間を更に効率よく使うことが出来る。
上に述べた呼吸周期を決定する方法の詳細は、第３図に
示す様なパルス順序を考えれば最も判り易い。この場合
、位相符号化勾配Ｇｙを印加せずに、例えば６４個のス
ピンエコーＳｉ （ｔ）、ｉ＝１，２・・・６４を観
察する。スピンエコーの大きさくスピンエコー指数Ｍ＋
＞は次の式に従って計算することが出来る。
示す様なパルス順序を考えれば最も判り易い。この場合
、位相符号化勾配Ｇｙを印加せずに、例えば６４個のス
ピンエコーＳｉ （ｔ）、ｉ＝１，２・・・６４を観
察する。スピンエコーの大きさくスピンエコー指数Ｍ＋
＞は次の式に従って計算することが出来る。
こ１でＴｓはスピンエコー信号を標本化する持続時間（
例えばＴｓ−８，１９２ミリ秒）であり、Ｒ（Ｓｉ）及
びＩ（３１は複素スピンエコーの実数及び虚数部分であ
る。
例えばＴｓ−８，１９２ミリ秒）であり、Ｒ（Ｓｉ）及
びＩ（３１は複素スピンエコーの実数及び虚数部分であ
る。
離散的な形では、Ｍｉ は次の様に表わすことが出来る
。
。
こ為でＥｌは離散的なスピンエコー信号のサンプルであ
り、Ｎｓは時間Ｔｓ内に求めたサンプルの数である。こ
の式は、走査データが実際にＮＭＲスピンエコー信号の
離散的なサンプルで構成されるので特に適している。
り、Ｎｓは時間Ｔｓ内に求めたサンプルの数である。こ
の式は、走査データが実際にＮＭＲスピンエコー信号の
離散的なサンプルで構成されるので特に適している。
スペクトル密度対図又は時間のグラフを求める為に必要
なエコー指数のフーリエ変換（スペクトル密度）は、次
の式に従って計算される。
なエコー指数のフーリエ変換（スペクトル密度）は、次
の式に従って計算される。
ＦＦＴ （Ｍ＋ ）ＦＫ （３
）こ）でに＝１．・−Ｎ図である。次に、ＦＫを縦軸に
とり、１／Ｋを横軸にとることにより、第６図に示す様
なグラフが得られる。
）こ）でに＝１．・−Ｎ図である。次に、ＦＫを縦軸に
とり、１／Ｋを横軸にとることにより、第６図に示す様
なグラフが得られる。
サンプル期間（８，１９２秒）にわたって積分した複素
スピンエコーの位相を使って呼吸周期を決定することが
出来るかどうかを判断しようとして、その為の解析も行
なった。図の番号に対する位相のグラフが第７図に示さ
れている。恐らく積分が全体のスピンエコ一時間に及ぶ
為、位相のグラフの解釈は大きざのグラフの解釈程容易
ではない。
スピンエコーの位相を使って呼吸周期を決定することが
出来るかどうかを判断しようとして、その為の解析も行
なった。図の番号に対する位相のグラフが第７図に示さ
れている。恐らく積分が全体のスピンエコ一時間に及ぶ
為、位相のグラフの解釈は大きざのグラフの解釈程容易
ではない。
フーリエ変換のグラフの解析も行なった。これは実際に
は、被検体を読出し軸に投影したことに相当する。こう
いうグラフは、積分が投影幅全体を含む時、第５図と非
常によく似ていた。投影の中心の１０ｃｍを構成するＦ
ＦＴサンプルだけを使った場合、変化の振幅は予想通り
に一層大きくなったが、目立つ様な余分の情報は得られ
なかった。
は、被検体を読出し軸に投影したことに相当する。こう
いうグラフは、積分が投影幅全体を含む時、第５図と非
常によく似ていた。投影の中心の１０ｃｍを構成するＦ
ＦＴサンプルだけを使った場合、変化の振幅は予想通り
に一層大きくなったが、目立つ様な余分の情報は得られ
なかった。
好ましい実施例では、スピンエコー指数（Ｍｉ ＞及
びＭｉ のフーリエ変換（ＦＫ）を計算機１１４で決定
した。これは前に第１図について説明したＮＭＲ装置の
一部分である。然しＭｌ及びＦＫを決定する為に専用の
回路装置を使うことが出来る。
びＭｉ のフーリエ変換（ＦＫ）を計算機１１４で決定
した。これは前に第１図について説明したＮＭＲ装置の
一部分である。然しＭｌ及びＦＫを決定する為に専用の
回路装置を使うことが出来る。
使うことの出来る装置の例を次に第８図について説明す
る。
る。
第８図には直角位相検波器１３４を示してあり、その入
力は前置増幅器１３２からのＮＭＲ信号を受取る様に結
合しである。検波器及び前置増幅器の両方は第１図の装
置と同じ装置１３４，１３２であってよい。直角出力■
及びＱが夫々実時間のアナログ積分回路２２０，２２２
に印加される。
力は前置増幅器１３２からのＮＭＲ信号を受取る様に結
合しである。検波器及び前置増幅器の両方は第１図の装
置と同じ装置１３４，１３２であってよい。直角出力■
及びＱが夫々実時間のアナログ積分回路２２０，２２２
に印加される。
積分器２２０．２２２の作用は、期間Ｔｓ （＝８゜
１９２ミリ秒）にわたって、パルス順序の各々の図の間
に標本化されたＮＭＲ信号を積分することである。この
目的の為、１つの図に対する積分動作が完了した時、（
例えば第１図のパルス制御モジュール１１２からのリセ
ット信号により）積分器をリセットし、次の図に対する
積分動作に備える。積分回路の出力を回路２２４，２２
６で自乗し、加算回路２２８で加算する。こうして得ら
れた和信号の平方根を平方根回路２３０で求める。
１９２ミリ秒）にわたって、パルス順序の各々の図の間
に標本化されたＮＭＲ信号を積分することである。この
目的の為、１つの図に対する積分動作が完了した時、（
例えば第１図のパルス制御モジュール１１２からのリセ
ット信号により）積分器をリセットし、次の図に対する
積分動作に備える。積分回路の出力を回路２２４，２２
６で自乗し、加算回路２２８で加算する。こうして得ら
れた和信号の平方根を平方根回路２３０で求める。
上に述べた動作により、式（１〉に従ってスピンエコー
指数Ｍｉの値が決定されることが判る。
指数Ｍｉの値が決定されることが判る。
平方根回路２３０に結合したモニタ２３２にスピンエコ
ー指数を表示することが出来る。この表示は、第５図に
示すグラフと同様であるが、次にこれを使って前に述べ
た様に呼吸周期を評価することが出来る。信号の平方根
を求める工程をなくしても、呼吸周期を評価する為のグ
ラフの使い方に大幅な変化はない。その影響は、第５図
に示すピークの振幅が変わることがあって、その時間に
対する位置は変わらない。
ー指数を表示することが出来る。この表示は、第５図に
示すグラフと同様であるが、次にこれを使って前に述べ
た様に呼吸周期を評価することが出来る。信号の平方根
を求める工程をなくしても、呼吸周期を評価する為のグ
ラフの使い方に大幅な変化はない。その影響は、第５図
に示すピークの振幅が変わることがあって、その時間に
対する位置は変わらない。
平方根回路２３０の出力をＡ／Ｄ変換器２３４１４に印
加する。表示装置２３６によって表示されるグラフは第
５図に示すものと同様であり、同じ目的の為に使うこと
が出来る。ディジタル化した信号を標本化し、ＦＦＴ装
置２３８に貯蔵する。
加する。表示装置２３６によって表示されるグラフは第
５図に示すものと同様であり、同じ目的の為に使うこと
が出来る。ディジタル化した信号を標本化し、ＦＦＴ装
置２３８に貯蔵する。
この装置２３８が式（３）に従ってスピンエコー指数の
フーリエ変換ＦＫを計算する。ＦＦＴ装置に結合した表
示装置２４０を使って、前に第６図について説明した様
に、縦軸にＦＫ、横軸に１／Ｋを表示する。計算機１１
４は、別個のＦＦＴ装置２３８の代りに、変換動作を行
なう様に利用することが出来る。
フーリエ変換ＦＫを計算する。ＦＦＴ装置に結合した表
示装置２４０を使って、前に第６図について説明した様
に、縦軸にＦＫ、横軸に１／Ｋを表示する。計算機１１
４は、別個のＦＦＴ装置２３８の代りに、変換動作を行
なう様に利用することが出来る。
この発明を特定の実施例及び例について説明したが、当
業者には、以上の説明からこの他の変更が考えられよう
。従って、この発明はこ１に具体的に説明した以外の形
でも実施することが出来ることを承知されたい。
業者には、以上の説明からこの他の変更が考えられよう
。従って、この発明はこ１に具体的に説明した以外の形
でも実施することが出来ることを承知されたい。
第１図はこの発明を実施するのに役立つＮＭＲ装置の一
部分の簡略ブロック図、第２図は選択励起過程によって
選択された作像スライスを含む仰臥患者の一部分を示す
図、第３図はこの発明を実施する時に利用することが出
来るスピン捩れ形パルス順序の１例を示すグラフ、第４
図は、平均の数を２に等しいとした場合の、呼吸速度Ｂ
Ｇに対するパルス順序繰返し時間ＴＲの比に対して、実
際の像からの人為効果のゴーストの変位を示すグラフ、
第５図は呼吸速度の情報を求めることが出来る様にする
、積分したスピンエコーの大きざを時間に対して示すグ
ラフ、第６図は第５図に示したスペクトル密度のフーリ
エ変換を求めることによって得られる、スピンエコー・
スペクトル密度対呼吸速度のグラフ、第７図は第６図と
同様なグラフであるが、積分したスピンエコー位相対時
間の関係を示す。第８図はこの発明に従って呼吸周期を
決定する為に用いることが出来る装置の１例を示すブロ
ック図である。 特許出願人ゼネラル・エレクトリックφカンノ−ニイ代
理人 （’１１６３０）生沼徳二 ＦＩＧ、１ ＦＩＧ、３ ｙｇＬ／ａζ ＦＩＧ、４
部分の簡略ブロック図、第２図は選択励起過程によって
選択された作像スライスを含む仰臥患者の一部分を示す
図、第３図はこの発明を実施する時に利用することが出
来るスピン捩れ形パルス順序の１例を示すグラフ、第４
図は、平均の数を２に等しいとした場合の、呼吸速度Ｂ
Ｇに対するパルス順序繰返し時間ＴＲの比に対して、実
際の像からの人為効果のゴーストの変位を示すグラフ、
第５図は呼吸速度の情報を求めることが出来る様にする
、積分したスピンエコーの大きざを時間に対して示すグ
ラフ、第６図は第５図に示したスペクトル密度のフーリ
エ変換を求めることによって得られる、スピンエコー・
スペクトル密度対呼吸速度のグラフ、第７図は第６図と
同様なグラフであるが、積分したスピンエコー位相対時
間の関係を示す。第８図はこの発明に従って呼吸周期を
決定する為に用いることが出来る装置の１例を示すブロ
ック図である。 特許出願人ゼネラル・エレクトリックφカンノ−ニイ代
理人 （’１１６３０）生沼徳二 ＦＩＧ、１ ＦＩＧ、３ ｙｇＬ／ａζ ＦＩＧ、４
Claims (1)
- 【特許請求の範囲】 １）フーリエ変換作像パルス順序を使う時の位相符号化
次元に於ける、動きによつて誘起された位相誤差による
ＮＭＲ像の人為効果を減少する方法に於て、 走査する被検体の呼吸サイクルを決定し、 パルス順序の繰返し時間ＴＲを呼吸サイクルＢＣの１／
４の大体奇数倍に選ぶ工程を含む方法。 ２）特許請求の範囲１）に記載した方法に於て、繰返し
時間ＴＲをＢＣ／４に選ぶ方法。 ３）特許請求の範囲１）に記載した方法に於て、繰返し
時間ＴＲを３ＢＣ／４に選ぶ方法。 ４）特許請求の範囲１）に記載した方法に於て、被検体
が動物であり、繰返し時間ＴＲを５ＢＣ／４に等しいか
又はそれより大きく選ぶ方法。 ５）特許請求の範囲１）に記載した方法に於て、フーリ
エ変換作像パルス順序がスピン捩れ形作像パルス順序で
ある方法。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US06/673,690 US4567893A (en) | 1984-11-21 | 1984-11-21 | Method of eliminating breathing artifacts in NMR imaging |
| US673690 | 1984-11-21 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPS61128953A true JPS61128953A (ja) | 1986-06-17 |
Family
ID=24703717
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP60252996A Pending JPS61128953A (ja) | 1984-11-21 | 1985-11-13 | Ｎｍｒ像の人為効果を減少する方法 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US4567893A (ja) |
| EP (1) | EP0182266A3 (ja) |
| JP (1) | JPS61128953A (ja) |
| IL (1) | IL76795A0 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4857418B2 (ja) * | 2005-02-16 | 2012-01-18 | アポロ メディカル イメージング テクノロジー ピーティーワイ リミテッド | 検体のモーションアーチファクト補正の方法およびシステム |
Families Citing this family (45)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB8417290D0 (en) * | 1984-07-06 | 1984-08-08 | Picker Int Ltd | Nuclear magnetic resonance method |
| EP0172345B1 (de) * | 1984-08-01 | 1989-05-10 | Siemens Aktiengesellschaft | Gerät zur Erzeugung von Bildern eines Untersuchungsobjektes |
| DE3430625A1 (de) * | 1984-08-20 | 1986-02-27 | Siemens AG, 1000 Berlin und 8000 München | Einrichtung fuer die kernspin-tomographie |
| US4616183A (en) * | 1984-10-22 | 1986-10-07 | General Electric Company | Method for reducing baseline error components in NMR signals |
| US4612504A (en) * | 1984-11-21 | 1986-09-16 | General Electric Company | Method for removing the effects of baseline error components in NMR imaging applications |
| NL8403627A (nl) * | 1984-11-29 | 1986-06-16 | Philips Nv | Werkwijze en inrichting voor het bepalen van een kernmagnetisatieverdeling in een deel van een lichaam. |
| DE3576408D1 (de) * | 1984-12-21 | 1990-04-12 | Philips Patentverwaltung | Verfahren zur erzeugung eines bewegungssignals und kernspintomograph fuer ein solches verfahren. |
| US4667159A (en) * | 1985-06-10 | 1987-05-19 | General Electric Company | Method of, and apparatus for, minimizing magnetic resonance imaging artifacts due to power line interference |
| US4724386A (en) * | 1985-09-30 | 1988-02-09 | Picker International, Inc. | Centrally ordered phase encoding |
| US4841248A (en) * | 1985-08-14 | 1989-06-20 | Picker International, Inc. | Transverse field limited localized coil for magnetic resonance imaging |
| US4663591A (en) * | 1985-08-16 | 1987-05-05 | General Electric Company | Method for reducing image artifacts due to periodic signal variations in NMR imaging |
| US4706026A (en) * | 1985-08-16 | 1987-11-10 | General Electric Company | Method for reducing image artifacts due to periodic signal variations in NMR imaging |
| US4728890A (en) * | 1985-08-16 | 1988-03-01 | Picker International, Inc. | Motion artifact suppression technique of magnetic resonance imaging |
| JPS62106755A (ja) * | 1985-11-02 | 1987-05-18 | 株式会社東芝 | 磁気共鳴イメ−ジング装置 |
| US4697149A (en) * | 1985-11-04 | 1987-09-29 | Wisconsin Alumni Research Foundation | NMR flow imaging using a composite excitation field and magnetic field gradient sequence |
| NL8503525A (nl) * | 1985-12-20 | 1987-07-16 | Philips Nv | Mri-werkwijze en inrichting ter reductie van beweginsartefacten. |
| JPS62231647A (ja) * | 1986-03-31 | 1987-10-12 | 株式会社東芝 | 磁気共鳴映像装置 |
| GB2190502B (en) * | 1986-05-16 | 1989-12-13 | Gen Electric Plc | Nuclear magnetic resonance methods and apparatus |
| NL8601845A (nl) * | 1986-07-15 | 1988-02-01 | Philips Nv | Mri-werkwijze en inrichting voor het reduceren van artefacten door middel van fasecodering. |
| US4903704A (en) * | 1986-08-14 | 1990-02-27 | North American Philips Corporation | Method and apparatus for MR imaging |
| JPS63164943A (ja) * | 1986-09-03 | 1988-07-08 | 株式会社日立製作所 | Ｎｍｒイメ−ジング方式 |
| US4710717A (en) * | 1986-12-29 | 1987-12-01 | General Electric Company | Method for fast scan cine NMR imaging |
| JP2644744B2 (ja) * | 1987-02-04 | 1997-08-25 | 株式会社日立製作所 | Ｎｍｒスペクトロスコピックイメージング装置 |
| IL82184A (en) * | 1987-04-10 | 1990-07-26 | Elscint Ltd | Reducing respiratory motion artifacts in nuclear magnetic resonance images |
| US4751462A (en) * | 1987-05-26 | 1988-06-14 | General Electric Company | Method for acquiring NMR data which is subject to periodic variations |
| FR2616936B1 (fr) * | 1987-06-19 | 1989-10-13 | Thomson Cgr | Procede de prise en compte, dans une image, des mouvements d'un objet |
| US4761613A (en) * | 1987-08-12 | 1988-08-02 | Picker International, Inc. | Monitored echo gating for the reduction of motion artifacts |
| JPH01148250A (ja) * | 1987-12-04 | 1989-06-09 | Hitachi Ltd | 磁気共鳴検査装置 |
| US5035244A (en) * | 1988-02-23 | 1991-07-30 | Elscint Ltd. | Motion artifact minimization |
| JPH01218439A (ja) * | 1988-02-26 | 1989-08-31 | Hitachi Ltd | 磁気共鳴イメージング装置 |
| US4961426A (en) * | 1988-08-19 | 1990-10-09 | Siemens Medical Systems, Inc. | Method for retrospectively gating NMR data |
| US4994965A (en) * | 1988-11-23 | 1991-02-19 | General Electric Company | Method for reducing motion induced image artifacts in projection imaging |
| US4982161A (en) * | 1989-08-24 | 1991-01-01 | North American Philips Corporation | Multimode magnetic resonance fast imaging method |
| US4994743A (en) * | 1989-10-27 | 1991-02-19 | General Electric Company | Method for monitoring respiration with acquired NMR data |
| US5251128A (en) * | 1990-11-19 | 1993-10-05 | General Electric Company | Motion artifact reduction in projection imaging |
| US5122747A (en) * | 1990-11-21 | 1992-06-16 | Mayo Foundation For Medical Education And Research | Spiral three-dimensional fourier transform NMR scan |
| DE4232897C2 (de) * | 1992-09-30 | 2002-10-31 | Bayer Ag | Flammwidrige thermoplastische Polycarbonat-Formmassen und Verfahren zu ihrer Herstellung |
| US5339035A (en) * | 1993-05-07 | 1994-08-16 | General Electric Company | MR imaging with rectangular magnetization transfer pulse |
| 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 |
| US6184682B1 (en) | 1998-04-15 | 2001-02-06 | Mayo Foundation For Medical Education And Research | Correction of MR images for motion artifacts using navigator echoes and autocorrection |
| US6647134B1 (en) * | 2000-03-30 | 2003-11-11 | Mayo Foundation For Medical Education And Research | Autocorrection of MR projection images |
| US7945305B2 (en) * | 2005-04-14 | 2011-05-17 | The Board Of Trustees Of The University Of Illinois | Adaptive acquisition and reconstruction of dynamic MR images |
| WO2010058732A1 (ja) * | 2008-11-18 | 2010-05-27 | 株式会社 日立メディコ | 磁気共鳴イメージング装置及び磁気共鳴イメージング方法 |
| TWI423656B (zh) * | 2010-12-31 | 2014-01-11 | Altek Corp | 數位影像之斜紋雜訊的檢測方法 |
| EP2713868B1 (en) * | 2011-05-23 | 2016-07-13 | Koninklijke Philips N.V. | Motion compensated magnetic resonance imaging method and system |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6125543A (ja) * | 1984-07-05 | 1986-02-04 | シーメンス、アクチエンゲゼルシヤフト | 臓器運動の無接触測定装置 |
| JPS6173642A (ja) * | 1984-09-19 | 1986-04-15 | 株式会社 日立メデイコ | デ−タ収集方法 |
| JPS6173643A (ja) * | 1984-09-19 | 1986-04-15 | 株式会社 日立メデイコ | デ−タ収集方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3871360A (en) * | 1973-07-30 | 1975-03-18 | Brattle Instr Corp | Timing biological imaging, measuring, and therapeutic timing systems |
| US3993995A (en) * | 1975-12-08 | 1976-11-23 | Rca Corporation | Respiration monitor |
| US4387722A (en) * | 1978-11-24 | 1983-06-14 | Kearns Kenneth L | Respiration monitor and x-ray triggering apparatus |
| US4443760A (en) * | 1982-07-01 | 1984-04-17 | General Electric Company | Use of phase alternated RF pulses to eliminate effects of spurious free induction decay caused by imperfect 180 degree RF pulses in NMR imaging |
| JPS59155239A (ja) * | 1983-02-23 | 1984-09-04 | 株式会社東芝 | 診断用核磁気共鳴装置 |
-
1984
- 1984-11-21 US US06/673,690 patent/US4567893A/en not_active Expired - Fee Related
-
1985
- 1985-10-23 IL IL76795A patent/IL76795A0/xx unknown
- 1985-11-12 EP EP85114392A patent/EP0182266A3/en not_active Withdrawn
- 1985-11-13 JP JP60252996A patent/JPS61128953A/ja active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6125543A (ja) * | 1984-07-05 | 1986-02-04 | シーメンス、アクチエンゲゼルシヤフト | 臓器運動の無接触測定装置 |
| JPS6173642A (ja) * | 1984-09-19 | 1986-04-15 | 株式会社 日立メデイコ | デ−タ収集方法 |
| JPS6173643A (ja) * | 1984-09-19 | 1986-04-15 | 株式会社 日立メデイコ | デ−タ収集方法 |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4857418B2 (ja) * | 2005-02-16 | 2012-01-18 | アポロ メディカル イメージング テクノロジー ピーティーワイ リミテッド | 検体のモーションアーチファクト補正の方法およびシステム |
Also Published As
| Publication number | Publication date |
|---|---|
| US4567893A (en) | 1986-02-04 |
| EP0182266A3 (en) | 1988-06-01 |
| IL76795A0 (en) | 1986-02-28 |
| EP0182266A2 (en) | 1986-05-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPS61128953A (ja) | Ｎｍｒ像の人為効果を減少する方法 | |
| US4564017A (en) | Method and apparatus for respiration monitoring with an NMR scanner | |
| US4614195A (en) | Method for reduction of motion artifacts in Fourier transform NMR imaging techniques | |
| US7425828B2 (en) | Frequency swept excitation for magnetic resonance | |
| JPS6279047A (ja) | 像の人為効果を減少する方法と装置 | |
| EP0121312B1 (en) | High speed nmr imaging system | |
| KR910003450B1 (ko) | 영상 아티팩트 저감 방법 | |
| EP1946137A1 (en) | Frequency swept excitation for magnetic resonance | |
| GB2206210A (en) | A method for acquiring nmr data which is subject to periodic variations | |
| US20130147481A1 (en) | Magnetic resonance imaging apparatus and vibrational error magnetic field reduction method | |
| US20150226821A1 (en) | System and method for gradient-modulated sweep imaging with fourier transformation magnetic resonance imaging | |
| EP0182267B1 (en) | A method for removing the effects of baseline error components in nmr imaging applications | |
| US4994743A (en) | Method for monitoring respiration with acquired NMR data | |
| JPS61100645A (ja) | サンプル物体のｎｍｒ検査を行なう方法 | |
| JP2018033691A (ja) | 磁気共鳴測定装置および画像処理方法 | |
| JPH10277010A (ja) | Ｍｒｉ装置 | |
| JPS63111847A (ja) | 信号を発生する方法と装置 | |
| EP0415683A2 (en) | NMR system | |
| JPH01274749A (ja) | 核磁気共鳴イメージング装置 | |
| JP2880716B2 (ja) | 磁気共鳴イメージング装置 | |
| EP0415682A2 (en) | NMR system | |
| JP2636649B2 (ja) | Ｍｒイメージング装置 | |
| JPH06114030A (ja) | Ｍｒイメージング装置 | |
| JPH10277005A (ja) | Ｍｒｉ装置 |