JP4498947B2 - 磁気共鳴スペクトルの定量方法 - Google Patents
磁気共鳴スペクトルの定量方法 Download PDFInfo
- Publication number
- JP4498947B2 JP4498947B2 JP2005046419A JP2005046419A JP4498947B2 JP 4498947 B2 JP4498947 B2 JP 4498947B2 JP 2005046419 A JP2005046419 A JP 2005046419A JP 2005046419 A JP2005046419 A JP 2005046419A JP 4498947 B2 JP4498947 B2 JP 4498947B2
- Authority
- JP
- Japan
- Prior art keywords
- spectrum
- magnetic resonance
- value
- peak
- resonance spectrum
- 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 - Fee Related
Links
- 238000001228 spectrum Methods 0.000 title claims description 141
- 238000000034 method Methods 0.000 title claims description 59
- 238000011002 quantification Methods 0.000 title claims description 15
- 230000005415 magnetization Effects 0.000 claims description 31
- 238000001514 detection method Methods 0.000 claims description 15
- 238000000655 nuclear magnetic resonance spectrum Methods 0.000 claims description 13
- 238000005259 measurement Methods 0.000 claims description 12
- 238000012417 linear regression Methods 0.000 claims description 9
- 230000010354 integration Effects 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000001362 electron spin resonance spectrum Methods 0.000 claims description 2
- 230000003595 spectral effect Effects 0.000 claims description 2
- 230000005389 magnetism Effects 0.000 claims 1
- 238000012937 correction Methods 0.000 description 11
- 238000005481 NMR spectroscopy Methods 0.000 description 10
- 238000002474 experimental method Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 238000001646 magnetic resonance method Methods 0.000 description 2
- 230000010363 phase shift Effects 0.000 description 2
- 238000004611 spectroscopical analysis Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000004435 EPR spectroscopy Methods 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 239000003480 eluent Substances 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000004316 pulse electron spin resonance spectroscopy Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000013214 routine measurement Methods 0.000 description 1
- 238000004448 titration Methods 0.000 description 1
Images
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/46—NMR spectroscopy
- G01R33/4625—Processing of acquired signals, e.g. elimination of phase errors, baseline fitting, chemometric analysis
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10T—TECHNICAL SUBJECTS COVERED BY FORMER US CLASSIFICATION
- Y10T436/00—Chemistry: analytical and immunological testing
- Y10T436/24—Nuclear magnetic resonance, electron spin resonance or other spin effects or mass spectrometry
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Description
本発明は、磁気共鳴装置で用いられる磁気共鳴スペクトルの定量方法に関する。
一般に、フーリエ変換核磁気共鳴（ＦＴ−ＮＭＲ）法やフーリエ変換電子スピン共鳴（ＦＴ−ＥＳＲ）法などのフーリエ変換磁気共鳴法においては、観測したスペクトルに実数部と虚数部が混合していて、スペクトル観測後に位相補正を行なって、実数部だけのスペクトルと虚数部だけのスペクトルに分離する作業が必要とされる。
これまで、このような磁気共鳴スペクトルを定量する際には、手動またはその他の自動位相補正法を用いて、スペクトルに対して適切な位相補正を行なった後に、スペクトルの実数部の所定の部分を切り出して、その範囲の積分値を求める方法が用いられてきた。
このような場合、積分前に、必ず位相補正を行なう必要があるが、大小複数のピークが含まれるスペクトルに対する自動位相補正は、非常に困難であり、手動での位相補正が必ず必要であった。
従来法では、定量解析手続き中に、手動での位相補正が必ず介在してくるために、解析処理の他の部分を自動化したとしても、完全な自動化はできないという問題があった。
本発明の目的は、上述した点に鑑み、磁気共鳴スペクトルの実数部と虚数部の理論的な関係に基づいて、各ピークの位相補正を行なうことなく、すべてのピークに対して定量が行なえるような、磁気共鳴スペクトルの定量方法を提供することにある。
この目的を達成するため、本発明にかかる磁気共鳴スペクトルの定量方法は、
磁気共鳴スペクトルにおけるスペクトルピークの積分値強度を求める方法であって、
（１）直角検波によって実数部スペクトルと虚数部スペクトルを求め、
（２）得られた実数部スペクトル及び虚数部スペクトルについて、積分対象領域に渡る積分値ａ、ｂをそれぞれ求め、
（３）前記積分値強度を±sqrt（ａ2＋ｂ2）として求めることを特徴としている。（ここで「sqrt」は平方根「square root」の略である。以下同様。）
また、前記請求項１記載の磁気共鳴スペクトルの定量方法により得られた磁気共鳴ピーク（磁化）の積分値強度を用いて、磁気共鳴スペクトル測定時の最適パルス幅の自動決定を行なうことを特徴としている。
磁気共鳴スペクトルにおけるスペクトルピークの積分値強度を求める方法であって、
（１）直角検波によって実数部スペクトルと虚数部スペクトルを求め、
（２）得られた実数部スペクトル及び虚数部スペクトルについて、積分対象領域に渡る積分値ａ、ｂをそれぞれ求め、
（３）前記積分値強度を±sqrt（ａ2＋ｂ2）として求めることを特徴としている。（ここで「sqrt」は平方根「square root」の略である。以下同様。）
また、前記請求項１記載の磁気共鳴スペクトルの定量方法により得られた磁気共鳴ピーク（磁化）の積分値強度を用いて、磁気共鳴スペクトル測定時の最適パルス幅の自動決定を行なうことを特徴としている。
また、磁気共鳴スペクトルにおけるスペクトルピークの積分値強度を求める方法であって、
（１）直角検波によって実数部と虚数部とから成る一次元磁気共鳴スペクトルを観測し収集する工程、
（２）得られたｎ個の一次元磁気共鳴スペクトルに対して、各スペクトルの所定の周波数範囲の実数部と虚数部のそれぞれの積分値を求める工程、
（３）（２）の工程で得られた実数部の積分値ｘiと虚数部の積分値ｙi（ただし、i＝１，２，……，ｎ）を複素平面上で直線ｙ＝ａｘ＋ｂに一次回帰させてフィッティングさせる（ただし、ｂはオフセット値）工程、
（４）実測データ点（ｘi，ｙi）をフィッティング直線ｙ＝ａｘ＋ｂへ投影し、投影した点を新規の（ｘi，ｙi）とする工程、
（５）ピーク強度がゼロとなる場合の投影データ点（ｘ0，ｙ0）から各投影データ点（ｘi，ｙi）までのデータ点間距離（sqrt((ｘi−ｘ0)2＋(ｙi−ｙ0)2)の値を求め、この値をスペクトルピーク（磁化）の積分値強度とする工程、
から成ることを特徴としている。
（１）直角検波によって実数部と虚数部とから成る一次元磁気共鳴スペクトルを観測し収集する工程、
（２）得られたｎ個の一次元磁気共鳴スペクトルに対して、各スペクトルの所定の周波数範囲の実数部と虚数部のそれぞれの積分値を求める工程、
（３）（２）の工程で得られた実数部の積分値ｘiと虚数部の積分値ｙi（ただし、i＝１，２，……，ｎ）を複素平面上で直線ｙ＝ａｘ＋ｂに一次回帰させてフィッティングさせる（ただし、ｂはオフセット値）工程、
（４）実測データ点（ｘi，ｙi）をフィッティング直線ｙ＝ａｘ＋ｂへ投影し、投影した点を新規の（ｘi，ｙi）とする工程、
（５）ピーク強度がゼロとなる場合の投影データ点（ｘ0，ｙ0）から各投影データ点（ｘi，ｙi）までのデータ点間距離（sqrt((ｘi−ｘ0)2＋(ｙi−ｙ0)2)の値を求め、この値をスペクトルピーク（磁化）の積分値強度とする工程、
から成ることを特徴としている。
また、磁気共鳴スペクトルにおけるスペクトルピークの積分値強度を求める方法であって、
（１）直角検波によって実数部と虚数部とから成る一次元磁気共鳴スペクトルを観測し収集する工程、
（２）得られたｎ個の一次元磁気共鳴スペクトルに対して、各スペクトルの所定の周波数範囲の実数部と虚数部のそれぞれの積分値を求める工程、
（３）（２）の工程で得られた実数部の積分値ｘiと虚数部の積分値ｙi（ただし、i＝１，２，……，ｎ）を複素平面上で直線ｙ＝ａｘ＋ｂに一次回帰させてフィッティングさせる（ただし、ｂはオフセット値）工程、
（４）実測データ点（ｘi，ｙi）をフィッティング直線ｙ＝ａｘ＋ｂへ投影し、投影した点を新規の（ｘi，ｙi）とする工程、
（５）ピーク強度がゼロとなる場合の投影データ点（ｘ0，ｙ0）から各投影データ点（ｘi，ｙi）までのデータ点間距離（sqrt((ｘi−ｘ0)2＋(ｙi−ｙ0)2)の値を求め、この値をスペクトルピーク（磁化）の積分値強度の絶対値とする工程、
（６）基準データを１点目の座標値（ｘ1，ｙ1）とし、ｉ番目のデータ（ｘi，ｙi）に対して、(ｘi−ｘ0)／(ｘ1−ｘ0)の符号（signi）を求め、この符号をスペクトルピーク（磁化）の符号とする工程、
（７）（５）の工程で得られた絶対値と（６）の工程で得られた符号とに基づいて、（signi×sqrt((ｘi−ｘ0)2＋(ｙi−ｙ0)2)の値を求め、この値をスペクトルピーク（磁化）の積分値強度とする工程、
から成ることを特徴としている。
（１）直角検波によって実数部と虚数部とから成る一次元磁気共鳴スペクトルを観測し収集する工程、
（２）得られたｎ個の一次元磁気共鳴スペクトルに対して、各スペクトルの所定の周波数範囲の実数部と虚数部のそれぞれの積分値を求める工程、
（３）（２）の工程で得られた実数部の積分値ｘiと虚数部の積分値ｙi（ただし、i＝１，２，……，ｎ）を複素平面上で直線ｙ＝ａｘ＋ｂに一次回帰させてフィッティングさせる（ただし、ｂはオフセット値）工程、
（４）実測データ点（ｘi，ｙi）をフィッティング直線ｙ＝ａｘ＋ｂへ投影し、投影した点を新規の（ｘi，ｙi）とする工程、
（５）ピーク強度がゼロとなる場合の投影データ点（ｘ0，ｙ0）から各投影データ点（ｘi，ｙi）までのデータ点間距離（sqrt((ｘi−ｘ0)2＋(ｙi−ｙ0)2)の値を求め、この値をスペクトルピーク（磁化）の積分値強度の絶対値とする工程、
（６）基準データを１点目の座標値（ｘ1，ｙ1）とし、ｉ番目のデータ（ｘi，ｙi）に対して、(ｘi−ｘ0)／(ｘ1−ｘ0)の符号（signi）を求め、この符号をスペクトルピーク（磁化）の符号とする工程、
（７）（５）の工程で得られた絶対値と（６）の工程で得られた符号とに基づいて、（signi×sqrt((ｘi−ｘ0)2＋(ｙi−ｙ0)2)の値を求め、この値をスペクトルピーク（磁化）の積分値強度とする工程、
から成ることを特徴としている。
また、前記請求項４記載の磁気共鳴スペクトルの定量方法により得られたスペクトルピーク（磁化）の積分値強度を用いて、磁気共鳴スペクトル測定時の最適パルス幅の自動決定を行なうことを特徴としている。
また、磁気共鳴スペクトルにおけるスペクトルピークの積分値強度を求める方法であって、
（１）直角検波によって実数部と虚数部とから成る一次元磁気共鳴スペクトルを観測し収集する工程、
（２）得られたｎ個の一次元磁気共鳴スペクトルに対して、各スペクトルの所定の周波数範囲の実数部と虚数部のそれぞれの積分値を求める工程、
（３）（２）の工程で得られた実数部の積分値ｘiと虚数部の積分値ｙi（ただし、i＝１，２，……，ｎ）を複素平面上で直線ｙ＝ａｘ＋ｂに一次回帰させてフィッティングさせる（ただし、ｂはオフセット値）工程、
（４）実測データ点（ｘi，ｙi）をフィッティング直線ｙ＝ａｘ＋ｂへ投影し、投影した点を新規の（ｘi，ｙi）とする工程、
（５）基準データを１点目の座標値（ｘ1，ｙ1）とし、ｉ番目のデータ（ｘi，ｙi）に対して、(ｘi−ｘ1)の符号（signi）を求め、この符号をスペクトルピーク（磁化）の符号とする工程、
（６）（５）の工程で得られた符号（signi）に基づいて、（signi×sqrt((ｘi−ｘ１)2＋(ｙi−ｙ1)2)の値を求め、その値をパルス幅計算用最小自乗フィッティング・プログラムに渡す工程、
（７）パルス幅計算用最小自乗フィッティング・プログラムに基づいて、スペクトルピーク（磁化）の大きさのパルス幅依存性を、所定のモデル関数にカーブ・フィッティングさせ、フィッティングされたモデル関数の諸パラメータから、スペクトルの観測に最適なパルス幅を決定する工程、
から成ることを特徴としている。
（１）直角検波によって実数部と虚数部とから成る一次元磁気共鳴スペクトルを観測し収集する工程、
（２）得られたｎ個の一次元磁気共鳴スペクトルに対して、各スペクトルの所定の周波数範囲の実数部と虚数部のそれぞれの積分値を求める工程、
（３）（２）の工程で得られた実数部の積分値ｘiと虚数部の積分値ｙi（ただし、i＝１，２，……，ｎ）を複素平面上で直線ｙ＝ａｘ＋ｂに一次回帰させてフィッティングさせる（ただし、ｂはオフセット値）工程、
（４）実測データ点（ｘi，ｙi）をフィッティング直線ｙ＝ａｘ＋ｂへ投影し、投影した点を新規の（ｘi，ｙi）とする工程、
（５）基準データを１点目の座標値（ｘ1，ｙ1）とし、ｉ番目のデータ（ｘi，ｙi）に対して、(ｘi−ｘ1)の符号（signi）を求め、この符号をスペクトルピーク（磁化）の符号とする工程、
（６）（５）の工程で得られた符号（signi）に基づいて、（signi×sqrt((ｘi−ｘ１)2＋(ｙi−ｙ1)2)の値を求め、その値をパルス幅計算用最小自乗フィッティング・プログラムに渡す工程、
（７）パルス幅計算用最小自乗フィッティング・プログラムに基づいて、スペクトルピーク（磁化）の大きさのパルス幅依存性を、所定のモデル関数にカーブ・フィッティングさせ、フィッティングされたモデル関数の諸パラメータから、スペクトルの観測に最適なパルス幅を決定する工程、
から成ることを特徴としている。
また、前記磁気共鳴スペクトルは、ＮＭＲスペクトルであることを特徴としている。
また、前記磁気共鳴スペクトルは、ＥＳＲスペクトルであることを特徴としている。
本発明の磁気共鳴スペクトルの定量方法によれば、磁気共鳴スペクトルにおけるスペクトルピークの積分値強度を求める方法であって、
（１）直角検波によって実数部スペクトルと虚数部スペクトルを求め、
（２）得られた実数部スペクトル及び虚数部スペクトルについて、積分対象領域に渡る積分値ａ、ｂをそれぞれ求め、
（３）前記積分値強度を±sqrt（ａ2＋ｂ2）として求めるようにしたので、
磁気共鳴スペクトルの実数部と虚数部の理論的な関係に基づいて、各ピークの位相補正を行なうことなく、すべてのピークに対して定量が行なえるようになった。
（１）直角検波によって実数部スペクトルと虚数部スペクトルを求め、
（２）得られた実数部スペクトル及び虚数部スペクトルについて、積分対象領域に渡る積分値ａ、ｂをそれぞれ求め、
（３）前記積分値強度を±sqrt（ａ2＋ｂ2）として求めるようにしたので、
磁気共鳴スペクトルの実数部と虚数部の理論的な関係に基づいて、各ピークの位相補正を行なうことなく、すべてのピークに対して定量が行なえるようになった。
一般に、分光法では、エネルギーに比例した量を横軸とし、それに対する吸収や共鳴の強度を縦軸として、スペクトルを得ることが多い。ＦＴ−ＮＭＲ法やＦＴ−ＥＳＲ法などのフーリエ変換磁気共鳴分光法の場合は、ある周波数に対する共鳴の強度をスペクトルとして得るが、このとき、横軸の共鳴周波数は、共鳴エネルギーに比例する量であることが知られている。したがって、共鳴エネルギーの分布は、そのまま共鳴周波数の分布になるため、その実数部の分布は、一般に偶関数で表わされる分布を取る。
ここで、偶関数の例として、ローレンツ型分布（コーシー分布）関数の場合を例に取ると、共鳴ピークの実数部の形状ｆ(ｘ)は、次のような式で表わされる。
ここで、πは円周率、Ｔ2は緩和時間に関する定数、ｘ0は共鳴ピークの中心周波数である。このとき、互いに９０゜位相の異なる、実数部のスペクトル（ｆ(ｘ)）と虚数部のスペクトル（ｇ(ｘ)）は、ヒルベルト変換の関係にあるので、共鳴ピークの虚数部の形状ｇ(ｘ)は、次のような式で表わされる。
したがって、虚数部スペクトルと実数部スペクトルとは、積分と微分の関係にあり、実数部スペクトルが偶関数のとき、虚数部スペクトルは奇関数になる。
この関係は、スペクトルのピーク形状がローレンツ型分布でないときにも、一般に得られたピークを偶関数の重ね合わせとして表わすことが可能なので、実数部と虚数部が、偶関数と奇関数である、という関係は、すべてのスペクトルに対して成り立つ関係である。
更に、このスペクトルを積分した値は、
となる。つまり、観測されたスペクトルの位相が、適切に補正されている場合には、虚数部スペクトルの積分値は、ゼロになる。また、実数部スペクトルの積分値は、実際には規格化されていないので、ピーク強度に比例した値を持つことになる。
次に、実際に直角検波により観測された、互いに９０゜位相の異なるスペクトル（ｆ'(ｘ)，g'(x)）の位相が、それぞれｐ0だけずれていた場合を考えると、
となる。このとき、直角検波により観測された、互いに９０゜位相の異なるスペクトル（ｆ'(ｘ)，g'(x)）の積分値は、
となる。
このとき、上式中の積分項は、観測対象の磁化に比例する量であるから、観測対象の磁化が同一種のときは、強度の異なる同一種の磁化に対する実数部と虚数部を積分した値を、複素平面上にプロットすれば、磁化の大きさを定量することができる。また、このとき、プロットされた点は、複素平面上では、同一直線上に配列するが、その直線の傾きは、位相ｐ0の正弦・余弦成分を表わしているので、０次の位相補正値も同時に決定することができる。
つまり、図１で示すように、実数部の積分値ａを横軸、虚数部の積分値ｂを縦軸としてプロットしたときに、座標（ａ，ｂ）で表わされる点の原点からの距離sqrt（ａ2＋ｂ2）が、対象としているピークの磁化の強度に比例する量であり、この値は、従来の方法で、手動による位相補正を行なって、実数部の積分値を取った場合の値と、同じ値を取ることになる。すなわち、図１の中の複素平面上に引かれた直線上の０％、２５％、５０％、７５％、１００％に当たる点が、それぞれ観測している磁化の０％、２５％、５０％、７５％、１００％の大きさを表わしており、角度ｐ0は、その磁化の位相のずれを表わしている。
この解析方法において、手動で行なう必要があるのは、定量対象となるピークが含まれる領域を指定するだけであり、後はすべて自動で行なわれるため、フーリエ変換磁気共鳴法でのパルス幅変調実験や滴定実験などのピーク強度の変調を伴う実験結果の自動解析にとって、きわめて有効である。
ここで、実際に観測されるスペクトルにおいては、装置的にオフセットのないスペクトルを観測することが困難な場合がある。その場合、図１のような、信号が検出されない条件の下で測定されたスペクトル（ピーク強度が０％の場合のスペクトル）を実際に取得した上で、データを解析する必要がある。具体的な実施例を３つ示す。
まず、最初は、すでに分かっている成分の含有量をルーチン測定で定量する場合、すなわち、試料管を交換しながら同一物質を定量したり、液体クロマトグラフ装置からの溶離液をフローさせながら同一物質を定量したりする場合である。この場合、観測対象は既知物質であるから、ＮＭＲ装置のパルス幅を固定した測定になり、ピークの絶対値のみを観測することになる。スペクトルピークは、常に正の側にのみ表われる。図２は、そのような場合に本発明を適用したときのアルゴリズムを示したものである。測定の工程を以下に示す。
（１）濃度がゼロの場合を含む濃度の異なる複数の試料に対して、直角検波によって実数部と虚数部とから成る一次元ＮＭＲスペクトルを観測し収集する。（Ｓ−１１）
（２）得られたｎ個の一次元ＮＭＲスペクトルに対して、各スペクトルの所定の周波数範囲の実数部と虚数部のそれぞれの積分値を求める。（Ｓ−１２）
（３）（２）の工程で得られた実数部の積分値ｘiと虚数部の積分値ｙi（ただし、i＝１，２，……，ｎ）を複素平面上で直線ｙ＝ａｘ＋ｂに一次回帰させてフィッティングさせる（ただし、ｂはオフセット値）。（Ｓ−１３）
（４）実測データ点（ｘi，ｙi）をフィッティング直線ｙ＝ａｘ＋ｂへ投影し、投影した点を新規の（ｘi，ｙi）とする。（Ｓ−１４）
（５）濃度がゼロの場合の投影データ点（ｘ0，ｙ0）から各投影データ点（ｘi，ｙi）までのデータ点間距離（sqrt((ｘi−ｘ0)2＋(ｙi−ｙ0)2)の値を求め、この値をＮＭＲピーク（磁化）の積分値強度とする。（Ｓ−１５）
（１）濃度がゼロの場合を含む濃度の異なる複数の試料に対して、直角検波によって実数部と虚数部とから成る一次元ＮＭＲスペクトルを観測し収集する。（Ｓ−１１）
（２）得られたｎ個の一次元ＮＭＲスペクトルに対して、各スペクトルの所定の周波数範囲の実数部と虚数部のそれぞれの積分値を求める。（Ｓ−１２）
（３）（２）の工程で得られた実数部の積分値ｘiと虚数部の積分値ｙi（ただし、i＝１，２，……，ｎ）を複素平面上で直線ｙ＝ａｘ＋ｂに一次回帰させてフィッティングさせる（ただし、ｂはオフセット値）。（Ｓ−１３）
（４）実測データ点（ｘi，ｙi）をフィッティング直線ｙ＝ａｘ＋ｂへ投影し、投影した点を新規の（ｘi，ｙi）とする。（Ｓ−１４）
（５）濃度がゼロの場合の投影データ点（ｘ0，ｙ0）から各投影データ点（ｘi，ｙi）までのデータ点間距離（sqrt((ｘi−ｘ0)2＋(ｙi−ｙ0)2)の値を求め、この値をＮＭＲピーク（磁化）の積分値強度とする。（Ｓ−１５）
次に、同一試料に対して、ＮＭＲ側の測定条件を変えながら測定する場合、すなわち、パルス幅を変化させて、観測されるピーク強度の変化を観測するような場合である。この場合、変化する対象はパルス幅であるから、ピーク値は、正側に出る場合と負側に出る場合とがあり得る。したがって、この場合は、ピークの絶対値とピーク値の符号との両方の決定が必要な場合である。図３は、そのような場合に本発明を適用したときのアルゴリズムを示したものである。測定の工程を以下に示す。
（１）試料に対して、パルス幅を変えながら、直角検波によって実数部と虚数部とから成る一次元ＮＭＲスペクトルを観測し収集する。（Ｓ−２１）
（２）得られたｎ個の一次元ＮＭＲスペクトルに対して、各スペクトルの所定の周波数範囲の実数部と虚数部のそれぞれの積分値を求める。（Ｓ−２２）
（３）（２）の工程で得られた実数部の積分値ｘiと虚数部の積分値ｙi（ただし、i＝１，２，……，ｎ）を複素平面上で直線ｙ＝ａｘ＋ｂに一次回帰させてフィッティングさせる（ただし、ｂはオフセット値）。（Ｓ−２３）
（４）実測データ点（ｘi，ｙi）をフィッティング直線ｙ＝ａｘ＋ｂへ投影し、投影した点を新規の（ｘi，ｙi）とする。（Ｓ−２４）
（５）ＮＭＲ信号が検出されないようなパルス幅のパルスを試料に対して与えた場合（たとえば、パルス幅をゼロに設定した場合や、同じ幅のパルスを正負連続して与えた場合など）の投影データ点（ｘ0，ｙ0）から各投影データ点（ｘi，ｙi）までのデータ点間距離（sqrt((ｘi−ｘ0)2＋(ｙi−ｙ0)2)の値を求め、この値をＮＭＲピーク（磁化）の積分値強度の絶対値とする。（Ｓ−２５）
（６）基準データを１点目の座標値（ｘ1，ｙ1）とし、ｉ番目のデータ（ｘi，ｙi）に対して、(ｘi−ｘ0)／(ｘ1−ｘ0)の符号（signi）を求め、この符号をＮＭＲピーク（磁化）の符号とする。（Ｓ−２６）
（７）（５）の工程で得られた絶対値と（６）の工程で得られた符号とに基づいて、（signi×sqrt((ｘi−ｘ0)2＋(ｙi−ｙ0)2)の値を求め、この値をＮＭＲピーク（磁化）の積分値強度とする。（Ｓ−２７）
図２のアルゴリズムと図３のアルゴリズムとの違いは、投影データ点（ｘ0，ｙ0）から各投影データ点（ｘi，ｙi）までの距離を求めた後に、その距離の符号（正負）に関する情報を付加するか否かという点のみに過ぎない。したがって、図２のアルゴリズムと図３のアルゴリズムは、「ピークはすべて正のピークから成るか？」という判断工程を加えることで、図４のような１つのアルゴリズムで表わすこともできる。
（１）試料に対して、パルス幅を変えながら、直角検波によって実数部と虚数部とから成る一次元ＮＭＲスペクトルを観測し収集する。（Ｓ−２１）
（２）得られたｎ個の一次元ＮＭＲスペクトルに対して、各スペクトルの所定の周波数範囲の実数部と虚数部のそれぞれの積分値を求める。（Ｓ−２２）
（３）（２）の工程で得られた実数部の積分値ｘiと虚数部の積分値ｙi（ただし、i＝１，２，……，ｎ）を複素平面上で直線ｙ＝ａｘ＋ｂに一次回帰させてフィッティングさせる（ただし、ｂはオフセット値）。（Ｓ−２３）
（４）実測データ点（ｘi，ｙi）をフィッティング直線ｙ＝ａｘ＋ｂへ投影し、投影した点を新規の（ｘi，ｙi）とする。（Ｓ−２４）
（５）ＮＭＲ信号が検出されないようなパルス幅のパルスを試料に対して与えた場合（たとえば、パルス幅をゼロに設定した場合や、同じ幅のパルスを正負連続して与えた場合など）の投影データ点（ｘ0，ｙ0）から各投影データ点（ｘi，ｙi）までのデータ点間距離（sqrt((ｘi−ｘ0)2＋(ｙi−ｙ0)2)の値を求め、この値をＮＭＲピーク（磁化）の積分値強度の絶対値とする。（Ｓ−２５）
（６）基準データを１点目の座標値（ｘ1，ｙ1）とし、ｉ番目のデータ（ｘi，ｙi）に対して、(ｘi−ｘ0)／(ｘ1−ｘ0)の符号（signi）を求め、この符号をＮＭＲピーク（磁化）の符号とする。（Ｓ−２６）
（７）（５）の工程で得られた絶対値と（６）の工程で得られた符号とに基づいて、（signi×sqrt((ｘi−ｘ0)2＋(ｙi−ｙ0)2)の値を求め、この値をＮＭＲピーク（磁化）の積分値強度とする。（Ｓ−２７）
図２のアルゴリズムと図３のアルゴリズムとの違いは、投影データ点（ｘ0，ｙ0）から各投影データ点（ｘi，ｙi）までの距離を求めた後に、その距離の符号（正負）に関する情報を付加するか否かという点のみに過ぎない。したがって、図２のアルゴリズムと図３のアルゴリズムは、「ピークはすべて正のピークから成るか？」という判断工程を加えることで、図４のような１つのアルゴリズムで表わすこともできる。
カーブ・フィッティングを行なうことで、ピーク強度がゼロとなる点の実数部と虚数部の積分値（ｘ0，ｙ0）を解析的に決定できる（オフセットの補正機能がある）プログラムを用いる場合（例えば、特許文献４に記載のプログラムを用いる場合）には、図５に示すアルゴリズムが適用され、ピーク強度がゼロとなる点での実数部の積分値ｘ0と虚数部の積分値ｙ0を実測で取得する必要はない。実際に、この手法をパルス幅較正の際のニューテーション実験に適用する際には、以下のような工程を行なう。
（１）試料に対して、パルス幅を変えながら、直角検波によって実数部と虚数部とから成る一次元ＮＭＲスペクトルを観測し収集する。（Ｓ−３１）
（２）得られたｎ個の一次元ＮＭＲスペクトルに対して、各スペクトルの所定の周波数範囲の実数部と虚数部のそれぞれの積分値を求める。（Ｓ−３２）
（３）（２）の工程で得られた実数部の積分値ｘiと虚数部の積分値ｙi（ただし、i＝１，２，……，ｎ）を複素平面上で直線ｙ＝ａｘ＋ｂに一次回帰させてフィッティングさせる（ただし、ｂはオフセット値）。（Ｓ−３３）
（４）実測データ点（ｘi，ｙi）をフィッティング直線ｙ＝ａｘ＋ｂへ投影し、投影した点を新規の（ｘi，ｙi）とする。（Ｓ−３４）
（５）基準データを１点目の座標値（ｘ1，ｙ1）とし、ｉ番目のデータ（ｘi，ｙi）に対して、(ｘi−ｘ1）の符号（signi）を求め、この符号をＮＭＲピーク（磁化）の符号とする。（Ｓ−３５）
（６）（５）の工程で得られた符号（signi）に基づいて、（signi×sqrt((ｘi−ｘ１)2＋(ｙi−ｙ1)2)の値を求め、その値を最適パルス幅の自動決定技術（特許文献４）に記載のパルス幅計算用最小自乗フィッティング・プログラムに渡す。（Ｓ−３６）
（７）特許文献４に記載の手法に基づいて、磁化の大きさのパルス幅依存性を、所定のモデル関数にカーブ・フィッティングさせ、フィッティングされたモデル関数の諸パラメータから、磁化のオフセット値（ｘ0，ｙ0）を割り出すとともに、スペクトルの観測に最適なパルス幅（たとえば９０°パルスの幅）を自動的に決定する。（Ｓ−３７）
実際に、パルス幅較正のためのニューテーション実験のデータに適用した例を示す。図６は、シングルパルスによる１１本のスペクトルを並べた図で、図６（ａ）が実数部のスペクトル、図６（ｂ）が虚数部のスペクトルである。このデータに対して、前記最適パルス幅の自動決定技術（特許文献４）に記載の手法で最適なパルス幅を決定するには、実数部スペクトルの積分値を観測時の励起パルス幅に対してプロットし、その値を、所定のモデル関数に、最小自乗法でカーブ・フィッティングさせることが必要である。スペクトルの観測に最適なパルス幅は、フィッティングされたときのモデル関数の諸パラメータから、自動的に決定することができる。
（１）試料に対して、パルス幅を変えながら、直角検波によって実数部と虚数部とから成る一次元ＮＭＲスペクトルを観測し収集する。（Ｓ−３１）
（２）得られたｎ個の一次元ＮＭＲスペクトルに対して、各スペクトルの所定の周波数範囲の実数部と虚数部のそれぞれの積分値を求める。（Ｓ−３２）
（３）（２）の工程で得られた実数部の積分値ｘiと虚数部の積分値ｙi（ただし、i＝１，２，……，ｎ）を複素平面上で直線ｙ＝ａｘ＋ｂに一次回帰させてフィッティングさせる（ただし、ｂはオフセット値）。（Ｓ−３３）
（４）実測データ点（ｘi，ｙi）をフィッティング直線ｙ＝ａｘ＋ｂへ投影し、投影した点を新規の（ｘi，ｙi）とする。（Ｓ−３４）
（５）基準データを１点目の座標値（ｘ1，ｙ1）とし、ｉ番目のデータ（ｘi，ｙi）に対して、(ｘi−ｘ1）の符号（signi）を求め、この符号をＮＭＲピーク（磁化）の符号とする。（Ｓ−３５）
（６）（５）の工程で得られた符号（signi）に基づいて、（signi×sqrt((ｘi−ｘ１)2＋(ｙi−ｙ1)2)の値を求め、その値を最適パルス幅の自動決定技術（特許文献４）に記載のパルス幅計算用最小自乗フィッティング・プログラムに渡す。（Ｓ−３６）
（７）特許文献４に記載の手法に基づいて、磁化の大きさのパルス幅依存性を、所定のモデル関数にカーブ・フィッティングさせ、フィッティングされたモデル関数の諸パラメータから、磁化のオフセット値（ｘ0，ｙ0）を割り出すとともに、スペクトルの観測に最適なパルス幅（たとえば９０°パルスの幅）を自動的に決定する。（Ｓ−３７）
実際に、パルス幅較正のためのニューテーション実験のデータに適用した例を示す。図６は、シングルパルスによる１１本のスペクトルを並べた図で、図６（ａ）が実数部のスペクトル、図６（ｂ）が虚数部のスペクトルである。このデータに対して、前記最適パルス幅の自動決定技術（特許文献４）に記載の手法で最適なパルス幅を決定するには、実数部スペクトルの積分値を観測時の励起パルス幅に対してプロットし、その値を、所定のモデル関数に、最小自乗法でカーブ・フィッティングさせることが必要である。スペクトルの観測に最適なパルス幅は、フィッティングされたときのモデル関数の諸パラメータから、自動的に決定することができる。
このような前記最適パルス幅の自動決定技術（特許文献４）に記載の手法で、このデータの最適パルス幅を決定した場合、図７（ａ）に示すように、９０゜パルス幅は、６．５１２３５μｓと決定することができる。
しかし、実際に実験で得られる磁気共鳴スペクトルの位相は、測定条件により異なるため、位相補正を行なわない自動処理の場合は、図６（ｂ）のような９０゜位相のずれたスペクトルに対して、処理を行なってしまう場合がある。実際に、図６（ａ）のスペクトルに対して行なった処理と同じ処理を、図６（ｂ）のスペクトルに対して行なうと、フィッティングした結果は、図７（ｂ）のようになり、用いたモデル関数にそぐわない、明らかに望まれないものとなってしまう。
ここで、このときの実数部と虚数部の積分値を、複素平面上にプロットしてみると、図８（ａ）のようになる。この図から、それぞれの積分値を表わした点は、ほぼ一直線上に並んでいることが分かる。したがって、このようなスペクトルに対しては、図８（ｂ）に示すように、すべてのデータ点が最も良くフィットするような所定の直線への回帰を行なわせ、更に、図８（ｃ）に示すように、回帰した直線上へのデータ点の投影（データ点の変換）を行なわせることにより、位相が適切に補正されている場合の実数部の積分値に対応する値を得ることができる。
なお、ここで、データの直線への投影を行なうとは、回帰した直線と、その直線に対してデータ点から引かれた垂線との交点の座標に、データ点の座標を変換することを言う。例えば、回帰した直線がｙ＝ａｘ＋ｂで表わされる直線のとき、ｉ番目のデータ点（ｘi，ｙi）は、
ｙ＝ａｘ＋ｂ
と、その垂線、
ｙ＝−１／ａ（ｘ−ｘi）＋ｙi
との交点の座標（(ａｙi＋ｘi−ａｂ)／(ａ2＋１)，(ａ2ｙi＋ｘi＋ｂ)／(ａ2＋１)）に変換される。
ｙ＝ａｘ＋ｂ
と、その垂線、
ｙ＝−１／ａ（ｘ−ｘi）＋ｙi
との交点の座標（(ａｙi＋ｘi−ａｂ)／(ａ2＋１)，(ａ2ｙi＋ｘi＋ｂ)／(ａ2＋１)）に変換される。
このようにして得られた値に対して、前記最適パルス幅の自動決定技術（特許文献４）に記載の手法により、図６（ａ）と図６（ｂ）に対して行なった処理と同じ処理を行なったところ、図９に示すように、最適な９０゜パルス幅の値は、６．５０７３４μｓ〜６．５７４４３μｓという結果が得られた。この値は、あらかじめ、手動でスペクトルの位相補正を行なうことによって得た値と、ほぼ一致する値である。したがって、本発明にかかる磁気共鳴スペクトルの定量方法によって定量された磁化の大きさを、特許文献４に記載のスペクトルの観測に最適なパルス幅の決定方法に当てはめれば、位相ずれの極端に大きな実数部スペクトルを用いても、スペクトルの観測に最適なパルス幅を、正確に決定することができるようになる。
尚、この方法は、パルスＮＭＲによる測定だけでなく、パルスＥＳＲによる測定にも適用することができる。
フーリエ変換磁気共鳴装置に、広く利用できる。
Claims (8)
- 磁気共鳴スペクトルにおけるスペクトルピークの積分値強度を求める方法であって、
（１）直角検波によって実数部スペクトルと虚数部スペクトルを求め、
（２）得られた実数部スペクトル及び虚数部スペクトルについて、積分対象領域に渡る積分値ａ、ｂをそれぞれ求め、
（３）前記積分値強度を±sqrt（ａ2＋ｂ2）として求めることを特徴とする磁気共鳴スペクトルの定量方法。（ここで「sqrt」は平方根「square root」の略である。以下同様。） - 前記請求項１記載の磁気共鳴スペクトルの定量方法により得られた磁気共鳴ピーク（磁化）の積分値強度を用いて、磁気共鳴スペクトル測定時の最適パルス幅の自動決定を行なうことを特徴とする磁気共鳴スペクトルの定量方法。
- 磁気共鳴スペクトルにおけるスペクトルピークの積分値強度を求める方法であって、
（１）直角検波によって実数部と虚数部とから成る一次元磁気共鳴スペクトルを観測し収集する工程、
（２）得られたｎ個の一次元磁気共鳴スペクトルに対して、各スペクトルの所定の周波数範囲の実数部と虚数部のそれぞれの積分値を求める工程、
（３）（２）の工程で得られた実数部の積分値ｘiと虚数部の積分値ｙi（ただし、i＝１，２，……，ｎ）を複素平面上で直線ｙ＝ａｘ＋ｂに一次回帰させてフィッティングさせる（ただし、ｂはオフセット値）工程、
（４）実測データ点（ｘi，ｙi）をフィッティング直線ｙ＝ａｘ＋ｂへ投影し、投影した点を新規の（ｘi，ｙi）とする工程、
（５）ピーク強度がゼロとなる場合の投影データ点（ｘ0，ｙ0）から各投影データ点（ｘi，ｙi）までのデータ点間距離（sqrt((ｘi−ｘ0)2＋(ｙi−ｙ0)2)の値を求め、この値をスペクトルピーク（磁化）の積分値強度とする工程、
から成ることを特徴とする磁気共鳴スペクトルの定量方法。 - 磁気共鳴スペクトルにおけるスペクトルピークの積分値強度を求める方法であって、
（１）直角検波によって実数部と虚数部とから成る一次元磁気共鳴スペクトルを観測し収集する工程、
（２）得られたｎ個の一次元磁気共鳴スペクトルに対して、各スペクトルの所定の周波数範囲の実数部と虚数部のそれぞれの積分値を求める工程、
（３）（２）の工程で得られた実数部の積分値ｘiと虚数部の積分値ｙi（ただし、i＝１，２，……，ｎ）を複素平面上で直線ｙ＝ａｘ＋ｂに一次回帰させてフィッティングさせる（ただし、ｂはオフセット値）工程、
（４）実測データ点（ｘi，ｙi）をフィッティング直線ｙ＝ａｘ＋ｂへ投影し、投影した点を新規の（ｘi，ｙi）とする工程、
（５）ピーク強度がゼロとなる場合の投影データ点（ｘ0，ｙ0）から各投影データ点（ｘi，ｙi）までのデータ点間距離（sqrt((ｘi−ｘ0)2＋(ｙi−ｙ0)2)の値を求め、この値をスペクトルピーク（磁化）の積分値強度の絶対値とする工程、
（６）基準データを１点目の座標値（ｘ1，ｙ1）とし、ｉ番目のデータ（ｘi，ｙi）に対して、(ｘi−ｘ0)／(ｘ1−ｘ0)の符号（signi）を求め、この符号をスペクトルピーク（磁化）の符号とする工程、
（７）（５）の工程で得られた絶対値と（６）の工程で得られた符号とに基づいて、（signi×sqrt((ｘi−ｘ0)2＋(ｙi−ｙ0)2)の値を求め、この値をスペクトルピーク（磁化）の積分値強度とする工程、
から成ることを特徴とする磁気共鳴スペクトルの定量方法。 - 前記請求項４記載の磁気共鳴スペクトルの定量方法により得られたスペクトルピーク（磁化）の積分値強度を用いて、磁気共鳴スペクトル測定時の最適パルス幅の自動決定を行なうことを特徴とする磁気共鳴スペクトルの定量方法。
- 磁気共鳴スペクトルにおけるスペクトルピークの積分値強度を求める方法であって、
（１）直角検波によって実数部と虚数部とから成る一次元磁気共鳴スペクトルを観測し収集する工程、
（２）得られたｎ個の一次元磁気共鳴スペクトルに対して、各スペクトルの所定の周波数範囲の実数部と虚数部のそれぞれの積分値を求める工程、
（３）（２）の工程で得られた実数部の積分値ｘiと虚数部の積分値ｙi（ただし、i＝１，２，……，ｎ）を複素平面上で直線ｙ＝ａｘ＋ｂに一次回帰させてフィッティングさせる（ただし、ｂはオフセット値）工程、
（４）実測データ点（ｘi，ｙi）をフィッティング直線ｙ＝ａｘ＋ｂへ投影し、投影した点を新規の（ｘi，ｙi）とする工程、
（５）基準データを１点目の座標値（ｘ1，ｙ1）とし、ｉ番目のデータ（ｘi，ｙi）に対して、(ｘi−ｘ1)の符号（signi）を求め、この符号をスペクトルピーク（磁化）の符号とする工程、
（６）（５）の工程で得られた符号（signi）に基づいて、（signi×sqrt((ｘi−ｘ１)2＋(ｙi−ｙ1)2)の値を求め、その値をパルス幅計算用最小自乗フィッティング・プログラムに渡す工程、
（７）パルス幅計算用最小自乗フィッティング・プログラムに基づいて、スペクトルピーク（磁化）の大きさのパルス幅依存性を、所定のモデル関数にカーブ・フィッティングさせ、フィッティングされたモデル関数の諸パラメータから、スペクトルの観測に最適なパルス幅を決定する工程、
から成ることを特徴とする磁気共鳴スペクトルの定量方法。 - 前記磁気共鳴スペクトルは、ＮＭＲスペクトルであることを特徴とする請求項１ないし６のいずれか１つに記載の磁気共鳴スペクトルの定量方法。
- 前記磁気共鳴スペクトルは、ＥＳＲスペクトルであることを特徴とする請求項１ないし６のいずれか１つに記載の磁気共鳴スペクトルの定量方法。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2005046419A JP4498947B2 (ja) | 2004-04-15 | 2005-02-23 | 磁気共鳴スペクトルの定量方法 |
| EP05252375A EP1586915A1 (en) | 2004-04-15 | 2005-04-15 | Method of quantifying a magnetic resonance spectrum |
| US11/107,356 US7106059B2 (en) | 2004-04-15 | 2005-04-15 | Method of quantifying magnetic resonance spectrum |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2004119750 | 2004-04-15 | ||
| JP2005046419A JP4498947B2 (ja) | 2004-04-15 | 2005-02-23 | 磁気共鳴スペクトルの定量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2005326397A JP2005326397A (ja) | 2005-11-24 |
| JP4498947B2 true JP4498947B2 (ja) | 2010-07-07 |
Family
ID=34940863
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2005046419A Expired - Fee Related JP4498947B2 (ja) | 2004-04-15 | 2005-02-23 | 磁気共鳴スペクトルの定量方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US7106059B2 (ja) |
| EP (1) | EP1586915A1 (ja) |
| JP (1) | JP4498947B2 (ja) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6349876B2 (ja) * | 2014-03-31 | 2018-07-04 | 住友大阪セメント株式会社 | 酸化ケイ素被覆酸化亜鉛とその製造方法及び酸化ケイ素被覆酸化亜鉛含有組成物並びに化粧料 |
| GB2529135A (en) * | 2014-06-21 | 2016-02-17 | Neil Francis Peirson | Method and apparatus for signal enhancement by reduction of noise |
| US10605751B2 (en) | 2015-02-09 | 2020-03-31 | Pepric Nv | System and method for determining a quantity of magnetic particles |
| CN106198600B (zh) * | 2015-05-05 | 2019-01-22 | 首都师范大学 | 一种磁共振检测溶液浓度的方法 |
| JP6787740B2 (ja) * | 2015-10-29 | 2020-11-18 | ローム株式会社 | オフセット導出装置、オフセット算出装置、及び、方位角センサ |
| CN111239657B (zh) * | 2020-01-20 | 2022-05-06 | 上海东软医疗科技有限公司 | 谱图的相位校正方法、装置及设备 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5614146A (en) * | 1979-07-16 | 1981-02-10 | Jeol Ltd | Sweep type nuclear magnetic resonator |
| JPS61210933A (ja) * | 1985-03-15 | 1986-09-19 | Jeol Ltd | ２次元核磁気共鳴測定方法 |
| JPS6373176A (ja) * | 1986-09-17 | 1988-04-02 | Hitachi Metals Ltd | 強磁性共鳴吸収の測定方法 |
Family Cites Families (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE1798079C3 (de) * | 1967-08-21 | 1975-03-06 | Varian Associates, Palo Alto, Calif. (V.St.A.) | Verfahren zur Messung der gyromagnetischen Resonanz und zur Durchführung des Verfahrens geeignetes Spinresonanzspektrometer |
| NL8203519A (nl) * | 1982-09-10 | 1984-04-02 | Philips Nv | Werkwijze en inrichting voor het bepalen van een kernmagnetisatieverdeling in een deel van een lichaam. |
| JPS59136642A (ja) | 1983-01-26 | 1984-08-06 | Jeol Ltd | 核磁気共鳴装置 |
| US4763282A (en) * | 1983-05-27 | 1988-08-09 | Larry Rosenberg | Programming format and apparatus for the improved coherent beam coupler system and method |
| FR2574551B1 (fr) * | 1984-12-12 | 1986-12-26 | Commissariat Energie Atomique | Procede de generation et de traitement de signaux pour l'obtention par resonance magnetique nucleaire d'une image exempte de distorsions a partir d'un champ de polarisation inhomogene |
| JPS63142249A (ja) | 1986-12-04 | 1988-06-14 | Jeol Ltd | 電子スピン共鳴装置 |
| DE3914301A1 (de) * | 1989-04-29 | 1990-10-31 | Bruker Medizintech | Verfahren zur aufnahme von spinresoneanzspektren und zur spinresonanz-bildgebung |
| US5572125A (en) * | 1991-03-25 | 1996-11-05 | Dunkel; Reinhard | Correction and automated analysis of spectral and imaging data |
| US5262723A (en) * | 1991-07-09 | 1993-11-16 | General Electric Company | Method and apparatus for obtaining pure-absorption two-dimensional lineshape data for multidimensional NMR spectroscopy using switched acquisition time gradients |
| US6472870B1 (en) * | 1999-02-23 | 2002-10-29 | M. Robin Bendall | Radiofrequency irradiation schemes and methods of design and display for use in performing nuclear magnetic resonance spectroscopy |
| JP2004032540A (ja) | 2002-06-27 | 2004-01-29 | Nippon Park:Kk | ホームコントロールシステム |
| CA2490305A1 (en) * | 2002-07-11 | 2004-01-22 | Thomas A. Szyperski | A method of using g-matrix fourier transformation nuclear magnetic resonance (gft nmr) spectroscopy for rapid chemical shift assignment and secondary structure determination of proteins |
| JP4006356B2 (ja) | 2003-04-28 | 2007-11-14 | 日本電子株式会社 | Ｎｍｒ測定条件の最適化方法 |
-
2005
- 2005-02-23 JP JP2005046419A patent/JP4498947B2/ja not_active Expired - Fee Related
- 2005-04-15 US US11/107,356 patent/US7106059B2/en active Active
- 2005-04-15 EP EP05252375A patent/EP1586915A1/en not_active Ceased
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5614146A (en) * | 1979-07-16 | 1981-02-10 | Jeol Ltd | Sweep type nuclear magnetic resonator |
| JPS61210933A (ja) * | 1985-03-15 | 1986-09-19 | Jeol Ltd | ２次元核磁気共鳴測定方法 |
| JPS6373176A (ja) * | 1986-09-17 | 1988-04-02 | Hitachi Metals Ltd | 強磁性共鳴吸収の測定方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| US20050237058A1 (en) | 2005-10-27 |
| US7106059B2 (en) | 2006-09-12 |
| EP1586915A1 (en) | 2005-10-19 |
| JP2005326397A (ja) | 2005-11-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4498947B2 (ja) | 磁気共鳴スペクトルの定量方法 | |
| US7191069B2 (en) | Automatic identification of compounds in a sample mixture by means of NMR spectroscopy | |
| US10775462B2 (en) | System and method for direct saturation-corrected chemical exchange saturation transfer (DISC-CEST) | |
| US6617167B2 (en) | Method of determining presence and concentration of lipoprotein X in blood plasma and serum | |
| Martineau et al. | Fast quantitative 2D NMR for metabolomics and lipidomics: A tutorial | |
| US9285444B2 (en) | Method for substance identification from NMR spectrum | |
| JP2006267111A (ja) | スペクトル、特にｎｍｒスペクトルのセットをプロセッシングする方法 | |
| Deelchand et al. | Sensitivity and specificity of human brain glutathione concentrations measured using short‐TE 1H MRS at 7 T | |
| US10613044B2 (en) | NMR quantification of TMAO | |
| EP2906956B1 (en) | Nmr quantification of branched chain amino acids | |
| Zu et al. | Chemical exchange rotation transfer imaging of phosphocreatine in muscle | |
| KR20190106305A (ko) | 와전류탐상검사를 이용하여 튜브 확관부의 결함 측정을 위한 대비시험편 및 이를 이용한 결함측정방법 | |
| Schoenberger et al. | Guideline for qNMR analysis | |
| Okada et al. | Repeatability of proton magnetic resonance spectroscopy of the brain at 7 T: effect of scan time on semi-localized by adiabatic selective refocusing and short-echo time stimulated echo acquisition mode scans and their comparison | |
| Mantel et al. | Study of liquid–liquid interfaces by an easily implemented localized NMR sequence | |
| Petrie et al. | Optimizing spherical navigator echoes for three‐dimensional rigid‐body motion detection | |
| US10036793B2 (en) | Method and apparatus for reconstructing magnetic resonance images with phase noise-dependent selection of raw data | |
| Riemann et al. | Assessment of measurement precision in single‐voxel spectroscopy at 7 T: Toward minimal detectable changes of metabolite concentrations in the human brain in vivo | |
| Perdigão et al. | Okapi-EM: A napari plugin for processing and analyzing cryogenic serial focused ion beam/scanning electron microscopy images | |
| US20240077439A1 (en) | Methods and Systems for Measuring Citrate and Creatinine Levels by NMR Spectroscopy | |
| US11835609B2 (en) | Double offsets and powers magnetization transfer ratio for specific magnetization transfer magnetic resonance imaging | |
| JP6346991B2 (ja) | Ｎｍｒ測定の結果に信号として含まれる情報を抽出するための方法 | |
| US9989486B2 (en) | Induction thermography method | |
| CN115079073B (zh) | 频率差分准静态磁感应成像方法、系统、设备以及介质 | |
| Pradhan et al. | Quantifying spectra in the frequency domain |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20071225 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20100406 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20100414 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130423 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 4498947 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130423 Year of fee payment: 3 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130423 Year of fee payment: 3 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130423 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140423 Year of fee payment: 4 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |