JPH02152447A - カラー超音波診断装置 - Google Patents

カラー超音波診断装置

Info

Publication number
JPH02152447A
JPH02152447A JP30754688A JP30754688A JPH02152447A JP H02152447 A JPH02152447 A JP H02152447A JP 30754688 A JP30754688 A JP 30754688A JP 30754688 A JP30754688 A JP 30754688A JP H02152447 A JPH02152447 A JP H02152447A
Authority
JP
Japan
Prior art keywords
blood flow
distribution image
frequency
color
data
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
JP30754688A
Other languages
English (en)
Other versions
JPH0716489B2 (ja
Inventor
Tatsuro Baba
達朗 馬場
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.)
Toshiba Corp
Original Assignee
Toshiba Corp
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 Toshiba Corp filed Critical Toshiba Corp
Priority to JP30754688A priority Critical patent/JPH0716489B2/ja
Publication of JPH02152447A publication Critical patent/JPH02152447A/ja
Publication of JPH0716489B2 publication Critical patent/JPH0716489B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

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

Description

【発明の詳細な説明】 [発明の目的] (産業上の利用分野) 本発明は、超音波のドプラ効果を利用して被検体内の血
流情報を求め、これを2次元表示するカラー超音波診断
装置に関する。
(従来の技術) 超音波ドプラ法とパルス反射法とを併用することによっ
て一つの超音波プローブで血流情報と断層像(Bモード
像)情報を得、断層像に重ねて血流情taをリアルタイ
ムでカラー表示するようにしたカラー超音波診断装置が
知られている。このような装置によって血流情報を1q
る場合の動作原理は次の通りでおる。
すなわち、被検体である生体内を流れている血流に対し
て超音波パルスを送波すると、この超音波ビームの中心
周波数fcは流動する血球によって散乱されドプラ偏移
を受けて周波数fdだけ変化して、この受波周波数fは
f=fc +fdとなる。このとき周波数fc、fdは
次式のように示される。
に こで、■=平均血流速 θ:超音波ビームと血管とのなす角度 C:音速 従って、ドプラ偏移f(lを検出することによって平均
血流速Vを得ることができる。
このようにして得られた平均血流速Vの2次元画像表示
は次のように行われる。先ず第14図のように超音波プ
ローブ1から被検体に対してA。
B、C,・・・方向に順次超音波パルスを送波してセク
タく又はリニア)スキャンを行うにあたり、その超音波
パルスのスキャン制御が行われる。
最初に六方向に数回超音波パルスが送波されると、被検
体内の血流でドプラ偏移されて反射されたエコー信号は
同一プローブ1によって受波され、電気信号に変換され
て受信回路に送られる。
次に位相検波回路によってドプラ偏移信号が検出される
。このドプラ偏移信号は超音波パルスの送波方向に設定
された例えば256個のサンプル点ごとにとらえられる
。各サンプル点でとらえられたドプラ偏移信号は周波数
分析器で周波数分析され、フレームメモリを備えたり、
S、C,(ディジタル・スキャン・コンバータ)に送ら
れここで走査変換された後に、表示部に送出されA方向
の血流分布像が2次元画像としてリアルタイムで表示さ
れる。
以下B、C,・・・の各方向に対しても同様な動作が繰
り返されて、各スキャン方向に対応した血流分布像(カ
ラーフローマツピング(CFM)像とも称される)が表
示されることになる。
第15図は平均血流速の角度表現1周波数表現とそれに
対応する表示色との関係を示すもので、角度表現の場合
には一π乃至十πに対応して、周波数表現の場合には+
fr/2乃至−fr/2に対応して、表示色が青色、黒
色、赤゛色となる。尚、frは超音波パルスの繰返し周
波数である。
ところで、従来装置において、CFM像の時間的平滑化
(残光効果)処理が行われ、この時間的平滑化処理には
フレーム相関が用いられている。
このフレーム相関による時間的平滑化処理は、第9図に
示すように1フレーム構成周明(T OF >毎に連続
する複数のCFM像が形成される場合において、同一座
標(x、y)のデータの相関を取る手法である。原理的
には第10図に示すように巡回型のディジタルフィルタ
であり、物理的には時間軸(1)のローパスフィルタで
ある。ここでCFM@のフレーム周波数をf OFで表
わし、隣接フレーム間でのデータ変化をfで表わした場
合、f / f OF≧172が成立すると、階調変化
の所謂「折返し現象」を生じてしまう(第11図参照〉
(発明が解決しようとする課題) Bモード像データ(白黒エコーデータ)のような符号無
しデータの相関による時間的平滑化処理は、第12図の
ように何ら問題無く行われる。
しかしCFM像では第13図において実線2で示すにう
に平滑化され、ドプラ周波数のパルス繰返し周期による
折返し現象により赤→青、青→赤へ移行する際に周波数
ゼロラインと交差するように血流速度V (t)が変化
する。このぜ0交差部分3に対して黒色が割当てられる
ため、そこでCFM像をフリーズ(静止画像表示)する
と、実際には高速血流部位が存在するのにもかかわらず
、表示画像上では当該部位が黒色となり血流が無いよう
に表示されてしまう。
本発明は上記欠点を除去するもので、その目的とすると
ころは、血流分布像のフレーム相関による時間的平滑化
処理において折返し現象を生じた場合の高血流部位の黒
色表示を排除することにより、適切な血流分布像表示を
可能とするカラー超音波診断装置を提供することにある
[発明の構成1 (課題を解決するための手段) 本発明は、被検体よりの超音波エコーに基づいて形成さ
れた血流分布像に対してフレーム相関による時間的平滑
化処理を施し、その処理結果をカラー表示するようにし
たカラー超音波診断装置において、血流分布像のフレー
ム周波数をf OFとし、隣接するフレーム間での血流
分布像データの変化をfとしたとき、f / f OF
≧172の場合に血流情報の周波数表現上の周波数ぜロ
ラインと交差ジノない方向に血流分布像のフレーム相関
による時間的平滑化処理を行う平滑化処理手段を具備す
るものである。
(作 用) 本発明ではf / f OF≧172の場合に血流情報
の周波数表現上の周波数ゼロラインと交差しない方向に
時間的平滑化処理を行うようにしているので、折返し現
象が生じても周波数ゼロラインと交差する部分は存在せ
ず、従って、折返し現象を生じた場合の高血流速部位の
黒色表示を排除できる。
(実施例) 以下、本発明を実施例により具体的に説明する。
第1図は本発明の一実施例を示している。
10は複数の超音波振動子をアレイ状に配列して成る超
音波プローブで必り、11はこの超音波プローブ10を
介して超音波パルスの送受信を行う送受信部であり、1
2は超音波エコーの包絡線検波を行い、それをディジタ
ル信号に変換して出力するレシーバであり、13はこの
レシーバ13の出力が書込まれる白/黒フレームメモリ
である。
また、14は超音波エコーの直交検波を行う直交検波部
であり、15はCFM(カラーフローマツピング)処理
を行うCFMユニットであり、16はこのCFM処理出
力が書込まれるカラーフレームメモリである。このカラ
ーフレームメモリ16及び白/黒フレームメモリ13の
後段には、両メモリからの読出しデータを合成する合成
処理部17が配置され、更にこの合成処理部17の後段
には、この合成処理出力をカラー表示するRGBモニタ
18が配置されている。
ここで前記CFMユニット15は次のように構成されて
いる。
19はA/D (アナログ・ディジタル)変換器であり
、前記直交検波部14の出力がこのA/D変換器19に
よりディジタル信号に変換されるようになっている。こ
のA/D変換器19の後段には、クラッタ成分を除去す
るMTIフィルタ20が配置され、またこのフィルタ2
0の後段には、平均血流速(V)、トータルパワー(P
)9分散(写)の算出を行うCFM演算部21が配置さ
れ、更にこのCFM演算部21の後段には、フレーム相
関により血流分布像の時間的平滑化処理を行う平滑化処
理手段22が配置されている。ここでこの平滑化処理手
段22においては、血流分布像のフレーム周波数をf 
OFとし、隣接するフレーム間での血流分布像データの
変化をfとしたとき、f / f OF≧172 の成立により折返しを生じた場合に血流情報の周波数表
現上の周波数ゼロラインと交差しない方向に血流分布像
の時間的平滑化処理が行われるようになっており、これ
が本実施例装置の特徴点の一つとなっている。
次に、この平滑化処理手段22の機能を実現する。ため
の構成例について説明する。
第2図は平滑化処理手段22の詳細な構成を示している
25は第1図のCFM演算部21からの血流速信号V 
(t)を(1−a)倍する係数ROM(+、J−ドオン
リメモリ)であり、この係数ROM25の出力は、後段
に配置された加算器26に取込まれ、この加算器26に
おいて、後述する係数ROM29の出力と加算された後
にMPX (マルチプレクサ)27を介して画像メモリ
28に取込まれるようになっている。ここでaはフレー
ム相関処理における係数を意味する。MPX27は加算
器26の出力と初期条件ROM30の出力とを択一的に
画像メモリ28に伝達する機能を有し、いずれを選択す
るかは初期条件ROM30より選択信号SELによる。
初期条件ROM3.0は画像メモリ28の出力と入力信
号V([)との比較を行い、その比較結果に応じてMP
X27に選択信号SELを送出すると共に初期条件を送
出する。尚、この初期条件及びMPX27の選択論理に
ついては後に詳述する。
また画像メモリ28は第5図に示すように方位方向にM
ラスタ分、距離方向にN画素分の記憶容量を持つ2次元
構造のメモリであり、この画仰メモリ28へのデータ書
込み動作とこの画像メモリ28らのデータ読出し動作と
の切換え、及びアドレス指定はメモリコントローラ31
によって行われる。そしてこのメモリコントローラ31
には、1フレーム構成周期信号OFO1CMFスキャン
ラスタ先頭信号INITO、CFM演算部からの転送ク
ロックDEPIXCにが入力されるようになっている。
更にこの画像メモリ28の後段には、この画像メモリ2
8の出力を8倍する係数ROM29が配置されている。
そして係数ROM25.29には第1図のホストCPt
J (中央処理装置)23からの係数選択信号が取込ま
れるようになっている。この信号によって係数aの値が
決定される。すなわち係数ROM25.29における係
数aの値はホストCPU23の制御下で決定されるので
ある。
またMPX27の出力がこの平滑化処理手段22の出力
として第1図のカラーフレームメモリ16に送出される
ようになっている。
第8図は前記MPX27の選択論理を示しており、以下
、これについて説明する。
現在のデータをV(n)とし、このV (n)より1フ
レーム前のデータをV (n+1)とする。このV(n
) 、 V(n+1)の符号が等しい場合、初期条件R
OM30より選択信号SELは出力されず、MPX27
は加算器26の出力を選択する(フレーム相関処理継続
)。また、V(n) 、 V(n+1)の符号が異なる
場合であって、それらの差の絶対値がスレッショルドレ
ベルTH以上である場合には、フレーム相関処理が一旦
停止され、初期条件が与えられる。すなわち、この場合
には初期条件ROM30より選択信号SEL及び初期条
件が出力され、MPX27は初期条件ROM30よりの
初期条件を選択してそれを画像メモリ2B及びカラーフ
レームメモリ16に送出する。更に、V(n) 、 V
(n+1)の符号が異なる場合であって、それらの差の
絶対値がスレッショルドレベルTHより小ざい場合には
、初期条件ROM30より選択信号SELは出力されず
、MPX27は加算器26の出力を選択する(フレーム
相関処理継続)。
ここで、スレッショルドレベルTHは通常fr /2 
(frは超音波パルスの繰返し周波数)に設定されてい
るが、Oからfrまで必要に応じて調整できるようにな
っている。また、初期条件は、V (n+1)が正の場
合には+fr /2であり、逆にV (n+1)が負の
場合には−fr /2である。
次に、上記構成による実施例装置の作用について説明す
る。
超音波プローブ10から被検体に向けて送波された超音
波パルスの該被検体よりの反射成分は再び超音波プロー
ブ10によって受波され、それが送受信部11を介して
レシーバ12に取込まれた包絡線快波され、それがディ
ジタル信号に変換された後に白/黒フレームメモリ13
に転送され、ここに書込まれる。この白/黒フレームメ
モリ13の記憶内容はモニタ系のタイミングで続出され
、合成処理部17を介してRGBモニタ18に伝達され
る。
一方、ドプラ情報を得るための超音波パルスの送受によ
って得られたパルスエコーは直交検波部14において直
交検波され、その検波出力がCFMユニット15に送出
される。そして上記の直交検波出力はA/D変換器19
によりディジタル信号に変換され、MTIフィルタ20
によりクラッタ成分が除去された後にCFM@R部21
に取込まれ、平均血流速(V)、トータルパワー(P)
2分散(σ)の算出に供される。そしてこのCFM演算
部21での算出結果が平滑化処理手段22に取込まれ、
この平滑化処理手段22においてフレーム相関による時
間的平滑化処理が行われ、その処理結果がカラーフレー
ムメモリ16及び合成処理部17を介してRGBモニタ
18に伝達され、白黒Bモード像に重畳して表示される
ここでこの時間的平滑化処理の詳細について説明する。
フレーム相関による時間的平滑化処理は第9図に示すよ
うにCFM像の1フレーム構成周期(OF周期毎)のフ
レーム間での処理でおる。
第3図及び第4図はCFM出力タイミングを示しており
、第3図は方位方向のデータ出力タイミング図、第4図
は距離方向のデータ出力タイミング図である。
CFM1フレーム構成周期(OFO信号周期)中にCF
Mスキャンラスク先頭信号(INITO)がM(Mは正
の整数)個発生する。またCFM像の1ラスタのデータ
を得るためにレー1〜周波数(PRF)周期で同一ラス
タを数回スキャンしている。このようにして得られるC
FM出力信号V (t)は第3図中の斜線部のタイミン
グで出力される。またCFM1ラスタ分の出力信号V(
t)における距離方向のデータ1’   2’   3
’・・・。
N′は第4図に示すようにル−ト内に順次出力される。
第6図は画素(m、n)のフレーム相関処理タイミング
を示し、第7図は第6図のTa部分を拡大して示してい
る。
第6図に示すように画像メモリ28のアドレスはメモリ
RD信号に同期してインクリメントされ、(m、1)、
(m、2)、−・・(m、n−1)、(m、n)、・・
・のように変化する。このメモリアドレスの確定後、そ
れまで画像メモリ28に蓄えられた1叶周期前のデータ
が読出され、それが現在取込まれたデータと合わせて処
理され、この処理結果がメモリWR信号により再び画像
メモリ28に書込まれる。
この画素(m、n)についての相関処理の詳細は次の通
りである。
すなわち第7図に示すように画像メモリ28よりアドレ
ス(m、n)の旧データが読出され、それが係数ROM
29で8倍されて加算器26に取込まれ、又現在の入力
データV (t)が係数ROM25で(a−1)倍され
て加算機26に取込まれると、両データがこの加算器2
6で加算され、そしてその加算結果がアドレス(m、n
)の新たなデータとしてMPX27を介して画像メモリ
28に書込まれる(データ更新)。そしてこの占込みと
ほぼ同時にfr(m、n)データがカラーフレームメモ
リ16(第1図)に転送される。このような処理がCF
M像の全画素について行われる。
ここで、画像メモリ28より読出されたデータ(旧デー
タ)は係数ROM29に転送されるとともに初期条件R
OM30に伝達される。このデータをV(n+1)とす
る。そしてこのとき、新たに取込まれたデータV(t)
  (これをV(n)とする)と旧データV (n+1
)とが同符号の場合及び、それらが異符号であって、し
かも が成立する場合には、MPX27により加算器26の出
力が選択されるため、上記のフレーム相関処理が継続さ
れる(第8図参照)。しかし、V (n)  、 V 
(n+1)が異符号であって、しかもが成立する場合に
は、初期条件ROM30より選択信号SFLとともに初
期条件が出力され、MPX27は、この選択信号SEI
入力により、初期条件ROM30の初期条件を選択する
(第8図参照)。ここで、初期条件はV(n+1)が正
の場合は+fr/2であり、逆にV(n+1)が負の場
合に−fr/2である。すなわち前記加算器26の出力
に代えてそれぞれの血流方向の最高階調に等しい値が画
像メモリ28及びカラ」フレームメモリ16に与えられ
るのである。
このような処理によれば、f / f OF≧172に
より折返し現像を生じた場合、血流情報の周波数表現上
の周波数ゼロラインと交差しない方向にフレーム相関に
よる時間的平滑化処理が行われる。すなわち、折返し現
象を生じた場合、従来装置では第13図において実線2
で示すように周波数ゼロラインと交差する方向に平滑化
されるのに対し、本実施例装置では破線4で示すように
周波数ぜロラインと交差しない方向に平滑化されるので
おる。
これによれば、ゼロ交差部分3は存在せず、従ってCF
M像をフリーズした場合でも、高速血流部位が黒色表示
されるといった不都合は生じない。
尚、本発明は上記実施例に限定されない。
上記実施例では一定条件下で初期条件ROM30より初
期条件(+fr/2及び−fr/2)を与えるようにし
たものについて説明したが、係数a=1としてリセット
する方式と等価な方法を採用することもできる。この場
合、画像メモリ28の出力データについてそのままの状
態で初期条件ROM30を通過させ、それをMPX27
で選択させればよい。
また上記実施例ではCFMユニット15内に平滑止処理
手段22を設けたものについて説明したが、フレーム相
関処理は第1図のカラーフレームメモリ16を利用して
行うこともできるから、このCFMユニット15の後段
に平滑化処理手段を配置することもできる。
[発明の効果] 以上詳述したように本発明によれば、血流分布(CFM
)像のフレーム相関による時間的平滑化処理において折
返し現象を生じた場合の高血流速部位の黒色表示を排除
することができ、適切な血流分布象表示を可能とするカ
ラー超音波診断装置を提供できる。
【図面の簡単な説明】
第1図は本発明の一実施例装置のブロック図第2図は同
上装置の主要部の詳細な構成ブロック図、第3図、第4
図、第6図、第7図は本実施例装置の動作タイミング図
、第5図は同上装置における画像メモリのフォーマット
図、第8図は同上装置におけるMPXのデータ選択論理
の説明図、第9図はフレーム相関処理におけるデータサ
ンプリング説明図、第10図はフレーム相関処理を実行
する巡回型フィルタのブロック図、第11図はフレーム
相関処理の特性図、第12図は符号無しデータのフレー
ム相関による時間的平滑化処理の説明図、爪13図はC
FM像データのフレーム相関による時間的平滑化処理の
説明図、第14図は超音波スキャンのパターン図、第1
5図は平均血流速の角度表現1周波数表現とそれに対応
する表示色との関係説明図である。 10・・・超音波プローブ、 15・・・CFMユニット、 18・・・RGBモニタ、 21・・・CFM演算部、 22・・・平滑化処理手段。 OFO 第 図 RF 第 図 第 図 第 ■ 図 (1=0.75 第 ■ ■ 図

Claims (1)

    【特許請求の範囲】
  1. 被検体よりの超音波エコーに基づいて形成された血流分
    布像に対してフレーム相関による時間的平滑化処理を施
    し、その処理結果をカラー表示するようにしたカラー超
    音波診断装置において、血流分布像のフレーム周波数を
    f_D_Fとし、隣接するフレーム間での血流分布像デ
    ータの変化をfとしたとき、f/f_D_F≧1/2の
    場合に血流情報の周波数表現上の周波数ゼロラインと交
    差しない方向に血流分布像のフレーム相関による時間的
    平滑化処理を行う平滑化処理手段を具備することを特徴
    とするカラー超音波診断装置。
JP30754688A 1988-12-05 1988-12-05 カラー超音波診断装置 Expired - Lifetime JPH0716489B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP30754688A JPH0716489B2 (ja) 1988-12-05 1988-12-05 カラー超音波診断装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP30754688A JPH0716489B2 (ja) 1988-12-05 1988-12-05 カラー超音波診断装置

Publications (2)

Publication Number Publication Date
JPH02152447A true JPH02152447A (ja) 1990-06-12
JPH0716489B2 JPH0716489B2 (ja) 1995-03-01

Family

ID=17970395

Family Applications (1)

Application Number Title Priority Date Filing Date
JP30754688A Expired - Lifetime JPH0716489B2 (ja) 1988-12-05 1988-12-05 カラー超音波診断装置

Country Status (1)

Country Link
JP (1) JPH0716489B2 (ja)

Also Published As

Publication number Publication date
JPH0716489B2 (ja) 1995-03-01

Similar Documents

Publication Publication Date Title
US5409007A (en) Filter to reduce speckle artifact in ultrasound imaging
US5429137A (en) Acoustic scan conversion method and apparatus for velocity flow
US6406430B1 (en) Ultrasound image display by combining enhanced flow imaging in B-mode and color flow mode
US6814703B2 (en) Apparatus and method for ultrasonic diagnostic imaging using a contrast medium
US5349525A (en) Color flow imaging system utilizing a frequency domain wall filter
US5301670A (en) Ultrasonic diagnosis apparatus
JPH062134B2 (ja) 超音波診断装置
JP3251696B2 (ja) 超音波診断装置
US4930514A (en) Ultrasonic imaging apparatus
US6059729A (en) Method and apparatus for edge enhancement in ultrasound imaging
US5865752A (en) Method and apparatus for ultrasound imaging using normalized difference between successive frames
JP4717978B2 (ja) セグメント化およびフロー・ダイナミクスを強化するためのカラー・フロー撮像
JP2763126B2 (ja) カラー超音波診断装置
JPH078492A (ja) 超音波診断装置
US6048313A (en) Method and apparatus for fractal-based enhancement of ultrasound imaging
JP2719707B2 (ja) カラードプラmti装置の画像表示方法
JPH02152447A (ja) カラー超音波診断装置
JP2639391B2 (ja) 流体マップの画像形成装置
JPH1189846A (ja) 超音波診断装置
JP2714329B2 (ja) 超音波診断装置
JP2938125B2 (ja) 超音波診断装置
JP2950600B2 (ja) カラードプラ超音波診断装置
JP2602898B2 (ja) 超音波血流イメージング装置
JP3332090B2 (ja) 超音波診断装置
JP2963193B2 (ja) 超音波診断装置

Legal Events

Date Code Title Description
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080301

Year of fee payment: 13

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090301

Year of fee payment: 14

EXPY Cancellation because of completion of term
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090301

Year of fee payment: 14