JP6088711B2 - 磁気共鳴イメージングのための息止め検出 - Google Patents

磁気共鳴イメージングのための息止め検出 Download PDF

Info

Publication number
JP6088711B2
JP6088711B2 JP2016541371A JP2016541371A JP6088711B2 JP 6088711 B2 JP6088711 B2 JP 6088711B2 JP 2016541371 A JP2016541371 A JP 2016541371A JP 2016541371 A JP2016541371 A JP 2016541371A JP 6088711 B2 JP6088711 B2 JP 6088711B2
Authority
JP
Japan
Prior art keywords
space
interest
area
subject
samples
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
Application number
JP2016541371A
Other languages
English (en)
Other versions
JP2017500128A (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 JP2017500128A publication Critical patent/JP2017500128A/ja
Application granted granted Critical
Publication of JP6088711B2 publication Critical patent/JP6088711B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/567Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution gated by physiological signals, i.e. synchronization of acquired MR data with periodical motion of an object of interest, e.g. monitoring or triggering system for cardiac or respiratory gating
    • G01R33/5676Gating or triggering based on an MR signal, e.g. involving one or more navigator echoes for motion monitoring and correction
    • 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/283Intercom or optical viewing arrangements, structurally associated with NMR apparatus
    • 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/4818MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • 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/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56509Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling
    • 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/4818MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • G01R33/4824MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory
    • G01R33/4826MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory in three dimensions

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Signal Processing (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biophysics (AREA)
  • Cardiology (AREA)
  • Physiology (AREA)
  • Power Engineering (AREA)
  • Pulmonology (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Description

本発明は磁気共鳴イメージングの分野に関する。
磁気共鳴イメージング(MRI)法は、関心被験者またはその少なくとも一エリアの2次元または3次元画像を形成する目的で、磁場と核スピンとの間の相互作用を用いる。これらの方法は、今日広く使われている。特に医療診断の分野においてはそうである。その理由は、イオン化放射線を必要とせず、通常は非侵襲的だからである。MRIは、例えば、身体の構造的異常、例えば腫瘍成長を可視化するイメージング手法として用いられる。
MRI装置は強力な磁場を用いて、身体内の一部の原子核の磁化を整列させ、無線周波数(RF)場がシステマチックにこの磁化のアライメントを修正する。これにより、原子核が、スキャナにより検出可能な回転磁場を生成する。スキャナが、サンプリングされた身体のエリアの画像を構成する。磁場グラジエントにより、位置が異なる原子核は、異なる速さで回転する。異なる方向のグラジエントを使うことにより、2次元画像または3次元ボリュームを任意の芳香で取得できる。
MRIは、一般に、検査される関心被験者が、強力かつ一様な磁場B0中に配置されることを必要とする。磁場B0の方向は同時に、測定が基づく座標系の軸、通常はz軸を確定する。印加されるRF場による磁化の時間的バリエーションは、受信RFアンテナにより検出され得る。受信RFアンテナは、MR装置の検査ボリューム内に構成及び向けられ、磁化の時間的バリエーションがz軸に垂直な方向で測定されるようになっている。
身体内の空間的解像度は、磁場グラジエントを切り替えることにより実現し得る。磁場グラジエントは3つの主軸に沿って延在し、一様な磁場に重畳され、スピン共鳴周波数の線形な空間的依存性を生じる。受信アンテナでピックアップされた信号は、身体/被験者の異なる位置にリンクし得る異なる周波数成分を含む。
受信アンテナにより取得される信号データは、空間的周波数ドメインにマッチし、K空間データと呼ばれる。K空間データは、一般的に異なる位相エンコードステップで取得された複数のラインを含む。各ラインは多数のサンプルを収集することによりデジタル化される。K空間データのサンプルのセットは、例えば逆フーリエ変換により、MR画像に変換される。
現在、MRIスキャンを行う2つの方法がある、一方は息止めをするものであり、他方は自由呼吸できるものである。第1のオプションでは、関心被験者はMRIシステムのオペレータにより、スキャンの時間中、呼吸を止めるように指示されねばならない。自由呼吸スキャンは関心被験者とのインターラクションを必要としないが、息止めスキャンより長い時間がかかる。また、画質の一貫性も自由呼吸スキャンにおける潜在的な問題である。よって、スキャン効率を向上しつつ、画質と一貫性を維持することが望ましい。
このコンテキストにおいて、特許文献1は磁気共鳴イメージングシステムに関する。このシステムでは、静磁場が患者に印加され、グラジエント磁場と、所定パルスシーケンスによる励起パルス信号が患者に印加され、患者の選択スライスにおいて磁気共鳴減少を生じる。磁気共鳴現象の磁気共鳴データが取得され、磁気共鳴データから磁気共鳴画像が得られる。本システムは、患者に断続的に動かないように促すアナウンスメントセクションを有する。このシステムでは、データ取得セクションは、アナウンスメントセクションのアナウンスに応じて患者がじっとしている間だけ、制御セクションの管理下で動作し、それにより所定ボリュームの単位で磁気共鳴データを間欠的に取得する。
さらに、特許文献2は磁気共鳴イメージング装置について記載している。パルスシーケンスは、被験者の身体の動きの位置を検出する身体動き検出シーケンスと、被験者の画像を取得する写真シーケンスとを有する。制御手段は、身体動き検出シーケンスを繰り返し実行し、身体動き検出シーケンス及び写真シーケンスにより検出された身体動き情報に基づき、呼吸停止指令が送信手段により被験者に送られ、呼吸再開指令も送信手段により被験者に送られ、呼吸停止指令後でも、身体動き検出シーケンスを継続することにより、身体動き位置が所定範囲内である時、写真シーケンスが実行される。
またさらに、特許文献3は呼吸モニタを有するNMRシステムを説明している。モニターは、患者の横隔膜が同じ基準点に位置している一連の息止めを、患者が行えるようにする、視覚的フィードバックを患者に提供する。これにより、ぼけや画像アーティファクトを生じずに、一連の息止めにわたりNMRデータを取得できるようになる。息止めの間、ナビゲータパルスシーケンスを用いて、横隔膜位置を測定できるNMRデータを収集し、各息止め中、パルスシーケンスを変えて、NMR画像データを収集する。
特許文献4は、NMR画像再構成の方法を説明している。画像再構成に必要なデータは、中心K空間ビューと周辺K空間ビューに分割される。スキャン中に患者の呼吸を示すNMRナビゲータ信号を取得して、呼吸が狭い取得ウィンドウ内にある時、第1のゲーティング信号を生成し、呼吸がより広い取得ウィンドウ内にある時、第2のゲーティング信号を生成する。第1のゲーティング信号が生成された時、中心K空間ビューが取得され、第2のゲーティング信号が生成された時、周辺K空間ビューが取得される。
非特許文献1は腹部イメージングにおける動きアーティファクトを低減するための確率されたストラテジとして、息止めを説明している。呼吸の始まりのナビゲータベースの検出と組み合わせて、圧縮センシングと並列イメージングを用いて、任意の時点まで取得されたアンダーサンプリングされたデータからの画像再構成をサポートするサンプリングパターンがデザインされる。それにより、スキャンの終了と、一貫性のあるデータのみからの再構成が可能になり、動きアーティファクトが抑制される。空間的解像度は、サンプリング密度の下限により制約され、スキャンごとに向上し、任意の息止め期間中の信号対ノイズ比及びアンダーサンプリングアーティファクトとの妥協点になる。
米国特許第4,878,499A明細書 特開第2007−029250A号 米国特許第5,363,844A明細書 米国特許第6,144,874A明細書
Nadine Gdaniecetal.著「robust abdominal imaging within completebreath−holds」(Magnetic resonance in medicine、vol. 71、no. 5、pages 1733−1742、ISSN0740−3194)
本発明の一目的は、磁気共鳴イメージング方法、磁気共鳴イメージングシステムをアップグレードするソフトウェアパッケージ、及びMRイメージングシステムであって、関心被験者の動きに基づいて効率的に画質が高くアーティファクトが少ない磁気共鳴イメージングを可能とし、特に呼吸により動くエリアにおける関心被験者のイメージングの場合に画質を高くかつアーティファクトを少なくできるものをを提供することである。
この目的は、関心被験者のエリアの磁気共鳴(MR)イメージングをする方法により実現される。該方法は、
前記関心被験者に息止め指令を発するステップと、
前記関心被験者の動き検出を行って、前記関心被験者のエリアにおける息止め状態を検出するステップと、
前記関心被験者のエリアのK空間サンプリングを行って、前記関心被験者のエリアにおける息止め状態の検出の前に、前記関心被験者のエリアをカバーしている前記K空間サンプルのセットからK空間サンプルの第1のサブセットを取得するステップと、
前記関心被験者のエリアにおける息止め状態を検出したとき、一解像度で前記関心被験者のエリアのK空間サンプリングを行うステップであって、K空間サンプリングを行って前記関心被験者のエリアをカバーするK空間サンプルの第2のサブセットを取得することを含むステップと、
前記関心被験者のエリアをカバーするK空間サンプルを処理して前記関心被験者のエリアのMR画像を取得するステップであって、前記K空間サンプルの第1と第2のサブセットを処理して前記関心被験者のエリアの前記MR画像を取得することを含むステップとを含み、
前記K空間サンプルの第1と第2のサブセットは共に前記一解像度で前記関心被験者のエリアカバーし、
前記関心被験者の動き検出を行って前記関心被験者のエリアにおける息止め状態を検出するステップは、K空間サンプルの前記第1の及び/又は第2のサブセットの取得中に少なくとも部分的に行われ、
前記関心被験者のエリアをカバーするK空間サンプルを処理して前記関心被験者のエリアのMR画像を取得するステップはさらに、検出された前記関心被験者のエリアの動きに応じてK空間サンプルの前記第1の及び/又は前記第2のサブセットのうちのK空間サンプルを排除及び/又は加重するステップを含む。
この目的は、コンピュータに上記の方法を実行させる一組の命令を含む、ストレージ装置、フロッピーディスク、コンパクトディスク、CD、デジタルバーサタイルディスク、DVD、ブルーレイディスク、またはランダムアクセスメモリ、RAMなどのコンピュータ読み取り可能媒体が提供される。
この目的はさらに、コンピュータ利用可能プログラムコードであって上記の方法を実行するよう適応されたものを含むコンピュータ利用可能媒体を含むコンピュータプログラム製品により達成できる。
さらに、この目的は、上記方法によるMRイメージングシステムを制御する命令を含む、MRイメージングシステムをアップグレードするソフトウェアパッケージにより達成できる。
さらにまた、この目的は、
MRイメージングシステムの検査空間に配置された関心被験者のエリアの画像表現を提供するMRイメージングシステムであって、
上記の磁気共鳴イメージング方法を実行するように適応されている、MRイメージングシステムにより達成される。
動きは、画質劣化と一貫性問題を生じる関心被験者の身体のMRIにおける第1の問題である。これは、特に胸部の呼吸の動きを指す。息止めをしてMRIスキャンを行う際に、画質を改善する技術水準のアプローチは、取得されたK空間データに基づいたデータ処理に基づき、アーティファクト、特にMR画像が生成される関心被験者のエリア内のものは、関心被験者の動きに基づく。
関心被験者のエリアにおける関心被験者の動き検出では、K空間サンプリングにより取得されるデータの信頼性は、K空間サンプルの生成の段階ですでに改善され得る。動きの制御、特に呼吸運動の制御は、改善できる。息止め指令と動き停止を発する間の相対的関係が、関心被験者のエリアにおける動き検出により確立でされ得るからである。関心被験者の呼吸状態のフィードバック制御は、息止め指令の発行に加えて実現し得る。したがって、関心被験者のエリアのK空間サンプリングは、モニターされた息止め状態で行え、取得されたK空間サンプルにより、改善された画質で、例えば低減されたアーティファクトでMR画像を生成できる。データ処理、例えばアーティファクトを減らすデータ処理の必要性を低減して、MR画像が少ない計算パワーで取得できるようにし得る。さらに、関心被験者のエリアのK空間サンプリングは、高時間効率で行える。K空間サンプリングは息止め状態を検出するとすぐに開始できるからである。それゆえ、関心被験者の息止め状態の長さを全部利用して、MR画像を取得するのに好適なK空間サンプルを取得できる。部分的K空間サンプルの再取得またはMR画像の反復的生成を減らすことができる。上記方法では、関心被験者の動き検出を行って関心被験者のエリアにおける息止め状態を検出するステップは、関心被験者に息止め指令を発するステップの前に開始できる。方法ステップにおけるこれらの差は、息止め状態の検出には影響がない。関心被験者への息止め指令の発行は、イメージングシーケンスが開始された時、MRイメージングシステムにより、またはMRイメージングシステムのオペレータにより、自動的に開始できる。息止め指令は、任意の好ましいコマンドであってもよく、可聴コマンド、視覚的コマンド、または触覚コマンドが含められても良い。したがって、MRイメージングシステムは、息止め指令のコマンドを提供する手段、例えばラウドスピーカシステム、コマンドを示す信号ライト、読める方法でコマンドを表示するディスプレイ装置その他を含み得る。与えられた解像度は、MR画像の解像度を指す。解像度に応じて、必要とされるK空間サンプルの最小数は変化する。関心被験者のエリア自体の動きは最も重要であるが、関心被験者のエリアにおける進行中のまたは次回の動きの表示を与える場合、これはMRイメージングしだいであるが、関心被験者の他のエリアの動きも用いて、関心被験者の息止め状態を検出してもよい。このましくは、動きは関心被験者のエリアで直接検出される。K空間サンプルの再取得後、サンプルを処理して、関心被験者のエリアのMR画像を取得する。2、3または4次元空間の関心被験者のエリアをK空間にマッピングできる。このマッピングに基づいて、K空間データを用いて、関心被験者のエリアのMR画像を取得できる。K空間サンプルの全セットの生成にかかる時間は、一般的には約12ないし18秒であるから、息止めはこのスキャン時間に実現するのは困難であることがある。一般的な息止めの時間は、関心被験者たちが実現できるものであると、約8ないし12秒である。したがって、息止め状態をモニターする時、貴重な時間を無駄にすることなく、K空間データ取得をすぐに開始できる。さらに、K空間データ取得は、息止め状態が実現した後にならないと開始されず、MR画像生成に対して、取得されたK空間サンプルの信頼性が高いようにされる。例えば、息止め指令が発せられた後に、一般的な待ち時間を加えて、K空間データ生成が始まる前に息止めが実現することを保証することにより、息止めの必要性の不要な延長を回避できる。
一般的なアイデアとして、息止め状態が検出される前に、サンプルの1のサブセットの取得を開始する。したがって、方法ステップの順序を、息止め状態の検出前に変更し得る。例えば、関心被験者のエリアのK空間サンプリングを実行してK空間サンプルの第1のサブセットを取得するステップは、息止め指令を発する前に、間に、又は直後に行い得る。一実施形態では、K空間サンプルの第1のサブセットの取得は、息止めコマンドが発行された時に、自動的に開始される。K空間サンプルは、選択されたK空間取得アプローチ/ストラテジに応じて、第1と第2のサブセットにアサインされてもよい。さらに、異なる要件または品質基準、例えば、K空間データの品質基準、アプリオリに決定されたタイミング基準を確定して、K空間サンプルの第1のサブセットの取得が息止め指令の発行前の一定時間に開始できるようにしてもよい。あるいは、第1のサブセットのK空間サンプルの再取得は、息止め指令発行時にすぐに開始することもできる。また、息止め状態への集束を用いて、K空間サンプルの第1のサブセットの取得を開始できる。したがって、第1のK空間サンプルは、第2のサブセットのK空間サンプリングが開始された時に、すでに得られている。K空間サンプルの第1のサブセットのK空間サンプリングにより、息止め状態の検出前に、全エリアのK空間サンプリング、例えば、K空間サンプルの第1と第2のサブセットの取得を効率的な方法で行い得る。息止め状態の時間は関心被験者により維持されており、第2のサブセットのK空間サンプルを取得するために効率的に用いられるからである。第1のサブセットのK空間サンプルは、息止め状態の検出前にすでに取得できている。
K空間サンプルの取得中の動き検出により、取得されたK空間サンプルの動きコンプライアンスの確認が可能である。例えば、十分な動きコンプライアンスが検出されない場合、動きコンプライアンスが再度確立されるまで、取得されたK空間サンプルを破棄でき、またはK空間サンプルの取得を中断できる。動きコンプライアンスは、関心被験者のエリアの動きの状態を指し、これにより、MR画像のMR画像生成に適したK空間サンプルの取得が、高画質かつ低アーティファクトで可能になる。動きコンプライアンスは、K空間サンプルを取得するサンプリングシーケンスに依存し得る。例えば、中心K空間サンプルは、全MR画像に影響を有し、周辺K空間サンプルの動きコンプライアンスよりも、中心K空間サンプルの動きコンプライアンスがより厳しく考慮される。したがって、動きコンプライアンスは、全動きの停止、例えば完全な息止め、または動きが制限された動き状態、例えば低速度の動きであり得る。
好ましい一実施形態では、K空間サンプリングを行って関心被験者のエリアをカバーするK空間サンプルの第2のサブセットを取得するステップは、検出された関心被験者のエリアの動きに基づいて、K空間サンプルの第1の及び/又は第2のサブセットからK空間サンプルを再取得するステップを含む。息止め状態の検出前に第1のサブセットのK空間サンプルについて、十分な動きコンプライアンスが検出されない場合、取得されたK空間サンプルを破棄できる。破棄されたK空間サンプルは、息止め状態の検出後に、再取得でき、またはデータ処理を行って、例えば、取得されたK空間サンプルから、破棄されたK空間サンプルを生成できる。また、第2のサブセットのK空間サンプルの再取得は、動きコンプライアンスに関する上記原理により適用可能である。
好ましい一実施形態では、関心被験者のエリアのK空間サンプリングを行って一解像度で関心被験者のエリアをカバーするK空間サンプルのセットからK空間サンプルの第1のサブセットを取得するステップは、K空間サンプルの第1のサブセットの反復的K空間サンプリングを行うステップを含む。したがって、第1のサブセットのK空間サンプルは、息止め状態が検出されるまで、継続的に再取得される。よって、第1のサブセットのK空間サンプルは、息止め状態に近い関心被験者のエリアの動きで、提供される。これにより、第1のサブセットのK空間サンプルの再取得の品質が改善される。息止め状態に近いので動きが少ないと仮定できるからである。
好ましい一実施形態では、関心被験者のエリアのK空間サンプリングを行ってK空間サンプルの第1のサブセットを取得するステップはK空間の周辺から始めてK空間サンプリングを行いK空間の中心までK空間サンプリングを伸ばすステップを含む。関心被験者のエリアの中心K空間サンプルは、一般的に、関心被験者のMR画像を取得するデータ処理にとって最も重要である。特に、中心K空間サンプルは、MR画像のコントラスト特性を決定するのに重要である。それゆえ、中心K空間サンプルの取得中の動きを回避することは、グロスの動き、ぼやけ、及び信号損失、特に位相の非干渉性により生じるものの効果を回避または低減するために重要である。したがって、K空間サンプルの第1のサブセットの最も重要なK空間サンプルが、息止め状態の検出の直前に取得されたとき、関心被験者のエリアの動きが停止またはほぼ停止し、そのK空間サンプルはMR画像生成のために信頼性が高いと仮定する。
好ましい一実施形態では、K空間サンプリングを行って関心被験者のエリアをカバーするK空間サンプルの第2のサブセットを取得するステップは、K空間サンプリングをK空間の中心から始めて、K空間の周辺まで伸ばすステップを含む。再び、関心被験者のエリアの中心K空間サンプルは、一般的に、関心被験者のMR画像を取得するデータ処理にとって最も重要である。特に、中心K空間サンプルは、MR画像のコントラスト特性を決定するのに重要である。それゆえ、中心K空間サンプルの取得中の動きを回避することは、グロスの動き、ぼやけ、及び信号損失、特に位相の非干渉性により生じるものの効果を回避または低減するために重要である。したがって、K空間サンプルの第2のサブセットの最も重要なK空間サンプルが、息止め状態の検出の直後に取得されたとき、関心被験者のエリアの動きが停止または少なくともほぼ停止し、そのK空間サンプルはMR画像生成のために信頼性が高いと仮定する。所与の解像度で関心被験者のすべてのK空間サンプルを取得するのに必要な時間の間、関心被験者により息止め状態が維持できない可能性もあるので、中心K空間サンプルが最初に取得され、信頼性の高い中心K空間サンプルのみが取得されるようにする。
好ましい一実施形態では、K空間サンプリングを行い、関心被験者のエリアをカバーするK空間サンプルの第2のサブセットを取得するステップは、関心被験者のエリアの動き検出を行うステップと、関心被験者のエリアにおける動きの検出時に、息止め状態が再び検出されるまで、K空間サンプルを破棄するステップとを含む。このステップは、十分な数のK空間サンプルが生成され、取得されたK空間サンプルからMR画像が取得されるまで、繰り返され得る。サンプルシーケンスは、関心被験者のエリアの動きの検出時に任意のポイントで停止でき、K空間サンプルが関心被験者のエリアの動きの検出時に破棄されたサンプルシーケンスの同じステップで継続できる。したがって、たとえ関心被験者が全エリアのK空間サンプルを提供するのに必要な時間だけ動き停止を維持できなくても、動き停止の検出に基づいて中断されたK空間サンプリングは、関心被験者のエリアの全K空間をカバーするサンプルのセットを提供できる。サンプルは、動きがMR画像中にアーティファクトを生じることなしに提供できる。動きのない好ましいK空間サンプルのみを画像生成に用いることができるからである。特に息止めの場合、前の息止めに対して良い一貫性を示すほぼ同じ位置の胸部のデータのみを受け入れ、K空間サンプリングの中断が、低画質の画像生成とならないようにする。この方法により、動きを止められない、すなわち必要なK空間サンプルの取得に十分な時間だけ息止め状態を維持できない関心被験者たちのMR画像生成も可能になる。このましくは、動きの検出時、他の息止め指令を関心被験者に発行する。
好ましい一実施形態では、関心被験者の動き検出を行って関心被験者のエリアにおける息止め状態を検出するステップは、呼吸ナビゲータ信号をモニターするステップを含む。呼吸ナビゲータ信号は、1D、2D、または3D画像の画像ベースの取得を指す。1D呼吸ナビゲータ信号は、「ペンシルビーム」とも呼ばれ、関心被験者のエリアにおける関心被験者の1Dプロファイルを取得する。呼吸ナビゲータ信号を用いて、身長方向の関心被験者の横隔膜をモニターし、その横隔膜の動きによる呼吸運動を特定する。呼吸ナビゲータ信号が安定している時は、関心被験者が呼吸を止めた、すなわち息止めを行ったものと仮定する。呼吸ナビゲータ信号に基づく動き検出は、追加的ハードウェアが必要ナノので、ハードウェア修正の必要無しに、既存のMRイメージングシステムに容易に実装できる。呼吸ナビゲータ信号は、適応的イメージング方法を用いる一般的なMRイメージングシステムによりモニターできる。
好ましい一実施形態では、関心被験者の動き検出を行い関心被験者のエリアにおける息止め状態を検出するステップは、少なくとも1つのK空間ナビゲータを確定するステップと、少なくとも1つのK空間ナビゲータのリアルタイムモニタリングを行うことにより、関心被験者のエリアにおける関心被験者の動き検出を行うステップとを含む。K空間ナビゲータ手法の一例は、非類似性マトリックスを用いてK空間データの一貫性をチェックするいわゆるオービタルナビゲータ(O−NAV)である。これにより、関心被験者のリアルタイム動き検出の実行が可能になる。K空間ナビゲータは、MR画像のK空間サンプルの取得に加えて行われる、インターリービングK空間データ取得ステップとして実装してもよい。K空間ナビゲータは、等距離でインターリーブされ、一連の複数のK空間取得はスピンエコーまたはグラジエントエコー測定値(readoouts)またはショット(shots)のセグメントにグループ分けされる。それにもかかわらず、画像及びK空間ナビゲータは、非等距離でインターリーブされることもできる。インターリーブされたK空間データ取得の周波数は、条件が違えば異なっても良い。例えば、動き停止コマンドを発した後、K空間ナビゲータのモニタリングを増やして、いつ息止め状態に到達したか精密に検出できるようにすれば有用である。一方、息止め状態が検出された後は、動き検出頻度を減らすこともできる。あるいは、K空間ナビゲータデータは、サンプルが反復的に生成されているとき、MR画像生成に用いられる固有の(inherent)イメージングデータにすなわち、MR画像のK空間サンプルに対応する。これは、息止め状態が検出されるまでK空間サンプルの第1のサブセットが繰り返し取得されるとき、息止め状態の検出の前にK空間サンプルの第1のサブセットの取得を指してもよい。一般的に、複数のK空間ナビゲータが、動き検出のために実装される。このましくは、円形K空間ナビゲータが用いられる。円形K空間ナビゲータは、K空間の中心の周りに円形に配置されてもよい。K空間ナビゲータに基づく動き検出の実装は、追加的ハードウェアが必要ナノので、ハードウェア修正の必要無しに、既存のMRイメージングシステムに容易に実装できる。K空間ナビゲータを用いて、適応されたイメージング方法を用いる一般的なMRイメージングシステムにより関心被験者の動きをモニターできる。
好ましい一実施形態では、関心被験者の動き検出を行って関心被験者のエリアにおける息止め状態を検出するステップは、K空間サンプルの第1のサブセットのK空間サンプルを繰り返し処理して、関心被験者のエリアの画像を取得するステップと、取得された画像から関心被験者の動き検出を行うステップとを含む。画像における動きの検出の既知の手法を用いることができる。K空間サンプルの繰り返し処理は、K空間サンプルが繰り返し取得されることを示唆する。動き検出のこの実装は、高速画像処理を要する。K空間サンプルに基づく画像生成に典型的な時間は約50msだからである。K空間サンプリングに基づく画像に基づく動き検出の実装は、追加的ハードウェアが必要ナノので、ハードウェア修正の必要無しに、既存のMRイメージングシステムに容易に実装できる。K空間サンプリングに基づく画像を用いて、適応されたイメージング方法を用いる一般的なMRイメージングシステムにより関心被験者の動きをモニターできる。取得された画像に基づく動き検出を同様に適用して、K空間サンプルの第2のサブセットのK空間サンプリングを行うことができる。
好ましい一実施形態では、関心被験者の動き検出を行って関心被験者のエリアにおける息止め状態を検出するステップは、K空間サンプルの第1のサブセットに基づいてK空間サンプルの生成を行うステップと、生成された前記サンプルと、K空間内の対応する取得されたサンプルとの一貫性を評価するステップとを含む。K空間の生成は、例えばCOCOAアプローチ(Data Convolution and Combination Operation(COCOA) for Motion Ghost Artifacts Reduction)などの既知の手法に基づいて行われ得る。この手法は、K空間データの後処理を行って、すでに取得したK空間サンプルに基づいてアーティファクトを低減するのに現在用いられている。それにもかかわらず、K空間サンプルの全量を取得する前に、K空間生成を用いて、関心被験者が動いているか否か推定することができる。このましくは、K空間カーネルエリアからのサンプルを、K空間サンプルの生成に用いる。取得されたK空間サンプルに基づいて、隣接K空間サンプルを推定できる。生成されたK空間サンプルを含む後続のK空間サンプル間の偏差が閾値より大きいとき、これは動きのインジケータである、すなわち、関心被験者はまだ動くのをやめておらず、特に呼吸を止めていない。K空間生成に基づく動き検出の実装は、追加的ハードウェアが必要ナノので、ハードウェア修正の必要無しに、既存のMRイメージングシステムに容易に実装できる。K空間生成からのデータを用いて、適応されたイメージング方法を用いる一般的なMRイメージングシステムにより関心被験者の動きをモニターできる。さらに、K空間生成を実行しても、既存のMRイメージングシステムへの影響は非常に小さい。それにもかかわらず、この実装には高速データ処理が必要である。このアプローチは、取得されたK空間サンプルを置き換えるとの利点もある。動き停止の質は、一般的には息止め状態の最初の4秒間に高くなるので、その後に取得されるK空間サンプルは、一般的にそれ以前のK空間サンプルよりも安定である。COCOAを用いて、以前に取得されたK空間サンプルのためにK空間サンプルを生成できる。生成されたK空間サンプルは、取得されたK空間サンプルと比較される。比較時、先に取得されたK空間サンプルは、生成または再取得されたデータにより置き換えられてもよい。
好ましい一実施形態では、本方法は、関心被験者のエリアにおける関心被験者の動きをモニターする動き検出装置を設ける追加ステップを含み、関心被験者の動き検出を行って関心被験者のエリアにおける息止め状態を検出するステップは、動き検出装置を用いて息止め検出を行うステップを含む。対応するMRイメージングシステムは、関心被験者のエリアにおける関心被験者の動きをモニターする動き検出装置を含み、MRイメージングシステムは、動き検出装置を用いて関心被験者のエリアにおける関心被験者の動き検出を行うように適応されている。異なる種類の動き検出装置を用いても良い。動き検出装置の使用は、具体的な関心被験者のエリアに依存し得る。一般的に、例えばカメラに基づく光学的動き検出装置を含む、異なる種類の動き検出装置が好ましい。動き検出装置は、このましくはMRイメージングシステムの制御部に接続され、動き検出装置を用いて行われる動き検出に基づいてK空間サンプリングを制御できるようになっている。動き検出装置は、例えば、MRイメージングシステムの制御部に接続された任意種類のセンサを用いて実装できる。動きを検出する信号処理は、例えば、そのセンサからのセンサ入力に基づいて制御部により行われる。動き検出装置を用いた動き検出の実装は、既存のMRイメージングシステムにおいて容易に実装され得る。動き検出に要する追加的計算パワーは非常に少ないからである。したがって、動き検出装置は、一般的なMRイメージングシステムであって、その動き検出装置に基づく動き検出に適応したイメージング方法を用いるものと共に容易に利用することができる。動き検出装置は、可視またはTHz周波数帯域の光センサ、または画像生成に用いられる超音波装置を含んでいてもよい。さらに好ましくは、動き検出装置は、ビデオ生成をして、生成されたビデオ信号に基づいて動き検出を行うように適応されていてもよい。動き検出装置の利用の一利点は、それが取得と干渉せず、余分なスキャン時間を必要としないことである。
好ましい一実施形態では、関心被験者のエリアにおける関心被験者の動きをモニターする動き検出装置を設けるステップは、呼吸ベルトを設けるステップを含み、関心被験者のエリアにおける関心被験者の動き検出を行うステップは、呼吸ベルトを用いて動き検出を行うステップを含む。呼吸ベルトは、特に、呼吸検出に適している。呼吸ベルトは、関心被験者の胸部に直接巻いて、関心被験者の呼吸と息止めを検出する信頼できる信号を提供することができる。呼吸ベルトは、このましくは空気式呼吸ベルトであり、MR環境で容易に利用することができるものである。空気式動き検出は、磁場及び/又はRF場に敏感なコンポーネント無しで設けることができる。
好ましい一実施形態では、関心被験者の動き検出を行って関心被験者のエリアにおける息止め状態を検出するステップは、関心被験者のエリアにおける動きをモニターして、息止め状態への集束を検出するステップを含む。息止め状態への集束(convergence)、例えば動きがほぼ停止した状態は、K空間サンプルの第1の及び/又は第2のサブセットのK空間サンプルの取得が開始できる表示を提供できる。動きの停止は、予想できる方法で生じ、一般的には短時間で生じると思われるからである。息止め状態への動き集束は、例えば、現時点における動きの速さを推定することにより、実現できる。この動き集束を用いて、取得、またはK空間サンプルの再取得、または最初のスキャンフェーズにおける周辺K空間サンプルのサンプリングを最適化できる。さらに、中心K空間サンプルの取得または再取得は、動きが息止め状態に集束する場合にのみ行い得る。速度は、動きをモニターする頻度の、考慮中の動き状態差分を比較することにより得られる。このましくは、第1のサブセットのK空間サンプルは、例えば、息止め状態の検出前の1と1.5秒タイミング基準内であれば受け入れられる。したがって、K空間サンプルの第2のサブセットは、画像再構成に対して適当なサンプルで生成できる。息止め状態の前の動きは、すでに減少しており、画像生成への影響は小さいからである。
好ましい一実施形態では、関心被験者のエリアにおける動きをモニターして息止め状態への集束を検出するステップは、関心被験者のエリアの動きの速さをモニターするステップを含み、関心被験者のエリアにおける息止め状態の検出をするステップは、関心被験者のエリアにおける動きの速さが誤差も含め0に近くなることの検出として、息止め状態を検出するステップを含む。これにより、息止め状態の表示として速さ自体を用いることができる。息止めが完全でなくても息止め状態の検出をできるようにするため、誤差は必要である。息止めをしているときであっても、関心被験者は何らかの動きを維持しているからである。速度は、動きをモニターする頻度の、考慮中の動き状態差分を比較することにより得られる。
好ましい一実施形態では、本方法は、関心被験者の息止め状態の動き検出及び/又は集束に応じて、関心被験者のエリアのK空間サンプリングを実行するサンプリングストラテジに適応する追加的ステップとを有する。よって、K空間サンプルは、画像生成の意味で最も適当な時に、高画質及び低アーティファクトで取得できる。一例として、MR画質に大きな影響がある領域のK空間サンプルは、息止め状態で取得され、他のサンプルは、息止め状態が満たされていない時に取得できる。さらに、動き集束情報を用いて、いつ息止め状態でK空間サンプルを取得できるか事前に推定して、息止め状態の予見を可能にする。これにより、サンプリングストラテジの事前適応ができ、K空間サンプルの取得を効率的な方法で行えるようになる。このましくは、K空間サンプリングストラテジのリアルタイム適応を、現在の動き状態、又は息止め状態へ向けた動きの集束に基づいて行う。サンプリングストラテジの適応は、K空間サンプルの全セット、及びK空間サンプルの第1と第2のサブセットを指す。サンプリングストラテジの適応は、K空間サンプルの第1の及び/又は第2のサブセットが中心K空間サンプル及び周辺K空間サンプルを含む時に、役に立ち、中心K空間サンプルが最も好ましい動き条件において取得でき、すなわち中心K空間サンプルが息止め状態で取得され、息止め状態が例えば十分長い時間達成されなくても、周辺K空間サンプルがこのましくは息止め状態で必要とされる。
好ましい一実施形態では、関心被験者のエリアのK空間サンプリングを実行するステップは、分散サンプリング方式でK空間サンプリングを行うステップを含み、取得されたK空間サンプルを処理するステップは、分散サンプリングスキームに基づいて取得されたK空間サンプルに基づいてK空間中の穴を再構成するステップを含む。このましくは、サンプリング方式はポアソンディスクサンプリング方式である。K空間におけるK空間サンプルの分布により、少ないK空間サンプルで信頼できるMR画像生成を可能にする。さらに、K空間サンプルの再構成は、例えば圧縮センシングアルゴリズムにより容易に行い得る。この状況では、反復的センスアルゴリズムを用いることもできる。
本技術分野の当業者には言うまでもないが、本発明の態様は、システム、方法、またはコンピュータプログラム製品として実施され得る。したがって、本発明の態様は、完全にハードウェアによる実施形態、(ファームウェア、常駐ソフトウェア、マイクロコードなどを含む)完全にソフトウェアによる実施形態、またはソフトウェアとハードウェアの態様を組み合わせた実施形態の形式を取り得る。これらはすべてここでは概して「回路」、「モジュール」または「システム」と呼ぶ。さらに、本発明の態様は、コンピュータ読み取り可能プログラムコードが記録された一以上のコンピュータ読み取り可能媒体に化体されたコンピュータプログラムの形式を取っても良い。
これらのコンピュータプログラム命令は、コンピュータ、他のプログラマブルデータ処理装置、又は他のデバイスに、フローチャート及び/又はブロック図のブロックに指定された機能/動作を実装する命令を含む、コンピュータ読み取り可能媒体に記憶された命令が製品となるように機能するように指示できるコンピュータ読み取り可能媒体に記憶されてもよい。
コンピュータプログラム命令は、コンピュータ、その他のプログラマブルデータ処理装置、又はその他のデバイスにロードされ、一連の動作ステップがそのコンピュータ、その他のプログラマブルデータ処理装置、又はその他のデバイスで実行され、コンピュータ実装プロセスとなり、コンピュータまたはその他のプログラマブルデータ処理装置上で実行された命令がフローチャート及び/又はブロック図のブロックに指定された機能/動作を実装するプロセスを提供するようにする。
本発明の上記その他の態様を、以下に説明する実施形態を参照して明らかにし、説明する。しかし、かかる実施形態は必ずしも本発明の全範囲を表すものではなく、それゆえ特許請求の範囲を参照して、本発明の範囲を解釈する。
図中、
呼吸ベルトを含む、本発明の好ましい実施形態による磁気共鳴イメージングシステムの一般的なセットアップを示す図である。 好ましい実施形態による、磁気共鳴画像生成の方法を示す図である。 1回、2回、及び3回の繰り返し後のサンプリングポイントの密度と対応サンプリングパターンとを示す図である。 5回、10回、及び15回の繰り返し後のサンプリングポイントの密度と対応サンプリングパターンとを示す図である。 図3及び図4によるサンプリングポイントの密度と、タイミンググラフ及び呼吸グラフを示す図である。
図1は、MRスキャナー112を有する磁気共鳴(MR)イメージングシステム110の一実施形態の一部を示す模式図である。MRイメージングシステム110は、静磁場を発生するように設けられた主磁石114を含む。主磁石114は中心ボア(centralbore)を有する。これは、その中に配置される通常は患者である関心被験者120の中心軸118の周りに検査空間116を提供する。詳しくは図示していないが、主磁石114は、2セットのコイル、1セットの内側コイル及び1セットの外側コイルを有し、これらはそれぞれ主磁石114の内側部分と外側部分とに円筒状に構成され、主磁石114のz軸に沿って離間されている。さらに、MRイメージングシステム110は、静磁場に重畳されるグラジエント磁場を発生するために設けられたグラジエント磁場コイル系122を含む。主磁石122は、本技術分野では既知にように、主磁石114のボア内に同心状に構成されている。
さらに、MRイメージングシステム110は、チューブ状のボディを有する全身コイルとしてデザインされた無線周波数(RF)アンテナ装置140を含む。RF送信フェーズ中に検査空間116にRF磁場を印加して、関心被験者120の核を励起するRFアンテナ装置140が設けられる。RFアンテナ装置140は、RF受信フェーズ中に、励起された核からMR信号を受信するためにも設けられる。MRイメージングシステム110の動作状態では、RF送信フェーズとRF受信フェーズは連続して行われる。RFアンテナ装置140は、主磁石114のボア内に同心状に構成される。本技術分野では既知であるように、円筒形金属RFスクリーン124が、主磁石122とRFアンテナ装置140との間に同心状に構成される。
さらに、MRイメージングシステム110は、本技術分野では周知であるように、取得されたMR信号からを再構成するために設けられたMR画像再構成部130と、MRスキャナー112の機能を制御するために設けられた、モニター部128を有するMRイメージングシステム制御部126とを含む。制御ライン132はMRイメージングシステム制御部126とRF送信部134との間にインストールされ、RF送信フェーズ中にRFスイッチ部136を介してRFアンテナ装置140にMR無線周波数のRFパワーを供給(feed)するように設けられている。次いで、RFスイッチ部136は、MRイメージングシステム制御部126によっても制御され、他の制御ライン138が、MRイメージングシステム制御部126とRFスイッチ部136との間にインストールされ、その目的を果たす。RF受信フェーズ中、RFスイッチ部136は、前置増幅後、MR信号をRFアンテナ装置140からMR画像再構成部130に送る。
ここで、図2を参照して、関心被験者120のエリア144のMRイメージングをする方法を、説明する実施形態によるMRイメージングシステム110を参照して説明する。本方法はステップS100において始まり、関心被験者120がMRイメージングシステム110の検査空間116に配置される。本実施形態における関心被験者120のエリア144は、関心被験者120の胸部である。さらに、呼吸を止める息止め指令が、音響コマンドすなわち発話コマンドとして、MRイメージングシステム110のオペレータから関心被験者120に発せられる。別の一実施形態では、例えば制御部126上で実行されている制御SWが、ラウドスピーカによる音響シーケンスの再生として、関心被験者120に息止め指令を発しても良い。
ステップS110において、MRイメージングシステム110は、関心被験者120の胸部144をカバーするK空間154サンプルの第1のサブセットの反復的K空間154サンプリングの実行を開始する。K空間154サンプルの第1のサブセットは主に周辺K空間154サンプルを含む。さらに、K空間154サンプルの第1のサブセットのK空間154サンプリングは、K空間154の周辺152で開始され、K空間154の中心150までK空間154サンプリングが伸ばされる。第1のサブセットのK空間154サンプルは、選択されたK空間154取得アプローチ/ストラテジに応じた解像度で、関心被験者120の胸部144をカバーするK空間154サンプルから選択される。この解像度は、MR画像の解像度を指す。解像度に応じて、必要とされるK空間154サンプルの最小数は変化する。別の一実施形態では、K空間154サンプルの第1のサブセットの取得は、息止め指令が発せられる前からすでに開始されている。
ステップS120において、MRイメージングシステム110は、関心被験者120の胸部144の動き検出を行う。ステップS120の動き検出の詳細は、後で説明する。別の一実施形態では、動き検出を行うステップS120は、関心被験者120に息止め指令を発するステップS100の前に開始される。
ステップS130において、MRイメージングシステム110は息止め状態をチェックする。したがって、MRイメージングシステム110関心被験者120が呼吸を止め、安定した息止めが確立されたか確認する。息止め状態への集束を検出すると、動き検出が実施される。息止め状態への集束(convergence)は、息止め状態の前に何時息止め状態が始まったか時間に関する情報を提供し、現時点における関心被験者120の胸部144の動き速度を用いてモニターされる。速度は、動きをモニターする頻度の、考慮中の運き状態差分を比較することにより得られる。したがって、息止め状態への集束は、関心被験者120の胸部144における動きの速度に基づいて検出され、それにより速度が0に近いことは息止め状態を示す。動きが検出された場合、方法はステップS110に戻り、K空間154サンプルの第1のサブセットのK空間154サンプルが繰り返し取得され、直近のK空間154サンプルが前のK空間154サンプルを置き換える。そうでなければ、本方法はステップS140に続く。別の一実施形態では、動き収束(motion convergence)を用いて、K空間154サンプルの第1のサブセットの取得または再取得を適応させる。
ステップS140において、MRイメージングシステム110は、K空間154サンプルの第1のサブセットからK空間154サンプルの第2のサブセットに、K空間154サンプルの取得を切り替え、関心被験者120の胸部144のK空間154サンプリングを完了している。したがって、K空間154サンプルの第2のサブセットは、まだK空間154サンプルの第1のサブセットの一部ではない、与えられた解像度、で関心被験者120の胸部144をカバーするK空間154サンプルを含む。K空間154サンプルの第2のサブセットの取得は、息止め状態の検出時にすぐに開始され、第1の中心K空間154サンプルが取得され、K空間154サンプルの第2のサブセットのK空間154サンプリングが周辺K空間154サンプルまで伸ばされる。こうして、の関心被験者120胸部144に関する全K空間154がカバーされる。
K空間154サンプルの第2のサブセットのK空間154サンプルの取得が行われている間に、ステップS150において、動き検出が行われる。後で詳しく説明するように、ステップS120において適用可能な原理は、ステップS150にも適用される。別の一実施形態では、関心被験者120の胸部144のK空間154サンプリングを行うサンプリングストラテジは、息止め状態の動き検出及び/又は集束に応じて、適応される。したがって、息止め状態が安定的に達成されると、中心K空間154サンプルが取得され、一方、他の場合には、周辺K空間154サンプルがその後のステップS140において取得される。
ステップS160において、MRイメージングシステム110は、K空間154サンプリングがMR画像生成のために十分なK空間154サンプルを提供したか確認し、すなわち、K空間154サンプルの第1と第2のサブセットの所得されたK空間154サンプルが、与えられた解像度のMR画像生成にとって十分か確認する。これは、関心被験者120の胸部144の検出された動きに基づいて破棄されたK空間154サンプルが再取得されたか確認することを含む。後のステップS180と比較されたい。したがって、K空間154サンプリングが終わると、本方法は、ステップS190におけるMR画像生成に続く。そうでなければ、本方法はステップS140に続く。あるいはステップS160は、図2の異なるポイントで行え、例えば動き検出ステップS150の前、またはステップS170の後に行える。これは後で説明する。
ステップS170において、MRイメージングシステム110は、関心被験者120が呼吸を止めて息止め状態が達成されたか、確認する。動きまたは不十分な動きコンプライアンス(compliance)が検出されると、本方法は、ステップS180において、最近のK空間154サンプルを破棄し、前のサンプリングステップによるK空間154サンプルの第2のサブセットのK空間154サンプルを取得するステップS140に続く。これは、破棄されたK空間154サンプルを置き換える有効なK空間154サンプルを取得するためである。したがって、サンプルシーケンスは、関心被験者120の胸部144の動きの任意の検出時点で停止され、息止め状態の検出のサンプルシーケンスの同じステップに続く。K空間154サンプルの取得中の動き検出により、取得されたK空間154サンプルの動きコンプライアンスの確認ができる。動きコンプライアンス(動きコンプライアンス)は、関心被験者のエリアの動きの状態を指す。これにより、MR画像のMR画像生成に適したK空間154サンプルの取得を可能とする。そうでなければ、本方法は、胸部144のさらなるK空間154サンプルを取得するため、ステップS140に進み、K空間154サンプリングする。
ステップS190において、MRイメージングシステム110は、K空間154サンプルの第1と第2のサブセットのK空間154サンプルを処理して、関心被験者120の胸部144のMR画像を取得する。K空間154データは、異なる位相エンコーディングステップで取得され、これは分散サンプリングスキーム(この実施形態ではポアソンディスクサンプリングスキームである)により多数のサンプルを収集することにより、デジタイズされたものである。K空間154サンプルの第1と第2のサブセットは共に、与えられた解像度で関心被験者120の胸部144をカバーする。K空間154サンプルの第1と第2のサブセットを含むK空間154サンプルのセットは、処理され、例えばフーリエ変換により、MR画像に変換される。このステップでは、K空間154サンプルは、K空間154サンプルが取得された時に、関心被験者120の胸部144の検出された動きに応じて、拒絶及び/又は加重される。よって、K空間154サンプルの、関心被験者120の胸部144のMR画像へのマッピングが行われ、関心被験者120の胸部144のMR画像を取得する。さらに、分散サンプリングスキームに基づいて、取得されたK空間154サンプルに基づいたK空間154の穴(holes)の再構成が、圧縮センシングアルゴリズムにより行われる。
図3と図4は、例として、1、2、3、5、10及び15回反復後のサンプリングポイントの密度と、対応するサンプリングパターンとを示す。これは、K空間154サンプルの第2のサブセットのK空間154サンプルの取得が、中心K空間154サンプルから周辺K空間154サンプルまで伸びるとき、ステップS140において行われる。
図5は、さらに各グラフにより示されたタイミングと呼吸とを示す。図3A、3B、3C、4A、4B及び4Cは、サンプリングポイントの密度を示す。図3A’、3B’、3C’、4A’、4B’及び4C’は、K空間154中のサンプリングパターンを示す。サンプリングされたポイントの密度の分布は、2つの条件の結果である。空間的解像度は時間とともに増加し、サンプリング密度は周辺K空間154に向けて単調に減少する。最大空間的解像度はオペレータにより確定される。K空間154の中心150に近い最大サンプリング密度と、K空間154の周辺152に向かう減衰率との両方は、取得の間に適応され、一方、平均減衰係数は成功裏に減少する。サンプリングパターンを生成するため、上記の通り、ポアソンディスクを用いる。ポアソンディスクサンプリングは、隣接するサンプルは所定距離未満になってはならないとの付加的制約を有するサンプルのランダム分布である。上記により、息止め指令を発すると、呼吸が止まり息止め状態が実現される。その時、中心K空間154サンプルが最初に取得され、呼気が始まる時、すなわち息止め状態が終わる時に向けて、周辺K空間154サンプルは後で含められる。図から分かるように、サンプリング密度は、K空間154のエッジでゼロに下がる。複数のK空間154サンプルのランダムな疑似ランダム選択が周辺K空間154にあると、K空間154の周辺152のほぼローカルな一様なサブサンプリングとなる。
以下、本発明の異なる実施形態による、動き検出のステップS120を詳細に説明する。第1の実施形態では、ステップS120において、呼吸ナビゲータ信号がモニターされる。呼吸ナビゲータ信号は、この実施形態における1D画像の画像ベースの取得を指す。1D呼吸ナビゲータ信号は「ペンシルビーム」とも呼ばれる。呼吸ナビゲータ信号は、身長方向の関心被験者120の横隔膜をモニターし、その横隔膜の動きによる呼吸運動を特定する。呼吸ナビゲータ信号が安定している時は、関心被験者120が呼吸を止めた、すなわち息止めを行ったものと仮定する。呼吸ナビゲータ信号は、ハードウェア修正無しに、MRイメージングシステム110によりモニターされる。
第2の実施形態では、ステップS120において複数のK空間154ナビゲータが規定される。関心被験者120の胸部144における動き検出は、K空間154ナビゲータデータのモニタリングにより行われる。この実施形態におけるK空間154ナビゲータは、K空間154の中心150の周りに円周上に配置される。K空間154ナビゲータは、K空間154サンプルの第1のサブセットのK空間154サンプルの取得とインターリーブされ、この実施形態においては、オービタルナビゲータ(O−NAV)と非類似行列を用いて、収集されたデータのK空間154一貫性をチェックする。K空間154ナビゲータは、インターリーブK空間154データ取得ステップとして実装され、K空間154サンプルの第1のサブセットのK空間154サンプルの取得に加えて行われる。K空間154ナビゲータは、一連の複数のK空間154取得、すなわち複数のK空間154サンプルの取得と等距離にインターリーブされる。これらは、スピンエコー(TSE)またはグラジエントエコー(TFE)読み出しまたはショットのセグメントにグループ化される。K空間154ナビゲータは、ハードウェア修正無しに、MRイメージングシステム110によりモニターされる。
第3の実施形態では、ステップS120において、動き検出は、K空間154サンプルの第1のサブセットのK空間154を繰り返し処理し、サンプル 関心被験者120の胸部144の画像を取得することを含む。胸部144の動き検出は、取得された画像により行われる。画像における動きの検出の既知の手法を用いることができる。取得された画像に基づく動き検出は、ハードウェア修正無しに、MRイメージングシステム110により行われる。
第4の実施形態では、ステップS120において、関心被験者120の胸部144の動き検出は、K空間154サンプルの第1のサブセットと、生成されたサンプルの評価一貫性と、K空間154内の対応する取得されたサンプルとに基づいて、K空間154サンプルの生成により行われる。K空間154の生成は、COCOAアプローチ(Data Convolution and Combination Operation(COCOA) for Motion Ghost Artifacts Reduction)に基づき行われる。この実施形態における、K空間154サンプルの生成のためK空間154カーネルエリアからのサンプルが用いられる。取得されたK空間154サンプルに基づいて、隣接するK空間154サンプルが推定される。生成されたK空間154サンプルを含む後続のK空間154サンプル間の偏差が閾値より大きいとき、これは動きすなわち関心被験者120が呼吸していることのインジケータである。K空間154生成に基づく動き検出は、ハードウェア修正無しに、MRイメージングシステム110により行われる。
上記の第1ないし第4の実施形態の場合、ステップS150において、適宜動き検出が行われる。例えばK空間154サンプルの(取得された)異なるサブセットによる差は、適宜処理される。
第5の実施形態では、胸部144における関心被験者の動き120をモニターする動き検出装置としてが呼吸ベルト142設けられる。呼吸ベルト142は、空気圧式呼吸ベルトであり、関心被験者120の胸部144の周りに取り付けられ、関心被験者120の胸部144の動きを示す圧力信号を提供する。呼吸ベルト142は、MRイメージングシステム110の制御部126に接続され、制御部126は、呼吸ベルト142を用いて行われる動き検出を用いてMR画像生成の方法を実行するように適応される。したがって、ステップS120において、動き検出は動き検出装置142を用いて行われる。
一実施形態では、MRイメージングシステムの制御部126に上記の方法を実行させる一組の命令を含む、ストレージ装置、フロッピーディスク、コンパクトディスク、CD、デジタルバーサタイルディスク、DVD、ブルーレイディスク、またはランダムアクセスメモリ、RAMなどのコンピュータ読み取り可能媒体が提供される。
したがって、コンピュータ利用可能プログラムコードであって上記の方法を実行するよう適応されたものを含むコンピュータ利用可能媒体を含むコンピュータプログラム製品が提供される。
さらに別の一実施形態では、MRイメージングシステム110を更新するソフトウェアパッケージであって、MRイメージングシステム10を制御して上記方法を実行する命令を含むものが提供される。
本発明を、図面と上記の説明に詳しく示し説明したが、かかる例示と説明は例であり限定ではなく、本発明は開示した実施形態には限定されない。請求項に記載した発明を実施する際、図面、本開示、及び添付した特許請求の範囲を研究して、開示した実施形態のその他のバリエーションを、当業者は理解して実施することができるであろう。 請求項において、「有する(comprising)」という用語は他の要素やステップを排除するものではなく、「1つの(「a」又は「an」)」という表現は複数ある場合を排除するものではない。相異なる従属クレームに手段が記載されているからといって、その手段を組み合わせて有利に使用することができないということではない。請求項に含まれる参照符号は、その請求項の範囲を限定するものと解してはならない。
110 磁気共鳴(MR)イメージングシステム
112 磁気共鳴(MR)スキャナー
114 主磁石
116 RF検査空間
118 中心軸
120 関心被験者
122 グラジエント磁場コイルシステム
124 RFスクリーン
126 MRイメージングシステム制御部
128 モニター部
130 MR画像再構成部
132 制御ライン
134 RF送信部
136 RFスイッチ部
138 制御ライン
140 無線周波数(RF)アンテナ装置
142 呼吸ベルト、動き検出装置
144 胸部、関心被験者のエリア
150 中心K空間
152 外部K空間、周辺K空間
154 K空間

Claims (13)

  1. 関心被験者のエリアの磁気共鳴(MR)イメージングをする方法であって、
    前記関心被験者に息止め指令を発するステップと、
    前記関心被験者の動き検出を行って、前記関心被験者のエリアにおける息止め状態を検出するステップと、
    前記関心被験者のエリアのK空間サンプリングを行って、前記関心被験者のエリアにおける息止め状態の検出の前に、前記関心被験者のエリアをカバーしている前記K空間サンプルのセットからK空間サンプルの第1のサブセットを取得するステップと、
    前記関心被験者のエリアにおける息止め状態を検出したとき、一解像度で前記関心被験者のエリアのK空間サンプリングを行うステップであって、K空間サンプリングを行って前記関心被験者のエリアをカバーするK空間サンプルの第2のサブセットを取得することを含むステップと、
    前記関心被験者のエリアをカバーするK空間サンプルを処理して前記関心被験者のエリアのMR画像を取得するステップであって、前記K空間サンプルの第1と第2のサブセットを処理して前記関心被験者のエリアの前記MR画像を取得することを含むステップとを含み、
    前記K空間サンプルの第1と第2のサブセットは共に前記一解像度で前記関心被験者のエリアカバーし、
    前記関心被験者の動き検出を行って前記関心被験者のエリアにおける息止め状態を検出するステップは、K空間サンプルの前記第1の及び/又は第2のサブセットの取得中に少なくとも部分的に行われ、
    前記関心被験者のエリアをカバーするK空間サンプルを処理して前記関心被験者のエリアのMR画像を取得するステップはさらに、検出された前記関心被験者のエリアの動きに応じてK空間サンプルの前記第1の及び/又は前記第2のサブセットのうちのK空間サンプルを排除及び/又は加重するステップを含む、
    磁気共鳴イメージング方法。
  2. K空間サンプリングを行って、前記関心被験者のエリアをカバーするK空間サンプルの第2のサブセットを取得するステップは、検出された前記関心被験者のエリアの動きに応じてK空間サンプルの前記第1の及び/又は前記第2のサブセットからK空間サンプルを再取得するステップを含む、
    請求項1に記載の磁気共鳴イメージング方法。
  3. 前記関心被験者のエリアのK空間サンプリングを行って一解像度で前記関心被験者のエリアをカバーする前記K空間サンプルのセットからK空間サンプルの前記第1のサブセットを取得するステップは、K空間サンプルの前記第1のサブセットの反復的K空間サンプリングを行うステップを含む、
    請求項1または2に記載の磁気共鳴イメージング方法。
  4. 前記関心被験者のエリアのK空間サンプリングを行ってK空間サンプルの第1のサブセットを取得するステップは前記K空間の周辺から始めてK空間サンプリングを行い前記K空間の中心までK空間サンプリングを伸ばすステップを含む、。
    請求項1ないし3いずれか一項に記載の磁気共鳴イメージング方法。
  5. K空間サンプリングを行って前記関心被験者のエリアをカバーするK空間サンプルの第2のサブセットを取得するステップは、K空間サンプリングを前記K空間の中心から始めて、前記K空間の周辺まで伸ばすステップを含む、
    請求項1ないし4いずれか一項に記載の磁気共鳴イメージング方法。
  6. 前記関心被験者の動き検出を行って前記関心被験者のエリアにおける息止め状態を検出するステップは、
    K空間サンプルの前記第1のサブセットに基づいてK空間サンプルの生成を行うステップと、
    生成された前記サンプルと、前記K空間内の対応する取得されたサンプルとの一貫性を評価するステップとを含む、
    請求項1ないし5いずれか一項に記載の磁気共鳴イメージング方法。
  7. 前記関心被験者のエリアにおける前記関心被験者の動きをモニターする動き検出装置を提供するステップをさらに含み、
    前記関心被験者の動き検出を行って前記関心被験者のエリアにおける息止め状態を検出するステップは、前記動き検出装置を用いて息止め検出を行うステップを含む、。
    請求項1ないし6いずれか一項に記載の磁気共鳴イメージング方法。
  8. 前記関心被験者の動き検出を行って前記関心被験者のエリアにおける息止め状態を検出するステップは、前記関心被験者のエリアにおける動きをモニターして前記息止め状態への集束(convergence)を検出するステップを含む、。
    請求項1ないし7いずれか一項に記載の磁気共鳴イメージング方法。
  9. 関心被験者のエリアにおける動きをモニターして前記息止め状態への集束を検出するステップは、前記関心被験者のエリアにおける運動の速さをモニターするステップを含み、
    前記関心被験者のエリアにおける息止め状態を検出するステップは、前記関心被験者のエリアにおける動きの速さが誤差範囲で0に近いことを検出するステップを含む、。
    請求項8に記載の磁気共鳴イメージング方法。
  10. 前記関心被験者の前記息止め状態の前記動き検出及び/又は集束に応じて、前記関心被験者のエリアのK空間サンプリングを実行するサンプリングストラテジに適応するステップとを有する、
    請求項1ないし9いずれか一項に記載の磁気共鳴イメージング方法。
  11. 前記関心被験者のエリアのK空間サンプリングを実行するステップは、分散サンプリングスキームでK空間サンプリングを実行するステップを含み、
    取得された前記K空間サンプルを処理するステップは、分散サンプリングスキームに基づいて取得された前記K空間サンプルに基づいて前記K空間中の穴を再構成するステップを含む、
    請求項1ないし10いずれか一項に記載の磁気共鳴イメージング方法。
  12. MRイメージングシステムをアップグレードするソフトウェアパッケージであって、請求項1ないし11の方法のうち何れか一項に記載の方法により、前記MRイメージングシステムを制御する命令を含む、ソフトウェアパッケージ。
  13. MRイメージングシステムの検査空間に配置された関心被験者のエリアの画像表現を提供するMRイメージングシステムであって、
    前記MRイメージングシステムは、請求項1乃至12いずれか一項に記載の磁気共鳴イメージング方法を実行するように適応されている、MRイメージングシステム。
JP2016541371A 2013-12-20 2014-12-20 磁気共鳴イメージングのための息止め検出 Expired - Fee Related JP6088711B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP13198693 2013-12-20
EP13198693.7 2013-12-20
PCT/EP2014/078935 WO2015092062A1 (en) 2013-12-20 2014-12-20 Breath-hold detection for magnetic resonance imaging

Publications (2)

Publication Number Publication Date
JP2017500128A JP2017500128A (ja) 2017-01-05
JP6088711B2 true JP6088711B2 (ja) 2017-03-01

Family

ID=49920026

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016541371A Expired - Fee Related JP6088711B2 (ja) 2013-12-20 2014-12-20 磁気共鳴イメージングのための息止め検出

Country Status (5)

Country Link
US (1) US10401463B2 (ja)
EP (1) EP3084462B1 (ja)
JP (1) JP6088711B2 (ja)
CN (1) CN106030330B (ja)
WO (1) WO2015092062A1 (ja)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102015204483A1 (de) * 2015-03-12 2016-09-15 Siemens Healthcare Gmbh Magnetresonanz-Vorschau-Abbildung
DE102016203940B4 (de) 2016-03-10 2018-05-30 Siemens Healthcare Gmbh MR-Bildgebung mit optimiertem Bildgebungsarbeitsablauf
CN110226099B (zh) * 2017-01-23 2022-04-29 皇家飞利浦有限公司 在对象运动期间采集四维磁共振数据
DE102017202535A1 (de) * 2017-02-16 2018-08-16 Siemens Healthcare Gmbh Erstellen von MR-Bildern
CN110573896B (zh) * 2017-06-29 2022-05-24 上海联影医疗科技股份有限公司 用于磁共振成像加速的系统和方法
EP3486672A1 (en) * 2017-11-16 2019-05-22 Koninklijke Philips N.V. Magnetic resonance imaging system with rf motion detection
DE102018209359A1 (de) * 2018-06-12 2019-12-12 Siemens Healthcare Gmbh Bestimmung eines Merkmals charakterisierend ein beabsichtigtes Anhalten einer Atmung
US10672151B1 (en) * 2019-01-07 2020-06-02 Uih America, Inc. Systems and methods for magnetic resonance image reconstruction
CN109917315B (zh) 2019-04-30 2021-09-28 上海联影医疗科技股份有限公司 磁共振成像扫描方法、装置、计算机设备和存储介质

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62106755A (ja) 1985-11-02 1987-05-18 株式会社東芝 磁気共鳴イメ−ジング装置
US5417213A (en) 1993-06-07 1995-05-23 Prince; Martin R. Magnetic resonance arteriography with dynamic intravenous contrast agents
US5363844A (en) 1993-08-13 1994-11-15 Mayo Foundation For Medical Education And Research Breath-hold monitor for MR imaging
US5830143A (en) 1997-01-21 1998-11-03 Wisconsin Alumnin Research Foundation Gated time-resolved contrast-enhanced 3D MR angiography
US7254437B2 (en) * 1998-04-17 2007-08-07 Kabushiki Kaisha Toshiba MR imaging providing tissue/blood contrast image
US6144874A (en) * 1998-10-15 2000-11-07 General Electric Company Respiratory gating method for MR imaging
EP1139114A3 (en) * 2000-03-27 2003-08-06 GE Medical Systems Global Technology Company LLC Slice ordering method for breath-hold abdominal MR imaging
JP4807825B2 (ja) * 2005-07-25 2011-11-02 株式会社日立メディコ 磁気共鳴イメージング装置
US8200311B2 (en) 2006-04-21 2012-06-12 Cornell Research Foundation, Inc. Cardiac motion artifact suppression using ECG ordering
WO2010038850A1 (ja) * 2008-10-03 2010-04-08 株式会社 日立メディコ 磁気共鳴イメージング装置及び息止め撮像方法
US8818487B2 (en) * 2008-10-15 2014-08-26 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP5551367B2 (ja) * 2009-02-05 2014-07-16 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 画像診断装置
DE102009017778A1 (de) * 2009-04-20 2010-11-18 Siemens Aktiengesellschaft Verfahren, Magnetresonanzgerät und Computerprogramm für eine Akquisition von Messdaten eines Untersuchungsbereichs eines Untersuchungsobjekts während einer kontinuierlichen Fahrt des Untersuchungsbereichs durch ein Magnetresonanzgerät
US8306299B2 (en) 2011-03-25 2012-11-06 Wisconsin Alumni Research Foundation Method for reconstructing motion-compensated magnetic resonance images from non-Cartesian k-space data
RU2014142029A (ru) * 2012-03-19 2016-05-20 Конинклейке Филипс Н.В. Способ восстановления магнитно-резонансного изображения с обнаружением дыхательного движения во время дискретизации центральной и переферийной областей k-пространства

Also Published As

Publication number Publication date
US20160313429A1 (en) 2016-10-27
CN106030330B (zh) 2020-02-14
CN106030330A (zh) 2016-10-12
WO2015092062A1 (en) 2015-06-25
JP2017500128A (ja) 2017-01-05
EP3084462B1 (en) 2019-05-01
US10401463B2 (en) 2019-09-03
EP3084462A1 (en) 2016-10-26

Similar Documents

Publication Publication Date Title
JP6088711B2 (ja) 磁気共鳴イメージングのための息止め検出
JP6232043B2 (ja) K空間領域の中央および周辺のサンプリング最中での呼吸運動の検出を用いた磁気共鳴画像の再構成
CN107850652B (zh) 具有运动检测的mr成像
US8768034B2 (en) Motion compensated MR imaging system
US9305376B2 (en) Magnetic resonance imaging apparatus and method of acquiring functional image
JP5599893B2 (ja) ナビゲータを使用するmrイメージング
US9664764B2 (en) Magnetic resonance imaging apparatus and susceptibility-weighted imaging method using the same
JP5138043B2 (ja) 磁気共鳴イメージング装置
JP2015506775A (ja) ディクソン法を使って取得されるナビゲーターを使う動き補正のあるmri
JP2012245350A (ja) 拡散強調エコープラナー撮像法において高次渦電流に誘発された歪みを予測補正するためのシステムおよび方法
US11029383B2 (en) Operating a magnetic resonance device
US10132902B2 (en) Intrinsic navigation from velocity-encoding gradients in phase-contrast MRI
KR101653589B1 (ko) 이동가능 검사 대상의 적어도 하나의 자기 공명 이미지 데이터 세트의 조정을 위한 방법 및 상응적으로 설계된 자기 공명 장치
WO2017119994A1 (en) Methods and systems for correcting k-space trajectories
JP3891667B2 (ja) 磁気共鳴イメージング装置
US10094900B2 (en) 4D velocity profile image recording with a magnetic resonance system
JP2010233907A (ja) 磁気共鳴イメージング装置及び感度補正方法
US20100045292A1 (en) Magnetic resonance angiography method and apparatus
JP4391214B2 (ja) 磁気共鳴イメージング装置
JP2005111276A (ja) 低信号領域を伴うsenseのための画質改善
JP5421600B2 (ja) 核磁気共鳴イメージング装置および核磁気共鳴イメージング装置の作動方法
JP3952310B2 (ja) 核磁気共鳴撮影装置
US20220230318A1 (en) Non-respiratory body movement detection in respiratory triggering
EP1685421B1 (en) Readout ordering in collection of radial magnetic resonance imaging data
US11112480B2 (en) Compressed sensing with regularization parameter

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20160803

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20160803

TRDD Decision of grant or rejection written
A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20161227

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20170110

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20170203

R150 Certificate of patent or registration of utility model

Ref document number: 6088711

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

LAPS Cancellation because of no payment of annual fees