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

超音波診断装置

Info

Publication number
JPH0866397A
JPH0866397A JP7021795A JP2179595A JPH0866397A JP H0866397 A JPH0866397 A JP H0866397A JP 7021795 A JP7021795 A JP 7021795A JP 2179595 A JP2179595 A JP 2179595A JP H0866397 A JPH0866397 A JP H0866397A
Authority
JP
Japan
Prior art keywords
speed
data
scale
velocity
diagnostic apparatus
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
JP7021795A
Other languages
English (en)
Other versions
JP3833282B2 (ja
Inventor
Nobuo Yamazaki
延夫 山崎
Fumiyasu Sakaguchi
文康 坂口
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 JP02179595A priority Critical patent/JP3833282B2/ja
Priority to US08/415,479 priority patent/US5513640A/en
Publication of JPH0866397A publication Critical patent/JPH0866397A/ja
Priority to US09/073,876 priority patent/USRE38209E1/en
Application granted granted Critical
Publication of JP3833282B2 publication Critical patent/JP3833282B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details 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/52053—Display arrangements
    • G01S7/52057—Cathode ray tube displays
    • G01S7/52071—Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
    • G01N29/0609—Display arrangements, e.g. colour displays
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8979—Combined Doppler and pulse-echo imaging systems
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/899—Combination of imaging systems with ancillary equipment
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details 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/52053—Display arrangements
    • G01S7/52057—Cathode ray tube displays
    • G01S7/5206—Two-dimensional coordinated display of distance and direction; B-scan display

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Chemical & Material Sciences (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

(57)【要約】 【目的】ＴＤＩで得られる低速領域の運動に対する計測
機能を活かし、低速度領域の表示能を向上させ、関心部
位の正常／異常の見分けを容易にする。 【構成】エコー信号から抽出したドプラ信号に基づいて
組織の運動の速度データを求める手段と、エコー信号か
ら断層像のデータを形成する手段と、超音波パルス信号
のパルス繰返し周波数で制限される測定可能なドプラ周
波数に対する速度測定範囲の速度データと当該速度デー
タに割り当てるカラー表示の階調データとをその速度測
定範囲内の低速度域をその残りの速度域よりも強調した
状態で対応付けるスケールを設定する手段と、スケール
に基づいて速度データを階調データに変換する手段と、
階調データ及び速度データの何れか一方が指定しきい値
を超えるサンプル点をブランク処理に付す手段と、ブラ
ンク処理を経た速度データを断層像データに重畳すると
共に速度データをカラー化して表示する手段とを備え
た。

Description

【発明の詳細な説明】
【０００１】
【産業上の利用分野】本発明は超音波診断装置に係り、
特に超音波パルスドプラ法を用いた組織ドプライメージ
ング（Tissue Doppler Imaging: ＴＤＩ）に適用した超
音波診断装置に関する。
【０００２】
【従来の技術】従来、組織ドプライメージングを行う超
音波診断装置としては例えば本出願人が既に提案してい
る特開平６−１１４０５９号（発明の名称：超音波カラ
ードプラ断層装置）がある。この公報記載の超音波診断
装置は、心筋や血管壁の組織の運動速度をパルスドプラ
法とローパスフィルタを用いて検出し、その運動速度か
ら運動に関する種々の物理量を演算して、演算結果を適
宜な態様でカラー表示する機構を備えている。この組織
の運動速度を検出にあたっては、組織の運動速度が血流
に比べて著しく低いことから、送信する超音波パルス
（レートパルス）のパルス繰返し周波数（Pulse Repeti
tion Frequency：ＰＲＦ）を下げて組織の超低速の運動
速度を計測可能にしている。
【０００３】この演算結果のカラー表示には種々の態様
があり、上記公開公報に係る発明の中でも二次元のカラ
ー表示を提案している。このカラー表示の階調に関して
は、組織ドプライメージングと共通の基盤に立つ、カラ
ードプラ装置による血流イメージングにおける手法を適
用することができる。
【０００４】上記血流イメージングでは、ドプラ偏移周
波数ｆｄが「−ｆｒ／２〜ｆｒ／２」（ｆｒ：超音波パ
ルス信号のパルス繰返し周波数：ＰＲＦ）の範囲を色の
輝度または色相を変えて、例えば図２７に示す如く３２
階調（１階調当りｆｒ／３２）で表示している。すなわ
ち、「−ｆｒ／２〜ｆｒ／２」の速度範囲全体に、速度
（ドプラ偏移周波数）に対する階調変化が一様なスケー
ルを割り充てて、赤（黄）−青（水色）に係るカラー表
示の階調を設定している。
【０００５】
【発明が解決しようとする課題】上述の如く、組織ドプ
ライメージングは高いパルス繰返し周波数（＝高フレー
ムレート）に設定して組織の超低速の運動速度を計測可
能にしたことを特徴の１つとしており、その「超低速」
という特性から、組織ドプライメージングで表示に特に
要求されるドプラ偏移周波数の範囲は血流イメージング
の場合よりも狭く、例えば「−ｆｒ／８〜ｆｒ／８」程
度である。
【０００６】これにも拘らず、現状では血流イメージン
グにおけるカラー階調の割当てをそのまま組織ドプライ
メージングに適用せざるを得ない。この結果、低速度領
域に割り当てられる階調数も極く僅かとなり、観察した
い心筋などの組織部位が例えば殆ど同じ色合いの赤また
は殆ど同じ輝度の赤で表示され、低速度領域における速
度差を視覚的に正確に評価することが極めて困難である
と共に、ｆｒ／３２以下の速度差を検出できても表示さ
れる色合いまたは輝度は変わらないことから、検出側の
高精度化が表示能の低さで相殺されてしまい、高検出精
度の機能を充分に発揮させることができないという未解
決の問題があった。
【０００７】さらに組織ドプライメージングを利用する
診断にあっては心筋などの組織の関心部位が正常か異常
かを早く見分けたい。これに対し従来の表示方法によれ
ば、カラー階調の割当てが低速度域から高速度域まで一
様であるため、前述の如く低速度域の表示能が低く、診
断部位の正常／異常の見分けが難しくなり、診断に長時
間を必要とし、また検査者に高度な熟練度を要求すると
いう未解決の問題があった。
【０００８】本発明は上述した未解決の問題に鑑みてな
されたもので、組織ドプライメージングで得られる低速
領域の運動に対する計測機能を活かし、その低速度領域
の表示能を向上させることを、その第１の目的とする。
【０００９】また、本発明は、上記第１の目的を達成し
つつ、関心部位の正常／異常の見分けが容易なＴＤＩ像
を提供することを、第２の目的とする。
【００１０】
【課題を解決するための手段】上記目的を達成するた
め、請求項１〜１３記載の発明は、被検体の断層面に含
まれる心筋，血管壁などの運動に関するカラー画像を表
示する超音波診断装置において、超音波パルス信号を前
記断層面に沿って走査して当該断層面からの超音波反射
信号に応じた電気量のエコー信号を得る走査手段と、前
記エコー信号から前記組織の運動に従うドプラ信号を抽
出する抽出手段と、前記ドプラ信号に基づいて前記組織
の運動の速度データを前記断層面のサンプル点毎に求め
る速度算出手段と、前記超音波パルス信号のパルス繰返
し周波数で制限される測定可能なドプラ周波数範囲の前
記速度データと当該速度データに割り当てるカラー表示
の階調データとをそのドプラ周波数範囲内の低速度域を
その残りの速度域よりも強調した状態で対応付けるスケ
ールを設定するスケール設定手段と、前記スケールに基
づいて前記速度データを前記階調データに変換する速度
変換手段と、この速度変換手段により変換された階調デ
ータをカラー化して表示する表示手段とを備えたことを
特徴とする。前記抽出手段は、前記組織によりドプラ偏
移を受けた前記ドプラ信号を選択的に抽出するローパス
フィルタを備える。前記スケールは、前記速度データの
変化に対する前記階調データの変化の割合が非線形に設
定されている。例えば、前記スケールは、前記低速度域
の前記割合を前記残りの速度域の前記割合よりも高く設
定されている。前記低速度領域は、例えば、前記超音波
パルス信号のパルス繰返し周波数をｆｒ，ドプラ偏移周
波数をｆｄとしたとき、−ｆｒ／８≦ｆｄ≦ｆｒ／８，
−ｆｒ／１２≦ｆｄ≦ｆｒ／１２，および−ｆｒ／１６
≦ｆｄ≦ｆｒ／１６のうちの１つである。また、前記階
調データは前記速度データの変化に応じて変わる連続的
なつながりを有する階調度のカラーデータを含むととも
に、当該階調データの少なくとも最高値はその低階調度
側とは階調度的に不連続なカラーデータである。
【００１１】また請求項１４記載の発明は、被検体の断
層面の断層像に当該断層面に含まれる組織の運動に関す
るカラー画像を重畳して表示する超音波診断装置におい
て、超音波パルス信号を前記断層面に沿って走査して当
該断層面からの超音波反射信号に応じた電気量のエコー
信号を得る走査手段と、前記エコー信号から前記組織の
運動に従うドプラ信号を抽出する抽出手段と、前記ドプ
ラ信号に基づいて前記組織の運動の速度データを前記断
層面のサンプル点毎に求める速度算出手段と、前記エコ
ー信号から前記断層像のデータを形成する断層像形成手
段と、前記速度データが速度しきい値を超えた場合、そ
の速度データをブランク処理するブランク処理手段と、
前記ブランク処理された前記速度データを前記断層像の
データに重畳するとともに当該速度データをカラー化し
て表示する表示手段とを備えたことを特徴とする。
【００１２】さらに請求項１５〜２７記載の発明は、被
検体の断層面の断層像に当該断層面に含まれる組織の運
動に関するカラー画像を重畳して表示する超音波診断装
置において、超音波パルス信号を前記断層面に沿って走
査して当該断層面からの超音波反射信号に応じた電気量
のエコー信号を得る走査手段と、前記エコー信号から前
記組織の運動に従うドプラ信号を抽出する抽出手段と、
前記ドプラ信号に基づいて前記組織の運動の速度データ
を前記断層面のサンプル点毎に求める速度算出手段と、
前記エコー信号から前記断層像のデータを形成する断層
像形成手段と、前記超音波パルス信号のパルス繰返し周
波数で制限される測定可能なドプラ周波数に対する速度
測定範囲の前記速度データと当該速度データに割り当て
るカラー表示の階調データとをその速度測定範囲内の低
速度域をその残りの速度域よりも強調した状態で対応付
けるスケールを設定するスケール設定手段と、前記スケ
ールに基づいて前記速度データを前記階調データに変換
する速度変換手段と、この速度変換手段により変換され
た階調データ及び前記速度算出手段により算出された前
記速度データの何れか一方の量が指定されたしきい値を
超える前記サンプル点をブランク処理に付すブランク処
理手段と、前記ブランク処理手段によるブランク処理を
経た前記速度データを前記断層像のデータに重畳すると
ともに当該速度データをカラー化して表示する表示手段
とを備えたことを特徴とする。前記スケール設定手段
は、前記ドプラ周波数の変化に対する前記階調データの
変化の割合を前記被検体の運動流体解析時の当該割合よ
りも高く、且つ前記階調データの最高値に対応する前記
速度データよりも大きい速度データには当該最高値を割
り当てる前記スケールを設定する手段である。さらに、
前記スケールの設定とは独立して前記しきい値を指定す
るしきい値指定手段を備えることができる。さらに、前
記スケール設定手段は、例えば、前記スケールの設定に
連動して前記しきい値を自動的に決める手段である。
【００１３】
【作用】本発明の一態様に係る超音波診断装置は、組織
ドプライメージング（ＴＤＩ）に適用される。このイメ
ージングの場合、心筋などの部位が超音波パルスドプラ
法で走査されてエコー信号が得られ、走査断層面のサン
プル点毎に運動速度が演算される。この運動速度は二次
元のカラー表示に付される。この表示のとき、超音波パ
ルスドプラ法で測定可能な速度範囲「−ｆｒ／２≦ｆｄ
≦ｆｒ／２」（ｆｒ：超音波パルスのパルス繰返し周波
数、ｆｄ：ドプラ偏移周波数）のうち、例えば低速域
「−ｆｒ／８≦ｆｄ≦ｆｒ／８」における「速度デー
タ」対「カラー表示の階調データ」のスケールの勾配が
その他の速度域よりも上げられる。したがって、かかる
低速域の表示分解能が上がるから、心筋の低速の運動の
僅かな変化が色の輝度または色合いの多階調の変化によ
って感度良く表示される。このため、高いパルス繰返し
周波数にして低速運動の計測機能を上げても、その計測
機能が損われることもない。また、例えば、最高階調度
の速度を有する画素は、それ以下の速度の画素とは不連
続な色相で表示されるので、その見分けが容易で、視認
性が高くなる。
【００１４】本発明のほかの態様に係る超音波診断装置
では、エコー信号から断層像のデータが形成される一方
で、超音波パルス信号のパルス繰返し周波数で制限され
る測定可能なドプラ周波数に対する速度測定範囲の速度
データと当該速度データに割り当てるカラー表示の階調
データとを、その速度測定範囲内の低速度域をその残り
の速度域よりも強調した状態で対応付けるスケールが設
定される。このスケールに基づいて速度データが階調デ
ータに変換され、この階調データ及び速度データの何れ
か一方の量が指定しきい値を超えるサンプル点がブラン
ク処理に付される。このブランク処理を経た速度データ
は断層像のデータに重畳されるとともに、当該速度デー
タがカラー化されて表示される。つまり、組織運動の速
度がしきい値より高いサンプル点のＴＤＩカラー像は表
示されず、背景のＢモード断層像のみが表示される。こ
のため、しきい値を適宜に指定することで、階調度の無
い速度カラー像を排除又は最小限に抑えることができ
る。
【００１５】
【実施例】
（第１実施例）本発明の第１実施例を図１〜図６に基づ
き説明する。この第１実施例に係る超音波診断装置は、
組織としての心筋（心臓壁）のＴＤＩ（組織ドプライメ
ージング）画像を得る診断装置である。
【００１６】図１には、超音波診断装置のブロック構成
を示す。図に示すように、この超音波診断装置１０は、
被検者との間で超音波信号の送受信を担う超音波プロー
ブ１１と、この超音波プローブ１１を駆動し且つ超音波
プローブ１１の受信信号を処理する装置本体１２と、こ
の装置本体１２に接続され且つ心電情報を検出するＥＣ
Ｇ（心電計）１３と、装置本体１２に接続され且つオペ
レータからの指示情報を装置本体に出力可能な操作パネ
ル１４とを備える。
【００１７】装置本体１２は、その扱う信号経路の種別
に拠り超音波プローブ系統、ＥＣＧ系統及び操作パネル
系統に大別することができる。超音波プローブ系統とし
ては、超音波プローブ１１に接続された超音波送受信部
１５を備え、この超音波送受信部１５の出力側に配置さ
れたＢモード用ＤＳＣ（デジタルスキャンコンバータ）
部１６、Ｂモード用フレームメモリ（ＦＭ）１７、メモ
リ合成部１８及び表示器１９を備える一方、同じく超音
波プローブ１１に接続された、カラーマッピングのため
の位相検波部２０、フィルタ部２１、周波数解析部２
２、符号化演算部２３、ＴＤＩ用ＤＳＣ部２４、及びＴ
ＤＩ用フレームメモリ２５を備えている。また、ＥＣＧ
系統としては、ＥＣＧ１３に接続されたＥＣＧ用アンプ
４０を備え、このアンプ４０の出力側に接続されたトリ
ガ信号発生器４１及び参照データメモリ４２を備える。
さらに、操作パネル系統としては、操作パネル１４から
の操作情報を入力するＣＰＵ（中央処理装置）４３と、
このＣＰＵ４３の管理下に置かれるタイミング信号発生
器４４とを備える。なお、ＣＰＵ４３は、オペレータが
操作パネル１４を介して指令したＲＯＩ（関心領域）の
設定信号を、ＲＯＩ設定に必要な各構成に供給できるよ
うになっている。
【００１８】この実施例にあっては、超音波プローブ１
１及び超音波送受信部１５が本発明の走査手段を形成
し、位相検波部２０が抽出手段を成し、フィルタ部２１
及び周波数解析部２２が本発明の速度演算手段を形成し
ている。また、ＴＤＩ用ＤＳＣ部２４、ＴＤＩ用フレー
ムメモリ２５、メモリ合成部１８及び表示器１９が本発
明の表示手段を形成している。符号化演算部２３がスケ
ール設定手段及び速度変換手段に対応している。
【００１９】超音波プローブ１１は、短冊状の複数の圧
電振動子を配列させたフェーズドアレイ形のトランスデ
ューサを内蔵している。各圧電振動子は、超音波送受信
部１５からの駆動信号によって励振される。各駆動信号
の遅延時間を制御することにより、スキャン方向を変更
してセクタ電子走査可能になっている。超音波送受信部
１５の遅延時間パターンは、後述するタイミング信号発
生器４４から送られてくる基準信号を基準時として、Ｃ
ＰＵ４３により制御される。超音波送受信部１５は、ス
キャン方向に対応して遅延時間パターンが制御された駆
動電圧信号を超音波プローブ１１に出力する。この駆動
電圧信号を受けた超音波プローブ１１は、そのトランス
デューサにおいて電圧信号を超音波信号に変換する。こ
の変換された超音波信号は、被検者の心臓に向けて送波
される。この送波された超音波信号は、心臓を含む各組
織で反射され、再び超音波プローブ１１に戻ってくる。
そこで、プローブ１１内のトランスデューサでは反射超
音波信号が再び電圧信号（エコー信号）に変換され、そ
のエコー信号は超音波送受信部１５に出力される。
【００２０】上記超音波送受信部１５の信号処理回路
は、送信時と同様に、入力したエコー信号に遅延をかけ
て整相加算し、スキャン方向に超音波ビームを絞ったと
等価なエコービーム信号を生成する。この整相加算され
たエコービーム信号は、検波された後、Ｂモード用ＤＳ
Ｃ部１６に出力される。このＤＳＣ部１６は超音波走査
のエコーデータを標準テレビ走査のデータに変換し、メ
モリ合成部１８に出力する。また、これと並行して、Ｂ
モード用ＤＳＣ部１６は、任意の心位相における複数枚
の画像データをＢモード用フレームメモリ１７に記憶さ
せる。
【００２１】一方、超音波送受信部１５で処理されたエ
コー信号は、位相検波部２０にも出力される。位相検波
部２０はミキサとローパスフィルタを備える。心筋のよ
うな運動をしている部位で反射したエコー信号は、ドプ
ラ効果によって、その周波数にドプラ偏移を受けてい
る。位相検波部２０はそのドプラ周波数について位相検
波を行い、ドプラ信号のみをフィルタ部２１に出力す
る。
【００２２】フィルタ部２１は、運動速度の大きさが
「心筋＜弁＜血流」の関係にあることを利用して（図２
参照）、位相検波されたドプラ信号から、心臓壁以外の
弁運動、血流などの不要なドプラ成分を除去し、超音波
ビーム方向の心筋のドプラ信号を効率良く検出する。こ
の場合、フィルタ部２１はローパスフィルタとして機能
させる。
【００２３】上記フィルタ部は既に実用化されている、
血流情報を得るためのカラードプラ断層装置にも搭載さ
れているものである。この血流情報を得るカラードプラ
断層装置の場合には、血流と心臓壁、弁運動とのドプラ
信号が混在したエコー信号に対してハイパスフィルタと
して機能させ、血流以外のドプラ信号を除去している。
このため、フィルタ部は装置の目的に応じてローパスフ
ィルタとハイパスフィルタとを切換可能にすることで汎
用性を高めることができる。
【００２４】フィルタ部２１でフィルタリングされたド
プラ信号は、次段の周波数解析部２２に出力される。周
波数解析部２２は、超音波ドプラ血流計測で用いられて
いる周波数分析法である、ＦＦＴ法及び自己相関法など
を応用するものであり、スキャンされる断層面内の個々
のサンプル点における観測時間（時間窓）内での平均速
度や最大速度を速度データとして演算する。具体的に
は、例えば、ＦＦＴ法又は自己相関法を用いてサンプル
各点の平均ドプラ周波数（即ち、その点での観測対象の
運動の平均速度）や分散値（ドプラスペクトラムの乱れ
度）を、さらにはＦＦＴ法を用いてドプラ周波数の最大
値（即ち、その点での観測対象の運動の最大速度）など
をほぼリアルタイムで演算する。このドプラ周波数の解
析結果は運動速度のカラードプラ情報として次段の符号
化演算部２３に出力される。
【００２５】符号化演算部２３は、ＣＰＵ機能を備え、
周波数解析部２２から送られてくる、断層面の各サンプ
ル点毎のドプラ偏移周波数ｆｄを、指定された速度変換
スケールを使って所定ビット数の速度表示データに符号
化するようになっている。この超音波パルスドプラ法で
は、超音波パルス信号のパルス繰返し周波数ｆｒはサン
プリング周波数に相当する。よってサンプリング定理か
ら、測定可能なドプラ偏移周波数の最大値ｆｄmax は、
【数１】ｆｄmax ＝ｆｒ／２ である。
【００２６】周波数解析部２２において折返り現象が発
生しない演算可能なドプラ偏移周波数ｆｄは、
【数２】−ｆｒ／２≦ｆｄ≦ｆｒ／２ であるが、この周波数（すなわち速度）範囲のうち、
【数３】−ｆｒ／８≦ｆｄ≦ｆｒ／８ の周波数範囲を「ｆｒ／１２８」の量子化率で量子化し
て、例えば５ビットのデータ長の速度表示コードに符号
化する。この場合、ドプラ偏移周波数ｆｄが、
【数４】−ｆｒ／２≦ｆｄ＜−ｆｒ／８ のときは、一方の運動方向（例えば超音波ビームから遠
ざかる方向）の最高階調度に対応する、「−ｆｒ／８」
のときと同じ符号化データに設定し、
【数５】ｆｒ／８＜ｆｄ≦ｆｒ／２ のときは、もう一方の運動方向（例えば超音波ビームに
近づく方向）の最高階調度に対応する、「ｆｒ／８」の
ときと同じ符号化データに設定する。
【００２７】この結果、ドプラ偏移周波数の範囲「−ｆ
ｒ／２≦ｆｄ≦ｆｒ／２」を横軸にとり、カラー表示色
の赤（黄）−青（水色）の色相変化を表す速度表示コー
ドを縦軸にとったときの速度変換スケールは図３に示す
ように絶対値が「ｆｒ／８」より大きい速度では飽和す
る形で表わされる。このように符号化された、データ長
が５ビットの各サンプル点毎の速度表示データは次段の
ＴＤＩ用ＤＳＣ部２４に出力される。以上の手順の概略
を、図４のステップＳ１〜Ｓ５に示す。
【００２８】ＴＤＩ用ＤＳＣ部２４は、走査方式変換用
のＤＳＣ２４ａとコードに変換された速度表示データを
カラー化するためにルックアップ用テーブルを備えたカ
ラー回路２４ｂとを備えている。このため、符号化演算
部２３から送られてきた速度データは、ＤＳＣ２４ａで
超音波走査信号が標準テレビ走査信号に変換されると共
に、カラー回路２４ｂでカラー表示用データに変換さ
れ、その変換信号が前記メモリ合成部１８に出力され
る。
【００２９】ここで、上記カラー回路２４ｂで処理され
る心筋速度のカラー表示方式について触れる。このカラ
ー表示を大別すると、（ｉ）速度の大きさ（絶対値）の
表示、（ｉｉ）運動の方向と速度の大きさの表示、（ｉ
ｉｉ）運動の方向の表示、に分けられる。（ｉ）の表示
法としては、ａ：単色で大きさに応じて輝度を変える、
ｂ：大きさに応じて色を変える、がある。（ｉｉ）の表
示法については、方向を色で示し、大きさを輝度で示す
やり方のほか、方向を色合いで示し、大きさをその色合
いの変化で示すやり法などがあり、この内、方向につい
ては、得られる速度情報の態様に応じて、適用可能な表
現法が制限される。ここでは、ＴＤＩ用ＤＳＣ部２４の
カラー回路２４ｂにおいて、図５に示したように、カラ
ーが決められる。即ち、従来知られている超音波ビーム
に近づく運動を赤、超音波ビームから遠ざかる運動を青
で示す方法に対応させて、心筋の収縮運動を赤（黄）、
心筋の拡張運動を青（水色）で示し、且つ、その絶対値
が大きくなるにしたがって黄色又は水色にグラデーショ
ンを変化させるようにしたものである。この結果、所望
の低速度領域「−ｆｒ／８≦ｆｄ≦ｆｒ／８」の速度表
示データは「水色→青→赤→黄」の一連のカラー表示色
に対して運動方向毎に３２階調の色情報に変換される。
【００３０】また、ＴＤＩ用ＤＳＣ部２４のＤＳＣ２４
ａはさらに、任意の心時相における複数枚のカラードプ
ラ画像をＴＤＩ用フレームメモリ２５に記憶させる。
【００３１】一方、前述したＥＣＧ１３は被検者の各心
時相の心電図情報を検出するようになっている。この検
出信号は、ＥＣＧ用アンプ４０を経てトリガ信号発生器
４１及び参照データメモリ４２に各々出力される。この
内、参照データメモリ４２は各心時相における心電図情
報を記憶しておき、必要に応じて必要な情報をメモリ合
成部１８に供給する。トリガ信号発生器４１は、各心時
相のタイミング情報を前記タイミング信号発生器４４に
知らせるようになっている。タイミング信号発生器４４
は、通常、操作パネル１４からの指示に応じて超音波送
受信部１５における遅延時間パターンを制御するＣＰＵ
４３のコントロール下にあるが、トリガ信号発生器４１
から各心時相のタイミングが告知されると、超音波送受
信部１５に対して超音波送受のための基準信号を発振す
る。
【００３２】上述したようにメモリ合成部１８には、Ｂ
モード用ＤＳＣ部１６から出力されたＢモード画像信
号、ＴＤＩ用ＤＳＣ部２４から出力されたＴＤＩモード
の画像信号、さらには必要に応じて前記参照データメモ
リ４２からの心電図情報が入力するようになっている。
メモリ合成部１８では、それらの入力信号データが重畳
され、その重畳データが表示器１９に出力される。表示
器１９はここではＣＲＴで成る。
【００３３】この結果、血流や弁のドプラ信号は既にフ
ィルタ部２１でカットされているから、表示器１９には
心臓のＢモード断層像（白黒階調）と、心筋の動きを図
３に示す速度変換スケールに従って、図５に示すカラー
バーで色分けしたカラー画像とを重畳させた断層像が、
例えば図６に示すように殆どリアルタイムに表示される
（同図においてハッチング部分が心筋ＨＭを示す）。つ
まり、図６に示す心筋ＨＭのカラーは収縮運動時には赤
（黄）、拡張運動時には青（水色）となり、その赤、青
が周期的に且つリアルタイムに繰り返される。しかも収
縮、拡張運動の最中における運動速度の変化は、赤もし
くは黄または青もしくは水色の色合い変化によって殆ど
リアルタイムに表現される。よって、心筋ＨＭの運動速
度をカラーでほぼリアルタイム且つ精度良く表示させる
ことができ、心臓の機能低下を定量的且つ高精度に評価
するための基礎画像を取得できる。
【００３４】特に、組織ドプライメージングにおいて関
心のある超低速度の範囲「−ｆｒ／８≦ｆｄ≦ｆｒ／
８」がカラーバーの全階調度範囲に相当する運動方向毎
の３２階調（量子化率５ビット）の色合いで表示される
ため、ここでは１階調当り「ｆｒ／１２８」の、刻みが
細かいドプラ偏移周波数（すなわち運動速度）が実質的
に割り当られることになる。これは、前述したフルスケ
ール「−ｆｒ／２≦ｆｄ≦ｆｒ／２」に３２階調を割り
当てる従来法の場合に比べて、超低速度範囲「−ｆｒ／
８≦ｆｄ≦ｆｒ／８」の階調表示能が４倍に向上したこ
とになる。これにより、高いパルス繰返し周波数にして
検出した心筋の超低速の運動速度が、従来よりも各段に
多い複数段の階調で表示されるから、超低速域内の僅か
な速度差が異なる色合いの表示色となり、速度差を視覚
的に容易に差別化して評価することができる。
【００３５】本出願人の臨床評価では、上記低速度域の
最大速度範囲は４cm／s 及び１０cm／s の間で選択し、
この最大速度範囲よりも大きい速度は飽和して最高輝度
（色相変化に代えて指定カラーの輝度変化をカラー表示
用階調データとして採用した）の赤または青の速度表示
データに符号化されるようにした。しかしながら、サン
プリング周波数から得られるエイリアシング速度はその
飽和した最大速度よりも４倍または８倍高く、３０cm／
s 及び４０cm／s の間であった。この条件のもとで心室
壁の運動速度を計測したが、エイリアシングは発生しな
いことが確認できた。
【００３６】よって、従来と同じパルス繰返し周波数の
超音波パルス信号によるスキャンであっても、組織の超
低速の運動速度に対する計測機能を損わないで済み、装
置の高機能化を図ることができる。
【００３７】なお、上記実施例における診断装置はＢモ
ード用とＴＤＩ用の２種類のフレームメモリ１７、２５
を備えているため、必要に応じて、スローモーション再
生、コマ送り再生などのシネループ再生や動画再生を行
ったり、心時相が異なる画像をＢモード用とＣＦＭ用と
で個別に或いは並列に表示させることができる。
【００３８】また、上記断層装置には、心筋の動きをド
プラ表示させるためのドプラフィルタやＦＦＴ（高速フ
ーリエ変換）周波数分析器を付加することもできる。
【００３９】さらに、上記実施例では心筋ＴＤＩ画像を
重畳させる画像がＢモード断層像であり、また診断対象
が心臓である構成について説明してきたが、この発明は
必ずしもそのような構成に限定されるものではない。例
えば、Ｂモード像の代わりに、Ｍモード像であってもよ
いし（この場合には、Ｂモード像取得のための各構成要
素をＭモード像のそれに置換すればよい）、心筋の代わ
りに血管壁を診断してもよい（この場合には、フィルタ
部２１のカットオフ周波数を血管壁用に合わせる）。ま
た、それらＢモード像やＭモード像を重畳しないで、Ｔ
ＤＩ（カラードプラ）像のみを単独で表示させてもよ
い。
【００４０】さらに、通常のＢモード断層装置及び血流
カラーフローマッピングで見られるように、心電図など
の生体信号との対応を明確にするため、生体信号波形の
同時表示や、心電図Ｒ波などからの時間差表示を行って
もよい。
【００４１】さらにまた、上記実施例における周波数解
析部２２と符号化演算部２３との間に、絶対速度演算部
を挿入して、心筋などの組織の運動の絶対速度（すなわ
ち、各サンプル点における組織の運動方向の速度それ自
体）を例えば特開平６−１１４０５９号に示す如く推定
演算させ、この絶対速度を二次元カラー表示することも
できる。
【００４２】一方、本発明に係る階調表示能を階調する
低速度領域については、前述した実施例の「−ｆｒ／８
≦ｆｄ≦ｆｒ／８」に限定されることなく、符号化演算
部のプログラムなどを組み換えることにより、例えば図
７中の一点鎖線の速度変換スケールで示す如く「−ｆｒ
／１２≦ｆｄ≦ｆｒ／１２」の範囲や、同図中二点鎖線
の速度変換スケールで示す如く、「−ｆｒ／１６≦ｆｄ
≦ｆｒ／１６」の範囲に設定するとしてもよい。また、
それらの周波数範囲「−ｆｒ／８≦ｆｄ≦ｆｒ／８」、
「−ｆｒ／１２≦ｆｄ≦ｆｒ／１２」、「−ｆｒ／１６
≦ｆｄ≦ｆｒ／１６」をオペレータからのマニュアル操
作信号に応じて切り換え、表示器の画面を見ながら適宜
なレンジを選択するようにしてもよい。これには、図１
において、操作パネル１４からのマニュアル操作信号を
受けたＣＰＵ４３が、その操作に係る切換信号を符号化
演算部２３に送るようにすればよい。
【００４３】さらに、本発明に係る低速度域強調表示に
あっては、例えば図８に示す如く、所望の低速度領域、
例えば「−ｆｒ／１２≦ｆｄ≦ｆｒ／１２」については
勾配のより大きい速度変換スケールを割り当て、その外
側の中速度領域についてはより緩やかな勾配の速度変換
スケールを割り当てるようにしてもよく、同図に示す複
数段の速度変換スケールの勾配の切替制御を符号化演算
部で行なうようにすればよい。これにより低速度領域と
その周辺部との速度関係にも若干、配慮した二次元カラ
ー像が得られる。
【００４４】さらに本発明に係る低速度域の強調表示法
にあっては図９の実線に示す速度変換スケールを使うこ
ともできる（同図の１点鎖線は従来の血流解析時のスケ
ールである）。この速度変換スケールは予め符号化演算
部２３に例えば記憶テーブルの形で記憶させているもの
である。低速度域として定めた、例えば−ｆｒ／８＜ｆ
d ＜＋ｆｒ／８の範囲では組織運動の方向毎に赤色から
黄色に及び青色から水色に徐々に（すなわち連続的に）
変化する色相の階調データに対応する速度表示コードが
割り当てられている。しかし、平均化されたドプラ偏移
周波数ｆｄ（すなわち組織の平均運動速度）がｆｄ≧±
ｆｒ／８に達すると、それまでの赤から黄または青から
水色への連続変化する色相に変え、色相が全く連続しな
い特殊な表示色ＣＬ１，ＣＬ２の速度表示コードを一律
に割り当てるものである。この特殊表示色ＣＬ１，ＣＬ
２は、例えば赤系，青系の色相に一定の色相を混ぜて生
成される色である。勿論、このときの速度しきい値は±
ｆｒ／８に限定されず、変更可能な値である。これによ
り、予め定めた低速度域を超える組織の運動速度は不連
続な色相に拠って一目で見分けが付くから、低速度域の
強調表示と合せて診断のための画像観察の容易化が一層
押し進められる。
【００４５】さらに、本発明のカラー表示に係る速度表
示コードの階調データは前述したように速度の大きさに
応じて色相を変えるもののほか、速度の大きさに応じて
赤（または青）の輝度を変化させるようにしてもよい。
【００４６】さらにまた、前記実施例及びその変形例は
必要に応じて従来の血流ドプライメージングを実施する
超音波診断装置に適用することもできる。
【００４７】（第２実施例）次に本発明の第２実施例を
図１０〜図１５に基づいて説明する。前述した第１実施
例では低速度領域の表示能を向上させることを目的とし
ていたが、この第２実施例の超音波診断装置は、かかる
表示能向上に加えて、関心領域の正常／異常を見分ける
診断を容易化させるものである。この目的を達成するた
め、この第２実施例の超音波診断装置は図１０に示すよ
うに構成されている。
【００４８】具体的には、この超音波診断装置は、被検
体との間で超音波信号の送受信を担う超音波プローブ１
００と、この超音波プローブ１００を駆動し且つ超音波
プローブ１００の受信信号を処理する装置本体１０１と
を備える。
【００４９】この内、超音波プローブ１００は第１実施
例と同様にフェーズド・アレイ形に形成されている。装
置本体１０１は、この超音波プローブ１００に接続され
た超音波送受信部１１０を備え、この出力側に位相検波
器１１１，Ａ／Ｄ変換器１１２，フィルタ１１３，運動
速度解析部１１４，ＤＳＣ１１５，カラー処理器１１
６，Ｄ／Ａ変換器１１７，及びカラーモニタ１１８が順
次接続されている。
【００５０】また、運動速度解析部１１４には、後述す
る組織運動のカラー画像の少なくとも一部をブランク処
理するためのブランク制御器１２１が接続されている、
このブランク制御器１２１には入力器１２２を介して検
査者から必要な情報が与えられる。
【００５１】上記超音波送受信部１１０は、与えられる
レートパルスの周期で超音波プローブ１００を駆動する
とともに、超音波プローブ１００からのエコー信号を整
相加算する送受信回路１１０ａ、この送受信回路１１０
ａにラスタアドレスなどの必要な情報を与えるＲＰＧ
（レートパルスジェネレータ）１１０ｂ、及びＢモード
用画像信号を生成する包絡線検波器１１０ｃとを備え、
第１実施例と同等に機能する。
【００５２】送受信回路１１０ａの出力側は位相検波器
１１１，Ａ／Ｄ変換器１１２，フィルタ１１３に順次接
続されており、これらのユニットは第１実施例における
位相検波部２０，フィルタ部２１（Ａ／Ｄ変換器を含
む）と同等に機能する。
【００５３】周波数解析部１１４は、自己相関法などに
拠り断層面の各サンプル点のドプラ周波数解析を行う周
波数解析器１１４ａを有するとともに、その解析結果に
基づいて、各サンプル点の平均ドプラ周波数（平均速
度）を演算する速度演算器１１４ｂ，分散値（スペクト
ラムの乱れ度）を演算する分散演算器１１４ｃ及び強さ
（パワー）を演算するパワー演算器１１４ｄを備える。
【００５４】速度演算器１１４ｂは本実施例では図１１
に示すように、ＣＰＵ１１４０及びメモリ１１４１など
を備えている。メモリ１１４１の所定記憶領域には図１
２に示す処理のプログラムが予め内蔵されており、速度
演算器１１４ｂの起動と共に、係る処理が自動的に実行
される。なお、この速度演算器１１４ｂは図１２と同等
に機能するようにアナログ，デジタルの電子回路などを
組み合せて構成してもよい。
【００５５】ブランク制御器１２１は入力器１２２から
供給される信号を解読するとともに、その解読値に対応
した、ブランク処理に必要な偏角値θのしきい値θthの
しきい値信号Ｓθth及びスケール変換係数Ｋの計数信号
ＳK を速度演算器１１ｂに供給する。
【００５６】ＤＳＣ１１５は、包絡線検出波器１１０ｃ
から出力された白黒のＢモード断層像の画像データと速
度演算器１１４ｂ，分散演算器１１４ｃ，及びパワー演
算器１１４ｄから出力されたＴＤＩ像の画像データとを
入力し、Ｂモード像にＴＤＩ像を重畳（合成）したフレ
ーム像データを形成する。このフレーム像データはカラ
ー処理器１１６に送られる。カラー処理器１１６はＴＤ
Ｉ像の画素には速度表示コードに対応した色付けを行
い、そのカラー化したフレーム像データをＤ／Ａ変換器
１１７を介してカラーモニタ１１８に送るようになって
いる。
【００５７】ここで図１２に基づいて速度演算器１１４
ｂの処理を説明する。
【００５８】速度演算器１１４ｂのＣＰＵ１１４０は、
まずステップ２００でブランク制御器１２１から供給さ
れるしきい値信号Ｓθthを読み込む。次いでステップ２
０１に移行してしきい値信号Ｓθthで指定されたしきい
値θthを記憶する。速度演算器１１４ｂは、後述するよ
うに、周波数解析器１１４ａの解析結果である平均ドプ
ラ周波数（平均速度）の例えば自己相関係数Ｒｅ，Ｉｍ
（複素数），Ｐｏ（パワー）の内のＲｅ，Ｉｍを入力し
て、組織の運動速度Ｖに対応する複素平面上の単位円上
の点を表わす偏角θを求めるので（図１３参照）、上記
しきい値θthは図１４に示す如く運動速度Ｖ（すなわち
ドプラ周波数）のしきい値Ｖthに相当する。
【００５９】ＣＰＵ１１４０は次いでステップ２０２に
移行し、ＴＤＩの表示能強調用のスケール変換係数Ｋ
（＞１）を指定する係数信号ＳK をブランク制御器１２
１から読み込む。そしてステップ２０３で係数信号ＳK
に対応するスケール変換係数Ｋを記憶する。この係数Ｋ
＝１のときは血流速度解析時である。
【００６０】次いで、ステップ２０４で、ＣＰＵ１１４
０は周波数解析器１１４ａで解析された例えば自己相関
係数Ｒｅ，Ｉｍを入力し、ステップ２０５で前述した組
織の運動速度Ｖに相当する偏角θを演算する。このとき
の偏角θは、スケール変換係数Ｋ＝１、すなわち血流速
度解析に相当する値であり、図１４の１点鎖線ｄで示す
直線（速度変換スケール）上を動く。同図は横軸に運動
速度Ｖ、縦軸に赤（プローブに近付く運動方向）及び青
（プローブから離れる運動方向）の階調データとしての
輝度階調を表わす速度表示コードＣＤV （例えば８ビッ
トの論理値データ）を各々とっている。同図上には、後
述するように本発明に係る種々の直線、すなわち速度変
換スケールを引けるが、前記直線ｄは血流解析時に相当
する速度変換スケールである。この速度変換スケールｄ
は周知の如く、折り返り速度±ｆｒ／２の範囲内で全て
の速度Ｖに連続的に変化する速度表示コードＣＤV が一
様に割り当てられる。
【００６１】さらにステップ２０６では、ＣＰＵ１１４
０により、ステップ２０５で演算した偏角θがステップ
２０１で設定したしきい値θthに対して、θ＞θth（す
なわち、Ｖ＞Ｖth）か否かが判断される。この判断でＹ
ＥＳ（θ＞θth）となる場合、ステップ２０７に移行し
てそのθ＞θth（Ｖ＞Ｖth）となった運動速度を有する
画素に対するブランク処理が指令される。このブランク
処理は、その画素の速度表示コードＣＤV を強制的にＶ
＝０のときの速度表示コードＣＤV ＝「０，０，…，
０」（ブランクコード）に設定する処理である。
【００６２】このブランク指令が終ったとき、又はステ
ップ２０６の判断でＮＯ（θ≦θth）となるときは、ス
テップ２０８に処理が進められる。このステップ２０８
ではステップ２０３で設定されたスケール変換係数Ｋ
（＞１）にステップ２０５で演算した偏角値θが掛けら
れる。
【００６３】この係数の乗算「Ｋ・θ」によって偏角値
θ、すなわち組織の運動速度ＶがＴＤＩ像の低速度域の
強調表示態様に変換されるようになる。例えばＫ＝Ｋ1
（＞１；例えば「４」）となる所定値にする速度変換ス
ケールは図１４上で血流解析時に相当する直線ｄからａ
に移る。既にブランキング処理されている偏角値θ（＝
０）に係数Ｋを乗じてもその値は零であるからＴＤＩ用
の一例を示す速度変換スケールａは図１４に示すよう
に、Ｖ＝±Ｖａ（例えば±ｆｒ／８）の間で表示色
「赤」及び「青」の設定輝度階調全部の速度表示コード
を使う直線で表わされ、Ｖａ≧｜Ｖ｜≧Ｖthの範囲では
最大輝度に対応する速度表示コード±CDV(MAX)で飽和し
た形となり、Ｖ＝±Ｖthで原点を通る横軸、すなわち黒
の表示色に対応まで立ち下がる。
【００６４】また別の例として係数Ｋ＝Ｋ2 （＞Ｋ1 ：
例えば「８」），Ｋ＝Ｋ3 （＞Ｋ2：例えば「１６」）
を偏角値θに乗じたときの速度変換スケールは図１４上
の特性線ｂ，ｃのようになる。速度変換スケールｂはＶ
＝±Ｖｂ（例えば±ｆｒ／１２）の間で直線的に増減
し、Ｖｂ≧｜Ｖ｜≧Ｖthの間で飽和して赤，青の最大輝
度の速度表示コード±ＣＤV (MAX) に設定される。もう
一方の速度変換スケールｃの場合、Ｖ＝±Ｖｃ（例えば
±ｆｒ／１６）の間で同様に直線的に増減し、Ｖｃ≧｜
Ｖ｜≧Ｖthの間で飽和して最大輝度の速度表示コード±
ＣＤV (MAX) に設定される。両スケールｂ，ｃ共に｜Ｖ
｜≧Ｖthの範囲ではブランク処理により立ち下がる。こ
のように、乗じる係数Ｋの値に比例して速度変換スケー
ルの直線部分が急峻になる。このためスケール変換係数
Ｋの値に比例してその低速度域の階調性のある部分の表
示能が強調されることになる。
【００６５】このようにＴＤＩ表示のためのスケール変
換係数Ｋの乗算が終ると、ＣＰＵ１１４０は図１２のス
テップ２０９に移って、Ｋ倍された偏角値θ（すなわち
速度Ｖ）の速度表示コードＣＤV への変換を行う。この
変換は予めメモリ１１４１に格納してある偏角値θを速
度表示コードＣＤV （例えば８ビット値）の対応を示す
記憶テーブルを参照することで行われる。次いでステッ
プ２１０では、変換した速度表示コードＣＤV をＤＳＣ
１１５に出力する。
【００６６】さらにステップ２１１で再びしきい値信号
Ｓθth及び係数信号ＳK のステップ２１２で偏角値θの
しきい値θth（Ｖth）又は／及び係数Ｋを変更するか否
かを判断し、変更する場合はステップ２０１又は２０３
の処理に戻る。しきい値θthや係数Ｋを変更しないでブ
ランク処理を続ける場合、ステップ２１３を経てステッ
プ２０４に戻る。ブランク制御は以上のようにしてなさ
れる。
【００６７】超音波プローブ１００より生体内に超音波
ビームを送信して得られたエコー信号は再びプローブ１
００により受信される。このエコー信号はプローブ１０
０で電気量の受信信号に変換され、送受信回路１１ａで
受信処理され、さらに位相検波器１１１で直交検波され
る。検波により検出されたドプラ信号はフィルタ１１３
により組織運動に相当する周波数成分を含む成分が抽出
され、クラッタ成分を含んだままの当該ドプラ信号が周
波数解析器１１４ａに送られる。この周波数解析結果
は、速度演算器１１４ｂ，分散演算器１１４ｃ，パワー
演算器１１４ｄに送られて、組織運動に関する目的とす
る量が演算される。パワー演算器１１４ｄでは「Ｋ・ l
og10Ｐｏ」の処理が行われる。
【００６８】今、図１４に示す如く、速度演算器１１４
ｂに、入力器１２２およびブランク制御器１２１を介し
て、偏角値θに対するしきい値θth（すなわち組織運動
Ｖに対する速度のしきい値Ｖth）が指令されており、Ｔ
ＤＩ表示に係るスケール変換係数Ｋ＝Ｋ2 （直線ｂ）が
指令されているとする。この場合、図１４上でのコード
変換特性全体は一点鎖線ｍｂで表わされる。
【００６９】このため、ＴＤＩの強調表示のために演算
された組織運動の平均速度（Ｋ・θの値）が｜Ｖ｜＜±
Ｖｂの低速度範囲に入るとき、許容される速度表示コー
ドＣＤV 全体が使われる。速度Ｖが±Ｖｂ≦｜Ｖ｜＜±
Ｖthの範囲に収まるときは±ＣＤV (MAX) の最高輝度の
コードで飽和し、±Ｖth≦｜Ｖ｜となる速度はブランク
コードが割り当てられる。
【００７０】このようにブランク処理された平均速度デ
ータを含む運動速度解析結果はＤＳＣ１１５に送られ
る。このＤＳＣ１１５には包絡線検波器１１０ｃからＢ
モード像データも送られており、このＢモード像に運動
速度情報が重畳される。この重畳されたフレーム像デー
タはカラー処理器１１６で色付けされ、カラーモニタ１
１８に表示される。この表示像は白黒のＢモード像を背
景として、これに組織運動のカラー像などの運動速度情
報が重畳したものであるが、指定速度Ｖth以上の速度値
を有する画素は図１２のステップ２０７によりブランク
処理されているので、その画素のカラー像は表示されな
い。すなわち、±Ｖth≦｜Ｖ｜となる高速度域の画素に
はカラー情報が上書きされず、背景像のＢモード像が検
査者に見えることになる。
【００７１】このため、高速度域をカットするために指
定する速度しきい値Ｖth（具体的には偏角θのしきい値
θth）を適宜な値に設定することで、心筋などの組織の
運動の階調性の無い最高輝度部分は必要最小限の面積
（画素）に抑えることができる。これにより、見易い画
面となり、組織運動の解析にとっては無用な最高輝度部
分が診断の邪魔になるという事態を殆ど回避でき、診断
の能率（診断時間，診断労力）向上に寄与できる。
【００７２】またそのような高輝度像に代えて背景像
（白黒Ｂモード像）が現われるので、かえって診断に有
効な情報が増えるという利点がある。
【００７３】当然に、低速度域を強調する速度変換スケ
ールａ，ｂ，ｃ（図１４参照）を用いているので、壊死
した組織の遅い運動速度を容易に見分けることができ、
関心部位の正常／異常の診断能が向上する。
【００７４】さらに、本実施例では運動速度Ｖのしきい
値Ｖthを速度変換スケールａ（…ｃ）とは独立して設定
でき、モニタ画面上に残す赤、，青の最高輝度±CD
V (MAX)の部分を、同一の速度変換スケールａ（…ｃ）
（すなわち、同一の低速度強調機能）に対して適宜に調
整できる。このため、速度しきい値Ｖth及びスケール変
換係数Ｋの選択具合によっては図１５の１点鎖線ｍ1 で
示すように、最高階調の輝度領域をモニタ画面上に残す
こともできるし、同図実線ｍ2 で示すように最高階調の
輝度まで達しない内にＴＤＩカラー像のみをブランクに
してしまうこともできる。
【００７５】本実施例ではさらに組織運動解析用に複数
の速度変換スケールａ…ｃを予め用意し、それらを適
宜、選択／切換できるようになっている（図１２中ステ
ップ２０２，２１１参照）。このため非常に汎用性が高
められる。なお、図１４上の仮想線ｍａ，ｍｃは速度し
きい値Ｖthを一定としたときの速度変換特性の別の例を
示している。
【００７６】（第３実施例）本発明の第３実施例を図１
６〜図２０に基づいて説明する。この第３実施例に係る
超音波診断装置は前述した組織運動の高速域のブランク
処理を、その低速域強調の速度表示コードを演算した後
に行うようにするものである。なお、この第３実施例以
降の実施例では第２実施例と同一又は同等の構成要素に
は同一符号を用いてその説明を省略又は簡略化する。
【００７７】図１６に第３実施例の超音波診断装置のブ
ロック構成を示す。運動速度解析部１１４とＤＳＣ１１
５との間にはブランク処理器１２５が介挿されている。
ブランク制御器１２１からは前述した係数信号ＳK のみ
が前述した速度演算器１１４ｂに供給される一方で、ブ
ランク処理器１２５に速度表示コードのコードしきい値
信号ＳCDthが供給されるようになっている。
【００７８】ブランク制御器１２１はＣＰＵ１２１０及
びメモリ１２１１を含むコンピュータ機能を有する。Ｃ
ＰＵ１２１０は図１７に示す一連の処理を行う。速度演
算器１１４ｂのＣＰＵ１１４０は本実施例では図１８に
示す一連の処理を行う。さらにブランク処理器１２５も
ＣＰＵ１２５０及びメモリ１２５１を含み、ＣＰＵ１２
５０は図１９に示す一連の処理を行う。
【００７９】最初にブランク制御器１２１の動作を図１
７で説明する。ＣＰＵ１２１０はステップ２５０，２５
１にて入力器１２２から読み込んだ操作信号からスケー
ル変換係数Ｋを算出する。そしてこの算出係数Ｋに対応
した係数信号ＳK を速度演算器１１４ｂに出力する（ス
テップ２５２）。この後ステップ２５３で、血流解析時
に相当する速度変換スケールをＫ倍したときの、最高輝
度の速度表示コード±ＣＤV (MAX) に対応する速度±Ｖ
MAX を演算する（図２９参照）。なお、この±ＶMAX の
演算は速度演算部１１４ｂで行わせて、その演算値を送
り返してもらうようにしてもよい。
【００８０】次いでＣＰＵはステップ２５４で入力器１
２２からの操作信号を読み込み、ステップ２５５で検査
者が欲している所望の速度しきい値Ｖthを算出する（図
２０参照）。ここでの速度しきい値は直接的には偏角θ
に関係しない。というのは、この速度しきい値Ｖthは速
度演算器１１４ｂから出力される、既に変換された速度
表示コードＣＤV に対するしきい値を与えようとするも
のだからである。
【００８１】上記速度しきい値Ｖthを用いて、次のステ
ップ２５６では｜Ｖth｜≦ＶMAX か否かが判断される。
この判断でＹＥＳ、すなわち−ＶMAX ≦Ｖth≦＋ＶMAX
になるときはそのしきい値Ｖthに対応する速度表示コー
ドのしきい値ＣＤth（図２０参照）を演算し（ステップ
２５７）、次いでこのしきい値ＣＤthに対応したコード
しきい値信号ＳCDthをブランク処理器１２５に出力する
（ステップ２５８）。これに対し、ステップ２５６の判
断でＮＯ、すなわち｜Ｖth｜＞ＶMAX となるときは、ス
テップ２５９で、ブランク処理が不可である旨の表示を
ＤＳＣ１１５を介してモニタ１１８に表示させる（図１
６の信号Ｓun参照）。
【００８２】ステップ２５９及び前記ステップ２５８の
処理後、ステップ２６０で再び入力器１２２からの操作
信号の読込みを試み、ステップ２６１で検査者がスケー
ル変換係数Ｋ又は／及び速度しきい値Ｖthの変更を欲し
ているか否かを判断する。この判断でＹＥＳ（変更）の
ときは、前記ステップ２５１又は２５５に戻り上述した
処理を繰り返す。さらにステップ２６１でＮＯ（変更し
ない）となるときは、ステップ２６２で処理終了か否か
を判断し、処理継続の場合ステップ２６０〜２６２の処
理を行いながら待機する。
【００８３】続いて、速度演算器１１４ｂの動作を図１
８に示す。この動作は前述した図１２の一連の処理の中
の対応する符号部分と同一又は同等であって、ステップ
２００，２０１，２０６及び２０７を除いたものであ
る。これにより、ブランク制御器１２１から指令された
スケール変換係数Ｋに応答してリアルタイムに、低速域
を強調した速度表示コードＣＤV がブランク処理器１２
５に出力される。
【００８４】さらにブランク処理器１２５の動作を図１
９に基づき説明する。この処理器１２５のＣＰＵ１２５
０はまず、ステップ２７０でブランク処理器１２１から
供給されるコードしきい値信号ＳCDthを読み込み、ステ
ップ２７１で速度表示コードのしきい値ＣＤthを設定す
る。次いでステップ２７２で速度演算器１１４ｂからの
速度表示コードＣＤV を読み込み、ステップ２７３で｜
ＣＤV ｜＞ＣＤthか否かを判断する。この判断がＹＥＳ
（送られてきたコードＣＤV の絶対値がしきい値ＣＤth
を超えている）のとき、ステップ２７４に移ってその速
度表示コードＣＤV にブランク処理を施す。これにより
｜ＣＤV ｜＞ＣＤthとなる画素の速度Ｖにはブランクコ
ードが強制的に割り当てられる。
【００８５】ステップ２７３でＮＯ、すなわち｜ＣＤV
｜≦ＣＤthとなるとき、及びステップ２７４でブランク
処理を終えた速度表示コード（ブランクコードを含むコ
ード）ＣＤV ＊は、次いで、ＤＳＣ１１５に出力される
（ステップ２７５）。
【００８６】ＣＰＵ１２５０は次いでステップ２７６に
てブランク処理器１２１から供給されているかもしれな
いコードしきい値信号ＳCDthの読込みを試み、ステップ
２７７にてしきい値ＣＤthの変更を欲しているか否かを
判断する。しきい値ＣＤthを変更する場合はステップ２
７１に戻り、変更しない場合はステップ２７８で処理終
了か否かを判断し、終了しないときはステップ２７２に
戻って上述した処理を繰り返す。
【００８７】なお、このブランク処理器１２５は特に図
示していないが、速度の分散値及び／又はパワー値を単
独で観察するときは、これらの値を直接ＤＳＣ１１５に
供給することともに、速度Ｖと併用して観察するとき
は、速度Ｖ（すなわち速度表示コードＣＤV ）を優先
し、｜ＣＤｖ＞ＣＤthとなる画素をブランク処理するよ
うになっている。それらの処理は図示しないコンソール
からの指令信号Ｓcon に基づいて行われる。
【００８８】ブランク制御器１２１，速度演算器１１４
ｂ，ブランク処理器１２５は以上のように動作するの
で、図２０に示す如く、検査者が指令したブランク処理
のしきい値Ｖthが、組織運動の低速度域を強調表示した
ときの最高階調度に対応する速度±Ｖthの間に入ってい
るときのみ、速度しきい値±Ｖthを超える速度の画素が
自動的にブランクされる。これにより前述した第２実施
例と同等の作用効果が得られるほか、スケール変換係数
Ｋを変更する場合でもブランク処理は固定できるという
利点がある。
【００８９】但し、本実施例では図２０で±ＶMAX 以降
の速度表示コードＣＶV には階調性が無いので、検査者
が±ＶMAX を超えるしきい値±Ｖthを指定したとして
も、検査者にその旨の告知があるだけでブランク処理は
なされない。これにより検査者は低速域強調表示におけ
る±ＶMAX が何処にあるかに特別な注意を払わなくても
自由に±Ｖthの値を指令でき、操作が容易になる。
【００９０】なお、このブランク処理器１２５の機能
は、図２１に示す如く、ＤＳＣ１３０に含めてもよい。
この場合、コードしきい値信号ＳCDth及びその他の制御
信号Ｓun，Ｓcon はＤＳＣ１３０に供給される。ＤＳＣ
１３０は前述と同様に速度表示コードＣＤV を速度しき
い値Ｖthに対応するコードしきい値ＣＤthでブランク処
理を行うとともに、組織運動情報の画像データとＢモー
ド像データとを重畳する。これにより、速度演算出力Ｃ
ＤV の階調性を持った低速域に適宜な速度しきい値Ｖth
（ＣＤth）を設定することができ、第３実施例と同等の
ブランク作用及び効果が得られる。
【００９１】（第４実施例）以上説明した第２〜第３実
施例は、何れも低速度強調表示のための速度変換スケー
ルに無関係にブランク処理のしきい値を設定できるもの
であったが、速度変換スケールの設定に連動してしきい
値（すなわちブランク域）を決めることもできる。以下
の実施例ではそのような超音波診断装置を例示する。
【００９２】本発明の第４実施例を図２２〜図２４を参
照して説明する。
【００９３】図２２に示す超音波診断装置は前述した速
度演算器１１４ｂを備えており、そのＣＰＵ１１４０は
図２３に示す一連の処理を行うようになっている。ブラ
ンク制御器１２１からは低速域強調表示用の係数信号Ｓ
K のみが供給されるようになっている。
【００９４】図２３において、速度演算器１１４ｂのＣ
ＰＵ１１４０はステップ３００、３０１の処理（図１２
のステップ２０２、２０３と同じ処理）を行った後、ス
テップ３０２で速度変換スケール（図２４の直線ａ…ｃ
参照）の最高階調度に相当するコードしきい値±ＣＤ
V(MAX)を演算する。その後、ステップ３０３、３０４の
処理（図１２のステップ２０４、２０５と同じ処理）を
行う。この図２３の処理ではステップ３０４で偏角θが
求まると、直ちにステップ３０５、３０６の処理（図１
２のステップ２０８、２０９と同じ処理）を実行する。
次いでＣＰＵ１１４０はステップ３０７に移行し、コー
ドしきい値±ＣＤV(MAX)を用いて、速度表示コードＣＤ
V が｜ＣＤV ｜＜ＣＤV(MAX)か否かを判断する。
【００９５】この判断は最高階調度に対応する速度以上
の組織運動はブランク処理しようとする趣旨に基づくも
ので、したがってＫ倍した後のスケールの最高階調度に
よってブランク処理する速度しきい値Ｖthが自動的に決
まってしまう。換言すれば、低速度強調用の速度変換ス
ケールに連動して速度しきい値が決まることになる。例
えば図２４の例ではスケール変換係数Ｋ＝Ｋ1 （＞１）
のときの速度変換スケールａの場合、速度しきい値±Ｖ
th-1に、Ｋ＝Ｋ2 （＞Ｋ1 ）のときのスケールｂの場
合、速度しきい値±Ｖth-2に、さらにＫ＝Ｋ3 （＞
Ｋ2 ）のときのスケールｃの場合、速度しきい値±Ｖ
th-3に各々、連動して決まる。
【００９６】このステップ３０７の判断でＹＥＳ（｜Ｃ
ＤV ｜＜ＣＤV(MAX)）となるときははステップ３０９に
移行して、そのままのコード値を有する速度表示コード
ＣＤV をＤＳＣ１１５に出力する。しかし、ＮＯ（｜Ｃ
ＤV ｜≧ＣＤV(MAX)）となるとき、すなわち求めた速度
表示コードＣＤV が赤又は青の輝度の最高階調度に達し
ている場合、ステップ３０８で速度表示コードＣＤV の
コード値＝ブランクコードに強制設定する。このように
ブランク処理されたコードＣＤV はステップ３０９でや
はりＤＳＣ１１５に出力される。この結果、例えば±Ｃ
ＤV (MAX) ＝±１２８のとき、ＣＤV ＝±１２７までは
ブランク処理されないが、ＣＤV ＝±１２８のときはブ
ランク処理される。これは図２４に示すように最高階調
度±ＣＤV (MAX) からＣＤV ＝０に立ち下がる特性にな
る。
【００９７】次いでＣＰＵ１１４０はステップ３１０〜
３１２の処理を、図１２のステップ２１１〜２１３と同
じに行う。
【００９８】以上のブランク処理によって、組織運動の
低速域強調表示のための速度変換スケールａ（…ｃ）を
任意に選択／切換できるとともに、各スケールａ（…
ｃ）に連動して決まる速度値Ｖth以上の速度の画素は自
動的にブランクされ、背景像が現われる。
【００９９】したがって、この実施例によっても前述し
た低速域強調表示及びブランク処理の利点を良好に享受
できる。また検査者が入力器１２２から指令するのはス
ケール変換係数Ｋだけで済む利点もある。
【０１００】なお、上記実施例ではブランク処理のしき
い値を±ＣＤV (MAX) の最高階調に設定したが、必要に
応じて、これより低い適宜なコード値に設定することも
でき、その場合にもそのコード値に連動して速度しきい
値が自動的に決まり、それ以上の速度のカラー画素のみ
は同様にブランクされる。
【０１０１】（第５実施例）続いて本発明の第５実施例
を図２５，２６を参照して説明する。
【０１０２】この実施例では図２５に示す如く、組織運
動解析部１１４とＤＳＣ１１５との間に図１６と同様に
別個のブランク処理器１３１を介挿し、図２６に示すス
テップを含む処理をそのＣＰＵ（図示せず）で実行させ
るようにしている。ブランク制御器１２１はここではコ
ードしきい値ＣＤthを表わすしきい値信号ＳCDthブラン
ク処理器１３１に与える。
【０１０３】ブランク処理器１３１は図２６に示すよう
に、コードしきい値信号ＳCDthを読み込み、対応するコ
ードしきい値ＣＤth（ここでは例えば前記実施例と同様
に±ＣＤV (MAX) ）を設定する（ステップ３２０，３２
１）。次いで、演算された速度表示コードＣＤV が読み
込まれ、｜ＣＤV ｜＜ＣＤthか否かが判断される（ステ
ップ３２３，３２４）。｜ＣＤV ｜＝ＣＤthとなるとき
はブランク処理される（ステップ３２５）。
【０１０４】したがって、ブランクする速度範囲を速度
変換スケールに連動させる場合、平均ドプラ速度の表示
コードを演算した後でも、コードしきい値を適宜に決め
てやることで容易に実施でき、図２４と同一の作用効果
が得られる。
【０１０５】なお、この実施例においてコードしきい値
ＳCDthは速度演算器１１４ｂで演算してブランク処理器
１３１に与えるようにしてもよい。
【０１０６】またなお、上述したブランク処理器１３１
の機能を表示系としてのＤＳＣ１１５の中に一体に含め
て実施することもできる（図２１参照；但し、この場合
しきい値設定の可否に関する制御信号Ｓunの供給は不要
である）。
【０１０７】さらに前記第２〜第５実施例及びその変形
例では速度表示コードＣＤV を組織運動の方向を識別す
る赤色及び青色の輝度の階調として表わしたが、速度表
示コードＣＤV はそのほかにも例えば色相のグラデーシ
ョンで表わすようにしてもよい。
【０１０８】
【発明の効果】以上説明したように、本発明の一態様の
超音波診断装置によれば、被検体内の組織の運動速度を
超音波パルスドプラ法を用いて計測し、二次元カラー表
示するためのスケールの勾配を、低速度域でその他の速
度域のそれよりも上げるようにしたため、その低速度域
の速度表示能が高められる。これによって、低いパルス
繰返し周波数にすることに依存する、心筋などの低速度
領域に対する高感度な検出能を損うことなく、僅かな速
度差であっても異なる階調（例えば色の輝度または色合
いの違い）で表示されることから、低速度領域における
速度差を視覚を通して容易に且つ高精度に評価すること
が可能となり、組織ドプライメージングにおける表示機
能を充実させることができる。また、また、例えば、最
高階調度の速度を有する画素は、それ以下の速度の画素
とは不連続な色相で表示されるので、その見分けが容易
で、その輪郭に対する視認性が高くなる。
【０１０９】本発明のほかの態様に係る超音波診断装置
では、心筋などの組織の運動速度がしきい値より高いサ
ンプル点のＴＤＩカラー像は表示されず、背景のＢモー
ド断層像のみが表示されるため、しきい値を適宜に指定
することで、階調度が無く、診断上有効な画像情報を提
供し得ない最高階調度（例えば最高輝度）の速度カラー
像を排除又は必要最小限に抑えることができ、その代わ
りに背景像であるＢモード断層像を表示できるので、そ
の背景像からより有益な画像情報が得られ、関心部位の
正常／異常の見分けが容易なＴＤＩ像を提供することが
できる。
【図面の簡単な説明】
【図１】本発明の第１実施例に係る超音波診断装置のブ
ロック図。
【図２】同実施例における組織、血流などのドプラ偏移
周波数（運動速度）の特性と組織ドプライメージングに
用いるローパスフィルタの特性を示すグラフ。
【図３】同実施例の速度変換スケールの一例を示すグラ
フ。
【図４】同実施例における符号化演算部の処理の概略フ
ローチャート。
【図５】同実施例における心筋の運動方向とカラーバー
の関係を示す図。
【図６】同実施例における心筋の表示例を示す図。
【図７】速度変換スケールの別の例を示す図。
【図８】速度変換スケールのさらに別の例を示す図。
【図９】速度変換スケールのさらに別の例を示す図。
【図１０】本発明の第２実施例に係る超音波診断装置の
ブロック図。
【図１１】同実施例における速度演算器のブロック図。
【図１２】速度演算器の処理例を示すフローチャート。
【図１３】周波数解析器の出力データを説明する図。
【図１４】速度変換スケールの例を示す図。
【図１５】速度しきい値の変化に対応するブランク域の
増減を説明する図。
【図１６】本発明の第３実施例に係る超音波診断装置の
ブロック図。
【図１７】同実施例におけるブランク制御器の処理例を
示すフローチャート。
【図１８】同実施例における速度演算器の処理例を示す
フローチャート。
【図１９】同実施例におけるブランク処理器の処理例を
示すフローチャート。
【図２０】同実施例における速度変換スケールの一例を
示す図。
【図２１】本発明の変形例に係る超音波診断装置のブロ
ック図。
【図２２】本発明の第４実施例に係る超音波診断装置の
ブロック図。
【図２３】同実施例における速度演算器の処理例を示す
フローチャート。
【図２４】速度変換スケールに連動するブランク域の速
度しきい値の例を示す図。
【図２５】本発明の第５実施例に係る超音波診断装置の
部分ブロック図。
【図２６】同実施例のブランク制御器の処理例を部分的
に示すフローチャート。
【図２７】従来技術の一例を示す血流速度解析時の速度
変換スケールの図。
【符号の説明】
１０ 超音波診断装置 １１ 超音波プローブ １４ 操作パネル １５ 超音波送受信部 １９ 表示器 ２０ 位相検波部 ２１ フィルタ部 ２２ 周波数解析部 ２３ 符号化演算部 ２４ ＴＤＩ用ＤＳＣ部 １００ 超音波プローブ １１０ 超音波送受信部 １１１ 位相検波器 １１２ Ａ／Ｄ変換器 １１３ フィルタ １１４ａ 周波数解析器 １１４ｂ 速度演算器 １１５ ＤＳＣ １１６ カラー処理器 １１７ Ｄ／Ａ変換器 １１８ カラーモニタ １２１ ブランク制御器 １２２ 入力器 １２５ ブランク処理器 １３０ ＤＳＣ １３１ ブランク処理器

Claims (27)

    【特許請求の範囲】
  1. 【請求項１】 被検体の断層面に含まれる心筋，血管壁
    などの運動に関するカラー画像を表示する超音波診断装
    置において、 超音波パルス信号を前記断層面に沿って走査して当該断
    層面からの超音波反射信号に応じた電気量のエコー信号
    を得る走査手段と、前記エコー信号から前記組織の運動
    に従うドプラ信号を抽出する抽出手段と、前記ドプラ信
    号に基づいて前記組織の運動の速度データを前記断層面
    のサンプル点毎に求める速度算出手段と、前記超音波パ
    ルス信号のパルス繰返し周波数で制限される測定可能な
    ドプラ周波数範囲の前記速度データと当該速度データに
    割り当てるカラー表示の階調データとをそのドプラ周波
    数範囲内の低速度域をその残りの速度域よりも強調した
    状態で対応付けるスケールを設定するスケール設定手段
    と、前記スケールに基づいて前記速度データを前記階調
    データに変換する速度変換手段と、この速度変換手段に
    より変換された階調データをカラー化して表示する表示
    手段とを備えたことを特徴とする超音波診断装置。
  2. 【請求項２】 前記抽出手段は、前記組織によりドプラ
    偏移を受けた前記ドプラ信号を選択的に抽出するローパ
    スフィルタを備える請求項１記載の超音波診断装置。
  3. 【請求項３】 前記スケールは、前記速度データの変化
    に対する前記階調データの変化の割合が非線形に設定さ
    れている請求項２記載の超音波診断装置。
  4. 【請求項４】 前記スケールは、前記低速度域の前記割
    合を前記残りの速度域の前記割合よりも高く設定されて
    いる請求項３記載の超音波診断装置。
  5. 【請求項５】 前記低速度領域は、前記超音波パルス信
    号のパルス繰返し周波数をｆｒ，ドプラ偏移周波数をｆ
    ｄとしたとき、−ｆｒ／８≦ｆｄ≦ｆｒ／８，−ｆｒ／
    １２≦ｆｄ≦ｆｒ／１２，および−ｆｒ／１６≦ｆｄ≦
    ｆｒ／１６のうちの１つである請求項４記載の超音波診
    断装置。
  6. 【請求項６】 前記スケールは、前記低速度域を前記階
    調データの全てに割り当て且つ前記残りの速度域を全て
    当該階調データの最高値に割り当てたスケールである請
    求項４記載の超音波診断装置。
  7. 【請求項７】 前記スケールは、前記低速度域に対する
    前記割合が直線的に変化する請求項６記載の超音波診断
    装置。
  8. 【請求項８】 前記スケールは、前記低速度域に対する
    前記割合が第１の値及びこの第１の値より低い第２の値
    を有する折れ曲がった直線状に変化するとともに、前記
    第１の値を有する直線部分の変化が前記ドプラ周波数範
    囲内の低速側に位置する請求項６記載の超音波診断装
    置。
  9. 【請求項９】 前記階調データは前記組織の運動方向を
    各別に表わす色の輝度および色合いのうちの一方を変化
    させた複数のカラーデータから成る請求項１記載の超音
    波診断装置。
  10. 【請求項１０】 前記階調データは前記速度データの変
    化に応じて変わる連続的なつながりを有する階調度のカ
    ラーデータを含むとともに、当該階調データの少なくと
    も最高値はその低階調度側とは階調度的に不連続なカラ
    ーデータである請求項９記載の超音波診断装置。
  11. 【請求項１１】 前記抽出手段は、前記組織によりドプ
    ラ偏移を受けた前記ドプラ信号を選択的に抽出するロー
    パスフィルタを備える請求項１０記載の超音波診断装
    置。
  12. 【請求項１２】 前記スケールは、前記速度データの変
    化に対する前記階調データの変化の割合が非線形に設定
    されている請求項１１記載の超音波診断装置。
  13. 【請求項１３】 前記スケールは、前記低速度域の前記
    割合を前記残りの速度域の前記割合よりも高く設定され
    ている請求項１２記載の超音波診断装置。
  14. 【請求項１４】 被検体の断層面の断層像に当該断層面
    に含まれる組織の運動に関するカラー画像を重畳して表
    示する超音波診断装置において、 超音波パルス信号を前記断層面に沿って走査して当該断
    層面からの超音波反射信号に応じた電気量のエコー信号
    を得る走査手段と、前記エコー信号から前記組織の運動
    に従うドプラ信号を抽出する抽出手段と、前記ドプラ信
    号に基づいて前記組織の運動の速度データを前記断層面
    のサンプル点毎に求める速度算出手段と、前記エコー信
    号から前記断層像のデータを形成する断層像形成手段
    と、前記速度データが速度しきい値を超えた場合、その
    速度データをブランク処理するブランク処理手段と、前
    記ブランク処理された前記速度データを前記断層像のデ
    ータに重畳するとともに当該速度データをカラー化して
    表示する表示手段とを備えたことを特徴とする超音波診
    断装置。
  15. 【請求項１５】 被検体の断層面の断層像に当該断層面
    に含まれる組織の運動に関するカラー画像を重畳して表
    示する超音波診断装置において、 超音波パルス信号を前記断層面に沿って走査して当該断
    層面からの超音波反射信号に応じた電気量のエコー信号
    を得る走査手段と、前記エコー信号から前記組織の運動
    に従うドプラ信号を抽出する抽出手段と、前記ドプラ信
    号に基づいて前記組織の運動の速度データを前記断層面
    のサンプル点毎に求める速度算出手段と、前記エコー信
    号から前記断層像のデータを形成する断層像形成手段
    と、前記超音波パルス信号のパルス繰返し周波数で制限
    される測定可能なドプラ周波数に対する速度測定範囲の
    前記速度データと当該速度データに割り当てるカラー表
    示の階調データとをその速度測定範囲内の低速度域をそ
    の残りの速度域よりも強調した状態で対応付けるスケー
    ルを設定するスケール設定手段と、前記スケールに基づ
    き前記速度データを前記階調データに変換する速度変換
    手段と、この速度変換手段により変換された階調データ
    及び前記速度算出手段により算出された前記速度データ
    の何れか一方の量が指定されたしきい値を超える前記サ
    ンプル点をブランク処理に付すブランク処理手段と、前
    記ブランク処理手段によるブランク処理を経た前記速度
    データを前記断層像のデータに重畳するとともに当該速
    度データをカラー化して表示する表示手段とを備えたこ
    とを特徴とする超音波診断装置。
  16. 【請求項１６】 前記スケール設定手段は、前記ドプラ
    周波数の変化に対する前記階調データの変化の割合を前
    記被検体の運動流体解析時の当該割合よりも高く、且つ
    前記階調データの最高値に対応する前記速度データより
    も大きい速度データには当該最高値を割り当てる前記ス
    ケールを設定する手段である請求項１５記載の超音波診
    断装置。
  17. 【請求項１７】 前記階調データは、前記組織の運動方
    向を各別に表わす指定カラーの輝度の階調データから成
    る請求項１６記載の超音波診断装置。
  18. 【請求項１８】 前記階調データは、前記組織の運動方
    向を各別に表わす指定カラーの色相の階調データから成
    る請求項１６記載の超音波診断装置。
  19. 【請求項１９】 前記スケールの設定とは独立して前記
    しきい値を指定するしきい値指定手段を備える請求項１
    ６記載の超音波診断装置。
  20. 【請求項２０】 前記一方の量は、前記速度算出手段に
    より算出された速度データである請求項１９記載の超音
    波診断装置。
  21. 【請求項２１】 前記速度算出手段は、前記ドプラ信号
    の周波数成分を解析する周波数解析手段と、この周波数
    解析後の解析結果に基づいて前記各サンプル点の前記速
    度データを演算する速度演算手段とを有し、この速度演
    算手段、前記速度変換手段及びブランク処理手段は１つ
    の演算ユニットに設けてある請求項２０記載の超音波診
    断装置。
  22. 【請求項２２】 前記しきい値は、前記スケールで決ま
    る前記階調データの最高値に対応した前記速度データの
    値である請求項２０記載の超音波診断装置。
  23. 【請求項２３】 前記一方の量は、前記速度変換手段に
    より変換された階調データであって、前記しきい値指定
    手段は、前記階調データの最高値よりも低い階調データ
    のしきい値を指定する手段である請求項１９記載の超音
    波診断装置。
  24. 【請求項２４】 前記ブランク処理手段は、少なくとも
    前記速度算出手段及び表示手段とは独立した処理器に設
    けてある請求項２３記載の超音波診断装置。
  25. 【請求項２５】 前記表示手段は、前記速度データを前
    記断層像のデータに重畳するデジタルスキャンコンバー
    タを有し、このデジタルスキャンコンバータは前記ブラ
    ンク処理手段を含む請求項２３記載の超音波診断装置。
  26. 【請求項２６】 前記スケール設定手段は、前記スケー
    ルの設定に連動して前記しきい値を自動的に決める手段
    である請求項１６記載の超音波診断装置。
  27. 【請求項２７】 前記一方の量は、前記速度変換手段に
    より変換された階調データである請求項２６記載の超音
    波診断装置。
JP02179595A 1994-06-24 1995-02-09 超音波診断装置 Expired - Lifetime JP3833282B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP02179595A JP3833282B2 (ja) 1994-06-24 1995-02-09 超音波診断装置
US08/415,479 US5513640A (en) 1994-06-24 1995-03-31 Diagnostic ultrasound system
US09/073,876 USRE38209E1 (en) 1994-06-24 1998-05-07 Diagnostic ultrasound system

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP6-143528 1994-06-24
JP14352894 1994-06-24
JP02179595A JP3833282B2 (ja) 1994-06-24 1995-02-09 超音波診断装置

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2006167401A Division JP2006314807A (ja) 1994-06-24 2006-06-16 超音波診断装置

Publications (2)

Publication Number Publication Date
JPH0866397A true JPH0866397A (ja) 1996-03-12
JP3833282B2 JP3833282B2 (ja) 2006-10-11

Family

ID=26358899

Family Applications (1)

Application Number Title Priority Date Filing Date
JP02179595A Expired - Lifetime JP3833282B2 (ja) 1994-06-24 1995-02-09 超音波診断装置

Country Status (2)

Country Link
US (2) US5513640A (ja)
JP (1) JP3833282B2 (ja)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0797106A2 (en) * 1996-03-22 1997-09-24 Advanced Technology Laboratories, Inc. Three dimensional medical ultrasonic diagnostic imaging of tissue texture and vasculature
JP2000189418A (ja) * 1998-12-24 2000-07-11 General Electric Co <Ge> パワ―・ドップラ―・イメ―ジの最適デ―タ・マッピングのための方法および装置
JP2003501177A (ja) * 1999-06-08 2003-01-14 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 同時的組織及びモーション超音波診断撮像
WO2004062503A1 (ja) * 2003-01-15 2004-07-29 Hitachi Medical Corporation 超音波診断装置
WO2005122906A1 (ja) * 2004-06-18 2005-12-29 Hitachi Medical Corporation 超音波診断装置
US7662100B2 (en) 2002-12-20 2010-02-16 Aloka Co., Ltd. Ultrasonic diagnostic device
JP2013135869A (ja) * 2013-02-18 2013-07-11 Toshiba Corp 医用画像処理装置
JP2014144370A (ja) * 2014-04-07 2014-08-14 Hitachi Medical Corp 超音波診断装置

Families Citing this family (53)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5622174A (en) * 1992-10-02 1997-04-22 Kabushiki Kaisha Toshiba Ultrasonic diagnosis apparatus and image displaying system
JPH0975348A (ja) * 1995-09-12 1997-03-25 Toshiba Corp 超音波診断装置
US5793883A (en) * 1995-09-29 1998-08-11 Siemens Medical Systems, Inc. Method for enhancing ultrasound image
FR2741521A1 (fr) * 1995-11-23 1997-05-30 Philips Electronique Lab Procede d'acquisition automatique de l'angle doppler, en echographie, et agencement pour la mise en oeuvre du procede
NO963175D0 (no) * 1996-07-30 1996-07-30 Vingmed Sound As Analyse- og målemetode
US5800356A (en) * 1997-05-29 1998-09-01 Advanced Technology Laboratories, Inc. Ultrasonic diagnostic imaging system with doppler assisted tracking of tissue motion
US5857973A (en) * 1997-09-30 1999-01-12 Siemens Medical Systems, Inc. Fuzzy logic tissue flow determination system
US5882315A (en) * 1997-12-23 1999-03-16 Acuson Corporation Ultrasonic imaging method and image for doppler tissue parameters
US6039692A (en) * 1998-06-19 2000-03-21 Vingmed Sound A/S Method and apparatus for processing ultrasound signals
US6196972B1 (en) 1998-11-11 2001-03-06 Spentech, Inc. Doppler ultrasound method and apparatus for monitoring blood flow
US6358208B1 (en) * 1998-11-21 2002-03-19 Philipp Lang Assessment of cardiovascular performance using ultrasound methods and devices that interrogate interstitial fluid
US6120451A (en) * 1998-12-31 2000-09-19 General Electric Company Ultrasound color flow display optimization by adjustment of threshold
US6176830B1 (en) * 1999-07-27 2001-01-23 Siemens Medical Systems, Inc. Method and system for pre-determining spectral doppler user parameters
US7678048B1 (en) 1999-09-14 2010-03-16 Siemens Medical Solutions Usa, Inc. Medical diagnostic ultrasound system and method
US6561979B1 (en) 1999-09-14 2003-05-13 Acuson Corporation Medical diagnostic ultrasound system and method
US6312381B1 (en) * 1999-09-14 2001-11-06 Acuson Corporation Medical diagnostic ultrasound system and method
US6497664B1 (en) 1999-09-14 2002-12-24 Ecton, Inc. Medical diagnostic ultrasound system and method
US6524244B1 (en) 1999-09-14 2003-02-25 Ecton Inc. Medical diagnostic ultrasound system and method
US6468213B1 (en) 1999-09-14 2002-10-22 Ecton, Inc. Medical diagnostic ultrasound system and method
US6436039B1 (en) 1999-09-14 2002-08-20 Ecton, Inc. Medicial diagnostic ultrasound system and method
US6508763B1 (en) 1999-09-14 2003-01-21 Ecton, Inc. Medical diagnostic ultrasound system and method
JP3495710B2 (ja) * 2001-02-01 2004-02-09 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 血流イメージング装置および超音波診断装置
US6599248B1 (en) * 2001-03-20 2003-07-29 Aloka Method and apparatus for ultrasound diagnostic imaging
US7044913B2 (en) * 2001-06-15 2006-05-16 Kabushiki Kaisha Toshiba Ultrasonic diagnosis apparatus
US6749571B2 (en) 2002-09-19 2004-06-15 Wisconsin Alumni Research Foundation Method and apparatus for cardiac elastography
CA2501647C (en) 2002-10-10 2013-06-18 Visualsonics Inc. High frequency high frame-rate ultrasound imaging system
JP3669588B2 (ja) * 2003-05-06 2005-07-06 学校法人慶應義塾 超音波流速分布計及び流量計、超音波流速分布及び流量測定方法並びに超音波流速分布及び流量測定処理プログラム
JP4800214B2 (ja) * 2003-05-30 2011-10-26 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ カラーフローバイプレーンの超音波撮像システム及び方法
US7128713B2 (en) * 2003-07-10 2006-10-31 Spentech, Inc. Doppler ultrasound method and apparatus for monitoring blood flow and hemodynamics
EP1515158B1 (en) 2003-09-09 2013-07-17 Esaote S.p.A. Ultrasound imaging method combined with the presence of contrast media in the body under examination
JP2007525299A (ja) * 2004-03-01 2007-09-06 サニーブルック アンド ウィメンズ カレッジ ヘルス サイエンシーズ センター Ｅｃｇトリガレトロスペクティブカラーフロー超音波イメージングのためのシステムおよび方法
US7261695B2 (en) * 2004-03-09 2007-08-28 General Electric Company Trigger extraction from ultrasound doppler signals
JP4127827B2 (ja) * 2004-06-30 2008-07-30 オリンパス株式会社 超音波診断装置
US7713210B2 (en) 2004-11-23 2010-05-11 St. Jude Medical, Atrial Fibrillation Division, Inc. Method and apparatus for localizing an ultrasound catheter
US20060122505A1 (en) * 2004-11-23 2006-06-08 Ep Medsystems, Inc. M-Mode presentation of an ultrasound scan
CN101141920B (zh) * 2005-03-15 2011-12-14 株式会社东芝 超声波诊断装置及其控制方法
US7771358B2 (en) * 2005-05-20 2010-08-10 Spentech, Inc. System and method for grading microemboli monitored by a multi-gate doppler ultrasound system
US8162837B2 (en) * 2005-06-13 2012-04-24 Spentech, Inc. Medical doppler ultrasound system for locating and tracking blood flow
JP4920302B2 (ja) * 2005-06-20 2012-04-18 株式会社東芝 超音波診断装置及び超音波計測方法
CN100525711C (zh) * 2005-08-29 2009-08-12 深圳迈瑞生物医疗电子股份有限公司 基于运动插值的解剖m型成像方法和装置
JP4694930B2 (ja) * 2005-09-21 2011-06-08 富士フイルム株式会社 超音波診断装置
JP2007105400A (ja) * 2005-10-17 2007-04-26 Toshiba Corp 超音波診断装置及び画像処理装置
JP4730125B2 (ja) * 2006-02-22 2011-07-20 株式会社日立製作所 血流画像表示装置
KR101055581B1 (ko) * 2008-03-20 2011-08-23 삼성메디슨 주식회사 초음파 시스템 및 클러터 신호 제거용 필터 설정 방법
US8366624B1 (en) * 2008-10-02 2013-02-05 Hitachi Aloka Medical, Ltd. Methods and apparatus for ultrasound imaging
JP2012030048A (ja) * 2010-07-02 2012-02-16 Toshiba Corp 超音波診断装置、画像処理装置、及び解析プログラム
JP2012176232A (ja) * 2011-02-04 2012-09-13 Toshiba Corp 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム
JP5788229B2 (ja) * 2011-06-06 2015-09-30 株式会社東芝 超音波診断装置
KR101484958B1 (ko) * 2012-03-28 2015-01-21 삼성메디슨 주식회사 초음파 시스템 및 초음파 영상 획득 방법
CN105324078B (zh) * 2013-06-28 2019-04-19 皇家飞利浦有限公司 对介入器械的不依赖扫描器的跟踪
EP2989986B1 (en) * 2014-09-01 2019-12-18 Samsung Medison Co., Ltd. Ultrasound diagnosis apparatus and method of operating the same
US20190269380A1 (en) * 2018-03-05 2019-09-05 Edan Instruments, Inc Blood flow displaying method, apparatus and device for ultrasonic system
JP7334494B2 (ja) * 2019-06-25 2023-08-29 コニカミノルタ株式会社 超音波診断装置、超音波診断装置の制御方法、及び、超音波診断装置の制御プログラム

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4800891A (en) * 1987-11-13 1989-01-31 Siemens Medical Laboratories, Inc. Doppler velocity processing method and apparatus
US5170792A (en) * 1989-11-27 1992-12-15 Acoustic Imaging Technologies Corporation Adaptive tissue velocity compensation for ultrasonic Doppler imaging
US5148809A (en) * 1990-02-28 1992-09-22 Asgard Medical Systems, Inc. Method and apparatus for detecting blood vessels and displaying an enhanced video image from an ultrasound scan
US5215094A (en) * 1990-11-14 1993-06-01 Advanced Technology Laboratories, Inc. Ultrasonic flow velocity imaging systems with velocity image presistence
JP2791255B2 (ja) * 1992-10-02 1998-08-27 株式会社東芝 超音波カラードプラ断層装置
US5285788A (en) * 1992-10-16 1994-02-15 Acuson Corporation Ultrasonic tissue imaging method and apparatus with doppler velocity and acceleration processing

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0797106A3 (en) * 1996-03-22 1999-09-15 Advanced Technology Laboratories, Inc. Three dimensional medical ultrasonic diagnostic imaging of tissue texture and vasculature
EP0797106A2 (en) * 1996-03-22 1997-09-24 Advanced Technology Laboratories, Inc. Three dimensional medical ultrasonic diagnostic imaging of tissue texture and vasculature
JP4712931B2 (ja) * 1998-12-24 2011-06-29 ゼネラル・エレクトリック・カンパニイ パワー・ドップラー・イメージの最適データ・マッピングのための方法および装置
JP2000189418A (ja) * 1998-12-24 2000-07-11 General Electric Co <Ge> パワ―・ドップラ―・イメ―ジの最適デ―タ・マッピングのための方法および装置
JP2003501177A (ja) * 1999-06-08 2003-01-14 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 同時的組織及びモーション超音波診断撮像
US7662100B2 (en) 2002-12-20 2010-02-16 Aloka Co., Ltd. Ultrasonic diagnostic device
US7628754B2 (en) 2003-01-15 2009-12-08 Hitachi Medical Corporation Ultrasonographic device
WO2004062503A1 (ja) * 2003-01-15 2004-07-29 Hitachi Medical Corporation 超音波診断装置
US8206298B2 (en) 2003-01-15 2012-06-26 Hitachi Medical Corporation Ultrasonographic elasticity imaging device
WO2005122906A1 (ja) * 2004-06-18 2005-12-29 Hitachi Medical Corporation 超音波診断装置
JPWO2005122906A1 (ja) * 2004-06-18 2008-04-10 株式会社日立メディコ 超音波診断装置
JP5203605B2 (ja) * 2004-06-18 2013-06-05 株式会社日立メディコ 超音波診断装置
JP2013135869A (ja) * 2013-02-18 2013-07-11 Toshiba Corp 医用画像処理装置
JP2014144370A (ja) * 2014-04-07 2014-08-14 Hitachi Medical Corp 超音波診断装置

Also Published As

Publication number Publication date
USRE38209E1 (en) 2003-08-05
US5513640A (en) 1996-05-07
JP3833282B2 (ja) 2006-10-11

Similar Documents

Publication Publication Date Title
JP3833282B2 (ja) 超音波診断装置
JP2791255B2 (ja) 超音波カラードプラ断層装置
US6884216B2 (en) Ultrasound diagnosis apparatus and ultrasound image display method and apparatus
JP5100343B2 (ja) 超音波診断装置、及び超音波診断装置の制御プログラム
US7871379B2 (en) Ultrasonic diagnostic apparatus and method of ultrasonic measurement
US6126605A (en) Ultrasound color flow display optimization by adjusting dynamic range
JP2006520619A (ja) 超音波カラーフロードプラ情報から合成されるフロースペクトログラム
US6814703B2 (en) Apparatus and method for ultrasonic diagnostic imaging using a contrast medium
US20060170714A1 (en) Ultrasound diagnosis apparatus and ultrasound data generating method
US4785402A (en) Ultrasonic imaging apparatus for color display of flow velocity
EP0349321B1 (en) Method and apparatus for determining and displaying the absolute value of quantitative backscatter
US5664571A (en) Ultrasonic diagnostic apparatus
JP3253409B2 (ja) 超音波ドプラ診断装置
JP2004313291A (ja) 超音波診断装置、医用画像解析装置及び医用画像解析方法
JP3187008B2 (ja) 超音波カラードプラ断層装置
US20120203111A1 (en) Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, and ultrasonic image acquisition method
US20070073152A1 (en) Systems and methods for acquiring images simultaneously
US6500125B1 (en) Ultrasound b/color priority threshold calculation
JP4304217B2 (ja) 超音波診断装置
US7371219B2 (en) Ultrasound diagnosis apparatus operable in doppler mode
JP5455567B2 (ja) 超音波診断装置
JP2006314807A (ja) 超音波診断装置
JP7343342B2 (ja) 超音波診断装置、及び画像処理装置
JP3263171B2 (ja) 超音波診断装置
JP3405578B2 (ja) 超音波診断装置

Legal Events

Date Code Title Description
A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060616

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20060719

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

Free format text: PAYMENT UNTIL: 20100728

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20100728

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20110728

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20120728

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20130728

Year of fee payment: 7

EXPY Cancellation because of completion of term