JP4047553B2 - 磁気共鳴イメージング装置 - Google Patents
磁気共鳴イメージング装置 Download PDFInfo
- Publication number
- JP4047553B2 JP4047553B2 JP2001122213A JP2001122213A JP4047553B2 JP 4047553 B2 JP4047553 B2 JP 4047553B2 JP 2001122213 A JP2001122213 A JP 2001122213A JP 2001122213 A JP2001122213 A JP 2001122213A JP 4047553 B2 JP4047553 B2 JP 4047553B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- sensitivity
- coil
- magnetic resonance
- calculation
- 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
Links
- 238000002595 magnetic resonance imaging Methods 0.000 title claims description 18
- 230000035945 sensitivity Effects 0.000 claims description 70
- 238000003384 imaging method Methods 0.000 claims description 46
- 238000009826 distribution Methods 0.000 claims description 45
- 238000004364 calculation method Methods 0.000 claims description 32
- 239000011159 matrix material Substances 0.000 claims description 31
- 238000005259 measurement Methods 0.000 claims description 13
- 238000005481 NMR spectroscopy Methods 0.000 claims description 8
- 230000002194 synthesizing effect Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000000034 method Methods 0.000 description 20
- 230000008569 process Effects 0.000 description 14
- 238000010586 diagram Methods 0.000 description 8
- 230000000877 morphologic effect Effects 0.000 description 8
- 238000001514 detection method Methods 0.000 description 5
- 239000000523 sample Substances 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 210000001015 abdomen Anatomy 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000013421 nuclear magnetic resonance imaging Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000004904 shortening Methods 0.000 description 2
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 1
- 239000012491 analyte Substances 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 239000013256 coordination polymer Substances 0.000 description 1
- 210000004351 coronary vessel Anatomy 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000005281 excited state Effects 0.000 description 1
- 238000002599 functional magnetic resonance imaging Methods 0.000 description 1
- 235000015220 hamburgers Nutrition 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000002610 neuroimaging Methods 0.000 description 1
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 1
- 229910052698 phosphorus Inorganic materials 0.000 description 1
- 239000011574 phosphorus Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
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/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)
- 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)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、被検体中の水素や燐等からの核磁気共鳴（以下、「NMR」という）信号を測定し、核の密度分布や緩和時間分布等を映像化する核磁気共鳴撮影（MRI）装置に関し、心臓等の高速撮影が要求される領域において、特に複数のRF受信コイルを用い、各RF受信コイルにて位相エンコードを間引いて取得した信号をRF受信コイルの感度分布を用いて行列演算により展開する撮影方法（以下、パラレルイメージング法）において、アーチファクトの無い画像を得ることを可能としたMRI装置に関する。
【０００２】
【従来の技術】
MRIでは、位相エンコード量を変えながらシーケンスを繰り返し実行し、１枚の画像再構成に必要なエコー信号を取得する。そのため、画像の取得時間は繰り返し回数が大きく影響する。高速撮影を行う場合、一般的には、１回の繰り返し内に複数のエコー信号を発生させるマルチエコータイプのシーケンスを用いたり、繰り返しの時間間隔を数〜数十msにまで短縮したシーケンスを用いる。しかし、マルチエコータイプのシーケンスでは、画像のコントラストが低下したり、画像歪みの原因となることがある。これは、画像のコントラストに寄与するエコー時間がマルチエコータイプのシーケンスでは各エコー信号で異なるため、コントラストの低下を招く。また、エコー時間が異なるため、各エコー信号間での位相変化が異なると、画像歪みとなって画像に現れる。
【０００３】
また、心臓領域（冠状動脈撮影など）を撮影する場合、より高速に画像を取得する必要があり、パラレルイメージング法と呼ばれる高速撮影方法も提案されている。パラレルイメージング法は、マルチプルRF受信コイルを用いて、位相エンコードを等間隔に間引いたシーケンスを実行し、繰り返し回数を低減して撮影時間を短縮する。通常、位相エンコードを等間隔に間引いて計測を行うと、画像には折り返しが発生するが、各RF受信コイルの感度分布をもとに行列演算を行うことで画像を展開し、折り返しを除去する。一般的にパラレルイメージング法では、撮影に用いたRF受信コイルの数分だけ撮影時間を短縮できる。
【０００４】
【発明が解決しようとする課題】
パラレルイメージング法では、マルチエコータイプのシーケンスを用いなくとも、撮影時間を短縮することができるため、コントラストの低下や画像歪みの影響を減らすことができる。しかしながら、マルチプルRF受信コイルの配置や感度分布の形状により結果の画像が大きく変化する。特に、背景等の低信号領域を含んだまま行列演算を行うと、ノイズの影響で折り返しの展開時に誤差が大きくなり、画像中に輝点のアーチファクトが発生することがある。また、各RF受信コイルの感度分布を精度良く算出するために、比較的感度分布の均一な全身用ボディコイルで計測した画像を用いて感度分布を算出するのが望ましいが、装置によっては全身用のボディコイルが無かったり、受信チャンネル数が少ないために同時に画像を取得できない問題が有る。
【０００５】
そこで、本発明の目的は、行列演算時の計算誤差を抑え、輝点アーチファクトの発生を無くし、また全身用ボディコイルの無い装置、或いはチャンネル数の少ない装置でもパラレルイメージングを可能とすることにある。あるいは、全身用ボディコイルを有する装置でも、装置構成及び信号処理フローを簡潔にすることにある。
【０００６】
【課題を解決するための手段】
上記目的を達成するために、本発明は、複数のＲＦ受信コイルを備えて被検体からの核磁気共鳴信号を受信する受信手段と、前記受信手段を用いてｋ空間のエンコードステップを間引いて前記核磁気共鳴信号を計測して、前記ＲＦ受信コイル毎のｋ空間データを取得する計測制御手段と、前記ＲＦ受信コイル毎の感度画像を用いて該ＲＦ受信コイル毎の感度分布データを取得する感度分布演算と、前記ＲＦ受信コイル毎のｋ空間データと感度分布データとを用いてパラレルイメージング法に基づく演算を適用して前記被検体の画像を再構成する画像再構成演算と、を行う演算手段と、を備えた磁気共鳴イメージング装置において、前記感度分布演算は、前記複数のＲＦ受信コイルの内の少なくとも２つの感度画像を合成して略均一感度画像を取得し、該略均一感度画像を用いて、前記複数のＲＦ受信コイルの内の少なくとも１つのＲＦ受信コイルの感度分布データを取得することを特徴とする。
【０００７】
また、少なくとも2個のRF受信コイルを備え、前記RF受信コイルで受信する核磁気共鳴信号の計測を計測空間のエンコードステップを間引くよう計測して各RF受信コイル毎に感度画像及び形態画像を取得し、前記RF受信コイルの感度画像に基づく感度分布から各形態画像の折り返し除去の演算を行なうと共に前記各形態画像を結合して１枚の形態画像を得る制御手段を有する磁気共鳴イメージング装置において、前記制御手段は、各RF受信コイルで取得した感度画像を結合した全体感度画像を作成し、前記全体感度画像を用いて各RF受信コイルの感度分布を算出し、前記全体感度画像を用いて画像の非画像領域と画像領域を分けるマスクを作成し、前記マスクに基づいて演算を行なうものである。
【０００８】
さらに、計測領域全体を受信する第1のRF受信コイルと、前記計測領域を少なくとも2つ以上分割した領域を受信する複数の第2のRF受信コイルを備え、前記各第2のRF受信コイルで受信する核磁気共鳴信号の計測を計測空間のエンコードステップを間引くよう計測して前記第2のRF受信コイル毎に感度画像及び形態画像を取得し、前記第1のRF受信コイルによる感度画像及び前記第2のRF受信コイルによる感度画像から求めた感度分布から前記形態画像の折り返し除去の演算を行なうと共に前記各形態画像を結合して１枚の形態画像を得る制御手段を有する磁気共鳴イメージング装置において、前記制御手段は、前記受信感度分布において非画像領域と画像領域を分けたマスクを作成し、前記マスクを用いて演算を行なうものである。
【０００９】
【発明の実施の形態】
以下、本発明の核磁気共鳴イメージング装置について、図面を参照して詳述する。図4は典型的な核磁気共鳴イメージング装置の構成である。被検体401の周囲に静磁場を発生する磁石402と、該空間に傾斜磁場を発生する傾斜磁場コイル403と、この領域に高周波磁場を発生するRFコイル404と被検体401が発生するMR信号を検出するRFプローブ405がある。傾斜磁場コイル403は、X、Y、Zの３方向の傾斜磁場コイルで構成され、傾斜磁場電源409からの信号に応じてそれぞれ傾斜磁場を発生する。RFコイル404はRF送信部410の信号に応じて高周波磁場を発生する。RFプローブ405の信号は、信号検出部406で検出され、信号処理部407で信号処理され、また計算により画像信号に変換される。画像は表示部408で表示される。傾斜磁場電源409、RF送信部410、信号検出部406は制御部411で制御され、制御のタイムチャートは一般にパルスシーケンスと呼ばれている。ベッド412は被検体が横たわるためのものである。
【００１０】
現在MRIの撮影対象は、臨床で普及しているものとしては、被検体の主たる構成物質、プロトンである。プロトン密度の空間分布や、励起状態の緩和現象の空間分布を画像化することで、人体頭部、腹部、四肢等の形態または、機能を2次元もしくは3次元的に撮影する。
【００１１】
RFコイル405の例として、複数の受信コイルを用いた「マルチプルRFコイル」もしくは「フェーズドアレイコイル」と呼ばれる技術が用いられている。マルチプルRFコイルとは、相対的に高感度な小型RF受信コイルを複数個並べて、各コイルで取得した信号を合成することより、RF受信コイルの高い感度を保ったまま視野を拡大し、高感度化を図る受信専用RFコイルである。水平磁場頭部用マルチプルRFコイルとしては、Array Head Coil for Improved Functional MRI(Christoph Leussler), 1996 ISMRM abstruct p.249がある。また、水平磁場頭部用QDマルチプルRFコイルとしては、Helmet and Cylindrical Shaped CP Array Coils for Brain Imaging: A Comparison of Signal-to-Noise Characteristics(H.A.Stark, E.M.Haacke), 1996 ISMRM abstract p.1412がある。水平磁場腹部用QDマルチプルRFコイルとしては、Four Channel Wrap-Around Coil with Inductive Decoupler for 1.5T Body Imaging(T.Takahashi et.al), 1995 ISMRM abstruct p.1418がある。
【００１２】
マルチプルRFコイルの信号検出部の一部を図3に示す。図3では、４個のRF受信コイル405が、それぞれプリアンプ302に接続されて一つのマルチプルコイル301を構成する。信号検出部406は、4個のAD変換・直交検波回路303が並列してなり、前記各プリアンプの出力が接続されている。信号処理部407は、各信号をフーリエ変換、バックプロジェクション法、ウエーブレット変換などにより、それぞれのRF受信コイル405で検出したMRI画像を求め、これらの画像信号を合成する演算304からなる。
【００１３】
次に、撮影方法を説明する。図6は一般的なグラディエントエコーシーケンスである。601は高周波パルス、602はスライス選択傾斜磁場パルス、603は位相エンコード傾斜磁場パルス、604は読み出し傾斜磁場パルス、605はサンプリングウインド、606はエコー信号、607は繰り返し時間（高周波パルス601の照射間隔）である。MRIでは、各繰り返し時間607毎に位相エンコード傾斜磁場パルス603の量を変え、異なる位相エンコードを与え、それぞれの位相エンコードで得られるエコー信号606を検出する。この操作を位相エンコードの数だけ繰り返し、画像取得時間608で1枚の画像再構成に必要なエコー信号を取得する。位相エンコードの数は通常1枚の画像あたり64、128、256、512等の値が選ばれる。各エコー信号は通常128、256、512、1024個のサンプリングデータからなる時系列信号として得られる。これらのデータを2次元フーリエ変換して1枚のMR画像を作成する。
【００１４】
マルチプルRF受信コイルを用いた高速撮影であるパラレルイメージング法の場合、位相エンコードステップ間隔を一定の割合で間引いて撮影の繰り返し回数を減らす。この間引き率を一般的に倍速数と呼ぶ。例えば、位相エンコードステップを2倍に間引く場合、倍速数は2である。以下、パラレルイメージング法の原理を、図2を用いて説明する。通常撮影の場合、図2(a)に示すように各位相エンコード量で取得した信号202nを配置し画像１枚分のデータ201を形成し、これをフーリエ変換することで図2(c)示すような画像を得る。次に、パラレルイメージング法において、例えば位相エンコードステップ間隔を2倍にし、データを間引いた場合、図2(b)に示すように、1ラインおきにデータ204mを計測し、205mの位置に相当するデータは計測しない。このとき、計測したデータ204mの量は通常撮影と比べ半分になるので、マトリクスを半分にして画像を作成するが、図2(d)のような折り返しの発生した画像が得られる。図2においては、y方向を位相エンコード方向としており、この折り返しは、位相エンコード方向の画像が折り返して発生する。つまり、図2(c)に示すような画像207内の被検体206の、上側領域2071内の被検体画像2061と下側領域2072内の被検体画像2062が重なり、図2(d)に示すような折り返しの有る画像208となる。このようにして発生した信号折り返しは、SENSEと呼ばれる信号処理法で除去する（SENSE:Sensitivity Encoding for Fast MRI(Klaas P.Pruessmann et.al), Magnetic Resonance in Medicine 42: 952-962 (1999)参照）。
【００１５】
以下、折り返し除去について説明する。x,y方向の画像マトリクスをそれぞれX,Yとした場合、画像内の座標(x,y)(x:1≦x≦X, y:1≦y≦Y)における画素値をsi (x,y)とする(ここで、iはコイル番号であり、2≦i≦Nである)。図2(d)の場合では、位相エンコードステップを2倍に間引いているので、間引き後の画像の位相エンコード方向のマトリクスは、Y'≡Y/2となる。図2(d)の画像の座標を(x,y')(1≦y'≦Y')としたとき、その画素値s'(x,y')208は、元画像207の2領域2071、2072が重なって、
【００１６】
【数１】
となる。ここでａは定数を表す。次に、RF受信コイルの感度分布と画像について述べる。i番目のRF受信コイルの2次元的な感度分布を、ci(x,y)とすると、受信した信号si (x,y)は、受信コイルの感度分布ci(x,y)と被検体のプロトン密度分布p(x,y)との積、
【００１７】
【数２】
で表される。(2)式を用いると、(1)式は、
【００１８】
【数３】
【００１９】
【数４】
とおくと、(3)式は、
【００２０】
【数５】
となる。これは、N行2列の行列として表すことができ、
【００２１】
【数６】
となる。これより、コイルの感度分布Cijが分かれば、逆行列を計算することで、被検体のプロトン密度分布Pjが分かる。
【００２２】
同様に、N個のコイルを用いて、M倍速で撮影を行った場合は、Y'≡Y/M、1≦y'≦Y'として、
【００２３】
【数７】
となる。ここでｂは定数を表す。逆行列演算によって折り返しを除去する事から、パラレルイメージング法のコイル数と倍速数との関係は、数学的にN≧Mである。
【００２４】
通常、パラレルイメージング法では、各RF受信コイルの感度分布Cijを前計測等で予め計測するなどして取得する。しかしながら、感度分布Cijを直接算出する事は難しく、一般的には比較的感度分布の均一な全身用ボディコイルの画像を用いて、各RF受信コイルで取得した画像をそれぞれ除算し、近似的なコイルの感度分布を求めて行列演算を行う。
【００２５】
図7に一般的なパラレルイメージング法の処理を示す。図では1つの全身用ボディコイルと、3つのマルチプルRF受信コイルの構成である。まず、全身用ボディコイルで取得した感度画像701と、各RF受信コイルで取得した感度画像7021〜7023を用いて、感度分布算出処理7041〜7043を行い、各RF受信コイルの感度分布7051〜7053を得る。感度分布算出処理704としては、全身用ボディコイルの画像をsc(x,y)として、例えば、
【００２６】
【数８】
がある。このようにして算出した感度分布7051〜7053と、各RF受信コイルで取得した折り返しの有る本計測画像7031〜7033を用いて、行列作成処理706により行列を作成後、逆行列演算処理707によって折り返しを除去した画像708を得る。
次に、本発明の第1の実施形態を図1，図5により説明する。従来の一般的なパラレルイメージング法では、背景等の低信号領域に関しては何ら考慮されていないものであるが、本実施形態では低信号領域を考慮した構成となっている。つまり、全身用ボディコイルの画像を用いて、マスク作成のステップ101により、マスク102を作成し、マスクの処理を行うステップ103により精度を向上した処理後のマスク104を作成し、マスクを用いた行列作成処理105を行うよう構成した。ここで、マスク処理について説明する。図5では3倍速の場合であり、画像501は通常画像(図5(b))の3領域5041、5042、5043が重なった結果である（ここで、Y'≡Y/3、1≦y'≦Y'である）。図5(a)中に、2箇所の注目点A(xA, y'A) 502、B(xB, y'B) 503を設ける。このとき、A、B点の画素値は、それぞれ、図の5021〜5023、5031〜5033の領域が重なっている事から、
【００２７】
【数９】
となる。ここでｄは定数を表す。しかし、5021は、被検体の無い背景部分の信号であり、これを用いて行列演算を行うと、逆行列演算時に行列が発散して、結果画像に輝点のアーチファクトを生じる場合が有った。そこで、このような背景の影響を無くすため、図5(c)に示すようなマスクm(x,y)102を用いる。図1中のマスク作成処理101としては、全身用ボディコイルの感度画像701に対して、例えば閾値を設けて図5（c）に示すような画素値が閾値以下の画素を背景領域506、画素値が閾値以上の画素を被検体領域505とする。また、閾値設定方法として、例えば画像内の最大画素値の1/10程度を閾値としたり、図5(d)に示すようなような画像内の画素値ヒストグラム507を作成し、ノイズの分布を求めて閾値508を設定する事もできる。次に、マスク作成101により、例えば被検体領域505の値を1に、背景領域506の値を0にしたマスク102を得る。そして、得られたマスク102に対して、さらにマスクの精度を向上させるため、マスクの処理103を行なう。そして、マスク処理103によって得られた処理後のマスク104を用いて行列作成処理105を行なう。これにより背景領域を除いて行列演算できるため、輝点アーチファクトの発生をなくすことができる。
【００２８】
ここで、マスク処理103としては、例えば被検体内部に構造が存在する場合はその領域の信号をそのまま使用するためマスク領域を内挿により埋めたり、本来計算に入れない領域に誤って取られた孤立点を除去する処理等を行なう。
また、マスクを用いた行列作成処理105では、行列の要素に背景領域を含まないように、各要素を0にする判別処理を行なう。例えば、
【００２９】
【数１０】
で表される。また、背景領域を0とした場合、画像が不自然となるような場合は、行列の要素を所定の定数にして処理してもよい。例えば
【００３０】
【数１１】
で表される。ここで、Const.は任意の定数を示す。
さらに、行列作成処理105内で(10)式に示すような各要素を0にする判別処理を用いることは、感度分布とマスクとの積、
【００３１】
【数１２】
とすることであり、判別処理を無くし処理を簡略化できる。
(10)式のように行列要素の一部を"0"とすることは、行列のマトリクスサイズを縮小することと等価である。すなわち、「各画素毎に折り返しの状況を調べ、必要十分な行列サイズを決定し、必要十分な演算を行う」ことになる。この結果、従来と比べ、不要なノイズの混入が無くなり、画質が向上する。又、演算時間が短縮するメリットも有る。
【００３２】
なお、行列作成処理105では処理後のマスク104を用いて演算したが、マスク102を用いて演算することもできる。これにより処理フローを簡略化することができる。また、マスク102は背景領域と被検体領域に分けて説明したが、背景領域とは被検体領域外の背景のみならず被検体領域中に存在する低信号領域（例えば、空隙等）も含む。
【００３３】
次に、本発明の第2の実施形態を図8を用いて説明する。この場合では、全身用ボディコイルの画像がない場合である。この場合、各RF受信コイルの感度分布を算出するために必要な全身用ボディコイルの画像は、信号結合処理801により作成する。801の具体的な処理として、例えば、各RF受信コイルで取得した画像si(x,y)と、si(x,y)にローパスフィルタを施した画像をwi(x,y)とするとき、合成後の擬似的なボディコイル画像s'c(x,y)は、
【００３４】
【数１３】
で求まる（ここで、*は複素共役、||は絶対値を表す）。しかし、マルチプルRF受信コイルの有感度領域は通常のコイルと比較して狭いので、この様に合成した擬似的なボディコイルの画像には、シェーディングの影響が大きい場合がある。そこで、合成した擬似的なボディコイルの画像に感度補正処理802を施し、シェーディングの無い擬似画像803を作成して、パラレルイメージング法の処理を行う。なお、感度補正処理802としては、例えば、画像からシェーディングを算出し補正する方法（特開平7-222724号公報）を用いる。
【００３５】
このように、全身用ボディコイルの信号を合成することで、例えば全身用ボディコイルが無い、或いは取りつけ不可能な装置や、使用可能なコイルのチャンネル数が少ない装置でも、パラレルイメージング法が可能となる。また、全身用ボディコイルを有する装置においても、マルチプルRF受信コイルの信号のみの処理となるため、装置構成やソフトウエア構成、データフローが簡略化できる点である。例えば、ハンバーガー型オープンMRIでは、照射コイルの感度が比較的低いため、照射コイルでの受信ができない場合がある。又、物理的に受信できても、パラレルイメージング法のためだけに受信用部品を付けると、MRI装置のコストアップになる。しかし、本実施例を適用すると、このようなデメリットを克服できる。感度分布の演算は、マルチスライス撮影や三次元撮影にも適用できる。マルチスライスではスライス毎に(13)式を計算する。三次元撮影でも、スライス毎に(13)式を計算する。
【００３６】
本発明は、以上の実施例で開示された内容にとどまらず、本発明の趣旨を踏まえた上で各種形態を取り得る。本実施例では、3つのマルチプルRF受信コイルを用いたパラレルイメージング法の例を示したが、RF受信コイルの数Nは2以上の任意数を使用可能である。また、倍速数M=2,3の例を示したが、倍速数はM≦Nの範囲内で選択可能である。さらに、本発明ではグラディエントエコーシーケンスについて記載したが、パラレルイメージング法では、シーケンスの形状には依存しない。例えば、SEシーケンス、FSEシーケンス、EPIシーケンス、スパイラルシーケンス、SSFPシーケンスなどに適用できる。また、三次元計測に本発明を適用する場合は、位相エンコード方向だけでなく、スライスエンコード方向にデータを間引いて高速化することもできる。或いは、位相エンコード方向、スライス方向を組み合わせてデータを間引き、高速化することもできる。更に、本アルゴリズムを適用して心臓イメージングを行うと、輝点アーチファクトの無い画像が高時間分解能で取得できる。
【００３７】
【発明の効果】
本発明は以上のように構成されたので、パラレルイメージング法で輝点アーチファクトを無くす事ができる。また、全身用ボディコイルが無い装置でもパラレルイメージング法が可能となる。さらに、全身用ボディコイルを有する装置においても、マルチプルRF受信コイルの信号のみでパラレルイメージング法が可能となるため、使用チャンネル数を低減でき、装置構成や処理・データフローを簡略化できる。
【図面の簡単な説明】
【図１】本発明の第1の実施形態における信号処理を説明する図。
【図２】パラレルイメージング法の画像の折り返しを説明する図。
【図３】本発明が適用されるRFコイルの受信部を示す図。
【図４】本発明が適用されるMRI装置の全体構成を示す図。
【図５】本発明における画像の折り返しを説明する図。
【図６】一般的なグラディエントエコーのシーケンスを説明する図。
【図７】従来のパラレルイメージング法の処理を説明する図。
【図８】本発明の第2の実施形態における信号処理を説明する図。
【符号の説明】
601 高周波パルス、602 スライス選択傾斜磁場、603 位相エンコード傾斜磁場パルス、604 読み出し傾斜磁場パルス、605 データサンプルウインド、606 エコー信号、607 繰り返し時間間隔、608 画像取得時間、401 被検体、402 磁石、403 傾斜磁場コイル、404 RFコイル、405 RFプローブ、406 信号検出部、407 信号処理部、408 表示部、409 傾斜磁場電源、410 RF送信部、411 制御部、412 ベッド
【発明の属する技術分野】
本発明は、被検体中の水素や燐等からの核磁気共鳴（以下、「NMR」という）信号を測定し、核の密度分布や緩和時間分布等を映像化する核磁気共鳴撮影（MRI）装置に関し、心臓等の高速撮影が要求される領域において、特に複数のRF受信コイルを用い、各RF受信コイルにて位相エンコードを間引いて取得した信号をRF受信コイルの感度分布を用いて行列演算により展開する撮影方法（以下、パラレルイメージング法）において、アーチファクトの無い画像を得ることを可能としたMRI装置に関する。
【０００２】
【従来の技術】
MRIでは、位相エンコード量を変えながらシーケンスを繰り返し実行し、１枚の画像再構成に必要なエコー信号を取得する。そのため、画像の取得時間は繰り返し回数が大きく影響する。高速撮影を行う場合、一般的には、１回の繰り返し内に複数のエコー信号を発生させるマルチエコータイプのシーケンスを用いたり、繰り返しの時間間隔を数〜数十msにまで短縮したシーケンスを用いる。しかし、マルチエコータイプのシーケンスでは、画像のコントラストが低下したり、画像歪みの原因となることがある。これは、画像のコントラストに寄与するエコー時間がマルチエコータイプのシーケンスでは各エコー信号で異なるため、コントラストの低下を招く。また、エコー時間が異なるため、各エコー信号間での位相変化が異なると、画像歪みとなって画像に現れる。
【０００３】
また、心臓領域（冠状動脈撮影など）を撮影する場合、より高速に画像を取得する必要があり、パラレルイメージング法と呼ばれる高速撮影方法も提案されている。パラレルイメージング法は、マルチプルRF受信コイルを用いて、位相エンコードを等間隔に間引いたシーケンスを実行し、繰り返し回数を低減して撮影時間を短縮する。通常、位相エンコードを等間隔に間引いて計測を行うと、画像には折り返しが発生するが、各RF受信コイルの感度分布をもとに行列演算を行うことで画像を展開し、折り返しを除去する。一般的にパラレルイメージング法では、撮影に用いたRF受信コイルの数分だけ撮影時間を短縮できる。
【０００４】
【発明が解決しようとする課題】
パラレルイメージング法では、マルチエコータイプのシーケンスを用いなくとも、撮影時間を短縮することができるため、コントラストの低下や画像歪みの影響を減らすことができる。しかしながら、マルチプルRF受信コイルの配置や感度分布の形状により結果の画像が大きく変化する。特に、背景等の低信号領域を含んだまま行列演算を行うと、ノイズの影響で折り返しの展開時に誤差が大きくなり、画像中に輝点のアーチファクトが発生することがある。また、各RF受信コイルの感度分布を精度良く算出するために、比較的感度分布の均一な全身用ボディコイルで計測した画像を用いて感度分布を算出するのが望ましいが、装置によっては全身用のボディコイルが無かったり、受信チャンネル数が少ないために同時に画像を取得できない問題が有る。
【０００５】
そこで、本発明の目的は、行列演算時の計算誤差を抑え、輝点アーチファクトの発生を無くし、また全身用ボディコイルの無い装置、或いはチャンネル数の少ない装置でもパラレルイメージングを可能とすることにある。あるいは、全身用ボディコイルを有する装置でも、装置構成及び信号処理フローを簡潔にすることにある。
【０００６】
【課題を解決するための手段】
上記目的を達成するために、本発明は、複数のＲＦ受信コイルを備えて被検体からの核磁気共鳴信号を受信する受信手段と、前記受信手段を用いてｋ空間のエンコードステップを間引いて前記核磁気共鳴信号を計測して、前記ＲＦ受信コイル毎のｋ空間データを取得する計測制御手段と、前記ＲＦ受信コイル毎の感度画像を用いて該ＲＦ受信コイル毎の感度分布データを取得する感度分布演算と、前記ＲＦ受信コイル毎のｋ空間データと感度分布データとを用いてパラレルイメージング法に基づく演算を適用して前記被検体の画像を再構成する画像再構成演算と、を行う演算手段と、を備えた磁気共鳴イメージング装置において、前記感度分布演算は、前記複数のＲＦ受信コイルの内の少なくとも２つの感度画像を合成して略均一感度画像を取得し、該略均一感度画像を用いて、前記複数のＲＦ受信コイルの内の少なくとも１つのＲＦ受信コイルの感度分布データを取得することを特徴とする。
【０００７】
また、少なくとも2個のRF受信コイルを備え、前記RF受信コイルで受信する核磁気共鳴信号の計測を計測空間のエンコードステップを間引くよう計測して各RF受信コイル毎に感度画像及び形態画像を取得し、前記RF受信コイルの感度画像に基づく感度分布から各形態画像の折り返し除去の演算を行なうと共に前記各形態画像を結合して１枚の形態画像を得る制御手段を有する磁気共鳴イメージング装置において、前記制御手段は、各RF受信コイルで取得した感度画像を結合した全体感度画像を作成し、前記全体感度画像を用いて各RF受信コイルの感度分布を算出し、前記全体感度画像を用いて画像の非画像領域と画像領域を分けるマスクを作成し、前記マスクに基づいて演算を行なうものである。
【０００８】
さらに、計測領域全体を受信する第1のRF受信コイルと、前記計測領域を少なくとも2つ以上分割した領域を受信する複数の第2のRF受信コイルを備え、前記各第2のRF受信コイルで受信する核磁気共鳴信号の計測を計測空間のエンコードステップを間引くよう計測して前記第2のRF受信コイル毎に感度画像及び形態画像を取得し、前記第1のRF受信コイルによる感度画像及び前記第2のRF受信コイルによる感度画像から求めた感度分布から前記形態画像の折り返し除去の演算を行なうと共に前記各形態画像を結合して１枚の形態画像を得る制御手段を有する磁気共鳴イメージング装置において、前記制御手段は、前記受信感度分布において非画像領域と画像領域を分けたマスクを作成し、前記マスクを用いて演算を行なうものである。
【０００９】
【発明の実施の形態】
以下、本発明の核磁気共鳴イメージング装置について、図面を参照して詳述する。図4は典型的な核磁気共鳴イメージング装置の構成である。被検体401の周囲に静磁場を発生する磁石402と、該空間に傾斜磁場を発生する傾斜磁場コイル403と、この領域に高周波磁場を発生するRFコイル404と被検体401が発生するMR信号を検出するRFプローブ405がある。傾斜磁場コイル403は、X、Y、Zの３方向の傾斜磁場コイルで構成され、傾斜磁場電源409からの信号に応じてそれぞれ傾斜磁場を発生する。RFコイル404はRF送信部410の信号に応じて高周波磁場を発生する。RFプローブ405の信号は、信号検出部406で検出され、信号処理部407で信号処理され、また計算により画像信号に変換される。画像は表示部408で表示される。傾斜磁場電源409、RF送信部410、信号検出部406は制御部411で制御され、制御のタイムチャートは一般にパルスシーケンスと呼ばれている。ベッド412は被検体が横たわるためのものである。
【００１０】
現在MRIの撮影対象は、臨床で普及しているものとしては、被検体の主たる構成物質、プロトンである。プロトン密度の空間分布や、励起状態の緩和現象の空間分布を画像化することで、人体頭部、腹部、四肢等の形態または、機能を2次元もしくは3次元的に撮影する。
【００１１】
RFコイル405の例として、複数の受信コイルを用いた「マルチプルRFコイル」もしくは「フェーズドアレイコイル」と呼ばれる技術が用いられている。マルチプルRFコイルとは、相対的に高感度な小型RF受信コイルを複数個並べて、各コイルで取得した信号を合成することより、RF受信コイルの高い感度を保ったまま視野を拡大し、高感度化を図る受信専用RFコイルである。水平磁場頭部用マルチプルRFコイルとしては、Array Head Coil for Improved Functional MRI(Christoph Leussler), 1996 ISMRM abstruct p.249がある。また、水平磁場頭部用QDマルチプルRFコイルとしては、Helmet and Cylindrical Shaped CP Array Coils for Brain Imaging: A Comparison of Signal-to-Noise Characteristics(H.A.Stark, E.M.Haacke), 1996 ISMRM abstract p.1412がある。水平磁場腹部用QDマルチプルRFコイルとしては、Four Channel Wrap-Around Coil with Inductive Decoupler for 1.5T Body Imaging(T.Takahashi et.al), 1995 ISMRM abstruct p.1418がある。
【００１２】
マルチプルRFコイルの信号検出部の一部を図3に示す。図3では、４個のRF受信コイル405が、それぞれプリアンプ302に接続されて一つのマルチプルコイル301を構成する。信号検出部406は、4個のAD変換・直交検波回路303が並列してなり、前記各プリアンプの出力が接続されている。信号処理部407は、各信号をフーリエ変換、バックプロジェクション法、ウエーブレット変換などにより、それぞれのRF受信コイル405で検出したMRI画像を求め、これらの画像信号を合成する演算304からなる。
【００１３】
次に、撮影方法を説明する。図6は一般的なグラディエントエコーシーケンスである。601は高周波パルス、602はスライス選択傾斜磁場パルス、603は位相エンコード傾斜磁場パルス、604は読み出し傾斜磁場パルス、605はサンプリングウインド、606はエコー信号、607は繰り返し時間（高周波パルス601の照射間隔）である。MRIでは、各繰り返し時間607毎に位相エンコード傾斜磁場パルス603の量を変え、異なる位相エンコードを与え、それぞれの位相エンコードで得られるエコー信号606を検出する。この操作を位相エンコードの数だけ繰り返し、画像取得時間608で1枚の画像再構成に必要なエコー信号を取得する。位相エンコードの数は通常1枚の画像あたり64、128、256、512等の値が選ばれる。各エコー信号は通常128、256、512、1024個のサンプリングデータからなる時系列信号として得られる。これらのデータを2次元フーリエ変換して1枚のMR画像を作成する。
【００１４】
マルチプルRF受信コイルを用いた高速撮影であるパラレルイメージング法の場合、位相エンコードステップ間隔を一定の割合で間引いて撮影の繰り返し回数を減らす。この間引き率を一般的に倍速数と呼ぶ。例えば、位相エンコードステップを2倍に間引く場合、倍速数は2である。以下、パラレルイメージング法の原理を、図2を用いて説明する。通常撮影の場合、図2(a)に示すように各位相エンコード量で取得した信号202nを配置し画像１枚分のデータ201を形成し、これをフーリエ変換することで図2(c)示すような画像を得る。次に、パラレルイメージング法において、例えば位相エンコードステップ間隔を2倍にし、データを間引いた場合、図2(b)に示すように、1ラインおきにデータ204mを計測し、205mの位置に相当するデータは計測しない。このとき、計測したデータ204mの量は通常撮影と比べ半分になるので、マトリクスを半分にして画像を作成するが、図2(d)のような折り返しの発生した画像が得られる。図2においては、y方向を位相エンコード方向としており、この折り返しは、位相エンコード方向の画像が折り返して発生する。つまり、図2(c)に示すような画像207内の被検体206の、上側領域2071内の被検体画像2061と下側領域2072内の被検体画像2062が重なり、図2(d)に示すような折り返しの有る画像208となる。このようにして発生した信号折り返しは、SENSEと呼ばれる信号処理法で除去する（SENSE:Sensitivity Encoding for Fast MRI(Klaas P.Pruessmann et.al), Magnetic Resonance in Medicine 42: 952-962 (1999)参照）。
【００１５】
以下、折り返し除去について説明する。x,y方向の画像マトリクスをそれぞれX,Yとした場合、画像内の座標(x,y)(x:1≦x≦X, y:1≦y≦Y)における画素値をsi (x,y)とする(ここで、iはコイル番号であり、2≦i≦Nである)。図2(d)の場合では、位相エンコードステップを2倍に間引いているので、間引き後の画像の位相エンコード方向のマトリクスは、Y'≡Y/2となる。図2(d)の画像の座標を(x,y')(1≦y'≦Y')としたとき、その画素値s'(x,y')208は、元画像207の2領域2071、2072が重なって、
【００１６】
【数１】
【００１７】
【数２】
【００１８】
【数３】
【数４】
【００２０】
【数５】
【００２１】
【数６】
【００２２】
同様に、N個のコイルを用いて、M倍速で撮影を行った場合は、Y'≡Y/M、1≦y'≦Y'として、
【００２３】
【数７】
【００２４】
通常、パラレルイメージング法では、各RF受信コイルの感度分布Cijを前計測等で予め計測するなどして取得する。しかしながら、感度分布Cijを直接算出する事は難しく、一般的には比較的感度分布の均一な全身用ボディコイルの画像を用いて、各RF受信コイルで取得した画像をそれぞれ除算し、近似的なコイルの感度分布を求めて行列演算を行う。
【００２５】
図7に一般的なパラレルイメージング法の処理を示す。図では1つの全身用ボディコイルと、3つのマルチプルRF受信コイルの構成である。まず、全身用ボディコイルで取得した感度画像701と、各RF受信コイルで取得した感度画像7021〜7023を用いて、感度分布算出処理7041〜7043を行い、各RF受信コイルの感度分布7051〜7053を得る。感度分布算出処理704としては、全身用ボディコイルの画像をsc(x,y)として、例えば、
【００２６】
【数８】
次に、本発明の第1の実施形態を図1，図5により説明する。従来の一般的なパラレルイメージング法では、背景等の低信号領域に関しては何ら考慮されていないものであるが、本実施形態では低信号領域を考慮した構成となっている。つまり、全身用ボディコイルの画像を用いて、マスク作成のステップ101により、マスク102を作成し、マスクの処理を行うステップ103により精度を向上した処理後のマスク104を作成し、マスクを用いた行列作成処理105を行うよう構成した。ここで、マスク処理について説明する。図5では3倍速の場合であり、画像501は通常画像(図5(b))の3領域5041、5042、5043が重なった結果である（ここで、Y'≡Y/3、1≦y'≦Y'である）。図5(a)中に、2箇所の注目点A(xA, y'A) 502、B(xB, y'B) 503を設ける。このとき、A、B点の画素値は、それぞれ、図の5021〜5023、5031〜5033の領域が重なっている事から、
【００２７】
【数９】
【００２８】
ここで、マスク処理103としては、例えば被検体内部に構造が存在する場合はその領域の信号をそのまま使用するためマスク領域を内挿により埋めたり、本来計算に入れない領域に誤って取られた孤立点を除去する処理等を行なう。
また、マスクを用いた行列作成処理105では、行列の要素に背景領域を含まないように、各要素を0にする判別処理を行なう。例えば、
【００２９】
【数１０】
【００３０】
【数１１】
さらに、行列作成処理105内で(10)式に示すような各要素を0にする判別処理を用いることは、感度分布とマスクとの積、
【００３１】
【数１２】
(10)式のように行列要素の一部を"0"とすることは、行列のマトリクスサイズを縮小することと等価である。すなわち、「各画素毎に折り返しの状況を調べ、必要十分な行列サイズを決定し、必要十分な演算を行う」ことになる。この結果、従来と比べ、不要なノイズの混入が無くなり、画質が向上する。又、演算時間が短縮するメリットも有る。
【００３２】
なお、行列作成処理105では処理後のマスク104を用いて演算したが、マスク102を用いて演算することもできる。これにより処理フローを簡略化することができる。また、マスク102は背景領域と被検体領域に分けて説明したが、背景領域とは被検体領域外の背景のみならず被検体領域中に存在する低信号領域（例えば、空隙等）も含む。
【００３３】
次に、本発明の第2の実施形態を図8を用いて説明する。この場合では、全身用ボディコイルの画像がない場合である。この場合、各RF受信コイルの感度分布を算出するために必要な全身用ボディコイルの画像は、信号結合処理801により作成する。801の具体的な処理として、例えば、各RF受信コイルで取得した画像si(x,y)と、si(x,y)にローパスフィルタを施した画像をwi(x,y)とするとき、合成後の擬似的なボディコイル画像s'c(x,y)は、
【００３４】
【数１３】
【００３５】
このように、全身用ボディコイルの信号を合成することで、例えば全身用ボディコイルが無い、或いは取りつけ不可能な装置や、使用可能なコイルのチャンネル数が少ない装置でも、パラレルイメージング法が可能となる。また、全身用ボディコイルを有する装置においても、マルチプルRF受信コイルの信号のみの処理となるため、装置構成やソフトウエア構成、データフローが簡略化できる点である。例えば、ハンバーガー型オープンMRIでは、照射コイルの感度が比較的低いため、照射コイルでの受信ができない場合がある。又、物理的に受信できても、パラレルイメージング法のためだけに受信用部品を付けると、MRI装置のコストアップになる。しかし、本実施例を適用すると、このようなデメリットを克服できる。感度分布の演算は、マルチスライス撮影や三次元撮影にも適用できる。マルチスライスではスライス毎に(13)式を計算する。三次元撮影でも、スライス毎に(13)式を計算する。
【００３６】
本発明は、以上の実施例で開示された内容にとどまらず、本発明の趣旨を踏まえた上で各種形態を取り得る。本実施例では、3つのマルチプルRF受信コイルを用いたパラレルイメージング法の例を示したが、RF受信コイルの数Nは2以上の任意数を使用可能である。また、倍速数M=2,3の例を示したが、倍速数はM≦Nの範囲内で選択可能である。さらに、本発明ではグラディエントエコーシーケンスについて記載したが、パラレルイメージング法では、シーケンスの形状には依存しない。例えば、SEシーケンス、FSEシーケンス、EPIシーケンス、スパイラルシーケンス、SSFPシーケンスなどに適用できる。また、三次元計測に本発明を適用する場合は、位相エンコード方向だけでなく、スライスエンコード方向にデータを間引いて高速化することもできる。或いは、位相エンコード方向、スライス方向を組み合わせてデータを間引き、高速化することもできる。更に、本アルゴリズムを適用して心臓イメージングを行うと、輝点アーチファクトの無い画像が高時間分解能で取得できる。
【００３７】
【発明の効果】
本発明は以上のように構成されたので、パラレルイメージング法で輝点アーチファクトを無くす事ができる。また、全身用ボディコイルが無い装置でもパラレルイメージング法が可能となる。さらに、全身用ボディコイルを有する装置においても、マルチプルRF受信コイルの信号のみでパラレルイメージング法が可能となるため、使用チャンネル数を低減でき、装置構成や処理・データフローを簡略化できる。
【図面の簡単な説明】
【図１】本発明の第1の実施形態における信号処理を説明する図。
【図２】パラレルイメージング法の画像の折り返しを説明する図。
【図３】本発明が適用されるRFコイルの受信部を示す図。
【図４】本発明が適用されるMRI装置の全体構成を示す図。
【図５】本発明における画像の折り返しを説明する図。
【図６】一般的なグラディエントエコーのシーケンスを説明する図。
【図７】従来のパラレルイメージング法の処理を説明する図。
【図８】本発明の第2の実施形態における信号処理を説明する図。
【符号の説明】
601 高周波パルス、602 スライス選択傾斜磁場、603 位相エンコード傾斜磁場パルス、604 読み出し傾斜磁場パルス、605 データサンプルウインド、606 エコー信号、607 繰り返し時間間隔、608 画像取得時間、401 被検体、402 磁石、403 傾斜磁場コイル、404 RFコイル、405 RFプローブ、406 信号検出部、407 信号処理部、408 表示部、409 傾斜磁場電源、410 RF送信部、411 制御部、412 ベッド
Claims (6)
- 複数のＲＦ受信コイルを備えて被検体からの核磁気共鳴信号を受信する受信手段と、
前記受信手段を用いてｋ空間のエンコードステップを間引いて前記核磁気共鳴信号を計測して、前記ＲＦ受信コイル毎のｋ空間データを取得する計測制御手段と、
前記ＲＦ受信コイル毎の感度画像を用いて該ＲＦ受信コイル毎の感度分布データを取得する感度分布演算と、前記ＲＦ受信コイル毎のｋ空間データと感度分布データとを用いてパラレルイメージング法に基づく演算を適用して前記被検体の画像を再構成する画像再構成演算と、を行う演算手段と、
を備えた磁気共鳴イメージング装置において、
前記感度分布演算は、前記複数のＲＦ受信コイルの内の少なくとも２つの感度画像を合成して略均一感度画像を取得し、該略均一感度画像を用いて、前記複数のＲＦ受信コイルの内の少なくとも２つのＲＦ受信コイルの感度分布データを取得し、所定の閾値を用いて前記略均一感度画像における低信号強度領域と高信号強度領域とを分けるマスクを作成し、該マスクを用いて前記感度分布データの変換を行うことを特徴とする磁気共鳴イメージング装置。 - 前記感度分布演算は、前記マスクを用いた感度分布データの変換において、マスクされるデータを所定の値に変換することを特徴とする請求項１記載の磁気共鳴イメージング装置。
- 前記画像再構成演算は、前記略均一感度画像における前記所定の閾値より大きい高信号強度領域のみに前記パラレルイメージング法に基づく演算を適用することを特徴とする請求項１又は２に記載の磁気共鳴イメージング装置。
- 前記画像再構成演算は、行列を用いて前記パラレルイメージング法に基づく演算を行い、
前記行列のマトリクスサイズは、前記高信号強度領域のみが反映されるように調整されていることを特徴とする請求項３記載の磁気共鳴イメージング装置。 - 前記感度分布演算は、前記閾値として、前記低信号強度領域が被検体画像の背景領域となるような値を用いることを特徴とする請求項１乃至４のいずれか一項に記載の磁気共鳴イメージング装置。
- 前記計測制御手段は、前記各コイルの感度画像用エコー信号の計測を行い、
前記演算手段は、前記コイル毎の感度画像用エコー信号を用いて、該コイル毎の感度画像を取得することを特徴とする請求項１乃至５のいずれか一項に記載の磁気共鳴イメージング装置。
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001122213A JP4047553B2 (ja) | 2001-04-20 | 2001-04-20 | 磁気共鳴イメージング装置 |
| PCT/JP2002/001653 WO2002065907A1 (fr) | 2001-02-23 | 2002-02-25 | Procede et appareil d'imagerie par resonance magnetique |
| US10/468,963 US6876201B2 (en) | 2001-02-23 | 2002-02-25 | Magnetic resonance imaging apparatus and method |
| EP09014853A EP2159590A1 (en) | 2001-02-23 | 2002-02-25 | Magnetic resonance imaging apparatus and method |
| EP02703889A EP1371327A4 (en) | 2001-02-23 | 2002-02-25 | METHOD AND APPARATUS FOR MAGNETIC RESONANCE IMAGING |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001122213A JP4047553B2 (ja) | 2001-04-20 | 2001-04-20 | 磁気共鳴イメージング装置 |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JP2002315731A JP2002315731A (ja) | 2002-10-29 |
| JP2002315731A5 JP2002315731A5 (ja) | 2005-09-22 |
| JP4047553B2 true JP4047553B2 (ja) | 2008-02-13 |
Family
ID=18971956
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001122213A Expired - Fee Related JP4047553B2 (ja) | 2001-02-23 | 2001-04-20 | 磁気共鳴イメージング装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4047553B2 (ja) |
Families Citing this family (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4443079B2 (ja) | 2001-09-13 | 2010-03-31 | 株式会社日立メディコ | 磁気共鳴イメージング装置及び磁気共鳴イメージング装置用ｒｆ受信コイル |
| JP4068954B2 (ja) * | 2002-12-26 | 2008-03-26 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ｍｒｉ装置 |
| JP3929047B2 (ja) * | 2003-04-24 | 2007-06-13 | 株式会社日立メディコ | 磁気共鳴イメージング装置 |
| US6919722B2 (en) * | 2003-10-09 | 2005-07-19 | Ge Medical Systems Global Technology Company, Llc | Image quality improvement for SENSE with low signal regions |
| JP4698231B2 (ja) * | 2004-06-11 | 2011-06-08 | 株式会社日立メディコ | 磁気共鳴診断装置 |
| JP4679158B2 (ja) * | 2004-10-08 | 2011-04-27 | 株式会社日立メディコ | 磁気共鳴イメージング装置 |
| JP4975614B2 (ja) * | 2005-03-30 | 2012-07-11 | 株式会社日立メディコ | 磁気共鳴イメージング装置及び方法 |
| US7714577B2 (en) * | 2005-10-18 | 2010-05-11 | Tursiop Technologies Llc | Method and apparatus for high-gain magnetic resonance imaging |
| JP4980662B2 (ja) * | 2006-06-30 | 2012-07-18 | 株式会社日立メディコ | 磁気共鳴イメージング装置 |
| JP4817381B2 (ja) * | 2006-10-17 | 2011-11-16 | 株式会社日立メディコ | 磁気共鳴イメージング装置 |
| JP5156958B2 (ja) * | 2007-07-17 | 2013-03-06 | 株式会社日立メディコ | 磁気共鳴イメージング装置 |
| JP5686533B2 (ja) * | 2009-07-10 | 2015-03-18 | 株式会社東芝 | 磁気共鳴イメージング装置、画像処理装置および画像処理方法 |
| US9864032B2 (en) * | 2010-01-05 | 2018-01-09 | National Health Research Institutes | Magnetic resonance imaging system |
| DE102011088828B4 (de) * | 2011-12-16 | 2013-08-29 | Siemens Aktiengesellschaft | Erstellung eines MR-Bildes eines Untersuchungsobjekts unter Verwendung einer für ein Empfangsspulenelement erstellten Maske |
| US10534060B2 (en) | 2013-06-10 | 2020-01-14 | Toshiba Medical Systems Corporation | Parallel MRI with spatially misregistered signal |
| JP6518559B2 (ja) * | 2015-09-18 | 2019-05-22 | 株式会社日立製作所 | 磁気共鳴イメージング装置、及び方法 |
| JP7608268B2 (ja) | 2021-05-20 | 2025-01-06 | 富士フイルム株式会社 | 磁気共鳴撮影装置及び感度分布算出プログラム |
-
2001
- 2001-04-20 JP JP2001122213A patent/JP4047553B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2002315731A (ja) | 2002-10-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US6876201B2 (en) | Magnetic resonance imaging apparatus and method | |
| JP4152381B2 (ja) | 磁気共鳴イメージング装置 | |
| EP1830198B1 (en) | Data correction apparatus and method | |
| US8526703B2 (en) | Making different corrections for different areas of magnetic resonance image data | |
| JP3952247B2 (ja) | 核磁気共鳴撮影装置 | |
| EP3004908B1 (en) | Parallel mri with b0 distortion correction and multi-echo dixon water-fat separation using regularised sense reconstruction | |
| JP4047553B2 (ja) | 磁気共鳴イメージング装置 | |
| CN102670201B (zh) | 用于产生磁共振图像数据的方法和磁共振断层造影系统 | |
| US7385397B2 (en) | Magnetic resonance imaging device and data processing method for magnetic resonance imaging device | |
| JP2003079595A (ja) | 磁気共鳴イメージング装置及び磁気共鳴イメージング装置用ｒｆ受信コイル | |
| CN102597795A (zh) | 使用导航器的mr成像 | |
| EP1372110B1 (en) | Method and system for image reconstruction | |
| US20060058629A1 (en) | Magnetic resonance method | |
| US9157978B2 (en) | Magnetic resonance imaging apparatus and method | |
| JP5536358B2 (ja) | 磁気共鳴イメージング装置及び感度補正方法 | |
| JP4679158B2 (ja) | 磁気共鳴イメージング装置 | |
| JP4975614B2 (ja) | 磁気共鳴イメージング装置及び方法 | |
| JP4980662B2 (ja) | 磁気共鳴イメージング装置 | |
| JP2009022319A (ja) | 磁気共鳴イメージング装置 | |
| JP4675936B2 (ja) | 核磁気共鳴撮影装置 | |
| WO2018001759A1 (en) | Diffusion weighted mr imaging using multi-shot epi with motion detection and modified sense reconstruction | |
| JP4817381B2 (ja) | 磁気共鳴イメージング装置 | |
| JP3983792B2 (ja) | 核磁気共鳴撮影装置 | |
| JP3952310B2 (ja) | 核磁気共鳴撮影装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050413 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050413 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20061121 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20061226 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070205 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070307 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070507 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070703 |
|
| 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: 20071119 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20071122 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101130 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 4047553 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101130 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111130 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111130 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121130 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131130 Year of fee payment: 6 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| LAPS | Cancellation because of no payment of annual fees |