JP6339299B1 - 主磁場の不均一性に対する感度が減少された磁気共鳴フィンガープリンティング - Google Patents

主磁場の不均一性に対する感度が減少された磁気共鳴フィンガープリンティング Download PDF

Info

Publication number
JP6339299B1
JP6339299B1 JP2017553237A JP2017553237A JP6339299B1 JP 6339299 B1 JP6339299 B1 JP 6339299B1 JP 2017553237 A JP2017553237 A JP 2017553237A JP 2017553237 A JP2017553237 A JP 2017553237A JP 6339299 B1 JP6339299 B1 JP 6339299B1
Authority
JP
Japan
Prior art keywords
magnetic resonance
gradient
pulse sequence
magnetic
pulse
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.)
Expired - Fee Related
Application number
JP2017553237A
Other languages
English (en)
Other versions
JP2018516622A (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
Application granted granted Critical
Publication of JP6339299B1 publication Critical patent/JP6339299B1/ja
Publication of JP2018516622A publication Critical patent/JP2018516622A/ja
Expired - Fee Related 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/445MR involving a non-standard magnetic field B0, e.g. of low magnitude as in the earth's magnetic field or in nanoTesla spectroscopy, comprising a polarizing magnetic field for pre-polarisation, B0 with a temporal variation of its magnitude or direction such as field cycling of B0 or rotation of the direction of B0, or spatially inhomogeneous B0 like in fringe-field MR or in stray-field 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/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/50NMR imaging systems based on the determination of relaxation times, e.g. T1 measurement by IR sequences; T2 measurement by multiple-echo sequences
    • 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/543Control 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
    • 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/546Interface between the MR system and the user, e.g. for controlling the operation of the MR system or for the design of pulse sequences
    • 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
    • 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
    • G01R33/5612Parallel RF transmission, i.e. RF pulse transmission using a plurality of independent transmission channels
    • 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/5615Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5617Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using RF refocusing, e.g. RARE
    • 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/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56563Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0

Landscapes

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

Abstract

本発明は、測定領域内に主磁場を生成する磁石104と、少なくとも1つの方向それぞれの傾斜磁場コイル112のセットに電流を供給することによって、少なくとも1つの方向において、測定領域内に傾斜磁場を生成する傾斜磁場システム110、112とを含む磁気共鳴システム100を提供する。命令は、プロセッサ130に磁気共鳴システムを制御させる。マシン実行可能命令の実行は、プロセッサに、パルスシーケンスコマンドを用いて磁気共鳴システムを制御させることによって、磁気共鳴データを取得させる(200)。パルスシーケンスコマンド140は、磁気共鳴システムに、磁気共鳴フィンガープリンティング技術に従って磁気共鳴データを取得させる。パルスシーケンスコマンドは、一連500ノパルスシーケンス反復502、504を特定し、それぞれ、固定反復時間302を有する。各反復は、パルスシーケンス反復の開始から固定遅延316において生じる、無線周波数パルスの分布から選択される無線周波数パルス310か又はサンプリングイベント404を含む。パルスシーケンスコマンドは、傾斜磁場コイルのセットへの供給電流を制御することによる少なくとも1つの方向における傾斜磁場308の印加を特定する。傾斜磁場コイルそれぞれに対して、供給電流の積分は、各固定反復時間に対して一定である。命令は更に、プロセッサに、磁気共鳴データを磁気共鳴フィンガープリンティングディクショナリ144と比較させることによって、所定物質のセットの各物質の存在量を計算させる(202)。

Description

本発明は、磁気共鳴撮像に関し、具体的には、磁気共鳴フィンガープリンティングを行う技術に関する。
磁気共鳴(MR)フィンガープリンティングは、様々な材料及び組織からの信号が測定MR信号に一意の寄与を有するように、時間的に分布した幾つかのRFパルスが印加される新しい技術である。一組の又は固定数の物質からの事前に計算した信号寄与の限定ディクショナリが測定MR信号と比較され、単一のボクセルにおける組成を決定することができる。例えばボクセルが、水分、脂肪及び筋肉組織しか含まないことが分かっていると、これらの3つの材料からの寄与のみを考慮すればよく、当該ボクセルの組成を正確に決定するためには、幾つかのRFパルスしか必要でない。より高い分解能を有するより大きいディクショナリが使用される場合、MRフィンガープリンティングを使用して、ボクセルの様々な組織パラメータ(例えばT1、T2、…)を同時にかつ定量的に決定することができる。
磁気共鳴フィンガープリンティング技術は、Ma他によるジャーナル論文「Magnetic Resonance Fingerprinting」(Nature、第495巻、187〜193頁、doi:10.1038/nature11971)において紹介された。磁気共鳴フィンガープリンティング技術は、米国特許出願公開第2013/0271132A1号及び第2013/0265047A1号においても説明されている。
Jiang他による会議議事録「MR Fingerprinting Using Spiral QUEST」(Proc、Intl.Soc.Mag.Reson.Med.21(2013)、0019頁)では、クイック・エコー・スプリット撮像技術(QUEST)を、磁気共鳴フィンガープリンティング(MRF)シーケンスの構成単位として使用するMRFが開示されている。
本発明は、独立請求項における磁気共鳴撮像システム、コンピュータプログラムプロダクト及び方法を提供する。実施形態は、従属請求項に提供される。
Ma他によるNatureの論文は、磁気共鳴フィンガープリンティングの基本概念と、本明細書では「事前計算した磁気共鳴フィンガープリンティングディクショナリ」、「磁気共鳴フィンガープリンティングディクショナリ」及び「ディクショナリ」と呼ぶディクショナリといったこの技術を説明するために使用される専門用語とを紹介している。
当業者には理解されるように、本発明の態様は、装置、方法又はコンピュータプログラムプロダクトとして具体化されてよい。したがって、本発明の態様は、完全にハードウェアの実施形態、完全にソフトウェアの実施形態(ファームウェア、常駐ソフトウェア、マイクロコード等を含む)又はソフトウェア態様とハードウェア態様とを組み合わせた実施形態の形を取ってよく、これらはすべて、本明細書では、一般的に、「回路」、「モジュール」又は「システム」と呼ぶ。更に、本発明の態様は、コンピュータ実行可能なコードが具体化されている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ポート、IEEE−488ポート、ブルートゥース(登録商標)接続、ワイヤレスローカルエリアネットワーク接続、TCP/IP接続、イーサネット(登録商標)接続、制御電圧インターフェース、MIDIインターフェース、アナログ入力インターフェース及びデジタル入力インターフェースが挙げられるが、これらに限定されない。
「ディスプレイ」又は「ディスプレイデバイス」は、本明細書において使用される場合、画像又はデータを表示する出力デバイス又はユーザインターフェースを包含する。ディスプレイは、視覚、音声及び/又は触知データを出力する。ディスプレイの例としては、コンピュータモニタ、テレビジョンスクリーン、タッチスクリーン、触知電子ディスプレイ、点字(Braille)スクリーン、陰極線管(CRT)、蓄積管、双安定ディスプレイ、電子ペーパー、ベクトルディスプレイ、フラットパネルディスプレイ、真空蛍光ディスプレイ(VF)、発光ダイオード(LED)ディスプレイ、エレクトロルミネッセンスディスプレイ(ELD)、プラズマディスプレイパネル(PDP)、液晶ディスプレイ(LCD)、有機発光ダイオードディスプレイ(OLED)、プロジェクタ及び頭部取付け式ディスプレイが挙げられるが、これらに限定されない。
磁気共鳴(MR)データは、本明細書において、磁気共鳴撮像スキャン中に、磁気共鳴装置のアンテナを使用して、原子スピンによって放射された無線周波数信号の記録された測定結果と定義される。磁気共鳴データは、医用画像データの一例である。磁気共鳴撮像(MRI)画像は、本明細書において、磁気共鳴撮像データ内に含まれる解剖学的データの再構成された2又は3次元視覚化と定義される。この視覚化は、コンピュータを使用して行われてよい。
1つの態様では、測定領域内の被験者からの磁気共鳴データを取得する磁気共鳴システムを提供する。様々な例において、磁気共鳴システムは、様々な形を取る。例えば一例では、磁気共鳴システムはNMRスペクトロメータである。別の例では、磁気共鳴システムは磁気共鳴撮像システムである。磁気共鳴システムがNMRスペクトロメータである場合、被験対象は、測定領域内に置かれた何らかの種類の容器内の化学的被験対象である。磁気共鳴システムは、測定領域内に主磁場を生成する磁石を含む。
磁気共鳴システムは更に、少なくとも1つの方向それぞれの傾斜磁場コイルのセットに電流を供給することによって、少なくとも1つの方向において、測定領域内に傾斜磁場を生成する傾斜磁場システムを含む。つまり、傾斜磁場が生成される各方向に対する傾斜磁場コイルのセットがある。磁気共鳴撮像システムでは、典型的に、傾斜磁場コイルの3つの直交セットがある。NMRスペクトロメータの場合、傾斜磁場コイルの1つ、2つ又は3つのセットがある。
磁気共鳴システムは更に、マシン実行可能命令を格納するメモリを含む。メモリは更に、パルスシーケンスコマンドを格納する。パルスシーケンスコマンドは、磁気共鳴システムに、磁気共鳴フィンガープリンティング技術に従って磁気共鳴データを取得させる。パルスシーケンスコマンドは、一連のパルスシーケンス反復を特定する。これは更に、パルスシーケンスコマンドは、磁気共鳴撮像システムによって行われるパルスシーケンス反復のシーケンス又は一連のパルスシーケンス反復を特定するとも解釈される。各パルスシーケンス反復は、固定反復時間を有する。つまり、各パルス反復は、同じ持続時間を有する。各パルスシーケンス反復は、パルスシーケンス反復の始まりから固定遅延において生じる無線周波数パルスか又はサンプリングイベントを含む。各パルスシーケンス反復内で、無線周波数パルス又はサンプリングイベントのどちらかが生じるが、両方は生じない。
ほとんどの場合、固定遅延は、パルスシーケンス反復の持続時間と等しく、RFパルス又はサンプリングイベントは、パルスシーケンス反復の終わりにおいて生じる。しかし、特定のパルスシーケンス反復がいつ開始又は終了するかを測定できる様々な方法がある。したがって、固定遅延は、パルスシーケンス反復がいつ始まるか又は開始するかを解釈するこれらの様々な方法を表現する手段である。無線周波数パルス又はサンプリングイベントは、毎回、パルスシーケンス反復内に同じ時間にあるように発生する。無線周波数パルスは、無線周波数パルスの分布から選択される。無線周波数パルスの分布は、フリップ角の分布に向けて磁気スピンを回転させる。パルスシーケンスコマンドは、傾斜磁場コイルのセットに供給される電流を制御することによる少なくとも1つの方向における傾斜磁場の印加を特定する。パルスシーケンス反復は更に、所定の(選択された)基本時間単位を有してよく、これは、RF励起パルス又は信号読出しを任意に含んでよい。基本時間単位は、RF励起及び信号読出し(取得)といったイベントが行われる最小時間単位である。連続する基本時間単位は、同じ又は異なるイベントを含んでよく、また、基本時間単位の連続は、(シーケンスの反復時間を規定する)自己繰り返し単位を含んでも含まなくてもよい。しかし、様々なイベントを有する基本時間単位の(疑似ランダム)非反復連続があってもよい。
傾斜磁場コイルのセットの各傾斜磁場コイルに対して、供給電流の積分は、各固定反復時間に対して一定である。無線周波数パルスは、無線周波数パルスの分布から選択され、したがって、被験者の様々な特性を、磁気共鳴フィンガープリンティング技術を使用して試験することができる。磁気共鳴システムは、当該磁気共鳴システムを制御するプロセッサを含む。本明細書では、プロセッサは、1つ以上のプロセッサと、更にはコントローラとを包含するように使用される。マシン実行可能命令の実行は、プロセッサに、パルスシーケンスコマンドを用いて磁気共鳴システムを制御することによって磁気共鳴データを取得させる。
マシン実行可能命令の実行は更に、プロセッサに、磁気共鳴データを磁気共鳴フィンガープリンティングディクショナリと比較することによって、所定物質のセットの各物質の存在量を計算させる。磁気共鳴フィンガープリンティングディクショナリは、所定物質のセットに対するパルスシーケンスコマンドの実行に反応して計算された磁気共鳴信号のリストを含む。典型的に、磁気共鳴フィンガープリンティング技術は、いつすべてのサンプリングイベントが生じ、また、いつすべての無線周波数パルスが生じるかを規定する特定のパルスシーケンスに対して行われる。無線周波数パルスの分布から選択される無線周波数パルスのシーケンスも規定される。非常に典型的に、磁気共鳴フィンガープリンティングディクショナリは、使用される特定のパルスシーケンスコマンドに合わせられている。磁気共鳴データを磁気共鳴フィンガープリンティングディクショナリと比較することによって、例えば所定物質のどの物質が又はどれくらいの量で、磁気共鳴システムによって測定される様々なボクセル又はボリューム内にあるかが推測される。
上記方法は、磁気共鳴フィンガープリンティング技術は、測定領域内のB0、即ち、主磁場の不均一性によってあまり影響を受けないという利点を有する。使用された特定のパルスシーケンスは、反復時間と傾斜磁場コイルの供給電流の積分の詳述との組み合わせを使用し、B0不均一性への測定結果の依存を減少させる方法がもたらされる。これは、特に磁気共鳴撮像にメリットがある。これは、1つ以上の傾斜磁場コイルを追加することによって、核磁気共鳴器具、即ち、NMRスペクトロメータにも有益である。NMRスペクトロメータの傾斜磁場コイルは、空間符号化には使用されないが、主磁場のB0不均一性の作用を減少させるために使用される。
磁気共鳴フィンガープリンティングディクショナリは、所定物質のセットに対するパルスシーケンスコマンドの励起に反応して計算した磁気共鳴信号のリストを含む。
パルスシーケンスコマンドが実行されると、パルスシーケンス反復は、1つずつ実行される。これは、サンプリング時間中に、各パルスシーケンス反復に対してデータ取得されることにつながる。磁気共鳴フィンガープリンティングディクショナリは、特定の物質に対する期待磁気共鳴信号を含む。測定信号は、1ボクセルに含まれる様々な物質からの磁気共鳴信号の一次結合である。ボクセルサイズに応じて、1ボクセル内に1つ以上の物質があってよい。磁気共鳴フィンガープリンティング技術では、様々な物質の可能な組成が検討される。各物質の可能なフィンガープリントが、実測物質と比較され、物質の組成が、磁気共鳴フィンガープリンティングディクショナリを使用して分解される。
通常は、磁気共鳴フィンガープリンティングディクショナリが計算される場合、磁場の不均一性を考慮する必要がある。ボクセルサイズが、空間場変化と比較して小さい場合、多数の様々な磁場に対する計算信号応答を含むディクショナリが、十分に優れたマッチを提供することができる。より大きいボクセルサイズは、所定物質のセットの各物質に対して基本的に不鮮明なフィンガープリントをもたらす場合がある。上記パルスシーケンスコマンドの使用は、磁気共鳴フィンガープリンティング技術へのB0不均一性の作用を減少することができる。
別の実施形態では、無線周波数パルス又はサンプリングイベントは、固定遅延の時間に中心がある。つまり、無線周波数パルスの中心又はサンプリングイベントの取得時間の中間点は、固定遅延の時間にその位置を有する。
なお、供給電流の積分が、傾斜磁場コイルのセットの各傾斜磁場コイルに対して各固定反復時間に対して一定であることは、特定の方向における傾斜磁場強度の積分が、各パルス反復時間に対して一定であるように設定することと同等である。これは、傾斜磁場コイル電流と磁場強度とは比例するからである。したがって、ある期間内の供給電流の一定積分は、特定の方向おける又は同じ期間内の当該特定のコイルの磁場強度の積分と同等である。
例えば、パルスシーケンスコマンドが、傾斜磁場コイルのセットへの供給電流を制御することによる少なくとも1つの方向における傾斜磁場の印加を特定し、傾斜磁場コイルのセットの各傾斜磁場コイルに対して、供給電流の積分は、各固定反復時間に対して一定であるとのテキストは、パルスシーケンスコマンドが、少なくとも1つの方向それぞれにおける経時的な傾斜磁場強度の積分が、各固定反復時間又は各基本時間単位に対して一定であるように、傾斜磁場コイルのセットへの供給電流を制御することによる少なくとも1つの方向における傾斜磁場の印加を特定するとのテキストと置き換えられてよい。
パルスシーケンスコマンドは、パルス反復毎に、可変反復時間、可変フリップ角及び可変測定時間で磁気共鳴データの測定を行う命令を含んでよい。これは、優れたサンプリングを提供し、様々なコンポーネントの磁気共鳴フィンガープリンティングディクショナリとのマッチングを可能にするパルス時間の有用な分布を提供する。より詳細には、磁気共鳴フィンガープリンティング技術には、磁気共鳴信号の時間的信号展開を記録することが含まれる。磁気共鳴信号は、RFパルスと傾斜磁場パルスとを含み、1つ又は幾つかの可変要素を有する取得シーケンスによって生成される。これらの可変要素は、取得シーケンスの時間進行に亘って変化する。記録された信号展開は、同じ取得シーケンスを使用してシミュレートされた信号展開と比較される。シミュレーションは、ブロッホ(Bloch)方程式に基づいて行われてよい。測定信号展開のシミュレートされた基準との比較は、ディクショナリアプローチに基づいて行われてよい。MRフィンガープリンティング技術は、組織及び材料パラメータの変化に対する感度が高い。MRフィンガープリンティング技術に関する見識は、様々な材料又は組織に対して、当該材料又は組織を表す一意の信号展開があることである。磁気共鳴フィンガープリンティング技術は、フリップ角、RF位相、反復時間及びk空間サンプリングパターンといった取得パラメータの疑似ランダム変化による時間的及び部分的インコヒーレンスを有する。磁気共鳴フィンガープリンティングは、信号レベルが有限時間量の後、定常状態レベルに達する従来の磁気共鳴撮像方法に比べて、より長い信号展開に沿ってより有益な点をサンプリングすることを可能にする。
RFパルス(フリップ角)のシーケンス、反復時間等は、ランダム又は疑似ランダムであってよい。RFパルスの疑似ランダムシーケンス又は可能なRFパルスの分布から選択されたRFパルスでは、RFパルスのシーケンスは、当該シーケンスが、様々な種に対する潜在MR反応間に最も高度な多様性を達成するようにその符号化パワーを最大にするように選択されてよい。重要な点は、パルスシーケンスが、単一値ではなく、ある範囲の反復時間及びフリップ角を含むことである。これは、結果として得られる磁気共鳴信号が様々な組織に対して異なり、フィンガープリントに似ているようなやり方で選択されてよい。
k空間サンプリングは様々であってよい。例えば1次元における均一k空間サンプリング、1次元における不均一k空間サンプリング及び1次元におけるランダムk空間サンプリングである。zスライス選択といった1次元スライス選択を使用し、x及びy傾斜磁場なしで(即ち、一度に1つのzスライス全体を)サンプリングする場合、k空間における単一点(原点)のみがサンプリングされると言える。z傾斜磁場を、スライス選択ではなく、ここでもx及びy傾斜磁場なしで、z方向におけるk空間のサンプリングに使用することができる。この場合、k空間は1次元であり、サンプリングは、k空間における点の均一又は不均一分布を使用して行われてよい。
別の実施形態では、パルスシーケンスは、一連のパルス反復を含む。一連のパルス反復の各パルス反復は、疑似ランダム分布、持続時間の分布から事前選択された持続時間又は疑似ランダム持続時間を有する。事前選択された持続時間は、結果として得られる一連のRFパルスがランダム又は疑似ランダムに見えるように上記分布から選択されてよいが、他の特性も最適にするように選択されてよい。例えば、上記されたように、RFパルスは、様々な種に対する潜在MR反応間に最も高度な多様性を達成するようにシーケンスの符号化パワーを最大にするように選択される。より具体的には、パルスシーケンスは、RF励起又は信号読出しを有する基本時間単位のランダムな入れ替わりを有する。RF励起を有する乱数の基本時間単位と、信号読出しを有する乱数の基本時間単位とは交互に現れる。別の実施形態では、磁気共鳴システムは、磁気共鳴撮像システムである。測定領域は撮像領域である。傾斜磁場システムは、3つの直交方向において傾斜磁場を生成する。この場合、傾斜磁場コイルの3つのセットが使用される。傾斜磁場システムは、サンプリングイベント中に、上記3つの方向において磁気共鳴データを空間符号化するために、測定領域内に位相符号化傾斜磁場を追加的に生成する。空間符号化は、磁気共鳴データを離散ボクセルに分割する。
主磁場は、しばしば、B0磁場とも呼ばれる。パルスシーケンスコマンドは更に、磁気共鳴データの取得中に、磁気共鳴データの空間符号化を行うように傾斜磁場システムを制御する命令を含む。空間符号化は、磁気共鳴データを離散ボクセルに分割する。この実施形態は、被験者の空間的な結果組成をより迅速に決定する手段を提供できるので有益である。
別の実施形態では、パルスシーケンスコマンドは、位相符号化傾斜磁場が、各サンプリングイベントに関して完全にバランスが取れていることを特定する。位相符号化傾斜磁場が、完全にバランスが取れていると述べることは、総傾斜磁場面積が0であると言っているに等しい。総傾斜磁場面積は、例えば特定の傾斜磁場コイルに供給される電流と解釈されてよい。総傾斜磁場面積は、0であるので、個々の反復の各固定時間又は基本時間単位に対して一定である供給電流の積分に影響を及ぼさない。
別の実施形態では、パルスシーケンスは、スポイル傾斜磁場エコー、任意選択的に疑似T1スポイルシーケンスとして構成される。スポイル傾斜磁場エコーMRFシーケンスは、明白な反復間位相蓄積を有する。すべての横方向磁化は、コヒーレント経路に沿って増大した信号、即ち、FID、励起エコー及び共役励起エコーだけがサンプリングされるように、各反復内でスポイルされる。疑似T1スポイルMRFシーケンスは、明確で有限の数の寄与コヒーレンスを有する。これは、主磁場不均一性に無感応であるサンプリングされた信号を増強し、したがって、MRFライブラリエントリが主磁場変化に依存するものとして列挙される必要がないため、MRF符号化空間は縮小される。
更に、傾斜磁場スポイルされ、任意選択の疑似T1スポイルパルスシーケンスとして形成される本発明のMRFシーケンスは、フリップ角及びシーケンス反復時間の変化を使用する。正味傾斜磁場面積(時間積分)は、シーケンス反復時間と比例して設定される。デフェージングは、連続デフェージング状態間に生じる。即ち、連続(セットの)反復間に十分なデフェージングがある。シーケンス反復時間は、個々の反復の(ベース)固定反復時間又は基本時間単位の整数として選択される。RF励起及び信号サンプリングのないダミー反復が、シーケンス反復時間を更に変更させるように挿入されてもよい。本発明のこれらの態様は、次のことを含む:
−ブロッホ(Bloch)シミュレーションの著しい単純化を可能にする「疑似T1−」スポイルMRFシーケンスの使用(MRF信号に非常に明確な方法で寄与する明白で有限の数のコヒーレンスのみが検出されること見なすことができる)。
−必要なMRF符号化空間を縮小することによって、MRF取得をより効率的にする(オフ共鳴を符号化する必要はない)。結果は、オフ共鳴とは無関係である。
−RFパルスの位相を使用して、専用「RFスポイリング」スキームを使用する追加の符号化要素となることを可能にし、効率的なMRFサンプリング及び信頼できる信号マッチングを可能にする。
別の実施形態では、空間符号化は、1次元である。離散ボクセルは、離散スライスのセットである。上記方法は更に、磁気共鳴データを、離散スライスのセットに分割するステップを含む。所定物質のセットの各物質の存在量は、離散スライスのセットの各セット内で、離散スライスのセットの各スライスの磁気共鳴データを、磁気共鳴フィンガープリンティングディクショナリと比較することによって計算される。
別の実施形態では、空間符号化は、パルスシーケンスの実行中に、所定方向において一定の傾斜磁場を生成するように、傾斜磁場システムを制御することによって行われる。
別の実施形態では、空間符号化は、少なくとも部分的にサンプリングイベント中に、1次元読出し傾斜磁場を生成するように、傾斜磁場システムを制御することによって行われる。
別の実施形態では、空間符号化は、3次元である。空間符号化は、少なくとも部分的にサンプリングイベント中に、3次元読出し傾斜磁場を生成するように、傾斜磁場システムを制御することによって行われる。
別の実施形態では、空間符号化は、非デカルト空間符号化として行われる。空間符号化は、非デカルト順序(non-Cartesian order)でk空間をサンプリングするサンプリングイベント中に、読出し傾斜磁場を生成するように、傾斜磁場システムを制御することによって行われる。
別の実施形態では、離散ボクセルそれぞれの磁気共鳴データを磁気共鳴フィンガープリンティングディクショナリと比較することによる離散ボクセルそれぞれ内の所定の組織タイプそれぞれの存在量の計算は、最初に、磁気共鳴データの各磁気共鳴信号を、所定物質のセットの各物質からの信号の一次結合として表し、次に、最少化手法を使用して一次結合を解くことによって、所定物質のセットの各物質の存在量を決定することによって行われる。
別の実施形態では、磁気共鳴システムは、核磁気共鳴スペクトロメータである。これは、NMRスペクトロメータとしても知られている。
別の実施形態では、マシン実行可能命令の実行は更に、プロセッサに、磁気共鳴フィンガープリンティングディクショナリを計算させる。磁気共鳴フィンガープリンティングディクショナリの実際の計算は、NMR信号をモデル化する様々な技術のどの技術で行われてもよい。例えばNMR信号は、いわゆるブロッホ(Bloch)方程式を使用して計算される多数の単一スピンを合計することによってモデル化されてよい。ディクショナリは、物質パラメータの具体的なセットとパルスシーケンスコマンドによって特定される特定のMRシーケンスとに対するボクセルからの期待NMR信号を計算することによって作成される。
別の実施形態では、磁気共鳴フィンガープリンティングディクショナリは、拡張位相グラフ(Extended phase graph)形式を使用して、所定物質それぞれをモデル化することによって計算される。拡張位相グラフ形式は、例えばWeigel,M.(2015)、Extended phase graphs: Dephasing, RF pulses, and echoes-pure and simple(J.Magn.Reson.Imaging、41:266−295、doi:10.1002/jmri.24619)に説明され、また、Scheffler,K.(1999)、A pictorial description of steady-states in rapid magnetic resonance imaging(Concepts Magn.Reson.、11:291−304、doi:10.1002/(SICI)1099−0534(1999)11:5<291::AID−CMR2>3.0.CO;2−J)でも説明されている。
別の実施形態では、パルスシーケンスコマンドは、固定遅延におけるk空間中心の読出しを特定する。
別の実施形態では、命令の実行は更に、プロセッサに、少なくとも1つの較正ファントムの磁気共鳴データの測定を繰り返させる。少なくとも1つの較正ファントムは、所定物質のセットの少なくとも1つの物質の既知のボリュームを含む。
1次元に沿って磁気共鳴データを測定するシステムと共に使用される場合、較正ファントムそれぞれが較正軸を有してよい。この場合、少なくとも1つの較正ファントムは、較正軸が所定方向と一直線になると、所定物質のセットの少なくとも1つの物質の既知のボリュームを含む。例えば較正ファントムが3次元又は2次元撮像が行われるシステム内で使用される他の場合では、所定物質は、較正ファントム内に既知の集結で均一に分布してよい。
別の態様では、本発明は、測定領域内の被験者からの磁気共鳴データを取得する磁気共鳴システムを制御するプロセッサによる実行のためのマシン実行可能命令を含むコンピュータプログラムプロダクトを提供する。磁気共鳴システムは、測定領域内に主磁場を生成する磁石を含む。磁気共鳴システムは更に、少なくとも1つの方向それぞれの傾斜磁場コイルのセットに電流を供給することによって、上記少なくとも1つの方向において、測定領域内に傾斜磁場を生成する傾斜磁場システムを含む。マシン実行可能命令の実行は、プロセッサに、パルスシーケンスコマンドを用いて磁気共鳴システムを制御することによって磁気共鳴データを取得させる。
パルスシーケンスコマンドは、磁気共鳴システムに、磁気共鳴フィンガープリンティング技術に従って磁気共鳴データを取得させる。パルスシーケンスコマンドは、一連のパルスシーケンス反復を特定する。各パルスシーケンス反復は、固定の反復時間を有する。各パルスシーケンス反復は、パルスシーケンス反復の始まりから固定遅延において生じる無線周波数パルスか又はサンプリングイベントを含む。無線周波数パルスは、無線周波数パルスの分布から選択される。無線周波数パルスの分布は、フリップ角の分布に向けて磁気スピンを回転させる。パルスシーケンスコマンドは、傾斜磁場コイルのセットに供給される電流を制御することによる少なくとも1つの方向における傾斜磁場の印加を特定する。傾斜磁場コイルのセットの各傾斜磁場コイルに対して、供給電流の積分は、各固定反復時間に対して一定である。
命令の実行は更に、プロセッサに、磁気共鳴データを磁気共鳴フィンガープリンティングディクショナリと比較することによって、所定物質のセットの各物質の存在量を計算させる。磁気共鳴フィンガープリンティングディクショナリは、計算した磁気共鳴信号のリストを含み、また、所定物質のセットに対するパルスシーケンスコマンドの実行に反応する。
別の態様では、本発明は、測定領域内の被験者からの磁気共鳴データを取得するように磁気共鳴システムを作動させる方法を提供する。磁気共鳴システムは、測定領域内に主磁場を生成する磁石を含む。磁気共鳴システムは更に、少なくとも1つの方向それぞれの傾斜磁場コイルのセットに電流を供給することによって、少なくとも1つの方向において、測定領域内に傾斜磁場を生成する傾斜磁場システムを含む。上記方法は、パルスシーケンスコマンドを用いて磁気共鳴システムを制御することによって磁気共鳴データを取得するステップを含む。
パルスシーケンスコマンドは、磁気共鳴システムに、磁気共鳴フィンガープリンティング技術に従って磁気共鳴データを取得させる。パルスシーケンスコマンドは、一連のパルスシーケンス反復を特定する。各パルスシーケンス反復は、固定反復時間を有する。各パルスシーケンス反復は、パルスシーケンス反復の始まりから固定遅延において生じる無線周波数パルスか又はサンプリングイベントを含む。無線周波数パルスは、無線周波数パルスの分布から選択される。無線周波数パルスの分布は、フリップ角の分布に向けて磁気スピンを回転させる。パルスシーケンスコマンドは、傾斜磁場コイルのセットに供給される電流を制御することによる少なくとも1つの方向における傾斜磁場の印加を特定する。傾斜磁場コイルのセットの各傾斜磁場コイルに対して、供給電流の積分は、各固定反復時間に対して一定である。
上記方法は更に、磁気共鳴データを磁気共鳴フィンガープリンティングディクショナリと比較することによって、所定物質のセットの各物質の存在量を計算するステップを含む。磁気共鳴フィンガープリンティングディクショナリは、所定物質のセットに対するパルスシーケンスコマンドの実行に反応して計算された磁気共鳴信号のリストを含む。
なお、本発明の上記実施形態の1つ以上は、組み合わせられた実施形態が相互に排他的ではない限り、組み合わせられてよい。
以下において、本発明の好適な実施形態が、ほんの一例として、図面を参照して説明される。
図1は、磁気共鳴撮像システムの一例を示す。 図2は、図1の磁気共鳴撮像システムを作動させる方法を説明するフローチャートを示す。 図3は、パルスシーケンスの一部を示す。 図4は、パルスシーケンスの更なる一部を示す。 図5は、一連のパルスシーケンス反復を示す。 図6は、図5に示されるパルスシーケンスの位相グラフを示す。 図7は、図5及び図6に示されるパルスシーケンスの代替表現を示す。 図8は、磁気共鳴フィンガープリンティングディクショナリの一例を示す。 図9は、磁気共鳴フィンガープリンティングディクショナリの一例を示す。
これらの図面内の同様の参照符号を有する要素は、同等の要素であるか又は同じ機能を行う。先に説明された要素は、機能が同等であるならば、後の図面において必ずしも説明されない。
図1は、磁石104を有する磁気共鳴撮像システム100の一例を示す。磁石104は、中にボア106を有する円筒型の超電導磁石104である。異なるタイプの磁石の使用も可能である。例えばスプリット円筒磁石といわゆるオープン磁石との両方を使用することも可能である。スプリット円筒磁石は、クライオスタットが二分されて、磁石のイソプレーン(iso-plane)へのアクセスが可能にされている点を除いて標準的な円筒磁石と同様であり、このような磁石は、例えば荷電粒子線療法と併せて使用される。オープン磁石は、被検者を受容するのに十分大きい空間を間に有して上下にある2つの磁石部を有する。これらの2つの磁石部の配置は、ヘルムホルツ(Helmholtz)コイルのそれに似ている。オープン磁石は、被検者の閉塞性がより少ないため、人気がある。円筒磁石のクライオスタット内に、超電導コイルの集合がある。円筒磁石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は、複数の送信/受信チャネルを有してもよい。
被験者支持体120が、被験者支持体及び被験者118を撮像領域108内で動かすことができるオプションのアクチュエータ122に取り付けられる。このようにすると、被験者118の大部分又は被験者118全体を撮像することができる。トランシーバ116、傾斜磁場電源112及びアクチュエータ122は、すべて、コンピュータシステム126のハードウェアインターフェース128に接続されていることが見て取れる。コンピュータストレージ134は、磁気共鳴フィンガープリンティング技術を行うためのパルスシーケンスコマンド140を含むものとして示されている。
パルスシーケンスコマンドは、磁気共鳴システムに、磁気共鳴フィンガープリンティング技術に従って磁気共鳴データを取得させる。パルスシーケンスコマンドは、一連のパルスシーケンス反復を特定する。各パルスシーケンス反復は、固定反復時間を有する。各パルスシーケンス反復は、パルスシーケンス反復の始まりから固定遅延において生じる無線周波数パルスか又はサンプリングイベントを含む。無線周波数パルスは、無線周波数パルスの分布から選択される。無線周波数パルスの分布は、フリップ角の分布に向けて磁気スピンを回転させる。パルスシーケンスコマンドは、傾斜磁場コイルのセットへの供給電流を制御することによる少なくとも1つの方向における傾斜磁場の印加を特定する。傾斜磁場コイルのセットの各傾斜磁場コイルに対して、供給電流の積分は、各固定反復時間に対して一定である。
コンピュータストレージ134は更に、磁気共鳴撮像システム100を制御するように、パルスシーケンスコマンド140を使用して取得された磁気共鳴データ142を含むものとして示されている。コンピュータストレージ134は更に、磁気共鳴フィンガープリンティングディクショナリ144を含むものとして示されている。コンピュータストレージは更に、磁気共鳴データ142と磁気共鳴フィンガープリンティングディクショナリ144とを使用して再構成された磁気共鳴画像146を含むものとして示されている。
コンピュータメモリ136は、プロセッサ130が磁気共鳴撮像システム100の動作及び機能を制御することができるようにするオペレーティングシステム又は他の命令のようなコードを含む制御モジュール150を含む。
コンピュータメモリ136は更に、磁気共鳴フィンガープリントディクショナリ生成モジュール152を含むものとして示されている。フィンガープリント生成モジュール152は、磁気共鳴フィンガープリンティングディクショナリ144を作成するために、各ボクセルに対しブロッホ(Bloch)方程式を使用して1つ以上のスピンをモデル化する。コンピュータメモリ136は更に、磁気共鳴データ142と磁気共鳴フィンガープリンティングディクショナリ144とを使用して磁気共鳴画像146を再構成する画像再構成モジュールを含むものとして示されている。例えば磁気共鳴画像146は、被験者118内の所定物質の1つ以上の物質の空間分布のレンダリングであってよい。
図1の例は、磁気共鳴撮像システム又は装置100が核磁気共鳴(NMR)スペクトロメータとなるように変更されてよい。傾斜磁場コイル110及び傾斜磁場コイル電源112を用いずに、装置100は撮像ゾーン108内で0次元測定を行う。
図2は、図1の磁気共鳴撮像システム100を作動させる方法を説明するフローチャートを示す。まず、ステップ200において、パルスシーケンスコマンド140を用いて磁気共鳴撮像システムを制御することによって、磁気共鳴データ142が取得される。次に、ステップ202において、磁気共鳴データ142を磁気共鳴フィンガープリンティングディクショナリ144と比較することによって、所定物質のセットの各物質の存在量が計算される。存在量は、例えば磁気共鳴画像146内にプロット又は表示される。
MRフィンガープリンティングは、定量的MRIへの将来有望な新しいアプローチである。本発明開示は、MRフィンガープリンティングに必要な柔軟性を可能にするが、ΔB0作用(B0、即ち、主磁場の不均一性による作用)とは無関係のMR信号をもたらす新規のMRシーケンスのクラスについて説明する。これらのシーケンスは、ボクセル内デフェージングによる信号損失といった問題や、MRフィンガープリンティングディクショナリのシミュレーションにΔB0を不必要な余剰の次元として含めなければならないことを回避することができる。説明されるシーケンスの大きな利点は、これらのシーケンスが、180度パルスを必要とせず、これによりSAR問題が回避される点にある。拡張位相グラフ(Extended phase graph:EPG)形式を使用して、期待される信号を容易にシミュレートすることができ、これは、適度な時間内に大きいディクショナリを計算することを可能にする。
MRフィンガープリンティング(MRF)は、定量的MRIへの新しいアプローチである。MR撮像への従来のアプローチは、1つの画像を再構成するために必要なすべてのデータを取得するために、様々な位相符号化ステップに対して同じ基本シーケンス構成単位を何回も繰り返すことに基づいている。対照的に、MRフィンガープリンティングは、多くの可変要素(フリップ角、TR、…)を含むシーケンスを使用して、MR信号の時間的展開を記録する。次に、取得信号を、同じシーケンスを使用してシミュレートされた信号展開と比較することによって、撮像物体の磁化特性(M0、T1、T2、…)及び信号展開に影響を及ぼすシステムパラメータが取得される。通常、シミュレーションは、大きいセットのパラメータ組み合わせの実験の前に行われ、ディクショナリに格納される。測定が行われると、物体の特性は、取得データにベストマッチする信号展開をディクショナリから見つけることによって取得される。
MRフィンガープリンティングのシーケンスは、マッチング処理における臨床的問題に関連する組織/材料のパラメータの精度を高めるために、当該シーケンスが、これらの変数の変化に対する感度が高いように選択される必要がある。その一方で、シーケンスは、MR信号に潜在的に影響を及ぼすすべての他の要素をディクショナリにおける追加の次元として含める(指数関数的に遅いマッチング/シミュレーションとマッチングにおける潜在的な曖昧さとにつながる)ことを回避するために、これらの変数に対する感度は高くあるべきではない。
通常、臨床上全く重要ではないが、MR信号に好ましくない影響を及ぼす1つのパラメータは、オフ共鳴(ΔB0)である。実施例は、MRフィンガープリンティングに必要な柔軟性を可能にするが、ΔB0作用のないMR信号をもたらす新しいMRシーケンスの新しいクラスを提供する。これらのシーケンスは、ボクセル内デフェージングによる信号損失といった問題や、MRフィンガープリンティングディクショナリのシミュレーションにΔB0を不必要な余剰の次元として含めなければならないことを回避する。
MRフィンガープリンティングシーケンスの例は、場合によっては、シーケンスオブジェクトに特定の制限を課すことによって、主磁場の変化に対して無感応であるようにされてよい。しかし、これらの制限は、MRフィンガープリンティングに必要な変化を達成するために、十分な柔軟性を残す。
Tを、シーケンスの選択された基本時間単位の持続時間とする。シーケンスは、次の特徴のうちの1つ以上を有する。
1.RFパルスはTの整数倍の時間にのみ置かれる。RFパルスのフリップ角及び位相は任意に選択される。
2.3つの主方向のすべてにおいて、Tの各間隔中に同じ総傾斜磁場面積が独立して蓄積されなければならない。(この要件は、位相符号化傾斜磁場が各データ取得近くで完全にバランスが取れていることを条件に、位相符号化傾斜磁場には適用されない。)
3.データ取得期間は、Tの整数倍近くに、その時間にRFパルスがない場合は置かれる。(k空間中心k=0は、常に、ちょうどTの整数倍において通過される。)
長さ及び傾斜磁場面積の要件は、シーケンスによって生成されたすべての磁化状態が、EPG形式で等間隔の状態の単一のセットによってモデル化可能であることを確実にし、これは、磁化の効率的なリフォーカシングを確実にする。これは更に、オフ共鳴による追加の位相が、Tの整数倍において、即ち、各読出しの中心において常にちょうどゼロであることも保証する。
Tのすべての整数倍においてエコーが形成されるので、RFパルスによって使用されていないTのすべての整数倍においてデータ取得期間を置いて、物体に関して取得される情報の割合を最大限にすることが有利である。即ち、フィンガープリンティングシーケンスを、励起(E)及び取得(A)のシーケンスによって特徴付けることができる。これらの2文字のすべてのシーケンスは、オフ共鳴に無感応である有効なフィンガープリンティングシーケンスを表現する(各Eにおいて、フリップ角及び位相を選択する追加の自由がある)。
以下において、可能なシーケンスの一例が、エコーの形成を視覚化する対応する位相グラフと共に示される。以下に説明される図3及び図4は、Eセグメント(図3)及びAセグメント(図4)の象徴的表現を示す。
図3は、パルスシーケンス300の一部を示す。横棒302は、固定反復時間302を表す。3本の線があり、第1の線304は、読出し傾斜磁場を特定するために使用される。線306は、位相符号化傾斜磁場306を特定するための空間を示し、線308は、スライス選択傾斜磁場を特定するための空間と、更に、RFパルス310の位置とを示す。図3に示される例では、パルスシーケンス300は、位相符号化を示さない。単一の読出し傾斜磁場312が示されている。上記例では、線308上に、固定反復時間302の始まり及び終わりの両方において生じるスライス選択傾斜磁場314とRFパルス310とがある。この例では、中心線316があり、スライス選択傾斜磁場314及びRFパルス310はこの線を中心とする。RFパルス310及びスライス選択傾斜磁場314は、固定反復時間302の始まり及び終わりによって部分に分割されているものとして示されている。しかし、これは少し人工的であり、というのは、固定反復時間302の始まり及び終わりはずらすことが可能であり、中心線316の位置を、パルスシーケンス反復の始まりからの固定遅延と解釈することができる。つまり、図3は、RFパルス310及びスライス選択傾斜磁場314全体が単一の固定反復時間302内にあるように再構成されてもよい。図3は、RFパルスが遅延316において印加されるパルスシーケンス反復を表す。
図4は、パルスシーケンス400の一部の更なる例を示す。図4は、サンプリングイベント404が生じるパルスシーケンス反復を示す。固定反復時間302の始まりにおいて、RFパルス310がある。固定反復時間302の終わりにはRFパルスがない。代わりに、固定遅延316’に関して対称に印加される読出し傾斜磁場312がある。位相符号化傾斜磁場402もあり、これも、固定遅延316’に関して対称である。同様に、固定遅延316’の近くで、スライス選択傾斜磁場314’が2つの対称部分に分割されている。図3と同様に、固定反復時間302の正確な位置は、特定のパルスシーケンス反復のすべてのコンポーネントが固定反復時間302内にあるように調整される。固定遅延316及び316’は共に、表示されている固定反復時間302の始まり及び終わりにおいて示されている。しかし、例えば固定反復時間302の始まりは、ちょうど316と316’との間となるようにずらされてもよい。この場合、傾斜磁場312、402及び314’をすべて同じ固定反復時間302内に含めることができる。
図3及び図4の両方において、1つのTの間のM304、P306及びS308における総面積は、それぞれ、Am、0、Asである。読出し傾斜磁場の総面積は2Amである。
これらの要素を組み合わせることによって(ここでもA、Eフラグメントを使用して)、より長いシーケンスを作成することができる。
図5は、一連のパルスシーケンス反復500を示す。この場合、タイムラインは、等しい持続時間を有する幾つかの部分に分割されている。これらの部分は、固定反復時間に相当する。固定反復時間は、参照符号502か又は504が付けられる。固定反復時間502は、固定遅延316に無線周波数パルスを有するパルスシーケンス反復に対応する。パルスシーケンス反復504は、固定遅延316を中心にサンプリングイベント404を有するパルスシーケンス反復に対応する。図5における期間は、図3及び図4とは異なるように分割されている。傾斜磁場312、402、314、314’はすべて、対応するパルスシーケンス反復502又は504内に含まれる特定のパルスシーケンス反復に対するものであることが見て取れる。図5は、磁気共鳴フィンガープリンティングを行うのに有用であるパルスシーケンスコマンドを形成するように、パルスシーケンス反復502又は504の基本構成単位をつなぎ合わせるために使用する方法を示す。
図6は、図5に示されるパルスシーケンス500の位相グラフを示す。RFパルス310及びサンプリングイベント404の位置が示されている。図5に示されるパルスシーケンスは、ランダムシーケンスを表す。各時間単位又は固定反復時間の間、同じ傾斜磁場面積が蓄積される。RFパルスは、任意のフリップ角及び位相を有し、固定反復時間の整数倍に置かれてよい。この例では、RFパルスは、固定遅延316に置かれている。次に、固定反復時間のすべての整数倍においてエコーが生成される。エコーは、固定遅延316において生成される。固定遅延316においてRFパルスがない場合、エコーが読み出されてよい。斜めの線600は、位相グラフベースでのスピン系の状態を表す。すべての状態が示されているわけではなく、一部の状態の展開は、図の上縁及び下縁において切り取られている。線の傾きは不平衡の傾斜磁場による位相の取得を表す。
図7は、図5及び図6に示されているパルスシーケンスの別の表現を示す。図7は、2つのプロットを示す。上のプロット700は、選択されたフリップ角の分布をプロットしたものであり、第2のプロット、即ち、下のプロットは、次の無線周波数パルスが印加されるまで待機するステップの数702を示す。x軸704は、パルスシーケンス反復の数704である。上のプロットにおけるy軸706は、選択されたフリップ角を示す。下のy軸708は、単位ブロックの数708を示す。単位ブロックの数は、パルスシーケンス反復に相当する。
このシーケンスは、50のステップからなる(x軸上に示される)。各ステップは、時間Tの単位ブロックの整数を含む(下のグラフに示される)。各ステップの始まりにおいて、ランダムフリップ角のRFパルスが印加される(上のグラフ)。
長さnTを有するステップは、1つのRFパルスと(n−1)の測定結果とを含む。したがって、結果として得られるフィンガープリント信号の長さは、シーケンスステップの数とは異なる。次のグラフは、様々なT1/T2組み合わせについて上記シーケンスから計算される2つの異なるフィンガープリント信号の一例を示す。このような計算信号のセットは、測定結果との比較のためのMRFディクショナリとして使用することができる。
図8は、磁気共鳴フィンガープリンティングディクショナリ800の一例を示す。第1のエントリ802と第2のエントリ804とがある。この例では、ディクショナリ800は、パルスシーケンスコマンドの特定のセットについて計算される。各エントリ802、804は、2つの異なる材料について測定される期待MR信号を表す。材料802は、400msのT1時間と100msのT2時間とを有する。材料804は、1000msのT1時間と500msのT2時間とを有する。実測MR信号を、2つのディクショナリエントリ802、804と比較することができ、例えば当該2つのエントリの一次結合を、測定MR信号に近似させるために追加されてよい。このようにすると、特定ボリューム内の第1の材料802と第2の材料804との相対比が推測される。
本発明の他の態様は、傾斜磁場及び任意選択のRFスポイリングを含むスポイル傾斜磁場エコーシーケンスを使用して、MRフィンガープリンティングを行う。スポイルシーケンスは、(オフ共鳴及び傾斜磁場切り替えによる)明白なTR間位相蓄積と、T1重み付けを達成するために任意選択的に使用される適切なRF信号スポイリングとによって特徴付けられる。すべての横方向磁化は、各TR内でスポイルされるので、離散コヒーレンス経路からの信号(基本的に、FID、スピンエコー及び励起エコー)のみがコヒーレントに重ね合わされ、測定MR信号に寄与する。具体的には、オフ共鳴作用は、選択された傾斜磁場エコー時間に固有のT2コントラストにまで減少される。更に、非常に多くの磁気モーメントからの寄与を合計するのではなく、可算個のコヒーレンスのみを追跡すればよいので、ディクショナリの計算は非常に簡易化される。
MRF取得において、一般に、シーケンスのフリップ角及びTRが、選択されたディクショナリによる意図的な選択によって変化するので、定常状態は増大しない。可変フリップ角αは、構造理論のフレームワークにおいて既に取り扱われているが、可変TRはすぐに可能ではない。各TR間隔の終わりに、単純な可変遅延を追加することは、切り替えられた傾斜磁場及び静的傾斜磁場(即ち、オフ共鳴)からの非同期位相寄与につながり、これは、計算(ブロッホ(Bloch)シミュレーション)を非常に困難にする。したがって、第1の要件として、切り替えられた傾斜磁場の正味傾斜磁場面積が、対応するTRに比例するように調整されなければならない。第2の要件は、様々なデフェージング状態を混合することが可能ではないという事実から生じる。つまり、隣接するデフェージング状態間には十分なデフェージングが常になければならず、これは、状態1=0のみが測定信号に寄与することが確実にされる。そうでなければ、ずれた「エコートップ」を有するエコーが重なり合い、深刻なアーチファクトを引き起こし、MRF信号受信を損なう。これは、可変TRを元のシーケンスのベースTRの整数倍として選択することによって達成される。実用的な実施態様では、両要件は、RF励起及び取得の両方がオフにされている「ダミーシーケンスモジュール」を取り入れることによって満たされる(図9を参照)。
本発明は、図面及び上記説明において詳細に例示及び説明されたが、このような例示及び説明は、例示的であって、限定と考えられるべきではない。本発明は、開示された実施形態に限定されない。
開示された実施形態の他の変形態様は、図面、開示内容及び添付の請求項の検討から、請求項に係る発明を実施する当業者によって理解され、実施される。請求項において、「含む」との用語は、他の要素又はステップを排除するものではなく、また、「a」又は「an」との不定冠詞も、複数形を排除するものではない。単一のプロセッサ又は他のユニットが、請求項に記載される幾つかのアイテムの機能を果たしてもよい。特定の手段が相互に異なる従属請求項に記載されることだけで、これらの手段の組み合わせを有利に使用することができないことを示すものではない。コンピュータプログラムは、他のハードウェアと共に又はその一部として供給される光学ストレージ媒体又は固体媒体といった適切な媒体上に記憶及び/又は分散されてもよいが、インターネット又は他の有線若しくは無線通信システムを介するといった他の形式で分配されてもよい。請求項における任意の参照符号は、範囲を限定するものと解釈されるべきではない。
100 磁気共鳴システム
104 磁石
106 磁石のボア
108 測定領域又は撮像領域
110 傾斜磁場コイル
112 傾斜磁場コイル電源
114 無線周波数コイル
116 トランシーバ
118 被験者
120 被験者支持体
122 アクチュエータ
124 所定方向
125 スライス
126 コンピュータシステム
128 ハードウェアインターフェース
130 プロセッサ
132 ユーザインターフェース
134 コンピュータストレージ
136 コンピュータメモリ
140 パルスシーケンスコマンド
142 磁気共鳴データ
144 磁気共鳴フィンガープリンティングディクショナリ
146 磁気共鳴画像
150 制御モジュール
152 磁気共鳴フィンガープリントディクショナリ生成モジュール
154 画像再構成モジュール
200 パルスシーケンスコマンドを用いて磁気共鳴撮像システムを制御することによって、磁気共鳴データを取得する
202 磁気共鳴データを磁気共鳴フィンガープリンティングディクショナリと比較することによって、所定物質のセットの各物質の存在量を計算する
300 パルスシーケンスの一部
302 固定反復時間
304 読出し
306 位相符号化
308 スライス選択
310 RFパルス
312 読出し傾斜磁場
314 スライス選択傾斜磁場
314’ スライス選択傾斜磁場
316 固定遅延
316’ 固定遅延
400 パルスシーケンスの一部
402 位相符号化傾斜磁場
404 サンプリングイベント
500 一連のパルスシーケンス反復
502 RFパルスを有するパルスシーケンス反復
504 サンプリングイベントを有するパルスシーケンス反復
600 状態の組み合わせ
700 選択されたフリップ角の分布
702 次のRFパルスまで待機するステップ数
704 パルスシーケンス反復の数
706 選択されたフリップ角
708 単位ブロックの数
800 磁気共鳴フィンガープリンティングディクショナリ
802 第1のエントリ
804 第2のエントリ

Claims (15)

  1. 測定領域内の被験者からの磁気共鳴データを取得する磁気共鳴システムであって、前記磁気共鳴システムは、
    前記測定領域内に主磁場を生成する磁石と、
    少なくとも1つの方向それぞれの傾斜磁場コイルのセットに電流を供給することによって、前記少なくとも1つの方向において、前記測定領域内に傾斜磁場を生成する傾斜磁場システムと、
    マシン実行可能命令とパルスシーケンスコマンドとを格納するメモリであって、前記パルスシーケンスコマンドは、前記磁気共鳴システムに、磁気共鳴フィンガープリンティング技術に従って前記磁気共鳴データを取得させ、前記パルスシーケンスコマンドは、一連のパルスシーケンス反復を特定し、各パルスシーケンス反復は、固定反復時間を有し、各パルスシーケンス反復は、前記パルスシーケンス反復の始まりから固定遅延において生じる無線周波数パルスか又はサンプリングイベントを含み、前記無線周波数パルスは、無線周波数パルスの分布から選択され、前記無線周波数パルスの分布は、フリップ角の分布に向けて磁気スピンを回転させ、前記パルスシーケンスコマンドは、前記傾斜磁場コイルのセットに供給される前記電流を制御することによる前記少なくとも1つの方向における傾斜磁場の印加を特定し、前記傾斜磁場コイルのセットの各傾斜磁場コイルに対して、供給電流の積分は、各固定反復時間に対して一定である、前記メモリと、
    前記磁気共鳴システムを制御するプロセッサと、
    を含み、
    前記マシン実行可能命令の実行は、前記プロセッサに、
    パルスシーケンスコマンドを用いて前記磁気共鳴システムを制御することによって、前記磁気共鳴データを取得させ、
    前記磁気共鳴データを磁気共鳴フィンガープリンティングディクショナリと比較することによって、所定物質のセットの各物質の存在量を計算させ、
    前記磁気共鳴フィンガープリンティングディクショナリは、所定物質のセットの前記パルスシーケンスコマンドの実行に反応して計算された磁気共鳴信号のリストを含む、磁気共鳴システム。
  2. 前記磁気共鳴システムは、磁気共鳴撮像システムであり、前記測定領域は、撮像領域であり、前記傾斜磁場システムは、3つの直交方向において前記傾斜磁場を生成し、前記傾斜磁場システムは、前記サンプリングイベント中に、前記3つの直交方向において前記磁気共鳴データを空間符号化するように、前記測定領域内に位相符号化傾斜磁場を追加的に生成し、前記空間符号化は、前記磁気共鳴データを離散ボクセルに分割する、請求項1に記載の磁気共鳴システム。
  3. 前記パルスシーケンスコマンドは、前記位相符号化傾斜磁場が、各サンプリングイベントに関して完全にバランスが取れていることを特定する、請求項2に記載の磁気共鳴システム。
  4. 前記空間符号化は、1次元であり、前記離散ボクセルは、離散スライスのセットであり、更に、前記磁気共鳴データを前記離散スライスのセットに分割し、所定物質のセットの各物質の存在量は、前記離散スライスのセットの各スライスにおいて、前記離散スライスのセットの各スライスの前記磁気共鳴データを、前記磁気共鳴フィンガープリンティングディクショナリと比較することによって計算される、請求項2又は3に記載の磁気共鳴システム。
  5. 前記空間符号化は、前記パルスシーケンスコマンドの実行中に、所定方向において一定傾斜磁場を生成するように前記傾斜磁場システムを制御することによって行われる、請求項2又は3に記載の磁気共鳴システム。
  6. 前記空間符号化は、少なくとも部分的に前記サンプリングイベント中に、1次元読出し傾斜磁場を生成するように前記傾斜磁場システムを制御することによって行われる、請求項2又は3に記載の磁気共鳴システム。
  7. 前記空間符号化は、3次元であり、前記空間符号化は、少なくとも部分的に前記サンプリングイベント中に、3次元読出し傾斜磁場を生成するように前記傾斜磁場システムを制御することによって行われる、請求項2又は3に記載の磁気共鳴システム。
  8. 前記空間符号化は、非デカルト空間符号化として行われ、前記空間符号化は、非デカルト順序でk空間をサンプリングする前記サンプリングイベント中に、読出し傾斜磁場を生成するように前記傾斜磁場システムを制御することによって行われる、請求項2又は3に記載の磁気共鳴システム。
  9. 前記離散ボクセルそれぞれの前記磁気共鳴データを、前記磁気共鳴フィンガープリンティングディクショナリと比較することによる前記離散ボクセルそれぞれ内の所定組織タイプそれぞれの存在量の計算は、
    前記磁気共鳴データの各磁気共鳴信号を、前記所定物質のセットの各物質からの前記磁気共鳴信号の一次結合として表し、
    最少化手法を使用して前記一次結合を解くことによって、前記所定物質のセットの各物質の存在量を決定することによって行われる、請求項2又は3に記載の磁気共鳴システム。
  10. 前記固定反復時間に等しい持続時間を有する1つ又は幾つかのダミーシーケンスモジュールが、前記一連のパルスシーケンス反復において適用され、各ダミーシーケンスは、RF励起及びサンプリングイベントがない、請求項1に記載の磁気共鳴システム。
  11. 前記一連のパルスシーケンス反復は、傾斜磁場スポイルされ、任意選択的に、疑似T1−スポイルされたシーケンスを形成する、請求項10に記載の磁気共鳴システム。
  12. 前記マシン実行可能命令の実行は更に、前記プロセッサに、前記磁気共鳴フィンガープリンティングディクショナリを計算させる、請求項1乃至11の何れか一項に記載の磁気共鳴システム。
  13. 前記パルスシーケンスコマンドは、前記固定遅延におけるk空間中心の読出しを特定する、請求項1乃至12の何れか一項に記載の磁気共鳴システム。
  14. 測定領域内の被験者からの磁気共鳴データを取得する磁気共鳴システムを制御するプロセッサによる実行のためのマシン実行可能命令を含むコンピュータプログラムであって、
    前記磁気共鳴システムは、前記測定領域内に主磁場を生成する磁石を含み、前記磁気共鳴システムは更に、少なくとも1つの方向それぞれの傾斜磁場コイルのセットに電流を供給することによって、前記少なくとも1つの方向において、前記測定領域内に傾斜磁場を生成する傾斜磁場システムを含み、
    前記マシン実行可能命令の実行は、前記プロセッサに、
    パルスシーケンスコマンドを用いて前記磁気共鳴システムを制御することによって、前記磁気共鳴データを取得させ、
    前記磁気共鳴データを磁気共鳴フィンガープリンティングディクショナリと比較することによって、所定物質のセットの各物質の存在量を計算させ、
    前記パルスシーケンスコマンドは、前記磁気共鳴システムに、磁気共鳴フィンガープリンティング技術に従って前記磁気共鳴データを取得させ、前記パルスシーケンスコマンドは、一連のパルスシーケンス反復を特定し、各パルスシーケンス反復は、固定反復時間を有し、各パルスシーケンス反復は、前記パルスシーケンス反復の始まりから固定遅延において生じる無線周波数パルスか又はサンプリングイベントを含み、前記無線周波数パルスは、無線周波数パルスの分布から選択され、前記無線周波数パルスの分布は、フリップ角の分布に向けて磁気スピンを回転させ、前記パルスシーケンスコマンドは、前記傾斜磁場コイルのセットに供給される前記電流を制御することによる前記少なくとも1つの方向における傾斜磁場の印加を特定し、前記傾斜磁場コイルのセットの各傾斜磁場コイルに対して、供給電流の積分は、各固定反復時間に対して一定であり、
    前記磁気共鳴フィンガープリンティングディクショナリは、所定物質のセットの前記パルスシーケンスコマンドの実行に反応して計算された磁気共鳴信号のリストを含む、コンピュータプログラム。
  15. 測定領域内の被験者からの磁気共鳴データを取得するように磁気共鳴システムを作動させる方法であって、前記磁気共鳴システムは、前記測定領域内に主磁場を生成する磁石を含み、前記磁気共鳴システムは更に、少なくとも1つの方向それぞれの傾斜磁場コイルのセットに電流を供給することによって、前記少なくとも1つの方向において、前記測定領域内に傾斜磁場を生成する傾斜磁場システムを含み、
    前記方法は、
    パルスシーケンスコマンドを用いて前記磁気共鳴システムを制御することによって前記磁気共鳴データを取得するステップと、
    前記磁気共鳴データを磁気共鳴フィンガープリンティングディクショナリと比較することによって、所定物質のセットの各物質の存在量を計算するステップと、
    を含み、
    前記パルスシーケンスコマンドは、前記磁気共鳴システムに、磁気共鳴フィンガープリンティング技術に従って前記磁気共鳴データを取得させ、前記パルスシーケンスコマンドは、一連のパルスシーケンス反復を特定し、各パルスシーケンス反復は、固定反復時間を有し、各パルスシーケンス反復は、前記パルスシーケンス反復の始まりから固定遅延において生じる無線周波数パルスか又はサンプリングイベントを含み、前記無線周波数パルスは、無線周波数パルスの分布から選択され、前記無線周波数パルスの分布は、フリップ角の分布に向けて磁気スピンを回転させ、前記パルスシーケンスコマンドは、前記傾斜磁場コイルのセットに供給される前記電流を制御することによる前記少なくとも1つの方向における傾斜磁場の印加を特定し、前記傾斜磁場コイルのセットの各傾斜磁場コイルに対して、供給電流の積分は、各固定反復時間に対して一定であり、
    前記磁気共鳴フィンガープリンティングディクショナリは、所定物質のセットの前記パルスシーケンスコマンドの実行に反応して計算された磁気共鳴信号のリストを含む、方法。
JP2017553237A 2015-04-14 2016-04-13 主磁場の不均一性に対する感度が減少された磁気共鳴フィンガープリンティング Expired - Fee Related JP6339299B1 (ja)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
EP15163505.9 2015-04-14
EP15163505 2015-04-14
US201562206327P 2015-08-18 2015-08-18
US62/206,327 2015-08-18
PCT/EP2016/058048 WO2016166119A1 (en) 2015-04-14 2016-04-13 Magnetic resonance fingerprinting with reduced sensitivity to inhomogeneities in the main magnetic field

Publications (2)

Publication Number Publication Date
JP6339299B1 true JP6339299B1 (ja) 2018-06-06
JP2018516622A JP2018516622A (ja) 2018-06-28

Family

ID=52829020

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017553237A Expired - Fee Related JP6339299B1 (ja) 2015-04-14 2016-04-13 主磁場の不均一性に対する感度が減少された磁気共鳴フィンガープリンティング

Country Status (5)

Country Link
US (1) US20180106876A1 (ja)
EP (1) EP3283897A1 (ja)
JP (1) JP6339299B1 (ja)
CN (1) CN107533121A (ja)
WO (1) WO2016166119A1 (ja)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016074946A1 (en) * 2014-11-14 2016-05-19 Koninklijke Philips N.V. Magnetic resonance fingerprinting in slices along a one-dimensional extension
US10241178B2 (en) * 2015-05-29 2019-03-26 Case Western Reserve University System and method for magnetic resonance fingerprinting at high field strengths
JP7027423B2 (ja) * 2016-11-22 2022-03-01 バイエル・ヘルスケア・エルエルシー 造影剤を用いた合成磁気共鳴イメージングおよび磁気共鳴指紋法、ならびに造影剤注入プロトコルの決定
CN106597337B (zh) * 2016-12-09 2019-05-07 深圳先进技术研究院 一种磁共振t2*加权快速成像方法及装置
CN107942272B (zh) * 2018-01-11 2020-08-21 郑州航空工业管理学院 主磁场不均匀下的分数域磁共振成像方法
US11373392B2 (en) * 2018-04-06 2022-06-28 Case Western Reserve University System and method for magnetic resonance fingerprinting using a plurality of pulse sequence types
EP3680680A1 (en) * 2019-01-11 2020-07-15 Siemens Healthcare GmbH Method for obtaining an optimised operating parameter, storage medium and magnetic resonance apparatus
US11131733B2 (en) * 2019-04-25 2021-09-28 Case Western Reserve University System and method for magnetic resonance fingerprinting with non-locally sequential sampling of k-space
DE102019206827A1 (de) 2019-05-10 2020-11-12 Siemens Healthcare Gmbh Magnetresonanz-Fingerprinting-Verfahren

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6570383B1 (en) * 2000-12-27 2003-05-27 Ge Medical Systems Global Technology Company, Llc Method and apparatus for a pulse sequence for magnetic resonance imaging in an inhomogeneous magnetic field
US7046004B2 (en) * 2004-05-07 2006-05-16 University Of Basel Generic eddy-current compensation in magnetic resonance imaging
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
CN104583799B (zh) * 2012-08-27 2018-12-18 皇家飞利浦有限公司 基于快速图像采集的运动跟踪
EP2897523B1 (en) * 2012-09-19 2024-05-22 Case Western Reserve University Nuclear magnetic resonance (nmr) fingerprinting
US9568579B2 (en) * 2012-12-19 2017-02-14 Case Western Reserve University Magnetic resonance fingerprinting (MRF) with echo splitting

Also Published As

Publication number Publication date
CN107533121A (zh) 2018-01-02
US20180106876A1 (en) 2018-04-19
WO2016166119A1 (en) 2016-10-20
EP3283897A1 (en) 2018-02-21
JP2018516622A (ja) 2018-06-28

Similar Documents

Publication Publication Date Title
JP6339299B1 (ja) 主磁場の不均一性に対する感度が減少された磁気共鳴フィンガープリンティング
US10321845B2 (en) Magnetic resonance fingerprinting in slices along a one-dimensional extension
CN107110938B (zh) 使用具有额外180度rf脉冲的自旋回波脉冲序列的磁共振指纹
JP6608564B1 (ja) 勾配インパルス応答関数マッピング
CN109791186B (zh) 对磁共振指纹期间的b0偏共振场的直接测量
CN110073232B (zh) 用于多状态磁共振指纹识别的磁共振成像系统、方法和计算机可读介质
CN107850650B (zh) 使用补充磁场线圈生成磁共振指纹识别词典的方法和仪器
JP2019518498A (ja) サイレント3d磁気共鳴フィンガープリンティング
JP2020515331A (ja) サブボクセル解像度磁気共鳴フィンガープリンティングイメージング
CN108603921B (zh) 稳态磁共振指纹
WO2019149769A1 (en) Mri with fat/water separation
JP6613161B2 (ja) Mriシミュレーションシステム

Legal Events

Date Code Title Description
A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20180316

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: 20180410

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20180509

R150 Certificate of patent or registration of utility model

Ref document number: 6339299

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees