JPH01242043A - 超音波診断装置 - Google Patents
超音波診断装置Info
- Publication number
- JPH01242043A JPH01242043A JP63068774A JP6877488A JPH01242043A JP H01242043 A JPH01242043 A JP H01242043A JP 63068774 A JP63068774 A JP 63068774A JP 6877488 A JP6877488 A JP 6877488A JP H01242043 A JPH01242043 A JP H01242043A
- Authority
- JP
- Japan
- Prior art keywords
- signals
- circuit
- signal
- weighting
- moving average
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 claims description 6
- 238000001514 detection method Methods 0.000 abstract description 14
- 238000013016 damping Methods 0.000 abstract description 6
- 230000002194 synthesizing effect Effects 0.000 abstract 2
- 238000010586 diagram Methods 0.000 description 5
- 210000004185 liver Anatomy 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 238000003786 synthesis reaction Methods 0.000 description 3
- 210000004204 blood vessel Anatomy 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 210000004087 cornea Anatomy 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000011888 foil Substances 0.000 description 1
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000003387 muscular Effects 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
〔概要〕
媒体から受信した超音波信号の減衰係数を表示する超音
波診断装置に関し、 検波信号Ｙ （ｉ）に含まれるスペキュラ成分のみを取
り除いて正しい減衰係数を幅広く取り出して表示するこ
とを目的とし、 検波信号Ｙ （ｉ）について、所定のウィンドウ幅Ｎで
求めた各深さに於ける平均値を当該ウィンドウの任意の
位置に対応づけた複数の移動平均値＜Ｙ （ｉ） ＞ｍ
を生成し、当該検波信号Ｙ （ｉ）から移動平均値＜Ｙ
（ｉ） ＞ｍをそれぞれ減算して得られた複数の信号
に対して、所定の閾値でスライスしてこれらの所定の閾
値を超える部分を重みＯ゛、超えない部分を重みｌ”と
した複数の重み信号を合成して重み付最小２乗法によっ
て減衰係数を算出し、表示するように構成する。
波診断装置に関し、 検波信号Ｙ （ｉ）に含まれるスペキュラ成分のみを取
り除いて正しい減衰係数を幅広く取り出して表示するこ
とを目的とし、 検波信号Ｙ （ｉ）について、所定のウィンドウ幅Ｎで
求めた各深さに於ける平均値を当該ウィンドウの任意の
位置に対応づけた複数の移動平均値＜Ｙ （ｉ） ＞ｍ
を生成し、当該検波信号Ｙ （ｉ）から移動平均値＜Ｙ
（ｉ） ＞ｍをそれぞれ減算して得られた複数の信号
に対して、所定の閾値でスライスしてこれらの所定の閾
値を超える部分を重みＯ゛、超えない部分を重みｌ”と
した複数の重み信号を合成して重み付最小２乗法によっ
て減衰係数を算出し、表示するように構成する。
本発明は、媒体から受信した超音波信号の減衰係数を表
示する超音波診断装置に関するものである。
示する超音波診断装置に関するものである。
〔従来の技術と発明が解決しようとする問題点〕近年、
生体組織の超音波音響特性を画像化し、音響特性と＆ｌ
ｌ時特性を対応づける試みが行われている。音響特性の
１つに組織の減衰特性がある。
生体組織の超音波音響特性を画像化し、音響特性と＆ｌ
ｌ時特性を対応づける試みが行われている。音響特性の
１つに組織の減衰特性がある。
この減衰特性を推定する手法として、超音＄ＲＦ波の解
析信号のパワスペクトラムのモーメントから減衰特性を
推定する手法が知られている。この手法は、検波信号を
Ｙ （ｘ）とすると、減衰特性β（Ｘ）は、近偵的に下
式（１）によって表される。
析信号のパワスペクトラムのモーメントから減衰特性を
推定する手法が知られている。この手法は、検波信号を
Ｙ （ｘ）とすると、減衰特性β（Ｘ）は、近偵的に下
式（１）によって表される。
β（ｘ）−（ｄｙ （ｘ）／ｄｘｌ ／ （−２（ｃ
）・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・
・ ・ ・ ・ ・ ・ ・ ・（１）この式（１）は
、減衰係数傾斜β（ｘ）が検波信号Ｙ （ｘ）の変化分
ｄｙ （ｘ）／ｄｘと、この時の中心周波数ｆｃとから
推定できることを表している。以下定性的に簡単に説明
する。
）・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・ ・
・ ・ ・ ・ ・ ・ ・ ・（１）この式（１）は
、減衰係数傾斜β（ｘ）が検波信号Ｙ （ｘ）の変化分
ｄｙ （ｘ）／ｄｘと、この時の中心周波数ｆｃとから
推定できることを表している。以下定性的に簡単に説明
する。
第７図ｆａｔに示す検波信号Ｙ （ｘ）を適当な長さの
データ（ウィンドウ幅内のデータ）に対して最１Ｊ１２
乗法を適用し、第７図（ｂｌに示すように傾きを各深さ
Ｘについて算出する。第７図ｔｅｌは、超音波反射波を
検波する前の信号から中心円ａ１故ｆ、を算出したもの
である。これら第７図ｆｂｌ、＋ＣＩを式（１１に代入
し、第７図ｆｄｌに示すような減衰係数傾斜β（ｘ）を
算出する。
データ（ウィンドウ幅内のデータ）に対して最１Ｊ１２
乗法を適用し、第７図（ｂｌに示すように傾きを各深さ
Ｘについて算出する。第７図ｔｅｌは、超音波反射波を
検波する前の信号から中心円ａ１故ｆ、を算出したもの
である。これら第７図ｆｂｌ、＋ＣＩを式（１１に代入
し、第７図ｆｄｌに示すような減衰係数傾斜β（ｘ）を
算出する。
この第７図ｆｄ＋に示す減衰係数傾斜β（ｘ）を求める
際、第８図ｔａｌに示すように媒体からの反射信号中に
スペキュラ成分を含む場合、例えば人の脂肪層と筋肉層
との境界、血管と肝臓実質部との境界、肝臓実質部と横
隅膜との境界などからの強烈な反射（スペクトラム成分
）を含む場合、第７図を用いて説明した手法を用いて減
衰係数傾斜β（Ｘ）を算出すると、第８図ｆｂ）−実線
に示すように、強烈な当該スペクトラム成分の影響を受
は正しく算出されず、第８図ｆｂ）点線を用いて示す真
の求めたい減衰係数傾斜β（Ｘ）から大幅に外れてしま
うという問題点があった。
際、第８図ｔａｌに示すように媒体からの反射信号中に
スペキュラ成分を含む場合、例えば人の脂肪層と筋肉層
との境界、血管と肝臓実質部との境界、肝臓実質部と横
隅膜との境界などからの強烈な反射（スペクトラム成分
）を含む場合、第７図を用いて説明した手法を用いて減
衰係数傾斜β（Ｘ）を算出すると、第８図ｆｂ）−実線
に示すように、強烈な当該スペクトラム成分の影響を受
は正しく算出されず、第８図ｆｂ）点線を用いて示す真
の求めたい減衰係数傾斜β（Ｘ）から大幅に外れてしま
うという問題点があった。
本発明は、検波信号Ｙ（ｉ）に含まれるスペキュラ成分
のみを取り除いて正しい減衰係数を可及的に幅広く取り
出して表示することを目的としている。
のみを取り除いて正しい減衰係数を可及的に幅広く取り
出して表示することを目的としている。
第１図を参照して問題点を解決するだめの手段を説明す
る。
る。
第１図において、ＣＦ Ａ Ｒ＋ （Ｃｏｎｓｔａｎｔ
ＦａｌｓｅＡｌａｒｍ Ｒａｔｅ）回路７、ＣＦ Ａ
Ｒｔ、回路８、ＣＦＡＲ１回路９は、検波信号Ｙ （
ｉ）について所定のウィンドウ幅Ｎで求めた各平均値を
当該検波信号Ｙ （ｉ）のウィンドウ幅の右端、左端、
および中央の近傍に対応づけた移動平均値〈Ｙｌ）〉８
、＜Ｙ （ｉ）＞ｍ 、＜ｙ （ｉ）＞ｃを計算し、検
波信号Ｙ （ｉ）から移動平均値＜Ｙ（ｉ）＞ｍ１、＜
Ｙ （ｉ） ＞ｔ 、＜’／ （ｉ） ＞ｍを夫々減算
した信号Ｙ＋＋ （ｉ） 、Ｙｔ （ｉ） 、Ｙｃ
（ｉ）を生成するものである。
ＦａｌｓｅＡｌａｒｍ Ｒａｔｅ）回路７、ＣＦ Ａ
Ｒｔ、回路８、ＣＦＡＲ１回路９は、検波信号Ｙ （
ｉ）について所定のウィンドウ幅Ｎで求めた各平均値を
当該検波信号Ｙ （ｉ）のウィンドウ幅の右端、左端、
および中央の近傍に対応づけた移動平均値〈Ｙｌ）〉８
、＜Ｙ （ｉ）＞ｍ 、＜ｙ （ｉ）＞ｃを計算し、検
波信号Ｙ （ｉ）から移動平均値＜Ｙ（ｉ）＞ｍ１、＜
Ｙ （ｉ） ＞ｔ 、＜’／ （ｉ） ＞ｍを夫々減算
した信号Ｙ＋＋ （ｉ） 、Ｙｔ （ｉ） 、Ｙｃ
（ｉ）を生成するものである。
ＴＬ、回路ｌ０１ＴＬ、回路１１．ＴＩ−＊回路１２は
、信号Ｙ＊ （ｉ） 、Ｙｔ に） 、ＹＣ（ｉ）
に対し、所定閾値でスライスしてこれらの所定閾値を超
える部分を重み°０”、超えない部分を重み“１″とし
た重み信号を夫々生成するものである。
、信号Ｙ＊ （ｉ） 、Ｙｔ に） 、ＹＣ（ｉ）
に対し、所定閾値でスライスしてこれらの所定閾値を超
える部分を重み°０”、超えない部分を重み“１″とし
た重み信号を夫々生成するものである。
合成回路１３は、生成された夫々の重み信号を合成する
ものでる。
ものでる。
重み付き最小２乗回路１４は、合成回路１３によって合
成された重み信号と、検波信号Ｙ （ｉ）とに基づいて
、重み付き最小２乗をし、検波信号Ｙい）の環さ方向の
傾斜を求めるものである。
成された重み信号と、検波信号Ｙ （ｉ）とに基づいて
、重み付き最小２乗をし、検波信号Ｙい）の環さ方向の
傾斜を求めるものである。
乗除算回路１６は、検波信号Ｙ （ｉ）の傾斜と、中心
周波数ｆｃとに基づいて、減衰係数傾斜β（Ｘ）を算出
するものである。
周波数ｆｃとに基づいて、減衰係数傾斜β（Ｘ）を算出
するものである。
本発明は、第１図に示すように、ＣＦＡＲ，回路７ない
し９が媒体２０から受信した検波１３号Ｙい）について
所定のウィンド、つ幅Ｎで求めた各平均値を当該検波信
号Ｙい）のウィンドウ幅の右端、左端、および中央の近
傍に対応づけた移動平均値＜Ｙ（ｉ）＞え、＜Ｙ （ｉ
） ＞ｔ 、＜Ｙ（ｉ）＞ｃを計算し、検波信号Ｙ （
ｉ）から移動平均値〈Ｙ （ｉ） 〉Ｒ、〈Ｙ （ｉ）
＞Ｌ 、＜Ｙ（ｉ）〉。を夫々減算した信号Ｙ、（ｉ
）、Ｙｔ（ｉ） 、Ｙｃ （ｉ）を生成し、ＴＬ、回
路１０ないし１２が信号Ｙａ （ｉ）％ Ｙｔ （
ｉ） 、Ｙｃ（ｉ）に対し、所定閾値でスライスしてこ
れらの所定閾値を超える部分を重み“０”、超えない部
分を重み“ｌ”とした重み信号を夫々生成し、合成回路
１３が生成された夫々の重み信号を合成し、重み付き最
小２乗回路１４がこの合成された重み信号と、検波信号
Ｙ （ｉ）とに基づいて、重み付き最小２乗をして検波
信号Ｙ（＋）の深さ方向の傾斜を求め、乗除算回路１６
が検波信号Ｙ（ｉ）の傾斜と、中心周波数ｆｃとに基づ
いて、減衰係数傾斜β（ｘ）を算出して表示するように
している。
し９が媒体２０から受信した検波１３号Ｙい）について
所定のウィンド、つ幅Ｎで求めた各平均値を当該検波信
号Ｙい）のウィンドウ幅の右端、左端、および中央の近
傍に対応づけた移動平均値＜Ｙ（ｉ）＞え、＜Ｙ （ｉ
） ＞ｔ 、＜Ｙ（ｉ）＞ｃを計算し、検波信号Ｙ （
ｉ）から移動平均値〈Ｙ （ｉ） 〉Ｒ、〈Ｙ （ｉ）
＞Ｌ 、＜Ｙ（ｉ）〉。を夫々減算した信号Ｙ、（ｉ
）、Ｙｔ（ｉ） 、Ｙｃ （ｉ）を生成し、ＴＬ、回
路１０ないし１２が信号Ｙａ （ｉ）％ Ｙｔ （
ｉ） 、Ｙｃ（ｉ）に対し、所定閾値でスライスしてこ
れらの所定閾値を超える部分を重み“０”、超えない部
分を重み“ｌ”とした重み信号を夫々生成し、合成回路
１３が生成された夫々の重み信号を合成し、重み付き最
小２乗回路１４がこの合成された重み信号と、検波信号
Ｙ （ｉ）とに基づいて、重み付き最小２乗をして検波
信号Ｙ（＋）の深さ方向の傾斜を求め、乗除算回路１６
が検波信号Ｙ（ｉ）の傾斜と、中心周波数ｆｃとに基づ
いて、減衰係数傾斜β（ｘ）を算出して表示するように
している。
従って、検波信号Ｙ （ｉ）に含まれるスペキュラ成分
を削除し、幅広く有意な減衰係数１１Ｊｉ斜β（Ｘ）を
表示することが可能となる。
を削除し、幅広く有意な減衰係数１１Ｊｉ斜β（Ｘ）を
表示することが可能となる。
次に、第１図から第６図を用いて本発明の１実施例の構
成および動作を順次詳細に説明する。
成および動作を順次詳細に説明する。
第１図において、媒体２０は、減衰係数傾斜β（Ｘ）を
求めて表示しようとするもの、例えば人の肝臓などであ
る。
求めて表示しようとするもの、例えば人の肝臓などであ
る。
トランスデユーサ１は、超音波を送受信するものである
。
。
プリアンプ２は、トラスデューサ１が受信した超音波信
号を増幅するものである。
号を増幅するものである。
ＴＧＣアンプ３は、受信した超音波信号の深さ方向に対
応してゲインを変えるものである。β（Ｘ）を推定する
場合には、スイッチ（ＳＷ）２１を使って、ＴＧＣアン
プ３を通さないようにする。
応してゲインを変えるものである。β（Ｘ）を推定する
場合には、スイッチ（ＳＷ）２１を使って、ＴＧＣアン
プ３を通さないようにする。
ログアンプ４は、積分型の増幅器である。
検波回路５は、受信した超音波信号を検波するものであ
る。
る。
ＬＰＦ Ｃローパスフィルタ）６は、高周波成分を遮断
するものである。
するものである。
ＺＣ（ゼロクロス検出）回路１５は、受信した超音波信
号のゼロ交差を計数し、中心周波数を推定するものであ
る。
号のゼロ交差を計数し、中心周波数を推定するものであ
る。
表示用メモリ１７は、表示部１８に表示しようとする減
衰係数傾斜β（Ｘ）の信号（画像）を格納するものであ
る。
衰係数傾斜β（Ｘ）の信号（画像）を格納するものであ
る。
表示部１Ｂは、ＣＩ’？Ｔなどの画像表示装置である。
まず、第２図を用いて本発明の全体の動作を簡単に説明
する。
する。
第２図１ｄｌは、スペキュラ成分を含む検波信号Ｙ（ｘ
）を示す。
）を示す。
第２図（ｂｌは、第２図ｆａｔのスペキュラ成分を含む
検波信号Ｙ （ｘ）の平均区間（ウィンドウ）の平均値
＜Ｙ （Ｘ） ＞を示す。
検波信号Ｙ （ｘ）の平均区間（ウィンドウ）の平均値
＜Ｙ （Ｘ） ＞を示す。
第２図１ｄｌは、第２図（ａｌのＹ （ｘ）から第２図
（ｂｌの＜Ｙ （ｘ） ＞を減算したものであって、第
１図のＣＦＡＲ，回路７ないし９の出力値を示す（後述
する）。
（ｂｌの＜Ｙ （ｘ） ＞を減算したものであって、第
１図のＣＦＡＲ，回路７ないし９の出力値を示す（後述
する）。
第２図１ｄｌは、第２図（Ｃ１のＹｃ （ｘ）について
所定の閾値を超えるものをＷ （Ｘ）＝Ｏ，超えないも
のをｗ （ｘ）＝１とした重み信号ｗ （ｘ）を示す。
所定の閾値を超えるものをＷ （Ｘ）＝Ｏ，超えないも
のをｗ （ｘ）＝１とした重み信号ｗ （ｘ）を示す。
第２図１ｄｌは、減衰係数傾斜β（Ｘ）を示す。これは
、第２図ｆａｔの検波（δ号Ｙ （ｘ）と、第２図１ｄ
ｌの重み信号ｗ （ｘ）とに基づいて、重み付き最小２
乗を計算し、中心周波数ｆｃで正規化することにより、
式（１）に対応する減衰係数傾斜β（Ｘ）を求めたもの
である。
、第２図ｆａｔの検波（δ号Ｙ （ｘ）と、第２図１ｄ
ｌの重み信号ｗ （ｘ）とに基づいて、重み付き最小２
乗を計算し、中心周波数ｆｃで正規化することにより、
式（１）に対応する減衰係数傾斜β（Ｘ）を求めたもの
である。
ここで、第２図（ｄｌの重み信号ｗ （ｘ）の閾値（ス
レッシュホールド値）、および平均区間（ウィンドウ幅
）の取り方によって、減衰係数傾斜β（Ｘ）の様子が変
わる。このため、第１図ＣＦＡＲ回路７などによって、
検波信号Ｙ （ｘ）の平均区間における平均値＜Ｙ （
ｘ） ＞を元の検波信号Ｙ （ｘ）から差し引いたＣＦ
ＡＲ信号ＹＣＦＡＩＩ（Ｘ）を下式によって求める。
レッシュホールド値）、および平均区間（ウィンドウ幅
）の取り方によって、減衰係数傾斜β（Ｘ）の様子が変
わる。このため、第１図ＣＦＡＲ回路７などによって、
検波信号Ｙ （ｘ）の平均区間における平均値＜Ｙ （
ｘ） ＞を元の検波信号Ｙ （ｘ）から差し引いたＣＦ
ＡＲ信号ＹＣＦＡＩＩ（Ｘ）を下式によって求める。
Ｙｃｉｓａ （Ｘ）＝Ｙ （ｘ） ＜Ｙ （Ｘ） ＞
・・（２１また、正負の２つのスレッシュホールド値
を第２図１ｄｌのように設定し、ＣＦＡＲ信号ＹＣＦＡ
ＩＩ（Ｘ）が２つのスレソシニホールド値の間にあると
きは、第２図ｆｄｌに示すように重み信号ｗ （ｘ）
−１、その他のときは重み信号ｗ （ｘ）＝Ｏとする。
・・（２１また、正負の２つのスレッシュホールド値
を第２図１ｄｌのように設定し、ＣＦＡＲ信号ＹＣＦＡ
ＩＩ（Ｘ）が２つのスレソシニホールド値の間にあると
きは、第２図ｆｄｌに示すように重み信号ｗ （ｘ）
−１、その他のときは重み信号ｗ （ｘ）＝Ｏとする。
この際、式（２）によって求めるＣＦＡＲ信号ＹＣＦ□
（ｘ）を、例えば■平均区間の前方位置、■平均区間の
後方位置、■平均区間の中心位置とする３通りについて
算出する。数式を用いて表すと下式となる。
（ｘ）を、例えば■平均区間の前方位置、■平均区間の
後方位置、■平均区間の中心位置とする３通りについて
算出する。数式を用いて表すと下式となる。
ｊ・Ｏ
ｊ・Ｏ
ｊ・０
ここで、Ｎは平均区間（ウィンドウ幅）の個数を表す。
■ないし■の第２項の移動平均値くＹ（ｉ）〉え、＜Ｙ
（ｉ） ＞Ｌ 、＜Ｙ （ｉ） ＞（について、第３
図検波体号Ｙ （ｉ）について夫々法めると、第４図１
８）ないしｔｃ＋のようになる。以下順次詳細に説明す
る。
（ｉ） ＞Ｌ 、＜Ｙ （ｉ） ＞（について、第３
図検波体号Ｙ （ｉ）について夫々法めると、第４図１
８）ないしｔｃ＋のようになる。以下順次詳細に説明す
る。
第３図は、検波信号Ｙ （ｉ）を示す。ここで、スペキ
ュラ成分は、図示のように脂肪や筋肉質の部分、肝臓の
実質部の血管の部分、横隔膜の部分で強い反射波として
検出されるものである。
ュラ成分は、図示のように脂肪や筋肉質の部分、肝臓の
実質部の血管の部分、横隔膜の部分で強い反射波として
検出されるものである。
第４図ｉ８＋は、式（３）を用いて第３図検波信号Ｙ（
ｉ）に対して算出した移動平均値＜Ｙ （ｉ） ＞ｍを
示す。これは、第４図ｆａ１図示平均区間Ｎの第３図検
波信号Ｙ（１）の平均値を、当該平均区間Ｎの１つ右側
の位置の値（移動平均値＜”ｔ’ （ｉ）＞ａ）とした
ものである。
ｉ）に対して算出した移動平均値＜Ｙ （ｉ） ＞ｍを
示す。これは、第４図ｆａ１図示平均区間Ｎの第３図検
波信号Ｙ（１）の平均値を、当該平均区間Ｎの１つ右側
の位置の値（移動平均値＜”ｔ’ （ｉ）＞ａ）とした
ものである。
第４図ｆｂｌは、式（４）を用いて第３図検波信号Ｙ（
ｉ）に対して算出した移動平均値＜Ｙ （ｉ） ＞ｍを
示す、これは、第４図（ｂ１図示平均区間Ｎの第３図検
波信号Ｙ （ｉ）の平均値を、当該平均区間Ｎの１つ左
側の位置の値（移動平均値＜Ｙ （＋）＞Ｌ）としたも
のである。
ｉ）に対して算出した移動平均値＜Ｙ （ｉ） ＞ｍを
示す、これは、第４図（ｂ１図示平均区間Ｎの第３図検
波信号Ｙ （ｉ）の平均値を、当該平均区間Ｎの１つ左
側の位置の値（移動平均値＜Ｙ （＋）＞Ｌ）としたも
のである。
第４図（Ｃ１は、式（５）を用いて第３Ｉ２Ｉ検波信号
Ｙい→に対して算出した移動平均値＜Ｙ （ｉ） ＞ｍ
を示す。これは、第４図ｔｅ１図示平均区間Ｎの第３図
検波信号Ｙ （ｉ）の平均値を、当該平均区間Ｎの中央
の位置の値（移動平均値〈７口）〉１）としたものであ
る。
Ｙい→に対して算出した移動平均値＜Ｙ （ｉ） ＞ｍ
を示す。これは、第４図ｔｅ１図示平均区間Ｎの第３図
検波信号Ｙ （ｉ）の平均値を、当該平均区間Ｎの中央
の位置の値（移動平均値〈７口）〉１）としたものであ
る。
第５図（ａ−１） 、（ｂ−１） 、（ｃ−１）は、第
３図検波信号Ｙ （ｉ）から第４図（ａｔ＜Ｙ （ｉ）
＞＊ 、第４図（ｂｌ＜Ｙ （ｉ） ＞ｍ 、第４図
ｆｃ）＜Ｙ（ｉ）＞。を夫々引算して、式（２）のＣＦ
ＡＲ信号ＹＣＦＡＩを夫々算出したものである。
３図検波信号Ｙ （ｉ）から第４図（ａｔ＜Ｙ （ｉ）
＞＊ 、第４図（ｂｌ＜Ｙ （ｉ） ＞ｍ 、第４図
ｆｃ）＜Ｙ（ｉ）＞。を夫々引算して、式（２）のＣＦ
ＡＲ信号ＹＣＦＡＩを夫々算出したものである。
第５図（ａ−２） 、（ｂ−２） 、（ｃ−２）は、第
５図（ａ−１） 、（ｂ−１） 、（ｃ−１）について
、閾値ＴＵ、、ＴＬ６、閾値Ｔ ｔＪｚ 、Ｔ Ｌｚ
、閾値ＴＵ、 、ＴＬ３を図示のように設定し、これら
閾値の間に挟まれるものを重み“ｌ”、それ以外を重み
“０′として夫々生成したものである。
５図（ａ−１） 、（ｂ−１） 、（ｃ−１）について
、閾値ＴＵ、、ＴＬ６、閾値Ｔ ｔＪｚ 、Ｔ Ｌｚ
、閾値ＴＵ、 、ＴＬ３を図示のように設定し、これら
閾値の間に挟まれるものを重み“ｌ”、それ以外を重み
“０′として夫々生成したものである。
第５図（ａ−３） 、（ｂ−３） 、（ｃ−３）は、第
５図（ａ−２） 、（ｂ−２） 、（ｃ−２）について
、第６図に示すように、重み“０”の部分を広げ、不要
な部分を更に取り除くようにしている。具体的には、例
えば第５図（ａ−２１の重みをある幅りで平均値を求め
、この平均値が所定のスレノソユホールドレヘル例えば
０．９５を超えるものを新たな重み“１′とし、超えな
いものを重み“０”とするは作を行う。同様に、第５図
（ｂ−２） 、（ｃ−２）についても幅を広げる操作を
行う。
５図（ａ−２） 、（ｂ−２） 、（ｃ−２）について
、第６図に示すように、重み“０”の部分を広げ、不要
な部分を更に取り除くようにしている。具体的には、例
えば第５図（ａ−２１の重みをある幅りで平均値を求め
、この平均値が所定のスレノソユホールドレヘル例えば
０．９５を超えるものを新たな重み“１′とし、超えな
いものを重み“０”とするは作を行う。同様に、第５図
（ｂ−２） 、（ｃ−２）についても幅を広げる操作を
行う。
第５図＋ｄ＋は、第５図（ａ−３） ＋ （ｂ−３）
、（ｃ−３）を合成したものである。この合成は、図示
矢印を用いて示すように、減衰係数伸行β（ｘ）を求め
ようとする第３図肝臓の実質部の幅が最大限拡大される
ように、第５図Ｆｄｌの左側の部分は第５図（ｂ−３）
の図示左部分を切り出し、第５図ｆｄｌの中央部分は図
（ｃ−３）の図示中央部分を切り出し、第５閃（ｄｌの
右側の部分は第５１２１（ａ−３）の図示右側の部分を
切り出し、１つの重み信号を生成するようにしている。
、（ｃ−３）を合成したものである。この合成は、図示
矢印を用いて示すように、減衰係数伸行β（ｘ）を求め
ようとする第３図肝臓の実質部の幅が最大限拡大される
ように、第５図Ｆｄｌの左側の部分は第５図（ｂ−３）
の図示左部分を切り出し、第５図ｆｄｌの中央部分は図
（ｃ−３）の図示中央部分を切り出し、第５閃（ｄｌの
右側の部分は第５１２１（ａ−３）の図示右側の部分を
切り出し、１つの重み信号を生成するようにしている。
このため、第５図（Ｃ〜３）に示すように、式（５）に
示す移動平均値のみを用いて算出した場合に比し、図中
に記載したように、 有効データが（ｃ−３）より、こ
の部分だけｆｄ＋の方が広い“こととなる。そして、こ
の第５図（ｄｌに示すように合成した後の重み信号ｗ
（ｘ）と、第３図検波信号Ｙ（ｉ）とに基づいて、重み
付きの最小２乗法番こまって減衰係数傾斜β（Ｘ）を算
出し、表示部１８に表示する。
示す移動平均値のみを用いて算出した場合に比し、図中
に記載したように、 有効データが（ｃ−３）より、こ
の部分だけｆｄ＋の方が広い“こととなる。そして、こ
の第５図（ｄｌに示すように合成した後の重み信号ｗ
（ｘ）と、第３図検波信号Ｙ（ｉ）とに基づいて、重み
付きの最小２乗法番こまって減衰係数傾斜β（Ｘ）を算
出し、表示部１８に表示する。
第６図は、第５図（ａ−２） 、（ｂ−２） 、（ｃ−
２）の重み信号の幅を広げる操作を記載したものである
。
２）の重み信号の幅を広げる操作を記載したものである
。
第６図（ａｌは、例えば第５図（ａ−２）につい０幅り
で走査し、この幅りの平均値（移動平均値）を求めて図
示点線として表示したものである。
で走査し、この幅りの平均値（移動平均値）を求めて図
示点線として表示したものである。
第６図（ｂｌは、第６図ｔａ＋の点線のうち、スレノシ
ュホールドレーヘル例えば０．９５を超える部分を新た
な重み“Ｉ゛とし、それ以外を重み“Ｏ′としたもので
ある。
ュホールドレーヘル例えば０．９５を超える部分を新た
な重み“Ｉ゛とし、それ以外を重み“Ｏ′としたもので
ある。
本構成では、タイプの違う３つのＣＦＡＲ回路を用いて
いるが、ＣＦＡＲ回路の個数としては、複数個以上であ
ればよい。また、移動平均信号〈Ｙ （ｉ） ＞ｍを下
式（６）と表現した場合、ｍとして、ｍ−−Ｎ、−Ｎ／
２、ｌに限らず、任意の整数とするようにしてもよい。
いるが、ＣＦＡＲ回路の個数としては、複数個以上であ
ればよい。また、移動平均信号〈Ｙ （ｉ） ＞ｍを下
式（６）と表現した場合、ｍとして、ｍ−−Ｎ、−Ｎ／
２、ｌに限らず、任意の整数とするようにしてもよい。
以上の操作により、重み信号の微細なパルス状のものが
削除される。
削除される。
以上説明したように、本発明によれば、検波１３号Ｙ（
ｉ）について所定のウィンドウ幅で求めた平均値を例え
ば当該ランインドウ幅の右端、中央、および左端の近傍
に位置づけた移動平均値を求め、これらの移動平均値と
検波信号Ｙ （ｉ）との差を閾値でスライスした重みを
夫々算出し、更にこれらの重みを合成して観察したい領
域に対応する減衰係数傾斜β（Ｘ）を最大限有効に取り
出して表示する構成を採用しているため、検波信号Ｙ（
ｉ）中に含まれるスペキュラ成分を削除し、幅広く有意
へ部分の減衰係数傾斜β（ｘ）情叩を精度良好に画面上
に表示することができる。
ｉ）について所定のウィンドウ幅で求めた平均値を例え
ば当該ランインドウ幅の右端、中央、および左端の近傍
に位置づけた移動平均値を求め、これらの移動平均値と
検波信号Ｙ （ｉ）との差を閾値でスライスした重みを
夫々算出し、更にこれらの重みを合成して観察したい領
域に対応する減衰係数傾斜β（Ｘ）を最大限有効に取り
出して表示する構成を採用しているため、検波信号Ｙ（
ｉ）中に含まれるスペキュラ成分を削除し、幅広く有意
へ部分の減衰係数傾斜β（ｘ）情叩を精度良好に画面上
に表示することができる。
第１図は本発明の１実施例構成図、第２図は本発明の詳
細な説明図、第３図は検波信号例、第４図は移動平均値
例、第５図は重み信号生成説明図、第６図は重み信号整
形説明図、第７図は減衰係数傾斜β（ｘ）の算出説明図
、第８図は従来技術の説明図を示す。 図中、ｌはトランスデユーサ、５は検波回路、７ないし
９はＣ，Ｆ Ａ Ｒ回路、ｌＯないし１２はＴＬ開回路
１３は合成回路、１４は重み付き最小２乗回路、Ｉ５は
ＺＣ（ゼロクロス）回路、】６は乗除算回路、１８は表
示部、２０は媒体を表す。 ｃ 人光明の１笑克傅ｊ構成層 沁 １ 閃 乃動乎刀値例 箔４図 も ７ 図
細な説明図、第３図は検波信号例、第４図は移動平均値
例、第５図は重み信号生成説明図、第６図は重み信号整
形説明図、第７図は減衰係数傾斜β（ｘ）の算出説明図
、第８図は従来技術の説明図を示す。 図中、ｌはトランスデユーサ、５は検波回路、７ないし
９はＣ，Ｆ Ａ Ｒ回路、ｌＯないし１２はＴＬ開回路
１３は合成回路、１４は重み付き最小２乗回路、Ｉ５は
ＺＣ（ゼロクロス）回路、】６は乗除算回路、１８は表
示部、２０は媒体を表す。 ｃ 人光明の１笑克傅ｊ構成層 沁 １ 閃 乃動乎刀値例 箔４図 も ７ 図
Claims (1)
- 【特許請求の範囲】 媒体から受信した超音波信号の減衰係数を表示する超音
波診断装置において、 検波信号Ｙ（ｉ）について、所定のウィンドウ幅Ｎで求
めた各深さに於ける平均値を当該ウィンドウの任意の位
置に対応づけた複数の移動平均値＜Ｙ（ｉ）＞＿ｍを生
成し、当該検波信号Ｙ（ｉ）から移動平均値＜Ｙ（ｉ）
＞＿ｍをそれぞれ減算して得られた複数の信号に対して
、所定の閾値でスライスしてこれらの所定の閾値を超え
る部分を重み“０”、超えない部分を重み“１”とした
複数の重み信号を合成して重み付最小２乗法によって減
衰係数を算出し、表示するように構成したことを特徴と
する超音波診断装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP63068774A JPH01242043A (ja) | 1988-03-23 | 1988-03-23 | 超音波診断装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP63068774A JPH01242043A (ja) | 1988-03-23 | 1988-03-23 | 超音波診断装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH01242043A true JPH01242043A (ja) | 1989-09-27 |
Family
ID=13383412
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP63068774A Pending JPH01242043A (ja) | 1988-03-23 | 1988-03-23 | 超音波診断装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH01242043A (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2570806A4 (en) * | 2010-05-10 | 2015-09-30 | Jfe Steel Corp | METHOD FOR ILLUMINATING A STRUCTURAL FORM TO BE WELDED AND A DEVICE THEREFOR |
| JP2017093913A (ja) * | 2015-11-26 | 2017-06-01 | 東芝メディカルシステムズ株式会社 | 超音波診断装置、信号処理装置及び解析プログラム |
-
1988
- 1988-03-23 JP JP63068774A patent/JPH01242043A/ja active Pending
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2570806A4 (en) * | 2010-05-10 | 2015-09-30 | Jfe Steel Corp | METHOD FOR ILLUMINATING A STRUCTURAL FORM TO BE WELDED AND A DEVICE THEREFOR |
| JP2017093913A (ja) * | 2015-11-26 | 2017-06-01 | 東芝メディカルシステムズ株式会社 | 超音波診断装置、信号処理装置及び解析プログラム |
| US11123044B2 (en) | 2015-11-26 | 2021-09-21 | Canon Medical Systems Corporation | Signal processing device, ultrasonic diagnostic apparatus, and method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Rabben et al. | Ultrasound-based vessel wall tracking: an auto-correlation technique with RF center frequency estimation | |
| US7421377B2 (en) | Method and apparatus for supressing noise in a doppler system | |
| US5544659A (en) | Ultrasonic doppler imager having a reduced hardware adaptive tissue rejection filter arrangement | |
| US6689064B2 (en) | Ultrasound clutter filter | |
| US7894874B2 (en) | Method and apparatus for enhancing the detecting and tracking of moving objects using ultrasound | |
| KR100742466B1 (ko) | 이미징 방법 및 장치와 그를 이용한 혈액 움직임 이미징 시스템 | |
| US4542657A (en) | Time domain technique to determine mean frequency | |
| Gastounioti et al. | Comparison of Kalman-filter-based approaches for block matching in arterial wall motion analysis from B-mode ultrasound | |
| US5653234A (en) | Method and apparatus for adaptive spatial image filtering | |
| JPS6034434A (ja) | 生体類の循環系の血流速の多次元流れマップを生成する装置 | |
| JP5256210B2 (ja) | 超音波画像処理方法および超音波画像処理装置 | |
| Van Leeuwen et al. | Simulation of real-time frequency estimators for pulsed Doppler systems | |
| US5494037A (en) | Ultrasonic doppler imager having a spatially smoothed control signal for an adaptive tissue rejection filter | |
| US8758249B2 (en) | Method and apparatus for tissue doppler imaging | |
| Srinivasan et al. | A zero-crossing strain estimator for elastography | |
| US20060030777A1 (en) | T-statistic method for suppressing artifacts in blood vessel ultrasonic imaging | |
| US8162838B2 (en) | Doppler signal spectrum calculating method and apparatus | |
| JPH07265304A (ja) | 動脈の分節を超音波エコーグラフィにより検出し特徴付ける方法および装置 | |
| US8801616B2 (en) | Ultrasonic diagnostic apparatus and ultrasonic image processing method | |
| Kristoffersen et al. | A time-shared ultrasound Doppler measurement and 2-D imaging system | |
| Hoeks et al. | Biophysical principles of vascular diagnosis | |
| JP3773779B2 (ja) | レーダ信号処理装置 | |
| JPH01242043A (ja) | 超音波診断装置 | |
| JP3061430B2 (ja) | 速度測定装置 | |
| JP3281435B2 (ja) | 超音波ドプラ診断装置 |