(57)【要約】【課題】バブル組織比の高い造影画像を生成すること。【解決手段】実施形態の超音波診断装置は、フィルタと、フィルタ処理部と、加減算部とを備える。フィルタは、造影剤非存在下において、振幅変調法又は振幅位相変調法のいずれかの変調法に基づいて同一走査線で複数回送信された超音波それぞれの複数の受信信号である複数の非造影受信信号を、前記変調法に応じた加減算により合成した非造影合成信号のエネルギーを最小化するように設計されたフィルタ係数が設定される。フィルタ処理部は、前記フィルタを有し、造影撮影時に、前記変調法に応じて同一走査線で複数回送信された超音波それぞれの複数の受信信号の少なくとも１つの受信信号に対して、前記フィルタをかける。加減算部は、前記フィルタ処理後の複数の受信信号を、前記変調法に応じた加減算により合成した合成信号を出力する。【選択図】図２
2015016251.tif 000002
【特許請求の範囲】
【請求項１】
造影剤非存在下において、振幅変調法又は振幅位相変調法のいずれかの変調法に基づいて同一走査線で複数回送信された超音波それぞれの複数の受信信号である複数の非造影受信信号を、前記変調法に応じた加減算により合成した非造影合成信号のエネルギーを最小化するように設計されたフィルタ係数が設定されたフィルタと、
前記フィルタを有し、造影撮影時に、前記変調法に応じて同一走査線で複数回送信された超音波それぞれの複数の受信信号の少なくとも１つの受信信号に対して、前記フィルタをかけるフィルタ処理部と、
前記フィルタ処理後の複数の受信信号を、前記変調法に応じた加減算により合成した合成信号を出力する加減算部と、
を備える超音波診断装置。
【請求項２】
前記複数の非造影受信信号それぞれは、信号レベルが飽和していない不飽和領域からの受信信号であることを特徴とする請求項１に記載の超音波診断装置。
【請求項３】
前記フィルタのフィルタ長は、送信超音波のパルス長の略２倍であることを特徴とする請求項１又は２に記載の超音波診断装置。
【請求項４】
前記複数の非造影受信信号それぞれは、非線形伝播の発生が低減される音圧の超音波送信により受信された受信信号であることを特徴とする請求項１〜３のいずれか１つに記載の超音波診断装置。
【請求項５】
前記複数の非造影受信信号それぞれは、ファントム、又は、生体から受信した受信信号であることを特徴とする請求項１〜４のいずれか１つに記載の超音波診断装置。
【請求項６】
複数の送信条件ごとに予め設計された複数のフィルタ係数を記憶するフィルタ係数テーブル、
を更に備え、
前記フィルタ処理部は、造影撮影時の送信条件に対応するフィルタ係数を前記フィルタ係数テーブルから取得して、取得したフィルタ係数を前記フィルタに設定することを特徴とする請求項１〜５のいずれか１つに記載の超音波診断装置。
【請求項７】
前記複数の非造影受信信号に基づいて、前記フィルタ係数を設計するフィルタ設計部、
を更に備え、
前記フィルタ処理部は、前記フィルタ設計部が設計した前記フィルタ係数を前記フィルタに設定することを特徴とする請求項１〜５のいずれか１つに記載の超音波診断装置。
【請求項８】
前記フィルタ設計部が用いる前記複数の非造影受信信号それぞれは、造影撮影と同一の送信条件により受信される受信信号であって、造影撮影が行なわれる被検体の撮影部位であり、造影剤が存在しない状態の撮影部位から受信した受信信号であることを特徴とする請求項７に記載の超音波診断装置。
【請求項９】
前記フィルタ設計部が用いる前記複数の非造影受信信号それぞれは、造影剤が到達していない前記撮影部位又は造影剤投与前の前記撮影部位から受信した受信信号であることを特徴とする請求項８に記載の超音波診断装置。
【請求項１０】
前記フィルタ設計部が用いる前記複数の非造影受信信号それぞれは、造影剤投与後に造影剤を破壊可能な音圧の超音波送信が行なわれた前記撮影部位から受信した受信信号であることを特徴とする請求項８又は９に記載の超音波診断装置。
【請求項１１】
前記不飽和領域は、操作者により設定されることを特徴とする請求項２〜１０のいずれか１つに記載の超音波診断装置。
【請求項１２】
前記フィルタ設計部は、前記複数の非造影受信信号の信号レベルに基づいて、当該複数の非造影受信信号が前記不飽和領域からの受信信号であるか否かを判定することを特徴とする請求項７〜１０のいずれか１つに記載の超音波診断装置。
【請求項１３】
前記フィルタ設計部は、操作者の入力を契機として、前記フィルタ係数を設計することを特徴とする請求項７〜１２のいずれか１つに記載の超音波診断装置。
【請求項１４】
前記フィルタ設計部は、造影剤投与タイマーに連動して、前記フィルタ係数を設定することを特徴とする請求項７〜１２のいずれか１つに記載の超音波診断装置。
【請求項１５】
前記フィルタ処理部は、前記フィルタをＩＱ信号又はＲＦ信号に対してかけることを特徴とする請求項１〜１４のいずれか１つに記載の超音波診断装置。
【請求項１６】
前記フィルタは、複素有限インパルス応答フィルタ、又は、実数有限インパルス応答フィルタであることを特徴とする請求項１５のいずれか１つに記載の超音波診断装置。
【請求項１７】
造影撮影時に、振幅変調法又は振幅位相変調法のいずれかの変調法に基づいて同一走査線で複数回送信された超音波それぞれの複数の受信信号の少なくとも１つの受信信号に対して、フィルタをかけるフィルタ処理手順と、
前記フィルタ処理後の複数の受信信号を、前記変調法に応じた加減算により合成した合成信号を出力する加減算手順と、
をコンピュータに実行させ、
前記フィルタは、造影剤非存在下において、前記変調法に応じて同一走査線で複数回送信された超音波それぞれの複数の受信信号である複数の非造影受信信号を、前記変調法に応じた加減算により合成した非造影合成信号のエネルギーを最小化するように設計されたフィルタ係数が設定される
ことを特徴とする超音波診断装置の制御プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、超音波診断装置及びその制御プログラムに関する。
【背景技術】
【０００２】
近年、静脈投与型の超音波造影剤が製品化され、超音波診断装置により、ＣＨＩ（Contrast  Harmonic  Imaging）と呼ばれる造影エコー法が行われている。造影エコー法は、例えば、心臓や肝臓等の検査において、静脈から超音波造影剤を注入して血流信号を増強し、血流動態の評価を行うことを目的としている。超音波造影剤の多くは、微小気泡（マイクロバブル）が反射源として機能するものである。しかし、気泡というデリケートな基材の性質上、通常の診断レベルの超音波照射であっても、超音波の機械的作用によって気泡が崩壊し、結果的にスキャン面からの信号強度が低下してしまう。
【０００３】
従って、還流の動的な様子をリアルタイムで観察するためには、低音圧の超音波送信によって画像化する等、スキャンによる気泡の崩壊を比較的低減させることが必要となる。この様な低音圧の超音波送信による画像化では、信号／雑音比（Ｓ／Ｎ比）も低下してしまう。これを補うため、位相変調法(ＰＭ：Phase  Modulation)や振幅変調法(ＡＭ：Amplitude  Modulation)、位相振幅変調法(ＡＭＰＭ)のように、様々な信号処理法が考案されている。これらの映像化手法により、リアルタイムで高いＳ／Ｎ比の造影画像を表示することが可能となっている。超音波造影は、リアルタイム性や高空間分解能から、Ｘ線ＣＴ装置やＭＲＩ装置では視覚化できない微小構造（例えば、微小血管構造）の精査に利用される。
【０００４】
例えば、振幅変調法は、バブル組織比や深部感度に優れた映像化手法である。振幅変調法は、造影剤の非線形応答を抽出する一方、組織からの線形信号をキャンセルし、造影剤を特異的に抽出する映像化手法である。振幅変調法の実現には、高精度の波形形成が求められる。しかし、超音波診断装置は、システム構成や開口制御・振幅制御などの実装方式、回路の非線形性により、組織由来の線形信号を完全にはキャンセルできない場合がある。かかる場合、組織由来の線形信号が残存し、造影画像におけるバブル組織比が低下する。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平１０−１４９２１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明が解決しようとする課題は、バブル組織比の高い造影画像を生成することできる超音波診断装置及びその制御プログラムを提供することである。
【課題を解決するための手段】
【０００７】
実施形態の超音波診断装置は、フィルタと、フィルタ処理部と、加減算部とを備える。フィルタは、造影剤非存在下において、振幅変調法又は振幅位相変調法のいずれかの変調法に基づいて同一走査線で複数回送信された超音波それぞれの複数の受信信号である複数の非造影受信信号を、前記変調法に応じた加減算により合成した非造影合成信号のエネルギーを最小化するように設計されたフィルタ係数が設定される。フィルタ処理部は、前記フィルタを有し、造影撮影時に、前記変調法に応じて同一走査線で複数回送信された超音波それぞれの複数の受信信号の少なくとも１つの受信信号に対して、前記フィルタをかける。加減算部は、前記フィルタ処理後の複数の受信信号を、前記変調法に応じた加減算により合成した合成信号を出力する。
【図面の簡単な説明】
【０００８】
【図１】図１は、第１の実施形態に係る超音波診断装置の構成例を示すブロック図である。
【図２】図２は、第１の実施形態に係るＢモード処理ユニットの構成例を示すブロック図である。
【図３】図３は、第１の実施形態に係るフィルタ係数を説明するための図である。
【図４】図４は、第１の実施形態に係るフィルタ係数設計部によるフィルタ係数設計処理を説明するための図である。
【図５】図５は、第１の実施形態に係るフィルタ係数設計部によるフィルタ係数設計処理を説明するための図である。
【図６】図６は、第１の実施形態に係る超音波診断装置による処理の手順を示すフローチャートである。
【図７】図７は、第１の実施形態に係る超音波診断装置による超音波の送受信処理の手順を示すフローチャートである。
【図８】図８は、第１の実施形態に係る超音波診断装置によるＢモードデータ生成処理の手順を示すフローチャートである。
【図９】図９は、第２の実施形態に係るＢモード処理ユニットの構成例を示すブロック図である。
【図１０】図１０は、第２の実施形態に係る超音波診断装置による処理の手順を示すフローチャートである。
【図１１】図１１は、第２の実施形態に係る超音波診断装置によるフィルタ設計処理の手順を示すフローチャートである。
【図１２】図１２は、第３の実施形態に係る超音波診断装置による処理の手順を示すフローチャートである。
【発明を実施するための形態】
【０００９】
（第１の実施形態）
まず、図１を用いて、第１の実施形態に係る超音波診断装置の構成について説明する。図１は、第１の実施形態に係る超音波診断装置１の構成例を示すブロック図である。かかる超音波診断装置１は、図１に例示するように、超音波プローブ１０と、入力装置２０と、モニタ３０と、装置本体１００とを有する。
【００１０】
超音波プローブ１０は、複数の圧電振動子を有し、これら複数の圧電振動子は、後述する装置本体１００が有する送受信ユニット１１０から供給される駆動信号に基づき超音波を発生する。また、超音波プローブ１０は、被検体Ｐからの反射波信号を受信して電気信号に変換する。また、超音波プローブ１０は、圧電振動子に設けられる整合層と、圧電振動子から後方への超音波の伝播を防止するバッキング材などを有する。なお、超音波プローブ１０は、装置本体１００と着脱自在に接続される。
【００１１】
超音波プローブ１０から被検体Ｐに超音波が送信されると、送信された超音波は、被検体Ｐの体内組織における音響インピーダンスの不連続面で次々と反射され、反射波信号として超音波プローブ１０が有する複数の圧電振動子にて受信される。受信される反射波信号の振幅は、超音波が反射される不連続面における音響インピーダンスの差に依存する。なお、送信された超音波パルスが、移動している血流や心臓壁などの表面で反射された場合の反射波信号は、ドプラ効果により、移動体の超音波送信方向に対する速度成分に依存して、周波数偏移を受ける。
【００１２】
入力装置２０は、装置本体１００と接続され、トラックボール２１、各種スイッチ２２、各種ボタン２３、マウス２４、キーボード２５などを有する。かかる入力装置２０は、操作者からの各種指示を装置本体１００に通知する。例えば、各種指示には、関心領域（ＲＯＩ：Region  of  Interest）の設定指示、超音波の送信条件を含む撮影条件の設定指示、被検体Ｐに造影剤を投与してからの経過時間を表示する指示等が含まれる。
【００１３】
モニタ３０は、超音波診断装置１の操作者が入力装置２０を用いて各種設定を行うためのＧＵＩ（Graphical  User  Interface）を表示したり、装置本体１００において生成された超音波画像などを表示したりする。具体的には、モニタ３０は、後述する画像生成部１４０から入力されるビデオ信号に基づいて、生体内の形態学的情報や血流情報を画像として表示する。また、モニタ３０は、被検体Ｐに造影剤を投与してからの経過時間を表示する指示を受付けた場合、造影剤を投与後の時間を表示する。
【００１４】
装置本体１００は、超音波プローブ１０が受信した反射波信号に基づいて超音波画像を生成する。かかる装置本体１００は、図１に例示するように、送受信ユニット１１０と、Ｂモード処理ユニット１２０と、ドプラ処理ユニット１３０と、画像生成部１４０と、画像メモリ１５０と、ソフトウェア格納部１６０と、インタフェース部１７０と、記憶部１８０と、制御部１９０とを有する。なお、装置本体１００に内蔵される送受信ユニット１１０、Ｂモード処理ユニット１２０、ドプラ処理ユニット１３０、画像生成部１４０等は、集積回路などのハードウェアにより実現されてもよいし、ソフトウェア的にモジュール化されたソフトウェアプログラムにより実現されてもよい。
【００１５】
送受信ユニット１１０は、図示しない遅延回路、パルサ回路、トリガ発生回路等を有し、超音波プローブ１０に駆動信号を供給する。パルス発生回路は、所定の繰り返し周波数（ＰＲＦ：Pulse  Repetition  Frequency）で、送信超音波を形成するためのレートパルスを繰り返し発生する。なお、繰り返し周波数は、レート周波数等とも呼ばれる。また、遅延回路は、超音波プローブ１０から発生される超音波をビーム状に集束し、かつ送信指向性を決定するために必要な圧電振動子ごとの遅延時間を、各レートパルスに対し与える。また、トリガ発生回路は、遅延回路により遅延時間が与えられた各レートパルスに基づくタイミングで、超音波プローブ１０に駆動信号（駆動パルス）を印加する。なお、送信方向に応じた遅延時間は記憶部１８０に記憶されており、遅延回路は、記憶部１８０を参照して遅延時間を与える。
【００１６】
また、送受信ユニット１１０は、図示しないアンプ回路、Ａ／Ｄ（Analog／Digital）変換器、加算器等を有し、超音波プローブ１０が受信した反射波信号に対して各種処理を行って、例えば、ＲＦ（Radio  Frequency）信号を反射波データとして生成する。アンプ回路は、反射波信号をチャネル毎に増幅する。Ａ／Ｄ変換器は、増幅された反射波信号をＡ／Ｄ変換し、受信指向性を決定するために必要な遅延時間を与える。加算器は、遅延時間が与えられた反射波信号の加算処理を行なって反射波データを生成する。加算器の加算処理により、反射波信号の受信指向性に応じた方向からの反射成分が強調され、受信指向性と送信指向性とにより超音波送受信の総合的なビームが形成される。なお、受信方向に応じた遅延時間は記憶部１８０に記憶されている。また、以下では、反射波データのことを「受信信号」と呼ぶ場合がある。
【００１７】
なお、送受信ユニット１１０は、制御部１９０からの指示に従って、遅延情報、送信周波数、送信駆動電圧、開口素子数等を瞬時に変更可能な機能を有している。特に、送信駆動電圧の変更は、瞬間にその値を切り替え可能なリニアアンプ型の発信回路、または、複数の電源ユニットを電気的に切り替える機構によって実現される。このように、送受信ユニット１１０は、超音波の送受信における送信指向性と受信指向性とを制御する。
【００１８】
Ｂモード処理ユニット１２０は、送受信ユニット１１０から反射波データを受け取り、対数増幅、包絡線検波処理などを行って、信号強度が輝度の明るさで表現されるデータ（Ｂモードデータ）を生成する。Ｂモードデータは、走査線上の各サンプル点に、信号強度に応じた輝度が割り当てられたデータである。ここで、Ｂモード処理ユニット１２０は、高調波成分を映像化するハーモニックイメージングを行なうための信号処理を行なうことができる。
【００１９】
例えば、ハーモニックイメージングとしては、コントラストハーモニックイメージング（ＣＨＩ：Contrast  Harmonic  Imaging）やティッシュハーモニックイメージング（ＴＨＩ：Tissue  Harmonic  Imaging）が知られている。また、ハーモニックイメージングには、スキャン方式として、振幅変調（ＡＭ：Amplitude  Modulation）法、位相変調（ＰＭ：Phase  Modulation）法、ＡＭ法とＰＭ法とを組み合わせることで、ＡＭ法の効果及びＰＭ法の効果の双方が得られるＡＭＰＭ法が知られている。
【００２０】
なお、上記の送受信ユニット１１０は、造影像を生成するための映像化手法であるＣＨＩを行なう場合、超音波の走査線毎に、異なる波形を複数回送信する。例えば、ＣＨＩをＡＭ法で行なう場合、送受信ユニット１１０は、１回目の送信波形に対して同じ位相極性で振幅の比率を１/２に変調した波形を２回目に送信し、反射波データをそれぞれ生成する。なお、ハーモニックイメージングを行なう場合、送受信ユニット１１０は、後述する制御部１９０によって設定されたスキャンシーケンスで超音波を送信する。Ｂモード処理ユニット１２０は、上記のＡＭ法の場合、マイクロバブルが注入された被検体Ｐに対する２つの反射波データを、送受信ユニット１１０から受信する。そして、Ｂモード処理ユニット１２０は、送受信ユニット１１０から受信した２つの反射波データの振幅の比率を補正した後、２つの反射波データを差分することにより、基本波成分が抑制され、高調波成分（非線形成分）が抽出された反射波データを生成する。なお、ＡＭ法では、例えば、同じ位相極性で、振幅の比率を「１：２：１」に変調させた超音波を各走査線で３回送信させ、３つの反射波データを受信してもよい。この場合の、ＡＭ法では、振幅の比率が「１」である反射波データを加算し、加算した反射波データと振幅の比率が「２」である反射波データと差分することにより、高調波成分（非線形成分）が抽出された反射波データを生成する。
【００２１】
また、例えば、ＣＨＩをＡＭＰＭ法で行なう場合、送受信ユニット１１０は、１回目の送信波形に対して位相極性を反転させ、かつ振幅の比率を１/２に変調した波形を２回目に送信し、反射波データをそれぞれ生成する。Ｂモード処理ユニット１２０は、上記のＡＭＰＭ法の場合、マイクロバブルが注入された被検体Ｐに対する２つの反射波データを、送受信ユニット１１０から受信する。そして、Ｂモード処理ユニット１２０は、送受信ユニット１１０から受信した２つの反射波データの振幅の比率を補正した後、２つの反射波データを加算することにより、基本波成分が抑制され、高調波成分（非線形成分）が抽出された反射波データを生成する。なお、ＡＭＰＭ法では、例えば、振幅の比率を「１：２：１」に変調させ、更に、１回目及び３回目の送信超音波の極性と２回目の送信超音波の極性とを反転させた超音波を各走査線で３回送信させ、受信した３つの反射波データを加算することにより、高調波成分（非線形成分）が抽出された反射波データを生成してもよい。
【００２２】
続いて、Ｂモード処理ユニット１２０は、高調波成分が抽出された反射波データに対して、包絡線検波処理等を行なって、造影画像を生成するためのＢモードデータを生成する。この結果、後述する画像生成部１４０は、被検体Ｐ内を流動する造影剤を高感度に映像化した造影画像、及び組織を映像化した組織画像を生成することができる。
【００２３】
ドプラ処理ユニット１３０は、送受信ユニット１１０から受け取った反射波データから速度情報を周波数解析し、ドプラ効果による血流や組織、造影剤エコー成分を抽出し、平均速度、分散、パワー等の血流情報（ドプラデータ）を多点について算出する。
【００２４】
画像生成部１４０は、Ｂモード処理ユニット１２０が生成したＢモードデータから、信号強度が輝度の明るさで表現されるＢモード画像を生成し、ドプラ処理ユニット１３０が生成した血流情報から、血流の速度、分散、血流量等を示すパワー成分等を色によって識別可能に表示するカラードプラ画像を生成する。なお、画像生成部１４０に入力される前のデータは「生データ」と呼ばれることがある。
【００２５】
具体的には、画像生成部１４０は、Ｂモードデータ及びドプラデータに対してフィルタリング処理を実行することにより、超音波スキャン走査線信号列からノイズ成分を除去し、フィルタリング処理後のデータを画像メモリ１５０に格納する。そして、画像生成部１４０は、フィルタリング処理を実行したデータの超音波スキャン走査線信号列をテレビなどの一般的なビデオフォーマットの走査線信号列に変換する。画像生成部１４０は、走査線信号列に対して、輝度やコントラストの調整処理や、空間フィルタ等の画像処理、もしくは種々の設定パラメータの文字情報やメモリ等を合成する合成処理を実行し、ビデオ信号としてモニタ３０に出力する。これにより、画像生成部１４０によって生成された被検体組織形状を表す断層像等の超音波画像は、モニタ３０に表示される。
【００２６】
画像メモリ１５０は、画像生成部１４０が生成した超音波画像や、超音波画像を画像処理することで生成した画像を記憶するメモリである。例えば診断の後に、操作者が検査中に記録された画像を画像メモリ１５０から呼び出すことが可能となっており、静止画像的に、あるいは複数枚を使って動画的に再生することが可能である。また、画像メモリ１５０は、送受信ユニット１１０を通過した後の画像輝度信号、その他の生データ、ネットワークを介して取得した画像等を必要に応じて記憶する。
【００２７】
ソフトウェア格納部１６０は、後述する制御部１９０によって各種装置制御プログラムが展開される記憶領域である。
【００２８】
インタフェース部１７０は、入力装置２０、外部装置（図示を省略）、ネットワークに関するインタフェースである。超音波診断装置１によって得られた超音波画像などのデータは、インタフェース部１７０によって、ネットワークを介して他の装置へ転送することができる。
【００２９】
記憶部１８０は、スキャンシーケンス、画像処理及び表示処理等を実行するための各種装置制御プログラムや、診断情報（例えば、患者ＩＤ、医師の所見等）、診断プロトコルや各種設定情報等の各種データ群を記憶する。なお、各種装置制御プログラムは、制御部１９０と同様の処理を実行する手順が記述されたプログラムを含む場合もある。また、記憶部１８０は、必要に応じて、画像メモリ１５０が記憶する超音波画像の保管などにも使用される。なお、記憶部１８０が記憶する各種データは、インタフェース部１７０を経由して、外部装置へ転送することができる。
【００３０】
制御部１９０は、情報処理装置（計算機）としての機能を実現する制御プロセッサ（ＣＰＵ：Central  Processing  Unit）であり、超音波診断装置１における処理全体を制御する。具体的には、制御部１９０は、入力装置２０を介して操作者から入力された各種指示や設定指示、記憶部１８０から読み込んだ各種装置制御プログラムをソフトウェア格納部１６０に展開し、各種設定情報に基づいて、送受信ユニット１１０、Ｂモード処理ユニット１２０、ドプラ処理ユニット１３０及び画像生成部１４０の処理を制御したり、画像メモリ１５０が記憶する超音波画像などをモニタ３０にて表示するように制御したりする。
【００３１】
以上、第１の実施形態に係る超音波診断装置１の全体構成について説明した。かかる構成のもと、第１の実施形態に係る超音波診断装置１は、例えば、ＡＭ法でＣＨＩを行なうことで、マイクロバブルが注入された被検体Ｐに対して、２次高調波成分がより強調された造影像を生成する。しかし、超音波診断装置１は、システム構成や開口制御・振幅制御などの実装方式、回路の非線形性により、組織由来の線形信号を完全にはキャンセルできない場合がある。かかる場合、組織由来の線形信号が残存し、造影画像におけるバブル組織比が低下する。
【００３２】
このようなことから、第１の実施形態に係る超音波診断装置１において、Ｂモード処理ユニット１２０は、送受信ユニット１１０から受信した反射波データを、変調法に応じたビーム加減算処理を行う前に波形整形を行い、組織由来の線形信号をキャンセルする。以下に、図２〜図６を用いて、第１の実施形態に係るＢモード処理ユニット１２０についてより詳細に説明する。
【００３３】
図２は、第１の実施形態に係るＢモード処理ユニット１２０の構成例を示すブロック図である。図２に例示するように、Ｂモード処理ユニット１２０は、直交検波部１２１ａと、フィルタ処理部１２１ｂと、フィルタ係数テーブル１２１ｃと、加減算部１２１ｄとを有する。なお、説明の便宜上、以下では、超音波診断装置１が、ＡＭ法でＣＨＩを行う場合に、例えば、（０．５，１．０）のように、同じ位相極性で振幅の比率を「１：２」に変調させた超音波を各走査線で２回送信させる例を仮定する。なお、振幅の比率が「２」である超音波を大振幅送信レートと呼び、振幅の比率が「１」である超音波を小振幅送信レートと呼ぶ。
【００３４】
直交検波部１２１ａは、反射波データとして送受信ユニット１１０から出力されたＲＦ信号を、ベースバンド帯域の同相信号（Ｉ信号、Ｉ：In-pahse）と直交信号（Ｑ信号、Ｑ：Quadrature-phase）とに変換する直交検波を行なう。そして、直交検波部１２１ａは、Ｉ信号及びＱ信号（以下、ＩＱ信号と記載する）を、反射波データ（受信信号）として後段の処理部に出力する。本実施形態では、直交検波部１２１ａは、小振幅送信レートのＩＱ信号を加減算部１２１ｄに出力し、大振幅送信レートのＩＱ信号をフィルタ処理部１２１ｂに出力する。なお、このＩＱ信号は、反射波データ（受信信号）の一種である。
【００３５】
フィルタ処理部１２１ｂは、図示しないフィルタを有する。この実施形態に係るフィルタは、ＦＩＲ（Finite  Impulse  Response）フィルタである。また、この実施形態では、ＩＱ信号に対してフィルタ処理を行うので、フィルタは、複素ＦＩＲフィルタである。このフィルタには、フィルタ係数が設定される。
【００３６】
フィルタ係数は、造影剤非存在下において、ＡＭ法又はＡＭＰＭ法のいずれかの変調法に基づいて同一走査線で複数回送信された超音波それぞれの複数の受信信号である複数の非造影受信信号を、変調法に応じた加減算により合成した非造影合成信号のエネルギーを最小化するように設計される。非造影受信信号は、造影剤非存在下において、受信される受信信号である。
【００３７】
ここで、フィルタ係数の設計上、複数の非造影受信信号それぞれは、信号レベルが飽和していない不飽和領域からの受信信号であることが望ましい。また、フィルタ係数の設計上、複数の非造影受信信号それぞれは、非線形伝播の発生が低減される音圧の超音波送信により受信された受信信号であることが望ましい。
【００３８】
フィルタ処理部１２１ｂは、造影撮影時に、変調法に応じて同一走査線で複数回送信された超音波それぞれの複数の受信信号の少なくとも１つの受信信号に対して、フィルタをかける。例えば、フィルタ処理部１２１ｂは、大振幅送信レートのＩＱ信号に対して、畳み込みを行う。言い換えると、フィルタ処理部１２１ｂは、フィルタ係数テーブル１２１ｃが記憶するフィルタ係数を用いて、大振幅送信レートのＩＱ信号を補正する。なお、フィルタ処理部１２１ｂは、ＡＭＰＭ法においても同様に、大振幅送信レートのＩＱ信号に対して、畳み込みを行う。そして、フィルタ処理部１２１ｂは、補正したＩＱ信号を加減算部１２１ｄに出力する。なお、フィルタ処理部１２１ｂは、造影撮影時の送信条件に対応するフィルタ係数をフィルタ係数テーブル１２１ｃから取得して、取得したフィルタ係数をフィルタに設定する。
【００３９】
フィルタ係数テーブル１２１ｃは、複数の送信条件（周波数、送信フォーカスの深さ方向の位置、及び音圧等）ごとに予め設計された複数のフィルタ係数を記憶する。例えば、フィルタ係数テーブル１２１ｃは、造影剤が投与されていない被検体Ｐから受信した各反射波データの加減算後のエネルギーを最小にするフィルタ係数を記憶する。言い換えると、フィルタ係数テーブル１２１ｃは、組織由来の線形信号を最小にするフィルタ係数を記憶する。このフィルタ係数テーブル１２１ｃは、例えば、他の超音波診断装置から取得したフィルタ係数を記憶する。なお、フィルタ係数の設計方法については、後述する。
【００４０】
図３は、第１の実施形態に係るフィルタ係数を説明するための図である。図３は、造影剤非存在下でのＡＭ法により得られた非造影合成信号に残存する線形成分を最小化させるように設計されたフィルタ係数の一例である。また、図３は、非造影合成信号のエネルギーを最小化するために大振幅レートの受信信号（ＩＱ信号）に対してフィルタ処理を行なうために設計されたフィルタ係数の一例である。図３の（Ａ）は、時間成分で表したフィルタ係数の一例を示す。図３の（Ａ）において、縦軸は、フィルタ係数を示し、横軸は、時間を示す。ここで、時間は、超音波を送信後の時間に対応した走査線上の深さ方向を示す。また、図３の（Ａ）において、ＩＱ信号の実部（Ｉ信号）に対するフィルタ係数を３ａで示し、ＩＱ信号の虚部（Ｑ信号）に対するフィルタ係数を３ｂで示す。図３の（Ａ）に示すように、各深さには、フィルタ係数が各々対応付けられる。フィルタ処理部１２１ｂは、時間領域でフィルタ処理を行う場合、ＩＱ信号の実部の波形を深さ毎に順次サンプリングし、サンプリングしたＩＱ信号の実部を深さ毎に予め設定されたフィルタ係数で成形する。また、フィルタ処理部１２１ｂは、時間領域でフィルタ処理を行う場合、ＩＱ信号の虚部の波形を深さ方向に順次サンプリングし、サンプリングしたＩＱ信号の虚部を深さ毎に予め設定されたフィルタ係数で成形する。
【００４１】
図３の（Ｂ）は、図３の（Ａ）に示す時間成分で表したフィルタ係数をフーリエ変換によって周波数成分で表した図である。図３（Ｂ）の横軸は、周波数を示し、縦軸は、振幅特性或いは位相特性を示す。また、図３の（Ｂ）において、振幅特性を３ｃで示し、位相特性を３ｄで示す。図３の（Ｂ）に示すように、各周波数には、振幅特性及び位相特性が各々対応付けられる。フィルタ処理部１２１ｂは、周波数領域でフィルタ処理を行う場合、ＩＱ信号の波形を周波数毎に順次サンプリングし、周波数毎に予め設定された振幅特性及び位相特性で成形する。図３の（Ｂ）に例示するフィルタ係数の振幅特性は、画像化される周波数帯域において凹状の形状となり、フィルタ処理では、Ｉ線形成分をキャンセルするために、ＩＱ信号の振幅は、変更されることになる。なお、フィルタ処理部１２１ｂは、周波数領域又は時間領域のいずれかにおいて、フィルタ処理を行なう。
【００４２】
加減算部１２１ｄは、フィルタ処理部１２１ｂによるフィルタ処理後の複数の受信信号を、変調法に応じた加減算により合成した合成信号を出力する。例えば、加減算部１２１ｄは、ＡＭ法でＣＨＩを行う場合、直交検波部１２１ａから入力された小振幅送信レートのＩＱ信号と、フィルタ処理部１２１ｂから入力された補正された大振幅送信レートのＩＱ信号とを加算して合成信号を生成する。なお、後述するように、フィルタ係数には、２つの反射波データの振幅の比率及び位相極性が加味されているので、加減算部１２１ｄは、ＡＭ法でＣＨＩを行う場合、２つの反射波データを加算して合成信号を生成する。また、例えば、加減算部１２１ｄは、ＡＭＰＭ法でＣＨＩを行う場合、直交検波部１２１ａから入力された小振幅送信レートのＩＱ信号と、フィルタ処理部１２１ｂから入力された補正された大振幅送信レートのＩＱ信号とを加算して合成信号を生成する。これにより、加減算部１２１ｄは、組織由来の線形信号をキャンセルし、造影剤由来の高調波成分を抽出することができる。加減算部１２１ｄにより生成される合成信号は、造影像を生成するための「Ｂモードデータ」となる。
【００４３】
そして、画像生成部１４０は、Ｂモードデータから、Ｂモード画像データ（造影像データ）を生成し、モニタ３０に造影像を表示する。
【００４４】
続いて、フィルタ係数の設計方法について説明する。第１の実施形態に係るフィルタ係数は、他の超音波診断装置によって、造影剤非存在下の他の被検体、造影剤非存在下の被検体Ｐ、及びファントムを用いて設計される。すなわち、第１の実施形態に係るフィルタ係数は、ファントム、又は、生体から受信した受信信号である複数の非造影受信信号を用いて設計される。なお、この他の超音波診断装置はＢモード処理ユニットの構成が異なる点を除いて、超音波診断装置１と同様の構成を有する。具体的には、他の超音波診断装置は、Ｂモード処理ユニット内に、フィルタ係数設計部を有する。以下では、このフィルタ係数設計部によるフィルタ係数設計処理について説明する。なお、フィルタ係数は、ワークステーションによって設計されてもよい。
【００４５】
図４は、第１の実施形態に係るフィルタ係数設計処理を説明するための図である。図４は、大振幅送信レートの送信信号のＩＱ信号を画像化した一例を示す。図４中に示す４ａは、大振幅送信レートの送信信号のＩＱ信号が飽和した領域の一例を示し、図４中に示す４ｂは、大振幅送信レートの送信信号のＩＱ信号が飽和していないが組織由来の線形信号が生じている領域（不飽和領域とも言う）の一例を示す。フィルタ係数を設計する場合、造影剤からの信号領域や信号レベルが飽和している部分を避けて反射波データをサンプリングする必要がある。このため、第１の実施形態に係るフィルタ係数設計部は、操作者により不飽和領域の設定を受付ける。例えば、フィルタ係数設計部は、操作者から図４中の４ｂで示す矩形領域をサンプリングする領域に設定する指示を受付ける。そして、フィルタ係数設計部は、造影剤が存在しない状態のビーム加算後の信号エネルギーが最小となるようなフィルタを設計する。
【００４６】
以下では、フィルタ係数設計部によるフィルタの設計について説明する。なお、説明の便宜上、送信信号の送信レートが２であるＡＭ法でＣＨＩを行う場合を例にフィルタ設計例を説明する。また、フィルタ係数は、線形成分のキャンセル性能のみを向上させるフィルタであるものとする。このため、フィルタ係数設計部は、非線形成分の影響を無視してフィルタを設定する。
【００４７】
送受信ユニット１１０が送信する小振幅送信レートの送信信号をｓＴｘＬ（ｔ）、大振幅送信レートの送信信号をｓＴｘＨ（ｔ）とする。これらの送信信号が、物体（反射係数：ρ（ｚ＝２ｔｃ０）、ｃ０は伝搬速度）で反射された場合に得られる、小振幅送信レートの送信信号に対する受信信号（ＲＦ信号又はＩＱ信号）をｒＴｘＬ（ｔ）とし、大振幅送信レートの送信信号に対する受信信号（ＲＦ信号又はＩＱ信号）をｒＴｘＨ（ｔ）とする。また、小振幅送信レートの送信信号の送信で受信されるノイズ(白色ノイズ)をｎＴｘＬ（ｔ）とし、大振幅送信レートの送信信号の送信で受信されるノイズ(白色ノイズ)を、ｎＴｘＨ（ｔ）とする。なお、「ｔ」は、時間を示し、深さ方向に沿った各サンプル点の位置は、「ｔ」により表される。
【００４８】
フィルタ処理後に加算して得られる合成信号（「ＡＭ信号」とも言う）ｒＡＭ（ｔ）は、式（１）で表される。ここで、ｈ（ｔ）は、フィルタのインパルス応答関数である。なお、ｈ（ｔ）には、２つの反射波データの振幅の比率及び位相極性が加味される。
【００４９】
【数１】
2015016251.tif 000003
【００５０】
この式（１）の両辺を周波数成分にフーリエ変換すると、式（２）が得られる。
【００５１】
【数２】
2015016251.tif 000004
【００５２】
この加算後の合成信号は、消え残りとノイズである。このため、フィルタ係数設計部は、加算後の合成信号のエネルギーを最小とするフィルタＨ（ω）を求める。加算後のＡＭ信号ＲＡＭ（ω）の絶対値２乗のアンサンブル平均であるε（ω）は、式（３）で表される。また、ここで白色ノイズの性質を式（４）で表す。また散乱係数の性質は、式（５）で表される。なお、「ω」は、周波数を示す。
【００５３】
【数３】
2015016251.tif 000005
【００５４】
ε（ω）は、ノイズとは無相関とする場合、式（６）で表される。
【００５５】
【数４】
2015016251.tif 000006
【００５６】
式（６）において、第２項以降は定数であるので、ε（ω）を最小とするＨ（ω）は第１項を最小とするＨＯｐｔ（ω）で表される（式（７））。
【００５７】
【数５】
2015016251.tif 000007
【００５８】
ここで、実際の送信信号スペクトルｓＴｘＬ（ω）やｓＴｘＨ（ω）が既知であれば、式（７）に従って最適なフィルタが得られる。しかし、状況によっては一様な散乱体群からの反射信号しか得られない場合が考えられる。この場合、準最適なフィルタを求める方法が必要である。なお、非線形伝播の発生が低減される音圧の超音波（例えば、低中音圧の超音波）を送信させることにより、受信信号には組織由来の非線形成分が含まれていないものとする。
【００５９】
大振幅送信レートの送信信号の送信で得られる受信信号ｒＴｘＨ（ｔ）は、式（８）で表される。この式（８）の両辺をフーリエ変換すると、式（９）が得られる。
【００６０】
【数６】
2015016251.tif 000008
【００６１】
この式（９）のパワースペクトルのアンサンブル平均を求めると、式（１０）が得られる。この式（１０）は、式（７）の分母である。
【００６２】
【数７】
2015016251.tif 000009
【００６３】
次に、小振幅送信レートの送信信号の送信時に得られる受信信号ｒＴｘＬ（ｔ）は、式（１１）で表される。この式（１１）の両辺をフーリエ変換すると、式（１２）が得られる。
【００６４】
【数８】
2015016251.tif 000010
【００６５】
ここで、大振幅送信レートの送信信号の送信時に得られる受信信号のスペクトル、式（９）の複素共役と、この小振幅送信レートの送信信号の送信時に得られる受信信号のスペクトル、式（１２）との積は、式（１３）で表される。そして、このアンサンブル平均は、式（４）〜式（５）の性質を用いて、式（１４）で表される。
【００６６】
【数９】
2015016251.tif 000011
【００６７】
この式（１４）は、式（７）の分子である。よって、ある走査線ｉの小振幅送信レートの送信信号に対する受信信号及び大振幅送信レートの送信信号に対する受信信号のサンプルスペクトルＲＴｘＬｉ（ω），ＲＴｘＨｉ（ω）を各々求めると、フィルタ係数は、式（１５）で与えられる。
【００６８】
【数１０】
2015016251.tif 000012
【００６９】
このようにして、フィルタ係数設計部は、一つの走査線に対してフィルタ係数を算出する。また、フィルタ係数設計部は、同様にして各走査線でフィルタ係数を算出する。そして、フィルタ係数設計部は、算出したフィルタ係数を平均処理する。これにより、超音波診断装置１は、ロバスト性を向上させることができる。なお、フィルタ係数設計部は、一つの走査線に対してだけで算出することにより、フィルタ係数を設計するようにしてもよい。
【００７０】
図５は、第１の実施形態に係るフィルタ係数設計部によるフィルタ係数設計処理を説明するための図である。図５では、フィルタ係数設計部によって設計されたフィルタ係数で補正した場合の大振幅送信レートの送信信号のＩＱ信号を画像化した一例を示す。なお、図５に示す撮影領域は、図４に示した撮影領域と同一である。図５に示すように、大振幅送信レートの送信信号のＩＱ信号にフィルタ係数で補正を行った場合、図４中の４ａに示した組織由来の線形信号がキャンセルされる。
【００７１】
ここで、フィルタのフィルタ長（カーネル長、フィルタ係数長、或いはタップ長とも言う）は、送信超音波のパルス長の略２倍となるように設計される。すなわち、フィルタ係数設計部は、フィルタのフィルタ長さをパルス長の約２倍程度に設計する。これにより、補正対象である受信信号の波形の先頭が、フィルタの中心に位置した場合、補正対象である受信信号の波形の末端がフィルタ内に収まるようになる。なお、フィルタ長さが長すぎると、組織由来の線形信号をキャンセルできるが、深さ方向の空間分解能が犠牲になってしまう場合がある。また、フィルタ長さが短すぎると、組織由来の線形信号をキャンセルできない。フィルタ係数設計部は、様々な送信条件ごとにフィルタ係数を設計する。
【００７２】
フィルタ係数テーブル１２１ｃは、様々な送信条件ごとに設計されたフィルタ係数を記憶する。そして、フィルタ処理部１２１ｂは、撮影時には、送信条件に合致するフィルタ係数或いは、送信条件と近似するフィルタ係数をフィルタ係数テーブル１２１ｃから取得して、フィルタに設定する。
【００７３】
次に、図６〜図８を用いて、超音波診断装置１による処理の手順について説明する。図６は、第１の実施形態に係る超音波診断装置１による処理の手順を示すフローチャートである。図６に示すように、制御部１９０は、造影剤の投与を受付ける（ステップＳ１０１）。そして、送受信ユニット１１０は、超音波の送受信処理を実行する（ステップＳ１０２）。なお、この超音波の送受信処理の詳細については、図７を用いて後述する。
【００７４】
続いて、Ｂモード処理ユニット１２０は、Ｂモードデータ生成処理を実行する（ステップＳ１０３）。なお、このＢモードデータ生成処理の詳細については、図８を用いて後述する。そして、画像生成部１４０は、造影画像データを生成し（ステップＳ１０４）、モニタ３０に造影画像を表示させる（ステップＳ１０５）。
【００７５】
図７は、第１の実施形態に係る超音波診断装置１による超音波の送受信処理の手順を示すフローチャートである。なお、この処理は、図６に示したステップＳ１０２の処理に対応する。図７に示すように、送受信ユニット１１０は、被検体Ｐに小振幅送信レートの超音波を送信する（ステップＳ２０１）。そして、送受信ユニット１１０は、小振幅送信レートの反射波信号を受信して、反射波データを生成する（ステップＳ２０２）。また、送受信ユニット１１０は、被検体Ｐに大振幅送信レートの超音波を送信する（ステップＳ２０３）。そして、送受信ユニット１１０は、大振幅送信レートの反射波信号を受信して、反射波データを生成する（ステップＳ２０４）。なお、送受信ユニット１１０は、ステップＳ２０３とステップＳ２０４とを実行してから、ステップＳ２０１とステップＳ２０２とを実行してもよい。図７に示すフローチャートは、１フレーム分の複数の走査線について行われる。
【００７６】
図８は、第１の実施形態に係る超音波診断装置１によるＢモードデータ生成処理の手順を示すフローチャートである。なお、この処理は、図６に示したステップＳ１０３の処理に対応する。Ｂモード処理ユニット１２０において、直交検波部１２１ａは、送受信ユニット１１０から入力された受信信号の直交検波を行う（ステップＳ３０１）。直交検波部１２１ａは、小振幅送信レートの反射波データ（ＩＱ信号）を加減算部１２１ｄに出力し、大振幅送信レートの反射波データ（ＩＱ信号）をフィルタ処理部１２１ｂに出力する。そして、フィルタ処理部１２１ｂは、フィルタ係数テーブル１２１ｃから送信条件と合致する或いは送信条件と近似するフィルタ係数を読出して、大振幅送信レートの反射波データ（ＩＱ信号）を補正する（ステップＳ３０２）。そして、加減算部１２１ｄは、小振幅送信レートの反射波データ（ＩＱ信号）と、補正された大振幅送信レートの反射波データ（ＩＱ信号）とを加減算処理を行う（ステップＳ３０３）。図８に示すフローチャートは、１フレーム分の複数の走査線について行われる。これにより、１フレーム分の造影像用のＢモードデータが生成される。
【００７７】
上述したように、第１の実施形態に係る超音波診断装置１は、受信信号に対して、組織由来の線形信号をキャンセルするフィルタ係数が設定されたフィルタを有し、ＣＨＩを行う場合、複数の受信信号の少なくとも１つの受信信号をこのフィルタで処理する。これにより、超音波診断装置１は、組織由来の線形信号をキャンセルし、造影剤由来の高調波成分を抽出することができる。この結果、超音波診断装置１は、バブル組織比の高い造影画像を生成することができる。
【００７８】
また、複数の非造影受信信号それぞれは、ファントム、又は、生体から受信した受信信号であって、信号レベルが飽和していない不飽和領域からの受信信号である。これにより、フィルタ係数設計部は、バブル組織比の高いフィルタ係数を設計することが可能となる。
【００７９】
第１の実施形態では、非線形伝播の発生が低減される音圧の超音波送信により受信された受信信号を用いてフィルタ係数を設計する。このため、第１の実施形態では、組織由来の非線形成分が含まれていない条件でフィルタ係数を設計するので、組織由来の線形成分をキャンセル可能であるとともに、造影剤由来の高調波成分だけを抽出することが可能なフィルタを設計することができる。
【００８０】
フィルタ係数テーブル１２１ｃは、複数の送信条件ごとに予め設計された複数のフィルタ係数を記憶する。これにより、操作者が造影撮影用の送信条件を入力するだけで、入力された送信条件に適合するフィルタ係数が自動的に選択される。この結果、バブル組織比の高い造影画像を容易に生成することができる。
【００８１】
（第２の実施形態）
第１の実施形態では、フィルタ係数が事前に設定されている場合について説明した。しかし、より高画質な造影像を生成するには、事前に設計されたフィルタ係数を用いてフィルタ処理を行うよりも、撮影部位に応じて適応的に設計されたフィルタ係数を用いてフィルタ処理を行う方が望ましい。このため、第２の実施形態では、超音波診断装置１が、スキャンする被検体Ｐからの受信信号に基づき適応的にフィルタ係数を設計する例を説明する。
【００８２】
第２の実施形態に係る超音波診断装置１の構成は、Ｂモード処理ユニットの構成が一部異なる点を除いて、図１に示した第１の実施形態に係る超音波診断装置１の構成と同様である。図９は、第２の実施形態に係るＢモード処理ユニット１２０の構成例を示すブロック図である。図９に示すように、第２の実施形態に係るＢモード処理ユニット１２０は、直交検波部１２１ａと、フィルタ処理部１２１ｂと、加減算部１２１ｄと、フィルタ係数設計部１２１ｅとを有する。
【００８３】
第２の実施形態に係るフィルタ係数設計部１２１ｅは、複数の非造影受信信号に基づいて、フィルタ係数を設計する。すなわち、第２の実施形態に係るフィルタ係数設計部１２１ｅは、第１の実施形態で説明したフィルタ係数設計部と同様の機能を有する。
【００８４】
また、フィルタ係数設計部１２１ｅが用いる複数の非造影受信信号それぞれは、造影撮影と同一の送信条件により受信される受信信号であって、造影撮影が行なわれる被検体Ｐの撮影部位であり、造影剤が存在しない状態の撮影部位から受信した受信信号である。具体的には、複数の非造影受信信号それぞれは、造影剤が到達していない撮影部位又造影剤投与前の撮影部位から受信した受信信号である。ここで、例えば、フィルタ係数設計部１２１ｅは、操作者の入力を契機として、フィルタ係数を設計する。具体的には、造影剤投与前に、撮像部位の超音波走査を開始したタイミングでフィルタ係数を設計する。
【００８５】
また、フィルタ係数設計部１２１ｅは、複数の非造影受信信号の信号レベルに基づいて、当該複数の非造影受信信号が不飽和領域からの受信信号であるか否かを判定する。例えば、フィルタ係数設計部１２１ｅは、入力された受信信号のエネルギーと、ダイナミックレンジとに基づいて、受信信号が飽和しているか否かを判定する。具体的には、フィルタ係数設計部１２１ｅは、ダイナミックレンジよりも大きなエネルギーの受信信号が入力された場合に、受信信号が飽和していると判定する。また、フィルタ係数の設計を行う際に、フィルタ係数設計部１２１ｅは、例えば、大振幅送信レートの受信信号を用いて、受信信号が飽和しているか否かを判定する。なお、フィルタ係数設計部１２１ｅは、例えば、小振幅送信レートの受信信号を用いて、受信信号が飽和しているか否かを判定してもよい。また、フィルタ係数設計部１２１ｅは、例えば、大振幅送信レートの受信信号と、小振幅送信レートの受信信号とを用いて、受信信号が飽和しているか否かを判定してもよい。
【００８６】
また、フィルタ係数設計部１２１ｅは、フレーム内の全ての走査線でフレーム内のフィルタ係数を設計するようにしてもよいし、フレーム内の一部の複数の走査線でフレーム内のフィルタ係数を設計するようにしてもよいし、フレーム内の一つの走査線でフレーム内のフィルタ係数を設計するようにしてもよい。また、これらの処理を複数のフレームで行ってもよい。また、フィルタ係数設計部１２１ｅは、操作者によりＲＯＩ（Region  Of  Interest）の設定を受付けて、ＲＯＩ内の全ての走査線でＲＯＩ内のフィルタ係数を設計するようにしてもよいし、ＲＯＩ内の一部の複数の走査線でＲＯＩ内のフィルタ係数を設計するようにしてもよいし、ＲＯＩ内の一つの走査線でＲＯＩ内のフィルタ係数を設計するようにしてもよい。また、これらの処理を複数のフレームで行ってもよい。
【００８７】
第２の実施形態に係るフィルタ処理部１２１ｂは、フィルタ係数設計部１２１ｅが設計したフィルタ係数をフィルタに設定する。そして、フィルタ処理部１２１ｂは、複数の受信信号の少なくとも１つの受信信号に対して、フィルタをかける。加減算部１２１ｄは、フィルタ処理後の複数の受信信号を、変調法に応じた加減算により合成した合成信号（Ｂモードデータ）を出力する。そして、画像生成部１４０は、Ｂモードデータから、信号強度が輝度の明るさで表現されるＢモード画像を生成する。
【００８８】
図１０は、第２の実施形態に係る超音波診断装置１による処理の手順を示すフローチャートである。図１０に示すように、フィルタ係数設計部１２１ｅは、フィルタ設計処理を実行する（ステップＳ４０１）。なお、このフィルタ設計処理の詳細については、図１１を用いて後述する。そして、制御部１９０は、造影剤の投与を受付ける（ステップＳ４０２）。そして、送受信ユニット１１０は、超音波の送受信処理を実行する（ステップＳ４０３）。なお、この超音波の送受信処理の手順は、図７に示した処理の手順と同様である。
【００８９】
続いて、Ｂモード処理ユニット１２０は、Ｂモードデータ生成処理を実行する（ステップＳ４０４）。なお、このＢモードデータ生成処理の手順は、図８に示した処理の手順と同様である。そして、画像生成部１４０は、造影画像データを生成し（ステップＳ４０５）、モニタ３０に造影画像を表示させる（ステップＳ４０６）。
【００９０】
なお、超音波診断装置１は、投与した造影剤が流出され切った後に、再度造影剤を投与して、異なる送信条件で造影撮影を行う場合には、ステップＳ４０１に移行し、ステップＳ４０１以降の処理を繰り返し実行するようにしてもよい。
【００９１】
図１１は、第２の実施形態に係る超音波診断装置１によるフィルタ設計処理の手順を示すフローチャートである。図１１に示すように、フィルタ係数設計部１２１ｅは、各レートの反射波データを受信する（ステップＳ５０１）。例えば、フィルタ係数設計部１２１ｅは、直交検波部１２１ａから小振幅送信レートの反射波データと、大振幅送信レートの反射波データとを受信する。
【００９２】
そして、フィルタ係数設計部１２１ｅは、フィルタ係数を設計可能な信号か否かを判定する（ステップＳ５０２）。ここで、フィルタ係数設計部１２１ｅは、フィルタ係数を設計可能な信号であると判定しなかった場合（ステップＳ５０２、Ｎｏ）、ステップＳ５０１に移行して、各レートの反射波データを受信する。なお、フィルタ係数設計部１２１ｅは、フィルタ係数を設計可能な信号であると判定しなかった場合、例えば、プローブの当接位置を変更する旨や送信条件を変更する旨をモニタ３０に表示させるように制御部１９０に通知してもよい。これにより、再度受信した受信信号によりフィルタ係数を設計可能な信号か否かを判定する。
【００９３】
一方、フィルタ係数設計部１２１ｅは、フィルタ係数を設計可能な信号であると判定した場合（ステップＳ５０２、Ｙｅｓ）、フィルタ係数を設計する（ステップＳ５０３）。例えば、フィルタ係数設計部１２１ｅは、小振幅送信レートの送信信号の反射波データ（ＩＱ信号）と、大振幅送信レートの送信信号の反射波データ（ＩＱ信号）とからフィルタ係数を算出する。そして、フィルタ係数設計部１２１ｅは、算出したフィルタ係数をフィルタ処理部１２１ｂに通知する。
【００９４】
上述したように、第２の実施形態に係る超音波診断装置１は、被検体Ｐと送信条件とに適応したフィルタ係数を設計する。ここで、複数の非造影受信信号それぞれは、造影撮影と同一の送信条件により受信される受信信号であって、造影撮影が行なわれる被検体の撮影部位であり、造影剤が存在しない状態の撮影部位から受信した受信信号である。このため、超音波診断装置１は、バブル組織比の高い造影画像を生成することができる。
【００９５】
なお、造影剤を投与後、撮影対象部位まで造影剤が届くまでには、約２０〜３０秒かかる。このため、第２の実施形態では、造影剤投与前ではなく造影剤投与後に、撮影対象部位まで造影剤が届くまでの時間を利用して、フィルタ係数を設計してもよい。かかる場合、この処理には、通常造影撮影に用いられる造影剤投与タイマーを利用する。制御部１９０は、この造影剤投与タイマーから造影剤投与後の経過時間を取得する。そして、例えば、造影剤投与から１０秒後に自動的に造影剤投与後の非造影受信信号群からフィルタ係数を設計する指示をフィルタ係数設計部１２１ｅに送出する。すなわち、フィルタ係数設計部１２１ｅは、造影剤投与の時間を計測する造影剤投与タイマーに連動して、フィルタ係数を設計する。これにより、操作者は、造影剤投与前にフィルタ係数の設計指示を入力しなくてもよく、自動的にフィルタ係数を設計することができる。
【００９６】
また、第２の実施形態において、Ｂモード処理ユニット１２０は、フィルタ係数テーブル１２１ｃを有するようにしてもよい。この場合、フィルタ係数設計部１２１ｅは、適応的に設計したフィルタ係数をフィルタ係数テーブル１２１ｃに記憶させる。これにより、フィルタ処理部１２１ｂは、同一の被検体かつ同一の送信条件で設計したフィルタ係数がフィルタ係数テーブル１２１ｃに記憶されている場合、フィルタ係数テーブル１２１ｃからフィルタ係数を読出して、フィルタに設定することができる。
【００９７】
なお、第１の実施形態で説明した内容は、造影剤到達前に、造影撮影に用いる条件（撮影条件）で得られた非造影受信信号群を用いて、適応型のフィルタ係数を設計する点以外、第２の実施形態においても適用可能である。
【００９８】
（第３の実施形態）
第２の実施形態では、スキャンする被検体Ｐからの受信信号に基づき適応的にフィルタを設計する例を説明した。ところで、造影剤投与後に、送信条件を変えて造影像を取得したい場合がある。かかる場合、変更後の送信条件に適したフィルタ係数を設計することが望ましい。しかし、造影剤投与後であるので、被検体Ｐ内には造影剤バブルが残ったままであり、フィルタ係数を設定することができない。このため、第３の実施形態では、超音波診断装置１は、残存する造影剤バブルを壊すために超音波を送信して造影剤バブルが存在しない状態で、組織由来の線形信号をキャンセルできるフィルタを設計する。以下では、造影剤を破壊可能な音圧の超音波のことを「フラッシュ」と呼ぶ。
【００９９】
第３の実施形態に係る制御部１９０は、「フラッシュ用のボタン」の選択を操作者から受付けた場合、送受信ユニット１１０に造影剤を破壊可能な音圧の超音波（フラッシュ）を送信させるように指示する。
【０１００】
第３の実施形態に係るフィルタ係数設計部１２１ｅは、第２の実施形態に係るフィルタ係数設計部１２１ｅと同様の機能を有することに加えて、更に以下の機能を有する。例えば、第３の実施形態に係るフィルタ係数設計部１２１ｅは、造影剤投与後に造影剤を破壊可能な音圧の超音波送信が行なわれた撮影部位から受信した受信信号である複数の非造影受信信号それぞれを用いてフィルタ係数を設計する。
【０１０１】
第３の実施形態に係るフィルタ処理部１２１ｂは、フィルタ係数設計部１２１ｅが設計したフィルタ係数をフィルタに設定する。そして、フィルタ処理部１２１ｂは、第２の実施形態と同様に、複数の受信信号の少なくとも１つの受信信号に対して、フィルタをかける。加減算部１２１ｄは、第２の実施形態と同様に、フィルタ処理後の複数の受信信号を、変調法に応じた加減算により合成した合成信号（Ｂモードデータ）を出力する。そして、画像生成部１４０は、第２の実施形態と同様に、Ｂモードデータから、信号強度が輝度の明るさで表現されるＢモード画像を生成する。
【０１０２】
図１２は、第３の実施形態に係る超音波診断装置１による処理の手順を示すフローチャートである。図１２に示すように、制御部１９０は、フラッシュを予め設定されたフレーム分送信する（ステップＳ６０１）。そして、フィルタ係数設計部１２１ｅは、フィルタ設計処理を実行する（ステップＳ６０２）。なお、このフィルタ設計処理の手順は、図１１に示した処理の手順と同様である。
【０１０３】
そして、送受信ユニット１１０は、超音波の送受信処理を実行する（ステップＳ６０３）。なお、この超音波の送受信処理の手順は、図７に示した処理の手順と同様である。続いて、Ｂモード処理ユニット１２０は、Ｂモードデータ生成処理を実行する（ステップＳ６０４）。なお、このＢモードデータ生成処理の手順は、図８に示した処理の手順と同様である。そして、画像生成部１４０は、造影画像データを生成し（ステップＳ６０５）、モニタ３０に造影画像を表示させる（ステップＳ６０６）。なお、超音波診断装置１は、更に送信条件を変更して造影撮影を行う場合、ステップＳ６０１に移行して、フラッシュを送信し、フィルタ係数を設計する処理を実行する。そして、超音波診断装置１は、フィルタ係数を設計した後、ステップＳ６０２以降の処理を実行する。
【０１０４】
上述したように、第３の実施形態に係る超音波診断装置１は、造影剤バブルを崩壊する超音波送信を行うことで、造影撮影用の送信条件を適宜変更する場合であっても、適応的なフィルタ係数を設計することが可能である。したがって、第３の実施形態では、造影撮影用の送信条件が変更された場合であっても、バブル組織比の高い造影画像を生成することができる。
【０１０５】
なお、第１の実施形態及び第２の実施形態で説明した内容は、フラッシュ送信によりフィルタ係数を再設計する点以外、第３の実施形態においても適用可能である。
【０１０６】
（その他の実施形態）
上述した実施形態においては、ＡＭ法でＣＨＩを行う場合を例に説明したが、これに限定されるものではない。例えば、ＡＭＰＭ法でＣＨＩを行う場合でも適用可能である。更に、ＡＭＰＭ法では、ＴＨＩを行う場合にも適用可能である。
【０１０７】
また、上述した実施形態においては、フィルタ処理部１２１ｂは、大振幅送信レートのＩＱ信号にフィルタ処理を行う場合を説明したが、これに限定されるものではない。例えば、フィルタ処理部１２１ｂは、小振幅送信レートのＩＱ信号にフィルタ処理を行ってもよい。或いは、フィルタ処理部１２１ｂは、大振幅送信レートのＩＱ信号と、小振幅送信レートのＩＱ信号とにフィルタ処理を行ってもよい。
【０１０８】
また、上述した実施形態においては、フィルタ処理部１２１ｂは、直交検波後のＩＱ信号に対してフィルタ処理を行う例を説明したが、これに限定されるものではない。例えば、フィルタ処理部１２１ｂは、ＲＦ信号に対してフィルタ処理を行うようにしてもよい。すなわち、フィルタ処理部１２１ｂは、フィルタをＩＱ信号又はＲＦ信号に対してかけるようにしてもよい。なお、ＩＱ信号を処理する場合のフィルタは、複素ＦＩＲ（有限インパルス応答）フィルタであり、ＲＦ信号を処理する場合のフィルタは、実数ＦＩＲ（有限インパルス応答）フィルタである。
【０１０９】
フィルタ係数設計部１２１ｅは、１画像（１フレーム）だけではなく、数フレーム分のフィルタ係数を算出して、各走査線で平均したフィルタ係数を設計するようにしてもよい。また、フィルタ係数設計部１２１ｅは、１フレームを短冊状の領域に分割し、分割した領域それぞれでフィルタ係数を設計するようにしてもよい。また、フィルタ係数設計部１２１ｅは、深さ方向に応じてフィルタ係数を設計するようにしてもよい。これにより、超音波診断装置１は、マルチフォーカスを行う場合に、組織由来の線形信号をキャンセルし、バブル組織比の高い造影画像を生成することができる。また、フィルタ係数設計部１２１ｅは、撮影領域をメッシュ状に分割し、分割した領域に含まれる複数のピクセル単位でフィルタ係数を算出するようにしてもよい。
【０１１０】
以上説明した少なくともひとつの実施形態によれば、バブル組織比の高い造影画像を生成することができる。
【０１１１】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【０１１２】
１  超音波診断装置
１００  装置本体
１１０  送受信ユニット
１２０  Ｂモード処理ユニット
１２１ｂ  フィルタ処理部
１２１ｄ  加減算部
【図１】
2015016251.tif 000013
【図２】
2015016251.tif 000014
【図３】
2015016251.tif 000015
【図４】
2015016251.tif 000016
【図５】
2015016251.tif 000017
【図６】
2015016251.tif 000018
【図７】
2015016251.tif 000019
【図８】
2015016251.tif 000020
【図９】
2015016251.tif 000021
【図１０】
2015016251.tif 000022
【図１１】
2015016251.tif 000023
【図１２】
2015016251.tif 000024