JPH0775639A - 超音波診断装置 - Google Patents

超音波診断装置

Info

Publication number
JPH0775639A
JPH0775639A JP22318093A JP22318093A JPH0775639A JP H0775639 A JPH0775639 A JP H0775639A JP 22318093 A JP22318093 A JP 22318093A JP 22318093 A JP22318093 A JP 22318093A JP H0775639 A JPH0775639 A JP H0775639A
Authority
JP
Japan
Prior art keywords
frequency band
circuit
frequency
blood flow
filter
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
JP22318093A
Other languages
English (en)
Other versions
JP3471860B2 (ja
Inventor
Takatoshi Okumura
貴敏 奥村
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
Original Assignee
Toshiba 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 Toshiba Corp filed Critical Toshiba Corp
Priority to JP22318093A priority Critical patent/JP3471860B2/ja
Publication of JPH0775639A publication Critical patent/JPH0775639A/ja
Application granted granted Critical
Publication of JP3471860B2 publication Critical patent/JP3471860B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

(57)【要約】 【目的】本発明は、フィルタ特性を維持しつつ折り返し
現象の発生を解消し、それに伴って結果物たる血流情報
のS/N比を向上させ得る超音波診断装置を提供するこ
とを目的とする。 【構成】本発明に係る超音波診断装置は、超音波を被検
体内に送波し、反射体からの反射波を受波し、受信信号
を得る手段と、前記受信信号について所定のフィルタ通
過周波数帯域に合致する成分を通過させるフィルタ回路
6と、前記フィルタ通過周波数帯域と略同一のサンプル
周波数帯域を有し、フィルタ回路6の出力をアナログ/
ディジタル変換するアナログ/ディジタル変換回路7
と、前記アナログ/ディジタル変換回路7の出力を周波
数分析することにより前記送波した超音波の血流による
反射波のドプラ偏移周波数を検出して血流速を求める周
波数分析回路8と、その血流速を表示する表示系12と
を具備することを特徴とする。

Description

【発明の詳細な説明】
【0001】
【産業上の利用分野】本発明は、超音波を生体内に送波
し該超音波の血球からの反射波のドプラ偏移周波数を検
出して血流速度を表示する超音波診断装置に関する。
【0002】
【従来の技術】超音波診断装置におけるドプラ血流検出
回路は、被検体内に超音波を送波したときに、超音波反
射体が移動体である場合に、その移動体からの反射波の
周波数が、送波時の周波数に対して偏移するドプラ現象
を利用して、血流情報を検出しようとするものである。
周波数の偏移量、即ち偏移周波数は、血流速度の超音波
送受波方向成分に対応するものであって、ドプラ血流検
出回路はこの偏移周波数(血流速)を検出するために、
反射波を受信して得たエコー信号をサイドローブ等の等
の緩動反射体からの不要信号成分(ノイズ成分)を除去
するLPF(ロー・パス・フィルタ)等のフィルタ回路
を介して後、A/D変換器によりディジタル信号に変換
し、そしてディジタル信号処理系である高速フーリエ変
換(周波数解析)処理を行う。なお、高速フーリエ変換
処理の特性上、解析可能な周波数帯域△fc は、A/D
変換器におけるサンプリング周波数fsampに依存し、そ
の関係は次式(1)に示す通りである。
【0003】 −fsamp/2 < △fc < +fsamp/2 −(1) 周波数解析回路(FFT)に入力した信号の内、周波数
帯域△fc の範囲外の周波数成分は、いわゆる折返し現
象を生じさせる。
【0004】ここで、A/D変換器のサンプリング周波
数帯域△fsampと、フィルタ回路を通過するフィルタ通
過周波数帯域△fLPF と、偏移周波数(血流速)を表示
する表示周波数帯域△fdispとは、一般に次式(2)に
示す関係で設定されている。なお、記号「…」は略同一
を意味する。
【0005】 △fdisp…△fsamp…△fLPF /2 −(2) 図7は、これらの関係および折返し現象が生じたFFT
からの出力波形(表示波形)について示した図であり、
図8は折返し現象が生じやすい場合の一例として心臓の
償帽弁逆流について示した図である。図8に示したよう
に、償帽弁MVにおいて逆流MRが発生するような場
合、偏移周波数(血流速)の最大変化幅は増大し、その
一部(図7の点線部分T2')がサンプリング周波数帯域
△fsampおよび表示周波数帯域△fdispを越えることと
なり、該部T2'は、折返し現象により表示周波数帯域△
fdispのゼロ位置P0 を挾んで逆側に現れてしまう。
【0006】また、上記式(2)の関係によると、フィ
ルタ通過周波数帯域△fLPF の内、サンプリング周波数
帯域△fsampに含まれない部分、いわゆる折返し部分N
1 (斜線部分)は、ノイズ成分に重畳され、結果的に得
られる血流情報のS/N比が低下してしまう。
【0007】なお、式(2)の関係を維持しつつ、上記
S/N比を向上させるには、LPF周波数帯域△fLPF
を広く設定し、それに応じてサンプリング周波数帯域△
fsampおよび表示周波数帯域△fdispも拡大させるれば
よい。しかし、上述したようにこのLPF周波数帯域△
fLPF はサイドローブ等の緩動反射体からのノイズ成分
を除去するフィルタ特性を考慮して設定されるものであ
り、広く設定すればするほどノイズ成分が増加してしま
う。
【0008】
【発明が解決しようとする課題】そこで本発明の目的
は、フィルタ特性を維持しつつ折り返し現象の発生を解
消し、それに伴って結果物たる血流情報のS/N比を向
上させ得る超音波診断装置を提供することである。
【0009】
【課題を解決するための手段】本発明に係る超音波診断
装置は、超音波を被検体内に送波し、反射体からの反射
波を受波し、受信信号を得る手段と、前記受信信号につ
いて所定のフィルタ通過周波数帯域に合致する成分を通
過させるフィルタ回路と、前記フィルタ通過周波数帯域
と略同一のサンプル周波数帯域を有し、前記フィルタ回
路の出力をアナログ/ディジタル変換するアナログ/デ
ィジタル変換回路と、前記アナログ/ディジタル変換回
路の出力を周波数分析することにより前記送波した超音
波の血流による反射波のドプラ偏移周波数を検出して血
流速を求める周波数分析回路と、その血流速を表示する
表示系とを具備することを特徴とする。
【0010】
【作用】本発明によれば、超音波を被検体内に送波し、
反射体からの反射波を受波し、受信信号を得、その受信
信号についてフィルタ回路においてフィルタ通過周波数
帯域に合致する成分だけを通過させ、前記フィルタ回路
の出力をアナログ/ディジタル変換回路で前記フィルタ
通過周波数帯域と略同一のサンプル周波数帯域を用いて
アナログ/ディジタル変換することによって、ノイズ成
分の増加を抑えることができる。
【0011】
【実施例】以下、図面を参照しながら実施例を説明す
る。図1は本発明の一実施例に係る超音波診断装置の構
成を示すブロック図であり、図2は図1に示したfmax
測定装置の一例を示した図である。
【0012】超音波診断装置は図示しないシステムコン
トローラをシステム全体の制御中枢として、超音波探触
子1、送波回路2、受波回路3、直交検波回路(90度D
D)4、バンドパスフィルタ(BPF)6、アナログ/
ディジタル変換器(ADC)7、高速フーリエ変換器
(FFT)8、fmax 測定回路9、ゼロシフト量設定回
路10、ディジタル・スキャン・コンバータ(DSC)
11および画像表示装置12を備えている。なお、本実
施例装置の構成は、説明の便宜上、本発明が主要に関わ
る連続波ドプラモード(CWDモード)処理系に関する
部分のみについて示し、この他の超音波断層像を得るB
モード処理系等については省略した。
【0013】超音波探触子1は、複数の圧電振動子を並
設してなり、これらの振動子により被検体に対して超音
波パルスを繰返し送信および受信する。送波回路2は、
図示しないが、クロック発生器、送信遅延回路、連続波
駆動回路を備えていて、各振動子を駆動する駆動信号を
送信側の超音波探触子1の各振動子に供給し、超音波探
触子1から所定の単一周波数で超音波パルスを送波させ
るものである。このような送波回路2により超音波探触
子1が送信駆動されると、超音波探触子1から被検体内
に超音波パルスが送波され、被検体で流動する血球で反
射し、受信側の超音波探触子の同一振動子で受波され、
ドプラ偏移をともなう受信信号に変換される。勿論、受
信側の超音波探触子は送信側の超音波探触子1と兼用す
ることが一般的である。
【0014】この超音波の送受波で得た受信信号を受け
る受波回路3は、図示しないが、プリアンプ、受信遅延
回路、および加算器を備えていて、受信信号を所定のレ
ベルまで増幅し、増幅された受信信号に前記送信遅延回
路で与えた遅延時間を基に戻すような遅延時間を各振動
子毎に与え、受信信号を加算合成し、出力する。
【0015】直交検波回路(90度DD)4は、ドプラ復
調回路の1つであって、パルスドプラ法には多くのドプ
ラ信号が存在するため、参照波周波数によって所望の送
信スペクトルに対するドプラ信号を抽出するものであ
る。なお、このドプラ復調回路は、直交検波であるの
で、受信信号を2チャンネルに分け、互いに90度位相
の異なる参照波と掛け合わせ(直交検波)、そして不要
な高周波成分を取り除くものである。
【0016】バンドパスフィルタ5は、所定のフィルタ
通過周波数帯域△fLPF を有するものであり、まず該出
力を低域通過フィルタにより波形形成し、さらに低周波
成分除去フィルタにより復調されたドプラ信号中の心臓
壁や弁運動の緩動体からの低周波成分を除去する。
【0017】ADC7は、バンドパスフィルタ6までの
アナログ信号処理系と、後流回路のディジタル信号処理
系である高速フーリエ変換器8およびDSC11等とを
つなぐため、アナログ信号をディジタル信号に変換する
ものである。なお、このADC7のサンプリング周波数
帯域△fsampは、バンドパスフィルタ6のフィルタ通過
周波数帯域△fLPF に対し次の式(3)の関係で設定さ
れる。なお、記号「…」は略同一を意味する。
【0018】 △fsamp…△fLPF −(3) 図3は、このフィルタ通過周波数帯域△fLPF とサンプ
リング周波数帯域△fsampと後述する画像表示装置12
の表示周波数帯域△fdispとの関係および高速フーリエ
変換器からの出力波形WFFT について周波数空間上で示
した図である。
【0019】このようにサンプリング周波数帯域△fsa
mpが、フィルタ通過周波数帯域△fLPF と略同一の関係
で設定されていることによって、サンプリング周波数帯
域△fsamp以外のノイズ成分は、既にバンドパスフィル
タ6で減衰されているために、折返しノイズは解消され
ることとなる。これにより、従来の場合に比し、S/N
比が約3dB向上する。
【0020】高速フーリエ変換器8は、ADC7の出力
に対し高速フーリエ変換処理することにより、反射波の
偏移周波数fn を得る。この偏移周波数fn は、血球の
流動状態に応じて刻々と変化するものであって、従来の
場合と比し、フィルタ通過周波数帯域△fLPF と略同一
程度にまでサンプリング周波数帯域△fsampが拡大され
て設定されていることにより、高速フーリエ変換器8の
処理範囲は、前記偏移周波数fn の変化の幅△fFFT を
十分捕らえることができ、従来のように折返しの無い出
力波形WFFT を得ることができる。
【0021】ここで、フィルタ通過周波数帯域△fLPF
とサンプリング周波数帯域△fsampとを略同一に設定し
たことにより、高速フーリエ変換器8からの出力波形W
FFTには折返し現象が生じないこととなったが、図3に
示すように、一般にそれらの周波数帯域より狭く設定さ
れる画像表示装置12における血流速を表示する表示周
波数帯域△fdispを越える周波数成分(斜線部分)が発
生することとなり、その部分については、図4に示すよ
うに、画像表示装置12で折返し表示(点線)されるこ
とはないが、一点鎖線で示した該部分の波形が表示され
なくなってしまうこととなる。なお、それらの帯域の関
係は、前記式(3)の関係を踏まえ、式(4)に示す関
係で設定される。
【0022】 △fdisp…△fsamp/2…△fLPF /2 −(4) そこで、fmax 測定回路9およびゼロシフト量設定回路
10で、前記偏移周波数fn の変化幅△fFFT をこの表
示周波数帯域△fdispに収めるように、図5に示したよ
うに、ゼロ位置P0 をP0'にずらすためのシフト量Sf
を自動測定する。
【0023】fmax 測定回路9は、高速フーリエ変換器
8から出力される偏移周波数fn の中の最大値fmax を
測定するものである。なお、この測定方法は、種々考え
られるが、ここでは、図2に示したように、血流速を可
聴音として処理するオーディオ回路に可変ハイパスフィ
ルタ91 と、この可変ハイパスフィルタ91 の出力をノ
イズと分離する比較器92 とを用いて、可変ハイパスフ
ィルタ91 の初期設定されている遮断レベルを越える信
号の有無を調べる方法を採用する。なお、可変ハイパス
フィルタ91 の遮断レベルは、比較器92 からのカット
オフ(CF)信号により変化する。これによって、連続
的に供給される偏移周波数fn の中の最大値fmax を測
定することができる。
【0024】ゼロシフト量設定回路10は、最大値fma
x および表示周波数帯域△fdispに基づいて、前記偏移
周波数fn の変化幅△fFFT をこの表示周波数帯域△f
dispに収めるように、ゼロ位置P0 をP0'にずらすため
のシフト量Sfを自動測定する。すなわち、最大値fma
x と、表示周波数帯域△fdispの最高位レベルとを比較
し、それらが同程度となるようにシフト量Sfを設定す
るものである。
【0025】このシフト量Sfは、DSC11に供給さ
れ、DSC11はこのシフト量Sfに応じてゼロ位置P
0 ををP0'に変更し、画像表示装置12のモニタ走査方
式に応じてCWDモード画像を出力する。画像表示装置
12は、そのCWDモード画像をアナログ信号に変換し
た後、表示する。
【0026】以上のように、バンドパスフィルタ6のフ
ィルタ通過周波数帯域△fLPF と、ADC7のサンプリ
ング周波数帯域△fsamp、画像表示装置12の血流速を
表示する表示周波数帯域△fdispとを、式(4)の関係
で設定することとしたことによって、高速フーリエ変換
器8で周波数解析した後の出力には、従来の場合のよう
に折返し現象が生じることはなく、さらに、その出力、
即ち偏移周波数fn の変化幅△fFFT をこの表示周波数
帯域△fdispに収めるようにゼロシフト量Sfを自動設
定することによって、偏移周波数fn の変化、即ち血流
速変化を全て表示することができる。
【0027】なお本発明は上記実施例に限定されること
なく種々変形して実施例可能である。例えば、上述した
ように偏移周波数fn の最大値fmax の測定装置は通常
用いられている他の装置を用いてもよい。
【0028】また、上記実施例では、偏移周波数fn の
最大値fmax のみを測定してゼロシフト量Sfを設定し
ていたが、偏移周波数fn の変化幅△fFFT が十分広
く、表示周波数帯域△fdispに収まらない場合には、こ
のゼロシフトによる対応だけでは全ての血流速を表示し
きれない場合が生じる。図6の斜線で示した部分は、こ
の場合の表示されない波形部分である。この場合には、
偏移周波数fn の最大値fmax のみではなく最小値fmi
n をも測定し、その最大値fmax および最小値fmin に
応じて表示周波数帯域△fdispの幅を△fdisp´に自動
調整することとしてもよい。なおfmax 測定回路と逆の
フィルタ特性を有する最小値測定回路をfmax 測定回路
とFFT8に対して並列に接続し、その最小値fmin お
よび最大値fmax を入力しそれらに規定される幅と表示
周波数帯域△fdispとを比較し表示周波数帯域△fdisp
を該幅に応じて△fdisp´に変化させる回路を付設すれ
ばよい。なお、式(4)△fdisp…△fsamp/2…△f
LPF /2の関係について、予め、式(5)のように設定
しておけば、表示周波数帯域を調整する必要をなくする
ことができる。 △fdisp…△fsamp…△fLPF −(5)
【0029】
【発明の効果】以上説明したように本発明によれば、超
音波を被検体内に送波し、反射体からの反射波を受波
し、受信信号を得、その受信信号についてフィルタ回路
においてフィルタ通過周波数帯域に合致する成分だけを
通過させ、前記フィルタ回路の出力をアナログ/ディジ
タル変換回路で前記フィルタ通過周波数帯域と略同一の
サンプル周波数帯域を用いてアナログ/ディジタル変換
することによって、ノイズ成分の増加を抑えることがで
きる超音波診断装置を提供することができる。
【0030】また、前記アナログ/ディジタル変換回路
の出力を周波数分析回路で周波数分析することにより前
記送波した超音波の血流による反射波のドプラ偏移周波
数を検出して血流速を求めると共に、前記ドプラ偏移周
波数の最大値を計測し、その最大値に基づいてゼロシフ
ト量を算出し、そのゼロシフト量に応じて前記ドプラ偏
移周波数の表示範囲を規定する表示周波数帯域のゼロ位
置を転移させ、その表示周波数帯域の範囲内で前記周波
数分析回路で求めた血流速を表示することによって、検
出された前記ドプラ偏移周波数(血流速)の全範囲を表
示することができる超音波診断装置を提供することがで
きる。
【図面の簡単な説明】
【図1】本発明一実施例に係る超音波診断装置の構成を
示すブロック図。
【図2】図1に示したfmax 測定回路の構成を示すブロ
ック図。
【図3】図1に示したバンドパスフィルタのフィルタ通
過周波数帯域△fLPF とADCのサンプリング周波数帯
域△fsampと画像表示装置の表示周波数帯域△fdispと
の関係および高速フーリエ変換器からの出力波形の一例
を示す図。
【図4】図3に示した出力波形をそのまま画像表示装置
で表示した場合の表示波形を示す図。
【図5】図1に示したゼロシフト量設定回路によるゼロ
シフト量に応じてゼロ位置をシフトした場合の図3に示
した出力波形を示す図。
【図6】本実施例の変形例として掲げた表示周波数帯域
△fdispの変化の様子について示す図。
【図7】従来装置のバンドパスフィルタのフィルタ通過
周波数帯域△fLPF とADCのサンプリング周波数帯域
△fsampと画像表示装置の表示周波数帯域△fdispとの
関係および高速フーリエ変換器からの出力波形の一例を
示す図。
【図8】折返り現象の生じやすい逆流の発生しやすい部
位について示す図。
【符号の説明】
1…超音波探触子、2…送波回路、3…受波回路、4…
直交検波回路、6…バンドパスフィルタ、7…アナログ
ディジタル変換器、8…高速フーリエ変換器、9…fma
x 測定回路、10…ゼロシフト設定回路、11…ディジ
タルスキャンコンバータ、12…画像表示装置。

Claims (2)

    【特許請求の範囲】
  1. 【請求項1】 超音波を被検体内に送波し、反射体から
    の反射波を受波し、受信信号を得る手段と、 前記受信信号について所定のフィルタ通過周波数帯域に
    合致する成分を通過させるフィルタ回路と、 前記フィルタ通過周波数帯域と略同一のサンプル周波数
    帯域を有し、前記フィルタ回路の出力をアナログ/ディ
    ジタル変換するアナログ/ディジタル変換回路と、 前記アナログ/ディジタル変換回路の出力を周波数分析
    することにより前記送波した超音波の血流による反射波
    のドプラ偏移周波数を検出して血流速を求める周波数分
    析回路と、 その血流速を表示する表示系とを具備することを特徴と
    する超音波診断装置。
  2. 【請求項2】 前記周波数分析回路から出力される前記
    ドプラ偏移周波数の最大値を計測する手段と、 その最大値に基づいてゼロシフト量を算出しする手段
    と、 そのゼロシフト量に応じて前記表示系における前記ドプ
    ラ偏移周波数の表示範囲を規定する表示周波数帯域のゼ
    ロ位置を転移させ、その表示周波数帯域の範囲内で前記
    周波数分析回路で求めた血流速を表示させる手段とをさ
    らに具備することを特徴とする請求項1に記載された超
    音波診断装置。
JP22318093A 1993-09-08 1993-09-08 超音波診断装置 Expired - Fee Related JP3471860B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP22318093A JP3471860B2 (ja) 1993-09-08 1993-09-08 超音波診断装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP22318093A JP3471860B2 (ja) 1993-09-08 1993-09-08 超音波診断装置

Publications (2)

Publication Number Publication Date
JPH0775639A true JPH0775639A (ja) 1995-03-20
JP3471860B2 JP3471860B2 (ja) 2003-12-02

Family

ID=16794064

Family Applications (1)

Application Number Title Priority Date Filing Date
JP22318093A Expired - Fee Related JP3471860B2 (ja) 1993-09-08 1993-09-08 超音波診断装置

Country Status (1)

Country Link
JP (1) JP3471860B2 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111048168A (zh) * 2019-11-25 2020-04-21 安徽名流健康管理有限公司 多普勒超声诊断质量追溯系统

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111048168A (zh) * 2019-11-25 2020-04-21 安徽名流健康管理有限公司 多普勒超声诊断质量追溯系统

Also Published As

Publication number Publication date
JP3471860B2 (ja) 2003-12-02

Similar Documents

Publication Publication Date Title
EP1609423B1 (en) Ultrasonic diagnostic apparatus
JP5462966B2 (ja) 超音波画像処理方法および超音波画像処理装置
US6447455B2 (en) Ultrasound diagnostic apparatus and method for measuring blood flow velocity using doppler effect
US6530890B2 (en) Ultrasound diagnostic apparatus and method for measuring blood flow velocity using doppler effect
US5706818A (en) Ultrasonic diagnosing apparatus
JP5256210B2 (ja) 超音波画像処理方法および超音波画像処理装置
US6773403B2 (en) Ultrasonic apparatus and method for measuring the velocities of human tissues using the doppler effects
JP3093823B2 (ja) 超音波ドプラ診断装置
KR100483783B1 (ko) 도플러 신호 처리 장치 및 초음파 진단 장치
US20030073905A1 (en) Apparatus and method for indicating mechanical stiffness properties of body tissue
US8758249B2 (en) Method and apparatus for tissue doppler imaging
JP3281435B2 (ja) 超音波ドプラ診断装置
JP3471860B2 (ja) 超音波診断装置
JP2001249057A (ja) 大気温度測定方法および装置
JP4086684B2 (ja) パルス圧縮法及びパルス圧縮装置
JP3314991B2 (ja) 超音波ドプラ装置
JP4698073B2 (ja) 超音波診断装置
JP2002017725A (ja) 超音波ドップラ診断装置
JPH05200024A (ja) 超音波ドプラ診断装置
JPH0947456A (ja) 超音波ドプラ診断装置
JPH0518943A (ja) 内部欠陥探知方法およびその装置
JPH04146738A (ja) 超音波診断装置
JPH02264644A (ja) 超音波診断装置
Tanaka et al. Estimation of Doppler shift frequency using selected phase information for high frame rate color flow mapping
JPS6284742A (ja) 非線形パラメ−タ測定装置

Legal Events

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

Free format text: PAYMENT UNTIL: 20070912

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20080912

Year of fee payment: 5

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

Year of fee payment: 5

Free format text: PAYMENT UNTIL: 20080912

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

Free format text: PAYMENT UNTIL: 20090912

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20090912

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20100912

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20100912

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20110912

Year of fee payment: 8

LAPS Cancellation because of no payment of annual fees