JP2008017958A - パラレルイメージング方法およびｍｒｉ装置 - Google Patents

パラレルイメージング方法およびｍｒｉ装置 Download PDF

Info

Publication number
JP2008017958A
JP2008017958A JP2006190999A JP2006190999A JP2008017958A JP 2008017958 A JP2008017958 A JP 2008017958A JP 2006190999 A JP2006190999 A JP 2006190999A JP 2006190999 A JP2006190999 A JP 2006190999A JP 2008017958 A JP2008017958 A JP 2008017958A
Authority
JP
Japan
Prior art keywords
channel
receiving coil
coil
image
row
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2006190999A
Other languages
English (en)
Other versions
JP4464373B2 (ja
Inventor
Yoshikazu Ikezaki
吉和 池崎
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.)
GE Medical Systems Global Technology Co LLC
Original Assignee
GE Medical Systems Global Technology Co LLC
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 GE Medical Systems Global Technology Co LLC filed Critical GE Medical Systems Global Technology Co LLC
Priority to JP2006190999A priority Critical patent/JP4464373B2/ja
Priority to US11/776,115 priority patent/US8054078B2/en
Priority to CN2007101287812A priority patent/CN101105522B/zh
Publication of JP2008017958A publication Critical patent/JP2008017958A/ja
Application granted granted Critical
Publication of JP4464373B2 publication Critical patent/JP4464373B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/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)

Abstract

【課題】３個以上の受信コイルを用いたパラレルイメージングにおいて、少ない演算量かつ短い処理時間で合成画像を得る。
【解決手段】３個以上の受信コイルを用い且つ位相エンコードステップを間引いたスキャンを行って各受信コイルのデータを収集し、各データから画像ａiをそれぞれ生成する（ステップＱ２）。各画像ａiの組合せの中から折返しを元に戻す演算に用いる画像の組合せを選び、選んだ画像の組合せ及び対応する受信コイルの感度係数の正方行列を利用した演算により合成画像を得る（ステップＱ３）。
【効果】無視してもほとんど影響を受けない画像の組合せについては演算しないので、演算量を減らすことが出来ると共に処理時間を短縮できる。
【選択図】図３

Description

本発明は、パラレルイメージング方法およびＭＲＩ（Magnetic Resonance Imaging）装置に関し、さらに詳しくは、３個以上の受信コイルを用いて、少ない演算量かつ短い処理時間で折返しを元に戻した合成画像を得ることが出来るパラレルイメージング方法およびＭＲＩ装置に関する。
従来、複数個の受信コイルを用いて位相エンコードステップを間引いたスキャンを行って各受信コイルのデータを収集し、各データから各画像を生成し、各受信コイルの感度分布の差を利用した折返しを元に戻す（unfolding）演算により折返しを元に戻した合成画像を得るパラレルイメージング方法が知られている（例えば、特許文献１、非特許文献１参照。）。
特開２００３−７９５９５号公報 Klaas P. Pruessmann et al."SENSE:Sensitivity Encoding for Fast MRI"Magnetic Resonance in Medicine 42(1999) pp.952-962
従来の折返しを元に戻す演算では、受信コイル数がＮ（≧３）個になると、受信コイル数が２個の場合の演算量に対して、演算量がＮ個から２個を取り出す組合せ数倍になり、演算量が急激に増加してしまう問題点があった。例えば、受信コイル数が３個になると演算量が３倍に増加し、受信コイル数が４個になると演算量が６倍に増加してしまう。
そこで、本発明の目的は、３個以上の受信コイルを用いて、少ない演算量かつ短い処理時間で折返しを元に戻した合成画像を得ることが出来るパラレルイメージング方法およびＭＲＩ装置を提供することにある。
第１の観点では、本発明は、３個以上の受信コイルを用い且つ位相エンコードステップを間引いたスキャンを行って前記各受信コイルのデータを収集し、前記各データから画像をそれぞれ生成し、前記各画像の組合せの中から折返しを元に戻す演算に用いる画像の組合せを選び、選んだ画像の組合せ及び対応する受信コイルの感度係数の正方行列を利用した演算により合成画像を得ることを特徴とするパラレルイメージング方法を提供する。
一般的に、パラレルイメージング方法では、Ｎ個の受信コイルの各校正データｃ(ｎ)から複素数画像Ｃ(ｎ)をそれぞれ生成し、ボディコイルの校正データｃ(０)から複素数画像Ｃ(０)を生成し、各受信コイルの複素数画像Ｃ(ｎ)をボディコイルの複素数画像Ｃ(０)で割り算することにより各受信コイルの感度マップｓ(ｎ)を得る。そして、各受信コイルの感度マップｓ(ｎ)を順に並べた感度マトリクスＳと各受信コイルのデータｈ(ｎ)から生成した複素数画像ａ(ｎ)を順に並べた画像マトリクスＡとから次式により合成画像Ｖを求める。
Ｖ＝（Ｓ*Ψ-1Ｓ）-1Ｓ*Ψ-1Ａ
上式で、Ｓ*は、Ｓの随伴行列（conjugate transpose）である。Ψは、noise correlation matrix である。この計算は、画素毎に行われる。
なお、上式は、非特許文献１（Klaas P. Pruessmann et al.“SENSE:Sensitivity Encoding for Fast MRI”Magnetic Resonance in Medicine 42(1999) pp.952-962）に記載されている。
説明の簡単化のために、noise correlation matrix を使用しないとすると、（数１）となる。
Figure 2008017958
例えば、受信コイル数Ｎ＝３とし、リダクションファクターＲ＝２とすると、（数２）となる。
Figure 2008017958
 ａiは、ｉチャンネルの受信コイルの画像上の点の画素値である。
ｓi1，ｓi2は、ｉチャンネルの受信コイルの画像上の点で折返しにより重なっている元の２点のｉチャンネルの受信コイルに対応する各感度係数である。
（数１）に（数２）を入れて展開すると、（数３）となる。
Figure 2008017958
 |Ｓkm|2は、ｋチャンネルの受信コイルの感度係数ｓk1，ｓk2を第１行とし、ｍチャンネルの受信コイルの感度係数ｓm1，ｓm2を第２行とする行列式の２乗である。
Ｓkmは、ｋチャンネルの受信コイルの感度係数ｓk1，ｓk2を第１行とし、ｍチャンネルの受信コイルの感度係数ｓm1，ｓm2を第２行とする正方行列である。
（数３）から、受信コイル数がＮ（≧３）個になると、受信コイル数が２個の場合の演算量に対して、演算量がＮ個から２個を取り出す組合せ数だけ倍増し、演算量が急激に増加してしまうことが判る。
ここで、（数３）における|Ｓ12|2，|Ｓ13|2，|Ｓ23|2に着目すると、もしも|Ｓ12|2が|Ｓ13|2や|Ｓ23|2に比べて小さいなら、|Ｓ12|2の項を無視してもほとんど影響を受けない。すなわち、（数４）とすることが出来る。
Figure 2008017958
（数４）は、各画像の組合せＡ12，Ａ13，Ａ23の中から折返しを元に戻す演算に用いる画像の組合せＡ13，Ａ23を選び、選んだ画像の組合せＡ13，Ａ23及び対応する受信コイルの感度係数の正方行列Ｓ13，Ｓ23を利用した演算により合成画像Ｖを得ることに外ならない。
上記第１の観点によるパラレルイメージング方法では、例えば（数４）により合成画像を得るが、これは（数３）に比べると、演算量が少なくて済む。よって、短い処理時間で折返しを元に戻した合成画像を得ることが出来る。
第２の観点では、本発明は、上記第１の観点によるパラレルイメージング方法において、前記受信コイルの各組合せに対応する受信コイルの感度係数の正方行列の行列式の２乗に基づいて、前記各画像の組合せの中から折返しを元に戻す演算に用いる画像の組合せを選ぶことを特徴とするパラレルイメージング方法。
上記第２の観点によるパラレルイメージング方法では、１回だけ受信コイルの各組合せに対応する受信コイルの感度係数の正方行列の行列式の２乗を計算する必要があるが、各画像の組合せの中から折返しを元に戻す演算に用いる画像の組合せを正確に選ぶことが出来る。
第３の観点では、本発明は、前記第１の観点によるパラレルイメージング方法において、前記受信コイルの配置および前記フェーズエンコード方向に基づいて、前記各画像の中から折返しを元に戻す演算に用いる画像の組合せを選ぶことを特徴とするパラレルイメージング方法を提供する。
２つの受信コイルの配置がフェーズエンコード方向に対向している場合のこれら受信コイルの感度係数の正方行列の行列式の２乗は、周波数エンコード方向に並んでいる場合のこれら受信コイルの感度係数の正方行列の行列式の２乗よりもずっと大きな値になる。つまり、受信コイルの配置がフェーズエンコード方向に対向している受信コイルの画像の組合せを選び、受信コイルの配置が周波数エンコード方向に並んでいる受信コイルの画像の組合せを選ばなければよい。
上記第３の観点によるパラレルイメージング方法では、各受信コイルの感度係数の正方行列の行列式の２乗を計算しなくても、受信コイルの配置およびフェーズエンコード方向に基づいて、演算に用いる画像の組合せを選ぶことが出来る。
第４の観点では、本発明は、前記第１から前記第３のいずれかの観点によるパラレルイメージング方法において、前記受信コイルが、第１の受信コイルと、周波数エンコード方向に前記第１の受信コイルと並んだ第２の受信コイルと、位相エンコード方向に前記第１の受信コイルと対向する第３の受信コイルとからなり、リダクションファクターＲ＝２とすることを特徴とするパラレルイメージング方法を提供する。
上記第４の観点によるパラレルイメージング方法では、３個の受信コイルを用いて、少ない演算量かつ短い処理時間で折返しを元に戻した合成画像を得ることが出来る。
第５の観点では、本発明は、前記第４の観点によるパラレルイメージング方法において、ｉチャンネルの受信コイルの画像上の点の画素値をａiとし、ｉチャンネルの受信コイルの画像上の点で折返しにより重なっている元の２点のｉチャンネルの受信コイルに対応する各感度係数をｓi1，ｓi2とし、ｋチャンネルの受信コイルの感度係数ｓk1，ｓk2を第１行とすると共にｍチャンネルの受信コイルの感度係数ｓm1，ｓm2を第２行とする行列式の２乗を|Ｓkm|2とし、ｋチャンネルの受信コイルの感度係数ｓk1，ｓk2を第１行とし、ｍチャンネルの受信コイルの感度係数ｓm1，ｓm2を第２行とする正方行列をＳkmとするとき、前記演算が、
Figure 2008017958
により表されることを特徴とするパラレルイメージング方法を提供する。
上記第５の観点によるパラレルイメージング方法では、従来方法に比べて演算量を２／３に低減できる。
第６の観点では、本発明は、前記第１から前記第３のいずれかの観点によるパラレルイメージング方法において、前記受信コイルが、第１の受信コイルと、周波数エンコード方向に前記第１の受信コイルと並んだ第２の受信コイルと、位相エンコード方向に前記第１の受信コイルと対向する第３の受信コイルと、周波数エンコード方向に前記第３の受信コイルと並んだ第４の受信コイルとからなり、リダクションファクターＲ＝２とすることを特徴とするパラレルイメージング方法を提供する。
上記第６の観点によるパラレルイメージング方法では、４個の受信コイルを用いて、少ない演算量かつ短い処理時間で折返しを元に戻した合成画像を得ることが出来る。
第７の観点では、本発明は、前記第６の観点によるパラレルイメージング方法において、ｉチャンネルの受信コイルの画像上の点の画素値をａiとし、ｉチャンネルの受信コイルの画像上の点で折返しにより重なっている元の２点のｉチャンネルの受信コイルに対応する各感度係数をｓi1，ｓi2とし、ｋチャンネルの受信コイルの感度係数ｓk1，ｓk2を第１行とすると共にｍチャンネルの受信コイルの感度係数ｓm1，ｓm2を第２行とする行列式の２乗を|Ｓkm|2とし、ｋチャンネルの受信コイルの感度係数ｓk1，ｓk2を第１行とし、ｍチャンネルの受信コイルの感度係数ｓm1，ｓm2を第２行とする正方行列をＳkmとするとき、前記演算が、
Figure 2008017958
により表されることを特徴とするパラレルイメージング方法を提供する。
上記第７の観点によるパラレルイメージング方法では、従来方法に比べて演算量を２／３に低減できる。
第８の観点では、本発明は、前記第１から前記第３のいずれかの観点によるパラレルイメージング方法において、前記選んだ画像の組合せに対応する受信コイルを第１チャンネルから第Ｌチャンネルまでとし、リダクションファクターをＲとするとき、第１チャンネルから第Ｌチャンネルまでの中からＲチャンネルを選んだ各組合せを一般的に第1Rチャネルから第RRチャネルとし、ｉチャンネルの受信コイルの画像上の点の画素値をａiとし、ｉチャンネルの受信コイルの画像上の点で折返しにより重なっている元のＲ点のｉチャンネルの受信コイルに対応する各感度係数をｓi1，…，ｓiRとし、ｋチャンネルの受信コイルの感度係数ｓk1，…，ｓkRを各行とする正方行列の行列式の２乗を|Ｓ1R-RR|2とし、ｋチャンネルの受信コイルの感度係数ｓk1，…，ｓkRを各行とする正方行列をＳ1R-RRとするとき、前記演算が、
Figure 2008017958
により表されることを特徴とするパラレルイメージング方法を提供する。
上記第８の観点によるパラレルイメージング方法では、従来方法に比べて演算量を低減できる。
第９の観点では、本発明は、３個以上の受信コイルと、前記３個以上の受信コイルを用い且つ位相エンコードステップを間引いたスキャンを行って前記各受信コイルのデータを収集するスキャン手段と、前記各データから画像をそれぞれ生成する画像生成手段と、前記各画像の組合せの中から折返しを元に戻す演算に用いる画像の組合せを選ぶ組合せ選択手段と、選んだ画像の組合せ及び前記各受信コイルの感度分布を利用して前記演算により合成画像を得る演算手段とを具備したことを特徴とするＭＲＩ装置を提供する。
上記第９の観点によるＭＲＩ装置では、前記第１の観点によるパラレルイメージング方法を好適に実施できる。
第１０の観点では、本発明は、前記第９の観点によるＭＲＩ装置において、前記組合せ選択手段は、前記受信コイルの各組合せに対応する受信コイルの感度係数の正方行列の行列式の２乗に基づいて、前記各画像の組合せの中から折返しを元に戻す演算に用いる画像の組合せを選ぶことを特徴とするＭＲＩ装置を提供する。
上記第１０の観点によるＭＲＩ装置では、前記第２の観点によるパラレルイメージング方法を好適に実施できる。
第１１の観点では、本発明は、前記第９の観点によるＭＲＩ装置において、前記組合せ選択手段は、前記受信コイルの配置および前記フェーズエンコード方向に基づいて、前記各画像の中から折返しを元に戻す演算に用いる画像の組合せを選ぶことを特徴とするＭＲＩ装置を提供する。
上記第１１の観点によるＭＲＩ装置では、前記第３の観点によるパラレルイメージング方法を好適に実施できる。
第１２の観点では、本発明は、前記第９から前記第１１のいずれかの観点によるＭＲＩ装置において、前記受信コイルが、第１の受信コイルと、周波数エンコード方向に前記第１の受信コイルと並んだ第２の受信コイルと、位相エンコード方向に前記第１の受信コイルと対向する第３の受信コイルとからなり、前記スキャン手段が、リダクションファクターＲ＝２とすることを特徴とするＭＲＩ装置を提供する。
上記第１２の観点によるＭＲＩ装置では、前記第４の観点によるパラレルイメージング方法を好適に実施できる。
第１３の観点では、本発明は、前記第１２の観点によるＭＲＩ装置において、ｉチャンネルの受信コイルの画像上の点の画素値をａiとし、ｉチャンネルの受信コイルの画像上の点で折返しにより重なっている元の２点のｉチャンネルの受信コイルに対応する各感度係数をｓi1，ｓi2とし、ｋチャンネルの受信コイルの感度係数ｓk1，ｓk2を第１行とすると共にｍチャンネルの受信コイルの感度係数ｓm1，ｓm2を第２行とする行列式の２乗を|Ｓkm|2とし、ｋチャンネルの受信コイルの感度係数ｓk1，ｓk2を第１行とし、ｍチャンネルの受信コイルの感度係数ｓm1，ｓm2を第２行とする正方行列をＳkmとするとき、前記演算手段が、
Figure 2008017958
により表される演算を行うことを特徴とするＭＲＩ装置を提供する。
上記第１３の観点によるＭＲＩ装置では、前記第５の観点によるパラレルイメージング方法を好適に実施できる。
第１４の観点では、本発明は、前記第９から前記第１１のいずれかの観点によるＭＲＩ装置において、前記受信コイルが、第１の受信コイルと、周波数エンコード方向に前記第１の受信コイルと並んだ第２の受信コイルと、位相エンコード方向に前記第１の受信コイルと対向する第３の受信コイルと、周波数エンコード方向に前記第３の受信コイルと並んだ第４の受信コイルとからなり、前記スキャン手段が、リダクションファクターＲ＝２とすることを特徴とするＭＲＩ装置を提供する。
上記第１４の観点によるＭＲＩ装置では、前記第６の観点によるパラレルイメージング方法を好適に実施できる。
第１５の観点では、本発明は、前記第１４の観点によるＭＲＩ装置において、ｉチャンネルの受信コイルの画像上の点の画素値をａiとし、ｉチャンネルの受信コイルの画像上の点で折返しにより重なっている元の２点のｉチャンネルの受信コイルに対応する各感度係数をｓi1，ｓi2とし、ｋチャンネルの受信コイルの感度係数ｓk1，ｓk2を第１行とすると共にｍチャンネルの受信コイルの感度係数ｓm1，ｓm2を第２行とする行列式の２乗を|Ｓkm|2とし、ｋチャンネルの受信コイルの感度係数ｓk1，ｓk2を第１行とし、ｍチャンネルの受信コイルの感度係数ｓm1，ｓm2を第２行とする正方行列をＳkmとするとき、前記演算手段が、
Figure 2008017958
により表される演算を行うことを特徴とするＭＲＩ装置を提供する。
上記第１５の観点によるＭＲＩ装置では、前記第７の観点によるパラレルイメージング方法を好適に実施できる。
第１６の観点では、本発明は、前記第９から前記第１１のいずれかの観点によるＭＲＩ装置において、前記組合せ選択手段が選んだ画像の組合せに対応する受信コイルを第１チャンネルから第Ｌチャンネルまでとし、リダクションファクターをＲとするとき、第１チャンネルから第Ｌチャンネルまでの中からＲチャンネルを選んだ各組合せを一般的に第1Rチャネルから第RRチャネルとし、ｉチャンネルの受信コイルの画像上の点の画素値をａiとし、ｉチャンネルの受信コイルの画像上の点で折返しにより重なっている元のＲ点のｉチャンネルの受信コイルに対応する各感度係数をｓi1，…，ｓiRとし、ｋチャンネルの受信コイルの感度係数ｓk1，…，ｓkRを各行とする正方行列の行列式の２乗を|Ｓ1R-RR|2とし、ｋチャンネルの受信コイルの感度係数ｓk1，…，ｓkRを各行とする正方行列をＳ1R-RRとするとき、前記演算手段が、
Figure 2008017958
により表される演算を行うことを特徴とするＭＲＩ装置を提供する。
上記第１６の観点によるＭＲＩ装置では、前記第８の観点によるパラレルイメージング方法を好適に実施できる。
本発明のパラレルイメージング方法およびＭＲＩ装置によれば、３個以上の受信コイルを用いたパラレルイメージングにおいて、少ない演算量かつ短い処理時間で折返しを元に戻した合成画像を得ることが出来る。
以下、図に示す実施の形態により本発明をさらに詳細に説明する。なお、これにより本発明が限定されるものではない。
図１は、実施例１にかかるＭＲＩ装置１００を示すブロック図である。
このＭＲＩ装置１００において、マグネットアセンブリ１０１は、内部に被検体を挿入するための空間部分（ボア）を有し、この空間部分を取りまくようにして、被検体に一定の静磁場を印加する静磁場コイル１０１Ｃと、Ｘ軸，Ｙ軸，Ｚ軸の勾配磁場を発生するための勾配コイル１０１Ｇと、被検体内の原子核のスピンを励起するためのＲＦパルスを与える送信コイル１０１Ｔと、被検体からのＮＭＲ信号を受信するためのボディコイル１０１(0)およびＩ（≧３）チャンネルの受信コイル１０１(1)，…，１０１(Ｉ)とが配置されている。
静磁場コイル１０１Ｃ，勾配コイル１０１Ｇ，送信コイル１０１Ｔは、それぞれ静磁場電源１０２，勾配コイル駆動回路１０３，ＲＦ電力増幅器１０４に接続されている。また、ボディコイル１０１(0)，受信コイル１０１(1)，…，１０１(Ｉ)は、それぞれ前置増幅器１０５(0)，１０５(1)，…，１０５(Ｉ)に接続されている。
なお、ボディコイル１０１(0)を送信コイル１０１Ｔとして使用する場合もある。
また、静磁場コイル１０１Ｃの代わりに永久磁石を用いてもよい。
シーケンス記憶回路１０８は、計算機１０７からの指令に従い、記憶しているパルスシーケンスに基づいて勾配コイル駆動回路１０３を操作し、勾配コイル１０１Ｇから勾配磁場を発生させると共に、ゲート変調回路１０９を操作し、ＲＦ発振回路１１０の搬送波出力信号を所定タイミング・所定包絡線形状・所定位相のパルス状信号に変調し、それをＲＦパルスとしてＲＦ電力増幅器１０４に加え、ＲＦ電力増幅器１０４でパワー増幅した後、送信コイル１０１Ｔに印加する。
セレクタ１１１は、ボディコイル１０１(0)，受信コイル１０１(1)，…，１０１(Ｉ)で受信され前置増幅器１０５(0)，１０５(1)，…，１０５(Ｉ)で増幅されたＮＭＲ信号をｍ個のレシーバ１１２(1)，１１２(2)，…，１１２(m)に伝達する。これは、ボディコイル１０１(0)，受信コイル１０１(1)，…，１０１(Ｉ)とレシーバ１１２(1)，１１２(2)，…，１１２(m)の対応を可変にするためである。
レシーバ１１２(1)，１１２(2)，…，１１２(m)は、ＮＭＲ信号をデジタル信号に変換し、計算機１０７に入力する。
計算機１０７は、レシーバ１１２からデジタル信号を読み込み、処理を施して、ＭＲ画像を生成する。また、計算機１０７は、操作卓１１３から入力された情報を受け取るなどの全体的な制御を受け持つ。
表示装置１０６は、画像やメッセージを表示する。
図２は、実施例１に係るボディコイル１０１(0)およびＩ＝３の場合の受信コイル１０１(1)，１０１(2)，１０１(3)の配置例を示す概念図である。
第１チャンネルＣＨ１の受信コイル受信コイル１０１(1)は被検体Ｈの上の左寄りに置かれ、第２チャンネルＣＨ２の受信コイル受信コイル１０１(2)は被検体Ｈの上の右寄りに置かれ、第３チャンネルＣＨ３の受信コイル受信コイル１０１(3)は被検体Ｈの下に置かれている。
被検体Ｈの上下方向が位相エンコード方向Ｐであり、被検体Ｈの左右方向が周波数エンコード方向Ｆである。
図３は、実施例１に係るパラレルイメージング処理を示すフロー図である。
ステップＱ１では、ボディコイル１０１(0)を用い且つ位相エンコードステップを間引かないスキャンを行って校正用データを収集する。また、受信コイル１０１(1)，１０１(2)，１０１(3)を用い且つリダクションファクターＲ＝２で位相エンコードステップを間引いたスキャンを行って校正用データを収集する。そして、それら校正用データを基に、受信コイル１０１(1)の感度係数ｓ11，ｓ12と、受信コイル１０１(2)の感度係数ｓ21，ｓ22と、受信コイル１０１(3)の感度係数ｓ31，ｓ32とを求める。
ステップＱ２では、受信コイル１０１(1)，１０１(2)，１０１(3)を用い且つリダクションファクターＲ＝２で位相エンコードステップを間引いたスキャンを行ってデータを収集する。そして、それらデータを基に、受信コイル１０１(1)の画像ａ1と、受信コイル１０１(2)の画像ａ2と、受信コイル１０１(3)の画像ａ3とを求める。
ステップＱ３では、
Figure 2008017958
により、合成画像Ｖを生成する。そして、処理を終了する。
実施例１のＭＲＩ装置１００によれば、従来方法に比べて演算量を２／３に低減でき、処理時間も短縮できる。
図４は、実施例２に係るボディコイル１０１(0)およびＩ＝３の場合の受信コイル１０１(1)，１０１(2)，１０１(3)，１０１(4)の配置例を示す概念図である。
第１チャンネルＣＨ１の受信コイル受信コイル１０１(1)は被検体Ｈの上の左寄りに置かれ、第２チャンネルＣＨ２の受信コイル受信コイル１０１(2)は被検体Ｈの上の右寄りに置かれ、第３チャンネルＣＨ３の受信コイル受信コイル１０１(3)は被検体Ｈの下の左寄りに置かれ、第４チャンネルＣＨ４の受信コイル受信コイル１０１(4)は被検体Ｈの下の右寄りに置かれている。
被検体Ｈの上下方向が位相エンコード方向Ｐであり、被検体Ｈの左右方向が周波数エンコード方向Ｆである。
図５は、実施例２に係るパラレルイメージング処理を示すフロー図である。
ステップＴ１では、ボディコイル１０１(0)を用い且つ位相エンコードステップを間引かないスキャンを行って校正用データを収集する。また、受信コイル１０１(1)，１０１(2)，１０１(3)，１０１(4)を用い且つリダクションファクターＲ＝２で位相エンコードステップを間引いたスキャンを行って校正用データを収集する。そして、それら校正用データを基に、受信コイル１０１(1)の感度係数ｓ11，ｓ12と、受信コイル１０１(2)の感度係数ｓ21，ｓ22と、受信コイル１０１(3)の感度係数ｓ31，ｓ32と、受信コイル１０１(4)の感度係数ｓ41，ｓ42とを求める。
ステップＴ２では、受信コイル１０１(1)，１０１(2)，１０１(3)，１０１(4)を用い且つリダクションファクターＲ＝２で位相エンコードステップを間引いたスキャンを行ってデータを収集する。そして、それらデータを基に、受信コイル１０１(1)の画像ａ1と、受信コイル１０１(2)の画像ａ2と、受信コイル１０１(3)の画像ａ3と、受信コイル１０１(4)の画像ａ4とを求める。
ステップＴ３では、
Figure 2008017958
により、合成画像Ｖを生成する。そして、処理を終了する。
実施例２のＭＲＩ装置によれば、従来方法に比べて演算量を２／３に低減でき、処理時間も短縮できる。
図６は、実施例３に係るパラレルイメージング処理を示すフロー図である。
ステップＧ１では、ボディコイル１０１(0)を用い且つ位相エンコードステップを間引かないスキャンを行って校正用データを収集する。また、受信コイル１０１(1)，…，１０１(I)を用い且つリダクションファクターＲで位相エンコードステップを間引いたスキャンを行って校正用データを収集する。そして、それら校正用データを基に、ｉ＝１，…，Ｉとするとき、受信コイル１０１(i)の感度係数ｓi1，…，ｓiRを求める。
ステップＧ２では、受信コイル１０１(1)，…，１０１(I)を用い且つリダクションファクターＲで位相エンコードステップを間引いたスキャンを行ってデータを収集する。そして、それらデータを基に、受信コイル１０１(i)の画像ａiを求める。
ステップＧ３では、受信コイル１０１(1)，…，１０１(I)の中からＲ個を選び、組合せを作る。作られる組合せの数は、組合せ数IＣR個になる。次いで、各組合せに属する受信コイルの感度係数の正方行列の行列式の２乗を計算する。例えば、ある組合せに受信コイル１０１(1R)，…，１０１(RR)が属するなら、次式のような感度係数の正方行列の行列式の２乗｜Ｓ1R-RR｜2を計算する。
Figure 2008017958
そして、値を比較し、例えば値の平均値の２０％以上の値を持つ組合せだけを選ぶ。
ステップＧ４では、選んだ組合せだけを用いて、
Figure 2008017958
により、合成画像Ｖを生成する。そして、処理を終了する。
実施例３のＭＲＩ装置によれば、従来方法に比べて演算量を低減でき、処理時間も短縮できる。
本発明のパラレルイメージング方法及びＭＲＩ装置は、高速ＭＲＩに利用できる。
実施例１に係るＭＲＩ装置を示す構成ブロック図である。 実施例１に係る受信コイルの配置を示す概念図である。 実施例１に係るパラレルイメージング処理を示すフロー図である。 実施例２に係る受信コイルの配置を示す概念図である。 実施例２に係るパラレルイメージング処理を示すフロー図である。 実施例３に係るパラレルイメージング処理を示すフロー図である。
符号の説明
８ シーケンス記憶回路
１００ ＭＲＩ装置
１０７ 計算機
１０１Ｔ 送信コイル
１０１(0) ボディコイル
１０１(1)〜１０１(Ｉ) 受信コイル

Claims (16)

  1. ３個以上の受信コイルを用い且つ位相エンコードステップを間引いたスキャンを行って前記各受信コイルのデータを収集し、前記各データから画像をそれぞれ生成し、前記各画像の組合せの中から折返しを元に戻す演算に用いる画像の組合せを選び、選んだ画像の組合せ及び前記各受信コイルの感度分布を利用して前記演算により合成画像を得ることを特徴とするパラレルイメージング方法。
  2. 請求項１に記載のパラレルイメージング方法において、前記受信コイルの各組合せに対応する受信コイルの感度係数の正方行列の行列式の２乗に基づいて、前記各画像の組合せの中から折返しを元に戻す演算に用いる画像の組合せを選ぶことを特徴とするパラレルイメージング方法。
  3. 請求項１に記載のパラレルイメージング方法において、前記受信コイルの配置および前記フェーズエンコード方向に基づいて、前記各画像の中から折返しを元に戻す演算に用いる画像の組合せを選ぶことを特徴とするパラレルイメージング方法。
  4. 請求項１から請求項３のいずれかに記載のパラレルイメージング方法において、前記受信コイルが、第１の受信コイルと、周波数エンコード方向に前記第１の受信コイルと並んだ第２の受信コイルと、位相エンコード方向に前記第１の受信コイルと対向する第３の受信コイルとからなり、リダクションファクターＲ＝２とすることを特徴とするパラレルイメージング方法。
  5. 請求項４に記載のパラレルイメージング方法において、ｉチャンネルの受信コイルの画像上の点の画素値をａiとし、ｉチャンネルの受信コイルの画像上の点で折返しにより重なっている元の２点のｉチャンネルの受信コイルに対応する各感度係数をｓi1，ｓi2とし、ｋチャンネルの受信コイルの感度係数ｓk1，ｓk2を第１行とすると共にｍチャンネルの受信コイルの感度係数ｓm1，ｓm2を第２行とする行列式の２乗を|Ｓkm|2とし、ｋチャンネルの受信コイルの感度係数ｓk1，ｓk2を第１行とし、ｍチャンネルの受信コイルの感度係数ｓm1，ｓm2を第２行とする正方行列をＳkmとするとき、前記演算が、
    Figure 2008017958
    により表されることを特徴とするパラレルイメージング方法。
  6. 請求項１から請求項３のいずれかに記載のパラレルイメージング方法において、前記受信コイルが、第１の受信コイルと、周波数エンコード方向に前記第１の受信コイルと並んだ第２の受信コイルと、位相エンコード方向に前記第１の受信コイルと対向する第３の受信コイルと、周波数エンコード方向に前記第３の受信コイルと並んだ第４の受信コイルとからなり、リダクションファクターＲ＝２とすることを特徴とするパラレルイメージング方法。
  7. 請求項６に記載のパラレルイメージング方法において、ｉチャンネルの受信コイルの画像上の点の画素値をａiとし、ｉチャンネルの受信コイルの画像上の点で折返しにより重なっている元の２点のｉチャンネルの受信コイルに対応する各感度係数をｓi1，ｓi2とし、ｋチャンネルの受信コイルの感度係数ｓk1，ｓk2を第１行とすると共にｍチャンネルの受信コイルの感度係数ｓm1，ｓm2を第２行とする行列式の２乗を|Ｓkm|2とし、ｋチャンネルの受信コイルの感度係数ｓk1，ｓk2を第１行とし、ｍチャンネルの受信コイルの感度係数ｓm1，ｓm2を第２行とする正方行列をＳkmとするとき、前記演算が、
    Figure 2008017958
    により表されることを特徴とするパラレルイメージング方法。
  8. 請求項１から請求項３のいずれかに記載のパラレルイメージング方法において、前記選んだ画像の組合せに対応する受信コイルを第１チャンネルから第Ｌチャンネルまでとし、リダクションファクターをＲとするとき、第１チャンネルから第Ｌチャンネルまでの中からＲチャンネルを選んだ各組合せを一般的に第1Rチャネルから第RRチャネルとし、ｉチャンネルの受信コイルの画像上の点の画素値をａiとし、ｉチャンネルの受信コイルの画像上の点で折返しにより重なっている元のＲ点のｉチャンネルの受信コイルに対応する各感度係数をｓi1，…，ｓiRとし、ｋチャンネルの受信コイルの感度係数ｓk1，…，ｓkRを各行とする正方行列の行列式の２乗を|Ｓ1R-RR|2とし、ｋチャンネルの受信コイルの感度係数ｓk1，…，ｓkRを各行とする正方行列をＳ1R-RRとするとき、前記演算が、
    Figure 2008017958
    により表されることを特徴とするパラレルイメージング方法。
  9. ３個以上の受信コイルと、前記３個以上の受信コイルを用い且つ位相エンコードステップを間引いたスキャンを行って前記各受信コイルのデータを収集するスキャン手段と、前記各データから画像をそれぞれ生成する画像生成手段と、前記各画像の組合せの中から折返しを元に戻す演算に用いる画像の組合せを選ぶ組合せ選択手段と、選んだ画像の組合せ及び前記各受信コイルの感度分布を利用して前記演算により合成画像を得る演算手段とを具備したことを特徴とするＭＲＩ装置。
  10. 請求項９に記載のＭＲＩ装置において、前記組合せ選択手段は、前記受信コイルの各組合せに対応する受信コイルの感度係数の正方行列の行列式の２乗に基づいて、前記各画像の組合せの中から折返しを元に戻す演算に用いる画像の組合せを選ぶことを特徴とするＭＲＩ装置。
  11. 請求項９に記載のＭＲＩ装置において、前記組合せ選択手段は、前記受信コイルの配置および前記フェーズエンコード方向に基づいて、前記各画像の中から折返しを元に戻す演算に用いる画像の組合せを選ぶことを特徴とするＭＲＩ装置。
  12. 請求項９から請求項１１のいずれかに記載のＭＲＩ装置において、前記受信コイルが、第１の受信コイルと、周波数エンコード方向に前記第１の受信コイルと並んだ第２の受信コイルと、位相エンコード方向に前記第１の受信コイルと対向する第３の受信コイルとからなり、前記スキャン手段が、リダクションファクターＲ＝２とすることを特徴とするＭＲＩ装置。
  13. 請求項１２に記載のＭＲＩ装置において、ｉチャンネルの受信コイルの画像上の点の画素値をａiとし、ｉチャンネルの受信コイルの画像上の点で折返しにより重なっている元の２点のｉチャンネルの受信コイルに対応する各感度係数をｓi1，ｓi2とし、ｋチャンネルの受信コイルの感度係数ｓk1，ｓk2を第１行とすると共にｍチャンネルの受信コイルの感度係数ｓm1，ｓm2を第２行とする行列式の２乗を|Ｓkm|2とし、ｋチャンネルの受信コイルの感度係数ｓk1，ｓk2を第１行とし、ｍチャンネルの受信コイルの感度係数ｓm1，ｓm2を第２行とする正方行列をＳkmとするとき、前記演算手段が、
    Figure 2008017958
    により表される演算を行うことを特徴とするＭＲＩ装置。
  14. 請求項９から請求項１１のいずれかに記載のＭＲＩ装置において、前記受信コイルが、第１の受信コイルと、周波数エンコード方向に前記第１の受信コイルと並んだ第２の受信コイルと、位相エンコード方向に前記第１の受信コイルと対向する第３の受信コイルと、周波数エンコード方向に前記第３の受信コイルと並んだ第４の受信コイルとからなり、前記スキャン手段が、リダクションファクターＲ＝２とすることを特徴とするＭＲＩ装置。
  15. 請求項１４に記載のＭＲＩ装置において、ｉチャンネルの受信コイルの画像上の点の画素値をａiとし、ｉチャンネルの受信コイルの画像上の点で折返しにより重なっている元の２点のｉチャンネルの受信コイルに対応する各感度係数をｓi1，ｓi2とし、ｋチャンネルの受信コイルの感度係数ｓk1，ｓk2を第１行とすると共にｍチャンネルの受信コイルの感度係数ｓm1，ｓm2を第２行とする行列式の２乗を|Ｓkm|2とし、ｋチャンネルの受信コイルの感度係数ｓk1，ｓk2を第１行とし、ｍチャンネルの受信コイルの感度係数ｓm1，ｓm2を第２行とする正方行列をＳkmとするとき、前記演算手段が、
    Figure 2008017958
    により表される演算を行うことを特徴とするＭＲＩ装置。
  16. 請求項９から請求項１１のいずれかに記載のＭＲＩ装置において、前記組合せ選択手段が選んだ画像の組合せに対応する受信コイルを第１チャンネルから第Ｌチャンネルまでとし、リダクションファクターをＲとするとき、第１チャンネルから第Ｌチャンネルまでの中からＲチャンネルを選んだ各組合せを一般的に第1Rチャネルから第RRチャネルとし、ｉチャンネルの受信コイルの画像上の点の画素値をａiとし、ｉチャンネルの受信コイルの画像上の点で折返しにより重なっている元のＲ点のｉチャンネルの受信コイルに対応する各感度係数をｓi1，…，ｓiRとし、ｋチャンネルの受信コイルの感度係数ｓk1，…，ｓkRを各行とする正方行列の行列式の２乗を|Ｓ1R-RR|2とし、ｋチャンネルの受信コイルの感度係数ｓk1，…，ｓkRを各行とする正方行列をＳ1R-RRとするとき、前記演算手段が、
    Figure 2008017958
    により表される演算を行うことを特徴とするＭＲＩ装置。
JP2006190999A 2006-07-12 2006-07-12 Ｍｒｉ装置 Active JP4464373B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2006190999A JP4464373B2 (ja) 2006-07-12 2006-07-12 Ｍｒｉ装置
US11/776,115 US8054078B2 (en) 2006-07-12 2007-07-11 Parallel imaging method and MRI apparatus
CN2007101287812A CN101105522B (zh) 2006-07-12 2007-07-12 并行成像方法和mri设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2006190999A JP4464373B2 (ja) 2006-07-12 2006-07-12 Ｍｒｉ装置

Publications (2)

Publication Number Publication Date
JP2008017958A true JP2008017958A (ja) 2008-01-31
JP4464373B2 JP4464373B2 (ja) 2010-05-19

Family

ID=38949306

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006190999A Active JP4464373B2 (ja) 2006-07-12 2006-07-12 Ｍｒｉ装置

Country Status (3)

Country Link
US (1) US8054078B2 (ja)
JP (1) JP4464373B2 (ja)
CN (1) CN101105522B (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105182263A (zh) * 2008-04-28 2015-12-23 康奈尔大学 分子mri中的磁敏度精确量化

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9864032B2 (en) * 2010-01-05 2018-01-09 National Health Research Institutes Magnetic resonance imaging system
US9396562B2 (en) * 2012-09-26 2016-07-19 Siemens Aktiengesellschaft MRI reconstruction with incoherent sampling and redundant haar wavelets
US10107884B2 (en) 2013-03-13 2018-10-23 Koninklijke Philips N.V. Automatic optimization of parallel imaging acceleration parameters
DE102014204706B4 (de) 2014-03-13 2019-06-06 Siemens Healthcare Gmbh Empfänger-Baugruppe eines bildgebenden Magnetresonanz-Systems und ein bildgebendes Magnetresonanz-System
US11047935B2 (en) 2015-05-14 2021-06-29 Ohio State Innovation Foundation Systems and methods for estimating complex B1+ fields of transmit coils of a magnetic resonance imaging (MRI) system
WO2018136705A1 (en) 2017-01-19 2018-07-26 Ohio State Innovation Foundation Estimating absolute phase of radio frequency fields of transmit and receive coils in a magnetic resonance
JP6820876B2 (ja) * 2018-03-20 2021-01-27 株式会社日立製作所 磁気共鳴撮影装置、感度分布算出方法、及び、感度分布算出プログラム

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4751643A (en) * 1986-08-04 1988-06-14 General Electric Company Method and apparatus for determining connected substructures within a body
US4984157A (en) * 1988-09-21 1991-01-08 General Electric Company System and method for displaying oblique planar cross sections of a solid body using tri-linear interpolation to determine pixel position dataes
US5332968A (en) * 1992-04-21 1994-07-26 University Of South Florida Magnetic resonance imaging color composites
JP3149097B2 (ja) 1992-08-28 2001-03-26 カシオ計算機株式会社 音響成分抽出装置、それを用いた電子楽器、及び周波数成分抽出装置
JP2641417B2 (ja) 1996-05-09 1997-08-13 安川商事株式会社 時空間微分法を用いた計測装置
JP4380867B2 (ja) 1999-12-28 2009-12-09 株式会社東芝 磁気共鳴スペクトロスコピー装置
US7006085B1 (en) * 2000-10-30 2006-02-28 Magic Earth, Inc. System and method for analyzing and imaging three-dimensional volume data sets
US6841998B1 (en) * 2001-04-06 2005-01-11 Mark Griswold Magnetic resonance imaging method and apparatus employing partial parallel acquisition, wherein each coil produces a complete k-space datasheet
JP4443079B2 (ja) 2001-09-13 2010-03-31 株式会社日立メディコ 磁気共鳴イメージング装置及び磁気共鳴イメージング装置用ｒｆ受信コイル
US6950543B2 (en) * 2002-05-02 2005-09-27 Ge Medical Systems Global Technology Company, Llc Method and system for image reconstruction
JP4607430B2 (ja) 2003-03-28 2011-01-05 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｒｆコイルおよび磁気共鳴撮影装置
JP4391214B2 (ja) 2003-12-12 2009-12-24 株式会社日立メディコ 磁気共鳴イメージング装置
JP4037385B2 (ja) * 2004-04-23 2008-01-23 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 画像生成方法およびｍｒｉ装置
US7511495B2 (en) * 2005-04-25 2009-03-31 University Of Utah Systems and methods for image reconstruction of sensitivity encoded MRI data
EP2069815A1 (en) * 2006-09-13 2009-06-17 Eidgenössische Technische Hochschule (ETH) Method and apparatus for mri coil array compression
US8219176B2 (en) * 2007-03-08 2012-07-10 Allegheny-Singer Research Institute Single coil parallel imaging
US20090285463A1 (en) * 2008-04-18 2009-11-19 Ricardo Otazo Superresolution parallel magnetic resonance imaging
US7782058B2 (en) * 2008-04-28 2010-08-24 General Electric Company System and method for accelerated MR imaging
US8354844B2 (en) * 2009-12-23 2013-01-15 Universitaetsklinikum Freiburg Method for data acquisition acceleration in magnetic resonance imaging (MRI) with N-dimensional spatial encoding using two or more receiver coil arrays and non-linear phase distributions

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105182263A (zh) * 2008-04-28 2015-12-23 康奈尔大学 分子mri中的磁敏度精确量化

Also Published As

Publication number Publication date
CN101105522A (zh) 2008-01-16
US20080013811A1 (en) 2008-01-17
CN101105522B (zh) 2011-08-10
JP4464373B2 (ja) 2010-05-19
US8054078B2 (en) 2011-11-08

Similar Documents

Publication Publication Date Title
KR100742458B1 (ko) Ｍｒ 화상 생성 방법 및 ｍｒｉ 장치
JP4464373B2 (ja) Ｍｒｉ装置
US7288937B2 (en) Image generating method and MRI apparatus
US6703834B2 (en) Systems and methods for correcting phase errors
US7187791B2 (en) Coil sensitivity map generating method, parallel imaging method, and MRI apparatus
JP2014083445A (ja) 磁気共鳴イメージングシステム及び磁気共鳴イメージング方法
JP3891799B2 (ja) Ｍｒｉ装置
US7538550B2 (en) Magnetic resonance imaging apparatus, magnetic resonance imaging method, scan apparatus, program and storage medium
JP2004275198A (ja) 核磁気共鳴イメージング装置、および核磁気共鳴イメージング方法
JP3117670B2 (ja) マルチスライスｍｒイメージング方法およびｍｒｉ装置
JP2001299720A (ja) 磁場ドリフト補正用データ収集方法、磁場ドリフト補正方法およびｍｒｉ装置
US7486072B2 (en) Parallel multiple coil MR scanning method and MRI apparatus with flow compensation in a calibration scan including a calibration of a coil sensitivity of each coil
US20080246474A1 (en) Mri apparatus and control method therefor
US6470203B2 (en) MR imaging method, phase error measuring method and MRI system
US20200088822A1 (en) Parallel multi-slice mr imaging
JPWO2006112497A1 (ja) 磁気共鳴イメージング装置及び方法
US10073152B2 (en) Magnetic resonance imaging device and magnetic resonance imaging method using the same device
JP6867256B2 (ja) 磁気共鳴撮像装置及び画像処理方法
JP3720752B2 (ja) ０次位相検出方法およびｍｒｉ装置
JP2015047240A (ja) 磁気共鳴装置およびプログラム
JP6850224B2 (ja) 磁気共鳴イメージング装置及び信号抑制方法
JPH06205754A (ja) Ｍｒ装置のプリパレーションパルス印加方法
JP4607928B2 (ja) 画像生成方法およびｍｒｉ装置
JP4678926B2 (ja) Ｍｒｉ装置
JP2022178683A (ja) 磁気共鳴撮影装置及び感度分布算出プロブラム

Legal Events

Date Code Title Description
A625 Written request for application examination (by other person)

Free format text: JAPANESE INTERMEDIATE CODE: A625

Effective date: 20080516

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20081027

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20081104

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090202

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090407

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090703

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090908

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090916

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

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20100218

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130226

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Ref document number: 4464373

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

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130226

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140226

Year of fee payment: 4

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250