JPS6058539A - Νmr断層像撮影装置 - Google Patents

Νmr断層像撮影装置

Info

Publication number
JPS6058539A
JPS6058539A JP58168040A JP16804083A JPS6058539A JP S6058539 A JPS6058539 A JP S6058539A JP 58168040 A JP58168040 A JP 58168040A JP 16804083 A JP16804083 A JP 16804083A JP S6058539 A JPS6058539 A JP S6058539A
Authority
JP
Japan
Prior art keywords
magnetic field
gradient
field gradient
nxn
signals
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
Application number
JP58168040A
Other languages
English (en)
Other versions
JPH0421492B2 (ja
Inventor
Akinori Fujita
明徳 藤田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shimadzu Corp
Shimazu Seisakusho KK
Original Assignee
Shimadzu Corp
Shimazu Seisakusho KK
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shimadzu Corp, Shimazu Seisakusho KK filed Critical Shimadzu Corp
Priority to JP58168040A priority Critical patent/JPS6058539A/ja
Publication of JPS6058539A publication Critical patent/JPS6058539A/ja
Publication of JPH0421492B2 publication Critical patent/JPH0421492B2/ja
Granted legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 (イ)産業上の利用分野 この発明はNMR断層像撮影装置に係り、qlに、同時
に採取される2組のN′個のサンプリング点よりなるN
行の核磁気共鳴伝号からフーリエ変換によりNXN“個
の画素を持つ萌IM像をi!する装置に関する。
(ロ)従来技術 第1図において体軸方向をZ方向、れ(方向をY方向、
措方向をX方向とした場合、jiL末のN M R断層
像撮影装置では、前記X方向Gこ磁界グラジェントGx
を加えるごとにより、被検体に1−)、え人口 ″定周
波数の高周波を周波数変調し、一方、i’+ii記Y方
向に磁界グラジェントGyを加えろことで前記高周波を
位相変調することにもとすいて、コンピュータ計算を行
って被検体の1IJi層像を4iJ影゛4゛イ)十法(
Kuma r、 We l L、Ijr n s LO
) ’、次ノcフーリエ変換法、あるいはこの改良法で
J〕るスピンワープ法等で、以後二次元フーリエ変換法
という。)がある(参考文献: Z、 11.cIIO
他著 ” I’ourier Transform N
uclear Magnetic Resonance
Tomographic Imaging PIIOC
EEDINGS OF THEIEEE 、VOL、7
0 、NO,Lo、1152 〜1173,0CTOr
lER1982)。
しかして、前記二次元フーリエ変換により断層像を計算
する場合、一つの断層面内のNXN個の画素像を得るた
めには各行が2N個のサンプリング点をもつ2N行の核
磁気共鳴信号が必要であった。
しかして、1行の核磁気共鳴信号を採取するに要する時
間としては核磁気モーメントの縦緩和時間が約1秒であ
ることからそれと略等しく約1秒が必要である。
したがって、必要な個数の核磁気共1!!X信号を(4
4,’るための信号採取時間が長くかかるという欠点が
あった。更に、前記断層像を得るための計q貼間もそれ
に応して長くなるため、信号を採取しはじめてから1f
7N像を写し出すまでの時間、ずなわぢ前記信号採取時
間と前記i’ll算時間の合計11.1間も長くなって
しまうという欠点があった。
かかる欠点を排除するため、本発明借は、Y方向に対し
ては所定の関係で定まるθ+に+’タラシ1ントGyを
加えることによりそれぞれ2N’個のサンプリング点よ
りなるN行の核ωヶ気共鳴(−号をもとめ、このN行の
核磁気共鳴信号をコンピー・、−タで二次元フーリエ変
換することによりN X N ’の画素からなる断面像
を得る手段を提案している。
(特願昭58−056880号)。
しかし、これによると、採取イ1jυ:、)か牛骨に4
1′るため断層像のS/Nは、(it来の1/仮になる
いう欠点をかある。
(ハ)目的 この発明は従来のNMR断bri IS:撮影装置に比
・咬して、S/Nを低下させることなく 、+ii+記
信゛・」採取時間及びコンピュータによる計算時間をそ
の半分にすることのできるN M R断1<・i +象
4:I!r ’4’+ユ装置を1是供することを目的と
している。
(ニ)構成 この発明はY方向(縦方向)にグ・1しては、通、・i
への2倍の磁界グラジエン1−Gyを−fノえ、かつ、
11いに90′位相の異なる二つのレファレンス(3号
を用いることにより、各行それぞれN“個のサンプリン
グ点よりなるN行の核磁気共鳴信号を2組得て、従来装
置と同じS/Nで、同じ分解能を持つNXN“の断面像
を得るように構成されている。
(ホ)実施例 前記第1図は本発明の説明にも用いられるもので、同図
(イ)において1は被検体、2は撮;;;ニずべき断面
をしめず。また、同図(ロ)は前記断面2そのものを取
り出したところを示しCいる。
ここで、核磁気共鳴信号を得る手lll、Hについて説
明する。
即ち、体軸(Z軸)方向の断面位置は90”RFパルス
(例えば6Mz)の印加と同時に、磁界グラジェントG
zを印加することにより決定される。
つき′に、X方1句について磁界グラジェントGXを与
えることにより周波数変調をおこなう(但し、X=Oと
χ−Lxの位置では共鳴周波数差がfとなるように磁界
グラジェントGx4よ定められる)。
しかして、Y方向の磁界グランフーンl−G yとして
は(1)式に示すように、各行(L=0、■、2、・・
・・N−1)ごとに加える。
Gy (L) −GOy(L−N/2) ・・・ (1
)ただし、ayは第2図にしめずように、中心Y=Ly
/2で零となるY方向の(j陽1’グラジエンl−であ
って、 y=oで、γGOy=−11 Y = ]、 yで、γQOy=IIであ、乙。
ここで、Tは磁気回転比である。ずなわし、上記したこ
とはY方向の磁界グラジェントを(ヌf末の倍にしたこ
とをは味している。
前記条件のもとで、発生した核(a?気共鳴信′・シを
X = L x / 2におLJる共鳴周波数を持し、
かり、90′位相の異なる二組のレソアランス(J3号
で後述の第6図のブロック図で示ずよ]Gこミキンンク
したあと、ローパスフィルタを通”」ことにより、0〜
f/2(上記例ではθ〜5KllZ)の信号に変換し、
さらにこれらをサンプリング定理を活たず1/fの時間
m1隔でN点ずつサンプリングする。
ごのような手1IurtをI−= 0、■、2、・・・
N−1として繰り返ずことでN(l!lIのザンプリン
グ点よりなるN行の信号か二月1得られる。
つぎに、上記のようにしてj5Jられた二組のNxNの
信号の一方を実部として、他方を虚部として、NxN複
素行列データをつくり、これを二次元フーリエ変換して
、その絶対値を計算することにより、NxN個の画素よ
りなる断面像を得る。
上記した原理を数式を用いて、以下説明する。
第3図は被検体の断面のうちで撮影すべき対象物(同図
では三角形)21を表している。被検体に加えたZ軸方
向の静磁、場をトI。とする。そして、X = L x
 / 2の位置での共鳴周波数(以下、中心共鳴周波数
という)をf′とすると、 M=211ω’ = 2 n r I−1(、・・・・
・ (2)さて、第3図の図示するようなρ<x、y 
)の密度を持つ対象に対し、前述のシーケンスを通用す
ると、第■3回目の繰り返しO4に冑られるべき核磁気
共鳴信号s <v、+、 >ば、 y となる。
コレラ・x = L x / 2の共鳴周波数を1.5
つレフアランス信号rl(L) rl (L)=CO3(2月Y′110)−・141で
ミキシングすると、 Ml (t、I−)=3 (t、L) ・r I CL
)211−−(1、−N/ 2) 4 UJ I31 して(5)i:Cの第1JAを除番ノば、(6)式の信
りS’l (t、L)を得る。すなわち、 x 1・y ・ ・ ・ ・ (6) これをサンプリング定理を病たず間1へλ、4、Δt=
1/fでN点ザンプリングすると、t=にIf (k=
0.1.2. ・−・・If−1)より、(7)式が成
立する。ずなゎjノ、+、y ・ ・ ・ ・ ・ (7) ここで、k=Q、1.2.・・・、l1lL=0. 1
. 2. ・・・、N−1である。
さらに、第3図のように、× (0〜1.x)および、
y (0〜Ly)をN分;(;す”4ろごとにより、対
象を1ζ)を軟化すれば、 x=n+Lx /N、、m−0,1,・・ N 1y=
nLy /N、、n=0.1. ・・・N−1・・ (
8) を用いて(9)式のようにあられされろ。
但し、 k=o、1.2. ・・・、N−1 L=0.1,2. ・・・、N−1 ・・・ (10) なお、一般的にはNは128または256程度の数をも
ちいるのが普通である。
さて、もう一つのレファレンス信号r2(L)として、
rl(L)と同じ周波数を持し、位相が90°遅れた信
号を用いる。
すなわち、 r2 (t)=CO3(2rlf ’ を十〇−II 
/ 2 )=S IN (21’l r ’−+−0)
 ・・・ (11)(3)式で表される核磁気共鳴信号
S(t、L)に対して、r2(+、)でミキソンクをi
Jい1.!−バスフィルタを通し、さらに△L−1/ 
l”ir N点→」ンプリングを行うと、+iii述の
、ll算さ同L’、’ 、I:、−)にして、 (12
)j:Cで表されるS2 (k、l、)を111ン〕。
・・・ (12) を作る。
s (k、L) −31(k、L) x j +S2 
(k、L) ・・・ (13) すると、s (k、、1.、)は(14)弐〇表されろ
・ ・ ・ (14) 但し、k=o、i ・・ ・、N−1 L=O,i ・ ・ ・、N−1 ・ ・ ・ (15) 以上のs (k 、1.)を得る手続きを第4図に示す
さて、離散フーリエ変換の周期性より、対象ρ(rrl
、n)は、05m<N、0≦n < N以外でば、m、
nに関し、それぞれNの周期を持つとしてよい。
ここで、同じ(、rn、nに関して周期Nを持つ関数ρ
′ (rn、n)を(■6)式で定義する。
・ ・ ・ (l 6) すると、s (k、L)は1次のように変形−(きる。
(L−−) −0)j) Exp (−11nj)Exp (−jO))X=N2
DI;l’、N Ilt ’ (m、n) xExp 
(−11n j) F、X[l (j D) 1いての
N点の811敗二次元フーリコー変換を表し、1)it
’s (k (k、L) )は、mおよび[1につい′
(のN点の離散二次元逆フーリエ変換をAI、ずものと
yる。
(17)式より、次式を得る。
DN、N fs (k、L)) =N2p ’ (m、
n)xExp (−nnj)EXI)(jθ)・・・ 
(18) (■8)式の絶対値をとれば、次式をf4jる。
IDNμfs (k、L) l l ・・・ (19) (19)式で得られる像を第5図に示す。
m、nは周期Nをもつもので、 m’=rn+ − ・・・ (20) n””rl+□ として、第5図を変換すれば、第3図の密度分布を持つ
対象の像が1縁られる。
第6図は上述したN M RyB漸撮影装置を実現して
、2組のN行N列の核磁気共鳴信号を取り出してこれを
コンピュータに取り込むための回路ブロック図を示して
いる。
ずなわち、90’RFパルス発仕回路20が基準クロッ
ク発生回路10により頼ス動され、90゜RFパルスを
発生し、RF増幅器30を介し一〇、送受切り替え回路
40に入り、送受切り替え回路40ばそれに連結された
R F二フィル50に00”RFパルスをおくる。
磁界グラジェントGzと前記(1)式でL = +)と
したときのGy(L)が加えられ、さらにGyが加えら
れることに基すいて、核磁気共鳴信号か1行求められ、
前記送受切り替え回路40 、+’+i+ii・W゛増
幅器60を介し−i1’ 2 iJIのミキサー回路7
0に人力する。
一方、前記基準クロック発η゛回路10からのノ、l、
準クロックはレファランス信号回路80に入り、ここで
互いに90°位相の異なる2つのシファフンス信号が発
生ずる。この2つの(11号は、にの1)の前記ミキサ
ー回路70に人力しミキサー回路70からの出力は1寸
個のサンプリング点を含むアナログ信号としてローパス
フィルタ()0に入力し、A/D回路100の出力は第
1行のデータとしてコンピュータ110に人力する。
上記と同様の操作を行うことにより、2組のし−】、2
・・・1″−J−1行のデータがコンピュータ110に
入力する。
コンピュータ110は2組のN個のサンプリング点から
なるN行のデータに基づいて、」−述の二次光フーリエ
変換を行い、第3図で図示したρ(x+y)の分布密度
をもった所望の映像を別設の再生装置に現出する。
また、上述の実施例では、得るべき画素としてNXHの
ように正方形の映像について、説明したが、実際の被検
体の検査するべき映像ば長方形のものもある。したがっ
て、−走査中に含まれるデータの個数は上記Nに限定さ
れず、一般的にはN′とし“ζ取り扱うことができる。
(へ)効果 この発明に係るN M R断層像撮影装置はnIYツ面
のX方向を周波数変調し、Y方向を位相変調して核磁気
共鳴信号を得る装置であって、Y方向に対してば所定の
関係で定まる磁界クラジエンl−G yを加え、また、
H−、いに90°位相の5’4!、 ’/、;・に)2
つのレファランス信号でミキシングする、二とよって、
2組のそれぞれN 117.1のサンプリング点よりな
るN行の核磁気共鳴信号を同時に取り込み、この(Ii
号に対し、2次元フーリエ変換を1」っでいるから、(
iL来のNMR断層像撮影装jjI4:に比・\、S/
Nを低下させることなく、信号採取時間を1−うjに’
p、j紬゛ζきる。
【図面の簡単な説明】
第1図は体軸方向をZ方向、糸υ(方向をY方向、横方
向をX方向とした場合の被検体とそのlすt)τ・)像
の位置関係を表す説明図、第21′;!目41本光明ご
あるNMR断層(象1最影装置におけるQ yクシ−/
上ン1のY方向分布をあられず説明図、第;)図は破4
1ミ体の断面のうらで撮影゛yへきり・1゛′ρ(1,
、l lス1“(はづ1]形)を表した説明図、第4図
は採取−1イ)核磁気共鳴信号の処理手順を説明するだ
めの図、第5図はI):j記(18)式の値を絶対値で
11、小し71図、第(3図は上述したNMR断層像撮
影装置を実現3′ろために2組のN行の核磁気共鳴信号
を取り出し、これをコンピュータに取り込むための回路
のフ゛口・ノクレ1をそれぞれ示している。 特許出11’I人 株式会社 島津製作所代理人 弁理
士 大 西 孝 治

Claims (1)

    【特許請求の範囲】
  1. (1)被検体のX方向に磁界グラジエン)Gxを加える
    ことにより、被検体に与えた一定周波数の高周波を周波
    数変調し、一方、前記Y方向に磁界グラジエンI−G 
    yを加えることにより前記高周波を位相変調することに
    基づいて被検体の断層面を得るNMR断層像撮影装置に
    おいて、Y方向に対しては所定の関係で定まる磁界グラ
    ジェントcyを加え、かつ、互いに90°位相の異なる
    2つのレファランス信号を用いることによって、それぞ
    れN“個のサンプリング点よりなるN行の核磁気共鳴信
    号を2組同時に採り込み、この−組を実部、信組を虚部
    とたNXN“の複素行列データをつくり、これをコンピ
    ュータで二次元フーリエ変換を行うことによりNXN“
    の画素からなる断面像を得ることを特徴とするNMR断
    fmり撮影装置。
JP58168040A 1983-09-12 1983-09-12 Νmr断層像撮影装置 Granted JPS6058539A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP58168040A JPS6058539A (ja) 1983-09-12 1983-09-12 Νmr断層像撮影装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP58168040A JPS6058539A (ja) 1983-09-12 1983-09-12 Νmr断層像撮影装置

Publications (2)

Publication Number Publication Date
JPS6058539A true JPS6058539A (ja) 1985-04-04
JPH0421492B2 JPH0421492B2 (ja) 1992-04-10

Family

ID=15860698

Family Applications (1)

Application Number Title Priority Date Filing Date
JP58168040A Granted JPS6058539A (ja) 1983-09-12 1983-09-12 Νmr断層像撮影装置

Country Status (1)

Country Link
JP (1) JPS6058539A (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6462142A (en) * 1987-09-03 1989-03-08 Toshiba Corp Magnetic resonance imaging method
JPH01146529A (ja) * 1987-12-02 1989-06-08 Toshiba Corp 磁気共鳴イメージング装置

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7588486B2 (ja) 2020-09-16 2024-11-22 株式会社きもと 農作物販売システム及び区画選択受付方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS57147042A (en) * 1981-01-26 1982-09-10 Nat Res Dev Inspection of sample

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS57147042A (en) * 1981-01-26 1982-09-10 Nat Res Dev Inspection of sample

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6462142A (en) * 1987-09-03 1989-03-08 Toshiba Corp Magnetic resonance imaging method
JPH01146529A (ja) * 1987-12-02 1989-06-08 Toshiba Corp 磁気共鳴イメージング装置

Also Published As

Publication number Publication date
JPH0421492B2 (ja) 1992-04-10

Similar Documents

Publication Publication Date Title
Holsinger et al. The importance of phase‐encoding order in ultra‐short TR snapshot MR imaging
Wigström et al. Temporally resolved 3D phase‐contrast imaging
FI87953C (fi) Foerfarande och anordning foer framstaellning av nmr-bilder, fria fraon aliasingfel
US6515476B1 (en) Magnetic field inhomogeneity measurement method and apparatus, phase correction method and apparatus, and magnetic resonance imaging apparatus
Barends et al. Logic gates at the surface code threshold: Superconducting qubits poised for fault-tolerant quantum computing
JPS6058539A (ja) Νmr断層像撮影装置
US4818941A (en) MR tomography method and apparatus for performing the method
Lee et al. Accelerated 3D myelin water imaging using joint spatio‐temporal reconstruction
JP2598038B2 (ja) 磁気共鳴映像装置
US5363043A (en) Producing dynamic images from motion ghosts
JPS6289176A (ja) 核磁気共鳴画像の補間拡大方法
KR850000679A (ko) 핵자기 공명을 영상화하는 방법
CN114325528B (zh) 一种磁共振成像方法及相关设备
JPS59182351A (ja) Nmr断層像撮影装置
JP3033851B2 (ja) 磁気共鳴イメージング装置
Rybicki et al. Reconstruction algorithm for novel ultrafast magnetic resonance imaging
Levin et al. A new class of sampling theorems for Fourier imaging of multiple regions
JPH01121041A (ja) 磁気共鳴映像法
Wang et al. A fast wavelet multigrid algorithm for solution of electromagnetic integral equations
Itami et al. A realization method of cyclic filter banks
JPH04183435A (ja) 磁気共鳴イメージング装置
JPH04150833A (ja) 核磁気共鳴イメージング装置
JPH04138134A (ja) 核磁気共鳴イメージング装置
EP0563651A1 (en) Method for improving the resolution of solid-state NMR multiple-pulse imaging systems using second averaging
JPH0271727A (ja) 核磁気共鳴イメージング装置