JPH0420141B2 - - Google Patents
Info
- Publication number
- JPH0420141B2 JPH0420141B2 JP58039907A JP3990783A JPH0420141B2 JP H0420141 B2 JPH0420141 B2 JP H0420141B2 JP 58039907 A JP58039907 A JP 58039907A JP 3990783 A JP3990783 A JP 3990783A JP H0420141 B2 JPH0420141 B2 JP H0420141B2
- Authority
- JP
- Japan
- Prior art keywords
- measurement
- ultrasonic
- continuous
- wave
- nonlinear parameter
- 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
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K11/00—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
- G01K11/22—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using measurement of acoustic effects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H3/00—Measuring characteristics of vibrations by using a detector in a fluid
- G01H3/10—Amplitude; Power
- G01H3/12—Amplitude; Power by electric means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/024—Mixtures
- G01N2291/02491—Materials with nonlinear acoustic properties
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/02881—Temperature
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/04—Wave modes and trajectories
- G01N2291/042—Wave modes
- G01N2291/0421—Longitudinal waves
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Description
【発明の詳細な説明】
(A) 発明の技術分野
本発明は、超音波媒体の等価非線形パラメータ
分布測定装置、特に生体組織等の超音波媒体の物
理特性の空間的分布を測定する方式に係り、音速
が音圧に対し、一次近似としては一定値である
が、二次近似としては音圧に比例するという非線
形性を示すのを利用し、この非線形パラメータの
空間的分布を媒体の特性値として測定し、更には
必要に応じてその空間的分布の映像化を高速に行
うようにした測定装置に関する。
分布測定装置、特に生体組織等の超音波媒体の物
理特性の空間的分布を測定する方式に係り、音速
が音圧に対し、一次近似としては一定値である
が、二次近似としては音圧に比例するという非線
形性を示すのを利用し、この非線形パラメータの
空間的分布を媒体の特性値として測定し、更には
必要に応じてその空間的分布の映像化を高速に行
うようにした測定装置に関する。
(B) 技術の背景と問題点
本発明が利用する非線形パラメータ映像法の原
理は特願昭57−167036号（特許第1561131号（特
公平１−47182号公報））に詳述されている。そこ
では、ポンピング波として波長を順次変化させた
連続波を用い、ポンピング波の各波長に対する受
信信号の位相変化をそれぞれ求め、その波長−位
相関係をフーリエ変換する事により各部位におけ
る（Ｂ／Ａ）eを求めている。しかし、上述の方法に は以下の様な問題がある。第１に、ポンピング波
の周波数をかなり広い範囲にわたつて変化させる
必要があるが、超音波振動子の帯域幅を充分に広
くした上に、ポンピング波の発生に必要とされる
だけの効率の良い（音圧の出せる）振動子を得る
のは必らずしも容易でない。第２に、（ＢＡ）eを 得る際にフーリエ変換を行うが、このため分解能
を上げようとすればする程多くの周波数のポンピ
ング波を用いる必要があり、データの集収に多く
の時間を必要とすると共にフーリエ変換にも時間
を要していた。
理は特願昭57−167036号（特許第1561131号（特
公平１−47182号公報））に詳述されている。そこ
では、ポンピング波として波長を順次変化させた
連続波を用い、ポンピング波の各波長に対する受
信信号の位相変化をそれぞれ求め、その波長−位
相関係をフーリエ変換する事により各部位におけ
る（Ｂ／Ａ）eを求めている。しかし、上述の方法に は以下の様な問題がある。第１に、ポンピング波
の周波数をかなり広い範囲にわたつて変化させる
必要があるが、超音波振動子の帯域幅を充分に広
くした上に、ポンピング波の発生に必要とされる
だけの効率の良い（音圧の出せる）振動子を得る
のは必らずしも容易でない。第２に、（ＢＡ）eを 得る際にフーリエ変換を行うが、このため分解能
を上げようとすればする程多くの周波数のポンピ
ング波を用いる必要があり、データの集収に多く
の時間を必要とすると共にフーリエ変換にも時間
を要していた。
(C) 発明の目的と構成
本発明は、ポンピング波として連続波でなくパ
ルス波を用い、かつポンピング波を、被観察部位
をカバーしてかつ測定用ビームに交叉する平面波
に近いものとするよう構成し、測定用ビームが被
観察部位を透過中にポンピング波によつて受けた
位相変調を復調し、その復調信号に簡単な処理を
施すだけで従来の手法に比べ比較的簡単で容易か
つ高速に被観察部位の等価非線形パラメータ
（Ｂ／Ａ）eの空間分布を得る事のできる超音波媒体の 非線形パラメータ測定装置を提供することを目的
としている。
ルス波を用い、かつポンピング波を、被観察部位
をカバーしてかつ測定用ビームに交叉する平面波
に近いものとするよう構成し、測定用ビームが被
観察部位を透過中にポンピング波によつて受けた
位相変調を復調し、その復調信号に簡単な処理を
施すだけで従来の手法に比べ比較的簡単で容易か
つ高速に被観察部位の等価非線形パラメータ
（Ｂ／Ａ）eの空間分布を得る事のできる超音波媒体の 非線形パラメータ測定装置を提供することを目的
としている。
本発明は、進行する測定波上の一つの点に注目
した場合、その点がポンピング波の影響を受ける
間に、その点が存在した領域の非線形パラメータ
（Ｂ／Ａ）e（但し場所の関数）と、ポンピング波の音 圧Ｐ（但し時間の、従つて場所の関数）とのたた
み込み積分により定まる位相変調を受ける事を利
用して、受信信号を復調して得た位相情報をデコ
ンボリユーシヨンする事によりＢ／Ａの分布を得よ うとするものである。そしてそのため、本発明は
特許請求の範囲記載の構成をもつことを特徴とし
ている。以下具体的に説明する。
した場合、その点がポンピング波の影響を受ける
間に、その点が存在した領域の非線形パラメータ
（Ｂ／Ａ）e（但し場所の関数）と、ポンピング波の音 圧Ｐ（但し時間の、従つて場所の関数）とのたた
み込み積分により定まる位相変調を受ける事を利
用して、受信信号を復調して得た位相情報をデコ
ンボリユーシヨンする事によりＢ／Ａの分布を得よ うとするものである。そしてそのため、本発明は
特許請求の範囲記載の構成をもつことを特徴とし
ている。以下具体的に説明する。
(D) 発明の実施例
超音波媒体内の音圧がゼロの時の音速をC0、
密度をρ0とすると測定波と直交する方向からＰな
る音圧が加えられたときの音速Ｃは Ｃ＝C0＋１／2ρ0C0（Ｂ／Ａ）eＰ ……(1) ここで（Ｂ／Ａ）eは次式で定義され、「等価非線形 パラメータ」と呼ばれるものであり、不均一媒質
においては場所により異つた値を取る。
密度をρ0とすると測定波と直交する方向からＰな
る音圧が加えられたときの音速Ｃは Ｃ＝C0＋１／2ρ0C0（Ｂ／Ａ）eＰ ……(1) ここで（Ｂ／Ａ）eは次式で定義され、「等価非線形 パラメータ」と呼ばれるものであり、不均一媒質
においては場所により異つた値を取る。
（Ｂ／Ａ）e＝2ρ0C0（∂Ci／∂Pj）s ……(2)
ｉとｊは直交する方向であり、Ｓは等エントロ
ピーであることを示す。
ピーであることを示す。
一方、一般に音響の分野で用いられている非線
形パラメータＢ／Ａは次のように定義されている。
形パラメータＢ／Ａは次のように定義されている。
Ｂ／Ａ＝2ρ0C0（∂Ci／∂Pi）s……(3)
等方体であればボアソン比をνとして（Ｂ／Ａ）e
とＢ／Ａとの間には
（Ｂ／Ａ）e＝（１＋2ν−１／１−ν）Ｂ／Ａ…
…(4) という関係があり、液体、生体組織ではν≒0.5
であるので（Ｂ／Ａ）eＢ／Ａとみなすことができる
。
…(4) という関係があり、液体、生体組織ではν≒0.5
であるので（Ｂ／Ａ）eＢ／Ａとみなすことができる
。
従つて、ポンピング波の音圧Ｐにより、音速Ｃ
は、 △Ｃ＝１／2ρ0C0（Ｂ／Ａ）e・Ｐ ……(5) だけ変化する事になる。
は、 △Ｃ＝１／2ρ0C0（Ｂ／Ａ）e・Ｐ ……(5) だけ変化する事になる。
今、第１図に示す様に、送信用超音波振動子
XAから測定用の連続波超音波ビームWnを図のＺ
軸方向に被測定媒体中に送り込み、受信用超音波
振動子XBでこのビームを受信する。この時、測
定用ビームWnと交叉する方向から、別の超音波
振動子XPにより平面パルス波WPをポンピング波
として図示の様にｘ軸方向に送り込む。以下簡単
のため交叉角は直角であるとする。また、Wnの
音圧はWPの音圧に比べて充分に小さく、式(1)〜
(5)のＰとしてはWPによるものだけを考慮すれば
良いものとする。
XAから測定用の連続波超音波ビームWnを図のＺ
軸方向に被測定媒体中に送り込み、受信用超音波
振動子XBでこのビームを受信する。この時、測
定用ビームWnと交叉する方向から、別の超音波
振動子XPにより平面パルス波WPをポンピング波
として図示の様にｘ軸方向に送り込む。以下簡単
のため交叉角は直角であるとする。また、Wnの
音圧はWPの音圧に比べて充分に小さく、式(1)〜
(5)のＰとしてはWPによるものだけを考慮すれば
良いものとする。
ポンピング波WPが測定用ビームWnと交叉して
いる間には、測定用ビームWnの各Ｚにおける波
形は、式(5)により、場所により異なる音速変化 △Ｃ（ｚ）＝１／2ρ0C0（Ｂ／Ａ）e（ｚ）・Ｐ……(
6) を受ける事になり、従つて場所（ｚ）により異つ
た△Ｃ（ｚ）に比例する位相変化を受ける事にな
る。従つて、受信振動子XBからの受信信号の位
相情報及び予め測定しておいたＰを用いて逆に
（Ｂ／Ａ）e（ｚ）を推定する事ができる。
いる間には、測定用ビームWnの各Ｚにおける波
形は、式(5)により、場所により異なる音速変化 △Ｃ（ｚ）＝１／2ρ0C0（Ｂ／Ａ）e（ｚ）・Ｐ……(
6) を受ける事になり、従つて場所（ｚ）により異つ
た△Ｃ（ｚ）に比例する位相変化を受ける事にな
る。従つて、受信振動子XBからの受信信号の位
相情報及び予め測定しておいたＰを用いて逆に
（Ｂ／Ａ）e（ｚ）を推定する事ができる。
以下では（Ｂ／Ａ）e（ｚ）の推定法を詳述する。
第２図は、測定用ビームWnが進行するにつれ
て、ポンピング波WPからどの様な音圧を受ける
かという事を模式的に示したものである。
て、ポンピング波WPからどの様な音圧を受ける
かという事を模式的に示したものである。
ポンピング波の音圧をＰ（ｔ）とする。時刻ｔ
＝０でｚ＝z0にあつて音圧Ｐ（ｏ）を受けていた
測定波（第２図Ａ参照）は、時刻ｔ＝△ｔにはｚ
＝z0＋△ｚ＝z0＋C0△ｔの位置にあり、この時に
はＰ（ｔ＋△ｔ）の音圧を受ける（第２図Ｂ参
照）。以下同様にして時刻ｔ＝２△ｔにおいて第
２図Ｃ図示、時刻ｔ＝３△ｔにおいて第２図Ｄ図
示の如き音圧を受ける事になる。これから判る様
に、一般に、ｔ＝０においてｚ＝z0にあつた測定
波が、ｚ＝ｚで受ける音速変化は △Ｃ（ｚ）＝１／2ρ0C0（Ｂ／Ａ）e（ｚ）Ｐ（△ｔ） ＝１／2ρ0C0（Ｂ／Ａ）e（ｚ）Ｐ（ｚ−z0／C0）…
…(7) 但し、ここで△Ｃ（ｚ）の値はC0に比べて小さ
く、△ｔは測定波がｚ＝z0からｚ＝ｚまで移動す
るのに要する時間で △ｔ＝ｚ−z0／C0 ……(8) で与えられると近似した。
＝０でｚ＝z0にあつて音圧Ｐ（ｏ）を受けていた
測定波（第２図Ａ参照）は、時刻ｔ＝△ｔにはｚ
＝z0＋△ｚ＝z0＋C0△ｔの位置にあり、この時に
はＰ（ｔ＋△ｔ）の音圧を受ける（第２図Ｂ参
照）。以下同様にして時刻ｔ＝２△ｔにおいて第
２図Ｃ図示、時刻ｔ＝３△ｔにおいて第２図Ｄ図
示の如き音圧を受ける事になる。これから判る様
に、一般に、ｔ＝０においてｚ＝z0にあつた測定
波が、ｚ＝ｚで受ける音速変化は △Ｃ（ｚ）＝１／2ρ0C0（Ｂ／Ａ）e（ｚ）Ｐ（△ｔ） ＝１／2ρ0C0（Ｂ／Ａ）e（ｚ）Ｐ（ｚ−z0／C0）…
…(7) 但し、ここで△Ｃ（ｚ）の値はC0に比べて小さ
く、△ｔは測定波がｚ＝z0からｚ＝ｚまで移動す
るのに要する時間で △ｔ＝ｚ−z0／C0 ……(8) で与えられると近似した。
従つて、ｔ＝０でｚ＝z0にあつた測定波がｚ＝
ｚで受ける位相の変化はＫを比例定数として、 △（ｚ）＝Ｋ・（Ｂ／Ａ）e（ｚ）・Ｐ（ｚ−z0／C0
） ……(9) となり、この測定波が送信振動子XAを出てから
受信振動子XBに受信されるまでの間に受ける位
相変化の総和は φ（z0）＝∫∞ -∞△（ｚ）dz ＝K∫∞ -∞（Ｂ／Ａ）e（ｚ）・Ｐ（ｚ−z0／C0）dz
……(10) となる。ここで ｇ（z0−ｚ）＝Ｐ（ｚ−z0／Ｃ）……（11） と置くと式(10)は φ（z0）＝K∫∞ -∞（Ｂ／Ａ）e（ｚ）ｇ（z0−ｚ）dy ＝Ｋ〔（Ｂ／Ａ）e（ｚ）＊ｇ（ｚ）〕 ……（12） となる。
ｚで受ける位相の変化はＫを比例定数として、 △（ｚ）＝Ｋ・（Ｂ／Ａ）e（ｚ）・Ｐ（ｚ−z0／C0
） ……(9) となり、この測定波が送信振動子XAを出てから
受信振動子XBに受信されるまでの間に受ける位
相変化の総和は φ（z0）＝∫∞ -∞△（ｚ）dz ＝K∫∞ -∞（Ｂ／Ａ）e（ｚ）・Ｐ（ｚ−z0／C0）dz
……(10) となる。ここで ｇ（z0−ｚ）＝Ｐ（ｚ−z0／Ｃ）……（11） と置くと式(10)は φ（z0）＝K∫∞ -∞（Ｂ／Ａ）e（ｚ）ｇ（z0−ｚ）dy ＝Ｋ〔（Ｂ／Ａ）e（ｚ）＊ｇ（ｚ）〕 ……（12） となる。
ここで、式（12）の「＊」はいわゆるたたみ込
み積分である。この式は、（Ｂ／Ａ）e（ｚ）に対し て、ｇ（ｚ）なる関数が重み関数となつてφ（z0）
を与える事を示している。
み積分である。この式は、（Ｂ／Ａ）e（ｚ）に対し て、ｇ（ｚ）なる関数が重み関数となつてφ（z0）
を与える事を示している。
式（12）をフーリエ変換すれば
Φ（ω）＝Ｋ・（Ｂ／Ａ）e（ω）・Ｇ（ω）……（13
） （但しΦ（ω）、（Ｂ／Ａ）e（ω）、Ｇ（ω）はそれ
ぞ れφ（ｚ）、（Ｂ／Ａ）e（ｚ）、ｇ（ｚ）のフーリエ
変 換、ωは座標軸Ｚに対する空間周波数） となり、これから （Ｂ／Ａ）e（ω）＝１／Ｋ１／Ｇ（ω）Φ（ω）……
（14） （Ｂ／Ａ）e（ｚ）＝１／ＫF-1〔１／Ｇ（ω）Φ（ω）
〕 ……（15） として（Ｂ／Ａ）e（ｚ）を求める事ができる。
） （但しΦ（ω）、（Ｂ／Ａ）e（ω）、Ｇ（ω）はそれ
ぞ れφ（ｚ）、（Ｂ／Ａ）e（ｚ）、ｇ（ｚ）のフーリエ
変 換、ωは座標軸Ｚに対する空間周波数） となり、これから （Ｂ／Ａ）e（ω）＝１／Ｋ１／Ｇ（ω）Φ（ω）……
（14） （Ｂ／Ａ）e（ｚ）＝１／ＫF-1〔１／Ｇ（ω）Φ（ω）
〕 ……（15） として（Ｂ／Ａ）e（ｚ）を求める事ができる。
式（14）は、Φ（ω）を、１／Ｋ１／Ｇ（ω）なる周
波 数特性のフイルタに通す事によつて（Ｂ／Ａ）e（ω） が得られる事を示している。従つて（Ｂ／Ａ）e（ｚ） を求める際には必ずしも式（15）の如きフーリエ
逆変換を行う必要はなく、周波数特性が１／Ｋ １／Ｇ（ω）であるフイルタを予め用意しておき、そ れにφ（ｚ）を入力する事により（Ｂ／Ａ）e（ｚ）が 出力として得られる。またＧ（ω）は式（11）を
用いてＰ（ｔ）から求めておく事ができる。この
様にして得られたＢ／Ａ（ｚ）を映像として表示す れば、Ｂ／Ａ（ｚ）の空間分布像を得る事もできる。
波 数特性のフイルタに通す事によつて（Ｂ／Ａ）e（ω） が得られる事を示している。従つて（Ｂ／Ａ）e（ｚ） を求める際には必ずしも式（15）の如きフーリエ
逆変換を行う必要はなく、周波数特性が１／Ｋ １／Ｇ（ω）であるフイルタを予め用意しておき、そ れにφ（ｚ）を入力する事により（Ｂ／Ａ）e（ｚ）が 出力として得られる。またＧ（ω）は式（11）を
用いてＰ（ｔ）から求めておく事ができる。この
様にして得られたＢ／Ａ（ｚ）を映像として表示す れば、Ｂ／Ａ（ｚ）の空間分布像を得る事もできる。
以上、理論的な説明を行なつたが、次に本発明
の１実施例について第３図を用いて説明する。
の１実施例について第３図を用いて説明する。
第３図において、１はポンピング波の送信のタ
イミングを発生するタイミング制御部、２は測定
用連続波のための発振器、３は発振器出力を受け
て振動子を駆動するためのドライバ、４は測定波
送信用の振動子XA、５は被測定超音波媒体Ｔ、
６は測定波受信用の振動子XB、７は受信増幅器、
８は位相検出器、９は式（14）の１／Ｋ１／Ｇ（ω）で
定 義されるフイルタ、１０はポンピング波用のドラ
イバ、１１はポンピング波発生用の振動子XPで
ある。またx1はXPから測定用ビームまでの距離、
z2は被測定区間の距離、z1は被測定区間とXBとの
距離である。第４図には第３図の主要部の時間波
形を示しており、第４図には第３図と同じ信号名
称を記してある。ポンピング波を送信してからポ
ンピング波が測定波の位置に到達するまでの時間
（x1／C0）と、ポンピング波による位相変調を受け た測定波が初めて振動子XBに到達するのに要す
る時間（z1／C0）との和をt1とすると、受信増幅器 ７の出力VRは、第４図に示す様にポンピング波
の送信から時間t1の後に位相変調を受け始め、ポ
ンピング波による位相変調を受けた測定波が被測
定区間を通過してしまうまでの時間t2＝ｚ／C0の間 位相変調出力を出し続ける。位相検出器８は、発
振器２の出力とVRの位相を比較し、VRの位相を
時間の関数φ（ｔ）として、従つて座標ｚの関数
として出力する。この出力が、フイルタ９に入力
され、等価非線形パラメータ（Ｂ／Ａ）eが時間の関 数Ｂ／Ａ（ｔ）として、従つて座標ｚの関数として 得られる。
イミングを発生するタイミング制御部、２は測定
用連続波のための発振器、３は発振器出力を受け
て振動子を駆動するためのドライバ、４は測定波
送信用の振動子XA、５は被測定超音波媒体Ｔ、
６は測定波受信用の振動子XB、７は受信増幅器、
８は位相検出器、９は式（14）の１／Ｋ１／Ｇ（ω）で
定 義されるフイルタ、１０はポンピング波用のドラ
イバ、１１はポンピング波発生用の振動子XPで
ある。またx1はXPから測定用ビームまでの距離、
z2は被測定区間の距離、z1は被測定区間とXBとの
距離である。第４図には第３図の主要部の時間波
形を示しており、第４図には第３図と同じ信号名
称を記してある。ポンピング波を送信してからポ
ンピング波が測定波の位置に到達するまでの時間
（x1／C0）と、ポンピング波による位相変調を受け た測定波が初めて振動子XBに到達するのに要す
る時間（z1／C0）との和をt1とすると、受信増幅器 ７の出力VRは、第４図に示す様にポンピング波
の送信から時間t1の後に位相変調を受け始め、ポ
ンピング波による位相変調を受けた測定波が被測
定区間を通過してしまうまでの時間t2＝ｚ／C0の間 位相変調出力を出し続ける。位相検出器８は、発
振器２の出力とVRの位相を比較し、VRの位相を
時間の関数φ（ｔ）として、従つて座標ｚの関数
として出力する。この出力が、フイルタ９に入力
され、等価非線形パラメータ（Ｂ／Ａ）eが時間の関 数Ｂ／Ａ（ｔ）として、従つて座標ｚの関数として 得られる。
以上の如く、本発明によれば、超音波媒体の等
価非線形パラメータの空間分布（Ｂ／Ａ）e（ｚ）を、 従来の様にフーリエ変換の様な操作を伴う事な
く、簡単な構成で高速に得る事ができる。
価非線形パラメータの空間分布（Ｂ／Ａ）e（ｚ）を、 従来の様にフーリエ変換の様な操作を伴う事な
く、簡単な構成で高速に得る事ができる。
以下、ポンピング波として正しい平面波が得ら
れない場合の対策について述べる。一般に、完全
な平面波を得るには無限に広い平板振動子を用い
る必要があり、これは実際上不可能である。近似
的に平面波を得る方法としては、振動子面から充
分に離れた所で、充分に曲率半径の大きくなつた
球面波の一部を使う事もあるが、ポンピング波と
して用いるには、単位面積当りのエネルギー密度
が小さくなり過ぎるので必らずしも適当でない。
従つて、充分な音圧を持つたポンピング波は、実
際には平面波とはならず、従つてその音圧は、時
間ｔだけでなく座標ｚによつても変化する事にな
り、平面波からのズレが無視できない位大きくな
つた場合にはその補正が必要となる。この補正は
以下の様にして行なう事ができる。
れない場合の対策について述べる。一般に、完全
な平面波を得るには無限に広い平板振動子を用い
る必要があり、これは実際上不可能である。近似
的に平面波を得る方法としては、振動子面から充
分に離れた所で、充分に曲率半径の大きくなつた
球面波の一部を使う事もあるが、ポンピング波と
して用いるには、単位面積当りのエネルギー密度
が小さくなり過ぎるので必らずしも適当でない。
従つて、充分な音圧を持つたポンピング波は、実
際には平面波とはならず、従つてその音圧は、時
間ｔだけでなく座標ｚによつても変化する事にな
り、平面波からのズレが無視できない位大きくな
つた場合にはその補正が必要となる。この補正は
以下の様にして行なう事ができる。
先ず、ポンピング波の音圧分布を、時間ｔ及び
座標の関数Ｐ（ｚ、ｔ）として測定波ビームの線
上で予め測定しておく。このＰ（ｚ、ｔ）から、
ｔ＝０でｚ＝z0にあつた測定波が、測定波ビーム
線上で受ける音圧変化をz0をパラメータとしてｚ
の関数gz0（ｚ）として求めておく事ができる。こ
のgz0（ｚ）が式(9)のＰ（ｚ−z0／C0）に相当するもの となる。これを用いて、ｔ＝０でｚ＝z0にあつた
測定波に対する重み関数（式（11）のｇ（ｚ−z0）
に相当する）を平面波の場合と同様にして求めて
おく事ができる。従つて、この重み関数の逆特性
を持つフイルタを平面波の場合と同様に予め用意
しておく事ができ、従つて（Ｂ／Ａ）e（z0）を得る 事ができる。但し、平面波の場合と異なり、測定
波の注目している部分がｔ＝０でｚ軸上のどの部
分にあるかによつて、つまりz0の値によつて、そ
の注目している部分がｚ軸上で受ける音圧分布が
異つてくる。その様子を第５図に模式的に示す。
なお、第５図Ａはｔ＝０の時点に球面波が到来し
ている場合の音圧分布を表わし、第５図Ｂはｔ＝
０においてｚ＝z1にあつた測定波に対して本来の
平面波のポンピング波を受けた場合のｚ軸上での
音圧分布、第５図Ｃはｔ＝０においてｚ＝z1にあ
つた測定波が球面波のポンピング波を受けた場合
のｚ軸上での音圧分布、第５図Ｄはｔ＝０におい
てｚ＝０であつた測定波が球面波のポンピング波
を受けた場合のｚ軸上での音圧分布を夫々表わし
ている。従つて、式（11）のｇ（ｚ）に相当する
重み関数も、更には式（14）の１／Ｇ（ω）に対応す るフイルタも、z0をパラメータとして異なつたも
のとなる。つまり、このフイルタは、どのｚ座標
のＢ／Ａ（ｚ）を求めるのかという事により変化す るものとなるが、第５図に示す如く、z0の各値に
対して注目している測定波がｚ軸上で受ける音圧
を求める事ができるので、これから各z0の値に対
するフイルタ特性を予め計算しておく事ができ
る。第４図における位相検出器８の出力を時間的
に見るとこれは、逐次異なるz0に対応する位相変
化が出力されているから、この出力に対し、各z0
に対する上述のフイルタを順次作用させれば、平
面ポンピング波の時と同様にＢ／Ａ（ｚ）の分布を 得る事ができる。この時、順次異なる特性のフイ
ルタを作用させる方法としては、例えば、位相検
出器を一旦Ａ／Ｄ変換した上でデイジタル・フイ
ルタリングを施しても良い。
座標の関数Ｐ（ｚ、ｔ）として測定波ビームの線
上で予め測定しておく。このＰ（ｚ、ｔ）から、
ｔ＝０でｚ＝z0にあつた測定波が、測定波ビーム
線上で受ける音圧変化をz0をパラメータとしてｚ
の関数gz0（ｚ）として求めておく事ができる。こ
のgz0（ｚ）が式(9)のＰ（ｚ−z0／C0）に相当するもの となる。これを用いて、ｔ＝０でｚ＝z0にあつた
測定波に対する重み関数（式（11）のｇ（ｚ−z0）
に相当する）を平面波の場合と同様にして求めて
おく事ができる。従つて、この重み関数の逆特性
を持つフイルタを平面波の場合と同様に予め用意
しておく事ができ、従つて（Ｂ／Ａ）e（z0）を得る 事ができる。但し、平面波の場合と異なり、測定
波の注目している部分がｔ＝０でｚ軸上のどの部
分にあるかによつて、つまりz0の値によつて、そ
の注目している部分がｚ軸上で受ける音圧分布が
異つてくる。その様子を第５図に模式的に示す。
なお、第５図Ａはｔ＝０の時点に球面波が到来し
ている場合の音圧分布を表わし、第５図Ｂはｔ＝
０においてｚ＝z1にあつた測定波に対して本来の
平面波のポンピング波を受けた場合のｚ軸上での
音圧分布、第５図Ｃはｔ＝０においてｚ＝z1にあ
つた測定波が球面波のポンピング波を受けた場合
のｚ軸上での音圧分布、第５図Ｄはｔ＝０におい
てｚ＝０であつた測定波が球面波のポンピング波
を受けた場合のｚ軸上での音圧分布を夫々表わし
ている。従つて、式（11）のｇ（ｚ）に相当する
重み関数も、更には式（14）の１／Ｇ（ω）に対応す るフイルタも、z0をパラメータとして異なつたも
のとなる。つまり、このフイルタは、どのｚ座標
のＢ／Ａ（ｚ）を求めるのかという事により変化す るものとなるが、第５図に示す如く、z0の各値に
対して注目している測定波がｚ軸上で受ける音圧
を求める事ができるので、これから各z0の値に対
するフイルタ特性を予め計算しておく事ができ
る。第４図における位相検出器８の出力を時間的
に見るとこれは、逐次異なるz0に対応する位相変
化が出力されているから、この出力に対し、各z0
に対する上述のフイルタを順次作用させれば、平
面ポンピング波の時と同様にＢ／Ａ（ｚ）の分布を 得る事ができる。この時、順次異なる特性のフイ
ルタを作用させる方法としては、例えば、位相検
出器を一旦Ａ／Ｄ変換した上でデイジタル・フイ
ルタリングを施しても良い。
特別な場合として、測定波ビームがポンピング
波の遠距離音場にある場合を考えると、ポンピン
グ波が球面波であるので、測定波ビーム上での音
圧分布はＰ（ωt−ｋ√2＋2＋2） と書くことができる。ここでｙ≫ｘ、ｙであり、
ポンプ波の十分遠距離音場であるとすると、ビー
ム上（ｙ＝y0）で Ｐ（ωt−ｋ√2＋2＋2） ≒Ｐ（ωt−ｋ（y0＋x2＋y2／2y0）） ＝Ｐ（（ωt−ky0）−ｋx2＋y2／2y0） ……（16） 第１項は（ｘ、ｙ、ｚ）＝（０、y0、０）なる映
像面の中心の位相であり、第２項は位置による位
相のずれを表わしている。従つて今 Ｚ＝ｚ＋ｋ（x2＋y2／2y0） ……（17） という変数変換を行えば、式(9)は △（Ｚ）＝Ｋ（Ｂ／Ａ）e（ｚ）Ｐ（Ｚ−Z0／C0） ……（18） となり、以下式(9)ないし（15）によつて、（Ｂ／Ａ）e （ｚ）の像が得られる。よつてその像を式（17）
の関係から変数をＺからｚに戻してやると（Ｂ／Ａ）e （ｚ）の像が得られる。このようにポンピング
波と測定波ビームとをある程度距離をもたせるこ
とにより信号処理が非常に容易になる。以上の如
くポンピング波は必ずしも測定波に直交する平面
波でなくともよく、任意の形状、任意の交差角の
ものが利用できる。
波の遠距離音場にある場合を考えると、ポンピン
グ波が球面波であるので、測定波ビーム上での音
圧分布はＰ（ωt−ｋ√2＋2＋2） と書くことができる。ここでｙ≫ｘ、ｙであり、
ポンプ波の十分遠距離音場であるとすると、ビー
ム上（ｙ＝y0）で Ｐ（ωt−ｋ√2＋2＋2） ≒Ｐ（ωt−ｋ（y0＋x2＋y2／2y0）） ＝Ｐ（（ωt−ky0）−ｋx2＋y2／2y0） ……（16） 第１項は（ｘ、ｙ、ｚ）＝（０、y0、０）なる映
像面の中心の位相であり、第２項は位置による位
相のずれを表わしている。従つて今 Ｚ＝ｚ＋ｋ（x2＋y2／2y0） ……（17） という変数変換を行えば、式(9)は △（Ｚ）＝Ｋ（Ｂ／Ａ）e（ｚ）Ｐ（Ｚ−Z0／C0） ……（18） となり、以下式(9)ないし（15）によつて、（Ｂ／Ａ）e （ｚ）の像が得られる。よつてその像を式（17）
の関係から変数をＺからｚに戻してやると（Ｂ／Ａ）e （ｚ）の像が得られる。このようにポンピング
波と測定波ビームとをある程度距離をもたせるこ
とにより信号処理が非常に容易になる。以上の如
くポンピング波は必ずしも測定波に直交する平面
波でなくともよく、任意の形状、任意の交差角の
ものが利用できる。
本発明は更に、測定した（Ｂ／Ａ）e（ｚ）分布の
Ｓ／Ｎ比を向上させる手段を提供している。等価
非線形パラメータと音圧との積による測定波の位
相変化量が小さい時には、本来の位相変化量φ
（ｚ）に対して回路内その他で発生する雑音を無
視できなくなる。このため、φ（ｚ）から求めた
（Ｂ／Ａ）e（ｚ）にも大きな雑音が含まれる事にな る。この様子を第６図に示す。第６図Ａ図示の無
雑音の位相検出器出力φ（ｚ）に第６図Ｂ図示の
雑音Ｎが加わる事により、第６図Ｃの図示の如き
信号が実際には出力される。この対策として、同
一測定部位をＫ回測定して第６図Ｃ図示のS1〜
SKの如き雑音に加わつた信号を得て、これらを
同一ｚ座標の点毎に加算すると、各点において信
号成分は振幅でＫ倍されるが、雑音成分は電力で
Ｋ倍されるに過ぎず、もし雑音が不規則雑音であ
れば各点における雑音振幅は√倍されるにとど
まる。従つてＳ／Ｎ比は√倍改善されて第６図
Ｄ図示の如き出力が得られる事になる。
非線形パラメータと音圧との積による測定波の位
相変化量が小さい時には、本来の位相変化量φ
（ｚ）に対して回路内その他で発生する雑音を無
視できなくなる。このため、φ（ｚ）から求めた
（Ｂ／Ａ）e（ｚ）にも大きな雑音が含まれる事にな る。この様子を第６図に示す。第６図Ａ図示の無
雑音の位相検出器出力φ（ｚ）に第６図Ｂ図示の
雑音Ｎが加わる事により、第６図Ｃの図示の如き
信号が実際には出力される。この対策として、同
一測定部位をＫ回測定して第６図Ｃ図示のS1〜
SKの如き雑音に加わつた信号を得て、これらを
同一ｚ座標の点毎に加算すると、各点において信
号成分は振幅でＫ倍されるが、雑音成分は電力で
Ｋ倍されるに過ぎず、もし雑音が不規則雑音であ
れば各点における雑音振幅は√倍されるにとど
まる。従つてＳ／Ｎ比は√倍改善されて第６図
Ｄ図示の如き出力が得られる事になる。
この方法を用いたシステム構成例を第７図に示
す。第７図において、第３図と同一の構成要素に
は同一の番号を付しており説明は省略する。第７
図において、１３はポンピング波の送信繰り返し
周期Ｔだけ信号を遅らせる遅延回路であり、例え
ばCCDやBBD等のアナログ的手段で実現しても
良いし、またＡ／Ｄ変換器とシフトレジスタ又は
メモリとを用いたデイジタル的手段で実現して良
い事は言うまでもない。１２は加算器であり、遅
延回路の種類に応じてアナログ又はデイジタルの
いずれのタイプでも良い。第７図図示の構成の場
合、位相検出器８の出力Siが同一のｚ軸座標の各
点毎に加算される事は明らかであり、いわゆる同
期加算によりＳ／Ｎ比の改善を行なう事ができ
る。
す。第７図において、第３図と同一の構成要素に
は同一の番号を付しており説明は省略する。第７
図において、１３はポンピング波の送信繰り返し
周期Ｔだけ信号を遅らせる遅延回路であり、例え
ばCCDやBBD等のアナログ的手段で実現しても
良いし、またＡ／Ｄ変換器とシフトレジスタ又は
メモリとを用いたデイジタル的手段で実現して良
い事は言うまでもない。１２は加算器であり、遅
延回路の種類に応じてアナログ又はデイジタルの
いずれのタイプでも良い。第７図図示の構成の場
合、位相検出器８の出力Siが同一のｚ軸座標の各
点毎に加算される事は明らかであり、いわゆる同
期加算によりＳ／Ｎ比の改善を行なう事ができ
る。
本発明においては、更に（Ｂ／Ａ）e（ｚ）の２次
元又は３次元分布を得る事が可能となる。これま
での説明で明らかな様に、特定の走査線上の
（Ｂ／Ａ）e（ｚ）の分布を得る事ができるから、送 信・受信振動子の対を相対位置を一定に保つたま
までｘまたは／およびｙ方向に移動させ、各ｘま
たは／およびｙ座標に対応したメモリアドレスに
（Ｂ／Ａ）e（ｚ）の値を記憶しておけば、（Ｂ／Ａ）e （ｚ）の２次元または／および３次元分布を得る
事ができる。
での説明で明らかな様に、特定の走査線上の
（Ｂ／Ａ）e（ｚ）の分布を得る事ができるから、送 信・受信振動子の対を相対位置を一定に保つたま
までｘまたは／およびｙ方向に移動させ、各ｘま
たは／およびｙ座標に対応したメモリアドレスに
（Ｂ／Ａ）e（ｚ）の値を記憶しておけば、（Ｂ／Ａ）e （ｚ）の２次元または／および３次元分布を得る
事ができる。
第８図はその構成例で、第７図と同じ構成要素
には同一番号を付しており、説明は省略する。第
８図において、駆動部１４は例えばステツピング
モータによつて振動子４，６の対を移動させる部
分であり、タイミング制御部１からの制御パルス
に応じて振動子の対の移動を行なう。位置検出部
１５は例えばステツピングモータの軸に取り付け
られたロータリ・エンコーダ等により振動子４，
６の対の位置を検出する部分である。メモリ制御
部１６は、タイミング制御部１からのポンピング
波送信同期信号やクロツク等を受けて位置検出部
１５からの出力に対応したメモリアドレスを発生
する。メモリ１７は、メモリ制御部１６からの書
込／読出等の制御信号及びアドレスに従つて、フ
イルタ９の出力である（Ｂ／Ａ）e（ｚ）をＡ／Ｄ変 換した後、記憶する。表示部１８は記憶された
（Ｂ／Ａ）e（ｚ）の分布をメモリから読み出して表示 する。
には同一番号を付しており、説明は省略する。第
８図において、駆動部１４は例えばステツピング
モータによつて振動子４，６の対を移動させる部
分であり、タイミング制御部１からの制御パルス
に応じて振動子の対の移動を行なう。位置検出部
１５は例えばステツピングモータの軸に取り付け
られたロータリ・エンコーダ等により振動子４，
６の対の位置を検出する部分である。メモリ制御
部１６は、タイミング制御部１からのポンピング
波送信同期信号やクロツク等を受けて位置検出部
１５からの出力に対応したメモリアドレスを発生
する。メモリ１７は、メモリ制御部１６からの書
込／読出等の制御信号及びアドレスに従つて、フ
イルタ９の出力である（Ｂ／Ａ）e（ｚ）をＡ／Ｄ変 換した後、記憶する。表示部１８は記憶された
（Ｂ／Ａ）e（ｚ）の分布をメモリから読み出して表示 する。
以上の如くすれば（Ｂ／Ａ）e（ｚ）の分布の２次
元又は３次元分布を得る事ができる事は明らかで
ある。
ある。
尚、第８図において、振動子４，６の対の移動
は機械的に行うものとしたが、振動子の対として
いわゆる電子スキヤン・プローブの対を用いる事
により、機械的走査を行なわなくても（Ｂ／Ａ）e （ｚ）の２次元分布が得られる様になる事は言う
までもない。又、第８図において、メモリ１７を
用いず、表示管の残光特性を用いて（Ｂ／Ａ）e（ｚ） の２次元分布像を得ても良い事は言うまでもな
い。
は機械的に行うものとしたが、振動子の対として
いわゆる電子スキヤン・プローブの対を用いる事
により、機械的走査を行なわなくても（Ｂ／Ａ）e （ｚ）の２次元分布が得られる様になる事は言う
までもない。又、第８図において、メモリ１７を
用いず、表示管の残光特性を用いて（Ｂ／Ａ）e（ｚ） の２次元分布像を得ても良い事は言うまでもな
い。
(E) 発明の効果
以上述べた如く、本発明によれば、以前に提案
した超音波媒体の非線形パラメータ（Ｂ／Ａ）e（ｚ） の測定に比べ、より簡単なハードウエアで、より
高速にその分布を測定する事ができる。
した超音波媒体の非線形パラメータ（Ｂ／Ａ）e（ｚ） の測定に比べ、より簡単なハードウエアで、より
高速にその分布を測定する事ができる。
第１図および第２図は本発明の概念を説明する
説明図、第３図は本発明の一実施例構成、第４図
は第３図図示構成の動作を説明する説明図、第５
図は球面波のポンピング波が到来する場合を説明
する説明図、第６図はＳ／Ｎ比を向上する態様を
説明する説明図、第７図は第６図に対応する一実
施例構成、第８図は二次元または三次元像を得る
場合の一実施例構成を示す。 図中、１はタイミング制御部、２は発振器、３
はドライバ、４は送信用振動子、５は被測定超音
波媒体、６は受信用振動子、７は受信増幅器、８
は位相検出器、９はフイルタ、１０はポンピング
波用のドライバ、１１はポンピング波発生用振動
子、１３は遅延回路、１４は駆動部を表わす。
説明図、第３図は本発明の一実施例構成、第４図
は第３図図示構成の動作を説明する説明図、第５
図は球面波のポンピング波が到来する場合を説明
する説明図、第６図はＳ／Ｎ比を向上する態様を
説明する説明図、第７図は第６図に対応する一実
施例構成、第８図は二次元または三次元像を得る
場合の一実施例構成を示す。 図中、１はタイミング制御部、２は発振器、３
はドライバ、４は送信用振動子、５は被測定超音
波媒体、６は受信用振動子、７は受信増幅器、８
は位相検出器、９はフイルタ、１０はポンピング
波用のドライバ、１１はポンピング波発生用振動
子、１３は遅延回路、１４は駆動部を表わす。
Claims (1)
- 【特許請求の範囲】 １ (イ) 超音波媒体を透過する測定用連続超音波
ビームを送信・受信する一対の超音波振動子
と、 (ロ) 該測定用連続超音波ビームと交差し、被観察
部位をカバーするだけのビーム幅を持つポンピ
ング・パルス用超音波を送信する超音波振動子
と、 (ハ) 該ポンピング・パルスによつて、各々の部位
における等価非線形パラメータ値に応じた位相
遷移を受けて透過してくる測定用連続超音波を
受信してその位相遷移量を時間の関数として検
出する手段と、 (ニ) 該時間関数としての位相遷移量を、上記ポン
ピング・パルスの周波数特性の逆特性にほぼ等
しいフイルタに通すことにより、上記測定用連
続超音波ビーム上の等価非線形パラメータ超音
波の空間分布を得る手段と、 を持つことを特徴とする超音波媒体の等価非線形
パラメータ分布測定装置。 ２ (イ) 超音波媒体を透過する測定用連続超音波
ビームを送信・受信する一対の超音波振動子
と、 (ロ) 該測定用連続超音波ビームと交差し、被観察
部位をカバーするだけのビーム幅を持つポンピ
ング・パルス用超音波を送信する超音波振動子
と、 (ハ) 該ポンピング・パルスによつて、各々の部位
における等価非線形パラメータ値に応じた位相
遷移を受けて透過してくる測定用連続超音波を
受信してその位相遷移量を時間の関数として検
出する手段と、 (ニ) 該時間関数としての位相遷移量を、上記ポン
ピング・パルスの周波数特性の逆特性にほぼ等
しいフイルタに通すことにより、上記測定用連
続超音波ビーム上の等価非線形パラメータ超音
波の空間分布を得る手段と、 (ホ) ポンピング波の音圧分布を場所及び時間の関
数として予め測定しておき、上記等価非線形パ
ラメータの測定線上の空間分布を得るに当たつ
て、上記予め測定された音圧分布を用いて補正
を行う手段と、 を持つことを特徴とする超音波媒体の等価非線形
パラメータ分布測定装置。 ３ (イ) 超音波媒体を透過する測定用連続超音波
ビームを送信・受信する一対の超音波振動子
と、 (ロ) 該測定用連続超音波ビームと交差し、被観察
部位をカバーするだけのビーム幅を持つポンピ
ング・パルス用超音波を送信する超音波振動子
と、 (ハ) 該ポンピング・パルスによつて、各々の部位
における等価非線形パラメータ値に応じた位相
遷移を受けて透過してくる測定用連続超音波を
受信してその位相遷移量を時間の関数として検
出する手段と、 (ニ) 該時間関数としての位相遷移量を、上記ポン
ピング・パルスの周波数特性の逆特性にほぼ等
しいフイルタに通すことにより、上記測定用連
続超音波ビーム上の等価非線形パラメータ超音
波の空間分布を得る手段と、 (ホ) 上記等価非線形パラメータの測定ビーム線上
の空間分布を得るに当たつて、同一の測定ビー
ム線上の同一位置に対応する位相遷移量を、複
数回のポンピング波の送信に対応して加算する
同期加算回路と、 を持つことを特徴とする超音波媒体の等価非線形
パラメータ分布測定装置。 ４ (イ) 超音波媒体を透過する測定用連続超音波
ビームを送信・受信する一対の超音波振動子
と、 (ロ) 該測定用連続超音波ビームと交差し、被観察
部位をカバーするだけのビーム幅を持つポンピ
ング・パルス用超音波を送信する超音波振動子
と、 (ハ) 該ポンピング・パルスによつて、各々の部位
における等価非線形パラメータ値に応じた位相
遷移を受けて透過してくる測定用連続超音波を
受信してその位相遷移量を時間の関数として検
出する手段と、 (ニ) 該時間関数としての位相遷移量を、上記ポン
ピング・パルスの周波数特性の逆特性にほぼ等
しいフイルタに通すことにより、上記測定用連
続超音波ビーム上の等価非線形パラメータ超音
波の空間分布を得る手段と、 (ホ) 上記送信・受信振動子の対を一次元的または
二次元的に走査する事により、二次元的または
三次元的な等価非線形パラメータ分布像を得る
ための手段と、 を持つことを特徴とする超音波媒体の等価非線形
パラメータ分布測定装置。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP58039907A JPS59164956A (ja) | 1983-03-10 | 1983-03-10 | 超音波媒体の等価非線形パラメ−タ分布測定装置 |
| US06/587,522 US4566460A (en) | 1983-03-10 | 1984-03-08 | Measuring method and apparatus for non-linear parameter of acoustic medium and its application |
| EP84301591A EP0128635B1 (en) | 1983-03-10 | 1984-03-09 | Measuring non-linear parameter of acoustic medium |
| DE8484301591T DE3475651D1 (en) | 1983-03-10 | 1984-03-09 | Measuring non-linear parameter of acoustic medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP58039907A JPS59164956A (ja) | 1983-03-10 | 1983-03-10 | 超音波媒体の等価非線形パラメ−タ分布測定装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS59164956A JPS59164956A (ja) | 1984-09-18 |
| JPH0420141B2 true JPH0420141B2 (ja) | 1992-03-31 |
Family
ID=12566021
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP58039907A Granted JPS59164956A (ja) | 1983-03-10 | 1983-03-10 | 超音波媒体の等価非線形パラメ−タ分布測定装置 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US4566460A (ja) |
| EP (1) | EP0128635B1 (ja) |
| JP (1) | JPS59164956A (ja) |
| DE (1) | DE3475651D1 (ja) |
Families Citing this family (28)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6060838A (ja) * | 1983-09-13 | 1985-04-08 | 富士通株式会社 | 超音波音場観測装置 |
| US4610255A (en) * | 1983-12-02 | 1986-09-09 | Fujitsu Limited | Ultrasonic non-linear parameter measuring system |
| JPS6111025A (ja) * | 1984-06-26 | 1986-01-18 | 株式会社東芝 | 超音波組織診断装置 |
| JPS61100237A (ja) * | 1984-10-23 | 1986-05-19 | 株式会社東芝 | 超音波診断装置 |
| US4802490A (en) * | 1984-11-01 | 1989-02-07 | Johnston G Gilbert | Catheter for performing volumetric flow rate determination in intravascular conduits |
| FR2577324B1 (fr) * | 1985-02-08 | 1987-03-06 | Labo Electronique Physique | Appareil d'exploration de milieux par echographie ultrasonore |
| US4807633A (en) * | 1986-05-21 | 1989-02-28 | Indianapolis Center For Advanced Research | Non-invasive tissue thermometry system and method |
| US4970656A (en) * | 1986-11-07 | 1990-11-13 | Alcon Laboratories, Inc. | Analog drive for ultrasonic probe with tunable phase angle |
| JPS63122923A (ja) * | 1986-11-13 | 1988-05-26 | Agency Of Ind Science & Technol | 超音波測温装置 |
| US5001649A (en) * | 1987-04-06 | 1991-03-19 | Alcon Laboratories, Inc. | Linear power control for ultrasonic probe with tuned reactance |
| JPH01299537A (ja) * | 1988-05-27 | 1989-12-04 | Agency Of Ind Science & Technol | 音響特性測定装置及び測温装置 |
| DE3829999A1 (de) | 1988-09-01 | 1990-03-15 | Schering Ag | Ultraschallverfahren und schaltungen zu deren durchfuehrung |
| US5143070A (en) * | 1989-11-17 | 1992-09-01 | The University Of Texas Systems Board Of Regents | Transaxial compression technique for sound velocity estimation |
| US5293870A (en) * | 1989-11-17 | 1994-03-15 | Board Of Regents The University Of Texas System | Method and apparatus for elastographic measurement and imaging |
| US5108369A (en) * | 1990-03-15 | 1992-04-28 | Diagnostic Devices Group, Limited | Dual-diameter multifunction catheter |
| US5086775A (en) * | 1990-11-02 | 1992-02-11 | University Of Rochester | Method and apparatus for using Doppler modulation parameters for estimation of vibration amplitude |
| DE4229817C2 (de) * | 1992-09-07 | 1996-09-12 | Siemens Ag | Verfahren zur zerstörungsfreien und/oder nichtinvasiven Messung einer Temperaturänderung im Inneren eines insbesondere lebenden Objektes |
| JP3094742B2 (ja) * | 1993-09-03 | 2000-10-03 | 松下電器産業株式会社 | 超音波診断装置 |
| US5482045A (en) * | 1994-10-12 | 1996-01-09 | Advanced Technology Laboratories, Inc. | Ultrasonic diagnostic system gain control |
| AU2833097A (en) * | 1996-05-10 | 1997-12-05 | Christian H. Constantinov | System and method for registration and reproduction of sensory fields |
| US5921928A (en) * | 1996-12-05 | 1999-07-13 | Mayo Foundation For Medical Education And Research | Acoustic force generation by amplitude modulating a sonic beam |
| KR100439656B1 (ko) * | 2002-01-04 | 2004-07-12 | (주)휴먼정보통신 | 초음파를 이용한 비접촉식 온도분포 계측장치 |
| US20050277835A1 (en) * | 2003-05-30 | 2005-12-15 | Angelsen Bjorn A | Ultrasound imaging by nonlinear low frequency manipulation of high frequency scattering and propagation properties |
| US8038616B2 (en) * | 2003-05-30 | 2011-10-18 | Surf Technology As | Acoustic imaging by nonlinear low frequency manipulation of high frequency scattering and propagation properties |
| US6984209B2 (en) * | 2003-07-02 | 2006-01-10 | The Brigham And Women's Hospital, Inc. | Harmonic motion imaging |
| JP5399192B2 (ja) * | 2009-09-30 | 2014-01-29 | 富士フイルム株式会社 | 超音波診断装置、および超音波診断装置の作動方法 |
| US20110092818A1 (en) * | 2009-10-16 | 2011-04-21 | Artann Laboratories, Inc. | Ultrasonometer for bone assessment in infants |
| CN115372952B (zh) * | 2022-07-08 | 2025-08-22 | 陕西师范大学 | 一种涡旋波束多基地水下目标检测装置与方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4431008A (en) * | 1982-06-24 | 1984-02-14 | Wanner James F | Ultrasonic measurement system using a perturbing field, multiple sense beams and receivers |
| JPS5955245A (ja) * | 1982-09-25 | 1984-03-30 | 富士通株式会社 | 超音波診断処理方式 |
-
1983
- 1983-03-10 JP JP58039907A patent/JPS59164956A/ja active Granted
-
1984
- 1984-03-08 US US06/587,522 patent/US4566460A/en not_active Expired - Lifetime
- 1984-03-09 EP EP84301591A patent/EP0128635B1/en not_active Expired
- 1984-03-09 DE DE8484301591T patent/DE3475651D1/de not_active Expired
Also Published As
| Publication number | Publication date |
|---|---|
| EP0128635A2 (en) | 1984-12-19 |
| EP0128635A3 (en) | 1986-01-15 |
| EP0128635B1 (en) | 1988-12-14 |
| DE3475651D1 (en) | 1989-01-19 |
| JPS59164956A (ja) | 1984-09-18 |
| US4566460A (en) | 1986-01-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH0420141B2 (ja) | ||
| US10222461B2 (en) | Phase inversion ultrasonic imaging | |
| US7785260B2 (en) | Ultrasonograph | |
| JP3510025B2 (ja) | 超音波撮像装置 | |
| JPH05237117A (ja) | コヒーレント振動エネルギビーム作像装置及び伝達媒質内の収差を補正する方法 | |
| CN1744858B (zh) | 超声波摄像装置 | |
| JP3093823B2 (ja) | 超音波ドプラ診断装置 | |
| WO2003032841A1 (en) | Apparatus and method for indicating mechanical stiffness properties of body tissue | |
| JP7235722B2 (ja) | 超音波診断装置および超音波診断装置の制御方法 | |
| US7481769B2 (en) | Ultrasonic diagnosing apparatus | |
| JPH03155843A (ja) | 超音波診断装置 | |
| JPH0556466B2 (ja) | ||
| JP3177677B2 (ja) | 超音波連続波ドプラ血流計 | |
| CN102240214B (zh) | 超声波诊断装置 | |
| JP3519111B2 (ja) | 超音波診断装置 | |
| JP2002034985A (ja) | 超音波診断装置 | |
| JP3492095B2 (ja) | 超音波撮像装置 | |
| JPH056142B2 (ja) | ||
| JPH0139057B2 (ja) | ||
| JPS60144660A (ja) | 超音波エコーグラフイによる物体の検査装置 | |
| JPH08173431A (ja) | 超音波診断装置 | |
| JPH08117227A (ja) | 超音波診断装置 | |
| JPH08257027A (ja) | 超音波診断装置 | |
| JPS63194644A (ja) | 超音波計測装置 | |
| JPS6279045A (ja) | 合成開口法を用いた超音波診断装置 |