JP2023180216A - ダイナミック2d心臓mriの交互交差取得を通したプロスペクティブスライストラッキング(pst) - Google Patents
ダイナミック2d心臓mriの交互交差取得を通したプロスペクティブスライストラッキング(pst) Download PDFInfo
- Publication number
- JP2023180216A JP2023180216A JP2023074620A JP2023074620A JP2023180216A JP 2023180216 A JP2023180216 A JP 2023180216A JP 2023074620 A JP2023074620 A JP 2023074620A JP 2023074620 A JP2023074620 A JP 2023074620A JP 2023180216 A JP2023180216 A JP 2023180216A
- Authority
- JP
- Japan
- Prior art keywords
- image
- images
- plane
- heart
- heartbeat
- 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
- 230000002526 effect on cardiovascular system Effects 0.000 title 1
- 238000000034 method Methods 0.000 claims abstract description 38
- 238000006073 displacement reaction Methods 0.000 claims abstract description 25
- 238000003384 imaging method Methods 0.000 claims description 27
- 230000008569 process Effects 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 7
- 210000003484 anatomy Anatomy 0.000 claims description 4
- 239000002872 contrast media Substances 0.000 claims description 4
- 230000004941 influx Effects 0.000 claims description 2
- 230000003287 optical effect Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 7
- 238000012937 correction Methods 0.000 description 5
- 230000000747 cardiac effect Effects 0.000 description 4
- 230000010412 perfusion Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000013184 cardiac magnetic resonance imaging Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 230000002107 myocardial effect Effects 0.000 description 2
- 230000000241 respiratory effect Effects 0.000 description 2
- 230000029058 respiratory gaseous exchange Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000001627 detrimental effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000002861 ventricular Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- 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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56509—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling
-
- 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/543—Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription
-
- 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/56308—Characterization of motion or flow; Dynamic imaging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4053—Scaling of whole images or parts thereof, e.g. expanding or contracting based on super-resolution, i.e. the output image resolution being higher than the sensor resolution
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
- G06T7/0014—Biomedical image inspection using an image reference approach
- G06T7/0016—Biomedical image inspection using an image reference approach involving temporal comparison
-
- 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/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
- G01R33/4833—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices
- G01R33/4835—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices of multiple slices
-
- 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/5601—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution involving use of a contrast agent for contrast manipulation, e.g. a paramagnetic, super-paramagnetic, ferromagnetic or hyperpolarised contrast agent
-
- 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/56366—Perfusion imaging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10088—Magnetic resonance imaging [MRI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30048—Heart; Cardiac
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Theoretical Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- High Energy & Nuclear Physics (AREA)
- Signal Processing (AREA)
- Medical Informatics (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Vascular Medicine (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
【課題】心臓のMR画像を生成する。【解決手段】本方法は、基準心拍中に心臓の2D基準画像の第1のセットを収集し再構成し、第1のセットは、心臓を通る第1の基準平面の位置に沿って収集される第1の基準画像を含む、過程と、第2の心拍中に心臓の2D画像の第2のセットを収集する過程と、を含み、第2のセットを収集する過程は、第1の基準平面に沿って収集される第2のセット中の第1の画像を収集する過程と、第1の基準画像と第2のセット中の第1の画像とに基づいて、第1の基準平面の位置の画像面における心臓の第1のインプレイン変位を決定する過程と、決定した第1のインプレイン変位に基づいて、第2のセット中の第1の画像に交差して収集する第2のセット中の第2の画像に関するシフトさせた平面の位置を決定する過程と、シフトさせた平面の位置に沿って第2の画像を収集し再構成する過程と、を含む。【選択図】図4
Description
本願は、心臓のMR画像を生成する方法、これに対応するMR撮像システム、プログラムコードからなるコンピュータプログラム、コンピュータプログラムを含んだ担体に関する。
ダイナミック2D心臓MRI(磁気共鳴撮像)は、動的時系列の収集(acquisitions)間で呼吸運動によって影響を受ける。どのような2D-MR撮像でも、呼吸運動の影響は2つのカテゴリーに分類することができる:1)撮像された2D平面の中での動態間の心臓の動き(in-plane/インプレイン)と、2)収集間の撮像された2D平面の内外への心臓の動き(through-plane/スループレイン)である。インプレイン(面内)及びスループレイン(面通過)運動の影響及び視認性は、解剖学的構造及び所定のスライス位置によって変わり得るが、いずれも通常、心臓MRI(特に短軸)で使用される標準スライス方向において重大である。さらに、両分類の重大な運動は、深呼吸時に断続的に、すなわち不規則で予測不能に、出現し得る。
インプレイン運動は画像ベースの運動補正によって遡及的に補正できるが、スループレイン運動には、収集自体が動作する部分をもはや含んでおらず、遡及的な画像ベースの運動補正が不可能であり、そして、何らかの形式のプロスペクティブ(予測)補正(例えば、スライストラッキング)を必要とする、という問題がある。心臓内の同じ箇所の長いダイナミックデータウィンドウを通例で含む灌流定量化には、スループレイン運動が特に有害である。
プロスペクティブ運動補正は、個々の収集それぞれの前に別々のナビゲータ収集を使用して達成されることが多い。このナビゲータからの情報を使用して、受容ウィンドウ外のデータを排斥し、及び/又は、ナビゲータによって示される運動に従って次の収集のスライス位置をシフトさせることができる。
セグメント化3D心臓撮像の場合、これを横隔膜の1Dコーン/ペンシルビームナビゲータとすることができ、各方向の運動を推定し、続く画像収集に予測的に影響させるか(すなわち、シフト)、遡及的にデータを補正する。
ダイナミック2D撮像の場合、1D横隔膜ナビゲータによるプロスペクティブスライストラッキングが可能であるが、いくつかの問題が残っている。
1)断続的な又は特異な深呼吸イベントをカバーしない、横隔膜運動と心筋運動との間の一定のトラッキング因子の想定、あるいは、潜在的に不正確な学習フェーズの必要性。
2)1Dナビゲータと、ダイナミック撮像で通常使用される逆位又は飽和先行パルスとの間の複雑な相互作用で、実際の診断画像収集に影響を及ぼす可能性のある追加のリストアパルスの使用が必要である。
3)横隔膜フィートヘッド運動を心筋短軸スライスに適用すると、その向きは個々の解剖学的構造に依存し、フィートヘッドが必ずしもスループレイン運動に対する唯一の寄与ではない。
1)断続的な又は特異な深呼吸イベントをカバーしない、横隔膜運動と心筋運動との間の一定のトラッキング因子の想定、あるいは、潜在的に不正確な学習フェーズの必要性。
2)1Dナビゲータと、ダイナミック撮像で通常使用される逆位又は飽和先行パルスとの間の複雑な相互作用で、実際の診断画像収集に影響を及ぼす可能性のある追加のリストアパルスの使用が必要である。
3)横隔膜フィートヘッド運動を心筋短軸スライスに適用すると、その向きは個々の解剖学的構造に依存し、フィートヘッドが必ずしもスループレイン運動に対する唯一の寄与ではない。
以上のとおり、上記の問題を克服すること、発生するスループレイン運動を効果的な方法で克服することのできる心臓のMR撮像方法を提供することが必要である。
この課題は、独立形式請求項の特徴によって解決される。さらなる態様が従属形式請求項に記載されている。
第1の態様によれば、心臓のMR画像を生成する方法が提供される。この方法は、基準心拍中に心臓の2D基準画像の第1のセットを収集して再構成する過程を含む。この第1のセットは、心臓を通る第1の基準平面の位置に沿って収集された第1の基準画像からなる。さらに、心臓の2D画像の第2のセットを第2の心拍中に収集する。この第2のセットの収集は、第1の基準平面に沿って収集される第2のセット中の第1の画像を収集する過程を含む。第1の基準画像と第2のセット中の第1の画像に基づいて、第1の基準平面の位置の画像面の心臓について第1のインプレイン変位を決定する。そして、第1の基準画像と第2のセット中の第1の画像とから決定された第1のインプレイン変位に基づいて、第2のセット中の第1の画像に交差して収集する第2のセット中の第2の画像について、シフトさせた平面の位置を決定する。さらに、シフトさせた平面の位置に沿った第2の画像を収集して再構成する。2D画像の第1のセットのすべての画像と2D画像の第2のセットのすべての画像とは、診断上関連する画像である。
第2のセット中の第2の画像はセットそれぞれの中の第1の画像と交差しているので、各セットの第1の画像におけるスループレイン運動を決定することができる。画像面が交差しているので、各セットの対応する第2の画像でインプレイン運動であればスループレイン運動になる。第2の心拍中に第2のセット中の第2の画像が決定される前に、インプレイン変位は、各種のマッチングアルゴリズムを用いて決定することができる。マッチングアルゴリズムにおいて、第1の基準画像又はこの画像中の何らかの解剖学的ランドマークが、同じ平面の位置に沿って収集される第2のセット中の第1の画像と比較される。交差する方向では、これが、第2の画像の画像の向きでスループレインの動きに相当し、したがって、第2のセット中の第2の画像に関する第2の又はシフトさせた平面の位置を決めるときに、当該変位を考慮することができる。
好ましくは、前記第2のセットがN個の画像を含んでいるとすると、第2のセット中の付加画像2~Nに関して、これら付加画像それぞれの平面の位置は、互いに平行であり、第1の基準平面の位置と交差する。この場合、付加画像それぞれの平面の位置は、第1のインプレイン変位に基づいて決定し得る。第2のセット中の付加画像はすべて、心臓を通る平行な平面であり、例えば、心臓を通る異なる位置で見た心臓の短軸像である。これらすべての画像について、対応する第1の画像から決定したインプレイン変位が、第2の心拍における画像面を設定し適合させるのに役立ち、したがって、実質的に同じ解剖学的領域が画像の第1のセットと対比した画像中に示される。
別の態様において、基準心拍中にN個の基準画像を収集し再構成することを含む方法で、第1のセットを収集し再構成することが可能である。この場合、画像の第2のセットもN個の画像からなる。両方のセットに関し、画像nの画像面の位置は、両方のセットにおいて画像n+1の画像面の位置と交差していることが可能で、nは、1からN-1までの整数である。この場合、第1のセット中の画像n及び第2のセット中の画像nで決定された対応するインプレイン変位に基づいて、第2のセット中の画像n+1に関してシフトさせた平面の位置を決定する。これにより、第2の心拍における第2のセットの各画像を、現在と直前の心拍における先行画像に基づいて補正することができる。このことは、各連続画像が直前の画像と交差する画像面を有することで可能になる。すなわち、第1の画像面が短軸像であれば後続の画像面は長軸像であり、あるいは、その逆である。
ここで、第1のセット中の第2の画像及び第2のセット中の第2の画像はいずれも、心臓の長軸又は心臓の短軸に沿うような、同じ画像面に沿った平面の位置を有することが可能である。この場合、その間の画像は交差した画像面を有する。
また、第1及び第2のセットそれぞれの中の奇数画像番号を有する画像は、心臓の長軸に沿った心臓の画像を表現する、ということが可能である。この場合、第1及び第2のセットそれぞれの中の偶数画像番号を有する画像は、心臓の短軸に沿った心臓の画像を表現する。同様に、奇数画像番号を有する画像は短軸像、偶数画像番号を有する画像は長軸像としてもよい。
好ましくは、第2のセット中の第2の画像及び付加基準画像は、実質的に心臓内の同じ解剖学的位置を有する。このことは、対応する第1の画像から決定されるインプレイン運動が、交差する平面における第2の画像に関する画像面の位置及びスループレイン運動を取得するのに役立つことから、可能になる。
また、上記方法は、第1及び第2のセットにおいて収集された画像の各々の最終版を再構成する過程を含むことができる。この場合、最終版は、それぞれ、インプレイン変位を決定するために使用した第1の画像及び第1の基準画像と比較して高い解像度を有する。各セット中の第1の画像の再構成とインプレイン変位の決定は第2の心拍中に起きなければならず、それにより第2の画像を画像の第2のセットにおいて収集することができるので、高速撮像シーケンスと再構成が必要である。一例を挙げると、撮像シーケンスには、各画像を再構成するために使用されるMR信号を1つの励起RFパルスの後に取得する撮像シーケンスを利用することができる。すなわち、撮像シーケンスは、シングルショット撮像シーケンスとすることができ、好ましくはグラジエントエコー式の撮像シーケンスとし得る。
また、マッチング処理を、第1の基準画像及び第1の画像に使用することができる。この処理において、第1の基準画像の少なくとも1つのサブセットを第1の画像の少なくとも1つのサブセットと比較し、これら2つの画像におけるインプレイン変位を決定する。このインプレイン変位は、次に続く画像のスループレイン変位に相当する。
第2の心拍中に2D画像の第2のセットを収集する過程とその一連の全過程は、継続する心拍に対しいくらでも繰り返すことができる。また、本方法の各過程は、MR造影剤の投与時に適用することができ、第1及び第2のセットの収集画像がMR造影剤の心臓への流入を表現するということが可能である。
また、相応のMR撮像システムが提供され、この撮像システムは、心臓のMR画像を生成するように構成される。この場合、撮像システムは、上述したように又は以下でさらに詳述するように作動可能な制御ユニットを含む。
また、プログラムコードを含むコンピュータプログラムが提供され、このプログラムコードは、MR撮像システムの制御ユニットによって実行されると、MR撮像システムに、上述した又は以下でさらに詳述する方法を実行させる。
また、上記コンピュータプログラムの担体が提供され、この担体は、電子信号、光信号、無線信号、及びコンピュータ可読記憶媒体のうちの1つである。
本発明の以上の及びその他の特徴及び効果は、次の図面を参照する以下の詳細な説明から明らかとなる。図中、類似の参照符号は類似の要素を表す。
心臓のMR撮像において生じるスループレインの動きを最適補正して心臓のMR画像を生成することができるMRシステムの概略図。
同じ心拍中に後続の画像面の位置を先行の画像面と交差するように決定する方法の概略図。
1回の長軸画像収集後の3回の短軸画像収集の第1の概略図。
スループレイン運動を考慮するために、異なる心拍の画像をどのように収集して処理するかを示す概略図。
スループレインアーチファクトを回避しながら心臓のMR画像を決定するために必要ないくつかの過程を含むフローチャートの概略図。
以下、本発明の実施の形態について、図面を参照して詳細に説明する。以下に述べる実施の形態の説明が限定の意図をもつと捉えるべきではないのはもちろんである。本発明の範囲は、ここに説明する実施の形態又は図面によって限定されるべきではなく、それらは単なる例示である。
図面は、模式的、概要的な表現であると解釈されるべきであり、図示された要素は必ずしも原寸に比例して(正確な縮尺で)示されているわけではない。むしろ、種々の要素は、それらの機能及び汎用性が当分野で通常の知識をもつ者に明らかであるように表されている。図面に示しそしてここに説明する機能ブロック、デバイス(装置)、物理的又は機能的ユニットのコンポーネントの間の接続や連結は、間接的な接続や連結でも実施することができる。コンポーネント間の連結は、有線又は無線接続を通して確立され得る。機能ブロックは、ハードウエア、ソフトウェア、ファームウェア、又はそれらの組み合わせで実施することができる。
図1は、分極場B0を生成する磁石10を備えたMRシステム1の概略図を示す。テーブル11上に横たわる被検体12をMRシステム1の中央へ移動させ、RF励起後のMR信号が受信コイル2によって検出される。受信コイル2は個々のコイルセクションからなり、コイルセクションのそれぞれは、対応する検出チャネル3と結びつけられている。RFパルス及び磁場勾配を印加することにより、被検体12において、特に受信コイル2内に位置する部位(本例では心臓)において、核スピンが励起され、位置コード化され、リラクゼーションによって誘導される電流を検出することができる。MR画像を生成する方法と、RFパルスのシーケンス及び磁場勾配のシーケンスを用いてMR信号を検出する方法とは、当技術分野において知られているので詳細な説明は省略する。
MRシステムは、MRシステムを制御するために使用される制御モジュール13を含む。制御モジュール13は、磁場勾配を制御しスイッチングする勾配制御ユニット14と、撮像シーケンスのためにRFパルスを制御し生成するRF制御ユニット15と、を含む。印加RFパルス及び磁場勾配のシーケンスを制御して、勾配制御ユニット14及びRF制御ユニット15を制御する、画像シーケンス制御ユニット16が設けられる。メモリ17において、MRシステムの作動に必要なコンピュータプログラムとMR画像を生成するために必要な撮像シーケンスとを、生成されたMR画像と共に記憶することができる。生成されたMR画像はディスプレイ18に表示することができ、入力ユニット19が提供され、MRシステムのユーザによってMRシステムの機能を制御するために使用される。処理ユニット20は、図1に示す各機能ユニットの作動を管理し、メモリ17に記憶された命令を実行可能な1つ以上のプロセッサを含む。メモリには、処理ユニット20によって実行されるプログラムコードが入っている。処理ユニットは、検出された画像に基づいて、MR画像を再構成することができる。
以下の開示において、運動が起こっていても心臓に関して同じ画像面を収集できるようにスループレイン運動を考慮することを可能にする方法について説明する。特に、制御モジュール13は、後述のとおりに作動するように構成される。
時間的に互いに前後した同じ向きの平行平面のダイナミック2D画像スキャンを実行する代わりに、短軸像のような対象の向きの収集と、この対象の向きと交差する診断的に関連した画像の向きの収集とを、交互に行うことが提案される。
図2は、連続した交差潅流(パーフルージョン)スライスの図を示す。画像セットの第1の画像の例として、画像22は、長軸に沿った2つ又は3つの心室、すなわち長軸像であり、これに短軸像を示す画像23が続き、そして、長軸スライス位置で4つの心室像の画像24及び短軸画像25が続き、これらは、同じ心拍中の1つの時間ウィンドウで互いに前後して収集される。各画像において、次に続く画像の平面の位置が示されており、これは、画像22の場合、平面の位置が線31で表され、画像23の画像面を示している、ということである。同様にして画像23の場合、第3の画像24の平面の位置を表す平面位置(線)32が示されている。画像24では平面位置33が示され、画像25では平面位置34が示され、これは、画像22で示された同様の平面で、次の心拍で収集される可能性をもつ。
このような処理で、後続の交差した向きの画像を再構成することにより、それを先行する心拍の同じ向きのいくつかの任意の基準と整合させることにより、そして、それにより得られたシフトを、次の交差収集にプロスペクティブに(前向きに/予測的に)適用することにより、スライス位置に関するプロスペクティブスライス補正を、後続の心拍について得ることができる。これについては以下に図4を参照してより詳細に説明する。図2において、第1の画像中の2つの矢印28,29は、後続の短軸スライス(画像面31、画像23のスライス位置を選択された基準にアジャストするために使用される)に交差する運動の投影を表す。この方法により、画像23に表示される短軸像は、検査全体を通して同じ解剖学的位置で収集され、すべての残りのインプレイン運動を遡及的に補正することができる。連続したスライス収集の間の時間は短いが、前のスライスの予備的で低解像度の再構成にかかる時間と、続くスライス位置を調整するためのその基準に対するマッチングとが、通常採用される飽和回復(SR)準備によって可能とされ、約50ms以上の次のデータ収集の開始前の時間を残す。交差の向きは、いくつかの異なる方法で適用することができる。
第1の例が図3に示され、他の例が図4に示されている。両方の例において、心拍ごとにN個の2Dスライス画像がダイナミック潅流実験において収集されている。図3は、2つの心室又は3つの心室又は4つの心室などの長軸像である第1の画像41を示し、これにN-1個の短軸スライスが続き、これらは、各心拍において同じ第1の長軸スライスに従ってプロスペクティブにシフトさせてある。このことは、画像42,43,44が収集される前に、画像41が、先行する心拍で取得された同じ画像と比較され、これら2つの画像の比較に基づいており、そしてインプレイン運動がこの画像面の位置に関して決定可能であることを意味する。線47によって示される次の平面が交差しており、スライス位置は対応して適合させることができる。
図4は、異なる画像がどのように収集され、後処理がどのように行われるかを説明する一種のフローチャートを示している。第1の心拍において、画像51~54が収集される。このときの画像スライス位置56~59は、後続の画像スライスが互いに交差するようになっている。図示の例で、画像51は長軸像であり、画像52は短軸像であり、画像53は長軸像であり、画像54は短軸像である。つまり、交差位置は交互方式で収集される。この第1の過程の後、第2の過程において、これら画像は、最初の位置の第1の心拍の基準画像として再構成もされる。第1の画像セット中の画像51~54は、基準画像セットの役割をもつ。
第3の過程において、すべての残りの心拍について、基準画像セットと同じ順番で同じスライス位置が使用され、この第3の過程の最初に、過程3-aにおいて、画像61が、患者の呼吸運動を考慮して、直前の心拍における前回の収集に対して潜在的にシフトさせたスライス位置で収集され、好ましくは低解像度で再構成される。次に続く画像62~64が収集される前に、過程3-bにおいて矢印80によって示すように、画像61は、画像51とマッチングさせる。このマッチングは、完全な画像51と画像61、又はこれら画像のサブセットに基づき得る。次の過程3-cにおいて、後続のスライス位置に関するスループレイン運動に相当するインプレイン運動変位を計算し、そして過程3-dにおいて、次に続くスライス位置、すなわち画像62の生成に使用する平面位置65、をシフトさせて過程3-b及び過程3-cの基準位置とマッチングさせる。以上に説明したのと同じ過程3-a~dが、後続の画像に適用される。すなわち、画像62を画像52と比較してこれらの画像に示されたインプレイン運動を決定し、この運動は、画像63に関するスループレイン運動に対応するので、その結果、画像63が収集される前に画像面の位置を適合させることができる。図4において画像62から画像64の間に示す各矢印についても、過程a~dが繰り返される。同じ処理を第3の心拍の画像スライスに適用することができ、このときには第2の心拍の画像が基準画像の役割をもつ。代案として、第1の心拍の画像が基準画像の役割をもつ。
図3に示す例に適用すると、次のようになり得る。画像51~54が図4に示すとおり収集されるか、又は、画像51のみが収集され、第2の心拍において、図3に示す画像41が決定され、異なる短軸像42、短軸像43、又は短軸像44が図3に象徴されるように収集される前に過程a~dを実行することができる。
さらなる過程として、第3の過程の過程a~dの後、第4の過程において、全収集サイクルの終了時に、それぞれのスライス及び心拍のすべてて全画像のフル解像度版を再構成することが可能である。最終的な診断画像再構成の後回しは、画像62を収集する前に必要となるフル解像度再構成及びマッチングがリアルタイムでは時間がかかりすぎることから、有益である。
横隔膜サロゲートからではなく、心臓自体の画像から2D画像スライスに交差して運動が推定され、1Dナビゲータで可能である精度よりも正確にスループレイン運動を決定することができる。先行する完全に診断用の収集がナビゲーションの基礎であるから、追加のデータ収集が必要ない。このことは、さらに、実際の灌流収集時のナビゲータリストアパルスの潜在的干渉を回避し、横隔膜運動と心臓運動との間の関係に関する学習された又は先験的な仮定を含まない。
図5は、上述した過程の一部をまとめたものである。第1の過程S81において、心臓の2D基準画像の第1のセットが収集され、画像51~54によって示したように再構成される。このとき、少なくとも第1の基準画像51が提供される。過程S82において、第2のセットが収集され、図3に示した画像41~44又は図4に示した画像61~64のいずれかである。第2のセットを収集するこの過程は、過程S83を含み、過程S83において、第1の基準平面に沿って第2のセット中の第1の画像が収集され、つまり、基準セット中の第1の画像と同じ画像面に沿って第2の心拍中の第1の画像が収集される。過程S84において、第1の基準画像及び第2のセット中の第1の画像に基づいて、すなわち、画像61及び画像51又はこれら画像のサブセットに基づいて、第1のインプレイン変位が決定される。この第1のインプレイン変位に基づいて、第2のセット中の第1の画像と交差して収集されるべき第2のセット中の第2の画像に関し、シフトさせた平面の位置を決定することが可能である。これが過程S85であり、次の過程S86において、シフトさせた平面の位置に沿って第2の画像を収集し再構成することができる。2D画像の第1のセット中のすべての画像と2D画像の第2のセット中のすべての画像は、診断的に関連する画像である。
最後に、上述した本発明は、適用分野の1つとして灌流ベースの画像に有用な、同じスライス位置で心臓の画像を生成する効果的な方法を提供する。
Claims (15)
- 心臓のMR画像を生成する方法であって、
基準心拍中に心臓の2D基準画像の第1のセットを収集し再構成し、前記第1のセットは、心臓を通る第1の基準平面の位置に沿って収集される第1の基準画像を含む、過程と、
第2の心拍中に心臓の2D画像の第2のセットを収集する過程と、を含み、
前記第2のセットを収集する過程は、
前記第1の基準平面に沿って収集される前記第2のセット中の第1の画像を収集する過程と、
前記第1の基準画像と前記第2のセット中の前記第1の画像とに基づいて、前記第1の基準平面の位置の画像面における心臓の第1のインプレイン変位を決定する過程と、
決定した前記第1のインプレイン変位に基づいて、前記第2のセット中の前記第1の画像に交差して収集する前記第2のセット中の第2の画像に関するシフトさせた平面の位置を決定する過程と、
前記シフトさせた平面の位置に沿って前記第2の画像を収集し再構成する過程と、を含み、
2D画像の前記第1のセットのすべての画像及び2D画像の前記第2のセットのすべての画像は、診断的に関連する画像である、方法。 - 前記第2のセットがN個の画像を含み、
前記第2のセット中の付加画像(2~N)は、これら付加画像それぞれの平面の位置が互いに平行であると共に前記第1の基準平面の位置と交差し、
前記付加画像それぞれの前記平面の位置は、前記第1のインプレイン変位に基づいて決定される、請求項1に記載の方法。 - 前記第1のセットを収集し再構成する過程は、前記基準心拍中にN個の基準画像を収集し再構成する過程を含み、
前記第2のセットがN個の画像を含んでおり、その画像nの画像面の位置が前記第1のセット及び前記第2のセット中の画像n+1の画像面の位置に交差し、この場合のnは1からN-1までの整数であり、
前記第2のセット中の画像n+1に関し決定される前記シフトさせた平面の位置は、前記第1のセット中の画像n及び前記第2のセット中の画像nで決定される、対応するインプレイン変位に基づいて、決定される、請求項1に記載の方法。 - 前記第1のセット及び前記第2のセット中のすべての前記第2の画像が、心臓の長軸に沿った平面の位置又は心臓の短軸に沿った平面の位置を有する、請求項3に記載の方法。
- 前記第1のセット及び前記第2のセットそれぞれの中の奇数画像番号をもつ画像が心臓の長軸に沿った心臓の画像を表現し、前記第1のセット及び前記第2のセットそれぞれの中の偶数画像番号をもつ画像が心臓の短軸に沿った心臓の画像を表現する、請求項4に記載の方法。
- 前記第2のセット中の前記第2の画像及び別の前記基準画像は、心臓内で同じ解剖学的位置を実質的に有する、請求項1~5のいずれか1項に記載の方法。
- 前記第1のセット及び前記第2のセット中の収集された画像それぞれの最終版を再構成する過程をさらに含み、
前記最終版は、それぞれ、前記第1の画像及び前記第1の基準画像と比較して高い解像度を有する、請求項1~6のいずれか1項に記載の方法。 - 前記第1のセット及び前記第2のセット中の各画像は、シングルショットグラジエントエコー撮像シーケンスで収集される、請求項1~7のいずれか1項に記載の方法。
- 前記第1の基準画像及び前記第1の画像に関しマッチング処理を使用し、
このマッチング処理において、前記第1の基準画像の少なくともサブセットと前記第1の画像の少なくともサブセットとを比較してインプレイン変位を決定する、請求項1~8のいずれか1項に記載の方法。 - 前記第2の心拍中に心臓の2D画像の前記第2のセットを収集する過程とその一連の全過程は、継続する心拍に対しいくらでも繰り返される、請求項1~9のいずれか1項に記載の方法。
- 当該方法が、少なくとも部分的にMR造影剤の投与時に適用され、収集される画像データが心臓への造影剤の流入を表現する、請求項1~10のいずれか1項に記載の方法。
- 制御モジュールを含み、心臓のMR画像を生成するように構成されたMR撮像システムであって、
前記制御モジュールは、MR画像を生成するために、
基準心拍中に心臓の2D基準画像の第1のセットを収集し再構成し、
第2の心拍中に心臓の2D画像の第2のセットを収集する、ように構成され、
前記第1のセットは、心臓を通る第1の基準平面の位置に沿って収集される第1の基準画像と、前記第1の基準平面の位置と交差する第2の基準平面の位置で収集される少なくとも1つの別の基準画像と、を含み、前記別の基準画像が心臓の所定の解剖学的領域を表現し、
前記第2の心拍中に心臓の2D画像の前記第2のセットを収集するときに、前記制御モジュールは、
前記第1の基準平面に沿って収集される前記第2のセット中の第1の画像を収集し、
前記第1の基準画像と前記第2のセット中の前記第1の画像とに基づいて、前記第1の基準平面の位置の画像面における心臓の第1のインプレイン変位を決定し、
決定した前記第1のインプレイン変位に基づいて、心臓の前記第2の基準平面の位置と平行に収集する前記第2のセット中の第2の画像に関するシフトさせた平面の位置を決定し、
前記シフトさせた平面の位置に沿って前記第2の画像を収集し再構成する、ように構成され、
2D画像の前記第1のセットのすべての画像及び2D画像の前記第2のセットのすべての画像は、診断的に関連する画像である、MR撮像システム。 - 請求項2~11のいずれか1項に記載の方法を実行するようにさらに構成されている、請求項12に記載のMR撮像システム。
- プログラムコードを含むコンピュータプログラムであって、前記プログラムコードがMR撮像システムの制御ユニットによって実行されると、該MR撮像システムが請求項1~11のいずれか1項に記載の方法を実行する、コンピュータプログラム。
- 請求項14に記載のコンピュータプログラムを含む担体であって、電子信号、光信号、無線信号、及びコンピュータ可読記憶媒体のうちの1つである、担体。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP22177768 | 2022-06-08 | ||
EP22177768.3A EP4290262A1 (en) | 2022-06-08 | 2022-06-08 | Prospective slice tracking through interleaved perpendicular acquisitions of dynamic 2d cardiovascular mri |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2023180216A true JP2023180216A (ja) | 2023-12-20 |
Family
ID=81984735
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2023074620A Pending JP2023180216A (ja) | 2022-06-08 | 2023-04-28 | ダイナミック2d心臓mriの交互交差取得を通したプロスペクティブスライストラッキング(pst) |
Country Status (4)
Country | Link |
---|---|
US (1) | US20230401767A1 (ja) |
EP (1) | EP4290262A1 (ja) |
JP (1) | JP2023180216A (ja) |
CN (1) | CN117197184A (ja) |
-
2022
- 2022-06-08 EP EP22177768.3A patent/EP4290262A1/en active Pending
-
2023
- 2023-04-28 JP JP2023074620A patent/JP2023180216A/ja active Pending
- 2023-06-06 US US18/329,655 patent/US20230401767A1/en active Pending
- 2023-06-06 CN CN202310665396.0A patent/CN117197184A/zh active Pending
Also Published As
Publication number | Publication date |
---|---|
EP4290262A1 (en) | 2023-12-13 |
US20230401767A1 (en) | 2023-12-14 |
CN117197184A (zh) | 2023-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5565777A (en) | Method/apparatus for NMR imaging using an imaging scheme sensitive to inhomogeneity and a scheme insensitive to inhomogeneity in a single imaging step | |
US8427153B2 (en) | Method for motion correction in magnetic resonance imaging using radio frequency coil arrays | |
US5427101A (en) | Diminishing variance process for real-time reduction of motion artifacts in MRI | |
US9804246B2 (en) | Magnetic resonance imaging apparatus | |
US9684979B2 (en) | MRI 3D cine imaging based on intersecting source and anchor slice data | |
KR101664138B1 (ko) | 멀티 슬라이스 자기 공명 영상을 사용하는 조직 특성 감별 시스템 및 방법 | |
CN103792502A (zh) | 用于成像中时序的自动系统 | |
WO2012060252A1 (ja) | 磁気共鳴イメージング装置及び磁気共鳴イメージング方法 | |
US10165960B2 (en) | Magnetic resonance 2D navigator technique | |
US6307369B1 (en) | Autocorrection of 3D MR images for motion artifacts | |
CN106456046A (zh) | 经改进的多时相动态对比增强磁共振成像的方法 | |
JP2001309903A (ja) | Mr投影画像の自動修正 | |
US20070285090A1 (en) | Phase cycling method and magnetic resonance imaging apparatus | |
US20180217218A1 (en) | Image reconstruction for mri using multiplexed sensitivity encoding | |
US11454692B2 (en) | Method of performing magnetic resonance imaging and a magnetic resonance apparatus | |
JP2003515404A (ja) | 素材の変形及びひずみを検査するための高速3次元磁気共鳴タッギング | |
JP2005040416A (ja) | 磁気共鳴イメージング装置 | |
JP4912808B2 (ja) | 磁気共鳴イメージング装置 | |
JP2641486B2 (ja) | Nmrイメージング装置 | |
JP4349647B2 (ja) | 磁気共鳴イメージング装置 | |
JP2023180216A (ja) | ダイナミック2d心臓mriの交互交差取得を通したプロスペクティブスライストラッキング(pst) | |
US20220338829A1 (en) | Automatic Determination of a Motion Parameter of the Heart | |
JP6356412B2 (ja) | 磁気共鳴装置およびプログラム | |
CN113466768B (zh) | 磁共振成像方法及磁共振成像系统 | |
JP4519827B2 (ja) | 磁気共鳴映像化装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20231027 |