JP5479427B2 - 磁気共鳴撮影装置 - Google Patents
磁気共鳴撮影装置 Download PDFInfo
- Publication number
- JP5479427B2 JP5479427B2 JP2011222934A JP2011222934A JP5479427B2 JP 5479427 B2 JP5479427 B2 JP 5479427B2 JP 2011222934 A JP2011222934 A JP 2011222934A JP 2011222934 A JP2011222934 A JP 2011222934A JP 5479427 B2 JP5479427 B2 JP 5479427B2
- Authority
- JP
- Japan
- Prior art keywords
- magnetic resonance
- signal
- resonance signal
- real space
- space
- 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.)
- Active
Links
- 238000002595 magnetic resonance imaging Methods 0.000 title claims description 44
- 238000005259 measurement Methods 0.000 claims description 185
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 113
- 238000012937 correction Methods 0.000 claims description 93
- 230000001629 suppression Effects 0.000 claims description 88
- 238000001228 spectrum Methods 0.000 claims description 68
- 230000015572 biosynthetic process Effects 0.000 claims description 45
- 238000003786 synthesis reaction Methods 0.000 claims description 45
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 claims description 26
- 230000003595 spectral effect Effects 0.000 claims description 23
- 238000004364 calculation method Methods 0.000 claims description 22
- 230000003068 static effect Effects 0.000 claims description 20
- 238000009826 distribution Methods 0.000 claims description 19
- 238000006243 chemical reaction Methods 0.000 claims description 15
- 238000001514 detection method Methods 0.000 claims description 15
- 230000035945 sensitivity Effects 0.000 claims description 12
- 238000005481 NMR spectroscopy Methods 0.000 claims description 11
- 230000036278 prepulse Effects 0.000 claims description 10
- 238000000034 method Methods 0.000 description 39
- 230000005415 magnetization Effects 0.000 description 16
- 230000005284 excitation Effects 0.000 description 15
- 238000012545 processing Methods 0.000 description 13
- 230000000694 effects Effects 0.000 description 12
- 239000002207 metabolite Substances 0.000 description 10
- 230000008569 process Effects 0.000 description 10
- 238000003384 imaging method Methods 0.000 description 9
- 230000005540 biological transmission Effects 0.000 description 8
- 238000003860 storage Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 239000000126 substance Substances 0.000 description 6
- 230000007423 decrease Effects 0.000 description 4
- 230000006872 improvement Effects 0.000 description 4
- 230000002194 synthesizing effect Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 230000005389 magnetism Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000007781 pre-processing Methods 0.000 description 3
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 2
- 239000007864 aqueous solution Substances 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000005281 excited state Effects 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 238000000655 nuclear magnetic resonance spectrum Methods 0.000 description 2
- 238000004904 shortening Methods 0.000 description 2
- OTCCIMWXFLJLIA-UHFFFAOYSA-N N-acetyl-DL-aspartic acid Natural products CC(=O)NC(C(O)=O)CC(O)=O OTCCIMWXFLJLIA-UHFFFAOYSA-N 0.000 description 1
- KTHDTJVBEPMMGL-VKHMYHEASA-N N-acetyl-L-alanine Chemical compound OC(=O)[C@H](C)NC(C)=O KTHDTJVBEPMMGL-VKHMYHEASA-N 0.000 description 1
- OTCCIMWXFLJLIA-BYPYZUCNSA-N N-acetyl-L-aspartic acid Chemical compound CC(=O)N[C@H](C(O)=O)CC(O)=O OTCCIMWXFLJLIA-BYPYZUCNSA-N 0.000 description 1
- 108091081062 Repeated sequence (DNA) Proteins 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 238000000701 chemical imaging Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000005429 filling process Methods 0.000 description 1
- 235000015220 hamburgers Nutrition 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000005236 sound signal Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
- 238000001308 synthesis method Methods 0.000 description 1
- 230000002123 temporal effect 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/46—NMR spectroscopy
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Radiology & Medical Imaging (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Description
以下、本発明を適用した第一の実施の形態を図面を参照して説明する。
φ(t)=tan−1(Im(F(t))/Re(F(t))) (式1)
ここで、tan−1はアークタンジェント関数、Im(X)は複素数Xの虚部、Re(X)は複素数Xの実部を表す。
φ(i)=tan−1(Im(S(i))/Re(S(i))) (式2)
P(i)=C1×i+C0 (式3)
位相補正は、(式2)で求められる信号位相φ(i)から、(式3)で求められる位相歪みP(i)分を減算することにより行う。
φw(L(i))(x(a),y(b),z(c))(t(N))=tan−1{Im(Sw2(L(i))(x(a),y(b),z(c))(t(N)))/Re(Sw2(L(i))(x(a),y(b),z(c))(t(N)))} (式4)
なお、求めた小型コイルL(i)毎の信号位相値φw(L(i))(x(a),y(b),z(c))(t(N))は、本計測(水抑圧計測)で計測する信号を補正する補正値、すなわち、MAC合成用補正値として記憶装置11に記憶する。この信号位相値φw(L(i))(x(a),y(b),z(c))(t(N))は、上述のデータの番号iに依存しない信号位相φ(i)のオフセット変動成分C0に相当する。
Re(Sm3(L(i))(x(a),y(b),z(c))(t(j)))=cos((−1)×π×φw(L(i))(x(a),y(b),z(c))(t(N))/180)×Re(Sm2(L(i))(x(a),y(b),z(c))(t(j)))−sin((−1)×π×φw(L(i))(x(a),y(b),z(c))(t(N))/180)×Im(Sm2(L(i))(x(a),y(b),z(c))(t(j)))
(式5)
Im(Sm3(L(i))(x(a),y(b),z(c))(t(j)))=sin((−1)×π×φw(L(i))(x(a),y(b),z(c))(t(N))/180)×Re(Sm2(L(i))(x(a),y(b),z(c))(t(j)))+cos((−1)×π×φw(L(i))(x(a),y(b),z(c))(t(N))/180)×Im(Sm2(L(i))(x(a),y(b),z(c))(t(j)))
(式6)
次に本発明の第二の実施の形態を説明する。本実施形態においても、第一の実施形態と同様に、非水抑圧計測で得られた小型コイル毎の位相補正値を用いて、その後に行われる水抑圧計測で得られる信号を小型コイル毎に位相補正する。ただし、本実施形態では、水抑圧計測(本計測)による信号を足し合わせる際(MAC合成の際)、各小型コイルのSNRに応じた異なる加算係数を乗算する。なお、第一の実施形態は、本実施形態の全ての小型コイルの加算係数を常に1にしたものといえる。以下、第一の実施形態と異なる点についてのみ説明する。
SNR(L(i))(x(a),y(b),z(c))=(Aw(L(i))(x(a),y(b),z(c))(N)−Na(L(i))(x(a),y(b),z(c)))/Nd(L(i))(x(a),y(b),z(c)) (式13)
SNR(L(1))(x(0),y(0),z(0))=e (式14)
SNR(L(2))(x(0),y(0),z(0))=f (式15)
SNR(L(3))(x(0),y(0),z(0))=g (式16)
SNR(L(4))(x(0),y(0),z(0))=h (式17)
R(L(1))(x(0),y(0),z(0))=e/max(e,f,g,h)=E (式18)
R(L(2))(x(0),y(0),z(0))=f/max(e,f,g,h)=F (式19)
R(L(3))(x(0),y(0),z(0))=g/max(e,f,g,h)=G (式20)
R(L(4))(x(0),y(0),z(0))=h/max(e,f,g,h)=H (式21)
なお、max(e,f,g,h)は、e,f,g,hのうちの最大の値を取る関数である。
W(L(1))(x(0),y(0),z(0))=E×P/(E+F+G+H) (式22)
W(L(2))(x(0),y(0),z(0))=F×P/(E+F+G+H) (式23)
W(L(3))(x(0),y(0),z(0))=G×P/(E+F+G+H) (式24)
W(L(4))(x(0),y(0),z(0))=H×P/(E+F+G+H) (式25)
Re(Sm3(L(i))(x(a),y(b),z(c))(t(j)))=cos((−1)×π×φw(L(i))(x(a),y(b),z(c))(t(N))/180)×Re(Sm2(L(i))(x(a),y(b),z(c))(t(j)))−sin((−1)×π×φw(L(i))(x(a),y(b),z(c))(t(N))/180)×Im(Sm2(L(i))(x(a),y(b),z(c))(t(j)))
(式26)
Im(Sm3(L(i))(x(a),y(b),z(c))(t(j)))=sin((−1)×π×φw(L(i))(x(a),y(b),z(c))(t(N))/180)×Re(Sm2(L(i))(x(a),y(b),z(c))(t(j)))+cos((−1)×π×φw(L(i))(x(a),y(b),z(c))(t(N))/180)×Im(Sm2(L(i))(x(a),y(b),z(c))(t(j)))
(式27)
Re(Sm4(L(i))(x(a),y(b),z(c))(t(j)))=Re(Sm3(L(i))(x(a),y(b),z(c))(t(j)))×W(L(i))(x(a),y(b),z(c))(t(N))
(式28)
Im(Sm4(L(i))(x(a),y(b),z(c))(t(j)))=Re(Sm3(L(i))(x(a),y(b),z(c))(t(j)))×W(L(i))(x(a),y(b),z(c))(t(N))
(式29)
Claims (7)
- 静磁場、高周波磁場および傾斜磁場をそれぞれ発生する磁場発生手段と、前記静磁場中に置かれた被検体から発生する磁気共鳴信号を検出する検出手段と、前記磁場発生手段および検出手段を制御する計測制御手段と、前記核磁気共鳴信号を用いて磁気共鳴スペクトル画像を作成し、表示する演算手段とを備え、前記検出手段が複数の小型コイルからなる受信コイルを有する磁気共鳴撮影装置であって、
前記受信コイルは複数の小型コイルL(i)(ここで、iはコイル番号を表す1以上I以下の整数:i=1,2,…,I)を備えるマルチアレイコイル(MAC)であり、
前記計測制御手段は、
非水抑圧時の、2次元k空間もしくは3次元k空間(kx(a),ky(b),kz(c))(ここで、aはkx軸方向の計測番号を表すA以下の整数、bはky軸方向の計測番号を表すB以下の整数、cはkz軸方向の計測番号を表すC以下の整数、2次元の場合はA,B,Cのいずれか1つが1で他の2つは2以上の整数、3次元の場合はA,B,Cのいずれもが2以上の整数)の前記小型コイル毎の磁気共鳴信号Sw1(L(i))(kx(a),ky(b),kz(c))(t(j))(ここで、jは時間軸(t軸)方向のデータ番号を表すJ以下の整数)を計測する第一の計測手段と、
水抑圧時の、前記2次元k空間もしくは3次元k空間(kx(a),ky(b),kz(c))の前記小型コイル毎の磁気共鳴信号Sm1(L(i))(kx(a),ky(b),kz(c))(t(j))を計測する第二の計測手段と、を備え、
前記演算手段は、
前記磁気共鳴信号Sw1(L(i))(kx(a),ky(b),kz(c))(t(j))に対して、各k空間軸(kx,ky,kz)方向に逆フーリエ変換を施し、2次元実空間もしくは3次元実空間(x(a),y(b),z(c))の磁気共鳴信号Sw2(L(i))(x(a),y(b),z(c))(t(j))を算出する第一のk空間−実空間変換手段と、
前記各小型コイルの実空間点毎に、前記磁気共鳴信号Sw2(L(i))(x(a),y(b),z(c))(t(j))の先頭からN番目の信号Sw2(L(i))(x(a),y(b),z(c))(t(N))の位相値φw(L(i))(x(a),y(b),z(c))(t(N))を算出し、算出結果から補正値を求める位相値算出手段と、
前記第二の計測手段が計測した前記小型コイル毎のk空間の磁気共鳴信号Sm1(L(i))(kx(a),ky(b),kz(c))(t(j))に対して、各k空間軸方向に逆フーリエ変換を施し、2次元実空間もしくは3次元実空間(x(a),y(b),z(c))の磁気共鳴信号Sm2(L(i))(x(a),y(b),z(c))(t(j))を算出する第二のk空間−実空間変換手段と、
各小型コイルの実空間点毎に、前記磁気共鳴信号Sm2(L(i))(x(a),y(b),z(c))(t(j))の全点の前記位相値φm(L(i))(x(a),y(b),z(c))(t(j))を前記補正値を用いて位相補正し、前記小型コイル毎の位相補正後の磁気共鳴信号Sm3(L(i))(x(a),y(b),z(c))(t(j))を得る位相値補正手段と、
前記磁気共鳴信号Sm3(L(i))(x(a),y(b),z(c))(t(j))を、実空間点毎に複素数加算し、信号Sm4(x(a),y(b),z(c))(t(j))を得るMAC合成手段と、
前記磁気共鳴信号Sm4(x(a),y(b),z(c))(t(j))に対して、時間軸方向にフーリエ変換を施し、周波数軸(f軸)方向の情報を有する2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像Sm5(x(a),y(b),z(c))(f(j))を算出する時間−周波数変換手段と、を備え、
前記第一の計測手段は、MRSIシーケンスを実行することにより、前記磁気共鳴信号Sw1を得、
前記第二の計測手段は、前記MRSIシーケンスと、前記MRSIシーケンスに先立ち、水信号を抑圧するプリパルスシーケンスとを実行することにより、前記磁気共鳴信号Sm1を得ること
を特徴とする磁気共鳴撮影装置。 - 請求項1記載の磁気共鳴撮影装置であって、
前記演算手段は、前記磁気共鳴スペクトル画像Sm5(x(a),y(b),z(c))(f(j))に対してスペクトル軸方向の位相補正を施し、2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像Sm6(x(a),y(b),z(c))(f(j))を算出する第二の位相値補正手段をさらに備えること
を特徴とする磁気共鳴撮影装置。 - 静磁場、高周波磁場および傾斜磁場をそれぞれ発生する磁場発生手段と、前記静磁場中に置かれた被検体から発生する磁気共鳴信号を検出する検出手段と、前記磁場発生手段および検出手段を制御する計測制御手段と、前記核磁気共鳴信号を用いて磁気共鳴スペクトル画像を作成し、表示する演算手段とを備え、前記検出手段が複数の小型コイルからなる受信コイルを有する磁気共鳴撮影装置であって、
前記受信コイルは複数の小型コイルL(i)(ここで、iはコイル番号を表す1以上I以下の整数:i=1,2,…,I)を備えるマルチアレイコイル(MAC)であり、
前記計測制御手段は
非水抑圧時の、2次元k空間もしくは3次元k空間(kx(a),ky(b),kz(c))(ここで、aはkx軸方向の計測番号を表すA以下の整数、bはky軸方向の計測番号を表すB以下の整数、cはkz軸方向の計測番号を表すC以下の整数、2次元の場合はA,B,Cのいずれか1つが1で他の2つは2以上の整数、3次元の場合はA,B,Cのいずれもが2以上の整数)の前記小型コイル毎の磁気共鳴信号Sw1(L(i))(kx(a),ky(b),kz(c))(t(j))(ここで、jは時間軸(t軸)方向のデータ番号を表すJ以下の整数)を計測する第一の計測手段と、
水抑圧時の、前記2次元k空間もしくは3次元k空間(kx(a),ky(b),kz(c))の前記小型コイル毎の磁気共鳴信号Sm1(L(i))(kx(a),ky(b),kz(c))(t(j))を計測する第二の計測手段と、を備え、
前記演算手段は、
前記Sw1(L(i))(kx(a),ky(b),kz(c))(t(j))に対して、各k空間軸(kx,ky,kz)方向に逆フーリエ演算を施し、2次元実空間もしくは3次元実空間(x(a),y(b),z(c))の磁気共鳴信号Sw2(L(i))(x(a),y(b),z(c))(t(j))を算出する第一のk空間−実空間変換手段と、
前記各小型コイルの実空間点毎に、前記磁気共鳴信号Sw2(L(i))(x(a),y(b),z(c))(t(j))の先頭からN番目の信号Sw2(L(i))(x(a),y(b),z(c))(t(N))の絶対値Aw(L(i))(x(a),y(b),z(c))(t(N))と位相値φw(L(i))(x(a),y(b),z(c))(t(N))とを算出し、両算出結果をそれぞれ補正用絶対値および補正用位相値とする絶対値位相値算出手段と、
前記各小型コイルの実空間点毎に、前記磁気共鳴信号Sw2(L(i))(x(a),y(b),z(c))(t(j))の所定数の点の絶対値信号から、ノイズ平均値Na(L(i))(x(a),y(b),z(c))とノイズ標準偏差Nd(L(i))(x(a),y(b),z(c))とを算出するノイズ算出手段と、
前記各小型コイルの実空間点毎に、前記補正用絶対値と、前記ノイズ平均値Na(L(i))(x(a),y(b),z(c))および前記ノイズ標準偏差Nd(L(i))(x(a),y(b),z(c))とから、信号雑音比SNR(L(i))(x(a),y(b),z(c))を算出する信号雑音比算出手段と、
前記磁気共鳴信号Sm1(L(i))(kx(a),ky(b),kz(c))(t(j))に対して、各k空間軸方向に逆フーリエ変換を施し、前記2次元実空間もしくは3次元実空間(x(a),y(b),z(c))の磁気共鳴信号Sm2(L(i))(x(a),y(b),z(c))(t(j))を算出する第二のk空間−実空間変換手段と、
前記各小型コイルの実空間点毎に、前記磁気共鳴信号Sm2(L(i))(x(a),y(b),z(c))(t(j))の全点の位相値φm(L(i))(x(a),y(b),z(c))(t(j))を前記補正用位相値を用いて位相補正し、位相補正後の磁気共鳴信号Sm3(L(i))(x(a),y(b),z(c))(t(j))を得る位相値補正手段と、
前記各小型コイルの実空間点毎に、前記磁気共鳴信号Sm3(L(i))(x(a),y(b),z(c))(t(j))に前記信号雑音比SNR(L(i))(x(a),y(b),z(c))を用いて算出した加算係数を乗算し、磁気共鳴信号Sm4(L(i))(x(a),y(b),z(c))(t(j))を得る信号雑音比乗算手段と、
前記磁気共鳴信号Sm4(L(i))(x(a),y(b),z(c))(t(j))を、実空間点毎に複素数加算し、信号Sm5(x(a),y(b),z(c))(t(j))を得るMAC合成手段と、
前記磁気共鳴信号Sm5(x(a),y(b),z(c))(t(j))に対して、時間軸方向にフーリエ変換を施し、周波数軸(f軸)方向の情報を有する前記2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像Sm6(x(a),y(b),z(c))(f(j))を算出する時間−周波数変換手段と、を備え、
前記第一の計測手段は、MRSIシーケンスを実行することにより、前記磁気共鳴信号Sw1を得、
前記第二の計測手段は、前記MRSIシーケンスと、前記MRSIシーケンスに先立ち、水信号を抑圧するプリパルスシーケンスとを実行することにより、前記磁気共鳴信号Sm1を得ること
を特徴とする磁気共鳴撮影装置。 - 請求項3記載の磁気共鳴撮影装置であって、
前記ノイズ算出手段は、前記磁気共鳴信号Sw2(L(i))(x(a),y(b),z(c))(t(j))の後部のM点分の絶対値信号を用いて前記ノイズ平均値Na(L(i))(x(a),y(b),z(c))と前記ノイズ標準偏差Nd(L(i))(x(a),y(b),z(c))とを算出すること
を特徴とする磁気共鳴撮影装置。 - 請求項3または4記載の磁気共鳴撮影装置であって、
前記演算手段は、前記磁気共鳴スペクトル画像Sm6(x(a),y(b),z(c))(f(j))に対してスペクトル軸方向の位相補正を施し、前記2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像Sm7(x(a),y(b),z(c))(f(j))を算出する第二の位相値補正手段をさらに備えること
を特徴とする磁気共鳴撮影装置。 - 静磁場、高周波磁場および傾斜磁場をそれぞれ発生する磁場発生手段と、前記静磁場中に置かれた被検体から発生する磁気共鳴信号を検出する検出手段と、前記磁場発生手段および検出手段を制御する計測制御手段と、前記核磁気共鳴信号を用いて磁気共鳴スペクトル画像を作成し、表示する演算手段とを備え、前記検出手段が複数の小型コイルからなる受信コイルを有する磁気共鳴撮影装置であって、
前記受信コイルは複数の小型コイルL(i)(ここで、iはコイル番号を表す1以上I以下の整数:i=1,2,…,I)を備えるマルチアレイコイル(MAC)であり、
前記計測制御手段は、
非水抑圧時の、2次元k空間もしくは3次元k空間(kx(a),ky(b),kz(c))(ここで、aはkx軸方向の計測番号を表すA以下の整数、bはky軸方向の計測番号を表すB以下の整数、cはkz軸方向の計測番号を表すC以下の整数、2次元の場合はA,B,Cのいずれか1つが1で他の2つは2以上の整数、3次元の場合はA,B,Cのいずれもが2以上の整数)の前記小型コイル毎の磁気共鳴信号Sw1(L(i))(kx(a),ky(b),kz(c))(t(j))(ここで、jは時間軸(t軸)方向のデータ番号を表すJ以下の整数)を計測する第一の計測手段と、
水抑圧時の、前記2次元k空間もしくは3次元k空間(kx(a),ky(b),kz(c))の前記小型コイル毎の磁気共鳴信号Sm1(L(i))(kx(a),ky(b),kz(c))(t(j))を計測する第二の計測手段と、を備え、
前記演算手段は、
前記磁気共鳴信号Sw1(L(i))(kx(a),ky(b),kz(c))(t(j))に対して、各k空間軸(kx,ky,kz)方向に逆フーリエ変換を施し、2次元実空間もしくは3次元実空間(x(a),y(b),z(c))の磁気共鳴信号Sw2(L(i))(x(a),y(b),z(c))(t(j))を算出する第一のk空間−実空間変換手段と、
前記各小型コイルの実空間点毎に、前記磁気共鳴信号Sw2(L(i))(x(a),y(b),z(c))(t(j))の先頭からN番目の信号Sw2(L(i))(x(a),y(b),z(c))(t(N))の位相値φw(L(i))(x(a),y(b),z(c))(t(N))を算出し、算出結果から補正値を求める位相値算出手段と、
前記第二の計測手段が計測した前記小型コイル毎のk空間の磁気共鳴信号Sm1(L(i))(kx(a),ky(b),kz(c))(t(j))に対して、各k空間軸方向に逆フーリエ変換を施し、2次元実空間もしくは3次元実空間(x(a),y(b),z(c))の磁気共鳴信号Sm2(L(i))(x(a),y(b),z(c))(t(j))を算出する第二のk空間−実空間変換手段と、
各小型コイルの実空間点毎に、前記磁気共鳴信号Sm2(L(i))(x(a),y(b),z(c))(t(j))の全点の前記位相値φm(L(i))(x(a),y(b),z(c))(t(j))を前記補正値を用いて位相補正し、前記小型コイル毎の位相補正後の磁気共鳴信号Sm3(L(i))(x(a),y(b),z(c))(t(j))を得る位相値補正手段と、
前記磁気共鳴信号Sw1(L(i))(kx(a),ky(b),kz(c))(t(j))を用いて前記小型コイル毎の感度分布を算出する感度分布算出手段と、
前記磁気共鳴信号Sm3(L(i))(x(a),y(b),z(c))(t(j))と、前記感度分布とを用いて、画像の折り返しを除去する演算を行ない、前記小型コイル毎の画像を合成する合成手段と、
前記合成後の画像に対して、時間軸方向にフーリエ変換を施し、周波数軸(f軸)方向の情報を有する2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像Sm5(x(a),y(b),z(c))(f(j))を算出する時間−周波数変換手段と、を備え、
前記第一の計測手段は、MRSIシーケンスを実行することにより、前記磁気共鳴信号Sw1を得、
前記第二の計測手段は、前記MRSIシーケンスと、前記MRSIシーケンスに先立ち、水信号を抑圧するプリパルスシーケンスとを実行することにより、前記磁気共鳴信号Sm1を得ること
を特徴とする磁気共鳴撮像装置。 - 静磁場、高周波磁場および傾斜磁場をそれぞれ発生する磁場発生手段と、前記静磁場中に置かれた被検体から発生する磁気共鳴信号を検出する検出手段と、前記磁場発生手段および検出手段を制御する計測制御手段と、前記核磁気共鳴信号を用いて磁気共鳴スペクトル画像を作成し、表示する演算手段とを備え、前記検出手段が複数の小型コイルからなる受信コイルを有する磁気共鳴撮影装置であって、
前記受信コイルは複数の小型コイルL(i)(ここで、iはコイル番号を表す1以上I以下の整数:i=1,2,…,I)を備えるマルチアレイコイル(MAC)であり、
前記計測制御手段は
非水抑圧時の、2次元k空間もしくは3次元k空間(kx(a),ky(b),kz(c))(ここで、aはkx軸方向の計測番号を表すA以下の整数、bはky軸方向の計測番号を表すB以下の整数、cはkz軸方向の計測番号を表すC以下の整数、2次元の場合はA,B,Cのいずれか1つが1で他の2つは2以上の整数、3次元の場合はA,B,Cのいずれもが2以上の整数)の前記小型コイル毎の磁気共鳴信号Sw1(L(i))(kx(a),ky(b),kz(c))(t(j))(ここで、jは時間軸(t軸)方向のデータ番号を表すJ以下の整数)を計測する第一の計測手段と、
水抑圧時の、前記2次元k空間もしくは3次元k空間(kx(a),ky(b),kz(c))の前記小型コイル毎の磁気共鳴信号Sm1(L(i))(kx(a),ky(b),kz(c))(t(j))を計測する第二の計測手段と、を備え、
前記演算手段は、
前記Sw1(L(i))(kx(a),ky(b),kz(c))(t(j))に対して、各k空間軸(kx,ky,kz)方向に逆フーリエ演算を施し、2次元実空間もしくは3次元実空間(x(a),y(b),z(c))の磁気共鳴信号Sw2(L(i))(x(a),y(b),z(c))(t(j))を算出する第一のk空間−実空間変換手段と、
前記各小型コイルの実空間点毎に、前記磁気共鳴信号Sw2(L(i))(x(a),y(b),z(c))(t(j))の先頭からN番目の信号Sw2(L(i))(x(a),y(b),z(c))(t(N))の絶対値Aw(L(i))(x(a),y(b),z(c))(t(N))と位相値φw(L(i))(x(a),y(b),z(c))(t(N))とを算出し、両算出結果をそれぞれ補正用絶対値および補正用位相値とする絶対値位相値算出手段と、
前記各小型コイルの実空間点毎に、前記磁気共鳴信号Sw2(L(i))(x(a),y(b),z(c))(t(j))の所定数の点の絶対値信号から、ノイズ平均値Na(L(i))(x(a),y(b),z(c))とノイズ標準偏差Nd(L(i))(x(a),y(b),z(c))とを算出するノイズ算出手段と、
前記各小型コイルの実空間点毎に、前記補正用絶対値と、前記ノイズ平均値Na(L(i))(x(a),y(b),z(c))および前記ノイズ標準偏差Nd(L(i))(x(a),y(b),z(c))とから、信号雑音比SNR(L(i))(x(a),y(b),z(c))を算出する信号雑音比算出手段と、
前記磁気共鳴信号Sm1(L(i))(kx(a),ky(b),kz(c))(t(j))に対して、各k空間軸方向に逆フーリエ変換を施し、前記2次元実空間もしくは3次元実空間(x(a),y(b),z(c))の磁気共鳴信号Sm2(L(i))(x(a),y(b),z(c))(t(j))を算出する第二のk空間−実空間変換手段と、
前記各小型コイルの実空間点毎に、前記磁気共鳴信号Sm2(L(i))(x(a),y(b),z(c))(t(j))の全点の位相値φm(L(i))(x(a),y(b),z(c))(t(j))を前記補正用位相値を用いて位相補正し、位相補正後の磁気共鳴信号Sm3(L(i))(x(a),y(b),z(c))(t(j))を得る位相値補正手段と、
前記各小型コイルの実空間点毎に、前記磁気共鳴信号Sm3(L(i))(x(a),y(b),z(c))(t(j))に前記信号雑音比SNR(L(i))(x(a),y(b),z(c))を用いて算出した加算係数を乗算し、磁気共鳴信号Sm4(L(i))(x(a),y(b),z(c))(t(j))を得る信号雑音比乗算手段と、
前記磁気共鳴信号Sw1(L(i))(kx(a),ky(b),kz(c))(t(j))を用いて前記小型コイル毎の感度分布を算出する感度分布算出手段と、
前記磁気共鳴信号Sm3(L(i))(x(a),y(b),z(c))(t(j))と、前記感度分布とを用いて、画像の折り返しを除去する演算を行ない、前記小型コイル毎の画像を合成する合成手段と、
前記合成後の画像に対して、時間軸方向にフーリエ変換を施し、周波数軸(f軸)方向の情報を有する2次元実空間もしくは3次元実空間の磁気共鳴スペクトル画像Sm5(x(a),y(b),z(c))(f(j))を算出する時間−周波数変換手段と、を備え、
前記第一の計測手段は、MRSIシーケンスを実行することにより、前記磁気共鳴信号Sw1を得、
前記第二の計測手段は、前記MRSIシーケンスと、前記MRSIシーケンスに先立ち、水信号を抑圧するプリパルスシーケンスとを実行することにより、前記磁気共鳴信号Sm1を得ること
を特徴とする磁気共鳴撮影装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011222934A JP5479427B2 (ja) | 2007-01-17 | 2011-10-07 | 磁気共鳴撮影装置 |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007008008 | 2007-01-17 | ||
JP2007008008 | 2007-01-17 | ||
JP2011222934A JP5479427B2 (ja) | 2007-01-17 | 2011-10-07 | 磁気共鳴撮影装置 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2008553988A Division JP4845972B2 (ja) | 2007-01-17 | 2007-12-17 | 磁気共鳴スペクトル画像算出方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2012035097A JP2012035097A (ja) | 2012-02-23 |
JP5479427B2 true JP5479427B2 (ja) | 2014-04-23 |
Family
ID=39635827
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2008553988A Active JP4845972B2 (ja) | 2007-01-17 | 2007-12-17 | 磁気共鳴スペクトル画像算出方法 |
JP2011222934A Active JP5479427B2 (ja) | 2007-01-17 | 2011-10-07 | 磁気共鳴撮影装置 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2008553988A Active JP4845972B2 (ja) | 2007-01-17 | 2007-12-17 | 磁気共鳴スペクトル画像算出方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US8258785B2 (ja) |
JP (2) | JP4845972B2 (ja) |
WO (1) | WO2008087822A1 (ja) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4845972B2 (ja) * | 2007-01-17 | 2011-12-28 | 株式会社日立メディコ | 磁気共鳴スペクトル画像算出方法 |
JP2010100474A (ja) * | 2008-10-23 | 2010-05-06 | Covalent Materials Corp | シリコン単結晶引上げ水平磁場の最適化方法およびシリコン単結晶の製造方法 |
CN101963657B (zh) * | 2009-07-24 | 2012-12-05 | 西门子(深圳)磁共振有限公司 | 边带抑制方法和边带抑制装置 |
US9389289B2 (en) | 2010-10-28 | 2016-07-12 | Hitachi Medical Corporation | Magnetic resonance imaging device |
WO2013027539A1 (ja) * | 2011-08-24 | 2013-02-28 | 株式会社日立メディコ | 磁気共鳴撮影装置、位相値補正方法およびプログラム |
US9766317B2 (en) * | 2011-10-25 | 2017-09-19 | New York University | Apparatus, systems and methods which are based on magnetic resonance imaging for evaluation(s) of radio frequency emitting device(s) |
CN104220873B (zh) * | 2012-02-10 | 2017-12-22 | 儿童医疗中心公司 | 基于nmr的代谢物筛选平台 |
JP6084420B2 (ja) * | 2012-10-02 | 2017-02-22 | 株式会社日立製作所 | 多チャンネル画像合成方法及び磁気共鳴イメージング装置 |
US9008399B2 (en) * | 2013-01-28 | 2015-04-14 | Wisconsin Alumni Research Foundation | System and method for joint reconstruction of spatial and spectral imaging data |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5709208A (en) * | 1994-04-08 | 1998-01-20 | The United States Of America As Represented By The Department Of Health And Human Services | Method and system for multidimensional localization and for rapid magnetic resonance spectroscopic imaging |
JP3598411B2 (ja) * | 1997-09-09 | 2004-12-08 | 株式会社日立メディコ | 磁気共鳴診断装置 |
JP3427059B2 (ja) | 2001-04-13 | 2003-07-14 | 株式会社東芝 | 磁気共鳴診断装置 |
JP2002315732A (ja) | 2001-04-23 | 2002-10-29 | Toshiba Corp | 磁気共鳴装置及び磁気共鳴装置における信号処理方法 |
DE10243830B4 (de) * | 2002-09-13 | 2006-11-16 | Universität Bremen | Spektroskopisches Bildgebungsverfahren sowie Verwendung desselben zur Materialcharakterisierung |
US6995560B2 (en) * | 2003-03-20 | 2006-02-07 | Duerk Jeffrey L | Chemical species suppression for MRI imaging using spiral trajectories with off-resonance correction |
CN100396239C (zh) * | 2005-02-28 | 2008-06-25 | 西门子(中国)有限公司 | 磁共振多通道成像水脂分离重建算法 |
US7298144B2 (en) * | 2005-05-06 | 2007-11-20 | The Board Of Trustee Of The Leland Stanford Junior University | Homodyne reconstruction of water and fat images based on iterative decomposition of MRI signals |
JP4845972B2 (ja) * | 2007-01-17 | 2011-12-28 | 株式会社日立メディコ | 磁気共鳴スペクトル画像算出方法 |
-
2007
- 2007-12-17 JP JP2008553988A patent/JP4845972B2/ja active Active
- 2007-12-17 US US12/523,259 patent/US8258785B2/en active Active
- 2007-12-17 WO PCT/JP2007/074205 patent/WO2008087822A1/ja active Application Filing
-
2011
- 2011-10-07 JP JP2011222934A patent/JP5479427B2/ja active Active
Also Published As
Publication number | Publication date |
---|---|
US8258785B2 (en) | 2012-09-04 |
US20100013481A1 (en) | 2010-01-21 |
JP4845972B2 (ja) | 2011-12-28 |
JPWO2008087822A1 (ja) | 2010-05-06 |
JP2012035097A (ja) | 2012-02-23 |
WO2008087822A1 (ja) | 2008-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5479427B2 (ja) | 磁気共鳴撮影装置 | |
US9766313B2 (en) | MR imaging using apt contrast enhancement and sampling at multiple echo times | |
JP5189203B2 (ja) | 磁気共鳴装置 | |
CN105308469B (zh) | 一种mr成像方法、mr设备以及相关数据载体 | |
US9223001B2 (en) | MR imaging using navigators | |
JP4889791B2 (ja) | 磁気共鳴撮影装置 | |
US9709641B2 (en) | Magnetic resonance imaging apparatus, image processing apparatus, and susceptibility map calculation method | |
JP4610611B2 (ja) | 磁気共鳴撮影装置 | |
JP6038654B2 (ja) | 磁気共鳴イメージング装置及び振動誤差磁場低減方法 | |
JP4564015B2 (ja) | 磁気共鳴撮影装置及び磁気共鳴撮影方法 | |
JP4781120B2 (ja) | 磁気共鳴撮影装置および磁気共鳴スペクトル計測方法 | |
JP5214209B2 (ja) | 磁気共鳴イメージング装置 | |
JP5536358B2 (ja) | 磁気共鳴イメージング装置及び感度補正方法 | |
JP4980662B2 (ja) | 磁気共鳴イメージング装置 | |
US11965945B2 (en) | Magnetic resonance system and shimming method and imaging method thereof | |
JP5156958B2 (ja) | 磁気共鳴イメージング装置 | |
JP5636058B2 (ja) | 磁気共鳴撮影装置 | |
US20150253402A1 (en) | Metabolite detection system and operation thereof | |
JP4832510B2 (ja) | 磁気共鳴撮影装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A132 Effective date: 20130226 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20130418 |
|
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: 20140114 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20140212 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5479427 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
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 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |