JP3401199B2 - 超音波映像化システムの超音波信号集束方法及び装置 - Google Patents

超音波映像化システムの超音波信号集束方法及び装置

Info

Publication number
JP3401199B2
JP3401199B2 JP31961998A JP31961998A JP3401199B2 JP 3401199 B2 JP3401199 B2 JP 3401199B2 JP 31961998 A JP31961998 A JP 31961998A JP 31961998 A JP31961998 A JP 31961998A JP 3401199 B2 JP3401199 B2 JP 3401199B2
Authority
JP
Japan
Prior art keywords
focusing
signal
ultrasonic
unit
apodization
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.)
Expired - Lifetime
Application number
JP31961998A
Other languages
English (en)
Other versions
JPH11221214A (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.)
Samsung Medison Co Ltd
Original Assignee
Samsung Medison Co Ltd
Medison Co Ltd
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 Samsung Medison Co Ltd, Medison Co Ltd filed Critical Samsung Medison Co Ltd
Publication of JPH11221214A publication Critical patent/JPH11221214A/ja
Application granted granted Critical
Publication of JP3401199B2 publication Critical patent/JP3401199B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K11/00Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/18Methods or devices for transmitting, conducting or directing sound
    • G10K11/26Sound-focusing or directing, e.g. scanning
    • G10K11/34Sound-focusing or directing, e.g. scanning using electrical steering of transducer arrays, e.g. beam steering
    • G10K11/341Circuits therefor
    • G10K11/346Circuits therefor using phase variation
    • 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/52046Techniques for image enhancement involving transmitter or receiver
    • 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/52046Techniques for image enhancement involving transmitter or receiver
    • G01S7/52049Techniques for image enhancement involving transmitter or receiver using correction of medium-induced phase aberration

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Medical Informatics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は超音波映像化システ
ムの超音波信号集束方法及び装置に係り、特に、対象体
に向って発射される超音波信号を送信及び/または受信
集束(focusing)する超音波信号集束方法及び装置に関
する。
【0002】
【従来の技術】一般に、超音波映像化システムは超音波
を利用して人体のような対象体の内部の断面構造を示
す。このシステムは超音波信号集束装置及び映像化装置
を備える。超音波信号集束装置は対象体に超音波信号を
発射し、対象体の音響インピーダンスの不連続面から反
射され戻ってくる超音波信号を電気信号に変換させる。
映像化装置はこの電気信号をディスプレーして対象体の
内部状態を見せてくれる。
【0003】このような超音波映像化システムの機能向
上のために求められる重要な要因の一つは超音波映像の
解像度であり、この解像度を改善するための開発がなさ
れつつある。超音波映像は軸方向(axial direction)
解像度よりは側方向(lateral)解像度が一層劣り、こ
の側方向解像度を決定するのには集束が最大の役割をは
たす。側方向解像度を向上させるため、超音波映像化シ
ステムは多数個の変換素子(transducer element)で形
成された配列型変換器(array transducer)をプローブ
(probe)として使用し、電気的な信号処理を用いた送
信集束(transmission focusing)及び受信集束(recei
pt focusing)を行なう。受信集束時には、超音波の進
行速度に鑑みて集束点(focal point)の位置を配列型
変換器に近い位置から遠い位置に連続的に移動させなが
ら集束を行なう動的集束(dynamic focusing)が使用可
能である。この動的集束は集束点を固定させた時より一
層良好な側面解像度を提供する。
【0004】一般の送信集束を図1に基づき説明すれば
次の通りである。図1において、x軸は配列型変換器1
0の超音波送信面に平行であり、z軸は超音波送信面に
垂直である。説明の便宜のため、配列型変換器10の特
定変換素子の表面と走査線(scanline)とが交わる点を
原点(0、0)に定め、原点に位置した変換素子を“変
換素子0"という。ここで、走査線はかじ取り角(steer
ing angle)が“θ"の場合に“変換素子0"の超音波送
信面の中央と送信集束点Fとを連結するラインである。
【0005】全ての変換素子から発射された超音波パル
スが送信集束点Fに同じ時点に到達するためには、原点
から遠く離れた変換素子ほど原点に近い変換素子より先
に超音波パルスを発射しなければならない。例えば、
“変換素子3"は“変換素子0"に比べてl3/vほど早
く超音波パルスを送信しなければならない。ここで、v
は超音波パルスの速度であり、rは原点から集束点Fま
での距離で、lは原点でない位置に置かれた変換素子の
超音波送信側面から集束点Fまでの距離(r+l)と距
離(r)間の距離差である。従って、z軸から反時計方
向に太線で示された走査線までのかじ取り角(θ)を利
用すれば、位置xlを有する“変換素子1"に対応する
遅延素子のための送信遅延時間td(x1)は次の数式
3で求められる。
【数3】
【0006】遅延素子に負数の送信遅延時間を与えると
いうことは不可能である。従って、実際は、全ての変換
素子に対応する送信遅延時間に同じ正(positive)の値
を足して全ての送信遅延時間が正数になるようにするこ
とが必要である。
【0007】しかし、説明の便宜のため、負数の送信遅
延時間が可能なことと仮定し、“変換素子0"が超音波
パルスを送信する瞬間の時間を“0"と定める。する
と、数式3におけるrをvt/2に置換することによっ
て、中心位置xを有する変換素子に対応する送信遅延時
間または受信遅延時間は次の数式4で表現される。
【数4】 ここで、tは超音波パルスが変換素子から集束点Fまで
を往復するのにかかる時間である。
【0008】パルス発生器(図示せず)から発生された
各電気パルスに前述した数式4の送信遅延時間を利用す
れば、変換素子が集束点Fに向かって送信する超音波パ
ルスは図2に示したようなリップル(ripples)を形成
する。図2において、各波(ripple)は、変換素子がパ
ルス発生器から出てくる一つのパルスを他の送信遅延時
間によって遅延させた後送信した結果である。
【0009】送信集束時には、全ての超音波パルスが同
じ時間に集束点Fに到達し、集束点Fに到達したこの超
音波パルスは全て同じ位相を持つ。従って、超音波の振
幅または強さは集束点Fで最大となる。しかし、集束点
F以外の位置(AまたはB)では、超音波パルスが同じ
時間に到達しなくなり、位相も相異なるため、超音波パ
ルスは互いに相殺されて超音波の強さは集束点Fに比べ
て小さくなる。この時、軸方向より側方向に遠ざかるほ
ど超音波パルスの強さはより一層小さくなる。
【0010】前述したような送受信集束の総合的な結果
により解像度が決定されるが、送信集束時には集束点を
固定させるしかないので、受信集束が理想的である。し
かし、受信集束の場合も、側方向解像度は集束点の近傍
で最良であり、その残り部分では比較的劣化する。
【0011】
【発明が解決しようとする課題】本発明は対象体のため
に用いられる超音波信号を走査線上の送受信可能な全て
の点で集束する超音波映像化システムの超音波信号集束
方法及び装置に関することであって、本発明の目的は超
音波映像化において、送信集束点から遠ざかるほど次第
に側方向解像度が低くなる既存の超音波映像化システム
の短所を克服するため、走査線上の全ての点に集束点が
置かれているように送受信集束を行って全体的な解像度
を向上させる超音波映像化システムの信号集束方法及び
装置を提供するにある。
【0012】
【課題を解決するための手段】前述したような目的を達
成するための本発明の特徴は、超音波信号を集束する装
置において、相異なる箇所に位置した変換素子によって
各々相異なる時間遅延を持つように作られて側方向の送
信集束を行なう配列型変換器と、前記変換素子に戻って
きた超音波信号を適切に増幅してデジタル信号に変換す
る受信部と、前記受信部からの信号を受けて送信する毎
に別途に格納(記憶)し、格納された超音波受信信号に
加えるべき時間遅延を計算して種々の送信走査線の信号
を一つの変換素子内で集束するイントラチャンネルプロ
セッシング部と、前記イントラチャンネルプロセッシン
グ部の出力を全て合算して全ての変換素子からの信号を
利用して集束を完成する受信集束部とを含む超音波信号
集束装置にある。
【0013】
【発明の実施の形態】以下、添付した図面に基づき本発
明の望ましい実施例を詳述する。図3は本発明に適用さ
れた超音波信号の集束を説明するための幾何学的な図面
である。
【0014】図3において、f0及びf2は集束点であ
り、P1、P1'、P2、P3及びP4は対象点であ
る。そして、Z1は“変換素子0"から走査線L0に沿
って集束点f0に向かって発射した超音波パルスが対象
点P1に到達する時までの経路の長さを、Z2は“変換
素子2"から走査線L2に沿って集束点f2に向かって
発射した超音波パルスが対象点P1に到達するのにかか
る時間間隔に対応する経路の長さを示したことである。
Zfは“変換素子2"から走査線L2に沿って集束点f
2に向かって発射した超音波パルスが集束点f2に到達
する時までの経路の長さを示し、Xfは“変換素子0"
と“変換素子2"との距離を表す。
【0015】図4は図3に示された走査線L0及びL2
に沿って送信された後、対象点P1から戻ってくる超音
波パルスを説明するための図面である。図4において、
(a)は走査線L0上の集束点f0に向かって送信した
超音波パルスが対象点P1から“変換素子0"に戻って
くる時間を示し、(b)は集束点f2に向かって送信し
た超音波パルスが対象点P1から“変換素子0"に戻っ
てくる時間を示す。
【0016】区間A内にある“変換素子−4"乃至“変
換素子4"は超音波パルスを、走査線L0上の集束点f
0に向かって送信集束されるように発射し、区間B内に
ある“変換素子−2"乃至“変換素子6"は超音波パルス
を、走査線L2上の集束点f2に向かって送信集束され
るように発射する。走査線L0について送信した超音波
パルスは“変換素子0"を出発してから距離Z1を進行
する時間が経過した後、対象点P1に到達する。その
後、距離Z1に該当する時間が過ぎた後“変換素子0"
に戻ってくる。走査線L2に対して送信した超音波パル
スは“変換素子2"を出発してから距離Z2を進行する
時間が経過した後対象点P1に到達する。対象点P1に
到達した超音波パルスは再び距離Z1に該当する時間が
過ぎた後“変換素子0"で受信される。“変換素子0"及
び“変換素子2"から集束点f0に向かって送信された
超音波パルスは送信集束により同じ時間に集束点f0に
到達する。従って、“変換素子0"から送信された超音
波パルスが対象点P1に到達する時間は“変換素子2"
から送信された超音波パルスが対象点P1'に到達する
時間と同一である。このような関係は、図3において集
束点f0を中心に有する円(点線で表示)で示した。
【0017】“変換素子0"で受信した二つの受信超音
波パルス即ち、走査線L0上の送信集束点f0に向かっ
て送信した超音波パルスから得られた受信超音波パルス
と走査線L2上の送信集束点f2に向かって送信した超
音波パルスから得られた受信超音波パルスは、対象点P
1から反射された受信超音波パルスが同じ時間に加えら
れるようにその遅延時間が調整される。こうすれば、送
信集束点f0に向かって送信して“変換素子0"で受信
した対象点P1に対応する受信超音波パルスは、送信集
束点f2に向かって送信して対象点P1から反射された
後“変換素子0"で受信される受信超音波パルスを利用
してその信号の強さが補強される。
【0018】このように全ての走査線に対して送信集束
を行なった信号を合成(synthesizing)すれば、結果的
に目標とする対象点から反射された超音波受信パルスは
全てエンベロープと位相が整列されて加えられるのでそ
の信号強度が大きくなるが、その他の点に対応する信号
強度はそうでなくなる。
【0019】図3を参照して信号強さ補強に必要な遅延
時間を説明すれば次の通りである。図3に示された“変
換素子0"、“変換素子2"及び対象点P1を考慮すれ
ば、次の数式5が満たされる。
【数5】
【0020】Z1=vt/2という関係を利用すれば、
数式5は次の数式6に置換される。
【数6】 この数式6は“変換素子2"から走査線L2上の集束点
f2に向かって送信した超音波パルスに対し、“変換素
子0"が対象点P1から受信する受信超音波パルスの遅
延時間になる。
【0021】超音波信号を送受信し、この信号を合成す
る時隣接した種々の走査線に送信集束して得られた信号
を全て同じ割合で合算する必要はない。側方向解像度を
向上させるためには色々な方法があるが、その中の一つ
として加重値を置く方法がある。例えば、走査線L0上
の送信集束点に向かって発射された信号から戻ってきた
信号のための加重値が‘1'とすれば、走査線L2上の
送信集束点に向かって発射された信号から戻ってくる信
号のための加重値は‘0.7'、それより一層離れてい
る走査線L3から戻ってきた信号のための加重値は
‘0.5'に定めて足す方法がある。このような方法を
ビームフォーミング技術ではアポダイゼーション(apod
ization)という。本発明の場合も適切にアポダイゼー
ションを使用して側方向解像度を改善できる。
【0022】今までは走査線L0のために“変換素子
0"だけ受信する場合を説明した。しかし、走査線L0
上の集束点f0に対する送信集束に参加できる全ての隣
接した変換素子(図3では“変換素子−4"乃至“変換
素子4")が走査線L0上の可能な全ての対象点に対す
る信号強度の補強に寄与できる。例えば、走査線L0に
沿って送信して“変換素子0"で受信した信号を‘基準
信号'と定める。すると、走査線L2に沿って送信した
信号を“変換素子1"が受信して格納(記憶)し、この
信号に対する適切な遅延時間を加えて基準信号と合せて
走査線L0上の対象点から得られた基準信号を補強でき
る。この対象点は勿論固定された点でなくても良く、超
音波パルスの進行を追って走査線L0上で動くことが一
層望ましい。こうすれば、走査線L0上の所望の全ての
対象点から得られる対象体からの信号は補強干渉により
大きくなるが、その他の点は相殺干渉により小さくなっ
て解像度が向上する。
【0023】図5は“変換素子0"に受信される超音波
パルスを補強するために隣接する“変換素子−1"に受
信される超音波パルスのために用いられる遅延時間を説
明するための図面である。説明の明瞭さのため、送信走
査線L2は支点(始点)(Xf、0)から集束点f2を
通り、集束点f2の座標は(Xf、Zf)である。ま
た、超音波パルスを受信する変換素子が支点(始点)
(Xe、0)にあり、最終的に集束する走査線は原点
(0、0)を通るL0と定める。すると、対象点P1か
ら“変換素子0"に受信される超音波パルスを補強する
ために“変換素子−1"に受信される超音波パルスに加
えなければならない遅延時間を数式で表現すれば次の数
式7の通りである。
【数7】
【0024】ここで、Ztとrrは次の数式8で表現で
きる。
【数8】 従って、遅延時間は次の数式9の通りである。
【数9】
【0025】P1をZp=vt/2に置き換えれば、数
式10のように変形できる。
【数10】
【0026】従って、支点(Xe、0)にある“変換素
子−1"で受信した超音波パルスをs(t、Xf、Z
f、Xe)とすれば、対象点P1に対する最終集束の結
果信号Sfocused(t)は、数式11の通りである。
【数11】 ここで、A(t、Xe、Xf)は前述したアポダイゼー
ション因子(apodization factor)であって、t、X
e、Xfにより各々適切な関数を使用する。
【0027】
【外1】 は対象点P1に対する超音波パルスの受信を許す全ての
送信スキャンランのための項であり、
【外2】 は対象点P1から反射される超音波パルスを受信する全
ての変換素子のための項である。この時、tもやはり前
述したように送信する毎に0に置く。また、Xe、Xf
は適切な間隔で適切な範囲内で変化させながら合算す
る。このような集束を行うためには最終集束された一本
の走査線のためにも多くのXe、Xfに対して一定の範
囲内の信号が必要である。しかし、最終集束された走査
線の位置が以前のものから全く移動しなければ、大部分
の信号はそのまま使用でき、一部が追加または削除され
るだけである。今までは説明の便宜上、最終集束する走
査線が固定されたように表現したが、実際は最終集束走
査線が移動可能であり、この時は今までの数式のx軸が
平行移動することと説明可能である。
【0028】図6は本発明の望ましい実施例にともなう
超音波映像化システムの集束装置を示す。送信部20か
ら発生された電気パルスは配列型変換器10の変換素子
に各々印加される。この時、各電気パルスは相異なる個
所に位置した変換素子によって各々相異なる時間遅延を
持つように作られ、側方向の送信集束になるようにす
る。変換素子は電気パルスの電圧によって超音波パルス
を発射する。超音波パルスは対象体内に伝播され、この
対象体内の支点から反射されて再び配列型変換器10の
変換素子に入射されて電気的な信号に変換される。受信
部40は変換素子からの超音波信号を適切に増幅してデ
ジタル信号に変換する。イントラチャンネルプロセッシ
ング部50は受信部40からのデジタル信号を受けて送
信する毎に別途に格納し、前述した数式11のXf項の
加算
【外3】 に該当する動作を遂行する。
【0029】受信集束部60は乗算器62、加算器64
及びチャンネル間アポダイゼーション値計算部66を具
備する。各乗算器62は対応するイントラチャンネルプ
ロセッシング部50から出力されるデータにチャンネル
間アポダイゼーション値計算部66から供給されるアポ
ダイゼーション値をかけ、その結果を加算器64に出力
する。加算器64は乗算器62から供給されるデータを
加算して、その結果である最終集束された値を表示部7
0へ供給する。表示部70は集束された超音波信号をモ
ニター80上にディスプレーさせる。なお、図6におい
て符号30は送信部20および表示部70の動作を制御
する制御部である。
【0030】図7は図6に示されたイントラチャンネル
プロセッシング部50の詳細な構成を示す。イントラチ
ャンネルプロセッシング部50は受信部40でデジタル
信号に変換されて印加されるRFデータを順序通り格納
(記憶)する多数のバッファメモリ51を備える。各バ
ッファメモリ51は対応する書込制御部53の制御下に
RFデータを格納し、読出制御部54の制御下に遅延補
間器52に格納されたデータを供給する。読出制御部5
4は多数のバッファメモリ51に対応するように構成さ
れ、集束遅延量計算部55から印加される遅延量によっ
てバッファメモリ51に格納されたデータを読出す。バ
ッファメモリ51は内部的に複数のSRAMで構成され
ており、交代に一つのSRAMにRFデータを記録し、
他の一つのSRAMで既に記録されたRFデータを出力
する。従って、記録及び読出の干渉がなく、読む間デー
タが変わらない。集束遅延量計算部55はバッファメモ
リ51に格納されたデータを集束するために前述した数
式10を利用して遅延時間を計算する。この数式10を
適用する過程自体は既存のデジタル集束における方法中
一つとほぼ同一である。例えば、数式10で示した集束
遅延時間(focusingdelay time)を受信部40でADC
のために用いられるサンプリング間隔(sampling inter
val)で分けて整数部と小数部とに分離する。集束遅延
量計算部55は整数データをバッファメモリ51の読出
アドレスに換算して読出制御部54に出力する。そし
て、集束遅延量計算部55は小数データを遅延補間器5
2に供給する。読出制御部54は該当読出アドレスを利
用してバッファメモリ51からデータを読み出し、読み
出したデータは遅延補間器52に出力される。遅延補間
器52はバッファメモリ51からのデータを集束遅延量
計算部55から印加される小数データによって時間遅延
及び補間したデータを発生する。このため、遅延補間器
52は集束遅延量計算部55からの小数データに対応す
る幾つかのデータをバッファメモリ51から供給される
データから決定し、決定されたデータを補間させる。そ
の結果、遅延補間器52は、バッファメモリ51からの
データに、集束遅延量計算部55において計算された遅
延量の小数部に該当する遅延時間を補間を通じて加え
る。補間は、広く知られた多相フィルタ(polyphase fi
lter)などを使用することができ、これらの動作は、前
後の幾つかのデータにそれぞれ適 切な係数を乗じて全て
を合わせたものである。このとき係数は、計算された遅
延量の小数部によって決定される。より詳しくは、遅延
補間器52は集束遅延量計算部55で計算された遅延量
の整数部に対応するデータを含んだいくつかのデータを
バッファメモリ51から読み出して小数部に対応する係
数を使ってデータを補間させる。その結果、遅延補間器
52はバッファメモリ51からのデータを集束遅延量計
算部55からの遅延量のうち小数部に該当する精密な時
間間隔だけ遅延させることができるようになる。。遅延
補間器52により発生されたデータはアポダイゼーショ
ン部90に出力される。アポダイゼーション部90内の
各乗算器92は対応する遅延補間器52から供給される
データにイントラチャンネルアポダイゼーション値計算
部96で計算された対応するアポダイゼーション値を掛
ける。前述した全ての過程は各送信走査線に該当する部
分で各々同時に行なわれ、この時発生された補間データ
は乗算器92を経て加算器94に供給される。加算器9
4は乗算器92から出力する信号を足して最終集束され
た信号を出力する。
【0031】今まで便宜上、全ての内容において線形配
列変換器におけるθ=0即ち、ビームかじ取り(beam s
teering)のない場合についてだけ説明したが、実際は
全ての集束について原理的に同じ方法を適用できる。例
えば、コンベックスアレイ(convex array)、かじ取り
(steering)、2次元アレイ(2-D array)などに前述
したような方法を適用できる。更に、本発明によると、
図2に示した焦点が配列型変換器から出射される超音波
パルスが発散する形態を有する波(ripple)様子を形成
するよう、配列型変換器の裏の任意の位置に置かれるよ
うにすることも可能である。
【0032】
【発明の効果】以上述べたように、本発明では集束点を
順次に変動させながら超音波の送信を行い、送信毎に全
ての変換素子が各々受信し、この受信信号を全て各々格
納し、格納された信号を利用して集束を行う。このよう
な方法で隣接した多くの走査線上に送信集束して得た信
号を全て活用すれば、受信しようとする走査線上の全て
の点について集束点を置いたことと同じ効果を得られ
る。
【図面の簡単な説明】
【図1】集束点Fに送信集束する既存の技術の原理を説
明するための概念図である。
【図2】従来の送信集束において、超音波パルスが集束
点Fに進行される姿を示す図面である。
【図3】本発明の望ましい一実施例にともなう超音波信
号集束を説明するための概念図である。
【図4】走査線L0及びL2に沿って送信した超音波パ
ルスが対象点に到達してから同じ変換素子に戻ってくる
ことを比較するための図である。
【図5】本発明の望ましい一実施例にともなう超音波信
号集束を説明するための概念図である。
【図6】本発明の超音波信号の集束装置を示した図であ
る。
【図7】図6の装置においてイントラチャンネルプロセ
ッシング部を示す図である。
【符号の説明】
50 イントラチャンネルプロセッシング部 51 バッファメモリ 52 遅延補間器 53 書込制御部 54 読出制御部 55 集束遅延量計算部 60 受信集束部 62 92 乗算器 64 94 加算器 66 96 アポダイゼーション値計算部 90 アポダイゼーション部
───────────────────────────────────────────────────── フロントページの続き (56)参考文献 特開 昭58−69538(JP,A) 特開 昭61−217143(JP,A) 特開 昭62−47348(JP,A) 特開 平6−14926(JP,A) 特開 平10−277042(JP,A) 特表 昭60−502141(JP,A) 国際公開96/3921(WO,A1) (58)調査した分野(Int.Cl.7,DB名) A61B 8/00

Claims (8)

    (57)【特許請求の範囲】
  1. 【請求項1】 超音波信号を集束する装置において、 相異なる箇所に位置した変換素子によって各々相異なる
    時間遅延を持つように作られて側方向の送信集束を行な
    う配列型変換器と、 前記変換素子に戻ってきた超音波信号を適切に増幅して
    デジタル信号に変換する受信部と、前記送信集束時毎に戻ってきた超音波信号を各送信走査
    線別に格納し、格納された超音波受信信号に加えるべき
    遅延時間を計算した後に前記遅延時間によって遅延させ
    た種々の送信走査線別受信信号を各変換素子別に 集束す
    るイントラチャンネルプロセッシング部と、 前記イントラチャンネルプロセッシング部の出力を全て
    合算して全ての変換素子からの信号を利用して集束を完
    成する受信集束部を含む超音波映像化システムの超音波
    信号集束装置。
  2. 【請求項2】 前記受信集束部は、 各チャンネル間アポダイゼーション値を計算するチャン
    ネル間アポダイゼーション値計算部と、 前記イントラチャンネルプロセッシング部から印加され
    る信号を前記チャンネル間アポダイゼーション値計算部
    から印加されるアポダイゼーション値によって乗算演算
    を行なう乗算器と、 前記乗算器から出力された信号を全て合算して最終集束
    された値を出力する加算器とを含むことを特徴とする
    求項1に記載の超音波映像化システムの超音波信号集束
    装置。
  3. 【請求項3】 前記イントラチャンネルプロセッシング
    部は、 前記各変換素子に受信された超音波信号をデジタル信号
    に変換したRFデータを走査線別に格納する多数のバッ
    ファメモリと、 前記バッファメモリに前記RFデータを格納するための
    制御信号を出力する書込制御部と、 前記バッファメモリに格納されたRFデータを集束する
    ための集束遅延時間を計算する集束遅延量計算部と、 前記集束遅延量計算部から印加される読出アドレスによ
    って前記バッファメモリに格納された超音波信号を読出
    すための制御信号を出力する読出制御部と、 前記集束遅延量計算部において計算された集束遅延量の
    小数部によって信号を補間して該信号を遅延量および小
    数部に該当する精密な量だけ遅延させる遅延補間器と、 前記遅延補間器によって遅延された信号に適切な加重値
    を加え、適切な加重値が加えられた信号を全て合算する
    アポダイゼーション部とを含むことを特徴とする請求項
    に記載の超音波信号集束装置。
  4. 【請求項4】 前記アポダイゼーション部は、 前記アポダイゼーション値を計算するイントラチャンネ
    ルアポダイゼーション値計部と、 前記イントラチャンネルプロセッシング部から印加され
    る信号を前記チャンネル間アポダイゼーション値計算部
    から印加されるアポダイゼーション値により適切な加重
    値に掛ける乗算演算を行う乗算器と、 前記乗算器から出力された信号を全て合算して一本の走
    査線についてアポダイゼーションが行なわれた値を出力
    する加算器とを含むことを特徴とする請求項3に記載の
    超音波信号集束装置。
  5. 【請求項5】 前記集束遅延量計算部は、送信集束して
    送信された信号が最終集束点に達した後に各々の変換素
    子に戻ってくる遅延時間を勘案して集束遅延時間を計算
    することを特徴とする請求項3に記載の超音波信号集束
    装置。
  6. 【請求項6】 前記集束遅延量計算部は、送信走査線が
    (Xf、0)を通り、その集束点の座標を(Xf、Z
    f)、受信する変換素子の座標を(Xe、0)、最終的
    に集束する走査線は原点(0、0)を通り、対象点で変
    換器の表面まで到達する経路の長さがvt/2である場
    合、前記送信走査線上の対象点からの信号を補強するた
    めに受信信号に加えるべき前記集束遅延時間td(t、
    xf、Zf、xe)を次式、 【数1】 により計算することを特徴とする請求項5に記載の超音
    波信号集束装置。
  7. 【請求項7】 前記イントラチャンネルプロセッシング
    部は全ての送信集束に対して個々の変換素子の超音波受
    信信号に前記集束遅延時間td(t、xf、Zf、x
    e)を加えた後に全て合算して集束することを特徴とす
    る請求項5に記載の超音波信号集束装置。
  8. 【請求項8】 前記イントラチャンネルプロセッシング
    部は、前記集束遅延時間td(t、xf、Zf、xe)
    を加えた後に全て合算して集束した最終集束信号を次式 【数2】 により計算し、このときのA(t、xe、xf)はt、
    xe、xfによって各々適切な関数を使用し、tは毎回
    送信毎に0に設定して計算する超音波信号集束装置であ
    ることを特徴とする請求項5に記載の超音波信号集束装
    置。
JP31961998A 1997-11-10 1998-11-10 超音波映像化システムの超音波信号集束方法及び装置 Expired - Lifetime JP3401199B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
KR97-59057 1997-11-10
KR1019970059057A KR100280197B1 (ko) 1997-11-10 1997-11-10 초음파영상화시스템의초음파신호집속방법및장치

Publications (2)

Publication Number Publication Date
JPH11221214A JPH11221214A (ja) 1999-08-17
JP3401199B2 true JP3401199B2 (ja) 2003-04-28

Family

ID=19524460

Family Applications (1)

Application Number Title Priority Date Filing Date
JP31961998A Expired - Lifetime JP3401199B2 (ja) 1997-11-10 1998-11-10 超音波映像化システムの超音波信号集束方法及び装置

Country Status (5)

Country Link
US (1) US6231511B1 (ja)
EP (1) EP0916966B1 (ja)
JP (1) JP3401199B2 (ja)
KR (1) KR100280197B1 (ja)
DE (1) DE69805484T2 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007072589A1 (ja) * 2005-12-19 2007-06-28 Jfe Steel Corporation 超音波による断面検査方法及び装置
KR101442221B1 (ko) * 2012-12-28 2014-09-19 서강대학교산학협력단 다상 보간필터를 이용한 멀티 빔 집속 장치 및 방법
US10197670B2 (en) 2014-10-30 2019-02-05 Seiko Epson Corporation Ultrasonic measurement apparatus, ultrasonic diagnostic apparatus, and ultrasonic measurement method

Families Citing this family (65)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20000038847A (ko) * 1998-12-09 2000-07-05 이민화 초음파영상화시스템에서의 집속방법
KR100330855B1 (ko) 1999-02-09 2002-04-03 이민화 초단수신신호의 저장과 저장된 데이터의 재사용이 가능한 디지털 초음파영상장치
US6514206B2 (en) * 2001-03-09 2003-02-04 Koninklijke Philips Electronics, N.V. Simultaneous fundamental and harmonic ultrasonic imaging
US6527723B2 (en) * 2001-06-26 2003-03-04 Koninklijke Philips Electronics N.V. Variable multi-dimensional apodization control for ultrasonic transducers
US6622561B2 (en) 2001-08-14 2003-09-23 Varco I/P, Inc. Tubular member flaw detection
US6578422B2 (en) 2001-08-14 2003-06-17 Varco I/P, Inc. Ultrasonic detection of flaws in tubular members
US6748808B2 (en) 2001-08-14 2004-06-15 Varco I/P, Inc. Flaw detection in tubular members
US6745136B2 (en) 2002-07-02 2004-06-01 Varco I/P, Inc. Pipe inspection systems and methods
ITSV20020040A1 (it) * 2002-08-13 2004-02-14 Esaote Spa Metodo e dispositivo per il rilevamento di immagini ecografiche
ES2277473B1 (es) * 2004-01-30 2008-07-16 Consejo Sup. Investig. Cientificas Composicion coherente de señales por correccion focal progresiva.
US7662114B2 (en) * 2004-03-02 2010-02-16 Focus Surgery, Inc. Ultrasound phased arrays
US8234923B2 (en) * 2004-09-20 2012-08-07 Innervision Medical Technologies Inc. Systems and methods for ultrasound imaging
US7850611B2 (en) * 2004-09-20 2010-12-14 Innervision Medical Technologies Inc. System and methods for improved ultrasound imaging
WO2007092054A2 (en) 2006-02-06 2007-08-16 Specht Donald F Method and apparatus to visualize the coronary arteries using ultrasound
US8473239B2 (en) 2009-04-14 2013-06-25 Maui Imaging, Inc. Multiple aperture ultrasound array alignment fixture
KR100947824B1 (ko) * 2006-03-29 2010-03-18 주식회사 메디슨 초음파 시스템에서 시그마-델타 아날로그-디지털 변환기를이용한 수신 집속 장치
WO2007133879A2 (en) * 2006-05-12 2007-11-22 Koninklijke Philips Electronics, N.V. Ultrasonic synthetic transmit focusing with motion compensation
US8317712B2 (en) * 2006-05-12 2012-11-27 Koninklijke Philips Electronics N.V. Eindhoven Retrospective dynamic transmit focusing for spatial compounding
EP2019626B1 (en) * 2006-05-12 2015-09-02 Koninklijke Philips Electronics, N.V. Incoherent retrospective dynamic transmit focusing
WO2008051639A2 (en) 2006-10-25 2008-05-02 Maui Imaging, Inc. Method and apparatus to produce ultrasonic images using multiple apertures
US8220334B2 (en) 2006-11-10 2012-07-17 Penrith Corporation Transducer array imaging system
US20080114246A1 (en) * 2006-11-10 2008-05-15 Penrith Corporation Transducer array imaging system
US7984651B2 (en) * 2006-11-10 2011-07-26 Penrith Corporation Transducer array imaging system
US9295444B2 (en) * 2006-11-10 2016-03-29 Siemens Medical Solutions Usa, Inc. Transducer array imaging system
KR100869496B1 (ko) * 2007-03-07 2008-11-21 주식회사 메디슨 초음파 영상 시스템 및 초음파 영상 형성 방법
US8968202B2 (en) 2007-04-13 2015-03-03 Medison Co., Ltd. System of forming ultrasound image and method of forming scan line data
KR100930569B1 (ko) 2007-04-13 2009-12-09 주식회사 메디슨 초음파 영상 시스템 및 스캔라인 데이터 형성 방법
CN101373181B (zh) * 2007-08-24 2012-03-21 深圳迈瑞生物医疗电子股份有限公司 实时计算逐点变迹系数的方法及装置
US9247926B2 (en) 2010-04-14 2016-02-02 Maui Imaging, Inc. Concave ultrasound transducers and 3D arrays
US9282945B2 (en) * 2009-04-14 2016-03-15 Maui Imaging, Inc. Calibration of ultrasound probes
EP2053420B1 (en) * 2007-10-25 2012-12-05 Samsung Medison Co., Ltd. Method of removing an effect of side lobes in forming an ultrasound synthetic image by motion estimation and compensation
US20100076314A1 (en) * 2008-03-25 2010-03-25 Robert Muratore System and method for creating virtual force field
JP5355924B2 (ja) * 2008-03-31 2013-11-27 株式会社東芝 超音波診断装置
DK200800633A (en) * 2008-05-02 2009-05-23 Bk Medical Aps Method and apparatus for processing ultrasonic signals
JP4776707B2 (ja) * 2009-03-30 2011-09-21 株式会社東芝 超音波画像化装置
CN101664321B (zh) * 2009-09-07 2011-06-15 无锡祥生科技有限公司 组织声速实时可调的超声诊断设备及其波束合成方法
US10069668B2 (en) * 2009-12-31 2018-09-04 Mediguide Ltd. Compensation of motion in a moving organ using an internal position reference sensor
KR101126181B1 (ko) * 2010-02-17 2012-03-22 서강대학교산학협력단 초음파 영상 시스템의 화소점에 수신 집속하는 방법 및 장치
JP6274724B2 (ja) 2010-02-18 2018-02-07 マウイ イマギング,インコーポレーテッド 多開口超音波撮像を用いた点音源送信及び音速補正
WO2012049124A2 (en) * 2010-10-11 2012-04-19 B-K Medical Aps Methods and systems for producing compounded ultrasound images
WO2012051305A2 (en) 2010-10-13 2012-04-19 Mau Imaging, Inc. Multiple aperture probe internal apparatus and cable assemblies
TW201336478A (zh) 2011-12-01 2013-09-16 Maui Imaging Inc 使用以回音為基及多孔徑都卜勒超音波之移動偵測
JP2015503404A (ja) 2011-12-29 2015-02-02 マウイ イマギング,インコーポレーテッド 任意経路のmモード超音波イメージング
KR101969305B1 (ko) * 2012-01-04 2019-08-13 삼성전자주식회사 초음파 영상 생성 장치 및 방법
EP2816958B1 (en) 2012-02-21 2020-03-25 Maui Imaging, Inc. Determining material stiffness using multiple aperture ultrasound
EP2833791B1 (en) 2012-03-26 2022-12-21 Maui Imaging, Inc. Methods for improving ultrasound image quality by applying weighting factors
WO2014026185A1 (en) 2012-08-10 2014-02-13 Maui Imaging, Inc. Calibration of multiple aperture ultrasound probes
US9986969B2 (en) 2012-08-21 2018-06-05 Maui Imaging, Inc. Ultrasound imaging system memory architecture
JP5873412B2 (ja) * 2012-09-28 2016-03-01 富士フイルム株式会社 超音波診断装置、音速決定方法およびプログラム
EP2730314B1 (en) * 2012-11-08 2016-02-03 Storz Medical Ag Endoscopic device for generating acoustic waves with variable focus
WO2014160291A1 (en) 2013-03-13 2014-10-02 Maui Imaging, Inc. Alignment of ultrasound transducer arrays and multiple aperture probe assembly
US8929177B2 (en) 2013-03-14 2015-01-06 Fujifilm Sonosite, Inc. System and method for performing progressive beamforming
JP5964774B2 (ja) * 2013-03-22 2016-08-03 富士フイルム株式会社 超音波診断装置、超音波診断装置の信号処理方法およびプログラム
US9883848B2 (en) 2013-09-13 2018-02-06 Maui Imaging, Inc. Ultrasound imaging using apparent point-source transmit transducer
KR102164456B1 (ko) * 2014-02-13 2020-10-12 삼성전자주식회사 초음파 측정 장치 및 초음파 측정 방법
JP6722656B2 (ja) 2014-08-18 2020-07-15 マウイ イマギング,インコーポレーテッド ネットワークベース超音波イメージングシステム
EP3260047A4 (en) 2015-02-18 2018-09-26 Hitachi, Ltd. Ultrasound image capturing device and method of processing ultrasound signal
CN113729764A (zh) 2016-01-27 2021-12-03 毛伊图像公司 具有稀疏阵列探测器的超声成像
JP6595360B2 (ja) * 2016-02-04 2019-10-23 株式会社日立製作所 超音波撮像装置および超音波撮像方法
JP6724797B2 (ja) * 2017-01-13 2020-07-15 コニカミノルタ株式会社 超音波診断装置、および、その制御方法
US11583254B2 (en) * 2017-06-28 2023-02-21 Duke University Systems and methods for intra-beam compounding for ultrasound imaging
US11751849B2 (en) 2017-09-27 2023-09-12 B-K Medical Aps High-resolution and/or high-contrast 3-D and/or 4-D ultrasound imaging with a 1-D transducer array
JP6848793B2 (ja) * 2017-09-28 2021-03-24 コニカミノルタ株式会社 超音波信号処理方法、及び超音波信号処理装置。
JP6526163B2 (ja) * 2017-11-28 2019-06-05 キヤノン株式会社 光音響トモグラフィの受信データ処理装置
CN110613476B (zh) * 2018-06-19 2022-09-23 青岛海信医疗设备股份有限公司 超声信号的处理方法和装置

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0208002B1 (de) * 1985-07-05 1990-01-03 Siemens Aktiengesellschaft Verfahren und Vorrichtung zur Ultraschallabtastung eines Objekts mit einem Ultraschallkopf
US4781199A (en) * 1985-01-07 1988-11-01 Kabushiki Kaisha Toshiba System and method for measuring sound velocity of internal tissue in an object being investigated
US4779623A (en) * 1985-08-09 1988-10-25 Kabushiki Kaisha Toshiba System for measuring sound velocity of tissue in an object and simultaneously providing a real-time B-mode image of the object
US4821574A (en) * 1986-07-22 1989-04-18 Kabushiki Kaisha Toshiba Method and apparatus for measuring ultrasonic velocity by crossed beam
DE3788757D1 (de) * 1986-08-20 1994-02-24 Siemens Ag Verfahren und Einrichtung zur adaptiven Fokussierung bei einem medizinischen Ultraschall-Bildgabegerät.
US4839652A (en) * 1987-06-01 1989-06-13 General Electric Company Method and apparatus for high speed digital phased array coherent imaging system
US4835689A (en) * 1987-09-21 1989-05-30 General Electric Company Adaptive coherent energy beam formation using phase conjugation
US5299577A (en) * 1989-04-20 1994-04-05 National Fertility Institute Apparatus and method for image processing including one-dimensional clean approximation
US5172343A (en) * 1991-12-06 1992-12-15 General Electric Company Aberration correction using beam data from a phased array ultrasonic scanner
US5230340A (en) * 1992-04-13 1993-07-27 General Electric Company Ultrasound imaging system with improved dynamic focusing
US5331964A (en) * 1993-05-14 1994-07-26 Duke University Ultrasonic phased array imaging system with high speed adaptive processing using selected elements
US5675554A (en) * 1994-08-05 1997-10-07 Acuson Corporation Method and apparatus for transmit beamformer
US5793701A (en) * 1995-04-07 1998-08-11 Acuson Corporation Method and apparatus for coherent image formation
US5605154A (en) * 1995-06-06 1997-02-25 Duke University Two-dimensional phase correction using a deformable ultrasonic transducer array

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007072589A1 (ja) * 2005-12-19 2007-06-28 Jfe Steel Corporation 超音波による断面検査方法及び装置
JP2007170871A (ja) * 2005-12-19 2007-07-05 Jfe Steel Kk 超音波による断面検査方法及び装置
KR101442221B1 (ko) * 2012-12-28 2014-09-19 서강대학교산학협력단 다상 보간필터를 이용한 멀티 빔 집속 장치 및 방법
US10197670B2 (en) 2014-10-30 2019-02-05 Seiko Epson Corporation Ultrasonic measurement apparatus, ultrasonic diagnostic apparatus, and ultrasonic measurement method

Also Published As

Publication number Publication date
KR100280197B1 (ko) 2001-02-01
KR19990039095A (ko) 1999-06-05
EP0916966A1 (en) 1999-05-19
JPH11221214A (ja) 1999-08-17
DE69805484T2 (de) 2003-01-09
US6231511B1 (en) 2001-05-15
DE69805484D1 (de) 2002-06-27
EP0916966B1 (en) 2002-05-22

Similar Documents

Publication Publication Date Title
JP3401199B2 (ja) 超音波映像化システムの超音波信号集束方法及び装置
JP3700990B2 (ja) ビーム形成器及び超音波イメージング・システム
JP2777197B2 (ja) 超音波診断装置
US8672846B2 (en) Continuous transmit focusing method and apparatus for ultrasound imaging system
EP1004894B1 (en) Method and apparatus for high-frame-rate high-resolution ultrasonic image data acquisition
US8545406B2 (en) Dynamic aperture control and normalization for apodization in beamforming
US20060173313A1 (en) Coherence factor adaptive ultrasound imaging
JPH10293170A (ja) 超音波ビーム形成装置
Hazard et al. Theoretical assessment of a synthetic aperture beamformer for real-time 3-D imaging
US5793703A (en) Digital time-delay acoustic imaging
JP2001187054A (ja) 超音波ビーム経路の数値的最適化方式
CN109363714B (zh) 一种超声成像设备及其超声成像方法
JPH11197151A (ja) 超音波イメージング・システム用のbモード処理装置および検出後像処理方法
CN108024798B (zh) 超声波诊断装置以及延迟数据生成方法
EP0834750A2 (en) Method and apparatus for forming ultrasonic three-dimensional images using cross array
US20140257103A1 (en) Three Dimensional (3D) Transverse Oscillation Vector Velocity Ultrasound Imaging
JP2588185B2 (ja) 超音波診断装置
US8926514B2 (en) Iterative time delay values for ultrasound beamforming
JP3492095B2 (ja) 超音波撮像装置
CN117825512B (zh) 面向钢轨轨头缺陷的高速超声动态成像检测方法及装置
JPH0560B2 (ja)
JPH09108223A (ja) 超音波診断装置
JP2000316854A (ja) 超音波装置
JP4593260B2 (ja) 超音波診断装置
JP2000300564A (ja) 超音波診断装置

Legal Events

Date Code Title Description
R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080221

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090221

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100221

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100221

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110221

Year of fee payment: 8

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120221

Year of fee payment: 9

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130221

Year of fee payment: 10

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140221

Year of fee payment: 11

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

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

EXPY Cancellation because of completion of term