JPH0357907A - 物体の2つの面を分離している距離eの測定方法 - Google Patents
物体の2つの面を分離している距離eの測定方法Info
- Publication number
- JPH0357907A JPH0357907A JP2190927A JP19092790A JPH0357907A JP H0357907 A JPH0357907 A JP H0357907A JP 2190927 A JP2190927 A JP 2190927A JP 19092790 A JP19092790 A JP 19092790A JP H0357907 A JPH0357907 A JP H0357907A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- sensor
- frequency
- thickness
- absolute value
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims description 55
- 238000001228 spectrum Methods 0.000 claims abstract description 13
- 230000009466 transformation Effects 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 13
- 239000002131 composite material Substances 0.000 claims description 11
- 238000002604 ultrasonography Methods 0.000 claims description 10
- 230000000737 periodic effect Effects 0.000 claims description 9
- 230000002123 temporal effect Effects 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 230000001902 propagating effect Effects 0.000 claims description 2
- 210000001367 artery Anatomy 0.000 abstract description 9
- 238000012545 processing Methods 0.000 abstract description 5
- 238000006243 chemical reaction Methods 0.000 abstract description 2
- 239000010409 thin film Substances 0.000 abstract description 2
- 238000002592 echocardiography Methods 0.000 description 11
- 238000012937 correction Methods 0.000 description 10
- 238000005259 measurement Methods 0.000 description 10
- 235000006693 Cassia laevigata Nutrition 0.000 description 4
- 241000522641 Senna Species 0.000 description 4
- 229940124513 senna glycoside Drugs 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 239000010410 layer Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 2
- 239000010931 gold Substances 0.000 description 2
- 229910052737 gold Inorganic materials 0.000 description 2
- 238000000338 in vitro Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 210000003323 beak Anatomy 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 150000001768 cations Chemical class 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000003973 paint Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0858—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B17/00—Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations
- G01B17/02—Measuring arrangements characterised by the use of infrasonic, sonic or ultrasonic vibrations for measuring thickness
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices Characterised By Use Of Acoustic Means (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〔発明の利用分野〕
本発明は、物体の2つの面を分離している距離eの測定
方法に関している。この方法は、物体の方向に速度Cで
伝播する超音波を発信することと、上記物体の面からの
上記超音波の反射によシ発生されたエコーを七ンサによ
り受信することを含んでいる。上記七ンサは、時間の関
数として変化する振幅を有する複合電気信号g(t)供
給する。
方法に関している。この方法は、物体の方向に速度Cで
伝播する超音波を発信することと、上記物体の面からの
上記超音波の反射によシ発生されたエコーを七ンサによ
り受信することを含んでいる。上記七ンサは、時間の関
数として変化する振幅を有する複合電気信号g(t)供
給する。
物体の2つの面の間の距離を測定するのに超音波を使用
している方法は、既にいくつか提供されている。2つの
面の間の距離が大きい場合には、エコーの形で受信され
る2つの信号金翻訳することは一般にあtb困難ではな
い。なぜならば、面は、時間的にはつきbと離間した2
つのピークを発生するので、音波が伝播する媒体中の音
波の伝播速度を知れば、2つの面の間の問題となってい
る距離を計算することができるからである。しかし、上
記面が非常に接近し、波長程度の距離しか離れていない
場合、それらが発生するエコーは一体化しているように
見える。
している方法は、既にいくつか提供されている。2つの
面の間の距離が大きい場合には、エコーの形で受信され
る2つの信号金翻訳することは一般にあtb困難ではな
い。なぜならば、面は、時間的にはつきbと離間した2
つのピークを発生するので、音波が伝播する媒体中の音
波の伝播速度を知れば、2つの面の間の問題となってい
る距離を計算することができるからである。しかし、上
記面が非常に接近し、波長程度の距離しか離れていない
場合、それらが発生するエコーは一体化しているように
見える。
基板にコードされた5〜20μmのコーティングの厚さ
の測定は、90〜5 1 0 MHzの周波数を用いた
超音波干渉によシ行なわれることが多かった。
の測定は、90〜5 1 0 MHzの周波数を用いた
超音波干渉によシ行なわれることが多かった。
このような方法は、1984年1月. J. Appl
. 55(1)の第194〜198頁に示されている
「超音波干渉による薄膜の厚さの測定」と題する論文に
示されている。しかし、この方法は、物体が安定してい
る場合、すなわち測定時に物体の次元が変化しない場合
にのみ満足のいく結果を得ることができる。
. 55(1)の第194〜198頁に示されている
「超音波干渉による薄膜の厚さの測定」と題する論文に
示されている。しかし、この方法は、物体が安定してい
る場合、すなわち測定時に物体の次元が変化しない場合
にのみ満足のいく結果を得ることができる。
IEEE超音波シンポジウム1988年の第1027〜
1029頁に示されている「多層媒体の厚さ測定に関す
るデイジタル信号処理方法」と題する論文は、自動車産
業に訃ける塗料の薄い層に適用される超音波測定につい
て述べている。この場合、層の厚さは20〜110μm
である。この場合の中心周波数は約75MHzで帯域幅
は100MHzである。
1029頁に示されている「多層媒体の厚さ測定に関す
るデイジタル信号処理方法」と題する論文は、自動車産
業に訃ける塗料の薄い層に適用される超音波測定につい
て述べている。この場合、層の厚さは20〜110μm
である。この場合の中心周波数は約75MHzで帯域幅
は100MHzである。
この測定方法は、無線周波数信号のパヮー●スペクトル
密度の対数の逆フーリエ変換と定義されるパワー・ケプ
ストラムの分析に基いている。人が表面層の目の前にい
る場合には、高周波で使用することができ、そのため優
れた解像度を得ることができる。しかし、この方法は、
超音波金減衰させる媒体中に位置する物体の厚さの測定
には適して釦らず、したがって、本発明の場合と同様に
、比較的低い周波数(xoMHz)において使用しなけ
ればならない。
密度の対数の逆フーリエ変換と定義されるパワー・ケプ
ストラムの分析に基いている。人が表面層の目の前にい
る場合には、高周波で使用することができ、そのため優
れた解像度を得ることができる。しかし、この方法は、
超音波金減衰させる媒体中に位置する物体の厚さの測定
には適して釦らず、したがって、本発明の場合と同様に
、比較的低い周波数(xoMHz)において使用しなけ
ればならない。
このケプストラム方法は、筐た、グルノーブルのナショ
ナル●ポリテクニカル学会においてジェイ●シー●パル
トによ#)1979年2月9日に提出された@Leg
operatcurs cepstres+appli
cationm a la separation d
’@chos rapproeha〆と題する論文にも
示されている。しかし、この論文の方法を非常に接近し
たエコーに対して使用した場合、翻訳が困難で、しかも
観測すべきピークは、見ることのできる限界に位置して
いる。
ナル●ポリテクニカル学会においてジェイ●シー●パル
トによ#)1979年2月9日に提出された@Leg
operatcurs cepstres+appli
cationm a la separation d
’@chos rapproeha〆と題する論文にも
示されている。しかし、この論文の方法を非常に接近し
たエコーに対して使用した場合、翻訳が困難で、しかも
観測すべきピークは、見ることのできる限界に位置して
いる。
超音波画像6 . 1984年,第435〜451頁の
「高速ヂコンボリューション方法によるレンジ解像度の
改善」と題する論文には、生きている組織である生体に
適用される超音波エコグラフイ方法が示されている。数
学的物理学の多くの問題は、y=h●Xの形のコンポリ
ューションの式で表すことができる。xkよびyは入力
釦よび出力の方式を示し、hはその方式のパルス・レス
ポンスで6る。デコンポリューションは、実験的に測定
された出力に基いて入力を再構成するという逆の問題で
ある。この論文に示された方法の釦かげで、眼底を精密
に調べることができ、さらに心臓の循環時にかける動脈
壁の厚さの変化を追うことができる。しかし、ここに示
された方法によれば、時間的に非常に接近しているエコ
ー金分離することはかなり困難である、という欠点を有
している。
「高速ヂコンボリューション方法によるレンジ解像度の
改善」と題する論文には、生きている組織である生体に
適用される超音波エコグラフイ方法が示されている。数
学的物理学の多くの問題は、y=h●Xの形のコンポリ
ューションの式で表すことができる。xkよびyは入力
釦よび出力の方式を示し、hはその方式のパルス・レス
ポンスで6る。デコンポリューションは、実験的に測定
された出力に基いて入力を再構成するという逆の問題で
ある。この論文に示された方法の釦かげで、眼底を精密
に調べることができ、さらに心臓の循環時にかける動脈
壁の厚さの変化を追うことができる。しかし、ここに示
された方法によれば、時間的に非常に接近しているエコ
ー金分離することはかなり困難である、という欠点を有
している。
米国特許第4.370.985号には、動脈の内径を測
定できる非侵入形のセンナが示されている。しかし、こ
こに示された方法では、動脈の実際の厚さを測定するこ
とはできない。
定できる非侵入形のセンナが示されている。しかし、こ
こに示された方法では、動脈の実際の厚さを測定するこ
とはできない。
前述したように、現在使用されている方法は、媒体中で
測定しなければならない場合には適しておらず、1た測
定する厚さが非常に薄い場合翻訳するのが困難である。
測定しなければならない場合には適しておらず、1た測
定する厚さが非常に薄い場合翻訳するのが困難である。
これら問題を解決するため、本発明は、センサにより受
信される複合電気信号に次のような信号処理過程を施す
。
信される複合電気信号に次のような信号処理過程を施す
。
a)センナからの複合電気信号g (t)を、周波数の
関数である振幅金有しかつ実数部R @ C G(f)
:) >よび虚数部1ml:G(f)”)を有する周波
数スペクトルG(f)に変換する過程、 b)周波数スペクトルl a(r) l = ( Re
C a(f) :]2十Irn( G(f)〕2}n
/2)n/2のn次の絶対値を計算する過程、 C)複合電気信号からエコーを分離するため、過程b)
にシいて計算された絶対値の対数導関数( IQ(f)
l )’/|G(f)l を計算する過程で、上記
対数導関数は、超音波周波数に実質的に集中した有効部
分と上記有効部分の両側に位置するノイズを有する部分
とから成り、上記有効部分は、実質的線形信号が重畳さ
れた周期信号から戒っている上記過程、 d)その有効部分を分離しかつノイズを有する部分をな
くすため、対数導関数に周波数の窓を供給する過程、 ●)少iくとも上記有効部分の周期信号に過程&)の変
換の逆の変換を行なう過程、 f)時間に関して離間した一連の可変振幅ピークを有す
る、逆変換の包結線を決定する過程、g)有効部分から
の周期信号により、最大振幅のピークから時間的原点1
での時間的拡がりJtを決定する過程。
関数である振幅金有しかつ実数部R @ C G(f)
:) >よび虚数部1ml:G(f)”)を有する周波
数スペクトルG(f)に変換する過程、 b)周波数スペクトルl a(r) l = ( Re
C a(f) :]2十Irn( G(f)〕2}n
/2)n/2のn次の絶対値を計算する過程、 C)複合電気信号からエコーを分離するため、過程b)
にシいて計算された絶対値の対数導関数( IQ(f)
l )’/|G(f)l を計算する過程で、上記
対数導関数は、超音波周波数に実質的に集中した有効部
分と上記有効部分の両側に位置するノイズを有する部分
とから成り、上記有効部分は、実質的線形信号が重畳さ
れた周期信号から戒っている上記過程、 d)その有効部分を分離しかつノイズを有する部分をな
くすため、対数導関数に周波数の窓を供給する過程、 ●)少iくとも上記有効部分の周期信号に過程&)の変
換の逆の変換を行なう過程、 f)時間に関して離間した一連の可変振幅ピークを有す
る、逆変換の包結線を決定する過程、g)有効部分から
の周期信号により、最大振幅のピークから時間的原点1
での時間的拡がりJtを決定する過程。
なお物体の2つの而の間の距離eは、e=c・At/2
により表わされる。
により表わされる。
以下、添付の図面に基いて、本発明の実施例に関し説明
する。
する。
tlcl図は、壁2の厚さを測定される物体1の断面を
示している。壁は、外而3と内面1を有している。この
測定を行なうため、超音波七ンサのエミッターレシーバ
6により超音波5が発信される。
示している。壁は、外而3と内面1を有している。この
測定を行なうため、超音波七ンサのエミッターレシーバ
6により超音波5が発信される。
各超音波は、速度Cで物体の方向に伝播される。
発信の際、センナは発振器7Kより励起され、その周波
数は10MHzのオーダである。本実施例では、センサ
6ぱ、物体の面3,4からの超音波の反射により生じた
エコーのレシーバとしても作用する。このセンサは、時
間tの関数として変化する振幅Aの複合電気信号g(t
)を供給する。この電気{i号は、本発明の方法にした
がって処理される。
数は10MHzのオーダである。本実施例では、センサ
6ぱ、物体の面3,4からの超音波の反射により生じた
エコーのレシーバとしても作用する。このセンサは、時
間tの関数として変化する振幅Aの複合電気信号g(t
)を供給する。この電気{i号は、本発明の方法にした
がって処理される。
全ての処理は、第1図に示したパーソナル●コンピュー
タ8(PC)によって行なわれる。この方法の各処理過
程において、ディスプレイ9により、実施された計算を
見ることができる。後述するように、本発明の方法は、
全て、8087形数理コブロセツサt加tたIBM A
Tパーソナル●コンピュータによって実施することがで
きる。ここではその一例を示したにすぎず、他の装置を
使用してもよい。
タ8(PC)によって行なわれる。この方法の各処理過
程において、ディスプレイ9により、実施された計算を
見ることができる。後述するように、本発明の方法は、
全て、8087形数理コブロセツサt加tたIBM A
Tパーソナル●コンピュータによって実施することがで
きる。ここではその一例を示したにすぎず、他の装置を
使用してもよい。
第2図の流れ図は、それぞれ、本発明の方法の1つの操
作を示した様々なブロックを示している。
作を示した様々なブロックを示している。
ここには完全な配置が示されているが、以下に説明する
ように、本発明の思想から離れることなく、いくつかの
ブロックを省略することができる。
ように、本発明の思想から離れることなく、いくつかの
ブロックを省略することができる。
発振器7は、約10MHzの中心周波数を有する超音波
を発信する圧電センサ6を励起する。この同じセンナで
、物体1の壁2の面3,4から反射して戻ってきた信号
g(t)が集められる。この信号は、時間tの関数とし
て変化する振幅Aを有し、センサ6の出力に釦いてはア
ナログ形態で示されている。増幅器10にかける増幅後
、アナログーデイジタル●コンバータ11により信号g
(t)をデイジタル化することは、必ずしも必要では
ないが都合が良い。ここに示されている例では、信号は
100MHzで8ビットにわたシサンプルされ、その後
、2048ポイントを有するメモリに記憶される。
を発信する圧電センサ6を励起する。この同じセンナで
、物体1の壁2の面3,4から反射して戻ってきた信号
g(t)が集められる。この信号は、時間tの関数とし
て変化する振幅Aを有し、センサ6の出力に釦いてはア
ナログ形態で示されている。増幅器10にかける増幅後
、アナログーデイジタル●コンバータ11により信号g
(t)をデイジタル化することは、必ずしも必要では
ないが都合が良い。ここに示されている例では、信号は
100MHzで8ビットにわたシサンプルされ、その後
、2048ポイントを有するメモリに記憶される。
本発明の方法を説明するためにここで使用されている物
体は、プラスチックの管で、その壁の厚さを測定する。
体は、プラスチックの管で、その壁の厚さを測定する。
この厚さを、たとえば、マイクロメータのような他の周
知の手段により測定することができるので、2つの方式
の測定結果を照らし合せて、本発明の方法により得られ
た結果を立証することができる。
知の手段により測定することができるので、2つの方式
の測定結果を照らし合せて、本発明の方法により得られ
た結果を立証することができる。
第3図,第4図,第5図のグラフは、それぞれ、3つの
数値の管理の厚さに関する、デイジタル化した後の信号
g (t)を示している。これらグラフは、縦座標に信
号の振幅(相対値)と横座標に時間(マイクロ秒)を示
している。
数値の管理の厚さに関する、デイジタル化した後の信号
g (t)を示している。これらグラフは、縦座標に信
号の振幅(相対値)と横座標に時間(マイクロ秒)を示
している。
第3図は、壁の厚さ@=0.7511の場合の信号g
(t)を示している。ここでは、壁の外面3と内面4か
ら戻ってきた2つのエコー20.21をはつきbと見分
けると七ができる。信号の2つの主なピークを分離して
いるグラフにおける時間Jt を測定すれば、at=i
μSであることがわかる。さらに、この場合、音波の伝
播速度Cが1. 5 uv’μmとわかっているので、
容易に壁の厚さ●=CIIAt/2を計算することがで
きる。2で割ることは、音波の往復を考慮している。式
のCとJtを1.5カよびlにそれぞれ置き換えれば、
厚さ●の値0.75nが得られる。このように、エコー
かはつきbと分離されている場合には、厚さの測定が簡
単にでき、本発明の方法は必要ない。
(t)を示している。ここでは、壁の外面3と内面4か
ら戻ってきた2つのエコー20.21をはつきbと見分
けると七ができる。信号の2つの主なピークを分離して
いるグラフにおける時間Jt を測定すれば、at=i
μSであることがわかる。さらに、この場合、音波の伝
播速度Cが1. 5 uv’μmとわかっているので、
容易に壁の厚さ●=CIIAt/2を計算することがで
きる。2で割ることは、音波の往復を考慮している。式
のCとJtを1.5カよびlにそれぞれ置き換えれば、
厚さ●の値0.75nが得られる。このように、エコー
かはつきbと分離されている場合には、厚さの測定が簡
単にでき、本発明の方法は必要ない。
第4図は、壁の厚さが0.3朋の場合の信号g(t)を
示している。ここでは、2つのエコーが互いに接近して
おシ、前と同様の操作を行なった場合、ある程度真実に
近い結果を、選択したピークから得ることができる。
示している。ここでは、2つのエコーが互いに接近して
おシ、前と同様の操作を行なった場合、ある程度真実に
近い結果を、選択したピークから得ることができる。
第5図は、壁の厚さが0.3朋の場合の信号g (t)
を示している。ここでは、選択すべき2つのピークを決
定することは困難で、実際には不可能である。エコーは
、一体化して見える。したがって、本発明の目的は、壁
の厚さを測定するため、考慮されるべきピークを分離す
ることができる計算方法を提供することである。この方
法が、第5図に示された初期信号g(t)を生じる管に
適用される。
を示している。ここでは、選択すべき2つのピークを決
定することは困難で、実際には不可能である。エコーは
、一体化して見える。したがって、本発明の目的は、壁
の厚さを測定するため、考慮されるべきピークを分離す
ることができる計算方法を提供することである。この方
法が、第5図に示された初期信号g(t)を生じる管に
適用される。
第5図乃至第13図に示された一連のグラフの他、第2
図に関連して、本発明の方法を説明する。
図に関連して、本発明の方法を説明する。
龜)ブロック11の出力に》ける第5図の複合電気信号
g(t)は、ブロック12において操作される。この操
作は、時間に関する信号から周波数スペクトルに関する
信号に変換する操作で、その振幅Aは、周波数tにより
決筐る。この操作に先だって、信号g (t)にオフセ
ットがある場合には、一般に、オフセット補正が行なわ
れる。極めて一般的には、これは、サーボ機構の研究論
文においてフルワインにより紹介された2に関する変換
に関係している。これに関しては、1947年、ニュー
ヨーク所在のマツクグロウーヒル●ブック社によシ出版
された「サーボ機構の理論」の第5章に記載されている
。我々が興味を示すのは、周知のフーリエ変換である2
変換の特殊なケースを有利に用いることができるという
ことである。ブロック12の出力にkける変換後の信号
G(f)は、第6図に示された実数部Re[G(f)〕
および第7図に示された虚数部ImCG(f)〕2}n
/の形である。
g(t)は、ブロック12において操作される。この操
作は、時間に関する信号から周波数スペクトルに関する
信号に変換する操作で、その振幅Aは、周波数tにより
決筐る。この操作に先だって、信号g (t)にオフセ
ットがある場合には、一般に、オフセット補正が行なわ
れる。極めて一般的には、これは、サーボ機構の研究論
文においてフルワインにより紹介された2に関する変換
に関係している。これに関しては、1947年、ニュー
ヨーク所在のマツクグロウーヒル●ブック社によシ出版
された「サーボ機構の理論」の第5章に記載されている
。我々が興味を示すのは、周知のフーリエ変換である2
変換の特殊なケースを有利に用いることができるという
ことである。ブロック12の出力にkける変換後の信号
G(f)は、第6図に示された実数部Re[G(f)〕
および第7図に示された虚数部ImCG(f)〕2}n
/の形である。
b)その後、周波数スペクトルG(f)の2つの部分、
すなわち実数部と虚数部は、ブロック13(第2図)に
供給される。ここでの目的は、式|G(f) l =
( Rs (G(f)〕2}n/ + Irn (G(
f)〕2}n/”)n/2KL&がって上記スペクトル
のn次の絶対値を計算することである。たとえば、n=
2とすることにより周波数スペクトルの二乗された絶対
値を有利に用いることができる。この二乗された絶対値
は第8図に示されており、ここには、超音波センf(約
10MH!)の中心発信周波数において非常にはつきシ
としたピークがある。
すなわち実数部と虚数部は、ブロック13(第2図)に
供給される。ここでの目的は、式|G(f) l =
( Rs (G(f)〕2}n/ + Irn (G(
f)〕2}n/”)n/2KL&がって上記スペクトル
のn次の絶対値を計算することである。たとえば、n=
2とすることにより周波数スペクトルの二乗された絶対
値を有利に用いることができる。この二乗された絶対値
は第8図に示されており、ここには、超音波センf(約
10MH!)の中心発信周波数において非常にはつきシ
としたピークがある。
C)次に、ブロック13の出力において得られた絶対値
は、ブロック15に供給される。ここでの目的は、実際
の絶対値|G(f)I に関する絶対値の導関数(|
G(f)l ) の商により定義される絶対値の対数
導関数を計算することである。この過程は、本発明の最
も重要な過程で、この過程によって得られる2つのエコ
ーから複合電気信号を分離することができる。第9図に
は、ブロック15の出力にかける対数導関数が示されて
いる。対数導関数は、はつきシした2つの部分を有して
いる。
は、ブロック15に供給される。ここでの目的は、実際
の絶対値|G(f)I に関する絶対値の導関数(|
G(f)l ) の商により定義される絶対値の対数
導関数を計算することである。この過程は、本発明の最
も重要な過程で、この過程によって得られる2つのエコ
ーから複合電気信号を分離することができる。第9図に
は、ブロック15の出力にかける対数導関数が示されて
いる。対数導関数は、はつきシした2つの部分を有して
いる。
すなわち、2つの部分は、超音波の周波数にほぼ集中し
ている有効部分25と、有効部分の両側にあるノイズ部
分26から成っている。関心のある有効部分25は、周
期信号27に1i畳された実質的線形信号28を有して
いる。有効部分のこれら2つの戒分は、次のように物理
的に説明される。
ている有効部分25と、有効部分の両側にあるノイズ部
分26から成っている。関心のある有効部分25は、周
期信号27に1i畳された実質的線形信号28を有して
いる。有効部分のこれら2つの戒分は、次のように物理
的に説明される。
最初の信号g (t)は、実際には、検出されるべき2
つのピークを有している。この2つのピークは、七ンサ
により発信された信号と混合され、ガウス定数によシ振
幅変調された正弦波信号に同化される。信号g (t)
のフーIJ 工変換G(f)f′i.、混合信号を示し
、そこから対数導関数を計算する場合、ピークは周期信
号27に変換され、センサによシ供給される信号は実質
的線形信号28に変換される。
つのピークを有している。この2つのピークは、七ンサ
により発信された信号と混合され、ガウス定数によシ振
幅変調された正弦波信号に同化される。信号g (t)
のフーIJ 工変換G(f)f′i.、混合信号を示し
、そこから対数導関数を計算する場合、ピークは周期信
号27に変換され、センサによシ供給される信号は実質
的線形信号28に変換される。
ここで注目すべきことは、この分離は、過程a)におい
て行なわれた変換にかける時間ドメインに関しては不可
能であるが、周波数ドメインに関しては簡単に得ること
ができる。
て行なわれた変換にかける時間ドメインに関しては不可
能であるが、周波数ドメインに関しては簡単に得ること
ができる。
本発明の望プしいが、不可欠とは言えない例では、その
後、第9図において傾斜の形で示されている線形信号2
8を抑制している。この傾斜の抑制は、後述するような
利点をもたらす。この傾斜すなわち線形信号を抑制する
には、いくつかの方法が使用されてかb1これら方法は
、第lO図に示された信号G!(f)を生じる。
後、第9図において傾斜の形で示されている線形信号2
8を抑制している。この傾斜の抑制は、後述するような
利点をもたらす。この傾斜すなわち線形信号を抑制する
には、いくつかの方法が使用されてかb1これら方法は
、第lO図に示された信号G!(f)を生じる。
第1の方法は、第9図の信号をスクリーン上にビジュア
ル化し、その後、第10図に示したようなグラフ的に満
足のいく結果を得る1で、試行咬たは連続的近似により
この信号を整流して線形信号の抑制を行なう方法である
。これは、傾斜28が周波数の軸と一体化される場合で
ある。このような操作は、第2図のブロック16に示さ
れている。
ル化し、その後、第10図に示したようなグラフ的に満
足のいく結果を得る1で、試行咬たは連続的近似により
この信号を整流して線形信号の抑制を行なう方法である
。これは、傾斜28が周波数の軸と一体化される場合で
ある。このような操作は、第2図のブロック16に示さ
れている。
第2の方法は、数学的方法を用いている。この方法では
、前の過程b)で得られた二乗された絶対値(n=2)
の式から、上記絶対値の1次のモーメントmlと2次の
モーメン}m2を計算する。
、前の過程b)で得られた二乗された絶対値(n=2)
の式から、上記絶対値の1次のモーメントmlと2次の
モーメン}m2を計算する。
これらモーメントは、m1に関してn = 1 >よび
m2に関してはn = 2とすれば次の式により得らl
=一ω 量=+ω 次に、f=mlを通過し、傾斜p=ksm2を有する直
線を決定する。この直線は、関係p(f−ml)によシ
表わされる。fは周波数,kは定数である。
m2に関してはn = 2とすれば次の式により得らl
=一ω 量=+ω 次に、f=mlを通過し、傾斜p=ksm2を有する直
線を決定する。この直線は、関係p(f−ml)によシ
表わされる。fは周波数,kは定数である。
最後に、過程C)において得られた対数導関数から上記
直線を減算して、次の式を得る。
直線を減算して、次の式を得る。
G+(f)=(|G(f)I )’/ |G(f)l
−p(f−n+1 )第2図の流れ図において、ブロ
ック14は、ブロック13で得られた絶対値のモーメン
}mi+m2の計算を表し、ブロック16は、前述した
減算によシ得られた傾斜すなわち線形信号の実際の抑制
を示している。
−p(f−n+1 )第2図の流れ図において、ブロ
ック14は、ブロック13で得られた絶対値のモーメン
}mi+m2の計算を表し、ブロック16は、前述した
減算によシ得られた傾斜すなわち線形信号の実際の抑制
を示している。
第3の方法は、基準測定により先験的に知られている実
質的線形部分の特性を用いている。
質的線形部分の特性を用いている。
d)有効部分25を分離しかつノイズ部分26を除去す
るため、過程C)にシいて得られた対数導関数に対して
周波数の窓を供給する。この操作は、前述したような傾
斜の補正を行なわない場合には、ブロック15に示した
過程C)のすぐ後に行なわれる。傾斜の補正が行なわれ
る(ブロック16)場合には、この補正の前後に関係な
く周波数の窓を供給することができる。傾斜補正が行な
われている第2図の流れ図に釦いては、上記補正(ブロ
ック16)の後に窓が供給(第2のブロック1T)され
ている。第11図のグラフは、窓の供給後の信号を示し
ている。この窓は、第11図に示すように、有効部分の
両側にあるノイズ部分を除去する帯域フィルタである。
るため、過程C)にシいて得られた対数導関数に対して
周波数の窓を供給する。この操作は、前述したような傾
斜の補正を行なわない場合には、ブロック15に示した
過程C)のすぐ後に行なわれる。傾斜の補正が行なわれ
る(ブロック16)場合には、この補正の前後に関係な
く周波数の窓を供給することができる。傾斜補正が行な
われている第2図の流れ図に釦いては、上記補正(ブロ
ック16)の後に窓が供給(第2のブロック1T)され
ている。第11図のグラフは、窓の供給後の信号を示し
ている。この窓は、第11図に示すように、有効部分の
両側にあるノイズ部分を除去する帯域フィルタである。
窓は、ハニングと呼ばれるベル形の窓であることが望1
しい(1958年,アール●ビー−ブラックマンおよび
ジエイ●ダプリュ●タッキーによる、ニューヨークのド
ーバ出版の「パワー●スペクトルの測定」と題する著作
参照)。なか、他の窓を使用してもよく、たとえば、1
981年ダノドにより出版されたエム●クントの「Tr
altement num≦rique dosalg
nauxJと題する著作も考えられる。
しい(1958年,アール●ビー−ブラックマンおよび
ジエイ●ダプリュ●タッキーによる、ニューヨークのド
ーバ出版の「パワー●スペクトルの測定」と題する著作
参照)。なか、他の窓を使用してもよく、たとえば、1
981年ダノドにより出版されたエム●クントの「Tr
altement num≦rique dosalg
nauxJと題する著作も考えられる。
e)窓の供給後の信号は、第2図のブロック18におい
て行なわれる時間ドメインに戻されるならば、使用に便
利である。時間ドメインに戻すのは、過程&)において
説明された変換と逆の変換を信号に対して行なえばよい
。上記説明から明らかなように、この操作は、2に関す
る逆変換により、壕たはここに示されているケースでは
、逆フーリ工変換によりきわめて簡単に行なわれる。変
換後の信号は、第12図に示されている。逆フーリエ変
換を゛rF により、傾斜の補正後の二乗された絶対
値の対数導関数をGK(f)により、および周波数の窓
をH(f)により表すならば、第12図の信号は、次の
式で表される。
て行なわれる時間ドメインに戻されるならば、使用に便
利である。時間ドメインに戻すのは、過程&)において
説明された変換と逆の変換を信号に対して行なえばよい
。上記説明から明らかなように、この操作は、2に関す
る逆変換により、壕たはここに示されているケースでは
、逆フーリ工変換によりきわめて簡単に行なわれる。変
換後の信号は、第12図に示されている。逆フーリエ変
換を゛rF により、傾斜の補正後の二乗された絶対
値の対数導関数をGK(f)により、および周波数の窓
をH(f)により表すならば、第12図の信号は、次の
式で表される。
A(t)=TF [G+(f)・H(f)]f)本発
明の次の過程は、過程e)において行なわれた逆変換の
包絡線を決定する過程である。
明の次の過程は、過程e)において行なわれた逆変換の
包絡線を決定する過程である。
この操作は、第2図のブロック19により行なわれ、そ
の結果は第13図に示されている。第13図に示すよう
に、包絡緩は、様々な振幅の時間的に離間した一連のピ
ーク30.31.32・・・を示している。この結果を
得るには、1984年マツクグロウ●ヒル社によシ出版
された、工−●パボリスによる著作「確率、ランダム変
数、ふ・よび推計プロセス」に釦いて述べられているヒ
ルベルト変換を用いれば有利である。ヒルベルト変換を
THにより表せば、第13図の包絡線は次の式で表され
る。
の結果は第13図に示されている。第13図に示すよう
に、包絡緩は、様々な振幅の時間的に離間した一連のピ
ーク30.31.32・・・を示している。この結果を
得るには、1984年マツクグロウ●ヒル社によシ出版
された、工−●パボリスによる著作「確率、ランダム変
数、ふ・よび推計プロセス」に釦いて述べられているヒ
ルベルト変換を用いれば有利である。ヒルベルト変換を
THにより表せば、第13図の包絡線は次の式で表され
る。
( C A(t)〕2}n/” + C TH −A(
t)〕2}n/” )”″A (t)の値は、前に得ら
れた値である。
t)〕2}n/” )”″A (t)の値は、前に得ら
れた値である。
g)前の過程で得られた包絡線から、最も大きい振幅の
ビーク30から時間原点筐での時間的拡がシltを決定
する。このピークは、前に定義された有効部分の周期信
号によるものである。前述したように、これは、求めら
れているエコーが発見される有効部分の周期信号である
。ここで、時間原点は、最初のエコーに対応する相対的
原点である。測定されたインターバルJt は、2つの
エコーを分離しているインターバルに非常に良く対応し
ている。この原点は、第5図に示した原点とは異なって
いる。これは、過程b)で計算された平方された絶対値
を使用すると、最初のエコーから第5図に示された時間
原点筐での時間を含んでいる情報が失われるからである
。
ビーク30から時間原点筐での時間的拡がシltを決定
する。このピークは、前に定義された有効部分の周期信
号によるものである。前述したように、これは、求めら
れているエコーが発見される有効部分の周期信号である
。ここで、時間原点は、最初のエコーに対応する相対的
原点である。測定されたインターバルJt は、2つの
エコーを分離しているインターバルに非常に良く対応し
ている。この原点は、第5図に示した原点とは異なって
いる。これは、過程b)で計算された平方された絶対値
を使用すると、最初のエコーから第5図に示された時間
原点筐での時間を含んでいる情報が失われるからである
。
特別なケースでは、有効信号の傾斜が抑制されているの
で、周期信号による最も大きい振幅を有するピークは、
最初に出会うビーク30で、管1の壁2の2つの面3.
4の間の距離eは、式e=C−zt/2によシ表される
。第13図のグラフにおいて、Jt=0.2μgと測定
される。c = 1. 5gi Bの場合、e = 1
. 5 ●0.2/ 2 = 0.1 5 ax が得
られ、これは、他の測定方法により得られた管の厚さ0
.15txに良く対応している。
で、周期信号による最も大きい振幅を有するピークは、
最初に出会うビーク30で、管1の壁2の2つの面3.
4の間の距離eは、式e=C−zt/2によシ表される
。第13図のグラフにおいて、Jt=0.2μgと測定
される。c = 1. 5gi Bの場合、e = 1
. 5 ●0.2/ 2 = 0.1 5 ax が得
られ、これは、他の測定方法により得られた管の厚さ0
.15txに良く対応している。
追加考察
本発明の方法を説明するためここで使用されているプラ
スチック管は、特に、測定される管が剛性を有している
ので、はつきシとした信号を生じる。しかし、多くの場
合には、厚さを測定される物体は、剛性金有しておらず
、さらに通常手が届く所にはなく、シかも機械的手段で
は測定できない。たとえば、これは、非侵入手段で壁の
厚さを測定しなければならない生体の動脈の場合である
。
スチック管は、特に、測定される管が剛性を有している
ので、はつきシとした信号を生じる。しかし、多くの場
合には、厚さを測定される物体は、剛性金有しておらず
、さらに通常手が届く所にはなく、シかも機械的手段で
は測定できない。たとえば、これは、非侵入手段で壁の
厚さを測定しなければならない生体の動脈の場合である
。
第14図は、試験管内の(イン・ビトロ)豚のの動脈か
ら、本発明による方法の最後の過程の1つ前の過程にお
いて得られた信号包絡lfsを示している。この場合、
過程C)に釦いて示されたような傾斜補正が行なわれて
いる。包絡線は、いくつかのピークを示し、かつ考察さ
れるべき最大振幅のピークは、0.5μsのピークであ
ることは明白である。
ら、本発明による方法の最後の過程の1つ前の過程にお
いて得られた信号包絡lfsを示している。この場合、
過程C)に釦いて示されたような傾斜補正が行なわれて
いる。包絡線は、いくつかのピークを示し、かつ考察さ
れるべき最大振幅のピークは、0.5μsのピークであ
ることは明白である。
第15図のグラフは、傾斜補正が行なわれていない他は
同じ動脈から得られた包結線を示している。求められて
いるピークに対応するピークは0,56μSにおいて容
易に発見できるが、このピークの前に、これよシ大きい
振幅の他のビーク40がある。このようなピークは、有
効信号の線形部分(すなわち傾斜)により発生され、こ
のピークが、ここでは優勢である。この方法は、測定さ
れる物体がほとんど安定しておらずlた手が届かない場
合などいずれの場合にも、傾斜補正を体系的に使用でき
るようにする。
同じ動脈から得られた包結線を示している。求められて
いるピークに対応するピークは0,56μSにおいて容
易に発見できるが、このピークの前に、これよシ大きい
振幅の他のビーク40がある。このようなピークは、有
効信号の線形部分(すなわち傾斜)により発生され、こ
のピークが、ここでは優勢である。この方法は、測定さ
れる物体がほとんど安定しておらずlた手が届かない場
合などいずれの場合にも、傾斜補正を体系的に使用でき
るようにする。
以上のように、本発明は、動脈の厚さを測定するのに使
用できる方法を提供する。また、この方法は、動脈の内
径を測定するのに使用しても有効的で、しかも薄い層を
有するどんな物体の厚さの測定にも使用できる。
用できる方法を提供する。また、この方法は、動脈の内
径を測定するのに使用しても有効的で、しかも薄い層を
有するどんな物体の厚さの測定にも使用できる。
第1図は、壁の厚さ金測定すべき物体と、コンピュータ
に複合電気信号を供給する七ンサと金示している測定配
置の概要図で、この信号は、本発明による処理過程を受
ける。 第2図は、本発明による訣よび信号の処理の一例による
方法の様々な過程の状態を示した流れ図である。 第3図は、壁の厚さが0.75y+mの場合、センサに
よシ供給される複合信号を示している。ここでは、壁は
、プラスチック材料でできている。 第4図は、0.3msの厚さの場合の信号を示している
。 第5図は、0.15soaの厚さの場合の信号を示して
いる。 第6図は、第5図の信号に関してセンナから生じた信号
を変換することによυ得られた周波数スペクトルの虚数
部を示したグラフである。 @7図は、第5図の信号に関してセンサから生じた信号
金変換することによシ得られた周波数スペクトルの実数
部金示したグラフである。 第8図は、第6図かよび第7図に示した周波数スペクト
ルの二乗された絶対値を示したグラフである。 第9図は、第8図に示した二乗された絶対値の対数導関
数を示したグラフである。 第10図は、曲線調整が行なわれた第9図のグラフと同
様のグラフである。 第11図は、周波数の窓金供給した後の第10図の信号
金示したグラフである。 第12図は、第1l図の悟号の逆変換を表している。 第13図は、拡大された横座標の、第l2図の信号の包
絡線を表したグラフで、との包絡線から管壁の厚さ金計
算することができる。 第14図は、動物の動脈の試験管内測定により得られた
包絡線を示している。 第15図は、第10図の過程が実施さtしていない以外
は、第14図と同じ条件で得られたグラフをボしている
。 1●●●●物体、2●●●●壁、3,4●●●●面、5
・・・・超音波、6●●●・センブ、T11 @ 11
@発振器、8.●●●パーソナル●コンビユータ、9
●●●●ディスプレイ、10●●●●増幅器、e・・・
・距離。
に複合電気信号を供給する七ンサと金示している測定配
置の概要図で、この信号は、本発明による処理過程を受
ける。 第2図は、本発明による訣よび信号の処理の一例による
方法の様々な過程の状態を示した流れ図である。 第3図は、壁の厚さが0.75y+mの場合、センサに
よシ供給される複合信号を示している。ここでは、壁は
、プラスチック材料でできている。 第4図は、0.3msの厚さの場合の信号を示している
。 第5図は、0.15soaの厚さの場合の信号を示して
いる。 第6図は、第5図の信号に関してセンナから生じた信号
を変換することによυ得られた周波数スペクトルの虚数
部を示したグラフである。 @7図は、第5図の信号に関してセンサから生じた信号
金変換することによシ得られた周波数スペクトルの実数
部金示したグラフである。 第8図は、第6図かよび第7図に示した周波数スペクト
ルの二乗された絶対値を示したグラフである。 第9図は、第8図に示した二乗された絶対値の対数導関
数を示したグラフである。 第10図は、曲線調整が行なわれた第9図のグラフと同
様のグラフである。 第11図は、周波数の窓金供給した後の第10図の信号
金示したグラフである。 第12図は、第1l図の悟号の逆変換を表している。 第13図は、拡大された横座標の、第l2図の信号の包
絡線を表したグラフで、との包絡線から管壁の厚さ金計
算することができる。 第14図は、動物の動脈の試験管内測定により得られた
包絡線を示している。 第15図は、第10図の過程が実施さtしていない以外
は、第14図と同じ条件で得られたグラフをボしている
。 1●●●●物体、2●●●●壁、3,4●●●●面、5
・・・・超音波、6●●●・センブ、T11 @ 11
@発振器、8.●●●パーソナル●コンビユータ、9
●●●●ディスプレイ、10●●●●増幅器、e・・・
・距離。
Claims (1)
- 【特許請求の範囲】 物体に向けて速度cで伝播する超音波の発信と、上記物
体の面から上記超音波の反射により発生されたエコーの
、時間の関数として変化する振幅の複合電気信号g(t
)を供給するセンサによる受信とを含む、物体の2つの
面を分離している距離eを測定する方法において、 a)センサからの複合電気信号g(t)を、周波数の関
数である振幅を有しかつ実数部R_e〔G(f)〕およ
び虚数部I_m〔G(f)〕を有する周波数スペクトル
G(f)に変換する過程と、 b)周波数スペクトル|G(f)|^n={R_e〔G
(f)〕+Im〔G(f)〕^2}^n^/^2のn次
の絶対値を計算する過程と、 c)複合電気信号からエコーを分離するため過程b)に
おいて計算された絶対値の対数導関数(|G(f)|^
n)′/|G(f)|^nを計算する過程であつて、上
記対数導関数は、超音波の周波数に実質的に集中した有
効部分と上記有効部分の両側に位置するノイズを有する
部分とから成り、上記有効部分は、実質的線形信号が重
畳された周期信号から成る上記過程と、 d)対数導関数の有効部分を分離しかつノイズを有する
部分を除去するため、対数導関数に周波数の窓を供給す
る過程と、 e)少なくとも上記有効部分の周期信号に過程a)の変
換とは逆の変換を行なう過程と、 f)時間に関して離間した一連の可変振幅ピークを有す
る、逆変換の包絡線を決定する過程と、g)有効部分か
らの周期信号により、最大振幅を有するピークから時間
原点までの時間的拡がりΔtを決定する過程と、 から成る一連の処理過程を上記複合電気信号に施して、
物体の2つの面を分離している距離eはe=c・Δt/
2により表されることを特徴とする物体の2つの面を分
離している距離eを測定する方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR8909926 | 1989-07-20 | ||
FR8909926A FR2650071B1 (fr) | 1989-07-20 | 1989-07-20 | Procede de traitement d'un signal electrique |
Publications (1)
Publication Number | Publication Date |
---|---|
JPH0357907A true JPH0357907A (ja) | 1991-03-13 |
Family
ID=9384066
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2190927A Pending JPH0357907A (ja) | 1989-07-20 | 1990-07-20 | 物体の2つの面を分離している距離eの測定方法 |
Country Status (7)
Country | Link |
---|---|
US (1) | US5197019A (ja) |
EP (1) | EP0409054B1 (ja) |
JP (1) | JPH0357907A (ja) |
CA (1) | CA2021509A1 (ja) |
DE (1) | DE69000081D1 (ja) |
FR (1) | FR2650071B1 (ja) |
NO (1) | NO903643L (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009103459A (ja) * | 2007-10-19 | 2009-05-14 | Ricoh Elemex Corp | 超音波板厚測定装置 |
JP2009257945A (ja) * | 2008-04-17 | 2009-11-05 | Jfe Steel Corp | 薄肉鋼管の肉厚測定方法 |
Families Citing this family (45)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CH685420A5 (fr) * | 1991-09-27 | 1995-07-14 | Asulab Sa | Procédé de mesure par ultrasons de la position d'au moins une paroi mobile et appareil pour la mise en oeuvre de ce procédé. |
US5349859A (en) * | 1991-11-15 | 1994-09-27 | Scientific Engineering Instruments, Inc. | Method and apparatus for measuring acoustic wave velocity using impulse response |
US5404754A (en) * | 1993-02-10 | 1995-04-11 | Shell Oil Company | Ultrasonic detection of high temperature hydrogen attack |
US5723791A (en) * | 1993-09-28 | 1998-03-03 | Defelsko Corporation | High resolution ultrasonic coating thickness gauge |
US5491668A (en) * | 1994-05-13 | 1996-02-13 | Western Atlas International, Inc. | Method for determining the thickness of a casing in a wellbore by signal processing pulse-echo data from an acoustic pulse-echo imaging tool |
DE19612925C1 (de) * | 1996-04-01 | 1997-05-22 | Fraunhofer Ges Forschung | Verfahren zur Werkstoffprüfung von beschichteten Werkstoffen |
US5608165A (en) * | 1996-05-06 | 1997-03-04 | Ford Motor Company | Ultrasonic thickness gauge for multilayer plastic fuel tanks |
ES2352340T3 (es) * | 1997-07-05 | 2011-02-17 | Hudson-Sharp Machine Company | Aparato para la aplicación de cierres resellables sobre una banda de película. |
JP2889568B1 (ja) * | 1998-05-18 | 1999-05-10 | 正男 伊藤 | 血管膜厚測定装置及び動脈硬化診断装置 |
EP1055133A1 (en) * | 1998-12-15 | 2000-11-29 | Koninklijke Philips Electronics N.V. | Ultrasound process and apparatus for the determination of the location of a parietal surface in a tissue and of the absolute radius of an artery |
IT1306874B1 (it) * | 1999-06-25 | 2001-10-11 | Consiglio Nazionale Ricerche | Metodo acustico non invasivo per la rivelazione e la mappatura deidistacchi nei dipinti ad affresco e nei dipinti murali in genere, e |
US6397680B1 (en) * | 2000-07-24 | 2002-06-04 | National Research Council Of Canada | Ultrasonic spectroscopy apparatus for determining thickness and other properties of multilayer structures |
US6494097B1 (en) | 2000-09-05 | 2002-12-17 | Elias Edmond Shihadeh | Method and apparatus for measuring thickness of a layer in a multi-layered object |
JP3643562B2 (ja) * | 2002-02-08 | 2005-04-27 | コーリンメディカルテクノロジー株式会社 | 脈波伝播速度測定装置 |
US8221322B2 (en) * | 2002-06-07 | 2012-07-17 | Verathon Inc. | Systems and methods to improve clarity in ultrasound images |
US20100036252A1 (en) * | 2002-06-07 | 2010-02-11 | Vikram Chalana | Ultrasound system and method for measuring bladder wall thickness and mass |
US20090062644A1 (en) * | 2002-06-07 | 2009-03-05 | Mcmorrow Gerald | System and method for ultrasound harmonic imaging |
US7520857B2 (en) * | 2002-06-07 | 2009-04-21 | Verathon Inc. | 3D ultrasound-based instrument for non-invasive measurement of amniotic fluid volume |
US7819806B2 (en) * | 2002-06-07 | 2010-10-26 | Verathon Inc. | System and method to identify and measure organ wall boundaries |
GB2391625A (en) * | 2002-08-09 | 2004-02-11 | Diagnostic Ultrasound Europ B | Instantaneous ultrasonic echo measurement of bladder urine volume with a limited number of ultrasound beams |
US8221321B2 (en) | 2002-06-07 | 2012-07-17 | Verathon Inc. | Systems and methods for quantification and classification of fluids in human cavities in ultrasound images |
US20060025689A1 (en) * | 2002-06-07 | 2006-02-02 | Vikram Chalana | System and method to measure cardiac ejection fraction |
US20040127797A1 (en) * | 2002-06-07 | 2004-07-01 | Bill Barnard | System and method for measuring bladder wall thickness and presenting a bladder virtual image |
US20080262356A1 (en) * | 2002-06-07 | 2008-10-23 | Vikram Chalana | Systems and methods for ultrasound imaging using an inertial reference unit |
US20090112089A1 (en) * | 2007-10-27 | 2009-04-30 | Bill Barnard | System and method for measuring bladder wall thickness and presenting a bladder virtual image |
US20040100282A1 (en) * | 2002-09-11 | 2004-05-27 | Henrik Christensen | Method and apparatus for determination of layer thickness in a multi-layer structure |
WO2004025221A1 (en) * | 2002-09-11 | 2004-03-25 | Beta Lasermike | Method and apparatus for determination of layer thickness in a multi-layer structure |
US6745632B1 (en) * | 2003-06-03 | 2004-06-08 | Joseph Ernest Dryer | Method for measuring ultrasonic transit times |
US20050038343A1 (en) * | 2003-07-10 | 2005-02-17 | Alfred E. Mann Institute For Biomedical Research At The University Of Southern California | Apparatus and method for locating a bifurcation in an artery |
KR20050014051A (ko) * | 2003-07-29 | 2005-02-07 | 안희태 | 초음파 신호에서 주파수 분리를 이용한 거리측정 방법 및장치 |
CA2598335A1 (en) * | 2004-02-17 | 2005-09-01 | Verathon Inc. | System and method for measuring bladder wall thickness and mass |
FR2886534A1 (fr) * | 2005-06-03 | 2006-12-08 | Theraclion Soc Par Actions Sim | Tete d'imagerie et de traitement d'organes d'etres vivants et procede de fabrication |
FR2886551B1 (fr) * | 2005-06-03 | 2007-09-07 | Theraclion Soc Par Actions Sim | Procede de determination de distance et appareil de traitement mettant en oeuvre ub tel procede |
FR2886533B1 (fr) | 2005-06-03 | 2007-09-14 | Theraclion Soc Par Actions Sim | Tete d'imagerie et de traitement d'organes d'etres vivants et procede de fabrication |
US7434988B1 (en) | 2006-05-17 | 2008-10-14 | Enertechnix, Inc. | Low pressure acoustic pyrometer signal generator |
US7954378B2 (en) * | 2007-02-07 | 2011-06-07 | Pepperl + Fuchs Gmbh | Ultrasonic sensor and method for determining a separation of an object from an ultrasonic sensor |
US8167803B2 (en) * | 2007-05-16 | 2012-05-01 | Verathon Inc. | System and method for bladder detection using harmonic imaging |
NO330292B1 (no) * | 2007-09-12 | 2011-03-21 | Det Norske Veritas As | Akustiske tykkelsesmalinger ved bruk av gass som et koblingsmedium |
US8225998B2 (en) * | 2008-07-11 | 2012-07-24 | Es&S Innovations Llc | Secure ballot box |
WO2010017508A1 (en) * | 2008-08-07 | 2010-02-11 | Verathon Inc. | Device, system, and method to measure abdominal aortic aneurysm diameter |
CN107635472A (zh) | 2015-06-19 | 2018-01-26 | 神经***分析公司 | 经颅多普勒探测器 |
CA2996221C (en) | 2015-08-24 | 2023-09-19 | Hifi Engineering Inc. | Method and system for determining the distance to an acoustically reflective object in a conduit |
JP2019504670A (ja) | 2016-01-05 | 2019-02-21 | ニューラル アナリティクス、インコーポレイテッド | 臨床指標を決定するためのシステム及び方法 |
US10617388B2 (en) | 2016-01-05 | 2020-04-14 | Neural Analytics, Inc. | Integrated probe structure |
US11589836B2 (en) | 2016-01-05 | 2023-02-28 | Novasignal Corp. | Systems and methods for detecting neurological conditions |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2308909A1 (fr) * | 1975-04-22 | 1976-11-19 | Commissariat Energie Atomique | Dispositif de mesure precise des dimensions d'un objet par ultra-sons |
BE840456A (fr) * | 1975-04-22 | 1976-10-07 | Dispositif de mesure precise des dimensions d'un objet par ultra-sons | |
US4063549A (en) * | 1975-12-22 | 1977-12-20 | Technicon Instruments Corporation | Ultrasonic method and apparatus for imaging and characterization of bodies |
US4534220A (en) * | 1984-01-27 | 1985-08-13 | Kleeper John R | Apparatus for measuring the characteristics of an ultrasonic wave medium |
US4648276A (en) * | 1984-01-27 | 1987-03-10 | Klepper John R | Apparatus for measuring the characteristics of an ultrasonic wave medium |
US4817432A (en) * | 1985-09-27 | 1989-04-04 | Design Team Partners | Digital ultrasonic instrument for ophthalmic use |
US4799177A (en) * | 1985-12-31 | 1989-01-17 | The Boeing Company | Ultrasonic instrumentation for examination of variable-thickness objects |
US4855911A (en) * | 1987-11-16 | 1989-08-08 | Massachusetts Institute Of Technology | Ultrasonic tissue characterization |
US4926870A (en) * | 1988-08-30 | 1990-05-22 | Osteo-Technology, Inc. | Method and apparatus for ultrasonic analysis of bone strength in vivo |
FR2637400B1 (fr) * | 1988-09-30 | 1990-11-09 | Labo Electronique Physique | Dispositif de traitement ameliore d'un signal echographique |
FR2644913A1 (fr) * | 1989-03-24 | 1990-09-28 | Labo Electronique Physique | Dispositif d'imagerie par echographie ultrasonore, utilisant un filtre adaptatif perfectionne |
US5038615A (en) * | 1990-05-11 | 1991-08-13 | General Motors Corporation | Ultrasonic multilayer paint thickness measurement |
-
1989
- 1989-07-20 FR FR8909926A patent/FR2650071B1/fr not_active Expired - Lifetime
-
1990
- 1990-07-11 DE DE9090113192T patent/DE69000081D1/de not_active Expired - Lifetime
- 1990-07-11 EP EP90113192A patent/EP0409054B1/fr not_active Expired - Lifetime
- 1990-07-19 US US07/554,420 patent/US5197019A/en not_active Expired - Fee Related
- 1990-07-19 CA CA002021509A patent/CA2021509A1/en not_active Abandoned
- 1990-07-20 JP JP2190927A patent/JPH0357907A/ja active Pending
- 1990-08-17 NO NO90903643A patent/NO903643L/no unknown
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009103459A (ja) * | 2007-10-19 | 2009-05-14 | Ricoh Elemex Corp | 超音波板厚測定装置 |
JP2009257945A (ja) * | 2008-04-17 | 2009-11-05 | Jfe Steel Corp | 薄肉鋼管の肉厚測定方法 |
Also Published As
Publication number | Publication date |
---|---|
US5197019A (en) | 1993-03-23 |
DE69000081D1 (de) | 1992-06-04 |
FR2650071A1 (fr) | 1991-01-25 |
NO903643L (no) | 1991-02-21 |
CA2021509A1 (en) | 1991-01-21 |
EP0409054A1 (fr) | 1991-01-23 |
NO903643D0 (no) | 1990-08-17 |
EP0409054B1 (fr) | 1992-04-29 |
FR2650071B1 (fr) | 1991-09-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JPH0357907A (ja) | 物体の2つの面を分離している距離eの測定方法 | |
US4754760A (en) | Ultrasonic pulse temperature determination method and apparatus | |
US7811233B2 (en) | Finite amplitude distortion-based inhomogeneous pulse echo ultrasonic imaging | |
US5099848A (en) | Method and apparatus for breast imaging and tumor detection using modal vibration analysis | |
US4545250A (en) | Measuring apparatus utilizing spectrum profile | |
Cobb | Finite amplitude method for the determination of the acoustic nonlinearity parameter B/A | |
JPH02280739A (ja) | 超音波エコーグラフ結像装置 | |
US9642600B2 (en) | Shear wave attenuation from k-space analysis system | |
US5103427A (en) | Method and apparatus generating high resolution data and echo identification | |
US4844082A (en) | Ultrasonic examination apparatus | |
Vecchio et al. | Prediction of ultrasonic field propagation through layered media using the extended angular spectrum method | |
JPH03261461A (ja) | 超音波パルスシーケンスの振幅を特徴付ける方法 | |
Sun et al. | Improving time resolution of ultrasonic signals with frequency-domain sparse blind deconvolution (FSBD) method | |
Bowler et al. | An investigation into the feasibility of internal strain measurement in solids by correlation of ultrasonic images | |
Bige et al. | Analysis of microstructural alterations of normal and pathological breast tissue in vivo using the AR cepstrum | |
Pedersen et al. | Ultrasonic monitoring of film condensation for applications in reduced gravity | |
US9918698B2 (en) | System and method for gradient-based k-space search for shear wave velocity dispersion estimation | |
US20020035438A1 (en) | Fast chirp transform | |
Isaev et al. | Measurement of sound reflection coefficients as a function of frequency in an undamped tank | |
CN113367664B (zh) | 一种基于人体浅表皮肤光声显微成像的声学补偿方法及*** | |
US4546772A (en) | Method and means for determining ultrasonic wave attenuation in tissue using phase locked loop | |
JP2669204B2 (ja) | 探査装置 | |
CA2314305A1 (en) | Ultrasonic spectroscopy apparatus and method for determining thickness and other properties of multilayer structures | |
JPH0222660B2 (ja) | ||
Jang et al. | Ultrasound attenuation estimation using the LMSE filters and the median filter |