JPH02224737A - Mr imaging method - Google Patents
Mr imaging methodInfo
- Publication number
- JPH02224737A JPH02224737A JP1049377A JP4937789A JPH02224737A JP H02224737 A JPH02224737 A JP H02224737A JP 1049377 A JP1049377 A JP 1049377A JP 4937789 A JP4937789 A JP 4937789A JP H02224737 A JPH02224737 A JP H02224737A
- Authority
- JP
- Japan
- Prior art keywords
- magnetic field
- phase
- matphase
- matrix
- inclined magnetic
- 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
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
Description
この発明は、核磁気共鳴(NMR)を利用してイメージ
ングを行なうMR撮像法に関する。The present invention relates to an MR imaging method that performs imaging using nuclear magnetic resonance (NMR).
通常のMR撮像法では、NXNのマトリクスの画像を得
ようとする場合、位相エンコード方向のエンコード数は
Nとし、N回の撮像シーケンスを繰り返してデータの収
集を行なう(これをフルエンコード法と呼ぶ)。
位相エンコード数を減少させることは撮像シーケンスの
繰り返し回数を減少させることになり、全体の撮像時間
を短縮することになる。
そこで、従来より、上記のフルエンコード法での位相エ
ンコード周辺領域をノイズレベルと見なしてデータ収集
を行なわず、そのデータの欠けた領域にゼロデータを補
充して正方生データマトリクスを生成し、これを2次元
フーリエ変換して画像再構成することが行なわれている
(これをエンコードリダクション法と呼ぶ)。In the normal MR imaging method, when trying to obtain an NXN matrix image, the number of encodes in the phase encoding direction is set to N, and data is collected by repeating the imaging sequence N times (this is called the full encoding method). ). Reducing the number of phase encodes will reduce the number of repetitions of the imaging sequence, which will shorten the overall imaging time. Therefore, in the past, the area around the phase encoding in the full encoding method described above is regarded as a noise level and data is not collected, and the area where the data is missing is filled with zero data to generate a square raw data matrix. Image reconstruction is performed by performing a two-dimensional Fourier transform on the image (this is called the encode reduction method).
しかしながら、エンコードリダクション法では、たしか
に全体の撮像時間は短縮されるものの、イメージング視
野を正方視野に保ったままデータ収集しているため、再
構成された画素が非等方的になり、再構成画像の空間分
解能が劣化する問題が生じる。
この発明は、フルエンコード法とまったく同一の空間分
解能を得ることができ、しがち全体の撮像時間をエンコ
ードリダクション法と同様に短縮できる、MR撮像法を
提供することを目的とする。However, although the encode reduction method certainly shortens the overall imaging time, since data is collected while keeping the imaging field of view square, the reconstructed pixels become anisotropic and the reconstructed image The problem arises that the spatial resolution of the image is degraded. An object of the present invention is to provide an MR imaging method that can obtain exactly the same spatial resolution as the full encoding method and shorten the overall imaging time as well as the encode reduction method.
上記目的を達成するため、この発明によるMR撮像法に
おいては、被写体の形状に応じて位相エンコード方向の
サイズを周波数エンコード方向のサイズよりも小さく定
めた非正方イメージング視野に対して、画素が等方的に
なるように位相エンコード方向のマトリクスと周波数エ
ンコード方向のマトリクスとを決め、位相エンコード数
を位相エンコード方向マトリクスと同数とするとともに
、上記両方向のマトリクスの比により各回の位相エンコ
ード量を定めるようにしている。In order to achieve the above object, in the MR imaging method according to the present invention, pixels are arranged isotropically with respect to a non-square imaging field in which the size in the phase encoding direction is determined to be smaller than the size in the frequency encoding direction according to the shape of the subject. The matrix in the phase encoding direction and the matrix in the frequency encoding direction are determined so that the matrix in the phase encoding direction is the same as the matrix in the phase encoding direction, and the amount of phase encoding each time is determined by the ratio of the matrices in both directions. ing.
イメージング視野を、その位相エンコード方向のサイズ
を周波数エンコード方向のサイズよりも小さく定めるこ
とにより、非正方とする。そしてこの非正方視野に対し
て、画素が等方的になるように位相エンコード方向のマ
トリクスと周波数エンコード方向のマトリクスとを決め
る。
この位相エンコード方向のマトリクスに対応した数の位
相エンコードを行なう、これにより、位相エンコード方
向のサイズを周波数エンコード方向のサイズよりも小さ
くした分だけ位相エンコード数を少なくでき、撮像時間
の短縮を図れる。さらに、その各回の位相エンコード量
は位相エンコード方向のマトリクスと周波数エンコード
方向のマトリクスとの比により定めるので、画素自体の
等方性は保たれ、再構成画像の空間分解能が劣化するこ
とを避けることができる。
すなわち、イメージング視野を位相エンコード方向に狭
くすることにより、画像の分解能の劣化を招かずに撮像
時間の短縮を達成することができる。The imaging field of view is made non-square by setting its size in the phase encoding direction to be smaller than its size in the frequency encoding direction. Then, for this non-square visual field, a matrix in the phase encoding direction and a matrix in the frequency encoding direction are determined so that the pixels are isotropic. The number of phase encodes corresponding to the matrix in the phase encode direction is performed, whereby the number of phase encodes can be reduced by the amount that the size in the phase encode direction is smaller than the size in the frequency encode direction, and the imaging time can be shortened. Furthermore, since the amount of phase encoding each time is determined by the ratio of the matrix in the phase encoding direction and the matrix in the frequency encoding direction, the isotropy of the pixel itself is maintained and the spatial resolution of the reconstructed image is prevented from deteriorating. Can be done. That is, by narrowing the imaging field of view in the phase encoding direction, it is possible to shorten the imaging time without causing deterioration in image resolution.
つぎにこの発明の一実施例について図面を参照しながら
説明する。第1図はこの発明にががるMR撮像法をフィ
ールドエコー法に適用した一実施例のパルスシーケンス
を示すもので、これはたとえば第3図のようなシステム
により行なわれる。
すなわちRF送受信装置34よりRFコイル35に第1
図Aに示すようなエンベロープの高周波信号を送り、撮
像対象たる被写体(図示しない)のスピンを励起する。
このとき測定制御装置31の制御のもとに傾斜磁場電源
32から第1図Bに示すような波形の電流が傾斜磁場発
生用コ、イル33に送られ、スライス選択用の傾斜磁場
(ここではGzとする)が印加される。その後、傾斜磁
場電源32から第1図C,Dのような波形の電流が傾斜
磁場発生用コイル33に送られて、周波数エンコード用
傾斜磁場(ここではGxとする)と位相エンコード用傾
斜磁場(ここではGyとする)が印加される。これらの
波形及びタイミングは測定制御装置31によって定めら
れる。被写体から発生したNMR信号はR,Fコイル3
5により受信され、RF送受信装置34を経てコンピュ
ータ36に取り込まれてデータの収集が行なわれる。
撮像対象たる被写体は通常人体であり、一般には第2図
に示すようにその断面像21の横・縦のサイズは必ずし
も同一でない部位が多く、したがって横のサイズLxを
縦のサイズLyよりも大きくした非正方のイメージング
視野22の方がよい場合が多い、そのため、イメージン
グ視野22をLx=Lyとすることは重要でなく、むし
ろ画素自体の等方性の方が重要である。したがって、こ
の非正方のイメージング視野22に対して画素の等方性
を損なわないように
L x / MAT freq= L y / MAT
phaseの関係式を満たす、横方向(周波数エンコー
ド方向)マトリクスMAT freqと縦方向(位相エ
ンコード方向)マトリクスMATρhaseとを定める
。
そして位相エンコード数はMATphaseと同数とし
て、第1図で示したパルスシーケンスの繰り返し回数を
MAT freq−MATphaseだけ減少させる。
コノことにより全体の撮像時間の短縮が図れる。
ここで、上記の任意マトリクスを規定するマトリクスウ
ェイト(VV+++at)として、W mat=MAT
phase/MAT freqで定義されたパラメータ
を導入する。第1図りの斜線部で示す毎回の位相エンコ
ード量epは次のようにして定める。ナイキスト周波数
にてサンプリングするのと等価な位相制御を要すること
から、位相エンコード波形Gyは
γLyiGydt=π
(γ;磁気回転比)
のrIA係を満たす必要がある。任意マトリクス時の撮
像においては、W matパラメータを含むよう上記の
3つの式より次の式を導く。
SGy dt= (π/γLx) (1/W mat
)この式の右辺のπ/γLxは通常のフルエンコード法
におけるエンコード量を示すから、任意マトリクス時の
第1図の1回の位相エンコード量Gpは、フルエンコー
ド法における位相エンコード量を1 / W mat倍
したものとすればよいことが分かる。このW matの
範囲は、
Q(Wmat≦1
であり、W mat=1はフルエンコード法による撮像
となる。
断面像21に対応してイメージング視野22のサイズL
x、Lyを定め、MAT freq、 MATphas
eを決めると、コンピュータ36においてW matの
計算がなされ、エンコード量Gpが求められる。このG
pは測定制御装置31に送られ、これに応じて位相エン
コード用傾斜磁場Gyの1回のエンコード量が制御され
る。こうして得られた生データはコンピュータ36で2
次元フーリエ変換されて第2図で示すような非正方のイ
メージング視野22での画像再構成がなされる。この再
構成画像では画素の等方性は保たれ、画像の空間分解能
は劣化しない。
たとえば任意マトリクス70%(W mat= 0 。
7)としたとき、フルエンコード法による撮像に比べて
撮像時間が70%(MATphase= 0 、7 X
MAT freqだから)、位相エンコード方向の視
野サイズは周波数エンコード方向の視野サイズの70%
(Ly=0.7Lx)となる。
なお、上記ではこの発明をフィールドエコー法に適用し
た一実施例について説明したが、位相エンコードを行な
うものであればスピンエコー法などの他の撮像法にも適
用できることは勿論である。
また、2次元イメージングのみならず3次元イメージン
グへも容易に適用できる。さらに3次元方向へのエンコ
ードにも使用できる。Next, an embodiment of the present invention will be described with reference to the drawings. FIG. 1 shows a pulse sequence of an embodiment in which the MR imaging method according to the present invention is applied to a field echo method, and this is performed, for example, by a system as shown in FIG. That is, the RF transmitter/receiver 34 sends the first
A high-frequency signal with an envelope as shown in Figure A is sent to excite the spins of a subject (not shown) to be imaged. At this time, under the control of the measurement control device 31, a current having a waveform as shown in FIG. Gz) is applied. Thereafter, currents with waveforms as shown in FIG. 1C and D are sent from the gradient magnetic field power supply 32 to the gradient magnetic field generating coil 33, and a gradient magnetic field for frequency encoding (herein referred to as Gx) and a gradient magnetic field for phase encoding ( Here, Gy) is applied. These waveforms and timings are determined by the measurement control device 31. The NMR signal generated from the object is sent to the R and F coils 3.
5, and is taken into the computer 36 via the RF transmitter/receiver 34, where the data is collected. The subject to be imaged is usually a human body, and generally, as shown in FIG. 2, the horizontal and vertical sizes of the cross-sectional image 21 are not always the same in many parts, so the horizontal size Lx is larger than the vertical size Ly. Therefore, it is not important that the imaging field 22 has Lx=Ly, but rather the isotropy of the pixel itself is more important. Therefore, for this non-square imaging field of view 22, L x / MAT freq = L y / MAT so as not to impair pixel isotropy.
A horizontal (frequency encode direction) matrix MAT freq and a vertical direction (phase encode direction) matrix MAT rhase that satisfy the relational expression of phase are determined. The number of phase encodes is set to be the same as MATphase, and the number of repetitions of the pulse sequence shown in FIG. 1 is decreased by MAT freq-MATphase. This allows the overall imaging time to be shortened. Here, as the matrix weight (VV+++at) that defines the above arbitrary matrix, W mat=MAT
Introduce parameters defined in phase/MAT freq. The phase encode amount ep each time indicated by the shaded area in the first diagram is determined as follows. Since phase control equivalent to sampling at the Nyquist frequency is required, the phase encode waveform Gy needs to satisfy the rIA coefficient of γLyiGydt=π (γ: gyromagnetic ratio). In imaging at the time of arbitrary matrix, the following equation is derived from the above three equations so as to include the W mat parameter. SGy dt= (π/γLx) (1/W mat
) Since π/γLx on the right side of this equation indicates the encoding amount in the normal full encoding method, the one-time phase encoding amount Gp in FIG. It can be seen that it is sufficient to multiply the value by mat. The range of this W mat is Q(W mat ≦ 1, and W mat = 1 means imaging by the full encoding method. The size L of the imaging field of view 22 corresponds to the cross-sectional image 21.
Define x, Ly, MAT freq, MATphas
Once e is determined, W mat is calculated in the computer 36, and the encoded amount Gp is determined. This G
p is sent to the measurement control device 31, and the amount of one encode of the phase encoding gradient magnetic field Gy is controlled accordingly. The raw data obtained in this way is
After dimensional Fourier transformation, the image is reconstructed in a non-square imaging field of view 22 as shown in FIG. In this reconstructed image, the isotropy of pixels is maintained, and the spatial resolution of the image is not degraded. For example, when the arbitrary matrix is set to 70% (W mat = 0.7), the imaging time is 70% (MATphase = 0, 7
MAT freq), the field of view size in the phase encoding direction is 70% of the field of view size in the frequency encoding direction.
(Ly=0.7Lx). Although an embodiment in which the present invention is applied to the field echo method has been described above, it is of course applicable to other imaging methods such as the spin echo method as long as phase encoding is performed. Moreover, it can be easily applied not only to two-dimensional imaging but also to three-dimensional imaging. Furthermore, it can also be used for encoding in three-dimensional directions.
この発明のMR撮像法によると、イメージング視野を非
正方なものとすることを許すような被写体形状であれば
、再構成画像の空間分解能をまったく損なうことなく撮
像時間の短縮を図ることができる。このように撮像時間
を短縮できることにより、患者の拘束時間の短縮、スル
ープットの向上はもとより、被写体の動きなどによるア
ーティファクト成分の混入確率を抑えることができる。According to the MR imaging method of the present invention, if the object shape allows the imaging field to be non-square, the imaging time can be shortened without any loss in the spatial resolution of the reconstructed image. By shortening the imaging time in this way, it is possible not only to shorten the patient's restraint time and improve throughput, but also to suppress the probability of contamination by artifact components due to subject movement.
第1図はこの発明をフィールドエコー法に適用した一実
施例のパルスシーケンスを表わすタイムチャート、第2
図は被写体の断面像とイメージング視野との関係を示す
図、第3図は上記の実施例にかかるシステム構成を示す
ブロック図である。
21・・・断面像、22・・・イメージング視野、31
・、。
測定制御装置、32・・・傾斜磁場電源、33・・・傾
斜磁場発生用コイル、34・・・RF送受信装置、35
・・・RFコイル、36・・・コンピュータ。FIG. 1 is a time chart showing a pulse sequence of an embodiment in which the present invention is applied to a field echo method, and FIG.
The figure is a diagram showing the relationship between a cross-sectional image of a subject and an imaging field of view, and FIG. 3 is a block diagram showing a system configuration according to the above embodiment. 21... Cross-sectional image, 22... Imaging field of view, 31
・、. Measurement control device, 32... Gradient magnetic field power supply, 33... Gradient magnetic field generation coil, 34... RF transmitting/receiving device, 35
...RF coil, 36...computer.
Claims (1)
ズを周波数エンコード方向のサイズよりも小さく定めた
非正方イメージング視野に対して、画素が等方的になる
ように位相エンコード方向のマトリクスと周波数エンコ
ード方向のマトリクスとを決め、位相エンコード数を位
相エンコード方向マトリクスと同数とするとともに、上
記両方向のマトリクスの比により各回の位相エンコード
量を定めるようにしたMR撮像法。(1) For a non-square imaging field of view where the size in the phase encoding direction is set smaller than the size in the frequency encoding direction according to the shape of the subject, the matrix in the phase encoding direction and the frequency encoding are set so that the pixels are isotropic. An MR imaging method in which a direction matrix is determined, the number of phase encodes is the same as the phase encode direction matrix, and the amount of phase encode each time is determined by the ratio of the matrices in both directions.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP1049377A JPH02224737A (en) | 1989-02-28 | 1989-02-28 | Mr imaging method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP1049377A JPH02224737A (en) | 1989-02-28 | 1989-02-28 | Mr imaging method |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH02224737A true JPH02224737A (en) | 1990-09-06 |
JPH0572814B2 JPH0572814B2 (en) | 1993-10-13 |
Family
ID=12829334
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP1049377A Granted JPH02224737A (en) | 1989-02-28 | 1989-02-28 | Mr imaging method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPH02224737A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10132895B2 (en) | 2013-04-25 | 2018-11-20 | General Electric Company | Scan condition determining device, magnetic resonance imaging system, scan condition determining method, and program |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62217950A (en) * | 1986-03-18 | 1987-09-25 | 横河メディカルシステム株式会社 | Nmr imaging apparatus |
-
1989
- 1989-02-28 JP JP1049377A patent/JPH02224737A/en active Granted
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS62217950A (en) * | 1986-03-18 | 1987-09-25 | 横河メディカルシステム株式会社 | Nmr imaging apparatus |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10132895B2 (en) | 2013-04-25 | 2018-11-20 | General Electric Company | Scan condition determining device, magnetic resonance imaging system, scan condition determining method, and program |
Also Published As
Publication number | Publication date |
---|---|
JPH0572814B2 (en) | 1993-10-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6513398B2 (en) | MR image reconstruction using prior information constrained regularization | |
JP2905804B2 (en) | Magnetic resonance imaging method and apparatus | |
US7015696B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
US5704357A (en) | Method of producing MR image and MRI system | |
JPH0732765B2 (en) | Method for correcting phase error caused by NMR apparatus | |
JPH0484937A (en) | MRI compensated for spurious and rapid changes in static magnetic fields during a single MRI sequence | |
US5923789A (en) | Band limited interpolation and projection of spatial 3-D images | |
WO2021196865A1 (en) | 3d oscillating gradient-prepared gradient spin-echo imaging method, and device | |
CN1190572A (en) | Method for Reducing Maxwell Term Artifacts in Fast Spin Echo Magnetic Resonance Images | |
EP3191862B1 (en) | Zero echo time mr imaging | |
US6046591A (en) | MRI system with fractional decimation of acquired data | |
US10175322B2 (en) | Zero echo time MR imaging with sampling of K-space center | |
JPH1033498A (en) | Batch type plural volume angiography using nmr imaging | |
US6975113B1 (en) | Method and system for moving table MRI with partial fourier imaging | |
CN106175765B (en) | Magnetic resonance imaging system and method | |
JP2010075747A (en) | Magnetic resonance imaging method involving sub-sampling and system | |
JP3992934B2 (en) | Magnetic resonance imaging apparatus and method | |
US5349292A (en) | Magnetic resonance imaging apparatus | |
JP3519270B2 (en) | Magnetic resonance imaging apparatus and signal processing method in magnetic resonance imaging apparatus | |
JP4698231B2 (en) | Magnetic resonance diagnostic equipment | |
JP6820876B2 (en) | Magnetic resonance imaging device, sensitivity distribution calculation method, and sensitivity distribution calculation program | |
US4966149A (en) | Reducing the effects of coherence in magnetic resonance imaging | |
JPH02224737A (en) | Mr imaging method | |
JP2002085376A (en) | Nuclear magnetic resonance imaging device and method | |
EP1461632B1 (en) | Parallel mr imaging method |