JP7059189B2 - 定常状態磁気共鳴フィンガープリンティング - Google Patents

定常状態磁気共鳴フィンガープリンティング Download PDF

Info

Publication number
JP7059189B2
JP7059189B2 JP2018541247A JP2018541247A JP7059189B2 JP 7059189 B2 JP7059189 B2 JP 7059189B2 JP 2018541247 A JP2018541247 A JP 2018541247A JP 2018541247 A JP2018541247 A JP 2018541247A JP 7059189 B2 JP7059189 B2 JP 7059189B2
Authority
JP
Japan
Prior art keywords
magnetic resonance
imaging system
resonance imaging
space
pulse sequence
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2018541247A
Other languages
English (en)
Other versions
JP2019512282A5 (ja
JP7059189B6 (ja
JP2019512282A (ja
Inventor
トーマス エリック アムソー
ピーター コーケン
カルステン ソマー
マリヤ イワノワ ドネヴァ
ペーター ボルネート
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
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 Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of JP2019512282A publication Critical patent/JP2019512282A/ja
Publication of JP2019512282A5 publication Critical patent/JP2019512282A5/ja
Application granted granted Critical
Publication of JP7059189B2 publication Critical patent/JP7059189B2/ja
Publication of JP7059189B6 publication Critical patent/JP7059189B6/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01R33/5608Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
    • 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/46NMR spectroscopy
    • G01R33/4625Processing of acquired signals, e.g. elimination of phase errors, baseline fitting, chemometric analysis
    • 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/4818MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/482MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a Cartesian trajectory
    • 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/4818MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/4824MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory
    • 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/4828Resolving the MR signals of different chemical species, e.g. water-fat imaging
    • 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
    • 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
    • G01R33/561Image 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/5611Parallel 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
    • 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
    • G01R33/561Image 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/5613Generating steady state signals, e.g. low flip angle sequences [FLASH]
    • 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/4818MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space

Landscapes

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

Description

本発明は、磁気共鳴イメージング、特に、磁気共鳴フィンガープリンティングに関する。
磁気共鳴フィンガープリンティング(MRF)は、異なる材料又は組織からの信号が、測定された磁気共鳴(MR)信号に独自の寄与を有するように、時間的に分散したいくつかのRFパルスが印加される新しい技術である。設定された又は固定された数の物質からの事前計算された信号寄与の限られたディクショナリが、測定されたMR信号と比較され、単一のボクセル内で組成が決定され得る。例えば、ボクセルが水、脂肪、及び筋組織しか含まないことが分かっている場合、これらの3つの物質からの寄与しか考慮する必要がなく、ボクセルの組成を正確に決定するのにごくわずかのRFパルスしか必要とされない。
磁気共鳴フィンガープリンティング技術は、Ma等による学術論文「Magnetic Resonance Fingerprinting」、Nature、495巻、187から193ページ、doi:10.1038/Nature11971で紹介された。磁気フィンガープリンティング技術は、米国特許出願第2013/0271132A1号及び米国特許出願第2013/0265047A1号にも記載されている。米国特許出願第2015/0346300号は、同時マルチボリューム取得を用いて適用される磁気フィンガープリンティング(MRF)法を開示している。
本発明は、独立請求項に磁気共鳴イメージングシステム、コンピュータプログラム製品及び方法を提供する。実施形態は従属請求項に記載される。
Ma等によるNatureの論文は、磁気共鳴フィンガープリンティングの基本的な考え方と、本明細書で「定常状態磁気共鳴フィンガープリンティングディクショナリ」又は「過渡共鳴フィンガープリンティングディクショナリ」と呼ばれるディクショナリなどのこの技術を説明するために使用される用語とを紹介している。「ディクショナリ」という用語は、ただ磁気共鳴フィンガープリンティングディクショナリだけを指すと理解されたい。
磁気共鳴(MR)データは、本明細書では、磁気共鳴イメージングスキャン中に磁気共鳴装置のアンテナを使用して原子スピンによって発せられた無線周波数信号の記録された測定結果として定義される。磁気共鳴データは、医療画像データの一例である。磁気共鳴イメージング(MRI)画像は、本明細書において、磁気共鳴イメージングデータ内に含まれる解剖学的データの復元された2次元又は3次元視覚化であるとして定義される。この視覚化は、コンピュータを使用して実行され得る。
一態様では、本発明は、測定ゾーン内の被検体から磁気共鳴データを取得するための磁気共鳴イメージングシステムを提供する。磁気共鳴イメージングシステムは、磁気共鳴イメージングシステムを制御するプロセッサを含む。磁気共鳴イメージングシステムは、マシン実行可能命令を記憶するメモリを更に含む。メモリは、パルスシーケンスコマンドを更に記憶するためのものである。パルスシーケンスコマンドは、磁気共鳴フィンガープリンティングプロトコルに従って磁気共鳴データを取得するように磁気共鳴イメージングシステムを制御するように構成される。本明細書で使用されるパルスシーケンスコマンドは、磁気共鳴イメージングプロトコルに従って磁気共鳴イメージングシステムを制御する命令又はコマンドに変換される命令又はデータのいずれかを含む。パルスシーケンスコマンドは、高周波パルス列を生成するように磁気共鳴イメージングシステムを制御するように構成される。パルスシーケンスコマンドは、磁気共鳴データを複数のk空間トレースとして取得するように磁気共鳴イメージングシステムを制御するように構成される。
磁気共鳴イメージングにおいて、データは、所謂k空間でサンプリングされる。k空間トレースは、サンプリングされるk空間を通る経路である。デカルト空間でサンプリングされる場合、k空間トレースは、サンプリングされるk空間の所謂ラインを指す。パルスシーケンスコマンドは、複数のk空間トレースの各々の取得のために繰り返されるべき高周波パルス列を制御するように構成される。
マシン実行可能命令の実行によって、プロセッサは、パルスシーケンスコマンドで磁気共鳴イメージングシステムを制御することによって、磁気共鳴データの複数のk空間トレースを順次取得する。マシン実行可能命令の実行によって、更に、プロセッサは、複数のk空間トレースのうちの初期に取得される所定の数のk空間トレースが取得された後に取得されるk空間トレースに対する所定の物質のセットの各々の存在量を計算する。所定の物質のセットの各々の存在量は、磁気共鳴データを定常状態磁気共鳴フィンガープリンティングディクショナリと比較することによって決定される。定常状態磁気共鳴フィンガープリンティングディクショナリは、所定の物質のセットに対してRFパルス列に応答して計算された磁気共鳴信号のリストを含む。
この実施形態では、初期に取得される所定の数のk空間トレースが取得された後にk空間トレースと定常状態磁気共鳴フィンガープリンティングディクショナリとを比較することにより、複数のk空間トレースの取得をより迅速に実行することが可能になる。RFパルス列が1回実行された後、測定ゾーン内のビンに残留磁気が存在する。所定の数のk空間トレースが取得された後、磁化は最終的に平衡点に達する。この点の後、定常状態磁気共鳴フィンガープリンティングディクショナリが使用される。これにより、例えば、磁気共鳴データを連続的に取得することが可能になる。これにより、更に、磁気共鳴フィンガープリンティングを使用してより複雑な磁気共鳴技術を実施することが可能になる。
一例として、SENSEなどの並列イメージング技術では、磁気共鳴フィンガープリンティングのためのRFパルス列が使用される。これを行う難しさは、より多くのデータが取得され得る前に、残留磁気が長い遅延を必要とすることである。一定の数のRFパルス列が実行された後でのみ使用される定常状態磁気共鳴フィンガープリンティングディクショナリを使用することによって、磁気共鳴データを使用して、定常状態磁気共鳴フィンガープリンティングディクショナリにおける所定の物質の存在量を正確に計算することが可能になる。
別の実施形態では、マシン実行可能命令の実行によって、プロセッサは、初期に取得した所定の数のk空間トレースを廃棄する。例えば、スピンを定常状態磁化にするために、パルスシーケンスが初期に所定の回数実行される。これは、単一の磁気共鳴フィンガープリンティングディクショナリしか必要でなく、この単一の磁気共鳴フィンガープリンティングディクショナリを使用して、磁気共鳴データの取得が速められるという技術的な効果を有する。1つの例では、最初の所定の数の取得が廃棄され、この後、使用されるべき真の又は実際のデータが取得される。
別の実施形態では、パルスシーケンス命令は、連続的に繰り返されるべきRFパルス列を制御するように構成される。これは、磁気共鳴データができるだけ迅速に取得されるという技術的な効果を有する。
別の実施形態では、パルスシーケンス命令は、各RFパルス列の間の最大遅延でRFパルス列を繰り返すように磁気共鳴イメージングシステムを制御するように構成される。最大遅延は、下記の5秒未満、1秒未満、0.5秒未満、0.1秒未満、0.05秒未満、及び0.00秒のうちのいずれかである。この実施形態は、取得の間の遅延を最小にするという技術的な効果を有する。
別の実施形態では、磁気共鳴イメージングシステムは、メイン磁場をもつ磁石を含む。メイン磁場は、多くの場合、所謂B0場と呼ばれる。メイン磁場は、測定ゾーン内に平均の磁場の大きさを有する。マシン実行可能命令の実行によって、プロセッサが、各RFパルス列を、下記の平均の磁場の大きさにおける所定の物質のセットの最大T1時間の5倍、平均の磁場の大きさにおける所定の物質のセットの最大T1時間の1倍、平均の磁場の大きさにおける所定の物質のセットの最大T1時間の0.5倍、及び平均の磁場の大きさにおける所定の物質のセットの最大T1時間の0.1倍のうちのいずれかよりも速く繰り返す。定常状態磁気共鳴フィンガープリンティングディクショナリにおける所定の物質のセットの各々は、特定のT1緩和時間を有する。これらの物質の全ての最大T1値は、更に、メイン磁場の平均強度に依存する。上述で指定した時間の各々は、定常状態磁気共鳴フィンガープリンティングディクショナリが使用されない場合よりも磁気共鳴データを迅速に取得することを可能にする。
別の実施形態では、マシン実行可能命令の実行によって、更に、プロセッサは、定常状態磁気共鳴フィンガープリンティングディクショナリを計算する。定常状態磁気共鳴フィンガープリンティングディクショナリは、定常状態磁気共鳴フィンガープリンティングディクショナリを計算する標準方法のうちの任意の1つを使用して計算される。例えば、ディクショナリは、ブロッホ方程式を解くことによって又は拡張位相グラフ計算を実行することによって計算される。
定常状態磁気共鳴フィンガープリンティングディクショナリは、例えば、RFパルス列を繰り返すことによって計算され、計算は定常状態の結果に収束する。磁気共鳴フィンガープリンティングでは、ディクショナリは特定のRFパルスに対して計算される。同じ技術が、定常状態ディクショナリを計算するために使用され得る。単一RFパルスが複数回繰り返され、同じ計算を使用してスピンの応答を計算する。計算が所定の回数繰り返された後、RFパルス列への応答は定常状態に達する。定常状態磁気共鳴フィンガープリンティングディクショナリをこのように計算するプロセスは、所定の数を決定するためにも使用される。
別の実施形態では、マシン実行可能命令の実行によって、更に、プロセッサは、離散的なボクセルの各々に対してブロッホ方程式を用いて、所定の物質の各々を単一スピンとしてモデリングすることによって定常状態磁気共鳴フィンガープリンティングディクショナリを計算する。例えば、離散的なボクセルの各々において、仮想スピンが、ブロッホ方程式と、パルスシーケンス命令を使用する磁気共鳴システムのシミュレーションとを使用してモデリングされる。そのとき、サンプリング時間の各々における計算された磁気共鳴データは、モデリングされた特定タイプのスピンの磁気共鳴フィンガープリンティングディクショナリである。これは、測定ゾーンが単に単一ボクセルに分割される場合に特に適切に機能する。それは、空間符号化のための勾配磁場がない場合にも当てはまる。例えば、磁気共鳴システムは、試料の化学分析を行うための所謂NMRシステムである。
別の実施形態では、マシン実行可能命令の実行によって、更に、プロセッサは、所定の物質のセットの各々へのRFパルス列の繰り返し印加をモデリングすることによって定常状態磁気共鳴フィンガープリンティングディクショナリを計算する。
別の実施形態では、マシン実行可能命令の実行によって、更に、プロセッサは、定常状態磁気共鳴フィンガープリンティングディクショナリの計算の間に所定の数を決定する。所定の数は、定常状態磁気共鳴フィンガープリンティングディクショナリの計算の間に生成されたシミュレート磁気共鳴データの収束を決定する収束基準を使用して決定される。収束基準は、例えば、複数個のデータの距離を測定しフィッティングを計算するために使用される最小二乗ルーチン又は他のアルゴリズムである。
別の実施形態では、パルスシーケンスコマンドは、下記のk空間において複数のk空間トレースを回転させること、ラジアルサンプリングを使用して複数のk空間トレースを取得すること、デカルトサンプリングを使用して複数のk空間トレースを取得すること、及び非デカルトサンプリングを使用して複数のk空間トレースを取得することのうちのいずれかを実行するように磁気共鳴イメージングシステムを制御するように構成される。
別の実施形態では、パルスシーケンス命令は、PROPELLER磁気共鳴イメージングプロトコルに従って磁気共鳴データを取得するように磁気共鳴イメージングシステムを制御するように構成される。PROPELLER磁気共鳴イメージングプロトコルは、空間におけるk空間ライン又はトレースを、それらを回転させることによって取得する。定常状態磁気共鳴フィンガープリンティングディクショナリの使用によって、磁気共鳴データをより迅速に取得し、それによって、被検体の動きの影響を低減することが可能になる。PROPELLER磁気共鳴イメージングプロトコルの使用は、被検体の動きを補正するためにも利用される。2つを組み合わせると、磁気共鳴データが取得されるとき被検体の動きに影響されにくいより迅速に取得された磁気共鳴データがもたらされる。
別の実施形態では、パルスシーケンス命令は、複数のk空間トレースをk空間におけるスパイラルとして取得するように磁気共鳴イメージングシステムを制御するように構成される。
別の実施形態では、パルスシーケンス命令は、複数のk空間トレースをデカルトk空間におけるラインとして取得するように磁気共鳴イメージングシステムを制御するように構成される。
別の実施形態では、パルスシーケンス命令は、磁気共鳴データを、下記の並列イメージング磁気共鳴イメージングプロトコル、SENSE磁気共鳴イメージングプロトコル、及びGRAPPA磁気共鳴イメージングプロトコルのうちのいずれかに従って取得するように磁気共鳴イメージングシステムを制御するように構成される。
別の実施形態では、パルスシーケンス命令は、SENSE磁気共鳴イメージングプロトコルに従って磁気共鳴データを取得するように磁気共鳴イメージングシステムを制御するように構成される。
別の態様では、本発明は、測定ゾーン内の被検体から磁気共鳴データを取得するように磁気共鳴イメージングシステムを制御するプロセッサによる実行のためのマシン実行可能命令を含むコンピュータプログラム製品を提供する。マシン実行可能命令によって、プロセッサは、磁気共鳴フィンガープリンティングプロトコルに従って磁気共鳴データを取得するように磁気共鳴イメージングシステムを制御するためのパルスシーケンスコマンドで磁気共鳴イメージングシステムを制御することによって、磁気共鳴データの複数のk空間トレースを順次取得する。パルスシーケンスコマンドは、RFパルス列を生成するように磁気共鳴イメージングシステムを制御するように構成される。パルスシーケンスコマンドは、磁気共鳴データを複数のk空間トレースとして取得するように磁気共鳴イメージングシステムを制御するように構成される。パルスシーケンスコマンドは、複数のk空間トレースの各々の取得のために繰り返されるべきRFパルス列を制御するように構成される。マシン実行可能命令の実行によって、更に、プロセッサは、複数のk空間トレースのうちの初期に取得される所定の数のk空間トレースが取得された後に取得されるk空間トレースに対する所定の物質のセットの各々の存在量を計算する。所定の物質のセットの各々の存在量は、磁気共鳴データを定常状態磁気共鳴フィンガープリンティングディクショナリと比較することによって決定される。定常状態磁気共鳴フィンガープリンティングディクショナリは、所定の物質のセットに対してRFパルス列に応答して計算された磁気共鳴信号のリストを含む。
別の態様では、本発明は、更に、測定ゾーン内の被検体から磁気共鳴データを取得するための磁気共鳴イメージングシステムを動作させる方法を提供する。この方法は、磁気共鳴フィンガープリンティングプロトコルに従って磁気共鳴データを取得するように磁気共鳴イメージングシステムを制御するためのパルスシーケンスコマンドで磁気共鳴イメージングシステムを制御することによって、磁気共鳴データの複数のk空間トレースを順次取得するステップを有する。パルスシーケンスコマンドは、更に、RFパルス列を生成するように磁気共鳴イメージングシステムを制御するように構成される。パルスシーケンスコマンドは、磁気共鳴データを複数のk空間トレースとして取得するように磁気共鳴イメージングシステムを制御するように構成される。パルスシーケンスコマンドは、複数のk空間トレースの各々の取得のために繰り返されるべきRFパルス列を制御するように構成される。この方法は、複数のk空間トレースのうちの初期に取得される所定の数のk空間トレースが取得された後に取得されるk空間トレースに対する所定の物質のセットの各々の存在量を計算するステップを更に有する。所定の物質のセットの各々の存在量は、磁気共鳴データを定常状態磁気共鳴フィンガープリンティングディクショナリと比較することによって決定される。定常状態磁気共鳴フィンガープリンティングディクショナリは、所定の物質のセットに対してRFパルス列に応答して計算された磁気共鳴信号のリストを含む。
本発明の前述の実施形態のうちの1つ又は複数は、組み合わされた実施形態が相互に排他的でない限り組み合わせることができることが理解されよう。
当業者には理解されるように、本発明の態様は、装置、方法又はコンピュータプログラムプロダクトとして具体化され得る。従って、本発明の態様は、全面的にハードウェア実施形態、全面的にソフトウェア実施形態(ファームウェア、常駐ソフトウェア、マイクロコード等を含む)又は本明細書において全て一般的に「回路」、「モジュール」若しくは「システム」と称され得るソフトウェア及びハードウェア態様を組み合わせた実施形態の形態をとり得る。更に、本発明の態様は、コンピュータ可読媒体上で具現化されたコンピュータ実行可能コードを有する1つ又は複数のコンピュータ可読媒体において具体化されたコンピュータプログラムプロダクトの形態をとり得る。
1つ又は複数のコンピュータ可読媒体の任意の組み合わせが利用されてもよい。コンピュータ可読媒体は、コンピュータ可読信号媒体又はコンピュータ可読ストレージ媒体でもよい。本明細書で使用される「コンピュータ可読ストレージ媒体」は、コンピューティングデバイスのプロセッサによって実行可能な命令を保存することができる任意の有形ストレージ媒体を包含する。コンピュータ可読ストレージ媒体は、コンピュータ可読非一時的ストレージ媒体と称される場合もある。コンピュータ可読ストレージ媒体はまた、有形コンピュータ可読媒体と称される場合もある。一部の実施形態では、コンピュータ可読ストレージ媒体はまた、コンピューティングデバイスのプロセッサによってアクセスされることが可能なデータを保存可能であってもよい。コンピュータ可読ストレージ媒体の例は、フロッピー(登録商標)ディスク、磁気ハードディスクドライブ、半導体ハードディスク、フラッシュメモリ、USBサムドライブ、ランダムアクセスメモリ(RAM)、読み取り専用メモリ(ROM)、光ディスク、磁気光学ディスク、及びプロセッサのレジスタファイルを含むが、これらに限定されない。光ディスクの例は、例えば、CD-ROM、CD-RW、CD-R、DVD-ROM、DVD-RW、又はDVD-Rディスクといったコンパクトディスク(CD)及びデジタル多用途ディスク(DVD)を含む。コンピュータ可読ストレージ媒体という用語は、ネットワーク又は通信リンクを介してコンピュータデバイスによってアクセスされることが可能な様々な種類の記録媒体も指す。例えば、データは、モデムによって、インターネットによって、又はローカルエリアネットワークによって読み出されてもよい。コンピュータ可読媒体上で具現化されたコンピュータ実行可能コードは、限定されることはないが、無線、有線、光ファイバケーブル、RF等を含む任意の適切な媒体、又は上記の任意の適切な組み合わせを用いて送信されてもよい。
コンピュータ可読信号媒体は、例えばベースバンドにおいて又は搬送波の一部として内部で具体化されたコンピュータ実行可能コードを備えた伝搬データ信号を含んでもよい。このような伝搬信号は、限定されることはないが電磁気、光学的、又はそれらの任意の適切な組み合わせを含む様々な形態のいずれかをとり得る。コンピュータ可読信号媒体は、コンピュータ可読ストレージ媒体ではない及び命令実行システム、装置、若しくはデバイスによって又はそれと関連して使用するためのプログラムを通信、伝搬、若しくは輸送できる任意のコンピュータ可読媒体でもよい。
「コンピュータメモリ」又は「メモリ」は、コンピュータ可読ストレージ媒体の一例である。コンピュータメモリは、プロセッサに直接アクセス可能な任意のメモリである。「コンピュータストレージ」又は「ストレージ」は、コンピュータ可読ストレージ媒体の更なる一例である。コンピュータストレージは、任意の不揮発性コンピュータ可読ストレージ媒体である。一部の実施形態では、コンピュータストレージは、コンピュータメモリであってもよい又はその逆でもよい。
本明細書で使用される「プロセッサ」は、プログラム、マシン実行可能命令、又はコンピュータ実行可能コードを実行可能な電子コンポーネントを包含する。「プロセッサ」を含むコンピューティングデバイスへの言及は、場合により、2つ以上のプロセッサ又は処理コアを含むと解釈されるべきである。プロセッサは、例えば、マルチコアプロセッサである。プロセッサは、また、単一のコンピュータシステム内の、又は複数のコンピュータシステムの中へ分配されたプロセッサの集合体も指す。コンピュータデバイスとの用語は、各々が1つ又は複数のプロセッサを有するコンピュータデバイスの集合体又はネットワークを指してもよいと理解されるべきである。コンピュータ実行可能コードは、同一のコンピュータデバイス内の、又は複数のコンピュータデバイス間に分配された複数のプロセッサによって実行される。
コンピュータ実行可能コードは、本発明の態様をプロセッサに行わせるマシン実行可能命令又はプログラムを含んでもよい。本発明の態様に関する動作を実施するためのコンピュータ実行可能コードは、Java(登録商標)、Smalltalk(登録商標)、又はC++等のオブジェクト指向プログラミング言語及び「C」プログラミング言語又は類似のプログラミング言語等の従来の手続きプログラミング言語を含む1つ又は複数のプログラミング言語の任意の組み合わせで書かれてもよい及びマシン実行可能命令にコンパイルされてもよい。場合によっては、コンピュータ実行可能コードは、高水準言語の形態又は事前コンパイル形態でもよい及び臨機応変にマシン実行可能命令を生成するインタプリタと共に使用されてもよい。
コンピュータ実行可能コードは、完全にユーザのコンピュータ上で、部分的にユーザのコンピュータ上で、スタンドアローンソフトウェアパッケージとして、部分的にユーザのコンピュータ上で及び部分的にリモートコンピュータ上で、又は完全にリモートコンピュータ若しくはサーバ上で実行することができる。後者の場合、リモートコンピュータは、ローカルエリアネットワーク(LAN)若しくは広域ネットワーク(WAN)を含む任意の種類のネットワークを通してユーザのコンピュータに接続されてもよい、又はこの接続は外部コンピュータに対して(例えば、インターネットサービスプロバイダを使用したインターネットを通して)行われてもよい。
本発明の態様は、本発明の実施形態による方法、装置(システム)及びコンピュータプログラムプロダクトのフローチャート、図及び/又はブロック図を参照して説明される。フローチャート、図、及び/又はブロック図の各ブロック又は複数のブロックの一部は、適用できる場合、コンピュータ実行可能コードの形態のコンピュータプログラム命令によって実施され得ることが理解されよう。相互排他的でなければ、異なるフローチャート、図、及び/又はブロック図におけるブロックの組み合わせが組み合わせられてもよいことが更に理解される。これらのコンピュータプログラム命令は、コンピュータ又は他のプログラム可能データ処理装置のプロセッサを介して実行する命令がフローチャート及び/又はブロック図の1つ又は複数のブロックにおいて指定された機能/行為を実施するための手段を生じさせるようにマシンを作るために、汎用コンピュータ、特定用途コンピュータ、又は他のプログラム可能データ処理装置のプロセッサへと提供されてもよい。
これらのコンピュータプログラム命令はまた、コンピュータ可読媒体に保存された命令がフローチャート及び/又はブロック図の1つ又は複数のブロックにおいて指定された機能/行為を実施する命令を含む製品を作るように、コンピュータ、他のプログラム可能データ処理装置、又は他のデバイスにある特定の方法で機能するように命令することができるコンピュータ可読媒体に保存されてもよい。
コンピュータプログラム命令はまた、コンピュータ又は他のプログラム可能装置上で実行する命令がフローチャート及び/又はブロック図の1つ又は複数のブロックにおいて指定された機能/行為を実施するためのプロセスを提供するように、一連の動作ステップがコンピュータ、他のプログラム可能装置又は他のデバイス上で行われるようにすることにより、コンピュータ実施プロセスを生じさせるために、コンピュータ、他のプログラム可能データ処理装置、又は他のデバイス上にロードされてもよい。
本明細書で使用される「ユーザインタフェース」は、ユーザ又はオペレータがコンピュータ又はコンピュータシステムとインタラクトすることを可能にするインタフェースである。「ユーザインタフェース」は、「ヒューマンインタフェースデバイス」と称される場合もある。ユーザインタフェースは、情報若しくはデータをオペレータに提供することができる及び/又は情報若しくはデータをオペレータから受信することができる。ユーザインタフェースは、オペレータからの入力がコンピュータによって受信されることを可能にしてもよい及びコンピュータからユーザへ出力を提供してもよい。つまり、ユーザインタフェースはオペレータがコンピュータを制御する又は操作することを可能にしてもよい、及びインタフェースはコンピュータがオペレータの制御又は操作の結果を示すことを可能にしてもよい。ディスプレイ又はグラフィカルユーザインタフェース上のデータ又は情報の表示は、情報をオペレータに提供する一例である。キーボード、マウス、トラックボール、タッチパッド、指示棒、グラフィックタブレット、ジョイスティック、ゲームパッド、ウェブコム、ヘッドセット、ギアスティック、ステアリングホイール、ペダル、有線グローブ、ダンスパッド、リモコン、及び加速度計を介したデータの受信は、すべて、オペレータから情報又はデータの受信を可能にするユーザインタフェース要素の例である。
本明細書で使用される「ハードウェアインタフェース」は、コンピュータシステムのプロセッサが外部コンピューティングデバイス及び/又は装置とインタラクトする及び/又はそれを制御することを可能にするインタフェースを包含する。ハードウェアインタフェースは、プロセッサが外部コンピューティングデバイス及び/又は装置へ制御信号又は命令を送ることを可能にしてもよい。ハードウェアインタフェースはまた、プロセッサが外部コンピューティングデバイス及び/又は装置とデータを交換することを可能にしてもよい。ハードウェアインタフェースの例は、ユニバーサルシリアルバス、IEEE1394ポート、パラレルポート、IEEE1284ポート、シリアルポート、RS-232ポート、IEEE488ポート、ブルートゥース(登録商標)接続、無線LAN接続、TCP/IP接続、イーサネット(登録商標)接続、制御電圧インタフェース、MIDIインタフェース、アナログ入力インタフェース、及びデジタル入力インタフェースを含むが、これらに限定されない。
本明細書で使用される「ディスプレイ」又は「ディスプレイデバイス」は、画像又はデータを表示するために構成された出力デバイス又はユーザインタフェースを包含する。ディスプレイは、視覚、音声、及び/又は触覚データを出力してもよい。ディスプレイの例は、コンピュータモニタ、テレビスクリーン、タッチスクリーン、触覚電子ディスプレイ、点字スクリーン、陰極線管(CRT)、蓄積管、双安定ディスプレイ、電子ペーパー、ベクターディスプレイ、平面パネルディスプレイ、真空蛍光ディスプレイ(VF)、発光ダイオード(LED)ディスプレイ、エレクトロルミネッセントディスプレイ(ELD)、プラズマディスプレイパネル(PDP)、液晶ディスプレイ(LCD)、有機発光ダイオードディスプレイ(OLED)、プロジェクタ、及びヘッドマウントディスプレイを含むが、これらに限定されない。
以下で、本発明の好ましい実施形態を、単なる例として図面を参照して説明する。
磁気共鳴イメージングシステムの一例を示す図である。 図1の磁気共鳴イメージングシステムを使用する方法を示すフローチャートである。 定常状態に達する磁気共鳴信号の概念を示す図である。 複数の既知の試料によるファントムの3つの磁気共鳴フィンガープリンティング画像を示す図である。
図において似通った参照番号を付された要素は、等価な要素であるか、同じ機能を実行するかのいずれかである。先に考察された要素は、機能が等価である場合は、後の図においては必ずしも考察されない。
図1は、磁石104を含む磁気共鳴イメージングシステム100の一例を示す。磁石104は、それを貫通するボア106を有した超伝導円筒型磁石である。異なるタイプの磁石の使用も可能であり、例えば、分割円筒磁石及び所謂開放磁石の両方を使用することも可能である。分割円筒磁石は、クライオスタットが2つの部分に分割されて、磁石の等平面へのアクセスを可能にする点を除いて、標準的な円筒磁石に類似し、斯かる磁石は、例えば荷電粒子ビーム療法と共に使用される。開放磁石は、2つの磁石部分を有し、被検体を受容するのに十分な大きさの空間を間に有して一方が他方の上にあり、2つの部分の領域の配置は、ヘルムホルツコイルのものと類似する。被検体が閉じ込められる度合いが低いことから、開放磁石が普及している。円筒磁石のクライオスタットの内部には、超伝導コイルの一群が存在する。円筒磁石104のボア106内部には、磁気共鳴イメージングを行うのに十分な程、磁場が強く均一であるイメージングゾーン108が存在する。
磁石のボア106内部には、磁気共鳴データの取得のために使用されて、磁石104のイメージングゾーン108内で磁気スピンを空間的に符号化する磁場勾配コイル110のセットも存在する。磁場勾配コイル110は、磁場勾配コイル電源112に接続される。磁場勾配コイル110は代表的なものであることが意図される。一般的に、磁場勾配コイル110は、3つの直交する空間方向に空間的に符号化を行うための3つの別々のコイルセットを含む。磁場勾配電源は、電流を磁場勾配コイルに供給する。磁場勾配コイル110に供給される電流は、時間の関数として制御され、傾斜がつけられる又はパルス出力される。
イメージングゾーン108に隣接するのは、イメージングゾーン108内の磁気スピンの配向を操作するため及び同じくイメージングゾーン108内のスピンから無線伝送を受信するための高周波コイル114である。高周波アンテナは、複数のコイル素子を含む。高周波アンテナは、チャネル又はアンテナとも呼ばれる。高周波コイル114は、高周波送受信器116に接続される。高周波コイル114及び高周波送受信器116は、別個の送信及び受信コイル並びに別個の送信器及び受信器によって置き換えられてもよい。高周波コイル114及び高周波送受信器116は、代表的なものであることが理解される。高周波コイル114は、専用送信アンテナ及び専用受信アンテナも表すように意図される。同様に、送受信器116は、別個の送信器及び複数の受信器も表す。高周波コイル114は、複数の受信/送信素子も有し、高周波送受信器116は、複数の受信/送信チャネルを有する。例えば、SENSEなどの並列イメージング技術が実行される場合、高周波コイル114は複数のコイル素子を有する。
送受信器116及び磁場勾配コイル電源112は全てコンピュータシステム126のハードウェアインタフェース128に接続される。
コンピュータストレージ134はパルスシーケンスコマンド140を含むものとして示される。パルスシーケンスコマンドは、磁気共鳴フィンガープリンティングプロトコルに従って磁気共鳴データ142を取得するように磁気共鳴イメージングシステム100を制御する。磁気共鳴データ142は、コンピュータストレージ134に記憶されるものとして示される。コンピュータストレージ134は、更に、定常状態磁気共鳴フィンガープリンティングディクショナリ144を含むものとして示される。コンピュータストレージ134はまた、更に、磁気共鳴フィンガープリンティングプロトコルに従って磁気共鳴データ142から再構成された磁気共鳴画像148を含むものとして示される。磁気共鳴画像148は、例えば、所定の物質のセットの様々な物質の濃度の2次元又は3次元マップである。
コンピュータメモリ136は、制御命令150を含むものとして示される。制御命令150は、プロセッサ130が磁気共鳴イメージングシステム100の動作及び機能を制御することを可能にする命令を含む。例えば、制御命令150は、プロセッサ130が、磁気共鳴イメージングシステムを制御するためにパルスシーケンスコマンド140を使用することを可能にするコマンドを含む。コンピュータメモリ136は、更に、磁気共鳴フィンガープリンティングディクショナリ生成命令152を含むものとして示される。磁気共鳴フィンガープリンティングディクショナリ生成命令152はオプションである。磁気共鳴フィンガープリンティングディクショナリ生成命令152は、例えば、定常状態磁気共鳴フィンガープリンティングディクショナリ144を生成するために使用されるモデルを含む。コンピュータメモリ136は、更に、画像再構成命令154を含むものとして示される。画像再構成命令154は、プロセッサ130が磁気共鳴データ142から磁気共鳴画像148を再構成することを可能にする命令を含む。例えば、磁気共鳴画像148は、被検体118内の所定の物質のうちの1つ又は複数の空間分布のレンダリングである。
コンピュータストレージ134及びコンピュータメモリ136の内容は、互いに複製してもよく、又はどちらかの内容が交換されてもよい。図1の例は、磁気共鳴イメージングシステム又は装置100が核磁気共鳴(NMR)分光計と等価となるように変更することができる。装置100は、勾配コイル110及び勾配コイル電源112がない場合、イメージングゾーン108において0次元測定を実行する。
図2は、図1の磁気共鳴イメージングシステム100を動作させる方法の一例を示すフローチャートを示す。最初に、ステップ200において、磁気共鳴データ142の複数のk空間トレースが、パルスシーケンスコマンド140で磁気共鳴イメージングシステム100を制御することによって順次取得される。パルスシーケンスコマンド140は、磁気共鳴フィンガープリンティングプロトコルに従って磁気共鳴データを取得するように磁気共鳴イメージングシステムを制御する。パルスシーケンスコマンド140は、RFパルス列を生成するように磁気共鳴イメージングシステム100を制御するように構成される。パルスシーケンスコマンド140は、磁気共鳴データ142を複数のk空間トレースとして取得するように磁気共鳴イメージングシステムを制御するように構成される。パルスシーケンスコマンド140は、複数のk空間トレースの各々の取得のために繰り返されるべきRFパルス列を制御するように構成される。
この方法は、複数のk空間トレースのうちの初期に取得される所定の数のk空間トレースが取得された後に取得されるk空間トレースに対する所定の物質のセットの各々の存在量を計算するステップ202を更に有する。所定の物質のセットの各々の存在量は、磁気共鳴データを定常状態磁気共鳴フィンガープリンティングディクショナリ144と比較することによって決定される。定常状態磁気共鳴フィンガープリンティングディクショナリ144は、所定の物質のセットに対してRFパルス列に応答して計算された磁気共鳴信号のリストを含む。いくつかの例では、初期に取得された所定の数のk空間トレースは廃棄される。
磁気共鳴フィンガープリンティング(MRF)は、定量的マルチパラメータ測定及び組織特性評価のための有望な技術である。RFパルスの列は、特定の組織タイプのためのフィンガープリントとして働くMR信号を生成する。このフィンガープリントは、後で、既知の信号のディクショナリと比較され得る。
多くの以前に報告された手法は、高速取得のために、長いフィンガープリントシーケンスをスパイラルサンプリングと組み合わせている。しかしながら、MRFシーケンスを短くし、取得を速めるためにデカルトサンプリングをSENSEと組み合わせて使用することも技術的に有利である。この場合、パルス列は、サンプリングされるべきk空間のトレースごとに1回印加されなければならない。高分解能ラジアル又はスパイラルサンプリングについても同様であり、複数のk空間軌道が引き続き測定される。本発明は、k空間を別々の部分でサンプリングするこの場合に焦点を当てる。
インタリーブ式に(デカルトサンプリング技術におけるように)k空間で複数のトレース又はセグメントを引き続きサンプリングするには、通常、フィンガープリント取得の間にスピン系が緩和するのを待つことが必要とされ、その結果、一般的な取得時間は非常に長い。
この欠点は、以下のやり方で克服される可能性がある。スピン系が完全に緩和するのを待つ代わりに、遅延が短くされ、フィンガープリントシーケンスの間のスピン系の時間展開がディクショナリの計算において考慮に入れられる。この目的のためにフィンガープリント信号の定常解を見いだすことができることを示す。
このようにして、デカルト又はさもなければ高密度サンプリングMRFスキャンが、照合精度を損なうことなく著しく速められ得る。
例は、次の特徴のうちの1つ又は複数を含む。
MR信号を取得するMRイメージングシステム。
擬似ランダムMRフィンガープリントシーケンスに従ってデータを取得するプログラム。
以下で説明するように定常状態MRF信号に基づいてMRFディクショナリを計算するソフトウェア。
測定された信号をディクショナリと照合する方法。
ディクショナリエントリを計算する1つのやり方は以下の通りである。
1.列間遅延tをもつ基本MRF列のいくつかの繰り返しからなるシーケンスを構成する。
2.遅延中のスピン系の時間展開を考慮に入れて、完全な合成シーケンスに対する組織パラメータのセットのための予想信号を計算する。計算のタイプは、使用されるMRシーケンスのタイプに依存する。それは、単一スピンのブロッホシミュレーション、スピンの集合のブロッホシミュレーション、拡張位相グラフ計算、又は系を十分に適切に説明する任意の他のスピンモデル計算である。
3.計算された信号から、最後のフィンガープリント繰り返しに対応する部分を切り取り、このフィンガープリント信号をディクショナリエントリとして使用する。
そのとき、MRF測定は、計算で仮定したものと同じ列間遅延tを使用して実行されなければならない。

この例で使用される磁気共鳴フィンガープリントシーケンスは、スポイルド勾配エコーシーケンスに基づく。それは200個のフリップ角からなる列からなり、その前に反転パルスがあり、200ステップに対して3秒の合計時間がある。シーケンスは数回繰り返され、完全なパルス列の予想信号が、拡張位相グラフ形式を使用して計算される。測定は、既知のゲル試料を備えたファントムを使用して実行された。MRFディクショナリは、異なる試料(227ms≦T1≦1646ms;48ms≦T2≦369ms)の18個のエントリ(そのうちの12個はファントムに存在する)と、バックグラウンドのための1つのエントリからなる。
図3は、定常状態に達する磁気共鳴信号の概念を示すために使用される。3つのグラフがある。グラフ300は、磁気共鳴フィンガープリンティングシーケンスで使用される4つの高周波パルスを示す。以下で、グラフ302及びグラフ304には、グラフ300のRFパルスに応答して測定された磁気共鳴信号が示される。グラフ302において、T1緩和時間及びT2緩和時間は、グラフ300のRFパルスの各々の全長と比較して短い。グラフ302では測定された信号の間に小さい差しかないことが分かる。以下で、グラフ304において、T1時間は、グラフ300のRFパルスと比較して長い。第1のMR信号306は、後続の磁気共鳴信号と根本的に異なることが分かる。第1のMR信号306が下げられるか又は異なる磁気共鳴フィンガープリンティングディクショナリに適合する場合、測定される後続の磁気共鳴信号は定常状態磁気共鳴フィンガープリンティングディクショナリに適合することが分かる。
図4は、複数の既知の試料によるファントムの3つの磁気共鳴フィンガープリンティング画像400、402、404を示す。これらの3つの画像では、既知のゲル試料を備えたファントムが、18個の異なるエントリをもつ磁気共鳴フィンガープリンティングディクショナリと比較されている。この例では、ファントムは、18個の異なる試料のうちの12個を含む。画像400の磁気共鳴データを取得するとき、RFパルス列の高周波パルスの各々の間に5秒の遅延があった。画像400では、12個の試料の各々が正しく識別された。画像402では、RFパルスの各々の間の遅延はわずか0.5秒であった。402から404において識別された試料を比較すると、画像402の左上の1、2、及び3とラベル付けされた試料の3つのみが、正しく識別されていることが分かる。他の試料は、間違って識別され、それが間違っていることを示すためにそれを通る対角線で印をつけられている。
画像404では、定常状態磁気共鳴フィンガープリンティングディクショナリが使用された。試料の全てが、画像400でどのように識別されているかと一致して識別されていることが分かる。画像404において、遅延は、同様に、わずか0.5秒であった。画像400の取得は8分33秒かかった。画像404の取得は3分45秒しかかからなかった。これは、定常状態磁気共鳴フィンガープリンティングディクショナリの使用が、磁気共鳴フィンガープリンティングプロトコルのための磁気共鳴データの取得を根本的に速めるためにどのように利用できるかを示している。
繰り返し間の遅延が短い場合、スピン系は次の反転パルスの前に完全には緩和していない。ショートT1系(図3、画像302)の信号の応答は、繰り返しごとにほとんど同一であるが、ロングT1系(図3、画像304)は異なる信号応答を示す。しかしながら、3つの繰り返しの後、この系でさえ定常フィンガープリントに達する。第1の測定では、MRFシーケンス繰り返し間の遅延は大きく(5s)、その結果、完全に緩和されたスピン系を毎回仮定することができる。この仮定で計算されたディクショナリを使用して、全ての12個のサンプル番号が正しく照合されている(図4、画像400)。遅延が0.5sに減少された場合、緩和されたスピン系の仮定がもはや有効でないので、同じディクショナリは間違った照合をもたらす(図4、画像402)。3つのシーケンス繰り返しからなる列を計算し、ディクショナリに対して定常状態フィンガープリント信号を使用する場合、全ての試料が再び正しく照合される(図4、404)。第1の測定は8分33秒かかり、第2の測定は3分45秒しかかからない。
別の例では、最初の取得がk空間の最下位部分で実行され、その結果、初期(過渡)MRF信号は、結果として生じる画像をあまり歪ませない。
別の例では、データの取得が始まる前に、MRFシーケンスの少数の「ダミーの」繰り返しが系に適用される。このようにして、系は、最初の取得時には既に定常状態にある。
別の例では、計算で使用する繰り返しの数は、決定されるべき個々のディクショナリエントリのT1値及びT2値に動的に適合される。繰り返しの数は、最後の繰り返しで系が定常状態に達しているように選ばれる。緩和定数がより小さい系は定常状態に早く近づくので、計算時間をここで節約することができる。
更なる例では、2~3のスポークによるラジアルサンプリング又はいくつかのインタリーブによるスパイラルサンプリングが、列間遅延t=0を用いて受容的に実行される。
本発明は、図面及び前述の記載において詳細に図示及び説明されたが、このような図示及び記載は、説明的又は例示的であって限定するものではないと見なされるべきである。すなわち本発明は、開示された実施形態に限定されるものではない。
開示された実施形態のその他の変形が、図面、本開示及び添付の請求項の検討から、請求項に係る発明を実施する当業者によって理解されて実現され得る。請求項において、「comprising(含む、備える)」という単語は、他の要素又はステップを除外するものではなく、不定冠詞「a」又は「an」は、複数を除外するものではない。単一のプロセッサ又は他のユニットが請求項に記載されたいくつかのアイテムの機能を果たす。特定の手段が相互に異なる従属請求項に列挙されているという単なる事実は、これらの手段の組み合わせが有利に用いられないことを示すものではない。コンピュータプログラムは、他のハードウェアと共に若しくは他のハードウェアの一部として供給される光記憶媒体又はソリッドステート媒体等の適当な媒体に保存/分配されてもよいが、インターネット又は他の有線若しくは無線の電気通信システムを介して等の他の形式で分配されてもよい。請求項における任意の参照符号は、本発明の範囲を限定するものと解釈されるべきではない。
1 ゲル試料1
2 ゲル試料2
3 ゲル試料3
4 ゲル試料4
6 ゲル試料6
7 ゲル試料7
8 ゲル試料8
9 ゲル試料9
10 ゲル試料10
11 ゲル試料11
12 ゲル試料12
13 ゲル試料13
14 ゲル試料14
100 磁気共鳴イメージングシステム
104 磁石
106 磁石のボア
108 測定ゾーン又はイメージングゾーン
110 磁場勾配コイル
112 磁場勾配コイル電源
114 高周波コイル
116 送受信器
118 被検体
120 被検体支持体
122 アクチュエータ
124 所定の方向
125 スライス
126 コンピュータシステム
128 ハードウェアインタフェース
130 プロセッサ
132 ユーザインタフェース
134 コンピュータストレージ
136 コンピュータメモリ
140 パルスシーケンスコマンド
142 磁気共鳴データ
144 定常状態磁気共鳴フィンガープリンティングディクショナリ
148 磁気共鳴画像
150 制御命令
152 磁気共鳴フィンガープリンティングディクショナリ生成命令
154 画像再構成命令
200 パルスシーケンスコマンドで磁気共鳴イメージングシステムを制御することによって、磁気共鳴データの複数のk空間トレースを順次取得するステップ
202 複数のk空間トレースのうちの所定の数のk空間トレースが取得された後に取得されるk空間トレースに対する所定の物質のセットの各々の存在量を計算するステップ
300 RFパルス
302 MR信号
304 MR信号
306 第1のMR信号

Claims (13)

  1. 測定ゾーン内の被検体から磁気共鳴データを取得するための磁気共鳴イメージングシステムであって、前記磁気共鳴イメージングシステムは、
    前記磁気共鳴イメージングシステムを制御するためのプロセッサと、
    マシン実行可能命令及びパルスシーケンスコマンドを記憶するためのメモリとを含み、
    前記パルスシーケンスコマンドが、磁気共鳴フィンガープリンティングプロトコルに従って前記磁気共鳴データを取得するように前記磁気共鳴イメージングシステムを制御し、前記パルスシーケンスコマンドが、RFパルス列を生成するように前記磁気共鳴イメージングシステムを制御し、前記パルスシーケンスコマンドが、前記磁気共鳴データを複数のk空間トレースとして取得するように前記磁気共鳴イメージングシステムを制御し、前記パルスシーケンスコマンドが、前記複数のk空間トレースの各々の前記取得のために繰り返されるべき前記RFパルス列を制御し、前記マシン実行可能命令によって、前記プロセッサが、
    パルスシーケンスコマンドで前記磁気共鳴イメージングシステムを制御することによって、磁気共鳴データの前記複数のk空間トレースを順次取得することと、
    初期に取得される所定の数のk空間トレースが取得された後でのみ、前記k空間トレースを定常状態磁気共鳴フィンガープリンティングディクショナリと比較することであって、前記所定の数のk空間トレースが取得された後、磁化は平衡状態に達する、当該比較することと、
    前記複数のk空間トレースのうちの前記所定の数のk空間トレースが取得された後に取得されるk空間トレースに対する所定の物質のセットの各々の存在量を計算することとを実行させ、所定の物質のセットの各々の前記存在量が、前記磁気共鳴データを前記定常状態磁気共鳴フィンガープリンティングディクショナリと比較することによって決定され、前記定常状態磁気共鳴フィンガープリンティングディクショナリが、所定の物質のセットに対して前記RFパルス列に応答して計算された磁気共鳴信号のリストを含み、
    前記所定の数のk空間トレースが、決定されるべき個々のディクショナリエントリのT1値及びT2値に動的に適合される、磁気共鳴イメージングシステムにおいて、
    前記マシン実行可能命令の実行によって、前記プロセッサは、前記定常状態磁気共鳴フィンガープリンティングディクショナリの計算の間に前記所定の数を決定し、前記所定の数が、前記定常状態磁気共鳴フィンガープリンティングディクショナリの前記計算の間に生成されたシミュレート磁気共鳴データの収束を決定する収束基準を使用して決定される、磁気共鳴イメージングシステム
  2. 前記マシン実行可能命令の実行によって、前記プロセッサが、前記所定の数のk空間トレースを廃棄する、請求項1に記載の磁気共鳴イメージングシステム。
  3. 前記パルスシーケンスコマンドが、各RFパルス列の間の最大遅延で前記RFパルス列を繰り返すように前記磁気共鳴イメージングシステムを制御し、前記最大遅延が、5秒未満、1秒未満、0.5秒未満、0.1秒未満、0.05秒未満、及び0.00秒未満のうちのいずれかである、請求項1又は2に記載の磁気共鳴イメージングシステム。
  4. 前記磁気共鳴イメージングシステムがメイン磁場をもつ磁石を含み、前記メイン磁場が、前記測定ゾーン内に平均の磁場の大きさを有し、前記マシン実行可能命令の実行によって、前記プロセッサが、各RFパルス列を、前記平均の磁場の大きさにおける前記所定の物質のセットの最大T1時間の5倍、前記平均の磁場の大きさにおける前記所定の物質のセットの最大T1時間の1倍、前記平均の磁場の大きさにおける前記所定の物質のセットの最大T1時間の0.5倍、及び前記平均の磁場の大きさにおける前記所定の物質のセットの最大T1時間の0.1倍のうちのいずれかよりも速く繰り返す、請求項1又は2に記載の磁気共鳴イメージングシステム。
  5. 前記マシン実行可能命令の実行によって、更に、前記プロセッサが、前記定常状態磁気共鳴フィンガープリンティングディクショナリを計算する、請求項1から4のいずれか一項に記載の磁気共鳴イメージングシステム。
  6. 前記マシン実行可能命令の実行によって、更に、前記プロセッサが、前記所定の物質のセットの各々への前記RFパルス列の繰り返し印加をモデリングすることによって前記定常状態磁気共鳴フィンガープリンティングディクショナリを計算する、請求項5に記載の磁気共鳴イメージングシステム。
  7. 前記パルスシーケンスコマンドが、k空間において前記複数のk空間トレースを回転させること、ラジアルサンプリングを使用して前記複数のk空間トレースを取得すること、デカルトサンプリングを使用して前記複数のk空間トレースを取得すること、及び非デカルトサンプリングを使用して前記複数のk空間トレースを取得することのうちのいずれかを実行するように、前記磁気共鳴イメージングシステムを制御する、請求項1からのいずれか一項に記載の磁気共鳴イメージングシステム。
  8. 前記パルスシーケンスコマンドが、PROPELLER磁気共鳴イメージングプロトコルに従って前記磁気共鳴データを取得するように前記磁気共鳴イメージングシステムを制御する、請求項1からのいずれか一項に記載の磁気共鳴イメージングシステム。
  9. 前記パルスシーケンスコマンドが、前記複数のk空間トレースをk空間におけるスパイラルとして取得するように前記磁気共鳴イメージングシステムを制御する、請求項1からのいずれか一項に記載の磁気共鳴イメージングシステム。
  10. 前記パルスシーケンスコマンドが、前記複数のk空間トレースをデカルトk空間におけるラインとして取得するように前記磁気共鳴イメージングシステムを制御する、請求項1からのいずれか一項に記載の磁気共鳴イメージングシステム。
  11. 前記パルスシーケンスコマンドが、前記磁気共鳴データを、並列イメージング磁気共鳴イメージングプロトコル、SENSE磁気共鳴イメージングプロトコル、及びGRAPPA磁気共鳴イメージングプロトコルのうちのいずれかに従って取得するように前記磁気共鳴イメージングシステムを制御する、請求項10に記載の磁気共鳴イメージングシステム。
  12. 測定ゾーン内の被検体から磁気共鳴データを取得するように磁気共鳴イメージングシステムを制御するプロセッサによる実行のためのマシン実行可能命令を含むコンピュータプログラムであって、前記マシン実行可能命令によって、前記プロセッサが、
    磁気共鳴フィンガープリンティングプロトコルに従って前記磁気共鳴データを取得するように前記磁気共鳴イメージングシステムを制御するためのパルスシーケンスコマンドで前記磁気共鳴イメージングシステムを制御することによって、磁気共鳴データの複数のk空間トレースを順次取得することであって、前記パルスシーケンスコマンドが、RFパルス列を生成するように前記磁気共鳴イメージングシステムを制御し、前記パルスシーケンスコマンドが、前記磁気共鳴データを複数のk空間トレースとして取得するように前記磁気共鳴イメージングシステムを制御し、前記パルスシーケンスコマンドが、前記複数のk空間トレースの各々の前記取得のために繰り返されるべき前記RFパルス列を制御する、当該取得することと、
    初期に取得される所定の数のk空間トレースが実行された後でのみ、前記k空間トレースを定常状態磁気共鳴フィンガープリンティングディクショナリと比較することであって、前記所定の数のk空間トレースが取得された後、磁化は平衡状態に達する、当該比較することと、
    前記複数のk空間トレースのうちの前記所定の数のk空間トレースが取得された後に取得されるk空間トレースに対する所定の物質のセットの各々の存在量を計算することであって、所定の物質のセットの各々の前記存在量が、前記磁気共鳴データを前記定常状態磁気共鳴フィンガープリンティングディクショナリと比較することによって決定され、前記定常状態磁気共鳴フィンガープリンティングディクショナリが、所定の物質のセットに対して前記RFパルス列に応答して計算された磁気共鳴信号のリストを含み、前記所定の数のk空間トレースが、決定されるべき個々のディクショナリエントリのT1値及びT2値に動的に適合される、当該計算することとを行う、コンピュータプログラムにおいて、
    前記マシン実行可能命令の実行によって、前記プロセッサは、前記定常状態磁気共鳴フィンガープリンティングディクショナリの計算の間に前記所定の数を決定し、前記所定の数が、前記定常状態磁気共鳴フィンガープリンティングディクショナリの前記計算の間に生成されたシミュレート磁気共鳴データの収束を決定する収束基準を使用して決定される、コンピュータプログラム
  13. 測定ゾーン内の被検体から磁気共鳴データを取得するための磁気共鳴イメージングシステムを動作させる方法であって、前記方法は、
    磁気共鳴フィンガープリンティングプロトコルに従って前記磁気共鳴データを取得するように前記磁気共鳴イメージングシステムを制御するためのパルスシーケンスコマンドで前記磁気共鳴イメージングシステムを制御することによって、磁気共鳴データの複数のk空間トレースを順次取得するステップであって、前記パルスシーケンスコマンドが、RFパルス列を生成するように前記磁気共鳴イメージングシステムを制御し、前記パルスシーケンスコマンドが、前記磁気共鳴データを複数のk空間トレースとして取得するように前記磁気共鳴イメージングシステムを制御し、前記パルスシーケンスコマンドが、前記複数のk空間トレースの各々の前記取得のために繰り返されるべき前記RFパルス列を制御する、当該取得するステップと、
    初期に取得される所定の数のk空間トレースが実行された後でのみ、前記k空間トレースを定常状態磁気共鳴フィンガープリンティングディクショナリと比較するステップであって、前記所定の数のk空間トレースが取得された後、磁化は平衡状態に達する、当該比較するステップと、
    前記複数のk空間トレースのうちの前記所定の数のk空間トレースが取得された後に取得されるk空間トレースに対する所定の物質のセットの各々の存在量を計算するステップであって、所定の物質のセットの各々の前記存在量が、前記磁気共鳴データを前記定常状態磁気共鳴フィンガープリンティングディクショナリと比較することによって決定され、前記定常状態磁気共鳴フィンガープリンティングディクショナリが、所定の物質のセットに対して前記RFパルス列に応答して計算された磁気共鳴信号のリストを含み、前記所定の数のk空間トレースが、決定されるべき個々のディクショナリエントリのT1値及びT2値に動的に適合される、当該計算するステップと
    を有する、方法において、
    前記定常状態磁気共鳴フィンガープリンティングディクショナリの計算の間に所定の数を決定するステップを更に有し、前記所定の数が、前記定常状態磁気共鳴フィンガープリンティングディクショナリの前記計算の間に生成されたシミュレート磁気共鳴データの収束を決定する収束基準を使用して決定される、方法
JP2018541247A 2016-02-09 2017-02-06 定常状態磁気共鳴フィンガープリンティング Active JP7059189B6 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP16154766.6 2016-02-09
EP16154766 2016-02-09
PCT/EP2017/052481 WO2017137333A1 (en) 2016-02-09 2017-02-06 Steady state magnetic resonance fingerprinting

Publications (4)

Publication Number Publication Date
JP2019512282A JP2019512282A (ja) 2019-05-16
JP2019512282A5 JP2019512282A5 (ja) 2020-03-12
JP7059189B2 true JP7059189B2 (ja) 2022-04-25
JP7059189B6 JP7059189B6 (ja) 2022-06-03

Family

ID=55349681

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018541247A Active JP7059189B6 (ja) 2016-02-09 2017-02-06 定常状態磁気共鳴フィンガープリンティング

Country Status (5)

Country Link
US (1) US10788556B2 (ja)
EP (1) EP3414584B1 (ja)
JP (1) JP7059189B6 (ja)
CN (1) CN108603921B (ja)
WO (1) WO2017137333A1 (ja)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10705168B2 (en) * 2017-01-17 2020-07-07 Case Western Reserve University System and method for low rank approximation of high resolution MRF through dictionary fitting
CN111492256B (zh) * 2017-12-18 2023-07-14 皇家飞利浦有限公司 运动补偿的磁共振成像
DE102019206827A1 (de) 2019-05-10 2020-11-12 Siemens Healthcare Gmbh Magnetresonanz-Fingerprinting-Verfahren

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015119866A (ja) 2013-12-24 2015-07-02 株式会社東芝 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
JP2015528384A (ja) 2012-09-19 2015-09-28 ケース ウェスタン リザーブ ユニバーシティCase Westernreserve University 核磁気共鳴(nmr)フィンガープリンティング

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7348776B1 (en) * 2006-09-01 2008-03-25 The Board Of Trustees Of The Leland Stanford Junior University Motion corrected magnetic resonance imaging
JP5132681B2 (ja) * 2007-07-02 2013-01-30 株式会社日立メディコ 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
US9645211B2 (en) * 2010-04-29 2017-05-09 The Board Of Trustees Of The University Of Illinois Magnetic resonance imaging using steering-propeller
US10627468B2 (en) 2011-03-18 2020-04-21 Case Western Reserve University Nuclear magnetic resonance (NMR) fingerprinting
US9097781B2 (en) 2012-04-12 2015-08-04 Mark Griswold Nuclear magnetic resonance (NMR) fingerprinting with parallel transmission
US9568579B2 (en) * 2012-12-19 2017-02-14 Case Western Reserve University Magnetic resonance fingerprinting (MRF) with echo splitting
DE102013201814B4 (de) 2013-02-05 2018-11-08 Siemens Healthcare Gmbh Verfahren zur Magnetresonanz-Bildgebung mit multidimensional ortsselektiven HF-Pulsen in einem Außenbereich
US9625547B2 (en) 2013-03-28 2017-04-18 Universitatsspital Basel Magnetic resonance imaging method for the quantification of the T1 and/or T2 relaxation times in a sample
WO2015175028A1 (en) 2014-02-11 2015-11-19 The General Hospital Corporation Systems and methods for acceleration magnetic resonance fingerprinting
WO2015138941A1 (en) 2014-03-14 2015-09-17 The General Hospital Corporation System and method for magnetic resonance fingerprinting
US10335049B2 (en) 2014-04-22 2019-07-02 Case Western Reserve University Magnetic resonance fingerprinting with steady state precession (MRF-FISP)
WO2015181091A2 (en) * 2014-05-27 2015-12-03 Koninklijke Philips N.V. Mr imaging using multi-echo k-space acquisition
US9897675B2 (en) * 2014-05-28 2018-02-20 Case Western Reserve University Magnetic resonance fingerprinting (MRF) with simultaneous multivolume acquisition
US10422845B2 (en) 2014-10-24 2019-09-24 The General Hospital Corporation Systems and methods for steady-state magnetic resonance fingerprinting
US10241178B2 (en) 2015-05-29 2019-03-26 Case Western Reserve University System and method for magnetic resonance fingerprinting at high field strengths

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015528384A (ja) 2012-09-19 2015-09-28 ケース ウェスタン リザーブ ユニバーシティCase Westernreserve University 核磁気共鳴(nmr)フィンガープリンティング
JP2015119866A (ja) 2013-12-24 2015-07-02 株式会社東芝 磁気共鳴イメージング装置及び磁気共鳴イメージング方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Guido Buonincontri et al.,MR Fingerprinting with Simultaneous B1 Estimation,Magn Reson Med,2015年10月28日,vol.76,pp.1127-1135
Yun Jiang et al.,MR Fingerprinting Using Fast Imaging with Steady State Precession (FISP) with Spiral Readout,Magn Reson Med,2014年12月09日,vol.74,pp.1621-1631

Also Published As

Publication number Publication date
EP3414584B1 (en) 2023-08-23
JP7059189B6 (ja) 2022-06-03
EP3414584A1 (en) 2018-12-19
US10788556B2 (en) 2020-09-29
WO2017137333A1 (en) 2017-08-17
CN108603921B (zh) 2021-08-10
CN108603921A (zh) 2018-09-28
US20190033413A1 (en) 2019-01-31
JP2019512282A (ja) 2019-05-16

Similar Documents

Publication Publication Date Title
US10321845B2 (en) Magnetic resonance fingerprinting in slices along a one-dimensional extension
RU2693837C2 (ru) Магнитно-резонансный метод пальцевых отпечатков
JP6339299B1 (ja) 主磁場の不均一性に対する感度が減少された磁気共鳴フィンガープリンティング
CN110073232B (zh) 用于多状态磁共振指纹识别的磁共振成像系统、方法和计算机可读介质
JP7041072B2 (ja) サイレント3d磁気共鳴フィンガープリンティング
CN109791186B (zh) 对磁共振指纹期间的b0偏共振场的直接测量
JP7090636B2 (ja) サブボクセル解像度磁気共鳴フィンガープリンティングイメージング
JP7080316B2 (ja) 磁気共鳴イメージングにおける緩和時間の定量的測定
JP2018522628A (ja) Mrf辞書生成
JP7059189B2 (ja) 定常状態磁気共鳴フィンガープリンティング
CN110114685B (zh) 使用最大似然重建的对心脏的t1映射
US11953571B2 (en) Quantatative magnetic resonance imaging in the presence of fat
JP2020525157A (ja) フーリエ空間における磁気共鳴フィンガプリンティング

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200203

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200203

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20210126

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210222

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210806

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20211102

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220131

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20220315

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220413

R150 Certificate of patent or registration of utility model

Ref document number: 7059189

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150