JP2006519051A - Movable table MRI using sub-sampling - Google Patents
Movable table MRI using sub-sampling Download PDFInfo
- Publication number
- JP2006519051A JP2006519051A JP2006502603A JP2006502603A JP2006519051A JP 2006519051 A JP2006519051 A JP 2006519051A JP 2006502603 A JP2006502603 A JP 2006502603A JP 2006502603 A JP2006502603 A JP 2006502603A JP 2006519051 A JP2006519051 A JP 2006519051A
- Authority
- JP
- Japan
- Prior art keywords
- patient
- region
- main magnet
- data
- profile
- 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.)
- Pending
Links
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/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
- G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
- G01R33/56375—Intentional motion of the sample during MR, e.g. moving table imaging
-
- 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
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
- G01R33/5611—Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (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)
- Vascular Medicine (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
磁気共鳴撮像装置の主磁場の制限された一様領域でサンプリングされた複数の信号から対象の画像を形成する新しい磁気共鳴撮像方法が記載されている。台の上に配置された患者は、主磁石の孔を通って連続的に移動され、前記患者の所定領域のスピンが送信アンテナからの励起パルスにより励起され、この結果、画像が前記制限された領域を大きく越える領域にわたり形成される。データは、前記患者の励起された領域内の維持された磁化及びプリサチュレーション条件を持つ時間の連続した測定のブロックとして定義される複数の受信状況において少なくとも1つの受信アンテナを用いて前記制限された領域においてアンダーサンプリングされる。前記アンダーサンプリングによる折り返しアーチファクトは、前記受信アンテナの既知の感度パターン及び/又は前記受信状況を決定する選択された要素の性質を用いて展開される。A new magnetic resonance imaging method is described that forms an image of an object from a plurality of signals sampled in a limited uniform region of the main magnetic field of a magnetic resonance imaging apparatus. The patient placed on the table is continuously moved through the hole of the main magnet, and the spin of the predetermined area of the patient is excited by the excitation pulse from the transmitting antenna, so that the image is limited. It is formed over a region that greatly exceeds the region. Data was limited using at least one receive antenna in multiple receive situations defined as a block of time continuous measurements with sustained magnetization and pre-saturation conditions in the excited region of the patient Undersampled in region. The aliasing artifacts due to the undersampling are developed using the known sensitivity pattern of the receiving antenna and / or the nature of the selected element that determines the reception situation.
Description
本発明は、磁気共鳴装置の主磁場の制限された一様領域(homogeneity region)内でサンプリングされた複数の信号から対象の画像を形成する磁気共鳴(MR)方法に関し、ここで台の上に配置された患者は、画像が前記磁気共鳴撮像装置の主磁場の制限された領域を大きく越える領域に対して形成されるように、主磁石の孔(bore)を通って連続的に移動される。 The present invention relates to a magnetic resonance (MR) method for forming an image of an object from a plurality of signals sampled in a restricted homogeneity region of a main magnetic field of a magnetic resonance apparatus, wherein the method is on a platform. The placed patient is continuously moved through the bore of the main magnet so that an image is formed for a region that greatly exceeds the limited area of the main magnetic field of the magnetic resonance imaging device. .
本発明は、前記方法を実行するMR装置及びコンピュータプログラムにも関する。 The invention also relates to an MR apparatus and a computer program for performing the method.
患者にとって、MR撮像装置の超電導主磁石(main superconductive magnet)の大きな孔の中に配置されるのは、非常に居心地が悪いので、より短い磁石を使用して、患者台又は寝台を磁石の孔を通して移動する傾向がある。このような装置は、例えば欧州特許出願公開公報EP-A-1024371に記載されている。MR撮像装置において、励起パルスは、磁場が一様である磁石の孔の制限された領域に印加される。収集されたデータサンプルは、フーリエ変換されて前記制限された領域の体積画像を形成する。モータは、関心領域が良い磁場の領域を通過するように患者寝台を連続的に移動する。前記収集されたデータサンプルは、動きを補償するように補正され、この結果、体積画像は、前記制限された領域の体積画像より長く形成される。 It is very uncomfortable for a patient to be placed in a large hole in the main superconductive magnet of an MR imaging device, so a shorter magnet is used to place the patient table or bed into the magnet hole. Tend to move through. Such a device is described, for example, in EP-A-1024371. In an MR imaging device, excitation pulses are applied to a limited area of a magnet hole where the magnetic field is uniform. The collected data samples are Fourier transformed to form a volume image of the restricted area. The motor continuously moves the patient bed so that the region of interest passes through a region of good magnetic field. The collected data samples are corrected to compensate for motion, so that a volumetric image is formed longer than the volumetric image of the restricted area.
上述のMR撮像は、従来の撮像シーケンスで実行され、連続して撮像された体積間の重複さえも中心領域の誤差により生成されたアーチファクトを切り捨てるために好まれるので、非常に長い時間を要する。明らかに、RF送信コイル及びRF受信コイルは、前記主磁石の基準フレーム内、即ち直接的な近傍に配置される。しかしながら、患者の異なる場所に配置されたボディコイルによる撮像は、この参考文献の論点ではなかった。いずれにしろ、ボディコイルの使用は、容易には抑制されることができない追加のアーチファクトを生じる可能性があることが既知である。 The MR imaging described above is performed in a conventional imaging sequence and takes a very long time because it is preferred to truncate artifacts generated by errors in the central region, even the overlap between successively imaged volumes. Obviously, the RF transmitter coil and the RF receiver coil are arranged in the reference frame of the main magnet, ie in the immediate vicinity. However, imaging with body coils located at different locations on the patient was not the issue of this reference. In any case, it is known that the use of body coils can result in additional artifacts that cannot be easily suppressed.
したがって、本発明の目的は、撮像が大幅に短い時間期間に実行されるように比較的短い主磁石及び可動台を用いてMR撮像方法を改良することである。本発明の他の目的は、この方法を実行するMR装置及びコンピュータプログラムを提供することである。 Accordingly, it is an object of the present invention to improve the MR imaging method using a relatively short main magnet and movable platform so that imaging is performed in a significantly shorter time period. Another object of the present invention is to provide an MR apparatus and a computer program for executing this method.
これらの目的は、請求項1に記載の方法により、請求項11に記載のMR装置により、及び請求項13に記載のコンピュータプログラムにより達成される。
These objects are achieved by the method according to claim 1, by the MR apparatus according to claim 11 and by the computer program according to
本発明は、制限された一様領域を持つ、より短い主磁石が使用されることができるが、データがアンダーサンプリングのために大幅に短い時間でサンプリングされることができ、したがって走査時間が、可動台又はベッドを有する従来のMR撮像に対して大幅に減少されるという主要な利点を有する。本発明の他の利点は、視野の方向が如何なる方向にも向けられることができ、例えば流れアーチファクトが適切な方向に向けられることができ、又は前記患者の組織構造に対する好適な向きに単一のスライスを向けることができることである。 The present invention can use a shorter main magnet with a limited uniform area, but the data can be sampled in a much shorter time due to undersampling, so the scan time is It has the major advantage of being greatly reduced over conventional MR imaging with a movable platform or bed. Another advantage of the present invention is that the field of view can be oriented in any direction, for example, flow artifacts can be oriented in a suitable direction, or a single orientation in a suitable orientation relative to the patient's tissue structure. The slices can be directed.
本発明のこれら及び他の利点は、従属請求項及び本発明の模範的実施例が添付図面に関して記載される以下の記載に開示される。 These and other advantages of the present invention are disclosed in the dependent claims and the following description in which exemplary embodiments of the invention are described with reference to the accompanying drawings.
基本原理及び記号
図1において、患者1は可動台又はベッド2に固定して配置される。局所コイルLC1ないしLC5のセットは患者1に対して固定して配置され、したがって配置“患者1”、“ベッド2”及び“局所コイルLC1ないしLC5”は、前記患者に関する座標系3において互いに固定された関係で配置される。ベクトル
は、この座標系又は基準系3に対する位置を示す。更に、主磁場4には、ベッドが矢印6の方向に前後に移動されることができる孔5が設けられる。磁石/勾配システム7は、主磁場4及び複数の受信コイル8−一般には1又は0−、並びに送信コイル(描かれていない)を有する。これらの受信コイル8及び送信コイルは、“大域的(global)”であると示される。磁石/勾配システム7の要素は全て、主磁石4に関する他の座標系9により定義される他の固定された関係で関連する。受信コイル8及び前記送信コイルは大きく、前記ベッドの運動の方向(矢印6)に滑らかな感度パターンを持つと更に仮定される。磁石座標系9は、患者座標系3に対する距離
によりオフセットされる。患者1の性質、例えば体の陽子密度は、患者座標系3に明らかに与えられ、
と表現される。理想的には、i=1〜5である局所コイルLCiの1つの受信感度パターン
についても同様である。これに対し、受信コイル7の感度パターンは、
と記載される。
Basic Principle and Symbol In FIG. 1, a patient 1 is fixedly placed on a movable table or
Indicates the position relative to this coordinate system or reference system 3. Furthermore, the main magnetic field 4 is provided with a
Offset by The nature of the patient 1, for example the proton density of the body, is clearly given in the patient coordinate system 3,
It is expressed. Ideally, one reception sensitivity pattern of the local coil LC i where i = 1 to 5
The same applies to. On the other hand, the sensitivity pattern of the receiving coil 7 is
It is described.
ここで表現“受信状況(receive situation)”が導入され、これは、時間の連続した測定のブロックを意味し、この間に台2が主磁石4の“一様サイズ”の比較的小さな距離(例えば半分若しくは3分の1又はその間)だけ移動する。単一の受信状況の間に、(プリサチュレーション・スラブ(presaturation slab)、体積選択等のような)全ての磁化プレパレーション(magnetization preparations)は、おおむね患者1と共に“移動”する。単一の受信状況の間に、全ての関心のあるk空間が、(潜在的に重要である)アンダーサンプリングを除き、カバーされる。k空間全体の完全走査は、Kの受信状況をカバーする。原理的に、MR収集シーケンスは、前記“磁化プレパレーション”が位置を変えることができ、収集されたk空間グリッドにオフセットが存在する可能性があることを除き、全ての受信状況に対して同一である。しかしながら、これは概念の改良である。
The expression “receive situation” is introduced here, which means a block of continuous measurement of time during which the
受信状況の間、コイルLC1ないしLC5は、患者1から情報を受信する。(I+J)個の情報セットが供給され、ここでIは局所受信コイルの数であり、Jは大域受信コイルの数である。これらのセットの幾つかは、例えば受信コイルLC1ないしLC5が完全に主磁石4の一様体積(homogeneity volume)の外である場合に、実際には0である可能性がある。K個の受信状況が存在する。したがって、全体で、情報はK・(I+J)個の“受信インスタンス(receive instances)”に集められる。 During the reception situation, the coils LC 1 to LC 5 receive information from the patient 1. (I + J) information sets are provided, where I is the number of local receive coils and J is the number of global receive coils. Some of these sets may actually be zero if, for example, the receiving coils LC 1 to LC 5 are completely outside the homogeneity volume of the main magnet 4. There are K reception situations. Thus, in total, information is collected in K · (I + J) “receive instances”.
本発明の基本アイデア
臨床医は、患者1の関心のある比較的大きな体積からの情報を望む。データは、台又はベッド2が主磁石4に対して移動している間に収集される。磁石4に対する変位
はいつでも既知である。MR収集の全体は、複数の受信状況に分割される。各状況の間に、励起又はプリサチュレーションプロファイルは、患者1と共に“移動”する。前記収集されたデータは、前記データが収集された台2のオフセットに対して位相補正される。受信状況の間に、患者1からの情報は、如何なる向きの如何なるMR撮像シーケンスを使用しても収集される。実際には、アーチファクト(例えば流れアーチファクト)を適切な方向に向けるか、又は単一のスライスを特定の人間組織構造に対する好適な向きに向けることは有用であるかもしれない。しかしながら、前記情報は、著しくアンダーサンプリングされている可能性がある。この方法の例は図2にスケッチされている。受信状況ごとに対して、ボックス11により示される体積を名目上収集する。この大多数は、1次元の行、2次元の行、又は例えば螺旋形収集の場合の折り返し位置(folding location)の連続体であることさえあり得る。
Basic Idea Clinicians of the Invention The clinician wants information from a relatively large volume of patient 1's interest. Data is collected while the platform or
Is always known. The entire MR acquisition is divided into a plurality of reception situations. During each situation, the excitation or pre-saturation profile “moves” with the patient 1. The collected data is phase-corrected with respect to the offset of the table 2 from which the data was collected. During the reception situation, information from patient 1 is collected using any MR imaging sequence in any orientation. In practice, it may be useful to direct artifacts (eg, flow artifacts) in the proper orientation, or to direct a single slice in a suitable orientation for a particular human tissue structure. However, the information may be significantly undersampled. An example of this method is sketched in FIG. For each reception situation, the volume indicated by box 11 is nominally collected. The majority can be a one-dimensional row, a two-dimensional row, or even a continuum of folding locations, for example in the case of a helical collection.
ここで、受信状況は同じサブサンプリングパターンを有する。異なる受信状況間の唯一の違いは、以下の選択された要素、即ち、
1.励起又は飽和プロファイル(患者1に対する)
2.前記送信コイルの送信パターン
3.固定された受信コイル8の受信パターン
4.受信器システムの周波数応答パターン
の性質に対する重み付けである。
Here, the reception situation has the same sub-sampling pattern. The only difference between the different reception situations is the following selected elements:
1. Excitation or saturation profile (for patient 1)
2. 2. Transmission pattern of the transmission coil 3. Reception pattern of fixed reception coil 8. A weighting for the nature of the frequency response pattern of the receiver system.
項目2及び3に関して、理論的には、可動台2を用いて、これらのパターンは1つの状況内でも変化するが、状況ごとの運動が甚だしく大きくない場合には、これは無視されることができることに注意する。
With respect to
加えて、受信コイル(局所又は大域)は、独自の感度パターンを有する。したがって受信コイル及び受信状況の組み合わせである受信インスタンスは、前記コイル感度パターンと前記受信状況の組み合わされたパターン(励起、飽和、送信及び周波数応答)との積である異なる“全体パターン”を“見る”。 In addition, the receiving coil (local or global) has its own sensitivity pattern. Thus, a receive instance that is a combination of a receive coil and a receive situation “sees” a different “overall pattern” that is the product of the coil sensitivity pattern and the combined pattern of the receive situation (excitation, saturation, transmission and frequency response). ”.
折り返し点(folding point)12の大多数は、受信インスタンスごとの全ての前記“全体パターン”を使用して展開される(unfolded)。これは、SENSE法(例えばK. Pruessmann他、Proc. ISMRM、1998、要約 pp.579、799、803及び2087)のように行われるが、形式的には、展開マトリクス(unfolding matrix)は、より大きい可能性がある。そうでなければ、式の系は適度に安定であることができ、即ち、例えば各状況に対して、多くの折り返し点12は、励起スラブの外になるので、実際には0パターンを“見る”かもしれない。前記台の運動はフレキシブルであり、即ち、これは線形、非線形且つ非単調又は2次元であることができることに更に注意すべきである。制限は、これが1つの受信状況の間に過剰に、例えば前記一様体積の半分も移動してはならないことである。
The majority of
パターンの測定
較正
SENSE法において、コイル感度プロファイルのマップは利用可能でなければならない。このようなマップは、基準又は明確なファントムを有する粗い較正測定値で与えられることができる。他の可能性は、アレイコイル素子からの信号と、前記台上の前記患者の基準コイルからの信号との比の低解像度マップを使用することである。本発明において、同様な基準測定が与えられる。データは一般に2次元の台運動を含む台運動パターンにおいて収集されることができるが、ここで線形の1次元運動が、データ収集に使用される。
Pattern measurement
In the calibration SENSE method, a map of the coil sensitivity profile must be available. Such a map can be given with coarse calibration measurements with a reference or a clear phantom. Another possibility is to use a low resolution map of the ratio of the signal from the array coil element to the signal from the patient's reference coil on the table. In the present invention, a similar reference measurement is given. Data can generally be collected in a table motion pattern that includes a two-dimensional table motion, where linear one-dimensional motion is used for data collection.
関連する原理のより良い理解のために、初めにレビューすると、
1.前記患者の各位置
において、幾つかの位置依存組織性質(陽子密度、T1、T2、流れ、拡散...に関する全てのデータ)が既知である。
2.前記患者の各位置
において、局所コイルLCiの感度パターン
が得られることができる。大域コイルに対する結合のために
に対する依存性が存在するかもしれないが、これはここでは無視される。
3.全てのコイル(大域又は局所、送信又は受信)のパターンは、特に3T以上の主磁場における誘電共振により影響を受ける可能性がある。この効果は、全てのコイルのコイル感度に対して純粋に乗算的であると見なされることができ、−前記患者及び前記患者における位置のみに依存する効果であり、したがって前記コイル又はコイル位置に依存しないと仮定される。1.5Tの主磁場に対して、前記効果は無視することができる。誘電共振効果は、
で示される。各コイルは実際には感度
を“見る”。
4.前記大域コイルは、磁石フレーム内の前記患者の実際の位置に独立であると近似される感度パターンをも有する(もちろん明示的に計上される誘電共振効果を除く)。この感度パターンは、
と表現される。
5.全体の受信パターンは、所定の読み出し勾配方向において
と示される受信器周波数応答の効果を含むこともできる。
For a better understanding of the relevant principles, first review
1. Each position of the patient
Several position-dependent tissue properties are known (all data on proton density, T1, T2, flow, diffusion ...).
2. Each position of the patient
Sensitivity pattern of local coil LC i
Can be obtained. For coupling to the global coil
There may be a dependency on, but this is ignored here.
3. The pattern of all coils (global or local, transmit or receive) can be affected by dielectric resonance, especially in main magnetic fields above 3T. This effect can be regarded as purely multiplicative with respect to the coil sensitivity of all coils, and is an effect that depends only on the patient and the position in the patient and is therefore dependent on the coil or coil position. It is assumed not to. For a main magnetic field of 1.5 T, the effect can be ignored. Dielectric resonance effect is
Indicated by Each coil is actually sensitive
I see".
4). The global coil also has a sensitivity pattern that is approximated to be independent of the patient's actual position within the magnet frame (except of the explicitly accounted for dielectric resonance effect). This sensitivity pattern is
It is expressed.
5. The overall reception pattern is determined in the predetermined readout gradient direction.
The effect of the receiver frequency response indicated as
所定の測定タイプに対して、横向き磁化(transverse magnetisation)は、送信パターン−即ち大域コイルが使用される場合の磁石座標系9における位置−、及び組織性質−即ち患者座標系3における位置−の両方に依存する非線形効果である。これは非常に扱うのが難しいので、近似が必要であり、これは非常に小さな先端角度(およそ5度以下)を持つFFE測定を使用することにより達成されることができる。この場合、前記横向き磁化は、
として近似されることができ、ここで
は、所定の測定に対する組織性質である“スピン密度”関数である(提案されている場合には、主に陽子密度に依存する)。αは
に比例するように近似された前記先端角度である。
For a given measurement type, the transverse magnetisation is both a transmission pattern-i.e. a position in the magnet coordinate system 9 when a global coil is used-and a tissue property-i.e. a position in the patient coordinate system 3-. This is a nonlinear effect that depends on. Since this is very difficult to handle, an approximation is necessary, which can be achieved by using FFE measurements with very small tip angles (approximately 5 degrees or less). In this case, the transverse magnetization is
Can be approximated as
Is a “spin density” function that is a tissue property for a given measurement (in the proposed case it depends mainly on the proton density). α is
The tip angle approximated to be proportional to.
問題点
後でより詳しく記載される提案された方法の再構成に対して、理想的には
の全ての情報を知りたい(しかしながら最後の1つは重要ではない)。
の値は事前に知られていると仮定され、
の値は(少なくとも原理的には)較正測定の関心のある結果ではない。
Against reconstruction of the proposed method problems described in more detail below, ideally
I want to know all the information (but the last one is not important).
Is assumed to be known in advance,
The value of (at least in principle) is not an interesting result of the calibration measurement.
この問題点は、これら全てのパラメータが測定されるわけではないことである。所定の台位置
に対し、以下の値、即ち、
−前記局所コイルからの信号:
−前記大域コイルからの信号:
が測定される。
これら2つの信号の単純な商は、
となり、これは、
が一様でない場合に、直接的には役に立たない。
The problem is that not all these parameters are measured. Predetermined platform position
For the following values:
The signal from the local coil:
The signal from the global coil:
Is measured.
The simple quotient of these two signals is
And this is
If is not uniform, it is not directly useful.
第1の解決策
前記大域コイルのパターンは、既知であると仮定される(少なくとも1つ)。実際には、前記関数はロード依存である可能性があるが、これは、感度関数の較正されたロード依存性が存在し、前記ロードが容易に測定されることができるという事実により克服されることができる。
First Solution The global coil pattern is assumed to be known (at least one). In practice, the function may be load dependent, but this is overcome by the fact that there is a calibrated load dependency of the sensitivity function and the load can be easily measured. be able to.
複数の離散的な基準走査セグメントが収集される。このようなセグメントの間に、前記台は静止している。各セグメントは結果として完全な“画像”を生じる。この場合、前記台は、複数のセグメントにわたり、したがって
の複数の異なる値に対して進められる。
の情報はマトリクス
を計算することを可能にする。数学的に、各セグメントは、
の全範囲にわたり結果を与えるが、実際には、各セグメントは、全ての
の部分範囲にわたってのみ正確な結果を提供する。これらの部分的に重複する範囲は、最小二乗フィッティングを使用して結合されることができ、結果として
の完全なマップを生じる。
A plurality of discrete reference scan segments are collected. During such a segment, the platform is stationary. Each segment results in a complete “image”. In this case, the platform spans multiple segments and is therefore
Is advanced for a number of different values.
Information is a matrix
Makes it possible to calculate Mathematically, each segment is
Results in the entire range, but in practice, each segment
Provide accurate results only over a sub-range. These partially overlapping ranges can be combined using least square fitting, and as a result
Produces a complete map of
第2の解決策
3次元測定は、周波数符号化が前記台運動方向−z方向に平行になるように適応される。サンプリング帯域幅は、この方向において非常に高い。他の2つの方向は位相符号化であり、一方は台運動と比較して“速い”。“遅い”位相符号化方向はkyと称される。患者座標系9において、主磁石4は前記患者に対してゆっくりと移動する。前記磁石の運動中に、プロファイルが、線形に増加するkyに対して収集される。図3において、ラインセグメント13は、主磁石4の一様体積により検出される患者1の部分を示す。
The second solution 3D measurement is adapted so that the frequency coding is parallel to the table motion direction-z direction. The sampling bandwidth is very high in this direction. The other two directions are phase encoding, and one is “fast” compared to table motion. "Slow" phase encoding direction is referred to as k y. In the patient coordinate system 9, the main magnet 4 moves slowly with respect to the patient. During movement of the magnet, the profile is collected for k y linearly increasing. In FIG. 3, the
Lは、主磁石4が位相符号化プロファイルの完全なセットに対してカバーする距離である。前記台の速度は、Lが最大で前記一様体積の半分−しかしそれ以下(例えば1/4)が好ましい−であるように整えられる。図3においてLは1/2に選択される。結果は、プロファイルが主磁石4の少なくとも2つの異なる位置(ここで無視される外縁領域を除く)に対してサンプリングされることになる。例えば、プロファイルkyは、変位
及び変位
の両方に対して測定されている。これは、主磁石4が位置
(この値が
と
との間にある限り)にある場合に測定された前記プロファイルを近似することを可能にする。これは、補間法(例えば、2つのプロファイルが利用可能な場合には線形補間、又は4つのプロファイルが存在する場合には三次スプライン)により行われることができる。これはプロファイルごとに行われることができ、磁石オフセット
における感度較正測定の完全な再構成を可能にする。
L is the distance that the main magnet 4 covers for a complete set of phase encoding profiles. The speed of the pedestal is adjusted so that L is maximum and half of the uniform volume-but less (eg 1/4) is preferred. In FIG. 3, L is selected to be 1/2. The result is that the profile is sampled for at least two different positions of the main magnet 4 (except for the outer edge region which is ignored here). For example, the profile k y, the displacement
And displacement
Measured against both. This is because the main magnet 4 is located
(This value is
When
It is possible to approximate the measured profile as long as it is between. This can be done by interpolation (eg, linear interpolation if two profiles are available, or cubic spline if there are four profiles). This can be done per profile, and magnet offset
Allows a complete reconstruction of the sensitivity calibration measurement.
しかしながら、これは、
の如何なる値に対しても行われることができる。この原理は、
−前記局所コイルからの信号:
と、
−前記大域コイルからの信号:
と、
を収集することを可能にする。
これは、“問題点”セクションとは反対に、
の1つの値だけでなく、全ての値に対して当てはまる。これは、
−前記局所コイルからの信号を
及び
に分離し、
−前記大域コイルからの信号を
及び
に分離する
ことを可能にする。
これは、
及び
の算出を可能にする。実際には、Fが既知である場合、
も算出されることができるが、不幸なことに誘電共振
は算出されることができないが、これは結果に対して重要ではない。
However, this
Can be done for any value of This principle is
The signal from the local coil:
When,
The signal from the global coil:
When,
Makes it possible to collect.
This is contrary to the “problem” section,
This applies to all values, not just one value. this is,
-The signal from the local coil
as well as
Separated into
-The signal from the global coil
as well as
Makes it possible to separate.
this is,
as well as
Enables the calculation of. In fact, if F is known,
Can also be calculated, but unfortunately dielectric resonance
Cannot be calculated, but this is not important to the result.
実際のデータの受信及び再構成
セクション“本発明の基本アイデア”において、データの受信及び再構成は、大いに説明されている。前記データは、プロファイル間の比較的大きなステップを使用して収集され、結果として(名目上)小さな“折り返し体積(folding volume)”、即ちたくさんの折り返し(folding)を生じる。非カーテシアン・シーケンス(non-cartesian sequences)に対して、これは、データが比較的まばらにサンプリングされることを意味する。k空間の全範囲が1走査状況(scan situation)の間に収集され、1走査状況の間の変位は前記一様体積の一部(例えば1/3)である。入力サンプルは、
の非常に独特な値で収集される。これは、入力データのサンプルごとに
を乗算することにより考慮に入れられる。ここで、Tは、患者座標から斜め走査(oblique scan)に関連する走査座標への座標変換である。この演算は、たとえ前記一様体積のかなり外であっても、前記収集されたデータを前記患者の中心に“移動”する。明らかに、プロファイルごとに固定された
で補正することにより単純化することができる。
In the actual data reception and reconstruction section “Basic Ideas of the Present Invention”, data reception and reconstruction is greatly explained. The data is collected using relatively large steps between the profiles, resulting in a (nominal) small “folding volume”, ie a lot of folding. For non-cartesian sequences, this means that the data is sampled relatively sparsely. The entire range of k-space is collected during a scan situation, and the displacement during a scan situation is part of the uniform volume (eg 1/3). The input sample is
Collected with very unique values. This is done for each sample of input data
Is taken into account by multiplying by. Here, T is a coordinate transformation from patient coordinates to scan coordinates associated with an oblique scan. This operation “moves” the collected data to the center of the patient, even well outside the uniform volume. Obviously fixed per profile
It can be simplified by correcting with.
前記走査状況は、平均オフセット
において収集された。これは、
−当該状況に対する送信プロファイル:
と、
−同様に大域コイル受信プロファイルと、
−同様に当該状況に対する周波数応答関数と、
を推定することに関連する。
The scanning situation is the average offset
Collected in this is,
-Transmission profile for the situation:
When,
-Similarly, the global coil reception profile,
-Similarly the frequency response function for the situation,
Related to estimating.
励起プロファイル及びプリサチュレーションプロファイルは、前記患者座標系に対して固定され、即ち主磁石4から見ると台2と共に移動する。したがって、これらは
において既知である。局所コイルプロファイルも
において既知である(較正セクション参照)。これらの全てのコイルプロファイルは、各コイルに対する全体的な受信パターンに増やされることができる。全体として、Kの受信状況が存在する場合、情報はK・(I+J)の受信インスタンスに集められている。各インスタンスは異なる全体的な受信パターンを持つ。これは、患者1全体に対してK・(I+J)より少ない折り返された点が存在する場合に、SENSEのような再構成を可能にする。原理的に、この収集は全ての受信状況に対して同じままである。改良として、位相符号化が少量だけオフセットされることができる。これは、他の全ての述べられた符号化の上に追加の位相符号化を与え、SENSE再構成の安定性を向上する。
The excitation profile and pre-saturation profile are fixed with respect to the patient coordinate system, i.e. move with the
Is known. Local coil profile
(See calibration section). All these coil profiles can be augmented to an overall reception pattern for each coil. Overall, if there are K reception situations, the information is collected in K · (I + J) reception instances. Each instance has a different overall reception pattern. This allows a reconstruction like SENSE when there are fewer folded points than K · (I + J) for the entire patient 1. In principle, this collection remains the same for all reception situations. As an improvement, the phase encoding can be offset by a small amount. This provides additional phase encoding over all other described encodings and improves the stability of the SENSE reconstruction.
幾何学的補正
上の推論の全てが完全に線形の勾配システムを仮定している。
の台変位は、(a)広範囲で上のように処理されている異なる重み付け、及び(b)サンプルに対する
の位相修正を生じると仮定する。不幸なことに、幾何学的歪は、他の重大な複雑にする要素である。オフセット
で収集された画像は、たとえ全ての重み付け及び変位
が計上されたとしても、オフセット
で収集された画像と比較される場合に幾何学的に歪んでいる。
All of the geometric correction inferences assume a perfectly linear gradient system.
Table displacements are: (a) different weightings that have been processed extensively as above, and (b) for samples
Is assumed to result in a phase correction of Unfortunately, geometric distortion is another significant complicating factor. offset
The images collected at are all weighted and displaced
Offset even if
Are geometrically distorted when compared to images collected at
この問題は、前記励起及び送信プロファイルが大きくなる傾向にある場合に前記歪が小さくなる傾向にあるので、解決されることができる。これは、(a)幾何学的に非一様な領域に多すぎる励起が存在しないようにシーケンスを作成することにより、及び(b)これらの領域において低い送信感度を持つシステムを設計することにより果たされることができる。この解決策は、反復的再構成であり、この説明は図4に視覚的に与えられる。主要なアイデアは、基本的に、初めに幾何学的歪の存在を無視することにより結果画像を再構成することである。この結果について、全ての画像領域は、走査状況が存在するだけ多くの方向に歪まされる。歪んでいないバージョンを減算した歪んだバージョンは、“折り返し歪誤差”の推定値を与える。最終画像に対する受信インスタンスの寄与(又は重み)が既知である場合、前記折り返し歪誤差の推定値はこの結果から抽出されることができる。 This problem can be solved because the distortion tends to decrease when the excitation and transmission profiles tend to increase. This is due to (a) creating sequences so that there are not too many excitations in geometrically non-uniform regions, and (b) designing systems with low transmission sensitivity in these regions. Can be fulfilled. This solution is an iterative reconstruction and this explanation is given visually in FIG. The main idea is basically to reconstruct the resulting image by first ignoring the presence of geometric distortion. For this result, all image areas are distorted in as many directions as there are scanning conditions. The distorted version obtained by subtracting the undistorted version gives an estimate of the “folding distortion error”. If the received instance's contribution (or weight) to the final image is known, an estimate of the aliasing error can be extracted from this result.
図4(a)において、概略的には異なる位置において2つの部分、即ちブロック16として示される第1部分及び心臓17として示される第2部分からなる非常に単純化された患者15が描かれている。図4(b)は、ブロック16が主として前記主磁石の中心にあり、心臓18がより離れており、したがって歪んでいるときの撮像を示す。図4(c)に示されるように、歪んだ心臓18は、アンダーサンプリングによりブロック16に折り重なる。まだ、心臓18は、より低い励起を見たので、より薄れている。図4(d)による他の受信状況において、心臓17は主磁石4の中心にあり、ブロック19はより離れており、したがって歪んでいる。アンダーサンプリングにより、これは、図4(e)に示されるような受信状況で測定される。したがって、図4(c)及び図4(e)で画像が実際に測定され、図4(f)は、図4(c)及び図4(e)の2つの線形結合からなるSENSE展開(unfolding)の後の結果を示す。これは不完全であり、例えば右の部分は、ブロック16と、減算されていない心臓18と、“間違って減算された”心臓17(破線)と、2次の誤差(点線)とからなる。図4(g)には、(f)が可能である条件下の(b)と同じ様式の意図的な歪が示される。図4(h)は、(f)の一部分が折り返された対の部分から減算され、(g)の一部が加算され、結果として2次アーチファクトのみを生じる最終結果を示す。
In FIG. 4 (a) a very
Claims (13)
主磁石と、
前記患者の所定領域のスピンを励起する送信アンテナと、
主磁場の制限された一様領域において複数の信号をサンプリングする少なくとも1つの受信アンテナと、
患者を支える台と、
前記主磁石の孔を通って前記台を連続的に移動する手段と、
基準点に対する連続的に移動する前記台の位置に前記送信アンテナの周波数プロファイル及び前記受信アンテナの周波数プロファイルを適合する手段と、
前記患者の励起された領域における維持された磁化及びプリサチュレーション条件を持つ時間の連続した測定のブロックとして定義される所定の受信状況において前記制限された領域内でデータをサンプリングする手段と、
前記受信アンテナの感度パターン及び/又は前記受信状況を決定する選択された要素の性質を用いてアンダーサンプリングによる折り返しアーチファクトを展開する手段と、
を有する装置。 In a magnetic resonance imaging apparatus for obtaining MR images from a plurality of signals,
A main magnet;
A transmitting antenna that excites spins of a predetermined region of the patient;
At least one receiving antenna for sampling a plurality of signals in a limited uniform region of the main magnetic field;
A platform to support the patient,
Means for continuously moving the platform through a hole in the main magnet;
Means for adapting the frequency profile of the transmitting antenna and the frequency profile of the receiving antenna to the position of the stage that moves continuously with respect to a reference point;
Means for sampling data within the limited region in a predetermined reception situation defined as a block of time continuous measurements with sustained magnetization and presaturation conditions in the excited region of the patient;
Means for developing aliasing artifacts due to undersampling using a sensitivity pattern of the receiving antenna and / or a property of a selected element that determines the reception situation;
Having a device.
主磁石により主磁場を作成するステップと、
送信アンテナにより前記患者の所定領域のスピンを励起するステップと、
少なくとも1つの受信アンテナにより前記主磁場の制限された一様領域で複数の信号をサンプリングするステップと、
患者を支える台を前記主磁石の孔を通って連続的に移動するステップと、
基準点に対する連続的に移動する前記台の位置に前記送信アンテナの周波数プロファイル及び前記受信アンテナの周波数プロファイルを適合するステップと、
前記患者の励起された領域における維持された磁化及びプリサチュレーション条件を持つ時間の連続した測定のブロックとして定義される所定の受信状況において前記制限された領域でデータをアンダーサンプリングするステップと、
前記受信アンテナの感度パターン及び/又は前記受信状況を決定する選択された要素の性質を用いて前記アンダーサンプリングによる折り返しアーチファクトを展開するステップと、
をコンピュータに制御させるコンピュータ読取可能プログラム手段を有するコンピュータプログラム。 In a computer program stored on a computer usable medium for forming an image using a magnetic resonance method,
Creating a main magnetic field with a main magnet;
Exciting a spin of a predetermined region of the patient with a transmitting antenna;
Sampling a plurality of signals in a limited uniform region of the main magnetic field with at least one receiving antenna;
Continuously moving the platform supporting the patient through the hole in the main magnet;
Adapting the frequency profile of the transmitting antenna and the frequency profile of the receiving antenna to the position of the platform that moves continuously with respect to a reference point;
Undersampling data in the limited region in a predetermined reception situation defined as a block of time continuous measurements with sustained magnetization and pre-saturation conditions in the excited region of the patient;
Unfolding aliasing artifacts due to the undersampling using a sensitivity pattern of the receiving antenna and / or a property of a selected element that determines the reception situation;
A computer program comprising computer readable program means for controlling a computer.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP03100506 | 2003-02-28 | ||
PCT/IB2004/050131 WO2004077086A1 (en) | 2003-02-28 | 2004-02-18 | Moving table mri with subsampling |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2006519051A true JP2006519051A (en) | 2006-08-24 |
Family
ID=32921620
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2006502603A Pending JP2006519051A (en) | 2003-02-28 | 2004-02-18 | Movable table MRI using sub-sampling |
Country Status (4)
Country | Link |
---|---|
US (1) | US7417429B2 (en) |
EP (1) | EP1599741A1 (en) |
JP (1) | JP2006519051A (en) |
WO (1) | WO2004077086A1 (en) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7417429B2 (en) * | 2003-02-28 | 2008-08-26 | Koninklijke Philips Electronics N.V. | Moving table MRI with subsampling in parallel |
WO2005047918A1 (en) * | 2003-11-14 | 2005-05-26 | Koninklijke Philips Electronics, N.V. | Fast continuous moving bed magnetic resonance imaging with multiple stationary receive coils for sense reconstruction |
EP1763681A1 (en) * | 2004-06-28 | 2007-03-21 | Koninklijke Philips Electronics N.V. | Parallel magnetic resonance imaging |
WO2006077936A1 (en) * | 2005-01-21 | 2006-07-27 | Hitachi Medical Corporation | Magnetic resonance imaging method and device |
DE102005018939B4 (en) * | 2005-04-22 | 2007-09-20 | Siemens Ag | Improved MRI imaging based on conventional PPA reconstruction techniques |
JP5243421B2 (en) | 2006-07-18 | 2013-07-24 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Multi-coil MRI artifact suppression |
WO2008075265A1 (en) * | 2006-12-19 | 2008-06-26 | Koninklijke Philips Electronics N.V. | Motion correction in a pet/mri hybrid imaging system |
US8320647B2 (en) | 2007-11-20 | 2012-11-27 | Olea Medical | Method and system for processing multiple series of biological images obtained from a patient |
WO2010001285A1 (en) * | 2008-07-04 | 2010-01-07 | Koninklijke Philips Electronics N.V. | Mr imaging with extended field of view |
DE102014212943B4 (en) * | 2014-07-03 | 2016-11-24 | Siemens Healthcare Gmbh | Magnetic resonance imaging taking into account different frequency coding patterns |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001190517A (en) * | 2000-01-06 | 2001-07-17 | Shimadzu Corp | Mr imaging device |
JP2002010992A (en) * | 2000-04-25 | 2002-01-15 | Toshiba Corp | Magnetic resonance imaging device and magnetic resonance imaging method |
JP2003038459A (en) * | 2001-07-13 | 2003-02-12 | Koninkl Philips Electronics Nv | Rf coil system for mr apparatus |
JP2004097826A (en) * | 2002-09-12 | 2004-04-02 | General Electric Co <Ge> | Method and system for imaging extended volume using mri with parallel reception |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2345139A (en) | 1998-12-24 | 2000-06-28 | Marconi Electronic Syst Ltd | MRI apparatus with continuous movement of patient |
JP2003500134A (en) | 1999-05-20 | 2003-01-07 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Magnetic resonance imaging with subsampling |
WO2000072050A1 (en) | 1999-05-24 | 2000-11-30 | Kyriakos Walid E | Method and apparatus for parallel data acquisition from a mri coil array |
US6946836B2 (en) | 2000-04-25 | 2005-09-20 | Kabushiki Kaisha Toshiba | Magnetic resonance imaging involving movement of patient's couch |
EP1307757B1 (en) | 2000-07-31 | 2012-09-12 | Koninklijke Philips Electronics N.V. | Magnetic resonance imaging method with sub-sampled acquisition |
JP4515616B2 (en) * | 2000-09-25 | 2010-08-04 | 株式会社東芝 | Magnetic resonance imaging system |
US6445181B1 (en) | 2000-11-09 | 2002-09-03 | The Board Of Trustees Of The Leland Stanford Junior University | MRI method apparatus for imaging a field of view which is larger than a magnetic field |
DE60127382T2 (en) | 2000-11-24 | 2007-12-06 | Koninklijke Philips Electronics N.V. | METHOD FOR OBTAINING IMAGES MAGNETIC RESONANCE BY UNDER-SAMPLE IN A VERTICAL FIELD MRI DEVICE |
DE10059772A1 (en) | 2000-11-30 | 2002-06-13 | Philips Corp Intellectual Pty | MR image reconstruction |
US6897655B2 (en) * | 2001-03-30 | 2005-05-24 | General Electric Company | Moving table MRI with frequency-encoding in the z-direction |
US6794869B2 (en) * | 2001-03-30 | 2004-09-21 | General Electric Company | Moving table MRI with frequency-encoding in the z-direction |
US6912415B2 (en) * | 2001-04-09 | 2005-06-28 | Mayo Foundation For Medical Education And Research | Method for acquiring MRI data from a large field of view using continuous table motion |
US6564082B2 (en) * | 2001-04-19 | 2003-05-13 | General Electric Company | Method for incremental field-of-view-MR imaging |
DE10119660B4 (en) * | 2001-04-20 | 2006-01-05 | Siemens Ag | Method for the rapid acquisition of a magnetic resonance image |
DE10128534C2 (en) * | 2001-06-13 | 2003-04-30 | Siemens Ag | Magnetic resonance imaging method using multiple receiving antennas |
US7417429B2 (en) * | 2003-02-28 | 2008-08-26 | Koninklijke Philips Electronics N.V. | Moving table MRI with subsampling in parallel |
-
2004
- 2004-02-18 US US10/547,196 patent/US7417429B2/en not_active Expired - Fee Related
- 2004-02-18 WO PCT/IB2004/050131 patent/WO2004077086A1/en active Application Filing
- 2004-02-18 JP JP2006502603A patent/JP2006519051A/en active Pending
- 2004-02-18 EP EP04712128A patent/EP1599741A1/en not_active Withdrawn
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001190517A (en) * | 2000-01-06 | 2001-07-17 | Shimadzu Corp | Mr imaging device |
JP2002010992A (en) * | 2000-04-25 | 2002-01-15 | Toshiba Corp | Magnetic resonance imaging device and magnetic resonance imaging method |
JP2003038459A (en) * | 2001-07-13 | 2003-02-12 | Koninkl Philips Electronics Nv | Rf coil system for mr apparatus |
JP2004097826A (en) * | 2002-09-12 | 2004-04-02 | General Electric Co <Ge> | Method and system for imaging extended volume using mri with parallel reception |
Also Published As
Publication number | Publication date |
---|---|
WO2004077086A1 (en) | 2004-09-10 |
EP1599741A1 (en) | 2005-11-30 |
US20060145696A1 (en) | 2006-07-06 |
US7417429B2 (en) | 2008-08-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3869337B2 (en) | Magnetic resonance imaging device | |
JP4443625B2 (en) | Imaging method and apparatus by magnetic resonance | |
US9588207B2 (en) | System for reconstructing MRI images acquired in parallel | |
RU2557334C2 (en) | Dynamic contrast enhanced mr imaging with compressed sensing reconstruction | |
JP3952247B2 (en) | Nuclear magnetic resonance imaging system | |
JP6333293B2 (en) | Metal resistant MR imaging | |
CN109814058B (en) | Magnetic resonance system and method for generating a magnetic resonance image of an examination subject therein | |
US8155419B2 (en) | MRI acquisition using sense and highly undersampled fourier space sampling | |
US20020060567A1 (en) | Magnetic resonance method for forming a fast dynamic image | |
JP4610611B2 (en) | Magnetic resonance imaging device | |
JPWO2002056767A1 (en) | Parallel MR imaging using high-precision coil sensitivity map | |
JP2004275198A (en) | Apparatus and method for nuclear magnetic resonance imaging | |
US20130207652A1 (en) | System for Accelerated Magnetic Resonance Imaging Using Parallel Coils | |
JPH06189936A (en) | Magnetic resonance imaging method | |
US6564082B2 (en) | Method for incremental field-of-view-MR imaging | |
US7042217B2 (en) | Magnetic resonance method | |
JP2006519051A (en) | Movable table MRI using sub-sampling | |
US6745064B2 (en) | Magnetic resonance method for forming a fast dynamic imaging | |
US20030004410A1 (en) | Method for MR imaging with an array of RF coils | |
JP4202855B2 (en) | Magnetic resonance imaging system | |
EP1365253A2 (en) | Automatic field of view optimization in MR imaging for maximization of resolution and elimination of aliasing artifacts | |
JP4607430B2 (en) | RF coil and magnetic resonance imaging apparatus | |
JP3938717B2 (en) | Magnetic resonance imaging device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070216 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20091119 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20091126 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20100225 |
|
A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20100304 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100521 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20100930 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20110301 |