(57)【特許請求の範囲】
【請求項１】
スキャン領域へ超音波を送信し、その反射波を受信することによりエコー信号を得る超音波プローブと、
前記エコー信号からドプラ信号を抽出するドプラ信号抽出手段と、
前記ドプラ信号に基づいて、複数心拍に亘る第１のドプラ波形画像を生成する第1のドプラ波形画像生成手段と、
前記ドプラ信号に基づいて、前記複数心拍より少ない心拍数に亘る第2のドプラ波形画像を生成する第2のドプラ波形画像生成手段と、
前記ドプラ信号に基づいて、ドプラ信号における心拍の安定を示す情報を検出する安定検出手段と、
前記安定検出手段により心拍の安定を示す情報が検出されると、前記第1のドプラ波形画像と前記第２のドプラ波形画像とを並べて表示すると共に、前記安定を示す情報も付帯して表示する表示手段とを有することを特徴とする超音波診断装置。
【請求項２】
少なくとも前記第２のドプラ波形画像を記憶する記憶部を更に備え、
前記安定検出手段にて検出した情報により心拍が安定していることを示すと、前記超音波プローブによる前記超音波の送信を停止し、前記記憶部に記憶された前記第２のドプラ波形画像を用いて計測を行う計測モードへ移行する計測部とを更に備えることを特徴とする請求項１に記載の超音波診断装置。
【請求項３】
予め取得した基準ドプラ信号を記憶する記憶部と、
前記基準ドプラ信号と前記ドプラ信号との心拍毎の比較に基づいて、前記ドプラ信号における心拍の安定を示す情報を検出する安定検出手段とを更に備えたことを特徴とする請求項１に記載の超音波診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、超音波のドプラ波形を用いて診断を行う超音波診断装置に関する。
【背景技術】
【０００２】
超音波診断装置は、超音波プローブに内蔵された超音波振動子から超音波を被検体に対して送信し、被検体内で生じた反射波（以下、単にエコー信号と記載する）を超音波振動子で受信して、エコー信号に基づく超音波画像を生成する装置である。
【０００３】
近年、被検体の任意の位置における血流速度を計測する方法として、ドプラスペクトラム法が普及している。このドプラスペクトラム法では、被検体の同一部位に対して一定間隔で複数回の超音波送受波を行い、血球などの移動反射体において反射した超音波反射波からドプラ信号を検出する。そして、このドプラ信号にスペクトル解析を施すことにより、ドプラ波形（ドプラスペクトラムの画像データ）を生成する。このドプラ波形の時間変化から、例えば任意の位置の血流速度の最大流速や平均流速といった血流に関する情報を計測することができる（例えば、特許文献１を参照。）
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１０−２４０１９８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
先述したような超音波診断装置を用いてドプラ波形に基づいた計測作業を行う場合には、まず被検体に超音波を送受信するプローブを当接させた状態でドプラ信号の収集を行い、リアルタイムで更新されるドプラ波形をモニタ上に表示する。一般的にはこのとき、心臓が５〜６心拍分拍動する期間に収集されたドプラ波形が１画面中に表示される。プローブがリアルタイムでドプラ信号を収集するのに合わせて、画面中のドプラ波形も順次最新の波形に更新される。
【０００６】
ここで、ドプラ波形を用いた計測は、通常被験体の拍動が安定した状態で行われる。被検体の拍動の安定は、モニタに表示されるドプラ波形が周期的に変化することによって読み取れる。超音波診断装置の操作者は、モニタに表示されたドプラ波形が周期的に変化していることを確認すると、プローブによるドプラ信号の収集を停止させて、モニタ上のドプラ波形の更新を停止させた（フリーズ表示した）状態で計測を行う。
【０００７】
ここで、操作者による計測は、ドプラ波形上にマーカなどを配置してマーカ間の距離を測定することにより行われる。ここで、モニタ上には５〜６心拍などの多数の心拍を１画面中に表示させているために、１心拍あたりのドプラ波形の表示領域は小さいものとなってしまう。モニタ上に小さく表示されたドプラ波形上にマーカを配置しなければならないため、波形の細部の観察が困難であり、正しい位置にマーカを配置しばらつきの少ない計測を行うことは操作者にとって困難であった。一方モニタ上に５心拍などの多数の心拍が表示された状態から、画面を拡大表示する操作を行って１心拍を表示することも考えられるが、計測を行うたびに拡大表示する操作を行うことは操作者にとって操作負担が大きく、また複数心拍のうちどの心拍について表示を行っているか把握することが困難であるという課題もある。
【０００８】
そこで本願においては、ドプラ波形を用いた計測をより簡便に行うことを目的とする。
【課題を解決するための手段】
【０００９】
上記課題を解決するため実施形態の超音波診断装置においては、スキャン領域へ超音波を送信し、その反射波を受信することによりエコー信号を得る超音波プローブと、前記エコー信号からドプラ信号を抽出するドプラ信号抽出手段と、前記ドプラ信号に基づいて、複数心拍に亘る第１のドプラ波形画像を生成する第1のドプラ波形画像生成手段と、前記ドプラ信号に基づいて、前記複数心拍より少ない心拍数に亘る第2のドプラ波形画像を生成する第2のドプラ波形画像生成手段と、前記ドプラ信号に基づいて、ドプラ信号における心拍の安定を示す情報を検出する安定検出手段と、前記安定検出手段により心拍の安定を示す情報が検出されると、前記第1のドプラ波形画像と前記第２のドプラ波形画像とを並べて表示すると共に、前記安定を示す情報も付帯して表示する表示手段とを有する。
【図面の簡単な説明】
【００１１】
【図１】実施例に係る超音波診断装置の内部構成を示すブロック図。
【図２】実施例に係るドプラ波形及び超音波画像の表示例を示す図。
【図３】実施例に係るドプラ波形の処理に関連する各構成要素についての、機能的な接続態様を示したブロック図。
【図４】実施例に係るリングバッファの構成について示す図。
【図５】実施例に係る拡大ドプラ波形を表示する時間領域の表示例を示す図。
【図６】実施例に係る平均ドプラ波形を拡大ドプラ波形と重畳表示した際の表示例を示す図。
【図７】実施例に係る心拍の安定を示す情報の表示例を示す図。
【図８】実施例に係る、拡大ドプラ波形を表示する処理に関連する処理の流れを示すフローチャート。
【発明を実施するための形態】
【００１２】
（超音波診断装置１の構成）
以下、本開示の実施例について図面を参照して説明する。図１は、各実施例に係る超音波診断装置１の構成を示したブロック図である。超音波診断装置１は、システム制御部１０、入力部２１、表示部２３、超音波プローブ３０、及び生体信号検出ユニット４０を組み合わせることにより構成される。なお、本発明における超音波診断装置１の構成はこれに限られるものではなく、適宣構成要素を追加しても構わないし、超音波診断装置１の構成要素の幾つかを省略して、超音波診断装置１の外部に接続した装置が省略された構成要素の役割を果たすものであっても構わない。例えば後述するドプラ波形処理部１７などのドプラ波形の生成や処理に用いる構成要素を超音波診断装置１の外部接続した外部処理装置に行われる構成としても構わない。あるいは本実施例に述べるドプラ波形の拡大表示機能やドプラ波形処理機能を実行するためのプログラムを任意のコンピュータにインストールして実行することで、コンピュータを実質的に超音波診断装置１として動作させるものであっても構わない。
【００１３】
超音波プローブ３０はシステム制御部１０に接続される、超音波を送受波するための部材である。超音波プローブ３０は超音波振動子３１を内蔵し、送受信部１３から出力された駆動信号に基づいて超音波振動子３１を駆動する。
【００１４】
超音波振動子３１は、送受信部１３を介して入力された駆動信号に基づいて超音波を発生し、また被検体から反射された超音波を受信して電気信号へと変換する部材である。超音波振動子３１は超音波を受信することにより生じた電気信号（以下、単にエコー信号と記載する）を送受信部１３へと出力する。なお、本実施例における超音波振動子３１は、リニア形、コンベックス形、セクタ形などの２次元領域をスキャンしうる種々の形状によって構成されていても構わない。
【００１５】
システム制御部１０は、記憶部１２、超音波画像生成部１５、ドプラ波形処理部１７、Bモード処理部１８、ドプラ処理部１９から構成される。システム制御部１０はシステム制御部１０に接続された各構成要素へ制御命令を出力し、超音波診断装置１を統括的に制御する。
【００１６】
送受信部１３は、超音波プローブ３０が受信したエコー信号を処理するアンプ回路、Ａ／Ｄ変換器、加算器などを備える。アンプ回路は、超音波プローブ３０が受信したエコー信号を増幅し、Ａ／Ｄ変換器へ出力する。Ａ／Ｄ変換器は、増幅されたエコー信号の受信指向性を決定するために必要な遅延時間をエコー信号に与え、加算器へと出力する。加算器は遅延時間を与えられたエコー信号を加算することで、超音波を送信するスキャンラインに対応したエコー信号を得る。送受信部２１はスキャンラインに対応したエコー信号をＢモード処理部１８あるいはドプラ処理部１９へと出力する。また送受信部１３は、超音波振動子３１を駆動し超音波を送信するための駆動信号を超音波超音波プローブ３０へと出力する。
【００１７】
Ｂモード処理部１８は、送受信部１３が出力したエコー信号の振幅強度に応じて変化するＢモード信号を生成する。Ｂモード処理部１８は生成したＢモード信号を記憶部１２へ出力する。
【００１８】
ドプラ処理部１９は、エコー信号に対し直交位相検波処理を施すことにより、エコー信号の周波数遷移を示すドプラ信号を生成する。このドプラ信号の収集位置は、超音波画像中に表示したサンプルゲートの位置を調節することによって設定される。図２にサンプルゲートの表示例を示す。超音波診断装置１の操作者は入力部２１を操作して、超音波画像１００中のドプラ信号を収集した領域にサンプルゲート１０２を位置させる。送受信部１３はサンプルゲート１０２を通る走査線１０１に対して超音波の送信を行い、ドプラ処理部１９へと出力する。ドプラ処理部１９はドプラ信号を生成すると、これを記憶部１２へ出力する。
【００１９】
超音波画像生成部１５は、Ｂモード処理部１８から出力されたＢモード信号、あるいはドプラ処理部１９から出力されたドプラ信号を記憶部１２から読み出して、超音波画像を生成する。超音波画像生成部１５がＢモード信号を用いて超音波画像（以下、単にＢモード画像と記載する）を生成する場合は、超音波画像生成部１５はＢモード信号の振幅値から算出された画素値を超音波の送受信方向に対応した２次元平面へマッピングすることでＢモード画像を生成する。一方超音波画像生成部１５がドプラ信号に基づいて超音波画像（以下、単にドプラ画像と記載する）を生成する場合には、超音波画像生成部１５はドプラ処理部１９から受けたドプラ信号にフィルタ処理を施すことにより、被検体中の移動していない組織から得られたクラッタ成分を除去し、血流によって得られたドプラ信号のみを抽出する。そして、抽出したドプラ信号に自己相関処理を施して、血流の平均流速値や分散値を色情報に変換したカラードプラ信号を得る。そして、超音波画像生成部１５はこのカラードプラ信号を超音波の送受信方向に対応した２次元平面へマッピングすることで、ドプラ画像を生成する。
【００２０】
ドプラ波形処理部１７は、ドプラ処理部１９から出力された所定の期間中のドプラ信号を記憶部１２から読み出して、ドプラ信号の時間的な変化を示すドプラ波形を生成する。このドプラ波形は、縦軸を流速として、横軸を時間軸とした２次元波形である。ドプラ波形処理部１７は記憶部１２からドプラ信号を循環的に読み出してドプラ波形の生成処理を行う。このドプラ波形の生成処理については後に詳しく述べる。
【００２１】
ドプラ波形処理部１７は、所定の期間に亘るドプラ信号の変化を標準ドプラ波形として生成し表示部２３へと出力する。標準ドプラ波形には、例えば被検体の約６心拍分のドプラ信号の変化が収まるような、比較的長い期間のドプラ信号の変化が表される。ドプラ波形処理部１７は標準ドプラ波形を生成する一方で、標準ドプラ波形から１心拍分のドプラ波形の変化を切り出した拡大ドプラ波形を生成して表示部２３へと出力する。なお、これ以降標準ドプラ波形と拡大ドプラ波形とを総称する場合には、単にドプラ波形とだけ記載する。またドプラ波形処理部１７は、ドプラ波形の形状をより分かりやすく表示するために、各時点におけるドプラ信号の最大値を抽出して最大値波形（以下、オートトレース波形と記載する）を出力するオートトレース機能を備える。表示部２３はドプラ波形上に、オートトレース波形を示す曲線２１１を表示する。
【００２２】
またドプラ波形処理部１７は、ドプラ処理部１９から出力されたドプラ信号と後述する生体信号検出ユニット４０から出力される生体信号とを組み合わせて、ドプラ信号に対する計測処理を行う。具体的には、ドプラ波形処理部１７は生体信号をもとにドプラ波形中の１心拍の期間、拡張期の時相、収縮期の時相を設定し、１心拍期間中のドプラ信号を用いて計測処理を行い、ドプラ計測値として最高流速値（Ｖｍａｘ）、平均流速値（Ｖｍｅａｎ）、抵抗指数（ＲＩ値）、拍動指数（ＰＩ値）を算出する。なお、最高流速値（Ｖｍａｘ）は、１心拍期間中に得られた最大のドプラ信号値を示す。平均流速値（Ｖｍｅａｎ）は、１心拍期間中の平均のドプラ信号値を示す。抵抗指数（ＲＩ値）は、（収縮期流速−拡張期流速）／収縮期流速を計算することにより求められる。拍動指数（ＰＩ値）は、（収縮期流速−拡張期流速）／平均流速を計算することにより求められる。なお、ドプラ波形処理部１７が算出するドプラ計測値の種類は、ここに挙げたものに限定されない。
【００２３】
またドプラ波形処理部１７は、過去の心拍におけるドプラ計測値と拡大ドプラ波形２１０におけるドプラ計測値とを比較して、最新の心拍が安定しているか否かを判断する安定検出機能を備える。心拍の安定検出方法については後に詳しく述べる。
【００２４】
記憶部１２は、例えばＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）やフラッシュメモリ、及びＨＤＤ（Ｈａｒｄ  Ｄｉｓｋ  Ｄｒｉｖｅ）などから構成される記憶媒体である。記憶部１２は超音波画像生成部１６が生成した超音波画像や、あるいはドプラ処理部１９が出力したドプラ信号、あるいはドプラ波形処理部１７が出力したドプラ計測値などを記憶する。また、記憶部１２は過去の検査データとして、過去に測定して得た被検体のドプラ波形のデータやドプラ計測値、あるいは標準の検査データとして、健常者のドプラ波形のデータやドプラ計測値を記憶する。
【００２５】
入力部２１は、例えば機械的なボタン、ダイヤル、トラックボール、ジョイスティック、スライダやホイールなどの種々の操作デバイスを用いて構成され、オペレータの入力操作を受け付けて、入力操作を電気信号に変換してシステム制御部１０へ出力する部材である。入力部２１は受け付けた入力操作に応じて、例えば送受信部１３にスキャンの開始・停止を指示する指示信号や、表示部２３に表示される標準ドプラ波形と拡大ドプラ波形の更新を停止させて計測モードへ移行するフリーズ信号をシステム制御部１０へ出力する。
【００２６】
表示部２３は、例えばＬＣＤ（Ｌｕｃｉｄ  Ｃｒｙｓｔａｌ  Ｄｉｓｐｌａｙ）ディスプレイや有機ＥＬ（Ｅｌｅｃｔｒｏ  ｌｕｍｉｎｅｓｃｅｎｃｅ）ディスプレイ、あるいはブラウン管ディスプレイなどの任意のディスプレイにより構成される。表示部２３は、参照画像やライブ像を表示する。あるいは、送受信部１３がスキャンを行う際のスキャンパラメータや、超音波画像生成部１５がライブ像を生成する際の画像化パラメータ、参照画像生成部１３が参照ボリュームデータから参照画像を生成する際の画像化パラメータ、種々のメッセージなどを表示する。
【００２７】
生体信号検出ユニット４０は、被検体に取り付けられ、被検体から定期的に発せられる生体信号を電気信号に変換して制御部１０へ出力する部材である。ドプラ波形処理部１７は、生体信号検出ユニット４０から出力される電気信号に基づいて被検体の心時相を検出して、ドプラ波形の生成や計測処理を行う。例えば生体信号検出ユニット４０は心電計によって構成され、被検体に取り付けた電極から心拍に同期して発せられる電気信号を検出して、ＥＣＧ（心電波形）として制御部１０へ出力する。また例えば生体信号検出ユニット４０は心音計によって構成され、被検体に取り付けた小型のスピーカから心拍に同期して発せられる心音信号を検出して、ＰＣＧ（心音波形）として制御部１０へ出力する。出力された生体信号は、図２中の生体信号波形２１２に示されるようにドプラ波形と並べて表示される。以降の実施例では、例として生体信号検出ユニット４０が心電計によって構成されるものとして述べる。なお、ドプラ波形処理部１７がドプラ信号の変化に基づいて心時相を検出する場合には、生体信号検出ユニット４０の構成及び生体信号波形２１２を省略して超音波診断装置１を構成しても構わない。
【００２８】
（拡大ドプラ波形の表示処理）
先に述べた各構成要素に基づいて、超音波診断装置１は標準ドプラ波形２００及び拡大ドプラ波形２１０の生成及び表示を行う。図２は、表示部２３に表示されるＢモード画像１００、標準ドプラ波形２００、及び拡大ドプラ波形２１０の表示例を示した図である。例えば表示部２３の右上の領域には、サンプルゲート１０２を設定するためにＢモード画像１００が表示される。操作者は入力部２１を用いてＢモード画像１００に表示されたサンプルゲート１０２の表示位置を変化させて、ドプラ信号を収集するための領域を設定する。Ｂモード画像１００上には、サンプルゲート１０２と合わせてサンプルゲート１０２を通過する走査線１０１の位置が表示される。
【００２９】
例えばＢモード画像１００の下には、ドプラ波形処理部１７が生成した標準ドプラ波形２００が表示される。標準ドプラ波形２００は先述したように、例えば被検体の約６心拍分を表示するような、長い時間領域のドプラ信号の変化が表示される。ここで、標準ドプラ波形２００と並べて心電波形を表示しても構わない。また、標準ドプラ波形２００上には、現在の時相の位置を示すムービングバー２０１が表示される。また、標準ドプラ波形２００上には、拡大ドプラ波形２１０が対応する領域を示す枠線２０２が表示される。
【００３０】
例えばＢモード画像１００の左には、ドプラ波形処理部１７が生成した拡大ドプラ波形２１０が表示される。拡大ドプラ波形２１０は、標準ドプラ波形２００のうち最新の１心拍分を拡大して表示した波形である。拡大ドプラ波形２１０上には、オートトレース波形２１１が重畳して表示される。なお、標準ドプラ波形２００と並べて心電波形を表示している場合には、拡大ドプラ波形２１０と並べて心電波形を表示しても構わないし、更に心電波形上に心電波形の平均値を示す心電オートトレース波形２１２を重畳して表示しても構わない。また、表示部２３は、例えば拡大ドプラ波形２１０の心拍に関するドプラ計測値を計測値表示領域２２０に表示する。
【００３１】
ここで図２に示されるように、拡大ドプラ波形２１０は、１心拍分あたりの表示領域が標準ドプラ波形２００よりも大きくなるように大きな表示領域に表示される。これにより、操作者は被検体の最新の心拍がどのように変化しているかを大きな表示領域で視認し、心拍が安定しているか否かをより簡便に判断することができる。また、ドプラ波形の表示をフリーズし、ドプラ波形上に手動でマーカを配置する計測操作を行う際にも、大きく表示された拡大ドプラ波形２１０を用いてより正確にマーカを配置することができる。また、大きな表示領域でドプラ波形を視認する際に操作者の拡大操作などが不要であるため、操作者の操作負担を軽減することができる。
【００３２】
なお、図２において示した各画像やドプラ波形の表示位置あるいは表示方法は単なる例であり、表示位置あるいは表示方法はここに挙げたものに限定されない。例えばＢモード画像１００の表示を省略して標準ドプラ波形２００と拡大ドプラ波形２１０を並べて表示しても構わないし、標準ドプラ波形２００の表示を省略してＢモード画像１００と拡大ドプラ波形２１０とを並べて表示しても構わない。あるいは、各画像やドプラ波形の表示位置を並べ替えても構わない。
【００３３】
（リングバッファを用いたドプラ波形の生成処理）
図３は、ドプラ波形の処理に関連する各構成要素についての、機能的な接続態様を示したブロック図である。ここで、記憶部１２は複数の記憶領域を繋げて構成したリングバッファ１２１によって構成される。ここでは、リングバッファ１２を用いたドプラ波形の生成処理及びドプラ計測値の算出処理について述べる。
【００３４】
まずドプラ処理部１９は、送受信部１３から出力されたエコー信号にドプラ処理を施し、ドプラ信号として記憶部１２へと出力する。リングバッファ１２１は、ドプラ信号が出力された時刻に合わせて記憶領域を順に変化させながら記憶を行う。リングバッファ１２１は、例えば時刻ｔ０〜ｔ１の期間中に出力されたドプラ信号を領域０に、時刻ｔ１〜ｔ２の期間中に出力されたドプラ信号を領域１に記憶する。リングバッファ１２は時刻が進むのに合わせて記憶領域を変化させながらドプラ信号の記憶を続け、ｔｎ−１〜ｔｎの期間中に出力された領域を最後の領域である領域ｎに記憶すると、次のｔｎ〜ｔｎ＋１の期間中に出力されたドプラ信号を、最初の領域０へ上書きする。更に次のｔｎ＋１〜ｔｎ＋２の期間中に出力されたドプラ信号は、次の領域１へ上書きされる。このように、リングバッファ１２はドプラ信号の記憶領域を時刻の変化と共に変化させながら記憶を行い、すべての記憶領域に記憶を行うと最初の記憶領域から順にドプラ信号を上書きしてゆく。
【００３５】
一方ドプラ波形処理部１７が標準ドプラ波形２００を生成する場合には、記憶領域を順に変化させながらドプラ信号を読み出して、波形生成処理を行う。まずドプラ波形処理部１７は領域０に記憶されたドプラ信号を読み出してこれを画像情報へ変換する。次にドプラ波形処理部１７は領域１に記憶されたドプラ信号を読み出して画像情報へ変換し、これを領域０における画像情報に付け加える。ドプラ波形処理部１７は記憶領域を変化させながら画像情報の変換処理を続け、ｍ個（ここで、ｍはｍ＜ｎを満たす数）のドプラ信号による画像情報を繋ぎ合わせることで、これを標準ドプラ波形２００として表示部２３へ出力する。そして、ドプラ信号を読み出す領域が最後の領域である領域ｎに到達すると、次は領域０に戻ってドプラ信号の読み出し及び画像情報の変換処理を続ける。ドプラ処理部１９によるリングバッファ１２への書き込み及び、ドプラ波形処理部１７によるドプラ信号の読み出しが循環的に行われることによって、ドプラ処理部１９によるドプラ信号の書き込みとドプラ波形処理部１７による読み出しと画像情報の生成処理がリアルタイムで行われる。
【００３６】
なお、リングバッファ１２の各領域へはドプラ信号と合わせて生体信号検出ユニット４０から検出された時相情報が合わせて記憶される。図４は、各領域に書き込まれるドプラ信号及び時相情報の様子を模式的に表した図である。生体信号検出ユニット４０は、心電波形上で特徴的な波形であるＲ波を検出すると、これをＲ波タグデータとしてドプラ処理部１９へと出力する。ドプラ処理部１９はＲ波タグデータを受け取ると、ドプラ信号と合わせてこれをリングバッファ１２へと出力する。図４に示すように、リングバッファ１２はドプラ信号を時系列順に記憶するとともに、Ｒ波タグデータが得られた時刻のドプラ信号と同じ領域にＲ波タグデータを記憶する。
【００３７】
ドプラ波形処理部１７は、リングバッファ１２に記憶されたＲ波タグデータを利用して、拡大ドプラ波形２１０中のドプラ計測値を算出するためのドプラ信号読み出し領域を設定する。即ち、ドプラ波形処理部１７はＲ波タグデータが記憶された領域からｐ個の領域分離れた領域を始点として、約１心拍分の領域数であるｑ個（ここで、ｑはｑ＜ｍを満たす数）のドプラ信号を抽出し、このｑ個のドプラ信号を用いてドプラ計測値を算出する。図４を用いて具体的に説明すると、ドプラ波形処理部１７は図４に示すように、Ｒ波タグデータが得られた領域から２領域分離れた領域を始点とした、Ｘ（ｋ）に示す６個分のドプラ信号を抽出する。そして、ドプラ波形処理部１７は抽出した６個のドプラ信号を基にドプラ計測値を算出する。リングバッファ１２上の各領域のドプラ信号はドプラ信号の収集に合わせて次々と上書きされてゆくため、ドプラ波形処理部１７は上書きされるドプラ信号に合わせてドプラ信号読み出し領域を変化させながらドプラ計測値の算出を行う。過去の時点において算出されたドプラ計測値（図４中でＸ（ｋ−１）で示される領域から得られたドプラ計測値）は、過去のドプラ計測値として記憶部１２に記憶される。
【００３８】
図５（ａ）は、拡大ドプラ波形２１０を表示する時間領域を標準ドプラ波形２００上に示した図である。心電波形上でＲ波が得られた時刻をＴｅｃｇとすると、ドプラ波形処理部１７はＴｅｃｇから得られた時刻からｄＴｏｆｆｃｅｔだけ離れた時間領域（即ち、リングバッファ１２上でｐ個分のドプラ信号を記憶する時間領域）を始点として、ｄＴｃｕｔの時間領域（即ち、リングバッファ１２上でｑ個分のドプラ信号を記憶する時間領域）内にあるドプラ信号を抽出する。
【００３９】
なお、ドプラ波形処理部１７はＲ波タグデータを用いる代わりに、オートトレース波形を用いてドプラ計測値を算出するためのドプラ信号読み出し領域を設定しても構わない。この場合、リングバッファ１２上にはオートトレース波形で最大値を検出した時相と対応する領域に、最大値タグデータが書き込まれる。最大値タグデータが得られた時刻をＴｍａｘとすると、ドプラ波形処理部１７はＴｍａｘが得られた時刻からｄＴｏｆｆｃｅｔだけ離れた時間領域を始点として、ｄＴｃｕｔの時間領域内にあるドプラ信号を抽出する。
【００４０】
ｄＴｏｆｆｃｅｔ及びｄＴｃｕｔの長さは、例えば入力部２１を用いて操作者が任意に設定するものであっても構わないし、例えば心電波形やオートトレース波形から検出された心拍の時間間隔に基づいて制御部１０が設定するものであっても構わない。
【００４１】
（拡大ドプラ波形の平均化機能）
ドプラ波形処理部１７は先述した拡大ドプラ波形２１０を表示する際に、最新の一心拍を拡大して表示する表示モードの他に、過去の数心拍を平均化した平均ドプラ波形２１１を拡大ドプラ波形２１０の表示領域に重畳して表示する機能を備える。例えば入力部２１の操作によって平均ドプラ波形２１１を表示することが指示される。すると、ドプラ波形処理部１７は過去の例えば５心拍に亘るドプラデータを平均化して、平均ドプラ波形２１１を生成する。具体的には、ドプラ波形処理部１７は平均ドプラ波形２１１を以下の数式（数１）によって生成する。
【００４２】
（数１）    Ｙ（Ｎ）  ＝ＳＵＭ（ａiＸ（Ｎ−ｉ））  ｉ＝０  ｔｏ  Ｎ  
上述の（数１）において、ＳＵＭ（）は括弧内の合計値を、Ｙ(Ｎ)は平均ドプラ波形２１１のドプラ値を、Ｎは平均化を行う心拍の数を、Ｘ（ｉ）は各心拍におけるドプラ値を、aiは平均化演算を行う際の各心拍に対する重みづけ係数をそれぞれ示している。例えば５心拍分のドプラデータを重み付けなしで平均化する際には、Ｎ＝５、a1〜a5＝１として（数１）の演算を行えばよい。平均化を行う心拍の数や重みづけ係数はここにあげたものに限られず、例えば２心拍について平均化を行っても構わないし、あるいは重みづけ係数を最新の心拍に近いドプラ波形ほど大きくなるように設定して、平均ドプラ波形２１１の表示がより最新の心拍の変化を反映するようにしても構わない。
【００４３】
図６に、平均ドプラ波形２１１を拡大ドプラ波形２１０と重畳表示した際の画面表示例を示す。図２に示した拡大ドプラ波形２１０の表示領域内に、平均ドプラ波形２１１が点線で重畳表示される。ここで、平均ドプラ波形２１１の表示はオートトレース波形として示しているが、単にオートトレース処理を行う前のドプラ値を平均化した波形を表示するものであっても構わない。この場合、拡大ドプラ波形２１０と平均ドプラ波形２１１の表示領域が重なるために、表示部２３は拡大ドプラ波形２１０と平均ドプラ波形２１１の表示色を変化させて、両者を区別して表示する。また、平均ドプラ波形２１１のオートトレース波形を表示する際には、過去の心拍のオートトレースによって得られた値をＸ（ｎ）として（数１）に述べた演算を行っても構わない。これにより、より少ない演算数で平均ドプラ波形２１１のオートトレース波形を算出することができる。また、ドプラ波形だけでなく心電波形を平均化した平均心電波形を更に表示しても構わない。
【００４４】
また、図２に示したときと同様に、Ｂモード画像１００の下には標準ドプラ波形２００が表示される。標準ドプラ波形２００内には、拡大ドプラ波形２１０および平均ドプラ波形２１１が対応する領域を示す枠線２０２が表示される。
【００４５】
この平均ドプラ波形２１１を表示する機能により、操作者は最新の心拍の様子を表す拡大ドプラ波形２１０と、過去の心拍変化の概要を表す平均ドプラ波形２１１とを同じ表示領域に重畳表示して、両者の違いを見比べながら診断を行うことができる。
【００４６】
（ドプラ波形の安定検出処理）
ドプラ波形処理部１７は先述した拡大ドプラ波形２１０を表示する機能の他に、被検体の心拍が安定しているか否かを検出する安定検出機能を更に備える。安定検出機能は、例えば被検体の最新の心拍におけるドプラ計測値と、過去の数心拍において得られたドプラ計測値の平均値とを比較することによって、最新の心拍が過去の心拍に比べて大きく変化していないか、即ち心拍が安定しているかを検出する機能である。
【００４７】
例えば入力部２１の操作によってドプラ波形処理部１７の安定検出機能を動作させるよう指示を受けると、ドプラ波形処理部１７は過去の数心拍に亘るドプラ計測値と、最新の心拍におけるドプラ計測値との差を比較する。このとき比較に用いるドプラ計測値は、先にドプラ計測値の例として挙げた最高流速値（Ｖｍａｘ）、平均流速値（Ｖｍｅａｎ）、抵抗指数（ＲＩ値）、拍動指数（ＰＩ値）のいずれかを用いる。いずれを用いるかは例えば入力部２１の操作によって定められるものであっても構わないし、複数の計測値を用いて比較を行っても構わない。
【００４８】
ドプラ波形処理部１７は、過去数心拍のドプラ計測値と最新の１心拍におけるドプラ計測値との標準偏差を比較し、偏差の値が予め定めた閾値以下の値に収まった場合には心拍が安定しているものと判断する。なお、心拍の安定を判断する閾値はここに挙げた値の他に、任意に設定されるものであって構わない。例えば過去数心拍のドプラ計測値の平均値と、最新の１心拍におけるドプラ計測値との値の差を求め、この差がドプラ計測値の平均値に対して２０％以下であった場合に心拍の安定を判断しても構わない。
【００４９】
また、ドプラ波形処理部１７はリアルタイムでドプラ値を取得して、最新の１心拍のドプラ計測値を算出するたびに上述したドプラ計測値の比較を行う。そして、最新の１心拍のドプラ計測値が安定していると連続で５回判断されたときに、現在の心拍が安定していると判断する。なお、心拍の安定を判断する回数は５回に限られず、入力部２１の操作によって任意の値に変更するものであっても構わない。
【００５０】
ドプラ波形処理部１７は現在の心拍が安定していると判断すると、表示部２３に心拍の安定を示す情報を表示させる。図７に心拍の安定を示す情報の表示例を示す。図７（ａ）は心拍が安定していない時の拡大ドプラ波形２１０の表示領域を、図７（ｂ）は心拍が安定した時の拡大ドプラ波形２１０の表示領域をそれぞれ示している。ドプラ波形処理部１７は表示部２３を用いて、図７（ｂ）に示すように計測値表示ウインドウ２２０の表示を白黒反転させ、また拡大ドプラ波形２１０の表示領域下部に心拍が安定したことを示すメッセージを表示させる。更に、拡大ドプラ波形２１０の表示領域を示す枠線を例えば太くしたり、色を変えたりして強調表示する。なお、心拍が安定したことを示す方法はここに挙げた動作に限られず、例えばオートトレース波形の色を変化させたり、あるいは音声やビープ音などを発することで心拍の安定を操作者に示しても構わない。
【００５１】
この動作によれば、操作者は過去の心拍におけるドプラ計測値と最新の心拍におけるドプラ計測値の比較に基づいて、被検体の心拍の安定を検出することができる。これにより、波形の目視に基づいて心拍の安定を判断する従来の手法に比べて、より定量的な診断を行うことが可能となる。
【００５２】
なお、本実施形態においては平均ドプラ計測値と最新の心拍におけるドプラ計測値とを比較すると述べたが、動作はこれに限られるものではない。例えばドプラ計測値の代わりに、過去の心拍におけるドプラ値と最新の心拍におけるドプラ値とを比較して心拍の安定を検出しても構わないし、あるいはドプラ計測値の平均化を行なわず、過去の心拍それぞれで得られたドプラ計測値と最新の心拍のドプラ計測値とをそれぞれ比較しても構わない。
【００５３】
またあるいは、ドプラ波形処理部１７が現在の波形が安定していることを検出すると、自動的に超音波プローブ３０による超音波の送受信を停止し、ドプラ波形の更新を停止するフリーズ動作を行い、更に標準ドプラ波形２００や拡大ドプラ波形２１０中にマーカを配置して計測値を得る計測モードへの移行を自動的に行っても構わない。これにより、操作者は心拍が安定した時に逐一フリーズ操作と計測モードへの移行動作を行う必要がなく、操作者の操作負担を軽減した迅速な診断を行うことができる。
【００５４】
またあるいは、記憶部１２に過去の検査で得られた被検体のドプラ計測値や、あるいは検査を行う被検体と同一の年代、性別の健常者の被検体を検査して得られたドプラ計測値を比較対象として記憶しておく。ドプラ波形処理部１７は、記憶部１２に記憶されたドプラ計測値と、最新のドプラ計測値とを比較して現在の心拍が異常値であるか否かを検出しても構わない。先述した実施例においては心拍が安定した場合に心拍が安定したことを示すメッセージを表示させていたが、例えば健常者のドプラ計測値との比較を行う場合、健常者のドプラ計測値と最新のドプラ計測値との差が大きい場合には現在の心拍が異常な値を示していることを示す警告メッセージを表示させても構わない。このとき、過去のドプラ計測値や健常者のドプラ計測値は記憶部１２に記憶される場合に限られず、例えば超音波診断装置１の外部に接続された記憶装置から読みだすものであっても構わないし、あるいは超音波診断装置１にネットワーク接続された他の端末からドプラ計測値を取得するものであっても構わない。
【００５５】
（ドプラ波形表示処理の流れ）
図８は、超音波診断装置１が拡大ドプラ波形２１０を表示した上で、心拍安定検出機能に基づいて安定を検出して計測モードへ移行する際の、一連の処理の流れを示したフローチャートである。以下、図８に沿って処理の流れについて述べる。
【００５６】
まず、超音波診断装置１の使用者がドプラ波形の計測を開始する（ステップ１０００）と、超音波プローブ３０は超音波の送受信を行う。そして、ドプラ処理部１９が出力したドプラデータに基づいて、ドプラ波形処理部１７が標準ドプラ波形２００を生成する（ステップ１００１）。ドプラ波形処理部１７は標準ドプラ波形２００を生成すると、ドプラ値を用いてオートトレース処理を行う（ステップ１００２）。ドプラ波形処理部１７は、生成した標準ドプラ波形２００及びオートトレース波形を表示部２３へ表示させる。
【００５７】
ドプラ波形処理部１７が標準ドプラ波形２００を生成しオートトレース処理を終えると、ドプラ波形処理部１７は最新の１心拍分のドプラデータを記憶部１２から抽出し、抽出したドプラデータから拡大ドプラ波形２１０を生成する（ステップ１００３）。ドプラ波形処理部１７は、生成した拡大ドプラ波形２１０を表示部２３へ表示させる。更に、ドプラ波形処理部１７は抽出したドプラデータに基づいて、最新の１心拍に関するドプラ計測値を算出する（ステップ１００４）。ドプラ波形処理部１７は、算出したドプラ計測値を心拍毎に区別して記憶部１２に記憶する。ここまでの動作により、表示部２３には標準ドプラ波形２００とは別に、最新の１心拍を表す拡大ドプラ波形２１０と拡大ドプラ波形２１０のドプラ計測値が表示部２３中に表示される。
【００５８】
次に、ドプラ波形処理部１７は記憶部１２に記憶された過去の心拍におけるドプラ計測値と、最新の１心拍におけるドプラ計測値との差を比較する。そして、ドプラ計測値同士の差異が所定の閾値以下に収まっているか否かを判断する。ドプラ波形処理部１７は差異が連続でＮ回以上閾値以下に収まっていると判断すると（ステップ１００５のＹｅｓ）、拡大ドプラ波形２１０の表示領域内に心拍が安定したことを示すメッセージを表示させる（ステップ１００６）。一方ドプラ波形処理部１７が、ドプラ計測値同士の差異が連続して閾値以下に収まっている回数がＮ回以下であると判断すると（ステップ１００５のＮｏ）、ドプラ波形処理部１７はステップ１００１に戻ってドプラ波形の生成を繰り返す。
【００５９】
表示部２３がステップ１００６にて心拍が安定したことを示すメッセージを表示させると、ドプラ波形処理部１７は記憶部１２に記憶された過去の被検体のドプラ計測値、あるいは健常者のドプラ計測値と最新の１心拍におけるドプラ計測値との差異を比較する（ステップ１００７）。そして、ドプラ波形処理部１７は差異が所定の閾値以上であると判断すると（ステップ１００７のＮｏ）、拡大ドプラ波形２１０の表示領域に現在の心拍におけるドプラ計測値が異常値であることを示すメッセージを表示させる（ステップ１００８）。
【００６０】
ドプラ波形処理部１７が、差異が所定の閾値以下であると判断する（ステップ１００７のＹｅｓ）か、あるいはドプラ計測値が異常値であることを示すメッセージを表示させる（ステップ１００８）と、超音波プローブ３０は超音波の送受信を停止する。更に、ドプラ波形処理部１７は標準ドプラ波形２００及び拡大ドプラ波形２１０の生成を停止し、表示部２３に表示させる標準ドプラ波形２００及び拡大ドプラ波形２１０を停止させる、いわゆるフリーズ動作を行う（ステップ１００９）。ドプラ波形処理部１７がフリーズ動作を行うと、制御部１０は表示部２３にマーカを表示させて、入力部２１の操作によってマーカを移動させてマーカの位置に基づく計測値を表示させる計測モードへと移行して（ステップ１０１０）、処理を終了する（ステップ１０１１）。
【００６１】
なお、図８に述べた処理ではステップ１００７に述べた過去のドプラ計測値あるいは健常者のドプラ計測値との比較を心拍安定後に行うものとして述べた。しかし、動作はこれに限られるものではなく、例えばステップ１００４で最新の１心拍の計測値を算出するたびごとに過去のドプラ計測値あるいは健常者のドプラ計測値との比較を行っても構わない。
【００６２】
以上に述べた処理により、超音波診断装置１は１心拍分あたりの表示領域を標準ドプラ波形２００よりも大きくした拡大ドプラ波形２１０を表示する。これにより、操作者は被検体の最新の心拍がどのように変化しているかを大きな表示領域で視認し、心拍が安定しているか否かをより簡便に判断することができる。また、ドプラ波形の表示をフリーズし、ドプラ波形上に手動でマーカを配置する計測操作を行う際にも、大きく表示された拡大ドプラ波形２１０を用いてより正確にマーカを配置することができる。また、大きな表示領域でドプラ波形を視認する際に操作者の拡大操作などが不要であるため、操作者の操作負担を軽減することができる。
【００６３】
また以上に述べた処理により、超音波診断装置１は拡大ドプラ波形２１０に重畳して、過去の心拍におけるドプラ波形を平均化した平均ドプラ波形２１１を表示する。これにより、操作者は拡大ドプラ波形２１０と平均ドプラ波形２１１とを比較して、現在の心拍が過去の心拍と比較して安定しているかを確認し、また心拍のどの時相上で差異が大きいかを容易に確認することができる。
【００６４】
また以上に述べた処理により、超音波診断装置１は現在の心拍が安定しているか否かを自動で検出し、検出結果を表示部２３に表示する。心拍の安定を機械的に検出することができるため、操作者の標準ドプラ波形２００の目視による安定検出に比べてより定量的な判断を行うことができる。また、操作者は標準ドプラ波形２００を観察し続ける必要がないため、診断に係る操作者の負担をより軽減することができる。
【００６５】
また以上に述べた処理により、超音波診断装置１は被検体の過去のドプラ計測値や健常者のドプラ計測値を比較して、現在の心拍が異常であるか否かを自動で検出し、検出結果を表示部２３に表示する。心拍の異常を機械的に検出することができ、操作者が心拍の異常を判断し損ねる事態を防止することができる。
【００６６】
また以上に述べた処理により、超音波診断装置１は心拍が安定したことを検出すると、フリーズ操作及び計測モードへの移行を自動的に行う。これにより、心拍が安定した際に素早く計測モードへ移行することができる。更に、フリーズ操作及び計測モードへの移行操作が不要になるために、操作者の操作負担をより軽減することができる。
【００６７】
なお、本実施例における超音波診断装置１の構成はこれに限定されず、適宣構成や動作を省略し、あるいは追加しても構わない。例えば図２の画面表示例では標準ドプラ波形２００と拡大ドプラ波形２１０とを並べて表示するものとして述べたが、標準ドプラ波形２００の表示を省略しても構わない。あるいはドプラ波形処理部１７が心拍の安定を検出したときに標準ドプラ波形２００の表示を消去し、代わりに拡大ドプラ波形２１０を表示させる動作を行なっても構わない。
【００６８】
以上に本発明のいくつかの実施例を説明したが、これらの実施例は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施例は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施例やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
【符号の説明】
【００６９】
１      超音波診断装置
１０    システム制御部
１２    記憶部
１３    送受信部
１５    超音波画像生成部
１７    ドプラ波形処理部
１８    Ｂモード処理部
１９    ドプラ処理部
２１    入力部
２３    表示部
３０    超音波プローブ
３１    超音波振動子
４０    生体信号検出ユニット
【図１】
0005820128.tif 000002
【図２】
0005820128.tif 000003
【図３】
0005820128.tif 000004
【図４】
0005820128.tif 000005
【図５】
0005820128.tif 000006
【図６】
0005820128.tif 000007
【図７】
0005820128.tif 000008
【図８】
0005820128.tif 000009