JPH0510620B2 - - Google Patents
Info
- Publication number
- JPH0510620B2 JPH0510620B2 JP59038854A JP3885484A JPH0510620B2 JP H0510620 B2 JPH0510620 B2 JP H0510620B2 JP 59038854 A JP59038854 A JP 59038854A JP 3885484 A JP3885484 A JP 3885484A JP H0510620 B2 JPH0510620 B2 JP H0510620B2
- Authority
- JP
- Japan
- Prior art keywords
- sensor
- ultrasonic
- sensors
- array sensor
- ultrasonic waves
- 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
Links
Classifications
-
- 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/22—Details, e.g. general constructional or apparatus details
- G01N29/26—Arrangements for orientation or scanning by relative movement of the head and the sensor
- G01N29/262—Arrangements for orientation or scanning by relative movement of the head and the sensor by electronic orientation or focusing, e.g. with phased arrays
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/18—Methods or devices for transmitting, conducting or directing sound
- G10K11/26—Sound-focusing or directing, e.g. scanning
- G10K11/34—Sound-focusing or directing, e.g. scanning using electrical steering of transducer arrays, e.g. beam steering
- G10K11/341—Circuits therefor
- G10K11/345—Circuits therefor using energy switching from one active element to another
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Description
【発明の詳細な説明】
〔発明の利用分野〕
本発明は、金属中における超音波反射体を表示
するのに好適とされたアレイセンサによる超音波
送受信方法に係り、特に不感帯を生じせしめるこ
となく、しかも遅延回路不要にして鮮明な画像表
示が可能とされたアレイセンサによる超音波送受
信方法に関するものである。
するのに好適とされたアレイセンサによる超音波
送受信方法に係り、特に不感帯を生じせしめるこ
となく、しかも遅延回路不要にして鮮明な画像表
示が可能とされたアレイセンサによる超音波送受
信方法に関するものである。
これまでのアレイセンサの電子走査方法として
は各センサへ発信信号を遅延回路を介し印加する
ことによつて超音波を発信させる一方、超音波の
受信にあたつては各センサからの受信信号を遅延
回路を介し加算することが知られている。遅延回
路の遅延時間を各センサ単位に変えることで、各
センサから発信された超音波が干渉し、この結果
斜角超音波ビームや集束超音波ビームが実現され
得るものである。遅延時間を制御する場合は、超
音波ビームの進行方向および超音波の焦点距離を
変え得、焦点位置可変として焦点位置を走査し得
るわけである。このようにして被検体内部を走査
する場合は、固定焦点の凹面探触子を機械的に走
査して反射体を映像化する場合に比し10倍以上走
査時間が短縮され得ることから、そのような走査
方法は動いている反射体、たとえば内臓などの映
像化を目的とした医療機器の分野で採り入れられ
急速に発展しているのが実状であり、この他金属
中の欠陥検出にも応用が進んでいる。
は各センサへ発信信号を遅延回路を介し印加する
ことによつて超音波を発信させる一方、超音波の
受信にあたつては各センサからの受信信号を遅延
回路を介し加算することが知られている。遅延回
路の遅延時間を各センサ単位に変えることで、各
センサから発信された超音波が干渉し、この結果
斜角超音波ビームや集束超音波ビームが実現され
得るものである。遅延時間を制御する場合は、超
音波ビームの進行方向および超音波の焦点距離を
変え得、焦点位置可変として焦点位置を走査し得
るわけである。このようにして被検体内部を走査
する場合は、固定焦点の凹面探触子を機械的に走
査して反射体を映像化する場合に比し10倍以上走
査時間が短縮され得ることから、そのような走査
方法は動いている反射体、たとえば内臓などの映
像化を目的とした医療機器の分野で採り入れられ
急速に発展しているのが実状であり、この他金属
中の欠陥検出にも応用が進んでいる。
ここで、遅延回路を有してなるアレイセンサに
よる走査方法を第１図により説明すれば、アレイ
センサ１自体はセンサAi（ｉ＝１〜ｍ）が等間隔
ｄに配列されたものとしてなり、被検体２を走査
する場合にはアレイセンサ１は被検体２表面に設
置されるものとなつている。ここで座標軸を図示
の如くにとり、各センサAi，…Ai，…Anの座標
（ｘ，ｚ）をそれぞれ（x0，０），…，（x0＋（ｉ−
１）ｄ，０），…，（x0＋（ｍ−１）ｄ，０）とす
れば、センサAiからの超音波が全て位置Ｆに同時
に到達するように超音波を発信すれば、位置Ｆで
超音波が集束することになることは明らかであ
る。いま、位置Ｆの座標を例えば（x0，z0）とす
れば、センサAiから超音波を発信した時刻は超音
波が点Ｆに到達した時刻より時間tiだけ前になり
時間tiは次式によつて与えられることになる。
よる走査方法を第１図により説明すれば、アレイ
センサ１自体はセンサAi（ｉ＝１〜ｍ）が等間隔
ｄに配列されたものとしてなり、被検体２を走査
する場合にはアレイセンサ１は被検体２表面に設
置されるものとなつている。ここで座標軸を図示
の如くにとり、各センサAi，…Ai，…Anの座標
（ｘ，ｚ）をそれぞれ（x0，０），…，（x0＋（ｉ−
１）ｄ，０），…，（x0＋（ｍ−１）ｄ，０）とす
れば、センサAiからの超音波が全て位置Ｆに同時
に到達するように超音波を発信すれば、位置Ｆで
超音波が集束することになることは明らかであ
る。いま、位置Ｆの座標を例えば（x0，z0）とす
れば、センサAiから超音波を発信した時刻は超音
波が点Ｆに到達した時刻より時間tiだけ前になり
時間tiは次式によつて与えられることになる。
ti＝li／ｖ
＝√〔0−｛0＋（−１）｝2＋（０−
z0）2／ｖ ＝√（−１）2 2＋2 0／ｖ ……(1) 但し、ｖは被検体中での音速である。
z0）2／ｖ ＝√（−１）2 2＋2 0／ｖ ……(1) 但し、ｖは被検体中での音速である。
ｉの値によつて時間tiも変化することは明らか
であるが、第２図はその変化の様子を示したもの
である。使用されるセンサの数をｍとすれば時間
tnが最大値をとることになるものである。したが
つて、ｍ個の各センサから発信された超音波が全
て同時に点Ｆに到達するためには、ｍ番目のセン
サAnから超音波が発信された時点よりｉ番目の
センサAiからはτi＝tn−ti後に超音波を発信させ
なければならない。即ち、遅延時間τiは次のよう
に与えられる。
であるが、第２図はその変化の様子を示したもの
である。使用されるセンサの数をｍとすれば時間
tnが最大値をとることになるものである。したが
つて、ｍ個の各センサから発信された超音波が全
て同時に点Ｆに到達するためには、ｍ番目のセン
サAnから超音波が発信された時点よりｉ番目の
センサAiからはτi＝tn−ti後に超音波を発信させ
なければならない。即ち、遅延時間τiは次のよう
に与えられる。
τi＝tn−ti＝｛√（−１）2 2＋2 0−√（−１
）
2d2＋z2 0｝／ｖ ＝｛√（−１）2 2＋2 0−√（−１）2 2＋
z2 0｝／ｖ ……(2) このように遅延時間を制御して超音波を集束さ
せる場合には、使用されるセンサのうち何れのセ
ンサから最も早く超音波を発信させなくてはなら
ないかを決定する必要があり、また式(2)よりその
他のセンサから超音波を発信させる遅延時間を求
めるのが複雑になる。もちろん、アレイセンサで
は、集束位置Ｆをｘ方向，ｚ方向にも移動させる
ため、各センサの遅延時間および、それらの遅延
時間の組合せが膨大になり、また、遅延時間を生
みだす遅延回路は少なくともセンサ数だけは必要
になる。ちなみに遅延回路の価格は１個当り30万
円程度であり、アレイセンサのセンサ数を多くで
きない一つの理由ともなつている。
）
2d2＋z2 0｝／ｖ ＝｛√（−１）2 2＋2 0−√（−１）2 2＋
z2 0｝／ｖ ……(2) このように遅延時間を制御して超音波を集束さ
せる場合には、使用されるセンサのうち何れのセ
ンサから最も早く超音波を発信させなくてはなら
ないかを決定する必要があり、また式(2)よりその
他のセンサから超音波を発信させる遅延時間を求
めるのが複雑になる。もちろん、アレイセンサで
は、集束位置Ｆをｘ方向，ｚ方向にも移動させる
ため、各センサの遅延時間および、それらの遅延
時間の組合せが膨大になり、また、遅延時間を生
みだす遅延回路は少なくともセンサ数だけは必要
になる。ちなみに遅延回路の価格は１個当り30万
円程度であり、アレイセンサのセンサ数を多くで
きない一つの理由ともなつている。
ところで、金属の探傷にアレイセンサを適用す
る場合の問題点として次の２点が挙げられる。
る場合の問題点として次の２点が挙げられる。
(i) 不感帯の存在
金属中を伝播する超音波の速度は、人体中を伝
播する速度の約４倍である。すなわち、金属中の
波長は人体中での超音波波長の４倍になる。この
結果センサから近い位置では、各センサから発信
された超音波が干渉しきれず集束ビームを形成し
得ないというわけである。この限界距離は超音波
の波長に比例する。したがつて集束ビームを形成
し得ない領域は、金属探傷では人体の映像化に比
し４倍以上広がることになる。このため、アレイ
センサが設置される試料表面近傍に欠陥が存在す
る場合には、その欠陥は映像化され得ないという
わけである。
播する速度の約４倍である。すなわち、金属中の
波長は人体中での超音波波長の４倍になる。この
結果センサから近い位置では、各センサから発信
された超音波が干渉しきれず集束ビームを形成し
得ないというわけである。この限界距離は超音波
の波長に比例する。したがつて集束ビームを形成
し得ない領域は、金属探傷では人体の映像化に比
し４倍以上広がることになる。このため、アレイ
センサが設置される試料表面近傍に欠陥が存在す
る場合には、その欠陥は映像化され得ないという
わけである。
(ii) 遅延回路の制御
金属中での音速は人体中でのそれの約４倍であ
らることは既に述べたところであるが、医療用ア
レイセンサで所定の位置に焦点を結ばせる遅延時
間に比し金属中で同じ位置に焦点を形成させる遅
延時間の値は1/4になる。すなわち金属探傷では、
遅延時間を４倍高い精度で制御することが必要で
ある。
らることは既に述べたところであるが、医療用ア
レイセンサで所定の位置に焦点を結ばせる遅延時
間に比し金属中で同じ位置に焦点を形成させる遅
延時間の値は1/4になる。すなわち金属探傷では、
遅延時間を４倍高い精度で制御することが必要で
ある。
また、遅延回路はセンサの数と同じ数だけ必要
であり、さらに遅延回路各々は全センサにスイツ
チを介し選択的に連結可とされるべく、最低でも
センサ数の二乗分の結線が必要となつている。セ
ンサ自体が小型でも結線は膨大であり、また、遅
延回路の実装に際して制御部の小型化は困難なも
のとなつている。遅延回路としては通常コイルを
使用するが、製作上高度の技術が必要であること
から、コストが高くなること（約30万円／個）も
大きな問題となつている。
であり、さらに遅延回路各々は全センサにスイツ
チを介し選択的に連結可とされるべく、最低でも
センサ数の二乗分の結線が必要となつている。セ
ンサ自体が小型でも結線は膨大であり、また、遅
延回路の実装に際して制御部の小型化は困難なも
のとなつている。遅延回路としては通常コイルを
使用するが、製作上高度の技術が必要であること
から、コストが高くなること（約30万円／個）も
大きな問題となつている。
よつて本発明の目的は、被検体が金属であつて
も不感帯を生じせしめることなく、しかも遅延回
路不要にして鮮明な画像が表示され得るアレイセ
ンサによる超音波送受信方法を供するにある。
も不感帯を生じせしめることなく、しかも遅延回
路不要にして鮮明な画像が表示され得るアレイセ
ンサによる超音波送受信方法を供するにある。
この目的のため本発明は、ｎの値が順次更新さ
れるものとしてｎ個間隔毎に作動状態におかれた
センサから超音波を同時に被検体に向けて発信す
る一方、ｎの値が更新される毎に少なくともそれ
らセンサからの反射超音波受信信号を遅延させる
ことなく加算したうえ加算結果の波高値を記録す
ることによつて空間的なセンサ選択の周期に依存
したスペクトルを作成するようにし、作成された
スペクトルをフーリエ変換することによつて超音
波反射体の位置およびその形状を求め表示するよ
うにしたものである。
れるものとしてｎ個間隔毎に作動状態におかれた
センサから超音波を同時に被検体に向けて発信す
る一方、ｎの値が更新される毎に少なくともそれ
らセンサからの反射超音波受信信号を遅延させる
ことなく加算したうえ加算結果の波高値を記録す
ることによつて空間的なセンサ選択の周期に依存
したスペクトルを作成するようにし、作成された
スペクトルをフーリエ変換することによつて超音
波反射体の位置およびその形状を求め表示するよ
うにしたものである。
以下、本発明を説明する。
先ず本発明の理論的背景について第３図から第
８図により説明する。第３図に示す如く被検体２
中に点状反射体Ｐが存在し、その反射体Ｐから各
センサAiに反射波が入射する場合を想定すれば、
点Ｐの座標を（xp，zp）、センサAiのｚ軸座標を
０として点Ｐから伝播される反射波でセンサAiで
受信される波φiは以下のように表わされる。
８図により説明する。第３図に示す如く被検体２
中に点状反射体Ｐが存在し、その反射体Ｐから各
センサAiに反射波が入射する場合を想定すれば、
点Ｐの座標を（xp，zp）、センサAiのｚ軸座標を
０として点Ｐから伝播される反射波でセンサAiで
受信される波φiは以下のように表わされる。
φi＝（ｃ／li）・exp｛jωt−jkli｝ ……(3)
但し、ｊは複素記号（√−１）を、ωは超音波
の角速度を、ｋは超音波の波数を、ｃは定数をそ
れぞれ示す。
の角速度を、ｋは超音波の波数を、ｃは定数をそ
れぞれ示す。
第４図は式(3)に示す波φiの位相がセンサAiのｘ
座標によつて如何に変化するかを示したものであ
る。図に示すようにｘ座標xp，xqの位置では互い
に位相はπラジアンだけ異なつており、したがつ
てｘ座標xp，xqでの波は互いに反転した振動状態
にある。また、ｘ座標xp，xrでの位相は同じ値を
示している。即ち、ｘ座標xp，xrでの波の振動は
それぞれ同じ時間で同じ振動をする。しかして、
ｘ座標xp，xqでの振動を加算すると双方打ち消し
合つてしまう。これとは逆にｘ座標xp，xrの振動
を加算すると同位相なので振動の振幅は強め合つ
た振動になる。ここで、ｘ座標xpの振動を基準と
して任意のｘ座標ｘでの振動との加算を想定すれ
ば、先ずｘ座標xpでの振動φxpは式(3)より以下の
ようになる。
座標によつて如何に変化するかを示したものであ
る。図に示すようにｘ座標xp，xqの位置では互い
に位相はπラジアンだけ異なつており、したがつ
てｘ座標xp，xqでの波は互いに反転した振動状態
にある。また、ｘ座標xp，xrでの位相は同じ値を
示している。即ち、ｘ座標xp，xrでの波の振動は
それぞれ同じ時間で同じ振動をする。しかして、
ｘ座標xp，xqでの振動を加算すると双方打ち消し
合つてしまう。これとは逆にｘ座標xp，xrの振動
を加算すると同位相なので振動の振幅は強め合つ
た振動になる。ここで、ｘ座標xpの振動を基準と
して任意のｘ座標ｘでの振動との加算を想定すれ
ば、先ずｘ座標xpでの振動φxpは式(3)より以下の
ようになる。
φxp＝（ｃ／zp）・exp｛jωt−jkzp｝ ……(4)
同様にしてｘ座標ｘでの振動φxは以下のよう
に表わされる。
に表わされる。
φx＝ｃ／√（p−）2＋p 2・exp｛jωt−jk
√（p−）2＋p 2｝……(5) ここで｜xp−ｘ｜≪zpであると仮定すれば、式
(5)は以下のように変形され得る。
√（p−）2＋p 2｝……(5) ここで｜xp−ｘ｜≪zpであると仮定すれば、式
(5)は以下のように変形され得る。
φx≒〔ｃ／｛zp＋（xp−ｘ）2／2zp｝〕・exp〔
jωt−jk｛zp＋（xp−ｘ）2／2zp｝〕……(6) 第５図は式(4)，(5)の振動を加算した場合での振
動の振幅を示すが、図に示すようにｘ座標xpから
√2p離れたｘ座標位置の波を加算した場合に
振幅が高くなることが判る。但し、λは超音波の
波長である。したがつて、ｘ座標xp位置から√
2zpλの整数倍離れたｘ座標位置での波を全て加
算していけば、波の振幅が徐々に大きな値になる
ことは明らかである。
jωt−jk｛zp＋（xp−ｘ）2／2zp｝〕……(6) 第５図は式(4)，(5)の振動を加算した場合での振
動の振幅を示すが、図に示すようにｘ座標xpから
√2p離れたｘ座標位置の波を加算した場合に
振幅が高くなることが判る。但し、λは超音波の
波長である。したがつて、ｘ座標xp位置から√
2zpλの整数倍離れたｘ座標位置での波を全て加
算していけば、波の振幅が徐々に大きな値になる
ことは明らかである。
このことを数式的に説明すると以下のようにな
る。
る。
即ち、第６図に示す如く点Ｐを反射体位置と
し、ｘ軸上にアレイセンサが位置すると想定する
場合、ｘ座標ｘでの受信波Ψxはｘ≪z0のとき、
次のように近似し得る。
し、ｘ軸上にアレイセンサが位置すると想定する
場合、ｘ座標ｘでの受信波Ψxはｘ≪z0のとき、
次のように近似し得る。
Ψx＝〔ｃ／（z0＋x2／2z0）〕・exp〔jωt−jk（z
0＋x2／2z0）｝……(7) ここでｘ＝０での受信波をΨ0と表わすとすれ
ば式(7)は更に以下のように変形される。
0＋x2／2z0）｝……(7) ここでｘ＝０での受信波をΨ0と表わすとすれ
ば式(7)は更に以下のように変形される。
Ψx＝Ψ0（１−x2／2z0）・exp（−jkx2／2z0）…
…(8) 式(8)で与えられる振幅（実数部の値）を第７図
に示すが、このスペクトルは第５図で示すスペク
トルをxpだけ平行移動したスペクトルである。し
たがつて、反射体Ｐのｘ座標が如何なる位置でも
第５図，第７図で示すスペクトルの形状は保存さ
れ、最大振幅とその次に高い振幅との間隔は√
2zpλであることは明らかである。
…(8) 式(8)で与えられる振幅（実数部の値）を第７図
に示すが、このスペクトルは第５図で示すスペク
トルをxpだけ平行移動したスペクトルである。し
たがつて、反射体Ｐのｘ座標が如何なる位置でも
第５図，第７図で示すスペクトルの形状は保存さ
れ、最大振幅とその次に高い振幅との間隔は√
2zpλであることは明らかである。
次に、第７図で示すスペクトルの周波数分析を
行なうが、これは式(8)をフーリエ変換、即ち、以
下の式(9)によつてスペクトルを得るようにするも
のである。
行なうが、これは式(8)をフーリエ変換、即ち、以
下の式(9)によつてスペクトルを得るようにするも
のである。
式(9)の結果を図示すると第８図のようになる。
第８図に示す如くFxは幅の狭い正規分布を示す。
式(9)のFxは式(8)のΨxのフーリエ変換関数である
から、ΨxはFxのフーリエ変換関数でもある。Fx
は幅の狭い正規分布であるので、第６図に示す点
のｘ軸上の存在確率密度として考えることもでき
る。即ち、第７図で得られるスペクトルをフーリ
エ変換して反射体Ｐのｘ軸上位置を得る。ｚ軸位
置については、第７図の振幅極大値の間隔が√
2zpλであることから、これよりzpを得るもので
ある。このようにして反射体位置（xp，zp）が正
確に求められるわけである。
第８図に示す如くFxは幅の狭い正規分布を示す。
式(9)のFxは式(8)のΨxのフーリエ変換関数である
から、ΨxはFxのフーリエ変換関数でもある。Fx
は幅の狭い正規分布であるので、第６図に示す点
のｘ軸上の存在確率密度として考えることもでき
る。即ち、第７図で得られるスペクトルをフーリ
エ変換して反射体Ｐのｘ軸上位置を得る。ｚ軸位
置については、第７図の振幅極大値の間隔が√
2zpλであることから、これよりzpを得るもので
ある。このようにして反射体位置（xp，zp）が正
確に求められるわけである。
それでは、アレイセンサを用いた如何に走査す
れば第７図のようなスペクトルが得られるかを次
に説明する。
れば第７図のようなスペクトルが得られるかを次
に説明する。
第７図のスペクトルは式(8)のΨxを図示したも
のであり、また、ΨxはFxのフーリエ変換関数で
あることは既に述べた通りである。更に第７図と
第５図は同じスペクトルであることも既に述べた
ところである。結局第５図のスペクトルを得る方
法としては、複数あるセンサのうちから適当に複
数のセンサを選択したうえ選択されたセンサの受
信信号を加算した結果の振幅を測定すればよい。
この場合加算すべきセンサの間隔は等間隔にしな
くてはならない。これは、ΨxがFxのフーリエ変
換関数であるから、センサの間隔はｘ軸に対して
空間的な周波数に対応した値にするためである。
即ち、使用するセンサの間隔を１個，２個……，
ｉ個といつた具合にして、それらのセンサを選択
したうえ選択されたセンサの信号振幅を記録すれ
ば、そのスペクトルが第５あるいは第７図に示す
スペクトルとして得られるものである。
のであり、また、ΨxはFxのフーリエ変換関数で
あることは既に述べた通りである。更に第７図と
第５図は同じスペクトルであることも既に述べた
ところである。結局第５図のスペクトルを得る方
法としては、複数あるセンサのうちから適当に複
数のセンサを選択したうえ選択されたセンサの受
信信号を加算した結果の振幅を測定すればよい。
この場合加算すべきセンサの間隔は等間隔にしな
くてはならない。これは、ΨxがFxのフーリエ変
換関数であるから、センサの間隔はｘ軸に対して
空間的な周波数に対応した値にするためである。
即ち、使用するセンサの間隔を１個，２個……，
ｉ個といつた具合にして、それらのセンサを選択
したうえ選択されたセンサの信号振幅を記録すれ
ば、そのスペクトルが第５あるいは第７図に示す
スペクトルとして得られるものである。
以上説明した如く本発明では、アレイセンサに
おける複数のセンサのうちから１個，２個，…ｉ
個毎のセンサを順次選択使用し、それら選択使用
に係るセンサ全ての受信信号を加算した信号の振
幅を記録することで空間的なセンサ選択の周期に
依存したスペクトルを作成し、そのスペクトルを
フーリエ変換することによつて反射体位置を求め
たうえその像を表示することが可能となるもので
ある。このように本発明による場合は何等超音波
の送受信に遅延回路が要されることなく反射体の
像を表示可なわけである。
おける複数のセンサのうちから１個，２個，…ｉ
個毎のセンサを順次選択使用し、それら選択使用
に係るセンサ全ての受信信号を加算した信号の振
幅を記録することで空間的なセンサ選択の周期に
依存したスペクトルを作成し、そのスペクトルを
フーリエ変換することによつて反射体位置を求め
たうえその像を表示することが可能となるもので
ある。このように本発明による場合は何等超音波
の送受信に遅延回路が要されることなく反射体の
像を表示可なわけである。
さて、本発明を第９図から第１６図により具体
的に説明すれば以下のようである。
的に説明すれば以下のようである。
先ず本発明に係る超音波送受信装置について説
明する。第９図はその装置の構成を概略的に示し
たものである。これによると複数のセンサよりな
るアレイセンサ１にはスイツチ群３を介しパルス
発振器５より高電圧パルスが直接的に印加可とな
つている。センサ対応のスイツチがオン状態にあ
る場合はそのスイツチ対応のセンサからは超音波
が被検体に向けて発射され、また、超音波反射体
からの反射超音波はそのセンサ，スイツチ、更に
はセンサ対応のアイソレータ４および前置増幅器
６を介し他のオン状態にあるセンサからの超音波
受信信号と加算器７で加算されるようになつてい
るものである。ところで、何れのスイツチがオン
状態におかれるかはスイツチ制御器１３からのス
イツチ制御データによつている。スイツチ制御器
１３はパルス発振器５からのパルス発振同期信号
をカウントすることによつてパルス発振回数をカ
ウントしており、パルス発振回数に応じスイツチ
を所定にオン状態におくようになつているもので
ある。スイツチ制御器１３はまた使用されるセン
サ間隔に応じ空間位相信号を発生するが、これは
加算器７の出力に乗じられるべく乗算器８に与え
られるようになつている。なお、アイソレータ４
ではパルス発振器５からの高電圧パルスはしや断
され、加算器７では超音波受信信号のみが加算さ
れる。
明する。第９図はその装置の構成を概略的に示し
たものである。これによると複数のセンサよりな
るアレイセンサ１にはスイツチ群３を介しパルス
発振器５より高電圧パルスが直接的に印加可とな
つている。センサ対応のスイツチがオン状態にあ
る場合はそのスイツチ対応のセンサからは超音波
が被検体に向けて発射され、また、超音波反射体
からの反射超音波はそのセンサ，スイツチ、更に
はセンサ対応のアイソレータ４および前置増幅器
６を介し他のオン状態にあるセンサからの超音波
受信信号と加算器７で加算されるようになつてい
るものである。ところで、何れのスイツチがオン
状態におかれるかはスイツチ制御器１３からのス
イツチ制御データによつている。スイツチ制御器
１３はパルス発振器５からのパルス発振同期信号
をカウントすることによつてパルス発振回数をカ
ウントしており、パルス発振回数に応じスイツチ
を所定にオン状態におくようになつているもので
ある。スイツチ制御器１３はまた使用されるセン
サ間隔に応じ空間位相信号を発生するが、これは
加算器７の出力に乗じられるべく乗算器８に与え
られるようになつている。なお、アイソレータ４
ではパルス発振器５からの高電圧パルスはしや断
され、加算器７では超音波受信信号のみが加算さ
れる。
さて、乗算器８では、加算器７からの加算超音
波受信信号とスイツチ制御器１３からの空間位相
信号Cωとが乗算されるが、この結果は振幅検出
器９に出力されるようになつている。この場合空
間位相信号Cωの値は使用されるセンサの間隔を
ｉ個として次式で与えられる。
波受信信号とスイツチ制御器１３からの空間位相
信号Cωとが乗算されるが、この結果は振幅検出
器９に出力されるようになつている。この場合空
間位相信号Cωの値は使用されるセンサの間隔を
ｉ個として次式で与えられる。
Cω＝１／ｉ・exp（2πi／ｍ） ……(10)
振幅検出器９では乗算結果信号の振幅値が検出
され、記憶装置１０ではスイツチ制御データに対
応し、例えば2πi／ｍの値に比例した番地に検出
された振幅値を格納するようになつているもので
ある。フーリエ変換器１１では、スイツチ制御デ
ータにもとづきスイツチ間隔がｍ／２に到達した
時点で憶装置１０にあるスペクトルデータを読み
出したうえフーリエ変換するが、変換されたスペ
クトルは座標データとしてデイスプレイ１２に出
力されることによつて反射体の像が表示されるよ
うになつているものである。
され、記憶装置１０ではスイツチ制御データに対
応し、例えば2πi／ｍの値に比例した番地に検出
された振幅値を格納するようになつているもので
ある。フーリエ変換器１１では、スイツチ制御デ
ータにもとづきスイツチ間隔がｍ／２に到達した
時点で憶装置１０にあるスペクトルデータを読み
出したうえフーリエ変換するが、変換されたスペ
クトルは座標データとしてデイスプレイ１２に出
力されることによつて反射体の像が表示されるよ
うになつているものである。
ここで、スイツチ制御器１３によるスイツチ３
の制御について説明すれば、制御方法としては４
通りある。
の制御について説明すれば、制御方法としては４
通りある。
先ずその１は超音波は全てのセンサより発生さ
れるようにするが、反射超音波の受信は１個，２
個，…，ｎ個おきのセンサによつて行なう方法で
ある。受信時にスイツチをオン状態とする組合せ
の例を第１０図に示す。図中、白丸はオン状態
を、黒丸はオフ状態を表わしている。
れるようにするが、反射超音波の受信は１個，２
個，…，ｎ個おきのセンサによつて行なう方法で
ある。受信時にスイツチをオン状態とする組合せ
の例を第１０図に示す。図中、白丸はオン状態
を、黒丸はオフ状態を表わしている。
その２は第１１図に示す如くである。この場
合、発信回数がｉ（ｉ＋１）／２から （ｉ＋１）（ｉ＋２）／２−１までの間では記憶装置１ ０で同番地に記録してある値に新たに得た振幅値
を加算していくようにするものである。
合、発信回数がｉ（ｉ＋１）／２から （ｉ＋１）（ｉ＋２）／２−１までの間では記憶装置１ ０で同番地に記録してある値に新たに得た振幅値
を加算していくようにするものである。
その３は第１２図に示すようである。この場合
にはセンサの数は２の整数倍とされ、発振回数ｌ
が（ｉ／２−１）ｉ／２＋１からｉ／２（ｉ／２＋ １）の範囲では式(10)の代わりに次式(11)が用いられ
るものとなつている。
にはセンサの数は２の整数倍とされ、発振回数ｌ
が（ｉ／２−１）ｉ／２＋１からｉ／２（ｉ／２＋ １）の範囲では式(10)の代わりに次式(11)が用いられ
るものとなつている。
Cω＝exp（2πΔl／2i） ……(11)
但し、Δl＝ｌ−（ｉ／２−１）ｉ／２−１であ
り、この場合でも記憶装置１０では同番地に記憶
してある値には新たに得た振幅値が加算される。
り、この場合でも記憶装置１０では同番地に記憶
してある値には新たに得た振幅値が加算される。
その４は第１３図に示す如くである。この場合
はセンサ数ｍは２の整数乗とされ、発信回数が2i
−１から2i+1−２までの範囲では式(10)の代わりに
次式(12)が用いられる。
はセンサ数ｍは２の整数乗とされ、発信回数が2i
−１から2i+1−２までの範囲では式(10)の代わりに
次式(12)が用いられる。
Cω＝exp（2πΔl／2i） ……(12)
但し、Δl＝ｌ−2i＋１であり、この場合にも同
様に記憶装置１０では同番地に記憶してある値に
新たな発信回数で得た振幅値が加算記録されるも
のとなつている。
様に記憶装置１０では同番地に記憶してある値に
新たな発信回数で得た振幅値が加算記録されるも
のとなつている。
以上第１０図から第１３図に挙げた方法におい
ては、全センサで超音波を発生させることなく超
音波を受信すべきセンサのみ超音波を発生させる
ようにすることも可能である。
ては、全センサで超音波を発生させることなく超
音波を受信すべきセンサのみ超音波を発生させる
ようにすることも可能である。
次に発信回数がこれまでに述べたものよりも少
なくて済まされる方法について説明する。第１４
図はその方法に係る超音波送受信装置の構成を示
したものである。図示の如く第９図に示すものと
異なるところは、スイツチ群３の代わりに異なつ
た機能をもつスイツチ群３′が設けられ、また、
アイソレータ４′、前置増幅器６′、加算器７′お
よび減算器１４が新たに追加されていることであ
る。この場合アイソレータ４′はアイソレータ４
と同一機能を有しているが、スイツチ群３′はス
イツチ群３と異なり単なるオンオフ動作ではなく
アレイセンサ１とアイソレータ４，４′の何れか
との接続を切替えるためのものである。即ち、ス
イツチ制御器１３からのスイツチ制御データに従
つてあるスイツチはアレイセンサ１とアイソレー
タ４とを接続し、それ以外のスイツチはアレイセ
ンサ１とアイソレータ４′とを接続するものとな
つている。前置増幅器６，６′ではそれぞれアイ
ソレータ４，４′からの超音波受信信号が増幅さ
れるが、前置増幅器６，６′それぞれの出力は更
に加算器７，７′で加算されるようになつている。
減算器１４では加算器７の出力と加算器７′の出
力との間で減算が行なわれその結果は乗算器８に
与えられるようになつている。乗算器８以降の処
理内容は第９図の場合と同一である。したがつ
て、この方法による場合は、発振回数は例えば第
１２図に示す方法による場合に比し第１５図に示
す如くその1/2で済まされるものである。但し、
第１５図中、白丸はアレイセンサ１とアイソレー
タ４とのスイツチ接続状態を、三角印はアレイセ
ンサ１とアイソレータ４′とのスイツチ接続状態
をそれぞれ示す。なお、この方法による場合発信
回数ｌが（ｉ／２−１）４＋１からｉ／４（ｉ／
２＋１）の範囲では式(11)の代わりに次式(13)が用い
られる。
なくて済まされる方法について説明する。第１４
図はその方法に係る超音波送受信装置の構成を示
したものである。図示の如く第９図に示すものと
異なるところは、スイツチ群３の代わりに異なつ
た機能をもつスイツチ群３′が設けられ、また、
アイソレータ４′、前置増幅器６′、加算器７′お
よび減算器１４が新たに追加されていることであ
る。この場合アイソレータ４′はアイソレータ４
と同一機能を有しているが、スイツチ群３′はス
イツチ群３と異なり単なるオンオフ動作ではなく
アレイセンサ１とアイソレータ４，４′の何れか
との接続を切替えるためのものである。即ち、ス
イツチ制御器１３からのスイツチ制御データに従
つてあるスイツチはアレイセンサ１とアイソレー
タ４とを接続し、それ以外のスイツチはアレイセ
ンサ１とアイソレータ４′とを接続するものとな
つている。前置増幅器６，６′ではそれぞれアイ
ソレータ４，４′からの超音波受信信号が増幅さ
れるが、前置増幅器６，６′それぞれの出力は更
に加算器７，７′で加算されるようになつている。
減算器１４では加算器７の出力と加算器７′の出
力との間で減算が行なわれその結果は乗算器８に
与えられるようになつている。乗算器８以降の処
理内容は第９図の場合と同一である。したがつ
て、この方法による場合は、発振回数は例えば第
１２図に示す方法による場合に比し第１５図に示
す如くその1/2で済まされるものである。但し、
第１５図中、白丸はアレイセンサ１とアイソレー
タ４とのスイツチ接続状態を、三角印はアレイセ
ンサ１とアイソレータ４′とのスイツチ接続状態
をそれぞれ示す。なお、この方法による場合発信
回数ｌが（ｉ／２−１）４＋１からｉ／４（ｉ／
２＋１）の範囲では式(11)の代わりに次式(13)が用い
られる。
Cω＝exp（πΔl／2i） ……(13)
但し、Δl＝ｌ−（ｉ／２−１）ｉ／４−１であ
る。
る。
同様に第１３図に示す方法による場合に比し発
振回数はその場合での1/2で済まされることにな
る。
振回数はその場合での1/2で済まされることにな
る。
第１６図は第１４図に示す構成よりさらに詳細
な振幅を検出するための装置構成を示す。アレイ
センサ１、パルス発振器５、スイツチ制御器１
３、乗算器８、振幅検出器９、記憶装置１０、フ
ーリエ変換器１１およびデイスプレイ１２の構成
とその機能は第１４図に示すものと変わらない
が、受信信号をアイソレートする組合せを変えた
構成にしているのが特徴である。これによるとア
イソレータ４，４′，４″，４は、スイツチ３′，
３″，３，３′′′′を介してアレイセンサ１に対し
それぞれセンサ１個分ずつシフトして連結されて
いる。アイソレータ群の数はセンサ個数より少な
い数で、各アイソレータ群の結線のシフト量も２
個以上でも構わない。アイソレータ群の出力はそ
れぞれ前置増幅器群６，６′，６″，６で増幅さ
れ、それら増幅信号は加算器７，７′，７″，７
で加算されるようになつている。これらの加算信
号は乗算器８′，８″，８で重みをつけられたう
え加算器７′′′′で更に加算されるものとなつてい
る。乗算器８′，８″，８での重みは、スイツチ
３′の結線に対するそれぞれのスイツチが結線し
ているセンサのシフト分に相当する量であり、シ
フト分をｎ、スイツチ３′で使用されるセンサの
間隔をｉとすると、重みWaは以下のようにな
る。
な振幅を検出するための装置構成を示す。アレイ
センサ１、パルス発振器５、スイツチ制御器１
３、乗算器８、振幅検出器９、記憶装置１０、フ
ーリエ変換器１１およびデイスプレイ１２の構成
とその機能は第１４図に示すものと変わらない
が、受信信号をアイソレートする組合せを変えた
構成にしているのが特徴である。これによるとア
イソレータ４，４′，４″，４は、スイツチ３′，
３″，３，３′′′′を介してアレイセンサ１に対し
それぞれセンサ１個分ずつシフトして連結されて
いる。アイソレータ群の数はセンサ個数より少な
い数で、各アイソレータ群の結線のシフト量も２
個以上でも構わない。アイソレータ群の出力はそ
れぞれ前置増幅器群６，６′，６″，６で増幅さ
れ、それら増幅信号は加算器７，７′，７″，７
で加算されるようになつている。これらの加算信
号は乗算器８′，８″，８で重みをつけられたう
え加算器７′′′′で更に加算されるものとなつてい
る。乗算器８′，８″，８での重みは、スイツチ
３′の結線に対するそれぞれのスイツチが結線し
ているセンサのシフト分に相当する量であり、シ
フト分をｎ、スイツチ３′で使用されるセンサの
間隔をｉとすると、重みWaは以下のようにな
る。
Wa＝cos（2πn／ｉ） ……(14)
シフト分ｎ、センサ間隔ｉはスイツチ制御器１
３から重み変換器２０に送信され、変換器２０は
式(14)にもとづき重みが算出されたうえそれぞれ乗
算器８′，８″，８に出力されるようになつてい
る。加算器７′′′′で重み付けの加算を行つた後の
信号処理は第１４図の場合と同じである。この構
成では、反射波の受信に関しアイソレータ群で、
１回の発信で並列に受信信号処理を行なうため、
その分発信回数を少なくできること、重み付けし
た信号の加算により加算信号の振幅が精度よく検
出できることになる。
３から重み変換器２０に送信され、変換器２０は
式(14)にもとづき重みが算出されたうえそれぞれ乗
算器８′，８″，８に出力されるようになつてい
る。加算器７′′′′で重み付けの加算を行つた後の
信号処理は第１４図の場合と同じである。この構
成では、反射波の受信に関しアイソレータ群で、
１回の発信で並列に受信信号処理を行なうため、
その分発信回数を少なくできること、重み付けし
た信号の加算により加算信号の振幅が精度よく検
出できることになる。
以上説明したように本発明は、複数のセンサか
らは超音波が同時に被検体に向けて発信されるよ
うにするが、被検体からの反射超音波を受信する
際には少なくとも受信すべきセンサからの反射超
音波受信信号は加算されたうえその波高値が記録
され、それにより作成された、空間的なセンサ選
択の周期に依存したスペクトルをフーリエ変換す
ることによつて超音波反射体位置を求め表示する
ようにしたものである。したがつて、本発明によ
る場合は、超音波送受信の際にこれまで要されて
いた高価な遅延回路が不要とされながらも鮮明な
画像表示が可能となるばかりか、装置の小型化や
経済性の面で有利となる。また、集束ビーム形成
の際アレイセンサ近傍領域においては超音波ビー
ムの音圧が不安定となり、その領域内に反射体が
存在しても十分な強度でそれを検出し得ないが、
本発明による場合は集束ビームは形成されなく、
したがつて、被検体が金属であつてもアレイセン
サに近い領域における超音波反射体を十分な強度
で表示し得ることになる。
らは超音波が同時に被検体に向けて発信されるよ
うにするが、被検体からの反射超音波を受信する
際には少なくとも受信すべきセンサからの反射超
音波受信信号は加算されたうえその波高値が記録
され、それにより作成された、空間的なセンサ選
択の周期に依存したスペクトルをフーリエ変換す
ることによつて超音波反射体位置を求め表示する
ようにしたものである。したがつて、本発明によ
る場合は、超音波送受信の際にこれまで要されて
いた高価な遅延回路が不要とされながらも鮮明な
画像表示が可能となるばかりか、装置の小型化や
経済性の面で有利となる。また、集束ビーム形成
の際アレイセンサ近傍領域においては超音波ビー
ムの音圧が不安定となり、その領域内に反射体が
存在しても十分な強度でそれを検出し得ないが、
本発明による場合は集束ビームは形成されなく、
したがつて、被検体が金属であつてもアレイセン
サに近い領域における超音波反射体を十分な強度
で表示し得ることになる。
第１図は、遅延回路を有してなるアレイセンサ
による被検体走査方法を説明するための図、第２
図は、特定位置に超音波をアレイセンサによつて
集束させる場合での超音波発生時点がセンサ位置
によつて如何に変化するかを示す図、第３図、第
４図、第５図、第６図、第７図、第８図は、本発
明の理論的背景あるいは原理を説明するための
図、第９図は、本発明に係る超音波送受信装置の
構成を概略的に示す図、第１０図、第１１図、第
１２図、第１３図は、その装置による超音波送受
信方法を説明するための図、第１４図は、本発明
に係る超音波送受信装置の他の例での構成を概略
的に示す図、第１５図は、その装置による超音波
送受信方法を説明するための図、第１６図は、詳
細な振幅検出が可とされた本発明に係る超音波送
受信装置の一例での構成を示す図である。 １…アレイセンサ、３，３′，３″，３，３
′′′′…スイツチ群、４，４′，４″，４…アイソ
レータ、５…パルス発振器、６，６′，６″，６
…前置増幅器、７，７′，７″，７…加算器、８
…乗算器、９…振幅検出器、１０…記憶装置、１
１…フーリエ変換器、１２…デイスプレイ、１３
…スイツチ制御器。
による被検体走査方法を説明するための図、第２
図は、特定位置に超音波をアレイセンサによつて
集束させる場合での超音波発生時点がセンサ位置
によつて如何に変化するかを示す図、第３図、第
４図、第５図、第６図、第７図、第８図は、本発
明の理論的背景あるいは原理を説明するための
図、第９図は、本発明に係る超音波送受信装置の
構成を概略的に示す図、第１０図、第１１図、第
１２図、第１３図は、その装置による超音波送受
信方法を説明するための図、第１４図は、本発明
に係る超音波送受信装置の他の例での構成を概略
的に示す図、第１５図は、その装置による超音波
送受信方法を説明するための図、第１６図は、詳
細な振幅検出が可とされた本発明に係る超音波送
受信装置の一例での構成を示す図である。 １…アレイセンサ、３，３′，３″，３，３
′′′′…スイツチ群、４，４′，４″，４…アイソ
レータ、５…パルス発振器、６，６′，６″，６
…前置増幅器、７，７′，７″，７…加算器、８
…乗算器、９…振幅検出器、１０…記憶装置、１
１…フーリエ変換器、１２…デイスプレイ、１３
…スイツチ制御器。
Claims (1)
- 【特許請求の範囲】 １ 複数のセンサが等間隔に配置されてなるアレ
イセンサによる超音波送受信方法にして、ｎの値
が更新される程にｎ個間隔毎に作動状態におかれ
るセンサより少なくとも超音波を同時に被検体に
向けて発信する一方、超音波発信の度に超音波発
信に係るセンサからの反射超音波受信信号を遅延
させることなく加算したうえ加算結果の波高値を
記録することによつて空間的なセンサ選択の周期
に依存した空間周波数スペクトルを作成し、該ス
ペクトルをフーリエ変換により実空間強度分布に
変換した上、被検体内における超音波反射体の位
置および形状を表示することを特徴とするアレイ
センサによる超音波送受信方法。 ２ ｎは2k（ｋ；正の整数）にしたがつて更新さ
れる特許請求の範囲第１項記載のアレイセンサに
よる超音波送受信方法。 ３ 超音波を送受信するセンサには該センサに対
してｊ（ｊ＜ｎ）個分間隔をずらしたセンサが含
まれ、ｎ個間隔のセンサ群と、ｊ個分センサ間隔
がずれたセンサ群との超音波受信信号は重みをつ
けて加算される特許請求の範囲第１項記載のアレ
イセンサによる超音波送受信方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59038854A JPS60183553A (ja) | 1984-03-02 | 1984-03-02 | アレイセンサによる超音波送受信方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59038854A JPS60183553A (ja) | 1984-03-02 | 1984-03-02 | アレイセンサによる超音波送受信方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS60183553A JPS60183553A (ja) | 1985-09-19 |
| JPH0510620B2 true JPH0510620B2 (ja) | 1993-02-10 |
Family
ID=12536786
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP59038854A Granted JPS60183553A (ja) | 1984-03-02 | 1984-03-02 | アレイセンサによる超音波送受信方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS60183553A (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH06100584B2 (ja) * | 1985-12-16 | 1994-12-12 | 株式会社日立メデイコ | 超音波撮像装置の送受波回路 |
| CN105078472B (zh) * | 2015-08-24 | 2018-10-12 | 联想(北京)有限公司 | 生理特征图像获取装置和方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5396285A (en) * | 1977-02-02 | 1978-08-23 | Hitachi Medical Corp | Ultrasonic diagnosing device |
| DE3024995A1 (de) * | 1980-07-02 | 1982-01-28 | Philips Patentverwaltung Gmbh, 2000 Hamburg | Ultraschall-untersuchungsanordnung |
-
1984
- 1984-03-02 JP JP59038854A patent/JPS60183553A/ja active Granted
Also Published As
| Publication number | Publication date |
|---|---|
| JPS60183553A (ja) | 1985-09-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH07303644A (ja) | 超音波診断装置 | |
| JPS61185259A (ja) | 超音波エコーグラフイによる物体の検差装置 | |
| KR0180057B1 (ko) | 초음파시스템의 3차원 영상 획득장치 | |
| US4702112A (en) | Ultrasonic phase reflectoscope | |
| JPH03500454A (ja) | 人為構造を除外した超音波反射伝送映像化方法および装置 | |
| JPS6111659A (ja) | 超音波検査装置 | |
| JPH0518942A (ja) | Ｖ（ｚ）特性による超音波音速測定装置およびこれを用いた超音波顕微鏡 | |
| JP6819184B2 (ja) | 超音波測定装置および超音波測定装置の制御方法 | |
| CN114305495B (zh) | 基于超声换能器的超声成像方法、超声换能器及系统 | |
| JP2005043164A (ja) | 音波伝播時間測定装置 | |
| JPH0510620B2 (ja) | ||
| CN113925528B (zh) | 多普勒成像方法和超声设备 | |
| WO1996024053A1 (en) | Acoustic probing apparatus | |
| JPH0592001A (ja) | 超音波診断装置 | |
| JPH0364831B2 (ja) | ||
| JPH0526655A (ja) | 膜厚測定方法及び装置 | |
| CN105548363A (zh) | 基于多途识别的超声检测成像方法 | |
| JPH05149931A (ja) | 音速・密度測定方法及び装置 | |
| JP3136635B2 (ja) | 探査方法およびその装置 | |
| JP3002065U (ja) | 多重反射波連続自動消滅装置 | |
| KR100543736B1 (ko) | 초음파 영상 표시 방법 | |
| JPH0467857A (ja) | 超音波診断装置 | |
| JPS597260A (ja) | 超音波探傷方法および装置 | |
| Kozuka et al. | Visualization of sound source positions using pinpoint nonlinear secondary emission by ultrasound focus scanning | |
| JP2854730B2 (ja) | 超音波可視化装置 |