(57)【要約】　　　（修正有）【課題】位相回転効果を補正し、真の生体内音響インピーダンス変化のイメージングを可能とする超音波撮像装置及び超音波撮像方法を提供する。【解決手段】検査対象から探触子１に戻った超音波信号は探触子１によって電気信号に変換され、送受切り替えスイッチ群５を介し複素受波ビームフォーマ２０に送られ、受信のタイミングに応じて遅延時間を調整するダイナミックフォーカスが実行される。実部と虚部のビームを出力する複素ビームフォーマ２０の出力を用いて、位相回転補正部２１で、周波数に依存する減衰に起因する位相回転の補正、または、ビームの方位方向の位相回転分の補正、もしくはこれらの両方の補正を行なう。位相回転の補正の語、音響インピーダンス変化量演算部２２で、音響インピーダンスの空間位置に関する微分を求めフィルタ処理がなされた信号は、スキャンコンバータ２３を経て表示部２４に画像表示される。【選択図】　図１
2004113364.tif 000002
【特許請求の範囲】
【請求項１】
超音波探触子を用いて検査対象に超音波を送受波して、前記検査対象を画像化する超音波撮像装置において、前記検査対象に送波信号を送信し、前記検査対象から返ってきた受波信号の前記送波信号を基準とした位相変化から、前記検査対象内の音響インピーダンスの分布に起因する前記受波信号の位相のずれと、それ以外の要因に起因する前記受波信号の位相ずれを区別して検出する手段を有することを特徴とする超音波撮像装置。
【請求項２】
請求項１に記載の超音波撮像装置において、前記受波信号に対して前記検出された前記位相ずれを前記受波信号に対して補正する手段を有することを特徴とする超音波撮像装置。
【請求項３】
請求項１に記載の超音波撮像装置において、前記音響インピーダンスの分布に起因する位相情報を表示する手段を有することを特徴とする超音波撮像装置。
【請求項４】
請求項１に記載の超音波撮像装置において、前記受波信号を直交検波して複素信号で表現する手段と、前記複素信号を２乗する手段と、ローパスフィルタと、周波数に依存する減衰に起因する位相回転を前記ローパスフィルタにかけた前記複素信号に基づいて補正する手段と、前記位相回転が補正された前記複素信号から前記検査対象の前記音響インピーダンスの分布もしくは前記音響インピーダンスを求める手段とを有することを特徴とする超音波撮像装置。
【請求項５】
請求項１に記載の超音波撮像装置において、基本波の偶数次の高調波を重畳した送波信号を送信する手段と、受信された前記第２高調波の直交成分を利用することにより、前記受波信号に対して前記超音波の伝播の過程で生じる周波数に依存する減衰による位相回転による位相ずれを前記受波信号に対して補正する手段と、前記位相回転が補正された複素信号から前記検査対象の音響インピーダンスを求める手段とを有することを特徴とする超音波撮像装置。
【請求項６】
請求項１に記載の超音波撮像装置において、前記超音波の送受信の条件によって定まる点応答関数の方位方向の位相回転を算出する手段と、前記位相回転が最小となるフィルタリング処理手段と、回折効果による位相回転に起因して生じる前記超音波のビームの方位方向の位相ずれを前記受波信号に対して補正する手段と、前記位相ずれが補正された複素信号から前記検査対象の音響インピーダンスを求める手段を有することを特徴とする超音波撮像装置。
【請求項７】
請求項６に記載の超音波撮像装置において、前記点応答関数の方位方向の位相回転を最小とするために最小２乗法に基づいたフィルタリング処理手段と、前記方位方向の位相回転分の影響を補正して前記検査対象内での前記音響インピーダンスの分布を抽出する手段と、前記音響インピーダンスの分布を画像化する手段を有することを特徴とする超音波撮像装置。
【請求項８】
超音波探触子を用いて検査対象に超音波を送受波して、前記検査対象を画像化する超音波撮像装置において、前記検査対象に送波信号を送信し、前記検査対象から返ってきた受波信号の前記送波信号を基準とした位相変化から、周波数に依存する減衰による位相回転、及び／又は、回折効果による位相回転、に起因する位相のずれを補正する手段と、前記検査対象中の共振物体による位相変化を抽出もしくは強調する手段を有することを特徴とする超音波撮像装置。
【請求項９】
請求項８に記載の超音波撮像装置において、前記超音波の送信周波数を変化させることにより、共振状態のスペクトロスコピーを反映したイメージングを行なう手段を有することを特徴とする超音波撮像装置。
【請求項１０】
請求項８に記載の超音波撮像装置において、前記位相変化を抽出もしくは強調する手段は、前記検査対象中の超音波造影剤による位相変化を抽出もしくは強調することを特徴とする超音波撮像装置。
【請求項１１】
請求項１０に記載の超音波撮像装置において、前記受波信号を直交検波して複素信号で表現する手段と、前記複素信号を２乗する手段と、ローパスフィルタと、周波数に依存する減衰に起因する位相回転を前記ローパスフィルタにかけた前記複素信号に基づいて補正する手段と、前記位相回転の補正の後に残る前記複素信号の位相変化部から前記検査対象の中の超音波造影剤の位置、量、移動速度を求める手段とを有することを特徴とする超音波撮像装置。
【請求項１２】
請求項１０に記載の超音波撮像装置において、基本波に対しその第２高調波を重畳した送波信号を送信することにより、前記送波信号の位相に対する同相成分の受信された前記第２高調波を利用することにより前記超音波の伝播の過程で生じる位相回転と造影剤の存在による位相変化を分離する手段と、前記検査対象の超音波造影剤の位置、量、移動速度を求める手段とを有することを特徴とする超音波撮像装置。
【請求項１３】
請求項１０に記載の超音波撮像装置において、送受信の条件によって定まる点応答関数の方位方向の位相回転を算出し、前記位相回転が最小となるフィルタリング処理によって、回折効果による位相回転に起因して生じる前記超音波のビームの方位方向の位相ずれの補正を行い、前記超音波のビームの回折で生じる位相回転と造影剤の存在による位相変化を分離する手段と、前記検査対象の超音波造影剤の位置、量、移動速度を求める手段とを有し、前記点応答関数の方位方向の位相回転を最小とするために最小２乗法に基づいたフィルタリング処理を行なうことを特徴とする超音波撮像装置。
【請求項１４】
超音波探触子を用いて検査対象に超音波を送受波して、前記検査対象を画像化する超音波撮像装置において、前記検査対象に送波信号を送信し、前記検査対象から返ってきた受波信号の前記送波信号を基準とした位相変化から、前記検査対象の音響インピーダンスの不連続性に起因する前記受波信号の位相のずれと、それ以外の要因に起因する前記受波信号の位相ずれを区別する手段と、前記検査対象の前記音響インピーダンスの時間変化を求め、前記音響インピーダンスの時間変化を表示する手段と有することを特徴とする超音波撮像装置。
【請求項１５】
請求項１４に記載の超音波撮像装置と、前記超音波撮像装置により得られた前記音響インピーダンスの不連続性の時間変化を、超音波を用いた治療装置による治療状態として、前記治療装置にフィードバックする手段を有することを特徴とする治療装置。
【請求項１６】
超音波探触子を用いて検査対象に超音波を送受波して、前記検査対象を画像化する超音波撮像装置において、前記検査対象に送波信号を送信し、前記検査対象から返ってきた受波信号の前記送波信号を基準とした位相変化から、周波数に依存する減衰による位相回転、及び／又は、回折効果による位相回転、に起因して生じる前記超音波のビームの方位方向の位相ずれを前記受波信号に対して補正する手段とを有することを特徴とする超音波撮像装置。
【請求項１７】
超音波探触子を用いて検査対象に超音波を送受波して、前記検査対象を画像化する超音波撮像装置において、前記検査対象に送波信号を送信する手段と、前記検査対象から返ってきた受波信号を直交検波する手段と、直交検波された前記受波信号から位相変化に対応する成分を求める手段と、前記位相変化に対応する成分から急峻な変化を除去するローパスフィルタと、前記ローパスフィルタの出力信号を用いて、前記検査対象から返ってきた受波信号の前記送波信号を基準とした位相変化から、周波数に依存する減衰による位相回転、及び／又は、回折効果による位相回転、に起因して生じる前記超音波のビームの方位方向の位相ずれを前記受波信号に対して補正する手段と、表示する位相を選択する入力手段と、前記入力手段によって選ばれた位相、もしくは位相の範囲の信号を抽出もしくは強調して表示する手段とを有することを特徴とする超音波撮像装置。
【請求項１８】
超音波探触子を用いて検査対象に超音波を送受波して、前記検査対象を画像化する超音波撮像装置において、前記検査対象に基本波の偶数次の高調波を重畳した送波信号を送信する手段と、前記検査対象から返ってきた受波信号を直交検波する手段と、直交検波された前記受波信号から前記偶数次の高調波に対する特定の位相成分を入力する手段と、前記入力手段によって選ばれた範囲の位相の信号を抽出もしくは強調する手段とこれを表示する手段とを有することを特徴とする超音波撮像装置。
【請求項１９】
超音波探触子を用いて検査対象に超音波を送受波して、前記検査対象を画像化する超音波撮像方法において、前記検査対象に送波信号を送信する工程と、前記検査対象から返ってきた受波信号の前記送波信号を基準とした位相変化から、周波数に依存する減衰による位相回転、及び／又は、回折効果による位相回転に起因して生じる前記超音波のビームの方位方向の位相ずれを補正する工程と、補正された前記受波信号から前記検査対象の音響インピーダンス画像、及び／又は、前記音響インピーダンスの空間位置に関する微分画像を求める工程と、前記音響インピーダンス画像、及び／又は、前記微分画像を表示する工程とを有することを有することを特徴とする超音波撮像方法。
【請求項２０】
超音波探触子を用いて検査対象に超音波を送受波して、前記検査対象を画像化する超音波撮像方法において、前記検査対象に送波信号を送信し、前記検査対象から返ってきた受波信号を直交検波する工程と、直交検波された前記受波信号から位相変化に対応する成分を求める工程と、前記位相変化に対応する成分から急峻な変化をローパスフィルタにより除去する工程と、前記ローパスフィルタの出力信号を用いて、前記検査対象から返ってきた受波信号の前記送波信号を基準とした位相変化から、周波数に依存する減衰による位相回転、及び／又は、回折効果による位相回転、に起因して生じる前記超音波のビームの方位方向の位相ずれを前記受波信号に対して補正して、位相が補正された前記受信信号に基づいて、前記検査対象の音響インピーダンス画像、及び／又は、前記音響インピーダンスの空間位置に関する微分画像を求める工程と、前記音響インピーダンス画像、及び／又は、前記微分画像を表示する工程とを有することを特徴とする超音波撮像方法。
【請求項２１】
請求項２０に記載の超音波撮像方法において、前記音響インピーダンスが前記検査対象に注入された造影剤に基づくことを特徴とする超音波撮像方法。
【請求項２２】
超音波探触子を用いて検査対象に超音波を送受波して、前記検査対象を画像化する超音波撮像方法において、前記検査対象に基本波の偶数次の高調波を重畳した送波信号を送信する工程と、前記検査対象から返ってきた受波信号を直交検波し、直交検波された前記受波信号から前記偶数次の高調波の直交成分を抽出する工程と、抽出された前記直交成分に基づいて、前記検査対象の音響インピーダンス画像、及び／又は、前記微分画像を表示する工程と、前記音響インピーダンス画像、及び／又は、前記微分画像を表示する手段とを有することを特徴とする超音波撮像方法。
【請求項２３】
請求項２２に記載の超音波撮像方法において、前記音響インピーダンスが前記検査対象に注入された造影剤に基づくことを特徴とする超音波撮像方法。
【請求項２４】
超音波探触子を用いて検査対象に超音波を送受波して、前記検査対象を画像化する超音波撮像方法において、前記検査対象に基本の偶数次の高調波を重畳した送波信号を送信する工程と、前記検査対象から返ってきた受波信号を直交検波し、直交検波された前記受波信号から前記偶数次の高調波の同相成分を抽出する工程と、抽出された前記同相成分に基づいて、前記検査対象の音響インピーダンス画像、及び／又は、前記微分画像を求める工程と、前記音響インピーダンス画像、及び／又は、前記微分画像を表示する工程とを有することを特徴とする超音波撮像方法。
【請求項２５】
請求項２４に記載の超音波撮像方法において、前記音響インピーダンスが前記検査対象に注入された造影剤に基づくことを特徴とする超音波撮像方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、超音波を用いて検査対象内の画像を撮像する超音波撮像装置及び超音波撮像方法に関する。
【０００２】
【従来の技術】
医療画像診断に用いられる超音波撮像装置は、超音波パルス反射法を用いて、生体の軟部組織の断層像や生体内を流れる血流像等をほぼリアルタイムでモニタに表示して観察でき、また、放射線を用いる画像診断装置のような放射線被爆を検査対象に与えないことから安全性も高いとされ、更に小型で安価なことも加わり、広く医療の分野で用いられている。
【０００３】
従来の超音波撮像装置の撮像メカニズムは、以下のように説明される。超音波撮像装置においては、圧電体をアレイ化した超音波探触子から超音波を検査対象に送波する。アレイの各素子から超音波を送波するタイミングをコントロールすることで、超音波が検査対象内で焦点を結ぶ位置をコントロールし、アレイの各素子が受波した信号を所望の位置と各素子位置の距離の差に応じた時間をずらして加算することで、所望の位置からの信号を増強する。超音波パルスを撮像領域全域に渡り走査することで、断層像が得られる。従来、このような超音波のパルスエコー法は、検査対象中の音響インピーダンスが変化する界面からのエコー（受信信号）を受信し、受信信号を包絡線検波することで、そのエコーの絶対値を断層像の画像化に用いてきた。
【０００４】
一方、エコーの符号を残したまま、もしくは位相情報を用いて画像化を行う試みに関する従来技術が提案されている（例えば、特許文献１、特許文献２参照。）。
【０００５】
【特許文献１】
特開昭５５−１３６０４３号公報
【特許文献２】
特開平１１−１３７５４６号公報
【０００６】
【発明が解決しようとする課題】
しかし、前記の従来技術の方法で、エコー（受信信号）の符号情報を利用しようとすると、音響インピーダンスが変化したことによる反射信号の位相変化と、その他の要因による反射信号の位相変化が区別出来ないため、実際の検査対象の音響インピーダンス分布を画像化することは困難であった。
【０００７】
上述した、その他の要因とは、以下の（１）、（２）の２つの要因による。
【０００８】
まず、（１）生体中の減衰によって、超音波パルスの中心周波数が変化するため、見かけ位相が回っている。もう１つは、（２）超音波超音波ビームの方位方向への位相まわりのために、注目する超音波ビーム近傍の反射体からの信号が位相が回った信号として受信される。
【０００９】
このような従来技術の方法における第１の課題は、上記（１）、（２）による受信信号の位相回転効果を補正し、真の生体内音響インピーダンス変化のイメージングとを可能とすることである。
【００１０】
また、第２の課題としては、以下に述べる加熱凝固治療の治療効果のモニタリングに関する課題である。加熱凝固治療とは、強力収束超音波照射もしくはＲＦ波やマイクロ波などの電磁波の照射によって、熱的に前立腺肥大、前立腺がんや肝臓がんなどの患部を治療する方法である。前述の位相情報の画像化は、近年発達して来た、加熱凝固を用いた低侵襲治療の発達を益々加速するものとして期待される。
【００１１】
これらの治療法は、体表の大きな切開を伴わずに治療出来るので、臨床での使用が拡大してきているが、この場合に効果判定方法が重要と成る。がんの治療もれがあると、結局そこから再びがんが増殖してしまうため、治療目的部位を完全に治療しきったかどうか判定することが重要である。現在は、造影ＣＴ、造影ＭＲＩ、造影超音波などの方法を用いているが、ＣＴやＭＲＩなど治療中に撮像することが困難であったり、造影剤を用いているため、一度治療をやめて効果判定をして、だめであれば再び治療をするという処理になっている。
【００１２】
これに対し、より効果的な治療を行なうために、治療しながら効果判定が出来ることが望まれていた。加圧法を用いた歪イメージングなどが超音波で造影剤を用いない、組織性状の画像化方法として検討されているが、これも、加圧が効果的に行なえる場所にしか適用出来ないと考えられる。上記の理由から簡便な組織の硬さのイメージング法が必要とされていた。
【００１３】
さらに、第３の課題としては、造影超音波に関する課題である。造影超音波は造影剤を注入するという行為をトリガとして、生体の機能の時間差を観測する方法である。これによって、例えば、肝臓の組織毎の血管系の違い（即ち、肝臓がんの場合は動脈支配であるのに対し、それ以外の正常組織部分は門脈支配）のコントラストが画像化可能となる。
【００１４】
この場合、造影剤からの信号と、生体組織からの信号を区別して画像化することが重要となる。従来は、造影剤の応答の非線形性を用いて、第２高調波を使う等反射信号の周波数変化を使う方法によって造影剤と生体組織からのエコー信号のコントラストをつけてきたが、生体組織も非線形性を持つため、このコントラストがつきにくく、このコントラストをさらに改善したいという要求があった。
【００１５】
本発明は、上記の点に鑑みてなされたものであり、位相回転効果を補正し、真の生体内音響インピーダンス変化のイメージングを可能とする超音波撮像装置及び超音波撮像方法を提供することを目的とする。
【００１６】
【課題を解決するための手段】
従来、音響インピーダンス撮像において、検査対象内の音響インピーダンス変化に起因するエコーの位相回転と、超音波ビームの方位方向の位相回転及び周波数シフトに起因する位相回転が混ざり真の検査対象内の音響インピーダンス変化を撮像することが困難であったが、本発明の超音波撮像装置では、超音波ビームの方位方向の位相回転及び周波数シフトに起因する位相回転を仮想的な第２高調波法又は偶数次の高調波法と、マトリックスフィルタによって補正することで、検査対象内の音響インピーダンス変化のみを抽出してイメージングすることが可能となる。
【００１７】
本発明の超音波撮像装置では、対象に超音波を送信し、該送信によって該対象から返ってきた受信波形について、送信波形を基準とした該受信波形の位相変化から、観測対象の音響インピーダンスの空間変化（あるいは、分布）に起因するエコーの位相のずれと、それ以外の要因に起因するエコーの位相ずれを区別する手段を備え、該対象の音響インピーダンスの空間変化を推定する手段を備えたことを特徴とする。なお、上述した、それ以外の要因とは、周波数依存減衰、及び／もしくは、回折効果による位相回転を扱うことが望ましい。
【００１８】
周波数依存減衰に起因する位相回転の除去方法の一つの手段としては、前記受信波形を直交検波して複素信号表現する手段を備え、これを二乗する手段を備え、音響インピーダンスの空間変化に起因する位相回転以外の効果を抽出する手段を備え、逆に周波数依存減衰に起因する位相回転を推定する手段を備え、これを除去する手段を備えることである。
【００１９】
周波数依存減衰に起因する位相回転の除去方法の別の一つの手段としては、前記送信波形として、基本波に対しその第二高調波もしくは偶数次高調波を重畳した波形を送波する手段を備え、受信第二高調波の直交成分を利用することにより伝播の過程で生じる位相補正する手段を備えることである。これらいずれかの方法によって、周波数依存減衰に起因する位相回転の効果を補正する手段を備えることが可能となる。
【００２０】
回折効果による位相回転に起因する位相のずれに対しては、送受信の条件によって定まる点応答関数の方位方向の位相回転を算出する手段を備え、この位相回転が最小となるフィルタリング処理手段を備えることで、除去することが可能となる。
【００２１】
本発明によれば、反射信号の位相回転を、観測対象の空間分布に起因するものと、伝播の過程で蓄積する位相回転や、回折効果に起因するものを分離し、反射信号のなかから特定の位相成分を抽出する手段を備えることで、観測対象の音響インピーダンスの変化や、音響インピーダンスの微分の空間変化や、もしくは共鳴状態、造影剤の分布などをイメージングする手段を備えることが可能となる。
【００２２】
また、本発明は、治療に用いることも可能である。超音波により観測対象の状態を画像化する超音波モニタリング機構を備えた治療装置において、対象に超音波を送信し、該送信によって該対象から返ってきた受信波形について、前記の周波数依存減衰に起因する位相回転、及び／もしくは、回折効果を補正して送信波形を基準とした該受信波形の位相変化から、観測対象の音響インピーダンスの時間変化を推定する手段を備え、表示する手段を備え、この結果を治療状態にフィードバックする手段を備えることが可能となる。
【００２３】
本発明の超音波撮像方法では、対象に超音波を送信し、該送信によって該対象から返ってきた受信波形について、送信波形を基準とした該受信波形の位相変化から、観測対象の音響インピーダンスの空間変化に起因するエコーの位相のずれと、それ以外の要因に起因するエコーの位相ずれを区別する方法であり、該対象の音響インピーダンスの空間変化（あるいは、分布）を推定することを特徴とする。なお、上述した、それ以外の要因とは、周波数依存減衰、及び／もしくは、回折効果による位相回転を扱うことが望ましい。
【００２４】
周波数依存減衰に起因する位相回転の除去方法の一つとしては、前記受信波形を直交検波して複素信号表現し、これを二乗し、音響インピーダンスの空間変化に起因する位相回転以外の効果を抽出し、逆に周波数依存減衰に起因する位相回転を推定し、これを除去する。
【００２５】
周波数依存減衰に起因する位相回転の除去方法の別の一つ方法としては、前記送信波形として、基本波に対しその第二高調波もしくは偶数次高調波を重畳した波形を送波し、受信第二高調波の直交成分を利用することにより伝播の過程で生じる位相補正することである。これらいずれかの方法によって、周波数依存減衰に起因する位相回転の効果を補正することが可能となる。
【００２６】
回折効果による位相回転に起因する位相のずれに対しては、送受信の条件によって定まる点応答関数の方位方向の位相回転を算出し、この位相回転が最小となるフィルタリング処理を行ない、回折効果に起因する位相回転の効果を除去することが可能となる。
【００２７】
本発明によれば、反射信号の位相回転を、観測対象の空間分布に起因するものと、伝播の過程で蓄積する位相回転や、回折効果に起因するものを分離し、反射信号のなかから特定の位相成分を抽出することで、観測対象の音響インピーダンスの変化や、音響インピーダンスの微分の空間変化や、もしくは共鳴状態、造影剤の分布などをイメージングすることが可能となる。
【００２８】
【発明の実施の形態】
以下、本発明の実施例を、図面を用いて詳細に説明する。
【００２９】
（実施例１）
図１は、本発明の実施例１において音響インピーダンスの微分画像を得る超音波撮像装置の構成例を示す図である。
【００３０】
図２は、本発明の実施例１において音響インピーダンス画像を得る超音波撮像装置の構成例を示す図である。
【００３１】
超音波探触子１は、検査対象（図示してない）に対し、超音波ビームを送信し、検査対象からのエコーを受信する。送波信号は、送波ビームフォーマ３により、制御系４の制御下で、信号に応じた送波焦点に合わせた遅延時間で、送受切り替えＳＷ（スイッチ）群５を介して超音波探触子１に送られる。送受切り替えＳＷ群５は制御系４の制御下にある。検査対象内で反射もしくは散乱されて超音波探触子１に戻った超音波信号は、超音波探触子１によって電気信号に変換され、再び送受切り替えスイッチ群５を介し複素受波ビームフォーマ２０に送られる。ここで、複素受波ビームフォーマ２０は、制御系４の制御の下で、受信のタイミングに応じて遅延時間を調整するダイナミックフォーカスを行う。複素受波ビームフォーマ２０で、９０度位相がずれた２つの信号のミキシングを行っており、実部と虚部のビームを出力する。複素受波ビームフォーマ２０の出力信号を用いて、位相回転補正部２１で、周波数に依存する減衰に起因する位相回転の補正、または、ビームの方位方向の位相回転分の補正、もしくはこれらの両方の補正を行なう。
【００３２】
位相回転補正部２１での位相回転の補正の後、図２に示す構成では、音響インピーダンス演算部２５で音響インピーダンスを計算する、また、図１に示す構成では、音響インピーダンス変化量演算部２２で、音響インピーダンスの空間位置に関する微分を計算する。音響インピーダンス変化量演算部２２（図１）、音響インピーダンス演算部２５（図２）の出力信号に対して、従来公知である図示しないフィルタ処理がなされた信号は、スキャンコンバータ２３を経て表示部２４に結果として画像表示される。
【００３３】
この画像表示は、従来のＢモード画像、図１で得られる音響インピーダンスの微分画像、図２で得られる音響インピーダンス画像の何れか、もしくはこれらから選択された複数の画像が重畳して表示可能なように構成される。音響インピーダンス画像の表示に関しては、平均的な硬さのところに対してより硬いところは赤、柔らかいところは青などのように硬さをカラーで表示する方法などが考えられる。
【００３４】
次に、生体中からのエコーの位相について詳細に説明する。
【００３５】
図３は、本発明の実施例１においてエコーの中心周波数の移動を説明する図である。
【００３６】
図４は、本発明の実施例１においてエコーの位相に対する中心周波数の移動の影響を説明する図である。
【００３７】
図４（ａ）は、音響インピーダンスの変化を模試的に示す図、図４（ｂ）は、位相回転が無い場合の、図４（ａ）に示す音響インピーダンスの変化に対する受波ＲＦ波形（エコー）を示す。
【００３８】
図４（ｂ）に示すように、エコーの位相は、硬い（音響インピーダンズが大きい）ものに変る界面では変化しないが、柔らか（音響インピーダンズが小さい）いものに変る界面では１８０度回転する。しかし、実際には、図３（ａ）に示すように、生体減衰は周波数が高くなるほど大きくなるので、より長い距離を伝播して来た生体深部からの反射信号ほど、高周波成分が失われる。よって、図３（ｂ）（送信帯域、探触子に対する３通りのエコーの伝播長に対するエコーの中心周波数の変化を示す。）に示すように、反射信号の中心周波数はエコーの反射してくる深さによって低周波数側へ変化してしまう。
【００３９】
このため、図４（ｃ）に示すように、生体減衰に起因する位相回転のために受波ＲＦ波形は位相回転を受ける。反射信号の生体による減衰に起因する中心周波数の変化は見せかけの位相回転を与えるので、このままでは硬さの変化に起因する位相回転なのか、中心周波数移動に起因するものか区別する事が出来ない。もし生体が完全に均質な媒体であれば、中心周波数の移動は補正が可能であるが、実際には音速も減衰係数も場所によって異なるため、これを補正するのは容易ではない。
【００４０】
そこで、本発明では、まず、この中心周波数移動を第２高調波を考えることで補正する。図４（ｄ）は、位相回転が無い場合の、図４（ａ）に示す音響インピーダンスの変化に対する第２高調波による受波ＲＦ波形を示す。図４（ｅ）は、位相回転がある場合の、図４（ａ）に示す音響インピーダンスの変化に対する第２高調波による受波ＲＦ波形を示す。図４（ｄ）と図４（ｅ）の比較から明らかなように、図４（ｅ）に示す受波ＲＦ波形は、硬さに起因する位相回転の影響を受けない。即ち、基本波にとっての１８０度の位相回転は、第２高調波にとっては、３６０度の位相回転になるので、０度の位相回転と区別されないので、中心周波数移動による位相回転の効果のみを見ることが出来る。
【００４１】
図５は、本発明の実施例１において仮想的な高調波を想定して音響インピーダンス像を求める信号処理プロセス（仮想高調波位相法信号処理プロセス）を説明する図である。
【００４２】
図５に示す処理は、第２高調波を用いた中心周波数の移動補正の例である。これは、受信後に仮想的な第２高調波を考えた処理によって生体減衰に起因する受信信号の中心周波数の変化によって引き起こされる位相変化をトラッキングする方法である。
【００４３】
複素ビームフォーマ２０による直交検波処理１００により、受信信号はｃ（ｔ）ｃｏｓ（ω０ｔ）＋ｓ（ｔ）ｓｉｎ（ω０ｔ）と展開される。
【００４４】
次に、位相成分演算処理１０１により位相変化に対する成分を求める。ａ（ｔ）＝ｓｑｒｔ（ｃ２（ｔ）＋ｓ２（ｔ））を求め、Ｃ（ｔ）＝ｃ（ｔ）／ａ（ｔ）、Ｓ（ｔ）＝ｓ（ｔ）／ａ（ｔ）を求める。
【００４５】
次に、仮想２乗成分演算処理１０２により、検波信号の実部Ｃ（ｔ）と虚部Ｓ（ｔ）に対する複素信号Ｃ（ｔ）＋ｊＳ（ｔ）の２乗信号Ｃ２（ｔ）＋ｊＳ２（ｔ）を求める。ここで、ｊは虚数単位である。Ｃ２（ｔ）はＣ２（ｔ）＝Ｃ２（ｔ）−Ｓ２（ｔ）から、Ｓ２（ｔ）はＳ２（ｔ）＝２Ｃ（ｔ）Ｓ（ｔ）から導かれる。
【００４６】
この複素２乗信号Ｃ２（ｔ）＋ｊＳ２（ｔ）に対し、ローパスフィルタ処理１０３でローパスフィルタ（ＬＰＦ）を通し、急峻な変化を取り除くことが好ましい。ローパスフィルタ処理１０３により、ローパスフィルタの出力信号として、Ｃ３（ｔ）＋ｊＳ３（ｔ）＝ＬＰＦ〈Ｃ２（ｔ）＋ｊＳ２（ｔ）〉を求める。
【００４７】
複素共役演算処理１０４により、Ｃ３（ｔ）−ｊＳ３（ｔ）を求める（省略可能）。
【００４８】
次に、開平処理１０５により、Ｃ４（ｔ）＝±ｓｑｒｔ（（１＋Ｃ３（ｔ））／２）、Ｓ４（ｔ）＝ｓｇｎ（−Ｓ３Ｃ４）ｓｑｒｔ（（１−Ｃ３（ｔ））／２）を求め、Ｃ４（ｔ）、Ｓ４（ｔ）を用いて、エコー位相回転補正処理１０６により、Ｃ５（ｔ）＋ｊＳ５（ｔ）＝（Ｃ（ｔ）＋ｊＳ（ｔ））（Ｃ４（ｔ）＋ｊＳ４（ｔ））を求め、エコーの位相を補正する。
【００４９】
符号の決定処理１０７により、例えば、レンズと生体界面での反射信号など、符号の既知なエコー信号に着目し、Ｃ４の符号を決定する。音響インピーダンス決定処理１０８により、Ｃ５（ｔ）をもとに音響インピーダンス像を求める。
【００５０】
超音波の生体の伝播中の減衰に起因する位相回転補正に関しては、以下のような方法も有力である。
【００５１】
図６は、本発明の実施例１において高調波を重畳して送信する超音波撮像装置の構成例を示す図である。
【００５２】
図７は、本発明の実施例１において高調波を重畳して送信して音響インピーダンス像を求める信号処理プロセス（実高調波位相法によるインピーダンス撮像信号処理プロセス）を説明する図である。
【００５３】
図６に示す超音波撮像装置の構成は、図１に示す超音波撮像装置の構成に送波波形（第２高調波重畳波形）メモリ２６を付加した構成である。図６に示す超音波撮像装置の構成において、音響インピーダンス変化量演算部２２を音響インピーダンス演算部２５に代えることもできる。図６に示す超音波撮像装置は、実際に高調波を重畳して送信する構成である。この場合、送信信号は基本波に対する第２高調波の位相に任意性があるが、以下ではこの位相差が無い場合を例にとって説明する（実際はこの位相差が定まっていれば、どのような位相差を用いても良い。）。
【００５４】
第２高調波重畳送信処理１１０により、第２高調波を重畳して送信する、即ち、送信パルス信号として、ｓ１（ｔ）ｓｉｎ（ωｔ）＋ｓ２（ｔ）ｓｉｎ（２ωｔ）を送信する。
【００５５】
次に、直交検波処理１１１により受信信号を直交検波して、受信信号をＣ１（ｔ）ｃｏｓ（ωｔ）＋Ｓ１（ｔ）ｓｉｎ（ωｔ）＋Ｃ２（ｔ）ｃｏｓ（２ωｔ）＋Ｓ２（ｔ）ｓｉｎ（２ωｔ）と展開する。
【００５６】
２乗処理１１２により、Ｃ１（ｔ）＋ｊＳ１（ｔ）を２乗して、実部としてＣ３＝Ｃ１２（ｔ）−Ｓ１２（ｔ）、虚部としてＳ３＝２Ｃ１（ｔ）Ｓ１（ｔ）を求める。
【００５７】
次に、直交成分抽出処理１１３により、受信した第２高調波Ｃ２（ｔ）＋ｊＳ２（ｔ）から、Ｃ３（ｔ）＋ｊＳ３（ｔ）に直交する成分Ｒ（ｔ）＝（Ｃ２（ｔ）Ｓ３（ｔ）＋Ｓ２（ｔ）Ｃ３（ｔ））／ｓｑｒｔ（Ｃ３２（ｔ）＋Ｓ３２（ｔ））を求める。音響インピーダンス決定処理１１４により、Ｒ（ｔ）をもとに音響インピーダンス像を作る。
【００５８】
高周波を重畳して送信する場合、前記のように周波数空間での処理法以外にも、時間軸上での処理、すなわち送波波形と受波信号の相関をとって位相反転を検出する方法も有効である。もちろん、減衰に対応して、受波信号に対して相関をとる波形を徐々に変形してもよい。
【００５９】
次に、ビームの回折効果の補正方法について述べる。
【００６０】
図８は、本発明の実施例１において超音波パルス波形、点応答関数を説明する図である。
【００６１】
まず、現状の装置の超音波ビームの回折パターンと、それが超音波エコーの位相にどのような影響を与えているか説明する。現状のほとんどの超音波診断機において、電気信号の超音波への変換、超音波の電気信号への変換には、ＰＺＴ（チタン酸ジルコン酸鉛）を用いている。これは、感度、即ち、電気信号を超音波信号に変換する効率（これは超音波信号を電気信号に変換する効率とも等しい）の良さと、均質な材料を安定して生産できること、比較的機械加工がし易いなどの理由による。
【００６２】
しかし、感度が良いといっても、電気信号を超音波信号に変換する効率は６０から８０％程度である。本来はパルスエコー法を用いるには、デルタ関数的なパルスを送波することが好ましいが、なるべく送信エネルギーを稼ぐために、共振現象を使って、ＰＺＴの共振周波数で電気的に励起するということを使っているため、そこから出ていく超音波は、図８（ａ）に示す様にｓｉｎ波的なパルスとなる。
【００６３】
このような波形をアレイの各素子から送信すると、焦点近傍の音場は、図８（ｂ）に示す様になる。超音波断層像は、この音場を点応答関数とみなして、この点応答関数と、生体中の反射源との２次元コンボリューションを取ったものと見なすことができる。図８（ｂ）に示す点応答関数の形状を見てみると、深さ方向にキャリアの位相が回転しているのと、同時に方位方向にも位相が回転しているのがわかる。このキャリアによる局所的な深さ方向の位相回転に関しては、複素ビームフォーマ出力では除去されているが、方位方向の位相回転も位相イメージングにおいて誤った情報を与えることになる。反射体が送受信口径の焦点位置
から方位方向にはずれた位置にある場合、ずれの量に応じて位相が回転して、エコーが受信されてしまう。
【００６４】
この効果を補正するために、本発明においては、以下に説明する最小２乗法を使ったマトリックスフィルタを用いる。送波波形をＢ、マトリックスフィルタをｆ、マトリックスフィルタの出力信号をｃとすると、マトリックスフィルタ出力後の信号ｃは、（数１）で与えられる。
【００６５】
以下の説明で、記号Ｂは行列を示しｃ、ｆ、ｄはベクトル、記号「Ｔ」は転置を示す。ｄはマトリックスフィルタの出力として望まれる形を示すベクトルである。
【００６６】
（数１）でｆ＝（ｆ１，ｆ２，…，ｆｍ）であり、Ｂは、第１行を（ｂ１，ｂ２，…，ｂｎ，０，…，０）、第２行を（０，ｂ１，ｂ２，…，ｂｎ，０，…，０）、…、第ｎ行を（０，…，０，ｂ１，ｂ２，，…，ｂｎ）とするｍ行（ｍ＋ｎ−１）列の行列である。ｍはマトリックスフィルタのタップ数（要素数）であり、ｎは送波波形をサンプリングしてベクトル化した時のタップ数である。
ｃ＝ｆＢ　　　　　　　　　　　　　　　　　‥‥‥‥‥‥‥‥‥‥（数１）
マトリックスフィルタ透過後の信号ｃとし、望ましい波形をｄとすると、ｃとｄの２乗誤差和Ｉは、（数２）となる。２乗誤差和Ｉを最小にするｆがマトリックスフィルタである。加算Σは、ｉ＝１，２，…，ｍ＋ｎ−１について行なう。
Ｉ＝Σ（ｃｉ−ｄｉ）２＝（ｆＢ−ｄ）（ｆＢ−ｄ）Ｔ＝ｆＢＢＴｆＴ−ｄＢＴｆＴ−ｆＢｄＴ＋ｄｄＴ　　　‥‥‥‥‥‥‥（数２）
（数３）の条件から（数４）を全てのｉ（ｉ＝１，２，…，ｍ）について求めると（数５）となるので、ｆは（数６）のように求まる。
∂Ｉ／∂ｆｉ＝０　　　　　　　　　　　　　　　　‥‥‥‥‥‥‥‥（数３）
∂Ｉ／∂ｆｉ＝ＢｉＢＴｆＴ＋ｆＢＢｉＴ−ｄＢｉＴ−ＢｉｄＴ＝２（ｆＢＢｉＴ−ｄＢｉＴ）＝０　　　　‥‥‥‥‥‥‥（数４）
（∂Ｉ／∂ｆ１，∂Ｉ／∂ｆ２，…，∂Ｉ／∂ｆｍ）Ｔ＝ｆＢＢＴ−ｄＢＴ＝０‥‥‥‥‥‥‥‥‥（数５）
ｆ＝ｄＢＴ（ＢＢＴ）−１　　　　　　　　　　　‥‥　‥‥‥‥‥‥‥（数６）
このような最小２乗法に基づいたフィルタ（最小２乗フィルタ）を設計することで、受信信号の周波数帯域を保ったまま、波形を変更することが出来る。このことによって、ビームの方位方向の位相回転を補正することが出来る。信号の帯域を変更しないのは、帯域を変えて、元々信号の無い帯域にまで信号を作り出す操作を行なうと、ノイズが入ったときに、このノイズを増強してしまうことがあり、実際のノイズがある系では使えなくなってしまうからである。
【００６７】
図９は、本発明の実施例１において最小２乗フイルタ透過後の超音波パルス波形、方位方向の位相回転補正後の点応答関数を説明する図である。ここでは、図９（ａ）に示す様に、送波パルスの帯域の上限と下限に近い帯域を持つガウシアンパルスの差分波形を、マトリックスフィルタの出力信号として望ましい波形Ｄに用いる。これを用いて設計したフィルタを、図８（ｂ）に示す点応答関数に対して用いた結果が図９の（ｂ）であり、方位方向の位相回転が無くなっていることが解る。この方式を用いた装置構成例を、図１０に示す。
【００６８】
図１０は、本発明の実施例１において方位方向の位相回転補正を行なう超音波撮像装置の構成例を示す図である。複素受波ビームフォーマ２０の出力信号がマトリックスフィルタ２７に入力され、反射体が送受信口径の焦点位置から方位方向にはずれた位置にあり、ずれの量に応じて位相が回転してエコーが受信された場合、点応答関数の方位方向の位相回転補正がなされ、位相回転補正がなされた信号が音響インピーダンス変化量演算部２２に入力される。
【００６９】
図１０に示す超音波撮像装置の構成において、音響インピーダンス変化量演算部２２を音響インピーダンス演算部２５に代えることもできる。点応答関数の方位方向の位相回転補正は、点応答関数が対象域の探触子からの距離に応じて変ることから、マトリックスフィルタ２７もエコーの受信時間の関数として変える様に構成されることが好ましい。
【００７０】
実施例１では、周波数減衰起因の位相回転補正、方位方向の位相回転の補正に関して説明したが、実際には、これら補正を同時に使うことはより精確な位相補正を行なう上ではより好ましい。
【００７１】
（実施例２）
次に、本発明による位相補正を応用して、共振物を強調して画像化する方法に関して述べる。進行波に対する厚みが波長の半分もしくは、この奇数倍の厚みを持つ対象物が生体中に有る場合、対象物に超音波が到達すると、この対象物は共振を起こす。対象物が共振を起こすと、対象物からの信号は位相が回転することが知られている。よって、本発明の方法によって、位相回転の補正が行なわれると、この共振状態を抽出してイメージングすることが可能になる。この場合は、超音波スペクトロスコピー的に入力超音波周波数を掃引することで、系の周波数応答をイメージングすることが可能となる。
【００７２】
図１２は、本発明の実施例２において造影剤からの信号を抽出する超音波撮像装置の構成例を示す図である。図１２に示す超音波撮像装置の構成は、図１の超音波撮像装置の構成において、音響インピーダンス変化量演算部２２を、位相の空間変化を抽出する位相変化抽出部２８に置き換えた構成である。位相変化抽出部２８は、位相回転補正部２１の出力信号を使用して位相の空間変化を抽出する。特に共鳴系として、最近は以下に述べる超音波造影剤の画像化が重要である。
【００７３】
本発明によって、生体からの反射信号の位相変化が意味ある情報として取得出来るようになると、造影剤からの信号と生体組織からの信号を区別する事にも使うことが出来るようになる。即ち、生体からの反射信号は前述のように、柔らかいところから硬いところに変る界面もしくはその逆の界面で、反射信号の位相が０度もしくは１８０度となるが、造影剤の場合、より複雑な応答をして全位相成分を含んでいる。
【００７４】
図１１は、本発明の実施例２において、例えば、Ｒａｙｌｅｉｇｈ−Ｐｌｅｓｓｅｔ方程式を解いて、造影剤の応答（造影剤からの受波信号）のシミュレーション計算結果例を示す図である。図１１に示すように、点線で示す励起波に対して、実線で示す泡の応答の様子は、周波数もずれているが、位相もずれた形であることが解る。つまり周波数変化を用いなくても、同じ周波数の信号に対しても造影剤とそうでない信号に対してコントラストをつけることが可能となる。即ち、０度と１８０度以外の位相の信号を抽出することで、造影剤からの信号のみを抽出することが可能となる。
【００７５】
図１３は、本発明の実施例２において仮想的な高調波を想定して造影剤分布像を求める信号処理プロセス（仮想高調波位相法信号処理プロセス）を説明する図である。図１３の信号処理プロセスの各処理のうち、直交検波処理１００、位相成分演算処理１０１、仮想２乗成分演算処理１０２、ローパスフィルタ処理１０３、複素共役演算処理１０４、開平処理１０５、エコー位相回転補正処理１０６、符号の決定処理１０７は、図５に示す信号処理プロセスの各処理と同じである。
【００７６】
図１３の信号処理プロセスでは、図５に示す音響インピーダンス決定処理１０８に代えて、造影剤信号抽出処理１０９により、Ｓ５（ｔ）をもとに造影剤像を求める。
【００７７】
図１４は、本発明の実施例２において高調波を重畳して送信して造影剤分布像を求める信号処理プロセスを説明する図である。図１４の信号処理プロセスの各処理のうち、第２高調波重畳送信処理１１０（図１４に図示せず）、直交検波処理１１１、２乗処理１１２は、図７に示す信号処理プロセスの各処理と同じである。
【００７８】
図１４の信号処理プロセスでは、図７に示す直交成分抽出処理１１３に代えて同相成分抽出処理１１５を行ない、音響インピーダンス決定処理１１４に代えて造影剤信号抽出処理１１６を行なう。同相成分抽出処理１１５では、受信した第２高調波Ｃ２（ｔ）＋ｊＳ２（ｔ）からＣ３（ｔ）＋ｊＳ３（ｔ）に同相な成分Ｒ（ｔ）＝（Ｃ２（ｔ）Ｃ３（ｔ）−Ｓ２（ｔ）Ｓ３（ｔ））／ｓｑｒｔ（Ｃ３２（ｔ）＋Ｓ３２（ｔ））を求める。造影剤信号抽出処理１１６により、Ｒ（ｔ）をもとに造影剤分布像を作る。
【００７９】
本発明の方法によれば、周波数変化していない信号をも造影剤イメージングに用いる事が出来る。これまでは造影剤イメージングにおいては、周波数変化した分のみを用いて、周波数変化していない信号は切り捨ててきたが、これは一番信号の多い分を無駄にして来た事に相当するので、本方法を用いることによって、信号対ノイズの比が改善することが期待される。
【００８０】
（実施例３）
更に、本発明は、以下に示す様に、超音波による加熱凝固治療法のモニタリングに用いる事も可能である。
【００８１】
図１５、図１６、図１７は、本発明の実施例３において加熱凝固治療法のモニタに使用する超音波撮像装置の構成例を示す図である。
【００８２】
図１５に示す超音波撮像装置の構成では、画像１フレーム分、もしくは関心領域分の治療前の複素ビームフォーマ２０の出力信号を蓄えるメモリ２９と、複素ビームフォーマ２０の出力信号が入力され、位相の時間変化を演算する位相の時間変化演算部３０を具備しており、図１５に示す他の構成要素は図１と同じである。
【００８３】
図１６に示す超音波撮像装置の構成では、図１に示す超音波撮像装置の構成のうち音響インピーダンス変化量演算部２２に代えて、メモリ２９と位相の時間変化演算部３０を使用する。位相回転補正部２１の出力信号は、位相の時間変化演算部３０に入力されると共に、メモリ２９に蓄えられる。メモリ２９に蓄えられた信号は位相の時間変化演算部３０に読み出される。
【００８４】
図１７に示す超音波撮像装置の構成では、音響インピーダンス演算部２５、メモリ２９、音響インピーダンスの時間変化を演算する音響インピーダンス時間変化演算部３１を具備しており、図１７に示す他の構成要素は、図１と同じである。位相回転補正部２１の出力信号は、音響インピーダンス演算部２５に入力され、音響インピーダンス演算部２５の出力信号は、音響インピーダンス時間変化演算部３１に入力されると共に、メモリ２９に蓄えられる。メモリ２９に蓄えられた信号は音響インピーダンス時間変化演算部３１に読み出される。
【００８５】
図１５、図１６、図１７に示すように、画像１フレーム分、もしくは関心領域分の治療前の複素ビームフォーマ出力をメモリに蓄えておき、これと治療中の複素ビームフォーマ出力と比較することで治療部位の位相変化を画像化出来る。
【００８６】
この治療による変化の検出は、図１６に示す様に位相回転補正後に行なっても良いし、図１７に示す様に位相回転補正後、音響インピーダンスの時間変化として画像化することでも良い。また、治療前後で比較するのでなく、治療中の適切な画面と比較することで、治療中の体動もしくは超音波探触子の動きによるアーチファクトの影響を最小限に押さえることも可能となる。
【００８７】
図１８は、本発明の実施例３において加熱凝固治療装置と超音波撮像装置の組合せ例を示す図である。図１８に示す例では、図５に示す超音波撮像装置と超音波治療装置（システム）４０と一体化している。この一体化により、位相の時間変化演算部３０で求められた信号を治療装置４０の制御信号としてフイードバックして使用する。この結果、この治療効果領域が治療前に定めた予定の大きさになったら、自動的に治療を終了する機能を持たせる事も可能となる。
【００８８】
（実施例４）
本発明においては、位相補正された複素受波信号から、任意の位相成分を表示することで、実部から虚部へ連続的に画質を変化させることが可能である。
【００８９】
図１９は、図１の構成にて音響インピーダンス変化量演算部２２を特定位相成分抽出部３２に置き換えて、位相成分選択入力部３３を付加した構成を示す。図１９に示す様に、位相選択範囲を入力するための入力部を本体に備えることで、造影剤を強調して表示したり、音響インピーダンス等の組織特性の空間変化を強調して表示したりすることが可能となる。
【００９０】
この方法には、これまでに述べてきた第二高調波を重畳して送波することは有効であり、構成としては、図２０に示すような形態となる。図２０は、図６の構成にて音響インピーダンス変化量演算部２２を特定位相成分抽出部３２に置き換え、位相成分選択入力部３３を付加した構成を示す。
【００９１】
以上、上述した実施例で説明したように、本発明は、位相回転効果を補正し、真の生体内音響インピーダンス変化のイメージングを可能とする超音波撮像装置及び超音波撮像方法を実現することができる。
【００９２】
また、音響インピーダンスイメージングにおいて、エコーの周波数帯域移動による位相回転及びビームの方位方向の位相回転を補正することによって、検査対象内の音響インピーダンスの変化による、エコー信号の位相回転分を表示する超音波撮像装置の実現が可能となる。
【００９３】
さらに、上記した位相回転の計算から、組織の硬さの変化の表示や、共振部位の表示、造影剤からの信号の表示、凝固治療による組織の硬さの変化を表示することが可能となる。
【００９４】
なお、本発明は、上述した特定の実施の形態に限定されるものでなく、その技術思想の範囲を逸脱しない範囲で様々な変形が可能である。
【００９５】
【発明の効果】
本発明によれば、位相回転効果を補正し、真の生体内音響インピーダンス変化のイメージングを可能とする超音波撮像装置及び超音波撮像方法を提供できる。
【図面の簡単な説明】
【図１】本発明の実施例１において音響インピーダンスの微分画像を得る超音波撮像装置の構成例を示す図。
【図２】本発明の実施例１において音響インピーダンス画像を得る超音波撮像装置の構成例を示す図。
【図３】本発明の実施例１においてエコーの中心周波数の移動を説明する図。
【図４】本発明の実施例１においてエコーの位相に対する中心周波数の移動の影響を説明する図。
【図５】本発明の実施例１において仮想的な高調波を想定して音響インピーダンス像を求める信号処理プロセスを説明する図。
【図６】本発明の実施例１において高調波を重畳して送信する超音波撮像装置の構成例を示す図。
【図７】本発明の実施例１において高調波を重畳して送信して音響インピーダンス像を求める信号処理プロセスを説明する図。
【図８】本発明の実施例１において超音波パルス波形、点応答関数を説明する図。
【図９】本発明の実施例１において最小２乗フイルタ透過後の超音波パルス波形、方位方向の位相回転補正後の点応答関数を説明する図。
【図１０】本発明の実施例１において方位方向の位相回転補正を行なう超音波撮像装置の構成例を示す図。
【図１１】本発明の実施例２において造影剤の応答を計算した結果例を示す図。
【図１２】本発明の実施例２において造影剤からの信号を抽出する超音波撮像装置の構成例を示す図。
【図１３】本発明の実施例２において仮想的な高調波を想定して造影剤分布像を求める信号処理プロセスを説明する図。
【図１４】本発明の実施例２において高調波を重畳して送信して造影剤分布像を求める信号処理プロセスを説明する図。
【図１５】本発明の実施例３において加熱凝固治療法のモニタに使用する超音波撮像装置の構成例を示す図。
【図１６】本発明の実施例３において加熱凝固治療法のモニタに使用する超音波撮像装置の構成例を示す図。
【図１７】本発明の実施例３において加熱凝固治療法のモニタに使用する超音波撮像装置の構成例を示す図。
【図１８】本発明の実施例３において加熱凝固治療装置と超音波撮像装置の組合せ例を示す図。
【図１９】本発明の実施例４の一例を説明する図。
【図２０】本発明の実施例４の別の例を説明する図。
【符号の説明】
１…超音波探触子、３…送波ビームフォーマ、４…制御系、５…送受切り替えスイッチ郡、２０…複素受波ビームフォーマ、２１…位相回転補正部、２２…音響インピーダンス変化量演算部、２３…スキャンコンバータ、２４…表示部、２５…音響インピーダンス演算部、２６…送波波形メモリ、２７…マトリックスフィルタ、２８…位相変化抽出部、２９…メモリ、３０…位相の時間変化演算部、３１…音響インピーダンス時間変化演算部、４０…治療装置、１００…直交検波処理、１０１…位相成分演算処理、１０２…仮想２乗成分演算処理、１０３…ローパスフィルタ処理、１０４…複素共役演算処理、１０５…開平処理、１０６…エコー位相回転補正処理、１０７…符号の決定処理、１０８…音響インピーダンス決定処理、１０９…造影剤信号抽出処理、１１０…第２高調波重畳送信処理、１１１…直交検波処理、１１２…２乗処理、１１３…直交成分抽出処理、１１４…音響インピーダンス決定処理、１１５…同相成分抽出処理、１１６…造影剤信号抽出処理。
【図１】
2004113364.tif 000003
【図２】
2004113364.tif 000004
【図３】
2004113364.tif 000005
【図４】
2004113364.tif 000006
【図５】
2004113364.tif 000007
【図６】
2004113364.tif 000008
【図７】
2004113364.tif 000009
【図８】
2004113364.tif 000010
【図９】
2004113364.tif 000011
【図１０】
2004113364.tif 000012
【図１１】
2004113364.tif 000013
【図１２】
2004113364.tif 000014
【図１３】
2004113364.tif 000015
【図１４】
2004113364.tif 000016
【図１５】
2004113364.tif 000017
【図１６】
2004113364.tif 000018
【図１７】
2004113364.tif 000019
【図１８】
2004113364.tif 000020
【図１９】
2004113364.tif 000021
【図２０】
2004113364.tif 000022