特許 7422099 超音波撮像装置、信号処理装置、および、信号処理方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2024-01-17

(45)【発行日】2024-01-25

(54)【発明の名称】超音波撮像装置、信号処理装置、および、信号処理方法

(51)【国際特許分類】

   A61B 8/14 20060101AFI20240118BHJP

【ＦＩ】

A61B8/14

【請求項の数】 14

(21)【出願番号】P 2021007305

(22)【出願日】2021-01-20

(65)【公開番号】P2022111699

(43)【公開日】2022-08-01

【審査請求日】2023-07-24

【早期審査対象出願】

(73)【特許権者】

【識別番号】320011683

【氏名又は名称】富士フイルムヘルスケア株式会社

(74)【代理人】

【識別番号】110000888

【氏名又は名称】弁理士法人山王坂特許事務所

(72)【発明者】

【氏名】池田 貞一郎

(72)【発明者】

【氏名】網野 和宏

(72)【発明者】

【氏名】久津 将則

【審査官】佐野 浩樹

(56)【参考文献】

【文献】国際公開第２０１６／１２５５０９（ＷＯ，Ａ１）

【文献】特開２０１８－０８２８３５（ＪＰ，Ａ）

【文献】特開２００８－１６１２０８（ＪＰ，Ａ）

【文献】特開２０１５－２２６７６２（ＪＰ，Ａ）

【文献】特開２０１６－１８５２１２（ＪＰ，Ａ）

【文献】国際公開第２０１５／１６６８６９（ＷＯ，Ａ１）

【文献】特開２０２０－０１０８７４（ＪＰ，Ａ）

【文献】特開２０１１－０８７７１０（ＪＰ，Ａ）

【文献】米国特許出願公開第２０１６／０１５１０４２（ＵＳ，Ａ１）

【文献】特開２０１０－１６２１８５（ＪＰ，Ａ）

【文献】特開平０２－０７５２８３（ＪＰ，Ａ）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ａ６１Ｂ ８／００ － ８／１５

(57)【特許請求の範囲】

【請求項1】

接続されている超音波素子アレイの複数の超音波素子からそれぞれ、所定の送信焦点に集束するように位相を遅延させた超音波を被検体に対して送信させる送信ビームフォーマと、
前記超音波の送信を受けた前記被検体から前記超音波素子アレイに戻った超音波を、前記超音波素子アレイの複数の超音波素子がそれぞれ受信した受信信号を受け取って、前記被検体の深度に応じた遅延時間によりそれぞれ遅延させて加算し、整相信号を生成する受信ビームフォーマと、
前記遅延時間が前記超音波素子ごとに格納された遅延時間格納部と、
前記遅延時間を演算により生成して前記遅延時間格納部に格納する遅延時間演算部と、
合成部と、
前記被検体の撮像対象の部位または臓器の選択を操作者から受け付けるコンソールとを含み、
前記遅延時間演算部は、被検体の複数の部位ごと、または、複数の臓器ごとに予め求めておいた、遅延時間の数を増加させることが解像度の向上に効果的となる深度範囲に基づいて、前記コンソールを介して操作者から受け付けた撮像対象の部位または臓器の、遅延時間の数を増加させる深度範囲を求め、当該深度範囲については、他の深度範囲より生成する遅延時間の数を増加させ、
前記受信ビームフォーマは、前記遅延時間格納部に前記遅延時間が２以上格納されている深度範囲では、同一の前記受信信号を２以上の前記遅延時間によってそれぞれ遅延させることにより、２以上の整相信号を生成し、
前記合成部は、前記整相信号が２以上生成された深度範囲について、前記２以上の整相信号を合成することを特徴とする超音波撮像装置。

【請求項2】

請求項１に記載の超音波撮像装置であって、操作者から送受信パラメータを受け付ける受付部をさらに有し、
前記遅延時間演算部は、前記遅延時間を前記送受信パラメータに基づいて演算により求めることを特徴とする超音波撮像装置。

【請求項3】

請求項２に記載の超音波撮像装置であって、前記遅延時間演算部は、送受信パラメータと、深度範囲ごとの遅延時間の数および値との関係を予め定めておいたテーブルを参照して、前記受付部が受け付けた前記送受信パラメータに対応する前記遅延時間の数および値を求めることを特徴とする超音波撮像装置。

【請求項4】

請求項１に記載の超音波撮像装置であって、
前記遅延時間演算部および前記遅延時間格納部は、前記受信ビームフォーマ内に配置され、
前記遅延時間演算部は、前記受信ビームフォーマによる前記受信信号の遅延加算処理に対応させて、前記深度範囲によって個数の異なる遅延時間をリアルタイムに算出し、前記遅延時間格納部内に格納されている遅延時間を順次上書きしていくことを特徴とする超音波撮像装置。

【請求項5】

請求項１に記載の超音波撮像装置であって、前記遅延時間演算部は、予め定めておいた基準となる遅延時間に対して、予め定めた調整係数α１，α２、α３・・αｎをそれぞれ掛けることにより、深度範囲ごとにｎ個の遅延時間の値を算出することを特徴とする超音波撮像装置。

【請求項6】

請求項１に記載の超音波撮像装置であって、前記遅延時間演算部は、前記超音波素子ごとに発せられ、前記超音波素子ごとに受信される超音波の伝搬経路を、予め定めた音波伝搬シミュレーションに基づいて求め、求めた前記伝搬経路を用いて、前記深度範囲によって個数の異なる遅延時間を算出することを特徴とする超音波撮像装置。

【請求項7】

請求項１に記載の超音波撮像装置であって、前記遅延時間演算部は、機械学習モデルを備え、前記機械学習モデルに、本撮像の前に撮像した超音波画像またはその受信信号に基づき、適切な遅延時間の数および／または深さに応じた遅延時間の値を深度ごとに算出させることを特徴とする超音波撮像装置。

【請求項8】

請求項１に記載の超音波撮像装置であって、前記受信ビームフォーマは、前記遅延時間格納部に格納されている前記遅延時間の最大個数と同数の遅延加算回路を備え、それぞれの遅延加算回路に同一の受信信号が入力され、異なる遅延時間によって遅延加算処理されることを特徴とする超音波撮像装置。

【請求項9】

請求項１に記載の超音波撮像装置であって、前記受信ビームフォーマは、前記遅延時間格納部に格納されている前記遅延時間の個数よりも少ない数の遅延加算回路と、受信信号記憶部と、整相信号記憶部を備え、
前記受信信号記憶部に格納された同一の受信信号を、前記遅延時間格納部内の遅延時間で遅延加算処理し、整相信号記憶部に格納する処理を、前記遅延時間格納部内の複数の遅延時間について順番に行うことを特徴とする超音波撮像装置。

【請求項10】

請求項１に記載の超音波撮像装置であって、前記合成部は、前記２以上の整相信号を、深さに応じて設定された重みにより、重みづけしてから加算することを特徴とする超音波撮像装置。

【請求項11】

請求項２に記載の超音波撮像装置であって、前記送受信パラメータは、開口径、周波数、および、送信焦点位置の少なくとも１つを含むことを特徴とする超音波撮像装置。

【請求項12】

請求項１に記載の超音波撮像装置であって、前記超音波素子アレイが搭載された超音波探触子をさらに含み、
前記受信ビームフォーマは、前記超音波探触子に搭載されていることを特徴とする超音波撮像装置。

【請求項13】

超音波の送信を受けた被検体から超音波素子アレイに到達した超音波を複数の超音波素子がそれぞれ受信した受信信号を受け取って処理する信号処理装置であって、
前記被検体の深度に応じた遅延時間により前記受信信号をそれぞれ遅延させて加算し、整相信号を生成する受信ビームフォーマと、
前記遅延時間が前記超音波素子ごとに格納された遅延時間格納部と、
前記遅延時間を演算により生成して前記遅延時間格納部に格納する遅延時間演算部と、
合成部と、
前記被検体の撮像対象の部位または臓器の選択を操作者から受け付けるコンソールとを含み、
遅延時間演算部は、被検体の撮像対象の複数の部位ごと、または、複数の臓器ごとに予め求めておいた、遅延時間の数を増加させることが解像度の向上に効果的となる深度範囲に基づいて、前記コンソールを介して操作者から受け付けた撮像対象の部位または臓器の、遅延時間の数を増加させる深度範囲を求め、当該深度範囲については、他の深度範囲より生成する遅延時間の数を増加させ、
前記受信ビームフォーマは、前記遅延時間格納部に前記遅延時間が２以上格納されている深度範囲では、同一の前記受信信号を２以上の前記遅延時間によってそれぞれ遅延させることにより、２以上の整相信号を生成し、
前記合成部は、前記整相信号が２以上生成された深度範囲について、前記２以上の整相信号を合成することを特徴とする信号処理装置。

【請求項14】

超音波の送信を受けた被検体から超音波素子アレイに到達した超音波を複数の超音波素子がそれぞれ受信した受信信号を受け取って処理する信号処理方法であって、
前記被検体の撮像対象の部位または臓器の選択を操作者から受け付け、
被検体の撮像対象の複数の部位ごと、または、複数の臓器ごとに予め求めておいた、遅延時間の数を増加させることが解像度の向上に効果的となる深度範囲に基づいて、操作者から受け付けた撮像対象の部位または臓器の、遅延時間の数を増加させる深度範囲を求め、
当該遅延時間の数を増加させる前記深度範囲および他の深度範囲について遅延時間を生成し、その際、前記遅延時間の数を増加させる前記深度範囲については、前記他の深度範囲よりも、多くの数の遅延時間を生成し、
生成された前記遅延時間によって、前記受信信号をそれぞれ遅延させて加算することにより、前記深度範囲によって１つ、または、２以上の整相信号を生成し、
前記整相信号が２以上生成された深度範囲について、前記２以上の整相信号を合成することを特徴とする信号処理方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、超音波を用いて被検体内の画像を撮像する超音波撮像技術に関する。

【背景技術】

【0002】

超音波撮像技術とは、超音波（聞くことを意図しない音波、一般的には２０ｋＨｚ以上の高周波数の音波）を用いて人体をはじめとする被検体の内部を非侵襲的に画像化する技術である。

【0003】

超音波撮像装置による超音波の送受信は、有限の開口径を持つ超音波素子アレイによって行われるため、開口部のエッジによる超音波の回折の影響を受け、方位角方向の分解能を向上させることが難しい。そのため、適応ビームフォーマや、開口合成などの新しい整相方式が提案されている。

【0004】

特許文献１には、集束型送信を行う超音波撮像技術において、仮想音源法を改良した方法を用いて、開口合成を行う技術が開示されている。具体的には、超音波ビームのエネルギが焦点に収束する領域（特許文献１の図２の領域A）では、焦点を仮想音源とみなして開口合成を行い、その周辺の超音波エネルギが拡散する領域（領域B,C）では探触子の端部から球面波が放射されたとみなして開口合成を行う。

【0005】

一方、特許文献２には、２種類以上の遅延時間によって受信信号をそれぞれ遅延加算する２以上の遅延加算部を備えた超音波撮像装置が開示されている。第１の遅延加算部は、超音波素子から送信された送信ビーム（干渉波）から生じた受信信号を整相する第１の遅延時間によって遅延加算する。第２の遅延加算部は、送信ビームとは位相の異なる回折波（球面波）から生じた受信信号を整相する第２の遅延時間によって受信信号を遅延加算する。第１および第２遅延加算部がそれぞれ遅延加算した後の信号は合成される。

【先行技術文献】

【特許文献】

【0006】

【文献】特開平１０－２７７０４２号公報

【文献】国際公開第２０１６／１２５５０９号

【発明の概要】

【発明が解決しようとする課題】

【0007】

特許文献１の技術は、送信間で開口合成を行う技術であり、１回の送信のみでは、高精度な画像を得ることはできない。

【0008】

特許文献２の技術は、送信ビーム（干渉波）と、それ以外の回折波（球面波）で異なる遅延時間を用いることを開示しているが、遅延加算部の数は、２個または３個で固定的であり、遅延加算部ごとに遅延時間が設定されるため、遅延時間の種類も２つまたは３つである。

【0009】

発明者らの研究によると、被検体内には、超音波素子アレイの複数の超音波素子からそれぞれ送信された超音波（球面波）が伝搬しており、これらは相互に干渉して送信ビーム（干渉波）を形成する他に、球面波のまま各方向に伝搬したり、一部のみが干渉しあったり、交差したりしながら深さ方向に複雑に伝搬していく。そのため、特許文献２のように、２個または３個の遅延加算部では、それらの波によってそれぞれ生じる受信信号の一部しか整相加算することができない。一方、遅延加算部の数を超音波を送信した素子数と同程度もしくはそれ以上に増やせば、複雑な波によって生じた受信信号をそれぞれ整相加算することが可能になるが、回路規模が増大してしまう。

【0010】

発明の目的は、複数の超音波素子からそれぞれ送信され、被検体の深さ方向に複雑に伝搬していく超音波によって生じた受信信号を、効率よく整相加算して、より高解像度な画像を生成することにある。

【課題を解決するための手段】

【0011】

本発明の超音波撮像装置は、接続されている超音波素子アレイの複数の超音波素子からそれぞれ、所定の送信焦点に集束するように位相を遅延させた超音波を被検体に対して送信させる送信ビームフォーマと、超音波の送信を受けた被検体から超音波素子アレイに戻った超音波を、超音波素子アレイの複数の超音波素子がそれぞれ受信した受信信号を受け取って、被検体の深度に応じた遅延時間によりそれぞれ遅延させて加算し、整相信号を生成する受信ビームフォーマと、遅延時間が超音波素子ごとに格納された遅延時間格納部と、合成部とを含む。遅延時間格納部に格納されている遅延時間の個数は、被検体の深度範囲によって異なる。受信ビームフォーマは、遅延時間格納部に格納されている遅延時間が２以上格納されている深度範囲では、同一の受信信号を２以上の遅延時間によってそれぞれ遅延させることにより、２以上の整相信号を生成する。合成部は、整相信号が２以上生成された深度範囲について、２以上の整相信号を合成する。

【発明の効果】

【0012】

本発明によれば、送信ビームだけでなく、複数の超音波素子からそれぞれ送信され、被検体の深さ方向に複雑に伝搬していく球面波によって生じた受信信号を、効率よく整相加算して、より高解像度な画像を生成することができる。

【図面の簡単な説明】

【0013】

【図1】第１実施形態の超音波撮像装置の構成を示すブロック図。

【図2】送信ビーム３１の形状および波面と、球面波とを示す説明図。

【図3】第１実施形態の深度範囲ごとの遅延曲線の一例を示すグラフ。

【図4】第１実施形態の深度範囲ごとの遅延曲線の一例を示すグラフ。

【図5】第１実施形態の深度範囲ごとの遅延曲線の一例を示すグラフ。

【図6】第２実施形態の超音波撮像装置の構成を示すブロック図。

【図7】第３実施形態の超音波撮像装置の構成を示すブロック図。

【図8】第３実施形態の超音波撮像装置の撮影時の各部の動作を示すフロー図。

【図9】（ａ）は、比較例の超音波撮像装置で撮像した超音波画像、（ｂ）は、本実施形態の超音波撮像装置で撮像した超音波画像、（ｃ）は、図（ａ）および（ｂ）の超音波画像の深度方向の輝度のプロファイルを示すグラフ、（ｄ）および（ｅ）は、図（ａ）および（ｂ）の超音波画像の方位方向の輝度のプロファイルを示すグラフである。

【発明を実施するための形態】

【0014】

本発明の一実施形態の超音波撮像装置について説明する。

【0015】

＜＜第１実施形態＞＞
まず、第１実施形態の超音波撮像装置１００について図１～図５を用いて説明する。図１は、超音波診断装置の構成を示す図であり、図２は、超音波素子アレイから送信された超音波の波面の例を示す図であり、図３～図５は、遅延時間格納部に格納されている深度ごとの遅延時間を示すグラフである。

【0016】

図１に示すように、本実施形態の超音波撮像装置１００は、超音波撮像装置本体１０２および超音波撮像装置１０２に接続された超音波素子アレイ１０１を備えた超音波探触子１１６で構成されている。超音波素子アレイ１０１は、複数の超音波素子をアレイ状に配列した構成である。超音波素子アレイ１０１において、送信に用いる複数の超音波素子（チャンネル）を送信チャンネルと呼び、受信に用いる複数の超音波チャンネルを受信チャンネルと呼ぶ。

【0017】

超音波撮像装置本体１０２は、送信ビームフォーマ１０４と、受信ビームフォーマ１３と、遅延時間格納部１７と、合成部１５と、遅延時間演算部１８と、合成重み演算部１９と、画像処理部１０９と、送受分離回路１０７と、アナログ／デジタルコンバータ（ＡＤＣ）１１と、全体を制御する制御部１１１とを備えて構成される。

【0018】

送信ビームフォーマ１０４は、超音波素子アレイ１０１の複数の送信チャンネル１０５に送信信号を出力し、複数の送信チャンネル１０５からそれぞれ、所定の送信焦点３０に集束するように位相を遅延させた超音波を被検体９０に対して送信させる（図２参照）。

【0019】

これにより、図２に示すように、被検体９０内には、送信チャンネル１０５からそれぞれ送信された超音波（球面波）３０１～３０６等が伝搬し、これらは相互に干渉して送信ビーム（干渉波）３１０を形成する他に、一部は球面波３０１～３０６等のまま各方向に伝搬し、さらにそれらの一部のみが干渉しあったりしながら深さ方向に複雑に伝搬していく。

【0020】

超音波３０１～３０６、３１０等の送信を受けた被検体９０から超音波素子アレイ１０１に戻った超音波は、超音波素子アレイ１０１の複数の受信チャンネル１０６によってそれぞれ受信され、受信信号に変換される。

【0021】

受信ビームフォーマ１３は、送受分離回路１０７およびアナログ／デジタルコンバータを介して、複数の受信チャンネル１０６から受信信号を受け取って、被検体９０の深度に応じた遅延時間によりそれぞれ遅延させて加算し、受信走査線３６上の複数の撮像対象点からの受信音波に対して位相（受信焦点）の合った整相信号を生成する。

【0022】

遅延時間格納部１７には、受信ビームフォーマ１３が受信信号を遅延させる際に用いる遅延時間が、超音波素子（受信チャンネル１０６）ごとに格納されている。遅延時間は、受信走査線３６上の複数の撮像対象点に位相（受信焦点）を合わせるため、例えば図３の破線の曲線３５０に示すように深度によって変化している。ここでは、深度ごとの遅延時間の変化を遅延曲線とも呼ぶ。

【0023】

本実施形態では、送信ビーム３１０により生じた受信信号のみならず、球面波３０１～３０６等のまま伝搬する波や、それらの一部が干渉した干渉波等から生じた受信信号からも整相信号を生成するために、図３～図５に例示したように複数の遅延時間（遅延曲線）が遅延時間格納部１７に格納されている。

【0024】

これにより、受信ビームフォーマ１３は、同一の受信信号を、それぞれの遅延時間で遅延させて加算することにより、同一の受信走査線３６について複数の整相信号を生成することができる。

【0025】

合成部１５は、生成された複数の整相信号を同一深度ごとに加算することにより、送信ビーム３１０により生じた受信信号から生成した整相信号のみならず、球面波３０１～３０６等の他の波の一部または全部により生じた受信信号から生成した整相信号を合成した整相信号を同一の受信走査線について得ることができる。よって、送信ビーム３１０のみについて得た整相信号よりも、高解像度の整相信号を得ることができる。

【0026】

このとき、すべての球面波３０１～３０６等や、それらの一部の干渉波等から整相信号を得るためには、遅延時間の個数（遅延曲線の本数）をそれらの球面波３０１～３０６等や一部の干渉波等の数の分だけ用意する必要があり、送信チャンネルと同数程度の遅延時間の個数が必要になる。その場合、受信ビームフォーマ１３の演算量も増加する。

【0027】

そこで、本実施形態では、被検体９０の深度範囲によって遅延時間の個数（遅延曲線の本数）を異ならせることにより、遅延時間の個数を増加させる深度範囲を制限し、演算量を抑制する。

【0028】

例えば、遅延時間格納部１７に格納する遅延時間の数と値を深度範囲ごとに演算により算出する遅延時間演算部１８をさらに備える構成とすることができる。

【0029】

遅延時間演算部１８は、例えば、コンソール（受付部）１１０を介して操作者が設定した撮像条件パラメータおよび／または撮像パラメータに基づいて設定される送受信パラメータおよび／または接続されている超音波探触子１１６の種類に基づいて、深度範囲ごとの遅延時間の数と値を、演算により求める。撮像条件パラメータとはコンソール１１０上で操作者に対して明示的に表示される入力条件であり、例えば撮像モード（カラーフロー撮像、ドプラ撮像、非線形撮像、開口合成、周波数・空間コンパウンド、コントラスト・エッジ強調撮像など各種撮像モード）などいくつかの選択肢の中から操作者が選ぶものや、音波強度、高・中・低などの周波数帯域、深さ毎の撮像ゲイン、画像拡大率などツマミやボタンなどで段階もしくは連続的に変化した値として操作者が入力可能な条件などがある。送受信パラメータは前記の撮像条件パラメータに基づいてあらかじめ装置内に用意されたテーブルや数式演算を介して変換された超音波送受信に関するパラメータである。送受信パラメータとしては、例えば、送信開口径、受信開口径、周波数（中心周波数、周波数帯域）、送信フォーカス位置、送信パルス波の形状（波数・振幅）などを用いる。

【0030】

具体的には、撮像条件パラメータおよび／または送受信パラメータおよび／または超音波探触子１１６の種類と、深度範囲ごとの遅延時間の数および値との関係を定めたテーブルを予め定めて遅延時間演算部１８内のメモリに格納しておき、遅延時間演算部１８はこのテーブルを参照することにより、遅延時間の数と値を求める構成とすることができる。また、テーブルに限られるものではなく、予め定めておいた数式に送信パラメータの値を代入して、深度範囲ごとの遅延時間の数と値を算出してもよい。

【0031】

遅延時間演算部１８は、図３に示すように予め定めておいた基準となる遅延時間に対して、調整係数α１，α２、α３・・αｎをそれぞれ掛ける等することにより、深度範囲ごとにｎ個の遅延時間の値（遅延曲線）を生成する構成としてもよい。α１，α２、α３・・αｎは、深度（ｄ）の関数である。

【0032】

遅延時間演算部１８は、遅延線の内部計算で遅延時間の値を深度ごとに算出する構成としてよい。遅延線の内部計算とは、深度ごとの遅延時間を信号処理中にリアルタイムに演算する方式である。この場合、遅延時間演算部１８および遅延時間格納部１７は受信ビームフォーマ１３の中に位置する構成となる。遅延線の内部計算を行う場合、遅延時間演算部１８はＡＤＣ１１から逐次入力される受信信号に対してサンプルごともしくは受信走査線ごと、もしくは音波送信ごとなどごく短時間の時間間隔で、遅延時間を計算し、同じく受信ビームフォーマ１３の中に位置する遅延時間格納部１７に即座に格納する。また遅延時間格納部１７中の遅延時間の値は、受信ビームフォーマ１３での遅延加算処理が遂行された直後に、次のサンプルもしくは受信走査線もしくは音波送信について使用される遅延時間の値にリフレッシュ（上書き）される。この方式によれば一度に遅延時間格納部１７に格納する遅延時間の値のデータ量を大幅に低減することができる。よって、遅延時間格納部１７として大容量のメモリを必要とせず、例えばＦＰＧＡやＡＳＩＣ内部のトランジエントメモリによって遅延時間格納部１７を代替することが可能となる。

【0033】

さらに、遅延時間演算部１８は、予め定めた音波伝搬シミュレーションに基づいて、演算する遅延時間を深度ごとに算出する構成としてもよい。一般的に遅延曲線は音波伝搬を簡易的にモデル化した数式で表されるが、送受信パラメータに基づいた音波伝搬シミュレーションを遅延時間演算部１８が直接実行することでより、正確な遅延時間を計算することができる。具体的には、送信開口、受信開口、周波数、送信フォーカス位置、送信パルス波の形状などを入力として、送信１チャンネルごとに発する超音波と受信１チャンネルごとに受信する超音波の音波伝搬の様態を１次元または２次元もしくは３次元的に微分方程式、差分式、グリーン関数などを用いた伝搬シミュレーションによって演算することにより、音波の伝搬経路を正確に求め、伝搬経路を用いて遅延曲線の遅延時間を算出する。音波伝搬シミュレーションに基づいて演算された正確な遅延時間を受信ビームフォーミングに用いることで、さらに高画質な超音波撮像が可能となる。音波伝搬シミュレーションは演算規模が大きいため、受信ビームフォーマ１３がＦＰＧＡやＡＳＩＣなどのロジックデバイスを利用した構造の場合、音波伝搬シミュレーションを実現するためデバイス規模を大きくすることが望ましい。受信ビームフォーマのデバイスとしてＣＰＵ、ＧＰＵを利用した構造の場合には、ソフトウェア演算の形態により音波伝搬シミュレーションを好適に実現することができる。

【0034】

また、例えば、予め求めた撮像対象の部位（腹部、循環器、胸部、脚部、血管部、消化器、妊婦検診等）や臓器（肝臓、心臓、腎臓、膵臓、胆嚢、卵巣、頸動脈、甲状腺等）ごとに、遅延時間の数を増加させることが解像度の向上に効果的となるような深度範囲を予め求めておき、遅延時間演算部１８は、その深度範囲の遅延時間の数を増加させる構成にすることができる。部位や臓器の選択は、制御部１１１に接続されているコンソール１１０を介して、操作者から受け付ける。

【0035】

また、遅延時間演算部１８は、コンソール１１０を介して操作者から、高解像度の撮像を望む深度範囲を受け付け、その深度範囲の遅延時間の数を増加させてもよい。もしくは、予め所定の深度範囲において遅延時間の数を増加させておいた遅延時間のパターンを用意しておき、操作者が高解像度の撮像を望む深度範囲に応じて用いることもできる。

【0036】

さらには、遅延時間演算部１８は、機械学習モデルを内蔵し、撮像条件パラメータ、送受信パラメータ、探触子の種類など超音波撮像に関係する各種パラメータや、本撮像の前に準備として撮像した超音波画像またはその受信信号を、機械学習モデルに入力して、適切な遅延時間の数、および／または深さに応じた遅延時間の値深度ごとに算出させることも可能である。機械学習モデルは、撮像条件や、撮像した画像または受信信号や、その撮像に用いた遅延時間の数を入力データとし、遅延時間の数を増加させて得た高精度の画像やその時の遅延時間、遅延時間の数を正解データとして、予め学習させたものを用いる。また、撮像の精度や分解能などは臓器の種類によっても異なるため、前期の機械学習モデルは、肝臓、腎臓、血管、乳腺などことなる臓器によって別々に学習されたものであってもよい。

【0037】

遅延曲線の具体例を図３から図５を用いて説明する。図３の遅延曲線は、送信焦点３０よりも深い深度範囲３７２には２本の遅延曲線３５０，３５１が用意され、送信焦点３０よりも浅い深度範囲３７１には、１本の遅延曲線３５０のみが用意されている例である。

【0038】

図４の遅延曲線は、送信焦点３０よりも深い深度範囲３７２には５本の遅延曲線３５０～３５４が用意され、送信焦点３０よりも浅い深度範囲３７１には、１本の遅延曲線３５０のみが用意されている例である。

【0039】

図５の遅延曲線は、送信焦点３０よりも深い深度範囲３７２には４本の遅延曲線３５０～３５４が用意され、送信焦点３０よりも浅い深度範囲３７１には、６本の遅延曲線３５０、３６１～３６５が用意されている例である。

【0040】

なお、図３から図５は、送信焦点３０を境に２つの深度範囲３７１，３７２で異なる本数の遅延曲線を設定しているが、本実施形態の深度範囲は、この２つの範囲に限定されない。任意の深度範囲を任意の数だけ設定して、その深度範囲ごとに所望の数の遅延時間（所望の本数の遅延曲線）を設定することが可能である。

【0041】

制御部１１１は、送信チャンネル１０５を移動させながら、画像生成に必要な数の受信走査線３６の整相信号が得られるまで送受信を繰り返すように各部を制御する。なお、この制御方法はリニア（線形）走査、コンベック走査などのスキャン（走査）方式の場合である。セクタ（フェーズドアレイ）型走査の場合には、送信と受信の口径は同じであるが、受信走査線３６を角度方向に傾けることで２次元平面上に複数の送信および受信走査線を設定し、その方向に沿って扇形領域の撮像を行うような形態をとる。例えば中心を探触子口径の中心とした±４５°や±６０°の扇形形状の中に５０－１３００本程度の走査線が用意されるような形態がある。この場合でも制御部１１１は、送信チャンネル１０５の代わりに送信の角度方向を移動させながら、画像生成に必要な数の受信走査線３６の整相信号が得られるまで送受信を繰り返すように各部を制御する。

【0042】

画像処理部１０９は、画像生成に必要な数の整相信号または整相信号からサンプル／ピクセルごとの信号強度・輝度値に変換された画像信号並べることにより、画像を生成し、接続されている画像表示部１０３に表示させる。

【0043】

このように、本実施形態では、深度に応じて遅延時間の数（遅延曲線の本数）を異ならせることにより、送信ビーム３１０により生じた受信信号のみならず、他の波（球面波３０１～３０６等）の一部または全部から生じた受信信号からも整相信号を得ることができるため、それらを深さごとに合成することにより、高解像度の整相信号を、演算量を抑制しながら生成することができる。

【0044】

また、本実施形態では少ない超音波送信において効率的に高解像度の整相信号を得ることができるため、送信開口合成方式、空間コンパウンド方式、符号化送受信方式、マルチビーム送信方式といった複数の送信ビームを使用した高解像度方式と比べて、大きく演算量を抑制することが可能である。またこれらの方式と異なり、１回の送信で高解像度の整相信号を得ることができるため、送信間開口合成などの複数送信を利用した方式で得た画像と比較して、時間分解能が高く、また運動体の運動に伴う画像のブレ（動き・体動アーチファクト）を伴わない。よって、本実施形態の超音波撮像装置は、時間分解能が要求される高速の運動体（心臓の拍動、心臓弁）や気泡の高速振動を利用した造影剤撮像等にも適している。

【0045】

ただし、本実施形態の超音波撮像装置は、１回の送信で１本の受信走査線について整相信号を得る構成に限定されるものではなく、１回の送信で複数本の受信走査線についてそれぞれ整相信号を生成することも可能である。これにより、１枚の画像を生成するのに必要な送信回数を低減して、高速撮像を行うことが可能である。また、時間分解能が要求されない撮像部位については、必要に応じて、送信開口合成を行ってもよい。また、また本実施形態は送信ビームフォーマ１３という超音波装置の信号処理最上流での処理であるため、送信開口合成のみならず他の超音波撮像方式、例えば非線形（ハーモニック）撮像、ドプラ撮像、カラーフロー撮像、コヒーレンス撮像、適応ビームフォーミングを利用した撮像など他の超音波撮像方式と組み合わせて用いることができる。

【0046】

なお、合成部１５は、受信ビームフォーマ１３が生成した深度範囲ごとに数の異なる整相信号を合成する際に、重みづけした後加算する構成としてもよい。その場合、合成重み演算部１９は、コンソール１１０から操作者が設定した撮像条件や、遅延時間の深度範囲ごとの数のパターン等に応じて、適切な重みを深さに応じて生成してもよい。

【0047】

また、受信ビームフォーマ１３は、複数の受信チャンネル１０６から受け取った受信信号を、複数の遅延時間でそれぞれ遅延させた後加算するために、複数の遅延加算回路を並列に備えた構成であってもよい。また、受信ビームフォーマ１３は、同一の受信信号を、遅延時間の個数よりも少ない数の遅延加算回路により、時系列に並んだ（サンプル点ごとや、深さに応じたブロックごとや、送信ごとなど）順番の遅延時間によって順番に遅延加算処理して整相信号を生成し、生成した整相信号を内蔵するメモリに順次格納していく構成であってもよい。後者の場合、遅延加算回路が最低一つでよいため、回路規模が小さくて済む。このため、超音波探触子の１１６の内部に受信ビームフォーマ１３を搭載することも可能になる。また、後者の場合、時分割処理により、複数の遅延時間で同一の受信信号を並列に遅延加算処理していくことも可能である。

【0048】

以下、第２、第３実施形態では、本実施形態の超音波撮像装置の構成及び動作を具体的に説明する。

【0049】

＜＜第２実施形態＞＞
第２実施形態の超音波撮像装置の受信ビームフォーマ１３の構成について具体的に説明する。第２実施形態では、受信ビームフォーマ１３は、並列に配置された２つの遅延加算回路１３－１、１３－２により構成されている。

【0050】

図６に示すように、受信ビームフォーマ１３が、並列な第１遅延加算回路１３－１と第２遅延加算回路１３－２を含む。第１および第２遅延加算回路１３－１、１３－２は、遅延回路と、加算回路とを含み、複数の受信チャンネル１０６からそれぞれ受け取った同一の受信信号（同一送信の超音波による撮像対象からの反射波を複数の描く受信チャンネル１０６がそれぞれ受け取った受信信号）を、受信チャンネル１０６ごとに遅延させた後、まとめて加算する。第１遅延加算回路１３－１の遅延回路は、受信信号を例えば図３の第１遅延時間３５０により遅延させた後加算することにより、第１整相信号を生成する。第１遅延時間３５０は、送信ビーム３１０を被検体９０が反射した反射波の受信信号を整相するように設定されている。第２遅延加算回路１３－２は、受信信号を図３の第２遅延時間３５１により遅延させることにより、第２整相信号を生成する。第２遅延時間は、送信ビーム３１０以外の波（例えば球面波３０１）の被検体９０が反射した反射波の受信信号を整相するように設定されている。

【0051】

図３から明らかなように、第１遅延時間３５０は、深さ方向の全範囲について設定されているため、第１遅延加算回路１３－１は、受信走査線３６のすべての深度範囲について第１整相信号を生成する。一方、第２遅延時間３５１は、送信焦点３０よりも深い深度範囲にのみ値が設定されているため、第２遅延加算回路１３－２は、受信走査線３６の送信焦点３０よりも深い深度範囲にのみ第２整相信号を生成する。

【0052】

合成部１５は、第１遅延加算回路１３－１が遅延により生成した第１整相信号と、第２遅延加算回路１３－２が遅延により生成した第２整相信号とを受け取って同一の深度の信号を合成重み演算部１９の算出した深度ごとの重みを掛けて重み付けした後加算する。

【0053】

加算後の整相信号を用いて画像を生成することにより、送信ビーム３１０のみならず、球面波（例えば球面波３０１）の情報を用いて画像を生成することができる。しかも、第２遅延時間３５１は、送信焦点３０よりも深い深度範囲にのみ設定されているため、第２遅延加算回路１３－２の演算量が低減され、受信ビームフォーマ１３全体の演算量を低減することができる。

【0054】

なお、第２実施形態では、図３の最大２つの遅延時間３５０，３５１を用いる構成であるため、受信ビームフォーマ１３は、２つの遅延加算回路１３－１、１３－２を備えていれば足りるが、図４や図５のように最大の遅延時間の数が５つまたは６つある場合には、遅延時間の最大数と同数以上の遅延加算回路を受信ビームフォーマ１３に配置する必要がある。

【0055】

超音波撮像装置の上述した以外の構成は、第１の実施形態の実施形態と同じであるので説明を省略する。

【0056】

＜＜第３実施形態＞＞
第３実施形態として、受信ビームフォーマ１３の別の構成について具体的に説明する。第３実施形態の受信ビームフォーマ１３は、遅延時間の最大数よりも少ない数の遅延加算回路を備え、時分割処理によりすべての遅延時間についての整相信号を生成する。

【0057】

具体的には、図７のように、受信ビームフォーマ１３は、１つの遅延加算回路１３－４と、受信信号記憶部１３－３と、整相信号記憶部１３－５とを備えて構成される。受信信号記憶部１３－３は、超音波素子アレイ１０１の複数の受信チャンネル１０６からそれぞれ受け取った受信信号を記憶する。整相信号記憶部１３－５は、遅延加算回路１３－４が生成した整相信号を記憶する。整相信号記憶部１３－５には、第１メモリ１３－５１～第ｎメモリ１３－５ｎまでｎ個のメモリまたは記憶領域が備えられており、遅延時間格納部１７に格納されている遅延曲線ごとに算出された整相信号は、それぞれ別々のメモリまたは記憶領域に格納される。よって、第ｎメモリ１３－５ｎのｎの数は、整相信号の生成に用いる遅延時間（遅延曲線）の最大値以上用意されている。

【0058】

遅延加算回路１３－４は、時分割処理によって、受信信号記憶部１３－３に格納した受信信号を、最大ｎ本の遅延曲線について、順に整相信号を生成し、第１メモリ１３－５１～第ｎメモリ１３－５ｎに格納していく。

【0059】

つぎに、本実施形態の超音波撮像装置の撮像時の各部の動作について、図８を用いて説明する。

【0060】

（ステップ１３１）
まず、制御部１１１は、コンソール１１０を介して撮像条件パラメータおよび／または撮像パラメータに基づいて設定される送受信パラメータおよび／または接続されている超音波探触子１１６の種類を受け付ける。送受信パラメータは、例えば、送信開口径、受信開口系、周波数（中心周波数、周波数帯域）、送信フォーカス位置、送信パルス波の形状（波数・振幅）等である。

【0061】

（ステップ１３２）
制御部１１１は、ステップ１３１で受け付けた条件等に基づいて、送信ビーム３１の形状を計算により求める。

【0062】

（ステップ１３３）
遅延時間演算部１８は、ステップ１３２で算出した送信ビーム３１の形状等を用いて、受信走査線３６の位置を設定し、各受信走査線３６上に複数の撮像対象点（受信焦点）を設定する。遅延時間演算部１８には、内蔵されているメモリに予め、送受信パラメータおよび超音波探触子１１６の種類と、深度範囲ごとの各受信チャンネル１０６の遅延時間（遅延曲線）の数と値との関係を定めたテーブルが格納されている。遅延時間演算部１８は、ステップ１３２で受け付けた送受信パラメータおよび超音波探触子１１６の種類に対応する、深度範囲ごとの遅延時間（遅延曲線）の数と各深度（整相対象点）ごとの遅延値を求める。

【0063】

例えば、図４のように、送信焦点３０よりも浅い深度範囲３７１においては、１本の遅延曲線３５０が設定され、深い深度範囲３７２においては、５本の遅延曲線３５０～３５４を設定する。

【0064】

（ステップ１３４）
遅延時間演算部１８は、求めた深度範囲ごとの各受信チャンネル１０６の遅延時間（遅延曲線）の数と値を遅延時間格納部１７に格納する。

【0065】

（ステップ１３５）
制御部１１１は、送信焦点３０の位置、送信周波数、送信回数等の送信条件を送信ビームフォーマ１０４に受け渡す。送信ビームフォーマ１０４は、送信信号を生成し、超音波素子アレイ１０１の送信チャンネル１０５の超音波素子に出力する。送信チャンネル１０５の超音波素子はそれぞれ、送信信号を超音波に変換して送信する。超音波素子アレイ１０１の受信チャンネル１０６は、ステップ１３５の送信により生じた被検体からの音波を受信し、受信信号を出力する。

【0066】

（ステップ１３６）
制御部１１１は、アナログ／デジタルコンバータ１１によりデジタル信号に変換された受信信号を受信信号記憶部１３－３に格納する。

【0067】

（ステップ１３７）
遅延加算回路１３－４は、遅延時間格納部１７から各受信チャンネル１０６についての１本目の遅延曲線とその深度範囲を読み込むとともに、受信信号記憶部１３－３から各受信チャンネル１０６の受信信号を読み込む。読み込んだ遅延曲線の深度範囲の各受信信号を、各受信チャンネル１０６の遅延曲線の示す遅延時間で遅延後、一致する深度について遅延後信号を各チャンネル分まとめて加算することにより整相信号を生成する。

【0068】

（ステップ１３８）
生成した整相信号を、整相信号記憶部１３－５の第１メモリ１３－５１に格納する。このとき、その整相信号の深度範囲を示す情報もあわせて格納する。

【0069】

（ステップ１３９）
遅延加算回路１３－４は、上記ステップ１３７、１３８を遅延時間格納部１７に格納されているすべての遅延曲線について遅延加算処理が終了するまで繰り返す。

【0070】

（ステップ１４０）
合成部１５は、整相信号記憶部１３－５の各メモリ１３－５１～１３－５ｎの整相信号を受け取って、合成重み演算部１９が演算により定めた整相信号ごとかつ深度ごとの重みを読み込んで、整相信号を重み付けした後加算する。このとき、整相信号によって存在する深度範囲が異なるため、合成部１５は、整相信号に添付されている深度範囲情報を参照して、一致する深度の整相信号を加算する。

【0071】

画像処理部１０９は、超音波の送信ごとに受け取った加算後の整相信号を、走査線の型（リニア型、コンベックス型、セクタ・フェーズドアレイ型）に従った座標変換・スキャンコンバージョンなどの処理をして２次元・３次元座標空間上に並べたり、対数圧縮などの信号ダイナミックレンジ変換によりピクセルごとの輝度データへの変換処理を施したり、リサンプル・補間・バンドパス理など線形フィルタ処理したり、等をして画像を生成し、接続されている画像表示部１０３に表示させる。

【0072】

なお、遅延加算回路１３－５による遅延加算処理は、上記説明では、遅延曲線ごとに順次行うように説明したが、遅延曲線をそれぞれ予め定めた深度範囲で複数の範囲に分割し、分割した深度範囲ごとに遅延加算処理する構成としてもよい。

【0073】

第３実施形態では、１つの遅延加算回路１３－４のみでありながら、送信焦点３０より浅い深度範囲３７１では、１本の遅延曲線３５０について整相信号を得ることができ、かつ、深い深度範囲３７２においては、複数（例えば５本）の遅延曲線３５０～３５４を用いて、送信ビーム３１０以外の球面波（例えば球面波３０１～３０４）についてそれぞれ整相を得ることができる。

【0074】

よって、第３実施形態では、少ない数の遅延加算回路の規模で、多くの種類の整相信号を１回の送信で得て、高解像度の画像を生成することができる。

【0075】

また、多くの遅延曲線３５０～３５４を演算に用いながら、その深度範囲が限られているため、演算量および処理後の遅延加算後の整相信号のデータ量も、全深度範囲に遅延曲線が設定されている場合よりも、低減することができる。

【0076】

このように、受信ビームフォーマの回路規模が小さく、かつ、演算量も整相信号のデータ量も低減されているため、受信ビームフォーマを超音波探触子１１６に搭載することや、さらには探触子から無線通信で合成部１５に整相信号を送信することも可能である。

【0077】

なお、上述の第１～第３実施形態において、受信ビームフォーマ１３および遅延時間演算部１８は、ハードウエアにより構成することができる。例えば、ＡＳＩＣ（Application Specific Integrated Circuit）のようなカスタムＩＣや、ＦＰＧＡ（Field-Programmable Gate Array）のようなプログラマブルＩＣを用いて、各部の機能を実現するように回路設計を行えばよい。なお、受信ビームフォーマ１３および遅延時間演算部１８は、その一部および全部をソフトウエアでその機能を実現することも可能である。その場合、ＣＰＵ（Central Processing Unit）やＧＰＵ（Graphics Processing Unit）等のプロセッサーと、メモリとを備えたコンピュータ等によって受信ビームフォーマ１３および遅延時間演算部１８を構成し、ＣＰＵが、メモリに格納されたプログラムを読み込んで実行することにより、それらの機能を実現する。

【0078】

図９（ａ）～（ｅ）を用いて、本実施形態の効果を説明する。図９（ａ）の超音波画像１５１は、１本のみの遅延線を用いて超音波ファントムを撮像した超音波画像（比較例）である。図９（ｂ）の超音波画像１５２は、本実施形態の２本以上の遅延線を用いて超音波ファントムを撮像した超音波画像である。比較例の超音波画像１５１と比べて、本実施形態の超音波画像１５２は、超音波ファントム内の血管模擬領域や嚢胞模擬領域（ともに黒い領域）においてそれ以外の組織領域とのコントラストが大きくなっており視認性が向上していることがわかる。また、超音波ファントム内の各反射体の撮像分解能も向上している。

【0079】

グラフ１５３－１５５はそれぞれ、超音波画像１５１中の関心領域（ROI: Region of Interest）１５１１－１５１３および超音波画像１５２中の関心領域１５２１－１５２３について、画像輝度のプロファイルを抽出して並べたグラフである。グラフ１５３は、関心領域１５１１、１５２１の深度方向の画像輝度のプロファイルを、グラフ１５４は、関心領域１５１２、１５２２の方位方向の画像輝度のプロファイルを、グラフ１５５は、関心領域１５１３、１５２３の方位方向画像輝度のプロファイルを示す。グラフ１５３－１５５中の点線のプロファイル１５１１，１５１２，１５１３は、１本のみの遅延線を用いた超音波画像１５１に対応し、実線はプロファイル１５２１，１５２２，１５２３は、本実施形態の２本以上の遅延線を用いた超音波画像１５２に対応している。

【0080】

グラフ１５３およびグラフ１５４をみると、血管領域および嚢胞領域において本実施形態の実線のプロファイル１５２１，１５２２の方が、比較例の点線のプロファイル１５１１，１５１２よりも信号輝度強度が顕著に低くなっており、余分な音響雑音を深さ方向および方位方向の両方において低減する効果を生んでいることが確認できる。またグラフ１５５を見ると、本実施形態の実線のプロファイル１５２３の方が、比較例の点線のプロファイル１５１３よりも、血管模擬領域中の点散乱体の方位方向の幅が顕著に狭くなっており、超音波撮像の分解能が本実施形態によって大きく向上していることが確認できる。

【0081】

以上のように図９（ａ）～（ｅ）より、本実施形態の超音波撮像装置は、得られる超音波画像のコントラストおよび分解能など画質描出能を向上させる効果を有し、これによってより視認性が高い高画質の超音波撮像を実現できることが確認できる。

【符号の説明】

【0082】

１３ 受信ビームフォーマ
１５ 合成部
１７ 遅延時間格納部
１８ 遅延時間演算部
１９ 合成重み演算部
３０ 送信焦点
３６ 受信走査線
１００ 超音波撮像装置
１０１ 超音波素子アレイ
１０２ 超音波撮像装置本体
１０３ 画像表示部
１０４ 送信ビームフォーマ
１０５ 送信チャンネル（超音波素子）
１０６ 受信チャンネル（超音波素子）
１０７ 送受信分離回路（Ｔ／Ｒ）
１０９ 画像処理部
１１０ コンソール
１１１ 制御部
１１６ 超音波探触子
３０１～３０６ 球面波
３１０ 送信ビーム
３５０～３５４ 遅延時間（遅延曲線）
３６１～３６５ 遅延時間（遅延曲線）
３７１、３７２ 深度範囲

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】