JPH02220635A - 基準マーカを使用したｎｍｒ信号位相および振幅補正方式 - Google Patents
基準マーカを使用したｎｍｒ信号位相および振幅補正方式Info
- Publication number
- JPH02220635A JPH02220635A JP1307835A JP30783589A JPH02220635A JP H02220635 A JPH02220635 A JP H02220635A JP 1307835 A JP1307835 A JP 1307835A JP 30783589 A JP30783589 A JP 30783589A JP H02220635 A JPH02220635 A JP H02220635A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- complex
- signals
- frequency
- correction function
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 54
- 239000003550 marker Substances 0.000 title 1
- 238000001914 filtration Methods 0.000 claims abstract description 6
- 238000001228 spectrum Methods 0.000 claims description 38
- 238000000034 method Methods 0.000 claims description 16
- 230000001131 transforming effect Effects 0.000 claims description 7
- 238000012545 processing Methods 0.000 abstract description 21
- 238000012935 Averaging Methods 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 abstract description 2
- 230000010287 polarization Effects 0.000 abstract description 2
- 230000003068 static effect Effects 0.000 abstract description 2
- 230000005540 biological transmission Effects 0.000 abstract 1
- 238000005481 NMR spectroscopy Methods 0.000 description 25
- 238000002474 experimental method Methods 0.000 description 17
- 238000003384 imaging method Methods 0.000 description 8
- 238000012565 NMR experiment Methods 0.000 description 6
- 230000002452 interceptive effect Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 150000001875 compounds Chemical class 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- CSCPPACGZOOCGX-UHFFFAOYSA-N Acetone Chemical compound CC(C)=O CSCPPACGZOOCGX-UHFFFAOYSA-N 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000004611 spectroscopical analysis Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- UYZGNJORQKVKGM-UHFFFAOYSA-N 1,2,3,4-tetramethylsiline Chemical compound CC1=C(C)C(C)=[Si](C)C=C1 UYZGNJORQKVKGM-UHFFFAOYSA-N 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000012306 spectroscopic technique Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000008646 thermal stress Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
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/58—Calibration of imaging systems, e.g. using test probes, Phantoms; Calibration objects or fiducial markers such as active or passive RF coils surrounding an MR active material
- G01R33/583—Calibration of signal excitation or detection systems, e.g. for optimal RF excitation power or frequency
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
発明の背景
本発明はＮＭＲ分光法およびＮＭＲイメージングに関す
るものであり、更に詳しくは相次ぐ実験または走査で取
得される２つのＮＭＲ信号の間の位相および振幅の変動
を補正するための方法に関するものである。
るものであり、更に詳しくは相次ぐ実験または走査で取
得される２つのＮＭＲ信号の間の位相および振幅の変動
を補正するための方法に関するものである。
相次ぐＮＭＲ実験の間に取得される信号の変動の発生源
は主として、印加される分極磁界Ｂ。の強さの変化の結
果である。Ｂｏの変化は振動または熱応力によるＮＭＲ
磁石の物理的変化から生じることがあり、あるいは磁界
上の外部の物体の影響で生じることもある。
は主として、印加される分極磁界Ｂ。の強さの変化の結
果である。Ｂｏの変化は振動または熱応力によるＮＭＲ
磁石の物理的変化から生じることがあり、あるいは磁界
上の外部の物体の影響で生じることもある。
相次ぐＮ Ｍ Ｒ実験相互の間の信号振幅の変動は各実
験の間に取得される数個の信号を平均してその実験に対
する合成信号とする付加的なステップにより生ずる。１
つの実験で取得されるこのような多数の信号の相互の間
の位相変移が平均された信号の振幅変動に変換される。
験の間に取得される数個の信号を平均してその実験に対
する合成信号とする付加的なステップにより生ずる。１
つの実験で取得されるこのような多数の信号の相互の間
の位相変移が平均された信号の振幅変動に変換される。
正確に再生できる分光信号の位相および振幅に対する要
求は、例えば水からの結合されていない陽子の共鳴のよ
うな大きな望ましくない信号が存在している状態である
元素の間の陽子結合から生じる小さな信号を調べること
を必要とするようなりラスの分光技術で生じる。このよ
うな大きな妨害信号の抑圧は相次いで取得される信号が
関心のある信号の位相差を除けば開−であるように実験
を設計することによって行なわれる。これらの相次ぐ信
号の減算によって、実験相互の間で一定に留まっている
望ましくない妨害信号が相殺される。
求は、例えば水からの結合されていない陽子の共鳴のよ
うな大きな望ましくない信号が存在している状態である
元素の間の陽子結合から生じる小さな信号を調べること
を必要とするようなりラスの分光技術で生じる。このよ
うな大きな妨害信号の抑圧は相次いで取得される信号が
関心のある信号の位相差を除けば開−であるように実験
を設計することによって行なわれる。これらの相次ぐ信
号の減算によって、実験相互の間で一定に留まっている
望ましくない妨害信号が相殺される。
位相が変化した関心のある信号は相殺されない。
妨害信号が関心のある信号に比べて数桁大きいことが多
い。したがって、妨害信号の位相および振幅が相次ぐ実
験相互の間で一定に留まっていることが重要であり、そ
うでないと、減算処理の後、望ましくない信号のかなり
大きな成分が相殺されないで残る。
い。したがって、妨害信号の位相および振幅が相次ぐ実
験相互の間で一定に留まっていることが重要であり、そ
うでないと、減算処理の後、望ましくない信号のかなり
大きな成分が相殺されないで残る。
位Ｉｌｌおよび振幅の安定度が重大であるこのような１
つの分光実験が米国特許出願第１８１，９５６号（米国
特許第４，８４３，３２１号）に述べられている。これ
には相次ぐ実験相互の間のこれらの信号の成る成分のみ
の位相を調節することにより望ましくないＮＭＲ信号を
相殺するための多数の技術が述べられている。
つの分光実験が米国特許出願第１８１，９５６号（米国
特許第４，８４３，３２１号）に述べられている。これ
には相次ぐ実験相互の間のこれらの信号の成る成分のみ
の位相を調節することにより望ましくないＮＭＲ信号を
相殺するための多数の技術が述べられている。
ＮＭＲイメージングの分野では、正確な位を口と振幅の
安定性によって、結果として得られる信号の信号対雑音
比が改善され、１つ以上の軸で位相情報によって左右さ
れる画像またはスペクトルの分解能が改善される。長い
取得遅延による「スピンエコー」のようなイメージング
技術は特に位相誤差の影響を受けやすい。
安定性によって、結果として得られる信号の信号対雑音
比が改善され、１つ以上の軸で位相情報によって左右さ
れる画像またはスペクトルの分解能が改善される。長い
取得遅延による「スピンエコー」のようなイメージング
技術は特に位相誤差の影響を受けやすい。
発明の要約
本発明は相次いで取得されたＮＭＲ信号の位相および振
幅の誤差であって、分光またはイメージング装置の不安
定性、特に印加される分極磁界Ｂｏの変動によって生じ
る誤差を補正するための計算手段に関するものである。
幅の誤差であって、分光またはイメージング装置の不安
定性、特に印加される分極磁界Ｂｏの変動によって生じ
る誤差を補正するための計算手段に関するものである。
更に詳しく述べると本発明の一実施例では、フーリエ変
換されたＮＭＲ信号の基準共鳴ピークが分離されて、相
次ぐ実験において位相および振幅の変動を補正するため
の基準として使用される。
換されたＮＭＲ信号の基準共鳴ピークが分離されて、相
次ぐ実験において位相および振幅の変動を補正するため
の基準として使用される。
分離されたピークは逆変換されて、帯域通過濾波された
複素表現の時間領域信号を発生する。相次ぐ実験から得
られるこのような２つの信号の比を用いて、複素補正係
数が作られる。この補正係数はＮＭＲ信号に印加したと
き、基準ピークに対してだけでなく同等の位相および振
幅の誤差がある信号スペクトルの他のピークに対しても
信号の位相および振幅を補正する。
複素表現の時間領域信号を発生する。相次ぐ実験から得
られるこのような２つの信号の比を用いて、複素補正係
数が作られる。この補正係数はＮＭＲ信号に印加したと
き、基準ピークに対してだけでなく同等の位相および振
幅の誤差がある信号スペクトルの他のピークに対しても
信号の位相および振幅を補正する。
代りの実施例では、基準として使用される共鳴ピークは
数学的に求めるか、または慎重に制御された実験により
取得され、相次いで取得される信号の位相を正規化する
ために使用される。
数学的に求めるか、または慎重に制御された実験により
取得され、相次いで取得される信号の位相を正規化する
ために使用される。
本発明の一般的な目的はイメージングまたは分光の用途
で取得される信号の信号対雑音比を改良することである
。信号対雑音比を改良するために多数の信号が平均され
る分光またはイメージングの用途に本発明を用いた場合
、このような信号相互の間の位相誤差のｈａ正によって
、付加的に信号を追加する毎に信号対雑音比を増大して
いくことができる。これは、位相誤差によって、信号が
完全に相加的になることが防止されないからである。
で取得される信号の信号対雑音比を改良することである
。信号対雑音比を改良するために多数の信号が平均され
る分光またはイメージングの用途に本発明を用いた場合
、このような信号相互の間の位相誤差のｈａ正によって
、付加的に信号を追加する毎に信号対雑音比を増大して
いくことができる。これは、位相誤差によって、信号が
完全に相加的になることが防止されないからである。
本発明のもう１つの目的は各実験における数個の信号の
平均化の結果として、相次ぐＮＭＲ実験で取得される信
号相互の間に生じる振幅変動を小さくすることである。
平均化の結果として、相次ぐＮＭＲ実験で取得される信
号相互の間に生じる振幅変動を小さくすることである。
このような平均された信号相互の間の位相誤差は平均さ
れた信号の振幅誤差に変換される。位相誤差が小さくな
ると、この振幅誤差が小さくなる。
れた信号の振幅誤差に変換される。位相誤差が小さくな
ると、この振幅誤差が小さくなる。
本発明の更にもう１つの目的は分光またはイメージング
装置の信号取得チェーンの利１Ｆ４変化から生じる振幅
誤差を小さくすることである。取得されたＮ Ｍ Ｒ信
号の相次ぐ基準ピークの相対振幅は、このような位相か
ら独立した振幅の変動に対する補正を行なう補正係数に
組込まれる。
装置の信号取得チェーンの利１Ｆ４変化から生じる振幅
誤差を小さくすることである。取得されたＮ Ｍ Ｒ信
号の相次ぐ基準ピークの相対振幅は、このような位相か
ら独立した振幅の変動に対する補正を行なう補正係数に
組込まれる。
本発明の付加的な目的はＮＭＲイメージングの用途で幾
何学的な歪を小さくすることである。ＮＭＲ信号の位相
は１つ以上の次元でＮＭＲ画像を規定する。位相誤差の
減少により、ＮＭＲ画像の空間的な正確さが改善される
。
何学的な歪を小さくすることである。ＮＭＲ信号の位相
は１つ以上の次元でＮＭＲ画像を規定する。位相誤差の
減少により、ＮＭＲ画像の空間的な正確さが改善される
。
本発明の上記および他の目的および利点は以下の説明か
ら明らかとなる。以下の説明では添付の図面を参照する
。図面は本発明の実施例を図示したものである。しかし
、このような実施例は必らずしも本発明の全範囲を表わ
すものではないので、本発明の範囲の解釈のためには特
許請求の範囲を参照しなければならない。
ら明らかとなる。以下の説明では添付の図面を参照する
。図面は本発明の実施例を図示したものである。しかし
、このような実施例は必らずしも本発明の全範囲を表わ
すものではないので、本発明の範囲の解釈のためには特
許請求の範囲を参照しなければならない。
発明の詳細な説明
まず第１図には、本発明で使用するのに適したＮＭＲ分
光システムの概略構成図が示されている。
光システムの概略構成図が示されている。
試料ｌはデカルト座標系の正の２軸方向を向いた物質の
静止分極磁界Ｂｏの中に保持されている。
静止分極磁界Ｂｏの中に保持されている。
座標系の原点は試料１の中心にあるものとする。
当業者には理解されるように、送信コイル（図示しない
）からの外部無線周波の場が試料１に印加されることに
より、試料１の核またはスピンが励起されてラーモア周
波数で共鳴する。この無線周波の場は単一または多重の
パルスで構成することができる。次に、励起されたスピ
ンの減衰振動によって生じるＮＭＲ信号が検出できるよ
うに外部無線周波の場が除去される。信号対雑音比を改
善するために、第１図にアンテナ２および３として示さ
れている直角に配置されている２個のループアンテナに
よってこのＮＭＲ無線周波信号が取得される。ループア
ンテナ２はその軸が上記座標系のＸ軸と揃うように配置
される。ループアンテナ３はその軸が座標系のｙ軸と揃
うように配置される；アンテナ２および３からのＮＭＲ
信号はそれぞれ前置増幅器４および５によって増幅され
る。
）からの外部無線周波の場が試料１に印加されることに
より、試料１の核またはスピンが励起されてラーモア周
波数で共鳴する。この無線周波の場は単一または多重の
パルスで構成することができる。次に、励起されたスピ
ンの減衰振動によって生じるＮＭＲ信号が検出できるよ
うに外部無線周波の場が除去される。信号対雑音比を改
善するために、第１図にアンテナ２および３として示さ
れている直角に配置されている２個のループアンテナに
よってこのＮＭＲ無線周波信号が取得される。ループア
ンテナ２はその軸が上記座標系のＸ軸と揃うように配置
される。ループアンテナ３はその軸が座標系のｙ軸と揃
うように配置される；アンテナ２および３からのＮＭＲ
信号はそれぞれ前置増幅器４および５によって増幅され
る。
このように増幅された信号は更に増幅して検出および一
波を行なうために、それぞれ受信器６および７に印加さ
れる。信号のサンプリングとディジタル化を容易にする
ため、受信器６および７は通常のヘテロダイン処理によ
って信号の周波数をメガヘルツ範囲からキロヘルツ範囲
へ変移させる。
波を行なうために、それぞれ受信器６および７に印加さ
れる。信号のサンプリングとディジタル化を容易にする
ため、受信器６および７は通常のヘテロダイン処理によ
って信号の周波数をメガヘルツ範囲からキロヘルツ範囲
へ変移させる。
次に信号はコンピュータ１１による処理のため、平均化
器８および９によってサンプリングされて、ディジタル
化される。コンピュータ１１はフーリエ変換のようなデ
ータ処理、ディスク記憶装置１２と協働してのデータ記
憶、ならびに他の論理および算術演算を行なう。通常、
これらの動作はすべてミニコンピユータによって行なわ
れるので、機能的にしか説明しない。
器８および９によってサンプリングされて、ディジタル
化される。コンピュータ１１はフーリエ変換のようなデ
ータ処理、ディスク記憶装置１２と協働してのデータ記
憶、ならびに他の論理および算術演算を行なう。通常、
これらの動作はすべてミニコンピユータによって行なわ
れるので、機能的にしか説明しない。
後で更に詳しく説明するコンピュータ１１による後続の
計算処理のため、一方のアンテナからの信号は複素無線
周波信号の虚数成分として扱われ、他方のアンテナから
の信号は実数成分として扱われる。したがって、１つの
ＮＭＲ実験で２つの信号が取得される。すなわち、これ
らの信号は実数信号Ａｒ（ｔ）および虚数信号Ａｔ
（ｔ）であり、これらは時間に対する単一の複素関数Ｖ
（ｔ）として扱うことができる。
計算処理のため、一方のアンテナからの信号は複素無線
周波信号の虚数成分として扱われ、他方のアンテナから
の信号は実数成分として扱われる。したがって、１つの
ＮＭＲ実験で２つの信号が取得される。すなわち、これ
らの信号は実数信号Ａｒ（ｔ）および虚数信号Ａｔ
（ｔ）であり、これらは時間に対する単一の複素関数Ｖ
（ｔ）として扱うことができる。
Ｖ （ｔ）＝Ａｒ （ｔ）＋ｊＡ＋ （ｔ）
（１）後で明らかとなる理由のため、Ｖ （ｔ）を単
一の周波数成分すなわち単一の回転ベクトルで近似でき
る場合には、次式が成立することに注意すべきである。
（１）後で明らかとなる理由のため、Ｖ （ｔ）を単
一の周波数成分すなわち単一の回転ベクトルで近似でき
る場合には、次式が成立することに注意すべきである。
Ｖ （ｔ） −Ａｃｏｓ （ωｔ） ＋ ｊＡｓｉｎ
（ωｔ）但し、 ω■２π／Ｔ Ｔ−（単一の周波数成分の周期） 上式のＣＯＳおよびｓｉｎの項は直交アンテナによって
取得される両信号の間の９０度の位相差の結果として生
じる。ｒＡＪはその周波数成分の振幅を表わす。
（ωｔ）但し、 ω■２π／Ｔ Ｔ−（単一の周波数成分の周期） 上式のＣＯＳおよびｓｉｎの項は直交アンテナによって
取得される両信号の間の９０度の位相差の結果として生
じる。ｒＡＪはその周波数成分の振幅を表わす。
上記の通り、信号Ｖ （ｔ）は平均化器８および９によ
ってサンプリングされてディジタル化されるので、１−
０からｉｍｎまでの値Ｖ （ｔｔ ）の配列と考えるこ
とができる。ここで指標値ｔ１ とｔ１＋１との差は平
均化器８および９のサンプリング周期であり、通常、キ
ロヘルツ範囲の信号に対して５０マイクロ秒のオーダー
である。代表的なＮＭＲ信号の取得の場合、サンプル数
ｒｎＪは１６０００のオーダーである。表記法を簡単に
するため、これらのサンプリングされた配列は普通の時
間関数、たとえばＶ （ｔ）と表わすべきである。
ってサンプリングされてディジタル化されるので、１−
０からｉｍｎまでの値Ｖ （ｔｔ ）の配列と考えるこ
とができる。ここで指標値ｔ１ とｔ１＋１との差は平
均化器８および９のサンプリング周期であり、通常、キ
ロヘルツ範囲の信号に対して５０マイクロ秒のオーダー
である。代表的なＮＭＲ信号の取得の場合、サンプル数
ｒｎＪは１６０００のオーダーである。表記法を簡単に
するため、これらのサンプリングされた配列は普通の時
間関数、たとえばＶ （ｔ）と表わすべきである。
これらの関数は実際にはコンピュータ１１内でサンプリ
ングされてディジタル化された配列として操作される。
ングされてディジタル化された配列として操作される。
次に２つの相次ぐＮＭＲ実験および２つの複素信号の取
得を考える。実験１における複素信号は上記のように実
数部と虚数部を含むＡ （ｔ）と名付けられる。実験２
における複素信号はＢ （ｔ）と名付けられる。
得を考える。実験１における複素信号は上記のように実
数部と虚数部を含むＡ （ｔ）と名付けられる。実験２
における複素信号はＢ （ｔ）と名付けられる。
第５図はコンピュータ１１により遂行する振幅補正プロ
グラム１９のブロック図を示す。このプログラムの第１
のステップは、処理ブロック２０で示すように、そして
次式で表わされるように、取得した各信号をフーリエ変
換することで構成される。
グラム１９のブロック図を示す。このプログラムの第１
のステップは、処理ブロック２０で示すように、そして
次式で表わされるように、取得した各信号をフーリエ変
換することで構成される。
Ｃ（ｆ）−デ［Ａ （ｔ）１ （３）
Ｄ （ｆ）−Ｆ［Ｂ（ｔ）］ （４）
このフーリエ変換は当業者に知られている多数の高速フ
ーリエ変換アルゴリズムのうちの１つで行なうことがで
きる。たとえば、アール・プレースウェル（Ｒｊ）ｒａ
ｃｅｗｅｌｌ ）の著書［高速フーリエ変換とその応用
（Ｔｈｅ Ｐｏｕｒｌｅｒ Ｔｒａｎｓｌ’ｏｒｍ ａ
ｎｄ Ｉｔｓ Ａｐｐｌｉｃａｔｉｏｎ ） Ｊ 、マ
グロ−ヒル発行（１９６５）、を参照されたい。
Ｄ （ｆ）−Ｆ［Ｂ（ｔ）］ （４）
このフーリエ変換は当業者に知られている多数の高速フ
ーリエ変換アルゴリズムのうちの１つで行なうことがで
きる。たとえば、アール・プレースウェル（Ｒｊ）ｒａ
ｃｅｗｅｌｌ ）の著書［高速フーリエ変換とその応用
（Ｔｈｅ Ｐｏｕｒｌｅｒ Ｔｒａｎｓｌ’ｏｒｍ ａ
ｎｄ Ｉｔｓ Ａｐｐｌｉｃａｔｉｏｎ ） Ｊ 、マ
グロ−ヒル発行（１９６５）、を参照されたい。
一般にフーリエ変換は変換される信号の周波数成分の振
幅を周波数の関数として示すと共に、その信号の周波数
成分の位相を周波数の関数として示す一組の組を作る。
幅を周波数の関数として示すと共に、その信号の周波数
成分の位相を周波数の関数として示す一組の組を作る。
フーリエ変換の代りの具現法は一組の実数値と虚数値を
作るものであり、この実数値と虚数値は一緒になってス
ペクトルの各周波数に対するベクトルの向きと長さを定
める。
作るものであり、この実数値と虚数値は一緒になってス
ペクトルの各周波数に対するベクトルの向きと長さを定
める。
これらの２つの結果は数学的には等価であり、本明細書
での説明のためにはフーリエ変換の後者の具現法を採用
するものとする。
での説明のためにはフーリエ変換の後者の具現法を採用
するものとする。
第２Ａ図はフーリエ変換処理前のＮＭＲ信号の代表的な
実数部を示す。第２Ａ図の水平軸は時間である。第２Ａ
図の波形は上記のように外部ＲＦの場による励起後の試
料内の多数の原子核の歳差運動から生ずる多重周波数減
衰信号によって作られるＮＭＲ信号を表わそうとしたも
のである。
実数部を示す。第２Ａ図の水平軸は時間である。第２Ａ
図の波形は上記のように外部ＲＦの場による励起後の試
料内の多数の原子核の歳差運動から生ずる多重周波数減
衰信号によって作られるＮＭＲ信号を表わそうとしたも
のである。
第３Ａ図は第２Ａ図の信号のフーリエ変換から生ずる代
表的な複素周波スペクトルの実数部を示し、２つの、す
なわち１つは正に向い且つ他の１つは負に向うピーク１
５と１群の妨害ピーク１６を有する。水平軸は変換の結
果として周波数である。ピーク１５は分光分析の目的の
ための関心のある信号であって、調べている化合物の共
鳴周波数ｆ１に近い信号を概略表示したものである。ピ
ーク１６は妨害信号であり、この信号は本明細書で述べ
る位相補正処理の後に適用される周知の多数の位相変移
（移相）技術のう、ちの１つによって最終的に相殺すべ
きものである。わかりやすくするためにピーク１５は誇
張して示されており、通常は、ピーク１５はピーク１６
によって不明瞭になっていて、明確に区別できない。ピ
ーク１７は基準ピークすなわち基準周波数成分であり、
試料１の中に含まれている化合物から自然に生じるもの
か或いは特定の物質が試料の近くに含有されていること
によって生じるものであってよい。「マカ」と名付けら
れたこの物質は、周波数スペクトルの中に関心のある信
号の周波数から離れて鋭い明確なピークを持っているこ
とが好ましい。アセトン、テトラメチルシリンのような
種々の物質が基準ピークを発生するのに適している。本
発明の実施例では、この目的のため水が用いられる。
表的な複素周波スペクトルの実数部を示し、２つの、す
なわち１つは正に向い且つ他の１つは負に向うピーク１
５と１群の妨害ピーク１６を有する。水平軸は変換の結
果として周波数である。ピーク１５は分光分析の目的の
ための関心のある信号であって、調べている化合物の共
鳴周波数ｆ１に近い信号を概略表示したものである。ピ
ーク１６は妨害信号であり、この信号は本明細書で述べ
る位相補正処理の後に適用される周知の多数の位相変移
（移相）技術のう、ちの１つによって最終的に相殺すべ
きものである。わかりやすくするためにピーク１５は誇
張して示されており、通常は、ピーク１５はピーク１６
によって不明瞭になっていて、明確に区別できない。ピ
ーク１７は基準ピークすなわち基準周波数成分であり、
試料１の中に含まれている化合物から自然に生じるもの
か或いは特定の物質が試料の近くに含有されていること
によって生じるものであってよい。「マカ」と名付けら
れたこの物質は、周波数スペクトルの中に関心のある信
号の周波数から離れて鋭い明確なピークを持っているこ
とが好ましい。アセトン、テトラメチルシリンのような
種々の物質が基準ピークを発生するのに適している。本
発明の実施例では、この目的のため水が用いられる。
再び第５図を参照して説明すると、処理ブロック２１で
示される補正プログラムの次のステップは変換された各
信号に、点毎に、ｆ、に中心があり次式で定義される窓
関数Ｕ （ｆ）を乗算することで構成される。
示される補正プログラムの次のステップは変換された各
信号に、点毎に、ｆ、に中心があり次式で定義される窓
関数Ｕ （ｆ）を乗算することで構成される。
Δｆは関心のあるＮＭＲ信号１５のどの部分も含まず、
かつスペクトル周波数の中の異なる周数の他のピークを
含むことなく、ｆ、の基準ピーク１７のまわりに対称な
窓を作るようにできる限り大きく選定される。ここで注
意しなければならない点は、ｆｏおよびｆｌがヘテロゲ
イン変換前の化合物の共鳴周波数を指すということであ
る。第３Ｂ図は窓関数Ｕ （ｆ）を示し、第３Ｃ図は次
式で示すように変換されたＮＭＲ信号に窓関数を印加し
た後の代表的なスペクトルを示す。
かつスペクトル周波数の中の異なる周数の他のピークを
含むことなく、ｆ、の基準ピーク１７のまわりに対称な
窓を作るようにできる限り大きく選定される。ここで注
意しなければならない点は、ｆｏおよびｆｌがヘテロゲ
イン変換前の化合物の共鳴周波数を指すということであ
る。第３Ｂ図は窓関数Ｕ （ｆ）を示し、第３Ｃ図は次
式で示すように変換されたＮＭＲ信号に窓関数を印加し
た後の代表的なスペクトルを示す。
Ｅ （ｆ） −Ｃ＜ｆ＞ ・Ｕ （ｆ）
（６）Ｆ （ｆ）−Ｄ （ｆ） ・Ｕ
（ｆ） （７）この処理の効果はｆ
ｏの基準ピークの近くの複素周波数スペクトルの一部に
窓を設けることにより基準ピークを分離することである
。これはｆ。
（６）Ｆ （ｆ）−Ｄ （ｆ） ・Ｕ
（ｆ） （７）この処理の効果はｆ
ｏの基準ピークの近くの複素周波数スペクトルの一部に
窓を設けることにより基準ピークを分離することである
。これはｆ。
を中心とするこの窓の外側のすべての点におけるスペク
トルの値をゼロに設定することにより行なわれる。上記
の式（５）で与えられる窓関数は簡単である利点はある
が、周波数（ｆｏ−Δｆ）および（ｆｏ＋Δｆ）の所で
Ｕ （ｆ）が不連続であるので切捨てによるアーティフ
ァクトが生ずる。
トルの値をゼロに設定することにより行なわれる。上記
の式（５）で与えられる窓関数は簡単である利点はある
が、周波数（ｆｏ−Δｆ）および（ｆｏ＋Δｆ）の所で
Ｕ （ｆ）が不連続であるので切捨てによるアーティフ
ァクトが生ずる。
このような切捨て誤差を避けるため、正弦曲線またはガ
ウス曲線に基く窓関数のような他の窓関数を選択しても
よい。窓はスペクトルＥ （ｆ）およびＦ （ｆ）に対
して同じである限り、任意の周波数区間の中で一定であ
る必要はない。
ウス曲線に基く窓関数のような他の窓関数を選択しても
よい。窓はスペクトルＥ （ｆ）およびＦ （ｆ）に対
して同じである限り、任意の周波数区間の中で一定であ
る必要はない。
再び第５図を参照すると、処理ブロック２２に示される
次のステップは窓を通したスペクトルを逆フーリエ変換
することである。その結果得られる信号は次式で表わさ
れる。
次のステップは窓を通したスペクトルを逆フーリエ変換
することである。その結果得られる信号は次式で表わさ
れる。
Ｇ （ｔ） −Ｅ″−’ ［Ｅ （ｆ）］
（８）Ｈ（ｔ）−Ｊ’−’ ［Ｆ
（ｆ）］ （９）このステップの結果得
られる代表的な信号が第２Ｂ図に示されている。第２Ｂ
図の信号は信号Ａ（１）またはＢ （ｔ）に帯域通過濾
波動作を加えることによって作成される信号と同等であ
る。窓Δｆが次第に小さくなるにつれて、信号Ｇ（ｔ）
およびＨ（ｔ）は単一の信号成分しか持たないと考えて
よく、式（２）に関して次のように表わすことができる
。
（８）Ｈ（ｔ）−Ｊ’−’ ［Ｆ
（ｆ）］ （９）このステップの結果得
られる代表的な信号が第２Ｂ図に示されている。第２Ｂ
図の信号は信号Ａ（１）またはＢ （ｔ）に帯域通過濾
波動作を加えることによって作成される信号と同等であ
る。窓Δｆが次第に小さくなるにつれて、信号Ｇ（ｔ）
およびＨ（ｔ）は単一の信号成分しか持たないと考えて
よく、式（２）に関して次のように表わすことができる
。
Ｇ （ｔ） ＝Ａ＋ ｃｏｓ （ωｔ） ＋Ａ＋
ｊｓｉｎ （ωｔ）Ｈ（ｔ） −Ａ２ ｃｏｓ
（ωｔ＋φ）＋Ａ２ ｊｓｌｎ （ωｔ＋φ）
（１１）但し、φは信号Ｇ （ｔ）とＨ（ｔ）との
間の位相差すなわち誤差である。
ｊｓｉｎ （ωｔ）Ｈ（ｔ） −Ａ２ ｃｏｓ
（ωｔ＋φ）＋Ａ２ ｊｓｌｎ （ωｔ＋φ）
（１１）但し、φは信号Ｇ （ｔ）とＨ（ｔ）との
間の位相差すなわち誤差である。
オイラー（Ｅｕｌｅｒ ）の式により、これらの波形は
次のように表わすことができる。
次のように表わすことができる。
Ｇ（ｔ）−Ａ１ ｊ″１ｔ （１２）Ｈ（
ｔ） 、ｐ、２ｊ （″）ｔ＋φ） （１３）再
び第５図を参照すると、処理ブロック２３で示される次
のステップは補正係数Ｊ （ｔ）を決定するものであり
、その詳細は第６図に示されている。処理ブロック２６
で示すように、まずＧ （ｔ）およびＨ（ｔ）がゼロに
等しいか否か試験する。
ｔ） 、ｐ、２ｊ （″）ｔ＋φ） （１３）再
び第５図を参照すると、処理ブロック２３で示される次
のステップは補正係数Ｊ （ｔ）を決定するものであり
、その詳細は第６図に示されている。処理ブロック２６
で示すように、まずＧ （ｔ）およびＨ（ｔ）がゼロに
等しいか否か試験する。
両方ともゼロでなければ、処理ブロック２７で、次式で
示されるような点毎の複素除算により中間の振幅および
位相の補正係数１ （ｔ）を求める。
示されるような点毎の複素除算により中間の振幅および
位相の補正係数１ （ｔ）を求める。
１ （ｔ） −Ｇ （ｔ） ／Ｈ（ｔ） （
１４）１ （ｔ） ＝Ｍ （ｔ） ｅ ｊφ（ｔ）
（１５）ここで、Ｍ （ｔ）はこの複素関数のベクト
ル表現の大きさであり、φ（１）はその角度である。
１４）１ （ｔ） ＝Ｍ （ｔ） ｅ ｊφ（ｔ）
（１５）ここで、Ｍ （ｔ）はこの複素関数のベクト
ル表現の大きさであり、φ（１）はその角度である。
関数１ （ｔ）は角度と位相が各サンプル時点に於け
るＧ （ｔ）と１（（ｔ）との間の位相と振幅の誤差を
表わすベクトルと考えてもよい。誤差関数の表現が第４
Ａ図にベクトル形式で示されている。
るＧ （ｔ）と１（（ｔ）との間の位相と振幅の誤差を
表わすベクトルと考えてもよい。誤差関数の表現が第４
Ａ図にベクトル形式で示されている。
第３Ａ図の基準ピーク１７に適用可能な上記の位相補正
は、ピーク１７．１６および１５の間の（すなわちｆｏ
と【１との間の）周波数差が平均周波数（ｆ＋ ＋ｆｏ
）／２に比べて小さい、すなわち ｆ＋−ｆ。
は、ピーク１７．１６および１５の間の（すなわちｆｏ
と【１との間の）周波数差が平均周波数（ｆ＋ ＋ｆｏ
）／２に比べて小さい、すなわち ｆ＋−ｆ。
＜＜ １ （１６）
（ｆ＋ 十ｆ。）／２
であるとき、第３Ａ図のピーク１５．１６にも適用可能
である。
である。
典型的なＭＲ大実験は、上記の値は百方分率（ｐｐｍ
）で表わすと数ｐｐ献のオーダーであり、式（１６）の
条件を満足する。したがって、位相および振幅の補正値
φ（１）およびＭ （ｔ）も関心のあるＮＭＲ信号を正
確に補正する。
）で表わすと数ｐｐ献のオーダーであり、式（１６）の
条件を満足する。したがって、位相および振幅の補正値
φ（１）およびＭ （ｔ）も関心のあるＮＭＲ信号を正
確に補正する。
分光装置の分解能限界、ディジタル化の処理、および元
の信号の中の雑音の存在によって、」１記で計算したよ
うに誤差関数１ （ｔ）に偶発的な不連続が生じるこ
とがある。これらの不連続は次の２つの条件のうちの１
つが生じているか試験することにより検出することがで
きる。すなわち（イ）補正係数の竹田Ｈ（ｔ）または分
子Ｇ （ｔ）の値がゼロである場合、あるいは（ロ）誤
差関数の大きさ１１（ｔ）１が第１の限界より大きいか
、または第２の限界より小さい場合である。実施例では
、第１の限界は０．２５であり、第２の限界は４である
。
の信号の中の雑音の存在によって、」１記で計算したよ
うに誤差関数１ （ｔ）に偶発的な不連続が生じるこ
とがある。これらの不連続は次の２つの条件のうちの１
つが生じているか試験することにより検出することがで
きる。すなわち（イ）補正係数の竹田Ｈ（ｔ）または分
子Ｇ （ｔ）の値がゼロである場合、あるいは（ロ）誤
差関数の大きさ１１（ｔ）１が第１の限界より大きいか
、または第２の限界より小さい場合である。実施例では
、第１の限界は０．２５であり、第２の限界は４である
。
この第１の条件（イ）は前に処理ブロック２６で試験さ
れており、第２の条件（ロ）は処理ブロック２８で試験
される。これらの条件のいずれかが真であれば、処理ブ
ロック３０でその時間の引数ｔに対する誤差関数１
（ｔ）は（１＋Ｏｊ）、すなわち（実数軸に沿った単位
ベクトル）に等しくなるように設定される。第４Ｂ図は
この置き換えられた補正係数のベクトル表現を示す。こ
の試験および補正（必要な場合）の後の補正係数を以後
Ｊ （ｔ）と呼ぶことにする。
れており、第２の条件（ロ）は処理ブロック２８で試験
される。これらの条件のいずれかが真であれば、処理ブ
ロック３０でその時間の引数ｔに対する誤差関数１
（ｔ）は（１＋Ｏｊ）、すなわち（実数軸に沿った単位
ベクトル）に等しくなるように設定される。第４Ｂ図は
この置き換えられた補正係数のベクトル表現を示す。こ
の試験および補正（必要な場合）の後の補正係数を以後
Ｊ （ｔ）と呼ぶことにする。
与えられた時間の引数に於ける誤差関数をこの値（１＋
ｏｊ）に置き換えることの効果は、修正された誤差関数
Ｊ （ｔ）を補正すべき信号に乗算したときに信号の振
幅または位相に変化が生じないということである。
ｏｊ）に置き換えることの効果は、修正された誤差関数
Ｊ （ｔ）を補正すべき信号に乗算したときに信号の振
幅または位相に変化が生じないということである。
要約すると、誤差関数の分子Ｇ（ｔ）および竹田Ｈ（ｔ
）のどちらもゼロに等しくない場合、Ｇ （ｔ）のＨ（
ｔ）による複素除算が処理ブロック２７に示すように行
なわれ、値Ｉ （ｔ）が求められる。判定ブロック２８
でこの値１ （ｔ）が０゜２５と４との間にあるか否か
を判定する。Ｉ （ｔ）がこれらの値の間にない場合、
Ｊ （ｔ）は再び（１＋Ｏｊ）に等しくなるように設定
される。
）のどちらもゼロに等しくない場合、Ｇ （ｔ）のＨ（
ｔ）による複素除算が処理ブロック２７に示すように行
なわれ、値Ｉ （ｔ）が求められる。判定ブロック２８
でこの値１ （ｔ）が０゜２５と４との間にあるか否か
を判定する。Ｉ （ｔ）がこれらの値の間にない場合、
Ｊ （ｔ）は再び（１＋Ｏｊ）に等しくなるように設定
される。
１ （ｔ）が０．２５と４との間にある場合は、処理
ブロック２９に示すようにＪ （ｔ）は！ （ｔ）に等
しくなるように設定される。この補正処理は判定ブロッ
ク３１で配列全体が処理されたと判定されるまでｔの各
指標値に対して繰返えされる。
ブロック２９に示すようにＪ （ｔ）は！ （ｔ）に等
しくなるように設定される。この補正処理は判定ブロッ
ク３１で配列全体が処理されたと判定されるまでｔの各
指標値に対して繰返えされる。
再び第５図を参照して説明すると、次に、処理ブロック
２４および次式で示されるように点毎の複素乗算により
、修正された誤差関数Ｊ （ｔ）が元のＮＭＲ信号Ｂ
（ｔ）に適用される。
２４および次式で示されるように点毎の複素乗算により
、修正された誤差関数Ｊ （ｔ）が元のＮＭＲ信号Ｂ
（ｔ）に適用される。
Ｂ’ （ｔ） −Ｂ （ｔ） ・Ｊ （ｔ）
（１７）Ｂ （ｔ）にＪ （ｔ）を乗算することの
効果は各サンプル点でＢ （ｔ）の振幅をＭ （ｔ）だ
け変え、またＢ （ｔ）の位相をφ（１）だけ変えるこ
とである。
（１７）Ｂ （ｔ）にＪ （ｔ）を乗算することの
効果は各サンプル点でＢ （ｔ）の振幅をＭ （ｔ）だ
け変え、またＢ （ｔ）の位相をφ（１）だけ変えるこ
とである。
第２Ｃ図は補正後の信号Ｂ’ （ｔ）を実線で、補正
前の信号Ｂ （ｔ）を破線で示す。Ｂ’ （ｔ）の位
相φおよび振幅Ｍが調整されていることに注目すべきで
ある。
前の信号Ｂ （ｔ）を破線で示す。Ｂ’ （ｔ）の位
相φおよび振幅Ｍが調整されていることに注目すべきで
ある。
上記の処理は２つの相次ぐＮＭＲ実験または走査からの
信号に対して遂行される。しかし、この処理を任意の数
の相次ぐ信号に対して各対ごとに遂行することにより、
これらの信号の振幅および位相を補正することができる
ことは明らかである。
信号に対して遂行される。しかし、この処理を任意の数
の相次ぐ信号に対して各対ごとに遂行することにより、
これらの信号の振幅および位相を補正することができる
ことは明らかである。
また、比較を行なう信号はデータ取得の目的で取得する
必要はなく、単に所望の特性のテンプレート（ｔｅｌｐ
ｌａｔｅ）信号であってもよい。このテンプレート信号
は制御された条件下で慎重な実験によって作成してもよ
いし、あるいは簡単な場合には数学的に導き出してもよ
い。
必要はなく、単に所望の特性のテンプレート（ｔｅｌｐ
ｌａｔｅ）信号であってもよい。このテンプレート信号
は制御された条件下で慎重な実験によって作成してもよ
いし、あるいは簡単な場合には数学的に導き出してもよ
い。
再び第５図を参照すると、補正方法の一部ではないが、
次の処理ブロック２５は信号Ａ （ｔ）および補正され
た信号Ｂ’ （ｔ）の代表的な使用例を示す。詳しく
説明すると、信号Ｂ’ （ｔ）が信号Ａ （ｔ）から
減算されて、各信号の同相である部分は相殺され、各信
号の同相でない部分は元のままになる。この減算の結果
束じた信号の代表的なフーリエ変換が第３Ｄ図に示され
ている。図かられかるようにピーク１５が顕著になる。
次の処理ブロック２５は信号Ａ （ｔ）および補正され
た信号Ｂ’ （ｔ）の代表的な使用例を示す。詳しく
説明すると、信号Ｂ’ （ｔ）が信号Ａ （ｔ）から
減算されて、各信号の同相である部分は相殺され、各信
号の同相でない部分は元のままになる。この減算の結果
束じた信号の代表的なフーリエ変換が第３Ｄ図に示され
ている。図かられかるようにピーク１５が顕著になる。
もちろん、上記の相殺は信号Ａ （ｔ）と信号Ｂ （ｔ
）との間に適当な位＋１＋変移を生じさせる実験を必要
とする。その説明は本発明の範囲外である。
）との間に適当な位＋１＋変移を生じさせる実験を必要
とする。その説明は本発明の範囲外である。
第７図には２つのスペクトルが示されており、各スペク
トルは処理ブロック２５で行なわれるような２つの信号
の組み合わせである。スペクトル３１は本発明による位
相および振幅の補正技術を適用しない場合の信号の組合
わせを示す。スペクトル３２は本発明による位相および
振幅の補正係数を適用した場合の同じ信号の組合わせを
示す。
トルは処理ブロック２５で行なわれるような２つの信号
の組み合わせである。スペクトル３１は本発明による位
相および振幅の補正技術を適用しない場合の信号の組合
わせを示す。スペクトル３２は本発明による位相および
振幅の補正係数を適用した場合の同じ信号の組合わせを
示す。
Claims (1)
- 【特許請求の範囲】 １、第１および第２の複素ＮＭＲ信号Ａ（ｔ）およびＢ
（ｔ）の間の位相および振幅の誤差を補正する方法であ
つて、 上記の両方のＮＭＲ信号に存在する第１の基準周波数成
分を識別するステップ、 中心周波数が上記第１の基準周波数成分の周波数とほぼ
同じである帯域通過フィルタで、上記第１および第２の
複素ＮＭＲ信号を帯域通過濾波することにより、対応す
る第１および第２の濾波後の複素信号Ｇ（ｔ）およびＨ
（ｔ）を発生するステップ、 上記第１の濾波後の複素信号Ｇ（ｔ）を上記第２の濾波
後の複素信号Ｈ（ｔ）で除算することにより、補正関数
Ｊ（ｔ）を計算するステップ、上記第２の複素ＮＭＲ信
号Ｂ（ｔ）に上記補正関数Ｊ（ｔ）を適用することによ
り、位相および振幅を調節した補正済みの第２の複素Ｎ
ＭＲ信号Ｂ′（ｔ）を発生するステップ、 を含むことを特徴とする補正方法。 ２、上記帯域通過濾波するステップが、 上記第１および第２の複素ＮＭＲ信号Ａ（ｔ）およびＢ
（ｔ）をフーリエ変換することにより第１および第２の
複素周波数スペクトルＣ（ｆ）およびＤ（ｆ）を作成す
るステップ、 上記基準周波数以外のすべての周波数で上記第１および
第２の複素周波数スペクトルＣ（ｆ）およびＤ（ｆ）を
ゼロに設定することにより、対応する第１および第２の
修正周波数スペクトルＥ（ｆ）およびＦ（ｆ）を作成す
るステップ、上記第１および第２の修正周波数スペクト
ルＥ（ｔ）およびＦ（ｆ）に対して逆フーリエ変換を行
なうことにより、対応する上記第１および第２の濾波後
の複素信号Ｇ（ｔ）およびＨ（ｔ）を作成するステップ
、で構成されている請求項１記載の補正方法。 ３、上記補正関数Ｊ（ｔ）を計算するステップが、 上記第１の濾波後の複素信号Ｇ（ｔ）を上記第２の濾波
後の複素信号Ｈ（ｔ）で除算することにより中間補正関
数Ｉ（ｔ）を作成するステップ、上記中間補正関数Ｉ（
ｔ）が第１の限界値を超えるかまたは第２の限界値より
小さいとき、上記補正関数Ｊ（ｔ）を所定のデフォルト
（ｄｅｆａｕｌｔ）値に設定するステップ、および 上記中間補正関数Ｉ（ｔ）の大きさが上記第１の限界値
と上記第２の限界値との間にあるとき、上記補正関数Ｊ
（ｔ）を上記中間補正関数Ｉ（ｔ）に等しく設定するス
テップ、で構成されている請求項１記載の補正方法。 ４、第１および第２の複素ＮＭＲ信号Ａ（ｔ）とＢ（ｔ
）との間の位相および振幅の誤差を補正する方法であっ
て、 上記複素信号Ａ（ｔ）およびＢ（ｔ）をフーリエ変換す
ることにより、周波数スペクトルＣ（ｔ）およびＤ（ｆ
）を作成するステップ、 各スペクトルＣ（ｆ）およびＤ（ｆ）において周波数ｆ
＿０の基準ピークを識別するステップ、上記周波数スペ
クトルＣ（ｆ）およびＤ（ｆ）に、ｆ＿０を中心とする
窓関数Ｕ（ｆ）を乗算することにより、スペクトルＥ（
ｆ）およびＦ（ｆ）を作成するステップ、 上記スペクトルＥ（ｆ）およびＦ（ｆ）を逆フーリエ変
換することにより、信号Ｇ（ｔ）およびＨ（ｔ）を作成
するステップ、 上記Ｇ（ｔ）を上記信号Ｈ（ｔ）で除算することにより
、補正関数Ｉ（ｔ）を作成するステップ、上記第２の複
素ＮＭＲ信号Ｂ（ｔ）に上記補正関数Ｉ（ｔ）を乗算す
ることにより、補正された複素ＮＭＲ信号Ｂ′（ｔ）を
作成するステップ、を含むことを特徴とする補正方法。 ５、第１および第２の複素ＮＭＲ信号Ａ（ｔ）とＢ（ｔ
）の間の位相および振幅の誤差を補正する方法であって
、 上記複素ＮＭＲ信号Ａ（ｔ）およびＢ（ｔ）をフーリエ
変換することにより、周波数スペクトルＣ（ｆ）および
Ｄ（ｆ）を作成するステップ、各周波数スペクトルＣ（
ｆ）およびＤ（ｆ）における周波数ｆ＿０の基準ピーク
を識別するステップ、 上記周波数スペクトルＣ（ｆ）およびＤ（ｆ）にｆ＿０
を中心とする窓関数Ｕ（ｆ）を乗算することにより、ス
ペクトルＥ（ｆ）およびＦ（ｆ）を作成するステップ、 上記信号Ｅ（ｆ）およびＦ（ｆ）を逆フーリエ変換する
ことにより、信号Ｇ（ｔ）およびＨ（ｔ）を作成するス
テップ、 上記信号Ｇ（ｔ）を上記Ｈ（ｔ）で除算することにより
、補正関数Ｉ（ｔ）を作成するステップ、上記補正関数
Ｉ（ｔ）が第１の限界値より大きいか、または第２の限
界値より小さいとき、引数ｔに対するＩ（ｔ）の値を所
定のデフォルト値に置き換えることにより、修正補正関
数Ｊ（ｔ）を作成するステップ、 上記複素ＮＭＲ信号Ｂ（ｔ）に上記修正補正関数Ｊ（ｔ
）を乗算することにより、位相および振幅を補正した複
素ＮＭＲ信号Ｂ′（ｔ）を作成するステップ、 を含むことを特徴とする補正方法。 ６、第１および第２の複素ＮＭＲ信号Ａ（ｔ）およびＢ
（ｔ）の間の位相および振幅の誤差を補正する装置であ
って、 上記の両方のＮＭＲ信号に存在する第１の基準周波数成
分を識別する手段と、 中心周波数が上記第１の基準周波数成分の周波数とほぼ
同じである帯域通過フィルタを用いて、上記第１および
第２の複素ＮＭＲ信号を帯域通過濾波することにより、
対応する第１および第２の濾波後の複素信号Ｇ（ｔ）お
よびＨ（ｔ）を発生する手段と、 上記第１の濾波後の複素信号Ｇ（ｔ）を上記第２の濾波
後の複素信号Ｈ（ｔ）で除算することにより、補正関数
Ｊ（ｔ）を計算する手段と、上記第２の複素ＮＭＲ信号
Ｂ（ｔ）に上記補正関数Ｊ（ｔ）を適用することにより
、位相および振幅を調節した補正済みの第２の複素ＮＭ
Ｒ信号Ｂ′（ｔ）を発生する手段と、を含むことを特徴
とする補正装置。 ７、上記帯域通過濾波する手段が、 上記第１および第２の複素ＮＭＲ信号Ａ（ｔ）およびＢ
（ｔ）をフーリエ変換することにより第１および第２の
複素周波数スペクトルＣ（ｆ）およびＤ（ｆ）を作成す
る手段と、 上記基準周波数以外のすべての周波数で上記第１および
第２の複素周波数スペクトルＣ（ｆ）およびＤ（ｆ）を
ゼロに設定することにより、対応する第１および第２の
修正周波数スペクトルＥ（ｆ）およびＦ（ｆ）を作成す
る手段と、上記第１および第２の修正周波数スペクトル
Ｅ（ｆ）およびＦ（ｆ）に対して逆フーリエ変換を行な
うことにより、対応する上記第１および第２の濾波後の
複素信号Ｇ（ｔ）およびＨ（ｔ）を作成する手段と、で
構成されている請求項１記載の補正方法。 ８、上記補正関数Ｊ（ｔ）を計算する手段が、上記第１
の濾波後の複素信号Ｇ（ｔ）を上記第２の濾波後の複素
信号Ｈ（ｔ）で除算することにより中間補正関数Ｉ（ｔ
）を作成する手段と、上記中間補正関数Ｉ（ｔ）が第１
の限界値を超えるかまたは第２の限界値より小さいとき
、上記補正関数Ｊ（ｔ）を所定のデフォルト（ｄｅｆａ
ｕｌｔ）値に設定する手段と、 上記中間補正関数Ｉ（ｔ）の大きさが上記第１の限界値
と上記第２の限界値との間にあるとき、上記補正関数Ｊ
（ｔ）を上記中間補正関数Ｉ（ｔ）に等しく設定する手
段と、で構成されている請求項１記載の補正装置。 ９、第１および第２の複素ＮＭＲ信号Ａ（ｔ）とＢ（ｔ
）との間の位相および振幅の誤差を補正する装置であっ
て、 上記複素信号Ａ（ｔ）およびＢ（ｔ）をフーリエ変換す
ることにより、周波数スペクトルＣ（ｆ）およびＤ（ｆ
）を作成する手段と、 各スペクトルＣ（ｆ）およびＤ（ｆ）において周波数ｆ
＿０の基準ピークを識別する手段と、上記周波数スペク
トルＣ（ｔ）およびＤ（ｆ）に、ｆ＿０を中心とする窓
関数Ｕ（ｆ）を乗算することにより、スペクトルＥ（ｆ
）およびＦ（ｆ）を作成する手段と、 上記スペクトルＥ（ｆ）およびＦ（ｆ）を逆フーリエ変
換することにより、信号Ｇ（ｔ）およびＨ（ｔ）を作成
する手段と、 上記信号Ｇ（ｔ）を上記信号Ｈ（ｔ）で除算することに
より、補正関数Ｉ（ｔ）を作成する手段と、 上記第２の複素ＮＭＲ信号Ｂ（ｔ）に上記補正関数Ｉ（
ｔ）を乗算することにより、補正された複素ＮＭＲ信号
Ｂ′（ｔ）を作成する手段と、を含むことを特徴とする
補正装置。 １０、第１および第２の複素ＮＭＲ信号Ａ（ｔ）とＢ（
ｔ）の間の位相および振幅の誤差を補正する装置であっ
て、 上記複素ＮＭＲ信号Ａ（ｔ）およびＢ（ｔ）をフーリエ
変換することにより、周波数スペクトルＣ（ｆ）および
Ｄ（ｆ）を作成する手段と、各周波数スペクトルＣ（ｆ
）およびＤ（ｆ）における周波数ｆ＿０の基準ピークを
識別する手段と、上記周波数スペクトルＣ（ｆ）および
Ｄ（ｆ）に、ｆ＿０を中心とする窓関数Ｕ（ｆ）を乗算
することにより、スペクトルＥ（ｆ）およびＦ（ｆ）を
作成する手段と、 上記信号Ｅ（ｆ）およびＦ（ｆ）を逆フーリエ変換する
ことにより、信号Ｇ（ｔ）およびＨ（ｔ）を作成する手
段と、 上記信号Ｇ（ｔ）を上記信号Ｈ（ｔ）で除算することに
より、補正関数Ｉ（ｔ）を作成する手段と、 上記補正関数Ｉ（ｔ）が第１の限界値より大きいか、ま
たは第２の限界値より小さいとき、引数ｔに対するＩ（
ｔ）の値を所定のデフォルト値に置き換えることにより
、修正補正関数Ｊ（ｔ）を作成する手段と、 上記複素ＮＭＲ信号Ｂ（ｔ）に上記修正補正関数Ｊ（ｔ
）を乗算することにより、位相および振幅を補正した複
素ＮＭＲ信号Ｂ′（ｔ）を作成する手段と、を含むこと
を特徴とする補正装置。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US07/278,121 US4885549A (en) | 1988-11-30 | 1988-11-30 | Method of phase and amplitude correction of NMR signals using a reference marker |
| US278121 | 1988-11-30 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH02220635A true JPH02220635A (ja) | 1990-09-03 |
| JPH0422575B2 JPH0422575B2 (ja) | 1992-04-17 |
Family
ID=23063758
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP1307835A Granted JPH02220635A (ja) | 1988-11-30 | 1989-11-29 | 基準マーカを使用したｎｍｒ信号位相および振幅補正方式 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US4885549A (ja) |
| EP (1) | EP0371349B1 (ja) |
| JP (1) | JPH02220635A (ja) |
| DE (1) | DE68925620T2 (ja) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4970457A (en) * | 1989-04-05 | 1990-11-13 | The Regents Of The University Of California | MRI compensated for spurious rapid variations in static magnetic field during a single MRI sequence |
| FR2735865B1 (fr) * | 1995-06-26 | 1997-08-22 | Centre Nat Rech Scient | Dispositif d'analyse quantitative par resonance magnetique nucleaire |
| DE10032345C2 (de) * | 2000-07-04 | 2002-05-29 | Bruker Medical Gmbh | Verfahren zum Korrigieren von Störeinflüssen auf die MR-Signale einer im Messvolumen einer MR-Apparatur angeordneten Substanz sowie MR-Apparatur zur Durchführung des Verfahrens und Rechnereinheit |
| AUPR005900A0 (en) * | 2000-09-11 | 2000-10-05 | Nmr Holdings No. 2 Pty Limited | Asymmetric radio frequency coils for magnetic resonance |
| US7323876B2 (en) * | 2003-09-10 | 2008-01-29 | Koninklijke Philips Electronics N.V. | Magnetic resonance imaging receive chain with dynamic gain and wireless receiver coil |
| US8970217B1 (en) | 2010-04-14 | 2015-03-03 | Hypres, Inc. | System and method for noise reduction in magnetic resonance imaging |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6078337A (ja) * | 1983-07-21 | 1985-05-04 | ザ・リ−ジエンツ・オブ・ザ・ユニバ−シテイ・オブ・カリフオルニア | 矢状又は冠状ｎｍｒ画像形成におけるエイリアシングを低減する方法及びその装置 |
| JPS6117054A (ja) * | 1984-07-03 | 1986-01-25 | Yokogawa Medical Syst Ltd | 核磁気共鳴断層撮影装置 |
| JPS63108254A (ja) * | 1986-10-24 | 1988-05-13 | Jeol Ltd | ２次元核磁気共鳴における位相補正方法 |
| JPS63302841A (ja) * | 1987-06-03 | 1988-12-09 | Fuji Electric Co Ltd | 核磁気共鳴映像装置の計測デ−タ誤差補償方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB1496886A (en) * | 1974-04-08 | 1978-01-05 | Nat Res Dev | Nuclear magnetic resonance apparatus |
| 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. |
| US4703267A (en) * | 1985-07-17 | 1987-10-27 | Advanced Nmr Systems, Inc. | High dynamic range in NMR data acquisition |
| US4716368A (en) * | 1985-08-09 | 1987-12-29 | Picker International, Inc. | Magnetic resonance reconstruction and scanning techniques using known information, constraints, and symmetry relations |
| JPS62231647A (ja) * | 1986-03-31 | 1987-10-12 | 株式会社東芝 | 磁気共鳴映像装置 |
| US4703270A (en) * | 1986-04-18 | 1987-10-27 | The University Of British Columbia | Zero quantum NMR imaging and spectroscopy in a low homogeneity magnetic field |
-
1988
- 1988-11-30 US US07/278,121 patent/US4885549A/en not_active Expired - Fee Related
-
1989
- 1989-11-17 EP EP89121267A patent/EP0371349B1/en not_active Expired - Lifetime
- 1989-11-17 DE DE68925620T patent/DE68925620T2/de not_active Expired - Fee Related
- 1989-11-29 JP JP1307835A patent/JPH02220635A/ja active Granted
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6078337A (ja) * | 1983-07-21 | 1985-05-04 | ザ・リ−ジエンツ・オブ・ザ・ユニバ−シテイ・オブ・カリフオルニア | 矢状又は冠状ｎｍｒ画像形成におけるエイリアシングを低減する方法及びその装置 |
| JPS6117054A (ja) * | 1984-07-03 | 1986-01-25 | Yokogawa Medical Syst Ltd | 核磁気共鳴断層撮影装置 |
| JPS63108254A (ja) * | 1986-10-24 | 1988-05-13 | Jeol Ltd | ２次元核磁気共鳴における位相補正方法 |
| JPS63302841A (ja) * | 1987-06-03 | 1988-12-09 | Fuji Electric Co Ltd | 核磁気共鳴映像装置の計測デ−タ誤差補償方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| DE68925620D1 (de) | 1996-03-21 |
| DE68925620T2 (de) | 1996-09-26 |
| EP0371349A2 (en) | 1990-06-06 |
| EP0371349B1 (en) | 1996-02-07 |
| EP0371349A3 (en) | 1991-02-06 |
| US4885549A (en) | 1989-12-05 |
| JPH0422575B2 (ja) | 1992-04-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5289127A (en) | Correction of signal distortion in an NMR apparatus | |
| EP0577770B1 (en) | Method for correcting spectral and imaging data and for using such corrected data in magnet shimming | |
| US4876507A (en) | Method of and device for automatic phase correction of complex NMR spectra | |
| Liao et al. | Modal parameter identification using the log decrement method and band-pass filters | |
| US4438400A (en) | Method for suppression of acoustic ringing in NMR measurements | |
| JP4071625B2 (ja) | デカップリングサイドバンド分解ｎｍｒ分光におけるパルスシーケンス法 | |
| JPS5938636A (ja) | 核磁気共鳴装置 | |
| US6066949A (en) | Gradient characterization using fourier-transform | |
| JPH0113054B2 (ja) | ||
| US5311132A (en) | Method of enhancing the focus of magnetic resonance images | |
| US4857844A (en) | Method of and device for automatic phase correction of complex NMR spectra | |
| JPH02220635A (ja) | 基準マーカを使用したｎｍｒ信号位相および振幅補正方式 | |
| JPS6297544A (ja) | Ｎｍｒ信号受信方法 | |
| Liu et al. | Adaptive pre-whiten filtering for the free induction decay transversal signal in weak magnetic detection | |
| JPH04262281A (ja) | 勾配増強ｎｍｒ相関分光法 | |
| Weiss et al. | A study of precision in the measurement of chemical shifts | |
| JPH1090384A (ja) | Ｎｒｍにおける放射減衰の抑制 | |
| JP3478924B2 (ja) | 核磁気共鳴スペクトルの自動位相補正装置 | |
| US7166998B2 (en) | Exact half pulse synthesis via the inverse scattering transform | |
| US5221899A (en) | Signal acquisition in magnetic resonance analysis | |
| US4780673A (en) | Acquisition and processing of spin-echo NMR spectra | |
| US4677383A (en) | Two-dimensional nuclear magnetic resonance spectrometry | |
| US4766382A (en) | Two-dimensional nuclear magnetic resonance spectrometry | |
| AU655079B2 (en) | Compensation of pulsed linear magnet field gradient inducted eddy currents in NMR systems | |
| US6268729B1 (en) | Correction of magnetic resonance signals in a magnetic resonance imaging method, and a magnetic resonance device |