JP2012521843A - 加速b1マッピング - Google Patents

加速b1マッピング Download PDF

Info

Publication number
JP2012521843A
JP2012521843A JP2012502841A JP2012502841A JP2012521843A JP 2012521843 A JP2012521843 A JP 2012521843A JP 2012502841 A JP2012502841 A JP 2012502841A JP 2012502841 A JP2012502841 A JP 2012502841A JP 2012521843 A JP2012521843 A JP 2012521843A
Authority
JP
Japan
Prior art keywords
coils
coil
field mapping
high frequency
frequency transmit
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2012502841A
Other languages
English (en)
Other versions
JP5525596B2 (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
Koninklijke Philips Electronics 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, Koninklijke Philips Electronics NV filed Critical Koninklijke Philips NV
Publication of JP2012521843A publication Critical patent/JP2012521843A/ja
Application granted granted Critical
Publication of JP5525596B2 publication Critical patent/JP5525596B2/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/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/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/5659Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the RF magnetic field, e.g. spatial inhomogeneities of the RF magnetic field
    • 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/58Calibration of imaging systems, e.g. using test probes, Phantoms; Calibration objects or fiducial markers such as active or passive RF coils surrounding an MR active material
    • G01R33/583Calibration of signal excitation or detection systems, e.g. for optimal RF excitation power or frequency
    • 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/24Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
    • G01R33/246Spatial mapping of the RF magnetic field B1
    • 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/58Calibration of imaging systems, e.g. using test probes, Phantoms; Calibration objects or fiducial markers such as active or passive RF coils surrounding an MR active material
    • G01R33/583Calibration of signal excitation or detection systems, e.g. for optimal RF excitation power or frequency
    • G01R33/586Calibration of signal excitation or detection systems, e.g. for optimal RF excitation power or frequency for optimal flip angle of RF pulses

Landscapes

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

Abstract

方法は、B磁場マッピングデータセットを収集するために高周波送信コイルのセット11を用いてある回数のB磁場マッピングシーケンス24を実行するステップであって、前記回数は高周波送信コイルのセットにおける高周波送信コイルの数よりも少ない、ステップと、収集されたB磁場マッピングデータセットに基づいて高周波送信コイルのセットに対するコイル感度30を決定するステップとを有する。一部の実施形態において、実行されるB磁場マッピングシーケンスは、(i)直交仮想高周波送信コイルのセット42を作るために高周波送信コイルのセットに線形変換40を実行するステップと、(ii)実行されるB磁場マッピングシーケンスを規定する直交仮想高周波送信コイルのセットのサブセット46を選択するステップ44とによって規定される。

Description

下記は磁気共鳴法、磁気共鳴画像法、磁気共鳴分光法、医用磁気共鳴画像法及び分光法、並びに関連技術に関する。
磁気共鳴(MR)画質はB送信磁場の均一性に依存する。このような目的で、磁場励起のために使用される高周波コイルは視野(FOV)全体に十分に均一なB送信磁場を生じるように設計される。これは複数の高周波送信コイル又はそのアレイを用いることによってなされることができ、これらは集合的にFOV全体に十分に均一なB送信磁場を生じるようにFOVに各々戦略的に配置される。しかしながらこうした設計は不完全である。さらに、無負荷状態において均一であるように設計されたB送信磁場は被験者の磁化率によって歪められる可能性がある。これは被験者負荷効果として知られ、この効果は約3テスラ以上の静磁場など、高いB静磁場において顕著になる。低い静磁場であっても、被験者負荷効果は無視できない可能性がある。
シミングは無負荷又は負荷状態においてB送信磁場の均一性を向上させるために使用されることができる。付加的に又は代替的に、B送信磁場不均一性は収集後MR画像再構成処理中に数学的に補正されることができる。B送信磁場の空間マップはBシミング及び収集後の数学的補正の両方に対する入力である。被験者負荷効果を考慮すると、B送信磁場は好適にはMRスキャナ及び撮像位置に乗せられる撮像中の特定の被験者とマップされる。別の方法では、B送信磁場はMRスキャナに乗せられる被験者に適切に類似するファントムとマップされることができる。
既存のB送信磁場マッピングMRシーケンスにおいて、測定は二次元又は三次元で個々の送信チャネル又はコイルの各々に対して実行される。Nコイルアレイ(Nはコイルの数を示す)の場合、従ってN回送信磁場マッピングシーケンスが実行される。あいにく、これらのB送信磁場マッピングMRシーケンスは比較的遅く、不都合に撮像セッション時間を延ばす。
本明細書に例として示され記載される特定の実施形態例によれば、方法は、B磁場マッピングデータセットを収集するために高周波送信コイルのセットを用いてある回数のB磁場マッピングシーケンスを実行するステップであって、前記回数は高周波送信コイルのセットにおける高周波送信コイルの数よりも少ない、ステップと、収集されたB磁場マッピングデータセットに基づいて高周波送信コイルのセットに対するコイル感度を決定するステップとを有する。
本明細書に例として示され記載される特定の実施形態例によれば、記憶媒体は直前の段落に記載の方法を実行するためにデジタルプロセッサによって実行可能な命令を保存する。本明細書に例として示され記載される特定の実施形態例によれば、デジタルプロセッサは直前の段落に記載の方法を実行するように構成される。本明細書に例として示され記載される特定の実施形態例によれば、磁気共鳴システムは、磁気共鳴スキャナと、高周波送信コイルのセットと、直前の段落に記載の方法を実行するように構成されるデジタルプロセッサとを有し、ある回数のB磁場マッピングシーケンスを実行するステップは、磁気共鳴スキャナに前記高周波送信コイルのセットを用いて前記B磁場マッピングシーケンスを実行させるステップを有する。
1つの利点はより効率的なBシミングを提供することにある。
別の利点はより効率的なB送信磁場マッピングを提供することにある。
さらなる利点は以下の詳細な説明を読んで理解することで当業者によって理解されるだろう。
図面は好適な実施形態を例示する目的に過ぎず、本発明を限定するものと解釈されない。対応する参照数字は様々な図面において使用されるとき、図中の対応する要素をあらわす。
高周波送信コイルのセット及びBシミングモジュールを含む又はそれらとともに機能する磁気共鳴システムを図示する。 マッピングシーケンスの減少した回数を決定するように構成される図1のBシミングモジュールのプロセッサによって適切に実行されるプロセスをフローチャートに描く。 磁場マッピングデータセットを収集し、収集されたB磁場マッピングデータセットに基づいてN個の高周波送信コイルのセットに対するコイル感度を決定するためにBマッピングシーケンスをMRシステムに実行させるように構成される図1のBシミングモジュールの1つ又は複数のプロセッサによって適切に実行されるプロセスをフローチャートに描く。 一部の例示的なコイル感度マップを示す。 一部の例示的なコイル感度マップを示す。 一部の例示的なコイル感度マップを示す。 一部の例示的なコイル感度マップを示す。
図1を参照すると、画像システムは、図示されたAchieva(登録商標)磁気共鳴スキャナ(Koninklijke Philips Electronics N.V.,Eindhoven,The Netherlandsから利用可能)、又はIntera(登録商標)若しくはPanorama(登録商標)磁気共鳴スキャナ(両方ともKoninklijke Philips Electronics N.V.,Eindhoven,The Netherlandsから利用可能)、又は別の市販の磁気共鳴スキャナ、又は非商用磁気共鳴スキャナなどといった磁気共鳴(MR)スキャナ10を含む。典型的な実施形態において、磁気共鳴スキャナは、静(B)磁場を生じる超伝導又は常伝導主磁石、及び静磁場に選択された傾斜磁場を重ねるための傾斜磁場コイル巻線のセットなどの内部部品(図示せず)を含む。
MRスキャナはまた、全部でN個の高周波送信コイルのセット11も含む、又はそれらとともに機能する。例えば、高周波送信コイルのセット11は8個の送信コイルのアレイ(この場合N=8)であってもよく、又は16個の送信コイルのアレイ(この場合N=16)、20個の送信コイルのアレイ(この場合N=20)などであってもよい。高周波送信コイルのセット11は、中で共通の筐体又は支持体がN個の高周波送信コイルを含む若しくは支持する一体コイルアレイユニットとして具体化されることができ、又はN個の高周波送信コイルは個々のコイルユニットとして別々に具体化されてもよく、別個の及びグループ化された筐体若しくは支持体のいくつかの組み合わせを持ってもよい。一部の実施形態において、高周波送信コイルのセット11は複数のコイル素子を含む頭部コイルユニット、全身コイルユニット、又は他の局部コイルアセンブリとして具体化されてもよい。MRスキャナ10はさらに1つ以上の高周波受信コイルを含む、又はそれらとともに機能し、これらは高周波送信コイルのセット11によって具体化され得る(この場合高周波送信コイルのセット11のコイルは送信/受信若しくはT/Rコイルとして構成される)か、又は1つ以上の別の受信コイルとして具体化され得る(図1に図示せず)。
磁気共鳴スキャナ10及び高周波送信コイルのセット11は、磁気共鳴励起を規定する磁気共鳴シーケンスを実行し、磁気共鳴励起によって生じる磁気共鳴信号を受信するように磁気共鳴制御モジュール12によって制御される。MR画像応用例において、高周波送信コイルのセット11は磁気共鳴を励起するために励起され、これはMRスキャナ10の傾斜磁場コイルによって選択的に印加される傾斜磁場によって空間エンコードされ、励起され空間エンコードされた磁気共鳴信号は1つ又は複数の受信コイルを用いて読み取られる。随意に、磁気共鳴励起は磁気共鳴励起位相中にMRスキャナ10の傾斜磁場コイルによって印加される傾斜磁場によってスライス又は他の励起領域に空間的に制限される。
この画像応用例の場合、再構成モジュール14は収集された磁気共鳴信号を再構成して磁気共鳴画像を生成し、これは磁気共鳴画像メモリ16に保存される。他の応用の場合、適切な収集後処理が画像再構成処理の代わりに又は組み合わせて適用される。例えば、MR分光法応用において、選択位置からのMR信号は周波数の関数としてプロットされ得るか、又はMR信号は周波数などに基づいてフィルタされ得る。一部の実施形態において、部品12、14、16は磁気共鳴スキャナ10の製造業者及び/又は1つ以上の第三者ベンダによって供給される汎用市販磁気共鳴画像機器である。代替的に、部品12、14、16のうちの1つ以上又は全ては特注若しくは顧客が変更した部品であってもよい。
引き続き図1を参照すると、Bシミングモジュール20はMRシステムと、例えばMR制御モジュール12と動作可能なように通信し、高周波送信コイルのセット11によって生じるB送信磁場の空間的均一性を高めるためにMRシステムにBシミングを実行させる。このような目的で、プロセッサ22はBマッピングを実行するためのBマッピングシーケンスの減少した回数を決定するように構成される。"減少した回数"とは、実行されるべきBマッピングシーケンスの回数が高周波送信コイルのセット11の高周波送信コイルの数Nよりも少ないことを意味する。プロセッサ22はB磁場マッピングシーケンス24を出力し、その回数はN未満である。プロセッサ26はB磁場マッピングデータセットを作り出すようにMRシステム10、11、12にB磁場マッピングシーケンス24を実行させるように構成される。プロセッサ28は収集されたB磁場マッピングデータセットに基づいてN個の高周波送信コイルのセット11に対するコイル感度30を決定するように構成される。
図1は別々に図示されたプロセッサ22、26、28を有するシミングモジュール20を図示するが、一部の実施形態においてプロセッサ22、26、28は随意にMR制御モジュール12を具体化する同じプロセッサである単一のプロセッサとして具体化され得る。例えば、図示されたコンピュータCは様々な部品12、22、26、28を具体化するために記憶媒体(ハードディスク若しくは他の磁気ディスク、光ディスク、フラッシュメモリ若しくは他の静電メモリ、ランダムアクセスメモリ、リードオンリーメモリなど)によってプログラムされるデジタルプロセッサ(図示せず)を、また随意に再構成モジュール14も適切に含む。コンピュータ例Cはまた、ディスプレイD及び図示されたキーボードKなどの1つ以上のユーザ入力装置も含む。例えば、ディスプレイDは画像メモリ16から読み出される再構成画像を表示するために使用されることができ、一方1つ以上のユーザ入力装置KはMRシステム10、11、12、再構成モジュール14、及びBシミングモジュール20を操作するためのコマンドを放射線科医若しくは他のユーザが入力することを可能にするために使用されることができる。
生成されたコイル感度30は、少なくとも視野(FOV)内で高周波送信コイルのセット11によって生じるB1送信磁場の空間的均一性を高めるために、MRシステム10、11、12によって実行される磁気共鳴シーケンスの送信位相中に高周波送信コイルのセット11のBシミングを実行するために適切に使用される。
シミングモジュール20の効率は、B磁場マッピングシーケンス24の回数を高周波送信コイルのセット11の送信コイルの数Nよりも十分に少なくすることによって促進される。一部の実施形態において、例えば、B磁場マッピングシーケンス24の回数はNの半分であり、一方他の実施形態ではB磁場マッピングシーケンス24の回数はNの半分より少ない。プロセッサ22は、減少した回数のB磁場マッピングシーケンス24が、1回のB磁場マッピングシーケンスが各送信コイルに対して実行される(Nシーケンス)従来のB磁場マッピングシーケンスのセットによって提供されるBマッピング情報の大部分を依然として提供することができるように、B磁場マッピングシーケンス24を選択するように構成される。
引き続き図1を参照し、さらに図2を参照すると、プロセッサ22によって適切に実行されるプロセス例が記載される。直交仮想高周波送信コイルのセット42を作るために、高周波送信コイルのセット11に対して線形変換40が実行される。線形変換40は低周波B送信磁場成分から高周波B送信磁場成分を分けるように適切に選ばれる。例えば、一部の実施形態において、線形変換40は適切には固有モード変換であり、その場合高次固有モードは典型的には高周波B送信磁場成分を具体化し、一方低次固有モードは典型的には低周波B送信磁場成分を具体化する。選択操作44において、直交仮想高周波送信コイルのセット42のサブセット46が選択される。サブセット46は、低周波B送信磁場成分を十分に保ちながら、主に高い空間周波数成分を持つB送信磁場を生じる仮想高周波送信コイルを選択されたサブセット46から除外するように適切に選択される。例えば、線形変換40が固有モード変換である実施形態において、選択操作44は、最低固有モード次数を持つ直交仮想高周波送信コイルのセット42の複数の仮想高周波送信コイルを適切に選択する。
例えば、高周波コイルのセット11がN=8コイルを含む実施例を考慮すると、直交仮想高周波送信コイルのセット42が固有モード0‐7に対応する仮想高周波コイルを含むように固有モード変換が使用される。B磁場マッピングシーケンスの回数がコイルの数の半分であることが望ましい場合(つまりこの実施例では4回のB磁場マッピングシーケンス)、選択操作44は最低固有モ−ド次数0‐3に対応する4個の仮想高周波コイルから構成されるようにサブセット46を適切に選択する。
いかなる特定の操作理論にも限定されることなく、Bシミングは空間的に平滑な操作として最適に実施されることが期待される。Bシミングは低い空間周波数に対応するB送信磁場における漸進的な又は大規模の空間的変動の補正に有効であることが期待される。他方で、Bシミングは高い空間周波数に対応するB送信磁場における突然の又は小規模の空間的変動の補正にはあまり有効でないことが予想される。従って、選択されたサブセット46において、主に低い空間周波数成分を持つB送信磁場を生じる仮想高周波送信コイルを保ちながら、主に高い空間周波数成分を持つB送信磁場を生じる仮想高周波送信コイルを選択されたサブセット46から除外することによって、結果として得られるサブセット46は、Bシミングが補正に有効である可能性が低い高い空間周波数情報を捨てながら、Bシミングによって容易に補正可能な低い空間周波数情報を提供することが期待される。
直交仮想高周波送信コイルのセット42のサブセット46は、好適にはB磁場マッピングシーケンス24によって収集される最も適切なBマッピング情報を特定する。しかしながら、サブセット46は線形変換40に基づく表現である。他方で、B磁場マッピングシーケンス24はN個の高周波送信コイルのセット11を用いて実施される。従って、サブセット46の各仮想高周波送信コイルは選択操作50において順に選択され、選択された仮想コイルに対してエンコード操作52が実行される。エンコード操作52はエンコード基準としてN個の高周波送信コイルのセットを用いて選択された仮想高周波送信コイルをエンコードする。エンコード操作52は、選択された振幅と位相において動作する動作中のN個の高周波送信コイル11の重ね合わせが、選択された仮想高周波送信コイルを再現するように、N個の高周波送信コイル11に対する振幅と位相を生成する(54)。このプロセス50、52、54はB磁場マッピングシーケンス24を作るためにサブセット46の各仮想コイルに対して繰り返される。要するに、直交仮想高周波送信コイルのセット42はエンコード操作50、52、54を介してB磁場マッピングシーケンス24を規定する。
引き続き図1を参照し、さらに図3を参照すると、プロセッサ26、28によって適切に実行されるプロセス例が記載される。プロセッサ26は、プロセッサ22によって決定されるB磁場マッピングシーケンス24をMRシステム10、11、12に実行させる操作60を実行する。収集操作60は収集されたB磁場マッピングデータセット62を生成し、その成分は仮想高周波送信コイルのサブセット46に対応するB送信磁場の重ね合わせをマップする。物理的な高周波送信コイルのセット11の高周波送信コイルに対するBマッピングを提供するために、収集されたB磁場マッピングデータセット62は、例えば線形変換40(図2参照)の逆である図示された逆線形変換64によって変換され、N個の高周波送信コイルのセット11に基づいて対応する収集されたB磁場マッピングデータセット66を生成する。B磁場マッピングデータセット66は、N個の物理的コイルを連続的に用いて実行されるN回のB磁場マッピングシーケンスを用いて従来収集されるようなB磁場マッピングデータセットに対応する。従って、計算操作68はB磁場マッピングデータセット66からN個の高周波送信コイルのセット11のN個のコイルに対するコイル感度30を計算する。一部の実施形態において、計算操作68はN個の高周波送信コイルのセット11のN個のコイルに対するコイル感度30を計算するために第一原理的に(ab initio)演算する。他の実施形態においてはコイル感度30についての先験的情報も利用される。例えばコイル感度30についての先験的情報は高周波送信コイルのセット11に対する保存されたコイル感度のセットの形をとってもよく、計算操作68はN個の高周波送信コイルのセット11に基づいてあらわされる収集されたB磁場マッピングデータセット66に基づいて保存されたコイル感度のセットを調節する。
開示されたB送信マッピングの一実施例が記載される。高周波送信コイルのセット11の各個々のコイルのコイル感度分布(Bマップ)は、行列Sによって記載される物理的コイル11を、仮想コイルが固有モードに対応する対応仮想固有モード系にマップする固有モード変換を用いて固有‐コイル系に変換される。行列Sの行は個々の送信コイルの空間感度を含むベクトルを成す。コイル感度SはS=U×Σ×Vとしてその特異値分解(SVD)表現で与えられることができ、行列Σは特異値を含み、これはコイルアレイの固有値に対応する。行列U、又はより正確にはUのエルミート(U)は、行列Sによってあらわされる物理的コイルを、E=U×Sによって与えられる固有モードEに対応する仮想コイルにマップする射影行列となる。
個々のコイルのBマップを測定する代わりに、固有モードを形成するそれらの適切な重ね合わせを測定することもできる。固有モードは行列Σにおいて与えられる特異値によって、全コイル系の特性に対するそれらの寄与に従ってランクされることができる。このランキングはBマッピングにとって重要な情報提供のための内容を失うことなく、定量的な方法で最も重要でない固有モードコイルを除外する機会を与える。むしろ、Bマッピングに最も強く寄与するごく少数の固有モードが測定される。
マッピング情報にほとんど寄与しない固有モードを無視することは、E'=R×U×Sに従って系に簡約演算子Rを適用することによって数学的に実現されることができ、E'は簡約化された固有‐コイル系を示す。演算子Rは単位行列と同様であるが、Bマッピング測定中に除外されるべき固有モードの位置でゼロを持つ。当然のことながら測定されるべき固有モードの数を制限することは一般的方法であり、測定結果を得るために使用されるBマッピングシーケンスから独立して利用されることができる。Bマッピングは二次元又は三次元のいずれかで実行されることができる。高周波送信コイルのセット11のコイル感度についての粗い先験的情報に基づき、Bマッピングにとってほとんどの情報に寄与する固有モードを得るための線形変換が推定されることができ、従ってコイルアレイは少ない数の仮想コイル素子から成るコイルとみなされることができる。
選択された固有モードの測定は随意に直接実行されない。この方法はBマッピングシーケンスがいくらかの不正確さを持つ可能性があることを認識する。例えば、Bマップは低感度の空間領域を含む高いダイナミックレンジを持ち得る。高次固有モードは低い信号エネルギーを含み、こうした低感度領域をより多く示し、それらを十分な精度で測定することを困難にする。従って一部の実施形態において、Bマッピングプロセスの精度を向上させるために、固有モードは、例えば物理的コイル11の1つを除く全てが使用される"1つを除く全ての"方法を用いてモードの重ね合わせによってエンコードされる。行列Aが、異なるBマッピングステップ中に適用される"1つを除く全ての"エンコード法に従って修正されたコイルモードZを作るための個々の固有モードの重ね合わせをあらわすとすると、次式はZ=A×E'=A×R×U×Sと書き換えられることができる。積(A×R×U)はBマッピングシーケンスに対する高周波送信コイルのセット11の個々の送信コイル素子に対する振幅と位相を与える。一番右の式におけるコイル感度情報Sは送信コイルアレイ11についての先験的知識から適切に導出される。
測定後、式Z=A×E'=A×R×U×Sから得られる適切な変換は、物理的な高周波送信コイルのセット11に対するコイル感度30の対応する表現を導出することを可能にし、そしてこれはBシミングのために、又は並列送信高周波パルスパラメータを決定するために適切に使用される。固有モードのサブセットに基づいて高周波パルスを設計することもまた便利である。固有モード空間において演算することによって、総計算時間は物理的コイル11に基づいて演算する計算と比較して削減されることができる。
二次元又は三次元で実行される前のコイル感度較正からの知識もまた、適切な学習手順を用いて利用されることができる。例えば、コイルのセット11の各物理的コイルに対する三次元B感度行列は行列S'に保存されることができる。この行列は新たなBマッピング測定に基づいて個々のコイル振幅及び位相を導出するために調節されるコイル感度の平均モデルとなる。新たなBマッピング測定の各々は二次元又は三次元で実行されるので、保存されたモデルS'は、例えば適切な重み付け平均化プロセスを用いて、新たに得られる情報を組み込むことによって適切に調節される。このプロセスは随意に異なる空間位置における既存のモデルS'の信頼性を考慮に入れる。モデルS'は、例えば物理座標系において(つまり高周波送信コイルのセット11に基づいて)、又は固有‐コイル系に基づいて、異なる形式で保存されることができる。後者の保存法は、無視される固有値の数に依存して、保存されたデータ精度の異なるレベルを許容する。この精度のレベルをユーザが選択することが考慮される。感度S'を保存し、保存された感度を各Bマッピング反復で調節することは、高周波コイルのセット11のコイル間の相対的空間的配置が共通の筐体又は支持体によって予め規定される、全身コイルなどの固定され取り付けられた送信アレイの場合に特に有用である。
保存されたS'において与えられる保存されたコイル感度情報の情報を組み込むことによって、関係式Z=A×E'=A×R×U×SからBマッピングシーケンスに対する比吸収率(SAR)が推定されることができることも有利である。
図4‐7を参照すると、上記Bマッピングに対するいくつかのBマップが提示される。ここで対象とする用途は斜め方向に二次元スキャンとして実行される全身MR画像である。こうした撮像におけるB変動を補正するために、Bシミングが有利に実行される。ここでは8チャネル送信システムが想定される(つまりこの実施例の場合N=8)。スライス幾何学情報(FOV、スライス方向)に基づいて、3D trainedコイル感度は内挿及び外挿法を用いて表現S'から再フォーマットされる。個々のコイルに対する再フォーマットされた感度は関係式E=U×Sに従って分析される。この実施例において、4個の固有モードがBマッピングのために使用され、これは8コイルシステムに対するBマッピング時間を1/2に減らす。"1つを除く全ての"測定法が選ばれる。対応するエンコード行列Aは簡単にするために実数1であるとみなされる。Aの全成分は1に等しいが、"1つを除く全ての"測定法を実施するために対角成分はゼロに設定された。(これは単に一実施例である。他のエンコード行列もまた考慮される。例えば他の実施形態においては、測定されるべき全モードに対するマッピング手順の精度を向上させるように最適エンコード行列Aを選択するために、異なる固有モードに対する対応する固有値が考慮される。)エンコード法を選択後、高速Bマッピングのためのチャネル駆動スケールが導出され、所望の空間分解能で測定が実行される。4個の測定された固有モードに基づいて、コイル感度、従ってBシミング係数が導出される。図4は8コイルシステムに対するコイル感度(Bマップ)を示す。図5は次数0‐7の固有モードを示す。図6は"1つを除く全ての"方法を用いてエンコードされる4個の最低次数固有モードのみの測定からあらわされるように計算される8コイルシステムに対するコイル感度(Bマップ)を示す。図4及び6の比較は、4個の最高次数固有モードの除外が、得られるコイル感度に対して比較的ほとんど影響を及ぼさないことを示す。最後に、図7の一番上の行は図6の一番上の行を再現し、一方図7の一番下の行は(i)高次の固有モードを除外することによって得られるこれらのコイル感度と、(ii)図4に示される実際のコイル感度との間の差を示す。再度、FOVの外側端付近を除いてほとんど差は見られない。
本発明は好適な実施形態に関して記載されている。修正と変更は上記の詳細な説明を読んで理解することで想到され得る。本発明はこうした修正と変更を添付のクレーム又はその均等物の範囲内にある限り全て含むものと解釈されることが意図される。クレームにおいて括弧の間に置かれる任意の参照符号はクレームを限定するものと解釈されてはならない。"有する"という語はクレームに列挙されたもの以外の要素又はステップの存在を除外しない。ある要素に先行する"a"又は"an"という語はこうした要素の複数の存在を除外しない。開示された方法は、複数の別個の要素を有するハードウェアを用いて、及び適切にプログラムされたコンピュータを用いて実施されることができる。複数の手段を列挙するシステムクレームにおいて、これらの手段のいくつかはコンピュータ可読ソフトウェア又はハードウェアの1つの同じ項目によって具体化されることができる。特定の手段が相互に異なる従属クレームにおいて列挙されているという単なる事実は、これらの手段の組み合わせが有利に使用されることができないことを示すものではない。

Claims (15)

  1. 磁場マッピングデータセットを収集するために高周波送信コイルのセットを用いてある回数のB磁場マッピングシーケンスを実行するステップであって、前記回数は前記高周波送信コイルのセットにおける高周波送信コイルの数よりも少ない、ステップと、
    前記収集されたB磁場マッピングデータセットに基づいて前記高周波送信コイルのセットに対するコイル感度を決定するステップとを有する、方法。
  2. 前記実行されるB磁場マッピングシーケンスが、(i)直交仮想高周波送信コイルのセットを作るために前記高周波送信コイルのセットに線形変換を実行するステップと、(ii)前記実行されるB磁場マッピングシーケンスを規定する前記直交仮想高周波送信コイルのセットのサブセットを選択するステップとによって規定される、請求項1に記載の方法。
  3. 前記線形変換が固有モード変換である、請求項2に記載の方法。
  4. 前記選択操作(ii)が、
    前記直交仮想高周波送信コイルのセットの最低固有モード次数を持つ複数の仮想高周波送信コイルを選択するステップを有する、請求項3に記載の方法。
  5. 前記実行されるB磁場マッピングシーケンスがさらに、(iii)前記高周波送信コイルのセットの前記高周波送信コイルに対する対応する振幅及び位相の値を決定するためにエンコード基準として前記高周波送信コイルを用いて前記選択されたサブセットのコイルをエンコードするステップによって規定される、請求項2に記載の方法。
  6. 前記選択操作(ii)が、主に高い空間周波数成分を持つB送信磁場を生じる仮想高周波送信コイルを、前記選択されたサブセットから除外するステップを有する、請求項2に記載の方法。
  7. 前記コイル感度を決定するステップが、前記高周波送信コイルのセットに基づいて収集されたB磁場マッピングデータセットを生成するために、前記収集されたB磁場マッピングデータセットに前記線形変換に対応する逆変換を適用するステップを有する、請求項2に記載の方法。
  8. 前記実行されるB磁場マッピングシーケンスの回数が、前記高周波送信コイルのセットにおける前記高周波送信コイルの数の半分以下である、請求項1に記載の方法。
  9. 各B磁場マッピングシーケンスが、前記高周波送信コイルのセットにおける前記高周波送信コイルの1つを除く全てを用いて実行される、請求項1に記載の方法。
  10. 実行されるB磁場マッピングシーケンスの各々が前記高周波送信コイルのセットにおける前記高周波送信コイルの異なる1つを除外する、請求項1に記載の方法。
  11. 前記収集されたB磁場マッピングデータセットに基づいて前記高周波送信コイルのセットに対するコイル感度を決定するステップが、
    前記収集されたB磁場マッピングデータセットに基づいて前記高周波送信コイルのセットに対する保存されたコイル感度のセットを調節するステップを有する、請求項1に記載の方法。
  12. 前記収集されたB磁場マッピングデータセットに基づいて前記高周波送信コイルのセットに対して決定された前記コイル感度を用いてBシミングを実行するステップをさらに有する、請求項1に記載の方法。
  13. 請求項1に記載の方法を実行するためにデジタルプロセッサによって実行可能な命令を保存する記憶媒体。
  14. 請求項1に記載の方法を実行するデジタルプロセッサ。
  15. 磁気共鳴システムであって、
    磁気共鳴スキャナと、
    高周波送信コイルのセットと、
    請求項1に記載の方法を実行するデジタルプロセッサとを有し、ある回数のB磁場マッピングシーケンスを実行するステップは、前記磁気共鳴スキャナに前記高周波送信コイルのセットを用いて前記B磁場マッピングシーケンスを実行させるステップを有する、磁気共鳴システム。
JP2012502841A 2009-03-31 2010-03-22 加速b1マッピング Active JP5525596B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP09156820 2009-03-31
EP09156820.4 2009-03-31
PCT/IB2010/051225 WO2010113062A1 (en) 2009-03-31 2010-03-22 Accelerated b1 mapping

Publications (2)

Publication Number Publication Date
JP2012521843A true JP2012521843A (ja) 2012-09-20
JP5525596B2 JP5525596B2 (ja) 2014-06-18

Family

ID=42232637

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012502841A Active JP5525596B2 (ja) 2009-03-31 2010-03-22 加速b1マッピング

Country Status (6)

Country Link
US (1) US8847593B2 (ja)
EP (1) EP2417470A1 (ja)
JP (1) JP5525596B2 (ja)
CN (1) CN102369451B (ja)
RU (1) RU2534724C2 (ja)
WO (1) WO2010113062A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014525812A (ja) * 2011-08-17 2014-10-02 コーニンクレッカ フィリップス エヌ ヴェ 磁気共鳴撮像中における既定の体積内の無線周波数伝送場の減少
JP2015054252A (ja) * 2013-09-11 2015-03-23 シーメンス アクチエンゲゼルシヤフトSiemens Aktiengesellschaft B1マップの算出

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2497692B (en) * 2010-08-20 2017-06-07 Koninklijke Philips Nv Virtual coil emulation in parallel transmission mri
US9157979B2 (en) 2011-05-02 2015-10-13 Toshiba Medical Systems Corporation Efficient multi-station MRI
JP5677226B2 (ja) * 2011-07-28 2015-02-25 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴装置およびプログラム
CN103874457B (zh) * 2011-09-29 2016-02-03 株式会社日立医疗器械 磁共振摄像装置以及高频磁场照射方法
WO2013069513A1 (ja) * 2011-11-08 2013-05-16 株式会社 日立メディコ 磁気共鳴イメージング装置および照射磁場分布計測方法
US20130289912A1 (en) * 2012-03-30 2013-10-31 Siemens Aktiengesellschaft Eigen-vector approach for coil sensitivity maps estimation
TW201351861A (zh) * 2012-06-08 2013-12-16 Novatek Microelectronics Corp 控制電源轉換裝置的方法及其相關電路
DE102013205785B4 (de) * 2013-04-02 2014-10-09 Siemens Aktiengesellschaft Ermittlung einer Magnetresonanzsystem-Ansteuersequenz auf Basis einer reduzierten Anzahl an Feldverteilungskarten
WO2016183572A1 (en) 2015-05-14 2016-11-17 Ohio State Innovation Foundation Systems and methods for estimating complex b1+ fields of transmit coils of a magnetic resonance imaging (mri) system
WO2018136705A1 (en) 2017-01-19 2018-07-26 Ohio State Innovation Foundation Estimating absolute phase of radio frequency fields of transmit and receive coils in a magnetic resonance
EP3511727A1 (en) * 2018-01-11 2019-07-17 Koninklijke Philips N.V. Active b1+ shimming of transmission coils
US10527699B1 (en) * 2018-08-01 2020-01-07 The Board Of Trustees Of The Leland Stanford Junior University Unsupervised deep learning for multi-channel MRI model estimation

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0722576B2 (ja) * 1988-04-11 1995-03-15 三菱電機株式会社 磁気共鳴装置
US6294972B1 (en) * 2000-08-03 2001-09-25 The Mcw Research Foundation, Inc. Method for shimming a static magnetic field in a local MRI coil
DE102004045691B4 (de) 2003-10-27 2009-10-01 Siemens Ag Verfahren zum Erzeugen eines homogenen hochfrequenten Magnetfelds in einem räumlichen Untersuchungsvolumen einer Magnetresonanzanlage
DE102004013422B4 (de) * 2004-03-18 2009-02-19 Siemens Ag Verfahren zur Homogenisierung eines B1-Felds, Magnetresonanzsystem und Computerprogrammprodukt
US20070170917A1 (en) * 2005-01-05 2007-07-26 Koninklijke Philips Electronics N.V. Magnetic resonance imaging with real-time magnetic filed mapping
DE102006017439B3 (de) * 2006-04-13 2007-10-11 Siemens Ag Verfahren und System zur Steuerung einer Magnetresonanzanlage
WO2007124245A1 (en) * 2006-04-21 2007-11-01 Koninklijke Philips Electronics, N.V. Magnetic resonance with time sequential spin excitation
US20090080750A1 (en) * 2006-04-21 2009-03-26 Koninklijke Philips Electronics N. V. Passive mr visualisation of interventional instruments

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014525812A (ja) * 2011-08-17 2014-10-02 コーニンクレッカ フィリップス エヌ ヴェ 磁気共鳴撮像中における既定の体積内の無線周波数伝送場の減少
JP2015054252A (ja) * 2013-09-11 2015-03-23 シーメンス アクチエンゲゼルシヤフトSiemens Aktiengesellschaft B1マップの算出

Also Published As

Publication number Publication date
JP5525596B2 (ja) 2014-06-18
CN102369451B (zh) 2015-02-25
WO2010113062A1 (en) 2010-10-07
CN102369451A (zh) 2012-03-07
EP2417470A1 (en) 2012-02-15
US20120007600A1 (en) 2012-01-12
RU2534724C2 (ru) 2014-12-10
US8847593B2 (en) 2014-09-30
RU2011143726A (ru) 2013-05-10

Similar Documents

Publication Publication Date Title
JP5525596B2 (ja) 加速b1マッピング
JP5731373B2 (ja) 患者の安全性、及び走査性能の改良のためのリアルタイム性の局所的及び大域的なsar推定
JP5706899B2 (ja) 磁気共鳴イメージングのための画像強度訂正
US9317917B2 (en) Method, reconstruction device, and magnetic resonance apparatus for reconstructing magnetic resonance raw data
JP5619136B2 (ja) k空間に依存するRFパルス選択によるパラレル送信におけるSAR低減
JP2013505046A (ja) Mriにおけるrfパワー及びrfフィールド均一性の同時最適化
US10203394B2 (en) Metal resistant MR imaging
US8686726B2 (en) RF shimming with RF power regularization using a multi-channel RF transmit system for MRI
US9229081B2 (en) Accelerated MRI with nonlinear spatial encoding gradients
CN107072586B (zh) 磁共振成像装置
US10185008B2 (en) Medical imaging examination apparatus and operating method
JP2017519576A (ja) マルチショット磁気共鳴(mr)イメージングシステム及びその操作方法
US10001538B2 (en) Method and magnetic resonance system that determines operational sample points in k-space of a random undersampling scheme when acquiring magnetic resonance data
US20130281822A1 (en) Image correction using multichannel blind deconvolution with homomorphic filtering
JP4975614B2 (ja) 磁気共鳴イメージング装置及び方法
KR101751974B1 (ko) 자기 공명 영상 장치 및 이를 이용한 자기 공명 영상 처리 방법
US11061097B2 (en) Data driven correction of phase depending artefacts in a magnetic resonance imaging system
JP7183048B2 (ja) 磁気共鳴イメージングシステム、磁気共鳴イメージング方法、及び磁気共鳴イメージングプログラム
US20160124065A1 (en) Method and apparatus for correction of magnetic resonance image recordings with the use of a converted field map
Giliyar Radhakrishna et al. Improving spreading projection algorithm for rapid k‐space sampling trajectories through minimized off‐resonance effects and gridding of low frequencies
US8396277B2 (en) Magnetic resonance imaging apparatus
US20240183921A1 (en) Magnetic resonance system, magnetic resonance image correction method, and magnetic resonance imaging method
KR101942854B1 (ko) 자기 공명 영상 장치 및 이를 이용한 자기 공명 영상 처리 방법
JPH0686763A (ja) 磁気共鳴診断装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20130318

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20140226

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140411

R150 Certificate of patent or registration of utility model

Ref document number: 5525596

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250