Motion compensated image-guided focused ultrasound therapy system

Abstract
An image-guide therapy system comprises a thermal treatment device (e.g., an ultrasound transducer) configured for transmitting a therapeutic energy beam, and The system further comprises an imaging device (e.g., a magnetic resonant imaging (MRI) device) configured for acquiring images of the target tissue mass and the thermal treatment device. The system further comprises a controller configured for controlling thermal dose properties of the thermal treatment device to focus the energy beam on a target tissue mass located in an internal body region of a patient, and a processor configured for tracking respective positions of the thermal treatment device and the target tissue mass in a common coordinate system based on the acquired images. The system may optionally comprise a display configured for displaying the acquired images.
Description
FIELD OF THE INVENTION

The present invention relates to thermal treatment systems, and more particularly, to image-guided focused ultrasound therapy systems.


BACKGROUND OF THE INVENTION

Thermal energy, such as generated by high intensity focused ultrasound (acoustic waves with a frequency greater than about 20 KHz, and more typically between 50 KHz and 5 MHz), may be used to therapeutically treat internal tissue within a patient. For example, ultrasonic waves may be used to ablate tumors (e.g., breast tumors and uterine fibroids), thereby obviating the need for invasive surgery. For this purpose, a phased transducer array having transducer elements driven by electrical signals to produce ultrasonic energy can be placed external to the patient, but in close proximity to the target tissue mass to be ablated.


The transducer is geometrically shaped and positioned, such that the ultrasonic energy is focused at a “focal zone” corresponding to the target tissue mass within the patient. During the wave propagation across tissue, a portion of the ultrasound energy is absorbed, leading to increased temperature and eventually to cellular necrosis, preferably at the target tissue mass at the focal zone. The focal zone of the transducer can be rapidly displaced by independently adjusting the amplitude and phase of the electrical signal input to each of the transducer elements. The focal size and focal length of the ultrasound transducer will depend on the ultrasound frequency, the focal depth and the aperture size of the transducer. Because the size of target tissue mass is often greater than the size of the focal zone, the transducer may be sequentially focused and activated at a number of target sites within the target tissue mass to fully coagulate the target volume. The sequential “sonications” are used to cause coagulation necrosis of an entire tissue structure, such as a tumor, of a desired size and shape.


Image-guided focused ultrasound therapy systems offer the benefit of target visualization and localization. In particular, before a focused ultrasound treatment procedure is performed, a patient may be initially imaged to localize the mass and/or to plan a trajectory of the ultrasound beam. For example, using displayed images of the internal body region, a treatment boundary can be defined around the target tissue mass, and obstacle boundaries can be defined around tissue that should not be exposed to the ultrasound energy beam. The ultrasound transducer can then be operated based on these defined boundaries. During treatment, the patient can be continuously imaged to ensure that the target tissue mass is treated without damaging surrounding healthy tissue.


Magnetic Resonance Imaging (MRI) guidance offers the additional benefit of temperature mapping in vivo, which can be used to verify that a sufficient temperature is reached during each application of ultrasonic energy (i.e., sonication) to kill the target tissue mass or portion thereof. Temperature mapping can be accomplished by measuring the temperature change (rise) of the portion of the tissue mass being heated during each sonication using conventional MR imaging techniques coupled with image processing to extract the temperature from the MRI data. Thus, accurate temperature measurements allow verification of the proper location of the focal zone and computation of the accumulated thermal dose during treatment for prediction of tissue ablation.


Current image-guided focused ultrasound therapy systems, such as those based on MRI, assume that the acoustic transducer is in a predefined and known position relative to the target tissue mass to be treated. In this case, the respective coordinate systems defined by the imaging device, therapeutic device, and patient remain registered with each other. If unplanned movement of either the transducer or the patient is detected, however, one or both of these coordinate systems become mis-registered with the other coordinate systems. As such, the treatment process must be stopped and the ultrasound beam trajectory re-planned. This introduces significant inefficiencies in the treatment process and may generate significant delays. Some image-guided focused ultrasound therapy systems use mechanically aligned imaging and ultrasound transducer arrangements or use the same transducer to perform both imaging and therapy tasks. In these cases, the coordinates of the therapy system remain registered with the coordinates of the imaging system if the transducer moves. However, if the patient moves, the patient coordinate system in which the target tissue region is defined, will become mis-registered with the imaging and therapeutic coordinate systems.


Current generation of systems require that all three coordinate systems—patient, imaging, and therapeutic, be locked or be continually registered with respect to each other, so that treatment of the target tissue region without damaging surrounding healthy tissue can be assured and verified. Conventionally, this has been accomplished by mechanically aligning the imaging and therapy transducers, as briefly discussed above, and immobilizing the target tissue region. However, electronic phasing changes or electrical impedance changes within the control circuitry of the transducer may still generate misalignment between the imaging and therapy beams; that is, cause mis-registration between the imaging and therapeutic coordinate systems. In addition, immobilization of the target tissue region of the patient may not always be practically accomplished.


For example, after delivery of a thermal dose, e.g., an ultrasound sonication, a cooling period is required to avoid harmful and painful heat build up in healthy tissue adjacent the target tissue mass. This cooling period may be significantly longer than the thermal dosing period. Since a large number of sonications may be required in order to fully treat the target tissue mass, the overall time required can be quite significant. This means that the patient must remain motionless in the imaging device for a significant period of time, which can be very stressful. At the same time, it may be critical that the entire target tissue mass be ablated (such as, e.g., in the case of a malignant cancer tumor), and that no short cuts be taken during the procedure just in the name of patient comfort. Thus, the use of image-guided focused therapy systems are limited to treating tissue masses with small motion amplitudes or to those that are easily immobilized.


Recently, it has been suggested that the movement of a target tissue mass relative to an MRI device and ultrasound transducer can be compensated for when the motion is periodic (e.g., motion caused by a physiological cycle, such as a cardiac or respiratory cycle). See Denis de Senneville, B., Mougenot, C., and Moonen, C., Real-Time Adaptive Methods for Treatment of Mobile Organs by MRI-Controlled High-Intensity Focused Ultrasound, Magnetic Resonance in Medicine 57:319-330 (2007). In particular, an MRI of a target tissue mass currently acquired during a treatment procedure can be compared to reference MRI images of the target tissue mass previously acquired during a periodic cycle before the treatment procedure, which information can then be used to predict the displacement of the target tissue mass relative to the MRI device and ultrasound transducer. The ultrasound beam generated by the transducer can then be electronically controlled in real-time to continually maintain its focal zone at the moving target tissue mass. While this technique may be successful when tracking a target tissue mass that moves in accordance with a periodic cycle, it does not address non-periodic movement of the target tissue mass; for example, when the patient moves within the imaging device.


There, thus, remains a need for an improved method and system that can correlate relative displacement between an imaging device, a therapeutic device, and a target tissue mass in real time.


SUMMARY OF THE INVENTION

In accordance with a first aspect of the present inventions, an image-guide therapy system is provided. The system comprises a thermal treatment device configured for transmitting a therapeutic energy beam. In one embodiment, the thermal treatment device is an external thermal treatment device configured for transcutaneously transmitting the energy beam into a patient. In another embodiment, the thermal treatment device is an ultrasound transducer, in which case, the energy beam will be an ultrasound energy beam. For example, the ultrasound transducer may comprise a multiplicity of independently controlled transducer elements, and the thermal dose properties may be controlled by adjusting the one or more of a phase, frequency, and amplitude of respective drive signals supplied to the respective transducer elements. The system further comprises an imaging device (e.g., a magnetic resonance imaging (MRI) device) configured for acquiring images of the target tissue mass and the thermal treatment device. In one embodiment, the thermal treatment device and imaging device are different devices, although the thermal treatment device and imaging device can be the same device.


The system further comprises a controller configured for controlling thermal dose properties of the thermal treatment device to focus the energy beam on a target tissue mass located in an internal body region of a patient, and a processor configured for tracking respective positions of the thermal treatment device and the target tissue mass in a common coordinate system (e.g., a three-dimensional coordinate system) based on the acquired images. The system may optionally comprise a display configured for displaying the acquired images.


In one embodiment, the processor is configured for tracking the respective positions of the thermal treatment device and the target tissue mass during non-periodic relative movement between the thermal treatment device and the target tissue mass. The respective positions of the thermal treatment device and the target tissue mass may be tracked in any one of a variety of manners. For example, the processor may be configured for tracking the respective positions of the thermal treatment device and the target tissue mass by analyzing fiducial markers associated with one or both of the thermal treatment device and the target tissue mass within the acquired images, by transmitting or receiving signals to or from a localization element associated with one or both of the thermal treatment device and the target tissue mass, or by comparing the acquired images with one or more reference images.


The tracked positions of the thermal treatment device and the target tissue mass can be used in any one of a variety of manners. For example, the controller may be configured for automatically adjusting one or more of the thermal dose properties of the thermal treatment device (e.g., by adjusting electrical and/or mechanical parameters of the thermal treatment device) based on the tracked positions of the thermal treatment device and the target tissue mass to compensate for a change in the relative position between the thermal treatment device and the target tissue mass. In this example, the system may optionally comprise a user interface configured for allowing a system operator to define a treatment boundary of the target tissue mass using the acquired images, in which case, the controller may be configured for automatically controlling the thermal dose properties of the thermal treatment device to focus the energy beam within the defined treatment boundary of the target tissue mass. The user interface may also be configured for allowing the system operator to define an obstacle boundary of tissue, in which case, the controller can be further configured for controlling the thermal dose properties of the thermal treatment device to avoid traversing the defined obstacle boundary with the energy beam. As another example, the processor may be configured for tracking treated regions and untreated regions of the target tissue mass based on the tracked positions of the thermal treatment device and the target tissue mass.


In accordance with a second aspect of the present inventions, another image-guide therapy system is provided. The system comprises a thermal treatment device configured for transmitting a therapeutic energy beam having a focal zone, and an imaging device configured for acquiring images of the target tissue mass. The thermal treatment device and imaging device may have the same features described above. The system further comprises a controller configured for controlling thermal dose properties of the thermal treatment device to locate the focal zone of the energy beam coincident with a target tissue mass located in an internal body region of a patient, and a processor configured for tracking relative non-periodic movement between the focal zone of the energy beam and the target tissue mass in real-time based on the acquired images. The system may optionally comprise a display configured for displaying the acquired images.


In one embodiment, the relative non-periodic movement between the focal zone of the energy beam and the target tissue mass is tracked in a common three-dimensional coordinate system. In another embodiment, the processor is configured for tracking the relative non-periodic movement between the focal zone of the energy and the target tissue mass in a common three-dimensional coordinate system. The relative non-periodic movement between the focal zone of the energy beam and the target tissue mass can be tracked in any one of a variety of manners. For example, the respective positions of the thermal treatment device and the target tissue mass may be tracked during non-periodic relative movement between the thermal treatment device and the target tissue mass using any one of the exemplary techniques described above, or if the acquired images are thermally sensitive images, the thermal signature of the focal zone of the energy beam within the thermally sensitive images can be analyzed.


The tracked relative non-periodic movement between the thermal treatment device and the target tissue mass can be used in any one of a variety of manners, including those described above. For example, the controller may be configured for automatically adjusting one or more of the thermal dose properties of the thermal treatment device (e.g., by adjusting electrical and/or mechanical parameters of the thermal treatment device) based on the tracked relative non-periodic movement between the focal zone of the energy beam and the target tissue mass to compensate for a change in the relative position between the thermal treatment device and the target tissue mass. In this example, the system may optionally comprise a user interface having the same functions described above. As another example, the processor may be configured for tracking treated regions and untreated regions of the target tissue mass based on the tracked relative non-periodic movement between the thermal treatment device and the target tissue mass.


Other and further aspects and features of the invention will be evident from reading the following detailed description of the preferred embodiments, which are intended to illustrate, not limit, the invention.





BRIEF DESCRIPTION OF THE DRAWINGS

The drawings illustrate the design and utility of preferred embodiments of the present invention, in which similar elements are referred to by common reference numerals. In order to better appreciate how the above-recited and other advantages and objects of the present inventions are obtained, a more particular description of the present inventions briefly described above will be rendered by reference to specific embodiments thereof, which are illustrated in the accompanying drawings.


Understanding that these drawings depict only typical embodiments of the invention and are not therefore to be considered limiting of its scope, the invention will be described and explained with additional specificity and detail through the use of the accompanying drawings in which:



FIG. 1 is a block diagram of an image-guided therapy system constructed in accordance with one embodiment of the present inventions;



FIG. 2 is a plan view of an ultrasound transducer used in the system of FIG. 1 to treat a target tissue mass located in the interior region of a patient;



FIG. 3 is a cross-section view of the ultrasound transducer of FIG. 2 being controlled by electrical drive signals;



FIG. 4A is a top view of one embodiment of the ultrasound transducer of FIG. 2



FIG. 4B is a top view of an alternative embodiment of the ultrasound transducer of FIG. 2;



FIG. 5 is a cross-sectional view of a tissue mass treated with a series of sonications performed by the system of FIG. 1;



FIGS. 6A-6C are two-dimensional pixel representations of untreated and treated regions of a target tissue mass during a series of sonications performed by the system of FIG. 1; and



FIG. 7 is a flow diagram illustrating a method of operating the system of FIG. 1 to treat a target tissue mass within the interior region of a patient.





DETAILED DESCRIPTION OF THE EMBODIMENTS

Referring generally to FIG. 1, an image-guided therapy system 10 arranged in accordance with one embodiment of the present inventions will now be described. The system 10 is designed to be operated by a system operator 12 to treat a target tissue mass 16 (e.g., a tumor) within an internal body region of a patient 14. The system 10 generally comprises a focused thermal treatment subsystem 18, a patient table 20, an imaging subsystem 22, a planner 24, a user interface 26, and a motion compensation processor 28. It should be noted that the elements illustrated in FIG. 1 are functional in nature, and are not meant to limit the structure that performs these functions in any manner. For example, several of the functional blocks can be embodied in a single device, or one of the functional blocks can be embodied in multiple devices. Also, the functions can be performed in hardware, software, or firmware.


Referring further to FIG. 2, the thermal treatment system 12 comprises a thermal treatment device 30 configured for delivering an energy beam 32 to the target tissue mass 16 within the patient 14. The energy beam 32 raises the temperature of the target tissue mass 16 by transferring the power of the energy beam as heat to the target tissue mass 16. The energy beam 32 is focused on the target tissue mass 16 in order to raise the temperature of the target tissue mass 16 to a point where it is destroyed. The heat distribution within the tissue is controlled by the energy density within a focal zone 34 of the energy beam 32, the tissue acoustic parameters and, to some extent, the duration of the energy beam application. In the illustrated embodiment, the thermal treatment device 30 is configured for being placed external to the patient 14, in which case, the energy beam 32 will be transcutaneously transmitted into the patient 14, so that the energy beam 32 is focused on the target tissue mass 16 some distance from a skin surface 36 of the patient 14. The distance from the skin surface 36 to the target tissue mass 16 is the near field, which contains healthy tissue. It is important that tissue in the near field is not damaged by the energy beam 32. Thus, the energy beam 32 is preferably focused within the target zone, where the energy is transferred as heat to the target tissue mass 16.


In one implementation illustrated in FIGS. 3 and 4A, the thermal treatment device 30 takes the form of a phased array ultrasound transducer, in which case, the energy beam 32 takes the form of an ultrasound energy beam. In one embodiment, the transducer 30 may have a concave or bowl shape, such as a “spherical cap” shape; that is, having a substantially constant radius of curvature, such that the transducer 30 has an inside surface 38 defining a portion of a sphere. Alternatively, as illustrated in FIG. 4B, the transducer 30 may have a substantially flat configuration and/or may include an outer perimeter that is generally rectangular. It should be appreciated that a variety of geometric designs for the transducer 30 may be employed. Additionally, alternative embodiments of the system 10 may use focused radiators, acoustic lenses, or acoustic reflectors in order to achieve optimal focus of the beam 32. In a preferred embodiment, the transducer 30 has an outer diameter of between about two and sixteen centimeters.


The ultrasound transducer 30 converts an electronic drive signal (described in further detail below) into acoustic power in the ultrasound beam 32. In particular, the ultrasound transducer 30 comprises a plurality of transducer elements 40, each composed of piezoelectric material or other materials like silicon-based transducers, such that, upon being driven with a sinusoidal signal near the resonant frequency of the piezoelectric material or silicon transducers, the transducer elements 40 vibrate according to the phase and amplitude of the exciting sinusoidal signal, thereby creating the desired ultrasound energy beam 32.


The exemplary embodiment of the transducer 30 is divided radially into six concentric rings 42 and circumferentially into eight sectors 44 disposed about a central axis 46, thereby dividing the transducer 30 into forty eight transducer elements 40. Alternatively, the transducer 30 may be divided into any desired number of rings 42 and/or sectors 44. In the illustrated embodiment, each transducer element 40 has an arcuate shape, although each transducer element 40 may have any one or more of a variety of geometric shapes, such as hexagons, triangles, squares, and the like. In the illustration FIG. 4B, a different implementation is described in which the array is flat with very large number of elements isotropically spread over the surface. The configuration of the transducer 14, however, is not important to the present invention, and any of a variety of known ultrasound transducers may be used, such as flat circular arrays, linear arrays, and the like. Additional information on the construction of an ultrasound transducer appropriate for use with the system 10 may be found, for example, in Cain, C. and Umemura, S., “Concentric-Ring and Sector-Vortex Phased-Array Applicators for Ultrasound Hyperthermia,” IEEE Transactions on Microwave Theory and Techniques, Vol. MTT-34, No. 5, pages 542-551 (May 1986); and Fjield, T. and Hynyen, K., “The Combined Concentric-Ring and Sector-Vortex Phased Array for MRI Guided Ultrasound Surgery, IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, Vol. 44, No. 5, pages 1157-1167 (September 1997), the disclosures of which are expressly incorporated herein by reference.


The thermal treatment subsystem 18 further comprises an electrical driver 48 to which each of the transducer elements 40 is individually coupled in a conventional manner. The electrical driver 48 is configured for providing electrical drive signals 50 to the transducer elements 40 at one or more frequencies, preferably at radio frequencies (RF), for example, from 0.25 MHz to 10 MHz, and more preferably from 0.5 MHz to 3.0 MHz. When electrical drive signals 50 are provided to the transducer elements 40, the transducer 30 emits ultrasonic energy from its inner surface 38 as is known to those skilled in the art.


Referring specifically now to FIG. 1, the thermal treatment subsystem 18 further comprises a mechanical driver 52 (or positioner) that is coupled to the transducer 30. In one embodiment, the mechanical driver 52 may be operated to translate with the transducer 30 in a plane that is perpendicular to the axis of the transducer 30, as well as controlling the roll and pitch of the transducer 30. A preferred mechanical driver (or positioner) is disclosed in U.S. Pat. No. 6,582,381, entitled “Mechanical Positioner for MRI Guided Ultrasound Therapy System,” which is expressly incorporated herein by reference.


The thermal treatment subsystem 18 further comprises a focus ultrasound (FUS) controller 54 that controls the thermal dose properties of the transducer 30 via the electrical driver 48 and the mechanical driver 52. Thermal dose properties may include the duration and power of the ultrasound beam 32, the frequency of the acoustic beam, and the position and size of the focal zone 34 of the ultrasound beam 32. In order to control the thermal dose properties of the transducer 30, the FUS controller 54 is coupled to the electrical driver 48 to dictate the electrical parameters (e.g., amplitude, frequency, and/or phase and duration) of the drive signals 50 (shown in FIG. 3), and thus the electrical parameters of the ultrasound energy emitted by the respective transducer elements 40. For example, the FUS controller 54 may control the amplitude of the drive signals 50 to control the intensity of ultrasonic energy delivered by the transducer 30, and the phase and frequency of the drive signals 50 to control the location, shape and size of the focal zone 34 and/or to adjust the focal distance (which in essence is a subset of location) of the focal plane (i.e., the distance from the face of the transducer 30 to the center of the focal zone 34). The frequency of the drive signals 50 can also be adjusted to optimize the energy density within the focal zone 34. To provide for precise positioning and dynamic movement and reshaping of the focal zone 34, it is desirable to be able to adjust the phase and/or amplitude of the individual transducer elements 40 relatively quickly, e.g., in the microsecond range. In order to further control the thermal dose properties of the transducer 30, and in particular, the position of the focal zone 34, the FUS controller 54 is also coupled to the mechanical driver 52, so that the mechanical position of the transducer 30 can be adjusted to displace the focal zone 34 to a desired location.


More advanced techniques for obtaining specific focal distances and shapes are disclosed in U.S. Pat. Nos. 6,506,171, entitled “Systems and Methods for Controlling Distribution of Acoustic Energy Around a Focal Point Using a Focused Ultrasound System,” U.S. Pat. No. 6,419,648, entitled “Systems and Methods for Reducing Secondary Hot Sports in a Phase Array Focused Ultrasound System,” and U.S. Pat. No. 6,613,004, entitled “Systems and Methods for Creating Longer Necrosed Volumes Using a Phased Array Focused Ultrasound System,” all of which are expressly incorporated herein by reference.


Significantly, as will be described in further detail below, the FUS controller 54, based on information supplied by the motion compensation processor 28, automatically controls the electrical and mechanical properties of the transducer 30 based on a tracked relative movement between the focal zone 34 of the ultrasound beam 32 and the target tissue mass 16 to maintain the location of the focal zone 34 of the ultrasound energy beam 32 inside the target tissue mass 16. In one embodiment, the FUS controller 54 automatically controls the electrical and mechanical properties of the transducer 30 based on tracked positions of the transducer 30 and the target tissue mass 16 to compensate for any change in the relative position between the therapy device and the target tissue mass.


Referring still to FIG. 1, the patient table 56 comprises a chamber 58 filled with degassed water or similar acoustically transmitting fluid in which the transducer 30 is disposed. Alternatively, the chamber 58 may be disposed in a fluid-filled bag mounted on a movable arm that may be placed against the patient's body (not shown). The patient table 56 further comprises a flexible membrane 60 on which the patient 14 may be disposed. That membrane 60 is composed of a material that is substantially transparent to ultrasound, such as Mylar®, polyvinyl chloride (PVC), or other suitable plastic material. A fluid-filled bag (not shown) may be provided on the membrane 60 that may conform easily to the contours of the patient 14 disposed on the table 56, thereby acoustically coupling the patient 14 to the transducer 30 within the chamber 58. In addition or alternatively, an acoustic gel, water, or other fluid (not shown) may be provided between the patient 14 and the membrane 60 to facilitate further acoustic coupling.


In the illustrated embodiment, the imaging subsystem 22 is configured for acquiring three-dimensional images of the target tissue mass 16 at a fast rate (i.e., in real-time. Notably, whether images are acquired in real-time will depend on the velocity at which the target tissue mass 16 is expected to move and the error that can be tolerated when tracking the movement of the target tissue mass 16, as will be described in further detail below. For example, if is expected that the target tissue mass 16 may move at 10 cm/sec at a maximum, and the tracking error is less than 1 mm, than it is preferred that the imaging subsystem 22 generate at least 100 images per second. If the target tissue mass 16 moves in accordance with a periodic cycle (e.g., liver motion caused by a cardiac cycle may be 4 cm/sec), the movement of the target tissue mass 16 is more predictable, whereas if the target tissue mass 16 moves in accordance with a non-periodic cycle (e.g., caused by general movement of the patient 14), the movement of the target tissue mass 16 may not be as predictable. Preferably, the rate at which the imaging subsystem 22 acquires images is at least 5 images/second, more preferably, greater than 10 images/second, and most preferably greater than 50 images/second.


In the illustrated embodiment, the imaging subsystem 22 takes the form of a Magnetic Resonance Imaging (MRI) subsystem, which generally comprises a high field magnet 62, a gradient field amplifier 64, a radio frequency (RF) transmitter 66, an RF receiver 68, an MRI controller 70, an MRI processor 72, and a display 74.


The magnet 62 may be integrated with the patient table 56 and includes a region for receiving the patient 14 therein. The magnet 62 provides a static, relatively homogenous magnetic field over the patient 14, as is well known in the art. The gradient field amplifier 64 generates magnetic field gradients that vary the static magnetic field generated by the magnet 62 in a known manner. The RF transmitter 66 generates and transmits RF pulse sequences or other signals over the patient 14 to cause the target tissue mass 16 to emit RF response signals, which may include free induction decay (FID) signals and/or echo signals. In a preferred embodiment, the RF transmitter 66 includes RF coils (not shown) in the magnet 62, and a pulse transmitter (also not shown), which may have a pulse transmitter frequency supplied by a synthesizer (not shown) and/or controlled by the MRI controller 70. The RF receiver 68 senses raw MR response signals. The RF receiver 68 may include a separate set of RF coils (not shown), from the RF transmitter 66. Alternatively, an RF transmitter/receiver (not shown) may be provided that is configured to operate alternately in a transmit mode to transmit RF pulse sequences and a receive mode to sense the MR response signals. Additional receive only or receive/transmit MR imaging coils (not shown) may be placed adjacent the target tissue volume to accomplish a better localized imaging quality.


The MRI controller 70 provides or otherwise controls the timing sequence used to operate the MRI subsystem 22. For example, the timing sequence may include one or more signals instructing the RF transmitter 66 to transmit RF pulse sequences, and/or instructing the RF receiver 68 to listen for MR response signals. For the alternative embodiment wherein an RF transmitter/receiver is provided, the MRI controller 70 may control this operation, for example, by switching the RF coils of the RF transmitter/receiver between the transmit and receive modes. The MRI processor 72, which may include an analog-to-digital converter and/or image processor (both not shown) receives the MR response signals from the RF receiver 68 and generates MR images therefrom in a conventional manner for display on the display 74.


In addition or alternatively, the MRI processor 72 may rapidly acquire temperature-sensitive images of the patient 14. Such temperature-sensitive images may be superimposed on other medical diagnostic images (such as conventional MR images) or provided on separate displays. Preferably, both the target tissue mass 16 and regions heated by the focused ultrasound subsystem 22 may be imaged simultaneously, thereby enabling an operator to verify that the heated region (i.e., the focal zone 34) corresponds to the target tissue mass 16. Additional information on systems and methods for obtaining temperature-sensitive MR images may be found in U.S. Pat. No. 6,559,644, entitled “MRI-Based Temperature Mapping With Error Resolution,” and U.S. Pat. No. 6,618,608, entitled “Thermal Imaging of Fat and Muscle Using a Simultaneous Phase and Magnitude Double Echo Sequence,” which are expressly incorporated herein by reference.


Notably, in order to quickly adjust the location of the focal zone 34 of the ultrasound energy beam 32 relative to the target tissue mass 16 in response to movements made by internal organs, such as the liver, kidney, spleen, heart, etc., it is desirable that the MR images generated by the MR imager 22 be generated as quickly as possible. To accomplish this, the MR imager 22 preferably acquires images of the target tissue mass 16 using a relatively small scanning window. As such, without adjusting the MR imager 22, the moving target tissue mass 16 may move out of the scanning window of the MR imager 22. To compensate for this, the motion compensation processor 28, which will be described in further detail below, transmits control signals to the MRI controller 70 to adjust the location of the scanning window to correspond with the location of the target tissue mass 16.


The planner 24 automatically constructs a treatment plan, which consists of a series of treatment sites represented by thermal dose properties. The purpose of the treatment plan is to ensure complete ablation of the target tissue mass 16 by planning a series of sonications (i.e., distinct applications of the ultrasound beam 32) that will apply a series of thermal doses 76 at various points within the target tissue mass 16, resulting in a composite thermal dose sufficient to treat the entire target tissue mass 16, as illustrated in FIG. 5. The plan will include the afore-described thermal dose properties determined based on default parameters, such as a thermal dose threshold, maximum allowed energy for each thermal dose, thermal dose duration for each treatment site, cooling time between ultrasound applications (sonications), sonication grid density (how much the sonications should overlap), the physical parameters of the transducer 30, the patient anatomy and anatomy constraints.


In order to construct the treatment plan, the planner 24 receives input from the user interface 26. For example, in one implementation, a user specifies the clinical application protocol, i.e., breast, pelvis, eye, prostrate, etc., via the user interface 26. Selection of the clinical application protocol may control at least some of the default parameters. In other implementations, some or all of these parameters are input through the user interface 26 as user specified parameters. Additionally, a user may edit any of the default parameters via the user interface 26. In one implementation, the user interface 26 comprises a Graphical User Interface (GUI). In this case, a user employs a mouse or touch screen to navigate through menus or choices as displayed on the display 74 or another display in order to make the appropriate selections and supply the required information.


To further facilitate construction of the treatment plan, the planner 24 uses images supplied by the MRI subsystem 22 in conjunction with input supplied by the user interface 26. In one implementation, the images are used to define a treatment boundary around the target tissue mass 16, e.g., by tracing a line on the image displayed on the display 74. In the case where the images are three-dimensional, the treatment boundary can be defined in three dimensions. The images may also be used to define obstacle boundaries around tissue that is acoustically reflective (such as gas) or bone and tissue that is otherwise sensitive to thermal exposure.


There are several methods that are used to change or update the treatment plan. For example, at the end of each sonication, there may be regions within the target tissue mass 16 that are not covered by accumulated treatment. After each sonication, the treated regions and untreated regions may be taken into account in generating an updated treatment plan. In order to accomplish the tracking of untreated regions, each target tissue mass 16 may be maintained as a two-dimensional linked list of pixel ranges sorted by (y) and then (x) coordinates, as illustrated in FIGS. 6A-6C. Notably, if the target tissue mass 16 is to be represented as a three-dimensional mass, the target tissue mass 16 may be maintained as a three-dimensional linked list of voxel ranges.


As can be seen in FIG. 6A, the target tissue mass 16 is represented as a continuous region 78 of lighter pixels. The pixel distribution is then represented by the data structure 80. This type of representation is called “run-length encoding.” The data structure 80 contains linked lists 82 by row indicating the pixels that contain the target tissue mass 16. Thus, for row 1, a list element 82a indicates that the pixel range from 1 to 2 contains a portion of the target tissue mass 16, a list element 82b indicates that the pixel range from 7 to 7 contains a portion of the target tissue mass 16, and a list element 82c indicates that the pixel range from 11 to 11 contains a portion of the target tissue mass 16. For row 2, a list element 82d indicates that the pixel range from 1 to 3 contains a portion of the target tissue mass 16, and a list element 82e indicates that the pixel range from 7 to 11 contains a portion of the target tissue mass 16. For row 3, a list element 82f indicates that the pixel range from 1 to 11 contains a portion of the target tissue mass 16. For row 4, a list element 82g indicates that the pixel range from 1 to 10 contains a portion of the target tissue mass 16. It can also be seen that rows 0 and 5 do not contain any portion of the target tissue mass 16. Therefore, these rows in the data structure 80 do not contain any pixel ranges.


Once a sonication is applied, the area of the target tissue mass 16 that is destroyed (the treated region) is presented in the same fashion as the untreated region shown in FIG. 6A, and the treated region is subtracted from the untreated region in order to define the remaining untreated region within the target tissue mass 16. After several sonications, the target tissue mass 16 may be represented as a continuous region 78 of lighter pixels, as illustrated in FIG. 6B. In this case, for row 1, a list element 82h indicates that the pixel range from 1 to 2 contains a portion of the target tissue mass 16, and a list element 82i indicates that the pixel range from 6 to 11 contains a portion of the target tissue mass 16. For row 2, a list element 82j indicates that the pixel range from 2 to 3 contains a portion of the target tissue mass 16, a list element 82k indicates that the pixel range from 7 to 7 contains a portion of the target tissue mass 16, and a list element 82l indicates that the pixel range from 11 to 11 contains a portion of the target tissue mass 16. For row 3, a list element 82m indicates that the pixel range from 1 to 11 contains a portion of the target tissue mass 16. For row 4, a list element 82n indicates that the pixel range from 4 to 9 contains a portion of the target tissue mass 16. After several more sonications, the target tissue mass 16 may be represented as discontinuous regions 78 of lighter pixels, as illustrated in FIG. 6C. In this case, for row 2, a list element 82o indicates that the pixel range from 1 to 1 contains a portion of the target tissue mass 16, and a list element 82p indicates that the pixel range from 8 to 10 contains a portion of the target tissue mass 16. For row 4, a list element 82q indicates that the pixel range from 1 to 3 contains a portion of the target tissue mass 16, and a list element 82r indicates that the pixel range from 10 to 10 contains a portion of the target tissue mass 16. Further details discussing tracking untreated regions in a target tissue mass are described in U.S. Pat. No. 6,618,620, entitled “Focused Ultrasound System With MRI Synchronization,” which is expressly incorporated herein by reference.


The treatment plan or revised treatment plans are passed to the motion compensation processor 28, which calculates the values of the electrical and/or mechanical parameters used by the FUS controller 54 to control the thermal dose properties of the transducer 30 via the electrical driver 48 and the mechanical driver 52. In the case where a treatment boundary and obstacle boundaries are defined, the FUS controller 54 will focus the ultrasound energy beam within the defined treatment boundary of the target tissue mass and/or avoid traversing the defined obstacle boundaries with the ultrasound beam 32.


Significantly, the motion compensation processor 28 is configured for tracking relative movement between the focal zone 34 of the ultrasound energy beam 32 and the target tissue mass 16 in real-time based on the images acquired by the MRI subsystem 22. The relative movement tracked by the motion compensation processor 28 can be non-periodic (for example, random movement of the patient within the MRI subsystem 22). The relative movement tracked by the motion compensation processor 28 can also be periodic (for example, movement of an organ caused by a physiological parameter, such as a respiratory cycle and/or cardiac cycle). Notably, there is a direct correlation between the location of the focal zone 34 of the ultrasound beam 32 and the location of the transducer 30 within the imaging coordinate system. Thus, the motion compensation processor 28 can indirectly track the relative movement between the focal zone 34 of the ultrasound energy beam 32 relative to the target tissue mass 16 by tracking the respective positions of the transducer 30 and the target tissue mass 16 in a common three-dimensional coordinate system (and in particular, the imaging coordinate system) based on the acquired images.


In one implementation, the transducer 30 and the target tissue mass 16 can be located by comparing the currently acquired images of the transducer 30 and the target tissue mass 16 with reference images. For example, the transducer 30 can be located by performing image template matching between the MR image of the transducer 30 and a predefined graphical template that is scaled based on the field of view used, and inferring therefrom the coordinates of the transducer 30 within the image coordinate system. The target tissue mass 16 may be located by performing an image-based maximum correlation algorithm on the currently acquired image of the target tissue mass 16 and a reference image of the target tissue mass 16. In an alternative implementation, image comparison is not performed, but rather the transducer 30 can be located by placing localization elements 31 on the transducer 30. For example, the localization elements 31 may be fiducial elements, in which case artifacts representing the fiducial elements can be located within the MR images in a conventional manner to locate the transducer 30. As another example, the localization elements 31 can be transmitting/emitting/sensing elements that transmit and/or receive signals. By way of non-limiting examples, the transmitting/sensing elements 31 can be MR tracking devices, such as microcoils filled with MR sensitive material, or gradient-based tracking devices, such as microcoils sensitive to gradient magnetic field changes, such that the microcoils (and thus the transducer 30) can be located relative to the MRI subsystem 22. In any event, the determined coordinates of the transducer 30 and the target tissue mass 16 can equal any number, but in one implementation, the coordinates can number six (e.g., x-, y-, and z-rectilinear coordinates, and yaw, pitch, and roll coordinates) or a subset (e.g., four) of these coordinates.


While the motion compensation processor 28 has been described as tracking the position of the focal zone 34 of the ultrasound energy beam 32 indirectly by tracking the location of the transducer 30, the motion compensation processor 28 may alternatively or optionally track the position of the focal zone 34 of the ultrasound energy beam 32 relative to the target tissue mass 16 directly from the images. For example, the focal zone 34 of the ultrasound energy beam can be located by analyzing the optional thermal sensitive images and identifying the “hot spot” or thermal signature in the thermal sensitive image corresponding to the focal zone 34 relative to the target tissue mass 16 in the conventional MR images. In the case where the focal zone 34 is displaced in accordance with a physiological cycle by constructing an atlas of motion with the MR images acquired prior to the treatment period and comparing the respective MR images currently acquired during the treatment period to the images stored in the atlas. Further details discussing this technique are disclosed in Denis de Senneville, B., Mougenot, C., and Moonen, C., Real-Time Adaptive Methods for Treatment of Mobile Organs by MRI-Controlled High-Intensity Focused Ultrasound, Magnetic Resonance in Medicine 57:319-330 (2007), which is expressly incorporated herein by reference.


The motion compensation processor 28 uses the tracked relative movement between the focal zone 34 of the ultrasound energy beam 32 and the target tissue mass 16 to continually recalculate or transform the nominal values of the electrical and/or mechanical parameters previously calculated based on the treatment plan generated by the planner 24. These adjusted electrical and/or mechanical parameter values are passed to the FUS controller 54, which automatically adjusts the thermal dose properties of the transducer 30 to maintain the location of the focal zone 34 of the ultrasound energy beam 32 at the target tissue mass 16 or otherwise compensate for the change in the relative position between the transducer 30 and the target tissue mass 16.


As previously discussed, the thermal dose properties of the transducer 30, and in particular the location of the focal zone 34 of the ultrasound energy beam 32, can be changed via control of the electrical driver 48 and mechanical driver 52. Because the thermal dose properties of the transducer 30 should be modified as quickly as possible in order to compensate for relative movement between the transducer 30 and the target tissue mass 16 in real-time, it is preferred that the location of the focal zone 34 of the ultrasound energy beam 32 be adjusted, at least in part, via control of the electrical driver 48, which is faster than the mechanical driver 52.


As an alternative to, or as an adjunct to, compensating for the relative movement between the transducer 30 and the target tissue mass 16, the adjusted electrical and/or mechanical parameter values are passed to the planner 24, such that the treated and untreated regions of the target tissue mass 16 can be more accurately tracked. That is, by knowing the relative position between the focal zone 34 of the ultrasound energy beam 30 during each sonication, the linked lists 82 can be more accurately generated.


Having described the structure and function of the system 10, its operation in treating the target tissue mass 16 will now be described with further reference to FIG. 7. Initially, the system operator 12 selects an appropriate clinical application protocol via the user interface 26, after which the planner 24 or the FUS controller 54 will select or change the default parameters (step 100). After the clinical application is selected, relevant MR images of the target tissue mass 16 are acquired by the MRI subsystem 22 (step 102). Next, the system operator 12 uses the acquired MR images to define the treatment and obstacle boundaries via the user interface 26 (step 104). The system operator 12 may then enter additional default parameters or modify previously defined default parameters via the user interface 26 (step 106). The planner 24 then automatically constructs the treatment plan based on the defined treatment and obstacle boundaries and the default parameters (step 108). Next, the system operator 12 may edit the treatment plan via the user interface 26, e.g., by adding or deleting treatment boundaries and/or obstacle boundaries, changing the location of some or al of the treatment sites, or changing thermal dose properties (step 110). If the treatment plan is edited, the process returns to step 108 where the treatment plan is reconstructed/recalculated. Once the treatment plan is set, proper registration of the system 10 with regard to the position of the focal zone 34 relative to the target tissue mass 16 is verified (step 112). This verification step can comprise performing a low energy thermal dose at a predefined site within the target tissue mass 16 and generating a thermally sensitive image in order to ensure that the focal zone 34 is located within the target tissue mass 16. If the system 10 is not properly registered, the FUS controller 54 may be operated via the user interface 26 to adjust the mechanical position of the transducer 30 via the mechanical driver 52, which mechanical position can then be set as a “home” position (step 114). The process can then return to step 112 to again verify the registration of the system 10. Once verification is complete, the treatment plan is implemented (step 116).


During implementation step of the treatment plan, a sonication can be performed via operation of the focused ultrasound subsystem 18 while MRI images are acquired by the MRI subsystem 22. The MRI images may be temperature sensitive images sequences as each step of the treatment plan is implemented. These images will illustrate the actual thermal dose distribution resulting from each successive thermal dose. Significantly, during implementation of the treatment plan, the motion compensation processor 28 tracks the relative position between the focal zone 34 of the ultrasound energy beam 30 and the target tissue mass 16 and calculates the electrical and/or mechanical parameters necessary to maintain the focal zone 34 at the desired site of the target tissue mass 16, and the FUS controller 54 controls the electrical driver 48 and/or the mechanical driver 52 to adjust the thermal dose parameters of the transducer 30. The motion compensation processor 28 may also track the treated and untreated regions of the target tissue mass 16 based on the tracked relative position between the focal zone 34 and the target tissue mass 16.


Although particular embodiments of the present inventions have been shown and described, it will be understood that it is not intended to limit the present inventions to the preferred embodiments, and it will be obvious to those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the present inventions. Thus, the present inventions are intended to cover alternatives, modifications, and equivalents, which may be included within the spirit and scope of the present inventions as defined by the claims.

Claims
  • 1. An image-guided therapy system for applying thermal treatment to a target tissue mass notwithstanding motion of the target tissue mass, the system comprising: a thermal treatment device configured for transmitting a therapeutic energy beam;an imaging device configured for acquiring images of a scanning window around the target tissue mass;a processor configured for tracking relative movement between a focal zone of the energy beam and the target tissue mass in a common coordinate system based on the acquired images; anda controller configured for controlling the thermal treatment device to focus the energy beam on the target tissue mass, the controller being responsive to the processor so as to automatically maintain the focus of the energy beam on the target tissue mass notwithstanding the motion,wherein the imaging device is responsive to the processor so as to adjust a location of the scanning window to correspond with a position of the target tissue mass notwithstanding the motion.
  • 2. The system of claim 1, wherein the processor is configured for tracking the relative movement between the focal zone and the target tissue mass directly from the images.
  • 3. The system of claim 1, wherein the processor is configured for tracking the relative movement between the focal zone and the target tissue mass indirectly, based on a direct correlation between a location of the focal zone and a location of the thermal treatment device, by tracking respective locations of the thermal treatment device and the target tissue mass in a common coordinate system based on the images.
  • 4. The system of claim 3, wherein the processor is configured for tracking the position of the thermal treatment device using a first technique and tracking the position of the target tissue mass using a second technique different from the first.
  • 5. The system of claim 4, wherein the first technique comprises one of template matching between the images and a pre-defined scalable graphical template of the treatment device or localizing the thermal treatment device based on a localization element placed thereon, and the second technique comprises image correlation between the images and reference images of the target tissue mass.
  • 6. The system of claim 1, further comprising a planner configured for automatically constructing, based on the images, a treatment plan comprising a series of treatment sites having associated thermal dose properties.
  • 7. The system of claim 6, wherein construction of the treatment plan is further based on user input.
  • 8. The system of claim 7, wherein the user input comprises at least one user-specified parameter selected from the group consisting of a thermal dose threshold, a maximum allowed energy for each thermal dose, a thermal dose duration for each treatment site, a cooling time between applications of the thermal doses, or a grid density of the treatment sites.
  • 9. The system of claim 7, wherein the user input comprises at least one user-specified physical parameter of the thermal treatment device.
  • 10. The system of claim 7, wherein the user input comprises at least one user-specified parameter related to at least one of an anatomy or an anatomy constraint.
  • 11. The system of claim 7, wherein the user input comprises a specification of a clinical application having at least one default parameter associated therewith.
  • 12. The system of claim 7, wherein the user input comprises, defined in the image, at least one tissue boundary.
  • 13. The system of claim 12, wherein the at least one tissue boundary comprises at least one of a boundary around a target tissue, an obstacle boundary around reflective tissue, or a boundary around tissue sensitive to thermal exposure.
  • 14. The system of claim 12, wherein the at least one tissue boundary is defined three-dimensionally.
  • 15. The system of claim 7, wherein the planner is configured for updating the treatment plan following treatment of any of the treatment sites, taking into account treated regions and untreated regions.
  • 16. The system of claim 15, wherein the planner is responsive to the processor and configured, when updating the treatment plan, to further take into account the relative movement.
  • 17. An image-guided therapy method for applying thermal treatment to a target tissue mass notwithstanding motion of the target tissue mass, the system comprising: focusing a therapeutic energy beam from a thermal treatment device on the target tissue mass;acquiring images of a scanning window around the target tissue mass;tracking relative movement between a focal zone of the energy beam and the target tissue mass in a common coordinate system based on the acquired images;controlling the thermal treatment device so as to automatically maintain the focus of the energy beam on the target tissue mass notwithstanding the motion; andadjusting a location of the scanning window to correspond with a position of the target tissue mass notwithstanding the motion.
  • 18. The method of claim 17, wherein the relative movement between the focal zone and the target tissue mass is tracked directly from the images.
  • 19. The method of claim 17, wherein the relative movement between the focal zone and the target tissue mass is tracked indirectly, based on a direct correlation between a location of the focal zone and a location of the thermal treatment device, from tracked respective locations of the thermal treatment device and the target tissue mass in a common coordinate system based on the images.
  • 20. The method of claim 19, wherein the position of the thermal treatment device is tracked using at least one of template matching between the images and a pre-defined scalable graphical template of the treatment device or localizing the thermal treatment device based on a localization element placed thereon, and the position of the target tissue mass is tracked using image correlation between the images and reference images of the target tissue mass.
  • 21. The method of claim 17, further comprising automatically constructing, based on the images, a treatment plan comprising a series of treatment sites having associated thermal dose properties.
CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a Continuation of U.S. patent application Ser. No. 11/865,662, filed Oct. 1, 2007, which is incorporated herein by reference.

US Referenced Citations (223)
Number Name Date Kind
2795709 Camp Jun 1957 A
3142035 Harris Jul 1964 A
3942150 Booth et al. Mar 1976 A
3974475 Burckhardt et al. Aug 1976 A
3992693 Martin et al. Nov 1976 A
4000493 Spaulding et al. Dec 1976 A
4074564 Anderson Feb 1978 A
4206653 LeMay Jun 1980 A
4339952 Foster Jul 1982 A
4454597 Sullivan Jun 1984 A
4478083 Hassler et al. Oct 1984 A
4505156 Questo Mar 1985 A
4526168 Hassler et al. Jul 1985 A
4537074 Dietz Aug 1985 A
4549533 Cain et al. Oct 1985 A
4554925 Young Nov 1985 A
4662222 Johnson May 1987 A
4858597 Kurtze et al. Aug 1989 A
4865042 Umemura et al. Sep 1989 A
4888746 Wurster et al. Dec 1989 A
4889122 Watmough et al. Dec 1989 A
4893284 Magrane Jan 1990 A
4893624 Lele Jan 1990 A
4937767 Reuschel et al. Jun 1990 A
5209221 Riedlinger May 1993 A
5211160 Talish et al. May 1993 A
5247935 Cline et al. Sep 1993 A
5271400 Dumoulin et al. Dec 1993 A
5275165 Ettinger et al. Jan 1994 A
5291890 Cline et al. Mar 1994 A
5307812 Hardy et al. May 1994 A
5307816 Hashimoto et al. May 1994 A
5318025 Dumoulin et al. Jun 1994 A
5323779 Hardy et al. Jun 1994 A
5327884 Hardy et al. Jul 1994 A
5329930 Thomas, III et al. Jul 1994 A
5339282 Kuhn et al. Aug 1994 A
5368031 Cline et al. Nov 1994 A
5368032 Cline et al. Nov 1994 A
5379642 Reckwerdt et al. Jan 1995 A
5391140 Schaetzle et al. Feb 1995 A
5413550 Castel May 1995 A
5435312 Spivey et al. Jul 1995 A
5443068 Cline et al. Aug 1995 A
5474071 Chapelon et al. Dec 1995 A
5485839 Aida et al. Jan 1996 A
5490840 Uzgiris et al. Feb 1996 A
5507790 Weiss Apr 1996 A
5520188 Hennige et al. May 1996 A
5526814 Cline et al. Jun 1996 A
5549638 Burdette Aug 1996 A
5553618 Suzuki et al. Sep 1996 A
5557438 Schwartz et al. Sep 1996 A
5573497 Chapelon Nov 1996 A
5582578 Zhong et al. Dec 1996 A
5590653 Aida et al. Jan 1997 A
5590657 Cain et al. Jan 1997 A
5601526 Chapelon et al. Feb 1997 A
5605154 Ries et al. Feb 1997 A
5617371 Williams Apr 1997 A
5617857 Chader et al. Apr 1997 A
5643179 Fujimoto Jul 1997 A
5662170 Donovan et al. Sep 1997 A
5665054 Dory Sep 1997 A
5666954 Chapelon et al. Sep 1997 A
5676673 Ferre et al. Oct 1997 A
5687729 Schaetzle Nov 1997 A
5694936 Fujimoto et al. Dec 1997 A
5711300 Schneider et al. Jan 1998 A
5722411 Suzuki et al. Mar 1998 A
5739625 Falcus Apr 1998 A
5743863 Chapelon Apr 1998 A
5752515 Jolesz et al. May 1998 A
5759162 Oppelt et al. Jun 1998 A
5762616 Talish Jun 1998 A
5769790 Watkins et al. Jun 1998 A
5810008 Dekel et al. Sep 1998 A
5831739 Ota Nov 1998 A
5873845 Cline et al. Feb 1999 A
5897495 Aida et al. Apr 1999 A
5938600 Van Vaals et al. Aug 1999 A
5938608 Bieger et al. Aug 1999 A
5947900 Derbyshire et al. Sep 1999 A
5984881 Ishibashi et al. Nov 1999 A
6004269 Crowley et al. Dec 1999 A
6023636 Wendt et al. Feb 2000 A
6042556 Beach et al. Mar 2000 A
6071239 Cribbs et al. Jun 2000 A
6113559 Klopotek Sep 2000 A
6128522 Acker et al. Oct 2000 A
6128958 Cain Oct 2000 A
6193659 Ramamurthy et al. Feb 2001 B1
6242915 Hurd Jun 2001 B1
6246896 Dumoulin et al. Jun 2001 B1
6263230 Haynor et al. Jul 2001 B1
6267734 Ishibashi et al. Jul 2001 B1
6289233 Dumoulin et al. Sep 2001 B1
6309355 Cain et al. Oct 2001 B1
6317619 Boernert et al. Nov 2001 B1
6322527 Talish Nov 2001 B1
6334846 Ishibashi et al. Jan 2002 B1
6374132 Acker et al. Apr 2002 B1
6392330 Zloter et al. May 2002 B1
6397094 Ludeke et al. May 2002 B1
6413216 Cain et al. Jul 2002 B1
6419648 Vitek et al. Jul 2002 B1
6424597 Bolomey et al. Jul 2002 B1
6425867 Vaezy et al. Jul 2002 B1
6428532 Doukas et al. Aug 2002 B1
6461314 Pant et al. Oct 2002 B1
6475150 Haddad Nov 2002 B2
6478739 Hong Nov 2002 B1
6506154 Ezion et al. Jan 2003 B1
6506171 Vitek et al. Jan 2003 B1
6511428 Azuma et al. Jan 2003 B1
6522142 Freundlich Feb 2003 B1
6523272 Morales Feb 2003 B1
6559644 Froundlich et al. May 2003 B2
6566878 Komura et al. May 2003 B1
6582381 Yehezkeli et al. Jun 2003 B1
6599256 Acker et al. Jul 2003 B1
6612988 Maor et al. Sep 2003 B2
6613004 Vitek et al. Sep 2003 B1
6613005 Friedman et al. Sep 2003 B1
6618608 Watkins et al. Sep 2003 B1
6618620 Freundlich et al. Sep 2003 B1
6626854 Friedman et al. Sep 2003 B2
6626855 Weng et al. Sep 2003 B1
6629929 Jago et al. Oct 2003 B1
6645162 Friedman et al. Nov 2003 B2
6652461 Levkovitz Nov 2003 B1
6666833 Friedman et al. Dec 2003 B1
6676601 Lacoste et al. Jan 2004 B1
6679855 Horn et al. Jan 2004 B2
6705994 Vortman et al. Mar 2004 B2
6719694 Weng et al. Apr 2004 B2
6733450 Alexandrov et al. May 2004 B1
6735461 Vitek et al. May 2004 B2
6761691 Tsuzuki Jul 2004 B2
6770031 Hynynen et al. Aug 2004 B2
6770039 Zhong et al. Aug 2004 B2
6788619 Calvert Sep 2004 B2
6790180 Vitek Sep 2004 B2
6824516 Batten et al. Nov 2004 B2
6951540 Ebbini et al. Oct 2005 B2
6961606 DeSilets et al. Nov 2005 B2
7001379 Behl et al. Feb 2006 B2
7077820 Kadziauskas et al. Jul 2006 B1
7094205 Marmarelis Aug 2006 B2
7128711 Medan et al. Oct 2006 B2
7155271 Halperin et al. Dec 2006 B2
7175596 Vitek et al. Feb 2007 B2
7175599 Hynynen et al. Feb 2007 B2
7264592 Shehada Sep 2007 B2
7264597 Cathignol Sep 2007 B2
7267650 Chow et al. Sep 2007 B2
7344509 Hynynen et al. Mar 2008 B2
7377900 Vitek et al. May 2008 B2
7452357 Vlegele et al. Nov 2008 B2
7462488 Madsen et al. Dec 2008 B2
7505805 Kuroda Mar 2009 B2
7505808 Anderson et al. Mar 2009 B2
7510536 Foley et al. Mar 2009 B2
7511501 Wexler Mar 2009 B2
7535794 Prus et al. May 2009 B2
7553284 Vaitekunas Jun 2009 B2
7603162 Danz et al. Oct 2009 B2
7611462 Vortman et al. Nov 2009 B2
7652410 Prus Jan 2010 B2
7699780 Vitek et al. Apr 2010 B2
20010031922 Weng et al. Oct 2001 A1
20020035779 Krieg et al. Mar 2002 A1
20020082589 Friedman et al. Jun 2002 A1
20020188229 Ryaby Dec 2002 A1
20030004439 Pant et al. Jan 2003 A1
20030060820 Maguire et al. Mar 2003 A1
20030187371 Vortman et al. Oct 2003 A1
20040030227 Littrup et al. Feb 2004 A1
20040030251 Ebbini et al. Feb 2004 A1
20040068186 Ishida et al. Apr 2004 A1
20040122323 Vortman et al. Jun 2004 A1
20040199068 Bucholz et al. Oct 2004 A1
20040210134 Hynynen et al. Oct 2004 A1
20040236253 Vortman et al. Nov 2004 A1
20040267126 Takeuchi Dec 2004 A1
20050033201 Takahashi et al. Feb 2005 A1
20050096542 Weng et al. May 2005 A1
20050131301 Peszynski et al. Jun 2005 A1
20050203444 Schonenberger et al. Sep 2005 A1
20050251046 Yamamoto et al. Nov 2005 A1
20060052661 Gannot et al. Mar 2006 A1
20060052701 Carter et al. Mar 2006 A1
20060052706 Hynynen et al. Mar 2006 A1
20060058671 Vitek et al. Mar 2006 A1
20060058678 Vitek et al. Mar 2006 A1
20060079773 Mourad et al. Apr 2006 A1
20060106300 Seppenwoolde et al. May 2006 A1
20060173385 Lidgren et al. Aug 2006 A1
20060184069 Vaitekunas Aug 2006 A1
20060206105 Chopra et al. Sep 2006 A1
20060229594 Francischelli et al. Oct 2006 A1
20070016013 Camus Jan 2007 A1
20070016039 Vortman et al. Jan 2007 A1
20070055140 Kuroda Mar 2007 A1
20070066897 Sekins et al. Mar 2007 A1
20070073135 Lee et al. Mar 2007 A1
20070098232 Matula et al. May 2007 A1
20070167781 Vortman et al. Jul 2007 A1
20070197918 Vitek et al. Aug 2007 A1
20070219470 Talish et al. Sep 2007 A1
20080027342 Rouw et al. Jan 2008 A1
20080031090 Prus et al. Feb 2008 A1
20080033278 Assif Feb 2008 A1
20080082026 Schmidt et al. Apr 2008 A1
20080108900 Lee et al. May 2008 A1
20080183077 Moreau-Gobard et al. Jul 2008 A1
20080228081 Becker et al. Sep 2008 A1
20080249408 Palmeri et al. Oct 2008 A1
20080292194 Schmidt et al. Nov 2008 A1
20080312562 Routh et al. Dec 2008 A1
20090088623 Vortman et al. Apr 2009 A1
20100056962 Vortman et al. Mar 2010 A1
20110034800 Vitek et al. Feb 2011 A1
Foreign Referenced Citations (32)
Number Date Country
4345308 Feb 2001 DE
1132054 Sep 2001 EP
0151073 Nov 2005 EP
1774920 Apr 2007 EP
1790384 May 2007 EP
2806611 Sep 2001 FR
11313833 Nov 1999 JP
WO-910005 Jan 1991 WO
WO-9852465 Nov 1998 WO
WO-0031614 Jun 2000 WO
WO-0158337 Aug 2001 WO
WO-0166189 Sep 2001 WO
WO-0180709 Nov 2001 WO
WO-02058791 Aug 2002 WO
WO-03013654 Feb 2003 WO
WO-03097162 Nov 2003 WO
WO-03098232 Nov 2003 WO
WO-2004093686 Nov 2004 WO
WO-2005-58029 Jun 2005 WO
WO-2006018837 Feb 2006 WO
WO-2006025001 Mar 2006 WO
WO-2006087649 Aug 2006 WO
WO-2006119572 Nov 2006 WO
WO-2007073551 Jun 2007 WO
WO-2008039449 Apr 2008 WO
WO-2008050278 May 2008 WO
WO-2008075203 Jun 2008 WO
WO-2008119054 Oct 2008 WO
WO-2009044276 Apr 2009 WO
WO-2009055587 Apr 2009 WO
WO-2009094554 Jul 2009 WO
WO-2011015949 Feb 2011 WO
Non-Patent Literature Citations (50)
Entry
Baudouin Denis de Senneville, Charles Mougenot, Chrit T.W. Moonn, Real-Time Adaptive Methods for Treatment of Mobile Organs by MRI-Controlled High-Intensity Focused Ultrasound. Magnetic Resonance in Medicine 57: 319-330 (2007), Wiley-Liss, Inc.
International Search Report and Written Opinion for PCT/IB2008/003069—actual completion date of Apr. 2, 2009 (11 pages).
Cline, H.E., et al., “MR Temperature Mapping of Focused Ultrasound Surgery”, Magnetic Resonance in Medicine, Academic Press, Duluth, MN, US, vol. 32, No. 6, 1 Kiem 1994 (Jan. 6, 1994), ISSN; 0740-3194, p. 626, figure 6.
Cain, C. and Umemura, S., “Concentric-Ring and Sector-Vortex Phased-Array Applicators for Ultrasound Hyperthermia, ” IEEE Transactions on Microwave Theory and Techniques, vol. MTT-34, No. 5, pp. 542-551 (May 1986).
Fjieid, T. and Hynyen, K., The Combined Concentric-Ring and Sector-Vortex Phased Array for MRI Guided Ultrasound Surgery, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, vol. 44, No. 5, pp. 1157-1167 (Sep. 1997).
Botros et al., “A hybrid computational model for ulrasound phased-array heating in presence of strongly scattering obstacles, ” IEEE Trans. On Biomed. Eng., vol. 44, No. 11, pp. 1039-1050 (Nov. 1997).
Chen et al., “MR Acoustic Radiation Force Imaging: Comparison of Encoding Gradients.”
Cline et al., “Focused US system for MR imaging-guide tumor ablation,” Radiology, vol. 194, No. 3, pp. 731-738 (Mar. 1995).
Cline et al., “Simultaneous magnetic resonanace phase and magnitude temperature maps in muscle,” Magnetic Resonance in Medicine, vol. 35, No. 3, pp. 309-315 (Mar. 1996).
Daum et al., “Design and evaluation of a feedback based phased array system for ultrasound surgery,” IEE Trans Ultrason. Ferroelec. Freq. Control, vol. 45, No. 2, pp. 431-434 (1998).
Herbert et al., “Energy-based adaptive focusing of waves: application to ultrasonic transcranial therapy.” 8th Intl. Symp. On therapeutic Ultrasound.
Huber et al., “A New Noninvasive Approach in Breast Cancer Therapy Using Magnetic Resonance Imaging-Guided Focussed Ultrasound Surgery.” Cancer Research 61, 8441-8447 (Dec. 2001).
International Preliminary Report on Patentability in International Patent Application No. PCT/IB2004/001512, mailed Dec. 8, 2005.
International Search Report and Written Opinion in International Patent Application No. PCT/IB2004/001498, dated Aug. 31, 2004.
International Search Report and Written Opinion in International Patent Application No. PCT/IB2005/002273, mailed Dec. 20, 2005.
International Search Report and Written Opinion in International Patent Application No. PCT/IB2005/002413, mailed Nov. 22, 2005.
International Search Report and Written Opinion in International Patent Application No. PCT/IB2006/001641, mailed Sep. 25, 2006.
International Search Report and Written Opinion in International Patent Application No. PCT/IB2006/003300, mailed Feb. 14, 2008.
International Search Report and Written Opinion in International Patent Application No. PCT/IB2007/001079, mailed Dec. 10, 2007.
International Search Report and Written Opinion in International Patent Application No. PCT/IB2007/002134, mailed Dec. 13, 2007.
International Search Report and Written Opinion in International Patent Application No. PCT/IB2007/002140, mailed Dec. 29, 2008.
Jolesz et al., “Integration of interventional MRI with computer-assisted surgery,” J. Magnetic Resonance Imaging, 12:69-77 (2001).
Kohler et al., “Volumetric HIFU Ablation guided by multiplane MRI thermometry,” 8th Intl. Symp. On Therapeutic Ultrasound, edited by E.S. Ebbini, U. of Minn. (Sep. 2009).
Kowalski et al., “Optimization of electromagnetic phased-arrays for hyperthermia via magnetic resonance temperature estimation,” IEEE Trans. On Biomed. Eng., vol. 49, No. 11, pp. 1229-1241 (Nov. 2002).
Maxwell et al. “Noninvasive thrombolysis using pulsed ultrsound cavitation therapy—Histotripsy,” Abstract U.S. Natl. Lib. Of Med., NIH, Ultrsound Med. Biol. (Oct. 23, 2009).
McDannold et al., “MRI evaluation of thermal ablation of tumors and focused ultrasounds,” JMRI vol. 8, No. 1, pp. 91-100 (1998).
McDannold et al., “Magnetic resonance acoustic radiation force imaging,” Med. Phys. vol. 35, No. 8, pp. 3748-3758 (Aug. 2008).
Medei et al., “Sonothrombolysis: An emerging modality for the managemnet of stroke,” Neurosurgery, vol. 65, No. 5, pp. 979-993.
Mougenot et al., “MR monitoring of the near-filed HIFU heating,” 8th Intl. Symp. On Therapeutic Ultrasound, edited by E.S. Ebbini, U. of Minn. (Sep. 2009).
Partial International Search Report and Written Opinion in International Patent Application No. PCT/IB2007/001079, dated Sep. 25, 2007.
Vimeus et al., “Real-time control of focused ultrasound heating based on rapid MR thermometry,” Investig. Radiology, vol. 43, No. 3, pp. 190-193.
Vykhodtseva et al., “MRS detection of the thermal effects of focused ultrasound on the brain,” Ultrasound in Med. & Biol., vol. 26, No. 5, pp. 871-880 (2000).
Written Opinion in International Patent Application No. PCT/IL01/00340, mailed Feb. 24, 2003.
Written Opinion in International Patent Application No. PCT/IL02/00477, mailed Feb. 25, 2003.
Written Opinion in International Patent Application No. PCT/IB03/05551, mailed Sep. 10, 2004.
“How is Ablatherm treatment performed?” http://www.edap-hifu.com/eng/physicians/hifu/3c—treatment—treat-description.htm, accessed Jan. 3, 2003.
“What is HIFU? HIFU: High Intensity Focused Ultrasound,” http://www.edap-hifu.com/eng/physicians/hifu2a—hifu—overview.htm, accessed Jan. 3, 2003.
“What are the physical principles?” http://www.edap-hifu.com/eng/physicians/hifu/2c—hifu—physical.htm, accessed Jan. 3, 2003.
“How does HIFU create a lesion?” http://www.edap-hifu.com/eng/physicians/hifu/2d—hifu—lesion.htm, accessed Jan. 3, 2003.
“Prostate Cancer Phase I Clinical Trials Using High Intensity Focused Ultrsound (HIFU),” Focus Surgery, http://www.focus-surgery.com/PCT%20Treatment%20with%20HIFU.htm, accessed Jan. 3, 2003.
“Abstract” Focus Surgery, http://www.focus-surgery.com/Sanghvi.htm, accessed Jan. 3, 2003.
Exablate 2000 Specificaiton, InSightec, Ltd. (2 pages).
FDA Approves Exablate 2000 as Non-Invasive surgery for Fibroids, Oct. 22, 2004.
International Search Report and Written Opinion in International Patent Application No. PCT/IB2010/000189, mailed Jun. 1, 2010.
International Search Report for PCT/IB03/05551 completion date Mar. 2, 2004 (5 pages).
International Search Report and Written Opinion in International Patent Application No. PCT/IB2010/000971, mailed Jul. 29, 2010 (9 pages).
McGough, et al., “Direct Computation of Ultrasound Phased-Array Driving Signals from a Specified Temperature Distribution for Hyperthermia,” IEEE Trans. On Biomedical Engineering, vol. 39, No. 8, pp. 825-835 (Aug. 1992).
International Preliminary Report on Patentbility mailed on Apr. 15, 2010 for International Application No. PCT/IB2008/003069 (10 pages).
McDannold, et al., “Quality Assurance and System Stability of a Clinical MRI-guided focused ulatrsound system: Four-year experience,” Medical Physics, vol. 33, No. 11, pp. 4307-4313 (Oct. 2006).
Office Action mailed Dec. 26, 2011 for Chinese Patent Application No. 200880109807.1 (translation) 13 pages.
Related Publications (1)
Number Date Country
20130030283 A1 Jan 2013 US
Continuations (1)
Number Date Country
Parent 11865662 Oct 2007 US
Child 13555338 US