[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

JPH07178095A - Ultrasonic diagnostic system - Google Patents

Ultrasonic diagnostic system

Info

Publication number
JPH07178095A
JPH07178095A JP5323484A JP32348493A JPH07178095A JP H07178095 A JPH07178095 A JP H07178095A JP 5323484 A JP5323484 A JP 5323484A JP 32348493 A JP32348493 A JP 32348493A JP H07178095 A JPH07178095 A JP H07178095A
Authority
JP
Japan
Prior art keywords
time
change
signal
level
luminance signal
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
Application number
JP5323484A
Other languages
Japanese (ja)
Other versions
JP3405578B2 (en
Inventor
Yoshitaka Watanabe
欣孝 渡辺
Takanobu Uchibori
孝信 内堀
Isao Uchiumi
勲 内海
Yasushi Ueki
康至 上木
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.)
Toshiba Corp
Canon Medical Systems Corp
Original Assignee
Toshiba Corp
Toshiba Medical Systems Engineering 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 Toshiba Corp, Toshiba Medical Systems Engineering Co Ltd filed Critical Toshiba Corp
Priority to JP32348493A priority Critical patent/JP3405578B2/en
Publication of JPH07178095A publication Critical patent/JPH07178095A/en
Application granted granted Critical
Publication of JP3405578B2 publication Critical patent/JP3405578B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

PURPOSE:To diagnose the blood streaming condition quickly in good performance when a contrast agent is injected in a blood vessel by calculating the change of a brightness signal with time, and displaying the obtained time change together with a tomographic image. CONSTITUTION:According to a delay control of a delay line 24, a drive voltage is fed from a pulser 22 in an analog part 2 and impressed on each piezo vibrator of an ultrasonic probe 1, and an ultrasonic beam is transmitted, and then the reflected waves from inside of a vital organism will be received by probe 11. The reception signal is sent to a wave detector 26 via a preamplifier 21, the delay line 24, and an adder 25, and therein the envelope of reception signals is detected, and the result is passed to a scan converting means 3 as a brightness signal via an A/D converter. Also, the brightness signal is fed to a signal level time-change extractor 6 to determine the time change of the average for all points within the interest region, and also the level change of the brightness signal is determined. Each time change is converted into color information according to its level, and display is made while it is overlapped on the tomographic image as a two-dimensional image.

Description

【発明の詳細な説明】Detailed Description of the Invention

【0001】[0001]

【産業上の利用分野】本発明は、被検体内を超音波ビー
ムで走査することにより得られた受信信号の包絡線を検
波し、この検波信号を用いて断層像を表示する超音波診
断装置に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an ultrasonic diagnostic apparatus for detecting an envelope of a received signal obtained by scanning the inside of a subject with an ultrasonic beam and displaying a tomographic image using this detected signal. Regarding

【0002】[0002]

【従来の技術】近年、超音波診断の分野では、X線診断
における造影剤と同様の機能としての超音波反射特性が
良好なコントラスト剤の開発が進み、実用域に入ってい
る。このコントラスト剤を例えば血管に注入すると、B
モード画像に血流の様子がその周辺組織と異なるコント
ラストで鮮明に表れる。
2. Description of the Related Art In recent years, in the field of ultrasonic diagnosis, a contrast agent having a good ultrasonic reflection characteristic as a function similar to that of a contrast agent in X-ray diagnosis has been developed and is in a practical range. If this contrast agent is injected into a blood vessel, for example, B
The appearance of blood flow clearly appears in the mode image with a contrast different from that of the surrounding tissue.

【0003】ところで、血流が関心部分に流入してから
流出するまでの時間的変化には貴重な情報が含まれてい
る。しかし従来の超音波診断装置では、上記時間的変化
を計算する機能を備えていないため、検査終了後に複数
枚のBモード画像データを時間的変化を計算する処理装
置に転送して、そこで計算するようにしていた。このた
め、診断結果が出されるまでに長時間必要であった。
By the way, valuable information is included in the temporal change from the inflow of blood flow to the outflow of the region of interest. However, since the conventional ultrasonic diagnostic apparatus does not have the function of calculating the temporal change, a plurality of B-mode image data are transferred to the processing device for calculating the temporal change after the inspection is completed, and the B-mode image data are calculated there. Was doing. For this reason, it took a long time before the diagnosis result was issued.

【0004】[0004]

【発明が解決しようとする課題】そこで、本発明は、上
述した事情に対処するべく、Bモード画像の生成に並行
して輝度の時間的変化をリアルタイムで計算する超音波
診断装置を提供することを目的とする。
SUMMARY OF THE INVENTION Therefore, the present invention provides an ultrasonic diagnostic apparatus for coping with the above-mentioned circumstances, in which a temporal change in luminance is calculated in real time in parallel with the generation of a B-mode image. With the goal.

【0005】[0005]

【課題を解決するための手段】本発明による被検体内を
超音波ビームで走査することにより得られた受信信号の
包絡線を検波し、この検波結果を輝度信号として断層像
を表示する超音波診断装置は、前記輝度信号を入力し、
前記輝度信号の時間的変化を計算する手段と、前記時間
的変化を前記断層像と共に表示する手段とを具備する。
An ultrasonic wave for detecting the envelope of a received signal obtained by scanning the inside of a subject with an ultrasonic beam according to the present invention and displaying a tomographic image using the detected result as a luminance signal. The diagnostic device inputs the luminance signal,
It is provided with means for calculating the temporal change of the luminance signal and means for displaying the temporal change together with the tomographic image.

【0006】本発明による被検体内を超音波ビームで走
査することにより得られた受信信号の包絡線を検波し、
この検波結果を輝度信号として断層像を表示する他の超
音波診断装置は、前記輝度信号を入力し、前記検波信号
の時間的変化を前記断層像の各点毎に計算する手段と、
前記各時間的変化をそのレベルに応じてカラー情報に変
換し2次元画像として前記断層像に重ねて表示する手段
とを具備する。
The envelope of the received signal obtained by scanning the inside of the subject with the ultrasonic beam according to the present invention is detected,
Another ultrasonic diagnostic apparatus for displaying a tomographic image as a luminance signal with this detection result is a means for inputting the luminance signal and calculating a temporal change of the detection signal for each point of the tomographic image,
A means for converting each temporal change into color information according to its level and displaying it as a two-dimensional image by superimposing it on the tomographic image.

【0007】[0007]

【作用】本発明による超音波診断装置によれば、輝度信
号を取り込んでその時間的変化を計算して、この時間的
変化を断層像と共にリアルタイムで表示することができ
るので、コントラスト剤を血管に注入した場合には血流
状態を短時間のうちに診断することができる。
According to the ultrasonic diagnostic apparatus of the present invention, it is possible to capture a luminance signal, calculate its temporal change, and display this temporal change together with a tomographic image in real time. When injected, the blood flow condition can be diagnosed in a short time.

【0008】本発明による他の超音波診断装置によれ
ば、輝度信号を取り込んでその時間的変化を各点毎に計
算して、各時間的変化をそのレベルに応じてカラー情報
に変換し2次元画像として断層像に重ねてリアルタイム
で表示することができるので、コントラスト剤を血管に
注入した場合には血流状態を短時間のうちに診断するこ
とができ、またその2次元的な状態も把握することがで
きる。
According to another ultrasonic diagnostic apparatus of the present invention, a luminance signal is taken in, its temporal change is calculated for each point, and each temporal change is converted into color information according to its level. Since it can be displayed in real time as a three-dimensional image overlaid on a tomographic image, when a contrast agent is injected into a blood vessel, the blood flow state can be diagnosed in a short time, and the two-dimensional state can also be obtained. You can figure it out.

【0009】[0009]

【実施例】図面を参照して本発明の一実施例について説
明する。図1は第1実施例のブロック図である。セクタ
電子走査型超音波プローブ1は複数の圧電振動子が一次
元に配列される。このプローブ1には送受信制御と受信
信号の検波とを兼ねたアナログ部2が接続される。この
アナログ部2にはパルサ22が含まれる。
An embodiment of the present invention will be described with reference to the drawings. FIG. 1 is a block diagram of the first embodiment. The sector electronic scanning ultrasonic probe 1 has a plurality of piezoelectric vibrators arranged one-dimensionally. The probe 1 is connected to an analog section 2 which has both transmission and reception control and detection of a received signal. The analog section 2 includes a pulser 22.

【0010】パルサ22はディレーライン24から供給
される各圧電振動子毎の遅延情報にしたがって異なるタ
イミングで各圧電振動子に駆動パルスを供給する。この
駆動パルスによりプローブ1から超音波ビームが上記遅
延情報に応じた方向に送信される。この遅延情報をラス
タ毎に変化させることにより超音波ビームは扇状に移動
する。超音波ビームは被検体内の音響インピーダンスの
境界で反射する。この反射波は、プローブ1の各振動子
で受信される。各受信信号はプリアンプ21を介してデ
ィレーライン24に個別に送られ、そこで所定の焦点深
度に応じた遅延時間を与えられた後、加算器25で加算
される。この加算信号は、検波器26に送られ、そこで
対数増幅および包絡線が検波される。この検波結果は輝
度信号として一般にはディジタルスキャンコンバータと
呼ばれる走査変換手段3と共に、信号レベル時間変化抽
出器6に送られる。
The pulsar 22 supplies a drive pulse to each piezoelectric vibrator at different timings according to the delay information for each piezoelectric vibrator supplied from the delay line 24. With this drive pulse, the ultrasonic beam is transmitted from the probe 1 in the direction according to the delay information. By changing this delay information for each raster, the ultrasonic beam moves in a fan shape. The ultrasonic beam is reflected at the boundary of the acoustic impedance inside the subject. This reflected wave is received by each transducer of the probe 1. Each reception signal is individually sent to the delay line 24 via the preamplifier 21, where it is given a delay time corresponding to a predetermined depth of focus and then added by the adder 25. This added signal is sent to the detector 26, where the logarithmic amplification and the envelope are detected. The detection result is sent as a luminance signal to the signal level time change extractor 6 together with the scan conversion means 3 generally called a digital scan converter.

【0011】信号レベル時間変化抽出器6には、超音波
ビームによる走査の繰り返しに応じて図示しないマウス
やトラックボール等の関心領域設定部により設定された
関心領域に相当する輝度信号が2次元走査面の走査毎に
検波器26から連続的に供給されてくる。信号レベル時
間変化抽出器6には縦軸を輝度レベル、横軸を時間とし
たグラフ領域を内部メモリに備えている。関心領域の輝
度信号が供給される毎に、その輝度信号の輝度レベルお
よび走査時刻に応じたグラフ領域の位置に点が分布され
る。このグラフ情報は、次の輝度信号が供給されるまで
に作成され、走査変換手段3に出力される。
The signal level temporal change extractor 6 two-dimensionally scans a luminance signal corresponding to a region of interest set by a region of interest setting unit such as a mouse or a trackball (not shown) in response to repeated scanning with an ultrasonic beam. It is continuously supplied from the detector 26 every time the surface is scanned. The signal level time change extractor 6 is provided with a graph area in which the vertical axis represents the brightness level and the horizontal axis represents time in the internal memory. Each time the luminance signal of the region of interest is supplied, points are distributed at the positions of the graph region according to the luminance level of the luminance signal and the scanning time. This graph information is created by the time the next luminance signal is supplied and is output to the scan conversion means 3.

【0012】走査変換手段3には表示手段5の表示画面
に対応する2次元メモリが備えられている。この2次元
メモリには検波器26からの輝度信号を記憶する領域
と、信号レベル時間変化抽出器6からのグラフ情報を記
憶する領域とが予め区別して割り当てられている。走査
変換手段3は輝度信号やグラフ情報を一旦ディジタル信
号として記憶する。走査変換手段3から表示手段5の走
査方式にしたがって読み出された輝度信号やグラフ情報
は、D/A変換手段4を介してモニタとしての表示手段
5に送られ、そこで関心領域に関する輝度信号の強度
(レベル)の時間的変化がグラフとしてBモード画像
(断層像)と別々に表示される。
The scan conversion means 3 is provided with a two-dimensional memory corresponding to the display screen of the display means 5. An area for storing the luminance signal from the detector 26 and an area for storing the graph information from the signal level time change extractor 6 are separately assigned to the two-dimensional memory in advance. The scan conversion means 3 temporarily stores the luminance signal and the graph information as a digital signal. The luminance signal and the graph information read from the scanning conversion means 3 according to the scanning method of the display means 5 are sent to the display means 5 as a monitor via the D / A conversion means 4, and there, of the luminance signal relating to the region of interest. The temporal change in intensity (level) is displayed as a graph separately from the B-mode image (tomographic image).

【0013】次に本実施例の動作について説明する。X
線診断における造影剤と同様の機能としての超音波反射
特性が良好なコントラスト剤を例えば血管に注入した
後、以下のような動作が開始される。
Next, the operation of this embodiment will be described. X
After injecting into a blood vessel, for example, a contrast agent having a good ultrasonic reflection characteristic as a function similar to that of a contrast agent in line diagnosis, the following operation is started.

【0014】アナログ部2のパルサ22からプローブ1
1の各圧電振動子にディレーライン24の遅延制御にし
たがって駆動電圧がされる。これにより超音波ビームが
被検体内に送信される。生体中で反射してきた反射波
は、プローブ11で受信される。この受信信号は、アナ
ログ部2のプリアンプ21、ディレーライン24、加算
器25を順に介して検波器26に送られる。受信信号は
その包絡線を検波器26で検波され、図示しないアナロ
グディジタル変換器を介して輝度信号として走査変換手
段3に送られる。走査変換手段3では輝度信号は当該ラ
スタの方位角にしたがってメモリ領域に記憶される。
From the pulsar 22 of the analog section 2 to the probe 1
A drive voltage is applied to each piezoelectric vibrator 1 according to the delay control of the delay line 24. Thereby, the ultrasonic beam is transmitted into the subject. The reflected wave reflected in the living body is received by the probe 11. This received signal is sent to the detector 26 through the preamplifier 21, the delay line 24, and the adder 25 of the analog section 2 in order. The envelope of the received signal is detected by the detector 26, and is sent to the scan conversion means 3 as a luminance signal via an analog-digital converter (not shown). In the scan conversion means 3, the luminance signal is stored in the memory area according to the azimuth angle of the raster.

【0015】検波器26からの輝度信号は信号レベル時
間変化抽出器6にも供給される。信号レベル時間変化抽
出器6では、関心領域設定部の操作により設定された少
なくとも1点(1画素)の関心領域の輝度信号の時間経
過に伴うレベル変化をグラフとして作成するが、複数点
(複数画素)を含む関心領域が設定された場合、本実施
例では、関心領域内全点の平均値の時間的変化を求める
こと、関心領域内の最大値及び/又は最小値の時間的変
化を求めること、関心領域内全点の平均値の所定の短時
間間隔における変化率の時間的変化を求めることのいず
れかにより時間経過に伴う輝度信号のレベル変化を求め
る。複数点を含む関心領域内全点の平均値を求める場
合、各点に対し窓関数による重み付けを行うこともあ
る。
The luminance signal from the detector 26 is also supplied to the signal level time change extractor 6. The signal level time change extractor 6 creates a level change over time of the luminance signal of at least one point (one pixel) of the region of interest set by the operation of the region of interest setting unit as a graph. When a region of interest including (pixels) is set, in this embodiment, the temporal change of the average value of all the points in the region of interest is obtained, and the temporal change of the maximum value and / or the minimum value in the region of interest is obtained. The change in the level of the luminance signal with the passage of time is obtained by either obtaining the change over time in the rate of change of the average value of all the points in the region of interest in a predetermined short time interval. When obtaining the average value of all the points in the region of interest including a plurality of points, each point may be weighted by a window function.

【0016】また、信号レベル時間変化抽出器6は、1
心拍期間のグラフから輝度レベルが最大値を示す時刻か
らその最大値の半値になる時刻までの時間幅をレベル変
化を示す数値情報として求める。
Further, the signal level time change extractor 6 is
From the graph of the heartbeat period, the time width from the time when the brightness level shows the maximum value to the time when it becomes half the maximum value is obtained as the numerical information showing the level change.

【0017】図2には関心領域内全点の平均値の時間的
変化を求める場合の信号レベル時間変化抽出器6による
2次元走査面の1回の走査動作に対応する動作手順が順
番に示されている。なお、この動作は2次元走査面を一
通り走査する毎に繰り返される。
FIG. 2 shows in sequence the operation procedure corresponding to one scanning operation of the two-dimensional scanning plane by the signal level time change extractor 6 when obtaining the time change of the average value of all the points in the region of interest. Has been done. This operation is repeated every time the two-dimensional scanning surface is scanned once.

【0018】ステップ(S1)で、1回の2次元走査面
分の輝度信号が検波器26から信号レベル時間変化抽出
器6に供給される。ステップ(S2)では、1回の2次
元走査面分の輝度信号の中の表示手段5に表示されたB
モード画像上に関心領域設定部の操作により設定された
少なくとも1点(1画素)の関心領域の輝度信号だけが
信号レベル時間変化抽出器6の内部メモリに格納され、
関心領域以外の輝度信号は破棄される。
In step (S1), the luminance signal for one two-dimensional scanning plane is supplied from the detector 26 to the signal level time change extractor 6. In step (S2), B displayed on the display means 5 in the luminance signal for one two-dimensional scanning plane is displayed.
Only the luminance signal of at least one point (one pixel) of the ROI set on the mode image by the operation of the ROI setting unit is stored in the internal memory of the signal level time change extractor 6,
Luminance signals other than the region of interest are discarded.

【0019】ステップ(S3)では、内部メモリに格納
されている関心領域内の輝度信号のレベルの平均値が以
下の(1)式にしたがって計算される。なお、I´を輝
度レベルの平均値とし、x,yそれぞれx方向,y方向
の画素の座標(関心領域x0〜x1 ,y0 〜y1 の範囲
内の座標)とし、I(x,y) は座標(x,y)の輝度レベ
ルとする。
In step (S3), the average value of the luminance signal levels in the region of interest stored in the internal memory is calculated according to the following equation (1). It should be noted that I'is the average value of the brightness level, x and y are the coordinates of the pixels in the x direction and the y direction (the coordinates within the region of interest x0 to x1, y0 to y1), and I (x, y) is The brightness level is at the coordinate (x, y).

【0020】[0020]

【数1】 [Equation 1]

【0021】この平均値はステップ(S4)で内部メモ
リに記憶された後、ステップ(S5)で、所定時間前か
ら直前までにすでに作成されているグラフ情報に、当該
平均値のレベル及び走査時刻に応じた位置に点が追加さ
れて、最新のグラフ情報が作成される。この最新のグラ
フ情報は、ステップ(S6)で走査変換手段3に送ら
れ、走査変換手段3のメモリ領域における輝度信号を記
憶する領域以外の専用領域に記憶される。
This average value is stored in the internal memory in step (S4), and then in step (S5), the level of the average value and the scanning time are added to the graph information already created from the time before the predetermined time to the time immediately before. The point is added to the position corresponding to and the latest graph information is created. This latest graph information is sent to the scan conversion means 3 in step (S6) and is stored in a dedicated area other than the area for storing the luminance signal in the memory area of the scan conversion means 3.

【0022】走査変換手段3では、同じ2次元走査によ
る受信信号から作成されたBモード画像とグラフ情報を
順次走査によりD/A変換手段4に読出し、図3(a)
に示すように表示手段5の1画面内にBモード画像とグ
ラフ情報を別々に表示させる。なお図3(b)には図3
(a)のグラフ情報の拡大図である。
In the scanning conversion means 3, the B-mode image and the graph information created from the received signals by the same two-dimensional scanning are sequentially scanned and read out to the D / A conversion means 4, as shown in FIG.
As shown in, the B mode image and the graph information are separately displayed in one screen of the display means 5. Note that FIG.
It is an enlarged view of the graph information of (a).

【0023】このように関心領域の輝度信号の時間変化
曲線がBモード画像と共にリアルタイムで表示される。
図4には関心領域内の最大値及び/又は最小値の時間的
変化を求める場合の信号レベル時間変化抽出器6による
2次元走査面の1回の走査動作に対応する動作手順が順
番に示されている。なお、この動作は2次元走査面を一
通り走査する毎に繰り返される。図4において図2と同
じステップは同符号を付して説明は省略する。ステップ
(S13)では、内部メモリに格納されている関心領域
内の全ての点の輝度信号の中の最大値及び最小値が選択
される。これらの最大値及び最小値はステップ(S
4)、ステップ(S5)、ステップ(S6)を経過し
て、最大値及び最小値の時間変化を現すグラフ情報とし
て走査変換手段3に送られ、図3(a)と同様にBモー
ド画像と共に同時表示される。図5には表示手段5に表
示されるグラフ情報を示している。
In this way, the time change curve of the luminance signal of the region of interest is displayed in real time together with the B-mode image.
FIG. 4 sequentially shows an operation procedure corresponding to one scanning operation of the two-dimensional scanning plane by the signal level time change extractor 6 when obtaining the time change of the maximum value and / or the minimum value in the region of interest. Has been done. This operation is repeated every time the two-dimensional scanning surface is scanned once. 4, the same steps as those in FIG. 2 are designated by the same reference numerals and the description thereof will be omitted. In step (S13), the maximum value and the minimum value of the luminance signals of all the points in the region of interest stored in the internal memory are selected. These maximum and minimum values are calculated in step (S
4), step (S5) and step (S6) are passed, and the graph information showing the time change of the maximum value and the minimum value is sent to the scan conversion means 3 and is sent together with the B mode image as in FIG. 3 (a). Displayed simultaneously. FIG. 5 shows the graph information displayed on the display means 5.

【0024】図6には関心領域内全点の平均値の所定の
短時間間隔における変化率の時間的変化を求める場合の
信号レベル時間変化抽出器6による2次元走査面の1回
の走査動作に対応する動作手順が順番に示されている。
なお、この動作は2次元走査面を一通り走査する毎に繰
り返される。図6において図2と同じステップは同符号
を付して説明は省略する。
FIG. 6 shows one scanning operation of the two-dimensional scanning plane by the signal level time change extractor 6 when the time change of the change rate of the average value of all the points in the region of interest in a predetermined short time interval is obtained. The operation procedure corresponding to is shown in order.
This operation is repeated every time the two-dimensional scanning surface is scanned once. 6, the same steps as those in FIG. 2 are designated by the same reference numerals, and the description thereof will be omitted.

【0025】ステップ(S25)では、以下の(2)に
したがって内部メモリに格納されている関心領域内の全
ての点の輝度信号の平均値と前回または所定時間前の走
査による輝度信号の平均値との差dIが求められ、この
dIが今回から前回または所定時間前の走査時刻までの
所定の短時間間隔dtで割り算されることにより、その
変化率△Iが計算される。
In step (S25), the average value of the luminance signals of all the points in the region of interest stored in the internal memory according to the following (2) and the average value of the luminance signals obtained by the previous scanning or a predetermined time before scanning. The difference dI is calculated and the change rate ΔI is calculated by dividing this dI by a predetermined short time interval dt from the present time to the scanning time of the previous time or a predetermined time.

【0026】 △I=dI/dt ・・・(2) この変化率はステップ(S4)と同じステップ(S2
6)、ステップ(S5)、ステップ(S6)を経過し
て、所定の短時間間隔における変化率の大きさに応じた
長さの線分で、その線分の上点が今回の平均値に一致す
るグラフ情報として、図7に示すように表示手段5に表
示される。勿論、変化率の時間変化も図3(a)と同様
にBモード画像と共に同時表示される。
ΔI = dI / dt (2) This rate of change is the same step (S2) as step (S4)
6), step (S5), and step (S6) have passed, and the line segment has a length corresponding to the magnitude of the change rate in a predetermined short time interval, and the upper point of the line segment is the average value of this time. The matching graph information is displayed on the display means 5 as shown in FIG. Of course, the change over time in the rate of change is also displayed simultaneously with the B-mode image, as in FIG.

【0027】また、信号レベル時間変化抽出器6は、1
心拍期間のグラフから輝度レベルが最大値を示す時刻か
らその最大値の半値になる時刻までの半値時間幅をレベ
ル変化を示す数値情報として求める。図8はこの半値時
間幅を説明する図であり、1心拍期間のグラフを示して
いる。信号レベル時間変化抽出器6は、1心拍期間が経
過する毎に、そのグラフの中の輝度レベルが最大値を示
す時刻t1 と、その最大値の半値になる時刻t2 を求
め、時刻t1 からt2 に至る半値時間幅tw を数値情報
として計算する。このように求めた半値時間幅tw は、
走査変換手段3を介して、Bモード画像及びグラフ情報
と共に表示手段5に表示される。
Further, the signal level time change extractor 6 is set to 1
From the graph of the heartbeat period, the half value time width from the time when the brightness level shows the maximum value to the time when the brightness level reaches the half value is obtained as the numerical information showing the level change. FIG. 8 is a diagram for explaining this half-value time width and shows a graph of one heartbeat period. The signal level time change extractor 6 obtains the time t1 at which the brightness level in the graph shows the maximum value and the time t2 at which the maximum value becomes half the value every time one heartbeat period elapses, and from the time t1 to t2. The half-value time width tw leading to is calculated as numerical information. The half-value time width tw thus obtained is
It is displayed on the display means 5 together with the B-mode image and the graph information via the scan conversion means 3.

【0028】このように本実施例によれば、Bモード画
像と同じ輝度信号を用いて、関心領域の輝度信号の輝度
レベルの時間的変化がグラフ情報として、また半値時間
幅が数値情報としてBモード画像と同じ受信信号から求
められ、Bモード画像と共にリアルタイムで表示され
る。したがって、超音波反射特性が良好なコントラスト
剤を血管に注入した場合、コントラスト剤が関心領域に
流入するに応じて輝度レベルが序々に高くなり、コント
ラスト剤が関心領域から流出するに応じて輝度レベルが
序々に低くなる様子を、リアルタイムで且つ定量的に診
断することができる。
As described above, according to this embodiment, using the same luminance signal as that of the B-mode image, the temporal change of the luminance level of the luminance signal of the region of interest is used as graph information and the half-value time width is used as numerical information. It is obtained from the same received signal as the mode image and displayed in real time together with the B mode image. Therefore, when a contrast agent having good ultrasonic reflection characteristics is injected into a blood vessel, the brightness level gradually increases as the contrast agent flows into the region of interest, and the brightness level increases as the contrast agent flows out of the region of interest. It can be diagnosed in real time and quantitatively how the value gradually decreases.

【0029】次に第2実施例について説明する。本実施
例は、輝度信号の時間的変化をBモード画像の各点毎に
計算し、この計算した結果(これを以下「インデック
ス」という)をそのレベルに応じてカラー情報に変換
し、カラーの2次元画像としてBモード画像に重ねて表
示することを特徴とする。インデックスを得る手段は前
述の通り、複数画素の平均をその中心画素のインデック
スとすることも可能である。
Next, a second embodiment will be described. In this embodiment, the temporal change of the luminance signal is calculated for each point of the B-mode image, and the calculation result (hereinafter referred to as “index”) is converted into color information according to the level, and the color information is calculated. It is characterized in that it is displayed as a two-dimensional image so as to be superimposed on the B-mode image. As described above, the means for obtaining the index can use the average of a plurality of pixels as the index of the central pixel.

【0030】図9は本実施例のブロック図を示す。な
お、図9において図1と同じ部分は同符号を付して説明
は省略する。検波器26からの輝度信号はディジタル信
号としてインデックス演算器10に供給される。インデ
ックス演算器10では、輝度レベルと前回または所定時
間前の走査による輝度信号の輝度レベルとの差が2次元
走査面の全ての点について各々求められ、この差が今回
から前回または所定時間前の走査時刻までの所定の短時
間間隔で割り算されることにより、その変化率が2次元
走査面の全ての点について各々計算される。
FIG. 9 shows a block diagram of this embodiment. In FIG. 9, the same parts as those in FIG. The brightness signal from the detector 26 is supplied to the index calculator 10 as a digital signal. The index calculator 10 obtains the difference between the luminance level and the luminance level of the luminance signal obtained by the previous scanning or the scanning before the predetermined time at each point on the two-dimensional scanning plane, and this difference is from the current time to the previous or the predetermined time before. The rate of change is calculated for each of all points on the two-dimensional scanning plane by being divided by a predetermined short time interval until the scanning time.

【0031】このインデックス演算器10で演算された
2次元走査面内の各点の変化率は、走査変換手段3と同
じ走査変換手段32に送られ、一旦記憶される。そし
て、走査変換手段3の読出しと同期して走査変換手段3
2から読み出され、カラー情報変換手段12に送られ、
そこで各点毎に又は微小領域毎にその変化率(微小領域
内の各点の変化率合計)の相違に応じて異なる色彩又は
彩度のカラー情報に変換された後、マルチプレクサ(M
PX)7に送られる。マルチプレクサ7には、走査変換
手段3の出力も供給されている。図示しない表示選択入
力手段の操作による変化率の2次元カラー画像とBモー
ド画像を重ねて表示するモードと、2次元カラー画像ま
たはBモード画像だけを表示するモードのいずれかのモ
ードにしたがって、変化率の2次元カラー画像とBモー
ド画像とが1枚の画像に合成されて出力されるか、また
は変化率の2次元カラー画像またはBモード画像だけが
出力され、D/A変換手段4を介して表示手段5に表示
される。図10は変化率の2次元カラー画像とBモード
画像との合成画像の表示例を示している。
The rate of change of each point in the two-dimensional scanning plane calculated by the index calculator 10 is sent to the same scan conversion means 32 as the scan conversion means 3 and temporarily stored therein. Then, the scanning conversion means 3 is synchronized with the reading of the scanning conversion means 3.
2 is read out and sent to the color information conversion means 12,
Therefore, after being converted into color information of different colors or saturations depending on the difference in the rate of change (sum of rate of change of each point in the minute area) for each point or for each minute area, the multiplexer (M
PX) 7. The output of the scan conversion means 3 is also supplied to the multiplexer 7. Change in the rate of change due to operation of a display selection / input means (not shown) according to one of a mode for displaying a two-dimensional color image and a B-mode image in an overlapping manner and a mode for displaying only a two-dimensional color image or a B-mode image Rate two-dimensional color image and B-mode image are combined and output as one image, or only the two-dimensional color image or B-mode image with a change rate is output and is output via the D / A conversion means 4. Is displayed on the display means 5. FIG. 10 shows a display example of a composite image of a two-dimensional color image with a change rate and a B-mode image.

【0032】このように本実施例によれば、Bモード画
像と同じ輝度信号を用いて、2次元走査面全域の輝度レ
ベルの時間的変化の空間的な様子が色彩又は彩度の相違
としてBモード画像と共にリアルタイムで表示される。
したがって、超音波反射特性が良好なコントラスト剤を
血管に注入した場合、コントラスト剤がある部分に流入
してから流出する様子をリアルタイムで且つ空間的に診
断することができる。インデックスをそのレベルに応じ
た輝度信号として画像化しても同様の効果が得られるこ
とは勿論である。
As described above, according to the present embodiment, the same luminance signal as that of the B-mode image is used, and the spatial state of the temporal change of the luminance level over the entire two-dimensional scanning plane is the difference in color or saturation as B. It is displayed in real time with the mode image.
Therefore, when a contrast agent having a good ultrasonic reflection characteristic is injected into a blood vessel, it is possible to make a real-time and spatial diagnosis of how a contrast agent flows into a portion where it exists and then flows out. Of course, the same effect can be obtained by imaging the index as a luminance signal corresponding to the level.

【0033】本発明による関心領域の輝度レベルを得る
とき、操作者が超音波プローブを誤って僅かに動かして
しまった場合や、被検体の呼吸や心拍動による臓器の移
動がある場合、Bモード画像又はそれに代わる2次元又
は3次元情報に基づいてパターンマッチング技法等によ
り、関心領域の動きの補正をすることが考えられる。こ
れにより精度向上を図ることができる。
When obtaining the brightness level of the region of interest according to the present invention, when the operator accidentally slightly moves the ultrasonic probe, or when there is movement of the organ due to respiration or heartbeat of the subject, the B mode is used. It is conceivable to correct the movement of the region of interest by a pattern matching technique or the like based on the image or two-dimensional or three-dimensional information that substitutes for it. This can improve accuracy.

【0034】本発明は上述した実施例に限定されず、種
々変形して実施可能である。例えば上述の説明では、輝
度レベルの時間的変化をグラフや数値情報として表示す
ることと説明したが、一般的な超音波ドプラ装置で計算
される周波数シフト(偏移周波数)のパワーを輝度レベ
ルに置換して、このパワーの時間的変化を求めるように
してもよい。また上述の説明では、検波結果としての輝
度信号を用いて時間的変化を求めているが、検波器に入
力される加算器の出力信号を用いて、その振幅値の時間
的変化を求めるようにしてもよい。
The present invention is not limited to the above-mentioned embodiments, but can be implemented with various modifications. For example, in the above description, it has been explained that the temporal change of the brightness level is displayed as a graph or numerical information, but the power of frequency shift (shift frequency) calculated by a general ultrasonic Doppler device is set to the brightness level. Alternatively, the change in power may be calculated with time. Further, in the above description, the temporal change is obtained by using the luminance signal as the detection result, but the temporal change of the amplitude value is obtained by using the output signal of the adder input to the detector. May be.

【0035】[0035]

【発明の効果】本発明による被検体内を超音波ビームで
走査することにより得られた受信信号の包絡線を検波
し、この検波結果を輝度信号として断層像を表示する超
音波診断装置は、前記輝度信号を入力し、前記輝度信号
の時間的変化を計算する手段と、前記時間的変化を前記
断層像と共に表示する手段とを具備するようにしたの
で、断層像としての輝度信号を取り込んでその時間的変
化を計算して、この時間的変化を断層像と共にリアルタ
イムで表示することができ、これによりコントラスト剤
を血管に注入した場合には血流状態を短時間のうちに且
つ良好に診断することができる。
The ultrasonic diagnostic apparatus for detecting the envelope of the received signal obtained by scanning the inside of the subject with the ultrasonic beam according to the present invention and displaying a tomographic image using the detection result as a luminance signal is The luminance signal is input, the means for calculating the temporal change of the luminance signal and the means for displaying the temporal change together with the tomographic image are provided, so that the luminance signal as a tomographic image is captured. The temporal change can be calculated, and this temporal change can be displayed in real time together with a tomographic image. With this, when a contrast agent is injected into a blood vessel, the blood flow state can be diagnosed well in a short time. can do.

【0036】また、本発明による被検体内を超音波ビー
ムで走査することにより得られた受信信号の包絡線を検
波し、この検波結果を輝度信号として断層像を表示する
他の超音波診断装置は、前記輝度信号を入力し、前記検
波信号の時間的変化を前記断層像の各点毎に計算する手
段と、前記各時間的変化をそのレベルに応じてカラー情
報に変換し2次元画像として前記断層像に重ねて表示す
る手段とを具備するようにしたので、輝度信号を取り込
んでその時間的変化を各点毎に計算して、各時間的変化
をそのレベルに応じてカラー情報に変換し2次元画像と
して断層像に重ねてリアルタイムで表示することがで
き、これによりコントラスト剤を血管に注入した場合に
は血流状態を短時間のうちに診断することができ、また
血流の空間的な様子も把握することができる。
Further, another ultrasonic diagnostic apparatus according to the present invention detects the envelope of the received signal obtained by scanning the inside of the subject with an ultrasonic beam and displays a tomographic image using the detection result as a luminance signal. Means for inputting the luminance signal, calculating time change of the detection signal for each point of the tomographic image, and converting each time change into color information according to the level thereof to obtain a two-dimensional image. Since a means for displaying the tomographic image in an overlapping manner is provided, the luminance signal is taken in, the temporal change thereof is calculated for each point, and the temporal change is converted into color information according to the level. As a two-dimensional image, it can be displayed in real time superimposed on a tomographic image. This allows the blood flow state to be diagnosed in a short time when a contrast agent is injected into a blood vessel, and the blood flow space State It is possible to grasp.

【図面の簡単な説明】[Brief description of drawings]

【図1】第1実施例のブロック図。FIG. 1 is a block diagram of a first embodiment.

【図2】関心領域内全点の平均値の時間的変化を求める
場合の図1の信号レベル時間変化抽出器による2次元走
査面の1回の走査動作に対応する動作手順を示す流れ
図。
FIG. 2 is a flowchart showing an operation procedure corresponding to one scanning operation of a two-dimensional scanning plane by the signal level time change extractor of FIG. 1 when obtaining a time change of an average value of all points in a region of interest.

【図3】図2の場合の表示手段の表示画面を示す図。FIG. 3 is a diagram showing a display screen of a display unit in the case of FIG.

【図4】関心領域内の最大値及び/又は最小値の時間的
変化を求める場合の図1の信号レベル時間変化抽出器に
よる2次元走査面の1回の走査動作に対応する動作手順
を示す流れ図。
FIG. 4 shows an operation procedure corresponding to one scanning operation of a two-dimensional scanning plane by the signal level time change extractor of FIG. 1 when obtaining a time change of a maximum value and / or a minimum value in a region of interest. flow diagram.

【図5】図4の場合の表示手段の表示画面を示す図。FIG. 5 is a diagram showing a display screen of a display unit in the case of FIG.

【図6】関心領域内全点の平均値の所定の短時間間隔に
おける変化率の時間的変化を求める場合の図1の信号レ
ベル時間変化抽出器による2次元走査面の1回の走査動
作に対応する動作手順を示す流れ図。
FIG. 6 shows a single scanning operation of a two-dimensional scanning plane by the signal level time change extractor of FIG. 1 in the case of obtaining the time change of the change rate of the average value of all the points in the region of interest in a predetermined short time interval. The flowchart which shows the corresponding operation procedure.

【図7】図6の場合の表示手段の表示画面を示す図。7 is a diagram showing a display screen of the display means in the case of FIG.

【図8】図1の信号レベル時間変化抽出器で計算される
半値時間幅を説明する図。
8 is a diagram illustrating a half value time width calculated by the signal level time change extractor of FIG. 1;

【図9】第2実施例のブロック図。FIG. 9 is a block diagram of a second embodiment.

【図10】第2実施例の表示画面を示す図。FIG. 10 is a diagram showing a display screen of the second embodiment.

【符号の説明】[Explanation of symbols]

1…超音波プローブ、2…アナログ部、3…走査変換手
段、4…D/A変換手段、5…表示手段、6…信号レベ
ル時間変化抽出器、21…プリアンプ、22…パルサ、
24…ディレーライン、25…加算器、26…検波器。
DESCRIPTION OF SYMBOLS 1 ... Ultrasonic probe, 2 ... Analog part, 3 ... Scan conversion means, 4 ... D / A conversion means, 5 ... Display means, 6 ... Signal level time change extractor, 21 ... Preamplifier, 22 ... Pulser,
24 ... delay line, 25 ... adder, 26 ... detector.

───────────────────────────────────────────────────── フロントページの続き (72)発明者 内海 勲 栃木県大田原市下石上1385番の1 株式会 社東芝那須工場内 (72)発明者 上木 康至 栃木県大田原市下石上1385番の1 株式会 社東芝那須工場内 ─────────────────────────────────────────────────── ─── Continuation of the front page (72) Inventor Isao Utsumi 1385-1 Shimoishigami, Otawara-shi, Tochigi Stock company Toshiba Nasu factory (72) Inventor Yasushi Ueki 1385-1,385 Shimoishi, Otawara, Tochigi Stock company Toshiba Nasu factory

Claims (2)

【特許請求の範囲】[Claims] 【請求項1】 被検体内を超音波ビームで走査すること
により得られた受信信号の包絡線を検波し、この検波結
果を輝度信号として断層像を表示する超音波診断装置に
おいて、 前記輝度信号を入力し、前記輝度信号の時間的変化を計
算する手段と、 前記時間的変化を前記断層像と共に表示する手段とを具
備することを特徴とする超音波診断装置。
1. An ultrasonic diagnostic apparatus for detecting an envelope of a received signal obtained by scanning the inside of a subject with an ultrasonic beam, and displaying a tomographic image using the detection result as a luminance signal. The ultrasonic diagnostic apparatus comprises: means for calculating the temporal change of the luminance signal; and means for displaying the temporal change together with the tomographic image.
【請求項2】 被検体内を超音波ビームで走査すること
により得られた受信信号の包絡線を検波し、この検波結
果を輝度信号として断層像を表示する超音波診断装置に
おいて、 前記輝度信号を入力し、前記検波信号の時間的変化を前
記断層像の各点毎に計算する手段と、 前記各時間的変化をそのレベルに応じてカラー情報に変
換し2次元画像として前記断層像に重ねて表示する手段
とを具備することを特徴とする超音波診断装置。
2. An ultrasonic diagnostic apparatus for detecting an envelope of a received signal obtained by scanning the inside of a subject with an ultrasonic beam, and displaying a tomographic image using the detection result as a luminance signal. Means for calculating the temporal change of the detection signal for each point of the tomographic image, and converting the temporal change into color information according to the level thereof and superimposing it on the tomographic image as a two-dimensional image. An ultrasonic diagnostic apparatus, comprising:
JP32348493A 1993-12-22 1993-12-22 Ultrasound diagnostic equipment Expired - Lifetime JP3405578B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP32348493A JP3405578B2 (en) 1993-12-22 1993-12-22 Ultrasound diagnostic equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP32348493A JP3405578B2 (en) 1993-12-22 1993-12-22 Ultrasound diagnostic equipment

Publications (2)

Publication Number Publication Date
JPH07178095A true JPH07178095A (en) 1995-07-18
JP3405578B2 JP3405578B2 (en) 2003-05-12

Family

ID=18155210

Family Applications (1)

Application Number Title Priority Date Filing Date
JP32348493A Expired - Lifetime JP3405578B2 (en) 1993-12-22 1993-12-22 Ultrasound diagnostic equipment

Country Status (1)

Country Link
JP (1) JP3405578B2 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006110028A (en) * 2004-10-14 2006-04-27 Ge Medical Systems Global Technology Co Llc Ultrasonic imaging apparatus
JP2006181058A (en) * 2004-12-27 2006-07-13 Toshiba Corp Ultrasonic diagnostic equipment and display method of ultrasonic diagnostic image
JP2009028194A (en) * 2007-07-26 2009-02-12 Ge Medical Systems Global Technology Co Llc Ultrasonic imaging system
WO2014142037A1 (en) * 2013-03-11 2014-09-18 株式会社東芝 Ultrasound diagnostic apparatus and image processing device

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01254151A (en) * 1988-04-04 1989-10-11 Fujitsu Ltd Ultrasonic diagnosing device
JPH04117950A (en) * 1990-09-07 1992-04-17 Ken Ishihara Ultrasonic diagnostic device
JPH04208143A (en) * 1990-10-12 1992-07-29 Hitachi Medical Corp Ultrasonic diagnostic apparatus
US5235984A (en) * 1992-03-30 1993-08-17 Hewlett-Packard Company On-line acoustic densitometry tool for use with an ultrasonic imaging system
JPH0759781A (en) * 1993-06-14 1995-03-07 Toshiba Corp Ultrasonic diagnostic system
JPH0775638A (en) * 1993-07-12 1995-03-20 Toshiba Corp Ultrasonic diagnostic device
JPH0779974A (en) * 1993-09-13 1995-03-28 Toshiba Corp Ultrasonic diagnostic apparatus

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01254151A (en) * 1988-04-04 1989-10-11 Fujitsu Ltd Ultrasonic diagnosing device
JPH04117950A (en) * 1990-09-07 1992-04-17 Ken Ishihara Ultrasonic diagnostic device
JPH04208143A (en) * 1990-10-12 1992-07-29 Hitachi Medical Corp Ultrasonic diagnostic apparatus
US5235984A (en) * 1992-03-30 1993-08-17 Hewlett-Packard Company On-line acoustic densitometry tool for use with an ultrasonic imaging system
JPH0759781A (en) * 1993-06-14 1995-03-07 Toshiba Corp Ultrasonic diagnostic system
JPH0775638A (en) * 1993-07-12 1995-03-20 Toshiba Corp Ultrasonic diagnostic device
JPH0779974A (en) * 1993-09-13 1995-03-28 Toshiba Corp Ultrasonic diagnostic apparatus

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006110028A (en) * 2004-10-14 2006-04-27 Ge Medical Systems Global Technology Co Llc Ultrasonic imaging apparatus
JP2006181058A (en) * 2004-12-27 2006-07-13 Toshiba Corp Ultrasonic diagnostic equipment and display method of ultrasonic diagnostic image
JP4679141B2 (en) * 2004-12-27 2011-04-27 株式会社東芝 Ultrasonic diagnostic apparatus and method for displaying ultrasonic diagnostic image
JP2009028194A (en) * 2007-07-26 2009-02-12 Ge Medical Systems Global Technology Co Llc Ultrasonic imaging system
WO2014142037A1 (en) * 2013-03-11 2014-09-18 株式会社東芝 Ultrasound diagnostic apparatus and image processing device
JP2014171755A (en) * 2013-03-11 2014-09-22 Toshiba Corp Ultrasound diagnostic device, image processing system and image processing program
CN105188555A (en) * 2013-03-11 2015-12-23 株式会社东芝 Ultrasound diagnostic apparatus and image processing device
CN105188555B (en) * 2013-03-11 2018-06-22 东芝医疗系统株式会社 Diagnostic ultrasound equipment and image processing apparatus
US11298109B2 (en) 2013-03-11 2022-04-12 Canon Medical Systems Corporation Ultrasonic diagnostic apparatus and image processing apparatus

Also Published As

Publication number Publication date
JP3405578B2 (en) 2003-05-12

Similar Documents

Publication Publication Date Title
US9241689B2 (en) Ultrasonic diagnostic equipment and imaging processing apparatus
JP4786150B2 (en) Ultrasonic diagnostic apparatus and image processing apparatus
US6824514B2 (en) System and method for visualizing scene shift in ultrasound scan sequence
JP3043873B2 (en) Ultrasonic aperture synthesis system
US7758507B2 (en) Blood flow imaging
KR100748178B1 (en) Ultrasound diagnostic system and method for displaying arbitrary m-mode images
JP5134787B2 (en) Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, and ultrasonic image processing program
JP2007513726A (en) Ultrasound imaging system with automatic control of penetration, resolution and frame rate
WO2007114375A1 (en) Ultrasound diagnostic device and control method for ultrasound diagnostic device
JPH08173422A (en) Ultrasonic diagnostic system
JP7313902B2 (en) ultrasound diagnostic equipment
JP2010017537A (en) Ultrasonic diagnostic apparatus
JPH11327A (en) Ultrasonograph
JP4764209B2 (en) Ultrasonic signal analysis apparatus, ultrasonic signal analysis method, ultrasonic analysis program, ultrasonic diagnostic apparatus, and control method of ultrasonic diagnostic apparatus
JPH10118061A (en) Ultrasonograph
US20120203111A1 (en) Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, and ultrasonic image acquisition method
US20220313220A1 (en) Ultrasound diagnostic apparatus
JP4758736B2 (en) Ultrasonic diagnostic equipment
JP5606025B2 (en) Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, and ultrasonic image processing program
JP3405578B2 (en) Ultrasound diagnostic equipment
CN113301854A (en) Image analysis device
JP2007130063A (en) Ultrasonographic apparatus
JP5823184B2 (en) Ultrasonic diagnostic apparatus, medical image processing apparatus, and medical image processing program
WO2003001240A1 (en) Ultrasonic diagnostic system for selectively developing ultrasound diagnostic data
JP3047265B2 (en) Three-dimensional image creation method for ultrasonic diagnostic equipment

Legal Events

Date Code Title Description
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080307

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20090307

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20100307

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20100307

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20110307

Year of fee payment: 8

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

Free format text: PAYMENT UNTIL: 20120307

Year of fee payment: 9

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

Free format text: PAYMENT UNTIL: 20130307

Year of fee payment: 10

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

Free format text: PAYMENT UNTIL: 20140307

Year of fee payment: 11

EXPY Cancellation because of completion of term