JPH02119851A - 超音波検査装置 - Google Patents

超音波検査装置

Info

Publication number
JPH02119851A
JPH02119851A JP63274201A JP27420188A JPH02119851A JP H02119851 A JPH02119851 A JP H02119851A JP 63274201 A JP63274201 A JP 63274201A JP 27420188 A JP27420188 A JP 27420188A JP H02119851 A JPH02119851 A JP H02119851A
Authority
JP
Japan
Prior art keywords
vibration
frequency
amplitude
spectrum
doppler effect
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
JP63274201A
Other languages
English (en)
Other versions
JP2699473B2 (ja
Inventor
Hirokuni Sato
佐藤 拓宋
Yoshiki Yamakoshi
芳樹 山越
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.)
Shimadzu Corp
Original Assignee
Shimadzu 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 Shimadzu Corp filed Critical Shimadzu Corp
Priority to JP63274201A priority Critical patent/JP2699473B2/ja
Publication of JPH02119851A publication Critical patent/JPH02119851A/ja
Application granted granted Critical
Publication of JP2699473B2 publication Critical patent/JP2699473B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

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

Description

【発明の詳細な説明】 【産業上の利用分野】
この発明は、超音波によって物体内部を検査する装置に
関する。
【従来の技術】
従来より、超音波検査装置は、とくに医学の分野におい
て診断の重要な手段として定着してきている。現在医用
診断の手段として用いられている超音波診断装置は、大
多数が超音波の反射を利用して生体内部の組織の不連続
性を映像化しようとするものに過ぎず、軟部組織の病変
(癌や肝硬変など)に対しては識別能力がない。 そこで、最近、生体内部の弾性的な特性に間する情報を
得るため、生体外部から周波数10HzないしI K 
Hz程度の低周波振動を加え、生体内部に伝わる振動変
位の分布を、同時に送波したプローブ用超音波により計
測することが提案されている(山越芳樹、鈴木正人、佐
藤拓宋:低周波振動とプローブ用超音波を用いた弾性的
特性の映像系、日超医論文集 50 301 1987
)。 これは、同時に送波した超音波の反射波の振幅の変化か
ら内部の振動振幅の大きさを間接的に推定し、その分布
より内部の弾性特性に関する情報を得ようというもので
ある。 また、生体内部の振動振幅を超音波を用いて計測するた
め、任意の2つの時点で得られた超音波画像などの相互
相関関数を求め、そのピーク位置より移動量を推定しよ
うという提案もなされている(稲葉克、へ木晋−1中山
淑:RFエコー相関方による変位測定、日超医論文集 
50 303一 1987)  。
【発明が解決しようとする課題】
しかしながら、前者にあっては、内部振動振幅の直接計
測ができず、さらに反射振幅の統計的な変動を取り除く
ために、得られた映像に移動平均処理を行なうことが必
要で、映像の分解能が劣るという問題点がある。 また、後者では、加振により内部の散乱体が超音波ビー
ムに対して回転してしまうと相互相関関数に明確なピー
クが現われなくなり、このときに正確な変位量の測定が
難しいという問題点がある。 この発明は、被検査体内部の振動の振幅や位相の正確な
測定ができる、超音波検査装置を提供することを目的と
する。
【課題を解決するための手段】
上記目的を達成するため、この発明による超音波検査装
置においては、被検査体に対し低周波の機械的振動を加
える手段と、該被検体に対して超音波を送ってその反射
波を受ける超音波送受手段と、該受波信号からドツプラ
ー効果による周波数変調成分を取り出す手段とが備えら
れる。
【作  用】
超音波の反射波は被検査体内部の振動によってドツプラ
ー効果を受け、その周波数がシフトする。 そこで、受波信号からドツプラー効果による周波数変調
成分を取り出せば、そのスペクトルから変調指数を求め
て内部の振動振幅を直接計測したり、振動の位相を計測
することが可能となる。 得られたデータの統計的な変動を押さえるための移動平
均処理などが不必要で、空間的に高分解能なデータを計
測できるため、高分解能な画像が得られる。
【実 施 例】
つぎにこの発明の一実施例について図面を参照しながら
説明する。第1図の実施例において、メカニカルバイブ
レータ12は低周波発振器14の出力を増幅する増幅器
15により駆動され、その先端に取り付けられた適当な
大きさの振動板(たとえばアクリル板でなる)13より
被検査体10の表面に−様な振動変位を与えるようにし
ている。 メ カニカルバイブレータ つるされており、被検査体10の表面への振動板13の
接触圧力が嚢えちれるようにされている。 超音波トランスデユーサアレイ21は、振動板13から
の低周波振動が伝わらないように、柔らかいコンタクト
ゲル22を間に挟んで被検査体10の表面に接触させら
れている。このトランスデユーサアレイ21にはパルス
発生器23からのこの実施例では3.5MHzの信号が
増幅器24を介して送られてきており、被検査体10の
内部に超音波を送波する。被検査体10からの反射波が
このトランスデユーサアレイ21によって受波され、受
波信号は増幅器26を経て、90°移相器31、乗算器
32、33、ローパスフィルタ34、35、サンプルホ
ールド回路36、37よりなる周波数変調検出部と、エ
ンベロープ検波器41、ローパスフィルタ42及σサン
プルホールド回路43よりなる振幅変調検出部とに送ら
れる。乗算器32、33はパルス発生器23からの周波
数3。 5MHzの参照信号と受波信号とを乗算するもので、ロ
ーパスフィルタ34、35は高域カット周波数が約50
0KHzのものでこれにより高調波成分が除去される。 このローパスフィルタ34、35の出力はサンプルホー
ルド回路36、う7によりサンプリングされ、あらかじ
め指定した測定点からの信号だけが時系列信号として取
り出される。 受波信号をそのまま及び90°移相して2系統で処理を
行なっているのは直交検波方式とするためであり、サン
プリングされて得た信号は一方が実部、他方が虚部の信
号としてA/Dコンバータ51を経てマイクロコンピュ
ータ52に取り込まれる。被検査体10を振動させる場
合、その内部の反射体の運動によって、反射してくる超
音波はドツプラー周波数シフトのほかに振幅変調も受け
ているので、振幅変調検出信号をA/Dコンバータ51
を経てマイクロコンピュータ52に入力して、その影響
を補正するようにしている。こうして各測定点ごとに低
周波振動の振幅が求められ、その2次元分布像がデイス
プレー装置53でカラ−表示される。なお、タイミング
信号発生器54は各部の動作タイミングを決めるための
信号を発生する。この実施例では、超音波送波パルスの
繰り返し周期は最高200μsであり、振動振幅測定の
ために加振周波数の第5次高調波成分まで用いるとする
と、約500 Hzまでの低周波振動を加えることがで
きる。 つぎに、上記の実施例の動作原理を説明する。 まず、振動板13により被検査体10に対し10Hz〜
IKHz程度の低周波振動を加えると、その内部の反射
体がその周波数で振動するので、反射体の振動速度は0
式のように表わすことができる。 V=Vosin ωbt            ■こ
こで、ωbは加えた振動の角周波数、Voは速度振幅で
ある。 こうして振動している反射体に向けて超音波を送波する
と、この反射体で反射し、その反射超音波はドツプラー
効果により周波数変調を受ける。 このとき周波数シフト量Δfは0式で与えられる。 Δf、=f  (2V/C) 一Δf wax Sin O1) t        
    ■ここでfは超音波の周波数、Cは音速であり
、Δf y+nxはドツプラー効果による最大周波数シ
フトで、つぎのようになる。 Δf 、、!= 2 f (V o/ C)     
     ■一方、超音波トランスデユーサアレイ21
で受波した反射波の受波信号r(t)は、 r(t)−r。5in(ωt+m(sinωbt+φ)
     ■と表わすことができる。ここで、roは反
射超音波の振幅であり、ωはドツプラー効果を受けてい
ないときの超音波の角周波数、φは伝播による位相差で
ある6また、meは周波数変調の変調指数であり、 m r= 2 πΔf 1111.X / Ot1  
         ■と表わせる。 いま、反射体の振動振幅をξ。とすると、この変調指数
は、 mr−(4πf/C)ξ0         ■と表わ
すこともできる。 さて、受波信号r(t)からドツプラー効果による周波
数変調成分のみを取り出すために、参照信号5(t)を
、 5(t)・50(cosωt+jsinωt〉    
    ■とじて、r(t)と5(t)とを乗算しく乗
算器32.33)、ローパスフィルタ34.35で周波
数2ω付近のスペクトラムを持つ信号成分を除去すれば
、ローバルフィルタ34.35の出力信号0Lpp(t
)は、 OLPF (t) = + r o s o ×((s
inφ+jcosφ)[Jo(mr)+2ΣJ2 n 
(m r )cos2nωbt]+(cosφ−jsi
nφ)[2ΣJzn+t(mr)sin(2n+1)ω
btll■ どなる。ここで、Jn(x)はn次のベッセル関数であ
る。0式をみると、0LPF(t)は直流分と角周波数
ωb、2ωb、3ωb、・・・の輝線スペクトラムを持
つことが分かるが、そのi次のスペクトラムの大きさA
、はつぎの0式で与えられる。 Ao=+ r O3(1+J、O(mr)A i=r 
as o J o (mr>  (14,2,・・、’
)   ■つまり各スペクトラムの大きさは変調指数m
fの関数になるため、振動振幅が変化するとそれに応じ
て各スペクトラムの大きさが変化するようになる。この
ことより、Ot、pp(t)のスペクトラムの大きさか
ら変調指数を推定すれば、その点での振動振幅が推定で
きることになる。ここでは、隣り合う次数のスペクトラ
ムの大きさの比から振動振幅を推定するようにした。隣
り合う次数のスペクトラムの比は、 Al/Ao= 2 J t (mr) / Jo (m
r)At+t/ At−2J tit (mr) / 
J i(mr)  (IQ)となる。各スペクトラムの
大きさと、それらの比の値を変調指数をパラメータとし
て示すと第2図(a)、(b)のようになる。このよう
にして振動変位を推定するには、まず0LPF(t)の
スペクトラムを計算し、隣り合うスペクトラムの大きさ
の比をとり第2図(b)を参照して変調指数m。 を求め、0式から導かれる振動振幅と変調指数との関係 ξ。−(C/4πf)m、         ■を使っ
て振動振幅を求めればよい。しがし、第2図(b)を見
ても分がるように一般にスペクトラム比の大きさに対し
てm(の値は多価となっており、このままでは振動変位
の一意な推定ができない。そこで、ここではスペクトラ
ムの大きさを比較し振動振幅が入る大まかなレンジを求
め、それをもとに推定のために用いるスペクトラム比を
適応的に変えることで内部振動変位を一意に推定できる
ようにしている。具体的には、まずoLPF(t)の直
流成分、1次、2次、3次スペクトラム成分のなかで、
もし直流成分が最大の値をとれば、振動振幅は小さいこ
とが予想されるからスペクトラムの直流成分と1次の成
分との比を用いて振動振幅を求める。直流成分が最大で
ない場合には、1次成分と2次成分との比や、2次成分
と3次成分との比を用いて求める。このような方法を高
次のスペクトラムにまで適応していくことにより振動振
幅を一意に推定することができる。 この第1図の実施例での実験結果について説明すると、
まず、観測されたドツプラー検出信号〇LPF(t )
は第3図、のようになった。この第3図(a)、(c)
で01.t、pp(t ) 、02.t、pp(t )
はそれぞれ0LPF(t)の実部と虚部を表わす。ここ
では、孔検査体10として平均粒径45μmのグラファ
イト粉末を重量比で1,5%混ぜたグラファイト混入カ
ンテンファントム(カンテン粉末の重量比は1.5%)
を用い、表面から40mm下の位置から反射してくる超
音波を捉えるようにした。第3図(a)は加振した低周
波振動の周波数を1゜OHzとし、メカニカルバイブレ
ータ12への印加電圧を2Vとしたときの01.LPF
(t ) 、02.LPF(t )を表わし、第3図(
b)は第3図(a)の信号をフーリエ変換して求めたス
ペクトラムである。また、第3図(c)は加振周波数は
100H2のまま、メカニカルバイブレーク12への印
加電圧を6■に上げたときの○1.LPF(t ) 、
02.LPF(t )を表わし、第3図(d)はその信
号をフーリエ変換して求めたスペクトラムである。この
第3図(b)、(d)から、さきに原理を示したように
、加振周波数とその高調波成分からなる輝線スペクトラ
ムになっていることが分かる。これら第3図(a)、(
b)と第3図(c)、(d)とを比較すると、メカニカ
ルバイブレータ12への印加電圧が高くなって振動振幅
を大きくすれば、ドツプラー検出信号が歪み、高調波成
分が大きく現われることがよく分かる。 このスペクトラムの大きさの比から振動振幅が求められ
るが、こうして求められた振動振幅を、直径約2.5m
mの領域内での振動振幅をサブミクロンの精度で非接触
に計測できる光変位計を用いて測定した結果と比較して
みたところ、非常に高い精度となっていることが確認で
きた。まず、第4図(a)、(b)に示すように、水槽
61の底部に直径約20+nmの円形の孔をあけ、そこ
にビニールフィルム62を張り水63を入れたものを用
意した。第4図(a)では、縦10mmX横10+nm
×高さ5mmのカンテンのみのファントム64をこのビ
ニールフィルム62上に接着し、ファントム64の真下
からメカニカルバイブレーク12及び振動板13で振動
を与え、水槽61の水面側に超音波トランスデユーサア
レイ21を配置してこのファントム64の表面における
振動振幅を測定した。第4図(b)では、このビニール
フィルム62の上に反射体としてアルミニウム箔65を
張りつけ、光変位計のプローブ66を上方に配置して振
動振幅を測定した。ここでは加振周波数を100Hzと
しており、メカニカルバイブレータ13に加える電圧を
変化させて振幅を変えて測定したところ、第5図のよう
な結果が得られた。この測定結果は、本実施例の場合も
光変位計を用いた場合も、表面上の5点で測定してその
平均値を求めたものである。この第5図から、この実施
例による測定値は、光変位計で得られた測定値とよく一
致しており、振動振幅の定量的な測定が精度高くできる
ことが分かる。 ここまでの計測では、スペクトラム解析のために約10
0 m s分の時系列データを用いて振動変位を推定し
たが、映像を実時間的に得ようとすると、さらに短い時
系列データからスペクトラム解析を行なう必要がある。 一般にここで用いているような高調波成分の解析には最
低、振動の1周期分のデータがあればよいが、データ長
が短いとスペクトラム分解能が低下してしまう。そこで
、さきに示したグラファイト混入カンテンファントム6
3(第4図〉を用いて、振動変位の違う2点についてス
ペクトラム解析のために用いるデータ長を変え、推定精
度への影響を実験的に評価した。 その結果は、第6図に示される。この実験では100H
zで加振したが、最短の1周期分のデータから推定した
振動変位と5周期分のデータから推定した振動変位とは
、2点とも、はぼ同じ値となっており、この場合は、1
周期分のデータでも正しい測定ができているといえる。 このため、以降の実験では映像を高速に得ることができ
る1周期分のデータを用いて振動変位の測定を行なうこ
ととしている。 つぎに、加振周波数について実験的に検討する。 同一のファントムに対して加振周波数を変えて内部の振
動振幅の分布を測定する。ファントムとして、比較的軟
らかいグラファイト混入カンテンファントム(カンテン
粉末の重量比;1.5%、グラファイト粉末の重量比;
1.5%)の中に縦50 mm X横50鵬×高さ50
+nmの比較的硬いファントム(カンテン粉末の重量比
;4%、グラファイト粉末の重量比;1.5%)を埋め
込んだものを使用した。第7図(a)、(b)、(c)
、(d)にこのファントムでのY方向(奥行き方向)の
振動変位の1次元分布を示す。図中ファントム上端から
40mmより90mmの点線で囲んだ部分に硬いカンテ
ンが埋め込まれているが、加振周波数が20Hz (第
7図(a→)と60Hz (第7図(b))とでは2つ
のカンテンの違いがこの図からはよく区別できていない
。これは周波数が低いとファントム全体が振動してしま
うためと考えられる。 また加振周波数が140Hz (第7図(d))では硬
いカンテン上面での弾性特性の違いがよく映像化されて
いるものの、周波数が高いため奥行き方向への浸透が悪
く硬いカンテンで振動が遮られてしまっている。この中
間の周波数100H2(第7図(C))では、2つのカ
ンテンの間の弾性的な特性の違いがよく映像化されてお
り、また硬いカンテンの下部に振動が入り込まない部分
が生じているものの、表面から約120mmを超える部
分で再び振動振幅が大きくなり、このようなファントム
では100Hz前後の加振周波数を用いればよいことが
分かる。 最後に、上記実施例で得られた内部の振動変位の映像例
を示す。実際には、ここで示す2次元画像はカラー画像
として表示されるが、ここではモノクロ画像として図示
している。第8図はカンテン粉末の重量比が1.5%で
グラファイト粉末の重量比が1,5%の−様なグラファ
イト混入カンテンファントムに対して得た映像である。 ここで、加振周波数は100Hzで、振動板13として
30mmX、20+nmの大きさのアクリル板を用いて
いる。 この図を見ると、ファントム内部に伝わる振動振幅の様
子がよく分かる。一方、第9図は、第7図で示した実験
で用いたのと同じ、の内部に縦50+nm×横50mm
X高さ50mmの比較的硬いファントムを埋め込んだグ
ラファイト混入カンテンファントムでの映像である。加
振周波数は100Hzとし、図中点線で囲まれた部分に
硬いカンテンが埋め込まれている。この像は、振動を硬
いカンテンの真上から加えたときの内部の振動変位分布
を示すもので、振動の伝播の様子がよく分かるが、硬い
カンテンの上面と右側面に比べて左側面と下面での違い
がよく映像かできていない。そこで、振動板の位置を左
に少し動かして映像をとり、両映像の各画素ごとに、振
動振幅の最大値をその画素の値として2つの映像を合成
して、第10図に示すような映像を得た。この映像では
奥にいくほど振動振幅を上げるような減衰補正を行なっ
ており、また50μm以上の振動振幅がある点はすべて
同色で示している。このような処理により、内部の振動
振幅から、ファントム外側の比較的軟らかい部分と、内
側の比較的硬い部分の弾性特性の違いをより明瞭に区別
して映像化できることが分かる。 第11図のように左上腕部にメカニカルバイブレータ1
2と振動板13と超音波トランスデユーサアレイ21と
を密着させ、IN  VIVOでの内部振動変位分布を
測定した。こうして、第12図(a)、(b)に示すよ
うな内部振動分布のY方向1次元像を得た。このとき、
加振周波数は100Hzで、第12図(a)は筋肉を弛
緩したとき、第12図(b)は緊張させたときである。 どちらも同じ部位を観測し、メカニカルバイブレータ1
2への印加電圧も同じ場合の結果であるが、筋肉を緊張
させると、弛緩時に比べて一般に組織が硬くなって振動
しにくくなる様子がよく分かる。 なお、本発明は上記の実施例に限定されるものではなく
、本発明の趣旨を逸脱しない範囲で種々に変更可能であ
る。たとえば、上記の実施例では映像化すべき2次元空
間の各点について1点ずつドツプラー効果による周波数
シフトを検出しているが、ドツプラー効果による周波数
シフトの検出を同時に複数点で行なうようにしてリアル
タイム性を実現することなども可能である。 また、振動振幅以外に、振動の伝搬速度に関連している
振動の位相分布を計測することも可能である。すなわち
、反射体の振動の位相φbを考慮すると、上記0式は、 ■−■。5in(ωbt÷φb)       ■′と
なる。このとき、反射波の受波信号はr(t)−r。5
in(ωt+mtsin(ωbt+φb)+φ) ■′
と表わすことができる。これを参照信号(0式)を用い
て直交検波すると、その複素検波信号(複素ドツプラー
検出信号)OLPF(t)の実部01LPF(1)、虚
部02LPF(t)はつぎのようになる。 01LPF (t) −+ r o s osinφ[
Jo (mr)+2ΣJ2n(Illr)cos(2n
ω、t+2nφb) ] ++rosocosφ[2Σ
’2n+1(u+r)sin((2n+1)ωbt+(
2n+1)φb)]      ■′02LPF (t
) −+  r o s  ocO3φ [Jo(mr
)+2ΣJzn(lllr)cos(2nωbt+2n
φb) ]  + r as oslnφ[2ΣJ2n
+1(lllr)sin((2n+1)ωbt+(2n
+1)φ5)1      ■″この■°式、■″式を
見ると、振動に位相差φbがあると、周波数nωbのス
ペクトラムに対してnφbの位相差を持つことが分かる
。つまり■′式、■゛式を各々フーリエ変換し、そのω
bのスペクトラムの偏角を導出すれば、φ、を求めるこ
とができる。そこで、先の振動振幅測定と組み合わせれ
ば、物体内部に伝わる振動の振幅と位相の分布を同時に
映像化できることになる。
【発明の効果】
この発明の超音波検査装置によれば、被検査体躯 内部の低周波振動の振幅を+μmの精度で測定したり、
振動の位相を測定できるので、生体軟部組織の硬さの違
いによる低周波振動加振に対する応答の差を検出でき、
癌の良性・悪性の区別などを容易に行える。
【図面の簡単な説明】
第1図はこの発明の一実施例のブロック図、第2図(a
)はドツプラーシフトの各スペクトラムの大きさと変調
指数との理論的関係を示すグラフ、第2図(b)はその
スペクトラムの比と変調指数との理論的関係を示すグラ
フ、第3図(a)、 (c)は実際のドツプラー検出信
号を示すグラフ、第3図(b)、(d)はそれらのスペ
クトラムをそれぞれ示すグラフ、第4図(a>、 (b
)は振動振幅の測定精度を確認するための実験例を示す
模式図、第5図は上記の実験結果として得られたデータ
を示すグラフ、第6図はデータ長と測定振動振幅との関
係を示すグラフ、第7図(a)、 (b)、 (c)、
 (d)は加振周波数を変えたときの測定振幅の奥行き
方向の1次元分布を表わすグラフ、第8図、第9図及び
第10図は映像例をそれぞれ示す図、第11図はIN 
 VIVO測定例を示す模式図、第12図(a)、 (
b)は第11図において測定した振幅のY方向1次元分
布を表わすグラフである。 10・・・被検査体、11・・・ばねばかり、12・・
・メカニカルバイブレータ、13・・・振動板、14・
・・低周波発振器、15.24.25・・・増幅器、2
1・・・超音波トランスデユーサアレイ、22・・・コ
ンタクトゲル、23・・・パルス発生器、31・・・9
0°移相器、32.33・・・乗算器、34.35.4
2・・・ローパスフィルタ、36.37.43・・・サ
ンプルホールド回路、41・・・エンベロープ検波器、
51・・・A/Dコンバータ、52・・・マイクロコン
ピュータ、53・・・デイスプレー装置、54・・・タ
イミング信号発生器、61・・・水槽、62・・・ビニ
ールフィルム、63・・・水、64・・・ファントム、
65・・・アルミニウム箔、 66・・・光変位計のプローブ、 71・・・上腕部。 (ト快 CO%

Claims (1)

    【特許請求の範囲】
  1. (1)被検査体に対し低周波の機械的振動を加える手段
    と、該被検査体に対して超音波を送ってその反射波を受
    ける超音波送受手段と、該受波信号からドップラー効果
    による周波数変調成分を取り出す手段とからなることを
    特徴とする超音波検査装置。
JP63274201A 1988-10-29 1988-10-29 超音波検査装置 Expired - Fee Related JP2699473B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP63274201A JP2699473B2 (ja) 1988-10-29 1988-10-29 超音波検査装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP63274201A JP2699473B2 (ja) 1988-10-29 1988-10-29 超音波検査装置

Publications (2)

Publication Number Publication Date
JPH02119851A true JPH02119851A (ja) 1990-05-07
JP2699473B2 JP2699473B2 (ja) 1998-01-19

Family

ID=17538442

Family Applications (1)

Application Number Title Priority Date Filing Date
JP63274201A Expired - Fee Related JP2699473B2 (ja) 1988-10-29 1988-10-29 超音波検査装置

Country Status (1)

Country Link
JP (1) JP2699473B2 (ja)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0556257A1 (en) * 1990-11-02 1993-08-25 The University Of Rochester Method and apparatus for breast imaging
JP2002243711A (ja) * 2001-02-14 2002-08-28 Aloka Co Ltd 振動計測方法及び超音波顕微鏡システム
JP2009531101A (ja) * 2006-03-29 2009-09-03 スーパー ソニック イマジン 粘弾性媒質を撮像する方法及び装置
CN103083045A (zh) * 2011-10-28 2013-05-08 Ge医疗系统环球技术有限公司 超声波诊断装置
JP2014030774A (ja) * 2013-11-22 2014-02-20 Super Sonic Imagine 関心領域の粘弾性の平均値を測定するための方法および装置
JP2014217605A (ja) * 2013-05-09 2014-11-20 国立大学法人群馬大学 映像装置及び映像形成方法
JP2015504755A (ja) * 2012-01-26 2015-02-16 パルティ、ヨーラム 肺振動中の経胸郭肺ドップラー超音波を使用した肺疾患の診断

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01151445A (ja) * 1987-12-08 1989-06-14 Toshiba Corp 超音波診断装置
JPH01151444A (ja) * 1987-12-08 1989-06-14 Toshiba Corp 超音波診断装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01151445A (ja) * 1987-12-08 1989-06-14 Toshiba Corp 超音波診断装置
JPH01151444A (ja) * 1987-12-08 1989-06-14 Toshiba Corp 超音波診断装置

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0556257A1 (en) * 1990-11-02 1993-08-25 The University Of Rochester Method and apparatus for breast imaging
EP0556257A4 (ja) * 1990-11-02 1994-08-31 The University Of Rochester
JP2002243711A (ja) * 2001-02-14 2002-08-28 Aloka Co Ltd 振動計測方法及び超音波顕微鏡システム
JP4672879B2 (ja) * 2001-02-14 2011-04-20 アロカ株式会社 振動計測方法及び超音波顕微鏡システム
JP2009531101A (ja) * 2006-03-29 2009-09-03 スーパー ソニック イマジン 粘弾性媒質を撮像する方法及び装置
CN103083045A (zh) * 2011-10-28 2013-05-08 Ge医疗系统环球技术有限公司 超声波诊断装置
JP2015504755A (ja) * 2012-01-26 2015-02-16 パルティ、ヨーラム 肺振動中の経胸郭肺ドップラー超音波を使用した肺疾患の診断
JP2014217605A (ja) * 2013-05-09 2014-11-20 国立大学法人群馬大学 映像装置及び映像形成方法
JP2014030774A (ja) * 2013-11-22 2014-02-20 Super Sonic Imagine 関心領域の粘弾性の平均値を測定するための方法および装置

Also Published As

Publication number Publication date
JP2699473B2 (ja) 1998-01-19

Similar Documents

Publication Publication Date Title
US5099848A (en) Method and apparatus for breast imaging and tumor detection using modal vibration analysis
Yamakoshi et al. Ultrasonic imaging of internal vibration of soft tissue under forced vibration
US5903516A (en) Acoustic force generator for detection, imaging and information transmission using the beat signal of multiple intersecting sonic beams
JPH06506055A (ja) 振幅推定用のドップラー変調パラメータの利用
US9717471B2 (en) Method and apparatus for multiple-wave doppler velocity meter
JPH0570769B2 (ja)
US7713201B2 (en) Method and apparatus for shear property characterization from resonance induced by oscillatory radiation force
Kelsey et al. Applications of ultrasound in speech research
KR20040066795A (ko) 신체 조직의 기계적 강도 특성을 나타내는 장치 및 방법
JPH02119851A (ja) 超音波検査装置
US10624610B2 (en) Ultrasonic imaging system
JP2629734B2 (ja) 超音波物体検査装置
Palanchon et al. New technique for emboli detection and discrimination based on nonlinear characteristics of gas bubbles
Nabavizadeh et al. Velocity measurement by vibro-acoustic Doppler
Hermawan et al. Color Doppler shear wave elastography using commercial ultrasound machine with compensated transducer scanning delay
JP2007309850A (ja) 軟質薄膜の物性値測定方法とそのための装置
JPS63194644A (ja) 超音波計測装置
JPH0147182B2 (ja)
JPH05168633A (ja) 超音波組織変位計測装置
JPS62227331A (ja) 超音波計測装置
JP6161028B2 (ja) 映像装置及び映像形成方法
JPS63134949A (ja) 超音波計測装置
JP2024073674A (ja) 体組織硬さの表示装置、表示方法
Baldwin et al. Influence of bright intramural echoes on estimates of ultrasonic attenuation from backscattered ultrasound in excised myocardium
JPH0428376B2 (ja)

Legal Events

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

Free format text: PAYMENT UNTIL: 20070926

Year of fee payment: 10

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

Free format text: PAYMENT UNTIL: 20080926

Year of fee payment: 11

LAPS Cancellation because of no payment of annual fees