JP7005206B2 - 超音波診断装置及び超音波イメージングプログラム - Google Patents

Description

本発明の実施形態は、超音波診断装置及び超音波イメージングプログラムに関する。

従来、超音波診断装置において、複数の振動素子が一列に配置された１次元超音波プローブ（「１Ｄアレイプローブ」とも表記する）や、複数の振動素子がマトリクス状に配置された２次元超音波プローブ（「２Ｄアレイプローブ」とも表記する）が利用されている。２Ｄアレイプローブでは、１Ｄアレイプローブと比べて多くのチャンネルを有するため、１Ｄアレイプローブと比べて空間分解能が高い画像を収集したり、プローブ位置を固定したままスキャン断面を回転させたりすることができる。

一般的に、２Ｄアレイプローブのチャンネル数は、装置本体側のチャンネル数よりも多い。このため、複数の振動素子をサブアレイと呼ばれるグループに分割し、収集した信号をサブアレイ単位で合成（遅延加算）することで、２Ｄアレイプローブから装置本体側へ送られるチャンネル数を減らすことが行われている。しかしながら、信号を合成する過程でノイズが発生してしまうことがある。

特開２０１５－０９７６５５号公報 特開２００５－３４２１９４号公報 特開２０１３－０００３５１号公報

本発明が解決しようとする課題は、ノイズを除去することができる超音波診断装置及び超音波イメージングプログラムを提供することである。

実施形態の超音波診断装置は、超音波プローブと、演算部と、画像生成部とを備える。超音波プローブは、第１の位相の第１超音波の送信及び前記第１の位相と実質的に９０度異なる第２の位相の第２超音波の送信のみからなる１セットの超音波送信を、複数の走査線に沿って順々に実行する。演算部は、前記超音波プローブを介して取得された、前記第１超音波に対応する第１エコー信号から前記第２超音波に対応する第２エコー信号を減算することにより減算信号を生成する。画像生成部は、複数の走査線について生成された前記減算信号に基づいて、超音波画像を生成する。

図１は、実施形態に係る超音波診断装置の構成例を示すブロック図である。 図２は、図１に示すＢモード処理回路の構成例を示すブロック図である。 図３Ａは、図１に示す超音波プローブの構成を説明するための図である。 図３Ｂは、図１に示す超音波プローブの構成を説明するための図である。 図４は、実施形態に係る超音波診断装置における撮像法を説明するための図である。 図５は、実施形態に係る超音波診断装置における処理手順を説明するためのフローチャートである。 図６は、実施形態に係る超音波診断装置による効果を説明するための図である。 図７は、その他の実施形態に係る超音波診断装置における処理手順を説明するためのフローチャートである。 図８は、その他の実施形態に係る超音波診断装置の構成例を示すブロック図である。 図９は、その他の実施形態に係る超音波診断装置の構成例を示すブロック図である。

以下、図面を参照して、実施形態に係る超音波診断装置及び超音波イメージングプログラムを説明する。

（実施形態）
図１は、実施形態に係る超音波診断装置１の構成例を示すブロック図である。図１に示すように、実施形態に係る超音波診断装置１は、装置本体１００と、超音波プローブ１０１と、入力装置１０２と、ディスプレイ１０３とを備える。超音波プローブ１０１、入力装置１０２、及びディスプレイ１０３は、それぞれ装置本体１００に接続される。

超音波プローブ１０１は、複数の振動素子（圧電振動子）を有する。超音波プローブ１０１は、被検体Ｐの体表面に接触され、超音波の送受信（超音波走査）を行う。複数の振動素子は、後述する装置本体１００が有する送信回路１１０から供給される駆動信号に基づいて、超音波を発生させる。発生した超音波は、被検体Ｐ内の音響インピーダンスの不整合面で反射され、組織内の散乱体によって散乱された成分等を含む反射波信号（受信エコー）として複数の振動素子にて受信される。超音波プローブ１０１は、複数の振動素子にて受信した反射波信号を、送信回路１１０へ送る。

なお、本実施形態では、超音波プローブ１０１が、マトリクス状（格子状）に配列された複数の振動素子を有する２次元超音波プローブ（「２Ｄアレイプローブ」とも表記する）である場合を説明するが、これに限定されるものではない。例えば、超音波プローブ１０１は、所定方向に１次元で配列された複数の振動素子を有する１次元超音波プローブ（「１Ｄアレイプローブ」とも表記する）であってもよい。

入力装置１０２は、マウス、キーボード、ボタン、パネルスイッチ、タッチコマンドスクリーン、フットスイッチ、トラックボール、ジョイスティック等を有し、超音波診断装置１の操作者からの各種設定要求を受け付け、装置本体１００に対して受け付けた各種設定要求を転送する。

ディスプレイ１０３は、超音波診断装置１の操作者が入力装置１０２を用いて各種設定要求を入力するためのＧＵＩ（Graphical User Interface）を表示したり、装置本体１００において生成された超音波画像データ等を表示したりする。

装置本体１００は、超音波プローブ１０１が受信した反射波信号に基づいて、超音波画像データを生成する装置である。図１に示すように、装置本体１００は、例えば、送信回路１１０と、受信回路１２０と、Ｂモード処理回路１３０と、ドプラ処理回路１４０と、画像生成回路１５０と、画像メモリ１６０と、記憶回路１７０と、制御回路１８０とを有する。送信回路１１０、受信回路１２０、Ｂモード処理回路１３０と、ドプラ処理回路１４０、画像生成回路１５０、画像メモリ１６０、記憶回路１７０、及び制御回路１８０は、互いに通信可能に接続される。

送信回路１１０は、パルサ回路等を有する。パルサ回路は、所定のレート周波数（ＰＲＦ：Pulse Repetition Frequency）で、送信超音波を形成するためのレートパルスを繰り返し発生し、発生したレートパルスを超音波プローブ１０１に出力する。

また、送信回路１１０は、制御回路１８０による制御を受けて、後述する送受信制御回路３１に、後述するパルサ３３が出力する駆動信号の振幅の値を出力する。また、送信回路１１０は、制御回路１８０による制御を受けて、後述する遅延加算回路３７における反射波信号の遅延量を出力する。

受信回路１２０は、Ａ／Ｄ変換器及び受信ビームフォーマを有する。受信回路１２０は、超音波プローブ１０１から出力された反射波信号を受信すると、まず、Ａ／Ｄ変換器が、反射波信号をデジタルデータに変換し、受信ビームフォーマが、これらの各チャンネルからのデジタルデータに対し整相加算処理を行って反射波データを生成し、生成した反射波データをＢモード処理回路１３０及びドプラ処理回路１４０に送信する。

Ｂモード処理回路１３０は、受信回路１２０から出力された反射波データを受信し、受信した反射波データに対して対数増幅、包絡線検波処理等を行って、信号強度が輝度の明るさで表現されるデータ（Ｂモードデータ）を生成する。

図２は、図１に示すＢモード処理回路の構成例を示すブロック図である。図２に示すように、Ｂモード処理回路１３０は、Ｂモードデータ生成機能１３１と、演算機能１３２と、フィルタ機能１３３とを備える。Ｂモードデータ生成機能１３１は、反射波データに対して、対数増幅処理、包絡線検波処理、対数圧縮処理等を行って、Ｂモードデータを生成する。通常のＢモード撮影が行われている場合には、演算機能１３２及びフィルタ機能１３３による処理は実行されず、Ｂモードデータ生成機能１３１は、受信回路１２０から受信した反射波データからＢモードデータを生成する。一方、演算機能１３２及びフィルタ機能１３３による処理が実行される場合には、Ｂモードデータ生成機能１３１は、フィルタ機能１３３が出力したデータからＢモードデータを生成する。なお、演算機能１３２及びフィルタ機能１３３の処理については、後述する。

ドプラ処理回路１４０は、受信回路１２０から出力された反射波データを受信し、受信した反射波データから速度情報を周波数解析し、ドプラ効果による血流や組織、造影剤エコー成分を抽出し、平均速度、分散、パワー等の移動体情報を多点について抽出したデータ（ドプラデータ）を生成する。

画像生成回路１５０は、Ｂモード処理回路１３０及びドプラ処理回路１４０が生成したデータから超音波画像データを生成する。画像生成回路１５０は、Ｂモード処理回路１３０が生成したＢモードデータから反射波の強度を輝度で表したＢモード画像データを生成する。また、画像生成回路１５０は、ドプラ処理回路１４０が生成したドプラデータから移動体情報を表すドプラ画像データを生成する。ドプラ画像データは、速度画像データ、分散画像データ、パワー画像データ、又は、これらを組み合わせた画像データである。

ここで、画像生成回路１５０は、一般的には、超音波走査の走査線信号列を、テレビ等に代表されるビデオフォーマットの走査線信号列に変換（スキャンコンバート）し、表示用の超音波画像データを生成する。具体的には、画像生成回路１５０は、超音波プローブ１０１による超音波の走査形態に応じて座標変換を行うことで、表示用の超音波画像データを生成する。また、画像生成回路１５０は、スキャンコンバート以外に、種々の画像処理として、例えば、スキャンコンバート後の複数の画像フレームを用いて、輝度の平均値画像を再生成する画像処理（平滑化処理）や、画像内で微分フィルタを用いる画像処理（エッジ強調処理）等を行う。また、画像生成回路１５０は、超音波画像データに、種々のパラメータの文字情報、目盛り、ボディーマーク等を合成する。

画像メモリ１６０は、画像生成回路１５０が生成した画像データ（Ｂモード画像データ、ドプラ画像データ等）を記憶するメモリである。また、画像メモリ１６０は、Ｂモード処理回路１３０及びドプラ処理回路１４０が生成したデータを記憶することも可能である。画像メモリ１６０が記憶するＢモードデータ及びドプラデータは、例えば、操作者が呼び出すことが可能となっており、画像生成回路１５０を経由して表示用の超音波画像データとなる。

記憶回路１７０は、超音波送受信、画像処理及び表示処理などを行うための制御プログラムや、診断情報（例えば、患者ＩＤ、医師の所見等）や、診断プロトコルや各種ボディーマーク等の各種データを記憶する。また、記憶回路１７０は、必要に応じて、画像メモリ１６０が記憶する画像データの保管等にも使用される。また、記憶回路１７０が記憶するデータは、図示しないインタフェース部を介して、外部装置へ転送することができる。

制御回路１８０は、超音波診断装置１の処理全体を制御する。具体的には、制御回路１８０は、入力装置１０２を介して操作者から入力された各種設定要求や、記憶回路１７０から読み込んだ各種制御プログラムおよび各種データに基づき、送信回路１１０、受信回路１２０、Ｂモード処理回路１３０、ドプラ処理回路１４０、画像生成回路１５０等の処理を制御する。また、制御回路１８０は、画像メモリ１６０が記憶する超音波画像データをディスプレイ１０３に表示させる。

なお、装置本体１００に内蔵される送信回路１１０、受信回路１２０、Ｂモード処理回路１３０、ドプラ処理回路１４０、画像生成回路１５０、及び制御回路１８０は、プロセッサ（ＣＰＵ（Central Processing Unit）、ＭＰＵ（Micro-Processing Unit）、集積回路等）のハードウェアにより構成される。

また、各プロセッサが実行する各処理機能は、例えば、コンピュータによって実行可能なプログラムの形態で記憶回路１７０に記録されている。つまり、各プロセッサは、各プログラムを記憶回路１７０から読み出し、実行することで各プログラムに対応する機能を実現する。例えば、図２に示したＢモードデータ生成機能１３１、演算機能１３２及びフィルタ機能１３３の各処理機能は、Ｂモード処理回路１３０が各処理機能に対応するプログラムを記憶回路１７０から読み出し実行することで、実現される機能である。換言すると、各プログラムを読み出した状態のＢモード処理回路１３０は、図２のＢモード処理回路１３０内に示された各機能を有することとなる。

なお、図１及び図２においては、単一のＢモード処理回路１３０にて、Ｂモードデータ生成機能１３１、演算機能１３２及びフィルタ機能１３３にて行われる処理機能が実現されるものとして説明するが、複数の独立したプロセッサを組み合わせて処理回路を構成し、各プロセッサがプログラムを実行することにより各機能を実現するものとしても構わない。

ここで、本実施形態において装置本体１００に接続される超音波プローブ１０１について説明する。本実施形態では、超音波プローブ１０１は、複数の振動素子がマトリクス状に２次元で配置された２Ｄアレイプローブである。

図３Ａ及び図３Ｂは、図１に示す超音波プローブ１０１の構成を説明するための図である。図３Ａには、２次元アレイプローブである超音波プローブ１０１が有する複数の振動素子の配置を例示する。図３Ｂには、超音波プローブ１０１の内部の構成例を例示する。なお、図３Ａ及び図３Ｂの内容はあくまで一例であり、従来の如何なる２次元アレイプローブが適用されてもよい。

図３Ａに示すように、超音波プローブ１０１は、ラテラル方向及びエレベーション方向に複数の振動素子２０（「振動素子群」とも表記する）が並んだ２Ｄアレイプローブである。２Ｄアレイプローブにおいて、超音波プローブ１０１が有する全ての振動素子２０は、メインアレイ２１と呼ばれる。メインアレイ２１は、ラテラル方向及びエレベーション方向に、複数のサブアレイ２２に分割される。サブアレイ２２とは、例えば、メインアレイ２１に対応する複数の振動素子２０が所定数の振動素子２０ごとに分割されたグループ分けを表す。図３Ａに示す例では、サブアレイ２２は、ラテラル方向に並んだ５個の振動素子２０が、エレベーション方向に５列並んだ２５個の振動素子２０を有する。言い換えると、サブアレイ２２は、「５×５」のマトリクス状に配置された２５個の振動素子２０で構成される。なお、図３Ａの例では、一つの振動素子のみに符号「２０」を付し、他の振動素子については符号「２０」を付していない。同様に、一つのサブアレイのみに符号「２２」を付し、他のサブアレイについては符号「２２」を付していない。

図３Ｂに示すように、超音波プローブ１０１は、特定用途向け集積回路（Application Specific Integrated Circuit：ＡＳＩＣ）３０を備える。また、超音波プローブ１０１は、ＡＳＩＣ３０上に、送受信制御回路３１、送信遅延回路３２、パルサ３３、送受信スイッチ（Ｔ／Ｒ ＳＷ（Transmission/Reception Switch））３４、低雑音増幅器（ＬＮＡ（Low Noise Amplifier））３５、タイムゲインコントローラ（ＴＧＣ（Time Gain Controller）３６、及び遅延加算回路３７の各電子回路を備える。ここで、超音波プローブ１０１において、１つの振動素子２０に対して１つのチャンネルが割り当てられる。超音波プローブ１０１は、チャンネルごとに、送信遅延回路３２、パルサ３３、送受信スイッチ３４、低雑音増幅器３５、及びタイムゲインコントローラ３６を有する。また、超音波プローブ１０１は、１つのサブアレイ２２に対して、送受信制御回路３１及び遅延加算回路３７を有する。つまり、超音波プローブ１０１は、メインアレイ２１に対して、図３Ｂに示す各電子回路を全てのサブアレイ２２（図３Ａの例では４２個のサブアレイ２２）分だけ備える。なお、ＡＳＩＣ３０は、超音波プローブ１０１に対して１つ又は数個備えられる。

送受信制御回路３１は、超音波の送受信を制御する。例えば、送受信制御回路３１は、送信回路１１０から出力されたレートパルスを受信して、受信したレートパルスを送信遅延回路３２へ送る。また、送受信制御回路３１は、送信回路１１０から出力された反射波信号の遅延時間を受信して、受信した反射波信号の遅延時間を後述する遅延加算回路３７に設定する。

送信遅延回路３２は、振動素子２０から発生される超音波をビーム状に集束して送信指向性を決定するために必要な振動素子２０ごとの遅延時間を、装置本体１００から供給されたレートパルスに対して与える。例えば、送信遅延回路３２は、送受信制御回路３１から出力されたレートパルスに対し、チャンネルごとに設定された遅延時間を与えてパルサ３３へ出力する。なお、レートパルスに与えられる遅延時間は、送受信制御回路３１により制御される。

パルサ３３は、所定の振幅値の駆動信号を発生させる。例えば、パルサ３３は、送信遅延回路３２から出力されたレートパルスに基づくタイミングで駆動信号を発生させ、振動素子２０へ出力する。なお、発生される駆動信号の振幅値は、送受信制御回路３１により制御される。

送受信スイッチ３４は、振動素子２０の接続先を、パルサ３３及び低雑音増幅器３５のいずれか一方に選択的に切り替える。送受信スイッチ３４がパルサ３３に接続される場合、送受信スイッチ３４は、パルサ３３から出力された駆動信号を振動素子２０に送信する。一方、送受信スイッチ３４が低雑音増幅器３５に接続される場合、送受信スイッチ３４は、振動素子２０から送信された反射波信号を低雑音増幅器３５に出力する。

ここで、パルサ３３に駆動信号を発生させるレートパルスは、送信回路１１０に由来する。また、低雑音増幅器３５へ出力される反射波信号は、後述するように、受信回路１２０に受信される。すなわち、送受信スイッチ３０は、超音波プローブ１０１に含まれる振動素子２０の接続先を、送信回路１１０及び受信回路１２０を含む選択肢の中から選択的に切り替える。なお、送受信スイッチ３４は、スイッチング回路の一例である。

低雑音増幅器３５は、送受信スイッチ３４を介して、振動素子２０から反射波信号を受信すると、予め設定されたゲインによって受信した反射波信号を増幅し、増幅した反射波信号をタイムゲインコントローラ３６へ出力する。

タイムゲインコントローラ３６は、低雑音増幅器３５から送信された反射波信号を受信すると、反射波信号を増幅する。そして、タイムゲインコントローラ３６は、増幅した反射波信号を遅延加算回路３７へ出力する。

遅延加算回路３７は、タイムゲインコントローラ３６から出力された各チャンネルの反射波信号を受信すると、各チャンネルの反射波信号に対して、受信指向性を決定するのに必要な遅延量を与える遅延処理を実行する。そして、遅延加算回路３７は、遅延処理後の各チャンネルの反射波信号を加算する加算処理を実行し、加算処理後の反射波信号を装置本体１００の受信回路１２０に出力する。この加算処理は、サブアレイ２２内のチャンネルに対して行われる。つまり、遅延加算回路３７は、サブアレイ２２内の各チャンネルの反射波信号をサブアレイ２２ごとに合成（遅延加算処理）する。

このように、実施形態に係る超音波プローブ１０１は、マトリクス状に配列された複数の振動素子２０を有する２Ｄアレイプローブである。超音波プローブ１０１は、複数の振動素子２０で構成されるサブアレイ２２と、複数のサブアレイ２２で構成されるメインアレイ２１と、サブアレイ２２毎に反射波信号の遅延加算処理を実行する遅延加算回路３７とを備える。

ところで、上述した２Ｄアレイプローブを用いる場合、決まったタイミングで生じるノイズが受信エコーに含まれてしまうことがある。例えば、遅延加算回路３７の処理において、遅延メモリに記憶された各チャンネルの遅延時間を繰り返し更新させるために、アナログスイッチが用いられる。このアナログスイッチは、スイッチの切り替え時にノイズ（スイッチノイズ）を発生させる。アナログスイッチは、遅延時間を更新させるために所定の更新レートで切り替わるため、アナログスイッチに由来するノイズも更新レートに従って生じる結果、周期的に生じるノイズ（「周期ノイズ」とも表記する）として受信エコーに含まれてしまうことがある。

また、固定ノイズとしては、疑似送信によるノイズ（「固定ノイズ」とも表記する）も受信エコーに含まれてしまうことがある。ここで、疑似送信とは、本来送信されるべきでないタイミングで超音波が送信されてしまうことである。例えば、送受信スイッチ３４における送受信切り替え時のスイッチノイズが振動素子２０に伝わると、本来送信されるべきでないタイミングで微弱な超音波が送信（疑似送信）される。疑似送信は本来の超音波同様に振動素子より送信波として出ていくため、被検体Ｐ内で反射される。反射された受信エコーが振動素子２０に受信されることにより、本来の受信エコーにノイズとして含まれてしまうことがある。具体的には、疑似送信による固定ノイズは、送信側から受信側へ切り替わる場合と、受信側から送信側へ切り替わる場合とに発生するため、本来の受信エコーの前後に出現して信号を三重化させてしまう。

このように、２Ｄアレイプローブを用いる場合には、周期ノイズや疑似送信による固定ノイズが受信エコーに含まれてしまうことがある。周期ノイズや固定ノイズは、大きさが決まっているため、本来の送信超音波の振幅が低い場合に現れやすいと言える。例えば、マイクロバブルを造影剤として用いるコントラストハーモニックイメージング（Contrast Harmonic Imaging：ＣＨＩ）法においては、送信超音波の振幅が低めに設定されるため、固定ノイズが現れやすい。

そこで、実施形態に係る超音波診断装置１は、周期ノイズや固定ノイズを除去するために、ＰＳ（Pulse Subtraction）－ＴＨＩ（Tissue Harmonic Imaging）法を応用した以下の撮像法を実行する。

すなわち、実施形態に係る超音波診断装置１において、超音波プローブ１０１は、第１の位相の第１超音波の送信及び第１の位相と９０度異なる第２の位相の第２超音波の送信のみからなる１セットの超音波送信を、複数の走査線に沿って順々に実行する。演算機能１３２は、超音波プローブ１０１を介して取得された、第１超音波に対応する第１エコー信号から第２超音波に対応する第２エコー信号を減算することにより減算信号を生成する。画像生成回路１５０は、複数の走査線について生成された減算信号に基づいて、超音波画像を生成する。以下、超音波診断装置１における撮像法について説明する。

なお、以下の実施形態では、２Ｄアレイプローブを用いる場合に生じる周期ノイズや固定ノイズを除去する処理を説明するが、実施形態はこれに限定されるものではない。例えば、上述した疑似送信によるノイズは、２Ｄアレイプローブ以外の超音波プローブ１０１（例えば、１Ｄアレイプローブ）であっても生じるノイズである。したがって、本実施形態は、２Ｄアレイプローブ以外の超音波プローブ１０１を用いる場合にも、周期ノイズや固定ノイズを除去する効果を奏するものである。

また、以下の実施形態では、２次高調波成分を抽出するＴＨＩ法を実行する場合を説明するが、実施形態はこれに限定されるものではない。例えば、以下の実施形態は、３次以上の高調波成分を抽出するＴＨＩ法に適用されても良い。また、以下の実施形態は、ＴＨＩ法に限らず、ＣＨＩ法に適用されても良いし、基本波を含む撮像法に広く適用されてもよい。

図４は、実施形態に係る超音波診断装置１における撮像法を説明するための図である。図４には、２回の超音波の送信からなる１セットの超音波送信により収集された２回分のエコー信号を減算する処理を例示する。

図４において、左側の波形は、１波目の送信により得られた１波目のエコー信号（第１エコー信号）の波形を表す。また、中央の波形は、２波目の送信により得られた２波目のエコー信号（第２エコー信号）の波形を表す。また、右側の波形は、減算信号の波形を表す。図４に示す例では、左側の波形のエコー信号から中央の波形のエコー信号を減算することで、右側の波形の減算信号を得る場合を例示する。また、各波形において、横軸は時間に対応し、縦軸は振幅に対応する。

例えば、送信回路１１０は、制御回路１８０により設定されたスキャンシーケンスに従って、超音波プローブ１０１に超音波走査を実行させる。具体的には、送信回路１１０は、第１の位相の第１超音波の送信及び第１の位相と実質的に９０度異なる第２の位相の第２超音波の送信のみからなる１セットの超音波送信を、複数の走査線に沿って順々に実行する超音波走査を、超音波プローブ１０１に実行させる。

つまり、送信回路１１０は、走査範囲に含まれる各走査線に対して、２回の超音波の送信からなる１セットの超音波送信を超音波プローブ１０１に実行させる。例えば、送信回路１１０は、１波目（１回目）の超音波（第１超音波）の送信として、正弦波の超音波を送信させ、２波目（２回目）の超音波（第２超音波）の送信として、１波目の超音波に対して位相を９０°回転させた超音波（余弦波）を送信させる。これにより、超音波プローブ１０１は、１波目の超音波の送信として、正弦波の超音波を送信し、２波目の超音波の送信として、１波目の超音波に対して位相を９０°回転させた超音波を送信する１セットの超音波送信を、各走査線に対して順々に実行する。

そして、受信回路１２０は、各走査線に対して、２回のエコー信号を含む反射波データを生成する。例えば、受信回路１２０は、１波目の送信超音波により得られる１波目のエコー信号と、２波目の送信超音波により得られる２波目のエコー信号とを含む反射波データを生成する（図４参照）。受信回路１２０は、生成した反射波データを、演算機能１３２へ出力する。

演算機能１３２は、超音波プローブ１０１を介して取得された、第１超音波に対応する第１エコー信号から第２超音波に対応する第２エコー信号を減算することにより減算信号を生成する。なお、演算機能１３２は、演算部の一例である。

ここで、上述した周期ノイズや固定ノイズは、超音波を送受信するごとに、決まったタイミングで含まれるノイズである。つまり、周期ノイズや固定ノイズは、１波目のエコー信号にも２波目のエコー信号にも、決まったタイミングで含まれている。そこで、演算機能１３２は、図４に示すように、１波目のエコー信号から２波目のエコー信号を１回減算する。これにより、演算機能１３２は、１波目のエコー信号及び２波目のエコー信号において決まったタイミングで含まれる周期ノイズや固定ノイズを除去することができる。そして、演算機能１３２は、１波目のエコー信号から２波目のエコー信号を１回減算することにより、減算信号を生成する。演算機能１３２は、複数の走査線それぞれについて減算信号を生成し、生成した各走査線の減算信号をフィルタ機能１３３へ出力する。

フィルタ機能１３３は、演算機能１３２から出力された減算信号に対して、ハイパスフィルタ又はバンドパスフィルタを適用する。例えば、フィルタ機能１３３は、各走査線の減算信号に含まれる基本波成分を取り除き、２次高調波成分を抽出するフィルタを適用する。なお、フィルタ機能１３３は、フィルタ部の一例である。

図４に示す例では、１波目のエコー信号に含まれる基本波成分の波形４０は、ｓｉｎθで表される。この場合、１波目のエコー信号に含まれる２次高調波成分の波形４１は、ｓｉｎ２θで表される。また、２波目のエコー信号に含まれる基本波成分の波形４２は、１波目の位相が９０°回転された位相であるため、ｃｏｓθで表される。この場合、２波目のエコー信号に含まれる２次高調波成分の波形４３は、－ｓｉｎ２θで表される。

ここで、２次高調波成分に着目すると、１波目の２次高調波成分の波形４１から２波目の２次高調波成分の波形４３が減算された結果、減算信号に含まれる２次高調波成分の波形４４の振幅は、減算前の振幅と比べて２倍になる。そこで、フィルタ機能１３３は、例えば、各走査線の減算信号に対してハイパスフィルタ（若しくはバンドパスフィルタ）を適用することで、波形４４で表される２次高調波成分をそれぞれ抽出する。フィルタ機能１３３は、抽出した２次高調波成分を含む各走査線の減算信号を、Ｂモードデータ生成機能１３１へ出力する。

これにより、Ｂモードデータ生成機能１３１は、各走査線の減算信号を用いてＢモードデータを生成する。そして、画像生成回路１５０は、Ｂモードデータ生成機能１３１により生成されたＢモードデータから超音波画像データを生成する。すなわち、画像生成回路１５０は、複数の走査線について生成された減算信号に基づいて、超音波画像データを生成する。具体的には、画像生成回路１５０は、フィルタ機能１３３により抽出された２次高調波成分を用いて、超音波画像データを生成する。

このように、実施形態に係る超音波診断装置１は、１回の減算により、固定ノイズを除去するとともに、２次高調波成分が２倍（６ｄＢ）となった減算信号を生成する。なお、図４の内容はあくまで一例であり、これに限定されるものではない。例えば、１波目として送信される超音波は、正弦波に限定されるものではなく、任意の初期位相を有する超音波であってよい。

また、１波目の送信超音波の位相と２波目の送信超音波の位相の差は、９０°に限定されるものではなく、例えば、２７０°であっても良い。言い換えると、超音波プローブ１０１は、第１の位相の第１超音波の送信及び第１の位相と実質的に９０度異なる第２の位相の第２超音波の送信のみからなる１セットの超音波送信を、複数の走査線に沿って順々に実行する。

図５は、実施形態に係る超音波診断装置１における処理手順を説明するためのフローチャートである。図５に示す処理手順は、例えば、撮像を開始する旨の指示を操作者から受け付けた場合に、開始される。

図５に示すように、例えば、入力装置１０２が撮像を開始する旨の指示を操作者から受け付けると（ステップＳ１０１肯定）、制御回路１８０は、ステップＳ１０２以降の処理を開始する。なお、撮像を開始する旨の指示を受け付けるまで（ステップＳ１０１否定）、制御回路１８０は、以下の処理を開始せず、待機状態である。

撮像が開始されると、超音波プローブ１０１は、第１の位相の第１超音波の送信及び第１の位相と９０度異なる第２の位相の第２超音波の送信のみからなる１セットの超音波送信を、複数の走査線に沿って順々に実行する（ステップＳ１０２）。

続いて、演算機能１３２は、超音波プローブ１０１を介して取得された、第１超音波に対応する第１エコー信号から第２超音波に対応する第２エコー信号を減算して、減算信号を生成する（ステップＳ１０３）。例えば、演算機能１３２は、１波目のエコー信号から２波目のエコー信号を１回減算することで、固定ノイズを除去した減算信号を生成する。

そして、フィルタ機能１３３は、複数の走査線ごとの減算信号に対して、ハイパスフィルタを適用する（ステップＳ１０４）。例えば、フィルタ機能１３３は、各走査線の減算信号に含まれる基本波成分をハイパスフィルタにより取り除き、２次高調波成分を抽出する。

そして、画像生成回路１５０は、複数の走査線について生成された減算信号に基づいて、超音波画像を生成する（ステップＳ１０５）。例えば、画像生成回路１５０は、Ｂモードデータ生成機能１３１によって減算信号から生成されたＢモードデータを用いて、超音波画像データを生成する。

ここで、略リアルタイムで撮像を行う場合には、超音波診断装置１は、ステップＳ１０２～ステップＳ１０５の処理を繰り返し実行することで、略リアルタイムの超音波画像データを生成し、表示する。そして、撮像を終了する旨の操作を操作者から受け付けた場合に、超音波診断装置１は、図５の処理手順を終了する。

なお、図５の内容はあくまで一例であり、これに限定されるものではない。例えば、ステップＳ１０４の処理は、必ずしも実行されなくてもよい。この場合、２次高調波成分の抽出が行われないので、基本波成分を含む減算信号が超音波画像データの生成に利用されることとなる。

また、例えば、ステップＳ１０４の処理は、ステップＳ１０３の処理の前に実行されてもよい。この場合、フィルタ機能１３３は、第１エコー信号及び第２エコー信号を含む反射波データに対してハイパスフィルタを適用して、第１エコー信号及び第２エコー信号に含まれる基本波成分を取り除き、２次高調波成分を抽出する。そして、演算機能１３２は、抽出された第１エコー信号及び第２エコー信号の２次高調波成分を用いて、減算信号を生成する。

このように、実施形態に係る超音波診断装置１は、ＰＳ－ＴＨＩ法を応用した撮像法を実行することで、ノイズを除去することができる。

図６は、実施形態に係る超音波診断装置１による効果を説明するための図である。図６の左図には、超音波診断装置１の撮像法を適用しない場合の超音波画像５０を例示し、図６の右図には、超音波診断装置１の撮像法を適用した場合の超音波画像５１を例示する。なお、超音波画像５０及び超音波画像５１は、ノイズレベルを同程度に調整した場合の画像である。

ここで、超音波診断装置１の撮像法を適用した場合の超音波画像５１において、減算信号に含まれる２次高調波成分の振幅は、「ｓｉｎ２θ－（－ｓｉｎ２θ）＝２ｓｉｎ２θ」となり、減算前の振幅と比べて２倍になる。この結果、図６に示すように、ノイズレベルを同程度に調整すると、超音波画像５１は、超音波画像５０と比べて明るく描出される。つまり、超音波診断装置１の撮像法により、Ｓ／Ｎ（Signal to Noise）比が６ｄＢ（２０×ｌｏｇ２＝６）改善したことがわかる。また、図示しないが、超音波画像５０及び超音波画像５１の信号レベルを同程度に調整すると、超音波画像５１では超音波画像５０よりも実際に６ｄＢ少ないゲインで画像が得られる。すなわち、超音波診断装置１の撮像法により、Ｓ／Ｎ（Signal to Noise）比が６ｄＢ改善したと言える。

（その他の実施形態）
上述した実施形態以外にも、種々の異なる形態にて実施されてもよい。

（振幅変調）
例えば、上記の実施形態では、位相変調を行った２回の超音波送信からなる１セットの超音波走査を行う場合を説明したが、実施形態はこれに限定されるものではない。例えば、振幅変調を行った２回の超音波送信からなる１セットの超音波走査を行った場合にも、同様の処理を実行することができる。

図７は、その他の実施形態に係る超音波診断装置１における処理手順を説明するためのフローチャートである。図７に示す処理手順は、例えば、撮像を開始する旨の指示を操作者から受け付けた場合に、開始される。

図７に示すように、例えば、入力装置１０２が撮像を開始する旨の指示を操作者から受け付けると（ステップＳ２０１肯定）、制御回路１８０は、ステップＳ２０２以降の処理を開始する。なお、撮像を開始する旨の指示を受け付けるまで（ステップＳ２０１否定）、制御回路１８０は、以下の処理を開始せず、待機状態である。

撮像が開始されると、超音波プローブ１０１は、第１の振幅の第１超音波と第１の振幅より小さい第２超音波のみからなる１セットの送信を、複数の走査線に沿って順々に実行する（ステップＳ２０２）。例えば、送信回路１１０は、振幅変調を行って、第１超音波、及び、第１超音波の振幅より小さい振幅の第２超音波のみからなる１セットの超音波送信を各走査線に対して実行する超音波走査を、超音波プローブ１０１に実行させる。

続いて、演算機能１３２は、第１超音波に対応する第１エコー信号の強度から第２超音波に対応する第２エコー信号の強度の相対的関係を維持したまま減算して、減算信号を生成する（ステップＳ２０３）。例えば、演算機能１３２は、第１エコー信号の強度及び第２エコー信号の強度のうちいずれか一方を変調させたり、両者を異なる倍率で変調させたりする通常の振幅変調法とは異なり、第１エコー信号の強度及び第２エコー信号の強度の相対的関係を維持したまま減算する。つまり、第１エコー信号の強度及び第２エコー信号の強度を変調させずに減算するか、変調させる場合には両者を同じ倍率で変調させた上で減算する。これにより、演算機能１３２は、固定ノイズを除去した減算信号を生成する。

そして、フィルタ機能１３３は、複数の走査線ごとの減算信号に対して、ハイパスフィルタを適用する（ステップＳ２０４）。例えば、フィルタ機能１３３は、各走査線の減算信号に含まれる基本波成分をハイパスフィルタにより取り除き、２次高調波成分を抽出する。

そして、画像生成回路１５０は、複数の走査線について生成された減算信号に基づいて、超音波画像を生成する（ステップＳ２０５）。例えば、画像生成回路１５０は、Ｂモードデータ生成機能１３１によって減算信号から生成されたＢモードデータを用いて、超音波画像データを生成する。

なお、図７の内容はあくまで一例であり、これに限定されるものではない。例えば、ステップＳ２０４の処理は、必ずしも実行されなくてもよい。また、例えば、ステップＳ２０４の処理は、ステップＳ２０３の処理の前に実行されてもよい。

（３次以上の高調波成分の利用）
また、例えば、上記の実施形態では、２次高調波成分を用いた撮像法を説明したが、これに限らず、３次以上の高調波成分が利用されても良い。

例えば、超音波プローブ１０１は、１波目の超音波の送信と、１波目の超音波に対して位相を６０°回転させた２波目の超音波の送信とを１セットとする超音波送信を、各走査線に対して順々に実行する。この場合、１波目のエコー信号に含まれる３次高調波成分の波形は、ｓｉｎ３θで表される。また、２波目のエコー信号に含まれる３次高調波成分の波形は、１波目の位相が６０°回転された位相であるため、－ｓｉｎ３θで表される。

そして、演算機能１３２は、１波目のエコー信号から２波目のエコー信号を減算することにより減算信号を生成する。この場合、３次高調波成分の振幅は、「ｓｉｎ３θ－（－ｓｉｎ３θ）＝２ｓｉｎ３θ」となり、減算前の振幅と比べて２倍になる。これにより、超音波診断装置１は、固定ノイズを除去しつつ、高調波成分の感度を得ることができる。すなわち、超音波診断装置１において、フィルタ機能１３３は、少なくとも２次以上の高調波成分を抽出する。

（ホワイトノイズによる減算処理）
また、例えば、上記の実施形態では、２回の超音波の送信を１セットの超音波送信として実行する場合を説明したが、これに限定されるものではない。例えば、２回目の超音波送信を行わずにホワイトノイズを収集し、減算処理を行っても良い。

すなわち、超音波プローブ１０１は、超音波の送信を行った後に、送信超音波によるエコー信号を受信する。次に、超音波プローブ１０１は、超音波の送信を行わずに、エコー信号を受信する処理と同様の受信処理を行って、ホワイトノイズのノイズ信号を受信する。このノイズ信号には、通常のホワイトノイズも含まれるものの、２Ｄアレイプローブにおける遅延加算処理により、周期ノイズも含まれることとなる。

そして、演算機能１３２は、超音波プローブ１０１により生成されたエコー信号から、ノイズ信号を減算して、減算信号を生成する。これにより、超音波診断装置１は、エコー信号から周期ノイズを除去することができる。

（高機能型超音波プローブを用いた超音波診断装置への適用）
近年、超音波の送受信に関する主要な機能を超音波プローブの筐体内に組み込み、この超音波プローブ（以下、「高機能型超音波プローブ」と称する）をパーソナルコンピュータやタブレット端末等の汎用的な情報処理装置に接続することで、超音波診断装置を実現する技術が知られている。上述した実施形態は、高機能型超音波プローブを用いた超音波診断装置にも適用することができる。

すなわち、高機能型超音波プローブは、超音波の送受信に関する主要な機能を実現するための各種回路（例えば、送信回路１１０及び受信回路１２０）を備えている。このため、これらの回路の駆動により周期ノイズや固定ノイズが生じる可能性が高い。そこで、上述した実施形態に係る構成を適用することで、高機能型超音波プローブにおいて生じる周期ノイズや固定ノイズを除去することが可能である。

図８は、その他の実施形態に係る超音波診断装置２の構成例を示すブロック図である。図８に示すように、その他の実施形態に係る超音波診断装置２は、高機能型超音波プローブ２００と、情報処理装置３００とを備える。高機能型超音波プローブ２００及び情報処理装置３００は、有線通信若しくは無線通信により接続される。

高機能型超音波プローブ２００は、複数の振動素子２０と、送受信回路２１０とを備える。この送受信回路２１０は、送受信制御回路３１、送信遅延回路３２、パルサ３３、送受信スイッチ３４、低雑音増幅器３５、タイムゲインコントローラ３６、遅延加算回路３７、送信回路１１０、及び受信回路１２０を備える。なお、図８に示した送受信制御回路３１、送信遅延回路３２、パルサ３３、送受信スイッチ３４、低雑音増幅器３５、タイムゲインコントローラ３６、及び遅延加算回路３７は、図３Ｂに示した送受信制御回路３１、送信遅延回路３２、パルサ３３、送受信スイッチ３４、低雑音増幅器３５、タイムゲインコントローラ３６、及び遅延加算回路３７と基本的に同様の処理を行うので、同一の符号を付して説明を省略する。また、図８に示した送信回路１１０及び受信回路１２０は、図１に示した送信回路１１０及び受信回路１２０と基本的に同様の処理を行うので、同一の符号を付して説明を省略する。

情報処理装置３００は、例えば、パーソナルコンピュータやタブレット端末等の汎用装置である。情報処理装置３００は、入力回路３０１、ディスプレイ３０２、記憶回路３１０、及び処理回路３２０を備える。入力回路３０１、ディスプレイ３０２、記憶回路３１０、及び処理回路３２０は、相互に通信可能に接続される。

入力回路３０１は、マウス、キーボード、タッチパネル等、操作者からの各種の指示や設定要求を受け付けるための入力装置である。ディスプレイ３０２は、医用画像を表示したり、操作者が入力回路３０１を用いて各種設定要求を入力するためのＧＵＩを表示したりする表示装置である。

記憶回路３１０は、例えば、ＮＡＮＤ（Not AND）型フラッシュメモリやＨＤＤ（Hard Disk Drive）であり、医用画像データやＧＵＩを表示するための各種のプログラムや、当該プログラムによって用いられる情報を記憶する。

処理回路３２０は、情報処理装置３００における処理全体を制御する電子機器（プロセッサ）である。処理回路３２０は、演算機能３２１及び画像生成機能３２２を実行する。演算機能３２１は、図２に示した演算機能１３２と基本的に同様の処理を行う。また、画像生成機能３２２は、図１に示した画像生成回路１５０と基本的に同様の処理を行う。演算機能３２１及び画像生成機能３２２は、例えば、コンピュータによって実行可能なプログラムの形態で記憶回路３１０内に記録されている。処理回路３２０は、各プログラムを読み出し、実行することで読み出した各プログラムに対応する機能（演算機能３２１及び画像生成機能３２２）を実現する。なお、図示しないが、情報処理装置３００では、図１に示したＢモード処理回路１３０及びドプラ処理回路１４０と同様の機能が実現される。

すなわち、高機能型超音波プローブ２００は、第１の位相の第１超音波の送信及び前記第１の位相と実質的に９０度異なる第２の位相の第２超音波の送信のみからなる１セットの超音波送信を、複数の走査線に沿って順々に実行する。そして、演算機能３２１は、高機能型超音波プローブ２００を介して取得された、第１超音波に対応する第１エコー信号から第２超音波に対応する第２エコー信号を減算することにより減算信号を生成する。そして、画像生成機能３２２は、複数の走査線について生成された減算信号に基づいて、超音波画像を生成する。

これによれば、超音波診断装置２は、上述した超音波診断装置１と同様にノイズを除去することができる。つまり、上述した実施形態は、２Ｄアレイプローブに限らず、例えば、高機能型超音波プローブ２００を用いた超音波診断装置２に適用される場合にも、周期ノイズや固定ノイズを除去することができる。

なお、図８は一例に過ぎず、図８に示した内容に限定されるものではない。例えば、図８では、高機能型超音波プローブ２００で収集された反射波データが、情報処理装置３００側で画像化される場合を例示したが、これに限定されるものではない。例えば、反射波データが高機能型超音波プローブ２００側で画像化され、画像データが高機能型超音波プローブ２００から情報処理装置３００へ送信される構成であってもよい。例えば、図９に示すように、高機能型超音波プローブ２００が、処理回路２２０を有する。処理回路２２０は、演算機能２２１及び画像生成機能２２２を有する。そして、高機能型超音波プローブ２００において、演算機能２２１は、第１超音波に対応する第１エコー信号から第２超音波に対応する第２エコー信号を減算することにより減算信号を生成する。また、高機能型超音波プローブ２００において、画像生成機能２２２は、複数の走査線について生成された減算信号に基づいて、超音波画像を生成する。そして、高機能型超音波プローブ２００は、ノイズが除去された超音波画像を情報処理装置３００へ送信する。図９に示す超音波診断装置２は、図８に示した構成と比較して、高機能型超音波プローブ２００から情報処理装置３００へ送信するデータ量を低減させることができる。このため、図９に示す超音波診断装置２は、高機能型超音波プローブ２００から情報処理装置３００への送信ができないというリスクを低減させることができる。

上記説明において用いた「プロセッサ」という文言は、例えば、ＣＰＵ（Central Processing Unit）、ＧＰＵ(Graphics Processing Unit)、或いは、特定用途向け集積回路（Application Specific Integrated Circuit：ＡＳＩＣ））、プログラマブル論理デバイス（例えば、単純プログラマブル論理デバイス（Simple Programmable Logic Device：ＳＰＬＤ）、複合プログラマブル論理デバイス（Complex Programmable Logic Device：ＣＰＬＤ）、及びフィールドプログラマブルゲートアレイ（Field Programmable Gate Array：ＦＰＧＡ））等の回路を意味する。プロセッサは記憶回路に保存されたプログラムを読み出し実行することで機能を実現する。なお、記憶回路１７０にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むよう構成しても構わない。この場合、プロセッサは回路内に組み込まれたプログラムを読み出し実行することで機能を実現する。なお、本実施形態の各プロセッサは、プロセッサごとに単一の回路として構成される場合に限らず、複数の独立した回路を組み合わせて１つのプロセッサとして構成し、その機能を実現するようにしてもよい。さらに、図１における複数の構成要素を１つのプロセッサへ統合してその機能を実現するようにしてもよい。

また、上記の実施形態において説明した各処理のうち、自動的に行われるものとして説明した処理の全部又は一部を手動的に行うこともでき、或いは、手動的に行われるものとして説明した処理の全部又は一部を公知の方法で自動的に行うこともできる。この他、上記文書中や図面中で示した処理手順、制御手順、具体的名称、各種のデータやパラメータを含む情報については、特記する場合を除いて任意に変更することができる。

また、上記の実施形態で説明した超音波イメージング方法は、予め用意された超音波イメージングプログラムをパーソナルコンピュータやワークステーション等のコンピュータで実行することによって実現することができる。この超音波イメージング方法は、インターネット等のネットワークを介して配布することができる。また、この超音波イメージング方法は、ハードディスク、フレキシブルディスク（ＦＤ）、ＣＤ－ＲＯＭ、ＭＯ、ＤＶＤ等のコンピュータで読み取り可能な記録媒体に記録され、コンピュータによって記録媒体から読み出されることによって実行することもできる。

以上説明した少なくともひとつの実施形態によれば、ノイズを除去することができる。

本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。

１ 超音波診断装置
１００ 装置本体
１０１ 超音波プローブ
１３０ Ｂモード処理回路
１３２ 演算機能
１５０ 画像生成回路

Claims (3)

1. 第１の位相の第１超音波の送信及び前記第１の位相と９０度又は２７０度異なる第２の位相の第２超音波の送信のみからなる１セットの超音波送信を、複数の走査線に沿って順々に実行する超音波プローブであって、マトリクス状に配列された複数の振動素子を有し、複数の振動素子で構成されるサブアレイと、複数のサブアレイで構成されるメインアレイと、前記サブアレイ毎にエコー信号の遅延加算処理を実行する遅延加算回路とを有する超音波プローブと、
   前記超音波プローブに含まれる振動素子の接続先を、送信回路及び受信回路を含む選択肢の中から選択的に切り替えるスイッチング回路と、
   前記遅延加算回路の処理において、遅延メモリに記憶された各チャンネルの遅延時間を繰り返し更新させるためのアナログスイッチに由来して周期的に発生する周期ノイズ、および、前記スイッチング回路における送受信切り替え時のスイッチノイズが前記振動素子に伝わることで本来送信されるべきでないタイミングで送信される微弱な超音波に由来する固定ノイズを除去しつつ高調波成分を増幅させるために、前記超音波プローブを介して取得された、前記第１超音波に対応する第１エコー信号から前記第２超音波に対応する第２エコー信号を減算することにより減算信号を生成する演算部と、
   複数の走査線について生成された前記減算信号に対してハイパスフィルタ又はバンドパスフィルタを適用し、前記ハイパスフィルタ又はバンドパスフィルタを適用後の信号に基づいて、超音波画像を生成する画像生成部と、
   を備えた超音波診断装置。
2. 前記画像生成部は、前記第１エコー信号及び前記第２エコー信号、又は前記減算信号に含まれる基本波成分を取り除き、少なくとも２次高調波成分を抽出するフィルタを適用し、
   前記画像生成部は、抽出された前記２次高調波成分を用いて前記超音波画像を生成する、
   請求項１に記載の超音波診断装置。
3. 第１の位相の第１超音波の送信及び前記第１の位相と９０度又は２７０度異なる第２の位相の第２超音波の送信のみからなる１セットの超音波送信を、複数の走査線に沿って順々に実行する超音波走査を超音波プローブであって、マトリクス状に配列された複数の振動素子を有し、複数の振動素子で構成されるサブアレイと、複数のサブアレイで構成されるメインアレイと、前記サブアレイ毎にエコー信号の遅延加算処理を実行する遅延加算回路とを有する超音波プローブに実行させ、
   前記遅延加算回路の処理において、遅延メモリに記憶された各チャンネルの遅延時間を繰り返し更新させるためのアナログスイッチに由来して周期的に発生する周期ノイズ、および、前記超音波プローブに含まれる振動素子の接続先を送信回路及び受信回路を含む選択肢の中から選択的に切り替えるスイッチング回路における送受信切り替え時のスイッチノイズが前記振動素子に伝わることで本来送信されるべきでないタイミングで送信される微弱な超音波に由来する固定ノイズを除去しつつ高調波成分を増幅させるために、前記超音波プローブを介して取得された、前記第１超音波に対応する第１エコー信号から前記第２超音波に対応する第２エコー信号を減算することにより減算信号を生成し、
   複数の走査線について生成された前記減算信号に対してハイパスフィルタ又はバンドパスフィルタを適用し、前記ハイパスフィルタ又はバンドパスフィルタを適用後の信号に基づいて、超音波画像を生成する、
   各処理をコンピュータに実行させる、超音波イメージングプログラム。

Images