(57)【要約】【課題】折り返しの発生の有無を精度よく判定し、折り返しが発生している場合にはその影響を好適に補正することができる超音波診断装置を提供する。【解決手段】折り返し判定部１０５は、最新のフレームの血流速度データＶｃｕｒｒｅｎｔと、最新より１つ前のフレームの血流速度データＶｏｕｔ−１と、に加えて、最新より以前のフレームにおいて判定された、当該フレームにおいて残像処理された血流速度データの折り返し判定結果Ｒａｌｉａｓｉｎｇ−１に基づいて折り返し判定を行う。パーシスタンス演算部１０８は、折り返しが発生した血流速度データに対して補正を行ってからパーシスタンス演算を行う。【選択図】図１
2018050925.tif 000002
【特許請求の範囲】
【請求項１】
探触子を繰り返し駆動し、前記探触子の駆動により送信された超音波が被検体において反射された反射エコーを前記探触子によって受信し、複数の受信信号を逐次生成する送受信部と、
前記複数の受信信号に基づき、各フレームにおける前記被検体中の血流部分の血流速度データを逐次生成するカラーフローマッピング信号処理部と、
前記各フレームにおける血流速度データに対して残像処理を行う残像処理部と、
を備え、
前記残像処理部は、
最新のフレームの血流速度データである第１血流速度データと、最新より以前のフレームにおいて残像処理された血流速度データである第２血流速度データと、最新より以前のフレームにおいて判定された、当該フレームにおいて残像処理された血流速度データの折り返しの有無を示す第１折り返し判定結果と、に基づいて、前記第１血流速度データと前記第２血流速度データに折り返しが発生したか否かの判定を行い、前記第１血流速度データの折り返し判定の結果である第２折り返し判定結果と、第２血流速度データの折り返し判定の結果である第３折り返し判定結果と、を出力する折り返し判定部と、
前記第２および第３折り返し判定結果と、残像効果を調整するための所定のパーシスタンス係数と、に基づいて、前記第１血流速度データおよび前記第２血流速度データを用いてパーシスタンス演算を行い、前記パーシスタンス演算の結果を最新のフレームにおいて残像処理された血流速度データである第３血流速度データとして出力し、かつ、前記第３血流速度データに折り返しが発生したか否かの判定を行い、その結果である第４折り返し判定結果を出力するパーシスタンス演算部と、
を備える、超音波診断装置。
【請求項２】
前記パーシスタンス演算部は、
前記第２および第３折り返し判定結果に基づいて、前記最新のフレームの血流速度データと前記最新より以前のフレームの残像処理された血流速度データを補正した後、前記パーシスタンス演算を行う、
請求項１に記載の超音波診断装置。
【請求項３】
前記残像処理部は、
前記第２および第３折り返し判定結果と、前記第１血流速度データと、に基づいて、前記パーシスタンス係数を決定するパーシスタンス係数決定部、
をさらに備える、
請求項１または２に記載の超音波診断装置。
【請求項４】
前記残像処理部は、
前記第１血流速度データを記憶する第１メモリ部と、
前記第２血流速度データと、前記第１折り返し判定結果と、を記憶する第２メモリ部と、
をさらに備え、
前記折り返し判定部は、前記第１メモリ部および前記第２メモリ部から、前記第１血流速度データと、前記第２血流速度データと、前記第１折り返し判定結果と、をそれぞれ読み出して前記折り返しの判定を行い、
前記パーシスタンス演算部は、前記第１メモリ部から前記第１血流速度データを、前記第２メモリ部から前記第２血流速度データと、前記第１折り返し判定結果と、を読み出して前記パーシスタンス演算を行い、前記パーシスタンス演算の結果である前記第３血流速度データと、前記第４折り返し判定結果と、を前記第２メモリ部に出力する、
請求項１から３のいずれか一項に記載の超音波診断装置。
【請求項５】
前記折り返し判定部は、前記第１血流速度データの大きさと、前記第２血流速度データの大きさと、所定のしきい値と、を比較した結果、および、前記第１折り返し判定結果に基づいて、前記第１血流速度データと前記第２血流速度データに折り返しが発生したか否かの判定を行う、
請求項１から４のいずれか一項に記載の超音波診断装置。
【請求項６】
前記残像処理部は、
前記第２および第３の折り返し判定結果と、前記第１血流速度データに基づいて第１パーシスタンス係数を決定する第１パーシスタンス係数決定部と、
前記第２および第３の折り返し判定結果と、前記第２血流速度データに基づいて第２パーシスタンス係数を決定する第２パーシスタンス係数決定部と、
前記第２および第３折り返し判定結果と、前記第１パーシスタンス係数と、に基づいて、前記第１血流速度データおよび前記第２血流速度データを用いてパーシスタンス演算を行う第１パーシスタンス演算部と、
前記第２および第３折り返し判定結果と、前記第２パーシスタンス係数と、に基づいて、前記第１血流速度データおよび前記第２血流速度データを用いてパーシスタンス演算を行う第２パーシスタンス演算部と、
前記第１パーシスタンス演算部から出力される演算結果の絶対値および前記第２パーシスタンス演算部から出力される演算結果の絶対値を比較し、大きい方の演算結果を前記第３の血流速度データとして出力し、かつ、前記第３血流速度データに折り返しが発生したか否かの判定を行い、その結果である第４折り返し判定結果を出力する最大値選択部と、
をさらに有する、請求項１または２に記載の超音波診断装置。
【請求項７】
前記残像処理部は、
前記第２および第３の折り返し判定結果と、前記第１血流速度データに基づいて第１パーシスタンス係数を決定する第１パーシスタンス係数決定部と、
前記第２および第３の折り返し判定結果と、前記第２血流速度データに基づいて第２パーシスタンス係数を決定する第２パーシスタンス係数決定部と、
前記第１血流速度データの絶対値と前記第２血流速度データの絶対値とを比較し、前記第１血流速度データの絶対値の方が大きい場合は前記第１パーシスタンス係数を、そうでない場合は前記第２パーシスタンス係数を選択するパーシスタンス係数選択部と、
をさらに有し、
前記パーシスタンス演算部は、前記第２および第３折り返し判定結果と、前記選択されたパーシスタンス係数と、に基づいて、前記第１血流速度データおよび前記第２血流速度データを用いてパーシスタンス演算を行う、
請求項１または２に記載の超音波診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、カラーフローマッピング時の残像処理を行う超音波診断装置に関する。
【背景技術】
【０００２】
超音波診断装置は、超音波を被検体に送信し、その反射エコーに含まれる情報を解析することにより、被検体内の画像を作成する。カラーフローマッピング（Color Flow Mapping、以下ＣＦＭと略す場合がある）と呼ばれる手法によって、被検体内の血流を画像化することが可能であり、医療分野全般において、血流状態を表示することのできる超音波診断装置が広く利用されている。
【０００３】
カラーフローマッピングはカラードップライメージング（Color Doppler Imaging、ＣＤＩ）とも呼ばれ、ドップラ効果を利用する手法である。超音波が血流を照射すると、ドップラ効果により反射エコーに血流速度に応じたドップラ偏移が生じる。このドップラ偏移の情報を直交検波によって検出し、ＭＴＩ（Moving Target Indicator）フィルタと呼ばれるハイパスフィルタ処理、自己相関処理およびノイズカット処理を施すことによって、血流速度に関する情報が得られる。得られた血流速度の情報を色情報に変換し、二次元的にＢモード断層画像（エコーの振幅を輝度として表した画像）に重畳することによって、被検体内の血流状態を好適に表示することができる。
【０００４】
血流から得られる反射エコーによる受信信号強度は、Ｂモード断層画像の生成に用いられる組織散乱体および組織境界から得られる反射エコーの受信信号強度と比較して小さい。このため、カラーフローマッピングにおける信号処理によって得られる血流速度および血流パワー（移動する血流量）は不安定になりやすい。
【０００５】
特に、観察したい部分の血流速度が遅い場合、あるいは、観察したい部分が末梢血管である場合、血流パワーが小さくなるため、本来、システムノイズおよび音響ノイズのみをカットするべきノイズカット処理において、血流速度あるいは血流パワーに関する情報が除去されてしまいやすい。その結果、血流画像において本来血流として表示される部分が黒く抜けてしまう現象が発生する。具体的には、例えば、１秒間に数フレームから数十フレームの割合で、被検体内の血流を画像化する場合、そのうちのいくつかのフレームにおいて、血流部分が黒く表示されてしまう。このため、断層画像中の血流部分が突然消滅し、画像が滑らかでなくなったり、違和感のあるものとなったりする。
【０００６】
このような問題を解決した超音波診断装置として、信号処理後段においてパーシスタンス処理（残像処理）と呼ばれる時間方向補間を行う超音波診断装置が例えば特許文献１に開示されている。以下では、パーシスタンス処理を行う超音波診断装置の具体的な構成例について簡単に説明する。
【０００７】
図８に示す従来の超音波診断装置において、超音波送受信部４０２は、探触子４０１を駆動し、超音波を被検体へ送信する。また、探触子４０１によって被検体で生じる反射エコーを受信し、受信信号を生成する。Ｂモード断層画像を生成する場合、超音波送受信部４０２は、Ｂモード断層画像の生成に適した送受信を行い、得られた受信信号を断層画像信号処理部４０９へ出力する。カラーフローマッピング断層画像を生成する場合は、カラーフローマッピング断層画像の生成に適した送受信を行い、得られた受信信号をカラーフローマッピング処理部４０３（以下、ＣＦＭ信号処理部と略す）へ出力する。一般的に、カラーフローマッピング断層画像を生成する場合、安定したカラーフローマッピング断層画像を得るために、超音波送受信部４０２は、同じ音響線上において超音波の送受信を複数回行う。
【０００８】
ＣＦＭ信号処理部４０３は、受信信号に対して直交検波処理、ＭＴＩフィルタ処理および自己相関処理を行い、血流速度および血流パワーを算出したのち、システムあるいは音響ノイズを排除するノイズカット処理を行い、フレームメモリ部４０４に血流速度および血流パワーを出力する。
【０００９】
フレームメモリ部４０４は、リングバッファで構成されており、現在の走査フレームからＮフレーム前（Ｎは、１以上の整数）までの血流速度および血流パワーをフレーム単位で保存している。ここで、フレームとは、１画面のＣＦＭ断層画像を構成する血流速度データおよび血流パワーデータ群を示す。
【００１０】
フレームメモリ選択部４０５は、フレームメモリ部４０４から、あらかじめ設定された複数のＣＦＭフレームデータを選択し、パーシスタンス演算部４０７へ出力するようフレームメモリ部４０４に指令を出力する。パーシスタンス演算部４０７は、フレームメモリ部４０４から読み出されたＣＦＭフレームデータとパーシスタンス係数設定部４０６から出力されるパーシスタンス係数に基づいてパーシスタンス演算を実施し、ＣＦＭ  ＤＳＣ（Digital Scan Converter）部４０８に出力する。パーシスタンス演算は、単純な重み付け演算であり、パーシスタンス係数設定部４０６から出力されるパーシスタンス係数は、あらかじめシステムから設定された固定係数である。パーシスタンス演算とは、現在およびそれに連続する複数フレームの画像に対し重み付け加算を行い、その結果を表示に利用することによって、残像効果を与える演算である。
【００１１】
ＣＦＭ  ＤＳＣ部４０８は、パーシスタンス演算部４０７から出力されるＣＦＭフレームデータの座標を変換し、画像合成部４１１へ出力する。
【００１２】
断層画像信号処理部４０９は受信信号に対してダイナミックフィルタ処理を施すことにより、不要なノイズをカットしたのち、包絡線検波処理およびダイナミックレンジ圧縮処理を施して断層画像ＤＳＣ部４１０に断層画像フレームデータを出力する。断層画像ＤＳＣ部４１０は、断層画像信号処理部４０９からの断層画像フレームデータの座標を変換し、画像合成部４１１へ出力する。
【００１３】
画像合成部４１１は、ＣＦＭ  ＤＳＣ部４０８および断層画像ＤＳＣ部４１０から出力される各フレームデータを画素ごとに合成し、合成画像フレームデータを生成する。具体的には、血流速度がゼロである場合は、断層画像フレームデータを表示し、そうでない場合は、ＣＦＭフレームデータを表示するように２つのデータを画素ごと、あるいは、対応する測定点のデータごとに合成する。また、血流速度や血流の方向に応じてデータを色情報に変換し、表示部４１２に出力する。表示部４１２は、画像合成部４１１から受け取ったデータを表示する。
【先行技術文献】
【特許文献】
【００１４】
【特許文献１】特開平２−２８６１４０号公報
【特許文献２】特許第５６５２３９５号公報
【発明の概要】
【発明が解決しようとする課題】
【００１５】
ところで、ドップラ偏移によって計測できる血流速度はパルス繰り返し周波数（ＰＲＦ）による制限を受ける。その結果、±ＰＲＦ／２を超える周波数変化に対応する血流速度が、反対向きの血流速度として観測される現象である折り返しが発生し、血流速度を正確に評価することが困難となる。
【００１６】
このように折り返しが発生する場合には、良好な画像を得るために、折り返しが発生しているか否かを判定し、その影響を補正する必要がある。折り返しの発生を検出し、補正を行う技術として、例えば特許文献２に開示された技術がある。特許文献２には、血流速度および折り返しの状態に基づいてパーシスタンス処理に用いるパーシスタンス係数を動的に決定した後にＣＦＭフレームデータに対してパーシスタンス演算を行うことにより、折り返しの影響を補正する技術が開示されている。しかしながら、折り返しが発生しているか否かを判定する精度をより向上させたいという要望がある。
【００１７】
本発明は、折り返しの発生の有無を精度よく判定し、折り返しが発生している場合にはその影響を好適に補正することができる超音波診断装置を提供することを目的とする。
【課題を解決するための手段】
【００１８】
本発明の超音波診断装置は、探触子を繰り返し駆動し、前記探触子の駆動により送信された超音波が被検体において反射された反射エコーを前記探触子によって受信し、複数の受信信号を逐次生成する送受信部と、前記複数の受信信号に基づき、各フレームにおける前記被検体中の血流部分の血流速度データを逐次生成するカラーフローマッピング信号処理部と、前記各フレームにおける血流速度データに対して残像処理を行う残像処理部と、を備え、前記残像処理部は、最新のフレームの血流速度データである第１血流速度データと、最新より以前のフレームにおいて残像処理された血流速度データである第２血流速度データと、最新より以前のフレームにおいて判定された、当該フレームにおいて残像処理された血流速度データの折り返しの有無を示す第１折り返し判定結果と、に基づいて、前記第１血流速度データと前記第２血流速度データに折り返しが発生したか否かの判定を行い、前記第１血流速度データの折り返し判定の結果である第２折り返し判定結果と、第２血流速度データの折り返し判定の結果である第３折り返し判定結果と、を出力する折り返し判定部と、前記第２および第３折り返し判定結果と、残像効果を調整するための所定のパーシスタンス係数と、に基づいて、前記第１血流速度データおよび前記第２血流速度データを用いてパーシスタンス演算を行い、前記パーシスタンス演算の結果を最新のフレームにおいて残像処理された血流速度データである第３血流速度データとして出力し、かつ、前記第３血流速度データに折り返しが発生したか否かの判定を行い、その結果である第４折り返し判定結果を出力するパーシスタンス演算部と、を備える。
【発明の効果】
【００１９】
本発明によれば、折り返しの発生の有無を精度よく判定し、折り返しが発生している場合にはその影響を好適に補正することができる超音波診断装置を提供することができる。
【図面の簡単な説明】
【００２０】
【図１】本発明の第１の実施の形態に係る超音波診断装置の構成の一例を示すブロック図
【図２】ドップラ効果に基づいて得られる血流の速度と血流の方向とを図示した概要図
【図３Ａ】条件（１０）のＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１を複素平面上に示した図
【図３Ｂ】条件（１１）のＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１を複素平面上に示した図
【図３Ｃ】条件（１２）のＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１を複素平面上に示した図
【図３Ｄ】条件（１３）のＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１を複素平面上に示した図
【図３Ｅ】条件（１４）のＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１を複素平面上に示した図
【図３Ｆ】条件（１５）のＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１を複素平面上に示した図
【図４】参照インデックスとパーシスタンス係数との対応関係の一例を示すグラフ
【図５】本発明の第２の実施の形態に係る超音波診断装置の構成の一例を示すブロック図
【図６Ａ】第１参照インデックスと第１パーシスタンス係数との対応関係の一例を示すグラフ
【図６Ｂ】第２参照インデックスと第２パーシスタンス係数との対応関係の一例を示すグラフ
【図７】本発明の第３実施の形態に係る超音波診断装置の構成の一例を示すブロック図
【図８】従来の超音波診断装置の構成の一例を示すブロック図
【発明を実施するための形態】
【００２１】
（第１の実施の形態）
＜超音波診断装置の構成例＞
以下、図面を参照しながら本発明の第１の実施の形態に係る超音波診断装置について説明する。図１は、本発明の第１の実施の形態に係る超音波診断装置の構成の一例を示すブロック図である。
【００２２】
図１に示す超音波診断装置１１は、探触子１０１と、超音波送受信部１０２と、ＣＦＭ信号処理部１０３と、残像処理部１１５と、断層画像信号処理部１１１と、ＣＦＭ  ＤＳＣ部１１０と、断層画像ＤＳＣ部１１２と、画像合成部１１３と、表示部１１４とを備える。これらの構成のうち、探触子１０１および表示部１１４には汎用の探触子および表示装置を用いることが可能であり、超音波診断装置１１は、探触子１０１および表示部１１４を備えていなくてもよい。
【００２３】
超音波送受信部１０２は、探触子１０１を駆動する駆動信号を生成して探触子１０１へ出力する。これに応じて、探触子１０１は被検体に向けて超音波を送信する。また、探触子１０１は送信された超音波が被検体において反射することによって得られた反射エコーを受信し、受信信号を生成する。超音波送受信部１０２は、本発明の送受信部に対応する。
【００２４】
より具体的には、探触子１０１は複数の圧電素子を含み、各圧電素子から超音波が送信される。複数の圧電素子から送信された超音波により超音波ビームが構成される。超音波送受信部１０２は、被検体を複数の超音波ビームによって走査するように、各圧電素子の遅延制御を行いながら探触子１０１を駆動する。
【００２５】
被検体による反射エコーは各圧電素子によって受信される。超音波送受信部１０２は、各圧電素子の遅延制御を行うことにより、送信された超音波ビームに対応する受信信号を探触子１０１に生成させる。超音波ビームが被検体を１回走査すると、１フレーム分のデータが得られる。探触子１０１は、１秒間に数回〜数十回繰り返して被検体を走査することにより、毎秒数フレームから数十フレームの受信信号を逐次生成することができる。
【００２６】
本第１の実施の形態に係る超音波診断装置１１は、Ｂモード断層画像とカラーフローマッピング画像と生成し、これらを合成して表示部１１４に表示する。このため、超音波送受信部１０２の上述した超音波の送受信は、Ｂモード断層画像の生成およびカラーフローマッピング画像の生成のそれぞれについて行われる。Ｂモード断層画像の１秒当たりのフレーム数およびカラーフローマッピング画像の１秒当たりのフレーム数は同じであってもよいし、異なっていてもよい。フレーム数が同じである場合には、Ｂモード断層画像生成のための超音波の送受信およびカラーフローマッピング画像生成のための超音波の送受信を交互に繰り返し行ってもよい。
【００２７】
Ｂモード断層画像を生成する場合、超音波送受信部１０２は、Ｂモード断層画像の生成に適した送受信を行い、得られた受信信号を断層画像信号処理部１１１へ出力する。一方、カラーフローマッピング断層画像を生成する場合、超音波送受信部１０２は、カラーフローマッピング断層画像の生成に適した送受信を行い、得られた受信信号をＣＦＭ信号処理部１０３へ出力する。一般的に、カラーフローマッピング断層画像を生成する場合、安定したカラーフローマッピング断層画像を得るために、超音波送受信部１０２は、同じ音響線上において超音波の送受信を複数回行う。
【００２８】
ＣＦＭ信号処理部１０３は、受信信号に対して直交検波処理、ＭＴＩフィルタ処理および自己相関処理を行った後、血流速度および血流パワーを算出したのち、システムあるいは音響ノイズを排除するノイズカット処理を行ってＣＦＭフレームデータを生成する。ＣＦＭフレームデータは少なくとも血流速度データを含んでおり、その他に血流パワーデータや血流速度の分散データを含んでいてもよい。ＣＦＭ信号処理部１０３はこの処理を、各フレームを構成する受信信号ごとに逐次繰り返して行う。ＣＦＭ信号処理部１０３が生成したＣＦＭフレームデータは、フレームごとに残像処理部１１５へ出力される。ＣＦＭ信号処理部１０３は、本発明のカラーフローマッピング信号処理部に対応する。
【００２９】
残像処理部１１５は、パーシスタンス係数を用いて、フレームごとにＣＦＭフレームデータに対して残像処理を行う。本第１の実施の形態に係る超音波診断装置１１では、残像処理部１１５は、最新のフレームの血流速度に応じてパーシスタンス係数を決定する。つまり、パーシスタンス係数は一定ではなく、最新のフレームの血流速度に基づく動的な値である。これにより、血流速度に応じてパーシスタンス係数を変化させ、残像効果を調整することができる。ただし、ドップラ偏移によって計測できる血流速度はパルス繰り返し周波数（ＰＲＦ）による制限を受ける。その結果、±ＰＲＦ／２を超える周波数変化に対応する血流速度が、反対向きの血流速度として観測される現象である折り返しが発生し、血流速度を正確に評価することが困難となる。このため、残像処理部１１５は、折り返しが発生しているか否かを判定し、折り返しが発生していると判定した場合には、血流速度データに対して補正を行い、折り返しの影響をキャンセルする。これにより、本第１の実施の形態に係る超音波診断装置１１では、折り返しの発生に拘わらず正確な血流速度を評価することができる。
【００３０】
本第１の実施の形態に係る超音波診断装置１１では、残像処理部１１５は、最新のフレームの血流速度データと、１つ前のフレームの血流速度データと、を用いて、折り返しが生じているか否かを判定する。このような処理を行うための構成として、残像処理部１１５は、フレームメモリ部（第１メモリ部）１０４と、折り返し判定部１０５と、パーシスタンス係数決定部１０６と、パーシスタンス係数参照メモリ部（第３メモリ部）１０７と、パーシスタンス演算部１０８と、パーシスタンスメモリ部（第２メモリ部）１０９と、を有する。
【００３１】
フレームメモリ部１０４は、最新の（現在走査中の）フレームのＣＦＭフレームデータを記憶する。パーシスタンスメモリ部１０９は、最新より１つ前のフレームのパーシスタンス演算部１０８の出力結果であるＣＦＭフレームデータと、パーシスタンス演算部１０８から出力された、最新より１つ前のフレームの折り返し判定結果と、を記憶する。パーシスタンスメモリ部１０９に記憶されるＣＦＭフレームデータには、パーシスタンス演算、すなわち残像処理が施されている。以下では、フレームメモリ部１０４に記憶されたＣＦＭフレームデータのうちの血流速度データをＶｃｕｒｒｅｎｔ、パーシスタンスメモリ部１０９に記憶されたＣＦＭフレームデータのうちの血流速度データをＶｏｕｔ−１、パーシスタンスメモリ部１０９に記憶された、最新より１つ前のフレームの折り返し判定結果をＲａｌｉａｓｉｎｇ−１と呼ぶ。
【００３２】
折り返し判定部１０５は、フレームメモリ部１０４からＣＦＭフレームデータのうちの血流速度データＶｃｕｒｒｅｎｔを、パーシスタンスメモリ部１０９からＣＦＭフレームデータのうちの血流速度データＶｏｕｔ−１をそれぞれ読み出し、それぞれの血流速度データに対して折り返し判定を行う。より具体的には、折り返し判定部１０５は、血流速度データＶｃｕｒｒｅｎｔおよび血流速度データＶｏｕｔ−１と複数のしきい値とを比較することによって、最新フレームの血流速度データＶｃｕｒｒｅｎｔと、最新より１つ前の血流速度データＶｏｕｔ−１のそれぞれに対して折り返しが発生しているか否かを判定し、判定結果をパーシスタンス係数決定部１０６およびパーシスタンス演算部１０８に出力する。
【００３３】
パーシスタンス係数決定部１０６は、折り返し判定部１０５からの２つの判定結果（Ｖｃｕｒｒｅｎｔの判定結果とＶｏｕｔ−１の判定結果）と、フレームメモリ部１０４から読み出した血流速度データＶｃｕｒｒｅｎｔと、に基づいて、パーシスタンス係数参照メモリ部１０７への参照インデックスを作成する。また、パーシスタンス係数決定部１０６は、パーシスタンス係数参照メモリ部１０７にアクセスし、参照インデックスに対応付けられたパーシスタンス係数を読み出して、パーシスタンス演算部１０８に設定する。パーシスタンス係数参照メモリ部１０７には、あらかじめ血流速度の値に対応付けられたパーシスタンス係数の参照テーブルが記憶されている。この参照テーブルは、血流速度の値に応じて対応付けられた異なる２つ以上の値のパーシスタンス係数を含む。
【００３４】
パーシスタンス演算部１０８は、パーシスタンス係数決定部１０６から設定されたパーシスタンス係数と、折り返し判定部１０５からの２つの判定結果と、に基づいて、最新の血流速度データＶｃｕｒｒｅｎｔに対してパーシスタンス演算を行う。パーシスタンス演算とは、現在およびそれに連続する複数フレームの画像に対し重み付け加算を行い、その結果を表示に利用することによって、残像効果を与える演算である。
【００３５】
なお、残像処理部１１５における、折り返し判定部１０５、パーシスタンス係数決定部１０６、およびパーシスタンス演算部１０８が行う処理の詳細は後述する。
【００３６】
ＣＦＭ  ＤＳＣ部１１０は、パーシスタンス演算部１０８から出力される血流速度データの座標を変換し、画像合成部１１３へ出力する。
【００３７】
断層画像信号処理部１１１は、受信信号に対してダイナミックフィルタ処理を施すことにより、不要なノイズをカットしたのち、包絡線検波処理およびダイナミックレンジ圧縮処理を施して断層画像ＤＳＣ部１１２に断層画像フレームデータを出力する。断層画像ＤＳＣ部１１２は、断層画像信号処理部１１１からの断層画像フレームデータの座標を変換し、画像合成部１１３へ出力する。
【００３８】
画像合成部１１３は、ＣＦＭ  ＤＳＣ部１１０および断層画像ＤＳＣ部１１２から出力される各フレームデータを画素ごとに、あるいは、対応する測定点のデータごとに合成し、合成画像フレームデータを生成する。具体的には、血流速度がゼロである場合は、断層画像フレームデータを表示し、そうでない場合は、ＣＦＭフレームデータを表示するように２つのデータを画素ごとに、あるいは、対応する測定点のデータごとに合成する。また、血流速度や血流の方向に応じてデータを色情報に変換し、表示部１１４に出力する。表示部１１４は、画像合成部１１３から受け取ったデータを表示する。
【００３９】
＜折り返し判定部の処理＞
次に、折り返し判定部１０５における折り返し判定方法と、パーシスタンス係数決定部１０６におけるパーシスタンス係数の決定方法と、について、より詳細に説明する。
【００４０】
折り返し判定部１０５は、フレームメモリ部１０４から最新のフレームのＣＦＭフレームデータに含まれる血流速度データＶｃｕｒｒｅｎｔを読み出し、パーシスタンスメモリ部１０９から最新の１フレーム前における残像処理された血流速度データＶｏｕｔ−１を読み出す。そして、折り返し判定部１０５は、読み出したＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１の値を用いて、Ｖｃｕｒｒｅｎｔに折り返しが発生しているかの判定と、Ｖｏｕｔ−１に折り返しが発生しているかの判定と、を行う。
【００４１】
図２は、ドップラ効果に基づいて得られる血流の速度と血流の方向とを図示した概要図である。図２は複素平面を示した図であって、血流速度は複素平面上において位相角を有するベクトルとして表すことができる。すなわち、血流速度が増減するということは、図２に示す複素平面上において血流速度に対応するベクトルが原点を中心に回転することを意味する。
【００４２】
図２に示すベクトルＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１は、ＣＦＭ信号処理部１０３が生成した、あるフレームにおけるＣＦＭフレームデータに含まれる血流速度データＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１が複素平面上に表されたものである。上述したように、図２に示す複素平面において、血流速度データＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１の大きさは、＋Ｘ軸からの回転角（すなわち位相角）に対応している。
【００４３】
図２に示す複素平面において、血流速度が正の値である場合、その血流速度に対応するベクトルは第１象限または第２象限に位置し、血流速度が負の値である場合、その血流速度に対応するベクトルは第３象限または第４象限に位置することになる。なお、血流速度が正の値であるとは、血流が探触子１０１に近づく方向に移動していることを意味し、血流速度が負の値であるとは、血流が探触子１０１から遠ざかる方向に移動していることを意味する。
【００４４】
複素平面において、ある血流速度に対応する速度ベクトルの大きさが０から徐々に増大すると、血流速度データに対応するベクトルは、＋Ｘ軸からの回転角が大きくなる方向に移動する。具体的には、＋Ｘ軸から−Ｘ軸に向かって第１象限、第２象限の順に移動する。そして、血流速度が所定の値＋Ｖｍａｘまたは−Ｖｍａｘ（±ＰＲＦ／２に相当する周波数変化に対応する速度）を超え、−Ｘ軸を超えて第２象限から第３象限へ移動すると、折り返しが生じ、正の値であるはずの血流速度が負の値として測定されてしまう。なお、速度ベクトルの大きさが０から徐々に減少する場合は、血流速度データに対応するベクトルは、上記と反対の方向に移動する。
【００４５】
以上のことを踏まえて、折り返し判定部１０５は、以下のように最新のフレームの血流速度データＶｃｕｒｒｅｎｔの大きさと、最新より１つ前のフレームの血流速度データＶｏｕｔ−１の大きさと、所定のしきい値Ｖｔｈを比較することで、ＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１に折り返しが発生しているか否かを判定する。
【００４６】
ここで、まず折り返し判定の方法として、上記した特許文献２に開示されている方法について説明し、その方法において改良すべき点（要改良点）を挙げた後、本発明が採用する折り返し判定の方法について説明する。
【００４７】
［従来の折り返し判定方法（特許文献２に開示された方法）］
以下の表１は、従来の折り返し方法を用いた場合の折り返し判定結果を示している。なお、表１において、しきい値Ｖｔｈは、隣接するフレームの時間間隔において想定される血流速度の変化の最大値である。
【００４８】
【表１】
2018050925.tif 000003
【００４９】
表１の条件（０）に示すように、最新より１つ前のフレームの血流速度データＶｏｕｔ−１が正であり、かつ、最新のフレームの血流速度データＶｃｕｒｒｅｎｔが−Ｖｔｈよりも小さい場合、仮に折り返しが発生していないとすると、最新より１つ前のフレームでは正方向（探触子１０１に近づく方向）であった血流が、最新のフレームでは逆の負方向（探触子１０１から遠ざかる方向）の高速な血流に急激に変化したことになるため不自然である。すなわち、最新より１つ前のフレームで正方向であった血流の速度が最新のフレームで大きくなり、＋ＰＲＦ／２に対応する血流速度、すなわち折り返しが発生する血流速度よりも大きな値となっていると推定される。このため、条件（０）が満たされるとき、Ｖｃｕｒｒｅｎｔに折り返しが発生していると判定される。条件（１）は、条件（０）の符号が逆転した場合である。なお、以下では、折り返しが発生する血流速度をＶｍａｘまたは−Ｖｍａｘと記載する。
【００５０】
条件（２）に示すように、Ｖｏｕｔ−１が−Ｖｔｈより小さい場合、Ｖｃｕｒｒｅｎｔが正の値となったとすると、仮に折り返しが発生していないとすると、最新より１つ前のフレームでは負方向であった高速な血流が、最新のフレームでは逆の正方向の血流に急激に変化したことになるため不自然である。すなわち、Ｖｏｕｔ−１は折り返しが発生していて正方向の＋ＰＲＦ／２より大きな値であり、最新のフレームでは同じ正方向で血流の速度が遅くなったと推定される。従って、折り返しが発生していると判定される。条件（３）は、条件（２）の符号が逆転した場合である。
【００５１】
条件（０）から（３）のいずれも満たさない場合には、折り返しが発生していないと判定される。
【００５２】
［従来の折り返し判定方法の要改良点］
しかしながら、上記した従来の折り返し判定方法には、以下のような要改良点が存在する。第１の要改良点は、折り返しが発生していると判定したとしても、ＶｃｕｒｒｅｎｔとＶｏｕｔ−１のどちらに折り返しが発生しているのかが特定されない、という点である。
【００５３】
具体的には、例えばＶｏｕｔ−１＞ＶｔｈかつＶｃｕｒｒｅｎｔ＜−Ｖｔｈである場合、表１では、条件（０）または条件（３）に該当し、折り返しが発生していると判定されるが、ＶｃｕｒｒｅｎｔとＶｏｕｔ−１のどちらが折り返しているのかを特定することができない場合がある。すなわち、条件（０）ならＶｃｕｒｒｅｎｔに折り返しが発生していることになり、条件（３）ならＶｏｕｔ−１に折り返しが発生していることになり、どちらが折り返しているのかを特定することができない。また、Ｖｏｕｔ−１＜−ＶｔｈかつＶｃｕｒｒｅｎｔ＜Ｖｔｈである場合も、条件（１）または条件（２）に該当し、同様にどちらが折り返しているのかを特定することができない場合がある。
【００５４】
第２の要改良点は、折り返しが発生していないと判定されたとしても、実際には折り返しが発生している場合がある、という点である。具体的には、例えばＶｏｕｔ−１＜−ＶｔｈかつＶｃｕｒｒｅｎｔ＜−Ｖｔｈである場合、上記表１では条件（４）に該当し、折り返しが発生していないと判定されるが、実際には、ＶｃｕｒｒｅｎｔとＶｏｕｔ−１の両方に折り返しが発生している場合にもＶｏｕｔ−１＜−ＶｔｈかつＶｃｕｒｒｅｎｔ＜−Ｖｔｈが満たされうる。Ｖｏｕｔ−１＞ＶｔｈかつＶｃｕｒｒｅｎｔ＞Ｖｔｈである場合も同様である。
【００５５】
以下説明する本発明の折り返し判定方法は、このような従来の折り返し判定方法を改良するものである。
【００５６】
［本発明の折り返し判定部による折り返し判定方法］
本発明においては、折り返し判定部１０５は、最新のフレームの血流速度データＶｃｕｒｒｅｎｔの大きさと、最新より１つ前のフレームの血流速度データＶｏｕｔ−１の大きさと、最新より１つ前のフレームの折り返し判定結果Ｒａｌｉａｓｉｎｇ−１と、に基づいて、折り返し判定を行う。以下の表２は、折り返し判定部１０５による折り返し判定結果を示している。なお、表２において、「折り返し有り」とは、最新より１つ前のフレームの折り返し判定結果Ｒａｌｉａｓｉｎｇ−１が、折り返しがあったという判定結果であることを意味し、「折り返し無し」とは、最新より１つ前のフレームの折り返し判定結果Ｒａｌｉａｓｉｎｇ−１が、折り返しがなかったという判定結果であることを意味する。
【００５７】
【表２】
2018050925.tif 000004
【００５８】
表２の条件（１０）に示すように、Ｖｏｕｔ−１＜−ＶｔｈかつＶｃｕｒｒｅｎｔ＜−Ｖｔｈである場合、上記したように、実際には、ＶｃｕｒｒｅｎｔとＶｏｕｔ−１の両方に折り返しが発生していない、あるいは、両方に折り返しが発生している。このため、折り返し判定部１０５は、Ｒａｌｉａｓｉｎｇ−１を参照し、最新より１フレーム前の折り返し判定結果Ｒａｌｉａｓｉｎｇ−１が条件（１０）に示すように「折り返し有り」であった場合には、ＶｃｕｒｒｅｎｔとＶｏｕｔ−１にも折り返しが発生したと判定する。その理由は以下の通りである。最新より１フレーム前の折り返し判定結果Ｒａｌｉａｓｉｎｇ−１が「折り返し有り」ということはＶｏｕｔ−１に折り返しが発生しているということを意味するため、Ｖｏｕｔ−１は＋ＰＲＦ／２より大きな値である。このとき、仮にＶｃｕｒｒｅｎｔに折り返しが発生していないとすると、最新より１つ前のフレームでは正方向の＋ＰＲＦ／２に相当する速度を超える非常に高速であった血流が、最新のフレームでは逆の負方向の高速な血流に急激に変化したことになるため不自然である。このため、最新のフレームでも正方向の＋ＰＲＦ／２に相当する速度を超える非常に高速な血流のままと推定される。条件（１１）に示す、Ｖｏｕｔ−１＜−ＶｔｈかつＶｃｕｒｒｅｎｔ＜−Ｖｔｈである場合も同様である。
【００５９】
図３Ａは、条件（１０）のＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１を図２と同様の複素平面上に示した図である。また、図３Ｂは、条件（１１）のＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１を図２と同様の複素平面上に示した図である。
【００６０】
また、条件（１２）に示すように、Ｖｏｕｔ−１＞０かつＶｃｕｒｒｅｎｔ＜−Ｖｔｈで、さらに折り返し判定結果Ｒａｌｉａｓｉｎｇ−１が「折り返し無し」であった場合には、Ｖｏｕｔ−１に折り返しが発生していないということであるから、表１の条件（０）と同様の理由によって、折り返し判定部１０５は、Ｖｃｕｒｒｅｎｔに折り返しが発生しており、Ｖｏｕｔ−１には折り返しが発生していないと判定する。条件（１３）は条件（１２）の符号が逆転した場合である。
【００６１】
図３Ｃは、条件（１２）のＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１を図２と同様の複素平面上に示した図である。また、図３Ｄは、条件（１３）のＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１を図２と同様の複素平面上に示した図である。
【００６２】
条件（１４）に示すように、Ｖｏｕｔ−１が−Ｖｔｈより小さい場合に、Ｖｃｕｒｒｅｎｔが正の値となり、かつ、折り返し判定結果Ｒａｌｉａｓｉｎｇ−１が「折り返し有り」であった場合には、Ｖｏｕｔ−１に折り返しが発生しているということであるから、表１の条件（２）と同様の理由によって、折り返し判定部１０５は、Ｖｃｕｒｒｅｎｔに折り返しが発生しておらず、Ｖｏｕｔ−１に折り返しが発生していると判定する。条件（１５）は条件（１４）の符号が逆転した場合である。
【００６３】
図３Ｅは、条件（１４）のＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１を図２と同様の複素平面上に示した図である。また、図３Ｆは、条件（１５）のＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１を図２と同様の複素平面上に示した図である。
【００６４】
ここで、Ｖｏｕｔ−１＞ＶｔｈかつＶｃｕｒｒｅｎｔ＜−Ｖｔｈである場合、折り返し判定結果Ｒａｌｉａｓｉｎｇ−１が「折り返し無し」であれば条件（１２）に該当し、折り返し判定結果Ｒａｌｉａｓｉｎｇ−１が「折り返し有り」であれば条件（１５）に該当するから、それぞれＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１の折り返しの有無を特定することができる。また、Ｖｏｕｔ−１＜−ＶｔｈかつＶｃｕｒｒｅｎｔ＞Ｖｔｈである場合も、同様に条件（１３）または条件（１４）に該当し、ＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１の折り返しの有無を特定することができる。
【００６５】
そして、条件（１０）から（１５）のいずれも満たさない場合には、折り返し判定部１０５は、折り返しが発生していないと判定する（条件（１６））。
【００６６】
このように、本発明の折り返し判定部１０５による折り返し判定方法によれば、従来の折り返し判定方法では折り返しが発生しているか否かを正確に判定できない場合でも、折り返しが発生しているか否かを正確に判定することができる。
【００６７】
＜パーシスタンス係数決定部の処理＞
次に、パーシスタンス係数決定部１０６がパーシスタンス計数をパーシスタンス演算部１０８に設定する方法について説明する。
【００６８】
パーシスタンス係数決定部１０６は、折り返し判定部１０５から出力された、最新のフレームの血流速度データＶｃｕｒｒｅｎｔの折り返し判定結果と、最新より１つ前のフレームの血流速度データＶｏｕｔ−１の折り返し判定結果と、フレームメモリ部１０４から読み出したＶｃｕｒｒｅｎｔの絶対値と、に基づいて、パーシスタンス係数参照メモリ部１０７への参照インデックスを作成する。以下の表３は、パーシスタンス係数決定部１０６が作成される参照インデックスを示す。
【００６９】
【表３】
2018050925.tif 000005
【００７０】
Ｖｃｕｒｒｅｎｔに折り返しが発生している場合、血流速度データＶｃｕｒｒｅｎｔの大きさは、実際には、Ｖｍａｘまたは−Ｖｍａｘを超える絶対値の大きな値と考えられる。このため、参照インデックスはＶｍａｘとなる。その他の場合には、Ｖｃｕｒｒｅｎｔの絶対値Abs(Ｖｃｕｒｒｅｎｔ)となる。
【００７１】
パーシスタンス係数参照メモリ部１０７には、参照インデックスに対応付けられたパーシスタンス係数によって構成される参照テーブルが記憶されている。パーシスタンス係数決定部１０６は、パーシスタンス係数参照メモリ部１０７にアクセスし、作成された参照インデックスに対応付けられたパーシスタンス係数を読み出し、パーシスタンス演算部１０８へ出力する。
【００７２】
図４は、参照インデックスとパーシスタンス係数との対応関係の一例を示すグラフである。図４において、横軸は参照インデックスを示し、縦軸は、パーシスタンス係数を示している。表３に示したように、参照インデックスは、ＶｍａｘまたはＶｃｕｒｒｅｎｔの絶対値Abs(Ｖｃｕｒｒｅｎｔ)である。Ｖｃｕｒｒｅｎｔの絶対値がしきい値Ｖｔｈ以下である場合には、Ｖｃｕｒｒｅｎｔの増大にともなって単調に減少するパーシスタンス係数Ｃｐｅｒｓｉｓｔｅｎｃｅが対応付けられている。つまり、Ｖｃｕｒｒｅｎｔの絶対値がしきい値Ｖｔｈ以下である場合には、最新のフレームの血流速度データＶｃｕｒｒｅｎｔの大きさに応じて異なるパーシスタンス係数Ｃｐｅｒｓｉｓｔｅｎｃｅが対応付けられる。
【００７３】
上記したように、パーシスタンス演算とは、現在およびそれに連続する複数フレームの画像に対し重み付け加算を行い、その結果を表示に利用することによって、残像効果を与える演算である。パーシスタンス演算の詳細については後述する。
【００７４】
＜パーシスタンス演算部１０８の処理＞
パーシスタンス演算部１０８は、パーシスタンス係数決定部１０６から設定されたパーシスタンス係数と、折り返し判定部１０５からの２つの判定結果と、に基づいて、最新の血流速度データＶｃｕｒｒｅｎｔに対してパーシスタンス演算を行う。
【００７５】
折り返し判定部１０５からの判定結果のうち、Ｖｃｕｒｒｅｎｔの判定結果が真であった場合、すなわち、Ｖｃｕｒｒｅｎｔに折り返しが発生していると判定された場合、パーシスタンス演算部１０８は、パーシスタンス演算を行う前に、Ｖｃｕｒｒｅｎｔに対して以下の補正を行う。具体的には、パーシスタンス演算部１０８は、Ｖｃｕｒｒｅｎｔ＜０の場合はＶｃｕｒｒｅｎｔに繰り返し周波数ＰＲＦを加えてＶｃｕｒｒｅｎｔ＋ＰＲＦに補正し、Ｖｃｕｒｒｅｎｔ＞０の場合はＶｃｕｒｒｅｎｔからＰＲＦを減じてＶｃｕｒｒｅｎｔ−ＰＲＦに補正する。
【００７６】
同様に、折り返し判定部１０５からの判定結果のうち、Ｖｏｕｔ−１の判定結果が真であった場合、すなわち、Ｖｏｕｔ−１に折り返しが発生していると判定された場合、パーシスタンス演算部１０８は、パーシスタンス演算を行う前に、Ｖｏｕｔ−１に対して上記Ｖｃｕｒｒｅｎｔと同様の補正を行う。具体的には、パーシスタンス演算部１０８は、Ｖｏｕｔ−１＜０の場合はＶｏｕｔ−１をＶｏｕｔ−１＋ＰＲＦに補正する。また、パーシスタンス演算部１０８は、Ｖｏｕｔ−１＞０の場合はＶｏｕｔ−１をＶｏｕｔ−１−ＰＲＦに補正する。
【００７７】
そして、パーシスタンス演算部１０８は、血流速度データＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１を用いて、以下の式（１）を用いてパーシスタンス演算を行う。
Ｖｏｕｔ＝（１−Ｃｐｅｒｓｉｓｔｅｎｃｅ）×Ｖｃｕｒｒｅｎｔ＋Ｃｐｅｒｓｉｓｔｅｎｃｅ×Ｖｏｕｔ−１・・（１）
【００７８】
式（１）において、Ｖｏｕｔは最新のフレームにおける、パーシスタンス演算された（残像処理された）血流速度データであり、Ｃｐｅｒｓｉｓｔｅｎｃｅ（０≦Ｃｐｅｒｓｉｓｔｅｎｃｅ≦１）は、パーシスタンス係数決定部１０６により設定されるパーシスタンス係数である。
【００７９】
さらに、パーシスタンス演算部１０８は、上記したように式（１）を用いてパーシスタンス演算を行った結果、Ｖｏｕｔ＞＋ＰＲＦ／２となった場合は、ＶｏｕｔをＶｏｕｔ−ＰＲＦに補正し、Ｖｏｕｔ＜−ＰＲＦ／２となった場合は、ＶｏｕｔをＶｏｕｔ＋ＰＲＦに補正する。そして、パーシスタンス演算部１０８は、補正したＶｏｕｔをＣＦＭ  ＤＳＣ部１１０およびパーシスタンスメモリ部１０９に出力する。
【００８０】
このような処理により、Ｖｃｕｒｒｅｎｔ、Ｖｏｕｔ−１、あるいはＶｏｕｔに折り返しが発生していたとしても、折り返しの影響をキャンセルしたＶｏｕｔを算出することができる。
【００８１】
なお、上記説明したパーシスタンス演算部１０８によるパーシスタンス演算は、１フレーム分の血流速度データ（ＣＦＭフレームデータ）の各画素または測定点ごとに行われる。
【００８２】
このようにパーシスタンス係数決定部１０６が決定したパーシスタンス係数を用いて、パーシスタンス演算部１０８がパーシスタンス演算を行うことにより、以下のような効果が得られる。
【００８３】
すなわち、最新のフレームの血流速度データＶｃｕｒｒｅｎｔの大きさが小さい場合には、パーシスタンス係数決定部１０６が決定するパーシスタンス係数Ｃｐｅｒｓｉｓｔｅｎｃｅは大きくなるため、パーシスタンス演算部１０８が式（１）のようなパーシスタンス演算を行うと、１つ前のフレームの血流速度データＶｏｕｔ−１の重みが大きくなる。その結果、Ｖｃｕｒｒｅｎｔが小さい場合には、Ｖｏｕｔ−１を大きく反映させたパーシスタンス演算の結果であるＶｏｕｔが算出される。このため、Ｖｏｕｔを用いてＣＦＭ  ＤＳＣ部１１０および画像合成部１１３が生成したカラーフローマッピング画像の変化が滑らかとなり、黒抜けが生じにくくなる。
【００８４】
また、Ｖｃｕｒｒｅｎｔが大きい場合には、パーシスタンス係数決定部１０６が決定するパーシスタンス係数Ｃｐｅｒｓｉｓｔｅｎｃｅは小さくなるため、パーシスタンス演算部１０８が式（１）のようなパーシスタンス演算を行うと、Ｖｏｕｔ−１の重みが小さくなる。その結果、Ｖｃｕｒｒｅｎｔが大きい場合には、Ｖｏｕｔ−１の影響が小さいＶｏｕｔが算出され、急激な血流速度の増大をリアルタイムで反映させたカラーフローマッピング画像を実現できる。
【００８５】
また、パーシスタンス係数決定部１０６は、Ｖｃｕｒｒｅｎｔの増大にともなってパーシスタンス係数Ｃｐｅｒｓｉｓｔｅｎｃｅを単調に減少させる。このため、時間の経過とともに血流速度が増大するような場合には、パーシスタンス係数決定部１０６は、時間の経過とともにパーシスタンス係数Ｃｐｅｒｓｉｓｔｅｎｃｅを減少させるので、パーシスタンス演算部１０８によるパーシスタンス演算の残像効果が小さくなり、カラーフローマッピング画像の変化が急激になる。反対に、時間の経過とともに血流速度が減少する場合には、パーシスタンス係数決定部１０６は、時間の経過とともにパーシスタンス係数Ｃｐｅｒｓｉｓｔｅｎｃｅを増大させるので、パーシスタンス演算部１０８によるパーシスタンス演算の残像効果が大きくなり、カラーフローマッピング画像の変化が緩やかになる。
【００８６】
パーシスタンス演算部１０８の処理の説明に戻る。パーシスタンス演算部１０８は、式（１）を用いてパーシスタンス演算を行った後、パーシスタンス演算の結果Ｖｏｕｔを用いて、最新のフレームにおける、残像処理された血流速度データＶｏｕｔに折り返しがあるか否かの判定を行う。ここで、パーシスタンス演算部１０８は、Ｖｏｕｔ＞＋ＰＲＦ／２またはＶｏｕｔ＜−ＰＲＦ／２となった場合はＶｏｕｔに折り返し有りと判定し、それ以外（−ＰＲＦ／２≦Ｖｏｕｔ≦ＰＲＦ／２）の場合は折り返し無しと判定する。
【００８７】
そして、パーシスタンス演算部１０８は、最新のフレームにおける血流速度データＶｏｕｔの折り返し判定結果Ｒａｌｉａｓｉｎｇをパーシスタンスメモリ部１０９に出力する。パーシスタンス演算部１０８から出力された折り返し判定結果Ｒａｌｉａｓｉｎｇは、パーシスタンスメモリ部１０９に記憶され、次のフレームにおける処理の際に、最新より１つ前のフレームの折り返し判定結果Ｒａｌｉａｓｉｎｇ−１として折り返し判定部１０５に出力される。
【００８８】
＜作用・効果＞
以上説明したように、本発明の第１の実施の形態に係る超音波診断装置１１では、折り返し判定部１０５が、最新のフレームの血流速度データＶｃｕｒｒｅｎｔと、最新より１つ前のフレームの血流速度データＶｏｕｔ−１と、に加えて、最新より以前のフレームにおいて判定された、当該フレームにおいて残像処理された血流速度データの折り返し判定結果Ｒａｌｉａｓｉｎｇ−１に基づいて折り返し判定を行う。このため、最新のフレームの血流速度データＶｃｕｒｒｅｎｔと、最新より１つ前のフレームの血流速度データＶｏｕｔ−１と、に折り返しが発生したか否かを正確に判定することができる。
【００８９】
また、本発明の第１の実施の形態に係る超音波診断装置１１では、このような折り返し判定方法による折り返し判定の結果に基づいて、パーシスタンス演算部１０８が、ＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１を用いてパーシスタンス演算を行い、その結果を出力する。このため、血流速度を正確に反映した断層画像を生成して表示することができる。
【００９０】
さらに、本発明の第１の実施の形態に係る超音波診断装置１１では、パーシスタンス演算の結果である、最新のフレームにおいて残像処理された血流速度データＶｏｕｔに基づいて、パーシスタンス演算部１０８がＶｏｕｔに折り返しが発生しているか否かの判定を行い、その結果であるＲａｌｉａｓｉｎｇをパーシスタンスメモリ部１０９に出力する。そして、次のフレームの処理において、パーシスタンスメモリ部１０９は、前のフレームの処理において入力されたＶｏｕｔを次のフレームにおけるＶｏｕｔ−１として、前のフレームの処理において入力されたＲａｌｉａｓｉｎｇをＲａｌｉａｓｉｎｇ−１として、それぞれ出力する。
【００９１】
そして、本発明の第１の実施の形態に係る超音波診断装置１１では、最新のフレームの血流速度データＶｃｕｒｒｅｎｔあるいは最新より１つ前のフレームの血流速度データＶｏｕｔ−１に折り返しが発生していると判定された場合、パーシスタンス演算部１０８は、パーシスタンス演算を行う前に、折り返しが発生した血流速度データに対して補正を行う。さらに、パーシスタンス演算部１０８は、パーシスタンス演算を行った結果である、最新のフレームにおいて残像処理された血流速度データＶｏｕｔに折り返しが発生しているか否かを判定し、折り返しが発生していると判定した場合は、Ｖｏｕｔに対して補正を行う。このため、Ｖｃｕｒｒｅｎｔ、Ｖｏｕｔ−１、あるいはＶｏｕｔに折り返しが発生していたとしても、折り返しの影響をキャンセルしたＶｏｕｔを算出することができる。
【００９２】
（第２の実施の形態）
以下、本発明の第２の実施の形態に係る超音波診断装置１２について説明する。図５は、本発明の第２の実施の形態に係る超音波診断装置の構成の一例を示すブロック図である。
【００９３】
図５に示す超音波診断装置１２は、残像処理部がパーシスタンス係数決定部、パーシスタンス係数参照メモリ部、パーシスタンス演算部を２つずつ有する点において上記第１の実施の形態の超音波診断装置１１と異なっている。以下では、第１の実施の形態との相違点についてのみ説明する。
【００９４】
図５に示すように、残像処理部１１５’は、フレームメモリ部１０４と、折り返し判定部１０５と、第１パーシスタンス係数決定部１０６Ａと、第２パーシスタンス係数決定部１０６Ｂと、第１パーシスタンス係数参照メモリ部１０７Ａと、第２パーシスタンス係数参照メモリ部１０７Ｂと、第１パーシスタンス演算部１０８Ａと、第２パーシスタンス演算部１０８Ｂと、パーシスタンスメモリ部１０９と、を有する。
【００９５】
第１パーシスタンス係数決定部１０６Ａは、折り返し判定部１０５から出力された２つの折り返し判定結果と、フレームメモリ部１０４から読み出したＶｃｕｒｒｅｎｔの絶対値と、に基づいて、第１パーシスタンス係数参照メモリ部１０７Ａへの第１参照インデックスを作成する。第１参照インデックスの作成については後述する。
【００９６】
また、第１パーシスタンス係数決定部１０６Ａは、第１パーシスタンス係数参照メモリ部１０７Ａにアクセスし、第１参照インデックスに対応付けられた第１パーシスタンス係数を読み出し、第１パーシスタンス演算部１０８Ａに設定する。第１パーシスタンス係数参照メモリ部１０７Ａには、あらかじめ血流速度データの大きさに対応付けられた第１パーシスタンス係数を含む第１参照テーブルが記憶されている。この第１参照テーブルは、血流速度データの大きさに対応付けられた異なる２つ以上の値のパーシスタンス係数を含む。
【００９７】
これに対し、第２パーシスタンス係数決定部１０６Ｂは、折り返し判定部１０５から出力された２つの折り返し判定結果と、パーシスタンスメモリ部１０９から読み出したＶｏｕｔ−１の絶対値と、に基づいて、第２パーシスタンス係数参照メモリ部１０７Ｂへの第２参照インデックスを作成する。第２パーシスタンス係数決定部１０６Ｂにおける第２参照インデックスの作成については後述する。
【００９８】
また、第２パーシスタンス係数決定部１０６Ｂは、第２パーシスタンス係数参照メモリ部１０７Ｂにアクセスし、第２参照インデックスに対応付けられた第２パーシスタンス係数を読み出し、第２パーシスタンス演算部１０８Ｂに設定する。第２パーシスタンス係数参照メモリ部１０７Ｂには、あらかじめ血流速度データの大きさに対応付けられた第２パーシスタンス係数を含む第２参照テーブルが記憶されている。第２参照テーブルも、血流速度の値に応じて対応付けられた異なる２つ以上の値のパーシスタンス係数を含むが、以下において詳細に説明するように、同じ血流速度値に対して対応付けられる第１パーシスタンス係数と第２パーシスタンス係数とは値が異なっている。
【００９９】
第１パーシスタンス演算部１０８Ａは、第１パーシスタンス係数決定部１０６Ａから設定されたパーシスタンス係数と、折り返し判定部１０５からの２つの判定結果と、に基づいて、第１の実施の形態のパーシスタンス演算部１０８と同様に、折り返しの補正を行ってから、式（１）を用いてＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１に対してパーシスタンス演算を行う。
【０１００】
同様に、第２パーシスタンス演算部１０８Ｂも、第２パーシスタンス係数決定部１０６Ｂから設定されたパーシスタンス係数と折り返し判定部１０５からの２つの判定結果と、に基づいて、第１の実施の形態のパーシスタンス演算部１０８と同様に、折り返しの補正を行ってから、式（１）を用いてＶｃｕｒｒｅｎｔおよびＶｏｕｔ−１に対してパーシスタンス演算を行う。
【０１０１】
第１パーシスタンス演算部１０８Ａおよび第２パーシスタンス演算部１０８Ｂにおける演算は、設定されたパーシスタンス係数が異なる点を除いて同じである。
【０１０２】
第１パーシスタンス演算部１０８Ａおよび第２パーシスタンス演算部１０８Ｂは、パーシスタンス演算の結果をそれぞれ最大値選択部１１６に出力する。
【０１０３】
最大値選択部１１６は、第１パーシスタンス演算部１０８Ａおよび第２パーシスタンス演算部１０８Ｂからそれぞれ入力されたパーシスタンス演算の結果に基づいて、画素ごとに、あるいは、対応する測定点のデータごとにパーシスタンス演算の結果である血流速度データの絶対値を比較し、大きい方を選択して、最新のフレームにおいて残像処理された血流速度データＶｏｕｔを生成する。
【０１０４】
また、最大値選択部１１６は、第１の実施の形態のパーシスタンス演算部１０８が行ったように、Ｖｏｕｔに対して折り返しの補正を行う。
【０１０５】
そして、最大値選択部１１６は、補正を行った補正したＶｏｕｔをＣＦＭ  ＤＳＣ部１１０およびパーシスタンスメモリ部１０９に出力する。
【０１０６】
さらに、最大値選択部１１６は、第１の実施の形態のパーシスタンス演算部１０８が行ったように、最新のフレームにおける、残像処理された血流速度データＶｏｕｔに折り返しがあるか否かの判定を行う。そして、最大値選択部１１６は、最新のフレームにおける血流速度データＶｏｕｔの折り返し判定結果Ｒａｌｉａｓｉｎｇをパーシスタンスメモリ部１０９に出力する。
【０１０７】
次に、第１および第２パーシスタンス係数の決定方法について詳細に説明する。上述したように、第１パーシスタンス係数決定部１０６Ａは、折り返し判定部１０５から出力された２つの折り返し判定結果と、フレームメモリ部１０４から読み出したＶｃｕｒｒｅｎｔの絶対値と、に基づいて、第１パーシスタンス係数参照メモリ部１０７Ａへの第１参照インデックスを作成する。
【０１０８】
表４に作成される第１参照インデックスを示す。
【表４】
2018050925.tif 000006
【０１０９】
Ｖｃｕｒｒｅｎｔに折り返しが発生している場合、Ｖｃｕｒｒｅｎｔは、実際には、Ｖｍａｘまたは−Ｖｍａｘを超える絶対値の大きな値であると考えられる。このため、Ｖｃｕｒｒｅｎｔに折り返しが発生している場合の第１参照インデックスはＶｍａｘとなる。その他の場合には、Ｖｃｕｒｒｅｎｔの絶対値Abs(Ｖｃｕｒｒｅｎｔ)となる。
【０１１０】
図６Ａは、第１参照インデックスと第１パーシスタンス係数との対応関係の一例を示すグラフである。図６Ａにおいて、横軸は第１参照インデックスを示し、縦軸は、第１パーシスタンス係数を示している。表４に示すように、第１参照インデックスは、ＶｍａｘまたはＶｃｕｒｒｅｎｔの絶対値Abs(Ｖｃｕｒｒｅｎｔ)である。Ｖｃｕｒｒｅｎｔの絶対値がしきい値Ｖｔｈ以下である場合には、Ｖｃｕｒｒｅｎｔの増大にともなって単調に減少する第１パーシスタンス係数Ｃｐｅｒｓｉｓｔｅｎｃｅが対応付けられている。つまり、Ｖｃｕｒｒｅｎｔの絶対値がしきい値Ｖｔｈ以下である場合には、最新のフレームの血流速度データＶｃｕｒｒｅｎｔの大きさに応じて異なるパーシスタンス係数Ｃｐｅｒｓｉｓｔｅｎｃｅが対応付けられる。
【０１１１】
一方、第２パーシスタンス係数決定部１０６Ｂは、折り返し判定部１０５から出力された２つの折り返し判定結果と、パーシスタンスメモリ部１０９から読み出したＶｏｕｔ−１の絶対値と、に基づいて、第２パーシスタンス係数参照メモリ部１０７Ｂへの第２参照インデックスを作成する。第２参照インデックスについても、表４に示す。
【０１１２】
Ｖｏｕｔ−１に折り返しが発生している場合、Ｖｏｕｔ−１は、実際には、Ｖｍａｘまたは−Ｖｍａｘを超える絶対値の大きな値であると考えられる。このため、Ｖｏｕｔ−１に折り返しが発生している場合の第２参照インデックスはＶｍａｘとなる。その他の場合には、Ｖｏｕｔ−１の絶対値Abs(Ｖｏｕｔ−１)となる。
【０１１３】
図６Ｂは、第２参照インデックスと第２パーシスタンス係数との対応関係の一例を示すグラフである。図６Ｂにおいて、横軸は第２参照インデックスを示し、縦軸は、第２パーシスタンス係数を示している。表４に示すように、第２参照インデックスは、ＶｍａｘまたはＶｏｕｔ−１の絶対値Abs(Ｖｏｕｔ−１)である。Ｖｏｕｔ−１の絶対値がしきい値Ｖｔｈ以下である場合には、Ｖｏｕｔ−１の増大にともなって単調に増加する第２パーシスタンス係数Ｃｐｅｒｓｉｓｔｅｎｃｅが対応付けられている。つまり、Ｖｏｕｔ−１の絶対値がしきい値Ｖｔｈ以下である場合には、最新より１つ前のフレームの血流速度データＶｏｕｔ−１に応じて異なるパーシスタンス係数Ｃｐｅｒｓｉｓｔｅｎｃｅが対応付けられる。
【０１１４】
図６Ａおよび図６Ｂに示すように、第１パーシスタンス係数は、最新のフレームの血流速度Ｖｃｕｒｒｅｎｔに対応付けられ、かつ、Ｖｃｕｒｒｅｎｔの絶対値が大きいほど小さな値である。第１パーシスタンス係数が大きくなれば、１つ前のフレームの血流速度をより考慮した演算となるため、第１パーシスタンス演算部１０８Ａは、残像効果を抑制し、すみやかに血流速度を変化させる演算を行う。これに対し、第２パーシスタンス係数は、１つ前のフレームの血流速度Ｖｏｕｔ−１に対応付けられ、かつ、Ｖｏｕｔ−１の絶対値が大きいほど大きな値であるため、第２パーシスタンス演算部１０８Ｂは、残像効果を高め、血流速度の変化を抑制する演算を行う。
【０１１５】
また、上述したように第１パーシスタンス演算部１０８Ａは、残像効果を抑制し、すみやかに血流速度を変化させる演算を行うため、第１参照インデックスが増加するにつれて第１パーシスタンス係数を単調に減少させ、血流速度が高くなるにつれて急激な血流速度の増大をリアルタイムで反映させたカラーフローマッピング画像を実現できる。
【０１１６】
また、第２パーシスタンス演算部１０８Ｂは、残像効果を高めた画像表示を行うため、血流速度が高くなるにつれて、血流画像に色調あるいは諧調の色付けを行っていると、血流速度が低い場合には、暗めの表示が必要以上に長い時間残像として表示されてしまう。例えば、探触子を動かした場合に、血流表示が尾引く印象を与えてしまう。このため、第２参照インデックスが増加するにつれて第２パーシスタンス係数を単調に増加させ、血流速度が低くなるにつれて残像効果を抑制することができる。
【０１１７】
（第３の実施の形態）
以下、本発明の第３の実施の形態に係る超音波診断装置１３について説明する。図７は、本発明の第３実施の形態に係る超音波診断装置の構成の一例を示すブロック図である。
【０１１８】
図７に示す超音波診断装置１３は、残像処理部が２つのパーシスタンス係数決定部、２つのパーシスタンス係数参照メモリ部を２つずつ有し、２つのパーシスタンス係数決定部が決定したパーシスタンス係数のうちいずれかを選択するパーシスタンス係数選択部を有する点において上記第１の実施の形態の超音波診断装置１１と異なっている。以下では、第１および第２の実施の形態との相違点についてのみ説明する。
【０１１９】
図７に示すように、残像処理部１１５’’は、フレームメモリ部１０４と、折り返し判定部１０５と、第１パーシスタンス係数決定部１０６Ａと、第２パーシスタンス係数決定部１０６Ｂと、パーシスタンス係数選択部１１７と、第１パーシスタンス係数参照メモリ部１０７Ａと、第２パーシスタンス係数参照メモリ部１０７Ｂと、パーシスタンス演算部１０８と、パーシスタンスメモリ部１０９と、を有する。
【０１２０】
パーシスタンス係数選択部１１７は、フレームメモリ部１０４から最新のフレームの血流速度データＶｃｕｒｒｅｎｔを、パーシスタンスメモリ部１０９から最新より１つ前のフレームの血流速度データＶｏｕｔ−１を、それぞれ読み出し、Ｖｃｕｒｒｅｎｔの絶対値とＶｏｕｔ−１の絶対値とを比較して、比較結果に基づいて、第１パーシスタンス係数決定部１０６Ａが決定した第１パーシスタンス係数と、第２パーシスタンス係数決定部１０６Ｂが決定した第２パーシスタンス係数のいずれをパーシスタンス演算部１０８に設定するかを選択する。
【０１２１】
具体的には、パーシスタンス係数選択部１１７は、Ｖｃｕｒｒｅｎｔの絶対値の方がＶｏｕｔ−１の絶対値より大きい場合は、第１パーシスタンス係数を選択し、そうでない場合は第２パーシスタンス係数を選択する。
【０１２２】
従って、第２の実施の形態と同様に、第１パーシスタンス係数は、最新のフレームの血流速度Ｖｃｕｒｒｅｎｔに対応付けられ、かつ、Ｖｃｕｒｒｅｎｔの絶対値が大きいほど小さな値となる。第１パーシスタンス係数が大きくなれば、１つ前のフレームの血流速度をより考慮した演算となるため、第１パーシスタンス係数を用いた場合のパーシスタンス演算部１０８は、残像効果を抑制し、すみやかに血流速度を変化させる演算を行うことができる。これに対し、第２パーシスタンス係数は、１つ前のフレームの血流速度Ｖｏｕｔ−１に対応付けられ、かつ、Ｖｏｕｔ−１の絶対値が大きいほど大きな値であるため、第２パーシスタンス係数を用いた場合のパーシスタンス演算部１０８は、残像効果を高め、血流速度の変化を抑制する演算を行うことができる。
【０１２３】
以上、図面を参照しながら各種実施の形態について説明したが、本発明はかかる例に限定されない。特許請求の範囲の記載範囲内において、当業者が想到できる各種の変更例または修正例についても、本発明の技術的範囲に含まれる。また、開示の趣旨を逸脱しない範囲において、上記実施の形態における各構成要素を任意に組み合わせてもよい。
【０１２４】
上記実施の形態では、ＣＦＭフレームデータに含まれる血流速度データを用いて処理を行っていたが、ＣＦＭフレームデータに血流速度データＶｃｕｒｒｅｎｔあるいはＶｏｕｔ−１以外のデータ（血流パワーデータや血流速度の分散データ等）が含まれる場合には、パーシスタンス演算部１０８は、血流速度データ以外のデータに対しても、血流速度データと同様に最新のフレームのデータおよび最新より１つ前のフレームのデータと求めたパーシスタンス係数Ｃｐｅｒｓｉｓｔｅｎｃｅを用いて上記式（１）と同様のパーシスタンス演算を行い、残像処理されたデータを得るようにしてもよい。
【０１２５】
また、上記実施の形態では、最新のフレームおよび１つ前のフレームの血流速度データを用いてパーシスタンス処理を行っていたが、２つ前、あるいは３つ以上前のフレームの血流速度データも用いてパーシスタンス処理を行ってもよい。また、式（１）に限られず、他の演算式を用いてパーシスタンス処理を行ってもよい。
【０１２６】
上記実施の形態の説明に用いた各機能ブロックは、典型的には、入力端子および出力端子を有する集積回路であるＬＳＩとして実現される。これらは個別に１チップ化されてもよいし、一部または全てを含むように１チップ化されてもよい。ここでは、ＬＳＩとしたが、集積度の違いにより、ＩＣ（Integrated Circuit）、システムＬＳＩ、スーパーＬＳＩ、ウルトラＬＳＩと呼称されることもある。
【０１２７】
また、集積回路化の手法はＬＳＩに限るものではなく、専用回路または汎用プロセッサを用いて実現してもよい。ＬＳＩ製造後に、プログラムすることが可能なＦＰＧＡ（Field Programmable Gate Array）、ＬＳＩ内部の回路セルの接続または設定を再構成可能なリコンフィギュラブル  プロセッサ（Reconfigurable Processor）を利用してもよい。
【０１２８】
さらには、半導体技術の進歩または派生する別技術により、ＬＳＩに置き換わる集積回路化の技術が登場すれば、当然、その技術を用いて機能ブロックを集積化してもよい。
【産業上の利用可能性】
【０１２９】
本発明は、血流速度を表示することができる超音波診断装置に好適である。
【符号の説明】
【０１３０】
１１，１２，１３  超音波診断装置
１０１  探触子
１０２  超音波送受信部
１０３  ＣＦＭ信号処理部
１０４  フレームメモリ部
１０５  折り返し判定部
１０６  パーシスタンス係数決定部
１０６Ａ  第１パーシスタンス係数決定部
１０６Ｂ  第２パーシスタンス係数決定部
１０７  パーシスタンス係数参照メモリ部
１０７Ａ  第１パーシスタンス係数参照メモリ部
１０７Ｂ  第２パーシスタンス係数参照メモリ部
１０８  パーシスタンス演算部
１０８Ａ  第１パーシスタンス演算部
１０８Ｂ  第２パーシスタンス演算部
１０９  パーシスタンスメモリ部
１１０  ＣＦＭ  ＤＳＣ部
１１１  断層画像信号処理部
１１２  断層画像ＤＳＣ部
１１３  画像合成部
１１４  表示部
１１５，１１５’，１１５’’  残像処理部
１１６  最大値選択部
１１７  パーシスタンス係数選択部
４０１  探触子
４０２  超音波送受信部
４０３  カラーフローマッピング処理部
４０３  ＣＦＭ信号処理部
４０４  フレームメモリ部
４０５  フレームメモリ選択部
４０６  パーシスタンス係数設定部
４０７  パーシスタンス演算部
４０８  ＣＦＭ  ＤＳＣ部
４０９  断層画像信号処理部
４１０  断層画像ＤＳＣ部
４１１  画像合成部
４１２  表示部
【図１】
2018050925.tif 000007
【図２】
2018050925.tif 000008
【図３Ａ】
2018050925.tif 000009
【図３Ｂ】
2018050925.tif 000010
【図３Ｃ】
2018050925.tif 000011
【図３Ｄ】
2018050925.tif 000012
【図３Ｅ】
2018050925.tif 000013
【図３Ｆ】
2018050925.tif 000014
【図４】
2018050925.tif 000015
【図５】
2018050925.tif 000016
【図６Ａ】
2018050925.tif 000017
【図６Ｂ】
2018050925.tif 000018
【図７】
2018050925.tif 000019
【図８】
2018050925.tif 000020