"超音波を被検体に送信して反射波を受信し超音波検出信号を出力する超音波探触子を有し、形状を表す形態画像を生成する超音波画像生成装置であって、前記形態画像の着目領域における音速値を算出する音速値算出部と、前記着目領域に対応する前記音速値に基づいて、該音速値の信頼度情報を生成する信頼度情報生成部と、を有することを特徴とする超音波画像生成装置。さらに、前記着目領域に対応する前記音速値に基づいて、音速画像を生成する音速画像生成部を有し、前記信頼度情報生成部は、前記音速画像の前記着目領域に対応する領域ごとに前記信頼度情報を生成することを特徴とする請求項1に記載の超音波画像生成装置。さらに、前記信頼度情報に基づき信頼度画像を生成する信頼度画像生成部を有することを特徴とする請求項2に記載の超音波画像生成装置。さらに、前記形態画像、前記音速画像、および前記信頼度画像のうち1以上から表示画像を生成する画像処理部と、前記表示画像を表示する表示部と、を有することを特徴とする請求項3に記載の超音波画像生成装置。前記信頼度情報は、前記音速画像の各画素および該各画素の周囲における所定範囲の画素から求めた情報であることを特徴とする請求項2~4のいずれかに記載の超音波画像生成装置。前記信頼度情報は、複数の前記音速画像間の各画素から求めた情報であることを特徴とする請求項2~4のいずれかに記載の超音波画像生成装置。前記信頼度情報は、前記音速画像に対してエッジ検出処理またはギャップ検出処理により不連続点を検出した場合の、測定エラー情報であることを特徴とする請求項2~4のいずれかに記載の超音波画像生成装置。前記音速値は環境音速値であり、前記音速画像は環境音速画像であることを特徴とする請求項2~7のいずれかに記載の超音波画像生成装置。前記音速値算出部は、環境音速値を算出する環境音速値算出部と、前記環境音速値に基づき局所音速値を算出する局所音速値算出部とを有し、前記音速値は局所音速値であり、前記音速画像は局所音速画像であることを特徴とする請求項2~7のいずれかに記載の超音波画像生成装置。前記音速値算出部は、前記形態画像の各画素および該各画素の周囲における所定範囲の画素から、該各画素のフォーカス指標値を設定音速ごとに求め、または該各画素における仮想受信波と仮想合成受信波との誤差を、仮の局所音速ごとに求め、前記信頼度情報生成部は、前記設定音速ごとの前記フォーカス指標値、または前記仮の局所音速ごとの前記誤差に基づいて信頼度情報を生成することを特徴とする請求項1に記載の超音波画像生成装置。さらに、音速画像生成部を有し、前記音速値算出部は、さらに、前記設定音速ごとの前記フォーカス指標値、または前記仮の局所音速ごとの前記誤差に基づき、前記形態画像の着目領域における環境音速値または局所音速値を算出し、前記音速画像生成部は、前記着目領域に対応する前記環境音速値または前記局所音速値に基づいて音速画像を生成し、前記信頼度情報生成部は、さらに、前記音速画像の前記着目領域に対応する領域ごとに信頼度情報を生成することを特徴とする請求項10に記載の超音波画像生成装置。さらに、前記信頼度情報に基づき信頼度画像を生成する信頼度画像生成部を有することを特徴とする請求項10または11に記載の超音波画像生成装置。さらに、前記形態画像、前記音速画像、および前記信頼度画像のうち1以上から表示画像を生成する画像処理部と、前記表示画像を表示する表示部と、を有することを特徴とする請求項12に記載の超音波画像生成装置。前記表示部に、前記形態画像と、前記音速画像または前記信頼度画像とをオーバーレイ表示することを特徴とする請求項4~9,13のいずれかに記載の超音波画像生成装置。前記信頼度画像の各画素のうち、前記信頼度情報が所定の値よりも低い画素をマスクして表示することを特徴とする請求項4または13に記載の超音波画像生成装置。超音波を被検体に送信して反射波を受信し、形状を表す形態画像を生成する超音波画像生成方法であって、前記形態画像の着目領域における音速値を算出する音速値算出ステップと、前記着目領域に対応する前記音速値に基づいて、該音速値の信頼度情報を生成する信頼度情報生成ステップと、を有することを特徴とする超音波画像生成方法。さらに、前記着目領域に対応する前記音速値に基づいて、音速画像を生成する音速画像生成ステップを有し、前記信頼度情報生成ステップは、前記音速画像の前記着目領域に対応する領域ごとに前記信頼度情報を生成することを特徴とする請求項16に記載の超音波画像生成方法。前記音速値算出ステップは、前記形態画像の各画素および該各画素の周囲における所定範囲の画素から、該各画素のフォーカス指標値を設定音速ごとに求め、または該各画素における仮想受信波と仮想合成受信波との誤差を、仮の局所音速ごとに求め、前記信頼度情報生成ステップは、前記設定音速ごとの前記フォーカス指標値、または前記仮の局所音速ごとの前記誤差に基づいて信頼度情報を生成することを特徴とする請求項16に記載の超音波画像生成方法。さらに、音速画像生成ステップを有し、前記音速値算出ステップは、さらに、前記設定音速ごとの前記フォーカス指標値、または前記仮の局所音速ごとの前記誤差に基づき、前記形態画像の着目領域における環境音速値または局所音速値を算出し、前記音速画像生成ステップは、前記着目領域に対応する前記環境音速値または前記局所音速値に基づいて音速画像を生成し、前記信頼度情報生成ステップは、さらに、前記音速画像の前記着目領域に対応する領域ごとに信頼度情報を生成することを特徴とする請求項18に記載の超音波画像生成方法。さらに、前記信頼度情報に基づき信頼度画像を生成する信頼度画像生成ステップを有することを特徴とする請求項17~19のいずれかに記載の超音波画像生成方法。請求項16~20のいずれかに記載の超音波画像生成方法の各ステップをコンピュータに実行させるためのプログラム。請求項21に記載のプログラムが記録されたコンピュータ読み取り可能な記録媒体。しかしながら、特許文献1では、Bモード画像のフォーカスは向上するが、画像中の任意の領域における音速値を求めることはできない。また、特許文献2では、着目領域の局所音速値を求めること、求めた音速値によってフォーカスを向上させ、Bモード画像を再構築して高画質化すること、およびBモード画像に局所音速値の判定結果を重畳して表示する例えば、局所音速値に応じて色分けまたは輝度を変化させることは開示されているが、音速値の信頼性を求めることは開示されていない。本発明の目的は、着目領域の環境音速値または局所音速値だけでなく、これらの信頼度を求めることで、信頼性の高い環境音速値または局所音速値のみをフォーカスに用いたり、これらの値を直接または画像化して表示することで、より被検体内の病変の検出が行いやすい超音波画像を生成する超音波画像生成装置、超音波画像生成方法およびプログラムを提供することにある。上記課題を解決するために、本発明は、超音波を被検体に送信して反射波を受信し超音波検出信号を出力する超音波探触子を有し、形状を表す形態画像を生成する超音波画像生成装置であって、前記形態画像の着目領域における音速値を算出する音速値算出部と、前記着目領域に対応する前記音速値に基づいて、該音速値の信頼度情報を生成する信頼度情報生成部と、を有することを特徴とする超音波画像生成装置を提供する。さらに、前記着目領域に対応する前記音速値に基づいて、音速画像を生成する音速画像生成部を有し、前記信頼度情報生成部は、前記音速画像の前記着目領域に対応する領域ごとに前記信頼度情報を生成することが好ましい。さらに、前記信頼度情報に基づき信頼度画像を生成する信頼度画像生成部を有することが好ましい。さらに、前記形態画像、前記音速画像、および前記信頼度画像のうち1以上から表示画像を生成する画像処理部と、前記表示画像を表示する表示部と、を有することが好ましい。また、前記信頼度情報は、前記音速画像の各画素および該各画素の周囲における所定範囲の画素から求めた情報であることが好ましい。また、前記信頼度情報は、複数の前記音速画像間の各画素から求めた情報であることが好ましい。また、前記信頼度情報は、前記音速画像に対してエッジ検出処理またはギャップ検出処理により不連続点を検出した場合の、測定エラー情報であることが好ましい。また、前記音速値は環境音速値であり、前記音速画像は環境音速画像であることが好ましい。また、前記音速値算出部は、環境音速値を算出する環境音速値算出部と、前記環境音速値に基づき局所音速値を算出する局所音速値算出部とを有し、前記音速値は局所音速値であり、前記音速画像は局所音速画像であることが好ましい。また、前記音速値算出部は、前記形態画像の各画素および該各画素の周囲における所定範囲の画素から、該各画素のフォーカス指標値を設定音速ごとに求め、または該各画素における仮想受信波と仮想合成受信波との誤差を、仮の局所音速ごとに求め、前記信頼度情報生成部は、前記設定音速ごとの前記フォーカス指標値、または前記仮の局所音速ごとの前記誤差に基づいて信頼度情報を生成することが好ましい。さらに、音速画像生成部を有し、前記音速値算出部は、さらに、前記設定音速ごとの前記フォーカス指標値、または前記仮の局所音速ごとの前記誤差に基づき、前記形態画像の着目領域における環境音速値または局所音速値を算出し、前記音速画像生成部は、前記着目領域に対応する前記環境音速値または前記局所音速値に基づいて音速画像を生成し、前記信頼度情報生成部は、さらに、前記音速画像の前記着目領域に対応する領域ごとに信頼度情報を生成することが好ましい。さらに、前記信頼度情報に基づき信頼度画像を生成する信頼度画像生成部を有することが好ましい。さらに、前記形態画像、前記音速画像、および前記信頼度画像のうち1以上から表示画像を生成する画像処理部と、前記表示画像を表示する表示部と、を有することが好ましい。また、前記表示部に、前記形態画像と、前記音速画像または前記信頼度画像とをオーバーレイ表示することが好ましい。また、前記信頼度画像の各画素のうち、前記信頼度情報が所定の値よりも低い画素をマスクして表示することが好ましい。また、上記課題を解決するために、本発明は、超音波を被検体に送信して反射波を受信し、形状を表す形態画像を生成する超音波画像生成方法であって、前記形態画像の着目領域における音速値を算出する音速値算出ステップと、前記着目領域に対応する前記音速値に基づいて、該音速値の信頼度情報を生成する信頼度情報生成ステップと、を有することを特徴とする超音波画像生成方法を提供する。さらに、前記着目領域に対応する前記音速値に基づいて、音速画像を生成する音速画像生成ステップを有し、前記信頼度情報生成ステップは、前記音速画像の前記着目領域に対応する領域ごとに前記信頼度情報を生成することが好ましい。また、前記音速値算出ステップは、前記形態画像の各画素および該各画素の周囲における所定範囲の画素から、該各画素のフォーカス指標値を設定音速ごとに求め、または該各画素における仮想受信波と仮想合成受信波との誤差を、仮の局所音速ごとに求め、前記信頼度情報生成ステップは、前記設定音速ごとの前記フォーカス指標値、または前記仮の局所音速ごとの前記誤差に基づいて信頼度情報を生成することが好ましい。さらに、音速画像生成ステップを有し、前記音速値算出ステップは、さらに、前記設定音速ごとの前記フォーカス指標値、または前記仮の局所音速ごとの前記誤差に基づき、前記形態画像の着目領域における環境音速値または局所音速値を算出し、前記音速画像生成ステップは、前記着目領域に対応する前記環境音速値または前記局所音速値に基づいて音速画像を生成し、前記信頼度情報生成ステップは、さらに、前記音速画像の前記着目領域に対応する領域ごとに信頼度情報を生成することが好ましい。さらに、前記信頼度情報に基づき信頼度画像を生成する信頼度画像生成ステップを有することが好ましい。また、上記課題を解決するために、本発明は、上記のいずれかに記載の超音波画像生成方法の各ステップをコンピュータに実行させるためのプログラムを提供する。また、上記課題を解決するために、本発明は、上記に記載のプログラムが記録されたコンピュータ読み取り可能な記録媒体を提供する。本発明によれば、環境音速値または局所音速値を測定するための、専用の超音波を送受信する構成を用いることなく、より被検体内の病変の検出が行いやすい、環境音速画像、または局所音速画像、および信頼度画像を提供することができる。本発明に係る超音波画像生成方法を実施する本発明に係る超音波画像生成装置を、添付の図面に示す好適実施形態に基づいて以下に詳細に説明する。なお、以下の実施形態では、形態画像としてBモード画像を用いる場合について説明する。図1は、本発明に係る超音波画像生成装置の構成を表す第1実施形態のブロック図である。第1実施形態の超音波画像生成装置では、音速値として環境音速値を用いて音速画像および信頼度画像が生成される。図1に示す超音波画像生成装置10は、操作部12、制御部14、超音波探触子16、送受信部20、信号処理部22、環境音速値算出部24、音速画像生成部32、第1信頼度画像生成部34、第2信頼度画像生成部36、画像処理部38、表示部40、およびRFデータ記録再生部42によって構成される。また、環境音速値算出部24は、フォーカス指標算出部26、環境音速プロファイル生成部28、および環境音速値決定部30によって構成される。操作部12は、超音波画像生成装置10の各種操作をオペレータが行うためのものであり、操作情報を出力する。操作部12の具体的な態様には特に限定はなく、キーボード、マウス、タッチパネルなど、公知の各種の操作機器を用いればよい。制御部14は、超音波画像生成装置10の各部の動作を制御するためのものである。また、操作情報に従って各種処理が実施されるように、各部に対して制御信号CTLが出力される。また、環境音速値を求めるための設定音速または受信遅延パターンを、後述する送受信部20に設定する。超音波探触子16は、被検体に当接させて用いるプローブであり、1次元または2次元のトランスデューサアレイを構成する複数の超音波トランスデューサ18を備えている。超音波トランスデューサ18では、送受信部20から印加される駆動信号に基づいて超音波ビームが被検体に送信されると共に、被検体から反射される超音波エコーが受信されて検出信号が出力される。超音波トランスデューサ18は、圧電性を有する材料圧電体の両端に電極が形成されて構成された振動子を含んでいる。上記振動子を構成する圧電体としては、例えば、PZTチタン酸ジルコン酸鉛:Pbzirconatetitanateのような圧電セラミック、PVDFポリフッ化ビニリデン:polyvinylidenedifluorideのような高分子圧電素子を用いることができる。上記振動子の電極に電気信号を送って電圧を印加すると圧電体が伸縮し、この圧電体の伸縮により各振動子において超音波が発生する。例えば、振動子の電極にパルス状の電気信号を送るとパルス状の超音波が発生し、振動子の電極に連続波の電気信号を送ると連続波の超音波が発生する。そして、各振動子において発生した超音波が合成されて超音波ビームが形成される。また、各振動子により超音波が受信されると、各振動子の圧電体が伸縮して電気信号を発生する。各振動子において発生した電気信号は、超音波の検出信号として送受信部20に出力される。なお、超音波トランスデューサ18としては、超音波変換方式の異なる複数種類の素子を用いることも可能である。例えば、超音波を送信する素子として上記圧電体により構成される振動子を用いて、超音波を受信する素子として光検出方式の超音波トランスデューサを用いるようにしてもよい。ここで、光検出方式の超音波トランスデューサとは、超音波信号を光信号に変換して検出するものであり、例えば、ファブリーペロー共振器又はファイバブラッググレーティングである。送受信部20は、送信回路、受信回路、およびA/D変換器を備えている。送信回路は、制御部14からの制御信号に応じて駆動信号を生成して、該駆動信号を超音波トランスデューサ18に印加する。このとき、送信回路は、制御部14によって選択された送信遅延パターンに基づいて、各超音波トランスデューサ18に印加する駆動信号を遅延させる。ここで、送信回路は、複数の超音波トランスデューサ18から送信される超音波が超音波ビームを形成するように、各超音波トランスデューサ18に駆動信号を印加するタイミングを調整する遅延させる。なお、複数の超音波トランスデューサ18から一度に送信される超音波が被検体の撮像領域全体に届くように、駆動信号を印加するタイミングを調節するようにしてもよい。受信回路は、各超音波トランスデューサ18から出力される超音波検出信号を受信して増幅する。各超音波トランスデューサ18と被検体内の超音波反射源との間の距離がそれぞれ異なるため、各超音波トランスデューサ18に反射波が到達する時間は異なる。受信回路は遅延回路を備えており、制御部14によって選択された音速以下、仮定音速という。または音速の分布に基づいて設定される受信遅延パターンに従って、反射波の到達時刻の差遅延時間に相当する分、各検出信号を遅延させる。次に、受信回路は、遅延時間を与えた検出信号を整合加算することにより受信フォーカス処理を行う。超音波反射源Xと異なる位置に別の超音波反射源がある場合には、別の超音波反射源からの超音波検出信号は到達時刻が異なるので、加算回路で加算することにより、別の超音波反射源からの超音波検出信号の位相が打ち消し合う。これにより、超音波反射源Xからの受信信号が最も大きくなり、フォーカスが合う。この受信フォーカス処理によって、超音波エコーの焦点が絞り込まれた音線信号以下、RF信号という。が形成される。A/D変換器では、受信回路から出力されるアナログのRF信号がデジタルRF信号以下、RFデータという。に変換され出力される。ここで、RFデータには、受信波搬送波の位相情報が含まれている。信号処理部22では、RFデータに対して、STCSensitivityTimegainControlによって、超音波の反射位置の深度に応じて距離による減衰の補正がされた後、包絡線検波処理が施され、Bモード画像データが生成され出力される。フォーカス指標算出部26には、設定音速ごとのBモード画像データおよびRFデータが入力され、Bモード画像の画素ごとに、環境音速値を求めるために必要な設定音速ごとのフォーカス指標が算出される。ある設定音速についてフォーカス指標が算出されると、設定音速を変更してフォーカス指標が算出される。つまり、全ての設定音速についてフォーカス指標が算出され出力される。フォーカス指標としては、例えば、Bモード画像データからは画像のコントラスト、シャープネス、あるいは各画素における超音波検出信号の周期または振幅が用いられ、RFデータからはビーム幅が用いられる。環境音速プロファイル生成部28には、Bモード画像の画素ごとに全ての設定音速についてのフォーカス指標が入力される。入力されたフォーカス指標は、横軸を設定音速、縦軸をフォーカス指標としたグラフにプロットされ、設定音速プロファイル以下、環境音速プロファイルというが生成され出力される。環境音速値決定部30には、Bモード画像の画素ごとに求められた環境音速プロファイルが入力される。入力された環境音速プロファイルに基づき、画素ごとの最適音速値以下、環境音速値というが決定され出力される。ここで、最適音速値環境音速値とは、画像のコントラスト、シャープネスが最も高くなる音速値、およびビーム幅が最も狭くなる音速値であり、各画素における実際の音速値局所音速値とは必ずしも一致しない。なお、環境音速値を求める方法は、他にも、例えばスキャン方向の空間周波数、分散などから判定する方法例えば、特開平8317926号公報を用いてもよい。音速画像生成部32には、Bモード画像の画素ごとに求められた環境音速値が入力される。音速画像生成部32では、環境音速値に対応した値、例えば、階調表現が可能なように、一定範囲の音速について対応付けられた一定範囲の数値が、Bモード画像の各画素に割り当てられた画像以下、音速画像というが生成され、音速画像データとして出力される。なお、階調表現が可能な一定範囲の数値に対して、擬似カラー画像化が行われてもよい。さらに、超音波探触子16が静止された状態で複数の音速画像が生成され、時間軸フレーム方向の当該複数の音速画像間の平均値を算出することで測定エラーを抑制した測定エラー抑制音速画像としてもよい。第1信頼度画像生成部34には、音速画像データが入力される。第1信頼度画像生成部34は信頼度情報生成部および信頼度画像生成部であり、音速画像の画素ごとに信頼度情報が生成され、信頼度情報が画像化された信頼度画像が生成され、信頼度画像データとして出力される。ここで、信頼度情報は、例えば、音速画像の各画素および該各画素の周囲における所定範囲例えば、3×3画素の画素における環境音速値のバラつきの標準偏差を求めて、中心画素の値とすることができる。全ての画素について求められた信頼度情報は画像化され、信頼度画像標準偏差画像データとして出力される。また、信頼度情報は、他にも取得時間の異なる同一着目領域の、複数の音速画像間における対応する画素の環境音速値のバラつきについて標準偏差を求めてもよく、同様に画像化され信頼度画像データとして出力される。さらに、これらを組み合わせて、つまり、画像のx軸、y軸、および時間軸フレームの3次元の所定範囲の画素における環境音速値のバラつきについて標準偏差を求めてもよい。また、各画素における環境音速が一定だとしても、超音波探触子16と着目領域の間に音速の異なる媒質が存在した場合、環境音速は深さ方向に変化する。そこで、この変化を信頼性の計算から除外するために、1次式または2次式で深さ方向に変化する平面で近似し、各画素における環境音速値から本近似による環境音速値を引いた後に標準偏差を求めるようにしてもよい。さらに、超音波探触子16が移動されて複数の音速画像が生成される場合に、同一着目画素の各音速画像各フレーム間における環境音速値のバラつきの標準偏差を求めるようにしてもよい。第2信頼度画像生成部36には、Bモード画像の画素ごとに求められた環境音速プロファイルが入力される。第2信頼度画像生成部36は信頼度情報生成部および信頼度画像生成部であり、Bモード画像の画素ごとに、入力された環境音速プロファイルの形状または音速画像に基づき、信頼度情報が生成され、信頼度情報が画像化された信頼度画像が生成され、信頼度画像データが出力される。なお、信頼度情報生成部および信頼度画像生成部は、第1信頼度画像生成部34と第2信頼度画像生成部36とが、単独で、または併せて構成されたものである。ここで、環境音速プロファイルの形状について説明すると、環境音速値の信頼性が高い場合、例えば、図2aに示すように、フォーカス指標は環境音速値に近い程単調に大きくなる単峰性の変化となる。また、環境音速値におけるフォーカス指標と、環境音速値から外れた設定音速値におけるフォーカス指標とでは、その差が大きくなる。つまり、コントラストが大きくなる。一方、環境音速値の信頼性が低い場合、例えば、図2bに示すように、フォーカス指標は設定音速によって複雑に変化し、2峰性や3峰性の変化となる場合もある。また、環境音速値におけるフォーカス指標と、環境音速値から外れた設定音速値におけるフォーカス指標との差が小さくなる。つまり、コントラストが小さくなる。このため、環境音速プロファイルの形状の歪さ、例えば、曲線フィッティングの誤差が小さくなる次数、最大値を頂点とした2次曲線フィッティング結果との誤差、設定音速値の分散2次モーメント、および設定音速値の3次モーメント等を用いて、信頼度情報を生成することができる。また、フォーカス指標のコントラスト、例えば、フォーカス指標の最大値と最小値との差分、下記式1で表されるフォーカス指標、およびフォーカス指標の微分絶対値の平均値を用いて信頼度情報を生成してもよい。フォーカス指標=最大フォーカス指標-フォーカス指標1400フォーカス指標1650/21また、フォーカス指標がBモード画像の輝度や周波数に依存する事を踏まえて、上記の誤差またはコントラストを、フォーカス指標の平均値や最大値で規格化した値を用いて信頼度情報を生成してもよい。画像処理部38には、Bモード画像データ、音速画像データ、および信頼度画像データが入力される。画像処理部38は、DSCDigitalScanConverter機能、並びに、エッジ検出、ギャップ検出、各種画像データBモード画像データ、音速画像データ、および信頼度画像データの重畳オーバーレイ表示、強調表示、およびマスク処理等の画像処理機能を有する。画像処理部38からは、DSCおよび画像処理が行われた表示画像データが出力される。なお、重畳表示等に用いるBモード画像データは、画像全体のフォーカスが最も良好な設定音速におけるデータを用いるのがよい。DSC機能では、Bモード画像データ、およびこれを基に生成される音速画像データおよび信頼度画像データは、通常のテレビジョン信号の走査方式と異なる走査方式であるため、後述する表示部40に表示可能なように、通常の画像データ例えば、テレビジョン信号の走査方式NTSC方式の画像データへの変換ラスター変換が行われる。画像処理機能では、例えば、各種画像データの重畳画像を生成したり、信頼度の高い画素を強調表示したり、信頼度の低い画素をマスクして非表示とすることで、医師の診断を補助する画像が生成される。例えば、信頼性画像の各画素のうち、信頼性情報が所定値以下の着目画素を測定エラーとして、測定エラー情報を取得して測定エラー画像を生成し、測定エラー画像を音速画像に重畳表示するようにしてもよい。また、信頼性画像として、例えば、音速画像についてエッジ検出処理またはギャップ検出処理を行う、すなわち、着目画素周辺の2次元断面内、または時間軸フレーム方向も含めて3次元的に複数の画素を設定し、隣り合う各画素間の環境音速値の差が所定以上の画素組を含む場合つまり、不連続箇所を含む場合、当該着目画素を測定エラーとして検出することで、測定エラー情報を取得して測定エラー画像を生成し、測定エラー画像を音速画像に重畳表示するようにしてもよい。このとき、深さ方向に1次式または2次式で変化する平面で近似し、各画素における環境音速値から、本近似による環境音速値を引いた後に不連続箇所を判定し、測定エラーとして検出してもよい。表示部40には、表示画像データが入力され表示される。表示部40は、液晶、プラズマ、有機ELElectroLuminescence等のFPDFlatPanelDisplay、またはCRTCathodeRayTube等により構成される。表示部40は、複数の画像を並べて表示可能なように、表示面積が大きいもの、および画素数の多いものを用いるのがよい。RFデータ記録再生部42には、RFデータ、フレームレートに関する情報例えば、超音波の反射位置の深度、走査線の密度、視野幅を示すパラメータが入力され、内部のシネメモリに記録される。RFデータ記録再生部42は、シネメモリ記録モードとシネメモリ再生モードの2つの動作モードを持ち、通常観察時ライブモードにはシネメモリ記録モードとして動作し、RFデータが記録されている。シネメモリ再生モードは、シネメモリに格納されているRFデータに基づいて超音波診断画像の表示、解析計測を行うモードである。シネメモリ再生モード時には、シネメモリに格納されているRFデータが、オペレータの操作に応じて、信号処理部22へと出力され、オペレータはRFデータ記録再生部42に記録されたRFデータに基づく、Bモード画像、音速画像、および信頼度画像を見ることができる。次に、本発明に係る超音波画像生成方法を実現する、本発明に係る超音波画像生成装置10の動作を説明する。図3および図4は、本発明に係る超音波画像生成方法の処理の流れの一例を示すフローチャートである。まず、環境音速値を求めるための設定音速が制御部14によって複数選択され、送受信部20に対して設定されるステップS8。ここで、設定音速は人体内の音速の範囲概ね、1400m/s~1650m/sの中から、複数の音速が選択され設定される。選択される音速の数Lは、横軸を設定音速、縦軸をフォーカス指標としたグラフにプロットしたときに、図2a,bに示すようなグラフが描ける程度であればよい。なお、オペレータによって設定音速が設定されるようにしてもよい。続いて、オペレータが超音波探触子16を被検体に接触させて超音波の送受信が行われ、超音波検出信号が超音波探触子16から出力されるステップS10。超音波検出信号は送受信部20に入力され、選択された複数の設定音速に対応する受信遅延パターンに基づき、それぞれの設定音速ごとに受信フォーカス処理または送受信フォーカス処理が行われ、A/DAnalog/Digital変換されて、設定音速ごとのRFデータとして出力されるステップS12。設定音速ごとのRFデータは、信号処理部22に入力され、STCにより超音波の反射位置の深度に応じて距離による減衰の補正がされた後、包絡線検波処理が施され、図5に示すようなBモード画像のBモード画像データが設定音速ごとに生成され出力されるステップS14。全ての設定音速ごとのBモード画像データおよびRFデータは、環境音速値算出部24に入力される。環境音速値算出部24では、Bモード画像のうち音速画像または信頼性画像を得たい領域について着目領域として指定されるステップS16。着目領域は、例えば、図6に示す着目領域44のように、Bモード画像全体が初期設定にて自動的に設定されるようにしてもよいし、Bモード画像の一部が自動的に着目領域として設定されてもよい。また、オペレータが操作部12を操作することによりBモード画像の一部が設定されてもよい。着目領域は、例えば、始点座標[xmin,ymin]、および終点座標[xmax,ymax]で指定される。図3のフローチャートでは、例としてxがn~N,yがm~Mと設定される。着目領域が設定されると、環境音速値の算出を開始する開始着目画素例えば、x=n,y=mが設定されステップS18、着目画素の環境音速値の算出が行われるステップS20。ここで、着目画素の環境音速値の算出の詳細を、図4のフローチャートにより説明する。まず、ステップS8で設定された環境音速値を求めるための選択される音速の数設定音速の個数が、初期値Ci=1,最大値Cmax=Lに設定されるステップS202。次に、設定音速Ciのフォーカス指標が算出され出力されるステップS204。フォーカス指標値としては、例えば、Bモード画像データのコントラスト、シャープネスの値が算出され出力される。なお、着目画素のRFデータのビーム幅から、所定の指数が算出されてフォーカス指標として出力されてもよい。Ci=1初期値についてフォーカス指標の算出が終了すると、CiとCmax最大値の値が比較されステップS206、Ciの値がCmax未満であるときはステップS206でN、Ciに1が足されてステップS208、ステップS204のフォーカス指標の算出に戻る。フォーカス指標の算出ステップS204は、Ci=Cmaxとなるまで繰り返され、着目画素の全ての設定音速についてフォーカス指標が算出され出力される。着目画素の全ての設定音速についてのフォーカス指標は、環境音速プロファイル生成部28に入力され、横軸を設定音速、縦軸をフォーカス指標としたグラフにプロットされ、環境音速プロファイルが生成され出力されるステップS210。環境音速プロファイルは、環境音速値決定部30に入力され、例えば、図2aに示すような環境音速プロファイルであれば、フォーカス指標の最大値の設定音速値が環境音速値として決定され出力されるステップS212。着目画素の環境音速値の算出が終了すると、つまりステップS20が終了すると、着目画素のy座標の値がymaxと比較されステップS22、yの値がymax未満であるときはステップS22でN、yに1が足されてステップS24、ステップS20の着目画素の環境音速値の算出に戻る。着目画素の環境音速値の算出ステップS20は、y=ymaxとなるまで繰り返される。y=ymaxとなるとステップS22でY、着目画素のx座標の値がxmaxと比較されステップS26、xの値がxmax未満であるときはステップS26でN、xに1が足されてステップS28、y座標の値がyminy=mに設定されステップS30、ステップS20の着目画素の環境音速値の算出に戻る。つまり、y座標方向をラインとすると、x座標がnである1ライン目の環境音速値が算出されると、x座標が1だけインクリメントされn1、2ライン目の環境音速値が算出される。着目画素の環境音速値の算出ステップS20は、着目領域全体についてx=xmaxx=N,y=ymaxy=Mまで、環境音速値が算出されるまで繰り返される。着目領域全体について環境音速値の算出が終了すると、環境音速値は音速画像生成部32に入力され、環境音速プロファイルは、第2信頼度画像生成部36に入力される。ここで、信頼度画像の選択情報が予め設定されるか、またはオペレータにより操作部12が操作され制御部14により選択されるステップS32。信頼度画像として標準偏差画像が選択された場合ステップS32で標準偏差、音速画像生成部32では、Bモード画像の各画素に対応した環境音速値が割り当てられ、環境音速画像が生成され環境音速画像データ出力されるステップS34。環境音速画像データは、第1信頼度画像生成部34、および画像処理部38に入力される。第1信頼度画像生成部34では、音速画像の画素ごとに信頼度情報が生成され、信頼度情報が画像化された信頼度画像が生成され、信頼度画像データとして出力されるステップS36。一方、ステップS32において、信頼度画像として環境音速プロファイル形状画像が選択された場合ステップS32でプロファイル形状、環境音速画像を用いずに、第2信頼度画像生成部36に入力された、Bモード画像の画素ごとの環境音速プロファイル形状に基づきBモード画像の画素ごとに信頼度情報が生成され、信頼度情報が画像化された信頼度画像が生成され、信頼度画像データが出力されるステップS36。画像処理部38には、Bモード画像データ、環境音速画像データ、および信頼度画像データが入力される。画像処理部38では、入力された各種画像データに対して通常の画像データへの変換、重畳表示、強調表示、およびマスク処理等が行われ、表示画像データが生成され、表示部40に対して出力され表示されるステップS38。なお、各種画像データに対してエッジ検出、ギャップ検出等が行われ、測定エラー情報を取得して測定エラー画像が生成されてもよい。ここで、表示部40に表示される表示画像データの例を挙げる。例えば、図5に示すようなBモード画像または図7に示すような環境音速画像と、図8に示すような信頼度画像標準偏差画像または図9に示すような信頼度画像音速プロファイル形状とを並べて表示することができる。また、例えば、図10に示すように、環境音速画像の色輝度、色相、彩度を変調して、単独または並べて表示してもよいし、環境音速画像の代わりに信頼性画像または測定エラー画像の色を変調して、単独または並べて表示してもよい。さらに、図11に示すように、図5のBモード画像と図7の環境音速画像とを重畳表示してもよい。また、Bモード画像、または環境音速画像の表示を、信頼性が低い画素または測定エラーが発生している画素の表示を行わないよう、あるいは、信頼性が高い画素のみ表示するように、限定して表示してもよい。例えば、図12に示すように、図7の環境音速画像を図8の信頼性画像でマスクして表示することができる。他にも、例えば、図13に示すように、図5に示すBモード画像と、図10に示す環境音速画像の色を変調した画像を重畳表示してもよいし、図14に示すように、図5に示すBモード画像と、図12に示す環境音速画像を信頼性画像でマスクした画像を重畳表示してもよい。また、上記の各種の表示画像をオペレータが操作部12を介して、任意に表示モードが切り替えられるようにしてもよい。以上のように、環境音速値を測定するための専用の超音波を送受信する構成を用いることなく、より被検体内の病変の検出が行いやすい、様々な環境音速画像、および信頼度画像を提供することができる。次に、第2実施形態として、環境音速値の代わりに局所音速値を用いた場合について説明する。図15は、上述の第1実施形態と異なる、本発明に係る超音波画像生成装置50の構成を表す第2実施形態のブロック図である。超音波画像生成装置50は、図1に示す本発明に係る第1実施形態の超音波画像生成装置10と比べ、環境音速値から局所音速値をさらに求め、当該局所音速値を用いて音速画像および信頼度画像を生成するものであり、基本的に同様の構成を有するものである。よって、同様の構成要素には同一の符号を付し、その詳細な説明は省略する。図15に示す超音波画像生成装置50は、環境音速値算出部24の出力である各画素の環境音速値が局所音速値算出部52に入力され、局所音速値算出部52から出力される局所音速値が音速画像生成部32および第2信頼度画像生成部36に入力され、同じく局所音速値算出部52から出力される誤差プロファイルが第2信頼度画像生成部36に入力される。局所音速値算出部52は、仮想受信波仮想合成受信波誤差算出部54、誤差プロファイル生成部56、および局所音速値決定部58で構成される。仮想受信波仮想合成受信波誤差算出部54には、Bモード画像データおよび環境音速値が入力され、Bモード画像の画素ごとに、局所音速値を求めるために必要な仮の局所音速ごとの、仮想受信波と仮想合成受信波との誤差が算出される。つまり、全ての仮の局所音速について仮想受信波と仮想合成受信波との誤差が算出され出力される。誤差プロファイル生成部56には、Bモード画像の画素ごとに全ての仮の局所音速について仮想受信波と仮想合成受信波との誤差が入力される。入力された仮想受信波と仮想合成受信波との誤差は、横軸を仮の局所音速、縦軸を仮想受信波と仮想合成受信波との誤差としたグラフにプロットされ、誤差プロファイルが生成され出力される。局所音速値決定部58には、Bモード画像の画素ごとに求められた誤差プロファイルが入力される。入力された誤差プロファイルに基づき、画素ごとの局所音速値が決定され出力される。ここで、局所音速値の演算処理について説明する。図16は、局所音速値の演算処理を模式的に示す図である。図16bに示すように、被検体OBJ内の着目領域ROIを代表する格子点をXROI、格子点XROIよりも浅い即ち、超音波トランスデューサ18に近い位置にXY方向に等間隔で配置された格子点をA1,A2,...とし、少なくとも格子点XROIと各格子点A1,A2,...との間の音速はそれぞれ一定と仮定する。本例では、格子点A1,A2,...からの受信波それぞれWA1,WA2,...のT及び遅延時間ΔTが既知として、格子点XROIと格子点A1,A2,...の位置関係から格子点XROIにおける局所音速値を求める。具体的には、ホイヘンスの原理により、格子点XROIからの受信波WXと格子点A1,A2,...からの受信波を仮想的に合成した受信波WSUMが一致することを利用する。ここで、格子点XROIにおける局所音速値を求めるときの演算に使用する格子点A1,A2,...の範囲及び個数は予め決めておく。ここで、局所音速値演算に使用する格子点の範囲が広いと局所音速値の誤差が大きくなり、狭いと仮想受信波との誤差が大きくなるため、格子点の範囲はこれらの兼ね合いで決める。格子点A1,A2,...のX方向の間隔は、分解能と処理時間の兼ね合いで決定される。格子点A1,A2,...のX方向の間隔は、一例で1mmから1cmである。格子点A1,A2,...のY方向の間隔が狭いと誤差計算における誤差が大きくなり、広いと局所音速値の誤差が大きくなる。格子点A1,A2,...のY方向の間隔は、超音波画像の画像分解能の設定に基づいて決定される。格子点A1,A2,...のY方向の間隔は、一例で1cmである。なお、格子点A1,A2,...の間隔が広い場合、合成波の演算が困難になるため、補間によって細かい格子点を生成するようにすればよい。第2実施形態の超音波画像生成装置50の動作は、第1実施形態の超音波画像生成装置10の場合とほぼ同様であるから、繰り返しの説明は省略し、異なる部分のみ図17に示すフローチャートを参照して説明する。図17は、図3のステップS26とステップS32との間に追加されるフローチャートである。ステップS26までによって着目領域の全ての画素について環境音速値が求められる。続いて、局所音速値の算出を開始する開始着目画素例えば、x=n,y=mが設定されステップS40、着目画素の局所音速値の算出が行われるステップS42。ここで、着目画素の局所音速値の算出の詳細を、図18のフローチャートにより説明する。まず、格子点XROIにおける環境音速値に基づいて、格子点XROIを反射点とした時の仮想的な受信波WXの波形が算出されるステップS302。次に、格子点XROIにおける仮定音速の初期値が設定されるステップS304。そして、仮定音速が1ステップ変更されてステップS306、仮想的な合成受信波WSUMが算出されるステップS308。格子点XROIにおける局所音速値をVと仮定すると、格子点XROIから伝播した超音波が格子点A1,A2,...に到達するまでの時間はXROIA1/V,XROIA2/V,...となる。ここで、XROIA1,XROIA2,...は、それぞれ格子点A1,A2,...と格子点XROIとの間の距離である。格子点A1,A2,...における環境音速値は図3のステップS26までにより既知のため、各格子点A1,A2,...からの受信波は予め求めることができる。従って、格子点A1,A2,...からそれぞれ遅延XROIA1/V,XROIA2/V,...で発した反射波超音波エコーを合成することにより、仮想合成受信波WSUMを求めることができる。なお、実際には、素子データRF信号上で上記処理を行うため、格子点XROIから格子点A1,A2,...に到達するまでの時間それぞれT1,T2,...は下記の式2により表される。ここで、XA1,XA2,...は、それぞれ格子点A1,A2,...と格子点Xとの間のスキャン方向X方向の距離である。また、Δtは格子点のY方向時間間隔である。上記T1,T2,...に、格子点XROIと同音線の格子点Anから格子点XROIに到達するまでの時間Δt/2を足した遅延で各格子点A1,A2,...からの受信波を合成することにより、仮想合成受信波WSUMを求めることができる。ここで、格子点をY方向に時間軸で等間隔Δtに設定する場合、空間上での間隔は必ずしも等間隔にはならない。従って、各格子点に超音波が到達するまでの時間を計算するときに、式2においてΔt/2の代わりに補正したΔt/2を用いてもよい。ここで、補正したΔt/2は、例えば、格子点XROIと同音線の格子点Anに比べたA1,A2,...の深さY方向の距離の差をVで除算した値をΔt/2から加算減算した値である。各格子点A1,A2,...の深さはそれより浅い格子点において局所音速値が既知であることから求められる。また、仮想合成受信波WSUMの算出は、実際に格子点A1,A2,...から遅延XROIA1/V,XROIA2/V,...で発した既定のパルス波それぞれWA1,WA2,...を重ね合わせることにより行う。次に、仮想受信波WXと仮想合成受信波WSUMの誤差が算出されるステップS310。仮想受信波WXと仮想合成受信波WSUMの誤差は、互いの相互相関をとる方法、仮想受信波WXに仮想合成受信波WSUMから得られる遅延を掛けて位相整合加算する方法、または逆に仮想合成受信波WSUMに仮想受信波WXから得られる遅延を掛けて位相整合加算する方法により算出される。ここで、仮想受信波WXから遅延を得るには、格子点XROIを反射点とし、音速Vで伝播した超音波が各素子に到着する時刻を遅延とすればよい。また、仮想合成受信波WSUMから遅延を得るには、隣り合う素子間での合成受信波の位相差から等位相線を抽出し、その等位相線を遅延とするか、または単に各素子の合成受信波の最大ピーク位置の位相差を遅延としてもよい。また、各素子からの合成受信波の相互相関ピーク位置を遅延としてもよい。位相整合加算時の誤差は、整合加算後の波形のpeaktopeakとする方法、または包絡線検波した後の振幅の最大値とする方法により求められる。次に、ステップS306からS310が繰り返されて、全ての仮定音速の値での演算が終了するとステップS312でY、格子点XROIにおける局所音速値が判定されるステップS314。ホイヘンスの原理を厳密に適用した場合、上記ステップS308において求めた仮想合成受信波WSUMの波形は、格子点XROIにおける局所音速値をVと仮定した場合の仮想受信波反射波WXの波形と等しくなる。ステップS314では、仮想受信波WXと仮想合成受信波WSUMとの差が最小になる仮定音速の値を格子点XROIにおける局所音速値と判定する。なお、上記の方法仮想合成受信波形算出、仮想受信波形との誤差算出、音速判定の代わりに、格子点XROIの環境音速値と格子点A1,A2,...の環境音速値を入力として格子点XROIにおける音速値を出力とするテーブルを利用してもよい。また、異なる間隔、異なる範囲の格子点を用いて、局所音速値の判定を複数回行うようにしてもよい。着目画素の局所音速値の算出が終了すると、つまりステップS42が終了すると、着目画素のy座標の値がymaxと比較されステップS44、yの値がymax未満であるときはステップS44でN、yに1が足されてステップS46、ステップS42の着目画素の局所音速値の算出に戻る。着目画素の局所音速値の算出ステップS42は、y=ymaxとなるまで繰り返される。y=ymaxとなるとステップS44でY、着目画素のx座標の値がxmaxと比較されステップS48、xの値がxmax未満であるときはステップS48でN、xに1が足されてステップS50、y座標の値がyminy=mに設定されステップS52、ステップS42の着目画素の局所音速値の算出に戻る。つまり、y座標方向をラインとすると、x座標がnである1ライン目の局所音速値が算出されると、x座標が1だけインクリメントされn1、2ライン目の局所音速値が算出される。着目画素の局所音速値の算出ステップS42は、着目領域全体についてx=xmaxx=N,y=ymaxy=Mまで、局所音速値が算出されるまで繰り返される。着目領域全体について局所音速値の算出が終了すると、第1実施形態のステップS32以降と同様に、局所音速画像データおよび信頼度画像データが生成され、表示部40に表示される。なお、表示部40に表示される表示画像データは、第1実施形態と同様に、各種の表示画像データとすることができる。このように、環境音速値だけでなく局所音速値を用いる場合であっても、局所音速値を測定するための、専用の超音波を送受信する構成を用いることなく、より被検体内の病変の検出が行いやすい、局所音速画像、および信頼度画像を提供することができる。なお、上記の各実施形態では、Bモード画像を生成してその画素を着目画素としたが、Bモード画像を生成せず、RFデータに基づいて着目領域を複数の画素に分解し、当該画素を着目画素として音速画像が生成されるようにしてもよい。また、上記の各実施形態では、環境音速値または局所音速値をBモード画像の各画素に割り当てて音速画像を生成したが、これに限定されず、音速画像の画素はBモード画像の画素と1対1で対応しなくてもよい。例えば、Bモード画像の4画素分を音速画像の1画素としてもよい。また、上記の各実施形態では、通常観察時ライブモードの動作について説明したが、RFデータ記録再生部42に記録されたRFデータに基づいて、Bモード画像、音速画像、および信頼度画像が生成されるようにしてもよい。なお、本発明においては、上述した超音波画像生成方法の各工程をコンピュータに実行させるための超音波画像生成プログラムとして構成しても良いし、また、コンピュータを、超音波画像生成方法の各工程を実施する各手段として、または、上述した超音波画像生成装置を構成する各手段として機能させる超音波画像生成プログラムとして構成しても良い。また、本発明を、上述した超音波画像生成プログラムをコンピュータにより読取可能な媒体またはコンピュータにより読取可能なメモリとして構成してもよい。以上、本発明に係る超音波画像生成装置、超音波画像生成方法およびプログラムについて詳細に説明したが、本発明は、上記実施形態に限定されるものではなく、本発明の要旨を逸脱しない範囲において、各種の改良や変更を行ってもよい。"