(57)【要約】  被検体媒質の不均質性に起因する画質劣化を補償することのできる超音波撮像装置を提供する。  受信ビームフォーマ１０８は、超音波素子アレイ１０５が受信した信号を、２以上のステアリング方向についてそれぞれ整相処理した後合成する。２以上のステアリング方向は、ステアリング方向指示部１１２が指示する。２以上のステアリング方向は、受信合焦点の方向に対して、超音波素子アレイの配列方向に沿って左右にそれぞれ所定の角度をなす２方向を含む。この所定の角度は、ヌル角度であることが好ましい。
2014109392.tif 000002
【特許請求の範囲】
【請求項１】
所定の方向に沿って複数の超音波素子を配列した超音波素子アレイと、前記超音波素子アレイが受信した信号を整相する受信ビームフォーマと、前記受信ビームフォーマの出力する整相出力を用いて画像データを生成する画像処理部とを有し、
前記受信ビームフォーマは、前記超音波素子アレイが受信した信号を、２以上のステアリング方向についてそれぞれ整相処理した後合成する整相合成部と、前記整相合成部に前記２以上のステアリング方向を指示するステアリング方向指示部とを備え、
前記２以上のステアリング方向は、受信合焦点の方向以外の少なくとも２方向を含むことを特徴とする超音波撮像装置。
【請求項２】
請求項１に記載の超音波撮像装置において、前記２方向はそれぞれ、前記超音波素子アレイの指向性プロファイルのヌル角度の方向であることを特徴とする超音波撮像装置。
【請求項３】
請求項１または２に記載の超音波撮像装置において、前記２以上のステアリング方向は、前記受信合焦点の方向をさらに含むことを特徴とする超音波撮像装置。
【請求項４】
請求項２に記載の超音波撮像装置において、前記ステアリング方向指示部は、予め求めておいた前記超音波素子アレイの２以上の前記ヌル角度の方向を前記整相合成部に指示することを特徴とする超音波撮像装置。
【請求項５】
請求項２に記載の超音波撮像装置において、前記ステアリング方向指示部は、前記超音波素子アレイの２以上の前記ヌル角度の方向を計算により求めるヌル角度計算部を有し、前記ヌル角度計算部が求めたヌル角度の方向を前記整相合成部に指示することを特徴とする超音波撮像装置。
【請求項６】
請求項２に記載の超音波撮像装置において、前記ステアリング方向指示部は、前記超音波素子アレイが受信した信号から、２以上の前記ヌル角度の検出するヌル角度検出部を有し、前記ヌル角度検出部が求めたヌル角度の方向を前記整相合成部に指示することを特徴とする超音波撮像装置。
【請求項７】
請求項６に記載の超音波撮像装置において、前記ステアリング方向指示部は、前記超音波素子アレイの２以上の前記ヌル角度の方向を計算により求めるヌル角度計算部をさらに有し、前記ヌル角度検出部は、前記ヌル角度計算部が求めたヌル角度の周辺の予め定めた範囲についてヌル角度の検出を行ってヌル角度を求めることを特徴とする超音波撮像装置。
【請求項８】
請求項６または７に記載の超音波撮像装置において、前記整相合成部は、複数のステアリング方向についてそれぞれ、前記超音波素子アレイの各素子が受信した信号を遅延した後加算した整相信号を生成する遅延加算部を有し、
前記ステアリング方向指示部は、前記遅延加算部の出力する整相信号を用いて、整相信号が極小となるステアリング方向を前記ヌル角度の方向として検出することを特徴とする超音波撮像装置。
【請求項９】
請求項６または７に記載の超音波撮像装置において、前記整相合成部は、前記適応ビームフォーミングを行うことにより、前記２以上のステアリング方向のそれぞれについて適応重みを求める適応整相部と、前記適応重みを用いて前記素子アレイの受信した信号を重み付けする重み付け演算部とを有し、
前記ステアリング方向指示部は、前記重み付け演算部の出力する信号が極小となるステアリング方向を前記ヌル角度の方向として検出することを特徴とする超音波撮像装置。
【請求項１０】
請求項１に記載の超音波撮像装置において、前記整相合成部は、前記超音波素子アレイの各素子が受信した信号を前記ステアリング方向に応じて遅延した後加算することにより整相信号を生成する遅延加算部と、前記２以上のステアリング方向それぞれについての前記整相信号を合成して前記整相出力を得る合成部とを有することを特徴とする超音波撮像装置。
【請求項１１】
請求項１に記載の超音波撮像装置において、前記整相合成部は、前記適応ビームフォーミングを行うことにより、前記２以上のステアリング方向のそれぞれについて適応重みを求める適応整相部と、前記２以上のステアリング方向についての適応重みを合成する重み合成部と、前記合成部の合成した適応重みを用いて前記素子アレイの受信した信号を重み付けした後合成する重み付け演算部とを有することを特徴とする超音波撮像装置。
【請求項１２】
請求項１に記載の超音波撮像装置において、前記整相合成部は、前記適応ビームフォーミングを行うことにより、前記２以上のステアリング方向のそれぞれについて適応重みを求める適応整相部と、前記適応重みを用いて前記素子アレイの受信した信号を重み付けした後合成する重み付け演算部と、前記２以上のステアリング方向それぞれについての前記重み付け演算部の出力を合成する合成部とを有することを特徴とする超音波撮像装置。
【請求項１３】
請求項１記載の超音波撮像装置において、
前記複数の超音波素子に対してアクティブチャンネルを時系列に異なる位置に順次設定し、前記アクティブチャンネルに含まれる前記超音波素子の受信信号を前記受信ビームフォーマに受け渡す制御部をさらに有し、
ある時点の前記アクティブチャンネルの前記２以上のステアリング方向の示す点は、異なる時点の前記アクティブチャンネルについて前記２以上のステアリング方向の示す点の一部と、位置が重なっていることを特徴とする超音波撮像装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波を用いて被検体内の画像を撮像する超音波撮像技術に関する。
【背景技術】
【０００２】
超音波撮像技術とは、超音波（聞くことを意図しない音波、一般的には２０ｋＨｚ以上の高周波数の音波）を用いて人体をはじめとする被検体の内部を非侵襲的に画像化する技術である。一つの例として、医用超音波撮像装置について簡単に説明する。超音波探触子は、超音波を患者の体内に向けて送信し、患者体内から反射されるエコー信号を受信する。受信信号は、超音波探触子および超音波撮像装置本体の一方またはその双方において信号処理を施された後、画像表示部に受け渡され、超音波画像が表示される。より詳しく説明すると、例えば超音波撮像装置本体中の送信ビームフォーマにおいて送信ビームの信号が生成され、送受信分離回路（Ｔ／Ｒ）を経た後に、超音波探触子に送られる。超音波探触子は、超音波を発信する。超音波探触子は、体内からのエコー信号を受信した後、撮像装置本体に信号を伝達する。撮像装置本体中で、受信信号は再び送受信分離回路を経て、受信ビームフォーマにおいて整相処理が施された後、画像処理部に伝達される。画像処理部では、各種フィルタ、スキャンコンバータなど様々な画像処理が実行される。最終的に画像表示部に超音波画像が表示される。【０００３】
このように一般的な超音波診断装置は，送波ビームフォーミング、受信ビームフォーミングおよびバックエンドの画像処理の３つの技術によって構成されている。特に送信時および受信時のビームフォーマは、ＲＦ（高周波）レベルでの信号処理を行うため、ビームフォーマのアルゴリズムや実装アーキテクチャが超音波画像の基本画質を決定する。このため、ビームフォーマは、装置の基幹部である。【０００４】
受信ビームフォーマは、超音波探触子を構成する複数素子の各受信信号（受信データ）に対して、焦点位置と素子の位置との関係に応じて凹面型に遅延量が分布する遅延時間を与え、仮想的に空間のある一点に焦点をあわせた後、受信信号データを加算する。この方法は、遅延加算方式による整相と呼ばれている。この遅延加算方式では、超音波診断装置の複数素子で受信された受信データと診断装置に蓄えられた固定の重みベクトルとを乗算し、重み付けしてから加算する。これは受信ビームフォーマのみならず、送信ビームフォーマでも同様である。【０００５】
一方、超音波撮像装置の基本的な問題として方位方向分解能の制約があることが知られている。超音波の送受信は、有限の開口径を持つアレイによって行われるため、開口部エッジの回折の影響が発生する。無限に長いアレイを用意すれば深さ方向と同様に無限の分解能向上の可能性があるが、現実的には送受アレイの長さ、という装置設計上の物理的な制約があるため、方位方向の分解能向上が妨げられてきた。近年になって、ビームフォーマの遅延加算時に遅延のために用いる前述の固定の重みベクトルを、時系列の送受信データの一つ一つに対して適応的に変化させることにより、より高精細な超音波画像を得る試みが注目を集めている。これにより、ビームフォーミング技術における本質的な課題の一つである方位方向の分解能が格段に向上する可能性がある。【０００６】
特に近年では、移動体通信の分野で発達してきたＭＶＤＲ法（Minimum  Variance  Distortionless  Response;  Capon法）をはじめとした適応信号処理の技術を、受信データのビームフォーマに適用した方位方向の分解能の改善技術が，報告され始めている。これらの適応的手法は，受信データの相関行列に基づき、遅延加算のための重みベクトルの複素成分を適応的に変化させることで実現される。すなわち、従来、重みベクトルは固定値であったが、適応的手法では、受信信号の時間方向のサンプル点ごとに、受信信号を用いて演算により重みベクトルを求め、これを受信信号に乗算する。【０００７】
適応信号処理においては、従来の遅延加算方式と同様に、媒質中の不均質音速分布や、微小体の乱散乱に起因する超音波波面の歪みが課題となる。適応信号処理においても、装置が設定する受信ビームフォーマの焦点は、音速一定かつホモジーニアス（均質）な媒質を仮定して設定されている。そのため、音波伝搬に歪みが存在する場合に、像がぼやけてしまうことや、実際の位置とは違う場所に結像してしまうことがある。従来遅延加算方式においては、古くより波面歪みの補正が課題となっており、相互相関処理を利用した収差補正技術が検討されており、適応ビームフォーマにおいても、媒質内の不均質による超音波画像の画質を大きく変えることが同様に課題となる。【０００８】
ビームフォーマの先行技術としては、例えば、固定のヌル方向を利用した遅延加算ビームフォーマ技術が特許文献１に記載されている。【０００９】
特許文献１の技術では、サイドローブによるアーチファクトを低減するために、情報を取得したい対象物の所定の方向にメインローブを向けて第１受信信号を生成するとともに、上記所定の方向にヌル方向を向けて第２受信信号を生成する。第２受信信号は、対象物の所定の方向の情報をほとんど含まないアーチファクト（ノイズ）の信号であるので、第１受信信号から第２受信信号を差し引くことにより、超音波画像における有用な信号成分の欠落を抑制しつつ、第１受信信号からアーチファクトを低減できる（特に特許文献１の段落００６０）。
【先行技術文献】
【特許文献】
【００１０】
【特許文献１】特開２０１０−１５８３７４
【発明の概要】
【発明が解決しようとする課題】
【００１１】
実際の被検体媒質は不均質であるため、送信した超音波信号が焦点で反射され、焦点の周囲に散乱する。このため、焦点の対象物の情報が焦点の周囲にも存在する。適応ビームフォーマは、指向性が鋭いため、焦点の周囲に存在する超音波信号を取得することができず、超音波画像の画質が劣化するという問題がある。【００１２】
また、被検体媒質が不均質であることに起因して、音速不均一等による波面ひずみが生じ、この波面歪みにより、焦点周囲の媒質からの反射波が受信音波に混入するという問題がある。この周囲の媒質からの反射波は、焦点からの信号と相関性がある雑音信号（相関性雑音）であるために、通常の雑音低減技術で除去することが難しい。この相関性雑音のために、超音波画像の画質が劣化する。【００１３】
本発明の目的は、被検体媒質の不均質性に起因する画質劣化を補償することのできる超音波撮像装置を提供することにある。
【課題を解決するための手段】
【００１４】
本発明の超音波撮像装置の受信ビームフォーマは、超音波素子アレイが受信した信号を、２以上のステアリング方向についてそれぞれ整相処理した後合成する整相合成部と、整相合成部に２以上のステアリング方向を指示するステアリング方向指示部とを備えている。２以上のステアリング方向は、受信合焦点の方向以外の少なくとも２方向を含んでいる。
【発明の効果】
【００１５】
本発明によれば、受信合焦点方向以外の２方向を含む２以上のステアリング方向について整相処理をした信号を合成することにより、相関性雑音を低減することができる。よって、波面歪みに起因する相関性雑音による超音波画像の画質劣化の補償とＳ／Ｎ比の向上を実現できる。
【図面の簡単な説明】
【００１６】
【図１】第１の実施形態の受信ビームフォーマの構成を示すブロック図。
【図２】第１の実施形態のステアリング方向を示す説明図。
【図３】第１の実施形態の相関性雑音低減の作用を説明する説明図。
【図４】（ａ）第１の実施形態の超音波撮像装置の概略構成を示す斜視図、（ｂ）ブロック図。
【図５】第２の実施形態の受信ビームフォーマの構成を示すブロック図。
【図６】図５の受信ビームフォーマの動作を示すフローチャート。
【図７】第３の実施形態の受信ビームフォーマの構成を示すブロック図。
【図８】第４の実施形態の受信ビームフォーマの構成を示すブロック図。
【図９】図８の受信ビームフォーマの動作を示すフローチャート。
【図１０】第５の実施形態の受信ビームフォーマの構成を示すブロック図。
【図１１】図１０の受信ビームフォーマの詳細な構成を示すブロック図。
【図１２】第６の実施形態の受信ビームフォーマの構成を示すブロック図。
【図１３】図１２の整相処理部２０４の詳細な構成を示すブロック図。
【図１４】第７の実施形態の受信ビームフォーマの構成を示すブロック図。
【図１５】第８の実施形態の受信ビームフォーマの構成を示すブロック図。
【図１６】第９の実施形態の受信ビームフォーマの一部構成を示すブロック図。
【図１７】本実施形態の超音波撮像装置のコンソールの斜視図。
【図１８】本実施形態の超音波撮像装置の別の具体例のコンソールと画像表示部の斜視図。
【図１９】本発明による遅延加算処理でのヌル角度検出の結果を示す図（２次元プロット）。
【図２０】本発明による遅延加算処理でのヌル角度検出の結果を示す図（ラインプロファイル）。
【図２１】本発明による適応ビームフォーマでのヌル角度検出の結果を示す図（２次元プロット）。
【図２２】本発明による適応ビームフォーマでのヌル角度検出の結果を示す図（ラインプロファイル）。
【図２３】（ａ）および（ｂ）本発明による異なる被検体について得た画像のコントラストとヌル位置を示す図。
【発明を実施するための形態】
【００１７】
本発明の一実施形態の超音波診断装置について説明する。
【００１８】
（第１の実施形態）
第１の実施形態の超音波診断装置について、図１〜図３を用いて説明する。図１のように、第１の実施形態の超音波診断装置は、所定の方向に沿って複数の超音波素子（超音波振動子）１０５を配列した超音波素子アレイ１０１と、超音波素子アレイ１０１が受信した信号を整相する受信ビームフォーマ１０８と、受信ビームフォーマ１０８の出力する整相出力を用いて画像データを生成する画像処理部１０９とを有する構成である。受信ビームフォーマ１０８は、整相合成部１１３と、ステアリング方向指示部１１２を備えている。整相合成部１１３は、超音波素子アレイ１０１が受信した信号を、２以上のステアリング方向についてそれぞれ整相処理した後合成する処理を行う。ステアリング方向指示部１１２は、整相合成部１１３に２以上のステアリング方向を指示する。ステアリング方向指示部１１２が指示する２以上のステアリング方向は、図２に示すように受信合焦点１０の方向２０以外の少なくとも２方向（図２では、受信合掌点１０の方向２０に対して超音波素子アレイ１０６の配列方向に沿って左右にそれぞれ所定の角度θＬ、θＲをなす２方向）２１，２２を含んでいる。ここでいうステアリング方向とは、所定の超音波素子１０５を中心に撮像対象（被検体）を向く方向であり、超音波素子アレイ１０１の超音波素子１０５の配列方向（すなわち長手方向）と、超音波素子１０５の超音波送受信面の法線とを含む面内の方向である。
【００１９】
２方向２１，２２はそれぞれ、受信合焦点１０の方向２０以外の方向であれば、被検体媒質の不均質性に起因する画質劣化を補償する一定の効果を得ることができる。また、２方向２１，２２がそれぞれ、図３に示すように、超音波素子アレイ１０１の指向性プロファイルのヌル角度θnullの方向である場合には、画質劣化補修の効果が大きく、さらに望ましい。なお、ここでいうヌル角度θnullは、受信信号がゼロまたは極小値を示す角度である。また、２以上のステアリング方向には、受信合焦点１０の方向が含まれていてもよいし、含まれていなくてもよい。
【００２０】
以下、第１の実施形態の超音波撮像装置をさらに具体的に説明する。
【００２１】
図４（ａ），（ｂ）を用いて超音波撮像装置の全体構成についてさらに説明する。図４（ａ）は装置の斜視図、図４（ｂ）は内部の概略構成を示すブロック図である。
【００２２】
図４（ａ）のように、超音波撮像装置は、超音波探触子１０６と、装置本体１０２と、画像表示部１０３と、コンソール１１０を備えている。装置本体１０２内には、図１（ｂ）のように送信ビームフォーマ１０４と、送受信分離回路（Ｔ／Ｒ）１０７と、受信ビームフォーマ１０８と、画像処理部１０９と、これらの動作を制御する制御部１１１が配置されている。受信ビームフォーマ１０８と画像処理部１０９は、上述の図１の構成である。超音波探触子１０６は、図１に示した超音波素子アレイ１０１を備えている。
【００２３】
送信ビームフォーマ１０４は、送信ビーム用の信号を生成する。送信ビーム用信号は、送受信分離回路１０７を経て、超音波探触子１０６に受け渡される。超音波探触子１０６は、超音波素子アレイ１０１の超音波素子１０５から超音波を被検体１００の体内に向けて送信する。体内で反射されたエコー信号は、超音波探触子１０６の超音波素子アレイ１０１により受信される。受信信号は、再び送受信分離回路１０７を経て受信ビームフォーマ１０８において整相演算処理等が施される。
【００２４】
受信ビームフォーマ１０８は、上述のように整相合成部１１３と、ステアリング方向指示部１１２を備えている。ステアリング方向指示部１１２は、２以上のステアリング方向を整相合成部１１３に指示する。２以上のステアリング方向には、図２のように、受信合焦点１０の方向２０に対して、受信合焦点１０の超音波素子アレイ１０６の配列方向に沿ってある所定の０°よりも大きい角度θＬ、θＲをなす２方向２１，２２が少なくとも含まれている。受信合焦点１０の方向２０がさらに含まれていてもよい。所定の角度θＬ、θＲが、図３の超音波素子アレイの左右のヌル角度θnullの方向であればさらに望ましい。なお、ここでは便宜上、２つのステアリング方向を、受信合焦点１０の方向に対して左右の1つずつの方向に設定した例について説明するが、２以上のステアリング方向は、受信合焦点１０の方向に対してどのような方向に設定してもよい。特に、超音波トランスデューサが２次元アレイの場合は、左右という概念はなく、受信合焦点以外のいかなる任意の２以上の方向に設定してもよい。
【００２５】
整相合成部１１３は、整相処理部２０４と合成部２０６とを備えている。整相処理部２０４は、超音波素子アレイ１０１の複数の超音波素子１０５の出力する受信信号をそれぞれ、ステアリング方向２１の点１１に対して整相した後、加算する。ステアリング方向２２の点２２についても同様に整相加算する。他のステアリング方向がある場合も同様に整相加算する。整相処理部２０４の整相処理は、どのような処理方法であってもよく、例えば遅延加算処理や適応整相処理を用いることができる。
【００２６】
合成部２０６は、整相処理部２０４から複数のステアリング方向２１，２２等についてのそれぞれの整相加算結果を受け取り、加算処理等によりこれらを合成し、整相信号ｙ(n)を出力する。
【００２７】
合成部２０６の出力（整相信号ｙ(n)）は、画像処理部１０９に受け渡され、各種フィルタ、スキャンコンバータなど様々な画像処理が実行され、超音波画像が生成される。超音波画像は、画像表示部１０３に受け渡され、表示される。
【００２８】
このように、本発明では、受信合焦点１０の方向２０に対して、超音波素子アレイ１０６の配列方向に沿って所定の角度θＬ、θＲをなす２方向２１，２２について整相処理した信号を合成することにより、受信合焦点１０（受信すべき点、受信焦点）からの信号Ｓに対して、相関性を示す雑音（相関性雑音Ｎｃ）を低減した整相信号ｙ(n)が得られる。よって、Ｓ／Ｎ比の大きい整相信号ｙ(n)を得ることができ、超音波画像の画質を向上させることができる。
【００２９】
これについてさらに説明する。被検体媒質は不均質であるため、送信した超音波信号が受信合焦点１０で反射され、焦点の周囲の点に散乱するため、受信合焦点１０の対象物の情報が受信合焦点１０の周囲にも存在する。また、被検体媒質が不均質であることに起因して、音速不均一等による波面ひずみが生じ、この波面歪みにより、受信合焦点１０周囲の媒質からの反射波が受信音波に混入する。この周囲の媒質からの反射波は、受信合焦点１０からの信号と相関性がある雑音信号（相関性雑音）である。このため、焦点の周囲に散らばってしまった焦点の情報を再現しつつ、焦点周囲の媒質からの不要な相関性雑音を低減する必要がある。
【００３０】
なお、反射波の伝搬過程だけでなく、音波が送信され受信合焦点１０にたどり着くまでの送信音波の伝搬過程においても、波面歪は発生するため、前述の不要な相関性雑音を増強させてしまう要因となることは明らかである。
【００３１】
そこで本発明では、受信合焦点１０の方向２０に対して、超音波素子アレイ１０６の配列方向に沿ってそれぞれ所定の角度θＬ、θＲをなす２方向２１，２２について整相処理した信号を合成する。例えば、超音波素子アレイ１０１が、図３のように角度方向の指向性プロファイルを有する場合、本発明においては、両サイドに存在するヌル角度θnullの方向についてそれぞれ整相処理した信号と、必要に応じて、中心部すなわちメインローブ（超音波伝搬方向、θ＝０°）の方向について整相処理した信号とを合成する。受信合焦点１０の情報を表す信号をＳ、相関性雑音をＮｃ、無相関性雑音（ホワイトノイズ）をＮｕと表すと、メインローブ方向２０にステアリング方向を向けた場合の整相処理部２０４の出力信号は、（Ｓ＋Ｎｃ＋Ｎｕ）で表される。この場合、信号のS／N比はＳＮＲ＝Ｓ／（Ｓ＋Ｎｃ＋Ｎｕ）となる。相関性雑音Ｎｃは、生体不均質によって進行方向を乱された相関性雑音であり、図２のように受信合焦点１０とは異なる方向から超音波素子アレイ１０１に到来する。また、無相関性雑音Ｎｕは、熱雑音や電気的な雑音などであり、超音波探触子アレイ１０１の受信信号に混在する方向性を持たない信号である。
【００３２】
また、図２のように所定の角度θＬ（例えば、左側の一つ目のヌル角度θ_left_null-1）の方向２１にステアリング方向を向けた場合の整相処理部２０４の出力信号は、（Ｓ_left＋Ｎｃ_left＋Ｎｕ）で表される。また、所定の角度θＲ（例えば、右側の1つめのヌル角度θ_right_null-1）の方向にステアリング方向を向けた場合の整相処理部２０４の出力は、（Ｓ_right＋Ｎｃ_right＋Ｎｕ）と表される。ただし、Ｓ_leftおよびＳ_rightは、受信合焦点１０の情報を表す信号である。
【００３３】
ここで、角度θＬと角度θＲとは、左右方向について異なる方向を向いているため、Ｓ_leftとＳ_right、Ｎｃ_leftとＮｃ_rightは、信号の位相・時相が異なるのみならず、方向性が逆向きの成分を有する。すなわち、Ｓ_leftとＳ_right、Ｎｃ_leftとＮｃ_rightは時間空間の双方に非対称性を有し、加算した場合にお互いを打ち消しあう信号となっている。よって、整相合成部１１３の合成部２０６において、（Ｓ_left＋Ｎｃ_left＋Ｎｕ）と（Ｓ_right＋Ｎｃ_right＋Ｎｕ）とを加算することにより、相関性雑音Ｎｃの成分をキャンセルすることができる。また、特に、角度θＬおよびθＲをヌル角度θnullに設定した場合には、整相処理部２０４の出力信号に含まれるＳ_leftとＳ_rightの信号は、信号全体に占める割合は少ないため、左右にそれぞれ所定の角度θＬ、θＲをなす２方向２１，２２について整相処理した信号を合成した場合、受信合焦点１０の情報を表す信号Ｓ_leftとＳ_rightの低減効果よりも、相関性雑音Ｎｃの低減効果の方が大きくなる。よって、結果的に、ＳＮＲの向上に寄与する割合をより大きくできる。
【００３４】
さらに、図２および図３に示すように、左右のヌルビーム方向に加えてメインローブ方向２０をステアリング方向に設定し、それぞれのステアリング方向についての整相処理部２０４の出力を合成することによって、相関性のノイズＮｃ’（＝Ｎｃ_left＋Ｎｃ_righ＋Ｎ）を大きく低減することができる。また、受信合焦点１０の情報を表す信号Ｓ’（＝Ｓ＋Ｓ_left＋Ｓ_right）は、メインローブ方向２０のＳと比較して、相関性雑音Ｎｃ’ほど大きく低減しないため、ＳＮＲを向上させる効果がさらに大きくなる。
【００３５】
一方、無相関性雑音Ｎｕは、方向性を持たず、ｎ回加算することによって、（１／√ｎ）倍になるため、複数のステアリング方向の整相処理部２０４の出力を加算することにより一意に減少する。
【００３６】
以上のように、もともとのメインローブ方向２０の整相処理部２０４の出力信号（Ｓ＋Ｎｃ＋Ｎｕ）と比べて、左右のヌルビーム方向２１，２２に加えてメインローブ方向２０の整相処理部２０４の出力を合成した信号（S'＋Nc'＋Nu'） = (1/3)×｛（Ｓ＋Ｎ＋Ｒ）+（Ｓ_left＋Ｎｃ_left＋Ｎｕ）＋（Ｓ_right＋Ｎｃ_right＋Ｎｕ）｝は、ＳＮＲが向上する（式（１））。なお、ここでは、３方向についての整相処理部２０４の出力をいずれも重み１／３として加算しているが、異なる重みで重み付けすることも可能である。
【００３７】
【数１】
2014109392.tif 000003
【００３８】
なお、整相合成部１１３の整相処理方法は、どのような方法でもよいが、適応整相処理を用いることにより、受信合焦点１０の情報を表す信号Sの成分により指向性を持たせることができる。すなわち、所望信号成分の合成結果(１／３)×（Ｓ＋Ｓ_left＋Ｓ_right）の中の所望信号成分Ｓをほとんど減少させることなく、相関性雑音Ｎｃの減少効果をより高めることが可能になるため好ましい。
【００３９】
このように第１の実施形態では、受信合焦点の方向以外の少なくとも２方向のステアリング方向の整相処理結果を合成することにより、相関性雑音の時空間での非対称性を利用して、波面ひずみに起因する相関性雑音を低減することができる。しかも、受信合焦点１０の周囲からの情報も収集できるため、焦点で反射されて周囲に散らばった情報も収集できる。よって、被検体の媒質が不均質であっても画質劣化を防止できるとともに、波面歪による相関性雑音の影響も受けにくく、頑健性が高く、しかも、比較的処理負荷の小さい超音波撮像装置が得られる。
【００４０】
（第２の実施形態）
図５および図６を用いて第２の実施形態の超音波撮像装置について説明する。図５は、受信ビームフォーマ１０８の構成を示すブロック図であり、図６は、その動作を示すフローチャートである。
【００４１】
第２の実施形態では、第１の実施形態の整相合成部１１３内の整相処理部２０４が、遅延加算処理により整相処理を行うためにＰ個の遅延加算部２０４−１〜２０４−Ｐを備えている。Ｐ個は、ステアリング方向指示部１１２が指定可能な最大のステアリング方向の数である。また、ステアリング方向指示部１１２は、超音波素子アレイ１０１の条件および超音波照射条件に基づきヌル角度を計算するヌル角度計算部５０１を備え、計算により求めたヌル角度の方向をステアリング方向として整相処理部２０４に指定する。ヌル角度の方向は、受信合焦点１０の方向以外の角度を、少なくとも２つ、を含むようにする。整相処理部２０４の遅延加算部２０４−１〜２０４−Ｐは、指定されたステアリング方向についてそれぞれ遅延加算による整相を行う。
【００４２】
また、第２の実施形態では、制御部１１１内のアクティブチャンネル設定部１１１ａが、超音波素子アレイ１０１の有限口径の一部にアクティブチャンネル２０１，２０２，２０３等を設定する。受信ビームフォーマ１０８はアクティブチャンネルの受信信号について受信ビームフォーミング処理を行う。すなわち、１本の送信超音波ビームに対する受信エコーを受信した超音波素子アレイ１０１を構成する超音波素子１０５のうち、所定の範囲のＫ個の素子１０５をアクティブチャンネル２０１とし、受信ビームフォーマ１０８は、アクティブチャンネル２０１の受信信号を用いて超音波伝搬方向の１本の画像データ（ラスタ：整相出力ｙ（ｎ））を生成する。図５のように、超音波素子アレイ１０１上で位置を少しずつずらしながらアクティブチャンネル２０２、アクティブチャンネル２０１、アクティブチャンネル２０３を順次構成し、それぞれのアクティブチャンネル２０２、２０１，２０３についてラスタを生成し、その結果を並べたものが超音波画像となる。
【００４３】
ヌル角度計算部５０１は、設定されているアクティブチャンネル２０１等についてヌル角度を計算する。
【００４４】
以下、図６のフローを用いてさらに具体的に説明する。
【００４５】
まず、制御部１１１は、そのときの撮像条件に応じた、プローブ条件・超音波照射条件・ヌル角度選択条件を示す情報を含む制御信号を生成し、ヌル角度計算部５０１に出力する（ステップ６１）。
【００４６】
ヌル角度計算部５０１は、制御信号１１１の出力を受け取り（ステップ６２）、制御信号に含まれる情報をパラメータとして、以下の式（３）などを用いてヌル角度を計算する（ステップ６３）。式（３）は送信のみのヌル角度を算出するものであるが、受信に関しても同様の計算によりヌル角度を算出することができる。また、これら送信／受信の重ね合わせとしてのヌル角度を計算することも可能である。
【００４７】
ヌル角度の計算は、超音波の周波数、送信波形などの撮像条件と、使用しているアクティブチャンネル２０１〜２０３に含まれる超音波素子１０５の素子数（アクティブチャンネル素子数）・送信／受信の口径、チャンネル間隔（超音波素子１０５の間隔）、周波数特性によって第１義的なヌル角度を確定できる。すなわち、制御部１１１から伝達された制御信号を用いてヌル角度を計算することができる。
【００４８】
例えば、超音波素子１０５が方形ピストン型である場合、超音波素子アレイ１０１のアクティブチャンネル２０１によって送受される超音波ビームの指向性は、以下の式（２）によって計算できる。Ｄ（θ）は、ステアリング角θにおける超音波ビームの強度である。
【数２】
2014109392.tif 000004
【００４９】
式（２）において、θは、ステアリング角、aは、アクティブチャンネル（２０１〜２０３）の口径、ｋは、波数である。ｋは、音速ｃを超音波の周波数ｆで除したものでありｋ＝ｃ／ｆで表される。
【００５０】
式（２）の大きさがゼロであるような角度θnullが、この場合のヌル角度であり、式（３）から解析的に求めることができる。
【数３】
2014109392.tif 000005
【００５１】
式（３）を見るとわかるように、ヌル角度θnullは、一つではなく複数であり、図３のように、θ＝０°（超音波照射方向＝送信焦点の方向）を中心に左右に対称にｎ個ずつ存在する、と言う特徴を持つ。そのため、ヌル角度計算部５０１には制御信号として、ヌル角度を選択する条件も制御部１１１から入力されている。ヌル角度選択条件は、左右にそれぞれ存在するヌル角度θnull-1、θnull-2・・・θnull-nのうち、整相合成部１１３が処理に使用するヌル角度を選択する条件である。その条件は例えば、n=＋１とn=−１の2本のヌルを使用する、もしくは、n=±１、±２の４本を利用する、のように選択を指示するものである。なお、ｎの符号が、＋である場合は、θ＞０°すなわち、θ＝０°よりも右側に位置するヌル角度であることを示し、ｎの符号が−である場合、θ＜０°すなわち、θ＝０°よりも左側に位置するヌル角度であることを示す。
【００５２】
なお、n=±１の時にもっとも効果が高いため、n=±１の処理部のみを有し、ヌル選択を行わず、必ずn=±１の２本を使用するようにしてもよい。
【００５３】
ヌル角度計算部５０１は、計算により求めたヌル角度のうちヌル角度選択条件で選択したヌル角度の方向を、ステアリング方向として、整相処理部２０４に出力する。ステアリング方向の指定する情報の形式としては、例えば、θ＝１０°、−１０°のような物理的な角度情報であっても良い。また、整相処理部２０４の遅延加算部２０４−１・・・２０４−Ｐをアクティブ（ON）にするための制御信号と上述の角度情報を兼用する信号を直接遅延加算部２０４−１・・・２０４−ｎに受け渡してもよい。この方法は実用上より好ましい形態である。例えば、遅延加算部２０４−１にθＬ＝１０°を指定して遅延加算処理を指示し、遅延加算部２０４−ｎにθＲ＝−１０°を指定して遅延加算処理を指示する。
【００５４】
整相処理部２０４の遅延加算部２０４−１・・・２０４−Ｐは、ヌル角度計算部５０１から指定されたヌル角度の方向をステアリング方向とした場合の整相信号を生成し、出力する（ステップ６５）。この時、正面方向（θ＝０°：受信合焦点の方向２０）についても、整相信号を生成することが望ましい。すなわち、θ＝１０°、−１０°がヌル角度である場合、遅延加算部２０４−１，２０４−２，２０４−３においては、正面方向θ＝０°とθ＝１０°とθ＝−１０°の３方向について遅延加算処理を行い、その結果を合成部２０６に出力する。
【００５５】
合成部２０６は、遅延加算部２０４−１〜２０４−Ｐからヌル方向と正面方向についての整相信号の入力を受ける（ステップ６６）。合成部２０６ではこれらの入力を加算することにより合成し出力する（ステップ６７）。加算に際しては、例えば、正面方向θ＝０°とθ＝１０°とθ＝−１０°の３方向それぞれの単純加算（加算結果には１／３を乗じることが好ましい）でもよいし、ステアリング方向毎に、整相信号に所定の重みを乗じたのちに加算するなどしてもよい。例えば正面方向の整相信号の比重を大きくしたり、逆にヌル方向の整相信号に大きな重みを乗算するなどして合成結果を調整したりすることができる。
【００５６】
合成部２０６によって演算された合成信号は、画像処理部１０９に出力される。画像処理部１０９は、合成信号を受け取って（ステップ６８）、従来の整相結果と同様に処理を行い、１本の画像データ（ラスタ）を生成する。そして、アクティブチャンネル２０１，２０２，２０３を少しずつずらしながら、それぞれラスタを得て、ラスタを並べることにより超音波画像に加工する。超音波画像は、画像表示部１０３に表示する。以上をもって本発明の第２の実施形態のフローが完了する。
【００５７】
上述してきた図６のフローの中で、制御部１１１のヌル角度計算部５０１への制御信号の情報に変化がない場合、すなわちプローブ条件・超音波照射条件・ヌル角度選択条件に変更が施されない限りは、遅延加算部２０４−１〜２０４−Ｐへのヌル角度情報の出力までのステップ（ステップ６１〜６５）の動作は同じである。よって、図６のように、受信ビームフォーマ１０８は、ステップ６７において、１サンプルの受信信号（１回のサンプリングの受信信号のセット）について、合成部２０６が整相信号を画像処理部１０９に出力したならば、ステップ６５に戻って、次の１サンプルの受信信号を処理する動作を１サンプルごとに繰り返す。なお、実用上は、ステップ６５において次の１サンプルの受信信号の処理は、前の１サンプルの合成信号の出力がステップ６７で終了する前に始まっていても良い。すなわち、上述のように図６のフローのプロセスは多段にわたっているため、あるサンプル時刻での処理が次のプロセスに移った時点で、次のサンプルの処理を始めても、構成上および実用上全く差し支えがない。
【００５８】
また、上述してきた説明では、制御部１１１からの制御信号の示す条件が変更になる度に、ヌル角度の計算をする構成であったが、この構成に限らず、ヌル角度計算部５０１にメモリ５０２を配置し、制御信号の示す条件（プローブ条件・超音波照射条件・ヌル角度選択条件）の組み合わせごとのヌル角度の計算結果を対応づけたテーブル（ヌル角度ＬＵＴ（LUT：Look  Up  Table））を予め格納しておくことも可能である。これにより、ヌル角度計算部５０１は、制御信号の示す条件に対応するヌル角度をヌル角度ＬＵＴから読み出して遅延加算部２０４−１〜２０４−Ｐに設定することができるため、演算量を低減することができる。
【００５９】
ヌル角度ＬＵＴは、設定可能なすべての条件（プローブ条件・超音波照射条件・ヌル角度選択条件）について、予めヌル角度を求めて作成したものをメモリ５０２に格納することも可能であるが、これに限らず、一部の条件についてのみヌル角度ＬＵＴを作成してメモリ５０２に格納しておき、ヌル角度ＬＵＴにない条件の場合には、ヌル角度計算部５０１が計算により求めるように構成してもよい。
【００６０】
また、ヌル角度計算部５０１が計算により求めたヌル角度を、その時の制御信号の示す条件に対応づけて、メモリ５０２に格納していくことにより、ヌル角度ＬＵＴを逐次生成していく構成にすることも可能である。
【００６１】
なお、図５に示したように、整相処理部２０４には、Ｐ個の遅延加算部２０４−１〜２０４−Ｐが用意されている。このとき、Ｐ個の遅延加算部２０４−１〜２０４−Ｐにそれぞれ、予めステアリング角度に対応する遅延時間のセットを割り当てておくことも可能である。例えば、Ｐ番目の遅延加算部２０４−Ｐには、最大のステアリング角度の遅延時間が割り当てておく。これにより、ヌル角度計算部５０１は、求めたステアリング角度に対応する遅延加算部を遅延加算部２０４−１〜２０４−Ｐの中から選択して、演算動作を指示するだけで、そのステアリング方向についての整相処理を指示することができる。この方法に限らず、Ｐ個の遅延加算部２０４−１〜２０４−Ｐが、ヌル角度計算部５０１に指定されたステアリング方向に基づいて、その都度遅延時間セットを計算する構成（On  the  Fly演算）であってもよい。
【００６２】
上述してきたように第２の実施形態では、遅延加算処理による整相結果を受信合焦点方向以外の少なくとも２方向のステアリング方向について得て、これらを合成することにより、相関性雑音の時空間の非対称性を利用して、波面ひずみに起因する相関性雑音を低減することができる。また、焦点で反射されて周囲に散らばった情報も収集できる。よって、被検体の媒質が不均質であっても画質劣化を防止できるとともに、波面歪による相関性雑音の影響も受けにくく、頑健性の高く、しかも、比較的処理負荷の小さい超音波撮像装置が得られる。
【００６３】
なお、第２の実施形態の超音波撮像装置の他の構成は、第１の実施形態と同様であるので説明を省略する。
【００６４】
（第３の実施形態）
第３の実施形態として、超音波撮像装置を図７を用いて説明する。
【００６５】
図７のように、第３の実施形態では、受信ビームフォーマ１０８にフレームメモリ７０１とフレーム加算部７０２を備え、開口合成処理を行う。なお、整相合成部１１３には、合成部２０６は備えられていない。他の構成およびその動作は、第２の実施形態と同様である。
【００６６】
本発明における遅延加算部２０４の出力信号は、公知の開口合成処理におけるマルチルック（多方向）受信データと同一である。よって、公知の開口合成処理における、複数送信間のマルチルック受信データの合成時に本発明を適用することにより、本発明の２以上のステアリング方向についての整相合成の処理を開口合成処理によって実現することができる。
【００６７】
具体的には、図７のフレームメモリ部７０１には、遅延加算部２０４−１〜２０４−Ｐで超音波送信ごとに生成されたＰ本のデータが逐次蓄えられる。つぎにフレーム加算部７０２ではこの送信ごとに蓄えられたデータを複数送信間の同一点を見込む角度ごとに加算処理することにより、開口合成処理を行う。このような開口合成処理を行うことにより、ある撮像点に対し、複数の送信方向から得られた受信データを重ね合わせることが可能になり、高解像度・高ＳＮ比・高フレームレートの超音波画像を得ることができる。
【００６８】
フレーム加算部７０２で加算処理するデータは、第２の実施形態の合成部２０６で合成されるデータ（複数方向の遅延加算後データ）と同様のデータである。すなわち、ステアリング方向指示部１１２からの２以上のステアリング方向を合成角度情報としてフレーム加算部７０２へ入力する。フレーム加算部７０２は、指定された２以上のステアリング方向についてのデータをフレームメモリ７０１から読み出して加算処理する。これにより、開口合成処理と、ヌル方向のデータ加算（２以上のステアリング方向の加算）を同時に行うことができる。
【００６９】
（第４の実施形態）
第４の実施形態の超音波撮像装置について図８、図９を用いて説明する。
【００７０】
第４の実施形態では、ステアリング方向指示部１１２がヌル角度検出部２０７を備え、超音波探触子アレイ１０１の実際の受信信号の整相信号を用いて、ヌル角度を検出する。これを可能にするために、整相合成部１１３は、整相処理部２０４にＱ個の遅延加算部２０４−１〜２０４−Ｑを備えるとともに、Ｑ個の遅延加算部２０４−１〜２０４−Ｑの整相出力を格納するメモリ部２０５を備えている。
【００７１】
さらに、ヌル角度検出部２０７は、第２の実施形態で説明したヌル角度計算部５０１を備え、計算により求めたヌル角度周辺についてヌル角度のスキャンを行い、ヌル角度を検出する構成である。
【００７２】
他の構成は、第２の実施形態と同様であるので説明を省略する。
【００７３】
以下、第４の実施形態の受信ビームフォーマ１０８の動作について図９のフローチャートを用いて説明する。図９のフローにおいて図６のフローと同様のステップについては、同じ符号を付している。
【００７４】
第２の実施形態では、プローブ条件等の超音波照射条件からヌル角度を計算により求める構成であったため、装置の設定や超音波照射条件など装置側の条件でヌル角度が決定されていた。しかしながら、実際の受信信号の指向性は、生体における不均質な音波伝搬の影響を受けており、受信ビームの指向性は、受信信号のサンプル毎によって変化する。そこで、第４の実施形態では、ヌル角度検出部２０７が受信信号のサンプルごとに、実際の受信信号に基づいてヌル角度を検出する。
【００７５】
制御部１１１では、第２の実施形態と同様に、そのときの撮像条件に応じた、プローブ条件・超音波照射条件・ヌル角度選択条件を示す制御信号を生成し出力する（ステップ６１）。ヌル角度計算部５０１は、これらの制御信号を受け取り、ヌル角度を計算により求める（ステップ６２，６３）。これらのステップ６１〜６３の詳細は、第２の実施形態で説明した通りである。ヌル角度計算部２０７は、ヌル角度計算部５０１がヌル角度を計算により求めたならば、メモリ部２０５に整相信号のデータを要求する。
【００７６】
一方、整相合成部の遅延加算部２０４−１〜２０４−Ｑは、予め定められたθminからθmaxの角度範囲における所定の角度ずつ異なるＱ個のステアリング方向について、それぞれ遅延加算処理を行い、整相信号を出力する（ステップ９２）。メモリ部２０５は、θminからθmaxの角度範囲のステアリング角度ごとの整相信号を格納する（ステップ９３）。θminからθmaxの角度範囲は、物理的には例えば、−６０°〜６０°のステアリング角度である。
【００７７】
メモリ部２０５は、ステップ９１におけるヌル角度検出部２０７の要求に応じて、θminからθmaxの角度範囲のＱ個のステアリング方向についての整相信号をヌル角度検出部２０７に出力する（ステップ９４）。
【００７８】
ヌル角度検出部２０７は、θminからθmaxの角度範囲のＱ個のステアリング方向についての整相信号を受け取って、θminからθmaxの角度範囲のうち、所定の角度範囲θ１〜θ２の整相信号をスキャンし、最小（極小）の整相出力をヌル角度θnull_scanとして抽出し、出力する（ステップ９５、９６）。これを（θ１，θ２）のセットの数だけ繰り返す。
【００７９】
ここでθ１〜θ２の角度範囲は、ヌル角度計算部５０１が計算により求めたヌル角度を含む予め定めた大きさの角度範囲であり、その角度の間に実際のヌル角度があると想定される角度範囲である。すなわち、（θ１，θ２）のセットは、ヌル角度計算部５０１が求めたθnullに対して、θ１＜θnull＜θ２となる予め定めた大きさの角度範囲である。また、（θ１，θ２）のセット数は、制御信号に含まれるヌル選択条件によって決定されたヌル角度θnullの数と同数である。これについては、第２の実施形態で説明した通りである。
【００８０】
これらの処理により、ヌル角度検出部２０７は、ヌル角度選択条件に基づいた数と同数の、実際のヌル角度θnull_scanを検出する。
【００８１】
ヌル角度検出部２０７は、検出したヌル角度θnull_scanを、メモリ部２０５に指示する。メモリ部２０５は、遅延加算部２０４に使用するヌル角度情報を伝達する。検出したヌル角度θnull_scanの伝達形式は、第2の実施形態と同様であり、例えば、θ＝１０°、−１０°のような物理的な角度情報で伝達する。
【００８２】
メモリ部２０５は、ヌル角度検出部２０７から伝達されたヌル角度θnull_scanと正面方向（θ＝０°）のステアリング方向の整相信号を選択し、合成部２０６に出力する。すなわち、θ＝１０°、−１０°がヌル角度である場合、メモリ部２０５は、正面方向θ＝０°とθ＝１０°とθ＝−１０°の３方向整相信号を出力する。なお、θ＝０°は含まなくても良い。
【００８３】
合成部２０６は、メモリ部２０５からヌル角度θnull_scanと正面方向（θ＝０°）のステアリング方向の整相信号を受け取り（ステップ９７）、加算することにより合成して、画像処理部１０９に出力する（ステップ９８）。加算に際しては、単純加算でもよいし、整相信号毎に所定の重みを乗じたのちに加算してもよい。
【００８４】
画像処理部１０９は、合成部２０６の出力する合成された整相信号を受け取り、これを用いて、第２の実施形態と同様に、超音波画像を生成し、画像表示部１０３に表示する（ステップ９９）。以上をもって本発明の第４の実施形態の処理フローが完了する。
【００８５】
なお、受信信号の一サンプルについて、ステップ９８において合成部２０６による合成信号の出力が終了したならば、ステップ９１，９２にそれぞれ戻って、次のサンプルの受信信号について処理を行う。
【００８６】
以上のように、第４の実施形態では、実際の受信信号の指向性のヌル角度を検出することができるため、生体における不均質な音波伝搬の影響によるヌル角度の変化に対応してヌル角度を設定することができる。よって、被検体の媒質が不均質であっても画質劣化を防止できるとともに、波面歪による相関性雑音の影響も受けにくく、頑健性の高い超音波撮像装置が得られる。
【００８７】
なお、本実施形態では、ヌル角度計算部により、ヌル角度を計算により求め、計算により求めたヌル角度を含む所定の角度範囲θ１〜θ２についてのみ、ヌル角度検出部２０７がヌル角度のスキャンを行えばよい。このため、ヌル角度検出部２０７の演算量を低減することができる。ただし、本発明は、この構成に限定されるものではなく、ヌル角度検出部２０７がヌル角度計算部５０１を備えず全角度範囲についてヌル角度のスキャンを行ってもよい。
【００８８】
また、ヌル角度計算部５０１は、第２の実施形態で説明したようにメモリ５０２を備え、メモリ５０２内のヌル角度LUTを参照してヌル角度を求める構成であっても構わない。
【００８９】
（第５の実施形態）
第５の実施形態の超音波撮像装置について、図１０、図１１を用いて説明する。
【００９０】
第５の実施形態の超音波撮像装置は、第２の実施形態と同様であるが、整相処理部２０４が、適応整相処理を行う点が第２の実施形態とは異なっている。すなわち、整相処理部は、図１０に示すように、遅延回路５１１と、適応整相部５１２を備えている。遅延回路５１１は、ｋ個の超音波素子１０５からなるアクティブチャンネル２０１等のｋ個の受信信号を受信合焦点１０の位置に応じて遅延する回路である。なお、受信合焦点１０は超音波送信方向（０°方向）だけでなく、送信方向から少しだけ傾いた角度方向に複数用意してもよい（受信複ビーム処理）。その場合、その傾いた角度ごとに整相合成部１１３を並列に用意しておく。
【００９１】
適応整相部５１２は、図１１のように、行列演算部と、適応ビームステアリング部３０１とを備え、ヌル角度計算部５０１により指定されたステアリング方向についてそれぞれ適応重みｗを求める。適応合成部５１６は、適応重みｗを用いて、遅延回路５１１の出力する遅延後受信信号を整相し、合成する。なお、ヌル角度計算部５０１の動作は、第２の実施形態と同様である。
【００９２】
このように、第５の実施形態では、受信ビームフォーマ１０８が適応整相処理を行うため、所望の信号成分Ｓをほとんど減少させることなく、相関性雑音Ｎｕをより効果的に低減することができるという効果が得られる。
【００９３】
以下、適応整相処理について具体的に説明する。アクティブチャンネル２０１を構成する超音波素子１０５の数（チャンネル数）がＫ個であれば、あるスナップショット時刻nにおける、遅延回路５１１の遅延後受信データは、超音波素子１０５の出力xｋ(n)を用いて、下記式（４）の左辺のベクトルx(n)として表すことができる。
【数４】
2014109392.tif 000006
【００９４】
この遅延後受信データx(n)は、適応処理部５１２−１〜５１２−Ｐにそれぞれ入力される。適応処理部５１２−１〜５１２−Ｐは、ヌル角度計算部５０１に指定されたステアリング方向ごとに、スナップショット時刻nにおけるＫ個のチャンネルごとの重み値wK(n)から構成される適応重みベクトルwp(n)=[wp_1(n)，wp_2(n)，…，wp_K(n)]Tを生成する。すなわち、ステアリング方向がＰ個である場合には、w1(n)、w2(n)・・・wP(n)が生成される。
【００９５】
適応合成部５１６は、適応重みベクトルw1(n)、w2(n)・・・wP(n)を合成した後、遅延後受信データx(n)と内積演算（重み付け演算）することにより、あるスナップショット時刻nにおける適応ビームフォーマ出力y(n)を得る。
【００９６】
適応整相部５１２の適応処理部５１２−１〜５１２−Ｐは、図１１のように行列演算部３００と適応ビームステアリング部３０１を備えている。適応処理部５１２−１〜５１２−Ｐは、それぞれ行列演算部３００と適応ビームステアリング部３０１を備えていてもよいし、ひと組の行列演算部３００と適応ビームステアリング部３０１を適応処理部５１２−１〜５１２−Ｐで共有していもよい。
【００９７】
遅延後受信データベクトルx(n)は、行列演算部３００において、下式（５）により空間共分散行列R(n)を求められる。R(n)を求める際には，x(n)をそのままの実信号として用いてもよいし、ヒルベルト変換やベースバンド変調などを施して、複素数のデータに変換されたものを用いてもよい。ここでは、より一般性のある形態として、x(n)を複素データξ(n)に変換し、式（５）のように空間共分散行列R(n)を求める場合を例に説明する。式（５）のR(n)は、式（６）であらわされる複素ベクトルξ(n)と、その（複素）共役転置ベクトルξＨ(n)との積のアンサンブル平均をとったものである。
【００９８】
【数５】
2014109392.tif 000007
【数６】
2014109392.tif 000008
【００９９】
なお、式（５）において、アンサンブル平均数Nは，対象スナップショットにおけるξ(n)の前後Sサンプルずつの合計N=2S+1点として、式（５）の最右辺のように一様平均とすることができる。また、時間方向の平均化の仕方はその他にも、台形重みなど、時間方向各サンプルに任意の重みを乗算してから加算平均を取る方法をとってもよい。行列演算部３００から出力された空間共分散行列R(n)は、次に適応ビームステアリング部３０１に入力される。
【０１００】
空間共分散行列R(n)を受け取った適応ビームステアリング部３０１では、MVDR法を用いて重みベクトルwp(n)を計算する。ヌル角度計算部５０１から指定されたステアリングベクトルapは式（７）で表される。
【数７】
2014109392.tif 000009
【０１０１】
上式（７）において、pはステアリングベクトルの本数であり、総数Ｐとして、０＜ｐ＜Ｐ＋１を満たす整数である。式（７）のように、ステアリングベクトルapは、アクティブチャンネルの個数と等しいＫ個のベクトル要素（０〜（Ｋ−１））をもつ方向ベクトルであり、受信周波数fp，および超音波素子１０５の表面の法線ベクトル方向とステアリングベクトルとのなす角（以下、ステアリング角度と呼ぶ）(θp,φp)の関数で表される。
【０１０２】
θpは、法線ベクトルからの開き角、φpは、超音波素子１０５の配列方向からの旋回角度である。超音波探触子１０６の超音波素子１０５の配列が、１次元の直線（リニア）アレイであった場合、ステアリング角度は式（７）の最終式で表される。ここで、λpは周波数fpに対応する音波の波長、dは超音波素子１０５の素子中心同士の間隔（素子ピッチ）である。
【０１０３】
上記のようなステアリングベクトルapの方向についての、ＭＶＤＲ法による重みベクトルwp(n)は、ここでは式（８）により算出される。よって、適応処理部５１２−１〜５１２−Ｐでそれぞれ異なるステアリングベクトルapについて重みベクトルwp(n)を演算することにより、ステアリングベクトルapの本数Ｐの分だけの適応重みベクトルw1(n)〜wP(n)を得ることができる。
【数８】
2014109392.tif 000010
【０１０４】
式（８）において、R(n)は，式（５）により生成した、時間方向のあるスナップショットnにおける空間共分散行列であり、上付き添え字の−１は逆行列を表す。
【０１０５】
なお、受信合焦点１０はアクティブチャンネル２０１の中心軸上、すなわち、アクティブチャンネル数が偶数の場合は、アクティブチャンネル２０１を構成するＫ個の超音波素子１０５のうちＫ／２番目と（Ｋ＋２）／２番目の素子の間の中点を通る法線（アクティブチャンネル面と直交する法線）上に位置するとする。また、アクティブチャンネル数が奇数の場合は、アクティブチャンネル２０１を構成するＫ個の超音波素子１０５のうち、（Ｋ＋１）／２番目の超音波素子１０５の中心を通るような法線（アクティブチャンネル面と直交する法線）上に位置するとする。
【０１０６】
図１１には、ステアリング方向の総数Ｐ＝３のケースを示しているため、適応重みベクトルw1(n)〜wP(n)が３つの例を示しているが、実際には、ヌル角度計算部５０１により指定された個数の適応重みベクトルが適応ビームステアリング部３０１から出力される。
【０１０７】
遅延回路２０４において遅延処理を施しているため、例えばリニアスキャンのある１時刻nの受信データベクトルx(n)は、法線ベクトルの方向であるθ＝0°方向に波面が揃ったデータとなっている。よって、式（７）において、θp＝0°とすると、ステアリングベクトルapは、ap=[1，1，….1]Tであり、焦点方向についての適応重みベクトルw(n)が求められる。焦点方向についての適応重みベクトルw(n)は、従来のＭＶＤＲ法による適応重みベクトルw(n)と一致する。
【０１０８】
適応ビームステアリング部３０１から出力されたP本の適応重みベクトルw1(n)〜wP(n)は、適応合成部５１６に入力される。図１１のように、適応合成部５１６は、P本の適応重みベクトルw1(n)〜wP(n)を重み合成部３０６において加算し、算術平均を取り、式（９）のように合成重みベクトルwsum(n)を算出する。例えば、図１１のようにＰ＝３個の適応重みベクトルの合成重みベクトルwsum(n)は、wsum(n)={w1(n)+w2(n)+w3(n)}/3，となる。
【数９】
2014109392.tif 000011
【０１０９】
なお、この合成重みベクトルwsum(n)を求める前段に、固定アポダイゼーション乗算部３０５を配置し、重みベクトルwp(n)のそれぞれに対して、固定のアポダイゼーションを乗算してもよい。たとえば、θ= 0°方向の適応重みベクトルwP(n)の値を大きくし、それ以外の方向は値が小さくなるような、分布を持たせた固定アポダイゼーションbpを乗算することもでき、これは式（１０）の演算によって実現される。
【数１０】
2014109392.tif 000012
【０１１０】
例えば、図１１の構成においては、固定アポダイゼーション乗算部３０５において、b1,b2,b3の固定アポダイゼーション重みが用意されている。この場合、重み合成部から出力される合成重みは、wsum(n)={b1w1(n)+b2w2(n)+b3w3(n)}/3となる。
【０１１１】
次に、合成重みベクトルwsum(n)は、適応合成部５１６内の重み付け演算部３０７に入力される。重み付け演算部３０７は、乗算部３０７１を含み、合成重みベクトルwsum(n)と、遅延回路２０４からバイパス線２０７を通して伝達された遅延後の受信データベクトルx(n)との間の内積演算を式（１１）のように行い、整相出力y(n)を得る。具体的には、乗算部３０７１で、重みベクトルおよび遅延後受信データベクトルのそれぞれのチャンネル要素間(1〜K)の積を計算する。加算部３０７２は、重み付け演算部３０７の乗算部３０７１で得られたK個の積の総和を計算し、適応合成部５１６の最終出力（整相出力y(n)、スカラー値）とする。
【数１１】
2014109392.tif 000013
【０１１２】
この重み付け演算部３０７の整相出力y(n)は、画像処理部１０９に出力される。以上の処理を１ラスタを構成するＮサンプルのうち最初のサンプルｎ＝１から最終サンプルｎ＝Ｎまで繰り返すことで重み付け演算部３０７は、順次整相出力y(1),y(2),…y(N)を順次画像処理部１０９に出力する。式（１１）によって得られる１ラスタの整相出力ｙ(n)は、アクティブチャンネル２０１をアクティブチャンネル２０２，２０３へと超音波素子アレイ１０１上をシフトさせながら、それぞれ取得される。画像処理部１０９においては、超音波探触子１０６の走査方式などに対応したスキャンコンバータにより、全ラスタを並べ２次元画像を生成する。また、各種フィルタ処理等の様々なバックエンド画像処理や計測アプリケーションの演算がなされる。最終的に、画像表示部１０３に超音波画像や計測アプリケーションの演算結果が表示される。
【０１１３】
このように、第５の実施形態では、適応整相処理を用いることにより、所望の信号成分Ｓの低減を抑制しつつ、相関性雑音Ｎｕのより効果的に低減することができるという効果が得られる。
【０１１４】
なお、第５の実施形態において、上述した行列演算部３００における演算アルゴリズムの別の手法の一つとして、サブアレイ行列を用いた空間平均演算を行うことも可能である。空間平均演算におけるサブアレイ行列R∧SUBlは、式（１２）のように、部分空間ベクトルξ∧l(n)（式（１３））の積で表される。部分空間ベクトルξ∧l(n)は、アクティブチャンネルＫ個に対する遅延後受信データ（ここでは一般化された複素信号ベクトルξ(n)で表記しているが、代わりに実信号ベクトルx(n)でもよい。）のうち、一部（Ｌ素子分）の成分を抜き出したベクトルである。よって、部分空間ベクトルの総数はＫ−Ｌ＋１個（０＜ｌ（小文字Ｌ）＜Ｋ−Ｌ＋１）である。
【０１１５】
【数１２】
2014109392.tif 000014
【数１３】
2014109392.tif 000015
【０１１６】
このサブアレイ行列の主対角成分を空間共分散行列R(n)の主対角成分に一致させて１サンプルずつずらしてゆくと，Ｋ−Ｌ＋１個のサブアレイ行列の空間平均処理となり、式（１４）のサブアレイ空間共分散行列R∧(n)が得られる。このサブアレイ空間共分散行列R∧(n)を適応ビームステアリング部３０１で演算する際に、上記式（８）のR(n)と置き換えて用いることにより、式（１５）のように重みベクトルwp(n)を演算することができる。なお、この場合行列演算部３００の出力はＬ×Ｌのサイズであり、また、重みベクトルwp(n)を構成する要素の数はＬとなる。
【０１１７】
【数１４】
2014109392.tif 000016
【数１５】
2014109392.tif 000017
【０１１８】
また、空間平均法のもう一つの例として、前向き後向き空間平均法を用いてもよい。この場合、後向きのサブアレイ行列R〜SUBl(n)は、式（１６）のように、後向き部分空間ベクトルξ〜l(n)の積で求めることができる。なお、後向き部分空間ベクトルは、式（１７）で表される。また、式（１８）に示すように、前向き空間平均と同様の演算によって、後向きのサブアレイ空間行列R〜(n)を式（１８）のように計算することができる。次に、式（１９）のように、前述の前向き空間平均の場合のサブアレイ空間行列R∧(n)と算術平均することによって、最終的に前向き／後向きのサブアレイ空間共分散行列RFB(n)を求めることができる。前向き空間平均の場合と同様、この前向き／後向きのサブアレイ空間共分散行列RFB(n)を、適応ビームステアリング部３０１での演算において式（８）のR(n)と置き換えて用いることにより、式（２０）のように重みベクトルwp(n)を演算することができる。なお、この場合もやはり、行列演算部３００の出力はＬ×Ｌのサイズであり、また、重みベクトルwp(n)を構成する要素の数はＬとなる。
【０１１９】
【数１６】
2014109392.tif 000018
【数１７】
2014109392.tif 000019
【数１８】
2014109392.tif 000020
【数１９】
2014109392.tif 000021
【数２０】
2014109392.tif 000022
【０１２０】
空間平均法を用いて式（１５）もしくは式（２０）で計算された複数個（Ｐ個）の重みベクトルwpは、空間平均法を用いない場合と同様、適応合成部５１６に伝達される。
【０１２１】
ここで、空間平均法を用いた場合に、重みベクトルwpの要素数はＬである。よって、最終的に内積演算部３０７での演算を行うためには、Ｋ個の成分を持つξ(n)からＬ個の成分を持つベクトルg(n)を作るような演算（式（２１））を行うブロックが追加で必要になる。例えば、図１１に示すように、バイパス線２０７の途中であって重み付け演算部３０７の前段に、次元圧縮部３０８を配置し、次元圧縮部３０８が、Ｋ個の要素を持つξ(n)からＬ個の要素を持つベクトルg(n)を生成するように構成することができる。なお、次元圧縮部３０８は、合成部２０６の内部に配置することも可能である。
【０１２２】
空間平均法を用いた場合の整相出力y(n)は式（２２）で表される。
【数２１】
2014109392.tif 000023
【数２２】
2014109392.tif 000024
【０１２３】
このように、行列演算部３００において、サブアレイ行列の空間平均処理を行うことにより、超音波受信信号に含まれる相関性の雑音を抑制する効果が得られる。よって、サブアレイ行列の空間平均処理を、焦点の周辺の複数の適応重みベクトルを合成する構成に組み合わせることによって、より雑音の影響の少ない超音波画像を得るができる。
【０１２４】
以上、第５の実施形態では、適応ビームフォーミング手法の一例としてＭＶＤＲを例にとって説明したが、適応ビームステアリング部３０１において適用されるアルゴリズムは、行列演算部３００で計算される空間共分散行列を用いるアルゴリズムであればよく、ＭＭＳＥ法やＡＰＥＳ法、空間共分散行列とその固有値・固有ベクトルを利用するＥｉｇｅｎｓｐａｃｅ−ＭＶ法（ＥＳＭＶ、ＥＩＢＭＶ）、ＥＳＰＲＩＴ法、ＭＵＳＩＣ法など他のいかなる手法でもよい。
【０１２５】
（第６の実施形態）
第６の実施形態の超音波撮像装置を図１２および図１３を用いて説明する。第６の実施形態の超音波撮像装置は、第４の実施形態と同様に、整相処理部２０４の整相信号を用いて実際のヌル角度をヌル角度検出部２０７が検出する構成であるが、整相処理部２０４が適応整相処理を行うという点で第４の実施形態とは異なっている。
【０１２６】
具体的には、図１２に示すように、整相処理部２０４に、遅延回路５１１と、適応整相部６１２を備えている。適応整相部６１２は、Ｑ個の適応処理部６１２−１〜６１２−Ｑを備えている。遅延回路５１１は、第５の実施形態と同様の構成および動作である。
【０１２７】
適応処理部６１２−１〜６１２−Ｑは、第４の実施形態と同様に、θminからθmaxの角度範囲のＱ個のステアリング方向について整相出力を求める。適応処理部６１２−１〜６１２−Ｑは、図１３に示すように第５の実施形態の適応処理部５１２−１〜５１２−Ｐと同様に、行列演算部３００と適応ビームステアリング部３０１により、それぞれ適応重みベクトルw1(n)〜wQ(n)を求めるだけでなく、適応処理部６１２−１〜６１２−Ｑは、それぞれ乗算部３０７１と加算部３０７２とを含む。適応処理部６１２−１の乗算部３０７１は、重みベクトルw1(n)と、遅延回路５１１からバイパス線２０７を通して伝達された遅延後の受信データベクトルx(n)との間の内積演算を行い、加算部３０７２は、それを加算して整相出力y1(n)を出力する。同様に、適応処理部６１２−２〜６１２−Ｑは、整相出力y2(n)〜yQ(n)を出力する。なお、固定アポダイゼーション乗算部３０５を配置し、重みベクトルw1(n)に対して、固定のアポダイゼーションを乗算してもよい。
【０１２８】
メモリ部２０５およびヌル角度検出部２０７の動作は、第４の実施形態と同様であるので説明を省略する。
【０１２９】
第６の実施形態では、ステアリング方向について適応整相処理により求めた整相信号を合成するため、波面ひずみに起因する相関性雑音を低減する効果が大きいことができる。また、適応整相処理により求めた整相信号を用いてヌル角度を検出するため、ヌル角度の検出精度も高いという効果がある。
【０１３０】
（第７の実施形態）
第７の実施形態の超音波撮像装置を図１４を用いて説明する。第７の実施形態は、第４の実施形態の図４の構成と、第６の実施形態の図１２の構成とを組み合わせた受信ビームフォーマ１０８を有する。この受信ビームフォーマ１０８は、図８の遅延加算によりＱ個のステアリング方向について整相信号を得る整相処理部２０４とその整相信号を格納するメモリ部２０５を備え、ヌル検出部２０７は、遅延加算により整相処理されたＱ個の整相信号に基づいてヌル角度を検出する。一方、この受信ビームフォーマ１０８は、図１２の適応処理により整相信号を得る整相処理部２０４’を備え、ヌル角度検出部２０７により指定されたＰ個のヌル角度のステアリング方向についてのそれぞれ整相信号y1(n)〜yP(n)を適応整相部６１２が適応整相処理により求める。合成部２０６は、整相信号y1(n)〜yP(n)を加算して合成し、整相信号y(n)を求める。
【０１３１】
本実施形態では、遅延加算処理により求めた整相信号により、ヌル角度を検出し、ヌル角度検出部２０７が指定したステアリング方向についてのみ適応整相処理を行えばよいため、第６の実施形態よりも少ない演算量で、ステアリング方向について適応整相処理により求めた整相信号を合成することができる。よって、少ない演算量で、波面ひずみに起因する相関性雑音の低減効果を高めることができる。
【０１３２】
なお、図１４の適応整相部６１２および合成部２０６を、図１０の適応整相部５１２および適応合成部５１６に置きかえることも可能である。
【０１３３】
（第８の実施形態）
第８の実施形態の超音波撮像装置について図１５を用いて説明する。
【０１３４】
第８の実施形態では、第５の実施形態の図１０の受信ビームフォーマ１０８において開口合成を行う構成である。
【０１３５】
図１５を用いて、本発明の第８の実施形態の超音波撮像装置について具体的に説明する。第３の実施形態と同様に、第８の実施形態では、図１５のフレームメモリ部７０１に、超音波送信ごとに生成されたＲ本のデータが逐次蓄えられている。フレーム加算部７０２は、この送信ごとに蓄えられたデータを複数送信間の同一点を見込む角度ごとに加算処理することにより、本発明の２以上のステアリング方向についての整相合成処理と開口合成処理とを同時に行う。開口合成処理を行うことにより、ある撮像点に対し、複数の送信方向から得られた受信データを重ね合わせることが可能になり、高解像度・高ＳＮ比・高フレームレートの超音波画像を得ることができる。
【０１３６】
第８の実施形態は、図１５のように、第３の実施形態の図７の構成とは、フレームメモリ７０１の前段の整相処理が異なる。第３の実施形態においては複数の遅延加算部２０４−１〜２０４−Ｐが整相合成部１１３に設けられたが、第８の実施形態においては、整相合成部１１３が、第５の実施形態の図５の遅延回路５１１と適応整相部５１２を備えた整相処理部２０４をＲ個備えている。
【０１３７】
フレームメモリ部７０１には、Ｒ個の整相処理部２０４が超音波送信ごとに生成した適応整相処理したＲ本の整相信号が逐次蓄えられる。フレーム加算部７０２は、ステアリング方向指示部１１２からの２以上のステアリング方向を合成角度情報として受け取り、指定された２以上のステアリング方向についてのデータをフレームメモリ７０１から読み出して加算処理する。これにより、開口合成処理と、ヌル方向のデータ加算（２以上のステアリング方向の加算）を同時に行うことができる。
【０１３８】
第８の実施形態では、ステアリング方向について適応整相処理により求めた整相信号を合成するため、第３の実施形態とくらべて、波面ひずみに起因する相関性雑音を低減する効果が大きい。また、適応整相処理により求めた整相信号を用いてヌル角度を計算することができるため、ヌル角度の検出精度も高いという効果がある。
【０１３９】
（第９の実施形態）
第９の実施形態の超音波撮像装置について図１６を用いて説明する。
【０１４０】
第９の実施形態は第５から第８の実施形態における適応合成部の別の実施態様を示すものである。第５から第８の実施形態を通じて、適応合成部の内部は図１１で示される第５の実施形態で説明した構成であっても、第９の実施形態の図１６の構成であっても良い。
【０１４１】
図１１に示した第５の実施形態においては、適応合成部５１６は、重み合成部３０６で複数の適応重みベクトルw1(n),w2(n),w3(n)を合成し、得られた合成重みwsum(n)によって遅延後受信データx(n)の整相加算処理を行う構成であった。この演算は、線形演算であるため、重み合成と整相加算処理の順序を逆にしてもよい。そこで、第９の実施形態の適応合成部２０６では、図１６に示すように、複数の乗算部（３０７１−１、３０７１−２、３０７１−３）と、乗算後の要素（遅延後データ）の和を計算して合成前整相出力を得る複数の加算部３０７２−１、３０７２−２、３０７２−３とを複数の適応重みベクトルw1(n),w2(n),w3(n)ごとに配置している。乗算部３０７１と加算部３０７２それぞれ１つずつのセットは、適応重みベクトルごとの内積演算部を一つ構成しており、これらのセットが複数用意されることで、内積演算部３０７を構成している。乗算部３０７１−１、３０７１−２、３０７１−３それぞれバイパス入力２０７１、２０７２、２０７３により遅延後受信データx(n)を入力する。内積演算部３０７の後段に、出力合成部５００を配置する。
【０１４２】
このような構成により、内積演算部３０７において、バイパス入力２０７１を重みベクトルw1(n)を用いて整相加算し、バイパス入力２０７２を重みベクトルw2(n)を用いて整相加算し、バイパス入力２０７３を重みベクトルw3(n)を用いて整相加算することにより、各ステアリングベクトルに対応した合成前の複数の整相出力y1(n),y2(n),y3(n)を算出する（下式（２３））。その後段で出力合成部５００において、各ステアリングベクトルに対応した整相出力y1(n),y2(n),y3(n)の算術平均をとり、式（２４）のように、合成整相出力ysum(n)を得る。
【０１４３】
【数２３】
2014109392.tif 000025
【数２４】
2014109392.tif 000026
【０１４４】
なお、図１６に示すように、本実施形態においても固定アポダイゼーション乗算部３０５を付加することができる。図１６においては、固定アポダイゼーション乗算部３０５は、内積演算部３０７の前段に配置している。しかしながら、固定アポダイゼーション乗算部３０５の位置は、図１６の位置には限定されない。線形演算であることから明らかなように、内積演算部３０７と出力合成部５００の間に固定アポダイゼーション乗算部３０５を付加しても、さらには，内積演算部中の乗算部３０７１−１〜３と加算部３０７２−１〜３の間に固定アポダイゼーション乗算部３０５を付加してもよく、どちらの場合も、最終的な合成整相出力ysum(n)は式（２５）で表すことができる。
【数２５】
2014109392.tif 000027
【０１４５】
以上のように、第５の実施形態と第８の実施形態は、重みを合成してその後に内積演算を行うか、内積演算後の整相出力を合成するかの違いであるかの違いで実施形態が異なる。その一方で最終的な出力は線形演算のため同一で有り、第５の実施形態の最終的な整相出力である（２２）式のy(n)と、第８の実施形態の最終的な合成整相出力ysum(n)は等しい値となる。
【０１４６】
（第９の実施形態）
図１７に第１〜第８の超音波撮像装置のコンソール１１０の一例を示す。上述してきた実施形態を実現するため、超音波診断装置のコンソール１１０に、ステアリング方向の本数などを変更するための操作部として、目盛りをともなったつまみ部１００１，１００２等を配置することができる。また、上述の各実施形態で説明したヌル角度を計算で求めるか、検出するか等のモードの切り替えなどを行うための操作部として、スイッチ部１００４をコンソール１１０に配置できる。これによって、操作者は実際の超音波画像を見ながら、整相合成処理の各種パラメータを変更し、被検体１００ごとに最適な条件で撮像や診断を行うことができる。また、設定した値を画像表示部１０３の表示領域の一部１１０４に表示する構成にすることもできる。
【０１４７】
また、コンソールのモードの切り替え部は、異なる探触子での切り替えや、撮像条件、撮像シーケンスに対応して切り替えられるようになっていても良い。例えば、リニア探触子、コンベックス探触子、セクタ探触子、２次元アレイ探触子、機械式３Ｄ撮像探触子、などの探触子に対応して異なるステアリング方向のセットを適用して画像を生成することが出来るようなモードを切り替えられるような切り替え部であっても、また、例えば、送信・受信フォーカス、送信／受信周波数、フレームレート、パラレルビームフォーミング本数、組織ハーモニックス、造影コントラストなどの撮像条件、撮像シーケンスごとに異なるステアリング方向のセットを適用して画像を生成することが出来るような切り替え部であってもよい。
【０１４８】
図１８には本発明の超音波診断装置の別の具体例のコンソール１１０と、画像表示部１０３の斜視図を示す。図１８の構成では、術者は通常の超音波画像（本発明の合成処理を行っていない画像）１０３を参照しながら、特定のＲＯＩ（関心領域）１１０１を設定するための用手的な操作部１１０３（例えばマウス）が用意されている。これにより、術者は特定のＲＯＩ１１０１についてのみ、本発明の各実施形態の処理を施した画像１１０２を生成することができる。また、本発明の処理を施した画像１１０２は、画像表示部１０３の別領域に表示することも可能である。
【０１４９】
また、術者が用手的な操作部１１０３により、撮像条件のパラメータを設定することができる構成とすることも可能である。
【０１５０】
（実施形態の効果の具体例）
図１９に、第４の実施形態の図８の遅延加算処理を行う整相処理部２０４の整相信号のプロファイルを示す。図１９は、ステアリング角度θを横軸に、深さを縦軸に取った２次元マップ１２０１である。図１９中の点線１２０２は、第２の実施形態の図５の構成により、ヌル角度計算部５０１が超音波素子１０５のピッチ・周波数等から算出した送受信ビームのヌル角度を示す。実線１２０３は、第４の実施形態の図８のヌル角度検出部２０７が検出したヌル角度（極小点）のトレースであり、実際に超音波を被検体に照射した場合のヌル角度を示している。
【０１５１】
図１９を見ると明らかなように、計算により決定されるヌル角度と、実際に超音波を照射して検出したヌル角度には、ずれが生じていることがわかる。
【０１５２】
図２０は、図１９の異なる３つの深さについて、第４の実施形態の整相処理部２０４の整相信号の出力の大きさをステアリング角度θ方向に示したものである。すなわち、図２０は、異なる３つの深さについての指向性を示している。図２０を見ると明らかなように、それぞれのラインプロファイル１３０１，１３０２、１３０３上に見られる極小値の位置１３０４，１３０５、１３０６は異なっている。
【０１５３】
図１９、図２０より、第５の実施形態等のように、ヌル角度検出部２０７によりヌル検出を行うことにより、正確なヌル角度の方向にステアリング方向を設定することができるため、相関性雑音の低減効果の向上が期待できることがわかる。
【０１５４】
図２１に、第６の実施形態の図１２の適応整相処理を行う整相処理部２０４の整相信号のプロファイルを示す。図２１中の点線１２０２は、第２の実施形態の図５の構成により、ヌル角度計算部５０１が超音波素子１０５のピッチ・周波数等から算出した送受信ビームのヌル角度を示す。実線１４０３は、第６の実施形態の図１２のヌル角度検出部２０７が検出したヌル角度（極小点）のトレースであり、実際に超音波を被検体に照射した場合のヌル角度を示している。
【０１５５】
図２２は、図２１の異なる３つの深さについて、第６の実施形態の整相処理部２０４の整相信号の出力の大きさをステアリング角度θ方向に示したものである。図２２を見ると明らかなように、それぞれのラインプロファイル１５０１，１５０２、１５０３上に見られる極小値の位置１５０４，１５０５、１５０６は異なっている。
【０１５６】
図２１、図２２より、適応整相処理を行う場合についても第６の実施形態等のように、ヌル角度検出部２０７によりヌル検出を行うことにより、正確なヌル角度の方向にステアリング方向を設定することができるため、相関性雑音の低減効果の向上が期待できることがわかる。
【０１５７】
最後に図２３を用いて、本発明のヌル検出の効果を示す。図２３（ａ）、（ｂ）は、第６の実施形態の適応整相処理によりヌル角度を検出する図１２の超音波撮像装置により、異なる被検体を対象に超音波を照射したときのある深さにおけるステアリング方向θについての画像コントラスト１６０１，１６０２のプロファイルを示している。図２３（ａ）、（ｂ）において、点線１６０３、１６０４は、そえぞれ、ヌル角度検出部２０７が検出したヌル角度を示している。ヌル角度１６０３と１６０４の角度は、異なっており、ヌル角度は、被検体によって異なることがわかる。
【０１５８】
図２３（ａ），（ｂ）のラインプロファイル１６０１、１６０２を見ると、ちょうどヌル角度の方向において画像コントラストが最大になっていることがわかる。ヌル角度の画像コントラストは、他のステアリング角度と比べて5-6ｄＢ高い画像コントラストが得られている。これらのことから、第６の実施形態のようにヌル角度を検出して、左右のヌル角度の方向をステアリング方向に設定して、整相信号を合成することにより、ノイズを低減し、画像コントラストを顕著に向上させることができることが示されている。
【符号の説明】
【０１５９】
１００  被検体
１０１  超音波素子アレイ
１０２  装置本体
１０３  画像表示部
１０４  送信ビームフォーマ
１０６  超音波探触子
１０７  送受信分離回路（Ｔ／Ｒ）
１０８  受信ビームフォーマ
１０９  画像処理部
１１０  コンソール
１１１  制御部
２０４  整相処理部
２０５  メモリ部
２０６  周辺情報合成部
２０７  ヌル角度検出部
３０１  適応ビームステアリング部
５０１  ヌル角度計算部
７０１  チャンネルメモリ
７０２  フレーム加算部
１００１、１００２  つまみ部
【図１】
2014109392.tif 000028
【図２】
2014109392.tif 000029
【図３】
2014109392.tif 000030
【図４】
2014109392.tif 000031
【図５】
2014109392.tif 000032
【図６】
2014109392.tif 000033
【図７】
2014109392.tif 000034
【図８】
2014109392.tif 000035
【図９】
2014109392.tif 000036
【図１０】
2014109392.tif 000037
【図１１】
2014109392.tif 000038
【図１２】
2014109392.tif 000039
【図１３】
2014109392.tif 000040
【図１４】
2014109392.tif 000041
【図１５】
2014109392.tif 000042
【図１６】
2014109392.tif 000043
【図１７】
2014109392.tif 000044
【図１８】
2014109392.tif 000045
【図１９】
2014109392.tif 000046
【図２０】
2014109392.tif 000047
【図２１】
2014109392.tif 000048
【図２２】
2014109392.tif 000049
【図２３】
2014109392.tif 000050
【国際調査報告】
2014109392.tif 000051
2014109392.tif 000052