JP2016509940A - パラレルイメージング加速パラメータの自動最適化 - Google Patents

パラレルイメージング加速パラメータの自動最適化 Download PDF

Info

Publication number
JP2016509940A
JP2016509940A JP2015562491A JP2015562491A JP2016509940A JP 2016509940 A JP2016509940 A JP 2016509940A JP 2015562491 A JP2015562491 A JP 2015562491A JP 2015562491 A JP2015562491 A JP 2015562491A JP 2016509940 A JP2016509940 A JP 2016509940A
Authority
JP
Japan
Prior art keywords
reduction
scan
reduction factor
coil
snr
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
JP2015562491A
Other languages
English (en)
Other versions
JP6430417B2 (ja
JP2016509940A5 (ja
Inventor
ペーター ボルネート
ペーター ボルネート
ミシェル ポール ジュリアーン ジュリッセン
ミシェル ポール ジュリアーン ジュリッセン
マリヤ イワノワ ドネヴァ
マリヤ イワノワ ドネヴァ
アドリアヌス ジョセフ ウィリブロルドゥス ドゥアインダム
アドリアヌス ジョセフ ウィリブロルドゥス ドゥアインダム
ケイ ネルケ
ケイ ネルケ
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of JP2016509940A publication Critical patent/JP2016509940A/ja
Publication of JP2016509940A5 publication Critical patent/JP2016509940A5/ja
Application granted granted Critical
Publication of JP6430417B2 publication Critical patent/JP6430417B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5611Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/32Excitation or detection systems, e.g. using radio frequency signals
    • G01R33/34Constructional details, e.g. resonators, specially adapted to MR
    • G01R33/34092RF coils specially adapted for NMR spectrometers
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5608Data 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)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Signal Processing (AREA)
  • General Health & Medical Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

パラレル磁気共鳴映像システム1は、複数のコイル素子を有する少なくとも1個の無線周波(RF)コイル10、12と、スマートセレクトユニット24と、パラレルイメージングパラメータユニット28と、シーケンスコントロール16とを含む。スマートセレクトユニット24は、少なくとも1個のRFコイルによる対象者のプリスキャン又は事前スキャンから、様々な組の低減係数に基づいて信号マップ及び複数の雑音マップを構築する60。パラレルイメージングパラメータユニット28は、最も高い信号対雑音比(SNR)を有する雑音マップに対応する1組の低減係数を選択する。シーケンスコントロール16は、選択された低減係数に基づいて対象者の磁気共鳴映像スキャンを行う。

Description

以下の内容は、全体として医用撮像に関する。以下の内容は、磁気共鳴イメージングに関連して特定の応用例を見出し、磁気共鳴イメージングに特に関連して説明される。しかし、以下の内容は他の使用シナリオにも応用され、必ずしも上記の応用例に限定されないことが理解されよう。
パラレルイメージングを用いる3次元(3D)磁気共鳴(MR)スキャン、又はスライス方向へのパラレルイメージングを用いるマルチスライススキャンでは、複数のコイル素子及び/又はコイルアレイを有する局所コイル及び/又は全身コイルが使用される。例えば全身コイルは、同時に動作する8個、16個、32個等のコイル素子を含むことができる。コイル素子は、画像を再構築するために使用される重複空間情報をもたらす。パラレルイメージングの性能は、スキャンのジオメトリ及びスキャンに使用されるコイルの設計に依存する方向によって異なる。
オペレータは、スキャンごとにスライス方向、位相エンコード方向、周波数エンコード方向等、1つ又は複数のエンコード方向に関する加速係数を選択することができる。加速係数は、k空間を完全にサンプリングするのに必要な共鳴励起及びサンプリングサイクルを減らすことによってスキャン時間を短縮する。加速量は多くの場合低減係数と呼ばれ、R=1は低減なしを意味し、R=2は2倍の低減等を意味する。エンコード方向は、例えば1つ又は複数の位相エンコード方向、読取方向、又はスライス方向を含み得る。低減係数は、小数又は整数とすることができる。例えば低減係数は、1、8/7、4/3、1.5、8/5、2、2.5、8/3、3、3.5、4等を含み得る。厳密でない又は正しくない選択は低い信号対雑音比(SNR)の原因となり、低いSNRは低い画質を引き起こす。こうした低減係数の厳密でない又は正しくない選択は、スキャンが完了し、その結果生じる画像が見られた後でしかオペレータによって認識されない。更に、多数の個々の加速係数は複雑なセットアップがスキャンごとに繰り返されることを引き起こし、このことはシステムを使いにくくし、オペレータのミスを生まれやすくする。
他の手法は、所与の関心領域のMRスキャンを行うための1組の「最適な」コイル素子を選択することを含む。1組の最適なコイル素子を選択する過程は、雑音情報を与え得る調査スキャン又はプリスキャンを含む。或いは、対象者の事前スキャンが雑音情報を与え得る。
以下の内容は、上記の問題及び他の問題に対処するパラレル加速パラメータの新規且つ改善された自動最適化を開示する。
一態様によれば、パラレル磁気共鳴イメージングシステムが、複数のコイル素子を有する少なくとも1個の無線周波(RF)コイルと、スマートセレクトユニットと、パラレルイメージングパラメータユニットと、シーケンスコントロールとを含む。スマートセレクトユニットは、少なくとも1個のRFコイルによる対象者のプリスキャン又は事前スキャンから、1つ又は複数の低減方向に沿った様々な組の低減係数に基づいて信号マップ及び複数の雑音マップを構築する。パラレルイメージングパラメータユニットは、信号対雑音比(SNR)及びスキャン速度を最適化する、雑音マップに対応する1組の低減係数を選択する。シーケンスコントロールは、選択された低減係数に基づいて対象者の磁気共鳴イメージングスキャンを行う。
別の態様によれば、パラレル磁気共鳴イメージングの方法が、対象者のプリスキャンから、受け取られた全体的な低減係数に基づいて信号マップ及び複数の雑音マップを構築するステップを含み、各雑音マップは1つ又は複数の低減方向に沿った様々な組の低減係数に基づいて算出される。信号対雑音比(SNR)基準及びスキャン速度基準を満たす雑音マップに対応する低減係数が選択される。選択された低減係数に基づき、対象者の磁気共鳴イメージングスキャンが行われる。
別の態様によれば、パラレル磁気共鳴イメージングシステムが、複数のコイル素子を有する少なくとも1個のRFコイルと、ユーザインターフェイスと、1個又は複数個のプロセッサとを含む。ユーザインターフェイスは、画像データを生成するための最低限の画質及び/又は最大時間の少なくとも1つを受け取る。1個又は複数個のプロセッサは、少なくとも1個のRFコイルによる対象者の事前スキャンから、1つ又は複数の低減方向に沿って複数の低減係数のそれぞれについて信号マップ及び複数の雑音マップを構築するように構成される。1個又は複数個のプロセッサは、少なくとも受け取られた最低限の画質で、及び/又は受け取られた最大時間内に画像データを生成する低減係数を選択するように更に構成され、選択された低減係数に基づいて対象者の磁気共鳴イメージングスキャンをイメージングシステムに行わせる。
1つの利点は、選択されるパラメータの数が減ることである。
別の利点は、全体的な低減係数や最小信号対雑音比(SNR)等の定量化可能なパラメータに基づき、最適化されたエンコード方向及び/又は低減係数が自動で選択されることにある。
別の利点は、エンコード方向及び/又は低減係数を最適化するために、コイル感度情報及び患者ローディングを含む従来のコイル調査スキャンから得られる情報を使用することにある。
別の利点は、ユーザインターフェイスの単純化にある。
以下の詳細な説明を読んで理解すれば、更なる利点が当業者に理解されよう。
本発明は、様々な構成要素及び構成要素の配置、並びに様々なステップ及びステップの配列の形をとり得る。図面は好ましい実施形態の例示目的に過ぎず、本発明を限定するものとして解釈されるべきではない。
パラレル加速パラメータシステムを自動で最適化する一実施形態を概略的に示す。 スマートセレクト処理の一実施形態を流れ図で示す。 パラレル加速パラメータを自動で最適化する一方法を流れ図で示す。
図1を参照すると、パラレル加速パラメータシステム1を自動で最適化する一実施形態が概略的に示されている。システム1は、水平ボアスキャナ、オープンシステム、cタイプスキャナ等のMRスキャナ2を含む。このスキャナは、イメージングのために対象者4がその内部に配置される検査領域を画定する開口部又はボアを含む。MRスキャナ2は、主磁石6、1個又は複数個の傾斜磁場コイル8、及びマルチコイル素子の全身無線周波(RF)コイル10を含む。主磁石6は、超電導コイル又は抵抗コイルを含むことができ、水平又は垂直の静的な場等の静的なB場を作り出す。傾斜磁場コイルは、対象者の組織内の共鳴を方向的にエンコードし、勾配エコーを発生させるために傾斜磁場を加える。この方向的エンコードは、位相エンコード方向、プレパレーション方向、スライスエンコード方向等を定めることができる。3次元(3D)スキャンでは、傾斜磁場コイルが、大量の2D画像を明示して再構築するのではなく、読取方向に対して直角な2つの直交方向に概して位相エンコードする。スキャナ2は、胴体コイル、頭部コイル、前立腺コイル、脚部コイル等、1個又は複数個の局所RFコイル12を含むことができる。全身RFコイル10及び/又は局所RFコイル12は、感度エンコーディング(SENSE:sensitivity encoding)等のパラレルイメージング技法に基づきパラレルイメージングを行う、複数のコイル素子又はコイルアレイを有する少なくとも1個のコイルを含む。全身RFコイル及び/又は局所RFコイルは、対象者4の組織内の磁気共鳴を引き起こし、操作するRFパルスを発生させる。全身RFコイル及び/又は局所コイルは、磁気共鳴データとして共鳴減衰に由来するRF信号を受け取り、それらのRF信号はRF受信機14に伝えられる。コイル素子又はアレイは、選択された送信及び/又は受信低減係数を与えるために、個々に又は選択可能なグループで送受信することができる。
システム1は、T1スキャン、T2スキャン等のイメージングスキャンや、プリスキャン及び/又はコイル調査スキャンの動作を制御するシーケンスコントローラ16を含む。コイル調査スキャン又はプリスキャンは、例えば勾配エコーや高速フィールドエコー(FFE:fast field echo)シーケンスを含み得る。コイル調査スキャン又はプリスキャンは、k空間の直交系サンプリング又は六角形サンプリングを使用する。一実施形態では、サンプリングがk空間の圧縮パラレル又はセンシングサンプリングを含み得る。シーケンスコントローラ16は、全身RFコイル10及び/又は局所RFコイル12の動作を制御するRF送信機ユニット18、並びに傾斜磁場コイル8の動作を制御するグラジェントコントローラ20を制御する。RF送信機ユニット18は、複数の送信機を使用して等、コイル素子又はアレイを独立に駆動するように構成される。制御ユニットと対応するコイルとの間の通信は無線でも有線でも良い。
システム1は、ユーザインターフェイス22、スマートセレクトユニット24、信号対雑音比(SNR)算出ユニット26、及びパラレルイメージングパラメータユニット28を含む。ユーザインターフェイスは、全体的な低減係数、最小SNR、最大スキャン時間等を受け取り、オペレータによる審査等のために加速係数を表示するために等、医療関係者又はスキャナ2のオペレータと対話するように構成されるワークステーション36内に適切に組み入れられる。ワークステーション36は、電子プロセッサ又は電子処理装置30、表示装置32、及び少なくとも1つの入力装置34を含む。表示装置32は、メニュー、パネル、及びユーザコントロールを表示し、入力装置34は、全体的な低減係数、加速係数、又は最小SNRを含む入力を医療関係者から受け取る。入力装置は、プリスキャン又はコイル調査スキャンを開始するためのコマンドも受け取る。ワークステーション36は、デスクトップコンピュータ、ラップトップ、タブレット、モバイルコンピューティング装置、スマートフォン等とすることができる。表示装置32は、コンピュータモニタ、タッチスクリーン、陰極線管(CRT)、蓄積管、フラットパネルディスプレイ、発光ダイオード(LED)ディスプレイ、エレクトロルミネセントディスプレイ(ELD)、プラズマディスプレイパネル(PDP)、液晶ディスプレイ(LCD)、有機発光ダイオードディスプレイ(OLED)、プロジェクタ等を含み得る。入力装置34は、キーボード、マウス、マイクロフォン等とすることができる。
スマートセレクトユニット24は、RF受信機14からコイル調査データを受け取り、その調査データをメモリ38内に記憶する。メモリは、ランダムアクセスメモリ(RAM)等のコンピュータメモリ、又はハードディスク、ディスクアレイ、ソリッドステートメモリ等のディスク記憶域とすることができる。スマートセレクトユニットは、受け取られた全体的な低減係数及び/又は最小SNRを満たす複数の低減係数を生成する。低減係数は、位相エンコーディングを含む。例えば、4の全体的な低減係数は、(4×1)、(3.5×8/7)、(2.5×8/5)、(2×2)、(1.5×8/3)、(1×4)等の低減係数を含み得る。エンコード方向の低減係数の積は、全体的な低減係数の値に近づく。それらの値は、反復によって更に洗練させることができる粗い格子に基づく離散値である。スマートセレクトユニットは、信号マップ40及び複数の雑音マップ42を構築する。
SNR算出ユニットは、様々な方向において、様々なパラレル低減係数のSNRを計算する。SNRは、信号マップを雑音マップの1つで割ることによって算出され、このことはボクセルごとのSNR値をもたらす。SNRは、

等、ボリュームの平均によって算出され、Nはボクセルの数であり、snr(r)はr番目のボクセルの信号対雑音比(SNR)である。平均SNRを計算し、又は雑音マップから基準を得る他の方法も考えられる。一実施形態では、平均が関心領域内のボクセルだけを含むことができる。
パラレルイメージングパラメータユニット28は、許容可能な最も高いSNRを有する受け取られた全体的な低減係数を満たす低減係数及び低減方向を選択する。低減係数は、最も高いSNRを有する雑音マップを構築するために使用されるエンコード方向のためのものである。パラレルイメージングパラメータユニットは、選択された低減係数前後のより小さい刻み又はより小さい増分離散値で更に洗練させるために、選択された低減係数を反復内でスマートセレクトユニットに供給することができる。別の実施形態では、低減係数が離散的刻みでエンコード方向間に広がり、最も高いSNR前後でより小さい刻みで更に洗練される。
ユーザインターフェイスは、コントラスト画像としてのSNR値、低減係数の組ごとの平均SNR値、及び/又は最も高いSNRを表示することができる。ユーザインターフェイスは、選択された低減係数をシーケンスコントロールに供給し、選択された低減係数及び/又はエンコード方向に基づいて対象者の磁気共鳴イメージングスキャンを行う。スキャンは、k空間の直交系サンプリング、六角形サンプリング、又は圧縮センシングサンプリングを含み得る。
様々なユニット又はモジュール22、24、26、及び28が、ワークステーション36の電子プロセッサや電子処理装置30等の電子データ処理装置によって、又はネットワークによってワークステーション36に動作可能に接続されるネットワークベースのサーバコンピュータ等によって適切に具体化される。更に、開示されるユーザインターフェイス、信号及び雑音マップ構築、並びにパラメータ選択技法が、開示される技法を実行するために電子データ処理装置によって読取可能であり、電子データ処理装置によって実行可能な命令(例えばソフトウェア)を記憶する非一時的記憶媒体を使用して適切に実施される。
図2を参照すると、スマートセレクト処理の一実施形態が流れ図で示される。簡潔に要約すると、各方向のSENSE低減係数が自動で計算される。プリスキャンが、各コイル素子又はコイル素子群の雑音度を求める。様々なSENSE低減係数を用いて様々な方向についてボクセルごとのSNRが計算される。オペレータは全体的な低減係数を指定することができ、システムは、選択されたSENSE係数のSNRペナルティをオペレータに知らせる。或いは、オペレータが許容可能な最大SNRを指定し、システムが低減係数を選択する。この方法をより詳細に見ると、ステップ又はユニット44が、複数Nのデータチャネルからのデータ及びQボディデータ(Q-body data)を受け取るプリスキャン又はコイル調査スキャンを行う。各チャネルは、1個(又は複数個の)コイル素子に対応する。受け取られたデータから雑音共分散行列Ψが構築される。
ステップ又はユニット46が、受け取られたデータに対して多断面再構成を行い、多断面再構成では、臨床スキャンのジオメトリを有するオブリークスライスへとデータを変換する。ステップ又はユニット48が、コイル感度及び正則化データを計算し、患者又は関心領域によって占有されるイメージングボリュームの一部の感度を算出する。
ステップ又はユニット50により、コイル感度データ及び正則化データから信号マップ及び雑音マップが構築される。信号マップは低減係数の各組について同じであるのに対し、雑音マップは低減係数に基づいて異なる。最適なSNRを有する方向を明らかにするために、雑音マップの構築が様々な低減係数を使って繰り返される。一実施形態では、行列H=(SΨ−1S+R−1−1Ψ−1を用いて雑音マップが計算され、Sはコイル感度マップであり、Rは患者又は関心領域によって占有されるイメージボリューム内で予期される最大MR信号のマップであり、Ψは共分散行列である。雑音(r)は、対角(diag)(HΨH)として計算される。
ステップ又はユニット52により、信号マップが雑音マップで割られ、SNRマップがもたらされる。SNRマップは、ボクセルごとのSNR値を含む。ステップ又はユニット54は、

等、SNRマップの平均又は他の基準を計算し、Nはボクセルの数であり、snr(r)はr番目のボクセルの信号対雑音比(SNR)である。図1に関して論じたように、各雑音マップの最も高いSNRを求めるために、雑音マップに対応する全ボクセルの平均SNRが使用される。この情報に基づき、方向ごとの低減係数を含む低減係数が求められ、オペレータは画質(SNR)とイメージング速度(低減係数)との間の最適なバランスを選択することができる。
図3は、1つ又は複数のプロセッサ又は個別ユニット若しくはモジュールによって実施され得る、低減方向に沿った加速係数等のパラレル加速パラメータを自動で最適化する一方法を流れ図で示す。ステップ56で、全体的な低減係数又はSNR最低値がユーザインターフェイスによって受け取られる。ステップ58で、高速フィールドエコー(FFE)スキャン等のコイル調査スキャンが行われる。コイル調査は、k空間の直交系サンプリング、六角形サンプリング、又は圧縮センシングサンプリングによって行われ得る。或いは事前スキャンが行われる。
ステップ59で、受け取られた全体的な低減係数を満たす複数の低減係数が生成される。全体的な低減係数は最小SNRに基づき得る。例えば、4の全体的な低減係数では、(4、1)、(3.5、8/7)、(3、4/3)等が4の全体的な低減係数を満たす。低減係数はスキャン時間を加速し、又はk空間のサンプリングを減らす。低減係数は、2つのエンコード方向の低減係数を含む。一実施形態では、生成される低減係数の組が粗い格子又は離散値を含み、これらは粗値の最初の最も高いSNRに基づいてより小さい増分で更に洗練させることができる。
ステップ60で、信号マップ及び複数の雑音マップが構築される。このステップは、図2に関して説明された処理を含む。生成された低減係数の組について、及び/又は最小SNRに基づき、ステップ60で信号マップ及び雑音マップが構築される。図2に関して説明されたスマートセレクト処理を続け、ステップ62で、1組の低減係数に基づいて各雑音マップのボクセルごとにSNRが算出される。SNRは、信号マップを雑音マップで割ることによって計算される。関心対象又は全イメージングボリュームにわたる平均SNRも計算される。
判定ステップ64で、最初の最も高いSNR及び対応する低減係数付近のより小さい低減係数の増分により、雑音マップの構築及びSNRの算出を繰り返す反復が実行される。
ステップ66で、対応する雑音マップの最も高いSNRに基づいて低減係数が選択される。全体的な低減係数未満の低減係数の探索は低減方向に沿う。ステップ68で、選択された低減係数を用いて磁気共鳴イメージングスキャンが行われる。スキャンは、例えば選択された低減係数によって圧縮されるk空間のSENSEサンプリングを用いて行われ得る。
本発明の特定の例示的実施形態によれば、この例示的MRシステムは、所与のイメージングジオメトリ、対象の範囲、受信コイル構成に関するSENSE係数及びサンプリング手法スキーム、並びに/又は比較的低分解能のG因子評価を実質的に行う誤り伝搬の高速評価によって初めに入手可能な追加情報に応じてSENSEの性能を最適化することができる。使用されるSENSEスキャンパラメータは、パラメータ空間を形成する。これらのパラメータは、例えば検査パラメータによって利用者に提示され得る。使用可能な情報に基づき、この例示的MRシステムは、最適な(又は最良の)パラメータセットを突き止めることができ、又は特定のSENSEパラメータを或る有用なものに変えることにより、何らかの既定の検査パラメータに固有の制約(例えばSNR損失、G因子の分布/統計、特定の関心領域(ROI)等)を満たすことができる。利用者に自身の選択の結果を知らせ、問題を解決するための選択肢を与えることも可能である。
最終的な画質が、最適な感度エンコードプロセス、及びアンフォールディング問題の良いコンディショニングによって得られる(G因子(1))。低減FOV内の単一位置に関するSENSE問題は、例えばSx=yと見なされ得る。
この例では、互いに重なる(均一なk空間のアンダーサンプリングによるフォールディング)様々な位置からのNの未知のボクセル信号寄与度がベクトルxに集約され得る。受け取られる信号が感度エンコードされ、M個の受信コイルを用いて測定される(M≧N)。それらの対応する信号がベクトルyにまとめられる。行列Sは、(N×M)の感度行列を示す。選択された低減係数R(2D SENSEの場合はRtotal=R1+R2)及び選択されたk空間固有のサブサンプリング方式は、フルFOVのうちのどのボクセルが、測定過程中に適切に感度エンコードされて互いに重なるのかを定める。つまり、Sx=yの実際の形は、ベクトルΔに集約されるべき選択された実験パラメータに依存する。
この行列をSx=yに反転するとベクトルxがもたらされ、ベクトルxは対応するNのボクセル固有信号を含む。例えば次式が成立する。
行列

はSの疑似反転であり、そのノルムが測定値から最終画像内への誤り伝搬を示す。1つのあり得るノルムであるg因子は、良い数値的条件のための正規化された測度を与えることができる。例えば次式が成立する。
g因子は行列Sに依存し、その行列Sは選択されたサンプリングパラメータΔに依存する。従って、選択されたSENSE検査パラメータΔ(R1、R2、ダイヤモンドサンプリングパターン等)が最適か、又はより良い結果を得るために異なるように選択されるべきかを判定するためにg因子が使用され得る。これを行うために、あり得る全ての設定Δの間でgを最小化する必要がある。しかし、

は低減FOV内の1つのボクセルしか示さず、低減FOV(2Dイメージングでは行列L×K、3Dでは更に高次元)内のほぼ全てのボクセルの妥協点が見つけられなければならない。これは数値的に困難であり、評価がバックグラウンドに隠され得るが、加速のための幾つかの手段が考えられる。
rFOV内の全てのL×KボクセルがΔパラメータを最適化するために使用されない場合、問題の大きさが低減され得る。そのために、低減FOVが均一又は無作為に「アンダーサンプリング」され(低分解能手法)、それにより検討すべきg因子要素の数を減らすことができる。このことは、コイル感度がどちらかと言えば平滑な空間内の関数であることも含意する。更に、臨床的に特に関心のある画像内の重要領域(ROI)に評価を更に集中させることができる。
問題の大きさを低減する第2の手法は、関与する全てのコイルを(SVD又は同様の手法を用いて)仮想コイルシステムに変換することである。従って、Nの有意のコイルだけがこの評価内で検討される必要がある。
しかし、疑似反転の基礎的計算の要素として全てのコイルを使用して位置ごとにS^H Sを直接計算するのも効率的であり、コイル情報の一部を軽視するのを防ぐ。
或いは、

の代わりに更に単純なノルムが評価に使用され得る。
統計的過程を適用して低減FOV内のg因子の分布を解析し、様々なパラメータセットΔの比較を可能にすることができる。更に、最適なΔの探索は、最適化の進行中に適合する/高めることができるΔ空間内及びrFOV空間内の粗分解能から始めて階層的方法で行われても良い。
SENSEアンフォールディングに概して適用され得る追加情報(雑音相関行列、正則化画像、又は最も単純な事例では、信号がそこから予期される領域にコイル感度を制限するバイナリマスク)も最適化におけるg因子の計算に含めることができる。
更に、エンコード問題の構造に関する情報が組み込まれ得る。このことは、特に例えばトレーニングデータセットがアンフォールディングに使用され得るk−t SENSE応用において特に興味深い場合がある。この情報が測定前に既に(例えばライブラリからの)概算として、又は測定モデルから入手可能である場合、この患者又はイメージング問題に専用の適切なk−tサンプリング方式及び低減係数の選択への追加入力としてこの情報が使用され得る。
例えば、固定加速係数(R total=4)による3D FFEスキャンが所与のアレイコイルを使用して行われる。これは、例えばスキャンが20秒未満であることを確実にするための臨床的要件の理由からであり得る。利用者は、スキャンジオメトリ、読出し方向、最終的なイメージング行列、及びR totalを定める。これらの制約及びSENSEリファレンススキャン情報を所与とし、スキャナは、あり得るどちらのアンダーサンプリング因子R1及びR2を使用すべきか、並びにどのサンプリング方式を使用すべきかを決定する。システムの助けなしにこの決定を行うのは複雑であり得る。対照的に、本発明の例示的実施形態によれば、この決定がシステムによって行われ、画質の改善を助けながら使い易さ及び単純さを高めることができる。
本明細書で示された特定の例示的実施形態に関連し、特定の構造上の及び/又は機能的な特徴が定められた要素及び/又は構成要素内に組み込まれているものとして記載されていることが理解される。しかし、これらの特徴は、適切な場合は同じ又は同様の利益を得るために他の要素及び/又は構成要素に同様に組み込まれても良いと考えられる。所望の応用に適した他の代替的実施形態を実現するために、例示的実施形態の様々な態様が必要に応じて選択的に採用され、それにより他の代替的実施形態が自らに組み込まれる態様のそれぞれの利点を実現し得ることも理解すべきである。
本明細書に記載した特定の要素又は構成要素がその機能をハードウェア、ソフトウェア、ファームウェア、又はそれらの組合せによって適切に実装し得ることも理解すべきである。加えて、一緒に組み込まれるものとして本明細書に記載された特定の要素は、適切な条件下では独立型要素とすることができ、又は他の方法で分けられても良い。同様に、或る特定の要素によって実行されるものとして記載された複数の特定の機能が、個々の機能を実行するために独立に動作する複数の別個の要素によって実行されても良く、又は特定の個別機能が分割され、提携して動作する複数の別個の要素によって実行されても良い。或いは、そうではなく互いに異なるものとして本明細書で説明され且つ/又は図示された一部の要素又は構成要素が、必要に応じて物理的に又は機能的に組み合わせられても良い。
要約すれば、本明細書が好ましい実施形態に関して記載されてきた。当然ながら、本明細書を読んで理解すれば修正形態及び改変形態に他者は気付かれよう。かかる修正形態及び改変形態が添付の特許請求の範囲及びその等価物の範囲に含まれる限り、本発明がかかる全ての修正形態及び改変形態を含むものとして解釈されるよう意図される。つまり、上記で開示された及び他の特徴及び機能又はその代替策が、他の多くの異なるシステム又は応用例に望ましく集約されても良いこと、及びここでは予測外の又は予期せぬ様々な代替形態、修正形態、改変形態、又はその中の改善点が当業者によって後でもたらされる場合があり、それらも同様に添付の特許請求の範囲によって包含されるよう意図される。

Claims (20)

  1. パラレル磁気共鳴イメージングシステムであって、
    複数のコイル素子を有する少なくとも1個の無線周波(RF)コイルと、
    前記少なくとも1個のRFコイルによる対象者のプリスキャン又は事前スキャンから、1つ又は複数の低減方向に沿った様々な組の低減係数に基づいて信号マップ及び複数の雑音マップを構築するスマートセレクトユニットと、
    信号対雑音比(SNR)及びスキャン速度を最適化する、雑音マップに対応する1組の低減係数を選択するパラレルイメージングパラメータユニットと、
    前記選択された低減係数に基づいて前記対象者の磁気共鳴イメージングスキャンを行うシーケンスコントロールと
    を含む、システム。
  2. 前記1組の低減係数が、2つの低減方向について2つの低減係数を含む、請求項1に記載のシステム。
  3. 最大イメージング時間又は全体的な低減係数の1つを受け取るユーザインターフェイス
    を更に含み、
    前記スマートセレクトユニットが、前記最大イメージング時間又は前記全体的な低減係数に基づいて低減係数の複数の組を生成する、
    請求項1又は2の何れか一項に記載のシステム。
  4. 前記1組の低減係数が複数の低減係数を含み、各低減係数がエンコード方向に対応する、請求項1乃至3の何れか一項に記載のシステム。
  5. 最小SNRを受け取るユーザインターフェイス
    を更に含み、
    前記スマートセレクトユニットが、受け取られた前記最小SNRを満たす全体的な低減係数を生成する、
    請求項1又は2に記載のシステム。
  6. 前記パラレルイメージングパラメータユニットが、少なくとも前記受け取られた最小SNRで画像を生成するためのイメージング時間を最小化する前記全体的な低減係数を選択する、請求項5に記載のシステム。
  7. 前記低減係数を選択することが、各エンコード方向に沿って探索することを含む、請求項1乃至6の何れか一項に記載のシステム。
  8. 前記スマートセレクトユニットが、
    前記プリスキャン又は前記事前スキャンから1組の3次元(3D)データを構築し、
    前記イメージングスキャンのジオメトリと共にチャネルデータ及びQボディデータを求めるために、前記3Dデータの多断面再構成を行い、
    前記求められたチャネルデータ及びQボディデータに基づいてコイル感度及び正則化データを計算し、
    計算されたコイル感度及び正則化データに基づいて信号マップ及び雑音マップを生成する、
    請求項1乃至7の何れか一項に記載のシステム。
  9. 前記最も高いSNRが全ボクセルの平均SNRを含む、請求項1に記載のシステム。
  10. パラレル磁気共鳴イメージングの方法であって、前記方法は、
    対象者の事前スキャンから、受け取られた全体的な低減係数に基づいて信号マップ及び複数の雑音マップを構築するステップであって、各雑音マップは1つ又は複数の低減方向に沿った様々な組の低減係数に基づいて求められる、構築するステップと、
    信号対雑音比(SNR)基準及びスキャン速度基準を満たす雑音マップに対応する低減係数を選択するステップと、
    前記選択された低減係数に基づき、前記対象者の磁気共鳴イメージングスキャンを行うステップと
    を含む、方法。
  11. 前記低減係数を選択するステップが、2つの低減係数を選択するステップを含み、各低減係数は低減方向に関する、請求項10に記載の方法。
  12. 最大イメージング時間又は全体的な低減係数の1つを受け取るステップと、
    前記最大イメージング時間又は前記全体的な低減係数に基づいて複数組の低減方向を生成するステップと
    を更に含む、請求項10又は11に記載の方法。
  13. 低減係数の各組が複数の低減係数を含み、各低減係数がエンコード方向に対応する、請求項10乃至12の何れか一項に記載の方法。
  14. 最小SNR基準を受け取るステップと、
    前記受け取られたSNR基準を満たす全体的な低減係数を生成するステップと
    を更に含む、請求項10又は11に記載の方法。
  15. 前記受け取られた最小SNR基準を少なくとも満たす画像を生成するためのイメージング時間を最小化する前記低減係数及び/又は低減方向を選択するステップ
    を更に含む、請求項14に記載の方法。
  16. 前記低減係数を各エンコード方向に沿って離散的刻みで探索するステップと、
    最も高いSNR前後の前記低減係数の探索をより小さい刻みで洗練させるステップと
    を更に含む、請求項10乃至15の何れか一項に記載の方法。
  17. 前記事前スキャンから1組の3次元データを構築するステップと、
    前記イメージングスキャンのジオメトリと共にチャネルデータ及びQボディデータを求めるために、前記3Dデータの多断面再構成を行うステップと、
    前記求められたチャネルデータ及びQボディデータに基づいてコイル感度及び正則化データを計算するステップと、
    求められたコイル感度及び正則化データに基づいて信号マップ及び雑音マップを生成するステップと
    を更に含む、請求項10乃至16の何れか一項に記載の方法。
  18. 請求項10乃至17の何れか一項に記載の方法を実行するために1つ又は複数の電子データ処理装置を制御するソフトウェアを担持する、非一時的コンピュータ可読記憶媒体。
  19. 請求項10乃至17の何れか一項に記載の方法を実行する、電子データ処理装置。
  20. パラレル磁気共鳴イメージングシステムであって、
    複数のコイル素子を有する少なくとも1個のRFコイルと、
    画像データを生成するための最低限の画質及び/又は最大時間の少なくとも1つを受け取るユーザインターフェイスと、
    1個又は複数個のプロセッサであって、
    前記少なくとも1個のRFコイルによる対象者の事前スキャンから、1つ又は複数の低減方向に沿って複数の低減係数のそれぞれについて信号マップ及び複数の雑音マップを構築し、
    少なくとも前記受け取られた最低限の画質で、及び/又は前記受け取られた最大時間内に画像データを生成する低減係数を選択し、
    前記選択された低減係数に基づいて前記対象者の磁気共鳴イメージングスキャンを前記イメージングシステムに行わせる、
    1個又は複数個のプロセッサと
    を含む、システム。
JP2015562491A 2013-03-13 2014-03-12 パラレルイメージング加速パラメータの自動最適化 Active JP6430417B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US201361778582P 2013-03-13 2013-03-13
US61/778,582 2013-03-13
PCT/IB2014/059667 WO2014141086A1 (en) 2013-03-13 2014-03-12 Automatic optimization of parallel imaging acceleration parameters

Publications (3)

Publication Number Publication Date
JP2016509940A true JP2016509940A (ja) 2016-04-04
JP2016509940A5 JP2016509940A5 (ja) 2017-04-13
JP6430417B2 JP6430417B2 (ja) 2018-11-28

Family

ID=50440712

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2015562491A Active JP6430417B2 (ja) 2013-03-13 2014-03-12 パラレルイメージング加速パラメータの自動最適化

Country Status (5)

Country Link
US (1) US10107884B2 (ja)
EP (1) EP2972450B1 (ja)
JP (1) JP6430417B2 (ja)
CN (1) CN105308473B (ja)
WO (1) WO2014141086A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2019530491A (ja) * 2016-08-25 2019-10-24 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 補正された感度符号化磁気共鳴イメージング

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9689945B2 (en) * 2010-08-20 2017-06-27 Koninklijke Philips Electronics N.V. Virtual coil emulation in parallel transmission MRI
WO2014165050A1 (en) * 2013-03-12 2014-10-09 The General Hospital Corporation Method for increasing signal-to-noise ratio in magnetic resonance imaging using per-voxel noise covariance regularization
EP3470868A1 (en) 2017-10-16 2019-04-17 Koninklijke Philips N.V. Quantitative measurement of relaxation times in magnetic resonance imaging
US11023785B2 (en) * 2018-07-23 2021-06-01 International Business Machines Corporation Sparse MRI data collection and classification using machine learning
CN111830449B (zh) * 2019-04-19 2022-10-18 西门子(深圳)磁共振有限公司 呼吸信号提取方法及装置、磁共振成像系统、存储介质
EP3910358A1 (en) * 2020-05-15 2021-11-17 Koninklijke Philips N.V. Automated adjustment of the undersampling factor
WO2022032133A1 (en) * 2020-08-06 2022-02-10 University Of Pittsburgh - Of The Commonwealth System Of Higher Education Optimizing magnetic resonance imaging protocols for increased speed
CN112630709A (zh) * 2020-11-03 2021-04-09 成都易检医疗科技有限公司 金属植入物磁共振成像方法、系统、终端及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002505613A (ja) * 1998-04-17 2002-02-19 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 磁気共鳴による撮像方法及び装置
JP2005525185A (ja) * 2002-05-13 2005-08-25 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 磁場強度の関数としてのパラレルイメージングにおけるリダクションファクタの固有の限界
WO2005091011A1 (en) * 2004-03-12 2005-09-29 Koninklijke Philips Electronics N.V. Prescan for optimization of mri scan parameters
WO2012028955A2 (en) * 2010-09-01 2012-03-08 Commissariat A L Energie Atomique Et Aux Énergies Alternatives Method for performing parallel magnetic resonance imaging
WO2012085810A2 (en) * 2010-12-22 2012-06-28 Koninklijke Philips Electronics N.V. Rapid parallel reconstruction for arbitrary k-space trajectories

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
US6714010B2 (en) * 2001-04-20 2004-03-30 Brigham And Women's Hospital, Inc. Combining unfold with parallel magnetic resonance imaging
JP4030805B2 (ja) * 2002-06-07 2008-01-09 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー コイル感度マップ作成方法およびmri装置
JP4607431B2 (ja) * 2003-05-08 2011-01-05 株式会社東芝 Mri装置
DE10339019B4 (de) * 2003-08-25 2007-02-08 Siemens Ag Verfahren zur Bestimmung eines Beschleunigungsfaktors einer parallelen Bildaufnahme
WO2006003553A1 (en) * 2004-06-28 2006-01-12 Koninklijke Philips Electronics N.V. Parallel magnetic resonance imaging
US7449888B1 (en) * 2005-07-27 2008-11-11 General Electric Company Method and apparatus for multi-dimensional parallel MR imaging
US7429862B2 (en) 2006-04-18 2008-09-30 Siemens Aktiengesellschaft PPA MRT method and apparatus with scaled contribution of k-space data to the reconstructed image
JP4464373B2 (ja) 2006-07-12 2010-05-19 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Mri装置
US7382127B2 (en) * 2006-09-15 2008-06-03 General Electric Company System and method of accelerated MR propeller imaging
US7394252B1 (en) 2007-05-03 2008-07-01 The General Hospital Corporation Regularized GRAPPA reconstruction
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
US8638096B2 (en) * 2010-10-19 2014-01-28 The Board Of Trustees Of The Leland Stanford Junior University Method of autocalibrating parallel imaging interpolation from arbitrary K-space sampling with noise correlations weighted to reduce noise of reconstructed images

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002505613A (ja) * 1998-04-17 2002-02-19 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 磁気共鳴による撮像方法及び装置
JP2005525185A (ja) * 2002-05-13 2005-08-25 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 磁場強度の関数としてのパラレルイメージングにおけるリダクションファクタの固有の限界
WO2005091011A1 (en) * 2004-03-12 2005-09-29 Koninklijke Philips Electronics N.V. Prescan for optimization of mri scan parameters
JP2007528767A (ja) * 2004-03-12 2007-10-18 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Mri走査パラメータの最適化のための前段走査
WO2012028955A2 (en) * 2010-09-01 2012-03-08 Commissariat A L Energie Atomique Et Aux Énergies Alternatives Method for performing parallel magnetic resonance imaging
WO2012085810A2 (en) * 2010-12-22 2012-06-28 Koninklijke Philips Electronics N.V. Rapid parallel reconstruction for arbitrary k-space trajectories

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
PAUL T. WEAVERS ET AL.: "Acceleration Apportionment: A Method of Improved 2D SENSE Acceleration Applied to 3D Contrast-Enhanc", MAGNETIC RESONANCE IN MEDICINE, vol. 71, no. 2, JPN6018004585, 28 February 2013 (2013-02-28), pages 672 - 680, XP055122491, DOI: doi:10.1002/mrm.24700 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2019530491A (ja) * 2016-08-25 2019-10-24 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. 補正された感度符号化磁気共鳴イメージング
JP7025410B2 (ja) 2016-08-25 2022-02-24 コーニンクレッカ フィリップス エヌ ヴェ 補正された感度符号化磁気共鳴イメージング

Also Published As

Publication number Publication date
WO2014141086A1 (en) 2014-09-18
EP2972450A1 (en) 2016-01-20
US10107884B2 (en) 2018-10-23
CN105308473A (zh) 2016-02-03
CN105308473B (zh) 2018-11-27
EP2972450B1 (en) 2021-02-17
JP6430417B2 (ja) 2018-11-28
US20160018498A1 (en) 2016-01-21

Similar Documents

Publication Publication Date Title
JP6430417B2 (ja) パラレルイメージング加速パラメータの自動最適化
CN107106039B (zh) 切片中的沿着一维延伸的磁共振指纹识别
CN102159965B (zh) 用于mri的b1映射和b1l匀场
CN104781685B (zh) 用于使用非相干采样和冗余haar子波的动态磁共振成像的图像重建
JP5525596B2 (ja) 加速b1マッピング
US9317917B2 (en) Method, reconstruction device, and magnetic resonance apparatus for reconstructing magnetic resonance raw data
US10353023B2 (en) Calculating MRI RF coil sensitivities using interpolation into an enlarged field of view
CN107072586B (zh) 磁共振成像装置
US11762049B2 (en) Medical data processing apparatus, medical data processing method, and magnetic resonance imaging apparatus
JP7272818B2 (ja) データ処理装置、方法及びプログラム
JP2013505046A (ja) Mriにおけるrfパワー及びrfフィールド均一性の同時最適化
US20170315193A1 (en) Magnetic resonance fingerprinting using a spin-echo pulse sequence with an additional 180 degree pulse
JP2016516502A (ja) k空間サンプリング方法
WO2015140277A1 (en) Control of magnetic resonance imaging acquisition using modeling
CN110542872B (zh) 磁共振成像方法及设备、相位信息的获取方法及装置
JP2016526459A (ja) コイル感度を用いた補正磁気共鳴画像法
US11061097B2 (en) Data driven correction of phase depending artefacts in a magnetic resonance imaging system
JP5211403B2 (ja) 磁気共鳴イメージング装置
JP2023525165A (ja) アンダーサンプリング係数の自動化された調整
JP2018196661A (ja) 磁気共鳴撮像装置および体動補正方法
JP7183048B2 (ja) 磁気共鳴イメージングシステム、磁気共鳴イメージング方法、及び磁気共鳴イメージングプログラム
US9069050B2 (en) Varying blipped trajectory in MRI
Aelterman et al. Automatic high-bandwidth calibration and reconstruction of arbitrarily sampled parallel MRI
US10067212B2 (en) Dynamic method and apparatus for radial acquisition of magnetic resonance data
JP2023180045A (ja) 医用データ処理装置、医用データ処理方法、医用データ処理プログラムおよび磁気共鳴イメージング装置

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170309

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20170309

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20180125

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180220

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180501

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20181031

R150 Certificate of patent or registration of utility model

Ref document number: 6430417

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

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