JP7023954B2 - プロペラmrイメージング - Google Patents

プロペラmrイメージング Download PDF

Info

Publication number
JP7023954B2
JP7023954B2 JP2019527912A JP2019527912A JP7023954B2 JP 7023954 B2 JP7023954 B2 JP 7023954B2 JP 2019527912 A JP2019527912 A JP 2019527912A JP 2019527912 A JP2019527912 A JP 2019527912A JP 7023954 B2 JP7023954 B2 JP 7023954B2
Authority
JP
Japan
Prior art keywords
space
signal
image
blades
relaxation
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.)
Active
Application number
JP2019527912A
Other languages
English (en)
Other versions
JP2019535435A (ja
JP2019535435A5 (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 JP2019535435A publication Critical patent/JP2019535435A/ja
Publication of JP2019535435A5 publication Critical patent/JP2019535435A5/ja
Application granted granted Critical
Publication of JP7023954B2 publication Critical patent/JP7023954B2/ja
Active 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/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
    • 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/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/4828Resolving the MR signals of different chemical species, e.g. water-fat imaging
    • 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/50NMR imaging systems based on the determination of relaxation times, e.g. T1 measurement by IR sequences; T2 measurement by multiple-echo sequences
    • 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/5602Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by filtering or weighting based on different relaxation times within the sample, e.g. T1 weighting using an inversion pulse
    • 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/5613Generating steady state signals, e.g. low flip angle sequences [FLASH]
    • 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/5615Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/006Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00Indexing scheme for image generation or computer graphics
    • G06T2210/41Medical

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Description

本発明は磁気共鳴(MR)イメージングの分野に関する。それは、MR装置の検査ボリューム内に配置される患者の身体のMRイメージングの方法に関する。本発明は更に、MR装置、及びMR装置上で実行されるコンピュータプログラムに関する。
2次元又は3次元画像を形成するために磁場と核スピンとの間の相互作用を利用する画像形成MR方法は、軟組織の撮像に関して、それらが多くの点で他の撮像方法より優れており、電離放射線を必要とせず、そして通常侵襲的ではないため、今日、特に医療診断の分野において広く使用されている。
一般にMR方法によれば、検査される患者の身体は、強く均一な磁場B内に配置され、その方向は、同時に、測定が関連付けられる座標系の軸(通常z軸)を規定する。磁場Bは、磁場強度に応じて、規定される周波数(いわゆるラーモア周波数又はMR周波数)の電磁交番磁場(RF場)の印加によって励起(スピン共鳴)されることができる個々の核スピンについて異なるエネルギーレベルを生成する。巨視的な観点から、個々の核スピンの分布は、全体的な磁化を生じさせ、全体的な磁化は、適切な周波数の電磁パルス(RFパルス)の印加によって平衡状態から外れるように偏らされ、その一方、このRFパルスの対応する磁場Bは、z軸に垂直に延在し、その結果、磁化は、z軸を中心に歳差運動を行う。歳差運動は、円錐の表面を表現し、その開口角はフリップ角と呼ばれる。フリップ角の大きさは、印加される電磁パルスの強度及び持続時間に依存する。いわゆる90°パルスの場合、磁化はz軸から横方向断面に偏らされる(フリップ角90°)。
RFパルスの終了後、磁化は元の平衡状態に戻るように緩和し、z方向の磁が、は第1の時定数T(スピン格子又は縦緩和時間)で再び増大し、z方向に垂直な方向の磁化は、第2のより短い時定数T(スピン-スピン又は横緩和時間)で緩和する。横方向磁化及びその変動は、磁化の変動がz軸に垂直な方向において測定されるように、MR装置の検査ボリューム内に配置され方向付けられる受信RFコイルによって検出されることができる。横磁化の減衰は、局所的な磁場の不均一性によって引き起こされるRF励起の後に起こるディフェージングを伴い、これは、同じ信号位相を有する秩序のある状態から全ての位相角が一様に分布する状態への遷移を促進する。ディフェージングは、リフォーカスパルス(例えば180°パルス)によって補償されることができる。これは、受信コイルにエコー信号(スピンエコー)を生じさせる。
体内での空間分解能を実現するために、3つの主軸に沿って延在する時間変動する磁場勾配が、一様な磁場Bに重ねられ、それによりスピン共鳴周波数の線形空間依存性を与える。受信コイルにおいて収集される信号は、体内のそれぞれ異なるロケーションに関連付けられる複数の異なる周波数成分を含む。受信コイルを介して得られる信号データは、空間周波数ドメインに対応し、k空間データと呼ばれる。k空間データは通常、それぞれ異なる位相エンコードから取得される複数のラインを含む。各ラインは、多数のサンプルを収集することによってデジタル化される。1組のk空間データが、フーリエ変換によってMR画像に変換される。
いわゆるPROPELLERイメージングは、患者の動きに対するロバストネス及び固有の動き補償能力のような臨床MRイメージングに対して有望な利点を提供する。PROPELLER概念(Periodically Rotated Overlapping ParallEL Lines、 James G. Pipe: "Motion Correction With PROPELLER MRI: Application to Head Motion and Free-Breathing Cardiac Imaging", Magnetic Resonance in Medicine, vol. 42, 1999, pages 963-969参照)において、MR信号データは、k空間においてN本のストリップで取得され、各ストリップは、Lの平行なk空間ラインを含み、かかるラインは、デカルトベースのk空間サンプリングスキームにおけるL本の最低周波数位相エンコードラインに対応する。k空間ブレードとも呼ばれる各ストリップは、MRデータセット全体がk空間内の円をスパンするように、例えば180°/Nの角度ずつk空間内で回転される。直径Mを有する全k空間データ行列が望ましい場合、L及びNは、L×N=M×π/2であるように選択されることができる。PROPELLERの本質的な特徴は、直径Lを有するk空間内の中央の円形部分が、各k空間ブレード毎に取得されることである。この中央部分は、各k空間ブレード毎に低解像度のMR画像を再構成するために使用されることができる。低解像度MR画像が互いに比較されることにより、患者の動きに起因する面内の変位及び位相誤差を除去することができる。PROPELLER技法は、MR信号取得中の検査対象の患者の動きに対してロバストなMR画像取得技法を得るために、k空間の中央部分でのオーバーサンプリングを利用する。更に、SENSEのようなパラレルイメージング技術がMRデータ収集のために使用される場合、PROPELLERは、k空間ブレードの(加重)平均により、例えばB不均一性又は不正確なコイル感度マップに起因する他のイメージングアーチファクトを「平均化」する。
PROPELLERイメージングは、単一のRF励起の後に完全なk空間ブレードが取得されるマルチエコー取得と組み合わせるのに非常に適していることが知られている。例えば、TSE PROPELLERイメージングにおいて、PROPELLER取得スキームは、ターボスピンエコー(TSE)イメージングシーケンスと組み合わされる。TSEシーケンスの単一ショットは、磁気共鳴を励起するための初期RFパルスと、それに続く、高速に印加される複数の(典型的には180°)リフォーカスRFパルスと、を有し、それにより複数の異なる位相エンコードされたスピンエコー信号のトレインを生成する。エコー信号が取得され、各エコー信号は、k空間ライン、すなわちk空間の1次元サンプルを表し、個々のk空間ブレード内のk空間ラインの位置は、シーケンスの適用される周波数エンコード及び位相エンコードにより切り替えられる磁場勾配によって決定される。いわゆるターボファクタ(TF)は、各励起の後に取得されるエコーの数である。PROPELLER TSEイメージングでは、マルチエコーシーケンスの複数のショットが、取得されるk空間プロファイルからMR画像を再構成することができるように、k空間を完全にサンプリングするために適用され、空間ブレードは、それぞれ単一のRF励起の後に取得される。従って、TSEシーケンスのショット数は、通常、取得されるk空間ブレードの数に対応し、ターボファクタは、通常、k空間ブレードあたりのk空間ラインの数に対応する。
このような利点にもかかわらず、PROPELLER TSEは、k空間ブレードの大きさ及び回転、並びに種々異なる量の緩和減衰を経験するエコー信号のPROPELLER特有のk空間分散取得によってコントラストブラーをもたらすことが知られている。特にプロトン密度(PD)及びTコントラストは、一般により短いTSEエコー時間及びエコートレイン長により、増大されたコントラストブラーをもたらす。従って、PD及びTコントラスト強調MRイメージングは、典型的には、PROPELLER TSEイメージングとの組み合わせにおいて今日使用されていない。
D. Huo他によるアブストラクト"IDEAL with Turbo-PROP", ISMRM-2008 (p.650)は、ターボPROPシーケンスからのマルチエコー画像の使用に関する。これらの画像は、位相差に基づいて水信号と脂肪信号を分離するために使用されることができる。
前述から、改善されたMRイメージング技術の必要があることが容易に理解される。従って、本発明の目的は、マルチエコー取得と組み合わせられるPROPELLERイメージングにおけるコントラストブラーを低減する方法を提供することである。
本発明によれば、MR装置の検査ボリューム内に配置される患者の身体のMRイメージングの方法が開示される。この方法は、
-身体の少なくとも一部を、複数のRFパルス及び切り替え磁場勾配を含むMRイメージングシーケンスに曝すことによってMR信号を生成するステップと、
-PROPELLERスキームに従って時間的に連続して複数のk空間ブレードとしてMR信号を取得するステップであって、各k空間ブレードは、複数の実質的に平行なk空間ラインを含み、k空間ブレードは、k空間の中心を中心に回転され、それにより、MR信号の収集されるデータセット全体が、k空間内の円の少なくとも一部をスパンし、k空間の共通の中央円形領域が、全てのk空間ブレードによってカバーされ、MR信号の強調が、個々異なるk空間ブレードの間で変化する、ステップと、
-MR信号の緩和強調(relaxation weighting)を推定するステップと、
-推定された緩和強調に従って、取得されたMR信号を補償するステップと、
-補償されたMR信号からMR画像を再構成するステップと、
を有する。
本発明は、それぞれ異なるk空間ブレードのk空間ラインがそれぞれ異なる緩和時間コントラストを有する、PROPELLERスキームによるk空間ブレードの取得を提案する。例えば、T強調(T1-weighting)及び/又はT強調(T2-weighting)及び/又は個別のk空間ラインのスペクトルに関して誘起される位相が、個々異なるk空間ブレードの間で変化しうる。このようにして、更に、個々異なるコントラストを有するMR信号の集合が、1回のPROPELLER取得中にk空間の中心から取得される。本発明の洞察は、これが、信号取得中にMR信号の緩和減衰及び/又はスペクトルに関して誘起される位相を推定することを可能にし、それに応じてMR信号が補償され、よって、「ブラー除去された」MR画像が、補償されたMR信号からの画像再構成により、取得されることができることである。
本発明の方法は、画像のコントラストに関する重要な情報を含むk空間の中心のPROPELLERスキャニングによる固有の実質的なオーバーサンプリングを利用する。従って、本発明は、取得中に適用される各緩和毎に、個々の緩和強調される及び/又はスペクトル強調されるMR画像を再構成することを可能にする。MR信号の緩和強調を推定するために、対応するTマップ及び/又はTマップ及び/又はBマップ及び/又は水マップ及び/又は脂肪マップが、取得から導出されることができ、取得されたMR信号は、導出されたそれぞれのマップに従って変動に関して補償される。最後に、このようにして補償されたMR信号から、「ブラー除去された」MR画像を再構成することができる。
特定の緩和時間に属するMR信号データは、一般に、少なくともk空間の周辺部においてアンダーサンプリングされる。従って、アンダーサンプリングされたMR信号データから、それぞれのTマップ及び/又はTマップ及び/又はBマップ及び/又は水マップ及び/又は脂肪マップを導出するために、圧縮センシング(CS)が有利に使用されることができる。特定の緩和時間のk空間ラインは、k空間において不規則に分布しうる。CSの理論は、大幅な信号データ削減の大きな可能性を有することが知られている。CS理論では、変換ドメインにおいてスパース表現を有する信号データセットが、適切な正則化アルゴリズムを適用することによりアンダーサンプリングされた測定から回復されることができる。信号のサンプリング及び再構成の数学的フレームワークとして、CSは、k空間のサンプリング密度がナイキスト基準をはるかに下回る場合にでも、信号データセットが正確に又は少なくとも高い画像品質を伴って再構成されることができる条件を規定し、そのような再構成の方法も提供する。
更に、純粋なPD及びTコントラスト強調されたMR画像は、本発明による単一PROPELLER取得から生成されることができる。更に、「混合された」(PD、T、T)コントラストをもつMR画像が、1回のPROPELLER取得のMR信号データから合成されることができる。更に、定量的緩和時間マップ及び/又はBマップ、水マップ及び/又は脂肪マップが、取得データから導出されることができる。
好適な実施形態において、MR信号の緩和時間は、撮像シーケンスのエコー時間の変動によって、及び/又は(反転回復)事前パルス後の遅延(たとえば反転)時間の変動によって、変わる。従って、k空間ブレードは、個々異なるT強調を得るために個々異なるエコー時間を使用することにより、及び/又は個々異なるT強調を得るために個々異なる反転遅延を使用することにより、それぞれに異なる。
更に好適な実施形態において、撮像シーケンスは、エコー信号のトレインが各ショットによって生成されるマルチエコー撮像シーケンスでありえ、各k空間ブレードの各k空間ラインは、エコー信号を表す。撮像シーケンスは、例えば、ターボスピンエコー(TSE)シーケンス又はターボフィールドエコー(TFE)シーケンス又はエコープラナーイメージング(EPI)シーケンスであり得る。これにより、撮像シーケンスの単一ショットでフルk空間ブレードを高速に取得することが可能である。この実施形態では、緩和時間は、本発明に従って、ブレード毎にk空間ラインの取得順序を変更することによって変えられることができる。例えば、k空間ブレード内の所与の位置におけるk空間ラインのT強調は、生成されたエコー信号のトレイン内の対応するエコー信号の位置に依存する。従って、緩和時間強調は、k空間ラインの取得順序をブレード毎に変えることによって容易に変えられることができる。代替として、緩和時間は、個々のエコー信号のエコーシフトを変えることによって変えられることができる。
本発明により、複数のPROPELLERブレードの取得中、k空間の中心において個々異なるコントラストを有するMR信号のk空間分布を一層最適化するために、ブレードの角度に関する順序付けは、黄金角(golden angle)スキームに従って選択されることができる。黄金角スキームでは、k空間ラインの角度は、ブレード毎に毎回ΔΦ=111.25°ずつインクリメントされ、これは、180°に黄金比を乗じた値に一致する。従って、その後にサンプリングされるスポークは、以前にサンプリングされたスポークの組の中での最大のギャップを埋めつつ、補足情報を常に追加する。その結果、取得されたスポークの順次の組は、ほぼ一様にk空間をカバーする。これにより、例えば最初の3つの黄金角ブレードは、緩和時間の第1の設定を使用して取得され、それに続いて、3つの他の黄金角ブレードが、別の緩和時間を使用して取得される。更に、各ブレード内のエコー順序が、コントラスト分布を一層強めるために変えられることができる。
本発明の他の好適な実施形態において、MRイメージングシーケンスは、RFパルスのフリップ角スイープを含み、各々の取得されたMR信号に属するフリップ角は、緩和強調の推定において考慮される。k空間ラインのT強調は、TSE取得において通常使用されるフリップ角スイープに影響されうる。フリップ角スイープに関する情報は、緩和時間強調の推定において、例えば対応するモデルを通じて組み込まれる必要がある。更に、フリップ角スイープを調整することにより、緩和強調の推定の精度を向上させることもできる。例えば、k空間ブレードの総数が少ない場合、TSEショットの最中の緩和時間強調の数は、フリップ角スイープを使用して平坦にされることができる。
従来のPROPELLERイメージングのように、本発明の方法は、更に、k空間ブレード内の動きにより引き起こされる変位及び位相誤差を推定し補正するステップを含むことができる。例えば、k空間ブレードの中心k空間データから再構成された低解像度MR画像は、患者の動きによって引き起こされる面内変位及び位相誤差を除去するために互いに比較される。これらのファクタは、MR信号データを更に処理する前に、本発明に従って各k空間ブレードにおいて補正されるべきである。これは、本発明の方法を、MR信号取得中の検査された患者の動きに対してロバストにする。緩和時間強調の変動は、各低解像度MR画像のコントラストをも変えるので、患者の動きを正確に決定し補正するためには、追加の方策が必要となりうる。相互相関及び正規化相互情報等の類似性方法が、取得されたマルチコントラスト低解像度画像をアラインするために使用されることができる。
本発明の好適な実施形態によれば、撮像シーケンスは、エコープラナーイメージング(EPI)シーケンスである。スペクトル強調の変動は、個々のk空間ブレード間のEPIエコートレインのシフトの対応する変動によって達成されることができる。個々異なるエコー時間における個々の画像が、そこから再構成されることができ、水マップ、脂肪マップ及びBマップは、既知のディクソンスキームに従って再構成されることができる。これらは、その後、Bマップを考慮に入れることによって歪みのないMR画像を提供するために使用されることができる。MR画像は、所望の水/脂肪の表現と合成されることができる。
これまで述べた本発明の方法は、検査ボリューム内に一様な定常磁場Bを生成する少なくとも1つの主磁石コイルと、検査ボリューム内のそれぞれ異なる空間方向に切り替え磁場勾配を生成する複数の勾配コイルと、検査ボリューム内にRFパルスを生成する及び/又は検査ボリューム内に位置する患者の身体からMR信号を受信する少なくとも1つの身体RFコイルと、時間的に連続するRFパルス及び切り替え磁場勾配を制御する制御ユニットと、受信されたMR信号からMR画像を再構成する再構成ユニットと、を有するMR装置によって実施されることができる。本発明の方法は、MR装置の再構成ユニット及び/又は制御ユニットの対応するプログラミングによって実施されることができる。
本発明の方法は、今日臨床的に使用されている多くのMR装置上で有利に実行されることができる。この目的のために、MR装置が本発明の上記方法ステップを実行するように該MR装置が制御されるコンピュータプログラムを利用することのみが必要である。コンピュータプログラムは、MR装置の制御ユニットにインストールするためにダウンロードされるように、データ担体上又はデータネットワーク内のいずれに存在するものでもよい。
本発明の方法を実行するMR装置を示す図。 本発明のプロペラ取得スキームを概略的に示す図。 本発明のプロシージャを示すブロック図。
添付の図面は、本発明の好適な実施形態を開示している。しかしながら、図面は例示の目的でのみ設計されており、本発明の範囲を規定するものとして設計されていないことを理解されたい。
図1を参照すると、MR装置1がブロック図として示されている。装置は、超電導又は抵抗主磁石コイル2を有し、それにより、実質的に均一で時間的に一定の主磁場Bが、検査ボリュームを通ってz軸に沿って生成される。装置は更に、1組の(1次、2次、及び適用可能な場合は3次の)シミングコイル2'を有し、組2'の個々のシミングコイルを通る電流は、検査ボリューム内の主磁場Bの偏りを最小にする目的で制御可能である。
磁気共鳴生成及び操作システムは、核磁気スピンを反転又は励起し、磁気共鳴を引き起こし、磁気共鳴をリフォーカスし、磁気共鳴を操作し、磁気共鳴を空間的に及びその他の態様でエンコードし、スピンを飽和させる等により、MRイメージングを実施するために、一連のRFパルス及び切り替え磁場勾配を印加する。
より具体的には、勾配増幅器3は、検査ボリュームのx、y及びz軸に沿って全身勾配コイル4、5及び6のうちの選択されたものに電流パルス又は波形を印加する。デジタルRF周波数送信器7は、RFパルスを検査ボリューム内に送信するために、送信/受信スイッチ8を介して身体RFコイル9にRFパルス又はパルスパケットを送信する。一般的なMRイメージングシーケンスは、短い持続時間のRFパルスセグメントのパケットで構成され、かかるパケットは印加される磁場勾配と共に、核磁気共鳴の選択された操作を達成する。RFパルスは、共鳴を飽和させ、共鳴を励起し、磁化を反転し、共鳴をリフォーカスし、共鳴を操作し、又は検査ボリューム内に位置付けられた身体10の一部を選択するために、使用される。MR信号は、身体RFコイル9によっても収集される。
身体10の限られた領域のMR画像を生成するために、又はパラレルイメージングによるスキャン加速のために、1組のローカルアレイRFコイル11、12、13が、イメージングのために選択される領域に隣接して配置される。アレイコイル11、12、13は、身体コイルRF送信によって引き起こされるMR信号を受信するために使用されることができる。
結果的に生じるMR信号は、身体RFコイル9及び/又はアレイRFコイル11、12、13によって収集され、好適には前置増幅器(図示せず)を有する受信器14によって復調される。受信器14は、送信/受信スイッチ8を介してRFコイル9、11、12、13に接続される。
ホストコンピュータ15は、シミングコイル2'並びに勾配パルス増幅器3及び送信器7を制御して、エコープラナーイメージング(EPI)、エコーボリュームイメージング、グラジエント及びスピンエコーイメージング、高速スピンエコーイメージング、その他の複数のMRイメージングシーケンスのうちのいずれかを生成する。選択されたシーケンスに関して、受信器14は、各RF励起パルスに続き、単一又は複数のMRデータラインを高速に連続して受信する。データ取得システム16は、受信信号のアナログ-デジタル変換を実施し、各MRデータラインを更なる処理に適したデジタル形式に変換する。今日のMR装置において、データ取得システム16は、生の画像データの取得に特化した別のコンピュータである。
最終的に、デジタル生画像データは、フーリエ変換又はSENSE若しくはGRAPPA等の他の適切な再構成アルゴリズムを適用する再構成プロセッサ17によって画像表現に再構成される。MR画像は、患者を通る平面スライス、平行な平面スライスのアレイ、3次元ボリューム、その他を表現することができる。画像は、画像メモリに記憶され、画像メモリは、画像表現のスライス、投影、又は他の部分を、例えばビデオモニタ18を介して視覚化するために適した形式に変換するためにアクセスされることができ、ビデオモニタは、結果的に得られるMR画像の人間可読表示を与える。
ホストコンピュータ15は、本明細書の上記及び下記の本発明の方法を実行するようにプログラムされる。
図2は、本発明によるPROPELLER MRイメージングのk空間サンプリングを示す。図2の左側の図に示されるように、6つのk空間ブレード21-26が取得される。各ブレード21-26は、k空間のそれぞれ異なる部分をカバーし、各ブレード21-26毎にk空間の中央円形領域が取得される。取得されるMRデータセット全体がk空間の円をスパンするように、ブレード21-26がk空間の中心を中心に回転される。図2の右側の図には、TSEシーケンスの単一ショットを使用して取得される単一のk空間ブレード21が示されている。ブレードの向きに対する位相エンコード方向及び読み出し方向の配向は、k空間ブレード21-26の全ての回転角度について維持される。ブレード21-26における動きにより引き起こされる変位及び位相誤差は、既知の類似性測定方法を使用して従来のPROPELLERイメージングのように検出され補正される。ブレード21内の個々異なるラインは、k空間ライン27、28、29の取得順序を示す。一例として、実線27は、単一のRF励起の後に取得される一連のエコー信号のうちの第1のエコーに対応する。破線28は、第2のエコーに対応し、点線29は、第3のエコーに対応する。図2の左の像に見られるように、MR信号の緩和時間強調(T強調)が個々の異なるブレード21-26の間で変化するように、エコー信号の取得順序が、ブレード毎に変えられる。
本発明は、図2に示されるように、ブレード毎に異なるTSEコントラスト又はTSE順序を提案する。TSEコントラストの変化は、それが非対称TSE技法において一般的に使用されているように、例えばTSEエコーシフトによって達成されることができる。ブレード毎に異なるTSE順序は、単一のPROPELLER取得中のk空間の中心におけるコントラスト変動を可能にする。これは、それぞれ異なる(実線、破線、点線)線27、28、29によって図2に示されている。個々の緩和時間強調されたMR画像は、例えば、オーバーサンプリングされたk空間中心における反復的SENSE、CS及び/又はコントラスト強調の使用によって再構成されることができる(Song HK et. al, Magnetic Resonance in Medicine, vol. 44, no. 6, p. 825-832, 2000)。個々のコントラストの再構成は、コントラスト減衰(T)を推定し、コントラストマップ(Tマップ)を導出するために本発明に従って使用される。Tマップは、例えば、個々異なるT強調MR画像において信号減衰をピクセル毎に決定することによって取得されることができる。次いで、TSEプロペラデータは、推定される緩和崩壊を用いてブラー除去される。純粋なプロトン密度強調及びT強調MR画像は、1回のPROPELLER取得から生成される。更に、ユーザの好みに基づく混合されたコントラスト(T、T、PD)が合成されることができ、定量的緩和時間マップが更に導出されることができる。
本発明の概念は、異なるシーケンスタイプ及びコントラストの変動に拡張されることができる。別の例として、本発明は、Tマップ及び種々異なるTコントラストを生成する可能性を提供するように、単一反転パルス後、コントラスト変動と組み合わせて、PROPELLER平衡TFE取得に適用されることができる。
図3は、本発明の方法を概略的に示す。ステップ31において、それぞれ異なる緩和時間をもつ複数のk空間ブレードを取得した後、ステップ32において、動きが、適切な類似性方法を使用して推定され補正される。次に、相応して異なるコントラストを有する個々のMR画像33a-eが再構成され、コントラスト強調が、ステップ34において、取得されたコントラストモデルにフィットされて、取得されたMR信号の緩和強調が推定される。次いで、ステップ35において、取得されたMR信号が、コントラストモデルに従って補正される。純粋なコントラストMR画像(例えば、T強調MR画像、T強調MR画像、プロトン密度強調MR画像)が、ステップ39において、補償されたMR信号データから再構成されることができる。フィッティングプロシージャの結果として、コントラストマップ(例えば、Tマップ、Tマップ、脂肪マップ、水マップ)が、ステップ36において付加的に取得される。所望のコントラストに関するユーザの好みは、ステップ37において入力されることができる。これに基づき、ステップ38において、相応して合成され組み合わせられるコントラストを有するMR画像が、MR信号から再構成されることができる。

Claims (15)

  1. MR装置の検査ボリューム内に配される患者の身体のMRイメージング方法であって、
    身体の少なくとも一部を、複数のRFパルス及び切り替え磁場勾配を含むMRイメージングシーケンスに曝すことによってMR信号を生成するステップと、
    PROPELLERスキームに従って時間的に連続して複数のk空間ブレードとしてMR信号を取得するステップであって、各k空間ブレードは、複数の実質的に平行なk空間ラインを含み、前記k空間ブレードは、k空間の中心を中心に回転され、それにより、MR信号の収集されるデータセット全体が、k空間の円の少なくとも一部をスパンし、k空間の共通の中央円形領域が、全てのk空間ブレードによってカバーされ、前記MR信号の緩和強調が、個々異なるk空間ブレードの間で変化し、特定の緩和強調のk空間ラインが、k空間に不規則に分布する、ステップと、
    前記MR信号の緩和強調を推定するステップと、
    前記推定された緩和強調に従って、前記取得されたMR信号を補償するステップと、
    前記補償されたMR信号から、MR画像を再構成するステップと、
    を実行する方法。
  2. 個々のk空間ラインのT1強調及び/又はT2強調が、個々異なる空間ブレードの間で変化する、請求項1に記載の方法。
  3. T1マップ及び/又はT2マップ及び/又はB0マップ及び/又は水マップ及び/又は脂肪マップが、前記取得されたMR信号から導出され、前記取得されたMR信号が、前記導出されたT1マップ、T2マップ、B0マップ、水マップ及び/又は脂肪マップに従って変動に関して補償される、請求項1又は2に記載の方法。
  4. 前記MR信号の前記緩和強調は、前記MRイメージングシーケンスのエコー時間を変化させることによって、及び/又は準備段階の事前パルス後の時間遅延を変化させることによって、変えられる、請求項1乃至3のいずれか1項に記載の方法。
  5. 前記MRイメージングシーケンスは、マルチエコーイメージングシーケンス、特にターボスピンエコー(TSE)シーケンス又はターボフィールドエコー(TFE)又はエコープラナーイメージング(EPI)シーケンス、及び各ショットによって生成されるエコー信号のトレインであり、各k空間ブレードの各k空間ラインは、エコー信号を表す、請求項1乃至4のいずれか1項に記載の方法。
  6. 前記MR信号の前記緩和強調は、個々のk空間ブレードの間で前記k空間ラインの取得順序を変更することによって変えられる、請求項5に記載の方法。
  7. 前記MR信号の前記緩和強調は、個々のk空間ブレードの間で前記エコー信号のエコーシフトを変更することによって変えられる、請求項5に記載の方法。
  8. 緩和強調された及び/又はスペクトル強調されたMR画像が、各緩和強調毎に前記取得されたMR信号から再構成される、請求項1乃至7のいずれか1項に記載の方法。
  9. 緩和時間マップが、前記緩和強調されて再構成されるMR画像から導出される、請求項8に記載の方法。
  10. 前記緩和強調されたMR画像は、圧縮センシング又はパラレルイメージング再構成アルゴリズムを使用して再構成される、請求項8又は9に記載の方法。
  11. T1強調MR画像、T2強調MR画像、プロトン密度強調MR画像、水MR画像、及び/又は脂肪MR画像が、前記補償されたMR信号から再構成される、請求項1乃至10のいずれか1項に記載の方法。
  12. T1強調、T2強調、プロトン密度強調、水及び/又は脂肪コントラストを組み合わせたMR画像が、前記補償されたMR信号から合成される、請求項1乃至11のいずれか1項に記載の方法。
  13. 前記MRイメージングシーケンスが、前記RFパルスのフリップ角スイープを含み、取得された各MR信号に関するフリップ角が、前記緩和強調の推定において考慮される、請求項1乃至12のいずれか1項に記載の方法。
  14. 請求項1乃至13のいずれか1項に記載の方法を実行するMR装置であって、前記MR装置は、検査ボリューム内に均一な定常磁場を生成する少なくとも1つの主磁石コイルと、前記検査ボリューム内の個々異なる空間方向に切り替え磁場勾配を生成する複数の勾配コイルと、前記検査ボリューム内にRFパルスを生成し及び/又は前記検査ボリューム内に位置する対象物からMR信号を受信する少なくとも1つのRFコイルと、時間的に連続する前記RFパルス及び前記切り替え磁場勾配を制御する制御ユニットと、前記受信されるMR信号からMR画像を再構成する再構成ユニットと、を有し、前記MR装置は、
    前記身体の少なくとも一部を、複数のRFパルス及び切り替え磁場勾配を含むMRイメージングシーケンスに曝すことによりMR信号を生成するステップと、
    ROPELLERスキームに従って時間的に連続して複数のk空間ブレードとしてMR信号を取得するステップであって、各k空間ブレードは、複数の実質的に平行なk空間ラインを有し、前記k空間ブレードは、k空間の中心を中心に回転され、それにより、前記MR信号の収集されたデータセット全体が、k空間の円の少なくとも一部をスパンし、k空間の共通の中央円形領域が、全てのk空間ブレードによってカバーされ、前記MR信号の緩和強調が、個々異なるk空間ブレードの間で変化し、特定の緩和強調のk空間ラインが、k空間に不規則に分布する、ステップと、
    前記MR信号の緩和強調を推定するステップと、
    前記推定された緩和強調に従って、前記取得されたMR信号を補償するステップと、
    前記補償されたMR信号からMR画像を再構成するステップと、
    を実行する、MR装置。
  15. MR装置上で実行されるコンピュータプログラムであって、
    複数のRFパルス及び切り替え磁場勾配を含むMRイメージングシーケンスを実行することによってMR信号を生成する命令と、
    PROPELLERスキームに従って時間的に連続して複数のk空間ブレードとしてMR信号を取得する命令であって、各k空間ブレードは、複数の実質的に平行なk空間ラインを含み、前記k空間ブレードは、k空間の中心を中心に回転され、それにより、前記MR信号の収集されたデータセット全体が、k空間の円の少なくとも一部をスパンし、k空間の共通の中央円形領域が、全てのk空間ブレードによってカバーされ、前記MR信号の緩和強調が、個々異なるk空間ブレードの間で変化し、特定の緩和強調のk空間ラインが、k空間に不規則に分布する、命令と、
    前記MR信号の緩和強調を推定する命令と、
    前記推定された緩和強調に従って、前記取得されたMR信号を補償する命令と、
    前記補償されたMR信号からMR画像を再構成する命令と、
    を有するコンピュータプログラム。
JP2019527912A 2016-11-24 2017-11-20 プロペラmrイメージング Active JP7023954B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP16200465.9 2016-11-24
EP16200465 2016-11-24
PCT/EP2017/079745 WO2018095851A1 (en) 2016-11-24 2017-11-20 Propeller mr imaging

Publications (3)

Publication Number Publication Date
JP2019535435A JP2019535435A (ja) 2019-12-12
JP2019535435A5 JP2019535435A5 (ja) 2020-12-24
JP7023954B2 true JP7023954B2 (ja) 2022-02-22

Family

ID=57394462

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019527912A Active JP7023954B2 (ja) 2016-11-24 2017-11-20 プロペラmrイメージング

Country Status (5)

Country Link
US (1) US10996305B2 (ja)
EP (1) EP3545324A1 (ja)
JP (1) JP7023954B2 (ja)
CN (1) CN109983358B (ja)
WO (1) WO2018095851A1 (ja)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3686621A1 (en) * 2019-01-23 2020-07-29 Siemens Healthcare GmbH Imaging method with variable density in k-space for multi-spin echo magnetic resonance imaging with different delay times after a preparation pulse and/or different echo times
EP3770625B1 (en) * 2019-07-24 2024-06-26 Siemens Healthineers AG Mri using compressed sensing with improved regularization parameter
EP4086649A1 (de) * 2021-05-06 2022-11-09 Siemens Healthcare GmbH Adaptive rekonstruktion von mr-daten

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006304955A (ja) 2005-04-27 2006-11-09 Ge Medical Systems Global Technology Co Llc 磁気共鳴イメージング装置
JP2013078562A (ja) 2011-09-22 2013-05-02 Toshiba Corp 磁気共鳴イメージング装置
JP2015531251A (ja) 2012-09-04 2015-11-02 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. ディクソン水脂肪分離を用いるプロペラ

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4684891A (en) * 1985-07-31 1987-08-04 The Regents Of The University Of California Rapid magnetic resonance imaging using multiple phase encoded spin echoes in each of plural measurement cycles
US6882148B2 (en) * 2003-07-09 2005-04-19 Catholic Healthcare West Split-blade data collection for propeller MRI
US20090227860A1 (en) * 2004-08-13 2009-09-10 Koninklijke Philips Electronics, N.V. Mr method for the quantitative determination of local relaxation time values
CN1327809C (zh) * 2005-03-28 2007-07-25 南方医科大学 在t1加权磁共振成像中用propeller采样方式消除运动伪影的方法
JP5781760B2 (ja) * 2007-04-27 2015-09-24 コーニンクレッカ フィリップス エヌ ヴェ プロペラmri用磁気共鳴装置及び方法
WO2008134891A1 (en) * 2007-05-03 2008-11-13 National Research Counsil Of Canada Rf based spatially selective excitation in mri
WO2008135885A1 (en) * 2007-05-03 2008-11-13 Koninklijke Philips Electronics N.V. Propeller mri with phase correction
SE0950363A1 (sv) * 2009-05-22 2010-06-29 Cr Dev Ab Metod och system för magnetisk resonanstomografi, samt användning därav.
CN102959425A (zh) * 2010-04-29 2013-03-06 伊利诺伊大学受托管理委员会 使用转向螺旋桨技术的磁共振成像
CN101995561B (zh) * 2010-10-22 2012-11-07 南方医科大学 基于图像域叠加的propeller磁共振数据重建方法
WO2012137606A1 (ja) 2011-04-06 2012-10-11 株式会社 日立メディコ 磁気共鳴イメージング装置及びラディアルサンプリング方法
US9322894B2 (en) * 2012-08-07 2016-04-26 General Electric Company Multiple excitation blade acquisition for motion correction in magnetic resonance imaging
CN103777162B (zh) * 2012-10-19 2018-06-05 上海联影医疗科技有限公司 磁共振成像k空间运动伪影矫正并行采集重建方法
DE102013205208B4 (de) * 2013-03-25 2015-02-12 Siemens Aktiengesellschaft Multiecho-Magnetresonanz-Messsequenz mit erhöhter Auflösung
JP2016539722A (ja) * 2013-12-12 2016-12-22 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. マルチエコーによる区分化されたk空間の取得を用いたMR撮像
US9983283B2 (en) * 2015-03-16 2018-05-29 Toshiba Medical Systems Corporation Accelerated MRI using radial strips and undersampling of k-space

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006304955A (ja) 2005-04-27 2006-11-09 Ge Medical Systems Global Technology Co Llc 磁気共鳴イメージング装置
JP2013078562A (ja) 2011-09-22 2013-05-02 Toshiba Corp 磁気共鳴イメージング装置
JP2015531251A (ja) 2012-09-04 2015-11-02 コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. ディクソン水脂肪分離を用いるプロペラ

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
E. Aboussouan et al.,Blade by Blade Compressed Sensing for PROPELLER,Proc. Intl. Soc. Mag. Reson. Med,2009年,4591

Also Published As

Publication number Publication date
EP3545324A1 (en) 2019-10-02
US20200081086A1 (en) 2020-03-12
US10996305B2 (en) 2021-05-04
JP2019535435A (ja) 2019-12-12
CN109983358A (zh) 2019-07-05
CN109983358B (zh) 2021-10-26
WO2018095851A1 (en) 2018-05-31

Similar Documents

Publication Publication Date Title
EP3322997B1 (en) Mr imaging with motion detection
US10401456B2 (en) Parallel MR imaging with Nyquist ghost correction for EPI
US10466322B2 (en) MR imaging using multi-echo k-space acquisition
EP3635425B1 (en) Mr imaging using a stack-of-stars acquisition with variable contrast
US10203394B2 (en) Metal resistant MR imaging
EP2798364B1 (en) Mr imaging with suppression of flow artefacts
US20140070805A1 (en) Mr imaging with b1 mapping
EP3044604B1 (en) Metal resistant mr imaging
WO2016180983A1 (en) Multi-echo mri with correction of concomitant gradient-induced phase errors
CN107810425B (zh) 消除非t2加权信号贡献的t2加权mr成像
JP7023954B2 (ja) プロペラmrイメージング
US20140121492A1 (en) Mri with separation of different chemical species using a spectral model
EP3447517A1 (en) Dixon-type water/fat separation mr imaging
US10859652B2 (en) MR imaging with dixon-type water/fat separation
EP3185029A1 (en) Mr imaging using propeller acquisition with t2 decay correction
EP3736593A1 (en) Dual echo steady state mr imaging using bipolar diffusion gradients
WO2020173688A1 (en) Epi mr imaging with distortion correction
US20220390538A1 (en) Dixon-type water/fat separation mr imaging

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20201111

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20201111

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20210826

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210914

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20211203

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220209

R150 Certificate of patent or registration of utility model

Ref document number: 7023954

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150