JPH0380013B2 - - Google Patents
Info
- Publication number
- JPH0380013B2 JPH0380013B2 JP60149736A JP14973685A JPH0380013B2 JP H0380013 B2 JPH0380013 B2 JP H0380013B2 JP 60149736 A JP60149736 A JP 60149736A JP 14973685 A JP14973685 A JP 14973685A JP H0380013 B2 JPH0380013 B2 JP H0380013B2
- Authority
- JP
- Japan
- Prior art keywords
- signal
- noise
- magnetic resonance
- time
- time window
- 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.)
- Expired - Lifetime
Links
- 238000001228 spectrum Methods 0.000 claims description 16
- 230000000737 periodic effect Effects 0.000 claims description 14
- 238000003384 imaging method Methods 0.000 claims description 7
- 238000002595 magnetic resonance imaging Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims 2
- 238000005481 NMR spectroscopy Methods 0.000 description 27
- 230000015654 memory Effects 0.000 description 16
- 238000010586 diagram Methods 0.000 description 9
- 230000001360 synchronised effect Effects 0.000 description 6
- 238000005070 sampling Methods 0.000 description 5
- 238000010183 spectrum analysis Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 102220642281 PTB domain-containing engulfment adapter protein 1_D13A_mutation Human genes 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 238000002247 constant time method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000000034 method Methods 0.000 description 2
- 150000002894 organic compounds Chemical class 0.000 description 2
- 238000012916 structural analysis Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 1
- 238000009987 spinning Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Description
【発明の詳細な説明】
〔産業上の利用分野〕
この発明は、該磁気共鳴（nuclear magnetic
resonance、以下NRMと略記する）を医学診断
のために応用したNMR映像装置に関するもので
ある。
resonance、以下NRMと略記する）を医学診断
のために応用したNMR映像装置に関するもので
ある。
NMRは有機化合物の構造解析や物性物理の研
究に有力な手段として使われていた。
究に有力な手段として使われていた。
原子核はスピン（spin）と呼ばれるコマのよう
な回転運動をしているので、磁石と同様の性質を
持つている。このような原子核を一様な静磁場
H0の中に置くと、ラーモア周波数（Larmor
Frequency）と呼ばれる原子核の種類に固有の角
周波数ωで歳差運動を行う。
な回転運動をしているので、磁石と同様の性質を
持つている。このような原子核を一様な静磁場
H0の中に置くと、ラーモア周波数（Larmor
Frequency）と呼ばれる原子核の種類に固有の角
周波数ωで歳差運動を行う。
ω0＝rH0 ……(1)
でありｒは核の種類に固有の定数である。こうい
う原子核の集団に電磁波を加えれば、ω0に等し
い角周波数の電磁波が吸収され、スピン状態間の
遷移が起こる。これをNMRという。NMRを使
つてスピン系の向きを空間の特定方向にそろえた
後で、印加した電磁波を切ればスピン系は元の熱
平衡の状態に戻る。そのために要する時間を緩和
時間という。
う原子核の集団に電磁波を加えれば、ω0に等し
い角周波数の電磁波が吸収され、スピン状態間の
遷移が起こる。これをNMRという。NMRを使
つてスピン系の向きを空間の特定方向にそろえた
後で、印加した電磁波を切ればスピン系は元の熱
平衡の状態に戻る。そのために要する時間を緩和
時間という。
NMRを用いて被検体中に存在する或る特定の
原子核のスピン密度または緩和時間（緩和時定数
の形で測定する）あるいはこれらに関連した情報
を得ることは、たとえば有機化合物の構造解析の
ための有力な手段であつたが、このNMRを医学
診断に応用することができる。NMRを医学診断
に用いる場合はいわゆるCT（computed
tomography、コンピユータ断面撮影法）手法を
用い、Ｘ線CTや超音波CTと同様にブラウン管表
示装置上に映像として表示する。CTの手法は従
来よく知られており、かつこの発明に対しては直
接の関係がないので説明を省略する。
原子核のスピン密度または緩和時間（緩和時定数
の形で測定する）あるいはこれらに関連した情報
を得ることは、たとえば有機化合物の構造解析の
ための有力な手段であつたが、このNMRを医学
診断に応用することができる。NMRを医学診断
に用いる場合はいわゆるCT（computed
tomography、コンピユータ断面撮影法）手法を
用い、Ｘ線CTや超音波CTと同様にブラウン管表
示装置上に映像として表示する。CTの手法は従
来よく知られており、かつこの発明に対しては直
接の関係がないので説明を省略する。
NMR信号を検出する場合、問題となるのは、
装置の電源として用いられている交流電源からの
ハム等がNMRに混入する雑音成分として検出さ
れ、これを除去することが困難な点である。検出
されたNMR信号は複素フーリエ変換されて、
CTの方法により画像を再構成するが、交流電源
からのハムのような低周波のノイズが混入してい
ると再構成画像の中心部に偽像を生じる結果にな
る。
装置の電源として用いられている交流電源からの
ハム等がNMRに混入する雑音成分として検出さ
れ、これを除去することが困難な点である。検出
されたNMR信号は複素フーリエ変換されて、
CTの方法により画像を再構成するが、交流電源
からのハムのような低周波のノイズが混入してい
ると再構成画像の中心部に偽像を生じる結果にな
る。
このような理由から、NMRの検出に際しては
電源ハムのよう周期的雑音の混入を除去すること
が重要の問題となる。特開昭59−166846号「該磁
気共鳴映像装置」の明細書に開示される装置（以
下先行技術という）はこのような雑音を除去する
手段を備えたNMR映像装置である。第６図は上
記先行技術の第２図と同一構成を示すブロツク図
で、図において１は基準発振器（式(1)のω0を発
振する）、２はRF（ラジオ周波数すなわちω0）パ
ルス信号部、３は同調回路、４は送受信コイル、
５はNMR信号を増幅するための増幅器、６は
180°分配器で角周波数ω0の電圧を互に180°位相差
のあるφ1，φ2の位相に分配する。７は装置全体
の制御を行い、かつCTの方法により画像の再構
成を行う計算機、８は切換器で、計算機７の制御
によりφ1とφ2のいずれかの位相を切変て同期検
波の同期電圧として出力する。９は角周波数ω0
に対する90°移相器、１０Ａ，１０Ｂはそれぞれ
位相検波器で、増幅器５の出力に対し同期検波を
行う。１１Ａ，１１Ｂはそれぞれ増幅器、１２
Ａ，１２Ｂはそれぞれローパスフイルタ（以下
LPFと略記する）、１３Ａ，１３Ｂはそれぞれア
ナログデイズタル変換器（以下Ａ／Ｄと略記す
る）、１４Ａ，１４Ｂはメモリである。また、１
５は計算機７で再構成された画像を表示するため
の画像表示装置である。
電源ハムのよう周期的雑音の混入を除去すること
が重要の問題となる。特開昭59−166846号「該磁
気共鳴映像装置」の明細書に開示される装置（以
下先行技術という）はこのような雑音を除去する
手段を備えたNMR映像装置である。第６図は上
記先行技術の第２図と同一構成を示すブロツク図
で、図において１は基準発振器（式(1)のω0を発
振する）、２はRF（ラジオ周波数すなわちω0）パ
ルス信号部、３は同調回路、４は送受信コイル、
５はNMR信号を増幅するための増幅器、６は
180°分配器で角周波数ω0の電圧を互に180°位相差
のあるφ1，φ2の位相に分配する。７は装置全体
の制御を行い、かつCTの方法により画像の再構
成を行う計算機、８は切換器で、計算機７の制御
によりφ1とφ2のいずれかの位相を切変て同期検
波の同期電圧として出力する。９は角周波数ω0
に対する90°移相器、１０Ａ，１０Ｂはそれぞれ
位相検波器で、増幅器５の出力に対し同期検波を
行う。１１Ａ，１１Ｂはそれぞれ増幅器、１２
Ａ，１２Ｂはそれぞれローパスフイルタ（以下
LPFと略記する）、１３Ａ，１３Ｂはそれぞれア
ナログデイズタル変換器（以下Ａ／Ｄと略記す
る）、１４Ａ，１４Ｂはメモリである。また、１
５は計算機７で再構成された画像を表示するため
の画像表示装置である。
次に動作について説明する。計算機７は除去し
ようとする周期性雑音（たとえば交流ハム）の周
期に同期してRFパルス信号部２を制御し、基準
発振器１からの角周波数ω0の電圧をパルス変調
して同調回路３を経てコイル４から送出する。
NMR現象が発生しコイル４、同調回路３を経て
増幅器５に角周波数ω0の電圧が入力される。
ようとする周期性雑音（たとえば交流ハム）の周
期に同期してRFパルス信号部２を制御し、基準
発振器１からの角周波数ω0の電圧をパルス変調
して同調回路３を経てコイル４から送出する。
NMR現象が発生しコイル４、同調回路３を経て
増幅器５に角周波数ω0の電圧が入力される。
増幅器５の出力のうち基準発振器１の出力と同
位相の成分は位相検波器１０Ａで同期検波され、
基準発振器１の出力と90°位相差のある成分は位
相検波器１０Ｂで同期検波された後、それぞれ増
幅器、LPF、Ａ／Ｄを経てデイジタル信号とし
てメモリ１４Ａ，１４Ｂに一時記憶されて計算機
７における複素フーリエ分析に用いられるが、こ
のとき電源ハムのような低周波の雑音は位相検波
器の出力端からＡ／Ｄの出力端までの間等から混
入する。この雑音を除去するのが切換器８におけ
る切換である。
位相の成分は位相検波器１０Ａで同期検波され、
基準発振器１の出力と90°位相差のある成分は位
相検波器１０Ｂで同期検波された後、それぞれ増
幅器、LPF、Ａ／Ｄを経てデイジタル信号とし
てメモリ１４Ａ，１４Ｂに一時記憶されて計算機
７における複素フーリエ分析に用いられるが、こ
のとき電源ハムのような低周波の雑音は位相検波
器の出力端からＡ／Ｄの出力端までの間等から混
入する。この雑音を除去するのが切換器８におけ
る切換である。
RFパルス送号部２の変調パルスは除去すべき
雑音の位相に同期しているから毎回の測定におい
て同一量の雑音が混入する。しかし、位相検波器
１０Ａ，１０Ｂに加えられる基準位相は毎回の測
定においてφ1からφ2へ、φ2からφへと180°反転す
るので、Ａ／Ｄの出力では毎回その極性が反転す
る。
雑音の位相に同期しているから毎回の測定におい
て同一量の雑音が混入する。しかし、位相検波器
１０Ａ，１０Ｂに加えられる基準位相は毎回の測
定においてφ1からφ2へ、φ2からφへと180°反転す
るので、Ａ／Ｄの出力では毎回その極性が反転す
る。
すなわち、切換器８が位相φ1を選択したとき、
Ａ／Ｄ１３Ａ，１３Ｂの出力中信号成分をX1
（ｔ）、雑音成分をn1（ｔ）とし、位相φ2を選択し
たとき、Ａ／Ｄ１３Ａ，１３Ｂの出力中信号成分
をX2（ｔ）、雑音成分をn2（ｔ）と仮定するとX1
（ｔ）＝−X2（ｔ）、n1（ｔ）＝n2X1となるのでX1
（ｔ）＋n1（ｔ）とX2（ｔ）＋n2（ｔ）から雑音成分
n1（ｔ）（＝n2（ｔ））を除去することができる。
Ａ／Ｄ１３Ａ，１３Ｂの出力中信号成分をX1
（ｔ）、雑音成分をn1（ｔ）とし、位相φ2を選択し
たとき、Ａ／Ｄ１３Ａ，１３Ｂの出力中信号成分
をX2（ｔ）、雑音成分をn2（ｔ）と仮定するとX1
（ｔ）＝−X2（ｔ）、n1（ｔ）＝n2X1となるのでX1
（ｔ）＋n1（ｔ）とX2（ｔ）＋n2（ｔ）から雑音成分
n1（ｔ）（＝n2（ｔ））を除去することができる。
先行技術によるNMR映像装置では、交流電源
のハムによる雑音混入を除去するため以上のよう
な処理を行つているので、１つのNMR信号を得
るために２回以上RFパルス送号部を動作させて
それぞれのRFパルス送信に対して入力した信号
の同期検波を行なわせねばならず、信号取得時間
が長くなるという問題点があつた。またRFパル
スの送信を同期性雑音に同期させるためには、そ
の雑音成分の検出等を行う回路を付加する必要が
あつた。更に、たとえば、心臓のNMR−CT画
像を構成するためにはECG（electrocardiograph）
信号に同期させてRFパルスを送信する必要があ
るため、このような場合には周期性雑音に同期し
てRFパルスを発生することが出来ないという問
題点があつた。
のハムによる雑音混入を除去するため以上のよう
な処理を行つているので、１つのNMR信号を得
るために２回以上RFパルス送号部を動作させて
それぞれのRFパルス送信に対して入力した信号
の同期検波を行なわせねばならず、信号取得時間
が長くなるという問題点があつた。またRFパル
スの送信を同期性雑音に同期させるためには、そ
の雑音成分の検出等を行う回路を付加する必要が
あつた。更に、たとえば、心臓のNMR−CT画
像を構成するためにはECG（electrocardiograph）
信号に同期させてRFパルスを送信する必要があ
るため、このような場合には周期性雑音に同期し
てRFパルスを発生することが出来ないという問
題点があつた。
この発明は上記のような問題点を解決するため
になされたもので、LPFで除去できない低周波
の周期性雑音を特別な回路や装置を付加すること
なく除去できるNMR映像装置を得ることを目的
としている。
になされたもので、LPFで除去できない低周波
の周期性雑音を特別な回路や装置を付加すること
なく除去できるNMR映像装置を得ることを目的
としている。
この発明のNMR映像装置では、雑音の重畳し
た信号を周波数分析したスペクトルから、雑音だ
けが存在する場合の雑音を周波数分析したスペク
トルを減算することによつて信号だけの周波数ス
ペクトルを算出した。すなわち、先行技術におけ
るように２回に分けて取得したNMR信号を記憶
するかわりに１回だけのNMR信号とNMR信号
の存在しない場合の雑音波形とを記憶した。
た信号を周波数分析したスペクトルから、雑音だ
けが存在する場合の雑音を周波数分析したスペク
トルを減算することによつて信号だけの周波数ス
ペクトルを算出した。すなわち、先行技術におけ
るように２回に分けて取得したNMR信号を記憶
するかわりに１回だけのNMR信号とNMR信号
の存在しない場合の雑音波形とを記憶した。
NMR映像装置では必ずフーリエ変換によりス
ペクトル分析が行われているので、雑音だけが存
在する場合の雑音のスペクトル分析も特別な装置
を付加することなく容易に行うことができ、周波
数ドメイン（frequency domain）における処理
により、容易に雑音を除去することができる。
ペクトル分析が行われているので、雑音だけが存
在する場合の雑音のスペクトル分析も特別な装置
を付加することなく容易に行うことができ、周波
数ドメイン（frequency domain）における処理
により、容易に雑音を除去することができる。
以下この発明の実施例を図面について説明す
る。第１図はこの発明の一実施例を示すブロツク
図であり、第６図と同一符号は同一又は相当部分
を示し、第６図の構成に比し180°分配器６と切換
器８が省略されており、メモリ１６Ａ，１６Ｂが
追加されている。
る。第１図はこの発明の一実施例を示すブロツク
図であり、第６図と同一符号は同一又は相当部分
を示し、第６図の構成に比し180°分配器６と切換
器８が省略されており、メモリ１６Ａ，１６Ｂが
追加されている。
また、第２図はメモリ１４Ａ，１６Ａへの書き
込みに関連する回路を示すブロツク図であつて、
第１図と同一符号は同一部分を示し、書き込み制
御回路１７Ａは計算機７の制御に従つてメモリ１
４Ａ，１６Ａへの書込みを実行し、第１図では書
き込み制御回路１７Ａが計算機７に内蔵されると
して表わされている。メモリ１４Ｂ，１６Ｂへの
書き込み制御回路（図示せず）も同様に構成され
る。
込みに関連する回路を示すブロツク図であつて、
第１図と同一符号は同一部分を示し、書き込み制
御回路１７Ａは計算機７の制御に従つてメモリ１
４Ａ，１６Ａへの書込みを実行し、第１図では書
き込み制御回路１７Ａが計算機７に内蔵されると
して表わされている。メモリ１４Ｂ，１６Ｂへの
書き込み制御回路（図示せず）も同様に構成され
る。
第３図はNMR信号と周期性雑音との関係を示
す波形図で、図において横軸は時間ｔ、縦軸は振
幅を示し、ｘ（ｔ）はNMR号、ｎ（ｔ）は周期性
雑音である。周期性雑音の一周期はt0からt2、又
はt1またt3の期間によつて表すことができる。
す波形図で、図において横軸は時間ｔ、縦軸は振
幅を示し、ｘ（ｔ）はNMR号、ｎ（ｔ）は周期性
雑音である。周期性雑音の一周期はt0からt2、又
はt1またt3の期間によつて表すことができる。
第４図は第３図に示す信号及び雑音をFFT
（Fast Fourier Tra sform）処理により周波数
スペクトルに変換した場合のスペクトル図を示し
図の横軸は周波数、縦軸はスペクトルの強さを
示し、Ｘ(f)はｘ（ｔ）に対応し、Ｎ（fn）はｎ（ｔ）
に対応する。またＮ(f)は第３図に示していない
が、白色雑音（ホワイトノイズ）のスペクトル
で、全周波数にわたり一様な強度分布を有するこ
とはよく知られている所である。またfn＝１／t2−t0 ＝１／t3−t1であり、電源ハムの場合、fnは50Hz又は 60Hzおよびその高調波である。
（Fast Fourier Tra sform）処理により周波数
スペクトルに変換した場合のスペクトル図を示し
図の横軸は周波数、縦軸はスペクトルの強さを
示し、Ｘ(f)はｘ（ｔ）に対応し、Ｎ（fn）はｎ（ｔ）
に対応する。またＮ(f)は第３図に示していない
が、白色雑音（ホワイトノイズ）のスペクトル
で、全周波数にわたり一様な強度分布を有するこ
とはよく知られている所である。またfn＝１／t2−t0 ＝１／t3−t1であり、電源ハムの場合、fnは50Hz又は 60Hzおよびその高調波である。
次の第１図に示す回路の動作について説明す
る。Ａ／Ｄ１３４，１３Ｂの入力点までの動作は
第１図に示す回路では第６図に示す回路に比し切
換器８における切換がない他は第６図に示す回路
と同様である。
る。Ａ／Ｄ１３４，１３Ｂの入力点までの動作は
第１図に示す回路では第６図に示す回路に比し切
換器８における切換がない他は第６図に示す回路
と同様である。
一般にFFT処理を行う場合のデータ数ＮはＮ
＝2k（ｋは正の整数）とすることが望ましい。サ
ンプリング周期をTsとし、第１番から第Ｎ番ま
でＮ個のサンプリング点を取つてＡ／Ｄ１３Ａ，
１３Ｂの入力をデイジタル信号化し、メモリ１４
Ａ，１４Ｂの第１番地乃至第Ｎ番地にそれぞれ書
き込むと（第２図参照）、この書き込みにTS（Ｎ
−１）時間を必要とする。第３図に示す例ではt0
時点から書き込みを開始しTS（Ｎ−１）時間の後
t1時点において書き込みを終つたとする。このこ
とは第３図に示す信号に対し時点t0からt1に到る
時間窓（window）によつて切り出した信号を
FFT処理したことになる。この処理結果が第４
図に示すとおりであつたとする。
＝2k（ｋは正の整数）とすることが望ましい。サ
ンプリング周期をTsとし、第１番から第Ｎ番ま
でＮ個のサンプリング点を取つてＡ／Ｄ１３Ａ，
１３Ｂの入力をデイジタル信号化し、メモリ１４
Ａ，１４Ｂの第１番地乃至第Ｎ番地にそれぞれ書
き込むと（第２図参照）、この書き込みにTS（Ｎ
−１）時間を必要とする。第３図に示す例ではt0
時点から書き込みを開始しTS（Ｎ−１）時間の後
t1時点において書き込みを終つたとする。このこ
とは第３図に示す信号に対し時点t0からt1に到る
時間窓（window）によつて切り出した信号を
FFT処理したことになる。この処理結果が第４
図に示すとおりであつたとする。
第３図のｎ（ｔ）により第４図に出現するＮ(f)
を除去するためにはｘ（ｔ）が存在しない部分に
おいてt0−t1の時間窓（仮に第１の時間窓とい
う）に対するｎ（ｔ）の位相関係と同様な位相関
係にある時間窓（仮に第２の時間窓という）によ
つて信号ｎ（ｔ）を切り出してFFT処理をするこ
とが必要である。ｎ（ｔ）と周期を測定すること
は容易であるので、t1時点の後でｎ（ｔ）がt0時
点と同位相になるt2時点のサンプリング点のデー
ダを第（Ｎ＋１）番のメモリアドレス（メモリ１
６Ａ）の第１番のアドレスと考えてもよい）に書
き込み、TS（Ｎ−１）時間後t3時点に到つて書き
込みを停止しメモリ１６Ａ，１６ＢのＮ個のデー
タについてFFT処理を施す。処理結果は第５図
に示すとおりとなりN′（fn）はほぼＮ（fn）に等
しく、またN′(f)はほぼＮ(f)に等しいので、第４
図に示す結果と第５図に示す結果とからｘ（ｔ）
に対応するｘ(f)だけを抽出することは容易であ
る。
を除去するためにはｘ（ｔ）が存在しない部分に
おいてt0−t1の時間窓（仮に第１の時間窓とい
う）に対するｎ（ｔ）の位相関係と同様な位相関
係にある時間窓（仮に第２の時間窓という）によ
つて信号ｎ（ｔ）を切り出してFFT処理をするこ
とが必要である。ｎ（ｔ）と周期を測定すること
は容易であるので、t1時点の後でｎ（ｔ）がt0時
点と同位相になるt2時点のサンプリング点のデー
ダを第（Ｎ＋１）番のメモリアドレス（メモリ１
６Ａ）の第１番のアドレスと考えてもよい）に書
き込み、TS（Ｎ−１）時間後t3時点に到つて書き
込みを停止しメモリ１６Ａ，１６ＢのＮ個のデー
タについてFFT処理を施す。処理結果は第５図
に示すとおりとなりN′（fn）はほぼＮ（fn）に等
しく、またN′(f)はほぼＮ(f)に等しいので、第４
図に示す結果と第５図に示す結果とからｘ（ｔ）
に対応するｘ(f)だけを抽出することは容易であ
る。
なお、時間窓の幅をｎ（ｔ）の周期に比し十分
に大きくしてもよい。
に大きくしてもよい。
さらに、上記実施例では第３図の時点t1とt2と
の間においてメモリへの書き込みを停止するとし
て説明したが、メモリへの書き込みはt0時点から
t3時点まで連続的に実行し、FFT処理のデータは
メモリに格納されているデータのうちからt0−t1
間に書き込まれたものと、t2−t3間に書き込まれ
たものと用いるようにしてもよい。
の間においてメモリへの書き込みを停止するとし
て説明したが、メモリへの書き込みはt0時点から
t3時点まで連続的に実行し、FFT処理のデータは
メモリに格納されているデータのうちからt0−t1
間に書き込まれたものと、t2−t3間に書き込まれ
たものと用いるようにしてもよい。
なお、第４図と第５図に示す例のように顕著な
雑音成分が周波数nの点だけに存在することが
予期されるときは第４図のスペクトルと第５図の
スペクトルの減算は周波数nの点だけについて
行なえばよい。
雑音成分が周波数nの点だけに存在することが
予期されるときは第４図のスペクトルと第５図の
スペクトルの減算は周波数nの点だけについて
行なえばよい。
また、周期性雑音成分が複数ある場合は、これ
らの雑音成分の周期の最小公倍数をTMとすると
きTM／TSが整数になるようなサンプリング周期
を選び、周期性雑音成分に相当する周波数におい
てだけスペクトルを減算する演算を行つてもよ
い。
らの雑音成分の周期の最小公倍数をTMとすると
きTM／TSが整数になるようなサンプリング周期
を選び、周期性雑音成分に相当する周波数におい
てだけスペクトルを減算する演算を行つてもよ
い。
なお、この発明は広く一般的な雑音除去手段と
して利用することができる。すなわち、信号の存
在する場合の雑音と信号の存在しない場合の雑音
とがほぼ同一に保たれるとき、雑音の重畳した信
号をスペクトル分析した結果から、雑音だけをス
ペクトル分析した結果を演算して信号だけのスペ
クトルを抽出し、この抽出した信号をフーリエ逆
変換すれば雑音の除去された信号波形を得ること
ができる。
して利用することができる。すなわち、信号の存
在する場合の雑音と信号の存在しない場合の雑音
とがほぼ同一に保たれるとき、雑音の重畳した信
号をスペクトル分析した結果から、雑音だけをス
ペクトル分析した結果を演算して信号だけのスペ
クトルを抽出し、この抽出した信号をフーリエ逆
変換すれば雑音の除去された信号波形を得ること
ができる。
以上のようにこの発明によれば、周波数ドメイ
ンにおいて雑音を除去するので、雑音除去が容易
となり、かつRFパルスのタイミングを低周波の
周期性雑音に同期する必要がなくなつたので、た
とえば心臓の周期的動きに同期して、かつ雑音に
より妨害されることのないNMR−CT画像を得
ることができる。
ンにおいて雑音を除去するので、雑音除去が容易
となり、かつRFパルスのタイミングを低周波の
周期性雑音に同期する必要がなくなつたので、た
とえば心臓の周期的動きに同期して、かつ雑音に
より妨害されることのないNMR−CT画像を得
ることができる。
第１図はこの発明の一実施例を示すブロツク
図、第２図は第１図のメモリへの書込みに関連す
る回路を示すブロツク図、第３図はNMR信号と
周期性雑音との関係を示す波形図、第４図及び第
５図は第３図の波形の第１の時間窓及び第２の時
間窓内のスペクトル分析結果を示すスペクトル
図、第６図は従来の装置を示すブロツク図であ
る。 １は基準発振器、７は決算機、９は移相器、１
０Ａ，１０Ｂはそれぞれ位相検波器、１１Ａ，１
１Ｂはそれぞれ増幅器、１２Ａ，１２Ｂはそれぞ
れLPF、１３Ａ，１３ＤはそれぞれＡ／Ｄ、１
４Ａ，１４Ｂ，１６Ａ，１６Ｂはそれぞれメモ
リ、１７Ａはメモリ書き込み制御回路。尚、各図
中同一符号は同一又は相当部分を示す。
図、第２図は第１図のメモリへの書込みに関連す
る回路を示すブロツク図、第３図はNMR信号と
周期性雑音との関係を示す波形図、第４図及び第
５図は第３図の波形の第１の時間窓及び第２の時
間窓内のスペクトル分析結果を示すスペクトル
図、第６図は従来の装置を示すブロツク図であ
る。 １は基準発振器、７は決算機、９は移相器、１
０Ａ，１０Ｂはそれぞれ位相検波器、１１Ａ，１
１Ｂはそれぞれ増幅器、１２Ａ，１２Ｂはそれぞ
れLPF、１３Ａ，１３ＤはそれぞれＡ／Ｄ、１
４Ａ，１４Ｂ，１６Ａ，１６Ｂはそれぞれメモ
リ、１７Ａはメモリ書き込み制御回路。尚、各図
中同一符号は同一又は相当部分を示す。
Claims (1)
- 【特許請求の範囲】 １ 該磁気共鳴現象により誘起される信号を用い
て被検体内の特定原子核のスピン密度または緩和
時定数あるいはこれらに関連した情報を得て画像
化する該磁気共鳴映像装置において、 該磁気共鳴現象により誘起される上記信号を含
む所定の時間幅の第１の時間窓内の信号をフーリ
エ変換して第１の周波数スペクトルを得る手段、 該磁気共鳴現象により誘起される上記信号は含
まず、かつ上記信号の出力に重畳する雑音だけが
残存する時間領域において上記所定の時間幅の第
２の時間窓内の信号をフーリエ変換して第２の周
波数スペクトルを得る手段、 上記第１の周波数スペクトルから上記第２の周
波数スペクトルを減算して該磁気共鳴現象により
誘起される上記信号に対する周波数スペクトルを
抽出する手段を備えたことを特徴とする該磁気共
鳴映像装置。 ２ 該磁気共鳴現象により誘起される信号の出力
に重畳する低周波の周期性雑音と上記第１の時間
窓の間の位相関係を、上記周期性雑音と上記第２
の時間窓の間の位相関係に合致させるよう上記第
１の時間窓と上記第２の時間窓の間隔を制御する
手段を備えたことを特徴とする特許請求の範囲第
１項記載の該磁気共鳴映像装置。
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 JPS6211438A (ja) | 1987-01-20 |
| JPH0380013B2 true 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) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0767442B2 (ja) * | 1985-11-22 | 1995-07-26 | 株式会社日立製作所 | 核磁気共鳴イメージング方法 |
| JPS63310149A (ja) * | 1987-06-12 | 1988-12-19 | Hitachi Cable Ltd | 高集積化ｉｃ用パッケ−ジおよびその製造方法 |
| JPS6462143A (en) * | 1987-08-31 | 1989-03-08 | Shimadzu Corp | Mri apparatus |
| JP2713530B2 (ja) * | 1992-08-21 | 1998-02-16 | 株式会社日立製作所 | Ａ／ｄ変換データ処理装置 |
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 | Ｍｒ−ｃｔ装置の受信装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPS6211438A (ja) | 1987-01-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP0132358B1 (en) | Imaging method and apparatus for obtaining images using nmr | |
| US10698053B2 (en) | System and method for gradient-modulated sweep imaging with fourier transformation magnetic resonance imaging | |
| WO1987005201A1 (en) | Method of selective excitation in nmr imaging | |
| JPS59166846A (ja) | 磁気共鳴イメージング装置 | |
| JPH0380013B2 (ja) | ||
| EP0213614B1 (en) | Method for obtaining nuclear magnetic resonance information data | |
| JPH04246327A (ja) | ダイナミックｍｒイメージング法 | |
| WO1987002569A1 (en) | Method and apparatus for nmr imaging | |
| JPH10305021A (ja) | 核スピントモグラフィ装置用パルスシーケンス生成方法及び核スピントモグラフィ装置 | |
| JP3274891B2 (ja) | 核磁気共鳴検査装置 | |
| JPH11216126A (ja) | 核磁気共鳴検査装置 | |
| JPH10243934A (ja) | 磁気共鳴診断装置 | |
| JP3104709B2 (ja) | Ｍｒイメージング装置 | |
| JPH0795971A (ja) | Ｍｒイメージング装置 | |
| JP3205061B2 (ja) | 磁気共鳴イメージング装置 | |
| JPS6044880A (ja) | チャ−プトラッキング方式を用いた核磁気共鳴イメ−ジング装置 | |
| JPH0564635A (ja) | Ｍｒｉ装置における水，脂肪スペクトルの分離表示方法 | |
| JPH01274749A (ja) | 核磁気共鳴イメージング装置 | |
| JPH0475638A (ja) | 磁気共鳴イメージング装置 | |
| JP2859264B2 (ja) | 磁気共鳴イメージング装置 | |
| JPH0531096A (ja) | Ｍｒイメージング装置 | |
| JPH04197241A (ja) | 磁気共鳴イメージング装置における血管造影像撮影方法 | |
| JPS61161441A (ja) | Ｎｍｒ断層撮像装置の位相補正方法 | |
| JPH042898B2 (ja) | ||
| Rodrigues | An automatic nuclear quadrupole resonance spectrometer |