JPS6211438A - 核磁気共鳴映像装置 - Google Patents
核磁気共鳴映像装置Info
- Publication number
- JPS6211438A JPS6211438A JP60149736A JP14973685A JPS6211438A JP S6211438 A JPS6211438 A JP S6211438A JP 60149736 A JP60149736 A JP 60149736A JP 14973685 A JP14973685 A JP 14973685A JP S6211438 A JPS6211438 A JP S6211438A
- Authority
- JP
- Japan
- Prior art keywords
- magnetic resonance
- signal
- nuclear magnetic
- time window
- noise
- 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
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〔産業上の利用分野〕
この発明は、核磁気共鳴（ｎｕｃｌｅａｒ ｍａｇｎｅ
ｔｉｃｒｅｓｏｎａｎｃｅ 、以下ＮＲＭと略記する）
を医学診断のために応用したＮＭＲ映像装置に関するも
のである。
ｔｉｃｒｅｓｏｎａｎｃｅ 、以下ＮＲＭと略記する）
を医学診断のために応用したＮＭＲ映像装置に関するも
のである。
ＮＭＲは有機化合物の構造解析や物性物理の研究に有力
な手段として使われていた。
な手段として使われていた。
原子核はスピン（５ｐｉｎ ） と呼ばれるコマのよ
うな回転運動金しているので、磁石と同様の性質金持っ
ている。このような原子核を一様な静磁場ＨＯの中に置
くと、ラーモア周波数（ＬＢｒｍｏｒＦｒｅｑｕｅｎｃ
ｙ ）と呼ばれる原子核の種類に固有の角周波数ω で
歳差運動を行う。
うな回転運動金しているので、磁石と同様の性質金持っ
ている。このような原子核を一様な静磁場ＨＯの中に置
くと、ラーモア周波数（ＬＢｒｍｏｒＦｒｅｑｕｅｎｃ
ｙ ）と呼ばれる原子核の種類に固有の角周波数ω で
歳差運動を行う。
ω０＝γＨＯ・・・・・・・・・（１）でありｒは核の
種類に固有の定数である。こういう原子核の集団に電磁
波を加えれば、ω０に等しい角周波数の電磁波が吸収さ
れ、スピン状態間の遷移が起こる。これをＮＭＲという
。ＮＭＲｔ−使ってスピン系の向きを空間の特定方向に
そろえ九後で、印加した電磁波を切ればスピン系は元の
熱平衡の状態に戻る。その友めに要する時間を緩和時間
という。
種類に固有の定数である。こういう原子核の集団に電磁
波を加えれば、ω０に等しい角周波数の電磁波が吸収さ
れ、スピン状態間の遷移が起こる。これをＮＭＲという
。ＮＭＲｔ−使ってスピン系の向きを空間の特定方向に
そろえ九後で、印加した電磁波を切ればスピン系は元の
熱平衡の状態に戻る。その友めに要する時間を緩和時間
という。
ＮＭＲｔ″用いて被検体中に存在する成る特定の原子核
のスピン密度または緩和時間（緩和時定数の形で測定す
る）あるいはこれらに関連した情報を得ることは、たと
えば有機化合物の構造解析のための有力な手段であった
が、とのＮＭＲｔ医学診断に応用することができる。Ｎ
ＭＲｔ’医学診断に用いる場合はいわゆるＣＴ（ｃｏｍ
ｐｕｔｅｄ ｔｏｍｏｇｒａｐｈｙ 、コンピュータ断
面撮影法）の手法を用い、ＸｍＣＴや超音波ＣＴと同様
にブラウン管表示装置上に映像として表示する。ＣＴの
手法は従来よく知られてお夛、かつこの発明に対しては
直接の関係がないので説明を省略する。
のスピン密度または緩和時間（緩和時定数の形で測定す
る）あるいはこれらに関連した情報を得ることは、たと
えば有機化合物の構造解析のための有力な手段であった
が、とのＮＭＲｔ医学診断に応用することができる。Ｎ
ＭＲｔ’医学診断に用いる場合はいわゆるＣＴ（ｃｏｍ
ｐｕｔｅｄ ｔｏｍｏｇｒａｐｈｙ 、コンピュータ断
面撮影法）の手法を用い、ＸｍＣＴや超音波ＣＴと同様
にブラウン管表示装置上に映像として表示する。ＣＴの
手法は従来よく知られてお夛、かつこの発明に対しては
直接の関係がないので説明を省略する。
ＮＭＲ信号を検出する場合、問題となるのは、装置の電
源として用いられている交流電源からのハム等がＮＭＲ
信号に混入する雑音成分として検出され、これを除去す
ることが困難な点である。
源として用いられている交流電源からのハム等がＮＭＲ
信号に混入する雑音成分として検出され、これを除去す
ることが困難な点である。
検出されたＮＭＲ信号は複素フーリエ変換されて、ＣＴ
の方法により画像を再構成するが、交流電源からのハ
ムのような低周波のノイズが混入していると再構成画像
の中心部に偽像を生じる結果になる。
の方法により画像を再構成するが、交流電源からのハ
ムのような低周波のノイズが混入していると再構成画像
の中心部に偽像を生じる結果になる。
このような理由から、ＮＭＲの検出に際しては電源ハム
のような周期的雑音の混入を除去することが重要な問題
となる。特開昭５９−１６６８４６号「核磁気共鳴映像
装置」の明細書に開示される装置（以下先行技術という
）はこのような雑音を除去する手段を備えたＮＭＲ映像
装置である。第６図は上記先行技術の第２図と同一構成
を示すブロック図で、図において（１）は基準発振器（
式（１）のω０を発振する）、＋２）はＲＦ（ラジオ周
波数すなわちω０）パルス送信部、（３）は同調回路、
（４）は送受信コイル、（５）はＮＭＲ信号を増幅する
ための増幅器、（６）は１８σ分配器で角周波数ω０の
電圧を互に１８０°位相差のあるφ□、φ２−の位相に
分配する。（７）は装置全体の制御を行い、かつＣＴの
方法により画像の再構成を行う計算機、（８）は切換器
で、計算機（７）の制御によ！φ□とφ２のいずれかの
位相を切換えて同期検波の同期電圧として出力する。（
９）は角周波数ω０ に対する９０°移相器、（ＩＯＡ
）、（ＩＯＢ）はそれぞれ位相検波器で、増幅器１５）
の出力に対し同期検波を行う。
のような周期的雑音の混入を除去することが重要な問題
となる。特開昭５９−１６６８４６号「核磁気共鳴映像
装置」の明細書に開示される装置（以下先行技術という
）はこのような雑音を除去する手段を備えたＮＭＲ映像
装置である。第６図は上記先行技術の第２図と同一構成
を示すブロック図で、図において（１）は基準発振器（
式（１）のω０を発振する）、＋２）はＲＦ（ラジオ周
波数すなわちω０）パルス送信部、（３）は同調回路、
（４）は送受信コイル、（５）はＮＭＲ信号を増幅する
ための増幅器、（６）は１８σ分配器で角周波数ω０の
電圧を互に１８０°位相差のあるφ□、φ２−の位相に
分配する。（７）は装置全体の制御を行い、かつＣＴの
方法により画像の再構成を行う計算機、（８）は切換器
で、計算機（７）の制御によ！φ□とφ２のいずれかの
位相を切換えて同期検波の同期電圧として出力する。（
９）は角周波数ω０ に対する９０°移相器、（ＩＯＡ
）、（ＩＯＢ）はそれぞれ位相検波器で、増幅器１５）
の出力に対し同期検波を行う。
（ＩＩＡ）、（ＩＩＢ）はそれぞれ増幅器、（１２Ａ）
、（１２Ｂ）はそれぞれローパスフィルタ（以下ＬＰＦ
と略記する） 、（１３Ａ）＝（１３Ｂ） ｆ′ｉそれ
ぞれアナログディジタル変換器（以下Ａ／Ｄと略記する
）、（１４Ａ） 、 （１４Ｂ）はメモリである。また
、（１５）は計算機（７）で再構成された画像を表示す
るための画像表示装置でお次に動作について説明する。
、（１２Ｂ）はそれぞれローパスフィルタ（以下ＬＰＦ
と略記する） 、（１３Ａ）＝（１３Ｂ） ｆ′ｉそれ
ぞれアナログディジタル変換器（以下Ａ／Ｄと略記する
）、（１４Ａ） 、 （１４Ｂ）はメモリである。また
、（１５）は計算機（７）で再構成された画像を表示す
るための画像表示装置でお次に動作について説明する。
計算機（７）は除去しようとする周期性雑音（たとえば
交流ハム）の周期に同期してＲＦパルス送信部（２）を
制御し、基準発振器（１）からの角周波数ω０の電圧を
パルス変調して同調回路（３）を経てコイル（４）から
送出する。
交流ハム）の周期に同期してＲＦパルス送信部（２）を
制御し、基準発振器（１）からの角周波数ω０の電圧を
パルス変調して同調回路（３）を経てコイル（４）から
送出する。
ＮＭＲ現象が発生しコイル（４）、同調回路（３）ヲ経
て増幅器（５）に角周波数ω０の電圧が入力される。
て増幅器（５）に角周波数ω０の電圧が入力される。
増幅器（５）の出力のうち基準発振器（１）の出力と同
位相の成分は位相検波器（ＩＯＡ）で同期検波され、基
準発振器（１）の出力と９０°位相差のある成分は位相
検波器（ＩＯＢ）で同期検波された後、それぞれ増幅器
、Ｌ Ｐ Ｆ ｓ ｈ／Ｄ ’に経てディジタル信号と
してメモリ（１４Ａ）、（１４Ｂ）に一時記憶されて計
算機（７）における複素フーリエ分析に用いられるが、
このとき電源ハムのような低周波の雑音は位相検波器の
出力端からル０の出力端までの間等から混入する。
位相の成分は位相検波器（ＩＯＡ）で同期検波され、基
準発振器（１）の出力と９０°位相差のある成分は位相
検波器（ＩＯＢ）で同期検波された後、それぞれ増幅器
、Ｌ Ｐ Ｆ ｓ ｈ／Ｄ ’に経てディジタル信号と
してメモリ（１４Ａ）、（１４Ｂ）に一時記憶されて計
算機（７）における複素フーリエ分析に用いられるが、
このとき電源ハムのような低周波の雑音は位相検波器の
出力端からル０の出力端までの間等から混入する。
この雑音を除去するのが切換器（８）における切換であ
る。
る。
ＲＦパルス送信部（２）の変調パルスは除去すべき雑音
の位相に同期しているから毎回の測定において同一量の
雑音が混入する。しかし、位相検波器（１０Ａ） 、
（ｌ０Ｂ）に加えられる基準位相は毎回の測定において
φ からφ へ、φ からφ１へと１８０°反転するの
で、ｋツの出力では毎回その極性が反転する。
の位相に同期しているから毎回の測定において同一量の
雑音が混入する。しかし、位相検波器（１０Ａ） 、
（ｌ０Ｂ）に加えられる基準位相は毎回の測定において
φ からφ へ、φ からφ１へと１８０°反転するの
で、ｋツの出力では毎回その極性が反転する。
すなわち、切換器（８）が位相φ、を選択したとき、１
’ＪＤ （１３Ａ）、（１３Ｂ）の出力中信号成分をＸ
、（ｔ）、雑音成分をｎ□（１）とし、位相φ２ を選
択したとき、）ｖ’Ｄ （１３Ａ＞ 、 （１３Ｂ ）
の出力中信号成分をＸ２（ｔ）、雑音成分ｔｎ２（ｔ）
と仮定するとＸ□（ｔ） ＝ −Ｘ２（ｔ） 。
’ＪＤ （１３Ａ）、（１３Ｂ）の出力中信号成分をＸ
、（ｔ）、雑音成分をｎ□（１）とし、位相φ２ を選
択したとき、）ｖ’Ｄ （１３Ａ＞ 、 （１３Ｂ ）
の出力中信号成分をＸ２（ｔ）、雑音成分ｔｎ２（ｔ）
と仮定するとＸ□（ｔ） ＝ −Ｘ２（ｔ） 。
ｎｌ（ｔ） ｆｎ２（ｔ）となるのでＸ、 （ｔ）＋ｎ
、（ｔ）とｘ２（ｔ）＋ｎ２（ｔ）から雑音成分ｎ１（
ｔ） （＝ ｎ２（ｔ） ） ｆｔ’ｌｉｍすること
ができる。
、（ｔ）とｘ２（ｔ）＋ｎ２（ｔ）から雑音成分ｎ１（
ｔ） （＝ ｎ２（ｔ） ） ｆｔ’ｌｉｍすること
ができる。
先行技術によるＮＭＲ映像装置では、交流電源のハムに
よる雑音混入を除去するため以上のような処理を行って
いるので、１つのＮ Ｍ Ｒ信号を得るために２回収上
Ｕパルス送信部を動作させてそれぞれのＲＦパルス送信
に対して入力した信号の同期検波を行なわねばならず、
信号取得時間が長くなるという問題点があった。またＲ
Ｆパルスの送信を周期性雑音に同期させるためには、そ
の雑音成分の検出等を行う回路を付加する必要があった
。更に、たとえば、心臓のＮＭＲ−ＣＴ両画像構成する
ためにはＥＣＧ （ｅｌｅｃｔｒｏｃａｒｄｉｏｇｒａ
ｐｈ）信号に同期させてＲＦパルスを送信する必要があ
るため、このような場合には周期性雑音に同期してＲＦ
パルスを発生することが出来ないという問題点があった
。
よる雑音混入を除去するため以上のような処理を行って
いるので、１つのＮ Ｍ Ｒ信号を得るために２回収上
Ｕパルス送信部を動作させてそれぞれのＲＦパルス送信
に対して入力した信号の同期検波を行なわねばならず、
信号取得時間が長くなるという問題点があった。またＲ
Ｆパルスの送信を周期性雑音に同期させるためには、そ
の雑音成分の検出等を行う回路を付加する必要があった
。更に、たとえば、心臓のＮＭＲ−ＣＴ両画像構成する
ためにはＥＣＧ （ｅｌｅｃｔｒｏｃａｒｄｉｏｇｒａ
ｐｈ）信号に同期させてＲＦパルスを送信する必要があ
るため、このような場合には周期性雑音に同期してＲＦ
パルスを発生することが出来ないという問題点があった
。
この発明は上記のような問題点を解決するためになされ
たもので、ＬＰＦで除去できない低周波の周期性雑音を
特別な回路や装置を付加することなく除去できるＮＩＶ
ＩＲ映像装置を得ることを目的としている。
たもので、ＬＰＦで除去できない低周波の周期性雑音を
特別な回路や装置を付加することなく除去できるＮＩＶ
ＩＲ映像装置を得ることを目的としている。
この発明のＮＭＲ映像装置では、雑音の重畳した信号を
周波数分析したスペクトルから、雑音だけが存在する場
合の雑音１ｉ５波数分析したスペクトルを減算すること
によって信号だけの周波数スペクトルを算出し友。すな
わち、先行技術におけるように２回に分けて取得し九Ｎ
ＭＲ信号を記憶するかわりに１回だけのＮＭＲ信号とＮ
ＭＲ信号の存在しない場合の雑音波形とを記憶した。
周波数分析したスペクトルから、雑音だけが存在する場
合の雑音１ｉ５波数分析したスペクトルを減算すること
によって信号だけの周波数スペクトルを算出し友。すな
わち、先行技術におけるように２回に分けて取得し九Ｎ
ＭＲ信号を記憶するかわりに１回だけのＮＭＲ信号とＮ
ＭＲ信号の存在しない場合の雑音波形とを記憶した。
ＮＭＲ映像装置では必ずフーリエ変換によるスペクトル
分析が行われているので、雑音だけが存在する場合の雑
音のスペクトル分析も特別な装置を付加することなく容
易に行うことができ、周波数ドメイン（ｆｒｅｑｕｅｎ
ｃｙ ｄｏｍａｉｎ ） における処理により、容易
に雑音を除去することができる。
分析が行われているので、雑音だけが存在する場合の雑
音のスペクトル分析も特別な装置を付加することなく容
易に行うことができ、周波数ドメイン（ｆｒｅｑｕｅｎ
ｃｙ ｄｏｍａｉｎ ） における処理により、容易
に雑音を除去することができる。
以下この発明の実施例を図面について説明する。
第１図はこの発明の一実施例を示すブロック図であり、
第６図と同一符号は同−又は相当部分を示し、第６図の
構成に比し１８０°分配器（６）と切換器（８）が省略
されており、メモリ（１６Ａ）、（１６Ｂ）が追加され
ている。
第６図と同一符号は同−又は相当部分を示し、第６図の
構成に比し１８０°分配器（６）と切換器（８）が省略
されており、メモリ（１６Ａ）、（１６Ｂ）が追加され
ている。
ま念、第２図はメモリ（１４Ａ）、（１６Ａ）への書き
込みに関連する回路を示すブロック図であって、第１図
と同一符号は同一部分を示し、書き込み制御回路（１７
Ａ）は計算機（７）の制御に従ってメモリ（１４Ａ）、
（１６Ａ）への書込みを実行し、第１図では書。
込みに関連する回路を示すブロック図であって、第１図
と同一符号は同一部分を示し、書き込み制御回路（１７
Ａ）は計算機（７）の制御に従ってメモリ（１４Ａ）、
（１６Ａ）への書込みを実行し、第１図では書。
き込み制御回路（１７Ａ）が計算機（７）に内蔵される
として表わされている。メモリ（１４Ｂ）、（１６Ｂ）
への書き込み制御回路（図示せず）も同様に構成される
。
として表わされている。メモリ（１４Ｂ）、（１６Ｂ）
への書き込み制御回路（図示せず）も同様に構成される
。
第３図はＮＭＲ信号と周期性雑音との関係を示す波形図
で、図において横軸は時間ｔ、縦軸は撮幅を示し、ｘ
（ｔ）ｔＩｉＮ Ｍ Ｒ信号、ｎ （ｔ）は周期性雑音
である。周期性雑音の一周期はｔ。からｔ２、又はｔ
からｔ３の期間によって表すことができる。
で、図において横軸は時間ｔ、縦軸は撮幅を示し、ｘ
（ｔ）ｔＩｉＮ Ｍ Ｒ信号、ｎ （ｔ）は周期性雑音
である。周期性雑音の一周期はｔ。からｔ２、又はｔ
からｔ３の期間によって表すことができる。
！
第４図は第３図に示す信号及び雑音をＦＦＴ（Ｆａｓｔ
ｐｏｕｒｉｅｒ Ｔｒａｎｓｆｏｒｍ ） 処理に
より周波数スペクトルに変換した場合のスペクトル図を
示し図の横軸は周波数ｆ１縦軸はスペクトルの強さを示
し、Ｘ Ｔｆ）はｘ （ｔｌに対応し、Ｎ（ｆｎ）はｎ
（ｔ）に対応する。またＮ（ｆ ）は第３図には示して
ないが、白色雑音（ホワイトノイズ）のスペクトルで、
全周波数にわたり一様な強度分布を有することはよであ
り、電源ハムの場合、ｆｎ ！ｆｉ５０ Ｈｚ 又は
６０）Ｉｚおよびその高調波である。
ｐｏｕｒｉｅｒ Ｔｒａｎｓｆｏｒｍ ） 処理に
より周波数スペクトルに変換した場合のスペクトル図を
示し図の横軸は周波数ｆ１縦軸はスペクトルの強さを示
し、Ｘ Ｔｆ）はｘ （ｔｌに対応し、Ｎ（ｆｎ）はｎ
（ｔ）に対応する。またＮ（ｆ ）は第３図には示して
ないが、白色雑音（ホワイトノイズ）のスペクトルで、
全周波数にわたり一様な強度分布を有することはよであ
り、電源ハムの場合、ｆｎ ！ｆｉ５０ Ｈｚ 又は
６０）Ｉｚおよびその高調波である。
次に第１図に示す回路の動作について説明する。
ＩＶＤ（１３４）、（１３Ｂ）の入力点までの動作は第
１図に示す回路では第６図に示す回路に比し切換器（８
）における切換がない他は第６図に示す回路と同様であ
る。
１図に示す回路では第６図に示す回路に比し切換器（８
）における切換がない他は第６図に示す回路と同様であ
る。
一般にＦ Ｂ’ Ｔ処理を行う場合のデータ数ＮはＮ＝
２ （ｋは正の整数）とすることが望ましい。
２ （ｋは正の整数）とすることが望ましい。
サンプリング周期をＴｓ とし、第１番から第Ｎ番箇
でＮ個のサンプリング点を取ってＡ／Ｄ （１３Ａ）。
でＮ個のサンプリング点を取ってＡ／Ｄ （１３Ａ）。
（１３Ｂ）の入力をディジタル信号化し、メモＩＪ （
１４Ａ）。
１４Ａ）。
（１４Ｂ）の第１番地乃至第Ｎ番地にそれぞれ書き込む
と（第２図参照）、この書き込みにＴｆｉ（Ｎ−１）時
間を必要とする。第３図に示す例ではｔＯ時点から薔き
込みを開始しＴｓ（Ｎ−１）時間の後ｔ１時点において
書き込みを終ったとする。このことは第３図に示す信号
に対し時点ｔ。からｔ□に到る時間窓（ｗｉｎｄｏｗ
） によって切や出した信号をＦＦＴ処理したことに
なる。この処理結果が第４図に示すとおりであったとす
る。
と（第２図参照）、この書き込みにＴｆｉ（Ｎ−１）時
間を必要とする。第３図に示す例ではｔＯ時点から薔き
込みを開始しＴｓ（Ｎ−１）時間の後ｔ１時点において
書き込みを終ったとする。このことは第３図に示す信号
に対し時点ｔ。からｔ□に到る時間窓（ｗｉｎｄｏｗ
） によって切や出した信号をＦＦＴ処理したことに
なる。この処理結果が第４図に示すとおりであったとす
る。
第３図のｎ ｆｔｌにより第４図に出現するＮ ｆｆｌ
を除去するためにはｘ ｆｔｌが存在しない部分におい
て１（、−ｔ□ の時間窓（仮に第１の時間窓という）
に対するｎ ［ｔｌの位相関係と同様な位相関係にある
時間窓（仮に第２の時間窓という）によって信号ｎ （
ｔ１金切り出してＦＦＴ処理をすることが必要である。
を除去するためにはｘ ｆｔｌが存在しない部分におい
て１（、−ｔ□ の時間窓（仮に第１の時間窓という）
に対するｎ ［ｔｌの位相関係と同様な位相関係にある
時間窓（仮に第２の時間窓という）によって信号ｎ （
ｔ１金切り出してＦＦＴ処理をすることが必要である。
ｎ （ｔ）の周期と位相を測定することは容易であるの
で、ｔ１時点の後でｎ ｔｔｌがｔＱ時点と同位相にな
るｔ２ 時点のサンプリング点のデータを第（Ｎｉ１
）番のメモリアドレス（メモリ（１６Ａ）の第１番のア
ドレスと考えてもよい）に沓き込み、Ｔｓ （Ｎ−１）
時間後ｔ３時点に到って書き込みを停止しメモ！Ｊ （
１６Ａ） 。
で、ｔ１時点の後でｎ ｔｔｌがｔＱ時点と同位相にな
るｔ２ 時点のサンプリング点のデータを第（Ｎｉ１
）番のメモリアドレス（メモリ（１６Ａ）の第１番のア
ドレスと考えてもよい）に沓き込み、Ｔｓ （Ｎ−１）
時間後ｔ３時点に到って書き込みを停止しメモ！Ｊ （
１６Ａ） 。
（１６Ｂ）のＮ個のデータについてＦ’ＦＴ処理を施す
。
。
処理結果は第５図に示すとおりとなＦ） Ｎ’（ｆｎ）
はほぼＮ（ｆｎ）に等しく、またＮ’ （ｆ ）はほぼ
Ｎ （ｆｌに等しいので、第４図に示す結果と第５図に
示す結果とからｘ ｌｔｌに対応するＸ ｆｆｌだけを
抽出することは容易である。
はほぼＮ（ｆｎ）に等しく、またＮ’ （ｆ ）はほぼ
Ｎ （ｆｌに等しいので、第４図に示す結果と第５図に
示す結果とからｘ ｌｔｌに対応するＸ ｆｆｌだけを
抽出することは容易である。
なお、ＦＦＴの処理結果には時間窓の影響が入るので、
時間窓の始点と終点とはなるべ（ｎｆｔｌの０位相の点
に合致するよう、第３図の例で言えばｔ□ −ｔ２間の
時間窓、ｔ２− ｔ４間の時間窓のように設定できれば
よい。
時間窓の始点と終点とはなるべ（ｎｆｔｌの０位相の点
に合致するよう、第３図の例で言えばｔ□ −ｔ２間の
時間窓、ｔ２− ｔ４間の時間窓のように設定できれば
よい。
マ念、時間窓の幅をｎ （ｔｌの周期に比し十分に大き
くしてもよい。
くしてもよい。
さらに、上記実施例では第３図の時点ｔ□とｔ２との間
においてメモリへの書き込みを停止するとして説明した
が、メモリへの書き込みはｔｇ時点からｔ３時点まで連
続的に実行し、ＦＦＴ処理のデータはメモリに格納され
ているデータのうちからｔ。
においてメモリへの書き込みを停止するとして説明した
が、メモリへの書き込みはｔｇ時点からｔ３時点まで連
続的に実行し、ＦＦＴ処理のデータはメモリに格納され
ているデータのうちからｔ。
−１１間に書き込まれたものと、１２−１３間に書き込
まれたものとを用いるようにしてもよい。
まれたものとを用いるようにしてもよい。
なお、第４図と＠５図に示す例のように顕著な雑音成分
が周波数ｆｎの点だけに存在することが予期されるとき
は第４図のスペクトルと第５図のスペクトルの減算は周
波数ｆｎの点だけについて行なえばよいが、一般には低
周波の周期的雑音の周波数は複数であり、たとえば５０
Ｈｚの電源のハムは５０ ｆ（ｚの雑音の他にその高
調波の雑音をも含み、更には時間窓の影響によりｆｎ
の前後数点にサイドローブ的雑音が出現することもある
ので、一般には全周波数帯にわたり第４図のスペクトル
から第５図のスペクトルを減算する演算を行わねばなら
ぬ。
が周波数ｆｎの点だけに存在することが予期されるとき
は第４図のスペクトルと第５図のスペクトルの減算は周
波数ｆｎの点だけについて行なえばよいが、一般には低
周波の周期的雑音の周波数は複数であり、たとえば５０
Ｈｚの電源のハムは５０ ｆ（ｚの雑音の他にその高
調波の雑音をも含み、更には時間窓の影響によりｆｎ
の前後数点にサイドローブ的雑音が出現することもある
ので、一般には全周波数帯にわたり第４図のスペクトル
から第５図のスペクトルを減算する演算を行わねばなら
ぬ。
また、周期性雑音成分が複数ある場合は、これらの雑音
成分の周期の最小公倍数をＴＭとするときＴＭ／Ｔ８が
整数になるようなサンプリング周期を選び、周期性雑音
成分に相当する周波数においてだけスペクトルを減算す
る演算を行ってもよい。
成分の周期の最小公倍数をＴＭとするときＴＭ／Ｔ８が
整数になるようなサンプリング周期を選び、周期性雑音
成分に相当する周波数においてだけスペクトルを減算す
る演算を行ってもよい。
なお、この発明は広く一般的な雑音除去手段として利用
することができる。すなわち、信号の存在する場合の雑
音と信号の存在しない場合の雑音とがほぼ同一に保たれ
るとき、雑音の重畳した信号をスペクトル分析した結果
から、雑音だけをスペクトル分析した結果を減算して信
号だけのスペクトルを抽出し、この抽出した信号音フー
リエ逆変換すれば雑音の除去された信号波形を得ること
ができる。
することができる。すなわち、信号の存在する場合の雑
音と信号の存在しない場合の雑音とがほぼ同一に保たれ
るとき、雑音の重畳した信号をスペクトル分析した結果
から、雑音だけをスペクトル分析した結果を減算して信
号だけのスペクトルを抽出し、この抽出した信号音フー
リエ逆変換すれば雑音の除去された信号波形を得ること
ができる。
以上のようにこの発明によれば、周波数ドメインにおい
て雑音を除去するので、雑音除去が容易となり、かつＲ
Ｆパルスのタイミングを低周波の周期性雑音に同期する
必要がなくなったので、たとえば心臓の周期的動きに同
期して、かつ雑音により妨害されることのないＮＭＲ−
ＣＴ両画像得ることができる。
て雑音を除去するので、雑音除去が容易となり、かつＲ
Ｆパルスのタイミングを低周波の周期性雑音に同期する
必要がなくなったので、たとえば心臓の周期的動きに同
期して、かつ雑音により妨害されることのないＮＭＲ−
ＣＴ両画像得ることができる。
第１図はこの発明の一実施例を示すブロック図、第２図
は第１図のメモリへの書込みに関連する回路を示すブロ
ック図、第３図はＮＭＲ信号と周期性雑音との関係を示
す波形図、第４図及び第５図は第３図の波形の第１の時
間窓及びｗｃ２の時間窓内のスペクトル分析結果を示す
スペクトル図、第６図は従来の装置ｎ示すブロック図で
ある。 （１）は基準発振器、（７）は計算機、（９）は移相器
、（ＩＯＡ）。 （１０Ｂ）はそれぞれ位相検波器、（ＩＩＡ） 、 （
１１Ｂ）はそれぞれ増幅器、（１２Ａ）、（１２Ｂ）は
それぞれＬＰＦ％（１３Ａン、（１３Ｄ）はそれぞれ、
＾、／１〕、（１４Ａ） 、 （１４Ｂ） 、 （１６
Ａ） 。 （１６Ｂ）はそれぞれメモリ、 （１７Ａ）はメモリ
書き込み￥１ＩＩ−回路。 尚、各図中同一符号は同−又は相当部分を示す。
は第１図のメモリへの書込みに関連する回路を示すブロ
ック図、第３図はＮＭＲ信号と周期性雑音との関係を示
す波形図、第４図及び第５図は第３図の波形の第１の時
間窓及びｗｃ２の時間窓内のスペクトル分析結果を示す
スペクトル図、第６図は従来の装置ｎ示すブロック図で
ある。 （１）は基準発振器、（７）は計算機、（９）は移相器
、（ＩＯＡ）。 （１０Ｂ）はそれぞれ位相検波器、（ＩＩＡ） 、 （
１１Ｂ）はそれぞれ増幅器、（１２Ａ）、（１２Ｂ）は
それぞれＬＰＦ％（１３Ａン、（１３Ｄ）はそれぞれ、
＾、／１〕、（１４Ａ） 、 （１４Ｂ） 、 （１６
Ａ） 。 （１６Ｂ）はそれぞれメモリ、 （１７Ａ）はメモリ
書き込み￥１ＩＩ−回路。 尚、各図中同一符号は同−又は相当部分を示す。
Claims (2)
- （１）核磁気共鳴現象により誘起される信号を用いて被
検体内の特定原子核のスピン密度または緩和時定数ある
いはこれらに関連した情報を得て画像化する核磁気共鳴
映像装置において、 核磁気共鳴現象により誘起される上記信号を含む所定の
時間幅の第１の時間窓内の信号をフーリエ変換して第１
の周波数スペクトルを得る手段、核磁気共鳴現象により
誘起される上記信号は含まず、かつ上記信号の出力に重
畳する雑音だけが残存する時間領域において上記所定の
時間幅の第２の時間窓内の信号をフーリエ変換して第２
の周波数スペクトルを得る手段、 上記第１の周波数スペクトルから上記第２の周波数スペ
クトルを減算して核磁気共鳴現象により誘起される上記
信号に対する周波数スペクトルを抽出する手段を備えた
ことを特徴とする核磁気共鳴映像装置。 - （２）核磁気共鳴現象により誘起される信号の出力に重
畳する低周波の周期性雑音と上記第１の時間窓の間の位
相関係を、上記周期性雑音と上記第２の時間窓の間の位
相関係に合致させるよう上記第１の時間窓と上記第２の
時間窓の間隔を制御する手段を備えたことを特徴とする
特許請求の範囲第１項記載の核磁気共鳴映像装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP60149736A JPS6211438A (ja) | 1985-07-08 | 1985-07-08 | 核磁気共鳴映像装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP60149736A JPS6211438A (ja) | 1985-07-08 | 1985-07-08 | 核磁気共鳴映像装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS6211438A true JPS6211438A (ja) | 1987-01-20 |
| JPH0380013B2 JPH0380013B2 (ja) | 1991-12-20 |
Family
ID=15481684
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP60149736A Granted JPS6211438A (ja) | 1985-07-08 | 1985-07-08 | 核磁気共鳴映像装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS6211438A (ja) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62122646A (ja) * | 1985-11-22 | 1987-06-03 | 株式会社日立製作所 | 核磁気共鳴イメージング方法 |
| JPS63310149A (ja) * | 1987-06-12 | 1988-12-19 | Hitachi Cable Ltd | 高集積化ｉｃ用パッケ−ジおよびその製造方法 |
| JPS6462143A (en) * | 1987-08-31 | 1989-03-08 | Shimadzu Corp | Mri apparatus |
| JPH0663028A (ja) * | 1992-08-21 | 1994-03-08 | Hitachi Ltd | Ａ／ｄ変換データ処理装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61144552A (ja) * | 1984-12-18 | 1986-07-02 | Yokogawa Electric Corp | 核磁気共鳴デ−タの処理方法 |
| JPS61198042A (ja) * | 1985-02-28 | 1986-09-02 | Shimadzu Corp | Ｍｒ−ｃｔ装置の受信装置 |
-
1985
- 1985-07-08 JP JP60149736A patent/JPS6211438A/ja active Granted
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61144552A (ja) * | 1984-12-18 | 1986-07-02 | Yokogawa Electric Corp | 核磁気共鳴デ−タの処理方法 |
| JPS61198042A (ja) * | 1985-02-28 | 1986-09-02 | Shimadzu Corp | Ｍｒ−ｃｔ装置の受信装置 |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62122646A (ja) * | 1985-11-22 | 1987-06-03 | 株式会社日立製作所 | 核磁気共鳴イメージング方法 |
| JPS63310149A (ja) * | 1987-06-12 | 1988-12-19 | Hitachi Cable Ltd | 高集積化ｉｃ用パッケ−ジおよびその製造方法 |
| JPS6462143A (en) * | 1987-08-31 | 1989-03-08 | Shimadzu Corp | Mri apparatus |
| JPH0663028A (ja) * | 1992-08-21 | 1994-03-08 | Hitachi Ltd | Ａ／ｄ変換データ処理装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0380013B2 (ja) | 1991-12-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH0440015B2 (ja) | ||
| JPS62186852A (ja) | 核磁気分布決定方法及び装置 | |
| US3461381A (en) | Phase sensitive analog fourier analyzer readout for stored impulse resonance spectral data | |
| WO1987005201A1 (en) | Method of selective excitation in nmr imaging | |
| JPS649859B2 (ja) | ||
| JPS6211438A (ja) | 核磁気共鳴映像装置 | |
| US4766379A (en) | Method for obtaining nuclear magnetic resonance information data | |
| JPS6241648A (ja) | 非調和核磁気共鳴スピンエコー画像形成装置 | |
| US4689561A (en) | Nuclear magnetic resonance spectroscopy | |
| JPH10243934A (ja) | 磁気共鳴診断装置 | |
| JPS62139641A (ja) | Ｎｍｒイメ−ジング装置 | |
| JPS58116344A (ja) | 診断用核磁気共鳴装置 | |
| JP3279378B2 (ja) | Ｍｒ装置 | |
| JPH11216126A (ja) | 核磁気共鳴検査装置 | |
| GB2234595A (en) | Nuclear magnetic resonance imaging of solids | |
| JPS6222639A (ja) | 診断用磁気共鳴装置 | |
| Kozerke et al. | Spatial encoding and image reconstruction | |
| JP2970704B2 (ja) | 多次元核磁気共鳴測定方法 | |
| JPH09117423A (ja) | 核磁気共鳴検査装置 | |
| JPS62246355A (ja) | 磁気共鳴映像装置 | |
| JPS6044880A (ja) | チャ−プトラッキング方式を用いた核磁気共鳴イメ−ジング装置 | |
| JPS62295651A (ja) | 核磁気共鳴断層撮像装置 | |
| JPH0564635A (ja) | Ｍｒｉ装置における水，脂肪スペクトルの分離表示方法 | |
| JPS63212337A (ja) | 水・脂肪分離撮像法 | |
| JPS61269051A (ja) | 核磁気共鳴映像装置 |