Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP5619136B2 - k空間に依存するRFパルス選択によるパラレル送信におけるSAR低減 - Google Patents
[go: Go Back, main page]

JP5619136B2 - k空間に依存するRFパルス選択によるパラレル送信におけるSAR低減 - Google Patents

k空間に依存するRFパルス選択によるパラレル送信におけるSAR低減 Download PDF

Info

Publication number
JP5619136B2
JP5619136B2 JP2012502867A JP2012502867A JP5619136B2 JP 5619136 B2 JP5619136 B2 JP 5619136B2 JP 2012502867 A JP2012502867 A JP 2012502867A JP 2012502867 A JP2012502867 A JP 2012502867A JP 5619136 B2 JP5619136 B2 JP 5619136B2
Authority
JP
Japan
Prior art keywords
pulse
space
sar
data
scan
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
JP2012502867A
Other languages
English (en)
Other versions
JP2012522563A (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 JP2012522563A publication Critical patent/JP2012522563A/ja
Application granted granted Critical
Publication of JP5619136B2 publication Critical patent/JP5619136B2/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/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/32Excitation or detection systems, e.g. using radio frequency signals
    • G01R33/34Constructional details, e.g. resonators, specially adapted to MR
    • G01R33/341Constructional details, e.g. resonators, specially adapted to MR comprising surface coils
    • G01R33/3415Constructional details, e.g. resonators, specially adapted to MR comprising surface coils comprising arrays of sub-coils, i.e. phased-array coils with flexible receiver 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/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/288Provisions within MR facilities for enhancing safety during MR, e.g. reduction of the specific absorption rate [SAR], detection of ferromagnetic objects in the scanner room
    • 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/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

Landscapes

  • Physics & Mathematics (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)
  • High Energy & Nuclear Physics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Description

本発明は、特にMRI(Magnetic Resonance Imaging)システムに適用される。しかしながら、開示された技術はまた、他のイメージングシステム、他の磁気共鳴シナリオ、他の画像データ収集技術などに適用されてもよい。
MRIシステムにおけるパラレル送信は、最近数年間において大きな関心のあるトピックになってきている。複数の個別の無線周波数(RF)送信コイルの利用は、B均一性の制限を解消し、期間を短縮することによって多次元RFパルスを改善するのに利用される。これらすべての送信用途において、過剰な患者の加熱を回避するためにあるリミット以下に維持される必要がある比吸収率(Specific Absorption Rate:SAR)が特に関心がある。パラレル送信においてSARを低減するための各種アプローチが議論されてきた。RFパルス設計における自由度は、例えば、低SARを実現するのに役立つ正則化(regularization)技術などを介し、最小のSARによる手段の選択を可能にする。さらに、k空間軌跡とRF波形との間の相互作用が、SARの低減に利用可能である(例えば、可変レート選択的励起やVERSE技術など)。このように取得される最適なRFパルスは、対応するMRスキャンのために利用される。これに関して、パラレル送信RFパルスは、MR信号サンプリング処理からほぼ独立して最適化される。SARの制約は、特に高磁場ではMRにおいて重大な問題である。
本出願は、上述した問題点などを解消するMR検査においてSARを低減するための新規かつ改良されたシステム及び方法を提供する。
一態様によると、低減された比吸収率(SAR)を有するMRデータを取得する方法は、2以上の送信コイルを用いてMR取得スキャン中にMRスキャンデータを生成するため、少なくとも1つのRFパルスが繰り返し印加される磁気共鳴シーケンスを適用するステップと、前記RFパルスがいくつかの繰り返しにおいて他の繰り返しより低いSARを生じさせるように、前記繰り返し印加されるRFパルスの構成を変更するステップと、k空間において前記MRスキャンデータを取得するステップとを有する。
他の態様によると、比吸収率(SAR)を低減する磁気共鳴システムは、2以上の送信RFコイルと、RFパルス生成手段と、少なくとも1つの繰り返し印加されるRFパルスによりMRスキャンシーケンスを適用するよう前記RFパルス生成手段を制御するプロセッサとを有する。前記プロセッサは、各バージョンが異なるSARを示し、異なる周波数コンポーネントを有する前記繰り返し印加されるRFパルスの少なくとも2つの予め生成されたバージョンから選択する。さらに、前記プロセッサは、k空間の第1領域からMRデータを取得する際には前記繰り返し印加されるRFパルスのより高い周波数のバージョンを、k空間の第2領域からMRデータを取得する際には前記繰り返し印加されるRFパルスのより低い周波数のバージョンを前記RFパルス生成手段に提供し、前記MRデータを処理する。本システムはさらに、前記処理されたMRデータを表示するディスプレイを有する。
1つの効果は、スキャン全体の合計のSARが低減されることである。
他の効果は、画質が向上又は少なくとも維持されることである。
さらなる効果は、以下の詳細な説明を読み、理解した当業者によって理解されるであろう。
図1は、k空間位置の関数として取得スキャン中にRFパルスを変更することによってマルチトランスミットMRI装置におけるSARを低減するシステムを示す。 図2は、RFパルスのSARとその全体的な精度及び/又はパフォーマンスとの間の関係のグラフを示す。 図3は、4チャネルパラレルトランスミットシステム上で2Dの空間選択的RFパルスを利用してローカルな励起のための励起リードアウトを示す。 図4は、空間精細さにおいて異なるRFパルスのクラスを用いて実行されるk空間の各ラインの信号励起の図を示す。 図5は、2Dサンプリングk空間における異なるエリアの循環的なセグメントにより形成される構成を示す。対応する位相符号化ステップk,kが測定される場合、8つの異なるパラレルトランスミットRFパルスが予め計算され、適用される。 図6は、8チャネル送信コイルにおける8つの横軸Bコイルのコイル感度を示す。 図7は、8チャネルコイル構成の対応する固有モードを示す。 図8は、ここに開示される1以上の態様によるマルチチャネルトランスミットコイルMRI装置においてSARを低減する方法を示す。
図1は、k空間位置の関数として取得スキャン中にRFパルスを変更することによってマルチトランスミット(multi−transmit)MRI装置におけるSARを低減するシステム10を示す。例えば、RFパルスは、予め選択された空間精細さ(spatial definition)により設計される。より高い空間精細さによるパルスは、より正確な画像を与えるが、より高いSARを有する。k空間のいくつかの領域におけるシーケンスに対して標準的な空間精細さを有するRFパルスと、他の領域において低減された空間精細さを有するRFパルスとを使用することによって、シーケンスの全体的なSARは低減される。例えば、より高い空間精細さのパルスは、k空間の低周波数部分に用いられ、より低い空間精細さのパルスが、k空間の高周波数領域に用いられるが、他の割当てもまた可能である。特に複数の独立に動作するRF送信機を有するMRスキャナでは、Bシミングパルス(shimming pulse)を設計するためのより大きなフレキシビリティがある。このフレキシビリティは、シーケンスの全体的なSARを低減するため、k空間領域に従ってBシミングを調整するために、開示されるシステム及び方法によって利用される。
システム10は、プロセッサ14、メモリ16(コンピュータ可読媒体など)及び情報をユーザに表示するディスプレイ18に結合されるMR装置又はスキャナ12を有する。メモリは、ここに開示される各種機能を実行するための1以上のコンピュータにより実行可能な命令を格納し、プロセッサがそれを実行する。メモリ、プロセッサ及び/又はディスプレイは、MR装置12から分離されてもよいし、又は一体化されていてもよいことが理解されるであろう。
MR装置12は、各パルスがプロセッサにより指示されるように構成される無線周波数(RF)Bパルスを生成するRFパルス生成手段20を有する。MR装置はまた、複数のチャネルを介しRFパルスを送信するマルチチャネル送信機(MCTx)22を有する。メモリ16は、k空間の特定位置に相互参照される複数の予め生成されたRFパルスパターンを有するRFパルスルックアップテーブル(LUT)24を格納する。LUT24は、k空間の各領域がスキャンされるとき、特定のシーケンス又はパターンにおいてRFPG20により生成され、MCTx22により送信される特定のRFパルスを特定するため、MR取得スキャン中にプロセッサによりアクセスされる。他の態様によると、RFパルスは、予め計算され、LUT24に格納され、必要に応じて呼び出される。スキャン中に検出されたk空間データ26は、メモリ16に格納される。マルチチャネルMRシステムでは、通常は複数の独立したRFコイルがあり、各コイルは、マルチチャネル送信機22の対応するRF送信機又はチャネルにより独立に駆動される。
メモリは、不定のSNRにおいて固定された最適なRFパルスを用いた取得スキャンに関して、1以上の可変的なRFパルスシーケンス又はスキームのパフォーマンスを推定するためプロセッサにより実行されるNRMSE(Normalized Root Mean Square Error)計算アルゴリズム28を格納する。所与のNRMSEに対して、プロセッサは、スキャンにおいて異なる又は可変的なRFパルスを利用する固有モード(SAR)と固定RFパルスモード(SAR)との平均SARを比較する。シミングアルゴリズム30がまたメモリ16に格納され、k空間の領域に応じてBシミングを調整するためプロセッサ14により実行される。例えば、1以上の固有モードは、k空間の高周波数(周辺)領域に対してスキャンデータを取得する際にBフィールドをシミング処理するのに利用され、単一の中心の固有値が、k空間の低周波数(中心)領域に対してデータを取得する際にBフィールドをシミング処理するのに利用される。しかしながら、他のスキームもまた想到しうることは理解されるであろう。
メモリは、MR装置12のマルチエレメントコイルの送信コイル感度データ32を格納し、当該情報はプロセッサ14により解析される。例えば、各コイルのコイル感度情報は、各固有モード34に分解される。所望の精度レベルに基づき、RFシミングに対して異なる個数の固有モードが考えられる。このように、プロセッサは、所与のk空間の領域において所望のSARを実現するため適切なシミングアルゴリズム30を選択する。
メモリはさらに、プロセッサが取得したk空間データ26をSARが低減されたMR画像38に再構成するのに実行する1以上の再構成アルゴリズム36(コンピュータにより実行可能な命令など)を格納する。
MCTx22により実現されるパラレル送信は、高磁場MRI(3Tなどより大きい)においてB送信フィールド均一性を向上させ、又は精緻化された多次元RFパルスを促進するための効率的なアプローチである。従来の高磁場アプリケーションに関する1つの問題は、限定的な比吸収率(SAR)である。システム10は、パラレル送信技術を利用するMR装置12においてSARを低減する。MRIスキャナ又はMRS(MR Spectroscopy)スキャナのサンプリングされたk空間における異なるエリアは、信号不完全性に対して異なる感度を示す。
マルチチャネル送信スキャナでは、各Bパルスは、各チャネルからのBコンポーネントの連結である。各チャネルの相対的な構成を調整することによって、異なって構成されるが、“同様の”又は類似したRFパルスが生成可能である。例えば、スピンエコーシーケンスでは、180°反転パルスが繰り返し印加される。画像取得全体に対して1つのタイプの又は固定されたRFパルスを利用する代わりに、システム10は、異なって構成される180°パルスなどのk空間に依存した方法により異なって構成されたRFパルスを利用する。これらの異なって構成されるが同様のRFパルスのそれぞれは、異なるパフォーマンスを示す可能性があり(例えば、精度、k空間軌跡、チャネル振幅、位相、使用されるBチャネルなど)、異なるSRA値を生じさせる。スキャンの一部のパルスについてより低いSAR RFパルスを利用することによって、異なって構成されるRFパルスを用いた1回のスキャンにおける平均SARは、有意な画質の犠牲なく同一に構成されたRFパルスを利用したスキャンと比較して低減可能である。このコンセプトはまた、RFシミング、RFパルスのリフォーカシング、信号サンプリング前の縦磁化を調整するのに用いられる他の磁化準備RFパルスに適用可能である。後者はまた、空間的に選択的でないが、例えば、化学シフト選択的パルスであるRFパルスを含む。この場合、利用されるRFパルスの対応するクラスは、異なるものとすることができる。
図2は、RFパルスのSARとそれの全体的な精度との間の関係のグラフ50を示す。高い精細さのRFパルスのパフォーマンスは、高いSARにSARに対応する。開示されるシステム及び方法では、RFパルスの設計は、信号サンプリング処理に関して検討される。このアプローチの1つの検討は、キーホールイメージング(例えば、van Vaals J.et al.JMR I1993;3:671−75などを参照)と動き適応化ゲーティング(例えば、Weiger M et al.MRM 1997;38:322−33などを参照)とから知られる、k空間の異なるエリアが信号不完全性に対して異なる感度を示すことである(例えば、Fuderer M.IEEE TMI 1988;7:368−80;van Vaals J.et al.JMR I1993;3:671−75;Weiger M et al.MRM 1997;38:322−33などを参照)。従って、MRスキャン全体について固定的なRFパルス構成を利用する代わりに、繰り返されるRFパルスの2以上の異なる構成が1回のスキャンにおいて利用される。これらのRFパルスのそれぞれは、異なるパフォーマンス及び/又は精度を示すものであってもよく、異なるRFパルスに固有のSAR値を生じさせる。従って、RFパルスは、実際の励起パターン、B波形及び/又はk空間軌跡などにおいて若干異なるものであってもよい。スキャンにおける平均SARは、画質を犠牲にすることなく、固定的なRFパルスを利用するスキャンと比較して低減可能である。
図3は、2D空間選択的RFパルスを利用したローカル励起の励起リードアウト60を示す。リードアウトは、RF(i),RF(i),RF(i),RF(i)とそれぞれラベル付けされた4つの個々の送信RFチャネルを用いた2Dズームスピンワープイメージングを可能にする。G及びGは直交磁場勾配の波形を表し、ADCはアナログ・デジタルコンバータであり、当該装置がシーケンス中にアクティブであるときを示す。各kについて、個々のRFパルスRF(i)は、スキャン全体の平均SARを最適化するのに利用可能である。
一例によると、デカルトスピンワープサンプリングスキームが検討される。2D RF励起パルスが、信号の発信元のエリアを制限するローカライズドMRについて利用される。利用されるRFパルスの励起フィールド(FOX)より小さい視界(FOV)においてサンプリングが実行される。
図4は、空間精細さの異なるRFパルスのクラスを用いて実行されるk空間における各ラインの信号励起の図70を示す。(左の)含まれるサンプリングFOV74(破線の四角形)を有するFOXに与えられるRFパルスパターンは、イメージングk空間76に対応する。カラーバー78は、異なる空間精細さを有するRFパルスがkサンプリングに利用されることを示す。空間精細さRF(1),RF(2),RF(3),RF(4)は、ガウスカーネルによるRFパルスターゲットの磁化をフィルタリングし(サンプリングFOVに関して0.1〜4.0のピクセルFWHM)、対応するSARを推定しながら(例えば、Graesslin I,et al.2008;ISMRM 621を参照)、Grissom et al.(MRM 2006;56:620−29)に従ってパラレル送信RFパルスを計算することによって実現される。繰り返し印加されるRFパルスは、k空間の中心近くの低精細さRF(4)から、k空間のエッジに隣接する徐々に高くなる異なる空間精細さのレベルまでの異なる空間精細さのレベル、すなわち、異なるSARレベルによって生成される。各RFパルスの対応するk空間位置へのマッピングは経験的に実行可能であるか、又は上述されたNRMSE計算アルゴリズムを利用することができる。MR信号サンプリングは、128×128のマトリックスを仮定して実行されてもよい。所与の位相符号化ステップkのための信号を生成するため何れのRFパルスが利用されるか決定するため、トライアルアンドエラーサーチが、励起パターンのk空間表現のエネルギーを考慮して実行されてもよい(図1のプロセッサ14などによって)。
1回の取得において異なって構成される同様のRFパルスの利用は、データの不整合を生じさせる可能性がある。従って、生じたエラーが通常の画像ノイズによって又はノイズフロアのアーチファクトピーキングアウト(すなわち、超過値)によって占められているか判断される。一例では、所定のSNR(15など)が仮定されてもよい。異なるRFパルス構成のパフォーマンスを判断するため、NRSMEが、不定のSNRにおいて固定された最適なRFパルス構成を利用した取得スキャンに関して計算される。
テーブル1は、わずかに増加した励起エラー(NRMSE)を犠牲にして可能なSARの減少の間の妥協を示す一例となるデータを提供する。
Figure 0005619136
所与のエラー(NRMSE)に対して、平均SAR値が、スキャンにおいて異なる又は可変的なRFパルス構成を用いたモード(SAR)と、スキャン全体に対して最良のRFパルスを利用した固定構成RFパルスモード(SAR)とについて比較される。より高いエラーが許容されるとき、より大きなSARの減少が実現され、それは、所与の及び/又は期待されるSNR(Signal−to−Noise Ratio)に依存する。より大きなNRMSE(テーブル1の0.25のケースを参照)に対応する低いSNRのケースでは、SARは33%減少される。しかしながら、上述された一例となるシナリオは、周波数符号化に関するものであるため、最適でないかもしれず、結果を不明りょうにしうる。純粋な2D位相符号化は、上述した技術の効果を改良するのに利用可能である。例えば、このコンセプトは、RFパルスのリフォーカシング及びすべてのタイプの磁化準備RFパルスに適用可能である。
図5は、k空間80の2Dサンプリングにおける循環リング構成を示す。各リング82において、RFパルス波、異なるSRA値を生じさせるため異なるコンポーネントから構成される(すなわち、異なって構成される)。RFシミングのケースでは、同じ基本RF波形とk空間軌跡とが各RFパルスについて利用可能であるが、異なる送信振幅及び位相のセットが、異なるk空間位置の各チャネルについて利用されてもよい。曲線84は、シミングに用いられる固有モードの個数を1次元(半径)的に概略的に示す。すなわち、k空間の中心で最大のシミングであり、周辺に向かってシミングが低減される。
一実施例では、2以上のパルス構成が、各タイプのパルスについて予め計算され(図1のプロセッサ14などによって)、格納される(図1のメモリ16などに)。例えば、実行されるk空間サンプリングステップに従って対応するRFパルス構成又はRFパルスを選択するコンピュータプログラム(メモリ16に格納され、プロセッサ14により実行されるなど)が提供される。異なるRFチャネルについて異なるB波形に加えて、k空間軌跡がまたスキャン中に変更される場合、対応する勾配波形を予め格納し、スキャン中にアクティブ化させることが可能である。
開示されるシステム及び方法の理解をさらに促進するため、図1のシステム10により実行可能な以下の具体例が与えられる。高フィールドシステムにおいて、3DMRIスキャンが実行される。Bシミングが、波動効果により生じるB変動を低減させるのに利用される。各RFチャネルについて、同じ基本RFパルス波形(選択勾配が存在する場合のスライス選択のsinc波形など)が伝搬されるが、RFシミングを促進するため、i個の異なるチャネルに対して振幅と位相(a,ψ)は異なる。いくつかの異なるRFシミングセット(a,ψ)がスキャン前に計算され、これらはすべて対応するSARとBシミングパフォーマンスにおいて若干異なる。
サンプリングスキームの対応する2D位相符号化空間では、異なるRFシミング設定が適用される異なる領域が規定される。本例では、k空間の中心部において、低いSAR設定が利用される。ここでは、低い空間周波数が取得され、高い空間周波数を無視するRFシミング設定の利用を正当化するためである。k空間におけるより高い位相符号化ステップについて、より良好な空間精細さを有すが、より高いSARとなるRFシミング設定が利用される。スキャン対象のオブジェクトに関する知識がない場合、図5に与えられ得るように、シンプルなRFパルスアプリケーションパターンが利用される。曲線84は、RFパルスのSARと空間精度に関する所望数のコイル固有モードを示す。開示される分布は例示的な目的のために与えられており、限定的な意味で解釈されるべきでなく、ここに開示される特徴及び態様に従って他の分布も想到しうることが理解されるであろう。
本例について、図6及び7は、RFパルスパフォーマンスとSARとを操作するための固有モードアプローチを示す。図6は、コイルがそれぞれcoil0−coil7とラベル付けされた8チャネル送信コイルにおける8つの横軸Bコイルのコイル感度90を示す。図7は、E−Coil0−E−Coil7とそれぞれラベル付けされた8チャネル送信コイル構成の固有モード100を示す。これらの固有モードはまた仮想的なコイルとみなすことが可能であり、“コイル”という用語は図7に関してその説明において使用されることが理解されるであろう。
マルチエレメントコイルの送信感度は、それらの固有モードに分解可能である。所望の精度レベルに基づき、異なる個数の固有モードがRFシミングのために検討される。各固有モードの実現は、所与のSARに対応する特定のRFパワーに関する。しかしながら、アレイの空間送信特性を取得し、SARに関連する他の仮想的なコイル分解が検討されてもよい。このような分解はまた、“モード”と呼ばれる。その方法が繰り返し実行されると、現在最も高い固有モードが省略可能であり、Bシミング問題は、残りの固有モードを用いて低減されたSARを示すシミング係数(a,ψ)を決定することによって解かれる。
他の具体例によると、図3に関して上述されたような選択的なローカル励起に対して2D RFパルスが利用される。さらに又はあるいは、異なる範囲のパルス周波数コンポーネントが、マルチチャネルトランしミットシステムにおいて磁化準備RFパルスを生成するのに利用される。
図8は、ここに開示される1以上の態様によるマルチチャネル送信コイルMRI装置12(図1)においてSARを低減する方法を示す。特定ステップ110において、MRI装置(12)により利用されるスキャンシーケンスが特定される。スキャンシーケンスは、周辺(高解像度)k空間をスキャンし、その後に中心(低解像度)k空間をスキャンしたり、中心(低解像度)k空間をスキャンし、その後に周辺(高解像度)k空間をスキャンするなど、k空間データ取得の順序を記述する情報を含む。ルックアップテーブルステップ112において、特定されたスキャンシーケンスに一致する予め生成されたパルスパターンを特定するため、テーブルルックアップが実行される。パルスパターンは、典型的には、RF準備パルス、励起パルス及び他の磁化操作パルスを含む。各RFパルスについて、低周波数コンポーネントと高周波数コンポーネントとを有する空間的に正確なパルス構成、低減された高周波数コンポーネントを有した又は有さない空間的に低い精度のパルス構成など、少なくとも2つの異なるパルス構成が格納される。他の順序がここに開示された発明の範囲の範囲内であると考えられるが、高周波数パルスが周辺のk空間データの取得中に利用され、低周波数パルスが中心のk空間データの取得中に利用される。例えば、特定されたスキャンシーケンスが周辺のk空間データを取得し、その後に中心のk空間データを取得する場合、高周波数RFパルスとその後に低周波数RFパルスを含むパルスパターンが特定される。
データ取得ステップ114において、特定されたスキャンシーケンスは、特定されたパルスパターンを用いてMRスキャンデータを取得するため実行される。画像出力ステップ116において、SARが低減されたMR画像がディスプレイ上に出力される(例えば、取得したスキャンデータの再構成後など)。このように、可変的なRFパルスが、MR画像におけるSARを低減するのに利用される。
他の実施例によると、Bフィールドシミングが、SARを低減するのに実行される。例えば、1以上の固有モードが周辺のk空間データを取得するサイにシミングアルゴリズムにおいて利用可能であり、1つの中心の固有モードが、中心のk空間データを取得する際にシミングアルゴリズムにおいて利用可能である。
所与のRFパルスのコンポーネントは、マルチトランスミットシステムの送信要素の相対的な貢献を変更し、周波数コンポーネントを変更するなどによって可変とすることができる。
本発明は、実施例を参照して説明された。上記の詳細な説明を読み、理解した者には改良及び変更が想到するかもしれない。本発明は、添付された請求項又はその均等の範囲内にある限り、このようなすべての改良及び変更を含むものとして解釈されることが意図される。請求項において、括弧の間に配置される参照符号は、請求項を限定するものとして解釈されるものでない。“有する”という用語は、請求項に列記された以外の要素又はステップの存在を排除するものでない。要素に前置される“ある”という単語は、当該要素が複数存在することを排除するものでない。開示された方法は、複数の異なる要素を有するハードウェアによって、また適切にプログラムされたコンピュータによって実現可能である。複数の手段を列記したシステムの請求項では、これらの手段のいくつかは、同一のコンピュータ可読ソフトウェア又はハードウェアにより実現可能である。特定の手段が互いに異なる従属クレームに記載されている事実は、これらの手段の組み合わせが効果的に利用可能でないことを示すものでない。

Claims (13)

  1. 低減された比吸収率(SAR)を有するMRデータを取得する方法であって、
    2以上の送信コイルを用いてMR取得スキャン中にMRスキャンデータを生成するため、少なくとも1つのRFパルスが繰り返し印加される磁気共鳴シーケンスを適用するステップと、
    前記RFパルスがいくつかの繰り返しにおいて他の繰り返しより低いSARを生じさせるように、前記繰り返し印加されるRFパルスの構成を変更するステップと、
    k空間において前記MRスキャンデータを取得するステップと、
    を有し、
    前記より低いSARのRFパルスは、より低い周波数コンポーネントを有し、
    前記その他のRFパルスは、より高い周波数コンポーネントを有する方法。
  2. 前記スキャンデータを画像に再構成するステップと、
    前記画像を表示又は格納するステップと、
    をさらに有する、請求項1記載の方法。
  3. 前記より高い周波数コンポーネントを有する繰り返し印加されるRFパルスは、k空間の周辺領域の前記MRスキャンデータを取得するのに利用され、前記より低い周波数コンポーネントを有するRFパルスは、k空間の中心領域の前記MRスキャンデータを取得するのに利用され、又は、
    前記より高い周波数コンポーネントを有する繰り返し印加されるRFパルスは、k空間の中心領域の前記MRスキャンデータを取得するのに利用され、前記より低い周波数コンポーネントを有するRFパルスは、k空間の周辺領域の前記MRスキャンデータを取得するのに利用される、請求項1記載の方法。
  4. 周辺のk空間データを取得する際に1以上の固有モードを用いてBフィールドをシミングするステップと、
    中心のk空間データを取得する際に単一の中心固有モードを用いて前記Bフィールドをシミングするステップと、
    をさらに有する、請求項1記載の方法。
  5. 前記繰り返し印加されるRFパルスのコンポーネントは、複数の送信チャネルのそれぞれからの貢献を含み、
    前記繰り返し印加されるRFパルスの構成を変更するステップは、前記貢献と前記RFパルスのk空間軌跡との間の重み付けを変更することを含む、請求項1記載の方法。
  6. 前記RFパルスは、マルチチャネルトランスミットシステムにおいて用いられる磁化準備RFパルスである、請求項1記載の方法。
  7. 前記より低いSARを有する繰り返し印加されるRFパルスは、より低い空間解像度を有し、
    前記より高いSARを有する繰り返し印加されるRFパルスは、より高い空間解像度を有する、請求項1記載の方法。
  8. 請求項1記載の方法を実行するため前記繰り返し印加されるRFパルスのバージョンを生成するRFパルス生成手段を有するマルチトランスミットチャネルMRI装置。
  9. 請求項1記載の方法を実行するようプロセッサを制御するためのソフトウェアを担持するコンピュータ可読媒体。
  10. 比吸収率(SAR)を低減する磁気共鳴システムであって、
    2以上の送信RFコイルと、
    RFパルス生成手段と、
    プロセッサと、
    ディスプレイと、
    を有し、
    前記プロセッサは、
    少なくとも1つの繰り返し印加されるRFパルスによりMRスキャンシーケンスを適用するよう前記RFパルス生成手段を制御し、
    各バージョンが異なるSARを示し、異なる周波数コンポーネントを有する前記繰り返し印加されるRFパルスの少なくとも2つの予め生成されたバージョンから選択し、
    k空間の第1領域からMRデータを取得する際には前記繰り返し印加されるRFパルスのより高い周波数のバージョンを、k空間の第2領域からMRデータを取得する際には前記繰り返し印加されるRFパルスのより低い周波数のバージョンを前記RFパルス生成手段に提供し、
    前記MRデータを処理するよう構成され、
    前記ディスプレイは、前記処理されたMRデータを表示するよう構成されるシステム。
  11. 前記RFコイルは、各RFコイルが対応するRFパルス生成手段と接続される複数のRF送信コイルを有する、請求項10記載のシステム
  12. 前記プロセッサは、
    周辺のk空間データを取得する際に1以上の固有モードを用いてBフィールドをシミングし、
    中心のk空間データを取得する際に単一の中心固有値を用いて前記Bフィールドをシミングする、
    ためのコンピュータにより実行可能な命令を実行する、請求項11記載のシステム
  13. 前記k空間の第1領域は、k空間の中心領域と周辺領域の一方であり、
    前記k空間の第2領域は、k空間の中心領域と周辺領域の他方である、請求項10記載のシステム
JP2012502867A 2009-04-02 2010-04-01 k空間に依存するRFパルス選択によるパラレル送信におけるSAR低減 Expired - Fee Related JP5619136B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP09157140 2009-04-02
EP09157140.6 2009-04-02
PCT/IB2010/051423 WO2010113131A1 (en) 2009-04-02 2010-04-01 Sar reduction in parallel transmission by k-space dependent rf pulse selection

Publications (2)

Publication Number Publication Date
JP2012522563A JP2012522563A (ja) 2012-09-27
JP5619136B2 true JP5619136B2 (ja) 2014-11-05

Family

ID=42237404

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012502867A Expired - Fee Related JP5619136B2 (ja) 2009-04-02 2010-04-01 k空間に依存するRFパルス選択によるパラレル送信におけるSAR低減

Country Status (6)

Country Link
US (1) US8890524B2 (ja)
EP (1) EP2414858A1 (ja)
JP (1) JP5619136B2 (ja)
CN (1) CN102428382B (ja)
RU (1) RU2522179C2 (ja)
WO (1) WO2010113131A1 (ja)

Families Citing this family (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5685576B2 (ja) * 2009-04-01 2015-03-18 コーニンクレッカ フィリップス エヌ ヴェ パラレルトランスミッションmriにおける時間平均によるsarホットスポットの低減
JP5619136B2 (ja) * 2009-04-02 2014-11-05 コーニンクレッカ フィリップス エヌ ヴェ k空間に依存するRFパルス選択によるパラレル送信におけるSAR低減
DE102010015044A1 (de) * 2010-04-15 2011-10-20 Siemens Aktiengesellschaft Verfahren und Einrichtung zur Ermittlung einer Magnetresonanzsystem-Ansteuersequenz und Verfahren zum Betrieb eines Magnetresonanzsystems
WO2012023098A1 (en) * 2010-08-20 2012-02-23 Koninklijke Philips Electronics N.V. Virtual coil emulation in parallel transmission mri
DE102011005433B4 (de) * 2011-03-11 2012-09-27 Siemens Aktiengesellschaft Verfahren und Hochfrequenz-Kontrolleinrichtung zur Kontrolle einer Hochfrequenzsendeeinrichtung eines Magnetresonanztomographiesystems
JP5677226B2 (ja) * 2011-07-28 2015-02-25 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴装置およびプログラム
US9304179B1 (en) * 2011-08-12 2016-04-05 University Of New Brunswick Method of magnetic resonance imaging combining phase and frequency encoding
DE102012212402B3 (de) * 2012-03-20 2013-10-17 Siemens Aktiengesellschaft Ermittlung einer MR-Sequenz unter Berücksichtigung der Energieverteilungsfunktion im k-Raum
DE102012207132B3 (de) * 2012-04-27 2013-09-05 Siemens Aktiengesellschaft Ansteuerung eines Magnetresonanzsystems unter Ermittlung von Sendeskalierungsfaktoren für eine vorgegebene Zielmagnetisierung
CN103513204B (zh) * 2012-06-29 2016-03-30 西门子(深圳)磁共振有限公司 一种磁共振成像中k空间数据的轨迹校正方法和装置
DE102012215255B3 (de) 2012-08-28 2014-02-27 Siemens Aktiengesellschaft Ansteuerung eines Magnetresonanzsystems unter Berücksichtigung aktueller komponentenbedingter B1-Feld-Maximalwerte
DE102013205868B4 (de) * 2013-04-03 2014-11-27 Albert-Ludwigs-Universität Freiburg Verfahren zum Zuordnen von K-Raum-Zeilen zu Echozügen und Verfahren zur Erfassung von MR-Daten sowie entsprechend ausgestaltete Magnetresonanzanlagen
DE102013205830B4 (de) 2013-04-03 2024-05-29 Siemens Healthineers Ag Verfahren und Bilddatenerzeugungseinrichtung zur Erzeugung von Bilddaten eines sich bewegenden Objekts, Magnetresonanzanlage und Computerprogrammprodukt
JP6001784B2 (ja) * 2013-07-09 2016-10-05 株式会社日立製作所 磁場調整支援装置、磁場調整支援方法、mri装置および磁石装置
US10761155B2 (en) * 2013-07-16 2020-09-01 Siemens Aktiengesellschaft Hybrid averaging method for specific absorption rate supervision
DE102013214867B4 (de) * 2013-07-30 2015-02-12 Siemens Aktiengesellschaft Ermittlung einer Magnetresonanz-Ansteuersequenz mit konzentrischen, kreisförmigen Sendetrajektorien
JP6101352B2 (ja) 2013-08-27 2017-03-22 株式会社日立製作所 磁気共鳴撮像装置および撮像パラメータ決定方法
CN110177499A (zh) 2016-11-21 2019-08-27 西门子医疗有限公司 通过磁共振设备记录头部的诊断测量数据的方法
US11129541B2 (en) * 2016-11-21 2021-09-28 Siemens Healthcare Gmbh Method for recording diagnostic measurement data of a head of an examination object in head imaging via a magnetic resonance device
WO2022236171A1 (en) * 2021-05-07 2022-11-10 The United States Of America, As Represented By The Secretary, Department Of Health And Human Services Method of improving uniformity of magnetization and magnetic resonance images at high static fields
DE102023202430A1 (de) * 2023-03-20 2024-09-26 Siemens Healthineers Ag Verfahren zur Ermittlung einer lokalen spezifischen Absorptionsrateninformation und zur Ermittlung einer Messpulsgruppe von Hochfrequenzpulsen, Magnetresonanzeinrichtung, Computerprogramm und elektronisch lesbarer Datenträger

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6249121B1 (en) * 1999-05-17 2001-06-19 General Electric Company RF body coil
DE10219528A1 (de) * 2002-05-02 2003-11-13 Philips Intellectual Property Schnelles Kernresonanz-Bildgebungsverfahren mit Gradienten-Echos
US6975114B1 (en) 2002-11-20 2005-12-13 Nova Medical, Inc. Methods for transmit excitation in magnetic resonance imaging using a transmit pulse with time varying spatial characteristics
US6989673B2 (en) * 2003-11-26 2006-01-24 General Electric Company Method and apparatus to reduce RF power deposition during MR data acquisition
US7053618B2 (en) 2003-11-26 2006-05-30 General Electric Company Method and apparatus to generate an RF excitation consistent with a desired excitation profile using a transmit coil array
US7619413B2 (en) * 2004-05-04 2009-11-17 The General Hospital Corporation Transmit-receive array for high field MRI
US7078900B2 (en) * 2004-07-23 2006-07-18 General Electric Company Method and system of determining parameters for MR data acquisition with real-time B1 optimization
US7436176B2 (en) * 2006-03-27 2008-10-14 Fonar Corporation Magnetic resonance imaging using adaptive phase encoding gradient
US7385396B2 (en) 2006-04-20 2008-06-10 General Electric Company SAR reduction in MR imaging with parallel RF transmission
US7336145B1 (en) * 2006-11-15 2008-02-26 Siemens Aktiengesellschaft Method for designing RF excitation pulses in magnetic resonance tomography
US7737690B2 (en) * 2007-05-18 2010-06-15 General Electric Company System and method for amplitude reduction in RF pulse design
US7619410B2 (en) * 2007-06-08 2009-11-17 General Electric Company System and method for accelerated magnetic resonance parallel imaging
DE102007043445A1 (de) * 2007-09-12 2009-04-09 Siemens Ag Verfahren zur Aufnahme einer Magnetresonanztomographiemesssequenz mit einer Magnetresonanztomographieeinrichtung sowie zugehörige Mangetresonanztomographieeinrichtung
US8148985B2 (en) * 2008-10-15 2012-04-03 Massachusetts Institute Of Technology Method for reducing maximum local specific absorption rate in magnetic resonance imaging
JP5619136B2 (ja) * 2009-04-02 2014-11-05 コーニンクレッカ フィリップス エヌ ヴェ k空間に依存するRFパルス選択によるパラレル送信におけるSAR低減

Also Published As

Publication number Publication date
US20120019247A1 (en) 2012-01-26
WO2010113131A1 (en) 2010-10-07
US8890524B2 (en) 2014-11-18
JP2012522563A (ja) 2012-09-27
EP2414858A1 (en) 2012-02-08
RU2011144368A (ru) 2013-05-10
CN102428382A (zh) 2012-04-25
CN102428382B (zh) 2015-02-25
RU2522179C2 (ru) 2014-07-10

Similar Documents

Publication Publication Date Title
JP5619136B2 (ja) k空間に依存するRFパルス選択によるパラレル送信におけるSAR低減
JP6513398B2 (ja) 事前情報に制約される正則化を使用するmr画像再構成
JP4657710B2 (ja) 前情報により向上される動的な磁気共鳴イメージング
JP5361160B2 (ja) 3次元磁気共鳴画像化用測定シーケンスおよび磁気共鳴装置
CN109814058B (zh) 磁共振系统和在其中生成检查对象的磁共振图像的方法
CN109239633B (zh) 磁共振成像方法和系统、校正扩散加权磁共振数据的方法
JP5984258B2 (ja) 磁気共鳴イメージング装置及び画像処理装置
EP3149508B1 (en) Mr imaging using multi-echo k-space acquisition
JP5525596B2 (ja) 加速b1マッピング
US20130088225A1 (en) System for Reconstructing MRI Images Acquired in Parallel
EP4071494A1 (en) Method for acquiring a three-dimensional magnetic resonance image dataset and for generating a motion-corrected image dataset
US9658304B2 (en) MRI method for retrospective motion correction with interleaved radial acquisition
US10429463B2 (en) Quiet MRI with spin echo (SE) or fast spin echo (FSE)
US20160033606A1 (en) Dixon-type water/fat separation mri using high-snr in-phse image and lower-snr at least partially out-of-phase image
JP4364789B2 (ja) 加速されたデータ収集を用いる磁気共鳴イメージング方法
US20170356972A1 (en) Quiet mri using alternating gradient sequence
US10895619B2 (en) MR imaging with Dixon-type water/fat separation
CN115598575B (zh) 磁共振成像方法、装置、计算机设备和存储介质
Chu et al. Coil compression in simultaneous multislice functional MRI with concentric ring slice‐GRAPPA and SENSE
CN101248366A (zh) 用于并行磁共振成像的设备和方法
US7285955B2 (en) System and method for improved data acquisition for medical imaging
Bammer et al. Parallel-Imaging Reconstruction of Arbitrary-k-Space-Sampling Data

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20130328

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20140328

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20140422

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20140722

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20140916

R150 Certificate of patent or registration of utility model

Ref document number: 5619136

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees