(57)【要約】【課題】画像の分解能を向上させること。【解決手段】実施形態に係る超音波診断装置は、超音波プローブと、取得部と、逆フィルタ部を備える。超音波プローブは、超音波を送受信して被検体の走査を行う。取得部は、前記超音波プローブが前記走査を行うことにより得られた第１のデータを取得する。逆フィルタ部は、前記第１のデータに対して、音場を示す情報に基づいて得られた伝達関数を用いてデコンボリューションを行い第２のデータを生成する。【選択図】図１
2019037661.tif 000002
【特許請求の範囲】
【請求項１】
超音波を送受信して被検体の走査を行う超音波プローブと、
前記超音波プローブが前記走査を行うことにより得られた第１のデータを取得する取得部と、
前記第１のデータに対して、音場を示す情報に基づいて得られた伝達関数を用いてデコンボリューションを行い第２のデータを生成する逆フィルタ部と
を備える超音波診断装置。
【請求項２】
前記逆フィルタ部は、前記第１のデータに対してウィーナフィルタを用いた前記デコンボリューションを行い前記第２のデータを生成する、請求項１に記載の超音波診断装置。
【請求項３】
前記逆フィルタ部は、３次元フィルタ、２次元フィルタまたは１次元フィルタを用いて、前記デコンボリューションを行い前記第２のデータを生成する、請求項１に記載の超音波診断装置。
【請求項４】
前記第１のデータは、対数リニア変換が行われた後のＢモードデータまたは対数リニア変換が行われた後のカラーモードデータである、請求項１に記載の超音波診断装置。
【請求項５】
前記伝達関数は、音場計算を行うことにより算出される、請求項１に記載の超音波診断装置。
【請求項６】
前記伝達関数は、点散乱体を用いた測定に基づいて算出される、請求項１に記載の超音波診断装置。
【請求項７】
前記第１のデータは、前記超音波プローブの位置情報を検出しながら前記超音波プローブを動かして収集を行う位置情報付き３Ｄモードによる収集により得られた３次元データ、前記超音波プローブの位置情報を検出せずに前記超音波プローブを動かして収集を行う位置情報なし３Ｄモードによる収集により得られた３次元データ、前記超音波プローブを機械的に搖動させながら収集を行うメカ４Ｄモードによる収集により得られた３次元データ、２次元アレープローブを用いて電子的に３次元走査する２次元アレー３Ｄモードによる収集により得られた３次元データ、または前記超音波プローブの位置を固定して収集を行う２Ｄモードにより得られた２次元データである、請求項１に記載の超音波診断装置。
【請求項８】
前記第１のデータは、３次元ＲＡＷデータ、３次元ボリュームデータ、または２次元ＲＡＷデータである、請求項１に記載の超音波診断装置。
【請求項９】
前記第１のデータは、Ｉ信号データまたはＱ信号データである、請求項１に記載の超音波診断装置。
【請求項１０】
前記第２のデータからボリューム画像、ＭＰＲ（multi-planar  reconstruction）画像または２次元画像のうち少なくとも一つの画像を生成する画像生成部と、
前記画像生成部が生成した前記少なくとも一つの画像をディスプレイに表示させる制御部とを更に備える、請求項１に記載の超音波診断装置。
【請求項１１】
超音波プローブが超音波を送受信して被検体の走査を行うことにより得られた第１のデータを取得する取得部と、
前記第１のデータに対して、前記走査の対象となる場所の音場を示す情報に基づいて得られた伝達関数を用いてデコンボリューションを行い第２のデータを生成する逆フィルタ部とを備えるデータ処理装置。
【請求項１２】
コンピュータに、超音波プローブが超音波を送受信して被検体の走査を行うことにより得られた第１のデータに対して、前記走査の対象となる場所の音場を示す情報に基づいて得られた伝達関数を用いてデコンボリューションを行い第２のデータを生成する処理を実行させるデータ処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、超音波診断装置、データ処理装置及びデータ処理プログラムに関する。
【背景技術】
【０００２】
被検体に超音波を送受信して被検体の組織像や血流像を２次元または３次元画像化して表示する超音波診断装置が実用化されている。信号強度が輝度の明るさで表現されるモードであるＢモード処理においては、被検体に超音波を送受信し、受信信号から組織の輝度情報を検出し、画像化して表示している。また、カラーモード処理（ドプラモード処理）においては、被検体に超音波を送受信し、受信信号から血流のパワー情報、速度情報、分散情報等を検出し、画像化して表示している。
【０００３】
しかしながら、被検体実体が組織像や血流像として表示される迄に処理系が介在しているため、処理系の影響を受けて、組織像や血流像は被検体実体よりも分解能が低下した状態で表示される場合があった。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１５−２００６２号明細書
【発明の概要】
【発明が解決しようとする課題】
【０００５】
本発明が解決しようとする課題は、分解能を向上させることである。
【課題を解決するための手段】
【０００６】
実施形態に係る超音波診断装置は、超音波プローブと、取得部と、逆フィルタ部を備える。超音波プローブは、超音波を送受信して被検体の走査を行う。取得部は、前記超音波プローブが前記走査を行うことにより得られた第１のデータを取得する。逆フィルタ部は、前記第１のデータに対して、音場を示す情報に基づいて得られた伝達関数を用いてデコンボリューションを行い第２のデータを生成する。
【図面の簡単な説明】
【０００７】
【図１】図１は、実施形態に係る超音波診断装置について示した図である。
【図２】図２は、実施形態に係る超音波診断装置の行う処理について説明した図である。
【図３】図３は、実施形態に係る超音波診断装置の行う処理の流れについて示したフローチャートである。
【図４】図４は、実施形態に係る超音波診断装置の行う処理について説明した図である。
【図５】図５は、実施形態に係る超音波診断装置の行う処理について説明した図である。
【図６】図６は、実施形態に係る超音波診断装置の行う処理について説明した図である。
【図７】図７は、比較例に係る超音波診断装置により得られた画像の一例である。
【図８】図８は、実施形態に係る超音波診断装置により得られた画像の一例である。
【図９】図９は、比較例に係る超音波診断装置により得られた画像の一例である。
【発明を実施するための形態】
【０００８】
以下、添付図面を参照して、実施形態に係る超音波診断装置を詳細に説明する。なお、以下の説明において、同様の構成要素には共通の符号を付与するとともに、重複する説明を省略する。
【０００９】
（実施形態）
図１は、第１の実施形態に係る超音波診断装置の構成例を示すブロック図である。図１に例示するように、本実施形態に係る超音波診断装置は、超音波プローブ１と、位置センサ１３と、装置本体１０と、ディスプレイ１３５と、出入力インタフェース１３４とを有する。ここで、装置本体１０は、超音波プローブ１が受信した反射波信号に基づいて超音波画像データを生成する装置であり、超音波送受信部１１と、メモリ１３２と、処理回路１５０とを備える。
【００１０】
超音波プローブ１は、超音波を送受信して被検体Ｐの走査を行う。超音波プローブ１は、例えば、複数の圧電振動子を有し、これら複数の圧電振動子は、装置本体１０が有する超音波送受信部１１から供給される駆動信号に基づき超音波を発生する。また、超音波プローブ１は、被検体Ｐからの反射波を受信して電気信号に変換する。また、超音波プローブ１は、圧電振動子に設けられる整合層と、圧電振動子から後方への超音波の伝播を防止するバッキング材等を有する。なお、超音波プローブ１は、装置本体１０と着脱自在に接続される。
【００１１】
超音波プローブ１から被検体Ｐに超音波が送信されると、送信された超音波は、被検体Ｐの体内組織における音響インピーダンスの不連続面で次々と反射され、反射波信号として超音波プローブ１が有する複数の圧電振動子にて受信される。受信される反射波信号の振幅は、超音波が反射される不連続面における音響インピーダンスの差に依存する。なお、送信された超音波パルスが、移動している血流や心臓壁等の表面で反射された場合の反射波信号は、ドプラ効果により、移動体の超音波送受信方向に対する速度成分に依存して、周波数偏移を受ける。
【００１２】
超音波送受信部１１は、パルス発生器、送信遅延回路、パルサ等を有し、超音波プローブ１に駆動信号を供給する。パルス発生器は、所定のレート周波数で、送信超音波を形成するためのレートパルスを繰り返し発生する。また、送信遅延回路は、超音波プローブ１から発生される超音波をビーム状に集束し、かつ送信指向性を決定するために必要な圧電振動子ごとの遅延時間を、パルス発生器が発生する各レートパルスに対し与える。また、パルサは、レートパルスに基づくタイミングで、超音波プローブ１に駆動信号（駆動パルス）を印加する。すなわち、送信遅延部は、各レートパルスに対し与える遅延時間を変化させることで、圧電振動子面から送信される超音波の送信方向を任意に調整する。
【００１３】
なお、超音波送受信部１１は、処理回路１５０の指示に基づいて、所定のスキャンシーケンスを実行するために、送信周波数、送信駆動電圧等を瞬時に変更可能な機能を有している。特に、送信駆動電圧の変更は、瞬間にその値を切り替え可能なリニアアンプ型の発信回路、又は、複数の電源ユニットを電気的に切り替える機構によって実現される。
【００１４】
また、超音波送受信部１１は、超音波プローブ１により行われた超音波の送受信に基づいて、受信信号Ｐを生成する。かかる受信信号Ｐは、直交位相検波が行われた複素信号であり、その実部はＩ（In  Phase）信号、虚部はQ（Quadrature  Phase）信号と呼ばれる。受信信号Ｐは、ＩＱ信号とも呼ばれる。ある１つのピクセルに着目した場合、I信号の値をI、Q信号の値をQで表わすと、受信信号はＰ＝Ｉ＋jＱと表わされる。
【００１５】
また、位置センサ１３は、超音波プローブ１に装着されるセンサであり、例えば磁気センサである。一例として、図示しない磁場発生器が超音波プローブ１の近傍に設置されると共に超音波プローブ１に位置センサ１３が装着されることで、位置センサ１３は、超音波プローブ１の位置情報、傾き情報をリアルタイムで取得することができる。位置センサ１３は、このようにして取得した超音波プローブ１の位置情報などを、処理回路１５０に送信する。また、処理回路１５０は、取得機能１５０ａにより、位置センサ１３から、超音波プローブ１の位置情報等を取得する。
【００１６】
超音波プローブ１の収集モードについては、様々な収集モードが可能である。超音波プローブ１は、一例として、超音波プローブ１の位置情報を位置センサ１３により検出しながら超音波プローブ１を動かしてデータ収集を行う３Ｄモードによりデータ収集を行い、３次元データを得る。かかる場合、処理回路１５０は、位置センサ１３が検出し、取得機能１５０ａが取得した位置情報を用いて、収集したデータの位置を補正して、３次元画像を生成する。また、別の例として、超音波プローブ１は、超音波プローブ１の位置情報を検出せずに超音波プローブ１を動かして収集を行う位置情報なし３Ｄモードによりデータ収集を行い、３次元データを得る。かかる場合、処理回路１５０は、フレーム（超音波走査面）間隔が一定、またはフレーム間の角度が一定であるという前提で、３次元画像を生成する。また、別の例として、超音波プローブ１は、超音波プローブ１を機械的に搖動させながら収集を行うメカ４Ｄモードによりデータ収集を行い、３次元データを得る。また、別の例として、超音波プローブ１は、２次元アレープローブを用いて電子的に３次元走査する２次元アレー３Ｄモードによりデータ収集を行い、３次元データを得る。また、別の例として、超音波プローブ１は、超音波プローブ１の位置を固定して収集を行う２Ｄモードによりデータ収集を行い、２次元データを得る。
【００１７】
なお、以下、例えば、一次元アレー振動子を用いて方位方向に電子スキャンを行う超音波プローブ１を用い、超音波プローブ１の位置情報を位置センサ１３により取得しながら超音波プローブ１を移動させて３次元画像データを収集し３次元画像を生成する場合について説明する。
【００１８】
処理回路１５０は、取得機能１５０ａ、制御機能１５０ｂ、逆対数変換機能１５０ｃ、対数変換機能１５０ｄ、２乗加算機能１５０ｅ、信号処理機能１５０ｆ、ＭＴＩフィルタ機能１５０ｇ、速度・分散演算機能１５０ｈ、ボリュームデータ生成機能１５０ｉ、画像生成機能１５０ｊ、音場算出機能１５０ｋ、逆フィルタ処理機能１５０ｌ等を有する。実施形態では、取得機能１５０ａ、制御機能１５０ｂ、逆対数変換機能１５０ｃ、対数変換機能１５０ｄ、２乗加算機能１５０ｅ、信号処理機能１５０ｆ、ＭＴＩフィルタ機能１５０ｇ、速度・分散演算機能１５０ｈ、ボリュームデータ生成機能１５０ｉ、画像生成機能１５０ｊ、音場算出機能１５０ｋ、逆フィルタ処理機能１５０ｌ等にて行われる各処理機能は、コンピュータによって実行可能なプログラムの形態でメモリ１３２へ記憶されている。処理回路１５０はプログラムをメモリ１３２から読み出し、実行することで各プログラムに対応する機能を実現するプロセッサである。換言すると、各プログラムを読み出した状態の処理回路１５０は、図１の処理回路１５０内に示された各機能を有することになる。
【００１９】
これらの機能に関しては、複数の独立したプロセッサを組み合わせて処理回路１５０を構成し、各プロセッサがプログラムを実行することにより機能を実現してもよいし、逆に
単一の処理回路１５０にて、これらの処理機能が実現されてもよい。換言すると、各機能ごとに特定の機能が専用の独立したプログラム実行回路に実装される場合であってもよい。例えば、逆対数変換機能１５０ｃ、対数変換機能１５０ｄ、２乗加算機能１５０ｅ、信号処理機能１５０ｆ、ＭＴＩフィルタ機能１５０ｇ、速度・分散演算機能１５０ｈ、逆フィルタ処理機能１５０ｌは、それぞれ逆対数変換器、対数変換器、２乗加算器、信号処理器、ＭＴＩフィルタ、速度・分散演算器、逆フィルタ処理器と呼ばれる専用の回路により実装されてもよい。また、逆に、上述のそれぞれの機能がプログラムとして構成され、１つの処理回路１５０が各プログラムを実行する場合であってもよい。なお、取得機能１５０ａ、制御機能１５０ｂ、画像生成機能１５０ｊ、逆フィルタ処理機能１５０ｌは、それぞれ取得部、制御部、画像生成部、逆フィルタ部の一例である。
【００２０】
上記説明において用いた「プロセッサ」という文言は、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）、ＧＰＵ（Ｇｒａｐｈｉｃａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）或いは、特定用途向け集積回路（Ａｐｐｌｉｃａｔｉｏｎ  Ｓｐｅｃｉｆｉｃ  Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ：ＡＳＩＣ）、プログラマブル論理デバイス（例えば、単純プログラマブル論理デバイス（Ｓｉｍｐｌｅ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｌｏｇｉｃ  Ｄｅｖｉｃｅ：ＳＰＬＤ）、複合プログラマブル論理デバイス（Ｃｏｍｐｌｅｘ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｌｏｇｉｃ  Ｄｅｖｉｃｅ：ＣＰＬＤ）、及びフィールドプログラマブルゲートアレイ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ：ＦＰＧＡ））等の回路を意味する。プロセッサはメモリ１３２に保存されたプログラムを読み出し実行することで機能を実現する。
【００２１】
また、メモリ１３２にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むよう構成しても構わない。この場合、プロセッサは回路内に組み込まれたプログラムを読み出し実行することで機能を実現する。
【００２２】
処理回路１５０は、取得機能１５０ａ、制御機能１５０ｂ、逆対数変換機能１５０ｃ、対数変換機能１５０ｄ、２乗加算機能１５０ｅ、信号処理機能１５０ｆ、ＭＩＴフィルタ機能１５０ｇ、速度・分散演算機能１５０ｈなどにより、超音波送受信部１１により生成された受信信号に基づいて、生データであるＲＡＷデータやボリュームデータ等のデータを生成し、生成されたデータに対する処理を行う。かかる処理について、Ｂモードデータを生成する場合及びカラーモード（ドプラモード）データを生成する場合のそれぞれについて、簡単に説明する。
【００２３】
まず、処理回路１５０が、信号強度が輝度の明るさで表現されるモードであるＢモードにおいて、受信信号からＲＡＷデータ等を生成する場合について説明する。
【００２４】
処理回路１５０は、取得機能１５０ａにより、超音波送受信部１１が生成した受信信号Ｐ＝I＋jQを取得する。続いて、処理回路１５０は、２乗加算機能１５０ｅにより、受信信号Ｐに対して２乗加算を行い、信号値が（I２＋Q２）となる信号を出力する。続いて、処理回路１５０は、対数変換機能１５０ｄにより、２乗加算機能１５０ｅにより例えば２乗加算器から出力された信号に対して対数変換処理を行い、被検体組織の輝度信号の対数圧縮値を出力する。出力された信号値は、例えば10log10(I2+Q2)＝  20log10(I2+Q2)1/2となる。続いて、処理回路１５０は、信号処理機能１５０ｆにより、エッジ強調等の信号処理を行い、ＲＡＷデータを生成する。処理回路１５０は、生成したＲＡＷデータを、メモリ１３２に格納する。
【００２５】
なお、例えば超音波プローブ１の位置情報を検出しながら超音波プローブ１を動かして３次元データを収集する位置情報付き３Ｄモードでは、方位方向に走査することにより１フレームのＲＡＷデータが得られ、走査しながら超音波プローブ１を移動することにより複数フレームから成る３次元ＲＡＷデータが得られる。ここで、手動により超音波プローブ１の移動が行われているので、収集された複数フレームは必ずしも等間隔なデータとはならない。従って、かかる場合、処理回路１５０は、ボリュームデータ生成機能１５０ｉにより、信号処理機能１５０ｆにより生成されたＲＡＷデータから、等間隔のボリュームデータを生成する。また、処理回路１５０は、得られたボリュームデータを、メモリ１３２に格納する。なお、超音波プローブ１を機械的に搖動させながら３次元データを収集する場合や２次元アレープローブを用いて電子的に３次元走査し３次元データを収集する場合は、通常、フレーム間隔は一定に制御される。従って、かかる場合、ボリュームデータを生成せずにＲＡＷデータから直接ＭＰＲ画像データ等を生成することができる。
【００２６】
図２は、実施形態に係る超音波診断装置の行う処理について説明した図である。図２において、処理回路１５０が、ボリュームデータ生成機能１５０ｉにより、ＲＡＷデータ２からボリュームデータ３を生成する様子が模式的に示されている。図２において、Ｘ’軸が距離方向に対応し、Ｙ’方向が方位方向に対応し、Ｚ’方向がスライス方向に対応する。図２において、Ｘ’―Ｙ’面がフレームを表し、Ｚ’軸はフレームの移動方向を表す。処理回路１５０は、ボリュームデータ生成機能１５０ｉにより、ＲＡＷデータ２に基づいて、等間隔に並ぶ複数のフレームから構成され、各フレームはＲＡＷデータ２と同様のフォーマットになっているデータとして、ボリュームデータ３を生成する。かかるボリュームデータ３は、ラスタ（超音波走査線）に対応したデータであるラスタデータをフレーム毎に持ち、位置の異なる複数のフレームで３次元データが構成されている。この意味でこの様なボリュームデータを３次元ＲＡＷデータと呼ぶこともある。このように、ボリュームデータ３のデータのフォーマットを、ＲＡＷデータ２のフレームの集まりとなるようなフォーマットとすることで、他の画像処理などにそのまま再利用することが可能となり、データの汎用性を高めることができる。
【００２７】
また、この時、処理回路１５０は、ボリュームデータ生成機能１５０ｉにより、ＲＡＷデータ２の各フレームを、位置センサ１３から得られた位置情報を用いて位置合わせして配置する。また、処理回路１５０は、ボリュームデータ生成機能１５０ｉにより、ボリュームデータ３のデータを、例えばボリュームデータ３が生成されるべきピクセルを挟む隣り合う２フレームのＲＡＷデータ２から補間により生成する。
【００２８】
また、画像生成機能１５０ｊを有する処理回路１５０は、ボリュームデータ生成機能１５０ｉにより生成されたボリュームデータ３に対して補間処理を行うことにより例えば、ＭＰＲ（multi-planar  reconstruction）画像等その他のフォーマットのデータまたは画像を生成する。処理回路１５０は、画像生成機能１５０ｊにより、当該ＭＰＲ画像として、例えば直交３断面画像を生成する。加えて、処理回路１５０は、画像生成機能１５０ｊにより、ボリュームデータ３に対してレンダリング処理等を行うことにより、ＶＲ（Volume  Rendering）画像やＭＩＰ（Maximum  Intensity  Projection）画像等の２次元画像等を更に生成してもよい。また、処理回路１５０は、制御機能１５０ｂにより、画像生成機能１５０ｊにより生成されたボリューム画像、ＭＰＲ画像、または２次元画像等を、ディスプレイ１３５に表示させてもよい。
【００２９】
次に、処理回路１５０が、カラーモード（ドプラモード）において、受信信号からＲＡＷデータ等を生成する場合について説明する。
【００３０】
カラーモード（ドプラモード）においては、超音波送受信部１１は、超音波プローブ１を通じて、被検体に対して同じ方向に複数回（Ｎ回）の超音波の送受信を行い、受信信号Ｐを生成する。処理回路１５０は、取得機能１５０ａにより、超音波送受信部１１が生成した複数の受信信号Ｐｉ（ｉ＝１，・・・，Ｎ）を取得する。このとき、超音波送受信部１１が生成する複数の受信信号Ｐｉのそれぞれは、直交位相検波が行われた複素信号である。受信信号Ｐの実部はＩ信号、受信信号Ｐの虚部は、Ｑ信号とも呼ばれる。また、受信信号Ｐは、ＩＱ信号とも呼ばれる。ある１つのピクセルに着目した場合、受信信号Ｐの実部の値をＩ、受信信号Ｐの虚部の値をＱで表わし、Ｎを送信の回数、ｉをｉ番目の送信とし、その時の受信信号をＰｉで表すと、ｉ番目の受信信号Ｐｉは、Ｐｉ＝Ｉｉ＋ｊＱｉ（ｉ＝１，・・・，Ｎ）と表わされる。
【００３１】
この受信信号Ｐｉには、目的とする血流からのドプラ信号Ib、Qbの他に、組織エコー等に起因するクラッタと呼ばれる不要信号Ic、Qcが含まれている。即ち、ｉ番目の受信信号Ｐｉは、以下の式（１）で表される。
【数１】
2019037661.tif 000003
【００３２】
そこで、処理回路１５０は、ＭＴＩフィルタ機能１５０ｇにより、取得機能１５０ａにより取得した一連の受信信号Ｐｉ（ｉ＝１，・・・，Ｎ）に対して、ＭＴＩ(Moving  Target  Indication)フィルタを適用し、クラッタの除去された血流信号を出力する。ＭＴＩフィルタ機能１５０ｇにより出力される信号は、例えば以下の式（２）で表される。
【数２】
2019037661.tif 000004
【００３３】
続いて、血流の絶対値を算出するモードであるパワーモードの場合、処理回路１５０は、２乗加算機能１５０ｅにより、以下の式（３）のように、複数の受信信号に対して２乗加算処理を行い、結果を出力する。
【数３】
2019037661.tif 000005
【００３４】
続いて、処理回路１５０は、対数変換機能１５０ｄにより、２乗加算機能１５０ｅにより例えば２乗加算器から出力された信号に対して対数変換処理を行い、被検体組織の輝度信号の対数圧縮値を出力する。処理回路１５０が出力する信号値は、例えば10log10(Power)＝10log10(Σi=1N(Ibi2+Qbi2))＝  20log10(Σi=1N(Ibi2+Qbi2))1/2となる。続いて、処理回路１５０は、信号処理機能１５０ｆにより、エッジ強調等の信号処理を行い、ＲＡＷデータを生成する。処理回路１５０は、生成したＲＡＷデータを、メモリ１３２に格納する。
【００３５】
また、血流速度や血流速度の分散を算出するモードである速度モードや速度-分散モードの場合、処理回路１５０は、ＭＴＩフィルタ機能１５０ｇにより出力されたクラッタの除去された血流信号Pbi  (i=1,・・・,N）に対して、速度・分散演算機能１５０ｈにより、血流速度や、血流速度の分散を算出する。具体的には、処理回路１５０は、ＭＴＩフィルタ機能１５０ｇにより出力されたクラッタの除去された血流信号Pbi  (i=1,・・・,N）に対して例えば自己相関法を適用することにより、血流速度や、血流速度の分散を算出し、血流速度を表すデータや、血流速度の分散を表すデータを生成する。続いて、処理回路１５０は、信号処理機能１５０ｆにより、エッジ強調等の補助的な信号処理を行い、血流速度を表すＲＡＷデータや、血流速度の分散を表すＲＡＷデータを生成する。処理回路１５０は、生成したＲＡＷデータを、メモリ１３２に格納する。
【００３６】
また、Ｂモードの場合と同様、カラーモードの場合においても、処理回路１５０は、ボリュームデータ生成機能１５０ｉにより、信号処理機能１５０ｆにより生成されたＲＡＷデータから、ボリュームデータを生成し、得られたボリュームデータを、メモリ１３２に格納する。また、処理回路１５０は、画像生成機能１５０ｊにより、生成されたボリュームデータを基に、例えば、ＭＰＲ（multi-planar  reconstruction）画像、ＶＲ（Volume  Rendering）画像やＭＩＰ（Maximum  Intensity  Projection）画像等の２次元画像等を更に生成し、制御機能１５０ｂにより、画像生成機能１５０ｊにより生成されたボリューム画像、ＭＰＲ画像、または２次元画像等を、ディスプレイ１３５に表示させてもよい。
【００３７】
また、処理回路１５０は、音場算出機能１５０ｋにより、音場を算出する。また、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、ＩＱデータ等の信号データ、ＲＡＷデータ、ボリュームデータ、ＭＰＲデータ等に対してデコンボリューションを行って、画質が向上されたデータを生成する。これらの処理の詳細については後述する。
【００３８】
メモリ１３２は、信号データを記憶する信号データメモリ、ＲＡＷデータを記憶するＲＡＷデータメモリ、ボリュームデータを記憶するボリュームデータメモリ等である。メモリ１３２は、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）、フラッシュメモリ等の半導体メモリ素子、ハードディスク、光ディスク等である。また、メモリ１３２は、処理回路１５０が有する各機能に対応するプログラムを記憶する。
【００３９】
ディスプレイ１３５は、超音波診断装置の操作者が出入力インタフェース１３４を用いて各種設定要求を入力するためのＧＵＩ（Graphical  User  Interface）を表示したり、装置本体１０において生成された各種画像データ等を表示したりするためのディスプレイである。ディスプレイ１３５は、例えば、液晶表示器等の表示デバイスである。
【００４０】
出入力インタフェース１３４は、データ出入力インタフェース及びネットワークインタフェースにより構成され、操作者からの各種指示や情報入力を受け付けるための出入力インタフェースである。出入力インタフェース１３４は、例えば、マウス、キーボード、ボタン、パネルスイッチ、タッチコマンドスクリーン、フットスイッチ、トラックボール、ジョイスティック等で構成される。
【００４１】
続いて、実施形態に係る超音波診断装置が行う処理の詳細について説明する。
【００４２】
実施形態に係る超音波診断装置において、処理回路１５０は、取得機能１５０ａにより、超音波プローブ１が被検体に対して超音波を送受信して走査を行うことにより得られた第１のデータを取得し、逆フィルタ処理機能１５０ｌにより、取得された第１のデータに対して、走査の対象となる場所の音場を示す情報に基づいて得られた伝達関数を用いてデコンボリューションを行い第２のデータを生成する。特に、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、取得された第１のデータに対してウィーナフィルタを用いたデコンボリューションを行い第２のデータを生成する。
【００４３】
かかる背景について、簡単に説明する。
【００４４】
まず、実施形態に係る超音波診断装置がデコンボリューションを行う背景及び実施形態に係る超音波診断装置がウィーナフィルタを用いたデコンボリューションを行う背景について簡単に説明する。
【００４５】
入力信号ｆに対して処理ｈが行われた時の出力信号gは、畳込み（コンボリューション）を用いて、g=f*hと書くことができる。例えば、超音波画像を取り扱う場合を想定し、入力信号ｆ及び出力信号ｇが２次元データである場合、x,yを空間座標とし、入力信号をf(x,y)とし、処理系を表す関数をh(x,y)とし、出力信号をg(x,y)とすると、出力信号は、以下の式（４）のように、入力信号f(x,y)と処理系を表す関数h(x,y)との畳込みで書くことができる。
【数４】
2019037661.tif 000006
【００４６】
畳み込み積分のフーリエ変換は、各々のフーリエ変換の積となる。従って、u,vを空間周波数とし、f(x,y)、ｇ(x,y)、h(x,y)のフーリエ変換をそれぞれF(u,v)、H(u,v)、G(u,v)とすると、式（４）をフーリエ変換すると、以下の式（５）が成り立つ。
【数５】
2019037661.tif 000007
【００４７】
H(u,v)は、伝達関数とも呼ばれる。
【００４８】
ここで、F(u,v)は、被検体実体の原画像のフーリエ変換であり、H(u,v)は、超音波処理系の伝達関数であり、G(u,v)は、表示される超音波画像のフーリエ変換に対応する。超音波画像のフーリエ変換であるG(u,v)は、処理系のH(u,v)の影響により一般にF(u,v)よりも分解能が低下する。
【００４９】
ここで、G(u,v)からF(u,v)を得るような処理Hｉｎｖ(u,v)を用いると、処理回路１５０は、被検体実体の原画像のフーリエ変換F(u,v)を、以下の式（６）のように、超音波画像のフーリエ変換G(u,v)に基づいて算出することができる。
【数６】
2019037661.tif 000008
【００５０】
上述の式（６）の処理により、処理回路１５０は、超音波画像のフーリエ変換G(u,v)から、より被検体実体を反映していると考えられるデータであるF(u,v)を再構成することができる。かかる操作をデコンボリューションと呼び、Hinv(u,v)を逆フィルタと呼ぶ。
【００５１】
換言すると、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、取得された第１のデータと、伝達関数H(u,v)に基づく値Hinv(u,v)とを用いてデコンボリューションを行い、画質が向上した第２のデータを生成する。
【００５２】
ここで、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、取得された第１のデータに対してウィーナフィルタを用いたデコンボリューションを行い第２のデータを生成する。具体的には、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、ステップＳ１００において得られた第１のデータと、伝達関数H(u,v)に基づくウィーナフィルタを用いたデコンボリューションを行い第２のデータを生成する。かかる理由について説明する。
【００５３】
式（５）と式（６）とを比較すると、逆フィルタＨｉｎｖ（u,v）は、以下の式（７）のように伝達関数H(u,v)の逆数で与えられると考えられる。
【数７】
2019037661.tif 000009
【００５４】
しかし、伝達関数H(u,v)が0または0に近い値をとるときはHｉｎｖ(u,v)が発散し、表示される超音波画像のフーリエ変換G(u,v)に含まれるノイズ成分が増幅されて、デコンボリューション後のデータF(u,v)に大きなノイズが現れる。
【００５５】
このため、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、取得された第１のデータに対してウィーナフィルタHw(u,v)を用いたデコンボリューションを行い第２のデータを生成する。
【００５６】
ここで、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、原画像F(u,v)に近い復元画像F'(u,v)を、以下の式（８）により計算する。
【数８】
2019037661.tif 000010
【００５７】
ここで、ウィーナフィルタHw(u,v)は、式（５）にノイズ成分N(u,v)を加味して、以下の式（９）とおき、原画像F(u,v)と復元画像F'  (u,v)との平均２乗誤差が最小となるという条件のもとで得られる関数であり、式（１０）で与えられる。
【数９】
2019037661.tif 000011
【数１０】
2019037661.tif 000012
【００５８】
これにより、処理回路１５０は、デコンボリューションを行うに当たって０または０に近い値で割るという問題を回避することができ、処理回路１５０が逆フィルタ処理機能１５０ｌにより生成する第２のデータの数値的安定性を保つことができる。なお、式（１０）において、以下の式（１１）のように、数値パラメータΓを導入した。
【数１１】
2019037661.tif 000013
【００５９】
このように、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、ステップＳ１００において得られた第１のデータに対して、伝達関数Ｈ(u,v)と、パラメータΓとにより定まるウィーナフィルタＨｗ(u,v)とに基づいてデコンボリューションを行い、第２のデータを生成する。なお、パラメータΓの値は、例えば０．００５等、０に近い値が用いられる。一例として、超音波画像に対してパラメータΓの値を変えて第２のデータを複数生成し、生成される第２のデータが妥当なデータとなるようなパラメータΓの値を、ウィーナフィルタＨｗ(u,v)に用いるΓの値に設定し、設定した値に基づき、デコンボリューションを行っても良い。
【００６０】
なお、パラメータΓ＝０の極限では、ウィーナフィルタHw(u,v)は、式（７）による逆フィルタHinv(u,v)に等しくなる。
【００６１】
また、フィルタの性能や適用範囲、パフォーマンス等を考慮し、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、ウィーナフィルタHw(u,v)以外のフィルタ、例えばルーシー・リチャードソンアルゴリズムなどの処理を行い、第２のデータを生成してもよい。
【００６２】
続いて、実施形態に係る超音波診断装置が、音場を示す情報に基づいて得られた伝達関数を用いてデコンボリューションを行う背景について説明する。
【００６３】
一般に、超音波診断装置の処理系は、複数の処理系が縦列につながっている。それらの処理系が例えば、音場に対応する応答関数（フィルタ係数）、信号処理に対応する応答関数（フィルタ係数）、画像処理に対応する応答関数（フィルタ係数）の３つである場合を考え、音場に対応する応答関数（フィルタ係数）をh１(x,y)、信号処理に対応する応答関数（フィルタ係数）をh２(x,y)、画像処理に対応する応答関数（フィルタ係数）をh３(x,y)、とすると、式（４）は、以下の式（１２）で与えられる。
【数１２】
2019037661.tif 000014
【００６４】
また、式（１２）をフーリエ変換すると、以下の式（１３）が成り立つ。
【数１３】
2019037661.tif 000015
【００６５】
ここで、H１(u,v)は音場に対応する伝達関数であり、H２(u,v)は例えばそれに続く信号処理に対応する伝達関数であり、H３(u,v)は画像処理に対応する伝達関数である。
【００６６】
ここで、音場に対応する伝達関数H1(u,v)、信号処理に対応する伝達関数H2(u,v)、画像処理に対応する伝達関数H3(u,v)の中で、音場に対応する伝達関数H1(u,v)が、最終的な画像の分解能の向上に最も寄与が大きいと考えられる。また、音場に対応する伝達関数H1(u,v)やその逆フーリエ変換であるフィルタ係数h1(x,y)は、計算や実測によって近似的な値を比較的容易に得ることができる。従って、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、音場に対応する伝達関数H1(u,v)を用いてデコンボリューションを行い、第２のデータを生成する。例えば、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、音場に対応する伝達関数H1(u,v)を用いたウィーナフィルタH１ｗ(u,v)を用いてデコンボリューションを行い、以下の式（１４）の左辺で与えられる第２のデータＦ''(u,v)を生成する。
【数１４】
2019037661.tif 000016
【００６７】
次に、図３〜６を用いて、実施形態に係る超音波診断装置が行う処理の詳細について説明する。図３は、実施形態に係る超音波診断装置の行う処理の流れについて示したフローチャートである。また、図４は、実施形態に係る超音波診断装置の行う処理について、模式的に説明した図である。
【００６８】
図３のフローチャートにおいて、はじめに、処理回路１５０は、取得機能１５０ａ等により、超音波プローブが被検体の走査を行うことにより得られた、信号データまたは画像データである第１のデータを取得する（ステップＳ１００）。ここで、第１のデータとは、例えばＢモードデータまたはカラーモードデータであり。より具体的には、第１のデータとは、例えば対数値をリニアな値に変換する変換である対数リニア変換が行われた後のＢモードデータや、対数値をリニアな値に変換する変換である対数リニア変換が行われた後のカラーモードデータである。より具体的には、第１のデータとは、例えば対数リニア変換が行われた後のＢモードのＲＡＷデータや、対数リニア変換が行われた後のカラーモードのＲＡＷデータである。
【００６９】
前述したように、Ｂモードの場合、処理回路１５０は、取得機能１５０ａにより、超音波送受信部１１が生成した受信信号Ｐを取得し、続いて、２乗加算機能１５０ｅにより、２乗加算を行い、対数変換機能１５０ｄにより、対数変換処理を行い、続いて、信号処理機能１５０ｆにより、エッジ強調等の信号処理を行うことによりＢモードのＲＡＷデータを生成し、生成したＢモードのＲＡＷデータを、メモリ１３２に格納する。
【００７０】
ここで、生成されるＲＡＷデータは、対数圧縮されているデータである。処理回路１５０は、逆フィルタ処理機能１５０ｌにより、対数圧縮されているデータであるＢモードのＲＡＷデータに対してそのままデコンボリューションを行っても分解能を向上させることは可能である。しかしながら、リニアデータに変換させたデータに対してデコンボリューションを行いその後対数データに戻した方が、処理回路１５０は、より分解能改善効果を高めることができる。従って、処理回路１５０は、逆対数変換機能１５０ｃにより、メモリ１３２に格納されている対数圧縮されているＲＡＷデータに対して、対数値（20log10(I2+Q2)1/2）をリニアな信号値（(I2+Q2)1/2）に変換する変換である対数リニア変換を行い、対数圧縮されたＲＡＷデータをリニアデータに変換する。
【００７１】
図４の（ａ）及び（ｂ）に、かかる状況が示されている。図４において、図４の上下方向が距離方向に対応し、左右方向が方位方向に対応し、前後方向がスライス方向を示す。また、図４は、処理回路１５０が、Ｂモードの３次元ＲＡＷデータに対して、後述のステップＳ１２０で、方位方向、スライス方向の２次元フィルタを逆フィルタとして用いて、デコンボリューションを行う場合の例について示している。なお、実施形態はこれに限られず、例えばカラーモードデータに対しても、実施形態は同様に適用可能である。
【００７２】
図４（ａ）において、ＲＡＷデータ２ａ、２ｂ、２ｃ、２ｄ、２ｅのそれぞれは、処理回路１５０が信号処理機能１５０ｆ等の処理により生成したＲＡＷデータである。別の例として、ＲＡＷデータ２ａ、２ｂ、２ｃ、２ｄ、２ｅのそれぞれは、処理回路１５０が信号処理機能１５０ｆ等の処理により生成し、メモリ１３２に蓄えられたデータであってもよい。ＲＡＷデータ２ａ、２ｂ、２ｃ、２ｄ、２ｅのそれぞれは、対数圧縮されているデータであることから、処理回路１５０は、ＲＡＷデータ２ａ、２ｂ、２ｃ、２ｄ、２ｅのそれぞれに対して対数リニア変換を実行することにより、リニアデータに変換されたデータ４ａ、４ｂ、４ｃ、４ｄ、４ｅを第１のデータとして取得し、後続のステップでデコンボリューションを行う。
【００７３】
図３に戻り、処理回路１５０は、音場算出機能１５０ｋにより、被検体の走査の対象となる場所の音場を示す情報に基づいて、音場に対する伝達関数Ｈ１(u,v)を算出する（ステップＳ１１０）。音場に対する伝達関数の求め方には、例えば以下の様な方法がある。なお、３次元の場合はh(x,y,z)、Ｈ１(u,v,w)等となるが、煩雑さを避ける為２次元の表記を用いて説明を行う。
【００７４】
式（４）において、入力信号f(x,y)をデルタ関数δ(x,y)とした場合、デルタ関数のフーリエ変換は1になるので、F(u,v)=1となり、式（４）をフーリエ変換すると、式（５）は、以下の式（１５）で表される。
【数１５】
2019037661.tif 000017
【００７５】
すなわち、出力信号G(u,v)は処理系の伝達関数H(u,v)そのものになる。すなわち、被検体として点散乱体δ(x,y)を備えたファントムを用いると、得られる画像は、伝達関数Ｈ(u,v)のフーリエ変換h(x,y)を表わす。このh(x,y)はPSF（Point  Spread  Function、点拡がり関数）と呼ばれる。
【００７６】
従って、一例として、処理回路１５０は、音場算出機能１５０ｋにより、例えば点散乱体を用いた測定に基づいて、ＰＳＦであるフィルタ係数h(x,  y)を算出する。例えば、処理回路１５０は、音場算出機能１５０ｋにより、点散乱体を３次元的に配置したファントムを用いて行われた超音波走査により得られた画像データに基づいて、フィルタ係数h(x,  y)を算出する。一例として、超音波診断装置は、代表的な超音波送受信条件等で、点散乱体を用いた測定を行い、その測定値をメモリ１３２に格納する。処理回路１５０は、メモリ１３２に格納した測定値を読み込み、それらを内挿、外挿（伸縮）、またはスケールさせることにより、フィルタ係数h(x,y)を算出する。
【００７７】
また、別の例として、処理回路１５０は、音場算出機能１５０ｋにより、各深さ（距離）ｘｉでの音場の形状とビーム幅を求め、算出した音場の形状とビーム幅等から、フィルタ係数h(x,y)を算出する。
【００７８】
続いて、処理回路１５０は、音場算出機能１５０ｋにより、算出されたフィルタ係数h(x,y)に基づいて、音場に対応するフィルタ係数h1(x,y)を算出する。一例として、処理回路１５０は、音場算出機能１５０ｋにより、処理系において最も影響が大きいのは音場であると考えられるので、音場に対応するフィルタ係数h1(x,y)とフィルタ係数h(x,y)とが等しいとして、算出されたフィルタ係数h(x,y)に基づいて、音場に対応するフィルタ係数h1(x,y)を算出する。続いて、処理回路１５０は、音場算出機能１５０ｋにより、音場に対応するフィルタ係数h1(x,y)をフーリエ変換することにより、音場に対する伝達関数Ｈ１（u,ｖ）を算出する。このようにして、処理回路１５０は、音場算出機能１５０ｋにより、点散乱体を用いた測定に基づいて、走査の対象となる場所の音場を示す情報に基づいて、伝達関数Ｈ１(u,v)を算出する。
【００７９】
また、別の例として、処理回路１５０は、音場算出機能１５０ｋにより、音場計算を行うことにより、音場に対応するフィルタ係数h1(x,y)を算出する。
【００８０】
音場計算を行うことにより、音場に対応するフィルタ係数h1(x,y)を算出する第１の方法として、処理回路１５０が、音場算出機能１５０ｋにより、例えば音場計算ソフトなどを用いて、超音波送信の送受信条件を用いて深さ毎の音場を計算し、計算された深さ毎の音場を用いてフィルタ係数h1(x,y)を算出する方法がある。第１の方法は、算出されるフィルタ係数の精度は高くなるが、負荷は大きくなる方法となる。
【００８１】
また、音場計算を行うことにより、音場に対応するフィルタ係数h1(x,y)を算出する第２の方法として、処理回路１５０が、音場算出機能１５０ｋにより、例えば音場計算ソフトなどを用いて、超音波送信の送受信条件を用いてある間隔の深さで音場を計算し、その間の深さの値は、それを挟む深さの音場値から補間により求め、計算された深さ毎の音場を用いてフィルタ係数h1(x,y)を算出する方法がある。第２の方法は、第１の方法に比べて計算の精度は低下するが、負荷は軽くなる方法となる。
【００８２】
また、音場計算を行うことにより、音場に対応するフィルタ係数h1(x,y)を算出する第３の方法として、処理回路１５０が、音場算出機能１５０ｋにより、代表的な送受信条件で、例えばFundamental音場とharmonic音場とについて、音場に対応するフィルタ係数h1(x,y)を算出しておき、算出した値を、メモリ１３２に格納する方法がある。超音波送受信条件などのパラメータ値が変わった場合、処理回路１５０は、音場算出機能１５０ｋにより、メモリ１３２に格納されている、代表的な送受信条件での音場に対応するフィルタ係数h1(x,y)を基に、それらを内挿、外挿（伸縮）、またはスケールさせることにより、超音波送受信条件などのパラメータ値が変わった場合における音場に対応するフィルタ係数h1(x,y)を算出する。一例として、音場の拡散はλ/Ｄ（λ：波長、Ｄ：開口幅）に比例する。即ち、例えばフォーカス点付近のビーム幅は、開口幅が２倍になれば１／２になり、周波数が２倍になれば（波長が半分になれば）１／２になる。よって、開口幅や周波数などのパラメータ値が変わっても、処理回路１５０は、音場算出機能１５０ｋにより、メモリ１３２に格納されている、代表的な送受信条件での音場に対応するフィルタ係数h1(x,y)を援用して、音場に対応するフィルタ係数h1(x,y)を算出することができる。
【００８３】
なお、上述の第１の方法〜第３の方法は、方位方向とスライス方向についての、音場に対応するフィルタ係数についての計算方法を示している。距離方向のフィルタ係数の方法については、例えば、処理回路１５０は、音場算出機能１５０ｋにより、超音波パルス形状を計算で取得し、取得した超音波パルス形状を基に、フィルタ係数を算出する。また、処理回路１５０は、音場算出機能１５０ｋにより、被検体の減衰効果等、例えば周波数依存性減衰の効果を繰込み、精度を向上してもよい。
【００８４】
また、上述の第１の方法〜第３の方法は、装置の能力に応じて方法を選択可能としてもよい。例えば、計算能力の面で高負荷な計算にも適しているデータ処理装置では精度の高い方法を採用し、逆に計算能力の面で高負荷な計算に適していない超音波診断装置では、リアルタイム性を損なわない方法を選択してもよい。
【００８５】
続いて、処理回路１５０は、音場算出機能１５０ｋにより、音場に対応するフィルタ係数h1(x,y)をフーリエ変換することにより、音場に対する伝達関数Ｈ１（u,ｖ）を算出する。このようにして、処理回路１５０は、音場算出機能１５０ｋにより、音場計算を行うことにより走査の対象となる場所の音場を示す情報を算出し、算出した情報に基づいて、伝達関数Ｈ１(u,v)を算出する。
【００８６】
図５及び図６を用いて、音場に対応するフィルタ係数h1(x,y)の算出について更に説明する。
【００８７】
まず、図５を用いて、音場と分解能との関係性を簡単に述べ、続いて、音場と分解能との関係性を踏まえて、音場に対応するフィルタ係数h1の次元の選択について述べる。
【００８８】
まず、図５を用いて、音場と分解能との関係性を説明する。音場は、分解能に大きく影響する。音場は、深さ方向である距離方向と、超音波を走査する方向である方位方向と、超音波走査面と垂直な方向であるスライス方向との３方向の分布を持つ。音場の３方向の分布形状は、一般に大きく異なる。
【００８９】
音場の距離方向の分布形状は送信される超音波パルスの形状で決まり、距離分解能はパルスの長さで決定される。通常は分解能を良くするため、パルスの長さは、１〜３波長程度となるように、短く設定されている。
【００９０】
図５（ａ）は、上下方向を距離方向とし、左右方向を方位方向とした場合の、送信される超音波ビームの模式図を表している。音場の方位方向の分布形状は送受信音場の方位方向の形状で決まり、方位分解能は方位方向の音場の幅２０で決まる。一般には、スライス方向よりも開口を大きく設定できるので、スライス方向より超音波ビームを狭くでき、スライス方向より分解能が良いが、距離方向よりは分解能は悪くなる。通常は、送信は１つのフォーカス点を設定して行われるが、受信は深さ毎にフォーカス点を移動するダイナミックフォーカスを行うので、方位分解能は受信音場の分解能が支配的である。
【００９１】
図５（ｂ）は、上下方向を距離方向とし、左右方向をスライス方向とした場合の、送信される超音波ビームの模式図を示している。音場のスライス方向の分布形状は送受信音場のスライス方向の形状で決まり、スライス分解能はスライス方向の音場の幅２１で決まる。通常は、超音波振動子のスライス方向の素子数は１素子で、音響レンズで焦点を形成しているので、送受信ともフォーカス点は１点で同じ位置にあり、スライス分解能は送受信音場の影響を同程度に受ける。また、１．５Ｄアレーと呼ばれるプローブもある。この場合は、スライス方向に振動子を３素子、５素子等に分割している。また、２次元アレーの場合には、スライス方向に振動子が多素子に分割されている。この様な場合、受信はダイナミックフォーカスを行うことができるので、１素子の場合よりも広い深さの範囲でスライス分解能が向上する。
【００９２】
よって、一般に分解能は、距離方向が最も良くなり、方位方向、スライス方向となるにつれて悪くなる。
【００９３】
次に、分解能を踏まえて、処理回路１５０がデコンボリューションを行う逆フィルタの次元について述べる。まず、第１のデータが３次元画像データの場合、計算負荷を問題にしない場合、フィルタの次元が３次元であるのが、最も分解能の改善効果が大きい。かかる場合、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、第１のデータに対して距離方向、方位方向及びスライス方向すべてに対してデコンボリューションを行う３次元フィルタを用いて第２のデータを生成する。計算負荷を低減させる場合、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、第１のデータに対して２次元フィルタを用いて第２のデータを生成するが、その場合、例えば分解能がもともと比較的良い距離方向を除いた、方位方向及びスライス方向の２次元フィルタを用いて第２のデータを生成してもよい。また、更に計算負荷を低減させる場合、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、第１のデータに対して１次元フィルタを用いて第２のデータを生成するが、その場合、例えば分解能の影響がもともと最も大きいスライス方向の１次元フィルタを用いて第２のデータを生成してもよい。
【００９４】
同様に、第１のデータが走査面が距離方向と方位方向から成る２次元画像データの場合、計算負荷を問題にしない場合、フィルタの次元が２次元であるのが、最も分解能の改善効果が大きい。かかる場合、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、第１のデータに対して距離方向、方位方向に対してデコンボリューションを行う２次元フィルタを用いて第２のデータを生成する。計算負荷を低減させる場合、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、第１のデータに対して１次元フィルタを用いて第２のデータを生成するが、その場合、例えば分解能の影響がもともと大きい方位方向の１次元フィルタを用いて第２のデータを生成してもよい。
【００９５】
次に、図６を用いて、音場分布から音場に対応するフィルタ係数h1の係数を算出する方法について更に詳細に説明する。図６は、実施形態に係る超音波診断装置の行う処理について説明した図である。具体的には、図６には、ｘ方向を距離方向、ｙ方向を方位方向、ｚ方向をスライス方向として、ｙｚ面での２次元フィルタに関するフィルタ係数ｈ１を算出する場合について、フィルタｈ１(y,z)のカーネルサイズが5*5の例を用いて説明する。
【００９６】
図６（ａ）は、ＲＡＷデータを上面から見た図であり、図６（ａ）の丸一つ一つはラスタを表している。また、図６（ｂ）は、ＲＡＷデータを正面から見た図である。方位方向のラスタピッチについて、リニア走査では、ラスタピッチ３１は一定であり、以下Δyとおく。また、コンベックス走査、セクタ走査の場合は、ラスタピッチは深さx=xi毎に異なるが、各深さでは一定であり、これをΔyとおく。ここで、議論の簡略化のため、添え字ｉは省略している。
【００９７】
また、スライス方向のラスタピッチについて、操作者が、スライス方向に均一でより高い分解能を得るために、プローブを遅い・均一な速度で一定方向に動かす場合（平行移動の場合）について考える。かかる場合、プローブを略スライス方向に手動で移動するので、スライス方向のラスタピッチ３０ａ、３０ｂ等は一般には異なる。これをΔz1、Δz2、Δz3、Δz4、等とおく。しかし、プローブが動く速度はほぼ均一な速度と考えて良いので、各フレームはほぼ平行で間隔はほぼ均一と考えて、Δz=mean(Δzk)    (k=1,2,・・・,(M-1)、Mはフレーム枚数)を用いてもよい。
【００９８】
また、スライス方向のラスタピッチについて、操作者が、スライス方向に均一でより高い分解能を得るために、プローブを遅い・均一な角度で煽る場合（煽る場合）について考える。この場合、プローブを略スライス方向に手動で移動するので、スライス方向のラスタピッチは一般には異なる。また、深さによっても異なるが、或る深さx=xi（実空間上では円弧を成す）に着目してこれをΔz1、Δz2、Δz3、Δz4、とおく。しかし、プローブは均一な角度で煽られているので、各フレーム間の角度はほぼ均一と考えて、Δz=mean(Δzk)    (k=1,2,・・・,(M-1)、Mはフレーム枚数)を用いてもよい。Δzは深さx=xiによって異なる。ここで、議論の簡略化のため、添え字ｉは省略している。
【００９９】
次に、図６（ｃ）〜（ｅ）を用いて、フィルタ係数h1の係数の算出について説明する。図６（ｃ）の左右方向は方位方向を示し、グラフ３５は方位方向の音場分布を示している。また、図６（ｄ）の左右方向はスライス方向を示し、グラフ３６はスライス方向の音場分布を示している。
【０１００】
図６（ｃ）に示されているように、処理回路１５０は、音場算出機能１５０ｋにより、深さx=xiでの方位方向音場分布に対し、音場の中心を中央にして間隔Δyで５点のフィルタ係数ha1、ha2、ha3、ha4、ha5を、それぞれ位置３２ａ、３２ｂ、３２ｃ、３２ｄ、３２ｅにおける音場分布に基づいて算出する。
【０１０１】
また、図６（ｄ）に示されているように、処理回路１５０は、音場算出機能１５０ｋにより、深さx=xiでのスライス方向音場分布に対し、音場の中心を中央にして間隔Δｚ（または間隔Δz1、Δz2、Δz3、Δz4）で５点のフィルタ係数he1、he2、he3、he4、he5を、それぞれ位置３３ａ、３３ｂ、３３ｃ、３３ｄ、３３ｅにおける音場分布に基づいて算出する。
【０１０２】
また、図６（ｅ）は、フィルタのカーネルサイズを５×５にしたときの、yz面での音場に対応したフィルタ係数ｈ１（ｘ、ｙ）の係数を示している。図６（ｅ）からわかるように、処理回路１５０は、音場算出機能１５０ｋにより、２次元平面のフィルタ係数ｈ１（ｘ、ｙ）の値を、方位方向及びスライス方向それぞれの方向のフィルタ係数の積として算出する。なお、２次元平面のフィルタ係数ｈ１（ｘ、ｙ）の値は、合計値が１になるように正規化される。
【０１０３】
なお、図６の例では、ＲＡＷデータは絶対値を取っているため常に非負であるため、音場分布も非負の値になるような値を用いたが、実施形態はこれに限られない。例えば受信信号P＝I＋jQに対してデコンボリューションを行う場合は、処理回路１５０が音場算出機能１５０ｋによりフィルタ係数を算出するのに使用される音場分布も、例えば符号を持った量として取り扱われても良い。
【０１０４】
図３に戻り、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、信号データまたは画像データである第１のデータに対して、算出された音場に対する伝達関数Ｈ１を用いてデコンボリューションを行い第２のデータを生成する（ステップＳ１２０）。一例として、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、ステップＳ１１０において算出された音場に対する伝達関数Ｈ１に基づいて、式（１０）と同様の計算により、ウィーナフィルタＨ１ｗを算出し、式（１４）により、第１のデータに対して算出したウィーナフィルタＨ１ｗを用いたデコンボリューションを行い第２のデータを生成する。
【０１０５】
ステップＳ１２０の処理が、図４（ｃ）〜（ｆ）に示されている。図４（ｃ）は、図４（ｂ）の第１のデータを平面で切り出したデータを画像で模式的に表わし、上下方向が距離方向、左右方向が方位方向、前後方向がスライス方向を表す。平面５ａ、５ｂ、５ｃ、５ｄは、ｘ方向（距離方向）の位置が異なる平面での画像を表す。画像がぼけた状態になっている。ここで、画像のぼけ方は方位方向の音場の形状やビーム幅wa、スライス方向の音場の形状やビーム幅weに依存する。方位方向の音場の形状やビーム幅wa、スライス方向の音場の形状やビーム幅weは距離（深さ）xに依存するので、画像のぼけ方は各平面で異なる。処理回路１５０は、逆フィルタ処理機能１５０ｌにより、各深さxiでの音場の形状とビーム幅wai、weiを求める。続いて、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、各深さにおいて、音場の形状とビーム幅wai、weiから、式（１２）の音場に対応したフィルタｈ１（x,y）の係数を求める。続いて、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、求めたｈ１(x,y)の係数をフーリエ変換することにより伝達関数Ｈ１を算出し、算出した音場に対する伝達関数Ｈ１に基づいて、式（１０）と同様の計算により、ウィーナフィルタＨ1ｗを算出する。ビーム幅wa及びweiは距離（深さ）に応じて決まることから、ウィーナフィルタＨ1ｗは、それぞれの距離（深さ）ごとに計算される。処理回路１５０は、逆フィルタ処理機能１５０ｌにより、式（１４）で示されているように、図４（ｃ）のデータに対してウィーナフィルタＨ１ｗを用いたデコンボリューションを行って式（１４）の左辺の復元画像F''(u,v)を生成し、さらに復元画像F''(u,v)に対して逆フーリエ変換を行い、図４（ｄ）の６ａ、６ｂ、６ｃ、６ｄの様な分解能が向上した画像を得る。即ち、図４（ｂ）に対して図４（ｅ）に示す分解能が向上したリニアデータ７ａ、７ｂ、７ｃ、７ｄ、７ｅを得る。処理回路１５０は、さらに対数変換機能１５０ｄによりデータを対数圧縮して、図４（ｆ）に示されているように、分解能が向上した３次元ＲＡＷデータ８ａ、８ｂ、８ｃ、８ｄ、８ｅを得る。
【０１０６】
また、処理回路１５０は、画像生成機能１５０ｊにより、第２のデータからボリューム画像、ＭＰＲ画像または２次元画像のうち少なくとも一つの画像を生成してもよい。かかる場合、制御機能１５０ｂを有する処理回路１５０は、画像生成機能１５０ｊにより生成されたそれらの画像をディスプレイ１３５に表示させてもよい。
【０１０７】
図７〜図９を用いて、実施形態に係る逆フィルタ処理機能１５０ｌによる処理の効果について説明する。図７は、比較例に係る超音波診断装置により得られたＢモードＭＰＲ画像の一例である。それに対して、図８は、図７のＢモードＭＰＲ画像に対して対数リニア変換を実行してリニアデータを生成したのち、ウィーナフィルタをかけた画像である。図８の画像は、図７の画像と比較して、分解能が改善している。また、図９は、図７のＢモードＭＰＲ画像に対して、ウィーナフィルタではなく、式（７）の単純な逆フィルタをかけた場合を示している。かかる場合、画像の信号値が発散してしまっている。
【０１０８】
なお、実施形態はこれに限られない。
【０１０９】
ステップＳ１２０においてデコンボリューションを行う第１のデータとしては上述の例に限られず、例えば、対数変換機能１５０ｄによる処理後であって信号処理機能１５０ｆによる処理前のデータ（対数変換器出力データ）や、２乗加算機能１５０ｅによる処理後であって対数変換機能１５０ｄによる処理前のデータ（２乗加算器出力データ）に対して、処理回路１５０はデコンボリューションを行ってもよい。
【０１１０】
処理回路１５０が、対数変換器出力データに対してデコンボリューションを行う場合、処理回路１５０は、対数変換器出力データをメモリ１３２に格納する。かかる場合、１ボリューム分のデータが格納された時点で、処理回路１５０は、メモリ１３２からデータを読出し、逆対数変換機能１５０ｃにより逆対数変換を行い、逆フィルタ処理機能１５０ｌにより逆フィルタ処理を行ったのちに、メモリ１３２に再度格納する。その後、処理回路１５０は、メモリ１３２からデータを読出し後段の処理を行う。
【０１１１】
２乗加算器出力データを用いる場合は、２乗加算器出力データを(I2+Q2)1/2とし、処理回路１５０は、２乗加算器出力データを、メモリ１３２に格納する。１ボリューム分のデータが格納された時点で、処理回路１５０は、メモリ１３２からデータを読出し、逆フィルタ処理機能１５０ｌにより逆フィルタ処理を行い、メモリ１３２に再度格納する。その後、処理回路１５０は、メモリ１３２からデータを読出し後段の処理を行う。なお、処理回路１５０は、対数変換機能１５０ｄにより、対数圧縮を、例えば変換後の信号値が20log10(I2+Q2)1/2となるように行う。
【０１１２】
このように、対数変換器データ出力や２乗加算器データ出力を用いてデコンボリューションを行う場合、メモリ１３２を用意する必要があるが、後段の信号処理器１５０ｆの影響を受けずに逆フィルタ処理ができるので精度を向上させることができる。
【０１１３】
また、第１のデータの種類は上述の例に限られず、ラスタに対応したデータであるラスタデータをフレーム毎に持ち、位置の異なる複数フレームで３次元データが構成された３次元ＲＡＷデータや、３次元座標に対応付けられた３次元ボリュームデータ、またはラスタデータをフレーム単位で持つ２次元ＲＡＷデータであっても、実施形態は同様に適用可能である。また、第１のデータの例として、２次元アレープローブで収集を行う場合、ラスタデータはボリューム単位で保持されてもよい。
【０１１４】
また、第１のデータとしては、ＭＰＲ画像データでもよい。この場合、ＭＰＲ画像は２次元データなので、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、２次元逆フィルタまたは１次元逆フィルタを用いてデコンボリューションを行う。また、処理回路１５０は、制御機能１５０ｂにより、逆フィルタ後のデータをディスプレイ１３５に表示してもよい。処理回路１５０がデコンボリューションを行う際、リニアデータに一旦変換するのは他の実施例と同様である。また、逆フィルタを用いる場合、ＲＡＷデータと比べ、ＭＰＲ画像データを用いると、処理負荷を軽減することができる。
【０１１５】
また、第１のデータとして、ＩＱ信号そのものを用いることができる。すなわち、第１のデータは、Ｉ信号データまたはＱ信号データである。この場合、処理回路１５０は、逆フィルタ処理機能１５０ｌにより、I、Qそれぞれの信号に対して逆フィルタをかけ、デコンボリューションを行う。また、処理回路１５０は、デコンボリューション後のI、Q信号をメモリ１３２に格納する。１ボリューム分のデータがメモリ１３２に格納された時点で、処理回路１５０は、取得機能１５０ａにより、メモリ１３２からデータを読出し、逆フィルタ処理機能１５０ｌにより逆フィルタ処理を行い、メモリ１３２に再度格納する。その後、処理回路１５０はメモリ１３２からデータを読出し後段の処理を行う。
【０１１６】
第１のデータとして受信信号を用いるとメモリの消費量が２倍になるが、２乗加算器出力以降の信号は絶対値をとっているのに対し、受信信号は振幅と位相を持った符号付きの信号なので、全体としての精度は向上する。また、Ｂモードデータに対して音場に対応する逆フィルタでデコンボリューションするので、Ｂモード画像の分解能が向上し、診断能が向上する。
【０１１７】
（プログラム）
上述した実施形態の中で示した処理手順に示された指示は、ソフトウェアであるプログラムに基づいて実行されることが可能である。汎用コンピュータが、このプログラムを予め記憶しておき、このプログラムを読み込むことにより、上述した実施形態の超音波診断装置による効果と同様の効果を得ることも可能である。例えば、実施形態に係るデータ処理プログラムは、コンピュータに、超音波プローブが超音波を送受信して被検体の走査を行うことにより得られた第１のデータに対して、走査の対象となる場所の音場を示す情報に基づいて得られた伝達関数を用いてデコンボリューションを行い第２のデータを生成する処理を実行させる。
【０１１８】
上述した実施形態で記述された指示は、コンピューターに実行させることのできるプログラムとして、磁気ディスク（フレキシブルディスク、ハードディスクなど）、光ディスク（ＣＤ−ＲＯＭ、ＣＤ−Ｒ、ＣＤ−ＲＷ、ＤＶＤ−ＲＯＭ、ＤＶＤ±Ｒ、ＤＶＤ±ＲＷなど）、半導体メモリ、又はこれに類する記録媒体に記録される。コンピューター又は組み込みシステムが読み取り可能な記憶媒体であれば、その記憶形式は何れの形態であってもよい。コンピューターは、この記録媒体からプログラムを読み込み、このプログラムに基づいてプログラムに記述されている指示をＣＰＵで実行させれば、上述した実施形態の超音波診断装置と同様の動作を実現することができる。また、コンピューターがプログラムを取得する場合又は読み込む場合は、ネットワークを通じて取得又は読み込んでもよい。
【０１１９】
また、記憶媒体からコンピューターや組み込みシステムにインストールされたプログラムの指示に基づきコンピューター上で稼働しているＯＳ（Operating  System）や、データベース管理ソフト、ネットワーク等のＭＷ（Middleware）等が、上述した実施形態を実現するための各処理の一部を実行してもよい。更に、記憶媒体は、コンピュータあるいは組み込みシステムと独立した媒体に限らず、ＬＡＮ（Local  Area  Network）やインターネット等により伝達されたプログラムをダウンロードして記憶又は一時記憶した記憶媒体も含まれる。また、記憶媒体は１つに限られず、複数の媒体から、上述した実施形態における処理が実行される場合も、実施形態における記憶媒体に含まれ、媒体の構成は何れの構成であってもよい。
【０１２０】
なお、実施形態におけるコンピュータ又は組み込みシステムは、記憶媒体に記憶されたプログラムに基づき、上述した実施形態における各処理を実行するためのものであって、パソコン、マイコン等の１つからなる装置、複数の装置がネットワーク接続されたシステム等の何れの構成であってもよい。また、実施形態におけるコンピュータとは、パソコンに限らず、情報処理機器に含まれる演算処理装置、マイコン等も含み、プログラムによって実施形態における機能を実現することが可能な機器、装置を総称している。
【０１２１】
（データ処理装置）
実施形態において、超音波診断装置にて実行されるデータ処理は、超音波診断装置以外の他の装置で実行されてもよい。例えば、実施形態で説明したデータ処理は、超音波診断装置とは独立に設置されたデータ処理装置や画像処理装置により行われる場合であってもよい。具体的には、取得機能１５０ａや逆フィルタ処理機能１５０ｌ、画像生成機能１５０ｊ等の機能を有するデータ処理装置が、超音波診断装置、又は、ＰＡＣＳのデータベースや、電子カルテシステムのデータベースから超音波画像データ群や受信信号等を受信して上述した画像処理を行なう場合であってもよい。かかるデータ処理装置は、例えば図１で示される超音波診断装置の制御機能１５０ｂ、取得機能１５０ａ、画像生成機能１５０ｊや逆フィルタ処理機能１５０ｌと同様の構成や機能等を備える。かかるデータ処理装置は、例えば、取得機能と、逆フィルタ機能とを備える。取得機能は、超音波プローブが超音波を送受信して被検体の走査を行うことにより得られた第１のデータを取得する。逆フィルタ機能は、第１のデータに対して、走査の対象となる場所の音場を示す情報に基づいて得られた伝達関数を用いてデコンボリューションを行い第２のデータを生成する。
【０１２２】
また、一例として、かかるデータ処理装置（超音波画像処理装置）では、例えば実施形態に係る超音波診断装置がＤＩＣＯＭフォーマットでメモリ１３２内に保存したＲＡＷデータ、または送受信条件データ等を直接に、例えばＵＳＢメモリ等のデータ入出力インタフェースを介して、またはＬＡＮ等のネットワークインタフェースを介して、データ処理装置に入力する。超音波画像処理装置では、超音波診断装置で行った様々なデコンボリューション処理を行う。超音波画像処理装置では、超音波診断装置に比べて処理能力と処理時間に余裕を持たせることができるので、性能の高い逆フィルタや精度の高い音場データ生成方法を用いて、デコンボリューションの精度を向上させることもできる。
【０１２３】
以上のように、実施形態に係る超音波診断装置によれば、画像の分解能を向上させることができる。
【０１２４】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【０１２５】
１  超音波プローブ
１５０ａ  取得機能
１５０ｌ  逆フィルタ処理機能
【図１】
2019037661.tif 000018
【図２】
2019037661.tif 000019
【図３】
2019037661.tif 000020
【図４】
2019037661.tif 000021
【図５】
2019037661.tif 000022
【図６】
2019037661.tif 000023
【図７】
2019037661.tif 000024
【図８】
2019037661.tif 000025
【図９】
2019037661.tif 000026