JP5437820B2 - 超音波診断装置、超音波画像処理方法 - Google Patents

超音波診断装置、超音波画像処理方法 Download PDF

Info

Publication number
JP5437820B2
JP5437820B2 JP2009554356A JP2009554356A JP5437820B2 JP 5437820 B2 JP5437820 B2 JP 5437820B2 JP 2009554356 A JP2009554356 A JP 2009554356A JP 2009554356 A JP2009554356 A JP 2009554356A JP 5437820 B2 JP5437820 B2 JP 5437820B2
Authority
JP
Japan
Prior art keywords
graph
image
ultrasonic
unit
compression
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
JP2009554356A
Other languages
English (en)
Other versions
JPWO2009104657A1 (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.)
Hitachi Healthcare Manufacturing Ltd
Original Assignee
Hitachi Medical Corp
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 Hitachi Medical Corp filed Critical Hitachi Medical Corp
Priority to JP2009554356A priority Critical patent/JP5437820B2/ja
Publication of JPWO2009104657A1 publication Critical patent/JPWO2009104657A1/ja
Application granted granted Critical
Publication of JP5437820B2 publication Critical patent/JP5437820B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0833Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/485Diagnostic techniques involving measuring strain or elastic properties
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023Details of receivers
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023Details of receivers
    • G01S7/52036Details of receivers using analysis of echo signal for target characterisation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023Details of receivers
    • G01S7/52036Details of receivers using analysis of echo signal for target characterisation
    • G01S7/52042Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52053Display arrangements
    • G01S7/52057Cathode ray tube displays
    • G01S7/52071Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52053Display arrangements
    • G01S7/52057Cathode ray tube displays
    • G01S7/52074Composite displays, e.g. split-screen displays; Combination of multiple images or of images and alphanumeric tabular information

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Medical Informatics (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Radiology & Medical Imaging (AREA)
  • Pathology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biophysics (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Description

本発明は、超音波診断装置に関し、具体的には、被検体の断層部位の組織の硬さ又は軟らかさを表す弾性画像を撮像して表示する超音波診断装置、超音波画像処理方法及び超音波画像処理プログラムに関する。
超音波診断装置は、超音波探触子により被検体内部に超音波を送信するとともに、被検体内部から生体組織の構造に応じた超音波の反射エコー信号を受信し、例えば超音波断層像等の画像を構成して診断用に表示する。
近年、手動又は機械的な方法により超音波探触子で被検体を圧迫して超音波受信信号を計測し、計測時間が異なる2つの超音波受信信号のフレームデータに基づいて圧迫により生じた生体各部の変位を求め、その変位データに基づいて生体組織の弾性を表す弾性画像を生成することが開示されている(例えば、特許文献1)。
ところで、被検体への圧迫操作は簡便性の観点から超音波探触子を用いることが多い。つまり、超音波探触子を被検体の体表に押し当ててある程度の初期圧迫を加え、この初期圧迫状態から、微小な加圧と減圧を繰り返し与えながら弾性画像が生成される。ところが、この手法では、操作者の手加減によって加圧状態が大きく変動するので、必ずしも適正な弾性画像が得られない場合がある。
そこで、従来の超音波診断装置では、計測時刻の異なる2つの断層像データに基づいて各部の変位を計測し、変位データに基づいて断層像の歪み率分布を求めて、歪み率分布が適正範囲内か否かにより加圧の適否を判定することが開示されている(例えば、特許文献2)。
また、超音波探触子の被検体圧迫面に圧力センサを設けて、圧力センサの計測値の推移をグラフにして弾性画像とともに表示することも開示されている(例えば、特許文献3)。
特開2000-060853号公報 特開2004-351062号公報 国際公開WO2005-120358号公報
しかしながら、特許文献2,3に記載された技術は、被検体の断層部位の関心組織に加えられている圧迫状態をより精度よく操作者に提示することに関して考慮されていない。
すなわち、特許文献2に記載の技術は、断層部位の領域内の組織の歪みの分布に基づいて圧迫の評価を行うものであり、関心組織以外の周辺組織の歪みデータも含まれるため、関心組織に加えられている圧迫状態を精度よく求めることができない場合がある。
また、特許文献3に記載されているように超音波探触子の圧迫面の圧力センサ値を用いる場合は、センサ表面の圧迫強度が計測されるため、被検体内の関心組織に加えられている圧迫状態を直接的に反映することができない場合がある。
そこで、本発明は、2点の計測点間の距離の推移のグラフ化することにより、関心組織に加えられている圧迫状態の推移をより精度よく操作者に対して提示して適切な圧迫操作をサポートすることを課題とする。
上記課題を解決するため、本発明の超音波診断装置は、被検体に圧迫を加えながら被検体との間で超音波を送受信する超音波探触子と、該超音波探触子で計測された反射エコー信号に基づいて前記被検体の断層部位のRF信号フレームデータを生成する整相加算部と、前記RF信号フレームデータに基づいて前記被検体の断層部位の超音波画像を生成する超音波画像生成部と、前記RF信号フレームデータ又は前記超音波画像に計測対象を設定する計測対象設定部と、設定された計測対象の位置変化を追跡する計測対象位置追跡部と、前記計測対象の位置変化に基づいて前記計測対象に加えられている圧迫状態の推移を表すグラフを生成する変位画像生成部と、前記生成されたグラフを表示する表示器と、を備え、前記計測対象設定部は、前記RF信号フレームデータ又は前記超音波画像の関心組織とその関心組織と異なる他の領域との境界に、前記圧迫が加えられる方向に沿った2点の計測点を設定し、前記変位画像生成部は、設定された2点の計測点間の距離変化の推移のグラフを生成し、前記超音波画像生成部によって生成された超音波画像と前記変位画像生成部によって生成されたグラフとを関連づけて前記表示器に表示される画像処理部を備えたことを特徴とする。
以上説明したように、本発明よれば、2点の計測点間の距離の推移のグラフ化することにより、関心組織に加えられている圧迫状態の推移をより精度よく操作者に対して提示して適切な圧迫操作をサポートすることが可能になる。

本実施形態の超音波診断装置の構成を示すブロック図。 超音波診断装置の特徴部の第1実施例を説明する図。 第1実施例の処理のフローチャート。 超音波診断装置の特徴部の第2実施例を説明する図。 第2実施例の処理のフローチャート。 超音波診断装置の特徴部の第3実施例を説明する図。 第3実施例の処理のフローチャート。 超音波診断装置の特徴部の第4実施例を説明する図。 第4実施例の処理のフローチャート。 初期圧迫が適切か否かを説明するための模式図。
符号の説明
1 超音波診断装置、10 被検体、12 超音波探触子、18 整相加算部、20 断層画像構成部、26 画像表示器、30 変位計測部、32 弾性情報演算部、34 弾性画像構成部、38 変位演算部、40 変位画像構成部、42 操作部、44 腫瘍、46 計測点、48,53,55 グラフ、50 振幅閾値ライン、52 周期ライン、54 参照体、60 閾値
発明を実施するための形態
以下、本発明を適用してなる超音波診断装置の実施形態を説明する。なお、以下の説明では、同一機能部品については同一符号を付して重複説明を省略する。
図1は、本実施形態の超音波診断装置の構成を示すブロック図である。図1に示すように、超音波診断装置1には、被検体10に当接させて用いる超音波探触子12と、超音波探触子12を介して被検体10に時間間隔をおいて超音波を繰り返し送信する送信部14と、被検体10から発生する時系列の反射エコー信号を受信する受信部16と、送信部14と受信部16を制御する送受信制御部17と、受信部16で受信された反射エコーを整相加算する整相加算部18と、が備えられている。
また、整相加算部18からのRF信号フレームデータに基づいて被検体の濃淡断層画像例えば白黒断層画像を構成する断層画像構成部20と、断層画像構成部20の出力信号を画像表示器26の表示に合うように変換する白黒スキャンコンバータ22とが備えられている。
また、整相加算部18から出力されるRF信号フレームデータを記憶し、少なくとも2枚のフレームデータを選択するRFフレームデータ選択部28と、被検体10の生体組織の変位を計測する変位計測部30と、変位計測部30で計測された変位情報から歪み又は弾性率を求める弾性情報演算部32と、弾性情報演算部32で演算した歪み又は弾性率からカラー弾性画像を構成する弾性画像構成部34と、弾性画像構成部34の出力信号を画像表示器26の表示に合うように変換するカラースキャンコンバータ36とが備えられている。
また、白黒断層画像とカラー弾性画像を重ね合わせたり、並列に表示させたり、切替を行う切替加算部24と、合成された合成画像を表示する画像表示器26とが備えられている。なお、断層画像構成部20と弾性画像構成部34などにより、RF信号フレームデータに基づいて被検体の断層部位の超音波画像を生成する超音波画像生成手段が構成される。
ここで、超音波診断装置1の各構成部について詳細に説明する。超音波探触子12は、複数の振動子を配設して形成されており、被検体10に振動子を介して超音波を送受信する機能を有している。送信部14は、超音波探触子12を駆動して超音波を発生させるための送波パルスを生成するとともに、送信される超音波の収束点をある深さに設定する機能を有している。また、受信部16は、超音波探触子12で受信した反射エコー信号について所定のゲインで増幅してRF信号すなわち受波信号を生成するものである。整相加算部18は、受信部16で増幅されたRF信号を入力して位相制御し、一点又は複数の収束点に対し超音波ビームを形成してRF信号フレームデータを生成するものである。
断層画像構成部20は、整相加算部18からのRF信号フレームデータを入力してゲイン補正、ログ圧縮、検波、輪郭強調、フィルタ処理等の信号処理を行い、断層画像データを得るものである。また、白黒スキャンコンバータ22は、断層画像構成部20からの断層画像データをデジタル信号に変換するA/D変換器と、変換された複数の断層画像データを時系列に記憶するフレームメモリと、制御コントローラを含んで構成されている。この白黒スキャンコンバータ22は、フレームメモリに格納された被検体内の断層フレームデータを1画像として取得し、取得された断像フレームデータをテレビ同期で読み出すものである。
RFフレームデータ選択部28は、整相加算部18からの複数のRF信号フレームデータを格納し、格納されたRF信号フレームデータ群から取得時刻の異なる一対すなわち2つのRF信号フレームデータを選択する。例えば、整相加算部18から時系列すなわち画像のフレームレートに基づいて生成されるRF信号フレームデータをRFフレームデータ選択部28に順次記憶し、記憶されたRF信号フレームデータ(N)を第1のデータとして選択すると同時に、時間的に過去に記憶されたRF信号フレームデータ群(N-1、N-2、N-3…N―M)の中から1つのRF信号フレームデータ(X)を選択する。なお、ここでN、M、XはRF信号フレームデータに付されたインデックス番号であり、自然数とする。
そして、変位計測部30は、選択された一対のデータすなわちRF信号フレームデータ(N)及びRF信号フレームデータ(X)から1次元或いは2次元相関処理を行って、断層画像の各点に対応する生体組織における変位や移動ベクトルすなわち変位の方向と大きさに関する1次元又は2次元変位分布を求める。ここで、移動ベクトルの検出にはブロックマッチング法を用いる。ブロックマッチング法とは、画像を例えばN×N画素からなるブロックに分け、関心領域内のブロックに着目し、着目しているブロックに最も近似しているブロックを前のフレームから探し、これを参照して予測符号化すなわち差分により標本値を決定する処理を行う。
弾性情報演算部32は、変位計測部30から出力される計測値、例えば変位や移動ベクトルから断層画像上の各点に対応する生体組織の歪みや弾性率を演算し、その歪みや弾性率に基づいて弾性画像信号すなわち弾性フレームデータを生成するものである。
このとき、歪みのデータは、生体組織の移動量例えば変位を空間微分することによって算出される。また、弾性率のデータは、超音波探触子の圧迫面に設けた図示していない圧力センサから出力される圧力の変化を歪みの変化で除することによって計算される。例えば、変位計測部30により計測された変位をL(X)、圧力センサにより計測された圧力をP(X)とすると、歪みΔS(X)は、L(X)を空間微分することによって算出することができるから、
ΔS(X)=ΔL(X)/ΔXという式を用いて求められる。
また、弾性率データのヤング率Ym(X)は、Ym=(ΔP(X))/ΔS(X)
という式によって算出される。このヤング率Ymから断層画像の各点に相当する生体組織の弾性率が求められるので、2次元の弾性画像データを連続的に得ることができる。なお、ヤング率とは、物体に加えられた単純引張り応力と、引張りに平行に生じるひずみに対する比である。
弾性画像構成部34は、フレームメモリと画像処理部とを含んで構成されており、弾性情報演算部32から時系列に出力される弾性フレームデータをフレームメモリに確保し、確保されたフレームデータに対し画像処理を行うものである。
カラースキャンコンバータ36は、弾性画像構成部34からの弾性フレームデータに色相情報を付与する機能を有したものである。つまり、弾性フレームデータに基づいて光の3原色すなわち赤(R)、緑(G)、青(B)に変換するものである。例えば、歪みが大きい弾性データを赤色コードに変換すると同時に、歪みが小さい弾性データを青色コードに変換する。
切替加算部24は、フレームメモリと、画像処理部と、画像選択部とを備えて構成されている。ここで、フレームメモリは、白黒スキャンコンバータ22からの断層画像データ、カラースキャンコンバータ36からの弾性画像データ、及び後述する変位画像構成部40からのグラフなどの画像を格納するものである。また、画像処理部は、フレームメモリに確保された断層画像データ、弾性画像データ、及び関心組織に加えられている圧迫状態の推移を表すグラフデータの合成割合を変更して合成するものである。合成画像の各画素の輝度情報及び色相情報は、白黒断層画像とカラー弾性画像の各情報を合成割合で加算したものとなる。さらに、画像選択部は、フレームメモリ内の断層画像データと弾性画像データ及び画像処理部の合成画像データのうちから画像表示器26に表示する画像を選択するものである。
以下、本実施形態の超音波診断装置の特徴部分について説明する。本実施形態の超音波診断装置1は、図1に示すように、変位演算部38と、変位画像構成部40と、操作部42とを備えて構成されている。
変位演算部38は、白黒スキャンコンバータ22から出力される断層画像データの設定された計測点の位置変化を追跡して変位を演算するものである。計測点の設定は、画像表示器26に表示された断層画像上の例えば腫瘍などの関心組織に操作者がキーボード、マウス、タッチセンサなどの入力インターフェースを介して設定することもできるし、超音波診断装置により自動で腫瘍などの関心組織を検出して、その関心組織に設定することもできる。つまり、操作者による入力設定、或いは超音波診断装置による自動設定などにより、RF信号フレームデータ又は超音波画像の関心組織に計測点を設定する手段が構成される。
例えば画像表示器26上で計測点が設定されたら、計測点の位置を示すデータが操作部42を介して変位演算部38に入力され、変位演算部38は入力されたデータに基づいて断層画像データの計測点の位置変化を追跡する。変位演算部38は、設定された計測点の位置変化を追跡する手段として機能する。計測点の追跡は、例えばトラッキングなどの手法を用いることができる。また、自動で腫瘍などの関心組織を検出するのは、例えば腫瘍における反射エコーの強度は小さいことから、RFフレーム信号を強度により2値化して、ローエコー部分を腫瘍とみなすなどの手法を採用することができる。
変位画像構成部40は、変位演算部38から出力される計測点の位置変化に基づいて腫瘍などの関心組織に加えられている圧迫状態の推移を表すグラフを生成するものである。
なお、超音波診断装置に、上述のような入力インターフェース、画像表示器26のような出力インターフェース、ソフトウェアプログラムを格納するメモリ、及びソフトウェアプログラムを読み出して実行可能な演算手段を設けることにより、RF信号フレームデータ又は超音波画像の関心組織に計測点を設定する手段、変位演算部38、変位画像構成部40などによる機能を、ソフトウェアプログラムにより構成して実行させることができる。
以下、計測点の設定態様、生成されるグラフの態様、及び操作者の適切な圧迫操作をサポートするための画像の生成態様などについて、具体的な実施例を用いて説明する。
図2は、本実施形態の超音波診断装置の特徴部の第1実施例を説明する図であり、計測点の設定例、及びこの場合に生成されるグラフを模擬的に示している。図3は、第1実施例の処理のフローチャートである。
図2(a)に示すように、計測点46は、腫瘍44の境界部の圧迫が加えられる方向に沿って2点設定されている。変位演算部38は、この2つの計測点をトラッキングなどによりフレームごとに追跡するとともに、2点の計測点間の距離を全フレームメモリのデータに対して演算し、その結果を逐次変位画像構成部40に出力する。変位画像構成部40は、2点の計測点間の距離変化の推移をグラフ48として生成する。
変位画像構成部40から出力されたグラフ48は、切替加算部24により加算合成され、図2(b)に示すように、断層画像と弾性画像が重畳された画像とともに、画像表示器26に表示される。
次に、図3のフローチャートを用いて、第1実施例の動作を説明する。まず、関心組織である腫瘍の境界に、白黒スキャンコンバータ22から得られたデータに対して2値化処理を行うことにより自動的に、もしくは操作部42を介して手動で、2点の計測点46を圧迫方向に沿って設定する(S11)。続いて、変位演算部38にて設定された計測点間の距離Aを検出する(S12)。さらに、圧迫により変化する計測点の位置をトラッキングにより検出し、全てのフレームメモリにおける計測点間の距離(A)iを検出する(S13)。続いて、変位画像構成部40にて(A)iから組織の変位を表すグラフ48を生成し、切替加算部24を介して画像表示器26に表示する(S14)。
本実施例によれば、被検体の断層部位の関心組織の圧迫が加えられる方向に沿った2点の計測点を設定し、その計測点間の距離変化の推移のグラフ、言い換えれば圧迫状態の推移を表すグラフを生成して表示しているので、このグラフには、被検体の断層部位の関心組織に加えられている圧迫状態の推移が直接的に精度よく反映される。その結果、操作者は、自らが行っている操作により関心組織に加えられている圧迫状態を視覚的に的確に把握することができるので、これに基づいて圧迫操作が適切か否かを判断することができる。
計測点を1点のみ設定して、この計測点を追跡した位置変化をグラフ化することも可能であるが、本実施例のように、2点の計測点間の距離の推移のグラフ化することにより、関心組織に加えられている圧迫状態の推移をより精度よく操作者に対して提示することができる。
図4は、本実施形態の超音波診断装置の特徴部の第2実施例を説明する図であり、操作者の適切な圧迫操作をサポートするために生成される画像を示している。図5は、第2実施例の処理のフローチャートである。
本実施例は、グラフ48に現れる複数の極大点及び極小点に基づいてグラフの中点位置を求め、グラフ48の中点位置から所定の振幅位置に振幅閾値のラインを重畳して表示するとともに、グラフ48上に、あらかじめ設定された周期のラインを重畳して表示するものである。
さらに、グラフ48の周期のライン及び振幅閾値のラインの少なくとも一方を逸脱した箇所の表示態様を変更して操作者に適切な圧迫操作を促すようにしてものである。
具体的には、図4(a)に示すように、例えば第1実施例のようにして生成されたグラフ48上に、あらかじめ設定された周期ライン52を重畳して表示して、操作者に視覚的に適切な周期を提示することができる。
一方、グラフ48に対して、計測点間の距離の各極大点を検出して各極大点間の周期Cを検出し、以下の式を用いて全体の周期(C)aveを算出することもできる。
(C)ave=(Σ(各極大点間の周期C))/(極大点の数-1)
このように、操作者の操作による(C)aveを演算で求めることにより、あらかじめ設定された推奨される周期との相違の程度を数値として検出することができ、これを操作者に提示することができる。
また、グラフ48において、計測点間の距離の各極大点、及び各極小点を検出して計測点間の距離の各中点位置を算出するとともに、以下の式を用いて圧迫全体の計測点間の距離の中間値(D)aveを算出することにより、振幅の閾値を表示することができる。
(D)ave=(Σ(各中点位置))/(極大点の数)
つまり、このようにして(D)aveが求められたら、この(D)aveを基準に、グラフ48の振幅上下方向に、あらかじめ定められた幅の振幅閾値ライン50を重畳表示することができ、操作者に視覚的に適切な振幅を提示することができる。なお、推奨される周期、及び振幅は、対象組織に応じて適宜変更して保存しておき、対象部位を選択可能にしておいてもよいし、操作者の任意に入力するように構成してもよい。
また、図4(b)に示すように、グラフ48に重畳された周期ライン52と振幅閾値ライン50を用いて、圧迫の適切性の判定を行い、操作者に提示することができる。つまり、振幅閾値ライン50から外れた圧迫、もしくは周期ライン52から外れた圧迫を不適切な圧迫と判断し、グラフ48の該当する箇所の色を変更して表示する。これにより、操作者は、圧迫の適切性を視覚的に容易に判断することが可能となる。
次に、第2実施例の動作について図5を用いて説明する。まず、グラフ48にあらかじめ定められた推奨される周期ライン52を重畳する(S21)。続いて、圧迫全体の計測点間の距離の中間値(D)aveを算出し、中間値から上下方向に、あらかじめ定められた幅の振幅閾値ライン50を重畳する(S22)。そして、振幅閾値ライン50から外れた圧迫、もしくは周期ライン52から外れた圧迫を不適切な圧迫と判断し、組織変位グラフの該当する箇所の色を変更する(S23)。
図6は、本実施形態の超音波診断装置の特徴部の第3実施例を説明する図であり、計測点の設定例、及びこの場合に生成されるグラフを模擬的に示している。図7は、第3実施例の処理のフローチャートである。
図6(a)に示すように、腫瘍44の境界部の圧迫が加えられる方向に沿った2点の計測点46の対が、圧迫が加えられる方向に直交する方向の異なる位置に複数設定されている。また、変位演算部38は、各計測点をトラッキングなどによりフレームごとに追跡するとともに、2点の計測点の対ごとに、対となる計測点間の距離を全フレームメモリのデータに対して演算し、その結果を逐次変位画像構成部40に出力する。変位画像構成部40は、それぞれの2点の計測点の対ごとに、計測点間の距離変化の推移をグラフ53として生成する。
変位画像構成部40から出力されたグラフ53は、切替加算部24により加算合成され、図6(b)に示すように、断層画像と弾性画像が重畳された画像とともに、画像表示器26に表示される。つまり、計測点46の対を図面左側から順に対1,対2,対3とすると、対1〜対3のそれぞれに対する計測点間の距離変化の推移が縦に並べて表示される。
次に、図7のフローチャートを用いて、第1実施例の動作を説明する。まず、関心組織である腫瘍の境界に、白黒スキャンコンバータ22から得られたデータに対して2値化処理を行うことにより自動的に、もしくは操作部42を介して手動で、2点の計測点46の対を複数(例えば対1〜対3)圧迫方向に沿って設定する(S31)。続いて、対1〜対3のそれぞれに対して、変位演算部38にて設定された計測点間の距離Aを検出する(S32)。さらに、圧迫により変化する計測点の位置をトラッキングにより検出し、対1〜対3のそれぞれに対して、全てのフレームメモリにおける計測点間の距離(A)iを検出する(S33)。続いて、対1〜対3のそれぞれに対して、変位画像構成部40にて(A)iから組織の変位を表すグラフ53を生成し、切替加算部24を介して画像表示器26に表示する(S34)。
本実施例によれば、操作者は、これら各グラフを参照することにより、圧迫が加えられる方向に直交する水平方向について、関心組織である腫瘍全体が均一の力で圧迫されているか否か、言い換えれば圧迫の偏りがないかどうかを判断することができる。
つまり、作成したグラフにおいて、腫瘍全体が均一の力で圧迫されている場合は、全ての計測点ペア間の距離が略同一の周期及び振幅で変化するため、圧迫の偏りがないことが判る。一方、例えば図6(b)のように、対1〜対3の振幅にバラツキがあれば、振幅の大きい対3の方に偏って圧迫が加えられており、圧迫方向の見直しが必要である、などの判断をすることが可能である。
なお、各計測点の対に対応するグラフ表示は、図6(b)のように計測点の対とグラフとの対応関係がわかるように名前をつけて表示をしてもよいし、対応関係を色分け、その他の方法により表示してもよい。
図8は、本実施形態の超音波診断装置の特徴部の第4実施例を説明する図であり、計測点の設定例、及びこの場合に生成されるグラフを模擬的に示している。図9は、第4実施例の処理のフローチャートである。
図8(a)に示すように、計測点46が、腫瘍44の境界部の圧迫が加えられる方向に沿って2点設定されている点では第1実施例と同様であり、本実施例では、さらに、圧迫に対する変位の基準となる参照体54の圧迫が加えられる方向に沿って2点の計測点46が設定されている。
なお、参照体54は、例えば被検体の脂肪層、或いは探触子の圧迫面にあらかじめ設けられた圧迫板などを対象とすることができる。つまり、脂肪層は被検体間の個体差が比較的小さく、また圧迫に対する平均的な変位量があらかじめ既知となっており、圧迫板も圧迫に対する変位量があらかじめ既知となっているので、これを、圧迫に対する変位の基準となる参照体54として適用することができる。
変位演算部38は、各計測点をトラッキングなどによりフレームごとに追跡するとともに、腫瘍44の2点の計測点の対と、参照体54の2点の計測点の対のごとに、計測点間の距離を全フレームメモリのデータに対して演算し、その結果を逐次変位画像構成部40に出力する。変位画像構成部40は、2点の計測点の対ごとに、計測点間の距離変化の推移をグラフ55として生成する。
変位画像構成部40から出力されたグラフ55は、切替加算部24により加算合成され、図8(b)に示すように、断層画像と弾性画像が重畳された画像とともに画像表示器26に表示される。つまり、腫瘍44の計測点対を対4,参照体54の計測点対を対5とすると、対4,5のそれぞれに対する計測点間の距離変化の推移が縦に並べて表示される。
次に、図9のフローチャートを用いて、第4実施例の動作を説明する。まず、関心組織である腫瘍の境界に、白黒スキャンコンバータ22から得られたデータに対して2値化処理を行うことにより自動的に、もしくは操作部42を介して手動で、2点の計測点46を圧迫方向に沿って設定し、さらに同一圧迫方向の参照体54にも計測点46を設定する(S41)。
続いて、変位演算部38にて、腫瘍に設定された計測点46間の距離Aを検出するとともに、参照体に設定された計測点46間の距離Bを検出する(S42)。さらに、圧迫により変化する計測点の位置をトラッキングにより検出し、腫瘍及び参照体のそれぞれに対して、全てのフレームメモリにおける計測点間の距離(A)i,(B)iを検出する(S43)。続いて、腫瘍及び参照体のそれぞれに対して、変位画像構成部40にて(A)i,(B)iから組織の変位を表すグラフ55を生成し、切替加算部24を介して画像表示器26に表示する(S44)。
また、圧迫による参照体と腫瘍の歪みの比から、参照体を基準とした腫瘍の硬さを検出することが可能となる。つまり、それぞれのフレームメモリにおける腫瘍と参照体の計測点間の距離の比(A)i/(B)iを歪み比として算出し、画像表示器26に結果を表示してもよい。これにより、圧迫に対する変位量があらかじめ既知となっている参照体と比較することができるので、腫瘍の良悪性の判断をより簡便に行うことが可能となる。
本実施例のように参照体54を設ける場合、初期圧迫が適正か否かを判断する上で特に有効である。すなわち、超音波探触子による圧迫操作は、超音波探触子を被検体の体表に押し当ててある程度の初期圧迫を加え、この初期圧迫状態から、微小な加圧と減圧を繰り返し行われるのが一般的であるが、適切に弾性画像を生成するためには、初期圧迫の程度も重要な要素となる。
図10は、初期圧迫が適切か否かを説明するための模式的な図である。図10(a)は初期圧迫が強すぎた場合の例を示している。初期圧迫で強く圧迫を加えた状態でエラストグラフィの計測を行うと、柔らかい組織が既に潰された状態で微弱な圧迫が行われることとなり、組織の硬さを正確に表示できないことがある。
これに対して、図10(b)のように、圧迫を開始する前に参照体に例えば計測点を2点設定して、圧迫による計測点間の距離の推移をグラフ表示するとともに、適切な初期圧迫の範囲の閾値60を設定して表示しておけば、操作者はこれを参照して適切に初期圧迫を行うことができる。なお、計測点を1点のみ設定して計測点の変位の推移をグラフ表示してもよい。この場合は、これに対応した適切な閾値が別に設定される。また、図10では、参照体に設定された計測点に基づくグラフのみを表示しているが、図8に示すように、腫瘍に設定された計測点に基づくグラフと併せて表示することもできる。
このように、圧迫に対する変位があらかじめある程度定められている参照体に計測点を設定して、この計測点の変位に基づくグラフを表示することにより、このグラフを判断基準として初期圧迫量が適切になるように操作を行うことができる。
以上、本実施形態の超音波診断装置について説明したが、これには限定されない。例えば、計測点が設定されたら、白黒スキャンコンバータ22出力されるデータを用いて計測点を追跡して関心組織に加えられている圧迫状態の推移をグラフ化しているが、RF信号フレームデータを用いて追跡及びグラフ化を行うことも可能である。また、変位計測部30を、弾性画像の生成のためと、計測点の位置変化を検出するために共通して用いて、変位計測部30の出力に基づいて関心組織に加えられている圧迫状態の推移をグラフ化することも可能である。
また、計測点の代わりに、計測線、面などを計測対象として設定してもよい。また、超音波画像とグラフを表示器に並置して表示する例を説明したが、これに限らず、少なくともグラフを表示するように構成されていればよい。例えば、グラフのみを表示して、圧迫の統計解析を行い、熟練者と未熟練者の違いを明示して未熟練者の教育用として使用してもよい。

Claims (7)

  1. 被検体に圧迫を加えながら被検体との間で超音波を送受信する超音波探触子と、
    該超音波探触子で計測された反射エコー信号に基づいて前記被検体の断層部位のRF信号フレームデータを生成する整相加算部と、
    前記RF信号フレームデータに基づいて前記被検体の断層部位の超音波画像を生成する超音波画像生成部と、
    前記RF信号フレームデータ又は前記超音波画像に計測対象を設定する計測対象設定部と、
    設定された計測対象の位置変化を追跡する計測対象位置追跡部と、
    前記計測対象の位置変化に基づいて前記計測対象に加えられている圧迫状態の推移を表すグラフを生成する変位画像生成部と、
    前記生成されたグラフを表示する表示器と、を備え、
    前記計測対象設定部は、前記RF信号フレームデータ又は前記超音波画像の関心組織とその関心組織と異なる他の領域との境界に、前記圧迫が加えられる方向に沿った2点の計測点を設定し、
    前記変位画像生成部は、設定された2点の計測点間の距離変化の推移のグラフを生成し、
    前記超音波画像生成部によって生成された超音波画像と前記変位画像生成部によって生成されたグラフとを関連づけて前記表示器に表示される画像処理部を備えたことを特徴とする超音波診断装置。
  2. 前記超音波画像生成部は、前記被検体の断層部位の断層画像を生成する断層画像生成部と、取得時刻の異なる一対のRF信号フレームデータに基づいて前記被検体の断層部位の組織の硬さ又は軟らかさを表す弾性画像を生成する弾性画像構成部と、を含む請求項1に記載の超音波診断装置。
  3. 前記変位画像生成部は、前記計測対象に加えられている圧迫状態の推移のグラフの複数の極大点及び極小点に基づいてグラフの中点位置を求め、該グラフの中点位置から所定の振幅位置に振幅閾値のラインを重畳して前記表示器に表示するとともに、前記グラフ上に、あらかじめ設定された周期のラインを重畳して前記表示器に表示する請求項1に記載の超音波診断装置。
  4. 前記変位画像生成部は、前記グラフの前記周期のライン及び振幅閾値のラインの少なくとも一方を逸脱した箇所について前記表示器への表示態様を変更する請求項3の超音波診断装置。
  5. 前記計測対象設定部は、前記RF信号フレームデータ又は前記超音波画像の関心組織とその関心組織と異なる他の領域との境界に、前記圧迫が加えられる方向に沿った2点の計測点の対を、圧迫が加えられる方向に直交する方向の異なる位置に複数設定し、
    前記変位画像生成部は、それぞれの2点の計測点の対ごとに、設定された2点の計測点間の距離変化の推移のグラフを生成する請求項1の超音波診断装置。
  6. 前記計測対象設定部は、前記RF信号フレームデータ又は前記超音波画像の前記圧迫に対する変位の基準となる参照体に計測点を設定し、
    前記変位画像生成部は、前記参照体に設定された計測点の変化の推移のグラフを生成する請求項1の超音波診断装置。
  7. 超音波探触子により被検体に圧迫を加えながら被検体との間で超音波を送受信して計測された反射エコー信号に基づく前記被検体の断層部位のRF信号フレームデータから超音波画像を生成するとともに、前記超音波探触子による被検体の圧迫を補助する画像を生成する超音波画像処理方法であって、
    前記RF信号フレームデータ又は前記超音波画像に計測対象を設定するステップと、
    設定された計測対象の位置変化を追跡するステップと、
    前記計測対象の位置変化に基づいて前記計測対象に加えられている圧迫状態の推移を表すグラフを生成するステップと、
    生成されたグラフを表示するステップと、を含み、
    前記計測対象を設定するステップは、前記RF信号フレームデータ又は前記超音波画像の関心組織とその関心組織と異なる他の領域との境界に、前記圧迫が加えられる方向に沿った2点の計測点を設定し、
    前記圧迫状態の推移を表すグラフを生成するステップは、設定された2点の計測点間の距離変化の推移のグラフを生成し、
    前記生成された超音波画像と前記生成されたグラフとを関連づけて表示器に表示するステップを含むことを特徴とする超音波画像処理方法。
JP2009554356A 2008-02-21 2009-02-19 超音波診断装置、超音波画像処理方法 Active JP5437820B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2009554356A JP5437820B2 (ja) 2008-02-21 2009-02-19 超音波診断装置、超音波画像処理方法

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2008040688 2008-02-21
JP2008040688 2008-02-21
JP2009554356A JP5437820B2 (ja) 2008-02-21 2009-02-19 超音波診断装置、超音波画像処理方法
PCT/JP2009/052829 WO2009104657A1 (ja) 2008-02-21 2009-02-19 超音波診断装置、超音波画像処理方法及び超音波画像処理プログラム

Publications (2)

Publication Number Publication Date
JPWO2009104657A1 JPWO2009104657A1 (ja) 2011-06-23
JP5437820B2 true JP5437820B2 (ja) 2014-03-12

Family

ID=40985535

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2009554356A Active JP5437820B2 (ja) 2008-02-21 2009-02-19 超音波診断装置、超音波画像処理方法

Country Status (4)

Country Link
US (1) US8968200B2 (ja)
JP (1) JP5437820B2 (ja)
CN (1) CN101951840B (ja)
WO (1) WO2009104657A1 (ja)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4465535B2 (ja) * 2004-06-09 2010-05-19 株式会社日立メディコ 弾性画像表示方法及び超音波診断装置
US8968200B2 (en) * 2008-02-21 2015-03-03 Hitachi Medical Corporation Ultrasonic elastography for cross sectional tissue measurement method and apparatus
JP5665040B2 (ja) * 2009-09-10 2015-02-04 学校法人上智学院 変位計測方法及び装置、並びに、超音波診断装置
JP5501813B2 (ja) * 2010-03-17 2014-05-28 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置及びその制御プログラム
CN103908289B (zh) * 2012-12-31 2019-11-12 Ge医疗系统环球技术有限公司 消除剪切波中的背景噪声的方法和相应的超声成像系统
JP6144990B2 (ja) * 2013-07-31 2017-06-07 株式会社日立製作所 超音波画像撮像装置及び超音波画像撮像方法
US9613122B2 (en) * 2014-05-02 2017-04-04 Facebook, Inc. Providing eventual consistency for multi-shard transactions
US10678445B2 (en) * 2015-06-10 2020-06-09 Microsoft Technology Licensing, Llc Recovery in data centers
WO2017150368A1 (ja) * 2016-02-29 2017-09-08 コニカミノルタ株式会社 超音波診断装置及び超音波情報処理方法
US10893849B2 (en) * 2017-02-28 2021-01-19 Canon Medical Systems Corporation Ultrasound image diagnosis apparatus, medical image processing apparatus, and computer program product
CN111093526B (zh) * 2017-07-04 2022-03-04 富士胶片株式会社 声波诊断装置及声波诊断装置的工作方法
JP6836652B2 (ja) * 2017-08-25 2021-03-03 富士フイルム株式会社 音響波診断装置および音響波診断装置の制御方法
EP3723618B1 (en) 2017-12-11 2023-11-29 Hologic, Inc. Ultrasound localization system with advanced biopsy site markers
EP3818944A4 (en) * 2018-07-02 2021-09-01 FUJIFILM Corporation ACOUSTIC DIAGNOSTIC DEVICE AND METHOD OF CONTROLLING AN ACOUSTIC DIAGNOSTIC DEVICE

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006013916A1 (ja) * 2004-08-05 2006-02-09 Hitachi Medical Corporation 弾性像表示方法及び超音波診断装置
JP2007090003A (ja) * 2005-09-30 2007-04-12 Matsushita Electric Ind Co Ltd 超音波診断装置及びその制御方法

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4201396B2 (ja) 1998-08-20 2008-12-24 株式会社日立メディコ 超音波診断装置
JP3932485B2 (ja) 2003-05-30 2007-06-20 株式会社日立メディコ 超音波診断装置
IL158379A0 (en) * 2003-10-13 2004-05-12 Volurine Israel Ltd Non invasive bladder distension monitoring apparatus to prevent enuresis, and method of operation therefor
EP1711108A1 (en) 2004-01-23 2006-10-18 Koninklijke Philips Electronics N.V. Ultrasound imaging system for calculating a compression ratio
JP4515799B2 (ja) 2004-03-24 2010-08-04 株式会社日立メディコ 超音波診断装置
JP4465535B2 (ja) * 2004-06-09 2010-05-19 株式会社日立メディコ 弾性画像表示方法及び超音波診断装置
EP2168495B1 (en) * 2004-07-28 2012-03-14 Panasonic Corporation Ultrasonograph and ultrasonograph control method
WO2006121031A1 (ja) * 2005-05-09 2006-11-16 Hitachi Medical Corporation 超音波診断装置及び超音波画像表示方法
US7650025B2 (en) * 2005-08-01 2010-01-19 Siemens Medical Solutions Usa, Inc. System and method for body extraction in medical image volumes
US20070167784A1 (en) * 2005-12-13 2007-07-19 Raj Shekhar Real-time Elastic Registration to Determine Temporal Evolution of Internal Tissues for Image-Guided Interventions
US8358818B2 (en) * 2006-11-16 2013-01-22 Vanderbilt University Apparatus and methods of compensating for organ deformation, registration of internal structures to images, and applications of same
JP5322945B2 (ja) * 2007-11-06 2013-10-23 株式会社日立メディコ 超音波診断装置
US8968200B2 (en) * 2008-02-21 2015-03-03 Hitachi Medical Corporation Ultrasonic elastography for cross sectional tissue measurement method and apparatus

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006013916A1 (ja) * 2004-08-05 2006-02-09 Hitachi Medical Corporation 弾性像表示方法及び超音波診断装置
JP2007090003A (ja) * 2005-09-30 2007-04-12 Matsushita Electric Ind Co Ltd 超音波診断装置及びその制御方法

Also Published As

Publication number Publication date
US8968200B2 (en) 2015-03-03
CN101951840A (zh) 2011-01-19
CN101951840B (zh) 2014-05-14
JPWO2009104657A1 (ja) 2011-06-23
WO2009104657A1 (ja) 2009-08-27
US20110004100A1 (en) 2011-01-06

Similar Documents

Publication Publication Date Title
JP5437820B2 (ja) 超音波診断装置、超音波画像処理方法
US8734351B2 (en) Method of displaying elastic image and diagnostic ultrasound system
JP5559788B2 (ja) 超音波診断装置
JP5479353B2 (ja) 超音波診断装置
JP4898809B2 (ja) 超音波診断装置
JP5436533B2 (ja) 超音波診断装置及び弾性画像表示方法
JP5789593B2 (ja) 弾性画像の画質評価方法及び超音波診断装置
US7766832B2 (en) Ultrasonic diagnostic device and image processing device
WO2010024168A1 (ja) 超音波診断装置
WO2009098961A1 (ja) 超音波診断装置
JP5156421B2 (ja) 超音波診断装置
US9311704B2 (en) Ultrasonic diagnosis apparatus and image display method
JPWO2012114670A1 (ja) 超音波診断装置及び画像表示方法
JP5455592B2 (ja) 超音波診断装置、及び超音波画像表示方法
JPWO2011129237A1 (ja) 超音波診断装置
JP6230801B2 (ja) 超音波撮像装置及び超音波画像表示方法
JP4615528B2 (ja) 超音波診断装置

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120216

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20120216

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130204

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130403

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130917

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20131031

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20131212

R150 Certificate of patent or registration of utility model

Ref document number: 5437820

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

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

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