JPH0199537A - 超音波診断装置 - Google Patents

超音波診断装置

Info

Publication number
JPH0199537A
JPH0199537A JP25886387A JP25886387A JPH0199537A JP H0199537 A JPH0199537 A JP H0199537A JP 25886387 A JP25886387 A JP 25886387A JP 25886387 A JP25886387 A JP 25886387A JP H0199537 A JPH0199537 A JP H0199537A
Authority
JP
Japan
Prior art keywords
frequency
complex
received signal
signal
velocity
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
JP25886387A
Other languages
English (en)
Other versions
JPH0679601B2 (ja
Inventor
Ikuo Sakai
郁夫 坂井
Yasuhiro Nakamura
恭大 中村
Masami Kawabuchi
川淵 正巳
Nobuaki Furuya
古谷 伸昭
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.)
Panasonic Holdings Corp
Original Assignee
Matsushita Electric Industrial 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 Matsushita Electric Industrial Co Ltd filed Critical Matsushita Electric Industrial Co Ltd
Priority to JP25886387A priority Critical patent/JPH0679601B2/ja
Priority to US07/254,834 priority patent/US4979513A/en
Priority to DE88117037T priority patent/DE3884547T2/de
Priority to EP88117037A priority patent/EP0312059B1/en
Publication of JPH0199537A publication Critical patent/JPH0199537A/ja
Priority to US07/543,890 priority patent/US5127418A/en
Publication of JPH0679601B2 publication Critical patent/JPH0679601B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

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

Description

【発明の詳細な説明】 産業上の利用分野 本発明は、生体内の運動部分の運動速度ベクトル分布を
正確に測定して表示する超音波診断装置に関するもので
ある。
従来の技術 生体内の運動部分の運動速度を測定し、二次元に表示す
ることのできる従来の超音波診断装置は、例えば、特開
昭58−188433号公報に記載の構成が知られてい
る。この方法は超音波のドツプラー効果による受波信号
の位相変化を自己相関関数から求め、運動速度を演算し
、測定都立を微小量ずらしながらこの測定を繰り返す事
により、表示装置に生体内の運動部分の速度分布像を二
次元的に表示している。
発明が解決しようとする問題点 しかしながら、従来の方法では超音波のドツプラー効果
を利用しているため、運動速度の超音波の進行方向の成
分のみの測定で、運動速度の超音波の進行方向に対して
直交する方向の成分を測定できないため、真の運動速度
を測定できないばかりでなく、本来ベクトル量としての
運動速度の方向さえ知ることができない。さらに、例え
ば従来の方法でセクター型の探触子を用いて、セクター
走査させて二次元の運動速度表示させた場合、第7図に
示すように左から右に一様に流れる流体の運動速度の表
示は、図の左側では超音波ビームの進行方向でトランス
ジューサに近ずく方向に5図の右側ではトランスジュー
サから遠ざかる方向に、中央部分ではドソ、プラ効果に
よる周波数シフトを検出することが不可能であるため運
動はまったくしていないよう表示してしまい、実際の運
動とは全くかけ離れた運動速度の表示をしてしまうとい
う問題を有していた。
本発明は従来技術の以上のような問題を解決するもので
、生体内部の運動部分の運動速度の超音波の進行方向に
対して直交する方向の運動速度成分と、運動速度の超音
波の進行方向の運動速度成分を同時に測定し、本来ベク
トル量である生体内の運動部分の運動速度をベクトル量
として二次元ベクトル表示することを可能にする技術を
提供することを目的とするものである。
問題点を解決するための手段 本発明は、複数のチャンネルで同時に受波整相する並列
受波装置と、交互に各チャンネルの受波信号の符号を反
転して加算する加算器と1周波数スペクトラムを求める
第一の周波数分析器と、第一の速度演算器とを具備し、
複素信号変換器と。
自己相関器もしくは第二の周波数分析器と、第二の速度
演算器とを備えることにより、上記目的を達成するもの
である。
作   用 本発明は上記構成により、並列受波回路の各チャンネル
の受波信号の符号を交互に反転させながら加算し、加算
された受波信号の周波数スペクトラムを第一の周波数分
析器で演算する。このことは、生体内の運動部分の超音
波ビームの進行方向に対し直交する方向(チャンネル方
向)の運動によって生じる各チャンネル間の受波信号の
変化分を周波数スペクトラムとして求めることを意味し
、この周波数スペクトラムより生体内の運動部分の超音
波ビームの進行方向に対し直交する方向の運動速度を求
める。と同時忙、受波信号のドツプラ効果による位相の
変化を受波信号の自己相関関数の位相から求めるか、ま
たは第二の周超数分析器で直接周波数分析することによ
り求め、これより生体内部の運動部分の超音波ビームの
進行方向の運動速度成分を演算することにより、生体内
部の運動をベクトル量として測定するこζができる。
実施例 以下1図面を参照しながら本発明の実施例について説明
する。なお、実施例を説明するための全図において、同
一機能を有するものは同一符号をつけその繰り返しの説
明は省略する。
第1図から第3図は1本発明の一実施例を説明するため
の図であり、第1図は、その超音波診断装置の概略構成
を示すブロック図、第2図及び第3図は、並列受波回路
の詳細な構成の一例を示すブロック図である。
第1図に於て、■は超音波ビームを送受するための探触
子であり、第2図に示すように、n個の短冊状の振動子
(以下、エレメントと呼ぶ)を配列状に並べることによ
りトランスジューサを構成したものである。この探触子
lの各エレメント#1〜IL nは、切換回路2に接続
されている。
この切換回路2は、送波時間の間は、n個のエレメント
#1〜#nのうち順次に個(第2図では10個)のエレ
メントを選択し、送波時間の間だlk個のエレメントを
送波回路3に接続し、受波時間の間は、n個のエレメン
ト#1〜#nのうち。
m個(第2図では4個)のエレメントからなる3個(第
2図では4個)の隣あったチャンネルを選択するように
受波増幅器4に接続する。
前記送波回路3では、送波パルスを発生するばかりでな
く、送波パルスの位相側間も行い、前記に個のエレメン
トから送波される超音波ビームを制御する。5〜8は受
波整相回路で、前記j個のチャンネルを構成する各エレ
メントからの受波信号の位相を制御することにより、受
波時における指向性を制御している。9は受波整相回路
で位相整合された受波信号の符号を交互に反転するよう
に加算する加算器、10は基準パルス信号を906位相
シフトする位相シフト器、1)は加算器9で加算された
受波信号と送波信号の基準パルス信号とをミキシングす
るミキサ、 12は90°位相シフトされた基準パルス
信号と受波信号をミキシングするミキサ、13.14は
前記ミキサ1),12の出力をろ波し位相検波信号とし
て出力するローパスフィルタ(LPF)、15.16は
LPFの位相検波信号をデジタル信号に変換するA/D
変換器、17゜18はデジタル信号に変換された位相検
波信号の低周波成分を除去するキャンセラ、19は位相
検波された受波信号の周波数スペクトラムを演算する周
波数分析装置、20は周波数分析装置19で求められた
周波数スペクトラムより、生体内部の運動部分の超音波
ビームの進行方向に対して直交する方向の運動速度成分
を演算する速度演算回路である。21は画像メモリで、
速度演算回路20で求められた運動速度を一時記憶する
。22はD/A変換器、23は切換回路、24は表示装
置。
25は受波整相回路で、Bモード画像を表示するために
受波信号を位相整合する0 26は検波器、27はA/
D変換器、28は画像メモリ、29はD/A変換器、3
0は切換回路である。101は移相器で、送波信号の基
準パルス信号を90’ 位相シフトする。102. 1
03はミキサで、102は受波整相回路27で位相整合
された受波信号をミキシングし、103は移相器101
で90°位相シフトされた基準パルス信号と受波信号を
ミキシングする。104゜105は前記ミキサー102
.103の出力をろ波し位相検波信号として出力するロ
ーパスフィルター(LPF)、106.107はA/D
変換器で、LPFで得られた位相検波信号をデジタル信
号に変換する。
108、 109はキャンセラで、デジタル信号に変換
された位相検波信号の低周波成分を除去する。
1)0、1)1は遅延器で1位相検波信号の一部を遅延
させる。1)2は自己相関器で位相検波された信号から
自己相関関数を演算する。1)3は第2の速度演算器で
自己相関関数から運動部分の超音波ビームの進行方向の
運動速度成分を演算する。
μ上のような構成において、以下その動作を説明する。
第2図は、4チヤンネルが並列に且つ同時に受波し、ま
た、各チャンネルは4エレメントから構成された並列受
波回路の一実施例である。受波時に於て、各チャンネル
は、2エレメントずつずれながら構成され、計10エレ
メントで同時に受波している。即ち受波整相回路5〜8
の出力がそれぞれ各チャンネルの出力に対応している。
ここで、同時に受波するチャンネル数、各チャンネルを
構成するエレメント数、及び隣あったチャンネルのずれ
ピッチは任意に敗ることができ、第3図に示すよ5に、
同時に受波するチャンネル数は4.各チャンネルを構成
するエレメント数を8.隣あったチャンネルのずれピッ
チを1エレメントにしてもよい。第2図に於て、切換回
路2によって選択された各エレメントは対応する受波増
幅器4に接続され、受波増幅器4の出力は受波整相回路
5〜8に入力されるが、受波時の指向性を制菌するため
に、受波増幅器4の出力は受波増幅器40入力となった
エレメントに対応した遅延時間を有する遅延器のタップ
に入力され、各チャンネルで位相制菌されて加算される
。隣あったチャンネルは2工レメント間隔のピッチで並
んでおり、空間的に2工レメント間隔の情報を並列にチ
ャンネルの数だけ、しかも同時に取り込むことができる
一方、送波時においては、少なくとも同時に受波するチ
ャンネルに接続されたすべてのエレメントが送波回路3
に接続され、各エレメントから、送波回路3で位相側−
された送波パルスに応じて指向性の制(財)された超音
波ビームが送波される。
第1図に於て、受波整相回路5〜8から出力された同一
時刻の各チャンネルの受波信号は、加算器によって符号
を反転させながら加算される。加算器で加算された受波
信号の一部は、ミキサ1)で送波信号の基準パルス信号
とミキシングされ、LPF13で位相検波信号に変換さ
れる。また他の一部は、ミキサー12で90’位相シフ
トされた送波信号の基準パルス信号とミキシングされ、
LPF14で位相検波信号に変換され、それぞれA/D
変換器15,16でデジタル信号に変換されるが。
このデジタルに変換された位相検波信号は、互いに90
°位相シフトした関係、即ち、複素共役の関係を持つ複
素位相検波信号である。キャンセラー17.18では、
位相検波信号の中に含まれる生体内組織の体動にともな
う低周波数成分、いわゆるクラッタ−成分を除去する。
周波数分析装置19に入力された。キャンセラー17.
18でクラッタ−成分の除去された複素位相検波信号の
周波数スペクトラムが演算される。この周波数分析の演
算において、生体内の運動にともなった超音波ビームの
各チャンネルにおける受波信号をf(x)、受波信号の
符号を反転させながら加算することを示す関数をh (
xlと表現する。h(x)は第5図(a)に示すように
並列受波装置のチャンネルピッチに等しい間隔で矩形波
状に変化する関数である。加算器において加算された受
波信号をg (X)とし、さらにGtllをg tX)
のフーリエ変換、またH(kl、 F(k)はそれぞれ
)1(XI、f(x)のフーリエ変換とすると以下の関
係が成立する。
ここでVは測定しようとしている生体内の運動部分の超
音波ビームの進行方向に対し直交する方向の運動速度で
あり、また、HlklはH(k)の複素共役を表してお
り、第5図(b)に示すようにに=2π/p、(pはチ
ャンネル間隔のピッチ)で鋭いピークを示す関数である
。一方、F(k)はブロードな関数であるから、その積
で表されるGK)は鋭いピークを示す関数となる。また
、ここで、周波数スペクトラムは、n個の超音波ビーム
のパルス繰り返し周期Tの時間間隔の複素位相検波信号
で演算されているが、この個数iは%Bモード画像を構
成するための走査線数N1超音波ビームのパルス繰り返
し周期T、及び画像のフレームレートFからに次の関係
式 %式%(2) によって決定される。速度演算器20において、周波数
スペクトラムのピークをあたえるωPはこの周波数スペ
クトラムより演算されるが、(2)式の制限により、超
音波診断画像にリアルタイム性を持たせるために、lを
十分に大きく取れず、周波数分解能が粗くなるため、周
波数スペクトラムモーメントとして求める。
ωP=ω・ IG(ω) l         (31
このωPから、生体内の運動部分の超音波ビームの進行
方向に対し直交する方向の運動速度成分Vは。
で求められ、生体運動部分の運動速度をベクトル量とし
て測定1表示が可能となる。
以上の説明から明らかなように本実施例によれば、隣合
う複数のチャンネルで並列に、しかも同時に受波し、加
算器で、各チャンネルの受波信号を交互に符号を反転し
ながら加算し、周波数分析器で求めた周波数スペクトラ
ムから、生体内運動部分の超音波ビームの進行方向と直
交する方向の運動速度成分を求めることができる。一方
、受波整相回路29で位相整合された受波信号の一部は
、ミキサー102で送波信号の基準パルス信号とミキシ
ングされ、 LPF104で位相検波信号に変換される
。また他の一部は、ミキサー103で90°位相シフト
された送波信号の基準パルス信号とミキシングされ、L
PFt05で位相検波信号に変換され、それぞれA/D
変換器106.107でデジタル信号に変換されるが、
このデジタルに変換された位相検波信号は、互いに90
°位相シフトした関係、即ち、複素共役の関係を持つ複
素位相検波信号である。
キャンセラー108 、109では、クラッタ−成分を
除去する。キャンセラ108.109でクラッタ−成分
の除去された複素位相検波信号の一部は、遅延器1)0
.1)1で遅延され、遅延された複素位相検波信号と、
遅延されていない複素位相検波信号が入力された自己相
関器1)2では複素量としての自己相関関数が求められ
る。いま、自己相関関数の実数部をCr%虚数部をCi
とすると、生体内の運動部分での超音波ビームのドツプ
ラ効果による周波数シフ) fdとの間には次の関係が
ある。
ここで、Tは超音波ビームのパルス繰り返し周期である
。ドツプラー効果は、生体内部の運動部分の超音波ビー
ムの進行方向の運動速度成分による超音波ビームの周波
数シフトであるから、速度演算器1)3において、自己
相関器1)2で計算された自己相関関数より(4)式の
関係よりfdを求め、これより生体内部の運動部分の超
音波ビームの進行方向の運動速度成分を求めることがで
き、ベクトル情報としての生体内部の運動速度成分を求
めることができる。
以上本実施例によれば、生体内部の運動部分の超音波ビ
ームの進行方向と、それに直交する方向の運動速度成分
が同時に測定でき、生体内部の運動部分の状態をベクト
ルとして測定表示ができる。
次に本発明の第2の実施例について説明する。
第4図は本発明の第2の実施例における超音波診断装置
の概略を示すブロック因である。
第4図において第1図と異なる点は、遅延器1)0゜1
)1および自己相関器1)20代わりに周波数分析器2
01を備えている点である。この周波数分析器2旧にお
いて、複素位相検波信号を直接周波数分析することによ
り、生体内運動部分による超音波ビームのドツプラー効
果による周波数シフトを直接演算し、生体内の運動部分
の超音波ビームの進行方向の運動速度成分を求めている
発明の効果 以上のように本発明は、複数のチャンネルで同時に受波
整相する並列受波整相回路と、各チャンネル間の受波信
号を交互に反転させながら加算する加算器と、加算され
た受波信号を複素位相検波信号に変換する回路と、周波
数スペクトラムを演算する周波数分析器と、周波数スペ
クトラムより生体内の運動部分の速度を演算する速度演
算器を具備することにより、生体内運動部分の超音波パ
ルスビームの進行する方向の運動速度成分を求めること
ができ、従来の装置に於て、測定、表示できなかった正
確な生体内の運動部分の状態をベクトル情報として測定
1表示することが可能となり、非常に分かりやすい表示
で、正確な診断を行うことができる超音波診断装置を提
供することができ。
その効果は大きい。
【図面の簡単な説明】
第1図は本発明の第1の実施例における超音波診断装置
の概略を示すブロック図、第2図は本発明の第1の実施
例における並列受波回路の一実施例の詳細を示すブロッ
ク図、第3図は本発明の第1の実施例における並列受波
回路の別の実施例の詳細を示すブロック図、第4図は本
発明の第2の実施例における超音波診断装置の概略を示
すブロック図、第5図(alは加算器の特性を表わす線
図、第5図tb)は加算器の特性を周波数領域で表した
線図%第6図は従来の超音波診断装置の問題点を説明す
るための図である。 1・・・探触子、2・・・切替回路、3・・・送波回路
、41.、受波増幅器、5〜8・・・受波整相回路、9
・・・加算器、10・・・90’立相器、1).12・
・ミキサ、13.14・・・LPF、15.16・・・
A/D変換器、  17.18・・・キャンセラ、19
・・・周波数分析器、20・・・速度演算回路、21・
・・画像メモリ、22・・・D/A変換器、23・・・
切換回路、24・・・表示装置、25・・・受波整相回
路、26・・検波器、27・・・A/D変換器、28・
・・画像メモリ、29・・D/A変換器、30・・・切
換回路、101・・・90°位相器、  102,10
3・・・ミキサ、104.105・・LPF。 106.107・・・A/D変換器k  108,10
9・・・キャンセラ。 1)0、1)1・・・遅延器、1)2・・・自己相関器
、1)3・・速度演算回路、201・・・周波数分析器
。 代理人の氏名 弁理士 中 尾 敏 男 ほか1名第 
2 図 第3図 第5図 h  (x) aπ 第6rW

Claims (2)

    【特許請求の範囲】
  1. (1)超音波パルスビームを一定の繰り返し周期で生体
    内に送波し、その反射波を受波し、この受波信号を増幅
    し、この増幅された受波信号を複数のチャンネルで同時
    に受波整相する並列受波回路と、前記並列受波回路の各
    チャンネルの受波信号を交互に符号を反転させて加算す
    る加算器と、送波繰り返し周波数の整数倍の周波数を有
    し互いに複素共役関係にある一組の複素基準信号と前記
    加算器に於て加算された受波信号とを混合して、受波信
    号を複素信号に変換する第一の複素信号変換器と、複素
    信号に変換された受波信号の周波数スペクトラムを演算
    する周波数分析器と、周波数スペクトラムから生体内の
    運動部分の超音波パルスビームの進行方向と直交する方
    向の運動速度を演算する第一の速度演算器と、送波繰り
    返し周波数の整数倍の周波数を有し互いに複素共役関係
    にある一組の複素基準信号と前記増幅された受波信号と
    を混合して、受波信号を複素信号に変換する第二の複素
    信号変換器と、前記複素信号の遅れ時間を設けて複素信
    号の自己相関関数を演算する自己相関器と、前記自己相
    関器の出力から生体運動の超音波パルスビームの進行方
    向の速度を演算する第二の速度演算器とを具備し、超音
    波パルスビームの進行方向及びそれと直交する方向の速
    度ベクトル分布の測定及び表示することを特徴とする超
    音波診断装置。
  2. (2)超音波パルスビームを一定の繰り返し周期で生体
    内に送波し、その反射波を受波し、この受波信号を増幅
    し、この増幅された受波信号を複数のチャンネルで同時
    に受波整相する並列受波回路と、前記並列受波回路の各
    チャンネルの受波信号を交互に符号を反転させて加算す
    る加算器と、送波繰り返し周波数の整数倍の周波数を有
    し互いに複素共役関係にある一組の複素基準信号と前記
    加算器に於て加算された受波信号とを混合して、受波信
    号を複素信号に変換する第一の複素信号変換器と、複素
    信号に変換された受波信号の周波数スペクトラムを演算
    する第一の周波数分析器と、周波数スペクトラムから生
    体内の運動部分の超音波パルスビームの進行方向と直交
    する方向の運動速度を演算する第一の速度演算器と、送
    波繰り返し周波数の整数倍の周波数を有し互いに複素共
    役関係にある一組の複素基準信号と前記増幅された受波
    信号とを混合して、受波信号を複素信号に変換する第二
    の複素信号変換器と、前記複素信号の遅れ時間を設けて
    複素信号の周波数スペクトラムを演算する第二の周波数
    分析器と、前記第二の周波数分析器の出力から生体運動
    の超音波パルスビームの進行方向の速度を演算する第二
    の速度演算器とを具備し、超音波パルスビームの進行方
    向及びそれと直交する方向の速度ベクトル分布の測定及
    び表示することを特徴とする超音波診断装置。
JP25886387A 1987-10-14 1987-10-14 超音波診断装置 Expired - Fee Related JPH0679601B2 (ja)

Priority Applications (5)

Application Number Priority Date Filing Date Title
JP25886387A JPH0679601B2 (ja) 1987-10-14 1987-10-14 超音波診断装置
US07/254,834 US4979513A (en) 1987-10-14 1988-10-07 Ultrasonic diagnostic apparatus
DE88117037T DE3884547T2 (de) 1987-10-14 1988-10-13 Ultraschalldiagnostikgerät.
EP88117037A EP0312059B1 (en) 1987-10-14 1988-10-13 Ultrasonic diagnostic apparatus
US07/543,890 US5127418A (en) 1987-10-14 1990-06-27 Ultrasonic diagnostic apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP25886387A JPH0679601B2 (ja) 1987-10-14 1987-10-14 超音波診断装置

Publications (2)

Publication Number Publication Date
JPH0199537A true JPH0199537A (ja) 1989-04-18
JPH0679601B2 JPH0679601B2 (ja) 1994-10-12

Family

ID=17326080

Family Applications (1)

Application Number Title Priority Date Filing Date
JP25886387A Expired - Fee Related JPH0679601B2 (ja) 1987-10-14 1987-10-14 超音波診断装置

Country Status (1)

Country Link
JP (1) JPH0679601B2 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10414050B2 (en) 2015-03-31 2019-09-17 Canon Kabushiki Kaisha Automated assembly method using automated assembly apparatus

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10414050B2 (en) 2015-03-31 2019-09-17 Canon Kabushiki Kaisha Automated assembly method using automated assembly apparatus

Also Published As

Publication number Publication date
JPH0679601B2 (ja) 1994-10-12

Similar Documents

Publication Publication Date Title
US4573477A (en) Ultrasonic diagnostic apparatus
EP0144968B1 (en) Ultrasonic diagnostic apparatus
US4799490A (en) Doppler ultrasonic diagnostic apparatus
JPH0614930B2 (ja) 超音波診断装置
EP0475100B1 (en) Ultrasonic doppler blood flowmeter
EP0312059B1 (en) Ultrasonic diagnostic apparatus
JPH0199537A (ja) 超音波診断装置
JPH09299369A (ja) 超音波連続波ドプラ血流計
JPH0199542A (ja) 超音波診断装置
JPH0576302B2 (ja)
JPH0199539A (ja) 超音波診断装置
JPH0199538A (ja) 超音波診断装置
JPH0199541A (ja) 超音波診断装置
US5127418A (en) Ultrasonic diagnostic apparatus
JP2594994B2 (ja) パルスドプラ計測装置
JP3391578B2 (ja) 相関装置および流れ情報表示装置
JPH0332653A (ja) 超音波診断装置
JPH01314551A (ja) 超音波診断装置
JPS6058131A (ja) 超音波ドプラ装置
JP2574820B2 (ja) 超音波パレスドプラ計測装置
JP2714042B2 (ja) パルスドプラ計測装置
JPH04197248A (ja) パルスドプラ計測装置
JP2714067B2 (ja) パルスドプラ計測装置
JPS62204734A (ja) 超音波ドプラ診断装置
JPH0199540A (ja) 超音波診断装置

Legal Events

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