JPH0417842A

JPH0417842A - 超音波診断装置

Info

Publication number: JPH0417842A
Application number: JP12236090A
Authority: JP
Inventors: Kiyoshi Nakayama; 中山　淑
Original assignee: Fujitsu Ltd
Current assignee: Fujitsu Ltd
Priority date: 1990-05-11
Filing date: 1990-05-11
Publication date: 1992-01-22
Anticipated expiration: 2014-01-20
Also published as: JP2849159B2

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】〔概要〕生体組織の弾性特性を算出表示するための手段を備えた
超音波診断装置に関し、パルスドプラ法によって得られた各深さにおけるｍ織の
変位速度を深さ方向に微分することによって、組織の伸
び縮みの度合いを示す「ずり速度」を算出表示すること
を目的とし、被検体内の同一方向に複数回超音波パルスを送信し、該
同一方向から受信された複数本を受信信号を使って、各
深さに於ける被検体内組織の変位速度を検出する速度検
出手段と、該検出した速度を深さ方向に微分する深さ方
向微分手段と、前記深さ方向微分手段によって得られた
ずり速度の実効値を算出する実効値算出手段と、前記深
さ方向微分手段シこよって得られたずり速度を時間と深
さの二次元関数、又は、空間的な二次元分布像として、
また、前記実効値算出手段によって得られた結果を数値
として表示するための表示手段を有する構成とする。

［産業上の利用分野］本発明は、パルスドプラ法によって得ろれた各深さ乙こ
おける組織の変位速度を深さ方向に微分すること乙こよ
って、組織の伸び縮みの度合いを示す弾性的徂識特性を
算出表示する超音波診断装置に関するものである。

〔従来の技術と発明が解決しようとする課題］従来、組
織の弾性的特性を表示する超音波診断装置として、第７
図のように生体外部から低周波振動を加え、生体内部を
伝播する振動波の振幅と位相分布を計測しようとするも
の（第５２回日本超音波医学会論文集５２−１４３　　
Ｐ２Ｓ５−２８８、及び、第５３回日本超音波医学会論
文集５３−８３．Ｐ２７１−２７２）や、第８図のよう
に解析信号の空間相関関数を用いた不拘−Ｍｉ織の微小
変位の計測をしようとするもの（第５４回日本超音波医
学会論文集５４−１１６．Ｐ２Ｓ５−３６０）が提案さ
れている。

第７図は、低周波振動を加えた時に、生体内部を伝播す
る振動波の振動振幅と位相分布の計測するためのプロー
ブ位置と機械的振動子の位置関係を示したものである。

機械的振動子２ｏは、平板２１を通して、媒体１９に周
波数１０ｋＨｚがら１ｋＨｚ程度を振動を加える。また
、プローブ１１から送波された超音波パルスは媒体内部
で反射し、同プローブ１１で受信される。この受信信号
を使って、振動波２２の振動振幅と位相分布を求めよう
とするのが、第７図である。２３は、プローブと体表と
の間のコンタクト剤、２４は、機械的振動子２０の振動
を吸収するためのスプリングである。

第８図は、解析信号の二次元相関関数を用いた組織変位
を計測する方法である。第８図において１１は、超音波
を送受信するだめの超音波プローブで、送受信回路１０
につながっている。１゜は、ビームフォーマ等を含む送
受信回路で、プロワ１１の構成要素である振動子群を駆
動するための１言号を生成したり、各振動子で受信され
た信号から特定方向にフォーカスした受信信号１を生成
する。受信信号Ｉは、次段の解析信号生成手段１２二こ
送られる。１２では、受信信号をｘｒ（ｔ）とすると、
ｘｒ（ｔ）と直交する信号ｘｉ（ｔ）を生成し、複素信
号Ｚ　（ｔ）　＝ｘｒ（ｔ）　−１−ｊ　ｘｉ（ｔ）を
メモリ１３乙こ格納する。１３は、全ての走査線乙こ対
応する受信複素信号ｚ（ｔ）を複数フレーム分層上格納
する。

ここで、あるフレームに於ける、走査線ｍ、深さｎ！こ
左コける二次元複素信号データをＡ　（ｍ、ｎ）＝Ａｒ
（ｍ、ｎ）−＝　ｊＡｉ（ｍ、ｎ）式（１）と表現する。また、別のフレームに於ける、走査線ｍ、
深さｎ４こ於ける二次元複素信号データをＢ　（ｍ、ｎ
）＝Ｂｒ（ｍ、ｎ）＋ｊＢｉ（ｍ、ｎ）式（２）と表現すると、（ｍ、ｎ）の周りでの二次元相互相関関
数Ｃ（τ、ρ１ｍ、ｎ）は次式で定義できる。

Ｃ（τ、ρ１ｍ、ｎ）＝＜Ａ（ｍ、ｎ）Ｂ”　　（ｍ＋ｒ、ｎ＋ρ）＞［＜　　
　ｌ　Ａ（ｍ、ｎ）　　１２　＞＜　　ｌ　Ｂ（ｍ＋ｒ
、ｎ＋　　　ρ　）　　　ｌ　　２　＞　　コ　１″式
（３）ここで、上式中〈　〉は、二次元空間的な平均操作を示
す。上式の最大値を計算し、複素空間で、（１，０）か
らの距離を計算すると変位ベクトルを計算できる。

１４は、（３）式の二次元相互相関関数を計算するため
のもので、１５は、任意の位置で算出した二次元相互相
関関数を格納するためのメモリである。微小変位検出手
段１６は、変位ベクトルを計算し、結果を表示手段１７
に送る。各位置で計算された変位ベクトルは、第８図（
ｂ）に示すごとく、ベクトルの傾きを移動方向、ベクト
ルの長さを移動量として表示する。

しか巳、第７図の方法においては、外部から振動を与え
るだめの機械系か必要になり、実用的でないこと、また
、第８図の方法では、二次元相互相関関数を計算するの
に十分に広い領域を必要とするため高い空間分解能を得
ることが難しく、計算処理Ｓ工時間がかかるという問題
があった。

ところで、生体内組織では、拍動等の動きによって組織
か変位している。この組織の変位速度は、場所Ｓこよっ
て硬さが異なるために、−様ではなく、ばらついている
。変位の速度の場所によるばらつき（よ、組織の伸び縮
みを意味する。例えば、ある深さ；こ注目し、この深さ
よりプローブ寄りのところで、プローブに向かう変位速
度が、その深さより４Ｆｌ−た位置での変位速度より速
ければ、その深さては、組織が伸びていることを意味す
る。

もし、正常組織よりも硬い異常組織が、正常組織の中に
存在した場合、硬い異常組織内部での変位速度の場所に
よるばらつきは、正常組織のそれに比べて小さいと考え
られる。

本発明は、パルスドプラ法を用いて各場所の組織の変位
速度を検出し、場所による変位速度の変化の度合いを算
出することによって、弾性的特性を表示することを目的
としている。

［課題を解決するための手段］上記課題を解決するための具体的手段を第２図と第３図
を使って説明する。

従来、パルスドプラ法を使って血流の速度を検出する方
法が知られている。本発明においては、組織の変位速度
を検出するのにこの手法を用いる。第３図中２は、この
検出方法を説明するブロンクである。

第３図において、２−１〜２−６は直交検波器を構成す
る。従って、ローパスフィルタ（ＬＰＦ）２−５の出力
は、受信信号１の直交検波出力の実部成分である。また
、ローパスフィルタ（ＬＰＦ）２−６の出力は、受信信
号１の直交検波出力の虚部成分である。

同一走査線方向に、周期Ｔの間隔で８回超音波パルスを
送波した時の受信信号の直交検波出力の実部成分を第２
図（ａ）中、Ｒ，、Ｒ２，・・・・・・、Ｒ９で示す。

同様に、虚部成分を第２図（ａ）中、ＩＩ２．・・・・
・・２１Ｎで示す。これらの信号系列は、それぞれメモ
リ２−９．２−１０に一旦格納される。こ二て、パルス
ドプラ法を使うと、ある深さ２４こだ：するｌ采さ方向
の変位速度は、以下のように算出することができる。

第２図（ａ）中の点線矢印で示すように、深さＺにおい
て時間方向に読みたりた複素信号系列Ｖ、（Ｚ）をｖ、、（ｚ　）　−（Ｒ＋（ｚ）、・・・・・・　ＲＮ
−、（ｚ））−Ｊ　　（１１（ｚ）、　・・・・・・、
ｌＮ−１（１））　　弐（４）と表現する。また、周期
Ｔだけずらして、深さＺにおいて時間方向に読みだした
複素信号系列Ｖ、。

（ｚ）をＶｒ、−１（ｚ　）　−（Ｒ２（Ｚ）、・・・・・・　
Ｒ＋（ｚ））Ｊ　　（ｌｚ（ｚ）、＝−、ＩＮ（ｚ）　
）　　式（５）と表現する。

第３図中、複素自己相関器２−１１は、上記２つの複素
信号系列の相関係数を算出するものである。即ち、相関
係数Ｙ、（ｚ）は、Ｙｎ　　（Ｚ　）　＝Ｖ、、、（ｚ　）　ｖ、　”　（
Ｚ）＝ｌＹｎ　　（ｚ）ｌｅｘｐ　（ｊΔθ、（２））
　　式（６）ミＹＩＩＲ（ｚ　）　＋　Ｊ　Ｙｎｒ　（
ｚ　）　　　　　　　　式（７）となる。ここで、°は
複素共役、Ｙ−＊（ｚ）は相関係数の実数成分、Ｙｎｌ
（Ｚ）は相関係数の虚数成分、Δθ、、（ｚ）は、ドプ
ラ効果によって生じた位相差である。

上式より、ドプラシフト周波数ｆｄ、、（Ｚ）は、ｆ　
ｄｎ　（ｚ　）　＝　ｔ　ａ　ｎ−’　［Ｙｎｌ（ｚ）
／Ｙｆｉ＊（ｚ）　］／Ｔ／（２π）＝Δθ、、（ｚ）／Ｔ／（２π）　　式（８）となる。

従って、組織の深さ方向の変位速度は、ｖ、、（ｚ）　
−Ｌｆｉ（ｚ）Ｃ／　（２ｆｏ　）　　　式（９）とな
る。但し、ｆｏは、送信超音波パルスの中心周波数、Ｃ
は音速である。

第３図中、速度計算部２−１２は、（９）式を計算する
ためのものである。全ての深さについて計算した組織の
深さ方向の変位速度の概念図を第２図（ｂ）に示す。

組織の伸び縮みの度合いを示す弾性的組織パラメータ：
よ、変位速度の場所による変化の度合いで表現できるか
ら、第２図（ｂ）の変位速度を深さ方向二こ微分するこ
とによって得ることができる。我々は、こｎをＭｉ礒の
伸び縮みを表す「ずり速度」５（ｚ）として以下のよう
に定義する。

Ｓ、、（ｚ）−ｄｖ、（ｚ）／ｄｚ　　　　　弐〇〇）
第２１ｍ（ｂ）の変位速度ｖ、（ｚ）を（１０）弐のよ
う二こ深さ方向に微分したずり速度Ｓ、、（ｚ）を概念
的二こ第２図（Ｃ）に示す。

同一走査線方向に送波する回数を増やしく即ち、Ｎ回収
玉にする）、ずり速度Ｓ、（ｚ）をｎｎ’−１，ｎ÷２
．、、、について繰り返し算出し、ある走査線；こつい
ての時間と深さの二次元関数として、ｓ、、（Ｚ）、Ｓ
ｒ＋。、（ｚ）、＝　　を表示することかできる。また
、走査線をスキャンすることによって、各走査線に対し
て得られたＳｎ　（ｚ）を空間的二次元分布像を表示す
ることができる。

更に、従来の超音波診断装置のＢモード像、Ｍモード像
と組み合わせて表示することによって、対象部位の伸び
縮みの様子を的確に判断することができるようになる。

また、ある走査線において、深さｚ＝ｊΔＺｊ＝０．１
，２．・・・・・・（ΔＺは深さ方向のサンプル間隔）
について計算された（Ｓｎ　（Ｊ）　、　　Ｓ、。

ｊ）、・・・　）について、深さの範囲ＣＺＩ　、Ｚｚ
コ、時間の範囲［ｔ＋、ｔｚ］を指定すると、ずり速度
の実効値Ｐは次式で定義できる。

ｌ　　　　　　　　ｎｚ　　　ＪｚＴ（ｒ＋ｚ−ｎ、＋１）ΔＺ　（ｊｚ−ｊ＋＝１）　ｎ
＝ｎ＋　　Ｊ”　Ｊ＋（Ｓｎ　（ｊ）　　　ａ　ｖ）　
２）　””　　　　式（１１）ａｖは、深さの範囲［ｚ
＋　、Ｚｚ　］、時間の範囲［ｔ、、ｔ２ｉ内の平均値
で、 ■ ｎｚ　　　Ｊ２ Σ　　ΣＳ、（ｊ）Ｔ（ｎｚ−ｎ＋＋１）ΔＺ　（ｊｚ−Ｊ＋＋１）ｎ＝ｎ
＋　　Ｊ＝　Ｌ載面となる。また、簡単のためのａ　Ｖ＝Ｏとして計算して
も良い。但し、ｔ＋　＝Ｔｎ＋　　　Ｌｘ　＝Ｔｎｚｚ
Ｉ−ＪｌΔＺ　％　Ｚ　２−Ｊ　ｚΔ２である。

（１１）弐のようなずり速度の実効値は、組織の弾性的
特性の１つのパラメータとなる。また、（１１）式にお
いて、特定の深さＺについてのみ、計算りでも良いこと
はいうまでもない。

二作用：本発明二こよると、組織の弾性的特性であるずり速度を
、ある走査線うこおける時間と深さの二次元的関数、あ
る関心領域におｉする空間的な二次元分布像、或いは、
実効値として、従来の超音波診断装置のＢモード像、Ｍ
モード像と組み合わせて表示されるようになる。これに
よって、Ｍｉ織の位置と硬さの度合いの情報を提供でき
るようになる。

〔実施例］次に、第１図〜第６図を使って本発明を順次説明する。

第１図に、本発明の構成を示す。１１は、媒体１９に対
して、超音波を送受信するための超音波プローブで、送
受信回路１０につながっている。

１０は、ビームフォーマ等を含む送受信回路で、プロー
ブ１１の構成要素である振動子群を駆動するための信号
杏生成したり、各振動子で受信された信号から特定方向
にフォーカスした受信信号１を生成する。受信信号１は
、速度検出手段２、Ｂモード像生成手段８、Ｍモード像
生成手段９の入力になる。１０．１１．２．８及び９は
、従来の超音波診断装置に備わっているものである。速
度検出手段２は、（９）式を算出するためのもので、直
交検波回路等から構成され、任意の深さに於ける深さ方
向の組織の変位速度を検出する。メモＩ７３は、指定し
た走査線における深さ方向の組織の変位速度、或いは、
ある関心領域内の複数走査線における組織の変位速度を
格納するためのものである。深さ方向微分手段４は、（
１０）式を算出するためのもので、手段２によって算出
した変位速度を深さ方向乙こ微分する。実効値算出手段
５は、（１１）式を算出するためのもので、手段４によ
って算出したずり速度の実効値を算出する。

表示手段６は、指定した走査線におけるずり速度の時間
と深さの二次元的関数、走査線をスキャンすることによ
って得られるある関心領域における空間的なずり速度の
二次元分布像、または、手段５によって算出したずり速
度の実効値を従来の超音波診断装置のＢモード像、Ｍモ
ード像と組み合わせて表示するためのものである。走査
線またはＲＯ［指定手段７は、ずり速度を計算するため
の走査線、関心領域を指定するためのＲｏｔ、ずり速度
の実効値を算出するための深さの範囲や時間区間を指定
するための手段である。　第２図は、第１図の本発明を
概念的に説明するものである。

受信信号１は、速度検出手段２内において先ず直交検波
される。ここで、第２図を使って、手段２の動作を説明
する。

第２図（ａ）は、同一方向にＮ回送受信した時の直交検
波された複素信号系列を示している。ここで、複素信号
系列の内、ある深さＺに於ける複素信号系列Ｖ、、（Ｚ
）　　（（４）式）を得る。次に、送波周期Ｔずれた、
同じ深さＺに於ける信号複素信号系列Ｖ、−１（ｚ　）
　　（（５）式）を得る。次に、（６）式に従って、複
素自己相関係数を計算し、（９）式から、深さＺにおけ
る組織の変位速度ｖ７（ｚ）を算出する。全深さについ
て、Ｖ、（ｚ）を計算したのが、図２図（ｂ）である。

次に、深さ方向微分手段４においは、−旦メモリ３に蓄
えられた組織の変位速度ｖ、、（ｚ）を（１０）式に従
って深さ方向に微分する。第２図（Ｃ）は、この時のず
り速度を概念的に示したものである。

第３図は、速度検出手段２及び深さ方向微分手段４をよ
り詳細に説明したものである。　第３図において、速度
検出手段２は、直交検波器２−１〜２−６、直交検波信
号の実部成分をＡＤ変換するためのＡＤＣ２−７、直交
検波信号の虚部成分をＡＤ変換するためのＡＤＣ２−８
、実数成分を格納するメモリ２−９、虚数成分を格納す
るメモＩＪ２−１０、複素自己相関器２−１１、速度計
算部２−１２から構成される。２−１は、受信信号１に
参照信号（ｃｏｓωＱ　ｔ、）　２’−３を乗算するた
めのミキサ、２−５は、前記ミキサ２−１の出力の上側
波帯成分をカットするためのローパスフィルタである。

同様↓こ、２−２は、受信信号１に参照１３号１ニーｓ
ｉｎω。Ｌ）２−４を乗算するためのミキサ、２　６Ｂ
；、前記ミキサ２−２の出力の上１１’ｌＪ　ｅ；、帯
成分をカットするだめのローパスフィルタである。２−
５及び２−６の出力は、直交検波信号の実数部と虚数部
に対応する。直交検波信号の実数部と虚数部は、−旦、
メモリ２−９．２−１０に格納される。２−１１は、同
−深さの複素信号系列に対して、複素自己相関係数を（
７）弐に従って計算するものである。また、２−１１に
おいて算出された複素自己相関係数から組織の変位速度
ｖ、（ｚ）が（９）式に従って算出される。

ここで、血流速度を検出する場合、複素自己相開蓋２−
１１の前段にＭＴＩフィルタを入れるのが通常であるが
、本発明では必要としない。

算出した組織の変位速度は一旦メモリ３に格納される。

−旦格納した変位速度は、深さ方向の時系列として読み
だされ、ＦＩＲフィルタ４−１によって、深さ方向の微
分処理が行われる。ＦＩＲフィルタの係数列の１例とし
て第３図（ｂ）のようにすると、（１０）式に相当する
演算をすることができる。この係数を第３図（Ｃ）のＦ
ＩＲフィルタの係数列４−１−２−０〜４−１−２−ｍ
とする。

第３図（Ｃ）において、データシフトレジスタ４−１１
は、人力の変位速度の信号系列ｖ、、（ｉ）　　（ｉ＝
０．１，２．＝・・・・）を図示しないシステムクロッ
クに従って１段ずつシフトするもので、シフトされた信
号系列（ｖ、（ｉ）、Ｖｆｉ　（ｉ−１）Ｖ、、（ｉ−
２：＋、・・・・・・）にそれぞれ係数（ａ。

”ｌ１１−１＋・・・・・・１を乗算した結果を加算器
４−１３にて加算する。加算結果は、変位速度Ｖ。（ｉ
）を深さ方向に微分した結果に相当し、ずり速度Ｓ、（
ｊ）となる。深さ方向のデータ位置を示すｊは例えばｊ
＝ｉ　−ｍ／２である。

第４図は、ある走査線について算出したずり速度の時間
と深さの二次元関数の表示例を示したものである。

第４図（ａ）は、従来の超音波診断装置で得られるＢモ
ード”像である。ここで、ずり速度を算出するために、
走査線またはＲＯＴ指定手段７によって指定した走査線
方向を示す線を表示手段６によってＢモード像上に重ね
て表示する。この走査線に於けるずり速度のデータセン
ト（Ｓｎ　（Ｚ）、Ｓ、、、（Ｚ）、・・・）を深さ２
方向、時間ｎ方向の関数として表示した例が、第４図（
ｂ）、及び（Ｃ）である。（ｂ）は、ずり速度の値を振
幅として、また、（Ｃ）はずり速度をグレースケール或
いはカラー画像として表示したもので、同図（ｄ）は、
（Ｃ）のずり速度値に対応するクレースケール、または
、カラーを表示したものである。また、第４図（ｅ）の
ように、対応するＭモード像を同時に表示してもよい。

ここで、第４図（ｂ）内に示すように、走査線またはＲ
ＯＩ指定手段７によって、深さの範囲［ｚ。

ｌ２１、時間の範囲［ｔ＋、ｔｚ］を指定する（１１）
弐に従って、ずり速度の実効値を手段５によって算出す
ることができる。

第５図は、肝臓のずり速度を第４図（ｂ）の如く、時間
と深さの二次元関数として表示した例である。第５図（
ａ）は、正常肝臓例、第５図中）は、深さへの位置に転
移癌（転移した腫瘍）がある場合である。正常な肝臓の
例（ａ）では、各深さにおいて、ずり速度が時間的に変
化していく様子が観測できる。また、転移癌（ｂ）の場
合は、腫瘍内部の深さＡにおけるずり速度の時間的変化
が非常に小さくなる。これは、腫瘍内部が硬く、伸び縮
みが小さいことを示しているもの考えられる。更に、腫
瘍と正常組織との境界に当たる深さＢでは、速度変化が
非常に大きいことを示している。このように、本発明に
よって、ずり速度を表示すると、正常組織と腫瘍組織の
違いを明瞭に判断することができる第６図は、従来の超
音波診断装置で得られるＢモード像上に、ずり速度の空
間的二次元分布像を重ねて表示したものである。ずり速
度を算出する範囲シよ、図中、点線のＲＣＩＩで示した
。第６図の例で２よ、Ｂモード像に重ねて表示したが、
Ｂモード像と、ずり速度の空間的二次元像を分けて表示
しても良いことは言うまでもない。

［発明の効果〕以上説明したように、本発明によれば、組織の弾ｉ生的
特性であるずり速度を、ある走査線における時間と深さ
の二次元的関数、ある関心領域における空間的な二次元
分布像、或いは、実効値として、従来の超音波診断装置
のＢモード像、Ｍモード像と組み合わせて表示されるよ
うになる。これ乙こよ−で、′Ｍｉ織の位置と硬さの度
合いの情報を提供できるようになる。

【図面の簡単な説明】

第１図は、本発明の１実施例、第２図は、本発明の１実
施例を概念的に説明する図、第３図は、第１図の構成を
より詳細に説明する図、第４図は、表示例、第５図は、
正常肝臓と転移肝癌に本手法を適用した場合の例、第６
図は、表示の別の例、第７図及び第８図は従来の実施例
である。受信信号速度検出手段メモリ深さ方向微分手段実効値算出手段表示手段Ｒｏｔ指定手段Ｂモード像生成手段Ｍモード像生成手段・超音波プローブ・送受信回路深さ２第２図第１図第３図第４図第６図（ｂ）第５図第７図第８図

Claims

【特許請求の範囲】

（１）被検体内の同一方向に複数回超音波パルスを送信
し、該同一方向から受信された複数本の受信信号から、
各深さに於ける被検体内組織の変位速度を検出する速度
検出手段（２）と、該検出した速度を深さ方向に微分し、ずり速度を算出す
る深さ方向微分手段（４）と、前記深さ方向微分手段（４）によって得られたずり速度
を時間と深さの二次元関数、又は、空間的な二次元分布
像として、表示するための表示手段（６）を有すること
を特徴とする超音波診断装置。
（２）前記深さ方向微分手段（４）によって得られたず
り速度の空間的な二次元分布像をＢモード像と同時に同
一画面上に、又は、Ｂモード像に重ねて表示することを
特徴とする特許請求の範囲第１項記載の超音波診断装置
。
（３）前記深さ方向微分手段（４）によって得られたず
り速度を時間と深さの二次元関数として、Ｂモード像、
或いは、Ｍモード像と同時に表示することを特徴とする
特許請求の範囲第一項記載の超音波診断装置。
（４）走査線またはＲＯＩ指定手段（７）によって、走
査線を指定し、前記深さ方向微分手段は、前記走査線又はＲＯＩ指定手
段（７）によって指定された超音波走査線上で、ずり速
度を算出し、前記算出されたずり速度を時間と深さの二
次元関数として表示することを特徴とする特許請求の範
囲第一項記載の超音波診断装置。
（５）走査線またはＲＯＩ指定手段（７）によって、二
次元的な関心領域（ＲＯＩ）を指定し、該関心領域内の
上記手段（４）によって算出されたずり速度を空間的な
二次元分布像として表示することを特徴とする特許請求
の範囲第一項記載の超音波診断装置。
（６）前記深さ方向微分手段（４）によって得られたず
り速度のうち、指定された時間範囲と深さ範囲のずり速
度において、実効値を算出する実効値算出手段（５）を
有することを特徴とする特許請求の範囲第一項記載の超
音波診断装置。