JP2007518485A - 正規化可変密度センス - Google Patents
正規化可変密度センス Download PDFInfo
- Publication number
- JP2007518485A JP2007518485A JP2006548510A JP2006548510A JP2007518485A JP 2007518485 A JP2007518485 A JP 2007518485A JP 2006548510 A JP2006548510 A JP 2006548510A JP 2006548510 A JP2006548510 A JP 2006548510A JP 2007518485 A JP2007518485 A JP 2007518485A
- Authority
- JP
- Japan
- Prior art keywords
- magnetic resonance
- resonance imaging
- space
- encoded data
- imaging method
- 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.)
- Withdrawn
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/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
- G01R33/4824—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory
-
- 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
-
- 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/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
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
磁気共鳴撮像装置は、k空間の中心から少なくとも離れたところでアンダーサンプルされた可変密度感度エンコードデータを取得する複数の高周波コイルを有する。各々のコイルのための再構成処理器(52)は、そのコイルにより取得されたk空間の中心又はその中心近傍に配置されている可変密度感度エンコードデータの高密度部分と、そのコイルにより取得された可変密度感度エンコードデータから再構成された折り畳み画像とを再構成する。展開処理器(66)は折り畳み画像を展開する。その展開長谷行きか画像により正規化される。
Description
本発明は高周波技術に関する。本発明は、特に、k空間の非直交座標系サンプリングを用いる磁気共鳴撮像及び可変密度感度エンコード(SENSE)を用いる磁気共鳴撮像におけるアプリケーションに適合し、その特定の見本について詳述する。しかしながら、本発明はまた、例えば、磁気共鳴分光法のような他の磁気共鳴アプリケーションにまた、適合する。
SENSE磁気共鳴撮像技術において、複数のコイルが、共通の撮像スライス又はボリュームの磁気共鳴撮像データを同時に取得するように用いられる。SENSE技術を用いて高速で撮像データを取得するように、各々のコイルはk空間をアンダーサンプルする。各々のコイルにより取得されたアンダーサンプルされれた撮像データは、スライス又はボリュームの対応する折り畳み画像に再構成される。折り畳み画像は、スライス又はボリュームの展開画像を生成するようにコイルの感度特性に基づいて展開される。SENSE撮像技術については、例えば、文献、Preessmann et al.,Magnetic Resonance in Medicine 42,pp.952−962(1999)及びHajnal等による米国特許第6,380,741号明細書に記載されている。
従来のSENSE撮像においては、k空間は一様にサンプリングされる。可変密度SENSEにおいては、サンプリングは一様ではなく、少数のk空間サンプルはk空間の中心において又はその中心近傍においてスキップされる一方、多数のk空間サンプルはk空間の中心から離れたところでスキップされる。従来のSENSEと同様に、各々のコイルにより取得されたアンダーサンプルされた撮像データは対応する折り畳み画像に再構成される。その折り畳み画像は、次の行列式を解くことによるような適切な展開アルゴリズムを用いて展開され、
m(ky)=S(ry,ky)p(ry) (1)
ここで、kyはフーリエエンコードステップの指数であり、ryは展開画像位置の指数であり、m(ky)は折り畳み画像データを有し、p(ry)は展開画像に対応し、S(ry,ky)は、列に沿ったコイルの感度及び対応するフーリエエンコードステップを含む感度マトリクスに対応する。展開画像項pは展開画像のN個の真の画素値を有する。折り畳み画像項mはCN/R個の測定値を有し、ここで、Cはコイル数であり、Rはアンダーサンプリングに対応する時間短縮係数である。式(1)は展開画像を計算するようにp(ry)について解かれる。
m(ky)=S(ry,ky)p(ry) (1)
ここで、kyはフーリエエンコードステップの指数であり、ryは展開画像位置の指数であり、m(ky)は折り畳み画像データを有し、p(ry)は展開画像に対応し、S(ry,ky)は、列に沿ったコイルの感度及び対応するフーリエエンコードステップを含む感度マトリクスに対応する。展開画像項pは展開画像のN個の真の画素値を有する。折り畳み画像項mはCN/R個の測定値を有し、ここで、Cはコイル数であり、Rはアンダーサンプリングに対応する時間短縮係数である。式(1)は展開画像を計算するようにp(ry)について解かれる。
可変密度SENSE技術を用いて撮像データを迅速に取得するように、アンダーサンプリングが増加される。しかしながら、k空間のサンプリングがまばら過ぎる場合、式(1)は、過少決定又は特異な場合に繋がる。このことは計算上の困難性に繋がり、折り畳み画像においてアーティファクトを生成する可能性がある。このような問題点に対処するように、展開が、演繹的知識を用いて正規化されることができる。正規化は、演繹的知識の性質及び品質に依存する。典型的には、この付加情報は、展開中に可変密度SENSEデータと結合された付加磁気共鳴信号情報を取得することにより与えられる。
そのような正規化技術にも拘わらず、減少されたアーティファクトを有する展開画像を効果的に生成するようにアンダーサンプルされた可変密度SENSE撮像データの処理は難しい問題点のまま残されている。演繹的知識を有することにより得られる画像品質の改善は、演繹的知識を取得するために必要な長い撮像期間によりオフセットされる。更に、対象物の動き、又は演繹的知識のコンテンツの取得と可変みつドSENSEデータとの間に生じる他の変化が、演繹的知識と折り畳み画像との間の空間的位置ずれ又は他の不一致をもたらす可能性がある。
米国特許第6,380,741号明細書
Preessmann et al.,Magnetic Resonance in Medicine 42,pp.952−962(1999)
本発明は、上記の欠点及び他の欠点を克服する改善された装置及び方法について検討したものである。
一特徴に従って、磁気共鳴撮像方法を提供する。可変密度感度エンコードデータが、k空間の中心及びその中心の近傍において高密度で、そして、k空間から離れたところでは低いアンダーサンプルされた密度で取得される。1つ又はそれ以上の正規化画像は、k空間の中心及びその中心の近傍において可変密度感度エンコードデータの高密度部分から構成される。可変密度感度エンコードデータは、展開された再構成画像に再構成される。その再構成は、高密度及び低密度の可変密度感度エンコードデータを複数の折り畳み画像に再構成することと、1つ又はそれ以上の正規化画像を用いて展開画像を生成するように折り畳み画像を展開することとを有する。
他の特徴に従って、磁気共鳴撮像装置について開示している。複数の高周波コイルは、k空間の中心から少なくとも離れたところでアンダーサンプルされた可変密度感度エンコードデータを取得する。各々のコイルのための再構成処理器は、コイルにより取得されたk空間の中心において又はその中心近傍において配置されている可変密度感度エンコードデータの高密度部分から再構成された正規化画像と、コイルにより取得された可変密度感度エンコードデータから再構成された折り畳み画像とを再構成する。展開処理器は折り畳み画像を展開する。その展開は正規化画像により正規化される。
一有利点は、改善された速度を有するアンダーアンプルされた可変密度SENSE撮像データの再構成にある。
他の有利点は、最小アーティファクトを有するアンダーアンプルされた可変密度SENSE撮像データの再構成にある。
他の有利点は、例えば、径方向又は螺旋サンプリングのような種々のk空間サンプリング軌道により取得されるサンダーサンプルされた可変密度SENSE撮像データを処理する能力にある。
多くの付加的有利点及び優位性については、以下の好適な実施形態についての詳細説明に従って読むときに、当業者に明らかに成るであろう。
本発明は、種々の構成要素及び構成要素の構成、種々の処理操作及び処理操作の構成において実施される。図面は、単に好適な実施形態を例示することを目的とするものであり、本発明を制限するように意図されるものではない。
図1を参照するに、磁気共鳴撮像スキャナ10は、一般に円筒状のスキャナボア14を規定するハウジング12を有し、そのスキャナボアの内側に、関連撮像対象物16が配置されている。主磁場コイル20はハウジング12の内側に備えられている。主磁場コイル20は、スキャナボア14の中心軸22に沿って方向付けられている主B0磁場を生成するように一般電磁構成で備えられている。主磁場コイル20は、典型的には、クライオシュラウディング24の内側に備えられた超伝導コイルであるが、抵抗性主磁石をまた、用いることが可能である。更に、スキャナ10は、撮像対象物16にアクセスするための円筒形スキャナボア14の端部以外の付加アクセス開口を有することが可能である。例えば、閉じた一般に“O字”型断面を有する閉じた電磁構成ではなく、更に開いた一般に“U字”型断面の磁石を用いることが可能である。
ハウジング12はまた、中心軸22に対して横断する面内方向に沿って又は他の選択された方向に沿って、ボア14の中心軸22に平行な磁場勾配を選択的に生成するために磁場勾配コイル30を収容する又は支持する。任意に、ハウジング12は、磁気共鳴を選択的に励起する及び/又は検出するためのバードケージ高周波全身コイル32を更に収容する又は支持する。ボア14の内側に備えられている感度エンコード(SENSE)コイルアレイ34は、複数のSENSEコイル、特に、実施例としてのコイルアレイ34における4つのSENSEコイルを有するが、他の数のコイルを用いることが可能である。コイルアレイ34のコイルは、磁気共鳴信号に対して異なる感度を有し、それ故、SENSEエンコードを可能にする。一実施形態においては、コイルアレイ34は、撮像対象物16の近傍に備えられている表面コイルのアレイである。コイルアレイではなく、例えば、高周波コイルの複数の空間的に別個の部分から受信される高周波信号を移植することにより規定されるアンテナのような、適切な空間感度プロファイルを有する他の高周波アンテナをSENSEデータを取得するために用いることができる。
SENSEコイルアレイ34は、任意のバードケージ又は他の全身コイル32により励起される磁気共鳴を受信するために用いられる、又は、磁気共鳴はSENSEコイルアレイ34はコイルアレイ34により励起及び受信される。ハウジング12は、典型的には、スキャナボア14を規定する整形内側ライナー36を有する。
主磁場コイル20は主磁場B0を生成する。磁気共鳴撮像制御器40は、磁場勾配コイル30を選択的に電源供給するように磁石制御器42を動作させ、高周波コイル又はコイルアレイ32、34を選択的に電力供給するように、図示している高周波コイル32に結合された又はSENSEコイルアレイ34に結合された高周波送信器44を動作させる。磁場勾配コイル30及び高周波コイル又はコイルアレイ32、34を選択的に動作させっることにより、磁気共鳴は生成され、撮像対象物16の注目すべき選択された領域の少なくとも一部において空間的にエンコードされる。勾配コイル30を介して選択された磁場勾配を適用することにより、直交軌道、複数の径方向軌道又は螺旋軌道等の選択されたk空間軌道は横断される。選択されたk空間軌道の横断中、磁気共鳴撮像制御器40は、SENSEコイルアレイ34のコイルを用いてk空間サンプルを読み出すようにSENSEコイルアレイ34に結合されている高周波送信器46を操作する。k空間サンプルはk空間メモリ50に記憶されている。
記憶されているk空間サンプルは、各々のサンプルを受信したアレイの受信コイルによりソーティングされる。k空間サンプリングは、k空間の中心における又はその中心近傍における高密度k空間サンプリング及びk空間の中心から離れたところでアンダーサンプルされる低密度k空間サンプリングと共に可変である。再構成処理器52は、折り畳み画像メモリ又はバッファ60に記憶されている対応する折り畳み画像に各々のコイルにより取得されたk空間の全てをスパンするアンダーサンプルされたk空間データを再構成する。更に、高サンプリング密度で各々のコイルにより取得されたk空間の中心における又はその中心近傍におけるk空間の一部は、正規化画像メモリ又はバッファ62に記憶されている対応する低分解能正規化画像に再構成処理器52により再構成される。
展開処理器66は、再構成画像62により正規化される展開アルゴリズムを用いて折り畳み画像60を折り畳む。1つの適切な方法においては、折り畳み画像60は、次式のようなエラーパラメータδを最小化することにより展開され、
ここで、kyはフーリエエンコードステップの指数であり、ryは指数yにおける画像位置を表し、m(ky)は折り畳み画像に対応し、p(ry)は展開画像に対応し、S(ry,ky)は、列及び対応するフーリエエンコードステップに沿ったコイルの感度を含む感度マトリクスに対応し、Ireg(ry)は、正規化画像メモリ62に記憶されている低分解能画像から構成される正規化画像に対応する。展開画像項p破綻界画像のN個の真の画素値を有する。折り畳み画像項mはCN/R個の測定値を有し、ここで、Cはコイル数であり、Rはアンダーサンプルに対応する時間短縮係数である。
一実施形態においては、正規化画像Ireg(ry)は、種々の受信コイルにより取得された低分解能画像を組み合わせることにより構成される。低分解能正規化画像は、合成画像を生成するように二乗和処理を用いて結合される。他の実施形態においては、受信コイルの一により取得された単一の低解像度正規化画像が正規化画像として用いられる。この実施形態においては、選択画像は、好適には、撮像領域に亘って高い均一の感度を有する受信コイルにより取得される。
式(2)における演算子
はxのユークリッドノルムを表す。式(2)は、折り畳み画像に展開画像の忠実性を示す展開項
と、1つ又はそれ以上の正規化画像62に折り畳み画像の忠実性を示す正規化項
とを有する。パラメータqは、展開項及び正規化項の相対的影響を決定する正規化重み付けパラメータである。q=0である場合、非正規化の式(1)が回復されることが理解できる。
展開画像を得るように式(2)を評価するための一方法においては、展開処理器66は、展開画像項p(ry)に対して式(2)を最小化させる。この最小化は、好適には、分析的に、次式に従ってなされ、
p=Ireg+SH(SSH+q2)−1(m−SIreg) (3)
ここで、上付きの添え字“H”はSの転置複素共役であり、式(2)に示されている式(3)のky及びry依存性は省かれている。式(3)は、例えば、正規化重み付けパラメータqに対して式(3)の画像項pの最小二乗法の最小化を実行することにより、正規化重み付けパラメータqに対して更に最適化されることができる。
p=Ireg+SH(SSH+q2)−1(m−SIreg) (3)
ここで、上付きの添え字“H”はSの転置複素共役であり、式(2)に示されている式(3)のky及びry依存性は省かれている。式(3)は、例えば、正規化重み付けパラメータqに対して式(3)の画像項pの最小二乗法の最小化を実行することにより、正規化重み付けパラメータqに対して更に最適化されることができる。
展開画像は画像メモリ70に記憶され、ユーザインターフェース72に表示され、不揮発性メモリに記憶され、ローカルイントラネット又はインターネットにおいて送信され、視覚化され、記憶され、操作される、等なされる。ユーザインターフェース72はまた、磁気共鳴撮像スキャナ10の放射線技師、技術者又は他のオペレータが、磁気共鳴撮像シーケンスを選択し、修正し、実行するように磁気共鳴撮像制御器40と通信することを可能にする。
図2を参照するに、各々のコイルについての適切な可変密度SENSE k空間サンプリングプロファイル80を示している。サンプリング密度80は、一定であり且つアンダーサンプルされていない、k空間の中心に又はその中心の近傍に配置されている高密度部分82を有する。このアンダーサンプルされていない高密度のk空間部分82は、折り畳みを有しない再構成画像に再構成される。k空間部分82の制限されていないk空間領域のために、結果として得られる再構成画像62は撮像領域の低分解能画像である。
k空間サンプリングプロファイル80はまた、k空間の中心から離れたところに配置され、−kmax及びkmaxそれぞれまで広がっているアンダーサンプル部分84、86を有する。一実施形態においては、アンダーサンプルされていないk空間部分82は、可変密度SENSEサンプリングプロファイル80のk空間範囲[−kmax,kmax]の約1/8の範囲である。この1/8の比は良好な正規化画像を与える一方、尚も、データ取得速度における実質的な改善を与えている。しかしながら、SENSEの中心、アンダーサンプルされたk空間部分80及び全k空間範囲[−kmax,kmax]間の他の比がまた、適切である。部分82は少なくともアンダーサンプルされていず、冗長データを与えるように、任意にオーバーサンプルされる。更に、部分82のオーバーサンプリングが実行される場合、そのオーバーサンプリングは、部分82の一部又は全てに亘って非一定のサンプリング密度を任意に有する。
k空間サンプリングプロファイル80は、アンダーサンプルされていない中心部分のk空間部分82とアンダーサンプルされている外側部分84、86の各々との間で一般的にガウス分布のサンプリング密度の変化を有する。破線で図2に示している直線状遷移のような他のサンプリング密度の遷移がまた、可能である。
上記の実施形態においては、直交k空間サンプリング軌道が用いられ、その軌道において、位相エンコードラインが読み取られ、そして中心(k=0近傍)ラインは十分にサンプリングされ、且つ更に周方向の位相エンコードラインは、選択された可変密度SENSEアンダーサンプリングに従って間隔を置かれている。他の実施形態においては、二次元SENSEが用いられ、その二次元SENSEにおいて、第2の次元は、典型的には、x−y面に対して横断する。他の実施形態においては、直交k空間サンプリング軌道が用いられるが、可変密度SENSEアンダーサンプリングが、位相エンコード方向ではなく読み取り方向で実行される。この実施形態においては、式(2)のパラメータky及びryは、読み取り方向に対応するkx及びrxで適切に置き換えられる。他の直交k空間サンプリング軌道において、SENSEアンダーサンプリングが位相エンコード及び読み取り要項の両方で実行される。
図3を参照するに、他の実施形態において、径方向のサンプリング軌道90が用いられる。径方向のサンプリング軌道90は、k空間の中心に又はその中心近傍に収束する。このような径方向のサンプリング軌道90の幾何学的構成は、径方向のサンプリング軌道90の収束のために、k空間中心の中心近傍で高サンプリング密度領域92を与える。高サンプリング密度領域92は、好適には、各々のコイルについてオーバーサンプルされ、又は少なくともアンダーサンプルされず、各々のコイルにより取得される撮像データの高サンプリング密度領域92は、再構成画像62を生成するように適切に再構成される。図3は二次元的な径方向のサンプリング軌道90を示し、その図はまた、三次元的に半球又は完全球状の立体角に広がった径方向軌道を採用することがまた、可能である。
図4を参照するに、他の実施形態において、k空間の中心から遠ざかって又はその中心に向かって螺旋を描く螺旋サンプリング軌道96が採用される。図示している螺旋軌道96は、好適には、一定の間隔Pを有するが、k空間の中心から距離が離れるにつれて変化して広くなる又は狭くなる間隔をまた、用いることが可能である。その螺旋は二次元又は三次元であることが可能である。三次元においては、その螺旋は連続的に螺旋を描く、又は、同心円の表面に存在する等が可能である。また、螺旋は、一般的に、四角形状螺旋、六角形状螺旋、十二角形状螺旋等を規定する複数の直線状部分を有することが可能である。高サンプリング密度領域98は、各々のコイルについてオーバーサンプルされ又は少なくともアンダーサンプルされず、各々のコイルにより取得される撮像データの高サンプリング密度領域98は、正規化画像62を生成するように適切に再構成される。k空間の中心の高サンプリング密度を伴う他のサンプリングスキームをまた、検討することが可能である。
本発明について、好適な実施形態に関連して詳述した。明らかに、上記説明を読み且つ理解することにより、当業者は、修正及び変形を考案することが可能である。本発明は、同時提出の特許請求の範囲における又はそれと同等の範囲内にあるそのような修正及び変形全てを網羅することが意図されている。
Claims (20)
- 磁気共鳴撮像方法であって:
k空間の中心及びk空間の中心近傍において高密度の、及びk空間の中心から離れたところでアンダーサンプルされた低密度の可変密度感度エンコードデータを取得する段階; k空間の中心及びk空間の中心近傍に配置されている前記可変密度感度エンコードデータの高密度部分から1つ又はそれ以上の再構成画像を構成する段階;並びに
前記可変密度感度エンコードデータを展開された再構成画像に再構成する段階;
を有する磁気共鳴撮像方法であり、前記再構成する段階は:
前記高密度及び低密度の可変密度感度エンコードデータを複数の折り畳み画像に再構成する手順と;
前記1つ又はそれ以上の再構成画像を用いて前記展開画像を生成するように前記折り畳み画像を展開する手順と;
を有する、ことを特徴とする磁気共鳴撮像方法。 - 請求項1に記載の磁気共鳴撮像方法であって、前記複数の折り畳み画像は、複数のアンテナの対応するアンテナにより各々、取得され、前記の1つ又はそれ以上の正規化画像を構成する段階は:
各々のアンテナにより取得された前記可変密度感度エンコードデータの前記高密度部分から低分解能画像を再構成する手順と;
前記展開する手順において用いられる前記正規化画像を得るように前記の再構成された低分解能画像を結合する手順と;
を有する、ことを特徴とする磁気共鳴撮像方法。 - 請求項2に記載の磁気共鳴撮像方法であって、各々のアンテナにより取得された前記可変密度感度エンコードデータの前記高密度部分はアンダーサンプルされていない、ことを特徴とする磁気共鳴撮像方法。
- 請求項3に記載の磁気共鳴撮像方法であって、前記可変密度感度エンコードデータの前記高密度部分は、前記可変密度感度エンコードデータのk空間領域の約1/8の範囲である、ことを特徴とする磁気共鳴撮像方法。
- 請求項2に記載の磁気共鳴撮像方法であって、各々のアンテナにより取得された前記可変密度感度エンコードデータの前記高密度部分はオーバーサンプルされ且つ冗長データを有する、ことを特徴とする磁気共鳴撮像方法。
- 請求項1に記載の磁気共鳴撮像方法であって、感度エンコード画像を展開することは: 前記折り畳み画像に対する前記展開画像の忠実度を表す展開項及び前記正規化画像に対する前記展開画像の忠実度を表す正規化項の重み付けされた組み合わせを有するペナルティ関数を最適化すること;
を有する、ことを特徴とする磁気共鳴撮像方法。 - 請求項1に記載の磁気共鳴撮像方法であって、可変密度感度エンコードデータを取得する段階は:
アンダーサンプリングされていない均一のk空間サンプリング密度を有するk空間の中心及びその中心近傍に配置された前記高密度部分を取得する手順と;
前記高密度部分から距離が離れるにつれて徐々に減少するk空間サンプリング密度を用いてk空間の中心から離れたところでアンダーサンプルされたk空間データを取得する手順と;
を有する、ことを特徴とする磁気共鳴撮像方法。 - 請求項7に記載の磁気共鳴撮像方法であって、前記高密度部分と前記のk空間の中心から離れたところでアンダーサンプルされたk空間データとの間のサンプリング密度遷移は直線的且つガウス分布形状のサンプリング密度遷移を有する、ことを特徴とする磁気共鳴撮像方法。
- 請求項1に記載の磁気共鳴撮像方法であって、可変密度感度エンコードデータを取得する段階は:
非直交軌道を用いて可変密度感度エンコードデータを取得する手順であって、前記可変密度感度エンコードデータの前記高密度部分は前記非直交軌道の幾何学的構成により規定されている、手順;
を有する、ことを特徴とする磁気共鳴撮像方法。 - 請求項9に記載の磁気共鳴撮像方法であって、非直交軌道を用いて可変密度感度エンコードデータを取得する前記段階は:
複数の径方向k空間サンプリング軌道を取得する手順であって、前記可変密度感度エンコードデータの前記高密度部分は前記のk空間の中心及びその中心近傍において前記の複数の径方向k空間サンプリング軌道の収束により規定される、手順;
を有する、ことを特徴とする磁気共鳴撮像方法。 - 請求項9に記載の磁気共鳴撮像方法であって、非直交軌道を用いて可変密度感度エンコードデータを取得する前記段階は:
螺旋状k空間サンプリング軌道を取得する手順であって、前記可変密度感度エンコードデータの前記高密度部分は前記k空間サンプリング軌道の中心領域により規定される、手順;
を有する、ことを特徴とする磁気共鳴撮像方法。 - 請求項11に記載の磁気共鳴撮像方法であって、前記螺旋状k空間サンプリング軌道は、k空間の中心から距離が離れるにつれて増加する、広がる螺旋間隔と均一の螺旋間隔の一を有する、ことを特徴とする磁気共鳴撮像方法。
- 請求項9に記載の磁気共鳴撮像方法であって、前記の1つ又はそれ以上の正規化画像を構成する段階は:
前記可変密度感度エンコードデータの前記高密度部分からの前記取得において用いられる高周波アンテナに対応する複数の低分解能画像を再構成する手順であって、前記展開において用いられる前記正規化画像は前記再構成された低分解能画像から構成される、手順;
を有する、ことを特徴とする磁気共鳴撮像方法。 - 請求項9に記載の磁気共鳴撮像方法であって、前記可変密度感度エンコードデータの前記高密度部分はアンダーサンプルされていない、ことを特徴とする磁気共鳴撮像方法。
- 磁気共鳴撮像装置であって:
主磁石;
磁場勾配コイル;
複数の高周波受信コイル;及び
請求項1の前記磁気共鳴撮像方法を実行する処理器;
を有することを特徴とする磁気共鳴撮像装置。 - 磁気共鳴撮像装置であって:
前記のk空間の中心から少なくとも離れたところでアンダーサンプルされた可変密度感度エンコードデータを取得する複数の高周波コイル;
各々のコイルについて、そのコイルにより取得されたk空間の中心又はその中心近傍に配置された前記可変密度感度エンコードデータの高密度部分から再構成された正規化画像と、そのコイルにより取得された前記可変密度感度エンコードデータから再構成された折り畳み画像と、を再構成する再構成処理器;及び
前記折り畳み画像を展開する展開処理器であって、前記展開は前記正規化画像により正規化される、展開処理器;
を有することを特徴とする磁気共鳴撮像装置。 - 請求項16に記載の磁気共鳴撮像装置であって:
前記可変密度感度エンコードデータの取得中、非直交k空間サンプリング軌道を生成する磁場勾配コイルであって、前記非直交k空間サンプリング軌道は前記可変密度感度エンコードデータの前記高密度部分を規定する幾何学的構成を有する、磁場勾配コイル;
を更に有する、ことを特徴とする磁気共鳴撮像装置。 - 請求項16に記載の磁気共鳴撮像装置であって:
複数の径方向k空間サンプリング軌道及び螺旋状k空間サンプリング軌道の一を生成する磁場勾配コイル;
を、前記可変密度感度エンコードデータの前記取得中に更に有する、ことを特徴とする磁気共鳴撮像装置。 - 請求項16に記載の磁気共鳴撮像装置であって、前記展開処理器は:
前記折り畳み画像に対する前記展開画像の忠実度を表す展開項;及び
前記正規化画像に対する前記展開画像の忠実度を表す正規化項;
の重み付けされた組み合わせを有するペナルティ関数を最適化する、ことを特徴とする磁気共鳴撮像方法。 - 請求項16に記載の磁気共鳴撮像装置であって、前記複数の高周波コイルは:
アンダーサンプルされていない均一なk空間サンプリング密度を有する前記のk空間の中心又はその中心近傍に配置された前記高密度部分;及び
前記高密度部分から距離が離れるにつれて徐々に減少するk空間サンプリング密度を用いて前記のk空間の中心から離れたアンダーサンプルされたk空間データ;
を取得する、ことを特徴とする磁気共鳴撮像方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US53639904P | 2004-01-14 | 2004-01-14 | |
PCT/IB2005/050051 WO2005069031A1 (en) | 2004-01-14 | 2005-01-05 | Regularized variable density sense |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2007518485A true JP2007518485A (ja) | 2007-07-12 |
Family
ID=34794401
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2006548510A Withdrawn JP2007518485A (ja) | 2004-01-14 | 2005-01-05 | 正規化可変密度センス |
Country Status (5)
Country | Link |
---|---|
US (1) | US7403005B2 (ja) |
EP (1) | EP1706755A1 (ja) |
JP (1) | JP2007518485A (ja) |
CN (1) | CN1910469A (ja) |
WO (1) | WO2005069031A1 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017000573A (ja) * | 2015-06-12 | 2017-01-05 | 株式会社日立製作所 | 磁気共鳴イメージング装置およびその制御方法 |
JP2021514710A (ja) * | 2018-02-22 | 2021-06-17 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | マルチグラジエントエコーシーケンスを使用したディクソンmrイメージング |
Families Citing this family (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8538099B2 (en) * | 2005-03-23 | 2013-09-17 | General Electric Company | Method and system for controlling image reconstruction |
EP1902424B1 (en) * | 2005-07-08 | 2009-11-11 | Wisconsin Alumni Research Foundation | Highly constrained image reconstruction method |
WO2007008530A1 (en) | 2005-07-08 | 2007-01-18 | Wisconsin Alumni Research Foundation | Backprojection reconstruction method for ct imaging |
WO2007038227A1 (en) | 2005-09-22 | 2007-04-05 | Wiconsin Alumni Research Foundation | Functional mri involving a highly constrained backprojection |
US7408347B2 (en) | 2005-09-22 | 2008-08-05 | Wisconsin Alumni Research Foundation | Highly constrained magnetic resonance spectroscopy image reconstruction method |
EP1927009A1 (en) | 2005-09-22 | 2008-06-04 | Wisconsin Alumni Research Foundation | Reconstruction of images of the beating heart using a highly constrained backprojection |
WO2007038206A1 (en) | 2005-09-22 | 2007-04-05 | Wisconsin Alumni Research Foundation | Highly constrained backprojection reconstruction in diffusion weighted mri |
EP1927010B1 (en) | 2005-09-22 | 2016-01-06 | Wisconsin Alumni Research Foundation | Highly constrained reconstruction of motion encoded mr images |
US7592809B1 (en) * | 2005-11-23 | 2009-09-22 | General Electric Company | Hybrid k-t method of dynamic imaging with improved spatiotemporal resolution |
US7609058B2 (en) * | 2006-11-17 | 2009-10-27 | Siemens Medical Solutions Usa, Inc. | Method and apparatus for generating a magnetic resonance data file |
CN101201392B (zh) * | 2006-11-17 | 2011-05-04 | 美国西门子医疗解决公司 | 产生磁共振数据文件的方法和设备 |
EP1959397B1 (en) | 2007-02-19 | 2019-08-07 | Wisconsin Alumni Research Foundation | Iterative HYPR medical image reconstruction |
EP1959396B1 (en) | 2007-02-19 | 2012-01-18 | Wisconsin Alumni Research Foundation | Localized and highly constrained image reconstruction method |
US8306289B1 (en) * | 2007-02-23 | 2012-11-06 | University Of Virginia Patent Foundation | Method and system for off-resonance correction for non-cartesian parallel image reconstruction |
DE102007009209B3 (de) * | 2007-02-26 | 2008-08-28 | Siemens Ag | Verfahren zur Erzeugung von zeitlich hochaufgelösten MR-Aufnahmen durch Multiplikation unterabgetasteter Bilder mit einem Mittelwertbild |
WO2008111416A1 (ja) * | 2007-03-09 | 2008-09-18 | Hitachi Medical Corporation | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 |
WO2008155703A1 (en) * | 2007-06-19 | 2008-12-24 | Koninklijke Philips Electronics N.V. | Mri radio frequency receiver comprising digital down converter |
ES2362776B1 (es) | 2009-11-24 | 2012-05-22 | Telefonica,S.A. | Procedimiento de comunicación de est�?mulos f�?sicos mediante dispositivos móviles. |
EP2503349A1 (en) | 2011-03-24 | 2012-09-26 | Koninklijke Philips Electronics N.V. | MR image reconstruction using prior information-constrained regularization |
GB2504642B (en) * | 2011-04-21 | 2018-02-21 | Koninl Philips Electronics Nv | Magnetic resonance imaging of object inmotion |
US8620404B2 (en) | 2011-07-26 | 2013-12-31 | Wisconsin Alumni Research Foundation | System and method of high-frame rate, time-resolved, three-dimensional magnetic resonance angiograpy |
CN102565737A (zh) * | 2011-12-12 | 2012-07-11 | 中国科学院深圳先进技术研究院 | 快速磁共振成像方法及系统 |
WO2013109965A1 (en) * | 2012-01-19 | 2013-07-25 | Brigham And Women's Hospital, Inc. | Data reconstruction for improved ultrasound imaging |
KR101630762B1 (ko) * | 2014-07-30 | 2016-06-15 | 삼성전자주식회사 | 자기 공명 영상 생성 장치 및 방법 |
JP6357594B1 (ja) * | 2015-06-08 | 2018-07-11 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 可変密度サンプリング方式mri |
JP6618786B2 (ja) * | 2015-11-17 | 2019-12-11 | キヤノンメディカルシステムズ株式会社 | 磁気共鳴イメージング装置及び画像処理装置 |
DE102016201458B4 (de) * | 2016-02-01 | 2018-10-31 | Siemens Healthcare Gmbh | Radiales Erfassen von MR-Daten zur Vermeidung von Artefakten |
CN107145885B (zh) * | 2017-05-03 | 2020-12-15 | 金蝶软件(中国)有限公司 | 一种基于卷积神经网络的单字图文字识别方法及装置 |
CN109143132B (zh) * | 2018-09-12 | 2020-11-27 | 上海联影医疗科技股份有限公司 | 磁共振信号采集方法、成像方法、系统和介质 |
US11486954B2 (en) * | 2020-08-24 | 2022-11-01 | Siemens Healthcare Gmbh | Systems and methods for extending reconstructions to non-uniform k-space sampling |
JP2022185902A (ja) * | 2021-06-03 | 2022-12-15 | 富士フイルムヘルスケア株式会社 | 磁気共鳴イメージング装置、及び、画像処理装置 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6487435B2 (en) | 1998-04-10 | 2002-11-26 | Wisconsin Alumni Research Foundation | Magnetic resonance angiography using undersampled 3D projection imaging |
US6353752B1 (en) | 1999-05-14 | 2002-03-05 | Board Of Trustees Of The Leland Standford Junior University | Reduced field-of-view method for cine magnetic resonance imaging |
GB9926923D0 (en) | 1999-11-15 | 2000-01-12 | Marconi Electronic Syst Ltd | Magnetic resonance imaging |
WO2001073463A1 (en) * | 2000-03-24 | 2001-10-04 | Koninklijke Philips Electronics N.V. | Magnetic resonance imaging method with sub-sampling |
JP5085832B2 (ja) | 2000-07-31 | 2012-11-28 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | サブサンプリングされた捕捉を用いた磁気共鳴撮像方法 |
EP1307758A1 (en) * | 2000-07-31 | 2003-05-07 | Koninklijke Philips Electronics N.V. | Magnetic resonance method for forming a fast dynamic image |
DE10059772A1 (de) | 2000-11-30 | 2002-06-13 | Philips Corp Intellectual Pty | MR-Bildrekonstruktion |
US6556009B2 (en) * | 2000-12-11 | 2003-04-29 | The United States Of America As Represented By The Department Of Health And Human Services | Accelerated magnetic resonance imaging using a parallel spatial filter |
DE10119660B4 (de) | 2001-04-20 | 2006-01-05 | Siemens Ag | Verfahren zur schnellen Gewinnung eines Magnetresonanzbildes |
US6559642B2 (en) | 2001-05-09 | 2003-05-06 | Ge Medical Systems Global Technology Company, Llc | Calibration method for use with sensitivity encoding MRI acquisition |
AU2003234310A1 (en) | 2002-05-01 | 2003-11-17 | The Brigham And Women's Hospital, Inc. | Variable k-space-density parallel magnetic resonance imaging |
-
2005
- 2005-01-05 EP EP05702581A patent/EP1706755A1/en not_active Ceased
- 2005-01-05 CN CNA2005800023854A patent/CN1910469A/zh active Pending
- 2005-01-05 US US10/597,060 patent/US7403005B2/en not_active Expired - Fee Related
- 2005-01-05 WO PCT/IB2005/050051 patent/WO2005069031A1/en not_active Application Discontinuation
- 2005-01-05 JP JP2006548510A patent/JP2007518485A/ja not_active Withdrawn
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017000573A (ja) * | 2015-06-12 | 2017-01-05 | 株式会社日立製作所 | 磁気共鳴イメージング装置およびその制御方法 |
JP2021514710A (ja) * | 2018-02-22 | 2021-06-17 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | マルチグラジエントエコーシーケンスを使用したディクソンmrイメージング |
Also Published As
Publication number | Publication date |
---|---|
US7403005B2 (en) | 2008-07-22 |
CN1910469A (zh) | 2007-02-07 |
EP1706755A1 (en) | 2006-10-04 |
WO2005069031A1 (en) | 2005-07-28 |
US20070182412A1 (en) | 2007-08-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2007518485A (ja) | 正規化可変密度センス | |
CN107750338B (zh) | 使用星形叠层采集的mr成像方法和设备 | |
JP6232043B2 (ja) | K空間領域の中央および周辺のサンプリング最中での呼吸運動の検出を用いた磁気共鳴画像の再構成 | |
EP1716430B1 (en) | Magnetic resonance imaging method | |
JP5926285B2 (ja) | 任意k空間トラジェクトリの高速並行再構成 | |
JP4657710B2 (ja) | 前情報により向上される動的な磁気共鳴イメージング | |
US8379951B2 (en) | Auto calibration parallel imaging reconstruction method from arbitrary k-space sampling | |
JP4443625B2 (ja) | 磁気共鳴による撮像方法及び装置 | |
US7439739B2 (en) | Anti-aliased magnetic resonance image reconstruction using partially parallel encoded data | |
US9964615B2 (en) | MR image reconstruction using compressed sensing | |
CN106471389B (zh) | 在多激发mri中由于激发间运动造成的伪影的降低 | |
CN104939828B (zh) | 磁共振成像方法 | |
US20070096732A1 (en) | Parallel magnetic resonance imaging method using a radial acquisition trajectory | |
JP6445470B2 (ja) | パラレルのマルチスライスmr画像法 | |
JP2014508622A (ja) | 事前情報に制約される正則化を使用するmr画像再構成 | |
US20240012078A1 (en) | Mri in an inhomogeneous field with no pulsed gradients | |
CN105022009B (zh) | 采集二维体积片段的磁共振数据的方法以及磁共振设备 | |
JP2008529643A (ja) | 短いエコー時間での磁気共鳴イメージング | |
JP2007508092A (ja) | 磁気共鳴画像化システム及び方法 | |
Lin et al. | GRAPPA operator for wider radial bands (GROWL) with optimally regularized self‐calibration | |
JP2006130285A (ja) | 磁気共鳴イメージング装置 | |
JP2006507071A (ja) | 磁気共鳴方法 | |
US7106062B1 (en) | Method and system for reducing coupling in magnetic resonance RF coil arrays | |
JP2005152175A (ja) | 磁気共鳴撮像装置および方法 | |
CN113567901B (zh) | 一种磁共振旋转坐标系下的自旋晶格弛豫成像方法和系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20080104 |
|
A761 | Written withdrawal of application |
Free format text: JAPANESE INTERMEDIATE CODE: A761 Effective date: 20090501 |