JP2020124315A - 超音波断層像生成方法、超音波断層像生成装置、およびプログラム - Google Patents

Abstract

【課題】超音波断層像の空間分解能およびコントラストを、エコー信号間の相関性に基づく手法よりも向上させる技術を提供する。【解決手段】Ｍ（２以上の自然数）個の超音波振動子から発せられた超音波のエコーを受信してエコー信号を出力する超音波プローブから出力されるＭチャネルのエコー信号における雑音を推定し、受信焦点からのエコーを強調する重み係数を前記Ｍチャネルのエコー信号における信号対雑音比に応じて算出する推定ステップＳＡ１００と、超音波断層像を表すビームフォーマを、推定ステップＳＡ１００にて算出された重み係数を用いて前記Ｍチャネルのエコー信号から生成する生成ステップＳＡ１１０と、を含む超音波断層像生成方法を提供する。【選択図】図２

Description

本発明は、超音波断層像の生成技術に関する。

超音波診断は、被検者の体内の断層像を非侵襲的に測定する手法であり，広く医療現場に普及している。超音波診断により得られる超音波断層像の空間分解能およびコントラストは，診断精度に直結する重要な要素である。このため、超音波断層像の空間分解能およびコントラストを向上させる技術が種々提案されており、その一例としては非特許文献１に開示の技術が挙げられる。非特許文献１に開示の技術では、複数の超音波振動子からなる配列型超音波振動子により受信された各エコー信号間の相関性（coherence）に基づき，超音波断層像の方位分解能およびコントラストを向上させている。

超音波受信信号間の相関性に基づく超音波イメージング手法 P.-C. Li and M.-L. Li, "Adaptive imaging using the generalized coherence factor," IEEE Trans. Ultrason. Ferroelectr. Freq. Control, vol. 50, no. 2, pp. 128-141, 2003. H. Hasegawa and H. Kanai, "Effect of element directivity on adaptive beamforming applied to high-frame-rate ultrasound," IEEE Trans. Ultrason. Ferroelectr. Freq. Control, vol. 62, no. 3, pp. 511-523, 2015.

前述したように、超音波断層像の空間分解能およびコントラストは，診断精度に直結するため、高ければ高い程好ましい。

本発明は上述した事情に鑑みてなされたものであり、超音波断層像の空間分解能およびコントラストを、エコー信号間の相関性に基づく手法よりも向上させる技術を提供することを目的とする。

上記課題を解決するために本発明は、Ｍ（２以上の自然数）個の超音波振動子から発せられた超音波のエコーを受信してエコー信号を出力する超音波プローブから出力されるＭチャネルのエコー信号における雑音を推定し、受信焦点からのエコーを強調する重み係数を前記Ｍチャネルのエコー信号における信号対雑音比に応じて算出する推定ステップと、
超音波断層像を表すビームフォーマを、前記推定ステップにて算出された重み係数を用いて前記Ｍチャネルのエコー信号から生成する生成ステップと、を含むことを特徴とする超音波断層像生成方法、を提供する。

詳細については後述するが、本発明によれば、超音波断層像の空間分解能およびコントラストを、エコー信号間の相関性に基づく手法よりも向上させることが可能になる。

より好ましい態様の超音波断層像生成方法においては、前記推定ステップでは、前記Ｍチャネルのエコー信号に遅延和ビームフォーミングにおける遅延を付与して得られるｍ番目のエコー信号ｓ_ｍまでを累算して得られる累積要素信号ｕ_ｍと、エコー信号ｓ_ｍに含まれる直流成分ｙと付加的ノイズに起因するバイアスｎとを用いてモデル化したモデル化要素信号Ｕ_ｍ＝ｍ×ｙ＋ｎと、の平均二乗差αが最小になるようにｙおよびｎの値を設定し、設定したｙおよびｎを用いて前記平均二乗差αの最小値を算出し、当該最小値と前記設定したｙとから前記重み係数を算出し、前記生成ステップでは、前記設定したｙに前記重み係数を乗算して、超音波断層像を表すビームフォーマを生成する、ことを特徴としてもよい。

より好ましい態様の超音波断層像生成方法においては、前記推定ステップでは、前記Ｍチャネルのエコー信号に遅延和ビームフォーミングにおける遅延を付与して得られるｍ番目のエコー信号ｓ_ｍまでに含まれる雑音の積分値ｎ_ｍの二乗平均と、遅延補償後のエコー信号ｓ_ｍの平均Ｙ_DASの二乗平均とから前記重み係数を算出し、前記生成ステップでは、前記平均Ｙ_DASに前記重み係数を乗算して、超音波断層像を表すビームフォーマを生成する、ことを特徴としてもよい。

また、上記課題を解決するために本発明は、Ｍ（２以上の自然数）個の超音波振動子を有し、各超音波振動子から発せられた超音波のエコーを受信してエコー信号を出力する超音波プローブから出力されるＭチャネルのエコー信号における信号対雑音比を推定し、受信焦点からのエコーを強調する重み係数を前記Ｍチャネルのエコー信号における信号対雑音比に応じて算出する推定手段と、超音波断層像を表すビームフォーマを、前記推定手段により算出された重み係数を用いて前記Ｍチャネルのエコー信号から生成する生成手段と、を有することを特徴とする超音波断層像生成装置、を提供する。

また、上記課題を解決するために本発明は、コンピュータを、Ｍ（２以上の自然数）個の超音波振動子を有し、各超音波振動子から発せられた超音波のエコーを受信してエコー信号を出力する超音波プローブから出力されるＭチャネルのエコー信号における信号対雑音比を推定し、受信焦点からのエコーを強調する重み係数を前記Ｍチャネルのエコー信号における信号対雑音比に応じて算出する推定手段と、超音波断層像を表すビームフォーマを、前記推定手段により算出された重み係数を用いて前記Ｍチャネルのエコー信号から生成する生成手段と、して機能させることを特徴とするプログラム、を提供する。

本発明の一実施形態による超音波断層像生成装置２０を含む超音波医用システム１の構成例を示すブロック図である。 超音波断層像生成装置２０の信号処理部２３０が実行する信号処理の流れを示すフローチャートである。 超音波断層像の空間分解能を評価するための点ターゲットをイメージングした図である。 超音波断層像のコントラストを評価するためのファントムを画像化した図である。

以下、図面を参照しつつ本発明の実施形態を説明する。
（Ａ．実施形態）
図１は、本発明の一実施形態によれる超音波断層像生成装置２０を含む超音波医用システム１の構成例を示す図である。超音波医用システム１は、医療現場において被検者の体内の超音波断層像を非侵襲的に撮像するためのシステムである。図１に示すように、超音波医用システム１は、超音波断層像生成装置２０の他に、信号線を介して各々超音波断層像生成装置２０に接続された超音波プローブ１０、操作装置３０および表示装置４０、を有する。

超音波プローブ１０は、複数の超音波振動子よりなる配列型超音波振動子を有する。本実施形態の超音波医用システム１では、超音波プローブ１０として、Ｍ（２以上の自然数）個の超音波振動子を０．１ｍｍ間隔で配列したリニアアレイプローブ（ＰＵ−０５５８：上田日本無線株式会社）が用いられている。複数の超音波素子の各々は、超音波断層像生成装置２０による制御の下、被検者の検査部位に向けて超音波を放射するとともに、当該超音波のエコーを受信してエコー信号を出力する。

超音波断層像生成装置２０は、超音波プローブ１０に超音波を送信させるとともに、超音波プローブ１０からの出力信号に信号処理を施して画像データを生成する。操作装置３０は、例えばマウスなどのポインティングデバイスやキーボードを含む。操作装置３０は、超音波医用システム１の利用者（例えば、超音波診断のための各種操作を行う検査技師）に超音波断層像生成装置２０に対する各種入力操作を行わせるための装置である。表示装置４０は、例えば液晶ディスプレイである。表示装置４０は、超音波断層像生成装置２０の出力する画像データに応じて画像を表示する。

超音波断層像生成装置２０は、図１に示すように、制御部２００と、送信部２１０と、受信部２２０と、信号処理部２３０と、を有する。図１では詳細な図示を省略したが、超音波断層像生成装置２０は、ＯＳ（Operating System）などの各種ソフトウェアを記憶した記憶部（例えば、ハードディスク）も有する。

制御部２００は例えばＣＰＵ（Central Processing Unit）である。制御部２００は、上記記憶部に記憶されているソフトウェアを実行することにより、超音波断層像生成装置２０の制御中枢として機能し、各部の作動制御を行う。より詳細に説明すると、制御部２００は、従来と同様のライン毎の取得シーケンスによって超音波断層像が生成されるように、各部の作動制御を行う。

送信部２１０および受信部２２０には、信号線を介して超音波プローブ１０が接続されている。送信部２１０は、制御部２００から与えられる送信データにＤ／Ａ変換を施して送信信号を生成し、超音波プローブ１０が備えるＭ個の超音波振動子の各々に与える。これにより、超音波プローブ１０が備えるＭ個の超音波振動子の各々は超音波を放射する。受信部２２０は、超音波プローブ１０の複数の超音波振動子の各々から出力されるエコー信号にＡ／Ｄ変換を施し、さらに遅延を付与して遅延補償し、信号処理部２３０に与える。なお、本実施形態において、受信部２２０がエコー信号に付与する遅延は、従来の超音波断層像生成手法である遅延和ビームフォーミング（以下、ＤＡＳビームフォーミング）に準拠した遅延である。

超音波プローブ１０のｍ（ｍ＝０、１、２・・・Ｍ−１）番目の超音波振動子から出力され、受信部２２０により遅延を付与されたエコー信号をｓ_ｍとすると、超音波プローブ１０の受信開口に含まれるＭ個の超音波振動子により得られるエコー信号は以下の数１に示すベクトルＳで表される。遅延補償の後、ベクトルＳに含まれる受信焦点からのエコーは受信開口を横切る直流（ＤＣ）成分となる。したがって、従来のＤＡＳビームフォーミングでは、受信焦点ｙからのエコーに相当するビームフォーマ（すなわち、超音波断層像を表すビームフォーマ）Ｙ_DASは、遅延補償後のエコー信号ｓ_ｍの平均として以下の数２のように求められていた。

これに対して、信号処理部２３０は、遅延補償後のエコー信号ｓ_ｍに対して、本実施形態の特徴を顕著に示す信号処理（信号対雑音比に基づくビームフォーミング処理）を施して超音波断層像を表すビームフォーマを生成して表示装置４０に与える。信号処理部２３０は、例えばＤＳＰ（Digital Signal Processor）であり、図１では詳細な図示を省略したが、信号処理部２３０には、信号対雑音比に基づくビームフォーミング処理を当該信号処理部２３０に実行させる信号処理プログラムが予めインストールされている。信号処理部２３０は、受信部２２０により遅延を付与された信号に対して、上記信号処理プログラムにしたがって、信号対雑音比ビームフォーミングまたは線形回帰ビームフォーミングを実行する。信号処理部２３０により実行されるビームフォーミング処理としては、信号対雑音比ビームフォーミングと線形回帰ビームフォーミングの２種類がある。信号対雑音比ビームフォーミングおよび線形回帰ビームフォーミングは、広義にはどちらも信号対雑音比に基づくものであるが、２つの手法を区別するために異なる名称とした。以下、本実施形態の特徴を顕著に示す信号対雑音比ビームフォーミングおよび線形回帰ビームフォーミングについて説明する。

図２は、信号対雑音比ビームフォーミングおよび線形回帰ビームフォーミングの流れを示すフローチャートである。図２に示すように、両手法には、推定ステップＳＡ１００と、推定ステップＳＡ１００に後続する生成ステップＳＡ１１０の２つのステップが含まれる。つまり、信号処理プログラムにしたがって作動している信号処理部２３０は、図１に示すように、推定ステップＳＡ１００を実行する推定手段２３０ａ、および生成ステップＳＡ１１０を実行する生成手段２３０ｂとして機能する。

信号対雑音比ビームフォーミングにおける推定ステップＳＡ１００では、信号処理部２３０は、受信部２２０から出力されるＭチャネルのエコー信号における信号対雑音比を推定し、受信焦点からのエコーを強調する重み係数（信号対雑音比に応じた重み係数）を算出する。上述したように、受信焦点からのエコーyは、遅延補償後のエコー信号ｓ_mの直流成分となる。信号対雑音比ビームフォーミングにおける推定ステップＳＡ１００では、信号処理部２３０は、信号成分と雑音成分を遅延補償後のエコー信号ｓ_mの平均値と分散で推定し、受信焦点からのエコーを強調する重み係数Ｗ_ＳＮＲを以下の数３にしたがって算出する。そして、信号対雑音比ビームフォーミングにおける生成ステップＳＡ１１０では、信号処理部２３０は、信号対雑音比ビームフォーミングの出力（すなわち、超音波断層像を表すビームフォーマ）Ｙ_SNRを以下の数４にしたがって算出して表示装置４０に与える。

遅延補償後のエコー信号ｓ_mの信号対雑音比が非常に高い場合は、数３の分母が非常に小さくなり、その結果Ｗ_ＳＮＲが極端に大きくなってビームフォーマ出力が不安定となる。それを避けるため、安定化パラメータβ(実数)を数５のように導入してもよい。
βが０に近いほど数５の分母が小さくなることが回避され、ビームフォーマ出力は安定するが、空間分解能などの改善効果は低下する。安定化パラメータβの値については、ビームフォーマ出力の安定度と空間分解能等の改善効果との兼ね合いで適宜好適な値に定めるようにすればよい。以上が信号対雑音比ビームフォーミングの内容である。

次いで、線形回帰ビームフォーミングについて説明する。
線形回帰ビームフォーミングにおける推定ステップＳＡ１００では、信号処理部２３０は、信号対雑音比ビームフォーミングにおける処理とは異なる処理で、受信部２２０から出力されるＭチャネルのエコー信号における雑音を推定し、受信焦点からのエコーを強調する重み係数を算出する。より詳細に説明すると、信号処理部２３０は、まず、累積要素信号ｕ_ｍを以下の数６にしたがって算出する（ただし、ｕ_０＝０）。なお、数６の右辺におけるｓ_ｉはｉ番目の超音波振動子からの遅延補償後のエコー信号である。

上述したように、受信焦点からのエコーyは、遅延補償後のエコー信号ｓ_mの直流成分となる。したがって、累積要素信号u_ｍは、以下の数７のような線形関数モデル化される。なお、数７におけるｎは付加的ノイズに起因するバイアスである。以下では、数７にしたがってモデル化した信号をモデル化要素信号と呼ぶ。測定された累積要素信号u_ｍとモデル化要素信号Ｕ_ｍとの平均二乗差αは以下の数８のように定義され、信号処理部２３０は、数８で定義される平均二乗差αが最小になるように、ｙおよびｎの値（以下、最小二乗推定値）を設定（すなわち、信号対雑音比を推定）する。なお、ｙおよびｎの最小二乗推定値は、数９に示すように、yとnに対するαの偏微分をゼロに設定することによって得られる。

次いで、信号処理部２３０は、まず、数９にしたがって算出した最小二乗推定値ＹおよびＮを、数８におけるｙおよびｎに代入して平均二乗差αの最小値α_ｍｉｎを算出する。そして、信号処理部２３０は、受信焦点からのエコーを強調する重み係数Ｗ_ＬＲを以下の数１０にしたがって算出し、線形回帰ビームフォーミングにおける推定ステップＳＡ１００を終了する。線形回帰ビームフォーミングにおける生成ステップＳＡ１１０では、線形回帰ビームフォーマの出力（すなわち、超音波断層像を表すビームフォーマ出力）Ｙ_ＬRを以下の数１１にしたがって算出して表示装置４０に与える。

信号対雑音比ビームフォーミングと同様に、遅延補償後のエコー信号ｓ_mの信号対雑音比が非常に高い場合は、数１０の分母が非常に小さくなり、その結果Ｗ_ＬＲが極端に大きくなってビームフォーマ出力が不安定となる。それを避けるため、安定化パラメータγ(実数)を数１２のように導入してもよい。
γの値を大きくするほど、ビームフォーマ出力は安定するが、空間分解能などの改善効果は低減する。安定化パラメータγの値についても、前述のβと同様に、ビームフォーマ出力の安定度と空間分解能等の改善効果との兼ね合いで適宜好適な値に定めるようにすればよい。以上が線形回帰ビームフォーミングの内容である。

上述した線形回帰ビームフォーミングの推定ステップＳＡ１００では、信号対雑音比の推定に最小二乗法が用いられているため、信号対雑音比ビームフォーミングに比べ計算負荷が高い。そこで、線形回帰ビームフォーミングの計算効率を向上させる（すなわち、計算負荷を低減させる）ため、以下のように変形してもよい。
計算効率を向上させた線形回帰ビームフォーミングの推定ステップＳＡ１００では、信号処理部２３０は、ｍ番目の素子による受信信号ｓ_mに含まれる雑音の積分値ｎ_mを数１３により算出する。
計算効率を向上させた線形回帰ビームフォーミングにおける重み係数Ｗ_ＬＲeは、数１３により得られる雑音成分の積分値ｎ_mを用いて、以下の数１４のように定義される。計算効率を向上させた線形回帰ビームフォーミングの推定ステップＳＡ１００では、信号処理部２３０は、数１４にしたがって重み係数Ｗ_ＬＲeを算出する。なお、数１４におけるγは数１２におけるものと同様に安定化パラメータである。
計算効率を向上させた線形回帰ビームフォーミングの生成ステップＳＡ１１０では、信号処理部２３０は、超音波断層像を表すビームフォーマ出力Ｙ_ＬReを以下の数１５にしたがって算出し、表示装置４０に与える。

また、信号対雑音比ビームフォーミング、線形回帰ビームフォーミング、計算効率を向上させた線形回帰ビームフォーミングの何れについても，非特許文献２に示された開口分割処理を組み合わせることで、計算量を低減させてもよい。

図３は、超音波断層像の空間分解能を評価するための点ターゲットをイメージングした結果である。より詳細には、図３（ａ）はＤＡＳビームフォーミングにより得られた画像，図３（ｂ）は受信信号間の相関性に基づく手法により得られた画像、図３（ｃ）と（ｄ）はそれぞれ、本実施形態の信号対雑音比ビームフォーミングおよび線形回帰ビームフォーミングにより得られた画像である。図３（ａ）〜（ｄ）の各図において、画像の輝度（白色の強さ）は超音波散乱波の強度を示す。図３（ａ）〜（ｄ）に示す各画像を比較すれば明らかなように、本実施形態の信号対雑音比ビームフォーミングおよび線形回帰ビームフォーミングにより得られた画像（図３（ｃ）、（ｄ））では、白の輝点のサイズが図３（ａ）および図３（ｂ）に比較して小さくなっている。このことから，本実施形態の信号対雑音比ビームフォーミングおよび線形回帰ビームフォーミングによれば、ＤＡＳビームフォーミングおよび受信信号管の相関性に基づく手法に比較して、高い空間分解能が得られることが判る。

図４は、超音波断層像のコントラストを評価するためのファントム（虚像）を画像化したものである。より詳細に説明すると、図４（ａ）は、従来のＤＡＳビームフォーミングにより得られた画像、図４（ｂ）は、受信信号間の相関性により得られた画像、図４（ｃ）と図４（ｄ）はそれぞれ、本実施形態の信号対雑音比ビームフォーミングおよび線形回帰ビームフォーミングにより得られた画像である。なお、図４（ａ）〜図４（ｄ）の各図においては、中央部の暗い部分は超音波散乱波が発生しない媒質（具体的には、嚢胞模擬部）であり，黒一色で描出されることが望ましい。図４（ａ）および図４（ｂ）の各画像では，嚢胞模擬部においても白い輝点が発生しており，これらの白い輝点は虚像である。図４（ｃ）では虚像が低減されていることがわかる。さらに、図４（ｄ）に示す画像では嚢胞模擬部において虚像が描出されていない。つまり、本実施形態の信号対雑音比ビームフォーミングおよび線形回帰ビームフォーミングによれば、従来のＤＡＳビームフォーミングおよび受信信号管の相関性に基づく手法に比較して、嚢胞模擬部における虚像の描出を抑制することができ、コントラストが向上していることが判る。

図４において、虚像の抑圧効果が信号対雑音比ビームフォーミングに比べ線形回帰ビームフォーミングの方が高かった理由は、数６の処理による効果、すなわち、遅延補償後のエコー信号ｓ_ｍの積分効果である。積分は低域通過フィルタに対応する。積分操作以外のフィルタを適用することで、線形回帰ビームフォーマの出力をさらに向上させることも可能である。なお、数１３の積分処理についても同様に他のフィルタ処理に適宜変更してもよい。

以上説明したように、本発明によれば、超音波断層像の空間分解能およびコントラストを、従来のＤＡＳビームフォーミングは勿論、超音波受信信号間の相関性に基づく手法と比較しても、さらに向上させることが可能になる。

（Ｂ．変形）
以上、本発明の一実施形態について説明したが、この実施形態に以下の変形を加えても勿論よい。
（１）上記実施形態では、超音波医用システムへの本発明の適用例を説明したが、物体の非破壊検査など、医療用以外の超音波断層像の生成への適用も可能である。医療用以外の技術分野においても、超音波断層像の空間分解能およびコントラストは，高ければ高い程好ましいからである。

（２）線形回帰ビームフォーミングの推定ステップＳＡ１００では、Ｍチャネルのエコー信号の信号対雑音比を最小二乗法により推定したが、尤度を用いる手法など他の手法により信号対雑音比を推定してもよい。要は、Ｍ（２以上の自然数）個の超音波振動子から発せられた超音波のエコーを受信してエコー信号を出力する超音波プローブから出力されるＭチャネルのエコー信号における雑音を推定し、受信焦点からのエコーを強調する重み係数を前記Ｍチャネルのエコー信号における信号対雑音比に応じて算出する推定ステップと、超音波断層像を表すビームフォーマを、前記推定ステップにて算出された重み係数を用いて前記Ｍチャネルのエコー信号から生成する生成ステップと、を含むことを特徴とする超音波断層像生成方法であればよい。

（３）上記実施形態では、超音波断層像生成装置２０の信号処理部２３０が推定手段２３０ａおよび生成手段２３０ｂとして機能したが、制御部２００を推定手段２３０ａおよび生成手段２３０ｂとして機能させてもよい。具体的には、受信部２２０の出力信号を制御部２００に与え、制御部２００に上記実施形態の信号処理プログラムを実行させるようにすればよい。

（４）上記実施形態では、本実施形態の特徴を顕著に示す超音波断層像生成方法を実現する信号処理プログラムが超音波断層像生成装置２０に予めインストールされていた。しかし、ＣＰＵなどのコンピュータを、Ｍ（２以上の自然数）個の超音波振動子を有し、各超音波振動子から発せられた超音波のエコーを受信してエコー信号を出力する超音波プローブから出力されるＭチャネルのエコー信号における雑音を推定し、受信焦点からのエコーを強調する重み係数を前記Ｍチャネルのエコー信号における信号対雑音比に応じて算出する推定手段と、超音波断層像を表すビームフォーマを、前記推定手段により算出された重み係数を用いて前記Ｍチャネルのエコー信号から生成する生成手段と、して機能させるプログラムを単体で製造し、販売等の配布をしてもよい。上記プログラムの配布態様の具体例としては、インターネットなどの電気通信回線経由のダウンロードにより配布する態様や、ＣＤ−ＲＯＭ（Compact Disk-Read Only Memory）やフラッシュＲＯＭ（Read Only Memory）などのコンピュータ読み取り可能な記録媒体に書き込んで配布する態様が挙げられる。このようにして配布される上記プログラムにしたがって、コンピュータを作動させることで、当該コンピュータに、本発明の超音波断層像生成方法を実行させることが可能になるからである。

（５）上記実施形態では、本実施形態の特徴を顕著に示す超音波断層像生成方法の各ステップを実行する推定手段２３０ａおよび生成手段２３０ｂがソフトウェアモジュールとして実現されていた。しかし、Ｍ（２以上の自然数）個の超音波振動子を有し、各超音波振動子から発せられた超音波のエコーを受信してエコー信号を出力する超音波プローブから出力されるＭチャネルのエコー信号における雑音を推定し、受信焦点からのエコーを強調する重み係数を前記Ｍチャネルのエコー信号における信号対雑音比に応じて算出する推定手段と、超音波断層像を表すビームフォーマを、前記推定手段により算出された重み係数を用いて前記Ｍチャネルのエコー信号から生成する生成手段の各々をＡＳＩＣなどの電子回路で構成し、これら電子回路を組み合わせて本発明の超音波断層像生成を構成してもよい。

１…超音波医用システム、１０…超音波プローブ、２０…超音波断層像生成装置、３０…操作装置、４０…表示装置、２００…制御部、２１０…送信部、２２０…受信部、２３０…信号処理部、２３０ａ…推定手段、２３０ｂ…生成手段。

Claims (5)

1. Ｍ（２以上の自然数）個の超音波振動子から発せられた超音波のエコーを受信してエコー信号を出力する超音波プローブから出力されるＭチャネルのエコー信号における雑音を推定し、受信焦点からのエコーを強調する重み係数を前記Ｍチャネルのエコー信号における信号対雑音比に応じて算出する推定ステップと、
   超音波断層像を表すビームフォーマを、前記推定ステップにて算出された重み係数を用いて前記Ｍチャネルのエコー信号から生成する生成ステップと、
   を含むことを特徴とする超音波断層像生成方法。
2. 前記推定ステップでは、
   前記Ｍチャネルのエコー信号に遅延和ビームフォーミングにおける遅延を付与して得られるｍ番目のエコー信号ｓ_ｍまでを累算して得られる累積要素信号ｕ_ｍと、エコー信号ｓ_ｍに含まれる直流成分ｙと付加的ノイズに起因するバイアスｎとを用いてモデル化したモデル化要素信号Ｕ_ｍ＝ｍ×ｙ＋ｎと、の平均二乗差αが最小になるようにｙおよびｎの値を設定し、設定したｙおよびｎを用いて前記平均二乗差αの最小値を算出し、当該最小値と前記設定したｙとから前記重み係数を算出し、
   前記生成ステップでは、
   前記設定したｙに前記重み係数を乗算して、超音波断層像を表すビームフォーマを生成する
   ことを特徴とする請求項１に記載の超音波断層像生成方法。
3. 前記推定ステップでは、
   前記Ｍチャネルのエコー信号に遅延和ビームフォーミングにおける遅延を付与して得られるｍ番目のエコー信号ｓ_ｍまでに含まれる雑音の積分値ｎ_ｍの二乗平均と、遅延補償後のエコー信号ｓ_ｍの平均Ｙ_DASの二乗平均とから前記重み係数を算出し、
   前記生成ステップでは、
   前記平均Ｙ_DASに前記重み係数を乗算して、超音波断層像を表すビームフォーマを生成する
   ことを特徴とする請求項１に記載の超音波断層像生成方法。
4. Ｍ（２以上の自然数）個の超音波振動子を有し、各超音波振動子から発せられた超音波のエコーを受信してエコー信号を出力する超音波プローブから出力されるＭチャネルのエコー信号における信号対雑音比を推定し、受信焦点からのエコーを強調する重み係数を前記Ｍチャネルのエコー信号における信号対雑音比に応じて算出する推定手段と、
   超音波断層像を表すビームフォーマを、前記推定手段により算出された重み係数を用いて前記Ｍチャネルのエコー信号から生成する生成手段と、
   を有することを特徴とする超音波断層像生成装置。
5. コンピュータを、
   Ｍ（２以上の自然数）個の超音波振動子を有し、各超音波振動子から発せられた超音波のエコーを受信してエコー信号を出力する超音波プローブから出力されるＭチャネルのエコー信号における信号対雑音比を推定し、受信焦点からのエコーを強調する重み係数を前記Ｍチャネルのエコー信号における信号対雑音比に応じて算出する推定手段と、
   超音波断層像を表すビームフォーマを、前記推定手段により算出された重み係数を用いて前記Ｍチャネルのエコー信号から生成する生成手段と、
   して機能させることを特徴とするプログラム。