(57)【要約】【課題】  超音波ビームの入射角度に起因する血流計測誤差を高精度で補正する。【解決手段】  当該被検体の画像データに表示された血管の関心部位から得られるドプラ信号に基づいてドプラスペクトラムデータを生成する際、第１の画像データの関心部位にサンプルゲート設定部９２が設定したサンプルゲートを含み第１の画像データと交叉する第２の画像データを生成する。次いで、方向ベクトル設定部９３は、サンプルゲートを起点として設定した血管走行の方向ベクトルの前記画像データの各々に対する投影成分がこれらの画像データにおける血管の走行方向と一致するように方向ベクトルを更新する。そして、入射角度計測部５１は、更新後の方向ベクトルに基づいて関心部位に対する超音波ビームの入射角度を計測し、スペクトラムデータ補正部５２は、この入射角度を用いて前記ドプラスペクトラムデータを補正する。【選択図】  図１
2008073422.tif 000002
【特許請求の範囲】
【請求項１】
被検体に対して超音波の送受信を行なう振動素子を備えた超音波プローブと、
前記振動素子に対して送受信を行なう送受信手段と、
この送受信手段によって得られた受信信号を処理し交叉する複数画像断面の各々に対応する画像データを生成する画像データ生成手段と、
前記受信信号におけるドプラ成分を検出してドプラスペクトラムデータを生成するスペクトラムデータ生成手段と、
前記画像データの関心部位に設定されたサンプルゲートを起点として血流方向あるいは反射体の移動方向を示す方向ベクトルを３次元的に設定する方向ベクトル設定手段と、
前記方向ベクトルの設定情報に基づいて前記血流方向あるいは反射体の移動方向に対する超音波の入射角度を計測する入射角度計測手段と、
前記入射角度に基づいて前記ドプラスペクトラムデータを補正するスペクトラムデータ補正手段とを
備えたことを特徴とする超音波診断装置。
【請求項２】
前記超音波プローブは、２次元配列された複数個の振動素子を備えたことを特徴とする請求項１記載の超音波診断装置。
【請求項３】
前記画像データ生成手段は、前記受信信号を処理して２次元Ｂモード画像データ及び２次元カラードプラ画像データの少なくとも何れかを生成することを特徴とする請求項１記載の超音波診断装置。
【請求項４】
前記スペクトラムデータ生成手段は、前記受信信号の直交位相検波によって検出した前記ドプラ成分をＦＦＴ分析し、得られた周波数スペクトラムを時系列的に配列して前記ドプラスペクトラムデータを生成することを特徴とする請求項１記載の超音波診断装置。
【請求項５】
前記複数画像断面を設定する画像断面設定手段を備え、前記画像断面設定手段は、所定の基準画像断面（第１の画像断面）において得られる前記画像データの前記関心部位に設定されたサンプルゲートを含み前記基準画像断面と交叉する１つあるいは複数の第２の画像断面を設定することを特徴とする請求項１記載の超音波診断装置。
【請求項６】
前記複数画像断面に対する前記方向ベクトルの投影成分を算出する投影成分演算手段を備え、前記方向ベクトル設定手段は、前記複数画像断面の各々において生成された画像データに重畳表示される前記方向ベクトルの投影成分に基づき、前記方向ベクトルを前記血流方向あるいは反射体の移動方向へ設定することを特徴とする請求項１記載の超音波診断装置。
【請求項７】
前記方向ベクトル設定手段は、前記複数画像断面に対する方向ベクトルの投影成分と前記サンプルゲートを含み前記複数画像断面の各々に対し所定角度で交叉した平面に対する前記方向ベクトルの投影成分に基づき、前記方向ベクトルを前記血流方向あるいは反射体の移動方向へ設定することを特徴とする請求項１記載の超音波診断装置。
【請求項８】
前記スペクトラムデータ補正手段は、前記超音波の入射角度に基づいて前記ドプラスペクトラムデータにおける血流速度あるいはドプラ周波数を補正することを特徴とする請求項１記載の超音波診断装置。
【請求項９】
表示手段を備え、前記表示手段は、前記画像データ生成手段が生成した複数画像断面の画像データをリアルタイム表示することを特徴とする請求項１記載の超音波診断装置。
【請求項１０】
表示手段を備え、前記表示手段は、前記複数画像断面の各々において得られた複数画像断面の画像データの少なくとも何れかと、前記スペクトラムデータ補正手段によって補正された前記ドプラスペクトラムデータを並列表示することを特徴とする請求項１記載の超音波診断装置。
【請求項１１】
表示手段を備え、前記表示手段は、前記画像データ生成手段が生成した複数画像断面の各々における画像データと方向ベクトルを２次元的あるいは３次元的に合成して表示することを特徴とする請求項１記載の超音波診断装置。
【請求項１２】
サンプルゲート設定手段が、基準画像断面（第１の画像断面）において生成された第１の画像データの関心部位に対してサンプルゲートを設定するステップと、
画像断面設定手段が、前記サンプルゲートを含み前記基準画像断面と交叉する１つあるいは複数の第２の画像断面を設定するステップと、
方向ベクトル設定手段が、前記サンプルゲートを起点とする方向ベクトルを設定するステップと、
投影成分演算手段が、前記第２の画像断面において生成された第２の画像データあるいは前記第１の画像データと前記第２の画像データに対する前記方向ベクトルの投影成分を算出するステップと、
前記方向ベクトル設定手段が、前記第２の画像データあるいは前記第１の画像データと前記第２の画像データに対して重畳表示された前記投影成分に基づいて前記方向ベクトルの方向を更新するステップと、
入射角度計測手段が、更新された前記方向ベクトルの情報に基づいて前記関心部位に対する超音波の入射角度を計測するステップと、
スペクトラムデータ生成手段が、前記関心部位からの受信超音波に基づく受信信号を処理してドプラスペクトラムデータを生成するステップと、
スペクトラムデータ補正手段が、前記ドプラスペクトラムデータを前記超音波の入射角度に基づいて補正するステップとを
有することを特徴とする超音波計測方法
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波のドプラ効果を利用して、生体内の血流の流速情報や組織の移動情報などの計測を行なう超音波診断装置及び超音波計測方法に関する。
【背景技術】
【０００２】
超音波診断装置は、超音波プローブに内蔵された振動素子から発生する超音波パルス（送信超音波）を被検体内に放射し、被検体組織の音響インピーダンスの差異によって生ずる反射波（受信超音波）を上記振動素子により電気信号に変換してモニタ上に表示するものである。この診断方法は、超音波プローブを体表に接触させるだけの簡単な操作でリアルタイムの２次元画像が容易に観察できるため、各種臓器の機能診断や形態診断に広く用いられている。生体内の組織あるいは血球からの反射波により生体情報を得る超音波診断法は、超音波パルス反射法と超音波ドプラ法の２つの大きな技術開発により急速な進歩を遂げ、上記技術を用いて得られるＢモード画像とカラードプラ画像は、今日の超音波画像診断において不可欠のものとなっている。
【０００３】
一方、被検体の関心部位における血流速度を定量的且つ正確に得る方法としてドプラスペクトラム法がある。このドプラスペクトラム法では、前記関心部位に対し所定時間間隔で複数回の超音波送受信を行ない、血球などの移動反射体において反射した受信超音波に対し送信超音波の中心周波数と略等しい周波数を有する基準信号を用いた直交位相検波を行なってドプラ信号を検出する。そして、このドプラ信号の中から上述の関心部位におけるドプラ信号をサンプルゲートによって抽出し、抽出した所定区間のドプラ信号をＦＦＴ分析することにより周波数スペクトラムを計測する。
【０００４】
更に、当該被検体の関心部位から得られる一連のドプラ信号に対し周波数スペクトラムの計測を所定時間間隔で繰り返し、得られた複数の周波数スペクトラムを時系列的に配列することによって、所謂ドプラスペクトラムデータを生成する。周波数スペクトラムの時間的変化を示すドプラスペクトラムデータは、通常、縦軸にドプラ周波数あるいは血流速度、横軸に時間が設定され、各周波数成分の大きさ（パワー）は輝度によって表現される。尚、サンプルゲートが被検体の関心部位に対して正確に設定されていることを確認するために、サンプルゲートの設定はＢモード画像データあるいはカラードプラ画像データの観測下にて行なわれ、このとき、サンプルゲートの位置はこれらの画像データにおいてサンプルゲートマーカとして表示される。
【０００５】
ところで、被検体の関心部位から得られる反射超音波のドプラ偏移周波数に基づいて血流速度を計測する場合、実際に計測される血流速度は、超音波の送受信方向に沿った速度成分であり、送受信方向と血流方向が一致しない場合には真の血流速度を捉えることはできない。
【０００６】
そこで、被検体に対する２次元超音波走査により２次元のＢモード画像データやカラードプラ画像データの生成を可能とする超音波診断装置では、２次元画像データ上の関心部位における血管の走行方向と超音波送信方向との交叉角度（以下では、超音波ビームの入射角度と呼ぶ。）を計測し、ＦＦＴ分析によって得られた血流速度を上述の入射角度によって補正する方法が行なわれている（例えば、特許文献１参照。）。
【０００７】
この場合、先ず、２次元画像データ上に表示された血管の関心部位に対するサンプルゲートの設定と、このサンプルゲートを起点とした血管の走行方向を示すベクトル（以下では、方向ベクトルと呼ぶ。）の設定が行なわれる。次いで、前記サンプルゲートを通過する超音波送信方向と血管走行の方向ベクトルとの交叉角度（即ち、血管走行に対する超音波ビームの入射角度）が算出され、この超音波ビームの入射角度に基づいて血流速度の補正が行なわれる。
【０００８】
しかしながら、上述の方法によって得られた血管走行の方向ベクトルは２次元画像データの断面内において設定されるものであり、必ずしも真の方向ベクトルを示していない。即ち、真の方向ベクトルが２次元画像データの断面内に存在しない場合には相変わらず血流速度を正確に計測することはできない。
【０００９】
一方、近年では、振動素子が２次元配列された超音波プローブを用いることによって生体内の３次元データ（ボリュームデータ）を収集することが可能となり、このボリュームデータを用い超音波ビームの入射角度に起因する誤差を補正する方法が提案されている（例えば、特許文献２参照。）。
【００１０】
ボリュームデータの収集が可能な超音波診断装置を用いて血流速度の計測を行なう場合、関心部位の縦断面情報が含まれる２次元データを前記ボリュームデータの中から抽出することにより好適な２次元画像データを生成することができる。そして、この２次元画像データに対して設定された方向ベクトル（真の方向ベクトル）に基づいて計測された超音波ビームの入射角度を用いることによりドプラスペクトラムデータを高精度で補正することができる。
【特許文献１】特開平０４−１０８４３５号公報
【特許文献２】特開２００５−９５２７８号公報
【発明の開示】
【発明が解決しようとする課題】
【００１１】
しかしながら、上述の特許文献２に記載された方法によれば、ボリュームデータの中から抽出した２次元データに基づいて方向ベクトルが設定されているため、ボリュームデータの収集や２次元データの抽出に多くの時間を要する。このため、拍動性あるいは呼吸性の移動を伴なう臓器の血流計測に上述の方法を適用することは困難である。
【００１２】
本発明は、上述の問題点に鑑みてなされたものであり、その目的は、超音波のドプラ効果を用いて生体組織の移動速度や血流速度を計測する際に、超音波ビームの入射角度に起因する計測誤差を短時間かつ高精度で補正することが可能な超音波診断装置及び超音波計測方法を提供することにある。
【課題を解決するための手段】
【００１３】
上記課題を解決するために、請求項１に係る本発明の超音波診断装置は、被検体に対して超音波の送受信を行なう振動素子を備えた超音波プローブと、前記振動素子に対して送受信を行なう送受信手段と、この送受信手段によって得られた受信信号を処理し交叉する複数画像断面の各々に対応する画像データを生成する画像データ生成手段と、前記受信信号におけるドプラ成分を検出してドプラスペクトラムデータを生成するスペクトラムデータ生成手段と、前記画像データの関心部位に設定されたサンプルゲートを起点として血流方向あるいは反射体の移動方向を示す方向ベクトルを３次元的に設定する方向ベクトル設定手段と、前記方向ベクトルの設定情報に基づいて前記血流方向あるいは反射体の移動方向に対する超音波の入射角度を計測する入射角度計測手段と、前記入射角度に基づいて前記ドプラスペクトラムデータを補正するスペクトラムデータ補正手段とを備えたことを特徴としている。
【００１４】
又、請求項１２に係る本発明の超音波計測方法は、サンプルゲート設定手段が、基準画像断面（第１の画像断面）において生成された第１の画像データの関心部位に対してサンプルゲートを設定するステップと、画像断面設定手段が、前記サンプルゲートを含み前記基準画像断面と交叉する１つあるいは複数の第２の画像断面を設定するステップと、方向ベクトル設定手段が、前記サンプルゲートを起点とする方向ベクトルを設定するステップと、投影成分演算手段が、前記第２の画像断面において生成された第２の画像データあるいは前記第１の画像データと前記第２の画像データに対する前記方向ベクトルの投影成分を算出するステップと、前記方向ベクトル設定手段が、前記第２の画像データあるいは前記第１の画像データと前記第２の画像データに対して重畳表示された前記投影成分に基づいて前記方向ベクトルの方向を更新するステップと、入射角度計測手段が、更新された前記方向ベクトルの情報に基づいて前記関心部位に対する超音波の入射角度を計測するステップと、スペクトラムデータ生成手段が、前記関心部位からの反射超音波に基づく受信信号を処理してドプラスペクトラムデータを生成するステップと、スペクトラムデータ補正手段が、前記ドプラスペクトラムデータを前記超音波の入射角度に基づいて補正するステップとを有することを特徴としている。
【発明の効果】
【００１５】
本発明によれば、超音波のドプラ効果を用いて生体組織の移動速度や血流速度を計測する際に、移動方向あるいは血流方向の３次元的な設定を容易に行なうことができるため超音波ビームの入射角度に起因する計測誤差を短時間かつ高精度で補正することが可能となる。
【発明を実施するための最良の形態】
【００１６】
以下、図面を参照して本発明の実施例を説明する。
【実施例１】
【００１７】
以下に述べる本発明の実施例では、当該被検体に対し複数枚の画像データを生成すると共に血流計測が要求される血管の関心部位から得られるドプラ信号に基づいてドプラスペクトラムデータを生成する際、第１の画像データにおける血管の関心部位に対しサンプルゲートを設定し、このサンプルゲートを含み第１の画像データと交叉する第２の画像データを生成する。次いで、血管走行の方向ベクトルを、サンプルゲートを起点として設定し、前記画像データの各々に対する方向ベクトルの投影成分がこれらの画像データにおける血管の走行方向と一致するように方向ベクトルの方向を更新する。そして、更新された方向ベクトルに基づいて関心部位に対する超音波ビームの入射角度を計測し、この入射角度を用いてドプラスペクトラムデータにおける血流速度あるいはドプラ周波数を補正する。
【００１８】
尚、以下では、振動素子が２次元配列された超音波プローブを有するセクタ走査方式の超音波診断装置を用いて交叉する複数枚のＢモード画像データを生成し、これらのＢモード画像データに基づいて方向ベクトルの設定と超音波ビーム入射角度の計測を行なう場合について述べるがこれに限定されるものではない。
【００１９】
（装置の構成）
以下では、本発明の実施例における超音波診断装置の構成と各ユニットの基本動作につき図１乃至図５を用いて説明する。尚、図１は、本実施例における超音波診断装置の全体構成を示すブロック図であり、図２は、この超音波診断装置を構成する送受信部及び超音波データ生成部のブロック図である。
【００２０】
図１に示す超音波診断装置１００は、交叉する複数画像断面の各々に沿って超音波の送受信を行なう複数個の振動素子を有した超音波プローブ１と、これら振動素子の各々に対し駆動パルスを供給すると共に前記振動素子から得られる複数チャンネルの受信信号を整相加算（所定方向からの受信超音波に基づく受信信号を位相合わせして加算合成）する送受信部２と、送受信部２から出力された整相加算後の受信信号を処理してＢモードデータやドプラスペクトラムデータ等の超音波データを生成する超音波データ生成部３と、超音波データ生成部３において生成されたＢモードデータを順次保存して複数画像断面におけるＢモード画像データを生成する画像データ生成部４と、超音波データ生成部３において生成されたドプラスペクトラムデータの血流速度に対し超音波ビームの入射角度に基づいた補正を行なう速度補正部５を備えている。
【００２１】
又、超音波診断装置１００は、後述の入力部９にて設定あるいは更新された方向ベクトルの複数画像断面への投影成分を算出する投影成分演算部６と、この投影成分が重畳された複数画像断面におけるＢモード画像データや速度補正部５においてその血流速度が補正されたドプラスペクトラムデータを表示する表示部７と、送信超音波の中心周波数と略等しい周波数を有する連続波あるいは矩形波を発生し送受信部２及び超音波データ生成部３へ供給する基準信号発生部８と、画像断面の設定、サンプルゲートの設定、方向ベクトルの設定及び更新、各種コマンド信号の入力等を行なう入力部９と、上述の各ユニットを統括的に制御するシステム制御部１０を備えている。
【００２２】
超音波プローブ１は、被検体の表面に対してその前面を接触させ超音波の送受信を行なうものであり、２次元配列された複数個（Ｎ個）の微小な振動素子をその先端部に有している。この振動素子は電気音響変換素子であり、送信時には電気パルスを超音波パルス（送信超音波）に変換し、又、受信時には超音波反射波（受信超音波）を電気信号（受信信号）に変換する機能を有している。そして、これらの振動素子は、Ｎチャンネルのケーブルを介して送受信部２の送信部２１及び受信部２２に接続されている。尚、超音波プローブ１にはセクタ走査対応、リニア走査対応、コンベックス走査対応等があり、診断部位に応じて任意に選択される。以下では、セクタ走査対応の超音波プローブ１を用いた場合について述べるがこれに限定されない。
【００２３】
次に、送受信部２及び超音波データ生成部３の具体的な構成につき図２を用いて説明する。図２に示した送受信部２は、被検体の所定方向に対し送信超音波を放射するための駆動パルスを超音波プローブ１の振動素子へ供給する送信部２１と、これらの振動素子から得られたＮチャンネルの受信信号に対し整相加算を行なう受信部２２を備え、送信部２１は、レートパルス発生器２１１と、送信遅延回路２１２と、駆動回路２１３を備えている。レートパルス発生器２１１は、被検体内に放射する送信超音波の繰り返し周期を決定するレートパルスを、基準信号発生部８から供給される連続波あるいは矩形波を分周することによって生成し、このレートパルスを送信遅延回路２１２へ供給する。
【００２４】
送信遅延回路２１２は、例えば、送信に使用される振動素子と同数（Ｎチャンネル）の独立な遅延回路から構成され、送信において細いビーム幅を得るために所定の深さに送信超音波を収束するための遅延時間と所定方向に送信超音波を放射するための遅延時間を前記レートパルスに与える。そして、Ｎチャンネルから構成される駆動回路２１３は、超音波プローブ１に内蔵された振動素子を駆動するための駆動パルスを送信遅延回路２１２から供給される上述のレートパルスに基づいて生成する。
【００２５】
一方、受信部２２は、プリアンプ２２１と、受信遅延回路２２２と、加算器２２３を備えている。プリアンプ２２１は、振動素子により受信超音波から電気信号に変換された微小な受信信号を増幅して十分なＳ／Ｎを確保する。又、受信遅延回路２２２は、所定の深さからの受信超音波を収束するための遅延時間と、所定方向からの受信超音波に対して強い受信指向性を設定するための遅延時間をプリアンプ２２１から出力される受信信号に与える。そして、受信遅延回路２２２において所定の遅延時間が与えられた受信信号は加算器２２３へ送られ、この加算器２２３において加算合成（整相加算）される。
【００２６】
次に、超音波データ生成部３は、受信部２２の加算器２２３から出力される受信信号に対し所定の信号処理を行なってＢモードデータを生成するＢモードデータ生成部３１と、前記受信信号から検出したドプラ信号を周波数分析してドプラスペクトラムデータを生成するスペクトラムデータ生成部３２を備えている。
【００２７】
Ｂモードデータ生成部３１は、対数変換器３１１と包絡線検波器３１２とＡ／Ｄ変換器３１３を備えている。対数変換器３１１は、受信部２２の加算器２２３から供給された受信信号の振幅を対数変換して弱い信号を相対的に強調し、包絡線検波器３１２は、対数変換された受信信号を包絡線検波して振幅情報のみを検出する。そして、Ａ／Ｄ変換器３１３は、包絡線検波器３１２の出力信号をＡ／Ｄ変換してＢモードデータを生成する。
【００２８】
一方、スペクトラムデータ生成部３２は、π／２移相器３２１、ミキサ３２２−１及び３２２−２、ＬＰＦ（ローパスフィルタ）３２３−１及び３２３−２、ＳＨ（サンプルホールド回路）３２４−１及び３２４−２を備えており、更に、ＢＰＦ（バンドパスフィルタ）３２５−１及び３２５−２、Ａ／Ｄ変換器３２６−１及び３２６−２、ＦＦＴ（Fast-Fourier-Transform）分析器３２７を備えている。そして、受信部２２の加算器２２３から供給された受信信号を直交位相検波してドプラ信号を検出し、得られたドプラ信号に対してＦＦＴ分析を行なう。
【００２９】
次に、スペクトラムデータ生成部３２の構成と基本動作につき図３のタイムチャートを参照して更に詳しく説明する。受信部２２の加算器２２３から出力された受信信号（図３の（ｃ））は、スペクトラムデータ生成部３２のミキサ３２２−１及び３２２−２の第１の入力端子に入力される。一方、この受信信号の中心周波数と略等しい周波数をもった基準信号発生部８の基準信号（図３の（ａ））はミキサ３２２−１の第２の入力端子に直接供給され、π／２移相器３２１を介することにより位相が９０度だけシフトした基準信号はミキサ３２２−２の第２の入力端子に送られる。そして、ミキサ３２２−１及び３２２−２の出力は、ＬＰＦ３２３−１及び３２３−２に送られ、加算器２２３から出力される受信信号の周波数と基準信号発生部８から供給される基準信号の周波数との和の成分が除去され、差の成分のみがドプラ信号として検出される（図３の（ｄ））。
【００３０】
次に、ＳＨ３２４−１及び３２４−２には、ＬＰＦ３２３−１及び３２３−２から出力されたドプラ信号と、システム制御部１０が基準信号発生部８の基準信号を分周して生成したサンプリングパルス（サンプルゲートパルス）が供給され（図３の（ｅ））、このサンプルゲートパルスの発生タイミングに対応した距離に位置する関心部位からのドプラ信号のみがサンプルホールドされる（図３の（ｆ））。尚、このサンプルゲートパルスは、送信超音波の放射タイミングを示す周期Ｔｒのレートパルス（図３の（ｂ））から所定時間（Ｔｓ）後に発生し、この発生タイミングは入力部９の操作者によって任意に設定される。
【００３１】
即ち、操作者は、遅延時間Ｔｓを変更することによって超音波プローブ１から距離Ｌｇだけ離れた関心部位からのドプラ信号を選択的に検出することが可能となる。尚、超音波プローブ１から関心部位までの距離Ｌｇと遅延時間Ｔｓは、生体組織の音速度をＣとすれば、２Ｌｇ／Ｃ＝Ｔｓの関係にある。
【００３２】
次に、ＳＨ３２４−１及び３２４−２から出力された関心部位のドプラ信号に重畳する階段状のノイズ成分（図３の（ｆ））は、ＢＰＦ３２５−１及び３２５−２によって除去され（図３の（ｇ））、平滑化された上記ドプラ信号は、Ａ／Ｄ変換器３２６−１及び３２６−２によってデジタル信号に変換された後ＦＦＴ分析器３２７へ供給される。尚、上述のＢＰＦ３２５−１及び３２５−２は、血流速度より遅い移動速度を有する血管や生体組織等に起因したドプラ信号（クラッタ信号）を排除する機能をも併せもっている。
【００３３】
一方、ＦＦＴ分析器３２７は、図示しない演算回路と記憶回路を備え、前記記憶回路は、Ａ／Ｄ変換器３２６−１及び３２６−２から出力されたドプラ信号の保存を行ない、前記演算回路は、前記記憶回路に保存された一連のドプラ信号の所定区間においてＦＦＴ分析を行なって周波数スペクトラムを算出する。更に、ＦＦＴ分析器３２７は、ＦＦＴ分析を上述のドプラ信号に対して所定時間間隔で繰り返し行ない、このとき得られた複数の周波数スペクトラムを時系列的に配列してドプラスペクトラムデータを生成する。
【００３４】
図１に戻って、画像データ生成部４は記憶回路を備え、画像データの生成と保存を行なう機能を有している。即ち、被検体の交叉する複数画像断面に沿って送受信方向を順次変更しながら超音波送受信を繰り返すことによって得られた一連のＢモードデータを前記記憶回路に保存することにより、複数画像断面の各々に対して２次元のＢモード画像データを生成する。
【００３５】
次に、速度補正部５は、サンプルゲートを起点として設定された方向ベクトルの設定情報と前記サンプルゲートに対する超音波送信方向に関する情報に基づいて血管の走行方向（血流方向）に対する超音波ビームの入射角度を計測する入射角度計測部５１と、得られた入射角度を用いてドプラスペクトラムデータにおける血流速度あるいはドプラ周波数を補正するスペクトラムデータ補正部５２を備えている。尚、上述の方向ベクトルの設定情報は、入力部９からシステム制御部１０を介して供給され、サンプルゲートに対する超音波送信方向の情報は、システム制御部１０から供給される。
【００３６】
一方、投影成分演算部６は、血管走行を示す方向ベクトルの設定情報を入力部９から受信し、複数画像断面の各々に対する方向ベクトルの投影成分を算出する。
【００３７】
表示部７は、図示しない表示用データメモリと変換回路とモニタを備え、Ｂモード画像データやドプラスペクトラムデータ、更には、入力部９から入力される被検体情報等を示す文字や数字などは上記表示用データメモリで合成され、変換回路においてＤ／Ａ変換とテレビフォーマット変換が行われた後、ＣＲＴあるいは液晶などのモニタに表示される。
【００３８】
又、本実施例では、画像データ生成部４において生成された複数画像断面の各々におけるＢモード画像データと速度補正部５のスペクトラムデータ補正部５２において補正されたドプラスペクトラムデータの表示が表示部７において行なわれる。このとき、複数画像断面において生成された複数のＢモード画像データの各々には、投影成分演算部６によって算出された方向ベクトルの当該画像断面に対する投影成分が重畳表示される。
【００３９】
次に、入力部９は、操作パネル上に表示パネルやキーボード、トラックボール、マウス等の入力デバイスを備え、被検体に対し画像断面を設定する画像断面設定部９１と、前記画像断面において得られたＢモード画像データの関心部位にサンプルゲートを設定するサンプルゲート設定部９２と、サンプルゲートを起点として方向ベクトルを設定する方向ベクトル設定部９３を備えている。更に、入力部９では、被検体情報の入力、画像データの選択、超音波データ生成条件の設定、各種コマンド信号の入力等も上述の入力デバイスや表示パネルを用いて行なわれる。
【００４０】
システム制御部１０は、図示しないＣＰＵと記憶回路を備え、装置の初期設定時に入力部９にて行なわれる被検体情報の入力、画像データの選択、超音波データ生成条件の設定等に関する情報は前記記憶回路に保存され、前記ＣＰＵは、これらの初期設定情報に基づいて超音波診断装置１００が備える各ユニットの制御やシステム全体の制御を統括して行なう。
【００４１】
又、システム制御部１０は、入力部９の画像断面設定部９１から供給される複数画像断面の設定情報に基づいて、前記複数画像断面におけるＢモード画像データの生成を目的とした超音波送受信方向の制御を送受信部２の送信遅延回路２１２及び受信遅延回路２２２に対して行ない、更に、入力部９のサンプルゲート設定部９２から供給されるサンプルゲートの設定情報に基づいて、サンプルゲートに対する超音波送信方向の情報を速度補正部５の入射角度計測部５１へ供給する。
【００４２】
次に、速度補正部５の入射角度計測部５１において行なわれる超音波ビームの入射角度計測方法につき図４を用いて更に詳しく説明する。尚、ここでは、互いに直交する２つの画像断面において得られたＢモード画像データを用い超音波ビームの入射角度を計測する場合について述べるが、３つ以上の画像断面において得られた画像データを用いてもよい。
【００４３】
超音波ビームの入射角度計測方法を示す図４の図４（ａ）は、直交する第１のＢモード画像データＰ１及び第２のＢモード画像データＰ２と血管ＢＶの関心部位に設定されたサンプルゲートＳＧ及び方向ベクトルＤＶとの関係を模式的に示したものであり、図４（ｂ）は、第１のＢモード画像データＰ１及び第２のＢモード画像データＰ２の各々において表示される血管ＢＶの血管断面ＢＶＰ１及びＢＶＰ２と方向ベクトルＤＶの投影成分ＤＶＰ１及びＤＶＰ２を示している。尚、方向ベクトルＤＶの起点位置は、入力部９のサンプルゲート設定部９２が血管ＢＶに対して設定したサンプルゲートＳＧに基づいて設定され、その方向は入力部９の方向ベクトル設定部９３からの設定／更新情報に基づいて制御される。
【００４４】
例えば、Ｂモード画像データの観測下にて超音波プローブ１の位置や画像断面の位置を調整することにより、血流計測が要求される血管ＢＶの関心部位を含む基準画像断面（第１の画像断面）を設定する。次いで、入力部９のサンプルゲート設定部９２は、第１の画像断面において生成され表示部７に表示された第１のＢモード画像データＰ１の関心部位にサンプルゲートＳＧを設定し、入力部９の画像断面設定部９１は、サンプルゲートＳＧを含み第１の画像断面と直交する第２の画像断面を設定する。そして、第１の画像断面にて得られた第１のＢモード画像データＰ１と第２の画像断面にて得られた第２のＢモード画像データＰ２は表示部７において図４（ｂ）に示すように並列表示され、Ｂモード画像データの各々にはサンプルゲートＳＧを起点として任意の方向に設定された方向ベクトルＤＶの投影成分ＤＶＰ１及びＤＶＰ２が重畳表示される。
【００４５】
次いで、入力部９の方向ベクトル設定部９３は、Ｂモード画像データＰ１及びＰ２に表示された血管ＢＶの断面ＢＶＰ１及びＢＶＰ２に基づく血管の走行方向と、投影成分ＤＶＰ１及びＤＶＰ２の方向とが一致するように設定された方向ベクトルの方向を更新し、更新された方向ベクトルの情報を速度補正部５の入射角度計測部５１へ供給する。
【００４６】
そして、入射角度計測部５１は、方向ベクトル設定部９３から供給された方向ベクトルの設定情報とシステム制御部１０から供給されたサンプルゲートに対する超音波送信方向の情報に基づいて血管走行に対する超音波ビームの入射角度を計測する。
【００４７】
次に、速度補正部５のスペクトラムデータ補正部５２によって行なわれるドプラスペクトラムデータの補正方法について述べる。送信超音波の中心周波数をｆｃ、被検体組織の音速をＣ，血流速度をＶ、血流方向（即ち、血管の走行方向）に対する超音波ビームの入射角度をθとした場合、血球からの反射超音波が有するドプラ偏移周波数ｆｄは次式（１）によって表すことができる。
【数１】
2008073422.tif 000003
【００４８】
従がって、血流速度Ｖは、次式（２）によって算出することができる。
【数２】
2008073422.tif 000004
【００４９】
この場合、従来のような１枚の２次元画像データを用いて超音波ビームの入射角度を算出する方法は、既に述べた理由により大きな誤差が発生し正確な血流速度Ｖを計測することが困難である。これに対して、本実施例では３次元領域における真の超音波入射角度θを計測することができるため、正確な血流速度Ｖを得ることができる。
【００５０】
次に、本実施例におけるＢモード画像データ及びドプラスペクトラムデータの表示例につき図５を用いて説明する。
【００５１】
図５は、表示部７のモニタに表示されるＢモード画像データ及びドプラスペクトラムデータの具体例を示したものであり、上部領域には互いに直交する第１のＢモード画像データＰ１と第２のＢモード画像データＰ２が並列表示され、夫々の画像データには、図４において示した血管ＢＶの血管断面ＢＶＰ１及びＢＶＰ２や方向ベクトルＤＶの投影成分ＤＶＰ１及びＤＶＰ２、更には、サンプルゲートＳＧが表示されている。
【００５２】
一方、図５の下部領域には、周波数スペクトラムの時系列的な変化を示すドプラスペクトラムデータＤＳＤが表示され、このドプラスペクトラムデータＤＳＤにおける縦軸の目盛りは、補正後の血流速度あるいはドプラ周波数に基づいて更新される。
【００５３】
（ドプラスペクトラムデータの補正手順）
次に、本実施例におけるドプラスペクトラムデータの補正手順につき図６のフローチャートに沿って説明する。尚、ここでも説明を簡単にするために、互いに直交する２つの画像断面に対してＢモード画像データを生成する場合について述べるが、３つ以上の画像断面に対してＢモード画像データを生成してもよい。又、画像断面の交叉角度は任意に設定しても構わない。
【００５４】
ドプラスペクトラムデータの生成に先立って超音波診断装置１００の操作者は、入力部９にて被検体情報の入力、画像データの選択、超音波データ生成条件の設定等を行ない、これらの初期設定情報をシステム制御部１０の記憶回路に保存する。本実施例では、画像データとしてＢモード画像データの選択を行なう（図６のステップＳ１）。
【００５５】
上述の初期設定が終了したならば、操作者は、超音波プローブ１の先端部（超音波送受信面）を当該被検体の体表面に配置し、任意に設定した第１の画像断面における送受信方向θ１に対して超音波送受信を行なう。即ち、図２に示した送信部２１のレートパルス発生器２１１は、基準信号発生部８から供給された基準信号を分周して周期Ｔｒのレートパルスを生成し、このレートパルスを送信遅延回路２１２に供給する。
【００５６】
送信遅延回路２１２は、送信超音波を所定距離に収束するための遅延時間と、送受信方向θ１に超音波を送信するための遅延時間をレートパルスに与え、駆動回路２１３は、送信遅延回路２１２から出力されたレートパルスに基づいて生成した駆動パルスを超音波プローブ１の振動素子に供給して被検体内に送信超音波を放射する。
【００５７】
被検体内に放射された送信超音波の一部は、音響インピーダンスの異なる臓器の境界面あるいは組織にて反射し、又、送信超音波が心臓壁や血球など動きのある反射体で反射する場合は、その超音波周波数はドプラ偏移を受ける。被検体組織にて反射した受信超音波は、前記振動素子によって受信されて電気的な受信信号に変換され、この受信信号は、プリアンプ２２１にて増幅された後受信遅延回路２２２に送られる。
【００５８】
受信遅延回路２２２は、所定の深さからの受信超音波を収束するための遅延時間と、送受信方向θ１に対し強い受信指向性を持たせるための遅延時間をプリアンプ２２１から出力された受信信号に与え、この受信信号を加算器２２３へ供給する。そして、加算器２２３は、受信遅延回路２２２から出力された複数チャンネルの受信信号を整相加算し超音波データ生成部３へ供給する。
【００５９】
超音波データ生成部３のＢモードデータ生成部３１へ送られた整相加算後の受信信号は、対数変換、包絡線検波、Ａ／Ｄ変換がなされた後、図１に示した画像データ生成部４の記憶回路に保存される。
【００６０】
送受信方向θ１に対するＢモードデータの生成と保存が終了したならば、同様の手順により送受信方向θ２（θ２＝θ１＋Δθ）、θ３（θ３＝θ１＋２Δθ）・・・・に対し超音波の送受信を繰り返し行ない、このとき得られたＢモードデータも送受信方向に対応させて画像データ生成部４の記憶回路に保存される。即ち、画像データ生成部４の記憶回路において基準画像断面（第１の画像断面）における第１のＢモード画像データが生成され、この第１のＢモード画像データは、表示部７のモニタにリアルタイム表示される。
【００６１】
次いで、操作者は、表示部７に表示された第１のＢモード画像データを観測しながら暫定的に設定された第１の画像断面の方向を順次更新し、血管の関心部位の表示を可能とする方向に第１の画像断面を再設定する（図６のステップＳ２）。そして、再設定された第１の画像断面における第１のＢモード画像データを上述と同様な方法によって生成し、得られた第１のＢモード画像データを表示部７にリアルタイム表示する（図６のステップＳ３）。
【００６２】
次に、操作者は、入力部９のサンプルゲート設定部９２を操作することによって第１のＢモード画像データに表示された血管の関心部位にサンプルゲートを設定し（図６のステップＳ４）、更に、入力部９の画像断面設定部９１を操作することによって前記サンプルゲートを含み第１の画像断面に直交する第２の画像断面を設定する（図６のステップＳ５）。
【００６３】
一方、画像断面設定部９１から第２の画像断面の位置情報を受信したシステム制御部１０は、送信部２１の送信遅延回路２１２及び受信部２２の受信遅延回路２２２に対して遅延時間制御信号を供給し、第１の画像断面の場合と同様にして、第２の画像断面における送受信方向φ１、φ２（φ２＝φ１＋Δφ）、φ３（φ３＝φ１＋２Δφ）・・・に対し超音波送受信を行なう、そして、このとき画像データ生成部４において生成された第２のＢモード画像データは、表示部７において第１のＢモード画像データと共にリアルタイム表示される（図６のステップＳ６）。
【００６４】
表示部７にて同時表示された第１のＢモード画像データ及び第２のＢモード画像データを観察した操作者は、これらのＢモード画像データを観察しながら入力部９の方向ベクトル設定部９３を操作することによりＢモード画像データに示されたサンプルゲートを起点とする方向ベクトルを任意の方向に設定する。この場合、入力部９に備えられたトラックボールを操作することにより方向ベクトルの方向設定は短時間かつ容易に行なうことができる。
【００６５】
方向ベクトル設定部９３より方向ベクトルの設定情報を受信した投影成分演算部６は、前記方向ベクトルの第１の画像断面及び第２の画像断面に対する投影成分を算出する。そして、得られた投影成分を表示部７へ供給し、第１のＢモード画像データ及び第２のＢモード画像データに重畳して表示する。
【００６６】
次に、操作者は、入力部９の方向ベクトル設定部９３を操作することにより、第１のＢモード画像データ及び第２のＢモード画像データに表示された血管の走行方向と上述の投影成分の方向とが一致するように方向ベクトルの方向を更新し、更新された方向ベクトルの情報を速度補正部５の入射角度計測部５１へ供給する（図６のステップＳ７）。
【００６７】
そして、入射角度計測部５１は、方向ベクトル設定部９３から供給された方向ベクトルの設定情報とシステム制御部１０から供給されたサンプルゲートに対する超音波送信方向の情報に基づいて血管走行に対する超音波ビームの入射角度を計測する（図６のステップＳ８）。
【００６８】
一方、上述のステップＳ４において第１のＢモード画像データに対しサンプルゲートが設定されたならば、Ｂモード画像データ生成用の超音波送受信と交互して前記サンプルゲートの方向に対する超音波送受信が行なわれ、このとき、送受信部２の受信部２２から出力された受信信号は、超音波データ生成部３のスペクトラムデータ生成部３２へ供給される。
【００６９】
スペクトラムデータ生成部３２は、受信部２２から供給された受信信号をミキサ３２２及びＬＰＦ３２３により直交位相検波してドプラ信号を検出し、ＳＨ３２４へ供給する。ＳＨ３２４にはサンプルゲート設定部９２にて設定されたサンプルゲートの位置（深さ）に対応するサンプルゲートパルスがシステム制御部１０から供給され、このサンプルゲートパルスに従ってドプラ信号がサンプルホールドされる。そして、ＳＨ３２４の出力は、ＢＰＦ３２５において平滑化された後、Ａ／Ｄ変換器３２６においてデジタル信号に変換されてＦＦＴ分析器３２７の記憶回路に一旦保存される。
【００７０】
ＦＦＴ分析器３２７の演算回路は、前記記憶回路に保存された所定区間のドプラ信号をＦＦＴ分析して周波数スペクトラムを算出し、このＦＦＴ分析を一連のドプラ信号に対して所定時間間隔で繰り返す。そして、このとき得られた複数の周波数スペクトラムを時系列的に配列しドプラスペクトラムデータを生成する（図６のステップＳ９）。
【００７１】
次に、速度補正部５のスペクトラムデータ補正部５２は、スペクトラムデータ生成部３２によって生成されたドプラスペクトラムデータの血流速度あるいはドプラ周波数の値を入射角度計測部５１によって計測された超音波ビームの入射角度に基づいて補正し（図６のステップＳ１０）、補正後のドプラスペクトラムデータを表示部７において第１のＢモード画像データ及び第２のＢモード画像データと共にリアルタイム表示する（図６のステップＳ１１）。但し、ドプラスペクトラムデータは第１のＢモード画像データや第２のＢモード画像データに対して独立に表示しても構わない。
【００７２】
以上述べた本実施例によれば、超音波のドプラ効果を用いて血流速度を計測する際に、血流方向を３次元的に把握することができるため超音波ビームの入射角度に起因する計測誤差を高精度で補正することが可能となる。
【００７３】
特に、上述の実施例では、交叉するＢモード画像データの各々において血流方向を示す方向ベクトルの投影成分と血管像が重畳表示されるため、正確な方向ベクトルの設定が可能となり、精度の高い補正を行なうことができる。
【００７４】
又、交叉する複数のＢモード画像データの生成と表示は略リアルタイムで行なうことができるため、呼吸性あるいは拍動性の動きがある血管や臓器の血流計測においても高精度の補正が可能となる。
【００７５】
更に、上述の補正はボリュームデータを収集する必要がないため特許文献２に記載された方法と比較して短時間かつ簡単な操作によって行なうことができ、検査における操作者の負担が軽減される。
【００７６】
以上、本発明の実施例について述べてきたが、本発明は上述の実施例に限定されるものではなく、変形して実施することが可能である。例えば、補正されたドプラスペクトラムデータの表示に際し、ドプラスペクトラムデータの縦軸に示された血流速度あるいはドプラ周波数の目盛りを超音波ビームの入射角度に基づいて補正する場合について述べたが、前記入射角度に基づいてドプラスペクトラムデータの振幅を伸長あるいは圧縮して表示してもよい。
【００７７】
又、上述の実施例では、交叉する複数の画像断面において生成されたＢモード画像データを用いてサンプルゲートの設定や方向ベクトルの設定／更新等を行なったが、カラードプラ画像データあるいはＢモード画像データとカラードプラ画像データを用いて上述の設定を行なってもよい。特に、カラードプラ画像データを用いることにより、血流方向の把握が更に容易になる。尚、カラードプラ画像データの生成方法については、例えば、特開２００５−８１０８２号公報等に記載されているため詳細な説明は省略する。
【００７８】
更に、複数の画像断面に対して投影された方向ベクトルの投影成分に基づいて方向ベクトルの更新を行ったが、例えば、図７に示すようにサンプルゲートＳＧを含み複数画像断面の各々に対して垂直あるいは所定角度で交叉する平面（紙面）を設定し、この平面上に投影された前記方向ベクトルの観測を追加することにより方向ベクトルの設定／更新は更に容易となる。尚、図７において、ＩＰ１及びＩＰ２は第１の画像断面及び第２の画像断面を示し、ＤＶＰ１及びＤＶＰ２は第１の画像断面及び第２の画像断面に対する方向ベクトルＤＶの投影成分を示している。
【００７９】
一方、上述の実施例では、複数の振動素子が２次元配列された超音波プローブ１を用い交叉する２つのＢモード画像データを生成したが、複数の振動素子が１次元配列された超音波プローブあるいは単一の振動素子を有した超音波プローブを機械的に高速移動させて複数の画像データを生成してもよい。又、被検体に対する超音波走査方式は、実施例に示したセクタ走査方式に限定されるものではなく、リニア走査方式やコンベックス走査方式等の他の走査方式であっても構わない。
【００８０】
更に、操作者が入力部９の方向ベクトル設定部９３を操作することにより方向ベクトルの設定／更新を行なう場合について述べたが、複数の画像データにおける血管断面データと前記画像データの各々に投影された方向ベクトルの投影成分を画像処理することによって方向ベクトルの設定／更新を自動的に行なってもよい。
【００８１】
又、上述の実施例では、超音波のドプラ効果を用いて得られる被検体の血流速度を補正する場合について述べたが、生体組織等の移動速度に対して同様の補正を行なってもよい。
【図面の簡単な説明】
【００８２】
【図１】本発明の実施例における超音波診断装置の全体構成を示すブロック図。
【図２】同実施例の超音波診断装置を構成する送受信部及び超音波データ生成部の構成を示すブロック図。
【図３】同実施例におけるスペクトラムデータ生成部の基本動作を示すタイムチャート。
【図４】同実施例における超音波ビームの入射角度計測方法を説明するための図。
【図５】同実施例におけるＢモード画像データ及びドプラスペクトラムデータの表示例を示す図。
【図６】同実施例におけるドプラスペクトラムデータの補正手順を示すフローチャート。
【図７】同実施例における方向ベクトル投影方法の変形例を示す図。
【符号の説明】
【００８３】
１…超音波プローブ
２…送受信部
２１…送信部
２２…受信部
３…超音波データ生成部
３１…Ｂモードデータ生成部
３２…スペクトラムデータ生成部
４…画像データ生成部
５…速度補正部
５１…入射角度計測部
５２…スペクトラムデータ補正部
６…投影成分演算部
７…表示部
８…基準信号発生部
９…入力部
９１…画像断面設定部
９２…サンプルゲート設定部
９３…方向ベクトル設定部
１０…システム制御部
１００…超音波診断装置
【図１】
2008073422.tif 000005
【図２】
2008073422.tif 000006
【図３】
2008073422.tif 000007
【図４】
2008073422.tif 000008
【図５】
2008073422.tif 000009
【図６】
2008073422.tif 000010
【図７】
2008073422.tif 000011