Patent Application
20220126120
Embodiments of the present invention relate to systems and methods for ultrasound focusing. In particular, various embodiments are directed to efficient methods of focusing a phased array of ultrasound transducer elements, using preparatory measurements to adjust the relative phases of the transducer elements.
Physics of Ultrasound Transmission: Ultrasounds are mechanical waves (with a frequency of more than 20 kHz) that travel with alternating compression and rarefaction, thereby transmitting energy by molecular movements. In contrast to diagnostic ultrasonography (with a frequency in the megahertz range), the frequencies used for the transcranial FUS are either mid-frequency (e.g., around 650 kHz) or low frequency (e.g., around 220 kHz). The speed of ultrasound transmission is medium dependent (eg, water and soft tissue are excellent conductors, but air and bone are not). Transmission necessitates a coupling medium (eg, degassed water) between the transducer and biological tissue. In addition, there is significant absorption and reflection of ultrasound while traveling through tissues with different “ultrasonic densities,” such as layers of the skull (inner table, marrow, and outer table). The extent of ultrasound reflection is also dependent on the incident angle (eg, large (25° to 30°) incident angles result in higher reflection). Besides altering the velocity, heterogeneous tissues also alter the phase of ultrasound waves.
The biological consequences of FUS are dependent on the intensity, frequency, and duration of exposure. At low intensity, high-frequency (2 MHz) ultrasound can produce a reversible conduction block in peripheral nerves. This conduction block is associated with a mild increase in local temperature (41° C. to 45° C.) and is mediated by the inactivation of sodium channels. It can produce transient clinical results that may last for a few minutes and are particularly appealing for target localization in functional neurosurgery. Under certain conditions, FUS can also reversibly open the blood-brain barrier without ablation. Blood-brain barrier opening can be reliably achieved at subthreshold intensities with the use of microbubbles. A recent proof-of-concept study demonstrated localized blood-brain barrier openings of approximately 1 cm3 with very low sonication power (5 W and 230 kHz transducer).
At high intensity, FUS creates tissue ablation, the mechanism of which is dependent on frequency. For example, the mid-frequency system (650 kHz) primarily produces thermal ablation, whereas the low-frequency system (220 kHz) achieves ablation via cavitation or histotripsy, in which ultrasound interacts with trapped gas bubbles within tissues that leads to the rapid oscillation and collapse of those bubbles. Therapeutic sonications (temperature greater than 55° C.) denature cellular protein and produce lesions with 3 separate zones on T2-weighted MRI: mixed-intensity core (zone 1 with necrotic center), hyperintense periphery (zone 2 with apoptosis), and surrounding vasogenic edema (zone 3). The histological consequences of FUS ablation are associated with the duration of exposure (thermal dose). The target tissue characteristics influence temperature elevation during sonication (eg, tissue perfusion may act as a heat sink, although this reaction is minimal for short-duration sonications).
As recently reviewed (Vibhor Krishna et al, JAMA Neurol. 2018; 75(2):246-254), the ability to transiently open blood-brain barrier in localized brain regions offers a promising route for targeted drug delivery (eg, chemotherapy for neuro-oncology; gene and neurotrophic factors for neurodegenerative diseases. Moreover, targeted drug delivery can be an attractive option for intravascular thrombolysis (eg, in distal intracranial thrombosis where endovascular therapy has met with challenges).
The nonthermal qualities of FUS are attractive for neuromodulation. Successful cortical stimulation using transcranial FUS has been reported in mice (at 650 kHz), patients with chronic pain (at 8 MHz), and healthy volunteers (at 500 kHz). Recently was reported the case of a comatose patient with brain injury who received FUS stimulation to the entire thalamus. The patient showed some clinical improvement without major adverse effects. Finally, subthreshold sonications can be a screening tool for target selection for therapeutic interventions, especially in situations where either conventional screening (eg, levodopa challenge) provides insufficient answers (eg, subthalamic nucleus vs VIM for patients with a dual diagnosis of PD and essential tremor) or the most efficacious target for neuromodulation is unclear (eg, subthalamic nucleus vs globus pallidus internus for patients with PD).
Thermal ablation, as may be accomplished using focused ultrasound, has particular appeal for treating tissue within the brain and other tissue regions deep within the body, because it generally does not disturb intervening or surrounding healthy tissue. Focused ultrasound may also be attractive, because acoustic energy generally penetrates well through soft tissues, and ultrasonic energy, in particular, may be focused towards focal zones having a cross-section of only a few millimeters due to relatively short wavelengths (e.g., as small as 1.5 millimeters (mm) in cross-section at one Megahertz (1 MHz)). Thus, ultrasonic energy may be focused at a region deep within the body, such as a cancerous tumor or other diseased tissue, to ablate the diseased tissue without significantly damaging surrounding healthy tissue.
To focus ultrasonic energy towards a desired target, a piezoelectric transducer may be used that includes a plurality of transducer elements. A controller may provide drive signals to each of the transducer elements, thereby causing the transducer elements to transmit acoustic energy such that constructive interference occurs at a “focal zone”. The focal zone is typically defined as the region of intensity higher than half maximum, and is commonly characterized by a “peak width” in a given direction. The peak width may be anisotropic. In fact, most realized instrumental systems produce an elliptical shaped peak cross-section at half maximum. At the focal zone, sufficient acoustic energy may be delivered to generate the desired tissue activation (e.g., heating, necrosis, neural stimulation, etc. . . . ) within the focal zone and for a sufficient period until tissue affects occurs. Preferably, tissue along the path through which the acoustic energy passes (“the pass zone”) outside the focal zone, is affected (e.g, heated) only minimally, if at all, thereby minimizing damaging tissue outside the focal zone.
Phased arrays of ultrasound transducers are well-known as a system for focusing ultrasound energy at target sites inside the body. Constructive and destructive interference of acoustic waves transmitted by multiple transducers can be used to deliver complex spatiotemporal patterns of acoustic waves. Generally, phased arrays use tens to hundreds or even thousands of ultrasound transducers distributed spatially on the surface of the body. For instance, a phased array placed on the head can be used to target an area deep in the brain. However, phased arrays have important limitations for delivering ultrasound transcranially for neuromodulation. Phased arrays use spatially distributed transducers, requiring a larger form factor. Moreover, large and generally unportable power and control components are required to manage the timing, intensity, phase, and other properties of the ultrasound waves transmitted by each of the transducers.
The known art of ultrasound focusing onto tissue in general and brain tissue in particular is exemplified in U.S. Pat. Nos. 8,932,237, 5,329,930, 4,817,614, 8,088,067, 6,128,958, 7,611,462, 5,984,881, and references therein, the entire disclosures of which are hereby incorporated by reference.
Focused ultrasound (i.e., acoustic waves having a frequency greater than about 20 kilohertz) can be used to image or therapeutically treat internal body tissues within a patient. For example, ultrasonic waves may be used to ablate tumors, eliminating the need for the patient to undergo invasive surgery. For this purpose, a piezo-ceramic transducer is placed externally to the patient, but in close proximity to the tissue to be ablated (the “target”). The transducer converts an electronic drive signal into mechanical vibrations, resulting in the emission of acoustic waves (a process hereinafter referred to as “sonication”). The transducer may be shaped so that the waves converge in a focal zone. Alternatively or additionally, the transducer may be formed of a plurality of individually driven transducer elements whose phases (and, optionally, amplitudes) can each be controlled independently from one another and, thus, can be set so as to result in constructive interference of the individual acoustic waves in the focal zone. Such a “phased-array” transducer facilitates steering the focal zone to different locations by adjusting the relative phases between the transducers, and generally provides the higher a focus quality and resolution, the greater the number of transducer elements. Magnetic resonance imaging (MRI) may be utilized to visualize the focus and target in order to guide the ultrasound beam.
While the transducer is located external to the patient, it must be in direct contact and tightly coupled with a media that efficiently transmits the high frequency ultrasound waves. For example, the transducer can be positioned in a liquid bath that is capable of efficient transmission of the ultrasound waves. The patient's body must also be wetted and tightly coupled to the transmission media in order to ensure an optimal acoustic wave transmission path from the transducer to the focal zone.
While system parameters are generally fixed for a given transducer array, tissue homogeneity may vary significantly from patient to patient, and even between different tissue regions within the same patient. Tissue inhomogeneity may decrease intensity of the acoustic energy at the focal zone and may even move the location of the focal zone within the patient's body. Specifically, because the speed of sound differs in different types of tissue, as portions of a beam of acoustic energy travel along different paths towards the focal zone, they may experience a relative phase shift or time delay, which may change the intensity at the focal zone and/or move the location of the focal zone.
For example, the speed of sound through fat is approximately 1460 meters per second (m/s), while the speed of sound through muscle is approximately 1600 meters per second (m/s). The speed of sound through bone tissue is much faster, for example, approximately 3000 meters per second (m/s) for skull bone tissue. The speed of sound also varies in different organs. For example, the speed of sound in brain tissue is approximately 1570 meters per second (m/s), approximately 1555 meters per second (m/s) in the liver, and approximately 1565 meters per second (m/s) in the kidney.
The relative phases (alternatively relative “time shift”) at which the transducer elements need to be driven to result in a focus at the target location depend on the relative location and orientation of the transducer surface and the target, as well as on the dimensions and acoustic material properties (e.g., sound velocities) of the tissue or tissues between them (i.e., the “target tissue”). Thus, to the extent the geometry and acoustic material properties are known, the relative phases (and, optionally, amplitudes) can be calculated, as described, for example, in U.S. Pat. Nos. 6,612,988, 6,770,031, and 7,344,509, the entire disclosures of which are hereby incorporated by reference. In practice, however, knowledge of these parameters is often too incomplete or imprecise to enable high-quality focusing based on computations of the relative phases alone. For example, when ultrasound is focused into the brain to treat a tumor, the skull in the acoustic path may cause aberrations that are not readily ascertainable. In such situations, treatment is typically preceded by an auto-focusing procedure in which, iteratively, an ultrasound focus is generated at or near the target, the quality of the focus is measured (using, e.g., thermal imaging or acoustic radiation force imaging (ARFI)), and experimental feedback is used to adjust the phases of the transducer elements to achieve sufficient focus quality.
The auto-focusing procedure may thus take a substantial amount of time, which may render it impracticable or, at the least, inconvenient for a patient. While the effect of pre-therapeutic sonications may be minimized by employing an imaging technique that requires only low acoustic intensity (e.g., ARFI), it is generally desirable to limit the number of sonications prior to treatment. Accordingly, there is a need for more efficient ways of focusing a phased array of transducer element to create a high-quality ultrasound focus.
Another common technique for focusing ultrasound is by using a shaped lens with an acoustic velocity (i.e. speed of sound) that differs from adjoining air, tissue, or material to bend acoustic waves. Most standard ultrasound focusing lenses employ a single concave lens. However, a single concave lens focusing system for ultrasound has limitations. Ultrasound lenses comprised of a single concave lens are limited with regard to the range of focal lengths that can be achieved with a lens of a particular cross-sectional area. Short focal lengths cannot be achieved with smaller cross-sectional areas appropriate for systems affixed to the head or skull. Neuromodulation of superficial brain regions with an appropriate transcranial ultrasound system would be advantageous due to the importance of such superficial brain regions (e.g. cerebral cortex) to sensory, motor, higher cognitive function, and other brain functions.
To affect brain function transcranial ultrasound neuromodulation requires appropriate ultrasound waveform parameters, including acoustic frequencies generally less than about 10 MHz, spatial-peak temporal-average intensity generally less than about 10 W/cm2 (e.g., between 0.5 and 10 W/cm2), and appropriate pulsing and other waveform characteristics to ensure that heating of a targeted brain region does not exceed about 2 degrees Celsius for more than about 5 seconds. Transcranial ultrasound neuromodulation induces neuromodulation primarily through vibrational or mechanical mechanisms. Noninvasive and nondestructive transcranial ultrasound neuromodulation is in contrast to other transcranial ultrasound based techniques that use a combination of parameters to disrupt, damage, destroy, or otherwise affect neuronal cell populations so that they do not function properly and/or cause heating to damage or ablate tissue.
As by in the article Lindsey (Lindsey B D, Smith S W. Refraction Correction in 3D Transcranial Ultrasound Imaging. Ultrasonic imaging. 2014; 36(1):35-54. doi:10.1177/0161734613510287), Image quality in transcranial ultrasound remains limited by the deleterious effects of the presence of the skull, including attenuation, aberration, refraction, and mode conversion. The effects of aberration is induced by spatially inhomogeneous layers having a different longitudinal wave velocity from that typically assumed by the ultrasound scanner (c=1540 m/s). An anatomical sources of aberration include layers of bone in the skull (c≈2800 m/s, commonly within 15% variation due to difference in bone porosity and thickness) or layers of fat (c≈1450).
As discussed in the journal article by Ding et. al., Phys. Med. Biol. 60 (2015) 3975-3998, ultrasound penetration through the skull is better, with less energy deposition within the skull bone otself, at around 0.5 MHz compared with higher frequencies.
For the present disclosure, (i) the term outer-facing skull surface is used interchangeably with skull outer surface (or outer skull surface)—the abbreviation SOS may be used and (i) the term inner-facing skull surface is used interchangeably with skull inner surface (or inner skull surface—the abbreviation SIS may be used.
As shown in
Preferably, the delivered energy has as ‘tight’ of a focus as possible—e.g. most energy is focused at target location 900 within the peak area half-width half-maximum (HWHM).
Unfortunately, as shown in
Some embodiments of the invention relate to providing “feature A’ and/or “feature B” and/or “feature C” and/or “feature D” and/or “feature E” and/or “feature F” and/or “feature G” and/or “feature H” and/or “feature I,” described below in the present section.
Embodiments of the present invention relate to systems and methods for ultrasound focusing. In particular, various embodiments are directed to efficient methods of focusing a phased array of ultrasound transducer elements, using preparatory measurements to adjust the relative phases of the transducer elements.
For the present application, when material (e.g. the skull or brain) is ‘irradiated’ with ultrasound, this means that ultrasound is delivered to the material (e.g. the skull or brain).
Embodiments of the invention relate to apparatus and methods for generating an ultrasound intensity-peak within a human subject brain around a target peak location in the brain. In particular and not wishing to be bound by theory, embodiments of the present invention overcome any defocusing attributable to: (i) differences in the speed of sound between the skull and underlying brain tissue; (ii) the fact that the skull is of non-uniform thickness; and (iii) the fact that the speed of sound within the skull is not uniform.
Some embodiments are disclosed and claimed in terms of generating an ultrasound intensity-peak within a human subject brain. However, in general terms the presently disclosed teaching may be employed to produce an ultrasound intensity-peak within a generally round object (e.g. having a representative radius curvature of at (i) at least 7 cm or at least 10 cm and/or (ii) at most 50 cm or at most 40 cm or at most 25 cm)(of having: (i) an outer crust (e.g. including but not only skull/bone) having a lower liquid content (e.g. of non-uniform thickness) and (ii) underlying material (e.g. including brain but not only brain) having a greater liquid content. In this case, the ultrasound is delivery from transducers outside of the crust, through the crust, and into the underlying material so as to produce the ultrasound intensity-peak within the underlying material.
In order to counteract defocusing derived from the properties crust (e.g. skull) and/or the underlying material (e.g. brain), transducers disposed outside of the head (e.g. outside of the generally round object) are operated at phase-shifts that are computed by analyzing ultrasound measurement data acquired during a previously-performed calibration stage.
During this calibration stage, the skull (or any other crust) is subjected to plurality of ultrasound skull-probe events (USP) (or an ultrasound crust-probe UCP event for the more general case) where during each USP (or UCP) ultrasound of a respective ultrasound test signal (i.e. produced by ultrasound transducer(s) disposed over the skull or crust).
Instead of only probing a single location on the skull, in order to accurately compute phase-shift data it may be advantageous to probe the skull at a variety of locations on the skull (or other crust), and not just at a single location (Feature A). Furthermore, it is noted that the skull is not a flat plane—the skull is a generally round object, and the local surface orientation of the skull (i.e. the vector perpendicular to a local plane of any given location on the skull or other crust).
In some embodiments, not only is it desired to probe (i.e. with ultrasound of a calibration measurement) the skull (or other crust) at a variety of locations (Feature A) which may be distancesd from each other (Feature B), but this entails probing the skull or a variety of different local skull-surface orientations (or orientations of any other crust) (Feature C).
For non-binding theoretical reasons explained below, when probing the skull with an ultrasound test signal emitted by one or more of the transducers, it may be advantageous to do so such that an energy flux vector of the ultrasound test signal is substantially perpendicularly incident on the skull (Feature D).
Another feature provided by embodiments of the invention (Feature E) is that ultrasound transducers (i.e. those which participate in the calibration stage as calibration emitters) are held stable for the entirety of the calibration stage.
Thus, in some embodiments, in order to provide Feature D while also providing previous Features A-C and E, it may be useful to:
(I) deliver ultrasound for USP (or UCP event) of the calibration stage so that a majority of power of the respective ultrasound test signal is provided from transducers in a relatively small (e.g. a sphere having a radius at most 0.5 cm locale) locale (Feature F); and/or
(ii) arranged so that these locales are relatively close to the skull (Feature G)—e.g. a center of the 0.5-cm radius spherical locale in which the transducer(s) that provide at least 50% of power of the respective ultrasound signal for the given USP (or UCP) even is displaced from the skull 150 (or other crust) by at most 1.5 or at most 1 cm.
The ultrasound signal of each USP (or UCP event) is received by one or more ultrasound transducers. For each USP (or USP event), a difference in echo times (discussed below) is computed—for example, (i) a time where a signal reflected from the outer surface 151 of skull 150 may be the ‘first echo time’ and (ii) a time where a signal reflected from the inner surface 151 of skull 150 may be the ‘second echo time.’ The second echo time may be longer than the first echo time.
In some embodiments of the invention, when measuring the echo time-differences, it is possible to (Feature H) to “re-use” an ultrasound transducer (or a location thereof) which supplies the ultrasound test signal of an USP (or UCP event) by receiving reflected ultrasound back into the same ultrasound transducer or back into another ultrasound transducer displaced therefrom by no more than 0.5 cm. This may be carried out for at least one or at least some or at least a majority of the USP (or UCP events). In some examples, this may increase the accuracy and facilitate computation of face-differences that are useful for achieving a better focus.
In some embodiments of the invention, when producing the ultrasound-intensity peak within the brain beneath the skull (or other material beneath the crust), it may be possible to (Feature I) to “re-use” an ultrasound transducer (or a location thereof) which supplies the ultrasound test signal of an USP (or UCP event) during calibration. Thus, when operating transducers according using the computed phase-shifts to produce this ultrasound-intensity peak, at least some of the transducers (or some of the locales) deliver may be same as those used to probe the skull during the USP (or UCP) event of calibration.
Thus, embodiments of the present invention relate to delivery of ultrasound (e.g. non-surgically—i.e. non-invasively).
A method for generating an ultrasound intensity-peak within a human subject brain around a target-peak-location, by delivering ultrasound through the skull, the method comprising:
In some embodiments, the electronic processing comprises for each USP-eventi having its respective ultrasound test signal UTSi and its respective dominant emission-locale DEL(USP-eventi), computing a respective measured echo time-difference ETD(USP-eventi) between:
A method for generating an ultrasound intensity-peak within a human subject brain around a target-peak-location 900, by delivering ultrasound through the skull, the method comprising:
In some embodiments, the respective skull-surface locations max_intensity_SOS_LOC(USP-eventi) on the skull outer-facing surface form a location set SOS_LOC_SET {max_intensity_SOS_LOC((USP-event1), max_intensity_SOS_LOC((USP-event2) . . . max_intensity_SOS_LOC((USP-eventL)}, and wherein the skull-surface locations of the location set SOS_LOC_SET are distributed in space so that no two skull-surface locations (max_intensity_SOS_LOC(USP-eventj), max_intensity_SOS_LOC(USP-events)) on the skull outer-facing surface are displaced from each other by less than disp_numb cm, wherein disp_numb is a positive number whose value is at least 1 or at least 1.5 or at least 2.
In some embodiments, for each event USP-eventi of the L events of {USP-event1, USP-event2 . . . USP-eventL}, a center of the respective dominant emission-locale DEL(USP-eventi) thereof is displaced from the skull by at most 2.5 cm.
In some embodiments, for each event USP-eventi of the L events of {USP-event1, USP-event2 . . . USP-eventL}, a center of the respective dominant emission-locale DEL(USP-eventi) thereof is displaced from the skull by at most 2 cm or at most 1.5 cm.
In some embodiments, a nearest distance between (i) the target-peak-location 900 in brain tissue 98 beneath the skull and (ii) the skull is at least 2 cm or at least 4 cm or at least 5 cm.
In some embodiments, wherein each USP-eventi is defined by a different respective dominant emission-locale DEL(USP-eventi) such that during each USP-eventi the at least 50% of power of the respective ultrasound test signal UTSi received at the respective skull-surface location max_intensity_SOS_LOC(USP-eventi) on the skull is supplied by the transmitter(s) of the AUT (i) whose center(s) is(are) each disposed within a dominant emission-locale DEL(USP-eventi) that is spherical in shape with a radius of at most 0.5 cm; and (ii) which each have a width of at most ww cm, wherein 0≤ww≤1.
In some embodiments, wherein 0≤ww≤0.75.
In some embodiments, wherein 0≤ww≤0.5.
In some embodiments, wherein L≥7 or L≥8 or L≥10.
In some embodiments, wherein for each event USP-eventi of the L events of {USP-event1, USP-event2 . . . USP-eventL}, the respective one of more transducer(s) which:
In some embodiments, wherein the dominant emission locales are distributed in space so that a center-center distance between a first DEL(USP-events) and a second DEL(USP-event) of the dominant emission locales (j≠k) (i.e. for two different events) is at least 4 cm or at least 5 cm.
In some embodiments, at least 30% of ultrasound power of the intensity peak are supplied by transducer(s) disposed with any of the dominant emission-locales {DEL(USP-event1), DEL(USP-event2) . . . DEL(USP-eventL)} defined by the at least L ultrasound skull-probe (USP) events {USP-event1, USP-event2 . . . USP-eventL}
In some embodiments, wherein at least 50% (or at least 75%) of ultrasound power of the intensity peak are supplied by transducer(s) disposed with any of the dominant emission-locales {DEL(USP-event1), DEL(USP-event2) . . . DEL(USP-eventL)} defined by the at least L ultrasound skull-probe (USP) events {USP-event1, USP-event2 . . . USP-eventL}
In some embodiments, the ultrasound intensity-peak produced by operating the at least some transducers in the relative phases produces the intensity-peak so that a FWHM full width half maximum thereof in a widest dimension (e.g. of an oval-shaped peak) is at most 5 cm or at most 4 cm.
In some embodiments, the ultrasound intensity-peak produced by operating the at least some transducers in the relative phases produces the intensity-peak so that a FWHM full width half maximum thereof in a most narrow dimension (e.g. of an oval-shaped peak) is at most 2.5 cm or at most 2 cm or at most 1.5 cm or at most 1 cm.
In some embodiments, at least a majority of or at least 75% of or all ultrasound transducers employed to transmit ultrasound during calibration mode remain stationary during a time period which begins upon commencement of the operation in the calibration mode and concludes upon generation of the ultrasound intensity-peak.
In some embodiments, at least a majority of or at least 75% of or all ultrasound transducers operated according to the relative phases to form the intensity-peak remain stationary during a time period which begins upon commencement of the operation in the calibration mode and concludes upon generation of the ultrasound intensity-peak.
A system for generating an ultrasound intensity-peak within a human subject brain around a target-peak-location, by delivering ultrasound through the skull, the method comprises:
A system for generating an ultrasound intensity-peak within a human subject brain around a target-peak-location, by delivering ultrasound through the skull, the method comprises:
In some embodiments, each USP-eventi respectively defines an event-specific direction of an energy flux vector EFV(USP-eventi) describing the directional energy flux of the USP signal incident upon the skull outer surface to thereby defined an EVF set {EVF(USP-event1), EVF(USP-event2) . . . EVF(USP-eventL)} such that all EVF vectors of a sub-set of the EVC set are rotated from each other by at least 15 degrees, a cardinality of the sub-set being at least two.
42. The method or system of claim 41 wherein the cardinality is at least 3 or at least 5.
43. The method or system of any of claims 41-42 wherein each of {EVF(USP-eventi), EVF(USP-event2) . . . EVF(USP-eventL)} is substantially aligned with a local orientation of max_intensity_SOS_LOC(USP-eventi) within at most 10 degrees or at most 5 degrees.
For the present disclosure, any “circuitry” or “module” or “unit” may be implemented using any combination of hardware (e.g. a general purpose CPU or any other analog or digital hardware) and/or software and/or firmware.
A mechanical arrangement may include, for example, a device housing or other components that the skilled artisan would employ after reading the present disclosure.
A presently preferred embodiment of the invention will be described in detail, in conjunction with the accompanying drawings, in which:
The invention is herein described, by way of example only, with reference to the accompanying drawings. With specific reference now to the drawings in detail, it is stressed that the particulars shown are by way of example and for purposes of illustrative discussion of the embodiments of the exemplary system only and are presented in the cause of providing what is believed to be a useful and readily understood description of the principles and conceptual aspects of the invention. In this regard, no attempt is made to show structural details of the invention in more detail than is necessary for a fundamental understanding of the invention, the description taken with the drawings making apparent to those skilled in the art how several forms of the invention may be embodied in practice and how to make and use the embodiments.
It is possible to emit an ultrasound signal in vicinity of the skull 150 using one or more transducers 1550. For example, the transducer 1550 emits an ultrasound test signal which propagates in space, reaches skull 150, and echoes from the skull. The echoed ultrasound may be subsequently detected by the same transducer or by a different transducer.
Not shown in
As shown in
The time delay (tE−tD) (i.e. the delay between receiving of reflected ultrasound in
(i) a normal vector of surface 151 where ultrasound is incident thereon at time tB;
(ii) a directional energy flux vector (analogous to a Poynting vector from electricity and magnetism) of the ultrasound test signal when is incident thereon at time tB.
In the example of
Reference is made to
Specifically, in order to compute the phase-differences, a calibration measurement procedure P1110 and a calibration analysis procedure P1150 are first performed.
Embodiments of the invention relate to features of the calibration measurement procedure P1110. In particular, in some embodiments, the measurement procedure is performed using (i) ultrasound transducers (i.e. of an ultrasound array) where the transducers are each very “close” to the surface of skull outer surface 151; and/or (ii) employed to produce a plurality of ultrasound-probe events (discussed below), where for each ultrasound-probe event In some embodiments, the term ‘very close’ to the skull (or another generally round object) may be relative to one or more of (i) the size/width of the array used in the calibration measurement procedure P1110 and/or (ii) the radius of curvature (or variations thereof) of the skull and/or of a target portion thereof.
For example, each ultrasound skull-probe (USP) event (described below) may be performed by delivering ultrasound primarily (i.e. on a per-event basis) from transducer(s) whose respective centers are:
(A) located within 3 cm or within 2 cm of the outer skull-surface;
(B) localized within a very small locale (referred to as a dominant emission locale and discussed below—e.g. each locale may be defined by a sphere having a radius of at most 1 cm or at most 0.5 cm.
Not wishing to be bound by theory, providing one or both of these features may be useful for delivering ultrasound such that (i) a normal vector at a location of surface 151 where the ultrasound test signal is incident is substantially aligned with (ii) a directional energy flux vector (analogous to a Poynting vector from electricity and magnetism) of the ultrasound test signal used to generate ultrasound skull-probe (USP) event (described below). This may be useful, in some examples, for achieving a situation similar to that of
In the non-limiting example illustrated in
A. a calibration MEASUREMENT procedure P1110. For example (and without limitation), this procedure P1110 (e.g. together with at least some analysis P1120) may measure an indication of a (i) a thickness of skull 150 as a function of location on the skull and/or (ii) a speed of sound within the skull 150. The physical data acquired during the measurement procedure P1110 is subsequently processed in the calibration ANALYSIS procedure P1150 to compute relative phases used in step P1200 (e.g. for overcoming inhomogeneities in skull thickness and/or in acoustic properties of the skull and/or brain).
As will be discussed below, the calibration MEASUREMENT procedure P1110 may be described with reference to ultrasound skull-probe (USP) events, discussed below. For each USP event, a respective ultrasound test-signal (e.g. see
B. a calibration ANALYSIS procedure P1150—for example, measurement data may be used to determine an echo-time difference. The echo time-difference may be related to the time delay (tE−tD) discussed above with reference to
The analysis procedure may also comprise computation of relative phases used in step P1200 (e.g. for overcoming inhomogeneities in skull thickness and/or in acoustic properties of the skull and/or brain).
C. a sub-surface focus procedure P1200. Once the relative phases are computed in P1150, they are employed in the focus procedure P1200 where a plurality of ultrasound transducers (i.e. either the same transducers used in P1110 or different transducers) emit ultrasound according to the computed relative phases. In embodiments of the invention, P1200 may be performed such that:
In embodiments of the invention, during P1110 a plurality of Ultrasound probe (USP) events {USP-event1, USP-event2 . . . USP-eventL} are performed. In some embodiments, L≥5. In some embodiments, L≥7. In some embodiments, L≥8. In some embodiments, L≥10.
Each USP event USP-eventi (1≤i≤L) is characterized by:
Ultrasound Test Signals—Each Ultrasound probe (USP) event is associated with a respective ultrasound test signal—i.e. the ultrasound signal emitted by one or more of the transducers which is propagates to the skull (e.g. and is reflected therefrom). USP-event1 is associated with Ultrasound_test_signal(USP-event1), USP-event2 is associated with Ultrasound_test_signal(USP-event2), and so on. In one example, some or all ultrasound test signals are the same—e.g. Ultrasound_test_signal(USP-event1)=Ultrasound_test_signal(USP-event2)= . . . Ultrasound_test_signal(USP-eventL). However, this is not a requirement. In another example, each Ultrasound probe (USP) event is associated with a different ultrasound test signal.
In embodiments of the invention, the events {USP-event1, USP-event2 . . . USP-eventL} are performed serially.
In
In
In
In
A Discussion of
For ease of reference the following numbers in the figures are meant to refer to as follows
There is a need for more efficient ways of focusing a phased array of transducer element to create a high-quality ultrasound focus without recourse to an imaging model of the brain by MRI or Ultrasound. The key concept to keep in mind is that an array of emitters forms a focus peak intensity by coherent interference of the ultrasound wave-fronts from the different emitters at the peak location. The peak is deformed or destroyed by dephasing of wave-fronts at the intended peak location. Therefore, to improve or restore a sub-optimal focus peak, phase shift corrections need to be found such that coherent interference of the ultrasound wave-fronts from the different emitters at the peak location is improved. In particular, it is well known in the art of ultrasound arrays how to produce an ultrasound peak at a desired location in a uniform medium. Yet, in order to produce a focus peak within a human skull when ultrasound is produced by an array of emitters outside of the skull, there is a need to overcome the dephasing effect of the non-uniform and unknown geometry of skull bone layer.
The goal of the present invention is to provide a method and apparatus for obtaining phase shift corrections to improve the focus peak within a human skull when ultrasound is produced by an array of emitters outside of the skull.
In the following discussion, we present a general scheme of the procedures and apparatus, elaborate on several preferred embodiments, and present a numerical simulation of a realization example of an embodiment according to the invention.
As illustrated in
Yet, as illustrated in
The beam of acoustic energy emitted from the focus-emitters array 117 has a relatively wide aperture where it enters the body. Therefore, different parts of the acoustic energy, such as 123 and 124, may pass through different intervening tissue layer thickness between, which may shift the effective relative time delay of acoustic energy transmitted from respective transducer elements upon arrival to the focal zone. This phase shifting may decrease the constructive interference of the acoustic energy at the focal zone, or may even move the focal zone in an unpredictable manner. For example, an intervening skull bone layer thickness difference of 1.5 mm may introduce a phase shift of 180° at an ultrasonic frequency of one Megahertz (1 MHz), which would change desired constructive interference at the focal zone into destructive interference.
In some embodiments, ultrasound frequency lower than 1 MHz, such as between 100 KHz and 500 KHz, is used in order to reduce the aberration effect of non-bone tissue inhomogeneity in the brain.
As illustrated in
Thus, a principal goal, according to embodiments, is to determine a calibration delay set {CT(n)} alternatively “correction phase shifts” from which a calibrated focusing delay set {FPS(n)} can be computed such as {FPS(n)}={FT0(n)+CT(n)}. Thereby, the ultrasound wave paths pass through an intervening skull-tissue intermediate-layer bounded by first (outer) and second (inner) boundary surfaces within otherwise same uniform medium 99 in the path from the transducer surface to the resulting focus peak location 131 with resulting focal zone area 139, such that the focal zone area 139 is smaller than the focal zone 129 created when the focusing transducer array 110 is driven with an uncalibrated signal delay set {FT0(n)}. Hence the focus is improved. The resulting focus peak location 131 is very close to the planned target focus peak location 111 in the sense that the planned target focus peak location is within the resulting focal zone area 139.
The key problem is how to determine the calibration set {CT(n)}. The present disclosure provides methods and systems for determining the calibration set {CT(n)} and thereby creating an improved focus when the focused ultrasound waves path pass through an intervening skull-tissue intermediate-layer 150 bounded by first boundary surface 151 and second boundary surface 152 within otherwise approximately uniform out-of-skull medium 99, and/or similar approximately uniform medium of under-skull brain tissue 98, in the path from the transducer surface to the focus peak.
We conceptualize the method procedure as composed of two major stages: (i) scanning measurements and analysis procedure “PROC-1”, and (ii) calibrated focusing irradiation application procedure “PROC-2”.
The scanning procedure PROC-1 can be sub-divided into two prominent sub-tasks: (a) physical measurement procedure “MEASURE”, and (b) computational analysis “ANALYSIS” from which the key outcome is the determination of the calibration set {CT(n)}.
In the known art, the measurement procedure is an imaging scan followed by an analysis procedure which involves constructing a geometrical 3D image of the skull and assigning phase shift properties to the constructed 3D geometry. In contrast, in the present disclosure the measurement procedure comprises non-imaging non-focusing MEASURE measurements and the calibrated phase shift set elements are computed directly from the MEASURE measurements without going through the construction of a geometrical image of the skull.
The calibrated focusing irradiation application procedure “PROC-2” can be sub-divided into two prominent sub-tasks: (a) Input parameters set-up, and (b) irradiation application process.
The calibrated focusing irradiation application procedure “PROC-2” is similar to known art. Given the calibration set {CT(n)} from PROC-1, the system is activating the focusing emitters array {FE(n)} with the corrected set {FPS(n)} of input parameter phases. Various other irradiation application process characteristics (intensity, duration, etc. . . . ) are determined by clinical goals (e.g., nerve stimulation, tissue ablation, etc. . . . ).
The innovation is primarily contained in the preparatory measurement MEASURE process in PROC-1 method and apparatus and from it the computation method ANALYSIS from which the core outcome is the values of the calibration set {CT(n)} of input parameters that is used to define the corrected set {FPS(n)} of phases.
In order to facilitate the presentation, we define several semantic functional sets. “semantic” meaning is distinguished from “physical”, in the sense that, for example, physically the same component can be operated as a transceiver which is capable of both emitter functionality and receiver functionality.
These sets comprising:
The method and apparatus, according to embodiments, comprise using arrays (alternately and interchangeably referred to as “sets”) of ultrasound elements comprising:
As illustrated in
The measurement MEASURE procedure, according to some embodiments, as for example illustrated in
a. positioning the calibration-emitters array (CEA) 112 in proximity of the skull, such that the distance D1 of the CEA to the skull is smaller than 50% of the width W1 of the CEA, and positioning the array 116 of and calibration-receivers array {CR(p)} relative to calibration-emitters array (CEA) 112;
b. performing a calibration measurement procedure P110 of subjecting the skull to at least 10 local ultrasound skull-probe (USP) events P120 {USP-event1, USP-event2 . . . USP-eventL}(L≥10) such that:
For example, in some preferred embodiments, during each USP-eventj an individual calibration emitter CE(j) is activated to emit a test signal while all other calibration emitters CE(j′≠j) are kept dormant, thereby each USP event is dominated by a different calibration emitter CE(j). Since the calibration-emitters array (CEA) 112 is positioned in proximity of the skull, as illustrated in
Yet in some other embodiments it is possible to activate some calibration emitters CE(j) in parallel while still maintaining the conditions of the invention method and apparatus. As illustrated in
Continuing with the elaboration of the steps of performing a calibration measurement procedure:
The analysis ANALYSIS procedure comprises:
In particular, in some embodiment, as indicated by the flow chart of
It is known in the art how to focus a focusing array in a uniform medium with a signal delay set {FT0(n)} for focusing ultrasound with the array of focus-emitters {FE(n)} in the uniform medium. Thus, in the non-uniform medium comprising the skull, the correction phase-shifts set {CT(n)} are used to adjust the relative phase shifts focusing ultrasound with the array of focus-emitters {FE(n)} to a the plurality of corrected phase-shifts {FPS(n)}={FPS(1), . . . , FPS(N)} as a function of the correction phase-shifts set {CT(n)} and uncalibrated signal delay set {FT0(n)}. Hence, in some embodiments, the analysis procedure is further comprising of
Preferred embodiments of specific computation formulas are elaborated in the text below.
The echo time difference METD(j,k1,k2) is measured between an outer-skull-surface reflection-time and an inner-skull-surface reflection-time of ultrasound test signal that is emitted from the calibration emitter CE(j) of the pair [CE(j), CR(k)], thereby creating a time delay set {TD(n)}, the set {TD(n)} is computed such that TD(n) is a weighted sum or average of METD(j,k1,k2) elements where CE(j) are geometrically related to calibration-emitters to focus-emitter FE(n).
The emitter CE(j) is conceptualized as an ultrasound source of limited extent in space of less 22 Cm area, whether made of one or more transducers. In most cost effective embodiments, the emitter CE(j) is consisting of a single transducer. Thereby, for each given calibration-emitter:calibration-receiver triplet [CE(j),CR(k1),CR(k2)] of {[CE(j),CR(k1),CR(k2)]} triplets, measuring is performed without generating a geometrical imaging of the skull.
The number of calibration emitters in the set {CE(m)} and the number of focus emitters in the set {FE(n)} may be different, and hence also the number of measured echo time difference METD(j,k1,k2) elements in the set {METD(j,k1,k2)} may be different than the number of focus emitters in the set {FE(n)}. Therefore, in computing the may incorporate weighted averages of several measured echo time difference METD(j,k1,k2) elements within a given neighborhood of a local on the skull.
The measurement MEASURE procedure, according to some embodiments, comprises:
The analysis ANALYSIS procedure comprises:
The total set of emitters {E(q)}={E(1), E(2) . . . E(Q)} comprises of the union of {FE(n)} and {CE(m)}. The arrays {FE(n)} and {CE(m)} may overlap spatially. The sets {FE(n)} and {CE(m)} may intersect (i.e., share elements). In fact, from cost, performance and simplicity considerations, an embodiment is one in which the arrays {FE(n)} and {CE(m)} are physically the same elements. i.e., the array 112 and array 117 are one and the same), or one is a subset of the other (e.g., {CE(m)} is a subset of {FE(m)}, such that only a portion of the emitters are used for the calibration process. Therefore, some of the figures and embodiment examples are drawn as such. i.e, there is one array (112,117) of physical ultrasound emitters that serves both for the purpose of MEASURE procedure (in PROC-1) and for focusing (in PROC-2). Yet this should not be understood as limiting.
For a given MEASURE procedure, defining a plurality of distinct calibration-emitter:calibration-receiver triplets {[CE(j),CR(k1),CR(k2)]}. for measuring echo differences; METD(j,k1,k2) between an outer-skull-surface reflection-time received at receiver CR(k1) and an inner-skull-surface reflection-time received at receiver CR(k2) of an ultrasound test signal that is emitted from the calibration emitter CE(j).
The triplet [CE(j),CR(k1),CR(k2)] may be regarded as composed of two calibration-emitter:calibration-receiver pairs: the pair [CE(j), CR(k1)] for measuring the outer-skull-surface reflection-time received at receiver CR(k1) of an ultrasound test signal that is emitted from the calibration emitter CE(j); and the pair [CE(j), CR(k2)] for measuring the inner-skull-surface reflection-time received at receiver CR(k2) of an ultrasound test signal that is emitted from the same calibration emitter CE(j).
As known in known art, there are multiple methods of determining the calibration set {CT(n)} if there is a given geometry of the intermediate skull layer (e.g., skull bone layer). The problem is how to determine the skull properties and geometry non-invasively. For concreteness we define the intermediate skull layer geometry by knowledge of first boundary surface 151 shape function Z1(x,y) and the second layer 152 shape function Z2(x,y).
In known art, for skull bone, the skull bone layer geometry is determined from MRI or CT imaging. Thus, in known art the MEASURE sub-task measurement procedure PROC-1 is an imaging procedure done by MRI or CT scanning, from which various skull parameters are determined by supplemental external information. e.g., from multiple slices of MRI images a full 3D skull bone section shape is reconstructed. In addition, external information concerning speed of sound in the bone is supplemented to predict and determine the supposed time shifts created by the skull bone on ultrasound. i.e., the MRI or CT scanning is NOT by itself directly measuring ultrasound phase shift or time shift due to passage through the skull layer.
In known art ANALYSIS sub-task, each phase shift correction CT(n) to be applied to individual emitter En is determined by going through a geometrical reconstruction of the intermediate skull layer shape and considering the particular path of the ultrasound from the emitter En to the intended focus location.
In contrast, the present invention: (i) uses ultrasound emitters array not only for the irradiation PROC-2 procedure, but also for the measurement PROC-1 procedure, thereby eliminating completely the need for non-ultrasound MRI or CT at any step of the full procedure (PROC-1 and PROC-2) method and apparatus; (ii) what is measured is not geometrical imaging of the skull shape, but the direct ultrasonic shift effects of the skull. It is conjectured that the measurement is somehow capturing a measure that is proportional to the real focusing path without actually being the focusing path or focused at all; and (iii) the ultrasound MEASURE procedure method in the present invention is different from what is commonly understood as “ultrasound scanning” in known art of medical ultrasound imaging. In particular, known art medical ultrasound scanning is an imaging procedure using focused ultrasound, thereby, skull scanning and/or imaging is done by focusing ultrasound array onto the skull. In contrast, in the present invention, only individual calibration-emitter elements CE(m) are activated, without focusing on the skull or any brain tissue, during the MEASURE process of PROC-1. Preferably individual, or at most a small fraction of the calibration-emitters array elements {CE(m)}, e.g., less than 20%, are activated simultaneously. The focusing emitters array elements {FE(q)} are activated simultaneously only during the irradiation-focusing PROC-2 procedure, and not on skull tissue. Moreover, in some embodiments, the calibration-emitters are not producing a focused beam in general and specifically not focused on the skull surface.
In known art, scanning with ultrasound phase arrays, what is conventionally understood to be ultrasound scanning: (i) all of the array elements are operated to radiate simultaneously, (ii) a focused beam is used, and (iii) target area scanning is performed by steering the beam focus by way of modifying the relative phase shifts between simultaneously activated array elements. In order to find and trace the skull outer surface, known art phased array scan is moving the focus over a volume within which the skull is assumed to be residing somewhere. In contrast, none of the above is conducted in the present invention ultrasound MEASURE procedure.
The present invention MEASURE sub-task, in some embodiments, is characterized by that: (i) the array ultrasound elements are operated serially in time, such that individual array element (or small groups of elements the majority of which consisting of less than 10% of the number array elements) are operated on after the other, and preferably after the previous element signal reflection have been measured; (ii) the ultrasound beam is not focused; (iii) target area Measurement is performed by NOT by steering a beam focus, but instead by way of each individual array element (or small group of elements) measuring the small section of the intermediate skull layer (e.g., skull bone) closest to it.
In addition, in embodiments of the present invention ANALYSIS sub-task, each phase shift correction CT(n) to be applied to individual emitter En is determined directly from the associated individual emitter En MEASURE step, without going through a geometrical reconstruction of the intermediate skull layer shape and without considering the particular path of the ultrasound from the emitter to the intended focus.
In embodiments, CT(n)=(V2/V1−1)*TD(n)/2
i.e.,
FPS(n)=FT0(n)+(V2/V1−1)*TD(n)/2, where V1 is an assumed pre-determined average speed of sound in the interior brain tissue (preferably within 10% accuracy), and V2 is an assumed pre-determined average speed of sound in the intermediate skull layer (e.g., skull bone), preferably within 20% accuracy (and more preferably within 10% accuracy). The time TD(n) is determined from the time difference between the reflected signals 171 from the intermediate skull layer first surface 151 and the reflected signals 172 from the intermediate skull layer second surface 152.
In other embodiments, ANALYSIS sub-task is more conventionally performed, such that each phase shift correction CT(n) to be applied to individual emitter En is computed by going through a geometrical reconstruction of the intermediate skull layer shape and considering the particular path of the ultrasound from the emitter En to the intended focus location. Yet, in the present invention the geometrical reconstruction of the intermediate skull layer shape and thickness are determined in a novel way.
For example, as illustrated in
We interpret and associate each measured TD(n) as an estimation of the local time difference of ultrasound to so twice across the local skull bone thickness nearest to the emitter element En.
In embodiments, the calibration elements CT(n) forming the calibration set {CT(n)} are each a function of TD(n), V1 and V2. Preferably, as noted before,
CT(n)=(V2/V1−1)*TD(n)/2.
In addition to the phase corrections, amplitude corrections can be determined from the reflected test signals from individual emitters.
In addition, to maximize the transmitted intensity, a step of frequency-test scan is added to be performed prior to treatment application PROC-2, or prior to the MEASURE procedure. In the frequency-test scan, the goal is to find the frequency of maximum transmission through the intermediate skull layer in order to minimize loss of intensity at the focus (due to reflection from the intermediate bone layer) and in order to minimize heat deposition within the intermediate bone layer.
In a frequency-scan, for a local region the reflected signal intensity is measured while changing the emitter activation frequency range around a central work frequency, e.g., within 10% deviation from the central work frequency. For example, if the work frequency is chosen to be 500 KHz, a scan of frequency range between 450 KHz and 550 KHz is performed. Minimum local of reflected signal intensity indicates maximum local transmission at the associated frequency.
In embodiments, the frequency scan is performed for the array activation as a whole (rather than for local regions).
It should be noted that, in contrast to conventional focused-beam scanning used in ultrasound imaging, the amplitude of reflected beams received by receiver array is not used in the procedures disclosed herein except for (a) using the timing of peak amplitudes in reflected signals in order to acquire the time delays to be used for phase corrections (or phase shifts) for each emitter, and (b) as described in the preceding paragraphs, for optimizing the measurement frequency in order to maximize the transmitted intensity of the focused ultrasound radiation during the subsequent PROC-2 procedure. In conventional focused-beam scanning used in ultrasound imaging, on the other hand the amplitude of reflected beams is an important element of producing images or digital representations of the scanned tissues.
Thus, as illustrated in
In preferred embodiments, each individual emitter of the focus-emitters array produces ultrasound radiation which is unfocussed, in the sense that it is does not generate on its own a focused peak maximum on the skull bone layer.
Stated alternatively, the invention introduces a method for generating an ultrasound intensity-peak within a human subject brain around target-peak-location 131, by delivering ultrasound through the skull 150, the skull having an outer surface 151 and inner-facing second boundary surface 152, the method comprising first performing an ultrasound-based calibration procedure and subsequently forming an ultrasound peak around the target-peak location 131 using phase-shift data acquired during the calibration procedure, wherein:
In some embodiments, the measurement procedure, for selected calibration-emitter:calibration-receiver triplets [CE(j), CR(k)] of {[CE(m), CR(p)]}, of measuring, without imaging and without focusing ultrasound, a measured echo time difference METD(j,k1,k2) between an outer-skull-surface reflection-time and an inner-skull-surface reflection-time of ultrasound test signal that is emitted from the calibration emitter CE(m) of the pair [CE(j), CR(k)] comprises:
In some embodiments, TD(n) as a function of the measured echo time difference set {METD(j,k1,k2)} is such that, for the majority of n, TD(n)=METD(j,k1,k2) where CE(m) is the nearest calibration-emitter to the focus-emitter FE(n). In particular, there are embodiments where the array of calibration-emitters is the same physical elements as the focusing-emitters, i.e., {CE(m)}={FE(n)}.
In some embodiments, the ultrasound calibration-emitters are transceivers which can be operated also as receivers. Therefore, in some embodiments the array of calibration-emitters is comprising the same physical elements as the calibration-receivers. Hence these elements activation is switched between emitter and receiver functionality.
In some embodiments, TD(n) as a function of the measured echo time difference set {METD(j,k1,k2)} is such that, for the majority of n, TD(n) is a weighted average of the subset {METD(j,k1,k2)} including the calibration-emitters CE(m) nearest and next-nearest to the focus-emitter FE(n). For example, TD(n) may be equal to the average of selected METD(j,k1,k2) associated with nearest neighbors CE(m) calibration emitters.
The time delay set {TD(n)} is computed such that TD(n) is a weighted sum of METD(j,k1,k2) elements associated with the subset of calibration-emitters {CE(j),n} which are geometrically related to focus-emitter FE(n).
The simplest case is where the focus-emitter array {FE(n)} is positioned at the same location or in great proximity to the location of the calibration-emitters {CE(j)}. For example, where the focus emitters and calibration emitters are physically the same transducer element and the transducer array is not moved much or at all between the calibration and focusing procedures. In such a case, the geometrical association is essentially an identity. Yet, in more complex situations, where the focus-emitters array {FE(n)} is positioned at a different location from the calibration-emitters {CE(j)}, the geometrical association is more involved, as discussed below with respect to
In some embodiments, the value of TD(n) is a function of the measured echo-difference set {METD(j,k1,k2)} such that contributions from METD(j,k1,k2)} values are weighed higher where CE(j) is a nearer calibration-emitter to the location to the focus-emitter FE(n).
In some embodiments, the value of TD(n) is a function of the measured echo-difference set {METD(j,k1,k2)} such that contributions from METD(j,k1,k2)} is weighed highest where CE(j) is a nearest calibration-emitter to the location to the focus-emitter FE(n). For example, the value of TD(n) is set equal to METD(j,k1,k2)}, where CE(j) is a nearest calibration-emitter to the location to the focus-emitter FE(n).
In some embodiments, the value of TD(n) is a weighted sum of the elements of the measured echo-difference set {METD(j,k1,k2)} such that contributions from METD(j,k1,k2)} is weighed highest where CE(j) is a nearest calibration-emitter to the location to the focus-emitter FE(n). For example, the value of TD(n) is set equal to the weighted average of two METD(j,k1,k2)} associated with the two nearest calibration-emitters to the location to the focus-emitter FE(n), with the weight proportional to their distance from the focus-emitter FE(n).
In some embodiments, FPS(n)=(FT0(n)+CT(n)), where CT(n) is a function of TD(n)/2 and where {FT0(n)} is a delay set used for a uniform-space focusing of said ultrasound array.
In some embodiments, the complete set of ultrasound emitters (CSUE) is the sum of ultrasound calibration-emitters {CE(m)} and the set of ultrasound focus-emitters {FE(n)}.
In some embodiments, a dormant emitter is considered to be an emitter whose simultaneous intensity is less than 80% of the intensity of the highest intensity active emitter.
In some embodiments, a dormant emitter is considered to be an emitter whose simultaneous intensity is less than 50% of the intensity of the highest intensity active emitter.
In some embodiments, a dormant emitter is considered to be an emitter whose simultaneous intensity is less than 20% of the intensity of the highest intensity active emitter.
In some embodiments, a dormant emitter is considered to be an emitter whose simultaneous intensity is less than 10% of the intensity of the highest intensity active emitter.
In some embodiments, the position of the array 117 of focus-emitters {FE(n)} and array 112 of calibration-emitters {CE(n)} during the focusing is the same as during the measurement procedure.
In some embodiments, the position of the array 112 of calibration-emitters {CE(n)} remains the same throughout the measurement procedure.
As illustrated in
In some embodiments, the system is such wherein the function CT(n)=(V2/V1−1)*TD(n)/2, where V2 is an estimated average speed of sound in the intermediate-layer within at least 20% accuracy and V1 is an average speed of sound in the target tissue within at least 10% accuracy.
In some embodiments, the system is such wherein the control module for selectively activating sub-sets of calibration-emitters {CE(m)} is selectively activating individual array elements CE(m)) while maintaining other array elements in a dormant state.
In some embodiments, the system is such wherein the detection system serially performing a measurement of such echo signals from a plurality of selected array calibration-emitters elements, is serially performing the measurement on a majority of array calibration-emitter elements,
In some embodiments, the system is further comprising a positioning module capable of maintaining the position of the emitter phased array with respect to the target tissue.
In some embodiments, the system is such wherein a transducer emitter element E(n) can physically serve also as a receiver element R(m), associated with same emitter n (m=n) and/or associated with a different emitter element n′ (m=n′).
In some embodiments, the system is having at least two modes of activation, (i) a focusing beam mode at which the majority of the emitter array elements are simultaneously emitting ultrasound to create a focus peak, and (ii) an calibration mode at which only a minority of emitter elements are simultaneously activated and at least one receiver element is sensing the reflected signals and transmitting them for analysis to the computation module. In some embodiments, the minority of emitter array elements is less than 10% of the array elements, or less than 1% of the array elements.
In some embodiments, the minority of emitter array elements is an individual single element.
It helps to highlight from the outset certain distinguishing feature of the present invention embodiments in comparison with known art.
It is well known in the art that a proper coupling medium 118 needs to fill the gap between the transducer surface 113 and skin surface, in order to reduce ultrasound reflection from the entry interface into the subject body. In embodiments, the coupling medium is deformable, such as a liquid or gel, in order to be able to conform to the skin contour. In embodiments, the coupling medium has speed of sound similar to the skin tissue.
Figs.
FIG. 2
Since any of the above and also combination of the above transducer array models can be used to realize the present invention, we use a generic representation, illustrated in
The calibration emitters array (CEA) has a characteristic size CS(CEA). To achieve good focus at depth of about 26 Cm, the calibration array needs to be of about the same size. Hence the CEA characteristic size is at least 5 cm, preferably between 26 Cm to 10 cm. For completeness, since stray elements are of little effect, we shall define the characteristic size CS(CEA) as a diameter of the smallest sphere which includes a 80% of the emitters calibrations emitters array.
Each member CE(j) of the calibration emitters array is displaced from the skull by DISPL[CE(j)]. For good calibration, the calibration emitters (hence, the calibration array) should be in proximity of the skull.
Preferably the majority of calibration emitters are within 23 Cm of the skull skin surface, better within 22 Cm of the skull skin surface, better within 1 cm of the skull skin surface. i.e., for a majority of the calibration emitters of the calibration emitters array, a ratio between DISPL(CE(j)) and CE(CEA) is at most ½, preferably at most ⅓, preferably at most ¼.
As illustrated in
As illustrated in
The key problem is how to determine the calibration set {CT(n)}. The present invention provides a method and system for determine the calibration set {CT(n)} and thereby creating an improved focus when the focused ultrasound waves path pass through an intervening skull-tissue intermediate-layer bounded by first and second boundary surfaces within otherwise approximately uniform medium 99 in the path from the transducer surface to the focus peak.
We conceptualize the method as a procedure composed of two major stages: (i) Measurement procedure “PROC-1”, and (ii) calibrated focusing irradiation application procedure “PROC-2”.
The measurement procedure PROC-1 can be sub-divided into two prominent sub-tasks: (a) physical measurement procedure “MEASURE”, and (b) computational analysis “ANALYSIS” from which the key outcome is the determination of a time-set {TD(n)} from which is determined the calibration set {CT(n)}, where each element CT(n) is a function of corresponding TD(n).
The calibrated focusing irradiation application procedure “PROC-2” can be sub-divided into two prominent sub-tasks: (a) Input parameters set-up, and (b) irradiation application process.
In the calibrated focusing irradiation application procedure “PROC-2”, the calibration set {CT(n)} from PROC-1, the system is activating the emitters array with the corrected set {FPS(n)} of input parameters phases. The irradiation application process itself is determined by clinical goals (e.g., nerve stimulation, tissue ablation, etc. . . . ).
The innovation is primarily contained in the preparatory measurement process PROC-1 method and apparatus and from it the core outcome is the values of the calibration set {CT(n)} of input parameters that is used to define the corrected set {FPS(n)} of phases.
As illustrated in
In the path to the focus peak 111 from the emitter 141m as calibration emitter CE(m), the sound wave passes once through the intermediate skull layer 150 width. In contrast, the reflected test signals, 172 and 174, from boundary 152, each passes twice through the intermediate skull layer 150 width (corresponding to the incident and reflected portions of the path within the layer 150). Therefore, half of the reflected test signals time difference TD(n), i.e., (TD(n)/2), is a good approximation to the time shift contribution of the intermediate skull layer 150 to the total travel time of the focused beam from the particular emitter 141m as calibration emitter CE(m) to the focus peak location 111. Therefore, in embodiments of the present invention, adjusting the relative time shift by subtracting (TD(n)/2) from FT0(n) substantially eliminates the phase shift contribution of the intermediate skull layer 150 to the total travel path from the emitter 141m as calibration emitter CE(m) to the focus peak 111. Hence, in embodiments the calibrated focusing delay set {FPS(n)}={(FT0(n)+CT(n))} is then use for irradiating the target tissue using the ultrasound transducer array emitter elements {E(n)} having a corrected-delay {F(n)}={FPS(n)} to create a focus peak 131 with an calibrated focusing focal zone 139 of smaller cross section area than if irradiated with uncalibrated focusing delay set {FT0(n)}.
It remains to be decided which receiver sensor measurement to use for the determination of TD(n). Ideally, one would like to have a reflected path which is exactly twice the length of the focused beam path within the intermediate skull layer 150.
One preferred approximation, as illustrated for the path 172, is a test signal path for which at least the incident portion of the test signal path to be the same as the focused beam path 170 to the focus peak 111. For such an embodiment, fast switching needs to be operated to switch the emitter transducer to receiver mode of operation within the duration of the reflection time. To a good approximation, for ray 370, defined to be the ultrasound propagation “ray” from emitter 341m as would be in the selected uniform medium 99 traversing the focus peak 111, the path of ray 370 through the intermediate skull layer 150 is very close to the real path of ultrasound ray from emitter 141m as calibration emitter CE(m) to the focus peak 111.
Another preferred embodiment approximation is one for which the same physical transducer array element is used for both test signal emitter and as receiver sensor. Another preferred embodiment approximation is one for which, as illustrated for the path 171, for a given emitter 141m as calibration emitter CE(m) the associated receiver sensor 161p as paired calibration receiver CR(p) is a neighboring (e.g., nearest neighbor) transducer sensor element 161p as paired calibration receiver CR(p). A question in such embodiments is which direction of neighbor to choose as sensor placement relative to the emitter element. For example, should it be one to the left or to the right of the emitter element. The better approximation depends on the relative curvature of the entry boundary 151 of the intermediate skull layer 150 compared with the curvature of the focus beam at that boundary surface. For example, as illustrated in
There is no need for unique exclusive one-to-one association of a sensor to an emitter. In some embodiments, the number of receiver sensors is smaller than the number of emitter elements in {E(n)}. For example, as illustrated in
In some embodiments, the emitter and receiver sets are physically distinct. In some other embodiments, the emitter and receiver sets are overlapping. For example, referring to
It is preferred to be able to selectively activate individual emitter and/or receiver elements from the physical arrays sets. For example, for test signal and calibration procedures, one is preferably activating the array emitter elements serially or in sub-sections (i.e., not all together as for focus creation), and respectively receiving and/or analyzing the echo received signals only of the associated receiver sensors. That is relevant both for embodiments where the emitter and receiver sets are physically distinct and for embodiments where they are overlapping. In embodiments, the switching control and activation of the emitter 112 and/or receiver 116 arrays are managed by a controller module 255.
A focusing transducer 110 comprising an emitter phased array {E(n)} 112 of ultrasound transducer elements for generating an ultrasound focus in the target tissue. At least most transducer elements having means connected thereto for variably setting a delay for that emitter transducer E(n), including a delay module for setting delay set {F(n)} of signals delay to associated emitter elements {E(n)}.
A receiver array 116 comprising a plurality of ultrasound receiver elements {R(n)} associated with emitter elements {E(n)}; The receiver array 116 is connected to and controlled by a receiver control module 170.
A transmit/receive controller “T/R controller” module 255 for directing the activation and switching of individual elements of the emitter and receiver arrays. The T/R module 255 comprising connections to the emitters array 112, to the receiver array 116, to the signal generator module 140 and to the control & computation module 160.
An emitter mux module 145 receives input signal from the signal generator 140, the associated activation delay set {TD(n)} from control module 160, and the array elements activation directed by T/R module 255. The emitter mux 145 transmit the activation signals to the emitter array 112 elements {E(n)}.
For focusing-mode all, or at least the majority, of the emitter array 112 elements {E(n)} are activated simultaneously, with associated delay set {F(n)} of a common action-signal, to generate a focus peak at a certain location.
For calibration-mode a minority, preferably one, of the emitter elements E(n) is activated at a time, preferably with a distinct “test-signal”, which may be preferably different from the action-signal. The test signal reflected echo-signal is received at an associated “paired” receiver element of the receiver set 116 and transmitted for analysis to receiver control module 170. It is expected that the received reflected signal would include multiple reflections from various material boundaries such as skin surface, fatty tissue, and the intermediate skull bone layer 150. As highlighted by the table of
After calibration mode procedure, the focusing mode is activated with the emitter array 112 driven with an calibrated focusing time-delay set {FPS(n)}={FT0(n)+CT(n)} to create an calibrated focusing focus peak.
For effective focus improvement by the invention calibration procedure. The positioning module 165 maintains the transducer module 110 at a fixed orientation relative to the target tissue for both the calibration procedure and the focusing mode activation.
As previously discussed, better focusing calibration would be obtained if the test signals path through the intermediate tissue layer is better matching (i.e., closer) to the path through the intermediate skull layer of focusing beam to the focus peak location. As illustrated in
In embodiments, the preferred associated calibration-receiver for optimal path matching is not fixed for all intermediate skull layers, but is dependent on the orientation of the intermediate tissue boundary layers relative to the transducer focusing-axis 310. For example, as illustrated in
As illustrated in
At different receiver sensor the reflection times may be different. Therefore, the selection of receiver sensors at which time of reflected beams is measured is affecting the time delay correction calibration. i.e., for each calibration-emitter CE(m) there is an associated “paired” calibration-receiver CR(p) at which the echo signal is measured, thus defining the pair [CE(j),CR(k1),CR(k2)]. Different pairing may lead to different measurement and hence different calibration outcome. Some of the difference between various embodiments relates to the method or protocol of selecting the pairing of the calibration-emitter CE(m) there is an associated “paired” calibration-receiver CR(p) at which the echo signal is measured for the calibration ANALYSIS.
The time difference between these reflection signals T2(n)−T1(n) is determined by the time-shift determiner module 172 to create the reflection time difference TD(n)=T2(n)−T1(n). Repeating the process serially for multiple, preferably most or preferably all, of the emitter elements, lead to obtaining a reflection time difference set {TD(n)}.
After calibration mode procedure, the focusing mode is activated with the emitter array 112 driven with a calibrated focusing time-delay set {FPS(n)}={FT0(n)+CT(n)} to create an calibrated focusing focus peak.
The previously discussed embodiment, with reference to
In embodiments, the sensor for timing T1 of the arrival of reflection 571 is timed at receiver 561p at which reflected beam 571 is detected at peak intensity. In embodiments, the sensor for timing T2 the arrival of reflection 572 is timed at receiver 562p at which reflected beam 572 is detected at peak intensity. Repeating the test process for each emitter E(n) we obtain the reflected signals arrival times T1(n) and T2(n).
The time difference between these reflection signals T2(n)−T1(n) is determined by the time-shift determiner module 172 to create the reflection time difference TD(n)=T2(n)−T1(n). Repeating the process serially for multiple, preferably most or preferably all, of the emitter elements, lead to obtaining a reflection time difference set {TD(n)}.
After calibration mode procedure, the focusing mode is activated with the emitter array 112 driven with an calibrated focusing time-delay set {FPS(n)}={FT0(n)+CT(n)} to create an calibrated focusing focus peak.
In addition to phase corrections (i.e., time delays) for individual emitters En, also amplitude corrections are determined in embodiments. For example, for a given emitter test signal, the fractional intensity which is collected at the receiver peak (and optionally including the intensity of the nearest neighbors receivers also) for both the first boundary reflection and second boundary layer reflected echo signal is indicative of the remaining transmitted intensity that reaches the focus peak. Higher reflected fraction RI(n) means lower transmitted intensity fraction “TI(n)” contributing to the focus peak. Hence, in embodiments, the amplitude correction A(n) for each emitter En of the set {En} is a function of TI(n), e.g., proportional to 1/TI(n). For example, in embodiments, in order to obtain more uniform contribution to the focus peak from each emitter En, the intensity emitted from emitter En is set to be [1/TI(n)]*A0, where A0 is the intended average intensity of the transducer emitters array.
In addition, to maximize the transmitted intensity, a step of frequency-test scan is performed prior to treatment application. Typical skull bone thickness ranges between 4 mm to 12 mm. Typical ultrasound treatment frequency ranges between 0.25 MHz to 2 MHz, which at average speed of sound in bone of 3000 m/s translates to a range of wavelength between 1.5 mm to 12 mm. Hence, % wavelengths range in size between 0.375 mm to 3 mm. Due to the % wavelength maximum transmission effect, for a given skull bone sample, for any given desired central frequency W0 for treatment (e.g., 1 MHz, with associated wavelength of 3 mm for bone speed of sound of 3000 m/s) there will be a particular optimal frequency W1 in the neighborhood of W0 for which the transmitted ultrasound total beam intensity across the skull bone layer is maximized. This can respectively be detected as a minimum in the reflected intensity.
Since the intermediate skull layer (e.g., skull bone) thickness is non-uniform, it is possible and even likely that, the optimal frequency W1 of maximum transmitted intensity is different for different sub-regions of the surface area 155 under the full transducer area. Therefore, in embodiments, as illustrated in
The optimal frequency W1 is determined independently for each sub-section. Then, the treatment is delivered with each sub-section driven at its own local optimal frequency, in parallel or serially with other sub-sections of the transducer array.
In all of the above examples, the key in the preparatory MEASURE process is to get separate measurements of individual emitters En (or small local group of emitters). In practice, the signals from sufficiently distant emitters do not mix or interfere at the small distance of the MEASURE reflection measurements. Therefore, in embodiments, the MEASURE process can be performed simultaneously on multiple emitters sub-set of the full array set {En} such that the distance between the array elements is larger than the distance between each array element and the intermediate skull layer 151.
In embodiments, of the MEASURE procedure the array 112 is positioned close to the skull, hence close to the skull boundary surface 151. In embodiments, the curvature of the array is very close to concentric with the generic curvature of the human skull layer, typically between 26 Cm and 122 Cm radius of curvature. Yet, different areas of the skull have some difference in curvature, and also there is a range of curvature differences between human individuals. Therefore, in some preferred embodiments, instead of fixed array curvature, the array curvature used in a particular instance can be modifiable to better match to the area of the skull that is irradiated in practice. i.e., prior to the start of PROC-1 as discussed above, there is an initial “curvature adjustment” step.
The MEASURE procedure in the simulation consisted of serially performing the same MEASURE step on all the 20 emitters forming the emitters array 112 set {En}.
Finally,
Importantly, we note that throughout the process (i.e., MEASURE, ANALYSIS, and irradiation), as demonstrated in the simulation, there was no need to recourse to a geometrical reconstruction of the macro shape of the intermediate skull layer 150. Instead, local measurements where directly translated to local time shifts and calibration factors.
In the above embodiments examples the PROC-1 MEASURE process and the PROC-2 irradiation process where both performed with the emitters array 112 at the same relative position with respect to the intermediate skull layer 150. But this is not necessarily the case.
Elements of the receivers array 116 are marked by dotted hexagons, and the remaining elements are emitters. In this embodiment, the same physical emitters serve for both the calibration emitters array 112 and the focusing emitters array 117. Activated emitter elements are marked by line-filled hexagon (as element 241a in
The method and system of the present invention can be employed either
In the following discussion, we elaborate the case when the PROC-1 MEASURE process and the PROC-2 irradiation process are performed with the array 112 of ultrasound calibration-emitters {CE(m)}={CE(1), CE(2) . . . CE(M)}, and an array 117 of ultrasound focus-emitters {FE(n)}={FE(1), FE(2) . . . FE(N)} are at different relative positions with respect to the intermediate skull layer 150.
The key outcome, of this modified embodiment of the calibration ANALYSIS procedure, is the construction of a fictitious-skull geometry and associated phase shifts distribution. The PROC-2 irradiation focusing process is then performed as if the irradiated intermediate skull layer object is of the properties of the fictitious-skull rather than the real patient skull.
For removal of any doubt, the array 112, of the emitters used for calibration-emitters function, can be both geometrically and physically different from the array 117 of the emitters used for the focus-emitters function. From the calibration procedure outcome, it is the data of the fictitious-skull geometry and phase shifts which is the input in to the focusing procedure. Hence, the focusing array configuration can in principle be independent of the configuration of the calibration array used for producing the phase shifts input data.
For the sake of simplicity of the drawings and discussion, we shall make the presentation in terms of an engineering embodiment where there is one physical system of emitters that serves for both the calibration-emitters {CE(m)} and the focus-emitters {FE(n)}. This is in reality the simplest to construct from an engineering realization point of view. Yet, functionally we continue to designate the array as 112 when the emitters are used as calibration-emitters function, and designate the array as 117 when the emitters are used as focus-emitters function.
In order to distinguish this more particular embodiment of the PROC-1 MEASURE process, we shall refer to is as a MEASURE-Trace process. The PROC-1 MEASURE-Trace process implements serial activation of an array of ultrasound probes, each placed at an associated probing-location relative to the skull. The PROC-1 MEASURE-Trace process implements the serial activation of an array of ultrasound probes to trace, at discrete locations, the shape of the external skull and measure an effective local ultrasound time delay phase shift directly. An effective fictitious skull model shape and thickness model is produced as an outcome of the PROC-1 MEASURE-Trace process. Subsequently, for the PROC-2 irradiation focus process, the focus-emitters array 117, which can be placed at any geometrical position and orientation relative to the traced external skull geometry, is activated for focusing at any location within the skull, while computing the focusing phase shifts with respect to the effective fictitious skull model, taking the effective fictitious skull model external surface as coinciding with the human subject skull, using known in the art methods for focusing through a known layer geometry. An assumed speed of sound properties of the skull layer is also taken into consideration as is known in the art.
Yet, as illustrated in
In embodiments, in order to enable preferred embodiment where the PROC-1 MEASURE process and the PROC-2 irradiation process where performed with the emitters array 112 at the different relative position with respect to the intermediate skull layer 150, we introduce in PROC-1 an optional supplemental ANALYSIS step—labeled VIRTUAL LAYER CONSTRUCTION—and introduce in PROC-2 a supplemental irradiation preparation optional step—labeled PATH ANALYSIS. The key additional step is the VIRTUAL LAYER CONSTRUCTION sub-process at the end of the ANALYSIS process.
The VIRTUAL LAYER CONSTRUCTION outcome is a geometrical construction which overlap a portion of the intermediate skull layer which is facing the emitters array 112 during the MEASURE process.
The virtual layer construction comprises of at least two properties: (A) a geometrical surface shape z=B1(x,y), and (B) an associated effective geometrical-surface-time-difference value function GSTD(x,y,z) as a function of the measured echo time difference set {METD(j,k1,k2)}.
The time delays set {TD(n)} is generated as a function of the geometrical-surface-time-difference value function GSTD(x,y,z). In preferred embodiments, as illustrated in
Generating a corrected focusing phase-shifts set {FPS(n)} from the measured time delays set {TD(n)} is following the same procedure as previously discussed.
In some embodiments the coordinate system (x,y,z) axis are fixed relative to the geometry of the emitters array {CE(j)}. Thereby, positioning of the focusing emitters array {FE(n)} can be determined by knowledge of the relative positioning of the focusing emitters array to the emitters array {CE(j)}.
In some embodiments the coordinate system (x,y,z) axis are fixed relative to the geometry of the skull first boundary layer 151. Thereby, positioning of the focusing emitters array {FE(n)} can be determined by knowledge of the relative positioning of the focusing emitters array to the skull first boundary layer 151.
In some embodiments, the value of GSTD(x,y,z) is a function of the measured echo-difference set {METD(j,k1,k2)} such that contributions from METD(j,k1,k2)} values are weighed higher where CE(j) is a nearer calibration-emitter to the location to the geometrical point (x,y,z).
In some embodiments, the value of GSTD(x,y,z) is a function of the measured echo-difference set {METD(j,k1,k2)} such that contributions from METD(j,k1,k2)} is weighed highest where CE(j) is a nearest calibration-emitter to the location to the geometrical point (x,y,z). For example, the value of GSTD(x,y,z) is set equal to METD(j,k1,k2)} is weighed highest where CE(j) is a nearest calibration-emitter to the location to the geometrical point (x,y,z).
In some embodiments, the value of GSTD(x,y,z) is a weighted sum of the elements of the measured echo-difference set {METD(j,k1,k2)} such that contributions from METD(j,k1,k2)} is weighed highest where CE(j) is a nearest calibration-emitter to the location to the geometrical point (x,y,z). For example, the value of GSTD(x,y,z) is set equal to the weighted average of two METD(j,k1,k2)} associated with the two nearest calibration-emitters to the location to the geometrical point (x,y,z), with the weight proportional to their distance from the geometrical point (x,y,z).
For better outcome, it is advantageous if the geometrical-surface-time-difference value function GSTD(x,y,z) geometry, specified by a shape function z=B1(x,y), is an approximation of the outer boundary surface 151 of the skull. For example in some embodiments z=B1(x,y)) is a curved surface with a curvature of a typical human skull facing the location of the calibration emitters array.
In some embodiments, the geometrical shape function z=B1(x,y) construction is derived from the set of parameters obtained in the invention SCAN process and ANALYSIS as deduced from the echos reflection from the skull surface. For example, the timing of first echo from a calibration emitter CE(j) serves to determine the geometrical distance of the outer skull surface from the geometrical location of the emitter CE(j).
Generation of a measured time delays set {TD(n)} as a function of the geometrical-surface-time-difference value function GSTD(x,y,z) is performed by the following procedure: Since {TD(n)} is associated with the focus-emitters array {FE(n)}, the pre-determined intended positioning of the focus-emitters array {FE(n)} with respect to the geometrical shape function z=B1(x,y) of the geometrical-surface-time-difference value function GSTD(x,y,z) needs to be assumed as given (even if the physical array of focusing emitters is not yet placed there).
As indicated in
The fictitious skull layer thickness at various locations, i.e., the local distance between Z1(x,y) surface and Z2(x,y) surface can be set in principle from the measured echo time difference S125 METD(j,k1,k2)=T2−T1, between detecting S123 at receiver CR(k1) and measure time T1 of first echo from the skull first boundary surface and detecting S124 at receiver CR(k2) and measure time T2 of second echo from the skull second boundary surface, of an unfocused ultrasound test signal that is emitted S122 from the calibration emitter CE(j).
All publications, patent applications, patents, and other references mentioned herein are incorporated by reference in their entirety. In case of conflict, the patent specification, including definitions, will prevail. In addition, the materials, methods, and examples are illustrative only and not intended to be limiting.
The present invention has been described using detailed descriptions of embodiments thereof that are provided by way of example and are not intended to limit the scope of the invention. The described embodiments comprise different features, not all of which are required in all embodiments of the invention. Some embodiments of the present invention utilize only some of the features or possible combinations of the features. Variations of embodiments of the present invention that are described and embodiments of the present invention comprising different combinations of features noted in the described embodiments will occur to persons of the art.
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/IB2019/000185 | 1/31/2019 | WO | 00 |
| Number | Date | Country | |
|---|---|---|---|
| 62462398 | Feb 2017 | US | |
| 62468473 | Mar 2017 | US | |
| 62474715 | Mar 2017 | US | |
| 62492193 | Apr 2017 | US |
| Number | Date | Country | |
|---|---|---|---|
| Parent | 15902361 | Feb 2018 | US |
| Child | 17425347 | US |
