JP3352211B2 - 超音波ドプラ診断装置 - Google Patents
超音波ドプラ診断装置Info
- Publication number
- JP3352211B2 JP3352211B2 JP3999794A JP3999794A JP3352211B2 JP 3352211 B2 JP3352211 B2 JP 3352211B2 JP 3999794 A JP3999794 A JP 3999794A JP 3999794 A JP3999794 A JP 3999794A JP 3352211 B2 JP3352211 B2 JP 3352211B2
- Authority
- JP
- Japan
- Prior art keywords
- doppler shift
- histogram
- frequency
- circuit
- doppler
- 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
- 238000002604 ultrasonography Methods 0.000 title claims 2
- 238000001514 detection method Methods 0.000 claims description 6
- 239000000523 sample Substances 0.000 description 38
- 238000001228 spectrum Methods 0.000 description 19
- 238000010586 diagram Methods 0.000 description 9
- 238000005259 measurement Methods 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 230000003111 delayed effect Effects 0.000 description 4
- 230000017531 blood circulation Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 238000000034 method Methods 0.000 description 1
- 230000004218 vascular function Effects 0.000 description 1
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Description
【０００１】
【産業上の利用分野】本発明は、超音波ドプラ診断装
置、特にドプラ偏移周波数について周波数解析を行う超
音波ドプラ診断装置に関する。
置、特にドプラ偏移周波数について周波数解析を行う超
音波ドプラ診断装置に関する。
【０００２】
【従来の技術】例えば、血管機能を診断する場合には、
いわゆるパルスドプラ法を用いて血流速度分布の計測が
行われる。すなわち、図６に示すように、超音波ビーム
１０を走査して形成される走査面１２内にサンプル点Ｓ
を設定し、かかるサンプル点Ｓの受信信号に対してドプ
ラ演算及び周波数解析演算を行うことにより、図７に示
すようなパワースペクトラムを得ている。ここで、図７
のグラフの横軸はドプラ偏移周波数ψｄであり、これは
血流速度に相当する。縦軸はパワーＰである。
いわゆるパルスドプラ法を用いて血流速度分布の計測が
行われる。すなわち、図６に示すように、超音波ビーム
１０を走査して形成される走査面１２内にサンプル点Ｓ
を設定し、かかるサンプル点Ｓの受信信号に対してドプ
ラ演算及び周波数解析演算を行うことにより、図７に示
すようなパワースペクトラムを得ている。ここで、図７
のグラフの横軸はドプラ偏移周波数ψｄであり、これは
血流速度に相当する。縦軸はパワーＰである。
【０００３】実際の表示においては、そのスペクトラム
の経時的変動を観察するため、図８に示すようなパワー
スペクトラム表示がなされる。ここで、横軸は時間ｔで
あり、縦軸はドプラ偏移周波数ψｄである。更に、画面
上の輝度軸が各周波数のパワーの大きさに対応してい
る。
の経時的変動を観察するため、図８に示すようなパワー
スペクトラム表示がなされる。ここで、横軸は時間ｔで
あり、縦軸はドプラ偏移周波数ψｄである。更に、画面
上の輝度軸が各周波数のパワーの大きさに対応してい
る。
【０００４】図９には、周波数解析を行うことのできる
従来の超音波ドプラ診断装置の概略構成が示されてい
る。
従来の超音波ドプラ診断装置の概略構成が示されてい
る。
【０００５】受信信号１００は、直交検波回路１６で直
交検波され、サンプルホールド回路１７によって対象部
位のみの信号が選択された後、フィルタ１９、２１によ
り波形整形される。それがＡ／Ｄ変換器１８でデジタル
信号に変換され、そのデジタル信号に対して、ＦＦＴ演
算器２０にて複素ＦＦＴ演算による周波数解析が行わ
れ、これにて上述のパワースペクトラムが求まる。そし
て、それが図７に示した表示形式等で表示器２２に表示
される。なお、二次元ドプラ画像を併せて表示する場合
には、例えば自己相関回路２４により、自己相関演算が
実行される。
交検波され、サンプルホールド回路１７によって対象部
位のみの信号が選択された後、フィルタ１９、２１によ
り波形整形される。それがＡ／Ｄ変換器１８でデジタル
信号に変換され、そのデジタル信号に対して、ＦＦＴ演
算器２０にて複素ＦＦＴ演算による周波数解析が行わ
れ、これにて上述のパワースペクトラムが求まる。そし
て、それが図７に示した表示形式等で表示器２２に表示
される。なお、二次元ドプラ画像を併せて表示する場合
には、例えば自己相関回路２４により、自己相関演算が
実行される。
【０００６】
【発明が解決しようとする課題】しかしながら、ＦＦＴ
による周波数解析は複雑な計算が必要となり、迅速な演
算が困難であると共に回路規模が大きくなるという問題
があった。よって、簡易な構成で周波数解析を実現する
手段が要望されていた。
による周波数解析は複雑な計算が必要となり、迅速な演
算が困難であると共に回路規模が大きくなるという問題
があった。よって、簡易な構成で周波数解析を実現する
手段が要望されていた。
【０００７】本発明は、上記従来の課題に鑑みなされた
ものであり、その目的は、簡易な構成でドプラ偏移周波
数の周波数解析を実現できる装置を提供することにあ
る。
ものであり、その目的は、簡易な構成でドプラ偏移周波
数の周波数解析を実現できる装置を提供することにあ
る。
【０００８】
【課題を解決するための手段】上記目的を達成するため
に、請求項１記載の発明は、検波後の受信信号間で位相
差分演算を行い、ドプラ偏移情報を求める位相差分演算
回路と、複数個の前記ドプラ偏移情報に基づき、ドプラ
偏移の大きさ毎の度数をグラフ化したヒストグラムを作
成するヒストグラム作成回路と、を含むことを特徴とす
る。
に、請求項１記載の発明は、検波後の受信信号間で位相
差分演算を行い、ドプラ偏移情報を求める位相差分演算
回路と、複数個の前記ドプラ偏移情報に基づき、ドプラ
偏移の大きさ毎の度数をグラフ化したヒストグラムを作
成するヒストグラム作成回路と、を含むことを特徴とす
る。
【０００９】また、請求項２記載の発明は、前記複数個
のドプラ偏移情報は、同一サンプル点に対する複数回の
送受波により得られることを特徴とする。
のドプラ偏移情報は、同一サンプル点に対する複数回の
送受波により得られることを特徴とする。
【００１０】また、請求項３記載の発明は、前記複数個
のドプラ偏移情報は、同一超音波ビーム上に隣接設定さ
れたサンプル点群に対する送受波毎に得られることを特
徴とする。
のドプラ偏移情報は、同一超音波ビーム上に隣接設定さ
れたサンプル点群に対する送受波毎に得られることを特
徴とする。
【００１１】更に、請求項４記載の発明は、前記受信信
号の振幅値を演算する振幅値演算回路を含み、前記ヒス
トグラム作成回路は、振幅値の大きさに基づき重み付け
を行いつつ前記ヒストグラムを作成することを特徴とす
る。
号の振幅値を演算する振幅値演算回路を含み、前記ヒス
トグラム作成回路は、振幅値の大きさに基づき重み付け
を行いつつ前記ヒストグラムを作成することを特徴とす
る。
【００１２】
【作用】本発明によれば、位相差分演算回路において位
相差分演算が行われ、これによりドプラ偏移情報（ドプ
ラシフト周波数）が求められる。すなわち、超音波パル
スの送信繰り返し周波数の１周期分遅延させた受信信号
と遅延させない受信信号の差分が演算される。
相差分演算が行われ、これによりドプラ偏移情報（ドプ
ラシフト周波数）が求められる。すなわち、超音波パル
スの送信繰り返し周波数の１周期分遅延させた受信信号
と遅延させない受信信号の差分が演算される。
【００１３】そして、ヒストグラム作成回路において、
計測対象（生体内運動体）についての複数個のドプラ偏
移情報から、各ドプラ偏移周波数毎にその個数をカウン
トすることによって頻度グラフ（ヒストグラム）が作成
される。
計測対象（生体内運動体）についての複数個のドプラ偏
移情報から、各ドプラ偏移周波数毎にその個数をカウン
トすることによって頻度グラフ（ヒストグラム）が作成
される。
【００１４】その場合、複数個のドプラ偏移情報は、計
測対象についてサンプル点を１点に設定した場合には、
同一サンプル点に対する複数回の超音波パルスの送受波
により得られ、サンプル点を複数個に設定した場合に
は、空間的分布情報も加味されたドプラ偏移情報が得ら
れる。すなわち、前者の場合には、位置的特定性を優先
して、複数回の送受波で１個のヒストグラムを形成で
き、後者の場合には、更に空間的な速度分布をも意味す
るヒストグラムを形成できる。
測対象についてサンプル点を１点に設定した場合には、
同一サンプル点に対する複数回の超音波パルスの送受波
により得られ、サンプル点を複数個に設定した場合に
は、空間的分布情報も加味されたドプラ偏移情報が得ら
れる。すなわち、前者の場合には、位置的特定性を優先
して、複数回の送受波で１個のヒストグラムを形成で
き、後者の場合には、更に空間的な速度分布をも意味す
るヒストグラムを形成できる。
【００１５】このように形成されたヒストグラムにおい
て、各ドプラ偏移の大きさ毎の頻度は、ドプラ偏移が取
り得る確率に相当するものであるため、結果として、ヒ
ストグラムはパワースペクトラムに相当するものといえ
る。
て、各ドプラ偏移の大きさ毎の頻度は、ドプラ偏移が取
り得る確率に相当するものであるため、結果として、ヒ
ストグラムはパワースペクトラムに相当するものといえ
る。
【００１６】ただし、その仮定は、ドプラ偏移情報の抽
出時における各受信信号の振幅値（レベル）が一定とみ
なせる場合に成り立ち、各受信信号の振幅値が一定とみ
なせない場合には振幅値の大きさに基づき、ヒストグラ
ム作成時の度数に重み付けを行う必要がある。
出時における各受信信号の振幅値（レベル）が一定とみ
なせる場合に成り立ち、各受信信号の振幅値が一定とみ
なせない場合には振幅値の大きさに基づき、ヒストグラ
ム作成時の度数に重み付けを行う必要がある。
【００１７】
【実施例】以下、本発明の好適な実施例を図面に基づい
て説明する。
て説明する。
【００１８】図１には、本発明に係る超音波ドプラ診断
装置の全体構成が示されている。図１において、超音波
探触子（図示せず）から送波された超音波パルスは、血
流等の運動体にて反射する際、ドプラシフトを受け、そ
の後、超音波探触子にて受波される。超音波探触子から
の受信信号１００は、受信回路２６にて増幅等が施され
た後、直交検波回路の一部を成す一対のミキサ２８に入
力される。そして、参照信号を用いた直交検波により、
受信信号は実部信号及び虚部信号から成る複素信号に変
換される。
装置の全体構成が示されている。図１において、超音波
探触子（図示せず）から送波された超音波パルスは、血
流等の運動体にて反射する際、ドプラシフトを受け、そ
の後、超音波探触子にて受波される。超音波探触子から
の受信信号１００は、受信回路２６にて増幅等が施され
た後、直交検波回路の一部を成す一対のミキサ２８に入
力される。そして、参照信号を用いた直交検波により、
受信信号は実部信号及び虚部信号から成る複素信号に変
換される。
【００１９】サンプルホールド回路３０は、サンプル点
の受信信号のみを抽出するゲート回路としての機能を有
し、具体的には、例えば１μｓ程度に設定されるゲート
期間内の受信信号を積分して、この積分値をホールド出
力している。なお、ゲート期間は超音波ビーム上のサン
プルボリューム（サンプル点の大きさ）を定めるもので
あり、サンプルボリュームが２ｍｍの場合、ゲート期間
は例えば１．３μｓに設定される。
の受信信号のみを抽出するゲート回路としての機能を有
し、具体的には、例えば１μｓ程度に設定されるゲート
期間内の受信信号を積分して、この積分値をホールド出
力している。なお、ゲート期間は超音波ビーム上のサン
プルボリューム（サンプル点の大きさ）を定めるもので
あり、サンプルボリュームが２ｍｍの場合、ゲート期間
は例えば１．３μｓに設定される。
【００２０】ウオールモーションフィルタ３２は、生体
内低速運動体（例えば、心臓壁）の信号を除去するもの
であり、一対のハイパスフィルタで構成される。ローパ
スフィルタ３４は、受信信号中のベースバンド領域の信
号を取り出す回路であり、その後、受信信号は、Ａ／Ｄ
変換器３６においてデジタル信号に変換される。ここま
での構成はほぼ従来装置の構成と同様である。
内低速運動体（例えば、心臓壁）の信号を除去するもの
であり、一対のハイパスフィルタで構成される。ローパ
スフィルタ３４は、受信信号中のベースバンド領域の信
号を取り出す回路であり、その後、受信信号は、Ａ／Ｄ
変換器３６においてデジタル信号に変換される。ここま
での構成はほぼ従来装置の構成と同様である。
【００２１】デジタル化された複素信号は、逆正接回路
（ＡＴＮ回路）３８に入力され、複素信号の逆正接（ｔ
ａｎ−１）を計算することにより、各受信信号の位相角
度ψが求められる。
（ＡＴＮ回路）３８に入力され、複素信号の逆正接（ｔ
ａｎ−１）を計算することにより、各受信信号の位相角
度ψが求められる。
【００２２】すなわち、逆正接回路３８の入力は、 ｓｉｎ（２πｆｄｎＴ） ｃｏｓ（２πｆｄｎＴ） であり（ただし、ｆｄはドプラ偏移周波数）、その出力
ψ（位相角度）は、 ψ＝２πｆｄｎＴ となる。なお、逆正接回路３８はＲＯＭなどで構成でき
る。
ψ（位相角度）は、 ψ＝２πｆｄｎＴ となる。なお、逆正接回路３８はＲＯＭなどで構成でき
る。
【００２３】そして、その出力ψに対して、位相差分演
算回路４０にて位相差分演算が実行される。ここで、位
相差分演算回路４０は、超音波パルスの送信繰り返し周
期の１周期分だけ受信信号を遅延させるディレーライン
４２と、遅延された受信信号ψ１と遅延されない受信信
号ψ０との位相差分を演算する加算器４４と、で構成さ
れている。なお、実際にはディレーライン４２はＲＡＭ
などで構成できる。
算回路４０にて位相差分演算が実行される。ここで、位
相差分演算回路４０は、超音波パルスの送信繰り返し周
期の１周期分だけ受信信号を遅延させるディレーライン
４２と、遅延された受信信号ψ１と遅延されない受信信
号ψ０との位相差分を演算する加算器４４と、で構成さ
れている。なお、実際にはディレーライン４２はＲＡＭ
などで構成できる。
【００２４】すなわち、ディレーライン４２は、 Ｔ＝１／ＰＲＦ （ＰＲＦ：送信繰り返し周波数） の遅延量を有し、差分器４４では位相の差分が演算さ
れ、つまり、 ψ１＝２πｆｄ（ｎ＋１）Ｔ ψ０＝２πｆｄｎＴ の差分ψｄが次のように演算される。
れ、つまり、 ψ１＝２πｆｄ（ｎ＋１）Ｔ ψ０＝２πｆｄｎＴ の差分ψｄが次のように演算される。
【００２５】ψｄ＝ψ１−ψ０ ＝２πｆｄＴ このψｄがドプラ偏移情報すなわちドプラ偏移周波数に
対応する信号である。このドプラ偏移情報ψｄは、次の
ヒストグラム作成回路４６に入力される。ここで、ヒス
トグラム作成回路４６は、アキュムレータ４８と、移動
平均回路５０と、で構成される。ここで、アキュムレー
タ４８は、各ドプラ偏移周波数毎に設けられた複数のカ
ウンタで構成される。
対応する信号である。このドプラ偏移情報ψｄは、次の
ヒストグラム作成回路４６に入力される。ここで、ヒス
トグラム作成回路４６は、アキュムレータ４８と、移動
平均回路５０と、で構成される。ここで、アキュムレー
タ４８は、各ドプラ偏移周波数毎に設けられた複数のカ
ウンタで構成される。
【００２６】図２には、ヒストグラム（頻度グラフ）作
成原理が示されている。横軸は微小ドプラ偏移周波数Δ
ψｄ毎に区分され、縦軸にその度数が計数されている。
すなわち、ヒストグラム作成回路４６は、ドプラ偏移情
報が入力されると、その偏移周波数の大きさを判断し、
次に、その大きさに該当する区間のカウント値（度数
値）を１つ増加させる。そして、それを所定個のドプラ
偏移情報について順次行う。
成原理が示されている。横軸は微小ドプラ偏移周波数Δ
ψｄ毎に区分され、縦軸にその度数が計数されている。
すなわち、ヒストグラム作成回路４６は、ドプラ偏移情
報が入力されると、その偏移周波数の大きさを判断し、
次に、その大きさに該当する区間のカウント値（度数
値）を１つ増加させる。そして、それを所定個のドプラ
偏移情報について順次行う。
【００２７】このように形成されるヒストグラムにおい
て、Δψｄ毎の各区間の度数は、受信信号中のドプラ偏
移情報がその区間内に入る確率を示すものであるといえ
る。換言すれば、その度数（頻度）は、受信信号中のド
プラ偏移周波数ｆｄが、ある区間ｆｄｉとｆｄｉ＋Δｆ
ｄとの間に入る確率を示し、頻度分布は確率分布を意味
するものとなる。
て、Δψｄ毎の各区間の度数は、受信信号中のドプラ偏
移情報がその区間内に入る確率を示すものであるといえ
る。換言すれば、その度数（頻度）は、受信信号中のド
プラ偏移周波数ｆｄが、ある区間ｆｄｉとｆｄｉ＋Δｆ
ｄとの間に入る確率を示し、頻度分布は確率分布を意味
するものとなる。
【００２８】例えば、ヒストグラム作成に用いられた全
データ数をｎとし、ある区間であるｆｄｉとｆｄｉ＋Δ
ｆｄの間の度数をｍとすれば、ｍ／ｎは受信信号中のド
プラ偏移周波数ｆｄが当該区間内に入る確率を意味す
る。
データ数をｎとし、ある区間であるｆｄｉとｆｄｉ＋Δ
ｆｄの間の度数をｍとすれば、ｍ／ｎは受信信号中のド
プラ偏移周波数ｆｄが当該区間内に入る確率を意味す
る。
【００２９】従って、頻度分布は、パワースペクトラム
の確率関数Ｐ（ｆｄ）と同じ表現形式となり、頻度分布
であるヒストグラムはパワースペクトラムに相当するも
のとして取り扱うことができ、従来のＦＦＴによるパワ
ースペクトル表示に代えて、本実施例の頻度グラフをパ
ワースペクトル又は擬似パワースペクトルとして表示す
ることができる。
の確率関数Ｐ（ｆｄ）と同じ表現形式となり、頻度分布
であるヒストグラムはパワースペクトラムに相当するも
のとして取り扱うことができ、従来のＦＦＴによるパワ
ースペクトル表示に代えて、本実施例の頻度グラフをパ
ワースペクトル又は擬似パワースペクトルとして表示す
ることができる。
【００３０】ただし、各受信信号の振幅値が一定とみな
せる場合に上述の仮定が成り立ち、一定とみなせない場
合には、受信信号の振幅値に応じた重み付けを行う必要
があるが、それについては後述する。
せる場合に上述の仮定が成り立ち、一定とみなせない場
合には、受信信号の振幅値に応じた重み付けを行う必要
があるが、それについては後述する。
【００３１】さて、本実施例の装置では、２つの計測モ
ードが採用されている。１つは単サンプル点モードであ
り、他の１つは複数サンプル点モードである。
ードが採用されている。１つは単サンプル点モードであ
り、他の１つは複数サンプル点モードである。
【００３２】単サンプル点モードが選択された場合、図
６に示すように、生体内に１つサンプル点Ｓのみが設定
され、そのサンプル点Ｓにつき複数回（例えば、５００
回）の超音波パルスの送受波が行われ、それにより得ら
れた同一サンプル点についての複数個のドプラ情報によ
りヒストグラムが作成される。
６に示すように、生体内に１つサンプル点Ｓのみが設定
され、そのサンプル点Ｓにつき複数回（例えば、５００
回）の超音波パルスの送受波が行われ、それにより得ら
れた同一サンプル点についての複数個のドプラ情報によ
りヒストグラムが作成される。
【００３３】なお、その単サンプル点モードの場合、サ
ンプルホールド回路３０のゲート期間は、当該サンプル
点の大きさに合わせて設定され、ＰＲＦ毎に１回のゲー
ト抽出動作が実行される。これと同様に、Ａ／Ｄ変換回
路３６でもＰＲＦ毎に１回Ａ／Ｄ変換が行われる。
ンプルホールド回路３０のゲート期間は、当該サンプル
点の大きさに合わせて設定され、ＰＲＦ毎に１回のゲー
ト抽出動作が実行される。これと同様に、Ａ／Ｄ変換回
路３６でもＰＲＦ毎に１回Ａ／Ｄ変換が行われる。
【００３４】複数サンプル点モードが選択された場合、
図６に示すように、超音波ビーム上に小領域が設定さ
れ、その小領域内に、互いに隣接する複数の（例えば２
０個の）サンプル点から成るサンプル点群Ｇが設定され
る。すなわち、サンプルホールド回路３０のゲート期間
は、当該サンプル点の大きさに合わせて設定され、ゲー
ト抽出動作はサンプル点の数に合わせて、所定の回数、
連続して実行される。すなわち、ＰＲＦ毎に複数回のゲ
ート抽出動作が実行される。これと同様に、Ａ／Ｄ変換
回路３６でもＰＲＦ毎に複数回のＡ／Ｄ変換が行われ
る。
図６に示すように、超音波ビーム上に小領域が設定さ
れ、その小領域内に、互いに隣接する複数の（例えば２
０個の）サンプル点から成るサンプル点群Ｇが設定され
る。すなわち、サンプルホールド回路３０のゲート期間
は、当該サンプル点の大きさに合わせて設定され、ゲー
ト抽出動作はサンプル点の数に合わせて、所定の回数、
連続して実行される。すなわち、ＰＲＦ毎に複数回のゲ
ート抽出動作が実行される。これと同様に、Ａ／Ｄ変換
回路３６でもＰＲＦ毎に複数回のＡ／Ｄ変換が行われ
る。
【００３５】上述した単サンプル点モードが選択された
場合、複数回の送受波により形成されるヒストグラムを
そのまま表示すると、滑らかにその変動を表示できな
い。そこで、本実施例の装置では、移動平均回路５０が
設けられており、例えば１０ｍｓに設定された時間窓を
スイープさせて、各時間窓内でヒストグラムが形成され
るように、移動平均化処理を行っている。
場合、複数回の送受波により形成されるヒストグラムを
そのまま表示すると、滑らかにその変動を表示できな
い。そこで、本実施例の装置では、移動平均回路５０が
設けられており、例えば１０ｍｓに設定された時間窓を
スイープさせて、各時間窓内でヒストグラムが形成され
るように、移動平均化処理を行っている。
【００３６】以上のように、ヒストグラムが作成された
後、Ｄ／Ａ変換器５２を介して、表示器５４に、本実施
例の頻度グラフがパワースペクトル又は擬似パワースペ
クトルとして表示される。その表示例が図３に示されて
いる。ここで、度数（パワーに相当）は輝度として表さ
れている。
後、Ｄ／Ａ変換器５２を介して、表示器５４に、本実施
例の頻度グラフがパワースペクトル又は擬似パワースペ
クトルとして表示される。その表示例が図３に示されて
いる。ここで、度数（パワーに相当）は輝度として表さ
れている。
【００３７】次に、図４及び図５を用いて、第２実施例
について説明する。
について説明する。
【００３８】上述のように複数個のドプラ偏移情報を用
いてヒストグラムを作成する場合において、それらのド
プラ偏移情報を含んでいた各受信信号の振幅値が一定で
あるとみなせなければ、ヒストグラムがパワースペクト
ラムに相当するとの仮定が成り立たなくなる。
いてヒストグラムを作成する場合において、それらのド
プラ偏移情報を含んでいた各受信信号の振幅値が一定で
あるとみなせなければ、ヒストグラムがパワースペクト
ラムに相当するとの仮定が成り立たなくなる。
【００３９】そこで、図４に示す第２実施例では、一対
のＡ／Ｄ変換器３６の出力信号から受信信号の振幅値を
求める振幅値演算回路５６と、振幅値に基づきヒストグ
ラムに対して重み付けを行う重み付け回路５８と、が設
けられている。
のＡ／Ｄ変換器３６の出力信号から受信信号の振幅値を
求める振幅値演算回路５６と、振幅値に基づきヒストグ
ラムに対して重み付けを行う重み付け回路５８と、が設
けられている。
【００４０】すなわち、アキュムレータ４８では、図２
に示した各区間毎に度数カウントを行っているが、重み
付け回路５８はそのカウント値に重み付けを行ってい
る。例えば、振幅値自体を加算してヒストグラムを形成
してもよい。その重み付けされたヒストグラムが図５に
示されている。ここで、縦軸は振幅値を累積加算したも
のである。
に示した各区間毎に度数カウントを行っているが、重み
付け回路５８はそのカウント値に重み付けを行ってい
る。例えば、振幅値自体を加算してヒストグラムを形成
してもよい。その重み付けされたヒストグラムが図５に
示されている。ここで、縦軸は振幅値を累積加算したも
のである。
【００４１】このような重み付けによれば、ドプラ偏移
情報をヒストグラム値（度数）に反映させる場合に、受
信信号のレベルに応じた値が加算されるので、最終的に
作成されるヒストグラムに各受信信号のパワーの大きさ
を反映させることができる。よって、そのように作成さ
れたヒストグラムは、パワースペクトラムに相当するも
のであるといえる。なお、この図４に示した第２実施例
は、上述の単サンプル点モード及び複数サンプル点モー
ドのいずれのモードにも適用させることができる。ま
た、図３に示したような表示を行うことができる。
情報をヒストグラム値（度数）に反映させる場合に、受
信信号のレベルに応じた値が加算されるので、最終的に
作成されるヒストグラムに各受信信号のパワーの大きさ
を反映させることができる。よって、そのように作成さ
れたヒストグラムは、パワースペクトラムに相当するも
のであるといえる。なお、この図４に示した第２実施例
は、上述の単サンプル点モード及び複数サンプル点モー
ドのいずれのモードにも適用させることができる。ま
た、図３に示したような表示を行うことができる。
【００４２】単サンプル点モードでは、計測に時間的な
幅が生じてしまうものの、計測対象の空間的な広がりを
最小限にしつつ、注目する点についてのヒストグラムを
作成できる。一方、複数サンプル点モードでは、超音波
ビームに沿った空間的な計測対象の広がりが生じるもの
の、より少ない送受波で空間的速度分布を作成できると
いう利点がある。
幅が生じてしまうものの、計測対象の空間的な広がりを
最小限にしつつ、注目する点についてのヒストグラムを
作成できる。一方、複数サンプル点モードでは、超音波
ビームに沿った空間的な計測対象の広がりが生じるもの
の、より少ない送受波で空間的速度分布を作成できると
いう利点がある。
【００４３】なお、パワースペクトラムと共に二次元ド
プラ画像を表示する場合には従来同様に自己相関回路を
別途設ければよい。広義の意味の位相差分演算には、自
己相関演算が含まれるが、本実施例の位相差分演算回路
は１つのＲＡＭと１つの加算器（減算器）で構成でき、
極めて構成が簡易化されている。
プラ画像を表示する場合には従来同様に自己相関回路を
別途設ければよい。広義の意味の位相差分演算には、自
己相関演算が含まれるが、本実施例の位相差分演算回路
は１つのＲＡＭと１つの加算器（減算器）で構成でき、
極めて構成が簡易化されている。
【００４４】
【発明の効果】以上説明したように、本発明によれば、
パワースペクトラムに相当するヒストグラムを簡易な構
成で迅速に作成できるという効果がある。
パワースペクトラムに相当するヒストグラムを簡易な構
成で迅速に作成できるという効果がある。
【図１】本発明に係る超音波ドプラ診断装置の全体構成
を示すブロック図である。
を示すブロック図である。
【図２】ドプラ偏移情報によって作成されるヒストグラ
ムを示す説明図である。
ムを示す説明図である。
【図３】パワースペクトルの表示例を示す図である。
【図４】重み付けヒストグラムを作成する構成を示す図
である。
である。
【図５】重み付けヒストグラムを示す図である。
【図６】単一のサンプル点とサンプル点群を示す図であ
る。
る。
【図７】パワースペクトラムを示す図である。
【図８】パワースペクトラムの表示例を示す図である。
【図９】従来例を示す図である。
３８ 逆正接回路 ４０ 位相差分演算回路 ４２ ディレーライン ４４ 加算器（減算器） ４６ ヒストグラム作成回路 ４８ アキュムレータ ５０ 移動平均回路 ５６ 振幅値演算回路 ５８ 重み付け回路
Claims (4)
- 【請求項１】 検波後の受信信号間で位相差分演算を行
い、ドプラ偏移情報を求める位相差分演算回路と、 複数個の前記ドプラ偏移情報に基づき、ドプラ偏移の大
きさ毎の度数をグラフ化したヒストグラムを作成するヒ
ストグラム作成回路と、 を含むことを特徴とする超音波ドプラ診断装置。 - 【請求項２】 請求項１記載の装置において、 前記複数個のドプラ偏移情報は、同一サンプル点に対す
る複数回の送受波により得られることを特徴とする超音
波ドプラ診断装置。 - 【請求項３】 請求項１記載の装置において、 前記複数個のドプラ偏移情報は、同一超音波ビーム上に
隣接設定されたサンプル点群に対する送受波毎に得られ
ることを特徴とする超音波ドプラ診断装置。 - 【請求項４】 請求項２又は３記載の装置において、 前記受信信号の振幅値を演算する振幅値演算回路を含
み、 前記ヒストグラム作成回路は、振幅値の大きさに基づき
重み付けを行いつつ前記ヒストグラムを作成することを
特徴とする超音波ドプラ診断装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP3999794A JP3352211B2 (ja) | 1994-03-10 | 1994-03-10 | 超音波ドプラ診断装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP3999794A JP3352211B2 (ja) | 1994-03-10 | 1994-03-10 | 超音波ドプラ診断装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH07246206A JPH07246206A (ja) | 1995-09-26 |
| JP3352211B2 true JP3352211B2 (ja) | 2002-12-03 |
Family
ID=12568577
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP3999794A Expired - Fee Related JP3352211B2 (ja) | 1994-03-10 | 1994-03-10 | 超音波ドプラ診断装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3352211B2 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20240078175A (ko) * | 2022-11-25 | 2024-06-03 | 주식회사 알앤에스랩 | 혈관 모니터링을 위한 초음파 센서 및 그 제어방법 |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH11309146A (ja) * | 1998-03-26 | 1999-11-09 | General Electric Co <Ge> | 超音波散乱体の流れをイメ―ジングするシステムおよび方法 |
| JP4711498B2 (ja) * | 2000-10-12 | 2011-06-29 | 株式会社日立メディコ | 超音波画像診断装置 |
| JP5325617B2 (ja) * | 2009-03-13 | 2013-10-23 | 富士フイルム株式会社 | 超音波プローブ及び超音波診断装置 |
-
1994
- 1994-03-10 JP JP3999794A patent/JP3352211B2/ja not_active Expired - Fee Related
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20240078175A (ko) * | 2022-11-25 | 2024-06-03 | 주식회사 알앤에스랩 | 혈관 모니터링을 위한 초음파 센서 및 그 제어방법 |
| KR102811568B1 (ko) | 2022-11-25 | 2025-05-26 | (주)알앤에스랩 | 혈관 모니터링을 위한 초음파 센서 및 그 제어방법 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH07246206A (ja) | 1995-09-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4722283B2 (ja) | 連続データ獲得を用いた超音波フロー・イメージングにおける運動の可視化のための方法および装置 | |
| US6095980A (en) | Pulse inversion doppler ultrasonic diagnostic imaging | |
| CN100387199C (zh) | 超声波诊断装置和超声波诊断装置的控制方法 | |
| JP2849827B2 (ja) | パルスドップラ信号のフーリエプロセッシング方法および装置 | |
| US5971927A (en) | Ultrasonic diagnostic apparatus for obtaining blood data | |
| EP3513735B1 (en) | Device and method for generating ultrasound vector doppler image using plane wave synthesis | |
| JP4627366B2 (ja) | パケット・データ獲得を用いた超音波フロー撮像における運動の可視化のための方法および装置 | |
| JPH0120899B2 (ja) | ||
| US7666142B2 (en) | Ultrasound doppler diagnostic apparatus and image data generating method | |
| JPH10165400A (ja) | 超音波診断装置 | |
| JP2004329609A (ja) | 超音波診断装置 | |
| JP3352211B2 (ja) | 超音波ドプラ診断装置 | |
| JP2823252B2 (ja) | 超音波診断装置 | |
| EP1627602B1 (en) | Ultrasonographic device and ultrasonographic device data processing method | |
| JPH0218094B2 (ja) | ||
| EP1006369B1 (en) | Method and apparatus for adaptive filtering by counting acoustic sample zeroes in ultrasound imaging | |
| JP2916219B2 (ja) | 超音波診断装置 | |
| JP2022154977A (ja) | 超音波診断装置 | |
| JP2563656B2 (ja) | 超音波ドプラ映像装置 | |
| JPH08280689A (ja) | 超音波ドプラ診断装置 | |
| JP3391578B2 (ja) | 相関装置および流れ情報表示装置 | |
| JP2703943B2 (ja) | 超音波診断装置 | |
| JPH07286999A (ja) | 超音波診断装置 | |
| JPH064073B2 (ja) | 超音波血流イメージング装置 | |
| JPH03277351A (ja) | 超音波診断装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |