(57)【要約】【課題】メモリや回路、また処理時間を要することなく最適音速を求める。【解決手段】被検者に向けて超音波を送信すると共に、被検者から反射された超音波信号を受信することにより受信信号を出力する複数の素子が配列された超音波プローブと、前記被検者に向けて送信する超音波の実音速に対して予め設定された仮定音速を変化させる手段と、前記仮定音速を変化させ、前記受信信号から前記仮定音速に基づく遅延でフォーカスして得られるＲＦ信号により微小構造物の判定を行い、該微小構造物と判定されたＲＦ信号の位相情報から、被検者の超音波音速である最適音速を判定する最適音速判定手段とを有することを特徴とする超音波診断装置を提供することにより前記課題を解決する。【選択図】図１
2009101145.tif 000002
【特許請求の範囲】
【請求項１】
被検者に向けて超音波を送信すると共に、被検者から反射された超音波信号を受信することにより受信信号を出力する複数の素子が配列された超音波プローブと、
前記被検者に向けて送信する超音波の実音速に対して予め設定された仮定音速を変化させる手段と、
前記仮定音速を変化させ、前記受信信号から前記仮定音速に基づく遅延でフォーカスして得られるＲＦ信号により微小構造物の判定を行い、該微小構造物と判定されたＲＦ信号の位相情報から、被検者の超音波音速である最適音速を判定する最適音速判定手段と、
を有することを特徴とする超音波診断装置。
【請求項２】
前記最適音速判定手段は、前記微小構造物と判定されたＲＦ信号の振幅情報から前記最適音速を判定することを特徴とする請求項１に記載の超音波診断装置。
【請求項３】
前記最適音速判定手段は、前記微小構造物と判定されたＲＦ信号の、前記超音波プローブの素子の配列方向の位相変化から前記最適音速を判定することを特徴とする請求項１に記載の超音波診断装置。
【請求項４】
前記最適音速判定手段は、前記微小構造物と判定されたＲＦ信号の、前記仮定音速に依る位相変化から前記最適音速を判定することを特徴とする請求項１に記載の超音波診断装置。
【請求項５】
前記最適音速判定手段は、前記微小構造物と判定されたＲＦ信号の、前記仮定音速による振幅変化から前記最適音速を判定することを特徴とする請求項１に記載の超音波診断装置。
【請求項６】
前記最適音速判定手段は、１回の送信から複数の前記仮定音速を変化させて生成した信号を利用することを特徴とする請求項１〜５のいずれかに記載の超音波診断装置。
【請求項７】
前記最適音速判定手段は、複数フレームを利用することを特徴とする請求項１〜６のいずれかに記載の超音波診断装置。
【請求項８】
前記複数フレームは、１回の送信でスキャン方向に２音線以上のＲＦデータが生成可能な装置で得られることを特徴とする請求項７に記載の超音波診断装置。
【請求項９】
前記最適音速判定手段は、前記超音波プローブの素子の配列方向に位相情報の分解能が該素子間隔以上のデータを利用することを特徴とする請求項１〜８のいずれかに記載の超音波診断装置。
【請求項１０】
前記最適音速判定手段は、複数の各前記微小構造物と判定されたＲＦ信号毎に最適音速を得ることを特徴とする請求項１〜９のいずれかに記載の超音波診断装置。
【請求項１１】
請求項１〜１０のいずれかに記載の超音波診断装置であって、さらに、前記複数の各前記微小構造物と判定されたＲＦ信号毎に得られた複数の最適音速により生成された画像を単独で、又は複数表示する表示手段を有することを特徴とする超音波診断装置。
【請求項１２】
前記表示手段は、前記複数の最適音速により生成された画像を合成した画像を表示することを特徴とする請求項１１に記載の超音波診断装置。
【請求項１３】
請求項１１または１２に記載の超音波診断装置であって、さらに、前記表示手段の表示モードを通常表示モードと、複数の画像を重ねて又は並べて表示し、あるいは単独で又は複数表示する表示モードとを切り替えるモード切替手段を有することを特徴とする超音波診断装置。
【請求項１４】
被検者に向けて超音波を送信すると共に、被検者から反射された超音波信号を受信することにより受信信号を出力する複数の素子が配列された超音波プローブと、
前記被検者に向けて送信する超音波の実音速に対して予め設定された仮定音速を変化させる手段と、
前記仮定音速を変化させた時、所定の対象領域における前記受信信号から前記仮定音速に基づく遅延でフォーカスして得られるＲＦ信号の前記素子の配列方向の位相変化から、前記被検者の超音波音速である最適音速を判定する最適音速判定手段と、
を有することを特徴とする超音波診断装置。
【請求項１５】
前記最適音速判定手段は、微小構造物における仮定音速に依る位相凹凸変化から最適音速を判定することを特徴とする請求項１４に記載の超音波診断装置。
【請求項１６】
前記最適音速判定手段は、１回の送信から複数の前記仮定音速を変化させて生成した信号を利用することを特徴とする請求項１４または１５に記載の超音波診断装置。
【請求項１７】
前記最適音速判定手段は、前記超音波プローブの素子の配列方向に位相情報の分解能が該素子間隔以上のデータを利用することを特徴とする請求項１４〜１６のいずれかに記載の超音波診断装置。
【請求項１８】
前記最適音速判定手段は、複数の各対象領域毎に最適音速を得ることを特徴とする請求項１４〜１７のいずれかに記載の超音波診断装置。
【請求項１９】
請求項１４〜１８のいずれかに記載の超音波診断装置であって、さらに、複数の前記各対象領域毎に得られた複数の最適音速により生成された画像を単独で、又は複数表示する表示手段を有することを特徴とする超音波診断装置。
【請求項２０】
前記表示手段は、前記複数の最適音速により生成された画像を合成した画像を表示することを特徴とする請求項１９に記載の超音波診断装置。
【請求項２１】
請求項１９または２０に記載の超音波診断装置であって、さらに、前記表示手段の表示モードを通常表示モードと、複数の画像を重ねて又は並べて表示し、あるいは単独で又は複数表示する表示モードとを切り替えるモード切替手段を有することを特徴とする超音波診断装置。
【請求項２２】
複数の素子が配列された超音波プローブから被検者に向けて超音波を送信すると共に、被検者から反射された超音波信号を受信し、
前記被検者に向けて送信する超音波の実音速に対して予め設定された仮定音速を変化させたときに、前記受信した信号から前記仮定音速に基づく遅延でフォーカスして得られるＲＦ信号により微小構造物の判定を行い、該微小構造物と判定されたＲＦ信号の位相情報から、被検者の超音波音速である最適音速を判定することを特徴とする超音波診断方法。
【請求項２３】
前記最適音速の判定は、前記微小構造物と判定されたＲＦ信号の振幅情報から前記最適音速を判定することによって行うことを特徴とする請求項２２に記載の超音波診断方法。
【請求項２４】
前記最適音速の判定は、前記微小構造物と判定されたＲＦ信号の、前記超音波プローブの素子の配列方向の位相変化から前記最適音速を判定することによって行うことを特徴とする請求項２２に記載の超音波診断方法。
【請求項２５】
前記最適音速の判定は、前記微小構造物と判定されたＲＦ信号の、前記仮定音速に依る位相変化から前記最適音速を判定することによって行うことを特徴とする請求項２２に記載の超音波診断方法。
【請求項２６】
前記最適音速の判定は、前記微小構造物と判定されたＲＦ信号の、前記仮定音速による振幅変化から前記最適音速を判定することによって行うことを特徴とする請求項２２に記載の超音波診断方法。
【請求項２７】
複数の素子が配列された超音波プローブから被検者に向けて超音波を送信すると共に、被検者から反射された超音波信号を受信し、
前記被検者に向けて送信する超音波の実音速に対して予め設定された仮定音速を変化させたときに、所定の対象領域から受信した信号から前記仮定音速に基づく遅延でフォーカスして得られるＲＦ信号の前記素子の配列方向の位相変化から、前記被検者の超音波音速である最適音速を判定することを特徴とする超音波診断方法。
【請求項２８】
前記所定の対象領域として特に微小構造物における仮定音速に依る位相凹凸変化から最適音速を判定することを特徴とする請求項２７に記載の超音波診断方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波診断方法及び装置に係り、特に、超音波を用いて被検者の断層画像を撮影し表示する超音波診断方法及び装置に関する。
【背景技術】
【０００２】
従来より、超音波を用いて被検者の断層画像を取得して医療診断に供することが行われているが、超音波診断装置では複数の配列する探触子素子からの複数の受信信号に遅延時間分布を与えて、所定の方向に指向性を持つ超音波ビームを形成することにより被検者の断層画像を取得するようにしている。
【０００３】
このとき対象となる診断部位の代表的な音速を想定して（これを設定音速、あるいは仮定音速という）遅延時間を設定している。しかし、被検者の生体音速（実音速、あるいは最適音速という）は一様ではなく組織によって異なり、仮定音速と最適音速が異なると画質が劣化するという問題があった。
【０００４】
そこで、従来より、最適音速を設定する様々な方法が知られている。例えば、装置に設定された超音波音速値を補正可能としフォーカスを向上することができる超音波断層装置が知られている（例えば、特許文献１等参照）。
【０００５】
これは、操作者が入力した任意の超音波音速値に基づいてフォーカスが計算され、そのフォーカスで超音波の送受信が行われ、操作者が入力する超音波音速値を変えると、それに対応してフォーカスが変わり超音波画像の画質が変化するので、操作者が画像を見ながら最もフォーカスの合う超音波音速値を選ぶようにするものである。
【０００６】
またこれは、超音波受信信号の振幅が最大となる場合の超音波音速値を求めることにより、あるいは超音波受信信号のビーム幅が最小となる場合の超音波音速値を求めることにより、さらには超音波受信信号の振幅について空間周波数の高周波成分又は分散が最大となる場合の超音波音速値を求めることによって、設定された超音波音速値を補正することでフォーカスを向上させるようにするものである。
【特許文献１】特開平８−３１７９２６号公報
【発明の開示】
【発明が解決しようとする課題】
【０００７】
しかしながら、上記従来技術では、何種類もの仮定音速（設定音速）の画像を作成する必要があり、その分のメモリや回路、また処理時間が必要となるという問題がある。
【０００８】
本発明は、このような事情に鑑みてなされたもので、メモリや回路、また処理時間を要することなく最適音速を求めることのできる超音波診断方法及び装置を提供することを目的とする。
【課題を解決するための手段】
【０００９】
前記目的を達成するために、請求項１に記載の発明は、被検者に向けて超音波を送信すると共に、被検者から反射された超音波信号を受信することにより受信信号を出力する複数の素子が配列された超音波プローブと、前記被検者に向けて送信する超音波の実音速に対して予め設定された仮定音速を変化させる手段と、前記仮定音速を変化させ、前記受信信号から前記仮定音速に基づく遅延でフォーカスして得られるＲＦ信号により微小構造物の判定を行い、該微小構造物と判定されたＲＦ信号の位相情報から、被検者の超音波音速である最適音速を判定する最適音速判定手段と、を有することを特徴とする超音波診断装置を提供する。
【００１０】
これにより、微小構造物信号の仮定音速に依る位相変化特性を利用して最適音速を得ることができる。
【００１１】
また、請求項２に示すように、前記最適音速判定手段は、前記微小構造物と判定されたＲＦ信号の振幅情報から前記最適音速を判定することを特徴とする。
【００１２】
これにより、微小構造物信号の仮定音速に依る振幅変化特性を利用して最適音速を得ることができる。
【００１３】
また、請求項３に示すように、前記最適音速判定手段は、前記微小構造物と判定されたＲＦ信号の、前記超音波プローブの素子の配列方向の位相変化から前記最適音速を判定することを特徴とする。
【００１４】
また、請求項４に示すように、前記最適音速判定手段は、前記微小構造物と判定されたＲＦ信号の、前記仮定音速に依る位相変化から前記最適音速を判定することを特徴とする。
【００１５】
また、請求項５に示すように、前記最適音速判定手段は、前記微小構造物と判定されたＲＦ信号の、前記仮定音速による振幅変化から前記最適音速を判定することを特徴とする。
【００１６】
このように、いろいろな方法で最適音速を判定することが可能となる。
【００１７】
また、請求項６に示すように、前記最適音速判定手段は、１回の送信から複数の前記仮定音速を変化させて生成した信号を利用することを特徴とする。
【００１８】
また、請求項７に示すように、前記最適音速判定手段は、複数フレームを利用することを特徴とする。
【００１９】
これにより、各クレーム間でずれがなくなり、高フレームレートでの処理が可能となった。
【００２０】
また、請求項８に示すように、前記複数フレームは、１回の送信でスキャン方向に２音線以上のＲＦデータが生成可能な装置で得られることを特徴とする。
【００２１】
また、請求項９に示すように、前記最適音速判定手段は、前記超音波プローブの素子の配列方向に位相情報の分解能が該素子間隔以上のデータを利用することを特徴とする。
【００２２】
これにより、高ＳＮの判定が可能となる。
【００２３】
また、請求項１０に示すように、前記最適音速判定手段は、複数の各前記微小構造物と判定されたＲＦ信号毎に最適音速を得ることを特徴とする。
【００２４】
また、請求項１１に示すように、請求項１〜１０のいずれかに記載の超音波診断装置であって、さらに、前記複数の各前記微小構造物と判定されたＲＦ信号毎に得られた複数の最適音速により生成された画像を単独で、又は複数表示する表示手段を有することを特徴とする。
【００２５】
また、請求項１２に示すように、前記表示手段は、前記複数の最適音速により生成された画像を合成した画像を表示することを特徴とする。
【００２６】
また、請求項１３に示すように、請求項１１または１２に記載の超音波診断装置であって、さらに、前記表示手段の表示モードを通常表示モードと、複数の画像を重ねて又は並べて表示し、あるいは単独で又は複数表示する表示モードとを切り替えるモード切替手段を有することを特徴とする。
【００２７】
また、同様に前記目的を達成するために、請求項１４に記載の発明は、被検者に向けて超音波を送信すると共に、被検者から反射された超音波信号を受信することにより受信信号を出力する複数の素子が配列された超音波プローブと、前記被検者に向けて送信する超音波の実音速に対して予め設定された仮定音速を変化させる手段と、前記仮定音速を変化させた時、所定の対象領域における前記受信信号から前記仮定音速に基づく遅延でフォーカスして得られるＲＦ信号の前記素子の配列方向の位相変化から、前記被検者の超音波音速である最適音速を判定する最適音速判定手段と、を有することを特徴とする超音波診断装置を提供する。
【００２８】
このように、微小構造物に限定することなく所定の対象領域を設定することによりそこの位相変化から最適音速を得ることができる。
【００２９】
また、請求項１５に示すように、前記最適音速判定手段は、微小構造物における仮定音速に依る位相凹凸変化から最適音速を判定することを特徴とする。
【００３０】
このように、所定の対象領域として特に微小構造物における位相凹凸変化から最適音速得ることができる。
【００３１】
また、請求項１６に示すように、前記最適音速判定手段は、１回の送信から複数の前記仮定音速を変化させて生成した信号を利用することを特徴とする。
【００３２】
また、請求項１７に示すように、前記最適音速判定手段は、前記超音波プローブの素子の配列方向に位相情報の分解能が該素子間隔以上のデータを利用することを特徴とする。
【００３３】
これにより、高ＳＮでの判定が可能となる。
【００３４】
また、請求項１８に示すように、前記最適音速判定手段は、複数の各対象領域毎に最適音速を得ることを特徴とする。
【００３５】
また、請求項１９に示すように、請求項１４〜１８のいずれかに記載の超音波診断装置であって、さらに、複数の前記各対象領域毎に得られた複数の最適音速により生成された画像を単独で、又は複数表示する表示手段を有することを特徴とする。
【００３６】
また、請求項２０に示すように、前記表示手段は、前記複数の最適音速により生成された画像を合成した画像を表示することを特徴とする。
【００３７】
また、請求項２１に示すように、請求項１９または２０に記載の超音波診断装置であって、さらに、前記表示手段の表示モードを通常表示モードと、複数の画像を重ねて又は並べて表示し、あるいは単独で又は複数表示する表示モードとを切り替えるモード切替手段を有することを特徴とする。
【００３８】
また、同様に前記目的を達成するために、請求項２２に記載の発明は、複数の素子が配列された超音波プローブから被検者に向けて超音波を送信すると共に、被検者から反射された超音波信号を受信し、前記被検者に向けて送信する超音波の実音速に対して予め設定された仮定音速を変化させたときに、前記受信した信号から前記仮定音速に基づく遅延でフォーカスして得られるＲＦ信号により微小構造物の判定を行い、該微小構造物と判定されたＲＦ信号の位相情報から、被検者の超音波音速である最適音速を判定することを特徴とする超音波診断方法を提供する。
【００３９】
これにより、微小構造物信号の仮定音速に依る位相変化特性を利用して最適音速を得ることができる。
【００４０】
また、請求項２３に示すように、前記最適音速の判定は、前記微小構造物と判定されたＲＦ信号の振幅情報から前記最適音速を判定することによって行うことを特徴とする。
【００４１】
これにより、微小構造物信号の仮定音速に依る振幅変化特性を利用して最適音速を得ることができる。
【００４２】
また、請求項２４に示すように、前記最適音速の判定は、前記微小構造物と判定されたＲＦ信号の、前記超音波プローブの素子の配列方向の位相変化から前記最適音速を判定することによって行うことを特徴とする。
【００４３】
また、請求項２５に示すように、前記最適音速の判定は、前記微小構造物と判定されたＲＦ信号の、前記仮定音速に依る位相変化から前記最適音速を判定することによって行うことを特徴とする。
【００４４】
また、請求項２６に示すように、前記最適音速の判定は、前記微小構造物と判定されたＲＦ信号の、前記仮定音速による振幅変化から前記最適音速を判定することによって行うことを特徴とする。
【００４５】
また、同様に前記目的を達成するために、請求項２７に記載の発明は、複数の素子が配列された超音波プローブから被検者に向けて超音波を送信すると共に、被検者から反射された超音波信号を受信し、前記被検者に向けて送信する超音波の実音速に対して予め設定された仮定音速を変化させたときに、所定の対象領域から受信した信号から前記仮定音速に基づく遅延でフォーカスして得られるＲＦ信号の前記素子の配列方向の位相変化から、前記被検者の超音波音速である最適音速を判定することを特徴とする超音波診断方法を提供する。
【００４６】
また、請求項２８に示すように、前記所定の対象領域として特に微小構造物における仮定音速に依る位相凹凸変化から最適音速を判定することを特徴とする。
【発明の効果】
【００４７】
以上説明したように、本発明によれば、微小構造物信号の仮定音速に依る位相変化特性を利用して最適音速を得ることができる。
【発明を実施するための最良の形態】
【００４８】
以下、添付図面を参照して、本発明に係る超音波診断方法及び装置について詳細に説明する。
【００４９】
本発明は、受信した超音波画像から判定画像を生成する際に設定する超音波音速（仮定音速）を変化させた時の信号の位相変化特性を利用して微小構造物を判定し、微小構造物と判定された信号の仮定音速に依る位相変化特性を利用して被検者の超音波音速である最適音速を求めるものである。
【００５０】
ここで、微小構造物とは、具体的には、例えば、被検者内の微小石灰化、気泡、あるいは造影剤などを言う。
【００５１】
なお、上記設定する超音波音速を所定量ずつステップ刻みで複数変化させて判定画像を生成するに当たり、被検者に送信する実際の超音波音速（実音速あるいは最適音速）に対して、複数変化させて設定する超音波音速を以下の実施形態では設定音速あるいは仮定音速と言う。
【００５２】
そして、本発明は、画像生成時の仮定音速を種々に変更して得られる信号の位相変化特性を利用して、最適音速を求めるものである。
【００５３】
具体的に、その位相変化特性とは、後で詳しく説明するが、簡単に言うと、以下のようなものである。
【００５４】
すなわち、仮定音速が最適音速と異なる場合、微小構造物信号はスキャン方向（超音波プローブの素子（振動子）の配列方向）に特徴的な位相凹凸変化を示し、仮定音速を変化させた場合の、微小構造物信号、連続面・線信号及びスペックルの振幅・位相変化特性は異なっている。
【００５５】
まず、微小構造物については、仮定音速が最適音速（実音速）より小さい（遅い）場合、スキャン方向に位相は凸型に変化し、その傾きは最適音速に近い程急峻になり、また仮定音速が最適音速より大きい（速い）場合、スキャン方向に位相は凹型に変化し、その傾きは最適音速に近い程急峻になる。なお、振幅については、最適音速に近い程、大きく、また形状は急峻になる。
【００５６】
次に、連続面・線の場合には、仮定音速によらず位相は一様であり、振幅については最適音速に近い程、大きくなる。
【００５７】
また、スペックルの場合には、仮定音速に依って振幅も位相もランダムに変化する。
【００５８】
図１は、本発明に係る超音波診断装置の一実施形態の概略構成を示すシステム構成図である。
【００５９】
図１に示すように、超音波診断装置１は、超音波を用いて被検者の診断部位について超音波画像を撮影して表示するものであり、超音波プローブ１０、送受信部１２、走査制御部１４、ＡＤ変換部１６、画像生成部１８、最適音速判定部２０、表示画像生成部２２、モニタ２４及びモード切替手段２６を有して構成されている。
【００６０】
超音波プローブ１０は、被検者の体内の診断部位に向けて超音波を送信するとともに体内で反射してきた超音波を受信するものである。本実施形態の超音波プローブ１０は、１次元の超音波トランスデューサアレイを構成する複数の超音波トランスデューサを備えており、各超音波トランスデューサは、例えばＰＺＴ等の圧電素子の両端に電極を形成した振動子によって構成されている。この電極は信号線によって送受信部１２と接続されている。各電極に電圧を印加すると振動子は超音波を発生する。また、振動子は反射してきた超音波を受信すると電気信号を発生し、これが受信信号として出力される。
【００６１】
送受信部１２は、超音波プローブ１０に超音波送信信号を与え振動子から超音波を発生させ、走査制御部１４から与えられた遅延に基づいて送波する。そして、反射した超音波を受信して超音波プローブ１０が出力した各素子の受信信号をそのまま（受波フォーカスをかけず）増幅する。
【００６２】
ＡＤ変換部１６は、送受信部１２から超音波受信信号を受け取りＡＤ変換して画像生成部１８に引き渡す。画像生成部１８はＡＤ変換部１６から受け取った受信データを保存する。画像生成部１８では、保存された各素子の受信データから、詳しくは後述するが、様々に設定される音速（これを上述したように被検者に送波する実際の音速（実音速）に対して仮定音速という。）に基づく遅延で受波フォーカスされ、各仮定音速に基づくＲＦデータが生成される。
【００６３】
最適画像判定部２０は、画像生成部１８で生成されたＲＦデータから微小構造物のスキャン方向の位相変化を見て最適音速を判定するものである。
【００６４】
また、表示画像生成部２２は、画像生成部１８で生成された画像と、最適音速判定部２０で生成された判定画像による判定結果からモニタ２４に表示するための表示画像を生成するものである。モード切替手段２６は、モニタ２４への画像の表示モードを切り替えるものである。
【００６５】
次に、上記超音波診断装置１の作用を説明する前に、仮定音速を変化させた時の位相変化特性について説明する。
図２〜９に、仮定音速を変化させた時の位相変化特性を表したグラフを示す。
【００６６】
各グラフは、それぞれ仮定音速を大体１４００［ｍ／ｓ］から１６２０［ｍ／ｓ］まで４０［ｍ／ｓ］あるいは２０［ｍ／ｓ］刻みで変化させたときの位相変化特性を、横軸をスキャン方向（Ｘ位置）、縦軸を位相として表示したものである。
【００６７】
図２は、仮定音速１４００［ｍ／ｓ］〜１５００［ｍ／ｓ］における微小構造物信号の仮定音速に依る位相変化特性を示すグラフであり、図３は、仮定音速１５００［ｍ／ｓ］〜１６２０［ｍ／ｓ］における微小構造物信号の仮定音速に依る位相変化特性を示すグラフである。
【００６８】
仮定音速１４００［ｍ／ｓ］〜１５００［ｍ／ｓ］の図２の場合には、Ｘ位置１００〜１２０付近において、仮定音速１５００［ｍ／ｓ］のグラフは正の傾きを有し、その他の、仮定音速が１５００［ｍ／ｓ］より小の（すなわち仮定音速がより遅い）グラフはいずれも右下がりで、仮定音速が１５００［ｍ／ｓ］に近い程その傾きが急峻であり、仮定音速が１５００［ｍ／ｓ］より遅くなるにつれて傾きが緩やかになっている。
【００６９】
また、仮定音速１５００［ｍ／ｓ］〜１６２０［ｍ／ｓ］の図３の場合には、Ｘ位置１００〜１２０付近において、いずれも右上がりのグラフとなっている。そして、仮定音速が１５００［ｍ／ｓ］のときが最も傾きが大きく、仮定音速が１５００［ｍ／ｓ］より大きくなる程傾きが緩やかになっている。
【００７０】
図２及び図３のこのようなグラフの形状は、Ｘ位置１００〜１２０付近に微小構造物が存在していることを示すものであると考えられる。
【００７１】
図４は、仮定音速１４００［ｍ／ｓ］〜１４８０［ｍ／ｓ］における面信号の仮定音速に依る位相変化特性を示すグラフであり、図５は、仮定音速１５４０［ｍ／ｓ］〜１６２０［ｍ／ｓ］における面信号の仮定音速に依る位相変化特性を示すグラフである。
【００７２】
図４及び図５からわかるように、いずれもＸ位置１００〜１３０及び１５０〜１８０付近において、仮定音速を変えても位相があまり変化していない。これはその部分に面（連続面）が存在することを示すものであると考えられる。
【００７３】
図６は、仮定音速１４００［ｍ／ｓ］〜１４８０［ｍ／ｓ］におけるスペックルの仮定音速に依る位相変化特性を示すグラフであり、図７は、仮定音速１５４０［ｍ／ｓ］〜１６２０［ｍ／ｓ］におけるスペックルの仮定音速に依る位相変化特性を示すグラフである。図６及び図７に示すように、スペックルの場合は仮定音速を変えると位相はランダムに変化する。
【００７４】
また、図８は、仮定音速１４００［ｍ／ｓ］〜１５００［ｍ／ｓ］における微小構造物の仮定音速に依る振幅変化特性を示すグラフであり、図９は、仮定音速１５００［ｍ／ｓ］〜１６２０［ｍ／ｓ］における微小構造物信号の仮定音速に依る振幅変化特性を示すグラフである。
【００７５】
図８及び図９からわかるように、微小構造物の振幅の変化を示すグラフは、いずれもＸ位置１１０付近に頂上（最大値）を有する山型（上に凸）のグラフで、仮定音速１５００［ｍ／ｓ］のとき振幅値が最大で、仮定音速が１５００［ｍ／ｓ］に近づく程最大振幅値が大きくなるとともに、形状も急峻となっている。
【００７６】
次に、仮定音速を変化させたときの微小構造物の位相変化が図２、図３にグラフで示したような特性を有する理由を説明する。
【００７７】
図１０において、点Ａ（０，ｚ０ )から音速Ｖ０で反射した超音波を、超音波プローブ１０における位置Ｘの素子（振動子）で、反射後の時刻ｔにおいて観測したとする。すると、この時刻ｔは次の（１）式のように算出される。
【００７８】
ｔ＝sqrt(ｚ０２ +Ｘ２）／Ｖ０            ・・・（１）
なお、式（１）においてsqrt(  ）は（  ）内の値の平方根をとることを意味する。
【００７９】
また、図１０において、点Ａ’（ｘ，ｚ）から音速Ｖで反射した超音波を、同じく超音波プローブ１０における位置Ｘの素子（振動子）で、反射後の時刻ｔにおいて観測したとする。上と同様にこの時刻ｔは、次の式（２）で表される。
【００８０】
ｔ＝sqrt｛ｚ２ +（Ｘ−ｘ）２｝／Ｖ       ・・・（２）
それぞれ式（１）と式（２）で与えられる曲線が（Ｘ，ｔ）平面において接するときの点Ａ’の軌跡は次の式（３）で与えられる。
【００８１】
ｚ２＝ｘ２×｛Ｖ２／（Ｖ０２−Ｖ２）｝＋ｚ０２Ｖ２／Ｖ０２      ・・・（３）
点Ａ’は、最適音速（実音速）に対して仮定音速Ｖとして位相を整合して加算した場合に信号が強くなる位置を示している。
【００８２】
式（３）より点Ａ’（ｘ，ｚ）の軌跡は、Ｖ＞Ｖ０のときは、原点を中心とした楕円となり、Ｖ＜Ｖ０のときは、原点を中心とした双曲線となる。なお、Ｂモード画像においてはｚ軸の下方向を上としているので、Ｖ＞Ｖ０の場合の原点を中心とした楕円の軌跡を凹型、Ｖ＜Ｖ０の場合の原点を中心とした双曲線の軌跡を凸型とする。
【００８３】
なお、図１１に、式（１）と式（２）で与えられる曲線が（Ｘ，ｔ）平面において接している様子を示す。図１１（ａ）はＶ＜Ｖ０の場合であり、実線Ｊは点Ａからの反射波を、破線Ｈ１は図１０のようにｚ軸より右側にある場合の点Ａ’からの反射波を表している。また、破線Ｈ２は図示は省略するが図１０のｚ軸より左側にある点からの反射波を表している。いまＶ＜Ｖ０の場合であるので、同じＸの位置に対しては点Ａ’からの反射波の方が時刻ｔが大きいため、破線Ｈ１（Ｈ２）が実線Ｊより上側に表れている。
【００８４】
また、図１１（ｂ）は、Ｖ＝Ｖ０の場合でり、図１１（ｃ）はＶ＞Ｖ０の場合である。ＶがＶ０に近づく場合には破線Ｈは実線Ｊに近づき、Ｖ＝Ｖ０の場合には、破線Ｈは実線Ｊと一致する。図１１（ｃ）の場合は、図１１（ａ）とは逆で破線が実線よりも下側に表れる。
【００８５】
これらの図から、上記破線が上記実線に接するような点Ａ’（ｘ，ｚ）の軌跡は上で述べたような傾向を有することが直感的に理解できる。
【００８６】
なお、ここで説明したモデルにおいては、観測された反射波を単純に線としたが、実際にはｔ方向に幅を持った波形である事や、Ｘ方向に強度差がある事なども考慮する必要がある。また、本モデルでは簡単のために、点Ａと点Ａ’からの反射を同時としたが、実際には、それぞれの点に超音波を送波してから反射する迄の時間も考慮する必要がある。
【００８７】
次に、図１の装置構成における画像生成部１８の作用を図１２のフローチャートに沿って説明する。
【００８８】
画像生成部１８は、仮定音速を変化させていろいろな音速で得られたデータから画像を生成するものである。
【００８９】
まず図１２のステップＳ１００において、いろいろ変化させる仮定音速の初期値を設定する。この値は特に限定されるものではなく、適宜決めればよい。例えば、前述した図２等の例のように、１４００［ｍ／ｓ］のように決めればよい
そして設定された初期値により、走査制御部１４によって制御された送受信部１２から超音波プローブ１０に信号が送られ、その仮定音速初期値によるデータが取得され画像生成部１８に送られる。
【００９０】
次にステップＳ１１０において、仮定音速を所定量１ステップ変更し、変更された仮定音速による超音波データが取得される。この１ステップの所定量は、特に限定されず、例えば図２等の例のように４０［ｍ／ｓ］でもよいし、１０［ｍ／ｓ］でも、２０［ｍ／ｓ］でもよく、所定量だけ仮定音速を変化させていく。
【００９１】
次にステップＳ１２０において、得られた各仮定音速によるデータを位相を整合して加算し、ＲＦ（Radio Frequency)データを生成する。このＲＦデータは、振幅情報と位相情報の両方を含むものである。このようにすべての仮定音速での画像でＲＦデータを作成する。
【００９２】
そしてステップＳ１３０において、画像生成が終了したか否か判断し、まだ終了していない場合にはステップＳ１１０に戻り、また仮定音速を１ステップ変更し画像生成を続行する。画像生成の終了は、すべての仮定音速についての処理が終了したか否かで判断する。それは例えば、仮定音速を何ステップ変更したら終了するかを予め決めておき、その回数をカウントして判断するようにすればよい。
【００９３】
次に、最適音速判定部２０の作用を説明する。
【００９４】
図１３は、最適音速判定部２０における最適音速判定処理の流れを示すフローチャートである。
【００９５】
まず図１３のステップＳ２００において、仮定音速の初期値を設定する。ただ、これはすでに画像生成部１８における処理で得られているデータを用いればよい。次に、ステップＳ２０２において、微小構造物判定画像及び音速判定画像の各画素の値を０として初期化する。
【００９６】
次に、ステップＳ２０４において、仮定音速を１ステップ変更し、その音速でのデータを取得するが、これも画像生成部１８における処理で得られているデータを使用すればよい。
【００９７】
次に、ステップＳ２０６において、その仮定音速でのデータから位相スキャン方向の２次微分値を算出する。スキャン方向とは、超音波プローブ１０の振動子（素子）の配列方向と一致している。
【００９８】
次に、ステップＳ２０８において、所定サイズのカーネルで上記２次微分値を積分し、積分値を算出する。カーネルのサイズは、特に限定はされず、解像度に応じて、９×４や１６×８等のものが用いられる。
【００９９】
次に微小構造物判定画像及び音速判定画像の生成が並行して同時に行われる。
【０１００】
ステップＳ２１０において、上で算出した積分値の絶対値をとり、次のステップＳ２１２において、この積分値の絶対値を微小構造物判定画像に加算して行く。
【０１０１】
このようにして、全ての仮定音速について２次微分値を積分した値を当初０に初期設定されていた微小構造物判定画像に足し合わせて行くことにより微小構造物判定画像が作成される。微小構造物の場合は、仮定音速が最適音速よりも速いときは２次微分値は正となり、仮定音速が最適音速よりも遅いときは２次微分値は負となるので、各仮定音速による２次微分値を所定のカーネルで積分すると微小構造物のところだけ信号が強くでる。
【０１０２】
従って、これらを加算して生成された微小構造物判定画像は、微小構造物のところだけ信号が強く出た画像となっており、これにより微小構造物の判定が可能となる。
【０１０３】
また、一方ステップＳ２１４において、指定された仮定音速用の音速判定画像に、所定サイズでのカーネルによる積分値を格納する。
【０１０４】
次のステップＳ２１６において、全ての仮定音速について上記処理が終了したか否か判断し、まだ終了していなければステップＳ２０４に戻り、仮定音速を１ステップ変更して次の仮定音速についての処理を繰り返す。また、全ての仮定音速についての処理が終了した場合には、次のステップＳ２１８において、微小構造物判定画像の所定以上の値となる位置を微小構造物位置として記録する。
【０１０５】
最後にステップＳ２２０において、各微小構造物位置について各仮定音速の音速判定画像の値の変化から最適音速を判定し、記録する。
【０１０６】
次に、最適音速の判定方法を詳しく説明する。
【０１０７】
最適音速の判定に利用する仮定音速は１種類以上である。また、図１３のフローチャートでは、ステップＳ２０４からステップＳ２１４において、同じ仮定音速の下で微小構造物判定画像の処理と音速判定画像の処理を行っているように示されているが、微小構造物判定に利用する仮定音速と、最適音速判定に利用する仮定音速の種類や数は、必ずしも一致しなくても良い。
【０１０８】
また、微小構造物の判定方法は、図１３のフローチャートのものには限定されず、他の方法でも良い。例えば、積分値の絶対値をとらずに、仮定音速が最適音速より遅い場合には符号を反転する方法でも良い。本発明はそもそも最適音速を求めるものであり、まだ最適音速はわからないのであるが、所定以上遅い仮定音速及び所定以上速い仮定音速を利用することでこの方法を利用することが可能となる。これにより、最適音速がわからなくてもそこから判定画像を生成することができる。
【０１０９】
その他、符号のみ積分する方法、２次微分値の分散や傾き等、一様性を数値化し指標にする方法、あるいは位相ずれに応じた凹凸形状フィルタを掛け、位相又は波形画像に対して相互相関をとって抽出する方法などを利用することもできる。
【０１１０】
また、微小構造物判定は、ユーザが指定するようにしても良い。
【０１１１】
図１４に、最適音速判定において、仮定音速と微小構造物位置における音速判定画像の値の関係を示す。
【０１１２】
図１４に示すように、仮定音速が最適音速より遅いとすると音速判定画像の値は負の値として出て、仮定音速が最適音速より速いとすると音速判定画像の値は正の値として出る。さらに、いずれの場合も仮定音速が最適音速に近づく程、その傾きが急になるのでそれを積分した結果としての音速判定画像の値は、図に矢印で示したように、負の場合にはより小さな値となり、正の場合にはより大きな値となる。このように、音速判定画像の値は仮定音速が最適音速より速いか遅いか、又ずれ量によって略一意的に求まる。
【０１１３】
従って、最低１種類の仮定音速における値から最適音速が判定できる。その最適音速からのずれ量とその音速判定画像の値の関係は、予め実験的に求めてテーブルに持っても良いし、式で与えるようにしても良い。
【０１１４】
また、２種類以上の仮定音速を利用する場合には、各仮定音速における値から求められる最適音速値の重み付け平均をとる方法でも良いし、最適音速値を仮定し、各仮定音速につきテーブルや式から求めた値と、実際の値の重み付け二乗誤差が最小となる最適音速値を採用する方法でも良い。
【０１１５】
仮定音速と判定画像の値の関係が、テーブルや式から求められる値から所定以上ばらついてしまう位置を微小構造物とは見做さずに外すようにしても良い。
【０１１６】
最適音速を判定する他の方法として、２次微分値のカーネルでの積分値の代わりに、最大値、最小値を利用する方法や、各仮定音速での位相や波形画像と予め用意した凹凸形状フィルタとの相関をとり最適凹凸判定する方法でも良い。
【０１１７】
また、図８、９に示したような振幅変化特性をも合わせて利用することで最適音速判定の精度を向上させることができる。
【０１１８】
図１５に、最適音速判定において、仮定音速と振幅を利用した音速判定画像の値の関係を示す。
【０１１９】
図１５に示すように、振幅情報を利用した判定画像値というのは、最適音速値をはさんで左右対称となる。この場合、前述したように、最適音速に近い程凸形状が急峻になるという特性から、上記のような位相利用と同様の方法をとることができる。ただし、振幅の凸形状変化は、最適音速より速くても遅くても同じであるため、単独では速いか遅いかの判定ができない。そこで、最適音速判定に複数フレームを利用するようにしてもよい。
【０１２０】
最近のソフトウエアベースの超音波装置やアナログベースでも高性能な回路構成によって、同じ１回の送信から得られた受信信号から種々の仮定音速で画像生成することが可能となってきている。本装置構成においては、種々の仮定音速での画像をフレームずれ無しに得られるため、その結果、高精度に最適音速を判定することが可能となる。
【０１２１】
また、スキャン方向に高分解能な位相情報を利用することにより、より高精度に最適音速を判定することが可能となる。また、±１８０°の反転も起き難い。また、複数フレームを利用することで、より高精度に最適音速を判定することができる。また、同じ１回の送信から得られた受信データを利用することで、より高精度に最適音速の判定が可能である。
【０１２２】
また、ここでは微小構造物の位相変化特性をスキャン方向の凹凸変化で表現しているが、同じ位置における仮定音速に依る位相変化としても表現することができ、最適音速判定もどちらの特性として利用しても良い。
【０１２３】
図１６は、表示画像生成部２２における処理の流れを示すフローチャートである。
【０１２４】
まず、図１６のステップＳ３００において、各微小構造物位置と各位置における最適音速値を取得する。
【０１２５】
次に、ステップＳ３１０において、各位置における各最適音速値に対応する振幅画像を取得する。
【０１２６】
次に、ステップＳ３２０において、各最適音速の振幅画像を各微小構造物位置を中心として合成する。合成方法は、各画像を各微小構造物位置を中心として他画像とオーバーラップするように長方形、円、あるいはユーザの指定や、輝度レベルで決められる任意形状の領域を設定し、オーバーラップする領域では各微小構造物位置からの距離に応じた割合で振幅画像を足し合わせるようにする。
【０１２７】
なお、ここでは複数の最適音速画像を合成する例を挙げたが、単に複数画像を並べても良いし、平均の音速の画像を選択しても良い。また、複数の微小構造物のうち、利用するものを限定するようにしても良い。
【０１２８】
次に、ステップＳ３３０において、その結果を対数圧縮し、ゲイン／ＤＲ（ダイナミックレンジ）／ＳＴＣ（深さ重み付け）／グレーマップ調整し、さらにスキャンコンバートして表示画像を生成する。
【０１２９】
なお、画像の表示モードは、複数の最適音速により生成された画像を単独で表示しても良いし、複数表示しても良い。また、複数の最適音速により生成された画像を合成した画像を表示しても良い。
【０１３０】
さらに、モード切替手段２６によって、画像の表示モードを通常表示モードとその他の表示モードとで切り替えるようにしても良い。
【０１３１】
以上説明したように、本実施形態によれば、微小構造物信号の仮定音速と位相の関係は、仮定音速が最適音速より速いか遅いか、あるいは、ずれ量によって略一意的に決まるため、最低１種類の仮定音速のデータがあれば、最適音速を求めることができ、その結果、少ないメモリや回路、また処理時間で最適音速を得ることができる。
【０１３２】
上で述べたように、本実施形態によれば、微小構造物の仮定音速による位相変化特性を利用して最適音速を得ることができる。
【０１３３】
完全なスペックルの場合は、仮定音速によって位相がランダムに変化し、連続面の場合は仮定音速によらず一定となる。一方、微小構造物の位相はスキャン方向に凹凸変化を示し、また仮定音速によって凹凸形状が変化する。生体内においては、完全なスペックルや、表面が滑らかで各部分の反射が同じ完全な連続面ではなく、局所的に強い反射を含むようなスペックルで略構成されていると考えられる。
【０１３４】
図１７にスペックルの仮定音速による位相変化を示す。図１７（ａ）は、最適音速より速い場合であり、図１７（ｂ）は最適音速より遅い場合である。図１７（ａ）に示すように、この場合には局所的に下に凸の凸形状３０（孤立点）が見られる。また図１７（ｂ）においては局所的に上に凸の凸形状３２（孤立点）が見られる。
【０１３５】
また、図１８は微小構造物における位相凹凸変化のヒストグラムをとったものであり、図１８（ａ）は仮定音速が最適音速より遅い場合、図１８（ｂ）は仮定音速が最適音速より速い場合である。またそれぞれ（１）は孤立点を、（２）は完全なスペックルを、（３）は局所的に高エコーが混在するスペックルを表している。
【０１３６】
図１８（ａ）の（１）に示すように仮定音速が最適音速よりも遅い場合には、孤立点の位相凹凸値は負の側に現れ、また図１８（ｂ）の（１）に示すように仮定音速が最適音速よりも速い場合には、孤立点の位相凹凸値は正の側に現れる。
【０１３７】
また、完全なスペックルの場合には、位相がランダムに変化するため、その位相凹凸値のヒストグラムをとると図１８（ａ）の（２）、及び図１８（ｂ）の（２）に示すように、略正規分布をなす。
【０１３８】
また、局所的に高エコーが混在するスペックルの場合には、仮定音速が最適音速よりも遅い場合は図１８（ａ）の（３）に示すように位相凹凸値のヒストグラムは負の側にその分布が偏り、また一方仮定音速が最適音速よりも速い場合は図１８（ｂ）の（３）に示すように位相凹凸値のヒストグラムは正の側に分布が偏っている。このように位相凹凸変化の分布の偏りから最適音速を求めることができる。
【０１３９】
以下、ある対象領域における仮定音速の位相凹凸変化から最適音速を判定する例について説明する。
【０１４０】
図１９に、ある対象領域における仮定音速の位相凹凸変化から最適音速を判定する手順を示す。
【０１４１】
まず図１９のステップＳ４００において、対象領域      （ＲＯＩ）を設定する。この各対象領域の設定はユーザが指定して行ってもよいし、全画面を所定サイズで分割した複数対象領域としてもよい。
【０１４２】
次にステップＳ４１０において、前述した図１２のステップＳ１００と同様に、いろいろ変化させる仮定音速の初期値を設定する。この値は特に限定されるものではなく、適宜決めればよい。例えば、図２に示した例のように、１４００［ｍ／ｓ］のように決めればよい。そして設定された初期値により、走査制御部１４によって制御された送受信部１２から超音波プローブ１０に信号が送られ、その仮定音速初期値によるデータが取得され画像生成部１８に送られる。
【０１４３】
次にステップＳ４２０において、仮定音速を所定量１ステップ変更し、変更された仮定音速による超音波データが取得される。この１ステップの所定量は、特に限定されず、例えば図２等の例のように４０［ｍ／ｓ］でもよいし、１０［ｍ／ｓ］等他の値でもよく、所定量だけ仮定音速を変化させていく。
【０１４４】
次にステップＳ４３０において、対象領域内のスキャン方向位相凹凸を数値化する。この位相凹凸の数値化も特に限定されるものではなく、例えば、２次微分値の積分値でもよいし、凹凸パターンとの相関値でもよく、あるいは凹凸に沿った波形加算でもよい。
【０１４５】
次にステップＳ４４０において、全ての仮定音速についての処理が終了したか否か判断する。まだ、全ての仮定音速についての処理が終了していない場合には、ステップＳ４２０に戻り、仮定音速を１ステップ変更しながら上記処理を繰り返す。
【０１４６】
また、全ての仮定音速についての処理が終了した場合には、次のステップＳ４５０において、各仮定音速において対象領域内の位相凹凸分布の偏りの変化から最適音速を判定する。
【０１４７】
分布の偏りの評価は、位相凹凸の最大値と最小値の平均や所定個数の平均、または絶対値の差分などを用いておこなうことができる。
【０１４８】
また、図２０に仮定音速による分布の偏りの変化を示す。図２０に示すように、偏りが最も０に近い仮定音速を最適音速と判定する。または、偏りが正負同程度の仮定音速の中間値としても良い。
【０１４９】
図２１は、前述した図１６と同様に表示画像生成部２２における処理の流れを示すフローチャートである。図２１に示すように、各対象領域の最適音速値から最適音速画像を生成しても良い。
【０１５０】
まず図２１のステップＳ５００において、各対象領域における最適音速値を取得する。次にステップＳ５１０において、各最適音速に対応する振幅画像を取得する。次にステップＳ５２０において、各最適音速の振幅画像を各対象領域位置を中心として合成する。
【０１５１】
次にステップＳ５３０において、その結果を対数圧縮し、ゲイン／ＤＲ（ダイナミックレンジ）／ＳＴＣ（深さ重み付け）／グレーマップ調整し、さらにスキャンコンバートして表示画像を生成する。
【０１５２】
また、画像の表示モードは複数の最適音速により生成された画像を単独で表示しても良いし、複数表示しても良い。また、複数の最適音速によって生成された画像を合成した画像を表示してもよい。
【０１５３】
さらに、モード切替手段２６によって、画像の表示モードを通常表示モードとその他の表示モードとで切り替えるようにしても良い。
【０１５４】
以上説明した例によれば、スペックルの位相凹凸変化の分布の偏りから、仮定音速が最適音速より速いか遅いか、またはどれだけずれているかが判定できるため、これから最適音速を判定することができる。
【０１５５】
以上、本発明の超音波診断方法及び装置について詳細に説明したが、本発明は、以上の例には限定されず、本発明の要旨を逸脱しない範囲において、各種の改良や変形を行ってもよいのはもちろんである。
【図面の簡単な説明】
【０１５６】
【図１】本発明に係る超音波診断装置の一実施形態の概略構成を示すシステム構成図である。
【図２】仮定音速１４００［ｍ／ｓ］〜１５００［ｍ／ｓ］における微小構造物信号の仮定音速に依る位相変化特性を示すグラフである。
【図３】仮定音速１５００［ｍ／ｓ］〜１６２０［ｍ／ｓ］における微小構造物信号の仮定音速に依る位相変化特性を示すグラフである。
【図４】仮定音速１４００［ｍ／ｓ］〜１４８０［ｍ／ｓ］における面信号の仮定音速に依る位相変化特性を示すグラフである。
【図５】仮定音速１５４０［ｍ／ｓ］〜１６２０［ｍ／ｓ］における面信号の仮定音速に依る位相変化特性を示すグラフである。
【図６】仮定音速１４００［ｍ／ｓ］〜１４８０［ｍ／ｓ］におけるスペックルの仮定音速に依る位相変化特性を示すグラフである。
【図７】仮定音速１５４０［ｍ／ｓ］〜１６２０［ｍ／ｓ］におけるスペックルの仮定音速に依る位相変化特性を示すグラフである。
【図８】仮定音速１４００［ｍ／ｓ］〜１５００［ｍ／ｓ］における微小構造物の仮定音速に依る振幅変化特性を示すグラフである。
【図９】仮定音速１５００［ｍ／ｓ］〜１６２０［ｍ／ｓ］における微小構造物信号の仮定音速に依る振幅変化特性を示すグラフである。
【図１０】超音波信号受信状態を示す説明図である。
【図１１】式（１）と式（２）で与えられる曲線が（Ｘ，ｔ）平面において接している様子を示す説明図であり、（ａ）はＶ＜Ｖ０の場合であり、（ｂ）は、Ｖ＝Ｖ０の場合であり、（ｃ）はＶ＞Ｖ０の場合である。
【図１２】画像生成部の作用を示すフローチャートである。
【図１３】最適音速判定部の作用を示すフローチャートである。
【図１４】音速判定画像値と仮定音速との関係を示す線図である。
【図１５】振幅を利用した音速判定画像値と仮定音速との関係を示す線図である。
【図１６】表示画像生成部における処理内容を示すフローチャートである。
【図１７】スペックルの仮定音速による位相変化を示す説明図であり、（ａ）は、最適音速より速い場合、（ｂ）は最適音速より遅い場合を示す。
【図１８】微小構造物における位相凹凸変化のヒストグラムを示すグラフであり、（ａ）は仮定音速が最適音速より遅い場合、（ｂ）は仮定音速が最適音速より速い場合を示す。
【図１９】ある対象領域における仮定音速の位相凹凸変化から最適音速を判定する手順を示すフローチャートである。
【図２０】仮定音速による分布の偏りの変化を示すグラフである。
【図２１】表示画像生成部における処理の流れを示すフローチャートである。
【符号の説明】
【０１５７】
１…超音波診断装置、１０…超音波プローブ、１２…送受信部、１４…走査制御部、１６…ＡＤ変換部、１８…画像生成部、２０…最適音速判定部、２２…表示画像生成部、２４…モニタ、２６…モード切替手段
【図１】
2009101145.tif 000003
【図２】
2009101145.tif 000004
【図３】
2009101145.tif 000005
【図４】
2009101145.tif 000006
【図５】
2009101145.tif 000007
【図６】
2009101145.tif 000008
【図７】
2009101145.tif 000009
【図８】
2009101145.tif 000010
【図９】
2009101145.tif 000011
【図１０】
2009101145.tif 000012
【図１１】
2009101145.tif 000013
【図１２】
2009101145.tif 000014
【図１３】
2009101145.tif 000015
【図１４】
2009101145.tif 000016
【図１５】
2009101145.tif 000017
【図１６】
2009101145.tif 000018
【図１７】
2009101145.tif 000019
【図１８】
2009101145.tif 000020
【図１９】
2009101145.tif 000021
【図２０】
2009101145.tif 000022
【図２１】
2009101145.tif 000023