JP6730995B2 - 環境内に動きがある物体のｍｒ画像を生成する方法及びシステム - Google Patents
環境内に動きがある物体のｍｒ画像を生成する方法及びシステム Download PDFInfo
- Publication number
- JP6730995B2 JP6730995B2 JP2017542280A JP2017542280A JP6730995B2 JP 6730995 B2 JP6730995 B2 JP 6730995B2 JP 2017542280 A JP2017542280 A JP 2017542280A JP 2017542280 A JP2017542280 A JP 2017542280A JP 6730995 B2 JP6730995 B2 JP 6730995B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- data set
- region
- dynamic
- interest
- 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
Links
- 238000000034 method Methods 0.000 title claims description 22
- 230000003068 static effect Effects 0.000 claims description 36
- 230000000747 cardiac effect Effects 0.000 claims description 21
- 238000002595 magnetic resonance imaging Methods 0.000 claims description 20
- 230000035945 sensitivity Effects 0.000 claims description 14
- 230000000873 masking effect Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 4
- 238000004590 computer program Methods 0.000 claims description 2
- 238000013459 approach Methods 0.000 description 6
- 238000013184 cardiac magnetic resonance imaging Methods 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 5
- 238000005070 sampling Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 210000000038 chest Anatomy 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000013480 data collection Methods 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 210000000115 thoracic cavity Anatomy 0.000 description 1
- 238000012549 training Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
- G01R33/5611—Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
- G01R33/56308—Characterization of motion or flow; Dynamic imaging
- G01R33/56325—Cine imaging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2576/00—Medical imaging apparatus involving image processing or analysis
- A61B2576/02—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
- A61B2576/023—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the heart
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
- G01R33/5619—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences by temporal sharing of data, e.g. keyhole, block regional interpolation scheme for k-Space [BRISK]
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Signal Processing (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Vascular Medicine (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Pathology (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Description
本発明は、一定の期間内に複数の動作フェーズ（運動フェーズ）を含む動きを実行する物体のＭＲ画像（ＭＲ：磁気共鳴）を生成する分野に関する。本発明は更に、一定の期間内に複数の動作フェーズからなる動きを実行する物体のＭＲ画像を生成する対応するＭＲＩシステムに関する。
心臓ＭＲＩ（ＣＭＲＩ）の適用可能性を低下させる２つの制限要因がある。第１の制限要因は、ＣＭＲＩ検査は実行すべきために必要なスキャナ時間が長く、関連情報が取得されるまでに多くのダミーの計画スキャンを必要とすることである。第２の制限要因は、検査助手が正しい心臓の向きで心臓検査を行うことができるためには、周囲に広めるには困難な高度な訓練が必要なことである。胸部全体をカバーする３Ｄ等方性ボリューム収集がこれらの２つの制限を回避するのに役立ち、心臓検査のワークフローを改善する。
Ｌ．Ｈ．Ｈａｍｉｌｔｏｎらによる科学論文「Ｔｉｍｅ−Ｒｅｓｏｌｖｅｄ ｐａｒａｌｌｅｌ Ｉｍａｇｉｎｇ ｗｉｔｈ ｒｅｄｕｃｅｄ Ｄｙｎａｍｉｃ Ｆｉｅｌｄ ｏｆ Ｖｉｅｗ」；Ｍａｇｎ．Ｒｅｓｏｎ．Ｍｅｄ．２０１１；６５：１０６２−７４．には、心臓ＭＲＩ（ＣＭＲＩ）用の加速された動的画像収集のための“Ｐａｒａｌｌｅｌ Ｉｍａｇｉｎｇ ａｎｄ Ｎｏｑｕｉｓｔ ｉｎ Ｔａｎｄｅｍ”（ＰＩＮＯＴ）という名前の方法が記載されている。この方法は、パラレルイメージングのＳＰＡＣＥ−ＲＩＰ実装とＮｏｑｕｉｓｔの低減された視野（ｒＦＯＶ）によるイメージング方法とを組み合わせたもので、どちらも画像再構成に直接反転モデルを使用する。
ＰＩＮＯＴ再構成では、画像再構成を可能にするために、次のような線形システム方程式を解くことが求められる。
ここでＳ０，．．，ＳＣ−１は感度マップを表し、ＭＳは静的ＦｏＶを埋める領域からのフーリエ変換を表し、ＭＤは位置決めのたびごとに更新される領域からのフーリエ変換を表す。３Ｄ心臓ＭＲ画像の場合は、このシステム行列のサイズを管理することが難しくなり、画像再構成を行うシステム反転はＣＰＵに極めて負担がかかる。
ＵＳ２０１４／０３０３４８０には、人の患者の体に発生する動的事象の加速されたＭＲＩスキャンからの信号を処理する方法が記載されている。患者は自分の体の関連部分を含むＭＲＩ検査を受ける。スキャン時間の経過と共に実質的な変化がないボクセルが特定され、スキャン信号全体から差し引かれる。
本発明の目的は上述の制限を克服する方法及びシステムを提供することである。
この目的は独立請求項１及び１１の特徴により達成される。従属請求項は本発明の詳細な有利な実施形態を述べている。
本発明の様々な実施形態によれば、物体のＭＲ画像を生成する方法は、
− 物体の１つの動作フェーズに関する第１のデータセットを提供するステップと、
− 第１のデータセットから関心領域（ＲＯＩ）の第１の画像を生成するステップと、
− 動的領域及び静的領域を第１の画像内で特定するステップであって、これらの領域は一定の期間内に主に動的又は静的である、ステップと、
− 動的領域をマスクして外すことにより第１の画像を編集するステップと、
− 残りの静的領域を示す編集された第１の画像の逆フーリエ変換を実行するステップと、
− 物体の１つの動作フェーズに関する第２のデータセットを提供するステップと、
− 残りの静的領域の逆フーリエ変換されたものを第２のデータセットから差し引くステップと、
− 差し引かれた第２のデータセットにフーリエ変換を実行するステップと、
− 第１の画像の関心領域について低減された関心領域の第２の画像を生成するステップであって、低減された関心領域は動的領域を含む、ステップと、を含む。
− 物体の１つの動作フェーズに関する第１のデータセットを提供するステップと、
− 第１のデータセットから関心領域（ＲＯＩ）の第１の画像を生成するステップと、
− 動的領域及び静的領域を第１の画像内で特定するステップであって、これらの領域は一定の期間内に主に動的又は静的である、ステップと、
− 動的領域をマスクして外すことにより第１の画像を編集するステップと、
− 残りの静的領域を示す編集された第１の画像の逆フーリエ変換を実行するステップと、
− 物体の１つの動作フェーズに関する第２のデータセットを提供するステップと、
− 残りの静的領域の逆フーリエ変換されたものを第２のデータセットから差し引くステップと、
− 差し引かれた第２のデータセットにフーリエ変換を実行するステップと、
− 第１の画像の関心領域について低減された関心領域の第２の画像を生成するステップであって、低減された関心領域は動的領域を含む、ステップと、を含む。
好ましくは、特定することは第１の画像上で静的及び動的領域を手動で選択することにより実行される。例えば心臓周囲のボックスは動的領域として選択することができ、これによってボックス外側の領域は静的領域として選択される。
前述のように、動きがある物体は総期間内に複数の動作フェーズからなる動きを実行する。提供された第１及び第２のデータセットのデータは、ＭＲＩスキャナを使用したデータ収集により予め収集される。ステップ（ｇ）において、残りの静的領域の逆フーリエ変換を第２のデータセットから差し引くために、（ステップ（ｅ）で生成されると共に第１のデータセットに基づく）このフーリエ変換された画像は、第２のデータセットと同じサンプリング方式に従ってアンダーサンプリングされる。
本発明の思想は、概ね静的な環境にある動きがある物体の（予め収集された）データを提供して、データの冗長性及び高いＳＥＮＳＥ加速係数を利用して等方性解像度に到達することに基づく。この方法は、３Ｄ収集が上述のＰＩＮＯＴ再構成について想定する再構成問題に対処する新しいアプローチを用いる。ＰＩＮＯＴ収集と同様に、２つの異なるデータセットが区別され得る。第１のデータセットでは、静的及び動的な情報を収集するｋ空間ラインが収集される。第２のデータセットでは、ｋ空間ラインのサブグループが、一部の情報を最終的な再構成を行うことができるように完全なデータセットから回復することができるという考えに基づいて物体の動作フェーズ（運動フェーズ）ごとに収集される。ＰＩＮＯＴとの相違点であるこのアプローチでは、式（１）に記載されたような行列式全体を構築する必要がない。
本発明の好適な実施形態によれば、全ての画像又は少なくとも第２の画像は、ＳＥＮＳＥ再構成を使用することにより生成される。感度エンコーディング（ＳＥＮＳＥ）アプローチでは、複数の同時に動作される受信器コイルのアレイが信号収集に使用される。アレイ素子は通常、不均一性の強い互いに異なる空間感度を示す表面コイルである。ＳＥＮＳＥアプローチの基本原理は、コイル感度の影響を勾配エンコーディングと類似のエンコーディング効果と見なすことである。実際、感度効果は勾配エンコーディングに数学的に概ね類似している。しかし、勾配概念に優る重要な利点は、感度機構によって、受信器アレイの複数の異なる感度での同時エンコーディングが可能になることである。従って、スキャン時間の相当な節約が、順次的な勾配切り替えを並列感度エンコーディングで部分的に置換することにより達成され得る。
感度エンコーディングは、受信器感度が一般に線形場勾配によりフーリエ前処理に対して相補的なエンコーディング効果を有するという事実に基づいている。従って、複数の受信器コイルを並列に使用することにより、フーリエ撮像におけるスキャン時間が相当短縮され得る。感度エンコーディングされたデータからの画像再構成の問題は一般的なやり方で定式化され、任意のコイル設定及びｋ空間サンプリングパターンについて解決される。
本発明の別の好適な実施形態によれば、動きを実行する物体は心臓であり、ＭＲ画像は心臓ＭＲ画像であり、物体の動作フェーズは心臓領域からの心フェーズである。概ね静的な環境における動きがある物体は、データの冗長性及び高いＳＥＮＳＥ加速係数を利用して一回の息止めで心臓シネ画像の等方性解像度に到達する角のない環状ボリュームを有する胸部全体である。この収集及び再構成手法によって、息止め期間に数分の再構成時間で３Ｄ等方性の心臓画像を取得することができる。
本発明の更に別の好適な実施形態によれば、少なくとも１つの別の画像が、９つのステップ（ａ−ｉ）のうちのステップ（ｆ）〜（ｉ）をそれぞれ頻繁に繰り返すことにより生成される。これらの他の画像は、別の動作フェーズ（運動フェーズ）ＰＨ２〜ＰＨＮに関連する。
好ましくは、第１の画像はフルサンプル画像を使うことにより生成され、第２の画像（及び次の全ての他の画像）は、好ましくはＳＥＮＳＥ再構成を用いることにより生成される。
本発明の別の好適な実施形態によれば、第１のデータセット及び第２のデータセットに提供されたデータは、ＭＲＩ走査ユニット（ＭＲＩスキャナ）を使用したデータ収集により生成される。好ましくは、関心領域（ＲＯＩ）はＭＲＩ走査ユニットの視野（ＦｏＶ）により決定される。従って、低減された関心領域は低減された視野（ｒＦｏＶ）と呼ばれる。
好ましくは、複数の受信器コイルがデータ収集に使用される。好ましくは、編集された第１の画像は、複数の画像に分割され、各画像はｋ空間領域で逆フーリエ変換される。各コイルに対するフーリエ変換された画像は、第２のデータセットのサンプリングと同じサンプリング方式に従ってアンダーサンプリングされる。
本発明の別の好適な実施形態によれば、複数の画像の各画像は、逆フーリエ変換を行う前に受信器コイルのうち対応する受信器コイルのコイル感度により重み付けされる。
本発明の様々な実施形態によれば、以下のステップを実行する、物体のＭＲ画像を生成するＭＲＩシステムが確立される：
− 物体の１つの動作フェーズに関する第１のデータセットを提供するステップ；
− 第１のデータセットから関心領域の第１の画像を生成するステップ；
− 動的領域及び静的領域を第１の画像内で特定し、これらの領域は一定の期間内に主に動的又は静的である、ステップ；
− 動的領域をマスクして外すことにより第１の画像を編集するステップ；
− 残りの静的領域を示す編集された第１の画像の逆フーリエ変換を実行するステップ；
− 物体の１つの動作フェーズに関する第２のデータセットを提供するステップ；
− 残りの静的領域を含む編集された第１の画像の逆フーリエ変換を第２のデータセットから差し引くステップ；
− 差し引かれた第２のデータセットにフーリエ変換を実行するステップ；及び
− 第１の画像のＲＯＩについて低減された関心領域の第２の画像を生成し、低減された関心領域は動的領域を含む、ステップ。
− 物体の１つの動作フェーズに関する第１のデータセットを提供するステップ；
− 第１のデータセットから関心領域の第１の画像を生成するステップ；
− 動的領域及び静的領域を第１の画像内で特定し、これらの領域は一定の期間内に主に動的又は静的である、ステップ；
− 動的領域をマスクして外すことにより第１の画像を編集するステップ；
− 残りの静的領域を示す編集された第１の画像の逆フーリエ変換を実行するステップ；
− 物体の１つの動作フェーズに関する第２のデータセットを提供するステップ；
− 残りの静的領域を含む編集された第１の画像の逆フーリエ変換を第２のデータセットから差し引くステップ；
− 差し引かれた第２のデータセットにフーリエ変換を実行するステップ；及び
− 第１の画像のＲＯＩについて低減された関心領域の第２の画像を生成し、低減された関心領域は動的領域を含む、ステップ。
動きがある物体は、一定の期間内に複数の動作フェーズからなる動きを実行し、この期間は総収集時間である。
本発明の好適な実施形態によれば、システムは、周囲に動きがある物体、特に胸部の心臓のＭＲ画像を生成する上述の方法を実行するために確立される。
本発明の別の好適な実施形態によれば、ＭＲＩシステムは、（ｉ）複数の同時に動作される受信器コイルのアレイを備えた、データ記録用のＭＲＩ走査ユニットと、（ｉｉ）磁気共鳴撮像の画像収集のためのコンピュータによるデータ処理ユニットとを備える。ステップ（ａ）〜（ｉ）は、ＭＲＩシステムのデータ処理ユニットにより実行される。ステップ（ａ）及び（ｆ）に続くデータ収集ステップは、ＭＲＩ走査ユニットを使用することにより実行される。
本発明の様々な他の実施形態は、上述の方法を特に上述のＭＲＩシステムを使用することにより実行するコンピュータプログラム製品に関する。
本発明のこれら及び他の態様は、以下に説明される実施形態から明らかとなり、以下に説明される実施形態を参照して理解されるであろう。
以下の考察では、撮像される物体として心臓について言及する。しかし本発明は、他の器官等の他の物体にも適用可能である。心臓は単に例として選択されている。
図１は、静的情報及びＳＥＮＳＥ収集を同時に使用した全ての再構成ステップの略図を示す。この略図は、ＲＯＩ内に角のない環状ボリュームを有する胸部の心臓ＭＲ（ＣＭＲ）画像１０、２０を総期間内に生成する手順を表す。
手順は次のステップを含む：
ステップ１（Ｓ１）： 心臓部から（第１の）心フェーズの第１のデータセットを提供すること；
ステップ２（Ｓ２）： 第１のデータセットからＲＯＩの第１の画像１０を生成すること；
ステップ３（Ｓ３）： それぞれ主に総期間内において動的又は静的な動的領域１２及び静的領域１４を第１の画像１０内で特定すること；
ステップ４（Ｓ４）：動的領域１２をマスクすることにより第１の画像１０を編集すること；
ステップ５（Ｓ５）： 残りの静的領域１４を示す編集された第１の画像１６の逆フーリエ変換（ＦＦＴ−１）を実行すること；
ステップ６（Ｓ６）： 心臓部から（第２の）心フェーズに関する第２のデータセットを提供すること；
ステップ７（Ｓ７）： 残りの静的領域１４を含む編集された第１の画像１６の逆フーリエ変換（ＦＦＴ−１）を第２のデータセットから差し引くこと；
ステップ８（Ｓ８）： 差し引かれた第２のデータセット１８にフーリエ変換を実行すること；及び
ステップ９（Ｓ９）： 少なくとも動的領域１２を含む低減された関心領域（ｒＲＯＩ）の第２の画像２０を生成すること。
ステップ１（Ｓ１）： 心臓部から（第１の）心フェーズの第１のデータセットを提供すること；
ステップ２（Ｓ２）： 第１のデータセットからＲＯＩの第１の画像１０を生成すること；
ステップ３（Ｓ３）： それぞれ主に総期間内において動的又は静的な動的領域１２及び静的領域１４を第１の画像１０内で特定すること；
ステップ４（Ｓ４）：動的領域１２をマスクすることにより第１の画像１０を編集すること；
ステップ５（Ｓ５）： 残りの静的領域１４を示す編集された第１の画像１６の逆フーリエ変換（ＦＦＴ−１）を実行すること；
ステップ６（Ｓ６）： 心臓部から（第２の）心フェーズに関する第２のデータセットを提供すること；
ステップ７（Ｓ７）： 残りの静的領域１４を含む編集された第１の画像１６の逆フーリエ変換（ＦＦＴ−１）を第２のデータセットから差し引くこと；
ステップ８（Ｓ８）： 差し引かれた第２のデータセット１８にフーリエ変換を実行すること；及び
ステップ９（Ｓ９）： 少なくとも動的領域１２を含む低減された関心領域（ｒＲＯＩ）の第２の画像２０を生成すること。
好ましくは、特定することは第１の画像上で静的及び動的領域を手動で選択することにより実行される。例えば心臓周囲のボックスは動的領域として選択することができ、これによってボックス外側の領域は静的領域として選択される。
提供された第１及び第２のデータセットのデータはデータ収集により予め収集される。上述の総期間は総収集時間である。ステップ１（Ｓ１）に備えて、例えばフルＦｏＶでの一回の収集が折り重なりアーチファクトを除去するために行われる。ステップ６（Ｓ６）に備えて、心フェーズごとに低減されたＦｏＶ（ｒＦｏＶ）を更新するためにｋｙ−ｋｚ空間のアンダーサンプリング収集が行われる。フルＦｏＶという用語は（完全な）ＲＯＩに対応し、ｒＦｏＶという用語は低減されたＲＯＩに対応する。
当業者は、第２のデータセット（ｒＦｏＶ）が第１のデータセットと比較してアンダーサンプリングされる場合は、好ましくはステップ５ａがステップ５及びステップ７の間に適用されることを理解するであろう。ステップ５ａでは、編集された第１の画像１６の逆フーリエ変換（ＦＦＴ−１）は、第２のデータセットを提供（収集）する際に使用されたアンダーサンプリングと同じやり方でアンダーサンプリングされる。従って手順は次のステップを含む：
ステップ１（Ｓ１）： 心臓部から（第１の）心フェーズの第１のデータセットを収集すること；
ステップ２（Ｓ２）： 第１のデータセットからＲＯＩの第１の画像１０を生成すること；
ステップ３（Ｓ３）： 動的領域１２及び静的領域１４を第１の画像１０内で特定することであって、これらの領域は総期間内に主に動的又は静的であること；
ステップ４（Ｓ４）： 動的領域１２をマスクして外すことにより第１の画像１０を編集すること；
ステップ５（Ｓ５）： 残りの静的領域１４を示す編集された第１の画像１６の逆フーリエ変換（ＦＦＴ−１）を実行すること；
ステップ５ａ（Ｓ５ａ）： 編集された第１の画像１６の逆フーリエ変換（ＦＦＴ−１）をアンダーサンプリング戦略を用いてアンダーサンプリングし、これによってアンダーサンプリングされた第１のデータセット３０を提供すること；
ステップ６（Ｓ６）： 心臓部から（第２の）心フェーズに関する第２のデータセットを収集すること。ここで第２のデータセットは、第２のデータセットの収集中にアンダーサンプリング戦略を用いることにより第１のデータセットと比較してアンダーサンプリングされる；
ステップ７（Ｓ７）： アンダーサンプリングされた第１のデータセット３０を第２のデータセットから差し引くこと；
ステップ８（Ｓ８）： 差し引かれた第２のデータセット１８にフーリエ変換（ＦＦＴ）を実行すること；及び
ステップ９（Ｓ９）： 低減された関心領域（ｒＲＯＩ）の第２の画像２０を生成し、低減された関心領域は少なくとも動的領域１２を含むこと。
ステップ１（Ｓ１）： 心臓部から（第１の）心フェーズの第１のデータセットを収集すること；
ステップ２（Ｓ２）： 第１のデータセットからＲＯＩの第１の画像１０を生成すること；
ステップ３（Ｓ３）： 動的領域１２及び静的領域１４を第１の画像１０内で特定することであって、これらの領域は総期間内に主に動的又は静的であること；
ステップ４（Ｓ４）： 動的領域１２をマスクして外すことにより第１の画像１０を編集すること；
ステップ５（Ｓ５）： 残りの静的領域１４を示す編集された第１の画像１６の逆フーリエ変換（ＦＦＴ−１）を実行すること；
ステップ５ａ（Ｓ５ａ）： 編集された第１の画像１６の逆フーリエ変換（ＦＦＴ−１）をアンダーサンプリング戦略を用いてアンダーサンプリングし、これによってアンダーサンプリングされた第１のデータセット３０を提供すること；
ステップ６（Ｓ６）： 心臓部から（第２の）心フェーズに関する第２のデータセットを収集すること。ここで第２のデータセットは、第２のデータセットの収集中にアンダーサンプリング戦略を用いることにより第１のデータセットと比較してアンダーサンプリングされる；
ステップ７（Ｓ７）： アンダーサンプリングされた第１のデータセット３０を第２のデータセットから差し引くこと；
ステップ８（Ｓ８）： 差し引かれた第２のデータセット１８にフーリエ変換（ＦＦＴ）を実行すること；及び
ステップ９（Ｓ９）： 低減された関心領域（ｒＲＯＩ）の第２の画像２０を生成し、低減された関心領域は少なくとも動的領域１２を含むこと。
ＰＩＮＯＴ収集と同様に、図示された手順は２つの異なるデータセットを区別することができる。第１のデータセットでは、静的及び動的な情報を収集するｋ空間ラインが収集される。第２のデータセットでは、ｋ空間ラインのサブグループが、一部の情報を最終的な再構成を行うことができるように完全なデータセットから回復することができるという考えに基づいて心フェーズ（心臓フェーズ）ごとに収集される。ＰＩＮＯＴとの相違点であるこのアプローチでは、式（１）に記載されたような行列式全体を構築する必要がない。対照的に、再構成は３つの異なる段階に分かれる。
第１の段階： この再構成段階では、フル画像１０が上記第１のデータセットを使用して１つの心フェーズに対して生成される（Ｓ２）。この画像１０は、従来のＳＥＮＳＥ再構成又はフルサンプル画像を使用して生成され、信号精度を改善することができる。
第２の段階： この再構成段階では、フル再構成画像内の規定の動的領域は、全ての心フェーズで等しい静的領域１４からのみ情報を得るために０に設定される（Ｓ４）。この再構成段階では、画像１６は、各コイルのコイル感度により重み付けされ、ｋ空間領域で逆フーリエ変換される（Ｓ５）。各コイルに対するフーリエ変換された画像は、上記第２のデータセットと同じサンプリング方式に従ってアンダーサンプリングされる（Ｓ５及びＳ７間の未表示の矢印により示される）。最後に、静的領域から生成されたｋ空間ラインは、心フェーズごとに更新されたｋ空間ラインから差し引かれる（Ｓ７）。
第３の段階： この再構成段階では、第２の段階で生成された画像は、画像空間へとフーリエ変換され（Ｓ８）、低減されたＦｏＶから情報を取り出すのみで従来のＳＥＮＳＥ再構成を使用して再構成される（Ｓ９）。
このアプローチに従えば、ＳＥＮＳＥ再構成及びＮｏＱＵＩＳＴ流の再構成はかなり分離される。ＳＥＮＳＥ情報は第３の段階でのみ使用されるのに対し、ＮｏＱＵＩＳＴ流の情報は再構成の第２の段階でのみ使用される。更に、再構成の第３の段階では低減された領域のみが再構成され、３Ｄケースについての従来のＳＥＮＳＥ再構成と比較して再構成速度が向上する。
本発明は、図面及び前述の記載において詳細に図示及び説明されたが、このような図示及び記載は、説明的又は例示的であって限定するものではないと見なされるべきである。すなわち本発明は、開示された実施形態に限定されるものではない。開示された実施形態のその他の変形が、図面、本開示及び添付の請求項の検討から、請求項に係る発明を実施する当業者によって理解されて実現され得る。請求項において、「ｃｏｍｐｒｉｓｉｎｇ（含む、備える）」という単語は、他の要素又はステップを除外するものではなく、不定冠詞「ａ」又は「ａｎ」は、複数を除外するものではない。特定の手段が相互に異なる従属請求項に列挙されるという単なる事実は、これらの手段の組み合わせが有利に用いられないことを示すものではない。請求項における任意の参照符号は、本発明の範囲を限定するものと解釈されるべきではない。
本発明に繋がる研究は、助成契約第２４２０３８の下で欧州連合第７フレームワークプログラム（ＦＰ−７−ＨＥＡＬＴＨ−２００９）からの資金提供を受けた。
Claims (8)
- 一定の期間内に複数の動作フェーズを有する動きを実行する物体のＭＲ画像を生成する方法であって、当該方法は、
ＳＥＮＳＥアンダーサンプリング方式を用いて、複数の受信器コイルを使用して前記物体の１つの前記動作フェーズに関する第１のデータセットを収集するステップＳ１と、
ＳＥＮＳＥ再構成によって、前記第１のデータセットから関心領域の第１の画像を再構成するステップＳ２と、
動的領域及び静的領域を前記第１の画像内で特定するステップＳ３であって、これらの領域は前記一定の期間内に主に動的又は静的である、ステップＳ３と、
前記動的領域をマスクして外すことにより前記第１の画像を編集するステップＳ４と、
編集された前記第１の画像を、前記複数の受信器コイルのコイル感度により重み付けするステップＳ５と、
残りの前記静的領域を示す編集され重み付けされた前記第１の画像の逆フーリエ変換を実行するステップＳ６と、
複数の受信器コイルを使用して前記物体の１つの前記動作フェーズに関する第２のデータセットを収集するステップＳ７であって、収集中にアンダーサンプリングによって前記第２のデータセットが提供され、前記第１のデータセットの視野と比較して低減された視野が生じ、更に、前記第１のデータセットの収集に対するのと同じＳＥＮＳＥアンダーサンプリング方式を用いる、ステップＳ７と、
前記残りの静的領域を含む前記編集された第１の画像の逆フーリエ変換されたものを前記第２のデータセットから差し引くステップＳ８と、
差し引かれた前記第２のデータセットにフーリエ変換を実行するステップＳ９と、
前記第１の画像の前記関心領域について低減された関心領域の第２の画像を生成するステップＳ１０であって、前記低減された関心領域は動的領域を含む、ステップＳ１０と、
を含む、方法。 - 前記動きを実行する物体は心臓であり、前記ＭＲ画像は心臓ＭＲ画像であり、前記物体の前記動作フェーズは心臓領域からの心フェーズである、請求項１に記載の方法。
- 少なくとも１つの別の第２の画像が、前記ステップＳ１−Ｓ１０の１０個のステップのうちの前記ステップＳ７−Ｓ１０をそれぞれ頻繁に繰り返すことにより生成される、請求項１又は２に記載の方法。
- 前記編集された第１の画像は、前記複数の受信器コイルに対応する複数の画像に分割され、前記複数の画像の各画像はｋ空間領域で逆フーリエ変換される、請求項１に記載の方法。
- 一定の期間内に複数の動作フェーズを有する動きを実行する物体のＭＲ画像を生成するＭＲＩシステムであって、当該ＭＲＩシステムは、
ＳＥＮＳＥアンダーサンプリング方式を用いて、複数の受信器コイルを使用して前記物体の１つの前記動作フェーズに関する第１のデータセットを収集することと、
ＳＥＮＳＥ再構成によって、前記第１のデータセットから関心領域の第１の画像を再構成することと、
動的領域及び静的領域を前記第１の画像内で特定することであって、これらの領域は前記一定の期間内に主に動的又は静的である、特定することと、
前記動的領域をマスクして外すことにより前記第１の画像を編集することと、
編集された前記第１の画像を、前記複数の受信器コイルのコイル感度により重み付けすることと、
残りの前記静的領域を示す編集され重み付けされた前記第１の画像の逆フーリエ変換を実行することと、
複数の受信器コイルを使用して前記物体の１つの前記動作フェーズに関する第２のデータセットを収集することであって、収集中にアンダーサンプリングによって前記第２のデータセットが提供され、前記第１のデータセットの視野と比較して低減された視野が生じ、更に、前記第１のデータセットの収集に対するのと同じＳＥＮＳＥアンダーサンプリング方式を用いる、収集することと、
前記残りの静的領域を含む前記編集された第１の画像の逆フーリエ変換されたものを前記第２のデータセットから差し引くことと、
差し引かれた前記第２のデータセットにフーリエ変換を実行することと、
前記第１の画像の前記関心領域について低減された関心領域の第２の画像を生成することであって、前記低減された関心領域は動的領域を含む、生成することと、
を実行する、ＭＲＩシステム。 - 請求項１乃至４の何れか一項に記載の方法を実行する、請求項５に記載のＭＲＩシステム。
- 磁気共鳴撮像の画像生成のためのコンピュータによるデータ処理ユニットを含む、請求項５又は６に記載のＭＲＩシステム。
- 請求項１乃至４の何れか一項に記載の方法を実行するための、コンピュータプログラム。
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EP14192290 | 2014-11-07 | ||
| EP14192290.6 | 2014-11-07 | ||
| PCT/EP2015/072756 WO2016071054A1 (en) | 2014-11-07 | 2015-10-01 | Method and system for generating mr images of a moving object in its environment |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JP2017533064A JP2017533064A (ja) | 2017-11-09 |
| JP2017533064A5 JP2017533064A5 (ja) | 2018-11-08 |
| JP6730995B2 true JP6730995B2 (ja) | 2020-07-29 |
Family
ID=51893873
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2017542280A Active JP6730995B2 (ja) | 2014-11-07 | 2015-10-01 | 環境内に動きがある物体のｍｒ画像を生成する方法及びシステム |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US10416265B2 (ja) |
| EP (1) | EP3215863B1 (ja) |
| JP (1) | JP6730995B2 (ja) |
| CN (1) | CN107209242B (ja) |
| WO (1) | WO2016071054A1 (ja) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE102016218715A1 (de) * | 2016-09-28 | 2018-03-29 | Siemens Healthcare Gmbh | Verbesserte Erzeugung kombinierter Magnetresonanzbilder |
| DE102018203507B4 (de) * | 2018-03-08 | 2020-10-22 | Siemens Healthcare Gmbh | Verfahren zur Aufnahme und Rekonstruktion eines vierdimensionalen Dynamikbilddatensatzes mit einer Magnetresonanzeinrichtung, Magnetresonanzeinrichtung, Computerprogramm und elektronisch lesbarer Datenträger |
| JP2022512949A (ja) * | 2018-10-31 | 2022-02-07 | セントロ ナショナル デ インベスティガシオネス カルディオバスキュラレス カルロス テルセーロ （エフェ．エセ．ペ．） | 無症候性アテローム性動脈硬化症のバイオマーカー |
Family Cites Families (21)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0769151A1 (en) * | 1995-05-02 | 1997-04-23 | Koninklijke Philips Electronics N.V. | Method of and device for magnetic resonance imaging of objects |
| US5653233A (en) | 1995-08-11 | 1997-08-05 | The Board Of Trustees Of The Leland Stanford Junior University | Method and apparatus for improved temporal resolution in dynamic MRI |
| US6353752B1 (en) | 1999-05-14 | 2002-03-05 | Board Of Trustees Of The Leland Standford Junior University | Reduced field-of-view method for cine magnetic resonance imaging |
| US6528995B1 (en) * | 2001-09-10 | 2003-03-04 | Schlumberger Technology Corporation | Methods and apparatus for measuring flow velocity in a wellbore using NMR and applications using same |
| JP2005508238A (ja) * | 2001-11-07 | 2005-03-31 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 高速ダイナミック画像を形成する磁気共鳴方法 |
| JP4657710B2 (ja) * | 2002-05-13 | 2011-03-23 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 前情報により向上される動的な磁気共鳴イメージング |
| AU2003300665A1 (en) * | 2003-01-20 | 2004-08-13 | Koninklijke Philips Electronics N.V. | Sub-sampled moving table mri |
| EP1642149A1 (en) * | 2003-05-06 | 2006-04-05 | Koninklijke Philips Electronics N.V. | Undersampled magnetic resonance imaging |
| WO2005047919A1 (en) * | 2003-11-12 | 2005-05-26 | Invivo Corporation | Method for generating fast magnetic resonance images |
| DE10353342B4 (de) * | 2003-11-14 | 2008-07-17 | Siemens Ag | Verbesserte MRT-Bildgebung auf Basis konventioneller PPA-Rekonstruktionsverfahren |
| US20080262340A1 (en) * | 2004-05-14 | 2008-10-23 | Koninklijke Philips Electronics N.V. | Moving Table Mri |
| WO2006106448A2 (en) * | 2005-04-06 | 2006-10-12 | Koninklijke Philips Electronics, N.V. | Sense mr parallel imaging with continuously moving bed |
| JP2009533128A (ja) * | 2006-04-13 | 2009-09-17 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 動き補償を含む、連続的に移動する対象物のｍｒｉ |
| JP5037866B2 (ja) * | 2006-06-28 | 2012-10-03 | 株式会社日立メディコ | 磁気共鳴イメージング装置 |
| US20080292167A1 (en) * | 2007-05-24 | 2008-11-27 | Nick Todd | Method and system for constrained reconstruction applied to magnetic resonance temperature mapping |
| US10067213B2 (en) * | 2010-12-22 | 2018-09-04 | Koninklijke Philips N.V. | Parallel MRI method using calibration scan, coil sensitivity maps and navigators for rigid motion compensation |
| WO2013103791A1 (en) * | 2012-01-06 | 2013-07-11 | Cincinnati Children's Hospital Medical Center | Correlation imaging for multi-scan mri with multi-channel data acquisition |
| JP6073627B2 (ja) * | 2012-10-01 | 2017-02-01 | 東芝メディカルシステムズ株式会社 | 磁気共鳴イメージング装置及び画像処理装置 |
| US20140303480A1 (en) | 2013-04-04 | 2014-10-09 | General Electric Company | Mri with static tissue removal |
| US10534060B2 (en) * | 2013-06-10 | 2020-01-14 | Toshiba Medical Systems Corporation | Parallel MRI with spatially misregistered signal |
| US9797970B2 (en) * | 2015-04-14 | 2017-10-24 | The Board Of Trustees Of The University Of Illinois | Image domain segmented echo planar magnetic resonance imaging using a 2D excitation radiofrequency pulse |
-
2015
- 2015-10-01 WO PCT/EP2015/072756 patent/WO2016071054A1/en active Application Filing
- 2015-10-01 CN CN201580060220.6A patent/CN107209242B/zh active Active
- 2015-10-01 EP EP15771635.8A patent/EP3215863B1/en active Active
- 2015-10-01 JP JP2017542280A patent/JP6730995B2/ja active Active
- 2015-10-01 US US15/523,701 patent/US10416265B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| WO2016071054A1 (en) | 2016-05-12 |
| EP3215863A1 (en) | 2017-09-13 |
| US20170328974A1 (en) | 2017-11-16 |
| CN107209242B (zh) | 2020-07-03 |
| EP3215863B1 (en) | 2020-12-09 |
| JP2017533064A (ja) | 2017-11-09 |
| CN107209242A (zh) | 2017-09-26 |
| US10416265B2 (en) | 2019-09-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bustin et al. | Five‐minute whole‐heart coronary MRA with sub‐millimeter isotropic resolution, 100% respiratory scan efficiency, and 3D‐PROST reconstruction | |
| Pang et al. | Whole‐heart coronary MRA with 100% respiratory gating efficiency: self‐navigated three‐dimensional retrospective image‐based motion correction (TRIM) | |
| JP6467344B2 (ja) | ディクソン水脂肪分離を用いるプロペラ | |
| EP3384307B1 (en) | Removal of image artifacts in sense-mri | |
| US10809337B2 (en) | Reconstructing magnetic resonance images with different contrasts | |
| US11002815B2 (en) | System and method for reducing artifacts in echo planar magnetic resonance imaging | |
| CN113470139B (zh) | 一种基于mri的ct图像重建方法 | |
| JP6560023B2 (ja) | 磁気共鳴イメージング装置 | |
| JP6510273B2 (ja) | 磁気共鳴イメージング装置、磁気共鳴イメージング方法及び磁気共鳴イメージングプログラム | |
| CN107209238A (zh) | 具有对边带伪迹的抑制的并行多切片mr成像 | |
| JP2012505709A (ja) | 移動テーブルｍｒｉ装置及び方法 | |
| JP2017529960A (ja) | アーチファクト抑制を有するプロペラｍｒイメージング | |
| JP6730995B2 (ja) | 環境内に動きがある物体のｍｒ画像を生成する方法及びシステム | |
| JP7350568B2 (ja) | 磁気共鳴データ収集方法及び磁気共鳴イメージング装置 | |
| CN108780134B (zh) | 移动对象的磁化率映射的系统、方法和装置 | |
| Hutter et al. | Dynamic field mapping and motion correction using interleaved double spin-echo diffusion MRI | |
| Haji-Valizadeh et al. | Comparison of complex k-space data and magnitude-only for training of deep learning–based artifact suppression for real-time cine MRI | |
| US9709651B2 (en) | Compensated magnetic resonance imaging system and method for improved magnetic resonance imaging and diffusion imaging | |
| JP5156958B2 (ja) | 磁気共鳴イメージング装置 | |
| JP2017533064A5 (ja) | ||
| WO2017167937A1 (en) | Dynamic mr imaging with increased temporal and spatial resolution | |
| EP3772660A1 (en) | Magnetic resonance image reconstruction using cluster analysis | |
| JP6873134B2 (ja) | Ｓｅｎｓｅイメージングにおける画像アーチファクトの除去 | |
| JP2019205542A (ja) | 磁気共鳴イメージング装置、処理装置、および医用画像処理方法 | |
| JP7507793B2 (ja) | オフ共鳴アーチファクト補正を伴うスパイラルｍｒ撮像 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20170711 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20180927 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20180927 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20190620 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20190802 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20191025 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20200122 |
|
| 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: 20200610 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20200703 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6730995 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |