(57)【要約】【課題】合成開口法において、空間解像度、Ｓ／Ｎ比、および、フレームレートの低下を抑制しながら画像化領域を素子列方向に拡大する。【解決手段】合成開口法において、対応するフレーム音響線信号が形成されるべき被検体内の領域として画像化主領域と、前記画像化主領域に対して列方向に隣接する追加領域とを設定し、送信イベントごとに、被検体内で超音波が集束するエリアに含まれる領域を含む主対象領域と、前記主対象領域に列方向に隣接する領域の一部または全部が前記追加領域内に存在する場合に前記追加領域内の当該部分である副対象領域とを、観測点を含む領域として設定し、画像化主領域内かつ主対象領域内の観測点と副対象領域に含まれる観測点との間で、送信された超音波が観測点に到達する送信時間の算出方法を異ならせる。【選択図】図１０
2019072358.tif 000002
【特許請求の範囲】
【請求項１】
超音波プローブに列設された複数の振動子を選択的に駆動して被検体に超音波送信する送信イベントを複数回繰り返すとともに、各送信イベントに同期して被検体から反射超音波を受波し、受波した反射超音波に基づいて生成される複数のサブフレーム音響線信号からフレーム音響線信号を合成する超音波信号処理装置であって、
１回の送信イベントにおいて前記複数の振動子から送信振動子列を選択し、前記送信振動子列から超音波が被検体中で集束するように送信し、各送信イベントに同期して、超音波を送信する送信振動子列が列方向に順次移動するよう選択する送信部と、
各送信イベントに同期して前記複数の振動子から受信振動子列を選択し、前記受信振動子列が被検体内から受波した反射超音波に基づいて、前記受信振動子列の振動子各々に対する受信信号列を生成する受信部と、
対応するフレーム音響線信号が形成されるべき被検体内の領域として画像化主領域と、前記画像化主領域に対して列方向に隣接する追加領域とを設定する画像化領域設定部と、
前記送信イベントごとに、各観測点から得られた反射超音波に基づく前記受信信号列を整相加算して前記サブフレーム音響線信号を生成する整相加算部と、
前記整相加算部が生成した前記複数のサブフレーム音響線信号に基づき、前記フレーム音響線信号を合成する合成部と
を備え、
前記整相加算部は、前記被検体内で超音波が集束するエリアに含まれる領域を含む主対象領域と、前記主対象領域に列方向に隣接する領域の一部または全部が前記追加領域内に存在する場合に前記追加領域内の当該部分である副対象領域とを、前記観測点を含む領域として設定し、前記画像化主領域内かつ前記主対象領域内の観測点と前記副対象領域に含まれる観測点との間で、送信された超音波が観測点に到達する送信時間の算出方法を異ならせる
ことを特徴とする超音波信号処理装置。
【請求項２】
前記整相加算部は、前記第１主対象領域に含まれる観測点について、前記観測点の深さが前記フォーカス深さより深い位置の場合は、送信された超音波が前記主対象領域内であって前記フォーカス深さにある基準点に到達するまでの第１時間に、超音波が前記基準点から前記観測点に到達するまでの第２時間を加算して算出した到達時間を、前記観測点の深さが前記フォーカス深さ未満の場合は、前記第１時間から前記第２時間を減算して算出した到達時間を、前記観測点についての前記送信時間とする
ことを特徴とする請求項１に記載の超音波信号処理装置。
【請求項３】
前記整相加算部は、前記副対象領域に含まれる観測点について、送信された超音波が前記フォーカス点に最近接する前記送信振動子列上の点から前記観測点に到達するまでの時間を、前記観測点についての前記送信時間とする
ことを特徴とする請求項１または２に記載の超音波信号処理装置。
【請求項４】
前記整相加算部は、前記副対象領域に含まれる観測点について、送信された超音波が前記フォーカス点に最近接する前記送信振動子列上の点から前記観測点に到達するまでの時間をＴ1とし、前記観測点の深さが前記フォーカス深さより深い位置の場合は、送信された超音波が前記主対象領域に含まれ前記フォーカス深さにある基準点に到達するまでの第３時間に、超音波が前記基準点から前記観測点に到達するまでの第４時間を加算して算出した時間、前記観測点の深さが前記フォーカス深さ未満の場合は、前記第３時間から前記第４時間を減算して算出した時間をＴ2としたとき、前記送信時間ＴMを、数式
ＴM＝αＴ1＋（１−α）Ｔ2
を用いて算出し、αの値は、０より大きく１以下である
ことを特徴とする請求項１または２に記載の超音波信号処理装置。
【請求項５】
前記整相加算部は、前記副対象領域に含まれる観測点について、送信された超音波が前記フォーカス点に最近接する前記送信振動子列上の点から前記観測点に到達するまでの時間をＴ1とし、前記観測点の深さが前記フォーカス深さより深い位置の場合は、送信された超音波が前記主対象領域に含まれ前記フォーカス深さにある基準点に到達するまでの第３時間に、超音波が前記基準点から前記観測点に到達するまでの第４時間を加算して算出した時間、前記観測点の深さが前記フォーカス深さ未満の場合は、前記第３時間から前記第４時間を減算して算出した時間をＴ2とし、送信された超音波が前記超音波振動子列の列中心から前記観測点と同じ深さで且つ前記超音波振動子列の列中心に最近接する第２基準点に到達するまでの時間とＴ3としたとき、前記送信時間ＴMを、数式
ＴM＝αＴ1＋（１−α）βＴ2＋（１−α）（１−β）Ｔ3
を用いて算出し、αの値は、０より大きく１以下であり、βの値は、０以上１以下である
ことを特徴とする請求項１または２に記載の超音波信号処理装置。
【請求項６】
前記βの値は、対応する観測点の深さと前記フォーカス深さとの差異が大きくなるほど大きくなる
ことを特徴とする請求項５に記載の超音波信号処理装置。
【請求項７】
前記αの値は、対応する観測点の深さと前記フォーカス深さとの差異が大きくなるほど小さくなる
ことを特徴とする請求項５または６に記載の超音波信号処理装置。
【請求項８】
前記主対象領域は、前記フォーカス点に最近接する前記送信振動子列上の点と、前記フォーカス点とを結ぶ直線に対して線対称の形状を有する
ことを特徴とする請求項１から７のいずれか１項に記載の超音波信号処理装置。
【請求項９】
前記送信振動子列が前記超音波プローブの端部の振動子を含まない送信イベントにおいて、前記主対象領域は、前記フォーカス点と前記送信振動子列の一端とを結ぶ直線と、前記フォーカス点と前記送信振動子列の他端とを結ぶ直線との間に位置する領域である
ことを特徴とする請求項８に記載の超音波信号処理装置。
【請求項１０】
前記送信振動子列が前記超音波プローブの端部の振動子を含む送信イベントにおいて、前記振動子列を延長した仮想線上に、前記送信振動子列のうち前記超音波プローブの端部とは異なる端部を一方の端として、前記フォーカス点に最近接する前記送信振動子列上の点を中点とする仮想送信振動子列を仮定したとき、前記主対象領域は、前記フォーカス点と前記仮想送信振動子列の一端とを結ぶ直線と、前記フォーカス点と前記仮想送信振動子列の他端とを結ぶ直線との間に位置する領域である
ことを特徴とする請求項８または９に記載の超音波信号処理装置。
【請求項１１】
前記追加領域は、列方向における前記フォーカス点の位置が前記追加領域に最も近づく送信イベントにおいて前記仮想送信振動子列のうち前記送信振動子列に含まれない部分を底辺とする領域である
ことを特徴とする請求項１０に記載の超音波信号処理装置。
【請求項１２】
前記追加領域と前記画像化主領域との境界線と、前記フォーカス点と前記送信振動子列の一端とを結ぶ直線との列方向の距離をｍとしたとき、前記第１領域は、前記追加領域を列方向に沿ってｍだけシフトした領域から前記主対象領域との重複領域を除いた領域である
ことを特徴とする請求項１から１１のいずれか１項に記載の超音波信号処理装置。
【請求項１３】
前記整相加算部は、前記画像化主領域と前記追加領域とのそれぞれの範囲と、前記フォーカス点の位置によって定まる前記第１領域の範囲とに基づいて前記副対象領域が空領域であるか否かを判定するための情報を保持し、前記情報を用いて前記副対象領域が空領域であると判定した送信イベントにおいて、前記副対象領域を設定せず前記主対象領域のみを設定する
ことを特徴とする請求項１から１２のいずれか１項に記載の超音波信号処理装置。
【請求項１４】
前記合成部は、前記画像化主領域に含まれる観測点について、前記主対象領域に対応する前記サブフレーム音響線信号を観測点の位置に基づいて合成することでフレーム音響線信号を生成する
ことを特徴とする請求項１から１３のいずれか１項に記載の超音波信号処理装置。
【請求項１５】
超音波プローブと、
請求項１から１４のいずれか１項に記載の超音波信号処理装置と、
を備えることを特徴とする超音波診断装置。
【請求項１６】
超音波プローブに列設された複数の振動子を選択的に駆動して被検体に超音波送信する送信イベントを複数回繰り返すとともに、各送信イベントに同期して被検体から反射超音波を受波し、受波した反射超音波に基づいて生成される複数のサブフレーム音響線信号からフレーム音響線信号を合成する超音波信号処理方法であって、
１回の送信イベントにおいて前記複数の振動子から送信振動子列を選択し、前記送信振動子列から超音波が被検体中で集束するように送信し、各送信イベントに同期して、超音波を送信する送信振動子列が列方向に順次移動するよう選択し、
各送信イベントに同期して前記複数の振動子から受信振動子列を選択し、前記受信振動子列が被検体内から受波した反射超音波に基づいて、前記受信振動子列の振動子各々に対する受信信号列を生成し、
対応するフレーム音響線信号が形成されるべき被検体内の領域として画像化主領域と、前記画像化主領域に対して列方向に隣接する追加領域とを設定し、
前記送信イベントごとに、各観測点から得られた反射超音波に基づく前記受信信号列を整相加算して前記サブフレーム音響線信号を生成し、
前記整相加算部が生成した前記複数のサブフレーム音響線信号に基づき、前記フレーム音響線信号を合成する
超音波信号処理方法であって、
前記サブフレーム音響線信号を生成する際に、前記被検体内で超音波が集束するエリアに含まれる領域を含む主対象領域と、前記主対象領域に列方向に隣接する領域の一部または全部が前記追加領域内に存在する場合に前記追加領域内の当該部分である副対象領域とを、前記観測点を含む領域として設定し、前記画像化主領域内かつ前記主対象領域内の観測点と前記副対象領域に含まれる観測点との間で、送信された超音波が観測点に到達する送信時間の算出方法を異ならせる
ことを特徴とする超音波信号処理方法。
【請求項１７】
表示部が接続可能であって、請求項１６に記載の超音波信号処理方法を実行する超音波信号処理装置における超音波画像表示方法であって、
前記フレーム音響線信号を超音波画像に変換して前記表示部に表示する
超音波画像表示方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、超音波信号処理装置、および、それを備えた超音波診断装置に関し、特に、超音波信号処理装置における受信ビームフォーミング処理方法、および、超音波画像表示方法に関する。
【背景技術】
【０００２】
超音波診断装置は、超音波プローブ（以下、「プローブ」とする）により被検体内部に超音波を送信し、被検体組織の音響インピーダンスの差異により生じる超音波反射波（エコー）を受信する。さらに、この受信から得た受信信号に基づいて、被検体の内部組織の構造を示す超音波断層画像を生成し、モニタ（以下、「表示部」とする）上に表示するものである。超音波診断装置は、被検体への侵襲が少なく、リアルタイムに体内組織の状態を断層画像などで観察できるため、生体の形態診断に広く用いられている。
【０００３】
従来の超音波診断装置では、受信した反射超音波に基づく信号の受信ビームフォーミング方法として、一般的に整相加算法と呼ばれる方法が使用されている（例えば、非特許文献１）。この方法では、複数の振動子によって行われる被検体への超音波送信が行われる際、被検体のある深さで超音波ビームがフォーカスを結ぶよう送信ビームフォーミングがなされる。また、この方法では、送信超音波ビームの中心軸上に観測点を設定する。そのため、１回の超音波送信イベントでは送信超音波ビームの中心軸近傍の少数本の音響線信号しか生成することができず、超音波の利用効率が悪い。また、観測点がフォーカス点から離れた位置にある場合には、音響線信号の空間分解能及び信号Ｓ／Ｎ比が低くなる課題も有している。
【０００４】
これに対して、合成開口法(Synthetic Aperture Method)により、送信フォーカス点近傍以外の領域においても空間分解能の高い、高画質な画像を得る受信ビームフォーミング方法が考案されている（例えば、非特許文献２）。この方法によれば、超音波送信波の伝播経路と、その伝播経路による反射波の振動子への到達時間の両方を加味した遅延制御を行うことで、送信フォーカス点近傍以外に位置する超音波主照射領域からの反射超音波も反映した受信ビームフォーミングを行うことができる。その結果、１回の超音波送信イベントから超音波主照射領域全体に対して音響線信号を生成することができる。なお、超音波主照射領域とは、領域内のすべての点において、送信振動子列を構成する各振動子から送波される超音波の伝搬する領域を指す。また、合成開口法では、複数の送信イベントから得た同一観測点に対する複数の受信信号をもとに仮想的に送信フォーカスを合わせることで、非特許文献１記載の受信ビームフォーミング方法と比較して、空間分解能及びＳ／Ｎ比の高い超音波画像を得ることが可能となる。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平４−６０６５３号公報
【非特許文献】
【０００６】
【非特許文献１】伊東正安、望月剛共著「超音波診断装置」コロナ社出版、２００２年８月２６日（Ｐ４２−Ｐ４５）
【非特許文献２】"Virtual ultrasound sources in high resolution ultrasound imaging", S.I.Nikolov and J.A.Jensen, in Proc, SPIE − Progress in biomedical optics and imaging, vol. 3, 2002, P. 395−405
【発明の概要】
【発明が解決しようとする課題】
【０００７】
超音波画像を生成する対象となる領域（以下、「画像化領域」とする）をプローブにおける振動子の並ぶ方向（以下、「素子列方向」とする）に拡大することが検討されている。画像化領域を素子列方向に拡大する方法としては、例えば、特許文献１に開示されているように、（１）送信イベント毎に超音波ビームの送信方向を変更し、超音波ビームを放射状に送信することにより、超音波ビームの通過領域を拡大する方法、と、（２）通常の送信イベントに加えて、超音波ビームの送信方向を外側に向けた送信イベントを追加する方法、とがある。しかしながら、（１）の方法では、超音波ビームの通過領域が素子列方向に拡大する代償として、素子列方向における空間解像度の低下と、Ｓ／Ｎ比の低下を招くこととなる。また、（２）の方法では、１枚の超音波画像を生成するために必要な送信イベント回数の増加を伴うため、フレームレートの低下を招くこととなる。
【０００８】
本発明は、上記課題に鑑みてなされたものであり、一般的な送信ビームフォーミングを用いた合成開口法において、空間解像度、Ｓ／Ｎ比、および、フレームレートの低下を抑制しながら画像化領域を素子列方向に拡大できる超音波信号処理装置、超音波診断装置、超音波信号処理方法、および、超音波画像表示方法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
本発明の一態様に係る超音波信号処理装置は、超音波プローブに列設された複数の振動子を選択的に駆動して被検体に超音波送信する送信イベントを複数回繰り返すとともに、各送信イベントに同期して被検体から反射超音波を受波し、受波した反射超音波に基づいて生成される複数のサブフレーム音響線信号からフレーム音響線信号を合成する超音波信号処理装置であって、１回の送信イベントにおいて前記複数の振動子から送信振動子列を選択し、前記送信振動子列から超音波が被検体中で集束するように送信し、各送信イベントに同期して、超音波を送信する送信振動子列が列方向に順次移動するよう選択する送信部と、各送信イベントに同期して前記複数の振動子から受信振動子列を選択し、前記受信振動子列が被検体内から受波した反射超音波に基づいて、前記受信振動子列の振動子各々に対する受信信号列を生成する受信部と、対応するフレーム音響線信号が形成されるべき被検体内の領域として画像化主領域と、前記画像化主領域に対して列方向に隣接する追加領域とを設定する画像化領域設定部と、前記送信イベントごとに、前記送信イベントごとに、各観測点から得られた反射超音波に基づく前記受信信号列を整相加算して前記サブフレーム音響線信号を生成する整相加算部と、前記整相加算部が生成した前記複数のサブフレーム音響線信号に基づき、前記フレーム音響線信号を合成する合成部とを備え、前記整相加算部は、前記被検体内で超音波が集束するエリアに含まれる領域を含む主対象領域と、前記主対象領域に列方向に隣接する領域の一部または全部が前記追加領域内に存在する場合に前記追加領域内の当該部分である副対象領域とを、前記観測点を含む領域として設定し、前記画像化主領域内かつ前記主対象領域内の観測点と前記副対象領域に含まれる観測点との間で、送信された超音波が観測点に到達する送信時間の算出方法を異ならせることを特徴とする。
【発明の効果】
【００１０】
本発明の一形態に係る超音波信号処理装置、超音波診断装置、超音波信号処理方法、および、超音波画像表示方法によれば、フレームレートを低下させることなく、かつ、空間解像度とＳ／Ｎ比を低下させることなく、追加領域の面積だけ画像化領域を拡大することができる。
【図面の簡単な説明】
【００１１】
【図１】実施の形態に係る超音波診断装置１００の構成を示すブロック図である。
【図２】実施の形態に係る送信ビームフォーマ部１０３による送信超音波ビームの伝播経路を示す図である。
【図３】実施の形態に係る受信ビームフォーマ部１０４の構成を示す機能ブロック図である。
【図４】実施の形態に係る整相加算部１０４１の構成を示す機能ブロック図である。
【図５】実施の形態に係る主領域Ｂｘ、副領域Ｃｘを示す図である。
【図６】実施の形態に係る仮想受信開口Ｒｖ、受信開口Ｒｘと仮想送信開口Ｔｖ、送信開口Ｔｘとの関係を示す模式図である。
【図７】実施の形態に係る、送信開口Ｔｘから観測点Ｐｉｊを経由して受信振動子Ｒｋに到達する超音波の伝播経路を示す模式図である。
【図８】実施の形態に係る、仮想送信開口Ｔｖから観測点Ｑｍｎに到達する超音波の伝播経路を示す模式図である。
【図９】実施の形態に係る合成部１１４０の構成を示す機能ブロック図である。
【図１０】実施の形態に係る、各送信イベントにおける主領域Ｂｘ、副領域Ｃｘを示す図である。
【図１１】実施の形態に係る、各送信イベントにおける主領域Ｂｘ、副領域Ｃｘを示す図である。
【図１２】実施の形態に係る加算処理部１１４０１における合成音響線信号を合成する処理を示す模式図である。
【図１３】実施の形態に係る、合成音響線信号における最大重畳数と、増幅処理部１１４０２における増幅処理の概要を示す模式図である。
【図１４】実施の形態に係る受信ビームフォーマ部１０４のビームフォーミング処理動作を示すフローチャートである。
【図１５】実施の形態に係る受信ビームフォーマ部１０４における観測点Ｐｉｊ、Ｑｍｎについての音響線信号生成動作を示すフローチャートである。
【図１６】実施の形態に係る受信ビームフォーマ部１０４における観測点Ｐｉｊ、Ｑｍｎについての音響線信号生成動作を説明するための模式図である。
【図１７】変形例１に係る受信開口設定部により設定された仮想受信開口Ｒｖ、受信開口Ｒｘと仮想送信開口Ｔｖ、送信開口Ｔｘとの関係を示す模式図である。
【図１８】変形例１に係る受信ビームフォーマ部のビームフォーミング処理動作を示すフローチャートである。
【図１９】変形例１に係る受信ビームフォーマ部における観測点Ｐｉｊ、Ｑｍｎについての音響線信号生成動作を説明するための模式図である。
【図２０】変形例２に係る主領域Ｂｘ、副領域Ｃｘを示す図である。
【図２１】従来の合成開口法における画像化領域を示す図である。
【発明を実施するための形態】
【００１２】
≪発明を実施するための形態に至った経緯≫
発明者は、合成開口法を用いる超音波診断装置において、画像化領域を拡大させるために各種の検討を行った。
【００１３】
超音波診断装置においては、一般に、複数の振動子によって行われる被検体への超音波送信が行われる際、被検体のある深さで超音波ビームがフォーカスを結ぶよう、波面を集束させる送信ビームフォーミングがなされる（以下、フォーカスが結ばれる深さを「フォーカス深さ」とする）。そのため、１度の超音波の送信（送信イベント）によって、超音波送信に用いられる複数の振動子（以下、「送信振動子列」とする）から、超音波主照射領域に主として超音波が照射される。ここで、「超音波主照射領域」とは、送信振動子列を構成する各振動子から送波される超音波の伝搬する領域のことである。送信フォーカス点が１点である場合には、超音波主照射領域は、送信振動子列を底辺とし、底辺の両端のそれぞれから送信フォーカス点を通る２つの直線で囲まれる砂時計形状の領域となり、波面は、送信フォーカス点を中心とした円弧状となる。なお、必ずしも超音波ビームが１点でフォーカスを結ぶとは限らず、例えば、１．５素子分から数素子分程度にフォーカスした領域（以下、「フォーカス領域」とする）に集束するだけの場合もあるが、この場合、超音波主照射領域はフォーカス深さまでは列方向の幅が狭まり、フォーカス深さでフォーカス領域の列方向の幅となり、フォーカス深さより深い領域では再び列方向が広がる形状となる。すなわち、超音波主照射領域は、フォーカス深さで列方向の幅が最も狭くなり、それ以外の深さでは、フォーカス深さまでの距離に応じて列方向の幅が広がる形状となる。
【００１４】
合成開口法では、１回の送信イベントにおいて、超音波主照射領域のいずれの点に対しても観測点を設定できるため、超音波主照射領域全域を、音響線信号が生成される領域（以下、「対象領域」とする）とするのが好ましい。１つの送信イベントでは超音波画像を生成する領域全体を対象領域とすることができないために、１フレームの超音波画像を生成するためには、対象領域の異なる複数の送信イベントを行う。超音波の利用効率の観点からは対象領域は大きいことが好ましく、また、連続する２つの送信イベントの対象領域の重複面積が大きいほうが、空間解像度や信号Ｓ／Ｎ比の向上のために好ましい。
【００１５】
したがって、砂時計形状の超音波主照射領域の全域を対象領域として送信イベントを行い、送信イベントに同期して、超音波主照射領域と対象領域とを１素子分ずつ移動させて音響線信号を生成し（以下、１回の送信イベントで生成される音響線信号を「サブフレーム音響線信号」とする）、複数のサブフレーム音響線信号を合成して１フレームの超音波画像を生成することが従来行われている。
【００１６】
一方で、設定上の送信振動子列（以下、「仮想送信振動子列」とする）の一部が振動子列からはみ出る場合には、物理的に振動子が存在しない部分を送信振動子列とすることができない。したがって、素子列方向における仮想送信振動子列の中心（以下、「送信開口中心」とする）の設定可能範囲は、超音波プローブの振動子列の幅を超えることができない。さらに、対象領域の素子列方向の幅はフォーカス深さで最小となるため、画像化領域の素子列方向の幅もフォーカス深さにおいて最小となる。以上の理由により、例えば、リニアプローブの場合、超音波プローブの振動子列を底辺とする長方形状、扇形状、台形状の範囲が、画像化領域の最大範囲となる。また、例えば、コンベックスプローブの場合、略扇形上となる画像化領域において、超音波プローブの振動子列がなす円弧の中心角が、画像化領域の中心角の最大値となる。
【００１７】
これに対して画像化領域を素子列方向に拡大する方法としては、上述したように、１つには、送信イベント毎に超音波ビームの送信方向を変更し、超音波ビームを放射状に送信することにより、超音波ビームの通過領域を拡大する方法がある。これは、図２１（ｂ）に示すように、超音波プローブの振動子列の中心から送信開口中心が離れるほど、超音波ビームの進行方向を振動子列の中心から離れるようにステアリングを行う方法である。具体的には、送信開口中心がｘの正の方向（図中右側）にあるほど、ｘの正の方向（図中右上側）に超音波ビームの進行方向を傾け、同様に、送信開口中心がｘの負の方向（図中左側）にあるほど、ｘの負の方向（図中左上側）に超音波ビームの進行方向を傾けるというものである。これにより、画像化領域は、リニアプローブの場合、振動子列の幅を上底とする台形状となり、下底と上底の長さの差だけ、素子列方向に画像化領域を拡大することができる。しかしながら、この方法では、対象領域が放射状に広がるため、深さが大きくなるほど連続する２つの送信イベントの対象領域の重複面積が小さくなり、合成数が小さくなって空間解像度や信号Ｓ／Ｎ比が低下することとなる。
【００１８】
また、もう一つの画像化領域を拡大する方法としては、上述したように、通常の送信イベントに加えて、超音波ビームの送信方向を外側に向けた送信イベントを追加する方法がある。すなわち、図２１（ｃ）に示すように、通常の送信イベントで図２１（ａ）と同様の画像を生成し、通常の画像化領域の素子列方向における外側に追加領域を設けて、追加領域については図２１（ｂ）と同様に、超音波ビームの進行方向を追加領域に向けてステアリングを行う方法である。これにより、通常の画像化領域については空間解像度や信号Ｓ／Ｎ比に影響を与えることなく、画像化領域を素子列方向に拡大することができる。しかしながら、この方法では、追加領域を画像化するための送信イベントを別途行わなければならないため、１枚の超音波画像を生成するのに必要な送信イベント数が増加し、フレームレートの向上が困難となる。
【００１９】
そこで、発明者は、上記課題に鑑み、波面を集束させる送信ビームフォーミングと合成開口法の組み合わせにおいて、空間解像度や信号Ｓ／Ｎ比、フレームレートに悪影響を与えることなく画像化領域を素子列方向に拡大させる技術について検討を行い、実施の形態に係る超音波信号処理装置、超音波診断装置、および、超音波信号処理方法に想到するに至ったものである。
【００２０】
以下、実施の形態に係る超音波画像処理方法及びそれを用いた超音波診断装置について図面を用いて詳細に説明する。
【００２１】
≪実施の形態≫
＜全体構成＞
以下、実施の形態に係る超音波診断装置１００について、図面を参照しながら説明する。
【００２２】
図１は、実施の形態に係る超音波診断システム１０００の機能ブロック図である。図１に示すように、超音波診断システム１０００は、被検体に向けて超音波を送信しその反射波の受信する複数の振動子１０１ａを有するプローブ１０１、プローブ１０１に超音波の送受信を行わせプローブ１０１からの出力信号に基づき超音波画像を生成する超音波診断装置１００、超音波画像を画面上に表示する表示部１０６を有する。プローブ１０１、表示部１０６は、それぞれ、超音波診断装置１００に各々接続可能に構成されている。図１は超音波診断装置１００に、プローブ１０１、表示部１０６が接続された状態を示している。なお、プローブ１０１と、表示部１０６とは、超音波診断装置１００の内部にあってもよい。
【００２３】
＜超音波診断装置１００の構成＞
超音波診断装置１００は、プローブ１０１の複数ある振動子１０１ａのうち、送信又は受信の際に用いる振動子を各々に選択し、選択された振動子に対する入出力を確保するマルチプレクサ部１０２、超音波の送信を行うためにプローブ１０１の各振動子１０１ａに対する高電圧印加のタイミングを制御する送信ビームフォーマ部１０３と、プローブ１０１で受信した超音波の反射波に基づき、複数の振動子１０１ａで得られた電気信号を増幅し、Ａ／Ｄ変換し、受信ビームフォーミングして音響線信号を生成する受信ビームフォーマ部１０４を有する。また、受信ビームフォーマ部１０４からの出力信号に基づいて超音波画像（Ｂモード画像）を生成する超音波画像生成部１０５、受信ビームフォーマ部１０４が出力する音響線信号及び超音波画像生成部１０５が出力する超音波画像を保存するデータ格納部１０７と、各構成要素を制御する制御部１０８を備える。
【００２４】
このうち、マルチプレクサ部１０２、送信ビームフォーマ部１０３、受信ビームフォーマ部１０４、超音波画像生成部１０５は、超音波信号処理装置１５０を構成する。
【００２５】
超音波診断装置１００を構成する各要素、例えば、マルチプレクサ部１０２、送信ビームフォーマ部１０３、受信ビームフォーマ部１０４、超音波画像生成部１０５、制御部１０８は、それぞれ、例えば、ＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ  Ｓｐｅｃｉｆｉｃ  Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ）などのハードウエア回路により実現される。あるいは、プロセッサなどのプログラマブルデバイスとソフトウェアにより実現される構成であってもよい。プロセッサとしてはＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）やＧＰＧＰＵを用いることができ、ＧＰＵを用いる構成はＧＰＧＰＵ（Ｇｅｎｅｒａｌ−Ｐｕｒｐｏｓｅ  ｃｏｍｐｕｔｉｎｇ  ｏｎ  Ｇｒａｐｈｉｃｓ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）と呼ばれる。これらの構成要素は一個の回路部品とすることができるし、複数の回路部品の集合体にすることもできる。また、複数の構成要素を組合せて一個の回路部品とすることができるし、複数の回路部品の集合体にすることもできる。
【００２６】
データ格納部１０７は、コンピュータ読み取り可能な記録媒体であり、例えば、フレキシブルディスク、ハードディスク、ＭＯ、ＤＶＤ、ＤＶＤ−ＲＡＭ、ＢＤ、半導体メモリ等を用いることができる。また、データ格納部１０７は、超音波診断装置１００に外部から接続された記憶装置であってもよい。
【００２７】
なお、本実施の形態に係る超音波診断装置１００は、図１で示した構成の超音波診断装置に限定されない。例えば、マルチプレクサ部１０２がなく、送信ビームフォーマ部１０３と受信ビームフォーマ部１０４とが直接、プローブ１０１の各振動子１０１ａに接続されていてもよい。また、プローブ１０１に送信ビームフォーマ部１０３や受信ビームフォーマ部１０４、またその一部などが内蔵される構成であってもよい。これは、本実施の形態に係る超音波診断装置１００に限られず、後に説明する他の実施の形態や変形例に係る超音波診断装置でも同様である。
【００２８】
＜超音波診断装置１００の主要部の構成＞
実施の形態１に係る超音波診断装置１００は、プローブ１０１の各振動子１０１ａから超音波送信を行わせる送信ビームフォーマ部１０３と、プローブ１０１での超音波反射波の受信から得た電気信号を演算して超音波画像を生成するための音響線信号を生成する受信ビームフォーマ部１０４に特徴を有する。そのため、本明細書では、主に、送信ビームフォーマ部１０３及び受信ビームフォーマ部１０４について、その構成及び機能を説明する。なお、送信ビームフォーマ部１０３及び受信ビームフォーマ部１０４以外の構成については、公知の超音波診断装置に使われるものと同じ構成を適用可能であり、公知の超音波診断装置のビームフォーマ部に本実施の形態に係るビームフォーマ部を置き換えて使用することが可能である。
【００２９】
以下、送信ビームフォーマ部１０３と、受信ビームフォーマ部１０４の構成について説明する。
【００３０】
１．送信ビームフォーマ部１０３
送信ビームフォーマ部１０３は、マルチプレクサ部１０２を介してプローブ１０１と接続され、プローブ１０１から超音波の送信を行うためにプローブ１０１に存する複数の振動子１０１ａの全てもしくは一部に当たる送信振動子列からなる送信開口Ｔｘに含まれる複数の振動子の各々に対する高電圧印加のタイミングを制御する。送信ビームフォーマ部１０３は送信部１０３１から構成される。
【００３１】
送信部１０３１は、制御部１０８からの送信制御信号に基づき、プローブ１０１に存する複数の振動子１０１ａ中、送信開口Ｔｘに含まれる各振動子に超音波ビームを送信させるためのパルス状の送信信号を供給する送信処理を行う。具体的には、送信部１０３１は、例えば、クロック発生回路、パルス発生回路、遅延回路を備えている。クロック発生回路は、超音波ビームの送信タイミングを決定するクロック信号を発生させる回路である。パルス発生回路は、各振動子を駆動するパルス信号を発生させるための回路である。遅延回路は、超音波ビームの送信タイミングを振動子毎に遅延時間を設定し、遅延時間だけ超音波ビームの送信を遅延させて超音波ビームのフォーカシングを行うための回路である。
【００３２】
送信部１０３１は、超音波送信ごとに仮想送信開口Ｔｖを列方向に、所定の移動ピッチＭｐだけ移動させながら超音波送信を繰り返し、プローブ１０１に存する全ての振動子１０１ａから超音波送信を行う。ここで、仮想送信開口Ｔｖとは、振動子の有無に関係なく設定される設定上の送信開口であり、その列長（以下、「送信開口長」と呼ぶ）は、送信イベントに関係なく一定である。仮想送信開口Ｔｖの範囲内全域に振動子が存在する場合には、仮想送信開口Ｔｖと送信開口Ｔｘは一致する。一方、仮想送信開口Ｔｖの範囲内に振動子が存在しない場所が存在する場合、仮想送信開口Ｔｖのうち、振動子の存在する部分が送信開口Ｔｘとなる。本実施の形態では、移動ピッチＭｐを振動子１つ分とし、超音波送信ごとに、仮想送信開口Ｔｖが振動子１つ分ずつ移動していく。なお、移動ピッチＭｐは振動子１つ分に限られず、例えば、振動子０．５個分としてもよい。仮想送信開口Ｔｖの位置を示す情報、および送信開口Ｔｘに含まれる振動子の位置を示す情報は制御部１０８を介してデータ格納部１０７に出力される。例えば、プローブ１０１に存する振動子１０１ａ全数を１９２としたとき、送信開口Ｔｘを構成する振動子列の数として、例えば２０〜１００を選択してもよく、超音波送信毎に移動ピッチＭｐだけ移動させる構成としてもよい。以後、送信部１０３１により同一の送信開口Ｔｘから行われる超音波送信を「送信イベント」と称呼する。
【００３３】
図２は、送信ビームフォーマ部１０３による超音波送信波の伝播経路を示す模式図である。ある送信イベントにおいて、超音波送信に寄与するアレイ状に配列された振動子１０１ａの列（送信振動子列）を送信開口Ｔｘとして図示している。
【００３４】
送信ビームフォーマ部１０３において、仮想送信開口Ｔｖの中心に位置する振動子ほど送信タイミングを遅らせるように各振動子の送信タイミングを制御する。これにより、送信開口Ｔｘ内の振動子列から送信された超音波送信波は、被検体のある深度（Ｆｏｃａｌ  ｄｅｐｔｈ）において、波面がある一点、すなわち送信フォーカス点Ｆ（Ｆｏｃａｌ  ｐｏｉｎｔ）で、フォーカスがあう（集束する）状態となる。送信フォーカス点Ｆの深さ（Ｆｏｃａｌ  ｄｅｐｔｈ）（以下、「フォーカス深さ」とする）は、任意に設定することができる。送信フォーカス点Ｆで合焦した波面は、再び拡散し、送信開口Ｔｘを底とし送信フォーカス点Ｆを節とする交差する２つの直線で区切られた砂時計型の空間内を超音波送信波が伝播する。すなわち、送信開口Ｔｘで放射された超音波は、次第にその空間上での幅（図中の横軸方向）を小さくし、送信フォーカス点Ｆでその幅を最小化し、それよりも深部（図中では上部）に進行するにしたがって、再び、その幅を大きくしながら拡散し、伝播することとなる。この砂時計型の領域が超音波主照射領域Ａｘである。なお、上述したように、超音波主照射領域Ａｘは、１点の送信フォーカス点Ｆに集束せず、フォーカス領域に集束してもよい。
【００３５】
なお、上述したように、超音波主照射領域Ａｘは、送信振動子列の各振動子から送波された超音波の位相が揃っている領域であり、超音波主照射領域Ａｘの外部にも超音波送信波は伝播している。もっとも、超音波主照射領域Ａｘの外部においては、送信振動子列の各振動子から送波された超音波の位相が揃っていないため、超音波主照射領域Ａｘ内部と比べて超音波送信波が劣化しており、特に、超音波主照射領域Ａｘから遠ざかるほど劣化が顕著となる。
【００３６】
２．受信ビームフォーマ部１０４の構成
受信ビームフォーマ部１０４は、プローブ１０１で受信した超音波の反射波に基づき、複数の振動子１０１ａで得られた電気信号から音響線信号を生成する。なお、「音響線信号」とは、ある観測点に対する、整相加算処理がされた後の信号である。整相加算処理については後述する。図３は、受信ビームフォーマ部１０４の構成を示す機能ブロック図である。図３に示すように、受信ビームフォーマ部１０４は、受信部１０４０、整相加算部１０４１、合成部１１４０を備える。
【００３７】
以下、受信ビームフォーマ部１０４を構成する各部の構成について説明する。
【００３８】
（１）受信部１０４０
受信部１０４０は、マルチプレクサ部１０２を介してプローブ１０１と接続され、送信イベントに同期してプローブ１０１での超音波反射波の受信から得た電気信号を増幅した後ＡＤ変換した受信信号（ＲＦ信号）を生成する回路である。送信イベントの順に時系列に受信信号を生成しデータ格納部１０７に出力し、データ格納部１０７に受信信号を保存する。
【００３９】
ここで、受信信号（ＲＦ信号）とは、各振動子にて受信された反射超音波から変換された電気信号をＡ／Ｄ変換したデジタル信号であり、各振動子にて受信された超音波の送信方向（被検体の深さ方向）に連なった信号の列を形成している。
【００４０】
送信イベントでは、上述のとおり、送信部１０３１は、プローブ１０１に存する複数の振動子１０１ａ中、送信開口Ｔｘに含まれる複数の振動子の各々に超音波ビームを送信させる。これに対し、受信部１０４０は、送信イベントに同期してプローブ１０１に存する複数の振動子１０１ａの一部又は全部にあたる振動子の各々が得た反射超音波に基づいて、各振動子に対する受信信号の列を生成する。ここで、反射超音波を受波する振動子を「受波振動子」と称呼する。受波振動子の数は、送信開口Ｔｘに含まれる振動子の数よりも多いことが好ましい。また、受波振動子の数はプローブ１０１に存する振動子１０１ａの全数としてもよい。
【００４１】
送信部１０３１は、送信イベントに同期して送信開口Ｔｘを列方向に移動ピッチＭｐだけ移動させながら超音波送信を繰り返し、プローブ１０１に存する複数の振動子１０１ａ全体から超音波送信を行う。受信部１０４０は、送信イベントに同期して各受波振動子に対する受信信号の列を生成し、生成された受信信号はデータ格納部１０７に保存される。
【００４２】
（２）整相加算部１０４１
整相加算部１０４１は、送信イベントに同期して、被検体内においてサブフレーム音響線信号の生成を行う対象領域である主領域Ｂｘ、副領域Ｃｘを設定する。次に、主領域Ｂｘ内に存在する複数の観測点Ｐｉｊと副領域Ｃｘ内に存在する複数の観測点Ｑｍｎの各々について、観測点から各受信振動子Ｒｋが受信した受信信号列を整相加算する。そして、各観測点における音響線信号の列を算出することによりサブフレーム音響線信号を生成する回路である。図４は、整相加算部１０４１の構成を示す機能ブロック図である。図４に示すように、整相加算部１０４１は、対象領域設定部１０４２、受信開口設定部１０４３、送信時間算出部１０４４、受信時間算出部１０４５、遅延量算出部１０４６、遅延処理部１０４７、重み算出部１０４８、及び加算部１０４９を備える。
【００４３】
以下、整相加算部１０４１を構成する各部の構成について説明する。
【００４４】
ｉ）対象領域設定部１０４２
対象領域設定部１０４２は、被検体内においてサブフレーム音響線信号の生成を行う対象領域である主領域Ｂｘ、副領域Ｃｘを設定する。「対象領域」とは、送信イベントに同期して被検体内においてサブフレーム音響線信号の生成が行われるべき信号上の領域であり、主領域Ｂｘ内の観測点Ｐｉｊおよび副領域Ｃｘ内の観測点Ｑｍｎについて音響線信号が生成される。主領域Ｂｘ、副領域Ｃｘは、音響線信号の生成が行われる観測対象点の集合として、１回の送信イベントに同期して計算の便宜上設定される。
【００４５】
ここで、「サブフレーム音響線信号」とは、１回の送信イベントから生成される対象領域内に存在する全ての観測点に対する音響線信号の集合である。サブフレーム音響線信号は、主領域Ｂｘ内に存在する全ての観測点Ｐｉｊに対する音響線信号の集合である「主領域音響線信号」と、副領域Ｃｘ内に存在する全ての観測点Ｑｍｎに対する音響線信号の集合である「副領域音響線信号」とからなる。なお、「サブフレーム」とは、１回の送信イベントで得られ、対象領域内に存在する全ての観測点に対応するまとまった信号を形成する単位をさす。取得時間の異なる複数のサブフレームを合成したものがフレームとなる。
【００４６】
対象領域設定部１０４２は、まず、画像化主領域Ｈｘと追加領域Ｈａ、Ｈｂからなる画像化領域を設定する。そして、送信イベントに同期して、送信ビームフォーマ部１０３から取得する送信開口Ｔｘおよび仮想送信開口Ｔｖの位置を示す情報に基づき対象領域を設定する。
【００４７】
図５は、画像化領域および対象領域を示す模式図である。図５（ａ）に示すように、画像化主領域Ｈｘは、振動子列１０１ａを底辺とする長方形状の領域として設定される。そして、追加領域Ｈａ、追加領域Ｈｂはそれぞれ、素子列方向において画像化主領域Ｈｘに隣接するように設定される。本実施の形態では、追加領域Ｈａ、追加領域Ｈｂの深さｙが０（被検体とプローブの界面）における列方向の幅は送信開口幅をＴｘとしたとき、その送信開口幅Ｔｘの１／２であり、深さｙが追加領域Ｈａ、Ｈｂの最深部における列方向の幅はそれぞれＷａ、Ｗｂであり、深さｙが大きくなるに従って素子列方向の幅が増加する。本実施の形態では、Ｗａ＝Ｗｂ＞Ｔｘ／２であるとする。
【００４８】
次に、対象領域について説明する。送信開口Ｔｘが画像化主領域Ｈｘの端から十分離れている場合には、図５（ａ）に示すように、超音波主照射領域Ａｘと画像化領域Ｈｘの重複する領域の全域を、主領域Ｂｘとして設定する。このとき、副領域Ｃｘは設定されず、空領域となる。一方、送信開口Ｔｘが画像化主領域Ｈｘの端に近い場合には、主領域Ｂｘと副領域Ｃｘとを設定する。具体的には、図５（ｂ）に示すように、仮想送信開口Ｔｖを送信開口とみなして、超音波主照射領域Ａｘと画像化領域Ｈｘまたは追加領域Ｈａ、Ｈｂとの重複する領域の全域を主領域Ｂｘとして設定する。すなわち、仮想送信開口Ｔｖを底辺とし、底辺の両端のそれぞれから送信フォーカス点を通る２つの直線で囲まれる砂時計形状の領域を主領域Ｂｘとする。さらに、以下の方法により副領域Ｃｘを決定する。まず、仮想送信開口Ｔｖの中心軸である送信開口中心軸Ｔｘｏと追加領域Ｈａの右端が一致するように追加領域Ｈａを素子列方向に沿って右に平行移動した領域から主領域Ｂｘを除く領域を第１領域Ｊａとする。すなわち、送信開口中心から左側にＴｘ／２だけ離れた点と、送信開口中心軸Ｔｘｏ上にあり深さが画像化主領域Ｈｘの最深部にある点Ｃから左側にＷａだけ離れた点Ｊａとを結ぶ直線と、送信開口中心軸Ｔｘｏとに挟まれた領域のうち、主領域Ｂｘに重複しない領域を第１領域Ｊａとする。同様に、送信開口中心軸Ｔｘｏと追加領域Ｈｂの左端が一致するように追加領域Ｈｂを素子列方向に沿って左に平行移動した領域から主領域Ｂｘを除く領域を第２領域Ｊｂとする。すなわち、送信開口中心から右側にＴｘ／２だけ離れた点と、点Ｃから右側にＷｂだけ離れた点Ｊｂとを結ぶ直線と送信開口中心軸Ｔｘｏとに挟まれた領域のうち、主領域Ｂｘに重複しない領域を第２領域Ｊｂとする。そして、第１領域Ｊａと追加領域Ｈａとの重複領域、および、第２領域Ｊｂと追加領域Ｈｂとの重複領域を、副領域Ｃｘとする。なお、第１領域Ｊａと追加領域Ｈａとが重複しない場合、または、第２領域Ｊｂと追加領域Ｈｂとが重複しない場合は、副領域Ｃｘは設定されない。
【００４９】
図１０、図１１は、各送信イベントにおける主領域Ｂｘと副領域Ｃｘとをより詳細に示す模式図である。仮想送信開口Ｔｖが左端、すなわち送信開口中心が振動子列１０１ａの左端の振動子である場合には、図１０（ａ）に示すように、仮想送信開口Ｔｖを底辺とし、底辺の両端のそれぞれから送信フォーカス点を通る２つの直線で囲まれる砂時計形状の領域を主領域Ｂｘとする。このとき、送信開口中心軸Ｔｘｏと追加領域Ｈａの右端は一致しているので、追加領域Ｈａから主領域Ｂｘを除いた領域が第１領域Ｊａとなる。したがって、第１領域Ｊａと追加領域Ｈａの重複領域である副領域Ｃｘは、必然的に第１領域Ｊａの全域、すなわち、追加領域Ｈａから主領域Ｂｘを除いた領域である。これに対し、仮想送信開口Ｔｖが右に移動したが左端が振動子列１０１ａの外にある場合には、図１０（ｂ）に示すように、仮想送信開口Ｔｖを底辺とし、底辺の両端のそれぞれから送信フォーカス点を通る２つの直線で囲まれる砂時計形状の領域を主領域Ｂｘとする。このとき、送信開口中心軸Ｔｘｏは追加領域Ｈａの右端より右側にあるので、追加領域Ｈａを右にシフトした領域から主領域Ｂｘを除いた領域が第１領域Ｊａとなる。したがって、第１領域Ｊａと追加領域Ｈａの重複領域である副領域Ｃｘは、図１０（ａ）で示す場合と比べ、右にシフトし、かつ、画像化主領域Ｈｘとの重複部分がなくなった形状となる。一方、送信開口中心が振動子列１０１ａの略中心にある場合は、図１０（ｃ）に示すように、送信開口Ｔｘが仮想送信開口Ｔｖと一致し、送信開口Ｔｘを底辺とし、底辺の両端のそれぞれから送信フォーカス点を通る２つの直線で囲まれる砂時計形状の領域を主領域Ｂｘとする。このとき、第１領域Ｊａ、第２領域Ｊｂはともに全域が画像化主領域Ｈｘ内に存在することとなるため、副領域Ｃｘは存在しない。また、仮想送信開口Ｔｖが右に移動し右端が振動子列１０１ａの外にある場合には、図１１（ａ）に示すように、仮想送信開口Ｔｖを底辺とし、底辺の両端のそれぞれから送信フォーカス点を通る２つの直線で囲まれる砂時計形状の領域を主領域Ｂｘとする。このとき、送信開口中心軸Ｔｘｏは追加領域Ｈｂの左端より左側にあるので、追加領域Ｈｂを左にシフトした領域から主領域Ｂｘを除いた領域が第１領域Ｊａとなる。したがって、第２領域Ｊｂと追加領域Ｈｂの重複領域である副領域Ｃｘは、図１０（ｂ）で示す場合と比べて、左右反転したような形状となる。さらに、仮想送信開口Ｔｖが右端、すなわち送信開口中心が振動子列１０１ａの右端の振動子である場合には、図１１（ｂ）に示すように、仮想送信開口Ｔｖを底辺とし、底辺の両端のそれぞれから送信フォーカス点を通る２つの直線で囲まれる砂時計形状の領域を主領域Ｂｘとする。このとき、送信開口中心軸Ｔｘｏと追加領域Ｈｂの左端は一致しているので、追加領域Ｈｂから主領域Ｂｘを除いた領域が第２領域Ｊｂとなる。したがって、第２領域Ｊｂと追加領域Ｈｂの重複領域である副領域Ｃｘは、必然的に第２領域Ｊｂの全域、すなわち、追加領域Ｈｂから主領域Ｂｘを除いた領域である。
【００５０】
なお、第１領域Ｊａと第２領域Ｊｂの範囲は、追加領域Ｈａ、Ｈｂそれぞれの形状と、送信開口中心の位置により一意に定まるので、画像化主領域Ｈｘと追加領域Ｈａ、Ｈｂそれぞれの位置及び大きさが既知である場合、送信開口中心の位置に基づいて副領域Ｃｘが存在するか否かを判定することができる。例えば、上述のように追加領域Ｈａ、Ｈｂを定める場合、点Ｃと画像化主領域Ｈｘの左端との距離がＷａ未満の場合、または、点Ｃと画像化主領域Ｈｘの右端との距離がＷｂ未満の場合、副領域Ｃｘが存在し、それ以外の場合には副領域Ｃｘは存在しない。したがって、対象領域設定部１０４２は、画像化主領域Ｈｘと追加領域Ｈａ、Ｈｂそれぞれの位置及び大きさに対して、副領域Ｃｘが存在する送信開口中心の位置の範囲を保持していてもよい。このようにすることで、副領域Ｃｘが空領域であるか否かを第１領域Ｊａ、第２領域Ｊｂの位置を算出しなくても判定することができるため、副領域Ｃｘが空領域である場合には主領域Ｂｘの設定のみを行えばよく、演算量の低減に奏功する。なお、対象領域設定部１０４２は、副領域Ｃｘが存在する送信開口中心の位置の範囲に替えて、副領域Ｃｘが存在する仮想送信開口Ｔｖの位置の範囲、または、副領域Ｃｘが存在する送信フォーカス点Ｆの位置の範囲を保持してもよい。
【００５１】
設定された対象領域は送信時間算出部１０４４、受信時間算出部１０４５、遅延処理部１０４７に出力される。
【００５２】
ｉｉ）受信開口設定部１０４３
受信開口設定部１０４３は、制御部１０８からの制御信号と、送信ビームフォーマ部１０３からの送信開口Ｔｘの位置を示す情報とに基づき、プローブ１０１に存する複数の振動子の一部に当たり、列中心が観測点に最も空間的に近接する振動子と合致する振動子列（受信振動子列）を受信振動子として選択して受信開口Ｒｘを設定する回路である。
【００５３】
受信開口設定部１０４３は、列中心が観測点ＰｉｊまたはＱｍｎに最も空間的に近接する振動子Ｘｋと合致するよう仮想受信開口Ｒｖを選択する。仮想送信開口Ｒｖに含まれる各振動子は、受信開口Ｒｘを構成する。送信開口と同様、仮想受信開口Ｒｖの範囲内全域に振動子１０１ａが存在する場合には、受信開口Ｒｘは仮想受信開口Ｒｖと一致する。一方、仮想受信開口Ｒｖの範囲内に振動子１０１ａが存在しない部分が存在する場合には、仮想受信開口Ｒｖのうち振動子１０１ａが存在している範囲が受信開口Ｒｘとなる。図６は、受信開口設定部１０４３により設定された仮想受信開口Ｒｖ、受信開口Ｒｘと仮想送信開口Ｔｖ、送信開口Ｔｘとの関係を示す模式図である。図６に示すように、仮想受信開口Ｒｖの中心が、観測点Ｐｉｊに最も空間的に近接する振動子Ｘｋと合致するように仮想受信開口Ｒｖが選択される。そのため、仮想受信開口Ｒｖおよび受信開口Ｒｘの位置は、観測点Ｐｉｊ、Ｑｍｎの位置によって定まり、送信イベントに同期して変動する送信開口Ｔｘの位置に基づいては変化しない。すなわち、異なる送信イベントであっても、同一位置にある観測点Ｐｉｊ、Ｑｍｎについての音響線信号を生成する処理においては、同一の受信開口Ｒｘ内の受信振動子Ｒｋによって取得された受信信号に基づき整相加算が行われる。
【００５４】
また、超音波主照射領域全体からの反射波を受信するために、受信開口Ｒｘに含まれる振動子の数は、対応する送信イベントにおける送信開口Ｔｘに含まれる振動子の数以上に設定することが好ましい。受信開口Ｒｘを構成する振動子列の数は、例えば３２、６４、９６、１２８、１９２等としてもよい。
【００５５】
仮想受信開口Ｒｖおよび受信開口Ｒｘの設定は、少なくとも列方向における観測点ＰｉｊおよびＱｍｎの最大数と同じ回数だけ行われる。また、仮想受信開口Ｒｖおよび受信開口Ｒｘの設定は、送信イベントに同期して漸次行われる構成であってもよく、あるいは、全ての送信イベントが終了した後に、各送信イベントに対応した仮想受信開口Ｒｖおよび受信開口Ｒｘの設定が送信イベントの回数分まとめて行われる構成であってもよい。
【００５６】
選択された受信開口Ｒｘの位置を示す情報は制御部１０８を介してデータ格納部１０７に出力される。
【００５７】
データ格納部１０７は、受信開口Ｒｘの位置を示す情報と受信振動子に対応する受信信号とを、送信時間算出部１０４４、受信時間算出部１０４５、遅延処理部１０４７、重み算出部１０４８に出力する。
【００５８】
ｉｉｉ）送信時間算出部１０４４
送信時間算出部１０４４は、送信された超音波が被検体中の観測点Ｐに到達する送信時間を算出する回路である。送信イベントに対応して、データ格納部１０７から取得した送信開口Ｔｘに含まれる振動子の位置を示す情報と、対象領域設定部１０４２から取得した対象領域の位置を示す情報とに基づき、対象領域内に存在する任意の観測点ＰｉｊおよびＱｍｎについて、送信された超音波が被検体中の観測点Ｐｉｊ、Ｑｍｎに到達する送信時間を算出する。
【００５９】
図７は、送信開口Ｔｘから放射され主領域Ｂｘ内の任意の位置にある観測点Ｐｉｊにおいて反射され受信開口Ｒｘ内に位置する受信振動子Ｒｋに到達する超音波の伝播経路を説明するための模式図である。なお、図７（ａ）は観測点Ｐｉｊがフォーカス深さより深い場合、図７（ｂ）は観測点Ｐｉｊの深さがフォーカス深さ以下である場合を示している。
【００６０】
送信開口Ｔｘから放射された送信波は、経路４０１を通って送信フォーカス点Ｆにて波面が集束し、再び、拡散する。送信波が集束または拡散する途中で観測点Ｐｉｊに到達し、観測点Ｐｉｊで音響インピーダンスに変化があれば反射波を生成し、その反射波がプローブ１０１における受信開口Ｒｘ内の受信振動子Ｒｋに戻っていく。送信フォーカス点Ｆは送信ビームフォーマ部１０３の設計値として規定されているので、送信フォーカス点Ｆと任意の観測点Ｐｉｊとの間の経路４０２の長さは幾何学的に算出することができる。
【００６１】
送信時間の算出方法を、以下、さらに詳細に説明する。
【００６２】
まず、観測点Ｐｉｊがフォーカス深さより深い場合について、図７（ａ）を用いて説明する。観測点Ｐｉｊがフォーカス深さより深い場合は、送信開口Ｔｘから放射された送信波が、経路４０１を通って送信フォーカス点Ｆに到達し、送信フォーカス点Ｆから経路４０２を通って観測点Ｐｉｊに到達したものとして算出する。したがって、送信波が経路４０１を通過する時間と、経路４０２を通過する時間を合算した値が、送信時間となる。具体的な算出方法としては、例えば、経路４０１の長さと経路４０２の長さとを加算した全経路長を、被検体内における超音波の伝搬速度で除算することで求められる。
【００６３】
一方、観測点Ｐｉｊが観測点Ｐｉｊの深さがフォーカス深さ以下である場合について、図７（ｂ）を用いて説明する。観測点Ｐｉｊの深さがフォーカス深さ以下である場合は、送信開口Ｔｘから放射された送信波が、経路４０１を通って送信フォーカス点Ｆに到達する時刻と、経路４０４を通って観測点Ｐｉｊに到達した後、観測点Ｐｉｊから経路４０２を通って送信フォーカス点Ｆに到達する時刻とが同一であるものとして算出する。つまり、送信波が経路４０１を通過する時間から、経路４０２を通過する時間を差し引いた値が、送信時間となる。具体的な算出方法としては、例えば、経路４０１の長さから経路４０２の長さを減算した経路長差を、被検体内における超音波の伝搬速度で除算することで求められる。
【００６４】
図８（ａ）は、送信開口Ｔｘから放射され副領域Ｃｘ内の任意の位置にある観測点Ｑｍｎにおいて反射され受信開口Ｒｘ内に位置する受信振動子Ｒｋに到達する超音波の伝播経路を説明するための模式図である。なお、図８（ａ）は観測点Ｑｍｎがフォーカス深さより深い場合について記載している。送信された超音波が被検体中の観測点Ｑｍｎに到達する送信時間としては、以下の３つの算出方法がある。１つは、上述したように、送信開口Ｔｘから放射された送信波が、経路４０１を通って送信フォーカス点Ｆに到達し、送信フォーカス点Ｆから経路４０２を通って観測点Ｑｍｎに到達したものとして算出する方法である。すなわち、送信波が経路４０１を通過する時間と、経路４０２を通過する時間を合算した値が、送信時間となる。具体的な算出方法としては、例えば、経路４０１の長さと経路４０２の長さとを加算した全経路長を、被検体内における超音波の伝搬速度で除算することで求められる。この送信時間を、以下、送信時間Ｔ2とする。もう１つは、送信開口Ｔｘから放射された送信波が、送信フォーカス点Ｆを経由せず、直接、経路４１１を通って観測点Ｑｍｎに到達したものとして算出する方法である。すなわち、送信波が経路４１１を通過する時間が、送信時間となる。具体的な算出方法としては、例えば、経路４１１の長さを、被検体内における超音波の伝搬速度で除算することで求められる。この送信時間を、以下、送信時間Ｔ1とする。最後の１つは、送信開口Ｔｘから放射された送信波が観測点Ｑｍｎと同じ深さにある参照点Ｒに到達した時刻と同時刻に、観測点Ｑｍｎにも到達したとみなして算出する方法である。すなわち、送信波が経路４１２を通過する時間が、送信時間となる。具体的な算出方法としては、例えば、経路４１２の長さを、被検体内における超音波の伝搬速度で除算することで求められる。この送信時間を、以下、送信時間Ｔ3とする。
【００６５】
上述した３つの送信時間Ｔ1、Ｔ2、Ｔ3のそれぞれについて、観測点Ｑｍｎの深さを横軸に、送信時間を縦軸に示した図が図８（ｂ）である。図８（ｂ）において、系統５０１は送信時間Ｔ2を示し、系統５０２は送信時間Ｔ1を示し、系統５０３は送信時間Ｔ3を示す。送信時間Ｔ2はフォーカス深さにおいて不連続となる。これは、経路４０２の長さの最小値がゼロではなく観測点Ｑｍｎと送信フォーカス点Ｆとの素子列方向における距離となるため、当該距離をｄ、被検体内の超音波速度をｃｓとすると、フォーカス深さを挟んで隣接する２つの観測点Ｑの送信時間差が２ｄ／ｃｓとなるためである。したがって、フォーカス深さ近傍については、送信時間Ｔ1またはＴ3を採用することが好ましい。また、観測点Ｑｍｎの深さが浅い場合には、送信開口Ｔｘ内の各振動子から直接観測点Ｑｍｎに超音波が伝搬していると考えられるため、送信時間Ｔ1を採用することが好ましい。一方、主領域Ｂｘと副領域Ｃｘとの界面付近では、主領域Ｂｘと副領域Ｃｘとで送信時間の算出方法を変えると不整合が目立つため、送信時間Ｔ2を採用することが好ましい。以上を総合すると、（１）フォーカス深さ近傍では、送信時間Ｔ1またはＴ3、（２）フォーカス深さより浅い領域では送信時間Ｔ1、（３）領域Ｂｘと副領域Ｃｘとの界面付近では、送信時間Ｔ2が、それぞれ好ましいこととなる。したがって、送信時間ＴMの算出方法として、例えば、以下のように算出することができる。
ＴM＝αＴ1＋（１−α）βＴ2＋（１−α）（１−β）Ｔ3
ここで、α、βは、それぞれ、０＜α≦１、０≦β≦１であり、α、βともに観測点Ｑｍｎの深さに対する関数であることが好ましい。具体的には、αはフォーカス深さ近傍において１であることが好ましく、主領域Ｂｘとの界面近傍においては小さいことが好ましい。また、βは、観測点Ｑｍｎの深さが大きいほど増加する値であることが好ましい。例えば、図８（ｂ）における系統５１１となるように、α、βを定めることができる。なお、送信時間ＴMは、上述した条件を満たし、且つ、深さの増加に対して単調増加するものであれば、ここで例示した場合に限られない。
【００６６】
送信時間算出部１０４４は、１回の送信イベントに対し、対象領域内の全ての観測点Ｐｉｊ、Ｑｍｎについて、送信された超音波が被検体中の観測点Ｐｉｊ、Ｑｍｎに到達する送信時間を算出して遅延量算出部１０４６に出力する。
【００６７】
ｉｖ）受信時間算出部１０４５
受信時間算出部１０４５は、観測点Ｐ、Ｑからの反射波が、受信開口Ｒｘに含まれる受信振動子Ｒｋの各々に到達する受信時間を算出する回路である。送信イベントに対応して、データ格納部１０７から取得した受信振動子Ｒｋの位置を示す情報と、対象領域設定部１０４２から取得した対象領域の位置を示す情報とに基づき対象領域内に存在する任意の観測点Ｐｉｊ、Ｑｍｎについて、送信された超音波が被検体中の観測点Ｐｉｊ、Ｑｍｎで反射され受信開口Ｒｘの各受信振動子Ｒｋに到達する受信時間を算出する。
【００６８】
上述のとおり、観測点Ｐｉｊ、Ｑｍｎに到達した送信波は、観測点Ｐｉｊ、Ｑｍｎで音響インピーダンスに変化があれば反射波を生成し、その反射波がプローブ１０１における受信開口Ｒｘ内の各受信振動子Ｒｋに戻っていく。受信開口Ｒｘ内の各受信振動子Ｒｋの位置情報はデータ格納部１０７から取得されるので、任意の観測点Ｐｉｊ、Ｑｍｎから各受信振動子Ｒｋまでの経路４０３の長さは幾何学的に算出することができる。
【００６９】
受信時間算出部１０４５は、１回の送信イベントに対し、対象領域内に存在する全ての観測点Ｐｉｊ、Ｑｍｎについて、送信された超音波が観測点Ｐｉｊ、Ｑｍｎで反射して各受信振動子Ｒｋに到達する受信時間を算出して遅延量算出部１０４６に出力する。
【００７０】
ｖ）遅延量算出部１０４６
遅延量算出部１０４６は、送信時間と受信時間とから受信開口Ｒｘ内の各受信振動子Ｒｋへの総伝播時間を算出し、当該総伝播時間に基づいて、各受信振動子Ｒｋに対する受信信号の列に適用する遅延量を算出する回路である。遅延量算出部１０４６は、送信時間算出部１０４４から送信された超音波が観測点Ｐｉｊ、Ｑｍｎに到達する送信時間と、観測点Ｐｉｊ、Ｑｍｎで反射して各受信振動子Ｒｋに到達する受信時間を取得する。そして、送信された超音波が各受信振動子Ｒｋへ到達するまでの総伝播時間を算出し、各受信振動子Ｒｋに対する総伝播時間の差異により、各受信振動子Ｒｋに対する遅延量を算出する。遅延量算出部１０４６は、対象領域内に存在する全ての観測点Ｐｉｊ、Ｑｍｎについて、各受信振動子Ｒｋに対する受信信号の列に適用する遅延量を算出して遅延処理部１０４７に出力する。
【００７１】
ｖｉ）遅延処理部１０４７
遅延処理部１０４７は、受信開口Ｒｘ内の受信振動子Ｒｋに対する受信信号の列から、各受信振動子Ｒｋに対する遅延量に相当する受信信号を、観測点Ｐｉｊ、Ｑｍｎからの反射超音波に基づく各受信振動子Ｒｋに対応する受信信号として同定する回路である。
【００７２】
遅延処理部１０４７は、送信イベントに対応して、受信開口設定部１０４３から受信振動子Ｒｋの位置を示す情報、データ格納部１０７から受信振動子Ｒｋに対応する受信信号、対象領域設定部１０４２から取得した対象領域の位置を示す情報、遅延量算出部１０４６から各受信振動子Ｒｋに対する受信信号の列に適用する遅延量を入力として取得する。そして、各受信振動子Ｒｋに対応する受信信号の列から、各受信振動子Ｒｋに対する遅延量を差引いた時間に対応する受信信号を観測点Ｐｉｊ、Ｑｍｎからの反射波に基づく受信信号として同定し、加算部１０４９に出力する。
【００７３】
ｖｉｉ）重み算出部１０４８
重み算出部１０４８は、受信開口Ｒｘの列方向の中心に位置する振動子に対する重みが最大となるよう各受信振動子Ｒｋに対する重み数列（受信アポダイゼーション）を算出する回路である。
【００７４】
図６に示すように、重み数列は受信開口Ｒｘ内の各振動子に対応する受信信号に適用される重み係数の数列である。重み数列は、送信フォーカス点Ｆを中心として対称な分布をなす。重み数列の分布の形状は、ハミング窓、ハニング窓、矩形窓などを用いることができ、分布の形状は特に限定されない。重み数列は、受信開口Ｒｘの列方向の中心に位置する振動子に対する重みが最大となるように設定され、重みの分布の中心軸は、受信開口中心軸Ｒｘｏと一致する。重み算出部１０４８は、受信開口設定部１０４３から出力される受信振動子Ｒｋの位置を示す情報を入力として、各受信振動子Ｒｋに対する重み数列を算出し加算部１０４９に出力する。
【００７５】
ｖｉｉｉ）加算部１０４９
加算部１０４９は、遅延処理部１０４７から出力される各受信振動子Ｒｋに対応して同定された受信信号を入力として、それらを加算して、観測点Ｐｉｊ、Ｑｍｎに対する整相加算された音響線信号を生成する回路である。あるいは、さらに、重み算出部１０４８から出力される各受信振動子Ｒｋに対する重み数列を入力として、各受信振動子Ｒｋに対応して同定された受信信号に、各受信振動子Ｒｋに対する重みを乗じて加算して、観測点Ｐｉｊ、Ｑｍｎに対する音響線信号を生成する構成としてもよい。遅延処理部１０４７において受信開口Ｒｘ内に位置する各受信振動子Ｒｋが検出した受信信号の位相を整えて加算部１０４９にて加算処理をすることにより、観測点Ｐｉｊ、Ｑｍｎからの反射波に基づいて各受信振動子Ｒｋで受信した受信信号を重ね合わせてその信号Ｓ／Ｎ比を増加し、観測点Ｐｉｊ、Ｑｍｎからの受信信号を抽出することができる。
【００７６】
１回の送信イベントとそれに伴う処理から、対象領域内の全ての観測点Ｐｉｊ、Ｑｍｎについて音響線信号を生成することができる。そして、図１０、図１１に示すように、送信イベントに同期して送信開口Ｔｘを列方向に移動ピッチＭｐだけ移動させながら超音波送信を繰り返し、プローブ１０１に存する全ての振動子１０１ａから超音波送信を行うことにより１フレームの合成された音響線信号であるフレーム音響線信号を生成する。
【００７７】
また、フレーム音響線信号を構成する観測点ごとの合成された音響線信号を、以後、「合成音響線信号」と称呼する。
【００７８】
加算部１０４９により、送信イベントに同期して対象領域内に存在する全ての観測点Ｐｉｊ、Ｑｍｎに対するサブフレーム音響線信号が生成される。以下、観測点Ｐｉｊに対するサブフレーム音響線信号を「主領域音響線信号」、観測点Ｑｍｎに対するサブフレーム音響線信号を「副領域音響線信号」と呼ぶ。生成されたサブフレーム音響線信号は、データ格納部１０７に出力され保存される。
【００７９】
（３）合成部１１４０
合成部１１４０は、送信イベントに同期して生成されるサブフレーム音響線信号からフレーム音響線信号を合成する回路である。図９は、合成部１１４０の構成を示す機能ブロック図である。図９に示すように、合成部１１４０は、加算処理部１１４０１−１、１１４０１−２と、増幅処理部１１４０２−１、１１４０２−２、結合部１１４０３を備える。加算処理部１１４０１−１、１１４０１−２は、それぞれ、フレーム音響線信号を合成するための一連のサブフレーム音響線信号の生成が終了したのち、データ格納部１０７に保持されている複数のサブフレーム音響線信号を読み出す。そして、各サブフレーム音響線信号に含まれる音響線信号が取得された観測点Ｐｉｊ、Ｑｍｎの位置を指標として複数のサブフレーム音響線信号を加算する。本実施の形態では、加算処理部１１４０１−１は複数の主領域音響線信号を加算対象として主合成音響線信号を合成し、また、加算処理部１１４０１−２は複数の副領域音響線信号を加算対象として副合成音響線信号を合成する。つまり、合成音響線信号は、主合成音響線信号と副合成音響線信号とで構成される。なお、観測点Ｐｉｊ、Ｑｍｎに対して主領域音響線信号と副領域音響線信号とが存在した場合においても、主領域音響線信号と副領域音響線信号との間の加算は行わない。
【００８０】
以下、合成部１１４０を構成する各部の構成について説明する。
【００８１】
ｉ）加算処理部１１４０１−１、１１４０１−２
加算処理部１１４０１−１は、フレーム音響線信号を合成するための一連のサブフレーム音響線信号の生成が終了したのち、データ格納部１０７に保持されている複数のサブフレーム音響線信号を読み出す。そして、各主領域音響線信号に含まれる音響線信号が取得された観測点Ｐｉｊの位置を指標として複数のサブフレーム音響線信号を加算することにより、各観測点に対する合成音響線信号を生成して主部分音響線信号を合成する。同様に、加算処理部１１４０１−２は、各副領域音響線信号に含まれる音響線信号が取得された観測点Ｑｍｎの位置を指標として複数のサブフレーム音響線信号を加算することにより、各観測点に対する合成音響線信号を生成して副部分音響線信号を合成する。そのため、複数のサブフレーム音響線信号に含まれる同一位置の観測点に対する音響線信号は加算されて合成音響線信号が生成される。なお、以下、加算処理部１１４０１−１と加算処理部１１４０１−２とを区別しないときは、加算処理部１１４０１と表記する。
【００８２】
図１２は、加算処理部１１４０１における合成音響線信号を合成する処理を示す模式図である。上述のとおり、送信イベントに同期して送信振動子列（送信開口Ｔｘ）に用いる振動子を振動子列方向に振動子１つ分だけ異ならせて超音波送信が順次行われる。そのため、異なる送信イベントに基づく対象領域Ｂｘ、Ｃｘも送信イベントごとに同一方向に振動子１つ分だけ位置が異なる。複数のサブフレーム音響線信号を、各サブフレーム音響線信号に含まれる音響線信号が取得された観測点ＰｉｊおよびＱｍｎの位置を指標として加算することにより、全ての対象領域を網羅したフレーム音響線信号が合成される。
【００８３】
また、位置の異なる複数の対象領域にまたがって存在する観測点Ｐｉｊ、Ｑｍｎについては、各サブフレーム音響線信号における音響線信号の値が加算されるので、合成音響線信号は、跨りの程度に応じて大きな値を示す。以後、観測点Ｐｉｊ、Ｑｍｎが異なる対象領域に含まれる回数を「重畳数」、振動子列方向における重畳数の最大値を「最大重畳数」と称する。
【００８４】
また、本実施の形態では、主領域Ｂｘは砂時計形状の領域内に存在する。そのため、図１３（ａ）に示すように、重畳数及び最大重畳数は被検体の深さ方向において変化するので、合成音響線信号の値も同様に深さ方向において変化する。なお、図示はしないが、副領域Ｃｘにおいても、重畳数及び最大重畳数は被検体の深さ方向において変化する。
【００８５】
なお、各サブフレーム音響線信号に含まれる音響線信号が取得された観測点Ｐｉｊ、Ｑｍｎの位置を指標として加算する際に、観測点Ｐｉｊ、Ｑｍｎの位置を指標として重みづけしながら加算してもよい。
【００８６】
合成されたフレーム音響線信号は増幅処理部１１４０２に出力される。
【００８７】
ｉｉ）増幅処理部１１４０２−１、１１４０２−２
上述のとおり、合成音響線信号の値は被検体の深さ方向において変化する。これを補うために、増幅処理部１１４０２は、フレーム音響線信号に含まれる合成音響線信号の合成において、加算が行われた回数に応じて決定した増幅率を各合成音響線信号に乗じる増幅処理を行う。
【００８８】
図１３（ｂ）は、増幅処理部１１４０２における増幅処理の概要を示す模式図である。図１３（ｂ）に示すように、最大重畳数は被検体の深さ方向において変化するので、この変化を補うように、最大重畳数に応じて決定された被検体深さ方向において変化する増幅率が合成音響線信号に乗じられる。これにより、深さ方向における重畳数の変化に伴う合成音響線信号の変動要因は解消され、増幅処理後の合成音響線信号の値は深さ方向において均一化が図られる。
【００８９】
また、重畳数に応じて決定された振動子列方向において変化する増幅率を合成音響線信号に乗じる処理を行ってもよい。振動子列方向において重畳数が変化する場合に、その変動要因を解消し、振動子列方向において増幅処理後の合成音響線信号の値の均一化が図られる。
【００９０】
なお、生成した各観測点に対する合成音響線信号に増幅処理を施した信号をフレーム音響線信号としてもよい。
【００９１】
ｉｉｉ）結合部１１４０３
結合部１１４０３は、フレーム音響線信号を合成するための主部分音響線信号と副部分音響線信号とを結合してフレーム音響線信号を生成する。
【００９２】
結合部１１４０３は、画像化主領域Ｈｘ内の観測点Ｐについては、主部分音響線信号の値を、観測点Ｐにおけるフレーム音響線信号の値とする。一方、追加領域Ｈａ、Ｈｂ内の観測点Ｑについては、以下のようにフレーム音響線信号を作成する。観測点Ｑについて、対応する主部分音響線信号が存在しない場合、副部分音響線信号の値を、観測点Ｑにおいけるフレーム音響線信号の値とする。一方、対応する主部分音響線信号と副部分音響線信号のいずれも存在する場合には、主部分音響線信号および副部分音響線信号のいずれか一方、または両方を用いて観測点Ｑにおけるフレーム音響線信号の値を算出する。本実施の形態においては、主部分音響線信号の値を、観測点Ｑにおけるフレーム音響線信号の値とする。このようにすることで、主領域Ｂｘ内の観測点Ｐとして取得された主領域音響線信号のみに基づいてフレーム音響線信号の値を算出することができ、空間分解能や信号Ｓ／Ｎ比を向上させることができる。ただし、フレーム音響線信号の値の算出方法は、主部分音響線信号と副部分音響線信号とに基づく値、例えば、相加平均、相乗平均や一次結合であってもよい。このようにすることで、フレーム音響線信号の値が急激に変化し、領域の界面において画質の向上が十分とならないことを避けることができる。
【００９３】
生成されたフレーム音響線信号は、データ格納部１０７に出力され保存される。
【００９４】
＜動作＞
以上の構成からなる超音波診断装置１００の動作について説明する。
【００９５】
図１４は、受信ビームフォーマ部１０４のビームフォーミング処理動作を示すフローチャートである。
【００９６】
先ず、ステップＳ１０１において、送信部１０３１は、プローブ１０１に存する複数の振動子１０１ａ中送信開口Ｔｘに含まれる各振動子に超音波ビームを送信させるための送信信号を供給する送信処理（送信イベント）を行う。
【００９７】
次に、ステップＳ１０２において、受信部１０４０は、プローブ１０１での超音波反射波の受信から得た電気信号に基づき受信信号を生成しデータ格納部１０７に出力し、データ格納部１０７に受信信号を保存する。プローブ１０１に存する全ての振動子１０１ａから超音波送信が完了したか否かを判定する（ステップＳ１０３）。そして、完了していない場合にはステップＳ１０１に戻り、送信開口Ｔｘを列方向に移動ピッチＭｐだけ移動させながら送信イベントを行い、完了している場合にはステップＳ２０１に進む。
【００９８】
次に、ステップＳ２１０において、対象領域設定部１０４２は、送信イベントに同期して、仮想送信開口Ｔｖの位置を示す情報に基づき主領域Ｂｘ、副領域Ｃｘを設定する。１回目のループでは、図１０（ａ）に示されるように、初回の送信イベントにおける仮想送信開口Ｔｖから求められる主領域Ｂｘ、副領域Ｃｘが設定される。
【００９９】
次に、観測点同期型ビームフォーミング処理（ステップＳ２２０（Ｓ２２１〜Ｓ２２８））に進む。ステップＳ２２０では、まず、観測点Ｐｉｊ、Ｑｍｎの位置を示す座標ｉｊ、ｍｎを主領域Ｂｘ、副領域Ｃｘ内の最小値に初期化し（ステップＳ２２１、Ｓ２２２）、受信開口設定部１０４３は、列中心が観測点Ｐｉｊ、Ｑｍｎに最も空間的に近接する振動子Ｘｋと合致するよう受信開口Ｒｘ振動子列を選択する（ステップＳ２２３）。
【０１００】
次に、観測点Ｐｉｊについて音響線信号を生成する（ステップＳ２２４）。
【０１０１】
ここで、ステップＳ２２４における、観測点Ｐｉｊについて音響線信号を生成する動作について説明する。図１５は、受信ビームフォーマ部１０４における観測点Ｐｉｊ、Ｑｍｎについての音響線信号生成動作を示すフローチャートである。図１６は、受信ビームフォーマ部１０４における観測点Ｐｉｊ、Ｑｍｎについての音響線信号生成動作を説明するための模式図である。
【０１０２】
まず、ステップＳ２２４１において、送信時間算出部１０４４は、主領域Ｂｘ内に存在する任意の観測点Ｐｉｊおよび副領域Ｃｘ内に存在する任意の観測点Ｑｍｎについて、送信された超音波が被検体中の観測点ＰｉｊまたはＱｍｎに到達する送信時間を算出する。送信時間は、上述したように、観測点Ｐｉｊについては、（１）観測点Ｐｉｊがフォーカス深さより深い場合は、幾何学的に定まる仮想送信開口Ｔｖの中心から送信フォーカス点Ｆを経由して観測点Ｐｉｊに至る経路（４０１＋４０２）の長さを超音波の音速ｃｓで除することにより、（２）観測点Ｐｉｊがフォーカス深さより浅い場合には、幾何学的に定まる、幾何学的に定まる仮想送信開口Ｔｖの中心から送信フォーカス点Ｆに至る経路と観測点Ｐｉｊからフォーカス点に至る経路との差分（４０１−４０２）の長さを、超音波の音速ｃｓで除することにより、算出できる。観測点Ｑｍｎについては、（ａ）幾何学的に定まる仮想送信開口Ｔｖの中心から観測点Ｑｍｎに至る経路（４１１）の長さを超音波の音速ｃｓで除することにより算出した時間、（ｂ）上述の（１）または（２）で算出した送信時間、（ｃ）幾何学的に定まる仮想送信開口Ｔｖの中心から観測点Ｑｍｎと同じ深さの参照点Ｒに至る経路（４１２）の長さを超音波の音速ｃｓで除することにより算出した時間、のうち、（ａ）と、（ｂ）と（ｃ）のうち少なくとも一方を用いた一次結合により算出できる。
【０１０３】
次に受信開口Ｒｘから求められる受信開口Ｒｘ内の受信振動子Ｒｋの位置を示す座標ｋを受信開口Ｒｘ内の最小値に初期化し（ステップＳ２２４２）、送信された超音波が被検体中の観測点Ｐｉｊで反射され受信開口Ｒｘの受信振動子Ｒｋに到達する受信時間を算出する（ステップＳ２２４３）。受信時間は、幾何学的に定まる観測点Ｐｉｊ、Ｑｍｎから受信振動子Ｒｋまでの経路４０３の長さを超音波の音速ｃｓで除することにより算出できる。さらに、送信時間と受信時間の合計から、送信開口Ｔｘから送信された超音波が観測点Ｐｉｊ、Ｑｍｎで反射して受信振動子Ｒｋに到達するまでの総伝播時間を算出し（ステップＳ２２４４）、受信開口Ｒｘ内の各受信振動子Ｒｋに対する総伝播時間の差異により、各受信振動子Ｒｋに対する遅延量を算出する（ステップＳ２２４５）。
【０１０４】
受信開口Ｒｘ内に存在する全ての受信振動子Ｒｋについて遅延量の算出を完了したか否かを判定し（ステップＳ２２４６）、完了していない場合には座標ｋをインクリメント（ステップＳ２２４７）して、更に受信振動子Ｒｋについて遅延量の算出を行い（ステップＳ２２４３）、完了している場合にはステップＳ２２４８に進む。この段階では、受信開口Ｒｘ内に存在する全ての受信振動子Ｒｋについて観測点Ｐｉｊ、Ｑｍｎからの反射波到達の遅延量が算出されている。
【０１０５】
ステップＳ２２４８において、遅延処理部１０４７は、受信開口Ｒｘ内の受信振動子Ｒｋに対応する受信信号の列から、各受信振動子Ｒｋに対する遅延量を差引いた時間に対応する受信信号を観測点Ｐｉｊ、Ｑｍｎからの反射波に基づく受信信号として同定する。
【０１０６】
次に、重み算出部１０４８は、受信開口Ｒｘの列方向の中心に位置する振動子に対する重みが最大となるよう各受信振動子Ｒｋに対する重み数列を算出する（ステップＳ２２４９）。加算部１０４９は、各受信振動子Ｒｋに対応して同定された受信信号に、各受信振動子Ｒｋに対する重みを乗じて加算して、観測点Ｐｉｊ、Ｑｍｎに対する音響線信号を生成し（ステップＳ２２５０）、生成された観測点Ｐｉｊ、Ｑｍｎに対する音響線信号はデータ格納部１０７に出力され保存される（ステップＳ２２５１）。
【０１０７】
次に、図１４に戻り、座標ｉｊ、ｍｎをインクリメントしてステップＳ２２３、Ｓ２２４を繰り返すことにより、主領域Ｂｘ内の座標ｉｊに位置する全ての観測点Ｐｉｊ（図１５中の「・」）、副領域Ｃｘ内の座標ｍｎに位置する全ての観測点Ｑｍｎについて音響線信号が生成される。対象領域内に存在する全ての観測点について音響線信号の生成を完了したか否かを判定し（ステップＳ２２５、Ｓ２２７）、完了していない場合には座標ｉｊ、ｍｎをインクリメント（ステップＳ２２６、Ｓ２２８）して、観測点Ｐｉｊ、Ｑｍｎについて音響線信号を生成し（ステップＳ２２４）、完了した場合にはステップＳ２３０に進む。この段階では、１回の送信イベントに伴う対象領域に存在する全ての観測点Ｐｉｊ、Ｑｍｎについてのサブフレームの音響線信号が生成され、データ格納部１０７に出力され保存されている。
【０１０８】
次に、全ての送信イベントについて、サブフレームの音響線信号の生成が終了したか否かを判定し（ステップＳ２３０）、終了していない場合には、ステップＳ２１０に戻り、観測点Ｐｉｊ、Ｑｍｎの位置を示す座標ｉｊ、ｍｎを、次の送信イベントでの仮想送信開口Ｔｖから求められる対象領域内の最小値に初期化し（ステップＳ２２１、Ｓ２２２）、受信開口Ｒｘを設定（ステップＳ２２３）、音響線信号の作成（ステップＳ２２４）を行い、終了している場合にはステップＳ３０１に進む。
【０１０９】
次に、ステップＳ３０１において、加算処理部１１４０１は、データ格納部１０７に保持されている複数のサブフレーム音響線信号を読み出し、観測点Ｐｉｊ、Ｑｍｎの位置を指標として複数のサブフレーム音響線信号を加算して各観測点Ｐｉｊ、Ｑｍｎに対する合成音響線信号を生成して主部分音響線信号、副部分音響線信号を合成する。次に、増幅処理部１１４０２は、主部分音響線信号、副部分音響線信号のそれぞれに含まれる各合成音響線信号の加算回数に応じて決定された増幅率を各合成音響線信号に乗じ（ステップＳ３０２）、増幅された主部分音響線信号、副部分音響線信号を結合部１１４０３に出力する。次に、結合部１１４０３は、主部分音響線信号と副部分音響線信号とを結合してフレーム音響線信号を生成し（ステップＳ３０３）、超音波画像生成部１０５及びデータ格納部１０７に出力して（ステップＳ３０４）、処理を終了する。
【０１１０】
＜まとめ＞
以上、説明したように本実施の形態に係る超音波診断装置１００によれば、合成開口法により、異なる送信イベントにより生成された同一位置にある観測点についての音響線信号に重ね合わせて合成する。これにより、複数の送信イベントに対して送信フォーカス点Ｆ以外の深度にある観測点においても、仮想的に送信フォーカスを行った効果が得られ空間分解能と信号Ｓ／Ｎ比を向上することができる。
【０１１１】
また、超音波診断装置１００では、追加領域においては、少なくとも幾何学的に定まる仮想送信開口Ｔｖの中心から観測点Ｑｍｎに至る経路（４１１）の長さを超音波の音速ｃｓで除することにより算出した時間に基づいて送信時間を算出する。これにより、フォーカス深さ近傍の観測点や画像化主領域Ｈｘから素子列方向に離れた観測点において遅延処理における誤差を低減し信号Ｓ／Ｎ比を向上させた音響線信号を得ることができる。一方で、主領域Ｂｘに近接する観測点においては主領域Ｂｘと同様の方法で算出した送信時間に基づいて送信時間を算出することにより、主領域Ｂｘと副領域Ｃｘとの境界が顕在化することを避けることができる。したがって、画像化主領域Ｈｘのみを画像化する場合と比べて送信イベント数を追加することなく、かつ、画像化主領域Ｈｘと追加領域Ｈａ、Ｈｂとの境界ならびに追加領域Ｈａ内、Ｈｂ内における画質の急激な変化を抑制しながら、画像化領域を素子列方向に拡大することができる。
【０１１２】
また、超音波診断装置１００では、受信開口設定部１０４３は、列中心が観測点に最も空間的に近接する振動子と合致するよう受信開口Ｒｘ振動子列を選択し、送信イベントに依存せず観測点の位置に基づいて、観測点を中心として対称な受信開口を用いて受信ビームフォーミングを行う。そのため、送信フォーカス点Ｆを横軸方向に変化（移動）させる送信イベントに同期せず、受信開口の位置が一定となり、異なる送信イベントにおいても同一の観測点に対して同一の受信開口にて整相加算を行うことができる。併せて、観測点からの反射波を、観測点から距離が小さい振動子ほど大きな重み数列が適用されることができるので、超音波が伝播距離に依存して減衰することを鑑みても、観測点に対して最も感度よく反射波を受信することができる。その結果、局所的に高い空間分解能と信号Ｓ／Ｎ比を実現できる。
【０１１３】
≪変形例１≫
実施の形態１に係る超音波診断装置１００では、受信開口設定部１０４３は、列中心が観測点に最も空間的に近接する振動子と合致するよう受信開口Ｒｘを選択する構成とした。しかしながら、受信開口Ｒｘの構成は、送信開口Ｔｘから送信された超音波が送信フォーカス点Ｆを経由して対象領域内の観測点で反射され受信開口Ｒｘの受信振動子Ｒｋに到達するまでの総伝播時間を算出して総伝播経路に基づく遅延制御を行なうことで、対象領域内の全ての観測点についての音響線信号を生成するものであればよく、受信開口Ｒｘの構成は適宜変更することができる。
【０１１４】
変形例１では、列中心が送信開口中心と合致する受信開口Ｒｘ振動子列を選択する送信同期型受信開口設定部（以後、「Ｔｘ受信開口設定部」）を備えた点で実施の形態１と相違する。Ｔｘ受信開口設定部以外の構成については、実施の形態１に示した各要素と同じであり、同じ部分については説明を省略する。
【０１１５】
図１７は、Ｔｘ受信開口設定部により設定された仮想受信開口Ｒｖ、受信開口Ｒｘと仮想送信開口Ｔｖ、送信開口Ｔｘとの関係を示す模式図である。変形例１では、仮想受信開口Ｒｖの列中心が送信開口中心と合致するように仮想受信開口Ｒｖが選択される。仮想受信開口Ｒｖの中心軸Ｒｘｏの位置は、仮想送信開口Ｔｖの中心軸Ｔｘｏの位置と同一であり、仮想受信開口Ｒｖは、送信フォーカス点Ｆを中心として対称な開口である。したがって、送信イベントにごとに列方向に移動する仮想送信開口Ｔｖの位置変化に同期して、仮想受信開口Ｒｖの位置も移動する。仮想送信開口Ｒｖのうち、振動子列１０１ａ上の振動子の存在する範囲が、受信開口Ｒｘとなる。
【０１１６】
また、仮想受信開口Ｒｖの中心軸Ｒｘｏ及び仮想送信開口Ｔｖの中心軸Ｔｘｏ上に位置する振動子に対する重みが最大となるよう受信開口Ｒｘの各受信振動子Ｒｋに対する重み数列（受信アポダイゼーション）は算出される。重み数列は、送信開口中心を中心として対称な分布をなす。重み数列の分布の形状は、ハミング窓、ハニング窓、矩形窓などを用いることができ、分布の形状は特に限定されない。
【０１１７】
＜動作＞
図１８は、変形例１に係る超音波診断装置の受信ビームフォーマ部のビームフォーミング処理動作を示すフローチャートである。本フローチャートでは、図１４における観測点同期型ビームフォーミング処理（ステップＳ２２０（Ｓ２２１〜Ｓ２２８））に替えて送信同期型ビームフォーミング処理（ステップＳ４２０（Ｓ４２１〜Ｓ４２８））を行う点にて相違する。ステップＳ４２０以外の処理については、図１４と同じであり、同じ部分については説明を省略する。
【０１１８】
ステップＳ４２０の処理において、先ず、ステップＳ４２１では、Ｔｘ受信開口設定部は、送信イベントに対応して列中心が送信開口中心の列中心と合致する振動子列を仮想受信振動子Ｒｖとして選択して受信開口Ｒｘを設定する。
【０１１９】
次に、ステップＳ２１０で算出した対象領域内の観測点Ｐｉｊ、Ｑｍｎの位置を示す座標ｉｊ、ｍｎを対象領域内の最小値に初期化し（ステップＳ４２２、Ｓ４２３）、観測点Ｐｉｊ、Ｑｍｎについて音響線信号を生成する（ステップＳ４２４）。図１８は、変形例１に係る受信ビームフォーマ部における観測点Ｐｉｊ、Ｑｍｎについての音響線信号生成動作を説明するための模式図である。実施の形態１に関する図１６とは、送信開口Ｔｘと受信開口Ｒｘとの位置関係が異なる。ステップＳ４２４における処理方法については、図１１におけるステップＳ２２４（図１４におけるステップＳ２２４１〜ステップＳ２２５１）と同じである。
【０１２０】
座標ｉｊ、ｍｎをインクリメントしてステップＳ４２４を繰り返すことにより、対象領域内の座標ｉｊ、ｍｎに位置する全ての観測点Ｐｉｊ（図１６中の「・」）、Ｑｍｎについて音響線信号が生成される。対象領域内に存在する全ての観測点Ｐｉｊ、Ｑｍｎについて音響線信号の生成を完了したか否かを判定し（ステップＳ４２５、Ｓ４２７）、完了していない場合には座標ｉｊ、ｍｎをインクリメント（ステップＳ４２６、Ｓ４２８）して、観測点Ｐｉｊ、Ｑｍｎについて音響線信号を生成し（ステップＳ４２４）、完了している場合にはステップＳ２３０に進む。この段階で、１回の送信イベントに伴う対象領域内に存在する全ての観測点Ｐｉｊ、Ｑｍｎのサブフレームの音響線信号が生成されデータ格納部１０７に出力され保存されている。
【０１２１】
＜効果＞
以上説明した、変形例１に係る超音波診断装置では、実施の形態１において示した効果のうち観測点同期型の受信開口に関する部分を除いた効果に変えて、以下の効果を奏する。すなわち、変形例１では、Ｔｘ受信開口設定部は送信イベントに対応して列中心が送信開口中心と合致する振動子列を仮想受信開口Ｒｖとして選択して受信開口Ｒｘを設定する。そのため、仮想受信開口Ｒｖの中心軸Ｒｘｏの位置は、仮想送信開口Ｔｖの中心軸Ｔｘｏの位置と同一であり、送信イベントごとに列方向に移動する送信開口中心の位置変化に同期して、受信開口Ｒｘの位置も変化（移動）する。よって、送信イベントに同期してそれぞれ異なる受信開口にて整相加算を行うことができ、複数の送信イベントにわたって受信時刻は異なるものの、結果としてより一層広い受信開口を用いた受信処理の効果が得られ、広い観測領域で空間分解能を均一にすることができる。また、観測点Ｑｍｎについては、仮想受信開口Ｒｖの中心軸Ｒｘｏの位置が振動子列１０１ａ上に存在することとなるため、最も大きな重み数列が適用される振動子が必ず存在することとなり、受信感度を向上させることができる。
【０１２２】
≪変形例２≫
実施の形態および変形例では、超音波プローブがリニアプローブである場合について説明した。
【０１２３】
しかしながら、本開示の内容は超音波プローブがリニアプローブである場合に限らず、任意の種類のプローブであってよい。
【０１２４】
変形例２に係る超音波診断装置は、超音波プローブとして、振動子が円弧上に配置された、いわゆるコンベックスプローブを用いることを特徴とする。
【０１２５】
変形例２に係る超音波診断装置では、図２０に示すように、画像化主領域Ｈｘは、振動子列の表面たる円弧と、振動子列の一端を通り円弧と直交する直線と、振動子列の他端を通り円弧と直交する直線とで３方向を囲まれた、扇形から中心を同じくする小さな扇形を除いた形状の領域となる。また、追加領域Ｈａ、Ｈｂは、振動子の表面の円弧を延長した仮想的な円弧を一辺とする領域となる。送信イベントごとの、主領域Ｂｘ、副領域Ｃｘについては、実施の形態と同様の方法により設定することができる。
【０１２６】
これにより、コンベックスプローブを用いる場合においても、画像化領域を素子列方向に拡大することができる。
【０１２７】
≪実施の形態に係るその他の変形例≫
（１）実施の形態および各変形例では、主領域Ｂｘ内の観測点Ｐについては、（ａ）観測点Ｐｉｊがフォーカス深さより深い場合は、仮想送信開口Ｔｖの中心から送信フォーカス点Ｆを経由して観測点Ｐｉｊに至る経路の長さを超音波の音速ｃｓで除することにより、（ｂ）観測点Ｐｉｊがフォーカス深さより浅い場合には、仮想送信開口Ｔｖの中心から送信フォーカス点Ｆに至る経路と観測点Ｐｉｊからフォーカス点に至る経路との差分の長さを、超音波の音速ｃｓで除することにより、算出するとした。しかしながら、例えば、主領域Ｂｘ内の観測点Ｐのうち、追加領域Ｈａ内または追加領域Ｈｂ内に存在し、かつ、フォーカス深さより浅い観測点については、副領域Ｃｘ内の観測点Ｑと同様の手法により算出するとしてもよい。また、この場合において、副領域Ｃｘと主領域Ｂｘとの境界のうちフォーカス深さより浅い領域については、境界近傍であっても送信時間Ｔ2の重み付け係数（１−α）を小さいまま、としてもよい。
【０１２８】
あるいは、副領域Ｃｘ内の観測点Ｑについては、受信時間Ｔ1をそのまま送信時間として適用してもよい。これにより、特にフォーカス点から遠い領域、深さの浅い領域について、遅延処理の精度を向上させることができる。
【０１２９】
（２）実施の形態および各変形例では、超音波プローブの振動子の配列を直線状（リニアプローブ）、円弧状（コンベックスプローブ）としたが、これに限られず、任意の配列であってよい。
【０１３０】
（３）実施の形態および各変形例では、画像化主領域と追加領域とからなる領域について超音波画像を生成するものとしたが、本発明はこの場合に限られない。例えば、素子列方向において追加領域よりさらに外側となる領域について、超音波ビームの出射方向を外側に変更して（ステアリング）画像化を行う、としてもよい。この場合においても、追加領域をステアリングにより画像化する場合と比べて送信イベントの回数を削減することが可能となるため、従来技術と比べて１フレームあたりの送信イベント回数を削減することができ、フレームレートの向上に奏功する。
【０１３１】
（４）なお、本発明を上記実施の形態に基づいて説明してきたが、本発明は、上記の実施の形態に限定されず、以下のような場合も本発明に含まれる。
【０１３２】
例えば、本発明は、マイクロプロセッサとメモリを備えたコンピュータシステムであって、上記メモリは、上記コンピュータプログラムを記憶しており、上記マイクロプロセッサは、上記コンピュータプログラムにしたがって動作するとしてもよい。例えば、本発明の超音波信号処理方法のコンピュータプログラムを有しており、このプログラムに従って動作する（又は接続された各部位に動作を指示する）コンピュータシステムであってもよい。
【０１３３】
また、上記超音波診断装置の全部、もしくは一部、また超音波信号処理装置の全部又は一部を、マイクロプロセッサ、ＲＯＭ、ＲＡＭ等の記録媒体、ハードディスクユニットなどから構成されるコンピュータシステムで構成した場合も本発明に含まれる。上記ＲＡＭ又はハードディスクユニットには、上記各装置と同様の動作を達成するコンピュータプログラムが記憶されている。上記マイクロプロセッサが、上記コンピュータプログラムにしたがって動作することにより、各装置はその機能を達成する。
【０１３４】
また、上記の各装置を構成する構成要素の一部又は全部は、１つのシステムＬＳＩ（Ｌａｒｇｅ  Ｓｃａｌｅ  Ｉｎｔｅｇｒａｔｉｏｎ（大規模集積回路））から構成されているとしてもよい。システムＬＳＩは、複数の構成部を１個のチップ上に集積して製造された超多機能ＬＳＩであり、具体的には、マイクロプロセッサ、ＲＯＭ、ＲＡＭなどを含んで構成されるコンピュータシステムである。これらは個別に１チップ化されてもよいし、一部又は全てを含むように１チップ化されてもよい。なお、ＬＳＩは、集積度の違いにより、ＩＣ、システムＬＳＩ、スーパーＬＳＩ、ウルトラＬＳＩと呼称されることもある。上記ＲＡＭには、上記各装置と同様の動作を達成するコンピュータプログラムが記憶されている。上記マイクロプロセッサが、上記コンピュータプログラムにしたがって動作することにより、システムＬＳＩは、その機能を達成する。例えば、本発明のビームフォーミング方法がＬＳＩのプログラムとして格納されており、このＬＳＩがコンピュータ内に挿入され、所定のプログラム（ビームフォーミング方法）を実施する場合も本発明に含まれる。
【０１３５】
なお、集積回路化の手法はＬＳＩに限るものではなく、専用回路または汎用プロセッサで実現してもよい。ＬＳＩ製造後に、プログラムすることが可能なＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）や、ＬＳＩ内部の回路セルの接続や設定を再構成可能なリコンフィギュラブル・プロセッサー（Ｒｅｃｏｎｆｉｇｕｒａｂｌｅ  Ｐｒｏｃｅｓｓｏｒ）を利用してもよい。
【０１３６】
さらには、半導体技術の進歩または派生する別技術によりＬＳＩに置き換わる集積回路化の技術が登場すれば、当然、その技術を用いて機能ブロックの集積化を行ってもよい。
【０１３７】
また、各実施の形態に係る、超音波診断装置の機能の一部又は全てを、ＣＰＵ等のプロセッサがプログラムを実行することにより実現してもよい。上記超音波診断装置の診断方法や、ビームフォーミング方法を実施させるプログラムが記録された非一時的なコンピュータ読み取り可能な記録媒体であってもよい。プログラムや信号を記録媒体に記録して移送することにより、プログラムを独立した他のコンピュータシステムにより実施するとしてもよい、また、上記プログラムは、インターネット等の伝送媒体を介して流通させることができるのは言うまでもない。
【０１３８】
上記実施形態に係る超音波診断装置では、記憶装置であるデータ格納部を超音波診断装置内に含む構成としたが、記憶装置はこれに限定されず、半導体メモリ、ハードディスクドライブ、光ディスクドライブ、磁気記憶装置、等が、超音波診断装置に外部から接続される構成であってもよい。
【０１３９】
また、ブロック図における機能ブロックの分割は一例であり、複数の機能ブロックを一つの機能ブロックとして実現したり、一つの機能ブロックを複数に分割したり、一部の機能を他の機能ブロックに移してもよい。また、類似する機能を有する複数の機能ブロックの機能を単一のハードウエア又はソフトウェアが並列又は時分割に処理してもよい。
【０１４０】
また、上記のステップが実行される順序は、本発明を具体的に説明するために例示するためのものであり、上記以外の順序であってもよい。また、上記ステップの一部が、他のステップと同時（並列）に実行されてもよい。
【０１４１】
また、超音波診断装置には、プローブ及び表示部が外部から接続される構成としたが、これらは、超音波診断装置内に一体的に具備されている構成としてもよい。
【０１４２】
また、上記実施の形態においては、プローブは、複数の圧電素子が一次元方向に配列されたプローブ構成を示した。しかしながら、プローブの構成は、これに限定されるものではなく、例えば、複数の圧電変換素子を二次元方向に配列した二次元配列振動子や、一次元方向に配列された複数の振動子を機械的に揺動させて三次元の断層画像を取得する揺動型プローブを用いてもよく、測定に応じて適宜使い分けることができる。例えば、２次元に配列されたプローブを用いた場合、圧電変換素子に電圧を与えるタイミングや電圧の値を個々に変化させることによって、送信する超音波ビームの照射位置や方向を制御することができる。
【０１４３】
また、プローブは、送受信部の一部の機能をプローブに含んでいてもよい。例えば、送受信部から出力された送信電気信号を生成するための制御信号に基づき、プローブ内で送信電気信号を生成し、この送信電気信号を超音波に変換する。併せて、受信した反射超音波を受信電気信号に変換し、プローブ内で受信電気信号に基づき受信信号を生成する構成を採ることができる。
【０１４４】
また、各実施の形態に係る超音波診断装置、及びその変形例の機能のうち少なくとも一部を組み合わせてもよい。更に上記で用いた数字は、全て本発明を具体的に説明するために例示するものであり、本発明は例示された数字に制限されない。
【０１４５】
さらに、本実施の形態に対して当業者が思いつく範囲内の変更を施した各種変形例も本発明に含まれる。
【０１４６】
≪まとめ≫
（１）実施の形態に係る超音波信号処理装置は、超音波プローブに列設された複数の振動子を選択的に駆動して被検体に超音波送信する送信イベントを複数回繰り返すとともに、各送信イベントに同期して被検体から反射超音波を受波し、受波した反射超音波に基づいて生成される複数のサブフレーム音響線信号からフレーム音響線信号を合成する超音波信号処理装置であって、１回の送信イベントにおいて前記複数の振動子から送信振動子列を選択し、前記送信振動子列から超音波が被検体中で集束するように送信し、各送信イベントに同期して、超音波を送信する送信振動子列が列方向に順次移動するよう選択する送信部と、各送信イベントに同期して前記複数の振動子から受信振動子列を選択し、前記受信振動子列が被検体内から受波した反射超音波に基づいて、前記受信振動子列の振動子各々に対する受信信号列を生成する受信部と、対応するフレーム音響線信号が形成されるべき被検体内の領域として画像化主領域と、前記画像化主領域に対して列方向に隣接する追加領域とを設定する画像化領域設定部と、前記送信イベントごとに、各観測点から得られた反射超音波に基づく前記受信信号列を整相加算して前記サブフレーム音響線信号を生成する整相加算部と、前記整相加算部が生成した前記複数のサブフレーム音響線信号に基づき、前記フレーム音響線信号を合成する合成部とを備え、前記整相加算部は、前記被検体内で超音波が集束するエリアに含まれる領域を含む主対象領域と、前記主対象領域に列方向に隣接する領域の一部または全部が前記追加領域内に存在する場合に前記追加領域内の当該部分である副対象領域とを、前記観測点を含む領域として設定し、前記画像化主領域内かつ前記主対象領域内の観測点と前記副対象領域に含まれる観測点との間で、送信された超音波が観測点に到達する送信時間の算出方法を異ならせることを特徴とする。
【０１４７】
また、実施の形態に係る超音波信号処理方法は、超音波プローブに列設された複数の振動子を選択的に駆動して被検体に超音波送信する送信イベントを複数回繰り返すとともに、各送信イベントに同期して被検体から反射超音波を受波し、受波した反射超音波に基づいて生成される複数のサブフレーム音響線信号からフレーム音響線信号を合成する超音波信号処理方法であって、１回の送信イベントにおいて前記複数の振動子から送信振動子列を選択し、前記送信振動子列から超音波が被検体中で集束するように送信し、各送信イベントに同期して、超音波を送信する送信振動子列が列方向に順次移動するよう選択し、各送信イベントに同期して前記複数の振動子から受信振動子列を選択し、前記受信振動子列が被検体内から受波した反射超音波に基づいて、前記受信振動子列の振動子各々に対する受信信号列を生成し、前記送信イベントごとに、各観測点から得られた反射超音波に基づく前記受信信号列を整相加算して前記サブフレーム音響線信号を生成し、前記整相加算部が生成した前記複数のサブフレーム音響線信号に基づき、前記フレーム音響線信号を合成する超音波信号処理方法であって、前記サブフレーム音響線信号を生成する際に、前記被検体内で超音波が集束するエリアに含まれる領域を含む主対象領域と、前記主対象領域に列方向に隣接する領域の一部または全部が前記追加領域内に存在する場合に前記追加領域内の当該部分である副対象領域とを、前記観測点を含む領域として設定し、前記画像化主領域内かつ前記主対象領域内の観測点と前記副対象領域に含まれる観測点との間で、送信された超音波が観測点に到達する送信時間の算出方法を異ならせることを特徴とする。
【０１４８】
上記構成又は方法によれば、フレームレートを低下させることなく、かつ、空間解像度とＳ／Ｎ比を低下させることなく、追加領域の面積だけ画像化領域を拡大することができる。
【０１４９】
（２）また、上記（１）の超音波信号処理装置は、前記整相加算部は、前記主対象領域に含まれる観測点について、前記観測点の深さが前記フォーカス深さより深い位置の場合は、送信された超音波が前記主対象領域内であって前記フォーカス深さにある基準点に到達するまでの第１時間に、超音波が前記基準点から前記観測点に到達するまでの第２時間を加算して算出した到達時間を、前記観測点の深さが前記フォーカス深さ未満の場合は、前記第１時間から前記第２時間を減算して算出した到達時間を、前記観測点についての前記送信時間とする、としてもよい。
【０１５０】
上記構成によれば、主対象領域内の観測点について、誤差の小さい遅延処理を行うことにより、音響線信号のＳ／Ｎ比を向上させることができる。
【０１５１】
（３）また、上記（１）または（２）の超音波信号処理装置は、前記整相加算部は、前記副対象領域に含まれる観測点について、送信された超音波が前記フォーカス点に最近接する前記送信振動子列上の点から前記観測点に到達するまでの時間を、前記観測点についての前記送信時間とする、としてもよい。
【０１５２】
上記構成によれば、副対象領域内の観測点について、特にフォーカス点から遠い領域において誤差の小さい遅延処理を行うことにより、音響線信号のＳ／Ｎ比を向上させることができる。
【０１５３】
（４）また、上記（１）または（２）の超音波信号処理装置は、前記整相加算部は、前記副対象領域に含まれる観測点について、送信された超音波が前記フォーカス点に最近接する前記送信振動子列上の点から前記観測点に到達するまでの時間をＴ1とし、前記観測点の深さが前記フォーカス深さより深い位置の場合は、送信された超音波が前記主対象領域に含まれ前記フォーカス深さにある基準点に到達するまでの第３時間に、超音波が前記基準点から前記観測点に到達するまでの第４時間を加算して算出した時間、前記観測点の深さが前記フォーカス深さ未満の場合は、前記第３時間から前記第４時間を減算して算出した時間をＴ2としたとき、前記送信時間ＴMを、数式ＴM＝αＴ1＋（１−α）Ｔ2を用いて算出し、αの値は、０より大きく１以下である、としてもよい。
【０１５４】
上記構成によれば、副対象領域内の観測点について、特にフォーカス点から遠い領域において誤差の小さい遅延処理を行うことにより音響線信号のＳ／Ｎ比を向上させることができるとともに、主対象領域と副対象領域との境界で送信時間が大きく変化することを抑止し、画質向上が不十分となることを抑止することができる。
【０１５５】
（５）また、上記（１）または（２）の超音波信号処理装置は、前記整相加算部は、前記副対象領域に含まれる観測点について、送信された超音波が前記フォーカス点に最近接する前記送信振動子列上の点から前記観測点に到達するまでの時間をＴ1とし、前記観測点の深さが前記フォーカス深さより深い位置の場合は、送信された超音波が前記主対象領域に含まれ前記フォーカス深さにある基準点に到達するまでの第３時間に、超音波が前記基準点から前記観測点に到達するまでの第４時間を加算して算出した時間、前記観測点の深さが前記フォーカス深さ未満の場合は、前記第３時間から前記第４時間を減算して算出した時間をＴ2とし、送信された超音波が前記超音波振動子列の列中心から前記観測点と同じ深さで且つ前記超音波振動子列の列中心に最近接する第２基準点に到達するまでの時間とＴ3としたとき、前記送信時間ＴMを、数式ＴM＝αＴ1＋（１−α）βＴ2＋（１−α）（１−β）Ｔ3を用いて算出し、αの値は、０より大きく１以下であり、βの値は、０以上１以下である、としてもよい。
【０１５６】
上記構成によれば、副対象領域内の観測点について、特にフォーカス点から遠い領域において誤差の小さい遅延処理を行うことにより音響線信号のＳ／Ｎ比を向上させることができるとともに、フォーカス深さ近傍、及び、主対象領域と副対象領域との境界で送信時間が大きく変化することを抑止し、画質向上が不十分となることを抑止することができる。
【０１５７】
（６）また、上記（５）の超音波信号処理装置は、前記βの値は、対応する観測点の深さと前記フォーカス深さとの差異が大きくなるほど大きくなる、としてもよい。
【０１５８】
上記構成により、フォーカス深さ近傍において送信時間が大きく変化することを抑止し、画質向上が不十分となることを抑止することができる。
【０１５９】
（７）また、上記（５）または（６）の超音波信号処理装置は、前記αの値は、対応する観測点の深さと前記フォーカス深さとの差異が大きくなるほど小さくなる、としてもよい。
【０１６０】
上記構成により、主対象領域と副対象領域との境界で送信時間が大きく変化することを抑止し、画質向上が不十分となることを抑止することができる。
【０１６１】
（８）また、上記（１）〜（７）の超音波信号処理装置は、前記主対象領域は、前記フォーカス点に最近接する前記送信振動子列上の点と、前記フォーカス点とを結ぶ直線に対して線対称の形状を有する、としてもよい。
【０１６２】
上記構成により、合成部における合成回数が素子列方向においてほぼ均等となるため、音響線信号のＳ／Ｎ比および空間解像度のばらつきを抑止することができる。
【０１６３】
（９）また、上記（８）の超音波信号処理装置は、前記送信振動子列が前記超音波プローブの端部の振動子を含まない送信イベントにおいて、前記主対象領域は、前記フォーカス点と前記送信振動子列の一端とを結ぶ直線と、前記フォーカス点と前記送信振動子列の他端とを結ぶ直線との間に位置する領域である、としてもよい。
【０１６４】
上記構成により、主対象領域における超音波の利用効率を向上させるとともに、合成部における合成回数を最大化させることで、音響線信号のＳ／Ｎ比および空間解像度を向上させることができる。
【０１６５】
（１０）また、上記（８）または（９）の超音波信号処理装置は、前記送信振動子列が前記超音波プローブの端部の振動子を含む送信イベントにおいて、前記振動子列を延長した仮想線上に、前記送信振動子列のうち前記超音波プローブの端部とは異なる端部を一方の端として、前記フォーカス点に最近接する前記送信振動子列上の点を中点とする仮想送信振動子列を仮定したとき、前記主対象領域は、前記フォーカス点と前記仮想送信振動子列の一端とを結ぶ直線と、前記フォーカス点と前記仮想送信振動子列の他端とを結ぶ直線との間に位置する領域である、としてもよい。
【０１６６】
上記構成により、追加領域内に設定される主対象領域において、超音波の利用効率を向上させるとともに、合成部における合成回数を最大化させることで、音響線信号のＳ／Ｎ比および空間解像度を向上させることができる。
【０１６７】
（１１）また、上記（１０）の超音波信号処理装置は、前記追加領域は、列方向における前記フォーカス点の位置が前記追加領域に最も近づく送信イベントにおいて前記仮想送信振動子列のうち前記送信振動子列に含まれない部分を底辺とする領域である、としてもよい。
【０１６８】
上記構成により、追加領域内に副対象領域のみならず主対象領域も設定されるため、音響線信号のＳ／Ｎ比および空間解像度を向上させることができる。
【０１６９】
（１２）また、上記（１）〜（１１）の超音波信号処理装置は、前記追加領域と前記画像化主領域との境界線と、前記フォーカス点と前記送信振動子列の一端とを結ぶ直線との列方向の距離をｍとしたとき、前記第１領域は、前記追加領域を列方向に沿ってｍだけシフトした領域から前記主対象領域との重複領域を除いた領域である、としてもよい。
【０１７０】
上記構成により、追加領域内においても、合成開口法による音響線信号のＳ／Ｎ比および空間解像度の向上効果を最大限享受することができる。
【０１７１】
（１３）また、上記（１）〜（１２）の超音波信号処理装置は、前記整相加算部は、前記画像化主領域と前記追加領域とのそれぞれの範囲と、前記フォーカス点の位置によって定まる前記第１領域の範囲とに基づいて前記副対象領域が空領域であるか否かを判定するための情報を保持し、前記情報を用いて前記副対象領域が空領域であると判定した送信イベントにおいて、前記副対象領域を設定せず前記主対象領域のみを設定する、としてもよい。
【０１７２】
上記構成により、副対象領域が空領域となる場合には第１領域を仮定する必要がなくなり、演算量の低減に奏功する。
【０１７３】
（１４）また、上記（１）〜（１３）の超音波信号処理装置は、前記合成部は、前記画像化主領域に含まれる観測点について、前記主対象領域に対応する前記サブフレーム音響線信号を観測点の位置に基づいて合成することでフレーム音響線信号を生成する、としてもよい。
【０１７４】
上記構成により、主対象領域内については、Ｓ／Ｎ比の高い主対象領域に係る音響線信号に基づいてフレーム音響線信号を生成することができる。
【０１７５】
（１５）実施の形態に係る超音波画像表示方法は、表示部が接続可能であって、上記（１）の超音波信号処理方法を実行する超音波信号処理装置における超音波画像表示方法であって、前記フレーム音響線信号を超音波画像に変換して前記表示部に表示する。
【０１７６】
上記方法により、ユーザに対し視野角の広い超音波診断装置を提供することが可能となる。
【産業上の利用可能性】
【０１７７】
本開示にかかる超音波信号処理装置、超音波診断装置、超音波信号処理方法、プログラム、及びコンピュータ読み取り可能な非一時的な記録媒体は、従来の超音波診断装置の画像化領域の拡大による使用上の視野拡大に有用である。
【符号の説明】
【０１７８】
１０００  超音波診断システム
１００  超音波診断装置
１５０  超音波信号処理装置
１０１  プローブ
１０２  マルチプレクサ部
１０３  送信ビームフォーマ部
１０３１  送信部
１０４  受信ビームフォーマ部
１０４０  受信部
１０４１  整相加算部
１０４２  対象領域設定部
１０４３  受信開口設定部
１０４４  送信時間算出部
１０４５  受信時間算出部
１０４６  遅延量算出部
１０４７  遅延処理部
１０４８  重み算出部
１０４９  加算部
１１４０  合成部
１１４０１  加算処理部
１１４０２  増幅処理部
１１４０３  結合部
１０５  超音波画像生成部
１０６  表示部
１０７  データ格納部
１０８  制御部
【図１】
2019072358.tif 000003
【図２】
2019072358.tif 000004
【図３】
2019072358.tif 000005
【図４】
2019072358.tif 000006
【図５】
2019072358.tif 000007
【図６】
2019072358.tif 000008
【図７】
2019072358.tif 000009
【図８】
2019072358.tif 000010
【図９】
2019072358.tif 000011
【図１０】
2019072358.tif 000012
【図１１】
2019072358.tif 000013
【図１２】
2019072358.tif 000014
【図１３】
2019072358.tif 000015
【図１４】
2019072358.tif 000016
【図１５】
2019072358.tif 000017
【図１６】
2019072358.tif 000018
【図１７】
2019072358.tif 000019
【図１８】
2019072358.tif 000020
【図１９】
2019072358.tif 000021
【図２０】
2019072358.tif 000022
【図２１】
2019072358.tif 000023