JPH0464693B2 - - Google Patents
Info
- Publication number
- JPH0464693B2 JPH0464693B2 JP59111913A JP11191384A JPH0464693B2 JP H0464693 B2 JPH0464693 B2 JP H0464693B2 JP 59111913 A JP59111913 A JP 59111913A JP 11191384 A JP11191384 A JP 11191384A JP H0464693 B2 JPH0464693 B2 JP H0464693B2
- Authority
- JP
- Japan
- Prior art keywords
- information
- hue
- blood flow
- power
- modulation
- 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 - Lifetime
Links
- 230000017531 blood circulation Effects 0.000 claims description 36
- 238000001514 detection method Methods 0.000 claims description 13
- 238000002604 ultrasonography Methods 0.000 claims description 11
- 238000000034 method Methods 0.000 description 9
- 238000006243 chemical reaction Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 239000000523 sample Substances 0.000 description 6
- 210000005240 left ventricle Anatomy 0.000 description 4
- 210000000601 blood cell Anatomy 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000001028 reflection method Methods 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 238000004040 coloring Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000000004 hemodynamic effect Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Description
【発明の詳細な説明】
〔発明の技術分野〕
本発明は、超音波反射法によつて得た被検体の
断層像中に超音波ドプラ法によつて得た血流情報
を２次元的に表示可能にした超音波診断装置に関
する。
断層像中に超音波ドプラ法によつて得た血流情報
を２次元的に表示可能にした超音波診断装置に関
する。
心臓内に流れる血流の動態を心臓の断層像と同
時に観測することは臨床的に価値がある。この血
流情報の検出方法としては、距離分類能がすぐれ
ていることで超音波（パルス変調）ドプラ法が主
に用いられている。そして、上記超音波ドプラ法
と超音波反射法とを同時に行なえるようにした超
音波診断装置が具体化され臨床で実用に供されて
いる。
時に観測することは臨床的に価値がある。この血
流情報の検出方法としては、距離分類能がすぐれ
ていることで超音波（パルス変調）ドプラ法が主
に用いられている。そして、上記超音波ドプラ法
と超音波反射法とを同時に行なえるようにした超
音波診断装置が具体化され臨床で実用に供されて
いる。
また近年に至つては、レーダ技術等で用いられ
ているMTI（Moving Target Indication；移動
目標指示装置）を応用して、被検体内の血流を２
次元画像として可視化する技術が提案されている
（例えば、Levy Gerzbeg，“Nonolithic Power
−Spectrum Centroid Detector”in Technical
Peport No.G557−２，Stanford University，
May 1979）。
ているMTI（Moving Target Indication；移動
目標指示装置）を応用して、被検体内の血流を２
次元画像として可視化する技術が提案されている
（例えば、Levy Gerzbeg，“Nonolithic Power
−Spectrum Centroid Detector”in Technical
Peport No.G557−２，Stanford University，
May 1979）。
上記において血流情報としては、平均血流速度
（ドプラ偏移周波数（d））とその分散σ2とを演算
して、これら演算結果を輝度変調として例えばカ
ラー対応させて、例えば心臓内左室においては第
１図に示すＢモード像（断層像）M1と血流情報
（平均血流速、分散）とを、第２図に示すカラー
スケールM2にもとづいた色を上記血流情報に対
応させて、第１図に示すＢモード像と重ね合わせ
て同一表示している。ところが、リアルタイムで
上記分散σ2の演算を正確に行うのは困難であり、
短時間で行なおうとすると多大の誤差をともない
実用的でない。また、例えば左室の駆出血流は上
記方法では流れるように見えない。以下この理由
の第３図及び第４図を参照して説明する。例えば
第３図に示すように左室の内部位Ａ，Ｂ，Ｃにお
いて、この方向に血液が流れていて、駆出開始時
間をｔ＝０とすれば、この左室内部位のＡ，Ｂ，
Ｃでは、第４図ａ，ｂ，ｃに示すように時間−速
度曲線が異なつている。この場合、心臓壁などの
遅い動きを除去するドツプラフイルタ（例えば、
カツトオフ周波数c＝50Hzのハイパスフイルタ
ー）が装置に組込まれているため（第４図中の点
線でcを示す）、tC＜tB＜tAの関係から駆出時には
Ｃ，Ｂ，Ａの順で血流が表現されてしまい、実際
とは異なつて、流れるように見えない。
（ドプラ偏移周波数（d））とその分散σ2とを演算
して、これら演算結果を輝度変調として例えばカ
ラー対応させて、例えば心臓内左室においては第
１図に示すＢモード像（断層像）M1と血流情報
（平均血流速、分散）とを、第２図に示すカラー
スケールM2にもとづいた色を上記血流情報に対
応させて、第１図に示すＢモード像と重ね合わせ
て同一表示している。ところが、リアルタイムで
上記分散σ2の演算を正確に行うのは困難であり、
短時間で行なおうとすると多大の誤差をともない
実用的でない。また、例えば左室の駆出血流は上
記方法では流れるように見えない。以下この理由
の第３図及び第４図を参照して説明する。例えば
第３図に示すように左室の内部位Ａ，Ｂ，Ｃにお
いて、この方向に血液が流れていて、駆出開始時
間をｔ＝０とすれば、この左室内部位のＡ，Ｂ，
Ｃでは、第４図ａ，ｂ，ｃに示すように時間−速
度曲線が異なつている。この場合、心臓壁などの
遅い動きを除去するドツプラフイルタ（例えば、
カツトオフ周波数c＝50Hzのハイパスフイルタ
ー）が装置に組込まれているため（第４図中の点
線でcを示す）、tC＜tB＜tAの関係から駆出時には
Ｃ，Ｂ，Ａの順で血流が表現されてしまい、実際
とは異なつて、流れるように見えない。
本発明は、上記事情に基づいてなされたもの
で、その目的とするところは、超音波断層像中に
血流情報として血流からの散乱パワーに比例した
平均パワー（パワー情報）を表示し、さらには平
均速度と上記平均パワーとを組み合わせて表示す
ることにより、血流の流れの可視化をより正確に
行うことを可能とした超音波診断装置を提供する
ことにある。
で、その目的とするところは、超音波断層像中に
血流情報として血流からの散乱パワーに比例した
平均パワー（パワー情報）を表示し、さらには平
均速度と上記平均パワーとを組み合わせて表示す
ることにより、血流の流れの可視化をより正確に
行うことを可能とした超音波診断装置を提供する
ことにある。
本発明は、被検体内に超音波を送受波して上記
被検体内の断層像情報を得ると共に、上記送受波
による散乱超音波に基づき上記被検体内の血流情
報を得る超音波診断装置において、上記散乱超音
波から血流情報を示すドプラ偏移信号を検出する
ドプラ信号検出手段と、このドプラ信号検出手段
で得たドプラ偏移信号に基づき測定対象の血流の
パワーを示す情報を演算する演算手段と、この演
算手段で演算されたパワー情報にそのパワーの大
きさに応じて色相変調をかける色相変調手段と、
この色相変調手段で色相変調されたパワー情報を
上記断層像に重畳して表示する表示手段とを備え
たことを特徴とし、超音波断層像中に血流情報と
して血流からの散乱パワーに比例したパワー情報
を表示するようにした超音波診断装置を提供する
ものである。
被検体内の断層像情報を得ると共に、上記送受波
による散乱超音波に基づき上記被検体内の血流情
報を得る超音波診断装置において、上記散乱超音
波から血流情報を示すドプラ偏移信号を検出する
ドプラ信号検出手段と、このドプラ信号検出手段
で得たドプラ偏移信号に基づき測定対象の血流の
パワーを示す情報を演算する演算手段と、この演
算手段で演算されたパワー情報にそのパワーの大
きさに応じて色相変調をかける色相変調手段と、
この色相変調手段で色相変調されたパワー情報を
上記断層像に重畳して表示する表示手段とを備え
たことを特徴とし、超音波断層像中に血流情報と
して血流からの散乱パワーに比例したパワー情報
を表示するようにした超音波診断装置を提供する
ものである。
以下本発明に係る超音波診断装置を第５図に示
す一実施例について説明する。尚、本実施例では
超音波送信系、及び受信系の一部は省略してあ
る。第５図において、１０ａ、１０ｂは超音波プ
ローブからの受信信号と、基準信号発生器１１か
らの基準パルス信号0とを混合するミキサーであ
り、ミキサー１０ｂには移相器１２により90°位
相の異ならせた基準パルス信号0が入力される。
１３ａ，１３ｂは夫々上記ミキサー１０ａ，１０
ｂの混合出力を入力し、それをろ波し位相検波信
号として出力するローパスフイルタ（L.P.F）で
ある。この位相検波信号は、被検体内部の種々の
距離におけるドプラ偏移成分を含んだ信号であ
る。１４ａ，１４ｂはローパスフイルタ１３ａ，
１３ｂからの出力信号をＡ／Ｄデイジタル信号に
変換するＡ／Ｄ変換器（Ａ／Ｄ）であり、１５
ａ，１５ｂは上記デイジタル化された位相検波出
力から心臓壁などによる固定反射信号を取除くフ
イルタ回路（MTI filter）である。
す一実施例について説明する。尚、本実施例では
超音波送信系、及び受信系の一部は省略してあ
る。第５図において、１０ａ、１０ｂは超音波プ
ローブからの受信信号と、基準信号発生器１１か
らの基準パルス信号0とを混合するミキサーであ
り、ミキサー１０ｂには移相器１２により90°位
相の異ならせた基準パルス信号0が入力される。
１３ａ，１３ｂは夫々上記ミキサー１０ａ，１０
ｂの混合出力を入力し、それをろ波し位相検波信
号として出力するローパスフイルタ（L.P.F）で
ある。この位相検波信号は、被検体内部の種々の
距離におけるドプラ偏移成分を含んだ信号であ
る。１４ａ，１４ｂはローパスフイルタ１３ａ，
１３ｂからの出力信号をＡ／Ｄデイジタル信号に
変換するＡ／Ｄ変換器（Ａ／Ｄ）であり、１５
ａ，１５ｂは上記デイジタル化された位相検波出
力から心臓壁などによる固定反射信号を取除くフ
イルタ回路（MTI filter）である。
上記Ａ／Ｄ変換器１４ａ，１４ｂ及びフイルタ
回路１５ａ，１５ｂにより、上記位相検波信号が
含有している心臓内の壁などの固定反射による直
流及びドプラ偏移周波数より高い周波数成分を除
去している。
回路１５ａ，１５ｂにより、上記位相検波信号が
含有している心臓内の壁などの固定反射による直
流及びドプラ偏移周波数より高い周波数成分を除
去している。
上記において、フイルタ回路１５ａ，１５ｂは
例えば第６図に示すように加算器３０，３１、シ
フトレジスタ３２及び乗算器３３から構成される
ものであり、シフトレジスタ３２は、１レートパ
ルス間隔に相当する時間分の遅延を行なわせるも
のである。ここで、レートパルス周波数rを5K
Hzとすれば、レートパルス間隔時間は200μsに相
当する。
例えば第６図に示すように加算器３０，３１、シ
フトレジスタ３２及び乗算器３３から構成される
ものであり、シフトレジスタ３２は、１レートパ
ルス間隔に相当する時間分の遅延を行なわせるも
のである。ここで、レートパルス周波数rを5K
Hzとすれば、レートパルス間隔時間は200μsに相
当する。
また、フイルタ回路１５ａ，１５ｂのフイルタ
特性は、横軸を周波数、縦軸を出力信号電圧とし
て表わすと第７図に示されるフイルタ特性であ
り、このようなフイルタ特性はレーダ技術におけ
るMTIを応用したものである。第５図のフイル
タ回路１５ａ，１５ｂを２段設けるとさらに良い
フイルタ特性が達成でき、理論的にはフイルタ回
路１５ａ，１５ｂでの演算を10レートパルス分
（200μs×10）繰り返すことによつて、心臓壁など
からの固定反射信号をとり除くことができる。
特性は、横軸を周波数、縦軸を出力信号電圧とし
て表わすと第７図に示されるフイルタ特性であ
り、このようなフイルタ特性はレーダ技術におけ
るMTIを応用したものである。第５図のフイル
タ回路１５ａ，１５ｂを２段設けるとさらに良い
フイルタ特性が達成でき、理論的にはフイルタ回
路１５ａ，１５ｂでの演算を10レートパルス分
（200μs×10）繰り返すことによつて、心臓壁など
からの固定反射信号をとり除くことができる。
ここで第５図における受信信号Ｓ(t)を下記式(1)
とすると、 Ｓ(t)＝｛〓aocos（ω0t＋ωdoｔ＋φdo）｝ ＋bcos（ω0t＋φb） ……(1) なお、〓aocos（ω0t＋ωdoｔ＋φdo）はドプラ偏
移信号であり、bcos（ω0t＋φb）は固定反射信号
を示している。上記受信信号Ｓ(t)を入力したミキ
サー１０ａ，１０ｂの出力Va″(t)、Vb″(t)は下記
式(2)で示される。
とすると、 Ｓ(t)＝｛〓aocos（ω0t＋ωdoｔ＋φdo）｝ ＋bcos（ω0t＋φb） ……(1) なお、〓aocos（ω0t＋ωdoｔ＋φdo）はドプラ偏
移信号であり、bcos（ω0t＋φb）は固定反射信号
を示している。上記受信信号Ｓ(t)を入力したミキ
サー１０ａ，１０ｂの出力Va″(t)、Vb″(t)は下記
式(2)で示される。
Va″(t)＝Ｓ(t)・cosω0t ……(2)
Vb″(t)＝Ｓ(t)・sinω0t ……(3)
そしてローパスフイルタ１３ａ，１３ｂの出力
Va′(t)、Va′(t)は(1)(2)(3)式から下記式(4)、(5)と
導
かれる。
Va′(t)、Va′(t)は(1)(2)(3)式から下記式(4)、(5)と
導
かれる。
Va′(t)＝〓ao／２cos（ωdoｔ＋φdo）
＋ｂ／２cosφb ……(4)
Vb′(t)＝−〓ao／２sin（ωdoｔ＋φdo）
−ｂ／２sinφb ……(5)
フイルタ１５ａ，１５ｂの出力Va(t)、Vb(t)
は、DC成分がカツトされて下記式(6)、(7)となる。
は、DC成分がカツトされて下記式(6)、(7)となる。
Va(t)＝〓ao／２cos（ωdoｔ＋φdo） ……(6)
Vb(t)＝〓ao／２sin（ωdoｔ＋φdo） ……(7)
第５図において１６は演算回路であり、この演
算回路１６は例えば第８図に示すように、フイル
タ回路１５ａ，１５ｂの出力Va(t)、Vb(t)を取り
込む１レートデイレーライン（レートパルス周波
数rを5KHzとすれば、１レートデイレイ時間は
200μsである）３４，３５、データセレクタ３６，
３７及びこれからの出力を積和する積和回路３８
から構成されている。そしてこの演算回路３８で
は上記出力Va(t)、Vb(t)により、血球からの散乱
パワーに比例した量（以下均等パワーと称する）
Ｐと、平均流速（ドプラ偏移周波数：ｄ＝
ωd／2πを演算する。即ち、式(6)、(7)式の出力Va(t)
、 Vb(t)より、平均パワーは、下記(8)式で求まる。
算回路１６は例えば第８図に示すように、フイル
タ回路１５ａ，１５ｂの出力Va(t)、Vb(t)を取り
込む１レートデイレーライン（レートパルス周波
数rを5KHzとすれば、１レートデイレイ時間は
200μsである）３４，３５、データセレクタ３６，
３７及びこれからの出力を積和する積和回路３８
から構成されている。そしてこの演算回路３８で
は上記出力Va(t)、Vb(t)により、血球からの散乱
パワーに比例した量（以下均等パワーと称する）
Ｐと、平均流速（ドプラ偏移周波数：ｄ＝
ωd／2πを演算する。即ち、式(6)、(7)式の出力Va(t)
、 Vb(t)より、平均パワーは、下記(8)式で求まる。
＝１／Nτ・N〓
〓t=0
｛Va2(t)＋Vb2(t)｝ ……(8)
ｔ＝０、τ、2τ、…、Nτ、τ＝１／n
ここで簡単のために、単一の移動物体からのド
プラ信号を考えて、ｎ＝１（単一ドプラ周波数）
として、 ao＝2a、ωdo＝ωd、ωdo＝０とおくと、上記
(6)、(7)式は、 Va(t)＝acosωdt ……(6)′ Vb(t)＝asinωdt ……(7)′ これにより、上記(8)式は ＝１／Nτ・N〓 〓t=0 a2（cos2ωdt ＋sinb2ωdt）＝a2 ……(8)′ で求まる。
プラ信号を考えて、ｎ＝１（単一ドプラ周波数）
として、 ao＝2a、ωdo＝ωd、ωdo＝０とおくと、上記
(6)、(7)式は、 Va(t)＝acosωdt ……(6)′ Vb(t)＝asinωdt ……(7)′ これにより、上記(8)式は ＝１／Nτ・N〓 〓t=0 a2（cos2ωdt ＋sinb2ωdt）＝a2 ……(8)′ で求まる。
平均流速（ドプラ偏移周波数：ｄ＝ωd／2π）
は、複素相関演算を行うことにより、
Va(t)、Vb(t)の複素自己相関は、
Ｃ（τ）＝１／NτN〓
〓
〓t=0
〔Va（ｔ＋τ）＋jVb（ｔ＋τ）〕〔Va(t)−jVb(t)
〕 ＝１／NτN〓 〓 〓t=0 〔｛Va（ｔ＋τ）Va(t)＋Vb（ｔ＋τ）Vb(t)｝＋ｊ
｛Vb（ｔ＋τ）Va(t) −Va（ｔ＋τ）Vb(t)｝〕＝CR＋jCI ……(9) ただし、ｔ＝０、τ、2τ、…、Nτ τ＝１／r (6)′、(7)′を代入して CR＝１／NτN〓 〓 〓t=0 a2〔cosωd（ｔ＋τ）cosωdt＋sinωd（ｔ＋τ）si
nωdt〕 ＝１／NτN〓 〓t=0 a2cosωdτ＝a2cosωdτ ……(10) CI＝１／NτN〓 〓 〓t=0 a2〔sinωd（ｔ＋τ）cosωdt−cosωd（ｔ＋τ）si
nωdt〕 ＝１／NτN〓 〓t=0 a2sinωdτ＝a2sinωdτ ……(11) (9)、(10)、(11)式から ｄ／r＝tan-1（CI／CR） ∴ｄ＝r／2π・tan-1（CI／CR） が求められる。ここで、CI，CRの符号から血流の
方向（プローブに近づく流れか、プローブから遠
ざかる流れか）の判定も可能である。
〕 ＝１／NτN〓 〓 〓t=0 〔｛Va（ｔ＋τ）Va(t)＋Vb（ｔ＋τ）Vb(t)｝＋ｊ
｛Vb（ｔ＋τ）Va(t) −Va（ｔ＋τ）Vb(t)｝〕＝CR＋jCI ……(9) ただし、ｔ＝０、τ、2τ、…、Nτ τ＝１／r (6)′、(7)′を代入して CR＝１／NτN〓 〓 〓t=0 a2〔cosωd（ｔ＋τ）cosωdt＋sinωd（ｔ＋τ）si
nωdt〕 ＝１／NτN〓 〓t=0 a2cosωdτ＝a2cosωdτ ……(10) CI＝１／NτN〓 〓 〓t=0 a2〔sinωd（ｔ＋τ）cosωdt−cosωd（ｔ＋τ）si
nωdt〕 ＝１／NτN〓 〓t=0 a2sinωdτ＝a2sinωdτ ……(11) (9)、(10)、(11)式から ｄ／r＝tan-1（CI／CR） ∴ｄ＝r／2π・tan-1（CI／CR） が求められる。ここで、CI，CRの符号から血流の
方向（プローブに近づく流れか、プローブから遠
ざかる流れか）の判定も可能である。
なお、上記式(8)乃至(12)の演算においては、最初
の演算でデータセレクタ３６，３７はVa(o)、Va
(o)を選択し、積和回路３８にてVa2(o)を演算す
る。次にデータセレクタ３６，３７は、Vb(o)、
Vb(o)を選択しVa2(o)＋Vb2(o)を演算することによ
り、順次演算を行ない(8)式の平均パワーＰを求め
ることができる。
の演算でデータセレクタ３６，３７はVa(o)、Va
(o)を選択し、積和回路３８にてVa2(o)を演算す
る。次にデータセレクタ３６，３７は、Vb(o)、
Vb(o)を選択しVa2(o)＋Vb2(o)を演算することによ
り、順次演算を行ない(8)式の平均パワーＰを求め
ることができる。
次にデータセレクタ３６，３７は、Va(o)Va
（１／r）を選択し、Va(o)Va（１／r）を演算
し、順次Vb(o)Vb（１／r）を演算し式(9)のリア
ルパートVa（１／r）Va(o)＋Vb（１／r）Vb(o)
を計算し、次にイマジナリーパートVb（１／r）
Va(o)−Va（１／r）Vb(o)を順次求めることがで
きる。
（１／r）を選択し、Va(o)Va（１／r）を演算
し、順次Vb(o)Vb（１／r）を演算し式(9)のリア
ルパートVa（１／r）Va(o)＋Vb（１／r）Vb(o)
を計算し、次にイマジナリーパートVb（１／r）
Va(o)−Va（１／r）Vb(o)を順次求めることがで
きる。
第５図において２１は受信信号を検波する検波
回路であり、２２はこの検波回路２１の出力をデ
ジタル量に変換するＡ／Ｄ変換回路であり、これ
ら信号処理系により超音波断層像（Ｂモード像）
が生成される。
回路であり、２２はこの検波回路２１の出力をデ
ジタル量に変換するＡ／Ｄ変換回路であり、これ
ら信号処理系により超音波断層像（Ｂモード像）
が生成される。
第５図において１７は上記断層像情報と、平均
パワー及び平均流速（ドプラ偏移周波数ｄ
＝ωd／2πとを蓄えるフレームメモリ（FM）である。
パワー及び平均流速（ドプラ偏移周波数ｄ
＝ωd／2πとを蓄えるフレームメモリ（FM）である。
１８はフレームメモリ１７からフレーム毎に平
均パワー及び平均流速の情報を取出してカラ
ー画像信号に変換するカラー変換回路である。
均パワー及び平均流速の情報を取出してカラ
ー画像信号に変換するカラー変換回路である。
このカラー変換回路１８では輝度変調として、
第９図に示すように、平均流速（ドプラ偏移周
波数）を例えばが正の場合（順流）赤系の
色で表示し、負の場合（逆流）には青系の色で表
示し、の大きさをその輝度変化に対応させる。
また、平均パワーは、例えばが小さいほど緑
色をまぜてその色相を変化に対応させる。
第９図に示すように、平均流速（ドプラ偏移周
波数）を例えばが正の場合（順流）赤系の
色で表示し、負の場合（逆流）には青系の色で表
示し、の大きさをその輝度変化に対応させる。
また、平均パワーは、例えばが小さいほど緑
色をまぜてその色相を変化に対応させる。
第５図において１９はカラー変換回路１８によ
りカラー処理された断層像情報と平均パワーＰ及
び平均流速ｖの情報とをアナログ量に変換する
Ｄ／Ａ変換器であり、２０はＤ／Ａ変換器１９か
らの断層像と平均パワーＰ及び平均流速ｖとを重
畳して表示するTVモニタである。
りカラー処理された断層像情報と平均パワーＰ及
び平均流速ｖの情報とをアナログ量に変換する
Ｄ／Ａ変換器であり、２０はＤ／Ａ変換器１９か
らの断層像と平均パワーＰ及び平均流速ｖとを重
畳して表示するTVモニタである。
上記の如く構成された本実施例の装置であれば
以下の如くの作用効果が得られる。
以下の如くの作用効果が得られる。
本実施例では、従来の装置における平均流速ｖ
（ドプラ偏移周波数d）情報に加えて、平均パワ
ーＰ情報をその度合に応じて色付け（輝度変調）
して断層像に重畳して表示するようにしている。
（ドプラ偏移周波数d）情報に加えて、平均パワ
ーＰ情報をその度合に応じて色付け（輝度変調）
して断層像に重畳して表示するようにしている。
この場合、心臓内の血流は、時間的にコヒーレ
ンシーがあり、ある１サンプルボリユーム内での
血球の密度は時々刻々変化しているものであり、
その血球からの平均パワーＰも時々刻々変化して
いる。従つて、その変化を表示していけば、血流
の流れを追うことができ、よつて血行動態を極め
て正確に把握することが可能となる。第９図に示
す以外の血流情報の表現法として第１０図に示す
ように順流（プローブに近づく流れ：dが正の
場合）は赤系の色で、逆流（プローブから遠ざか
る流れ：dが負の場合）は青系の色で表現し、
そのカラーの輝度を平均パワーＰの大きさに対応
させても良い。この他に本発明はその要旨を逸脱
しない範囲で種々変形して実施できる。
ンシーがあり、ある１サンプルボリユーム内での
血球の密度は時々刻々変化しているものであり、
その血球からの平均パワーＰも時々刻々変化して
いる。従つて、その変化を表示していけば、血流
の流れを追うことができ、よつて血行動態を極め
て正確に把握することが可能となる。第９図に示
す以外の血流情報の表現法として第１０図に示す
ように順流（プローブに近づく流れ：dが正の
場合）は赤系の色で、逆流（プローブから遠ざか
る流れ：dが負の場合）は青系の色で表現し、
そのカラーの輝度を平均パワーＰの大きさに対応
させても良い。この他に本発明はその要旨を逸脱
しない範囲で種々変形して実施できる。
以上述べたように本発明によれば、被検体内に
超音波を送受波して上記被検体内の断層像情報を
得ると共に、上記送受波による散乱超音波に基づ
き上記被検体内の血流情報を得る超音波診断装置
において、上記散乱超音波から血流情報を示すド
プラ偏移信号を検出するドプラ信号検出手段と、
このドプラ信号検出手段で得たドプラ偏移信号に
基づき測定対象の血流のパワーを示す情報を演算
する演算手段と、この演算手段で演算されたパワ
ー情報にそのパワーの大きさに応じて色相変調を
かける色相変調手段と、この色相変調手段で色相
変調されたパワー情報を上記断層像に重畳して表
示する表示手段とを備えたことを特徴とし、超音
波断層像中に血流情報として血流からの散乱パワ
ーに比例したパワー情報を表示するようにしたの
で、血流の流れの可視化をより正確に行うことを
可能とした超音波診断装置を提供することができ
る。
超音波を送受波して上記被検体内の断層像情報を
得ると共に、上記送受波による散乱超音波に基づ
き上記被検体内の血流情報を得る超音波診断装置
において、上記散乱超音波から血流情報を示すド
プラ偏移信号を検出するドプラ信号検出手段と、
このドプラ信号検出手段で得たドプラ偏移信号に
基づき測定対象の血流のパワーを示す情報を演算
する演算手段と、この演算手段で演算されたパワ
ー情報にそのパワーの大きさに応じて色相変調を
かける色相変調手段と、この色相変調手段で色相
変調されたパワー情報を上記断層像に重畳して表
示する表示手段とを備えたことを特徴とし、超音
波断層像中に血流情報として血流からの散乱パワ
ーに比例したパワー情報を表示するようにしたの
で、血流の流れの可視化をより正確に行うことを
可能とした超音波診断装置を提供することができ
る。
第１図は超音波パルス反射波法による断層像
（Ｂモード像）の一例を示す図、第２図は超音波
パルスドプラ法による血流情報の表示の一例を示
す図、第３図及び第４図は血行状態を説明するた
めの図、第５図は本発明の一実施例を示すブロツ
ク図、第６図は同実施例におけるフイルタ回路の
詳細を示すブロツク図、第７図は同実施例におけ
るフイルタ回路の周波数特性を示す図、第８図は
同実施例における演算回路の詳細を示すブロツク
図、第９図、第１０図は夫々同実施例におけるカ
ラー変換回路の各情報と色付けとの反応を示す図
である。 １０ａ，１０ｂ……ミキサー、１１……基準信
号発生器、１２……移相器、１３ａ，１３ｂ……
ローパスフイルタ、１４ａ，１４ｂ……Ａ／Ｄ変
換器、１５ａ，１５ｂ……（MTI）フイルタ回
路、１６……演算回路、１７……フレームメモ
リ、１８……カラー変換回路、１９……Ｄ／Ａ変
換器、２０……TVモニタ。
（Ｂモード像）の一例を示す図、第２図は超音波
パルスドプラ法による血流情報の表示の一例を示
す図、第３図及び第４図は血行状態を説明するた
めの図、第５図は本発明の一実施例を示すブロツ
ク図、第６図は同実施例におけるフイルタ回路の
詳細を示すブロツク図、第７図は同実施例におけ
るフイルタ回路の周波数特性を示す図、第８図は
同実施例における演算回路の詳細を示すブロツク
図、第９図、第１０図は夫々同実施例におけるカ
ラー変換回路の各情報と色付けとの反応を示す図
である。 １０ａ，１０ｂ……ミキサー、１１……基準信
号発生器、１２……移相器、１３ａ，１３ｂ……
ローパスフイルタ、１４ａ，１４ｂ……Ａ／Ｄ変
換器、１５ａ，１５ｂ……（MTI）フイルタ回
路、１６……演算回路、１７……フレームメモ
リ、１８……カラー変換回路、１９……Ｄ／Ａ変
換器、２０……TVモニタ。
Claims (1)
- 【特許請求の範囲】 １ 被検体内に超音波を送受波して上記被検体内
の断層像情報を得ると共に、上記送受波による散
乱超音波に基づき上記被検体内の血流情報を得る
超音波診断装置において、上記散乱超音波から血
流情報を示すドプラ偏移信号を検出するドプラ信
号検出手段と、このドプラ信号検出手段で得たド
プラ偏移信号に基づき測定対象の血流のパワーを
示す情報を演算する演算手段と、この演算手段で
演算されたパワー情報にそのパワーの大きさに応
じて色相変調をかける色相変調手段と、この色相
変調手段で色相変調されたパワー情報を上記断層
像に重畳して表示する表示手段とを備えたことを
特徴とする超音波診断装置。 ２ 前記ドプラ信号検出手段で得たドプラ偏移信
号に基づき測定対象の血流の平均速度情報を演算
する手段と、この演算された平均速度情報にその
速度に応じて輝度変調をかける輝度変調手段とを
備え、前記表示手段は前記輝度変調手段で輝度変
調された平均速度情報を前記色相変調されたパワ
ー情報とともに前記断層像に重畳して表示するこ
とを特徴とする請求項１記載の超音波診断装置。 ３ 前記色相変調手段は前記平均速度情報にその
速度方向に応じて色相変調をかけるとともに、前
記表示手段は輝度変調及び色相変調された平均速
度情報を前記色相変調されたパワーとともに前記
断層像に重畳して表示することを特徴とする請求
項２記載の超音波診断装置。 ４ 前記色相変調手段は前記パワー情報にその測
定対象の血流方向に応じて色相変調をかけること
を特徴とする請求項１または２記載の超音波診断
装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP11191384A JPS60256440A (ja) | 1984-05-31 | 1984-05-31 | 超音波診断装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP11191384A JPS60256440A (ja) | 1984-05-31 | 1984-05-31 | 超音波診断装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS60256440A JPS60256440A (ja) | 1985-12-18 |
| JPH0464693B2 true JPH0464693B2 (ja) | 1992-10-15 |
Family
ID=14573257
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP11191384A Granted JPS60256440A (ja) | 1984-05-31 | 1984-05-31 | 超音波診断装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS60256440A (ja) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6359940A (ja) * | 1986-08-29 | 1988-03-15 | 横河メディカルシステム株式会社 | 超音波診断装置 |
| JPS6359938A (ja) * | 1986-08-29 | 1988-03-15 | 横河メディカルシステム株式会社 | 超音波診断装置 |
| JPS63105744A (ja) * | 1986-10-24 | 1988-05-11 | 富士通株式会社 | 超音波診断装置 |
| JPH0284944A (ja) * | 1988-09-22 | 1990-03-26 | Yokogawa Medical Syst Ltd | 超音波診断装置 |
| JPH02121642A (ja) * | 1988-10-31 | 1990-05-09 | Yokogawa Medical Syst Ltd | 超音波診断装置 |
| JP6801337B2 (ja) * | 2016-09-28 | 2020-12-16 | コニカミノルタ株式会社 | 超音波診断装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS51136443A (en) * | 1975-05-21 | 1976-11-25 | Yokogawa Hewlett Packard Ltd | Device for detection and display of moving objects |
| JPS5920820A (ja) * | 1982-07-28 | 1984-02-02 | Aloka Co Ltd | 超音波血流画像形成装置 |
-
1984
- 1984-05-31 JP JP11191384A patent/JPS60256440A/ja active Granted
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS51136443A (en) * | 1975-05-21 | 1976-11-25 | Yokogawa Hewlett Packard Ltd | Device for detection and display of moving objects |
| JPS5920820A (ja) * | 1982-07-28 | 1984-02-02 | Aloka Co Ltd | 超音波血流画像形成装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPS60256440A (ja) | 1985-12-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JPH01310648A (ja) | 超音波血流イメージング装置 | |
| JPH0613031B2 (ja) | 超音波血流イメ−ジング装置 | |
| JPH0216139B2 (ja) | ||
| JPS63317137A (ja) | 超音波血流イメ−ジング装置 | |
| CA1267968A (en) | Ultrasonic device for measuring acceleration of moving reflective member | |
| JPH0464693B2 (ja) | ||
| Brandestini et al. | Blood flow imaging using a discrete-time frequency meter | |
| JPH0418859B2 (ja) | ||
| JPH02193650A (ja) | 超音波ドプラ診断装置 | |
| Kasai et al. | Real-time two-dimensional blood flow imaging using ultrasound Doppler | |
| JP3388043B2 (ja) | 超音波イメージング表示方法及び超音波イメージング装置 | |
| JP2563656B2 (ja) | 超音波ドプラ映像装置 | |
| JPH05212032A (ja) | 超音波ドプラ診断装置 | |
| JPH0741038B2 (ja) | 超音波診断装置における血流動態表示方法 | |
| JP2703943B2 (ja) | 超音波診断装置 | |
| JPH01244738A (ja) | 超音波診断装置 | |
| JP3408647B2 (ja) | 超音波イメージング装置 | |
| JP3406096B2 (ja) | 超音波診断装置 | |
| JPH03272751A (ja) | 超音波診断装置 | |
| JPH02193649A (ja) | 超音波ドプラ診断装置 | |
| JPH0543381B2 (ja) | ||
| JPH0741039B2 (ja) | カラー超音波診断装置 | |
| JPH0199537A (ja) | 超音波診断装置 | |
| JPH0444745A (ja) | 超音波画像形成装置 | |
| JPH06225880A (ja) | 超音波血流イメージング装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EXPY | Cancellation because of completion of term |