JPH057588A - 超音波ドプラ診断装置 - Google Patents
超音波ドプラ診断装置Info
- Publication number
- JPH057588A JPH057588A JP3161913A JP16191391A JPH057588A JP H057588 A JPH057588 A JP H057588A JP 3161913 A JP3161913 A JP 3161913A JP 16191391 A JP16191391 A JP 16191391A JP H057588 A JPH057588 A JP H057588A
- Authority
- JP
- Japan
- Prior art keywords
- clutter
- signal
- information
- moving speed
- received 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/50—Systems of measurement, based on relative movement of the target
- G01S15/52—Discriminating between fixed and moving objects or between objects moving at different speeds
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8979—Combined Doppler and pulse-echo imaging systems
- G01S15/8981—Discriminating between fixed and moving objects or between objects moving at different speeds, e.g. wall clutter filter
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
(57)【要約】
【目的】本発明は、超音波のドプラ効果を利用して被検
体内の血流情報を得る超音波ドプラ診断装置に関し、非
常に遅い速度をもつ血流であってもそのドプラ偏移周波
数がクラッタ成分のドプラ偏移周波数と異なっている限
り、クラッタ成分の情報のみを選択的に除去する。 【構成】受信信号に基づいてクラッタ成分の移動速度及
び／又はその分散を求めるクラッタ情報検出手段と、こ
の求められた移動速度及び／又は分散に基づいて、受信
信号から、該受信信号が担持するクラッタ成分の情報が
選択的に除去されたクラッタ除去信号を求めるクラッタ
情報除去手段と、このクラッタ除去信号に基づいて血流
情報を求める血流情報検出手段とを備えた。
体内の血流情報を得る超音波ドプラ診断装置に関し、非
常に遅い速度をもつ血流であってもそのドプラ偏移周波
数がクラッタ成分のドプラ偏移周波数と異なっている限
り、クラッタ成分の情報のみを選択的に除去する。 【構成】受信信号に基づいてクラッタ成分の移動速度及
び／又はその分散を求めるクラッタ情報検出手段と、こ
の求められた移動速度及び／又は分散に基づいて、受信
信号から、該受信信号が担持するクラッタ成分の情報が
選択的に除去されたクラッタ除去信号を求めるクラッタ
情報除去手段と、このクラッタ除去信号に基づいて血流
情報を求める血流情報検出手段とを備えた。
Description
【０００１】
【産業上の利用分野】本発明は、超音波のドプラ効果を
利用して被検体内の血流情報を得る超音波ドプラ診断装
置に関し、詳細には低速の血流成分の情報をクラッタ成
分の情報から分離して検出するようにした超音波ドプラ
診断装置に関する。
利用して被検体内の血流情報を得る超音波ドプラ診断装
置に関し、詳細には低速の血流成分の情報をクラッタ成
分の情報から分離して検出するようにした超音波ドプラ
診断装置に関する。
【０００２】
【従来の技術】人体内に超音波ビームを送信し、人体内
の組織で反射されて戻ってくる超音波を受信して人体の
内蔵等の疾患の診断を行う超音波診断装置が従来より用
いられており、この超音波診断装置の一態様として、も
しくはＢモード（断層像）表示を行う超音波診断装置の
オプションとして、体内を流れる血球で反射された超音
波を受信して血流の速度、分散、パワー等の血流情報を
得る超音波ドプラ診断装置が用いられている。
の組織で反射されて戻ってくる超音波を受信して人体の
内蔵等の疾患の診断を行う超音波診断装置が従来より用
いられており、この超音波診断装置の一態様として、も
しくはＢモード（断層像）表示を行う超音波診断装置の
オプションとして、体内を流れる血球で反射された超音
波を受信して血流の速度、分散、パワー等の血流情報を
得る超音波ドプラ診断装置が用いられている。
【０００３】図１０は従来の超音波診断装置の一例の概
略構成図である。送信制御部１１から超音波プローブ１
２を構成する多数の振動子（図示せず）にそれぞれ所定
のタイミングでパルス信号Ｔｐが送信され、これにより
被検体（図示せず）内に超音波パルスビームが送信され
る。ここでは例えばセクタ走査がなされ、各方向に例え
ば８回のパルス状の超音波ビームが送信される。この超
音波パルスビームは、被検体内を流れる血球やその他の
組織で反射され、この超音波プローブ１２内の多数の振
動子で受信される。この各振動子で受信された受信信号
Ａｐは、遅延加算部１３に入力され、この遅延加算部１
３内で受信ダイナミックフォーカスを満足するように遅
延加算される。この遅延加算後の受信信号Ｓは、Ｂモー
ド検出部１４と血流情報検出部１５に入力される。
略構成図である。送信制御部１１から超音波プローブ１
２を構成する多数の振動子（図示せず）にそれぞれ所定
のタイミングでパルス信号Ｔｐが送信され、これにより
被検体（図示せず）内に超音波パルスビームが送信され
る。ここでは例えばセクタ走査がなされ、各方向に例え
ば８回のパルス状の超音波ビームが送信される。この超
音波パルスビームは、被検体内を流れる血球やその他の
組織で反射され、この超音波プローブ１２内の多数の振
動子で受信される。この各振動子で受信された受信信号
Ａｐは、遅延加算部１３に入力され、この遅延加算部１
３内で受信ダイナミックフォーカスを満足するように遅
延加算される。この遅延加算後の受信信号Ｓは、Ｂモー
ド検出部１４と血流情報検出部１５に入力される。
【０００４】Ｂモード検出部１４では、入力された受信
信号Ｓに基づいてＢモード表示（断層像）を担う信号Ｓ
A が生成され、この信号ＳA がＣＲＴディスプレイ装置
等からなる表示部１６に入力され、この表示部１６にお
いて断層像が表示され診断に供される。また血流情報検
出部１５では、入力された受信信号Ｓに基づいて、以下
に示すドプラ効果が利用されて血流情報が検出される。
信号Ｓに基づいてＢモード表示（断層像）を担う信号Ｓ
A が生成され、この信号ＳA がＣＲＴディスプレイ装置
等からなる表示部１６に入力され、この表示部１６にお
いて断層像が表示され診断に供される。また血流情報検
出部１５では、入力された受信信号Ｓに基づいて、以下
に示すドプラ効果が利用されて血流情報が検出される。
【０００５】即ち、血流内の血球によって反射された反
射超音波は、血球の移動により周波数偏移を受け、この
偏移量（ドプラ偏移周波数）ｆｄは、体内を流れる血流
の速度をＶ、この血流の方向と送信された超音波ビーム
の方向とのなす角をθ、送信超音波の中心周波数をｆ
ｃ、体内を伝播する超音波の速度をＣとしたとき、 ｆｄ ＝ （２Ｖｃｏｓθ／Ｃ）・ｆｃ …（１） で表わされ、反射超音波を受信した受信信号の中心周波
数ｆは ｆ ＝ ｆｃ＋ｆｄ …（２） となる。したがって、このドプラ偏移周波数ｆｄを検出
し、また上記断層像を担う信号ＳA に基づいて血管の延
びる方向を検出することにより、血流速度Ｖを検出する
ことができる。このドプラ偏移周波数ｆｄは、例えば広
く用いられている自己相関法、ＦＦＴ法や、そのほか微
小変位計測法（相互相関法（例えば、日本超音波医学会
第５４回研究発表会講演論文集第３５９頁〜第３６０頁
「解析信号の空間相関を用いた不均一組織の微小変位計
測」八木他参照）、位相追跡法（例えば、日本超音波医
学会第５７回研究発表会講演論文集第４４５頁〜第４４
６頁「位相追跡処理による生体内組織変位測定」新木
他、平成２年１０月１２日特許出願「超音波診断装置」
特願平２−２７３９１０号参照）、観測データにオリエ
ンテッドな方法（日本超音波医学会第５６回研究発表会
講演論文集第２３３頁〜第２３４頁「散乱体のランダム
な構造に依存しない反射型での微小変動の推定法」山越
他、平成２年４月３日特許出願「超音波診断装置」特願
平２−０８８５５３号参照））等の各種の方法により求
めることができる。
射超音波は、血球の移動により周波数偏移を受け、この
偏移量（ドプラ偏移周波数）ｆｄは、体内を流れる血流
の速度をＶ、この血流の方向と送信された超音波ビーム
の方向とのなす角をθ、送信超音波の中心周波数をｆ
ｃ、体内を伝播する超音波の速度をＣとしたとき、 ｆｄ ＝ （２Ｖｃｏｓθ／Ｃ）・ｆｃ …（１） で表わされ、反射超音波を受信した受信信号の中心周波
数ｆは ｆ ＝ ｆｃ＋ｆｄ …（２） となる。したがって、このドプラ偏移周波数ｆｄを検出
し、また上記断層像を担う信号ＳA に基づいて血管の延
びる方向を検出することにより、血流速度Ｖを検出する
ことができる。このドプラ偏移周波数ｆｄは、例えば広
く用いられている自己相関法、ＦＦＴ法や、そのほか微
小変位計測法（相互相関法（例えば、日本超音波医学会
第５４回研究発表会講演論文集第３５９頁〜第３６０頁
「解析信号の空間相関を用いた不均一組織の微小変位計
測」八木他参照）、位相追跡法（例えば、日本超音波医
学会第５７回研究発表会講演論文集第４４５頁〜第４４
６頁「位相追跡処理による生体内組織変位測定」新木
他、平成２年１０月１２日特許出願「超音波診断装置」
特願平２−２７３９１０号参照）、観測データにオリエ
ンテッドな方法（日本超音波医学会第５６回研究発表会
講演論文集第２３３頁〜第２３４頁「散乱体のランダム
な構造に依存しない反射型での微小変動の推定法」山越
他、平成２年４月３日特許出願「超音波診断装置」特願
平２−０８８５５３号参照））等の各種の方法により求
めることができる。
【０００６】上記のようにして求められた血流情報を表
わす信号ＳB は表示部１６に入力され、例えば上記断層
像に重畳されて、超音波プローブに近づく方向の血流は
例えば赤色、遠ざかる血流は例えば青色で表示される。
以下、ここでの主題である血流情報の検出に限って説明
する。図１１は、従来の超音波ドプラ診断装置の、図１
０に示す血流情報検出部１５に相当する部分の一例を表
わしたブロック図である。
わす信号ＳB は表示部１６に入力され、例えば上記断層
像に重畳されて、超音波プローブに近づく方向の血流は
例えば赤色、遠ざかる血流は例えば青色で表示される。
以下、ここでの主題である血流情報の検出に限って説明
する。図１１は、従来の超音波ドプラ診断装置の、図１
０に示す血流情報検出部１５に相当する部分の一例を表
わしたブロック図である。
【０００７】遅延加算された受信信号Ｓは直交検波回路
１７に入力されて直交検波される。図１２は、この直交
検波回路１７の内部構成を表わすブロック図である。こ
の直交検波回路１７に入力された受信信号Ｓは２系統に
分割され、それぞれが乗算器１７ａ，１７ｂに入力され
る。またこれらの乗算器１７ａ，１７ｂには、互いに９
０°位相の異なる２つの正弦波信号（キャリア信号）Ｓ
ＩＮ，ＣＯＳがそれぞれ入力され、これらの乗算器１７
ａ，１７ｂでは入力された各２つの信号（受信信号Ｓと
キャリア信号ＳＩＮ、又は受信信号Ｓとキャリア信号Ｃ
ＯＳ）がそれぞれ乗算され、これにより乗算前の各２つ
の信号の和と差の周波数の双方をもつ２つの信号ＳＭ
I ，ＳＭQ が生成される。これらの各信号ＳＭI ，ＳＭ
Qはそれぞれ各低域通過フィルタ１７ｃ，１７ｄを経由
することにより上記各２つの信号の差の周波数の信号の
みを担持する、直交検波後の受信信号ＳのＩ成分Ｓ I ，
Ｑ成分ＳQ が生成されこの直交検波回路１７から出力さ
れる。この直交検波された受信信号ＳのＩ成分ＳI ，Ｑ
成分ＳQ は、この直交検波回路１７から出力された後、
図１１に示す各Ａ／Ｄ変換器１８ａ，１８ｂに入力され
てＡ／Ｄ変換され、その後各ＭＴＩ（Ｍｏｖｉｎｇ Ｔ
ａｒｇｅｔ Ｉｎｄｉｃａｔｏｒ）フィルタ１９ａ，１
９ｂに入力される。このＭＴＩフィルタ１９ａ，１９ｂ
は、レーダーにおけるＭＴＩフィルタと同様のものであ
り、通常はパルス信号の繰返し周期に相当する遅延時間
をもつ遅延回路と積和器とで構成される、低周波信号を
遮断するディジタルフィルタであって、この超音波ドプ
ラ診断装置の分野において広く用いられているものであ
る。このＭＴＩフィルタ１９ａ，１９ｂは、受信信号Ｓ
には血流成分の情報とともに血流以外の被検体の動き等
に起因する比較的ゆるやかな動きの情報を担う、血流成
分と比較し例えば１００倍程度のパワーをもつクラッタ
成分の情報が強大な雑音成分として混入しているため、
これを取り除く目的で用いられるものである。このＭＴ
Ｉフィルタ１９ａ，１９ｂから出力されたクラッタ除去
信号ＳＣのＩ成分ＳＣI ，Ｑ成分ＳＣQ は相関器２０に
入力され、この相関器２０における自己相関演算により
血流の速度、速度分散、パワー等の血流情報を担持する
信号ＳB が生成される。
１７に入力されて直交検波される。図１２は、この直交
検波回路１７の内部構成を表わすブロック図である。こ
の直交検波回路１７に入力された受信信号Ｓは２系統に
分割され、それぞれが乗算器１７ａ，１７ｂに入力され
る。またこれらの乗算器１７ａ，１７ｂには、互いに９
０°位相の異なる２つの正弦波信号（キャリア信号）Ｓ
ＩＮ，ＣＯＳがそれぞれ入力され、これらの乗算器１７
ａ，１７ｂでは入力された各２つの信号（受信信号Ｓと
キャリア信号ＳＩＮ、又は受信信号Ｓとキャリア信号Ｃ
ＯＳ）がそれぞれ乗算され、これにより乗算前の各２つ
の信号の和と差の周波数の双方をもつ２つの信号ＳＭ
I ，ＳＭQ が生成される。これらの各信号ＳＭI ，ＳＭ
Qはそれぞれ各低域通過フィルタ１７ｃ，１７ｄを経由
することにより上記各２つの信号の差の周波数の信号の
みを担持する、直交検波後の受信信号ＳのＩ成分Ｓ I ，
Ｑ成分ＳQ が生成されこの直交検波回路１７から出力さ
れる。この直交検波された受信信号ＳのＩ成分ＳI ，Ｑ
成分ＳQ は、この直交検波回路１７から出力された後、
図１１に示す各Ａ／Ｄ変換器１８ａ，１８ｂに入力され
てＡ／Ｄ変換され、その後各ＭＴＩ（Ｍｏｖｉｎｇ Ｔ
ａｒｇｅｔ Ｉｎｄｉｃａｔｏｒ）フィルタ１９ａ，１
９ｂに入力される。このＭＴＩフィルタ１９ａ，１９ｂ
は、レーダーにおけるＭＴＩフィルタと同様のものであ
り、通常はパルス信号の繰返し周期に相当する遅延時間
をもつ遅延回路と積和器とで構成される、低周波信号を
遮断するディジタルフィルタであって、この超音波ドプ
ラ診断装置の分野において広く用いられているものであ
る。このＭＴＩフィルタ１９ａ，１９ｂは、受信信号Ｓ
には血流成分の情報とともに血流以外の被検体の動き等
に起因する比較的ゆるやかな動きの情報を担う、血流成
分と比較し例えば１００倍程度のパワーをもつクラッタ
成分の情報が強大な雑音成分として混入しているため、
これを取り除く目的で用いられるものである。このＭＴ
Ｉフィルタ１９ａ，１９ｂから出力されたクラッタ除去
信号ＳＣのＩ成分ＳＣI ，Ｑ成分ＳＣQ は相関器２０に
入力され、この相関器２０における自己相関演算により
血流の速度、速度分散、パワー等の血流情報を担持する
信号ＳB が生成される。
【０００８】
【発明が解決しようとする課題】図１３は、ＭＴＩフィ
ルタの特性の例を表わした図である。横軸はドプラ偏移
周波数ｆｄ（式（１）参照）を表わしている。また折れ
線３１はＭＴＩフィルタの特性を示しており、このＭＴ
Ｉフィルタはｆｄ＝０を中心とした｜ｆｄ｜≦Ｔｈの周
波数帯域の信号をカットし、｜ｆｄ｜＞Ｔｈの周波数帯
域の信号を通過させる特性を有している。また山３２，
３３は受信信号Ｓが担持するそれぞれクラッタ成分、血
流成分のドプラ偏移周波数分布を表わしている。
ルタの特性の例を表わした図である。横軸はドプラ偏移
周波数ｆｄ（式（１）参照）を表わしている。また折れ
線３１はＭＴＩフィルタの特性を示しており、このＭＴ
Ｉフィルタはｆｄ＝０を中心とした｜ｆｄ｜≦Ｔｈの周
波数帯域の信号をカットし、｜ｆｄ｜＞Ｔｈの周波数帯
域の信号を通過させる特性を有している。また山３２，
３３は受信信号Ｓが担持するそれぞれクラッタ成分、血
流成分のドプラ偏移周波数分布を表わしている。
【０００９】図１３（ａ）に示すように、血流の流れが
早く山３３が山３２と大きく離れている場合には、山３
２に対応するクラッタ成分の情報をカットするとともに
山３３に対応する血流成分の情報を通過させるようにＭ
ＴＩフィルタの信号遮断帯域を定めることにより、クラ
ッタ成分の情報が選択的に除去されたクラッタ除去信号
ＳＣ（ＳＣI ，ＳＣQ ）が出力されるが、血流の速度が
非常に遅い場合は、例えば図１３（ｂ）に示すように山
３３が山３２に近づき、これらの山３２，３３が互いに
分離している場合であっても、クラッタ成分の情報を除
去するようにＭＴＩフィルタの特性を定めると血流成分
の情報まで除去されてしまい、クラッタ成分の情報のみ
を選択的に除去することができないという問題がある。
特に近年では例えば肝臓等腹部の血流情報を検出する要
請が高まり、このためクラッタ成分と同程度のドプラ偏
移周波数をもつ非常に速度の遅い血流まで検出すること
が重要となってきている。
早く山３３が山３２と大きく離れている場合には、山３
２に対応するクラッタ成分の情報をカットするとともに
山３３に対応する血流成分の情報を通過させるようにＭ
ＴＩフィルタの信号遮断帯域を定めることにより、クラ
ッタ成分の情報が選択的に除去されたクラッタ除去信号
ＳＣ（ＳＣI ，ＳＣQ ）が出力されるが、血流の速度が
非常に遅い場合は、例えば図１３（ｂ）に示すように山
３３が山３２に近づき、これらの山３２，３３が互いに
分離している場合であっても、クラッタ成分の情報を除
去するようにＭＴＩフィルタの特性を定めると血流成分
の情報まで除去されてしまい、クラッタ成分の情報のみ
を選択的に除去することができないという問題がある。
特に近年では例えば肝臓等腹部の血流情報を検出する要
請が高まり、このためクラッタ成分と同程度のドプラ偏
移周波数をもつ非常に速度の遅い血流まで検出すること
が重要となってきている。
【００１０】本発明は、上記事情に鑑み、非常に遅い速
度をもつ血流であってもそのドプラ偏移周波数がクラッ
タ成分のドプラ変移周波数と異なっている限り、クラッ
タ成分の情報のみを選択的に除去することのできる超音
波ドプラ診断装置を提供することを目的とする。
度をもつ血流であってもそのドプラ偏移周波数がクラッ
タ成分のドプラ変移周波数と異なっている限り、クラッ
タ成分の情報のみを選択的に除去することのできる超音
波ドプラ診断装置を提供することを目的とする。
【００１１】
【課題を解決するための手段】図１は、本発明の超音波
ドプラ診断装置の特徴部分の原理ブロック図である。以
後、簡単のため、特徴部分のみを指標する場合も超音波
ドプラ診断装置と称する。前述のようにして遅延加算さ
れた受信信号Ｓは、クラッタ情報検出手段１とクラッタ
情報除去手段２との双方に入力される。
ドプラ診断装置の特徴部分の原理ブロック図である。以
後、簡単のため、特徴部分のみを指標する場合も超音波
ドプラ診断装置と称する。前述のようにして遅延加算さ
れた受信信号Ｓは、クラッタ情報検出手段１とクラッタ
情報除去手段２との双方に入力される。
【００１２】このクラッタ情報検出手段１では入力され
た受信信号Ｓに基づいてクラッタ成分の移動速度Ｖｃと
その移動速度の分散σc 2との双方又はこれら移動速度Ｖ
ｃと分散σc 2のいずれか一方が求められる。ここでクラ
ッタの移動速度Ｖｃ、分散σc 2はクラッタ成分のそれぞ
れ平均的な移動速度、移動速度の変動の大きさを総称す
るものであり、数学的な意味における平均速度、分散で
ある必要はない。例えば移動速度Ｖｃとして多数回測定
した際の中央値、（最大値＋最小値）／２等を採用して
もよく、分散σc 2として、標準偏差、（最大値−最少
値）等を採用してもよい。
た受信信号Ｓに基づいてクラッタ成分の移動速度Ｖｃと
その移動速度の分散σc 2との双方又はこれら移動速度Ｖ
ｃと分散σc 2のいずれか一方が求められる。ここでクラ
ッタの移動速度Ｖｃ、分散σc 2はクラッタ成分のそれぞ
れ平均的な移動速度、移動速度の変動の大きさを総称す
るものであり、数学的な意味における平均速度、分散で
ある必要はない。例えば移動速度Ｖｃとして多数回測定
した際の中央値、（最大値＋最小値）／２等を採用して
もよく、分散σc 2として、標準偏差、（最大値−最少
値）等を採用してもよい。
【００１３】この移動速度Ｖｃ及び／又は分散σc 2の求
め方は、特定の求め方に限定されるものではなく、例え
ば前述した相互相関法、位相追跡法、観測データにオリ
エンテッドな方法等の各種の微小変位計測法や、自己相
関法、ＦＦＴ法のいずれかあるいはこれらの組合わせ等
を用いることができ、例えば被検体内の一点の血流情報
を求めるいわゆるポイントドプラ法や２次元的な断層面
内の血流情報を求めるいわゆるカラードプラ法等に応
じ、それぞれ適切な求め方が選択される。
め方は、特定の求め方に限定されるものではなく、例え
ば前述した相互相関法、位相追跡法、観測データにオリ
エンテッドな方法等の各種の微小変位計測法や、自己相
関法、ＦＦＴ法のいずれかあるいはこれらの組合わせ等
を用いることができ、例えば被検体内の一点の血流情報
を求めるいわゆるポイントドプラ法や２次元的な断層面
内の血流情報を求めるいわゆるカラードプラ法等に応
じ、それぞれ適切な求め方が選択される。
【００１４】クラッタ情報検出手段１で求められた移動
速度Ｖｃ及び／又は分散σc 2は、クラッタ情報除去手段
２に入力される。このクラッタ情報除去手段２は、上記
クラッタ情報検出手段１で求められた移動速度Ｖｃ及び
／又は分散σc 2に基づいて、このクラッタ情報除去手段
２に入力された受信信号Ｓからこの受信信号Ｓに含まれ
るクラッタ成分の情報を選択的に除去するものであり、
これにより血流成分の情報のみを含むクラッタ除去信号
ＳＣが生成される。
速度Ｖｃ及び／又は分散σc 2は、クラッタ情報除去手段
２に入力される。このクラッタ情報除去手段２は、上記
クラッタ情報検出手段１で求められた移動速度Ｖｃ及び
／又は分散σc 2に基づいて、このクラッタ情報除去手段
２に入力された受信信号Ｓからこの受信信号Ｓに含まれ
るクラッタ成分の情報を選択的に除去するものであり、
これにより血流成分の情報のみを含むクラッタ除去信号
ＳＣが生成される。
【００１５】このクラッタ情報除去手段２も特定の構成
に限定されるものではない。例えば後述する実施例に示
すように、クラッタ情報検出手段１でクラッタ成分の移
動速度Ｖｃのみ、もしくは移動速度Ｖｃと分散σc 2との
双方を求め、求められた移動速度Ｖｃ、もしくは移動速
度Ｖｃと分散σc 2とに基づいて複素ＭＴＩフィルタの特
性を定めるようにしてもよく、又はクラッタ情報検出手
段１でクラッタ成分の移動速度Ｖｃのみ、もしくは移動
速度Ｖｃと分散σc 2との双方を求め、この移動速度Ｖｃ
に基づいて直交検波回路から出力される受信信号の直流
成分がクラッタ成分の移動速度の情報を担持するように
直交検波のキャリア信号の位相を制御し、この直交検波
回路から出力された受信信号を、信号遮断帯域の固定さ
れたＭＴＩフィルタもしくは分散σc 2に基づいてその帯
域幅が調整されたＭＴＩフィルタに入力するように構成
してもよい。また、クラッタ成分の移動速度ＶｃがＶｃ
≒０であることはあらかじめわかっているが、その分散
σc 2が大きくなったり小さくなったり変動する場合に
は、この分散σc 2が小さくなったときに血流成分の情報
を取出すために、クラッタ情報検出手段１でクラッタ成
分の移動速度の分散σ c 2を求め、この分散σc 2に基づい
てＭＴＩフィルタの信号遮断帯域の帯域幅を定めるよう
に構成してもよい。またＭＴＩフィルタや複素ＭＴＩフ
ィルタを用いる各種構成のほか、例えばポイントドップ
ラ法の場合等に、クラッタ成分と血流成分との双方の情
報を含む受信信号Ｓに基づいてクラッタ成分の移動速度
Ｖｃとその分散σc 2とを求め（クラッタ情報検出手段１
の機能）、その後、この移動速度Ｖｃと分散σc 2とに基
づいて周波数空間上でクラッタ成分の情報を除去するよ
うに構成してもよい。
に限定されるものではない。例えば後述する実施例に示
すように、クラッタ情報検出手段１でクラッタ成分の移
動速度Ｖｃのみ、もしくは移動速度Ｖｃと分散σc 2との
双方を求め、求められた移動速度Ｖｃ、もしくは移動速
度Ｖｃと分散σc 2とに基づいて複素ＭＴＩフィルタの特
性を定めるようにしてもよく、又はクラッタ情報検出手
段１でクラッタ成分の移動速度Ｖｃのみ、もしくは移動
速度Ｖｃと分散σc 2との双方を求め、この移動速度Ｖｃ
に基づいて直交検波回路から出力される受信信号の直流
成分がクラッタ成分の移動速度の情報を担持するように
直交検波のキャリア信号の位相を制御し、この直交検波
回路から出力された受信信号を、信号遮断帯域の固定さ
れたＭＴＩフィルタもしくは分散σc 2に基づいてその帯
域幅が調整されたＭＴＩフィルタに入力するように構成
してもよい。また、クラッタ成分の移動速度ＶｃがＶｃ
≒０であることはあらかじめわかっているが、その分散
σc 2が大きくなったり小さくなったり変動する場合に
は、この分散σc 2が小さくなったときに血流成分の情報
を取出すために、クラッタ情報検出手段１でクラッタ成
分の移動速度の分散σ c 2を求め、この分散σc 2に基づい
てＭＴＩフィルタの信号遮断帯域の帯域幅を定めるよう
に構成してもよい。またＭＴＩフィルタや複素ＭＴＩフ
ィルタを用いる各種構成のほか、例えばポイントドップ
ラ法の場合等に、クラッタ成分と血流成分との双方の情
報を含む受信信号Ｓに基づいてクラッタ成分の移動速度
Ｖｃとその分散σc 2とを求め（クラッタ情報検出手段１
の機能）、その後、この移動速度Ｖｃと分散σc 2とに基
づいて周波数空間上でクラッタ成分の情報を除去するよ
うに構成してもよい。
【００１６】このようにクラッタ情報除去手段２として
種々の構成を採用することができる。このクラッタ情報
除去手段２から出力されたクラック除去信号ＳＣは、血
流情報検出手段３に入力される。この血流情報検出手段
３では入力されたクラッタ除去信号ＳＣに基づいて血流
情報を表わす信号ＳB が求められる。この血流情報検出
手段３における血流情報の求め方も特定の手法に限定さ
れるものではなく、クラッタ除去信号ＳＣに基づいて、
前述した各種の方法、即ち相互相関法、位相追跡法、観
測データにオリエンテッドな方法等の各種の微小変位計
測法や、自己相関法、ＦＦＴ法のいずれかあるいはこれ
らの組合わせ等を用いて血流情報を求めることができ
る。
種々の構成を採用することができる。このクラッタ情報
除去手段２から出力されたクラック除去信号ＳＣは、血
流情報検出手段３に入力される。この血流情報検出手段
３では入力されたクラッタ除去信号ＳＣに基づいて血流
情報を表わす信号ＳB が求められる。この血流情報検出
手段３における血流情報の求め方も特定の手法に限定さ
れるものではなく、クラッタ除去信号ＳＣに基づいて、
前述した各種の方法、即ち相互相関法、位相追跡法、観
測データにオリエンテッドな方法等の各種の微小変位計
測法や、自己相関法、ＦＦＴ法のいずれかあるいはこれ
らの組合わせ等を用いて血流情報を求めることができ
る。
【００１７】
【作用】クラッタ情報除去手段１に入力される受信信号
Ｓにはクラッタ成分の情報のみならず血流成分の情報も
含まれているが、前述したようにクラッタ成分は血流成
分よりも極端に大きなパワーを有している。従って受信
信号Ｓには血流成分の情報も含まれてはいるものの、こ
の受信信号Ｓに基づいてこの受信信号Ｓが担持する全成
分の移動速度や分散等を求めても、その結果は実質上ク
ラッタ成分の移動速度Ｖｃ、その分散σc 2と見なして差
しつかえない。
Ｓにはクラッタ成分の情報のみならず血流成分の情報も
含まれているが、前述したようにクラッタ成分は血流成
分よりも極端に大きなパワーを有している。従って受信
信号Ｓには血流成分の情報も含まれてはいるものの、こ
の受信信号Ｓに基づいてこの受信信号Ｓが担持する全成
分の移動速度や分散等を求めても、その結果は実質上ク
ラッタ成分の移動速度Ｖｃ、その分散σc 2と見なして差
しつかえない。
【００１８】本発明はこの点に着目して完成されたもの
であり、受信信号Ｓに基づいて先ずクラッタ成分の移動
速度Ｖｃ及び／又はその分散σc 2を求め、この移動速度
Ｖｃ及び／又は分散σc 2を用いて受信信号Ｓからクラッ
タ成分の情報を選択的に除去するように構成したため、
クラッタ成分と血流成分のドプラ偏移周波数ｆｄが互い
に異なっている限りクラッタ成分のみを選択的に除去す
ることができ、したがって例えば腹部等速度の遅い血流
の情報を得る場合に高性能の超音波ドプラ診断装置が実
現される。
であり、受信信号Ｓに基づいて先ずクラッタ成分の移動
速度Ｖｃ及び／又はその分散σc 2を求め、この移動速度
Ｖｃ及び／又は分散σc 2を用いて受信信号Ｓからクラッ
タ成分の情報を選択的に除去するように構成したため、
クラッタ成分と血流成分のドプラ偏移周波数ｆｄが互い
に異なっている限りクラッタ成分のみを選択的に除去す
ることができ、したがって例えば腹部等速度の遅い血流
の情報を得る場合に高性能の超音波ドプラ診断装置が実
現される。
【００１９】
【実施例】以下、本発明の実施例について説明する。図
２は、本発明の超音波ドプラ診断装置の第一の実施例を
表わすブロック図、図３は図２に示す第一の実施例にお
ける複素ＭＴＩフィルタの特性と受信信号との関係を表
わした図である。これらの図において図１１〜図１３に
示した従来例の構成要素と対応する要素には図１１〜図
１３に付した番号、記号と同一の番号、記号を付し、説
明は省略する。
２は、本発明の超音波ドプラ診断装置の第一の実施例を
表わすブロック図、図３は図２に示す第一の実施例にお
ける複素ＭＴＩフィルタの特性と受信信号との関係を表
わした図である。これらの図において図１１〜図１３に
示した従来例の構成要素と対応する要素には図１１〜図
１３に付した番号、記号と同一の番号、記号を付し、説
明は省略する。
【００２０】複素ＭＴＩフィルタ２２は、通常用いられ
ているＭＴＩフィルタ（図１１、図１３参照）と比べフ
ィルタを構成する素子数が多く回路規模が非常に大きい
という欠点を有するが、一方、信号遮断帯域の中心周波
数ｆｏをｆｏ＝０以外の周波数とすることができるとい
う特性を備えている。この実施例はその特性を生かした
ものである。
ているＭＴＩフィルタ（図１１、図１３参照）と比べフ
ィルタを構成する素子数が多く回路規模が非常に大きい
という欠点を有するが、一方、信号遮断帯域の中心周波
数ｆｏをｆｏ＝０以外の周波数とすることができるとい
う特性を備えている。この実施例はその特性を生かした
ものである。
【００２１】２つのＡ／Ｄ変換器１８ａ，１８ｂでディ
ジタル化された受信信号ＳのＩ成分ＳI ，Ｑ成分ＳQ は
相関器１９に入力され、この相関器１９で自己相関演算
がなされ、これによりクラッタ成分の移動速度Ｖｃとそ
の分散σc 2が求められる。ここで、この相関器１９に入
力された受信信号Ｓにはクラッタ成分と血流成分の双方
の情報が含まれるが、クラッタ成分は血流成分と比べ１
００倍（４０ｄＢ）程度のパワーを有しているため、血
流成分の情報を含んだまま自己相関演算を行ってもほと
んど問題なくクラッタ成分の移動速度Ｖｃ、分散σc 2が
求められる。この求められたクラッタ成分の移動速度Ｖ
ｃ、分散σc 2を表わす信号はメモリ２０に入力される。
このメモリ２０は、本発明にいうフィルタ制御回路に相
当するものであり、メモリ２０内には移動速度Ｖｃと複
素ＭＴＩフィルタ２２の信号遮断帯域の中心周波数ｆｏ
を定める係数との対応表、および分散σc 2と複素ＭＴＩ
フィルタ２２の信号遮断帯域の帯域幅Ｗを定める係数と
の対応表がルックアップテーブルの形式であらかじめ記
憶されており、相関器１９で求められた移動速度Ｖｃ及
び分散σc 2がメモリ２０に入力されるとこのメモリ２０
から複素ＭＴＩフィルタ２２の信号遮断帯域の中心周波
数ｆｏを定める係数と帯域幅Ｗを定める係数が出力さ
れ、これらの係数が複素ＭＴＩフィルタ２２に入力され
る。
ジタル化された受信信号ＳのＩ成分ＳI ，Ｑ成分ＳQ は
相関器１９に入力され、この相関器１９で自己相関演算
がなされ、これによりクラッタ成分の移動速度Ｖｃとそ
の分散σc 2が求められる。ここで、この相関器１９に入
力された受信信号Ｓにはクラッタ成分と血流成分の双方
の情報が含まれるが、クラッタ成分は血流成分と比べ１
００倍（４０ｄＢ）程度のパワーを有しているため、血
流成分の情報を含んだまま自己相関演算を行ってもほと
んど問題なくクラッタ成分の移動速度Ｖｃ、分散σc 2が
求められる。この求められたクラッタ成分の移動速度Ｖ
ｃ、分散σc 2を表わす信号はメモリ２０に入力される。
このメモリ２０は、本発明にいうフィルタ制御回路に相
当するものであり、メモリ２０内には移動速度Ｖｃと複
素ＭＴＩフィルタ２２の信号遮断帯域の中心周波数ｆｏ
を定める係数との対応表、および分散σc 2と複素ＭＴＩ
フィルタ２２の信号遮断帯域の帯域幅Ｗを定める係数と
の対応表がルックアップテーブルの形式であらかじめ記
憶されており、相関器１９で求められた移動速度Ｖｃ及
び分散σc 2がメモリ２０に入力されるとこのメモリ２０
から複素ＭＴＩフィルタ２２の信号遮断帯域の中心周波
数ｆｏを定める係数と帯域幅Ｗを定める係数が出力さ
れ、これらの係数が複素ＭＴＩフィルタ２２に入力され
る。
【００２２】また、Ａ／Ｄ変換器１８ａ，１８ｂでディ
ジタル化された受信信号ＳのＩ成分ＳI 、Ｑ成分ＳQ
は、各遅延器２１ａ，２１ｂにも入力され、相関器１９
で移動速度Ｖｃ、分散σc 2が求められメモリ２０で上記
係数が求められてその係数が複素ＭＴＩフィルタ２２に
入力されるまでの時間遅延された後、複素ＭＴＩフィル
タ２２に入力される。この複素ＭＴＩフィルタ２２は、
図３に示すように、クラッタ成分（山３２）のみを選択
的に除去するように上記係数によりその特性が定められ
ており、したがって血流成分（山３３）が有効に取り出
される。
ジタル化された受信信号ＳのＩ成分ＳI 、Ｑ成分ＳQ
は、各遅延器２１ａ，２１ｂにも入力され、相関器１９
で移動速度Ｖｃ、分散σc 2が求められメモリ２０で上記
係数が求められてその係数が複素ＭＴＩフィルタ２２に
入力されるまでの時間遅延された後、複素ＭＴＩフィル
タ２２に入力される。この複素ＭＴＩフィルタ２２は、
図３に示すように、クラッタ成分（山３２）のみを選択
的に除去するように上記係数によりその特性が定められ
ており、したがって血流成分（山３３）が有効に取り出
される。
【００２３】尚、クラッタ成分の移動速度の分散σc 2に
ついては、比較的実験的、経験的に求めやすいため、複
素ＭＴＩフィルタ２２の信号遮断帯域の帯域幅Ｗはあら
かじめ実験的、経験的に定められた値に固定し、相関器
１９ではクラッタ成分の移動速度Ｖｃのみを求め、これ
に基づいて複素ＭＴＩフィルタ２２の信号遮断帯域の中
心周波数ｆｏを定めてもよい。ただしこの固定された帯
域幅Ｗは被検体の測定部位に応じて変更することが好ま
しい。
ついては、比較的実験的、経験的に求めやすいため、複
素ＭＴＩフィルタ２２の信号遮断帯域の帯域幅Ｗはあら
かじめ実験的、経験的に定められた値に固定し、相関器
１９ではクラッタ成分の移動速度Ｖｃのみを求め、これ
に基づいて複素ＭＴＩフィルタ２２の信号遮断帯域の中
心周波数ｆｏを定めてもよい。ただしこの固定された帯
域幅Ｗは被検体の測定部位に応じて変更することが好ま
しい。
【００２４】また複素ＭＴＩフィルタ２２の信号遮断帯
域の中心周波数ｆｏ、帯域幅Ｗの定め方によっては、直
流成分（ｆｄ＝０）がこの信号遮断帯域から外れてしま
うこともあり得るが、例えば直交検波回路１７を構成す
る乗算器１７ａ，１７ｂ（図１２参照）のオフセットや
図示しない増幅器のオフセット等に起因する不要な直流
成分をもつ信号が複数ＭＴＩフィルタ２２に入力される
ことも有り得るため、この複素ＭＴＩフィルタ２２で入
力される信号の直流成分を遮断するようにその信号遮断
帯域にｆｄ＝０を含ませることが好ましい。
域の中心周波数ｆｏ、帯域幅Ｗの定め方によっては、直
流成分（ｆｄ＝０）がこの信号遮断帯域から外れてしま
うこともあり得るが、例えば直交検波回路１７を構成す
る乗算器１７ａ，１７ｂ（図１２参照）のオフセットや
図示しない増幅器のオフセット等に起因する不要な直流
成分をもつ信号が複数ＭＴＩフィルタ２２に入力される
ことも有り得るため、この複素ＭＴＩフィルタ２２で入
力される信号の直流成分を遮断するようにその信号遮断
帯域にｆｄ＝０を含ませることが好ましい。
【００２５】図４は本発明の超音波ドプラ診断装置の第
二の実施例を表わすブロック図、図５は、図４に示す第
二の実施例におけるＭＴＩフィルタの特性と受信信号と
の関係を表わした図である。これらの図において、図１
１〜図１３および図２、図３に示した従来例、実施例の
構成要素と対応する要素には、それらの図に付した番
号、記号と同一の番号、記号を付し、説明は省略する。
二の実施例を表わすブロック図、図５は、図４に示す第
二の実施例におけるＭＴＩフィルタの特性と受信信号と
の関係を表わした図である。これらの図において、図１
１〜図１３および図２、図３に示した従来例、実施例の
構成要素と対応する要素には、それらの図に付した番
号、記号と同一の番号、記号を付し、説明は省略する。
【００２６】この実施例では直交検波回路１７に入力さ
れるキャリア信号ＳＩＮ、ＣＯＳ（図１２参照）の位相
制御を行う位相制御回路２４を備えることにより、複素
ＭＴＩフィルタは用いずに通常のＭＴＩフィルタ２５
ａ，２５ｂを用いるように構成したものである。２次元
的な断層像に重畳して、血流を青、赤等で表示するいわ
ゆるカラードプラ法の場合は、前述したように被検体内
の各方向に超音波パルスが例えば８回送信されるが、こ
こではそのうちの前半４回と後半４回とに分け、前半４
回の測定の間は相関器１９でクラッタ成分の移動速度Ｖ
ｃとその分散σc 2が求められ、後半４回の測定の間はそ
の受信信号ＳがＭＴＩフィルタ２５ａ，２５ｂを通過し
て相関器２３に入力され、この相関器２３により血流情
報が求められる。
れるキャリア信号ＳＩＮ、ＣＯＳ（図１２参照）の位相
制御を行う位相制御回路２４を備えることにより、複素
ＭＴＩフィルタは用いずに通常のＭＴＩフィルタ２５
ａ，２５ｂを用いるように構成したものである。２次元
的な断層像に重畳して、血流を青、赤等で表示するいわ
ゆるカラードプラ法の場合は、前述したように被検体内
の各方向に超音波パルスが例えば８回送信されるが、こ
こではそのうちの前半４回と後半４回とに分け、前半４
回の測定の間は相関器１９でクラッタ成分の移動速度Ｖ
ｃとその分散σc 2が求められ、後半４回の測定の間はそ
の受信信号ＳがＭＴＩフィルタ２５ａ，２５ｂを通過し
て相関器２３に入力され、この相関器２３により血流情
報が求められる。
【００２７】図６は、前半４回、即ち位相調整前のキャ
リア信号ＳＩＮと受信信号Ｓとの位相関係を表わした模
式図である。図に示す（１），（２），（３），（４）
はこの順で行われた４回の測定を表わしている。例えば
被検体内のある反射点で反射されて戻った超音波を受信
して受信信号を得た場合にこの反射点が所定の方向に移
動していたと仮定し、この受信信号Ｓを超音波パルスの
各送信のタイミングを起点として表示すると、（１）〜
（４）の各受信信号Ｓの波形に示すように、反射点の移
動に伴ってその位相が順次変化する。一方、キャリア信
号ＳＩＮはその位相が固定されており、したがってキャ
リア信号ＳＩＮと受信信号Ｓとの位相関係は各測定毎に
異なることとなり、この場合、上記反射点によるドプラ
偏移周波数ｆｄはｆｄ≠０として検出されることとな
る。
リア信号ＳＩＮと受信信号Ｓとの位相関係を表わした模
式図である。図に示す（１），（２），（３），（４）
はこの順で行われた４回の測定を表わしている。例えば
被検体内のある反射点で反射されて戻った超音波を受信
して受信信号を得た場合にこの反射点が所定の方向に移
動していたと仮定し、この受信信号Ｓを超音波パルスの
各送信のタイミングを起点として表示すると、（１）〜
（４）の各受信信号Ｓの波形に示すように、反射点の移
動に伴ってその位相が順次変化する。一方、キャリア信
号ＳＩＮはその位相が固定されており、したがってキャ
リア信号ＳＩＮと受信信号Ｓとの位相関係は各測定毎に
異なることとなり、この場合、上記反射点によるドプラ
偏移周波数ｆｄはｆｄ≠０として検出されることとな
る。
【００２８】図７は、後半４回、即ち位相調整後のキャ
リア信号ＳＩＮと受信信号Ｓとの位相関係を表わした模
式図である。この図に示す（５），（６），（７），
（８）は、図６に示す（１）〜（４）の４回の測定の後
の、この順で行われた４回の測定を表わしている。図６
の場合と同様に、受信信号Ｓは反射点の移動に伴ってそ
の位相が順次変化しているが、キャリア信号ＳＩＮと受
信信号Ｓとの相対的な位相関係は常に一定となるように
キャリア信号ＳＩＮの位相も変化される。この場合上記
反射点によるドプラ偏移周波数ｆｄはｆｄ＝０として検
出されることとなる。
リア信号ＳＩＮと受信信号Ｓとの位相関係を表わした模
式図である。この図に示す（５），（６），（７），
（８）は、図６に示す（１）〜（４）の４回の測定の後
の、この順で行われた４回の測定を表わしている。図６
の場合と同様に、受信信号Ｓは反射点の移動に伴ってそ
の位相が順次変化しているが、キャリア信号ＳＩＮと受
信信号Ｓとの相対的な位相関係は常に一定となるように
キャリア信号ＳＩＮの位相も変化される。この場合上記
反射点によるドプラ偏移周波数ｆｄはｆｄ＝０として検
出されることとなる。
【００２９】図４に示す本発明の第二の実施例は上記原
理を応用したものであり、先ず前半の４回の測定に基づ
いて相関器１９でクラッタ成分の移動速度Ｖｃとその分
散σ c 2とが求められ、移動速度Ｖｃは位相制御回路２４
に入力され、分散σc 2はメモリ２０に入力される。位相
制御回路２４は、前半４回の測定の間は、図６に示すよ
うに位相の固定されたキャリア信号ＳＩＮ，ＣＯＳを生
成して直交検波回路１７に入力するが、後半の４回につ
いては、入力された移動速度Ｖｃに基づいて、この移動
速度Ｖｃがドプラ偏移周波数ｆｄ＝０にあらわれるよう
に位相の調整されたキャリア信号ＳＩＮ，ＣＯＳが生成
されて直交検波回路１７に入力される。またメモリ２０
には、分散σc 2とＭＴＩフィルタ２５ａ，２６ｂの信号
遮断帯域の帯域幅Ｗを定める係数との対応表がルックア
ップテーブルの形式であらかじめ記憶されており、分散
σc 2が入力されると、ＭＴＩフィルタ２５ａ，２５ｂの
信号遮断帯域の帯域幅Ｗを定める係数が求められ、この
係数がＭＴＩフィルタ２５ａ，２５ｂに入力される。こ
のＭＴＩフィルタ２５ａ，２５ｂは複素ＭＴＩフィルタ
とは異なり、その信号遮断帯域の中心周波数ｆｏはｆｏ
＝０に固定されこれを変更することはできないが、後半
４回の測定においては、直交検波回路１７から出力され
る受信信号Ｓ（ＳI ，ＳQ ）は、クラッタ成分の情報を
ドプラ偏移周波数ｆｄ＝０を中心として分布する情報と
して担持する信号であるため、クラッタ成分の本来の移
動速度ＶｃがＶｃ≠０であっても通常のＭＴＩフィルタ
２５ａ，２５ｂによりクラッタ成分を選択的に除去する
ことが可能となる。この後半４回の測定において、ＭＴ
Ｉフィルタ２５ａ，２５ｂから出力されたクラッタ除去
信号ＳＣ（ＳＣI ，ＳＣQ ）が相関器２３に入力され、
この相関器２３における自己相関演算により血流情報を
表わす信号ＳB が生成される。
理を応用したものであり、先ず前半の４回の測定に基づ
いて相関器１９でクラッタ成分の移動速度Ｖｃとその分
散σ c 2とが求められ、移動速度Ｖｃは位相制御回路２４
に入力され、分散σc 2はメモリ２０に入力される。位相
制御回路２４は、前半４回の測定の間は、図６に示すよ
うに位相の固定されたキャリア信号ＳＩＮ，ＣＯＳを生
成して直交検波回路１７に入力するが、後半の４回につ
いては、入力された移動速度Ｖｃに基づいて、この移動
速度Ｖｃがドプラ偏移周波数ｆｄ＝０にあらわれるよう
に位相の調整されたキャリア信号ＳＩＮ，ＣＯＳが生成
されて直交検波回路１７に入力される。またメモリ２０
には、分散σc 2とＭＴＩフィルタ２５ａ，２６ｂの信号
遮断帯域の帯域幅Ｗを定める係数との対応表がルックア
ップテーブルの形式であらかじめ記憶されており、分散
σc 2が入力されると、ＭＴＩフィルタ２５ａ，２５ｂの
信号遮断帯域の帯域幅Ｗを定める係数が求められ、この
係数がＭＴＩフィルタ２５ａ，２５ｂに入力される。こ
のＭＴＩフィルタ２５ａ，２５ｂは複素ＭＴＩフィルタ
とは異なり、その信号遮断帯域の中心周波数ｆｏはｆｏ
＝０に固定されこれを変更することはできないが、後半
４回の測定においては、直交検波回路１７から出力され
る受信信号Ｓ（ＳI ，ＳQ ）は、クラッタ成分の情報を
ドプラ偏移周波数ｆｄ＝０を中心として分布する情報と
して担持する信号であるため、クラッタ成分の本来の移
動速度ＶｃがＶｃ≠０であっても通常のＭＴＩフィルタ
２５ａ，２５ｂによりクラッタ成分を選択的に除去する
ことが可能となる。この後半４回の測定において、ＭＴ
Ｉフィルタ２５ａ，２５ｂから出力されたクラッタ除去
信号ＳＣ（ＳＣI ，ＳＣQ ）が相関器２３に入力され、
この相関器２３における自己相関演算により血流情報を
表わす信号ＳB が生成される。
【００３０】この第二の実施例は、前述した第一の実施
例と比較し、回路規模の大きな複素ＭＴＩフィルタを用
いる必要がないという長所はあるが、その一方で前半の
測定におけるクラッタ成分の移動速度Ｖｃに基づいて後
半の測定におけるキャリア信号ＳＩＮ，ＣＯＳの位相を
調整するため、前半の測定におけるクラッタ成分の移動
速度と後半の測定の際のクラッタ成分の移動速度が異な
ることも考えられ、この場合クラッタ成分を選択的に除
去する精度が低下するおそれがあること、血流情報は後
半の測定の際に求めるため、血流情報を前述した第一の
実施例と同一の精度で求めるには超音波の送信回数を増
やす必要があること等の短所を併せもっている。
例と比較し、回路規模の大きな複素ＭＴＩフィルタを用
いる必要がないという長所はあるが、その一方で前半の
測定におけるクラッタ成分の移動速度Ｖｃに基づいて後
半の測定におけるキャリア信号ＳＩＮ，ＣＯＳの位相を
調整するため、前半の測定におけるクラッタ成分の移動
速度と後半の測定の際のクラッタ成分の移動速度が異な
ることも考えられ、この場合クラッタ成分を選択的に除
去する精度が低下するおそれがあること、血流情報は後
半の測定の際に求めるため、血流情報を前述した第一の
実施例と同一の精度で求めるには超音波の送信回数を増
やす必要があること等の短所を併せもっている。
【００３１】尚、この第二の実施例においても、前述し
た第一の実施例と同様に、相関器１９でクラッタ成分の
移動速度Ｖｃのみを求め、ＭＴＩフィルタ２５ａ，２５
ｂの信号遮断帯域の帯域幅Ｗはあらかじめ実験的、経験
的に求められた値に固定してもよく、この場合、この第
二の実施例では、分散σc 2と帯域幅Ｗを定める係数との
対応表を記憶しておくメモリ２０は不要となる。
た第一の実施例と同様に、相関器１９でクラッタ成分の
移動速度Ｖｃのみを求め、ＭＴＩフィルタ２５ａ，２５
ｂの信号遮断帯域の帯域幅Ｗはあらかじめ実験的、経験
的に求められた値に固定してもよく、この場合、この第
二の実施例では、分散σc 2と帯域幅Ｗを定める係数との
対応表を記憶しておくメモリ２０は不要となる。
【００３２】図８は、本発明の超音波ドプラ診断装置の
第三の実施例を表わすブロック図、図９は図８に示す第
三の実施例におけるＭＴＩフィルタの特性と受信信号と
の関係を表わした図である。これらの図において、前述
した図１１〜図１３に示した従来例、図２〜図７に示し
た各実施例の構成要素と対応する要素にはそれらの各図
に付した番号、記号と同一の番号、記号を付し、説明は
省略する。
第三の実施例を表わすブロック図、図９は図８に示す第
三の実施例におけるＭＴＩフィルタの特性と受信信号と
の関係を表わした図である。これらの図において、前述
した図１１〜図１３に示した従来例、図２〜図７に示し
た各実施例の構成要素と対応する要素にはそれらの各図
に付した番号、記号と同一の番号、記号を付し、説明は
省略する。
【００３３】この第三の実施例においては、図８に示す
ように、信号遮断帯域の中心周波数ｆｏがｆｏ＝０に固
定された通常のＭＴＩフィルタ２５ａ，２５ｂが備えら
れており、相関器１９ではクラッタの移動速度の分散σ
c 2が求められてこれに基づいてＭＴＩフィルタ２５ａ，
２５ｂの信号遮断帯域の帯域幅Ｗが制御される。この第
三の実施例は、例えば被検体内のある臓器が震えたり停
止したりしている場合等、クラッタ成分（臓器）は全体
としては移動しない（移動速度Ｖｃ≒０）が、図９
（ａ），（ｂ）に示すようにその分散σc 2（山３２の
幅）が大きくなったり小さくなったりしている場合等に
有用であり、ＭＴＩフィルタ２５ａ，２５ｂの信号遮断
帯域の帯域幅Ｗを臓器の震え（分散σc 2の大小）に追随
して変化させることにより、図９（ａ）に示すようにそ
の震え（分散σc 2）が小さくなった瞬間に血流情報（山
３３に対応する情報）を取出すことができる。
ように、信号遮断帯域の中心周波数ｆｏがｆｏ＝０に固
定された通常のＭＴＩフィルタ２５ａ，２５ｂが備えら
れており、相関器１９ではクラッタの移動速度の分散σ
c 2が求められてこれに基づいてＭＴＩフィルタ２５ａ，
２５ｂの信号遮断帯域の帯域幅Ｗが制御される。この第
三の実施例は、例えば被検体内のある臓器が震えたり停
止したりしている場合等、クラッタ成分（臓器）は全体
としては移動しない（移動速度Ｖｃ≒０）が、図９
（ａ），（ｂ）に示すようにその分散σc 2（山３２の
幅）が大きくなったり小さくなったりしている場合等に
有用であり、ＭＴＩフィルタ２５ａ，２５ｂの信号遮断
帯域の帯域幅Ｗを臓器の震え（分散σc 2の大小）に追随
して変化させることにより、図９（ａ）に示すようにそ
の震え（分散σc 2）が小さくなった瞬間に血流情報（山
３３に対応する情報）を取出すことができる。
【００３４】以上の実施例に示すように、クラッタ成分
と血流成分の双方の情報を含む受信信号Ｓに基づいてク
ラッタ成分の移動速度Ｖｃ及び／又はその分散σc 2を求
め、この求められた移動速度Ｖｃ及び／又は分散σc 2に
基づいてクラッタ成分の情報を選択的に除去するように
したため、非常に遅い速度の血流であっても、この血流
によるドプラ偏移周波数とクラッタ成分によるドプラ偏
移周波数が互いに異なっている限り、血流成分の情報を
クラッタ成分の情報と分離して取り出すことができる。
と血流成分の双方の情報を含む受信信号Ｓに基づいてク
ラッタ成分の移動速度Ｖｃ及び／又はその分散σc 2を求
め、この求められた移動速度Ｖｃ及び／又は分散σc 2に
基づいてクラッタ成分の情報を選択的に除去するように
したため、非常に遅い速度の血流であっても、この血流
によるドプラ偏移周波数とクラッタ成分によるドプラ偏
移周波数が互いに異なっている限り、血流成分の情報を
クラッタ成分の情報と分離して取り出すことができる。
【００３５】尚、上記各実施例にみるように、クラッタ
成分の移動速度Ｖｃとその分散σc 2の双方を求めること
は必ずしも必要ではなく、クラッタ成分を除去する精度
やその測定の情況等に応じ、移動速度Ｖｃ又は分散σc 2
のいずれか一方のみを求めた場合でもクラッタ成分を選
択的に除去することができる。また、本発明は上記各実
施例に限定されるものではない。例えば上記各実施例で
はアナログ的に直交検波を行った後ディジタル化してい
るが、直交検波を行う前にディジタル化し、その後の演
算で直交検波と等価な処理を行ってもよく、また例えば
被検体内の一点の血流情報を得るいわゆるポイントドプ
ラ法の場合等において、受信信号をディジタル化した後
フーリエ変換し、周波数領域でクラッタ成分の移動速度
Ｖｃやその分散σc 2の検出、クラッタ成分の選択的な除
去、血流情報の検出を行ってもよく、その他種々に構成
されるものである。
成分の移動速度Ｖｃとその分散σc 2の双方を求めること
は必ずしも必要ではなく、クラッタ成分を除去する精度
やその測定の情況等に応じ、移動速度Ｖｃ又は分散σc 2
のいずれか一方のみを求めた場合でもクラッタ成分を選
択的に除去することができる。また、本発明は上記各実
施例に限定されるものではない。例えば上記各実施例で
はアナログ的に直交検波を行った後ディジタル化してい
るが、直交検波を行う前にディジタル化し、その後の演
算で直交検波と等価な処理を行ってもよく、また例えば
被検体内の一点の血流情報を得るいわゆるポイントドプ
ラ法の場合等において、受信信号をディジタル化した後
フーリエ変換し、周波数領域でクラッタ成分の移動速度
Ｖｃやその分散σc 2の検出、クラッタ成分の選択的な除
去、血流情報の検出を行ってもよく、その他種々に構成
されるものである。
【００３６】
【発明の効果】以上詳細に説明したように、本発明の超
音波ドプラ診断装置は、受信信号に基づいてクラッタ成
分の移動速度及び／又は該移動速度の分散を求めるクラ
ッタ情報検出手段と、この求められた移動速度及び／又
は分散に基づいて、受信信号が担持するクラッタ成分の
情報が選択的に除去されたクラッタ除去信号を求めるク
ラッタ情報除去手段と、このクラッタ除去信号に基づい
て血流情報を求める血流情報検出手段とを備えたため、
血流成分のドプラ偏移周波数分布とクラッタ成分のドプ
ラ偏移周波数分布が互いに異なっている限りクラッタ成
分のみを選択的に除去し血流情報を得ることができる。
音波ドプラ診断装置は、受信信号に基づいてクラッタ成
分の移動速度及び／又は該移動速度の分散を求めるクラ
ッタ情報検出手段と、この求められた移動速度及び／又
は分散に基づいて、受信信号が担持するクラッタ成分の
情報が選択的に除去されたクラッタ除去信号を求めるク
ラッタ情報除去手段と、このクラッタ除去信号に基づい
て血流情報を求める血流情報検出手段とを備えたため、
血流成分のドプラ偏移周波数分布とクラッタ成分のドプ
ラ偏移周波数分布が互いに異なっている限りクラッタ成
分のみを選択的に除去し血流情報を得ることができる。
【図１】本発明の超音波ドプラ診断装置の特徴部分の原
理ブロック図である。
理ブロック図である。
【図２】本発明の超音波ドプラ診断装置の第一の実施例
を表わすブロック図である。
を表わすブロック図である。
【図３】図２に示す第一の実施例における複素ＭＴＩフ
ィルタの特性と受信信号との関係を表わした図である。
ィルタの特性と受信信号との関係を表わした図である。
【図４】本発明の超音波ドプラ診断装置の第二の実施例
を表わすブロック図である。
を表わすブロック図である。
【図５】図４に示す第二の実施例におけるＭＴＩフィル
タの特性と受信信号との関係を表わした図である。
タの特性と受信信号との関係を表わした図である。
【図６】位相調整前のキャリア信号と受信信号との位相
関係を表わした模式図である。
関係を表わした模式図である。
【図７】位相調整後のキャリア信号と受信信号との位相
関係を表わした模式図である。
関係を表わした模式図である。
【図８】本発明の超音波ドプラ診断装置の第三の実施例
を表わすブロック図である。
を表わすブロック図である。
【図９】図８に示す第三の実施例におけるＭＴＩフィル
タの特性と受信信号との関係を表わした図である。
タの特性と受信信号との関係を表わした図である。
【図１０】超音波診断装置の一例の概略構成図である。
【図１１】従来の超音波ドプラ診断装置の一例を表わし
たブロック図である。
たブロック図である。
【図１２】直交検波回路の内部構成を表わすブロック図
である。
である。
【図１３】ＭＴＩフィルタの特性の例を表わした図であ
る。
る。
１ クラッタ情報検出手段
２ クラッタ情報除去手段
３ 血流情報検出手段
１７ 直交検波回路
１８ａ，１８ｂ Ａ／Ｄ変換器
１９ 相関器
２０ メモリ
２２ 複素ＭＴＩフィルタ
２３ 相関器
２４ 位相制御回路
２５ａ，２５ｂ ＭＴＩフィルタ
３１ ＭＴＩフィルタ（複素ＭＴＩフィルタを含む）の
特性を表わす折れ線 ３２ クラッタ成分のドプラ偏移周波数分布を表わす山 ３３ 血流成分のドプラ偏移周波数分布を表わす山
特性を表わす折れ線 ３２ クラッタ成分のドプラ偏移周波数分布を表わす山 ３３ 血流成分のドプラ偏移周波数分布を表わす山
Claims (7)
- 【請求項１】 被検体内に送信され該被検体内で反射さ
れた超音波を受信することにより得られた受信信号に基
づいて前記被検体内の血流情報を得る超音波ドプラ診断
装置において、 前記受信信号に基づいてクラッタ成分の移動速度及び／
又は該移動速度の分散を求めるクラッタ情報検出手段
と、 前記クラッタ情報検出手段で求められた前記移動速度及
び／又は前記分散に基づいて、前記受信信号から、該受
信信号が担持するクラッタ成分の情報が選択的に除去さ
れたクラッタ除去信号を求めるクラッタ情報除去手段
と、 前記クラッタ情報除去手段で求められた前記クラッタ除
去信号に基づいて前記血流情報を求める血流情報検出手
段とを備えたことを特徴とする超音波ドプラ診断装置。 - 【請求項２】 前記クラッタ情報検出手段が、少なくと
も前記移動速度を求めるものであり、 前記受信信号を直交検波する直交検波回路を備え、 前記クラッタ情報除去手段が、前記直交検波回路で直交
検波された前記受信信号を入力とし前記クラッタ除去信
号を出力とする複素ＭＴＩフィルタと、前記クラッタ情
報検出手段で求められた前記移動速度に基づいて前記複
素ＭＴＩフィルタの信号遮断帯域の中心周波数を求める
フィルタ制御回路とを備えたことを特徴とする請求項１
記載の超音波ドプラ診断装置。 - 【請求項３】 前記クラッタ情報検出手段が、前記移動
速度とともに前記分散を求めるものであり、 前記フィルタ制御回路が、前記中心周波数を求めるとと
もに、前記クラッタ情報検出手段で求められた前記分散
に基づいて前記複素ＭＴＩフィルタの信号遮断帯域の帯
域幅を求めるものであることを特徴とする請求項２記載
の超音波ドプラ診断装置。 - 【請求項４】 前記フィルタ制御回路が、前記複素ＭＴ
Ｉフィルタに入力された前記受信信号の直流成分が遮断
されるように前記複素ＭＴＩフィルタの特性を定めるも
のであることを特徴とする請求項２又は３記載の超音波
ドプラ診断装置。 - 【請求項５】 前記クラッタ情報検出手段が、少なくと
も前記移動速度を求めるものであり、 前記受信信号を直交検波する直交検波回路と、 前記クラッタ情報検出手段で求められた前記移動速度に
基づいて、前記直交検波回路から出力される前記受信信
号の直流成分がクラッタ成分の移動速度の情報を担持す
るように前記直交検波回路に入力されるキャリア信号の
位相を制御する位相制御回路とを備え、 前記クラッタ情報除去手段が、前記直交検波回路で直交
検波された前記受信信号を入力とし前記クラッタ除去信
号を出力とするＭＴＩフィルタを備えたことを特徴とす
る請求項１記載の超音波ドプラ診断装置。 - 【請求項６】 前記クラッタ情報検出手段が、前記移動
速度とともに前記分散を求めるものであり、 前記クラッタ情報除去手段が、前記クラッタ情報検出手
段で求められた前記分散に基づいて前記ＭＴＩフィルタ
の信号遮断帯域の帯域幅を定めるフィルタ制御回路を備
えたことを特徴とする請求項５記載の超音波ドプラ診断
装置。 - 【請求項７】 前記クラッタ情報検出手段が、少なくと
も前記分散を求めるものであり、 前記受信信号を直交検波する直交検波回路を備え、 前記クラッタ情報除去手段が、前記直交検波回路で直交
検波された前記受信信号を入力とし前記クラッタ除去信
号を出力とするＭＴＩフィルタと、前記クラッタ情報検
出手段で求められた前記分散に基づいて前記ＭＴＩフィ
ルタの信号遮断帯域の帯域幅を求めるフィルタ制御回路
とを備えたことを特徴とする請求項１記載の超音波ドプ
ラ診断装置。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP03161913A JP3093823B2 (ja) | 1991-07-02 | 1991-07-02 | 超音波ドプラ診断装置 |
| US07/904,647 US5383464A (en) | 1991-07-02 | 1992-06-26 | Ultrasonic doppler diagnostic system |
| DE69214023T DE69214023T2 (de) | 1991-07-02 | 1992-07-02 | Filter zur Festechounterdrückung für ein Ultraschalldopplersystem |
| EP92111189A EP0521498B1 (en) | 1991-07-02 | 1992-07-02 | Elutter rejection filter for an ultrasonic doppler system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP03161913A JP3093823B2 (ja) | 1991-07-02 | 1991-07-02 | 超音波ドプラ診断装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH057588A true JPH057588A (ja) | 1993-01-19 |
| JP3093823B2 JP3093823B2 (ja) | 2000-10-03 |
Family
ID=15744406
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP03161913A Expired - Lifetime JP3093823B2 (ja) | 1991-07-02 | 1991-07-02 | 超音波ドプラ診断装置 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US5383464A (ja) |
| EP (1) | EP0521498B1 (ja) |
| JP (1) | JP3093823B2 (ja) |
| DE (1) | DE69214023T2 (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH07178092A (ja) * | 1993-12-24 | 1995-07-18 | Aloka Co Ltd | 超音波ドプラ診断装置 |
| JP2010259790A (ja) * | 2009-04-30 | 2010-11-18 | Medison Co Ltd | 適応クラッターフィルタリングを行う超音波システムおよび方法 |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5349524A (en) * | 1993-01-08 | 1994-09-20 | General Electric Company | Color flow imaging system utilizing a time domain adaptive wall filter |
| US5327894A (en) * | 1993-03-01 | 1994-07-12 | General Electric Company | Wall filter using circular convolution for a color flow imaging system |
| US6029116A (en) * | 1994-08-05 | 2000-02-22 | Acuson Corporation | Method and apparatus for a baseband processor of a receive beamformer system |
| US5544659A (en) * | 1994-12-29 | 1996-08-13 | Siemens Medical Systems, Inc. | Ultrasonic doppler imager having a reduced hardware adaptive tissue rejection filter arrangement |
| US5664575A (en) * | 1994-12-29 | 1997-09-09 | Siemens Medical Systems, Inc. | Ultrasonic doppler imager having an adaptive tissue rejection filter with variable parameters |
| JP3946288B2 (ja) * | 1996-10-01 | 2007-07-18 | 東芝医用システムエンジニアリング株式会社 | 超音波カラードプラ診断装置および超音波カラードプライメージングの信号処理方法 |
| US6039692A (en) * | 1998-06-19 | 2000-03-21 | Vingmed Sound A/S | Method and apparatus for processing ultrasound signals |
| US6296612B1 (en) * | 1999-07-09 | 2001-10-02 | General Electric Company | Method and apparatus for adaptive wall filtering in spectral Doppler ultrasound imaging |
| US6685645B1 (en) | 2001-10-20 | 2004-02-03 | Zonare Medical Systems, Inc. | Broad-beam imaging |
| US6733455B2 (en) * | 1999-08-20 | 2004-05-11 | Zonare Medical Systems, Inc. | System and method for adaptive clutter filtering in ultrasound color flow imaging |
| US7749166B2 (en) * | 2004-04-26 | 2010-07-06 | General Electric Company | System and method for filtering in imaging systems |
| US7643862B2 (en) * | 2005-09-15 | 2010-01-05 | Biomet Manufacturing Corporation | Virtual mouse for use in surgical navigation |
| KR101055581B1 (ko) * | 2008-03-20 | 2011-08-23 | 삼성메디슨 주식회사 | 초음파 시스템 및 클러터 신호 제거용 필터 설정 방법 |
| KR102311699B1 (ko) * | 2020-04-09 | 2021-10-12 | 한화시스템 주식회사 | 레이다 장치 및 이를 이용한 표적 탐지 방법 |
| CN115778434B (zh) * | 2022-12-02 | 2024-12-03 | 业成光电（深圳）有限公司 | 讯号量测方法与都卜勒超音波系统 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4450446A (en) * | 1981-05-29 | 1984-05-22 | Westinghouse Electric Corp. | Method and system for tracking targets in a pulse doppler radar system |
| JPH01110351A (ja) * | 1987-10-23 | 1989-04-27 | Aloka Co Ltd | 超音波ドプラ診断装置 |
| JPH03188841A (ja) * | 1989-09-20 | 1991-08-16 | Toshiba Corp | 超音波診断装置 |
| US5170792A (en) * | 1989-11-27 | 1992-12-15 | Acoustic Imaging Technologies Corporation | Adaptive tissue velocity compensation for ultrasonic Doppler imaging |
| US5107841A (en) * | 1989-11-27 | 1992-04-28 | Acoustic Imaging Technologies Corporation | Maximum entropy velocity estimator for ultrasonic flow imaging system |
| JPH03228752A (ja) * | 1990-02-05 | 1991-10-09 | Hitachi Ltd | 高限界速パルスドプラ計測装置 |
| US5190044A (en) * | 1990-03-30 | 1993-03-02 | Kabushiki Kaisha Toshiba | Ultrasonic blood flow imaging apparatus |
-
1991
- 1991-07-02 JP JP03161913A patent/JP3093823B2/ja not_active Expired - Lifetime
-
1992
- 1992-06-26 US US07/904,647 patent/US5383464A/en not_active Expired - Fee Related
- 1992-07-02 DE DE69214023T patent/DE69214023T2/de not_active Expired - Fee Related
- 1992-07-02 EP EP92111189A patent/EP0521498B1/en not_active Expired - Lifetime
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH07178092A (ja) * | 1993-12-24 | 1995-07-18 | Aloka Co Ltd | 超音波ドプラ診断装置 |
| JP2010259790A (ja) * | 2009-04-30 | 2010-11-18 | Medison Co Ltd | 適応クラッターフィルタリングを行う超音波システムおよび方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| EP0521498A1 (en) | 1993-01-07 |
| EP0521498B1 (en) | 1996-09-25 |
| US5383464A (en) | 1995-01-24 |
| JP3093823B2 (ja) | 2000-10-03 |
| DE69214023T2 (de) | 1997-02-06 |
| DE69214023D1 (de) | 1996-10-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5706818A (en) | Ultrasonic diagnosing apparatus | |
| US5971927A (en) | Ultrasonic diagnostic apparatus for obtaining blood data | |
| US4993417A (en) | Method and system for controlling ultrasound scanning sequence | |
| Evans et al. | Ultrasonic colour Doppler imaging | |
| JP3462584B2 (ja) | 超音波診断装置 | |
| JP3093823B2 (ja) | 超音波ドプラ診断装置 | |
| JP2003501177A (ja) | 同時的組織及びモーション超音波診断撮像 | |
| JPS626810B2 (ja) | ||
| JPH0120899B2 (ja) | ||
| JPH11342128A (ja) | 高フレ―ム速度表示超音波ドップラ―画像処理装置 | |
| EP1354556B1 (en) | Ultrasonic apparatus and method for measuring the velocities of human tissues using the doppler effects | |
| US5501224A (en) | Ultrasonic diagnostic apparatus | |
| US5058593A (en) | Apparatus for processing and displaying ultrasonic data | |
| JPH10165400A (ja) | 超音波診断装置 | |
| KR101809358B1 (ko) | 새로운 평면파 합성을 이용한 초음파 도플러 영상 장치 및 그 제어 방법 | |
| JP2004195228A (ja) | 最適化された送信シーケンスを用いた超音波撮像における適応フィルタリング | |
| JPH03155843A (ja) | 超音波診断装置 | |
| JPH0298344A (ja) | 超音波診断装置 | |
| Gronningsaeter et al. | Blood noise reduction in intravascular ultrasound imaging | |
| JP3281435B2 (ja) | 超音波ドプラ診断装置 | |
| JP7715521B2 (ja) | 超音波診断装置 | |
| US20070073152A1 (en) | Systems and methods for acquiring images simultaneously | |
| Salles et al. | Clutter filter wave imaging: a new way to visualize and detect mechanical waves propagation | |
| JPH03261466A (ja) | 超音波診断装置 | |
| JPH0647043A (ja) | 超音波診断装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20000718 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090728 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110728 Year of fee payment: 11 |
|
| EXPY | Cancellation because of completion of term |