(57)【要約】【課題】流体中の固体粒子と気泡とを互いに区別して識別できる超音波診断装置を提供すること。【解決手段】この発明の超音波診断装置（１）は、超音波を流体に対して照射する照射部（１０）と、超音波に由来した流体からのエコー信号を受信し、この受信した受信信号を周波数空間に変換して周波数スペクトルを取得する測定部（３６）と、周波数スペクトルに含まれた、粒子を表す特定のドプラ周波数を算出する第１算出部（４２）と、流体の平均流速を表すドプラ周波数と特定のドプラ周波数との間の周波数差を算出する第２算出部（４３）と、周波数差に基づいて、粒子が固体粒子であるか気泡であるかを判別する判定部（４３）とを備える。【選択図】図１
2020185281.tif 000002
【特許請求の範囲】
【請求項１】
流体中の粒子を識別する超音波診断装置であって、
超音波を上記流体に対して照射する照射部と、
上記超音波に由来した上記流体からのエコー信号を受信し、この受信した受信信号を周波数空間に変換して周波数スペクトルを取得する測定部と、
上記周波数スペクトルに含まれた、粒子を表す特定のドプラ周波数を算出する第１算出部と、
上記流体の平均流速を表すドプラ周波数と上記特定のドプラ周波数との間の周波数差を算出する第２算出部と、
上記周波数差に基づいて、上記粒子が固体粒子であるか気泡であるかを判別する判定部とを備えることを特徴とする超音波診断装置。
【請求項２】
請求項１の超音波診断装置において、
上記判定部は、上記周波数差に加えて、上記超音波の送信電力で正規化した上記受信信号の強度に基づいて、上記粒子が固体粒子であるか気泡であるかを判別することを特徴とする超音波診断装置。
【請求項３】
請求項２の超音波診断装置において、
上記判定部は、上記受信信号の強度として、上記周波数スペクトルの遅れ時間に伴う自己相関関数のピーク値を用いることを特徴とする超音波診断装置。
【請求項４】
請求項３の超音波診断装置において、
上記判定部は、上記自己相関関数の上記ピーク値として、上記遅れ時間が長くなるのに伴って現れる第３のピーク値を用いることを特徴とする超音波診断装置。
【請求項５】
流体中の粒子を識別する超音波診断方法であって、
超音波を上記流体に対して照射し、
上記超音波に由来した上記流体からのエコー信号を受信し、この受信した受信信号を周波数空間に変換して周波数スペクトルを取得し、
上記周波数スペクトルに含まれた、粒子を表す特定のドプラ周波数を算出し、
上記流体の平均流速を表すドプラ周波数と上記特定のドプラ周波数との間の周波数差を算出し、
上記周波数差に基づいて、上記粒子が固体粒子であるか気泡であるかを判別することを特徴とする超音波診断方法。
【請求項６】
請求項５に記載の超音波診断方法をコンピュータに実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、超音波診断装置及び超音波診断方法に関し、より詳しくは、流体中に存在する固体粒子などを検出する超音波診断装置及び超音波診断方法に関する。また、この発明は、そのような超音波診断方法をコンピュータに実行させるためのプログラムに関する。
【背景技術】
【０００２】
従来、この種の超音波診断装置としては、例えば特許文献１（特開２００１−１３７２４２号公報）に開示されているように、血流に対して超音波を照射し、血流からのエコー信号（受信信号）に基づいて血流の速度に比例して変化する超音波の周波数（ドプラ効果による周波数変化、ドプラ信号とも呼ばれる）を抽出し、そのドプラ周波数成分のパワーを演算することによって、パワースペクトラムを求めるものが知られている。パワースペクトラムには、時折、栓子（固体粒子）を表す一過性高輝度信号（High Intensity Transient Signals: HITS）が含まれる。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００１−１３７２４２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ところで、例えば、臓器移植手術を待つ患者が左心補助人工心臓（Left Ventricular Assist Device: LVASD）を装着した状態で長期間経過すると、ポンプがＬＶＡＳＤの拍動に伴って駆動することによって血流中に気泡が発生することがある。その他、カテーテルや手術器具の操作でも血流中に気泡が混入する。ここで、超音波診断装置で診断すると、これらの気泡も、固体粒子の場合と同様に一過性高輝度信号を発生する。このため、流体（血流）中に存在する固体粒子と気泡とを互いに区別して識別することが求められる。
【０００５】
そこで、この発明の課題は、流体中の固体粒子と気泡とを互いに区別して識別できる超音波診断装置及び超音波診断方法を提供することにある。また、この発明の課題は、そのような超音波診断方法をコンピュータに実行させるためのプログラムを提供することにある。
【課題を解決するための手段】
【０００６】
上記課題を解決するため、この開示の超音波診断装置は、
流体中の粒子を識別する超音波診断装置であって、
超音波を上記流体に対して照射する照射部と、
上記超音波に由来した上記流体からのエコー信号を受信し、この受信した受信信号を周波数空間に変換して周波数スペクトルを取得する測定部と、
上記周波数スペクトルに含まれた、粒子を表す特定のドプラ周波数を算出する第１算出部と、
上記流体の平均流速を表すドプラ周波数と上記特定のドプラ周波数との間の周波数差を算出する第２算出部と、
上記周波数差に基づいて、上記粒子が固体粒子であるか気泡であるかを判別する判定部とを備えることを特徴とする。
【０００７】
本明細書で、「流体」とは、主に空間を移動し得る液体を意味し、粒子を含んでいても良い。「粒子」とは、例えば、固体粒子、気泡（気体粒子）を総称する。
【０００８】
また、流体の「平均流速」とは、上記流体について、或る時点で、流速を空間にわたって平均的に見た速度を意味する。
【０００９】
また、「特定のドプラ周波数」とは、例えば、固体または気体からなる粒子の移動速度を表すドプラ周波数を指す。
【００１０】
また、流体の平均流速を表すドプラ周波数と特定のドプラ周波数との間の「周波数差」とは、例えば、超音波照射によって粒子に音響放射力が作用して、その結果粒子が加速した後の流速変化分に相当するドプラ周波数を指す。これは、音響放射力が媒質と粒子の音響インピーダンスの差に依存して強く作用するため、気泡と血液の音響インピーダンスの差が固体粒子と血液の音響インピーダンスの差に対して極めて大きいことから、気泡には音響放射力が強く作用するため、当該周波数差が大きくなる、という事実に基づいている。
【００１１】
この開示の超音波診断装置では、照射部は超音波を流体中の被検体に対して照射する。測定部は上記超音波に由来した上記流体からのエコー信号を受信し、この受信した受信信号を周波数空間に変換して周波数スペクトルを取得する。第１算出部は上記周波数スペクトルに含まれた、粒子を表す特定のドプラ周波数を算出する。第２算出部は上記流体の平均流速を表すドプラ周波数と上記特定のドプラ周波数との間の周波数差を算出する。判定部は上記周波数差に基づいて、上記粒子が固体粒子であるか気泡であるかを判別する。上述のように、上記周波数差は気泡のように音響インピーダンスの差が大きいと、上記周波数差が大きくなる。したがって、上記周波数差に基づいて、流体中の固体粒子と気泡とを互いに区別して識別することが可能となる。
【００１２】
一実施形態の超音波診断装置では、上記判定部は、上記周波数差に加えて、上記超音波の送信電力で正規化した上記受信信号の強度に基づいて、上記粒子が固体粒子であるか気泡であるかを判別することを特徴とする。
【００１３】
この一実施形態の超音波診断装置では、上記判定部は、上記周波数差に加えて、上記超音波の送信電力で正規化した上記受信信号の強度に基づいて、上記粒子が固体粒子であるか気泡であるかを判別する。その結果、流体中の固体粒子と気泡とを高精度で判別することが可能となる。
【００１４】
一実施形態の超音波診断装置では、上記判定部は、上記受信信号の強度として、上記周波数スペクトルの遅れ時間に伴う自己相関関数のピーク値を用いることを特徴とする。
【００１５】
この一実施形態の超音波診断装置では、上記判定部は、上記受信信号の強度として、上記周波数スペクトルの遅れ時間に伴う自己相関関数のピーク値を用いる。その結果、流体中の固体粒子と気泡とをより高精度で判別することが可能となる。
【００１６】
一実施形態の超音波診断装置では、上記判定部は、上記自己相関関数の上記ピーク値として、上記遅れ時間が長くなるのに伴って現れる第３のピーク値を用いることを特徴とする。
【００１７】
この一実施形態の超音波診断装置では、上記周波数差（Δｆ）の演算結果を流用でき、判別に用いる計算を簡易にすることが可能となる。
【００１８】
別の局面では、この開示の超音波診断方法は、
流体中の粒子を識別する超音波診断方法であって、
超音波を上記流体に対して照射し、
上記超音波に由来した上記流体からのエコー信号を受信し、この受信した受信信号を周波数空間に変換して周波数スペクトルを取得し、
上記周波数スペクトルに含まれた、粒子を表す特定のドプラ周波数を算出し、
上記流体の平均流速を表すドプラ周波数と上記特定のドプラ周波数との間の周波数差を算出し、
上記周波数差に基づいて、上記粒子が固体粒子であるか気泡であるかを判別することを特徴とする。
【００１９】
この一実施形態の超音波診断方法では、上記周波数差に基づいて、流体中の固体と気体とを判別することが可能となる。
【００２０】
さらに別の局面では、この開示のプログラムは、上記超音波診断方法をコンピュータに実行させるためのプログラムである。
【００２１】
この開示のプログラムをコンピュータに実行させることによって、上記超音波診断方法を実施することができる。
【発明の効果】
【００２２】
以上より明らかなように、この開示の超音波診断装置及び超音波診断方法は、流体中の固体粒子と気泡とを互いに区別して識別することが可能となる。また、この開示のプログラムをコンピュータに実行させることによって、上記超音波診断方法を実施することができる。
【図面の簡単な説明】
【００２３】
【図１】この発明の一実施形態の超音波診断装置のブロック構成を模式的に示す図である。
【図２】上記超音波診断装置の超音波プローブでの超音波照射による粒子の移動状態を模式的に示す図である。
【図３】エコー信号からドプラ周波数を推定するための直交検波部のブロック構成を模式的に示す図である。
【図４】上記超音波診断処理（超音波診断方法）のメインフローを示す図である。
【図５】自己相関のピーク算出部が入力する直交検波出力信号を示す図である。
【図６】ドプラ信号の周波数スペクトルの、遅れ時間に伴う自己相関関数を示す図である。
【図７】時間経過に伴う心臓の拍動によって脈動する血流の周波数スペトログラムを示す図である。
【図８】ｆＢａｃｋ，ｆＭＥＳ算出部が入力する直交検波出力信号を示す図である。
【図９】血流のドプラ周波数に対するパワースペクトルの強度を示す図である。
【図１０】粒子についてドプラ周波数差に対する自己相関関数の第３ピーク値を表す図である。
【発明を実施するための形態】
【００２４】
以下、この発明の実施の形態を、図面を参照しながら詳細に説明する。
【００２５】
（超音波診断装置の構成）
図１は、この発明の一実施形態の超音波診断装置１のブロック構成を模式的に示している。この超音波診断装置１は、照射部である超音波プローブ１０と超音波診断装置本体２０を備え、超音波プローブ１０を用いて血流中に存在する栓子と気泡とを互いに区別して識別する。
【００２６】
超音波診断装置本体２０は、制御部３０と、操作部３１と、記憶部３２と、送受信部３３と、直交検波部３４と、短時間フーリエ変換部３５と、測定部であるドプラ周波数スペクトル取得部３６と、表示部３７と、演算部４０とを備えている。演算部４０は、自己相関のピーク値算出部４１と、第１算出部であるｆＢａｃｋ，ｆＭＥＳ算出部４２と、第２算出部である判定部４３とを備えている。
【００２７】
超音波プローブ１０は、測定対象である生体組織内の血管を流れる血流へ超音波を照射し、照射した超音波が血流に含まれる赤血球、粒子（栓子、気泡）などから反射することにより得られるエコー信号を受信する。この例では、超音波の送信周波数＝２ＭＨｚであり、送信波のパルス幅＝１μ秒であった。
【００２８】
図２に示すように、生体組織表面２００に設置された超音波プローブ１０から、生体組織内の血管１００に超音波が照射され、血流からのエコー信号が、超音波プローブ１０で受信される。これら照射される超音波と受信されるエコー信号は、超音波ビームを形成している。
【００２９】
図２によって分かるように、この例では、流体としての血流Ｂは、血液の成分としての赤血球Ｒに加えて、粒子（栓子、気泡）Ｐを含んでいるものとする。赤血球Ｒの円盤状の直径は７〜８μｍであり、栓子はこれよりも大きいものとする。血流に含まれる赤血球Ｒ、粒子（栓子、気泡）Ｐは、その速度は一定ではない。さらに、すべての赤血球Ｒが同じ方向に移動するわけではなく、それぞれが様々な方向へ移動している。これを総合的に見ると、血流Ｂは、或る時点では、血管１００内を或る平均流速Ｖによって流れていると考えられる。この例では、粒子Ｐが超音波ビーム内に１つ含まれ、超音波の照射によって速度Ｖｐで移動する様子を模式的に示している。
【００３０】
上述のように、気泡と血液の音響インピーダンスの差が栓子と血液の音響インピーダンスの差に対して極めて大きいことから、粒子Ｐに超音波が照射されると、気泡は、超音波が照射されている時間で加速して移動速度が大きくなる。これは音響放射力が気泡に対して強く作用し、栓子に対して音響放射力が弱く作用するためである。この超音波診断装置１は、この事実に基づいて、血流中の粒子が栓子であるか気泡であるかを判別する（詳しくは後述する。）。
【００３１】
図３に示すように、直交検波部３４は、エコー信号ｃｏｓ２π（ｆ０＋ｆｄ）ｔからドプラ周波数ｆｄを推定するために直交検波を実施し、Ｉ（ｔ）＝ｃｏｓ２πｆｄｔと、Ｑ（ｔ）＝−ｓｉｎ２πｆｄｔを得る。ここで、ドプラ周波数ｆｄは、次の式Ｅｑ１によって表される。
【００３２】
【数１】
2020185281.tif 000003
・・・（Ｅｑ１）
【００３３】
図１に示すように、短時間フーリエ変換部３５は、Ｉ信号及びＱ信号から窓関数によって切り出し、その短時間フーリエ関数を演算する。
【００３４】
ドプラ周波数スペクトル取得部３６は、短時間フーリエ変換部３５によって計算されたスペクトルを垂直な周波数軸に１次元分布として、時間経過とともに取得する。
【００３５】
演算部４０は、ソフトウェア（コンピュータプログラム）によって動作するＣＰＵ（中央演算処理ユニット）を含み、後述の超音波診断方法に従う処理や、その他各種処理を実行する。
【００３６】
表示部３７は、この例では、ＬＣＤ（液晶ディスプレイ）からなり、制御部３０からの制御信号によって、ドプラ周波数スペクトルを時間経過とともに表示する。さらに、演算部４０からの制御信号によって、判定部４３で判別された粒子Ｐについて、ディスプレイ上に栓子と気泡とを互いに区別して表示する。この例では、血流の周波数スペクトログラム画面上に、「栓子」であるか「気泡」であるかを重畳させて表示させることも好適である。
【００３７】
操作部３１は、この例では、キー入力スイッチ及び設定スイッチからなり、ユーザ（操作者）からの指示及びデータなどを入力または設定するために用いられる。
【００３８】
制御部３０は、操作部３１と、送受信部３３と、直交検波部３４と、短時間フーリエ変換部３５と、ドプラ周波数スペクトル取得部３６と、表示部３７と、演算部４０を制御し、この制御情報は、記憶部３２に記憶される。
【００３９】
図１中に示すように、この例では、プログラムされた演算部４０によって、自己相関のピーク値算出部４１と、ｆＢａｃｋ，ｆＭＥＳ算出部４２と、判定部４３とが構成されている。これらの動作については、次に述べる図３の超音波診断処理のフローの中で説明する。
【００４０】
（超音波診断方法）
図４は、超音波診断装置１の制御部３０及び演算部４０が実行する超音波診断処理のフローを示している。ユーザ（操作者）による、図１中に示した操作部３１を介した開始指示を受けると制御部３０及び演算部４０は、この超音波診断処理のフローを開始する。
【００４１】
まず、図４のステップＳ１に示すように、制御部３０は、送受信部３３を制御して、超音波プローブ１０から超音波を生体組織内の血管を流れる血流へ照射する。
【００４２】
次に、図４のステップＳ２に示すように、制御部３０は、送受信部３３を制御して、超音波プローブ１０から超音波に由来した血流からのエコー信号を受信する。
【００４３】
次に、図４のステップＳ３に示すように、制御部３０は、ドプラ周波数スぺクトル取得部３６を制御して周波数スペクトルを取得する。より詳しくは、制御部３０は、直交検波部３４を制御して、直交検波部３４からＩ信号とＱ信号を検波する。次いで、制御部３０は、短時間フーリエ変換部３５へＩ信号とＱ信号を入力し、出力をドプラ周波数スペクトル取得部３６に出力する。これによって、制御部３０は、周波数スペクトルを取得することが可能となる。
【００４４】
次に、図４のステップＳ４に示すように、演算部４０は自己相関のピーク値算出部４１として働いて、周波数スペクトルから、自己相関関数のピーク値を算出する。より詳しくは、制御部３０は、直交検波出力信号を短時間フーリエ変換して周波数スペクトルを取得し、次の式Ｅｑ２によって、自己相関関数Ａ（ｔ）を算出する。ここで、直交検波出力信号Ｉ（ｔ）とＱ（ｔ）をそれぞれ実数部と虚数部とした複素信号をｄ（τ）とし、遅れ時間をｔとする。この例では、図５に、自己相関のピーク算出部が入力する直交検波出力信号の実数部Ｉ（ｔ）を示す。
【００４５】
【数２】
2020185281.tif 000004
・・・（Ｅｑ２）
【００４６】
図６は、ドプラ信号の周波数スペクトルの、遅れ時間に対する自己相関関数を示す。図６によって分かるように、遅れ時間０秒に対する自己相関関数のピーク値が第１ピークＰ１として表される。次いで、自己相関関数は、遅れ時間に伴って順に、第２ピーク値Ｐ２、第３ピーク値Ｐ３が表される。これは、図３に示すように、直交検波出力信号が、上述の周期的な２つの信号ｃｏｓ２πｆｄｔとｓｉｎ２πｆｄｔという信号で表されるため、自己相関関数も周期的にピーク値を示すためである。
【００４７】
自己相関関数のピーク値は、第３ピーク値Ｐ３を用いると、ピークの遅れ時間からＨＩＴＳのドプラ周波数の近傍を推定できるため、後述する周波数差（Δｆ）の演算に流用でき、判別に用いる計算を簡易にすることが可能となる。ただし、自己相関関数のピーク値は、第３ピーク値Ｐ３に、限られるものではない。第１ピーク値Ｐ１、第２ピーク値Ｐ２を用いてもよい。
【００４８】
なお、自己相関関数のピーク値に代えて、超音波の送信電力で正規化した受信信号強度を用いてもよい。これによれば、血流中の栓子と気泡とを高精度で判別することが可能となる。
【００４９】
上述のように、気泡と血液の音響インピーダンスの差が栓子と血液の音響インピーダンスの差に対して極めて大きいことから、粒子に超音波が照射されると、気泡の方が、栓子よりも超音波の反射強度が大きくなる。したがって、粒子に超音波が照射された反射強度、すなわち、超音波の送信電力で正規化した受信信号強度（または、自己相関関数のピーク値）は、栓子に比して、気泡の方が大きい特性を有する。この結果、流体中の栓子と気泡とを互いに区別して識別することが可能となる。このように、超音波の送信電力で正規化した受信信号強度（または、自己相関関数のピーク値）に基づいて、流体中の栓子と気泡とを互いに区別して識別することが可能である。しかしながら、後述する周波数差（Δｆ）を考慮するとより精度よく識別することが可能となる。
【００５０】
次に、図４のステップＳ５に示すように、演算部４０はｆＢａｃｋ，ｆＭＥＳ算出部４２として働いて、ｆＢａｃｋ，ｆＭＥＳを算出する。より詳しくは、図９に、血流のドプラ周波数に対するパワースペクトルの強度を示す。図９によって分かるように、ｆＢａｃｋ，ｆＭＥＳ算出部４２は、ｆＢａｃｋを、血流の平均流速を表すドプラ周波数として、ドプラスペクトルの強度が４０ｄＢ以上の最大周波数の値を算出する。ｆＭＥＳを、粒子のドプラ周波数として、ドプラスペクトルの強度が最大値を示す周波数の値を算出する。この例では、制御部３０は、図８に示すような、直交検波出力信号を短時間フーリエ変換して周波数スペクトルを取得する。
【００５１】
ここで、演算部４０は判定部４３として働いて、ドプラ周波数差Δｆを、次の式Ｅｑ３によって算出する。
【００５２】
【数３】
2020185281.tif 000005
・・・（Ｅｑ３）
【００５３】
図７は、時間経過に伴う心臓の拍動によって脈動する血流の周波数スペクトログラムを示す。図７中（両矢印）によって分かるように、ドプラ周波数差Δｆは、周波数軸上に、粒子（栓子，気泡）に超音波が照射され加速された移動速度変化として表される。
【００５４】
次に、図４のステップＳ６に示すように、演算部４０は判定部４３として働いて、自己相関関数の第３ピーク値Ｐ３が閾値Ｔｈｐより小さいか否か、かつ、ｆＢａｃｋ−ｆＭＥＳ、すなわち、Δｆが閾値Ｔｈｇより小さいか否か判断する（ステップＳ６）、自己相関関数の第３ピーク値Ｐ３が閾値Ｔｈｐより小さく、かつ、Δｆが閾値Ｔｈｇより小さい場合（ステップＳ６でＹＥＳ）、ステップＳ７に進み、粒子が「栓子」であると判別して、表示部３７において「栓子」と表示し処理を終了する。
【００５５】
一方、自己相関関数の第３ピーク値Ｐ３が閾値Ｔｈｐ以上であるか、または、Δｆが閾値Ｔｈｇ以上である場合（ステップＳ６でＮＯ）、ステップＳ７に進み、粒子が「気泡」であると判別して、表示部３７において「気泡」と表示し処理を終了する。
【００５６】
この例では、自己相関関数の第３ピーク値Ｐ３に対する閾値Ｔｈｐ＝０．０５５であり、ドプラ周波数差Δｆに対する閾値Ｔｈｇ＝５５０（Ｈｚ）であるものとする。ただし、これに限られるものではない。ユーザは、適切な閾値Ｔｈｐ，Ｔｈｇを任意に設定することが可能である。
【００５７】
（検証）
図１０は、粒子について周波数差Δｆに対する自己相関関数の第３ピーク値Ｐ３を表している。図１０中、黒丸は、栓子を、白丸は気泡を表している。斜線の丸（LVAS）は、補助人工心臓を装着した際の栓子及び気泡の分布を表している。
【００５８】
図１０によって分かるように、栓子は、自己相関関数の第３ピーク値が閾値Ｔｈｐ＝０．０５５より小さく、かつ、Δｆが閾値Ｔｈｇ＝５５０（Ｈｚ）より小さい領域に、分布している。一方、気泡は、自己相関関数の第３ピーク値が閾値Ｔｈｐ以上であるか、または、Δｆが閾値Ｔｈｇ以上である領域に分布している。したがって、自己相関関数の第３ピーク値Ｐ３及びドプラ周波数差Δｆという２つのパラメータによって、流体中の栓子と気泡とを互いに区別して精度よく識別することが可能となる。
【００５９】
以上のように、気泡と血液の音響インピーダンスの差が栓子と血液の音響インピーダンスの差に対して極めて大きいことから、粒子に超音波が照射されると、気泡の方が、栓子よりも超音波の反射強度、及び音響放射力による移動速度変化が大きくなる。よって、粒子に超音波が照射された反射強度、すなわち、自己相関関数のピーク値、及び粒子に超音波が照射され加速された移動速度変化、すなわち、ドプラ周波数差Δｆという２つのパラメータを用いると、流体中の栓子と気泡とを互いに区別してより精度よく識別することが可能となる。
【００６０】
なお、上述の超音波診断方法を、ソフトウェア（コンピュータプログラム）として、ＣＤ（コンパクトディスク）、ＤＶＤ（デジタル多用途ディスク）、フラッシュメモリなどの非一時的（non-transitory）にデータを記憶可能な記録媒体に記録してもよい。このような記録媒体に記録されたソフトウェアを、パーソナルコンピュータ、ＰＤＡ（パーソナル・デジタル・アシスタンツ）、スマートフォン、ＰＬＣ（プログラマブルロジックコントローラ）などの実質的なコンピュータ装置にインストールすることによって、それらのコンピュータ装置に、上述の超音波診断方法を実行させることができる。
【００６１】
以上の実施形態は例示であり、この発明の範囲から離れることなく様々な変形が可能である。上述した複数の実施の形態は、それぞれ単独で成立し得るものであるが、実施の形態同士の組みあわせも可能である。また、異なる実施の形態の中の種々の特徴も、それぞれ単独で成立し得るものであるが、異なる実施の形態の中の特徴同士の組みあわせも可能である。
【符号の説明】
【００６２】
１    超音波診断装置
１０  超音波プローブ
２０  超音波診断装置本体
３５  短時間フーリエ変換部
３６  ドプラ周波数スペクトル取得部
４０  演算部
４１  自己相関のピーク値算出部
４２  ｆＢａｃｋ，ｆＭＥＳ算出部
４３  判定部
【図１】
2020185281.tif 000006
【図２】
2020185281.tif 000007
【図３】
2020185281.tif 000008
【図４】
2020185281.tif 000009
【図５】
2020185281.tif 000010
【図６】
2020185281.tif 000011
【図７】
2020185281.tif 000012
【図８】
2020185281.tif 000013
【図９】
2020185281.tif 000014
【図１０】
2020185281.tif 000015