JP4448921B2 - 骨を境界壁として囲まれた実質的に均質な媒体内において所定の音響波場を形成するための非侵襲的な方法および撮影方法ならびにそれら方法を実施するための装置 - Google Patents

骨を境界壁として囲まれた実質的に均質な媒体内において所定の音響波場を形成するための非侵襲的な方法および撮影方法ならびにそれら方法を実施するための装置 Download PDF

Info

Publication number
JP4448921B2
JP4448921B2 JP2004532220A JP2004532220A JP4448921B2 JP 4448921 B2 JP4448921 B2 JP 4448921B2 JP 2004532220 A JP2004532220 A JP 2004532220A JP 2004532220 A JP2004532220 A JP 2004532220A JP 4448921 B2 JP4448921 B2 JP 4448921B2
Authority
JP
Japan
Prior art keywords
transducer array
transducer
boundary wall
acoustic
bone boundary
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
JP2004532220A
Other languages
English (en)
Other versions
JP2005537071A (ja
Inventor
ジャン−フランソワ・オーブリ
マシア・ファンク
ミカエル・タントル
Original Assignee
スーパー・ソニック・イマジン
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by スーパー・ソニック・イマジン filed Critical スーパー・ソニック・イマジン
Publication of JP2005537071A publication Critical patent/JP2005537071A/ja
Application granted granted Critical
Publication of JP4448921B2 publication Critical patent/JP4448921B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0808Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the brain
    • A61B8/0816Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the brain using echo-encephalography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/42Details of probe positioning or probe attachment to the patient
    • A61B8/4209Details of probe positioning or probe attachment to the patient by using holders, e.g. positioning frames
    • A61B8/4218Details of probe positioning or probe attachment to the patient by using holders, e.g. positioning frames characterised by articulated arms
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4444Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
    • A61B8/4455Features of the external shape of the probe, e.g. ergonomic aspects
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52046Techniques for image enhancement involving transmitter or receiver
    • G01S7/52049Techniques for image enhancement involving transmitter or receiver using correction of medium-induced phase aberration
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/107Measuring physical dimensions, e.g. size of the entire body or parts thereof
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7253Details of waveform analysis characterised by using transforms
    • A61B5/7257Details of waveform analysis characterised by using transforms using Fourier transforms
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/42Details of probe positioning or probe attachment to the patient
    • A61B8/4272Details of probe positioning or probe attachment to the patient involving the acoustic interface between the transducer and the tissue
    • A61B8/4281Details of probe positioning or probe attachment to the patient involving the acoustic interface between the transducer and the tissue characterised by sound-transmitting media or devices for coupling the transducer to the tissue
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61NELECTROTHERAPY; MAGNETOTHERAPY; RADIATION THERAPY; ULTRASOUND THERAPY
    • A61N7/00Ultrasound therapy
    • A61N7/02Localised ultrasound hyperthermia

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Veterinary Medicine (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Neurology (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Thermotherapy And Cooling Therapy Devices (AREA)

Description

本発明は、骨を境界壁として囲まれた実質的に均質な媒体内において所定の音響波場を形成するための非侵襲的な方法に関するものであり、また、そのような方法を使用した医療用の撮影方法に関するものであり、さらに、それら方法を実施するための装置に関するものである。
より詳細には、本発明は、複数のトランスデューサからなる少なくとも1つのアレイ(1つまたは複数のトランスデューサを備えてなり、可能であれば、複数のトランスデューサからなる複数のサブアレイを構成する)によって骨境界壁を介して音響信号を伝達することにより、骨境界壁によって(全体的にあるいは部分的に)囲まれた実質的に均質な媒体内において少なくとも1つの所定のターゲット音響波場を形成するための非侵襲的な方法に関するものである。
『実質的に均質な媒体』という用語は、本明細書においては、音響波の伝搬に関して実質的に均質な特性を有した媒体を意味するものとして、使用される。そのような実質的に均質な媒体は、例えば、脳によって構成することができ、その場合、骨がなす境界壁は、頭蓋骨によって構成される。適切である場合には、実質的に均質な媒体は、心臓によって、あるいは、心肺システムによって、構成される。この場合には、骨境界壁は、胸郭によって構成される。
脳の1つまたは複数のポイントに焦点合わせしたパルス波によって、あるいは、空間的および時間的な観点からより複雑な場によって、そのようなターゲット音響波場を構成し得ることが観測された。
国際公開第02/32316号パンフレットには、このような方法の一例が開示されており、開示された方法は、結果という観点からは満足のいくものである。
本発明の格別の目的は、そのタイプの方法を、さらに改良することであり、特に、音響波の焦点合わせを、より正確なものとすることである。
国際公開第02/32316号パンフレット
この目的のために、本発明においては、
−トレーニングステージを行い、このトレーニングステージにおいては、
1a)少なくとも部分的にX線を使用することによって、骨境界壁の3次元画像を撮影し、これにより、様々なポイントにおける骨境界壁の空隙度合いを代理するパラメータを獲得し;
1b)3次元画像に基づいて、骨障壁内における少なくとも密度と音響波の速度と音響波の吸収とに関する3次元マップを決定し;
1c)骨境界壁に対してのトランスデューサアレイの特定位置を決定し;
1d)少なくとも実質的に均質な媒体の少なくとも1つのポイントと、トランスデューサアレイのうちの少なくともいくつかのトランスデューサと、の間にわたっての音響波の少なくとも1つの伝搬を、伝搬に関する機械的モデルと、密度と音響波の速度と音響波の吸収とに関する3次元マップと、に基づいて、シミュレートし;
1e)このシミュレーションに基づいて、ターゲット音響波場が得られるよう、トランスデューサアレイのうちの少なくともいくつかのトランスデューサから放出すべき個々の音響信号を計算し;
−骨境界壁上におけるトランスデューサアレイの正確な位置決めを実際に行うというステージを行い、このステージにおいては、
2a)骨境界壁上において、トランスデューサアレイを、特定位置に対応するよう初期的に大まかに位置決めし;
2b)トランスデューサアレイのうちの少なくともいくつかのトランスデューサを使用して、エコー像の撮影を行い、これにより、骨境界壁に対してのトランスデューサアレイの位置を求め;
2c)このステップ2b)において得られた位置に基づいて、骨境界壁に対してのトランスデューサアレイの精細な位置決めを行う;
ことを特徴としている。
これにより、施術者がその後に用途に応じて使用し得るような個々の音響信号が得られ、これにより、実質的に均質な媒体内において施術者が決定した領域に対して実際に音響波を正確に焦点合わせすることができる(あるいは、より一般的には、所望の音響波場を非常に正確に形成することができる)。これにより、例えば、固定型のあるいは機能的な医療用の撮影を行ったり、および/または、温熱療法を行ったりすることができる。
本発明による方法の好ましい実施態様においては、付加的に、1つまたは複数の以下の特徴点を有することができる。
−ステップ2b)の際に、骨境界壁に対してのトランスデューサアレイの相対位置を、トランスデューサアレイのうちの少なくとも一部を使用してエコー像の撮影を行うことより骨境界壁の少なくとも一部の外側形状を決定し;その後、得られた外側形状と、骨境界壁に関する3次元画像と、を比較する;ことにより求めるという特徴点。
−ステップ1a)の際に、骨境界壁に対して位置決めデバイスを堅固に固定し、この位置決めデバイスを、X線の吸収に関して適合したもの(少なくとも部分的に)とし、これにより、初期的ステップを実施するという特徴点。
−ステップ1a)の際に、3次元画像を、骨境界壁上の位置決めデバイスの位置をも撮影したものとするという特徴点。
−ステップ2b)の際に、トランスデューサアレイと位置決めデバイスとの間の相対位置を、エコー像によって認識するという特徴点。
−トランスデューサアレイを、少なくとも1つのフレキシブル壁を有した流体充填タンク(流体は、液体またはゲル等とすることができる)の中に設置し、フレキシブル壁を、骨境界壁に対して押圧するという特徴点。
実質的に均質な媒体が、脳の少なくとも一部から構成され、骨境界壁が、脳を囲んでいる頭蓋骨の一部から構成されているという特徴点。
−音響波の周波数を、0.5MHz〜3MHzという範囲とするという特徴点。
−ステップ1e)の次に、ステップ1f)を行い、このステップ1f)においては、トランスデューサアレイから放出される音響信号を、シミュレートし;信号を、所望の音響波場が得られるよう個々の音響信号に基づいて決定し;そのような放出によって形成される音響波の伝搬を、シミュレートし;伝搬が、ある種の所定基準(特に、焦点合わせ特性、音響波による加熱によって局所的に到達可能な温度範囲、キャビテーションが存在しないこと、等)を満たすかどうかを検証する;という特徴点。
−ステップ1d)の際に、実質的に均質な媒体内の少なくとも1つのポイントから、トランスデューサアレイのうちの少なくともいくつかのトランスデューサに向けての、音響波の伝搬をシミュレートし(特に、実質的に均質な媒体中のポイントからの音響パルスの放出をシミュレートする);iを1〜nの中のある整数としかつnをトランスデューサアレイ内におけるトランスデューサの総数としたときに、トランスデューサアレイ内における各トランスデューサ(i)の位置へと到達して受領されたシミュレーション音響信号Ri(t)を決定し;ステップ1e)の際に、対象をなす各トランスデューサ(i)から放出すべき個々の音響信号Ei(t)を、ステップ1d)において決定された受領シミュレーション音響信号Ri(t)の時間反転Ri(−t)に比例するものとして、決定するという特徴点。
−ステップ1e)の際に、放出される音響信号Ei(t)を、Ei(t) = Gi×Ri(−t)によって決定され、ここで、Giは、ゲイン因子であって、トランスデューサ(i)ごとに異なるものとされ、これにより、骨境界壁内の損失を補償し得るものとされるという特徴点。
−少なくともいくつかのトランスデューサに対応したゲイン因子(Gi)を、受領シミュレーション音響信号Ri(t)の振幅の2乗に反比例した係数としたものとするという特徴点。
−ステップ1d)の際に、ステップ(1b)において決定された実際の吸収係数(τ)とは正反対であるような吸収係数(−τ)を骨境界壁の各ポイントにおいて有しているような音響波の吸収の関しての仮想的3次元マップを使用することによって、シミュレートを行い;ステップ1e)の際に、放出すべき個々の音響信号Ei(t)を、時間反転Ri(−t)に等しいものとして決定するという特徴点。
−トランスデューサアレイを、骨境界壁に対して押圧することを意図した少なくとも1つのフレキシブル壁を有した流体充填タンク(液体またはゲル等によって充填されている)の中に設置し、ステップ1c)において決定したトランスデューサの特定位置を、骨境界壁に対して接触しないものとし、
ステップ1d)の際に、
−実質的に均質な媒体の少なくとも1つのポイントから、トランスデューサアレイのうちの少なくともいくつかのトランスデューサに向けての、音響波の伝搬を決定し;iを1〜nの中のある整数としかつnをトランスデューサアレイ内におけるトランスデューサの総数としたときに、トランスデューサアレイ内における各トランスデューサ(i)の位置へと到達して受領されたシミュレーション音響信号Ri(t)を決定し;
−その後、各トランスデューサ(i)からの、信号Ri(t)の時間反転に対応した音響信号Ri(−t)の放出をシミュレートし、骨境界壁に接触した各ポイントのところにおいて各トランスデューサ(i)に対応した仮想的トランスデューサ(i)に対しての、流体中の伝搬をシミュレートし、z仮想的トランスデューサ(i)の位置へと到達して受領されたシミュレーション音響信号R’i(t)を決定し;
−その後、各仮想的トランスデューサ(i)からの、音響信号G’i×R’i(−t)の放出を、シミュレートし、ここで、R’i(−t)を、信号R’i(t)の時間反転とし、G’iを、少なくともいくつかの仮想的トランスデューサ(i)に関しての、信号R’i(−t)の振幅の2乗に反比例した係数とし、
−その後、トランスデューサ(i)に向けての流体中の伝搬をシミュレートし、トランスデューサ(i)の位置へと到達して受領されたシミュレーション音響信号R”i(t)を決定し;
ステップ1e)の際に、
−放出されるべき個々の音響信号Ei(t)を、受領シミュレーション音響信号R”i(t)の時間反転R”i(−t)に等しいものとして、決定する;
という特徴点。
−ステップ1d)の際に、トランスデューサアレイ内の少なくともいくつかのトランスデューサ(i)からの音響波パルスの放出をシミュレートし、ここで、iを、トランスデューサアレイ内のトランスデューサを表す指数とし、nを、トランスデューサの総数を表す正の整数とし、rを、1〜mの中のある整数とし、mを、参照ポイントの総数を表す正の整数としたときに、各トランスデューサ(i)から、実質的に均質な媒体内に位置した複数の参照ポイント(r)に向けての、音響波の伝搬をシミュレートし;実質的に均質な媒体内の各参照ポイント(r)に到達するものとしてシミュレートされたインパルス応答hri(t)を決定し;ステップ1e)においては、以下のサブステップを行い、
1e1)それぞれ周波数(ωk)を有したシミュレートされた複数のインパルス応答の各々に関して、周波数成分の数(p)を決定し、ここで、kを、周波数成分を表す1〜pの間の指数とし;
1e2)iを、1〜nの中の整数とし、rを、1〜mの中の整数とし、Hri(ωk)を、インパルス応答Hri(t)のフーリエ形式における周波数(ωk)における値としたときに、p伝達行列H(ωk)[=Hri(ωk)]を決定し;
1e3)各参照ポイント(r)に関し、n個の成分Ei(ωk,r)を、F(ωk,r)=H(ωk)×E(ωk,r)として決定し、ここで、E(ωk,r)=[Ei(ωk,r)]を、n個の成分Ei(ωk,r)を有したベクトルとし、F(ωk,r)=[Fl(ωk,r)]を、m個の成分Fl(ωk,r)を有したベクトルとし、lを、1〜mの中の整数とし、m個の成分Fl(ωk,r)を、ポイント(r)におけるかつ周波数(ωk)における所望のターゲット音響波場の生成に対応したものとする;
という特徴点。
−サブステップ(1e3)においては、少なくとも伝達行列H(ωk)の逆行列を求めることにより、p行列H−1(ωk)を計算し、実質的に均質な媒体内の各参照ポイント(r)に関して、E(ωk,r) = H−1(ωk)×F(ωk,r)という式を使用して、ベクトルE(ωk,r)を計算するという特徴点。
−ステップ1d)の際に、実質的に均質な媒体内の複数の参照ポイント(r)と、トランスデューサアレイ内の少なくともいくつかのトランスデューサ(i)と、の間において、インパルス応答hri(t)を決定し、ここで、iを、1〜nの中のトランスデューサを表す指数とし、nを、トランスデューサの数を表す整数とし、rを、1〜mの中の整数とし、mを、参照ポイントの数を表す整数とし、ステップ1e)の際に、エコー像を形成し得るよう、各参照ポイント(r)において受領されるトランスデューサアレイの少なくとも一部の焦点合わせ方法を決定する;という特徴点。
本発明は、また、エコー像によって医療用の撮影を行うための方法を提供するものであって、この方法においては、上述したようなターゲット音響波場を形成するための非侵襲的方法を行い、さらに撮影ステージを行い、この撮影ステージにおいては、トレーニングステージの時点で決定したような個々の音響信号を使用することにより、トランスデューサアレイの少なくとも一部を使用して、実質的に均質な媒体に関する少なくとも1つのエコー像を撮影する。
最後に、本発明は、さらに、上述した方法を実施し得るよう構成された装置を提供するものであって、この装置は、少なくとも、
−複数のトランスデューサからなるアレイであるとともに、実質的に均質な媒体を囲む骨境界壁の外側に配置されるものとされた、トランスデューサアレイと;
−様々なポイントにおける骨境界壁の空隙度合いを与えるものとしてX線によって撮影された3次元画像に基づいて、骨障壁内における少なくとも密度と音響波の速度と音響波の吸収とに関する3次元マップを決定するためのマップ形成手段と;
−実質的に均質な媒体の少なくとも1つのポイントと、トランスデューサアレイのうちの少なくともいくつかのトランスデューサと、の間にわたっての音響波の少なくとも1つの伝搬を、伝搬に関する機械的モデルと、密度と音響波の速度と音響波の吸収とに関する3次元マップと、に基づいて、骨境界壁に対してのトランスデューサアレイの特定位置の関数として、シミュレートするシミュレート手段と;
−このシミュレーションに基づいて、実質的に均質な媒体内においてターゲット音響波場が得られるよう、トランスデューサアレイのうちの少なくともいくつかのトランスデューサから放出すべき個々の音響信号を計算するための計算手段と;
−トランスデューサアレイのうちの少なくともいくつかのトランスデューサを使用して、エコー像の撮影を行い、これにより、骨境界壁に対してのトランスデューサアレイの位置決めを行うための、位置決め手段と;
−骨境界壁に対してトランスデューサアレイが特定位置に対応するようにして骨境界壁に対して配置されたときに、トランスデューサアレイの位置の関数として、骨境界壁に対してのトランスデューサアレイの初期位置を、精密に位置決めするための精密位置決め手段と;
を具備している。
装置の好ましい実施形態においては、付加的に、1つまたは複数の以下の特徴点を有することができる。
−骨境界壁に対してのトランスデューサアレイの位置決めを行うための位置決め手段が、トランスデューサアレイのうちの少なくとも一部を使用してエコー像の撮影を行うことより骨境界壁の少なくとも一部の外側形状を決定し、その後、得られた外側形状と、骨境界壁に関する3次元画像と、を比較することにより、位置決めを行うものとされているという特徴点。
−位置決めデバイスと、この位置決めデバイスを骨境界壁上に堅固に固定し得る固定手段と、を具備し、位置決めデバイスは、X線を吸収し得るものとされ、これにより、骨境界壁の3次元画像内において可視のものとされ、骨境界壁に対してのトランスデューサアレイの位置決めを行うための位置決め手段が、エコー像の撮影によって、位置決めデバイスに対してのトランスデューサアレイの位置を位置決めするという特徴点。
−トランスデューサアレイが、骨境界壁に対して押圧されることとなる少なくとも1つのフレキシブル壁を有した流体充填タンク(液体またはゲル等によって充填されている)の中に、設置されるという特徴点。
トランスデューサアレイが、撮影のためのサブアレイと、温熱療法のためのサブアレイと、を備え、これら2つのサブアレイは、互いに異なるタイプのトランスデューサを有しているという特徴点。
本発明の他の特徴点や利点は、添付図面を参照しつつ、本発明を何ら限定するものではなく単なる例示としての好ましい実施形態に関する以下の詳細な説明を読むことにより、明瞭となるであろう。
添付図面においては、同じ部材または同様の部材に対しては、同じ符号が付されている。
音響波を生成するための装置(1)が、図1に示されている。この装置(1)は、患者(P)の身体の実質的に均質な媒体の中に所定の音響波場を生成し得るよう構成されている。実質的に均質な媒体は、骨がなす境界壁によって、少なくとも部分的に囲まれている。
図示の例においては、実質的に均質な媒体は、患者の脳によって構成されており、骨がなす境界壁は、患者の頭蓋によって構成されている。
より一般的には、実質的に均質な媒体は、ヒトの患者のあるいは任意の他の脊椎動物の、任意の組織媒体とすることができ、この組織媒体は、実質的に均質であるような音響波伝搬特性を有している。組織媒体は、可能であれば、患者(P)の心臓によって、あるいは、患者(P)の心肺システムによって、構成される。この場合には、骨境界壁は、胸郭によって構成される。
装置(1)は、患者(P)の脳内において超音波を発生させるように作用する。例えば、頭蓋骨の外部から、0.5MHz〜3MHzという程度の周波数の音響波を発生させるように作用する。
音響波の生成は、例えば、以下のことを意図している。
−固定型のあるいは機能的な撮影器具を使用して特にドップラー式撮影器具を使用して脳のエコー像あるいは一連をなすエコー像を撮影しこれにより血液の流れを可視化すること、あるいは、熱的な撮影を行いこれにより温熱療法によって引き起こされた加熱状態を可視化すること。
−温熱療法を行うことによって、特に、
−1つであるか複数であるかにかかわらず、良性または有害の腫瘍を破壊すること。
−出血を凝固させること(上述したドップラー式機能的撮影器具によって位置決めされる)。
−熱的に活性化可能な薬物を局所的に活性化すること。
−血液脳関門を局所的に破壊し、これにより、静脈内を経由して予め注入しておいた薬物を、局所的に拡散させること。
すべての場合において、可能な最大の正確さでもって、患者(P)の脳内において1つまたは複数の所定ターゲット音響波場を形成し得ることが必要である。これにより、例えば、複数のトランスデューサからなるアレイによって放出された音響波を脳内の1つまたは複数のポイントに焦点合わせすることができる、また、より複雑な音響波場を生成することができる。
複数のトランスデューサからなるアレイ(図1には示されていない)は、例えば、キャップ(2)内に組み込むことができる。キャップ(2)は、患者(P)の頭部(3)に対して押圧されるものであり、ロボットアーム(4)(あるいは、他の任意の配置システム)によって支持されている。ロボットアーム(4)は、互いに関節結合された複数のレバーアーム(5)を備えており、床に対して固定されたスタンド(6)によって支持されている。ロボットアームの基準座標系に対して患者の頭部(3)を正確に位置決めするために、スタンド(6)に対して固定された支持体(6a)を設けることができる。支持体(6a)上には、剛直な定位固定フレーム(8)が固定されている。定位固定フレーム(8)は、患者(P)がテーブル(9)上に位置しているときには、患者(P)の頭蓋骨に対して堅固に固定される。
ロボットアーム(6)あるいは他の位置決めデバイスは、有利には、マイクロコンピュータ(7)あるいは類似物によって制御される。マイクロコンピュータ(7)には、キーボード(7a)のような入力インターフェースと、スクリーン(7b)のような少なくとも1つの出力インターフェースと、が設けられている。マイクロコンピュータ(7)は、さらに、キャップ(2)内に組み込まれているトランスデューサアレイの動作を、直接的に、あるいは好ましくは信号処理のためのエレクトロニクスラック(10)(B)を介して、制御する。
図2に示すように、キャップ(2)は、例えば、剛直なカップ(11)を備えることができる。このカップ(11)は、フレキシブル壁(12)と協働することによって、タンク(13)を形成することができる。タンク(13)は、例えばゲルや液体といったような流体によって、特に水によって、充填されている。トランスデューサアレイ(17)は、カップ(11)の内面上においてこの流体に浸漬されている。
キャップ(2)が、ロボットアーム(4)によって患者の頭部(3)の頂部に対して配置された場合、フレキシブル壁(12)は、患者の頭蓋骨(14)の形状に対して密接に適合し、可能であれば、定位固定フレーム(8)の剛直保持アーム(15)の上端部をカバーする。剛直保持アーム(15)は、ネジ(16)によって頭蓋骨(14)に対して固定することができる。
複数のトランスデューサからなるアレイ(17)は、画像撮影の目的のためにのみ使用される総数n個の超音波トランスデューサを備えることができる、あるいは、温熱療法の目的のためにのみ使用されるトランスデューサを備えることができる、あるいは、実際に、画像撮影の目的と温熱療法の目的との双方のために使用されるトランスデューサを備えることができる。
付加的には、複数のトランスデューサからなるアレイ(17)は、以下のような2つのサブアレイを備えることができる。
−温熱療法のためのサブアレイ(17a)。このサブアレイ(17a)は、1以上とされた個数(n1)の、例えば100個以上や200個以上といったような、トランスデューサ(T1,T2,T3,…,Tn1)を備えている。これらトランスデューサは、流体を介して、患者の頭蓋骨(14)に対して音響的に連通した状態とされている。
−画像撮影のためのサブアレイ(17b)。このサブアレイ(17b)は、1以上とされた個数(n2)の(よって、n1+n2=n)、例えば100個以上や200個以上といったような、トランスデューサ(T’1,T’2,T’3,…,T’n2)を備えている。これらトランスデューサは、流体を介して患者の頭蓋骨(14)に対して音響的に連通した状態とされているとともに、例えば、患者の頭部(3)の矢状面内に実質的に位置している中央ストリップとして集合させることができる。
サブアレイ(17a,17b)の各トランスデューサは、有利には、互いに異なるタイプのものとされる。すなわち、サブアレイ(17a)のトランスデューサは、サブアレイ(17b)のトランスデューサよりも、大きなものとすることができ(例えば、サブアレイ(17a)のトランスデューサの直径は、約8mmとされ、サブアレイ(17b)のトランスデューサの直径は、約1mmとされる)、これにより、より大きなパワーレベルでもって、音響波を伝達することができる。
図3に示すように、エレクトロニクスラック(10)は、トランスデューサ(T1〜Tn1,T’1〜T’n2)を制御するものであって、以下のものを備えることができる。 −ロボットアーム(4)(POS)をも制御するマイクロコンピュータ(7)によって制御される中央演算処理装置(CPU)。
−中央演算処理装置(CPU)に対して接続された少なくとも1つの主メモリ(M)。
−トランスデューサ(T1〜Tn1,T’1〜T’n2)に対してそれぞれ接続されたサンプラ(E1,E2,…,En1,E’1,E’2,…,E’n2)。
−サンプラ(E1〜En1,E’1〜E’n2)に連通した複数のプロセッサあるいは中央演算処理装置(C1,C2,C3,…,Cn1,C’1,C’2,C’3,…,C’n2)。
−プロセッサ(C1〜Cn1,C’1〜C’n2)に対してそれぞれ接続されたメモリ(M1,M2,M3,…,Mn1,M’1,M’2,M’3,…,M’n2)。
上述した装置を動作させるには、まず最初に、定位固定フレーム(8)を、患者の頭蓋骨(14)の外側に固定する。
患者の頭蓋骨(14)に対して定位固定フレーム(8)を堅固に保持した後に、患者の頭蓋骨(14)の3次元画像を撮影する。この撮影は、例えば、コンピュータ断層撮影用スキャナ(図示せず)や他のタイプのスキャナを使用して行う。これにより、頭蓋骨(14)の各ポイントにおいて、空隙度合い(Φ)の3次元マップ(あるいは、空隙度合いを代理するパラメータに関しての3次元マップ)を形成することができる。この画像は、有利には、微細な解像度で撮影される。例えば、約0.2mmという解像度で撮影される。
この空隙度合いマップを得るために、断層スキャナによって得られた生のデータは、次式によって、ハウンズフィールド値へと変換することができる。
Figure 0004448921
ここで、μ、μbone、μwaterは、走査された組織と骨と水のそれぞれに関しての、X線に対しての光電子線減弱係数を示している。
X線の吸収の観点からは、頭蓋骨は、空隙部分が水によって充填された骨から形成された組織と見なすことができ、頭蓋骨の空隙度合い(Φ)は、次式によって係数(μ )に対して関連しているものと見なすことができる。
(2) μ = Φμwater + (1−Φ)μbone
よって、頭蓋骨(14)の各ポイントにおける空隙度合いは、次式によって、値(H)に対して直接的に関連している。
Figure 0004448921
よって、スキャナによって与えられた画像から、頭蓋骨(14)の空隙度合い、あるいは、それを代理するパラメータ(例えば、密度。密度は、場合によっては、ある種のスキャナの場合には、スキャナ画像から直接的に得られる)、を決定することができる。
頭蓋骨(14)の各ポイントに関する空隙度合いの値(Φ)に基づいて、頭蓋骨(14)のすべてのポイントにおいて、密度(ρ)と、音響波の速度(c)と、音響波の吸収(τ)と、に関しての3次元マップを決定することができる。
密度(ρ)は、次式を使用して計算することができる。
(4) ρ = Φ×ρwater + (1−Φ)×ρbone
ここで、ρwater は、水の密度(1000kg/m )であり、ρboneは、皮質骨の最大密度であって、約210kg/m と見積もることができる。
考慮している周波数における音響波の速度は、次式によって評価することができる。
(5) C = Cmin + (Cmax−Cmin)×(1−Φ)
ここで、Cmin は、水中における音響波の速度に対応しており(約1.5mm/μs)、Cmax は、皮質骨内における音響波の速度である(約2.9mm/μs)。
頭蓋骨(14)の各ポイントにおける吸収は、次式により与えることができる。
Figure 0004448921
ここで、τmin は、水中における音響波の最小吸収に対応しており(例えば、およそ0.1dB/mm〜0.5dB/mm、特に、約0.2dB/mm)、τmax は、皮質骨内における音響波の最大吸収に対応しており(例えば、およそ0.5dB/mm〜15dB/mm、特に、約8dB/mm)、βは、一定の係数であって、例えば、0.3〜0.7という範囲であり、特に、約0.5である。
密度や速度や吸収に関するこれら様々な3次元マップは、例えば、上記のマイクロコンピュータ(7)内においてあるいは他の任意の計算ユニット内において、コンピュータ断層撮影や他の走査によって得られた3次元画像に基づいて、計算することができる。コンピュータ断層撮影や他の走査によって得られた画像が、また、頭蓋骨(14)に対しての定位固定フレーム(8)の正確な位置を与えることに、注意されたい。
上記様々なマップ間における差に基づいて、マイクロコンピュータ(7)あるいは他の計算ユニットは、脳(19)を均質な媒体と見なすことによって、頭蓋骨(14)および脳(19)を通しての音響波の伝搬をシミュレートすることができる。脳(19)は、第1近似においては、超音波に対する振舞いに関しては、水と見なすことができる。
このシミュレーションを実行するために、使用者は、頭蓋骨(14)の表面上におけるトランスデューサアレイ(17)の位置を、マイクロコンピュータに対して入力する。これにより、マイクロコンピュータは、参照ポイントと称されるような脳(19)内の様々なポイント(18)と、アレイ(17)の様々なトランスデューサ(T1〜Tn1,T’1〜T’n2)によって占められた位置と、の間にわたっての超音波の伝搬をシミュレートすることができる。
このシミュレーションは、次式(7)のような波動方程式を使用することにより、実行することができる。
Figure 0004448921
ここで、rは、対象をなすポイントに関する位置ベクトルを示しており、pは、圧力を示しており、Sは、対象をなすポイントに存在するかもしれない音源によって生成された音響信号を示している。
頭蓋骨(14)および脳(19)内における超音波の伝搬は、上記方程式(7)を分離形式とすることにより、差分によってマイクロコンピュータ(7)の中でシミュレートすることができる。このシミュレーションは、また、有限要素法や、あるいは、パルス回折法、によっても実行することができる。
有利には、波動伝搬の計算は、脳(19)内における放射線追跡だけによって、および、頭蓋骨(14)およびその近傍における差分だけによって、計算を実行することにより軽減化することができる。その場合、計算の2つの手法に対応した領域間の界面のところにおいて、球状の音響波の放出が、シミュレートされる(球形波は、脳に向けて放出され、これにより、トランスデューサから頭蓋骨に向かうような波の伝搬をシミュレートすることができ、これにより、脳から来る波の伝搬をシミュレートすることができる)。
このシミュレーションにより、非侵襲的な手法で、つまり、患者の脳(19)内に必ずしも侵入する必要なく(そのような侵入は、必然的に排除されるわけではない。例えば生検を行う必要がある場合などは、そのような侵入を行うこととなる)、アレイ(17)の様々なトランスデューサ(T1〜Tn)によって放出されるする個々の音響波を決定することができ、これにより、脳(19)内に所定のターゲット音響波場を形成することができる。
例示するならば、ターゲット音響波場は、脳(19)の1つまたは複数の参照ポイント(18)上へと焦点合わせすることができる、および/または、より複雑な音響波場によってこうせいすることができる。すべての場合において、放出されるべき個々の信号は、脳内の様々な参照ポイント(18)上へと音響波を正確に焦点合わせし得るようにして、決定される。これにより、その後、様々な個々の音響信号を組み合わせることによって、より複雑な音響波場を形成することができる。
この決定は、様々な方法によって行うことができ、それらの方法においては、
−脳内の1つまたは複数の参照ポイント(18)と、アレイ(17)のトランスデューサ(T1〜Tn1,T’1〜T’n2)のうちの少なくともいくつかと、の間にわたっての音響波の伝搬をシミュレートする(脳(19)からトランスデューサアレイ(17)に向けて、および/または、アレイ(17)から脳(19)に向けて)ステップと、すなわち、この例においては、好ましくは、第1に、参照ポイント(18)と処理サブアレイ(17a)の各トランスデューサとの間において、および、第2に、参照ポイント(18)と撮影サブアレイ(17b)の各トランスデューサとの間において、シミュレートするステップと;
−上記シミュレーションを基礎として使用することにより、脳(19)内においてターゲット音響波場が得られるよう、想定している状況下において、アレイ(17)の各トランスデューサが放出すべき音響信号を計算するステップと;
を備えている。
すべての状況下において、脳(19)内の1つまたは複数のポイント(18)においてターゲット音響波場を得るために各トランスデューサ(T1〜Tn)によって放出すべき個々の音響信号の決定は、利用可能なコンピュータの計算能力に依存して、および、考慮している参照ポイント(18)の数に依存して、数十分〜数時間という時間範囲内に得ることができる。このようにして得られた音響信号を使用することにより、特に、
例えば約1mmという直径を有した焦点スポット上へと焦点合わせされるような非常に正確な音響波を得ることができる。これにより、エコー像に関して優秀な撮影精度を得ることができるとともに、優秀な温熱療法を行うことができる(例えば、放出線療法によって処置し得る腫瘍の精度は、数cmの程度である、すなわち、正確さがかなり悪いものである。また、放射線療法の場合には、温熱療法とは違って、副作用を引き起こす)。
また、各トランスデューサ(T1〜Tn1,T’1〜T’n2)から放出すべき個々の音響信号を決定するというステップの後に、本発明による方法および装置によれば、脳(19)内における1つまたは複数のターゲット音響波場に対応した音響信号をトランスデューサアレイ(17)によって仮想的に放出するという試験を行うことができることに注意されたい。これにより、頭蓋骨(14)および脳(19)内における音響波の伝搬が、ある種の所定の基準を満たすかどうかを、計算によって、検証することができる。特に、焦点合わせの精度に関する基準や、音響波による加熱に基づいて特定の温度範囲が得られるかどうかに関する基準や、傷害を引き起こしかねないようなキャビテーションが発生しないために基準、等を満たすかどうかを、計算によって、検証することができる。
超音波によって供給された熱に基づく温度に関しての時間的かつ空間的な変動をシミュレートするに際しては、例えば、次の拡散方程式を使用することができる。
Figure 0004448921
ここで、C は、媒体の熱容量を示しており、Kは、媒体の熱伝導度を示しており、μは、媒体の吸収係数を示している。拡散方程式において発生源に相当する項は、音響波に基づく圧力場によって媒体内へと供給された熱に対応している。
媒体(頭蓋骨+脳)内の圧力p(r,t)の3次元的分散が、デジタル的にシミュレートした音響波伝搬によって、既にモデル化されていることにより、このデータは、差分によって、上記熱拡散方程式のデジタルシミュレーション内に導入することができる。これにより、処理を実際に行うよりも前に、処理アレイ(17a)による温熱療法時にどのように温度分配が得られるかを予測することができ、例えば、骨内においておよび維持される必要のあるすべての敏感な領域において、温度が、過度に上昇しないことをチェックすることができる。
ターゲット音響波場を得るための上記個々の音響信号は、様々な方法で、決定することができる。
第1の方法においては、
−脳の少なくとも1つの参照ポイント(18)から、トランスデューササブアレイの各トランスデューサ(T1〜Tn1、あるいは、T’1〜T’n2)への、考慮している条件下においての、音響波の伝搬をシミュレートし、iを1〜nの中のある整数としかつnをトランスデューサアレイ(17)内におけるトランスデューサの総数としたときに、トランスデューサアレイ内において対象をなすサブアレイの各トランスデューサ(i)の位置へと到達して受領されたシミュレーション音響信号Ri(t)を決定し;
−放出される個々の音響信号Ei(t)を、上記のようにして決定された受領音響信号Ri(t)の時間反転Ri(−t)に比例するものとして、決定する。
よって、放出される個々の音響信号Ei(t)は、次式によって計算される。
Ei(t) = Gi×Ri(−t)
ここで、Giは、ゲイン因子であって、
−すべてのトランスデューサに関して同一のもの(可能であれば、1に等しい)とされる、あるいは、
−トランスデューサごとに異なるものとされ、頭蓋骨内の損失を補償し得るものとされる。
ゲイン因子(Gi)が、トランスデューサごとに異なる場合には、それらゲインは、以下の式を使用して計算することができる。
Figure 0004448921
ここで、aは、トランスデューサアレイ内のすべてのトランスデューサに共通の実数であり、Max(|Ri(t)|)は、信号Ri(t)の最大振幅の値である。有利には、ゲイン(Gi)は、十分な振幅(例えば、Max(|Ri(t)|)が、対象をなすすべてのトランスデューサiに関して、最大値Max(|Ri(t)|)の10%引きよりも大きなものである場合)を受領するようなトランスデューサに関してのみ上記の値となる。そうでなければ1に等しい。
第2の方法においては、Gi=1として、上記の手順を行う。しかしながら、仮想的な3次元音響波吸収マップを使用することにより、音響波の伝搬をシミュレートする。この場合、骨境界壁は、ステップ1b)において決定された実際の吸収係数(τ)とは正反対であるような吸収係数(−τ)を有している。これにより、媒体内における損失を自動的に補償される。
第3の方法においては、手順は、初期的には上記第1方法と同じであり、上記信号Ri(t)の計算に関しては、上記第1方法と同じである。
−その後、考慮している条件下において、各トランスデューサ(i)からの、信号Ri(t)の時間反転に対応した音響信号Ri(−t)の放出をシミュレートし、頭蓋骨に接触した各ポイントのところにおいて各トランスデューサ(i)に対応した仮想的トランスデューサ(Fi)(図2参照)に対しての、タンク(13)内の流体中の伝搬をシミュレートし、仮想的トランスデューサ(i)の位置へと到達して受領されたシミュレーション音響信号R’i(t)を決定する。
−R’i(−t)を、信号Ri(t)の時間反転とし、かつ、G’iを、少なくともいくつかの仮想的トランスデューサ(i)に関しての、信号R’i(−t)の振幅の2乗に反比例した係数としたときに、各仮想的トランスデューサ(i)からの音響信号G’i×R’i(−t)の放出をシミュレートする。ここで、G’iは、例えば、次式によって表される。
Figure 0004448921
ここで、aは、トランスデューサアレイ内のすべてのトランスデューサに共通の実数であり、Max(|Ri(t)|)は、信号Ri(t)の最大振幅の値である。有利には、ゲイン(Gi)は、十分な振幅(例えば、Max(|Ri(t)|)が、対象をなすすべてのトランスデューサiに関して、最大値Max(|Ri(t)|)の10%引きよりも大きなものである場合)を受領するようなトランスデューサに関してのみ上記の値となる。そうでなければ1に等しい。
−その後、トランスデューサ(i)に向けての流体中の伝搬をシミュレートし、トランスデューサ(i)の位置へと到達して受領されたシミュレーション音響信号R”i(t)を決定する。
−放出されるべき個々の音響信号Ei(t)を、受領音響信号R”i(t)の時間反転R”i(−t)に等しいものとして、決定する。
第4の方法においては、
−アレイ(17)内のあるいはサブアレイ(17a,17b)内の各トランスデューサ(i)からの音響波パルスの放出をシミュレートし、ここで、iを1〜nの中のある整数としかつnをトランスデューサアレイ(17)内におけるトランスデューサの総数を表す正の整数とし、rを1〜mの中のある整数とし、mを参照ポイント(18)の総数を表す正の整数としたときに、音響波は、考慮している条件下で、各トランスデューサ(i)から、脳内に位置した複数の参照ポイント(r)に向けて、伝搬するものとされ、
−脳内の各参照ポイント(r)に到達するものとしてシミュレートされたインパルス応答hri(t)を決定し、
−以下のサブステップによって、個々の音響信号を計算し、
1e1)それぞれ周波数(ωk)を有したシミュレートされた複数のインパルス応答の各々に関して、周波数成分の数(p)を決定し、ここで、kを、周波数成分を表す1〜pの間の指数(インデックス)とし;
1e2)iを、1〜nの中の整数とし、rを、1〜mの中の整数とし、Hri(ωk)を、インパルス応答Hri(t)のフーリエ形式における周波数(ωk)における値としたときに、p伝達行列H(ωk)[=Hri(ωk)]を決定し;
1e3)各参照ポイント(r)に関し、n個の成分Ei(ωk,r)を、F(ωk,r)=H(ωk)×E(ωk,r)として決定し、ここで、E(ωk,r)=[Ei(ωk,r)]を、n個の成分Ei(ωk,r)を有したベクトルとし、F(ωk,r)=[Fl(ωk,r)]を、m個の成分Fl(ωk,r)を有したベクトルとし、lを、1〜mの中の整数とし、m個の成分Fl(ωk,r)を、ポイント(r)におけるかつ周波数(ωk)における所望のターゲット音響波場の生成に対応したものとする。
サブステップ(1e3)においては、例えば、少なくとも伝達行列H(ωk)の逆行列を求めることにより、p行列H−1(ωk)を計算することができ、脳内の各参照ポイント(r)に関して、以下の式を使用して、ベクトルE(ωk,r)を計算する。
E(ωk,r) = H−1(ωk)×F(ωk,r)
上記トレーニングステージにおいては、方法の使用に関係なく、脳の複数の参照ポイント(r)(18)と、撮影アレイ(17b)の各トランスデューサと、の間において(つまり、脳から、トランスデューサアレイ(17)に対しての、音響波の伝搬に関し)、インパルス応答hir(t)を決定することが好ましく、また、エコー像を形成し得るよう、各参照ポイント(18)において受領される撮影アレイ(17b)の焦点合わせ方法を決定する。
よって、脳(19)によって後方散乱された音響信号を逆重畳することにより、エコー像の撮影目的の装置(1)を使用することができる。この逆重畳は、任意の公知手段によって実行することができる。例えば、単一周波数成分(ωk)の数(p)に関して、フーリエ変換後に周波数ドメイン内で動作させることにより、また、伝達行列H(ωk)=[Hir(ωk)]の逆行列を求めることにより、実行することができる。その場合、成分Hir(ωk)は、周波数(ωk)におけるインパルス応答hir(t)のフーリエ形式の周波数成分とされる。これにより、脳の各ポイントにおいて後方散乱された信号を得ることができ、特に、これらの信号の振幅(例えば、各振幅の最大値)により、脳(19)の正確なエコー像を形成することができる。
シミュレーションによって、上記個々の音響信号が決定された後に、音響トランスデューサの実際のアレイ(17)を、デジタルシミュレーションの際に使用されたのと正確に同じ位置において、患者の頭蓋骨(14)の外部において、位置決めする。
複数のトランスデューサからなるアレイ(17)の位置が、デジタルシミュレーションの際に使用され位置と同じ位置であることを保証するために、トランスデューサアレイを付帯したキャップ(2)を、初期的に、頭蓋骨(14)に対してほぼ正確な位置に配置する。次に、頭蓋骨に対するキャップ(2)の本当の位置を、エコー撮影を行う撮影サブアレイ(17b)によって位置決めする。
この位置決めは、以下のようにして実行することができる。
−エコー検査を使用して、複数のトランスデューサからなるアレイと、定位固定フレーム(8)のうちの、キャップのフレキシブル壁(12)と接触している複数の部分と、の間の相対的位置決めを行うことにより、実行することができる、および/または、
−トランスデューサアレイと頭蓋骨(14)との間の相対的位置決めを直接的に行い、その後、エコー像によって、頭蓋骨(14)の少なくとも一部に関する外側形状を決定し、その後、頭蓋骨の外側形状と、頭蓋骨の3次元画像とを、比較することにより、実行することができる(後者の方法は、特に定位固定フレーム(8)を使用しない場合に、それ自体で行うことができる)。
この相対的位置決めから開始して、マイクロコンピュータ(7)は、ロボットアームによってキャップ(2)を移動させることができる。これにより、キャップの位置を、シミュレーションのために使用した位置と、正確に対応させることができる。
複数のトランスデューサからなるアレイ(17)を頭蓋骨(14)に対して位置決めした後に、施術者は、例えば画像撮影といった目的や例えば温熱療法といった目的に応じて、患者の脳(19)内において、任意に、ターゲット音響波場を形成することができる。
施術者が、脳の1つまたは複数の領域において温熱療法を行おうとした場合には、音響波を、被処理領域に対して、非常に正確に焦点合わせすることができる(複数の領域を被処理領域とすることができ、焦点合わせは、連続的にあるいは並列的に行うことができる)。
すべての場合において、温熱療法は、有利には、不連続な順次的な期間について行われる。これにより、破壊を目的とした領域の外部にまで傷害を引き起こしかねないような熱拡散をもたらしてしまうような長時間の加熱を、回避することができる。温熱療法の期間どうしの間においては、上述したように、脳の順次的なエコー像を撮影することができる。
エコー像は、有利には、マイクロコンピュータ(7)によって処理することができる。これにより、脳の温度マップを形成することができる。特に複合撮影と称される方法を使用すれば、温熱療法の有効性を可視化することができる。この方法は、特に、M. Pernot 氏他による“Improvement of ultrasound based temperature estimation by compound
imaging”, Proceedings of the International Symposium on Therapeutic Ultrasound,
Seattle 2002 という文献、および、Entrekin 氏他による“Real time spatial imaging compound in breast ultrasound: technology and early clinical experiment”,Medica
Mundi, Vol. 43, Sept. 1999, p. 35-43 という文献、に記載されている。
本発明の一実施形態において超音波を発生するための装置を全体的に示す斜視図である。 図1の装置におけるキャップ形成部材を、患者の頭部上に適用した状態で、詳細に示す部分断面図である。 図1の装置を示すブロック図である。
符号の説明
4 ロボットアーム(精密位置決め手段)
7 マイクロコンピュータ(マップ形成手段、シミュレート手段、計算手段、位置決め手段、精密位置決め手段)
8 定位固定フレーム(位置決めデバイス)
12 フレキシブル壁
13 流体充填タンク
14 頭蓋骨(骨境界壁)
16 ネジ(固定手段)
17 トランスデューサアレイ
17a 温熱療法のためのサブアレイ
17b 撮影のためのサブアレイ
19 実質的に均質な媒体

Claims (32)

  1. 骨境界壁(14)により囲まれた実質的に均質な媒体(19)内に、少なくとも1つの所定のターゲット音響波場を形成するための非侵襲的方法であって、
    少なくとも1つのトランスデューサアレイ(17)を具備した装置が、この方法を実施するものとされ、
    この方法においては、前記骨境界壁(14)を介して前記少なくとも1つのトランスデューサアレイ(17)が、音響信号を放出し、
    さらに、この方法においては、
    −トレーニングステージを行い、このトレーニングステージにおいては、
    1a)少なくとも部分的にX線を使用することによって、前記骨境界壁(14)の3次元画像を撮影し、これにより、様々なポイントにおける前記骨境界壁の空隙度合いを代理するパラメータを獲得し;
    1b)前記3次元画像に基づいて、前記骨境界壁内における少なくとも密度と音響波の速度と音響波の吸収とに関する3次元マップを決定し;
    1c)前記骨境界壁(14)に対しての前記トランスデューサアレイ(17)の特定位置を決定し;
    1d)前記少なくとも実質的に均質な媒体の少なくとも1つのポイントと、前記トランスデューサアレイのうちの少なくともいくつかのトランスデューサと、の間にわたっての音響波の少なくとも1つの伝搬を、伝搬に関する機械的モデルと、密度と音響波の速度と音響波の吸収とに関する前記3次元マップと、に基づいて、シミュレートし;
    1e)このシミュレーションに基づいて、前記ターゲット音響波場が得られるよう、前記トランスデューサアレイのうちの前記少なくともいくつかのトランスデューサから放出すべき個々の音響信号を計算し;
    −前記骨境界壁上における前記トランスデューサアレイの正確な位置決めを実際に行うというステージを行い、このステージにおいては、
    2a)前記骨境界壁(14)上において、前記トランスデューサアレイ(17)を、前記特定位置に対応するよう初期的に大まかに位置決めし;
    2b)前記トランスデューサアレイのうちの前記少なくともいくつかのトランスデューサを使用して、エコー像の撮影を行い、これにより、前記骨境界壁(14)に対しての前記トランスデューサアレイ(17)の位置を求め;
    2c)このステップ2b)において得られた位置に基づいて、前記骨境界壁(14)に対しての前記トランスデューサアレイ(17)の精細な位置決めを行う;
    ことを特徴とする方法。
  2. 請求項1記載の方法において、
    前記ステップ2b)の際には、
    前記骨境界壁(14)に対しての前記トランスデューサアレイ(17)の相対位置を、
    前記トランスデューサアレイのうちの少なくとも一部を使用してエコー像の撮影を行うことより前記骨境界壁(14)の少なくとも一部の外側形状を決定し;
    その後、得られた外側形状と、前記骨境界壁に関する前記3次元画像と、を比較する;
    ことにより求めることを特徴とする方法。
  3. 請求項1記載の方法において、
    −前記ステップ1a)の際には、前記骨境界壁(14)に対して位置決めデバイス(8)を堅固に固定し、この位置決めデバイス(8)を、X線の吸収に関して適合したものとし、これにより、初期的ステップを実施し;
    −前記ステップ1a)の際には、前記3次元画像を、前記骨境界壁(14)上の前記位置決めデバイス(8)の位置をも撮影したものとし;
    −前記ステップ2b)の際には、前記トランスデューサアレイ(17)と前記位置決めデバイス(8)との間の相対位置を、エコー像によって認識する;
    ことにより求めることを特徴とする方法。
  4. 請求項1〜3のいずれか1項に記載の方法において、
    前記トランスデューサアレイを、少なくとも1つのフレキシブル壁(12)を有した流体充填タンク(13)の中に設置し、
    前記フレキシブル壁を、前記骨境界壁(14)に対して押圧することを特徴とする方法。
  5. 骨境界壁(14)を介して少なくとも1つのトランスデューサアレイ(17)から放出された音響信号によって、脳を囲んでいる頭蓋骨の一部から前記骨境界壁(14)が形成されている場合に、脳のうちの、前記骨境界壁(14)により囲まれている少なくとも一部からなる実質的に均質な媒体(19)内に、少なくとも1つの所定のターゲット音響波場を形成するための非侵襲的方法であって、
    この方法においては、
    −トレーニングステージを行い、このトレーニングステージにおいては、
    1a)少なくとも部分的にX線を使用することによって、前記骨境界壁(14)の3次元画像を撮影し、これにより、様々なポイントにおける前記骨境界壁の空隙度合いを代理するパラメータを獲得し;
    1b)前記3次元画像に基づいて、前記骨境界壁内における少なくとも密度と音響波の速度と音響波の吸収とに関する3次元マップを決定し;
    1c)前記骨境界壁(14)に対しての前記トランスデューサアレイ(17)の特定位置を決定し;
    1d)前記少なくとも実質的に均質な媒体の少なくとも1つのポイントと、前記トランスデューサアレイのうちの少なくともいくつかのトランスデューサと、の間にわたっての音響波の少なくとも1つの伝搬を、伝搬に関する機械的モデルと、密度と音響波の速度と音響波の吸収とに関する前記3次元マップと、に基づいて、シミュレートし;
    1e)このシミュレーションに基づいて、前記ターゲット音響波場が得られるよう、前記トランスデューサアレイのうちの前記少なくともいくつかのトランスデューサから放出すべき個々の音響信号を計算し;
    −前記骨境界壁上における前記トランスデューサアレイの正確な位置決めを実際に行うというステージを行い、このステージにおいては、
    2a)前記骨境界壁(14)上において、前記トランスデューサアレイ(17)を、前記特定位置に対応するよう初期的に大まかに位置決めし;
    2b)前記トランスデューサアレイのうちの前記少なくともいくつかのトランスデューサを使用して、エコー像の撮影を行い、これにより、前記骨境界壁(14)に対しての前記トランスデューサアレイ(17)の位置を求め;
    2c)このステップ2b)において得られた位置に基づいて、前記骨境界壁(14)に対しての前記トランスデューサアレイ(17)の精細な位置決めを行う;
    ことを特徴とする方法。
  6. 請求項1〜5のいずれか1項に記載の方法において、
    前記音響波の周波数を、0.5MHz〜3MHzとすることを特徴とする方法。
  7. 請求項1〜6のいずれか1項に記載の方法において、
    前記ステップ1e)の次に、ステップ1f)を行い、
    このステップ1f)においては、
    前記トランスデューサアレイ(17)から放出される音響信号を、シミュレートし;
    前記信号を、所望の音響波場が得られるよう前記個々の音響信号に基づいて決定し;
    そのような放出によって形成される音響波の伝搬を、シミュレートし;
    前記伝搬が、ある種の所定基準を満たすかどうかを検証する;
    ことを特徴とする方法。
  8. 請求項1〜7のいずれか1項に記載の方法において、
    −前記ステップ1d)の際には、
    前記実質的に均質な媒体内の少なくとも1つのポイント(18)から、前記トランスデューサアレイのうちの少なくともいくつかのトランスデューサに向けての、音響波の伝搬をシミュレートし;
    iを1〜nの中のある整数としかつnを前記トランスデューサアレイ内における前記トランスデューサの総数としたときに、前記トランスデューサアレイ内における各トランスデューサ(i)の位置へと到達して受領されたシミュレーション音響信号Ri(t)を決定し;
    −前記ステップ1e)の際には、
    対象をなす各トランスデューサ(i)から放出すべき個々の音響信号Ei(t)を、前記ステップ1d)において決定された受領シミュレーション音響信号Ri(t)の時間反転Ri(−t)に比例するものとして、決定する。
    ことを特徴とする方法。
  9. 請求項8記載の方法において、
    前記ステップ1e)の際には、放出される音響信号Ei(t)を、Ei(t) = Gi×Ri(−t)によって決定し、
    ここで、Giは、ゲイン因子であって、トランスデューサ(i)ごとに異なるものとされ、これにより、前記骨境界壁内の損失を補償し得るものとされていることを特徴とする方法。
  10. 請求項9記載の方法において、
    少なくともいくつかのトランスデューサに対応した前記ゲイン因子(Gi)を、前記受領シミュレーション音響信号Ri(t)の振幅の2乗に反比例した係数としたものとすることを特徴とする方法。
  11. 請求項8記載の方法において、
    −前記ステップ1d)の際には、
    前記ステップ(1b)において決定された実際の吸収係数(τ)とは正反対であるような吸収係数(−τ)を前記骨境界壁の各ポイントにおいて有しているような音響波の吸収の関しての仮想的3次元マップを使用することによって、前記シミュレートを行い;
    −前記ステップ1e)の際には、放出すべき前記個々の音響信号Ei(t)を、前記時間反転Ri(−t)に等しいものとして決定する;
    ことを特徴とする方法。
  12. 請求項1〜7のいずれか1項に記載の方法において、
    前記トランスデューサアレイ(17)を、前記骨境界壁(14)に対して押圧することを意図した少なくとも1つのフレキシブル壁(12)を有した流体充填タンク(13)の中に設置し、前記ステップ1c)において決定した前記トランスデューサの前記特定位置を、前記骨境界壁に対して接触しないものとし、
    前記ステップ1d)の際には、
    −前記実質的に均質な媒体(19)の少なくとも1つのポイント(18)から、前記トランスデューサアレイのうちの少なくともいくつかのトランスデューサに向けての、音響波の伝搬を決定し;
    iを1〜nの中のある整数としかつnを前記トランスデューサアレイ内における前記トランスデューサの総数としたときに、前記トランスデューサアレイ内における各トランスデューサ(i)の位置へと到達して受領されたシミュレーション音響信号Ri(t)を決定し;
    −その後、各トランスデューサ(i)からの、信号Ri(t)の時間反転に対応した音響信号Ri(−t)の放出をシミュレートし、前記骨境界壁(14)に接触した各ポイントのところにおいて各トランスデューサ(i)に対応した仮想的トランスデューサ(i)に対しての、前記流体中の伝搬をシミュレートし、z仮想的トランスデューサ(i)の位置へと到達して受領されたシミュレーション音響信号R’i(t)を決定し;
    −その後、各仮想的トランスデューサ(i)からの、音響信号G’i×R’i(−t)の放出を、シミュレートし、ここで、R’i(−t)を、前記信号R’i(t)の時間反転とし、G’iを、少なくともいくつかの前記仮想的トランスデューサ(i)に関しての、前記信号R’i(−t)の振幅の2乗に反比例した係数とし、
    −その後、前記トランスデューサ(i)に向けての前記流体中の伝搬をシミュレートし、前記トランスデューサ(i)の位置へと到達して受領されたシミュレーション音響信号R”i(t)を決定し;
    前記ステップ1e)の際には、
    −放出されるべき個々の音響信号Ei(t)を、前記受領シミュレーション音響信号R”i(t)の時間反転R”i(−t)に等しいものとして、決定する;
    ことを特徴とする方法。
  13. 請求項1〜7のいずれか1項に記載の方法において、
    前記ステップ1d)の際には、
    −前記トランスデューサアレイ(17)内の少なくともいくつかのトランスデューサ(i)からの音響波パルスの放出をシミュレートし、ここで、iを、トランスデューサアレイ内のトランスデューサを表す指数とし、nを、トランスデューサの総数を表す正の整数とし、rを、1〜mの中のある整数とし、mを、参照ポイントの総数を表す正の整数としたときに、各トランスデューサ(i)から、前記実質的に均質な媒体(19)内に位置した複数の参照ポイント(r)に向けての、音響波の伝搬をシミュレートし;
    −前記実質的に均質な媒体内の前記各参照ポイント(r)に到達するものとしてシミュレートされたインパルス応答hri(t)を決定し;
    前記ステップ1e)においては、以下のサブステップを行い、
    1e1)それぞれ周波数(ωk)を有したシミュレートされた複数のインパルス応答の各々に関して、周波数成分の数(p)を決定し、ここで、kを、周波数成分を表す1〜pの間の指数とし;
    1e2)iを、1〜nの中の整数とし、rを、1〜mの中の整数とし、Hri(ωk)を、インパルス応答Hri(t)のフーリエ形式における周波数(ωk)における値としたときに、p伝達行列H(ωk)[=Hri(ωk)]を決定し;
    1e3)各参照ポイント(r)に関し、n個の成分Ei(ωk,r)を、F(ωk,r)=H(ωk)×E(ωk,r)として決定し、ここで、E(ωk,r)=[Ei(ωk,r)]を、n個の成分Ei(ωk,r)を有したベクトルとし、F(ωk,r)=[Fl(ωk,r)]を、m個の成分Fl(ωk,r)を有したベクトルとし、lを、1〜mの中の整数とし、m個の成分Fl(ωk,r)を、ポイント(r)におけるかつ周波数(ωk)における所望のターゲット音響波場の生成に対応したものとする;
    ことを特徴とする方法。
  14. 請求項13記載の方法において、
    前記サブステップ(1e3)においては、少なくとも伝達行列H(ωk)の逆行列を求めることにより、p行列H−1(ωk)を計算し、前記実質的に均質な媒体内の各参照ポイント(r)に関して、E(ωk,r) = H−1(ωk)×F(ωk,r)という式を使用して、ベクトルE(ωk,r)を計算することを特徴とする方法。
  15. 請求項1〜14のいずれか1項に記載の方法において、
    前記ステップ1d)の際には、
    −前記実質的に均質な媒体内の複数の参照ポイント(r)と、前記トランスデューサアレイ(17)内の少なくともいくつかのトランスデューサ(i)と、の間において、インパルス応答hri(t)を決定し、ここで、iを、1〜nの中のトランスデューサを表す指数とし、nを、トランスデューサの数を表す整数とし、rを、1〜mの中の整数とし、mを、参照ポイントの数を表す整数とし、
    前記ステップ1e)の際には、
    −エコー像を形成し得るよう、各参照ポイント(r)において受領されるトランスデューサアレイ(17)の少なくとも一部の焦点合わせ方法を決定する;
    ことを特徴とする方法。
  16. エコー像によって医療用の撮影を行うための方法であって、
    請求項1〜15のいずれか1項に記載されたようなターゲット音響波場を形成するための非侵襲的方法を行い、さらに撮影ステージを行い、
    この撮影ステージにおいては、
    −前記トレーニングステージの時点で決定したような前記個々の音響信号を使用することにより、前記トランスデューサアレイ(17)の少なくとも一部を使用して、前記実質的に均質な媒体(19)に関する少なくとも1つのエコー像を撮影する;
    ことを特徴とする方法。
  17. 請求項1〜16のいずれか1項に記載されたような方法を実施し得るよう構成された装置であって、
    −複数のトランスデューサからなるアレイ(17)であるとともに、実質的に均質な媒体(19)を囲む骨境界壁(14)の外側に配置されるものとされた、トランスデューサアレイ(17)と;
    −様々なポイントにおける前記骨境界壁の空隙度合いを与えるものとしてX線によって撮影された3次元画像に基づいて、前記骨境界壁内における少なくとも密度と音響波の速度と音響波の吸収とに関する3次元マップを決定するためのマップ形成手段(7)と;
    −前記実質的に均質な媒体の少なくとも1つのポイントと、前記トランスデューサアレイのうちの少なくともいくつかのトランスデューサと、の間にわたっての音響波の少なくとも1つの伝搬を、伝搬に関する機械的モデルと、密度と音響波の速度と音響波の吸収とに関する前記3次元マップと、に基づいて、前記骨境界壁(14)に対しての前記トランスデューサアレイ(17)の特定位置の関数として、シミュレートするシミュレート手段(7)と;
    −このシミュレーションに基づいて、前記実質的に均質な媒体内において前記ターゲット音響波場が得られるよう、前記トランスデューサアレイ(17)のうちの前記少なくともいくつかのトランスデューサから放出すべき個々の音響信号を計算するための計算手段(7)と;
    −前記トランスデューサアレイ(17)のうちの少なくともいくつかのトランスデューサを使用して、エコー像の撮影を行い、これにより、前記骨境界壁(14)に対しての前記トランスデューサアレイの位置決めを行うための、位置決め手段(7)と;
    −前記骨境界壁に対して前記トランスデューサアレイが前記特定位置に対応するようにして前記骨境界壁に対して配置されたときに、前記トランスデューサアレイの位置の関数として、前記骨境界壁(14)に対しての前記トランスデューサアレイ(17)の初期位置を、精密に位置決めするための精密位置決め手段(7,4)と;
    を具備していることを特徴とする装置。
  18. 請求項17記載の装置において、
    前記骨境界壁(14)に対しての前記トランスデューサアレイ(17)の位置決めを行うための前記位置決め手段(7)が、前記トランスデューサアレイのうちの少なくとも一部を使用してエコー像の撮影を行うことより前記骨境界壁(14)の少なくとも一部の外側形状を決定し、その後、得られた外側形状と、前記骨境界壁に関する前記3次元画像と、を比較することにより、位置決めを行うものとされていることを特徴とする装置。
  19. 請求項17または18記載の装置において、
    位置決めデバイス(8)と、この位置決めデバイス(8)を前記骨境界壁(14)上に堅固に固定し得る固定手段(16)と、を具備し、
    前記位置決めデバイスは、X線を吸収し得るものとされ、これにより、前記骨境界壁(14)の前記3次元画像内において可視のものとされ、
    前記骨境界壁(14)に対しての前記トランスデューサアレイ(17)の位置決めを行うための前記位置決め手段(7)が、エコー像の撮影によって、前記位置決めデバイス(8)に対しての前記トランスデューサアレイ(17)の位置を位置決めすることを特徴とする装置。
  20. 請求項17〜19のいずれか1項に記載の装置において、
    前記トランスデューサアレイ(17)が、前記骨境界壁(14)に対して押圧されることとなる少なくとも1つのフレキシブル壁(12)を有した流体充填タンク(13)の中に、設置されることを特徴とする装置。
  21. 請求項17〜20のいずれか1項に記載の装置において、
    前記トランスデューサアレイ(17)が、撮影のためのサブアレイ(17b)と、温熱療法のためのサブアレイ(17a)と、を備え、
    これら2つのサブアレイは、互いに異なるタイプのトランスデューサを有していることを特徴とする装置。
  22. 請求項17〜21のいずれか1項に記載の装置において、
    前記トランスデューサアレイ(17)が、0.5MHz〜3MHzという範囲の周波数の音響波を放出するものとされていることを特徴とする装置。
  23. 請求項17〜22のいずれか1項に記載の装置において、
    前記シミュレート手段(7)が、
    −所望の音響波場が得られるよう、少なくとも1つのトランスデューサアレイ(17)からの、前記個々の音響信号に基づいて決定された音響信号の放出をシミュレートすることができ;
    −そのような放出によって形成される音響波の伝搬を、シミュレートすることができ;
    −前記伝搬が、ある種の所定基準を満たすかどうかを検証することができる;
    ことを特徴とする装置。
  24. 請求項17〜23のいずれか1項に記載の装置において、
    前記シミュレート手段(7)が、前記実質的に均質な媒体内の少なくとも1つのポイント(18)から、前記トランスデューサアレイのうちの少なくともいくつかのトランスデューサに向けての、音響波の伝搬をシミュレートすることができ、これにより、iを1〜nの中のある整数としかつnを前記トランスデューサアレイ内における前記トランスデューサの総数としたときに、前記トランスデューサアレイ内における各トランスデューサ(i)の位置へと到達して受領されたシミュレーション音響信号Ri(t)を決定することができ;
    前記計算手段(7)が、対象をなす各トランスデューサ(i)から放出すべき個々の音響信号Ei(t)を、前記ステップ1d)において決定された受領シミュレーション音響信号Ri(t)の時間反転Ri(−t)に比例するものとして、決定することができる。;
    ことを特徴とする装置。
  25. 請求項24記載の装置において、
    前記計算手段(7)が、放出される音響信号Ei(t)を、Ei(t) = Gi×Ri(−t)によって決定し、ここで、Giは、ゲイン因子であって、トランスデューサ(i)ごとに異なるものとされ、これにより、前記骨境界壁内の損失を補償し得るものとされていることを特徴とする装置。
  26. 請求項25記載の装置において、
    少なくともいくつかのトランスデューサに対応した前記ゲイン因子(Gi)が、前記受領シミュレーション音響信号Ri(t)の振幅の2乗に反比例した係数とされていることを特徴とする装置。
  27. 請求項24記載の装置において、
    前記シミュレート手段(7)が、前記ステップ(1b)において決定された実際の吸収係数(τ)とは正反対であるような吸収係数(−τ)を前記骨境界壁の各ポイントにおいて有しているような音響波の吸収の関しての仮想的3次元マップを使用することによって、前記音響波伝搬に関する前記シミュレートを行うことができ;
    前記計算手段(7)が、放出すべき前記個々の音響信号Ei(t)を、前記時間反転Ri(−t)に等しいものとして決定することができる;
    ことを特徴とする装置。
  28. 請求項17〜23のいずれか1項に記載の装置において、
    前記トランスデューサアレイ(17)が、前記骨境界壁(14)に対して押圧することを意図した少なくとも1つのフレキシブル壁(12)を有した流体充填タンク(13)の中に設置去れ、
    前記シミュレート手段(7)によって考慮された前記トランスデューサの前記特定位置が、前記骨境界壁に対して接触しないものとされ、
    前記シミュレート手段(7)が、
    −前記実質的に均質な媒体(19)の少なくとも1つのポイント(18)から、前記トランスデューサアレイのうちの少なくともいくつかのトランスデューサに向けての、音響波の伝搬をシミュレートし得るとともに;
    iを1〜nの中のある整数としかつnを前記トランスデューサアレイ内における前記トランスデューサの総数としたときに、前記トランスデューサアレイ内における各トランスデューサ(i)の位置へと到達して受領されたシミュレーション音響信号Ri(t)を決定することができ;
    −各トランスデューサ(i)からの、信号Ri(t)の時間反転に対応した音響信号Ri(−t)の放出をシミュレートし得るとともに、前記骨境界壁(14)に接触した各ポイントのところにおいて各トランスデューサ(i)に対応した仮想的トランスデューサ(i)に対しての、前記流体中の伝搬をシミュレートすることができ、さらに、仮想的トランスデューサ(i)の位置へと到達して受領されたシミュレーション音響信号R’i(t)を決定することができ;
    −各仮想的トランスデューサ(i)からの、音響信号G’i×R’i(−t)の放出を、シミュレートすることができ、ここで、R’i(−t)は、前記信号R’i(t)の時間反転とされ、G’iは、少なくともいくつかの前記仮想的トランスデューサ(i)に関しての、前記信号R’i(−t)の振幅の2乗に反比例した係数とされ;
    −前記トランスデューサ(i)に向けての前記流体中の伝搬をシミュレートし得るとともに、前記トランスデューサ(i)の位置へと到達して受領されたシミュレーション音響信号R”i(t)を決定することができ;
    前記計算手段(7)が、
    −放出されるべき個々の音響信号Ei(t)を、前記受領シミュレーション音響信号R”i(t)の時間反転R”i(−t)に等しいものとして、決定することができる;
    ことを特徴とする装置。
  29. 請求項17〜23のいずれか1項に記載の装置において、
    前記シミュレート手段(7)が、
    −前記トランスデューサアレイ(17)内の少なくともいくつかのトランスデューサ(i)からの音響波パルスの放出をシミュレートし得るとともに、ここで、iが、トランスデューサアレイ内のトランスデューサを表す指数とされ、nが、トランスデューサの総数を表す正の整数とされ、rが、1〜mの中のある整数とされ、mが、参照ポイントの総数を表す正の整数とされたときに、各トランスデューサ(i)から、前記実質的に均質な媒体(19)内に位置した複数の参照ポイント(r)に向けての、音響波の伝搬をシミュレートすることができ;
    −前記実質的に均質な媒体内の前記各参照ポイント(r)に到達するものとしてシミュレートされたインパルス応答hri(t)を決定することができ;
    前記計算手段(7)が、
    −それぞれ周波数(ωk)を有したシミュレートされた複数のインパルス応答の各々に関して、周波数成分の数(p)を決定することができ、ここで、kは、周波数成分を表す1〜pの間の指数とされ;
    −iを、1〜nの中の整数とし、rを、1〜mの中の整数とし、Hri(ωk)を、インパルス応答Hri(t)のフーリエ形式における周波数(ωk)における値としたときに、p伝達行列H(ωk)[=Hri(ωk)]を決定することができ;
    −各参照ポイント(r)に関し、n個の成分Ei(ωk,r)を、F(ωk,r)=H(ωk)×E(ωk,r)として決定することができ、ここで、E(ωk,r)=[Ei(ωk,r)]が、n個の成分Ei(ωk,r)を有したベクトルとされ、F(ωk,r)=[Fl(ωk,r)]が、m個の成分Fl(ωk,r)を有したベクトルとされ、lが、1〜mの中の整数とされ、m個の成分Fl(ωk,r)が、ポイント(r)におけるかつ周波数(ωk)における所望のターゲット音響波場の生成に対応したものとされる;
    ことを特徴とする装置。
  30. 請求項29記載の装置において、
    前記計算手段(7)が、
    −少なくとも伝達行列H(ωk)の逆行列を求めることにより、p行列H−1(ωk)を計算することができ;
    −前記実質的に均質な媒体内の各参照ポイント(r)に関して、E(ωk,r) = H−1(ωk)×F(ωk,r)という式を使用して、ベクトルE(ωk,r)を計算することができる;
    ことを特徴とする装置。
  31. 請求項17〜30のいずれか1項に記載の装置において、
    前記シミュレート手段(7)が、
    −前記実質的に均質な媒体内の複数の参照ポイント(r)と、前記トランスデューサアレイ(17)内の少なくともいくつかのトランスデューサ(i)と、の間において、インパルス応答hri(t)を決定することができ、ここで、iが、1〜nの中のトランスデューサを表す指数とされ、nは、トランスデューサの数を表す整数とされ、rは、1〜mの中の整数とされ、mは、参照ポイントの数を表す整数とされ;
    前記計算手段(7)が、
    −エコー像を形成し得るよう、各参照ポイント(r)において受領されるトランスデューサアレイ(17)の少なくとも一部の焦点合わせ方法を決定することができる;
    ことを特徴とする装置。
  32. 請求項17〜31のいずれか1項に記載の装置において、
    前記トランスデューサアレイ(17)の少なくとも一部を使用し、かつ、前記計算手段(7)によって決定された前記個々の音響信号を使用して、前記実質的に均質な媒体(19)の少なくとも1つのエコー像を撮影し得るものとされた撮影手段を具備している
    ことを特徴とする装置。
JP2004532220A 2002-08-28 2003-08-20 骨を境界壁として囲まれた実質的に均質な媒体内において所定の音響波場を形成するための非侵襲的な方法および撮影方法ならびにそれら方法を実施するための装置 Expired - Lifetime JP4448921B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR0210682A FR2843874B1 (fr) 2002-08-28 2002-08-28 Procede non invasif pour obtenir un champ predetermine d'ondes acoustiques dans un milieu sensiblement homogene masque par une barriere osseuse, procede d'imagerie, et dispositif pour la mise en oeuvre de ces procedes
PCT/FR2003/002554 WO2004019784A1 (fr) 2002-08-28 2003-08-20 Procede non invasif pour obtenir un champ predetermine d'ondes acoustiques dans un milieu sensiblement homogene masque par une barriere osseuse, procede d'imagerie, et dispositif pour la mise en oeuvre de procedes

Publications (2)

Publication Number Publication Date
JP2005537071A JP2005537071A (ja) 2005-12-08
JP4448921B2 true JP4448921B2 (ja) 2010-04-14

Family

ID=31502958

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2004532220A Expired - Lifetime JP4448921B2 (ja) 2002-08-28 2003-08-20 骨を境界壁として囲まれた実質的に均質な媒体内において所定の音響波場を形成するための非侵襲的な方法および撮影方法ならびにそれら方法を実施するための装置

Country Status (7)

Country Link
US (1) US7837623B2 (ja)
EP (1) EP1531734B1 (ja)
JP (1) JP4448921B2 (ja)
AU (1) AU2003274268A1 (ja)
FR (1) FR2843874B1 (ja)
IL (1) IL167103A (ja)
WO (1) WO2004019784A1 (ja)

Families Citing this family (49)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8328420B2 (en) 2003-04-22 2012-12-11 Marcio Marc Abreu Apparatus and method for measuring biologic parameters
US10123732B2 (en) 2002-04-22 2018-11-13 Geelux Holdings, Ltd. Apparatus and method for measuring biologic parameters
MXPA04010518A (es) 2002-04-22 2005-07-14 Marc Aurielo Martins Ab Marcio Metodo y aparato para medir parametros biologicos.
US10227063B2 (en) 2004-02-26 2019-03-12 Geelux Holdings, Ltd. Method and apparatus for biological evaluation
US10687785B2 (en) 2005-05-12 2020-06-23 The Trustees Of Columbia Univeristy In The City Of New York System and method for electromechanical activation of arrhythmias
CA2944239C (en) 2005-10-24 2018-03-20 Marcio Marc Abreu Apparatus and method for measuring biologic parameters
FR2912817B1 (fr) * 2007-02-21 2009-05-22 Super Sonic Imagine Sa Procede d'optimisation de la focalisation d'ondes au travers d'un element introducteur d'aberations.
US8545405B2 (en) * 2008-04-23 2013-10-01 Therataxis, Llc Device, methods, and control for sonic guidance of molecules and other material utilizing time-reversal acoustics
CA2725453C (en) * 2008-04-30 2021-09-14 Milux Holding S.A. Brain stimulation
EP2310094B1 (en) 2008-07-14 2014-10-22 Arizona Board Regents For And On Behalf Of Arizona State University Devices for modulating cellular activity using ultrasound
WO2010014977A1 (en) 2008-08-01 2010-02-04 The Trustees Of Columbia University In The City Of New York Systems and methods for matching and imaging tissue characteristics
WO2010030819A1 (en) 2008-09-10 2010-03-18 The Trustees Of Columbia University In The City Of New York Systems and methods for opening a tissue
US8582397B2 (en) * 2009-01-06 2013-11-12 Therataxis, Llc Creating, directing and steering regions of intensity of wave propagation in inhomogeneous media
EP2480144B1 (en) 2009-09-21 2024-03-06 The Trustees of Columbia University in the City of New York Systems for opening of a tissue barrier
JP6061678B2 (ja) * 2009-11-04 2017-01-18 アリゾナ・ボード・オブ・リージェンツ・オン・ビハーフ・オブ・アリゾナ・ステイト・ユニバーシティーArizona Board of Regents on behalf of Arizona State University 脳調節インターフェース装置
US20110190668A1 (en) * 2010-02-03 2011-08-04 Mishelevich David J Ultrasound neuromodulation of the sphenopalatine ganglion
WO2011079177A1 (en) * 2009-12-22 2011-06-30 The Trustees Of Columbia University In The City Of New York A planning system for targeting tissue structures with ultrasound
US8997572B2 (en) 2011-02-11 2015-04-07 Washington University Multi-focus optical-resolution photoacoustic microscopy with ultrasonic array detection
WO2012162664A1 (en) 2011-05-26 2012-11-29 The Trustees Of Columbia University In The City Of New York Systems and methods for opening of a tissue barrier in primates
US9042201B2 (en) 2011-10-21 2015-05-26 Thync, Inc. Method and system for direct communication
KR20130054003A (ko) * 2011-11-16 2013-05-24 삼성전자주식회사 해부학적 특성을 기반으로 초음파 조사 계획을 수립하는 방법 및 장치
WO2014036170A1 (en) 2012-08-29 2014-03-06 Thync, Inc. Systems and devices for coupling ultrasound energy to a body
WO2014059170A1 (en) 2012-10-10 2014-04-17 The Trustees Of Columbia University In The City Of New York Systems and methods for mechanical mapping of cardiac rhythm
US11020006B2 (en) * 2012-10-18 2021-06-01 California Institute Of Technology Transcranial photoacoustic/thermoacoustic tomography brain imaging informed by adjunct image data
US9247921B2 (en) 2013-06-07 2016-02-02 The Trustees Of Columbia University In The City Of New York Systems and methods of high frame rate streaming for treatment monitoring
US10322178B2 (en) 2013-08-09 2019-06-18 The Trustees Of Columbia University In The City Of New York Systems and methods for targeted drug delivery
US10028723B2 (en) 2013-09-03 2018-07-24 The Trustees Of Columbia University In The City Of New York Systems and methods for real-time, transcranial monitoring of blood-brain barrier opening
CN105814419A (zh) 2013-10-11 2016-07-27 马尔西奥·马克·阿布雷乌 用于生物学评估的方法和设备
WO2015077355A1 (en) 2013-11-19 2015-05-28 Washington University Systems and methods of grueneisen-relaxation photoacoustic microscopy and photoacoustic wavefront shaping
EP3091895A4 (en) 2014-01-10 2017-08-23 Marcio Marc Abreu Device for measuring the infrared output of the abreu brain thermal tunnel
WO2015106180A1 (en) 2014-01-10 2015-07-16 Marcio Marc Abreu Devices to monitor and provide treatment at an abreu brain tunnel
JP2017503604A (ja) 2014-01-22 2017-02-02 マーシオ マーク アブリュー Abreu脳熱トンネルでの治療を提供するように構成されたデバイス
CN104548390B (zh) * 2014-12-26 2018-03-23 中国科学院深圳先进技术研究院 一种获得用于发射穿颅聚焦超声的超声发射序列的方法及系统
WO2016145215A1 (en) 2015-03-10 2016-09-15 Marcio Marc Abreu Devices, apparatuses, systems, and methods for measuring temperature of an abtt terminus
FR3038217B1 (fr) * 2015-07-01 2017-07-21 Centre Nat De La Rech Scient - Cnrs - Procede d'insonification pour obtenir un champ predetermine d'ondes ultrasonores, et procede de fabrication pour realiser une lentille ultrasonore a ces fins
EP3179313B1 (en) 2015-12-11 2021-11-10 Max-Planck-Gesellschaft zur Förderung der Wissenschaften e.V. Apparatus and method for creating a holographic ultrasound field in an object
US11672426B2 (en) 2017-05-10 2023-06-13 California Institute Of Technology Snapshot photoacoustic photography using an ergodic relay
JP6993886B2 (ja) * 2018-01-19 2022-01-14 富士フイルムヘルスケア株式会社 超音波送受信装置
EP3836831A4 (en) 2018-08-14 2022-05-18 California Institute of Technology MULTIFOCAL PHOTOACOUSTIC MICROSCOPY THROUGH AN ERGODIC RELAY
WO2020051246A1 (en) 2018-09-04 2020-03-12 California Institute Of Technology Enhanced-resolution infrared photoacoustic microscopy and spectroscopy
CA3128067A1 (en) 2019-02-13 2020-08-20 Alpheus Medical, Inc. Non-invasive sonodynamic therapy
US11369280B2 (en) 2019-03-01 2022-06-28 California Institute Of Technology Velocity-matched ultrasonic tagging in photoacoustic flowgraphy
WO2020193667A2 (en) * 2019-03-25 2020-10-01 Creaholic S.A. Treatment parameters for acoustic wave stimulation
US11654635B2 (en) 2019-04-18 2023-05-23 The Research Foundation For Suny Enhanced non-destructive testing in directed energy material processing
WO2021092250A1 (en) 2019-11-05 2021-05-14 California Institute Of Technology Spatiotemporal antialiasing in photoacoustic computed tomography
WO2021154730A1 (en) * 2020-01-27 2021-08-05 Cordance Medical Inc. Ultrasound transducer assembly
US20220087643A1 (en) * 2020-09-23 2022-03-24 3Dintegrated Aps Patient bearing system, a robotic system
KR102637479B1 (ko) * 2020-11-18 2024-02-16 주식회사 뉴로소나 저강도 집속 초음파 치료 장치
CN112915409A (zh) * 2021-03-01 2021-06-08 北京儒奥医疗科技有限公司 可智能靶向定位的超声球面相控治疗装置

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4869247A (en) 1988-03-11 1989-09-26 The University Of Virginia Alumni Patents Foundation Video tumor fighting system
FR2631707B1 (fr) * 1988-05-20 1991-11-29 Labo Electronique Physique Echographe ultrasonore a coherence de phase controlable
FR2683915A1 (fr) * 1991-11-18 1993-05-21 Philips Electronique Lab Appareil d'examen de milieux par echographie ultrasonore.
US6122341A (en) * 1992-06-12 2000-09-19 Butler; William E. System for determining target positions in the body observed in CT image data
US5855582A (en) * 1995-12-19 1999-01-05 Gildenberg; Philip L. Noninvasive stereotactic apparatus and method for relating data between medical devices
FR2791136B1 (fr) * 1999-03-15 2001-06-08 Mathias Fink Procede et dispositif d'imagerie utilisant les ondes de cisaillement
FR2815717B1 (fr) * 2000-10-20 2003-01-10 Centre Nat Rech Scient Procede et dispositif non invasif de focalisation d'ondes acoustiques

Also Published As

Publication number Publication date
FR2843874A1 (fr) 2004-03-05
WO2004019784A1 (fr) 2004-03-11
EP1531734A1 (fr) 2005-05-25
JP2005537071A (ja) 2005-12-08
EP1531734B1 (fr) 2014-07-02
IL167103A (en) 2011-07-31
US7837623B2 (en) 2010-11-23
AU2003274268A1 (en) 2004-03-19
US20050277824A1 (en) 2005-12-15
FR2843874B1 (fr) 2004-11-05
AU2003274268A8 (en) 2004-03-19

Similar Documents

Publication Publication Date Title
JP4448921B2 (ja) 骨を境界壁として囲まれた実質的に均質な媒体内において所定の音響波場を形成するための非侵襲的な方法および撮影方法ならびにそれら方法を実施するための装置
Marsac et al. Ex vivo optimisation of a heterogeneous speed of sound model of the human skull for non-invasive transcranial focused ultrasound at 1 MHz
US10973461B2 (en) Mapping of intra-body cavity using a distributed ultrasound array on basket catheter
CN110248606B (zh) 气穴定位
US7175599B2 (en) Shear mode diagnostic ultrasound
JP5530685B2 (ja) 剛性が変化した領域を検出するためのシステム及び方法
US8155725B2 (en) Method for optimising the focussing of waves through an aberration-inducing element
JP4629034B2 (ja) 剪断モード治療用超音波
JP3843256B2 (ja) 音波を集束するための方法と非侵襲性デバイス
CN108135565A (zh) 用于配准使用各种成像模态获得的图像并验证图像配准的系统和方法
Arthur et al. Temperature dependence of ultrasonic backscattered energy in motion compensated images
KR20140102995A (ko) 의료 영상을 이용하여 관심 영역 내에 다중 초점들을 형성하는 초음파를 생성하는 방법, 장치 및 hifu 시스템
US20210196229A1 (en) Real-time ultrasound monitoring for ablation therapy
Mougenot et al. High intensity focused ultrasound with large aperture transducers: a MRI based focal point correction for tissue heterogeneity
JP6853197B2 (ja) 超音波の所定の場を得るための照射法およびそのための超音波レンズの製造方法
US10765892B1 (en) Systems and methods for optimizing transcranial ultrasound focusing
Soler López et al. Application of ultrasound in medicine Part II: the ultrasonic transducer and its associated electronics
Seip Feedback for ultrasound thermotherapy
Owen et al. In vivo evaluation of a mechanically oscillating dual-mode applicator for ultrasound imaging and thermal ablation
Amanatiadis et al. Transcranial ultrasonic propagation and enhanced brain imaging exploiting the focusing effect of the skull
Gentry et al. Finite-element analysis of temperature rise and lesion formation from catheter ultrasound ablation transducers
Ebbini et al. Temperature imaging using diagnostic ultrasound: methods for guidance and monitoring of thermal treatments of tissue
Gentry et al. Finite element analysis of temperature rise from an integrated 3-D intracardiac echo and ultrasound ablation transducer
Pedrosa Exploring tissue characterization by temperature induced changes on ultrasound backscattered energy
Santo Pedrosa Exploring Tissue Characterization by Temperature Induced Changes on Ultrasound Backscattered Energy

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20060529

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090421

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20090701

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20091117

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20091216

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20091217

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20091217

R150 Certificate of patent or registration of utility model

Ref document number: 4448921

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130205

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140205

Year of fee payment: 4

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S633 Written request for registration of reclamation of name

Free format text: JAPANESE INTERMEDIATE CODE: R313633

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

EXPY Cancellation because of completion of term