APPARATUS AND METHOD FOR ULTRASOUND BEAM SHAPING
20210187330 · 2021-06-24
Assignee
Inventors
- Michael R. Bailey (Seattle, WA, US)
- Adam D. Maxwell (Seattle, WA, US)
- Akshay Purushottamji Randad (Seattle, WA, US)
- Mohamed Abdalla Ghanem (Seattle, WA, US)
Cpc classification
A61B17/225
HUMAN NECESSITIES
B33Y80/00
PERFORMING OPERATIONS; TRANSPORTING
G16H20/40
PHYSICS
International classification
Abstract
Apparatus and method for ultrasound beam shaping are disclosed herein. In one embodiment, an ultrasonic therapy system is configured to apply ultrasound to a target in a body. The system includes: an ultrasonic transducer configured to generate the ultrasound; and a customizable lens configured to focus the ultrasound onto a focal area of the target. The target is an object or a portion of the object in the body. The customizable lens is designed and produced based on at least one acquired image of the target.
Claims
1. An ultrasonic therapy system configured to apply ultrasound to a target in a body, comprising: an ultrasonic transducer configured to generate the ultrasound; and a customizable lens configured to direct the ultrasound onto an area of a target, wherein the target is an object or a portion of the object in the body, and wherein the customizable lens is designed and produced based on a determined size, shape or location of at least one target being a design input for the customizable lens, wherein the customizable lens is designed by an iterative numerical computational method based on the at least one target such that a focal area of the customizable lens corresponds-to or exceeds a size and a shape of the object.
2. The system of claim 1, wherein the target is defined by at least one acquired image of the object in the body.
3. The system of claim 1, wherein the customizable lens is a three-dimensional (3D) printed lens.
4. The system of claim 1, wherein the iterative numerical computational method is based on an iterative angular spectrum approach (IASA).
5-6. (canceled)
7. The system of claim 1, wherein ultrasound pressure phases are focused onto the focal area.
8. The system of claim 3, wherein the customizable lens is configured to produce a plurality of target distributions of the ultrasound at a corresponding plurality of focal distances from the customizable lens.
9. The system of claim 3, wherein the customizable lens is configured to produce the plurality of target distributions of the ultrasound at a corresponding plurality of ultrasound frequencies.
10. The system of claim 2, wherein the at least one acquired image of the object is modified to introduce an asymmetry in a target ultrasound field of the focal area of the object.
11. The system of claim 10, wherein the customizable lens produces multiple high-pressure areas within the target ultrasound field.
12. The system of claim 1, further comprising: a mechanism configured to mate the customizable lens with the ultrasonic transducer; and an interface material configured to temporarily attach the customizable lens with the ultrasonic transducer.
13. The system of claim 12, wherein the mechanism is selected from a group consisting of a quick-change clamp, a hinge and a bolt.
14. The system of claim 1, wherein the ultrasound transducer is a phased array transducer comprising a plurality of ultrasound sources.
15. The system of claim 1, wherein the plurality of ultrasound sources of the phased array transducer is arranged along a curved surface.
16. A method for applying an ultrasound to a target in a body, comprising: defining a customizable lens based on a determined size, shape or location of at least one target, wherein the customizable lens is designed by an iterative numerical computational method, and wherein the target is an object or a portion of the object in the body; generating the ultrasound by an ultrasonic transducer; and directing the ultrasound onto an area of the object by the customizable lens such that a focal area of the customizable lens corresponds to or exceeds a size and a shape of the object.
17. The method of claim 16, further comprising: acquiring an image of the object, wherein the iterative numerical computational method is based on an iterative angular spectrum approach (IASA).
18. The method of claim 17, further comprising: acquiring additional images of the object while the target is being treated; and based on acquiring the additional images of the object, modifying the customizable lens.
19. The method of claim 16, further comprising manufacturing the customizable lens by three-dimensional (3D) additive-printing.
20. The method of claim 16, further comprising: applying an interface material to a surface of the customizable lens; mating the customizable lens with the ultrasonic transducer via the interface material; and after directing the ultrasound onto the area of the body, removing the customizable lens from the ultrasonic transducer.
21. The method of claim 16, wherein the ultrasound transducer is a phased array transducer comprising a plurality of ultrasound sources.
22. The method of claim 16, wherein directing the ultrasound onto the area of the body includes focusing ultrasound pressure amplitude distribution or ultrasound pressure phase distribution over the focal area.
23. The method of claim 16, further comprising: prior to defining the customizable lens, introducing an asymmetry in a target ultrasound field of the focal area; and in response to introducing the asymmetry, generating multiple high-pressure areas within the target ultrasound field by directing the ultrasound onto the focal area by the customizable lens.
24. The method of claim 16, further comprising generating a plurality of target distributions of the ultrasound at a corresponding plurality of focal distances from the customizable lens.
25. A non-transitory computer readable medium having computer executable instructions stored thereon that, in response to execution by one or more processors of one or more computing devices, cause the one or more computing device to perform actions comprising: acquiring an image of a size, shape or location of an object in a body; and determining a shape of a customizable lens based on the acquired image of the object in the body, wherein the customizable lens is designed by an iterative numerical computational method, wherein the customizable lens is configured for mating with an ultrasound transducer, and wherein the customizable lens is configured to direct the ultrasound transducer onto an area at the object such that a focal area of the customizable lens corresponds-to or exceeds a size and a shape of the object.
Description
DESCRIPTION OF THE DRAWINGS
[0028] The foregoing aspects and many of the attendant advantages of the inventive technology will become more readily appreciated as the same are understood with reference to the following detailed description, when taken in conjunction with the accompanying drawings, wherein:
[0029]
[0030]
[0031]
[0032]
[0033]
[0034]
[0035]
[0036]
[0037]
[0038]
[0039]
[0040]
[0041]
[0042]
[0043]
[0044]
[0045]
[0046]
[0047]
[0048]
[0049]
[0050]
[0051]
[0052]
[0053]
[0054]
[0055]
[0056]
[0057]
DETAILED DESCRIPTION
[0058] While several embodiments have been illustrated and described, it will be appreciated that various changes can be made therein without departing from the spirit and scope of the claimed subject matter.
[0059]
[0060]
[0061] In some embodiments, the customizable lens may be designed using the iterative angular spectrum approach (IASA). The method is described with reference to designing a customizable lens, but the method can also determine distribution and operation of the elements of an ultrasonic phased array transducer.
[0062] In some embodiments, an algorithm implements IASA numerically by iteratively comparing simulated conditions at the target focal surface against the target conditions at the focal surface. In some embodiments, an algorithm implements IASA numerically by iteratively comparing simulated conditions at the focal surface against the target conditions at the focal surface; and the complex pressure distribution at the source to the results from the previous iterative step.
[0063] In a first step, the algorithm introduces lens geometry, propagating wave front, and target focal surface in a given medium. The target focal surface may be defined by its pressure pattern (p), made up of an amplitude map (A) and a phase map (Φ). The target focal surface is located some known distance from the lens.
[0064] The pressure wave equation includes amplitude and phase functions describing pressure at a given position in Euclidean space:
p(x,y,z)={circumflex over (p)}(x,y,z)e.sup.jΔΦ(x,y,z) (1)
where {circumflex over (p)}(x, y, z) and ΔΦ(x, y, z) are the amplitude and phase functions, respectively.
[0065] The IASA method uses fast Fourier transform (FFT) and inverse fast Fourier transform (IFFT) methods to converge to an optimum error criterion, calculated as an error between the target focal surface and conditions at the focal surface. The general form of the FFT equation in Euclidean space is shown in Equation 2:
P(k.sub.x,k.sub.y)=∫∫.sub.−∞.sup.+∞p(x,y,0)e.sup.−j(k.sup.
The output of the FFT equation, P(k.sub.x,k.sub.y), gives an angular spectrum, where k.sub.i is the wavenumber in i space. The IFFT equation, excluding the evanescent wave components, is shown in Equation 3:
which provides the conditions at the focal surface in Euclidean space from the angular spectrum (P).
[0066] In the initial iteration of the loop shown in
P(k.sub.x,k.sub.u,z)=P(k.sub.x,k.sub.y,0)e.sup.jz√{square root over (k.sup.2−k.sub.x.sup.2−k.sub.y.sup.2)} (4)
which is used to calculate both propagation and backpropagation through the given medium between the focal surface and the lens. The propagating wave front then propagates through the lens and the given medium to produce an angular spectrum for a propagated wave front at the focal surface (the conditions at the focal surface).
[0067] As shown in
[0068] To account for the near field effects, the IASA incorporates a back-propagation of the propagated wave front from the focal surface to the lens, shown as a clockwise lower arrow in
[0069] In addition to conventional IASA method, the method uses the multiple checks in the convergence criterion to meet our desired goals. The algorithm iteratively compares the convergence of the simulated conditions to the target image specified at each target location. Second, after the first iteration step and in parallel to the previous check, the algorithm compares the complex pressure distribution at the source to that of the previous step as well to speed up and improve the convergence calculation criterion. The comparisons in the previous two checks are specified to be within a specific error tolerance below which convergence to the optimal hologram solution is achieved. Finally, a maximum number of iterations is determined for each run, such that when it is exceeded the method terminates and saves the optimal hologram solution. The error tolerance and maximum number of iterations is determined based on the complexity of the hologram, such as, the number of target locations for phase and or amplitude at different frequencies. These checks of convergence are checked at each iteration step to yield the optimal solution.
[0070] Incorporating back propagation into an iterated forward propagating wave equation permits a more precise calculation of the conditions at the focal surface for subsequent adjustment of the lens geometry. With each cycle of forward propagation and back propagation the conditions at the focal surface and the conditions at the lens converge to an optimal solution.
[0071] An output of the IASA algorithm is the lens geometry. As described in Equation 5, a spatial thickness parameter describes the lens geometry by taking into account the transmission coefficient (α) of the system, including acoustic impedances (Z) of the lens material (h), the given medium (m), and a transducer (t), a source of acoustic waves:
The thickness of the lens (T) can be calculated from the angular spectrum of the converged solution by creating a phase map for the surface of the lens. The lens creates constructive and destructive interference in the near-field by introducing phase offsets (ΔΦ) in the propagating wave front as it passes through the holographic lens. The thickness of the lens is described as follows in Equations 6-7:
Δω(x,y)=(k.sub.m−k.sub.h)ΔT(x,y) (6)
where T(x,y)=T.sub.0−ΔT(x,y). (7)
[0072] The IASA algorithm is capable of designing a lens that produces multiple target focal surfaces at as many distances from the lens in a given medium. To accomplish this, the IASA algorithm separately incorporates the backpropagation from the wave equations of each of the target focal surfaces when modulating the propagating wave equation.
[0073] In a similar manner, the Euclidean coordinate space of the solution permits a phased array element distribution to produce one or more target focal surfaces, by calculating IASA converged solutions for multiple propagating waves from an array of transducers. The IASA method can be used with different transducer geometries. For instance, for a focused transducer, the exact pressure field can be simulated and verified through holographic scanning in a plane. Next, the pressure field at the transducer aperture (obtained by back-projection) is used as the initial boundary condition over which we can impose the required phase to obtain the desired beam shape.
[0074] When compared to the conventional lens design methods, the IASA-based design method maximizes the power of the beam while producing an arbitrary pressure distribution in the plane of interest. Furthermore, the method can be extended to constrain the amplitude distribution in several different planes of propagation. Analogously, the method can be extended to produce different beam patterns using ultrasound transducers at different frequencies. The method can also be used to constrain the phase distribution in one or more planes, or both amplitude and phase distributions simultaneously. The desired target field may be binary or continuously varying in amplitude and/or phase over the focal plane of interest.
[0075] The IASA method can also be used with different transducer geometries. For instance, for a focused transducer, the exact pressure field can be simulated and verified through holographic scanning in a plane. Then the pressure field at the transducer aperture (obtained by back-projection) is used as the initial boundary condition over which we can impose the required phase to obtain the desired beam shape.
[0076] The above-described method for defining the thickness and shape of the lens uses the IASA. However, in different embodiments, other iterative methods for defining the thickness and shape of the lens are also possible. Sample results obtained with the IASA-designed lenses are described below.
[0077]
[0078] In the illustrated embodiment, the pressure amplitude field (i.e., distribution) that is generated at 32 mm distance from the source is shown in
[0079] In some embodiments, the customizable lens produces phase patterns, whereby the phase of the propagating wave front varies with position on the focal surface. This phase front of the propagating wave may permit non-invasive repositioning of a target located on a focal surface. Such targets may include kidney stones, bladder stones, calcifications, and other endogenous materials lodged in an anatomical vessel. In some embodiments, the customizable lenses create a pressure well around the target, pushing the target toward an area of relative low pressure.
[0080]
[0081]
[0082]
[0083] In some embodiments, a time-varying signal alters the beam pattern of a single-lensed transducer. For example, a customizable lens may generate multiple patterns at different frequencies simultaneously or may generate a single pressure pattern for a finite temporal period. In one embodiment, the frequency of a sinusoidal ultrasound signal may be varied over time to change the pattern, either continuously as a frequency chirp, or discretely in intervals. In another embodiment, a short signal pulse may be generated by the transducer to produce a temporary holographic image for a therapy such as shock wave lithotripsy, burst wave lithotripsy, or histotripsy. In other embodiments, the customizable lens may be designed to produce a target distribution of ultrasound phases that, for example, push the target in a desired direction.
[0084] Comparison of
[0085]
[0086]
[0087]
[0088] In some embodiments, the simulated pressure amplitude and phase distributions are obtainable using a single lens designed by the IASA-based methods. In general, the match is good between the target and simulated distribution for both the amplitudes and phases, and at both distances of interest. Therefore, it is possible to obtain different distributions of different parameters (e.g., pressure amplitude and phase) at different target distances from the source of ultrasound.
[0089] As explained above, the customizable lens can be designed based on the target amplitude/phase distributions that are shown in
[0090] In some embodiments, the transducer may be a phased array having transducer elements that are electronically controlled to generate the amplitude and/or phase at proper frequency.
[0091] The illustrated phased array elements 12-i are arranged in a plane, but, in other embodiments, the elements 12-i may be arranged along a curved surface. For example, the phased array elements 12-i may be angled towards a focus, and may be activated at the phase that takes into account this spatial distribution of the phased array elements 12-i. In some embodiments, such a transducer is more efficient than a planar transducer with a lens. Illustrative results obtained with a phased array are discussed with reference to
[0092]
[0093]
[0094]
[0095] As shown in
[0096]
[0097] In general, many of the sample pressure or phase distributions shown in
[0098]
[0099] In some embodiments, the holding mechanism 170 is a quick-change clamp assembly, including a clamp 170-1, anchored to one location of the transducer 12-T, and a quick-change bracket 170-2, fixed to the customizable lens 160 at another location, for example, opposite to the location of the clamp 170-1. In some embodiments, the clamp 170-1 is attached to the customizable lens 160. The clamp 170-1 may include a recess or a receptacle for inserting one end of the customizable lens 160, after the temporary interface 140 has been mounted to the customizable lens 160. Once aligned, the quick-change 170-2 may fit the customizable lens 160 and temporary interface 140 conformably to the surface of the transducer 12-T.
[0100] In some embodiments, the holding mechanism 170 is a threaded retaining ring that reversibly mounts the customizable lens 160 and the temporary interface 140 to the transducer 12-T using a threaded junction. In some embodiments, the holding mechanism 170 includes a plurality of clip-in fasteners attached in part to the customizable lens 160.
[0101]
[0102] The receive transducer 12-R may be placed at an oblique angle relative to the direction of the propagating wave front 160-T emitted by the transmit transducer 12-T. In some embodiments, the receive transducer 12-R converts ultrasound waves into an electronic signal and sends the signal to the controller 40 for further processing. In some embodiments, the receive transducer 12-R is a sensor such as an ultrasonic microphone, a laser interferometer, etc. In other embodiments, the receiver transducer 12-R may have similar structure as the transmit transducer 12-T, except for being configured to receive and process the ultrasound, and not to transmit the ultrasound. In some embodiments, the transducer 12-T may fulfil both transmit and receive functions.
[0103] In some embodiments, the controller 40 provides information about the condition and position of the target 22. For example, the propagating wave front 160-T may fragment or move the target 22, whereupon the receiver transducer 12-R may measure the change in the target 22 based on the reflected wave front 160-R, and then provide data to the controller 40.
[0104] In some embodiments, the receiver 12-R is, e.g., Computed Tomography (CT), magnetic resonance imaging (MRI) or other imaging system. Based on the ultrasound, CT, and/or MRI imaging, a 3D reconstruction of the stone may be obtained. In operation, the image may be used for designing the customizable lens and/or for monitoring the treatment process. In some embodiments, the customizable lens may be further modified by, for example, machining, based on the observed progress of the ultrasound treatment of the target 22.
[0105] Many embodiments of the technology described above may take the form of computer- or controller-executable instructions, including routines executed by a programmable computer or controller. Those skilled in the relevant art will appreciate that the technology can be practiced on computer/controller systems other than those shown and described above. The technology can be embodied in a special-purpose computer, controller or data processor that is specifically programmed, configured or constructed to perform one or more of the computer-executable instructions described above. Accordingly, the terms “computer” and “controller” as generally used herein refer to any data processor and can include Internet appliances and hand-held devices (including palm-top computers, wearable computers, cellular or mobile phones, multi-processor systems, processor-based or programmable consumer electronics, network computers, mini computers and the like).
[0106] From the foregoing, it will be appreciated that specific embodiments of the technology have been described herein for purposes of illustration, but that various modifications may be made without deviating from the disclosure. Moreover, while various advantages and features associated with certain embodiments have been described above in the context of those embodiments, other embodiments may also exhibit such advantages and/or features, and not all embodiments need necessarily exhibit such advantages and/or features to fall within the scope of the technology. Accordingly, the disclosure can encompass other embodiments not expressly shown or described herein.