JP2009240667A - 超音波撮像装置及び超音波撮像方法 - Google Patents

Abstract

【課題】超音波ビームの方向を順次切り換えながら超音波の送受信を何度も繰り返すことなく超音波断層画像を得ることができる超音波撮像装置を提供する。
【解決手段】この超音波撮像装置は、複数の超音波トランスデューサからそれぞれ出力される複数の受信信号を処理する送受信部と、被検体の断面における複数のサンプリング点からの超音波エコーを表す受信信号行列を生成するサンプリング部と、被検体の断面における複数の位置からの超音波エコーが合成されるときの重み付け情報を表す変換行列を計算する変換行列計算手段と、変換行列計算手段によって計算された変換行列の逆行列を計算する逆行列計算手段と、逆行列と受信信号行列との積を計算することによって、被検体内の反射音源に関する情報を表す音源信号行列を求め、被検体内の超音波画像を表す画像信号を生成する音源情報計算手段とを具備する。
【選択図】図１

Description

本発明は、超音波を送受信して生体内臓器の診断や非破壊検査を行うための超音波撮像装置及び超音波撮像方法に関する。

一般に、超音波撮像装置においては、超音波の送受信機能を有する複数の超音波トランスデューサを含む超音波探触子（プローブ）が用いられる。このような超音波探触子から、複数の超音波の合波によって形成される超音波ビームを被検体に向けて送信すると、超音波ビームは、被検体内部における音響インピーダンスが異なる領域、即ち、組織の境界において反射される。そのようにして生じた超音波エコー信号を受信し、超音波エコー信号の強度に基づいて画像を構成することにより、被検体内部の様子を表示画面に再現することができる。

従来の超音波撮像装置においては、被検体内に形成された超音波ビームが反射して得られる超音波エコーから超音波ビーム上の反射信号情報を推定し、超音波ビームの方向を順次切り換えることによって、関心領域全体の情報を得る方法が用いられている。これを通常のモニターに表示するためには、いわゆる走査変換が行われる。

図２０は、送信ビームフォーミングを説明するための図である。電子走査と呼ばれる方法では、配列された振動子群が用いられる。振動子群に含まれている個々の振動子に高周波の電圧（駆動信号）を与えることによって、超音波が被検体内に導入される。このとき、図２０に示すように、個々の振動子を駆動する時刻を変化させることで、被検体に導入される音にずれが生じる。各振動子が十分に小さければ、個々の振動子が発生する超音波は球面状に生体内を伝搬して行くが、このずれによって、超音波の位相が揃って超音波が強め合う領域が形成される。

図２０に示すように、ある位置に焦点を形成する場合には、その点から同心円上に各振動子からの波面が位置するように送信時刻を調整すれば良い。焦点位置においては、全ての振動子から送信された超音波の位相が揃うので、超音波が最も強め合うことになる。同様に、別の位置においても、一部の振動子からの超音波は強め合い、また別の位置においては、位相の反転によって超音波が弱め合う領域も形成される。その結果、いわゆるサイドローブが形成されることになる。このように、ビーム状の超音波音場のみが存在する訳ではなく、全領域に何らかの音場が形成される。

図２１は、受信ビームフォーミングを説明するための図である。受信ビームフォーミングの仕組みも、送信ビームフォーミングと全く同様である。超音波エコーの受信時においては、反射体から超音波エコーが球面状に広がって反射して来るものと想定される。従って、受信時においても、複数の振動子から出力される受信信号に送信時と同様の遅延を与えた後にそれらの受信信号を加算すれば、反射体からの受信信号が最も強め合うことになる。

ここで、反射体からの超音波エコーが各振動子に到達するのと同じ時刻に別の音源からの超音波エコーも振動子に到達すれば、別の音源からの超音波エコーも強め合う。即ち、加算された受信信号に含まれるのは、必ずしも関心ある反射体からの超音波エコーだけではない。これがノイズとなるが、送信超音波ビームが振動子と反射体とを結ぶ直線上にのみ存在する場合には、このノイズも存在しない。しかしながら、前述のように、従来の送信方式では、このような送信超音波ビームを形成することは不可能であるから、音場に起因したノイズの混入も不可避であると言える。

図２０及び図２１においては、振動子面に直交する方向の超音波ビームを示しているが、前述のように振動子が十分に小さい場合には、送信される超音波が点音源として四方に広がり、また受信時においても四方からの超音波エコーを受信することができる。図２２に示すように、遅延量を増大させて左右（図中の上下）で非対称にすれば、超音波ビームを傾けることができる。図２２は受信時を示しているが、送信時も全く同様である。一般に、セクター走査と呼ばれる方法においては、振動子幅をかなり小さくして、このような送受信方法が用いられている。

前述のように音場は全領域に形成されているが、送受信される超音波ビームは細いので、主として得られる受信信号も、超音波ビーム上の限られた領域からの情報を表している。従って、断層像を得るためには関心領域全体に超音波ビームの焦点を順次合わせなければならないので、超音波ビームを移動させる方法が取られている。

図２３は、リニア走査の場合における超音波ビームの移動を示す図である。図２３に示すように、送受信に寄与する振動子群を切り換えることによって、超音波ビームが移動する。しかしながら、この方法によれば、超音波ビームの間隔が振動子の間隔と等しくなり、断層画像の分解能がこのビーム間隔に依存してしまうので、十分な画像分解能が得られない。そこで、より多くの超音波ビームを得るために様々な方法が考案されている。図２４は、その一例を示しており、微小角法と呼ばれるものである。微小角法によれば、同じ振動子群を用いながら、遅延量の左右対称性をなくして、僅かに方向の異なるビームを得ることができる。また、セクター走査の場合には、図２２における遅延量を変えることによってビームの向きを切り換えて行く方法が取られている。

こうして得られた超音波ビーム上の情報は、初期の超音波撮像装置においては、Ｘ−Ｙモニター上に直接描画されていた。しかしながら、データの記録等の観点から、現在では、走査変換という手法で一般のモニター上に描画される。

図２５は、リニア走査の場合の走査変換の例を示す図である。リニア走査の場合の走査変換においては、超音波ビーム座標における菱形印で示す点上のデータを、モニターのＸ−Ｙ座標における丸印で示す点（Ｘ，Ｙ）上のデータに変換する処理が行われる。図２６は、セクター走査の場合の走査変換の例を示す図である。セクター走査の場合の走査変換においては、超音波ビーム座標における菱形印で示す点上のデータを、極座標における丸印で示す点（ｒ，φ）、即ち、モニターのＸ−Ｙ座標における丸印で示す点（Ｘ，Ｙ）上のデータに変換する処理が行われる。

微小な領域を見たときに、超音波ビームのデータ数がＸ−Ｙ座標のデータ数よりも多い場合には、データの間引きが必要になり、逆に、超音波ビームのデータ数がＸ−Ｙ座標のデータ数よりも少ない場合には、不足分のデータの追加が必要になる。ここで、不用意な間引きは、エイリアシングを引き起こす。また、データの追加においては、周囲の超音波ビームのデータから補間処理によって不足分のデータを作成する方法が主に取られているが、補間処理は画像処理の１つであり、周辺の情報に基づいて未知の情報を推定しているに他ならない。

このように、従来の超音波撮像装置においては、断層画像を得るために、ビームの方向を変えて超音波の送受信を何度も繰り返す必要があった。また、走査変換処理は、必ずしも生体の反射音源情報を忠実に再現するものではなく、ある種の人為的処理によって良好な画質を得る手段であるとも言える。さらに、超音波ビームの構成は、幾何学的な関係に依存するので、設定する生体内音速によって、得られる画像が異なるという欠点も有する。

関連する従来技術として、特許文献１には、幅が広いビームを用いて、１回の超音波送信によって画像を得るシステム及び方法が開示されている。しかしながら、特許文献１には、断層画像を推定する具体的な方法は、必ずしも開示されていない。

また、特許文献２には、ＳＮ比を向上させて高画質化を図ることを目的とする超音波撮像装置が開示されている。この超音波撮像装置においては、受信信号が直交検波器によって複素信号に変換され、実数部及び虚数部のそれぞれについてフィルタリングが行われ、実数部及び虚数部に基づいて振幅が計算されて、その振幅に基づいて画像が表示される。しかしながら、受信信号を得るためには、１つの方向について複数回の超音波の送受信が実行される。
米国特許第６６８５６４５号明細書（第１頁、図５） 特開平１１−８９８４６号公報（第１頁、図１）

そこで、上記の点に鑑み、本発明は、超音波ビームの方向を順次切り換えながら超音波の送受信を何度も繰り返すことなく超音波断層画像を得ることができる超音波撮像装置及び超音波撮像方法を提供することを目的とする。

上記課題を解決するため、本発明の１つの観点に係る超音波撮像装置は、少なくとも１つの超音波トランスデューサに駆動信号を供給して被検体に超音波を送信すると共に、被検体から超音波エコーを受信した複数の超音波トランスデューサからそれぞれ出力される複数の受信信号を処理する送受信部と、送受信部によって処理された複数の受信信号をサンプリングすることにより、被検体の断面における複数のサンプリング点からの超音波エコーを表す受信信号行列を生成するサンプリング部と、被検体の断面における複数の位置からの超音波エコーが合成されるときの重み付け情報を表す変換行列を計算する変換行列計算手段と、変換行列計算手段によって計算された変換行列の逆行列を計算する逆行列計算手段と、逆行列計算手段によって計算された逆行列とサンプリング部によって生成された受信信号行列との積を計算することによって、被検体内の反射音源に関する情報を表す音源信号行列を求め、該音源信号行列に基づいて、被検体内の超音波画像を表す画像信号を生成する音源情報計算手段とを具備する。

また、本発明の１つの観点に係る超音波撮像方法は、被検体の断面における複数の位置からの超音波エコーが合成されるときの重み付け情報を表す変換行列を計算するステップ（ａ）と、ステップ（ａ）において計算された変換行列の逆行列を計算するステップ（ｂ）と、少なくとも１つの超音波トランスデューサに駆動信号を供給して被検体に超音波を送信すると共に、被検体から超音波エコーを受信した複数の超音波トランスデューサからそれぞれ出力される複数の受信信号を処理するステップ（ｃ）と、ステップ（ｃ）において処理された複数の受信信号をサンプリングすることにより、被検体の断面における複数のサンプリング点からの超音波エコーを表す受信信号行列を生成するステップ（ｄ）と、ステップ（ｂ）において計算された逆行列とステップ（ｄ）において生成された受信信号行列との積を計算することによって、被検体内の反射音源に関する情報を表す音源信号行列を求め、該音源信号行列に基づいて、被検体内の超音波画像を表す画像信号を生成するステップ（ｅ）とを具備する。

本発明によれば、被検体の断面における複数の位置からの超音波エコーが合成されるときの重み付け情報を表す変換行列の逆行列と、被検体の断面における複数のサンプリング点からの超音波エコーを表す受信信号行列との積を計算することによって、被検体内の反射音源に関する情報を表す音源信号行列を求め、それに基づいて被検体内の超音波画像を表す画像信号を生成するので、超音波ビームの方向を順次切り換えながら超音波の送受信を何度も繰り返すことなく超音波断層画像を得ることができる。

以下、本発明の実施形態について、図面を参照しながら詳しく説明する。なお、同一の構成要素には、同一の参照番号を付して説明を省略する。
図１は、本発明の第１の実施形態に係る超音波撮像装置の構成を示すブロック図である。この超音波撮像装置は、超音波探触子１０と、駆動信号生成部２０と、受信信号処理部３０と、サンプリング部４１と、データ格納部４２と、画像信号生成部５０と、Ｄ／Ａ変換器６１と、表示部６２と、操作卓７１と、制御部７２と、格納部７３とを有している。ここで、駆動信号生成部２０及び受信信号処理部３０は、送受信部を構成する。

超音波探触子１０は、１次元のトランスデューサアレイを構成する複数の超音波トランスデューサ１０ａを備えている。各々の超音波トランスデューサ１０ａは、印加される駆動信号に基づいて超音波を送信すると共に、伝搬する超音波エコーを受信して受信信号を出力する。

各々の超音波トランスデューサは、例えば、ＰＺＴ（チタン酸ジルコン酸鉛：Pb(lead) zirconate titanate）に代表される圧電セラミックや、ＰＶＤＦ（ポリフッ化ビニリデン：polyvinylidene difluoride）に代表される高分子圧電素子等の圧電性を有する材料（圧電体）の両端に電極を形成した振動子によって構成される。そのような振動子の電極に、パルス状又は連続波の電気信号を印加すると、圧電体が伸縮する。また、振動子は、伝搬する超音波を受信することによって伸縮し、電気信号を発生する。その電気信号は、超音波の受信信号として出力される。

駆動信号生成部１４は、例えば、複数の超音波トランスデューサ１０ａに対応する複数のパルサによって構成されている。駆動信号生成部２０は、１つ又は複数の超音波トランスデューサを駆動するための駆動信号を超音波探触子１０に供給するか、あるいは、制御部７２によって選択された送信遅延パターンに基づいて、複数の超音波トランスデューサ１０ａから送信される超音波が所定の波面又はビームを形成するように複数の駆動信号の遅延量を調節して超音波探触子１０に供給する。

受信信号処理部３０は、複数の超音波トランスデューサ１０ａに対応して、複数の増幅器（プリアンプ）３１と、複数のＡＡ（アンチ・エイリアシング）フィルタ３２とを含んでいる。超音波トランスデューサ１０ａから出力される受信信号は、増幅器３１において増幅され、増幅器３１から出力される受信信号は、ＡＡフィルタ４１によって、サンプリングにおけるエイリアシングを防止するためのフィルタ処理が施される。また、図には示されていないが、受信信号は、検波処理によって、その包絡線情報に変換されている。

サンプリング部４１は、受信信号処理部３０から出力される複数の受信信号をサンプリングしてＡ／Ｄ変換することにより、被検体の断面におけるｍ×ｎ個のサンプリング点からの超音波エコーを表す受信信号行列（受信データ）を生成する。例えば、サンプリング部４１は、被検体の断面における深さ方向を第１の座標軸とし、複数の超音波トランスデューサ１０ａの配列方向を第２の座標軸として、受信信号処理部３０によって処理された複数の受信信号を２次元的にサンプリングする。

このときのサンプリング周波数は、サンプリング定理で通常規定される周波数よりも低い。そのため、先に説明したように、エイリアシングを防止するためのＡＡフィルター３２が必要となる。ＡＡフィルター３２におけるフィルタ設定は、サンプリング周波数に連動し、演算すべき断面における計算領域（ｍ×ｎ個のサンプリング点）の指定に依存する。サンプリング部４１によってサンプリングされた受信データは、ハードディスク又はメモリ等によって構成されるデータ格納部４２に一旦格納される。

画像信号生成部５０は、断面設定部５１と、計算領域指定部５２と、変換行列計算部５３と、逆行列計算部５４と、音源情報計算部５５とを含んでいる。断面設定部５１は、制御部７２の制御の下で、超音波画像を表示すべき被検体の断面を設定する。計算領域指定部５２は、断面設定部５１によって設定された断面において、計算領域（ｍ×ｎ個のサンプリング点）を指定する。

計算領域が指定されると、変換行列計算部５３は、被検体の断面における複数の位置からの超音波エコーが合成されるときの重み付け情報を表す変換行列を計算する。例えば、変換行列計算部５３は、被検体の断面における深さ方向を第１の座標軸とし、被検体の断面内で第１の座標軸に直交する方向を第２の座標軸として、被検体の断面において２次元的に配置された複数の位置からの超音波エコーが合成されるときの重み付け情報を表す変換行列を計算する。逆行列計算部５４は、変換行列計算部５３によって計算された変換行列の逆行列を計算する。

音源情報計算部５５は、逆行列計算部５４によって計算された変換行列の逆行列と、サンプリング部４１によって生成されてデータ格納部４２に格納されている受信信号行列との積を計算することによって、被検体内の反射音源に関する情報を表す音源信号行列を求め、該音源信号行列に基づいて、被検体内の超音波画像を表す画像信号を生成する。

Ｄ／Ａ変換器６１は、画像信号生成部５０から出力されるディジタルの画像信号を、アナログの画像信号に変換する。表示部６２は、例えば、ＣＲＴやＬＣＤ等のディスプレイ装置を含んでおり、アナログの画像信号に基づいて診断画像を表示する。

操作卓７１を用いたオペレータの操作に従って、制御部７２が、駆動信号生成部２０、受信信号処理部３０、画像信号生成部５０等を制御する。画像信号生成部５０及び制御部７２は、ＣＰＵとソフトウェア（プログラム）によって実現することができる。ソフトウェア（プログラム）は、格納部７３に格納される。格納部７３における記録媒体としては、内蔵のハードディスクの他に、フレキシブルディスク、ＭＯ、ＭＴ、ＲＡＭ、ＣＤ−ＲＯＭ、又は、ＤＶＤ−ＲＯＭ等を用いることができる。なお、画像信号生成部５０を、ディジタル回路によって構成しても良い。

次に、図１に示す超音波撮像装置の動作原理について詳しく説明する。
本発明は、従来技術におけるように受信信号を走査変換することなく、全く異なる演算方法によって断層情報を推定するものである。被検体に送信された超音波は、被検体内において、固有音響インピーダンスが異なるいわゆる音響境界面によって反射される。音響理論によれば、反射源が鏡面境界である場合には、音が入射角度と同じ反射角度をもって反射される。即ち、音が送信方向に反射されるためには、超音波ビームと垂直に反射面が存在しなければならない。

しかしながら、例えば、生体における反射源は超音波の波長よりも十分に小さいので、音は散乱を起こして全方位的に反射される。従来技術におけるように、細い超音波ビームを送信しながらも、そのビーム方向からの超音波エコーが受信波として得られて断層像を構築できるという事実は、このことの紛れも無い証拠である。そうだとすると、反射源からの反射信号は、配列された振動子のいずれにも略同等に帰って来るはずである。本発明は、１つの振動子が受ける受信信号は、被検体の断面内の全ての場所からの反射信号の総和であるという考えに基づくものである。

即ち、本発明は、受信信号が送信超音波ビーム上にあることを想定せず、超音波音場に関わりなく断層情報を得るものである。被検体内の反射源からの超音波エコーは、点音源からの超音波として振動子に戻って来るものと想定する。１つの振動子が受信する超音波エコーは、求める断面内の全領域から反射したものであり、振動子から同一距離にある反射源が、特定時刻における受信信号に寄与していると考える。そうすると、複数の振動子から得られる受信信号に基づいて、断面内の反射音源情報を変数とした連立１次方程式を立てることができ、この方程式を解くことによって断層画像が得られる。この方法によれば、十分に広い音場を与えることにより、１回の送受信で断層像を得ることが可能となる。また、生体内音速は均一であるという仮定は残るものの、音源情報の推定に音速の設定は不要である。

図２は、配列された振動子群と被検体内における超音波反射源との位置関係を示す図である。ここでは、配列された振動子群に含まれている個々の振動子に対応して、被検体内に有限個の超音波反射源が存在すると考える。図２に示すように、ｊ番目の振動子（以下、「振動子ｊ」ともいう）におけるｉ番目の深さ（以下、「深さｉ」ともいう）の超音波反射源からの超音波エコーに基づく受信信号をｒ_ｉｊとする。

配列された振動子の一部あるいは全部を用いて、パルス状の超音波が被検体内に送信されるものとする。パルス状の超音波が送信されてから次のパルス状の超音波が送信されるまでの間において超音波エコーが受信されることは、従来法と同じである。ある時刻において、この送信音場は、断層面内で波面を形成する。断層面内にある超音波反射源にこの波面が到達すると、それが音源となって振動子に向かう反射波が発生する。この断層面をｘｙ平面で表すことにすると、断層面内のある位置における音源信号をｓ_ｙｘとすることができる。

図３は、断層面の座標系と受信信号の座標系とを同じ座標系で示す図である。ｘ，ｙは離散的な値をとるものとして、図３においては、行列表現に合わせて直感的に理解しやすいように、深さｙが行に相当し、横方向の距離ｘが列に相当するようにした。

音源信号ｓ_ｙｘは、点音源性を有し、全方位に反射音が広がるものと考えれば、ｊ番目の振動子には、断面上の全ての音源信号が到達するものと考えられる。受信信号ｒ_ｉｊは、特定の時刻における受信信号であり、音場の音速ｃが一定であれば、時刻２ｉ／ｃに振動子ｊに到達したものになる。これは、特定の深さからの反射信号であり、この時刻に振動子ｊに到達する信号は、音速に関わらず、振動子ｊから深さｉと等距離にある音源が寄与していることになる。そのため、音源信号ｓ_ｙｘには、それぞれａ_ｙｘの重み付けがされるものとし、振動子ｊからの距離が受信信号ｒ_ｉｊと同一でない音源についてはａ_ｙｘ＝０とすれば良い。ただし、ａ_ｙｘはｉ，ｊに依存する重み付け情報であるから、これを改めてａ_{ｉｊ−ｙｘ}と記載する。

ここで、音源信号として、ｍ行ｎ列のマトリックスを考える。これは、求める断面の空間解像度を、深さ方向ｍ、横方向ｎと設定することと等価である。図３においては、断層面の座標と振動子による受信信号の座標とを同じにしている。受信信号ｒ_ｉｊは、次式（１）で表されるように、図３に示す円弧上にあるｓ_ｙｘの和となる。

ただし、ある時刻における受信信号には、等距離に位置する音源のみが寄与するから以下の条件が得られる。
（ｊ−ｘ）^２＋ｙ^２＝ｉ^２のときに、ａ_{ｉｊ−ｙｘ}＝１
（ｊ−ｘ）^２＋ｙ^２≠ｉ^２のときに、ａ_{ｉｊ−ｙｘ}＝０
・・・（２）

式（１）を展開して行列で表現すると、次式（３）となる。

振動子及び受信信号の位置を音源の位置と一致させると、ｉ_ｍａｘ＝ｍ、ｊ_ｍａｘ＝ｎであるから、式（３）から全ての受信信号についての行列を求めると、次式（４）となる。

それぞれの行列を、Ｒ，Ａ，Ｓで書くことにすれば、次式（５）となる。
ＡＳ＝Ｒ ・・・（５）
即ち、音源信号行列Ｓは、式（４）の連立１次方程式を解くことによって与えられる。連立１次方程式は、掃き出し法あるいはガウスの消去法と呼ばれる方法で解くことができる。また、変換行列Ａは、サイズがｍ×ｎであるから正方行列となり、逆行列Ａ⁻¹を求めることができる。音源信号行列Ｓは、次式（６）で示すように、変換行列の逆行列Ａ⁻¹と受信信号行列Ｒとの積で与えられる。
Ｓ＝Ａ⁻¹Ｒ ・・・（６）

変換行列Ａは、式（２）で示されるように、幾何学的な位置関係のみで与えられるので、信号を受信する前に事前に算出しておくことができる。従って、逆行列Ａ⁻¹も事前に演算しておくことができるから、実際に受信期間中に行わなければならない処理は、式（６）の行列の掛け算だけである。しかも、音源信号行列Ｓの内で求めたい要素のみを演算することが可能である。

図４は、断層面の座標系と受信信号の座標系とを別の座標系で示す図である。図３においては、求める断層面の座標系と受信信号の座標系とを同一にしたが、これは必ずしも必要ではない。ｍ×ｎ個の変数を有する音源信号行列Ｓに対して、ｍ×ｎ個の連立方程式が立てられれば良いのであって、例えば、図４に示すような座標設定であってもかまわない。この場合には、変換行列Ａの条件として、次式（７）を用いれば良い。
（ｘ_０＋ｄｊ−ｘ）^２＋（ｙ−ｙ_０）^２＝（ｔｉ）^２のときに、ａ_{ｉｊ−ｙｘ}＝１
（ｘ_０＋ｄｊ−ｘ）^２＋（ｙ−ｙ_０）^２≠（ｔｉ）^２のときに、ａ_{ｉｊ−ｙｘ}＝０
・・・（７）

即ち、求める断面座標における配列振動子の一端の位置（ｘ_０，ｙ_０）と、振動子の間隔ｄとが与えられれば、変換行列Ａを算出することができる。式（７）において、ｔはサンプリング間隔によって決まる生体内距離を意味するが、これは、ｘ，ｙを生体内距離に換算するために必要なだけであって、行列の算出においてはｔ＝１としても差し支えない。

図５は、コンベックス探触子を用いる場合の変換行列の算出方法を説明するための図である。コンベックス探触子のように、複数の振動子が凸型に配置されているような場合には、断面座標のｙ軸に直交する振動子の初期座標（ｘ_０，ｙ_０）が与えられ、振動子配列の仮想中心距離Ｄが与えられれば、次式（８）に従って変換行列Ａを求めることができる。
（ｘ_０−Ｄｓｉｎθ−ｘ）^２＋（ｙ−ｙ_０＋Ｄ−Ｄｃｏｓθ）^２＝（ｔｉ）^２のときに
ａ_{ｉｊ−ｙｘ}＝１
（ｘ_０−Ｄｓｉｎθ−ｘ）^２＋（ｙ−ｙ_０＋Ｄ−Ｄｃｏｓθ）^２≠（ｔｉ）^２のときに
ａ_{ｉｊ−ｙｘ}＝０
・・・（８）
ここで、振動子間隔ｄについて、次式（９）が成立する。
θ＝ｉ・ｓｉｎ^−１（ｄ／Ｄ） ・・・（９）
ただし、このときのｊは、（ｘ_０，ｙ_０）に位置する振動子からの配置順である。

以上に示したように、本発明によれば、全ての振動子による１回の送受信によって得られる受信信号に基づいて、断面内の全ての音源情報を得ることが可能である。即ち、送信される超音波音場が連続した波面を形成し、その波面が被検体内の関心領域の全てに行き届けば、１回の送受信で断層像を得ることができる。

次に、本発明の第２の実施形態について説明する。第２の実施形態においては、音場を考慮したマトリックス演算の重み付けが行われる。
一般に、均質な媒質に球面波音場が形成されていない限り、音場強度は空間内で一定にはならない。また、音の伝搬に従って音響強度は低下するので、反射位置が深いほど反射音源強度も弱くなる。これらを考慮して、反射音源が送信音場の強度ｗによって重み付けされるものと考え、ａ_{ｉｊ−ｙｘ}の替わりにｗ_{ｉｊ−ｙｘ}を用いることができる。従って、式（２）は、次式（１０）で置き換えることができる。
（ｊ−ｘ）^２＋ｙ^２＝ｉ^２のときに、ａ_{ｉｊ−ｙｘ}＝ｗ_{ｉｊ−ｙｘ}
（ｊ−ｘ）^２＋ｙ^２≠ｉ^２のときに、ａ_{ｉｊ−ｙｘ}＝０
・・・（１０）

図６は、本発明の第２の実施形態に係る超音波撮像装置の構成を示すブロック図である。第２の実施形態においては、図１に示す第１の実施形態に対して、画像信号生成部５０ａにおいて、音場推定部５６及び重み付け処理部５７が追加されている。音場推定部５６は、制御部７２によって設定された条件に基づいて、音場を推定する。また、重み付け処理部５７は、音場推定部５６によって行われた音場推定に基づいて、変換行列計算部５３によって計算された変換行列に含まれている各要素に重み付けを与える。逆行列計算部５４は、重み付け処理部５７によって重み付け処理された変換行列の逆行列を計算する。

次に、本発明の第３の実施形態について説明する。第３の実施形態においては、複素領域におけるマトリックス演算が行われる。
図７は、本発明の第３の実施形態に係る超音波撮像装置の構成を示すブロック図である。図７に示すように、受信信号処理部３０ａにおいて、複数の超音波トランスデューサ１０ａに対応して、複数の増幅器（プリアンプ）３１と、それぞれの受信信号をベースバンド周波数の信号に変換すると共に、実部と虚部とを検出する直交検波器３３と、それぞれの受信信号に対応して実部用のＡＡフィルタ３２ａ及び虚部用のＡＡフィルタ３２ｂとが設けられている。

図８は、直交検波器の構成例を示す図である。この直交検波器は、掛算器３４ａ及び３４ｂと、低域通過フィルタ３５ａ及び３５ｂと、互いに位相が直交する第１の信号（ｃｏｓωｔ）及び第２の信号（ｓｉｎωｔ）を生成する局部発振器３６とを含んでいる。掛算器３４ａは、受信信号と第１の信号（ｃｏｓωｔ）とを掛け合わせることによって受信信号の実部（Ｉ軸成分）を求め、また、受信信号と第２の信号（ｓｉｎωｔ）とを掛け合わせることによって受信信号の虚部（Ｑ軸成分）を求める。その後の受信信号の処理は、実部と虚部とについて同様に行われる。低域通過フィルタ３５ａ及び３５ｂは、受信信号の実部及び虚部から、不要な広域成分をそれぞれ除去する。

再び図７を参照すると、この超音波撮像装置においては、受信信号処理部３０ａから出力される複数の受信信号の実部をサンプリングして実部の受信データを生成するサンプリング部４１ａと、受信信号処理部３０ａから出力される複数の受信信号の虚部をサンプリングして虚部の受信データを生成するサンプリング部４１ｂと、実部の受信データを格納するデータ格納部４２ａと、虚部の受信データを格納するデータ格納部４２ｂとが設けられている。

また、画像信号生成部５０ｂにおいて、実部の受信データに基づいて音源信号の実部を計算する音源情報計算部５５ａと、虚部の受信データに基づいて音源信号の虚部を計算する音源情報計算部５５ｂと、音源信号の実部及び虚部に基づいて音源信号の振幅を求める振幅演算部５８とが設けられている。

音源情報計算部５５ａは、逆行列計算部４３によって計算された逆行列と、サンプリング部４１ａによって生成されてデータ格納部４２ａに格納されている受信信号行列との積を計算することによって、被検体内の反射音源に関する情報を表す音源信号行列（実部）を求める。また、音源情報計算部５５ｂは、逆行列計算部４３によって計算された逆行列と、サンプリング部４１ｂによって生成されてデータ格納部４２ｂに格納されている受信信号行列との積を計算することによって、被検体内の反射音源に関する情報を表す音源信号行列（虚部）を求める。

振幅演算部５８は、音源信号行列（実部）に含まれている音源信号の実部の自乗と、音源信号行列（虚部）に含まれている音源信号の虚部の自乗とを加算して、その平方根を計算することにより音源信号の振幅を求め、該振幅に基づいて、被検体内の超音波画像を表す画像信号を生成する。

第３の実施形態によれば、受信信号を直交検波することによって受信信号の位相情報を保持し、受信信号を複素情報として演算を行うことによって、ＳＮ比を向上させて高画質化を図ることができる。なお、図７においては、図６に示す第２の実施形態と同様に、音場推定部５６及び重み付け処理部５７が設けられているが、図１に示す第１の実施形態と同様に、それらを設けなくても良い。

次に、本発明の第４の実施形態について説明する。第４の実施形態においては、素子数よりも多い超音波分解能を得る方法が用いられる。超音波撮像装置の基本的な構成は、第１〜第３の実施形態のいずれかと同様である。

音源情報の算出に必要なことは、必要数の連立１次方程式を立てることに帰結する。即ち、求める断面座標データ数（ｍ×ｎ個）と同数の受信信号が得られれば良く、変換行列Ａの算出条件は、受信信号の座標系、即ち、探触子の形状と位置とによって定まる。この考え方によれば、振動子数が断層面の方位分解能（ｘ方向）に満たないような場合であっても、受信信号のサンプル間隔を短くすることによってｍ×ｎ個のサンプルデータを確保すれば良い。

図９は、本発明の第４の実施形態における配列された振動子群と受信信号のサンプル間隔との関係を示す図である。図９に示すのは、振動子数が方位分解能の半分の場合であり、距離分解能（ｙ方向）の倍の受信サンプルデータを得れば良い。ただし、標本化定理によれば、受信信号の最大周波数に対して２倍以上の周波数でサンプリングしても、実効的に独立な情報は得られない。従って、逆行列Ａ⁻¹は存在しなくなる。即ち、標本化定理に従うサンプルデータ数の上限が求める断層像の分解能の上限になる。

次に、本発明の第５の実施形態について説明する。第５の実施形態においては、部分送信の重ね合わせによる断層像構築が行われる。超音波撮像装置の基本的な構成は、第１〜第３の実施形態のいずれかと同様である。

例えば、図１に示す音源情報計算部５５が、被検体の断面を複数の領域に分けたときの各領域について、逆行列計算部５４によって計算された逆行列と、サンプリング部４１によって生成された受信信号行列との積を計算する。

図１０は、被検体の断面を小さい領域に分けて断層像を作成する例を示す図である。本発明による音源情報の計算は、連立１次方程式を解くことに帰結するが、断面座標データが多いほど計算は複雑になる。ｍ×ｎ個の方程式が成立することになる訳であるから、この数が減るほど演算は簡単になる。即ち、被検体の断面全体を一度に演算するよりも、複数の比較的小さい領域に分けて演算する方が、計算時間が軽減される。このとき、分けられた領域の幾何学的条件を一定に保てば、連立方程式における変換行列Ａは同じであるから、計算の事前準備が節約されるだけでなく、演算アルゴリズムも簡略化することができる。

次に、本発明の第６の実施形態について説明する。第６の実施形態においては、送信フォーカスを形成して関心領域のみの演算を行うことにより、部分送信の重ね合わせによる断層像構築が行われ、高速な撮像とＳＮ比の向上が図られる。超音波撮像装置の基本的な構成は、第１〜第３の実施形態のいずれかと同様である。

例えば、図１に示す駆動信号生成部２０が、被検体内において順次選択された断面領域に送信フォーカスを形成するように、複数の超音波トランスデューサ１０ａに複数の駆動信号をそれぞれ供給し、音源情報計算部５５が、順次選択された断面領域について、逆行列計算部５４によって計算された逆行列とサンプリング部４１によって生成された受信信号行列との積を計算する。

本発明によれば１回の送受信で断層像を取得することが可能であるが、セクター走査のように狭い超音波開口から扇状に広がる全空間に対して１回の送信で十分な音圧を発生することは困難である。従って、セクター用探触子においては、従来技術と同様の位相制御によって超音波ビーム方向を制御することで、生体内の音圧を部分的に高める必要がある。この場合においても、全空間で定義された逆行列Ａ⁻¹から断層像を求めることは可能であるが、音圧が低い領域ではＳＮ比が低下して良好な画像は得られない。

そこで、図１１に示すように、位相制御されて十分なＳＮ比が得られる領域のみの音源情報Ｓを算出する方が、計算時間も短く効率的である。そのために、送信方向を順次切り換えながら、被検体内において順次選択された断面領域に送信フォーカスを形成し、算出する音源情報の領域も送信方向に合わせて変更する。これを繰り返すことによって、関心領域全体の音源情報を得ることが可能である。本実施形態においては、超音波の送信については従来と同様な方法をとるのではあるが、画像構築において超音波ビームの太さや焦点などの影響によって分解能が左右される訳ではなく、良好なＳＮ比を得るために十分な被検体内の音圧を得ることが目的である。

次に、本発明の第７の実施形態について説明する。第７の実施形態においては、多重反射を考慮して変換行列及び逆行列を作成することにより、点音源に関する条件が拡張される。超音波撮像装置の基本的な構成は、第１〜第３の実施形態のいずれかと同様である。

振動子の表面においては、受信した超音波が大きく反射する。そのため、振動子から近距離に位置する強い音源から振動子に到達した超音波は、振動子の表面で反射して再び被検体内に導かれる。振動子の表面において反射した超音波が強い音源から発した超音波であるとすれば、そこには強い反射源があることが予想される。そのため、この超音波は、再びその反射源で反射し、振動子に向かうことが予想される。

図１２は、多重反射した超音波が振動子に到達する状態を示す図である。強い反射源によって多重反射した超音波は、図１２に示されるように、元の受信信号の伝播距離の３倍の距離を伝播して受信されることになる。これは、即ち、深さｉにおける受信信号には深さｉ／３にある音源情報が加算されている可能性を示すものである。

そこで、変換行列Ａにおける値として、式（１１）に示すようにｉ／３の深さからの受信信号を音源情報として加えることが有効である。例えば、図１に示す変換行列計算部５３が、ある深さｉからの受信信号に対して、その１／３の深さ（ｉ／３）にある音源情報が加算されているとして変換行列を計算する。この音源情報がどの程度寄与するかは一概に定めることはできないが、ｉ／３の深さからの受信信号の強度ｒ_{（ｉ／３）ｊ}に適当な係数ｂを掛けることが妥当である。
（ｊ−ｘ）^２＋ｙ^２＝ｉ^２のときに、ａ_{ｉｊ−ｙｘ}＝ｗ_{ｉｊ−ｙｘ}
（ｊ−ｘ）^２＋ｙ^２＝（ｉ／３）^２のときに、ａ_{ｉｊ−ｙｘ}＝ｂｒ_{（ｉ／３）ｊ−ｙｘ}
・・・（１１）

次に、本発明の第８の実施形態について説明する。第８の実施形態においては、拡大率に応じてサンプル周期を変え、限定された部分のみを演算することによって拡大表示が行われる。超音波撮像装置の基本的な構成は、第１〜第３の実施形態のいずれかと同様である。

例えば、図１に示すサンプリング部４１が、受信信号処理部３０によって処理された複数の受信信号のサンプリング周波数を高めることにより、音源情報計算部５５が、拡大された超音波画像を表す画像信号を生成する。

本発明において、連立１次方程式の数は、ｍ×ｎで与えられる。即ち、モニター上の画像表示解像度が一定であれば、断層面の被検体距離には依存しない。これは、拡大縮小表示をしても演算難易度が変わらないことを意味する。必要なことは、拡大縮小を行ってもｍ×ｎの数を保つことである。

図１３は、拡大表示の例を示す図である。一般に、図１３（ａ）に示す領域の一部を、図１３（ｂ）に示すように拡大表示を行う場合には、受信に関与する振動子の数は減少し、対象となる受信信号の時間範囲も減少する。即ち、そのままではデータ数が減少するが、サンプリング周波数を高めることで同じデータ数を確保することが可能である。

サンプリング周波数をどの程度高めればデータ数を確保できるかは、拡大の仕方にも依存するので、一概に決めることはできない。求める音源領域に寄与する受信信号の範囲において、ｍ×ｎ個のデータ数を確保するにはサンプリング周波数を何倍にする必要があるかという観点で決定される。ただし、サンプリング定理によって定まる最大サンプリング周波数、即ち、ナイキスト周波数以上にすることは意味がない。ナイキスト周波数以上でサンプルしても、独立した方程式が得られないので、結局、連立方程式の解を得ることができないからである。逆に言えば、ナイキスト周波数でサンプルしたときにｍ×ｎ個のデータ確保が可能な範囲が、最大拡大率である。

次に、本発明の第９の実施形態について説明する。第９の実施形態においては、任意の選択位置のみを演算することによって、任意方向の画像情報の時間的変化が観測される。超音波撮像装置の基本的な構成は、第１〜第３の実施形態のいずれかと同様である。

例えば、図１に示す音源情報計算部５５が、選択された断面領域について、逆行列計算部５３によって計算された逆行列とサンプリング部４１によって生成された受信信号行列との積を繰返し計算することにより、選択された断面領域における画像情報の時間的変化を表す画像信号を生成する。

図１４は、断面における任意方向の画像情報を選択する例を示す図である。行列による連立方程式の解法では、未知数Ｓをすべて同時に求める必要はない。所望するＳ_ｙｘのみを計算することができる。例えば、断面における任意方向の画像情報の時間的変化を求めることは、臨床上有効である。これは、任意方向Ｍモードと呼ばれるが、このとき、図１４に示す点のみを選択的に計算すれば良く、演算数が少なくなる。また、任意の音源情報を選択的に計算することが可能であるから、求める場所は任意であり、直線上の場所である必要もない。たとえば、図１４に示されるような心臓の場合に、心筋に沿った範囲を計算対象として心筋の変化を追従することも可能になる。

次に、本発明の第１０の実施形態について説明する。第１０の実施形態においては、振動子の大きさを考慮して逆行列を算出することにより、点音源に拡がりを持たせた重み付けが行われる。超音波撮像装置の基本的な構成は、第１〜第３の実施形態のいずれかと同様である。

例えば、図１に示す変換行列計算部５３が、ある深さからの受信信号に対して、各超音波トランスデューサ（振動子）１０ａの幅に対応する時間差を伴って受信される受信信号が加算されているとして変換行列を計算する。

図１５は、複数の振動子から送信された超音波が被検体内において反射されて生じる反射波を示す図である。振動子はある程度の大きさを有するので、同じ時刻において１つの振動子によって受信される信号は、図１５に示すように、振動子の大きさに相当する空間からの反射波の総和になる。深さｉに相当する時刻において、振動子の中心部分Ｃが受信する信号は、振動子中心から半径ｉの同心円上の点からの反射波である。

振動子の端部Ｌ又はＲにおいては、端部から半径ｉの同心円上の点を考慮しなければならず、これは、左の端部Ｌと右の端部Ｒとにおいては、相当な角度に広がる。全体として、円弧Ｌ_ｉ及び円弧Ｒ_ｉ、そして、場所によっては円弧Ｃ_ｉによって囲まれた領域からの反射波が同時刻に受信されることになる。変換行列Ａの算出において、式（１２）に示すように、振動子幅を考慮することにより、さらに正確に音源情報を推定することができる。
（ｊ＋Ｈ−ｘ）^２＋ｙ^２＝ｉ^２ 〜 （ｊ−Ｈ−ｘ）^２＋ｙ^２＝ｉ^２のときに、
ａ_{ｉｊ−ｙｘ}＝ｗ_{ｉｊ−ｙｘ} ・・・（１２）
ここで、Ｈは、振動子の幅を表している。

次に、本発明の第１１の実施形態について説明する。第１１の実施形態においては、複数の振動子に単一の駆動信号を供給することによって、送信回路（パルサ等）の規模が低減される。その他の構成は、第１〜第３の実施形態のいずれかと同様である。

図１６は、本発明の第１１の実施形態に係る超音波撮像装置の送受信部の一部を示す図である。本発明においては、送信される超音波ビームを位相制御することは、必ずしも必要ではない。従って、駆動信号は、複数の振動子について同位相及び同波形でもかまわない。ただし、受信においては、複数の振動子毎の受信信号を得る必要がある。

そこで、図１６に示すように、駆動信号生成部において共通の送信回路（パルサ等）を用いて、複数の振動子に共通の駆動信号を供給し、受信回路（プリアンプ等）のみを複数の振動子毎に設けるようにしても良い。一般に、駆動信号は十分に大きな電圧であるのに対して、受信信号は極めて微弱な電圧である。従って、受信系回路において、複数の振動子と基準電位（本実施形態においては接地電位）とを、逆方向に並列接続されたダイオードＤ１１及びＤ１２等を介して接続することによって、大きな電圧はダイオードＤ１１及びＤ１２等を通過し、微弱な電圧の受信信号は受信回路に入力されるようにすることができる。本実施形態によれば、送信回路の数を減らすことが可能であり、送信位相制御回路も不要になるので、大幅な回路規模の縮小を図ることができる。

次に、本発明の第１２の実施形態について説明する。第１２の実施形態においては、複数の振動子から送信される超音波を合成することによって球面波が形成される。その他の構成は、第１〜第３の実施形態のいずれかと同様である。

図１７は、複数の振動子から送信される超音波の合成を示す図である。本発明によれば、第１１の実施形態において説明したように、複数の振動子について駆動信号を共通化することが可能である。駆動信号が共通化された送信においては、配列振動子であっても、あたかも大きな平板振動子であるかのような挙動をし、図１７に示すように、送信音場として平面波が形成される。リニア走査やコンベックス走査においてはこれで問題ないが、振動子直下の領域以外の領域が関心領域として設定されるセクター走査においては、波面波を形成することができないことは問題となる。

そこで、図１８に示すように、送信音場が点音源からの音場となるように、中央の振動子に対して両端の振動子については駆動信号に遅延が与えられる。このとき、振動子の開口が十分に狭いものとすれば、振動子下の波面は無視して両脇に広がる波面のみにおいて関心領域を考えれば、駆動信号に与えられる遅延量は、振動子の配列位置で決まる固定量として差し支えない。

図１９は、本発明の第１２の実施形態に係る超音波撮像装置の送受信部の一部を示す図である。図１９に示すように、複数の振動子と送信回路とが、それぞれの遅延回路を介して接続される。これにより、複数の振動子にそれぞれ供給される駆動信号に固定の遅延量が与えられ、複数の振動子から送信される超音波が球面波を形成する。その他の点に関しては、図１６に示す第１１の実施形態と同様である。

以上説明した実施形態によれば、１回の送受信でも断層像を得ることができる。また、画像の生成方法が音速情報に依存しないので、音速設定によって画像が変わることがない。走査変換においては、超音波探触子の形状毎に変換方法及び補間方法を設定しなければならないが、これらの実施形態によれば、あらかじめ逆行列を求めておくことにより、全く同じ計算で断層像が得られる。送信時に超音波をフォーカスさせる必要がないので、個々の振動子に個別の駆動信号を与えなくても良く、送信回路を簡略化することができる。また、原理的には、３Ｄボリューム情報も算出可能である。さらに、複素領域で演算を行う場合には、原理上２回の送受信でもドプラ断層像が得られる。

本発明は、超音波を送受信して生体内臓器の診断や非破壊検査を行うための超音波撮像装置において利用することが可能である。

本発明の第１の実施形態に係る超音波撮像装置の構成を示すブロック図である。 配列された振動子群と被検体内における超音波反射源との位置関係を示す図である。 断層面の座標系と受信信号の座標系とを同じ座標系で示す図である。 断層面の座標系と受信信号の座標系とを別の座標系で示す図である。 コンベックス探触子を用いる場合の変換行列の算出方法を説明するための図である。 本発明の第２の実施形態に係る超音波撮像装置の構成を示すブロック図である。 本発明の第３の実施形態に係る超音波撮像装置の構成を示すブロック図である。 直交検波器の構成例を示す図である。 本発明の第４の実施形態における配列された振動子群と受信信号のサンプル間隔との関係を示す図である。 被検体の断面を小さい領域に分けて断層像を作成する例を示す図である。 順次選択された断面領域に送信フォーカスを形成して音源情報を計算する例を示す図である。 多重反射した超音波が振動子に到達する状態を示す図である。 拡大表示の例を示す図である。 断面における任意方向の画像情報を選択する例を示す図である。 複数の振動子から送信された超音波が被検体内において反射されて生じる反射波を示す図である。 本発明の第１１の実施形態に係る超音波撮像装置の送受信部の一部を示す図である。 複数の振動子から送信される超音波の合成を示す図である。 両端の振動子について駆動信号に遅延を与える例を示す図である。 本発明の第１２の実施形態に係る超音波撮像装置の送受信部の一部を示す図である。 送信ビームフォーミングを説明するための図である。 受信ビームフォーミングを説明するための図である。 超音波ビームを傾ける例を示す図である。 リニア走査の場合における超音波ビームの移動を示す図である。 微小角法を説明するための図である。 リニア走査の場合の走査変換の例を示す図である。 セクター走査の場合の走査変換の例を示す図である。

符号の説明

１０ 超音波探触子
１０ａ 超音波トランスデューサ
２０ 駆動信号生成部
３０、３０ａ 受信信号処理部
３１ 増幅器
３２、３２ａ、３２ｂ ＡＡフィルタ
３３ 直交検波器
４１、４１ａ、４１ｂ サンプリング部
４２、４２ａ、４２ｂ データ格納部
５０、５０ａ、５０ｂ 画像信号生成部
５１ 断面設定部
５２ 計算領域指定部
５３ 変換行列計算部
５４ 逆行列計算部
５５、５５ａ、５５ｂ 音源情報計算部
５６ 音場推定部
５７ 重み付け処理部
６１ Ｄ／Ａ変換器
６２ 表示部
７１ 操作卓
７２ 制御部
７３ 格納部

Claims (15)

1. 少なくとも１つの超音波トランスデューサに駆動信号を供給して被検体に超音波を送信すると共に、被検体から超音波エコーを受信した複数の超音波トランスデューサからそれぞれ出力される複数の受信信号を処理する送受信部と、
   前記送受信部によって処理された複数の受信信号をサンプリングすることにより、被検体の断面における複数のサンプリング点からの超音波エコーを表す受信信号行列を生成するサンプリング部と、
   被検体の断面における複数の位置からの超音波エコーが合成されるときの重み付け情報を表す変換行列を計算する変換行列計算手段と、
   前記変換行列計算手段によって計算された変換行列の逆行列を計算する逆行列計算手段と、
   前記逆行列計算手段によって計算された逆行列と前記サンプリング部によって生成された受信信号行列との積を計算することによって、被検体内の反射音源に関する情報を表す音源信号行列を求め、該音源信号行列に基づいて、被検体内の超音波画像を表す画像信号を生成する音源情報計算手段と、
   を具備する超音波撮像装置。
2. 前記送受信部が、複数の受信信号に対して、サンプリングにおけるエイリアシングを防止するためのフィルタ処理をそれぞれ施す複数のアンチ・エイリアシング・フィルタを含む、請求項１記載の超音波撮像装置。
3. 前記サンプリング部が、前記複数の超音波トランスデューサによって複数の時点で受信された超音波エコーに基づく複数の受信信号をサンプリングする、請求項１又は２記載の超音波撮像装置。
4. 前記変換行列計算手段が、被検体の断面において２次元的に配置された複数の位置からの超音波エコーが合成されるときの重み付け情報を表す変換行列を計算する、請求項１〜３のいずれか１項記載の超音波撮像装置。
5. 音場推定に基づいて、前記変換行列計算手段によって計算された変換行列に含まれている各要素に重み付けを与える重み付け処理部をさらに具備し、前記逆行列計算手段が、前記重み付け処理部によって重み付け処理された変換行列の逆行列を計算する、請求項１〜４のいずれか１項記載の超音波撮像装置。
6. 前記送受信部が、複数の受信信号に対して直交検波処理を施すことにより、それぞれの受信信号の実部及び虚部を求める複数の直交検波器を含み、
   前記サンプリング部が、受信信号の実部をサンプリングすることにより、被検体の断面における複数のサンプリング点からの超音波エコーの実部を表す第１の受信信号行列を生成する第１のサンプリング部と、受信信号の虚部をサンプリングすることにより、被検体の断面における複数のサンプリング点からの超音波エコーの虚部を表す第２の受信信号行列を生成する第２のサンプリング部とを含み、
   前記音源情報計算手段が、前記逆行列計算手段によって計算された逆行列と前記第１のサンプリング部によって生成された第１の受信信号行列との積を計算することによって、被検体内の反射音源に関する情報を表す第１の音源信号行列を求める第１の音源情報計算手段と、前記逆行列計算手段によって計算された逆行列と前記第２のサンプリング部によって生成された第２の受信信号行列との積を計算することによって、被検体内の反射音源に関する情報を表す第２の音源信号行列を求める第２の音源情報計算手段とを含み、
   第１の音源信号行列及び第２の音源信号行列に含まれている音源信号の実部及び虚部に基づいて音源信号の振幅を求め、音源信号の振幅に基づいて、被検体内の超音波画像を表す画像信号を生成する振幅演算手段をさらに具備する、請求項１〜５のいずれか１項記載の超音波撮像装置。
7. 前記音源情報計算手段が、被検体の断面を複数の領域に分けたときの各領域について、前記逆行列計算手段によって計算された逆行列と前記サンプリング部によって生成された受信信号行列との積を計算する、請求項１〜６のいずれか１項記載の超音波撮像装置。
8. 前記送受信部が、被検体内において順次選択された断面領域に送信フォーカスを形成するように、複数の超音波トランスデューサに複数の駆動信号をそれぞれ供給し、
   前記音源情報計算手段が、順次選択された断面領域について、前記逆行列計算手段によって計算された逆行列と前記サンプリング部によって生成された受信信号行列との積を計算する、請求項１〜６のいずれか１項記載の超音波撮像装置。
9. 前記変換行列計算手段が、ある深さからの受信信号に対して、その１／３の深さにある音源情報が加算されているとして変換行列を計算する、請求項１〜８のいずれか１項記載の超音波撮像装置。
10. 前記サンプリング部が、前記送受信部によって処理された複数の受信信号のサンプリング周波数を高めることにより、前記音源情報計算手段が、拡大された超音波画像を表す画像信号を生成する、請求項１〜９のいずれか１項記載の超音波撮像装置。
11. 前記音源情報計算手段が、選択された断面領域について、前記逆行列計算手段によって計算された逆行列と前記サンプリング部によって生成された受信信号行列との積を繰返し計算することにより、選択された断面領域における画像情報の時間的変化を表す画像信号を生成する、請求項１〜１０のいずれか１項記載の超音波撮像装置。
12. 前記変換行列計算手段が、ある深さからの受信信号に対して、各超音波トランスデューサの幅に対応する時間差を伴って受信される受信信号が加算されているとして変換行列を計算する、請求項１〜１１のいずれか１項記載の超音波撮像装置。
13. 前記送受信部が、複数の超音波トランスデューサに共通の駆動信号を供給する送信回路を含む、請求項１〜１２のいずれか１項記載の超音波撮像装置。
14. 前記送受信部が、複数の超音波トランスデューサに固有の遅延時間をそれぞれ与える複数の遅延回路と、前記複数の超音波トランスデューサに前記複数の遅延回路を介して駆動信号を供給する送信回路とを含む、請求項１〜１２のいずれか１項記載の超音波撮像装置。
15. 被検体の断面における複数の位置からの超音波エコーが合成されるときの重み付け情報を表す変換行列を計算するステップ（ａ）と、
    ステップ（ａ）において計算された変換行列の逆行列を計算するステップ（ｂ）と、
    少なくとも１つの超音波トランスデューサに駆動信号を供給して被検体に超音波を送信すると共に、被検体から超音波エコーを受信した複数の超音波トランスデューサからそれぞれ出力される複数の受信信号を処理するステップ（ｃ）と、
    ステップ（ｃ）において処理された複数の受信信号をサンプリングすることにより、被検体の断面における複数のサンプリング点からの超音波エコーを表す受信信号行列を生成するステップ（ｄ）と、
    ステップ（ｂ）において計算された逆行列とステップ（ｄ）において生成された受信信号行列との積を計算することによって、被検体内の反射音源に関する情報を表す音源信号行列を求め、該音源信号行列に基づいて、被検体内の超音波画像を表す画像信号を生成するステップ（ｅ）と、
    を具備する超音波撮像方法。

Images