(57)【要約】  本発明の超音波診断装置は、探触子を繰り返し駆動し、前記探触子を駆動することにより送信された超音波が被検体において反射することにより得られた反射エコーを前記探触子によって受信し、複数の受信信号を逐次生成する送受信部と、前記複数の受信信号に基づき、各フレームにおける前記被検体中の血流部分の血流速度データを逐次生成するカラーフローマッピング信号処理部と、前記各フレームにおける血流速度データに対して残像処理を行う残像処理部と、前記受信信号に基づき、Ｂモード断層画像フレームデータを生成する断層画像信号処理部と、前記残像処理された血流速度データおよび前記Ｂモード断層画像フレームデータを合成する画像合成部とを備え、前記残像処理部は、最新のフレームの血流速度データと最新より以前のフレームの残像処理された血流速度データに基づいて折り返し判定を行い、前記折り返し判定の結果と最新および最新より以前のフレームの血流速度データに基づいて動的にパーシスタンス係数を変更する。
2011036891.tif 000002
【特許請求の範囲】
【請求項１】
探触子を繰り返し駆動し、前記探触子を駆動することにより送信された超音波が被検体において反射することにより得られた反射エコーを前記探触子によって受信し、複数の受信信号を逐次生成する送受信部と、
前記複数の受信信号に基づき、各フレームにおける前記被検体中の血流部分の血流速度データを逐次生成するカラーフローマッピング信号処理部と、
前記各フレームにおける血流速度データに対して残像処理を行う残像処理部と、
前記受信信号に基づき、Ｂモード断層画像フレームデータを生成する断層画像信号処理部と、
前記残像処理された血流速度データおよび前記Ｂモード断層画像フレームデータを合成する画像合成部と
を備え、
前記残像処理部は、
最新のフレームの血流速度データと最新より以前のフレームの残像処理された血流速度データに基づいて折り返し判定を行い、前記折り返し判定の結果と最新および最新より以前のフレームの血流速度データに基づいて動的にパーシスタンス係数を変更する超音波診断装置。
【請求項２】
前記残像処理部は、
前記最新のフレームの血流速度データを記憶する第１のメモリ部と、
前記最新より以前のフレームの残像処理された血流速度データを記憶する第２のメモリ部と、
前記第１のメモリ部および前記第２のメモリ部から血流速度データをそれぞれ読み出し、前記折り返し判定を行う折り返し判定部と、
前記折り返し判定の結果および前記第１のメモリ部に記憶された血流速度データに応じてパーシスタンス係数を決定するパーシスタンス係数決定部と、
前記パーシスタンス係数および前記折り返し判定の結果に基づいて、前記第１のメモリ部に記憶された血流速度データに対しパーシスタンス演算を行い、演算結果を前記残像処理された血流速度データとして出力するパーシスタンス演算部とを含む請求項１に記載の超音波診断装置。
【請求項３】
前記折り返し判定部は、前記第１のメモリ部に記憶された血流速度データおよび前記第２のメモリ部に記憶された血流速度データと複数の閾値とを比較することによって、折り返しが発生しているかどうか、および、前記最新のフレームの血流速度データが折り返し領域にあるかどうかを判定する請求項２に記載の超音波診断装置。
【請求項４】
前記残像処理部は、前記血流速度の値に応じて対応付けられた異なる２つ以上の値のパーシスタンス係数を含む参照テーブルを記憶している第３のメモリ部をさらに含む、請求項２または３に記載の超音波診断装置。
【請求項５】
前記参照テーブルは、所定の値以上の血流速度に対して一定の値のパーシスタンス係数が対応付けられている請求項４に記載の超音波診断装置。
【請求項６】
前記残像処理部は、
前記最新のフレームの血流速度データを記憶する第１のメモリ部と、
前記最新より以前のフレームの残像処理された血流速度データを記憶する第２のメモリ部と、
前記第１のメモリ部および前記第２のメモリ部から血流速度データをそれぞれ読み出し、前記折り返し判定を行う折り返し判定部と、
前記折り返し判定の結果および前記第１のメモリ部に記憶された血流速度データに応じて第１のパーシスタンス係数を決定する第１のパーシスタンス係数決定部と、
前記第１のパーシスタンス係数および前記折り返し判定の結果に基づいて前記第１のメモリ部に記憶された血流速度データに対しパーシスタンス演算を行う第１のパーシスタンス演算部と、
前記折り返し判定の結果および前記第２のメモリ部に記憶された血流速度データに応じて第２のパーシスタンス係数を決定する第２のパーシスタンス係数決定部と、
前記第２のパーシスタンス係数および前記折り返し判定の結果に基づいて前記第１のメモリ部に記憶された血流速度データに対しパーシスタンス演算を行う第２のパーシスタンス演算部と、
前記第１のパーシスタンス演算部から出力される演算結果の絶対値および前記第２パーシスタンス演算部から出力される演算結果の絶対値を比較し、大きい方の演算結果を前記残像処理された血流速度データとして出力する最大値選択部と、
を含む請求項１に記載の超音波診断装置。
【請求項７】
前記折り返し判定部は、前記第１のメモリ部に記憶された血流速度データおよび前記第２のメモリ部に記憶された血流速度データと複数の閾値とを比較することによって、折り返しが発生しているかどうか、および、前記最新のフレームの血流速度データが折り返し領域にあるかどうかを判定する請求項６に記載の超音波診断装置。
【請求項８】
前記残像処理部は、
前記血流速度の値に応じて対応付けられた異なる２つ以上の値の第１のパーシスタンス係数を含む第１の参照テーブルを記憶している第３のメモリ部と、
前記血流速度の値に応じて対応付けられた異なる２つ以上の値の第２のパーシスタンス係数を含む第２の参照テーブルを記憶している第４のメモリ部と、
をさらに含む、請求項６または７に記載の超音波診断装置。
【請求項９】
前記第１の参照テーブルと前記第２の参照テーブルにおいて、同じ血流速度の値に応じて対応付けられた第１のパーシスタンス係数と第２のパーシスタンス係数とは互いに異なる値である請求項８に記載の超音波診断装置。
【請求項１０】
前記第１の参照テーブルは、所定の値以上の血流速度に対して一定の値のパーシスタンス係数が対応付けられている請求項８に記載の超音波診断装置。
【請求項１１】
前記最新より以前のフレームの残像処理された血流速度データは、前記最新より１つ前のフレームの残像処理された血流速度データである請求項１から１０のいずれかに記載の超音波診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波診断装置に関し、特にカラーフローマッピング時の残像処理方法に関する。
【背景技術】
【０００２】
超音波診断装置は、超音波を被検体に送信し、その反射エコーに含まれる情報を解析することにより、被検体内の画像を作成する。カラーフローマッピング（Color Flow Mapping、以下ＣＦＭと略す場合がある）と呼ばれる手法によって、被検体内の血流を画像化することも可能であり、医療分野全般において、血流状態を表示することのできる超音波診断装置が広く利用されている。【０００３】
カラーフローマッピングはカラードップライメージング（Color Doppler Imaging、ＣＤＩ）とも呼ばれ、ドップラ効果を利用する。超音波が血流を照射すると、ドプラ効果により反射エコーに血流速度に応じたドップラ偏移が生じる。このドップラ偏移の情報を直交検波によって検出し、ＭＴＩ（Moving Target Indicator）フィルタと呼ばれるハイパスフィルタ処理、自己相関処理およびノイズカット処理を施すことによって、血流速度に関する情報が得られる。得られた血流速度の情報を色情報に変換し、二次元的にＢモード断層画像に重畳して表示することによって、被検体内の血流状態をユーザに認識させることができる。【０００４】
血流から得られる反射エコーによる受信信号強度は、Ｂモード断層画像生成に用いる組織散乱体および組織境界から得られる反射エコーによる受信信号強度に比べ、かなり小さい。このため、カラーフローマッピングにおける信号処理によって得られる血流速度および血流パワー（移動する血流量）は不安定になりやすい。【０００５】
特に、観察したい部分の血流速度が遅い場合、あるいは、観察したい部分が末梢血管である場合、血流パワーが小さくなるため、本来、システムノイズおよび音響ノイズのみをカットするべきノイズカット処理において、血流速度あるいは血流パワーに関する情報が除去されてしまい易い。その結果、血流画像において本来血流として表示される部分が黒く抜けてしまう現象が発生する。例えば、１秒間に数フレームから数十フレームの割合で、被検体内の血流を画像化する場合、そのうちのいくつかのフレームにおいて、血流部分が黒く示されてしまう。このため、断層画像中の血流部分が突然消滅し、画像が滑らかでなくなったり、違和感のあるものとなったりする。【０００６】
この問題を解決するために、従来のカラーフローマッピングを行う超音波診断装置では、信号処理後段においてパーシスタンス処理（残像処理）と呼ばれる時間方向補間を行うのが通例である。以下、特許文献１に示される従来のカラーフローマッピングにおけるパーシスタンス処理を説明する。【０００７】
図７に示す従来の超音波診断装置において、超音波送受信部４０２は、探触子４０１を駆動し、超音波を被検体へ送信する。また、探触子４０１によって被検体で生じる反射エコーを受信し、受信信号を生成する。Ｂモード断層画像を生成する場合、超音波送受信部４０２は、Ｂモード断層画像の生成に適した送受信を行い、得られた受信信号を断層画像信号処理部４０９へ出力する。カラーフローマッピング断層画像を生成する場合は、カラーフローマッピング断層画像の生成に適した送受信を行い、得られた受信信号をカラーフローマッピング処理部４０３（以下、ＣＦＭ信号処理部と略す）へ出力する。一般的に、カラーフローマッピング断層画像を生成する場合、安定したカラーフローマッピング断層画像を得るために、超音波送受信部４０２は、同じ音響線上において超音波の送受信を複数回行う。【０００８】
ＣＦＭ信号処理部４０３は、受信信号に対して直交検波処理、ＭＴＩフィルタ処理および自己相関処理を行い、血流速度および血流パワーを算出したのち、システムあるいは音響ノイズを排除するノイズカット処理を行い、フレームメモリ部４０４に血流速度および血流パワーを出力する。【０００９】
フレームメモリ部４０４は、リングバッファで構成されており、現在の走査フレームからＮフレーム前（Ｎは、１以上の整数）までの血流速度および血流パワーをフレーム単位で保存している。ここで、フレームとは、１画面のＣＦＭ断層画像を構成する血流速度データおよび血流パワーデータ群を示す。【００１０】
フレームメモリ選択部４０５は、フレームメモリ部４０４から、あらかじめ設定された複数のＣＦＭフレームデータを選択し、パーシスタンス演算部４０７へ出力するようフレームメモリ部４０４に指令を出力する。パーシスタンス演算部４０７は、フレームメモリ部４０４から読み出されたＣＦＭフレームデータとパーシスタンス係数設定部４０６から出力されるパーシスタンス係数に基づいてパーシスタンス演算を実施し、ＣＦＭ  ＤＳＣ（Digital Scan Converter）部４０８に出力する。パーシスタンス演算は、単純な重み付け演算であり、パーシスタンス係数設定部４０６から出力されるパーシスタンス係数は、あらかじめシステムから設定された固定係数である。【００１１】
ＣＦＭ  ＤＳＣ部４０８は、パーシスタンス演算部４０７から出力されるＣＦＭフレームデータの座標を変換し、画像合成部４１１へ出力する。【００１２】
断層画像信号処理部４０９は受信信号に対してダイナミックフィルタ処理を施すことにより、不要なノイズをカットしたのち、包絡線検波処理およびダイナミックレンジ圧縮処理を施して断層画像ＤＳＣ部４１０に断層画像フレームデータを出力する。断層画像ＤＳＣ部４１０は、断層画像信号処理部４０９からの断層画像フレームデータの座標を変換し、画像合成部４１１へ出力する。【００１３】
画像合成部４１１は、ＣＦＭ  ＤＳＣ部４１０および断層画像ＤＳＣ部４１０から出力される各フレームデータを画素ごとに合成し、合成画像フレームデータを生成する。具体的には、血流速度がゼロである場合は、断層画像フレームデータを表示し、そうでない場合は、ＣＦＭフレームデータを表示するように２つのデータを画素ごと、あるいは、対応する測定点のデータごとに合成する。また、血流速度や血流の方向に応じてデータを色情報に変換し、表示部４１２に出力する。表示部４１２は、画像合成部４１１から受け取ったデータを表示する。
【先行技術文献】
【特許文献】
【００１４】
【特許文献１】特開平２-２８６１４０号公報
【発明の概要】
【発明が解決しようとする課題】
【００１５】
従来の超音波診断装置におけるパーシスタンス処理は、血流速度が遅く、あるいは、血流パワーが小さく、これらが不安定な状態にあることに起因してＣＦＭ信号処理部４０３の出力結果が不安定になり、表示血流に黒抜けが発生することを防止する。具体的には、現在走査中のフレームデータよりも過去のフレームデータに重みを置いたパースタンス係数を用いることによって、残像効果を生じさせ、画像に黒抜けが発生するのを抑制する。【００１６】
しかし、このような特徴は、血流速度が激しく変化する動脈の診断には適していない場合がある。例えば、頚動脈では、心臓の拡張収縮に応じて血流が激しく変化し、収縮期には、心周期に対してごく短い時間において血流が速くなり、拡張期では、血流が遅くなる。また、血流速度の最大値と最小値との差も他の診断部位に比べて大きい。その結果、拡張期における頚動脈の血流速度は、心周期に対して比較的長い時間、小さな値で一定となり、ＣＦＭ信号処理部４０３の出力が安定しない状態に陥る。【００１７】
表示する血流の画像に黒抜けが生じるのを抑制するためには、パーシスタンス処理によって、残像効果が高くなるようにパーシスタンス係数を設定することが好ましい。これにより、血流速度が低速となる場合でも、黒抜けが発生しない滑らかな動画表示が可能となるからである。しかし、この場合、収縮期における高速な血流速度の表示ができなくなる。【００１８】
また、従来の超音波診断装置におけるパーシスタンス処理は、末梢血管の診断には適していない場合がある。例えば、甲状腺、肝臓、腎臓などの臓器には、主流血管から派生する末梢血管が存在する。これらの臓器の診断にあたり、末梢血管構造を把握することは非常に重要である。【００１９】
末梢血管における時間的な血流変化は、比較的安定しているのであるが、血管が物理的に細いために、血流パワーが頚動脈や心臓の場合と比べ、極端に小さくなる傾向にある。従って、血流パワーが小さいことにより、ドップラ偏移の検出が不安定となり、ＣＦＭ信号処理部４０３の出力が結果として安定しない状態に陥る。【００２０】
このため、パーシスタンス処理を行わない場合、断層画像中の末梢血管は、時間的に点滅して表示されることになり、動画として見づらくなる。逆にパーシスタンス処理を行う場合、時間的方向に対して断層画像中の末梢血管の血流が平滑化されるため、パーシスタンス処理により末梢血管が消滅する可能性がある。この場合、末梢血管の検出率が著しく低下することになる。【００２１】
本発明は、このような従来技術の課題を解決し、血流速度が激しく変化する頚動脈などのような診断部位において、血流変化が明確に認識でき、かつ、低速な血流速度でも黒抜けが発生しない滑らかな血流動画を表示することができる超音波診断装置を提供することを目的とする。また、末梢血管など血流パワーが小さい血管部分でも見やすい動画表示を行うことのできる超音波診断装置を提供することを目的とする。
【課題を解決するための手段】
【００２２】
本発明の超音波診断装置は、探触子を繰り返し駆動し、前記探触子を駆動することにより送信された超音波が被検体において反射することにより得られた反射エコーを前記探触子によって受信し、複数の受信信号を逐次生成する送受信部と、前記複数の受信信号に基づき、各フレームにおける前記被検体中の血流部分の血流速度データを逐次生成するカラーフローマッピング信号処理部と、前記各フレームにおける血流速度データに対して残像処理を行う残像処理部と、前記受信信号に基づき、Ｂモード断層画像フレームデータを生成する断層画像信号処理部と、前記残像処理された血流速度データおよび前記Ｂモード断層画像フレームデータを合成する画像合成部とを備え、前記残像処理部は、最新のフレームの血流速度データと最新より以前のフレームの残像処理された血流速度データに基づいて折り返し判定を行い、前記折り返し判定の結果と最新および最新より以前のフレームの血流速度データに基づいて動的にパーシスタンス係数を変更する。【００２３】
ある好ましい実施形態の超音波診断装置は、探触子を繰り返し駆動し、前記探触子を駆動することにより送信された超音波が被検体において反射することにより得られた反射エコーを前記探触子によって受信し、複数の受信信号を逐次生成する送受信部と、前記複数の受信信号に基づき、各フレームにおける前記被検体中の血流部分の血流速度データを逐次生成するカラーフローマッピング信号処理部と、前記各フレームにおける血流速度データに対して残像処理を行う残像処理部と、前記受信信号に基づき、Ｂモード断層画像フレームデータを生成する断層画像信号処理部と、前記残像処理された血流速度データおよび前記Ｂモード断層画像フレームデータを合成する画像合成部とを備え前記残像処理部は、最新のフレームの血流速度データを記憶する第１のメモリ部と、最新より１つ前のフレームの残像処理された血流速度データを記憶する第２のメモリ部と、前記第１のメモリ部および前記第２のメモリ部から血流速度データをそれぞれ読み出し、折り返し判定を行う折り返し判定部と、前記折り返し判定の結果および前記第１のメモリ部に記憶された血流速度データに応じてパーシスタンス係数を決定するパーシスタンス係数決定部と、前記パーシスタンス係数および前記折り返し判定の結果に基づいて、前記第１のメモリ部に記憶された血流速度データに対しパーシスタンス演算を行い、演算結果を前記残像処理された血流速度データとして出力するパーシスタンス演算部とを含む。【００２４】
ある好ましい実施形態において、前記折り返し判定部は、前記第１のメモリ部に記憶された血流速度データおよび前記第２のメモリ部に記憶された血流速度データと複数の閾値とを比較することによって、折り返しが発生しているかどうか、および、前記最新のフレームの血流速度データが折り返し領域にあるかどうかを判定する。【００２５】
ある好ましい実施形態において、前記残像処理部は、前記血流速度の値に応じて対応付けられた異なる２つ以上の値のパーシスタンス係数を含む参照テーブルを記憶している第３のメモリ部をさらに含む。【００２６】
ある好ましい実施形態において、前記参照テーブルは、所定の値以上の血流速度に対して一定の値のパーシスタンス係数が対応付けられている。【００２７】
ある好ましい実施形態の超音波診断装置は、探触子を繰り返し駆動し、前記探触子を駆動することにより送信された超音波が被検体において反射することにより得られた反射エコーを前記探触子によって受信し、複数の受信信号を逐次生成する送受信部と、前記複数の受信信号に基づき、各フレームにおける前記被検体中の血流部分の血流速度データを逐次生成するカラーフローマッピング信号処理部と、前記各フレームにおける血流速度データに対して残像処理を行う残像処理部と、前記受信信号に基づき、Ｂモード断層画像フレームデータを生成する断層画像信号処理部と、前記残像処理された血流速度データおよび前記Ｂモード断層画像フレームデータを合成する画像合成部とを備え、前記残像処理部は、最新のフレームの血流速度データを記憶する第１のメモリ部と、最新より１つ前のフレームの残像処理された血流速度データを記憶する第２のメモリ部と、前記第１のメモリ部および前記第２のメモリ部から血流速度データをそれぞれ読み出し、折り返し判定を行う折り返し判定部と、前記折り返し判定の結果および前記第１のメモリ部に記憶された血流速度データに応じて第１のパーシスタンス係数を決定する第１のパーシスタンス係数決定部と、前記第１のパーシスタンス係数および前記折り返し判定の結果に基づいて前記第１のメモリ部に記憶された血流速度データに対しパーシスタンス演算を行う第１のパーシスタンス演算部と、前記折り返し判定の結果および前記第２のメモリ部に記憶された血流速度データに応じて第２のパーシスタンス係数を決定する第２のパーシスタンス係数決定部と、前記第２のパーシスタンス係数および前記折り返し判定の結果に基づいて前記第１のメモリ部に記憶された血流速度データに対しパーシスタンス演算を行う第２のパーシスタンス演算部と、前記第１のパーシスタンス演算部から出力される演算結果の絶対値および前記第２パーシスタンス演算部から出力される演算結果の絶対値を比較し、大きい方の演算結果を前記残像処理された血流速度データとして出力する最大値選択部とを含む。【００２８】
ある好ましい実施形態において、前記折り返し判定部は、前記第１のメモリ部に記憶された血流速度データおよび前記第２のメモリ部に記憶された血流速度データと複数の閾値とを比較することによって、折り返しが発生しているかどうか、および、前記最新のフレームの血流速度データが折り返し領域にあるかどうかを判定する。【００２９】
ある好ましい実施形態において、前記残像処理部は、前記血流速度の値に応じて対応付けられた異なる２つ以上の値の第１のパーシスタンス係数を含む第１の参照テーブルを記憶している第３のメモリ部と、前記血流速度の値に応じて対応付けられた異なる２つ以上の値の第２のパーシスタンス係数を含む第２の参照テーブルを記憶している第４のメモリ部とをさらに含む。【００３０】
ある好ましい実施形態において、前記第１の参照テーブルと前記第２の参照テーブルにおいて、同じ血流速度の値に応じて対応付けられた第１のパーシスタンス係数と第２のパーシスタンス係数とは互いに異なる値である。【００３１】
ある好ましい実施形態において、前記第１の参照テーブルは、所定の値以上の血流速度に対して一定の値のパーシスタンス係数が対応付けられている。【００３２】
ある好ましい実施形態において、前記最新より以前のフレームの残像処理された血流速度データは、前記最新より１つ前のフレームの残像処理された血流速度データである。
【発明の効果】
【００３３】
本発明によれば、最新のフレームの血流速度データおよび最新より前のフレームの血流速度データに基づき、折り返しの判定を行い、折り返し判定の結果および最新のフレームの血流速度データに基づき、パーシスタンス係数を動的に変更する。このため、血流変化が明確に認識でき、かつ、低速な血流速度でも黒抜けが発生しない滑らかな血流動画を表示することが可能な超音波診断装置が実現する。【００３４】
また、本発明によれば、最新のフレームの血流速度データに基づき決定されたパーシスタンス係数を用いて残像処理された血流速度、および、最新より前のフレームの血流速度データに基づき決定されたパーシスタンス係数を用いて残像処理された血流速度を求め、絶対値の大きいほうを選択して血流画像表示に用いる。このため、血流パワーが不安定となる末梢血管の血流表示を点滅させずに、かつ、平滑化に起因する末梢血管の血流消滅を起こさない血流動画を表示することが可能となる。
【図面の簡単な説明】
【００３５】
【図１】本発明による超音波診断装置の第１の実施形態を示すブロック図である。
【図２】（ａ）および（ｂ）は、第１の実施形態において、パーシスタンス係数を用いて血流速度データのパーシスタンス演算を行う場合に考慮する折り返し演算を説明する模式図である。
【図３】（ａ）は、第１の実施形態における折り返し判定を説明するための模式図であり、（ｂ）は、参照テーブルのデータが満たす関係を示すグラフである。
【図４】本発明による超音波診断装置の第２の実施形態を示すブロック図である。
【図５】（ａ）および（ｂ）は、第２の実施形態において、パーシスタンス係数を用いて血流速度データのパーシスタンス演算を行う場合に考慮する折り返し演算を説明する模式図である。
【図６】（ａ）は、第２の実施形態における折り返し判定を説明するための模式図であり、（ｂ）および（ｃ）は、第１および第２の参照テーブルのデータが満たす関係を示すグラフである。
【図７】従来の超音波診断装置を示すブロック図である。
【発明を実施するための形態】
【００３６】
（第１の実施形態）
以下、図面を参照しながら本発明による超音波診断装置の第１の実施形態を説明する。図１は、本発明による超音波診断装置の第１の実施形態を示すブロック図である。図１に示す超音波診断装置１１は、探触子１０１と、超音波送受信部１０２と、ＣＦＭ信号処理部１０３と、残像処理部１１５と、断層画像信号処理部１１１と、ＣＦＭ  ＤＳＣ部１１０と、断層画像ＤＳＣ部１１２と、画像合成部１１３と、表示部１１４とを備える。これらの構成のうち、探触子１０１および表示部１１４には汎用の探触子および表示装置を用いることが可能であり、探触子１０１および表示部１１４を超音波診断装置１１が備えていなくてもよい。
【００３７】
超音波送受信部１０２は、探触子１０１を駆動する駆動信号を生成し、探触子１０１へ出力することによって、探触子１０１から被検体に向けて超音波を送信する。また、送信された超音波が被検体において反射することによって得られた反射エコーを探触子１０１によって受信し、受信信号を生成する。より具体的には、探触子１０１は複数の圧電素子を含み、各圧電素子から送信される超音波が超音波ビームを構成し、被検体を複数の超音波ビームによって走査するように、超音波送受信部１０２が各圧電素子の遅延制御を行いながら探触子１０１を駆動する。反射エコーは各圧電素子によって受信され、超音波送受信部１０２が、各圧電素子の遅延制御を行うことにより、送信された超音波ビームに対応する受信信号が生成する。被検体を超音波ビームが１回走査することによって１フレーム分のデータが得られる。１秒間に数回〜数十回繰り返して超音波の送受信を行うことにより、毎秒数フレームから数十フレームの受信信号が逐次生成する。
【００３８】
本実施形態の超音波診断装置１１は、Ｂモード断層画像とカラーフローマッピング画像と生成し、これらを合成して表示部１１４に表示する。このため、超音波送受信部１０２の上述した超音波の送受信は、Ｂモード断層画像の生成およびカラーフローマッピング画像の生成のそれぞれについて行われる。Ｂモード断層画像の１秒当たりのフレーム数およびカラーフローマッピング画像の１秒当たりのフレーム数は同じであってもよいし、異なっていてもよい。フレーム数が同じである場合には、Ｂモード断層画像生成のための超音波の送受信およびカラーフローマッピング画像生成のための超音波の送受信を交互に繰り返し行ってもよい。
【００３９】
Ｂモード断層画像を生成する場合、超音波送受信部１０２は、Ｂモード断層画像の生成に適した送受信を行い、得られた受信信号を断層画像信号処理部１１１へ出力する。カラーフローマッピング断層画像を生成する場合は、カラーフローマッピング断層画像の生成に適した送受信を行い、得られた受信信号をＣＦＭ信号処理部１０３へ出力する。一般的に、カラーフローマッピング断層画像を生成する場合、安定したカラーフローマッピング断層画像を得るために、超音波送受信部１０２は、同じ音響線上において超音波の送受信を複数回行う。
【００４０】
ＣＦＭ信号処理部１０３は、受信信号に対して直交検波処理、ＭＴＩフィルタ処理および自己相関処理を行い、血流速度および血流パワーを算出したのち、システムあるいは音響ノイズを排除するノイズカット処理を行う。ＣＦＭフレームデータは少なくとも血流速度データを含む。この他、血流パワーデータや血流速度の分散データを含んでいてもよい。ＣＦＭ信号処理部１０３はこの処理を、各フレームを構成する受信信号ごとに逐次繰り返して行う。ＣＦＭ信号処理部１０３で生成したＣＦＭフレームデータはフレームごとに残像処理部１１５へ出力される。
【００４１】
残像処理部１１５は、パーシスタンス係数を用いて、フレームごとにＣＦＭフレームデータに対して残像処理を行う。本実施形態の超音波診断装置１１は、パーシスタンス係数を最新のフレームの血流速度に応じて決定する。つまり、パーシスタンス係数は一定ではなく、最新のフレームの血流速度に基づく動的な値である。これにより、血流速度に応じてパーシスタンス係数を変化させ、残像効果を調整することができる。ただし、血流を動画で表示するためには、パルスドプラ法を用いて超音波の送受信を行う必要があり、このため、計測できる血流速度はパルス繰り返し周波数（ＰＲＦ）による制限を受ける。その結果、血流速度に折り返しが生じ、血流速度を正確に評価することが困難となる。
【００４２】
本実施形態の超音波診断装置１１は折り返しが生じているかどうかを判定するため、最新のフレームの血流速度データと１つ前のフレームの血流速度データを用いる。このために、残像処理部１１５は、フレームメモリ部（第１のメモリ部）１０４と、折り返し判定部１０５と、パーシスタンス係数決定部１０６と、パーシスタンス係数参照メモリ部（第３のメモリ部）１０７とパーシスタンス演算部１０８とパーシスタンスメモリ部（第２のメモリ部）１０９とを含む。
【００４３】
フレームメモリ部１０４は、最新のフレームの（現在走査中の）ＣＦＭフレームデータを記憶する。パーシスタンスメモリ部１０９は、最新より１つ前のフレームのパーシスタンス演算部１０８の出力結果であるＣＦＭフレームデータを記憶する。パーシスタンスメモリ部１０９のＣＦＭフレームデータは残像処理が施されている。以下、フレームメモリ部１０４およびパーシスタンスメモリ部１０９にそれぞれ記憶されたＣＦＭフレームデータのうち、血流速度データをＶcurrentおよび血流速度データＶout-1と呼ぶ。
【００４４】
折り返し判定部１０５は、フレームメモリ部１０４からＣＦＭフレームデータのうちの血流速度データＶcurrentとパーシスタンスメモリ部１０９からＣＦＭフレームデータのうちの血流速度データＶout-1を読み出し、折り返し判定を行う。より具体的には、血流速度データＶcurrentおよび血流速度データＶout-1と複数の閾値とを比較することによって、折り返しが発生しているかどうか、および、血流速度データＶcurrentが折り返し領域にあるかどうかを判定し、結果をパーシスタンス係数決定部１０６およびパーシスタンス演算部１０８に出力する。
【００４５】
パーシスタンス係数決定部１０６は、折り返し判定部１０５からの２つの判定結果とフレームメモリ部１０４から読み出した血流速度データＶcurrentに基づいて、パーシスタンス係数参照メモリ部１０７への参照インデックスを作成する。また、パーシスタンス係数参照メモリ部１０７にアクセスし、参照インデックスに対応付けられたパーシスタンス係数を読み出し、パーシスタンス演算部１０８に設定する。パーシスタンス係数参照メモリ部１０７には、あらかじめ血流速度の値に対応付けられたパーシスタンス係数の参照テーブルが記憶されている。この参照テーブルは、血流速度の値に応じて対応付けられた異なる２つ以上の値のパーシスタンス係数を含む。
【００４６】
パーシスタンス演算部１０８は、パーシスタンス係数決定部１０６から設定されたパーシスタンス係数と折り返し判定部１０５からの折り返し発生の判定結果とに基づいて、以下に示す式（１）で血流速度データに対してパーシスタンス演算を行う。パーシスタンス演算によって求められる残像処理された血流速度データをＶoutとし、パーシスタンス係数をＣpersistence（０＜Ｃpersistence＜１）とすると、残像処理された血流速度データは以下の式（１）で求められる。
Ｖout＝（１−Ｃpersistence）×Ｖcurrent＋Ｃpersistence×Ｖout-1 ・・（１)
【００４７】
ＣＦＭフレームデータが、血流速度データ以外のデータを含む場合には、同様に最新のフレームのデータおよび最新より１つ前のフレームのデータと求めたパーシスタンス係数Ｃpersistenceを用いてパーシスタンス演算を行い、残像処理されたデータを得る。
【００４８】
折り返し判定部１０５からの折り返し発生の判定結果が真であった場合は、式（１）の演算式を符号なし演算で扱い、偽であった場合は、符号付演算で扱う。
【００４９】
上述したように測定にパルス波を用いるため、ドップラ偏移によって直接測定できる血流速度はパルス波の繰り返し周波数（ＰＲＦ）の制限を受ける。具体的には、±ＰＲＦ／２を超える周波数変化に対応する血流速度は反対向きの血流として観測される折り返しが発生する。
【００５０】
図２（ａ）および（ｂ）は、残像処理された血流速度データＶout、最新のフレームの血流速度データＶcurrentおよび最新より１つ前のフレームのパーシスタンス演算部１０８の出力結果である血流速度データＶout-1の大きさの関係を示している。図２（ａ）および（ｂ）において、横軸の第１象限部分は速度Ｖがゼロを意味し、横軸の第２象限部分は、＋Ｖまたは−Ｖを意味する。速度Ｖが正である場合、第１または第２象限に位置し、速度Ｖが負である場合には、第３または第４象限に位置する。
【００５１】
例えば、図２（ａ）に示すように、Ｖcurrentが第２象限にあり、Ｖout-1が第３象限にあり、折り返しが発生していると判定された場合、Ｖout-1は実際には、＋ＰＲＦ／２に対応する血流速度よりも大きな値であったことなるから、ゼロを通らない、つまり、符号変化をともなわない演算である。したがって、ＶcurrentおよびＶout-1の符号（プラスまたはマイナス）をとって、式（１）にこれらの値を代入し演算を行う。
【００５２】
一方、例えば、図２（ｂ）に示すように、Ｖcurrentが第１象限にあり、Ｖout-1が第４象限にあり、折り返しが発生していないと判定された場合、式（１）の演算は、ゼロを通り、符号変化が生じる演算となる。したがって、符号をつけたＶcurrentおよびＶout-1を式（１）に代入し演算を行う。この演算は、１フレーム分の血流速度データの各画素または測定点ごとに行われる。また、演算結果であるＶoutは、折り返しが発生している場合、符号なしの値となる。この場合、血流速度データＶoutの最上位ビットを符号として扱うことで符号あり値としてＣＦＭ  ＤＳＣ部１１０およびパーシスタンスメモリ部１０９に出力される。
【００５３】
ＣＦＭ  ＤＳＣ部１１０は、パーシスタンス演算部１０８から出力される血流速度データの座標を変換し、画像合成部１１３へ出力する。
【００５４】
断層画像信号処理部４０９は受信信号に対してダイナミックフィルタ処理を施すことにより、不要なノイズをカットしたのち、包絡線検波処理およびダイナミックレンジ圧縮処理を施して断層画像ＤＳＣ部４１０に断層画像フレームデータを出力する。断層画像ＤＳＣ部４１０は、断層画像信号処理部４０９からの断層画像フレームデータの座標を変換し、画像合成部４１１へ出力する。
【００５５】
画像合成部４１１は、ＣＦＭ  ＤＳＣ部４１０および断層画像ＤＳＣ部４１０から出力される各フレームデータを画素ごとに、あるいは、対応する測定点のデータごとに合成し、合成画像フレームデータを生成する。具体的には、血流速度がゼロである場合は、断層画像フレームデータを表示し、そうでない場合は、ＣＦＭフレームデータを表示するように２つのデータを画素ごとに、あるいは、対応する測定点のデータごとに合成する。また、血流速度や血流の方向に応じてデータを色情報に変換し、表示部４１２に出力する。表示部４１２は、画像合成部４１１から受け取ったデータを表示する。
【００５６】
次に、パーシスタンス係数の決定についてより詳細に説明する。パーシスタンス係数を決定するために、まず、血流速度に折り返しが発生しているかどうかを折り返し判定部１０５において判定する。
【００５７】
折り返し判定部１０５は、フレームメモリ部１０４から最新のＣＦＭフレームデータに含まれる血流速度データＶcurrentとパーシスタンスメモリ部１０９から１フレーム前のパーシスタンス演算部１０８の出力結果であるＣＦＭフレームデータに含まれる血流速度Ｖout-1とを読み出し、ＶcurrentおよびＶout-1の値から以下の２つの判定を行う。
１．折り返しが発生しているかどうか。
２．Ｖcurrentが折り返し領域にあるかどうか。
【００５８】
この２つの状態の判定には、あらかじめ定められた閾値とＶcurrentおよびＶout-1とを比較することにより行う。具体的には、閾値Ｖthおよび血流ゼロ速度ＶzeroとＶcurrentおよびＶout-1とを比較する。
【００５９】
【表１】
2011036891.tif 000003
【００６０】
図３（ａ）は、閾値Ｖth、血流ゼロ速度Ｖzero、ＶcurrentおよびＶout-1の大きさの関係を示している。図３（ａ）において、横軸の第１象限部分は血流ゼロ速度Ｖzeroを意味し、横軸の第２象限部分は、Ｖmaxまたは−Ｖmaxを意味する。速度Ｖが正である場合、第１または第２象限に位置し、速度Ｖが負である場合には、第３または第４象限に位置する。
【００６１】
ここで、Ｖthおよび−Ｖthには、例えば、隣接するフレームの時間間隔において、想定される血流速度の変化の最大値を設定する。
【００６２】
表１は、折り返し判定部１０５において判定される条件と判定結果とを示している。
【００６３】
条件（０）に示すように、Ｖout-1が正である場合、想定される血流速度の変化の最大値はＶthまたは−Ｖthであるから、Ｖcurrentが−Ｖthよりも小さくなることはない。したがって、Ｖcurrent＜−Ｖthを満たしていれば、Ｖcurrentは、実際には、＋ＰＲＦ／２に対応する最大血流速度Ｖmaxよりも大きな値となっており、折り返しが発生しているとともに、Ｖcurrentが折り返し領域にあると判定される。条件（１）は、条件（０）の符号が逆転した場合である。
【００６４】
条件（２）に示すように、Ｖout-1が−Ｖthより小さい場合、Ｖcurrentが正の値となることは、想定される血流速度の変化の最大値を超える変化であるから、折り返しが発生している。また、ＶcurrentがＶzeroを挟む±Ｖthの範囲にあるため、Ｖcurrentは折り返し領域ではない。条件（３）は、条件（２）の符号が逆転した場合である。
【００６５】
条件（０）から（３）のいずれも満たさない場合には、折り返しが発生しておらず、また、Ｖcurrentが折り返し領域にないと判定する。
【００６６】
パーシスタンス係数決定部１０６は、折り返し判定部１０５から出力される２つの判定結果とフレームメモリ部１０４から読み出した血流速度データＶcurrentの絶対値に基づいて、パーシスタンス係数参照メモリ部１０７への参照インデックスを作成する。表２に作成される参照インデックスを示す。
【００６７】
【表２】
2011036891.tif 000004
【００６８】
折り返しが発生しており、かつ、Ｖcurrentが折り返し領域にある場合、血流速度Ｖcurrentは、実際には、Ｖmaxをまたは−Ｖmaxを超える大きな値と考えられる。このため、参照インデックスはＶmaxとなる。その他の場合には、Ｖcurrentの絶対値Abs(Ｖcurrent)となる。
【００６９】
パーシスタンス係数参照メモリ部１０７には、参照インデックスに対応付けられたパーシスタンス係数によって構成される参照テーブルが記憶されている。パーシスタンス係数決定部１０６は、パーシスタンス係数参照メモリ部１０７にアクセスし、作成された参照インデックスに対応付けられたパーシスタンス係数を読み出し、パーシスタンス演算部１０８へ出力する。
【００７０】
図３（ｂ）は、参照インデックスとパーシスタンス係数との対応関係の一例を示すグラフである。図３（ｂ）において、横軸は参照インデックスを示し、縦軸は、パーシスタンス係数を示している。表２に示したように、参照インデックスは、ＶmaxまたはＶcurrentの絶対値Abs(Ｖcurrent)である。Ｖcurrentの絶対値が閾値Ｖth以下である場合には、Ｖcurrentの増大にともなって単調に減少するパーシスタンス係数Ｃpersistenceが対応付けられている。つまり、Ｖcurrentの絶対値が閾値Ｖth以下である場合には、最新のフレームの血流速度Ｖcurrentに応じて異なるパーシスタンス係数Ｃpersistenceが対応付けられる。これにより、最新のフレームの血流速度Ｖcurrentが小さい場合には、パーシスタンス係数Ｃpersistenceが大きくなる。つまり、１つ前のフレームの血流速度Ｖout-1の重みが大きくなる。その結果、最新のフレームの血流速度Ｖcurrentが小さい場合には、１つ前のフレームの血流速度Ｖout-1を大きく反映させた血流速度Ｖoutが決定され表示部１１４に表示される。このため、カラーフローマッピング画像の変化が滑らかとなり、黒抜けが生じにくくなる。
【００７１】
また、最新のフレームの血流速度Ｖcurrentが大きい場合には、パーシスタンス係数Ｃpersistenceが小さくなる。つまり、１つ前のフレームの血流速度Ｖout-1の重みが小さくなる。その結果、最新のフレームの血流速度Ｖcurrentが大きい場合には、１つ前のフレームの血流速度Ｖout-1の影響が小さくなり、急激な血流速度の増大をリアルタイムで反映させたカラーフローマッピング画像を実現できる。
【００７２】
また、Ｖcurrentの増大にともなってパーシスタンス係数Ｃpersistenceが単調に減少するため、時間の経過とともに血流速度が増大する場合には、パーシスタンス係数Ｃpersistenceが減少し、残像効果が小さくなり、カラーフローマッピング画像の変化が急激になる。時間の経過とともに血流速度が減少する場合には、パーシスタンス係数Ｃpersistenceが増大し、残像効果が大きくなり、カラーフローマッピング画像の変化が緩やかになる
。
【００７３】
なお、表１および表２から分かるように、Ｖcurrent＜−Ｖthでも、Ｖout-1＞０であれば、参照インデックスはＶmaxとなる（条件（０））一方、Ｖout-1＜０であれば、参照インデックスはＶcurrentの絶対値Abs(Ｖcurrent)となる（条件（４））。このため、Ｖcurrent＜−Ｖthを満たしていても、Ｖout-1が正か負であるかによって、参照インデックスが異なり、パーシスタンス係数Ｃpersistenceも異なる。その結果、Ｖcurrent＜−Ｖthを満たす隣接する領域であっても、Ｖout-1が正か負であるかによってカラーマッピングフロー画像として表示される色が異なり、画像に不連続な色調部分が生じてしまう。
【００７４】
このような不自然な表示を抑制するため、Ｖcurrentの絶対値が閾値Ｖth以上である場合には、参照インデックスに同じ値のパーシスタンス係数Ｃpersistenceが対応付けられていることが好ましい。これにより、折り返しが発生する血流領域、あるいはその境界近傍において、自然な表示を行うことができる。
【００７５】
このように本実施形態の超音波診断装置によれば、ＣＦＭフレームデータに対して血流速度および折り返し状態に基づいてパーシスタンス係数を動的に決定したのちにパーシスタンス演算を実施することにより、頚動脈のような血流変化が激しい診断部位においても血流変化が明確に認識でき、かつ、低速な血流速度でも黒抜けが発生しない滑らかな血流動画を表示することが可能となる。
【００７６】
なお、上記実施形態では、ＣＦＭフレームデータの血流速度に基づいてパーシスタンス係数を動的に決定し、血流速度にパーシスタンス演算を行っているが、上述したようにＣＦＭフレームデータの他のデータ、例えば、血流パワーデータに対してパーシスタンス演算を行ってもよいし、Ｂモード断層画像データに対してパーシスタンス演算を行ってもよい。
【００７７】
また、上記実施形態では、最新のフレームおよび１つ前のフレームの血流速度データを用いてパーシスタンス処理を行っていたが、２つ前、あるいは３つ以上前のフレームの血流速度データも用いてパーシスタンス処理を行ってもよい。また、式（１）に限られず、他の演算式を用いてパーシスタンス処理を行ってもよい。
【００７８】
（第２の実施形態）
以下、図面を参照しながら本発明による超音波診断装置の第２の実施形態を説明する。図４は、本発明による超音波診断装置の一実施形態を示すブロック図である。図４に示す超音波診断装置１２は、探触子１０１と、超音波送受信部１０２と、ＣＦＭ信号処理部１０３と、残像処理部１１５’と、断層画像信号処理部１１１と、ＣＦＭ  ＤＳＣ部１１０と、断層画像ＤＳＣ部１１２と、画像合成部１１３と、表示部１１４とを備える。これらの構成のうち、探触子１０１および表示部１１４には汎用の探触子および表示装置を用いることが可能であり、探触子１０１および表示部１１４を超音波診断装置１２が備えていなくてもよい。
【００７９】
第１の実施形態で説明したように、超音波送受信部１０２は、探触子１０１を駆動する駆動信号を生成し、探触子１０１へ出力することによって、探触子１０１から被検体に向けて超音波を送信する。また、送信された超音波が被検体において反射することによって得られた反射エコーを探触子１０１によって受信し、受信信号を生成する。より具体的には、探触子１０１は複数の圧電素子を含み、各圧電素子から送信される超音波が超音波ビームを構成し、被検体を複数の超音波ビームによって走査するように、超音波送受信部１０２が各圧電素子の遅延制御を行いながら探触子１０１を駆動する。反射エコーは各圧電素子によって受信され、超音波送受信部１０２が、各圧電素子の遅延制御を行うことにより、送信された超音波ビームに対応する受信信号が生成する。被検体を超音波ビームが１回走査することによって１フレーム分のデータが得られる。１秒間に数回〜数十回繰り返して超音波の送受信を行うことにより、毎秒数フレームから数十フレームの受信信号が逐次生成する。
【００８０】
本実施形態の超音波診断装置１２は、Ｂモード断層画像とカラーフローマッピング画像と生成し、これらを合成して表示部１１４に表示する。このため、超音波送受信部１０２の上述した超音波の送受信は、Ｂモード断層画像の生成およびカラーフローマッピング画像の生成のそれぞれについて行われる。Ｂモード断層画像の１秒当たりのフレーム数およびカラーフローマッピング画像の１秒当たりのフレーム数は同じであってもよいし、異なっていてもよい。フレーム数が同じである場合には、Ｂモード断層画像生成のための超音波の送受信およびカラーフローマッピング画像生成のための超音波の送受信を交互に繰り返し行ってもよい。
【００８１】
Ｂモード断層画像を生成する場合、超音波送受信部１０２は、Ｂモード断層画像の生成に適した送受信を行い、得られた受信信号を断層画像信号処理部１１１へ出力する。カラーフローマッピング断層画像を生成する場合は、カラーフローマッピング断層画像の生成に適した送受信を行い、得られた受信信号をＣＦＭ信号処理部１０３へ出力する。一般的に、カラーフローマッピング断層画像を生成する場合、安定したカラーフローマッピング断層画像を得るために、超音波送受信部１０２は、同じ音響線上において超音波の送受信を複数回行う。
【００８２】
ＣＦＭ信号処理部１０３は、受信信号に対して直交検波処理、ＭＴＩフィルタ処理および自己相関処理を行い、血流速度および血流パワーを算出したのち、システムあるいは音響ノイズを排除するノイズカット処理を行う。ＣＦＭフレームデータは少なくとも血流速度データを含む。この他、血流パワーデータや血流速度の分散データを含んでいてもよい。ＣＦＭ信号処理部１０３はこの処理を、各フレームを構成する受信信号ごとに逐次繰り返して行う。ＣＦＭ信号処理部１０３で生成したＣＦＭフレームデータはフレームごとに残像処理部１１５’へ出力される。
【００８３】
残像処理部１１５’は、パーシスタンス係数を用いてフレームごとにＣＦＭフレームデータに対して残像処理を行う。本実施形態の超音波診断装置１２は、パーシスタンス係数を血流速度に応じてパーシスタンス係数を決定する。つまり、パーシスタンス係数は一定ではなく、血流速度に応じた動的な値である。これにより、血流速度に応じてパーシスタンス係数を変化させ、残像効果を調整することができる。ただし、血流を動画で表示するためには、パルスドプラ法を用いて超音波の送受信を行う必要があり、このため、計測できる血流速度はパルス繰り返し周波数（ＰＲＦ）による制限を受ける。その結果、血流速度に折り返しが生じ、血流速度を正確に評価することが困難となる。
【００８４】
本実施形態の超音波診断装置１２は折り返しが生じているかどうかを判定するため、最新のフレームの血流速度データと１つまえのフレームの血流速度データを用いる。また、残像処理部１１５’は２つパーシスタンス演算部を備え、残像効果をあまり与えずに速やかに血流速度を変化させる第１のパーシスタンス演算と残像効果を強力に与え血流速度変化をなるべく維持する第２のパーシスタンス演算を同時に行う。このようにして生成した残像効果の異なる２つの血流速度データのうち絶対値の大きい方を用いて血流画像を構成する。これにより、血流パワーが小さい末梢血管でも、血流表示を点滅させずに、かつ、平滑化に起因する末梢血管の血流消滅を起こさない動画を表示することが可能となる。
【００８５】
このために、残像処理部１１５’は、フレームメモリ部（第１のメモリ部）１０４と、折り返し判定部１０５と、第１のパーシスタンス係数決定部１０６Ａと、第１のパーシスタンス係数参照メモリ部（第３のメモリ部）１０７Ａと第１のパーシスタンス演算部１０８Ａと、第２のパーシスタンス係数決定部１０６Ｂと、第２のパーシスタンス係数参照メモリ部（第４のメモリ部）１０７Ｂと第２のパーシスタンス演算部１０８Ｂと、最大値選択部１１６と、パーシスタンスメモリ部（第２のメモリ部）１０９とを含む。
【００８６】
フレームメモリ部１０４は、最新のフレームの（現在走査中の）ＣＦＭフレームデータを記憶する。パーシスタンスメモリ部１０９は、最新より１つ前の最大値選択部の出力結果であるＣＦＭフレームデータを記憶する。パーシスタンスメモリ部１０９のＣＦＭフレームデータは残像処理が施されている。第１の実施形態と同様、フレームメモリ部１０４およびパーシスタンスメモリ部１０９にそれぞれ記憶されたＣＦＭフレームデータのうち、血流速度データをＶcurrentおよび血流速度データＶout-1と呼ぶ。
【００８７】
折り返し判定部１０５は、フレームメモリ部１０４からＣＦＭフレームデータのうちの血流速度データＶcurrentとパーシスタンスメモリ部１０９からＣＦＭフレームデータのうちの血流速度データＶout-1を読み出し、折り返し判定を行う。より具体的には、血流速度データＶcurrentおよび血流速度データＶout-1と複数の閾値とを比較することによって、折り返しが発生しているかどうか、および、血流速度データＶcurrentが折り返し領域にあるかどうかを判定し、結果を第１のパーシスタンス係数決定部１０６Ａ、第２のパーシスタンス係数決定部１０６Ｂ、第１のパーシスタンス演算部１０８Ａおよび第２のパーシスタンス演算部１０８Ｂに出力する。
【００８８】
第１のパーシスタンス係数決定部１０６Ａは、折り返し判定部１０５からの２つの判定結果とフレームメモリ部１０４から読み出した血流速度データＶcurrentに基づいて、第１のパーシスタンス係数参照メモリ部１０７Ａへの参照インデックスを作成する。また、第１のパーシスタンス係数参照メモリ部１０７Ａにアクセスし、参照インデックスに対応付けられた第１のパーシスタンス係数を読み出し、第１のパーシスタンス演算部１０８Ａに設定する。第１のパーシスタンス係数参照メモリ部１０７Ａには、あらかじめ血流速度の値に対応付けられた第１のパーシスタンス係数を含む第１の参照テーブルが記憶されている。この第１の参照テーブルは、血流速度の値に応じて対応付けられた異なる２つ以上の値のパーシスタンス係数を含む。
【００８９】
これに対し、第２のパーシスタンス係数決定部１０６Ｂは、折り返し判定部１０５からの２つの判定結果とパーシスタンスメモリ部１０９から読み出した血流速度データＶout-1に基づいて、第２のパーシスタンス係数参照メモリ部１０７Ｂへの参照インデックスを作成する。また、第２のパーシスタンス係数参照メモリ部１０７Ｂにアクセスし、参照インデックスに対応付けられた第２のパーシスタンス係数を読み出し、第２のパーシスタンス演算部１０８Ｂに設定する。第２のパーシスタンス係数参照メモリ部１０７Ｂには、あらかじめ血流速度の値に対応付けられた第２のパーシスタンス係数を含む第２の参照テーブルが記憶されている。第２の参照テーブルも、血流速度の値に応じて対応付けられた異なる２つ以上の値のパーシスタンス係数を含むが、以下において詳細に説明するように、同じ血流速度値に対して対応付けられる第１のパーシスタンス係数値と第２のパーシスタンス係数値とは値が異なっている。
【００９０】
第１のパーシスタンス演算部１０８Ａは、第１のパーシスタンス係数決定部１０６Ａから設定されたパーシスタンス係数と折り返し判定部１０５からの折り返し発生の判定結果とに基づいて、以下に示す式（１）で血流速度データに対してパーシスタンス演算を行う。
【００９１】
パーシスタンス演算によって求められる残像処理された血流速度データをＶoutとし、パーシスタンス係数をＣpersistence（０＜Ｃpersistence＜１）とすると、残像処理された血流速度データは以下の式（１）で求められる。
Ｖout＝（１−Ｃpersistence）×Ｖcurrent＋Ｃpersistence×Ｖout-1 ・・（１)
【００９２】
同様に第２のパーシスタンス演算部１０８Ｂも、第２のパーシスタンス係数決定部１０６Ｂから設定されたパーシスタンス係数と折り返し判定部１０５からの折り返し発生の判定結果とに基づいて、式（１）で血流速度データに対してパーシスタンス演算を行う。
【００９３】
第１のパーシスタンス演算部１０８Ａおよび第２のパーシスタンス演算部１０８Ｂにおける演算は、決定したパーシスタンス係数が互い異なる点を除いて同じである。ＣＦＭフレームデータが、血流速度データ以外のデータを含む場合には、同様に最新のフレームのデータおよび最新より１つ前のフレームのデータと求めたパーシスタンス係数Ｃpersistenceを用いて第１のパーシスタンス演算部１０８Ａおよび第２のパーシスタンス演算部１０８Ｂでパーシスタンス演算を行い、残像処理されたデータをそれぞれ得る。
【００９４】
折り返し判定部１０５からの折り返し発生の判定結果が真であった場合は、式（１）の演算式を符号なし演算で扱い、偽であった場合は、符号付演算で扱う。
【００９５】
上述したように測定にパルス波を用いるため、ドップラ偏移によって直接測定できる血流速度はパルス波の繰り返し周波数（ＰＲＦ）の制限を受ける。具体的には、±ＰＲＦ／２を超える周波数変化に対応する血流速度は反対向きの血流として観測される折り返しが発生する。
【００９６】
図５（ａ）および（ｂ）は、残像処理された血流速度データＶout、最新のフレームの血流速度データＶcurrentおよび最新より１つ前のフレームのパーシスタンス演算部１０８の出力結果である血流速度データＶout-1の大きさの関係を示している。図５（ａ）および（ｂ）において、横軸の第１象限部分は速度Ｖがゼロを意味し、横軸の第２象限部分は、＋Ｖまたは−Ｖを意味する。速度Ｖが正である場合、第１または第２象限に位置し、速度Ｖが負である場合には、第３または第４象限に位置する。
【００９７】
例えば、図５（ａ）に示すように、Ｖcurrentが第２象限にあり、Ｖout-1が第３象限にあり、折り返しが発生していると判定された場合、Ｖout-1は実際には、＋ＰＲＦ／２に対応する血流速度よりも大きな値であったことなるから、ゼロをとらない、つまり、符号変化を伴わない演算である。したがって、ＶcurrentおよびＶout-1の符号（プラスまたはマイナス）をとって、式（１）にこれらの値を代入し演算を行う。
【００９８】
一方、例えば、図５（ｂ）に示すように、Ｖcurrentが第１象限にあり、Ｖout-1が第４象限にあり、折り返しが発生していないと判定された場合、式（１）の演算は、ゼロを通り、符号変化が生じる演算となる。したがって、符号をつけたＶcurrentおよびＶout-1を式（１）に代入し演算を行う。この演算は、１フレーム分の血流速度データの各画素または測定点ごとに行われる。また、演算結果であるＶoutは、折り返しが発生している場合、符号なしの値となる。この場合、血流速度データＶoutの最上位ビットを符号として扱うことで符号あり値として最大値選択部１１６にそれぞれ出力される。
【００９９】
最大値選択部１１６は、第１のパーシスタンス演算部１０８Ａおよび第２のパーシスタンス演算部１０８Ｂからそれぞれ演算結果、つまり、残像処理された血流速度データを受け取り、画素ごとに、あるいは、対応する測定点のデータごとに血流速度の絶対値を比較し、大きい方の血流速度を選択して、最新のフレームの残像処理された血流速度データを構成し、これをＣＦＭ  ＤＳＣ部１１０およびパーシスタンスメモリ部１０９に出力する。ＣＦＭ  ＤＳＣ部１１０は、選択した血流速度データの座標を変換し、画像合成部１１３へ出力する。
【０１００】
断層画像信号処理部４０９は受信信号に対してダイナミックフィルタ処理を施すことにより、不要なノイズをカットしたのち、包絡線検波処理およびダイナミックレンジ圧縮処理を施して断層画像ＤＳＣ部４１０に断層画像フレームデータを出力する。断層画像ＤＳＣ部４１０は、断層画像信号処理部４０９からの断層画像フレームデータの座標を変換し、画像合成部４１１へ出力する。
【０１０１】
画像合成部４１１は、ＣＦＭ  ＤＳＣ部４１０および断層画像ＤＳＣ部４１０から出力される各フレームデータを画素ごとに、あるいは、対応する測定点のデータごとに合成し、合成画像フレームデータを生成する。具体的には、血流速度がゼロである場合は、断層画像フレームデータを表示し、そうでない場合は、ＣＦＭフレームデータを表示するように２つのデータを画素ごとに、あるいは、対応する測定点のデータごとに合成する。また、血流速度や血流の方向に応じてデータを色情報に変換し、表示部４１２に出力する。表示部４１２は、画像合成部４１１から受け取ったデータを表示する。
【０１０２】
次に、第１および第２のパーシスタンス係数の決定についてより詳細に説明する。第１および第２のパーシスタンス係数を決定するために、まず、血流速度に折り返しが発生しているかどうかを折り返し判定部１０５において判定する。
【０１０３】
折り返し判定部１０５は、フレームメモリ部１０４から最新のＣＦＭフレームデータに含まれる血流速度データＶcurrentとパーシスタンスメモリ部１０９から１フレーム前のパーシスタンス演算部１０８の出力結果であるＣＦＭフレームデータに含まれる血流速度Ｖout-1とを読み出し、ＶcurrentおよびＶout-1の値から以下の２つの判定を行う。
１．折り返しが発生しているかどうか。
２．Ｖcurrentが折り返し領域にあるかどうか。
【０１０４】
この２つの状態の判定には、あらかじめ定められた閾値とＶcurrentおよびＶout-1とを比較することにより行う。具体的には、閾値Ｖthおよび血流ゼロ速度ＶzeroとＶcurrentおよびＶout-1とを比較する。
【０１０５】
【表３】
2011036891.tif 000005
【０１０６】
図６（ａ）は、閾値Ｖth、血流ゼロ速度Ｖzero、ＶcurrentおよびＶout-1の大きさの関係を示している。図６（ａ）において、横軸の第１象限部分は血流ゼロ速度Ｖzeroを意味し、横軸の第２象限部分は、Ｖmaxまたは−Ｖmaxを意味する。速度Ｖが正である場合、第１または第２象限に位置し、速度Ｖが負である場合には、第３または第４象限に位置する。
【０１０７】
ここで、Ｖthおよび−Ｖthには、例えば、隣接するフレームの時間間隔において、想定される血流速度の変化の最大値を設定する。
【０１０８】
表３は、折り返し判定部１０５において判定される条件と判定結果とを示している。
【０１０９】
条件（０）に示すように、Ｖout-1が正である場合、想定される血流速度の変化の最大値はＶthまたは−Ｖthであるから、Ｖcurrentが−Ｖthよりも小さくなることはない。したがって、Ｖcurrent＜−Ｖthを満たしていれば、Ｖcurrentは、実際には、＋ＰＲＦ／２に対応する最大血流速度Ｖmaxよりも大きな値となっており、折り返しが発生しているとともに、Ｖcurrentが折り返し領域にあると判定される。条件（１）は、条件（０）の符号が逆転した場合である。
【０１１０】
条件（２）に示すように、Ｖout-1が−Ｖthより小さい場合、Ｖcurrentが正の値となることは、想定される血流速度の変化の最大値を超える変化であるから、折り返しが発生している。また、ＶcurrentがＶzeroを挟む±Ｖthの範囲にあるため、Ｖcurrentは折り返し領域ではない。条件（３）は、条件（２）の符号が逆転した場合である。
【０１１１】
条件（０）から（３）のいずれも満たさない場合には、折り返しが発生しておらず、また、Ｖcurrentが折り返し領域にないと判定する。
【０１１２】
第１のパーシスタンス係数決定部１０６Ａは、折り返し判定部１０５から出力される２つの判定結果とフレームメモリ部１０４から読み出した血流速度データＶcurrentの絶対値に基づいて、第１のパーシスタンス係数参照メモリ部１０７への参照インデックスを作成する。表４に作成される参照インデックスを示す。
【０１１３】
【表４】
2011036891.tif 000006
【０１１４】
折り返しが発生しており、かつ、Ｖcurrentが折り返し領域にある場合、血流速度Ｖcurrentは、実際には、Ｖmaxをまたは−Ｖmaxを超える大きな値と考えられる。このため、参照インデックスはＶmaxとなる。その他の場合には、Ｖcurrentの絶対値Abs(Ｖcurrent)となる。
【０１１５】
第１のパーシスタンス係数参照メモリ部１０７Ａには、参照インデックスに対応付けられた第１のパーシスタンス係数によって構成される第１の参照テーブルが記憶されている。第１のパーシスタンス係数決定部１０６Ａは、第１のパーシスタンス係数参照メモリ部１０７Ａにアクセスし、作成された参照インデックスに対応付けられた第１のパーシスタンス係数を読み出し、第１のパーシスタンス演算部１０８Ａへ出力する。
【０１１６】
図６（ｂ）は、参照インデックスと第１のパーシスタンス係数との対応関係の一例を示すグラフである。図６（ｂ）において、横軸は参照インデックスを示し、縦軸は、パーシスタンス係数を示している。表４に示したように、参照インデックスは、ＶmaxまたはＶcurrentの絶対値Abs(Ｖcurrent)である。Ｖcurrentの絶対値が閾値Ｖth以下である場合には、Ｖcurrentの増大にともなって単調に増加する第１のパーシスタンス係数Ｃpersistenceが対応付けられている。つまり、Ｖcurrentの絶対値が閾値Ｖth以下である場合には、最新のフレームの血流速度Ｖcurrentに応じて異なるパーシスタンス係数Ｃpersistenceが対応付けられる。
【０１１７】
これに対し、第２のパーシスタンス係数決定部１０６Ｂは、折り返し判定部１０５から出力される２つの判定結果とパーシスタンスメモリ部１０９から読み出した血流速度データＶout-1の絶対値に基づいて、第１のパーシスタンス係数参照メモリ部１０７への参照インデックスを作成する。表５に作成される参照インデックスを示す。
【０１１８】
【表５】
2011036891.tif 000007
【０１１９】
第２のパーシスタンス係数決定部１０６Ｂは、条件（２）から（４）の場合に、パーシスタンスメモリ部１０９から読み出した血流速度データＶout-1の絶対値を参照インデックスとして生成する点で第１のパーシスタンス係数決定部１０６Ａと異なっている。
【０１２０】
第２のパーシスタンス係数参照メモリ部１０７Ｂには、参照インデックスに対応付けられた第２のパーシスタンス係数によって構成される第２の参照テーブルが記憶されている。第２のパーシスタンス係数決定部１０６Ｂは、第２のパーシスタンス係数参照メモリ部１０７Ｂにアクセスし、作成された参照インデックスに対応付けられた第２のパーシスタンス係数を読み出し、第２のパースタンス演算部１０８Ｂへ出力する。
【０１２１】
図６（ｃ）は、参照インデックスと第２のパーシスタンス係数との対応関係の一例を示すグラフである。図６（ｃ）において、横軸は参照インデックスを示し、縦軸は、パーシスタンス係数を示している。表５に示したように、参照インデックスは、ＶmaxまたはＶout-1の絶対値Abs(Ｖout-1)である。Ｖout-1の絶対値が閾値Ｖth以下である場合には、Ｖout-1の増大にともなって単調に増加する第２のパーシスタンス係数Ｃpersistenceが対応付けられている。つまり、Ｖout-1の絶対値が閾値Ｖth以下である場合には、１つ前の血流速度Ｖout-1に応じて異なる第２のパーシスタンス係数Ｃpersistenceが対応付けられる。
【０１２２】
図６（ｂ）および図６（ｃ）に示すように、参照インデックスがいずれの値をとっても第２のパーシスタンス係数のほうが第１のパーシスタンス係数より大きい。つまり、第１のパーススタンス係数は、最新のフレームの血流速度に対応づけられ、かつ、小さな値である。第１のパーススタンス係数が大きくなれば、１つ前のフレームの血流速度をより考慮した演算となるため、第１のパーシスタンス演算部１０８Ａは、残像効果を抑制し、すみやかに血流速度を変化させる演算を行う。これに対し、第２のパーススタンス係数は、１つ前のフレームの血流速度に対応付けられ、かつ、大きな値であるため、第２のパーシスタンス演算部１０８Ｂは、残像効果を高め、血流速度の変化を抑制する演算を行う。
【０１２３】
また、上述したように第１のパーシスタンス演算部１０８Ａは、残像効果を抑制し、速やかに血流速度を変化させる演算を行うため、血流速度は高いものの、血流パワーが小さいために血流を正しく検出できない場合、突然血流速度がゼロとなる可能性がある。この場合、血流速度が高くなるにつれて、血流画像に色調あるいは諧調の色付けを行っていると、血流画像が突然暗い諧調で色づけされ、画像が点滅するように表示される。このため、参照インデックスが増加するにつれて第１のパーシスタンス係数を単調に増加させ、血流速度が高くなるにつれて残像効果を高め、血流画像の点滅を抑制することができる。
【０１２４】
また、第２のパーシスタンス演算部１０８Ｂは、残像効果を高めた画像表示を行うため、血流速度が高くなるにつれて、血流画像に色調あるいは諧調の色付けを行っていると、血流速度が低い場合には、暗めの表示が必要以上に長い時間残像として表示されてしまう。例えば、探触子を動かした場合に、血流表示が尾引く印象を与えてしまう。このため、参照インデックスが増加するにつれて第２のパーシスタンス係数を単調に増加させ、血流速度が低くなるにつれて残像効果を抑制することができる。従って、血流速度の絶対値に基づく参照インデックスと第１および第２のパーシスタンス係数との間に適切な単調増加する関係を設定することにより、高品位な血流表示を実現することができる。
【０１２５】
なお、表３、表４および表５から分かるように、Ｖcurrent＜−Ｖthでも、Ｖout-1＞０であれば、参照インデックスはＶmaxとなる（条件（０））一方、Ｖout-1＜０であれば、参照インデックスはＶcurrentの絶対値Abs(Ｖcurrent)となる（条件（４））。このため、Ｖcurrent＜−Ｖthを満たしていても、Ｖout-1が正か負であるかによって、参照インデックスが異なり、パーシスタンス係数Ｃpersistenceも異なる。その結果、Ｖcurrent＜−Ｖthを満たす隣接する領域であっても、Ｖout-1が正か負であるかによってカラーマッピングフロー画像として表示される色が異なり、画像に不連続な色調部分が生じてしまう。
【０１２６】
このような不自然な表示を抑制するため、Ｖcurrentの絶対値が閾値Ｖth以上である場合には、参照インデックスに同じ値のパーシスタンス係数Ｃpersistenceが対応付けられていることが好ましい。これにより、折り返しが発生する血流領域、あるいはその境界近傍において、自然な表示を行うことができる。
【０１２７】
このようにして決定した第１および第２のパーシスタンス係数を用いて第１のパーシスタンス演算部１０８Ａおよび第２のパーシスタンス演算部１０８Ｂがそれぞれパーシスタンス処理の施された最新の血流速度データを生成する。
【０１２８】
最大値選択部１１６は、２つの血流速度データの絶対値の大きいほうを選択し、選択した血流速度データを残像処理された血流速度データとして出力する。つまり、絶対値の大きな血流速度データが得られるように、２つのパーシスタンス処理の結果を選択するため、甲状腺、肝臓、腎臓などに存在するような血流パワーが不安定となる末梢血管の血流表示を点滅させずに、かつ、平滑化に起因する末梢血管の血流消滅を起こさない血流動画を表示することが可能となる。
【０１２９】
なお、上記実施形態では、ＣＦＭフレームデータの血流速度に基づいてパーシスタンス係数を動的に決定し、血流速度にパーシスタンス演算を行っているが、上述したようにＣＦＭフレームデータの他のデータ、例えば、血流パワーデータに対してパーシスタンス演算を行ってもよいし、Ｂモード断層画像データに対してパーシスタンス演算を行ってもよい。
【０１３０】
また、上記実施形態では、最新のフレームおよび１つ前のフレームの血流速度データを用いてパーシスタンス処理を行っていたが、２つ前、あるいは３つ以上前のフレームの血流速度データも用いてパーシスタンス処理を行ってもよい。また、式（１）に限られず、他の演算式を用いてパーシスタンス処理を行ってもよい。
【産業上の利用可能性】
【０１３１】
本発明は、被検体の血流状態を表示することのできる超音波診断装置に好適に用いることができる。
【符号の説明】
【０１３２】
１０１、４０１  探触子
１０２、４０２  超音波送受部
１０３、４０３  ＣＦＭ信号処理部
１０４、４０４  フレームメモリ部
１０５  折り返し判定部
１０６  パーシスタンス係数決定部
１０６Ａ  第１のパーシスタンス係数決定部
１０６Ｂ  第２のパーシスタンス係数決定部
１０７  パーシスタンス係数参照メモリ部
１０７Ａ  第１のパーシスタンス係数参照メモリ部
１０７Ｂ  第２のパーシスタンス係数参照メモリ部
１０８、４０７  パーシスタンス演算部
１０８Ａ  第１のパーシスタンス演算部
１０８Ｂ  第２のパーシスタンス演算部
１０９  パーシスタンスメモリ部
１１０、４０８  ＣＦＭ  ＤＳＣ部
１１１、４０９  断層画像信号処理部
１１２、４１０  断層画像  ＤＳＣ部
１１３、４１１  画像合成部
１１４、４１２  表示部
１１５、１１５’  残像処理部
１１６  最大値選択部
４０５  フレームメモリ選択部
４０６  パーシスタンス係数設定部
【図１】
2011036891.tif 000008
【図２】
2011036891.tif 000009
【図３】
2011036891.tif 000010
【図４】
2011036891.tif 000011
【図５】
2011036891.tif 000012
【図６】
2011036891.tif 000013
【図７】
2011036891.tif 000014
【手続補正書】
【提出日】平成23年4月21日(2011.4.21)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波診断装置に関し、特にカラーフローマッピング時の残像処理方法に関する。
【背景技術】
【０００２】
超音波診断装置は、超音波を被検体に送信し、その反射エコーに含まれる情報を解析することにより、被検体内の画像を作成する。カラーフローマッピング（Color Flow Mapping、以下ＣＦＭと略す場合がある）と呼ばれる手法によって、被検体内の血流を画像化することも可能であり、医療分野全般において、血流状態を表示することのできる超音波診断装置が広く利用されている。
【０００３】
カラーフローマッピングはカラードップライメージング（Color Doppler Imaging、ＣＤＩ）とも呼ばれ、ドップラ効果を利用する。超音波が血流を照射すると、ドプラ効果により反射エコーに血流速度に応じたドップラ偏移が生じる。このドップラ偏移の情報を直交検波によって検出し、ＭＴＩ（Moving Target Indicator）フィルタと呼ばれるハイパスフィルタ処理、自己相関処理およびノイズカット処理を施すことによって、血流速度に関する情報が得られる。得られた血流速度の情報を色情報に変換し、二次元的にＢモード断層画像に重畳して表示することによって、被検体内の血流状態をユーザに認識させることができる。
【０００４】
血流から得られる反射エコーによる受信信号強度は、Ｂモード断層画像生成に用いる組織散乱体および組織境界から得られる反射エコーによる受信信号強度に比べ、かなり小さい。このため、カラーフローマッピングにおける信号処理によって得られる血流速度および血流パワー（移動する血流量）は不安定になりやすい。
【０００５】
特に、観察したい部分の血流速度が遅い場合、あるいは、観察したい部分が末梢血管である場合、血流パワーが小さくなるため、本来、システムノイズおよび音響ノイズのみをカットするべきノイズカット処理において、血流速度あるいは血流パワーに関する情報が除去されてしまい易い。その結果、血流画像において本来血流として表示される部分が黒く抜けてしまう現象が発生する。例えば、１秒間に数フレームから数十フレームの割合で、被検体内の血流を画像化する場合、そのうちのいくつかのフレームにおいて、血流部分が黒く示されてしまう。このため、断層画像中の血流部分が突然消滅し、画像が滑らかでなくなったり、違和感のあるものとなったりする。
【０００６】
この問題を解決するために、従来のカラーフローマッピングを行う超音波診断装置では、信号処理後段においてパーシスタンス処理（残像処理）と呼ばれる時間方向補間を行うのが通例である。以下、特許文献１に示される従来のカラーフローマッピングにおけるパーシスタンス処理を説明する。
【０００７】
図７に示す従来の超音波診断装置において、超音波送受信部４０２は、探触子４０１を駆動し、超音波を被検体へ送信する。また、探触子４０１によって被検体で生じる反射エコーを受信し、受信信号を生成する。Ｂモード断層画像を生成する場合、超音波送受信部４０２は、Ｂモード断層画像の生成に適した送受信を行い、得られた受信信号を断層画像信号処理部４０９へ出力する。カラーフローマッピング断層画像を生成する場合は、カラーフローマッピング断層画像の生成に適した送受信を行い、得られた受信信号をカラーフローマッピング処理部４０３（以下、ＣＦＭ信号処理部と略す）へ出力する。一般的に、カラーフローマッピング断層画像を生成する場合、安定したカラーフローマッピング断層画像を得るために、超音波送受信部４０２は、同じ音響線上において超音波の送受信を複数回行う。
【０００８】
ＣＦＭ信号処理部４０３は、受信信号に対して直交検波処理、ＭＴＩフィルタ処理および自己相関処理を行い、血流速度および血流パワーを算出したのち、システムあるいは音響ノイズを排除するノイズカット処理を行い、フレームメモリ部４０４に血流速度および血流パワーを出力する。
【０００９】
フレームメモリ部４０４は、リングバッファで構成されており、現在の走査フレームからＮフレーム前（Ｎは、１以上の整数）までの血流速度および血流パワーをフレーム単位で保存している。ここで、フレームとは、１画面のＣＦＭ断層画像を構成する血流速度データおよび血流パワーデータ群を示す。
【００１０】
フレームメモリ選択部４０５は、フレームメモリ部４０４から、あらかじめ設定された複数のＣＦＭフレームデータを選択し、パーシスタンス演算部４０７へ出力するようフレームメモリ部４０４に指令を出力する。パーシスタンス演算部４０７は、フレームメモリ部４０４から読み出されたＣＦＭフレームデータとパーシスタンス係数設定部４０６から出力されるパーシスタンス係数に基づいてパーシスタンス演算を実施し、ＣＦＭ  ＤＳＣ（Digital Scan Converter）部４０８に出力する。パーシスタンス演算は、単純な重み付け演算であり、パーシスタンス係数設定部４０６から出力されるパーシスタンス係数は、あらかじめシステムから設定された固定係数である。
【００１１】
ＣＦＭ  ＤＳＣ部４０８は、パーシスタンス演算部４０７から出力されるＣＦＭフレームデータの座標を変換し、画像合成部４１１へ出力する。
【００１２】
断層画像信号処理部４０９は受信信号に対してダイナミックフィルタ処理を施すことにより、不要なノイズをカットしたのち、包絡線検波処理およびダイナミックレンジ圧縮処理を施して断層画像ＤＳＣ部４１０に断層画像フレームデータを出力する。断層画像ＤＳＣ部４１０は、断層画像信号処理部４０９からの断層画像フレームデータの座標を変換し、画像合成部４１１へ出力する。
【００１３】
画像合成部４１１は、ＣＦＭ  ＤＳＣ部４１０および断層画像ＤＳＣ部４１０から出力される各フレームデータを画素ごとに合成し、合成画像フレームデータを生成する。具体的には、血流速度がゼロである場合は、断層画像フレームデータを表示し、そうでない場合は、ＣＦＭフレームデータを表示するように２つのデータを画素ごと、あるいは、対応する測定点のデータごとに合成する。また、血流速度や血流の方向に応じてデータを色情報に変換し、表示部４１２に出力する。表示部４１２は、画像合成部４１１から受け取ったデータを表示する。
【先行技術文献】
【特許文献】
【００１４】
【特許文献１】特開平２-２８６１４０号公報
【発明の概要】
【発明が解決しようとする課題】
【００１５】
従来の超音波診断装置におけるパーシスタンス処理は、血流速度が遅く、あるいは、血流パワーが小さく、これらが不安定な状態にあることに起因してＣＦＭ信号処理部４０３の出力結果が不安定になり、表示血流に黒抜けが発生することを防止する。具体的には、現在走査中のフレームデータよりも過去のフレームデータに重みを置いたパースタンス係数を用いることによって、残像効果を生じさせ、画像に黒抜けが発生するのを抑制する。
【００１６】
しかし、このような特徴は、血流速度が激しく変化する動脈の診断には適していない場合がある。例えば、頚動脈では、心臓の拡張収縮に応じて血流が激しく変化し、収縮期には、心周期に対してごく短い時間において血流が速くなり、拡張期では、血流が遅くなる。また、血流速度の最大値と最小値との差も他の診断部位に比べて大きい。その結果、拡張期における頚動脈の血流速度は、心周期に対して比較的長い時間、小さな値で一定となり、ＣＦＭ信号処理部４０３の出力が安定しない状態に陥る。
【００１７】
表示する血流の画像に黒抜けが生じるのを抑制するためには、パーシスタンス処理によって、残像効果が高くなるようにパーシスタンス係数を設定することが好ましい。これにより、血流速度が低速となる場合でも、黒抜けが発生しない滑らかな動画表示が可能となるからである。しかし、この場合、収縮期における高速な血流速度の表示ができなくなる。
【００１８】
また、従来の超音波診断装置におけるパーシスタンス処理は、末梢血管の診断には適していない場合がある。例えば、甲状腺、肝臓、腎臓などの臓器には、主流血管から派生する末梢血管が存在する。これらの臓器の診断にあたり、末梢血管構造を把握することは非常に重要である。
【００１９】
末梢血管における時間的な血流変化は、比較的安定しているのであるが、血管が物理的に細いために、血流パワーが頚動脈や心臓の場合と比べ、極端に小さくなる傾向にある。従って、血流パワーが小さいことにより、ドップラ偏移の検出が不安定となり、ＣＦＭ信号処理部４０３の出力が結果として安定しない状態に陥る。
【００２０】
このため、パーシスタンス処理を行わない場合、断層画像中の末梢血管は、時間的に点滅して表示されることになり、動画として見づらくなる。逆にパーシスタンス処理を行う場合、時間的方向に対して断層画像中の末梢血管の血流が平滑化されるため、パーシスタンス処理により末梢血管が消滅する可能性がある。この場合、末梢血管の検出率が著しく低下することになる。
【００２１】
本発明は、このような従来技術の課題を解決し、血流速度が激しく変化する頚動脈などのような診断部位において、血流変化が明確に認識でき、かつ、低速な血流速度でも黒抜けが発生しない滑らかな血流動画を表示することができる超音波診断装置を提供することを目的とする。また、末梢血管など血流パワーが小さい血管部分でも見やすい動画表示を行うことのできる超音波診断装置を提供することを目的とする。
【課題を解決するための手段】
【００２２】
本発明の超音波診断装置は、探触子を繰り返し駆動し、前記探触子を駆動することにより送信された超音波が被検体において反射することにより得られた反射エコーを前記探触子によって受信し、複数の受信信号を逐次生成する送受信部と、前記複数の受信信号に基づき、各フレームにおける前記被検体中の血流部分の血流速度データを逐次生成するカラーフローマッピング信号処理部と、前記各フレームにおける血流速度データに対して残像処理を行う残像処理部と、前記受信信号に基づき、Ｂモード断層画像フレームデータを生成する断層画像信号処理部と、前記残像処理された血流速度データおよび前記Ｂモード断層画像フレームデータを合成する画像合成部とを備え、前記残像処理部は、最新のフレームの血流速度データと最新より以前のフレームの残像処理された血流速度データに基づいて折り返し判定を行い、前記折り返し判定の結果と最新および最新より以前のフレームの血流速度データに基づいて動的にパーシスタンス係数を変更する。
【００２３】
ある好ましい実施形態の超音波診断装置は、探触子を繰り返し駆動し、前記探触子を駆動することにより送信された超音波が被検体において反射することにより得られた反射エコーを前記探触子によって受信し、複数の受信信号を逐次生成する送受信部と、前記複数の受信信号に基づき、各フレームにおける前記被検体中の血流部分の血流速度データを逐次生成するカラーフローマッピング信号処理部と、前記各フレームにおける血流速度データに対して残像処理を行う残像処理部と、前記受信信号に基づき、Ｂモード断層画像フレームデータを生成する断層画像信号処理部と、前記残像処理された血流速度データおよび前記Ｂモード断層画像フレームデータを合成する画像合成部とを備え前記残像処理部は、最新のフレームの血流速度データを記憶する第１のメモリ部と、最新より１つ前のフレームの残像処理された血流速度データを記憶する第２のメモリ部と、前記第１のメモリ部および前記第２のメモリ部から血流速度データをそれぞれ読み出し、折り返し判定を行う折り返し判定部と、前記折り返し判定の結果および前記第１のメモリ部に記憶された血流速度データに応じてパーシスタンス係数を決定するパーシスタンス係数決定部と、前記パーシスタンス係数および前記折り返し判定の結果に基づいて、前記第１のメモリ部に記憶された血流速度データに対しパーシスタンス演算を行い、演算結果を前記残像処理された血流速度データとして出力するパーシスタンス演算部とを含む。
【００２４】
ある好ましい実施形態において、前記折り返し判定部は、前記第１のメモリ部に記憶された血流速度データおよび前記第２のメモリ部に記憶された血流速度データと複数の閾値とを比較することによって、折り返しが発生しているかどうか、および、前記最新のフレームの血流速度データが折り返し領域にあるかどうかを判定する。
【００２５】
ある好ましい実施形態において、前記残像処理部は、前記血流速度の値に応じて対応付けられた異なる２つ以上の値のパーシスタンス係数を含む参照テーブルを記憶している第３のメモリ部をさらに含む。
【００２６】
ある好ましい実施形態において、前記参照テーブルは、所定の値以上の血流速度に対して一定の値のパーシスタンス係数が対応付けられている。
【００２７】
ある好ましい実施形態の超音波診断装置は、探触子を繰り返し駆動し、前記探触子を駆動することにより送信された超音波が被検体において反射することにより得られた反射エコーを前記探触子によって受信し、複数の受信信号を逐次生成する送受信部と、前記複数の受信信号に基づき、各フレームにおける前記被検体中の血流部分の血流速度データを逐次生成するカラーフローマッピング信号処理部と、前記各フレームにおける血流速度データに対して残像処理を行う残像処理部と、前記受信信号に基づき、Ｂモード断層画像フレームデータを生成する断層画像信号処理部と、前記残像処理された血流速度データおよび前記Ｂモード断層画像フレームデータを合成する画像合成部とを備え、前記残像処理部は、最新のフレームの血流速度データを記憶する第１のメモリ部と、最新より１つ前のフレームの残像処理された血流速度データを記憶する第２のメモリ部と、前記第１のメモリ部および前記第２のメモリ部から血流速度データをそれぞれ読み出し、折り返し判定を行う折り返し判定部と、前記折り返し判定の結果および前記第１のメモリ部に記憶された血流速度データに応じて第１のパーシスタンス係数を決定する第１のパーシスタンス係数決定部と、前記第１のパーシスタンス係数および前記折り返し判定の結果に基づいて前記第１のメモリ部に記憶された血流速度データに対しパーシスタンス演算を行う第１のパーシスタンス演算部と、前記折り返し判定の結果および前記第２のメモリ部に記憶された血流速度データに応じて第２のパーシスタンス係数を決定する第２のパーシスタンス係数決定部と、前記第２のパーシスタンス係数および前記折り返し判定の結果に基づいて前記第１のメモリ部に記憶された血流速度データに対しパーシスタンス演算を行う第２のパーシスタンス演算部と、前記第１のパーシスタンス演算部から出力される演算結果の絶対値および前記第２パーシスタンス演算部から出力される演算結果の絶対値を比較し、大きい方の演算結果を前記残像処理された血流速度データとして出力する最大値選択部とを含む。
【００２８】
ある好ましい実施形態において、前記折り返し判定部は、前記第１のメモリ部に記憶された血流速度データおよび前記第２のメモリ部に記憶された血流速度データと複数の閾値とを比較することによって、折り返しが発生しているかどうか、および、前記最新のフレームの血流速度データが折り返し領域にあるかどうかを判定する。
【００２９】
ある好ましい実施形態において、前記残像処理部は、前記血流速度の値に応じて対応付けられた異なる２つ以上の値の第１のパーシスタンス係数を含む第１の参照テーブルを記憶している第３のメモリ部と、前記血流速度の値に応じて対応付けられた異なる２つ以上の値の第２のパーシスタンス係数を含む第２の参照テーブルを記憶している第４のメモリ部とをさらに含む。
【００３０】
ある好ましい実施形態において、前記第１の参照テーブルと前記第２の参照テーブルにおいて、同じ血流速度の値に応じて対応付けられた第１のパーシスタンス係数と第２のパーシスタンス係数とは互いに異なる値である。
【００３１】
ある好ましい実施形態において、前記第１の参照テーブルは、所定の値以上の血流速度に対して一定の値のパーシスタンス係数が対応付けられている。
【００３２】
ある好ましい実施形態において、前記最新より以前のフレームの残像処理された血流速度データは、前記最新より１つ前のフレームの残像処理された血流速度データである。
【発明の効果】
【００３３】
本発明によれば、最新のフレームの血流速度データおよび最新より前のフレームの血流速度データに基づき、折り返しの判定を行い、折り返し判定の結果および最新のフレームの血流速度データに基づき、パーシスタンス係数を動的に変更する。このため、血流変化が明確に認識でき、かつ、低速な血流速度でも黒抜けが発生しない滑らかな血流動画を表示することが可能な超音波診断装置が実現する。
【００３４】
また、本発明によれば、最新のフレームの血流速度データに基づき決定されたパーシスタンス係数を用いて残像処理された血流速度、および、最新より前のフレームの血流速度データに基づき決定されたパーシスタンス係数を用いて残像処理された血流速度を求め、絶対値の大きいほうを選択して血流画像表示に用いる。このため、血流パワーが不安定となる末梢血管の血流表示を点滅させずに、かつ、平滑化に起因する末梢血管の血流消滅を起こさない血流動画を表示することが可能となる。
【図面の簡単な説明】
【００３５】
【図１】本発明による超音波診断装置の第１の実施形態を示すブロック図である。
【図２】（ａ）および（ｂ）は、第１の実施形態において、パーシスタンス係数を用いて血流速度データのパーシスタンス演算を行う場合に考慮する折り返し演算を説明する模式図である。
【図３】（ａ）は、第１の実施形態における折り返し判定を説明するための模式図であり、（ｂ）は、参照テーブルのデータが満たす関係を示すグラフである。
【図４】本発明による超音波診断装置の第２の実施形態を示すブロック図である。
【図５】（ａ）および（ｂ）は、第２の実施形態において、パーシスタンス係数を用いて血流速度データのパーシスタンス演算を行う場合に考慮する折り返し演算を説明する模式図である。
【図６】（ａ）は、第２の実施形態における折り返し判定を説明するための模式図であり、（ｂ）および（ｃ）は、第１および第２の参照テーブルのデータが満たす関係を示すグラフである。
【図７】従来の超音波診断装置を示すブロック図である。
【発明を実施するための形態】
【００３６】
（第１の実施形態）
以下、図面を参照しながら本発明による超音波診断装置の第１の実施形態を説明する。図１は、本発明による超音波診断装置の第１の実施形態を示すブロック図である。図１に示す超音波診断装置１１は、探触子１０１と、超音波送受信部１０２と、ＣＦＭ信号処理部１０３と、残像処理部１１５と、断層画像信号処理部１１１と、ＣＦＭ  ＤＳＣ部１１０と、断層画像ＤＳＣ部１１２と、画像合成部１１３と、表示部１１４とを備える。これらの構成のうち、探触子１０１および表示部１１４には汎用の探触子および表示装置を用いることが可能であり、探触子１０１および表示部１１４を超音波診断装置１１が備えていなくてもよい。
【００３７】
超音波送受信部１０２は、探触子１０１を駆動する駆動信号を生成し、探触子１０１へ出力することによって、探触子１０１から被検体に向けて超音波を送信する。また、送信された超音波が被検体において反射することによって得られた反射エコーを探触子１０１によって受信し、受信信号を生成する。より具体的には、探触子１０１は複数の圧電素子を含み、各圧電素子から送信される超音波が超音波ビームを構成し、被検体を複数の超音波ビームによって走査するように、超音波送受信部１０２が各圧電素子の遅延制御を行いながら探触子１０１を駆動する。反射エコーは各圧電素子によって受信され、超音波送受信部１０２が、各圧電素子の遅延制御を行うことにより、送信された超音波ビームに対応する受信信号が生成する。被検体を超音波ビームが１回走査することによって１フレーム分のデータが得られる。１秒間に数回〜数十回繰り返して超音波の送受信を行うことにより、毎秒数フレームから数十フレームの受信信号が逐次生成する。
【００３８】
本実施形態の超音波診断装置１１は、Ｂモード断層画像とカラーフローマッピング画像と生成し、これらを合成して表示部１１４に表示する。このため、超音波送受信部１０２の上述した超音波の送受信は、Ｂモード断層画像の生成およびカラーフローマッピング画像の生成のそれぞれについて行われる。Ｂモード断層画像の１秒当たりのフレーム数およびカラーフローマッピング画像の１秒当たりのフレーム数は同じであってもよいし、異なっていてもよい。フレーム数が同じである場合には、Ｂモード断層画像生成のための超音波の送受信およびカラーフローマッピング画像生成のための超音波の送受信を交互に繰り返し行ってもよい。
【００３９】
Ｂモード断層画像を生成する場合、超音波送受信部１０２は、Ｂモード断層画像の生成に適した送受信を行い、得られた受信信号を断層画像信号処理部１１１へ出力する。カラーフローマッピング断層画像を生成する場合は、カラーフローマッピング断層画像の生成に適した送受信を行い、得られた受信信号をＣＦＭ信号処理部１０３へ出力する。一般的に、カラーフローマッピング断層画像を生成する場合、安定したカラーフローマッピング断層画像を得るために、超音波送受信部１０２は、同じ音響線上において超音波の送受信を複数回行う。
【００４０】
ＣＦＭ信号処理部１０３は、受信信号に対して直交検波処理、ＭＴＩフィルタ処理および自己相関処理を行い、血流速度および血流パワーを算出したのち、システムあるいは音響ノイズを排除するノイズカット処理を行う。ＣＦＭフレームデータは少なくとも血流速度データを含む。この他、血流パワーデータや血流速度の分散データを含んでいてもよい。ＣＦＭ信号処理部１０３はこの処理を、各フレームを構成する受信信号ごとに逐次繰り返して行う。ＣＦＭ信号処理部１０３で生成したＣＦＭフレームデータはフレームごとに残像処理部１１５へ出力される。
【００４１】
残像処理部１１５は、パーシスタンス係数を用いて、フレームごとにＣＦＭフレームデータに対して残像処理を行う。本実施形態の超音波診断装置１１は、パーシスタンス係数を最新のフレームの血流速度に応じて決定する。つまり、パーシスタンス係数は一定ではなく、最新のフレームの血流速度に基づく動的な値である。これにより、血流速度に応じてパーシスタンス係数を変化させ、残像効果を調整することができる。ただし、血流を動画で表示するためには、パルスドプラ法を用いて超音波の送受信を行う必要があり、このため、計測できる血流速度はパルス繰り返し周波数（ＰＲＦ）による制限を受ける。その結果、血流速度に折り返しが生じ、血流速度を正確に評価することが困難となる。
【００４２】
本実施形態の超音波診断装置１１は折り返しが生じているかどうかを判定するため、最新のフレームの血流速度データと１つ前のフレームの血流速度データを用いる。このために、残像処理部１１５は、フレームメモリ部（第１のメモリ部）１０４と、折り返し判定部１０５と、パーシスタンス係数決定部１０６と、パーシスタンス係数参照メモリ部（第３のメモリ部）１０７とパーシスタンス演算部１０８とパーシスタンスメモリ部（第２のメモリ部）１０９とを含む。
【００４３】
フレームメモリ部１０４は、最新のフレームの（現在走査中の）ＣＦＭフレームデータを記憶する。パーシスタンスメモリ部１０９は、最新より１つ前のフレームのパーシスタンス演算部１０８の出力結果であるＣＦＭフレームデータを記憶する。パーシスタンスメモリ部１０９のＣＦＭフレームデータは残像処理が施されている。以下、フレームメモリ部１０４およびパーシスタンスメモリ部１０９にそれぞれ記憶されたＣＦＭフレームデータのうち、血流速度データをＶcurrentおよび血流速度データＶout-1と呼ぶ。
【００４４】
折り返し判定部１０５は、フレームメモリ部１０４からＣＦＭフレームデータのうちの血流速度データＶcurrentとパーシスタンスメモリ部１０９からＣＦＭフレームデータのうちの血流速度データＶout-1を読み出し、折り返し判定を行う。より具体的には、血流速度データＶcurrentおよび血流速度データＶout-1と複数の閾値とを比較することによって、折り返しが発生しているかどうか、および、血流速度データＶcurrentが折り返し領域にあるかどうかを判定し、結果をパーシスタンス係数決定部１０６およびパーシスタンス演算部１０８に出力する。
【００４５】
パーシスタンス係数決定部１０６は、折り返し判定部１０５からの２つの判定結果とフレームメモリ部１０４から読み出した血流速度データＶcurrentに基づいて、パーシスタンス係数参照メモリ部１０７への参照インデックスを作成する。また、パーシスタンス係数参照メモリ部１０７にアクセスし、参照インデックスに対応付けられたパーシスタンス係数を読み出し、パーシスタンス演算部１０８に設定する。パーシスタンス係数参照メモリ部１０７には、あらかじめ血流速度の値に対応付けられたパーシスタンス係数の参照テーブルが記憶されている。この参照テーブルは、血流速度の値に応じて対応付けられた異なる２つ以上の値のパーシスタンス係数を含む。
【００４６】
パーシスタンス演算部１０８は、パーシスタンス係数決定部１０６から設定されたパーシスタンス係数と折り返し判定部１０５からの折り返し発生の判定結果とに基づいて、以下に示す式（１）で血流速度データに対してパーシスタンス演算を行う。パーシスタンス演算によって求められる残像処理された血流速度データをＶoutとし、パーシスタンス係数をＣpersistence（０＜Ｃpersistence＜１）とすると、残像処理された血流速度データは以下の式（１）で求められる。
Ｖout＝（１−Ｃpersistence）×Ｖcurrent＋Ｃpersistence×Ｖout-1 ・・（１)
【００４７】
ＣＦＭフレームデータが、血流速度データ以外のデータを含む場合には、同様に最新のフレームのデータおよび最新より１つ前のフレームのデータと求めたパーシスタンス係数Ｃpersistenceを用いてパーシスタンス演算を行い、残像処理されたデータを得る。
【００４８】
折り返し判定部１０５からの折り返し発生の判定結果が真であった場合は、式（１）の演算式を符号なし演算で扱い、偽であった場合は、符号付演算で扱う。
【００４９】
上述したように測定にパルス波を用いるため、ドップラ偏移によって直接測定できる血流速度はパルス波の繰り返し周波数（ＰＲＦ）の制限を受ける。具体的には、±ＰＲＦ／２を超える周波数変化に対応する血流速度は反対向きの血流として観測される折り返しが発生する。
【００５０】
図２（ａ）および（ｂ）は、残像処理された血流速度データＶout、最新のフレームの血流速度データＶcurrentおよび最新より１つ前のフレームのパーシスタンス演算部１０８の出力結果である血流速度データＶout-1の大きさの関係を示している。図２（ａ）および（ｂ）において、横軸の第１象限部分は速度Ｖがゼロを意味し、横軸の第２象限部分は、＋Ｖまたは−Ｖを意味する。速度Ｖが正である場合、第１または第２象限に位置し、速度Ｖが負である場合には、第３または第４象限に位置する。
【００５１】
例えば、図２（ａ）に示すように、Ｖcurrentが第２象限にあり、Ｖout-1が第３象限にあり、折り返しが発生していると判定された場合、Ｖout-1は実際には、＋ＰＲＦ／２に対応する血流速度よりも大きな値であったことなるから、ゼロを通らない、つまり、符号変化をともなわない演算である。したがって、ＶcurrentおよびＶout-1の符号（プラスまたはマイナス）をとって、式（１）にこれらの値を代入し演算を行う。
【００５２】
一方、例えば、図２（ｂ）に示すように、Ｖcurrentが第１象限にあり、Ｖout-1が第４象限にあり、折り返しが発生していないと判定された場合、式（１）の演算は、ゼロを通り、符号変化が生じる演算となる。したがって、符号をつけたＶcurrentおよびＶout-1を式（１）に代入し演算を行う。この演算は、１フレーム分の血流速度データの各画素または測定点ごとに行われる。また、演算結果であるＶoutは、折り返しが発生している場合、符号なしの値となる。この場合、血流速度データＶoutの最上位ビットを符号として扱うことで符号あり値としてＣＦＭ  ＤＳＣ部１１０およびパーシスタンスメモリ部１０９に出力される。
【００５３】
ＣＦＭ  ＤＳＣ部１１０は、パーシスタンス演算部１０８から出力される血流速度データの座標を変換し、画像合成部１１３へ出力する。
【００５４】
断層画像信号処理部４０９は受信信号に対してダイナミックフィルタ処理を施すことにより、不要なノイズをカットしたのち、包絡線検波処理およびダイナミックレンジ圧縮処理を施して断層画像ＤＳＣ部４１０に断層画像フレームデータを出力する。断層画像ＤＳＣ部４１０は、断層画像信号処理部４０９からの断層画像フレームデータの座標を変換し、画像合成部４１１へ出力する。
【００５５】
画像合成部４１１は、ＣＦＭ  ＤＳＣ部４１０および断層画像ＤＳＣ部４１０から出力される各フレームデータを画素ごとに、あるいは、対応する測定点のデータごとに合成し、合成画像フレームデータを生成する。具体的には、血流速度がゼロである場合は、断層画像フレームデータを表示し、そうでない場合は、ＣＦＭフレームデータを表示するように２つのデータを画素ごとに、あるいは、対応する測定点のデータごとに合成する。また、血流速度や血流の方向に応じてデータを色情報に変換し、表示部４１２に出力する。表示部４１２は、画像合成部４１１から受け取ったデータを表示する。
【００５６】
次に、パーシスタンス係数の決定についてより詳細に説明する。パーシスタンス係数を決定するために、まず、血流速度に折り返しが発生しているかどうかを折り返し判定部１０５において判定する。
【００５７】
折り返し判定部１０５は、フレームメモリ部１０４から最新のＣＦＭフレームデータに含まれる血流速度データＶcurrentとパーシスタンスメモリ部１０９から１フレーム前のパーシスタンス演算部１０８の出力結果であるＣＦＭフレームデータに含まれる血流速度Ｖout-1とを読み出し、ＶcurrentおよびＶout-1の値から以下の２つの判定を行う。
１．折り返しが発生しているかどうか。
２．Ｖcurrentが折り返し領域にあるかどうか。
【００５８】
この２つの状態の判定には、あらかじめ定められた閾値とＶcurrentおよびＶout-1とを比較することにより行う。具体的には、閾値Ｖthおよび血流ゼロ速度ＶzeroとＶcurrentおよびＶout-1とを比較する。
【００５９】
【表１】
2011036891.tif 000015
【００６０】
図３（ａ）は、閾値Ｖth、血流ゼロ速度Ｖzero、ＶcurrentおよびＶout-1の大きさの関係を示している。図３（ａ）において、横軸の第１象限部分は血流ゼロ速度Ｖzeroを意味し、横軸の第２象限部分は、Ｖmaxまたは−Ｖmaxを意味する。速度Ｖが正である場合、第１または第２象限に位置し、速度Ｖが負である場合には、第３または第４象限に位置する。
【００６１】
ここで、Ｖthおよび−Ｖthには、例えば、隣接するフレームの時間間隔において、想定される血流速度の変化の最大値を設定する。
【００６２】
表１は、折り返し判定部１０５において判定される条件と判定結果とを示している。
【００６３】
条件（０）に示すように、Ｖout-1が正である場合、想定される血流速度の変化の最大値はＶthまたは−Ｖthであるから、Ｖcurrentが−Ｖthよりも小さくなることはない。したがって、Ｖcurrent＜−Ｖthを満たしていれば、Ｖcurrentは、実際には、＋ＰＲＦ／２に対応する最大血流速度Ｖmaxよりも大きな値となっており、折り返しが発生しているとともに、Ｖcurrentが折り返し領域にあると判定される。条件（１）は、条件（０）の符号が逆転した場合である。
【００６４】
条件（２）に示すように、Ｖout-1が−Ｖthより小さい場合、Ｖcurrentが正の値となることは、想定される血流速度の変化の最大値を超える変化であるから、折り返しが発生している。また、ＶcurrentがＶzeroを挟む±Ｖthの範囲にあるため、Ｖcurrentは折り返し領域ではない。条件（３）は、条件（２）の符号が逆転した場合である。
【００６５】
条件（０）から（３）のいずれも満たさない場合には、折り返しが発生しておらず、また、Ｖcurrentが折り返し領域にないと判定する。
【００６６】
パーシスタンス係数決定部１０６は、折り返し判定部１０５から出力される２つの判定結果とフレームメモリ部１０４から読み出した血流速度データＶcurrentの絶対値に基づいて、パーシスタンス係数参照メモリ部１０７への参照インデックスを作成する。表２に作成される参照インデックスを示す。
【００６７】
【表２】
2011036891.tif 000016
【００６８】
折り返しが発生しており、かつ、Ｖcurrentが折り返し領域にある場合、血流速度Ｖcurrentは、実際には、Ｖmaxをまたは−Ｖmaxを超える大きな値と考えられる。このため、参照インデックスはＶmaxとなる。その他の場合には、Ｖcurrentの絶対値Abs(Ｖcurrent)となる。
【００６９】
パーシスタンス係数参照メモリ部１０７には、参照インデックスに対応付けられたパーシスタンス係数によって構成される参照テーブルが記憶されている。パーシスタンス係数決定部１０６は、パーシスタンス係数参照メモリ部１０７にアクセスし、作成された参照インデックスに対応付けられたパーシスタンス係数を読み出し、パーシスタンス演算部１０８へ出力する。
【００７０】
図３（ｂ）は、参照インデックスとパーシスタンス係数との対応関係の一例を示すグラフである。図３（ｂ）において、横軸は参照インデックスを示し、縦軸は、パーシスタンス係数を示している。表２に示したように、参照インデックスは、ＶmaxまたはＶcurrentの絶対値Abs(Ｖcurrent)である。Ｖcurrentの絶対値が閾値Ｖth以下である場合には、Ｖcurrentの増大にともなって単調に減少するパーシスタンス係数Ｃpersistenceが対応付けられている。つまり、Ｖcurrentの絶対値が閾値Ｖth以下である場合には、最新のフレームの血流速度Ｖcurrentに応じて異なるパーシスタンス係数Ｃpersistenceが対応付けられる。これにより、最新のフレームの血流速度Ｖcurrentが小さい場合には、パーシスタンス係数Ｃpersistenceが大きくなる。つまり、１つ前のフレームの血流速度Ｖout-1の重みが大きくなる。その結果、最新のフレームの血流速度Ｖcurrentが小さい場合には、１つ前のフレームの血流速度Ｖout-1を大きく反映させた血流速度Ｖoutが決定され表示部１１４に表示される。このため、カラーフローマッピング画像の変化が滑らかとなり、黒抜けが生じにくくなる。
【００７１】
また、最新のフレームの血流速度Ｖcurrentが大きい場合には、パーシスタンス係数Ｃpersistenceが小さくなる。つまり、１つ前のフレームの血流速度Ｖout-1の重みが小さくなる。その結果、最新のフレームの血流速度Ｖcurrentが大きい場合には、１つ前のフレームの血流速度Ｖout-1の影響が小さくなり、急激な血流速度の増大をリアルタイムで反映させたカラーフローマッピング画像を実現できる。
【００７２】
また、Ｖcurrentの増大にともなってパーシスタンス係数Ｃpersistenceが単調に減少するため、時間の経過とともに血流速度が増大する場合には、パーシスタンス係数Ｃpersistenceが減少し、残像効果が小さくなり、カラーフローマッピング画像の変化が急激になる。時間の経過とともに血流速度が減少する場合には、パーシスタンス係数Ｃpersistenceが増大し、残像効果が大きくなり、カラーフローマッピング画像の変化が緩やかになる
。
【００７３】
なお、表１および表２から分かるように、Ｖcurrent＜−Ｖthでも、Ｖout-1＞０であれば、参照インデックスはＶmaxとなる（条件（０））一方、Ｖout-1＜０であれば、参照インデックスはＶcurrentの絶対値Abs(Ｖcurrent)となる（条件（４））。このため、Ｖcurrent＜−Ｖthを満たしていても、Ｖout-1が正か負であるかによって、参照インデックスが異なり、パーシスタンス係数Ｃpersistenceも異なる。その結果、Ｖcurrent＜−Ｖthを満たす隣接する領域であっても、Ｖout-1が正か負であるかによってカラーマッピングフロー画像として表示される色が異なり、画像に不連続な色調部分が生じてしまう。
【００７４】
このような不自然な表示を抑制するため、Ｖcurrentの絶対値が閾値Ｖth以上である場合には、参照インデックスに同じ値のパーシスタンス係数Ｃpersistenceが対応付けられていることが好ましい。これにより、折り返しが発生する血流領域、あるいはその境界近傍において、自然な表示を行うことができる。
【００７５】
このように本実施形態の超音波診断装置によれば、ＣＦＭフレームデータに対して血流速度および折り返し状態に基づいてパーシスタンス係数を動的に決定したのちにパーシスタンス演算を実施することにより、頚動脈のような血流変化が激しい診断部位においても血流変化が明確に認識でき、かつ、低速な血流速度でも黒抜けが発生しない滑らかな血流動画を表示することが可能となる。
【００７６】
なお、上記実施形態では、ＣＦＭフレームデータの血流速度に基づいてパーシスタンス係数を動的に決定し、血流速度にパーシスタンス演算を行っているが、上述したようにＣＦＭフレームデータの他のデータ、例えば、血流パワーデータに対してパーシスタンス演算を行ってもよいし、Ｂモード断層画像データに対してパーシスタンス演算を行ってもよい。
【００７７】
また、上記実施形態では、最新のフレームおよび１つ前のフレームの血流速度データを用いてパーシスタンス処理を行っていたが、２つ前、あるいは３つ以上前のフレームの血流速度データも用いてパーシスタンス処理を行ってもよい。また、式（１）に限られず、他の演算式を用いてパーシスタンス処理を行ってもよい。
【００７８】
（第２の実施形態）
以下、図面を参照しながら本発明による超音波診断装置の第２の実施形態を説明する。図４は、本発明による超音波診断装置の一実施形態を示すブロック図である。図４に示す超音波診断装置１２は、探触子１０１と、超音波送受信部１０２と、ＣＦＭ信号処理部１０３と、残像処理部１１５’と、断層画像信号処理部１１１と、ＣＦＭ  ＤＳＣ部１１０と、断層画像ＤＳＣ部１１２と、画像合成部１１３と、表示部１１４とを備える。これらの構成のうち、探触子１０１および表示部１１４には汎用の探触子および表示装置を用いることが可能であり、探触子１０１および表示部１１４を超音波診断装置１２が備えていなくてもよい。
【００７９】
第１の実施形態で説明したように、超音波送受信部１０２は、探触子１０１を駆動する駆動信号を生成し、探触子１０１へ出力することによって、探触子１０１から被検体に向けて超音波を送信する。また、送信された超音波が被検体において反射することによって得られた反射エコーを探触子１０１によって受信し、受信信号を生成する。より具体的には、探触子１０１は複数の圧電素子を含み、各圧電素子から送信される超音波が超音波ビームを構成し、被検体を複数の超音波ビームによって走査するように、超音波送受信部１０２が各圧電素子の遅延制御を行いながら探触子１０１を駆動する。反射エコーは各圧電素子によって受信され、超音波送受信部１０２が、各圧電素子の遅延制御を行うことにより、送信された超音波ビームに対応する受信信号が生成する。被検体を超音波ビームが１回走査することによって１フレーム分のデータが得られる。１秒間に数回〜数十回繰り返して超音波の送受信を行うことにより、毎秒数フレームから数十フレームの受信信号が逐次生成する。
【００８０】
本実施形態の超音波診断装置１２は、Ｂモード断層画像とカラーフローマッピング画像と生成し、これらを合成して表示部１１４に表示する。このため、超音波送受信部１０２の上述した超音波の送受信は、Ｂモード断層画像の生成およびカラーフローマッピング画像の生成のそれぞれについて行われる。Ｂモード断層画像の１秒当たりのフレーム数およびカラーフローマッピング画像の１秒当たりのフレーム数は同じであってもよいし、異なっていてもよい。フレーム数が同じである場合には、Ｂモード断層画像生成のための超音波の送受信およびカラーフローマッピング画像生成のための超音波の送受信を交互に繰り返し行ってもよい。
【００８１】
Ｂモード断層画像を生成する場合、超音波送受信部１０２は、Ｂモード断層画像の生成に適した送受信を行い、得られた受信信号を断層画像信号処理部１１１へ出力する。カラーフローマッピング断層画像を生成する場合は、カラーフローマッピング断層画像の生成に適した送受信を行い、得られた受信信号をＣＦＭ信号処理部１０３へ出力する。一般的に、カラーフローマッピング断層画像を生成する場合、安定したカラーフローマッピング断層画像を得るために、超音波送受信部１０２は、同じ音響線上において超音波の送受信を複数回行う。
【００８２】
ＣＦＭ信号処理部１０３は、受信信号に対して直交検波処理、ＭＴＩフィルタ処理および自己相関処理を行い、血流速度および血流パワーを算出したのち、システムあるいは音響ノイズを排除するノイズカット処理を行う。ＣＦＭフレームデータは少なくとも血流速度データを含む。この他、血流パワーデータや血流速度の分散データを含んでいてもよい。ＣＦＭ信号処理部１０３はこの処理を、各フレームを構成する受信信号ごとに逐次繰り返して行う。ＣＦＭ信号処理部１０３で生成したＣＦＭフレームデータはフレームごとに残像処理部１１５’へ出力される。
【００８３】
残像処理部１１５’は、パーシスタンス係数を用いてフレームごとにＣＦＭフレームデータに対して残像処理を行う。本実施形態の超音波診断装置１２は、パーシスタンス係数を血流速度に応じてパーシスタンス係数を決定する。つまり、パーシスタンス係数は一定ではなく、血流速度に応じた動的な値である。これにより、血流速度に応じてパーシスタンス係数を変化させ、残像効果を調整することができる。ただし、血流を動画で表示するためには、パルスドプラ法を用いて超音波の送受信を行う必要があり、このため、計測できる血流速度はパルス繰り返し周波数（ＰＲＦ）による制限を受ける。その結果、血流速度に折り返しが生じ、血流速度を正確に評価することが困難となる。
【００８４】
本実施形態の超音波診断装置１２は折り返しが生じているかどうかを判定するため、最新のフレームの血流速度データと１つまえのフレームの血流速度データを用いる。また、残像処理部１１５’は２つパーシスタンス演算部を備え、残像効果をあまり与えずに速やかに血流速度を変化させる第１のパーシスタンス演算と残像効果を強力に与え血流速度変化をなるべく維持する第２のパーシスタンス演算を同時に行う。このようにして生成した残像効果の異なる２つの血流速度データのうち絶対値の大きい方を用いて血流画像を構成する。これにより、血流パワーが小さい末梢血管でも、血流表示を点滅させずに、かつ、平滑化に起因する末梢血管の血流消滅を起こさない動画を表示することが可能となる。
【００８５】
このために、残像処理部１１５’は、フレームメモリ部（第１のメモリ部）１０４と、折り返し判定部１０５と、第１のパーシスタンス係数決定部１０６Ａと、第１のパーシスタンス係数参照メモリ部（第３のメモリ部）１０７Ａと第１のパーシスタンス演算部１０８Ａと、第２のパーシスタンス係数決定部１０６Ｂと、第２のパーシスタンス係数参照メモリ部（第４のメモリ部）１０７Ｂと第２のパーシスタンス演算部１０８Ｂと、最大値選択部１１６と、パーシスタンスメモリ部（第２のメモリ部）１０９とを含む。
【００８６】
フレームメモリ部１０４は、最新のフレームの（現在走査中の）ＣＦＭフレームデータを記憶する。パーシスタンスメモリ部１０９は、最新より１つ前の最大値選択部の出力結果であるＣＦＭフレームデータを記憶する。パーシスタンスメモリ部１０９のＣＦＭフレームデータは残像処理が施されている。第１の実施形態と同様、フレームメモリ部１０４およびパーシスタンスメモリ部１０９にそれぞれ記憶されたＣＦＭフレームデータのうち、血流速度データをＶcurrentおよび血流速度データＶout-1と呼ぶ。
【００８７】
折り返し判定部１０５は、フレームメモリ部１０４からＣＦＭフレームデータのうちの血流速度データＶcurrentとパーシスタンスメモリ部１０９からＣＦＭフレームデータのうちの血流速度データＶout-1を読み出し、折り返し判定を行う。より具体的には、血流速度データＶcurrentおよび血流速度データＶout-1と複数の閾値とを比較することによって、折り返しが発生しているかどうか、および、血流速度データＶcurrentが折り返し領域にあるかどうかを判定し、結果を第１のパーシスタンス係数決定部１０６Ａ、第２のパーシスタンス係数決定部１０６Ｂ、第１のパーシスタンス演算部１０８Ａおよび第２のパーシスタンス演算部１０８Ｂに出力する。
【００８８】
第１のパーシスタンス係数決定部１０６Ａは、折り返し判定部１０５からの２つの判定結果とフレームメモリ部１０４から読み出した血流速度データＶcurrentに基づいて、第１のパーシスタンス係数参照メモリ部１０７Ａへの参照インデックスを作成する。また、第１のパーシスタンス係数参照メモリ部１０７Ａにアクセスし、参照インデックスに対応付けられた第１のパーシスタンス係数を読み出し、第１のパーシスタンス演算部１０８Ａに設定する。第１のパーシスタンス係数参照メモリ部１０７Ａには、あらかじめ血流速度の値に対応付けられた第１のパーシスタンス係数を含む第１の参照テーブルが記憶されている。この第１の参照テーブルは、血流速度の値に応じて対応付けられた異なる２つ以上の値のパーシスタンス係数を含む。
【００８９】
これに対し、第２のパーシスタンス係数決定部１０６Ｂは、折り返し判定部１０５からの２つの判定結果とパーシスタンスメモリ部１０９から読み出した血流速度データＶout-1に基づいて、第２のパーシスタンス係数参照メモリ部１０７Ｂへの参照インデックスを作成する。また、第２のパーシスタンス係数参照メモリ部１０７Ｂにアクセスし、参照インデックスに対応付けられた第２のパーシスタンス係数を読み出し、第２のパーシスタンス演算部１０８Ｂに設定する。第２のパーシスタンス係数参照メモリ部１０７Ｂには、あらかじめ血流速度の値に対応付けられた第２のパーシスタンス係数を含む第２の参照テーブルが記憶されている。第２の参照テーブルも、血流速度の値に応じて対応付けられた異なる２つ以上の値のパーシスタンス係数を含むが、以下において詳細に説明するように、同じ血流速度値に対して対応付けられる第１のパーシスタンス係数値と第２のパーシスタンス係数値とは値が異なっている。
【００９０】
第１のパーシスタンス演算部１０８Ａは、第１のパーシスタンス係数決定部１０６Ａから設定されたパーシスタンス係数と折り返し判定部１０５からの折り返し発生の判定結果とに基づいて、以下に示す式（１）で血流速度データに対してパーシスタンス演算を行う。
【００９１】
パーシスタンス演算によって求められる残像処理された血流速度データをＶoutとし、パーシスタンス係数をＣpersistence（０＜Ｃpersistence＜１）とすると、残像処理された血流速度データは以下の式（１）で求められる。
Ｖout＝（１−Ｃpersistence）×Ｖcurrent＋Ｃpersistence×Ｖout-1 ・・（１)
【００９２】
同様に第２のパーシスタンス演算部１０８Ｂも、第２のパーシスタンス係数決定部１０６Ｂから設定されたパーシスタンス係数と折り返し判定部１０５からの折り返し発生の判定結果とに基づいて、式（１）で血流速度データに対してパーシスタンス演算を行う。
【００９３】
第１のパーシスタンス演算部１０８Ａおよび第２のパーシスタンス演算部１０８Ｂにおける演算は、決定したパーシスタンス係数が互い異なる点を除いて同じである。ＣＦＭフレームデータが、血流速度データ以外のデータを含む場合には、同様に最新のフレームのデータおよび最新より１つ前のフレームのデータと求めたパーシスタンス係数Ｃpersistenceを用いて第１のパーシスタンス演算部１０８Ａおよび第２のパーシスタンス演算部１０８Ｂでパーシスタンス演算を行い、残像処理されたデータをそれぞれ得る。
【００９４】
折り返し判定部１０５からの折り返し発生の判定結果が真であった場合は、式（１）の演算式を符号なし演算で扱い、偽であった場合は、符号付演算で扱う。
【００９５】
上述したように測定にパルス波を用いるため、ドップラ偏移によって直接測定できる血流速度はパルス波の繰り返し周波数（ＰＲＦ）の制限を受ける。具体的には、±ＰＲＦ／２を超える周波数変化に対応する血流速度は反対向きの血流として観測される折り返しが発生する。
【００９６】
図５（ａ）および（ｂ）は、残像処理された血流速度データＶout、最新のフレームの血流速度データＶcurrentおよび最新より１つ前のフレームのパーシスタンス演算部１０８の出力結果である血流速度データＶout-1の大きさの関係を示している。図５（ａ）および（ｂ）において、横軸の第１象限部分は速度Ｖがゼロを意味し、横軸の第２象限部分は、＋Ｖまたは−Ｖを意味する。速度Ｖが正である場合、第１または第２象限に位置し、速度Ｖが負である場合には、第３または第４象限に位置する。
【００９７】
例えば、図５（ａ）に示すように、Ｖcurrentが第２象限にあり、Ｖout-1が第３象限にあり、折り返しが発生していると判定された場合、Ｖout-1は実際には、＋ＰＲＦ／２に対応する血流速度よりも大きな値であったことなるから、ゼロをとらない、つまり、符号変化を伴わない演算である。したがって、ＶcurrentおよびＶout-1の符号（プラスまたはマイナス）をとって、式（１）にこれらの値を代入し演算を行う。
【００９８】
一方、例えば、図５（ｂ）に示すように、Ｖcurrentが第１象限にあり、Ｖout-1が第４象限にあり、折り返しが発生していないと判定された場合、式（１）の演算は、ゼロを通り、符号変化が生じる演算となる。したがって、符号をつけたＶcurrentおよびＶout-1を式（１）に代入し演算を行う。この演算は、１フレーム分の血流速度データの各画素または測定点ごとに行われる。また、演算結果であるＶoutは、折り返しが発生している場合、符号なしの値となる。この場合、血流速度データＶoutの最上位ビットを符号として扱うことで符号あり値として最大値選択部１１６にそれぞれ出力される。
【００９９】
最大値選択部１１６は、第１のパーシスタンス演算部１０８Ａおよび第２のパーシスタンス演算部１０８Ｂからそれぞれ演算結果、つまり、残像処理された血流速度データを受け取り、画素ごとに、あるいは、対応する測定点のデータごとに血流速度の絶対値を比較し、大きい方の血流速度を選択して、最新のフレームの残像処理された血流速度データを構成し、これをＣＦＭ  ＤＳＣ部１１０およびパーシスタンスメモリ部１０９に出力する。ＣＦＭ  ＤＳＣ部１１０は、選択した血流速度データの座標を変換し、画像合成部１１３へ出力する。
【０１００】
断層画像信号処理部４０９は受信信号に対してダイナミックフィルタ処理を施すことにより、不要なノイズをカットしたのち、包絡線検波処理およびダイナミックレンジ圧縮処理を施して断層画像ＤＳＣ部４１０に断層画像フレームデータを出力する。断層画像ＤＳＣ部４１０は、断層画像信号処理部４０９からの断層画像フレームデータの座標を変換し、画像合成部４１１へ出力する。
【０１０１】
画像合成部４１１は、ＣＦＭ  ＤＳＣ部４１０および断層画像ＤＳＣ部４１０から出力される各フレームデータを画素ごとに、あるいは、対応する測定点のデータごとに合成し、合成画像フレームデータを生成する。具体的には、血流速度がゼロである場合は、断層画像フレームデータを表示し、そうでない場合は、ＣＦＭフレームデータを表示するように２つのデータを画素ごとに、あるいは、対応する測定点のデータごとに合成する。また、血流速度や血流の方向に応じてデータを色情報に変換し、表示部４１２に出力する。表示部４１２は、画像合成部４１１から受け取ったデータを表示する。
【０１０２】
次に、第１および第２のパーシスタンス係数の決定についてより詳細に説明する。第１および第２のパーシスタンス係数を決定するために、まず、血流速度に折り返しが発生しているかどうかを折り返し判定部１０５において判定する。
【０１０３】
折り返し判定部１０５は、フレームメモリ部１０４から最新のＣＦＭフレームデータに含まれる血流速度データＶcurrentとパーシスタンスメモリ部１０９から１フレーム前のパーシスタンス演算部１０８の出力結果であるＣＦＭフレームデータに含まれる血流速度Ｖout-1とを読み出し、ＶcurrentおよびＶout-1の値から以下の２つの判定を行う。
１．折り返しが発生しているかどうか。
２．Ｖcurrentが折り返し領域にあるかどうか。
【０１０４】
この２つの状態の判定には、あらかじめ定められた閾値とＶcurrentおよびＶout-1とを比較することにより行う。具体的には、閾値Ｖthおよび血流ゼロ速度ＶzeroとＶcurrentおよびＶout-1とを比較する。
【０１０５】
【表３】
2011036891.tif 000017
【０１０６】
図６（ａ）は、閾値Ｖth、血流ゼロ速度Ｖzero、ＶcurrentおよびＶout-1の大きさの関係を示している。図６（ａ）において、横軸の第１象限部分は血流ゼロ速度Ｖzeroを意味し、横軸の第２象限部分は、Ｖmaxまたは−Ｖmaxを意味する。速度Ｖが正である場合、第１または第２象限に位置し、速度Ｖが負である場合には、第３または第４象限に位置する。
【０１０７】
ここで、Ｖthおよび−Ｖthには、例えば、隣接するフレームの時間間隔において、想定される血流速度の変化の最大値を設定する。
【０１０８】
表３は、折り返し判定部１０５において判定される条件と判定結果とを示している。
【０１０９】
条件（０）に示すように、Ｖout-1が正である場合、想定される血流速度の変化の最大値はＶthまたは−Ｖthであるから、Ｖcurrentが−Ｖthよりも小さくなることはない。したがって、Ｖcurrent＜−Ｖthを満たしていれば、Ｖcurrentは、実際には、＋ＰＲＦ／２に対応する最大血流速度Ｖmaxよりも大きな値となっており、折り返しが発生しているとともに、Ｖcurrentが折り返し領域にあると判定される。条件（１）は、条件（０）の符号が逆転した場合である。
【０１１０】
条件（２）に示すように、Ｖout-1が−Ｖthより小さい場合、Ｖcurrentが正の値となることは、想定される血流速度の変化の最大値を超える変化であるから、折り返しが発生している。また、ＶcurrentがＶzeroを挟む±Ｖthの範囲にあるため、Ｖcurrentは折り返し領域ではない。条件（３）は、条件（２）の符号が逆転した場合である。
【０１１１】
条件（０）から（３）のいずれも満たさない場合には、折り返しが発生しておらず、また、Ｖcurrentが折り返し領域にないと判定する。
【０１１２】
第１のパーシスタンス係数決定部１０６Ａは、折り返し判定部１０５から出力される２つの判定結果とフレームメモリ部１０４から読み出した血流速度データＶcurrentの絶対値に基づいて、第１のパーシスタンス係数参照メモリ部１０７への参照インデックスを作成する。表４に作成される参照インデックスを示す。
【０１１３】
【表４】
2011036891.tif 000018
【０１１４】
折り返しが発生しており、かつ、Ｖcurrentが折り返し領域にある場合、血流速度Ｖcurrentは、実際には、Ｖmaxをまたは−Ｖmaxを超える大きな値と考えられる。このため、参照インデックスはＶmaxとなる。その他の場合には、Ｖcurrentの絶対値Abs(Ｖcurrent)となる。
【０１１５】
第１のパーシスタンス係数参照メモリ部１０７Ａには、参照インデックスに対応付けられた第１のパーシスタンス係数によって構成される第１の参照テーブルが記憶されている。第１のパーシスタンス係数決定部１０６Ａは、第１のパーシスタンス係数参照メモリ部１０７Ａにアクセスし、作成された参照インデックスに対応付けられた第１のパーシスタンス係数を読み出し、第１のパーシスタンス演算部１０８Ａへ出力する。
【０１１６】
図６（ｂ）は、参照インデックスと第１のパーシスタンス係数との対応関係の一例を示すグラフである。図６（ｂ）において、横軸は参照インデックスを示し、縦軸は、パーシスタンス係数を示している。表４に示したように、参照インデックスは、ＶmaxまたはＶcurrentの絶対値Abs(Ｖcurrent)である。Ｖcurrentの絶対値が閾値Ｖth以下である場合には、Ｖcurrentの増大にともなって単調に増加する第１のパーシスタンス係数Ｃpersistenceが対応付けられている。つまり、Ｖcurrentの絶対値が閾値Ｖth以下である場合には、最新のフレームの血流速度Ｖcurrentに応じて異なるパーシスタンス係数Ｃpersistenceが対応付けられる。
【０１１７】
これに対し、第２のパーシスタンス係数決定部１０６Ｂは、折り返し判定部１０５から出力される２つの判定結果とパーシスタンスメモリ部１０９から読み出した血流速度データＶout-1の絶対値に基づいて、第１のパーシスタンス係数参照メモリ部１０７への参照インデックスを作成する。表５に作成される参照インデックスを示す。
【０１１８】
【表５】
2011036891.tif 000019
【０１１９】
第２のパーシスタンス係数決定部１０６Ｂは、条件（２）から（４）の場合に、パーシスタンスメモリ部１０９から読み出した血流速度データＶout-1の絶対値を参照インデックスとして生成する点で第１のパーシスタンス係数決定部１０６Ａと異なっている。
【０１２０】
第２のパーシスタンス係数参照メモリ部１０７Ｂには、参照インデックスに対応付けられた第２のパーシスタンス係数によって構成される第２の参照テーブルが記憶されている。第２のパーシスタンス係数決定部１０６Ｂは、第２のパーシスタンス係数参照メモリ部１０７Ｂにアクセスし、作成された参照インデックスに対応付けられた第２のパーシスタンス係数を読み出し、第２のパースタンス演算部１０８Ｂへ出力する。
【０１２１】
図６（ｃ）は、参照インデックスと第２のパーシスタンス係数との対応関係の一例を示すグラフである。図６（ｃ）において、横軸は参照インデックスを示し、縦軸は、パーシスタンス係数を示している。表５に示したように、参照インデックスは、ＶmaxまたはＶout-1の絶対値Abs(Ｖout-1)である。Ｖout-1の絶対値が閾値Ｖth以下である場合には、Ｖout-1の増大にともなって単調に増加する第２のパーシスタンス係数Ｃpersistenceが対応付けられている。つまり、Ｖout-1の絶対値が閾値Ｖth以下である場合には、１つ前の血流速度Ｖout-1に応じて異なる第２のパーシスタンス係数Ｃpersistenceが対応付けられる。
【０１２２】
図６（ｂ）および図６（ｃ）に示すように、参照インデックスがいずれの値をとっても第２のパーシスタンス係数のほうが第１のパーシスタンス係数より大きい。つまり、第１のパーススタンス係数は、最新のフレームの血流速度に対応づけられ、かつ、小さな値である。第１のパーススタンス係数が大きくなれば、１つ前のフレームの血流速度をより考慮した演算となるため、第１のパーシスタンス演算部１０８Ａは、残像効果を抑制し、すみやかに血流速度を変化させる演算を行う。これに対し、第２のパーススタンス係数は、１つ前のフレームの血流速度に対応付けられ、かつ、大きな値であるため、第２のパーシスタンス演算部１０８Ｂは、残像効果を高め、血流速度の変化を抑制する演算を行う。
【０１２３】
また、上述したように第１のパーシスタンス演算部１０８Ａは、残像効果を抑制し、速やかに血流速度を変化させる演算を行うため、血流速度は高いものの、血流パワーが小さいために血流を正しく検出できない場合、突然血流速度がゼロとなる可能性がある。この場合、血流速度が高くなるにつれて、血流画像に色調あるいは諧調の色付けを行っていると、血流画像が突然暗い諧調で色づけされ、画像が点滅するように表示される。このため、参照インデックスが増加するにつれて第１のパーシスタンス係数を単調に増加させ、血流速度が高くなるにつれて残像効果を高め、血流画像の点滅を抑制することができる。
【０１２４】
また、第２のパーシスタンス演算部１０８Ｂは、残像効果を高めた画像表示を行うため、血流速度が高くなるにつれて、血流画像に色調あるいは諧調の色付けを行っていると、血流速度が低い場合には、暗めの表示が必要以上に長い時間残像として表示されてしまう。例えば、探触子を動かした場合に、血流表示が尾引く印象を与えてしまう。このため、参照インデックスが増加するにつれて第２のパーシスタンス係数を単調に増加させ、血流速度が低くなるにつれて残像効果を抑制することができる。従って、血流速度の絶対値に基づく参照インデックスと第１および第２のパーシスタンス係数との間に適切な単調増加する関係を設定することにより、高品位な血流表示を実現することができる。
【０１２５】
なお、表３、表４および表５から分かるように、Ｖcurrent＜−Ｖthでも、Ｖout-1＞０であれば、参照インデックスはＶmaxとなる（条件（０））一方、Ｖout-1＜０であれば、参照インデックスはＶcurrentの絶対値Abs(Ｖcurrent)となる（条件（４））。このため、Ｖcurrent＜−Ｖthを満たしていても、Ｖout-1が正か負であるかによって、参照インデックスが異なり、パーシスタンス係数Ｃpersistenceも異なる。その結果、Ｖcurrent＜−Ｖthを満たす隣接する領域であっても、Ｖout-1が正か負であるかによってカラーマッピングフロー画像として表示される色が異なり、画像に不連続な色調部分が生じてしまう。
【０１２６】
このような不自然な表示を抑制するため、Ｖcurrentの絶対値が閾値Ｖth以上である場合には、参照インデックスに同じ値のパーシスタンス係数Ｃpersistenceが対応付けられていることが好ましい。これにより、折り返しが発生する血流領域、あるいはその境界近傍において、自然な表示を行うことができる。
【０１２７】
このようにして決定した第１および第２のパーシスタンス係数を用いて第１のパーシスタンス演算部１０８Ａおよび第２のパーシスタンス演算部１０８Ｂがそれぞれパーシスタンス処理の施された最新の血流速度データを生成する。
【０１２８】
最大値選択部１１６は、２つの血流速度データの絶対値の大きいほうを選択し、選択した血流速度データを残像処理された血流速度データとして出力する。つまり、絶対値の大きな血流速度データが得られるように、２つのパーシスタンス処理の結果を選択するため、甲状腺、肝臓、腎臓などに存在するような血流パワーが不安定となる末梢血管の血流表示を点滅させずに、かつ、平滑化に起因する末梢血管の血流消滅を起こさない血流動画を表示することが可能となる。
【０１２９】
なお、上記実施形態では、ＣＦＭフレームデータの血流速度に基づいてパーシスタンス係数を動的に決定し、血流速度にパーシスタンス演算を行っているが、上述したようにＣＦＭフレームデータの他のデータ、例えば、血流パワーデータに対してパーシスタンス演算を行ってもよいし、Ｂモード断層画像データに対してパーシスタンス演算を行ってもよい。
【０１３０】
また、上記実施形態では、最新のフレームおよび１つ前のフレームの血流速度データを用いてパーシスタンス処理を行っていたが、２つ前、あるいは３つ以上前のフレームの血流速度データも用いてパーシスタンス処理を行ってもよい。また、式（１）に限られず、他の演算式を用いてパーシスタンス処理を行ってもよい。
【産業上の利用可能性】
【０１３１】
本発明は、被検体の血流状態を表示することのできる超音波診断装置に好適に用いることができる。
【符号の説明】
【０１３２】
１０１、４０１  探触子
１０２、４０２  超音波送受部
１０３、４０３  ＣＦＭ信号処理部
１０４、４０４  フレームメモリ部
１０５  折り返し判定部
１０６  パーシスタンス係数決定部
１０６Ａ  第１のパーシスタンス係数決定部
１０６Ｂ  第２のパーシスタンス係数決定部
１０７  パーシスタンス係数参照メモリ部
１０７Ａ  第１のパーシスタンス係数参照メモリ部
１０７Ｂ  第２のパーシスタンス係数参照メモリ部
１０８、４０７  パーシスタンス演算部
１０８Ａ  第１のパーシスタンス演算部
１０８Ｂ  第２のパーシスタンス演算部
１０９  パーシスタンスメモリ部
１１０、４０８  ＣＦＭ  ＤＳＣ部
１１１、４０９  断層画像信号処理部
１１２、４１０  断層画像  ＤＳＣ部
１１３、４１１  画像合成部
１１４、４１２  表示部
１１５、１１５’  残像処理部
１１６  最大値選択部
４０５  フレームメモリ選択部
４０６  パーシスタンス係数設定部
【手続補正２】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
探触子を繰り返し駆動し、前記探触子を駆動することにより送信された超音波が被検体において反射することにより得られた反射エコーを前記探触子によって受信し、複数の受信信号を逐次生成する送受信部と、
前記複数の受信信号に基づき、各フレームにおける前記被検体中の血流部分の血流速度データを逐次生成するカラーフローマッピング信号処理部と、
前記各フレームにおける血流速度データに対して残像処理を行う残像処理部と、
前記受信信号に基づき、Ｂモード断層画像フレームデータを生成する断層画像信号処理部と、
前記残像処理された血流速度データおよび前記Ｂモード断層画像フレームデータを合成する画像合成部と
を備え、
前記残像処理部は、
最新のフレームの血流速度データと最新より以前のフレームの残像処理された血流速度データに基づいて折り返し判定を行い、前記折り返し判定の結果と最新および最新より以前のフレームの血流速度データに基づいて動的にパーシスタンス係数を変更する超音波診断装置。
【請求項２】
前記残像処理部は、
前記最新のフレームの血流速度データを記憶する第１のメモリ部と、
前記最新より以前のフレームの残像処理された血流速度データを記憶する第２のメモリ部と、
前記第１のメモリ部および前記第２のメモリ部から血流速度データをそれぞれ読み出し、前記折り返し判定を行う折り返し判定部と、
前記折り返し判定の結果および前記第１のメモリ部に記憶された血流速度データに応じてパーシスタンス係数を決定するパーシスタンス係数決定部と、
前記パーシスタンス係数および前記折り返し判定の結果に基づいて、前記第１のメモリ部に記憶された血流速度データに対しパーシスタンス演算を行い、演算結果を前記残像処理された血流速度データとして出力するパーシスタンス演算部とを含む請求項１に記載の超音波診断装置。
【請求項３】
前記折り返し判定部は、前記第１のメモリ部に記憶された血流速度データおよび前記第２のメモリ部に記憶された血流速度データと複数の閾値とを比較することによって、折り返しが発生しているかどうか、および、前記最新のフレームの血流速度データが折り返し領域にあるかどうかを判定する請求項２に記載の超音波診断装置。
【請求項４】
前記残像処理部は、前記血流速度の値に応じて対応付けられた異なる２つ以上の値のパーシスタンス係数を含む参照テーブルを記憶している第３のメモリ部をさらに含む、請求項２または３に記載の超音波診断装置。
【請求項５】
前記参照テーブルは、所定の値以上の血流速度に対して一定の値のパーシスタンス係数が対応付けられている請求項４に記載の超音波診断装置。
【請求項６】
前記残像処理部は、
前記最新のフレームの血流速度データを記憶する第１のメモリ部と、
前記最新より以前のフレームの残像処理された血流速度データを記憶する第２のメモリ部と、
前記第１のメモリ部および前記第２のメモリ部から血流速度データをそれぞれ読み出し、前記折り返し判定を行う折り返し判定部と、
前記折り返し判定の結果および前記第１のメモリ部に記憶された血流速度データに応じて第１のパーシスタンス係数を決定する第１のパーシスタンス係数決定部と、
前記第１のパーシスタンス係数および前記折り返し判定の結果に基づいて前記第１のメモリ部に記憶された血流速度データに対しパーシスタンス演算を行う第１のパーシスタンス演算部と、
前記折り返し判定の結果および前記第２のメモリ部に記憶された血流速度データに応じて第２のパーシスタンス係数を決定する第２のパーシスタンス係数決定部と、
前記第２のパーシスタンス係数および前記折り返し判定の結果に基づいて前記第１のメモリ部に記憶された血流速度データに対しパーシスタンス演算を行う第２のパーシスタンス演算部と、
前記第１のパーシスタンス演算部から出力される演算結果の絶対値および前記第２パーシスタンス演算部から出力される演算結果の絶対値を比較し、大きい方の演算結果を前記残像処理された血流速度データとして出力する最大値選択部と、
を含む請求項１に記載の超音波診断装置。
【請求項７】
前記折り返し判定部は、前記第１のメモリ部に記憶された血流速度データおよび前記第２のメモリ部に記憶された血流速度データと複数の閾値とを比較することによって、折り返しが発生しているかどうか、および、前記最新のフレームの血流速度データが折り返し領域にあるかどうかを判定する請求項６に記載の超音波診断装置。
【請求項８】
前記残像処理部は、
前記血流速度の値に応じて対応付けられた異なる２つ以上の値の第１のパーシスタンス係数を含む第１の参照テーブルを記憶している第３のメモリ部と、
前記血流速度の値に応じて対応付けられた異なる２つ以上の値の第２のパーシスタンス係数を含む第２の参照テーブルを記憶している第４のメモリ部と、
をさらに含む、請求項６または７に記載の超音波診断装置。
【請求項９】
前記第１の参照テーブルと前記第２の参照テーブルにおいて、同じ血流速度の値に応じて対応付けられた第１のパーシスタンス係数と第２のパーシスタンス係数とは互いに異なる値である請求項８に記載の超音波診断装置。
【請求項１０】
前記第１の参照テーブルは、所定の値以上の血流速度に対して一定の値のパーシスタンス係数が対応付けられている請求項８に記載の超音波診断装置。
【請求項１１】
前記最新より以前のフレームの残像処理された血流速度データは、前記最新より１つ前のフレームの残像処理された血流速度データである請求項１から１０のいずれかに記載の超音波診断装置。
【手続補正３】
【補正対象書類名】図面
【補正対象項目名】全図
【補正方法】変更
【補正の内容】
【図１】
2011036891.tif 000020
【図２】
2011036891.tif 000021
【図３】
2011036891.tif 000022
【図４】
2011036891.tif 000023
【図５】
2011036891.tif 000024
【図６】
2011036891.tif 000025
【図７】
2011036891.tif 000026
【国際調査報告】
2011036891.tif 000027
2011036891.tif 000028
2011036891.tif 000029
2011036891.tif 000030