(57)【特許請求の範囲】
【請求項１】
（ａ）測定ステップであって、該ステップの間、振動子Ｔｉのアレイ（２）が、撮像される媒体（１）の視野内に、前記視野内に焦点のない異なる波面を有する２以上の数Ｎの入射超音波ｌを放射し、各入射超音波は、前記振動子のアレイ（２）における複数の振動子Ｔｉによって放射され、Ｎは、２以上で、かつ１００未満であり、前記入射波ｌの伝搬の間、撮像される前記媒体と接触した少なくとも１つの電気センサ（Ｅｌ）が、生の電気信号Ｅｒａｗｌ（ｔ）のそれぞれを捕捉する、測定ステップと、
（ｂ）画像形成ステップであって、該ステップの間、電流のマップを含む前記媒体の画像が、ステップ（ａ）において得られた前記生の電気信号Ｅｒａｗｌ（ｔ）から決定される、画像形成ステップと
を含み、
前記生の電気信号は、前記媒体（１）中の、入射超音波と電流との相互作用から生じることを特徴とする、
音響電気撮像のための方法。
【請求項２】
ステップ（ｂ）の間、少なくとも前記Ｎの生の電気信号Ｅｒａｗｌ（ｔ）から、前記視野内のＭ個の仮想焦点Ｐｋに対して、電気値Ｅｃｏｈｅｒｅｎｔｋが決定され、各々が、仮に点Ｐｋに収束された超音波が前記振動子によって放射されたとしたら、補足されたであろう前記電気信号に対応する、請求項１に記載の方法。
【請求項３】
ステップ（ｂ）の間、逆ウェーブレット変換ＷＴ−１ 、次いで、逆ラドン変換Ｒ−１ が、前記生の電気信号Ｅｒａｗｌ（ｔ）に適用される、請求項１に記載の方法。
【請求項４】
ステップ（ｂ）の間、前記振動子のアレイ（２）を用いて得られた前記媒体の超音波画像が、前記電流のマップ上に重ね合わされる、請求項１または請求項２または請求項３に記載の方法。
【請求項５】
ステップ（ａ）の間、前記振動子Ｔｉが、それぞれ、前記入射波ｌから前記媒体によって反射された超音波を表す音響信号ＲＦｒａｗｌ、ｉ（ｔ）を捕捉し、
ステップ（ｂ）の間、捕捉された信号ＲＦｒａｗｌ、ｉ（ｔ）のＮの組から、仮に点Ｐｋに収束された超音波が前記振動子によって放射されたとしたら、前記振動子Ｔｉによって受け取られたであろう前記音響信号に対応するＭのコヒーレント音響信号ＲＦｃｏｈｅｒｅｎｔｋ、ｉ（ｔ）が決定され、前記媒体の前記超音波画像が、前記コヒーレント音響信号から計算される、請求項４に記載の方法。
【請求項６】
ステップ（ｂ）の間、前記超音波画像が、前記コヒーレント音響信号に基づくビームフォーミングによって決定される、請求項５に記載の方法。
【請求項７】
撮像される前記媒体が、ヒトまたは動物の組織である、請求項１から６のいずれか一項に記載の方法。
【請求項８】
振動子Ｔｉのアレイと、少なくとも１つの電気センサ（Ｅｌ）と、制御および処理手段（８、１０、４）とを備え、
（ａ）振動子Ｔｉの前記アレイ（２）に、異なる波面を有する数Ｎの非収束入射超音波ｌを、撮像される媒体（１）内に放射させ、各入射超音波は、振動子の前記アレイ（２）の中の複数の振動子Ｔｉによって放射され、Ｎは２以上で、かつ１００未満であり、および撮像される前記媒体と接触した少なくとも１つの電気センサに、前記入射波ｌの伝搬の間、生の電気信号Ｅｒａｗｌ（ｔ）のそれぞれを捕捉させ、
（ｂ）前記生の電気信号Ｅｒａｗｌ（ｔ）から電流のマップを含む前記媒体の画像を決定するように適応されており、
前記生の電気信号は、前記媒体（１）中の、入射超音波と電流との相互作用から生じることを特徴とする、請求項１から７のいずれか一項に記載の音響電気撮像のための方法を実現するデバイス。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、音響電気撮像のための方法およびデバイスに関する。
【背景技術】
【０００２】
心臓、骨格筋、および脳などの器官は、ニューロン内で情報を運ぶ、または筋肉もしくは心筋収縮を引き起こす電気インパルスによって連続的に通行される。多くの疾患を診断し、脳の機能診査を通じて脳のメカニズムを理解するために、これらのインパルスの伝搬を撮像できることは極めて重要である。
【０００３】
音響電気撮像は、超音波と電流との間の相互作用を利用して、超音波と組織との間の相互作用の点における、典型的には、収束超音波の焦点における、電流の値を決定する。
【０００４】
米国特許第8,057,390号は、収束超音波が線ごとに電流の画像を形成するように放射される音響電気撮像方法の例を開示する。この取得プロセスは遅く、結果として得られる電気信号が非常に弱いので高レベルの平均化が必要となるだけになおさらそうである。したがって、低フレームレートが得られる。
【０００５】
Kuchmentらは、「Synthetic focusing in ultrasound modulated tomography」Inverse problem and imaging、2009-10-01、1〜9頁、XP055116447、において、振動子が球面波を1つずつ放射する合成音響電気撮像のための方法を提案した。その結果は遅い処理である。さらに、入射超音波は、あまりに小さい振幅を有する。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】米国特許第8,057,390号
【特許文献２】EP2101191
【非特許文献】
【０００７】
【非特許文献１】Kuchmentら、「Synthetic focusing in ultrasound modulated tomography」Inverse problem and imaging、2009-10-01、1〜9頁、XP055116447
【非特許文献２】Montaldoら、「Coherent plane-wave compounding for very high frame rate ultrasonography and transient elastography」(IEEE Trans Ultrasound Ferroelectr Freq Control 2009 March、56(3)、489〜506頁)
【非特許文献３】E. J. Candes、「Ridgelets: theory and applications」Stanford University、1998
【発明の概要】
【発明が解決しようとする課題】
【０００８】
本発明はこの不利点を克服することが意図されている。
【課題を解決するための手段】
【０００９】
このために、本発明は、音響電気撮像のための方法を提案し、方法は、
（ａ）測定ステップであって、そのステップの間、振動子Ｔｉのアレイが、撮像される媒体の視野内に、視野内で収束されず、異なる波面を有する２以上の数Ｎの入射超音波ｌを放射し、各入射超音波が、振動子のアレイの中の複数の振動子Ｔｉによって放射され、ここでＮは、２以上で、かつ１００未満であり、撮像される媒体と接触した少なくとも１つの電気センサが、それぞれ、入射波ｌの伝搬の間、生の電気信号Ｅｒａｗｌ（ｔ）を捕捉する、測定ステップと、
（ｂ）画像形成ステップであって、そのステップの間、電流のマップ（言い換えれば、媒体の各点における局所電流密度を表す電気値のマップ）を含む媒体の画像が、ステップ（ａ）において得られた生の電気信号Ｅｒａｗｌ（ｔ）から決定される、画像形成ステップと
を含む。
【００１０】
これらの構成を用いて、人は観察される媒体における電気インパルスの超高速撮像を得ることができ、場合により、実時間においておよび1ミリメートルの分解能において組織内深くの電気インパルスの伝搬をフィルム撮影することができる。
【００１１】
本発明による方法の様々な実施形態において、以下の構成のうちの１つまたは複数を場合により使用することができる。
−  ステップ（ｂ）の間、少なくともＮの生の電気信号Ｅｒａｗｌ（ｔ）から、視野内のＭ個の仮想焦点Ｐｋに対して、各々、仮に点Ｐｋに収束された超音波が前記振動子によって放射されたとしたら、補足されたであろう電気信号に対応する、電気値Ｅｃｏｈｅｒｅｎｔｋが決定される。
−  ステップ（ｂ）の間、逆ウェーブレット変換ＷＴ−１ 、次いで、逆ラドン変換Ｒ−１ が生の電気信号Ｅｒａｗｌ（ｔ）に適用される（生の電気信号Ｅｒａｗｌ（ｔ）は、もちろん、予備処理を受けてから逆ラドン変換Ｒ−１ を受けることができる）。
−  ステップ（ｂ）の間、振動子アレイを用いて作成された媒体の超音波画像が、電流のマップ上に重ね合わせられる。
−  ステップ（ａ）の間、振動子Ｔｉは、それぞれ、入射波ｌからの媒体によって反射された超音波を表す音響信号ＲＦｒａｗｌ、ｉ（ｔ）を捕捉する。
−  ステップ（ｂ）の間、捕捉された信号ＲＦｒａｗｌ、ｉ（ｔ）のＮの組から、仮に点Ｐｋに収束された超音波が前記振動子によって放射されたとしたら、振動子Ｔｉによって受け取られたであろう音響信号に対応するＭのコヒーレント音響信号ＲＦｃｏｈｅｒｅｎｔｋ、ｉ（ｔ）が決定され、媒体の超音波画像がコヒーレント音響信号から計算される。
−  ステップ（ｂ）の間、超音波画像が、コヒーレント音響信号に基づくビームフォーミングによって決定される。
−撮像される媒体はヒトまたは動物の組織である。
【００１２】
本発明は、振動子Ｔｉのアレイと、少なくとも１つの電気センサと、制御および処理手段とを備える、音響電気撮像のための方法を実現するデバイスにも関し、
（ａ）振動子Ｔｉのアレイに、撮像される媒体内において、異なる波面を有する数Ｎの非収束入射超音波ｌを放射させるステップであって、各入射超音波が振動子のアレイの中の複数の振動子Ｔｉによって放射され、Ｎは２以上で、かつ１００未満である、放射させるステップ、および撮像される媒体と接触した少なくとも１つの電気センサに、それぞれ、入射波ｌの伝搬の間、生の電気信号Ｅｒａｗｌ（ｔ）を捕捉させるステップと、
（ｂ）生の電気信号Ｅｒａｗｌ（ｔ）から電流のマップを含む、媒体の画像を決定するステップとに適合される。
【００１３】
本発明の他の特徴および利点は、非限定例によっておよび添付の図面を参照して示される、その実施形態のうちの1つの以下の説明から明らかとなろう。
【図面の簡単な説明】
【００１４】
【図１】本発明の実施形態による方法を実現するためのデバイスの概略図である。
【図２】図1のデバイスの一部分のブロック図である。
【発明を実施するための形態】
【００１５】
様々な図において、同じ参照番号は同じまたは同様の要素を示す。
【００１６】
図1は、電気値の同時測定とともに、超音波圧縮波(例えば、0.2〜40MHzの周波数における)を放射し受け取ることによって、媒体1を撮像するステップに適合された音響電気撮像デバイスの例を示す。
【００１７】
撮像される媒体1は、具体的には患者または動物の組織、具体的には筋肉(心筋または他の)または脳から構成することができる。
【００１８】
撮像デバイスは、例えば、以下を備える。
-  例えば、数百の振動子を備え、撮像される媒体1内の視野(超音波によって走査される対象の領域)の2次元(2D)画像を得るステップに適合されたnの超音波振動子のアレイ2、
-  振動子アレイ2を制御するための、および振動子アレイによって捕捉された信号を取得するステップに適合された電子ベイ3または同様のもの、
-  電子ベイ3を制御し、前記捕捉された信号から得られた超音波画像を表示するためのコンピュータ4または同様のもの。
【００１９】
振動子アレイ2は、例えば、軸Xに沿って互いに隣接して配置された振動子の組によって形成された線形アレイでよく、X軸と垂直のZ軸は、視野内の深さ方向を表す。それに続いて、振動子は、Tiと表し、ここでiは軸Xに沿った各振動子の位置を示す添え字である。以下の説明は、このタイプの振動子アレイ2をその例として用いるが、2次元アレイを含む他の形の振動子アレイも本発明の範囲内で可能である。
【００２０】
デバイスは、例えば、電位の差を測定する2つの電極によって形成された少なくとも1つの電気センサEl(図2)をさらに備える。この電気センサは、有利には、振動子アレイ2に取り付けることができ、振動子アレイ2と同時に撮像される媒体1に接触するように適合することができる。
【００２１】
使用される電気センサElの数は、相対的に少なく、一般に10未満、好ましくは5未満、通常は1である。
【００２２】
図2に表すように、電子ベイ3は、例えば、以下を備えることができる。
-  個々に振動子アレイ2内のnの振動子Tiと、電気センサElとに接続されたn+1のアナログデジタルコンバータ5(A/Di〜A/De)、
-  それぞれ、nのアナログデジタルコンバータ5に接続されたn+1のバッファ6(Bi〜Be)、
-  バッファ6およびコンピュータ4と通信する中央処理装置8(CPU)、
-  中央処理装置8に接続されたメモリ9(MEM)、
-  中央処理装置8に接続されたデジタル信号プロセッサ10(DSP)。
【００２３】
n+1のアナログデジタルコンバータ5(A/Di〜A/De)は同じでよく、それはn+1のバッファ6(Bi〜Be)の場合もそうであり、したがって、使用されるデバイスは、単純に超高速の音響撮像に従来使用されているデバイスでよいことに留意されたい。
【００２４】
このデバイスにより、媒体1の音響電気撮像の方法を実現することが可能になり、方法は、具体的には、デジタル信号プロセッサ10によって補助される中央処理装置8によって実行される以下のステップ、
(a)測定するステップ(生データの送信/受信と保存)と、
(b)電気値のマップを含む媒体の画像を決定するステップと
を含む。
【００２５】
ステップ(a):測定するステップ(生データの送信/受信と保存):
振動子アレイ2および電気センサElは、媒体1と接触して配置され、数Nの入射超音波が振動子Tiによって媒体1中に放射される(Nは、例えば2〜100、具体的には5〜10でよい)。当該入射波は、非収束(より詳しくは、視野内で収束されない)であり、異なるそれぞれの波面、すなわち、異なる形状および/または異なる向きの波面を有する。有利には、入射波は、それらの伝搬Vの方向とZ軸との間で測定されたそれらのそれぞれの傾角θの角度、または空間の異なる点から生じたかのように放射された発散波を特徴とする、そのそれぞれの波面F(単一波の波面Fを図1に表す)が様々な異なる傾角を有する平面または発散波である。平面波の例は以下で検討する。
【００２６】
入射波は、一般に、1マイクロ秒未満のパルスであり、典型的には、中心周波数が約1〜10サイクルの超音波である。入射波の発射は、例えば、約50〜200マイクロ秒だけ離隔され得る。
【００２７】
各入射波は、媒体1内で入射波を反射する反射体に衝突する。反射された超音波は、アレイの振動子Tiによって捕捉される。したがって、各振動子Tiによって捕捉された信号は、入射波が放射時に収束されていないので、媒体1から全体として届く。同様に、電気センサElは、入射超音波の伝搬の間、電気信号E(t)を捕捉し、この電気信号は、各測定時に波面によって表される線全体に沿って、入射波と撮像される媒体1との相互作用から生じる。
【００２８】
nの振動子Tiによって捕捉された反射信号は、次いで、対応するアナログデジタルコンバータA/Diによってデジタル化され、対応するバッファBiに格納され、一方、電気信号はアナログデジタルコンバータA/Deによってデジタル化され、対応するバッファBeに格納される。各入射発射後バッファに格納されたこれらの信号は、以下、生データと称する。これらの生データは、入射超音波の発射l後、それぞれ、振動子Tiおよび電気センサElによって捕捉されたn+1の生の時間信号RFrawl,i(t)およびErawl(t)から構成される。
【００２９】
入射波の各発射1の後、バッファBi〜Beに格納された信号は、信号プロセッサ10のメモリ9に前記プロセッサによる処理のために転送される。したがって、ステップ(a)の終わりに、メモリ9は、n+1の生信号のNのアレイ(ベクトル)を含む。
【００３０】
ステップ(a)が500Hz以上などの高速で繰り返され、それは画像を得るのに使用される入射波の少ない数Nによって可能となる。
【００３１】
ステップ(b):電気値のマップを含む媒体の画像を決定するステップ:
このステップ(b)を実行するための2つの方法を以下に説明する。
【００３２】
ｂ１）第１の方法：コヒーレントデータの合成：
生データのＮのアレイから、それぞれ、視野のＭの点Ｐｋ（ｘ、ｚ）における（ｋは１〜Ｍの整数であり、ｘ、ｚはＸ、Ｚ軸上の点Ｐｋの座標である）合成コヒーレントデータのアレイ（ベクトル）の数Ｍがプロセッサ８によって計算される。これらの合成コヒーレントデータのＭのベクトルの各々は、それぞれ、仮に振動子が点Ｐｋに収束した入射波を放射していたとしたら、振動子Ｔｉによって補足されたであろう信号に対応するｎの時間信号ＲＦｃｏｈｅｒｅｎｔｋ、ｉ（ｔ）を含む。
【００３３】
コヒーレントデータのアレイは、例えば、具体的には文献EP2101191またはMontaldoらによる論文「Coherent plane-wave compounding for very high frame rate ultrasonography and transient elastography」(IEEE Trans Ultrasound Ferroelectr Freq Control 2009 March、56(3)、489〜506頁)に説明された原理により、媒体1全体にわたる超音波圧縮波に均一な伝搬速度cを想定することによって、得ることができる。
【００３４】
各発射lに対応する平面波の伝搬の方向が既知であり、伝搬速度vが既知であるので、プロセッサ8は、各点Pkに対する、入射波lから点Pkまでの伝搬時間τec(l,k)および点Pkから振動子Tiまでの反射波の伝搬時間τrec(l,k,i)を計算することができ、したがって、合計往復移動時間τ(l,k,i)=τec(l,Pk)+τrec(l,Pk,i)である。
【００３５】
次いで、虚焦点Pkに対応する、振動子Tiの空間的コヒーレント音響信号が、次式を用いて計算される。
【００３６】
【数１】
0006415555.tif 000002
【００３７】
ここで、B(l)は、入射波の各発射1の寄与の重み付け関数である(今回の場合、B(l)の値はすべて1に等しくてよい)。この信号RFcoherentkijは、各点Pkの単一の値を表示する。
【００３８】
同様に、人はコヒーレント電気信号Ecoherentkを計算することができる。
【００３９】
【数２】
0006415555.tif 000003
【００４０】
この電気値は、仮にＰｋに収束した入射超音波が放射されたとしたら、特に、音響電気画像を得るために、仮に十分な数の入射波が放射されたとしたら、例えば、高分解能画像画像を得るために、仮に４０〜１００の入射波が放射されたとしたら、電気センサＥｌによって測定されたであろうものである。
【００４１】
これらの値Ecoherentkは、上記の周知の音響電気撮像方法において捕捉された電気値と同様に、点Pkにおける電流を表し、したがって、視野内の電流のマップを提供する。
【００４２】
コヒーレントデータRFcoherentkのアレイおよび、場合により、値Ecoherentkは、次いで、例えば、特許EP2101191またはMontaldoらによる文献「Coherent plane-wave compounding for very high frame rate ultrasonography and transient elastography」(IEEE Trans Ultrasound Ferroelectr Freq Control 2009 March、56(3)、489〜506頁)に、例えば、説明されるように、場合により媒体1内の収差の作用を補正することによって改良することができる。
【００４３】
電流マップは、コンピュータ4の画面上に表示させることができ、場合により媒体1のBモード超音波画像上にまたは前記媒体1の何らかの他の画像上に、具体的には、例えば、前記文献EP2101191に説明されるように、受信モードでビームフォーミングによってアレイEcoherentkから得られた超音波画像上に重ね合わせることができる。
【００４４】
b2)第2の方法:ラドンおよびウェーブレット変換:
以下に説明するように、電気信号Erawk(t)から、点Pkにおける電流の局所値を直接送ることも可能である。
【００４５】
生の電気信号Erawk(t)を以下のようにモデル化することができる。
Erawk=∫volumeKρJ(x,y,z)ΔP(x,y,z)dxdydz         (2)
ここで、
Kは10-9Pa-1程度の相互作用定数であり、
ρは媒体の抵抗であり、
ΔPは圧力変動であり、
yは平面(X,Z)に垂直のY軸に沿った座標であり、
Jは検出された電流の密度分布、言い換えれば、電気センサElの電極の電極感度ベクトルの電流密度ベクトル倍のスカラー積である。
【００４６】
放射された超音波がパルス平面波であるので、ΔP(x,y,z)は、放射角度θと時間tとの関数として構成することができる。方向Yを無視することにより、本発明者らは以下を得る。
ΔP(x,z)=ΔP(-qsinθ+ctcosθ,qcosθ+ctsinθ)
ここで、qおよびctは、それぞれ、波面Fの方向と伝搬Vの方向との座標である。
【００４７】
放射超音波をディラックインパルスとして、すなわち、無限に短いインパルスとして考えることにより、音響電気信号は以下のようになる。
【００４８】
【数３】
0006415555.tif 000004
【００４９】
または同等に、
【００５０】
【数４】
0006415555.tif 000005
【００５１】
ここで、R[J]はラドン変換である。
【００５２】
実際には、入射波はディラックインパルスではなく、有限周波数帯のインパルス信号であり、それは結果としてラドン変換の変数ctに対して畳み込みとなる。
【００５３】
【数５】
0006415555.tif 000006
【００５４】
ここで、W(ct)は放射された波形であり、
【００５５】
【数６】
0006415555.tif 000007
【００５６】
は畳み込みの結果である。
【００５７】
例えば、典型的な超音波放射は、以下の畳み込みカーネルを生じる。
【００５８】
【数７】
0006415555.tif 000008
【００５９】
ここで、nおよびmは、振動子の周波数帯内で調整することができる。この畳み込みカーネルは、電流密度分布のリッジレット変換(E. J. Candes、「Ridgelets: theory and applications」Stanford University、1998)と同等である。
【００６０】
実際には、m=nであり、この畳み込みカーネルは、以下のパラメータ、すなわちa=nλ、b=ctおよびθを用いるリッジレット分解となる。
【００６１】
リッジレット分解は、パーセバル-プランシュレル関係式、再構成公式、線形不連続性から程遠いゆっくりと変動する物体のスパース表現など、いくつかの数学的特性を有し、ウェーブレット変換およびラドン変換の構成(composition)として表現することができる。
【００６２】
より詳しくは、ウェーブレットおよびリッジレット変換をそれぞれ、WT[.]およびRT[.]として表すことによって、以下を実証することができる。
【００６３】
【数８】
0006415555.tif 000009
【００６４】
ウェーブレット変換およびラドン変換の反転は、よく知られた問題である。正確な反転は、これらの2つの変換の、それぞれ、WT-1およびR-1として存在し、したがって、本発明者らは以下を得る。
【００６５】
【数９】
0006415555.tif 000010
【００６６】
実際には、反転は2つのステップで行われる。すなわち、まず、ウェーブレット変換WTの反転、次いでラドン変換Rの反転である。
【００６７】
これによって、撮像される媒体内の視野全体(入射波によって掃引される領域)にわたって電流密度マップが得られ、非常に高速の取得の後、そのようになり、電気インパルスの伝搬の実際の動画(movie)を得ることによって非常に速い電気現象の実時間監視が可能になる。
【００６８】
信号対雑音比(SNR)、分解能、およびフレームレートを最大化することも望ましい。
【００６９】
1つの取り組みは、入射波を可能な限り最短のインパルスの形で放射し、それによって、分解能を最適化することである。しかし、これは、非常に低いエネルギーを放射することに相当し、したがって、SNRが低くなる。
【００７０】
周波数帯をより長い放射(したがって、より多いエネルギー)に対応するサブバンドに分割することも可能である。理論的には、結果としてSNRの増加があるが、フレームレートが減少する(画像を形成するのに複数の放射を必要とするからである)。
【００７１】
最後に、第3の取り組みには、インパルス圧縮を実施するのに使用することができる「チャープ」を放射することが関与する。この取り組みは、フレームレートを維持しながらSNRを最大化する。
【００７２】
SNRは雑音の影響を制限することによって向上させることもできる。リッジレット変換が少数の大きい係数と多数の小さい係数とを有する電流密度分布を表すスパースベース(sparse basis)であるので、単純に閾値を得られた信号に適用することによって雑音を除去することができる。第1の取り組みは、「小さい」係数を除去する閾値から構成される。他の方法では、問題の物理特性を使用することも可能である。例えば、雑音を主に含む係数は、2つの対向する極性の放射に対して受け取られた信号の窓の間の相互相関によって識別することができる。さらに、これらの信号を、体系的アーチファクト(systemic artifact)を除去するために減算することができる。
【００７３】
ラドン変換の反転にはいくつかの技法が存在する。最も一般的なのは、おそらく、フィルタ逆投影であり、それには、逆投影(ビームフォーミングに対応する)の前にランプフィルタの適用が関与する。雑音レベルを増大させるこのステップを避けるために、入射超音波を、このフィルタを含む適切なインパルスの形で放射することも可能である。圧縮感知などの他の戦略も適切である。
【００７４】
さらに、視野の2次元(Bモード)超音波画像を、例えば、前記文献EP2101191に説明されているように、受信モードでビームフォーミングによって、さらに形成するために、アレイRFcoherentkを(方法b1で上記に説明したように)計算することができる。
【００７５】
この視野のBモード超音波画像(または他の何らかの画像、場合によっては超音波)を、場合により、以前に決定された電気値のマップ上に重ね合わせることができ、媒体の超音波画像および電流マップの両方をコンピュータ4の画面上に表示することができる。
【符号の説明】
【００７６】
1  媒体
2  超音波振動子のアレイ
3  電子ベイ
4  コンピュータ
5  アナログデジタルコンバータ
6  バッファ
8  中央処理装置
9  メモリ
10  デジタル信号プロセッサ
F  波面
El  電気センサ
l  入射超音波
Pk  点
Ti  振動子
V  伝搬の方向
X  軸
Y  軸
Z  軸
θ  放射角度
【図１】
0006415555.tif 000011
【図２】
0006415555.tif 000012