(57)【要約】【課題】血流からのドップラスペクトラムに弁や腱索の影響が混入している場合であっても、Ｅ波の減速時間ＤＣＴを高い信頼性で自動計測することができる超音波診断装置を提供する。【解決手段】本発明に係る超音波診断装置は、心臓内血流のドップラスペクトラムを生成するドップラ処理部と、ドップラスペクトラムの時間的変化をトレース波形として生成するトレース波形生成部と、トレース波形から所定のピーク波形を検出するピーク波形検出部と、ピーク波形の減速領域の波形に血流以外の影響による変動があると判断される場合は、この変動を取り除いて評価用減速波形を生成し、それ以外の場合は検出したピーク波形の減速領域の波形を評価用減速波形とする評価波形生成部と、評価用減速波形から、前記所定のピーク波形の減速時間を算出する減速時間算出部とを備える。【選択図】  図１
2010046284.tif 000002
【特許請求の範囲】
【請求項１】
心臓内の所定の部位から反射される超音波信号からドップラスペクトラムを生成するドップラ処理部と、
前記ドップラスペクトラムにおける所定のスペクトラム成分の時間的変化をトレース波形として生成するトレース波形生成部と、
生成された前記トレース波形から所定のピーク波形を検出するピーク波形検出部と、
前記ピーク波形の減速領域の波形に血流以外の影響による変動があると判断される場合は、この変動を取り除いて評価用減速波形を生成し、それ以外の場合は検出した前記ピーク波形の減速領域の波形を前記評価用減速波形とする評価波形生成部と、
前記評価用減速波形から、前記所定のピーク波形の減速時間を算出する減速時間算出部と、
を備えたことを特徴とする超音波診断装置。
【請求項２】
前記所定のピーク波形は、心臓拡張期のＥ波の波形であり、前記減速時間は、Ｅ波の減速時間ＤＣＴである、
ことを特徴とする請求項１に記載の超音波診断装置。
【請求項３】
心電波形を入力する心電波形入力部と、
複数の心電波形を取得する一方、血流以外の影響を受けていないピーク波形を各心電波形に同期させて取得し、取得した前記複数の心電波形と前記複数の純正血流のピーク波形とから、パラメトリックシステム同定モデルのパラメータを予め求めて記憶するパラメータ記憶部と、
をさらに備え、
前記評価波形生成部は、
実時間で入力される前記心電波形と、前記パラメータ記憶部に記憶されている前記パラメトリックシステム同定モデルのパラメータとから推定ピーク波形を実時間で生成し、
実時間で得られる前記ピーク波形と前記推定ピーク波形との振幅差が所定の閾値を超えている場合は、前記血流以外の影響による変動があると判断し、
前記閾値を超えている期間の波形を前記ピーク波形から除去することによって前記評価用減速波形を生成する、
ことを特徴とする請求項１に記載の超音波診断装置。
【請求項４】
前記評価波形生成部は、
前記ピーク波形から除去された期間の波形を、除去後に残っている前記ピーク波形によって補間し、前記評価用減速波形を生成する、
ことを特徴とする請求項３に記載の超音波診断装置。
【請求項５】
前記評価波形生成部は、
前記ピーク波形から除去された期間の波形を、該当する期間の前記推定ピーク波形と置換して、前記評価用減速波形を生成する、
ことを特徴とする請求項３に記載の超音波診断装置。
【請求項６】
前記評価波形生成部は、
抽出した前記ピーク波形のピーク値と、このピーク値よりも所定の比率だけ小さな中間値との間に極小値があるか否かを判定し、
前記極小値がある場合には、前記血流以外の影響による変動があると判断し、
前記ピーク波形の前記ピーク値から前記極小値までの範囲の波形を前記評価用減速波形として生成する、
ことを特徴とする請求項１に記載の超音波診断装置。
【請求項７】
前記減速時間算出部は、
前記評価用減速波形のピーク値と、このピーク値よりも所定の比率だけ小さな中間値との間の波形を近似する近似直線を求め、前記近似直線がゼロクロスする時刻と前記評価用減速波形のピーク値に対応する時刻との差から前記減速時間を算出する、
ことを特徴とする請求項１に記載の超音波診断装置。
【請求項８】
前記減速時間算出部は、
前記ピーク値と前記中間値との間の波形データに基づく最小自乗法によって前記近似直線を求める、
ことを特徴とする請求項７に記載の超音波診断装置。
【請求項９】
前記減速時間算出部は、
前記評価用減速波形のピークの座標が時間ゼロ、振幅ゼロの原点となるように前記評価用減速波形を平行移動させ、移動後の前記ピーク値と前記中間値との間の波形データに基づく最小自乗法によって、１次項のみからなる原点通過近似直線を求め、求めた前記原点通過近似直線を先に移動させた量だけ逆に平行移動させて前記近似直線を求め、これにより前記最小自乗法の演算量を低減させる、
ことを特徴とする請求項８に記載の超音波診断装置。
【請求項１０】
心臓内の所定の部位から反射される超音波信号からドップラスペクトラムを生成し、
前記ドップラスペクトラムにおける所定のスペクトラム成分の時間的変化をトレース波形として生成し、
生成された前記トレース波形から所定のピーク波形を検出し、
前記ピーク波形の減速領域の波形に血流以外の影響による変動があると判断される場合は、この変動を取り除いて評価用減速波形を生成し、それ以外の場合は抽出した前記ピーク波形を前記評価用減速波形とし、
前記評価用減速波形から、前記所定のピーク波形の減速時間を算出する、
ステップを備えたことを特徴とする診断パラメータ自動計測方法。
【請求項１１】
前記所定のピーク波形は、心臓拡張期のＥ波の波形であり、前記減速時間は、Ｅ波の減速時間ＤＣＴである、
ことを特徴とする請求項１０に記載の診断パラメータ自動計測方法。
【請求項１２】
心電波形を入力し、
複数の心電波形を取得する一方、血流以外の影響を受けていないピーク波形を各心電波形に同期させて取得し、取得した前記複数の心電波形と前記複数のピーク波形とから、パラメトリックシステム同定モデルのパラメータを予め求めて記憶する、
ステップをさらに備え、
前記評価用減速波形を生成するステップでは、
実時間で入力される前記心電波形と、前記パラメトリックシステム同定モデルのパラメータとから推定ピーク波形を実時間で生成し、
実時間で得られる前記ピーク波形と前記推定ピーク波形との振幅差が所定の閾値を超えている場合は、前記血流以外の影響による変動があると判断し、
前記閾値を超えている期間の波形を前記ピーク波形から除去することによって前記評価用減速波形を生成する、
ことを特徴とする請求項１０に記載の診断パラメータ自動計測方法。
【請求項１３】
前記評価用減速波形を生成するステップでは、
前記ピーク波形から除去された期間の波形を、除去後に残っている前記ピーク波形によって補間し、前記評価用減速波形を生成する、
ことを特徴とする請求項１２に記載の診断パラメータ自動計測方法。
【請求項１４】
前記評価用減速波形を生成するステップでは、
前記ピーク波形から除去された期間の波形を、該当する期間の前記推定ピーク波形と置換して、前記評価用減速波形を生成する、
ことを特徴とする請求項１２に記載の診断パラメータ自動計測方法。
【請求項１５】
前記評価用減速波形を生成するステップでは、
抽出した前記ピーク波形のピーク値と、このピーク値よりも所定の比率だけ小さな中間値との間に極小値があるか否かを判定し、
前記極小値がある場合には、前記血流以外の影響による変動があると判断し、
前記ピーク波形の前記ピーク値から前記極小値までの範囲の波形を前記評価用減速波形として生成する、
ことを特徴とする請求項１０に記載の診断パラメータ自動計測方法。
【請求項１６】
前記減速時間を算出するステップでは、
前記評価用減速波形のピーク値と、このピーク値よりも所定の比率だけ小さな中間値との間の波形を近似する近似直線を求め、前記近似直線がゼロクロスする時刻と前記評価用減速波形のピーク値に対応する時刻との差から前記減速時間を算出する、
ことを特徴とする請求項１０に記載の診断パラメータ自動計測方法。
【請求項１７】
前記減速時間を算出するステップでは、
前記ピーク値と前記中間値との間の波形データに基づく最小自乗法によって前記近似直線を求める、
ことを特徴とする請求項１６に記載の診断パラメータ自動計測方法。
【請求項１８】
前記減速時間を算出するステップでは、
前記評価用減速波形のピークの座標が時間ゼロ、振幅ゼロの原点となるように前記評価用減速波形を平行移動させ、移動後の前記ピーク値と前記中間値との間の波形データに基づく最小自乗法によって、１次項のみからなる原点通過近似直線を求め、求めた前記原点通過近似直線を先に移動させた量だけ逆に平行移動させて前記近似直線を求め、これにより前記最小自乗法の演算量を低減させる、
ことを特徴とする請求項１７に記載の診断パラメータ自動計測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波診断装置及び診断パラメータ自動計測方法に係り、特に、生体内の血流速度をドップラ情報に基づいて計測することができる超音波診断装置及び診断パラメータ自動計測方法に関する。
【背景技術】
【０００２】
超音波診断装置は、超音波プローブに内蔵された圧電振動子から発生する超音波信号を被検体内に送信し、被検体内からの反射信号を画像化して各種の診断を行うための装置である。診断用の画像としてＢモード画像やカラードップラ画像が広く用いられている。
【０００３】
一方、被検体の血管内や心臓内の血流情報を定量的にかつ精度よく測定する方法として、ドップラスペクトラム法がある。ドップラスペクトラム法では、血管内や心臓内の所望の方向に超音波ビームを向けて送信し、受信信号から所望の深度の信号をレンジゲートによって切り出す。切り出した受信信号の系列に対してＦＦＴ等の演算を施してドップラスペクトラムを求め、さらに求めたドップラスペクトラム（ドップラ周波数）から血流速度や速度成分の強さ等の血流情報を得ている。
【０００４】
心臓や血管の内部のどの位置の血流情報を得るかは予め指定する必要があるが、この指定は、通常Ｂモード画像を用いて行う。Ｂモード画像上で特定の方向やレンジを指定し、指定した方向とレンジからの受信信号を抽出してその部位のドップラスペクトラムを求めている。
【０００５】
超音波診断装置によって得られるドップラスペクトラムは、通常、横軸に時間、縦軸に周波数をとって画像化し、また周波数成分毎の強度を画像の輝度として表現している。この画像を、ドップラ画像或いはドップラ波形画像と呼んでいる。
【０００６】
また、ドップラスペクトラムから血流の最高流速Ｖｐや平均流速Ｖｍを求め、横軸に時間、縦軸に最高流速Ｖｐや平均流速Ｖｍをとった波形も求めている。この波形をトレース波形と呼んでいる。トレース波形はハードディスク等の適宜の記憶手段に保存され、必要に応じて超音波診断装置の表示部に表示される。
【０００７】
心臓の機能を評価し或いは診断する場合に、上記のトレース波形から得られる診断パラメータが非常に有用となる。
【０００８】
例えば、心機能を評価する場合、左室流入血流（LV-Inflow）流速や、肺静脈血流（PV）流速等に対してトレース波形を求め、このトレース波形から各種診断パラメータを求めている。
【０００９】
左室流入血流波形に対しては、Ｅ波（拡張早期波）やＡ波（心房収縮期波）の振幅（振幅が血流速度に対応する）、Ｅ波とＡ波の振幅比Ｅ／Ａ、ＤＣＴ（Ｅ波の減速時間（Deceleration  Time））等が診断パラメータとして広く用いられている。
【００１０】
また、肺静脈血流波形に対しては、Ｓ波（収縮期順行性血流）、Ｄ波（拡張期順行性血流）、ＡＲ波（収縮期逆行性血流）の振幅等が診断パラメータとして用いられている。
【００１１】
従来、これらの診断パラメータは、医師や技師等のユーザによるマニュアル操作によってトレース波形から読み取られていた。例えば、左室流入血流波形から診断パラメータを読み取る場合には、表示部に静止表示したトレース波形に対して２つのタイムカーソルを配置し、２つのタイムカーサの位置をＥ波とＡ波のピークに夫々ユーザが一致させる。そして、２つのタイムカーサの位置にあるＥ波の振幅とＡ波の振幅を装置が読み取り、診断パラメータとして表示部に表示する。また、読み取ったＥ波とＡ波の振幅から、振幅比Ｅ／Ａを算出し、診断パラメータとして表示する。
【００１２】
また、減速時間ＤＣＴを求める場合は、一方のタイムカーサをＥ波のピークに合わせ、他方のタイムカーサをＥ波が減速してほぼゼロとなる位置に合わせる。そして、装置が２つのタイムカーサの時間差を読みとって、診断パラメータＤＣＴを表示する。
【００１３】
このように、従来は、ユーザのマニュアル操作によって診断パラメータを求めていたため、操作が煩雑となるばかりでなく、診断パラメータを得るまでに多くの時間を要していた。
【００１４】
さらに、ドップラスペクトラム自体はほぼリアルタイムで得られるにもかかわらず、診断パラメータの取得は上記のようにマニュアル操作が介在するため、診断パラメータをリアルタイムで取得することは事実上不可能であった。
【００１５】
このような不都合を解消するため、診断パラメータを自動計測する技術が特許文献１等に開示されている。
【００１６】
特許文献１では、トレース波形から極大点、極小点のペアを複数求め、これらのペアの中からＥ波とＡ波のピーク時刻とピーク値を夫々自動的に求める技術等が開示されている。また、Ｅ波のピーク値からの下降曲線に対して接線を設定し、この接線とベースラインとが交わる時刻とＥ波のピーク時刻との間隔から減速時間ＤＣＴを求める技術も開示されている。
【特許文献１】特開２００６−１０２４８９号公報
【発明の開示】
【発明が解決しようとする課題】
【００１７】
しかしながら、超音波診断装置によって得られるドップラスペクトラムは、必ずしも必要な情報のみが含まれているとは限らない。
【００１８】
例えば、左室流入血流の速度情報を得ようとする場合、左心室入り口近傍の反射信号からドップラスペクトラムを取得することが多いが、左心室と左心房の間には僧帽弁と呼ばれる弁や僧帽弁を支持する多数の腱索が存在する。このため、計測するドップラスペクトラムには、評価対象とする左室流入血流のドップラ信号の他に、僧帽弁や腱索の動きに対応するドップラ信号が不要ドップラ信号として混在することになる。
【００１９】
例えば、左室流入血流のＥ波の減速時間ＤＣＴを測定しようとした場合、その測定対象範囲に、僧帽弁や腱索の影響が混入している割合は、３０％程度もあると言われている。
【００２０】
このため、特許文献１等に開示されている従来の方法でＤＣＴを自動計測することは、現実問題としてはほぼ困難な状況となっていた。
【００２１】
Ｅ波の減速時間ＤＣＴは、心機能を評価する上で重要な診断パラメータの１つであり、信頼性の高いＤＣＴの自動計測技術が強く求められている。
【００２２】
本発明は、上記事情に鑑みてなされたものであり、血流からのドップラスペクトラムに弁や腱索の影響が混入している場合であっても、Ｅ波の減速時間ＤＣＴを高い信頼性で自動計測することができる超音波診断装置及び診断パラメータ自動計測方法を提供することを目的とする。
【課題を解決するための手段】
【００２３】
上記課題を解決するため、本発明に係る超音波診断装置は、心臓内の所定の部位から反射される超音波信号からドップラスペクトラムを生成するドップラ処理部と、前記ドップラスペクトラムにおける所定のスペクトラム成分の時間的変化をトレース波形として生成するトレース波形生成部と、生成された前記トレース波形から所定のピーク波形を検出するピーク波形検出部と、前記ピーク波形の減速領域の波形に血流以外の影響による変動があると判断される場合は、この変動を取り除いて評価用減速波形を生成し、それ以外の場合は検出した前記ピーク波形の減速領域の波形を前記評価用減速波形とする評価波形生成部と、前記評価用減速波形から、前記所定のピーク波形の減速時間を算出する減速時間算出部と、を備えたことを特徴とする。
【００２４】
また、本発明に係る診断パラメータ自動計測方法は、心臓内の所定の部位から反射される超音波信号からドップラスペクトラムを生成し、前記ドップラスペクトラムにおける所定のスペクトラム成分の時間的変化をトレース波形として生成し、生成された前記トレース波形から所定のピーク波形を検出し、前記ピーク波形の減速領域の波形に血流以外の影響による変動があると判断される場合は、この変動を取り除いて評価用減速波形を生成し、それ以外の場合は抽出した前記ピーク波形を前記評価用減速波形とし、前記評価用減速波形から、前記所定のピーク波形の減速時間を算出する、ステップを備えたことを特徴とする。
【発明の効果】
【００２５】
本発明に係る超音波診断装置、及び診断パラメータ自動計測方法によれば、血流からのドップラスペクトラムに弁や腱索の影響が混入している場合であっても、Ｅ波の減速時間ＤＣＴを高い信頼性で自動計測することができる。
【発明を実施するための最良の形態】
【００２６】
本発明に係る超音波診断装置、及び診断パラメータ自動計測方法の実施形態について添付図面を参照して説明する。
【００２７】
（１）構成
図１は、本実施形態に係る超音波診断装置１の構成例を示すブロック図である。超音波診断装置１は、超音波プローブ２、基準信号生成部３、送受信部４、システム制御部５、操作部６、ドップラ信号検出部７、ドップラ処理部８、カラードップラデータ生成部９、Ｂモードデータ生成部１０、表示部１１、記憶部１２、血流評価部２０等を備えて構成されている。
【００２８】
超音波プローブ２は、被検体の体表面に先端部を接触させて超音波の送受信を行うデバイスであり、１次元或いは２次元に配列された複数の微小な圧電振動子をその先端部に有している。圧電振動子は電気音響変換素子であり、送信時には電気信号を超音波パルスに変換し、受信時には被検体から反射された超音波信号を電気信号に変換する。超音波プローブ２には、リニア走査対応、セクタ走査対応、コンベックス走査対応等のものがあり、診断部位に応じて選択される。また、圧電振動子が２次元配列されている超音波プローブ２では、いわゆる３次元走査が可能である。超音波プローブ２は、ケーブルを介して送受信部４に接続されている。
【００２９】
送受信部４は、その内部構成として送信ユニットと受信ユニット（いずれも図示を省略）を有している。送信ユニットは、基準信号生成部３から供給される基準信号（連続波或いは基準クロック）に基づいて、送信パルスの繰り返し周期を決定するレートパルスを生成する。さらに、超音波ビームの走査角（送信走査角）に対応する遅延時間をレートパルスに与え、その後パルサ（図示せず）によってレートパルスを駆動パルスに変換し、この駆動パルスを超音波プローブ２の各圧電振動子に供給している。
【００３０】
一方、受信ユニットでは、超音波プローブ２の各圧電振動子から出力される電気信号をデジタル信号に変換した後、超音波ビームの走査角（受信走査角）に対応する遅延量を印加して加算することによって、ビームフォーミングされた受信信号を得ている。
【００３１】
システム制御部５は、ユーザが操作部６に対して設定し、或いは入力した情報に基づいて、ビーム走査や送受信に関する諸元を決定し、これらの諸元を送受信部４に対して設定している。
【００３２】
送受信部４から出力される受信信号は、Ｂモードデータ生成部１０やドップラ信号検出部７に出力される。
【００３３】
Ｂモードデータ生成部では、受信信号を例えば振幅検波して、Ｂモード画像の表示に必要となる走査角ごとのレンジ対振幅情報のデータ（即ちＢモードデータ）を生成して記憶部１２に一時的に保存している。このＢモードデータは記憶部１２から読み出され、表示部１１にＢモード画像として表示される。ちなみに、図２は、Ｂモード画像の表示例として、心臓のＢモード画像を模式的に示した図である。
【００３４】
一方、ドップラ信号検出部７では、受信信号に対して位相検波を行っており、位相検波信号をカラードップラデータ生成部９とドップラ処理部８に供給している。
【００３５】
カラードップラデータ生成部９では、位相検波信号から自己相関値を求め、自己相関値から平均ドップラ周波数（平均血流速度に対応する量）や周波数の分散値を算出し、これらのカラードップラデータを対応する走査角データと共に記憶部１２に一時的に記憶している。カラードップラ表示を行う場合には、記憶部１２からカラードップラデータが読み出され、Ｂモード画像と同様な２次元画像の上に平均血流速度に対応するカラーの輝度表示がなされる。
【００３６】
ドップラ処理部８では、血流速度を詳細に評価するため、位相検波信号に対してＦＦＴ等の周波数解析処理を行ってドップラスペクトラムを生成している。ドップラスペクトラムは、被検体内の特定の部位を走査角とレンジによって指定し、この特定の部位からの反射信号を所定の期間収集したデータから求めている。
【００３７】
ドップラスペクトラムは、通常横軸に時間、縦軸に周波数をとって画像化し、また周波数成分毎の強度を画像の輝度として表現している。この画像を、ドップラ画像或いはドップラ波形画像と呼んでいる。ドップラ処理部８から出力されたドップラスペクトラムは記憶部１２に記憶され、表示部１１にドップラ画像として表示される。
【００３８】
また、ドップラスペクトラムから血流の最高流速Ｖｐや平均流速Ｖｍを求め、横軸に時間、縦軸に最高流速Ｖｐや平均流速Ｖｍをとった波形も生成している。この波形をトレース波形と呼んでいる。血流評価部２０のトレース波形生成部２１は、記憶部１２からドップラスペクトラムを読み出し、このドップラスペクトラムからトレース波形を生成する。生成されたトレース波形は再び記憶部１２に保存され、必要に応じて表示部１１に表示される。
【００３９】
本実施形態に係る超音波診断装置１では、血流のトレース波形から診断上有用な診断パラメータを自動的に計測する機能を有しており、この自動計測機能を血流評価部２０で実現している。血流評価部２０は、上述したトレース波形生成部２１の他、ピーク波形検出部２２、評価波形生成部２３、減速時間算出部２４、パラメータ記憶部２５、心電波形入力部２６等を有しており、これらの各部の機能の詳細については後述する。
【００４０】
（２）診断パラメータの自動計測
図２は、心臓の模式的なＢモード画像と、Ｂモード画像上で指定される血流の評価部位とを例示している図である。図２では、左心房から左心室へ流入する近傍が走査角とレンジゲートによって指定されており、左室流入血流を評価対象としている例を示している。
【００４１】
本実施形態に係る超音波診断装置１の血流評価対象はこの左室流入血流に限定されるものではないが、左室流入血流の評価は心機能の評価上非常に重要なものであるため、以下の説明では主として左室流入血流を評価対象の例としてとりあげて説明する。
【００４２】
心臓は収縮と拡張を繰り返しているが、拡張期には左心房から左心室に血液が流れる。この血流流速を左心室の入り口近傍を観測対象部位として観測したドップラスペクトラムが左室流入血流のドップラスペクトラムであり、このドップラスペクトラムの例えばピーク値をトレースした波形が左室流入血流のトレース波形である。
【００４３】
図３は、左室流入血流のトレース波形の実例を心臓の拡張期と収縮期を示す信号と併記して示している。拡張期と収縮期を示す信号は、心電計１００（図１参照）から出力される心電波形（ＥＣＧ波形ともいう）や心音波形から求めることができる。
【００４４】
左室流入血流のトレース波形（以下、単にトレース波形という）の縦軸は、ドップラ周波数から求めた左室流入血流の血流速度（図３に示す単位はｃｍ／ｓ）を示しており、横軸は時間である。
【００４５】
図３に示したように、拡張期のトレース波形には、Ｅ波とＡ波の２つの特徴的なピーク波形が存在する。拡張期の早期には、左心室が拡張することによって左心室の圧力が低下し、血液が左心房から左心室に流入する。このときの血流に対応する波形がＥ波（拡張早期波）である。拡張期の後期には、左心房が収縮して左心房に残っている血液を左心室に押し込む。このときの血流に対応する波形がＡ波（心房収縮期波）である。つまり、Ｅ波は左心室拡張による血液の左室流入を、Ａ波は左心房収縮による血液の左室流入を表している。
【００４６】
拡張期左室流入波形に基づく診断パラメータとしては、Ｅ波のピーク値（Ｅｐ）、Ａ波のピーク値（Ａｐ）、Ｅ波のピーク値とＡ波のピーク値との比（Ｅ／Ａ）、Ｅ波の減速時間ＤＣＴ（Deceleration  Time）等がある。Ｅ波の減速時間ＤＣＴは、Ｅ波のピークから血流速度がゼロとなるまでの時間として定義されるものである。これらの診断パラメータの値によって心機能を評価することができる。
【００４７】
ところで、図２に示したように、左室流入血流の測定対象部位の近傍（左心房と左心室の間）には、僧帽弁と呼ばれる弁があり、また僧帽弁を支持するための多数の腱索が存在する。僧帽弁や腱索は動きを伴うため、ドップラスペクトラムやトレース波形には血流以外の影響として弁や腱索の影響による変動波形（不要波形）が混入する。
【００４８】
図３に示した４つの拡張期のトレース波形のうち、左から３つの拡張期波形は血流のみによる波形（Ｅ波とＡ波）が支配的であるが、右端の拡張期のトレース波形には、Ｅ波の減速期間中に、弁や腱索の影響と見られる顕著な不要波形が観測されている。
【００４９】
多数の実測データを検証した結果、拡張期のトレース波形のうち３０％程度のトレース波形に図３の右端に示したように弁や腱索の影響を受けたトレース波形が含まれることが判明している。
【００５０】
弁や腱索の影響を受けた波形に基づいて診断パラメータを自動計測すると、誤った診断パラメータを出力することになる。特に、Ｅ波の減速時間ＤＣＴの自動計測に関してはその誤差は大きなものとなる。以下の説明では、診断パラメータのうち、Ｅ波の減速時間ＤＣＴの自動計測に焦点を絞り、どのようして弁や腱索の影響を排除し、どのようにして信頼性の高い減速時間ＤＣＴを自動計測によって得るかについて順次説明する。
【００５１】
図４は、本実施形態に係る評価パラメータ（特にＥ波の減速時間ＤＣＴ）の自動計測処理の一例を示すフローチャートである。
【００５２】
ステップＳＴ１では、心臓内の所定の部位（例えば、左心室入り口近傍）に超音波ビームの走査角とレンジゲートを設定し、その部位から反射される超音波信号からドップラスペクトラムを生成する。ドップラスペクトラムの生成はドップラ処理部８で行う。
【００５３】
ステップＳＴ２では、ドップラスペクトラムからトレース波形を生成する。具体的には、ドップラスペクトラムのピークを検出して時間軸方向にトレースした波形を生成する。トレース波形の生成は、血流評価部２０のトレース波形生成部２１で行う。
【００５４】
ステップＳＴ３では、トレース波形（この場合、左室流入血流のトレース波形）から所定のピーク波形、即ちＥ波とＡ波を検出する。
【００５５】
Ｅ波とＡ波の検出は、例えば特許文献１に開示されている技術を利用することができる。この技術では、まず、心電信号（ＥＣＧ信号）のＲ波から心拍周期の４０％を基準位置とする。この基準位置を収縮期から拡張期への移行時刻とみなす。次に、トレース波形の極大／極小ペアを検索して、最大の極大値と、２番目に大きな極大値とを有する２つの座標を検出し、前記基準位置に後続する極大の座標をＥ波のピークとして選定し、このＥ波のピークに後続する極大の座標をＡ波のピークとして選定する。特許文献１には、この他にも心音波形（ＰＣＧ波形）を用いて上記の基準位置を求める手法や、トレース波形のみからＥ波とＡ波の各ピークの座標を求める手法等が開示されている。ステップＳＴ３では、これらの手法を用いてＥ波のピーク座標とＡ波のピーク座標を検出する。また、ステップＳＴ３の処理、即ちとレース波形からピーク波形を検出し、Ｅ波のピーク座標とＡ波のピーク座標の検出する処理は、血流評価部２０のピーク波形検出部２２で行っている。
【００５６】
ステップＳＴ４では、Ｅ波の減速領域の波形に血流以外の影響（即ち、弁や腱索の影響）による変動があるか否かを判定する。この判定の具体的な方法については後述する。ステップＳＴ４で、血流以外の影響がＥ波の減速領域の波形には存在しないと判定された場合にはステップＳＴ５へ進む。
【００５７】
ステップＳＴ５では、Ｅ波の減速波形の所定の領域を抽出して評価用減速波形とする。具体的には、Ｅ波のピーク値（Ｅｐ）よりも所定の比率だけ小さい、例えばＥ波のピーク値の１／２の値の中間値（Ｅｐ／２）を求め、ピーク値（Ｅｐ）と中間値（Ｅｐ／２）の間にある減速波形を抽出して評価用減速波形とする。
【００５８】
ステップＳＴ７では、ステップＳＴ５で求めた評価用減速波形のデータ（ピーク値（Ｅｐ）と中間値（Ｅｐ／２）の間のデータ）に基づいて、評価用減速波形の近似直線を求める。この近似直線は、例えば最小自乗法等の手法で求める。さらに、近似直線がゼロクロスする時刻を求め、ゼロクロスの時刻とＥ波のピークの時刻との差をＥ波の減速時間ＤＣＴとして算出する。上述したステップＳＴ４及びステップＳＴ５の処理は評価波形生成部２３で行っており、ステップＳＴ７の処理は減速時間算出部２４で行っている。
【００５９】
図５の中央部に示す説明図は、弁や腱索の影響が無い場合の波形に対して上述したステップＳＴ５及びステップＳＴ７の処理を行い、Ｅ波の減速時間ＤＣＴが良好に求まることを示している。
【００６０】
これに対して、図５の右側に示す説明図は、弁や腱索の影響が有る波形に対して上述したステップＳＴ５及びステップＳＴ７の処理をそのまま行うと、減速時間ＤＣＴに大きな誤差が発生することを示している。この場合、Ｅ波のピーク値（Ｅｐ）と中間値（Ｅｐ／２）の間には、弁や腱索の影響に起因する不要なピーク波形が存在し、最小自乗法に使用するデータにはこの不要なピーク波形のデータが含まれる。このため、求めた近似直線は弁や腱索の影響を受けたものとなり、この近似直線から算出される減速時間ＤＣＴは大きな誤差をもつことになる。
【００６１】
そこで、本実施形態に係る超音波診断装置１では、弁や腱索の影響を排除する方法を２つの方法で実現している。第１の排除方法はパラメトリックシステム同定モデルを利用する排除方法であり、第２の排除方法は極小値検出を利用する排除方法である。
【００６２】
なお、パラメトリックシステム同定モデルにはいくつかの種類が存在するが、以下の説明ではパラメトリックシステム同定モデルの中のＡＲＸモデルを例にとって説明する。
【００６３】
（３）第１の排除方法
図６は、ＡＲＸモデルの利用によって弁や腱索の影響を排除する第１の排除方法の概念を示す図である。
【００６４】
この方法では、ＡＲＸモデルに基づく判定処理によって弁や腱索の影響を受けているか否かを判定し、影響をうけている場合には影響を受けている領域を特定する（ステップＳＴ４）。そして、影響を受けている領域をＥ波の減速波形から除去し、さらに除去した領域を線形補間する、或いはＡＲＸモデルに基づく推定トレース波形で置換することによって評価用減速波形を生成する（ステップＳＴ６）。なお、第１の排除方法における上記ステップＳＴ４、ステップＳＴ６の処理も評価波形生成部２３で行っている。第１の排除方法によるステップＳＴ４とステップＳＴ６のより詳細な処理を以下に説明する。
【００６５】
第１の排除方法では、図７に示したように、複数の心電波形（ＥＣＧ波形）ｕ（ｔ）と、血流以外の影響（即ち、弁や腱索の影響）を受けていない複数の実測トレース波形Ｙi(ｔ)を前記心電波形ｕ(t)と同期させて予め収集しておく。血流以外の影響を受けていないトレース波形Ｙi(ｔ)は、例えば、弁や腱索の存在しない心臓の部位で血流測定を行って収集することができる。或いは、弁や腱索の影響を受けているトレース波形から、ベテランの医師や技師が手動でその影響部位を除去することによって得ることができる。
【００６６】
収集したこれらのトレース波形Ｙi(t)と心電波形ｕ(t)を基にＡＲＸモデルによるシステム同定を行う。
【００６７】
ここで、ＡＲＸモデルは、次の（式１）から（式３）で表される。  
［数１］
Ａ(q)*Ｙi(t) ＝ Ｂ(q)*ｕ(t-nk)＋ｅ(t)              （式１）
Ａ(q) ＝ １＋a１q-1＋・・・anaq-na                    （式２）
Ｂ(q) ＝ b１＋b２q-1＋・・・bnbq-nb+1                  （式３）
【００６８】
（式１）から（式３）において、ｅ(t)は残差、即ち期待値と実測値との差、Ａ(q)及びＢ(q)は、既約なシフトオペレータqの多項式、nkは、心電波形ｕ(t)に対応する血流以外の影響を受けていないトレース波形Ｙi(t)との時間の遅れ、及びna、nbは整数の引数である。また、（式２）及び（式３）における
［数２］
ai ＝ (a１, a２, ・・・・ana)                        （式４）
bj ＝ (b１, b２, ・・・・bnb)                        （式５）
が血流以外の影響を受けていないときのＡＲＸモデルパラメータである。
【００６９】
収集した多数の血流以外の影響を受けていないトレース波形Ｙi(ｔ)、いわば理想的なトレース波形Ｙi(ｔ)と心電波形ｕ(t)によって（式４）、（式５）に示すＡＲＸモデルパラメータを予め求めておく。そして予め求めたＡＲＸモデルパラメータを血流評価部２０のパラメータ記憶部２５に保存しておく。
【００７０】
図８は、ステップＳＴ４における判定の概念を説明する図である。ステップＳＴ４における判定では、リアルタイムで得られる患者の実測トレース波形（これをＹp（t）と記述する）とＡＲＸモデルによって推定される推定トレース波形（これをＹe(t)と記述する）とを比較する。推定トレース波形Ｙe(t)は、リアルタイムで得られる患者の実測トレース波形Ｙp（t）と、これと同時に得られる実測心電波形ｕ(t)、及びパラメータ記憶部２５に保存されているＡＲＸモデルパラメータとから、次の（式６）によって求められる。  
［数３］
Ｙe(t) ＝ −Σ{ ai*Ｙp（t-i）} +Σ{ bj*ｕ(t-j）}    （式６）
（式６）で使用するＡＲＸモデルパラメータは、弁や腱索の影響を受けていない理想的な状態で収集されたトレース波形Ｙi(t)によって決定されたものであり、（式６）で得られる推定トレース波形Ｙe(t)も、実測トレース波形Ｙp（t）から弁や腱索の影響が取り除かれたものとして推定される波形となる。
【００７１】
従って、この推定トレース波形Ｙe(t)と実測トレース波形Ｙp（t）との差分をとり、差分の絶対値ＡＢＳ（Ｙe(t)−Ｙp（t））を所定の閾値と比較することによって実測トレース波形Ｙp（t）における弁や腱索の影響を判定することができる。即ち、ＡＢＳ（Ｙe(t)−Ｙp（t））が閾値より小さい場合は弁や腱索の影響を受けていないと判定し、ＡＢＳ（Ｙe(t)−Ｙp（t））が閾値より小さい場合は、その期間を弁や腱索の影響を受けている期間であると判定する。以上が第１の排除方法を用いた場合のステップＳＴ４の処理である。
【００７２】
次にステップＳＴ６では、弁や腱索の影響を排除した評価用減速波形を生成する。図９は、ステップＳＴ６の処理の概念を示す図である。
【００７３】
すでにステップＳＴ４で弁や腱索の影響を受けている期間が求められているため、実測トレース波形Ｙp（t）からこの期間の波形を除去する（図９の中央の図）。そして、除去された期間の波形を直線で補間して評価用減速波形を生成する（図９の右側上部の図）。
【００７４】
或いは、除去された期間を、該当する期間の推定トレース波形Ｙe(t)で置換して評価用減速波形を生成してもよい（図９の右側下部の図）。
【００７５】
このようにして弁や腱索の影響が排除された評価用減速波形が生成された後は、前述した方法と同様の方法によって減速期間の近似直線を生成して減速時間ＤＣＴを算出する（ステップＳＴ７）。
【００７６】
上記の説明では、ＡＲＸモデルを例として説明したが、このモデルに限定されるものではなく、ＡＲＸモデル以外のパラメトリックシステム同定モデルを用いた方法でも良い。
【００７７】
（４）第２の排除方法
図１０は、極小値検出を利用して弁や腱索の影響を排除する第２の排除方法の概念を示す図である。
【００７８】
第２の排除方法を用いた場合、ステップＳＴ４の処理では、Ｅ波のピーク値Ｅp2と、このピーク値Ｅp2よりも所定の比率だけ小さな中間値（例えばピーク値Ｅp2の半分のＥp2／２）との間に極小値が存在するか否かを判定する。ピーク値Ｅp2と中間値Ｅp2／２の間に極小値が存在しない場合は、血流以外の影響（弁や腱索の影響）は受けていないと判定し、ステップＳＴ５へ進む。
【００７９】
一方、ピーク値Ｅp2と中間値Ｅp2／２の間に極小値が存在する場合は、弁や腱索の影響を受けていると判定し、ステップＳＴ６へ進む。
【００８０】
極小値の存在の判定方法は特に限定するものではないが、例えば図１１に示したように、トレース波形の２次微分波形を求め、判定対象期間において２次微分波形が所定の閾値を超えるか否かによって判定することができる。
【００８１】
ステップＳＴ６では、図１０に示したように、実測トレース波形からピーク値Ｅp2と極小値との間の波形だけを抽出する。この期間の波形は弁や腱索の影響を受けていない波形であり、これを評価用減速波形とする。なお、第２の排除方法におけるステップＳＴ４やステップＳＴ６の処理も、評価波形生成部２３で行われる。
【００８２】
ステップＳＴ７では、この評価用減速波形（ピーク値Ｅp2と極小値との間の波形）に対して最小自乗法によって近似直線を求め、この近似直線からＥ波の減速時間ＤＣＴを算出する。
【００８３】
上述したように、第１の排除方法、或いは第２の排除方法のいずれの方法でも、Ｅ波の減速期間の波形に血流以外の影響（弁や腱索の影響）の有無を判定することが可能であり、影響のある場合にはこれを排除した評価用減速波形を生成することができる。
【００８４】
（５）簡略化最小自乗法
ステップＳＴ７では、生成した評価用減速波形から最小自乗法を用いて近似直線を求める演算を行うが、本実施形態では、通常の最小自乗法を簡略化して少ない演算量で近似直線を求める手法を採用している。
【００８５】
図１２は、従来から通常行われている最小自乗法によって近似直線を求める方法（図１２の左の図）と、本実施形態で用いている簡略化最小自乗法によって近似直線を求める方法（図１２の右の図）を示す図である。
【００８６】
前述したように、評価用減速波形から近似直線を求める場合、Ｅ波のピーク値Ｅpから減速期間の中間値Ｅp／２までのデータに対して最小自乗法を適用している。
【００８７】
従来の方法では、図１２の左図に示したように、求める近似直線を、
［数４］
ｙ ＝ ａ１*ｘ＋ａ0                            （式７）
としている。この場合、求める未知数は１次項の係数ａ１と定数項ａ0の２つとなり、次の連立方程式を解く必要がある。  
［数５］
ａ１*(Σxi)＋ａ0*Ｎ ＝ Σyi                   （式８）
ａ１*(Σxi2)＋ａ0*Σxi ＝ Σyi*xi              （式９）
ここで、xi, yi (i=1からＮ)は、ピーク値Ｅpから中間値Ｅp／２までのデータである。
【００８８】
これに対して、本実施形態では図１２右図に示したように、まず、Ｅ波のピーク位置が原点（０，０）となるように評価用減速波形を平行移動させる。この結果求める近似式（原点通過近似直線）は、
［数６］
ｙ ＝ ａ１*ｘ                                  （式１０）
となり、未知数は１次項の係数ａ１のみとなる。この係数ａ１は次のスカラー式から簡単に求めることができ、従来の方法に比べると積和演算の数は大幅に削減される。
【００８９】
［数７］
ａ１＝（Σyi*xi）／(Σxi2)                      （式１１）
ここで、（式１１）におけるxi, yi (i=1からＮ)は、評価用減速波形を平行移動させた後のピーク値Ｅpから中間値Ｅp／２までのデータである。求めた原点通過近似直線を、逆方向に平行移動させれば、従来の方法と全く同じ近似曲線が得られる。本実施形態では２回の平行移動を行っているが、平行移動に伴う演算負荷の増加はわずかなものであり、それよりも積和演算の低減による演算負荷の削減効果の方がはるかに大きい。
【００９０】
以上説明してきたように、本実施形態に係る超音波診断装置１及び診断パラメータ自動計測方法によれば、血流からのドップラスペクトラムに弁や腱索の影響が混入している場合であっても、Ｅ波の減速時間ＤＣＴを高い信頼性で自動計測することができる。
【００９１】
なお、本発明は上記の実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせても良い。
【図面の簡単な説明】
【００９２】
【図１】本発明の一実施形態に係る超音波診断装置の構成例を示すブロック図。
【図２】超音波診断装置の測定対象の一例を示す図。
【図３】左室流入血流のトレース波形の一例を示す図。
【図４】Ｅ波の減速時間ＤＣＴの自動計測処理の一例を示すフローチャート。
【図５】弁・腱索の影響を受けていないトレース波形と影響を受けているトレース波形を比較する図。
【図６】弁や腱索の影響をトレース波形から排除する第１の排除方法の概念説明図。
【図７】ＡＲＸモデルパラメータを求める方法の概念説明図。
【図８】ＡＲＸモデルパラメータを利用して弁や腱索の影響を受けている期間を判定する方法の概念説明図。
【図９】第１の排除方法における評価用減速波形の生成方法の概念説明図。
【図１０】弁や腱索の影響をトレース波形から排除する第２の排除方法の概念説明図。
【図１１】第２の排除方法において極小値を求める方法の概念説明図。
【図１２】本実施形態に係る簡略化最小自乗法の概念説明図。
【符号の説明】
【００９３】
１  超音波診断装置
２  超音波プローブ
８  ドップラ処理部
２１  トレース波形生成部
２２  ピーク波形検出部
２３  評価波形生成部
２４  減速時間算出部
２５  パラメータ記憶部
２６  心電波形入力部
【図１】
2010046284.tif 000003
【図２】
2010046284.tif 000004
【図３】
2010046284.tif 000005
【図４】
2010046284.tif 000006
【図５】
2010046284.tif 000007
【図６】
2010046284.tif 000008
【図７】
2010046284.tif 000009
【図８】
2010046284.tif 000010
【図９】
2010046284.tif 000011
【図１０】
2010046284.tif 000012
【図１１】
2010046284.tif 000013
【図１２】
2010046284.tif 000014