JPH0421493B2 - - Google Patents

Info

Publication number
JPH0421493B2
JPH0421493B2 JP16859983A JP16859983A JPH0421493B2 JP H0421493 B2 JPH0421493 B2 JP H0421493B2 JP 16859983 A JP16859983 A JP 16859983A JP 16859983 A JP16859983 A JP 16859983A JP H0421493 B2 JPH0421493 B2 JP H0421493B2
Authority
JP
Japan
Prior art keywords
wave
ultrasonic
measurement
probe
burst
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
JP16859983A
Other languages
Japanese (ja)
Other versions
JPS6060837A (en
Inventor
Hirokuni Sato
Nobuyuki Ichida
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.)
Fujitsu Ltd
Original Assignee
Fujitsu 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 Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to JP16859983A priority Critical patent/JPS6060837A/en
Publication of JPS6060837A publication Critical patent/JPS6060837A/en
Publication of JPH0421493B2 publication Critical patent/JPH0421493B2/ja
Granted legal-status Critical Current

Links

Description

【発明の詳細な説明】 本発明は、反射型非線形パラメータ分布測定装
置に係り、測定用超音波バースト波とポンピング
用超音波連続波平面波とを、超音波媒体内で交差
するように送波し、その交差領域から反射される
超音波バースト波の特定の周波数成分から、該交
差領域の等価非線形パラメータの変化量を検知
し、この該非線形パラメータの分布を二次元的あ
るいは三次元的に表示し得る装置に関する。
DETAILED DESCRIPTION OF THE INVENTION The present invention relates to a reflection type nonlinear parameter distribution measuring device, which transmits an ultrasonic burst wave for measurement and an ultrasonic continuous wave plane wave for pumping so as to intersect within an ultrasonic medium. , the amount of change in the equivalent nonlinear parameter of the intersection area is detected from a specific frequency component of the ultrasonic burst wave reflected from the intersection area, and the distribution of this nonlinear parameter is displayed two-dimensionally or three-dimensionally. Regarding the apparatus for obtaining.

従来の反射エコー法による超音波診断装置で
は、生体組織内の音響インピーダンスの変化に関
する情報、つまり肝臓等の組織の輪郭は描出でき
るが、組織の性質、例えば正常な肝臓なのか、悪
性腫瘍に侵された肝臓なのかという区別を行なう
のは困難であつた。このような組織の性質を識別
しようとする方法を組織特性弁別(Tissue
Characterization)と総称し、例えば超音波の非
線形パラメータが組織の性質により異なることか
ら、この非線形パラメータを測定することによ
り、組織の性質を識別しようという方法が考案さ
れている(特願昭58−39907号)。
Conventional ultrasonic diagnostic equipment using reflection echo methods can provide information about changes in acoustic impedance within living tissues, in other words the contours of tissues such as the liver, but they cannot provide information on the nature of the tissue, such as whether it is a normal liver or whether it has been invaded by a malignant tumor. It was difficult to distinguish whether the liver had been treated or not. Tissue characteristic discrimination is a method that attempts to identify the properties of such tissues.
For example, since the nonlinear parameters of ultrasound differ depending on the nature of the tissue, a method has been devised to identify the nature of the tissue by measuring these nonlinear parameters (Japanese Patent Application No. 58-39907 issue).

しかし、従来の非線形パラメータ測定装置は、
透過型のものはあるが、透過型は、送波用の探触
子と受波用の探触子をビーム軸を合わせて対向さ
せねばならず、超音波を透過しにくい骨や体腔内
の空隙また組織内の屈折などにより制限を受ける
という欠点を持つている。
However, conventional nonlinear parameter measurement devices
There is a transmission type, but the transmission type requires a transmitting probe and a receiving probe to face each other with their beam axes aligned. It has the disadvantage of being limited by voids and refraction within the tissue.

本発明の目的は、測定用超音波バースト波とポ
ンピング用超音波連続波を用い、かつこのポンピ
ング波を、被観察部位をカバーしてかつ測定用バ
ースト波に交差する平面波に近いものとするよう
構成し、ポンピング波による測定用バースト波の
反射への非線形効果を観測し、被観察部位の等価
非線形パラメータ(B/A)eの変化量の分布を得る ことのできる超音波媒体の非線形パラメータ分布
測定装置を提供するにある。
An object of the present invention is to use an ultrasonic burst wave for measurement and a continuous ultrasonic wave for pumping, and to make the pumping wave close to a plane wave that covers the area to be observed and intersects the burst wave for measurement. A nonlinear parameter distribution of the ultrasonic medium that can be configured, observe the nonlinear effect on the reflection of the measurement burst wave by the pumping wave, and obtain the distribution of the amount of change in the equivalent nonlinear parameter (B/A) e of the observed region. To provide measuring equipment.

本発明は、超音波媒体内で交差する測定用超音
波バースト波とポンピング用超音波連続波を利用
し、その交差領域から反射される超音波バースト
波の特定の周波数成分を分離して測定することに
より、該交差領域の等価非線形パラメータの変化
量を検知し、さらにこの交差領域を移動させて測
定を繰り返すことにより、該非線形パラメータの
分布を二次元的あるいは三次元的に表示できるよ
うにしたものである。
The present invention utilizes a measuring ultrasonic burst wave and a pumping ultrasonic continuous wave that intersect in an ultrasonic medium, and separates and measures a specific frequency component of the ultrasonic burst wave reflected from the intersecting region. By this, the amount of change in the equivalent nonlinear parameter of the intersection area is detected, and by further moving this intersection area and repeating the measurement, it is possible to display the distribution of the nonlinear parameter two-dimensionally or three-dimensionally. It is something.

超音波媒体内の音圧がゼロの時の音速をC0
密度をρ0とすると測定波と直交する方向からPな
る音圧が加えられたときの音速Cは C=C0+1/2ρ0C0(B/A)eP……(1) ここで(B/A)eは次式で定義され、「等価非線形 パラメータ」と呼ばれるものであり、不均一媒質
においては場所により異つた値を取る。尚、カツ
コ右下のeは等価(equivalent)の意味である。
The speed of sound when the sound pressure in the ultrasonic medium is zero is C 0 ,
If the density is ρ 0 , the sound speed C when a sound pressure P is applied from the direction perpendicular to the measurement wave is C=C 0 +1/2ρ 0 C 0 (B/A) e P...(1) Here (B/A) e is defined by the following equation and is called an "equivalent nonlinear parameter" and takes different values depending on the location in a non-uniform medium. Note that the letter e at the bottom right of the katsuko means equivalent.

(B/A)e=2ρ0C0(∂Ci/∂Pj)s ……(2) iとjは直交する方向であり、Sは等エントロ
ピーであることを示す。
(B/A) e = 2ρ 0 C 0 (∂Ci/∂Pj) s ...(2) i and j are orthogonal directions, and S indicates isentropy.

一方、一般に音響の分野で用いられている非線
形パラメータB/Aは次のように定義されている。
On the other hand, the nonlinear parameter B/A, which is generally used in the field of acoustics, is defined as follows.

B/A=2ρ0C0(∂Ci/∂Ci)s ……(3) 等方体であればポアソン比をνとして(B/A)eと B/Aとの間には (B/A)e=(1+2ν−1/1−ν B/A)……(
4) という関係があり、液体、生体組織ではν≒0.5
であるので(B/A)eB/Aとみなすことができる
B/A=2ρ 0 C 0 (∂Ci/∂Ci) s ...(3) If it is an isotropic body, Poisson's ratio is ν and (B/A) between e and B/A is (B/ A) e = (1+2ν-1/1-ν B/A)...(
4) There is a relationship, ν≒0.5 in liquids and biological tissues.
Therefore, it can be regarded as (B/A) e B/A.

従つて、ポンピング波の音圧Pにより、音速C
は、 △C=1/2ρ0C0(B/A)・P ……(5) だけ変化する事になる。
Therefore, due to the sound pressure P of the pumping wave, the sound speed C
will change by △C=1/2ρ 0 C 0 (B/A)・P...(5).

第1図に反射法による非線形パラメータの測定
原理を示す。Xs,Xpはそれぞれ測定用バースト
波およびポンピング用連続波を放射する探触子、
Mは超音波媒体を示す。一般に生体軟組織のよう
な超音波媒体Mの密度・音速の不連続面での超音
波の反射係数r0は次式で与えられる。
Figure 1 shows the principle of measuring nonlinear parameters using the reflection method. Xs and Xp are probes that emit measurement burst waves and pumping continuous waves, respectively;
M indicates an ultrasonic medium. Generally, the reflection coefficient r 0 of an ultrasound wave at a discontinuous surface of density and sound velocity of an ultrasound medium M such as a biological soft tissue is given by the following equation.

r0=ρ1C1−ρ2C2/ρ1C1+ρ2C2 ……(6) 但しρi、Ci(i=1、2)はそれぞれ領域iで
の組織の密度と音速である。
r 0 = ρ 1 C 1ρ 2 C 2 / ρ 1 C 1 + ρ 2 C 2 ...(6) However, ρi and Ci (i=1, 2) are the tissue density and sound speed in region i, respectively. .

今、この不連続面のポンピング用超音波連続波
(ポンプ波)が照射され、組織内部に次式のよう
な音圧変化が与えられたとする。
Now, suppose that this discontinuous surface is irradiated with a pumping ultrasonic continuous wave (pump wave), and a sound pressure change as shown in the following equation is given inside the tissue.

△P=△P0ejpt ……(7) 但しωp、△P0はそれぞれポンプ波の周波数、
音圧振幅である。この音圧変化によつて領域1,
2の音速は変化し、各々の領域での非線形パラメ
ータ(B/A)1,(B/A)2を使つて(5)式から次
のように表わされる。
△P=△P 0 e jpt ……(7) However, ωp and △P 0 are the frequency of the pump wave, respectively,
It is the sound pressure amplitude. Due to this sound pressure change, area 1,
The sound speed of 2 changes and can be expressed as follows from equation (5) using the nonlinear parameters (B/A) 1 and (B/A) 2 in each region.

C1=C10+(B/A)1/2ρ1C10△P C2=C20+(B/A)2/2ρ2C20△P ……(8) 但しC10,C20はポンプ波による音圧変化が生ず
る前の音速である。
C 1 = C 10 + (B/A) 1 /2ρ 1 C 10 △P C 2 = C 20 + (B/A) 2 /2ρ 2 C 20 △P ...(8) However, C 10 and C 20 are This is the speed of sound before the sound pressure changes due to pump waves.

ここで、測定用バースト波の音圧は、ポンプ波
の音圧に比べて充分小さく超音波媒体内にもたら
される非線形パラメータの影響はポンプ波による
ものだけとした。
Here, the sound pressure of the measurement burst wave is sufficiently small compared to the sound pressure of the pump wave, and the influence of nonlinear parameters brought into the ultrasonic medium is only caused by the pump wave.

そこで、(6)、(8)式を使い、音圧変化△Pが加え
られた時の反射係数r(△P)を求めると、 r(△P)r10+ρ1k1−ρ2k2/ρ1C10+ρ2C20
P……(9) となる。但しk1=(B/A)1/2ρ1C10、 k2=(B/A)2/2ρ2C20、r10=ρ1C10−ρ2C20
ρ1C10+ρ2C20 でかつk1,k2<<1とした。
Therefore, using equations (6) and (8) to find the reflection coefficient r (△P) when the sound pressure change △P is added, we get r (△P) r 10 + ρ 1 k 1 − ρ 2 k 21 C 102 C 20
P...(9) becomes. However, k 1 = (B/A) 1 /2ρ 1 C 10 , k 2 = (B/A) 2 /2ρ 2 C 20 , r 10 = ρ 1 C 10 −ρ 2 C 20 /
ρ 1 C 102 C 20 and k 1 , k 2 <<1.

さらに生体軟組織においては、密度・音速の変
化は小さいから、組織内部の平均的密度・音速を
ρ00とすると、(9)式は r(△P) r0+(B/A)1−(B/A)2/4ρ0 C0 2△P……(1
0) となる。
Furthermore, in biological soft tissue, changes in density and sound speed are small, so if the average density and sound speed inside the tissue are ρ 0 and 0 , equation (9) is r(△P) r 0 + (B/A) 1 −(B/A) 2 /4ρ 0 C 0 2 △P……(1
0).

よつて今、有限時間幅で継続する周波数ωsの
超音波バースト波Aを組織に入射し、この不連続
面からの反射バースト波を測定すると、その受波
信号A′は A′∝r(△P)・A ……(11) で与えられる。但しA=A0ejstでA0は音圧振幅
であるとすると(11)式は A′∝A0{r0ejst+(B/A)1−(B/
A)2/4ρ0 C0 2△P0ej(s+p)t}……(12) となる。
Therefore, if we inject an ultrasonic burst wave A with a frequency ωs that continues over a finite time width into the tissue and measure the reflected burst wave from this discontinuous surface, the received signal A' will be A'∝r(△ P)・A ……(11) is given by. However, if A = A 0 e jst and A 0 is the sound pressure amplitude, then equation (11) becomes A'∝A 0 {r 0 e jst + (B/A) 1 − (B/
A) 2 /4ρ 0 C 0 2 △P 0 e j(s+p)t }...(12).

この受波信号A′において周波数ωsの成分はr0
周波数(ωs+ωp)の成分は
(B/A)1−(B/A)2/4ρ0 C0 2△P0にそれぞれ比
例する。
In this received signal A′, the component of frequency ωs is r 0 ,
The frequency (ωs+ωp) component is proportional to (B/A) 1 −(B/A) 2 /4ρ 0 C 0 2 ΔP 0 .

即ちこの反射信号A′のスペクトルを図示すると
第2図のようになる。
That is, the spectrum of this reflected signal A' is shown in FIG. 2.

即ち、受波信号A′に対し、周波数ωs付近での
バンドパスフイルタ出力をとれば従来のパルスエ
コー法において得られる情報と対応するものが得
られるし、周波数ωs+ωp付近でのバンドパスフ
イルタ出力をとれば、非線形パラメータの変化量
に対応するものが得られる。この測定を一平面上
で行なうと、従来のBモード断層像に対応する画
像と、非線形パラメータの分布に対応する画像が
得られる。
That is, if we take the bandpass filter output around the frequency ωs for the received signal A', we can obtain information that corresponds to the information obtained in the conventional pulse echo method, and we can obtain the bandpass filter output around the frequency ωs + ωp. If we take the value, we can obtain the value corresponding to the amount of change in the nonlinear parameter. When this measurement is performed on one plane, an image corresponding to a conventional B-mode tomographic image and an image corresponding to the distribution of nonlinear parameters are obtained.

次に、装置全体の構成を第3図に、信号処理の
タイミングチヤートを第4図に示す。システムコ
ントローラ1はトリガー発生器2に指令を出しト
トリガーパルスを発生させる。トリガパルスに応
答してバースト波アンプ3からバースト波が出力
され、送信プローブ4に与えられる。また遅延回
路6に遅延幅を指定させ、受波プローブ5で受け
た信号への時間ゲート回路7の処理を制御する。
この時間ゲートは信号波バースト波とポンプ波が
交差する領域からの反射バースト波だけが通過す
るように設定される。さらにシステムコントロー
ラはポンプ波・信号波の交差領域を指定し、移動
するよう、ビームコントローラ8に指令を送る。
時間ゲートを通つた信号は周波数ωs付近または
ωs+ωpで透過するバンドパスフイルタ9に通
し、それぞれ組織の反射係数、非線形パラメータ
の変化に対応する信号に回路10,11にて分離
される。この信号を画像メモリ12に蓄え、さら
に一画面分の信号を蓄えた後、システムコントロ
ーラの指示により、表示制御部13の制御のもと
にBモード像または非線形パラメータ分布像をデ
イスプレイ14,15にて表示するようになつて
いる。
Next, the configuration of the entire device is shown in FIG. 3, and the timing chart of signal processing is shown in FIG. 4. The system controller 1 issues a command to the trigger generator 2 to generate a trigger pulse. In response to the trigger pulse, a burst wave is output from the burst wave amplifier 3 and applied to the transmission probe 4. It also causes the delay circuit 6 to specify a delay width and controls the processing of the time gate circuit 7 on the signal received by the receiving probe 5.
This time gate is set so that only the reflected burst wave from the region where the signal wave burst wave intersects with the pump wave passes through. Furthermore, the system controller specifies the intersection area of the pump wave and signal wave and sends a command to the beam controller 8 to move.
The signal that has passed through the time gate is passed through a bandpass filter 9 that transmits the signal at a frequency near ωs or ωs+ωp, and is separated by circuits 10 and 11 into signals corresponding to changes in the tissue reflection coefficient and nonlinear parameters, respectively. After storing this signal in the image memory 12 and further storing signals for one screen, a B-mode image or a nonlinear parameter distribution image is displayed on the displays 14 and 15 under the control of the display control section 13 according to instructions from the system controller. It is now displayed as follows.

尚、16はポンプ波用プローブである。 Note that 16 is a pump wave probe.

測定点の走査は、時間ゲートの送延量を変える
ことにより、測定用ビームに沿つた一次元的走査
が可能であり、またビーム方向と直角方向への走
査は送・受波プローブ4,5をリニアアレイ型探
触子とするか、又はセクタスキヤン型の探触子と
すれば容易に可能である。或いは測定用ビームと
ポンプ波とをともに細く絞り、かつ両探触子の位
置関係を固定としたうえで、両探触子の組全体を
2次元的に移動させることにより、両ビームの交
点を2次元的に走査することも可能である。
The measurement point can be scanned one-dimensionally along the measurement beam by changing the amount of time gate spread, and scanning in the direction perpendicular to the beam direction can be done using the transmitting/receiving probes 4 and 5. This can easily be achieved by using a linear array type probe or a sector scan type probe. Alternatively, by narrowing down both the measurement beam and the pump wave, fixing the positional relationship between both probes, and then moving the entire pair of probes two-dimensionally, the intersection point of both beams can be adjusted. It is also possible to scan two-dimensionally.

測定用超音波バースト波は有限の継続時間をも
つため測定には空間的分解能の劣化がもたらされ
る。これを避けるためには第5図に示す様にでき
るだけ細く絞られたポンプ波を用い、これとバー
スト波が直交するよう配置し、交差領域を小さく
とればよい。このことにより、空間的分解能の劣
化が改善される。但し、図中の記号は第1図のも
のと同様である。もちろん、ポンプ波を細く絞ぼ
るということは、ポンプ波が平面波であるという
仮定からはずれ、ポンプ波ビーム軸上で、距離と
ともに音圧分布が変化する。そこで前もつてポン
プ波の音圧分布を測定し、測定用バースト波とポ
ンプ波との交差領域の位置によつて、このポンプ
波の音圧を補正しておけば問題はない。
Since the measurement ultrasonic burst wave has a finite duration, the measurement results in a degraded spatial resolution. In order to avoid this, as shown in FIG. 5, a pump wave narrowed as narrow as possible is used, and the burst wave is arranged so as to be orthogonal to the pump wave, so that the crossing area is made small. This improves the deterioration of spatial resolution. However, the symbols in the figure are the same as those in FIG. Of course, narrowing down the pump wave deviates from the assumption that the pump wave is a plane wave, and the sound pressure distribution changes with distance on the pump wave beam axis. Therefore, there is no problem if the sound pressure distribution of the pump wave is measured in advance and the sound pressure of the pump wave is corrected based on the position of the intersection area of the measurement burst wave and the pump wave.

本発明によれば、生体組織など超音波媒体の非
線形パラメータの分布像が表示できるので、今ま
でのBモード像では見えない特徴が見いだせる。
また反射型映像系のため扱いやすく、かつ従来の
Bモード像も表示できる。
According to the present invention, a distribution image of nonlinear parameters of an ultrasound medium such as a living tissue can be displayed, so that features that cannot be seen in conventional B-mode images can be found.
Furthermore, since it is a reflection type image system, it is easy to handle and can also display conventional B-mode images.

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

第1図は反射型非線形パラメータ測定原理を説
明する図、第2図は受波信号A′のスペクトルを
示す図、第3図は一実施例装置の全構成を示す
図、第4図は信号処理のタイミングチヤートであ
る。また第5図はポンピング用超音波連続波を細
く絞つた場合の反射型非線形パラメータ測定を示
す図である。 第3図において、4,5は送・受波用プロー
プ、16はポンプ波用プローブ、9はバンドパス
フイタ、10は反射係数を求める回路、11は非
線形パラメータを求める回路である。
Fig. 1 is a diagram explaining the principle of reflection type nonlinear parameter measurement, Fig. 2 is a diagram showing the spectrum of the received signal A', Fig. 3 is a diagram showing the entire configuration of an embodiment of the device, and Fig. 4 is a diagram showing the signal This is a timing chart of processing. FIG. 5 is a diagram showing reflection type nonlinear parameter measurement when the pumping ultrasonic continuous wave is narrowed down. In FIG. 3, 4 and 5 are transmitting and receiving probes, 16 is a pump wave probe, 9 is a bandpass filter, 10 is a circuit for determining a reflection coefficient, and 11 is a circuit for determining nonlinear parameters.

Claims (1)

【特許請求の範囲】 1 (イ) 超音波媒体内に測定用の超音波バースト
波を送・受波する1個の兼用の、または個別の
超音波探触子と、 (ロ) ポンピング用の超音波連続波の平面波を送波
する超音波探触子と、 (ハ) 該測定用超音波バースト波と該ポンピング用
超音波連続波が該超音波媒体内で交差する配置
と、 (ニ) 該配置において、該測定用超音波探触子で受
波した該交差領域上での特定の場所からの反射
超音波バースト波の特定の周波数成分だけを取
り出し、該交差領域の特定の場所での非線形パ
ラメータの変化量を求める信号処理回路と、 (ホ) 該交差領域又は該交差領域上の該特定場所を
ずらし、上記非線形パラメータの変化量の分布
を求める手段とを有することを特徴とする超音
波媒体の反射型非線形パラメータ分布測定装
置。
[Claims] 1. (a) A dual-purpose or individual ultrasonic probe that transmits and receives ultrasonic burst waves for measurement into an ultrasonic medium; (b) A pumping probe. an ultrasonic probe that transmits a plane wave of a continuous ultrasonic wave; (c) an arrangement in which the measurement ultrasonic burst wave and the pumping ultrasonic continuous wave intersect within the ultrasonic medium; In this arrangement, only a specific frequency component of the reflected ultrasonic burst wave received by the measuring ultrasonic probe from a specific location on the intersection area is extracted, and (e) means for shifting the intersection region or the specific location on the intersection region to obtain a distribution of the variation in the nonlinear parameter; Reflection type nonlinear parameter distribution measuring device for acoustic wave media.
JP16859983A 1983-09-13 1983-09-13 Reflective type non-linear parameter distribution apparatus of ultrasonic medium Granted JPS6060837A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP16859983A JPS6060837A (en) 1983-09-13 1983-09-13 Reflective type non-linear parameter distribution apparatus of ultrasonic medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP16859983A JPS6060837A (en) 1983-09-13 1983-09-13 Reflective type non-linear parameter distribution apparatus of ultrasonic medium

Publications (2)

Publication Number Publication Date
JPS6060837A JPS6060837A (en) 1985-04-08
JPH0421493B2 true JPH0421493B2 (en) 1992-04-10

Family

ID=15871037

Family Applications (1)

Application Number Title Priority Date Filing Date
JP16859983A Granted JPS6060837A (en) 1983-09-13 1983-09-13 Reflective type non-linear parameter distribution apparatus of ultrasonic medium

Country Status (1)

Country Link
JP (1) JPS6060837A (en)

Also Published As

Publication number Publication date
JPS6060837A (en) 1985-04-08

Similar Documents

Publication Publication Date Title
US9439620B2 (en) Ultrasound diagnosis apparatus
US5361767A (en) Tissue characterization method and apparatus
US7785259B2 (en) Detection of motion in vibro-acoustography
US20140046173A1 (en) Ultrasound vibrometry with unfocused ultrasound
EP0139235B1 (en) Ultrasonic measurement method and apparatus therefor
US4511984A (en) Ultrasound diagnostic apparatus
JPH0570769B2 (en)
JP2010526626A (en) Method and apparatus for measuring the average value of viscoelasticity of a region of interest
JPH0467149B2 (en)
Kelsey et al. Applications of ultrasound in speech research
JP2012239546A (en) Ultrasound diagnostic apparatus
JP2003230560A (en) Ultrasonograph
CN110494082B (en) Ultrasound elastography method and system
JPH0421493B2 (en)
JP3280473B2 (en) Ultrasound diagnostic equipment
JP4077913B2 (en) Ultrasonic diagnostic equipment
JPH0548130B2 (en)
JPS61106137A (en) Ultrasonic diagnostic apparatus
Mobley et al. Evaluation of ultrasound techniques for brain injury detection
JPS6268441A (en) Ultrasonic tissue diagnostic apparatus
JPS6272336A (en) Ultrasonic tissue diagnostic apparatus
JPH0390139A (en) Ultrasonic measuring device with echo of blood flow as reference
Le Croissette Signal Processing In Ultrasound
JPH0512934B2 (en)
JPH0339146A (en) Ultrasonic measuring apparatus