JP2004313551A - 心臓の運動表示方法及び画像診断装置 - Google Patents
心臓の運動表示方法及び画像診断装置 Download PDFInfo
- Publication number
- JP2004313551A JP2004313551A JP2003113390A JP2003113390A JP2004313551A JP 2004313551 A JP2004313551 A JP 2004313551A JP 2003113390 A JP2003113390 A JP 2003113390A JP 2003113390 A JP2003113390 A JP 2003113390A JP 2004313551 A JP2004313551 A JP 2004313551A
- Authority
- JP
- Japan
- Prior art keywords
- image
- display
- wall thickness
- myocardium
- marks
- 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
- 230000033001 locomotion Effects 0.000 title claims abstract description 110
- 238000000034 method Methods 0.000 title claims abstract description 42
- 230000000747 cardiac effect Effects 0.000 title abstract description 4
- 210000004165 myocardium Anatomy 0.000 claims abstract description 68
- 238000003384 imaging method Methods 0.000 claims description 5
- 238000002059 diagnostic imaging Methods 0.000 claims description 4
- 238000003745 diagnosis Methods 0.000 abstract description 10
- 230000002861 ventricular Effects 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 description 21
- 210000001519 tissue Anatomy 0.000 description 18
- 238000010586 diagram Methods 0.000 description 15
- 238000012545 processing Methods 0.000 description 14
- 238000004364 calculation method Methods 0.000 description 11
- 230000006870 function Effects 0.000 description 8
- 230000002107 myocardial effect Effects 0.000 description 6
- 238000011156 evaluation Methods 0.000 description 4
- 230000017531 blood circulation Effects 0.000 description 3
- 230000008602 contraction Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 210000000056 organ Anatomy 0.000 description 3
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 238000004040 coloring Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 208000028867 ischemia Diseases 0.000 description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 230000007659 motor function Effects 0.000 description 2
- 238000013077 scoring method Methods 0.000 description 2
- 230000008054 signal transmission Effects 0.000 description 2
- 208000007536 Thrombosis Diseases 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 210000004351 coronary vessel Anatomy 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 210000002837 heart atrium Anatomy 0.000 description 1
- 230000004217 heart function Effects 0.000 description 1
- 230000000302 ischemic effect Effects 0.000 description 1
- 210000004115 mitral valve Anatomy 0.000 description 1
- 208000031225 myocardial ischemia Diseases 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000008929 regeneration Effects 0.000 description 1
- 238000011069 regeneration method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Abstract
【解決手段】負荷状態が異なる心臓の断層を撮像してなる少なくとも2つの動画像が記憶された画像記憶部から一フレーム画像を読み出して表示部に表示する第1ステップと、該表示された前記一フレーム画像上において心筋を挟んで対向する少なくとも一対のマークを設定する第2ステップと、前記各マークが設定された部位の心筋の動きを前記動画像上で追跡する第3ステップと、前記心筋の動きに合わせて前記各マークの位置を変更し、変更後の前記一対のマーク間の距離により前記心筋の壁厚を求める第4ステップと、前記壁厚と前記壁厚の変化率と前記壁厚の変化速度の少なくとも1つの運動量を求める第5ステップと、前記各対のマーク位置における前記心筋の運動量を、該運動量に応じて輝度と色の少なくとも一方を変えた運動表示画像を前記表示部に表示する第6ステップとからなり、前記負荷状態が異なる心臓の断層像について前記第1ステップから第6ステップを実行することを特徴とする。
【選択図】 図1
Description
【発明の属する技術分野】
本発明は、超音波診断画像、磁気共鳴画像又はX線CT画像に適用される心臓の運動表示方法、その方法を用いた画像診断装置に関する。
【0002】
【従来の技術】
超音波診断装置、磁気共鳴イメージング(MRI)装置、及びX線CT装置等の画像診断装置は、いずれも被検体の検査部位に係る断層像などをモニタに表示して診断に供するものである。例えば、心臓や血管等の循環器系及びその他の動きのある臓器の場合、それらを構成する生体組織(以下、組織と総称する)の動きを断層像により観察して、それら臓器等の機能を診断することが行なわれている。
【0003】
特に、心臓などの運動機能を定量的に評価できれば、診断の精度が一層向上することが期待されている。例えば、従来、超音波診断装置により得られた画像から心壁の輪郭を抽出し、その心壁輪郭に基づいて心室等の面積、容積、それらの変化率等から心機能(心臓ポンプ機能)を評価したり、局所の壁運動を観察して診断することが試みられている(特許文献1)。また、ドプラ信号等の計測信号に基づいて組織の変位を計測して、例えば局所的な収縮又は弛緩の分布を撮像し、これに基づいて心室の運動が活性化している場所を正確に決定したり、あるいは収縮期の心室壁の厚さ(壁厚)を計測する等、組織の運動を定量的に計測する方法が提案されている(特許文献2)。さらに、時々刻々変化する心房や心室の輪郭を抽出して、その輪郭を画像に重ねて表示するとともに、これに基づいて心室等の容量を求める技術が提案されている(特許文献3)。
【0004】
【特許文献1】
特開平9−13145号公報
【特許文献2】
特表2001−518342号公報
【特許文献3】
米国特許第5322067号公報(USP5,322,067)
【0005】
【発明が解決しようとする課題】
しかしながら、上記の従来技術は、いずれも心臓の全体的な機能を評価するための手法にとどまり、心筋などの各組織の動きを定量的に評価することについては配慮されていない。特に、心壁の輪郭を画像処理により抽出し、その輪郭に基づいて心室壁の厚み(壁厚)を計測したり、壁厚の変化を計測する従来技術は、必ずしも十分な精度を得るまでには至っていない。
【0006】
一般に、例えば、血栓等によって心筋に血が通わなくなると、心筋の動きが低下するといわれている。したがって、心室を構成する心筋の動きや壁厚の変化など、心臓の各組織の運動を定量的に計測できれば、治療法などを決定する際の有効な診断情報を提供できる。例えば、虚血の程度がわかれば、冠動脈再生術などの心臓の治療法選択及び治療部位を特定する指標として有効である。
【0007】
特に、心臓は収縮と拡張を繰返して血液を循環させることから、心筋壁厚の変化の度合いに点数をつけて診断するスコアリング法が提案されている。このスコアリング法は、例えば、心壁の延在方向に心壁を6分割し、動画像を観察しながら各分割領域ごとの壁厚の動きの程度を5段階に分けて評価する方法である。また、このスコアリングを、通常の状態と、患者に運動させて心臓に負荷をかけた負荷I状態と、さらに大きな負荷をかけた負荷II状態と、負荷を取り除いて回復した状態の4つの態様について実施して、心臓の診断を行なうことも行われている。
【0008】
しかし、従来のスコアリングは、断層像の動画像により心臓の動きを観察した観者の主観的な判断で評価せざるを得ないから、スコアリング結果が必ずしも客観的でなく、かつ的確な診断を行なうには相当の経験が要求されるという問題がある。
【0009】
そこで、本発明の発明者らは、心臓の断層像が表示されたモニタ上で、心筋壁の内外に一対の追跡点(マーク)を設定し、動画像上で追跡点を含む切り出し画像の移動先を相関法などの画像処理によって検出し、心筋壁の動きに合わせて追跡点を移動させて重畳表示する一方、その追跡結果に基づいて一対の追跡点間の距離(壁厚)等を定量的に計測することを提案している(特願2002−266864号)。
【0010】
しかし、定量的に計測した心室壁厚のデータを表示する具体的な方法についての配慮が十分でないことから、観者が的確な診断を行なう上で改善する点があった。
【0011】
本発明の課題は、定量的に計測した心室壁厚の表示法を改善して、的確なスコアリングを行えるようにすることにある。
【0012】
【課題を解決するための手段】
上記課題を解決する本発明の心臓の運動表示方法は、負荷状態が異なる心臓の断層を撮像してなる少なくとも2つの動画像が記憶された画像記憶部から一フレーム画像を読み出して表示部に表示する第1ステップと、該表示された前記一フレーム画像上において心筋を挟んで対向する少なくとも一対のマークを設定する第2ステップと、前記各マークが設定された部位の心筋の動きを前記動画像上で追跡する第3ステップと、前記心筋の動きに合わせて前記各マークの位置を変更し、変更後の前記一対のマーク間の距離により前記心筋の壁厚を求める第4ステップと、前記壁厚と前記壁厚の変化率と前記壁厚の変化速度の少なくとも1つの運動量を求める第5ステップと、前記各対のマーク位置における前記心筋の運動量を、該運動量に応じて輝度と色の少なくとも一方を変えた運動表示画像を前記表示部に表示する第6ステップとからなり、前記負荷状態が異なる心臓の断層像について前記第1ステップから第6ステップを実行することを特徴とする。
【0013】
このような本発明の心臓の運動表示方法によれば、心筋の動き、すなわち壁厚の変化等が、数値でなく輝度や色の変化としてマークの位置に対応付けて表示されるから、輝度や色に基づいて一目で心筋の動きに関する各分割領域のスコアを評価することができる。したがって、観測者ごとに主観的に定性評価する場合に比べて、評価のバラツキを低減できる。
ここで、各マークが設定された部位の心筋の動きをそれぞれ動画像上で追跡する上記のステップは、前記マークを含むサイズの切出し画像を前記一フレーム画像に設定するステップと、前記動画像の他のフレーム画像を検索して前記切出し画像と画像の一致度が最も高い同一サイズの局所画像を抽出するステップと、該一致度が最も高い局所画像と前記切出し画像の座標差に基づいて前記マークの移動先を求めるステップとを含んで構成することができる。この場合において、抽出された局所画像を新たな切出し画像とし、動画像のさらに他のフレーム画像に対して一致が最も高い局所画像を抽出するステップと、一致度が最も高い局所画像と切出し画像の座標差に基づいてマークの移動先を求めるステップとを繰り返し実行して、前記マークの移動先を順次求めることにより、マークを心筋の動きに追従させて移動表示することもできる。
【0014】
また、本発明の画像診断装置は、負荷状態が異なる心臓の断層を撮像して得られる少なくとも2つの動画像データが記憶される画像記憶手段と、前記心臓の断層像を表示する表示部と、前記画像記憶手段に記憶された前記動画像データを読み出して動画像又は静止画像を前記表示部に表示させる表示制御手段と、前記表示部に表示された前記静止画像上の心筋を挟んで対向する少なくとも一対のマークを設定する操作部と、前記表示部に表示された静止画像に前記各マークを含む切出し画像をそれぞれ設定し、前記表示部に表示される前記動画像に基づいて前記切出し画像の移動先を追跡する切り出し画像追跡手段と、該切出し画像の移動先に基づいて前記各マークの移動量を求める移動量演算手段と、該手段により求めた各マークの移動量に基づいて前記一対のマーク間の距離を求めて前記心筋の壁厚と該壁厚の変化率と該壁厚の変化速度の少なくとも1つの運動量を求める運動算出手段と、該手段により求めた前記運動量に応じて輝度と色の少なくとも一方の画素値を変えて前記表示部の表示領域に表示する運動表示画像を生成する画像生成手段と、前記負荷状態が異なる心臓の各断層像について前記運動表示画像を生成させて前記表示部に表示させる制御手段とを備えて構成することができる。
【0015】
本発明において、各負荷状態の運動表示画像を表示画面に並べて表示することいが好ましいが、各負荷状態の運動表示画像を時系列的に順次表示するようにしてもよい。また、心臓の1周期を複数の時相に分割し、負荷状態が異なる心臓の運動表示画像の時相を合わせて並べて表示することが望ましい。さらに、本発明の運動表示画像に時間軸を合わせて、心電計により計測される心電図や、心音計により計測される心音波形を並べて表示することが好ましい。
【0016】
具体的な運動表示画像は、矩形領域の一方の軸を前記心筋に沿って配置された複数対のマークの設定位置に対応させ、他方の軸を時間軸に対応させて表示することができる。この場合、矩形領域の時間と複数対のマークの設定位置の少なくとも一方を指定する線を入力設定することにより、指定線に沿った運動量の変化を表すグラフを表示部に表示するステップを設ける。
【0017】
また、運動表示画像の他の具体例としては、各対のマーク位置における動画像の心筋に対応する領域の輝度と色の少なくとも一方を当該位置における運動量に応じて変えた画像、あるいは心筋の内壁と外壁を模擬してなる模式図形に、該模式図形における各対のマーク位置に対応する領域の輝度と色の少なくとも一方を当該位置における運動量に応じて変えた画像を重ねた画像とすることができる。特に、心筋の模式図形に重ねて運動表示画像を表示することにより、心筋の動きが悪い部位を直視的に診断できる。
【0018】
【発明の実施の形態】
(実施の形態1)
本発明の心臓の運動表示方法を適用してなる一実施の形態の画像診断装置について、図1〜図4を用いて説明する。図1は本実施形態の心臓の運動表示方法の手順を示し、図2は図1の心臓の運動表示方法を適用してなる画像診断装置のブロック構成図である。図2に示すように、画像診断装置は、被検体である生体の断層を撮影してなる動画像が格納される画像記憶部1と、画像を表示する表示部2と、各種の指令を入力する操作卓3と、表示部2に表示される動画像の心臓のの動きを追跡する自動追跡部4と、自動追跡部4の追跡結果に基づいて又は運動算出部6の計測情報に基づいて運動表示画像を生成して記憶する画像生成記憶部5と、自動追跡部4の追跡結果に基づいて各種の計測情報を算出する運動算出部6と、これらを接続してなる信号伝送路7を含んで構成されている。
【0019】
画像記憶部1には、破線で示した診断画像撮像装置8から被検体の断層像を撮影してなる動画像がオンライン又はオフラインで格納されるようになっている。診断画像撮像装置8としては、超音波診断装置、磁気共鳴イメージング(MRI)装置及びX線CT装置等の診断装置が適用可能である。また、本実施の形態では、被検体の異なる運動負荷における心臓の断層像を撮像した例えば4状態の動画像が格納される。つまり、通常の負荷状態Aと、患者に運動をさせて心臓に負荷をかけた負荷I状態Bと、さらに大きな負荷をかけた負荷II状態Cと、負荷を取り除いて回復した状態Dの4つの態様についての動画像が格納される。また、それらの断層像は、図示していない心電計から出力される心電波形に同期させて撮像し、その心電波形を断層像に同期させて記憶しておくことが好ましいが、これに限らず時間軸を同期させて記憶しておくようにする。
【0020】
操作卓3は、表示部2に動画像の一フレーム画像(静止画像)を表示させる指令、表示部2に表示された静止画像上で動きを追跡したい生体組織の部位にマーク(目印)を重畳表示させる指令、表示部2に表示する画像の種類を選択する指令、等の各種指令を入力可能に形成されている。マークを設定する場合、4つの異なる負荷状態A〜Dの静止画像を4つ並べて表示してもよく、時間をずらして別々に表示してもよい。
【0021】
自動追跡部4は、画像診断装置全体を制御する制御手段10と、表示部2に表示する画像を切り替え制御する表示制御手段11と、表示部2に表示された一フレーム画像のマークの位置に対応する追跡部位を含むサイズの切出し画像を設定する切出し画像設定手段12と、画像記憶部1から動画像の他のフレーム画像を読み出して、切出し画像と画像の一致度が最も高い同一サイズの局所画像を抽出する切出し画像追跡手段13と、一致度が最も高い局所画像と切出し画像の座標差を求める移動量演算手段14と、その座標差に基づいて追跡部位の移動先座標を求める移動追跡手段15とを備えて構成されている。この移動追跡手段15により求められた追跡部位(マーク)の移動先は、動画像に合わせて順次記憶して動画像の追跡画像を生成するようにしてもよい。
【0022】
一方、運動算出部5は、自動追跡部4で求められた追跡部位の移動先座標に基づいて、追跡部位の移動量、移動速度、移動方向等の動きに関する物理量である計測情報を定量的に求めるとともに、これらの計測情報の変化を線図で表示部2に表示させる機能を有して構成されている。特に、本実施の形態の運動算出部5は、追跡部位の移動量に基づいて心筋の内外壁に対向させて設定された一対のマーク間の距離を心筋の壁厚として求めるとともに、心筋の壁厚の変化と、壁厚の変化率と、壁厚の変化速度等の運動量を求める機能を有して構成されている。
【0023】
また、画像生成記憶部9は、運動算出部5により求められた心筋の運動量を、その運動量に応じて輝度と色の少なくとも一方を変えた運動表示画像を生成するとともに記憶する機能を有して構成されている。特に、運動表示画像は4つの異なる負荷状態A〜Dごとに作成する。また、心電計で計測された心電波形を運動表示画像の同期させて表示するように構成されている。
【0024】
表示制御手段11は、操作卓3から入力される指令に応じて、画像記憶部1に記憶された動画像データを読み出して動画像又は静止画像と、画像生成記憶部5に記憶されている運動表示画像を切り替えて、又は重ねて表示する機能を有している。
【0025】
次に、本実施形態の画像診断装置の詳細な機能構成について、図1に示した処理手順に従って動作とともに説明する。まず、心臓の心筋の動きの追跡動作は、操作卓3から組織の動き追跡モードを選択する指令が入力されることによって開始する(S1)。表示制御手段11は、画像記憶部1から動画像の最初のフレーム画像ft(t=0)を読み出して表示部2に表示させる(S2)。例えば、最初のフレーム画像f0として図3に示す心臓の心室21の断層像が表示されたものとする。図3において、操作者が動きを追跡したい生体組織の追跡部位として、心筋22の特定の部位を選択したい場合、操作者は操作部3のマウスなどを操作してフレーム画像f0に重ねて追跡部位を設定するためのマークである追跡点23を表示させる。そして、その追跡点23を移動操作して所望の追跡部位に重畳表示させて追跡部位を入力設定する。本発明の特徴である心筋の壁厚を計測する場合は、図6(a)に示すように、心筋23を挟んで対向する一対の追跡点23を入力設定する。なお、図3において、符号24は僧帽弁である。
【0026】
追跡点23が入力設定されると、制御手段10はフレーム画像f0上の追跡点23の座標を取込み、切出し画像設定手段12に送る(S3)。切出し画像設定手段12は、図4(a)に示す様に、各追跡点23の画像を中心として、縦横2(A+1)画素(但しAは自然数)のサイズの矩形領域を切出し画像25として設定する(S4)。ここで、切出し画像25のサイズは、追跡点23の生体組織とは異なる生体組織を含む大きさの領域に設定することが好ましい。
【0027】
切出し画像追跡手段13は、画像記憶部1から動画像の次のフレーム画像f1を読み出し、切出し画像25と画像の一致度が最も高い同一サイズの局所画像を抽出する(S5)。この抽出処理は、いわゆるブロックマッチング法又は相関法と称される画像処理である。この抽出処理をフレーム画像f1の全領域について行なうと、処理時間がかかり過ぎる。そこで、抽出処理時間を短縮するため、本実施形態では、フレーム画像f1よりも十分に小さい、図4(b)に示す検索領域26について行なうようにしている。つまり、検索領域26は、切出し画像25に対して上下左右に一定の振り幅の画素数Bを付加した矩形領域とする。この画素数Bは、追跡部位に係る組織の移動量よりも大きく、例えば3〜10画素に設定する。これは、心臓などの循環器系の動く範囲は、通常の視野において、狭い領域に限られるからである。このようにして、検索領域26内の同一サイズの局所画像27を順次ずらして切出し画像25との画像の一致度を求める。図5に、相関法による画像追跡処理の具体例を示す。同例は、説明を簡単にするために、切出し画像25のサイズを矩形の9画素領域とし、検索領域26についても矩形の25画素領域として説明する。つまり、同図(a)に示す切出し画像25は、追跡点23の画素を中心としてA=1画素に設定した例であり、同図(b)に示す検索領域26はB=1画素に設定した例である。これによれば、同図(b)に示す様に、9個の局所領域27について画像の一致度を求めることになる。画像の一致度を求める方法については、公知の種々の方法を適用できる。
【0028】
次に、検索した複数の局所画像27の内で画像の一致度が最も高い局所画像27maxを抽出し、局所画像27maxを切出し画像25の移動先とし、局所画像27maxの座標を求める(S6)。これらの画像の座標は、中心画素の座標、あるいは矩形領域の何れかの角の座標で代表する。そして、局所画像27maxと切出し画像25の座標差を求め、これに基づいて追跡点23の移動先座標を求めて記憶する。必要に応じて、表示制御手段11は、追跡点23の移動先座標に基づいて表示部2のフレーム画像f1に追跡点23のマークを重ねて表示する(S7)。なお、局所画像27maxと切出し画像25における追跡点23の相対位置は変化しないものとして扱っている。
【0029】
運動算出部6は、S7で求められた追跡点23の移動先座標に基づいて、追跡点23の動き、つまり追跡部位の組織の動きに関する各種の計測情報を算出して記憶する(S8)。すなわち、移動前後の追跡点23の座標に基づいて、移動方向及び移動量を定量的に計測することができる。本実施の形態では一対の追跡点23間の距離から心筋23の壁厚を算出し、壁厚の変化、変化率、変化速度を指令に応じて算出し、記憶する(S9)。これに代えて、又は加えて、追跡部位の移動量、移動速度、移動方向等の動きに関する物理量である計測情報を定量的に求めることができる。そして、求められた追跡点23の移動に関する各種の計測情報、及びその変化をグラフで表示部に表示させることができる。これにより、観者は、追跡部位の動きを容易に観察することができる。
【0030】
次に、ステップS10に進み、動画像の全てのフレーム画像について追跡点23の追跡が終了したか否か判断し、未処理のフレーム画像があれば、ステップS5に戻ってS5〜S10の処理を繰り返す。全てのフレーム画像について追跡点23の追跡が終了した場合は、追跡処理動作を終了し、追跡点23の移動履歴のデータを記憶する。次いで、画像生成記憶部5は、指定された壁厚の変化、変化率、変化速度の運動表示画像を生成して表示部に表示させる。この表示画像の具体的な例は、後述する。このように、本実施形態によれば、追跡点23の移動先の座標を順次求めることができるから、追跡部位の動きを定量的に、かつ精度よく簡単に計測することができ、診断の情報を的確に提供することができる。
【0031】
次に、本発明の特徴部である心臓の運動表示方法の具体例について図6〜図11を用いて説明する。まず、操作卓3から指令を入力して表示部2に心臓の断層像を表示させる。このとき、図6(a)に示すように、運動負荷が異なる4つの断層像の静止画を表示させる。この静止画上で心筋22に沿って、心壁を挟んで複数対(1〜n)の追跡点23を設定し、動画像に切り替えて追跡処理を実行させる。なお、図6(a)では1つの断層像の静止画のみを示しているが、必要な場合は図7(a)に示すように4状態A〜Dの断層像の同一時相の静止画を並べて表示し、それぞれの状態A〜Dにそれぞれ複数対(1〜n、但しnは自然数)の追跡点23を設定することもできる。このようにして複数対(1〜n)の追跡点23の設定が終了すると、自動追跡部4が動作して、状態A〜Dごとに各追跡点23の移動を追跡する。運動算出部6は、各追跡点23の移動量の計測結果に基づいて、各対(1〜n)の追跡点23間の距離(壁厚)を計測するとともに、図6(b)、(c)に示すように、任意の一対の追跡点23における壁厚の変化と変化率を計測し、画像生成記憶部5において壁厚の変化と変化率のグラフを作成して、表示部2に表示させるようにすることができる。これにより、心臓の伸縮及び拡張に伴う心筋の動きを定量的に把握することができる。
【0032】
ところが、複数対(1〜n)の追跡点23の相互間における壁厚の変化や変化率を対比するには、n対の追跡点23ごとに計測される多数のグラフを並べて表示しなければならないから、煩雑になり実用的でない。つまり、心臓の運動機能を6分割した領域についてスコアリングを実施する場合、心筋22の各分割領域における壁厚の変化の平均値をグラフで表示するとしても、6つのグラフを並べて表示することになるから、対比観察が煩雑になる。
【0033】
そこで、図8(a)に示すように、各状態A〜Dに対応させて壁厚の変化を3次元的にグラフ表示をすることが考えられる。なお、同図において、横軸は時間を表し、縦軸は壁厚変化率を表し、奥行軸は心筋22に沿った計測部位を表しており、心臓が最も拡張したときに計測を開始して収縮期から再び拡張するまでの一周期について計測した結果である。このように、3次元画像で表示できれば、ある程度見やすくなるが、さらに見やすくして診断を支援するため、本実施の形態では、心筋の壁厚の変化(運動量)を時間軸に沿って輝度変化ないしカラー化した運動表示画像を生成して表示することを特徴とする。
【0034】
運動表示画像の一例としては、図8(b)に示すように、表示部2の表示画面に矩形の表示領域30を設定し、その横軸を時間軸とし、縦軸を計測部位として、その表示領域に各計測部位の心筋の壁厚の変化(運動量)を時間軸に沿って輝度変化ないしカラー化して表示するようにする。同図は模式的に示したものであり、図において、補助線は壁厚の等厚線に相当し、各補助線に囲まれる領域r1は低レベルの壁厚変化域、領域r2は中レベルの壁厚変化域、r3は高レベルの壁厚変化域を示し、それらの各領域は輝度(濃淡)又は色を連続的な壁厚変化に応じて連続的に変えて表示される。同図から明らかなように、心筋の真中付近の壁厚変化が他の域に比べて遅れている一方、心筋両端部の壁厚変化が激しいことが一目でわかる。例えば、壁厚変化が小さいときは輝度を最小値(例えば、黒)とし、壁厚変化が激しくなるにつれて輝度を増加するようにする。この運動表示画像は、運動算出部6の計測結果に基づいて、画像生成記憶部5によって作成されて記憶される。
【0035】
図8(b)の表示に代えて、図9〜図11に示す運動表示画像を形成することができる。すなわち、図9は、心筋の模擬画像であるシェーマを表示部に表示し、その画像上に図8(b)のような運動表示画像を重ねて表示するようにすることができる。この例によれば、心筋の形状に対応する部位の壁厚変化が表示されるから、心筋の運動状態を極めて認識しやすい。特に、4状態A〜Dにおける壁厚変化の激しい領域r3の変化を観察して、例えば負荷II状態に対応する状態Cになると2ヶ所に分かれていたr3がつながって来ることが観察できる。これは、例えば、激しい運動をすると血流が改善されて心筋が動き出すことを意味する。したがって、この部分の心筋は回復する余地があることを診断できる。なお、計測部位は心筋の断層像に対応させ、運動量の時間変化は動画像的に表示するようにする。
【0036】
また、図10は、図8(b)と同様の運動表示画像の4つの状態A〜Dを並べて表示した例であり、横軸は時間を、縦軸は計測部位に対応する。また、図11は、特定の心筋の位置における壁厚変化をグラフにして4つの状態A〜Dを並べて示した例であり、横軸は時間を、縦軸は壁厚変化を表している。図11の場合は、操作卓3を操作して、表示領域30の計測部位に指定線を入力することにより、その指定線上の計測部位における壁厚変化がグラフ化して表示される。
【0037】
以上の図8(b)〜図11は、心筋の壁厚変化の運動表示画像例であるが、壁厚の変化率又は壁厚の変化速度を同様に画像化して表示することもできる。例えば、壁厚の変化率を画像化する場合は、壁厚が厚く変化する方向を「赤」とし、薄くなる方向を「青」とする壁厚の変化率の運動表示画像を生成して表示することができる。また、壁厚の変化速度は、運動算出部6の計測結果に基づいて壁厚の変化速度のデータを求め、これに基づいて運動表示画像を生成して、表示領域に表示することができる。
【0038】
以上説明したように、上述の実施形態によれば、次のような効果が得られる。まず、心臓の手術後に心臓の動きが片側だけに偏る場合や、心筋壁の壁厚の変化が小さい部位の有無を観察して血の巡りが悪い部位を検査する場合、追跡点であるマークを心筋の動きに追従させて移動表示させても、心壁全体の揺れに惑わされて壁厚の変化の微妙な差異が観察しにくくなる。この点、上記の実施の形態によれば、マークの位置に対応させた心筋の各計測部位の壁厚の変化等の運動量を、輝度や色などの違いによって表示するようにしたから、一目で心臓の動きが偏る症状や、壁厚の変化が小さい部位である血の巡りが悪い部位を検知することができる。
【0039】
特に、心臓の機能をスコアリングによって評価する場合、一般に心臓の心壁に沿って心壁を6分割し、各分割領域ごとに壁厚の変化の度合いを5段階に分け、普通の動きの場合は「5」、動きがない場合は「1」、その中間は動きの度合いに応じて「2」〜「4」の点数をつける。さらに、通常状態Aと、患者に運動をさせて心臓に負荷をかけた負荷I状態Bと、さらに大きな負荷をかけた負荷II状態Cと、負荷を取り除いて回復した状態Dの4つの態様についてスコアリングを実施して、心臓を診断することが行なわれている。この場合、観者の主観的な判断で動きの度合いを決めて点数をつけると、スコアリング結果がバラツクという問題がある。この点、本実施の形態によれば、壁厚の変化、つまり心筋の動きを輝度の違い又は色の違いによって定量的に判別できるから、観者によるスコアリング結果のバラツキを低減して信頼性の高い評価を行うことができ、また経験が少ない観者でも同一の評価結果を得ることができる。特に4状態を対比して観察できるから、一層信頼性の高い評価を行うことができる。
【0040】
また、心筋の運動を定量的に認識できるから、虚血性心疾患において例えば虚血部位を特定したり、虚血の程度がわかるので、冠動脈再生術などの治療法選択および治療部位を特定する際の指標にできる。
【0041】
また、上述の実施形態において、心筋の動きを木目細かく観察するために、追跡部位のマークである追跡点23を操作卓3から入力設定する数が多くなり、設定作業が煩雑である。そこで、心筋壁に沿って追跡点23を設定する場合、操作者の判断で組織形状が緩やかに変化するような部位については適宜間隔を空けて、組織形状が大きく変化するような部位については間隔を狭めて設定するようにしてもよい。この場合、制御手段10により追跡点23を密な間隔に自動的に補完設定するようにすることが好ましい。
【0042】
また、上述の実施形態は、オフラインで行なう例について説明したが、画像追跡処理に係る速度を向上すれば、オンラインあるいはリアルタイムの動画像にも適用できる。また、2次元の断層像を例に説明したが、3次元断層像にも適用できることはいうまでもない。
【0043】
【発明の効果】
以上述べたように、本発明によれば、定量的に計測した心室壁厚の表示法を改善したことから、的確な診断を行うことができる。
【図面の簡単な説明】
【図1】図1は、本発明の生体組織の動き追跡表示方法の一実施形態の処理手順を示す図である。
【図2】図2は、図1の生体組織の動き追跡表示方法を適用してなる画像診断装置のブロック構成図である。
【図3】図3は、本発明の生体組織の動き追跡を、心臓の断層像に適用して説明するための図である。
【図4】図4は、本発明に係る画像追跡処理法の一実施形態を説明する図であり、(a)は切出し画像の一例を、(b)は検索領域の一例を示す図である。
【図5】図5は、相関法による画像追跡処理を、具体例を用いて説明する図である。
【図6】図6は、心壁を挟んで設定された複数対のマークの距離と、その距離の変化を計測してグラフにして表示する例である。
【図7】図7は、心臓のスコアリングのための運動負荷が異なる4つの心筋断層像を並べて表示した一例図である。
【図8】図8は、心筋の壁厚変化を示すグラフ、及び壁厚変化を時間軸に沿って輝度変化ないしカラー化して表示する運動表示画像の一例を示す図である。
【図9】図9は、心筋の模擬画像であるシェーマを表示部に表示し、その画像上に本発明の4状態の運動表示画像を重ねて表示する例を示した図である。
【図10】図10は、壁厚変化を時間軸に沿って輝度変化ないしカラー化して表示する運動表示画像を4状態について並べて表示した例を示した図である。
【図11】図11は、図10の特定の計測部位における壁厚変化をグラフ化して表示した一例図である。
【符号の説明】
1 画像記憶部
2 表示部
3 操作卓
4 自動追跡部
5 画像生成記憶部
6 運動算出部
7 信号伝送路
8 診断画像撮像装置
10 制御手段
11 表示制御手段
12 切出し画像設定手段
13 切出し画像追跡手段
14 移動量演算手段
15 移動追跡手段
Claims (8)
- 負荷状態が異なる心臓の断層を撮像してなる少なくとも2つの動画像が記憶された画像記憶部から一フレーム画像を読み出して表示部に表示する第1ステップと、該表示された前記一フレーム画像上において心筋を挟んで対向する少なくとも一対のマークを設定する第2ステップと、前記各マークが設定された部位の心筋の動きを前記動画像上で追跡する第3ステップと、前記心筋の動きに合わせて前記各マークの位置を変更し、変更後の前記一対のマーク間の距離により前記心筋の壁厚を求める第4ステップと、前記壁厚と前記壁厚の変化率と前記壁厚の変化速度の少なくとも1つの運動量を求める第5ステップと、前記各対のマーク位置における前記心筋の運動量を、該運動量に応じて輝度と色の少なくとも一方を変えた運動表示画像を前記表示部に表示する第6ステップとからなり、前記負荷状態が異なる心臓の断層像について前記第1ステップから第6ステップを実行する心臓の運動表示方法。
- 前記各負荷状態の前記運動表示画像を表示画面に並べて表示することを特徴とする請求項1に記載の心臓の運動表示方法。
- 前記運動表示画像は、矩形領域の一方の軸を前記心筋に沿って配置された複数対のマークの設定位置に対応させ、他方の軸を時間軸に対応させてなることを特徴とする請求項1又は2に記載の心臓の運動表示方法。
- 前記矩形領域の時間と複数対のマークの設定位置の少なくとも一方を指定する線を入力設定することにより、該指定線に沿った前記運動量の変化を表すグラフを前記表示部に表示する第7ステップを有することを特徴とする請求項3に記載の心臓の運動表示方法。
- 前記運動表示画像は、前記各対のマーク位置における前記動画像の心筋に対応する領域の輝度と色の少なくとも一方を当該位置における前記運動量に応じて変えた画像であることを特徴とする請求項1又は2に記載の心臓の運動表示方法。
- 前記運動表示画像は、前記心筋の内壁と外壁を模擬してなる模式図形に、該模式図形における前記各対のマーク位置に対応する領域の輝度と色の少なくとも一方を当該位置における前記運動量に応じて変えた画像を重ねて生成されるものであることを特徴とする請求項1又は2に記載の心臓の運動表示方法。
- 負荷状態が異なる心臓の断層を撮像して得られる少なくとも2つの動画像データが記憶される画像記憶手段と、前記心臓の断層像を表示する表示部と、前記画像記憶手段に記憶された前記動画像データを読み出して動画像又は静止画像を前記表示部に表示させる表示制御手段と、前記表示部に表示された前記静止画像上の心筋を挟んで対向する少なくとも一対のマークを設定する操作部と、前記表示部に表示された静止画像に前記各マークを含む切出し画像をそれぞれ設定し、前記表示部に表示される前記動画像に基づいて前記切出し画像の移動先を追跡する切り出し画像追跡手段と、該切出し画像の移動先に基づいて前記各マークの移動量を求める移動量演算手段と、該手段により求めた各マークの移動量に基づいて前記一対のマーク間の距離を求めて前記心筋の壁厚と該壁厚の変化率と該壁厚の変化速度の少なくとも1つの運動量を求める運動算出手段と、該手段により求めた前記運動量に応じて輝度と色の少なくとも一方の画素値を変えて前記表示部の表示領域に表示する運動表示画像を生成する画像生成手段と、前記負荷状態が異なる心臓の各断層像について前記運動表示画像を生成させて前記表示部に表示させる制御手段とを備えてなる画像診断装置。
- 前記表示制御手段は、前記負荷状態が異なる心臓の前記運動表示画像を表示画面に並べて表示させることを特徴とする請求項7に記載の画像診断装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003113390A JP4359749B2 (ja) | 2003-04-17 | 2003-04-17 | 生体組織の運動表示方法及び画像診断装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2003113390A JP4359749B2 (ja) | 2003-04-17 | 2003-04-17 | 生体組織の運動表示方法及び画像診断装置 |
Publications (3)
Publication Number | Publication Date |
---|---|
JP2004313551A true JP2004313551A (ja) | 2004-11-11 |
JP2004313551A5 JP2004313551A5 (ja) | 2006-06-15 |
JP4359749B2 JP4359749B2 (ja) | 2009-11-04 |
Family
ID=33473349
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2003113390A Expired - Lifetime JP4359749B2 (ja) | 2003-04-17 | 2003-04-17 | 生体組織の運動表示方法及び画像診断装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4359749B2 (ja) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007013321A1 (ja) * | 2005-07-27 | 2007-02-01 | Hitachi Medical Corporation | 画像診断装置及び計測点設定方法並びにプログラム |
JP2009028362A (ja) * | 2007-07-27 | 2009-02-12 | Toshiba Corp | 医用画像処理装置及び医用画像診断装置 |
WO2009078297A1 (ja) * | 2007-12-19 | 2009-06-25 | Konica Minolta Medical & Graphic, Inc. | 動態画像処理システム |
JP2009148336A (ja) * | 2007-12-19 | 2009-07-09 | Konica Minolta Medical & Graphic Inc | 動態画像診断支援システム |
JP2009160221A (ja) * | 2008-01-07 | 2009-07-23 | Toshiba Corp | X線コンピュータ断層撮影装置及び3次元画像処理装置 |
JP2009160314A (ja) * | 2008-01-09 | 2009-07-23 | Canon Inc | 画像処理装置、画像処理方法、及び、コンピュータプログラム |
JP2009538172A (ja) * | 2006-05-25 | 2009-11-05 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 心室壁濃厚部の定量化及び表示 |
JP2010075403A (ja) * | 2008-09-25 | 2010-04-08 | Canon Inc | 情報処理装置およびその制御方法、データ処理システム |
JP2010101811A (ja) * | 2008-10-24 | 2010-05-06 | Fujifilm Ri Pharma Co Ltd | 心電図同期画像処理装置、方法及びコンピュータプログラム |
WO2012056662A1 (ja) * | 2010-10-25 | 2012-05-03 | 株式会社 東芝 | 医用画像処理装置、医用画像撮影装置及び医用画像処理プログラム |
JP2012090722A (ja) * | 2010-10-26 | 2012-05-17 | Ge Medical Systems Global Technology Co Llc | 画像処理装置およびプログラム並びに画像診断装置 |
JP2013059695A (ja) * | 2013-01-10 | 2013-04-04 | Toshiba Corp | X線コンピュータ断層撮影装置、画像処理装置及びプログラム |
US8727992B2 (en) | 2006-06-02 | 2014-05-20 | Kabushiki Kaisha Toshiba | Ultrasonic doppler diagnostic apparatus, and method of controlling ultrasonic doppler diagnostic apparatus |
WO2016039100A1 (ja) * | 2014-09-11 | 2016-03-17 | 日立アロカメディカル株式会社 | 超音波診断装置 |
-
2003
- 2003-04-17 JP JP2003113390A patent/JP4359749B2/ja not_active Expired - Lifetime
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007013321A1 (ja) * | 2005-07-27 | 2007-02-01 | Hitachi Medical Corporation | 画像診断装置及び計測点設定方法並びにプログラム |
JP2009538172A (ja) * | 2006-05-25 | 2009-11-05 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 心室壁濃厚部の定量化及び表示 |
US8727992B2 (en) | 2006-06-02 | 2014-05-20 | Kabushiki Kaisha Toshiba | Ultrasonic doppler diagnostic apparatus, and method of controlling ultrasonic doppler diagnostic apparatus |
JP2009028362A (ja) * | 2007-07-27 | 2009-02-12 | Toshiba Corp | 医用画像処理装置及び医用画像診断装置 |
US8483456B2 (en) | 2007-12-19 | 2013-07-09 | Konica Minolta Medical & Graphic, Inc. | Dynamic image processing system |
WO2009078297A1 (ja) * | 2007-12-19 | 2009-06-25 | Konica Minolta Medical & Graphic, Inc. | 動態画像処理システム |
JP2009148336A (ja) * | 2007-12-19 | 2009-07-09 | Konica Minolta Medical & Graphic Inc | 動態画像診断支援システム |
US9125621B2 (en) | 2007-12-19 | 2015-09-08 | Konica Minolta, Inc. | Dynamic image processing system |
JP2009160221A (ja) * | 2008-01-07 | 2009-07-23 | Toshiba Corp | X線コンピュータ断層撮影装置及び3次元画像処理装置 |
JP2009160314A (ja) * | 2008-01-09 | 2009-07-23 | Canon Inc | 画像処理装置、画像処理方法、及び、コンピュータプログラム |
US9147098B2 (en) | 2008-01-09 | 2015-09-29 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method, and computer program storage medium |
JP2010075403A (ja) * | 2008-09-25 | 2010-04-08 | Canon Inc | 情報処理装置およびその制御方法、データ処理システム |
JP2010101811A (ja) * | 2008-10-24 | 2010-05-06 | Fujifilm Ri Pharma Co Ltd | 心電図同期画像処理装置、方法及びコンピュータプログラム |
CN102665566A (zh) * | 2010-10-25 | 2012-09-12 | 株式会社东芝 | 医用图像处理装置、医用图像拍摄装置和医用图像处理程序 |
WO2012056662A1 (ja) * | 2010-10-25 | 2012-05-03 | 株式会社 東芝 | 医用画像処理装置、医用画像撮影装置及び医用画像処理プログラム |
CN102665566B (zh) * | 2010-10-25 | 2015-10-21 | 株式会社东芝 | 医用图像处理装置、医用图像拍摄装置和医用图像处理方法 |
JP2012090722A (ja) * | 2010-10-26 | 2012-05-17 | Ge Medical Systems Global Technology Co Llc | 画像処理装置およびプログラム並びに画像診断装置 |
JP2013059695A (ja) * | 2013-01-10 | 2013-04-04 | Toshiba Corp | X線コンピュータ断層撮影装置、画像処理装置及びプログラム |
WO2016039100A1 (ja) * | 2014-09-11 | 2016-03-17 | 日立アロカメディカル株式会社 | 超音波診断装置 |
JP2016055040A (ja) * | 2014-09-11 | 2016-04-21 | 日立アロカメディカル株式会社 | 超音波診断装置 |
Also Published As
Publication number | Publication date |
---|---|
JP4359749B2 (ja) | 2009-11-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1543773B1 (en) | Biological tissue motion trace method and image diagnosis device using the trace method | |
JP4389081B2 (ja) | 関心領域の追従制御方法、その方法を用いた画像診断装置及び関心領域の追従制御プラグラム | |
AU2004273587B2 (en) | Method and device for visually supporting an electrophysiology catheter application in the heart | |
CN100525714C (zh) | 图像处理设备和超声诊断设备 | |
JP4596759B2 (ja) | 画像診断装置、及び画像診断装置の作動方法 | |
JP5581511B2 (ja) | 中心点軌跡を図表化するためのコンピュータプログラムが記録されたコンピュータ可読記憶媒体 | |
JP4359749B2 (ja) | 生体組織の運動表示方法及び画像診断装置 | |
JP5215036B2 (ja) | 医用画像処理装置、及び医用画像処理プログラム | |
JP2012205899A (ja) | 3次元的モデルを利用した身体臓器の映像生成方法及び装置並びにコンピュータ読み取り可能な記録媒体 | |
US20090220050A1 (en) | Method for Determining and Displaying at Least One Piece of Information on a Target Volume | |
JP4321121B2 (ja) | 診断画像における生体組織の動きの追跡方法及びその方法を用いた画像診断装置 | |
JP4503238B2 (ja) | 生体組織の運動表示方法及び画像診断装置 | |
JP3844799B2 (ja) | 超音波診断装置、超音波診断用のデータ処理方法、及びプログラム | |
US8224047B2 (en) | System and method for measuring left ventricular torsion | |
CN110477952A (zh) | 超声波诊断装置、医用图像诊断装置及存储介质 | |
JP4847684B2 (ja) | 画像診断装置 | |
JP2005087594A (ja) | 画像診断装置 | |
US20120008833A1 (en) | System and method for center curve displacement mapping | |
JP2008104695A (ja) | 超音波診断装置、画像処理装置、及び画像処理プログラム | |
JP5508035B2 (ja) | 画像診断装置 | |
JP2010505494A (ja) | 医用イメージングシステム | |
JP2004283583A (ja) | 画像形成医療検査装置の作動方法 | |
JP4352212B2 (ja) | 生体組織の動き追跡表示方法及び画像診断装置 | |
WO2022064445A1 (en) | Automated method for identifying and indicating to an operator pathological risk regions in at least one part of a patient's cardiovascular system by means of the reconstruction of an augmented reality of morphology and hemodynamics | |
JP2005160688A (ja) | 画像診断装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060413 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060413 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20081215 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090106 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090304 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090407 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090605 |
|
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: 20090714 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20090728 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 Ref document number: 4359749 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120821 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120821 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130821 Year of fee payment: 4 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
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 |
|
EXPY | Cancellation because of completion of term |