JP2016041384A - 磁気共鳴イメージング装置、時系列画像作成方法及びプログラム - Google Patents
磁気共鳴イメージング装置、時系列画像作成方法及びプログラム Download PDFInfo
- Publication number
- JP2016041384A JP2016041384A JP2016000049A JP2016000049A JP2016041384A JP 2016041384 A JP2016041384 A JP 2016041384A JP 2016000049 A JP2016000049 A JP 2016000049A JP 2016000049 A JP2016000049 A JP 2016000049A JP 2016041384 A JP2016041384 A JP 2016041384A
- Authority
- JP
- Japan
- Prior art keywords
- time
- magnetic field
- phase
- series
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims description 232
- 238000002595 magnetic resonance imaging Methods 0.000 title claims description 62
- 238000012545 processing Methods 0.000 claims abstract description 197
- 238000003384 imaging method Methods 0.000 claims abstract description 187
- 230000003068 static effect Effects 0.000 claims abstract description 168
- 238000012937 correction Methods 0.000 claims abstract description 101
- 230000008569 process Effects 0.000 claims description 130
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 90
- RGCLLPNLLBQHPF-HJWRWDBZSA-N phosphamidon Chemical compound CCN(CC)C(=O)C(\Cl)=C(/C)OP(=O)(OC)OC RGCLLPNLLBQHPF-HJWRWDBZSA-N 0.000 claims description 54
- 238000004364 calculation method Methods 0.000 claims description 27
- 238000005259 measurement Methods 0.000 claims description 10
- 238000000926 separation method Methods 0.000 abstract description 8
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 31
- 230000005540 biological transmission Effects 0.000 description 19
- 230000008859 change Effects 0.000 description 16
- 239000002872 contrast media Substances 0.000 description 15
- 238000010586 diagram Methods 0.000 description 12
- 230000006870 function Effects 0.000 description 8
- 238000001514 detection method Methods 0.000 description 7
- 230000014509 gene expression Effects 0.000 description 6
- 230000029058 respiratory gaseous exchange Effects 0.000 description 6
- 210000004185 liver Anatomy 0.000 description 4
- 230000036278 prepulse Effects 0.000 description 4
- 239000000126 substance Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000003252 repetitive effect Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 244000228957 Ferula foetida Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000004043 dyeing Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Images
Classifications
-
- 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/4828—Resolving the MR signals of different chemical species, e.g. water-fat 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/5601—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution involving use of a contrast agent for contrast manipulation, e.g. a paramagnetic, super-paramagnetic, ferromagnetic or hyperpolarised contrast agent
-
- 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
-
- 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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56545—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by finite or discrete sampling, e.g. Gibbs ringing, truncation artefacts, phase aliasing artefacts
-
- 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/546—Interface between the MR system and the user, e.g. for controlling the operation of the MR system or for the design of pulse sequences
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Signal Processing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Vascular Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biomedical Technology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Nonlinear Science (AREA)
Abstract
【解決手段】時系列で得られるTEの異なる複数のエコー信号から静磁場不均一マップを作成し、静磁場不均一マップを用いて水・脂肪分離画像を得る。この際、静磁場不均一マップに、位相又は周波数の空間方向の不連続を補正する処理(空間方向不連続補正処理)を行うとともに、時系列方向の不連続を補正する処理(時系列方向不連続補正処理)を行う。時系列方向の処理を行う際に、その位相または周波数オフセットの大きさを用いて、時系列方向不連続補正処理の要否を判断し、要と判断されたときに、時系列方向の不連続補正処理を部分的或いは全体的に行う。
【選択図】図4
Description
第1の実施形態は、2点Dixon法を用いた実施形態である。なお本実施形態の装置の構成は、上述した実施形態と同様であり、以下、適宜図1〜図3を参照する。
静磁場不均一マップ作成部1111は、画像データベース1104に格納された、TE1の画像S1とTE2の画像S2から静磁場不均一マップを算出する。静磁場不均一マップをPとし、TE2の画像S2からTE1の画像S1の位相を減算した位相減算画像をTとすると、静磁場不均一マップPは次式(1)、(2)により求まる。
前回の撮像番号(t-1)の撮像との間で、共鳴周波数のずれなどによる位相のオフセットを求める。静磁場不均一マップPを、前回の撮像のものと位相差分し、位相のオフセットαを算出する。ただし、最初の撮像番号の場合は本ステップはスキップする。
判定部1116により、ステップS902で算出した位相のオフセットαを用いて、時系列の位相アンラップ処理を実施するか否かの判定を行う。位相のオフセットαが大きい場合は、時系列の位相アンラップ処理を行っても、時系列方向の位相ラップが解消できないリスクが高くなるため、時系列の位相アンラップ処理は行わないほうがよい。そこで、たとえば位相のオフセットαが−π/2〜π/2ならば時系列の位相アンラップ処理を行う判定とし、ステップS904へ移行する。位相のオフセットαが−π/2未満、π/2を超える場合、即ちαの絶対値が大きい場合は、時系列の位相アンラップ処理を行わない判定とし、ステップS907へ移行する。ただし、最初の撮像番号の場合は、時系列の位相アンラップ処理が行えないため、ステップS907へ移行する。
不連続補正処理部1113(時系列不連続補正部1115)が、前回の撮像番号の静磁場不均一マップを用いて時系列の位相アンラップ処理を行う。図10はダイナミックイメージングの静磁場不均一マップのデータを示した図であり、ここでは3Dの静磁場不均一マップ701、702を示してる。図10のθは、位相アンラップ済みの静磁場不均一マップであり、位相が角度の情報で表されている。時系列方向の位相アンラップ処理は、時系列方向の座標(t-1)の位相アンラップ済みのピクセル(たとえばA)から、時系列方向の座標(t)の空間座標(xyz)が同じ未処理ピクセル(たとえばB)を隣接するピクセルとして位相アンラップ処理する。
時系列方向の位相アンラップ処理の後に、空間不連続補正部1114が空間的な位相アンラップ処理を実行することで残った未処理ピクセル(時系列方向の位相アンラップ処理をしていないピクセル)を位相アンラップする。未処理ピクセルは、前の撮像番号の隣接するピクセルが位相アンラップ済みでない場合や、時系列方向に隣接するピクセルの位相差が大きい場合(−π/2〜π/2の範囲外)に発生する。空間的な位相アンラップ処理は、時系列方向の位相アンラップ済みのピクセルと空間的に隣接する未処理ピクセルがある場合に、未処理ピクセルを位相アンラップ処理し、位相アンラップ処理が可能な未処理ピクセルが無くなるまで行う。
分離画像演算部1117は、最初に、位相アンラップ済みの静磁場不均一マップθ(x,y,z,t)を用いて、TE2の画像S2(x,y,z,t)を位相補正する。このとき、位相アンラップ済みの静磁場不均一マップθ(x,y,z,t)には、位相アンラップ処理されていないピクセルも含まれているため、位相アンラップ処理していないピクセルは外挿によって値を求める。また、ステップS901で位相の値は2倍にしているため、TE2の画像S2(x,y,z,t)の位相補正は、画像S2(x,y,z,t)から静磁場不均一マップθ(x,y,z,t)の半分の位相を引算し、位相補正画像S2’(x,y,z,t)を求める。即ち、虚数をi、指数関数をeで表すと、位相補正画像S2’(x,y,z,t)は次式(6)で求められる。
時系列の位相アンラップを行わない判定のときは、空間不連続補正部1114が、今回の撮像番号の静磁場不均一マップP(x,y,z,t)を用いて、空間的な位相アンラップ処理を行うことによって、位相アンラップ済みの静磁場不均一マップθ(x,y,z,t)を求める。このときの位相アンラップ処理は、処理を開始するピクセル(開始ピクセル)を決めて、x,y,z方向のいずれかに隣接するピクセルを順に位相アンラップ処理していく。開始ピクセルは、例えば静磁場不均一マップP(x,y,z,t)の絶対値が最大となるピクセルとし、開始ピクセルの位相アンラップ済みの静磁場不均一マップθ(x,y,z,t)の値には、同座標の静磁場不均一マップP(x,y,z,t)の位相を設定する。位相アンラップ処理するピクセルの順番は、位相アンラップ済みピクセルに隣接(空間的な隣接)する未処理ピクセルの内、信号値(静磁場不均一マップPの絶対値)の高い順に処理するなどの手法を組み合わせることによって位相ラップを低減できる。位相アンラップの方法はステップS905と同様である。
分離画像演算部1117が、位相アンラップ済みの静磁場不均一マップθ(x,y,z,t)、TE1の画像S1(x,y,z,t)、TE2の画像S2(x,y,z,t)を用いて、水画像W(x,y,z,t)と脂肪画像F(x,y,z,t)を作成する。処理内容はステップS906と同様である。
画像比較部1118が、時系列の画像を比較し、時系列間で画像を揃える処理を行うか否かの判定をする。ダイナミックイメージングの最初の撮像番号は、画像を揃える処理を行わない判定とし、ステップS911へ移行する。その他の撮像番号は画像を揃える処理を行う判定とし、ステップS910へ移行する。
前回の撮像の水画像W(x,y,z,t-1)、脂肪画像F(x,y,z,t-1)と、今回の撮像のステップS908で作成した水画像W(x,y,z,t)、脂肪画像F(x,y,z,t)との間で一致の度合いを比較する。比較の結果が一致の場合は、水画像W(x,y,z,t)、脂肪画像F(x,y,z,t)をそのまま画像送信部1107へ出力する。不一致の場合は、水画像W(x,y,z,t)と脂肪画像F(x,y,z,t)のデータを入れ替えた後、画像送信部1107へ出力する。また、不一致の場合は、静磁場不均一マップθにオフセットの位相ラップが発生しているため、次の撮像番号の位相アンラップ処理のために、静磁場不均一マップθの全ピクセルに2πを加算もしくは減算し位相ラップを除去する。
ステップS909で、時系列間で画像を揃える処理を行わない判定となった、最初の撮像番号の画像について、水・脂肪判定処理を行う。水・脂肪判定処理には、いくつかの方法がある。たとえば、TE1の画像S1とTE2の画像S2とのピクセル毎の信号比を求め、信号比のバラつきが大きい領域を占める画像を水画像とする方法がある。水画像の領域は分離した水画像W、脂肪画像Fを用いて、閾値処理などによって求める。その他の方法として、水画像W(x,y,z,t)と脂肪画像F(x,y,z,t)のヒストグラムを作成しその傾向から判定する(たとえばヒストグラムの分散が大きい方を脂肪画像とする)方法がある。
第2の実施形態は、3点Dixon法を用いた実施形態である。また、不連続補正処理として、第1の実施形態では位相回転量の静磁場不均一マップの位相アンラップ処理を行ったのに対し、最小二乗推定による反復演算を用いて周波数の静磁場不均一マップの不連続を補正する処理を行うことが特徴である。以下、第1の実施形態と異なる点を中心に本実施形態を説明する。なお本実施形態の装置の構成は、上述した実施形態と同様であり、以下、適宜図1〜図3を参照する。
静磁場不均一マップ作成部1111が、画像データベース1104に格納された、TE1の画像s1(x,y,z,t)とTE2の画像s2(x,y,z,t)および、前回の撮像番号(t-1)のTE1の画像s1(x,y,z,t-1)とTE2の画像s2(x,y,z,t-1の撮像との間で、共鳴周波数のずれなどによる位相もしくは周波数オフセットΨoffsetを式(12)により求める。位相と周波数は関連した物理量(位相=2π×周波数×時間)であるため、どちらを用いても良いが本実施形態では周波数を用いて説明する。
判定部1116が、周波数オフセットの値によって時系列処理を行うか否かを判断する。周波数オフセットΨoffsetが小さい場合は、時系列処理をすると判断しステップS1203へ移行し、周波数オフセットΨoffsetが大きい場合は、時系列処理をいないと判断しステップS1206へ移行する。大小の判断は、脂肪のケミカルシフトによる周波数(1.5Tのとき223Hz)を基準にして、たとえば、周波数オフセットΨoffsetが(−223/4〜223/4)のとき小さいとする。この判断は、第1の実施形態におけるステップS903と同様に、時系列処理により不連続補正処理が不正確になるのを回避するために行う。
不連続補正処理部1113(時系列不連続補正部1115)が、前回の撮像番号の静磁場不均一マップΨを用いて、時系列の処理を行う。この処理は最小二乗推定の反復演算により静磁場不均一マップを求める処理であり、ここでは領域拡張処理という。図13はダイナミックイメージングの静磁場不均一マップ1301,1302のデータを示した図である。図13の静磁場不均一マップは周波数Ψの情報で表されている。時系列方向の領域拡張処理は、時系列方向の座標(t-1)の領域拡張済みピクセル(たとえば1301のA)から、時系列方向の座標(t)の空間座標が同じ未処理ピクセル(たとえば1302のB)を隣接するピクセルとして領域拡張する。
時系列方向の領域拡張処理の後に、不連続補正処理部1113(空間不連続補正部1114)が、空間的な領域拡張処理を実行し、時系列方向に処理されていない未処理ピクセルを領域拡張処理する。未処理ピクセルは、前の撮像番号の隣接するピクセルが領域拡張済みでない場合に発生する。領域拡張済みのピクセルと空間的に隣接する未処理ピクセルがある場合に未処理ピクセルを領域拡張処理し、領域拡張処理が可能な未処理ピクセルが無くなるまで行う。
分離画像演算部1117は、領域拡張済みの静磁場不均一マップΨ(x,y,z,t)、TE1の画像s1(x,y,z,t)、TE2の画像s2(x,y,z,t)、TE3の画像s3(x,y,z,t)を用いて、水画像W(x,y,z,t)と脂肪画像F(x,y,z,t)を作成し、画像送信部1107へ出力する。このとき、領域拡張済みの静磁場不均一マップΨ(x,y,z,t)は、領域拡張処理されていないピクセルも含まれているため、領域拡張処理していないピクセルは外挿によって値を求める。水画像W(x,y,z,t)と脂肪画像F(x,y,z,t)は、ステップS1204で求めた静磁場不均一マップΨ(x,y,z,t)を式(13)に代入して、未知数W(x,y,z,t)、F(x,y,z,t)の線形方程式を最小二乗推定することによって求めることができる。
判定部1116による判定ステップS1202において、時系列の領域拡張処理を行わない判定のときは、空間的な領域拡張処理を行うことによって、領域拡張済みの静磁場不均一マップΨ(x,y,z,t)を求める。このときの領域拡張処理は、開始のピクセル決めて、x,y,z方向のいずれかに隣接するピクセルを順に領域拡張処理していく。開始ピクセルは例えば、TE1の画像s1、TE2の画像s2、TE3の画像s3の信号変化が大きいピクセルとする方法がある。これは、水と脂肪が混在するピクセルは信号変化が大きく、また水と脂肪が混在するピクセルは局所収束しにくいためこのような手法を用いる。領域拡張処理するピクセルの順番は、領域拡張済みピクセルに隣接(空間的な隣接)する未処理ピクセルの内、信号値(s1×s2×s3の絶対値)の高い順に処理するなどの手法を組み合わせることで、局所収束を低減することができる。
分離画像演算部1117が、領域拡張済みの静磁場不均一マップΨ(x,y,z,t)、TE1の画像s1(x,y,z,t)、TE2の画像s2(x,y,z,t)、TE3の画像s3(x,y,z,t)を用いて、水画像W(x,y,z,t)と脂肪画像F(x,y,z,t)を作成する。処理内容はステップS1205と同様である。
この処理は、第1の実施形態のステップS909と同様であり、画像比較部1118が、時系列の画像を比較し、時系列間で画像を揃える処理を行うか否かの判定をする。この判定はイナミックイメージングの最初の撮像番号は画像を揃える処理を行わない判定とし終了する。その他の撮像番号は画像を揃える処理を行う判定とし、ステップS1209へ移行する。
前回の撮像の水画像W(x,y,z,t-1)、脂肪画像F(x,y,z,t-1)と、今回の撮像のステップS1207で作成した水画像W(x,y,z,t)、脂肪画像F(x,y,z,t)との間で一致の度合いを比較する。比較の結果が一致の場合は、水画像W(x,y,z,t)、脂肪画像F(x,y,z,t)をそのまま画像送信部1107へ出力する。不一致の場合は、水画像W(x,y,z,t)と脂肪画像F(x,y,z,t)のデータを入れ替えた後、画像送信部1107へ出力する。また、不一致の場合は、静磁場不均一マップΨが局所収束しているため、次の撮像番号の領域拡張処理のために、静磁場不均一マップΨの全ピクセルに脂肪のケミカルシフトの周波数を加算もしくは減算する。比較の方法は、第1の実施形態で示したステップS910の方法と同じである。
第1の実施形態及び第2の実施形態のMRI装置では、時系列の位相アンラップ処理もしくは時系列の領域拡張処理をするか否かの判定を、静磁場不均一マップの位相オフセットもしくは周波数オフセットの値を用いて判定部1116が自動判定したが、本実施形態のMRI装置は、これらの判定を操作者が行うための手段を備える。第1及び第2の実施形態では、画像比較部1118により自動で画像の比較を行ったが、本実施形態のMRI装置は、これらの比較を操作者が行うための手段を備える。具体的には、MRI装置の操作部を介して、操作者の判定結果や選択を受け付ける。
本実施形態においても、装置の構成は前述の実施形態と共通するので、異なる特徴を中心に本実施形態を説明する。以下の説明では、前述の実施形態に用いた図面を適宜参照する。
図14に示した「4D Unwrapping」、「Start number」、「End number」パラメータを用いて時系列の位相アンラップの判定を行う。「4D Unwrapping」がONで撮像番号が「Start number」以上「End number」未満のとき本撮像番号の処理を終了する。「4D Unwrapping」がONで撮像番号が「End number」のとき、ステップS1502へ移行する。その他の場合は、ステップS1505へ移行する。
「Start number」から「End number」までのTE1の画像S1(x,y,z,t)とTE2の画像S2(x,y,z,t)から静磁場不均一マップP(x,y,z,t)を求める。計算式は式(1)と式(2)と同様である。
「Start number」から「End number」までの静磁場不均一マップP(x,y,z,t)を用いて、時系列と空間を合わせた4次元の位相アンラップ処理を行う。最初に開始ピクセルを決め、x,y,z,t方向のいずれかに隣接するピクセルを順に位相アンラップ処理し、位相アンラップ済みの静磁場不均一マップθ(x,y,z,t)を作成する。
「Start number」から「End number」までの、4次元の位相アンラップ済みの静磁場不均一マップθ(x,y,z,t)、TE1の画像S1(x,y,z,t)、TE2の画像S2(x,y,z,t)を用いて、「Start number」から「End number」までの水画像W(x,y,z,t)と脂肪画像F(x,y,z,t)を作成する。この処理内容は第1の実施形態のステップS906と同様であり、式(6)及び式(7−1)、(7−2)を用いて求める。
ステップS1505〜S1507は、4次元位相アンラップ処理もしくは4次元領域拡張を行わない場合であり、第1の実施形態のフロー(図9)のステップS901、S907、S908と同様である。すなわち、ステップS1505では、今回の撮像番号の、TE1の画像S1(x,y,z,t)とTE2の画像S2(x,y,z,t)を用いて、静磁場不均一マップP(x,y,z,t)を作成する。ステップS1506では、静磁場不均一マップP(x,y,z,t)を空間的な位相アンラップ処理を行うことによって、位相アンラップ済みの静磁場不均一マップθ(x,y,z,t)を求める。さらにステップS1507で、位相アンラップ済みの静磁場不均一マップθ(x,y,z,t)、TE1の画像S1(x,y,z,t)、TE2の画像S2(x,y,z,t)を用いて、水画像W(x,y,z,t)と脂肪画像F(x,y,z,t)を作成する。
一方、図14のUIにおいて、「Series matching」がONで、かつ「Maching number」に今回の撮像番号があるときは、ステップS1509へ移行し、その他の場合はステップS1510へ移行する。
図14のUI「Standard number」に入力された撮像番号(図示する例ではt=1)の水画像W(x,y,z,1)、脂肪画像F(x,y,z,1)と、今回の撮像のステップS1207で作成した水画像W(x,y,z,t)、脂肪画像F(x,y,z,t)を用いて、画像の比較と時系列方向の画像合わせを行う。この処理は第1の実施形態のステップS910と同様であり、画像の一致度から画像の入れ替えの要否を判断し、判断結果に従い入れ替えを実施或いは不実施とする。画像の入れ替えは今回の撮像番号が「End number」の場合には、「Start number」から「End number」までの画像をすべて入れ替える。
図14のUI「Series matching」がOFFの場合(或いはONが選択されない場合)は、時系列の画像の比較は行わずに、水・脂肪判定処理を行う。この処理は第1の実施形態のステップS911と同じである。
本実施形態のMRI装置によれば、操作者に時系列方向不連続補正処理を行うか否かの判定や選択を行わせる手段を備えたことによって、以下述べるように、さらに画像の精度と利便性をあげることができる。
本実施形態も、操作者が4次元位相アンラップ処理もしくは4次元領域拡張処理や時系列の画像比較を選択する手段を備えることは第3の実施形態と同様である。本実施形態は、操作者による選択をさらに利便化するため、操作者の判断の指標となる情報を提供する手段を備える。このような情報として、例えば共鳴周波数の時系列の変化がある。
一方、図14のUIにおいて、「Series matching」がONで、かつ「Matching number」に今回の撮像番号があるときは、ステップS1509へ移行し、その他の場合はステップS1510へ移行する。
図14のUI「Standard number」に入力された撮像番号(図示する例ではt=1)の水画像W(x,y,z,1)、脂肪画像F(x,y,z,1)と、今回の撮像のステップS1207で作成した水画像W(x,y,z,t)、脂肪画像F(x,y,z,t)を用いて、画像の比較と時系列方向の画像合わせを行う。この処理は第1の実施形態のステップS910と同様であり、画像の一致度から画像の入れ替えの要否を判断し、判断結果に従い入れ替えを実施或いは不実施とする。画像の入れ替えは今回の撮像番号が「End number」の場合には、「Start number」から「End number」までの画像をすべて入れ替える。
Claims (15)
- 静磁場におかれた被検体に対し高周波磁場及び傾斜磁場を印加し、1回の撮像においてエコー時間を異ならせた複数のエコー信号を収集する撮像部と、
前記エコー時間が異なる複数のエコー信号を用いて、前記エコー時間が異なる複数の画像を作成する信号処理部と、を備え、
前記信号処理部は、
前記複数の画像を用いて、静磁場不均一に起因して前記画像に生じる位相回転又は周波数のずれを表す静磁場不均一マップを作成する静磁場不均一マップ作成部と、
前記1回の撮像で取得した静磁場不均一マップに対し、空間方向に位相又は周波数の不連続を補正する不連続補正処理を行う空間方向不連続補正部と、
時系列の撮像で取得した複数の静磁場不均一マップに対し、時系列方向に位相又は周波数の不連続を補正する不連続補正処理を行う時系列方向不連続補正部と、を備え、
前記エコー時間が異なる複数の画像と、空間方向及び/又は時系列方向に不連続補正処理された前記静磁場不均一マップとを用いて、異なる組織コントラストを持つ複数種の時系列画像を作成することを特徴とする磁気共鳴イメージング装置。 - 請求項1に記載の磁気共鳴イメージング装置であって、
前記信号処理部は、時系列方向における静磁場不均一マップの位相又は周波数のオフセットが、所定の閾値の範囲内か否かを判定し、前記オフセットが所定の閾値の範囲内と判定したとき、前記時系列方向不連続補正部による時系列方向の不連続補正処理と前記空間方向不連続補正部による空間方向の不連続補正処理とを行うことを特徴とする磁気共鳴イメージング装置。 - 請求項2に記載の磁気共鳴イメージング装置であって、
前記信号処理部は、時系列方向における静磁場不均一マップの位相又は周波数のオフセットが、所定の閾値の範囲内か否かを判定し、前記オフセットが所定の閾値の範囲外と判定したときは、時系列方向の不連続補正処理を行うことなく、前記空間方向不連続補正部による空間方向の不連続補正処理を行うことを特徴とする磁気共鳴イメージング装置。 - 請求項1乃至3のいずれか一項に記載の磁気共鳴イメージング装置であって、
前記信号処理部は、前記複数種の時系列画像を、時系列方向に隣接する画像間で比較し、画像の一致度を判定する画像比較部をさらに備え、
前記画像比較部の判定結果に従い、同種の画像が揃った時系列画像を作成することを特徴とする磁気共鳴イメージング装置。 - 請求項1乃至4のいずれか一項に記載の磁気共鳴イメージング装置であって、
前記時系列方向不連続補正部が行う不連続補正処理は、時系列方向の位相アンラップ処理であることを特徴とする磁気共鳴イメージング装置。 - 請求項1乃至4のいずれか一項に記載の磁気共鳴イメージング装置であって、
前記時系列方向不連続補正部が行う不連続補正処理は、最小二乗推定の反復演算の初期値により周波数の不連続を補正する処理であることを特徴とする磁気共鳴イメージング装置。 - 請求項5又は6のいずれか一項に記載の磁気共鳴イメージング装置であって、
前記時系列方向不連続補正部は、基準となるピクセルを開始点として領域拡張処理により不連続補正処理を行うことを特徴とする磁気共鳴イメージング装置。 - 請求項1乃至7のいずれか一項に記載の磁気共鳴イメージング装置であって、
前記撮像部が行う撮像は、第1及び第2のエコー時間でそれぞれエコー信号を取得する2点ディクソン法による撮像であり、前記第1のエコー時間は、高周波磁場印加後に変化する水プロトンの位相と脂肪プロトンの位相とが揃う時間であり、前記第2のエコー時間は高周波磁場印加後に前記水プロトンの位相と前記脂肪プロトンの位相とが180度ずれる時間であることを特徴とする磁気共鳴イメージング装置。 - 請求項1乃至7のいずれか一項に記載の磁気共鳴イメージング装置であって、
前記撮像部が行う撮像は、第1、第2及び第3のエコー時間でそれぞれエコー信号を取得する3点ディクソン法による撮像であり、各エコー時間は、高周波磁場印加後に変化する水プロトンの位相と脂肪プロトンの位相との位相差が異なる時間に設定されていることを特徴とする磁気共鳴イメージング装置。 - 請求項1に記載の磁気共鳴イメージング装置であって、
操作者による時系列方向の不連続補正処理の選択を受け付ける操作部をさらに備え、
前記信号処理部は、前記操作部が時系列方向の不連続補正処理の選択を受け付けたときに、前記時系列方向不連続補正部による不連続補正処理を行うことを特徴とする磁気共鳴イメージング装置。 - 請求項10に記載の磁気共鳴イメージング装置であって、
操作者による時系列方向の不連続補正処理の範囲の指定を受け付ける操作部をさらに備え、
前記信号処理部は、前記操作部が受け付けた時系列方向の不連続補正処理の範囲で、前記時系列方向不連続補正部による不連続補正処理を行うことを特徴とする磁気共鳴イメージング装置。 - 請求項10に記載の磁気共鳴イメージング装置であって、
前記複数種の時系列画像を、時系列方向に隣接する画像間で比較し、画像の一致度を判定する画像比較部と、操作者による時系列方向の画像の比較の指示を受け付ける操作部とをさらに備えることを特徴とする磁気共鳴イメージング装置。 - 請求項10に記載の磁気共鳴イメージング装置であって、
前記撮像部は、1回の撮像毎に、共鳴周波数測定を行い、
前記操作部は、前記撮像部が測定した撮像毎の共鳴周波数を表示する表示部を備えることを特徴とする磁気共鳴イメージング装置。 - 磁気共鳴イメージング装置により時系列の撮像でそれぞれ取得した、エコー時間の異なるエコー信号から再構成した複数の原画像を用いて、原画像間の演算により、複数種の組織コントラスト画像を時系列画像として形成する時系列画像作成方法であって、
時系列の撮像毎に、前記複数の原画像から、静磁場不均一に起因する位相回転又は周波数のずれを表す静磁場不均一マップを作成する処理と、
時系列の撮像毎に得られる前記静磁場不均一マップに対し、位相又は周波数の空間方向の不連続を補正する処理と、
時系列の撮像毎に得られる前記静磁場不均一マップに対し、位相又は周波数の時系列方向のオフセットを所定の閾値と比較する処理と、
前記時系列方向の不連続が所定の閾値範囲であるときに、前記時系列方向の不連続を補正する処理と、
補正された前記静磁場不均一マップと前記複数の原画像とを用いて、複数種の組織コントラスト画像を作成する処理と、を含む時系列画像作成方法。 - コンピュータに、
時系列の撮像毎に、前記複数の原画像から、静磁場不均一に起因する位相回転又は周波数のずれを表す静磁場不均一マップを作成するステップと、
時系列の撮像毎に得られる前記静磁場不均一マップに対し、位相又は周波数の空間方向の不連続を補正するステップと、
時系列の撮像毎に得られる前記静磁場不均一マップに対し、位相又は周波数の時系列方向のオフセットを所定の閾値と比較するステップと、
前記時系列方向のオフセットが所定の閾値範囲であるときに、前記時系列方向の不連続を補正するステップと、
補正された前記静磁場不均一マップと前記複数の原画像とを用いて、複数種の組織コントラスト画像を作成するステップと、を実行させるプログラム。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016000049A JP5979327B2 (ja) | 2016-01-04 | 2016-01-04 | 磁気共鳴イメージング装置、その作動方法及び時系列画像作成プログラム |
US14/993,639 US9714999B1 (en) | 2016-01-04 | 2016-01-12 | Magnetic resonance imaging apparatus, time-series image generation method, and program |
CN201611106858.1A CN107019511B (zh) | 2016-01-04 | 2016-12-05 | 磁共振成像装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2016000049A JP5979327B2 (ja) | 2016-01-04 | 2016-01-04 | 磁気共鳴イメージング装置、その作動方法及び時系列画像作成プログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2016041384A true JP2016041384A (ja) | 2016-03-31 |
JP5979327B2 JP5979327B2 (ja) | 2016-08-24 |
Family
ID=55591502
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2016000049A Active JP5979327B2 (ja) | 2016-01-04 | 2016-01-04 | 磁気共鳴イメージング装置、その作動方法及び時系列画像作成プログラム |
Country Status (3)
Country | Link |
---|---|
US (1) | US9714999B1 (ja) |
JP (1) | JP5979327B2 (ja) |
CN (1) | CN107019511B (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10234521B2 (en) * | 2014-04-24 | 2019-03-19 | Junmin LIU | Systems and methods for field mapping in magnetic resonance imaging |
JP2019524303A (ja) * | 2016-08-15 | 2019-09-05 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | ディクソン型の水/脂肪分離を用いるmr撮像 |
JP2020022572A (ja) * | 2018-08-06 | 2020-02-13 | キヤノンメディカルシステムズ株式会社 | 磁気共鳴イメージング装置 |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6595393B2 (ja) * | 2016-04-04 | 2019-10-23 | 株式会社日立製作所 | 磁気共鳴イメージング装置、及び、画像処理方法 |
JP7220374B2 (ja) * | 2018-08-24 | 2023-02-10 | パナソニックIpマネジメント株式会社 | 通信端末および通信方法 |
JP7167564B2 (ja) * | 2018-09-05 | 2022-11-09 | 株式会社島津製作所 | X線撮影装置およびx線撮影装置の作動方法 |
JP7223619B2 (ja) * | 2019-03-29 | 2023-02-16 | キヤノンメディカルシステムズ株式会社 | 磁気共鳴イメージング装置、画像処理装置及び画像処理方法 |
CN110074786B (zh) * | 2019-04-30 | 2022-12-06 | 上海东软医疗科技有限公司 | 核磁共振匀场方法、装置、计算设备及核磁共振成像*** |
US12007458B2 (en) | 2019-07-08 | 2024-06-11 | Canon Medical Systems Corporation | Magnetic resonance imaging apparatus and magnetic resonance imaging method |
CN113466765A (zh) * | 2020-03-31 | 2021-10-01 | 通用电气精准医疗有限责任公司 | 磁共振扫描方法及***、计算机可读存储介质 |
DE102020208186A1 (de) * | 2020-06-30 | 2021-12-30 | Siemens Healthcare Gmbh | Erzeugung eines Homogenisierungsfeldes für Magnetresonanz-Bilddaten |
CN114115388B (zh) * | 2022-01-25 | 2022-06-24 | 康达洲际医疗器械有限公司 | 一种基于图像分析自适应的磁共振分时调控方法与*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002200056A (ja) * | 2000-12-28 | 2002-07-16 | Hitachi Medical Corp | 磁気共鳴イメージング方法及び装置 |
JP2012090867A (ja) * | 2010-10-28 | 2012-05-17 | Hitachi Medical Corp | 磁気共鳴イメージング装置及びケミカル画像取得方法 |
JP2013183967A (ja) * | 2012-03-09 | 2013-09-19 | Hitachi Medical Corp | 磁気共鳴イメージング装置及び位相補正方法 |
WO2013157299A1 (ja) * | 2012-04-19 | 2013-10-24 | 国立大学法人東京工業大学 | 信号処理装置、信号処理方法およびプログラム |
JP2015503432A (ja) * | 2012-01-12 | 2015-02-02 | コーニンクレッカ フィリップス エヌ ヴェ | B1マッピングのあるmr撮像 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH01303137A (ja) * | 1988-05-31 | 1989-12-07 | Yokogawa Medical Syst Ltd | 水,脂肪分離mri装置 |
JP2767258B2 (ja) * | 1988-09-22 | 1998-06-18 | ジーイー横河メディカルシステム株式会社 | 核磁気共鳴装置 |
JP2000342552A (ja) * | 1999-06-08 | 2000-12-12 | Hitachi Medical Corp | 磁気共鳴イメージング装置における静磁場不均一の補正方法 |
WO2002053031A1 (fr) * | 2000-12-28 | 2002-07-11 | Hitachi Medical Corporation | Appareil et procede d'imagerie par resonance magnetique |
DE10122874B4 (de) * | 2001-05-11 | 2004-09-23 | Siemens Ag | Verfahren zur Extraktion von Spinkollektiven mit unterschiedlicher chemischer Verschiebung aus phasenkodierten Einzelbildern unter Berücksichtigung von Feldinhomogenitäten sowie Vorrichtung dazu |
US7141972B2 (en) * | 2003-11-17 | 2006-11-28 | Toshiba America Mri, Inc. | Water-fat separation for fast spin echo imaging in an inhomogeneous field with progressive encoding |
US9430854B2 (en) * | 2012-06-23 | 2016-08-30 | Wisconsin Alumni Research Foundation | System and method for model consistency constrained medical image reconstruction |
DE102014222964A1 (de) * | 2014-11-11 | 2016-05-12 | Siemens Aktiengesellschaft | Dynamische Verzeichnungskorrektur von Magnetresonanz-Bildaufnahmen |
DE102015209838B4 (de) * | 2015-05-28 | 2017-04-13 | Siemens Healthcare Gmbh | Verfahren zu einer Korrektur von Magnetresonanz-Messdaten |
-
2016
- 2016-01-04 JP JP2016000049A patent/JP5979327B2/ja active Active
- 2016-01-12 US US14/993,639 patent/US9714999B1/en active Active
- 2016-12-05 CN CN201611106858.1A patent/CN107019511B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002200056A (ja) * | 2000-12-28 | 2002-07-16 | Hitachi Medical Corp | 磁気共鳴イメージング方法及び装置 |
JP2012090867A (ja) * | 2010-10-28 | 2012-05-17 | Hitachi Medical Corp | 磁気共鳴イメージング装置及びケミカル画像取得方法 |
JP2015503432A (ja) * | 2012-01-12 | 2015-02-02 | コーニンクレッカ フィリップス エヌ ヴェ | B1マッピングのあるmr撮像 |
JP2013183967A (ja) * | 2012-03-09 | 2013-09-19 | Hitachi Medical Corp | 磁気共鳴イメージング装置及び位相補正方法 |
WO2013157299A1 (ja) * | 2012-04-19 | 2013-10-24 | 国立大学法人東京工業大学 | 信号処理装置、信号処理方法およびプログラム |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10234521B2 (en) * | 2014-04-24 | 2019-03-19 | Junmin LIU | Systems and methods for field mapping in magnetic resonance imaging |
JP2019524303A (ja) * | 2016-08-15 | 2019-09-05 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | ディクソン型の水/脂肪分離を用いるmr撮像 |
JP2020022572A (ja) * | 2018-08-06 | 2020-02-13 | キヤノンメディカルシステムズ株式会社 | 磁気共鳴イメージング装置 |
JP7353735B2 (ja) | 2018-08-06 | 2023-10-02 | キヤノンメディカルシステムズ株式会社 | 磁気共鳴イメージング装置 |
Also Published As
Publication number | Publication date |
---|---|
JP5979327B2 (ja) | 2016-08-24 |
US9714999B1 (en) | 2017-07-25 |
CN107019511B (zh) | 2020-01-21 |
US20170192075A1 (en) | 2017-07-06 |
CN107019511A (zh) | 2017-08-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5979327B2 (ja) | 磁気共鳴イメージング装置、その作動方法及び時系列画像作成プログラム | |
JP5420206B2 (ja) | 撮像パラメータの決定方法、撮像パラメータの調整装置、コンピュータ読取可能な媒体および電子的に読取可能なデータ媒体 | |
US9329254B2 (en) | Magnetic resonance method and system to correct phase information in MR images | |
JP5735793B2 (ja) | Mr撮像において種信号を定量的に分離するシステム | |
RU2605524C2 (ru) | Магнитно-резонансная визуализация с подавлением артефактов потока | |
US10132899B2 (en) | Magnetic resonance imaging apparatus with voxel size control based on imaging parameters and magnetic susceptibility | |
US20140043022A1 (en) | Contrast enhanced magnetic resonance angiography with chemical shift encoding for fat suppression | |
JP2008508070A (ja) | 位相マッピングと、位相基準として用いる基準媒体が関係するmri温度測定 | |
US20170018080A1 (en) | Medical image diagnosis assistance device, magnetic resonance imaging apparatus and medical image diagnosis assistance method | |
JP2015019958A (ja) | 磁気共鳴イメージング装置、画像処理装置、および画像処理方法 | |
EP2762910A1 (en) | Susceptibility-weighted magnetic resonance imaging | |
KR101685377B1 (ko) | 화상해석장치, 화상해석방법 및 화상해석프로그램을 기록한 기록매체 | |
EP2495579B1 (en) | Method for generating 2D or 3D maps of MRI T1 and T2 relaxation times | |
Chung et al. | An improved PSF mapping method for EPI distortion correction in human brain at ultra high field (7T) | |
US10302726B2 (en) | Image reconstruction for MRI using multiplexed sensitivity encoding | |
JP2017528276A (ja) | 血流アーチファクトを抑制させたDixonMR撮像 | |
US10132902B2 (en) | Intrinsic navigation from velocity-encoding gradients in phase-contrast MRI | |
KR20100103533A (ko) | Pcmr 영상의 프로토콜을 자동으로 생성하는 방법 및 시스템 | |
JP6782681B2 (ja) | 磁気共鳴イメージング装置、撮像パラメータセット生成用演算装置及び撮像パラメータセット生成プログラム | |
JP4230875B2 (ja) | 磁気共鳴イメージング装置 | |
JP7128972B2 (ja) | 歪み補正を伴うepi mr画像 | |
JP6230882B2 (ja) | 磁気共鳴イメージング装置及びエコー時間設定方法 | |
JP5959888B2 (ja) | 磁気共鳴イメージング装置及び位相補正方法 | |
JP3907944B2 (ja) | 磁気共鳴イメージング方法及び装置 | |
Testud et al. | Direct magnetic field estimation based on echo planar raw data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20160118 |
|
A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20160118 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20160202 |
|
A975 | Report on accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A971005 Effective date: 20160330 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20160405 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20160418 |
|
TRDD | Decision of grant or rejection written | ||
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20160624 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20160628 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20160711 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5979327 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
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 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |