JP2741885B2 - 磁気共鳴を用いた検査装置におけるデータ処理方法 - Google Patents
磁気共鳴を用いた検査装置におけるデータ処理方法Info
- Publication number
- JP2741885B2 JP2741885B2 JP1036127A JP3612789A JP2741885B2 JP 2741885 B2 JP2741885 B2 JP 2741885B2 JP 1036127 A JP1036127 A JP 1036127A JP 3612789 A JP3612789 A JP 3612789A JP 2741885 B2 JP2741885 B2 JP 2741885B2
- Authority
- JP
- Japan
- Prior art keywords
- value
- estimated value
- signal
- estimated
- inspection apparatus
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Development Economics (AREA)
- High Energy & Nuclear Physics (AREA)
- General Health & Medical Sciences (AREA)
- Entrepreneurship & Innovation (AREA)
- Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Marketing (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Radiology & Medical Imaging (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Description
本発明な医療用の磁気共鳴を用いた検査装置に係り特
に得られたスペクトルの処理に好適なデータ処理方法に
関する。
に得られたスペクトルの処理に好適なデータ処理方法に
関する。
医療用の磁気共鳴検査装置において、NMRスペクトル
を計測することは診断に有効な情報を提供するものと期
待されている。しかしこの場合、対象とする特質の天然
存在比が小さい、計測にかけられる時間に限界がある、
などの理由でSN比の改善には限界がある。一方、医療用
のNMRスペクトロスコピーにおいては通常含まれるスペ
クトル成分の数、およびそれらの共鳴周波数はあらかじ
め知られている場合が多い。 したがって、これらの先験的な情報を利用することに
より、SN比の悪いスペクトルデータからでも診断に有効
な情報を取り出すことが考えられる。このような考え方
によるスペクトルデータの処理法がマグネティック・レ
ゾナンス・イン・メディスン3,97−104,1986（Magnetic
Resonance in Medicine3,97−104,1986）に論じられて
いる。 この方法によれば、含まれるスペクトル成分の数をK,
k番目のスペクトル成分の共鳴周波数をωｋ、振幅を
Ak、減衰乗数をbk、位相項をｋとすると、得られるFI
D Ｆ（ｔ）は、 と定式化できる。上論論文ではbk,ωｋがあらかじめ知
られているとして、ｋ及び診断に有効な情報を含んで
いると考えられるAkをFIDの実際の測定データＧ（ｔ）
から最小２乗法により決めている。 すなわち、k,Akに対する推定値を として、これらより計算したFIDを（ｔ）として、 となるようにAk,ｋを決める。ここでbk,ωｋは既知で
あるとする。この方法ではbk,ωｋを既知であるとして
いるので未知数として残るのはAk,ｋであり、十分な
測定点が得られれば、これを精度良く求めることができ
る。 この方法の問題点はωk,bkを既知としていることであ
る。共鳴周波数ωｋは知られている場合も多いが、局所
的なPHの影響をうける場合もあり、取り得る値の範囲の
みしか知り得ない場合もある。さらに、bkについては、
正確な値があらかじめわかるのはまれであり、多くの場
合、おおまかな値、あるいは値の範囲を知ることができ
るだけである。
を計測することは診断に有効な情報を提供するものと期
待されている。しかしこの場合、対象とする特質の天然
存在比が小さい、計測にかけられる時間に限界がある、
などの理由でSN比の改善には限界がある。一方、医療用
のNMRスペクトロスコピーにおいては通常含まれるスペ
クトル成分の数、およびそれらの共鳴周波数はあらかじ
め知られている場合が多い。 したがって、これらの先験的な情報を利用することに
より、SN比の悪いスペクトルデータからでも診断に有効
な情報を取り出すことが考えられる。このような考え方
によるスペクトルデータの処理法がマグネティック・レ
ゾナンス・イン・メディスン3,97−104,1986（Magnetic
Resonance in Medicine3,97−104,1986）に論じられて
いる。 この方法によれば、含まれるスペクトル成分の数をK,
k番目のスペクトル成分の共鳴周波数をωｋ、振幅を
Ak、減衰乗数をbk、位相項をｋとすると、得られるFI
D Ｆ（ｔ）は、 と定式化できる。上論論文ではbk,ωｋがあらかじめ知
られているとして、ｋ及び診断に有効な情報を含んで
いると考えられるAkをFIDの実際の測定データＧ（ｔ）
から最小２乗法により決めている。 すなわち、k,Akに対する推定値を として、これらより計算したFIDを（ｔ）として、 となるようにAk,ｋを決める。ここでbk,ωｋは既知で
あるとする。この方法ではbk,ωｋを既知であるとして
いるので未知数として残るのはAk,ｋであり、十分な
測定点が得られれば、これを精度良く求めることができ
る。 この方法の問題点はωk,bkを既知としていることであ
る。共鳴周波数ωｋは知られている場合も多いが、局所
的なPHの影響をうける場合もあり、取り得る値の範囲の
みしか知り得ない場合もある。さらに、bkについては、
正確な値があらかじめわかるのはまれであり、多くの場
合、おおまかな値、あるいは値の範囲を知ることができ
るだけである。
本発明の目的は前記bk,ωｋなどの値が正確に知られ
ていない場合でも、それらの取り得る値の範囲が知られ
ていれば、これを先験的な情報として利用することによ
り、SN比の悪いFIDからでも診断に有効な情報を取り出
すためのデータ処理法に関する。
ていない場合でも、それらの取り得る値の範囲が知られ
ていれば、これを先験的な情報として利用することによ
り、SN比の悪いFIDからでも診断に有効な情報を取り出
すためのデータ処理法に関する。
【課題を解決するための手段】 上記目的は（２）式で示された最適化演算にシミュレ
ーテッド アニーリングと呼ばれる手法を用い、しかも
そのコスト関数を工夫することにより達成される。以下
本発明を説明する。 まず、（１）式を次のように書きなおす。 ここで、Xkは推定値として求める変数であり、 振幅 Ak（ｋ＝1,…,K）はXk（ｋ＝1,…,K）で、 減衰乗数bk（ｋ＝1,…,K）はXk（ｋ＝Ｋ＋1,…,2K）
で、 周波数 ωｋ（ｋ＝1,…,K）はXk（ｋ＝2K＋1,…,3
K）で、 位相項 ｋ（ｋ＝1,…,K）はXk（ｋ＝3K＋1,…,4
K）で、 それぞれ表わす。おのおのの推定値は符号∧をつけて表
わす。 さて問題は、推定値ｋの取り得る値の範囲 Xk ν＜Xk＜Xk μ において を最小とするようなｋを求めることに帰着される。こ
こで（ｔ）は推定値ｋから（３）式を用いて、計算
により作り出したFIDであり、Ｇ（ｔ）はFIDの実測値で
ある。このようなXkはサイエンス,Vol.220,pp.671−68
0,1983（Science,Vol.220,pp.671−680,1983）に論じら
れているシミュレーティッド・アニーリングと呼ばれる
方法を注意深く行なうことにより求めることができる。
この方法を以下に説明する。 Xkに微小変位量ΔXkを与え、Ｅを（３），（４）式に
より求め、ΔXkを与える前のＥとの差ΔＥを計算する。
すなわち、 ΔＥ＝Ｅ（X1,…,Xk＋ΔXk,…,X4k） −Ｅ（X1,…,Xk,…,X4k） である。ここでΔＥ＜０であれば変位ΔXkを受け入れ
る。すなわち、Xk＋ΔXkを新しいXkとする。また、ΔＥ
＞０である場合は確立分布Ｐ（ΔＥ）＝exp（−ΔE/T）
に従がって変位ΔXkを受け入れるか、拒絶するかを決め
る。ここでＴは仮想的な温度と呼ばれる量であり、最初
十分高く取っておいて、上述のプロセスをｋ＝１→ｋ＝
4Kと何回も繰り返しながらゆっくりとＴを下げて行く。
もし十分ゆっくりと温度Ｔを下げるなら、Ｅの最小値と
して、局所的な最小値（local minimum）ではなく大域
的な最小値（global minimum）を与えるXk（ｋ＝1,…,4
Kのいずれか）を求めることができる。
ーテッド アニーリングと呼ばれる手法を用い、しかも
そのコスト関数を工夫することにより達成される。以下
本発明を説明する。 まず、（１）式を次のように書きなおす。 ここで、Xkは推定値として求める変数であり、 振幅 Ak（ｋ＝1,…,K）はXk（ｋ＝1,…,K）で、 減衰乗数bk（ｋ＝1,…,K）はXk（ｋ＝Ｋ＋1,…,2K）
で、 周波数 ωｋ（ｋ＝1,…,K）はXk（ｋ＝2K＋1,…,3
K）で、 位相項 ｋ（ｋ＝1,…,K）はXk（ｋ＝3K＋1,…,4
K）で、 それぞれ表わす。おのおのの推定値は符号∧をつけて表
わす。 さて問題は、推定値ｋの取り得る値の範囲 Xk ν＜Xk＜Xk μ において を最小とするようなｋを求めることに帰着される。こ
こで（ｔ）は推定値ｋから（３）式を用いて、計算
により作り出したFIDであり、Ｇ（ｔ）はFIDの実測値で
ある。このようなXkはサイエンス,Vol.220,pp.671−68
0,1983（Science,Vol.220,pp.671−680,1983）に論じら
れているシミュレーティッド・アニーリングと呼ばれる
方法を注意深く行なうことにより求めることができる。
この方法を以下に説明する。 Xkに微小変位量ΔXkを与え、Ｅを（３），（４）式に
より求め、ΔXkを与える前のＥとの差ΔＥを計算する。
すなわち、 ΔＥ＝Ｅ（X1,…,Xk＋ΔXk,…,X4k） −Ｅ（X1,…,Xk,…,X4k） である。ここでΔＥ＜０であれば変位ΔXkを受け入れ
る。すなわち、Xk＋ΔXkを新しいXkとする。また、ΔＥ
＞０である場合は確立分布Ｐ（ΔＥ）＝exp（−ΔE/T）
に従がって変位ΔXkを受け入れるか、拒絶するかを決め
る。ここでＴは仮想的な温度と呼ばれる量であり、最初
十分高く取っておいて、上述のプロセスをｋ＝１→ｋ＝
4Kと何回も繰り返しながらゆっくりとＴを下げて行く。
もし十分ゆっくりと温度Ｔを下げるなら、Ｅの最小値と
して、局所的な最小値（local minimum）ではなく大域
的な最小値（global minimum）を与えるXk（ｋ＝1,…,4
Kのいずれか）を求めることができる。
以下、本発明の一実施例を、図面により説明する。 第１図は、本発明の一実施例による核磁共鳴を用いた
検査装置（以下、単に「検査装置」と呼ぶ）の概略構成
図である。 第１図において、１は静磁場Hoを発生させる電磁石、
２は対象物体、３は高周波磁場を発生させると同時に、
対象物体２から生ずる信号を検出する信号を検出するた
めのコイル、4x,4y、および５は、それぞれＸ方向,Y方
向,Z方向の傾斜磁場を発生させるための傾斜磁場発生コ
イルである。傾斜磁場発生コイル５としては、互いに逆
格納するメモリーである。 さて、第１図に示す装置構成において、NMRスペクト
ルを測定する手法はいくつか提案されている。代表的な
ものはISISと呼ばれる手法であり、ジャーナル オブ
マグネティック レゾナンス 66,283−294,1986（Jour
nal of Magnetic Reconance 66,283−294,1986）に記載
されている。この手法を実現するためのRFフィールドと
傾斜磁場の印加シークエンスを第２図に示す。 このISISにおいては、90゜パルス前に印加するX,Y,Z
の選択照射は必ずしも毎回すべて行なうものではなく、
表１に示すような８通りの組み合わせで行なう。さて、
表１に示す８通りの選択により得られた90゜パルス直後
のFIDを表１に示す奇数番目のものについてはそのま
ま、表１に示す偶数番目のものについては−１倍して加
え合わせると希望する選択領域からの信号のみを取り出
すことができる。 さて、このようにして取り出された複素信号Ｇ（ｔ）
スペクトル情報を前述のシミュレーテッド アニーリン
グを行ない取り出す。第ｋ番目のスペクトル成分の振
幅，減衰乗数，共鳴周波数，位相項の推定値をそれぞれ とする。スペクトル成分の数をＫとして、これらの推定
値から計算された推定FIDを（ｔ）と書くと、 と表わされる。ここで、 なるＥの最小値を与える （ｋ＝1,…,K）を以下の様に求める。まず、 振幅 Ak（ｋ＝1,…,K）はXk（ｋ＝1,…,K）で、 減衰乗数bk（ｋ＝1,…,K）はXk（ｋ＝Ｋ＋1,…,2K）
で、 周波数 ωｋ（ｋ＝1,…,K）はXk（ｋ＝2K＋1,…,3
K）で、 位相項 ｋ（ｋ＝1,…,K）はXk（ｋ＝3K＋1,…,4
K）で、 それぞれ表わす。おのおののXk＝1,…,K）に対しその変
位量ΔXkを設定する。これはXkの推定値に対し、1/100
〜1/500程度に設定する。Xkがおおよそどの程度の値を
取るか全く分からない場合十分に小さく設定する。 まず十分大きなＴを設定し、第３図のフローチャート
にしたがってアニーリングを行なう。Ｔが十分大きいか
否かは、試しにΔＥを適当に与えた初期値により計算
し、ｐ＝ｅ−ΔE/T0.9〜0.95である事を確認すればよ
い。仮想的なエネルギーＥは以下の式により計算する。 ここで、（７）式の右辺第２項EL（X1,…,Xk）は各Xk
があらかじめ設定されている範囲 Xk νXkXk μ をはずれた時非常に大きな値を取るとする。すなわち、 EL（X,…,Xk）＝∞ Xk ν＞XkorXk＞Xk μ ＝０ Xk νXkXk（ｋ＝１μ，…,K） である。 ある温度Ｔにおいて、熱平衡となった場合に次のＴへ
移る。熱平衡の条件はN1＝N2であるが実際には|N1−N2|
/N1＜0.1〜0.05であれば熱平衡とみなす。次の温度へ移
るときの係数ξはなるべく１に近いことがのぞましい。
１に非常に近いと計算時間を要するのでこれとのかねあ
いであろう。前記サイエンスの文献などではξ＝0.9を
推奨している。第３図でＫが小さい場合、熱平衡である
か否かがよく分からない。したがってＭのループの部分
をMmax回繰り返す。K.Mmax〜100から200であれば十分で
あろう。N1＝０となったらアニーリング終了とする。 以上、FIDのモデル化に（１）式を用いた場合につい
ては説明したが、位相誤差ｋは時間ずれτと全成分に
共通な位相回りp0に集約される。したがってFID Ｆ
（ｔ）は とも書ける。この場合は推定すべき量は である。 さらにこの方法はNMRスペクトルスコピックイメージ
ングにも適用できる。第４図にこの場合の展型的なシー
クエンスを示す。図に示される通り、このスペクトロス
コピックイメージングの場合にはデータ取得はRFパルス
による励起から位相エンコードなどに必要な時間tCだけ
待ってから行なわなければならない。したがって従来の
処理法ではT2のみじかい成分がTC時間内に減衰してしま
い、小さな振幅となると言う欠点があった。本手法を用
いることにより、この欠点を克服することができる。 すなわち、空間周波数成分について２次元フーリエ変
換後の時間波形を とモデル化することにより求めたAk（x,y）の推定値はt
C時間内の減衰を含まないものとなる。
検査装置（以下、単に「検査装置」と呼ぶ）の概略構成
図である。 第１図において、１は静磁場Hoを発生させる電磁石、
２は対象物体、３は高周波磁場を発生させると同時に、
対象物体２から生ずる信号を検出する信号を検出するた
めのコイル、4x,4y、および５は、それぞれＸ方向,Y方
向,Z方向の傾斜磁場を発生させるための傾斜磁場発生コ
イルである。傾斜磁場発生コイル５としては、互いに逆
格納するメモリーである。 さて、第１図に示す装置構成において、NMRスペクト
ルを測定する手法はいくつか提案されている。代表的な
ものはISISと呼ばれる手法であり、ジャーナル オブ
マグネティック レゾナンス 66,283−294,1986（Jour
nal of Magnetic Reconance 66,283−294,1986）に記載
されている。この手法を実現するためのRFフィールドと
傾斜磁場の印加シークエンスを第２図に示す。 このISISにおいては、90゜パルス前に印加するX,Y,Z
の選択照射は必ずしも毎回すべて行なうものではなく、
表１に示すような８通りの組み合わせで行なう。さて、
表１に示す８通りの選択により得られた90゜パルス直後
のFIDを表１に示す奇数番目のものについてはそのま
ま、表１に示す偶数番目のものについては−１倍して加
え合わせると希望する選択領域からの信号のみを取り出
すことができる。 さて、このようにして取り出された複素信号Ｇ（ｔ）
スペクトル情報を前述のシミュレーテッド アニーリン
グを行ない取り出す。第ｋ番目のスペクトル成分の振
幅，減衰乗数，共鳴周波数，位相項の推定値をそれぞれ とする。スペクトル成分の数をＫとして、これらの推定
値から計算された推定FIDを（ｔ）と書くと、 と表わされる。ここで、 なるＥの最小値を与える （ｋ＝1,…,K）を以下の様に求める。まず、 振幅 Ak（ｋ＝1,…,K）はXk（ｋ＝1,…,K）で、 減衰乗数bk（ｋ＝1,…,K）はXk（ｋ＝Ｋ＋1,…,2K）
で、 周波数 ωｋ（ｋ＝1,…,K）はXk（ｋ＝2K＋1,…,3
K）で、 位相項 ｋ（ｋ＝1,…,K）はXk（ｋ＝3K＋1,…,4
K）で、 それぞれ表わす。おのおののXk＝1,…,K）に対しその変
位量ΔXkを設定する。これはXkの推定値に対し、1/100
〜1/500程度に設定する。Xkがおおよそどの程度の値を
取るか全く分からない場合十分に小さく設定する。 まず十分大きなＴを設定し、第３図のフローチャート
にしたがってアニーリングを行なう。Ｔが十分大きいか
否かは、試しにΔＥを適当に与えた初期値により計算
し、ｐ＝ｅ−ΔE/T0.9〜0.95である事を確認すればよ
い。仮想的なエネルギーＥは以下の式により計算する。 ここで、（７）式の右辺第２項EL（X1,…,Xk）は各Xk
があらかじめ設定されている範囲 Xk νXkXk μ をはずれた時非常に大きな値を取るとする。すなわち、 EL（X,…,Xk）＝∞ Xk ν＞XkorXk＞Xk μ ＝０ Xk νXkXk（ｋ＝１μ，…,K） である。 ある温度Ｔにおいて、熱平衡となった場合に次のＴへ
移る。熱平衡の条件はN1＝N2であるが実際には|N1−N2|
/N1＜0.1〜0.05であれば熱平衡とみなす。次の温度へ移
るときの係数ξはなるべく１に近いことがのぞましい。
１に非常に近いと計算時間を要するのでこれとのかねあ
いであろう。前記サイエンスの文献などではξ＝0.9を
推奨している。第３図でＫが小さい場合、熱平衡である
か否かがよく分からない。したがってＭのループの部分
をMmax回繰り返す。K.Mmax〜100から200であれば十分で
あろう。N1＝０となったらアニーリング終了とする。 以上、FIDのモデル化に（１）式を用いた場合につい
ては説明したが、位相誤差ｋは時間ずれτと全成分に
共通な位相回りp0に集約される。したがってFID Ｆ
（ｔ）は とも書ける。この場合は推定すべき量は である。 さらにこの方法はNMRスペクトルスコピックイメージ
ングにも適用できる。第４図にこの場合の展型的なシー
クエンスを示す。図に示される通り、このスペクトロス
コピックイメージングの場合にはデータ取得はRFパルス
による励起から位相エンコードなどに必要な時間tCだけ
待ってから行なわなければならない。したがって従来の
処理法ではT2のみじかい成分がTC時間内に減衰してしま
い、小さな振幅となると言う欠点があった。本手法を用
いることにより、この欠点を克服することができる。 すなわち、空間周波数成分について２次元フーリエ変
換後の時間波形を とモデル化することにより求めたAk（x,y）の推定値はt
C時間内の減衰を含まないものとなる。
以上のように本発明によれば、核磁気共鳴信号の各ス
ペクトルの共鳴周波数，減衰乗数が正確に特定できない
場合にも、SN比の悪いFIDから診断に有効な情報を取り
出すことができる。
ペクトルの共鳴周波数，減衰乗数が正確に特定できない
場合にも、SN比の悪いFIDから診断に有効な情報を取り
出すことができる。
【図面の簡単な説明】 第１図は本発明を実施する際の装置構成の一例を示す。 第２図はある特定の領域からのみNMRスペクトルを観測
するISISと呼ばれる方法を実施するためのパルスシーケ
ンスの一例を示す。 第３図はシミュレーテッドアニーリングを実施する場合
のフローチャートを示す。 第４図はスペクトロスコピックイメージングを実施する
際の展型的なパルスシーケンスを示す。
するISISと呼ばれる方法を実施するためのパルスシーケ
ンスの一例を示す。 第３図はシミュレーテッドアニーリングを実施する場合
のフローチャートを示す。 第４図はスペクトロスコピックイメージングを実施する
際の展型的なパルスシーケンスを示す。
Claims (6)
- 【請求項１】静磁場、傾斜磁場、及び高周波磁場発生手
段と、検査対象からの核磁気共鳴信号を検出する信号検
出手段と、該信号検出手段による検出信号の演算を行な
う演算処理手段と、該演算処理手段による処理結果を出
力する出力手段とを有する核磁気共鳴を用いた検査装置
におけるデータ処理方法において、高周波磁場印加後に
観測された磁化信号を、該磁化信号に含まれる各スペク
トル成分の振幅、減衰乗数、共鳴周波数、及び位相項の
推定値から計算により求めた推定磁化信号と比較し、前
記観測された磁化信号と前記推定磁化信号との差の絶対
値、又は該絶対値の２乗の時間軸方向での和Ｅが最小と
なるように、前記振幅、前記減衰乗数、前記共鳴周波
数、及び前記位相項を求めることを特徴とする核磁気共
鳴を用いた検査装置におけるデータ処理方法。 - 【請求項２】請求項１に記載の方法において、前記各ス
ペクトル成分の前記推定値に順次変位を与え新しい推定
値毎に前記Ｅを計算して、前記推定値に前記変位を与え
る以前に計算された前記Ｅとの差ΔＥが、負でありかつ
前記新しい推定値が予め設定された領域内にある場合
に、推定値を前記変位を変えた後の値で置き換えること
を特徴とする核磁気共鳴を用いた検査装置におけるデー
タ処理方法。 - 【請求項３】請求項１に記載の方法において、前記各ス
ペクトル成分の前記推定値に順次変位を与え新しい推定
値毎に前記Ｅを計算して、前記推定値に前記変位を与え
る以前に計算された前記Ｅとの差ΔＥが、正でありかつ
前記新しい推定値が予め設定された領域内にある場合
に、前記変位を前記ΔＥの関数である確率Ｐ（ΔＥ）で
推定値を前記変位を与えた後の値で置き換えることを特
徴とする核磁気共鳴を用いた検査装置におけるデータ処
理方法。 - 【請求項４】請求項１に記載の方法において、前記Ｐ
（ΔＥ）はexp（−ΔE/T）で与えられ、Ｔを順次減少さ
せていくことを特徴とする核磁気共鳴を用いた検査装置
におけるデータ処理方法。 - 【請求項５】請求項３に記載の方法において、前記新し
い推定値が予め設定された領域外にある場合に、推定値
を前記変位を与えた後の値で置き換えないことを特徴と
する核磁気共鳴を用いた検査装置におけるデータ処理方
法。 - 【請求項６】静磁場、傾斜磁場、及び高周波磁場発生手
段と、検査対象からの核磁気共鳴信号を検出する信号検
出手段と、該信号検出手段による検出信号の演算を行な
う演算処理手段と、該演算処理手段による処理結果を出
力する出力手段とを有する核磁気共鳴を用いた検査装置
におけるデータ処理方法において、高周波磁場印加後に
観測された磁化信号を、該磁化信号に含まれる各スペク
トル成分の振幅、減衰乗数、共鳴周波数、及び位相項の
推定値から計算により求めた推定磁化信号と比較し、前
記観測された磁化信号と前記推定磁化信号との差の絶対
値、又は該絶対値の２乗の時間軸方向での和Ｅが最小と
なるように、前記各スペクトル成分の前記推定値に順次
変位を与え新しい推定値毎に前記Ｅを計算して、前記推
定値に前記変位を与える以前に計算された前記Ｅとの差
ΔＥが、負でありかつ前記新しい推定値が予め設定され
た領域内にある場合に、推定値を前記変位を与えた後の
値で置き換え、前記ΔＥが、正でありかつ前記新しい推
定値が予め設定された領域内にある場合に、前記変位を
前記ΔＥの関数で与えられる確率Ｐ（ΔＥ）＝exp（−
ΔE/T）で推定値を前記変位を与えた後の値で置き換え
Ｔを順次減少させていき、前記新しい推定値が予め設定
された領域外にある場合に、推定値を前記変位を与えた
後の値で置き換えないようにして、前記振幅、前記減衰
乗数、前記共鳴周波数、及び前記位相項を求めることを
特徴とする核磁気共鳴を用いた検査装置におけるデータ
処理方法。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP1036127A JP2741885B2 (ja) | 1989-02-17 | 1989-02-17 | 磁気共鳴を用いた検査装置におけるデータ処理方法 |
| US07/480,946 US5218531A (en) | 1989-02-17 | 1990-02-16 | NMR spectroscopic analyzing method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP1036127A JP2741885B2 (ja) | 1989-02-17 | 1989-02-17 | 磁気共鳴を用いた検査装置におけるデータ処理方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH02215441A JPH02215441A (ja) | 1990-08-28 |
| JP2741885B2 true JP2741885B2 (ja) | 1998-04-22 |
Family
ID=12461118
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP1036127A Expired - Lifetime JP2741885B2 (ja) | 1989-02-17 | 1989-02-17 | 磁気共鳴を用いた検査装置におけるデータ処理方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US5218531A (ja) |
| JP (1) | JP2741885B2 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4617444B2 (ja) * | 2004-03-31 | 2011-01-26 | 株式会社国際電気通信基礎技術研究所 | 脳内電流源推定方法、脳内電流源推定プログラムおよび脳内電流源推定装置 |
| US10564240B2 (en) * | 2015-02-06 | 2020-02-18 | The University Of Melbourne | Method and system for magnetic resonance imaging using slice selective pulses |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS59200947A (ja) * | 1983-04-30 | 1984-11-14 | Toshiba Corp | Ｍｒｉ装置における静磁場強度調整方法 |
| JPS6321048A (ja) * | 1986-07-16 | 1988-01-28 | 横河メディカルシステム株式会社 | 核磁気共鳴画像の歪み補正方法 |
-
1989
- 1989-02-17 JP JP1036127A patent/JP2741885B2/ja not_active Expired - Lifetime
-
1990
- 1990-02-16 US US07/480,946 patent/US5218531A/en not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JPH02215441A (ja) | 1990-08-28 |
| US5218531A (en) | 1993-06-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4444413B2 (ja) | ４重フィールドエコーシーケンスを用いて水と脂肪を定量的にｍｒ撮影する装置 | |
| EP2539729B1 (en) | Magnetic resonance imaging of chemical species with a spectral model | |
| CN104422916B (zh) | 借助两点迪克松技术进行磁共振测量的方法和磁共振设备 | |
| WO2016184779A1 (en) | Time-domain mri | |
| EP0145276B1 (en) | Complex quotient nuclear magnetic resonance imaging | |
| US4654595A (en) | Method of analyzing properties of a matter or a magnetic field | |
| JPH074350B2 (ja) | Ｎｍｒ信号受信方法 | |
| JP6581584B2 (ja) | 水／脂肪分離を用いた位相感応型反転回復ｍｒｉ | |
| EP2402780A1 (en) | Method of reconstructing a magnetic resonance image of an object considering higher-order dynamic fields | |
| WO2009070695A1 (en) | System and method for phase offset and time delay correction in magnetic resonance spectroscopy data | |
| US20100305424A1 (en) | K-space sample density compensation for magnetic resonance image reconstruction | |
| EP0208483A2 (en) | Method and apparatus for rapid NMR imaging | |
| CN103885016A (zh) | 相位差图的确定 | |
| EP0100183A2 (en) | Nuclear magnetic resonance method and apparatus | |
| US20050137476A1 (en) | Method for evaluating magnetic resonance spectroscopy data using a baseline model | |
| Horsfield et al. | Two-species chemical-shift imaging using prior knowledge and estimation theory. Application to rock cores | |
| US4706027A (en) | Method for correcting phase errors in magnetic resonance imaging data | |
| JPH03264046A (ja) | 核磁気共鳴映像法及び装置 | |
| NL8801731A (nl) | Werkwijze en kernspinresonantie-inrichting om de dwarsrelaxatietijdkonstante t2 snel te bepalen. | |
| JP2741885B2 (ja) | 磁気共鳴を用いた検査装置におけるデータ処理方法 | |
| EP0213614A2 (en) | Method for obtaining nuclear magnetic resonance information data | |
| GB2139765A (en) | Nuclear magnetic resconance spectrum at a point with a sample | |
| US11385311B2 (en) | System and method for improved magnetic resonance fingerprinting using inner product space | |
| JP4187290B2 (ja) | 容積磁気共鳴分光方法及び装置 | |
| JP2001281315A (ja) | 磁気共鳴装置の運転方法 |