JP2009533163A - コイル感度の波動伝搬に基づく評価 - Google Patents
コイル感度の波動伝搬に基づく評価 Download PDFInfo
- Publication number
- JP2009533163A JP2009533163A JP2009505528A JP2009505528A JP2009533163A JP 2009533163 A JP2009533163 A JP 2009533163A JP 2009505528 A JP2009505528 A JP 2009505528A JP 2009505528 A JP2009505528 A JP 2009505528A JP 2009533163 A JP2009533163 A JP 2009533163A
- Authority
- JP
- Japan
- Prior art keywords
- coil
- sensitivity map
- low resolution
- data
- sensitivity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
- G01R33/32—Excitation or detection systems, e.g. using radio frequency signals
- G01R33/34—Constructional details, e.g. resonators, specially adapted to MR
- G01R33/341—Constructional details, e.g. resonators, specially adapted to MR comprising surface coils
- G01R33/3415—Constructional 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image 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/5611—Parallel 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
Landscapes
- Physics & Mathematics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
全身コイル18及びパラレルイメージングコイルの各コイル素子201,202,...,20nからの低解像度画像データが、メモリ又はバッファ34に受け取られる。再構成プロセッサ36は、低解像度全身コイルデータ及び各コイル素子からの低解像度データを、対応する低解像度画像38に再構成する。各コイル素子からの低解像度画像は、各コイル素子について対応する感度マップ441,442,...,44nを生成するために、全身コイルから低解像度画像で除算される(42)。低解像度ボディコイル画像が、ほぼゼロの値を有する領域において、又はボディコイル若しくは受信コイル画像の値が非常に急に変化する領域において、感度マップは、欠陥を有する。感度マップ補正回路又はアルゴリズム50は、感度マップの欠陥のある領域を決定し、(a)コイルジオメトリマップ56及び(b)波動伝搬モデル58に従って、感度マップの隣接する部分を内挿/外挿することにより、欠陥のある領域を補正し、視野の外側の領域に感度マップを拡げ、又は測定された感度マップを完全に置き換えて、各コイル素子について補正された感度マップを生成する。
Description
本発明は、パラレルイメージング技法に関する。本発明は、SENSEパラレルイメージング技法を使用する医用画像診断に関して特定のアプリケーションを見い出し、それを特に参照して記述される。しかしながら、本発明は、他のパラレルイメージング技法及び医用診断目的以外のイメージングにも適用できることが理解されるべきである。
磁気共鳴画像をより迅速に生成するためのさまざまなパラレルイメージング技法が、知られている。これらのパラレルイメージング技法は、SENSE、SMASH、その他を含む。磁気共鳴イメージングの間、それぞれ異なる受信コイル又は(複数の)受信コイルのそれぞれ異なるグループが、k空間のそれぞれ異なる部分を同時にサンプリングする。診断画像を再構成する際、各コイル(又はグループ)からのデータは、その感度に従って変換され又は「展開(unfolded)」される。最終画像の精度は、コイル感度を正確に決定することに依存する。
感度情報を得るために、全身受信コイル及び各受信コイル素子について大きい視野の画像データを取得する低解像度スキャンが、実施される。単一コイル画像の各々は、基準として利用できるボディコイル画像で除算される。この除算の結果は、対応する受信コイル素子又は協働する素子のグループの感度マップと考えられることができる。最終画像の精度は、感度マップの精度に依存する。決定された感度のエラーは、画像の不完全な展開に起因するいわゆる「SENSEアーチファクト」にいたることがあり、残りの信号部分は、画像アーチファクトとして現れ持続する。この問題は、SENSE加速ファクタとしてより顕著になっており、すなわち、k空間をサブサンプリングする程度が増える。
感度の不正確さにつながりうる1つの問題は、低信号密度の領域に起因する。これらの低信号密度領域において、基準スキャンの信号は、ノイズが多く、これは、不安定な又は不正確なコイル感度評価につながる。不正確さの別の原因は、基準スキャンの低解像度に起因する。原則として、コイル感度マップは、視野内の空間平滑化関数であり、それは数立方センチメートルの低解像度ボクセルサイズによって、正確にサンプリングされることができる。しかしながら、コイル素子の近傍で、コイル感度は、急勾配で上昇する。受信コイルが、患者上に又は患者の非常に近くに配置される表面コイルであるとき、低解像度の基準スキャンは、コイル近傍の急勾配の増加の領域において感度を正確に投影するには十分でない。
低信号レベルの領域に対処する1つの技法は、信号対雑音比を改善するために、信号を数回平均化することである。しかしながら、スキャン中の信号アーチファクトは、信号コンテントを支配しうる。生理的動態に関する領域において、平均化された画像が、すべての領域において正確に位置合わせされることができるというわけでない。例えば肺の領域では、心臓の動き及び血流が、ゴースト及びスメアリング(smearing)アーチファクトを生成する。線形又は他の数学的補間方法が、低信号領域におけるコイル感度を補正するために提案されている。それらは更に、外挿をある程度実施することができる。しかしながら、このような補間技法は、コイルジオメトリを考慮に入れていないので、それらは更に、内挿又は外挿された低信号領域及びコイル近傍の高感度領域の両方における不正確さに苦しむ。
本発明は、適用される内補又は外挿技法に電磁気制約を適用することによって、これら及び他の問題を解決することを目的とする。
1つの見地により、画像診断システムが提供される。補間器は、複数のパラレルイメージングコイル素子の各々について感度マップを受け取る。感度マップは、識別可能な領域に欠陥を有する。補間器は、(a)予めロードされたコイルジオメトリ及び(b)波動伝搬モデルに従って、各感度マップからのデータ、又は感度マップが生成される基礎をなすデータを補間することにより、欠陥のある領域を補正し、各コイル素子について感度マップを生成し、又は視野全体について(b)によってそれを完全に予測する。
感度マップは、
−低信号強度、
−基準スキャン及び実際のSENSEスキャンの合間又は最中の患者の動き、
により、情報が入手可能でなかった領域において受け取られる。
−低信号強度、
−基準スキャン及び実際のSENSEスキャンの合間又は最中の患者の動き、
により、情報が入手可能でなかった領域において受け取られる。
1つの利点は、パラレルイメージング技法における一層正確な展開にある。
別の利点は、SENSEアーチファクトの低減、及び非加速スキャンにおける改善されたコイル信号の組み合わせにある。
別の利点は、多数のパラレルイメージングチャネルを有するイメージングを容易にすることにある。
本発明の更なる利点が、当業者には、以下の詳細な説明を読み理解することにより認められるであろう。
本発明は、さまざまな構成要素及び構成要素の取り合わせ並びにさまざまなステップ及びステップの取り合わせの形をとりうる。図面は、好適な実施例を示す目的のみであり、本発明を制限するものとして解釈されるべきではない。
図1を参照して、MRIイメージングシステム10は、イメージング領域14に主磁界又はB0磁界を生成する1又は複数の主磁界コイル12を有する。主磁界コイルは、超電導、電気抵抗、永久磁石等でありうる。勾配コイル16は、空間符号化、周波数符号化及び位相符号化を提供するために、B0磁界を横切って勾配磁界Gx、Gy、Gzを印加する。全身送信/受信コイル18は、イメージング領域14に共鳴励起及び操作RFパルスを送信し、イメージング領域から磁気共鳴信号を受信する。
局所パラレルイメージングコイル20は、イメージング領域14に被検体と隣り合って配置される。パラレルイメージングコイルは、k空間のそれぞれ異なるサブ領域からイメージングデータを同時に生成するために、以下でコイル素子201,202,...,20nと呼ばれる、独立して又は小さいグループで機能する複数の素子又はループを有する。
シーケンスコントローラ22は、勾配磁界パルスを印加するために勾配コイルを制御する勾配増幅器24と、全身コイル18に磁気共鳴励起パルスを供給する送信器/受信器(T/R)ユニット26と、を制御する。シーケンスコントローラは更に、複数のn個のT/Rユニットを制御するために一連の送信器/受信器(T/R)ユニット281,282,...,28nを制御し、各ユニットは、パラレルイメージングRFコイル20のn個の独立して駆動可能なコイル素子の1つに関連付けられる。シーケンスコントローラ22は、他のオペレーションの中の被検体との最初のセットアップ較正中に、勾配増幅器、全身コイルT/Rユニット26及びパラレルイメージングT/Rユニット281,282,...,28nを制御して、低解像度イメージングシーケンスを実行し、感度マップを生成するためのデータを取得する。
感度マップの生成中、T/Rユニット26は、全身コイル18から共鳴信号を受け取り、復調する。T/Rユニット281,282,...,28nは、パラレルイメージングコイル20の各々の独立したコイル素子から、共鳴信号を受け取り、復調する。受け取られた共鳴信号は、イメージングデータメモリ34の個別のバッファ又は適当な部分にダウンロードされる。1又は複数の再構成プロセッサ36は、低解像度画像データを、画像メモリ38の個別の又は対応するセクションに記憶される対応する一連の低解像度画像に再構成する。平滑化関数40が、画像を平滑化する。除算器42は、ピクセルごとに、パラレルイメージングコイル20の各々のコイル素子からの低解像度画像を、全身画像で除算して、感度マップメモリ46の適当な部分に記憶されるnコイル素子についての対応する感度マップ441,442,...,44nを生成する。こうして生成されたこれらの感度マップは、以下、「ゴールド標準」マップとも呼ばれる。
感度マップ補正回路又はアルゴリズム50は、信号対雑音比が許容できないほど低い領域又は感度勾配変化のレートが予め選択されたレートよりも高い領域を決定するために、低解像度画像又はデータを調べるアルゴリズム又はプロセッサ52を有する。各々の識別された領域について、補間器54は、全身コイル及びパラレルイメージングコイル20の各々の独立したコイル素子のコイル構成及び電流特性を予めロードされている幾何学的パラメータメモリ56からの対応するコイルのコイルジオメトリパラメータに従って、許容できる信号対雑音比又は変化レートの領域内の隣りのボクセル又はピクセルから、感度又は画像値を補間する。補間器は更に、マクスウェル方程式又は波動モデルメモリ58からの他の波動伝搬モデルに従って、コイル感度を補間する。この幾何学的パラメータ及び電磁気制約を伴うマクスウェル方程式に基づく内挿及び外挿の結果は、画像の対応する低い信号対雑音比又は急に変化する感度勾配の領域を、補間された値と置き換えるために、低共鳴画像メモリ38に戻される。代替として、補間された領域は、比率で示され、感度マップにおいて直接置き換えられることができ、又は完全な感度マップが、波動伝搬結果と完全に置き換えられることができる。更に別の代替例として、より詳しく以下に説明されるように、ジオメトリ及び波動伝搬に基づく感度と元の感度とが、互いに調和(適合)され又は融合される。内挿/外挿の詳細は、以下に記述される。更なる改良として、感度マップは、コイル素子結合を反映するように調整されることができる。補正された感度マップは、患者ごとに調整される汎用的なコイル感度を記述するために一組の基底関数を規定することによって、生成されることができる。
パラレルイメージングシーケンスの間、パラレルイメージングコイル20の各コイル素子から受け取られるデータは、対応する高解像度サブ画像を生成するために、再構成プロセッサ又はコンピュータアルゴリズム60(36と同じでありうる)によって再構成され、展開プロセッサ62は、コイル感度マップメモリ44からの対応する補正されたコイル感度マップにより、各コイル素子からのサブ画像を展開し又は変換し、それらを合わせて、画像メモリ64に記憶される画像を生成する。この技法は、可動テーブルイメージング技法の場合及びパラレルイメージング加速が実施されない場合の簡単なコイル信号組み合わせについても有用である。ビデオプロセッサ66は、再構成された1又は複数の画像の部分を選択し、後処理エンハンスメント等を実施し、モニタ68又は他の人間可読のディスプレイ上の表示の生成を制御する。ビデオプロセッサは更に、将来の取り出しのために、患者レコードデータベースへの再構成されたイメージの転送を制御する。
一般に、それぞれ異なる受信コイルについての異なるすべての画像コンテントは、そのコイル感度に関連し、それは複合関数である。同様に複雑な基礎をなす解剖学的情報は、受信コイル及びボディコイルに関して同一である。これは、コイル感度が解剖学的情報から独立していることを意味するが、それは、より一般的な態様で患者によって影響を及ぼされうる。
コイル素子iによって取得される空間ドメインのボクセルの信号は、寄与、
ci=SiTQBCρ (1)
を含み、ここで、ciは、コイル素子iの総信号を表し、Siは、ボクセル位置における対応する感度であり、TQBCは、直交ボディコイル18のRF送信によって導かれるすべての影響を合わせたものである。基礎をなす解剖学的情報は、ボクセル密度ρによって指定される。これらの異なる構成要素を分離するために、同じスキャンが、信号受信のためにボディコイル18を使用してほぼ同時に取得される。コイル感度を得るために、信号ciは、除算器回路又はアルゴリズム42において、ボディコイル信号cQBCで除算される:
これは、
にいたる。
ci=SiTQBCρ (1)
を含み、ここで、ciは、コイル素子iの総信号を表し、Siは、ボクセル位置における対応する感度であり、TQBCは、直交ボディコイル18のRF送信によって導かれるすべての影響を合わせたものである。基礎をなす解剖学的情報は、ボクセル密度ρによって指定される。これらの異なる構成要素を分離するために、同じスキャンが、信号受信のためにボディコイル18を使用してほぼ同時に取得される。コイル感度を得るために、信号ciは、除算器回路又はアルゴリズム42において、ボディコイル信号cQBCで除算される:
これは、
にいたる。
これは、適用されるコイル感度、
が、ボディコイル感度の逆数によって重み付けされることを意味する。これは、一般に重要でない。ボディコイルの感度は、大きさ及び位相の両方において一定であると考えられ、アーチファクトのないSENSE再構成を可能にするからである。上述したように、一般に、コイル感度は、好適な基準スキャンのために低解像度のみを必要とする平滑化関数である。単一コイル画像及び全身コイル基準画像の両方の空間的な平滑化は、除算演算の前に適用されるcos2フィルタ40によって実施されることができる。これは、調整によって更に安定にされることができる。全身コイルからの低解像度画像で各コイル素子からの低解像度画像を除算することによって生成される感度マップは、「ゴールド標準」感度と呼ばれる。
が、ボディコイル感度の逆数によって重み付けされることを意味する。これは、一般に重要でない。ボディコイルの感度は、大きさ及び位相の両方において一定であると考えられ、アーチファクトのないSENSE再構成を可能にするからである。上述したように、一般に、コイル感度は、好適な基準スキャンのために低解像度のみを必要とする平滑化関数である。単一コイル画像及び全身コイル基準画像の両方の空間的な平滑化は、除算演算の前に適用されるcos2フィルタ40によって実施されることができる。これは、調整によって更に安定にされることができる。全身コイルからの低解像度画像で各コイル素子からの低解像度画像を除算することによって生成される感度マップは、「ゴールド標準」感度と呼ばれる。
コイル感度は、より理論的且つ汎用的なやり方で記述されることもできる。受信コイル素子が、特定の空間位置でどれくらい感受性があるかは、受信コイル素子の特性である。相反定理を使用して、受信コイルの感度は、コイル素子の単位電流によって生成されるその横断方向のH磁界に比例する:
これらの横断方向の磁界は、コイルジオメトリ、及びこのケースでは被検体の身体である指定された媒体中の波動伝播に依存する。説明を簡潔にするため、磁界は、ベクトルポテンシャル、
の回転ベクトルによって記述されることができる:
の回転ベクトルによって記述されることができる:
このベクトルポテンシャル、
は、コイルの導体の電流密度である磁界のソースによって誘起される。導体の極小の幅及び高さを仮定すると、この密度は、導体位置、
に沿って複素電流Iとして規定されることができる:
ここで、
は、磁界が計算される視野内のポイントをあらわす。これは、ベクトルポテンシャルの寄与が、コイル導体に沿って積分されることを意味する。
は、コイルの導体の電流密度である磁界のソースによって誘起される。導体の極小の幅及び高さを仮定すると、この密度は、導体位置、
に沿って複素電流Iとして規定されることができる:
ここで、
は、磁界が計算される視野内のポイントをあらわす。これは、ベクトルポテンシャルの寄与が、コイル導体に沿って積分されることを意味する。
式(7)と組み合わされて、これは、複素感度成分に関する最終的な式に至る:
ここで、
であり、ダイポールnの長さ
は、Δxn、Δyn、Δznに分解される。波数kは、複素数の性質をもつ。これらの一般的な式は、コイル感度を評価するために以下において使用される。数値最適化技法(例えばシンプレックスアルゴリズム)を使用して、式13(各々11−12)によって議論されるモデルに与えられる自由パラメータは、測定されたコイル感度に調和される。適当なペナルティ関数は、このモデルを使用して以後計算される磁界を、いくつかの補間ポイントにおいてゴールド標準感度と比較する。コイル位置及び装荷のようないくつかのパラメータを調整することによってこれらの補間ポイントの変量を最小限にすることにより、正確な感度評価が見つけられることができる。システム全体にわたる内挿及び外挿されたコイル感度マップが、生成されることができる。同様に、リジッドなアレイを有する対称コイルの場合、コイル感度マップ計算が、簡略化されうる。
ここで、
であり、ダイポールnの長さ
は、Δxn、Δyn、Δznに分解される。波数kは、複素数の性質をもつ。これらの一般的な式は、コイル感度を評価するために以下において使用される。数値最適化技法(例えばシンプレックスアルゴリズム)を使用して、式13(各々11−12)によって議論されるモデルに与えられる自由パラメータは、測定されたコイル感度に調和される。適当なペナルティ関数は、このモデルを使用して以後計算される磁界を、いくつかの補間ポイントにおいてゴールド標準感度と比較する。コイル位置及び装荷のようないくつかのパラメータを調整することによってこれらの補間ポイントの変量を最小限にすることにより、正確な感度評価が見つけられることができる。システム全体にわたる内挿及び外挿されたコイル感度マップが、生成されることができる。同様に、リジッドなアレイを有する対称コイルの場合、コイル感度マップ計算が、簡略化されうる。
式(4)、(10)及び(11)は、その形状及び幾何学的セットアップに基づいて受信コイル素子のコイル感度を規定するが、人体の物体特性によっても影響される。複素波数kは、人体の波動伝搬経路に依存する伝搬特性を含む。より具体的には、導電率σ及び誘電率εは、以下の相関関係によってカバーされる:
k 2=εμω2−jσμω (12)
k 2=εμω2−jσμω (12)
人体の10−5以下の透磁性率μの変量は、波動伝播に関して無視できるほどごくわずかであり、μ0によって置き換えられることができる。回転フレーム周波数を表すωは、既知の量であり、定数と考えられることができる。結果的に、人体に関する特性は、複素波数kによって包括的に考慮に入れられる。しかしながら、より複雑なモデルが、考えられる。人体の場合、水の波動伝搬数に基づく複素波数kが、数値最適化ルーチンの開始値として使用されることができる。
包括的なスケーリングパラメータとして考えられることにより、受信コイル素子に誘起される複素電流は、除去されることができる。モデルにより規定される電流が、微小な幅の導体を流れる。電流が小さい導体に沿って一定であると考えることは、未知数の数を1つの複素パラメータI 0にする近似である。小さいコイル素子の場合は特に、この状態は、明確に規定され、その場合、コイル電流の位相は、小さい受信ループに沿って一定である。より大きいコイルの場合、改善された電流分布モデルが、より正確になりうる。
図2を参照して、各々のコイル素子レイアウトの形状(互いに対するダイポールの相対位置)は、アプリオリ知識から知られる。電流を運ぶ導体の座標(xn,yn,zn)は、3つの並進パラメータ及び3つの回転パラメータによって記述される:質量中心x0,y0,z0及び角度φx,φy,φz。
メモリ56に記憶される受信コイル素子の幾何学的機構は、コイル感度分布に対して支配的な影響をもつ。パラレルイメージングRFコイル20が、ボア内の既知の位置に組み込まれており又は固定的に配置される場合、感度評価は、直接的で簡単であり、計算努力が大幅に低減される。しかしながら、自由に位置付け可能なコイルの場合でも、感度分布は、なお決定可能である。
各々の素子について別個に、図2のパラメータが、最適化される:
(a)その中央点の絶対位置は、x0,y0,z0によって記述される;
(b)その回転は、3つの軸φx,φx,φxに沿って規定される。
(a)その中央点の絶対位置は、x0,y0,z0によって記述される;
(b)その回転は、3つの軸φx,φx,φxに沿って規定される。
これは、パラレルイメージングコイル20の各コイル素子についてリジッドなジオメトリアプローチを示す。コイル素子の形状は知られており、その位置及び向きパラメータは、上記(a)及び(b)によって評価される。リジッドなアレイのコイル素子相互依存を考慮に入れることは、パラメータ評価プロシージャを改善し、簡略化することができる。
要するに、各コイル素子について評価されるパラメータは、6つの幾何学的パラメータ位置及び角度(C1−C6)と、包括的な身体特性を表す複素波動伝搬数kに関する2つのパラメータ(C7、C8)と、電流I 0に起因するコイル感度の包括的な振幅及び位相に関する2つの包括的なスケーリングパラメータ(C9、C10)と、それぞれ異なる指数距離タームを切り離す2つの付加のパラメータ(C11、C12)とを有する。これまでに言及されていない最後の2つのパラメータ(C11、C12)は、理論的な方程式(4)、(10)及び(11)と比較して、基準データとのより良い呼応を提供する。パラメータ(C11、C12)は、より局所化された効果を補償するようにみえる。実際のパラメータ化された方程式は、次のように実現される:
ここで、(xn,yn,zn)=f(n,C1,C2,C3,C4,C5,C6)であり、fは、受信コイル素子の形状に依存する。
ここで、(xn,yn,zn)=f(n,C1,C2,C3,C4,C5,C6)であり、fは、受信コイル素子の形状に依存する。
一実施例において、上述されたゴールド標準モデル及び上述の波動伝搬アプローチは、最適化プロセスにおいて組み合わせられる。ゴールド標準感度は、或る領域において不適当でありうる。主として、この問題は、ボディコイルの低い信号に関連し、その一方で、高いボディコイル信号の領域は、安定であり正確な感度評価を示す。結果的に、高信号レベルを有する低解像度の基準スキャンのポイントは、感度評価のための補間ポイントとして使用される。
感度評価は、波動伝搬アプローチを使用してこれらの補間ポイントにおいて計算される。上述したパラメータは、測定された感度と評価された感度との間の相違を最小限にするために適当な最適化戦略であるシンプレックス法を使用して調整される。
いかなる最適化に関しても、開始値が重要である。事前に知られていない場合、コイル位置についての開始値は、各受信コイル素子の最大補間ポイントによって記述される。角度についての開始値は、アプリケーション依存でありうる一般のコイル配置から得られる。最終的に、評価は、スケーリング及びローディングパラメータ(例えばεr=10,σ=0.4S/m)についての開始点として使用される。
ゴールド標準方法は、補間ポイントにおいて基準値を得るために使用される。受信コイル素子の低解像度画像は、ボディコイル基準で除算される(42)。高いボディコイル信号を有するボクセルのみを使用することにより、安定した値のみが、補間ポイントとして使用される。結果として得られるコイル感度評価が、計算される。
図3に示される画像を生成するために、直径110mmの2つのループコイルが、ファントム上のスキャナに固定され、第3のコイルが、ファントム下に配置された。基準スキャンの後、3つの同一の高解像度ボリュームが、頭部方向の0、10及び40mmのテーブルオフセット位置において取得された。これらのオフセットは、テーブルをこれらの対応する距離だけ移動させることによって達成された。これらのボリュームを再構成するために、コイル素子の受信感度は、ファントムと共にシフトされ、他方、第1及び第2のコイルの受信感度は、固定のままだった。スキャナに固定される第2のコイル素子の感度は、図3の3つの異なる位置に示されている。図3a−図3cは、ゴールド標準コイル感度を示し、図3d−図3fは、波動伝搬モデルによって評価される感度を示している。ゴールド標準感度について上述した問題は、図3cにおいて最も良く見られることができ、この図は、非常に不安定なコイル感度を示している。これは、スライスのオフセット位置から説明されることができる。基準スキャンの間、非常に低い信号のみが位置Δz=40mmにおいて受け取られた。これは、この不安定な感度規定に至る。
図3d−図3fに示されるモデルに基づいて評価される感度は、あらゆる位置において、ファントムの外側においても、安定した感度を生成することが可能である。図3e及び図3fに示される暗いドットは、現実の物理的なバックグラウンドを有する。コイル素子の近傍の領域において、このコイルは、横方向磁化に対して感受性がなく、その結果、ほぼゼロの感度を生じさせる。このポイントは、単一コイル画像においても見られることができるが、通常は隣り合った素子によって補償され、それにより、最終の再構成された画像では見えなくされる。
それぞれ異なるオフセットに関して取得される対応するスライスの再構成が、図4に示されている。ゴールド標準感度は、図4a−図4cに示される再構成に関して使用された。評価された感度を適用する対応する再構成が、図4d−図4fに示されている。基準スキャンとSENSE加速スキャンとの間の小さいオフセットは、図4bに示されるゴールド標準アプローチによって、なお補償されることができる。しかしながら、基準スキャンのホールにおける欠如する情報は、基準スキャンと画像スキャンとの間のより大きいオフセットを許さず、これは、図4cの深刻な再構成アーチファクトにいたる。評価されたコイル感度で全領域をカバーすることにより、問題は、図4fには存在せず、これは、患者テーブルの動きにかかわらず高品質な画像再構成を可能にする。
上述のプロセスは、さまざまな回路、構成要素、プロセッサアルゴリズム、全てが集合的にコンピュータプログラマブル媒体と示されるディスク又は他の電子記録媒体等に記憶されるコンピュータプログラムによって実施されることができる。
本発明は、好適な実施例に関して記述された。当業者であれば、変形及び変更が、前述の詳細な説明を読み理解することにより思い付くであろう。本発明は、すべてのこのような変形及び変更が添付の特許請求の範囲又はそれと等価なものの範囲内にある限り、それらを含むものとして構成されることが意図される。
Claims (21)
- 視野内の関心領域をイメージングする画像診断システムであって、
複数のパラレルイメージングコイル素子の各々について感度マップを受け取る補間器を有し、前記感度マップは、識別可能な領域に欠陥を有し、前記補間器は、コイルジオメトリ及び波動伝搬モデルに従って、各感度マップからのデータ又は該感度マップが生成される基礎をなすデータを内挿/外挿することにより、前記欠陥のある領域を補正し又は前記視野の外側の領域に前記感度マップを拡げて、各コイル素子について補正された感度マップを生成する、画像診断システム。 - 全身コイルから低解像度画像データを受け取るメモリ又はメモリ部分と、
各コイル素子から低解像度画像データを受け取るメモリ又はメモリ部分と、
前記低解像度全身コイルデータを低解像度全身コイル画像表現に再構成するとともに、各コイル素子からの前記低解像度データを各コイル素子について対応する低解像度画像表現に再構成する再構成プロセッサ又はアルゴリズムと、
前記欠陥のある感度マップを生成するために、前記低解像度コイル素子画像を前記低解像度全身コイル低解像度画像で除算する除算器であって、1つの欠陥のある感度マップが、各コイル素子に対応する、除算器と、
を更に有する、請求項1に記載のシステム。 - 各コイル素子から、k空間のそれぞれ異なるサブ領域を表す高解像度データを受け取り、各素子についての前記高解像度データを、対応する部分画像に再構成する再構成プロセッサ又はアルゴリズムと、
対応する前記補正された感度マップに従って各々の前記部分画像を展開し、前記展開された部分画像を合わせて高解像度画像表現にする展開プロセッサ又はアルゴリズムと、
を更に有する、請求項2に記載のシステム。 - 検査領域に主磁界を生成する主磁石と、
前記主磁界に勾配磁界を生成する勾配磁界コイルと、
少なくとも全身コイル低解像度画像データを生成する全身ラジオ周波数コイルと、
対応する低解像度データ及び対応する高解像度データを各々が生成する複数のコイル素子を有するパラレルイメージングコイルと、
を更に有する、請求項3に記載のシステム。 - 検査領域に主磁界を生成する主磁石と、
前記主磁界に勾配磁界を生成する勾配磁界コイルと、
少なくとも全身コイル低解像度画像データを生成する全身ラジオ周波数コイルと、
k空間のそれぞれ異なるサブ領域を表す対応する低解像度画像データ及び対応する高解像度画像データを各々が生成する複数のコイル素子を有するパラレルイメージングコイルと、
を更に有する、請求項1に記載のシステム。 - 各コイル素子から高解像度データを受け取り、各素子についての前記高解像度データを、対応する部分画像に再構成し、対応するコイル素子についての前記補正された感度マップに従って、各々の前記部分画像を展開し、前記展開された部分画像を合わせて高解像度画像表現にするプロセッサ又はアルゴリズムを更に有する、請求項5に記載のシステム。
- 前記波動伝搬モデルは、均質組織の波動伝搬に関するマックスウェル方程式に基づく、請求項1に記載のシステム。
- 前記補間器は、前記欠陥のある感度マップと前記波動伝搬モデル及び前記コイルジオメトリマップに基づいて生成される感度マップとの間の最小適合解析を実施する、請求項1に記載のシステム。
- 複数の撮像コイル素子の各々について感度マップを受け取る手段であって、前記感度マップが、識別可能な領域に欠陥を有する、手段と、
コイルジオメトリ及び波動伝搬モデルに従って、各感度マップからのデータ又は該感度マップが生成される基礎をなすデータを内挿/外挿することにより、前記欠陥のある領域を補正して、各コイル素子について補正された感度マップを生成する手段と、
を有する画像診断システム。 - 複数の撮像コイル素子の各々について感度マップを受け取るステップであって、前記感度マップが、識別可能な領域に欠陥を有する、ステップと、
コイルジオメトリ及び波動伝搬モデルに従って、各感度マップからのデータ又は該感度マップが生成される基礎をなすデータを内挿/外挿することにより、前記欠陥のある領域を補正して、各コイル素子について補正された感度マップを生成する内挿/外挿ステップと、
を含む画像診断方法。 - 前記内挿/外挿ステップは、前記波動伝搬モデルに従って、感度データを視野の周辺領域に拡げることを含む、請求項10に記載の方法。
- 全身コイルから低解像度画像データを受け取るステップと、
各コイル素子から低解像度画像データを受け取るステップと、
前記低解像度全身コイルデータを、低解像度全身コイル画像表現に再構成するステップと、
各コイル素子からの前記低解像度データを、各コイル素子について対応する低解像度画像表現に再構成するステップと、
前記欠陥のある感度マップを生成するために、前記低解像度コイル素子画像を前記低解像度全身コイル画像で除算するステップであって、1つの欠陥のある感度マップが、各コイル素子に対応する、ステップと、
を更に含む、請求項10に記載の方法。 - 各コイル素子からの高解像度データを対応する部分画像に再構成するステップと、
対応する前記補正された感度マップにより、各部分画像を展開するステップと、
前記展開された部分画像を組み合わせて高解像度画像表現にするステップと、
を更に含む、請求項10に記載の方法。 - 検査領域に主磁界を生成するステップと、
前記主磁界に勾配磁界を生成するステップと、
全身低解像度画像データを生成するステップと、
低解像度データ及び前記高解像度データを生成するステップであって、各パラレルイメージングコイル素子からのデータが、k空間のそれぞれ異なるサブ領域を表す、ステップと、
を更に含む、請求項13に記載の方法。 - 検査領域に主磁界を生成するステップと、
前記検査領域において前記主磁界に勾配磁界を生成するステップと、
全身コイルによって全身コイル低解像度画像データを生成するステップと、
複数のコイル素子から低解像度画像データを生成するステップと、
前記複数のコイル素子により高解像度画像データを生成するステップであって、各コイル素子からの高解像度画像データが、k空間のそれぞれ異なるサブ領域を表す、ステップと、
を更に含む、請求項10に記載の方法。 - 各コイル素子についての前記高解像度画像データを、対応する部分画像に再構成するステップと、
対応するコイル素子についての前記補正された感度マップに従って、各部分画像を展開するステップと、
前記展開された部分画像を組み合わせて高解像度画像表現にするステップと、
を更に含む、請求項15に記載の方法。 - 前記波動伝搬モデルは、マックスウェル方程式に基づく、請求項10に記載の方法。
- 前記欠陥のある感度マップはそれぞれ、前記波動伝搬モデル及び前記コイル感度マップに基づく各コイル素子に対応する感度マップと調和される、請求項10に記載の方法。
- 前記感度マップの前記欠陥のある領域は、前記感度マップが低解像度画像データから生成される該低解像度画像データの信号強度が、非常に低く、又は前記感度マップが生成される基礎をなすデータの値が、急に変化する、領域に対応する、請求項10に記載の方法。
- 請求項10に記載の方法を実施するようにプログラムされたコンピュータ媒体。
- コイルジオメトリ及び波動伝搬モデルに従って、波動伝搬感度マップを拡げるステップと、
低解像度画像データから欠陥のある感度マップを生成するステップと、
前記波動伝搬感度マップ及び前記欠陥のある感度マップを調和させて、補正された感度マップを生成するステップと、
を含む画像診断方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US74475506P | 2006-04-13 | 2006-04-13 | |
PCT/US2007/064194 WO2007121023A1 (en) | 2006-04-13 | 2007-03-16 | Wave-propagation based estimation of coil sensitivities |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2009533163A true JP2009533163A (ja) | 2009-09-17 |
Family
ID=38293329
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2009505528A Pending JP2009533163A (ja) | 2006-04-13 | 2007-03-16 | コイル感度の波動伝搬に基づく評価 |
Country Status (7)
Country | Link |
---|---|
US (1) | US20090278536A1 (ja) |
EP (1) | EP2010930A1 (ja) |
JP (1) | JP2009533163A (ja) |
CN (1) | CN101484821A (ja) |
RU (1) | RU2008144711A (ja) |
TW (1) | TW200745586A (ja) |
WO (1) | WO2007121023A1 (ja) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9036884B2 (en) * | 2009-09-17 | 2015-05-19 | Koninklijke Philips N.V. | Image intensity correction for magnetic resonance imaging |
JP5698554B2 (ja) * | 2010-03-05 | 2015-04-08 | 株式会社東芝 | 磁気共鳴イメージング装置 |
US8515696B2 (en) * | 2011-01-20 | 2013-08-20 | General Electric Company | Method and system for acquiring and analyzing coil data |
BR112013024733A2 (pt) * | 2011-03-28 | 2016-12-27 | Koninkl Philips Nv | método para processamento de uma ou mais imagens reconstruídas e sistema |
US9146293B2 (en) | 2012-02-27 | 2015-09-29 | Ohio State Innovation Foundation | Methods and apparatus for accurate characterization of signal coil receiver sensitivity in magnetic resonance imaging (MRI) |
BR112015000875A2 (pt) * | 2012-07-16 | 2017-06-27 | Koninklijke Philips Nv | meio de armazenamento não temporário que armazena instruções executáveis por um dispositivo de processamento de dados eletrônicos para realizar um método, e método |
DE102015203306A1 (de) * | 2015-02-24 | 2016-08-25 | Siemens Healthcare Gmbh | Verfahren zur Ermittlung von absoluten Empfangssensitivitätskarten für Empfangsspulen in einer Magnetresonanzeinrichtung und Magnetresonanzeinrichtung |
WO2019070848A1 (en) | 2017-10-06 | 2019-04-11 | University Of Cincinnati | SYSTEMS AND METHODS FOR ESTIMATING COMPLEX RADIO FREQUENCY FIELDS IN MAGNETIC RESONANCE IMAGING |
US10859646B2 (en) * | 2018-05-31 | 2020-12-08 | General Electric Company | Method and systems for coil selection in magnetic resonance imaging to reduce annefact artifact |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002056767A1 (fr) * | 2001-01-19 | 2002-07-25 | Kabushiki Kaisha Toshiba | Imagerie parallele par resonance magnetique utilisant une carte de sensibilite de bobine haute precision |
JP2004049907A (ja) * | 2002-07-01 | 2004-02-19 | Ge Medical Systems Global Technology Co Llc | 感度エンコーディング磁気共鳴イメージングの方法及びシステム |
JP2005131411A (ja) * | 2003-10-31 | 2005-05-26 | Ge Medical Systems Global Technology Co Llc | コイル感度プロフィールを較正するためのシステム及び方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7385397B2 (en) * | 2004-02-26 | 2008-06-10 | Kabushiki Kaisha Toshiba | Magnetic resonance imaging device and data processing method for magnetic resonance imaging device |
US7397242B2 (en) * | 2005-10-27 | 2008-07-08 | Wisconsin Alumni Research Foundation | Parallel magnetic resonance imaging method using a radial acquisition trajectory |
-
2007
- 2007-03-16 US US12/296,937 patent/US20090278536A1/en not_active Abandoned
- 2007-03-16 WO PCT/US2007/064194 patent/WO2007121023A1/en active Application Filing
- 2007-03-16 CN CNA2007800131364A patent/CN101484821A/zh active Pending
- 2007-03-16 EP EP07758715A patent/EP2010930A1/en not_active Withdrawn
- 2007-03-16 JP JP2009505528A patent/JP2009533163A/ja active Pending
- 2007-03-16 RU RU2008144711/28A patent/RU2008144711A/ru not_active Application Discontinuation
- 2007-04-10 TW TW096112553A patent/TW200745586A/zh unknown
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002056767A1 (fr) * | 2001-01-19 | 2002-07-25 | Kabushiki Kaisha Toshiba | Imagerie parallele par resonance magnetique utilisant une carte de sensibilite de bobine haute precision |
JP2004049907A (ja) * | 2002-07-01 | 2004-02-19 | Ge Medical Systems Global Technology Co Llc | 感度エンコーディング磁気共鳴イメージングの方法及びシステム |
JP2005131411A (ja) * | 2003-10-31 | 2005-05-26 | Ge Medical Systems Global Technology Co Llc | コイル感度プロフィールを較正するためのシステム及び方法 |
Non-Patent Citations (1)
Title |
---|
JPN6012026855; WINKELMANN R., ET AL.: '"Wave-propagation based estimation of coil sensitivities"' PROCEEDINGS OF THE INTERNATIONAL SOCIETY FOR MAGNETIC RESONANCE IN MEDICINE, 14TH MEETING PROCEEDING , 20060506, p2473 * |
Also Published As
Publication number | Publication date |
---|---|
CN101484821A (zh) | 2009-07-15 |
TW200745586A (en) | 2007-12-16 |
EP2010930A1 (en) | 2009-01-07 |
WO2007121023A1 (en) | 2007-10-25 |
RU2008144711A (ru) | 2010-05-20 |
US20090278536A1 (en) | 2009-11-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2009533163A (ja) | コイル感度の波動伝搬に基づく評価 | |
US7358732B2 (en) | System, method, software arrangement and computer-accessible medium for providing real-time motion correction by utilizing clover leaf navigators | |
JP6283009B2 (ja) | 磁気共鳴イメージングシステム | |
JP4034654B2 (ja) | 高精度コイル感度マップを用いたパラレルmrイメージング | |
JP5865262B2 (ja) | 電気特性断層画像化方法及びシステム | |
JP3976684B2 (ja) | 画像における動きの影響を低減する方法および装置 | |
US9797974B2 (en) | Nonrigid motion correction in 3D using autofocusing with localized linear translations | |
JP4981896B2 (ja) | 電気特性トモグラフィのための電界シミング | |
Manke et al. | Respiratory motion in coronary magnetic resonance angiography: a comparison of different motion models | |
CN105143906B (zh) | 金属抗mr成像 | |
EP1385018B1 (en) | Correcting geometry and intensity distortions in MR data | |
CN112596013A (zh) | 用于从校准图像生成定位器扫描设置的系统和方法 | |
WO2006106448A2 (en) | Sense mr parallel imaging with continuously moving bed | |
KR20030064641A (ko) | 측정 필드가 자동 조절되는 자기 공명 이미징 방법 | |
US7218107B2 (en) | Adaptive image homogeneity correction for high field magnetic resonance imaging | |
JP6912603B2 (ja) | デュアルエコーディクソン型水/脂肪分離mr撮像 | |
US9470772B2 (en) | Magnetic resonance image distortion correction with a reference volume | |
Honal et al. | Compensation of breathing motion artifacts for MRI with continuously moving table | |
Noterdaeme et al. | Intensity correction with a pair of spoiled gradient recalled echo images | |
WO2018028855A1 (en) | Phase offset determination in magnetic resonance imaging | |
JP6783642B2 (ja) | 磁気共鳴イメージング装置及び受信コイルの感度マップ算出方法 | |
JP2007511269A (ja) | 径方向磁気共鳴イメージングデータの収集における読み出し順序付け装置及び方法 | |
KR20230078217A (ko) | 자기 공명 영상으로부터 뇌 전기전도도 지도를 재구성하는 방법 및 장치 | |
CN115886780A (zh) | 在磁共振成像中确定外轮廓 | |
WO2019232518A1 (en) | Methods for uniform reconstruction of multi-channel surface-coil magnetic resonance data without use of a reference scan |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20100311 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20120529 |
|
A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20121023 |