(57)【要約】【課題】従来よりも生体組織の弾性を正確に反映し、画質が向上した弾性画像を得ることができる超音波診断装置又はその制御プログラムを提供する。【解決手段】時間的に異なる二つのフレームのいずれか一方に属する一の音線上のエコー信号への相関ウィンドウの設定にあたり、この相関ウィンドウと前記一の音線上において隣り合っていて直前に相関演算を行なった直前相関ウィンドウの相関演算で得られた算出値と、前記一の音線とは異なる他の音線上における他音線相関ウィンドウの相関演算で得られた算出値との平均演算を行なう平均演算部５３を備え、この平均演算部５３で得られた平均値に基づいて前記相関ウィンドウの設定を行なうことを特徴とする。【選択図】図２
2011024654.tif 000002
【特許請求の範囲】
【請求項１】
生体組織に対し超音波の送受信を行なって得られたエコー信号であって、時間的に異なる二つのフレームに属する同一音線上における二つのエコー信号に相関ウィンドウを設定し、該相関ウィンドウ間で相関演算を行なって生体組織の弾性に関する物理量を算出する物理量算出部と、
前記物理量に基づいて生体組織の弾性画像データを作成する弾性画像データ作成部と、
を備え、
前記物理量算出部は、前記二つのフレームのいずれか一方に属する一の音線上のエコー信号への相関ウィンドウの設定にあたり、該相関ウィンドウと前記一の音線上において隣り合っていて直前に相関演算を行なった直前相関ウィンドウの相関演算と、前記一の音線とは異なる他の音線上における他音線相関ウィンドウの相関演算とに従って、前記相関ウィンドウの設定を行なう
ことを特徴とする超音波診断装置。
【請求項２】
前記一の音線上における相関ウィンドウの設定にあたり、前記直前相関ウィンドウの相関演算で得られた算出値と、前記他音線相関ウィンドウの相関演算で得られた算出値との平均演算を行なう平均演算部を備え、
前記物理量算出部は、前記平均演算部で得られた平均値に基づいて前記相関ウィンドウの設定を行なうことを特徴とする請求項１に記載の超音波診断装置。
【請求項３】
前記平均演算部は、複数の前記他の音線についての前記他音線相関ウィンドウの相関演算の算出値と、前記直前相関ウィンドウの相関演算で得られた算出値との平均演算を行なうことを特徴とする請求項２に記載の超音波診断装置。
【請求項４】
前記他音線相関ウィンドウの相関演算で得られた算出値について、エラーか否の判定を行なうエラー判定部を備え、
前記平均演算部は、エラーの算出値を除いて前記平均演算を行なう
ことを特徴とする請求項３に記載の超音波診断装置。
【請求項５】
前記直前相関ウィンドウの相関演算で得られた算出値について、エラーか否かの判定を行なうエラー判定部と、
前記一の音線上における相関ウィンドウの設定にあたり、前記エラー判定部により、前記直前相関ウィンドウの相関演算で得られた算出値についてエラーと判定された場合、複数の前記他の音線についての前記他音線相関ウィンドウの相関演算で得られた算出値の平均演算を行なう平均演算部と、
をさらに備え、
前記物理量算出部は、前記平均演算部で得られた平均値に基づいて前記相関ウィンドウの設定を行なう
ことを特徴とする請求項１に記載の超音波診断装置。
【請求項６】
前記エラー判定部は、判定対象の算出値が得られた相関演算における相関係数に基づいてエラーか否かを判定することを特徴とする請求項４又は５に記載の超音波診断装置。
【請求項７】
前記エラー判定部は、判定対象の算出値が予め設定された所定の範囲内にない場合にエラーと判定することを特徴とする請求項４又は５に記載の超音波診断装置。
【請求項８】
前記エラー判定部は、判定対象の算出値について、他の判定対象の算出値の分布に基づいて、エラーか否かを判定することを特徴とする請求項４又は５に記載の超音波診断装置。
【請求項９】
前記平均演算部は、前記平均演算の対象となる算出値に対し、相関係数に応じた重み付けを行なうことを特徴とする請求項２〜８のいずれか一項に記載の超音波診断装置。
【請求項１０】
前記一の音線上における相関ウィンドウの設定にあたり、前記直前相関ウィンドウの相関演算で得られた算出値と、前記他音線相関ウィンドウの相関演算で得られた算出値の中から、前記相関ウィンドウの設定に適した算出値を選択する選択部を備え、
前記物理量算出部は、前記選択部によって選択された算出値に基づいて前記相関ウィンドウの設定を行なうことを特徴とする請求項１に記載の超音波診断装置。
【請求項１１】
前記選択部は、前記相関ウィンドウの設定に適した算出値の選択を、該算出値が得られた相関演算の相関係数又は前記算出値が所定の範囲内にあるか否かに基づいて行なうことを特徴とする請求項１０に記載の超音波診断装置。
【請求項１２】
前記選択部は、前記直前相関ウィンドウの相関演算における相関係数が所定の閾値を超える場合、該直前相関ウィンドウの相関演算によって得られた算出値を前記相関ウィンドウの設定に適した算出値として選択し、一方で前記直前相関ウィンドウの相関演算における相関係数が前記所定の閾値以下である場合、前記他音線相関ウィンドウの相関演算で得られた算出値であって、前記直前相関ウィンドウの相関演算における相関係数よりも高い相関係数の相関演算で得られた算出値又は所定の範囲内の算出値を、前記相関ウィンドウの設定に適した算出値として選択することを特徴とする請求項１０又は１１に記載の超音波診断装置。
【請求項１３】
前記選択部は、前記直前相関ウィンドウの相関演算における相関係数及び前記他音線相関ウィンドウの相関演算における相関係数のうち、相関係数が最も高い相関演算によって得られた算出値を、前記相関ウィンドウの設定に適した算出値として選択することを特徴とする請求項１０又は１１に記載の超音波診断装置。
【請求項１４】
前記選択部は、前記直前相関ウィンドウの相関演算によって得られた算出値が所定の範囲内である場合、該直前相関ウィンドウの相関演算によって得られた算出値を前記相関ウィンドウの設定に適した算出値として選択し、一方で前記直前相関ウィンドウの相関演算によって得られた算出値が所定の範囲外である場合、前記他音線相関ウィンドウの相関演算で得られた算出値であって、前記直前相関ウィンドウの相関演算における相関係数よりも高い相関係数の相関演算で得られた算出値又は前記所定の範囲内の算出値を、前記相関ウィンドウの設定に適した算出値として選択することを特徴とする請求項１０又は１１に記載の超音波診断装置。
【請求項１５】
前記選択部は、複数の前記他の音線についての前記他音線相関ウィンドウの相関演算で得られた算出値の分布に基づいて、前記直前相関ウィンドウの相関演算によって得られた算出値が前記相関ウィンドウの設定に適するか否かを判断し、該相関ウィンドウの設定に適していると判断した場合には前記直前相関ウィンドウの相関演算によって得られた算出値を前記相関ウィンドウの設定に適した算出値として選択し、一方で、前記相関ウィンドウの設定に適していないと判断した場合には、いずれかの前記他音線相関ウィンドウの相関演算で得られた算出値であって、前記直前相関ウィンドウの相関演算における相関係数よりも高い相関係数の相関演算で得られた算出値又は前記所定の範囲内の算出値を、前記相関ウィンドウの設定に適した算出値として選択することを特徴とする請求項１０又は１１に記載の超音波診断装置。
【請求項１６】
前記他音線相関ウィンドウは、前記直前相関ウィンドウと生体組織において同じ深さに位置することを特徴とする請求項１〜１５のいずれか一項に記載の超音波診断装置。
【請求項１７】
コンピュータに、
生体組織に対し超音波の送受信を行なって得られたエコー信号であって、時間的に異なる二つのフレームに属する同一音線上における二つのエコー信号に、相関ウィンドウを設定し、該相関ウィンドウ間で相関演算を行なって生体組織の弾性に関する物理量を算出する物理量算出機能と、
前記物理量に基づいて生体組織の弾性画像データを作成する弾性画像データ作成機能と、を実行させ、
前記物理量算出機能にあっては、前記二つのフレームのいずれか一方に属する一の音線上のエコー信号への相関ウィンドウの設定にあたり、該相関ウィンドウと前記一の音線上において隣り合っていて直前に相関演算を行なった直前相関ウィンドウの相関演算と、前記一の音線とは異なる他の音線上における他音線相関ウィンドウの相関演算とに従って、前記相関ウィンドウの設定を行なう
ことを特徴とする超音波診断装置の制御プログラム。
【請求項１８】
前記一の音線上における相関ウィンドウの設定にあたり、前記直前相関ウィンドウの相関演算で得られた算出値と、前記他音線相関ウィンドウの相関演算で得られた算出値との平均演算を行なう平均演算機能をコンピュータにさらに実行させ、
前記物理量算出機能にあっては、前記平均演算機能で得られた平均値に基づいて前記相関ウィンドウの設定を行なう
ことを特徴とする請求項１７に記載の超音波診断装置の制御プログラム。
【請求項１９】
前記直前相関ウィンドウの相関演算で得られた算出値について、エラーか否かの判定を行なうエラー判定機能と、
前記一の音線上における相関ウィンドウの設定にあたり、前記エラー判定機能により、前記直前相関ウィンドウの相関演算で得られた算出値についてエラーと判定された場合、複数の前記他の音線についての前記他音線相関ウィンドウの相関演算で得られた算出値の平均演算を行なう平均演算機能と、
をコンピュータにさらに実行させ、
前記物理量算出機能にあっては、前記平均演算機能で得られた平均値に基づいて前記相関ウィンドウの設定を行なう
ことを特徴とする請求項１７に記載の超音波診断装置の制御プログラム。
【請求項２０】
前記一の音線上における相関ウィンドウの設定にあたり、前記直前相関ウィンドウの相関演算で得られた算出値と、前記他音線相関ウィンドウの相関演算で得られた算出値の中から、前記相関ウィンドウの設定に適した算出値を選択する選択機能をコンピュータにさらに実行させ、
前記物理量算出機能は、前記選択機能によって選択された算出値に基づいて前記相関ウィンドウの設定を行なう
ことを特徴とする請求項１７に記載の超音波診断装置の制御プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波診断装置に関し、特に生体組織の硬さ又は軟らかさを表す弾性画像を表示する超音波診断装置及びその制御プログラムに関する。
【背景技術】
【０００２】
通常のＢモード画像と、生体組織の硬さ又は軟らかさを表す弾性画像とを合成して表示させる超音波診断装置が、例えば特許文献１などに開示されている。この種の超音波診断装置において、弾性画像は次のようにして作成される。先ず、生体組織に対し、超音波プローブによって体表面からの圧迫とその弛緩を繰り返しながら超音波の送受信を行い、エコー信号を取得する。そして、得られたエコー信号に基づいて、生体組織の弾性に関する物理量を算出し、この物理量を色相情報に変換してカラーの弾性画像を作成する。ちなみに、生体組織の弾性に関する物理量としては、例えば生体組織の変形による変位（以下、単に「変位」と云う）などを算出している。
【０００３】
前記物理量の算出手法の一例についてもう少し説明すると、先ず時間的に異なる二つのフレームに属する同一音線上の二つのエコー信号に、所定のデータ数分の幅を有する相関ウィンドウをそれぞれ設定し、この相関ウィンドウ間で相関演算を行なって前記物理量を算出する。例えば特許文献２では、相関ウィンドウ間で相関演算を行なうことによって、両エコー信号の波形のずれを算出し、この波形のずれを変位とみなしている。
【０００４】
前記相関ウィンドウは音線方向に順次設定され、各相関ウィンドウ毎に相関演算が行なわれて前記物理量が算出される。ここで、前記特許文献２では、異なるフレームに属する二つのエコー信号のうち、一のフレームに属するエコー信号については、同一音線上において直前に相関演算を行なった隣り合う直前相関ウィンドウから一定量だけ移動させて前記相関ウィンドウの設定を行なう。一方で、他のフレームに属するエコー信号については、前記直前相関ウィンドウからの移動量を、この直前相関ウィンドウを対象とする相関演算によって得られた前記物理量を用いて決定し、前記相関ウィンドウを設定する（特許文献２の［００４４］段落参照）。このようにして二つのフレームに属するエコー信号に設定された相関ウィンドウ間のマッチング度合いが高いほど、相関演算における相関係数が高くなり、得られる算出値が、生体組織の弾性をより正確に反映した算出値となる。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００５−１１８１５２号公報
【特許文献２】特開２００８−１２６０７９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
ところで、エコー信号の質が悪い場合には、相関ウィンドウ間のマッチングの度合いが低くなって相関係数が低い相関演算となり、生体組織の弾性を正確に反映した物理量を得ることができない。例えば、生体組織の中に石灰化などによって局所的に硬い部分がある場合、この硬い部分が前記超音波プローブによる圧迫とその弛緩によって横方向にずれることがある。この場合、ずれた部分のエコー信号の信号波形は、二つのフレームのそれぞれにおいて異なる波形になるため、相関演算を行なう相関ウィンドウ間のマッチングの度合いが低くなり、相関係数が低くなる。また、信号強度が小さい部分に設定された相関ウィンドウについて相関演算を行なった場合も、上記と同様にマッチングの度合いが低くなり相関係数が低くなる。このように、エコー信号の質が悪い場合、相関演算における相関係数が低くなり、算出される物理量が生体組織の弾性を正確に反映したものにならない。
【０００７】
ここで、前記直前相関ウィンドウのマッチング度合いが低く、その相関演算によって得られた物理量が、生体組織の弾性を正確に反映したものではない場合、この物理量に基づいて設定される次の相関ウィンドウについての相関演算の相関係数も低くなり、その演算結果として得られる物理量も生体組織の弾性を正確に反映したものとはならなくなる。従って、ある音線上において、生体組織の弾性を正確に反映した物理量が得られない相関ウィンドウが存在すると、信号波形によってはこの相関ウィンドウ以降に設定される相関ウィンドウ間のマッチング度合いが低い状態が続き、弾性画像に線状のアーチファクトが現れることがある。
【０００８】
本発明が解決しようとする課題は、従来よりも生体組織の弾性を正確に反映した弾性画像を得ることができる超音波診断装置又はその制御プログラムを提供することである。
【課題を解決するための手段】
【０００９】
この発明は、前記課題を解決するためになされたもので、第１の観点の発明は、生体組織に対し超音波の送受信を行なって得られたエコー信号であって、時間的に異なる二つのフレームに属する同一音線上における二つのエコー信号に、相関ウィンドウを設定し、該相関ウィンドウ間で相関演算を行なって生体組織の弾性に関する物理量を算出する物理量算出部と、前記物理量に基づいて生体組織の弾性画像データを作成する弾性画像データ作成部と、を備え、前記物理量算出部は、前記二つのフレームのいずれか一方に属する一の音線上のエコー信号への相関ウィンドウの設定にあたり、該相関ウィンドウと前記一の音線上において隣り合っていて直前に相関演算を行なった直前相関ウィンドウの相関演算と、前記一の音線とは異なる他の音線上における他音線相関ウィンドウの相関演算とに従って、前記相関ウィンドウの設定を行なうことを特徴とする超音波診断装置である。
【００１０】
第２の観点の発明は、第１の観点の発明において、前記一の音線上における相関ウィンドウの設定にあたり、前記直前相関ウィンドウの相関演算で得られた算出値と、前記他音線相関ウィンドウの相関演算で得られた算出値との平均演算を行なう平均演算部を備え、前記物理量算出部は、前記平均演算部で得られた平均値に基づいて前記相関ウィンドウの設定を行なうことを特徴とする超音波診断装置である。
【００１１】
第３の観点の発明は、第２の観点の発明において、前記平均演算部は、複数の前記他の音線についての前記他音線相関ウィンドウの相関演算の算出値と、前記直前相関ウィンドウの相関演算で得られた算出値との平均演算を行なうことを特徴とする超音波診断装置である。
【００１２】
第４の観点の発明は、第３の観点の発明において、前記他音線相関ウィンドウの相関演算で得られた算出値について、エラーか否かの判定を行なうエラー判定部を備え、前記平均演算部は、エラーの算出値を除いて前記平均演算を行なうことを特徴とする超音波診断装置である。
【００１３】
第５の観点の発明は、第１の観点の発明において、前記直前相関ウィンドウの相関演算で得られた算出値について、エラーか否かの判定を行なうエラー判定部と、前記一の音線上における相関ウィンドウの設定にあたり、前記エラー判定部により、前記直前相関ウィンドウの相関演算で得られた算出値についてエラーと判定された場合、複数の前記他の音線についての前記他音線相関ウィンドウの相関演算で得られた算出値の平均演算を行なう平均演算部と、をさらに備え、前記物理量算出部は、前記平均演算部で得られた平均値に基づいて前記相関ウィンドウの設定を行なうことを特徴とする超音波診断装置である。
【００１４】
第６の観点の発明は、第４，５の観点の発明において、前記エラー判定部は、判定対象の算出値が得られた相関演算における相関係数に基づいてエラーか否かを判定することを特徴とする超音波診断装置である。
【００１５】
第７の観点の発明は、第４，５の観点の発明において、前記エラー判定部は、判定対象の算出値が予め設定された所定の範囲内にない場合にエラーと判定することを特徴とする超音波診断装置である。
【００１６】
第８の観点の発明は、第４，５の観点の発明において、前記エラー判定部は、判定対象の算出値について、他の判定対象の算出値の分布に基づいて、エラーか否かを判定することを特徴とする超音波診断装置である。
【００１７】
第９の観点の発明は、第２〜８のいずれか一の観点の発明において、前記平均演算部は、前記平均演算の対象となる算出値に対し、相関係数に応じた重み付けを行なうことを特徴とする超音波診断装置である。
【００１８】
第１０の観点の発明は、第１の観点の発明において、前記一の音線上における相関ウィンドウの設定にあたり、前記直前相関ウィンドウの相関演算で得られた算出値と、前記他音線相関ウィンドウの相関演算で得られた算出値の中から、前記相関ウィンドウの設定に適した算出値を選択する選択部を備え、前記物理量算出部は、前記選択部によって選択された算出値に基づいて前記相関ウィンドウの設定を行なうことを特徴とする超音波診断装置である。
【００１９】
第１１の観点の発明は、第１０の観点の発明において、前記選択部は、前記相関ウィンドウの設定に適した算出値の選択を、該算出値が得られた相関演算の相関係数又は前記算出値が所定の範囲内にあるか否かに基づいて行なうことを特徴とする超音波診断装置である。
【００２０】
第１２の観点の発明は、第１０，１１の観点の発明において、前記選択部は、前記直前相関ウィンドウの相関演算における相関係数が所定の閾値を超える場合、該直前相関ウィンドウの相関演算によって得られた算出値を前記相関ウィンドウの設定に適した算出値として選択し、一方で前記直前相関ウィンドウの相関演算における相関係数が前記所定の閾値以下である場合、前記他音線相関ウィンドウの相関演算で得られた算出値であって、前記直前相関ウィンドウの相関演算における相関係数よりも高い相関係数の相関演算で得られた算出値又は所定の範囲内の算出値を、前記相関ウィンドウの設定に適した算出値として選択することを特徴とする超音波診断装置である。
【００２１】
第１３の観点の発明は、第１０，１１の観点の発明において、前記選択部は、前記直前相関ウィンドウの相関演算における相関係数及び前記他音線相関ウィンドウの相関演算における相関係数のうち、相関係数が最も高い相関演算によって得られた算出値を、前記相関ウィンドウの設定に適した算出値として選択することを特徴とする超音波診断装置である。
【００２２】
第１４の観点の発明は、第１０，１１の観点の発明において、前記選択部は、前記直前相関ウィンドウの相関演算によって得られた算出値が所定の範囲内である場合、該直前相関ウィンドウの相関演算によって得られた算出値を前記相関ウィンドウの設定に適した算出値として選択し、一方で前記直前相関ウィンドウの相関演算によって得られた算出値が所定の範囲外である場合、前記他音線相関ウィンドウの相関演算で得られた算出値であって、前記直前相関ウィンドウの相関演算における相関係数よりも高い相関係数の相関演算で得られた算出値又は前記所定の範囲内の算出値を、前記相関ウィンドウの設定に適した算出値として選択することを特徴とする超音波診断装置である。
【００２３】
第１５の観点の発明は、第１０，１１の観点の発明において、前記選択部は、複数の前記他の音線についての前記他音線相関ウィンドウの相関演算で得られた算出値の分布に基づいて、前記直前相関ウィンドウの相関演算によって得られた算出値が前記相関ウィンドウの設定に適するか否かを判断し、該相関ウィンドウの設定に適していると判断した場合には前記直前相関ウィンドウの相関演算によって得られた算出値を前記相関ウィンドウの設定に適した算出値として選択し、一方で、前記相関ウィンドウの設定に適していないと判断した場合には、いずれかの前記他音線相関ウィンドウの相関演算で得られた算出値であって、前記直前相関ウィンドウの相関演算における相関係数よりも高い相関係数の相関演算で得られた算出値又は前記所定の範囲内の算出値を、前記相関ウィンドウの設定に適した算出値として選択することを特徴とする超音波診断装置である。
【００２４】
第１６の観点の発明は、第１〜１５のいずれか一の観点の発明において、前記他音線相関ウィンドウは、前記直前相関ウィンドウと生体組織において同じ深さに位置することを特徴とする超音波診断装置である。
【００２５】
第１７の観点の発明は、コンピュータに、生体組織に対し超音波の送受信を行なって得られたエコー信号であって、時間的に異なる二つのフレームに属する同一音線上における二つのエコー信号に、相関ウィンドウを設定し、該相関ウィンドウ間で相関演算を行なって生体組織の弾性に関する物理量を算出する物理量算出機能と、前記物理量に基づいて生体組織の弾性画像データを作成する弾性画像データ作成機能と、を実行させ、前記物理量算出機能にあっては、前記二つのフレームのいずれか一方に属する一の音線上のエコー信号への相関ウィンドウの設定にあたり、該相関ウィンドウと前記一の音線上において隣り合っていて直前に相関演算を行なった直前相関ウィンドウの相関演算と、前記一の音線とは異なる他の音線上における他音線相関ウィンドウの相関演算とに従って、前記相関ウィンドウの設定を行なうことを特徴とする超音波診断装置の制御プログラムである。
【００２６】
第１８の観点の発明は、第１７の観点の発明において、前記一の音線上における相関ウィンドウの設定にあたり、前記直前相関ウィンドウの相関演算で得られた算出値と、前記他音線相関ウィンドウの相関演算で得られた算出値との平均演算を行なう平均演算機能をコンピュータにさらに実行させ、前記物理量算出機能にあっては、前記平均演算機能で得られた平均値に基づいて前記相関ウィンドウの設定を行なうことを特徴とする超音波診断装置の制御プログラムである。
【００２７】
第１９の観点の発明は、第１７の観点の発明において、前記直前相関ウィンドウの相関演算で得られた算出値について、エラーか否かの判定を行なうエラー判定機能と、前記一の音線上における相関ウィンドウの設定にあたり、前記エラー判定機能により、前記直前相関ウィンドウの相関演算で得られた算出値についてエラーと判定された場合、複数の前記他の音線についての前記他音線相関ウィンドウの相関演算で得られた算出値の平均演算を行なう平均演算機能と、をコンピュータにさらに実行させ、前記物理量算出機能にあっては、前記平均演算機能で得られた平均値に基づいて前記相関ウィンドウの設定を行なうことを特徴とする超音波診断装置の制御プログラムである。
【００２８】
第２０の観点の発明は、第１７の観点の発明において、前記一の音線上における相関ウィンドウの設定にあたり、前記直前相関ウィンドウの相関演算で得られた算出値と、前記他音線相関ウィンドウの相関演算で得られた算出値の中から、前記相関ウィンドウの設定に適した算出値を選択する選択機能をコンピュータにさらに実行させ、前記物理量算出機能は、前記選択機能によって選択された算出値に基づいて前記相関ウィンドウの設定を行なうことを特徴とする超音波診断装置の制御プログラムである。
【発明の効果】
【００２９】
本発明によれば、前記物理量算出部は、前記一の音線上における相関ウィンドウの設定にあたり、この相関ウィンドウと前記一の音線上において隣り合う前記直前相関ウィンドウの相関演算と、前記一の音線とは異なる他の音線上における前記他音線相関ウィンドウの相関演算とに従って前記相関ウィンドウの設定を行なう。従って、前記直前相関ウィンドウの相関演算の結果得られる算出値が、生体組織の弾性を正確に反映したものではない場合でも、前記他音線相関ウィンドウの相関演算の算出値が、生体組織の弾性をより正確に反映したものであれば、相関演算を行なう相関ウィンドウ間のマッチングの度合いが従来よりも高くなるように、前記一の音線上における前記相関ウィンドウの設定を行なうことができる。従って、弾性画像に線状のアーチファクトが現れることを抑制することができ、従来よりも生体組織の弾性を正確に反映した弾性画像を得ることができる。
【００３０】
また、本発明によれば、前記一の音線上における相関ウィンドウの設定にあたり、前記平均演算部が、前記直前相関ウィンドウの相関演算で得られた算出値と、前記他音線相関ウィンドウの相関演算で得られた算出値との平均演算を行ない、前記物理量算出部は、得られた平均値に基づいて前記相関ウィンドウの設定を行なう。従って、前記直前相関ウィンドウの相関演算の結果得られる算出値が、生体組織の弾性を正確に反映したものではない場合でも、前記他音線相関ウィンドウの相関演算の算出値が、生体組織の弾性をより正確に反映したものであれば、相関演算を行なう相関ウィンドウ間のマッチングの度合いが従来よりも高くなるように、前記一の音線上における相関ウィンドウの設定を行なうことができる。従って、弾性画像に線状のアーチファクトが現れることを抑制することができ、従来よりも生体組織の弾性を正確に反映した弾性画像を得ることができる。
【００３１】
また、本発明によれば、前記一の音線上における相関ウィンドウの設定にあたり、前記エラー判定部により、前記直前相関ウィンドウの相関演算で得られた算出値についてエラーと判定された場合、前記平均演算部が、複数の前記他の音線についての前記他音線相関ウィンドウの相関演算で得られた算出値の平均演算を行ない、前記物理量算出部は、得られた平均値に基づいて前記相関ウィンドウの設定を行なう。従って、前記直前相関ウィンドウの相関演算の結果得られる算出値が、生体組織の弾性を正確に反映したものではない場合でも、前記他音線相関ウィンドウの相関演算の算出値が、生体組織の弾性をより正確に反映したものであれば、相関演算を行なう相関ウィンドウ間のマッチングの度合いが従来よりも高くなるように、前記一の音線上における相関ウィンドウの設定を行なうことができる。従って、弾性画像に線状のアーチファクトが現れることを抑制することができ、従来よりも生体組織の弾性を正確に反映した弾性画像を得ることができる。
【００３２】
また、本発明によれば、前記選択部は、前記一の音線上における相関ウィンドウの設定にあたり、前記直前相関ウィンドウの相関演算で得られた算出値と、前記他音線相関ウィンドウの相関演算で得られた算出値の中から、前記相関ウィンドウの設定に適した算出値を選択し、選択された算出値に基づいて前記物理量算出部が前記相関ウィンドウの設定を行なう。従って、前記直前相関ウィンドウの相関演算の結果得られる算出値が、生体組織の弾性を正確に反映したものではない場合でも、前記他音線相関ウィンドウの相関演算の算出値が、生体組織の弾性をより正確に反映したものであれば、この他音線相関ウィンドウの相関演算の算出値が、前記相関ウィンドウの設定に適した算出値として選択される。これにより、相関演算を行なう相関ウィンドウ間のマッチングの度合いが従来よりも高くなるように、前記一の音線上における相関ウィンドウの設定を行なうことができる。従って、弾性画像に線状のアーチファクトが現れることを抑制することができ、従来よりも生体組織の弾性を正確に反映した弾性画像を得ることができる。
【図面の簡単な説明】
【００３３】
【図１】本発明に係る超音波診断装置の一実施形態の概略構成を示すブロック図である。
【図２】図１に示す超音波診断装置の弾性画像処理部の構成を示すブロック図である。
【図３】表示部に表示された超音波画像を示す図である。
【図４】弾性画像データを作成する際における相関ウィンドウの設定を説明するための図である。
【図５】弾性画像データを作成する際における相関ウィンドウの設定を説明するための図である。
【図６】弾性画像データを作成する際における相関ウィンドウの設定を説明するための図である。
【図７】第一実施形態の変形例における弾性画像処理部の構成を示すブロック図である。
【図８】第二実施形態における弾性画像処理部の構成を示すブロック図である。
【発明を実施するための形態】
【００３４】
以下、本発明の実施形態について図面に基づいて詳細に説明する。
（第一実施形態）
先ず、第一実施形態について図１〜図６に基づいて説明する。図１に示す超音波診断装置１は、超音波プローブ２、送受信部３、Ｂモード画像処理部４、弾性画像処理部５、合成部６、表示部７を備え、さらに制御部８及び操作部９を備える。
【００３５】
前記超音波プローブ２は、生体組織に対して超音波を送信しそのエコーを受信する。この超音波プローブ２を生体組織の表面に当接させた状態で圧迫と弛緩を繰り返しながら超音波の送受信を行なって取得されたエコー信号に基づいて、後述のように弾性画像が作成される。
【００３６】
前記送受信部３は、前記超音波プローブ２を所定の走査条件で駆動させて音線毎の超音波の走査を行なう。また、前記超音波プローブ２で受信したエコー信号について、整相加算処理等の信号処理を行なう。
【００３７】
ちなみに、前記送受信部３は、Ｂモード画像を作成するための走査と、弾性画像を作成するための走査とを別に行なう。弾性画像を作成するための走査としては、被検体における弾性画像を作成する領域において、同一音線上に二回の走査を行なう。
【００３８】
前記Ｂモード画像処理部４は、前記送受信部３から出力されたエコー信号に対し、対数圧縮処理、包絡線検波処理等のＢモード処理を行い、Ｂモード画像データを作成する。
【００３９】
前記弾性画像処理部５は、前記送受信部３から出力されたエコー信号に基づいて、弾性画像データを作成する。この弾性画像処理部５は、図２に示すように、物理量算出部５１と弾性画像データ作成部５２とを有し、さらに平均演算部５３を有している。
【００４０】
前記物理量算出部５１は、生体組織における各部の弾性に関する物理量として、前記超音波プローブ２による圧迫とその弛緩によって生じた生体組織における各部の変形による変位（以下、単に「変位」と云う）を算出する（変位算出機能）。前記物理量算出部５１は、後述の説明で用いる図４に示すように時間的に異なる二つのフレーム（ｉ），（ｉｉ）に属する同一音線上における二つのエコー信号に基づいて変位を算出する。詳細は後述する。前記物理量算出部５１は本発明における物理量算出部の実施の形態の一例であり、また前記変位算出機能は本発明における物理量算出機能の実施の形態の一例である。
【００４１】
前記弾性画像データ作成部５２は、前記物理量算出部５１によって算出された変位を色相情報に変換し、弾性画像を作成する領域における弾性画像データを作成する（弾性画像データ作成機能）。前記弾性画像データ作成部５２は本発明における弾性画像データ作成部の実施の形態の一例であり、また前記弾性画像データ作成機能は本発明における弾性画像データ作成機能の実施の形態の一例である。
【００４２】
ここで、本例では、図３に示すように前記表示部７に表示されたＢモード画像ＢＧ上に関心領域（ＲＯＩ：Ｒｅｇｉｏｎ  Ｏｆ  Ｉｎｔｅｒｅｓｔ）Ｒが設定され、この関心領域Ｒについて前記弾性画像データが作成される。ただし、本発明は、このように前記Ｂモード画像ＢＧの一部について弾性画像を作成する場合に限られるものではなく、前記Ｂモード画像ＢＧの全体について前記弾性画像データを作成してもよい。
【００４３】
前記平均演算部５３は、直前相関ウィンドウの相関演算で得られた変位（後述の変位Ｘｃ）と、前記他音線相関ウィンドウの相関演算で得られた変位（後述の変位Ｘａ，Ｘｂ，Ｘｄ，Ｘｅ）との平均演算を行なう（平均演算機能）。詳細については後述する。前記平均演算部５３は、本発明における平均演算部の実施の形態の一例である。
【００４４】
前記Ｂモード画像処理部４で作成されたＢモード画像データと、前記弾性画像処理部５で作成された弾性画像データは、前記合成部６で合成される。具体的には、この合成部６は、一フレーム分の前記Ｂモード画像データと前記弾性画像データとを加算処理し、前記表示部７に表示する一フレーム分の超音波画像データを作成する。そして、前記合成部６で得られた超音波画像データは、図３に示すように白黒のＢモード画像ＢＧとカラーの弾性画像ＥＧとが合成された超音波画像Ｇとして前記表示部７に表示される。本例では、前記弾性画像ＥＧは、前記関心領域Ｒ内に半透明で（背景のＢモード画像が透けた状態で）表示される。
【００４５】
前記制御部８は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）で構成され、図示しない記憶部に記憶された制御プログラムを読み出し、前記変位算出機能、前記弾性画像データ作成機能及び前記平均演算機能を実行させ、その他前記超音波診断装置１の各部における機能を実行させる。また、前記操作部９は、操作者が指示や情報を入力するためのキーボード及びポインティングデバイス（図示省略）などを含んで構成されている。
【００４６】
さて、本例の超音波診断装置１の作用について説明する。先ず、前記送受信部３は、前記超音波プローブ２から被検体の生体組織へ超音波を送信させ、そのエコー信号を取得する。このとき、前記超音波プローブ２により、被検体への圧迫とその弛緩を繰り返しながら超音波の送受信を行う。
【００４７】
そして、エコー信号が取得されると、前記Ｂモード画像処理部４は、前記送受信部３からのエコー信号に基づいてＢモード画像データを作成する。また、前記弾性画像処理部５は、前記送受信部３からのエコー信号に基づいて弾性画像データを作成する。前記Ｂモード画像データと前記弾性画像データは、前記合成部６で合成され、図３に示すようにＢモード画像ＢＧと弾性画像ＥＧとが合成された超音波画像Ｇが前記表示部７に表示される。
【００４８】
前記弾性画像処理部５における弾性画像データの作成について詳細に説明する。図４に示す前記フレーム（ｉ），（ｉｉ）は、複数本の音線上において取得されたエコー信号からなる。図４では、前記フレーム（ｉ）における複数本の音線の一部として、五本の音線Ｌ１ａ，Ｌ１ｂ，Ｌ１ｃ，Ｌ１ｄ，Ｌ１ｅが示され、また前記フレーム（ｉｉ）において前記音線Ｌ１ａ〜Ｌ１ｅに対応する音線として、音線Ｌ２ａ，Ｌ２ｂ，Ｌ２ｃ，Ｌ２ｄ，Ｌ２ｅが示されている。すなわち、前記音線Ｌ１ａ及び前記音線Ｌ２ａ、前記音線Ｌ１ｂ及び前記音線Ｌ２ｂ、前記音線Ｌ１ｃ及び前記音線Ｌ２ｃ、前記音線Ｌ１ｄ及び前記音線Ｌ２ｄ、前記音線Ｌ１ｅ及び前記音線Ｌ２ｅは、異なる二つのフレームに属する同一音線に該当する。また、図４においてＲ（ｉ），Ｒ（ｉｉ）は、前記関心領域Ｒに対応する領域を示している。
【００４９】
ここでは、前記音線Ｌ１ｃ，Ｌ２ｃ上における弾性画像データを作成する場合を例にして説明する。前記物理量算出部５１は、前記音線Ｌ１ｃ，Ｌ２ｃ上のエコー信号Ｓ１，Ｓ２（図示省略）のそれぞれに、相関ウィンドウＷ１，Ｗ２を設定し、この相関ウィンドウＷ１，Ｗ２間で相関演算を行なって前記変位を算出する。そして、この変位に基づいて、前記弾性画像データ作成部５２が一画素分の弾性画像データを作成する。
【００５０】
前記相関ウィンドウＷ１，Ｗ２の設定について説明すると、前記物理量算出部５１は、前記音線Ｌ１ｃ上における前記領域Ｒ（ｉ）の上端部１００から下端部１０１に向かって、相関ウィンドウＷ１を順次設定し、また前記音線Ｌ２ｃ上における前記領域Ｒ（ｉｉ）の上端部１００から下端部１０１に向かって、相関ウィンドウＷ２を順次設定する。
【００５１】
さらに詳細に説明する。ここでは、図４に示すように、音線Ｌ１ｃ，Ｌ２ｃ上のエコー信号Ｓ１，Ｓ２に相関ウィンドウＷ１１ｃ，Ｗ２１ｃを設定する場合を例に挙げて説明する。ちなみに、前記フレーム（ｉ）に属するエコー信号については、前記音線Ｌ１ｃ上のエコー信号Ｓ１への設定と同様にして相関ウィンドウが設定され、また前記フレーム（ｉｉ）に属するエコー信号については、前記音線Ｌ２ｃ上のエコー信号Ｓ２への設定と同様にして相関ウィンドウが設定される。
【００５２】
図４において、Ｗ１０ｃ，Ｗ２０ｃは、前記相関ウィンドウＷ１１ｃ，Ｗ２１ｃと前記音線Ｌ１ｃ，Ｌ２ｃ上において隣り合っていて直前に相関演算を行なった相関ウィンドウである。ここで、「隣り合う相関ウィンドウ」とは、相関ウィンドウの設定の順番が連続している相関ウィンドウをいう。また、設定しようとする相関ウィンドウ（ここでは前記相関ウィンドウＷ１１ｃ，Ｗ２１ｃ）と同一音線上において隣り合っていて直前に相関演算を行なった相関ウィンドウを直前相関ウィンドウとする。ここでは、直前相関ウィンドウは前記相関ウィンドウＷ１０ｃ，Ｗ２０ｃである。すなわち、この相関ウィンドウＷ１０ｃ，Ｗ２０ｃは、本発明における直前相関ウィンドウの実施の形態の一例である。
【００５３】
先ず、前記物理量算出部５１は、前記直前相関ウィンドウＷ１０ｃから予め定められた所定のデータ数分だけ移動させて前記相関ウィンドウＷ１１ｃを設定する。ここでは、前記相関ウィンドウＷ１０ｃのウィンドウ幅分のデータ数だけ移動させて前記相関ウィンドウＷ１１ｃを設定する。この相関ウィンドウＷ１１ｃのウィンドウ幅は、前記相関ウィンドウＷ１０ｃのウィンドウ幅と同じである。従って、前記音線Ｌ１ｃ上のエコー信号Ｓ１には、前記上端部１００から前記下端部１０１へ向かって、同じウィンドウ幅の相関ウィンドウＷ１が順次設定される。
【００５４】
次に、前記相関ウィンドウＷ２１ｃの設定について説明する。前記相関ウィンドウＷ２１ｃを設定するにあたっては、図５に示すように先ず前記音線Ｌ１ｃ，Ｌ２ｃとは異なる前記音線Ｌ１ａ，Ｌ２ａ、前記音線Ｌ１ｂ，Ｌ２ｂ、前記音線Ｌ１ｄ，Ｌ２ｄ及び前記音線Ｌ１ｅ，Ｌ２ｅに、前記相関ウィンドウＷ１０ｃ，Ｗ２０ｃと生体組織において同じ深さとなる位置に相関ウィンドウを設定する。具体的には、前記音線Ｌ１ａ，Ｌ２ａには相関ウィンドウＷ１０ａ，Ｗ２０ａを設定し、前記音線Ｌ１ｂ，Ｌ２ｂには相関ウィンドウＷ１０ｂ，Ｗ２０ｂを設定し、前記音線Ｌ１ｄ，Ｌ２ｄには相関ウィンドウＷ１０ｄ，Ｗ２０ｄを設定し、前記音線Ｌ１ｅ，Ｌ２ｅには相関ウィンドウＷ１０ｅ，Ｗ２０ｅを設定する。
【００５５】
ここで、相関ウィンドウを設定する音線とは異なる他の音線に設定される相関ウィンドウを他音線相関ウィンドウとする。ここでは、他音線相関ウィンドウは、前記相関ウィンドウＷ１０ａ，Ｗ２０ａ、前記相関ウィンドウＷ１０ｂ，Ｗ２０ｂ、前記相関ウィンドウＷ１０ｄ，Ｗ２０ｄ、前記相関ウィンドウＷ１０ｅ，Ｗ２０ｅであり、前記直前相関ウィンドウと生体組織において同じ深さに設定される。前記相関ウィンドウＷ１０ａ，Ｗ２０ａ、前記相関ウィンドウＷ１０ｂ，Ｗ２０ｂ、前記相関ウィンドウＷ１０ｄ，Ｗ２０ｄ、前記相関ウィンドウＷ１０ｅ，Ｗ２０ｅは、本発明における他音線相関ウィンドウの実施の形態の一例である。
【００５６】
ここで、他音線相関ウィンドウである前記相関ウィンドウＷ１０ａ，Ｗ２０ａ、前記相関ウィンドウＷ１０ｂ，Ｗ２０ｂ、前記相関ウィンドウＷ１０ｄ，Ｗ２０ｄ、前記相関ウィンドウＷ１０ｅ，Ｗ２０ｅは、直前相関ウィンドウである前記相関ウィンドウＷ１０ｃ，Ｗ２０ｃと生体組織において同じ深さ位置に設定されるものに限られるものではない。ただし、深さ方向の位置が離れすぎないことが望ましい。
【００５７】
ちなみに、ここでは前記音線Ｌ１ｃ，Ｌ２ｃは本発明における一の音線の実施の形態の一例である。また、前記音線Ｌ１ａ，Ｌ２ａ、音線Ｌ１ｂ，Ｌ２ｂ、音線Ｌ１ｄ，Ｌ２ｄ、音線Ｌ１ｅ，Ｌ２ｅは、本発明における他の音線の実施の形態の一例である。
【００５８】
次に、前記物理量算出部５１は、前記相関ウィンドウＷ１０ａ，Ｗ２０ａ間の相関演算と、前記相関ウィンドウＷ１０ｂ，Ｗ２０ｂ間の相関演算と、前記相関ウィンドウＷ１０ｄ，Ｗ２０ｄ間の相関演算と、前記相関ウィンドウＷ１０ｅ，Ｗ２０ｅ間の相関演算とを行ない、それぞれについて変位Ｘａ，Ｘｂ，Ｘｄ，Ｘｅを算出する。そして、前記平均演算部５３は、これら変位Ｘａ，Ｘｂ，Ｘｄ，Ｘｅと前記相関ウィンドウＷ１０ｃ，Ｗ２０ｃ間の相関演算で得られた変位Ｘｃとの平均演算を行ない、各変位の平均値ＸＡＶを算出する。
【００５９】
なお、他音線相関ウィンドウについての前記変位Ｘａ，Ｘｂ，Ｘｄ，Ｘｅは、前記平均演算を行なうにあたり、前記物理量算出部５１によって算出してもよいが、すでに他音線相関ウィンドウについて弾性画像データを作成する時に前記変位Ｘａ，Ｘｂ，Ｘｄ，Ｘｅが算出されている場合には、改めて算出することなく、この時に算出された値を前記平均演算に用いればよい。
【００６０】
ここで、前記平均演算を行なう際には、前記変位Ｘａ，Ｘｂ，Ｘｃ，Ｘｄ，Ｘｅに対し、相関係数に応じた重み付けを行ってもよい。すなわち、前記相関ウィンドウＷ１０ａ，Ｗ２０ａ間の相関演算における相関係数をＣａ、前記相関ウィンドウＷ１０ｂ，Ｗ２０ｂ間の相関演算における相関係数をＣｂ、前記相関ウィンドウＷ１０ｃ，Ｗ２０ｃ間の相関演算における相関係数をＣｃ、前記相関ウィンドウＷ１０ｄ，Ｗ２０ｄ間の相関演算における相関係数をＣｄ、前記相関ウィンドウＷ１０ｅ，Ｗ２０ｅ間の相関演算における相関係数をＣｅとすると、前記相関係数Ｃａに応じた重み係数を前記変位Ｘａに乗算し、前記相関係数Ｃｂに応じた重み係数を前記変位Ｘｂに乗算し、前記相関係数Ｃｃに応じた重み係数を前記変位Ｘｃに乗算し、前記相関係数Ｃｄに応じた重み係数を前記変位Ｘｄに乗算し、前記相関係数Ｃｅに応じた重み係数を前記変位Ｘｅに乗算して平均演算を行なう。重み係数は、相関係数が高くなるほど大きくなるような係数とする。
【００６１】
ちなみに、前記弾性画像データ作成部５２は、前記相関ウィンドウＷ１０ｃ，Ｗ２０ｃの弾性画像データの作成にあっては、前記変位Ｘｃに基づいて弾性画像データを作成するのではなく、前記平均値ＸＡＶに基づいて弾性画像データを作成してもよい。
【００６２】
前記平均演算部５３により前記平均値ＸＡＶが算出されると、前記物理量算出部５１は、この平均値ＸＡＶに基づいて、前記相関ウィンドウＷ２０ｃからの移動量を決定して図６に示すように前記相関ウィンドウＷ２１ｃを設定する。このようにして設定される前記相関ウィンドウＷ２１ｃは、前記相関ウィンドウＷ２０ｃと重複する部分を有し、またこの相関ウィンドウＷ２０ｃのウィンドウ幅とは異なるウィンドウ幅になっている。
【００６３】
ここで、前記相関ウィンドウＷ１１ｃは、前記相関ウィンドウＷ２１ｃが設定される前に設定されてもよいし、前記相関ウィンドウＷ２１ｃが設定された後に設定されてもよい。そして、前記物理量算出部５１は、前記相関ウィンドウＷ１１ｃ，Ｗ２１ｃを設定すると、これら相関ウィンドウＷ１１ｃ，Ｗ２１ｃ間で相関演算を行ない、変位を算出する。
【００６４】
ちなみに、前記音線Ｌ２ｃ上においては、前記上端部１００から前記下端部１０１に向かって、以下同様にして相関ウィンドウＷ２が順次設定される。すなわち、音線Ｌ２ｃ上において相関ウィンドウを設定するにあたり、直前に相関演算を行なった直前相関ウィンドウ間で得られた変位と、この直前相関ウィンドウと同じ深さに位置する他音線（前記音線Ｌ１ａ，Ｌ２ａ、音線Ｌ１ｂ，Ｌ２ｂ、音線Ｌ１ｄ，Ｌ２ｄ、音線Ｌ１ｅ，Ｌ２ｅ）の相関ウィンドウ間の相関演算で得られた変位との平均演算を行ない、得られた平均の変位に基づいて相関ウィンドウＷ２が順次設定される。このようにして前記上端部１００から前記下端部１０１へ向かって前記音線Ｌ２ｃ上のエコー信号Ｓ２に順次設定された相関ウィンドウＷ２は、必ずしも同じウィンドウ幅とならず、また互いに重複する部分を有する場合がある。
【００６５】
以上説明した本例の超音波診断装置１によれば、前記直前相関ウィンドウの相関演算と前記他音線相関ウィンドウの相関演算とに従って前記フレーム（ｉｉ）における相関ウィンドウの設定が行なわれる。例えば、前記相関ウィンドウＷ２１ｃは、前記相関ウィンドウＷ１０ｃ，Ｗ２０ｃ間で得られた変位Ｘｃと、前記相関ウィンドウＷ１０ａ，Ｗ２０ａ間で得られた変位Ｘａと、前記相関ウィンドウＷ１０ｂ，Ｗ２０ｂ間で得られた変位Ｘｂと、前記相関ウィンドウＷ１０ｄ，Ｗ２０ｄ間で得られた変位Ｘｄと、前記相関ウィンドウＷ１０ｅ，Ｗ２０ｅ間で得られた変位Ｘｅとの平均値ＸＡＶに基づいて設定される。従って、前記相関ウィンドウＷ１０ｃ，Ｗ２０ｃ間で得られた変位Ｘｃが生体組織の弾性を正確に反映したものではない場合でも、前記変位Ｘａ，Ｘｂ，Ｘｄ，Ｘｅが生体組織の弾性をより正確に反映したものであれば、従来のように前記変位Ｘｃのみに基づいて前記相関ウィンドウＸ２１の設定を行なう場合と比べれば、前記相関ウィンドウＷ１１ｃとのマッチング度合いが高くなるように前記相関ウィンドウＷ２１ｃを設定することができる。従って、弾性画像に線状のアーチファクトが現れることを抑制することができ、従来よりも生体組織の弾性を正確に反映した弾性画像を得ることができる。
【００６６】
次に、第一実施形態の変形例について説明する。この変形例では、図７に示すように、前記弾性画像処理部５は、エラー判定部５４をさらに有している。このエラー判定部５４は、前記平均演算部５３による前記平均演算を行なうにあたり、この平均演算に用いる変位を対象としてこれら変位がエラーか否かを判定する（エラー判定機能）。そして、前記平均演算部５３は、エラーとされた変位を除いて平均演算を行なう。前記エラー判定部５４は、本発明におけるエラー判定部の実施の形態の一例である。
【００６７】
例えば、前記エラー判定部５４は、前記変位Ｘａ，Ｘｂ，Ｘｃ，Ｘｄ，Ｘｅの平均演算を行なうにあたり、これら変位Ｘａ〜Ｘｅについてエラーか否かを判定する。そして、前記平均演算部５３は、前記変位Ｘａ〜Ｘｅのうち前記エラー判定部５４でエラーと判定された変位を除いて平均演算を行なう。すなわち、前記平均演算部５３は、前記他音線相関ウィンドウの相関演算で得られた変位であるＸａ，Ｘｂ，Ｘｄ，Ｘｅのいずれかがエラーと判定された場合、エラーと判定されたものを除く前記他音線相関ウィンドウにおける変位と、前記直前相関ウィンドウにおける変位Ｘｃとの平均演算を行なう。また、前記平均演算部５３は、前記変位Ｘｃがエラーと判定された場合、前記変位Ｘａ，Ｘｂ，Ｘｄ，Ｘｅの平均演算を行なう。
【００６８】
ちなみに、前記変位Ｘｃのエラー判定は、前記相関ウィンドウＷ１０ｃ，Ｗ２０ｃについての弾性画像データの作成前に行ってもよい。この場合において、前記変位Ｘｃがエラーと判定された場合、前記弾性画像データ作成部５２は、前記相関ウィンドウＷ１０ｃ，Ｗ２０ｃの弾性画像データの作成にあっては、前記変位Ｘｃに基づいて弾性画像データを作成するのではなく、前記前記変位Ｘａ，Ｘｂ，Ｘｄ，Ｘｅの平均値に基づいて弾性画像データを作成してもよい。
【００６９】
前記エラー判定部５４による判定手法の一例としては、判定対象の変位が得られた相関演算における相関係数Ｃ（０＜Ｃ＜１）に基づいてエラーか否かを判定する手法が挙げられる。この場合、相関係数Ｃについて所定の閾値ＣＴＨを予め設定しておき、この閾値ＣＴＨよりも相関係数Ｃが小さい場合、前記エラー判定部５４はエラーと判定する。例えば、前記エラー判定部５４は、前記変位Ｘａについてエラーか否かを判定する場合、前記相関ウィンドウＷ１０ａ，Ｗ２０ａ間で行なわれた相関演算の相関係数Ｃａを前記閾値ＣＴＨと比較し、前記相関係数Ｃａが閾値ＣＴＨよりも低かった場合、前記変位Ｘａをエラーと判定する。
【００７０】
また、前記エラー判定部５４による判定手法の他例としては、判定対象の変位が、予め設定された所定の範囲内にない場合にエラーと判定する手法が挙げられる。ここで、所定の範囲は、例えば操作者によって設定され、操作者が、通常得られると思料する変位の範囲である。
【００７１】
さらに、前記エラー判定部５４による判定手法の他例としては、判定対象の変位について、他の判定対象の変位の分布に基づいてエラーか否かを判定する手法が挙げられる。具体的には、前記エラー判定部５４は、判定対象の変位が、他の判定対象の変位の分布に対して著しく異なる場合、エラーと判定する。著しく異なると判定される範囲は、操作者の判断により予め設定される。
【００７２】
例えば、前記エラー判定部５４は、前記変位Ｘａについてエラーか否かを判定する場合、前記変位Ｘｂ，Ｘｃ，Ｘｄ，Ｘｅの平均値を求め、この平均値に対し操作者によって設定された±ｎ％の範囲内に前記変位Ｘａが入っていない場合にエラーと判定する。このように、前記エラー判定部５４は、前記変位Ｘｂ，Ｘｃ，Ｘｄ，Ｘｅの分布に基づいて、エラーか否かを判定する。
【００７３】
（第二実施形態）
次に、第二実施形態について説明する。この第二実施形態は第一実施形態と基本的構成を同じくし、第一実施形態と同様の事項については説明を省略する。
【００７４】
本例においても、前記物理量算出部５１は、第一実施形態と同様に、前記直前相関ウィンドウの相関演算と前記他音線相関ウィンドウの相関演算とに従って、前記フレーム（ｉｉ）における相関ウィンドウの設定を行なうが、具体的な設定手法が異なっている。
【００７５】
具体的に説明すると、本例では、図８に示すように、前記弾性画像処理部５は、前記物理量算出部５１及び前記弾性画像データ作成部５２のほか、選択部５５を有している。この選択部５５は、前記直前相関ウィンドウの変位と、前記他音線相関ウィンドウの相関演算で得られた変位の中から、前記相関ウィンドウの設定に適した変位を選択する（選択機能）。この選択部５５は、本発明における選択部の実施の形態の一例である。そして、前記物理量算出部５１は、前記選択部５５によって選択された変位に基づいて、前記フレーム（ｉｉ）における相関ウィンドウの設定を行なう。なお、前記フレーム（ｉ）における相関ウィンドウの設定は第一実施形態と同様である。
【００７６】
例えば、前記選択部５５は、前記直前相関ウィンドウの相関演算における相関係数が所定の閾値を超えている場合、この直前相関ウィンドウの相関演算で得られた変位を、前記相関ウィンドウの設定に適した変位として選択する。一方、前記選択部５５は、前記直前相関ウィンドウの相関演算における相関係数が所定の閾値以下である場合、前記他音線相関ウィンドウにおいて、前記直前相関ウィンドウの相関演算における相関係数よりも高い相関係数の相関演算によって得られた変位を、前記相関ウィンドウの設定に適した変位として選択する。
【００７７】
具体的に説明すると、前記相関ウィンドウＷ２１ｃを設定する場合、前記相関ウィンドウＷ１０ｃ，Ｗ２０ｃの相関演算における相関係数Ｃｃが前記閾値ＣＴＨを超えている場合、前記選択部５５により選択される変位は、前記相関ウィンドウＷ１０ｃ，Ｗ２０ｃ間で得られた変位Ｘｃであり、前記物理量算出部５１は、この変位Ｘｃに基づいて前記相関ウィンドウＷ２１を設定する。
【００７８】
一方、前記相関係数Ｃｃが前記閾値ＣＴＨ以下である場合、前記選択部５５は、前記音線Ｌ１ｃ，Ｌ２ｃの左隣の音線である前記音線Ｌ１ｂ，Ｌ２ｂにおける前記相関ウィンドウＷ１０ｂ，Ｗ２０ｂ間で行なわれた相関演算の相関係数Ｃｂを前記相関係数Ｃｃと比較する。そして、前記相関係数Ｃｂが前記相関係数Ｃｃよりも高い場合、前記選択部５５は、前記相関ウィンドウＷ１０ｂ，Ｗ２０ｂ間で得られた変位Ｘｂを前記相関ウィンドウＷ２１ｃを設定するための変位として選択する。一方、前記相関係数Ｃｂが前記相関係数Ｃｃ以下である場合、前記選択部５５は、前記音線Ｌ１ｃ，Ｌ２ｃの右隣の音線である前記音線Ｌ１ｄ，Ｌ２ｄにおける前記相関ウィンドウＷ１０ｄ，Ｗ２０ｄ間で行なわれた相関演算の相関係数Ｃｄを前記相関係数Ｃｃと比較する。そして、前記相関係数Ｃｄが前記相関係数Ｃｃよりも高い場合、前記選択部５５は、前記相関ウィンドウＷ１０ｄ，Ｗ２０ｄ間で得られた変位Ｘｄを前記相関ウィンドウＷ２１ｃを設定するための変位として選択する。一方、前記相関係数Ｃｄが前記相関係数Ｃｃ以下である場合、前記選択部５５は、さらに異なる音線について、上述の処理を繰り返し、前記相関係数Ｃｃよりも高い相関係数となった相関演算で得られた変位を、前記相関ウィンドウＷ２１を設定するための変位として選択する。
【００７９】
ただし、前記選択部５５は、前記相関係数Ｃｃが前記閾値ＣＴＨ以下である場合、前記相関係数Ｃａ，Ｃｂ，Ｃｄ，Ｃｅの中で、最も高い相関係数（相関係数Ｃｃよりも高い）の相関演算によって得られた変位を前記相関ウィンドウＷ２１ｃを設定するための変位として選択してもよい。
【００８０】
また、前記選択部５５は、前記相関係数Ｃｃが前記閾値ＣＴＨ以下である場合、前記変位Ｘａ，Ｘｂ、Ｘｄ，Ｘｅの中で、予め設定された所定の範囲内にあるいずれかの変位を、前記相関ウィンドウＷ２１ｃを設定するための変位として選択してもよい。
【００８１】
以上説明した本例によれば、例えば、前記相関ウィンドウＷ２１ｃは、前記相関ウィンドウＷ１０ｃ，Ｗ２１ｃの相関演算で得られた変位Ｘｃと他音線相関ウィンドウ（前記相関ウィンドウＷ１０ｂ，Ｗ２０ｂ、前記相関ウィンドウＷ１０ｄ，Ｗ２０ｄ等）の相関演算で得られた変位の中から、前記選択部５５によって選択された変位に基づいて設定される。従って、前記変位Ｘｃが生体組織の弾性を正確に反映したものではない場合でも、前記他音線相関ウィンドウの相関演算で得られた変位が生体組織の弾性をより正確に反映したものであれば、この他音線相関ウィンドウの相関演算の変位が前記相関ウィンドウＷ２１ｃの設定に適した変位として選択される。これにより、前記相関ウィンドウＷ１１ｃとのマッチング度合いが従来よりも高くなるように前記前記相関ウィンドウＷ２１ｃを設定することができる。従って、弾性画像に線状のアーチファクトが現れることを抑制することができ、従来よりも生体組織の弾性を正確に反映した弾性画像を得ることができる。
【００８２】
次に、第二実施形態の変形例について説明する。先ず、第一変形例について説明すると、前記選択部５５は、前記直前相関ウィンドウの相関演算における相関係数及び前記他音線相関ウィンドウの相関演算における相関係数のうち、相関係数が最も高い相関演算によって得られた変位を、相関ウィンドウの設定に適した変位として選択してもよい。例えば、前記選択部５５は、前記相関ウィンドウＷ１０ａ，Ｗ２０ａの相関演算における相関係数Ｃａ、前記相関ウィンドウＷ１０ｂ，Ｗ２０ｂの相関演算における相関係数Ｃｂ、前記相関ウィンドウＷ１０ｃ，Ｗ１０ｃの相関演算における相関係数Ｃｃ，前記相関ウィンドウＷ１０ｄ，Ｗ２０ｄの相関演算における相関係数Ｃｄ、前記相関ウィンドウＷ１０ｅ，Ｗ２０ｅの相関演算における相関係数Ｃｅのうち、相関係数が最も高い相関演算によって得られた変位を、前記相関ウィンドウＷ２１ｃを設定する変位として選択する。
【００８３】
次に第二変形例について説明する。この第二変形例では、前記選択部５５は、前記直前相関ウィンドウの相関演算によって得られた変位が所定の範囲内の変位である場合、この直前相関ウィンドウの相関演算によって得られた変位を、前記相関ウィンドウの設定に適した変位として選択する。一方、前記選択部５５は、前記直前相関ウィンドウの相関演算によって得られた変位が所定の範囲外の変位である場合、前記他音線相関ウィンドウの相関演算によって得られた変位であって、前記所定の範囲内の変位を、前記相関ウィンドウの設定に適した変位として選択する。
【００８４】
例えば、前記相関ウィンドウＷ２１ｃを設定する場合、前記相関ウィンドウＷ１０ｃ，Ｗ２０ｃの相関演算によって得られた変位Ｘｃが、予め設定された所定の範囲内である場合、前記選択部５５は、この変位Ｘｃを前記相関ウィンドウＷ２１ｃを設定する変位として選択する。
【００８５】
一方、前記相関ウィンドウＷ１０ｃ，Ｗ２０ｃの相関演算によって得られた変位Ｘｃが、前記所定の範囲内にない場合、前記選択部５５は、前記変位Ｘａ，Ｘｂ，Ｘｄ，Ｘｅのうち、前記所定の範囲内にあるいずれかの変位を、前記相関ウィンドウＷ２１ｃを設定する変位として選択する。この場合において、前記変位Ｘａ，Ｘｂ，Ｘｄ，Ｘｅのうち、相関係数が前記相関係数Ｃｃよりも高い相関演算によって得られたいずれかの変位を任意に選択してもよい。また、前記変位Ｘａ，Ｘｂ，Ｘｄ，Ｘｅのうち、前記相関係数Ｃｃよりも高い相関係数の相関演算の中で最も高い相関係数の相関演算によって得られた変位を選択してもよい。
【００８６】
ちなみに、前記所定の範囲は、本例においても、例えば操作者によって設定され、操作者が、通常得られると思料する変位の範囲である。
【００８７】
次に第三変形例について説明する。前記選択部５５は、複数の前記他の音線についての他音線相関ウィンドウの相関演算によって得られた変位の分布に基づいて、前記直前相関ウィンドウの相関演算によって得られた変位が前記相関ウィンドウの設定に適しているか否かを判断する。そして、前記直前相関ウィンドウの相関演算によって得られた変位が、前記他音線相関ウィンドウの相関演算によって得られた変位の分布に対して著しく異なるものでない場合、前記選択部５５は、前記直前相関ウィンドウの相関演算によって得られた変位を、前記相関ウィンドウの設定に適した変位として選択する。一方、前記直前相関ウィンドウの相関演算によって得られた変位が、前記他音線相関ウィンドウの相関演算によって得られた変位の分布に対して著しく異なる場合、前記選択部５５は、前記他音線相関ウィンドウの相関演算によって得られた変位であって、所定の範囲内の変位を、前記相関ウィンドウの設定に適した変位として選択する。ちなみに、著しく異なると判定される範囲は、操作者の判断により予め設定される。
【００８８】
例えば、前記相関ウィンドウＷ２１ｃを設定する場合、前記相関ウィンドウＷ１０ｃ，Ｗ２０ｃの相関演算によって得られた変位Ｘｃが、前記変位Ｘａ，Ｘｂ，Ｘｄ，Ｘｅの平均値に対し操作者によって設定された±ｎ％の範囲内である場合、前記変位Ｘｃが前記相関ウィンドウＷ２１ｃの設定に適した相関ウィンドウとして選択される。
【００８９】
一方、前記変位Ｘｃが、前記変位Ｘａ，Ｘｂ，Ｘｄ，Ｘｅの平均値に対し操作者によって設定された±ｎ％の範囲内に入っていない場合、前記変位Ｘａ，Ｘｂ，Ｘｄ，Ｘｅのいずれかの変位であって予め設定された所定の範囲内の変位が、前記相関ウィンドウＷ２１ｃの設定に適した相関ウィンドウとして選択される。この場合、前記変位Ｘａ，Ｘｂ，Ｘｄ，Ｘｅのうち、前記相関係数Ｃｃよりも高い相関係数又は所定の閾値ＣＴＨを超える相関係数の相関演算で得られたいずれかの変位を任意に選択してもよいし、前記相関係数Ｃｃよりも高い相関係数の相関演算の中で最も相関係数が高い相関演算によって得られた変位を選択してもよい。また、前記変位Ｘａ，Ｘｂ，Ｘｄ，Ｘｅのうち、予め設定された所定の範囲内にあるいずれかの変位を選択してもよい。
【００９０】
以上、本発明を前記各実施形態によって説明したが、本発明はその主旨を変更しない範囲で種々変更実施可能なことはもちろんである。例えば、前記物理量算出部５１は、生体組織の弾性に関する物理量として、生体組織の変形による変位の代わりに生体組織の歪みや弾性率を算出してもよい。
【００９１】
また、第二実施形態において、前記相関ウィンドウＷ２１ｃの設定に適した変位として前記Ｘｃが選択されなかった場合、前記弾性画像データ作成部５２は、前記相関ウィンドウＷ１０ｃ，Ｗ２０ｃの弾性画像データの作成にあたり、前記変位Ｘｃに基づいて弾性画像を作成するのではなく、前記相関ウィンドウＷ２１ｃの設定に適するものとして選択された変位に基づいて弾性画像データの作成を行ってもよい。
【符号の説明】
【００９２】
１  超音波診断装置
５１  物理量算出部
５２  弾性画像データ作成部
５３  平均演算部
５４  エラー判定部
５５  選択部
Ｗ１０ｃ，Ｗ２０ｃ  相関ウィンドウ（直前相関ウィンドウ）
Ｗ１０ａ，Ｗ２０ａ，Ｗ１０ｂ，Ｗ２０ｂ，Ｗ１０ｄ，Ｗ２０ｄ，Ｗ１０ｅ，Ｗ２０ｅ  相関ウィンドウ（他音線相関ウィンドウ）
【図１】
2011024654.tif 000003
【図２】
2011024654.tif 000004
【図３】
2011024654.tif 000005
【図４】
2011024654.tif 000006
【図５】
2011024654.tif 000007
【図６】
2011024654.tif 000008
【図７】
2011024654.tif 000009
【図８】
2011024654.tif 000010