JP3408647B2 - 超音波イメージング装置 - Google Patents
超音波イメージング装置Info
- Publication number
- JP3408647B2 JP3408647B2 JP29976694A JP29976694A JP3408647B2 JP 3408647 B2 JP3408647 B2 JP 3408647B2 JP 29976694 A JP29976694 A JP 29976694A JP 29976694 A JP29976694 A JP 29976694A JP 3408647 B2 JP3408647 B2 JP 3408647B2
- Authority
- JP
- Japan
- Prior art keywords
- data
- velocity
- signal
- moving object
- processing unit
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
- G01N29/0609—Display arrangements, e.g. colour displays
Landscapes
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Indicating Or Recording The Presence, Absence, Or Direction Of Movement (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Description
【０００１】
【産業上の利用分野】本発明は超音波イメージング表示
方法及び超音波イメージング装置に関し、特に、血流に
関して拍動流か定常流かを識別可能な状態で表示するこ
とが可能な超音波イメージング表示方法及び超音波イメ
ージング装置に関する。
方法及び超音波イメージング装置に関し、特に、血流に
関して拍動流か定常流かを識別可能な状態で表示するこ
とが可能な超音波イメージング表示方法及び超音波イメ
ージング装置に関する。
【０００２】
【従来の技術】超音波を被検体に照射すると、生体組織
を媒体として超音波が伝達されるが、臓器等の組織や病
変部のような周囲の組織との音響インピーダンスの差の
ある所から反射されて、また、血液における赤血球など
からは散乱によって、送波した超音波の一部が戻ってく
る。この反射体や散乱体が視線方向に運動または移動す
る物体であった場合、その反射波の周波数はドプラ効果
によって送信周波数から偏移する。超音波診断装置はこ
の周波数偏移量を測定して運動物体の速度及び移動方向
を表示観察して診断の用に供する装置である。
を媒体として超音波が伝達されるが、臓器等の組織や病
変部のような周囲の組織との音響インピーダンスの差の
ある所から反射されて、また、血液における赤血球など
からは散乱によって、送波した超音波の一部が戻ってく
る。この反射体や散乱体が視線方向に運動または移動す
る物体であった場合、その反射波の周波数はドプラ効果
によって送信周波数から偏移する。超音波診断装置はこ
の周波数偏移量を測定して運動物体の速度及び移動方向
を表示観察して診断の用に供する装置である。
【０００３】反射体が探触子の方に向かって動いている
時は受信周波数は送信周波数より高く、逆に遠ざかる時
は低くなる。そしてその偏移周波数は反射体の運動速度
に比例する。このドプラ効果を利用して、例えば、心臓
や血管内を流れる血流の方向と速度を知ることができ
る。
時は受信周波数は送信周波数より高く、逆に遠ざかる時
は低くなる。そしてその偏移周波数は反射体の運動速度
に比例する。このドプラ効果を利用して、例えば、心臓
や血管内を流れる血流の方向と速度を知ることができ
る。
【０００４】この超音波診断装置において、運動物体の
移動方向を例えば近付く場合に赤，遠ざかる場合に青の
ように色で表示し、その移動速度を輝度で表示するＣＦ
Ｍ（Color Flow Mapping）モードを有する装置がある。
移動方向を例えば近付く場合に赤，遠ざかる場合に青の
ように色で表示し、その移動速度を輝度で表示するＣＦ
Ｍ（Color Flow Mapping）モードを有する装置がある。
【０００５】ところで、超音波パルスの反射には固定物
体によるものと運動物体によるものとが混在しており、
特に運動物体に注目する時は固定物体による反射信号を
消去して運動物体の反射信号のみを表示したほうが分り
易いので、その目的のために、運動物体の反射信号のみ
を取り出すＭＴＩ(Moving Target Indicator) が用いら
れている。
体によるものと運動物体によるものとが混在しており、
特に運動物体に注目する時は固定物体による反射信号を
消去して運動物体の反射信号のみを表示したほうが分り
易いので、その目的のために、運動物体の反射信号のみ
を取り出すＭＴＩ(Moving Target Indicator) が用いら
れている。
【０００６】しかし、固定した反射物体に関する情報も
医用超音波診断装置では重要な意味も持っている。それ
は、上記の情報が無いと表示画面中のどこが何を表わし
ているのか特定できないし、また、単に運動物体を表示
しただけではそれが何であるかを判別することが困難で
ある。
医用超音波診断装置では重要な意味も持っている。それ
は、上記の情報が無いと表示画面中のどこが何を表わし
ているのか特定できないし、また、単に運動物体を表示
しただけではそれが何であるかを判別することが困難で
ある。
【０００７】このため、固定物体の表示と運動物体の表
示とを同一画面上に共通に表示し、相互に識別するため
に、血球からの散乱エコーから血流の方向，平均速度と
速度分散値及びパワー値などを演算し、色でそれらの諸
量をリアルタイムでＢモード画面上に２次元表示してい
る。このモードをＣＦＭモードといい、これらの機能を
有する装置がＣＦＭ超音波診断装置である。
示とを同一画面上に共通に表示し、相互に識別するため
に、血球からの散乱エコーから血流の方向，平均速度と
速度分散値及びパワー値などを演算し、色でそれらの諸
量をリアルタイムでＢモード画面上に２次元表示してい
る。このモードをＣＦＭモードといい、これらの機能を
有する装置がＣＦＭ超音波診断装置である。
【０００８】図６に従来のＣＦＭモードを有する超音波
診断装置の構成図を示す。この図６の構成において、超
音波探触子１は送波電気信号を超音波に変換して被検体
内に送波し、被検体内から反射されて戻って来た超音波
信号を電気信号に変換する。送受信部２は送信信号を増
幅して超音波探触子１に送り、受波された受信信号を直
交検波等により復調した後に各処理部に送る。すなわ
ち、送受信部２の出力（復調出力（ｉ，ｑ信号等））は
Ｂモード処理部３及びＣＦＭ演算部４にそれぞれ供給さ
れる。
診断装置の構成図を示す。この図６の構成において、超
音波探触子１は送波電気信号を超音波に変換して被検体
内に送波し、被検体内から反射されて戻って来た超音波
信号を電気信号に変換する。送受信部２は送信信号を増
幅して超音波探触子１に送り、受波された受信信号を直
交検波等により復調した後に各処理部に送る。すなわ
ち、送受信部２の出力（復調出力（ｉ，ｑ信号等））は
Ｂモード処理部３及びＣＦＭ演算部４にそれぞれ供給さ
れる。
【０００９】Ｂモード処理部３では増幅，対数圧縮，検
波等を行ってＢモード表示のための信号を出力して、画
像処理部５に供給する。また、送受信部２の復調出力は
ＣＦＭ演算部４に入力される。このＣＦＭ演算部４で
は、ＭＴＩフィルタ等により移動目標のみの反射波を抽
出し、例えば血流からの反射波を抽出して演算を行い血
流速度を検出する。また、血流速度のばらつきにより、
分散（速度分散）を検出する。
波等を行ってＢモード表示のための信号を出力して、画
像処理部５に供給する。また、送受信部２の復調出力は
ＣＦＭ演算部４に入力される。このＣＦＭ演算部４で
は、ＭＴＩフィルタ等により移動目標のみの反射波を抽
出し、例えば血流からの反射波を抽出して演算を行い血
流速度を検出する。また、血流速度のばらつきにより、
分散（速度分散）を検出する。
【００１０】これらＢモード処理部３及びＣＦＭ演算部
４のそれぞれの出力は画像処理部５に供給され、Ｂモー
ド表示上に運動物体（血流速度，速度分散）に応じて着
色されたイメージを生成し、このイメージの表示をＣＲ
Ｔ表示部７に行う。
４のそれぞれの出力は画像処理部５に供給され、Ｂモー
ド表示上に運動物体（血流速度，速度分散）に応じて着
色されたイメージを生成し、このイメージの表示をＣＲ
Ｔ表示部７に行う。
【００１１】例えば、移動方向で色（色相）を変え（赤
／青）、移動速度に応じて赤／青の輝度を変えるように
して表示を行う。また、速度分散に応じて緑を混ぜるよ
うにした表示を行う。通常、赤は近づく流れを、青は遠
ざかる流れを表し、分散の大きさは緑の輝度の違いよっ
て表している。従って、部位によって、青と緑との混合
が行われ、また、赤と緑との混合が行われる。この場
合、赤，青，緑は光の三原色であるので、混合された色
の変化により、方向，速度，分散の大きさを読みとれ
る。
／青）、移動速度に応じて赤／青の輝度を変えるように
して表示を行う。また、速度分散に応じて緑を混ぜるよ
うにした表示を行う。通常、赤は近づく流れを、青は遠
ざかる流れを表し、分散の大きさは緑の輝度の違いよっ
て表している。従って、部位によって、青と緑との混合
が行われ、また、赤と緑との混合が行われる。この場
合、赤，青，緑は光の三原色であるので、混合された色
の変化により、方向，速度，分散の大きさを読みとれ
る。
【００１２】図７にＣＦＭ演算部４の内部の回路構成の
一例を示す。この図７において、前述の図６の送受信部
２の復調出力が位相が直交する信号（ｉ信号，ｑ信号）
であった場合、ｉ，ｑ信号それぞれがＭＴＩフィルタ４
００，４０１によりフィルタリングされて、運動物体か
らの信号成分のみの信号が抽出される。
一例を示す。この図７において、前述の図６の送受信部
２の復調出力が位相が直交する信号（ｉ信号，ｑ信号）
であった場合、ｉ，ｑ信号それぞれがＭＴＩフィルタ４
００，４０１によりフィルタリングされて、運動物体か
らの信号成分のみの信号が抽出される。
【００１３】ＭＴＩフィルタ４００，４０１を通過後の
運動物体からの信号をそれぞれＩ，Ｑ信号とする。ここ
で、Ｉ，Ｑ信号それぞれは遅延部４０２，４０３に供給
される。この遅延部４０２，４０３の遅延時間はΔｔで
あるとする。また、この遅延時間Δｔは、送受信部２で
同一音線の信号を受信する間隔（若しくはその整数倍）
の間隔とする。そして、遅延部４０２，４０３により遅
延された信号と遅延されていない信号とがそれぞれ掛算
器４０５〜４０８で掛算され、更に加算器４０９と減算
器４１０で加減算される。
運動物体からの信号をそれぞれＩ，Ｑ信号とする。ここ
で、Ｉ，Ｑ信号それぞれは遅延部４０２，４０３に供給
される。この遅延部４０２，４０３の遅延時間はΔｔで
あるとする。また、この遅延時間Δｔは、送受信部２で
同一音線の信号を受信する間隔（若しくはその整数倍）
の間隔とする。そして、遅延部４０２，４０３により遅
延された信号と遅延されていない信号とがそれぞれ掛算
器４０５〜４０８で掛算され、更に加算器４０９と減算
器４１０で加減算される。
【００１４】ここで、Ｉ信号として、Ｉ1 ，Ｉ2 ，Ｉ
3 ，Ｉ4 とΔｔ毎の４信号を考える。同様にして、Ｑ信
号として、Ｑ1 ，Ｑ2 ，Ｑ3 ，Ｑ4 のΔｔ毎の４信号を
考える。このような信号が与えられた場合、加算器４０
９の出力にはΔｔ毎に、Ｉ1Ｉ 2 ＋Ｑ1 Ｑ2 ，Ｉ2 Ｉ3
＋Ｑ2 Ｑ3 ，Ｉ3 Ｉ4 ＋Ｑ3 Ｑ4 が得られる。この信
号をＤＥＮと呼ぶことにする。
3 ，Ｉ4 とΔｔ毎の４信号を考える。同様にして、Ｑ信
号として、Ｑ1 ，Ｑ2 ，Ｑ3 ，Ｑ4 のΔｔ毎の４信号を
考える。このような信号が与えられた場合、加算器４０
９の出力にはΔｔ毎に、Ｉ1Ｉ 2 ＋Ｑ1 Ｑ2 ，Ｉ2 Ｉ3
＋Ｑ2 Ｑ3 ，Ｉ3 Ｉ4 ＋Ｑ3 Ｑ4 が得られる。この信
号をＤＥＮと呼ぶことにする。
【００１５】また、同様に減算器４１０の出力にはΔｔ
毎に、Ｉ1 Ｑ2 −Ｑ1 Ｉ2 ，Ｉ2 Ｑ3 −Ｑ2 Ｉ3 ，Ｉ3
Ｑ4 −Ｑ3 Ｉ4 が得られる。この信号をＮＵＭと呼ぶこ
とにする。
毎に、Ｉ1 Ｑ2 −Ｑ1 Ｉ2 ，Ｉ2 Ｑ3 −Ｑ2 Ｉ3 ，Ｉ3
Ｑ4 −Ｑ3 Ｉ4 が得られる。この信号をＮＵＭと呼ぶこ
とにする。
【００１６】また、パワー算出部４０４では、２乗加算
処理が行なわれて、Δｔ毎に、Ｉ12 ＋Ｑ1 2 ，Ｉ2 2 ＋
Ｑ2 2 ，Ｉ3 2 ＋Ｑ3 2 ，Ｉ4 2 ＋Ｑ4 2 が得られる。
この信号をパワー値と呼び、Ｐと略することにする。
処理が行なわれて、Δｔ毎に、Ｉ12 ＋Ｑ1 2 ，Ｉ2 2 ＋
Ｑ2 2 ，Ｉ3 2 ＋Ｑ3 2 ，Ｉ4 2 ＋Ｑ4 2 が得られる。
この信号をパワー値と呼び、Ｐと略することにする。
【００１７】そして、遅延加算器４１１にて、Δｔ毎の
遅延加算（累加平均）が行なわれて、（（Ｉ1 Ｉ2 ＋Ｑ
1 Ｑ2 ）＋（Ｉ2 Ｉ3 ＋Ｑ2 Ｑ3 ）＋（Ｉ3 Ｉ4 ＋Ｑ3
Ｑ4））／３が得られる。これをここでは、＊ＤＥＮ＊
として略する。
遅延加算（累加平均）が行なわれて、（（Ｉ1 Ｉ2 ＋Ｑ
1 Ｑ2 ）＋（Ｉ2 Ｉ3 ＋Ｑ2 Ｑ3 ）＋（Ｉ3 Ｉ4 ＋Ｑ3
Ｑ4））／３が得られる。これをここでは、＊ＤＥＮ＊
として略する。
【００１８】同様にして、遅延加算器４１２にて、Δｔ
毎の遅延加算（累加平均）が行なわれて、（（Ｉ1 Ｑ2
−Ｑ1 Ｉ2 ）＋（Ｉ2 Ｑ3 −Ｑ2 Ｉ3 ）＋（Ｉ3 Ｑ4 −
Ｑ3Ｉ4 ））／３が得られる。同様にこれを＊ＮＵＭ＊
と略する。
毎の遅延加算（累加平均）が行なわれて、（（Ｉ1 Ｑ2
−Ｑ1 Ｉ2 ）＋（Ｉ2 Ｑ3 −Ｑ2 Ｉ3 ）＋（Ｉ3 Ｑ4 −
Ｑ3Ｉ4 ））／３が得られる。同様にこれを＊ＮＵＭ＊
と略する。
【００１９】また、遅延加算器４１３にて、Δｔ毎の遅
延加算（累加平均）が行なわれて、（（Ｉ1 2 ＋Ｑ1
2 ）＋（Ｉ2 2 ＋Ｑ2 2 ）＋（Ｉ3 2 ＋Ｑ3 2 ）＋
（Ｉ4 2 ＋Ｑ4 2 ））／４が得られる。これを＊Ｐ＊と
略する。
延加算（累加平均）が行なわれて、（（Ｉ1 2 ＋Ｑ1
2 ）＋（Ｉ2 2 ＋Ｑ2 2 ）＋（Ｉ3 2 ＋Ｑ3 2 ）＋
（Ｉ4 2 ＋Ｑ4 2 ））／４が得られる。これを＊Ｐ＊と
略する。
【００２０】このような＊ＤＥＮ＊，＊ＮＵＭ＊，＊Ｐ
＊が得られると、平均速度算出部４１４において、ｔａ
ｎ-1（（＊ＮＵＭ＊）／（＊ＤＥＮ＊））という処理が
実行されて、平均速度＊ｖ＊が求められる。
＊が得られると、平均速度算出部４１４において、ｔａ
ｎ-1（（＊ＮＵＭ＊）／（＊ＤＥＮ＊））という処理が
実行されて、平均速度＊ｖ＊が求められる。
【００２１】また、分散算出部４１５において、（１−
（（＊ＤＥＮ＊2 ＋＊ＮＵＭ＊2 ） 1/2）／＊Ｐ＊）と
いう処理が実行されて、分散＊ｖａｒ＊が求められる。
このようにして得られた平均速度と分散とを受けた画像
処理部５が、それぞれの信号値に応じて赤／青，緑など
の輝度を変えてカラー表示を行う。
（（＊ＤＥＮ＊2 ＋＊ＮＵＭ＊2 ） 1/2）／＊Ｐ＊）と
いう処理が実行されて、分散＊ｖａｒ＊が求められる。
このようにして得られた平均速度と分散とを受けた画像
処理部５が、それぞれの信号値に応じて赤／青，緑など
の輝度を変えてカラー表示を行う。
【００２２】以上のような平均速度，分散と血流との関
係を図８に示す。図８（イ）は血流速度を示しており、
横軸は時間ｔを示し、縦軸は血流速度ｖを示している。
尚、ハッチングで示した部分はバラツキを意味してい
る。ここでは、１心拍に相当する時間を示しており、こ
こでは時相Ｉ，ＩＩ，ＩＩＩ，ＩＶ，Ｖのそれぞれのタ
イミングで上述の平均速度と分散とを求める処理を実行
している。
係を図８に示す。図８（イ）は血流速度を示しており、
横軸は時間ｔを示し、縦軸は血流速度ｖを示している。
尚、ハッチングで示した部分はバラツキを意味してい
る。ここでは、１心拍に相当する時間を示しており、こ
こでは時相Ｉ，ＩＩ，ＩＩＩ，ＩＶ，Ｖのそれぞれのタ
イミングで上述の平均速度と分散とを求める処理を実行
している。
【００２３】従って、各時相Ｉ〜Ｖに応じて平均速度と
分散とが求められ、それぞれ図８（ロ）〜（ヘ）のよう
な速度スペクトラムになる。この各図では、山のピーク
の位置が平均速度、山の広がりが分散を示し、縦軸は強
度を示すスペクトラムを示している。
分散とが求められ、それぞれ図８（ロ）〜（ヘ）のよう
な速度スペクトラムになる。この各図では、山のピーク
の位置が平均速度、山の広がりが分散を示し、縦軸は強
度を示すスペクトラムを示している。
【００２４】尚、ここでは、各時相の平均速度と分散と
をグラフで示したが、実際には画像処理部５の処理によ
りカラー表示が行なわれる。従って、実際の装置では、
図８（ロ）〜（ヘ）に応じて、平均速度と分散とのカラ
ー表示が刻々と変換するように示される。
をグラフで示したが、実際には画像処理部５の処理によ
りカラー表示が行なわれる。従って、実際の装置では、
図８（ロ）〜（ヘ）に応じて、平均速度と分散とのカラ
ー表示が刻々と変換するように示される。
【００２５】
【発明が解決しようとする課題】以上のようにして平均
速度と分散とを表示する超音波診断装置は、異なった時
相毎に血流の平均速度と分散とを表示するものである。
速度と分散とを表示する超音波診断装置は、異なった時
相毎に血流の平均速度と分散とを表示するものである。
【００２６】従って、時間的な拍動性を表示することは
できず、血流が拍動流（血流速度や分散が周期的に変化
する流れ）なのか定常流（血流速度や分散がほぼ一定な
ような流れ）なのかを瞬時に識別することは困難であっ
た。
できず、血流が拍動流（血流速度や分散が周期的に変化
する流れ）なのか定常流（血流速度や分散がほぼ一定な
ような流れ）なのかを瞬時に識別することは困難であっ
た。
【００２７】このような場合に、拍動性を識別可能なよ
うに表示するためには、パルスドプラなどによりＦＦＴ
で周波数分析を行なって検証する必要が有り、装置が複
雑化する欠点がある。
うに表示するためには、パルスドプラなどによりＦＦＴ
で周波数分析を行なって検証する必要が有り、装置が複
雑化する欠点がある。
【００２８】また、ＣＦＭ演算部４の出力を長時間トレ
ースするような回路を設けることによっても識別が可能
であるが、大規模な回路付加が必要になる問題を有して
いる。
ースするような回路を設けることによっても識別が可能
であるが、大規模な回路付加が必要になる問題を有して
いる。
【００２９】以上のような問題のために、簡単な構成で
拍動性を識別可能なように表示することは極めて困難で
あった。本発明は上記の点に鑑みてなされたもので、第
１の目的は、血流に関して拍動流か定常流かを識別可能
な状態で表示することが可能な超音波イメージング表示
方法を実現することである。
拍動性を識別可能なように表示することは極めて困難で
あった。本発明は上記の点に鑑みてなされたもので、第
１の目的は、血流に関して拍動流か定常流かを識別可能
な状態で表示することが可能な超音波イメージング表示
方法を実現することである。
【００３０】また、第２の目的は、血流に関して拍動流
か定常流かを識別可能な状態で表示することが可能な超
音波イメージング装置を実現することである。
か定常流かを識別可能な状態で表示することが可能な超
音波イメージング装置を実現することである。
【００３１】
【課題を解決するための手段】前記の課題を解決する第
１の手段は、受波信号より得られた運動物体反射信号中
の同一音線データのサンプリング間隔Δｔより大きい間
隔ΔＴ毎に所定個数の各運動物体反射信号の遅延加算を
行い、遅延加算された運動物体反射信号から流速データ
を求めてカラーフローマッピングによる超音波イメージ
ング表示を行うことを特徴とする超音波イメージング表
示方法である。
１の手段は、受波信号より得られた運動物体反射信号中
の同一音線データのサンプリング間隔Δｔより大きい間
隔ΔＴ毎に所定個数の各運動物体反射信号の遅延加算を
行い、遅延加算された運動物体反射信号から流速データ
を求めてカラーフローマッピングによる超音波イメージ
ング表示を行うことを特徴とする超音波イメージング表
示方法である。
【００３２】ここで、運動物体反射信号とは受波信号か
ら検出された運動物体による反射波成分の信号であり、
同一音線のデータはΔｔ間隔で得られる。また、Δｔよ
り大きい間隔ΔＴとは異なる時相となるような時間間隔
を意味する。また、流速データとは、血流等の速度等を
示すデータである。
ら検出された運動物体による反射波成分の信号であり、
同一音線のデータはΔｔ間隔で得られる。また、Δｔよ
り大きい間隔ΔＴとは異なる時相となるような時間間隔
を意味する。また、流速データとは、血流等の速度等を
示すデータである。
【００３３】前記の課題を解決する第２の手段は、受波
信号を分析してＢモードデータを得るＢモード処理部
と、受波信号を分析して運動物体反射信号を求め、同一
音線データのサンプリング間隔Δｔより大きい間隔ΔＴ
毎に所定個数の運動物体反射信号の遅延加算を行って流
速データを求めるＣＦＭ演算部と、前記Ｂモード処理部
で得られたＢモードデータと前記ＣＦＭ演算部で遅延加
算されて得られた流速データとを用いてカラーフローマ
ッピングによる超音波イメージング表示用のイメージデ
ータを生成する画像処理部とを備えたことを特徴とする
超音波イメージング装置である。
信号を分析してＢモードデータを得るＢモード処理部
と、受波信号を分析して運動物体反射信号を求め、同一
音線データのサンプリング間隔Δｔより大きい間隔ΔＴ
毎に所定個数の運動物体反射信号の遅延加算を行って流
速データを求めるＣＦＭ演算部と、前記Ｂモード処理部
で得られたＢモードデータと前記ＣＦＭ演算部で遅延加
算されて得られた流速データとを用いてカラーフローマ
ッピングによる超音波イメージング表示用のイメージデ
ータを生成する画像処理部とを備えたことを特徴とする
超音波イメージング装置である。
【００３４】ここで、流速データとは受波信号から得ら
れる血流等の速度や分散等を示すデータのことである。
前記の課題を解決する第３の手段は、受波信号を分析し
てＢモードデータを得るＢモード処理部と、受波信号を
分析して運動物体反射信号を求め、同一音線データのサ
ンプリング間隔Δｔより大きい間隔ΔＴ毎に所定個数の
運動物体反射信号の遅延加算を行って速度データ及び分
散データを求めるＣＦＭ演算部と、前記Ｂモード処理部
で得られたＢモードデータ，前記ＣＦＭ演算部で遅延加
算されて得られた速度データ及び分散データを用いてカ
ラーフローマッピングによる超音波イメージング表示用
のイメージデータを生成する画像処理部とを備えたこと
を特徴とする超音波イメージング装置である。
れる血流等の速度や分散等を示すデータのことである。
前記の課題を解決する第３の手段は、受波信号を分析し
てＢモードデータを得るＢモード処理部と、受波信号を
分析して運動物体反射信号を求め、同一音線データのサ
ンプリング間隔Δｔより大きい間隔ΔＴ毎に所定個数の
運動物体反射信号の遅延加算を行って速度データ及び分
散データを求めるＣＦＭ演算部と、前記Ｂモード処理部
で得られたＢモードデータ，前記ＣＦＭ演算部で遅延加
算されて得られた速度データ及び分散データを用いてカ
ラーフローマッピングによる超音波イメージング表示用
のイメージデータを生成する画像処理部とを備えたこと
を特徴とする超音波イメージング装置である。
【００３５】ここで、速度データは血流等の速度を示
し、分散データは血流等の速度のばらつきを示してい
る。
し、分散データは血流等の速度のばらつきを示してい
る。
【００３６】
【作用】課題を解決する第１の手段である超音波イメー
ジング表示方法において、Δｔ間隔で得られた運動物体
反射信号から流速データが求められるが、この際に、Δ
ｔより大きい間隔であって異なる時相のΔＴ毎の運動物
体反射信号が遅延加算されたものにより流速データが求
められる。そして、このようにして求められた流速デー
タによりＣＦＭの超音波イメージング表示がなされる。
このようにすることで、拍動流か定常流かを識別可能な
状態で表示することが可能になる。
ジング表示方法において、Δｔ間隔で得られた運動物体
反射信号から流速データが求められるが、この際に、Δ
ｔより大きい間隔であって異なる時相のΔＴ毎の運動物
体反射信号が遅延加算されたものにより流速データが求
められる。そして、このようにして求められた流速デー
タによりＣＦＭの超音波イメージング表示がなされる。
このようにすることで、拍動流か定常流かを識別可能な
状態で表示することが可能になる。
【００３７】課題を解決する第２の手段である超音波イ
メージング装置において、Δｔ間隔で得られた運動物体
反射信号から流速データが求められるが、この際に、Δ
ｔより大きい間隔であって異なる時相のΔＴ毎の運動物
体反射信号が遅延加算されたものにより流速データが求
められる。そして、このようにして求められた流速デー
タによりＣＦＭの超音波イメージング表示がなされる。
このようにすることで、拍動流か定常流かを識別可能な
状態で表示することが可能になる。
メージング装置において、Δｔ間隔で得られた運動物体
反射信号から流速データが求められるが、この際に、Δ
ｔより大きい間隔であって異なる時相のΔＴ毎の運動物
体反射信号が遅延加算されたものにより流速データが求
められる。そして、このようにして求められた流速デー
タによりＣＦＭの超音波イメージング表示がなされる。
このようにすることで、拍動流か定常流かを識別可能な
状態で表示することが可能になる。
【００３８】課題を解決する第３の手段である超音波イ
メージング装置において、Δｔ間隔で得られた運動物体
反射信号から速度データ及び分散データが求められる
が、この際に、Δｔより大きい間隔であって異なる時相
のΔＴ毎の運動物体反射信号が遅延加算されたものによ
り速度データ及び分散データが求められる。そして、こ
のようにして求められた速度データ及び分散データによ
りＣＦＭの超音波イメージング表示がなされる。このよ
うにすることで、拍動流か定常流かを識別可能な状態で
表示することが可能になる。
メージング装置において、Δｔ間隔で得られた運動物体
反射信号から速度データ及び分散データが求められる
が、この際に、Δｔより大きい間隔であって異なる時相
のΔＴ毎の運動物体反射信号が遅延加算されたものによ
り速度データ及び分散データが求められる。そして、こ
のようにして求められた速度データ及び分散データによ
りＣＦＭの超音波イメージング表示がなされる。このよ
うにすることで、拍動流か定常流かを識別可能な状態で
表示することが可能になる。
【００３９】
【実施例】以下、図面を参照して本発明の実施例を詳細
に説明する。図１は本発明の超音波イメージング表示方
法の一実施例を実現するための装置であり、また、本発
明の超音波イメージング装置の一実施例の概略構成を示
す構成図である。まず、この図１により超音波イメージ
ング装置の概要を説明する。
に説明する。図１は本発明の超音波イメージング表示方
法の一実施例を実現するための装置であり、また、本発
明の超音波イメージング装置の一実施例の概略構成を示
す構成図である。まず、この図１により超音波イメージ
ング装置の概要を説明する。
【００４０】この図１は前述のＣＦＭ演算部４の内部の
回路構成の一例を示す構成図である。この図１におい
て、前述の図６の送受信部２の復調出力が位相が直交す
る信号（Ｉ信号，Ｑ信号）であった場合、Ｉ，Ｑ信号そ
れぞれがＭＴＩフィルタ４００，４０１によりフィルタ
リングされて、運動物体からの信号成分のみの信号（運
動物体反射信号）が抽出される。
回路構成の一例を示す構成図である。この図１におい
て、前述の図６の送受信部２の復調出力が位相が直交す
る信号（Ｉ信号，Ｑ信号）であった場合、Ｉ，Ｑ信号そ
れぞれがＭＴＩフィルタ４００，４０１によりフィルタ
リングされて、運動物体からの信号成分のみの信号（運
動物体反射信号）が抽出される。
【００４１】ＭＴＩフィルタ４００，４０１を通過後の
運動物体反射信号をそれぞれＩ，Ｑ信号とする。ここ
で、Ｉ，Ｑ信号それぞれは遅延部４０２，４０３に供給
される。この遅延部４０２，４０３の遅延時間はΔｔで
あるとする。また、この遅延時間Δｔは、送受信部２で
同一音線の信号を受信する間隔（若しくはその整数倍）
の間隔とする。そして、遅延部４０２，４０３により遅
延された信号と遅延されていない信号とがそれぞれ掛算
器４０５〜４０８で掛算され、更に加算器４０９と減算
器４１０で加減算される。
運動物体反射信号をそれぞれＩ，Ｑ信号とする。ここ
で、Ｉ，Ｑ信号それぞれは遅延部４０２，４０３に供給
される。この遅延部４０２，４０３の遅延時間はΔｔで
あるとする。また、この遅延時間Δｔは、送受信部２で
同一音線の信号を受信する間隔（若しくはその整数倍）
の間隔とする。そして、遅延部４０２，４０３により遅
延された信号と遅延されていない信号とがそれぞれ掛算
器４０５〜４０８で掛算され、更に加算器４０９と減算
器４１０で加減算される。
【００４２】ここで、Ｉ信号として、Ｉ1 ，Ｉ2 ，Ｉ
3 ，Ｉ4 とΔｔ毎の４信号を考える。同様にして、Ｑ信
号として、Ｑ1 ，Ｑ2 ，Ｑ3 ，Ｑ4 のΔｔ毎の４信号を
考える。このような信号が与えられた場合、加算器４０
９の出力にはΔｔ毎に、Ｉ1Ｉ 2 ＋Ｑ1 Ｑ2 ，Ｉ2 Ｉ3
＋Ｑ2 Ｑ3 ，Ｉ3 Ｉ4 ＋Ｑ3 Ｑ4 が得られる。この信
号をＤＥＮと呼ぶことにする。
3 ，Ｉ4 とΔｔ毎の４信号を考える。同様にして、Ｑ信
号として、Ｑ1 ，Ｑ2 ，Ｑ3 ，Ｑ4 のΔｔ毎の４信号を
考える。このような信号が与えられた場合、加算器４０
９の出力にはΔｔ毎に、Ｉ1Ｉ 2 ＋Ｑ1 Ｑ2 ，Ｉ2 Ｉ3
＋Ｑ2 Ｑ3 ，Ｉ3 Ｉ4 ＋Ｑ3 Ｑ4 が得られる。この信
号をＤＥＮと呼ぶことにする。
【００４３】また、同様に減算器４１０の出力にはΔｔ
毎に、Ｉ1 Ｑ2 −Ｑ1 Ｉ2 ，Ｉ2 Ｑ3 −Ｑ2 Ｉ3 ，Ｉ3
Ｑ4 −Ｑ3 Ｉ4 が得られる。この信号をＮＵＭと呼ぶこ
とにする。
毎に、Ｉ1 Ｑ2 −Ｑ1 Ｉ2 ，Ｉ2 Ｑ3 −Ｑ2 Ｉ3 ，Ｉ3
Ｑ4 −Ｑ3 Ｉ4 が得られる。この信号をＮＵＭと呼ぶこ
とにする。
【００４４】また、パワー算出部４０４では、２乗加算
処理が行なわれて、Δｔ毎に、Ｉ1 2 ＋Ｑ1 2 ，Ｉ2 2
＋Ｑ2 2 ，Ｉ3 2 ＋Ｑ3 2 ，Ｉ4 2 ＋Ｑ4 2 が得られ
る。この信号をパワー値と呼び、Ｐと略することにす
る。
処理が行なわれて、Δｔ毎に、Ｉ1 2 ＋Ｑ1 2 ，Ｉ2 2
＋Ｑ2 2 ，Ｉ3 2 ＋Ｑ3 2 ，Ｉ4 2 ＋Ｑ4 2 が得られ
る。この信号をパワー値と呼び、Ｐと略することにす
る。
【００４５】そして、遅延加算器４１１にて、Δｔ毎の
遅延加算（累加平均）が行なわれて、（（Ｉ1 Ｉ2 ＋Ｑ
1 Ｑ2 ）＋（Ｉ2 Ｉ3 ＋Ｑ2 Ｑ3 ）＋（Ｉ3 Ｉ4 ＋Ｑ3
Ｑ4））／３が得られる。同様にこれを＊ＤＥＮ＊と略
する。そして、図８において説明した時相Ｉ〜Ｖに対応
して、各時相の＊ＤＥＮ＊を、＊ＤＥＮ＊Ｉ，＊ＤＥＮ
＊ＩＩ，＊ＤＥＮ＊ＩＩＩ，＊ＤＥＮ＊ＩＶ，＊ＤＥＮ
＊Ｖのように示す。
遅延加算（累加平均）が行なわれて、（（Ｉ1 Ｉ2 ＋Ｑ
1 Ｑ2 ）＋（Ｉ2 Ｉ3 ＋Ｑ2 Ｑ3 ）＋（Ｉ3 Ｉ4 ＋Ｑ3
Ｑ4））／３が得られる。同様にこれを＊ＤＥＮ＊と略
する。そして、図８において説明した時相Ｉ〜Ｖに対応
して、各時相の＊ＤＥＮ＊を、＊ＤＥＮ＊Ｉ，＊ＤＥＮ
＊ＩＩ，＊ＤＥＮ＊ＩＩＩ，＊ＤＥＮ＊ＩＶ，＊ＤＥＮ
＊Ｖのように示す。
【００４６】同様にして、遅延加算器４１２にて、Δｔ
毎の遅延加算（累加平均）が行なわれて、（（Ｉ1 Ｑ2
−Ｑ1 Ｉ2 ）＋（Ｉ2 Ｑ3 −Ｑ2 Ｉ3 ）＋（Ｉ3 Ｑ4 −
Ｑ3 Ｉ4 ））／３が得られる。これをここでは、＊ＮＵ
Ｍ＊として略する。
毎の遅延加算（累加平均）が行なわれて、（（Ｉ1 Ｑ2
−Ｑ1 Ｉ2 ）＋（Ｉ2 Ｑ3 −Ｑ2 Ｉ3 ）＋（Ｉ3 Ｑ4 −
Ｑ3 Ｉ4 ））／３が得られる。これをここでは、＊ＮＵ
Ｍ＊として略する。
【００４７】そして、図８において説明した時相Ｉ〜Ｖ
に対応して、各時相の＊ＮＵＭ＊を、＊ＮＵＭ＊Ｉ，＊
ＮＵＭ＊ＩＩ，＊ＮＵＭ＊ＩＩＩ，＊ＮＵＭ＊ＩＶ，＊
ＮＵＭ＊Ｖのように示す。
に対応して、各時相の＊ＮＵＭ＊を、＊ＮＵＭ＊Ｉ，＊
ＮＵＭ＊ＩＩ，＊ＮＵＭ＊ＩＩＩ，＊ＮＵＭ＊ＩＶ，＊
ＮＵＭ＊Ｖのように示す。
【００４８】また、遅延加算器４１３にて、Δｔ毎の遅
延加算（累加平均）が行なわれて、（（Ｉ1 2 ＋Ｑ1
2 ）＋（Ｉ2 2 ＋Ｑ2 2 ）＋（Ｉ3 2 ＋Ｑ3 2 ）＋
（Ｉ4 2 ＋Ｑ4 2 ））／４が得られる。これを＊Ｐ＊と
略する。そして、図８において説明した時相Ｉ〜Ｖに対
応して、各時相の＊Ｐ＊を、＊Ｐ＊Ｉ，＊Ｐ＊ＩＩ，＊
Ｐ＊ＩＩＩ，＊Ｐ＊ＩＶ，＊Ｐ＊Ｖのように示す。
延加算（累加平均）が行なわれて、（（Ｉ1 2 ＋Ｑ1
2 ）＋（Ｉ2 2 ＋Ｑ2 2 ）＋（Ｉ3 2 ＋Ｑ3 2 ）＋
（Ｉ4 2 ＋Ｑ4 2 ））／４が得られる。これを＊Ｐ＊と
略する。そして、図８において説明した時相Ｉ〜Ｖに対
応して、各時相の＊Ｐ＊を、＊Ｐ＊Ｉ，＊Ｐ＊ＩＩ，＊
Ｐ＊ＩＩＩ，＊Ｐ＊ＩＶ，＊Ｐ＊Ｖのように示す。
【００４９】ここで、遅延加算器４１６にて、Δｔより
大きい間隔ΔＴ（異なる時相となるような間隔：ここで
は、ΔＴを時相Ｉ，ＩＩ，ＩＩＩ，ＩＶ，Ｖの間隔に合
わせた場合を示す）毎の遅延加算（累加平均）が行われ
る。この結果、（（＊ＤＥＮ＊Ｉ）＋（＊ＤＥＮ＊Ｉ
Ｉ）＋（＊ＤＥＮ＊ＩＩＩ）＋（＊ＤＥＮ＊ＩＶ）＋
（＊ＤＥＮ＊Ｖ）／５）のような演算が行われる。この
演算結果を、＊ＤＥＮ＊Ｔと呼ぶ。この＊ＤＥＮ＊Ｔ
は、異なる時相Ｉ〜ＶのデータＤＥＮの平均値である。
大きい間隔ΔＴ（異なる時相となるような間隔：ここで
は、ΔＴを時相Ｉ，ＩＩ，ＩＩＩ，ＩＶ，Ｖの間隔に合
わせた場合を示す）毎の遅延加算（累加平均）が行われ
る。この結果、（（＊ＤＥＮ＊Ｉ）＋（＊ＤＥＮ＊Ｉ
Ｉ）＋（＊ＤＥＮ＊ＩＩＩ）＋（＊ＤＥＮ＊ＩＶ）＋
（＊ＤＥＮ＊Ｖ）／５）のような演算が行われる。この
演算結果を、＊ＤＥＮ＊Ｔと呼ぶ。この＊ＤＥＮ＊Ｔ
は、異なる時相Ｉ〜ＶのデータＤＥＮの平均値である。
【００５０】また同様に、遅延加算器４１７にて、Δｔ
より大きい間隔ΔＴ（異なる時相となるような間隔）毎
の遅延加算（累加平均）が行われる。この結果、（（＊
ＮＵＭ＊Ｉ）＋（＊ＮＵＭ＊ＩＩ）＋（＊ＮＵＭ＊ＩＩ
Ｉ）＋（＊ＮＵＭ＊ＩＶ）＋（＊ＮＵＭ＊Ｖ）／５）の
ような演算が行われる。この演算結果を、＊ＮＵＭ＊Ｔ
と呼ぶ。この＊ＮＵＭ＊Ｔは、異なる時相Ｉ〜Ｖのデー
タＮＵＭの平均値である。
より大きい間隔ΔＴ（異なる時相となるような間隔）毎
の遅延加算（累加平均）が行われる。この結果、（（＊
ＮＵＭ＊Ｉ）＋（＊ＮＵＭ＊ＩＩ）＋（＊ＮＵＭ＊ＩＩ
Ｉ）＋（＊ＮＵＭ＊ＩＶ）＋（＊ＮＵＭ＊Ｖ）／５）の
ような演算が行われる。この演算結果を、＊ＮＵＭ＊Ｔ
と呼ぶ。この＊ＮＵＭ＊Ｔは、異なる時相Ｉ〜Ｖのデー
タＮＵＭの平均値である。
【００５１】更に同様に、遅延加算器４１８にて、Δｔ
より大きい間隔ΔＴ（異なる時相となるような間隔）毎
の遅延加算（累加平均）が行われる。この結果、（（＊
Ｐ＊Ｉ）＋（＊Ｐ＊ＩＩ）＋（＊Ｐ＊ＩＩＩ）＋（＊Ｐ
＊ＩＶ）＋（＊Ｐ＊Ｖ）／５）のような演算が行われ
る。この演算結果を、＊Ｐ＊Ｔと呼ぶ。この＊Ｐ＊Ｔ
は、異なる時相Ｉ〜ＶのデータＰの平均値である。
より大きい間隔ΔＴ（異なる時相となるような間隔）毎
の遅延加算（累加平均）が行われる。この結果、（（＊
Ｐ＊Ｉ）＋（＊Ｐ＊ＩＩ）＋（＊Ｐ＊ＩＩＩ）＋（＊Ｐ
＊ＩＶ）＋（＊Ｐ＊Ｖ）／５）のような演算が行われ
る。この演算結果を、＊Ｐ＊Ｔと呼ぶ。この＊Ｐ＊Ｔ
は、異なる時相Ｉ〜ＶのデータＰの平均値である。
【００５２】このような＊ＤＥＮ＊Ｔ，＊ＮＵＭ＊Ｔ，
＊Ｐ＊Ｔが得られると、平均速度算出部４１４におい
て、ｔａｎ-1（（＊ＮＵＭ＊Ｔ）／（＊ＤＥＮ＊Ｔ））
という処理が実行されて、平均速度＊ｖ＊Ｔが求められ
る。
＊Ｐ＊Ｔが得られると、平均速度算出部４１４におい
て、ｔａｎ-1（（＊ＮＵＭ＊Ｔ）／（＊ＤＥＮ＊Ｔ））
という処理が実行されて、平均速度＊ｖ＊Ｔが求められ
る。
【００５３】また、分散算出部４１５において、（１−
（（＊ＤＥＮ＊Ｔ2 ＋＊ＮＵＭ＊Ｔ 2 ）1/2 ）／＊Ｐ＊
Ｔ）という処理が実行されて、分散＊ｖａｒ＊Ｔが求め
られる。
（（＊ＤＥＮ＊Ｔ2 ＋＊ＮＵＭ＊Ｔ 2 ）1/2 ）／＊Ｐ＊
Ｔ）という処理が実行されて、分散＊ｖａｒ＊Ｔが求め
られる。
【００５４】このようにして得られた平均速度と分散と
を受けた画像処理部５が、それぞれの信号値に応じて赤
／青の色相に輝度を変え、また、この赤／青に緑などの
輝度を変えてた表示を重畳してカラー表示を行う。
を受けた画像処理部５が、それぞれの信号値に応じて赤
／青の色相に輝度を変え、また、この赤／青に緑などの
輝度を変えてた表示を重畳してカラー表示を行う。
【００５５】従って、本実施例では、異なった時相（Ｉ
〜Ｖ）にわたる血流の平均速度と分散との平均値を表示
するものである。このため、時間的な拍動性を表示する
ことができるようになり、血流が拍動流（血流速度や分
散が周期的に変化する流れ）なのか定常流（血流速度や
分散がほぼ一定なような流れ）なのかを瞬時に識別する
ことが可能になる。
〜Ｖ）にわたる血流の平均速度と分散との平均値を表示
するものである。このため、時間的な拍動性を表示する
ことができるようになり、血流が拍動流（血流速度や分
散が周期的に変化する流れ）なのか定常流（血流速度や
分散がほぼ一定なような流れ）なのかを瞬時に識別する
ことが可能になる。
【００５６】尚、上述の実施例では図８に示した時相Ｉ
〜ＶのΔＴ間隔の５時相のデータを遅延加算したが、こ
れ以外の個数のデータを遅延加算するようにしても構わ
ない。
〜ＶのΔＴ間隔の５時相のデータを遅延加算したが、こ
れ以外の個数のデータを遅延加算するようにしても構わ
ない。
【００５７】図２は拍動性の血流速度の時間変化を示し
ており、横軸は時間ｔを示し、縦軸は血流速度ｖを示し
ている。尚、ハッチングで示した部分はバラツキを意味
している。ここでは、１心拍に相当する時間を示してお
り、例えば上述の時相Ｉ，ＩＩ，ＩＩＩ，ＩＶ，Ｖに相
当するタイミングで上述の平均速度と分散とを求める処
理を実行している。
ており、横軸は時間ｔを示し、縦軸は血流速度ｖを示し
ている。尚、ハッチングで示した部分はバラツキを意味
している。ここでは、１心拍に相当する時間を示してお
り、例えば上述の時相Ｉ，ＩＩ，ＩＩＩ，ＩＶ，Ｖに相
当するタイミングで上述の平均速度と分散とを求める処
理を実行している。
【００５８】図３は定常性の血流速度の時間変化を示し
ており、横軸は時間ｔを示し、縦軸は血流速度ｖを示し
ている。尚、ハッチングで示した部分はバラツキを意味
している。ここでは、図２と同じ時間を示しており、例
えば上述の時相Ｉ，ＩＩ，ＩＩＩ，ＩＶ，Ｖに相当する
タイミングで上述の平均速度と分散とを求める処理を実
行している。
ており、横軸は時間ｔを示し、縦軸は血流速度ｖを示し
ている。尚、ハッチングで示した部分はバラツキを意味
している。ここでは、図２と同じ時間を示しており、例
えば上述の時相Ｉ，ＩＩ，ＩＩＩ，ＩＶ，Ｖに相当する
タイミングで上述の平均速度と分散とを求める処理を実
行している。
【００５９】このような図２，図３のような血流速度に
ついて、上述の演算を実行した場合の平均速度及び分散
はそれぞれ図４，図５に示すようになる。すなわち、拍
動流では、平均速度が広い範囲に分布している。一方、
定常流では、平均速度がほぼ一定の値に集中している。
ついて、上述の演算を実行した場合の平均速度及び分散
はそれぞれ図４，図５に示すようになる。すなわち、拍
動流では、平均速度が広い範囲に分布している。一方、
定常流では、平均速度がほぼ一定の値に集中している。
【００６０】以上の実施例によれば、図４，図５のよう
な特性に応じて着色された表示が行われるので、拍動流
であるか定常流であるかを瞬時に識別することができ
る。尚、図１の構成では、平均速度データと分散データ
との双方でΔＴ毎の遅延加算を行うようにしたが、いず
れか一方のみで遅延加算を行うような構成も可能であ
る。この場合にも、遅延加算を行ったいずれか一方の流
速データによりＣＦＭ表示を行うことで定常流か拍動流
かを識別することが可能であり、回路構成が簡略化でき
る利点も有している。
な特性に応じて着色された表示が行われるので、拍動流
であるか定常流であるかを瞬時に識別することができ
る。尚、図１の構成では、平均速度データと分散データ
との双方でΔＴ毎の遅延加算を行うようにしたが、いず
れか一方のみで遅延加算を行うような構成も可能であ
る。この場合にも、遅延加算を行ったいずれか一方の流
速データによりＣＦＭ表示を行うことで定常流か拍動流
かを識別することが可能であり、回路構成が簡略化でき
る利点も有している。
【００６１】以上詳細に説明したように、超音波イメー
ジング表示方法において、Δｔ間隔で得られた運動物体
反射信号から流速データを求める際に、Δｔより大きい
間隔であって異なる時相のΔＴ毎の運動物体反射信号が
遅延加算されたものにより流速データを求め、このよう
にして求められた流速データによりＣＦＭの超音波イメ
ージング表示を行うようにすることで、拍動流か定常流
かを識別可能な状態で表示することが可能になる。
ジング表示方法において、Δｔ間隔で得られた運動物体
反射信号から流速データを求める際に、Δｔより大きい
間隔であって異なる時相のΔＴ毎の運動物体反射信号が
遅延加算されたものにより流速データを求め、このよう
にして求められた流速データによりＣＦＭの超音波イメ
ージング表示を行うようにすることで、拍動流か定常流
かを識別可能な状態で表示することが可能になる。
【００６２】また、超音波イメージング装置において、
Δｔ間隔で得られた運動物体反射信号から流速データを
求める際に、Δｔより大きい間隔であって異なる時相の
ΔＴ毎の運動物体反射信号が遅延加算されたものにより
流速データを求め、このようにして求められた流速デー
タによりＣＦＭの超音波イメージング表示を行うように
することで、拍動流か定常流かを識別可能な状態で表示
することが可能になる。
Δｔ間隔で得られた運動物体反射信号から流速データを
求める際に、Δｔより大きい間隔であって異なる時相の
ΔＴ毎の運動物体反射信号が遅延加算されたものにより
流速データを求め、このようにして求められた流速デー
タによりＣＦＭの超音波イメージング表示を行うように
することで、拍動流か定常流かを識別可能な状態で表示
することが可能になる。
【００６３】そして、超音波イメージング装置におい
て、Δｔ間隔で得られた運動物体反射信号から速度デー
タ及び分散データを求める際に、Δｔより大きい間隔で
あって異なる時相のΔＴ毎の運動物体反射信号が遅延加
算されたものにより速度データ及び分散データを求め、
このようにして求められた速度データ及び分散データに
よりＣＦＭの超音波イメージング表示をすることで、拍
動流か定常流かを識別可能な状態で表示することが可能
になる。
て、Δｔ間隔で得られた運動物体反射信号から速度デー
タ及び分散データを求める際に、Δｔより大きい間隔で
あって異なる時相のΔＴ毎の運動物体反射信号が遅延加
算されたものにより速度データ及び分散データを求め、
このようにして求められた速度データ及び分散データに
よりＣＦＭの超音波イメージング表示をすることで、拍
動流か定常流かを識別可能な状態で表示することが可能
になる。
【００６４】
【発明の効果】以上詳細に説明したように、超音波イメ
ージング表示方法において、Δｔ間隔で得られた運動物
体反射信号から流速データを求める際に、Δｔより大き
い間隔であって異なる時相のΔＴ毎の運動物体反射信号
が遅延加算されたものにより流速データを求め、このよ
うにして求められた流速データによりＣＦＭの超音波イ
メージング表示を行うようにすることで、拍動流か定常
流かを識別可能な状態で表示することが可能になる。
ージング表示方法において、Δｔ間隔で得られた運動物
体反射信号から流速データを求める際に、Δｔより大き
い間隔であって異なる時相のΔＴ毎の運動物体反射信号
が遅延加算されたものにより流速データを求め、このよ
うにして求められた流速データによりＣＦＭの超音波イ
メージング表示を行うようにすることで、拍動流か定常
流かを識別可能な状態で表示することが可能になる。
【００６５】また、超音波イメージング装置において、
Δｔ間隔で得られた運動物体反射信号から流速データを
求める際に、Δｔより大きい間隔であって異なる時相の
ΔＴ毎の運動物体反射信号が遅延加算されたものにより
流速データを求め、このようにして求められた流速デー
タによりＣＦＭの超音波イメージング表示を行うように
することで、拍動流か定常流かを識別可能な状態で表示
することが可能になる。
Δｔ間隔で得られた運動物体反射信号から流速データを
求める際に、Δｔより大きい間隔であって異なる時相の
ΔＴ毎の運動物体反射信号が遅延加算されたものにより
流速データを求め、このようにして求められた流速デー
タによりＣＦＭの超音波イメージング表示を行うように
することで、拍動流か定常流かを識別可能な状態で表示
することが可能になる。
【００６６】そして、超音波イメージング装置におい
て、Δｔ間隔で得られた運動物体反射信号から速度デー
タ及び分散データを求める際に、Δｔより大きい間隔で
あって異なる時相のΔＴ毎の運動物体反射信号が遅延加
算されたものにより速度データ及び分散データを求め、
このようにして求められた速度データ及び分散データに
よりＣＦＭの超音波イメージング表示をすることで、拍
動流か定常流かを識別可能な状態で表示することが可能
になる。
て、Δｔ間隔で得られた運動物体反射信号から速度デー
タ及び分散データを求める際に、Δｔより大きい間隔で
あって異なる時相のΔＴ毎の運動物体反射信号が遅延加
算されたものにより速度データ及び分散データを求め、
このようにして求められた速度データ及び分散データに
よりＣＦＭの超音波イメージング表示をすることで、拍
動流か定常流かを識別可能な状態で表示することが可能
になる。
【図１】本発明の一実施例の超音波イメージング装置の
ブロック図である。
ブロック図である。
【図２】拍動流の速度スペクトラムを示す特性図であ
る。
る。
【図３】定常流の速度スペクトラムを示す特性図であ
る。
る。
【図４】拍動流の平均速度と分散を示す特性図である。
【図５】定常流の平均速度と分散を示す特性図である。
【図６】超音波イメージング装置の全体構成を示す構成
図である。
図である。
【図７】従来のＣＦＭ演算部の回路構成を示す構成図で
ある。
ある。
【図８】拍動流の速度スペクトラムと平均速度，分散の
様子を示す特性図である。
様子を示す特性図である。
４００，４０１ ＭＴＩフィルタ
４０２〜４０４ 遅延部
４０５〜４０８ 乗算器
４０９ 加算器
４１０ 減算器
４１１〜４１３ 遅延加算器
４１４ 平均速度算出部
４１５ 分散算出部
４１６〜４１８ 遅延加算部
Claims (2)
- 【請求項１】 受波信号からＢモードデータを得るＢモ
ード処理部と、 受波信号から運動物体反射信号を求め、同一音線データ
のサンプリング間隔Δｔより大きい間隔ΔＴ毎に所定個
数の運動物体反射信号の累加平均演算を行って速度デー
タ及び分散データを求めるＣＦＭ演算部と、 前記Ｂモード処理部で得られたＢモードデータ並びに前
記ＣＦＭ演算部で累加平均演算をして得られた速度デー
タ及び分散データを用いてカローフローマッピングによ
る超音波イメージング表示用のイメージデータを生成す
る画像処理部とを備えたことを特徴とする超音波イメー
ジング装置。 - 【請求項２】 受波信号を分析してＢモードデータを得
るＢモード処理部と、 受波信号を分析して運動物体反射信号を求め、同一音線
データのサンプリング間隔Δｔより大きい間隔ΔＴ毎に
所定個数の運動物体反射信号の遅延加算を行って速度デ
ータ及び分散データを求めるＣＦＭ演算部と、 前記Ｂモード処理部で得られたＢモードデータ、前記Ｃ
ＦＭ演算部で遅延加算して得られた速度データ及び分散
データを用いてカローフローマッピングによる超音波イ
メージング表示用のイメージデータを生成する画像処理
部とを備えたことを特徴とする超音波イメージング装
置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP29976694A JP3408647B2 (ja) | 1994-12-02 | 1994-12-02 | 超音波イメージング装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP29976694A JP3408647B2 (ja) | 1994-12-02 | 1994-12-02 | 超音波イメージング装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH08154939A JPH08154939A (ja) | 1996-06-18 |
| JP3408647B2 true JP3408647B2 (ja) | 2003-05-19 |
Family
ID=17876715
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP29976694A Expired - Fee Related JP3408647B2 (ja) | 1994-12-02 | 1994-12-02 | 超音波イメージング装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3408647B2 (ja) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP7313841B2 (ja) * | 2019-02-25 | 2023-07-25 | キヤノンメディカルシステムズ株式会社 | 超音波画像診断装置、医用画像処理装置及び医用画像処理プログラム |
-
1994
- 1994-12-02 JP JP29976694A patent/JP3408647B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JPH08154939A (ja) | 1996-06-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| TW565694B (en) | Ultrasound imaging system | |
| US5445156A (en) | Method for adaptively filtering doppler signals using a complex time domain filter | |
| US20050054931A1 (en) | Tracking clutter filter for spectral & audio doppler | |
| JPH0613031B2 (ja) | 超音波血流イメ−ジング装置 | |
| WO1996011413A1 (en) | Quantitative color flow | |
| JP2001178721A (ja) | パケット・データ獲得を用いた超音波フロー撮像における運動の可視化のための方法および装置 | |
| JPH0824678B2 (ja) | 超音波診断装置 | |
| JP2003529422A (ja) | Ｂカラー優先閾値計算を用いた超音波システム | |
| JP3408647B2 (ja) | 超音波イメージング装置 | |
| JP3281435B2 (ja) | 超音波ドプラ診断装置 | |
| JP3388043B2 (ja) | 超音波イメージング表示方法及び超音波イメージング装置 | |
| JP2807131B2 (ja) | 超音波診断装置 | |
| JP3847907B2 (ja) | 超音波ドプラ診断装置 | |
| JP3329589B2 (ja) | 超音波イメージング表示方法及び超音波診断装置 | |
| JP2563656B2 (ja) | 超音波ドプラ映像装置 | |
| JP2703943B2 (ja) | 超音波診断装置 | |
| JP2782229B2 (ja) | 超音波画像診断装置 | |
| JP3465197B2 (ja) | 超音波診断装置 | |
| JP3391578B2 (ja) | 相関装置および流れ情報表示装置 | |
| JPH02215449A (ja) | 超音波診断装置 | |
| JP3187169B2 (ja) | 超音波診断装置 | |
| JPH08206113A (ja) | 超音波イメージング処理方法及び超音波イメージング装置 | |
| JP2567246B2 (ja) | カラ−表示ｍｔｉ装置 | |
| JP3332090B2 (ja) | 超音波診断装置 | |
| JPH0292347A (ja) | カラー超音波診断装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| 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: 20080314 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090314 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100314 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100314 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100314 Year of fee payment: 7 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100314 Year of fee payment: 7 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100314 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110314 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110314 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120314 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120314 Year of fee payment: 9 |
|
| LAPS | Cancellation because of no payment of annual fees |