JPH0331059B2 - - Google Patents
Info
- Publication number
- JPH0331059B2 JPH0331059B2 JP59109831A JP10983184A JPH0331059B2 JP H0331059 B2 JPH0331059 B2 JP H0331059B2 JP 59109831 A JP59109831 A JP 59109831A JP 10983184 A JP10983184 A JP 10983184A JP H0331059 B2 JPH0331059 B2 JP H0331059B2
- Authority
- JP
- Japan
- Prior art keywords
- ultrasonic
- pumping
- pulse
- sound pressure
- measurement
- 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
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Description
【発明の詳細な説明】
(A) 発明の技術分野
本発明は超音波媒体の非線形パラメータ分布測
定装置、特に、生体組織等の超音波媒体の物理特
性の空間分布を測定する方式に係り、音速が音圧
に対し、一次近似としては一定値であるが、二次
近似としては音圧に比例するという非線形性を示
すのを利用し、この非線形性パラメータの空間的
分布を媒体の特性値として測定し、更には、必要
に応じてその空間的分布の映像化を高速且つ容易
に行う様にした測定装置に関する。更には、特
に、いわゆる透過法でなく反射法によつて上記非
線形パラメータを測定する事ができる装置に関す
る。
定装置、特に、生体組織等の超音波媒体の物理特
性の空間分布を測定する方式に係り、音速が音圧
に対し、一次近似としては一定値であるが、二次
近似としては音圧に比例するという非線形性を示
すのを利用し、この非線形性パラメータの空間的
分布を媒体の特性値として測定し、更には、必要
に応じてその空間的分布の映像化を高速且つ容易
に行う様にした測定装置に関する。更には、特
に、いわゆる透過法でなく反射法によつて上記非
線形パラメータを測定する事ができる装置に関す
る。
(B) 技術の背景と問題点
本発明が利用する非線形パラメータ映像法の原
理は、特願昭57−167036号又は特願昭58−119100
号に詳述されている。前者においては、比較的高
周波の測定用連続超音波ビームに直交する方向か
ら、比較的低周波のポンピング用平面パルス波を
交差させ、ポンピングパルスによつて位相変調さ
れた測定波を位相復調する事によつて、測定波ビ
ーム走査線上の非線形パラメータＢ／Ａを高速に
求めていた。又、後者においては、ポンピングパ
ルスを、測定用ビームに直交する様に与える代わ
りに、測定用ビームと同軸状で且つ進行方向が測
定用ビームと逆方向になる様に与える事により、
測定用ビームが殆どあらゆる所でほぼ同一形状の
ポンピングパルスを受ける様にして、機械的に大
きな構造をなくす等の効果を得ていた。
理は、特願昭57−167036号又は特願昭58−119100
号に詳述されている。前者においては、比較的高
周波の測定用連続超音波ビームに直交する方向か
ら、比較的低周波のポンピング用平面パルス波を
交差させ、ポンピングパルスによつて位相変調さ
れた測定波を位相復調する事によつて、測定波ビ
ーム走査線上の非線形パラメータＢ／Ａを高速に
求めていた。又、後者においては、ポンピングパ
ルスを、測定用ビームに直交する様に与える代わ
りに、測定用ビームと同軸状で且つ進行方向が測
定用ビームと逆方向になる様に与える事により、
測定用ビームが殆どあらゆる所でほぼ同一形状の
ポンピングパルスを受ける様にして、機械的に大
きな構造をなくす等の効果を得ていた。
しかし、これらの方法はいずれも、測定波送信
用振動子から放射された後、被測定媒体内を通過
して測定波受信用振動子で受信された信号を利用
するという、いわゆる透過法を用いているため、
実際の臨床応用において断層画像を得ようとする
と、その適用部位が測定用振動子対の間に被観察
部位を設ける事ができる様な乳房等に限られる
か、あるいは被観察部位を水枕で囲つて測定用振
動子対を走査する必要がある等、操作上の問題が
大きかつた。
用振動子から放射された後、被測定媒体内を通過
して測定波受信用振動子で受信された信号を利用
するという、いわゆる透過法を用いているため、
実際の臨床応用において断層画像を得ようとする
と、その適用部位が測定用振動子対の間に被観察
部位を設ける事ができる様な乳房等に限られる
か、あるいは被観察部位を水枕で囲つて測定用振
動子対を走査する必要がある等、操作上の問題が
大きかつた。
(C) 発明の目的と構成
本発明は、測定波として連続波を用いる代わり
にバースト状のパルス波を送受信兼用の測定波用
振動子から送信すると共に、該測定用超音波パル
スと殆ど同一場所から殆ど同一方向にポンピング
用の比較的低周波の超音波パルスを被観察媒体内
に送出する様にしておき、且つ、第１図に示す様
に、ポンピング波の音圧の高い部分（又は低い部
分）に測定用のバースト波が重畳する様に振動子
の駆動タイミングを調整しておき、ポンピングパ
ルスと測定用パルスの両方を送信した時に反射さ
れて帰つて来た測定用パルスの受信信号の位相
と、測定用パルスだけを送信した時に反射されて
帰つてきた測定用パルスの受信信号の位相との差
を求める事により、ポンピングパルスの影響だけ
による測定用パルスの位相変調を反射法で検出し
て、被観察媒体内の超音波非線形パラメータＢ／
Ａを求めようとするものである。
にバースト状のパルス波を送受信兼用の測定波用
振動子から送信すると共に、該測定用超音波パル
スと殆ど同一場所から殆ど同一方向にポンピング
用の比較的低周波の超音波パルスを被観察媒体内
に送出する様にしておき、且つ、第１図に示す様
に、ポンピング波の音圧の高い部分（又は低い部
分）に測定用のバースト波が重畳する様に振動子
の駆動タイミングを調整しておき、ポンピングパ
ルスと測定用パルスの両方を送信した時に反射さ
れて帰つて来た測定用パルスの受信信号の位相
と、測定用パルスだけを送信した時に反射されて
帰つてきた測定用パルスの受信信号の位相との差
を求める事により、ポンピングパルスの影響だけ
による測定用パルスの位相変調を反射法で検出し
て、被観察媒体内の超音波非線形パラメータＢ／
Ａを求めようとするものである。
本発明は、進行する測定波上の１つの点に注目
した場合、その点が反射体に至るまでの間に、測
定用パルスが通過した領域の非線形パラメータ
Ｂ／Ａ（但し場所の関数）とポンピング波の音圧
Ｐ（減衰の影響により場所の関数）との積の積分
値により定まる位相変調を受ける事を利用して、
受信信号を復調して得た位相信号を微分する事に
より、Ｂ／Ａの分布を得ようとするものである。
そしてそのため、本発明は特許請求の範囲記載の
構成を持つ事を特徴としている。以下具体的に説
明する。
した場合、その点が反射体に至るまでの間に、測
定用パルスが通過した領域の非線形パラメータ
Ｂ／Ａ（但し場所の関数）とポンピング波の音圧
Ｐ（減衰の影響により場所の関数）との積の積分
値により定まる位相変調を受ける事を利用して、
受信信号を復調して得た位相信号を微分する事に
より、Ｂ／Ａの分布を得ようとするものである。
そしてそのため、本発明は特許請求の範囲記載の
構成を持つ事を特徴としている。以下具体的に説
明する。
(D) 発明の実施例
超音波媒体内の音圧がゼロの時の音速をC0、
密度をρ0とすると、Ｐなる音圧が加えられた時の
音速Ｃは Ｃ＝C0＋１／2ρ0C0・Ｂ／Ａ・Ｐ …(1) （但しＢ／Ａ＝2ρ0C0（aC／aP）s （但し添字Ｓは等エントロピーである事を示
す。））となる。従つて、ポンピングパルスの音圧
Ｐにより、音速Ｃは、 ΔC＝１／2ρ0C0 Ｂ／ＡＰ …(2) だけ変化する事になる。
密度をρ0とすると、Ｐなる音圧が加えられた時の
音速Ｃは Ｃ＝C0＋１／2ρ0C0・Ｂ／Ａ・Ｐ …(1) （但しＢ／Ａ＝2ρ0C0（aC／aP）s （但し添字Ｓは等エントロピーである事を示
す。））となる。従つて、ポンピングパルスの音圧
Ｐにより、音速Ｃは、 ΔC＝１／2ρ0C0 Ｂ／ＡＰ …(2) だけ変化する事になる。
今、第２図に示す様に、Ｚ＝０に置かれた測定
波用超音波振動子Xnからバースト状の超音波パ
ルスWnを図のＺ軸方向に被測定媒体中に送り込
み、反射体Ｍからの反射波を受信する場合を考え
る。この時、測定パルスWnと殆ど同一場所に設
置されたポンピングパルス送信用振動子Xpから、
ポンピングパルスWpを図示の様にＺ軸方向に送
り込む。以下簡単のためにWpとWnとは重畳して
いると考えることとする。又、音圧Wnは音圧Wp
に比べて充分に小さく、式(1)(2)のＰとしてはWp
によるものだけを考慮すれば良いものとする。測
定用パルスWnがＺ＝Ｚにある反射体Ｍに到達す
るまでの間に、各Ｚにおいて、式(2)により、場所
により異なる音速変化 ΔC（Ｚ）＝１／2ρ0C0・Ｂ／Ａ（Ｚ）・Ｐ（Ｚ）…(
3) を受ける事になり、従つて場所（Ｚ）により異な
つたΔC（Ｚ）に比例する位相変化を受ける事にな
る。従つて、反射体Ｍに到達した時には、第３式
で示される音速変化の積分値に比例した位相変化 Φ（Ｚ）＝K∫z 0ΔC（Ｚ）dZ（Ｋ：比例定数） ＝K∫z 0１／2ρ0C0・Ｂ／Ａ（Ｚ）・Ｐ（Ｚ）dZ…(4) を受けている事になる。反射体Ｍで反射された
後、振動子Xnで受信されるまでの間は、ポンピ
ングパルスが存在しない（付確には、ポンピング
パルスの反射波が、測定用パルスの反射波と共に
戻つて来るが、往路での音圧に比べれば充分に小
さくなつている）ので、振動子Xnで受信された
測定用パルスの反射波を位相復調すれば、第４式
で示されるΦ（Ｚ）が得られる事になる。
波用超音波振動子Xnからバースト状の超音波パ
ルスWnを図のＺ軸方向に被測定媒体中に送り込
み、反射体Ｍからの反射波を受信する場合を考え
る。この時、測定パルスWnと殆ど同一場所に設
置されたポンピングパルス送信用振動子Xpから、
ポンピングパルスWpを図示の様にＺ軸方向に送
り込む。以下簡単のためにWpとWnとは重畳して
いると考えることとする。又、音圧Wnは音圧Wp
に比べて充分に小さく、式(1)(2)のＰとしてはWp
によるものだけを考慮すれば良いものとする。測
定用パルスWnがＺ＝Ｚにある反射体Ｍに到達す
るまでの間に、各Ｚにおいて、式(2)により、場所
により異なる音速変化 ΔC（Ｚ）＝１／2ρ0C0・Ｂ／Ａ（Ｚ）・Ｐ（Ｚ）…(
3) を受ける事になり、従つて場所（Ｚ）により異な
つたΔC（Ｚ）に比例する位相変化を受ける事にな
る。従つて、反射体Ｍに到達した時には、第３式
で示される音速変化の積分値に比例した位相変化 Φ（Ｚ）＝K∫z 0ΔC（Ｚ）dZ（Ｋ：比例定数） ＝K∫z 0１／2ρ0C0・Ｂ／Ａ（Ｚ）・Ｐ（Ｚ）dZ…(4) を受けている事になる。反射体Ｍで反射された
後、振動子Xnで受信されるまでの間は、ポンピ
ングパルスが存在しない（付確には、ポンピング
パルスの反射波が、測定用パルスの反射波と共に
戻つて来るが、往路での音圧に比べれば充分に小
さくなつている）ので、振動子Xnで受信された
測定用パルスの反射波を位相復調すれば、第４式
で示されるΦ（Ｚ）が得られる事になる。
Φ（Ｚ）をＺで微分する事により
dΦ（Ｚ）／dZ＝Ｋ１／2ρ0C0・Ｂ／Ａ（Ｚ）・Ｐ（Ｚ
）…(5) が得られ、第５式の左辺は実測で得られる値、右
辺のＫ、ρ0、C0は定数であるから、Ｐ（Ｚ）を知
る事ができればＢ／Ａ（Ｚ）を求める事ができる。
）…(5) が得られ、第５式の左辺は実測で得られる値、右
辺のＫ、ρ0、C0は定数であるから、Ｐ（Ｚ）を知
る事ができればＢ／Ａ（Ｚ）を求める事ができる。
超音波周波数の例として、ポンピングパルスWp
としては生体組織内での減衰が余り大きくない
500KHz程度のものが用いられ、測定波パルスWn
としてはそれよりも１桁程度高い5MHz程度のも
のが用いられる。生体組織中では、超音波はほぼ
1dB／ＭHz／cmの減衰を受ける事が良く知られて
おり、500KHzのWpの場合、ほぼ0.5dB／cmの減
衰を受ける事になる。従つて、第５式のＰ（Ｚ）
として、例えば という式を用いてＰ（Ｚ）を推定しても良い。こ
の場合、第５式より、 但しK′＝2ρ0C0／KP（０） α＝0.5／20loge10 つまり、位相復調出力Φ（Ｚ）を微分したものに、
距離Ｚと共に増大する様な係数K′e+〓Zを乗ずる事
によりＢ／Ａ（Ｚ）を得る事ができる。
としては生体組織内での減衰が余り大きくない
500KHz程度のものが用いられ、測定波パルスWn
としてはそれよりも１桁程度高い5MHz程度のも
のが用いられる。生体組織中では、超音波はほぼ
1dB／ＭHz／cmの減衰を受ける事が良く知られて
おり、500KHzのWpの場合、ほぼ0.5dB／cmの減
衰を受ける事になる。従つて、第５式のＰ（Ｚ）
として、例えば という式を用いてＰ（Ｚ）を推定しても良い。こ
の場合、第５式より、 但しK′＝2ρ0C0／KP（０） α＝0.5／20loge10 つまり、位相復調出力Φ（Ｚ）を微分したものに、
距離Ｚと共に増大する様な係数K′e+〓Zを乗ずる事
によりＢ／Ａ（Ｚ）を得る事ができる。
生体組織中での減衰が１・dB／ＭHz／cmとい
う様に場所によらず一定、という様な仮定が成立
しない場合は、通常のＢモード断層撮像装置でよ
く用いられている様な、距離Ｚ毎に第３図の如き
自由なゲインを与える事のできるいわゆるTGC
（Time Gain Control）を用いても良い事は言う
までもない。
う様に場所によらず一定、という様な仮定が成立
しない場合は、通常のＢモード断層撮像装置でよ
く用いられている様な、距離Ｚ毎に第３図の如き
自由なゲインを与える事のできるいわゆるTGC
（Time Gain Control）を用いても良い事は言う
までもない。
尚、媒体中の平均音速をC0とした時、Ｚ≒Cot
（ｔ：超音波が送信されてからの時間）であるか
ら、dZ＝Co dt ∴dΦ（Ｚ）／dZ＝１／C0 dΦ（Ｚ）／dt …(8) となり、第５式のＺによる微分は、時間微分に置
き換えられるのは言うまでもない。
（ｔ：超音波が送信されてからの時間）であるか
ら、dZ＝Co dt ∴dΦ（Ｚ）／dZ＝１／C0 dΦ（Ｚ）／dt …(8) となり、第５式のＺによる微分は、時間微分に置
き換えられるのは言うまでもない。
以上の説明では、距離Ｚにある反射体からの反
射波の位相Φ（Ｚ）が求まると仮定していたが、
以下の理由で、これは容易ではない。
射波の位相Φ（Ｚ）が求まると仮定していたが、
以下の理由で、これは容易ではない。
第１に、位相を検出するためには第４図ａに示
す如き基準信号と受信信号との位相比較を行うよ
うにされている。しかし反射法の場合、第４図
ｂ，ｃに示す如く、振動子と反射体との間の距離
が異なるだけで、基準信号に対する受信信号の位
相が大きく変化し、Ｂ／Ａ（Ｚ）・Ｐ（Ｚ）の影響
による位相シフトを隠してしまう。
す如き基準信号と受信信号との位相比較を行うよ
うにされている。しかし反射法の場合、第４図
ｂ，ｃに示す如く、振動子と反射体との間の距離
が異なるだけで、基準信号に対する受信信号の位
相が大きく変化し、Ｂ／Ａ（Ｚ）・Ｐ（Ｚ）の影響
による位相シフトを隠してしまう。
第２に、生体組織内からの反射波は、第５図に
示す様に、反射波どうしの重畳が頻繁に起こつて
おり、第５図ｃに示す如く合成された受信信号は
元のいずれの信号ａ，ｂとも大きく異つた位相を
持つ事になつて、やはりＢ／Ａ（Ｚ）・Ｐ（Ｚ）の影 響による位相シフトを隠してしまう。
示す様に、反射波どうしの重畳が頻繁に起こつて
おり、第５図ｃに示す如く合成された受信信号は
元のいずれの信号ａ，ｂとも大きく異つた位相を
持つ事になつて、やはりＢ／Ａ（Ｚ）・Ｐ（Ｚ）の影 響による位相シフトを隠してしまう。
この様な現象を除くため、先に本出願人が出願
した特願昭58−227949号においてはポンピングパ
ルスと共に測定用パルスを送信した時の受信信号
の位相と、ポンピングパルスを送らずに測定用パ
ルスだけを送信した時の受信信号の位相との差を
求める様にしている。上記第４図又は第５図に示
した如き、伝搬又はパルス重畳による位相の変化
は、ポンピングパルスの存在する時も存在しない
時も同様に現れるから、それぞれの場合の位相復
調出力を記憶しておいてその差を得る事により、
伝搬又はパルス重畳の影響を除いた、Ｂ／Ａ（Ｚ）・ Ｐ（Ｚ）の影響のみによる位相変化Φ（Ｚ）を抽出
する事ができる。
した特願昭58−227949号においてはポンピングパ
ルスと共に測定用パルスを送信した時の受信信号
の位相と、ポンピングパルスを送らずに測定用パ
ルスだけを送信した時の受信信号の位相との差を
求める様にしている。上記第４図又は第５図に示
した如き、伝搬又はパルス重畳による位相の変化
は、ポンピングパルスの存在する時も存在しない
時も同様に現れるから、それぞれの場合の位相復
調出力を記憶しておいてその差を得る事により、
伝搬又はパルス重畳の影響を除いた、Ｂ／Ａ（Ｚ）・ Ｐ（Ｚ）の影響のみによる位相変化Φ（Ｚ）を抽出
する事ができる。
これに対して本発明においては、第６図に示し
た如く、ポンピング波の高音圧部分（正音圧）に
測定波を重畳させた場合と、ポンピング波の低音
圧部分（負音圧）に測定波を重畳させた場合との
測定波の位相変化量を比較する事により、特願昭
58−227949号で述べたポンピング波の有無による
測定波の位相変化量の比較の場合よりも、高感度
にＢ／Ａの分布が検出できる様にしている。
た如く、ポンピング波の高音圧部分（正音圧）に
測定波を重畳させた場合と、ポンピング波の低音
圧部分（負音圧）に測定波を重畳させた場合との
測定波の位相変化量を比較する事により、特願昭
58−227949号で述べたポンピング波の有無による
測定波の位相変化量の比較の場合よりも、高感度
にＢ／Ａの分布が検出できる様にしている。
以下図を用いて実施例の構成例を、説明する。
第７図において、１はポンピングパルスの送信の
タイミングを発生するタイミング制御部、２は測
定波バーストパルスのための連続波発振器、３は
ポンピングパルス用のドライバ、４はポンピング
パルスの特定の位相（通常は音圧が最大又は最小
の位相）に測定用バーストパルスが重畳される様
に発振器２の出力を切り出すためのゲート回路、
５は測定用振動子７を駆動するためのドライバ、
６はポンピングパルス発生用の振動子、７は測定
用パルス発生用の振動子、８は被測定超音波媒
体、９は受信増幅器、１０は位相検出回路、１１
は位相検出器１０の出力を一時記憶するためのメ
モリ、１２はメモリ１１の出力から位相検出器１
０の出力を引算するための引算回路、１３は微分
回路、１４はいわゆるTime Gain Control回路
である。第７図ｂ，ｃには、振動子６，７の構成
例を、見取り図（第７図ｂ）及び断面図（同図
ｃ）をもつて示してある。第８図には第７図図示
の主要部における時間波形を示してある。第８図
ａのパルスが第７図のポンピングパルス用ドライ
バ３に与えられると、第８図ｃの波形が第７図の
ポンピングパルス発生用振動子６に印加される。
他方、第７図の連続波発振器２の出力は、ゲート
回路４によつて、例えばポンピングパルスの高音
圧部分にだけ測定用パルスが送出される様なタイ
ミング（第８図ｂ左端）だけ切り出され、ドライ
バ５を通つて、第８図ｄの如き駆動信号が第７図
の測定用振動子７に印加される。この結果、ポン
ピングパルスと測定用パルスとは、例えば第６図
ａに示した如きタイミング関係を保ちながら第７
図の被測定超音波媒体８の中を進行する事にな
る。被測定媒体中からの反射波は振動子７で受信
され、受信増幅器９で第８図ｅの如く増幅された
後、位相検出回路１０に入力される。位相検出回
路のもう一方の入力には連続波発振器２の出力が
基準信号として与えられており、第８図ｆの如く
この２つの入力間の位相差が位相検出回路１０か
ら出力されて記憶回路１１に送られる。
第７図において、１はポンピングパルスの送信の
タイミングを発生するタイミング制御部、２は測
定波バーストパルスのための連続波発振器、３は
ポンピングパルス用のドライバ、４はポンピング
パルスの特定の位相（通常は音圧が最大又は最小
の位相）に測定用バーストパルスが重畳される様
に発振器２の出力を切り出すためのゲート回路、
５は測定用振動子７を駆動するためのドライバ、
６はポンピングパルス発生用の振動子、７は測定
用パルス発生用の振動子、８は被測定超音波媒
体、９は受信増幅器、１０は位相検出回路、１１
は位相検出器１０の出力を一時記憶するためのメ
モリ、１２はメモリ１１の出力から位相検出器１
０の出力を引算するための引算回路、１３は微分
回路、１４はいわゆるTime Gain Control回路
である。第７図ｂ，ｃには、振動子６，７の構成
例を、見取り図（第７図ｂ）及び断面図（同図
ｃ）をもつて示してある。第８図には第７図図示
の主要部における時間波形を示してある。第８図
ａのパルスが第７図のポンピングパルス用ドライ
バ３に与えられると、第８図ｃの波形が第７図の
ポンピングパルス発生用振動子６に印加される。
他方、第７図の連続波発振器２の出力は、ゲート
回路４によつて、例えばポンピングパルスの高音
圧部分にだけ測定用パルスが送出される様なタイ
ミング（第８図ｂ左端）だけ切り出され、ドライ
バ５を通つて、第８図ｄの如き駆動信号が第７図
の測定用振動子７に印加される。この結果、ポン
ピングパルスと測定用パルスとは、例えば第６図
ａに示した如きタイミング関係を保ちながら第７
図の被測定超音波媒体８の中を進行する事にな
る。被測定媒体中からの反射波は振動子７で受信
され、受信増幅器９で第８図ｅの如く増幅された
後、位相検出回路１０に入力される。位相検出回
路のもう一方の入力には連続波発振器２の出力が
基準信号として与えられており、第８図ｆの如く
この２つの入力間の位相差が位相検出回路１０か
ら出力されて記憶回路１１に送られる。
記憶回路１１には、トリガ信号として、第８図
ａ図示のポンピング波駆動タイミングと、第８図
ｂ図示のゲートタイミングとが与えられており、
ポンピング波駆動タイミングの直後のゲートタイ
ミングから一走査線の時間（第８図ｆのＡの期
間）の間信号が記憶され、次のゲートタイミング
パルスで記憶内容を出力して引算回路１２に送り
出す。引算回路１２では、記憶回路１１の出力
（第８図ｆのＡ）から、位相検出回路１０の出力
（第８図ｆのＢ）を引き、ポンピング波の高音圧
部分に測定波が重畳した場合と、低音圧部分に重
畳した場合との位相検出回路の出力の差を第８図
ｇの如く得て、微分回路１３を経て、Time
Gain Control回路１４に送り出す。Time Gain
Control回路１４には、トリガ信号として第８図
ｂのゲートタイミング信号が与えられており、測
定用パルスの送信に同期して、第７図の微分回路
１３の出力（第８図ｈ）に対して、例えば、第３
図に示す如き時間的に変化する増幅が行われ、最
終出力Ｂ／Ａ（Ｚ）が得られる。尚、記憶回路１１ は、BBDやCCD等のデイレイラインによるアナ
ログ的なものでも、Ａ／Ｄ変換器とメモリ又はシ
フトレジスタを組合わせたデイジタル的なもので
も、どちらでも良い事は言うまでもない。
ａ図示のポンピング波駆動タイミングと、第８図
ｂ図示のゲートタイミングとが与えられており、
ポンピング波駆動タイミングの直後のゲートタイ
ミングから一走査線の時間（第８図ｆのＡの期
間）の間信号が記憶され、次のゲートタイミング
パルスで記憶内容を出力して引算回路１２に送り
出す。引算回路１２では、記憶回路１１の出力
（第８図ｆのＡ）から、位相検出回路１０の出力
（第８図ｆのＢ）を引き、ポンピング波の高音圧
部分に測定波が重畳した場合と、低音圧部分に重
畳した場合との位相検出回路の出力の差を第８図
ｇの如く得て、微分回路１３を経て、Time
Gain Control回路１４に送り出す。Time Gain
Control回路１４には、トリガ信号として第８図
ｂのゲートタイミング信号が与えられており、測
定用パルスの送信に同期して、第７図の微分回路
１３の出力（第８図ｈ）に対して、例えば、第３
図に示す如き時間的に変化する増幅が行われ、最
終出力Ｂ／Ａ（Ｚ）が得られる。尚、記憶回路１１ は、BBDやCCD等のデイレイラインによるアナ
ログ的なものでも、Ａ／Ｄ変換器とメモリ又はシ
フトレジスタを組合わせたデイジタル的なもので
も、どちらでも良い事は言うまでもない。
以上の如く、本発明によれば、超音波媒体の非
線形パラメータＢ／Ａ（Ｚ）の空間分布を、従来の 様な透過法でなく反射法で得る事ができ、透過法
の場合に必要な水槽等の大きな機械的構造を全く
必要としなくなると共に、種々の部位から対象を
観察できる様になり、操作性を大幅に向上させる
事ができ、しかも、特願昭58−227949号に示すも
のに比べて、Ｂ／Ａの検出感度を約２倍向上する
事ができる。
線形パラメータＢ／Ａ（Ｚ）の空間分布を、従来の 様な透過法でなく反射法で得る事ができ、透過法
の場合に必要な水槽等の大きな機械的構造を全く
必要としなくなると共に、種々の部位から対象を
観察できる様になり、操作性を大幅に向上させる
事ができ、しかも、特願昭58−227949号に示すも
のに比べて、Ｂ／Ａの検出感度を約２倍向上する
事ができる。
以上の実施例の説明は、特許請求の範囲第１項
および第３項および第４項に関するものであつ
た。特許請求の範囲第１項が、位相検出結果の間
の引算によつて、ポンピングパルスの高音圧部分
に測定波を重畳させた場合、及び、低音圧部分に
重畳させた場合との位相差を得ようとしていたの
に対して、特許請求の範囲第２項は、ポンピング
波の高音圧部分に測定波を重畳させた場合の受信
RF信号と、ポンピング波の低音圧部分に測定波
を重畳させた場合の受信RF信号とを、直接に位
相比較する事によつて、ポンピング波の高音圧部
分に重畳させた場合と低音圧部分に重畳させた場
合との、受信信号の位相差を求めようとするもの
である。構成例を第９図に示す。第９図におい
て、第７図と同じ構成要素には同一の番号を付し
てあり、説明を省略する。第７図と異なるのは、
位相検出回路１０と記憶回路１１の順序が逆にな
り、引算回路１２がなくなつた事である。第９図
の構成を用いる場合には、記憶回路１１として
は、RF受信信号を直接記憶できるだけの充分に
高速のアナログ（CCD等による）又はデイジタ
ル（Ａ／Ｄ変換器とメモリ又はシフトレジスタに
よる）的な記憶回路が必要である。第９図のシス
テムの動作は、第７図のシステムの動作から容易
に類推できるので説明は省略する。
および第３項および第４項に関するものであつ
た。特許請求の範囲第１項が、位相検出結果の間
の引算によつて、ポンピングパルスの高音圧部分
に測定波を重畳させた場合、及び、低音圧部分に
重畳させた場合との位相差を得ようとしていたの
に対して、特許請求の範囲第２項は、ポンピング
波の高音圧部分に測定波を重畳させた場合の受信
RF信号と、ポンピング波の低音圧部分に測定波
を重畳させた場合の受信RF信号とを、直接に位
相比較する事によつて、ポンピング波の高音圧部
分に重畳させた場合と低音圧部分に重畳させた場
合との、受信信号の位相差を求めようとするもの
である。構成例を第９図に示す。第９図におい
て、第７図と同じ構成要素には同一の番号を付し
てあり、説明を省略する。第７図と異なるのは、
位相検出回路１０と記憶回路１１の順序が逆にな
り、引算回路１２がなくなつた事である。第９図
の構成を用いる場合には、記憶回路１１として
は、RF受信信号を直接記憶できるだけの充分に
高速のアナログ（CCD等による）又はデイジタ
ル（Ａ／Ｄ変換器とメモリ又はシフトレジスタに
よる）的な記憶回路が必要である。第９図のシス
テムの動作は、第７図のシステムの動作から容易
に類推できるので説明は省略する。
特許請求の範囲第５項は、測定したＢ／Ａ（Ｚ）
分布のＳ／Ｎ比を向上させる手段に関するもので
ある。非線形パラメータと音圧との積による測定
用パルスの位相変化量が小さい時には、本来の位
相変化量Φ（Ｚ）に対して回路内その他で発生す
る雑音を無視できなくなる。更に、Ｂ／Ａ（Ｚ）を 求めるため、Φ（Ｚ）を微分するが、この微分動
作も雑音を大きくする原因となる。従つて得られ
たＢ／Ａ（Ｚ）には雑音が含まれ易い事になる。こ の様子を第１０図に示す。第１０図Ａ図示の無雑
音のＢ／Ａ（Ｚ）出力に、第１０図Ｂ図示の雑音Ｎ が加わる事により、第１０図Ｃ図示の如き信号が
実際には出力される。この対策として、同一部位
をＫ回測定して、第１０図Ｃ図示のS1ないしSK
の如き雑音の加わつた信号を行つて、こられを同
一Ｚ座標の点毎に加算すると、各点において信号
成分は振幅でＫ倍されるが、雑音成分は電力でＫ
倍されるに過ぎず、もし雑音が不規則雑音であれ
ば、各点における雑音振幅は√倍されるにとど
まる。従つて、Ｓ／Ｎ比は√倍改善されて、第
１０図Ｄ図示の如き出力が得られる。
ある。非線形パラメータと音圧との積による測定
用パルスの位相変化量が小さい時には、本来の位
相変化量Φ（Ｚ）に対して回路内その他で発生す
る雑音を無視できなくなる。更に、Ｂ／Ａ（Ｚ）を 求めるため、Φ（Ｚ）を微分するが、この微分動
作も雑音を大きくする原因となる。従つて得られ
たＢ／Ａ（Ｚ）には雑音が含まれ易い事になる。こ の様子を第１０図に示す。第１０図Ａ図示の無雑
音のＢ／Ａ（Ｚ）出力に、第１０図Ｂ図示の雑音Ｎ が加わる事により、第１０図Ｃ図示の如き信号が
実際には出力される。この対策として、同一部位
をＫ回測定して、第１０図Ｃ図示のS1ないしSK
の如き雑音の加わつた信号を行つて、こられを同
一Ｚ座標の点毎に加算すると、各点において信号
成分は振幅でＫ倍されるが、雑音成分は電力でＫ
倍されるに過ぎず、もし雑音が不規則雑音であれ
ば、各点における雑音振幅は√倍されるにとど
まる。従つて、Ｓ／Ｎ比は√倍改善されて、第
１０図Ｄ図示の如き出力が得られる。
この方法を用いたシステムの構成例を第１１図
に示す。第１１図において、第７図と同一の構成
要素には同一の番号を付してあり、説明は省略す
る。第１１図において、１６は、ポンピングパル
スの送信繰り返し周期Ｔだけ信号を遅らせる遅延
回路であり、例えばBBDやCCD等のアナログ的
手段で実現しても良いし、又、Ａ／Ｄ変換器とシ
フトレジスタ又はメモリとを用いたデイジタル的
手段で実現しても良い事は言うまでもない。１５
は加算器であり、遅延回路の種類に応じてアナロ
グ又はデイジタルのいずれのタイプでも良い。第
１１図図示の構成の場合、Time Gain Control
回路１４の出力が同一のＺ軸座標の各点毎に加算
される事は明らかであり、いわゆる同期加算によ
りＳ／Ｎ比の改善を行う事ができる。
に示す。第１１図において、第７図と同一の構成
要素には同一の番号を付してあり、説明は省略す
る。第１１図において、１６は、ポンピングパル
スの送信繰り返し周期Ｔだけ信号を遅らせる遅延
回路であり、例えばBBDやCCD等のアナログ的
手段で実現しても良いし、又、Ａ／Ｄ変換器とシ
フトレジスタ又はメモリとを用いたデイジタル的
手段で実現しても良い事は言うまでもない。１５
は加算器であり、遅延回路の種類に応じてアナロ
グ又はデイジタルのいずれのタイプでも良い。第
１１図図示の構成の場合、Time Gain Control
回路１４の出力が同一のＺ軸座標の各点毎に加算
される事は明らかであり、いわゆる同期加算によ
りＳ／Ｎ比の改善を行う事ができる。
本発明においては、更に特許請求の範囲第６項
に示す如く、Ｂ／Ａ（Ｚ）の二次元又は三次元分布 を得る事が可能となる。これまでの説明で明らか
なな様に、特定の走査線上のＢ／Ａ（Ｚ）の分布を 得る事ができるから、ポンピングパルス用振動子
６及び測定パルス用振動子７の相対位置を一定に
保つたままでｘまたは／およびｙ方向に移動さ
せ、各ｘまたは／およびｙ座標に対応したメモリ
アドレスにＢ／Ａ（Ｚ）の値を記憶しておけば、Ｂ／Ａ （Ｚ）の二次元または／および３次元分布を得る
事ができる。
に示す如く、Ｂ／Ａ（Ｚ）の二次元又は三次元分布 を得る事が可能となる。これまでの説明で明らか
なな様に、特定の走査線上のＢ／Ａ（Ｚ）の分布を 得る事ができるから、ポンピングパルス用振動子
６及び測定パルス用振動子７の相対位置を一定に
保つたままでｘまたは／およびｙ方向に移動さ
せ、各ｘまたは／およびｙ座標に対応したメモリ
アドレスにＢ／Ａ（Ｚ）の値を記憶しておけば、Ｂ／Ａ （Ｚ）の二次元または／および３次元分布を得る
事ができる。
例えば第１２図に示す如く、ポンピングパルス
用振動子Xpと測定パルス用振動子Xnとをそれぞ
れアレイ振動子としたものを使い、いわゆる電子
スキヤンを行えば、容易にＢ／Ａ（Ｚ）の二次元分 布を得る事ができるのは明らかである。又第１２
図図示の振動子アレイを図のｙ軸方向に拡張し
て、例えば第１３図の様な構成の振動子アレイを
用いれば、容易にＢ／Ａ（Ｚ）の三次元分布を得る 事ができるのは明らかである。
用振動子Xpと測定パルス用振動子Xnとをそれぞ
れアレイ振動子としたものを使い、いわゆる電子
スキヤンを行えば、容易にＢ／Ａ（Ｚ）の二次元分 布を得る事ができるのは明らかである。又第１２
図図示の振動子アレイを図のｙ軸方向に拡張し
て、例えば第１３図の様な構成の振動子アレイを
用いれば、容易にＢ／Ａ（Ｚ）の三次元分布を得る 事ができるのは明らかである。
以上の説明では、二次元又は三次元のＢ／Ａ（Ｚ）
の分布を得るのに、第１２図または第１３図の如
き振動子アレイを用いる場合について説明した
が、第７図ｂの如き振動子を用いて機械的に一次
元又は二次元に走査する事により、二次元又は三
次元のＢ／Ａ（Ｚ）の分布を得ても良い事は言うま でもない。
き振動子アレイを用いる場合について説明した
が、第７図ｂの如き振動子を用いて機械的に一次
元又は二次元に走査する事により、二次元又は三
次元のＢ／Ａ（Ｚ）の分布を得ても良い事は言うま でもない。
又、以上の説明では、ポンピングパルス用振動
子と測定パルス用振動子とは、横方向にずれた位
置に隣合わせに設けられていたが、特許請求の範
囲第７項に示した様に、例えば第１４図に示す如
く、ポンピングパルス用の振動子Xpと測定パル
ス用の振動子Xnとを重合わせても良い。この場
合、第１４図の測定パルス用振動子Xnとしては
例えば柔軟で音響インピーダンスが生体組織のそ
れに近いいわゆるPVDFを、ポンピングパルス用
振動子Xpとしては例えば硬くて音響インピーダ
ンスの大きないわゆるPZTを用いると、ポンピ
ングパルスはPVDF層を殆ど減衰なく通過でき、
しかもPZT層はPVDF層に対してパツキング層
BKとして作用する等、極めて具合の良い振動子
構造を実現できる。
子と測定パルス用振動子とは、横方向にずれた位
置に隣合わせに設けられていたが、特許請求の範
囲第７項に示した様に、例えば第１４図に示す如
く、ポンピングパルス用の振動子Xpと測定パル
ス用の振動子Xnとを重合わせても良い。この場
合、第１４図の測定パルス用振動子Xnとしては
例えば柔軟で音響インピーダンスが生体組織のそ
れに近いいわゆるPVDFを、ポンピングパルス用
振動子Xpとしては例えば硬くて音響インピーダ
ンスの大きないわゆるPZTを用いると、ポンピ
ングパルスはPVDF層を殆ど減衰なく通過でき、
しかもPZT層はPVDF層に対してパツキング層
BKとして作用する等、極めて具合の良い振動子
構造を実現できる。
(E) 発明の効果
以上述べた如く、本発明によれば、以前に提案
した超音波媒体の非線形パラメータＢ／Ａ（Ｚ）の 測定法に比べて、よりコンパクトな装置構成で、
しかもより操作性の優れた反射法によつて、Ｂ／Ａ （Ｚ）の分布を測定する事ができ、更に、特願昭
58−227949号に示すものよりも、２倍程度高感度
な測定が可能となる。
した超音波媒体の非線形パラメータＢ／Ａ（Ｚ）の 測定法に比べて、よりコンパクトな装置構成で、
しかもより操作性の優れた反射法によつて、Ｂ／Ａ （Ｚ）の分布を測定する事ができ、更に、特願昭
58−227949号に示すものよりも、２倍程度高感度
な測定が可能となる。
第１図および第２図は本発明の前提概念を説明
する説明図、第３図はいわゆるTime Gain
Controlの概念を説明する説明図、第４図および
第５図は、反射波の位相の概念を説明する説明
図、第６図は本発明において用いられる測定用超
音波パルスとポンピング用超音波パルスとの送出
タイミングを説明する説明図、第７図は本発明の
一実施例構成、第８図は第７図図示構成の動作を
説明する説明図、第９図は第７図に対応する一実
施例構成、第１０図は同期加算の概念を説明する
ための説明図、第１１図は同期加算を行うための
一実施例構成、第１２図ないし第１４図は振動子
の種々の実施形態例を示す。図中、１はタイミン
グ制御部、２は連続波発振器、３はポンピングパ
ルス用ドライバ、４はゲート回路、５は測定パル
ス用ドライバ、６はポンピングパルス用振動子、
７は測定パルス用振動子、８は被測定超音波媒
体、９は受信増幅器、１０は位相検出回路、１１
は記憶回路、１２は引算回路、１３は微分回路、
１４はTime Gain Control回路、１５は加算回
路、１６は遅延回路を表す。
する説明図、第３図はいわゆるTime Gain
Controlの概念を説明する説明図、第４図および
第５図は、反射波の位相の概念を説明する説明
図、第６図は本発明において用いられる測定用超
音波パルスとポンピング用超音波パルスとの送出
タイミングを説明する説明図、第７図は本発明の
一実施例構成、第８図は第７図図示構成の動作を
説明する説明図、第９図は第７図に対応する一実
施例構成、第１０図は同期加算の概念を説明する
ための説明図、第１１図は同期加算を行うための
一実施例構成、第１２図ないし第１４図は振動子
の種々の実施形態例を示す。図中、１はタイミン
グ制御部、２は連続波発振器、３はポンピングパ
ルス用ドライバ、４はゲート回路、５は測定パル
ス用ドライバ、６はポンピングパルス用振動子、
７は測定パルス用振動子、８は被測定超音波媒
体、９は受信増幅器、１０は位相検出回路、１１
は記憶回路、１２は引算回路、１３は微分回路、
１４はTime Gain Control回路、１５は加算回
路、１６は遅延回路を表す。
Claims (1)
- 【特許請求の範囲】 １ (イ) 超音波媒体内に測定用の比較的高周波且
つ低音圧の超音波パルスを送出する送受信兼用
の超音波振動子と、 (ロ) 該測定用超音波パルスと殆ど同一場所から殆
ど同一方向にポンピング用の比較的低周波且つ
高音圧の超音波パルスを該超音波媒体内に送出
する超音波振動子と、 (ハ) 該測定用超音波パルスが該ポンピング用超音
波パルスの予め定められた複数の特定の位相の
うちのいずれか１つの部分に重畳する様に送出
されるべく、振動子の駆動タイミングを制御す
る手段と、 (ニ) 該測定用超音波パルスの該超音波媒体内から
の反射波、から成る受信信号の位相を検出する
手段と、 (ホ) ポンピング用超音波パルスの上記特定の位相
のうちのいずれか１つに重畳する様に測定用超
音波パルスを送信した時に検出した受信信号の
位相（時間の関数）とポンピング用超音波パル
スの別の特定位相の部分に重畳する様に測定用
超音波パルスを送出した時に検出した受信信号
の位相（時間の関数）との差（時間の関数）を
求める手段と、 (ヘ) (ホ)で求めた、時間の関数としての位相差を微
分する事により、該超音波媒体の等価非線形パ
ラメータの、該測定用ビーム上の空間的分布を
求める手段 を有する事を特徴とする超音波媒体の非線形パラ
メータ分布測定装置。 ２ 上記予め定められた複数の特定の位相は、ポ
ンピング波の高音圧部（正音圧）及び低音圧部
（負音圧）である事を特徴とする特許請求の範囲
第１項記載の超音波媒体の非線形パラメータ分布
測定装置。 ３ ポンピング用超音波の音圧が超音波媒体内部
に行くと共に小さくなる効果を補償する手段とし
て、第１項(ヘ)の微分出力の後段に、時間と共に利
得の増大する増幅器を設けた事を特徴とする特許
請求範囲第１項または第２項記載の超音波媒体の
非線形パラメータ分布測定装置。 ４ 同一の測定ビーム線上の同一位置に対応する
非線形パラメータ値を、複数回の超音波パルスの
送受信毎に加算するいわゆる同期加算回路を持つ
事を特徴とする第１項ないし第３項いずれか記載
の超音波媒体の非線形パラメータ分布測定装置。 ５ 測定用およびポンピング用の振動子の組を一
次元的又は二次元的に走査する事により、二次元
的又は三次元的な等価非線形パラメータ分布像を
得るための走査部・表示部を持つ事を特徴とする
第１項ないし第４項いずれか記載の超音波媒体の
非線形パラメータ分布測定装置。 ６ ポンピング用の振動子は材料としていわゆる
PZT等の圧電セラミツクを使用した単板または
アレイ振動子であり、測定用の振動子は材料とし
ていわゆるPVDFの如き有機圧電フイルムを使用
した単板又はアレイ振動子であつて、且つ、後者
は前者の前面を被覆する様に配置された振動子を
用いる事を特徴とする特許請求の範囲第１項ない
し第５項いずれか記載の超音波媒体の非線形パラ
メータ分布測定装置。 ７ (イ) 超音波媒体内に測定用の比較的高周波且
つ低音圧の超音波パルスを送出する送受信兼用
の超音波振動子と、 (ロ) 該測定用超音波パルスと殆ど同一場所から殆
ど同一方向にポンピング用の比較的低周波且つ
高音圧の超音波パルスを該超音波媒体内に送出
する超音波振動子と、 (ハ) 該測定用超音波パルスが該ポンピング用超音
波パルスの予め定められた複数の特定の位相の
うちのいずれか１つの部分に重畳する様に送出
されるべく、振動子の駆動タイミングを制御す
る手段と、 (ニ) ポンピング用超音波パルスの上記特定位相の
いずれか１つに重畳する様に測定用超音波パル
スを送信した時の受信信号（いわゆるRF信号）
とポンピング用超音波パルスの別の特定位相の
部分に重畳する様に測定用超音波パルスを送出
した時の受信信号（RF信号）との位相差をRF
信号間の位相比較を行つて求める手段と、 (ホ) (ニ)で求めた、時間の関数としての位相差を微
分する事により、該超音波媒体の等価非線形パ
ラメータの、該測定用ビーム上の空間的分布を
求める手段 を有することを特徴とする超音波媒体の非線形パ
ラメータ分布測定装置。 ８ 上記予め定められた複数の特定の位相は、ポ
ンピング波の高音圧部（正音圧）及び低音圧部
（負音圧）である事を特徴とする特許請求の範囲
第７項記載の超音波媒体の非線形パラメータ分布
測定装置。 ９ ポンピング用超音波の音圧が超音波媒体内部
に行くと共に小さくなる効果を補償する手段とし
て、第７項(ホ)の微分出力の後段に、時間と共に利
得の増大する増幅器を設けた事を特徴とする特許
請求範囲第７項または第８項記載の超音波媒体の
非線形パラメータ分布測定装置。 １０ 同一の測定ビーム線上の同一位置に対応す
る非線形パラメータ値を、複数回の超音波パルス
の送受信毎に加算するいわゆる同期加算回路を持
つ事を特徴とする第７項ないし第９項いずれか記
載の超音波媒体の非線形パラメータ分布測定装
置。 １１ 測定用およびポンピング用の振動子の組を
一次元的又は二次元的に走査する事により、二次
元的又は三次元的な等価非線形パラメータ分布像
を得るための走査部・表示部を持つ事を特徴とす
る第７項ないし第１０項いずれか記載の超音波媒
体の非線形パラメータ分布測定装置。 １２ ポンピング用の振動子は材料としていわゆ
るPZT等の圧電セラミツクを使用した単板また
はアレイ振動子であり、測定用の振動子は材料と
していわゆるPVDFの如き有機圧電フイルムを使
用した単板又はアレイ振動子であつて、且つ、後
者は前者の前面を被覆する様に配置された振動子
を用いる事を特徴とする特許請求の範囲第７項な
いし第１１項いずれか記載の超音波媒体の非線形
パラメータ分布測定装置。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59109831A JPS60253434A (ja) | 1984-05-30 | 1984-05-30 | 超音波媒体の非線形パラメ−タ分布測定装置 |
| US06/675,612 US4610255A (en) | 1983-12-02 | 1984-11-28 | Ultrasonic non-linear parameter measuring system |
| DE8484308324T DE3474211D1 (en) | 1983-12-02 | 1984-11-30 | Ultrasonic non-linear parameter measuring system |
| EP84308324A EP0147955B1 (en) | 1983-12-02 | 1984-11-30 | Ultrasonic non-linear parameter measuring system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP59109831A JPS60253434A (ja) | 1984-05-30 | 1984-05-30 | 超音波媒体の非線形パラメ−タ分布測定装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS60253434A JPS60253434A (ja) | 1985-12-14 |
| JPH0331059B2 true JPH0331059B2 (ja) | 1991-05-02 |
Family
ID=14520308
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP59109831A Granted JPS60253434A (ja) | 1983-12-02 | 1984-05-30 | 超音波媒体の非線形パラメ−タ分布測定装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS60253434A (ja) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61228841A (ja) * | 1985-04-04 | 1986-10-13 | アロカ株式会社 | 超音波測定装置 |
| JPH066124B2 (ja) * | 1986-03-28 | 1994-01-26 | 松下電器産業株式会社 | 超音波計測装置 |
| JP2629734B2 (ja) * | 1987-09-24 | 1997-07-16 | 株式会社島津製作所 | 超音波物体検査装置 |
| JPH02191441A (ja) * | 1989-01-19 | 1990-07-27 | Gijutsu Kenkyu Kumiai Iryo Fukushi Kiki Kenkyusho | 超音波計測装置 |
| WO2006009469A2 (en) * | 2004-07-23 | 2006-01-26 | Angelsen Bjoern A J | Ultrasound imaging using non-linear manipulation of forward propagation |
| FR2905180B1 (fr) * | 2006-08-22 | 2008-12-05 | Super Sonic Imagine | Sonde d'imagerie ultrasonore pour imager une modification transitoire d'un milieu |
| US9291493B2 (en) * | 2011-10-03 | 2016-03-22 | Surf Technology As | Nonlinear imaging with dual band pulse complexes |
-
1984
- 1984-05-30 JP JP59109831A patent/JPS60253434A/ja active Granted
Also Published As
| Publication number | Publication date |
|---|---|
| JPS60253434A (ja) | 1985-12-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP0147955B1 (en) | Ultrasonic non-linear parameter measuring system | |
| EP0128635B1 (en) | Measuring non-linear parameter of acoustic medium | |
| US5632277A (en) | Ultrasound imaging system employing phase inversion subtraction to enhance the image | |
| US7785259B2 (en) | Detection of motion in vibro-acoustography | |
| US5903516A (en) | Acoustic force generator for detection, imaging and information transmission using the beat signal of multiple intersecting sonic beams | |
| US10222461B2 (en) | Phase inversion ultrasonic imaging | |
| JP4150866B2 (ja) | 超音波イメージング・システムの作動方法 | |
| IL298685A (en) | Ultrasonic imaging device with programmable anatomy and flow imaging | |
| US20070038095A1 (en) | Ultrasound vibrometry | |
| JP3510025B2 (ja) | 超音波撮像装置 | |
| JP2002113003A (ja) | スペックル・サイズ推定を用いた超音波ベースの定量的運動測定 | |
| US4679565A (en) | Ultrasonic diagnostic apparatus using non-linear parameters of an organ | |
| JPS62123355A (ja) | 媒体の超音波エコーグラフィック検査方法及び装置 | |
| JP2005506124A (ja) | 身体の組織の力学的剛性を示す装置と方法 | |
| JPH0331059B2 (ja) | ||
| JPH0222660B2 (ja) | ||
| US8118745B2 (en) | Ultrasound imaging apparatus | |
| Tabaru et al. | Measurement of elastic properties of tissue by shear wave propagation generated by acoustic radiation force | |
| JP5075503B2 (ja) | 超音波診断装置 | |
| Sandrin et al. | 2D transient elastography | |
| JPS62179437A (ja) | 超音波非線形パラメ−タ分布測定装置 | |
| JP2702983B2 (ja) | 超音波診断装置 | |
| JP2002301071A (ja) | 超音波撮像方法及び装置 | |
| JP2004041324A (ja) | 超音波診断装置 | |
| JP2002301072A (ja) | 超音波撮像方法及び装置 |