JPH0239253B2 - - Google Patents
Info
- Publication number
- JPH0239253B2 JPH0239253B2 JP59146926A JP14692684A JPH0239253B2 JP H0239253 B2 JPH0239253 B2 JP H0239253B2 JP 59146926 A JP59146926 A JP 59146926A JP 14692684 A JP14692684 A JP 14692684A JP H0239253 B2 JPH0239253 B2 JP H0239253B2
- Authority
- JP
- Japan
- Prior art keywords
- gaussian
- signal
- circuit
- spectrum
- echo signal
- 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
- 238000001228 spectrum Methods 0.000 claims description 36
- 230000003111 delayed effect Effects 0.000 claims description 16
- 238000000034 method Methods 0.000 claims description 11
- 239000000523 sample Substances 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000001934 delay Effects 0.000 claims description 3
- 238000002405 diagnostic procedure Methods 0.000 claims 1
- 238000002604 ultrasonography Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000003745 diagnosis Methods 0.000 description 3
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- NIXOWILDQLNWCW-UHFFFAOYSA-N acrylic acid group Chemical group C(C=C)(=O)O NIXOWILDQLNWCW-UHFFFAOYSA-N 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Description
【発明の詳細な説明】
［産業上の利用分野］
本発明は超音波パルスを生体内に放射し、生体
内の組織から反射してくるエコー信号を解析し生
体内の任意の組織の位置あるいは当該組織の超音
波に及ぼす減衰係数を測定し、疾病等の診断を行
なう超音波診断装置に関する。
内の組織から反射してくるエコー信号を解析し生
体内の任意の組織の位置あるいは当該組織の超音
波に及ぼす減衰係数を測定し、疾病等の診断を行
なう超音波診断装置に関する。
［従来技術］
従来の一般的な超音波エコー法を使用した診断
装置は生体内の組織から反射してくるエコー信号
の位置とその信号の振幅情報より主に組織の形態
情報を映像化するものであり、組織自体の固有特
性を詳細に知ることはできなかつた。
装置は生体内の組織から反射してくるエコー信号
の位置とその信号の振幅情報より主に組織の形態
情報を映像化するものであり、組織自体の固有特
性を詳細に知ることはできなかつた。
一方、生体中を伝搬する超音波はその伝搬中の
組織により組織固有の周波数依存性に基づく減衰
を受けることが知られている。すなわち、一般に
生体では高い周波数に対する減衰率が、低い周波
数のそれに比べ大きくなる特性を有する。そし
て、減衰特性α（ｆ）は次式のごとく周波数ｆに
比例することが知られている。
組織により組織固有の周波数依存性に基づく減衰
を受けることが知られている。すなわち、一般に
生体では高い周波数に対する減衰率が、低い周波
数のそれに比べ大きくなる特性を有する。そし
て、減衰特性α（ｆ）は次式のごとく周波数ｆに
比例することが知られている。
α（ｆ）＝βf
ここにβは組織により異なる減衰係数である。
また、この減衰係数は異常組織、正常組織の間で
異なることも報告されており、この減衰係数を正
確に測定することが疾病等の診断に強く望まれて
いた。
また、この減衰係数は異常組織、正常組織の間で
異なることも報告されており、この減衰係数を正
確に測定することが疾病等の診断に強く望まれて
いた。
この減衰係数を求める方法として従来提案され
ている方法の１つは、被測定部位の前後より得ら
れた各々のエコー信号のスペクトラムの差からそ
の間で受けた減衰係数を求める方法である。この
方法の欠点は測定距離内の組織構造に起因して生
ずる細かいスペクトラムが減衰係数の測定を困難
にしている点である。
ている方法の１つは、被測定部位の前後より得ら
れた各々のエコー信号のスペクトラムの差からそ
の間で受けた減衰係数を求める方法である。この
方法の欠点は測定距離内の組織構造に起因して生
ずる細かいスペクトラムが減衰係数の測定を困難
にしている点である。
そこで、直接にはスペクトラムの差を求めず、
各々のスペクトラムの中心周波数を求め、その中
心周波数の差から減衰係数を求める方法が提案さ
れている。この方法によれば、上記の細かい組織
構造によるスペクトラムの影響は受けにくい利点
があるが、中心周波数の差の値を減衰係数値に変
換する必要があり、この変換は探触子のスペクト
ラムに関係するため、一般には困難となる。とこ
ろが、そのスペクトラムがガウス型のスペクトラ
ムであると仮定すると、数式的な処理が著しく簡
単となり、その応用が期待されていた。しかしな
がら、現実にはこのようなガウス型のスペクトラ
ムを有する探触子を作ることは極めて困難であつ
た。
各々のスペクトラムの中心周波数を求め、その中
心周波数の差から減衰係数を求める方法が提案さ
れている。この方法によれば、上記の細かい組織
構造によるスペクトラムの影響は受けにくい利点
があるが、中心周波数の差の値を減衰係数値に変
換する必要があり、この変換は探触子のスペクト
ラムに関係するため、一般には困難となる。とこ
ろが、そのスペクトラムがガウス型のスペクトラ
ムであると仮定すると、数式的な処理が著しく簡
単となり、その応用が期待されていた。しかしな
がら、現実にはこのようなガウス型のスペクトラ
ムを有する探触子を作ることは極めて困難であつ
た。
［発明の目的］
本発明は上記従来の課題に鑑みなされたもので
あり、その目的は、通常の非ガウス型スペクトラ
ムを有する探触子を用いながらガウス型のスペク
トラムと同様の極めて簡単な演算処理によつて生
体組織の減衰係数を正確に測定し、これによつて
組織診断を容易に行うことのできる改良された超
音波診断装置を提供することにある。
あり、その目的は、通常の非ガウス型スペクトラ
ムを有する探触子を用いながらガウス型のスペク
トラムと同様の極めて簡単な演算処理によつて生
体組織の減衰係数を正確に測定し、これによつて
組織診断を容易に行うことのできる改良された超
音波診断装置を提供することにある。
［発明の構成］
上記目的を達成するために、本発明は、非ガウ
ス型超音波振動子から超音波パルスを生体内に放
射し生体内組織からのエコー信号をデイコンボリ
ユーシヨン処理で探触子固有のスペクトラムを除
去して解析し所望組織の減衰係数を測定する超音
波診断装置において、前記エコー信号を所定の異
なる遅延時間で遅延処理し複数の遅延信号として
出力する遅延回路と、予め定められ前記各遅延信
号に対応させた所定の重み付け信号と該各遅延信
号とを乗算して該遅延信号にガウス型の重み付け
処理を行なうガウス重み付け回路と、前記重み付
け処理された各遅延信号を加算処理してガウス型
パルス波形信号を演算算出する加算回路と、から
成り該エコー信号のガウス型スペクトラムを得る
デイコンボリユーシヨン処理回路を有することを
特徴とする。
ス型超音波振動子から超音波パルスを生体内に放
射し生体内組織からのエコー信号をデイコンボリ
ユーシヨン処理で探触子固有のスペクトラムを除
去して解析し所望組織の減衰係数を測定する超音
波診断装置において、前記エコー信号を所定の異
なる遅延時間で遅延処理し複数の遅延信号として
出力する遅延回路と、予め定められ前記各遅延信
号に対応させた所定の重み付け信号と該各遅延信
号とを乗算して該遅延信号にガウス型の重み付け
処理を行なうガウス重み付け回路と、前記重み付
け処理された各遅延信号を加算処理してガウス型
パルス波形信号を演算算出する加算回路と、から
成り該エコー信号のガウス型スペクトラムを得る
デイコンボリユーシヨン処理回路を有することを
特徴とする。
これにより、前記エコー信号のスペクトラムを
ガウス型のスペクトラムに変換し該ガウス型のス
ペクトラムの中心周波数の移動量により減衰係数
を正確にかつ容易に測定することができる。
ガウス型のスペクトラムに変換し該ガウス型のス
ペクトラムの中心周波数の移動量により減衰係数
を正確にかつ容易に測定することができる。
［実施例の説明］
以下図面に基づいて本発明を詳細に説明する。
一般に超音波エコー法では、第１図に示される
ように、振動子１０から放射され音響媒体１２中
を伝搬するパルス波Ｓ（ｔ）は音響インピーダン
スの異なる境界面で反射し、振動子１０にて受波
される。被検体の被観測区間での各境界面に対応
する反射系列をｒ（ｔ）とし、またその区間内で
パルス波Ｓ（ｔ）が変化しないとすると、得られ
るエコー信号ｅ（ｔ）は ｅ（ｔ）＝ｒ（ｔ）*Ｓ（ｔ） …(1) と表わされる。ただし＊はコンボリユーシヨンを
表わす。そして、これをフーリエ変換した後の各
関数を大文字で表わすと(1)式は Ｅ（ｆ）＝Ｒ（ｆ）・Ｓ（ｆ） …(2) と積の形で表わされ、故にＲ（ｆ）は Ｒ（ｆ）＝Ｅ（ｆ）／Ｓ（ｆ） …(3) として求めることができる。これは一般にデイコ
ンボリユーシヨンと呼ばれ、公知である。これに
ガウス型のスペクトラムＧ（ｆ）を掛けて EG（ｆ）＝Ｒ（ｆ）・Ｇ（ｆ） …(4) (4)式のEG（ｆ）を求める。｜Ｒ（ｆ）｜が｜Ｇ
（ｆ）｜より十分小さいとき、Ｇ（ｆ）の中心周波
数はほぼEG（ｆ）のものと一致する。
ように、振動子１０から放射され音響媒体１２中
を伝搬するパルス波Ｓ（ｔ）は音響インピーダン
スの異なる境界面で反射し、振動子１０にて受波
される。被検体の被観測区間での各境界面に対応
する反射系列をｒ（ｔ）とし、またその区間内で
パルス波Ｓ（ｔ）が変化しないとすると、得られ
るエコー信号ｅ（ｔ）は ｅ（ｔ）＝ｒ（ｔ）*Ｓ（ｔ） …(1) と表わされる。ただし＊はコンボリユーシヨンを
表わす。そして、これをフーリエ変換した後の各
関数を大文字で表わすと(1)式は Ｅ（ｆ）＝Ｒ（ｆ）・Ｓ（ｆ） …(2) と積の形で表わされ、故にＲ（ｆ）は Ｒ（ｆ）＝Ｅ（ｆ）／Ｓ（ｆ） …(3) として求めることができる。これは一般にデイコ
ンボリユーシヨンと呼ばれ、公知である。これに
ガウス型のスペクトラムＧ（ｆ）を掛けて EG（ｆ）＝Ｒ（ｆ）・Ｇ（ｆ） …(4) (4)式のEG（ｆ）を求める。｜Ｒ（ｆ）｜が｜Ｇ
（ｆ）｜より十分小さいとき、Ｇ（ｆ）の中心周波
数はほぼEG（ｆ）のものと一致する。
次に、生体のような減衰のある媒質中からのエ
コー信号のスペクトラムをEF（ｆ）と表わすと EF（ｆ）＝Ｒ（ｆ）・Ｓ（ｆ）e-〓fl …(5) のように表わせる。ただし、ｌは超音波パルスが
反射体から戻つてくるまでの往復距離とする。(5)
式を(3)式のＥ（ｆ）に代入し、更に(4)式のごとく
ガウス型スペクトラムＧ（ｆ）を掛けて(6)式を得
る。
コー信号のスペクトラムをEF（ｆ）と表わすと EF（ｆ）＝Ｒ（ｆ）・Ｓ（ｆ）e-〓fl …(5) のように表わせる。ただし、ｌは超音波パルスが
反射体から戻つてくるまでの往復距離とする。(5)
式を(3)式のＥ（ｆ）に代入し、更に(4)式のごとく
ガウス型スペクトラムＧ（ｆ）を掛けて(6)式を得
る。
EFG（ｆ）＝Ｒ（ｆ）e-〓flＧ（ｆ） …(6)
Ｇ（ｆ）はガウス型であるため(7)式のごとく仮
定する。
定する。
Ｇ（ｆ）＝exp｛（ｆ−fo）2／2σ2｝ …(7)
ただし、foは中心周波数、σは統計学で使用さ
れているガウス分布の標準偏差に相当するものと
する。(7)式を(6)式に代入し(8)式を得る。
れているガウス分布の標準偏差に相当するものと
する。(7)式を(6)式に代入し(8)式を得る。
EFG（ｆ）
＝KR（ｆ）exp［｛ｆ−fo＋Δf）｝2／2σ］ …(8)
ただしＫは定数とする。またΔfは
Δf＝αflσ2 …(9)
とした。(4)式のＧ（ｆ）に(7)式を代入すると次式
が得られる。
が得られる。
EG（ｆ）＝Ｒ（ｆ）exp｛（ｆ−fo）2／2σ2｝…(10)
(10)式を(8)式と比較すると、減衰特性の影響を受
けたEFG（ｆ）も、ガウス型スペクトラムを維持し
ていることが分かる。ただし、スペクトラムの中
心周波数がΔf移動している点が異なつている。
つまり、この中心周波数の移動分を求めることに
より(9)式を用いて減衰係数αを得ることができ
る。
けたEFG（ｆ）も、ガウス型スペクトラムを維持し
ていることが分かる。ただし、スペクトラムの中
心周波数がΔf移動している点が異なつている。
つまり、この中心周波数の移動分を求めることに
より(9)式を用いて減衰係数αを得ることができ
る。
本発明において特徴的なことは、非ガウス型の
探触子から送受したエコー信号に対しデイコンボ
リユーシヨン処理によりガウス型の探触子から得
たエコー信号と等価の信号を得ることができるこ
とである。
探触子から送受したエコー信号に対しデイコンボ
リユーシヨン処理によりガウス型の探触子から得
たエコー信号と等価の信号を得ることができるこ
とである。
次に、本発明のデイコンボリユーシヨン処理を
コンピユータで行つた例を示す。
コンピユータで行つた例を示す。
第２図は水中に置かれたアクリル板の表面より
得られたエコー波形であり、第３図はそのスペク
トラムである。第２図の波形から明らかなよう
に、このエコーは時間軸に沿つて後方へ減衰して
行く形を示し、非ガウス型のパルス波形である。
一方、デイコンボリユーシヨン処理を行いガウス
型の重み付けを施した第４図の波形は、そのエン
ベロープがガウス型で重み付けされたガウス型パ
ルスを示している。
得られたエコー波形であり、第３図はそのスペク
トラムである。第２図の波形から明らかなよう
に、このエコーは時間軸に沿つて後方へ減衰して
行く形を示し、非ガウス型のパルス波形である。
一方、デイコンボリユーシヨン処理を行いガウス
型の重み付けを施した第４図の波形は、そのエン
ベロープがガウス型で重み付けされたガウス型パ
ルスを示している。
第５図は第４図のスペクトラムであり、ガウス
型パルスのスペクトラムも同様にガウス型とな
り、第４図は実際の波形を処理して得られたデー
タを基としているため、形に少し凹凸があるもの
の、ガウス型スペクトラムを示していることが理
解される。つまり、これはガウス型パルスを送受
信して得られる波形とスペクトラムを示し、現実
には困難なガウス型パルスによるエコー信号を得
ることができたことを示す。
型パルスのスペクトラムも同様にガウス型とな
り、第４図は実際の波形を処理して得られたデー
タを基としているため、形に少し凹凸があるもの
の、ガウス型スペクトラムを示していることが理
解される。つまり、これはガウス型パルスを送受
信して得られる波形とスペクトラムを示し、現実
には困難なガウス型パルスによるエコー信号を得
ることができたことを示す。
第６図には、本発明の超音波診断装置の好適な
実施例が示されている。
実施例が示されている。
実施例の装置は、振動子１０のフオーカシング
調整を行うフオーカス回路及び振動子１０の電子
切替走査を行う切替回路を含むスキヤナ１６を有
し、オシレータ回路１８からスキヤナ１６に超音
波パルスの送信繰返し周期を定めるトリガパルス
が入力されると、このスキヤナ１６は該トリガパ
ルスの入力に同期して振動子１０を励振制御し、
被検体１２に向け超音波パルスＳ（ｔ）の放射を
行う。
調整を行うフオーカス回路及び振動子１０の電子
切替走査を行う切替回路を含むスキヤナ１６を有
し、オシレータ回路１８からスキヤナ１６に超音
波パルスの送信繰返し周期を定めるトリガパルス
が入力されると、このスキヤナ１６は該トリガパ
ルスの入力に同期して振動子１０を励振制御し、
被検体１２に向け超音波パルスＳ（ｔ）の放射を
行う。
このようにして放射された超音波パルスＳ（ｔ）
は、被検体１２内部を伝搬する途中で音響インピ
ーダンスの異なる各境界面で次々と反射される。
そして、その反射エコーｅ（ｔ）は振動子１０で
受波され、ここで電気信号に変換され、更にスキ
ヤナ１６内に設けられたフオーカス回路を介して
本発明の特徴的構成要素であるデイコンボリユー
シヨン回路及びガウス重み付け回路２０に入力さ
れる。
は、被検体１２内部を伝搬する途中で音響インピ
ーダンスの異なる各境界面で次々と反射される。
そして、その反射エコーｅ（ｔ）は振動子１０で
受波され、ここで電気信号に変換され、更にスキ
ヤナ１６内に設けられたフオーカス回路を介して
本発明の特徴的構成要素であるデイコンボリユー
シヨン回路及びガウス重み付け回路２０に入力さ
れる。
本発明の装置において、この回路２０は、エコ
ー信号を遅延処理しこれを一定の間隔で遅延され
た複数の遅延信号として出力する遅延回路と、各
遅延回路に所定の重み付け処理を行なう重み付け
回路と、重み付けされた各遅延信号を加熱出力す
る加算器とを含み、これら回路を用いて入力され
たエコー信号ｅ（ｔ）をデイコンボリユーシヨン
処理し、更にガウス型の重み付け処理を施すこと
により得られたガウス型エコー信号を増幅器２２
及び中心周波数検出器３０に向け出力している。
ー信号を遅延処理しこれを一定の間隔で遅延され
た複数の遅延信号として出力する遅延回路と、各
遅延回路に所定の重み付け処理を行なう重み付け
回路と、重み付けされた各遅延信号を加熱出力す
る加算器とを含み、これら回路を用いて入力され
たエコー信号ｅ（ｔ）をデイコンボリユーシヨン
処理し、更にガウス型の重み付け処理を施すこと
により得られたガウス型エコー信号を増幅器２２
及び中心周波数検出器３０に向け出力している。
このようにして増幅器２２に入力されたガウス
型信号は、ここで所定の増幅率をもつて増幅され
た後検波器２４を介してメモリ２６に入力され、
ここでビデオ信号として記憶された後CRT２８
上に診断画像として表示されることとなる。
型信号は、ここで所定の増幅率をもつて増幅され
た後検波器２４を介してメモリ２６に入力され、
ここでビデオ信号として記憶された後CRT２８
上に診断画像として表示されることとなる。
一方、中心周波数検出器３０にてガウス型信号
の中心周波数が指示された深さに相当する時点の
ガウス型信号より測定され、減衰係数演算器３２
に送られる。そして演算器３２では減衰のない状
態で測定された中心周波数との差から減衰係数が
計算され、メモリ２６に書き込まれ、CRTに表
示される。以下に本発明のデイコンボリユーシヨ
ン回路及びガウス重み付け回路２０の好適な実施
例を説明する。
の中心周波数が指示された深さに相当する時点の
ガウス型信号より測定され、減衰係数演算器３２
に送られる。そして演算器３２では減衰のない状
態で測定された中心周波数との差から減衰係数が
計算され、メモリ２６に書き込まれ、CRTに表
示される。以下に本発明のデイコンボリユーシヨ
ン回路及びガウス重み付け回路２０の好適な実施
例を説明する。
第７図には、非再帰型のデイコンボリユーシヨ
ン処理回路２０の好適な実施例が示されており、
実施例の回路は、スキヤナ１６から出力されるエ
コー信号ｅ（ｔ）をそのまま遅延回路３４に入力
し、エコー信号ｅ（ｔ）の非再帰型デイコンボリ
ユーシヨン処理を行つている。
ン処理回路２０の好適な実施例が示されており、
実施例の回路は、スキヤナ１６から出力されるエ
コー信号ｅ（ｔ）をそのまま遅延回路３４に入力
し、エコー信号ｅ（ｔ）の非再帰型デイコンボリ
ユーシヨン処理を行つている。
ここにおいて、前記遅延回路３４は、入力され
たエコー信号ｅ（ｔ）に所定の遅延処理を施し、
エコー信号ｅ（ｔ）をそれぞれ一定時間間隔τで
遅延した複数の遅延信号ｅ（ｔ−τ）、ｅ（ｔ−2τ）
………ｅ（ｔ−nτ）として重み付け回路３６に向
け出力している。
たエコー信号ｅ（ｔ）に所定の遅延処理を施し、
エコー信号ｅ（ｔ）をそれぞれ一定時間間隔τで
遅延した複数の遅延信号ｅ（ｔ−τ）、ｅ（ｔ−2τ）
………ｅ（ｔ−nτ）として重み付け回路３６に向
け出力している。
重み付け回路３６は、図示しない制御回路から
の信号に基づき各遅延信号に対応した複数の重み
付け信号a1，a2…aoを出力する重み付け器３８
と、これら重み付け信号と遅延回路３４の出力す
る各遅延信号とを乗算する複数の乗算器４０―
１，４０―２…４０―ｎとを含み、遅延回路３４
から出力される遅延信号に所定の重み付け処理を
行つている。
の信号に基づき各遅延信号に対応した複数の重み
付け信号a1，a2…aoを出力する重み付け器３８
と、これら重み付け信号と遅延回路３４の出力す
る各遅延信号とを乗算する複数の乗算器４０―
１，４０―２…４０―ｎとを含み、遅延回路３４
から出力される遅延信号に所定の重み付け処理を
行つている。
そして、このようにして重み付け回路３６から
出力される重み付けされた各遅延信号は加算器４
２にて加算され、デイコンボリユーシヨンされた
画像信号ｒ（ｔ）となる。
出力される重み付けされた各遅延信号は加算器４
２にて加算され、デイコンボリユーシヨンされた
画像信号ｒ（ｔ）となる。
従つて、加算器４２から出力される信号ｒ（ｔ）
は次式 ｒ（ｔ）＝o 〓i=1 aiｅ（ｔ―iτ） …(11) をもつて表わされることになる。
は次式 ｒ（ｔ）＝o 〓i=1 aiｅ（ｔ―iτ） …(11) をもつて表わされることになる。
ここにおいて、前記(11)式の両辺をフーリエ変換
すると、 Ｒ（ｆ）＝∫∞ -∞∞ 〓i=1 aiｅ（ｔ−iτ）e-j2〓ftdt ＝Ｅ（ｆ）∞ 〓i=1 aie-j2〓fi〓 …(12) となり、これはＥ（ｆ）と∞ 〓i=1 aie-j2〓fi〓のコンボリ
ユーシヨンの形を表わしている。従つて、この(12)
式を前記(3)式と比較すれば、 ∞ 〓i=1 aie-j2〓fi〓＝１／Ｓ（ｆ） …(13) の関係が得られ、この(13)式で表わされる１／Ｓ
（ｆ）を逆フーリエ変換することにより、aiを求
めることができる。
すると、 Ｒ（ｆ）＝∫∞ -∞∞ 〓i=1 aiｅ（ｔ−iτ）e-j2〓ftdt ＝Ｅ（ｆ）∞ 〓i=1 aie-j2〓fi〓 …(12) となり、これはＥ（ｆ）と∞ 〓i=1 aie-j2〓fi〓のコンボリ
ユーシヨンの形を表わしている。従つて、この(12)
式を前記(3)式と比較すれば、 ∞ 〓i=1 aie-j2〓fi〓＝１／Ｓ（ｆ） …(13) の関係が得られ、この(13)式で表わされる１／Ｓ
（ｆ）を逆フーリエ変換することにより、aiを求
めることができる。
従つて、１／Ｓ（ｆ）を求めておきこれを逆フ
ーリエ変換することによりデイコンボリユーシヨ
ンに必要なaiを予め求めることができ、このよう
にして求めたaiを重み付け器３８から出力するこ
とにより加算器４２から出力される前記(11)式の信
号ｒ（ｔ）は、エコー信号ｅ（ｔ）をデイコンボリ
ユーシヨン処理したものとすることができる。
ーリエ変換することによりデイコンボリユーシヨ
ンに必要なaiを予め求めることができ、このよう
にして求めたaiを重み付け器３８から出力するこ
とにより加算器４２から出力される前記(11)式の信
号ｒ（ｔ）は、エコー信号ｅ（ｔ）をデイコンボリ
ユーシヨン処理したものとすることができる。
次に、このデイコンボリユーシヨン処理された
信号をガウス型信号に変換する方法について述べ
る。ガウス型信号のスペクトラムをEG（ｆ）とす
ると、これは前述の(4)式のようにガウス型スペク
トラムＧ（ｆ）を掛けて得られる。従つて、 EG（ｆ）＝Ｅ（ｆ）・Ｇ（ｆ）／Ｓ（ｆ） …(14) と表わすことができる。
信号をガウス型信号に変換する方法について述べ
る。ガウス型信号のスペクトラムをEG（ｆ）とす
ると、これは前述の(4)式のようにガウス型スペク
トラムＧ（ｆ）を掛けて得られる。従つて、 EG（ｆ）＝Ｅ（ｆ）・Ｇ（ｆ）／Ｓ（ｆ） …(14) と表わすことができる。
つまり(13)式で重み付け定数aiを求める際に
∞
〓i=1
aie-j2〓fi〓＝Ｇ（ｆ）／Ｓ（ｆ） …(15)
(15)式よりaiを求めることによりデイコンボリ
ユーシヨンとガウス型の重み付けの両方を一度に
行うことが可能となる。
ユーシヨンとガウス型の重み付けの両方を一度に
行うことが可能となる。
［発明の効果］
以上説明したように、本発明によれば、非ガウ
ス型のスペクトラムを有する超音波探触子を用い
ながら、被検体から得られたエコー信号に対して
デイコンボリユーシヨン処理を行ない探触子固有
のスペクトラムを取り除いた後にガウス型のスペ
クトラムへの変換処理を行い、この結果、ガウス
型のスペクトラムを有する探触子と同等のエコー
信号を得ることができ、これによつて、減衰係数
を容易に求めることを可能とし、超音波診断によ
り得られる情報を著しく増大することが可能とな
る。
ス型のスペクトラムを有する超音波探触子を用い
ながら、被検体から得られたエコー信号に対して
デイコンボリユーシヨン処理を行ない探触子固有
のスペクトラムを取り除いた後にガウス型のスペ
クトラムへの変換処理を行い、この結果、ガウス
型のスペクトラムを有する探触子と同等のエコー
信号を得ることができ、これによつて、減衰係数
を容易に求めることを可能とし、超音波診断によ
り得られる情報を著しく増大することが可能とな
る。
第１図は一般的な超音波診断装置における放射
された超音波及び反射エコー信号を示す説明図、
第２，３，４，５図はそれぞれ本発明の原理を示
す説明図、第６図は本発明に係る超音波診断装置
の好適な実施例を示す要部概略図、第７図は本発
明に用いられるデイコンボリユーシヨン回路及び
重み付け回路の好適な実施例を示すブロツク図で
ある。 １０…振動子、１２…被検体、１６…スキヤ
ナ、２０…デイコンボリユーシヨン処理回路、２
４…検波器、２６…メモリ、２８…CRT、３０
…中心周波数検出器、３２…減衰係数演算器、３
６…重み付け回路、３８…重み付け器。
された超音波及び反射エコー信号を示す説明図、
第２，３，４，５図はそれぞれ本発明の原理を示
す説明図、第６図は本発明に係る超音波診断装置
の好適な実施例を示す要部概略図、第７図は本発
明に用いられるデイコンボリユーシヨン回路及び
重み付け回路の好適な実施例を示すブロツク図で
ある。 １０…振動子、１２…被検体、１６…スキヤ
ナ、２０…デイコンボリユーシヨン処理回路、２
４…検波器、２６…メモリ、２８…CRT、３０
…中心周波数検出器、３２…減衰係数演算器、３
６…重み付け回路、３８…重み付け器。
Claims (1)
- 【特許請求の範囲】 １ 非ガウス型超音波振動子から超音波パルスを
生体内に放射し生体内組織からのエコー信号をデ
イコンボリユーシヨン処理で探触子固有のスペク
トラムを除去して解析し所望組織の減衰係数を測
定する超音波診断装置において、 前記エコー信号を所定の異なる遅延時間で遅延
処理し複数の遅延信号として出力する遅延回路
と、 予め定められ前記各遅延信号に対応させた所定
の重み付け信号と該各遅延信号とを乗算して該遅
延信号にガウス型の重み付け処理を行なうガウス
重み付け回路と、前記重み付け処理された各遅延
信号を加算処理してガウス型パルス波形信号を演
算算出する加算回路と、から成るデイコンボリユ
ーシヨン処理回路を有し、 前記エコー信号のスペクトラムをガウス型のス
ペクトラムに変換し該ガウス型のスペクトラムの
中心周波数の移動量により所望組織の減衰係数を
正確に測定することを特徴とする超音波診断装
置。 ２ 特許請求の範囲１記載の装置において、前記
ガウス型パルス波形信号への変換は、前記エコー
信号を直接前記遅延回路に入力し該エコー信号の
非再帰型又は再帰型デイコンボリユーシヨン処理
を行う際に同時に行うことを特徴とする超音波診
断装置。 ３ 特許請求の範囲１、２のいずれかに記載され
た装置において、ガウス型に変換されたエコー信
号の中心周波数の移動量から所望組織の減衰係数
を求め、断層画像と同時に表示することを特徴と
する超音波診断装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP14692684A JPS6125535A (ja) | 1984-07-17 | 1984-07-17 | 超音波診断装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP14692684A JPS6125535A (ja) | 1984-07-17 | 1984-07-17 | 超音波診断装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS6125535A JPS6125535A (ja) | 1986-02-04 |
| JPH0239253B2 true JPH0239253B2 (ja) | 1990-09-04 |
Family
ID=15418695
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP14692684A Granted JPS6125535A (ja) | 1984-07-17 | 1984-07-17 | 超音波診断装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS6125535A (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2537183Y2 (ja) * | 1990-06-08 | 1997-05-28 | 門倉技研工業株式会社 | 鋳造機の射出シリンダ |
| JP4492556B2 (ja) | 2006-02-17 | 2010-06-30 | トヨタ自動車株式会社 | 内燃機関の吸気ポート構造 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS53133481A (en) * | 1977-04-25 | 1978-11-21 | Krautkramer Branson | Improved means in resolving power in real time supersonic shadowing apparatus |
| JPS5749439A (en) * | 1980-07-02 | 1982-03-23 | Philips Nv | Apparatus for inspecting inside structure of object to be inspected by ultrasonic |
| JPS60127457A (ja) * | 1983-12-14 | 1985-07-08 | Fujitsu Ltd | 超音波媒体特性値測定装置 |
| JPS60252258A (ja) * | 1984-05-29 | 1985-12-12 | Fujitsu Ltd | 超音波媒体特性値測定装置 |
-
1984
- 1984-07-17 JP JP14692684A patent/JPS6125535A/ja active Granted
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS53133481A (en) * | 1977-04-25 | 1978-11-21 | Krautkramer Branson | Improved means in resolving power in real time supersonic shadowing apparatus |
| JPS5749439A (en) * | 1980-07-02 | 1982-03-23 | Philips Nv | Apparatus for inspecting inside structure of object to be inspected by ultrasonic |
| JPS60127457A (ja) * | 1983-12-14 | 1985-07-08 | Fujitsu Ltd | 超音波媒体特性値測定装置 |
| JPS60252258A (ja) * | 1984-05-29 | 1985-12-12 | Fujitsu Ltd | 超音波媒体特性値測定装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPS6125535A (ja) | 1986-02-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US4452082A (en) | Ultrasonic measuring method | |
| Walker et al. | A fundamental limit on the performance of correlation based phase correction and flow estimation techniques | |
| US4936308A (en) | Method and apparatus for measuring acoustic characteristics and temperature | |
| US5197019A (en) | Method of measuring distance using ultrasonic waves | |
| US4771786A (en) | Ultrasonic method and apparatus for tissue characterization and imaging of nonlinear parameter | |
| EP0123427B1 (en) | Ultrasonic medium characterization | |
| EP0041403A1 (en) | Method and apparatus for ultrasonic wave measurement of characteristics of the internal structure of an object | |
| JPH0246213B2 (ja) | ||
| US4511984A (en) | Ultrasound diagnostic apparatus | |
| JPS6377437A (ja) | 超音波エコーグラフィによる可動物体検査装置 | |
| JPS60122367A (ja) | 超音波測定方法およびその装置 | |
| EP2903530B1 (en) | Shear wave attenuation from k-space analysis system | |
| CN110840488B (zh) | 一种基于剪切波的成像方法、系统及装置 | |
| JPH0713631B2 (ja) | 媒体の超音波エコーグラフィック検査方法及び装置 | |
| US20070238991A1 (en) | Ultrasound method and apparatus for characterizing and identifying biological tissues | |
| JPS61100237A (ja) | 超音波診断装置 | |
| JP3281435B2 (ja) | 超音波ドプラ診断装置 | |
| JPH0239253B2 (ja) | ||
| JP2669204B2 (ja) | 探査装置 | |
| US9918698B2 (en) | System and method for gradient-based k-space search for shear wave velocity dispersion estimation | |
| Lacefield et al. | Spatial coherence analysis applied to aberration correction using a two-dimensional array system | |
| JP2563656B2 (ja) | 超音波ドプラ映像装置 | |
| US4794546A (en) | Method of and apparatus for scanning objects by means of ultrasound echography | |
| JPH0221259B2 (ja) | ||
| JP2815622B2 (ja) | 超音波診断装置 |