JPH0614929B2 - 超音波計測装置 - Google Patents
超音波計測装置Info
- Publication number
- JPH0614929B2 JPH0614929B2 JP1174370A JP17437089A JPH0614929B2 JP H0614929 B2 JPH0614929 B2 JP H0614929B2 JP 1174370 A JP1174370 A JP 1174370A JP 17437089 A JP17437089 A JP 17437089A JP H0614929 B2 JPH0614929 B2 JP H0614929B2
- Authority
- JP
- Japan
- Prior art keywords
- unit
- ultrasonic
- subject
- signal
- pulse
- 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 - Fee Related
Links
- 239000000523 sample Substances 0.000 claims description 52
- 238000001228 spectrum Methods 0.000 claims description 40
- 239000000463 material Substances 0.000 claims description 26
- 230000001902 propagating effect Effects 0.000 claims description 15
- 238000006243 chemical reaction Methods 0.000 claims description 13
- 238000001514 detection method Methods 0.000 claims description 9
- 239000000284 extract Substances 0.000 claims description 5
- 238000010586 diagram Methods 0.000 description 16
- 239000002245 particle Substances 0.000 description 14
- 238000012545 processing Methods 0.000 description 12
- 230000001133 acceleration Effects 0.000 description 11
- 238000000034 method Methods 0.000 description 10
- 210000001519 tissue Anatomy 0.000 description 9
- 238000005259 measurement Methods 0.000 description 7
- 238000012546 transfer Methods 0.000 description 6
- 238000001727 in vivo Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 2
- 230000009022 nonlinear effect Effects 0.000 description 2
- 238000001028 reflection method Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 210000004303 peritoneum Anatomy 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Description
【発明の詳細な説明】 産業上の利用分野 本発明は、被検体内に超音波を送信し、被検体内からの
反射波を受信し、この反射波に包含する被検体内の音響
特性を計測する超音波計測装置に関するものである。
反射波を受信し、この反射波に包含する被検体内の音響
特性を計測する超音波計測装置に関するものである。
従来の技術 従来、被検体内の音響的情報を超音波を用いて得る方式
の例としては、超音波診断装置がある。この超音波診断
装置としては、被検体である生体内に超音波を送信し、
生体内からの反射波を受信し、この反射波に包含する生
体内の情報を得るようにしたパルス反射法を用いるもの
が主流を占めている。このパルス反射法は、通常、生体
内の音響インピーダンス差のある界面からの反射エコー
強度、すなわち振幅値と超音波の伝搬時間とから生体内
の情報を2次元的に集めて表示することにより断層像を
得るようになっている。しかし、近年、主に生体内組織
の形状判断を行う超音波診断装置に対し、生体内組織の
形状のみならず、質の情報をも得たいという要望が高ま
っている。このような生体組織の質に関する情報は、例
えば、生体内の各種臓器で特有の値を有する超音波の減
衰の大きさ、音速等を計測することにより得ることがで
きる。この超音波の減衰係数を計測する超音波計測装置
として、例えばウルトラソニック・イメージング:UL
TRASONIC IMAGING,VOL.5,N
O.2,1983,117〜135頁に記載されている構成が知ら
れている。以下、第7図を参照しながら上記従来の超音
波計測装置について説明する。
の例としては、超音波診断装置がある。この超音波診断
装置としては、被検体である生体内に超音波を送信し、
生体内からの反射波を受信し、この反射波に包含する生
体内の情報を得るようにしたパルス反射法を用いるもの
が主流を占めている。このパルス反射法は、通常、生体
内の音響インピーダンス差のある界面からの反射エコー
強度、すなわち振幅値と超音波の伝搬時間とから生体内
の情報を2次元的に集めて表示することにより断層像を
得るようになっている。しかし、近年、主に生体内組織
の形状判断を行う超音波診断装置に対し、生体内組織の
形状のみならず、質の情報をも得たいという要望が高ま
っている。このような生体組織の質に関する情報は、例
えば、生体内の各種臓器で特有の値を有する超音波の減
衰の大きさ、音速等を計測することにより得ることがで
きる。この超音波の減衰係数を計測する超音波計測装置
として、例えばウルトラソニック・イメージング:UL
TRASONIC IMAGING,VOL.5,N
O.2,1983,117〜135頁に記載されている構成が知ら
れている。以下、第7図を参照しながら上記従来の超音
波計測装置について説明する。
第7図において、101は超音波探触子であり、被検体A
に対して超音波を送信し、被検体Aからの反射波を受信
する超音波変換器102と、超音波変換器102から出力され
る超音波を伝搬させる伝搬材103と、超音波変換器102を
2次元に走査する駆動モータ104等から構成されてい
る。105は超音波変換器102を駆動するパルス駆動器、10
6は超音波変換器102の受信信号を増幅する受信回路、10
7は受信回路106の出力から周波数分析を行う周波数分析
器、108は周波数分析器107の出力から演算を行う演算
部、109は演算部108の出力を表示するTVモニタであ
る。
に対して超音波を送信し、被検体Aからの反射波を受信
する超音波変換器102と、超音波変換器102から出力され
る超音波を伝搬させる伝搬材103と、超音波変換器102を
2次元に走査する駆動モータ104等から構成されてい
る。105は超音波変換器102を駆動するパルス駆動器、10
6は超音波変換器102の受信信号を増幅する受信回路、10
7は受信回路106の出力から周波数分析を行う周波数分析
器、108は周波数分析器107の出力から演算を行う演算
部、109は演算部108の出力を表示するTVモニタであ
る。
以上の構成において、以下、その動作について説明す
る。
る。
まず、パルス駆動器105より送出された駆動パルスが超
音波変換器102に加えられる。超音波変換器102は、この
駆動パルスを超音波パルスに変換し、伝搬材103を介し
て被検体A内に送出する。送出された超音波パルスは被
検体A内を伝搬しながら組織の音響的性質に対応して次
々に散乱され、その一部は伝搬径路、すなわち音響走査
線上を逆行して超音波変換器102へ到達し、受信信号に
変換される。この過程で超音波パルスは生体組織の超音
波減衰特性や超音波散乱特性の影響を受ける。受信信号
は受信回路106で増幅され、周波数分析器107で中心周波
数が求められる。中心周波数の求め方としては、ゼロク
ロスカウンタ、あるいは高速フーリエ変換手法(FF
T)等を利用した方法がある。この中心周波数は、伝搬
距離や超音波の減衰係数に依存する値である。演算部10
8では周波数分析器107の出力である中心周波数から超音
波の減衰係数を求める。そして、超音波変換器102を駆
動モータ104で機械的に走査させることにより、被検体
A内の2次元の中心周波数や減衰係数の分布を求めるこ
とが可能となり、これをTVモニタ109に表示する。
音波変換器102に加えられる。超音波変換器102は、この
駆動パルスを超音波パルスに変換し、伝搬材103を介し
て被検体A内に送出する。送出された超音波パルスは被
検体A内を伝搬しながら組織の音響的性質に対応して次
々に散乱され、その一部は伝搬径路、すなわち音響走査
線上を逆行して超音波変換器102へ到達し、受信信号に
変換される。この過程で超音波パルスは生体組織の超音
波減衰特性や超音波散乱特性の影響を受ける。受信信号
は受信回路106で増幅され、周波数分析器107で中心周波
数が求められる。中心周波数の求め方としては、ゼロク
ロスカウンタ、あるいは高速フーリエ変換手法(FF
T)等を利用した方法がある。この中心周波数は、伝搬
距離や超音波の減衰係数に依存する値である。演算部10
8では周波数分析器107の出力である中心周波数から超音
波の減衰係数を求める。そして、超音波変換器102を駆
動モータ104で機械的に走査させることにより、被検体
A内の2次元の中心周波数や減衰係数の分布を求めるこ
とが可能となり、これをTVモニタ109に表示する。
発明が解決しようとする課題 しかし、上記のような従来例の構成では、被検体A内の
2次元に中心周波数や減衰係数の分布を求める際、伝搬
材103部分に相当する位置や、伝搬材103と被検体Aの境
界からの強い反射波を含む部分、または被検体Aが生体
であった場合、腹膜部分の複雑な組織からの多重反射を
含む部分も周波数分析器107で周波数分析してしまい、
誤った情報をTVモニタ109に表示してしまうという課
題があった。
2次元に中心周波数や減衰係数の分布を求める際、伝搬
材103部分に相当する位置や、伝搬材103と被検体Aの境
界からの強い反射波を含む部分、または被検体Aが生体
であった場合、腹膜部分の複雑な組織からの多重反射を
含む部分も周波数分析器107で周波数分析してしまい、
誤った情報をTVモニタ109に表示してしまうという課
題があった。
本発明は、以上のような従来技術の課題を解決するもの
であり、精度良く被検体内の音響特性を計測することが
できるようにした超音波計測装置を提供することを目的
とするものである。
であり、精度良く被検体内の音響特性を計測することが
できるようにした超音波計測装置を提供することを目的
とするものである。
課題を解決するための手段 上記課題を解決するための本発明の技術的手段は、 超音波を被検体に送信するとともに、当該被検体からの
反射超音波信号を受信する超音波変換器と、 前記超音波変換器を内包し当該超音波変換器から送受波
される超音波を伝搬させる伝搬材と前記超音波変換器を
2次元走査させる駆動モータからなる超音波探触子と、 前記超音波変換器を電気的に駆動するパルス駆動器と、 前記超音波変換器を介し、被検体からの反射超音波信号
を増幅する受信回路部と、 前記受信回路部の出力信号をアナログ/デジタル変換す
るアナログ/デジタル変換部と、 前記アナログ/デジタル変換部の出力信号である反射超
音波信号を記憶するバッファメモリ部と、 前記バッファメモリ部に記憶されている反射超音波信号
をパワースペクトラム値として周波数分析する周波数分
析部と、 前記周波数分析部の解析結果であるパワースペクトラム
値を、前記伝搬材と被検体との境界位置を示す予め設定
した値と比較し、当該予め設定した値が前記パワースペ
クトラム値より小さいことを検出すると、当該検出結果
が前記伝搬材における反射超音波信号であったとしてア
ドレス算出部へ検出信号を出力するパワースペクトラム
値比較部と、 前記パワースペクトラム値比較部の検出信号が入力され
ると、前記伝搬材における反射超音波信号であるとして
アドレス算出を中止するとともに、一方、前記パワース
ペクトラム値比較部の検出信号が検出されないと、被検
体からの反射超音波信号であるとしてアドレス算出を行
うことにより、被検体からの反射超音波信号のみを抜き
出してアドレス生成するアドレス算出部と、 前記周波数分析部の出力を、前記アドレス算出部がアド
レス生成したアドレスに基づき、当該アドレスに記憶す
るX−Y表示メモリ部と、 前記X−Y表示メモリ部に記憶されている周波数結果を
表示する表示部とを設けたものである。
反射超音波信号を受信する超音波変換器と、 前記超音波変換器を内包し当該超音波変換器から送受波
される超音波を伝搬させる伝搬材と前記超音波変換器を
2次元走査させる駆動モータからなる超音波探触子と、 前記超音波変換器を電気的に駆動するパルス駆動器と、 前記超音波変換器を介し、被検体からの反射超音波信号
を増幅する受信回路部と、 前記受信回路部の出力信号をアナログ/デジタル変換す
るアナログ/デジタル変換部と、 前記アナログ/デジタル変換部の出力信号である反射超
音波信号を記憶するバッファメモリ部と、 前記バッファメモリ部に記憶されている反射超音波信号
をパワースペクトラム値として周波数分析する周波数分
析部と、 前記周波数分析部の解析結果であるパワースペクトラム
値を、前記伝搬材と被検体との境界位置を示す予め設定
した値と比較し、当該予め設定した値が前記パワースペ
クトラム値より小さいことを検出すると、当該検出結果
が前記伝搬材における反射超音波信号であったとしてア
ドレス算出部へ検出信号を出力するパワースペクトラム
値比較部と、 前記パワースペクトラム値比較部の検出信号が入力され
ると、前記伝搬材における反射超音波信号であるとして
アドレス算出を中止するとともに、一方、前記パワース
ペクトラム値比較部の検出信号が検出されないと、被検
体からの反射超音波信号であるとしてアドレス算出を行
うことにより、被検体からの反射超音波信号のみを抜き
出してアドレス生成するアドレス算出部と、 前記周波数分析部の出力を、前記アドレス算出部がアド
レス生成したアドレスに基づき、当該アドレスに記憶す
るX−Y表示メモリ部と、 前記X−Y表示メモリ部に記憶されている周波数結果を
表示する表示部とを設けたものである。
作用 本発明は上記構成によれば、前記周波数分析部の解析結
果であるパワースペクトラム値を、パワースペクトラム
値比較部で、前記伝搬材と被検体との境界位置を示す予
め設定した値と比較し、当該パワースペクトラム値より
予め設定した値が小さいことを検出すると、当該検出結
果が前記伝搬材における反射超音波信号であつたとして
当該反射超音波信号を表示させないことで、被検体から
の反射超音波信号のみを抜き出して表示させることがで
きる。
果であるパワースペクトラム値を、パワースペクトラム
値比較部で、前記伝搬材と被検体との境界位置を示す予
め設定した値と比較し、当該パワースペクトラム値より
予め設定した値が小さいことを検出すると、当該検出結
果が前記伝搬材における反射超音波信号であつたとして
当該反射超音波信号を表示させないことで、被検体から
の反射超音波信号のみを抜き出して表示させることがで
きる。
実施例 以下、本発明の実施例を図面に基づいて詳細に説明す
る。
る。
まず、本発明の第1の実施例について説明する。
第1図は本発明の第1の実施例における超音波計測装置
を示すブロック図である。第1図において、1は超音波
探触子であり、被検体Aに対して超音波を送信し、被検
体Aからの超音波を受信する超音波変換器2と、超音波
変換器2から出力される超音波を伝搬させる伝搬材3
と、超音波変換器2を2次元に走査する駆動モータ4等
から構成されている。5は超音波変換器2をパルス駆動
するパルス駆動器、6は超音波変換器2の受信信号を増
幅する受信回路部、7は受信回路部6で増幅された受信
信号をデジタル信号に変換するA/D変換部、8はA/
D変換部7でデジタル信号に変換された受信信号から伝
搬材3と被検体Aとの境界を算出する比較部、9はA/
D変換部7でデジタル信号に変換された受信信号を記憶
するバッファメモリ部、10は比較部8の出力である境界
信号よりバッファメモリ部9に記憶されている被検体A
内からの反射信号を抜き出すゲート位置発生部、11はバ
ッファメモリ部9から転送された時系列データに対して
周波数分析を行い、中心周波数等を算出し、被検体A内
の減衰特性や減衰特性に依存するパラメータを求める周
波数分析部、12はゲート位置発生部10から出力されるゲ
ート位置信号より位置アドレスを算出して出力するアド
レス算出部、13はアドレス算出部12から出力される位置
アドレス信号に基づき周波数分析部11からの出力を記憶
するX−Y表示メモリ部、14はX−Y表示メモリ部13の
記憶内容を表示するTVモニタ、15は超音波の送信タイ
ミング等を発生するタイミング発生部である。
を示すブロック図である。第1図において、1は超音波
探触子であり、被検体Aに対して超音波を送信し、被検
体Aからの超音波を受信する超音波変換器2と、超音波
変換器2から出力される超音波を伝搬させる伝搬材3
と、超音波変換器2を2次元に走査する駆動モータ4等
から構成されている。5は超音波変換器2をパルス駆動
するパルス駆動器、6は超音波変換器2の受信信号を増
幅する受信回路部、7は受信回路部6で増幅された受信
信号をデジタル信号に変換するA/D変換部、8はA/
D変換部7でデジタル信号に変換された受信信号から伝
搬材3と被検体Aとの境界を算出する比較部、9はA/
D変換部7でデジタル信号に変換された受信信号を記憶
するバッファメモリ部、10は比較部8の出力である境界
信号よりバッファメモリ部9に記憶されている被検体A
内からの反射信号を抜き出すゲート位置発生部、11はバ
ッファメモリ部9から転送された時系列データに対して
周波数分析を行い、中心周波数等を算出し、被検体A内
の減衰特性や減衰特性に依存するパラメータを求める周
波数分析部、12はゲート位置発生部10から出力されるゲ
ート位置信号より位置アドレスを算出して出力するアド
レス算出部、13はアドレス算出部12から出力される位置
アドレス信号に基づき周波数分析部11からの出力を記憶
するX−Y表示メモリ部、14はX−Y表示メモリ部13の
記憶内容を表示するTVモニタ、15は超音波の送信タイ
ミング等を発生するタイミング発生部である。
以上のような構成において、以下、その動作について説
明する。
明する。
パルス駆動器5より送出された超音波駆動パルスが超音
波変換器2に加えられる。超音波変換器2はこの超音波
駆動パルスを超音波パルスに変換し、伝搬材3を介して
被検体Aに送出する。被検体Aに送出された超音波パル
スは被検体A内を伝搬しながら組織の音響的性質に対応
して次々に散乱され、その一部は伝搬径路、すなわち音
響走査線上を逆行して超音波変換器2へ到達し、受信信
号に変換される。受信信号は受信回路部6で増幅され、
A/D変換部7でデジタル信号に変換される。デジタル
信号に変換された受信信号は順次バッファメモリ部9に
記憶される。また、A/D変換部7でデジタル信号に変
換された受信信号は比較部8に入力され、比較部8にお
いて、上記受信信号から伝搬材3と被検体Aとの境界を
算出する。すなわち、伝搬材3からの反射はほとんどな
く、被検体Aとの境界からは強い反射があり、A/D変
換部7の入力は大きなものとなる。このA/D変換部7
の出力が大きく変化する部分を比較部8であらかじめ設
定しておいた値と比較し、その設定値より大きくなった
位置で境界を判断することができる。ゲート位置発生部
10は、比較部8の出力である境界位置信号よりバッファ
メモリ部9に記憶している被検体A内からの反射信号成
分だけを抜き出し、周波数分析部11に転送する。周波数
分析部11は転送された時系列データに対し、高速フーリ
エ変換手法(FFT)等により周波数分析を行い、中心
周波数等を算出し、被検体A内の減衰特性や減衰特性に
依存するパラメータを求め、X−Y表示メモリ部13に書
き込む。X−Y表示メモリ部13の書き込む位置は、ゲー
ト位置発生部10の出力であるゲート位置よりアドレス算
出部12がX−Y表示メモリ部13の格子配列に整合するよ
うに算出された位置アドレスに基づく。そして、超音波
変換器2を駆動モータ4で2次元に走査させることによ
り、2次元の周波数分析の結果をX−Y表示メモリ部13
に記憶させ、TVモニタ14に表示させることができる。
波変換器2に加えられる。超音波変換器2はこの超音波
駆動パルスを超音波パルスに変換し、伝搬材3を介して
被検体Aに送出する。被検体Aに送出された超音波パル
スは被検体A内を伝搬しながら組織の音響的性質に対応
して次々に散乱され、その一部は伝搬径路、すなわち音
響走査線上を逆行して超音波変換器2へ到達し、受信信
号に変換される。受信信号は受信回路部6で増幅され、
A/D変換部7でデジタル信号に変換される。デジタル
信号に変換された受信信号は順次バッファメモリ部9に
記憶される。また、A/D変換部7でデジタル信号に変
換された受信信号は比較部8に入力され、比較部8にお
いて、上記受信信号から伝搬材3と被検体Aとの境界を
算出する。すなわち、伝搬材3からの反射はほとんどな
く、被検体Aとの境界からは強い反射があり、A/D変
換部7の入力は大きなものとなる。このA/D変換部7
の出力が大きく変化する部分を比較部8であらかじめ設
定しておいた値と比較し、その設定値より大きくなった
位置で境界を判断することができる。ゲート位置発生部
10は、比較部8の出力である境界位置信号よりバッファ
メモリ部9に記憶している被検体A内からの反射信号成
分だけを抜き出し、周波数分析部11に転送する。周波数
分析部11は転送された時系列データに対し、高速フーリ
エ変換手法(FFT)等により周波数分析を行い、中心
周波数等を算出し、被検体A内の減衰特性や減衰特性に
依存するパラメータを求め、X−Y表示メモリ部13に書
き込む。X−Y表示メモリ部13の書き込む位置は、ゲー
ト位置発生部10の出力であるゲート位置よりアドレス算
出部12がX−Y表示メモリ部13の格子配列に整合するよ
うに算出された位置アドレスに基づく。そして、超音波
変換器2を駆動モータ4で2次元に走査させることによ
り、2次元の周波数分析の結果をX−Y表示メモリ部13
に記憶させ、TVモニタ14に表示させることができる。
次に本発明の第2の実施例について説明する。
第2図は本発明の第2の実施例における超音波計測装置
を示すブロック図である。本実施例においては、上記第
1の実施例と同一部分については同一符号を付してその
説明を省略し、異なる構成について説明する。本実施例
の特徴とするところは、第2図から明らかなように上記
第1の実施例における比較部8に替えて周波数分析部11
とアドレス算出部12との間にパワースペクトル値比較部
16が接続され、このパワースペクトル値比較部16は周波
数分析部11から転送されたある周波数、例えば超音波駆
動パルスの中心周波数でのパワースペクトル値をあらか
じめ設定された値と比較し、パワースペクトル値が設定
値より小さい場合、アドレス算出部12へ書き込みアドレ
スの出力を禁止する制御信号を出力する。
を示すブロック図である。本実施例においては、上記第
1の実施例と同一部分については同一符号を付してその
説明を省略し、異なる構成について説明する。本実施例
の特徴とするところは、第2図から明らかなように上記
第1の実施例における比較部8に替えて周波数分析部11
とアドレス算出部12との間にパワースペクトル値比較部
16が接続され、このパワースペクトル値比較部16は周波
数分析部11から転送されたある周波数、例えば超音波駆
動パルスの中心周波数でのパワースペクトル値をあらか
じめ設定された値と比較し、パワースペクトル値が設定
値より小さい場合、アドレス算出部12へ書き込みアドレ
スの出力を禁止する制御信号を出力する。
以上の構成において、以下、その動作について説明す
る。
る。
パルス駆動器5より送出された超音波駆動パルスが超音
波変換器2に加えられる。超音波変換器2はこの超音波
駆動パルスを超音波パルスに変換し、伝搬材3を介して
被検体Aに送出する。被検体Aに送出された超音波パル
スは被検体A内を伝搬しながら組織の音響的性質に対応
して次々に散乱され、その一部は伝搬径路、すなわち音
響走査線上を逆行して超音波変換器2へ到達し、受信信
号に変換される。受信信号は受信回路部6で増幅され、
A/D変換器7でデジタル信号に変換される。デジタル
信号に変換された受信信号は順次バッファメモリ部9に
記憶される。ゲート位置発生部10はバッファメモリ部に
記憶した受信信号を順次ある区間抜き出し、周波数分析
部11に転送する。周波数分析部11は転送された時系列デ
ータに対し、高速フーリエ変換手法(FFT)等により
周波数分析を行い、中心周波数等を算出し、被検体A内
の減衰特性や減衰特性に依存するパラメータを求め、X
−Y表示メモリ部13に書き込む。このとき、周波数分析
部11は周波数分析結果よりある周波数、例えば超音波駆
動パルスの中心周波数でのパワースペクトル値をパワー
スペクトル値比較部16へ転送する。パワースペクトル比
較部16は転送されてきた値があらかじめ設定された値よ
り小さい場合、アドレス算出部12へX−Y表示メモリ部
13に対する書き込みアドレスの出力を禁止する制御信号
を出力する。すなわち、伝搬材3からの反射信号等のパ
ワースペクトル値は小さいので、パワースペクトル値を
比較することにより被検体Aと伝搬材3の境界位置を指
し示すことができ、被検体A内からの反射信号成分だけ
を抜き出すことができる。そして、超音波変換器2を駆
動モータ4で2次元に走査させることにより、2次元の
周波数分析の結果をX−Y表示メモリ部13に記憶させ、
TVモニタ14に表示させることができる。
波変換器2に加えられる。超音波変換器2はこの超音波
駆動パルスを超音波パルスに変換し、伝搬材3を介して
被検体Aに送出する。被検体Aに送出された超音波パル
スは被検体A内を伝搬しながら組織の音響的性質に対応
して次々に散乱され、その一部は伝搬径路、すなわち音
響走査線上を逆行して超音波変換器2へ到達し、受信信
号に変換される。受信信号は受信回路部6で増幅され、
A/D変換器7でデジタル信号に変換される。デジタル
信号に変換された受信信号は順次バッファメモリ部9に
記憶される。ゲート位置発生部10はバッファメモリ部に
記憶した受信信号を順次ある区間抜き出し、周波数分析
部11に転送する。周波数分析部11は転送された時系列デ
ータに対し、高速フーリエ変換手法(FFT)等により
周波数分析を行い、中心周波数等を算出し、被検体A内
の減衰特性や減衰特性に依存するパラメータを求め、X
−Y表示メモリ部13に書き込む。このとき、周波数分析
部11は周波数分析結果よりある周波数、例えば超音波駆
動パルスの中心周波数でのパワースペクトル値をパワー
スペクトル値比較部16へ転送する。パワースペクトル比
較部16は転送されてきた値があらかじめ設定された値よ
り小さい場合、アドレス算出部12へX−Y表示メモリ部
13に対する書き込みアドレスの出力を禁止する制御信号
を出力する。すなわち、伝搬材3からの反射信号等のパ
ワースペクトル値は小さいので、パワースペクトル値を
比較することにより被検体Aと伝搬材3の境界位置を指
し示すことができ、被検体A内からの反射信号成分だけ
を抜き出すことができる。そして、超音波変換器2を駆
動モータ4で2次元に走査させることにより、2次元の
周波数分析の結果をX−Y表示メモリ部13に記憶させ、
TVモニタ14に表示させることができる。
次に、本発明の第3の実施例について説明する。第3図
は本発明の第3の実施例における超音波計測装置を示す
ブロック図である。本実施例においては、上記第1の実
施例と同一部分については同一符号を付してその説明を
省略し、異なる構成について説明する。本実施例の特徴
とするところは、第3図から明らかなように上記第1の
実施例における比較部8に替えてゲート位置発生部10と
アドレス算出部12にオペレータ部17が接続され、このオ
ペレータ部17はあらかじめ時間を設定し、ゲート位置発
生部10がこの設定した時間に対応する位置より深い区間
の受信信号をバッファメモリ部9から周波数分析部11に
転送するようになっている。
は本発明の第3の実施例における超音波計測装置を示す
ブロック図である。本実施例においては、上記第1の実
施例と同一部分については同一符号を付してその説明を
省略し、異なる構成について説明する。本実施例の特徴
とするところは、第3図から明らかなように上記第1の
実施例における比較部8に替えてゲート位置発生部10と
アドレス算出部12にオペレータ部17が接続され、このオ
ペレータ部17はあらかじめ時間を設定し、ゲート位置発
生部10がこの設定した時間に対応する位置より深い区間
の受信信号をバッファメモリ部9から周波数分析部11に
転送するようになっている。
以上の構成において、以下、その動作について説明す
る。
る。
パルス駆動器5より送出された超音波駆動パルスが超音
波変換器2に加えられる。超音波変換器2はこの超音波
駆動パルスを超音波パルスに変換し、伝搬材3を介して
被検体Aに送出する。被検体Aに送出された超音波パル
スは被検体A内を伝搬しながら組織の音響的性質に対応
して次々に散乱され、その一部は伝搬径路、すなわち音
響走査線上を逆行して超音波変換器2に到達し、受信信
号に変換される。受信信号は受信回路部6で増幅され、
A/D変換部7でデジタル信号に変換される。デジタル
信号に変換された受信信号は順次バッファメモリ部9に
記憶される。バッファメモリ部9に記憶された受信信号
は、ゲート位置発生部10によってある区間抜き出し、周
波数分析部11に転送する。このとき、ゲート位置発生部
10はあらかじめオペレータ部17において設定した時間に
対応する位置より深い区間の受信信号をバッファメモリ
部9から周波数分析部11に順次転送する。例えば、伝搬
材3と被検体Aの音速が約1500m/secで、ゲート区間発
生時間を超音波駆動パルスから130μsecとすると、約10
cmの深さに対応する位置データから順次周波数分析部11
に転送する。伝搬材3での超音波伝搬距離が10cm以内で
あれば、周波数分析部11に転送される受信信号は、被検
体A内からの受信信号成分だけとなる。周波数分析部11
は転送された時系列データに対し、高速フーリエ変換手
法(FFT)等により周波数分析を行い、中心周波数等
を算出し、被検体A内の減衰特性や減衰特性に依存する
パラメータを求め、X−Y表示メモリ部13においてアド
レス算出部12が出力した位置に書き込む。そして、超音
波変換器2を駆動モータ4で2次元に走査させることに
より、2次元の周波数分析の結果をX−Y表示メモリ部
13に記憶させ、TVモニタ14に表示させることができ
る。
波変換器2に加えられる。超音波変換器2はこの超音波
駆動パルスを超音波パルスに変換し、伝搬材3を介して
被検体Aに送出する。被検体Aに送出された超音波パル
スは被検体A内を伝搬しながら組織の音響的性質に対応
して次々に散乱され、その一部は伝搬径路、すなわち音
響走査線上を逆行して超音波変換器2に到達し、受信信
号に変換される。受信信号は受信回路部6で増幅され、
A/D変換部7でデジタル信号に変換される。デジタル
信号に変換された受信信号は順次バッファメモリ部9に
記憶される。バッファメモリ部9に記憶された受信信号
は、ゲート位置発生部10によってある区間抜き出し、周
波数分析部11に転送する。このとき、ゲート位置発生部
10はあらかじめオペレータ部17において設定した時間に
対応する位置より深い区間の受信信号をバッファメモリ
部9から周波数分析部11に順次転送する。例えば、伝搬
材3と被検体Aの音速が約1500m/secで、ゲート区間発
生時間を超音波駆動パルスから130μsecとすると、約10
cmの深さに対応する位置データから順次周波数分析部11
に転送する。伝搬材3での超音波伝搬距離が10cm以内で
あれば、周波数分析部11に転送される受信信号は、被検
体A内からの受信信号成分だけとなる。周波数分析部11
は転送された時系列データに対し、高速フーリエ変換手
法(FFT)等により周波数分析を行い、中心周波数等
を算出し、被検体A内の減衰特性や減衰特性に依存する
パラメータを求め、X−Y表示メモリ部13においてアド
レス算出部12が出力した位置に書き込む。そして、超音
波変換器2を駆動モータ4で2次元に走査させることに
より、2次元の周波数分析の結果をX−Y表示メモリ部
13に記憶させ、TVモニタ14に表示させることができ
る。
次に本発明の第4の実施例について説明する。第4(a)
ないし(c)は本発明の第4の実施例における超音波計測
装置を示し、第4図(a)は全体のブロック図、第4図(b)
は制御駆動部のブロック図、第4図(c)は信号処理部の
ブロック図である。第4図(a)において、18は低周波帯
の第1の超音波パルスであるポンプ波パルスを送出する
ポンプ波用超音波変換器、19は第1の超音波パルスより
周波数の高い第2の超音波パルスであるプローブ波パル
スを送出するプローブ波用超音波変換器、20はポンプ波
用超音波変換器18およびプローブ波用超音波変換器19を
位相制御駆動する制御駆動部であり、この制御駆動部20
は第4図(b)から明らかなようにポンプ波用超音波変換
器18を駆動するパルス駆動器31と、プローブ波用超音波
変換器19を駆動するパルス駆動器32と、パルス駆動器31
および32の相互のパルス発生タイミング差を制御する遅
延制御部33とから構成されている。第4図(a)におい
て、21はプローブ波用超音波変換器19の出力を増幅する
受信回路、22は受信回路21の出力に信号処理を行う信号
処理部であり、この信号処理部22は第4図(c)に示すよ
うに構成されている。第4図(c)において、41は受信回
路21から出力された受信信号をデジタル信号に変換する
A/D変換部、42はA/D変換部41の出力を記憶するバ
ッファメモリ部、43はバッファメモリ部42の書き込みア
ドレスを発生させる書き込みアドレス発生部、44はバッ
ファメモリ部42の読み出しアドレスを発生させる読み出
しアドレス発生部44、45、46はそれぞれバッファメモリ
部42の出力を任意の区画抜き出すウィンドウ特性部、4
7、48はそれぞれウィンドウ特性部45、46で抜き出され
たデータより受信信号のスペクトルを算出する周波数分
析部、49は周波数分析部47、48の出力により被検体Aの
音響特性を演算する演算部である。第4図(a)におい
て、23は制御駆動部20の遅延制御部33(第4図(b)参
照)、信号処理部22のA/D変換部41(第4図(c)参
照)、後述の主制御部25へクロックを供給するクロック
源、24は信号処理部22の出力を表示する表示部、25はシ
ステム全体の制御を行う主制御部である。
ないし(c)は本発明の第4の実施例における超音波計測
装置を示し、第4図(a)は全体のブロック図、第4図(b)
は制御駆動部のブロック図、第4図(c)は信号処理部の
ブロック図である。第4図(a)において、18は低周波帯
の第1の超音波パルスであるポンプ波パルスを送出する
ポンプ波用超音波変換器、19は第1の超音波パルスより
周波数の高い第2の超音波パルスであるプローブ波パル
スを送出するプローブ波用超音波変換器、20はポンプ波
用超音波変換器18およびプローブ波用超音波変換器19を
位相制御駆動する制御駆動部であり、この制御駆動部20
は第4図(b)から明らかなようにポンプ波用超音波変換
器18を駆動するパルス駆動器31と、プローブ波用超音波
変換器19を駆動するパルス駆動器32と、パルス駆動器31
および32の相互のパルス発生タイミング差を制御する遅
延制御部33とから構成されている。第4図(a)におい
て、21はプローブ波用超音波変換器19の出力を増幅する
受信回路、22は受信回路21の出力に信号処理を行う信号
処理部であり、この信号処理部22は第4図(c)に示すよ
うに構成されている。第4図(c)において、41は受信回
路21から出力された受信信号をデジタル信号に変換する
A/D変換部、42はA/D変換部41の出力を記憶するバ
ッファメモリ部、43はバッファメモリ部42の書き込みア
ドレスを発生させる書き込みアドレス発生部、44はバッ
ファメモリ部42の読み出しアドレスを発生させる読み出
しアドレス発生部44、45、46はそれぞれバッファメモリ
部42の出力を任意の区画抜き出すウィンドウ特性部、4
7、48はそれぞれウィンドウ特性部45、46で抜き出され
たデータより受信信号のスペクトルを算出する周波数分
析部、49は周波数分析部47、48の出力により被検体Aの
音響特性を演算する演算部である。第4図(a)におい
て、23は制御駆動部20の遅延制御部33(第4図(b)参
照)、信号処理部22のA/D変換部41(第4図(c)参
照)、後述の主制御部25へクロックを供給するクロック
源、24は信号処理部22の出力を表示する表示部、25はシ
ステム全体の制御を行う主制御部である。
以上の構成において、以下、その動作について説明す
る。
る。
まず、プローブ波パルスとポンプ波パルスの関係につい
て説明する。プローブ波パルスの一例を第5図(a)に示
し、ポンプ波パルスの一例を第5図(b)に示す。第5図
(c)はプローブ波パルスとポンプ波パルスが重畳された
例を示す。これらの波形は制御駆動部20の遅延制御部33
の制御によるものである。ポンプ波パルスの中心周波数
は、例えば0.3MHz、プローブ波パルスの中心周波数は、
例えば3MHzであり、これらの中心周波数は大幅に異なる
値が選ばれている。第5図(c)においては、プローブ波
パルスの波形の重心はポンプ波パルスの粒子速度が零近
傍で、かつその粒子加速度が正のピークを示すタイミン
グで重畳されている。ポンプ波パルスの波長をT、プロ
ーブ波パルスのパルス長をtとしたとき、次の(1)式の
関係に設定するのが望ましい。
て説明する。プローブ波パルスの一例を第5図(a)に示
し、ポンプ波パルスの一例を第5図(b)に示す。第5図
(c)はプローブ波パルスとポンプ波パルスが重畳された
例を示す。これらの波形は制御駆動部20の遅延制御部33
の制御によるものである。ポンプ波パルスの中心周波数
は、例えば0.3MHz、プローブ波パルスの中心周波数は、
例えば3MHzであり、これらの中心周波数は大幅に異なる
値が選ばれている。第5図(c)においては、プローブ波
パルスの波形の重心はポンプ波パルスの粒子速度が零近
傍で、かつその粒子加速度が正のピークを示すタイミン
グで重畳されている。ポンプ波パルスの波長をT、プロ
ーブ波パルスのパルス長をtとしたとき、次の(1)式の
関係に設定するのが望ましい。
2t<T……(1) 上記(1)式の関係を満たすと、プローブ波パルスの変調
特性の解析を容易に行うことができる。
特性の解析を容易に行うことができる。
次に第5図に示した各パルスが被検体A内を伝搬する様
子を詳しく説明する。通常の超音波診断装置で用いられ
る程度のピーク超音波出力レベルの場合でも、伝搬の非
線形現象により超音波の波形が歪むことが知られてい
る。その原因は次の(2)式により簡単に説明することが
できる。
子を詳しく説明する。通常の超音波診断装置で用いられ
る程度のピーク超音波出力レベルの場合でも、伝搬の非
線形現象により超音波の波形が歪むことが知られてい
る。その原因は次の(2)式により簡単に説明することが
できる。
△C=(1+B/2・A)・U……(2) ここで、Uは音波の粒子速度、B/Aは媒体の音響的非
線形パラメータ、△Cは非線形効果に基づく音速の変化
分である。上記(2)式は、音波の粒子速度Uの方向が音
波進行方向と一致するとき、音速変化△Cは増大し、逆
方向の場合には減少し、結果的に音波の波形が歪むこと
を示す。この伝搬の非線形現象が超音波パルスの波形に
与える影響を第6図に示す。第6図(a)は伝搬に伴うポ
ンプ波パルスの歪を示したものであり、第6図(b)は第
6図(a)に示すポンプ波パルスの歪により重畳されたプ
ローブ波パルスの中心周波数が伝搬に伴い圧縮され、高
周波側へシフトしている状態を示したものである。ま
た、逆にポンプ波パルスの粒子加速度が負のピーク位置
にプローブ波パルスを重畳させた場合には、伝搬に伴い
伸張され、低周波側にシフトする。したがって、中心周
波数が高周波側にシフトする場合と低周波側にシフトす
る場合の差をとれば、非線形効果に基づく中心周波数の
変化量を大きく得ることが可能になる。ポンプ波パルス
の粒子加速度が正のピーク位置にプローブ波パルスを重
畳させたときのプローブ波パルスの受信信号のスペクト
ルをP+(ω)、ポンプ波パルスの粒子加速度が負のピ
ーク位置にプローブ波パルスを重畳させたときのプロー
ブ波パルスの受信信号のスペクトルをP−(ω)とする
と、P±(ω)は次の(3)式で表わされる。
線形パラメータ、△Cは非線形効果に基づく音速の変化
分である。上記(2)式は、音波の粒子速度Uの方向が音
波進行方向と一致するとき、音速変化△Cは増大し、逆
方向の場合には減少し、結果的に音波の波形が歪むこと
を示す。この伝搬の非線形現象が超音波パルスの波形に
与える影響を第6図に示す。第6図(a)は伝搬に伴うポ
ンプ波パルスの歪を示したものであり、第6図(b)は第
6図(a)に示すポンプ波パルスの歪により重畳されたプ
ローブ波パルスの中心周波数が伝搬に伴い圧縮され、高
周波側へシフトしている状態を示したものである。ま
た、逆にポンプ波パルスの粒子加速度が負のピーク位置
にプローブ波パルスを重畳させた場合には、伝搬に伴い
伸張され、低周波側にシフトする。したがって、中心周
波数が高周波側にシフトする場合と低周波側にシフトす
る場合の差をとれば、非線形効果に基づく中心周波数の
変化量を大きく得ることが可能になる。ポンプ波パルス
の粒子加速度が正のピーク位置にプローブ波パルスを重
畳させたときのプローブ波パルスの受信信号のスペクト
ルをP+(ω)、ポンプ波パルスの粒子加速度が負のピ
ーク位置にプローブ波パルスを重畳させたときのプロー
ブ波パルスの受信信号のスペクトルをP−(ω)とする
と、P±(ω)は次の(3)式で表わされる。
P±(ω)=H±(ω):S±(ω)・G(ω)・T(ω)……
(3) 上記(3)式において、H±(ω)は変調されたプローブ波
パルス、S±(ω)は被検体A内の散乱特性、G(ω)は散
乱後の伝搬に伴う減衰特性、T(ω)はプローブ波パルス
用振動子の特性である。変調された2つの受信信号のス
ペクトルP±(ω)の比をとると、 次の(4)式となる。
(3) 上記(3)式において、H±(ω)は変調されたプローブ波
パルス、S±(ω)は被検体A内の散乱特性、G(ω)は散
乱後の伝搬に伴う減衰特性、T(ω)はプローブ波パルス
用振動子の特性である。変調された2つの受信信号のス
ペクトルP±(ω)の比をとると、 次の(4)式となる。
変調されたプローブ波パルスのスペクトルH±(ω)に対
し、被検体A内の散乱特性S±(ω)があまり変化しない
とするならば、上記(4)式は次の(5)式のように近似され
る。
し、被検体A内の散乱特性S±(ω)があまり変化しない
とするならば、上記(4)式は次の(5)式のように近似され
る。
上記(5)式において、P0(ω)の比は、被検体A内の複
雑な散乱特性S±(ω)を含まず、単純に変調されたプロ
ーブ波パルスのスペクトルH±(ω)の比となる。このス
ペクトルH±(ω)の比が1となる周波数をfxとする
と、この周波数fxは被検体の減衰特性や非線形パラメ
ータB/A値等に依存する値であることが解析的に求め
られる。以上より変調されたプローブ波パルスの受信信
号スペクトルP±(ω)の比が1となる周波数を算出する
ことにより、被検体A内の音響的性質を求めることがで
きる。以上のようにポンプ波パルスの粒子加速度の正、
負のピーク位置にプローブ波パルスを重畳した時の受信
信号のスペクトルを比較することで被検体A内の音響特
性を得ることができ、この計測を被検体内の複数の部位
で行うことにより音響特性の分布を求めることができ
る。
雑な散乱特性S±(ω)を含まず、単純に変調されたプロ
ーブ波パルスのスペクトルH±(ω)の比となる。このス
ペクトルH±(ω)の比が1となる周波数をfxとする
と、この周波数fxは被検体の減衰特性や非線形パラメ
ータB/A値等に依存する値であることが解析的に求め
られる。以上より変調されたプローブ波パルスの受信信
号スペクトルP±(ω)の比が1となる周波数を算出する
ことにより、被検体A内の音響的性質を求めることがで
きる。以上のようにポンプ波パルスの粒子加速度の正、
負のピーク位置にプローブ波パルスを重畳した時の受信
信号のスペクトルを比較することで被検体A内の音響特
性を得ることができ、この計測を被検体内の複数の部位
で行うことにより音響特性の分布を求めることができ
る。
次に、これらの受信信号が処理される過程を説明する。
第4図(a)に示すポンプ波用超音波変換器18およびプロ
ーブ波用超音波変換器19は、被検体Aに対し機械的に、
または電子的に扇形走査される。そして、まず、第4図
(b)に示す制御駆動部20の遅延制御部33は、プローブ波
パルスがポンプ波パルスの粒子加速度正のピーク位置に
重畳されるようにパルス駆動器31および32を制御する。
ポンプ波用超音波変換器18とプローブ波用超音波変換器
19から送出されたポンプ波パルスとプローブ波パルス
は、被検体Aの中を伝搬しながら、プローブ波パルスは
非線形現象に基づく伝搬歪を生じて変調される。同時に
被検体A内の音響的質の変化に対応して次々に散乱さ
れ、その一部はプローブ波用超音波変換器19に到達し、
受信信号に変換される。この受信信号は受信回路21でS
/N良く増幅された後、信号処理部22に入力される。受
信回路21の出力は、第4図(c)に示す信号処理部22内の
A/D変換部41でデジタル信号に変換され、バッファメ
モリ部42における書き込みアドレス発生部43の示す位置
に記憶される。次に、第4図(b)に示す遅延制御部33
は、プローブ波パルスがポンプ波パルスの粒子加速度負
のピーク位置に重畳されるようにパルス駆動器31および
32を制御する。ポンプ波用超音波変換器18とプローブ波
用超音波変換器19から送出されたポンプ波パルスとプロ
ーブ波パルスは、上記と同様に被検体A内の音響的質の
変化に対応して次々に散乱され、その一部はプローブ波
用超音波変換器19に到達し、受信信号に変換され、受信
回路21を通り、信号処理部22のA/D変換部41でデジタ
ル信号に変換され、バッファメモリ部42における書き込
みアドレス発生部43の示す位置に記憶され、この動作が
繰り返される。そして、第4図(c)に示すバッファメモ
リ部42から読み出しアドレス発生部44のアドレスにより
読み出し、ウィンドウ特性部45において、算出しようと
する部位の被検深さに対応するデータを、例えばハミン
グウィンドウで抜き出し、周波数分析部47で受信信号の
スペクトルP+(ω)を高速フーリエ変換手法(FF
T)等により算出する。同様に、バッファメモリ部42の
出力はウィンドウ特性部46で対象とする領域を抜き出
し、周波数分析部48で受信信号のスペクトルP−(ω)を
算出する。周波数分析部47および48で算出された2つの
スペクトルP±(ω)は、演算部49において、P±(ω)の
比が1となる周波数fxを求め、この周波数fxよりそ
れに依存する音響特性を求め、表示部24に転送して表示
する。なお、ポンプ波用超音波変換器18とプローブ波用
超音波変換器19は扇形走査することなく、同一部位にお
いて、プローブ波パルスをポンプ波パルスの粒子加速度
正のピーク位置に重畳させた時の受信信号と、プローブ
波パルスをポンプ波パルスの粒子加速度負のピーク位置
に重畳させた時の受信信号をバッファメモリ部42に記憶
させ、次にポンプ波用超音波変換器18とプローブ波用超
音波変換器19を微少角度移動させ、同様な測定を繰り返
しても良い。また、扇形走査する時にプローブ波パルス
をポンプ波パルスの粒子加速度正のピーク位置に重畳さ
せた時の受信信号だけをバッファメモリ部42に記憶し、
次に同じ位置を扇形走査する時にプローブ波パルスをポ
ンプ波パルスの粒子加速度負のピーク位置に重畳させた
時の受信信号をバッファメモリ部42に記憶させても良
い。更に、超音波変換器18、19は扇形走査に替えてリニ
アに走査させても良い。
第4図(a)に示すポンプ波用超音波変換器18およびプロ
ーブ波用超音波変換器19は、被検体Aに対し機械的に、
または電子的に扇形走査される。そして、まず、第4図
(b)に示す制御駆動部20の遅延制御部33は、プローブ波
パルスがポンプ波パルスの粒子加速度正のピーク位置に
重畳されるようにパルス駆動器31および32を制御する。
ポンプ波用超音波変換器18とプローブ波用超音波変換器
19から送出されたポンプ波パルスとプローブ波パルス
は、被検体Aの中を伝搬しながら、プローブ波パルスは
非線形現象に基づく伝搬歪を生じて変調される。同時に
被検体A内の音響的質の変化に対応して次々に散乱さ
れ、その一部はプローブ波用超音波変換器19に到達し、
受信信号に変換される。この受信信号は受信回路21でS
/N良く増幅された後、信号処理部22に入力される。受
信回路21の出力は、第4図(c)に示す信号処理部22内の
A/D変換部41でデジタル信号に変換され、バッファメ
モリ部42における書き込みアドレス発生部43の示す位置
に記憶される。次に、第4図(b)に示す遅延制御部33
は、プローブ波パルスがポンプ波パルスの粒子加速度負
のピーク位置に重畳されるようにパルス駆動器31および
32を制御する。ポンプ波用超音波変換器18とプローブ波
用超音波変換器19から送出されたポンプ波パルスとプロ
ーブ波パルスは、上記と同様に被検体A内の音響的質の
変化に対応して次々に散乱され、その一部はプローブ波
用超音波変換器19に到達し、受信信号に変換され、受信
回路21を通り、信号処理部22のA/D変換部41でデジタ
ル信号に変換され、バッファメモリ部42における書き込
みアドレス発生部43の示す位置に記憶され、この動作が
繰り返される。そして、第4図(c)に示すバッファメモ
リ部42から読み出しアドレス発生部44のアドレスにより
読み出し、ウィンドウ特性部45において、算出しようと
する部位の被検深さに対応するデータを、例えばハミン
グウィンドウで抜き出し、周波数分析部47で受信信号の
スペクトルP+(ω)を高速フーリエ変換手法(FF
T)等により算出する。同様に、バッファメモリ部42の
出力はウィンドウ特性部46で対象とする領域を抜き出
し、周波数分析部48で受信信号のスペクトルP−(ω)を
算出する。周波数分析部47および48で算出された2つの
スペクトルP±(ω)は、演算部49において、P±(ω)の
比が1となる周波数fxを求め、この周波数fxよりそ
れに依存する音響特性を求め、表示部24に転送して表示
する。なお、ポンプ波用超音波変換器18とプローブ波用
超音波変換器19は扇形走査することなく、同一部位にお
いて、プローブ波パルスをポンプ波パルスの粒子加速度
正のピーク位置に重畳させた時の受信信号と、プローブ
波パルスをポンプ波パルスの粒子加速度負のピーク位置
に重畳させた時の受信信号をバッファメモリ部42に記憶
させ、次にポンプ波用超音波変換器18とプローブ波用超
音波変換器19を微少角度移動させ、同様な測定を繰り返
しても良い。また、扇形走査する時にプローブ波パルス
をポンプ波パルスの粒子加速度正のピーク位置に重畳さ
せた時の受信信号だけをバッファメモリ部42に記憶し、
次に同じ位置を扇形走査する時にプローブ波パルスをポ
ンプ波パルスの粒子加速度負のピーク位置に重畳させた
時の受信信号をバッファメモリ部42に記憶させても良
い。更に、超音波変換器18、19は扇形走査に替えてリニ
アに走査させても良い。
そして、本実施例にあっては、第4図中には図示してい
ないが、上記第1の実施例と同様に比較部8とゲート位
置発生部10を備え、または上記第2の実施例と同様にゲ
ート位置発生部10とパワースペクトル値比較部16を備
え、または上記第3の実施例と同様にゲート位置発生部
10とオペレータ部17を備えたものである。
ないが、上記第1の実施例と同様に比較部8とゲート位
置発生部10を備え、または上記第2の実施例と同様にゲ
ート位置発生部10とパワースペクトル値比較部16を備
え、または上記第3の実施例と同様にゲート位置発生部
10とオペレータ部17を備えたものである。
なお、上記第1ないし第3の実施例における超音波探触
子1は、機械走査型のものを示しているが電子走査型で
あってもよい。
子1は、機械走査型のものを示しているが電子走査型で
あってもよい。
発明の効果 以上述べたように本発明によれば、周波数分析部による
周波数分析結果のパワースペクトル値をパワースペクト
ル値比較部において、設定値と比較し、被検体内からの
反射信号成分だけを抜き出している。このように、被検
体内からの反射信号成分だけを抜き出してその音響特性
を算出するので、精度良く被検体内の音響特性を計測す
ることができ、また操作者に対して誤診を招くようなま
ぎらしい結果を表示しないで済む等、その工業的価値は
大なるものがある。
周波数分析結果のパワースペクトル値をパワースペクト
ル値比較部において、設定値と比較し、被検体内からの
反射信号成分だけを抜き出している。このように、被検
体内からの反射信号成分だけを抜き出してその音響特性
を算出するので、精度良く被検体内の音響特性を計測す
ることができ、また操作者に対して誤診を招くようなま
ぎらしい結果を表示しないで済む等、その工業的価値は
大なるものがある。
第1図は本発明の第1の実施例における超音波計測装置
を示すブロック図、第2図は本発明の第2の実施例にお
ける超音波計測装置を示すブロック図、第3図は本発明
の第3の実施例における超音波計測装置を示すブロック
図、第4図(a)ないし(c)は本発明の第4の実施例におけ
る超音波計測装置を示し、同図(a)は全体のブロック
図、同図(b)は制御駆動部のブロック図、同図(c)は信号
処理部のブロック図、第5図(a)はプローブ波の一例を
示す図、第5図(b)はポンプ波パルスの一例を示す図、
第5図(c)はプローブ波パルスとポンプ波パルスが重畳
された例を示す図、第6図(a)は伝搬に伴うポンプ波パ
ルスの歪を示す図、第6図(b)はポンプ波パルスの歪に
より重畳されたプローブ波パルスの中心周波数が伝搬に
伴い圧縮され、高周波側へシフトしている状態を示す
図、第7図は従来の超音波計測装置を示すブロック図で
ある。 1……超音波探触子、2……超音波変換器、3……伝搬
材、4……駆動モータ、5……パルス駆動器、6……受
信回路部、7……A/D変換部、8……比較部、9……
バッファメモリ部、10……ゲート位置発生部、11……周
波数分析部、12……アドレス発生部、13……X−Y表示
メモリ部、14……TVモニタ、15……タイミング発生
部、16……パワースペクトル値比較部、17……オペレー
タ部、18……ポンプ波用超音波変換器、19……プローブ
波用超音波変換器、20……制御駆動部、21……受信回
路、22……信号処理部、23……クロック源、24……表示
部、25……主制御部、31、32……パルス駆動器、33……
遅延制御部、41……A/D変換部、42……バッファメモ
リ部、43…………書き込みアドレス発生部、44……読み
出しアドレス発生部、45、46……ウィンドウ特性部、4
7、48……周波数分析部、49……演算部、A……被検
体。
を示すブロック図、第2図は本発明の第2の実施例にお
ける超音波計測装置を示すブロック図、第3図は本発明
の第3の実施例における超音波計測装置を示すブロック
図、第4図(a)ないし(c)は本発明の第4の実施例におけ
る超音波計測装置を示し、同図(a)は全体のブロック
図、同図(b)は制御駆動部のブロック図、同図(c)は信号
処理部のブロック図、第5図(a)はプローブ波の一例を
示す図、第5図(b)はポンプ波パルスの一例を示す図、
第5図(c)はプローブ波パルスとポンプ波パルスが重畳
された例を示す図、第6図(a)は伝搬に伴うポンプ波パ
ルスの歪を示す図、第6図(b)はポンプ波パルスの歪に
より重畳されたプローブ波パルスの中心周波数が伝搬に
伴い圧縮され、高周波側へシフトしている状態を示す
図、第7図は従来の超音波計測装置を示すブロック図で
ある。 1……超音波探触子、2……超音波変換器、3……伝搬
材、4……駆動モータ、5……パルス駆動器、6……受
信回路部、7……A/D変換部、8……比較部、9……
バッファメモリ部、10……ゲート位置発生部、11……周
波数分析部、12……アドレス発生部、13……X−Y表示
メモリ部、14……TVモニタ、15……タイミング発生
部、16……パワースペクトル値比較部、17……オペレー
タ部、18……ポンプ波用超音波変換器、19……プローブ
波用超音波変換器、20……制御駆動部、21……受信回
路、22……信号処理部、23……クロック源、24……表示
部、25……主制御部、31、32……パルス駆動器、33……
遅延制御部、41……A/D変換部、42……バッファメモ
リ部、43…………書き込みアドレス発生部、44……読み
出しアドレス発生部、45、46……ウィンドウ特性部、4
7、48……周波数分析部、49……演算部、A……被検
体。
───────────────────────────────────────────────────── フロントページの続き (72)発明者 屋野 勉 神奈川県川崎市多摩区東三田3丁目10番1 号 松下技研株式会社内 (56)参考文献 特開 昭60−4857(JP,A) 特開 昭63−122923(JP,A) 特開 昭63−194644(JP,A) 特開 昭63−240829(JP,A) 実開 昭49−86688(JP,U)
Claims (1)
- 【請求項1】超音波を被検体に送信するとともに、当該
被検体からの反射超音波信号を受信する超音波変換器
と、 前記超音波変換器を内包し当該超音波変換器から送受波
される超音波を伝搬させる伝搬材と前記超音波変換器を
2次元走査させる駆動モータからなる超音波探触子と、 前記超音波変換器を電気的に駆動するパルス駆動器と、 前記超音波変換器を介し、被検体からの反射超音波信号
を増幅する受信回路部と、 前記受信回路部の出力信号をアナログ/デジタル変換す
るアナログ/デジタル変換部と、 前記アナログ/デジタル変換部の出力信号である反射超
音波信号を記憶するバッファメモリ部と、 前記バッファメモリ部に記憶されている反射超音波信号
をパワースペクトラム値として周波数分析する周波数分
析部と、 前記周波数分析部の解析結果であるパワースペクトラム
値を、前記伝搬材と被検体との境界位置を示す予め設定
した値と比較し、当該予め設定した値が前記パワースペ
クトラム値より小さいことを検出すると、当該検出結果
が前記伝搬材における反射超音波信号であったとしてア
ドレス算出部へ検出信号を出力するパワースペクトラム
値比較部と、 前記パワースペクトラム値比較部の検出信号が入力され
ると、前記伝搬材における反射超音波信号であるとして
アドレス算出を中止するとともに、一方、前記パワース
ペクトラム値比較部の検出信号が検出されないと、被検
体からの反射超音波信号であるとしてアドレス算出を行
うことにより、被検体からの反射超音波信号のみを抜き
出してアドレス生成するアドレス算出部と、 前記周波数分析部の出力を、前記アドレス算出部がアド
レス生成したアドレスに基づき、当該アドレスに記憶す
るX−Y表示メモリ部と、 前記X−Y表示メモリ部に記憶されている周波数結果を
表示する表示部と を具備する超音波計測装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP1174370A JPH0614929B2 (ja) | 1989-07-05 | 1989-07-05 | 超音波計測装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP1174370A JPH0614929B2 (ja) | 1989-07-05 | 1989-07-05 | 超音波計測装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0339146A JPH0339146A (ja) | 1991-02-20 |
| JPH0614929B2 true JPH0614929B2 (ja) | 1994-03-02 |
Family
ID=15977429
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP1174370A Expired - Fee Related JPH0614929B2 (ja) | 1989-07-05 | 1989-07-05 | 超音波計測装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH0614929B2 (ja) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61286483A (ja) * | 1985-06-10 | 1986-12-17 | 中央発條株式会社 | 窓の開閉操作装置 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS5213186Y2 (ja) * | 1972-11-14 | 1977-03-24 | ||
| JPS604857A (ja) * | 1983-06-22 | 1985-01-11 | Mitsubishi Electric Corp | 超音波送受信装置 |
| JPS63122923A (ja) * | 1986-11-13 | 1988-05-26 | Agency Of Ind Science & Technol | 超音波測温装置 |
| JPS63194644A (ja) * | 1987-02-09 | 1988-08-11 | 松下電器産業株式会社 | 超音波計測装置 |
| JPS63240829A (ja) * | 1987-03-30 | 1988-10-06 | 工業技術院長 | 超音波測温装置 |
-
1989
- 1989-07-05 JP JP1174370A patent/JPH0614929B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0339146A (ja) | 1991-02-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US4754760A (en) | Ultrasonic pulse temperature determination method and apparatus | |
| JPH0428375B2 (ja) | ||
| US4511984A (en) | Ultrasound diagnostic apparatus | |
| EP0952462A2 (en) | Method and apparatus for improving visualization of biopsy needle in ultrasound imaging | |
| JPS627856B2 (ja) | ||
| JPH06504453A (ja) | 乳房組織の画像を生成するための方法および装置 | |
| JPH0467149B2 (ja) | ||
| JPH043222B2 (ja) | ||
| Kelsey et al. | Applications of ultrasound in speech research | |
| JPS61100237A (ja) | 超音波診断装置 | |
| KR20040066795A (ko) | 신체 조직의 기계적 강도 특성을 나타내는 장치 및 방법 | |
| US9918698B2 (en) | System and method for gradient-based k-space search for shear wave velocity dispersion estimation | |
| JPH0614929B2 (ja) | 超音波計測装置 | |
| JPH0548130B2 (ja) | ||
| JP3280473B2 (ja) | 超音波診断装置 | |
| JPS63194644A (ja) | 超音波計測装置 | |
| JPS6053133A (ja) | 超音波診断装置 | |
| JPH0428376B2 (ja) | ||
| García et al. | Design of a continuous wave blood flow bi-directional Doppler system | |
| JPH0564059B2 (ja) | ||
| JPS6053130A (ja) | 超音波診断装置 | |
| JPH024344A (ja) | 超音波計測装置 | |
| JPS63134949A (ja) | 超音波計測装置 | |
| JPH0249103B2 (ja) | ||
| JPH066124B2 (ja) | 超音波計測装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313115 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| LAPS | Cancellation because of no payment of annual fees |