(57)【要約】【課題】強い反射体がある場合にも高画質な撮像を行うこと、及び、高い音圧を当ててはいけない部位がある場合にも全体の音圧を低下さないで撮像や治療を行うこと。【解決手段】目標とするビームを設定するためのデータを取得し、そのデータから回避する部位の位置と強度を検出し、位置と強度を目標ビーム形状に換算し、目標ビーム形状に沿ったビームを形成するフォーカスデータを算出し、そのフォーカスデータを用いて画像生成あるいは治療を行う。【選択図】図２
2010104584.tif 000002
【特許請求の範囲】
【請求項１】
超音波を送信あるいは受信する複数の素子を具備する探触子と、
前記複数の素子が被検体へ超音波信号を送信する際の指向性を付与する送信ビームフォーマと、
前記複数の素子が受信した超音波信号を指向性をともなって加算する受信ビームフォーマと、
前記受信ビームフォーマの出力した信号を信号処理して画像化する信号処理部と、
前記信号処理部の出力した画像を表示する表示手段とを備え、
前記信号処理部は、目標ビーム形状を設定する目標ビーム形状設定部と、前記目標ビーム形状が入力され、当該目標ビーム形状に沿ったビームを生成するフォーカスデータを算出するフォーカスデータ生成部と、画像を生成する画像生成部とを含むことを特徴とする超音波装置。
【請求項２】
請求項１に記載の超音波装置において、
前記目標ビーム形状設定部は、目標とする音場を設定するためのデータを取得する目標設定用データ取得部と、前記目標設定用データから回避する部位の位置と強度を検出する回避部検出部と、前記回避する部位の位置と強度を目標ビーム形状に換算する目標ビーム形状換算部とを具備し、
前記送信ビームフォーマと前記受信ビームフォーマの少なくとも一方は、前記フォーカスデータ生成部が出力したフォーカスデータを用いてビームを生成することを特徴とする超音波装置。
【請求項３】
請求項１に記載の超音波装置において、
さらに、情報を格納するメモリを具備し、前記メモリは前記複数の素子ごとの受信信号を格納する受信信号格納部を具備し、
前記目標ビーム形状設定部は、目標とする音場を設定するためのデータを取得する目標設定用データ取得部と、前記目標設定用データから回避する部位の位置と強度を検出する回避部検出部と、前記回避する部位の位置と強度を目標ビーム形状に換算する目標ビーム形状換算部とを具備し、
前記目標とする音場を設定するためのデータは、前記受信信号格納部に格納された、前記複数の素子ごとの受信信号であり、
前記画像生成部は、前記受信信号格納部から前記複数の素子ごとの受信信号を読み出し、前記フォーカスデータ生成部が出力したフォーカスデータを用いて画像を再構成することを特徴とする超音波装置。
【請求項４】
請求項１〜３のいずれか１項に記載の超音波装置において、
前記フォーカスデータ生成部は、フォーカスデータから音場への変換を表す演算子をＴとし、前記目標ビーム形状を表す関数をＷとして、演算子Ｔ-1Ｗ-1ＷＴの固有関数φに比例したフォーカスデータを出力することを特徴とする超音波装置。
【請求項５】
請求項１〜３のいずれか１項に記載の超音波装置において、
前記フォーカスデータ生成部は、フォーカスデータから音場への変換を表す演算子をＴとし、前記目標ビーム形状を表す関数をＷとして、演算子Ｔ+Ｗ+ＷＴＴ+Ｗ+ＷＴの固有関数φTに比例した送信用フォーカスデータとＴ+Ｗ+ＷＴφTに比例した受信用フォーカスデータφTを出力することを特徴とする超音波装置。
【請求項６】
請求項４又は５に記載の超音波装置において、
前記表示手段は、前記目標音場決定用データ、前記目標ビーム形状、前記フォーカスデータのうちの少なくとも１つを表示し、
撮像者の操作によって、前記目標ビーム形状、前記フォーカスデータのうち少なくとも１つを入力する入力手段を有することを特徴とする超音波装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波装置に関し、特に医療用に適した超音波装置に関する。
【背景技術】
【０００２】
医用超音波装置において、高画質化と医療上の安全性の確保は重要な課題である。特に、強い反射体があるとそれが雑音源となり、実像に虚像が重なるという画質の低下が頻繁にみられるという問題や、ＨＩＦＵ（高熱焼灼療法）において、血管など、強い音圧を当ててはいけない部分がある場合、全体の音圧を下げなければならず、焼灼部位の音圧と範囲が小さくなり、焼灼時間が長くなるという問題の改善は、大きな課題である。
【０００３】
強い反射体の例としては肋骨、横隔膜、ＨＩＦＵにおける金属性のプローブなどがある。強い反射体による画質の劣化という問題は日常的なものであり、これを解決する方法としては、プローブ位置を変えると実像と虚像の相対位置が変わることを用いて虚像を判別したり、プローブ位置と時間を変えて平均をとり、消去するというマニュアル的な方法や、コンパウンド技術などの平均処理によって消去するという技術的な方法が実施されている。特許文献１には、予め強反射体の位置がわかっている場合に、画像処理で輝度を局所的に調整して、画像全体で生体の構造が判別できるように画質を回復する方法が開示されている。
【０００４】
【特許文献１】特開２０００−３７３９３号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかし、上記のいずれの従来方法でも、ＳＮ比は改善されておらず、通常行われているマニュアル的及び技術的な平均処理方法ではさらに、時間分解能が低下するという問題もある。また、指定した位置のみ音圧を小さくする方法は開示されていない。
【０００６】
本発明の目的は、強い反射体がある場合にも時間分解能を低下させずにＳＮ比を改善して高画質な撮像を行うこと、及び、高い音圧を当ててはいけない部位がある場合にも全体の音圧を低下させずに撮像や治療を行うことのできる超音波装置を提供することにある。
【課題を解決するための手段】
【０００７】
本発明では、音場設計により高画質化と安全性の確保を図る。本発明の超音波装置は、超音波を送信あるいは受信する複数の素子を具備する探触子と、複数の素子が被検体へ超音波信号を送信する際の指向性を付与する送信ビームフォーマと、複数の素子が受信した超音波信号を指向性をともなって加算する受信ビームフォーマと、受信ビームフォーマの出力した信号を信号処理して画像化する信号処理部と、信号処理部の出力した画像を表示する表示手段とを備える。
【０００８】
信号処理部は、目標ビーム形状を設定する目標ビーム形状設定部と、目標ビーム形状を入力し、その目標ビーム形状に沿ったビームを生成するフォーカスデータを算出するフォーカスデータ生成部とを含む。送信ビームフォーマと受信ビームフォーマの少なくとも一方は、フォーカスデータ生成部が出力したフォーカスデータを用いてビームを生成する。あるいは、メモリに探触子の素子毎の受信信号を記憶しておき、その受信信号にフォーカスデータ生成部が出力したフォーカスデータを適用して画像を再構成する。
【発明の効果】
【０００９】
本発明の超音波装置によれば、意図した位置の音圧を抑制した音場を形成することができ、ノイズ源となる強い反射体がある場合にも時間分解能を低下させずにＳＮを向上させて高画質な画像を得ることができる。また、高い音圧を当ててはいけない部位がある場合にも全体の音圧を低下さないで撮像や治療を行うことができる。
【発明を実施するための最良の形態】
【００１０】
最初に、図１Ａから図９を用いて本発明の第１の実施例を説明する。
図１Ａは、本発明の第１の実施例の超音波装置の構成例を表すブロック図である。探触子１は複数の素子を備える。装置本体２は、送信ビームフォーマ３、増幅手段４、受信ビームフォーマ５、信号処理部６、メモリ７、表示手段８、入力手段９、制御部１０を有する。信号処理部６は、目標ビーム形状設定部６１、フォーカスデータ生成部６２、画像生成部６３を含む。さらに、目標ビーム形状設定部６１は、目標音場決定用データ取得部６１１、回避部検出部６１２、目標ビーム形状換算部６１３を含み、フォーカスデータ生成部６２は、目標ビーム形状入力部６２１、フォーカスデータ算出部６２２を含む。
【００１１】
送信ビームフォーマ３で生成された超音波パルスが探触子１から生体に送信され、生体から反射した超音波を探触子１が受信する。受信信号は増幅手段４に入力されて増幅され、受信ビームフォーマ５が整相加算する。信号処理部６の画像生成部６３には受信ビームフォーマ５が出力した受信信号が入力され、画像化が行われる。作成された画像はメモリ７に格納された後、読み出し・補間されて、表示手段８に表示される。なお、これらの処理は制御部１０によって制御される。
【００１２】
送信ビームフォーマ３及び受信ビームフォーマ５におけるビームフォーミングに用いるフォーカスデータは、予めメモリ７に格納されており、撮像時にメモリ７から読み出されてもよい。あるいは、メモリ７に格納されたデータである目標音場決定用データから計算してもよい。ここで、フォーカスデータとは、超音波信号を送受信する際の指向性を付与するために複数の素子に与える時間遅延とパルス強度、数値的な表現では複素振幅を意味する。目標音場決定用データは、たとえば前記の撮像を行った結果をメモリ７に一旦格納したものでもよい。この目標音場決定用データは、プレスキャン結果と考えるとその意味が理解しやすい。フォーカスデータを目標音場決定用データから計算する場合、目標音場決定用データがメモリ７に一旦格納されていることと、計算したフォーカスデータもメモリ７に一旦格納することを仮定すると、信号の流れは図１Ａ中に太矢印線で示したようになる。具体的には、メモリ７から目標ビーム形状設定部６１が目標音場決定用データを読み出して目標ビーム形状に換算してフォーカスデータ生成部６２に出力し、フォーカスデータ生成部６２が入力した目標ビーム形状をフォーカスデータに換算してメモリ７に出力する、という流れになる。
【００１３】
図１Ｂは、目標音場決定用データからフォーカスデータを計算し、フォーカスデータを使用する処理の流れを説明するフローチャートである。
【００１４】
まず、撮像が開始されると（開始）、フレーム数を表す整数ｉを０に設定し（Ｓ１１１）、メモリに予め格納してあるフォーカスデータを用いて撮像を行って、フレーム０の画像を得る（Ｓ１１２）。次に、電源が落とされるなどして撮影が終了しない限り（Ｓ１１６）、ｉをｉ＋１とし（Ｓ１１３）、フレームi−１の画像を目標音場決定用データとしてフォーカスデータを計算し（Ｓ１１４）、計算したフォーカスデータを送信ビームフォーマ、受信ビームフォーマの片方あるいは両方で用いてフレームiの撮像を行う（Ｓ１１５）という処理を続ける。
【００１５】
以下、図２から図８を用いて詳細に説明する。まず、図２から図４Ａのフローチャートを用いて信号処理部６の処理を説明する。
【００１６】
図２は、信号処理部６の処理を説明するフローチャートである。信号処理部６が起動されると（開始）、目標ビーム形状設定部６１が目標とするビームの形状を設定し（Ｓ２０１）、フォーカスデータ生成部６２が目標ビーム形状に近いビームを生成するフォーカスデータを算出し（Ｓ２０２）、算出したフォーカスデータを用いて撮像対象に対して超音波を送受信して受信データを画像生成部６３が入力し、画像を生成する（Ｓ２０３）。
【００１７】
図３は、信号処理部６における処理のうち、目標とするビームの形状を設定する処理（Ｓ２０１）を説明するフローチャートである。目標とするビームの形状を設定する処理（Ｓ２０１）が開始されると（Ｓ２０１開始）、目標音場決定用データ取得部６１１が目標とする音場を決定するためのデータを取得し（Ｓ２０１１）、回避部検出部６１２が回避する部位の位置と強度を検出し（Ｓ２０１２）、目標ビーム形状換算部６１３が回避する部位の位置と強度を目標とするビームの形状に換算し、目標ビーム形状を設定する（Ｓ２０１３）。
【００１８】
図４Ａは、信号処理部６における処理のうち、フォーカスデータ生成部６２の処理（Ｓ２０２）の一例を説明するフローチャートである。フォーカスデータ生成部６２の処理が開始されると（Ｓ２０２開始）、フォーカスデータから音場への変換を表す演算子Ｔ、及び、目標ビーム形状を表す関数Ｗを入力し（Ｓ２０２１）、演算子Ｔ+Ｗ+ＷＴを計算し（Ｓ２０２２）、演算子Ｔ+Ｗ+ＷＴの固有関数φnを計算し（Ｓ２０２３）、最大固有値をもつφnをフォーカスデータに設定する（Ｓ２０２４）。なお、すべての演算子は離散（行列表現あるいはテンソル表現）でもよいし、連続（関数表現）でもよい。
【００１９】
以下に、図４Ａにおけるフォーカスデータから音場への変換を表す演算子Ｔ、目標ビーム形状を表す関数Ｗ、及び、演算子Ｔ+Ｗ+ＷＴの固有関数φnの一例を示す。添え字の混乱を避けるため、ここでは固有関数φnをφと記した。
【００２０】
【数１】
2010104584.tif 000003
【００２１】
本例では、すべての演算子を離散演算子とし、ビーム形状を表す空間を１次元とし、探触子１の振動素子の並びを１次元とし、ビーム形状を表す空間の離散化数をｍ、振動素子数をｎとし、ビーム形状を表す空間の位置座標をｘ’1＝（ｘ’1，ｚ）〜ｘ’m＝（ｘ’m，ｚ）、振動素子の位置座標を（ｘ1，０）〜（ｘn，０）とした。このとき、フォーカスデータから音場への変換を表す演算子Ｔはｍ行ｎ列の行列となる。行列Ｔのi行j列の要素はj番目の振動素子がｉ番目の空間位置に作る音場の強度である。なお、本例の演算子Ｔはｘ’1／ｚ→０の極限でフーリエ変換を表す行列となるが、これは、近軸近似では、ビームの指向性は、振動素子に与えるフォーカスデータのフーリエ変換になることに相当する。目標ビーム形状を表す関数Ｗは、ｍ行ｍ列の対角行列となる。対角要素は、目標とするビーム形状がたとえば図４Ｂに示す４０１のような形、具体的には、座標（ｘ’2，ｚ）にメインビーム４０１１を持ち、座標（ｘ’m-1，ｚ）に抑制部４０１２を持つ形とする場合、ｘ’2行ｘ’2列要素がたとえば２と大きく、ｘ’m-1行ｘ’m-1列の要素がたとえば０と小さく、他の要素はたとえば１とその中間である行列となる。このとき、演算子Ｔ+Ｗ+ＷＴの固有関数φはｎ行複素ベクトルとなる。各要素は一般的に複素数であり、探触子４０２の振動素子４０２１〜４０２ｎのフォーカスデータに相当する。具体的には、ｉ番目の要素の振幅ａi及び位相θiは、振動素子４０２ｉがビームフォーミングする時の音響強度、及び、（遅延時間）・（音速）に相当する。
【００２２】
下記の数式を用いて、形状Ｗのビームを生成したい場合、フォーカスデータが演算子Ｔ+Ｗ+ＷＴの固有関数φとして求まるしくみを説明する。
【００２３】
【数２】
2010104584.tif 000004
【００２４】
いま、フォーカスデータφでビームを作るとする。この場合、フォーカスデータから音場への変換を表す演算子がＴで表されているため、ビームの形状は、式（１）に示すＢのようになる。ここで、フォーカスデータφによって作られるビーム形状Ｂと、目標とするビーム形状Ｗの一致度は、演算子の積ＷＢの自乗||ＷＢ||2で表される。このことは、今、演算子Ｗ及びＴが行列であり、２つのベクトルの一致度がベクトルの内積の自乗で表されることを考えると理解される。一方、探触子から出力される音圧の合計はフォーカスデータφのノルムの自乗||φ||2で表される。探触子から出力される音圧の合計は、近似的に、フォーカス点での音響強度を表すとみなせる。
【００２５】
ここで、フォーカス点での音響強度を一定にした条件下で可能な限り目標とするビーム形状Ｗに一致するビームを作るフォーカスデータを求めることを目標とする。するとこの問題を変分法で定式化すると、評価関数として式（２）の第１行が得られる。λはラグランジュの未定係数であり、係数のかかった項は制約条件である。評価関数（２）の第１行は、式（１）を代入して整理すると、式（２）の第３行のように整理され、変分をとると式（３）のようになる。微小変化δφ+は任意であるので、変分が０という条件式（４）の第１行を変形することで、問題は、式（４）の第３行に得られた、φに関する固有値問題に帰着することがわかる。すなわち、演算子Ｔ+Ｗ+ＷＴの固有関数φを求めてフォーカスデータとして用いれば、フォーカス点での音響強度を一定にした条件下で、可能な限り形状Ｗに近いビームが生成される。
【００２６】
以下、図５Ａ〜５Ｃの概念図及び図６の説明図を用いて、信号処理部６の処理手順を説明する。図５Ａ〜５Ｃは、本発明の超音波装置による撮像手順の例を説明する概念図である。
【００２７】
図５Ａは、目標とする音場を決定するためのデータを取得する処理（Ｓ２０１１）の一例を説明する概念図である。５０１は探触子、５０２は人体などの被検体、５０３は撮像領域、５０３１から５０３４はそれぞれ１回の送受信におけるビームの偏向方向、５０４１から５０４４は前記ビームの偏向方向５０３１から５０３４に超音波を送受する際のビームの指向性である。横軸が方位角、縦軸が指向性を表し、本図ではメインビームとサイドビームを有する。
【００２８】
目標とする音場を決定するためのデータを取得する処理（Ｓ２０１１）が開始されると、たとえば、目標音場決定用データ取得部６１１は、探触子５０１から撮像領域５０３内に偏向方向５０３１から５０３４を設定し、それぞれの偏向方向に対して、指向性５０４１から５０４４をもつビームを生成して超音波信号を取得して画像化し、目標音場決定用データとする。ここで、本例では説明のわかりやすさのために、超音波信号を送受する例を挙げたが、目標音場決定用データ取得部６１１は、図１Ｂに例を挙げたように、１フレーム前の撮像結果をメモリに一旦格納しておき、メモリから１フレーム前の撮像結果を読み出して目標音場決定用データとしてもよい。
【００２９】
図５Ｂは回避する部位の位置と強度を検出し、目標ビーム形状を設定する処理（Ｓ２０１２，Ｓ２０１３）の一例を説明する概念図である。撮像領域５０３内の画像は目標音場決定用データで、画像に写っている５０５１は画像化したい被検体、５０５２は肋骨などの、被検体の撮像にとってノイズ源となる強反射体である。５０６１から５０６４は、図５Ａの偏向方向５０３１から５０３４に超音波を送受する際に目標とするビーム形状である。くわしくは、５０６１１から５０６４１は、偏向方向を反映した、メインビームと呼ばれる音圧の大きい領域、５０６１２から５０６４２は、強反射体５０５２の位置を反映した、音圧の抑制部である。
【００３０】
回避する部位の位置と強度を検出する処理（Ｓ２０１２）が開始されると、回避部検出部６１２が、目標音場決定用データから、ノイズ源となる強反射体５０５２の位置と強度を検出する。次に、目標ビーム形状換算部６１３が回避する部位の位置と強度を目標ビーム形状に換算して、偏向方向５０３１から５０３４によって位置が変化するメインビーム５０６１１から５０６４１、及び、偏向方向５０３１から５０３４によって位置が変化しない抑制領域５０６１２から５０６４２を有する目標ビーム形状５０６１から５０６４を設定する（Ｓ２０１２）。
【００３１】
なお、以上の処理（Ｓ２０１２，Ｓ２０１３）は、ユーザによってマニュアル的に設定されてもよいし、装置の内部処理によって自動的に設定されてもよい。具体的には、目標音場決定用データを表示手段８が表示し、ユーザが強反射体５０５２の位置と強度を視認し、ユーザの知識に基づいて目標ビーム形状を判断し、入力手段９を用いて手動で入力してもよい。あるいは、目標音場決定用データを信号処理部６が読み込み、信号処理部がパターン認識などの画像処理を行って強反射体５０５２の位置と強度を自働的に検出し、回避する部位の影響を抑制する目標ビーム形状５０６１から５０６４を自動的に判断し、設定してもよい。
【００３２】
図５Ｃは、画像生成部６３が画像を生成する処理（Ｓ２０３）の一例を説明する概念図である。５０７１から５０７４は、偏向方向５０３１から５０３４に超音波を送受する際のビームの指向性である。
【００３３】
画像を生成する処理（Ｓ２０３）が開始されると、フォーカスデータ生成部６２が目標ビーム形状５０６１から５０６４を入力して計算したフォーカスデータを、送信ビームフォーマ３、及び／又は、受信ビームフォーマ５に入力し、探触子１を介して、偏向方向に音圧の大きいメインビーム５０７１１から５０７４１を持ち、強反射体５０５２の方向に抑制部５０７１２から５０７４２を持つ前記指向性５０７１から５０７４を有する超音波ビームを被検体に送受信して、画像を生成する。
【００３４】
図６は、目標ビーム形状換算部６１３が、回避する部位の位置及び強度を目標とするビームの形状に換算する際の、回避する部位の位置及び強度とビーム形状との数値的な関係の一例を説明する図である。
【００３５】
図６（ａ）は超音波装置の被検体を構成する散乱体のインピーダンス、各種散乱体が隣接した場合の超音波の反射率及び画像上の輝度差を示す表である。代表的な散乱体として、生体内の点反射体、血管壁などの連続した反射体、骨・金属などを考えた。生体内の点反射体のインピーダンスはおよそ１．３〜１．６MRayl、血管壁などの連続した反射体はおよそ１．６〜１．７MRayl、骨・金属などはおよそ７〜１２MRaylである。よって点反射体と連続した反射体が隣接した場合の反射率はおよそ〜２％、連続した反射体と骨・金属などが隣接した場合の反射率はおよそ７０〜８０％である。また、画像上の輝度差はそれぞれおよそ２０〜３０ｄＢ、およそ３０〜４０ｄＢとなる。
【００３６】
よって、生体内の点反射体を、骨や金属に起因するノイズに埋もれさせずに撮像するためには、図６（ｂ）に示したように、偏向方向５０３１に超音波を送受する際のビームの指向性５０７１のメインビーム５０７１１と抑制部５０７１２の強度差が（２０〜３０ｄＢ）＋（３０〜４０）ｄＢ、すなわちおよそ５０〜７０ｄＢあればよい。なお、抑制部を設定しない平均的な指向性において、メインビーム５０７１１とサイドローブレベル（指向性の低い方の集束値）の音圧差は４０ｄＢ程度である。
【００３７】
上記を本発明の装置の処理にあてはめて説明すると、回避部検出部６１２が、目標音場決定用データから、連続した反射体に比べて３０ｄＢ高い強反射体５０５２を検出した場合、目標ビーム形状換算部６１３が、強反射体の位置に、メインビームより５０から６０ｄＢ音圧が低い抑制部を有する目標ビーム形状５０６１から５０６４を設定すれば、強反射体に起因するノイズに埋もれさせずに生体内の点反射体を撮像できる。
【００３８】
図７は、フォーカスデータ生成部の入出力の計算例を示す図である。上段が、フォーカスデータ生成部６２における、目標ビーム形状を表す関数Ｗ、下段が、フォーカスデータ生成部が図４に示した処理によって計算したフォーカスデータを用いて形成される指向性の一例である。上段の縦軸は線形表示、下段の縦軸は対数表示、横軸はともに方位角である。指定した部分５０６１２のみがその周辺より有意に音圧が小さくなっており、メインビーム５０７１１と抑制部５０７１２の音圧差がおよそ５０ｄＢになっている。この値は図６を用いて説明したように、生体内の点反射体を、骨や金属に起因するノイズに埋もれさせずに撮像するために十分な音圧差である。なお本例では、図４に示した処理Ｓ２０２１における、フォーカスデータから音場への変換を表す演算子Ｔを、離散フーリエ変換とした。
【００３９】
図８は、本発明の超音波装置の効果を検証するシミュレーション画像である。図８（ａ）はファントム、図８（ｂ１）／（ｂ２）は従来法／本発明の方法によるシミュレーションに使ったビームの指向性、図８（ｃ１）／（ｃ２）は従来法／本発明の方法によるシミュレーション画像である。
【００４０】
ファントムは、薄板状の強反射体がある、輝度の低い１次元シストファントムで、平均反射強度が−３０ｄＢの背景８０１１の中に、０ｄＢの強反射体８０１２、−３４．０ｄＢ、−４０．０ｄＢの２つのシスト８０１３，８０１４を含む。シストの直径はビーム幅の約２倍である。図８（ａ）に示したファントム及び図８（ｂ）に示したビームは、計算の簡易化のため１次元としたが、本発明の方法は２次元、３次元の撮像対象及びビームに対しても適用可能である。また、図８（ａ）のシスト８０１３，８０１４を見やすくするため、表示レンジを狭くした抜き出し図を右側に記載した。図８（ｂ２）に示した本発明の方法によるビームのうち、強反射体の影響を抑制するための抑制部８０２１は効果を保障するために広くとり、本例では従来法のビームより約１０ｄＢ音圧を小さくしてある。図８（ａ）と図８（ｂ１）（ｂ２）の横スケールは同じで、図８（ｃ１）（ｃ２）は拡大図である。
【００４１】
シミュレーション画像では、図８（ｃ１）に示した従来法では、強反射体が生じさせたかぶりが見られるのに対して、図８（ｃ２）に示した本発明の方法ではかぶりが抑制されているのが見て取れる。かぶりが抑制された結果、特に８０１３では、形状の認識も容易になっている。なお、本シミュレーションでは、強反射体は８０１２に示すように薄板状、ビーム（ｂ２）は１次元連続波であるため、シストにのっているかぶりはサイドローブの形状そのものであり、ノイズであるかぶりと信号である円形の輝度変化が視認でも分離・区別できる。しかし、通常、強反射体の形状は単純なものとは限らず、かぶりも複雑なものになるのが一般的である。その結果、複雑な形状をしたかぶりのノイズによって信号である構造の形状が判別できなくなる。すなわち、本発明の方法によれば、ほとんどの場合に、図８（ｃ２）に示す以上の形状の視認性の向上が得られる。
【００４２】
本発明の方法では、ノイズであるかぶりが抑制されて信号Ｓ／Ｎ比が向上し、従来の方法では識別できない形状８０１３が識別できるようになっている。なお、従来の方法による画像（ｃ１）は、本発明における、目標音場決定用データとして用いることもできる。その場合、目標音場決定用データ取得部６１１が目標音場決定用データ８０１１を読み込み、回避部検出部６１２が周辺に比べて３０ｄＢ輝度の高い強反射体８０１２を方位角８０２１に検出し、フォーカスデータ生成部６２が強反射体８０１２の存在する方位角でメインビームより５０ｄＢ小さい抑制部８０２１をもつ受信ビーム（ｂ２）を作るフォーカスデータを計算して画像生成部６３が画像を生成する。その結果、かぶりが抑制されて信号Ｓ／Ｎ比が向上し、従来の方法では識別できない形状８０１３が識別できる画像（ｃ２）を得ることができる。
【００４３】
以上のような構成によれば、被検体からの信号はメインビーム５０７１１から５０７４１によって得られ、強反射体５０５２からのノイズは抑制部５０７１２から５０７４２によって除去されるため、ノイズ源となる強反射体がある場合にも、ＳＮ比の高い、高画質な画像が得られる。
【００４４】
より詳しくは、以上のような構成、特に図１Ｂの構成によれば、目標音場決定用データからフォーカスデータをリアルタイムに計算してノイズ源の影響を抑制した画像を得ることができ、画像のフレームレートを落とさずに、血流や体動によって動くノイズ源の影響をも抑制できる。
【００４５】
また、詳しくは、自動処理によりノイズ源を検出するため、ノイズ源に関する知識のない、予期しないノイズ源に対しても、手間をかけずに、その影響を抑制できる。
【００４６】
さらに、以上のような構成によれば、上述の例で強反射体としたところを血管などの、強い音圧を当ててはいけない領域、偏向方向をＨＩＦＵ治療の焼灼部位に置き換えて考えると、血管部での音圧が抑制され、血管を損傷せずに従来の方法より焼灼部位での音圧を強くすることができるため、短時間で、かつ、医療上の安全性を確保した治療を行うことができる。ただしこの場合は、図１ＢのステップＳ１１５において、計算したフォーカスデータを送信ビームフォーマ、あるいは、送信ビームフォーマ及び受信ビームフォーマで用いる。
【００４７】
次に、図９Ａ及び図９Ｂを用いて、本発明の第２の実施例を説明する。第１の実施例では、図１Ｂに示したように、リアルタイム処理で、i番目のフレームの画像を目標音場決定用データとしてフォーカスデータを計算し、i＋１番目のフレームの送信ビームフォーマ、受信ビームフォーマのフォーカスデータの片方あるいは両方で用いる構成を説明したが、目標音場決定用データから計算されたフォーカスデータは、送信ビームフォーマ３のみが用いてもよいし、受信ビームフォーマ５のみが用いてもよいし、送信ビームフォーマ３及び受信ビームフォーマ５の両方が用いてもよい。また、時間的にも、リアルタイムに撮像に反映されてもよいし、オフラインで非リアルタイムの画像再構成に用いられてもよい。ただし、オフライン処理を行う後者の場合、メモリ７は送信ビームフォーマ３及び受信ビームフォーマ５のフォーカスデータ及び送信信号、受信信号を格納しているとする。
【００４８】
第１の実施例で説明した以外の処理の例を、図９Ａ及び図９Ｂに示す。図９Ａは、リアルタイム処理で、i番目のフレームの画像を目標音場決定用データとしてフォーカスデータを計算し（Ｓ９１４’）、同じi番目のフレームの受信ビームフォーマで用いて画像を再構成する（Ｓ９１５’）場合である。この場合、フレームiの、振動素子ごとの受信信号、及び、フォーカスデータをメモリに格納する必要がある（Ｓ９１３、９１３’）。なお、Ｓ９１０（Ｓ９１１〜Ｓ９１５）は最初の（０番目の）フレームの処理、Ｓ９１０’（Ｓ９１１’〜Ｓ９１５’）はそれ以降のフレームの処理であり、対応する処理が多いため、対応する処理にはダッシュつきの番号を付けた。
【００４９】
図９Ｂは、オフライン処理で、i番目のフレームの画像を目標音場決定用データとしてフォーカスデータを計算し、同じi番目のフレームの受信ビームフォーマで用いて画像を再構成する場合である。オフライン処理では、計算したフォーカスデータは受信ビームフォーマにしか使えない。ステップＳ９２１が物理的に撮像対象に音を送受信している過程で、ステップＳ９２２（Ｓ９２３〜Ｓ９２７）がオフライン処理による画像再構成の過程である。オフライン処理Ｓ９２２では、フレーム数iを更新しつつ、フレームごとにメモリから振動素子ごとの受信信号、フォーカスデータ、画像を読み込み（Ｓ９２４）、フォーカスデータを計算し（Ｓ９２５）、それを用いて同じｉ番目のフレームの画像を再構成している（Ｓ９２６）。
【００５０】
以上のような構成、特に図９Ａに示した構成によれば、フレームレートは第1の実施例に劣るが、動きの早いノイズ源に対して高い抑制効果が得られ、ＳＮ比の高い、高画質な画像が得られる。また、以上のような構成、特に図９Ｂに示した構成によれば、撮像装置の演算能力が不十分な場合にも、ノイズ源の影響を抑制でき、ＳＮ比の高い、高画質な画像が得られる。
【００５１】
次に、図１０を用いて、本発明の第３の実施例を説明する。
図１０は、本発明の第３の実施例における、目標ビーム形状設定部６１の処理の一例を示す概念図である。８は表示部、１００１及び１００２は撮像者の指で、入力手段９の一例を示している。すなわち、本実施例は入力手段が入力ボタン表示１００３及び１００４を有するタッチパネル形式である例を示している。表示部８に目標音場決定用データとして、撮像結果５０３が表示されている。撮像結果５０３内にはノイズ源となる強反射体５０５２及び描出したい部位５０５３がある。撮像者はタッチパネル形式で回避部、目標ビーム形状のうち１つ以上を入力する。
【００５２】
回避部を入力する例を説明すると、たとえば、撮像者は目標音場決定用データとして表示された撮像結果５０３を目視し、強反射体５０５２を視認して回避部とみなし、タッチパネル上の強反射体５０５２の位置に接触し（１００１）、続いて回避部を指定したことを示す入力ボタン表示１００３に接触（１００１’）する。すると、回避部検出部６１２が内部処理を行って接触部位の位置及び信号強度をたとえば偏向角３０±５度の位置に０ｄＢと算出し、回避部の位置及び信号強度とする。さらに撮像者は、目標音場決定用データとして表示された撮像結果５０３を目視し、描出しようとする部位５０５３を視認して描出部とみなし、タッチパネル上の描出部５０５３の位置に接触し（１００２）、続いて抽出部を指定したことを示す入力ボタン表示１００４に接触してもよい。この場合、回避部検出部６１２が内部処理を行って描出部の位置及び信号強度をたとえば偏向角−２５±２度の位置に−５０ｄＢと算出し、描出部の位置及び信号強度とする。なお、タッチパネル上への接触は５０５３の例に示すように点的でもよいし、５０５２の例に示すように線的でもよい。
【００５３】
回避部検出部６１２は、接触を入力し、指定部の範囲と信号強度を出力する検出部を有するとする。たとえば前記検出部は、タッチパネル上への接触画素の位置を検出して、たとえば接触画素の平均値を中心、接触画素の分散の２倍を近傍距離として、目標音場決定用データの輝度の変化が少ない、たとえば画像全体５０３の最大輝度と最低輝度の１／１０以下のばらつきをもつ部位を同一の構造物に属すると判定して検出して、指定範囲、たとえば上記の例では接触１００１に対しては偏向角３０±５度の範囲を出力し、さらに、前記部位の代表輝度、たとえば平均輝度、あるいは平均輝度から輝度の分散の整数倍をさし引いた輝度を指定部の信号強度として出力する。
【００５４】
次に、以上のように、回避部検出部６１２が検出した、回避部、描出部の位置及び信号強度（それぞれ偏向角３０±５度の位置に０ｄＢ、偏向角-２５±２度の位置に−５０ｄＢ）を目標ビーム形状換算部６１３が入力して目標ビーム形状に換算する。上述した例では、メインビームの方向によらず、偏向角３０±５度の位置にメインビームより−５０ｄＢ小さい抑制部を持つ指向性５０６１から５０６４（図５）を目標ビーム形状とする。
【００５５】
なお、図１０ではタッチパネル形式の例を示したが、入力手段を限定するものではない。その例としては、ペンなどの器物、キーボード、マウスのうちの１つ以上でもよい。
【００５６】
以上のような構成によれば、ＨＩＦＵの治療用探触子など、予めわかっているノイズ源の影響を確実に回避することができる。さらに、撮像者の意図をより柔軟に反映することができる。
【００５７】
次に、図１１を用いて、本発明の第４の実施例を説明する。
図１１は、本発明の第４の実施例における、フォーカスデータ生成部６２の処理（Ｓ２０２）の一例を説明するフローチャートである。フォーカスデータ生成部６２の処理が開始されると（Ｓ２０２開始）、フォーカスデータから音場への変換を表す演算子Ｔ、及び、目標ビーム形状を表す関数Ｗを入力し（Ｓ２０２１’）、演算子Ｔ+Ｗ+ＷＴＴ+Ｗ+ＷＴを計算し（Ｓ２０２２’）、演算子Ｔ+Ｗ+ＷＴＴ+Ｗ+ＷＴの固有関数φnを計算し（Ｓ２０２３’）、最大固有値をもつφnを送信用フォーカスデータφTに設定し（Ｓ２０２４’）、Ｔ+Ｗ+ＷＴφTを受信用フォーカスデータφRに設定する（Ｓ２０２５’）。
【００５８】
下記の数式を用いて、図１１Ａにおいて、形状Ｗのビームを生成したい場合、送信用及び受信用のフォーカスデータが、演算子Ｔ+Ｗ+ＷＴＴ+Ｗ+ＷＴの固有関数φ、及びＴ+Ｗ+ＷＴφとして求まるしくみを説明する。
【００５９】
【数３】
2010104584.tif 000005
【００６０】
上記数式（５）〜（８）は、第１の実施例で説明した数式（１）〜（４）と対応している。送信用フォーカスデータφTと受信用フォーカスデータφRでビームを作るとする。この場合、フォーカスデータから音場への変換を表す演算子がＴで表されているため、式（５）に示すように、送信ビームの形状はＢT、受信ビームの形状はＢRのようになる。ここで、フォーカスデータφT及びφRによって作られる送受のビーム形状ＢR，ＢTと、目標とするビーム形状Ｗの一致度は、それぞれのビーム形状と目標の形状の積、ＷＢR及びＷＢTの積で表される。一方、探触子から出力される送信音圧の合計は送信用フォーカスデータφTのノルムの自乗||φT||2、受信音圧の合計は受信用フォーカスデータφRのノルムの自乗||φR||2で表される。
【００６１】
ここで、フォーカス点での音響強度を一定にした条件下で可能な限り目標とするビーム形状Ｗに一致するビームを作るフォーカスデータを求めることを目標とする。すると、この問題を変分法で定式化すると、評価関数として式（６）の第１行が得られる。λT及びλRはラグランジュの未定係数であり、係数のかかった項は制約条件である。評価関数（６）の第１行は、式（５）を代入して整理すると、式（６）の第３行のように整理され、変分をとると式（７）のようになる。微小変化δφT+及びδφR+は任意であるので、変分が０という条件式（８）の第１行を変形することで、問題は、式（８）の第３行に得られた、それぞれがφT及びφRに関する固有値問題である、２元連立方程式に帰着することがわかる。すなわち、演算子Ｔ+Ｗ+ＷＴＴ+Ｗ+ＷＴの固有関数φを求めて受信用あるいは送信用フォーカスデータとし、Ｔ+Ｗ+ＷＴφを送信用あるいは受信用フォーカスデータとして用いれば、フォーカス点での音響強度を一定にした条件下で、可能な限り形状Ｗに近いビームが生成される。
【００６２】
ここで、第１の実施例の式（４）では、１本の固有値方程式が得られたが、第４の実施例の式（８）では２元連立方程式である点が異なる。２元連立方程式であるため、送信用と受信用で異なるフォーカスデータを求めることができる。
【００６３】
以上のような構成によれば、送信用のフォーカスデータと、送信用とは異なる、受信用のフォーカスデータがそれぞれ求まり、送受での指向性が、サイドローブのみでなくグレーティングローブも抑制したビームを形成することが可能になる。その結果、セクタ探触子で偏向角を大きくしても高画質な撮像が可能になり、広い撮像範囲を設定できる。また、リニア探触子において、大きなオブリークをつけてスキャンした場合、たとえば台形型の広い撮像範囲を持つ場合、あるいは、コンパウンド処理を行う各場合にも高画質な撮像が可能になる。特にコンパウンド処理の場合は、加算フレーム数が増えるため、１スキャンにおけるＳＮ比の向上という直接的な効果のみでなく、加算フレーム数が増えたためにコントラストがより向上するという間接的な効果をも得られる。さらに、高周波探触子や２次元探触子で、素子ピッチが波長の半分以上になることから必然的に生じるグレーティングの影響を低減することができ、高画質な撮像が可能になる。
【００６４】
以下、図１２、図１３を用いて、本発明の第５の実施例を説明する。以下、特に断らない限り、演算子は離散、φは空間１次元、時間０次元、Ｔは空間１次元、時間０次元、Ｗは空間１次元、時間０次元の例を基本として説明する。
【００６５】
本実施例の超音波装置は、情報を格納するメモリを具備し、メモリは複数の素子ごとの受信信号を格納する受信信号格納部を備える。目標音場決定用データ取得部で取得するデータは、受信信号格納部に格納された、複数の素子ごとの受信信号であり、画像生成部は、受信信号格納部から複数の素子ごとの受信信号を読み出し、フォーカスデータ生成部が出力したフォーカスデータを用いて画像を再構成する。フォーカスデータ生成部は、フォーカスデータから音場への変換を表す演算子をＴとし、目標ビーム形状を表す関数をＷとして、演算子Ｔ-1Ｗ-1ＷＴの固有関数φに比例したフォーカスデータを出力する、あるいは、演算子Ｔ+Ｗ+ＷＴＴ+Ｗ+ＷＴの固有関数φTに比例した送信用フォーカスデータとＴ+Ｗ+ＷＴφTに比例した受信用フォーカスデータφTを出力する。
【００６６】
フォーカスデータ生成部の処理において、φは空間次元が１次元でも２次元でもよい。φの空間次元が１次元の場合は１次元探触子に相当し、Ｔは２次元行列、Ｗは２次元行列になる。φの空間次元が２次元の場合は２次元探触子に相当し、Ｔは３階テンソル、Ｗは２次元行列になる。以上のような構成によれば、探触子が１次元の場合も２次元の場合もフォーカスデータを求めることができる。
【００６７】
あるいは、フォーカスデータ生成部の処理において、φは周波数空間の次元あるいは時間次元が０次元でも１次元でもよい。φの周波数空間の次元あるいは時間次元が０次元の場合、Ｔは２次元行列、Ｗは２次元行列になる。φの周波数空間の次元あるいは時間次元が１次元の場合、Ｔは３階テンソル、Ｗは２次元行列になる。以上のような構成によれば、周波数空間の次元あるいは時間次元が０次元の場合、探触子のフォーカスデータを求めることができる。周波数空間の次元あるいは時間次元が１次元の場合、探触子のフォーカスデータを求めることができるのみでなく、素子ごとの最適なパルスを設計することができる。
【００６８】
あるいは、フォーカスデータ生成部の処理において、Ｗは空間次元が１次元でも２次元でも３次元でもよい。Ｗの空間次元が１次元の場合、Ｔは２次元行列、Ｗは２次元行列になる。Ｗの空間次元が２次元の場合、Ｔは３階テンソル、Ｗは４階テンソルになる。Ｗの空間次元が３次元の場合、Ｔは４階テンソル、Ｗは６階テンソルになる。
【００６９】
図１２は、Ｗの空間次元が２次元の場合の、目標ビーム形状とフォーカスデータ生成部の計算したフォーカスデータによって形成される音場を示した図である。１００１が目標ビーム形状、１２０２が形成される音場である。目標ビーム形状１００１によく一致して、深度方向の長さが通常より数倍長い焦域を持つ音場が得られている。このような長い焦域が得られると、広い深度領域で高画質な画像が得られる。さらに、造影剤を均質に破壊できる。さらに、従来より深い深度で、従来より強い高周波成分を発生させることができ、深部でのペネトレーションと画質を両立させた画像が得られる。
【００７０】
以上のような構成によれば、１次元空間のビーム設計も２次元空間のビーム設計も、３次元空間のビーム設計も行うことができる。
【００７１】
あるいは、フォーカスデータ生成部の処理において、Ｗは周波数空間の次元あるいは時間次元が０次元でも１次元でもよい。Ｗの周波数空間の次元あるいは時間次元が０次元の場合、Ｔは２次元行列、Ｗは２次元行列になる。Ｗの周波数空間の次元あるいは時間次元が１次元の場合、Ｔは３階テンソル、Ｗは４階テンソルになる。
【００７２】
以上のような構成によれば、時間を伴ったビーム設計が可能になる。
【００７３】
あるいは、フォーカスデータ生成部の処理において、Ｔは離散化する演算Ｇを含んでもよい（Ｔ’＝ＴＧ）。この場合、Ｔ’は２次元行列、Ｗは２次元行列になる。
【００７４】
図１３は、８個の素子からなる１次元探触子を２素子ずつ連結し、偏向角４５度の方向にメインビームを作る場合の、フォーカスデータ生成部の計算した連結パターン及びフォーカスデータ、及びそれらによって形成される音場を示した図である。
【００７５】
下式は、離散化する演算子Ｇの例である
【００７６】
【数４】
2010104584.tif 000006
【００７７】
散化する演算子Ｇは（素子数）行（素子数）列、すなわち８行８列の２次元行列で、素子に番号１から８をつけ、素子１と素子２、素子３と素子ｊを連結する場合を考えると、次式を満たす行列になる。
【００７８】
Ｇ11＝Ｇ12＝Ｇ21＝Ｇ22＝１／２，
Ｇ33＝Ｇ3j＝Ｇj3＝Ｇjj＝１／２，
．．．
【００７９】
フォーカスデータ生成部は、あり得るすべての連結パターン、すなわち8Ｃ2・6Ｃ2・4Ｃ2・2Ｃ2２通りの連結パターンに対して演算子Ｇを計算し、それぞれのＧに対してＴ’＝ＴＧを計算し、演算子Ｔ’ +Ｗ+ＷＴ’を計算し、演算子Ｔ’ +Ｗ+ＷＴ’の固有関数φnを計算し、すべてのＧに対する固有関数の中で最大の固有値をもつφnをフォーカスデータに設定する。
【００８０】
その結果出力された、連結パターンを表す行列Ｇが１３０２ａ、フォーカスデータの絶対値が１３０２ｂ、フォーカスデータの位相が１３０２ｃ、これらの連結パターン１３０２ａ、フォーカスデータ１３０２ｂ，１３０２ｃによって形成された音場が１３０２ｄである。連結がない場合にフォーカスデータ生成部が出力するフォーカスデータの絶対値を１３０３ｂ、フォーカスデータの位相を１３０３ｃ、それらを、本発明のように数値的に最適化せずに単純に考えて、フォーカスデータの位相の値が近い素子同士を連結した場合に形成される音場を１３０３ｄに、比較のために示した。なお、グラフ１３０２ａ、１３０２ｂ、１３０２ｃ、１３０３ｂ、１３０３ｃの横軸及び１３０２ａの縦軸は素子番号、グラフ１３０２ｄ、１３０３ｄの横軸は方位角、縦軸は音響強度（線形表示）である。
【００８１】
フォーカスデータの強度は、本発明の場合（１３０２ｂ）の方が単純に考えた従来例（１３０３ｂ）より大きいことから、音圧を確保できており、ペネトレーションのよい画像が得られている。また、形成される音場は、本発明の場合（１１０２ｄ）には従来例（１３０３ｄ）よりサイドローブが抑制されており、ＳＮ比のよい画像が得られている。
【００８２】
以上のような構成によれば、たとえば素子数が信号線数より格段に多く、複数の素子を１本の信号線に連結しなくてはいけない場合に、フォーカスデータのみでなく、素子のグループ分けのパターンをも目標とするビームに対して最適化することができる。
【図面の簡単な説明】
【００８３】
【図１Ａ】本発明の超音波装置の構成例を表すブロック図。
【図１Ｂ】処理の流れを説明するフローチャート。
【図２】信号処理部の処理を説明するフローチャート。
【図３】目標ビーム形状設定部の処理を説明するフローチャート。
【図４Ａ】フォーカスデータ生成部の処理を説明するフローチャート。
【図４Ｂ】関数Ｗと固有関数φの説明図。
【図５Ａ】目標とする音場を決定するためのデータを取得する処理の例を説明する概念図。
【図５Ｂ】目標ビーム形状を設定する処理の例を説明する概念図。
【図５Ｃ】画像を生成する処理の例を説明する概念図。
【図６】回避する部位の位置及び強度と目標とするビームの形状との数値的な関係の一例を説明する図
【図７】フォーカスデータ生成部の入出力の計算例を示す図。
【図８】本発明の超音波装置の効果を検証するシミュレーション画像。
【図９Ａ】リアルタイム処理で目標音場決定用データからフォーカスデータを計算し、フォーカスデータを使用する処理の流れを説明するフローチャート。
【図９Ｂ】オフライン処理で目標音場決定用データからフォーカスデータを計算し、フォーカスデータを使用する処理の流れを説明するフローチャート。
【図１０】目標ビーム形状設定部の処理の一例を示した図。
【図１１】フォーカスデータ生成部の処理の例を説明するフローチャート。
【図１２】本発明のフォーカスデータ生成部の、Ｗの空間次元が２次元の場合の、目標ビーム形状と、計算したフォーカスデータによって形成される音場を示した図。
【図１３】本発明のフォーカスデータ生成部の、離散化する演算Ｇ場合の、計算した連結パターン及びフォーカスデータ、及びそれらによって形成される音場を示した図。
【符号の説明】
【００８４】
１…探触子、２…装置本体、３…送信ビームフォーマ、４…増幅手段、５…受信ビームフォーマ、６…信号処理部、７…メモリ、８…表示部、９…入力手段、１０…制御部、６１…目標設定用データ取得部、６２…回避部検出部、６３…目標ビーム形状設定部、６４…フォーカスデータ生成部
【図１Ａ】
2010104584.tif 000007
【図１Ｂ】
2010104584.tif 000008
【図２】
2010104584.tif 000009
【図３】
2010104584.tif 000010
【図４Ａ】
2010104584.tif 000011
【図４Ｂ】
2010104584.tif 000012
【図５Ａ】
2010104584.tif 000013
【図５Ｂ】
2010104584.tif 000014
【図５Ｃ】
2010104584.tif 000015
【図６】
2010104584.tif 000016
【図７】
2010104584.tif 000017
【図８】
2010104584.tif 000018
【図９Ａ】
2010104584.tif 000019
【図９Ｂ】
2010104584.tif 000020
【図１０】
2010104584.tif 000021
【図１１】
2010104584.tif 000022
【図１２】
2010104584.tif 000023
【図１３】
2010104584.tif 000024