JPS62152437A - 超音波ドプラ診断装置 - Google Patents

超音波ドプラ診断装置

Info

Publication number
JPS62152437A
JPS62152437A JP29211285A JP29211285A JPS62152437A JP S62152437 A JPS62152437 A JP S62152437A JP 29211285 A JP29211285 A JP 29211285A JP 29211285 A JP29211285 A JP 29211285A JP S62152437 A JPS62152437 A JP S62152437A
Authority
JP
Japan
Prior art keywords
signal
velocity
vector
autocorrelation
calculator
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
Application number
JP29211285A
Other languages
English (en)
Other versions
JPH0254738B2 (ja
Inventor
滑川 孝六
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Aloka Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Aloka Co Ltd filed Critical Aloka Co Ltd
Priority to JP29211285A priority Critical patent/JPS62152437A/ja
Priority to US06/947,027 priority patent/US4759375A/en
Priority to CA000526091A priority patent/CA1262958A/en
Priority to DE8686117936T priority patent/DE3681268D1/de
Priority to EP86117936A priority patent/EP0228070B1/en
Publication of JPS62152437A publication Critical patent/JPS62152437A/ja
Publication of JPH0254738B2 publication Critical patent/JPH0254738B2/ja
Granted legal-status Critical Current

Links

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 [産業上の利用分野] 本発明は超音波ドプラ装置、特に運動反射体の速度を演
算し被検体内の動ぎを正確に表示することのできる超音
波ドプラ装置に関する。
[従来の技術] 被検体内の運動部、例えば心臓等の臓器、循環器及び血
管内の血流、体液流又は心筋などの運動反射体の速度を
測定するため、従来より超音波パルストアラ法が実用化
されており、被検体内運動反則体からの反射エコーの周
波数偏移によって運動速度を電気的に検出することがで
きる。
本出願人はこのようなドプラ装置として自己相関による
速度演算装置を、特開昭58−188433にて既に出
願している。
この装置によれば、被検体内から得られたドプラ受信信
号を複索信号に変換し、複索信号の共役積及び複累積を
自己相関法にて求め、この自己相関値から周波数偏移を
演算している。このようにして得られた周波数偏移は運
動反射体の速度に対応しており、この偏移の大きざから
正確な速度が求められる。
そして、これら運動反射体の速度はBモードあるいはM
ll−一ドにてCR7画面−Lに表示され、画像診断に
有効な情報を提供することが行われている。
[発明が解決しようとする問題点] 【1五五五二l直 ところで、超音波ドプラ装置において単一ど一ムにて得
られる速度は超音波ビーム方向に動くラジアル速度のみ
で速度の運動方向を正確に表示することができないとい
う問題がある。
この1=め、従来ではある程度離れた異なる位置から複
数のビームを被検体内の同一位置に放射し、得られた複
数の速度信号を合成して運動方向を含んだ運動反射体の
速度を求めている。
しかしながら、この方法では装置が複雑となるばかりで
なく、特に心臓等のようにビームを挿入する位はや角度
が限定される部位に用いる場合にはこの方法を適用する
ことができないという欠点があった。
発明の目的 本発明は前記従来の問題点に鑑みなされたものであり、
その目的は、異なる位置から超音波ビームをM射づるこ
となく、簡単な走査を行う装置にて運動方向を含む速度
を精iJ:<求めることのできる超音波ドプラ診断装置
を提供することにある。
E問題点を解決するための手段1 前記目的を達成するために、本発明は、一定の繰返し周
波数の超音波を被検体内に放射し送信信号と受信信号と
を比較して運動反射体の距離及び速度を検出する超音波
ドプラ装置において、複索信号変換器と自己相関器とを
協え、まず被検体内から得られたドプラ受信信号を、送
信繰返し周波数の整数倍の周波数を有し互いに複索関係
にある一組の複素基準信号と混合して複素信号に変換り
゛る。
そして、この複素信号の共役積及び複累積を自己相関法
にて求め、この自己相関値から周波数偏移を演算する。
本発明は、同時に任意方向の速度を記憶するメモリと、
接線速度演算器と、を備えている。前記メモリは任意方
向からの第1の受信信号に基づいて得られた第1の自己
相関器出力を記憶し、接線速度演算器は前記第1の受信
信号の受信方向に対して微小の偏向角差を有する第2の
受信信号から得られた第2の自己相関器出力を前記第1
の自己相関器出力と比較してセクタ走査円弧の法線速度
を演算することを特徴とする。
[作用1 以上の構成によれば、メモリに記憶されている第1の自
己相関信号と微小の偏向角差を有する第2の自己相関信
号とが比較されて接線達磨が求められるが、この接線速
度は被検体内の同一距離における速度の変化分を示して
おり、この接線速度とラジアル速度とからベクトル速度
、つまり速度の絶対値及びベクトル角度を求めることが
できる。
そして、このベクトル角度は運動反射体の運動方向を示
すものであり、運動方向を含んだベクトル速1ηは同一
ビーム内の速度分布の全てにおいて、かつ超音波送受波
領域の全てについて求められるので、正確な運動方向に
基づいた運動反射体の速度を画像表示することができる
[実施例] 以下、図面に基づいて本発明の好適な実施例を説明する
第1図、第2図には、超音波ドプラ装2の回路構成が示
されており、まず第2図に基づいてラジアル速度を求め
るための自己相関信号の演算を説明する。
自己相関信号部枠 第2図において、安定な6周波信号を発生する水晶発振
器10のの出力は分周同期回路12に供給され、この分
周同期回路12によって所望周波数の各種出力信号が得
られる。これらの出力信号は超音波パルスビーム送信用
の送信繰返し周波数信号100.複素変換のための複素
基準信号102.104 、超音波診断結果の表示を行
うための部用同期信号106及び装置各部の同期作用を
行うクロック信号108を含む。本発明において、前記
複索基準信号102.104は送信繰返し周波数信号1
00の整数倍の周波数を有しかつ互い   ゛に複索関
係となるようにしており、実施例においては、90°の
位相差を有する。
前記送信信号100は、駆動回路14及び送受切替回路
16を介して探触子18に供給されこの探触子18を励
賑し、超音波パルスビームが被検体20内に送信される
そして、被検体20からの反射エコーは探触子18によ
って電気信号に変換され、送受切替回路16から高周波
増幅器22へ送られて所望の増幅作用が施された後、ぞ
の一方の出力が通常の8モードあるいはMモード表示信
号として表示部に供給される。
通常のBモードあるいはMモード表示を行うための出力
信号は検波器24及びビデオ増幅器26から切替器28
を介してCRT表示器30に供給され、CRT表示器3
0の表示面を輝度変調する。
前記探触子18の超音波パルスビームを機械的あるいは
電気的な角度偏向などによって走査させ超音波パルスビ
ームで被検体20を周期的に走査し、あるいは所望の偏
向角にて走査を停止するために走査制御器32が設けら
番フており、この走査制御器32の走査位置信号及び前
記分周同11J1回路12から117られる部用同期信
号106は掃引トリガ発生器34に供給され、CRT表
示器30の171引制御が行われる。
前記高周波増幅器22の他方の出力は本発明において、
自己相関信号に供され、まず高周波増幅器22から得ら
れるドプラ受信信号は複素信号変換器36に供給され複
ズτ信号に変換される。
すなわち、実施例に、13いては、複素信号変換器36
は位相検波2;を含む一組のミキ膏す38a、38bを
イjし、各ミキサ゛38において前記受信信号は、それ
ぞれ前記複索基準信号102.104と混合され、〜索
基準信号102.104は前述したように互いに90°
位相の異なる複素関係にあるため、ミキサ38から高周
波信号に対応した複索信号を出力することができる。従
って、各ミキサ38は混合検波によって入力された受信
信号と?X/累基準信丹との画周波数の和と差の周波数
の信号を出力し、これら両信号が低域フィルタ40a、
40bに供給され、差の周波数成分のみが取り出される
前記ミキ+138の混合検波作用において、複素基準信
号102.104は単一周波数の連続波であるが、他方
の入力信号である受信信号はドプラ情報を含むパルス波
なので、前記低域フィルタ40の出力には多数のスペク
トル成分が現れることとなる。以下にこの複索変換を演
算式によって説明する。
一方の複索基準信号102は送信用の繰返し周波数fr
の整数倍の周波数fOを有し、その娠幅を1と覆れば、 sin  2 π fo  t           
             ・  (1)なる正弦波電
圧信号にて示される。他方、探触子18で受信される受
信信号は送信周波数をfoとすれば、 sin  (2πfo t、+ 2π■’d t)  
  −(2)にて示される、ただし、fdはドプラ偏移
周波数である。
なお、この受信信号には、一般に 5in(2π(fo±n fr )t+2πfd・(1
±nfr /fo )tl のスペクトル(frは送信繰返し周波数、nはOll、
2・・・なる自然数である)が含まれるが、以下に説明
を簡略化するために、n=oのときの(2)式に示され
るスペクトルについてのみ説明する。
ミキサ38aでは一方の複素基準信号102と受信信号
との積がとられるので、(1)式と(2)式の積の2倍
である次式が得られる。
cos 2πfdt−cos(4πfat+2πfdt
) そして、この出力は低域フィルタ40aで2fo+fd
の周波数が除去されているので、その出力信号は cos  2 π fdl−・・・ (3)となる。
一方、他方の複索1!準信号104は前記信号102と
90’位相が異なるので、 cos 2π1’0’t            ・・
・(4)なる余弦波電圧信号で示され、ミキサ38bの
温合検波及び低域フィルタ40bのフィルタ作用によっ
て、 sin 2yr fd t            −
(5)なる信号に変換され、前記(3)式を実数部、そ
して(5)式を虚数部とする複素信号に変換されたこと
となり、これら両信号は次の複素式によって示すことが
できる。
Zl−CO32πfd t+i sin 2πfdt・
・・(6) 以上のようにして複素変換された信号Z1は△D変換1
42a、42bによってデジタル信号に変換され、次段
の複素ディレーラインキャンセラ44に入力される。そ
して、前記AD変1[42にはクロック信号108が供
給されこのクロック信号によるサンプリングが行われて
いる。
実施例においては、前述した複素ディレーラインキャン
セラ44が設けられているので、被検体内の静止部ある
いは低速運動部からの受信信号を除去して運動部のみの
速度信号を取り出づ“ことができ、画像信号の品質を著
しく向上させることができる。
この複素ディレーラインキャンセラ44は繰返し信号の
1周期(T)に一致する遅延時間を有するディレーライ
ン46a146bを有し、このディレーラインは例えば
1周期の中に含まれるクロックパルスの数に等しい記憶
素子から成るメモリまたはシフトレジスタから形成する
ことができる。
そして、これらディレーライン46には、それぞれ差演
算器48a、48bが接続されており、差演算器48に
よってディレーライン46の入力すなわち現時刻の信号
と1周期前の信号とを同一深度において逐次比較して信
号の1周期間の差を演算する。従って、静止あるい(,
1低速度の生体組織からの縦用(t+ ”J (、l現
It) ’;]の信号と1周期前の信号との間に変化が
なく、あるいは変化が小さいため差rA算器48の差出
力は零に近くなり、一方速度の速い例えば血流信号の差
出力は大きな伯として検出され、これによって生体組織
からの反射信号つまりクラッタを抑制することができる
前記複索ディレーラインキセンセラ44の作用を以下に
演算式で説明する。なお、第2図においては、複素ディ
レーラインキ1!ンセラ44への入力はデジタル信号で
あるが、演算式では説明を簡単にするために、(6)式
のアナログ信号にで説明を行う。ディレーライン46の
入力Z1を(6)式で示すと、1周期遅延された出力Z
2はz2=cos 2πfd  (t−T)+−1si
n2πfd(t−T)  =(7)で示され、この結果
、差部n器48の差出力はZ  =Z  −72=−2
sin 2πfd  (T/2)−sin 2πfd 
 (t −(T/2) )+i 2sin 2πfd 
 (T/2>・cos 2πfd  (t −(T/2
>)となり、ここで差出力Z3を Z= y、3+ i  y3 にて示せば、各x 、y3は次式となる。
X3=−2sin 2zrfd  (T/2)−sin
  2yrfd  (t  −(’r、/2))・・・
 (8) y3=2sin 2πfd  (T/2  )・cos
  2πfd  (t  −(T/2>)・・・ (9
) このようにして、各差演算器48a、48bには、それ
ぞれx 1 y なる信号が出力されることどなる。
以上のようにして低速度信号が除去された複素信号は、
次に自己相関器50にによっC演樟処理され、近延出を
TとするZ3の自己相関が求められる。
まず、入力信号Z3はディレーライン52a、52bに
より1周期分遅延されてZ4が得られる。
この出カフ4は以下の式で表される。
Z  −X4+ +  y4 X4 =−2sin 2πfd  (T/2)・sin
  2πfd  (t  −(3/2T))・・・ (
10) 1/4 =2sin 2πfd  (T/2)−cos
 2πfd  (t −(3/2丁))・・・ (11
) そして、信号Z の共役信号はZ4”=x4−i y 
であり、Z とZ 廁との共役積を、以下の式ににっで
求めることにより自己相関が演算される。
Z3Z4”=(x3+i  y3)(x4−i  y4
)″ x3×4+y3y4 +1(Xa  ’13   Xa  V4)そして、こ
の相関を求めるため、自己相関=50には4個の掛算器
54a、54b、56a156b、そして、加減算器5
8a、58bが設けられ、前記相関演算が行われる。
加減算器58aの出力をRとすれば、前記(8)、(9
)、(10)、(11)の各関係式%式% が得られ、また加減算器58bの出力をIとすれば、同
様に 1=  X  V  −Xa V4=4sin227r
fd(T/2)sin2πfdT   −(13)が得
られ、両加減算器58の出力を合わせて自己相関信号は
次式にて示される。
s=R+r  I           ・・・(14
)そして、この出力Sは信号の変動成分や装置から発生
する雑音成分を含むので、これら雑音成分を除去するた
めに平均回路によって平均が求められ、この平均は’g
−4トi1で表され、複索相関が演鈴される。
前記平均回路はディレーライン608.60bにて1周
期荏延した出力を現時刻の入力信号に加0362a、6
2bにて加算し、再びこの出力をfイー−ライン60に
供給する操作を繰り返し、この加算を、例えばデジタル
回路で構成Jる場合には、その加締出力のL位ピッ1−
を出力すれば、平均値を得ることができる。しかし、単
にこの操作を繰り返し行っていくと、加算回数の増加に
伴い、出力伯が逐次増大し、ついには飽和する。そこで
、実施例にJ3いては、重み付回路64a、64bが設
【プられ、出力を減衰させて入力と加算している。プな
わら、減衰聞をαとすれば、現時刻の信号より例えば1
0周期前の信号はα”だけ減衰して現時刻の信号と加剪
されるので、出力に与える影響度が小さくなり、低域フ
ィルタや移動平均回路と同様の平均機能を果たすことが
可能となる。また、重み付回路64の重み付はを変える
ことにより、平均化の度合いを変更することが可能とな
る。
以上のJこうにして、本実施例においては、複素信号の
共役積を求めることによって自己相関が得られ、この自
己相関、出力は第1図に示すように速度演算器66に供
給され、この速度演算器66によって相関出力Sの偏角
θが求められる。すなわち、偏角0は(12)、(13
)式から、θ−tan−1< T/n > = 2 π
亙T・・・(15) として求められ、この結果、ドプラ偏移周波教祖は、 罹=θ/2π丁         ・・・(16)とし
て前記−角θから極めて容易に求められることとなる。
なお、前記自己相関は複索信号の共役積にて求めている
が、複素信号の複素積にて求めることもでき、この複累
積によれば、近距離にある運動反射体の速度を精度良く
求めることができる。
接線速度及びベクトル速度の演算 本発明において特徴的なことは、微小の偏向角差のある
2個の受信信号の自己相関信号を比較して同一距離にお
ける接線速度を求めることであり、実施例では、この接
線速度からベタ1ヘル速度を演算する。すなわち、前記
自己相関器50等に加えて、比較の対象となる第1の受
信信号の第1の自己相関信号と第2の受信信号の第2の
自己相関信号の1周期分の信号を記憶・読出しするライ
ンメモリ(又はディレーライン)72a、72bと、第
1の自己相関信号と第2の自己相関信号とを比較して接
線速度を演算する接線速度演算器76と、この接線速度
演c1器76の出力に基づいてベクトル速度演算を行う
ベクトル速度演算器84とが設けられている。ここで、
前記ラインメモリは第2の速度分子liと比較するため
に第1の速磨分イ1iを所定時間遅らせるものであり、
このような働きをするものであればラインメモリに限ら
ず各種の遅延線を用いることができる。
また、接線速度演算器76は比較器78、逆正接演算器
80.−1/Δψを乗専する乗篩器82とから成り、前
記比較器78は、実施例では自己相関器50内で共役積
を演算している掛算器54゜56及び加減算器58から
構成され、これらと同様の動作を行う。
本発明は以上の構成から成り、以下に接線速度演算から
ベクトル速度演算までの作用を説明する。
前記自己相関器50の出力であるn、Tは、次式で表わ
され、それぞれラインメモリ72a、72bに供給され
ている。
R= l Rl cosθ        ・・・(1
7)1 = l Rl sin  θ        
   ・(18)このR,Tは自己相関器50の最終値
であり任意方向で数回の超音波パルスを放射して得られ
た信号の平均値である。また、θは速度の信号が含まれ
、超音波ビーム方向と運動方向との成す角度をψとする
と、次式にて表すことができる。
θ+kVcO3ψ         −(19’1第3
図は、Bモード走査における接線速度の求め方を図式化
したものであり、任意方向0→aとこのO→a方向と微
小の偏向角度Δψ離れたO→b方向に超音波パルスを送
受信する場合を考える。
従って、本発明ではO−+a力方向らの反射エコーがパ
ルス列から成る第1の受信信号となり、o−+b力方向
らの反射エコーのパルス列が第2の受信信号となる。
まず、第1の受信信号で得られた第1の自己相関信号の
うちに信号はラインメモリ72aに、1信号はラインメ
モリ72bに記憶され、クロックパルス、走査アドレス
信号を受けたメ[り制御器74にてメモリの書込み・読
出しを行っている。
そして、第2の受信信号にて得られた第2の自己相関信
号はラインメモリ72を経ずに直接比較器78に供給さ
れるとともに、次の計算に備えてラインメモリ72に書
き込まれる。そして、この比較器78では第1と第2の
自己相関信号が比較されるが、例えば第1の自己相関信
号を9a、第2の自己相関信号をSbとすると、IRI
を相関の絶対値として次式にて表わすことができる。
’Sa= l RIcos  (kv cosφ)十j
 lR15in (kv cosψ) ・(20)gb
= l R1cos  [kv cos(ψ−Δφ)1
+j 1R1sin  [kv cos(Φ−Δψ)]
・・・(21) この両信号は比較器78にて、5aとHbの共役積が演
算÷→され、次のにうな接線速度成分の信号を得ること
ができる。
X7= l RI  CO3(kv [cosψ−CO
S  (ψ−Δ〉1) = l Rl 2CO3(−2kv  5in(Φ−Δ
ψ/ 2 ) sinΔψ/2) −11で12cos(kvΔψ sinψ)・・・ (
22) 2  ・ y7=−IRI  Sln  (kvΔψsinψ)・
・・ (23) そして、この信号X  、  V7は逆正接演算器80
で次式にて角度αが演算される。
α=  jan   y7/ X7 =kvΔΦ sinψ        ・(24)この
場合、比較器78と逆正接演算器80との間に、第2図
に示す前記ディレーライン60.加算器629重み付は
回路64とから成る積分器88a、88bを挿入するこ
とにより、雑音信号を除去して精度の更に良い測定を行
うこともできる。
前記角度αは乗算器82にて微小の走査角Δφの逆数の
角の値−1/Δψが型筒される。従って、乗算器82の
出力は、 Vt=kvsinψ         ・(25)とな
り、このVtが接線速度を表すこととなる。
ここで、速度演算器66の出力をVrとすると、vj’
 =kvcosψ         ・(26>となり
、これはラジアル速度を表すものであるから、この接線
速度Vtどラジアル速度Vrとからベク(−ル速I良V
をFtri算することができる。
ベタ1ヘル速瓜演r3器84は次式にてベクトル速度の
絶対値を演C1する。
+ v + =J覆t2+Vr’ =klVl・・・(
27)また、ベタ1−ル角度を次式にて求める。
ψ−tan−’(Vt/Vr)       =−(2
8>前記絶対値IVIは0−1a方向の運動反射体の真
の速度を、またベクトル角は運動方向を表しており、こ
のムク1−ル速度によって運動状況の正確な情報を得る
ことができる。
そして、実施例では、第3図に示されるように、走査角
の中心線OPに対するベクトル角度γに変換しており、
走査角βを前記Φから減算器86にて減算して、次式に
てムク1〜ル角度γを演算している。
γ−ψ−β            ・・・(29)こ
の走査角βは走査制御部32から供給される走査アドレ
スデータである。このアドレスデータを読み出して演算
すれば、ベクトル角度は中心線を基準としたベクトル角
度に変換される。また、この旦準は走査の開始線0−Q
に取ることも可能である。
このようにして、微小の偏向角Δψを順次2Δψ、3Δ
ψ、・・・に変えてそれぞれについて接線速度演算が行
われており、この結果、セクタ画面全体おける運動部の
ベクトル速度が得られる。
このベクトル速度演算器87Iの出力は表示器30に供
給され、ベクトル速度は表示器内のDSCに古き込まれ
C表示器−30上に表示される。
次に、第4図に基づいて接線速度演算器76の他の実施
例を説明する。
この他の実施例は偏向角ΔΦが約1°で前記(22)式
の×7が小さな値になることを考慮して、(23)式の
×7から接線速度を求めている。
すなわち、接線速度演算器76は近似演算器88゜比較
器90及び除算器92を備えており、近似演t5器88
は前記(17)、(18)式から次式を演算する。
−2−2,、、R,2 Rト  ■ また、比較器90は前記比較器78においてV5.V6
を処理する回路で構成される。従って、比較器90の出
力は前記(23)式と同じy7となるがΔψが小さいの
で、次式となる。
y7=−IRI  kv△ψ sinψまた、除篩器9
2は近似病t)器88出力と偏向角Δψとから次式によ
り接線速度Vt′を演算する。
Vt′=−’P2kvΔψsinΦ÷(−R2Δψ)−
kvsinψ このようにして求められた接線速度から前述のようにベ
クトル速度の絶対値とベクトル角とを求めることができ
る。この実施例によれば、演偉精1良を問題にしない場
合に簡略化された装置にでベクトル速度を求めることが
できるという利点を有する。
このような超音波ドプラ3A首はM−し−ドの場合、2
次元にステップ走査するBモードの場合、あるいは移動
目的検出装置(MTr)などに適用でさ、また肢検体内
の断FJ像に重ねて表示Jることができる。
以上のJ:うにして1りられたベクトル速゛度信号はD
A変換器68にJ:ってアナログ電圧信号に変換され、
切替器70を介してCRT表示器30に供給され、C1
マ工表示器30上に輝度g調信号としてB七−ドあるい
はMモードの運動速度分布画像が表示される。
実施例によれば、CRT表示器30は切替器70.28
によりビデオ増幅器26からの通常の画像信号と前述し
たドプラ信号どの両者を選択的に、あるいは同時に表示
可能であり、いずれか任意の画像、あるいはこれら両画
像を重ね合Uた表示を行うことができる。
[発明の効果] 以上説明したように、本発明によれば、微小の偏向角差
を有する2個の受信信号の自己相関信号から接線速度を
求めるようにしたので、この接線速度に基づいてベクト
ル速度を求める口とができ、離れた場所から超音波ビー
ムを放射することなく正確な運動方向に基づいた運動反
射体の速度を容易にVIることができる。
この結果、従来の超音波エコー法にJ:る超嘉波診1f
Ji装置からの6断情報に加えて、血流速亀、血流速度
分布を同時に画像表示することができ、実用上極めて多
くの診断情報をu供できる超音波診断装置を19ること
が可能となる。
【図面の簡単な説明】
第1図は、本発明に係る超音波ドプラ装置の実施例を示
す回路ブロック図、 第2図は、自己相関信号を得るまでの回路構成を示すブ
ロック図、 第3図は、接線速度及びラジアル速度とベクトル速度と
の関係を示す説明図、 第4図は、接線速度部i器の他の実施例の回路を示すブ
ロック図である。 10 ・・・ 水晶発成器 12 ・・・ 分周同期回路 30 ・・・ CRT表示器 32 ・・・ 走査制御器 36 ・・・ 複素信号変換器 38a、38 b  ・−ミキ1) 40a、40b  ・・・ 低域フィルタ44 ・・・
 複素ディレーラインキャンセラ50 ・・・ 自己相
関器 528.52b  ・・・ ディレーライン54a  
、 54b  、 56a  、  56b    −
f卦粋器58a、58b  ・・・ 加減0器 66 ・・・ 速度演算器 72a、72b  ・・・ ラインメモリ76 ・・・
 接線速度演鋒器 78 ・・・ 比較器 80 ・・・ 逆正接aIi算器 82 ・・・ 乗綽器 84 ・・・ ベクトル速;1演C軍器88 ・・・ 
近似演わ器 92 ・・・ 除口器 100  ・・・ 送信繰返し周波数信号102.10
4 ・・・ Nふ1.!準信号106 ・・・ (、+
1引同朋信8 108  ・・・ クロック13月。 第3図 第4図

Claims (6)

    【特許請求の範囲】
  1. (1)一定の繰返し周波数の超音波を被検体内に放射し
    送信信号と受信信号とを比較して運動反射体の距離及び
    速度を検出する超音波ドプラ装置において、受信信号を
    複素信号に変換する複素信号変換器と、送信繰返し周期
    の整数倍の遅れ時間を設けて前記複素信号の自己相関を
    演算する自己相関器と、この自己相関器から出力され任
    意方向から得られた第1の受信信号に対する第1の自己
    相関信号を記憶するメモリと、前記第1の受信信号の受
    信方向に対して微小の偏向角差を有する第2の受信信号
    から得られた第2の自己相関信号を前記第1の自己相関
    信号と比較してセクタ走査円弧の接線方向の接線速度分
    布を演算する接線速度演算器と、を備えたことを特徴と
    する超音波ドプラ診断装置。
  2. (2)特許請求の範囲(1)記載の装置において、前記
    接線速度演算器は自己相関信号の共役積を演算すること
    を特徴とする超音波ドプラ診断装置。
  3. (3)特許請求の範囲(1)又は(2)記載の装置にお
    いて、前記接線速度演算器の出力に基づいて運動反射体
    のベクトル速度を演算するベクトル速度演算器を備えた
    ことを特徴とする超音波ドプラ診断装置。
  4. (4)特許請求の範囲(3)記載の装置において、前記
    ベクトル速度演算器にて求められるベクトル角度は走査
    角により補正してセクタ走査の基準方向に対する角度と
    して演算することを特徴とする超音波ドプラ診断装置。
  5. (5)特許請求の範囲(3)記載の装置において、前記
    ベクトル速度演算器は接線速度演算器出力である接線速
    度の2乗と受信信号から得られたラジアル速度の2乗と
    の和の平方根からベクトル速度の絶対値を演算すること
    を特徴とする超音波ドプラ診断装置。
  6. (6)特許請求の範囲(3)又は(4)又は(5)記載
    の装置において、前記ベクトル速度演算器は接線速度演
    算器出力である接線速度と前記受信信号から得られたラ
    ジアル速度の比の逆正接からベクトル角度を演算するこ
    とを特徴とする超音波ドプラ診断装置。
JP29211285A 1985-12-26 1985-12-26 超音波ドプラ診断装置 Granted JPS62152437A (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
JP29211285A JPS62152437A (ja) 1985-12-26 1985-12-26 超音波ドプラ診断装置
US06/947,027 US4759375A (en) 1985-12-26 1986-12-22 Ultrasonic doppler diagnostic apparatus
CA000526091A CA1262958A (en) 1985-12-26 1986-12-22 Ultrasonic doppler diagnostic apparatus
DE8686117936T DE3681268D1 (de) 1985-12-26 1986-12-23 Mit doppler-effekt arbeitendes ultraschall-diagnosegeraet.
EP86117936A EP0228070B1 (en) 1985-12-26 1986-12-23 Ultrasonic doppler diagnostic apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP29211285A JPS62152437A (ja) 1985-12-26 1985-12-26 超音波ドプラ診断装置

Publications (2)

Publication Number Publication Date
JPS62152437A true JPS62152437A (ja) 1987-07-07
JPH0254738B2 JPH0254738B2 (ja) 1990-11-22

Family

ID=17777701

Family Applications (1)

Application Number Title Priority Date Filing Date
JP29211285A Granted JPS62152437A (ja) 1985-12-26 1985-12-26 超音波ドプラ診断装置

Country Status (1)

Country Link
JP (1) JPS62152437A (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6415033A (en) * 1987-07-09 1989-01-19 Aloka Co Ltd Ultrasonic doppler diagnostic apparatus
JPH02193649A (ja) * 1989-01-24 1990-07-31 Aloka Co Ltd 超音波ドプラ診断装置
JPH03244444A (ja) * 1990-02-21 1991-10-31 Aloka Co Ltd 超音波ドプラ診断装置
WO1992017114A1 (fr) * 1991-04-05 1992-10-15 Yokogawa Medical Systems, Ltd. Dispositif de cartographie d'ecoulement utilisant la correlation croisee d'un signal doppler
US11931202B2 (en) 2018-09-03 2024-03-19 Canon Medical Systems Corporation Ultrasound automatic scanning system, ultrasound diagnostic apparatus, ultrasound scanning support apparatus

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6226051A (ja) * 1985-07-24 1987-02-04 株式会社東芝 超音波血流量自動測定装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6226051A (ja) * 1985-07-24 1987-02-04 株式会社東芝 超音波血流量自動測定装置

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6415033A (en) * 1987-07-09 1989-01-19 Aloka Co Ltd Ultrasonic doppler diagnostic apparatus
JPH02193649A (ja) * 1989-01-24 1990-07-31 Aloka Co Ltd 超音波ドプラ診断装置
JPH03244444A (ja) * 1990-02-21 1991-10-31 Aloka Co Ltd 超音波ドプラ診断装置
WO1992017114A1 (fr) * 1991-04-05 1992-10-15 Yokogawa Medical Systems, Ltd. Dispositif de cartographie d'ecoulement utilisant la correlation croisee d'un signal doppler
US11931202B2 (en) 2018-09-03 2024-03-19 Canon Medical Systems Corporation Ultrasound automatic scanning system, ultrasound diagnostic apparatus, ultrasound scanning support apparatus

Also Published As

Publication number Publication date
JPH0254738B2 (ja) 1990-11-22

Similar Documents

Publication Publication Date Title
US5971927A (en) Ultrasonic diagnostic apparatus for obtaining blood data
US4622977A (en) Ultrasonic diagnostic apparatus
US4961427A (en) Ultrasonic doppler diagnostic apparatus
US4759375A (en) Ultrasonic doppler diagnostic apparatus
JPS6244494B2 (ja)
US4799490A (en) Doppler ultrasonic diagnostic apparatus
US4771789A (en) Ultrasonic device for measuring acceleration of moving reflective member
JPH03162837A (ja) 医用超音波装置
JPS62152437A (ja) 超音波ドプラ診断装置
JP3281435B2 (ja) 超音波ドプラ診断装置
JPH02193650A (ja) 超音波ドプラ診断装置
JPS6096232A (ja) 超音波血流測定装置
JPS62204734A (ja) 超音波ドプラ診断装置
Kasai et al. Real-time two-dimensional blood flow imaging using ultrasound Doppler
JPH02193649A (ja) 超音波ドプラ診断装置
JPH0254737B2 (ja)
JPS6096233A (ja) 超音波血流測定装置
JPS61154545A (ja) 運動反射体の超音波加速度測定装置
JPH03176036A (ja) 超音波ドプラ診断装置
JP2653141B2 (ja) 超音波診断装置
JPS63117739A (ja) 超音波ドプラ診断装置
JPH0479590B2 (ja)
JPH0431264B2 (ja)
JPH0669456B2 (ja) 超音波ドプラ診断装置
JPH03228751A (ja) 超音波血流測定装置

Legal Events

Date Code Title Description
LAPS Cancellation because of no payment of annual fees