JP5162923B2 - 超音波撮像装置 - Google Patents

Description

本発明は、超音波断層像を表示する装置に関する。

Ｘ線ＣＴ、ＭＲＩ、超音波診断装置などの画像診断モダリティが医療現場で必須のツールになって久しい。これらは生体内でのＣＴ値、スピン緩和時間、音響インピーダンスの違いをそれぞれ画像化したものであり、これら物理的性質の違いが専ら生体の構造（かたち）を反映することから、「形態イメージング」と呼ばれる。これに対し、構造的には同じ組織であっても機能的に異なる状態にある部位の画像化を行うものを「機能イメージング」と呼ぶ。この機能イメージングの内、特にタンパクなどの生体構成分子の存在状態の可視化を行うものが「分子イメージング」と呼ばれることが多い。分子イメージングは発生・分化といった生命現象の解明や疾病の診断・治療への応用が期待されることから、現在最も注目を浴びている研究領域のひとつである。分子イメージングでは、生体構成分子に選択性を有する構造を有する物質である「分子プローブ」を用いることが多く、この場合には分子プローブをなんらかの物理的手段で検出可能とする構造を付加し、体内での分子プローブの分布を可視化する。

上記分子プローブの検出用の物理的手段として医療用の画像診断装置を用いることができれば、従来の形態イメージングと分子イメージングとを融合することが可能となり、高次の診断が可能になるものと期待される。超音波を用いた分子イメージングを行う上で重要な技術として、分子プローブの体内分布を超音波画像上で可視化する技術が挙げられる。一般的には、分子プローブそのものを超音波画像診断機上で感度を持たせる設計を行なうのは困難である。そのため、分子プローブと化学的あるいは物理的に結合し、分子プローブの存在状態を画像診断装置上で可視化するための造影剤と分子プローブを組み合わせて用いる方が、分子プローブ、造影剤をそれぞれ最適設計が可能となり、開発の自由度が大きくなる。現在もっぱら用いられる超音波造影剤は数ミクロンの直径を有する気泡である。超音波診断装置では、物質の音響インピーダンスすなわち密度に音速を乗じた物性値の異なる部位が可視化されるという特徴を有しており、生体の音響インピーダンス約１．５×１０^６ｋｇ／ｍ^２・ｓに比べてはるかに小さい値約０．００４×１０^６ｋｇ／ｍ^２・ｓ（空気）をもつ気泡は生体中での可視化が容易である。さらには、数ミクロンの気泡は診断用に用いられる数ＭＨｚ帯域の超音波に共振し、照射した超音波の高調波成分を生成するという性質を有することから気泡からの信号のみを選択的に可視化することが可能となり、より高感度な造影が可能となっている。

超音波を用いた分子イメージングを実現するために分子プローブと共に用いる造影剤として、上記のように既に血流診断用に用いられている気泡に分子プローブを結合した構造の薬剤を用いることは極めて自然であり、例えば〔特許文献１〕にあるように血栓選択性造影剤として開発されている。また、例えば〔非特許文献１〕にあるように新生血管選択性造影剤も開発されている。しかしながら、これらミクロンサイズの気泡を用いる造影剤では、血管以外の造影を行うことは難しく、適用範囲が限られるという欠点がある。さらに、マイクロバブルは診断レベルの超音波強度でも破壊されてしまうため一度造影されてしまうと破壊され連続して造影することが困難である。また、バブル（気体）であることから、肺におけるガス交換により体外へ排出され、血液中に十数分間しか留まれないという欠点を有する。

これに対し、例えば〔非特許文献２〕に示されるように、音響インピーダンスが生体と異なる液体をカプセル化しサイズをサブミクロンまで小さくすることにより、血管から組織への移行を可能とした造影剤に関する研究も行われている。本手法によれば、血管以外の部位への適用も可能になるものと考えられる。また、血中滞留時間も気体より長いことが期待される。ただし、微小な液体を用いるために、上述のマイクロバブル造影剤における共振現象は生じず、感度の点で不利である。

これらの手法の長所を併せた形の造影剤として、〔特許文献２〕に示されるように、液体を界面活性剤により微粒子化した化合物を投与し、体内にて超音波照射により気化させるというタイプのものが考えられている。このタイプの造影剤によれば、体内滞留時間および適用部位に関する制約が少なく、また気体の状態で可視化するため、共振を利用した高感度な造影が可能であると考えられる。

液体を界面活性剤により微粒子化した化合物を投与し、体内にて超音波照射により気化させるというタイプの造影剤場合に、気化した後の数ミクロンのマイクロバブルからの信号を特異的に画像化する方法としては、送信周波数と異なる周波数の受信エコーを用いる方法が、例えば〔特許文献３〕，〔特許文献４〕に開示されている。

米国特許第６、５２１、２１１号 米国特許第５、７１６、５９７号 特開平１１−１３７５４７号公報 特開２００３−１３５４６７号公報 Ｅｌｌｅｇａｌａら（２００３）Ｃｉｒｃｕｌａｔｉｏｎ １０８：３３６−３４１ Ｌａｎｚａら（１９９６）Ｃｉｒｃｕｌａｔｉｏｎ ９４：３３３４−３３４０

形態・分子融合イメージングを考える上で、超音波を用いたシステムを実現する意義は大きい。一つには、リアルタイム性の高さが挙げられる。形態・分子融合イメージングは疾病の極早期での検出を可能にすると考えられるが、そのような極早期においては、疾病領域は小さく、一枚の画像では疾病領域かも知れないということしか判断できない状況がままありえる。かかる状況においては、関心領域を様々な角度から観察することが重要であると考えられる。このような目的にはリアルタイム性に優れる超音波は優れたツールとなりうる。二つ目の理由としては、超音波は照射条件を変えることにより、診断と治療の両方を行うことができ、しかもいずれにおいても生体への侵襲性が極めて低いことが挙げられる。上述の如く、形態・分子融合イメージングにおいては極早期の疾病の診断を行うことが可能になると考えられ、治療する対象も極小さな限られた領域になると予想される。このような状況においては侵襲性の高い治療法はリスクが高すぎるため、低侵襲治療が望ましい。超音波を用いることにより、診断し、その場で低侵襲治療を行うことも可能になると考えられる。

しかし、超音波を用いたシステムでは、超音波探触子の帯域が有限である場合、送信波形と、造影剤からの受信波形の帯域を分けようとすると、送信もしくは受信すくなくとも一方の帯域が狭くなり、空間分解能が低下する。すなわち、造影剤からの信号と他の信号の弁別比（造影剤信号の特異度）と空間分解能を両立するのが困難であった。

本発明においては、微粒子化した液体を、一度気化させて造影剤に変えてから、超音波パルスを送波して、可視化するのではなく、微粒子化した液体が超音波によって、気化するときに放出される、造影剤特異的な信号を検出し、その空間分布を画像化する。微粒子化した液体を気化する場合、エネルギーの蓄積が必要なため、超音波波形はある程度長い必要があり、周波数帯域では狭い信号となる。この長い波形に対するエコー信号は、基本波、高調波ともに、帯域が狭い。一方、気化するときに放出される超音波は、一気に気化するために、時間軸で短い、すなわち周波数空間では広帯域な波形となる。本発明では、送波には狭い帯域を使い、受波信号では広い帯域を使うことが可能になるため、受波信号の空間分解能が高い状態で、送波信号と受波信号を区別することが可能になる。

本発明に係る超音波撮像装置は、一例として、被検体の対象領域へ第１超音波信号を送信する送信手段と、前記被検体から前記第１超音波信号の照射によって生じる第２超音波信号を受信する受信手段と、前記第２超音波信号から前記第１超音波信号の基本波パルスと高調波パルスのうちの少なくとも１つを除去し検出信号を算出する計算部とを有する。

本発明に係る超音波撮像装置は、他の例として、被検体の相変化造影剤投与領域へ第１超音波信号と前記第１超音波信号よりも波数の多い第２超音波信号とを送信する送信手段と、前記被検体から、前記第１超音波信号及び前記第２超音波信号の各々によって生じる第３超音波信号及び第４超音波信号を受信する受信手段と、前記第3超音波信号について波数を前記第2超音波信号と同一にする波数処理を行って第５超音波信号とし、前記第3超音波信号と前記第5超音波信号との差分を検出信号として算出する計算部とを有する。

造影剤特異的な信号と他の信号の弁別比（造影剤信号の特異度）と空間分解能が両立することが出来る。

まず図１を用いて、超音波診断装置における、画像化のための信号処理の流れを説明する。ここには図示しない被検体の表面に設置された超音波探触子１に対して、送受切替スイッチ２を介して、制御系４の制御のもと送信ビームフォーマ３から送波電気パルスが送られる。送受切替スイッチ２は、超音波の送信手段及び受信手段（図示せず）と超音波探触子１との接続を切り替えるものである。このとき送信ビームフォーマは、所望の走査線上に超音波ビームが進むように、探触子１の各ｃｈ間の遅延時間が適した状態になるように制御している。この送信ビームフォーマ３から電気信号を受けて、超音波探触子１において電気信号は超音波信号に変換され、被検体内に超音波パルスが送波される。被検体内において、散乱された超音波パルスは一部がエコー信号として再び超音波探触子１によって、受信され、超音波信号から電気信号に変換される。この受信された信号は送受切替スイッチ２を介して、受信ビームフォーマ２０によって、所望の走査線上の所望の深さからのエコー信号が選択的に増強された、ある走査線上のデータとなる。この走査線上のＲＦ信号は造影剤エミッションパルス抽出部２１において、造影剤特異的な信号が強調され、包絡線検波部２２において、包絡線信号に変換され、Log圧縮され、スキャンコンバータ２３に送られる。スキャンコンバータ２３上では、ビデオ信号に、スキャンコンバージョンが行われる。このスキャンコンバージョン後のデータが表示部２４に送られ、超音波断層像として、表示される。送波波形と造影剤エミッションパルス抽出部以外は、ここでは詳しい説明は省略する。

図２および３を用いて、従来方法と比較しながら、超音波照射の流れを説明する。図２（ａ）に示す従来方法では、まず被検体に造影剤を投与する。つぎに、造影剤がある程度循環したところで、関心領域における造影剤の気泡化を行なうための超音波照射を行なう。気泡化の後で、一般的なマイクロバブルを造影するシーケンスによって、撮像を行い、造影剤画像を得る。この一連のシーケンスで生成する波形を図２（ｂ）に示す。超音波探触子から送波された基本波が放射されると、気泡の存在する領域で、気泡からのエコー、周囲の組織によって反射される基本波エコー、やはり周囲の組織によって反射される高調波エコーの３種類の信号が超音波探触子に返される。この場合には、この３つの種類の受波信号から、気泡からのエコーだけを選択的に抽出する必要がある。

一方、図３（ａ）に示す本発明では、造影剤を投与したあと、気泡化と画像化を同時に行なう。この場合は、図３（ｂ）に示すように、相変化（気泡化）パルスが放射されると、造影剤の存在する領域で、気泡化時に放出される超音波（以下エミッションパルスと呼ぶ）、周囲の組織によって反射される相変化パルスに対する基本波エコー、やはり周囲の組織によって反射される高調波エコーの３種類の信号が超音波探触子に返される。エミッションパルスとは、造影剤の気泡化の際に、造影剤そのものが周囲に放射する超音波をさす。図では模式的に、この気泡化時に超音波が照射されることを矢印で示している。この３つの信号成分の区別は周波数帯域による。図４を用いて、上記３つのパルスの周波数軸上での関係を説明する。相変化パルスは、時間軸上で長めに設定される。このため、周波数空間上では狭くなる。送波パルスの帯域が狭い場合は、高調波パルスの周波数帯域も狭くなる。一方、エミッションパルスは、これと無関係に広帯域で放射される。相変化パルスの基本波パルスと、二倍高調波パルスなどの高調波パルスが、それぞれ、超音波探触子の帯域の両端（両端やその近傍）にくるように設定すると、帯域の真ん中の広い部分を使って、エミッションパルスを扱うことが出来る。

この一連の処理をフローチャートにまとめると、図５のようになる。この図５の処理が図１の造影剤エミッションパルス抽出部（計算に関わる計算部）で行なわれる。バンドパスフィルタの実現方法はいくつかあるが、周波数移動とローパスフィルタであるハニング関数を組み合わせて実現する方法を説明する。ローパスフィルタとして、ハニング関数をもちいると、図８に示すように、ローパスフィルタでの所望の帯域幅に対して、ハニング関数の幅が定まる。周波数移動は受信信号に対して、移動周波数ｆのｓｉｎ波をかけることで実現できる。これは受信信号をｓｉｎ（ｆ１ｘｔ）とすると、これにｓｉｎ（ｆｘｔ）を掛けると、三角関数の積和定理より、ｓｉｎ（ｆ１ｘｔ）ｘｓｉｎ（ｆｘｔ）＝−０．５ｘ｛ｃｏｓ（（ｆ１＋ｆ）ｘｔ）−ｃｏｓ（（ｆ１−ｆ）ｘｔ）｝となる。この信号に前記のローパスフィルタをかけると、高周波成分ｃｏｓ（（ｆ１＋ｆ）ｘｔ）は除去されるので、ｆ１−ｆの周波数の信号になる。つまりｆだけ周波数が下がった信号になる。 このように、バンドパスフィルタの中心周波数ｆと比帯域幅を決めると、受信信号に対して設定した周波数帯域上の窓部分のみの信号を抽出することが出来るので、気泡化に用いた相変化パルス及びこの高調波成分と、造影剤気化時に送波されるエミッションパルスを弁別することが出来る。

以下、本構成の優位性を図６に示す従来方法と比較して説明する。（A）のケースは送受信ともにハニング関数を包絡線とするｓｉｎ３波の場合で、図8から明らかなように比帯域幅６０％となる。この場合、現在の超音波探触子の比帯域幅は７０％程度であるので、送信と受信で帯域がほとんど重なってしまう。例えば探触子の帯域が2ＭＨｚから4.2ＭＨｚの場合、送波パルスに比帯域幅６０％を使おうとすると、２ＭＨｚから３．７ＭＨｚとなり、探触子の帯域のほぼすべてを使い切ってしまう。そのため受波を探触子の帯域内で使おうとすると、２．２５ＭＨｚから４．２ＭＨｚで比帯域幅６０％となるので、2.25ＭＨｚから3.7ＭＨｚの範囲で送波と受波の帯域が重なってしまう。

一方、重ならないように、受波の帯域を設定しようとすると、探触子の帯域からはみでてしまうので、結局受波の感度が悪くなってしまう。（B）のケースでは送信、受信の比帯域幅を25％とした場合で、図８からｓｉｎ８波となる。比帯域幅の合計は５０％なので、探触子の帯域７０％内に収まるが、パルス波形が8波もあると、通常の診断装置に比べ、空間分解能が3倍程度悪くなってしまう。（現状の診断装置においては、ｓｉｎ２から3波程度）図6の（A）のケースでは、空間分解能は良いが、造影剤からの信号と組織からの信号の弁別比が悪くなる。

一方、（B）のケースでは、造影剤からの信号と組織からの信号の弁別比は良いが、空間分解能が悪くなる。超音波探触子の-6dB比帯域幅は通常６０から８０％と有限である。（尚、周波数空間において、最大の感度から-6dB感度の範囲の下限周波数をｆ１、上限周波数をｆ２とした場合（図１１）、-6dB比帯域幅は、（ｆ２−ｆ１）/（（ｆ１＋ｆ２）/２）が-6dB比帯域幅である。）（既に気泡化した造影剤からのエコー信号は、単純な数式で表現するのは難しいが、励起信号（すなわち送波パルス）が時間軸上で長いと、造影剤からの信号も長くなるので、（B）の図のように、送波と受波で帯域を区別するために、送波波形を時間軸上で長くすると、受波波形も時間軸上で長くなってしまい、空間分解能が劣化する。

一方（A）の図に示すように、受波信号を時間軸で、短くすることを優先して、送波パルスも時間軸で短くすると、探触子の帯域幅が有限のために、送波と受波の周波数帯域が重なり、区別が困難となってしまう。このような方法で、気泡化した後に造影画像を取得した場合では、造影剤からの信号と組織からの信号の弁別比と、空間分解能の間にトレードオフが存在し、両立させることが困難である。なお、信号ＡのＢに対する弁別比とは、「二つの信号の帯域の重なりの積分」／「信号Ａの帯域の積分」として定義できる。よって、周波数帯域での信号の重なりを小さくすることが弁別比の向上に必要である。本発明においては、気泡化する瞬間に放出されるエミッションパルスは、送波波形のパルスの長さに関わらず、時間軸では短くなるために、弁別比と、空間分解能との双方を高めることができる。このエミッションパルスが時間軸で短くなるのは、気化が起きる時間の長さが、超音波の周期（周波数の逆数）に比べ十分に短いからである。

これまで、パルス長について、長パルス、短パルスと説明してきた。以下、具体的なデータを用いて、定量的な長パルス、短パルスの説明を行なう。

まず、気泡化に必要なエネルギーの観点から、実験データを用いて説明する。図７は、アクリルアミドゲル中に固定した、パーフロロペンタンをリン脂質（フォスファチジルコリン）及びコレステロールでエマルション化した直径数１００ｎｍ程度の液滴である造影剤に対して超音波を照射して、Ｂ−ｍｏｄｅ画像を取得し、当該画像を用いて気泡化による送波音圧を評価した結果である。(a)は、2.226MHｚの超音波を照射し、光学顕微鏡画像によって、気泡化の閾値電圧を判定したときのパルス波数としきい値の関係を表すものである。（ｂ）は、6.7MHｚもしくは2.226MHｚの超音波を照射し、光学顕微鏡画像によって、気泡化の閾値電圧を判定したときのパルス波数としきい値の関係を表すものである。

菱形のドットは造影剤の構成要素である、パーフルオロペンタンとパーフルオロヘプタンの配分比率が１００：０の場合の結果を、四角のドットは造影剤の構成要素である、パーフルオロペンタンとパーフルオロヘプタンの配分比率が７５：２５の場合の結果を、三角のドットは造影剤の構成要素である、パーフルオロペンタンとパーフルオロヘプタンの配分比率が５０：５０の場合の結果を各々表す。（ｂ）の四角ドットは超音波の周波数が6.7MHzであり、造影剤の構成要素である、パーフルオロペンタンとパーフルオロヘプタンの配分比率が７５：２５の場合の結果を、三角ドットは、超音波の周波数が2.226MHzであり、やはり、造影剤の構成要素である、パーフルオロペンタンとパーフルオロヘプタンの配分比率が７５：２５の場合の結果を示している。図７にしめされるように、気泡化閾値はパルス波長をながくしていくとある変曲点を境に閾値が下がる。例えば、この実験に用いた造影剤の場合は、周波数に関わらず５μ秒を超えると、（２．２２６ＭＨｚで１０パルス、６．７ＭＨｚで１００パルス程度）気泡化の閾値が急激に低下する閾値があることが実験的に示されている。このような造影剤ごとに、あらかじめ気泡化に必要なエネルギーを調べておいて、送波音圧を決めると、必要なパルス長が定まるので、そのパルス長をこれまでに説明してきた長パルスとして、用いればよい。

次に帯域幅の観点から、長パルスに関する規定を行なう。図８はハニング関数を包絡線に用いた場合の、波数と比帯域幅（中心周波数に対する帯域幅）の関係を表したグラフである。ここでは、ハニング関数を包絡線とするパルス波形をフーリエ変換し、周波数空間上で、比帯域幅を算出した。ここでは、ハニング関数の例で説明するが、包絡線はどのような関数でも構わなく、包絡線に応じて、フーリエ変換を行なうことで、図８に類似のグラフを求めることが出来る。このように一度、パルスの波数と帯域幅の関係を調べておくことで、周波数帯域上で、長パルスとエミッションパルス受信用の窓を設定することが出来る。ここでいう受信用の窓とは、図４で示す探触子の帯域内で、相変化パルス（中心周波数fa）と相変化パルスの高調波成分（中心周波数fb）の間に挟まれた領域のことである。

例えば受信用の窓内における、これらの相変化パルス及びその高調波成分それぞれからのノイズ成分をそれぞれの中心周波数での信号の大きさに対して-6dB以下と設定する場合には、それぞれの-6dB帯域幅Δfa、Δfbを用いて表記すると、受信用の窓はfa+Δfa〜fb-Δfbの範囲となる。より具体的には、例えば、探触子の-6dB帯域が、2MHzから4MHz、比帯域幅６６％のような場合で検討を行なってみる。例えば長パルスを１６波とすると、図８に示すように比帯域幅は１２％になるので、探触子帯域の下限に長パルスを設定すると、中心周波数2.12MHzと設定する場合、長パルスは、2〜2.25MHzになり、高調波は帯域幅がほぼ同じである場合には、4〜4.5MHzとなるので、エミッションパルスの受信に使える帯域は、2.25MHzから4MHzと、56%の比帯域幅（３から４波数）と比較的良好な空間分解能が確保出来る。特に、送波パルスが16波とかなり長いことに比べれば、受信パルスが4波程度（図８のグラフより）とかなり短く、かつ弁別比も良好であることがはっきりわかる。ここでは-6dB幅で議論したが、求める造影剤からの信号と組織からの信号の弁別比に応じて、より低い（例えば-20dBや、-30dBなど）感度にまで帯域を検討することも可能である。また相変化パルスからのノイズの量とその高調波からのノイズの量を別々の値に設定することも可能である。

ここまでの議論と逆に、エミッションパルスに対する受信用の窓を最適化すると、現状の超音波診断装置での空間分解能に比べて、あまり大きく劣化することは出来ないので、せいぜいエミッションパルスの波数は３波程度までが限界である。その場合比帯域幅は６０％程度であるので、探触子の比帯域幅から６０％を差し引いた比帯域幅が長パルスに用いることが出来る。例えば探触子の比帯域幅が６６％の場合、下の６％を長パルスに用い、残りの６０％を受信用の窓にあて、長パルスの高調波は帯域外に持ってくるように設計出来る。このように、受信に要求する比帯域幅から、気泡化用の送波パルスを設計することも出来る。（尚、ここで比帯域幅で足し算、引き算が出来るのは、全体の帯域ｆ１からｆ３を、ｆ１からｆ２の帯域とｆ２からｆ３の帯域に分割した場合に、分割された帯域幅それぞれが、（ｆ２−ｆ１）/（ｆ２＋ｆ１）＊２、（ｆ３−ｆ２）/（ｆ３＋ｆ２）＊２となり、その合計、（ｆ２−ｆ１）/（ｆ２＋ｆ１）＊２＋（ｆ３−ｆ２）/（ｆ３＋ｆ２）＊２＝（ｆ３−ｆ１）/（ｆ３＋ｆ１）＊２となるからである。

このような、長いパルスを送波した結果生成するエミッションパルスを用いる場合、受波のダイナミックフォーカスを何処に合わせるのか、自由度がある。パルス長が非常に長い場合、エミッションのタイミングとダイナミックフォーカスのフォーカスのタイミングがずれてしまうと、常にフォーカスがずれてしまう恐れがある。パルス長が長い場合は、パルス長の中心に対応する時間にダイナミックフォーカスを合わせることが考えられるが、造影剤の気泡化の閾値寸前のエネルギーを投与している場合、送波パルスの一番後ろのあたりに、受信ダイナミックフォーカスのタイミングを合わせておくことが、エミッションパルスの送波にフォーカスをあわせることになる。ここでは、送波パルスによって一定のエネルギーが蓄積されてから造影剤の気化が起こることを想定しており、送波パルスの長さがこの気泡化に必要なエネルギーを少し超す程度であれば、送波パルスの最後の数波のところで気泡化が起きる。ここに受信ダイナミックフォーカスをあわせるとエミッションパルスにフォーカスをあわせることになる。ただし、超音波は伝播に伴い、減衰し徐々に音圧が下がっていくので、浅いところから、深いところに向かって、気泡化の起きるタイミングが長パルスの中で後ろにずれていく傾向がある。ダイナミックフォーカスの焦点もこの点を考慮して、若干、深さ方向に応じてずらして行くことも可能である。例えば送波パルス長が10波の場合、プローブ表面では6波目、深部の限界では10波目、その中間では距離に応じて、線形補間により6波と10波の中間に来るように設定するなどの方法がある。

尚、今回の説明には、分子プローブと造影剤が結合した造影剤の例で説明を行なったが、液体を界面活性剤により微粒子化した化合物を造影剤として用いることもできる。この場合、その大きさによって、腫瘍選択性があるEPR効果を有するものであって、例えば直径１００から数１００ｎｍ程度のものを用いれば、分子プローブをつけなくても、分子イメージング用の造影剤としてつかうことが出来る。本発明で使用する造影剤は、気泡化する物質を含むものであればよい。ここで、分子プローブがついていてもよくまたついていなくてもよい。

本実施例では、造影剤からの信号のみを抽出する方法に関して説明を行なってきた。しかし、造影が、生体の構造のどの部分に集中しているかを表示するには、通常のＢモード像の上に、造影像を重ねて表示した方が、好ましい場合もある。その場合は、長パルスと、通常のＢモード撮像用の超音波の送受波を交互に繰り返して、Bモード画像と、造影画像を重畳して表示することも有用である。

なお、本発明のように、エミッションパルスに対する選択性を向上することで、気泡化した薬剤の量を推定することが容易になる。さらには、ドラッグデリバリーシステムのように、不活性化した薬剤を体内に投与し、あるタイミングで活性化して治療するような場合に、活性化した薬剤の量を推定することもできる。

第二の実施例においては、周波数帯域で区別するのではなく、組織からのエコーと造影剤気泡化の際のエミッションパルスを区別する方法を説明する。この実施例では、造影剤が気泡化しない程度の短い短パルスと、気泡化する長パルスを用いる。後者は組織エコーとエミッションパルスを含み、前者は組織エコーのみを含むため、この短パルスと長パルスの差をとると、エミッションパルスのみを抽出することが出来る。しかし短パルスと長パルスは波形が異なるため、そのままでは差をとることが出来ない。そこでこの実施例では、短パルスに対して、受信後にパルスを長くする処理を行う。短パルスを送受信が終わってから長くして、エミッションパルスを含んだ長パルスから差し引くことでエミッションパルスのみが抽出出来るわけである。シーケンスのフローチャートを図９に示す。まず、通常のＢモード撮像に用いられるものと同様の１から2程度の波数の短パルスの送受波を行なう。一つのラスターに関して送受信が終わったら、次に実施例１で説明したような長パルスの超音波の送受波を行う。この長パルス送波時の受信エコーから、組織エコーとエミッションパルスを以下のようにして区別する。通常時間軸上のコンボリューション演算で、パルスを短くすることは困難であるが、長くすることは容易である。コンボリューション演算は次のような演算である。ベクトルuとｖのコンボリューションをｗとした場合、下記の通りに現せる。

〔数１〕

本実施例では、図９のフローチャートの１番目の項目の短パルスの送波に対して、２番目の項目の通り受信したエコー信号について、コンボリューション演算によって長波長化する。例えば、短パルスがｓｉｎ波の2波で、長パルスがｓｉｎ波の１２波の場合について説明する。この場合、大体二つの波数の差の波数、すなわちｓｉｎ波の１０波の波形と、もとの短パルスの間でコンボリューションを取ると、ほぼ、ｓｉｎ１２波の波形になる。（コンボリューション演算は前述の式のように、二つの波形をずらしながら、重なりの積分を求めていく操作になるので、幅Ｔａと幅Ｔｂの波形を前後からずらしていくと、重なりが生じる区間はＴａ＋Ｔｂとなり、ｓｉｎ2波とｓｉｎ10波のコンボリューションの場合はｓｉｎ12波の波形となる。）このような長波長化処理の例を図１０に示す。まず、短パルスの送波パルスに対する受信信号の第1の受波パルスに対して、一律に長波長化の処理を行う。

次に、長パルスの送波パルスに対する受信信号の第２の受波パルスを取得する。先の長波長化処理では、第1の受波パルスの波数を、第2の受波パルスの波数と同一にする波数処理を行うものとする。

第1の受波パルスは、組織エコーに対応する信号である。一方で、第２の受波パルスは、組織エコーとエミッションパルスの重ね合わせに対応する信号である。そこで、第２の受波パルスから第１の受波パルスを差し引くことで、エミッションパルスのみを抽出することが出来る。

装置のブロック図で説明すると、図１２のようになる。図１と共通の部分は説明を省略する。まずパルス照射制御部２５によって、長パルスと短パルスの切り替えの制御を、送波波形メモリと長波長化処理部に対して行う。受信ビームフォーマの出力を長波長化処理部２６によって、前述のような長波長化処理を行う。この長波長化処理を通った短パルスの受信信号と、長パルスの受信信号の差を計算部２７において、行い、エミッションパルスの抽出が行われる。

なお、長波長化処理は、より厳密には、長パルスを短パルスでデコンボリューションした関数でも良いし、次に示すような最小自乗法を用いて、フィルタ設計しても良い。なお、ベクトルAをベクトルBでデコンボリューションした結果がベクトルCとすると、BとCのコンボリューションがAとなるようなベクトルCを見つけることが、デコンボリューションである。

以下，不整合フィルタを用いた長パルス化用フィルタの設計方法について説明する。短パルス波形をＢ，長パルス化フィルタをｆ，長パルス化後の信号をｃとすると，長パルス化後の信号ｃは，（数２）の式で与えられる。以下の説明で，記号Ｂ，Ｗは行列を示しｃ，ｆ，ｄ，ｗはベクトル、記号「Ｔ」は転置を示す。

〔数２〕

長パルス化の信号Ｃとして望ましい波形をＤとすると，ＣとＤの自乗誤差和Ｉは（数３）となる。自乗誤差和Ｉを最小にするＦが不整合フィルタである。

〔数３〕
Ｉ＝Σ（ｃ_ｉ−ｄ_ｉ）^２＝（ｆＢ−ｄ）（ｆＢ−ｄ）^Ｔ
＝ｆＢＢ^Ｔｆ^Ｔ−ｄＢ^Ｔｆ^Ｔ−ｆＢｄ^Ｔ＋ｄｄ^Ｔ
（数４）の条件から（数５）を全てのｉ（ｉ＝１，２，…，ｍ）について求めると（数６）となるので，ｆは（数７）のように求まる。

〔数４〕
∂Ｉ／∂ｆ_ｉ＝０

〔数５〕
∂Ｉ／∂ｆ_ｉ＝Ｂ_ｉＢ^Ｔｆ^Ｔ＋ｆＢＢ_ｉ ^Ｔ−ｄＢ_ｉ ^Ｔ−Ｂ_ｉｄ^Ｔ
＝２（ｆＢＢ_ｉ ^Ｔ−ｄＢ_ｉ ^Ｔ）＝０

〔数６〕

〔数７〕
ｆ＝ｄＢ^Ｔ（ＢＢ^Ｔ）^−１
このように求めた不整合フィルタの概念によった長波長化フィルタを用いることで、より厳密にエミッションパルス成分を抽出することが出来る。

ここで、なるべく気泡が存在しない状態から、気泡化を行った方が、エミッションパルスを選択的に造影しやすい。そのため、ラスターの順番は端から順に一方方向に走査するよりは、なるべく時間的に隣接する送受波のラスターが空間的に離れるように設定した方が好ましい場合もある。例えば、ラスター数が100本で、左から、順にラスター１，２，３、、、１００と番号がついている場合に、１、５１、２５、７６、１３、３８、６３，８９、２、、、のような順番などが考えられる。このようにすると、隣接ラスター間の時間間隔を、順にラスター走査する場合に比べ、８倍に長くすることが出来る。

本発明の装置例のブロック図。 従来の考え方にたつ気泡化された造影剤のイメージング法の説明図。 本発明の気泡化と造影画像の取得を同時に行う方法の説明図。 本発明での送波パルスとエミッションパルスの時間軸波形と周波数帯域の説明図。 造影剤エミッションパルス信号の抽出シーケンスのフローチャート。 従来方法における組織エコーと造影信号の弁別比と空間分解能のトレードオフを説明する図。 パルス波長に対する気泡化閾値の実験データ。 ハニング関数を包絡線とした場合のパルスの波数と比帯域幅の関係を示すグラフ。 第二の実施例に関するフローチャート。 第二の実施例に関する波形を説明する図。 探触子帯域を説明する図。 本発明の装置例のブロック図。

符号の説明

１…超音波探触子、２…送受切替スイッチ、３…送波ビームフォーマ、４…制御系、５…送波波形メモリ、２０…受波ビームフォーマ、２１…造影剤エミッションパルス抽出部、２３…スキャンコンバータ、２４…表示部、２５…パルス照射制御部、２６…長波長化処理部、２７…計算部。

Claims (9)

1. 相変化造影剤投与領域である被検体の対象領域へ、前記対象領域に投与される造影剤の相変化に必要なパルス波数を有する第１超音波信号と、前記第１超音波信号よりも波数の少ない第３超音波信号とを送信する送信手段と、
   前記被検体から前記第１超音波信号の照射によって生じる第２超音波信号と、前記第３超音波信号によって生じる第４超音波信号とを受信する受信手段と、
   前記第４超音波信号について、波数を前記第２超音波信号と同一にする波数処理を行って、第５超音波信号とし、第２超音波信号と第５超音波信号との差分を、前記第２超音波信号から前記第１超音波信号の基本波パルスと高調波パルスを除去した、エミッションパルスとして算出する計算部とを有する超音波撮像装置。
2. 前記計算部はバンドパスフィルタを有することを特徴とする請求項１記載の超音波撮像装置。
3. 前記計算部の計算結果に基づいて画像を表示する表示部をさらに有し、前記表示部は、
   Ｂモード像に前記計算部で算出した検出信号に基づく造影像を重ねて表示することを特徴とする請求項１記載の超音波撮像装置。
4. 前記送信手段は前記第1超音波信号とＢモード像を取得するための超音波信号とを交互に送信することを特徴とする請求項１記載の超音波撮像装置。
5. 前記計算部は波長処理部を有し、前記波長処理部は、前記第４の超音波信号の波数と、
   前記第４の超音波信号の波数と前記第１の超音波信号との波数の差の波数を有するパルスとのコンボリューション演算を行うことを特徴とする請求項１記載の超音波撮像装置。
6. 前記計算部は波長処理部を有し、前記波長処理部は、前記第１の超音波信号を前記第３の超音波信号でデコンボリューションした関数を用いて処理することを特徴とする請求項１記載の超音波撮像装置。
7. 前記計算部は波長処理部を有し、前記波長処理部は、前記第４の超音波信号とコンボリューションをとった結果前記第１の超音波信号との自乗誤差を最小とする関数を用いて処理することを特徴とする請求項１記載の超音波撮像装置。
8. 前記計算部は波長処理部を有し、前記波長処理部は、不整合フィルタを有することを特徴とする請求項１に記載の超音波撮像装置。
9. 前記第３の超音波信号は組織エコーに該当する信号であり、前記第２の超音波信号は前記組織エコーとエミッションパルスとの重ね合わせに該当する信号であることを特徴とする請求項１に記載の超音波撮像装置。

Images